Oil & Gas Geology ›› 2025, Vol. 46 ›› Issue (3): 995-1005.doi: 10.11743/ogg20250320

• Methods and Technologies • Previous Articles     Next Articles

Molecular simulation of methane adsorption/desorption characteristics of deformed kerogen in shales

Liang HUANG1,2,3(), Zhenyao XU2, Qiujie CHEN2, Baohua TIAN2, Qin YANG2, Xinni FENG2, Zishuo QU2, Zhe YANG2, Sirun AN2   

  1. 1.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation,Chengdu University of Technology,Chengdu,Sichuan 610059,China
    2.College of Energy,Chengdu University of Technology,Chengdu,Sichuan 610059,China
    3.College of Chemistry,University of California,Berkeley,California,94720,USA
  • Received:2024-10-16 Revised:2025-05-15 Online:2025-06-30 Published:2025-06-26

Abstract:

Current studies on the methane adsorption/desorption behavior in kerogen are primarily conducted based on fixed pore structures. Consequently, the influence patterns and mechanisms of the pore deformation of kerogen on methane adsorption/desorption remain poorly understood. Using molecular simulation, we explore both the methane adsorption/desorption behavior in shale kerogen and the coupled pore deformation characteristics. By constructing and characterizing a nanopore model of kerogen, we investigate the influence of pressure on kerogen deformation. Using a coupling method combining Grand Canonical Monte Carlo (GCMC) and molecular dynamics (also referred to as the GCMC-MD method), we analyze the coupling characteristics of methane adsorption and kerogen deformation. Additionally, through methane desorption simulation using a kerogen model under high-pressure adsorption equilibrium, we investigate the characteristics and causes of methane desorption hysteresis. The results indicate that ultramicropores in Type ⅡD kerogen serve as the primary migration pathways, while isolated ultramicropores also have the potential to act as migration pathways due to the influence of pore deformation. The pore structures of kerogen generally are reversible during the pressure increase and decrease cycles. The methane adsorption-induced kerogen deformation can be divided into two stages: volume expansion and pore reorganization. The structural deformation of kerogen enhances the methane adsorption capacity, while the selective adsorption of methane onto high-affinity heterogeneous sites in kerogen leads to pronounced changes in pore space. Small desorption hysteresis loops are observed within the deformed structures of kerogen. The methane desorption hysteresis is primarily caused by pore deformation rather than the thermodynamic differences between adsorption and desorption and the capillary condensation effect.

Key words: adsorption, desorption hysteresis, deformation, molecular simulation, kerogen, shale gas

CLC Number: