# Nonlinear Seepage Model of Gas Transport in Multiscale Shale Gas Reservoirs and Productivity Analysis of Fractured Well

Authors: Ting Huang; Xiao Guo; Kun Wang
Journal: Journal of Chemistry
 (2015)

---

## Abstract

Shale is  abundant in nanoscale pores, so gas flow in shales cannot be simply represented by Darcy formula anymore. It is crucial to figure out the influence of gas flow in nano/micro pores on actual productivity, which can provide basic theories for optimizing parameters and improving the gas production from engineering perspective. This paper considers the effects of slippage and diffusion in nanoscale based on Beskok-Karniadakis (BK) equation, which can be applicable for different flow regimes including continuum flow, slip flow, transition flow, and free-molecule flow. A new non-Darcy equation was developed based on the analysis of effects of high order terms of BK equation on permeability correction factor. By using the conformal transformation principle and pressure coupling method, we established the productivity formula of fractured well (infinite and limited conductivity) satisfying mass variable seepage flowing in fractures. The simulation results have been compared with field data and influencing parameters are analyzed thoroughly. It is concluded that slippage effect affects gas production of fractured well when wellbore pressure is less than 15 MPa, and the effects of slippage and diffusion have greater influence on gas production of fractured well for reservoir with smaller permeability, especially when permeability is at nano-Darcy scale.


---

## Body

## 1. Introduction

Shale has been found to contain numerous well-developed nanoscale pores in recent researches: pore diameter of Haynesville basin in North America mainly distributes from 2 nm to 20 nm [1]; pore diameter of Mississippian basin is between 5 nm and 750 nm [2]; pore diameter of Sichuan basin in China is around 100 nm [3]. Core analysis shows that the permeability of shale is in a range of 1 × 10−9–1 × 10−3μm2 [4], and flow in extremely low permeability shales undergoes a transition from Darcy regime to other regimes owing to the significant effect of collisions between molecules and pore walls on gas transport. Darcy Formula cannot describe all the flow regimes in shale gas reservoirs, so it is necessary to establish a transport equation valid for all flow regimes in multiscale shale gas reservoirs.Gas flow in shales cannot be simply represented by Darcy equation anymore owing to pore throat diameters on the order of nanometers [4–11]. Javadpour et al. [5] employed a group of Gaussian distribution curves to model nonlinear flow in a combination of nano/micropore networks and was validated with the results of canister desorption tests. Wang and Reed [6] proposed that gas flow in shale matrix can be non-Darcy as a result of the slippage effect but a Darcy type in fractures. Huang et al. [7, 8] used apparent permeability which combines Knudsen diffusion coefficient and slippage factor to describe gas flow in nanoscale pores within shale matrix. Michel et al. [9] established a model to correct the permeability to account for all flow regimes including continuum, such as slip, transition, and molecular flow in shales based on Beskok and Karniadakis equation [10, 11], but the slip coefficient b is assumed to be constant. Deng et al. [4] developed a seepage model in regard to different slip coefficients on the basis of Michel’s model, and diffusion is considered as well.However, the slippage factor is only assigned to some different constants in previous articles which cannot really simulate the entire flow regimes for gas flow in shale gas reservoirs, but actually slippage factor is changing with pressure, pore radius, and other parameters. This paper considers the effects of slippage and diffusion based on Beskok-Karniadakis (BK) equation, and slippage factor and diffusion coefficient are varying with pressure and shale permeability in the process of simulation. A new non-Darcy equation of motion is obtained based on the analysis of the effects of high order terms of BK equation on permeability correction factor. By employing conformal transformation principle and pressure coupling method, a second order differential equation satisfying variable mass flowing in fractures is derived, and productivity formula of fractured well (infinite and limited conductivity) for shale gas reservoir is established. The simulation results are compared with field data, and effects of slippage and diffusion on productivity of fractured well have been analyzed thoroughly as well as other parameters.


## 2. Flow Regimes in Shale Gas Reservoirs

Knudsen number Kn is defined as the ratio of the fluid mean free pathλ and pore throat diameter r, which is a widely recognized dimensionless parameter to determine the degree of appropriateness of continuum model:(1)Kn=λ¯r.Gas flow regimes can be classified into four categories [12] depending on the Knudsen number (Table 1): (1) continuum flow, (2) slip flow, (3) transition flow, (4) and free-molecule flow. In continuum flow regime, no-slip boundary condition is valid and gas flow is linear. As Kn increases, the rarefaction effects become more pronounced and the continuum assumption breaks down eventually. So in flow regimes other than continuum flow, Darcy law is not applicable.Table 1
Classification of gas flow regimes based on Knudsen number.
Knudsen number
Kn
≤
0.001
0.001
<
Kn
≤
0.1
0.1
<
Kn
≤
10
Kn
>
10
Flow regime
Continuum flow
Slip flow
Transition flow
Free-molecule flowTo analyze the flow regimes in shale gas reservoirs, Figure1 presents the Knudsen number under different pore sizes between 1 nm and 100 μm and different pressure ranging from atmosphere to 100 MPa where gas flow from formation to surface undergoes. Knudsen number increases when pressure drops and pore throat diameter decreases. The horizontal plane in Figure 1 represents that Kn equals 0.001. Gas flow regime in shales is continuum flow when Knudsen number is beneath the horizontal plane, but it changes into other flow regimes when Knudsen number is above the horizontal plane.Figure 1
The relationship between Knudsen number, pressure, and pore radius.For example, when formation pressure is 100 MPa, gas flow in shale gas reservoirs is slip flow when pore diameter is between 1 nm and 63 nm and is transition flow when pore diameter is smaller than 1 nm, but when formation pressure drops to 10 MPa, gas flow in shale gas reservoirs is slip flow when pore diameter is between 10 nm and 631 nm and is transition flow when pore diameter is smaller than 10 nm. Gas flow in shale gas reservoirs is free-molecule flow when pore diameter is smaller than 6.3 nm and pressure is below 8 MPa. So we can conclude that gas flow in shale gas reservoirs is a multiscale flow process which includes continuum flow, slip flow, transition flow, and free-molecule flow.


## 3. Multiscale Non-Darcy Seepage Model in Shale Matrix

Volume flow rate of gas flow in reservoir can be represented by Darcy Law:(2)vg=-Kμdpdx.According to the Beskok-Karniadakis equation which accounts for different flow regimes of the continuum flow, slip flow, transition flow, and free-molecule flow, the relationship between flow rate and pressure gradient can be described as(3)vg=-K0μ1+αKn1+4Kn1-bKndpdx.Combining (2) and (3), permeability correction factor can be written as(4)K=K0ξξ=1+αKn1+4Kn1-bKn.We can see from (4) that the apparent permeability is close to the value of absolute permeability when Kn approaches 0. When Knudsen number gets bigger, which means that the flow in microtube is no longer Darcy flow, the apparent permeability should be corrected.The rarefaction coefficient is defined by Beskok and Karniadakis (1999) as(5)α=12815π2tan-14Kn0.4.To simplify Beskok-Karniadakis equation [4, 9, 13], permeability correction factor can be rewritten as(6)ξ=1+αKn+4Kn1-bKn+4αKn21-bKn.The first two terms of (6) represent the first order term to no-slip flow in Beskok-Karniadakis equation. When Kn<0.1, the following approximations are valid [4, 9]:(7)α≈12815π24Kn0.4-134Kn0.43(8)Kn1-bKn≈Kn1+bKn+b2Kn2.WhenKn>0.1, the value of rarefaction coefficient calculated by (7) is negative, which is very different from rarefaction coefficient calculated by (5), so we use another approximation instead:(9)α≈12815π24Kn0.4.By combining (8) and (9), the permeability correction factor can be rewritten as(10)ξ=1+4Kn+51215π2Kn1.4+4bKn2+204815π2Kn2.4+4b2Kn3+oKn3.Here, the influence of high order terms of (10) on permeability correction factor will be discussed. In Figures 2(a) and 2(b), the term of Kn2.4 has little influence on permeability correction factor whether the value of slippage factor b is big or small; the effect of term Kn2 on permeability correction factor can be neglected when the value of slippage factor is small, but this effect becomes greater when slippage factor gets bigger. For slip flow and continuum flow, the high order (>2) correction term of (10) can be neglected, and (3) can be derived as(11)vg=-K0μ1+4Kn+51215π2Kn1.4+4bKn2dpdx.Figure 2
Relationship between permeability correction factor and Knudsen number.
(a)
Slippage factorb=0.1
(b)
Slippage factorb=5The definition of mean path is given as [14](12)λ¯=πzRT2Mμp.Knudsen diffusion coefficient is mainly changing with pore diameter, which is defined as [5, 12, 15]:(13)DK=2r38RTπM0.5.Slippage factorb is changing with pressure and pore diameter, which is given by [16–18](14)b=8πRTM0.5μpavgr2f-1.Knudsen number can be represented by Knudsen diffusion coefficient according to (12) and (13):(15)Kn=3π8r2μpDK.Taking (15) into (11), we can get(16)vg=-K0μ1+43π8r2μpDK+51215π23π8r2μpDK1.4+4b3π8r2μpDK2dpdx.According to Hagen-Poiseuille model, relationship between permeability and pore diameter can be written as(17)K0=r28.Taking (17) into (16), we can derive(18)vg=-K0μ1+3π16K0μpDK+51215π23π64K0μpDK1.4+4b3π64K0μpDK2dpdx.Considering both sides of (18) multiplied by seepage area and divided by volume factor, we can get the volume flow rate at standard condition:(19)qv=vABg=AK0Bgμg1+3π16K0μpDK+51215π23π64K0μpDK1.4+4b3π64K0μpDK2dpdx,where the gas volume factor Bg can be derived according to the equation of state:(20)Bg=nzRT/pnRTsc/psc=pscpzTTsc.


## 4. Productivity Model of Fractured Well

The natural productivity of shale gas reservoirs is very low, so formation has to be artificially stimulated (hydraulically fractured) to render commercial production. So this paper establishes steady seepage model of fractured well for shale gas reservoirs by conformal transformation. According to the principle of conformal transformation, the gas production and pressure of the boundary remains the same, only the length of line and flow forms change [19, 20].


### 4.1. Physical Assumptions

To make this mathematical model more tractable, the following simplified assumptions are made.(1)
The reservoir is isotropic and gas flow in shale gas reservoir is planar (steady state).(2)
Flow in nanopores of shale matrix is assumed to be non-Darcy, and fluid in shale gas reservoirs is single phase.(3)
Fractured vertical well is produced at a constant wellbore pressure and is penetrated fully in shale gas reservoir. The height of hydraulic fractures is equal to the thickness of formation.(4)
The fractures are vertical and symmetrically located at two sides of wellbore. Fracture half-length isLf; fracture width is negligible when the vertical fractures are assumed to have infinite conductivity and is Wf when vertical fractures are assumed to have finite conductivity.


### 4.2. Fractured Well with Infinite Conductivity

Conformal transformation function is defined as(21)Z=LfchW.TakingZ=x+iy, W=u+iv into (21), we can get(22)x=Lfchucos⁡vy=Lfshusinv.We can see from Table2 that the upper half plane formation in Z plane (Figure 3(a)) is transformed to a half infinite formation with width of π at the right side of fractured well in W plane (Figure 3(b)); elliptical constant pressure boundary with radius of Re is transformed to a linear one with length of π; fractured well with length of 2Lf is transformed to a line sink in W plane.Table 2
Conformal transformations of the fractured well.
Location of point
W plane
Z plane
①
u
=
0
v
=
0
x
=
L
f
y
=
0
②
u
=
0
v
=
π
/
2
x
=
0
y
=
0
③
u
=
0
v
=
π
x
=
-
L
f
y
=
0
④
u
=
u
0
v
=
0
x
=
L
f
ch
u
0
y
=
0
⑤
u
=
u
0
v
=
π
/
2
x
=
0
y
=
L
f
sh
u
0
⑥
u
=
u
0
v
=
π
x
=
-
L
f
ch
u
0
y
=
0Figure 3
Sketch map of conformal transformation to fractured well.Equipotential line equation of gas flow inZ plane can be obtained through (22)(23)x2Lf2ch2u+y2Lf2sh2u=cos⁡2v+sin2v=1.The boundary can be viewed as circular when the boundary is far from fractured well, so(24)chu0≈shu0≈12eu0.Then the equipotential line equation of boundary can be rewritten as(25)x2+y2=Lf212eu02=Re2.Solving (25), we can get the relationship between the length of drainage area u0 in W plane and the radius of boundary Re in Z plane:(26)u0=ln⁡2ReLf.Combined with (19), production rate of drainage area in W plane can be derived:(27)Qv=2πK0hTscpscTμz¯u0pe2-pw22+3πμDK16K0pe-pw+5129π23πμDK64K01.4pe0.6-pw0.6+4b3πμDK64K02ln⁡pepw.Gas production of fractured well inZ plane can be obtained according to (26)(28)Qv=2πK0hTscpscTμz¯ln⁡⁡(2Re/Lf)·pe2-pw22+3πμDK16K0pe-pw+5129π23πμDK64K01.4pe0.6-pw0.6+4b3πμDK64K02ln⁡pepw.


### 4.3. Fractured Well with Finite Conductivity

Compared with hydraulic fracture with infinite conductivity, fractured vertical well with finite conductivity is more close to real condition, and the effects of parameters of hydraulic fractures on well performance can be analyzed as well.Considering an element of hydraulic fracture inW plane in Figure 3(b), we assume that gas only flows in v direction (flows to wellbore) in hydraulic fracture (Figure 4). So flow rate in fracture vv is a constant value in u direction and a variable value in v direction; flow rate from reservoir to the element vu equals a constant value.Figure 4
Sketch map of hydraulic fracture element.The mass balance equation of gas flow in hydraulic fracture can be obtained by principle of mass conservation:(29)-vvv+Δv-vvv·12Wfh+vu·Δvh=0.Considering both sides of (29) divided by Δv, the partial difference equation can be derived when Δv takes an infinitesimal:(30)dvvdv·12Wf=vu.Pseudopressure function of gas flow in hydraulic fracture is defined as(31)mp=2∫pep1+3π16K0μpDK+51215π23π64K0μpDK1.4+4b3π64K0μpDK2pμzdp.Equation (30) can be(32)∂2m∂v2-K0(1/2)kfWf1ln⁡⁡(2Re/Lf)m=-K01/2kfWf1ln⁡⁡2Re/Lfme.The boundary conditions of hydraulic fracture are given as follows:(33)dmdv=0;v=0m=mw;v=π2.Combining with (33), we can get the solution of (32): (34)mp=c1eλv+c2e-λv+me,where the parameters c1, c2, and λ are defined as follows:(35)λ=2K0kfWf1ln⁡2Re/Lfc1=c2=mw-mee(π/2)λ+e-(π/2)λ.The gas production of fractured vertical well can be written as(36)Qf=-kfWfhTscpscTdmdvv=π/2=kfWfhTscpscTλme-mweπ/2λ-e-π/2λeπ/2λ+e-π/2λ.Combined with (31), we can get the productivity formula of fractured vertical well with finite conductivity:(37)Qf=2kfWfhλTscpscTμz¯tanhπλ2·pe2-pw22+3πμDK16K0pe-pw+5129π23πμDK64K01.4pe0.6-pw0.6+4b3πμDK64K02ln⁡pepw,where the first term in the bracket is the gas production calculated by Darcy formula, which is defined as QD.


## 4.1. Physical Assumptions

To make this mathematical model more tractable, the following simplified assumptions are made.(1)
The reservoir is isotropic and gas flow in shale gas reservoir is planar (steady state).(2)
Flow in nanopores of shale matrix is assumed to be non-Darcy, and fluid in shale gas reservoirs is single phase.(3)
Fractured vertical well is produced at a constant wellbore pressure and is penetrated fully in shale gas reservoir. The height of hydraulic fractures is equal to the thickness of formation.(4)
The fractures are vertical and symmetrically located at two sides of wellbore. Fracture half-length isLf; fracture width is negligible when the vertical fractures are assumed to have infinite conductivity and is Wf when vertical fractures are assumed to have finite conductivity.


## 4.2. Fractured Well with Infinite Conductivity

Conformal transformation function is defined as(21)Z=LfchW.TakingZ=x+iy, W=u+iv into (21), we can get(22)x=Lfchucos⁡vy=Lfshusinv.We can see from Table2 that the upper half plane formation in Z plane (Figure 3(a)) is transformed to a half infinite formation with width of π at the right side of fractured well in W plane (Figure 3(b)); elliptical constant pressure boundary with radius of Re is transformed to a linear one with length of π; fractured well with length of 2Lf is transformed to a line sink in W plane.Table 2
Conformal transformations of the fractured well.
Location of point
W plane
Z plane
①
u
=
0
v
=
0
x
=
L
f
y
=
0
②
u
=
0
v
=
π
/
2
x
=
0
y
=
0
③
u
=
0
v
=
π
x
=
-
L
f
y
=
0
④
u
=
u
0
v
=
0
x
=
L
f
ch
u
0
y
=
0
⑤
u
=
u
0
v
=
π
/
2
x
=
0
y
=
L
f
sh
u
0
⑥
u
=
u
0
v
=
π
x
=
-
L
f
ch
u
0
y
=
0Figure 3
Sketch map of conformal transformation to fractured well.Equipotential line equation of gas flow inZ plane can be obtained through (22)(23)x2Lf2ch2u+y2Lf2sh2u=cos⁡2v+sin2v=1.The boundary can be viewed as circular when the boundary is far from fractured well, so(24)chu0≈shu0≈12eu0.Then the equipotential line equation of boundary can be rewritten as(25)x2+y2=Lf212eu02=Re2.Solving (25), we can get the relationship between the length of drainage area u0 in W plane and the radius of boundary Re in Z plane:(26)u0=ln⁡2ReLf.Combined with (19), production rate of drainage area in W plane can be derived:(27)Qv=2πK0hTscpscTμz¯u0pe2-pw22+3πμDK16K0pe-pw+5129π23πμDK64K01.4pe0.6-pw0.6+4b3πμDK64K02ln⁡pepw.Gas production of fractured well inZ plane can be obtained according to (26)(28)Qv=2πK0hTscpscTμz¯ln⁡⁡(2Re/Lf)·pe2-pw22+3πμDK16K0pe-pw+5129π23πμDK64K01.4pe0.6-pw0.6+4b3πμDK64K02ln⁡pepw.


## 4.3. Fractured Well with Finite Conductivity

Compared with hydraulic fracture with infinite conductivity, fractured vertical well with finite conductivity is more close to real condition, and the effects of parameters of hydraulic fractures on well performance can be analyzed as well.Considering an element of hydraulic fracture inW plane in Figure 3(b), we assume that gas only flows in v direction (flows to wellbore) in hydraulic fracture (Figure 4). So flow rate in fracture vv is a constant value in u direction and a variable value in v direction; flow rate from reservoir to the element vu equals a constant value.Figure 4
Sketch map of hydraulic fracture element.The mass balance equation of gas flow in hydraulic fracture can be obtained by principle of mass conservation:(29)-vvv+Δv-vvv·12Wfh+vu·Δvh=0.Considering both sides of (29) divided by Δv, the partial difference equation can be derived when Δv takes an infinitesimal:(30)dvvdv·12Wf=vu.Pseudopressure function of gas flow in hydraulic fracture is defined as(31)mp=2∫pep1+3π16K0μpDK+51215π23π64K0μpDK1.4+4b3π64K0μpDK2pμzdp.Equation (30) can be(32)∂2m∂v2-K0(1/2)kfWf1ln⁡⁡(2Re/Lf)m=-K01/2kfWf1ln⁡⁡2Re/Lfme.The boundary conditions of hydraulic fracture are given as follows:(33)dmdv=0;v=0m=mw;v=π2.Combining with (33), we can get the solution of (32): (34)mp=c1eλv+c2e-λv+me,where the parameters c1, c2, and λ are defined as follows:(35)λ=2K0kfWf1ln⁡2Re/Lfc1=c2=mw-mee(π/2)λ+e-(π/2)λ.The gas production of fractured vertical well can be written as(36)Qf=-kfWfhTscpscTdmdvv=π/2=kfWfhTscpscTλme-mweπ/2λ-e-π/2λeπ/2λ+e-π/2λ.Combined with (31), we can get the productivity formula of fractured vertical well with finite conductivity:(37)Qf=2kfWfhλTscpscTμz¯tanhπλ2·pe2-pw22+3πμDK16K0pe-pw+5129π23πμDK64K01.4pe0.6-pw0.6+4b3πμDK64K02ln⁡pepw,where the first term in the bracket is the gas production calculated by Darcy formula, which is defined as QD.


## 5. Results and Discussion

According to the productivity formula of fractured well with finite conductivity, the effects of permeability correction factor, slippage factor, Knudsen diffusion coefficient, matrix permeability, fracture half-length, and fracture conductivity on productivity of fractured well are analyzed by combining the data from a single well of shale gas reservoirs in China. The data for production simulation are presented in Table3.Table 3
Data of shale gas reservoir for production simulation.
Parameter
Value
Units
Permeability,K
0.005
mD
Porosity,ϕ
0.07
—
Formation temperature,T
366.15
K
Surface temperature,Tsc
293
K
Formation pressure,pe
28
MPa
Formation thickness,h
30.5
m
Pressure relief radius,Re
400
m
Gas viscosity,μ
0.014
mPa·s
Gas compressibility factor,z
0.89
—
Radius of wellbore,rw
0.1
m


### 5.1. Effect of Permeability Correction Factor on Productivity

Figure5 presents the effect of permeability correction factor on productivity of fractured well, and we can see from this figure that permeability correction factor has a great influence on productivity. Production of fractured well is very small when our model only considers Darcy flow in nano/microscale pores (ξ=1), but gas production increases a lot when matrix apparent permeability is corrected by permeability correction factor ξ. Field data shows that gas production of a fractured vertical well is 1.2 × 104 m3 when production pressure drop keeps 7 MPa and fracture half-length is 200 m [4]. The calculation results of multiscale non-Darcy seepage model are much more close to the field data compared with the production rate calculated by Darcy formula and production rate calculated by productivity formula only considering first order term of permeability correction factor.Figure 5
Effect of permeability correction factorξ on productivity.


### 5.2. Effects of Formation and Fracture Parameters on Productivity

The effect of matrix permeabilityK on productivity of fractured well is illustrated in Figure 6. Gas production of fractured well increases with the matrix permeability, and gas production increases to 1.71 times, 2.32 times, and 2.87 times when matrix permeability increases to 2 times, 3 times, and 4 times, respectively. Figure 7 reflects productivity of fractured well affected by fracture half-length. Gas production increases to 1.54 times, 2.06 times, and 2.65 times when fracture half-length increases to 3 times, 4 times, and 5 times, respectively. We can find that the increase of gas production rate affected by matrix permeability is still more than production rate affected by fracture half-length even if the increase ratio of fracture half-length is bigger than matrix permeability, which means that longer fracture half-length only can improve local seepage ability of formation adjacent to hydraulic fractures, but improvement of flow capacity in the whole formation can increase the gas production substantially.Figure 6
Effect of matrix permeabilityK on productivity.Figure 7
Effect of fracture half-lengthLf on productivity.Effect of fracture conductivity on productivity of fractured well has been analyzed under different fracture penetration ratio(Lf/Re). We can see from Figure 8 that gas production increases with the fracture conductivity and fracture penetration ratio. But gas production increase becomes slow when the fracture conductivity increases to a certain value, and this value becomes bigger when fracture penetration ratio increases. For examples, this value will be 2 μm2·cm when fracture penetration ratio is 0.1, but this value will increase to 4 μm2·cm when fracture penetration ratio is 0.3. This chart can help to optimize fracture conductivity under different fracture penetration ratio.Figure 8
Effect of fracture conductivitykf·Wf on productivity.


### 5.3. Effects of Slippage and Knudsen Diffusion on Productivity

The relationship between wellbore pressure of fractured well with finite conductivity and gas production under different slippage factorb has been shown in Figure 9. Slippage effect has negligible influence on productivity of fractured well when wellbore pressure is high, but slippage effect begins to affect gas production when wellbore pressure is lower than 15 MPa and productivity of the fractured well increases with the slippage factor. So decreasing the wellbore pressure of fractured well properly can improve gas production when exploiting shale gas reservoirs. We can see from Figure 10 that gas production rate increases with Knudsen diffusion coefficient, which has a greater impact on gas production than slippage effect.Figure 9
Effect of slippage factorb on productivity.Figure 10
Effect of Knudsen diffusion coefficientDK on productivity.As we can see in Figure11, the productivity of fractured well derived in this paper (Qf) has been compared with production calculated by Darcy formula (QD) under different matrix permeability and wellbore pressure. Figure 11(a) shows that bigger matrix permeability and a lower wellbore pressure lead to a higher productivity of fractured well, and gas production Qf is much bigger than Darcy production rate QD. The effects of Knudsen diffusion and slip flow on productivity of fractured well can be observed by taking the ratio of Qf and QD in Figure 11(b). The effects on gas production increase when matrix permeability becomes smaller and wellbore pressure declines, which become most significant under the smallest matrix permeability (1 × 10−5 mD) and the lowest wellbore pressure (1 MPa) in Figure 11(b).Figure 11
Effects of Knudsen diffusion coefficientDK and slip factor b on productivity.
(a)
Comparison betweenQD and Qf (at different permeability and wellbore pressure)
(b)
The value ofQf/QD at different matrix permeability and wellbore pressureFigure12 represents the relationship of gas production rate Qf and matrix permeability K at different wellbore pressure. For conventional reservoirs, production rate decreases with formation permeability; for shale gas reservoirs, gas production rate decreases with matrix permeability first but then increases with permeability owing to the effects of Knudsen diffusion and slip flow after a certain permeability value, which is called turning point in this paper. The value of turning point decreases when wellbore pressure rises, which is because Knudsen diffusion and slippage effect are more likely to occur in formation with smaller permeability and lower pressure, so a smaller turning point would be required when pressure rises.Figure 12
The relationship between gas production and matrix permeability at different wellbore pressure.Table4 shows the ratio of gas production rate at minimum permeability QKmin⁡ in Figure 12 and production at turning point under different wellbore pressures QTP. The ratios decrease when wellbore pressure increases, which means that Knudsen diffusion and slippage effects have a greater influence on productivity of fractured well under a lower wellbore pressure. But the gas production rate Qf (K = 1 × 10−5 mD, Pw = 1 MPa) is only 1.6 times the production Qf at turning point (Table 4) even though the production Qf is almost 150 times the Darcy production rate QD in Figure 11(b). Knudsen diffusion and slippage effects can improve the productivity at extremely low permeability significantly compared with gas production calculated by Darcy formula (without considering effects of diffusion and slip flow), but the improvement of productivity is still small compared with higher permeability when Knudsen diffusion and slippage effects are considered in both situations.Table 4
Comparison of production at turning point and production at minimum permeability under different wellbore pressure (the gas production rate unit: ×104 m3/d).
P
w/MPa
1
3
5
7
9
11
13
15
17
19
21
23
25
27
Q
K
min
⁡
0.505
0.333
0.250
0.197
0.158
0.127
0.103
0.082
0.065
0.050
0.037
0.025
0.014
0.005
Q
TP
0.315
0.238
0.193
0.160
0.134
0.112
0.093
0.076
0.061
0.048
0.035
0.024
0.014
0.004
Q
K
min
⁡/QTP
1.601
1.398
1.296
1.228
1.178
1.138
1.109
1.082
1.064
1.046
1.034
1.024
1.014
1.010


## 5.1. Effect of Permeability Correction Factor on Productivity

Figure5 presents the effect of permeability correction factor on productivity of fractured well, and we can see from this figure that permeability correction factor has a great influence on productivity. Production of fractured well is very small when our model only considers Darcy flow in nano/microscale pores (ξ=1), but gas production increases a lot when matrix apparent permeability is corrected by permeability correction factor ξ. Field data shows that gas production of a fractured vertical well is 1.2 × 104 m3 when production pressure drop keeps 7 MPa and fracture half-length is 200 m [4]. The calculation results of multiscale non-Darcy seepage model are much more close to the field data compared with the production rate calculated by Darcy formula and production rate calculated by productivity formula only considering first order term of permeability correction factor.Figure 5
Effect of permeability correction factorξ on productivity.


## 5.2. Effects of Formation and Fracture Parameters on Productivity

The effect of matrix permeabilityK on productivity of fractured well is illustrated in Figure 6. Gas production of fractured well increases with the matrix permeability, and gas production increases to 1.71 times, 2.32 times, and 2.87 times when matrix permeability increases to 2 times, 3 times, and 4 times, respectively. Figure 7 reflects productivity of fractured well affected by fracture half-length. Gas production increases to 1.54 times, 2.06 times, and 2.65 times when fracture half-length increases to 3 times, 4 times, and 5 times, respectively. We can find that the increase of gas production rate affected by matrix permeability is still more than production rate affected by fracture half-length even if the increase ratio of fracture half-length is bigger than matrix permeability, which means that longer fracture half-length only can improve local seepage ability of formation adjacent to hydraulic fractures, but improvement of flow capacity in the whole formation can increase the gas production substantially.Figure 6
Effect of matrix permeabilityK on productivity.Figure 7
Effect of fracture half-lengthLf on productivity.Effect of fracture conductivity on productivity of fractured well has been analyzed under different fracture penetration ratio(Lf/Re). We can see from Figure 8 that gas production increases with the fracture conductivity and fracture penetration ratio. But gas production increase becomes slow when the fracture conductivity increases to a certain value, and this value becomes bigger when fracture penetration ratio increases. For examples, this value will be 2 μm2·cm when fracture penetration ratio is 0.1, but this value will increase to 4 μm2·cm when fracture penetration ratio is 0.3. This chart can help to optimize fracture conductivity under different fracture penetration ratio.Figure 8
Effect of fracture conductivitykf·Wf on productivity.


## 5.3. Effects of Slippage and Knudsen Diffusion on Productivity

The relationship between wellbore pressure of fractured well with finite conductivity and gas production under different slippage factorb has been shown in Figure 9. Slippage effect has negligible influence on productivity of fractured well when wellbore pressure is high, but slippage effect begins to affect gas production when wellbore pressure is lower than 15 MPa and productivity of the fractured well increases with the slippage factor. So decreasing the wellbore pressure of fractured well properly can improve gas production when exploiting shale gas reservoirs. We can see from Figure 10 that gas production rate increases with Knudsen diffusion coefficient, which has a greater impact on gas production than slippage effect.Figure 9
Effect of slippage factorb on productivity.Figure 10
Effect of Knudsen diffusion coefficientDK on productivity.As we can see in Figure11, the productivity of fractured well derived in this paper (Qf) has been compared with production calculated by Darcy formula (QD) under different matrix permeability and wellbore pressure. Figure 11(a) shows that bigger matrix permeability and a lower wellbore pressure lead to a higher productivity of fractured well, and gas production Qf is much bigger than Darcy production rate QD. The effects of Knudsen diffusion and slip flow on productivity of fractured well can be observed by taking the ratio of Qf and QD in Figure 11(b). The effects on gas production increase when matrix permeability becomes smaller and wellbore pressure declines, which become most significant under the smallest matrix permeability (1 × 10−5 mD) and the lowest wellbore pressure (1 MPa) in Figure 11(b).Figure 11
Effects of Knudsen diffusion coefficientDK and slip factor b on productivity.
(a)
Comparison betweenQD and Qf (at different permeability and wellbore pressure)
(b)
The value ofQf/QD at different matrix permeability and wellbore pressureFigure12 represents the relationship of gas production rate Qf and matrix permeability K at different wellbore pressure. For conventional reservoirs, production rate decreases with formation permeability; for shale gas reservoirs, gas production rate decreases with matrix permeability first but then increases with permeability owing to the effects of Knudsen diffusion and slip flow after a certain permeability value, which is called turning point in this paper. The value of turning point decreases when wellbore pressure rises, which is because Knudsen diffusion and slippage effect are more likely to occur in formation with smaller permeability and lower pressure, so a smaller turning point would be required when pressure rises.Figure 12
The relationship between gas production and matrix permeability at different wellbore pressure.Table4 shows the ratio of gas production rate at minimum permeability QKmin⁡ in Figure 12 and production at turning point under different wellbore pressures QTP. The ratios decrease when wellbore pressure increases, which means that Knudsen diffusion and slippage effects have a greater influence on productivity of fractured well under a lower wellbore pressure. But the gas production rate Qf (K = 1 × 10−5 mD, Pw = 1 MPa) is only 1.6 times the production Qf at turning point (Table 4) even though the production Qf is almost 150 times the Darcy production rate QD in Figure 11(b). Knudsen diffusion and slippage effects can improve the productivity at extremely low permeability significantly compared with gas production calculated by Darcy formula (without considering effects of diffusion and slip flow), but the improvement of productivity is still small compared with higher permeability when Knudsen diffusion and slippage effects are considered in both situations.Table 4
Comparison of production at turning point and production at minimum permeability under different wellbore pressure (the gas production rate unit: ×104 m3/d).
P
w/MPa
1
3
5
7
9
11
13
15
17
19
21
23
25
27
Q
K
min
⁡
0.505
0.333
0.250
0.197
0.158
0.127
0.103
0.082
0.065
0.050
0.037
0.025
0.014
0.005
Q
TP
0.315
0.238
0.193
0.160
0.134
0.112
0.093
0.076
0.061
0.048
0.035
0.024
0.014
0.004
Q
K
min
⁡/QTP
1.601
1.398
1.296
1.228
1.178
1.138
1.109
1.082
1.064
1.046
1.034
1.024
1.014
1.010


## 6. Conclusions

In this paper, a multiscale comprehensive mathematical model had been established, which can simulate different flow regimes including continuum flow, slip flow, transition flow, and free-molecule flow. We deduced the productivity formula for fractured well with limited conductivity for shale gas reservoirs, and the influencing parameters were analyzed thoroughly. The following conclusions can be drawn.(1)
Flow regimes in shale gas reservoir were analyzed by calculating Knudsen number under different pressure and different pore throat diameter. Gas flow in shale gas reservoirs with nano/microscale pores is a combination of continuum flow, slip flow, transition flow, and free-molecule flow, which cannot simply be represented by Darcy formula anymore.(2)
A new non-Darcy equation was developed based on Beskok-Karniadakis equation which is applicable for all flow regimes, and effects of high order terms of BK equation on permeability correction factor were analyzed under different Knudsen number. Permeability correction factor is increasingly deviating from a Darcy type as Knudsen number increases.(3)
Productivity formula of fractured well satisfying variable mass flowing in fractures was obtained, and the simulated results showed that production rate calculated by non-Darcy seepage model with higher order terms of permeability correction factor is more close to the field production data, compared with the production rate calculated by Darcy formula and production rate calculated by productivity formula only considering first order term of permeability correction factor.(4)
Effects of slippage and diffusion which are changing with pressure and pore sizes in nanoscale pores were considered in simulation. Simulation results showed that slippage effect affects gas production of fractured well only when wellbore pressure less than 15 MPa, and the effects of slippage and diffusion have a greater influence on gas production of reservoirs with smaller matrix permeability, especially when permeability is at nano-Darcy scale. Matrix permeability, fracture half length, fracture penetration ratio, and fracture conductivity can improve gas production at different degrees as well.



---

*Source: 349507-2015-09-27.xml*