Molecular Simulation of Methane Adsorption Behavior in Kerogen Nanopores for Shale Gas Resource Assessment
- Jinrong Cao (The University of Tokyo) | Yunfeng Liang (The University of Tokyo) | Yoshihiro Masuda (The University of Tokyo) | Hiroaki Koga (Japan Oil, Gas and Metals National Corporation) | Hiroyuki Tanaka (Japan Oil, Gas and Metals National Corporation) | Kohei Tamura (Japan Oil, Gas and Metals National Corporation) | Sunao Takagi (Japan Oil, Gas and Metals National Corporation) | Toshifumi Matsuoka (Fukada Geological Institute)
- Document ID
- International Petroleum Technology Conference
- International Petroleum Technology Conference, 26-28 March, Beijing, China
- Publication Date
- Document Type
- Conference Paper
- 2019. International Petroleum Technology Conference
- 4.6 Natural Gas, 5.1 Reservoir Characterisation, 5.8 Unconventional and Complex Reservoirs, 5.5 Reservoir Simulation, 5.1.1 Exploration, Development, Structural Geology, 5.1 Reservoir Characterisation, 4.6 Natural Gas, 1.2.3 Rock properties, 5.8.2 Shale Gas, 5 Reservoir Desciption & Dynamics
- Molecular simulation, Methane adsorption, Kerogen nanopores, Shale gas
- 9 in the last 30 days
- 359 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 7.00|
|SPE Non-Member Price:||USD 23.00|
In this paper, we present an improved method to predict the methane adsorption isotherm for a real shale sample using molecular dynamics (MD) simulation with a realistic kerogen model. We compare our simulation results both to the experiment and to the simulation results on the basis of a simple graphite model, and show how our procedure leads to the creation of more accurate adsorption isotherms of a shale sample at a wide range of pressure. A Marcellus shale sample was chosen as an example to demonstrate how to calculate the adsorption isotherms using MD simulations. Type II kerogen molecular model was selected for the dry gas window. The constructed bulk kerogen model contains mesopores (> 2 nm) and micropores (≤ 2 nm) inside. Ten different mesopore sizes of kerogen nanopore systems were constructed. According to the characteristics of methane density distribution in the simulation system, three regions can be clearly distinguished, free gas, adsorbed gas, and absorbed gas. We show that the adsorbed gas per unit pore volume increases with the pore size decreased. This is similar to previous molecular simulations with graphite model. For predicting the total adsorption isotherm of a real shale sample, both adsorbed and absorbed gas were considered. For the adsorption amount, the calculated adsorption isotherms were averaged based on pore size distribution of that Marcellus Shale sample. For nanopores smaller than 5 nm, we used total organic carbon (TOC) data to weight the absorption contribution in the kerogen bulk (i.e. inside the micropores). The total adsorption isotherm thus obtained from our simulations reproduced experiments very well. Importantly, kerogen model has overcome the difficulties of prediction using graphite models (i.e. an underestimation of adsorption under high pressure conditions) as documented in previous studies. Furthermore, we predicted the adsorption isotherms for higher temperatures. With the temperature increased, lower adsorption amount is predicted. The novelty of our improved method is that it is able to predict methane adsorption isotherm at a wide range of pressure for a shale sample by considering both adsorption in kerogen mesopores and absorption in kerogen bulk. It can be readily used for any shale sample, where the pore size distribution, porosity, and TOC are known. We remark that the above results and conclusion resulted from our simple assumption. Further discussion might be necessary.
|File Size||2 MB||Number of Pages||19|
Chalmers GR, Bustin RM, Power IM. Characterization of gas shale pore systems by porosimetry, pycnometry, surface area, and field emission scanning electron microscopy/transmission electron microscopy image analyses: Examples from the Barnett, Woodford, Haynesville, Marcellus, and Doig Units. AAPG Bull. 2012; 96: 1099–1119.
Charpentier RR, de Witt W, Claypool, GE, Harris LD, Mast RF, Megeath JD, Roen JB, Schmoker JW. Estimates of unconverntional natural gas resources of Devonian shales of the Appalachian basin. In Petroleum Geology of Devonian and Mississippian Black Shale of Eastern North America; Roen JB, Kepferle RC, Eds.; United States Geological Survey (USGS): Reston, VA 1993; USGS Bullet. 1909B, pp. N1–N20.
Jin L, Pu H, Wang Y, Li Y. The consideration of pore size distribution in organic-rich unconventional formations may increase oil production and reserve by 25%, Eagle Ford case study. Paper URTEC 2148314 presented at the Unconventional Resources Technology Conference, San Antonio, Texas, USA, 20-22 July 2015.
Lemmon, E. W.; McLinden, M. O.; Friend, D. G. Thermophysical Properties of Fluid Systems, NIST Chemistry WebBook: NIST Standard Reference Database Number 69; Linstrom, P. J., Mallard, W. G., Eds.; http://webbook.nist.gov/chemistry/. 2017.
Passey QR, Bohacs KM, Esch WL, Klimentidis R, Sinha S. From oil-prone source rock to gas-producing shale reservoir – Geologic and petrophysical characterization of unconventional shale-gas reservoirs. Paper SPE-131350 presented at the CPS/SPE International Oil & Gas Conference and Exhibition, Beijing, China, 8-10 June 2010.
Saidian M, Kuila U, Rivera S, Godinez LJ, Prasad M. Porosity and pore size distribution in mudrocks: A comparative study for Haynesville, Niobrara, Monterey and Eastern European Silurian formations. Paper URTEC 1922745 presented at the Unconventional Resources Technology Conference, Denver, Colorado, USA, 25-27 August 2014.
Sawa Y, Liang Y, Murata S, Matsuoka T, Akai T, Takagi S. Pore-filling nature of CH4 adsorption behavior in kerogen nano pores: A molecular view based on full-atom kerogen models. Paper SPE 176999 presented at the SPE Asia Pacific Unconventional Resources Conference, Brisbane, Australia, 9-11 November 2015.