Understanding pure H2 and H2/CH4 adsorption and diffusion in earth materials is one vital step toward a successful and safe H2 storage in depleted gas reservoirs. Despite recent research efforts such understanding is far from complete. In this work we first use Nuclear Magnetic Resonance (NMR) experiments to study the NMR response of injected H2 into Duvernay shale and Berea sandstone samples, representing materials in confining and storage zones. Then we use molecular simulations to investigate H2/CH4 competitive adsorption and diffusion in kerogen, a common component of shale. Our results indicate that in shale there are two H2 populations, i.e., free H2 and adsorbed H2, that yield very distinct NMR responses. However, only free gas presents in sandstone that yields a H2 NMR response similar to that of bulk H2. About 10 % of injected H2 can be lost due to adsorption/desorption hysteresis in shale, and no H2 loss (no hysteresis) is observed in sandstone. Our molecular simulation results support our NMR results that there are two H2 populations in nanoporous materials (kerogen). The simulation results also indicate that CH4 outcompetes H2 in adsorption onto kerogen, due to stronger CH4-kerogen interactions than H2-kerogen interactions. Nevertheless, in a depleted gas reservoir with low CH4 gas pressure, about ∼30 % of residual CH4 can be desorbed upon H2 injection. The simulation results also predict that H2 diffusion in porous kerogen is about one order of magnitude higher than that of CH4 and CO2. This work provides an understanding of H2/CH4 behaviors in deleted gas reservoirs upon H2 injection and predictions of H2 loss and CH4 desorption in H2 storage.
Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO2 preferentially into their interlayers and therefore potentially act as a material for CO2 separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO2 into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO2 intercalation. Our free energy calculations indicate that CO2 favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO2/H2O ratio in clay interlayers as compared to that in bulk water. CO2 molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO2 intercalation. The structural charge distribution significantly affects the amount of CO2 intercalated in the interlayers. These results provide a mechanistic understanding of CO2 intercalation in clays for CO2 separation, capture, and storage.
Understanding the formation of H2CO3 in water from CO2 is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO2 to H2CO3 in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO2 conversion to H2CO3 in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Here our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO2 concentration.