Kerogens are extracted from deep shales to study pyrolysis of deep shale samples. The 2D molecular models of kerogens are obtained by a series of physical and chemical experiments by which the macromolecular models of kerogens are constructed. Then, the reasonable 3D macromolecular models are established by molecular mechanics and global energy minimization. The effects of temperature and heating rate on the chemical kinetics of kerogen pyrolysis are studied using reactive force field (ReaxFF). The hybrid molecular dynamics/force-biased Monte Carlo (MD/fbMC) approach is used to simulate the pyrolytic process at the experimental temperature, which is lower than the conventional one. The gaseous products and residues obtained by the simulations agree with the experimental results, which means a reliable simulation method for pyrolysis at experimental temperature is provided. This study constructs the rational macromolecular models of kerogen by experiments, and proposes the mechanisms of typical reactions of kerogen pyrolysis, which may help in understanding the formation of shale oil and gas.