Next Article in Journal
Life Cycle Analysis of the Bioethanol Production from Food Waste—A Review
Previous Article in Journal
A Power Hardware-in-the-Loop Based Method for FAPR Compliance Testing of the Wind Turbine Converters Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pressure-Transient Performances of Fractured Horizontal Wells in the Compartmentalized Heterogeneous Unconventional Reservoirs

1
School of Petroleum Engineering, China University of Petroleum, Qingdao 266580, China
2
The No.1 Gas Production Plant, Sinopec Southwest Oil and Gas Company, Deyang 618000, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(19), 5204; https://doi.org/10.3390/en13195204
Submission received: 31 August 2020 / Revised: 28 September 2020 / Accepted: 29 September 2020 / Published: 6 October 2020

Abstract

:
In order to investigate pressure performance of multiple fractured horizontal wells (MFHWs) penetrating heterogeneous unconventional reservoir and avoid the high computational cost of numerical simulation, a semi-analytical model for MFHWs combining Green function solution and boundary element method has been obtained, where the reservoir is divided into different homogeneous substructures and coupled at interface boundaries by plane source function in a closed rectangular parallelepiped. Hydraulic fractures are assumed uniform flux and dual porosity model is used for natural fractures system. Then the model is validated by compared with analytical solution of MFHWs in a homogeneous reservoir and trilinear flow model, which shows that this model can achieve high accuracy even with a small interface discretization number, and it can consider the radial flow around each hydraulic fractures. Finally, the pressure responses with heterogeneous parameters of reservoirs are discussed including heterogeneous permeability, non-uniform block-length and fracture half-length distribution as well as dual porosity parameters like elastic storage ratio and crossflow ratio.

Graphical Abstract

1. Introduction

With the rapid development of unconventional reservoirs, the horizontal well intercepted by hydraulic fractures has become the main technology [1,2], which can significantly increase contact areas between fractures and formation zone. Thus, the transient performances of Multiple Fractured Horizontal Wells (MFHWs) have attracted more and more scholars’ attention. There have been extensive studies [3,4,5] on the MFHWs of unconventional reservoir, which generally can be classified into analytical model, numerical model and semi-analytical model.
To consider the effect of formation heterogeneity, one of the most common analytical models is the trilinear flow model with three different flow areas of outer region, inner region and hydraulic fracture, which is proposed by Brown [6]. To describe the hydraulic fractures surrounding by the stimulated area, Stalgorova [7,8] subdivided the trilinear flow regions into five flow regions model. Guo [9] extended the flow area into seven parts including various fracture distribution along the horizontal wellbore. Later, on account of the fractal nature of fractured media, the new formulations are implemented in the trilinear flow model using fractal theory [3,10,11]. The fundamental petro-physical characteristic of unconventional reservoirs is considered, but all multi-linear models are based on two critical assumptions. One assumption is that fluid flow obeys linear flow in all regions. As a result, some other flow regimes are ignored such as first radial flow around hydraulic fractures and later elliptical flow [12]. The other assumption is that the formation properties around each fracture are assumed to be the same, and the interference between different fractures is not investigated. However, Zawilla [13] and Liu [14] indicated that heterogeneous properties are commonly existed in the unconventional reservoir formation and different fracturing stages.
The numerical simulation is an effective method to model the flow performance of MFHWs, which can deal with the formation heterogeneity easily. Especially, Discrete Fracture Model (DFM) [4,15,16] is used for complex fracture network with high computational costs, which accounts the effects of individual fractures explicitly. In order to avoid lots of meshing of unstructured grid around each fracture, Embedded Discrete Fracture Method (EDFM) is proposed [17] by orthogonal coarse grid instead of fine small unstructured grid, but the transient performance of MFHWs needs detailed pressure solution by small space scale and time scale.
Since Gringarten [18] (1973) introduced the source/Green’s function solution into the fluid flow in porous media, this semi-analytical method has greatly expanded our ability to solve different well problems. The transient performances of MFHWs are calculated by many authors [19,20]. Then, Ozkan and Raghavan [21,22] (1991) used Laplace transformation into the source function, which can easily extend to the naturally fractured media and variable rate production condition. After that, relevant transient solutions of MFHWs are obtained [23,24,25]. Most of their works focus on the different fracture distributions and different reservoir types, while reservoir heterogeneity is less studied. To consider the Stimulated Reservoir Volume (SRV) in the unconventional reservoir, a composite model [5,26,27] is built by deriving point source solution from circular composite reservoir. However, the structure of MFHWs in composite circular reservoirs cannot match the linear flow around the SRV area and the reservoir heterogeneity around each fracture is not considered either.
For the heterogeneous reservoirs, Kikani [28] (1993) modeled the pressure performance of different sectional homogeneous reservoir by BEM (boundary element method). Idorenyin [29] (2016) considered the arbitrary-shaped composite reservoirs by extending Kikani’s model. Based on the element boundary method and Green’s function, Medeiros [30] (2010) presented a semi-analytical model for the pressure transient analysis of horizontal wells in heterogeneous reservoirs. Based on the fact that the main contribution for hydrocarbon flow rate of MFHWs is the SRV area [31], so this paper is to build a semi-analytical model for MFHWs in heterogeneous stimulated reservoir volume of naturally fractured media. To avoid fine grid of numerical reservoir simulator around fractures and wells, the plane source functions in the closed domain and boundary element method are applied, where stimulated reservoir volume is divided into compartmentalized heterogeneous blocks. In addition, the result is verified with other analytical methods. Finally, heterogeneous effects on pressure and production performance are presented.

2. Modeling of MFHWs

Firstly, we will present basic the physical model and assumption for hydraulic fractures and stimulated reservoir volume (SRV). Then, the Green’s function formulation of the pressure-transient solution for a locally homogeneous reservoir substructure will be given with inner and outer boundaries.

2.1. Physical Model and Assumptions

A Multiple fractured horizontal well passing through a heterogeneous reservoir with different properties is shown in Figure 1, where each block is the subsections of the reservoir characterized by uniform average properties. Except for the interfaces between the blocks, other boundaries of the blocks are assumed to be impermeable. The basic assumptions of the model are as follows:
(1)
The reservoir is heterogeneous along the horizontal wellbore with nf blocks of different length xei (i = 1, 2, 3, …, nf), and there is one hydraulic fracture distributed in each block. The thickness of the stimulated reservoir volume is uniform h and width is ye. The initial pressure is assumed pi.
(2)
The artificial hydraulic fractures are vertical to the horizontal wellbore, with fracture half-length Lfi, fracture height hfi, (i = 1, 2, 3, …, nf), and uniform flux distribution assumed in each hydraulic fracture.
(3)
Each block of unconventional reservoir was established by Warren and Root’s (1963) dual porosity model, where natural fractures and matrix follow the Darcy’s flow model. For each blocks, different natural fracture permeability kfi, crossflow coefficient λi, and elastic storage ratio ωi (i = 1, 2, 3, …, nf) can be defined.
(4)
The fluid assumed slightly compressible with constant viscosity and compressibility.
(5)
Gravity and capillary effect are negligible.

2.2. Pressure Solution for a Reservoir Substructure

Based on the Green’s function, the pressure transient solution of dimensionless space and time in porous media can be written as follows [18,31,32], where the domain boundary is divided into the inner and outer boundaries, represented by BwD, and BeD, respectively:
Δ p ( M D , t D ) = 0 t D B wD [ G ( M D , M D , t D τ ) p ( M D , τ ) n B w D p ( M D , τ ) G ( M D , M D , t D τ ) n B w D ] M D B w D d M D d τ 0 t D B eD [ G ( M D , M D , t D τ ) p ( M D , τ ) n B e D p ( M D , τ ) G ( M D , M D , t D τ ) n B e D ] M D B e D d M D d τ
where Δ p ( M D , t D ) is the pressure drop, G ( M D , M D , t D τ ) is the Green’s function, p ( M D , τ ) is the pressure. The B w D corresponds to well surface and B e D is the outer surface of the dimensionless domain, the points M D and M D represent the observation and source locations in dimensionless space, respectively, and tD is dimensionless time, and defined as:
t D = k ϕ C t μ 2 t
where, μ is the fluid viscosity, k is the permeability, is the characteristic reference length, C t is the total compressibility, ϕ is the porosity and t is time.
Then, the homogeneous Neumann boundary condition is used to eliminate partial derivative of Green function. The point flux term is obtained by Darcy’s law to represent the normal pressure derivatives on the boundary surfaces. Therefore, we can rewrite Equation (1) as follows:
Δ p ( M D , t D ) = 0 t D B wD q ˜ w ( M D , τ ) G ( M D , M D , t D τ ) d M D d τ       + 0 t D B eD q ˜ e ( M D , τ ) G ( M D , M D , t D τ ) d M D d τ
where,   q ˜ w ( M D , τ ) , q ˜ e ( M D , τ ) are the point flux of well surface and outer surface, respectively. We divide the inner and outer boundary surface into n and m segments of uniform flux, and the source function are defined as
B w D = j = 1 n B w D j j = 1 , 2 , 3 , , n
B e D = k = 1 m B e D k k = 1 , 2 , 3 , , m
S α ( M D , t D ) = B α G ( M D , M D , t D τ ) d M D for   α = w D j   or   e D k
Therefore, the Equation (3) can be rewritten as follows:
Δ p ( M D , t D ) = j = 1 n 0 t D q ˜ w j ( τ ) S w D j ( M D , t D τ ) d τ       + k = 1 m 0 t D q ˜ e k ( τ ) S e D k ( M D , t D τ ) d τ
Ozkan and Raghavan [21] applied the Laplace transformation into Equation (7), the variable flow rate condition and naturally fractured media can be easily implemented by Green’s function. Convert the convolution integrals into algebraic expressions. Then, Equation (7) can be rewritten as
Δ p ¯ ( M D , s ) = j = 1 n s × q ˜ ¯ w j ( s ) S ¯ w D j ( M D , s ) + k = 1 m s × q ˜ ¯ e k ( s ) S ¯ e D k ( M D , s )
The function with bar signs in Equation (8) indicates their function in Laplace transforms with respect to the time tD, and s is the Laplace-transform variable.

2.3. Source Functions for the Inner and Outer Boundaries

To obtain the transient pressure solution of fractured horizontal well in a rectangular parallelepiped, the source functions of inner and outer boundaries are evaluated by the Green’s function. For the partial and full penetration of the hydraulic fractures in the reservoir, 3D plane source and 2D line source should be considered.
Corresponding to the point source in the rectangular parallelepiped, the Green’s function in dual porous media can be expressed as [21]
G ¯ ( M D , M D , s ) = μ 2 k f y e D h D { ch ( u x D 1 ) + ch ( u x D 2 ) u sh ( u x e D )     + 2 k = 1 cos k π y D y e D cos k π y D y e D ch ( ε k x D 1 ) + ch ( ε k x D 2 ) ε k sh ( ε k x e D )     + 2 n = 1 cos n π z D h D cos n π z w D h D [ ch ( ε n x D 1 ) + ch ( ε n x D 2 ) ε n sh ( ε n x e D )     + 2 k = 1 cos k π y D y e D cos k π y D y e D ch ε k , n x D 1 + ch ε k , n x D 2 ε k , n sh ( ε k , n x e D ) ] }
In Equation (9), we have defined   M D = ( x D , y D , z D ) ,   M D = ( x D , y D , z D ) ,   x eD = x e ,   y eD = y e   h D = h , x D 1 = x e D | x D x D | ,   x D 2 = x e D | x D + x D | ,   ε n = u + ( n π / h D ) 2 ,   ε k = u + ( k π / y e D ) 2 ,   ε k , n = u + ( k π / y e D ) 2 + ( n π / h D ) 2 and u = s f ( s ) , f ( s ) = s ω ( 1 ω ) + λ λ + s ( 1 ω ) for naturally fractured reservoir.
The inner and outer boundaries consist of the planar surfaces in the rectangular parallelepiped, which represented hydraulic fractures and boundary planes. The source function for a plane is perpendicular to the x-axis, so a planar source in the y-z plane can be calculated by integration of Equation (9) (shown in Figure 2)
S ¯ ( M D , s ) = h w D / 2 h w D / 2 L f D L f D G ¯ ( M D , M D , s ) d y D d z D = μ L f D h w D k h y e D { ch ( u x D 1 ) - ch ( u x D 2 ) u s h ( u x e D ) + 2 y e D π L f D k = 1 1 k cos ( k π y D y e D ) cos ( k π y w D y e D ) sin ( k π L f D y e D ) ch ε k x D 1 - ch ε k x D 2 ε k s h ε k x e D } + μ π k h × 4 h D L f D y e D s n = 1 1 n cos ( n π z D h D ) sin ( n π z w D h D ) cos ( n π h w D 2 h D ) { ch ε n x D 1 - ch ε n x D 2 ε n s h ε n x e D + 2 y e D π L f D k = 1 1 k sin ( k π y D y e D ) cos ( k π y w D y e D ) sin ( k π L f D y e D ) ch ε n , k x D 1 ch ε n , k x D 2 ε n , k s h ε n , k x e D }
When the reservoir is fully penetrated by hydraulic fractures, the boundary surface is at full penetration, namely h w D = h D ; thus, Equation (10) can simplify to (the second item of Equation (10) is equal to 0)
S ¯ ( M D , s ) = μ L f D h w D k h y e D { ch ( u x D 1 ) - ch ( u x D 2 ) u s h ( u x e D ) + 2 y e D π L f D k = 1 1 k cos ( k π y D y e D ) cos ( k π y w D y e D ) sin ( k π L f D y e D ) ch ε k x D 1 - ch ε k x D 2 ε k s h ε k x e D }

3. Solution of Coupling Multiple Blocks

To obtain the pressure drop solution of MFHW in multiple reservoir blocks, we will demonstrate the coupling process of nf rectangular reservoir blocks in series in the x direction as shown in Figure 1. Then, the pressure solution process of the linear system of equations is presented.
The fractured horizontal well extends through a heterogeneous reservoir with rectangular parallelepiped blocks. All block boundary surfaces are assumed to be impermeable, except for the interface between different blocks. We will assume that the hydraulic fractures are uniform flux distribution. For the infinite conductivity condition, the pressure response can be achieved by subdividing each fracture or assumed by 0.732 point pressure of fracture, and finite conductivity fractures can be obtained by incorporating pressure drop with hydraulic fractures [24]. Here, each hydraulic fracture as one segment represents inner boundary of each block, while boundary surfaces are discretized into me segments. From Equation (8), the pressure drop at the center of the ith fracture can be given by
Δ p ¯ ( M D w i , s ) = { q ˜ ¯ w 1 ( s ) S ¯ w D 1 ( M D w 1 , s ) + k = 1 m e q ˜ ¯ e 1 , 2 k ( s ) S ¯ e D 1 , 2 k ( M D w 1 , s )    i = 1 q ˜ ¯ w i ( s ) S ¯ w D i ( M D w i , s ) + k = 1 m e q ˜ ¯ e i , 1 k ( s ) S ¯ e D i , 1 k ( M D w i , s ) + k = 1 m e q ˜ ¯ e i , 2 k ( s ) S ¯ e D i , 2 k ( M D w i , s ) i = 2 , 3 , , n f 1 q ˜ ¯ w n f ( s ) S ¯ w D n f ( M D w n f , s ) + k = 1 m e q ˜ ¯ e n f , 1 k ( s ) S ¯ e D n f , 1 k ( M D w n f , s ) i = n f
where M D w i indicate the center point of ith fracture, so q ˜ ¯ w i is the point M D w i flux, S ¯ e D i , 1 k ( M D w i , s ) and S ¯ e D i , 2 k ( M D w i , s ) are the source function of left and right surface at point M D w i in ith block, respectively.
The pressure drop at the kth segment of two interfaces of ith block, where the left face marked 1, and right face marked 2, can be given by
Δ p ¯ ( M D e i , 1 k , s ) = { q ˜ ¯ w i ( s ) S ¯ w D i ( M D e i , 1 k , s ) + k = 1 m e q ˜ ¯ e i , 1 k ( s ) S ¯ e D i , 1 k ( M D e i , 1 k , s ) + + k = 1 m e q ˜ ¯ e i , 2 k ( s ) S ¯ e D i , 2 k ( M D e i , 1 k , s ) ,     i = 2 , 3 , , n f 1 q ˜ ¯ w n f ( s ) S ¯ w D n f ( M D e n f , 1 k , s ) + k = 1 m e q ˜ ¯ e n f , 1 k ( s ) S ¯ e D n f , 1 k ( M D e n f , 1 k , s ) , i = n f
Δ p ¯ ( M D e i , 2 k , s ) = { q ˜ ¯ w 1 ( s ) S ¯ w D 1 ( M D e 1 , 2 k , s ) + k = 1 m e q ˜ ¯ e 1 , 2 k ( s ) S ¯ e D 1 , 2 k ( M D e 1 , 2 k , s ) ,     i = 1 q ˜ ¯ w i ( s ) S ¯ w D i ( M D e i , 2 k , s ) + k = 1 m e q ˜ ¯ e i , 1 k ( s ) S ¯ e D i , 1 k ( M D e i , 2 k , s ) + + k = 1 m e q ˜ ¯ e i , 2 k ( s ) S ¯ e D i , 2 k ( M D e i , 2 k , s ) ,     i = 2 , 3 , , n f 1  
where k = 1 ,   2 ,   3 , . , m e , M D e i , 1 k , M D e i , 2 k indicate the center point of kth segment of left and right surface in ith block, so q ˜ ¯ e i , 1 k   q ˜ ¯ e i , 2 k is the according flux of point M D e i , 1 k and   M D e i , 2 k , respectively.
Except the first and last block of SRV, each block has 1 + 2me points with unknown parameters of pressure and flux, therefore there are ( 1 + 2 m e ) × 2 × n f 4 m e unknown. Based on Equations (12)–(14), we yield n f + 2 ( n f 1 ) × m e linear equations. In order to get the same number of equations as the unknowns, the pressure and flux continuity condition can be applied.
(1)
Continuity conditions of pressure and flux at the interface
The pressure drop of the same interface should be equal for each segment k of different blocks, and the flux of interfaces should be continuous, so we have
{ Δ p ¯ ( M D e i , 2 k , s ) = Δ p ¯ ( M D e i + 1 , 1 k , s ) q ˜ ¯ ( M D e i , 2 k , s ) = q ˜ ¯ ( M D e i + 1 , 1 k , s ) i = 1 , 2 , 3 , , n f 1 k = 1 , 2 , 3 , , m e
(2)
Horizontal wellbore conditions
Assuming that the horizontal wellbore is infinite conductivity, the pressure drop of each fracture mid-point M D w i should be equal to the pressure drop of horizontal wellbore Δ p ¯ ( M D w , s ) , so we can have nf − 1 additional equations
Δ p ¯ ( M D w i , s ) = Δ p ¯ ( M D w , s ) i = 1 , 2 , 3 , , n f
Moreover, the sum of the fluxes from all hydraulic fractures should be equal to the total production rate qtotal.
i = 1 n f q ˜ ¯ w i 2 L f i h e i = q t o t a l / s
The linear system defined by Equations (12)–(17) has ( 1 + 2 m e ) × 2 × n f 4 m e equations for ( 1 + 2 m e ) × 2 × n f 4 m e unknowns, which can be written in the matrix-vector form AX = B. Here, take a multiple fractured horizontal well with nf fractures and boundary segment me = 1; for example, the components of coefficient matrix A is
A = [ s × S ¯ w 1 w 1 s × S ¯ e 1 , 2 w 1 0 0 0 0 0 1 s × S ¯ w 1 e 1 , 2 s × ( S ¯ e 1 , 2 e 1 , 2 + S ¯ e 2 , 1 e 2 , 1 ) s × S ¯ w 2 e 2 , 1 s × S ¯ e 2 , 2 e 2 , 1 0 0 0 0 0 s × S ¯ e 2 , 1 w 2 s × S ¯ w 2 w 2 s × S ¯ e 2 , 2 w 2 0 0 0 1 0 s × S ¯ e 2 , 1 e 2 , 2 s × S ¯ w 2 e 2 , 2 s × ( S ¯ e 2 , 2 e 2 , 2 + S ¯ e 3 , 1 e 3 , 1 ) s × S ¯ w 3 e 3 , 1 0 0 0 0 0 0 s × S ¯ e 3 , 1 w 3 s × S ¯ w 3 w 3 0 0 1 0 0 0 0 0 s × S ¯ e n f , 1 w n f s × S ¯ w n f w n f 1 2 L f 1 h 1 s 0 2 L f 2 h 3 s 0 2 L f 3 h 3 s 2 L f n f h e n f s 0 0 ] ( 2 n f ) × ( 2 n f )
where to simplify writing, S ¯ D w i D w i   =   S ¯ D w i ( M D w i , s ) , S ¯ e D i , 2 k D w i   =   S ¯ e D i , 2 k ( M D w i , s ) . The solution vector, X, has the following components:
X = [ q ˜ ¯ w 1 , q ˜ ¯ e 1 , 2 , q ˜ ¯ w 2 , q ˜ ¯ e 2 , 2 , q ˜ ¯ w 3 , q ˜ ¯ e 3 , 2 , , q ˜ ¯ w ( n f 1 ) , q ˜ ¯ e ( n f 1 ) , 2 , q ˜ ¯ w n f , Δ p ¯ w D ] 2 n f × 1 T
and the components of the right-hand side vector, B, are
B = [ 0 , 0 , 0 , 0 , , 0 , q t o t a l ] ( 2 n f ) × 1 T
The pressure and flux solution of this linear system can be solved. Then, the solution in the Laplace transform domain can be inverted into real-time domain by the Stehfest Algorithm [33]. The dimensionless pressure and dimensionless fluxes of ith fracture are defined as:
p D = 2 π k h μ q t o t a l Δ p w D ;   q f i = 2 L f D i h w D i q ˜ ¯ w i ;   q f D i = 2 L f D i h w D i q ˜ ¯ w i q t o t a l , i = 1 , 2 , , n f .

4. Comparison and Verification

In order to verify the accuracy of our model, a relatively simple model of MFHW is discussed, which located in the center of a closed, homogenous- and single-porosity media, rectangular reservoir with three blocks. The transient pressure performance of MFHW in homogeneous reservoirs is calculated by [24] using the source function and the principle of superposition method. The properties of the hydraulic fractures and reservoir are given in Table 1, and the results with different discretized number me (Figure 3) are shown in Table 2 and plotted in Figure 4.
Table 2 presents the comparisons of our models with Zerzar’s model about pressure and pressure derivatives. Four cases of different interface discretization between the blocks are shown in Figure 3. Table 2 demonstrates that there are almost no differences among the four cases. All of them can match the Zerzar’s model very well. The maximum error (me = 1) between the Zerzar’s model and our model is 0.0664% of dimensionless pressure drop, 0.2043% of dimensionless pressure derivative, which is too small to be distinguished on the log-log plot of Figure 4, which shows high correctness of our model. To simplify the calculation process, we assumed me = 1 for the following discussion.
In Figure 4, the flow stages of MFHW can be identified including first linear flow, first radial flow and late pseudo-steady boundary flow, which are different from traditional MFHWs flow regimes. As hydraulic fractures are close to the impermeable boundaries in our model, the second linear and radial flows do not appear. In addition, we compare our model with trilinear flow model [6], which is widely used for unconventional reservoir because of its concise formula and convenient calculation. In Figure 4, it shows that the pressure and pressure derivative of our model can be matched perfectly with trilinear model at first linear flow and late pseudo-steady boundary flow period. However, radical flow around each hydraulic fracture is missing in trilinear flow model, which is not in line with real situation. Therefore, our model is more accurate than trilinear flow model.

5. Results and Discussion

The semi-analytical pressure solution for MFHWs within some rectangular reservoir blocks of different parameters is obtained. This section presents pressure behavior and flux distribution of each hydraulic fracture with some heterogeneous parameters. We assume that there are three compartments of the reservoir, the horizontal well penetrates these three compartments and there are three hydraulic fractures in the middle of each block. Besides the basic properties given in Table 1, the parameters in dual porosity media includes elastic storage ratio ω f = 0.1 and crossflow coefficient λ f = 2 . The heterogeneous parameter ratio values of three blocks are shown in Table 3.

5.1. Heterogeneous Permeability

We consider five cases of heterogeneous permeability distribution, and the natural fracture permeability ratio of these three blocks is shown in Table 3. In cases 2 and 4, the permeability increases in second and third compartment, while in cases 3 and 5, the permeability decreases in second and third compartment. Figure 5 shows the pressure and pressure derivative responses as well as the flux distribution of three different hydraulic fractures.
In Figure 5a, for the pressure-transient responses, we can see that the heterogeneous permeability distribution of different drainage blocks affects both arriving time and the pressure drop of the whole MFHW flow regimes including first linear flow, crossflow regime, first radial flow and late pseudo-steady boundary flow. The higher the permeability of cases 2 and 4, the earlier different regimes arrive, and the less the pressure drop compared with reference case 1 (homogeneous permeability), which is beneficial to production. On the contrary, with the lower permeability of cases 3 and 5, the higher-pressure transient curve means more pressure drop demands for constant flow rate condition of horizontal wellbore. Figure 5b shows the effect of heterogeneous permeability on the flux of each fracture. As shown in Figure 5b, there are great differences between different fracture fluxes with unequal permeability distribution. In homogeneous permeability case 1, the flux of each fracture is equal to qtotal/3. For the other heterogeneous permeability cases, the flux of fractures in higher permeability blocks is larger during the early linear flow period. Then, the difference will decrease because of crossflow between each block. Finally, the differences will become stable in the last pseudo-steady boundary period. Simultaneously, for lower-permeability case 5, the crossflow period appears earlier compared to the higher-permeability case 4, which is the same as the pressure response in Figure 5a.

5.2. Non-Uniform Length of Blocks

To consider the non-uniform length distribution effect, we use a heterogeneous permeability model with kf1:kf2:kf3 = 1:2:1. Three different cases with unequal length distributions of each compartment are considered in Table 3 under the constant whole drainage length condition. Then, the pressure response and flux distribution are shown in Figure 6.
From Figure 6a, we can see that if the whole drainage length is constant, the dimensionless pressure drop and pressure derivative curves have no change when the drainage length of each fracture is different. However, in Figure 6b, the flux distribution of each fracture is different after considering the block length heterogeneity. Non-uniform block length only affects the final flux proportion of each fracture. It has no effect on the arriving time to the crossflow and last pseudo-steady flow period. In other words, there is no obvious effect on the wellbore pressure drop when the partial section area is damaged or stimulated in the same drainage area, which only affect the flux proportion of each fracture.

5.3. Different Elastic Storage Ratio

To consider the different elastic storage ratio effect of each block in drainage area, three cases are considered in Table 3, and corresponding pressure and flux curves are shown in Figure 7.
Figure 7a illustrates that the elastic storage ratio between different drainage blocks mainly affects the first linear flow, first radial flow and crossflow regimes. When the elastic storage ratio increases from case 1 to case 2 and case 3, which means the storage capacity of natural fracture system is greater and more fluid is stored in the natural fractures, the dimensionless pressure derivative curves show smaller groove because of less fluid transfer between natural fractures and matrix. From Figure 7b, it shows the same performance as pressure response. The flux distribution is only affected in the early stages, while the flux of each fracture converges to the same qtotal/3 at the later pseudo-steady flow period. Similarly, the less natural fracture storage capacity ω, the less corresponding flux of hydraulic fracture.

5.4. Different Crossflow Ratio

Another important parameter of dual porosity for naturally fracture system is crossflow ratio. Three cases are investigated in Table 3, and the pressure responses and flux distribution are shown in Figure 8.
The crossflow ratio represents the fluid flowing ability from the matrix to the natural fracture system. As shown in Figure 8, the different crossflow ratio only affects crossflow stage in the pressure responses and flux distribution. The larger λ value of the block is, the earlier the crossflow appears between different cases and the higher the flux is during the crossflow period.

5.5. Non-Uniform Half-Length Distribution of Fractures

With the constant total fracture length, we consider three cases of non-uniform fracture half-length distribution. The ratio of three hydraulic fracture half-lengths is shown in Table 3. Figure 9 gives the pressure response and flux distribution.
From Figure 9a, we can see that if the total fracture half-length is constant, pressure response curves have nothing to do with different half-length distribution. However, In Figure 9b, the flux distribution of each fracture is changed by half-length heterogeneity. The flux difference between each fracture becomes smaller as the flow time increases, and becomes stable gradually after crossflow stages.

6. Conclusions

This work investigates the pressure transient solutions of MFHWs in heterogeneous systems where the sections of the naturally fractured reservoir display different characteristics. A semi-analytical method is presented which combines source function and boundary element method. Major conclusions are summarized below:
(1)
This model only discretizes hydraulic fractures and interface boundaries. With a small number of segments interface discretization, it can match verified model very well on the log–log plots. As a result, this model can greatly reduce fine gridding required by numerical approach.
(2)
Compared to the trilinear flow model, this model can consider the radial flow around the hydraulic fracture, and the interference effect between hydraulic fractures with different formation properties.
(3)
As the permeability increases along the horizontal well direction, the pressure drop decreases and flow regimes become earlier. The fractures flux fluctuates with flow regimes, where the difference between high- and low-permeability blocks decreases firstly, then increases gradually, and tends to be stable finally. When the entire drainage length keeps constant, the pressure responses have nothing to do with non-uniform heterogeneous block length while flux distribution has been changed obviously.
(4)
The elastic storage ratio mainly affects the early flow regimes including first linear flow, first radial flow and crossflow periods. The bigger the elastic storage ratio, the smaller the groove and the greater the flux. However, the crossflow ratio only affects the crossflow period. The bigger the crossflow ratio, the earlier the crossflow arrives.

Author Contributions

D.F. and J.Y. have conceived the analytical method. D.F. and H.S. have created the heterogeneous model by the Green method. D.F. and H.Z. have written the original manuscript. X.Y. has validated the outcome with numerical model. Z.S. has provided analysis tools. D.F., H.S. and H.Z. have analyzed the results and have made the revisions. Z.S., J.Y. and H.S. have provided the funding. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are thankful for financial support from the National Natural Science Foundation of China (Grant No. 51774308, No. 51774317), and the Science and Technology Support Plan for Youth Innovation of University in Shandong Province (Grant No. 2019KJH002).

Acknowledgments

The authors would also like to appreciate the reviewers and editors whose critical suggestions have been very helpful to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Acoefficient matrix
Bright-hand side vector
BWDThe dimensionless inner surface domain
BeDThe dimensionless outer surface domain
Ctthe total compressibility, MPa−1
GGreen’s function
hthe reservoir thickness, m
kpermeability, mD
kfithe ith fracture permeability, mD
Lfithe ith fracture half-length, m
lThe characteristic reference length, m
MDThe observation locations in dimensionless space
M D The source locations in dimensionless space
meinterface segment number
nfthe hydraulic fractures number
piinitial reservoir pressure, MPa
ΔpPressure drop, MPa
pDthe dimensionless pressure value
qtotalthe total wellbore flow rate, m3/d
tDthe dimensionless time
XThe solution vector
xethe length of blocks
yethe width of reservoir
λiThe ith block crossflow coefficient
ωiThe ith block elastic storage ratio
μViscosity
ϕPorosity

Subscripts and Superscripts

Ddimensionless
eboundary
iserial number of hydraulic fractures
ωwellbore

References

  1. Soliman, M.; Daal, J.; East, L. Fracturing unconventional formations to enhance productivity. J. Nat. Gas Sci. Eng. 2012, 8, 52–67. [Google Scholar] [CrossRef]
  2. Nguyen-Le, V.; Shin, H. Development of reservoir economic indicator for Barnett Shale gas potential evaluation based on the reservoir and hydraulic fracturing parameters. J. Nat. Gas Sci. Eng. 2019, 66, 159–167. [Google Scholar] [CrossRef]
  3. Dian, F.; Amin, E. Semi-analytical modeling of shale gas flow through fractal induced fracture networks with micro-seismic data. Fuel 2017, 193, 444–459. [Google Scholar]
  4. Fan, D.Y.; Yao, J.; Sun, H.; Zeng, H.; Wang, W. A composite model of hydraulic fractured horizontal well with stimulated reservoir volume in tight oil & gas reservoir. J. Nat. Gas Sci. Eng. 2015, 24, 115–123. [Google Scholar]
  5. Zhao, Y.; Zhang, L.-H.; Luo, J.-X.; Zhang, B.-N. Performance of fractured horizontal well with stimulated reservoir volume in unconventional gas reservoir. J. Hydrol. 2014, 512, 447–456. [Google Scholar] [CrossRef]
  6. Brown, M.; Ozkan, E.; Ragahavan, R.; Kazemi, H. Practical solutions for pressure transient responses of fractured horizontal wells in unconventional reservoirs. SPE Res. Eval. Eng. 2011, 14, 663–676. [Google Scholar] [CrossRef]
  7. Stalgorova, E.; Mattar, L. Practical analytical model to simulate production of horizontal wells with branch fractures. In Proceedings of the SPE Paper 162515 Presented at SPE Canadian Unconventional Resources Conference, Calgary, AB, Canada, 30 October–1 November 2012. [Google Scholar]
  8. Stalgorova, K.; Mattar, L. Analytical model for unconventional multi-fractured composite systems. SPE Reserv. Eval. Eng. 2013, 16, 246–256. [Google Scholar] [CrossRef]
  9. Guo, J.; Zeng, J.; Wang, X. Analytical model for multi-fractured horizontal wells in heterogeneous shale reservoirs. In Proceedings of the SPE Asia Pacific Oil & Gas Conference and Exhibition, Society of Petroleum Engineers (SPE), Perth, Australia, 25–27 October 2016. [Google Scholar]
  10. Ozcan, O.; Sarak, H.; Ozkan, E.; Raghavan, R. A trilinear flow model for a fractured horizontal well in a fractal unconventional reservoir. In Proceedings of the SPE Paper 170971 Presented at SPE Annual Technical Conference and Exhibition, Amsterdam, The Netherlands, 27–29 October 2014. [Google Scholar]
  11. Wang, W.; Su, Y.; Sheng, G.; Cossio, M.; Shang, Y. A mathematical model considering complex fractures and fractal flow for pressure transient analysis of fractured horizontal wells in unconventional reservoirs. J. Nat. Gas Sci. Eng. 2015, 23, 139–147. [Google Scholar] [CrossRef]
  12. Salam, A.R.; Tiab, D. Partially penetrating hydraulic fractures: Pressure responses and flow dynamics. In Proceedings of the SPE Paper 164500 Presented at SPE Production and Operation Symposium, Oklahoma City, OK, USA, 23–26 March 2013. [Google Scholar]
  13. Zawila, J.; Fluckiger, S.; Hughes, G.; Kerr, P.; Hennes, A.; Hofmann, M. An integrated, multi-disciplinary approach utilizing stratigraphy, petro-physics, and geophysics to predict reservoir properties of tight unconventional sandstones in the Powder River basin, Wyoming, USA. In Proceedings of the SEG Paper 5852423 Presented at SEG Annual Meeting, New Orleans, LA, USA, 19–20 August 2015. [Google Scholar]
  14. Liu, P.; Ju, Y.; Ranjith, P.G.; Zheng, Z.; Chen, J. Experimental investigation of the effects of heterogeneity and geostress difference on the 3D growth and distribution of hydro-fracturing cracks in unconventional reservoir rocks. J. Nat. Gas Sci. Eng. 2016, 35, 541–554. [Google Scholar] [CrossRef]
  15. Karimi-Fard, M.; Durlofsky, L.J.; Aziz, K. An efficient Discrete-Fracture model applicable for general-purpose reservoir simulators. In Proceedings of the SPE Paper 79699 Presented at SPE Reservoir Simulation Symposium, Houston, TX, USA, 3–5 February 2004. [Google Scholar]
  16. Warpinski, N.; Mayerhofer, M.; Vincent, M.; Cipolla, C.; Lolon, E. Stimulating Unconventional Reservoirs: Maximizing Network Growth While Optimizing Fracture Conductivity. J. Can. Pet. Technol. 2009, 48, 39–51. [Google Scholar] [CrossRef]
  17. Yan, X.; Huang, Z.; Yao, J.; Li, Y.; Fan, D. An efficient embedded discrete fracture model based on mimetic finite difference method. J. Pet. Sci. Eng. 2016, 145, 11–21. [Google Scholar] [CrossRef]
  18. Gringarten, A.C.; Ramey, H.J. The Use of Source and Green’s Functions in Solving Unsteady-Flow Problems in Reservoirs. Soc. Pet. Eng. J. 1973, 13, 285–296. [Google Scholar] [CrossRef]
  19. Chen, C.-C.; Rajagopal, R. A Multiply-Fractured Horizontal Well in a Rectangular Drainage Region. SPE J. 1997, 2, 455–465. [Google Scholar] [CrossRef]
  20. Guo, G.; Evans, R.D. Pressure-transient behavior and inflow performance of horizontal wells intersecting discrete fractures. In Proceedings of the SPE paper 26446 Presented at SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 3–6 October 1993. [Google Scholar]
  21. Ozkan, E.; Raghavan, R. New Solutions for Well-Test-Analysis Problems: Part 1-Analytical Considerations (includes associated papers 28666 and 29213 ). SPE Form. Eval. 1991, 6, 359–368. [Google Scholar] [CrossRef] [Green Version]
  22. Özkan, E.; Raghavan, R. New Solutions for Well-Test-Analysis Problems: Part 2 Computational Considerations and Applications. SPE Form. Eval. 1991, 6, 369–378. [Google Scholar] [CrossRef]
  23. Jia, P.; Cheng, L.; Huang, S.; Wu, Y. A semi-analytical model for the flow behavior of naturally fractured formations with multi-scale fracture networks. J. Hydrol. 2016, 537, 208–220. [Google Scholar] [CrossRef]
  24. Zerzar, A.; Bettam, Y. Interpretation of multiple hydraulically fractured horizontal wells in closed systems. In Proceedings of the SPE Paper 84888 Presented at SPE International Improved Oil Recovery Conference in Asia Pacific, Kuala Lumpur, Malaysia, 8–10 June 2003. [Google Scholar]
  25. Zhao, Y.; Zhang, L.-H.; Zhao, J.-Z.; Luo, J.-X.; Zhang, B.-N. “Triple porosity” modeling of transient well test and rate decline analysis for multi-fractured horizontal well in shale gas reservoirs. J. Pet. Sci. Eng. 2013, 110, 253–262. [Google Scholar] [CrossRef]
  26. Xu, J.; Guo, C.; Wei, M.; Jiang, R. Production performance analysis for composite shale gas reservoir considering multiple transport mechanisms. J. Nat. Gas Sci. Eng. 2015, 26, 382–395. [Google Scholar] [CrossRef]
  27. Zeng, H.; Fan, D.; Yao, J.; Sun, H. Pressure and rate transient analysis of composite shale gas reservoirs considering multiple mechanisms. J. Nat. Gas Sci. Eng. 2015, 27, 914–925. [Google Scholar] [CrossRef]
  28. Kikani, J.; Horne, R.N. Modeling Pressure-Transient Behavior of Sectionally Homogeneous Reservoirs by the Boundary-Element Method. SPE Form. Eval. 1993, 8, 145–152. [Google Scholar] [CrossRef]
  29. Idorenyin, E.H.; Shirif, E. Transient response in arbitrary-shaped composite reservoirs. SPE Res. Eval. Eng. 2016, 9, 1–13. [Google Scholar] [CrossRef]
  30. Medeiros, F.; Özkan, E.; Kazemi, H. A Semianalytical Approach to Model Pressure Transients in Heterogeneous Reservoirs. SPE Reserv. Eval. Eng. 2010, 13, 341–358. [Google Scholar] [CrossRef]
  31. Mayerhofer, M.J.; Lolon, E.; Warpinski, N.R.; Cipolla, C.L.; Walser, D.W.; Rightmire, C.M. What Is Stimulated Reservoir Volume? SPE Prod. Oper. 2010, 25, 89–98. [Google Scholar] [CrossRef]
  32. Carslaw, H.S.; Yaeger, J.C. Conduction of Heat in Solids; Oxford University Press: Oxford, UK, 1959. [Google Scholar]
  33. Stehfest, H. Numerical Inversion of Laplace Transforms. Commun. ACM 1970, 13, 47–49. [Google Scholar] [CrossRef]
Figure 1. Multiple fractured horizontal wells (MFHWs) through a compartmentalized heterogeneous reservoir.
Figure 1. Multiple fractured horizontal wells (MFHWs) through a compartmentalized heterogeneous reservoir.
Energies 13 05204 g001
Figure 2. A plane perpendicular to the x-axis in the rectangular parallelepiped. (a) One block with hydraulic fracture; (b) the fracture plane.
Figure 2. A plane perpendicular to the x-axis in the rectangular parallelepiped. (a) One block with hydraulic fracture; (b) the fracture plane.
Energies 13 05204 g002
Figure 3. Different discretization schemes at the block interfaces.
Figure 3. Different discretization schemes at the block interfaces.
Energies 13 05204 g003
Figure 4. Model verification of Zerzar model and our model (nf = 3, me = 1).
Figure 4. Model verification of Zerzar model and our model (nf = 3, me = 1).
Energies 13 05204 g004
Figure 5. Pressure response and flux distribution of heterogenous permeability. (a) Pressure and pressure derivative; (b) flux distribution of each hydraulic fracture.
Figure 5. Pressure response and flux distribution of heterogenous permeability. (a) Pressure and pressure derivative; (b) flux distribution of each hydraulic fracture.
Energies 13 05204 g005
Figure 6. Pressure response and flux distribution of non-uniform length distribution. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Figure 6. Pressure response and flux distribution of non-uniform length distribution. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Energies 13 05204 g006
Figure 7. Pressure response and flux distribution of different elastic storage ratio. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Figure 7. Pressure response and flux distribution of different elastic storage ratio. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Energies 13 05204 g007
Figure 8. Pressure response and flux distribution of different crossflow ratio. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Figure 8. Pressure response and flux distribution of different crossflow ratio. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Energies 13 05204 g008
Figure 9. Pressure response and flux distribution of different half-length distribution. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Figure 9. Pressure response and flux distribution of different half-length distribution. (a) Pressure and pressure derivative; (b) Flux distribution of each hydraulic fracture.
Energies 13 05204 g009
Table 1. Values of well and reservoir parameters.
Table 1. Values of well and reservoir parameters.
ParametersValueUnit
Reference length, 100m
Hydraulic fracture number3-
Reservoir size in x, y-direction xe, ye1200m
Each block size in x, y direction, xei, yei400m
Half-length of each hydraulic fracture, Lfi60m
Formation thickness, h100m
hydraulic fracture height, hfi 100m
Natural fracture permeability, kfi100mD
Matrix Permeability, km1mD
Initial pressure35MPa
elastic storage ratio ωfi0.1-
crossflow coefficient λfi2-
Table 2. Comparison of the model results with different segments of block interfaces.
Table 2. Comparison of the model results with different segments of block interfaces.
tDZerar 2004me = 1me = 2me = 3me = 4
pDdpD/dlntDpDdpD/dlntDpDdpD/dlntDpDdpD/dlntDpDdpD/dlntD
0.0010.0311390.0164650.0311390.0164650.0311390.0164650.0311390.0164650.0311390.016465
0.010.0984650.0520510.0984650.0520510.0984650.0520510.0984650.0520510.0984650.052051
0.10.2987810.1388950.2987810.1388950.2987810.1388950.2987810.1388950.2987810.138895
10.6442860.1840870.6440490.1832290.6440490.1832290.6440490.1832290.6440490.183229
101.2890510.509751.2855850.5096781.2855850.5096781.2855850.5096781.2855850.509678
1005.2373774.363735.2376414.3662735.2376414.3662735.2376414.3662735.2376414.366273
100044.519843.7094144.5338343.7808844.5338343.7808844.5338343.7808844.5338343.78088
10,000437.4155437.3679437.1711435.6028437.1711435.6028437.1711435.6028437.1711435.6028
Errors--0.0664%0.2043%0.0654%0.2003%0.065%0.2001%0.065%0.2001%
Table 3. The heterogeneous parameter ratio of three blocks.
Table 3. The heterogeneous parameter ratio of three blocks.
ParametersCase 1Case 2Case 3Case 4Case 5
k f 1 : k f 2 : k f 3 1:1:11:2:11:1/2:11:2:31:1/2:1/4
x e 1 : x e 2 : x e 3 1:1:11:2:11:1/2:1--
ω f 1 : ω f 2 : ω f 3 1:1:11:2:11:2:3--
λ f 1 : λ f 2 : λ f 3 1:1:11:2:11:2:3--
L f 1 : L f 2 : L f 3 1:1:11:2:12:1:2--

Share and Cite

MDPI and ACS Style

Fan, D.; Sun, H.; Yao, J.; Zeng, H.; Yan, X.; Sun, Z. Pressure-Transient Performances of Fractured Horizontal Wells in the Compartmentalized Heterogeneous Unconventional Reservoirs. Energies 2020, 13, 5204. https://doi.org/10.3390/en13195204

AMA Style

Fan D, Sun H, Yao J, Zeng H, Yan X, Sun Z. Pressure-Transient Performances of Fractured Horizontal Wells in the Compartmentalized Heterogeneous Unconventional Reservoirs. Energies. 2020; 13(19):5204. https://doi.org/10.3390/en13195204

Chicago/Turabian Style

Fan, Dongyan, Hai Sun, Jun Yao, Hui Zeng, Xia Yan, and Zhixue Sun. 2020. "Pressure-Transient Performances of Fractured Horizontal Wells in the Compartmentalized Heterogeneous Unconventional Reservoirs" Energies 13, no. 19: 5204. https://doi.org/10.3390/en13195204

APA Style

Fan, D., Sun, H., Yao, J., Zeng, H., Yan, X., & Sun, Z. (2020). Pressure-Transient Performances of Fractured Horizontal Wells in the Compartmentalized Heterogeneous Unconventional Reservoirs. Energies, 13(19), 5204. https://doi.org/10.3390/en13195204

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

Article Metrics

Back to TopTop