A system of thermal decompositions is presented providing an adequate description of the mechanism of oil shale pyrolysis. Specific rate constants for systems of thermal decompositions studied are estimated via a least-squares procedure. The estimates thus obtained provide a mathematical model of the corresponding thermal decomposition system. A model that generates a good fit to published Bureau of Mines experimental data, in comparison with other models studied, is based on a thermal decomposition system consisting of two first-order reactions and three second-order reactions. Values computed for the specific rate constants corresponding to this system are shown to exhibit a temperature dependence that is in good agreement with the Arrhenius equation.