This study explores the potential of organic-rich formations for geo-storage of hydrogen and carbon dioxide using a novel molecular simulation workflow combing Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) to delineate the interplay between pore structure, kerogen maturity, and the resulting adsorption–desorption behavior. The investigation utilizes two representative kerogen models – immature (II-A) and overmature (II-D) – to account for the inherent heterogeneity of organic-rich formations. The simulations demonstrate a direct correlation between kerogen maturity and gas sorption capacity. The overmature model, with its larger intermolecular spaces, exhibits a higher capacity for both H2 and CO2 compared to the immature one. Notably, CO2 displays a significantly stronger affinity towards the organic material, resulting in a much higher sorption capacity compared to H2 for both kerogen types. The key finding lies in the contrasting desorption profiles for H2 and CO2. H2 desorption exhibits minimal hysteresis, suggesting readily recoverable gas with minimal trapping within the kerogen structure. This signifies the potential for efficient H2 storage with minimal losses during retrieval stages. Conversely, CO2 desorption displays a pronounced hysteresis loop, indicating significant gas retention within the formation even at lower pressure. This behavior aligns perfectly with the objective of CO2 sequestration, where minimizing leakage risk is paramount. These contrasting effects of hysteresis on H2 storage and CO2 sequestration highlight the unique suitability of organic-rich formations for both applications. Additionally, the study suggests a potential link between kerogen maturity and the magnitude of the hysteresis loop. Such insights could be instrumental in selecting optimal formations for targeted geo-storage operations.
Read full abstract