Mechanism analysis of CO 2 enhanced shale gas exploitation : Competitive Adsorption between CO 2 and CH 4 in Shale Kerogen NanoPores

Kecheng Zeng,Ruina Xu,Bo Zhou,Peixue Jiang

2018-01-01

Abstract:Using carbon dioxide as a displacing fluid to increase the recovery ratio, in the enhanced oil recovery (EOR) technique and in shale gas development, is attractive since such potential technique also contributes to carbon dioxide geological storage. Supercritical methane in shale gas reservoirs is either adsorbed on solid shale matrix pore surfaces or stored freely in the pore void spaces. All the shale matrix pores are in nano-scale, in which the organic kerogen pore radii are mostly less than 10 nm. The nano-scale kerogen pores provide huge specific surface area for shale reservoir. The physical property of supercritical methane and carbon dioxide in the nano-pores, such as density and viscosity, is dominated by the influence of complicated solid surface. Considering the complex geometry and chemical components of shale organic nano-pores, it’s nearly impossible to conduct experiment and theoretical analysis. Numerical simulation is a feasible method to research the physical behavior of fluid in such a small scale confined space under high pressure and temperature. Comparing with the methane fluid stored in shale kerogen organic pores, carbon dioxide shows stronger adsorption affinity to shale organic matter because of its intermolecular interaction with the pore surfaces. However, the adsorption/desorption and diffusion behaviors related to real shale gas productions are not well understood in such nano-confined spaces. In this work, we use Molecular Dynamics and Monte Carlo methods to simulate the nano-scale methane/carbon dioxide/kerogen system. Based on a novel molecular simulation workflow, the kerogen pores are constructed using Molecular Dynamics under typical reservoir condition (up to 20 MPa), as illustrated in Fig.1, the geometry and scale of shale pore could be controlled based on cutter atoms method. As the adsorbed state and free state of fluid could be identified directly by density statistical results, the adsorption isothermal is calculated using the grand canonical Monte Carlo simulation results in a single kerogen pore and fitted using competitive adsorption models. The adsorption patterns of methane, carbon dioxide or the mixture of methane and CO2 in the organic pores could be visualized to reveal the adsorption mechanism in shale kerogen, which is shown in Fig.2. As reported in previous research, carbon dioxide shows different adsorption affinity in nano-pores, such as clay split pores, under different reservior pressure and temperature compared with methane. The competitive isothermal adsorption curve of methane and carbon dioxide in shale organic pores is illustrated in Fig.3, which shows that carbon dioxide takes larger proportion of adsorption state fluid than methane under different reservoir pressure and temperature. The mixed fluid of carbon dioxide and methane saturated in shale organic nano-pores shows different physical behavior, carbon dioxide is mainly distributed on the surface of solid and methane occupy the center of nanppores, which also means that carbon dioxide shows stronger adsorption affinity than methane. Replacement ratio is defined to represent how much methane is replaced by carbon dioxide,