The lattice Boltzmann method (LBM) has experienced tremendous advances and been well accepted as a popular method of simulation of various fluid flow mechanisms on pore scale in tight formations. With the introduction of an effective relaxation time and slip boundary conditions, the LBM has been successfully extended to solve micro-gaseous related transport and phenomena. As gas flow in shale matrix is mostly in the slip flow and transition flow regimes, given the difficulties of experimental techniques to determine extremely low permeability, it appears that the computational methods especially the LBM can be an attractive choice for simulation of these micro-gaseous flows. In this paper an extensive overview on a number of relaxation time and boundary conditions used in LBM-like models for micro-gaseous flow are carried out and their advantages and disadvantages are discussed. Furthermore, potential application of the LBM in flow simulation in shale gas reservoirs on pore scale and representative elementary volume(REV) scale is evaluated and summarised. Our review indicates that the LBM is capable of capturing gas flow in continuum to slip flow regimes which cover significant proportion of the pores in shale gas reservoirs and identifies opportunities for future research.

a b s t r a c tGas flow in shale is associated with both organic matter (OM) and inorganic matter (IOM) which contain nano-pores ranging in size from a few to hundreds of nano-meters. In addition to the non-continuum effect which leads to an apparent permeability of gas higher than the intrinsic permeability, the surface diffusion of adsorbed gas in organic pores also can influence the apparent permeability through its own transport mechanism. In this study, a generalized lattice Boltzmann model (GLBM) is employed for gas flow through the reconstructed shale matrix consisting of OM and IOM. The ExpectationMaximization (EM) algorithm is used to assign the pore size distribution to each component, and the dusty gas model (DGM) and generalized Maxwell-Stefan model (GMS) are adopted to calculate the apparent permeability accounting for multiple transport mechanisms including viscous flow, Knudsen diffusion and surface diffusion. Effects of pore radius and pressure on permeability of both IOM and OM as well as effects of Langmuir parameters on OM are investigated. The effect of total organic content and distribution on the apparent permeability of the reconstructed shale matrix at different surface diffusivity is also studied. It is found that the influence of pore size and pressure on the apparent permeability of organic matter is affected by the surface diffusion of adsorbed gas. Moreover, surface diffusion plays a significant role in determining apparent permeability and the velocity distribution of shale matrix.

Recent studies have shown that adsorbed gas and its surface diffusion have profound influence on micro-gaseous flow through organic pores in shale gas reservoirs. In this paper, a multiple-relaxation-time (MRT) LB model is adopted to estimate the apparent permeability of organic shale and a new boundary condition, which combines Langmuir adsorption theory with Maxwellian diffusive reflection boundary condition, is proposed to capture gas slip and surface diffusion of adsorbed gas. The simulation results match well with previous studies carried out using Molecular Dynamics (MD) and show that Maxwell slip boundary condition fails to characterize gas transport in the near wall region under the influence of the adsorbed gas. The total molar flux can be either enhanced or reduced depending on variations in adsorbed gas coverage and surface diffusion velocity. The effects of pore width, pressure as well as Langmuir properties on apparent permeability of methane transport in organic pores are further studied. It is found that the surface transport plays a significant role in determining the apparent permeability, and the variation of apparent permeability with pore size and pressure is affected by the adsorption and surface diffusion.