Next Article in Journal
Non-Classical Correlations in n-Cycle Setting
Next Article in Special Issue
Mesoscopic Simulation of the (2 + 1)-Dimensional Wave Equation with Nonlinear Damping and Source Terms Using the Lattice Boltzmann BGK Model
Previous Article in Journal
Application of Entropy Spectral Method for Streamflow Forecasting in Northwest China
Previous Article in Special Issue
Evaporation Boundary Conditions for the Linear R13 Equations Based on the Onsager Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lattice Boltzmann Model for Gas Flow through Tight Porous Media with Multiple Mechanisms

1
School of Sciences, Southwest Petroleum University, Chengdu 610500, China
2
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(2), 133; https://doi.org/10.3390/e21020133
Submission received: 14 October 2018 / Revised: 27 January 2019 / Accepted: 28 January 2019 / Published: 1 February 2019
(This article belongs to the Special Issue Thermodynamics of Non-Equilibrium Gas Flows)

Abstract

:
In the development of tight gas reservoirs, gas flow through porous media usually takes place deep underground with multiple mechanisms, including gas slippage and stress sensitivity of permeability and porosity. However, little work has been done to simultaneously incorporate these mechanisms in the lattice Boltzmann model for simulating gas flow through porous media. This paper presents a lattice Boltzmann model for gas flow through porous media with a consideration of these effects. The apparent permeability and porosity are calculated based on the intrinsic permeability, intrinsic porosity, permeability modulus, porosity sensitivity exponent, and pressure. Gas flow in a two-dimensional channel filled with a homogeneous porous medium is simulated to validate the present model. Simulation results reveal that gas slippage can enhance the flow rate in tight porous media, while stress sensitivity of permeability and porosity reduces the flow rate. The simulation results of gas flow in a porous medium with different mineral components show that the gas slippage and stress sensitivity of permeability and porosity not only affect the global velocity magnitude, but also have an effect on the flow field. In addition, gas flow in a porous medium with fractures is also investigated. It is found that the fractures along the pressure-gradient direction significantly enhance the total flow rate, while the fractures perpendicular to the pressure-gradient direction have little effect on the global permeability of the porous medium. For the porous medium without fractures, the gas-slippage effect is a major influence factor on the global permeability, especially under low pressure; for the porous medium with fractures, the stress-sensitivity effect plays a more important role in gas flow.

1. Introduction

Gas flows in tight porous media have recently attracted much attention because they are widely applied in engineering fields, such as oil-gas field development, chemical processes, new energy development, and so on. The wide applications of gas flows in tight porous media also stimulate great interest in experimental and theoretical studies [1,2,3,4,5,6,7,8,9,10,11,12,13]. Gas flow in porous media usually involves three scales: Pore scale, representative elementary volume (REV) scale, and domain scale [14]. An REV in a porous medium is the smallest volume at which the scale characteristics of the porous flow hold, and thus the flow in a porous medium can be simulated at the REV scale without detailed information of the pore structures. Compared with pore-scale models, REV-scale models have a higher computational efficiency. Therefore, recently, many REV-scale models have been proposed in the literature, such as the Darcy model, extended Darcy models (including Brinkman-extended Darcy and Forchheimer–Darcy models), and generalized model based on the generalized Navier–Stokes equations. Because the generalized Navier–Stokes equations [15] can consider the combined influences of the fluid and solid drag forces, they have received much attention recently. In general, these REV-scale models can be solved by the conventional numerical methods, such as the finite difference method, finite volume method, and so on [16].
The lattice Boltzmann method (LBM) is considered a promising method for simulating fluid flow and has been successfully applied to the flow in tight porous media [17,18,19,20,21]. According to the spatial scale in simulations, the lattice Boltzmann (LB) model for porous flow can be classified into two categories, i.e., the pore-scale LB model and REV-scale LB model, which are used to simulate the flow in porous media at the pore and REV scales, respectively. Owing to the advantage of the LBM in the simulation of fluid dynamics, some REV-scale LB models have been proposed [22,23,24]. However, the macroscopic equations derived from these REV-scale LB models are the Darcy equation and extended Darcy equations, which cannot consider the effects of all the fluid and solid drag forces. Guo and Zhao [14] presented a generalized REV-scale LB model for incompressible flow through porous media, which can recover the generalized Navier–Stokes equations [15] in the incompressible limit. Later, Guo and Zhao [25] extended the generalized REV-scale LB model [14] to thermal flow in porous media, which can be used to simulate both the velocity and temperature fields of the porous flow. Rong et al. [26] proposed an REV-scale LB model for axisymmetric thermal flow in porous media. Recently, Chen et al. [27] proposed a generalized REV-scale LB model with the Klinkenberg’s effect based on the model proposed by Guo and Zhao [14]. It was found that the Klinkenberg’s effect plays an important role in gas flow in tight porous media, especially when the reservoir pressure decreases.
However, in practice, gas flow through porous media in tight gas reservoirs usually takes place deep underground, where multiple mechanisms, including gas slippage and stress sensitivity of permeability and porosity, have proven to appear as the gas flows through porous media under gas-reservoir conditions [28,29,30]. To our knowledge, no work has been done to simultaneously incorporate these mechanisms in the LB model for simulating gas flow through porous media. Therefore, the aim of this paper is to propose an REV-scale LB model for gas flow through porous media with the consideration of these multiple mechanisms. The effects of these mechanisms on gas flow in tight porous media with (or without) fractures will be studied in detail.

2. Generalized Model for Gas Flow in Tight Porous Media

2.1. Generalized Navier–Stokes Equations

Owing to the nano/micron-sized porous media considered in the present work, the pressure difference in these porous media under gas-reservoir conditions is so small that the flow can be considered incompressible flow in simulations. Recently, the generalized Navier–Stokes equations proposed by Nithiarasu et al. [15] have been widely used to study the isotheral incompressible flow in porous media, which can be expressed as follows:
u = 0 ,
t u + ( u ) ( u ϕ ) = 1 ρ p + ν e 2 u + F ,
where u is the volume-averaged velocity, t is time, ϕ is the porosity, ρ is the fluid density, p is the pressure, ν e is the effective viscosity, F is the total body force, including the porous-medium resistance and other external forces expressed as follows:
F = ϕ v K u ϕ F ϕ K | u | u + ϕ G ,
where v is the shear viscosity that is not necessarily the same as ν e , K is the permeability, G is the external force, F ϕ is the geometric function that can be expressed as [31]:
F ϕ = 1.75 150 ϕ 3 .
Note that the pressure, p , in Equation (1b) is the fluid pressure in pores instead of the volume-averaged pressure [15]. The first and second terms on the right side of Equation (2) are the linear and nonlinear drag forces due to the presence of the porous media, respectively. Strictly speaking, the nonlinear term in Equation (2) should be considered for tight matrix and fractures. For the low-speed flow in tight matrix, the value of u 2 is usually negligibly small, and thus the effect of the nonlinear drag force can be neglected for the sake of simplification, while for the high-speed flow in fractures, the nonlinear drag force should be considered owing to the large value of u 2 [32,33]. In this work, the nonlinear drag force is considered for tight matrix and fractures.

2.2. Gas Slippage Effect

Owing to the small pore size in tight porous media, the Knudsen number, Kn = λ / r , where λ represents the mean free path of the gas and r is the characteristic pore size, is usually so large that the gas-slippage effect takes place and the apparent permeability increases. It should be noted that the apparent permeability reflects the real gas-transport capacity in porous media, which is dependent on the characteristics of the porous media and fluid. The initial intrinsic permeability, which only depends on the initial porous structures, is the measured liquid permeability when the pressure is the initial reservoir pressure. The relationship between the apparent permeability and initial intrinsic permeability can be given as follows:
K a = K i f c ,
where K a is the apparent permeability, K i is the initial intrinsic permeability, f c = f p f s is the total correction factor, f p is the correction factor due to the stress-sensitivity porous structures, f s is the correction factor due to the gas-slippage effect.
Note that the porous media are located deep underground, and thus the pore structures are sensitive to the pore pressure (or effective stress), and the permeability and porosity are functions of the pore pressure (or effective stress).
A detailed description of the stress sensitivity of permeability and porosity will be presented in the next section. Here, the stress-dependent measured liquid permeability, K d , can be expressed as follows:
K d = K i f p .
Note that K d only depends on the stress-dependent porous structures, and is not affected by the fluid properties.
Similarly, the measured gas permeability without the stress sensitivity of permeability and porosity is given as follows:
K s = K i f s ,
where K s is the Klinkenberg’s corrected permeability without the effect of the stress-dependent porous structures.
To consider the influence of the gas-slippage effect, the correction factor, f s , is introduced to correct the permeability. Currently, various expressions for f s are proposed to describe the gas-slippage effect. Klinkenberg [29] proposed a widely used expression as follows:
f s = 1 + b k p ,
where b k is the Klinkenberg’s slippage factor that usually depends on the porous structures. Nevertheless, the Klinkenberg’s correlation is only a first-order correlation, which may not be suitable for the tight porous media with a large Knudsen number.
Tang et al. [34] proposed a second-order correlation, which can be considered a direct extension of the Klinkenberg’s correlation. The correction factor, f s , presented by Tang et al. [34] is written as:
f s = 1 + A p + B p 2 ,
where A and B are the slippage factors that depend on the mean free path, characteristic pore size, pressure, and so on.
In order to study the gas transport in all flow regimes, i.e., continuum flow ( Kn 0.001 ), slip flow ( 0.001 < Kn 0.1 ), transition flow ( 0.1 < Kn 10 ), and free molecular flow ( Kn > 10 ), some scholars recommended the correction factor, f s , as follows [35]:
f s = [ 1 + β ( Kn ) Kn ] [ 1 + 4 Kn 1 b s Kn ] ,
where b s is the slip coefficient and is equal to 1 for slip flow, β ( Kn ) is the rarefaction coefficient, which is given as [36]:
β ( Kn ) = 1.358 1 + 0.170 Kn 0.4348 .
Equations (9) and (10) are derived from flows in a single pipe at micro and nano scales and have been widely applied to tight porous media. According to the definition of Kn , the mean free path, λ , and the characteristic pore size, r , should be determined to obtain the value of Kn . Based on the gas kinetic theory, the mean free path, λ , can be given by [37]:
λ = μ ρ π 2 R T ,
where μ = ρ ν is the dynamic viscosity, R is the gas constant, T is the temperature. Following the suggestion of Ziarani and Aguilera [38], the characteristic pore size, r , can be evaluated by [39]:
r = 8.886 × 10 2 ( K d ϕ d ) 0.5 ,
where ϕ d is the stress-dependent porosity, and the units of r and K d are micron ( μ m ) and millidarcy ( mD ), respectively.

2.3. Stress Sensitivity of Permeability and Porosity

The porous media located deep underground usually suffer from overburden pressure due to the presence of overlying rocks. In the development of oil and gas reservoirs, the pore pressure decreases as the fluid in the porous media is exploited. Therefore, the effective stress, which is approximately equal to the overburden pressure minus the pore pressure, will increase. The increase of the effective stress will result in the deformation of the reservoir rocks and then reduce the permeability and porosity. This phenomenon is called stress sensitivity of permeability and porosity. Many studies [28,30] have pointed out that these stress-sensitivity characteristics in low permeable reservoirs are especially obvious, and thus the effect of the stress sensitivity on the permeability and porosity should be considered in the simulation of gas flow in tight porous media.
The relationship between the permeability and the pore pressure is usually described by [40]:
K d = K i exp [ γ ( p i p ) ] ,
where p i is the initial reservoir pressure, γ is the permeability modulus.
The permeability varies with the porosity, which follows the following expression [40]:
K d K i = ( ϕ d ϕ i ) α ,
where α is the porosity sensitivity exponent, ϕ i is the porosity under the initial reservoir pressure.
The pore structures in real porous media are so complex that it is difficult to obtain a generalized porosity-permeability relationship. Ergun [31] proposed an expression to describe the relationship between the permeability, K d , and porosity, ϕ d :
K d = ϕ d 3 d p 2 150 ( 1 ϕ d ) 2 ,
where d p is the solid-particle diameter. Equation (15) can be used to simulate gas flows in porous media with an approximately uniform solid-particle diameter. Note that owing to the stress sensitivity of the reservoir rocks, the solid-particle diameter, d p , is not a constant, but a function of the pore pressure, p . However, it is difficult to determine the relationship between d p and p , and there is little work to study this. In this work, we only use Equation (15) to approximatively describe the relationship between the permeability, K i , and porosity, ϕ i , under the initial reservoir pressure, p i , and the relationship between the permeability, K d , and porosity, ϕd, at arbitrary reservoir pressure, p , is defined by Equation (14). The solid-particle diameter under the initial reservoir pressure is set to be d pi . Therefore, in simulations, stress sensitivity could only affect the permeability and porosity, but not affect the solid-particle diameter and the structures of the porous media.

2.4. Generalized Navier–Stokes Equations for Tight Porous Media

To study gas flow through tight porous media in gas reservoirs, Equations (2) and (3) should be modified to consider the effects of the gas slippage and stress-sensitivity pore structures:
F = ϕ d v K a u ϕ d F ϕ K a | u | u + ϕ d G ,
F ϕ = 1.75 150 ϕ d 3 .
In summary, Equations (1), (16) and (17), together with the apparent permeability and porosity calculated by Equations (4), (9), (10), (13)–(15), compose the generalized Navier–Stokes equations for gas flow in tight porous media with multiple mechanisms.

3. Lattice Boltzmann Model

There has been considerable debate in the literature about the appropriateness of the LB model with single relaxation time for simulating Darcy flows in granular porous domains and permeability estimates [41], and the LB model with multiple relaxation times is superior to the LB model with a single relaxation time in stability and accuracy. However, the LB model with a single relaxation time has a higher computational efficiency than the LB model with multiple relaxation times and can be used to accurately simulate the porous flow at the REV scale [14,25,26,27]. Therefore, to simulate the gas flow in tight porous media with multiple mechanisms, we develop an LB model with a single relaxation time based on the work of Guo and Zhao [14] for gas flow through tight porous media with multiple mechanisms, which is given as follows:
f i ( r + c i δ t , t + δ t ) f i ( r , t ) = 1 τ [ f i ( r , t ) f i ( eq ) ( r , t ) ] + δ t J i ,
where:
f i ( eq ) = ω i ρ [ 1 + c i u c s 2 + ( c i u ) 2 2 ϕ d c s 4 u u 2 ϕ d c s 2 ] ,
J i = ( 1 1 2 τ ) ρ ω i [ c i c s 2 + ( c i u ) c i c s 2 u ϕ d c s 4 ] F ,
f i is the density distribution function associated with the discrete velocity, c i , at the site, r , and time t , f i ( eq ) is the equilibrium distribution function, δ t is the time step, τ is the dimensionless relaxation time. For the two-dimensional nine-velocity (D2Q9) model, the weight coefficients, ωi, are ω 0 = 4 / 9 , ω 1 4 = 1 / 9 , ω 5 8 = 1 / 36 ; c s = c / 3 is the lattice sound speed; c i is the discrete velocity, which is given by:
c i = { ( 0 , 0 ) , i = 0 , ( cos [ ( i 1 ) π / 2 ] , sin [ ( i 1 ) π / 2 ] ) c , i = 1 4 , ( cos [ ( 2 i 9 ) π / 4 ] , sin [ ( 2 i 9 ) π / 4 ] ) 2 c , i = 5 8 ,
where c = 3 R T .
The density and velocity are defined, respectively:
ρ = i f i ,
ρ u = i c i f i + δ t 2 ρ F .
Through the Chapman–Enskog expansion, we can derive the generalized Navier–Stokes equations presented above in the limit of a small Mach number. The pressure, p , and effective viscosity, ν e , are defined, respectively:
p = ρ c s 2 ,
ν e = c s 2 ( τ 1 2 ) δ t .
It should be noted that F includes the velocity, u , and thus Equation (22b) is a nonlinear equation for the velocity, u . According to the suggestion proposed by Guo and Zhao [14], the velocity, u, can be obtained by:
u = υ c 0 + c 0 2 + c 1 | υ | ,
where υ is a temporal velocity given as:
ρ υ = i c i f i + δ t 2 ρ ϕ d G ,
and the parameters, c 0 and c 1 , are obtained by:
c 0 = 1 2 ( 1 + ϕ d δ t 2 v K a ) ,
c 1 = ϕ d δ t 2 F ϕ K a .
The present LB model, including Equations (18)–(27), can be used to study the gas flow in the low permeable reservoirs at the REV scale with the effects of the gas slippage and stress-sensitivity pore structures.

4. Simulation Results and Discussions

4.1. Flow in a Channel Filled with a Porous Medium

First, gas flow in a channel filled with a homogeneous porous medium is simulated to validate the present LB model. The flow is driven by pressure difference, Δ p , at the inlet and outlet, which can be used to calculate the external force by G = Δ p / ( ρ L ) with L representing the length of the channel. The streamwise velocity, u x (i.e., x -direction velocity), can be described by the following equation [14]:
v e ϕ d 2 u x y 2 + G v K a u x F ϕ K a u x 2 = 0 .
with u x ( x , 0 ) = u x ( x , H ) = 0 , where H is the height of the channel. Owing to the extremely low velocity in tight porous media without fractures, the nonlinear inertial effect (i.e., the fourth term on the left-hand side of Equation (28)) can be neglected. Therefore, Equation (28) is simplified as:
v e ϕ d 2 u x y 2 + G v K a u x = 0 .
The analytical solution of Equation (29) is as follows:
u x = G K a v ( 1 cosh [ A ( y H / 2 ) ] cosh ( A H / 2 ) ) ,
where A = v ϕ d / K a v e . Note that if some fractures exist in the tight porous media, the nonlinear inertial effect needs to be considered in simulations because of the large flow rate in the fractures. In the following simulations, the effective viscosity, ν e , is assumed to equal the shear viscosity, v . The initial pressure, p i , porosity sensitivity exponent, α , initial solid-particle diameter, d pi , temperature, T, and gas constant, R , are set as 30   MPa , 3, 50 nm, 333.15 K, and 519.6545 J/(K·kg), respectively. The dynamic viscosity, μ, is calculated by the method proposed by Dempsey [42].
Figure 1 shows the effect of the permeability modulus on Ka/Ki under different pressures. It can be seen that Ka/Ki first decreases as the pressure decreases from 30 MPa, and then Ka/Ki increases as the pressure continues to decrease, which results from the combined influences of the gas slippage and stress-sensitivity pore structures. With the decrease of the pressure, the gas-slippage effect makes the permeability increase, while the stress-sensitivity effect makes the permeability decrease. The stress-sensitivity effect plays a more important role in gas transport with the increase of the permeability modulus. For a fixed value of the permeability modulus, the gas-slippage effect becomes gradually significant as the pressure decreases.
Figure 2 shows the effect of the initial porosity on K a / K i under different pressures. It is clear that the gas-slippage effect becomes more obvious with the decrease of the initial porosity. For a fixed value of the initial porosity, the gas-slippage effect increases with the decrease of the pressure.
Figure 3 shows the effect of the permeability modulus on ϕ d / ϕ i under different pressures. For a given pressure, the value of ϕ d / ϕ i decreases with the increase of the permeability modulus. When γ = 0 , the stress-sensitivity effect is neglected, and thus the porosity remains constant under different pressures. When the value of γ is small (e.g., γ 0.06   MPa 1 ), the porosity decreases linearly as the pressure decreases. However, with the increase of the value of γ , a nonlinear relationship between the porosity and pressure will appear.
In the simulations, the present LB model in lattice units is implemented, and then the simulation results are converted from lattice units to physical units [43]. Both the length and height of the channel are set as 1   μ m (i.e., L × H = 1   μ m × 1   μ m ), the initial porosity, ϕ i , is fixed at 0.3, and the pressure gradient between the inlet and outlet is set to be 1.0 MPa / m . The lattice size is chosen as 200 × 200 . The dimensionless relaxation time, τ , is set as 1.0. Periodic boundary conditions are applied to the inlet and outlet, and the nonequilibrium extrapolation conditions [44] are applied to the top and bottom walls. The initial velocities ( u x and u y ) are set to be 0, and the initial density, ρ , is set as 1.0. The convergence criterion is given as follows:
i , j { [ u x ( i , j , t + 100 δ t ) u x ( i , j , t ) ] 2 + [ u y ( i , j , t + 100 δ t ) u y ( i , j , t ) ] 2 } i , j [ u x ( i , j , t + 100 δ t ) 2 + u y ( i , j , t + 100 δ t ) 2 ] < 10 6 .
The simulation results are shown in Figure 4 and Figure 5, from which it is seen that the simulation results by the present LB model are in excellent agreement with the analytic solutions in all cases with different permeability moduli and pore pressures. Therefore, the present LB model has the capability to simulate gas flow through tight porous media with multiple mechanisms, including gas slippage and stress sensitivity of permeability and porosity.

4.2. Flow in Tight Porous Media without Fractures

In practice, the real porous media are not homogeneous systems with uniform permeability and porosity, but heterogeneous systems that consist of different mineral components. Taking shale matrix as an example, shale matrix mainly consists of organic material and inorganic material. In general, the pores in shale matrix are nanoscale pores, and thus gas transport in both organic and inorganic materials is extremely slow. Compared with the inorganic material, organic material has a relatively lower permeability and porosity. Therefore, it is essential to study gas flow in porous media with different mineral components. In this section, the structure of a heterogeneous porous medium with two mineral components is considered, which is shown in Figure 6. The randomly distributed organic material is generated by the quartet structure generation set (QSGS) method [45]. The volume fractions of the inorganic material and organic material are 0.58 and 0.42, respectively. The domain size of the porous medium is L × H = 2   μ m × 2   μ m . The initial porosities of inorganic material and organic material are 0.3 and 0.1, respectively. The permeability moduli of inorganic material and organic material are 0.06   MPa 1 and 0.12   MPa 1 , respectively. Periodic boundary conditions are applied to all the boundaries. The pressure gradient between the inlet and outlet is set to be 1.0   MPa / m . The dimensionless relaxation time, τ , is set as 1.0. The initial conditions and convergence criterion are the same as that in Section 4.1. The values of parameters in physical units and lattice units for simulations of gas flow in tight porous media without fractures are listed in Table 1 and Table 2, respectively.
The grid dependency for the simulation results is tested based on the structure of the porous medium shown in Figure 6. The two different grid sizes, namely 200 × 200 and, 400 × 400 are employed to simulate the gas flow. The effect of the grid number on the velocity magnitude distributions in the porous medium is shown in Figure 7. The velocity magnitude, | u | , is defined as | u | = u x 2 + u y 2 . There is no obvious difference between the velocity magnitude distributions based on the two grid system. The global permeability is introduced to quantify the seepage capability of the porous medium. Note that the global permeability defined by Chen et al. [27] is dependent on the flow field in the entire domain rather than the local characteristic of a porous medium. The definition of the global permeability is given as follows:
K g = v 0 H | u ( L , y ) | d y H G ,
where K g is the global permeability. It is seen from Figure 8 that no obvious change of the global permeability is observed when the grid is refined, and thus the 200 × 200 lattice is sufficient to obtain grid-independent results.
Figure 9 shows the effects of the gas slippage and stress-sensitivity pore structures on velocity magnitude distributions in the porous medium. The gas-slippage effect enhances the global velocity magnitude, while the stress-sensitivity effect reduces the global velocity magnitude as shown in Figure 9. It is seen that under a fixed pressure, the velocity magnitude distributions are different for the cases with and without the gas-slippage and stress-sensitivity effects, which is because these influencing factors have different effects on gas transport in different materials. Figure 10 shows the effects of the gas slippage and stress-sensitivity pore structures on the velocity magnitude profiles at the outlet of the porous medium. It is clear that the velocity magnitude profiles at the outlet are also different for the cases with and without the gas-slippage and stress-sensitivity effects. Therefore, multiple mechanisms, including gas slippage and stress sensitivity of permeability and porosity, affect not only the global velocity magnitude, but also the flow field. To obtain a correct prediction of gas flow in porous media, these mechanisms should be considered in simulations.
We will focus on the effects of the multiple mechanisms on the global permeability. Figure 11 shows the effects of the gas slippage and stress-sensitivity pore structures on global permeability under different pressures. It is observed that the gas slippage and stress-sensitivity pore structures significantly affect the global permeability. The gas-slippage effect increases the global permeability and becomes stronger with a decrease of the pressure. The stress-sensitivity effect decreases the global permeability as the pressure is reduced. In general, the global permeability of a porous medium deep underground is affected by the combined influences of the gas slippage and stress sensitivity. The quantitative comparison between the global permeability with and without the effects of the gas slippage and stress-sensitivity pore structures is listed in Table 3, from which it can be seen that compared with the stress-sensitivity effect, the gas-slippage effect becomes more significant in matrix, especially under low pressure.

4.3. Flow in Tight Porous Media with Fractures

In general, natural fractures may exist in a porous medium. Compared with the matrix with organic and inorganic materials, natural fractures can significantly enhance the gas-transport capability of a porous medium. In this section, gas flow in tight porous media with fractures will be studied. The effects of the fracture distributions, gas slippage, and stress sensitivity on gas transport are analyzed in detail. The distributions of the organic and inorganic materials in a porous medium are the same as that in Figure 6. The fractures with different lengths and directions are embedded in the porous medium as shown in Figure 12. The domain size of the porous medium is 20   μ m × 20   μ m , and the width of the fracture is set as 1   μ m . Generally, gas flow in fractures can be approximately considered as the seepage in a porous medium with extremely high porosity and permeability. The initial porosity and permeability modulus of the fractures are 0.99 and 0.03   MPa 1 , respectively. Periodic boundary conditions are applied to all the boundaries. The pressure gradient between the inlet and outlet is set to be 1.0   MPa / m . The dimensionless relaxation time, τ , is set as 1.0. The initial conditions and convergence criterion are the same as that in Section 4.1. The values of parameters in physical units and lattice units for simulations of gas flow in tight porous media with fractures are listed in Table 4.
Figure 13 shows the effects of the fracture distributions on velocity magnitude distributions in a porous medium. It is seen that the fracture distributions have an important effect on the flow field. When the fracture is perpendicular to the pressure-gradient direction, there is a small change of the flow field compared with a porous medium without fractures. However, when the fracture is along the pressure-gradient direction, a large change of the flow field will appear. Furthermore, the average velocity in the fracture along the pressure-gradient direction increases observably compared with that in the fracture perpendicular to the pressure-gradient direction. The average velocity in the fracture along the pressure-gradient direction increases as the length of the fracture increases. The comparison between the global permeability of a porous medium with and without fractures under the pressure of 0.5   MPa is shown in Table 5. It is observed that the fracture along the pressure-gradient direction enhances the global permeability more significantly compared with the fracture perpendicular to the pressure-gradient direction. For the fracture along the pressure-gradient direction, the penetration ratio of the fracture in the porous medium is also an important factor on the global permeability. It is found that the global permeability with completely penetrated fracture is obviously higher than that with partly penetrated fracture, and the global permeability increases when the penetration ratio increases.
The effects of the gas slippage and stress-sensitivity pore structures on global permeability of a porous medium with fractures are studied. The quantitative results are presented in Table 6, from which it is seen that the gas-slippage effect has little effect on the global permeability, while the stress-sensitivity effect significantly affects it. According to the comparison between Table 3 and Table 6, it is found that for the porous medium without fractures, the gas-slippage effect is a major influence factor on the global permeability, especially under low pressure; for the porous medium with fractures, the stress-sensitivity effect plays a more important role in gas flow.

5. Conclusions

In this paper, the generalized Navier–Stokes equations were proposed for simulating gas flow through porous media with the effects of gas slippage and stress-sensitivity porous structures. The apparent permeability and porosity were obtained based on the intrinsic permeability, intrinsic porosity, permeability modulus, porosity sensitivity exponent, and pressure. An LB model was developed to solve the generalized Navier–Stokes equations. Gas flow in a two-dimensional channel filled with a homogeneous porous medium was simulated to validate the present LB model. The simulation of gas flow in a porous medium with different mineral components was carried out. The simulation results showed that the gas-slippage effect increases the global permeability of the porous medium and becomes stronger with the decrease of the pressure, while the stress-sensitivity effect decreases the global permeability as the pressure decreases. Furthermore, it was found that these effect factors affect not only the global permeability, but also the flow field. Gas flow in a porous medium with fractures was also investigated. It was found that the fractures along the pressure-gradient direction significantly enhance the total flow rate of the porous medium, while the fractures perpendicular to the pressure-gradient direction have little effect on the global permeability. For the porous medium without fractures, the gas-slippage effect is a major influence factor on the global permeability, especially under low pressure; for the porous medium with fractures, the stress-sensitivity plays a more important role in gas flow.

Author Contributions

Data curation, Q.Z.; Formal analysis, P.G.; Funding acquisition, J.R.; Investigation, J.R.; Methodology, J.R.; Project administration, J.R.; Visualization, Q.Z.; Writing—original draft, J.R.; Writing—review & editing, C.Z.

Funding

We acknowledge the support provided by Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (PLN201627), the Natural Science Project of Sichuan Province Department of Education (No. 17ZA0425), the Scientific Research Starting Project of SWPU (No. 2017QHZ031).

Acknowledgments

Junjie Ren would like to thank Zhaoli Guo (Huazhong University of Science and Technology) for his help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, X.; Liu, L.; Wu, Y. The Uniqueness of Positive Solution for a Fractional Order Model of Turbulent Flow in a Porous Medium. Appl. Math. Lett. 2014, 37, 26–33. [Google Scholar] [CrossRef]
  2. Liu, W.; Cui, J.; Xin, J. A Block-centered Finite Difference Method for an Unsteady Asymptotic Coupled Model in Fractured Media Aquifer System. J. Comput. Appl. Math. 2018, 337, 319–340. [Google Scholar] [CrossRef]
  3. Ren, J.; Guo, P. A Novel Semi-analytical Model for Finite-conductivity Multiple Fractured Horizontal Wells in Shale Gas Reservoirs. J. Nat. Gas. Sci. Eng. 2015, 24, 35–51. [Google Scholar] [CrossRef]
  4. Ren, J.; Guo, P. Anomalous Diffusion Performance of Multiple Fractured Horizontal Wells in Shale Gas Reservoirs. J. Nat. Gas. Sci. Eng. 2015, 26, 642–651. [Google Scholar] [CrossRef]
  5. Ren, J.; Guo, P. Performance of Vertical Fractured Wells with Multiple Finite-conductivity Fractures. J. Geophys. Eng. 2015, 12, 978–987. [Google Scholar] [CrossRef]
  6. Ren, J.; Guo, P. Nonlinear Flow Model of Multiple Fractured Horizontal Wells with Stimulated Reservoir Volume including the Quadratic Gradient Term. J. Hydrol. 2017, 554, 155–172. [Google Scholar] [CrossRef]
  7. Ren, J.; Guo, P. Nonlinear Seepage Model for Multiple Fractured Horizontal Wells with the Effect of the Quadratic Gradient Term. J. Porous Media 2018, 21, 223–239. [Google Scholar] [CrossRef]
  8. Ren, J.; Guo, P. A General Analytical Method for Transient Flow Rate with the Stress-sensitive Effect. J. Hydrol. 2018, 565, 262–275. [Google Scholar] [CrossRef]
  9. Ren, J.; Guo, P. Analytical Method for Transient Flow Rate with the Effect of the Quadratic Gradient Term. J. Pet. Sci. Eng. 2018, 162, 774–784. [Google Scholar] [CrossRef]
  10. Ren, J.; Guo, P.; Peng, S.; Ma, Z. Performance of Multi-stage Fractured Horizontal Wells with Stimulated Reservoir Volume in Tight Gas Reservoirs Considering Anomalous Diffusion. Environ. Earth Sci. 2018, 77, 768. [Google Scholar] [CrossRef]
  11. Wang, Z.; Tu, H.; Guo, P.; Zeng, F.; Sang, T.; Wang, Z. Fluid Behavior of Gas Condensate System with Water Vapor. Fluid Phase Equilib. 2017, 438, 67–75. [Google Scholar] [CrossRef]
  12. Wang, Z.; Li, Y.; Liu, H.; Zeng, F.; Guo, P.; Jiang, W. Study on the Adsorption, Diffusion and Permeation Selectivity of Shale Gas in Organics. Energies 2017, 10, 142. [Google Scholar] [CrossRef]
  13. Wang, Z.; Wang, Z.; Zeng, F.; Guo, P.; Xu, X.; Deng, D. The Material Balance Equation for Fractured Vuggy Gas Reservoirs with Bottom Water-drive Combining Stress and Gravity Effects. J. Nat. Gas. Sci. Eng. 2017, 44, 96–108. [Google Scholar] [CrossRef]
  14. Guo, Z.L.; Zhao, T.S. Lattice Boltzmann Model for Incompressible Flows through Porous Media. Phys. Rev. E 2002, 66, 036304. [Google Scholar] [CrossRef] [PubMed]
  15. Nithiarasu, P.; Seetharamu, K.N.; Sundararajan, T. Natural Convective Heat Transfer in a Fluid Saturated Variable Porosity Medium. Int. J. Heat Mass Transf. 1997, 40, 3955–3967. [Google Scholar] [CrossRef]
  16. Ferziger, J.H.; Perić, M. Computational Methods for Fluid Dynamics. Phys. Today 1997, 50, 80–84. [Google Scholar] [CrossRef]
  17. Ren, J.; Guo, P.; Guo, Z.; Wang, Z. A Lattice Boltzmann Model for simulating Gas Flow in Kerogen Pores. Transp. Porous Med. 2015, 106, 285–301. [Google Scholar] [CrossRef]
  18. Ren, J.; Guo, P.; Peng, S.; Yang, C. Investigation on Permeability of Shale Matrix using the Lattice Boltzmann Method. J. Nat. Gas Sci. Eng. 2016, 29, 169–175. [Google Scholar] [CrossRef]
  19. Ren, J.; Guo, P.; Guo, Z. Rectangular Lattice Boltzmann Equation for Gaseous Microscale Flow. Adv. Appl. Math. Mech. 2016, 8, 306–330. [Google Scholar] [CrossRef]
  20. Ren, J.; Guo, P. Lattice Boltzmann Simulation of Steady Flow in a Semi-Elliptical Cavity. Commun. Comput. Phys. 2017, 21, 692–717. [Google Scholar] [CrossRef]
  21. Ren, J.; Zheng, Q.; Guo, P.; Peng, S.; Wang, Z.; Du, J. Pore-scale Lattice Boltzmann Simulation of Two-component Shale Gas Flow. J. Nat. Gas. Sci. Eng. 2019, 61, 46–70. [Google Scholar] [CrossRef]
  22. Spaid, M.A.A.; Phelan, F.R. Lattice Boltzmann Methods for Modeling Microscale Flow in Fibrous Porous Media. Phys. Fluids 1997, 9, 2468–2474. [Google Scholar] [CrossRef]
  23. Dardis, O.; McCloskey, J. Lattice Boltzmann Scheme with Real Numbered Solid Density for the Simulation of Flow in Porous Media. Phys. Rev. E 1998, 57, 4834–4837. [Google Scholar] [CrossRef]
  24. Martys, N.S. Improved Approximation of the Brinkman Equation using a Lattice Boltzmann Method. Phys. Fluids 2001, 13, 1807–1810. [Google Scholar] [CrossRef]
  25. Guo, Z.L.; Zhao, T.S. A Lattice Boltzmann Model for Convection Heat Transfer in Porous Media. Numer. Heat Transf. B Fundam. 2005, 47, 157–177. [Google Scholar] [CrossRef]
  26. Rong, F.M.; Guo, Z.L.; Chai, Z.H.; Shi, B.C. A Lattice Boltzmann Model for Axisymmetric Thermal Flows through Porous Media. Int. J. Heat Mass Transf. 2010, 53, 5519–5527. [Google Scholar] [CrossRef]
  27. Chen, L.; Fang, W.Z.; Kang, Q.J.; Hyman, J.D.; Viswanathan, H.S.; Tao, W.Q. Generalized Lattice Boltzmann Model for Flow through Tight Porous Media with Klinkenberg’s Effect. Phys. Rev. E 2015, 91, 033004. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, H.-X.; Wang, G.; Chen, Z.-X.; Wong, R.C.K. Deformational Characteristics of Rock in Low Permeable Reservoir and their Effect on Permeability. J. Pet. Sci. Eng. 2010, 75, 240–243. [Google Scholar] [CrossRef]
  29. Klinkenberg, L.J. The Permeability of Porous Media to Liquids and Gases. API Drill. Prod. Pract. 1941, 2, 200–213. [Google Scholar] [CrossRef]
  30. Vairogs, J.; Hearn, C.L.; Dareing, D.; Rhoades, V.W. Effect of Rock Stress on Gas Production from Low-permeability Reservoirs. J. Pet. Technol. 1971, 23, 1161–1167. [Google Scholar] [CrossRef]
  31. Ergun, S. Fluid Flow through Packed Columns. Chem. Eng. Prog. 1952, 48, 89–94. [Google Scholar]
  32. Luo, W.; Wang, X.; Feng, Y.; Tang, C.; Zhou, Y. Productivity Analysis for a Vertically Fractured Well under Non-Darcy Flow Condition. J. Pet. Sci. Eng. 2016, 146, 714–725. [Google Scholar] [CrossRef]
  33. Guppy, K.H.; Cinco-Ley, H.; Ramey, H.J., Jr. Pressure Buildup Analysis of Fractured Wells Producing at High Flow Rates. J. Pet. Technol. 1982, 34, 2656–2666. [Google Scholar] [CrossRef]
  34. Tang, G.H.; Tao, W.Q.; He, Y.L. Gas Slippage Effect on Microscale Porous Flow using the Lattice Boltzmann Method. Phys. Rev. E 2005, 72, 056301. [Google Scholar] [CrossRef] [PubMed]
  35. Beskok, A.; Karniadakis, G.E. A Model for Flows in Channels, Pipes, and Ducts at Micro and Nano Scales. Microsc. Therm. Eng. 1999, 3, 43–77. [Google Scholar]
  36. Civan, F. Effective Correlation of Apparent Gas Permeability in Tight Porous Media. Transp. Porous Med. 2010, 82, 375–384. [Google Scholar] [CrossRef]
  37. Cercignani, C. Mathematical Methods in Kinetic Theory; Plenum Press: New York, NY, USA, 1990. [Google Scholar]
  38. Ziarani, A.S.; Aguilera, R. Knudsen’s Permeability Correction for Tight Porous Media. Transp. Porous Med. 2012, 91, 239–260. [Google Scholar] [CrossRef]
  39. Heid, J.G.; McMahon, J.J.; Nielsen, R.F.; Yuster, S.T. Study of the Permeability of Rocks to Homogeneous Fluids; American Petroleum Institute: Washington, DC, USA, 1950. [Google Scholar]
  40. David, C.; Wong, T.-F.; Zhu, W.L.; Zhang, J.X. Laboratory Measurement of Compaction-induced Permeability Change in Porous Rocks: Implications for the Generation and Maintenance of Pore Pressure Excess in the Crust. Pure Appl. Geophys. 1994, 143, 425–456. [Google Scholar] [CrossRef]
  41. Prestininzi, P.; Montessori, A.; La Rocca, M.; Succi, S. Reassessing the Single Relaxation Time Lattice Boltzmann Method for the Simulation of Darcy’s Flows. Int. J. Mod. Phys. C 2016, 27, 1650037. [Google Scholar] [CrossRef]
  42. Dempsey, J.R. Computer Routine Treats Gas Viscosity as a Variable. Oil Gas J. 1965, 63, 141–143. [Google Scholar]
  43. Succi, S. The Lattice Boltzmann Equation for Fluid Dynamics and Beyond; Oxford University Press: New York, NY, USA, 2001. [Google Scholar]
  44. Guo, Z.L.; Zheng, C.G.; Shi, B.C. Non-equilibrium Extrapolation Method for Velocity and Boundary Conditions in the Lattice Boltzman Method. Chin. Phys. 2002, 11, 366–374. [Google Scholar]
  45. Wang, M.; Wang, J.; Pan, N.; Chen, S. Mesoscopic Predictions of the Effective Thermal Conductivity for Microscale Random Porous Media. Phys. Rev. E 2007, 75, 036702. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The effect of the permeability modulus on K a / K i under different pressures. The pressure varies from 1 to 30 MPa, and the initial porosity is set as ϕ i = 0.1 .
Figure 1. The effect of the permeability modulus on K a / K i under different pressures. The pressure varies from 1 to 30 MPa, and the initial porosity is set as ϕ i = 0.1 .
Entropy 21 00133 g001
Figure 2. The effect of the initial porosity on K a / K i under different pressures. The pressure varies from 1 to 30 MPa, and the permeability modulus is set as γ = 0.09   MPa 1 .
Figure 2. The effect of the initial porosity on K a / K i under different pressures. The pressure varies from 1 to 30 MPa, and the permeability modulus is set as γ = 0.09   MPa 1 .
Entropy 21 00133 g002
Figure 3. The effect of the permeability modulus on ϕd/ϕi under different pressures. The pressure varies from 1 to 30 MPa, and the initial porosity is set as ϕ i = 0.1 .
Figure 3. The effect of the permeability modulus on ϕd/ϕi under different pressures. The pressure varies from 1 to 30 MPa, and the initial porosity is set as ϕ i = 0.1 .
Entropy 21 00133 g003
Figure 4. Velocity profiles in a channel filled with a porous medium under different pressures. Solid lines represent the analytical solutions and symbols represent the simulation results by the present model. The permeability modulus is set as γ = 0.06   MPa 1 .
Figure 4. Velocity profiles in a channel filled with a porous medium under different pressures. Solid lines represent the analytical solutions and symbols represent the simulation results by the present model. The permeability modulus is set as γ = 0.06   MPa 1 .
Entropy 21 00133 g004
Figure 5. Velocity profiles in a channel filled with a porous medium under different permeability moduli. Solid lines represent the analytical solutions and symbols represent the simulation results by the present model. The pressure is set as p = 25   MPa .
Figure 5. Velocity profiles in a channel filled with a porous medium under different permeability moduli. Solid lines represent the analytical solutions and symbols represent the simulation results by the present model. The pressure is set as p = 25   MPa .
Entropy 21 00133 g005
Figure 6. Structure of a heterogeneous porous medium with two mineral components: organic material with relatively lower permeability and porosity (white), and inorganic material with relatively higher permeability and porosity (black). The randomly distributed organic material is generated by the QSGS method. Volume fractions of the inorganic material and organic material are 0.58 and 0.42, respectively.
Figure 6. Structure of a heterogeneous porous medium with two mineral components: organic material with relatively lower permeability and porosity (white), and inorganic material with relatively higher permeability and porosity (black). The randomly distributed organic material is generated by the QSGS method. Volume fractions of the inorganic material and organic material are 0.58 and 0.42, respectively.
Entropy 21 00133 g006
Figure 7. Velocity magnitude distributions in the porous medium based on different grid numbers. The pressure is set as p = 10   MPa . (a) 200 × 200; (b) 400 × 400.
Figure 7. Velocity magnitude distributions in the porous medium based on different grid numbers. The pressure is set as p = 10   MPa . (a) 200 × 200; (b) 400 × 400.
Entropy 21 00133 g007
Figure 8. Global permeability based on different grid numbers. The pressure is set as p = 10   MPa .
Figure 8. Global permeability based on different grid numbers. The pressure is set as p = 10   MPa .
Entropy 21 00133 g008
Figure 9. The effects of the gas slippage and stress-sensitivity pore structures on velocity magnitude distributions in the porous medium. (a) p = 0.5   MPa , (b) p = 10   MPa , (c) p = 25   MPa . The first row: without the effects of gas slippage and stress-sensitivity pore structures; the second row: with only the effect of gas slippage; the third row: with only the effect of stress-sensitivity pore structures; the fourth row: with gas slippage and stress-sensitivity pore structures.
Figure 9. The effects of the gas slippage and stress-sensitivity pore structures on velocity magnitude distributions in the porous medium. (a) p = 0.5   MPa , (b) p = 10   MPa , (c) p = 25   MPa . The first row: without the effects of gas slippage and stress-sensitivity pore structures; the second row: with only the effect of gas slippage; the third row: with only the effect of stress-sensitivity pore structures; the fourth row: with gas slippage and stress-sensitivity pore structures.
Entropy 21 00133 g009
Figure 10. The effects of the gas slippage and stress-sensitivity pore structures on velocity magnitude profiles at the outlet of the porous medium. (a) p = 0.5   MPa , (b) p = 10   MPa , (c) p = 25   MPa . The first row: without the effects of gas slippage and stress-sensitivity pore structures; the second row: with only the effect of gas slippage; the third row: with only the effect of stress-sensitivity pore structures; the fourth row: with gas slippage and stress-sensitivity pore structures.
Figure 10. The effects of the gas slippage and stress-sensitivity pore structures on velocity magnitude profiles at the outlet of the porous medium. (a) p = 0.5   MPa , (b) p = 10   MPa , (c) p = 25   MPa . The first row: without the effects of gas slippage and stress-sensitivity pore structures; the second row: with only the effect of gas slippage; the third row: with only the effect of stress-sensitivity pore structures; the fourth row: with gas slippage and stress-sensitivity pore structures.
Entropy 21 00133 g010
Figure 11. The effects of the gas slippage and stress-sensitivity pore structures on the global permeability under different pressures. (b) is a subfigure of (a).
Figure 11. The effects of the gas slippage and stress-sensitivity pore structures on the global permeability under different pressures. (b) is a subfigure of (a).
Entropy 21 00133 g011
Figure 12. Structures of a heterogeneous porous medium with no fracture (a), a completely penetrated fracture along the vertical centerline (b), a completely penetrated fracture along the horizontal centerline (c), and partly penetrated fractures along the horizontal centerline with the penetration ratio being 1/4 (d), 1/2 (e), 3/4 (f), 7/8 (g), 15/16 (h), and 31/32 (i). White represents the fractures, black represents the inorganic matrix, and gray represents the organic matrix.
Figure 12. Structures of a heterogeneous porous medium with no fracture (a), a completely penetrated fracture along the vertical centerline (b), a completely penetrated fracture along the horizontal centerline (c), and partly penetrated fractures along the horizontal centerline with the penetration ratio being 1/4 (d), 1/2 (e), 3/4 (f), 7/8 (g), 15/16 (h), and 31/32 (i). White represents the fractures, black represents the inorganic matrix, and gray represents the organic matrix.
Entropy 21 00133 g012
Figure 13. The effects of the fracture distributions on velocity magnitude distributions in a porous medium under the pressure of 0.5 MPa. (ai) are the simulation results based on the structures of the porous medium shown in Figure 12a–i, respectively.
Figure 13. The effects of the fracture distributions on velocity magnitude distributions in a porous medium under the pressure of 0.5 MPa. (ai) are the simulation results based on the structures of the porous medium shown in Figure 12a–i, respectively.
Entropy 21 00133 g013
Table 1. Basic parameters in physical units used for simulations of gas flow in tight porous media without fractures.
Table 1. Basic parameters in physical units used for simulations of gas flow in tight porous media without fractures.
Parameter 0.5   MPa 10   MPa 25   MPa
L   ( m ) 2 × 10 6 2 × 10 6 2 × 10 6
H   ( m ) 2 × 10 6 2 × 10 6 2 × 10 6
G   ( Pa m 2 / kg ) 3.4625 × 10 5 1.7312 × 10 4 6.9249 × 10 3
v   ( Pa s m 3 / kg ) 3.8651 × 10 6 2.4740 × 10 7 1.3805 × 10 7
Inorganic material K a   ( m 2 ) 4.1931 × 10 18 5.9209 × 10 19 9.9148 × 10 19
Inorganic material K d   ( m 2 ) 1.5643 × 10 19 2.7661 × 10 19 6.8034 × 10 19
Inorganic material K i   ( m 2 ) 9.1837 × 10 19 9.1837 × 10 19 9.1837 × 10 19
Inorganic material ϕ d 0.16630.20110.2715
Inorganic material ϕ i 0.30.30.3
Organic material K a   ( m 2 ) 1.1587 × 10 19 1.5430 × 10 20 3.4592 × 10 20
Organic material K d   ( m 2 ) 5.9698 × 10 22 1.8666 × 10 21 1.1292 × 10 20
Organic material K i   ( m 2 ) 2.0576 × 10 20 2.0576 × 10 20 2.0576 × 10 20
Organic material ϕ d 0.03070.04490.0819
Organic material ϕ i 0.10.10.1
Table 2. Basic parameters in lattice units used for simulations of gas flow in tight porous media without fractures.
Table 2. Basic parameters in lattice units used for simulations of gas flow in tight porous media without fractures.
Parameter 0.5   MPa 10   MPa 25   MPa
L   ( lattice   unit ) 2 2 2
H   ( lattice   unit ) 2 2 2
G   ( lattice   unit ) 2.3178 × 10 14 2.8284 × 10 13 3.6310 × 10 13
v   ( lattice   unit ) 1.0 × 10 6 1.0 × 10 6 1.0 × 10 6
Inorganic material K a   ( lattice   unit ) 4.1931 × 10 6 5.9209 × 10 7 9.9148 × 10 7
Inorganic material K d   ( lattice   unit ) 1.5643 × 10 7 2.7661 × 10 7 6.8034 × 10 7
Inorganic material K i   ( lattice   unit ) 9.1837 × 10 7 9.1837 × 10 7 9.1837 × 10 7
Inorganic material ϕ d 0.16630.20110.2715
Inorganic material ϕ i 0.30.30.3
Organic material K a   ( lattice   unit ) 1.1587 × 10 7 1.5430 × 10 8 3.4592 × 10 8
Organic material K d   ( lattice   unit ) 5.9698 × 10 10 1.8666 × 10 9 1.1292 × 10 8
Organic material K i   ( lattice   unit ) 2.0576 × 10 8 2.0576 × 10 8 2.0576 × 10 8
Organic material ϕ d 0.03070.04490.0819
Organic material ϕ i 0.10.10.1
Table 3. Comparison between the global permeability with and without the effects of the gas slippage and stress-sensitivity pore structures in a porous medium without fractures. The variable with the subscript, g , represents the corresponding global permeability.
Table 3. Comparison between the global permeability with and without the effects of the gas slippage and stress-sensitivity pore structures in a porous medium without fractures. The variable with the subscript, g , represents the corresponding global permeability.
Pressure K gs / K gi K gd / K gi K ga / K gi
0.5 MPa16.07340.16214.4781
10 MPa1.80650.28970.6788
25 MPa1.45220.73191.1345
Table 4. Basic parameters in physical units and lattice units used for simulations of gas flow in tight porous media with fractures under the pressure of 0.5   MPa .
Table 4. Basic parameters in physical units and lattice units used for simulations of gas flow in tight porous media with fractures under the pressure of 0.5   MPa .
Parameter Physical UnitLattice Unit
L 20 × 10 6   ( m ) 20
H 20 × 10 6 ( m ) 20
G 3.4625 × 10 5   ( Pa m 2 / kg ) 2.3178 × 10 14
v 3.8651 × 10 6   ( Pa s m 3 / kg ) 1.0 × 10 6
Inorganic material K a 4.1931 × 10 18   ( m 2 ) 4.1931 × 10 6
Inorganic material K d 1.5643 × 10 19   ( m 2 ) 1.5643 × 10 7
Inorganic material K i 9.1837 × 10 19   ( m 2 ) 9.1837 × 10 7
Inorganic material ϕ d 0.16630.1663
Inorganic material ϕ i 0.30.3
Organic material K a 1.1587 × 10 19   ( m 2 ) 1.1587 × 10 7
Organic material K d 5.9698 × 10 22   ( m 2 ) 5.9698 × 10 10
Organic material K i 2.0576 × 10 20   ( m 2 ) 2.0576 × 10 8
Organic material ϕ d 0.03070.0307
Organic material ϕ i 0.10.1
Natural fracture K a 7.0996 × 10 14   ( m 2 ) 7.0996 × 10 2
Natural fracture K d 6.6743 × 10 14   ( m 2 ) 6.6743 × 10 2
Natural fracture K i 1.6172 × 10 13   ( m 2 ) 1.6172 × 10 1
Natural fracture ϕ d 0.73710.7371
Natural fracture ϕ i 0.990.99
Table 5. Comparison between the global permeability of a porous medium with and without the fractures under the pressure of 0.5   MPa .
Table 5. Comparison between the global permeability of a porous medium with and without the fractures under the pressure of 0.5   MPa .
Structures of a Porous Medium with FractureThe Ratio of Global Permeability with Fracture and Global Permeability without Fracture
Figure 12b1.4037
Figure 12c938.2724
Figure 12d2.6933
Figure 12e5.0843
Figure 12f7.6190
Figure 12g11.1979
Figure 12h13.4764
Figure 12i18.1392
Table 6. Comparison between the global permeability with and without the effects of the gas slippage and stress-sensitivity pore structures in a porous medium with fractures. The structure of the porous medium is shown in Figure 12c. The variable with the subscript, g , represents the corresponding global permeability.
Table 6. Comparison between the global permeability with and without the effects of the gas slippage and stress-sensitivity pore structures in a porous medium with fractures. The structure of the porous medium is shown in Figure 12c. The variable with the subscript, g , represents the corresponding global permeability.
Pressure K gs / K gi K gd / K gi K ga / K gi
0.5 MPa1.02260.55940.5797
10 MPa1.00140.68070.6820
25 MPa1.00080.91160.9123

Share and Cite

MDPI and ACS Style

Ren, J.; Zheng, Q.; Guo, P.; Zhao, C. Lattice Boltzmann Model for Gas Flow through Tight Porous Media with Multiple Mechanisms. Entropy 2019, 21, 133. https://doi.org/10.3390/e21020133

AMA Style

Ren J, Zheng Q, Guo P, Zhao C. Lattice Boltzmann Model for Gas Flow through Tight Porous Media with Multiple Mechanisms. Entropy. 2019; 21(2):133. https://doi.org/10.3390/e21020133

Chicago/Turabian Style

Ren, Junjie, Qiao Zheng, Ping Guo, and Chunlan Zhao. 2019. "Lattice Boltzmann Model for Gas Flow through Tight Porous Media with Multiple Mechanisms" Entropy 21, no. 2: 133. https://doi.org/10.3390/e21020133

APA Style

Ren, J., Zheng, Q., Guo, P., & Zhao, C. (2019). Lattice Boltzmann Model for Gas Flow through Tight Porous Media with Multiple Mechanisms. Entropy, 21(2), 133. https://doi.org/10.3390/e21020133

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop