Abstract: Using the Forcite module in MS (Materials Studio 2017) software, the energy minimization molecular dynamics simulation was performed on the self-constructed two-dimensional structural model of Fushun oil shale kerogen, and the initial optimized structure of kerogen was obtained through the energy optimization process. Then, molecular dynamics annealing simulations were performed to obtain a global energy optimization configuration, ie a three-dimensional structural model of oil shale kerogen molecules. Based on the density functional theory of quantum mechanics simulation method, a three-dimensional structural model of kerogen dynamics, bond energy, bond level, charge density and other parameters were calculated, and the chemical active sites were analyzed. The microchemical evolution mechanism of kerogen pyrolysis was discussed. And then, the reactivity was predicted.