A kinetic model for the catalytic reforming process was developed at the molecular level. The structure-oriented lumping (SOL) method was modified to fit catalytic reforming, and a structural increment was created to distinguish between the xylene isomers; the xylene isomers are important reforming products. To model the governing chemistry, a reaction network involving 140 molecules and 697 reactions was generated in terms of 16 reaction rules. Then, a mathematical model of the naphtha reforming process was developed and included a reaction kinetic model and a reformer reactor model. The linear free energy relation (LFER) method was used to reduce the kinetic parameters. Experiments under different conditions were conducted to optimize the parameters. The results showed that the calculated values of the product yield and molecular composition were in good agreement with the experimental data. The conversion of hydrocarbon molecules was revealed by calculating the distribution of the molecules in the reactor. Sensitivity analysis was performed to investigate the effect of the reaction temperature on the molecular transformations.