Computational Fluid Dynamics Modeling of Slip Flow Coupled with Gas Adsorption/Desorption Kinetics in Complex Pore Space
- Shanshan Yao (University of Regina) | Ronny Pini (Imperial College London) | Xiangzeng Wang (Yanchang Petroleum Group) | Fanhua Zeng (University of Regina) | Ning Ju (University of Regina)
- Document ID
- Society of Petroleum Engineers
- SPE Annual Technical Conference and Exhibition, 24-26 September, Dallas, Texas, USA
- Publication Date
- Document Type
- Conference Paper
- 2018. Society of Petroleum Engineers
- 1.10 Drilling Equipment, 5.8.2 Shale Gas, 5.8 Unconventional and Complex Reservoirs, 5 Reservoir Desciption & Dynamics, 1.10 Drilling Equipment
- Pore Space Reconstruction, Slip Flow, CFD, Adsorption/Desorption Kinetics
- 2 in the last 30 days
- 199 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 8.50|
|SPE Non-Member Price:||USD 25.00|
At reservoir conditions, gas flow confined in submicron pores of shale falls within slip flow and transition flow regimes. Beyond the common instant equilibrium assumption, we believe that gas adsorption/desorption on rough pore surfaces could be in non-equilibrium status when gas pressure keeps decreasing during production. We investigate the interplay of gas slip flow inside complex submicron-scale pores and gas adsorption/desorption kinetics on pore surfaces with computational fluid dynamics (CFD) under unsteady-state flow conditions.
Different from previous studies, the gas adsorption/desorption is in non-equilibrium state, which is closer to real reservoir conditions. Given pore pressure Pp at time t, linear driving force model with gas desorption rate coefficient kd is applied to describe the difference between the equilibrium adsorption amount (calculated with adsorption isotherms) and the actual adsorption amount per unit pore surface area. Free gas flow inside 3D reconstructions of shale pore space is modeled by Navier-Stokes equations with Maxwell's first-order slip boundary conditions. To include gas contributions from desorption, extra source with strength equal to the gas desorption rate is added to the slip boundaries.
Any type of adsorption isotherms can be incorporated into our CFD modeling. We investigate the coupling of slip flow and Langmuir adsorption isotherms for methane in 3D reconstructed pore space. We observe that not all of adsorbed gas measured in adsorption isotherms contribute to gas production. In our study the pore pressure, Pp, decreases along with time t. One significant finding is that there exists a key time point, tk, after which adsorbed gas starts desorbing off pore surfaces and the decreasing rate of pore pressure becomes smaller. The higher the gas desorption rate coefficient, kd, is, the earlier tk occurs. But the decreasing rate of pore pressure is no longer sensitive to the coefficient, kd, when kd is larger than 0.0005. Another significant finding is that optimum gas production rate should be chosen according to kd values. It is always better to have high production rate when kd is larger than 0.0005. If the kd is smaller than 0.0005, simulations are necessary to determine optimum gas production rate by balancing the development efficiency and final amount of gas desorption.
Gas adsorption/desorption is always regarded as an instant equilibrium process in shale reservoir simulations. This study considers the non-equilibrium gas adsorption/desorption process, which is closer to real reservoir conditions. No studies in the literature have considered the influence of gas adsorption/desorption kinetics when choosing optimum production rates. CFD simulations in this study provide insight and guidelines on optimizing shale gas development with evaluating slip flow as well as gas adsorption/desorption characteristics.
|File Size||1 MB||Number of Pages||12|
Chalmers, G. R., Bustin, R. M.Power, I. M.: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 Foig units. AAPG Bulletin 56(6),1099–1119 (2012)
Gasparik, M., Rexer, T.F.T., Aplin, A.C., Billemont, P., Weireld, G. De., Gensterblum, Y., Henry, M., Krooss, B. M., Liu, S., Ma, X., Sakurovs, R., Song, Z., Staib, G., Thomas, K. M., Wang, S., Zhang, T.: First international inter-laboratory comparison of high-pressure CH4, CO2 and C2H6 sorption isotherms on carbonaceous shales. International Journal of Coal Geology 132, 131–146 (2014)
Sondergeld, C.H., Ambrose, R. J., Rai, C.S., Moncrieff, J. M.: Mirco-structural studies of gas shales. SPE. https://doi.org/10.2118/131771-MS (2010)
Yao, S., Wang, X., Yuan, Q., Zeng, F.: Estimation of shale intrinsic permeability with process-based pore network modeling approach. Transp. Porous Med (2018) 10.1007/s11242-018-1091-5