Modified LB Model for Simulation of Gas Flow in Shale Pore Systems by Introducing End Effects and Local Effective Mean Free Path
Xinchi Jia,Lingjie Yu,Zhenmeng Sun,Ming Fan,Jian Cao,Xiancai Lu
DOI: https://doi.org/10.1016/j.petrol.2022.110285
2022-01-01
Abstract:Natural gas flow in shale pore systems determines the accumulation and production of shale gas. Under the conditions of reservoir, gas flow in such shale pores is significantly different from that in conventional reservoirs. Studies of gas flow under such conditions are usually limited to simple pore models with tube/slit geometries or bundle of tubes/slits, which constructs important theoretical basis, however, cannot represent the pore networks of shales. For directly simulating gas flow in complex media, we proposed a modified microscale lattice Boltzmann (MM-LB) model, in which the local effective mean free path (MFP) at any location of the pore system with complex geometry corrected by a custom two-dimensional (2D) eight-direction wall function was considered to capture the effective relaxation time for each lattice node in LB model. What's more, to consider slip velocity at solid boundaries, a combined bounce-back specular-reflection (BSR) scheme was adopted. The MM-LB model was firstly validated in 2D pore systems with simple slit/channel geometries and complex geometries (square and triangular cylinder flows) and in a 3D Sierpinski carpet, where good agreements with linearized Boltzmann solutions, molecular dynamics (MD) simulation and the direct simulation BGK (DSBGK) method results were found. Using the MM-LB model, we quantified the end effect on gas flow through pore throat and found that the end effect is caused by not only the streaming bending but also the variance of MFP (i.e., viscosity) throughout the pore throat. Finally, we applied this model to simulate shale gas flow in three digital reconstructions of kerogen pore systems and compared the results of apparent permeability prediction with Klinkenberg model and Beskok-Karniadakis (B–K) model. It is showed that the apparent permeability of gas flow in shales decreases with the decrease of temperature or the increase of pressure, and the pressure has a greater effect than temperature. For similar pore structures, the apparent permeability increases with the width of pore body and pore throat, especially the pore throat, which dominates the overall flow velocity of the entire flow field. According to the comparison and analysis of the results of MM-LB model with Klinkenberg and B–K model, it can be inferred that permeability prediction based on models of tube/slit (or bundle of tube/slit) pores misestimate the apparent permeability due to ignoring the connectivity and the end effect, especially for the conditions of lower pressures and nanoscale pore spaces.