# Effect of Matrix-Wellbore Flow and Porosity on Pressure Transient Response in Shale Formation Modeling by Dual Porosity and Dual Permeability System

Authors: Daolun Li; Wenshu Zha; Dewen Zheng; Longjun Zhang; Detang Lu
Journal: Journal of Chemistry
 (2015)

---

## Abstract

A mathematical dual porosity and dual permeability numerical model based on perpendicular bisection (PEBI) grid is developed to describe gas flow behaviors in shale-gas reservoirs by incorporating slippage corrected permeability and adsorbed gas effect. Parametric studies are conducted for a horizontal well with multiple infinite conductivity hydraulic fractures in shale-gas reservoir to investigate effect of matrix-wellbore flow, natural fracture porosity, and matrix porosity. We find that the ratio of fracture permeability to matrix permeability approximately decides the bottom hole pressure (BHP) error caused by omitting the flow between matrix and wellbore and that the effect of matrix porosity on BHP is related to adsorption gas content. When adsorbed gas accounts for large proportion of the total gas storage in shale formation, matrix porosity only has a very small effect on BHP. Otherwise, it has obvious influence. This paper can help us understand the complex pressure transient response due to existence of the adsorbed gas and help petroleum engineers to interpret the field data better.


---

## Body

## 1. Introduction

Shale formation has both free gas and adsorbed gas [1], which usually is described by Langmuir isotherm equation and affects pressure behaviors and production performance significantly. Another characteristic of shale formation is the ultralow permeability. Because of that, flow in shale-gas reservoir leads to several transport mechanisms, such as slip, transition, and molecular diffusion, which depends on the value of Knudsen number [2–7].Due to the incomplete development of the natural fracture network in shale, natural fracture network is not well interconnected and cannot form effective flow channels. In most cases, economic production is possible only if a very complex hydraulic fracture network is created that effectively connects a huge reservoir surface area to the wellbore [8]. Therefore, the shale formation consists of matrix, natural fracture, and hydraulic fracture. Many researchers use analytical or semianalytical solution of dual porosity model to study the pressure transient behaviors or production data analysis for fractured horizontal wells [9–13]. Clarkson et al. [14] and Yan et al. [15] use numerical simulation for flow behaviors in single porosity or three porosities model to study flow in shale gas. In order to decrease capillary pressure to improve gas flow, wettability alteration by chemical treatments is proposed [16].In the above studies, the dual porosity and dual permeability (DPDP) model, in which there are flow in matrix blocks and fracture blocks and flow between them, is seldom mentioned.In this paper, we first establish a mathematical dual porosity and dual permeability shale-gas flow model incorporating the gas retention and slip flow mechanism. Then the nonlinear pressure equation is solved by numerical way based on PBEI grid. Lastly, we use the numerical model to study the effect of flow between matrix and fracture, effect of natural fracture porosity, and effect of matrix porosity on pressure transient response for wells in shale-gas reservoir.


## 2. Flow Model Scheme of Shale Gas

### 2.1. Flow Model

Taking the adsorbed gas as accumulation term and using apparent permeability to describe the slip flow, we can obtain the following equation for matrix system and fracture system, respectively:(1)∇·kmμgBg∇pm-γg∇Z-αkmpm-pfμgBg=∂∂tϕmBg+∂∂tρsVads+qgsc,m,(2)∇·kfμgBg∇pf-γg∇Z+αkmpm-pfμgBg=∂∂tϕfBg+qgsc,f,where subscript m denotes matrix system, subscript f denotes natural fracture system, ρs is shale density in Kg/m3, Bg is formation volume factor of gas in fraction, ϕ is porosity in fraction, Vads is adsorption content in m3/Kg, and qgsc,m and qgsc,f are volume flow rate of the well divided by the volume of the well grid in 1/s.Many attempts have been made in the past to define the shape factor. The classical formula is given by Kazemi [17] as follows:(3)α=41Lx2+1Ly2+1Lz2,where Lx, Ly, and Lz are the natural fracture spacing in x, y, and z directions in meter, respectively, and α is shape factor to describe the interporosity flow between matrix and natural fracture. In this paper, for simplicity, we suppose that Lx, Ly, and Lz are equal and just use Lx to express the natural fracture spacing.LetQ in m3/s be the production of the well at surface condition, and considering wellbore storage, the flow equation in the wellbore is given by(4)2πkhμgBgln⁡re/rw+Spi-pwf-CΔtpwfn+1-pwfn=Q,where pwf is BHP in Pa, re is equivalent well radius in m and rw is actual well radius in m, C is wellbore storage in m3/Pa, and n is time step.


### 2.2. Apparent Permeability

In order to describe the flow mechanism caused by the ultralow permeability, the intrinsic permeability in (1) and (2) is replaced by apparent permeability, which can be described with Knudsen number for the slip flow [18] as follows:(5)k=k01+4Kn1+Kn.According to the definition of mean free path of gas molecules, Knudsen number is expressed as(6)Kn=μdpπRT2M,where μ is the gas viscosity in Pa·s, R=8314 J/kmol/K is the universal gas constant, T is the absolute temperature in K, and M is the molecular mass of gas in kg/kmol. Civan [4] estimates the equivalent hydraulic radius:(7)d=22τkϕ,where τ is the tortuosity and ϕ is the porosity of porous media.According to Schaaf and Chambre [19], free molecular flow happens if Kn>10, when collision between the gas molecules and the pore wall dominates. For 0.1≤Kn≤10, it is considered as transition regime. For 10-3<Kn<0.1, the flow at the wall cannot be neglected, and slip flow exists. For Kn<10-3, collision between molecules and the pore wall can be neglected, and the flow behaviors can be described by Darcy’s law. In this paper, the flow regime is limited to slip flow.


### 2.3. Isotherm Sorption of Gas

The amount of adsorption in (1) can be modeled by the Langmuir isotherm equation:(8)Vads=VLppL+p,where Vads is the standard volume of gas adsorbed per unit shale mass in m3/Kg under the pressure p, VL denotes the ultimate adsorption amount in m3/Kg, p is the gas pressure in Pa, and pL is the pressure when adsorption reaches half of the ultimate adsorption amount in Pa.


## 2.1. Flow Model

Taking the adsorbed gas as accumulation term and using apparent permeability to describe the slip flow, we can obtain the following equation for matrix system and fracture system, respectively:(1)∇·kmμgBg∇pm-γg∇Z-αkmpm-pfμgBg=∂∂tϕmBg+∂∂tρsVads+qgsc,m,(2)∇·kfμgBg∇pf-γg∇Z+αkmpm-pfμgBg=∂∂tϕfBg+qgsc,f,where subscript m denotes matrix system, subscript f denotes natural fracture system, ρs is shale density in Kg/m3, Bg is formation volume factor of gas in fraction, ϕ is porosity in fraction, Vads is adsorption content in m3/Kg, and qgsc,m and qgsc,f are volume flow rate of the well divided by the volume of the well grid in 1/s.Many attempts have been made in the past to define the shape factor. The classical formula is given by Kazemi [17] as follows:(3)α=41Lx2+1Ly2+1Lz2,where Lx, Ly, and Lz are the natural fracture spacing in x, y, and z directions in meter, respectively, and α is shape factor to describe the interporosity flow between matrix and natural fracture. In this paper, for simplicity, we suppose that Lx, Ly, and Lz are equal and just use Lx to express the natural fracture spacing.LetQ in m3/s be the production of the well at surface condition, and considering wellbore storage, the flow equation in the wellbore is given by(4)2πkhμgBgln⁡re/rw+Spi-pwf-CΔtpwfn+1-pwfn=Q,where pwf is BHP in Pa, re is equivalent well radius in m and rw is actual well radius in m, C is wellbore storage in m3/Pa, and n is time step.


## 2.2. Apparent Permeability

In order to describe the flow mechanism caused by the ultralow permeability, the intrinsic permeability in (1) and (2) is replaced by apparent permeability, which can be described with Knudsen number for the slip flow [18] as follows:(5)k=k01+4Kn1+Kn.According to the definition of mean free path of gas molecules, Knudsen number is expressed as(6)Kn=μdpπRT2M,where μ is the gas viscosity in Pa·s, R=8314 J/kmol/K is the universal gas constant, T is the absolute temperature in K, and M is the molecular mass of gas in kg/kmol. Civan [4] estimates the equivalent hydraulic radius:(7)d=22τkϕ,where τ is the tortuosity and ϕ is the porosity of porous media.According to Schaaf and Chambre [19], free molecular flow happens if Kn>10, when collision between the gas molecules and the pore wall dominates. For 0.1≤Kn≤10, it is considered as transition regime. For 10-3<Kn<0.1, the flow at the wall cannot be neglected, and slip flow exists. For Kn<10-3, collision between molecules and the pore wall can be neglected, and the flow behaviors can be described by Darcy’s law. In this paper, the flow regime is limited to slip flow.


## 2.3. Isotherm Sorption of Gas

The amount of adsorption in (1) can be modeled by the Langmuir isotherm equation:(8)Vads=VLppL+p,where Vads is the standard volume of gas adsorbed per unit shale mass in m3/Kg under the pressure p, VL denotes the ultimate adsorption amount in m3/Kg, p is the gas pressure in Pa, and pL is the pressure when adsorption reaches half of the ultimate adsorption amount in Pa.


## 3. Numerical Calculation Scheme

PEBI grids can overcome structure grid disadvantages, such as inflexibility in description of geologic features, inaccurate representation of multiwell locations, and grid orientation, and are especially widely used in numerical well test [20–23]. PEBI grid is actually Voronoi block, which is defined as the region of space that is closer to its grid point than to any other grid points, and its boundaries are normal to lines joining the nodes on the two sides of each boundary.Figure1(a) gives the grids of horizontal well with 5 hydraulic fractures. Figure 1(b) gives the fracture grids and horizontal well grids in detail. In order to capture the characteristics of the transient flow, the grid points are distributed by radial growth near the well and fractures. The well and fractures are modeled as grids and treated as infinite conductivity [24].Figure 1
Hydraulically fractured horizontal well grid without cell points in 2D form.
(a)
Grid
(b)
Part of the horizontal well and fracturesEquations (1), (2), and (4) form an equation system. There are a total of 2 variables of pf, pm for a grid and one variable for one well producing under constant flow rate.Nonlinear pressure equations (7) and (8) are implicitly solved by using the matrix solver GMRES (generalized minimal residual method) [25].


## 4. Simulation and Analysis

### 4.1. Validation

We use the IMEX model of STARS software to validate our core codes. STARS is a commercial software developed by CMG and is widely accepted in petroleum industry. The code we used here to compare with STARS (2012 version) is only the conventional mass balance equation containing gas phase flowing in DPDP and without adsorption and apparent permeability correction. The conventional balance equation is the core part of our model. Our code will be convincing if the core part of our model is validated.The horizontal wellbore only connects the fractures in this paper; thus, horizontal well with one fracture is equivalent to fractured vertical well. In our numerical code, the well and hydraulic fracture are treated as inner boundary.The grid number is81×81 in IMEX model. The size of grid is 10 m. The blocks including the hydraulic fracture are subdivided into 11×3 grids. The main parameters used in the validation are shown in Table 1.Table 1
Data used in validation.
Name
Value
Unit
Reservoir size
810 × 810
m
Boundary condition
Closed
Initial pressure
20
MPa
Reservoir thickness
30
m
Rock compressibility
0.0000145
1/MPa
Temperature
122
°F
Gas gravity
0.7
Half of hydraulic fracture
75
m
Natural fracture space
10
m
Fracture
Permeability
700
mD
Porosity
0.05
Matrix
Permeability
0.1
mD
Porosity
0.25
Horizontal well
Skin factor
0
Well storage
0
m3/MPa
Constant rate (STC)
12000
m3/day
Production time
200
Day
Deviation factorZ
Dranchuk-Abou-Kassem equation
Gas viscosity
Lee equationThe well is in the middle of reservoir and produces gas at a constant flow rate of 12000m3/day in a closed boundary reservoir. The BHP is compared between our code and STARS as shown in Figure 2. The almost overlap of the two output curves validates our code.Figure 2
Comparison of bottom hole flowing pressure.


### 4.2. Reservoir Parameters and Grid

Table2 gives the main parameters used for the simulation.Table 2
Simulation parameters (if not specified additionally).
Name
Value
Unit
Reservoir size
80000 × 80000
m
Initial pressure
20
MPa
Reservoir thickness
30
m
Rock compressibility
0.00015
1/MPa
Langmuir adsorptionVL
0.02
Std m3/Kg
Langmuir pressurepL
7.5
MPa
Horizontal well
Skin factor
0
Well storage
10
m3/MPa
Constant rate
10000
m3/day
Production time
10000
Day
Shut time
10000
Day
Deviation factorZ
Dranchuk-Abou-Kassem equation
Gas viscosity
Lee equationTable3 gives parameters about the DPDP and hydraulic fractures. The total simulation time is 20000 days (about 54.8 years). Gas viscosity and volume formation factor are calculated with PVT [26, 27] and are changing with pressure and obtained through iterative process.Table 3
DPDP data (if not specified additionally).
Matrix data
Natural fracture data
Hydraulic fracture data
k
m = 0.001 mD, ϕm = 0.1
k
f = 0.1 mD, ϕf = 0.001
Three hydraulic fractures with 40 meters long and fracture spacing is 280 meters


### 4.3. Effect of Flow between Matrix and Wellbore with Adsorption

In Figure3, ultimate adsorption capacity (VL) is 0.02 m3/Kg, and constant flow rate (Q) is 10000 m3/day. The other parameters are shown in the figures.Figure 3
Effect of flow during matrix and wellbore on BHP underLx=1 m and =10000 m3/day.
(a)
BHP comparison forkf=10km
(b)
BHP comparison forkf=100km
(c)
Derivative comparison forkf=10km
(d)
Derivative comparison forkf=100kmFigure3(a) shows that when fracture permeability is 10 times of the matrix permeability, the BHP difference at 10000th day is about 0.11 MPa. The pressure change is about 1.157 MPa under dual well model. If well-matrix flow is omitted, the BHP error is 9.03% (Table 4).Table 4
Error percentage of BHP without flow between wellbore and matrix.
Figure
Q
Parameter
Error percentage
Figure3(a)
10000 m3/day
L
x
=
1 m, kf=10km
9.03%
Figure4(a)
10000 m3/day
L
x
=
10 m, kf=10km
9.26%
Figure5(a)
50000 m3/day
L
x
=
10 m, kf=10km
12.6%
Figure3(b)
10000 m3/day
L
x
=
1 m, kf=100km
1.004%
Figure4(b)
10000 m3/day
L
x
=
1 m, kf=100km
1.028%
Figure5(b)
50000 m3/day
L
x
=
10 m, kf=100km
1.19%When the fracture permeability is 100 times of the matrix permeability (Figure3(b)), the BHP difference at 10000th day is 0.0146 MPa. If matrix-wellbore flow is omitted, the BHP error is 1.004% (Table 4).By comparison of the BHP error with the ratio ofkm/kf, we know that the ratio of the matrix permeability to the fracture permeability approximately decides the BHP error caused by omitting the matrix-wellbore flow.The corresponding pressure change and pressure derivative comparison is shown in Figures3(c) and 3(d).In Figure4, except the natural fracture spacing of Lx=10 m, other data are the same as data in Figure 3, from which we can see the influence of interporosity flow ability on the BHP error. The BHP error for kf=10km in Figure 4(a) is 9.26% and BHP error for kf=100km in Figure 4(b) is 1.028%, which is shown in Table 4. The BHP errors in Figure 4 are correspondingly larger than that in Figure 3, which is caused by smaller interporosity ability.Figure 4
Effect of flow during matrix and wellbore on BHP underLx=10 m and =10000 m3/day.
(a)
BHP comparison forkf=10km
(b)
BHP comparison forkf=100kmWhen interporosity ability is small, due to the large permeability of fracture system, some gas needs to flow from matrix system to fracture system. If the interporosity ability is not large enough, larger pressure difference between matrix and fracture is required, which causes larger error if the matrix-wellbore flow is omitted. Therefore, the low interporosity flow ability can enlarge the error when we omit the matrix-wellbore flow.In Figure5, the constant flow rate increases from 10000 m3/day to 50000 m3/day compared to the parameters in Figure 4. In Figure 5(a), the fracture permeability is 10 times of the matrix permeability, and the BHP difference at 10000th day is 0.77247 MPa. The BHP error is 12.6% if matrix-wellbore flow is neglected. In Figure 5(b), the fracture permeability is 100 times of the matrix permeability, the BHP difference at 10000th day is about 0.09457 MPa, and the BHP error is 1.1924% if the flow between wellbore and matrix is omitted.Figure 5
Effect of flow during matrix and wellbore underLx=10 m and =50000 m3/day.
(a)
BHP comparison forkf=10km
(b)
BHP comparison forkf=100kmCompared to the value of BHP error in Figure4 (see Table 4), the BHP error under Q=50000 m3/day increases, which shows that large flow rate leads to the large BHP error. We also find that the increase of BHP error under kf=10km is larger than the increase of BHP error under kf=100km, which shows that large ratio of kf/km decreases BHP error if matrix-wellbore flow is omitted.Table4 gives the BHP error comparison for Figures 3, 4, and 5. When natural fracture spacing increases from 1 m to 10 m under =10000 m3/day, the BHP error increases from 9.03% to 9.26% under kf=10km, which means that BHP error increases with decrease of interporosity flow ability.When constant flow rate increases 5 times, the BHP error increases from 9.26% to 12.6% underkf=10km, while the BHP error only increases from 1.028% to 1.19% under kf=100km, which has a smaller increase rate than that under kf=10km. Therefore, when ratio of kf/km is large, neglecting of the matrix-wellbore flow only causes very small BHP error. However, if ratio of kf/km is relatively large, the matrix-wellbore flow should not be omitted.


### 4.4. Effect of Porosity on Transient Pressure Response with Adsorption Gas

#### 4.4.1. Effect of Fracture Porosity on Transient Pressure Response

Figure6 shows the effect of natural fracture porosity on pressure transient responses for a horizontal well with 3 hydraulic fractures with ultimate adsorption capacity of VL=0.1 m3/Kg and natural fracture spacing of Lx=1 m.Figure 6
Effect of porosities of fracture system on pressure response forVL=0.1m3/Kg, φm=0.05, and Lx=1 m.
(a)
BHP for different fracture porosities
(b)
Pressure derivative of the pressure drawdown periodThe effect of fracture porosity on BHP is so small that it is not easy to identify the difference in Figure6(a). From pressure derivative curves in Figure 6(b), we can see the difference of the early linear flow regime. However, the time of the difference only lasts 0.01 days. When fracture porosity is 0.001, the early linear flow regime is obvious, while the early linear flow regime is not clear under fracture porosity of 0.0001. This is because the natural fracture system is the main flow channel. Larger fracture porosity can make the linear flow in nature fracture last longer and linear flow behavior becomes obvious.


#### 4.4.2. Effect of Matrix Porosity on Transient Pressure Response

Figure7 shows the effect of matrix porosity on pressure transient responses of a horizontal well with 3 hydraulic fractures for ultimate adsorption capacity of VL=1 m3/Kg, natural fracture spacing of Lx=1 m.Figure 7
Effect of matrix porosities on pressure response forVL=1 m3/Kg, Lx=1 m, and φf=0.0001.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porositiesFigure7 shows that the matrix porosity has minimal effect on BHP when the gas ultimate adsorption capacity is 1 m3/Kg. Matrix system is the storage space and its porosity should affect behaviors of BHP. However, in shale gas, the situation is different. When pressure drops, the adsorbed gas desorbed. If the ultimate adsorption capacity is big, the magnitude of the free gas obtained from the adsorption gas is much bigger than the free gas from the matrix porosity, which reduces the effect of the matrix porosity.When the ultimate adsorption capacity is small, the desorbed gas will have less effect on the BHP, and the effect of the matrix porosity on BHP will appear. Figure8 gives effect of matrix porosity on pressure response with ultimate adsorption capacity of 0.05 m3/Kg. With small adsorption content, the BHP difference in Figure 8(a) is 0.28 MPa at the end of the production period due to the porosity difference between matrix systems. When adsorption content decreases, the free gas will account for more proportion of the flow rate. Therefore, matrix porosity has more influence on the pressure transient response.Figure 8
Effect of matrix porosities on pressure response withVL=0.05 m3/Kg, Lx=1 m, and φf=0.0001.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porositiesFigure9 gives the case of no adsorption gas. The BHP difference in Figure 9(a) is 0.46 MPa, which is much larger than the BHP difference under VL=0.05 m3/Kg in Figure 8(a). The pressure change and pressure derivative in Figures 7, 8, and 9 together clearly show the effect of matrix porosity under different ultimate adsorption capacities.Figure 9
Effect of matrix porosities on pressure response withVL=0 m3/Kg.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porosities


## 4.1. Validation

We use the IMEX model of STARS software to validate our core codes. STARS is a commercial software developed by CMG and is widely accepted in petroleum industry. The code we used here to compare with STARS (2012 version) is only the conventional mass balance equation containing gas phase flowing in DPDP and without adsorption and apparent permeability correction. The conventional balance equation is the core part of our model. Our code will be convincing if the core part of our model is validated.The horizontal wellbore only connects the fractures in this paper; thus, horizontal well with one fracture is equivalent to fractured vertical well. In our numerical code, the well and hydraulic fracture are treated as inner boundary.The grid number is81×81 in IMEX model. The size of grid is 10 m. The blocks including the hydraulic fracture are subdivided into 11×3 grids. The main parameters used in the validation are shown in Table 1.Table 1
Data used in validation.
Name
Value
Unit
Reservoir size
810 × 810
m
Boundary condition
Closed
Initial pressure
20
MPa
Reservoir thickness
30
m
Rock compressibility
0.0000145
1/MPa
Temperature
122
°F
Gas gravity
0.7
Half of hydraulic fracture
75
m
Natural fracture space
10
m
Fracture
Permeability
700
mD
Porosity
0.05
Matrix
Permeability
0.1
mD
Porosity
0.25
Horizontal well
Skin factor
0
Well storage
0
m3/MPa
Constant rate (STC)
12000
m3/day
Production time
200
Day
Deviation factorZ
Dranchuk-Abou-Kassem equation
Gas viscosity
Lee equationThe well is in the middle of reservoir and produces gas at a constant flow rate of 12000m3/day in a closed boundary reservoir. The BHP is compared between our code and STARS as shown in Figure 2. The almost overlap of the two output curves validates our code.Figure 2
Comparison of bottom hole flowing pressure.


## 4.2. Reservoir Parameters and Grid

Table2 gives the main parameters used for the simulation.Table 2
Simulation parameters (if not specified additionally).
Name
Value
Unit
Reservoir size
80000 × 80000
m
Initial pressure
20
MPa
Reservoir thickness
30
m
Rock compressibility
0.00015
1/MPa
Langmuir adsorptionVL
0.02
Std m3/Kg
Langmuir pressurepL
7.5
MPa
Horizontal well
Skin factor
0
Well storage
10
m3/MPa
Constant rate
10000
m3/day
Production time
10000
Day
Shut time
10000
Day
Deviation factorZ
Dranchuk-Abou-Kassem equation
Gas viscosity
Lee equationTable3 gives parameters about the DPDP and hydraulic fractures. The total simulation time is 20000 days (about 54.8 years). Gas viscosity and volume formation factor are calculated with PVT [26, 27] and are changing with pressure and obtained through iterative process.Table 3
DPDP data (if not specified additionally).
Matrix data
Natural fracture data
Hydraulic fracture data
k
m = 0.001 mD, ϕm = 0.1
k
f = 0.1 mD, ϕf = 0.001
Three hydraulic fractures with 40 meters long and fracture spacing is 280 meters


## 4.3. Effect of Flow between Matrix and Wellbore with Adsorption

In Figure3, ultimate adsorption capacity (VL) is 0.02 m3/Kg, and constant flow rate (Q) is 10000 m3/day. The other parameters are shown in the figures.Figure 3
Effect of flow during matrix and wellbore on BHP underLx=1 m and =10000 m3/day.
(a)
BHP comparison forkf=10km
(b)
BHP comparison forkf=100km
(c)
Derivative comparison forkf=10km
(d)
Derivative comparison forkf=100kmFigure3(a) shows that when fracture permeability is 10 times of the matrix permeability, the BHP difference at 10000th day is about 0.11 MPa. The pressure change is about 1.157 MPa under dual well model. If well-matrix flow is omitted, the BHP error is 9.03% (Table 4).Table 4
Error percentage of BHP without flow between wellbore and matrix.
Figure
Q
Parameter
Error percentage
Figure3(a)
10000 m3/day
L
x
=
1 m, kf=10km
9.03%
Figure4(a)
10000 m3/day
L
x
=
10 m, kf=10km
9.26%
Figure5(a)
50000 m3/day
L
x
=
10 m, kf=10km
12.6%
Figure3(b)
10000 m3/day
L
x
=
1 m, kf=100km
1.004%
Figure4(b)
10000 m3/day
L
x
=
1 m, kf=100km
1.028%
Figure5(b)
50000 m3/day
L
x
=
10 m, kf=100km
1.19%When the fracture permeability is 100 times of the matrix permeability (Figure3(b)), the BHP difference at 10000th day is 0.0146 MPa. If matrix-wellbore flow is omitted, the BHP error is 1.004% (Table 4).By comparison of the BHP error with the ratio ofkm/kf, we know that the ratio of the matrix permeability to the fracture permeability approximately decides the BHP error caused by omitting the matrix-wellbore flow.The corresponding pressure change and pressure derivative comparison is shown in Figures3(c) and 3(d).In Figure4, except the natural fracture spacing of Lx=10 m, other data are the same as data in Figure 3, from which we can see the influence of interporosity flow ability on the BHP error. The BHP error for kf=10km in Figure 4(a) is 9.26% and BHP error for kf=100km in Figure 4(b) is 1.028%, which is shown in Table 4. The BHP errors in Figure 4 are correspondingly larger than that in Figure 3, which is caused by smaller interporosity ability.Figure 4
Effect of flow during matrix and wellbore on BHP underLx=10 m and =10000 m3/day.
(a)
BHP comparison forkf=10km
(b)
BHP comparison forkf=100kmWhen interporosity ability is small, due to the large permeability of fracture system, some gas needs to flow from matrix system to fracture system. If the interporosity ability is not large enough, larger pressure difference between matrix and fracture is required, which causes larger error if the matrix-wellbore flow is omitted. Therefore, the low interporosity flow ability can enlarge the error when we omit the matrix-wellbore flow.In Figure5, the constant flow rate increases from 10000 m3/day to 50000 m3/day compared to the parameters in Figure 4. In Figure 5(a), the fracture permeability is 10 times of the matrix permeability, and the BHP difference at 10000th day is 0.77247 MPa. The BHP error is 12.6% if matrix-wellbore flow is neglected. In Figure 5(b), the fracture permeability is 100 times of the matrix permeability, the BHP difference at 10000th day is about 0.09457 MPa, and the BHP error is 1.1924% if the flow between wellbore and matrix is omitted.Figure 5
Effect of flow during matrix and wellbore underLx=10 m and =50000 m3/day.
(a)
BHP comparison forkf=10km
(b)
BHP comparison forkf=100kmCompared to the value of BHP error in Figure4 (see Table 4), the BHP error under Q=50000 m3/day increases, which shows that large flow rate leads to the large BHP error. We also find that the increase of BHP error under kf=10km is larger than the increase of BHP error under kf=100km, which shows that large ratio of kf/km decreases BHP error if matrix-wellbore flow is omitted.Table4 gives the BHP error comparison for Figures 3, 4, and 5. When natural fracture spacing increases from 1 m to 10 m under =10000 m3/day, the BHP error increases from 9.03% to 9.26% under kf=10km, which means that BHP error increases with decrease of interporosity flow ability.When constant flow rate increases 5 times, the BHP error increases from 9.26% to 12.6% underkf=10km, while the BHP error only increases from 1.028% to 1.19% under kf=100km, which has a smaller increase rate than that under kf=10km. Therefore, when ratio of kf/km is large, neglecting of the matrix-wellbore flow only causes very small BHP error. However, if ratio of kf/km is relatively large, the matrix-wellbore flow should not be omitted.


## 4.4. Effect of Porosity on Transient Pressure Response with Adsorption Gas

### 4.4.1. Effect of Fracture Porosity on Transient Pressure Response

Figure6 shows the effect of natural fracture porosity on pressure transient responses for a horizontal well with 3 hydraulic fractures with ultimate adsorption capacity of VL=0.1 m3/Kg and natural fracture spacing of Lx=1 m.Figure 6
Effect of porosities of fracture system on pressure response forVL=0.1m3/Kg, φm=0.05, and Lx=1 m.
(a)
BHP for different fracture porosities
(b)
Pressure derivative of the pressure drawdown periodThe effect of fracture porosity on BHP is so small that it is not easy to identify the difference in Figure6(a). From pressure derivative curves in Figure 6(b), we can see the difference of the early linear flow regime. However, the time of the difference only lasts 0.01 days. When fracture porosity is 0.001, the early linear flow regime is obvious, while the early linear flow regime is not clear under fracture porosity of 0.0001. This is because the natural fracture system is the main flow channel. Larger fracture porosity can make the linear flow in nature fracture last longer and linear flow behavior becomes obvious.


### 4.4.2. Effect of Matrix Porosity on Transient Pressure Response

Figure7 shows the effect of matrix porosity on pressure transient responses of a horizontal well with 3 hydraulic fractures for ultimate adsorption capacity of VL=1 m3/Kg, natural fracture spacing of Lx=1 m.Figure 7
Effect of matrix porosities on pressure response forVL=1 m3/Kg, Lx=1 m, and φf=0.0001.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porositiesFigure7 shows that the matrix porosity has minimal effect on BHP when the gas ultimate adsorption capacity is 1 m3/Kg. Matrix system is the storage space and its porosity should affect behaviors of BHP. However, in shale gas, the situation is different. When pressure drops, the adsorbed gas desorbed. If the ultimate adsorption capacity is big, the magnitude of the free gas obtained from the adsorption gas is much bigger than the free gas from the matrix porosity, which reduces the effect of the matrix porosity.When the ultimate adsorption capacity is small, the desorbed gas will have less effect on the BHP, and the effect of the matrix porosity on BHP will appear. Figure8 gives effect of matrix porosity on pressure response with ultimate adsorption capacity of 0.05 m3/Kg. With small adsorption content, the BHP difference in Figure 8(a) is 0.28 MPa at the end of the production period due to the porosity difference between matrix systems. When adsorption content decreases, the free gas will account for more proportion of the flow rate. Therefore, matrix porosity has more influence on the pressure transient response.Figure 8
Effect of matrix porosities on pressure response withVL=0.05 m3/Kg, Lx=1 m, and φf=0.0001.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porositiesFigure9 gives the case of no adsorption gas. The BHP difference in Figure 9(a) is 0.46 MPa, which is much larger than the BHP difference under VL=0.05 m3/Kg in Figure 8(a). The pressure change and pressure derivative in Figures 7, 8, and 9 together clearly show the effect of matrix porosity under different ultimate adsorption capacities.Figure 9
Effect of matrix porosities on pressure response withVL=0 m3/Kg.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porosities


## 4.4.1. Effect of Fracture Porosity on Transient Pressure Response

Figure6 shows the effect of natural fracture porosity on pressure transient responses for a horizontal well with 3 hydraulic fractures with ultimate adsorption capacity of VL=0.1 m3/Kg and natural fracture spacing of Lx=1 m.Figure 6
Effect of porosities of fracture system on pressure response forVL=0.1m3/Kg, φm=0.05, and Lx=1 m.
(a)
BHP for different fracture porosities
(b)
Pressure derivative of the pressure drawdown periodThe effect of fracture porosity on BHP is so small that it is not easy to identify the difference in Figure6(a). From pressure derivative curves in Figure 6(b), we can see the difference of the early linear flow regime. However, the time of the difference only lasts 0.01 days. When fracture porosity is 0.001, the early linear flow regime is obvious, while the early linear flow regime is not clear under fracture porosity of 0.0001. This is because the natural fracture system is the main flow channel. Larger fracture porosity can make the linear flow in nature fracture last longer and linear flow behavior becomes obvious.


## 4.4.2. Effect of Matrix Porosity on Transient Pressure Response

Figure7 shows the effect of matrix porosity on pressure transient responses of a horizontal well with 3 hydraulic fractures for ultimate adsorption capacity of VL=1 m3/Kg, natural fracture spacing of Lx=1 m.Figure 7
Effect of matrix porosities on pressure response forVL=1 m3/Kg, Lx=1 m, and φf=0.0001.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porositiesFigure7 shows that the matrix porosity has minimal effect on BHP when the gas ultimate adsorption capacity is 1 m3/Kg. Matrix system is the storage space and its porosity should affect behaviors of BHP. However, in shale gas, the situation is different. When pressure drops, the adsorbed gas desorbed. If the ultimate adsorption capacity is big, the magnitude of the free gas obtained from the adsorption gas is much bigger than the free gas from the matrix porosity, which reduces the effect of the matrix porosity.When the ultimate adsorption capacity is small, the desorbed gas will have less effect on the BHP, and the effect of the matrix porosity on BHP will appear. Figure8 gives effect of matrix porosity on pressure response with ultimate adsorption capacity of 0.05 m3/Kg. With small adsorption content, the BHP difference in Figure 8(a) is 0.28 MPa at the end of the production period due to the porosity difference between matrix systems. When adsorption content decreases, the free gas will account for more proportion of the flow rate. Therefore, matrix porosity has more influence on the pressure transient response.Figure 8
Effect of matrix porosities on pressure response withVL=0.05 m3/Kg, Lx=1 m, and φf=0.0001.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porositiesFigure9 gives the case of no adsorption gas. The BHP difference in Figure 9(a) is 0.46 MPa, which is much larger than the BHP difference under VL=0.05 m3/Kg in Figure 8(a). The pressure change and pressure derivative in Figures 7, 8, and 9 together clearly show the effect of matrix porosity under different ultimate adsorption capacities.Figure 9
Effect of matrix porosities on pressure response withVL=0 m3/Kg.
(a)
BHP for different matrix porosities
(b)
Pressure derivative for different matrix porosities


## 5. Conclusion

In this paper, we first developed a dual porosity and dual permeability model for a multistage fractured horizontal well in shale-gas reservoirs incorporating the slippage corrected gas permeability and gas adsorption. Then we developed fully implicit numerical simulation model using PEBI grid and finally conducted parametric studies to quantify the transient pressure behaviors for flow in shale-gas formation. Our conclusions could be summarized as below.(1)
The ratio of fracture permeability to matrix permeability approximately decides the BHP error caused by omitting the flow between the matrix and the wellbore. The BHP error is also related to the interporosity flow ability. The smaller interporosity flow ability is, the larger BHP error is.(2)
In shale-gas formation with adsorption gas, the effect of matrix porosity on BHP is related to adsorption gas content. The desorbed free gas in matrix can offset the effect of the matrix porosity. When adsorbed gas accounts for a big part of total gas storage in shale formation, the change of the matrix porosity only has a very small effect of the total gas, which leads to a very small effect on BHP.(3)
Although effect of natural fracture system porosity on BHP is not discernable, it influences the flow behaviors during the early stage significantly.This paper can help us understand the complex pressure transient response due to existence of adsorbed gas and help petroleum engineers to interpret the field data better.



---

*Source: 426860-2015-09-28.xml*