A Review of Molecular Models for Gas Adsorption in Shale Nanopores and Experimental Characterization of Shale Properties

Shale gas, as a promising alternative energy source, has received considerable attention because of its broad resource base and wide distribution. The establishment of shale models that can accurately describe the composition and structure of shale is essential to perform molecular simulations of gas adsorption in shale reservoirs. This Review provides an overview of shale models, which include organic matter models, inorganic mineral models, and composite shale models. Molecular simulations of gas adsorption performed on these models are also reviewed to provide a more comprehensive understanding of the behaviors and mechanisms of gas adsorption on shales. To accurately understand the gas adsorption behaviors in shale reservoirs, it is necessary to be aware of the pore structure characteristics of shale reservoirs. Thus, we also present experimental studies on shale microstructure analysis, including direct imaging methods and indirect measurements. The advantages, disadvantages, and applications of these methods are also well summarized. This Review is useful for understanding molecular models of gas adsorption in shales and provides guidance for selecting experimental characterization of shale structure and composition.

Release of methane from nanochannels through displacement using CO2

In this work, we investigate the release of methane in quartz nanochannels through the method of displacement using carbon dioxide. Molecular dynamics (MD) simulations and theoretical analysis are performed to obtain the release percentage of methane for nanochannels of various diameters. It is found that both the pressure of CO2 and the channel size affect the release percentage of methane, which increases with increasing pressure of CO2 and channel diameter. Without CO2, the majority of methane molecules are adsorbed by the channel surface. When CO2 is injected into the channel, CO2 molecules replace many methane molecules due to the relatively strong molecular interactions between CO2 and the channel, which leads to the desorption of methane, reduces the energy barrier for the transport of methane, and consequently increases the release rate. Theoretical predictions using the kinetic energy of methane and the energy barrier inside the channel are also conducted, which are in good agreement with the MD simulations.

Molecular Investigation on the Displacement Characteristics of CH4 by CO2, N2 and Their Mixture in a Composite Shale Model

The rapid growth in energy consumption and environmental pollution have greatly stimulated the exploration and utilization of shale gas. The injection of gases such as CO2, N2, and their mixture is currently regarded as one of the most effective ways to enhance gas recovery from shale reservoirs. In this study, molecular simulations were conducted on a kaolinite–kerogen IID composite shale matrix to explore the displacement characteristics of CH4 using different injection gases, including CO2, N2, and their mixture. The results show that when the injection pressure was lower than 10 MPa, increasing the injection pressure improved the displacement capacity of CH4 by CO2. Correspondingly, an increase of formation temperature also increased the displacement efficiency of CH4, but an increase of pore size slightly increased this displacement efficiency. Moreover, it was found that when the proportion of CO2 and N2 was 1:1, the displacement efficiency of CH4 was the highest, which proved that the simultaneous injection of CO2 and N2 had a synergistic effect on shale gas production. The results of this paper will provide guidance and reference for the displacement exploitation of shale gas by injection gases.

Functionalized multiscale visual models to unravel flow and transport physics in porous structures

The fluid flow, species transport, and chemical reactions in geological formations are the chief mechanisms in engineering the exploitation of fossil fuels and geothermal energy, the geological storage of carbon dioxide (CO₂), and the disposal of hazardous materials. Porous rock is characterized by a wide surface area, where the physicochemical fluid-solid interactions dominate the multiphase flow behavior. A variety of visual models with differences in dimensions, patterns, surface properties, and fabrication techniques have been widely utilized to simulate and directly visualize such interactions in porous media. This review discusses the six categories of visual models used in geological flow applications, including packed beds, Hele-Shaw cells, synthesized microchips (also known as microfluidic chips or micromodels), geomaterial-dominated microchips, three-dimensional (3D) microchips, and nanofluidics. For each category, critical technical points (such as surface chemistry and geometry) and practical applications are summarized. Finally, we discuss opportunities and provide a framework for the development of custom-built visual models.

In situ study on interactions between hydroxyl groups in kaolinite and re-adsorption water

The interactions between O–H groups in kaolinite and re-adsorption water is an important aspect that should be considered in the hydraulic fracturing method for the production of shale gas, because the external water adsorbed by kaolinite in shale would significantly affect the desorption of methane. In this study, the interactions were investigated via changing the amount of O–H groups and re-adsorption water in kaolinite by heating treatment and water re-adsorption. To overcome the overlap of IR vibration bands of the O–H functional groups in H₂O and those in parent kaolinite, kaolinite samples with D₂O re-adsorption were prepared by drying the H₂O from raw kaolinite and soaking the dried kaolinite in D₂O. The interactions between O–H groups in kaolinite and D₂O molecules were investigated by in situ DRIFT and TG-MS. The results demonstrated that the vibration at 3670 ± 4 cm⁻¹ in the DRIFT spectra could be due to the outer O–H groups of the octahedral sheet on the upper surface of the kaolinite microcrystal structure, rather than a type of inner-surface O–H group. All types of O–H groups, including the inner O–H groups in kaolinite, could be transformed into O–D groups after D₂O re-adsorption at room temperature. The inner-surface O–H groups in kaolinite are the most preferred sites for D₂O re-adsorption; thus, they would be the key factor for studying the effect of re-adsorption water on methane desorption. When the temperature increased from 100 °C to 300 °C, two layers of kaolinite slipped away from each other, resulting in the transformation of inner-surface O–H groups into outer O–H groups. Thus, the temperature range of 100 to 300 °C was suggested for the heat treatment of kaolinite to decrease the content of inner-surface O–H groups; thereby, the amount of re-adsorption water was reduced. However, to thoroughly remove the re-adsorption water, a temperature higher than 650 °C should be used.

Because two layers slipped away from each other, inner-surface O–H was transformed into outer O–H during heating from 100–300 °C. Re-adsorption water could be thoroughly removed at 650 °C.

A new method to determine varying adsorbed density based on Gibbs isotherm of supercritical gas adsorption

In the calculation of the absolute adsorption of supercritical gas adsorbed on the microporous materials, most existing methods regard the adsorbed density as a constant, which is very unreasonable. In this study, an extended pressure point method combined with Langmuir adsorption model is proposed in which the varying adsorbed density under different pressures is considered at the same time. The utility of the proposed method to correlate accurately the experimental data for supercritical gas adsorption system is demonstrated by high-pressure methane adsorption measurements on two groups of shale samples. Taking advantage of the proposed method, we can obtain the adsorbed density and the adsorbed volume corresponding to different pressures. Compared with the conventional methods under the assumption of fixed and parameterized adsorbed density, the proposed method yields better fitting results with the experimental data. Our work should provide important fundamental understandings and insights into the supercritical gas adsorption system.

Investigation of dynamical properties of methane in slit-like quartz pores using molecular simulation

The dynamical properties of adsorption media confined in micropores play an important role in the adsorptive separation of fluids. However, a problem is that it is difficult to directly use approaches based on experimental measurements. Molecular simulation has been an effective tool for investigating the diffusion of fluids on the microscale in recent years. In this work, the diffusion properties of methane in quartz were mainly investigated from a microscale viewpoint using MD (molecular dynamics) methods, and this paper primarily discusses the influence of parameters such as pressure, temperature, pore size and water content on the diffusion and thermodynamic parameters of methane in slit-like quartz pores. The results demonstrate that the transport ability of quartz pores decreases with an increase in pressure in pores of a fixed size at a certain temperature and increases with an increase in pore size or temperature at a fixed pressure, which is related to changes in the interaction between methane molecules and quartz. In the pressure range used in the simulation, the average isosteric heat of adsorption of methane increases with an increase in pressure and is in the range of 6.52–10.794 kJ mol⁻¹. Therefore, the gas adsorption behavior is classed as physical adsorption because the heat of adsorption is significantly lower than the minimum heat of gas adsorption for chemisorption. The increase in the total adsorption entropy is caused by an increase in temperature due to an increase in internal energy, which brings about a reduction in the interactions between gas molecules and walls of quartz. However, with an increase in pore size the total adsorption entropy increases, for which an explanation may be that in pores of a larger size methane molecules are adsorbed at higher-energy sites and generate a higher isosteric heat, which causes a reduction in interactions between the adsorbate and adsorbent. Regarding the influence of different water contents on the diffusion of methane, it was demonstrated that with an increase in moisture the mobility of methane molecules initially increases and then decreases, which is related to the distance between gas molecules.

The dynamical properties of adsorption media confined in micropores play an important role in the adsorptive separation of fluids.

Geo-material surface modification of microchips using layer-by-layer (LbL) assembly for subsurface energy and environmental applications

A key constraint in the application of microfluidic technology to subsurface flow and transport processes is the surface discrepancy between microchips and the actual rocks/soils. This research employs a novel layer-by-layer (LbL) assembly technology to produce rock-forming mineral coatings on microchip surfaces. The outcome of the work is a series of 'surface-mimetic micro-reservoirs (SMMR)' that represent multi-scales and multi-types of natural rocks/soils. For demonstration, the clay pores of sandstones and mudrocks are reconstructed by representatively coating montmorillonite and kaolinite in polydimethylsiloxane (PDMS) microchips in a wide range of channel sizes (width of 10-250 μm, depth of 40-100 μm) and on glass substrates. The morphological and structural properties of mineral coatings are characterized using a scanning electron microscope (SEM), optical microscope and profilometer. The coating stability is tested by dynamic flooding experiments. The surface wettability is characterized by measuring mineral oil-water contact angles. The results demonstrate the formation of nano- to micro-scale, fully-covered and stable mineral surfaces with varying wetting properties. There is an opportunity to use this work in the development of microfluidic technology-based applications for subsurface energy and environmental research.

Molecular insights into competitive adsorption of CO2/CH4 mixture in shale nanopores

In the present study, competitive adsorption behaviour of supercritical carbon dioxide and methane binary mixture in shale organic nanopores was investigated by using grand canonical Monte Carlo (GCMC) simulations. The model was firstly validated by comparing with experimental data and a satisfactory agreement was obtained. Then the effects of temperature (298–388 K), pressure (up to 60 MPa), pore size (1–4 nm) and moisture content (0–2.4 wt%) on competitive adsorption behaviour of the binary mixture were examined and discussed in depth. It is found that the adsorption capacity of carbon dioxide in shale organic nanopores is much higher than that of methane under various conditions. The mechanism of competitive adsorption was discussed in detail. In addition, the results show that a lower temperature is favorable to both the adsorption amount and selectivity of CO₂/CH₄ binary mixture in shale organic nanopores. However, an appropriate CO₂ injection pressure should be considered to take into account the CO₂ sequestration amount and the exploitation efficiency of shale gas. As for moisture content, different influences on CO₂/CH₄ adsorption selectivity have been observed at low and high moisture conditions. Therefore, different simulation technologies for shale gas production and CO₂ sequestration should be applied depending on the actual moisture conditions of the shale reservoirs. It is expected that the findings in this work could be helpful to estimate and enhance shale gas resource recovery and also evaluate CO₂ sequestration efficiency in shale reservoirs.

Competitive adsorption behaviour of CO₂/CH₄ mixture in shale slit nanopores under various geological conditions was explored by molecular simulations.