Molecular Modeling of the Volumetric and Thermodynamic Properties of Kerogen: Influence of Organic Type and Maturity

Abstract

Molecular modeling is applied to a representative array of kerogens for the purpose of obtaining quantitative predictions of thermodynamic properties from quantum mechanics and volumetric properties from molecular dynamics. The kerogen model units (175−260 carbon atoms) have been built in the MedeA environment from the sole consideration of the elemental analysis and functional group analysis documented in the work of Exxon and IFP-EN scientists [Kelemen, S. R., et al., Energy Fuels, 2007, 21 (3), pp 1548−1561. The density results are in good agreement with the well-documented trends of kerogen density with thermal maturity and organic type. The heat capacity in the ideal gas state is predicted to increase as a function of temperature, as obtained from quantum mechanics at the semiempirical level (MOPAC-PM7). This result is in quantitative agreement with experimental heat capacity data on type I kerogen and on coal. This behavior appears clearly as a nonclassical feature, because of the quantization of energy levels in molecular vibrations. Also, the residual heat capacity estimated from molecular dynamics appears subordinate, compared with the ideal heat capacity evaluated from quantum mechanics. The change from negative to positive standard enthalpy of formation when changing from low-maturity kerogen to high-maturity kerogen is also predicted in agreement with correlative methods based on numerous experimental data from coals and fossil fuels. Kerogen model units are available for download free of charge in .xyz or .sci formats from www.materialsdesign. com/science/structures/kerogens_and_coals.

Nanostructural control of methane release in kerogen and its implications to wellbore production decline

“…It is representative of over-mature kerogen (Type II-D) found in the Duvernay organic-rich marine shale formation and similar to that found in Barnett shale. In the work of Ungerer et al 293031,. the kerogen molecule was modeled by implementing the PCFF+ force field32, which describes atomic dispersion-repulsion interactions using the Lennard Jones 6–9 potential.…”

Fast mass transport of oil and supercritical carbon dioxide through organic nanopores in shale

“…On the basis of the documented results for elemental analysis and functional group analysis, Ungerer et al [67] constructed molecular model units for different types of kerogen at various maturity levels. Recent molecular simulation shows that kerogen density, specific adsorption, and other properties estimated from their model are consistent with the experimental data [68].…”

“…As suggested by the developer of the kerogen model, we used all-atom pcff + force field to describe the interatomic potentials [67,68]. The LJ interactions between different species were determined through the Waldman-Hagler combining rules and the terms beyond the cutoff distance of 9.5 Å were also included via long-range tail correction.…”

Study on the Adsorption, Diffusion and Permeation Selectivity of Shale Gas in Organics

“…Based on the C/H and C/O ration of I, II and III kerogen, Ungerer et al applied a molecular simulation method to establish a molecular structure model. The thermodynamic properties of different types of organics such as maturity, density, compressibility, swelling and adsorption isotherm were studied [20]. Hu [49] studied the change of gas and water distribution in organic pore by molecular dynamics simulation.…”