Abstract
A method to judge shale gas flow regimes based on digital core analysis is proposed in this work. Firstly, three-dimensional shale digital cores in an anonymous shale formation in the Sichuan Basin are reconstructed by a Markov Chain Monte Carlo (MCMC) algorithm based on two-dimensional Scanning Electron Microscope (SEM) images. Then a voxel-based method is proposed to calculate the characteristic length of the three-dimensional shale digital core. The Knudsen number for three-dimensional shale digital cores is calculated by the ratio of the molecular mean free path to the characteristic length and is used to judge the flow regimes under different reservoir conditions. The results indicate that shale gas flow regimes are mainly located at the slip flow and transition flow region. Furthermore, adsorption has no obvious influence on the free gas flow regimes. Because adsorption only exists in organic pores, three-dimensional inorganic pores and organic pores in the Haynesville shale formation are reconstructed by a MCMC algorithm based on two-dimensional SEM images. The characteristic lengths of the three-dimensional inorganic pores and three-dimensional organic pores are both calculated and gas flow regimes in organic pores and inorganic pores are judged.
1 Introduction
With the rapid decline in conventional reserves, unconventional resources such as tight gas and shale gas reservoirs play an increasingly important role in the North American energy industry and have gradually become a key component of the world energy supply [1, 2, 3]. In shale gas reservoirs, gas is stored in organic pores and inorganic pores with nm-scale pore size and extremely low matrix intrinsic permeability [4, 5, 6, 7]. Zou et al. [8] reported that pore sizes in shale gas reservoirs range from 1 to 200 nm and that shale matrix intrinsic permeability ranges from 10−9 to 10−3 × 10−15 m2 [9, 10]. Gas flows in the forms of free gas and adsorbed gas in an organic pore system, while gas flows in the form of free gas in an inorganic pore system [11, 12, 13]. The size of shale nanopores approach the molecular mean free path. Therefore, the continuity assumption becomes invalid [14, 15, 16]. The Knudsen number is defined as the ratio of the molecular mean free path to the pore radius [17, 18]. At different Knudsen numbers, the corresponding gas flow regime in Figure 1 differs.
Lin et al. [19] reviewed recent advances on understanding gas flow processes in unconventional porous rocks. Javadpour [20] described gas flows in the forms of Knudsen diffusion and slip flow regime and proposed a model for gas flow in a nanopore duct. Darabi and Javadpour et al. [21] incorporated Knudsen diffusion and surface roughness into the gas flow model by Maxwell theory. However, the Javadpour model [20] includes one empirical coefficient known as the tangential momentum accommodation coefficient (TMAC) and accurate prediction of TMAC for different shale samples is not available. Beskok and Karniadakis [22] developed a unified Hagen–Poiseuille-type model valid in all flow regimes and this model has been adopted by Civan et al. [18, 23, 24] to consider the intrinsic permeability, porosity, and tortuosity of porous media. Freeman et al. [25] claimed that shale gas flows in the transitional flow regime. Landry et al. [26] deemed that gas flow in nanopores at reservoir pressure – temperature conditions falls within the slip flow and early transition flow regime. Wu et al. [27] considered shale gas flows in the forms of continuum flow, slip flow and transition flow. Sun et al. [28] concluded that Knudsen diffusion can be ignored when pressure is larger than 1 MPa. Although there are different viewpoints on shale gas flow in transition flow regime, most of the current studies suggest that shale gas flows in the slip flow regime at typical reservoir conditions [3, 29, 30].
The estimated ultimate recovery (EUR) of gas reservoirs is highly dependent on the pore structure characteristics [31, 32, 33, 34, 35, 36]. Characteristics of pore systems in Sichuan Basin shale formation in China have been studied in literature [37, 38, 39, 40]. However, the current analysis and results are mostly based on two-dimensional SEM images and the three dimensional pore structure is not fully understood. Though three dimensional organic matter and inorganic matter can be observed by Focused Ion Beam-Scanning Electron Microscopes (FIB-SEM) [41, 42, 43, 44], three dimensional organic pores and inorganic pores cannot be directly imaged and divided due to the limited resolution. On the other hand, the MCMC method is widely used for image processing research [45]. The MCMC method was applied to generate pore space of real heterogeneous porous media [46] and was also used to analyze shale pore structures [47] in our previous study.
Most of the current studies on gas flow in shale needs a predefined gas flow regime based on the Knudsen number to select the corresponding gas flow model. The Knudsen number for three-dimensional porous media is generally defined by the mean free path divided by the characteristic length of porous media. Previous work established a characteristic length analytical solution based on the two dimensional conceptual porous media [48, 49]. However, the characteristic length has not been calculated in the realistic subsurface porous media. With the recent advances in imaging techniques, detailed pore structure can be captured based on digital images, for example SEM and FIB-SEM. In this study we propose an algorithm to accurately calculate the characteristic length based on a digital image and Knudsen number for three-dimensional porous media to determine gas flow regimes and flow models. Three-dimensional shale digital cores in an anonymous shale formation in Sichuan Basin are first reconstructed by the MCMC algorithm based on two-dimensional SEM images. Then a voxel-based method is proposed to calculate the characteristic length of the shale digital core. The Knudsen number for the three-dimensional shale digital core is calculated by the ratio of the mean free path to the characteristic length and is used to judge the flow regimes under different reservoir conditions. Because adsorption only exists in organic pores, three-dimensional inorganic pores and three-dimensional organic pores are reconstructed respectively by the MCMC algorithm based on two-dimensional SEM images. Characteristic lengths of three-dimensional inorganic pores and three-dimensional organic pores are both calculated and free gas flow regimes in organic pores and inorganic pores are judged. Shale gas permeability at different reservoir conditions is also analyzed.
2 MCMC digital core reconstruction method
A Markov chain can be used for describing systems that follow a chain of linked events. Subsequent events depend only on the current state of the system. Suppose X(t) is a known random process describing the state of a process at time t0. If the state of random process X(t) in the future time t(t > t0) only depends on the state at time t0 and is uncorrelated with the state before time t0, this random process X(t) is called a Markov process. During a Markov process, the relationship between “PAST” and “FUTURE” can be described as “FUTURE” linked to the “PAST” through “NOW”, and if “NOW” is confirmed, “FUTURE” is irrelevant to the “PAST”. Suppose a system is composed of n voxels (i = 1,…,n), and X = X1,…,Xn indicates the state of the voxels. Therefore, Xi means that site i has a state xi. Specifically, for a certain voxel s, Λ_s represents all of the other points except s. The neighbor of voxel s, Ns must exist and can be described as:
Assume VLMN = {(l,m,n):0<l≤ L,0<m≤ M,0<n≤ N} represents L rows, M columns, and N layers of a rectangular grid filled with cube voxels. (i,j,k) represents the intersected voxel of row i, column j and layer k, and its associated state is expressed as Xijk. Vijk represents a rectangular parallelepiped array of voxels, and its associated state is expressed as vector X(Vijk). Any (i,j,k)∈ VLMN has the joint probability function:
The conditional probability of each voxel for its random Markov field is:
In the absence of three-dimensional information, three mutually perpendicular independent two-dimensional images can be used to reconstruct a three-dimensional chain. The core slice data from three mutually perpendicular planes xy, yz and xz are extracted and the corresponding binary images are obtained. The double voxel method is applied to simultaneously produce two new voxels, i.e. voxel (i, j, k) and (i, j+1, k), during the modeling process. Therefore, the reconstruction steps for the three-dimensional digital core are as follows. The porosity of the horizontal core slices is used as the conditional probability and the voxels on the first line of the first layer are simulated along the y-direction. The second voxel of the first line is simulated using the 2-neighborhood model, and after that, the 3-neighborhood model is used at the beginning of the simulation of the third voxel. Its conditional probability is calculated by the two-dimensional core slice on the xy plane (Figure 3(b)). The first layer is simulated along the x direction line by line. For double elements (i,j) and (i,j+1), 3 and 4 neighborhoods are used when simulating the edge voxels, and 5 and 6 neighborhoods are used when simulating the internal voxels. The conditional probability is derived by the two-dimensional core slice on the xy plane (Figure 3(c)). Each layer of voxels is simulated along the Z-direction and the three-dimensional model is reconstructed (Figure 3(d)).
3 Reconstruction of Sichuan Basin shale digital core
Three-dimensional shale digital cores in Figure 5 and Figure 6 are reconstructed based on binary SEM images (Figure 4) using the above mentioned MCMC algorithm. The blue colour represents pore phase and the red color represents matrix phase. In three-dimensional shale digital cores, pore phase is represented by number 0 and matrix is represented by number 1. Detailed parameters of shale digital cores are given in Table 1.
Core number | Geometrical size | Voxel size | Resolution | Porosity |
---|---|---|---|---|
1 | 1096nm×1096nm×1096nm | 100×100×100 | 10.96nm | 0.23 |
2 | 448nm×448nm×448nm | 100×100×100 | 4.48nm | 0.15 |
3 | 282nm×282nm×282nm | 100×100×100 | 2.82nm | 0.12 |
4 | 282nm×282nm×282nm | 100×100×100 | 2.82nm | 0.16 |
5 | 448nm×448nm×448nm | 100×100×100 | 4.48nm | 0.07 |
4 Voxel based characteristic length calculation method
For single tubes, the analysis on gas flow regimes has been well studied [20, 50]. The Knudsen number for a single tube and a three-dimensional porous media can be written respectively as:
lc is referred to as the characteristic length of porous media and can be calculated based on the chord length [48]. The chord lengths are those discontinuous lines distributed in pore phase when a line is used to cut porous media. The porous media characteristic length indicates the degree of gas diffusion in porous media and its physical meaning is the free path when gas transports within the porous media. The porous media characteristic length has been used to analyze gas transport ability in loose packings of spheres [51, 52] and can be calculated by [48, 49]:
Chord lengths are the distributions of lengths between intersections of lines with the interface (Figure 7). A chord is then a special line segment with its end points on the interface and all other points in one of the two phases. Here chord lengths in pore phase are calculated. p(z) is the chord length distribution function in pore phases. p(z)dz is the probability of finding a chord of length between z and z+dz in pore phase.
Previous studies have not applied Equation (5) to real rock samples. Here we propose a method to calculate the characteristic length based on digital image. The chord length in the binary digital image can be represented as the length of continuous number 0 in the pore phase. A regional example in two dimensions is shown in Figure 8, the voxel-based chord length is represented by the number of continuous 0 on red line. In three-dimensional pore phase, the first chord length is measured in the bottom slice on the XY plane and then the measuring window moves up in the Z direction layer by layer. Finally, chord length for three-dimensional pore phase can be obtained and Equation (8) is used to calculate the characteristic length of the three-dimensional pore phase. A validation of the proposed voxel based characteristic length calculation method can be found in the next section.
Gas adsorption is generally assumed to be a monolayer adsorption [53] and can be characterized with the Langmuir isotherm. Considering the influence of an adsorption layer, first adsorbed gas coverage on the pore wall is calculated by Equation (9) and the effective characteristic length of the three-dimensional pore phase is given in Equation (10).
According to Equation (5), Equation (6), Equation (9) and Equation (10), the Knudsen number without adsorption and with adsorption can be given respectively as:
5 Validation of the proposed voxel based characteristic length calculation method
The three-dimensional microscale lattice Boltzmann (LB) model with the D3Q19 discrete velocity model developed in our previous study [54] is adopted to simulate gas flow in a shale digital core. The developed LB model can be applied to simulate gas flow in slip flow and transition flow regimes. The basic evolution equation with the Bhatnagar-GrossKrook (BGK) collision approximation is shown as follows:
Where fα is the density distribution function of α direction; α = 0, 1, 2, … , 18; rs is the spatial location of the particles; eα is the velocity of α direction; t is time; δt is time step; τ is the relaxation time; Fα is the force term;
Where cs is the lattice sound speed; wα is the weight factor of α direction. The force term in Equation (13) can be obtained by Hermite expansion. For the D3Q19 model, it can be expressed as [55]:
Where a is the acceleration of the force. For gas flow in shale nanopores, the relaxation time should be determined by the Knudsen number. Considering the microscale effect and the effect of the Knudsen layer, the relaxation time can be expressed as [56]:
ψ(Kn) is the modification function and equals 1/(1+2Kn); N is the number of lattices occupied by the characteristic length. As the solid boundaries in shale are rough, the diffuse reflection boundary condition is very appropriate for such boundaries. The discrete format of the diffuse reflection boundary condition is expressed as [57]:
Where n is the inward unit normal vector; uw is the wall velocity; the subscript w represents the solid walls;
Detailed parameters are given in Table 2. The studied shale digital core is shown in Figure 10 and pressure drop distribution on the shale digital core is shown in Figure 11. The calculated permeability k1 is 2.7831×10−7 μm2.
Properties | Value |
---|---|
Pressure Pg/MPa | 30 |
Temperature/K | 358 |
Voxel size | 100×100×100 |
Porosity | 0.14 |
Tortuosity | 2.2 |
Resolution γ/m | 2×10−9 |
Pressure gradient/MPa/m | 0.1 |
Pressure drop on the shale digital core/Pa | 0.02 |
The Civan et al. model [18] in Equation (19) is applied to calculate gas permeability in slip flow regime and transition flow regime. The characteristic length is calculated based on the proposed voxel based method in Equation (8) and the Knudsen number is calculated using Equation (11).
The flow condition function f(Kn) is given by [22]:
The parameter α in Equation (20) is a dimensionless rarefaction coefficient and can be written as:
Though the slip coefficient β = −1 is initially considered only to be applicable to slip flow condition, evidence from the DSMC simulations and Boltzmann solutions [58] showed that β = −1 is valid within the full range of flow regimes. The calculated characteristic length in Equation (8) is 11.72nm and the calculated permeability k2 in Equation (19) is 2.9784×10−7 μm2. k2 is very close to k1 which suggests the proposed characteristic length calculation algorithm is correct.
6 Gas transport mechanism judgement and permeability analysis
6.1 Gas flow regimes judgement of Sichuan Basin shale
Firstly characteristic lengths are calculated by Equation (8) based on three-dimensional pore phase of shale digital cores shown in Figure 6. From the calculated results listed in Table 3, the average characteristic length for Sichuan Basin shale is 14.87 nm. Subsequently, the average characteristic length is used in Equation (11) and Equation (12) to study the change in the Knudsen number under different formation pressure and temperature without adsorption effect and with adsorption effect. The relative degree of deviation is defined to study the free gas transport differences in the above two cases:
Core number | Characteristic lengths /nm |
---|---|
1 | 25.2127 |
2 | 11.9197 |
3 | 10.5176 |
4 | 11.2995 |
5 | 15.3914 |
Average | 14.87 |
The Knudsen number values in Figure 12, Figure 13 are compared with the Knudsen number values in Figure 1. It is can be seen that shale gas flow regimes mainly locate at the slip flow and transition flow region. The value of the relative degree of deviation is larger in high pore pressure than that in low pore pressure. However, the maximum value of the relative degree of deviation is 5%. Therefore, the adsorption effect has no obvious influence on the free gas transport mechanism in the Sichuan Basin shale gas reservoir.
6.2 Gas flow regimes judgement based on three dimensional organic pores system and inorganic pores system
Based on the aforementioned MCMC method, SEM images (Figure 15) in Haynesville shale [59] are used to reconstruct three-dimensional organic pores and inorganic pores. First we identify organic pores and inorganic pores by the pixel value of red and green colour. Then, the corresponding binary images can be obtained using the maximum class separation distance method proposed by Otsu [60]. Taking a pixel as the basic unit, the pixel value is set to 1 when the pixel is located in the matrix phase and the pixel value is set to 0 when the pixel is located in the pore phase (Figure 16, Figure 17). Finally, the MCMC method is applied to reconstruct three-dimensional organic pores and inorganic pores based on binary images (Figure 18). The voxel sizes of the three-dimensional organic pores and inorganic pores are 800 × 800 × 800.
Based on the statement that gas transports in the forms of free gas and adsorbed gas in organic pores and gas transports only in the form of free gas in inorganic pores, Equation (8) is applied to calculate the characteristic length without adsorption effect in three-dimensional inorganic pores and organic pores. Our calculation results indicate characteristic lengths without adsorption effect in three-dimensional inorganic pores and organic pores are 49.5 nm and 35 nm respectively. Then, the value of the characteristic length without adsorption effect in three-dimensional organic pores of 35 nm is revised by Equation (10) to consider the influence of the adsorption layer. Finally, Equation (11) and Equation (12) are used respectively to calculate the Knudsen number in three-dimensional inorganic pores and organic pores. In Figure 19 and Figure 20, the Knudsen number differs a little in low pore pressure (< 5 MPa). But for shale gas reservoirs, pore pressure usually ranges from 10 MPa to 40 MPa [61] and the Knudsen number is almost the same in three-dimensional organic pores and inorganic pores under this condition. Therefore, free gas flow regimes can be deemed as the same in inorganic pores and organic pores. This may indicate that in macro scale gas flow numerical simulation, inorganic pores and organic pores can be viewed as a continuous system and a dual medium (matrix-micro fracture) mathematical model is suitable to simulate the shale gas production process.
6.3 Shale gas permeability at different reservoir conditions
According to Equation (10), Equation (12) and Equation (19), free gas permeability in consideration of adsorption can be given as:
Surface diffusion takes place within the adsorbed gas to enhance the transport of gas molecules along molecular concentration gradients [62]. Based on Hwang and Kammermeyer’s model [63], combined with methane adsorption experimental data, the surface diffusion coefficient for methane Ds0 when gas coverage is zero can be expressed as:
The isosteric adsorption heat ΔH is a function of gas coverage. According to Equation (25), the isosteric adsorption heat and gas coverage have a linear relationship and can be given as [64]:
The surface diffusion coefficient in Equation (24) is obtained under a low pressure condition by theory and experiments, and is a function of gas molecular weight, temperature, and gas activation energy, isosteric adsorption heat and is independent of pressure [63]. In order to describe the gas surface diffusion in nanopores of shale gas reservoirs under a high pressure condition, the influence of gas coverage on surface diffusion is considered. Chen et al. [65] used the kinetic method to calculate the surface diffusion coefficient:
Where H(1-κ) is Heaviside function. The surface diffusion coefficient at different pressure and temperature conditions is calculated according to Equation (26). Surface diffusion of adsorbed gas molecules can be modelled as the general diffusion process, using the molar flow rate per unit area of the concentration gradient within the adsorbed monolayer as developed in [62]:
Ca is calculated assuming Langmuir adsorption and is given by:
Camax can be expressed as [66]:
Combining Equation (28) and Equation (29), molar flow rate in the adsorbed layer is then expressed below:
From Equation (31), volumetric flow rate VA is:
According to Equation (26) and the Darcy law, adsorbed gas permeability can be written as:
Shale gas permeability in Equation (34) can be derived by combining free gas permeability in consideration of adsorption by Equation (23) with surface diffusion of adsorbed gas by Equation (33). Typical shale gas reservoir parameters are given in Table 4. Figure 21 shows that gas permeability increases with the decrease of pressure and the increase of temperature. The gas permeability value ranges from 4.3685 × 10−7 μm2 to 2.2451 × 10−6 μm2 for different reservoir conditions.
Properties | Value |
---|---|
Porosity | 0.1 |
Tortuosity | 2.3 |
Characteristic length lc /m | 1.487×10−8 |
Total organic grain volume per total grain volume ϵks | 0.01 [66] |
Langmuir pressure pL /MPa | 13.789514 [66] |
Maximum adsorbed gas concentration Cmax /mol/m3 | 328.7 [66] |
Isosteric adsorption heat at zero gas coverage Δ H(0) /J/mol | 16000 [67] |
Fitting coefficients of isosteric adsorption heat γ /J/mol | −4186 [67] |
Ideal gas constant R /J/(mol⋅k) | 8.314 |
Ratio of the rate constant for blockage to the rate constant for forward migration κ, dimensionless | 0.5 [67] |
Molecular weight M /kg/mol | 0.016 |
7 Conclusions
A method to judge shale gas flow regimes based on digital core analysis is studied. A Voxel-based method is proposed to the calculate characteristic length of a three-dimensional shale digital core. The Knudsen number for MCMC reconstructed three-dimensional shale digital core is calculated by the ratio of the molecular mean free path to the characteristic length and is used to judge the gas flow regime under different reservoir conditions. The results indicate that shale gas flow regimes mainly locate at the slip flow and transition flow region. Furthermore, adsorption has no obvious influence on the free gas flow regimes. Moreover, free gas flow regimes in organic pores and inorganic pores are judged. Our analysis results show that free gas flow regimes can be deemed as the same in two different pore systems and inorganic-organic pores can be viewed as the continuous system in a macro scale gas flow simulation. Gas permeability increases with the decrease of pressure and the increase of temperature.
Acknowledgement
This project was supported by the National Natural Science Foundation of China (No. 51504276, No.51490654), the Fundamental Research Funds for the Central Universities (No.16CX05018A, No. 18CX06008A), the Major Projects of the National Science and Technology (2016ZX05061), Graduate School Innovation Program of China University of Petroleum (YCX2017019).
Conflict of Interest
Conflict of Interests: The authors declare no conflict of interest.
Nomenclature
- ϵks
total organic grain volume per total grain volume, dimensionless
- κ
ratio of the rate constant for blockage to the rate constant for forward migration
- ΔH(0)
isosteric adsorption heat at zero gas coverage (J/mol)
- η
fitting coefficients of isosteric adsorption heat (J/mol)
- γ
resolution(m)
- λ
mean free path of molecules(m)
- μ
viscosity(Pa⋅s)
- Φ
porosity, dimensionless
- ρ
gas density (kg/m3)
- τ
tortuosity, dimensionless
- θ
gas coverage on the pore wall, dimensionless
- A
cross section area of shale digital core (m2)
- Camax
maximum adsorbed gas concentration inside the total core sample (mol/m3)
- Ca
adsorbed gas concentration (mol/m3)
- Cmax
maximum adsorbed gas concentration inside the organic matter (mol/m3)
- dm
gas molecular diameter (m) 4 × 10−10
- Ds
surface diffusion coefficient (m2/s)
- Ds0
surface diffusion coefficient when gas coverage is zero (m2/s)
- dev
relative degree of deviation, dimensionless
- Ja
molar flow rate per unit area (mol/(m2⋅s))
- kB
Boltzmann constant (J/K),1.3805 × 10−23
- Kn
Knudsen number, dimensionless
- Kn_3D
Knudsen number without adsorption effect
- Kn_3D_ads
Knudsen number with adsorption effect
- k
permeability (μm2)
- lc
characteristic length (m)
- lc_3D
characteristic length of three dimensional pore phases (m)
- lc(Pg)
characteristic length of three dimensional pore phases with adsorption effect(m)
- Mg
gas molecular weight (g/mol)
- N(ijk)
neighbor of (i,j,k)
- n
voxel number
- Pg
pore pressure(MPa)
- pL
Langmuir pressure(MPa)
- p(z)
chord length distribution function in pore phase, dimensionless
- Q
gas flux (m3/s)
- r
single tube radius (m)
- s
certain voxel, dimensionless
- t
time(s)
- vg
gas velocity (m/s)
- Vijk
rectangular parallelepiped array of voxels
- VLMN
rectangular grid filled with cube voxels, dimensionless
- X(t)
Markov process
- z_vol
chord length based on voxel unit in three dimensional digital image
- z
chord length (m)
References
[1] Zhang K., Ma X., Li Y, et al., Parameter prediction of hydraulic fracture for tight reservoir based on micro-seismic and history 65 matching, Fractals, 2018: 1840009.10.1142/S0218348X18400091Search in Google Scholar
[2] Zhang K., Zhang X., Zhang L, et al., Assisted history matching for the inversion of fractures based on discrete fracture-matrix model with different combinations of inversion parameters, Computational Geosciences, 2017, 21(5-6): 1365-1383.10.1007/s10596-017-9690-8Search in Google Scholar
[3] Wang L., Wang S., Zhang R, et al., Review of multi-scale and multi-physical simulation technologies for shale and tight gas reservoirs, Journal of Natural Gas Science and Engineering, 2017, 37: 560-578.10.1016/j.jngse.2016.11.051Search in Google Scholar
[4] Singh H., Javadpour F. Langmuir slip-Langmuir sorption permeability model of shale, Fuel, 2016, 164: 28-37.10.1016/j.fuel.2015.09.073Search in Google Scholar
[5] Ma J., Couples G. D., Jiang Z, et al., Predicting the Effective Gas Flow Properties of Gas Shale, from Nano-Pores to Shale Parasequences, (12-14 August, Denver 2013, Colorado, USA), Proceedings of URTEC: Unconventional Resources Technology Conference, Society of Petroleum Engineers, 2013.10.1190/urtec2013-237Search in Google Scholar
[6] Song W., Yao J., Ma J, et al., Pore-scale numerical investigation into the impacts of the spatial and pore-size distributions of organic matter on shale gas flow and their implications on multiscale characterisation, Fuel, 2018, 216: 707-721.10.1016/j.fuel.2017.11.114Search in Google Scholar
[7] Zhu H., Ju Y., Lu W, et al., The Characteristics and Evolution of Micro-Nano Scale Pores in Shales and Coals, Journal of Nanoscience and Nanotechnology, 2017, 17(9): 6124-6138.10.1166/jnn.2017.14529Search in Google Scholar
[8] Zou C., Zhu R., Wu S, et al., Types, characteristics, genesis and prospects of conventional and unconventional hydrocarbon accumulations: Taking tight oil and tight gas in China as an instance, Acta Petrolei Sinica, 2012, 33(2): 173-187.Search in Google Scholar
[9] Loucks R. G., Reed R. M., Ruppel S. C, et al., Morphology, genesis, and distribution of nanometer-scale pores in siliceous mudstones of the Mississippian Barnett Shale, Journal of sedimentary research, 2009, 79(12): 848-861.10.2110/jsr.2009.092Search in Google Scholar
[10] Wang F. P., Reed R. M., Pore networks and fluid flowin gas shales, (4-7 October, 2009, New Orleans, Louisiana, USA), Proceedings of the SPE Annual Technical Conference and Exhibition, Society 14 Ë W. Song et al. of Petroleum Engineers, 2009.10.2118/124253-MSSearch in Google Scholar
[11] Ma J., Sanchez J. P., Wu K, et al., A pore network model for simulating non-ideal gas flow in micro- and nano-porous materials, Fuel, 2014, 116(Supplement C): 498-508.10.1016/j.fuel.2013.08.041Search in Google Scholar
[12] Ma J., Couples G. D., Zhao X, et al., Assessing Impact of Shale Gas Adsorption on Free-Gas Permeability via a Pore Network Flow Mode, (20-22 July, 2015, San Antonio, Texas, USA), Proceedings of, URTEC: Unconventional Resources Technology Conference, Society of Petroleum Engineers, 2015.10.2118/178552-MSSearch in Google Scholar
[13] Riewchotisakul S., Akkutlu I. Y. Adsorption-Enhanced Transport of Hydrocarbons in Organic Nanopores, SPE Journal, 2016, 21(07): 1960-1969.10.2118/175107-PASearch in Google Scholar
[14] Ho C.-M., Tai Y.-C. Micro-electro-mechanical-systems (MEMS) and fluid flows, Annual Review of Fluid Mechanics, 1998, 30(1): 579-612.10.1146/annurev.fluid.30.1.579Search in Google Scholar
[15] Song W., Yao J., Li Y, et al., Apparent gas permeability in an organic-rich shale reservoir, Fuel, 2016, 181: 973-984.10.1016/j.fuel.2016.05.011Search in Google Scholar
[16] Song W., Yao J., Ma J, et al., Assessing relative contributions of transport mechanisms and real gas properties to gas flow in nanoscale organic pores in shales by pore network modelling, International Journal of Heat and Mass Transfer, 2017, 113: 524-537.10.1016/j.ijheatmasstransfer.2017.05.109Search in Google Scholar
[17] Ziarani A. S., Aguilera R. Knudsen’s permeability correction for tight porous media, Transport in Porous Media, 2012, 91(1): 239-260.10.1007/s11242-011-9842-6Search in Google Scholar
[18] Civan F. Effective correlation of apparent gas permeability in tight porous media, Transport in Porous Media, 2010, 82(2): 375-384.10.1007/s11242-009-9432-zSearch in Google Scholar
[19] Lin D., Wang J., Yuan B, et al., Review on gas flow and recovery in unconventional porous rocks, Adv Geo-Energy Res, 2017, 1: 39-53.10.26804/ager.2017.01.04Search in Google Scholar
[20] Javadpour F., Nanopores and apparent permeability of gas flow in mudrocks (shales and siltstone), Journal of Canadian Petroleum Technology, 2009, 48(08): 16-21.10.2118/09-08-16-DASearch in Google Scholar
[21] Darabi H., Ettehad A., Javadpour F, et al., Gas flow in ultra-tight shale strata, Journal of Fluid Mechanics, 2012, 710: 641-658.10.1017/jfm.2012.424Search in Google Scholar
[22] Beskok A., Karniadakis G. E., Report: a model for flows in channels, pipes, and ducts at micro and nano scales, Microscale Thermophysical Engineering, 1999, 3(1): 43-77.10.1080/108939599199864Search in Google Scholar
[23] Civan F., Devegowda D., Sigal R. F., Critical evaluation and improvement of methods for determination of matrix permeability of shale,(1-2 September October, 2013, New Orleans, Louisiana, USA), Proceedings of the SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, 2013.10.2118/166473-MSSearch in Google Scholar
[24] Civan F., Rai C. S., Sondergeld C. H., Shale-gas permeability and diffusivity inferred by improved formulation of relevant retention and transport mechanisms, Transport in Porous Media, 2011, 86(3): 925-94410.1007/s11242-010-9665-xSearch in Google Scholar
[25] Freeman C., Moridis G., Blasingame T., A numerical study of mi50 croscale flow behavior in tight gas and shale gas reservoir systems, Transport in Porous Media, 2011, 90(1): 253.10.1007/s11242-011-9761-6Search in Google Scholar
[26] Landry C. J., Prodanović M., Eichhubl P., Direct simulation of supercritical gas flow in complex nanoporous media and prediction of apparent permeability, International Journal of Coal Ge55 ology, 2016, 159: 120-134.10.1016/j.coal.2016.03.015Search in Google Scholar
[27] Wu K., Chen Z., Li X, et al., A model for multiple transport mechanisms through nanopores of shale gas reservoirs with real gas effect–adsorption-mechanic coupling, International Journal of Heat and Mass Transfer, 2016, 93: 408-426.10.1016/j.ijheatmasstransfer.2015.10.003Search in Google Scholar
[28] Sun H., Yao J., Cao Y.-c, et al., Characterization of gas transport behaviors in shale gas and tight gas reservoirs by digital rock analysis, International Journal of Heat and Mass Transfer, 2017, 104: 227-239.10.1016/j.ijheatmasstransfer.2016.07.083Search in Google Scholar
[29] Ilgen A. G., Heath J. E., Akkutlu I. Y, et al., Shales at all scales: Exploring coupled processes in mudrocks, Earth-Science Reviews, 2017, 166: 132-152.10.1016/j.earscirev.2016.12.013Search in Google Scholar
[30] Song W., Yao J., Li Y, et al., Fractal models for gas slippage factor in porous media considering second-order slip and surface adsorption, International Journal of Heat and Mass Transfer, 2018, 118: 948-960.10.1016/j.ijheatmasstransfer.2017.11.072Search in Google Scholar
[31] Wang C., Yao J., Wu K, et al., Pore-scale characterization of carbonate rock heterogeneity and prediction of petrophysical properties, (27-30 August, 2012, Aberdeen, Scotland, UK), Proceedings of the International Symposium of the Society of Core Analysts, 2012.Search in Google Scholar
[32] Jian-Chao C., A fractal approach to low velocity non-Darcy flow in a low permeability porous medium, Chinese Physics B, 2014, 23(4): 044701.10.1088/1674-1056/23/4/044701Search in Google Scholar
[33] Yao J., Wang C., Yang Y, et al., The construction of carbonate digital rock with hybrid superposition method, Journal of Petroleum Science and Engineering, 2013, 110: 263-267.10.1016/j.petrol.2013.10.005Search in Google Scholar
[34] Cai J., Luo L., Ye R, et al., Recent advances on fractal modeling of permeability for fibrous porous media, Fractals, 2015, 23(01): 1540006.10.1142/S0218348X1540006XSearch in Google Scholar
[35] Chen D., Yao Y., Fu G, et al., A new model for predicting liquid loading in deviated gas wells, Journal of Natural Gas Science and Engineering, 2016, 34: 178-184.10.1016/j.jngse.2016.06.063Search in Google Scholar
[36] Xia Y., Cai J., Wei W, et al., A new method for calculating fractal dimensions of porous media based on pore size distribution, Fractals, 2018, 26(01): 1850006.10.1142/S0218348X18500068Search in Google Scholar
[37] Guo X., Li Y., Liu R, et al., Characteristics and controlling factors of micropore structures of the Longmaxi Shale in the Jiaoshiba area, Sichuan Basin, Natural Gas Industry B, 2014, 1(2): 165-171.10.1016/j.ngib.2014.11.007Search in Google Scholar
[38] Jiao K., Yao S., Liu C, et al., The characterization and quantitative analysis of nanopores in unconventional gas reservoirs utilizing FESEM–FIB and image processing: An example from the lower Silurian Longmaxi Shale, upper Yangtze region, China, International Journal of Coal Geology, 2014, 128: 1-11.10.1016/j.coal.2014.03.004Search in Google Scholar
[39] Liang C., Jiang Z., Zhang C, et al., The shale characteristics and shale gas exploration prospects of the Lower Silurian Longmaxi shale, Sichuan Basin, South China, Journal of Natural Gas Science and Engineering, 2014, 21: 636-648.10.1016/j.jngse.2014.09.034Search in Google Scholar
[40] Shi M., Yu B., Xue Z, et al., Pore characteristics of organic-rich shales with high thermal maturity: A case study of the Longmaxi gas shale reservoirs from well Yuye-1 in southeastern Chongqing, China, Journal of Natural Gas Science and Engineering, 2015, 26: 948-959.10.1016/j.jngse.2015.07.042Search in Google Scholar
[41] Song W., Yao J., Li Y, et al., New pore size distribution calculation model based on chord length and digital image, Journal of Natural Gas Science and Engineering, 2017, 48, 111-118.10.1016/j.jngse.2016.12.041Search in Google Scholar
[42] Milliken K. L., Rudnicki M., Awwiller D. N, et al., Organic matter–hosted pore system, Marcellus formation (Devonian), Pennsylvania, AAPG bulletin, 2013, 97(2): 177-200.10.1306/07231212048Search in Google Scholar
[43] Bai B., Elgmati M., Zhang H, et al., Rock characterization of Fayetteville shale gas plays, Fuel, 2013, 105: 645-652.10.1016/j.fuel.2012.09.043Search in Google Scholar
[44] Song W., Yao J., Ma J, et al., Numerical simulation of multiphase flow in nanoporous organic matter with application to coal and gas shale systems, Water Resources Research, 2018, 54(2): 1077-1092.10.1002/2017WR021500Search in Google Scholar
[45] Geman S., Geman D., Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images, Pattern Analysis and Machine Intelligence, IEEE Transactions on, 1984, (6): 721-741.10.1109/TPAMI.1984.4767596Search in Google Scholar
[46] Wu K., Van Dijke M. I., Couples G. D, et al., 3D stochastic modelling of heterogeneous porous media–applications to reservoir rocks, Transport in Porous Media, 2006, 65(3): 443-467.10.1007/s11242-006-0006-zSearch in Google Scholar
[47] Yang Y., Yao J., Wang C, et al., New pore space characterization method of shale matrix formation by considering organic and inorganic pores, Journal of Natural Gas Science and Engineering, 2015, 27: 496-503.10.1016/j.jngse.2015.08.017Search in Google Scholar
[48] Lu B., Torquato S., Chord-length and free-path distribution functions for many-body systems, The Journal of chemical physics, 1993, 98(8): 6472-6482.10.1063/1.464812Search in Google Scholar
[49] Torquato S., Random heterogeneous materials: microstructure and macroscopic properties[M]. Springer Science & Business Media, 2013Search in Google Scholar
[50] Javadpour F., Fisher D., Unsworth M., Nanoscale gas flow in shale gas sediments, Journal of Canadian Petroleum Technology, 2007, 46(10):55-61.10.2118/07-10-06Search in Google Scholar
[51] Zalc J. M., Reyes S. C., Iglesia E., The effects of diffusion mechanism and void structure on transport rates and tortuosity factors in complex porous structures, Chemical Engineering Science, 2004, 59(14): 2947-2960.10.1016/j.ces.2004.04.028Search in Google Scholar
[52] Levitz P., Knudsen diffusion and excitation transfer in random porous media, The Journal of Physical Chemistry, 1993, 97(15): 3813-3818.10.1021/j100117a030Search in Google Scholar
[53] Heller R., Zoback M. Adsorption of methane and carbon dioxide on gas shale and pure mineral samples, Journal of Unconventional Oil and Gas Resources, 2014, 8: 14-24.10.1016/j.juogr.2014.06.001Search in Google Scholar
[54] Zhao J., Yao J., Zhang M, et al., Study of gas flow characteristics in tight porous media with a microscale lattice boltzmann model, Scientific Reports, 2016, 6: 32393.10.1038/srep32393Search in Google Scholar
[55] Suga K., Takenaka S., Ito T, et al., Evaluation of a lattice Boltzmann method in a complex nanoflow, Physical Review E, 2010, 82(1): 016701.10.1103/PhysRevE.82.016701Search in Google Scholar
[56] Li Q., He Y., Tang G, et al., Lattice Boltzmann modeling of microchannel flows in the transition flow regime, Microfluidics and nanofluidics, 2011, 10(3): 607-618.10.1007/s10404-010-0693-1Search in Google Scholar
[57] Tang G., TaoW., He Y., Lattice Boltzmann method for gaseous microflows using kinetic theory boundary conditions, Physics of 40 Fluids, 2005, 17(5): 058101.10.1063/1.1897010Search in Google Scholar
[58] Karniadakis G., Beskok A., Aluru N., Microflows and nanoflows: fundamentals and simulation[M]. Springer Science & Business Media, 2006.Search in Google Scholar
[59] Rine J., Dorsey W., Floyd M, et al., A Comparative SEM Study of Pore Types and Porosity Distribution in High to Low Porosity Samples from Selected Gas-Shale Formations, Gulf Coast Association of Geological Societies Transactions, 2010, 60:825.Search in Google Scholar
[60] Otsu N., A threshold selection method from gray-level histograms, Automatica, 1975, 11(285-296): 23-27.10.1109/TSMC.1979.4310076Search in Google Scholar
[61] Freeman C., Moridis G., Blasingame T., A numerical study of microscale flow behavior in tight gas and shale gas reservoir systems, Transport in Porous Media, 2011, 90(1): 253-268.10.1007/s11242-011-9761-6Search in Google Scholar
[62] Cunningham R. E., Williams R., Diffusion in gases and porous media[M]. Springer, 1980.10.1007/978-1-4757-4983-0Search in Google Scholar
[63] Hwang S. T., Kammermeyer K., Surface diffusion in microporous media, The Canadian Journal of Chemical Engineering, 1966, 44(2): 82-89.10.1002/cjce.5450440206Search in Google Scholar
[64] Nodzeński A., Sorption and desorption of gases (CH4, CO2) on hard coal and active carbon at elevated pressures, Fuel, 1998, 77(11): 1243-124610.1016/S0016-2361(98)00022-2Search in Google Scholar
[65] Chen Y. D., Yang R. T., Concentration dependence of surface diffusion and zeolitic diffusion, AIChE Journal, 1991, 37(10): 1579-1582.10.1002/aic.690371015Search in Google Scholar
[66] Wasaki A., Akkutlu I. Y., Permeability of organic-rich shale, SPE Journal, 2015, 20(06): 1,384-381,396.10.2118/170830-MSSearch in Google Scholar
[67] Wu K., Li X., Wang C, et al., Model for surface diffusion of adsorbed gas in nanopores of shale gas reservoirs, Industrial & Engineering Chemistry Research, 2015, 54(12): 3225-3236.10.1021/ie504030vSearch in Google Scholar
© 2018 W. Song et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.