Study of subcritical and supercritical gas adsorption behavior in different nanopore systems in shale using lattice Boltzmann method

Adsorption has been studied for decades in either petroleum or chemical industry. However, the detailed adsorption behavior is yet not well understood in nanoporous media with complex pore structures and surface properties, especially for subcritical gases when capillary condensation might occur. Shale (mudrock) has complex pore systems as a result of heterogeneous mineralogy and diagenesis, and interparticle (interP) and intraparticle (intraP) pores are two main components. Understanding gas adsorption and phase transition behavior in different pore systems is essential for petrophysical characterization, reserve estimation and production forecast. In this study, we use lattice Boltzmann method (LBM) to study nitrogen (subcritical) and methane (supercritical) adsorption behavior in reconstructed shale nanopore systems. The model accounts for the van der Walls (vdW) forces between fluid molecules by incorporating a modified pseudopotential model based on Peng-Robinson equation of state (PR-EOS). Rock-fluid interaction is also considered by introducing phenomenological surface forces, which take two different forms for subcritical and supercritical gases respectively. The correlation between parameters defined in LBM and physical properties (i.e. density, pressure, isosteric heat of adsorption) is established systematically. Using the developed model, we observe quantitative agreement between LBM and lattice density functional theory (LDFT)/experimental data in terms of adsorption isotherms and density distributions. We then compare adsorption behavior in two 'model' interP and intraP pore systems and a synthesized polyethylene porous medium with more complex pore shapes. We observe three adsorption stages for subcritical gas adsorption: a) early capillary condensation near grain contacts or sharp corners, forming isolated pendular rings or fluid clusters; b) growing, merging, and spreading of the isolated condensed phase, forming circular gas/liquid meniscus; c) the condensed phase becomes fully-connected across the geometry, leaving isolated gas bubbles. We observe a signature behavior for the adsorption isotherm of intraP pore system, where a stepped feature occurs around relative pressure of 0.8, and the adsorption uptake for the polyethylene porous medium is bounded by the upper and lower limits given by the interP and intraP pore systems. On the other hand, we find that supercritical methane adsorption in mesoporous media is less sensitive to pore shape, but rather controlled by the specific surface area. The existence of micropores also affects the adsorption uptake.