The kinetics of complex reaction systems were studied on molecular level with the combined method of Monte Carlo simulation and Structural Oriented Lumping by focusing on deep catalytic cracking (DCC) process, the model parameters were optimized by means of routine analytic data of a DCC unit. A model was established to transform the feedstock of the complex reaction systems such as DCC to 1000-10000 pseudo-molecules with the Monte Carlo simulation and every molecule was expressed by 19 attributes. The results of model simulation showed that these pseudo-molecules reflected the characteristics of feedstock very well and their average properties gave a good agreement with the plant data.