# Productivity Analysis of Volume Fractured Wells under Different Working Systems

Authors: Hongfei Ma; Wenqi Zhao; Meng Sun; Xiaodong Wang; Lun Zhao; Chunmei Zou; Bo Wang
Journal: Geofluids
 (2021)

---

## Abstract

The volume fracturing technique has been widely used to improve the productivity of ultralow-permeability reservoirs. This paper presents a new semianalytical model to simulate the pressure transient and production behaviour of finite conductivity vertical fractured wells with stimulated reservoir volume (SRV) in heterogeneous reservoirs. The model is based on the five-linear flow model, the Warren-Root model, and fracture conductivity influence function. The model is validated by comparing its results with a numerical model. One novelty of this model is its consideration of three different kinds of production prediction models. Constant rate, constant pressure, and compound working systems are taken into account. This paper illustrates the effects of the SRV size and shape, mobility ratio, initial flow rate, limiting wellbore pressure, and hydraulic fracture parameters under different working systems. Results show that the SRV and parameters of fractures have a significant influence on long-term well performance. Moreover, the initial rate can extend the constant rate period by 418%, and limiting wellbore pressure can effectively improve the cumulative recovery rate by 23%. Therefore, this model can predict long-term wells’ behaviour and provide practical guiding significance for hydraulic fracturing design.


---

## Body

## 1. Introduction

With the increasing demand of the energy market and the incredible processing of hydraulic fracturing technology, the development of ultralow-permeability oil and gas reservoirs becomes increasingly critical. The hydraulic fracturing technique has been widely used to improve the productivity of tight reservoirs. Induced fractures will generate around hydraulic fractures, and the stimulated reservoir volume (SRV) will reduce the seepage resistance, change the flow model, and increase single well production [1]. However, the problems of unsteady flow caused by vertical fractures with finite conductivity are very complicated. It is usually necessary to deal with the coupled integral equation of formation flow and fracture flow. The compound linear flow model can avoid this problem and get an approximate solution, especially for the closed reservoir of narrow channel type.Gringarten et al. [2] studied the pressure distribution of vertically fractured wells with infinite fracture conductivity in infinite homogeneous reservoirs, ignoring the influence of finite conductivity fractures. Lee and Brockenbrough [3] further developed an approximate solution for the trilinear flow of finite conductivity vertical fractured wells in infinite homogeneous reservoirs. Cinco et al. [4–7] simplified the trilinear flow model to the bilinear flow model, suitable for finite conductivity fractures in dual-porosity reservoirs. Houze et al. [8] introduced a solution for a reservoir with a fracture in a dual-porosity reservoir, which applies to the uniform flux and the infinite conductivity fracture. Brown et al. [9] proposed a new trilinear flow model where the SRV region was modelled by the Warren-Root [10] and Kazemi [11] models, but the enhanced region occupies the whole space between the fractures. Du et al. [12] described the SRV area based on Kazemi’s dual-porosity medium model. Brohi et al. [13] proposed an analytical trilinear flow model to describe the system’s pressure distribution. The SRV region occupies the entire space near the fractures in a heterogeneous reservoir. Stalgorova and Mattar [14, 15] improved the trilinear flow model to a five-region flow by simplifying the SRV into an enhanced area with limited width. However, the effect of the dual-porosity medium has not been considered. In recent years, Liu et al. [16, 17] used the fractal permeability composite reservoir model to evaluate tight reservoir productivity. Wang et al. [18] presented a semianalytical model by only taking the hydraulic fractures and SRV into account. Yuan et al. [19] proposed a new analytical solution model to evaluate the transient linear flow in heterogeneous reservoirs. However, few studies consider the production behaviour of finite conductivity vertical fractured wells with SRV in heterogeneous reservoirs despite the immense research efforts.In recent years, there has been an increasing amount of literature on finite conductivity fractures. Al-Kobaisi et al. [20] employed a numerical method to obtain the pressure dynamics of finite conductivity fractures. However, the solutions are complex and need a large amount of computation unfavourable for wide application. Blasingame and Poe [21] regarded fracture conductivity as a kind of skin that can cause an extra pressure drop and so proposed the hybrid solution of a trilinear flow model (equivalent to the finite conductivity function) and an infinite conductivity model. This, however, is challenging to simulate the formation linear flow stage of finite conductivity fracture transits from bilinear flow to pseudoradial flow. Wilkinsons [22] solved the early flow characteristic function (fracture linear flow and bilinear flow) of low conductivity fractures by employing the Fourier transform method to deal with this problem. Wang et al. [23] gave the conductivity influence function in the Laplace domain, which simplified the calculation and was well applied. Based on the research results of Wang et al. [24], Wang et al. [25] calculated the semianalytical pressure solution of finite conductivity fracture by using the boundary element numerical method (BEM) and gave the finite conductivity influence function of quasi-steady state by data regression. Wang et al. [26] calculated the productivity of horizontal gas wells based on the finite conductivity function. Luo et al. [27] made a detailed analysis of the application range and calculation error of the conductivity influence function, but it was calculated in the Laplace domain. In this paper, an improved conductivity influence function in the time domain is derived for studying the hydraulic fracture behaviour, which possesses a higher degree of accuracy and can be easily applied.The object of this paper is to establish a semianalytical model of ultralow-permeability reservoir with SRV. The reservoir is divided into several subsystems. The SRV region will be represented by Warren-Root dual-porosity model. The flow in the hydraulic fractures is calculated by the conductivity influence function. Expressions for both the bottom hole pressure function and rate function in the Laplace domain will be provided. The influence of sensitive parameters, such as the SRV size and shape, the permeability of SRV region, fracture conductivity, initial flow rate, and conversion pressure, is analysed under three different working systems, including the compound system, the constant flow rate system, and the constant flowing pressure system.The compound working system developed in this work possesses a high degree of accuracy and enables long-term well performance prediction to be made quickly. The compound working system implies that production begins at a constant flow rate, and then, the production rate drops as well produces against constant pressure [28]. Numerous mathematical models have been adopted to investigate the flow in a single working system. However, the oilfield’s actual production is usually under compound working system, and the dynamic model under a single working system cannot correctly predict well’s production performance.


## 2. Physical Model and Mathematical Model

### 2.1. Physical Model

The fracture network formation near the wellbore is mainly related to reservoir pressure, mechanical rock properties, and natural fractures. The weak structural surfaces between rocks cause volume fractures such as natural fractures, joints, and bedding or weak points in the formation. Furthermore, the rock brittleness index is the physical property basis for forming a fracture network. Modern hydraulic fracturing technology can form a particular enhanced region with induced fractures near hydraulic fractures by controlling the injection and reverse discharge flow [29]. In this way, the near-wellbore zone is a fracture network formed by the interweaving of primary fractures and multistage secondary fractures. The area of the secondary fracture is also known as the stimulated reservoir volume (SRV). The SRV region is simulated through a Warren-Root dual-porosity medium model.Figure1 shows the schematic diagram of the five-linear flow physical model. In the centre of a closed rectangular drainage area, there is a production well with hydraulic fractures that vertically penetrate through the reservoir. Only a quarter of the flow area will be considered because of the symmetry of the formation. The flow areas are divided into Regions 1, 2, 3, 4, and 5. Flow in the model is regarded as a combination of five linear flows within contiguous regions. Each region comprises one-dimensional linear flow, and the boundaries between regions are continuous in flux and pressure. To keep the solution general, each region has various parameters (permeability, porosity, and total compressibility). Region 4 represents the SRV region near the hydraulic fracture, while Region 5 represents the finite conductivity hydraulic fractures.Figure 1
(a) Schematic of the five-linear flow model representing five continuous flow regions for a vertical well with an SRV region. (b) The matrix of SRV regions. (c) Classic Warren-Root dual-porosity medium model.
(a)(b)(c)The basic assumptions of the model are as follows:(i)
The hydraulic fracturing vertical well is located in the centre of a rectangular low permeability oil reservoir with closed boundaries(ii)
The thickness of the reservoir ish, and the uniform initial pressure is pi(iii)
The hydraulic fractures thoroughly penetrate the formation(iv)
The effect of gravity and capillary forces on fluid flow in each region is negligible(v)
Isothermal single-phase fluid flow is assumed


### 2.2. Mathematical Model

The pressure distribution in each region is controlled by the flow equation given by(1)∇2p=ϕμck∂p∂t.To simplify the solution, this equation for each region was rewritten in dimensionless terms. The definitions of dimensionless parameters are given in AppendixA, and the derivations of equations and solutions are given in Appendix B.


#### 2.2.1. Governing Equations in Region 1

Region 1 produces one-dimensional flow in they-direction. This visualization of the flow behaviour is reasonable for fluid flow before the pseudo-steady-state flow period. Considering the symmetry, the flow governing differential equations and boundary conditions are as follows:
(2)∂2p1D∂yD2=1η1D∂p1D∂tD,xwD+xmD≤xD≤xeD,ywD+1≤yD≤yeD.The initial condition is(3)p1DxD,yD,0=0.The outer boundary condition is(4)∂p1DxD,yeD,tD∂yD=0.The interface conditions of Regions 1 and 3 are(5)p1DxD,ywD+1,tD=p3DxD,ywD+1,tD,(6)∂p1DxD,ywD+1,tD∂yD=∂p3DxD,ywD+1,tD∂yD.


#### 2.2.2. Governing Equations in Region 2

In Region 2, the flow is assumed to be one-dimensional and in they-direction. Considering the symmetry, the dimensionless governing equations are as follows xwmD=xwD+xmD:
(7)∂2p2D∂yD2=1η2D∂p2D∂tD,xwD≤xD≤xwmD,ywD+1≤yD≤yeD.The initial condition is(8)p2DxD,yD,0=0.The outer boundary condition is(9)∂p2DxD,yeD,tD∂yD=0.The interface conditions of Regions 2 and 4 are as follows:(10)p1DxD,ywD+1,tD=p3DxD,ywD+1,tD,(11)∂p1DxD,ywD+1,tD∂yD=∂p3DxD,ywD+1,tD∂yD.


#### 2.2.3. Governing Equations in Region 3

In Region 3, the one-dimensional flow in thex-direction is dominant. Take the integral average along the y-direction in Region 3. The flow is simplified as a linear flow perpendicular to the hydraulic fracture plane given by
(12)∂2p3D∂xD2+∫ywDywD+1∂∂yD∂pD∂yDdα=1η3D∂p3D∂tD.Since the normal plane of the centre of the hydraulic fracture plane has distributivity, combined with Equation (5), a partial differential equation can be attained as follows:
(13)∂2p3D∂xD2+∂p1DywD+1,tD∂yD=1η3D∂p3D∂tD.The initial condition is(14)p3DxD,yD,0=0.The outer boundary condition is(15)∂p3DxeD,yD,tD∂xD=0.The interface conditions of Regions 3 and 4 are as follows:(16)p3DxwmD,yD,tD=p4DxwmD,yD,tD,(17)kμ∂p3DxwmD,ywD,tD∂xD=k4μ∂pmDxwmD,ywD,tD∂xD.


#### 2.2.4. Governing Equations in Region 4

In Region 4, the matrix’s permeability is much lower than that of fracture, so the pressure drop rate is slower. There is a pressure difference between the matrix and the fracture, resulting in cross-flow. According to the law of conservation of mass, Darcy law, and the classical Warren-Root model, SRV region’s flow is assumed to be one-dimensional and in thex-direction perpendicular to the hydraulic fracture plane. The associated equations are given by
(18)∂2p4fD∂xD2+1yeDCRD∂p2DywD+1,s∂yD+λp4mD−p4fD=ω∂p4fD∂tD,(19)λp4mD−p4fD=1−ω∂p4mD∂tD.The initial condition is(20)p4mDxD,yD,0=p4fDxD,yD,0=0.The inner boundary condition is(21)∂p4fDxwD,tD∂xD=−π2constantrate.


### 2.3. Model Solution

We derive the solutions for Regions 1, 2, and 3 and the SRV region and then couple the solutions by using the flux and pressure continuity conditions on the interfaces between regions. The detailed derivation of the preceding equations with boundary conditions is provided in AppendixB.
(22)spwDinf~0,s=π2·1C2s·1+kDC1s/C2stanhxmDC2stanhxeD−xwmDC1stanhxmDC2s+kDC1s/C2stanhxeD−xwmDC1s,C1s=sη3D+sη3DtanhyeD−ywD−1sη3D,C2s=us+kDsη2D·tanhyeD−ywD−1sη2D.According to the convolution theorem, the wellbore flow rate under constant pressure is(23)sq~wDs=2C2sπtanhxmDC2s+kDC1s/C2stanhxeD−xwmDC1s1+kDC1s/C2stanhxmDC2stanhxeD−xwmDC1s.This solution can be inverted from the Laplace to the time domain with the numerical algorithm given by Stehfest [30].


#### 2.3.1. Conductivity Influence Function

Wilkinsons [22] placed the fracture in the formation with impermeable boundaries at its endpoint and analytically solved the early flow characteristic function of low conductivity fracture by the Fourier transform method. Based on the asymptotic fitting analysis method, the modified transformation is
(24)sf~1cfD=∑n=1∞2πn2π2cfD+2n2π2+s.Based on the numerical simulation results, the correction functionΔf~cfD is introduced to improve the accuracy of the hybrid solution to yield
(25)Δf~cfD=πΔEπ+ΔE·4s,Δf~cfD≈ΔEs→0,π4ss→∞,(26)ΔE=c11/cfDc2c11/cfD+c2,c1=4π−2π−π3,c2=lnπ2.Finally, a more accurate solution at the wellbore is obtained as(27)s~pwDs=sp~wDinf0,s+sf~1cfD+Δf~cfD.At the same time, the pseudo-steady-state flow asymptotic analysis (s approaches 0) of Equation (25) is done, and the pseudo-steady-state expression of the finite conductivity function in the time domain can be written as follows:
(28)fCfD=∑n=1∞2nnπcfD+2+ΔE.Wang et al. [25] gave the finite conductivity factor in the pseudo-steady-state through progressive analysis and multiple regression methods as follows:
(29)fCfD=0.95−0.56ψ+0.16ψ2−0.028ψ3+0.0028ψ4−0.00011ψ51+0.094ψ+0.093ψ2+0.0084ψ3+0.001ψ4+0.00036ψ5,where Ψ=lnCfD.Equation (29) is a fitting function and has no specific physical meaning but can simplify the calculation. Comparing and analysing Equations (28) and (29), the results show that the two methods provide good fitting results when the value of fracture conductivity is greater than 0.0631 (see Table 1 and Figure 2), indicating that this paper’s finite conductivity factor is reliable.Table 1
Fracture conductivity influence function.
CfDfCfD Equation (29)fCfDRelative error (%)CfDfCfD Equation (29)fCfDRelative error (%)0.017.83045.1652651.59740.316231.848651.835790.70050.012596.880144.9379339.33250.398111.648161.637750.63560.015856.126714.7105730.06300.501191.45681.448180.59520.019955.511794.4832522.94180.630961.275751.268710.55490.025124.99744.2560717.41820.794331.106381.100890.49870.031624.557654.0291513.116920.566090.56731-0.21500.039814.174243.802639.772430.40260.40658-0.97890.050123.833813.576677.189440.311410.31707-1.78510.06313.526373.351515.217350.253190.25993-2.59300.079433.244353.127433.738560.212780.22025-3.39160.12.981952.90482.655970.183120.19109-4.17080.125892.734782.684111.887880.160430.16875-4.93040.158492.499622.465981.364290.142530.15109-5.66550.199532.274272.251191.0252100.128050.13676-6.36880.251192.057432.040720.8188Figure 2
A validation analysis of the finite conductivity function.It deserves noting that when the SRV (Region 4) occupies all the space between fractures and boundaries, the five-region model is simplified to the trilinear model of Brown et al. [9]. When the contributions of Regions 1 and 2 are neglected, the five-region model reduces to the enhanced fracture region model of Stalgorova and Mattar [13]. For this reason, the five-region model is a generalization that covers both the five-linear model and Warren-Root model.


### 2.4. Possible Model Modifications

#### 2.4.1. Gas Flow

This model was derived for an oil system. However, it can also be applied to gas flow. It should be noted that the diffusivity term,η, is not constant but varies with pressure. Therefore, the gas reservoir time should be expressed in terms of pseudotime, and the pressure should be handled through a pseudopressure transformation [14, 15, 31].
(30)mp=2∫prefpξμξzξdξ,ta=∫0tμgicgiμξcξdξ.


#### 2.4.2. Wellbore Storage

According to the recommendations of Brown et al. [9], wellbore storage effects can also be included. To consider wellbore storage, the wellbore pressure given in Equation (22) can be modified to
(31)p~wD,storage=p~wD1+CfDs2p~wD.


### 2.5. Model Validation

To validate the developed model, comparisons between the analytical and numerical simulation results under different working systems were conducted. Isothermal single-phase fluid flow is assumed. The different working systems include the constant flow rate system and compound system. The semianalytical model is inverted with the Stehfest algorithm [30]. For numerical simulation, the 2D finite-difference method is used to simulate flow in the reservoir. The formation was modeled with 150×30 grid cells. A local grid fine area around the hydraulic fracture was defined to ensure the difference in model’s accuracy. The reservoir parameters and well parameters are shown in Table 2. The numerical and analytical solution results are shown in Figure 3. For these parameters, the numerical and analytical solutions in the different working systems are in excellent agreement. Minor divergences are believed because numerical simulator has more considerations. Therefore, the semianalytical model can be applied to predict long-term well performance.Table 2
Model parameters.
ParametersSymbolValueInitial pressure (MPa)pi28.5Reservoir length (m)xe1500Reservoir width (m)ye300Oil viscosity (mPa·s)μo8Oil volume factor (m3/m3)Bo1.2Porosityφ0.2Formation permeability (mD)k3SRV region permeability (mD)k480Fracture half-length (m)yf50SRV region width (m)2x4150Reservoir compressibility (MPa-1)c0.00373Figure 3
(a) Comparing analytical and numerical solutions in the constant rate system. (b) Comparing analytical and numerical solutions in the compound working system.
(a)(b)


## 2.1. Physical Model

The fracture network formation near the wellbore is mainly related to reservoir pressure, mechanical rock properties, and natural fractures. The weak structural surfaces between rocks cause volume fractures such as natural fractures, joints, and bedding or weak points in the formation. Furthermore, the rock brittleness index is the physical property basis for forming a fracture network. Modern hydraulic fracturing technology can form a particular enhanced region with induced fractures near hydraulic fractures by controlling the injection and reverse discharge flow [29]. In this way, the near-wellbore zone is a fracture network formed by the interweaving of primary fractures and multistage secondary fractures. The area of the secondary fracture is also known as the stimulated reservoir volume (SRV). The SRV region is simulated through a Warren-Root dual-porosity medium model.Figure1 shows the schematic diagram of the five-linear flow physical model. In the centre of a closed rectangular drainage area, there is a production well with hydraulic fractures that vertically penetrate through the reservoir. Only a quarter of the flow area will be considered because of the symmetry of the formation. The flow areas are divided into Regions 1, 2, 3, 4, and 5. Flow in the model is regarded as a combination of five linear flows within contiguous regions. Each region comprises one-dimensional linear flow, and the boundaries between regions are continuous in flux and pressure. To keep the solution general, each region has various parameters (permeability, porosity, and total compressibility). Region 4 represents the SRV region near the hydraulic fracture, while Region 5 represents the finite conductivity hydraulic fractures.Figure 1
(a) Schematic of the five-linear flow model representing five continuous flow regions for a vertical well with an SRV region. (b) The matrix of SRV regions. (c) Classic Warren-Root dual-porosity medium model.
(a)(b)(c)The basic assumptions of the model are as follows:(i)
The hydraulic fracturing vertical well is located in the centre of a rectangular low permeability oil reservoir with closed boundaries(ii)
The thickness of the reservoir ish, and the uniform initial pressure is pi(iii)
The hydraulic fractures thoroughly penetrate the formation(iv)
The effect of gravity and capillary forces on fluid flow in each region is negligible(v)
Isothermal single-phase fluid flow is assumed


## 2.2. Mathematical Model

The pressure distribution in each region is controlled by the flow equation given by(1)∇2p=ϕμck∂p∂t.To simplify the solution, this equation for each region was rewritten in dimensionless terms. The definitions of dimensionless parameters are given in AppendixA, and the derivations of equations and solutions are given in Appendix B.


### 2.2.1. Governing Equations in Region 1

Region 1 produces one-dimensional flow in they-direction. This visualization of the flow behaviour is reasonable for fluid flow before the pseudo-steady-state flow period. Considering the symmetry, the flow governing differential equations and boundary conditions are as follows:
(2)∂2p1D∂yD2=1η1D∂p1D∂tD,xwD+xmD≤xD≤xeD,ywD+1≤yD≤yeD.The initial condition is(3)p1DxD,yD,0=0.The outer boundary condition is(4)∂p1DxD,yeD,tD∂yD=0.The interface conditions of Regions 1 and 3 are(5)p1DxD,ywD+1,tD=p3DxD,ywD+1,tD,(6)∂p1DxD,ywD+1,tD∂yD=∂p3DxD,ywD+1,tD∂yD.


### 2.2.2. Governing Equations in Region 2

In Region 2, the flow is assumed to be one-dimensional and in they-direction. Considering the symmetry, the dimensionless governing equations are as follows xwmD=xwD+xmD:
(7)∂2p2D∂yD2=1η2D∂p2D∂tD,xwD≤xD≤xwmD,ywD+1≤yD≤yeD.The initial condition is(8)p2DxD,yD,0=0.The outer boundary condition is(9)∂p2DxD,yeD,tD∂yD=0.The interface conditions of Regions 2 and 4 are as follows:(10)p1DxD,ywD+1,tD=p3DxD,ywD+1,tD,(11)∂p1DxD,ywD+1,tD∂yD=∂p3DxD,ywD+1,tD∂yD.


### 2.2.3. Governing Equations in Region 3

In Region 3, the one-dimensional flow in thex-direction is dominant. Take the integral average along the y-direction in Region 3. The flow is simplified as a linear flow perpendicular to the hydraulic fracture plane given by
(12)∂2p3D∂xD2+∫ywDywD+1∂∂yD∂pD∂yDdα=1η3D∂p3D∂tD.Since the normal plane of the centre of the hydraulic fracture plane has distributivity, combined with Equation (5), a partial differential equation can be attained as follows:
(13)∂2p3D∂xD2+∂p1DywD+1,tD∂yD=1η3D∂p3D∂tD.The initial condition is(14)p3DxD,yD,0=0.The outer boundary condition is(15)∂p3DxeD,yD,tD∂xD=0.The interface conditions of Regions 3 and 4 are as follows:(16)p3DxwmD,yD,tD=p4DxwmD,yD,tD,(17)kμ∂p3DxwmD,ywD,tD∂xD=k4μ∂pmDxwmD,ywD,tD∂xD.


### 2.2.4. Governing Equations in Region 4

In Region 4, the matrix’s permeability is much lower than that of fracture, so the pressure drop rate is slower. There is a pressure difference between the matrix and the fracture, resulting in cross-flow. According to the law of conservation of mass, Darcy law, and the classical Warren-Root model, SRV region’s flow is assumed to be one-dimensional and in thex-direction perpendicular to the hydraulic fracture plane. The associated equations are given by
(18)∂2p4fD∂xD2+1yeDCRD∂p2DywD+1,s∂yD+λp4mD−p4fD=ω∂p4fD∂tD,(19)λp4mD−p4fD=1−ω∂p4mD∂tD.The initial condition is(20)p4mDxD,yD,0=p4fDxD,yD,0=0.The inner boundary condition is(21)∂p4fDxwD,tD∂xD=−π2constantrate.


## 2.2.1. Governing Equations in Region 1

Region 1 produces one-dimensional flow in they-direction. This visualization of the flow behaviour is reasonable for fluid flow before the pseudo-steady-state flow period. Considering the symmetry, the flow governing differential equations and boundary conditions are as follows:
(2)∂2p1D∂yD2=1η1D∂p1D∂tD,xwD+xmD≤xD≤xeD,ywD+1≤yD≤yeD.The initial condition is(3)p1DxD,yD,0=0.The outer boundary condition is(4)∂p1DxD,yeD,tD∂yD=0.The interface conditions of Regions 1 and 3 are(5)p1DxD,ywD+1,tD=p3DxD,ywD+1,tD,(6)∂p1DxD,ywD+1,tD∂yD=∂p3DxD,ywD+1,tD∂yD.


## 2.2.2. Governing Equations in Region 2

In Region 2, the flow is assumed to be one-dimensional and in they-direction. Considering the symmetry, the dimensionless governing equations are as follows xwmD=xwD+xmD:
(7)∂2p2D∂yD2=1η2D∂p2D∂tD,xwD≤xD≤xwmD,ywD+1≤yD≤yeD.The initial condition is(8)p2DxD,yD,0=0.The outer boundary condition is(9)∂p2DxD,yeD,tD∂yD=0.The interface conditions of Regions 2 and 4 are as follows:(10)p1DxD,ywD+1,tD=p3DxD,ywD+1,tD,(11)∂p1DxD,ywD+1,tD∂yD=∂p3DxD,ywD+1,tD∂yD.


## 2.2.3. Governing Equations in Region 3

In Region 3, the one-dimensional flow in thex-direction is dominant. Take the integral average along the y-direction in Region 3. The flow is simplified as a linear flow perpendicular to the hydraulic fracture plane given by
(12)∂2p3D∂xD2+∫ywDywD+1∂∂yD∂pD∂yDdα=1η3D∂p3D∂tD.Since the normal plane of the centre of the hydraulic fracture plane has distributivity, combined with Equation (5), a partial differential equation can be attained as follows:
(13)∂2p3D∂xD2+∂p1DywD+1,tD∂yD=1η3D∂p3D∂tD.The initial condition is(14)p3DxD,yD,0=0.The outer boundary condition is(15)∂p3DxeD,yD,tD∂xD=0.The interface conditions of Regions 3 and 4 are as follows:(16)p3DxwmD,yD,tD=p4DxwmD,yD,tD,(17)kμ∂p3DxwmD,ywD,tD∂xD=k4μ∂pmDxwmD,ywD,tD∂xD.


## 2.2.4. Governing Equations in Region 4

In Region 4, the matrix’s permeability is much lower than that of fracture, so the pressure drop rate is slower. There is a pressure difference between the matrix and the fracture, resulting in cross-flow. According to the law of conservation of mass, Darcy law, and the classical Warren-Root model, SRV region’s flow is assumed to be one-dimensional and in thex-direction perpendicular to the hydraulic fracture plane. The associated equations are given by
(18)∂2p4fD∂xD2+1yeDCRD∂p2DywD+1,s∂yD+λp4mD−p4fD=ω∂p4fD∂tD,(19)λp4mD−p4fD=1−ω∂p4mD∂tD.The initial condition is(20)p4mDxD,yD,0=p4fDxD,yD,0=0.The inner boundary condition is(21)∂p4fDxwD,tD∂xD=−π2constantrate.


## 2.3. Model Solution

We derive the solutions for Regions 1, 2, and 3 and the SRV region and then couple the solutions by using the flux and pressure continuity conditions on the interfaces between regions. The detailed derivation of the preceding equations with boundary conditions is provided in AppendixB.
(22)spwDinf~0,s=π2·1C2s·1+kDC1s/C2stanhxmDC2stanhxeD−xwmDC1stanhxmDC2s+kDC1s/C2stanhxeD−xwmDC1s,C1s=sη3D+sη3DtanhyeD−ywD−1sη3D,C2s=us+kDsη2D·tanhyeD−ywD−1sη2D.According to the convolution theorem, the wellbore flow rate under constant pressure is(23)sq~wDs=2C2sπtanhxmDC2s+kDC1s/C2stanhxeD−xwmDC1s1+kDC1s/C2stanhxmDC2stanhxeD−xwmDC1s.This solution can be inverted from the Laplace to the time domain with the numerical algorithm given by Stehfest [30].


### 2.3.1. Conductivity Influence Function

Wilkinsons [22] placed the fracture in the formation with impermeable boundaries at its endpoint and analytically solved the early flow characteristic function of low conductivity fracture by the Fourier transform method. Based on the asymptotic fitting analysis method, the modified transformation is
(24)sf~1cfD=∑n=1∞2πn2π2cfD+2n2π2+s.Based on the numerical simulation results, the correction functionΔf~cfD is introduced to improve the accuracy of the hybrid solution to yield
(25)Δf~cfD=πΔEπ+ΔE·4s,Δf~cfD≈ΔEs→0,π4ss→∞,(26)ΔE=c11/cfDc2c11/cfD+c2,c1=4π−2π−π3,c2=lnπ2.Finally, a more accurate solution at the wellbore is obtained as(27)s~pwDs=sp~wDinf0,s+sf~1cfD+Δf~cfD.At the same time, the pseudo-steady-state flow asymptotic analysis (s approaches 0) of Equation (25) is done, and the pseudo-steady-state expression of the finite conductivity function in the time domain can be written as follows:
(28)fCfD=∑n=1∞2nnπcfD+2+ΔE.Wang et al. [25] gave the finite conductivity factor in the pseudo-steady-state through progressive analysis and multiple regression methods as follows:
(29)fCfD=0.95−0.56ψ+0.16ψ2−0.028ψ3+0.0028ψ4−0.00011ψ51+0.094ψ+0.093ψ2+0.0084ψ3+0.001ψ4+0.00036ψ5,where Ψ=lnCfD.Equation (29) is a fitting function and has no specific physical meaning but can simplify the calculation. Comparing and analysing Equations (28) and (29), the results show that the two methods provide good fitting results when the value of fracture conductivity is greater than 0.0631 (see Table 1 and Figure 2), indicating that this paper’s finite conductivity factor is reliable.Table 1
Fracture conductivity influence function.
CfDfCfD Equation (29)fCfDRelative error (%)CfDfCfD Equation (29)fCfDRelative error (%)0.017.83045.1652651.59740.316231.848651.835790.70050.012596.880144.9379339.33250.398111.648161.637750.63560.015856.126714.7105730.06300.501191.45681.448180.59520.019955.511794.4832522.94180.630961.275751.268710.55490.025124.99744.2560717.41820.794331.106381.100890.49870.031624.557654.0291513.116920.566090.56731-0.21500.039814.174243.802639.772430.40260.40658-0.97890.050123.833813.576677.189440.311410.31707-1.78510.06313.526373.351515.217350.253190.25993-2.59300.079433.244353.127433.738560.212780.22025-3.39160.12.981952.90482.655970.183120.19109-4.17080.125892.734782.684111.887880.160430.16875-4.93040.158492.499622.465981.364290.142530.15109-5.66550.199532.274272.251191.0252100.128050.13676-6.36880.251192.057432.040720.8188Figure 2
A validation analysis of the finite conductivity function.It deserves noting that when the SRV (Region 4) occupies all the space between fractures and boundaries, the five-region model is simplified to the trilinear model of Brown et al. [9]. When the contributions of Regions 1 and 2 are neglected, the five-region model reduces to the enhanced fracture region model of Stalgorova and Mattar [13]. For this reason, the five-region model is a generalization that covers both the five-linear model and Warren-Root model.


## 2.3.1. Conductivity Influence Function

Wilkinsons [22] placed the fracture in the formation with impermeable boundaries at its endpoint and analytically solved the early flow characteristic function of low conductivity fracture by the Fourier transform method. Based on the asymptotic fitting analysis method, the modified transformation is
(24)sf~1cfD=∑n=1∞2πn2π2cfD+2n2π2+s.Based on the numerical simulation results, the correction functionΔf~cfD is introduced to improve the accuracy of the hybrid solution to yield
(25)Δf~cfD=πΔEπ+ΔE·4s,Δf~cfD≈ΔEs→0,π4ss→∞,(26)ΔE=c11/cfDc2c11/cfD+c2,c1=4π−2π−π3,c2=lnπ2.Finally, a more accurate solution at the wellbore is obtained as(27)s~pwDs=sp~wDinf0,s+sf~1cfD+Δf~cfD.At the same time, the pseudo-steady-state flow asymptotic analysis (s approaches 0) of Equation (25) is done, and the pseudo-steady-state expression of the finite conductivity function in the time domain can be written as follows:
(28)fCfD=∑n=1∞2nnπcfD+2+ΔE.Wang et al. [25] gave the finite conductivity factor in the pseudo-steady-state through progressive analysis and multiple regression methods as follows:
(29)fCfD=0.95−0.56ψ+0.16ψ2−0.028ψ3+0.0028ψ4−0.00011ψ51+0.094ψ+0.093ψ2+0.0084ψ3+0.001ψ4+0.00036ψ5,where Ψ=lnCfD.Equation (29) is a fitting function and has no specific physical meaning but can simplify the calculation. Comparing and analysing Equations (28) and (29), the results show that the two methods provide good fitting results when the value of fracture conductivity is greater than 0.0631 (see Table 1 and Figure 2), indicating that this paper’s finite conductivity factor is reliable.Table 1
Fracture conductivity influence function.
CfDfCfD Equation (29)fCfDRelative error (%)CfDfCfD Equation (29)fCfDRelative error (%)0.017.83045.1652651.59740.316231.848651.835790.70050.012596.880144.9379339.33250.398111.648161.637750.63560.015856.126714.7105730.06300.501191.45681.448180.59520.019955.511794.4832522.94180.630961.275751.268710.55490.025124.99744.2560717.41820.794331.106381.100890.49870.031624.557654.0291513.116920.566090.56731-0.21500.039814.174243.802639.772430.40260.40658-0.97890.050123.833813.576677.189440.311410.31707-1.78510.06313.526373.351515.217350.253190.25993-2.59300.079433.244353.127433.738560.212780.22025-3.39160.12.981952.90482.655970.183120.19109-4.17080.125892.734782.684111.887880.160430.16875-4.93040.158492.499622.465981.364290.142530.15109-5.66550.199532.274272.251191.0252100.128050.13676-6.36880.251192.057432.040720.8188Figure 2
A validation analysis of the finite conductivity function.It deserves noting that when the SRV (Region 4) occupies all the space between fractures and boundaries, the five-region model is simplified to the trilinear model of Brown et al. [9]. When the contributions of Regions 1 and 2 are neglected, the five-region model reduces to the enhanced fracture region model of Stalgorova and Mattar [13]. For this reason, the five-region model is a generalization that covers both the five-linear model and Warren-Root model.


## 2.4. Possible Model Modifications

### 2.4.1. Gas Flow

This model was derived for an oil system. However, it can also be applied to gas flow. It should be noted that the diffusivity term,η, is not constant but varies with pressure. Therefore, the gas reservoir time should be expressed in terms of pseudotime, and the pressure should be handled through a pseudopressure transformation [14, 15, 31].
(30)mp=2∫prefpξμξzξdξ,ta=∫0tμgicgiμξcξdξ.


### 2.4.2. Wellbore Storage

According to the recommendations of Brown et al. [9], wellbore storage effects can also be included. To consider wellbore storage, the wellbore pressure given in Equation (22) can be modified to
(31)p~wD,storage=p~wD1+CfDs2p~wD.


## 2.4.1. Gas Flow

This model was derived for an oil system. However, it can also be applied to gas flow. It should be noted that the diffusivity term,η, is not constant but varies with pressure. Therefore, the gas reservoir time should be expressed in terms of pseudotime, and the pressure should be handled through a pseudopressure transformation [14, 15, 31].
(30)mp=2∫prefpξμξzξdξ,ta=∫0tμgicgiμξcξdξ.


## 2.4.2. Wellbore Storage

According to the recommendations of Brown et al. [9], wellbore storage effects can also be included. To consider wellbore storage, the wellbore pressure given in Equation (22) can be modified to
(31)p~wD,storage=p~wD1+CfDs2p~wD.


## 2.5. Model Validation

To validate the developed model, comparisons between the analytical and numerical simulation results under different working systems were conducted. Isothermal single-phase fluid flow is assumed. The different working systems include the constant flow rate system and compound system. The semianalytical model is inverted with the Stehfest algorithm [30]. For numerical simulation, the 2D finite-difference method is used to simulate flow in the reservoir. The formation was modeled with 150×30 grid cells. A local grid fine area around the hydraulic fracture was defined to ensure the difference in model’s accuracy. The reservoir parameters and well parameters are shown in Table 2. The numerical and analytical solution results are shown in Figure 3. For these parameters, the numerical and analytical solutions in the different working systems are in excellent agreement. Minor divergences are believed because numerical simulator has more considerations. Therefore, the semianalytical model can be applied to predict long-term well performance.Table 2
Model parameters.
ParametersSymbolValueInitial pressure (MPa)pi28.5Reservoir length (m)xe1500Reservoir width (m)ye300Oil viscosity (mPa·s)μo8Oil volume factor (m3/m3)Bo1.2Porosityφ0.2Formation permeability (mD)k3SRV region permeability (mD)k480Fracture half-length (m)yf50SRV region width (m)2x4150Reservoir compressibility (MPa-1)c0.00373Figure 3
(a) Comparing analytical and numerical solutions in the constant rate system. (b) Comparing analytical and numerical solutions in the compound working system.
(a)(b)


## 3. Results and Discussion

Sensitivity studies illustrate the essential parameters affecting oil production performance. The parameters discussed include the SRV size and shape, fracture conductivity, mobility ratio, initial flowing rate, and limiting wellbore pressure under different working systems. Constant rate, constant pressure, and compound working systems are taken into account. Parameters for Regions 1, 2, and 3 can be considered identical for practical purposes. The most important influencing factors in the figures are displayed in a dimensional form to make the results intuitive.


### 3.1. Constant Production Working System

#### 3.1.1. Shape and Size of Stimulated Reservoir Volume

Figures4 and 5 show the effect of the SRV region’s different sizes and shapes on the wellbore pressure. In Figure 4, the SRV length remains fixed at 150 m, and the SRV width is varied (100 m, 150 m, 180 m, 200 m, and 240 m). As the SRV area increases, the rate of pressure at the wellbore decline slows down by 54%. This is mainly because, in the early stages, flow only occurs in the hydraulic fractures and the stimulated region near the well, which is not affected by the size of SRV. In the exploitation process, the production wave spreads to the SRV region edges, where reservoirs with larger effective SRVs have a lower pressure decline rate. Therefore, an appropriate SRV size should be ensured in the fracturing design.Figure 4
The effect of SRV size on the curves (constant rate).Figure 5
(a) The effect of SRV shape on the flow rate and cumulative yield. (b) Schematic diagram of five different SRV shapes (constant rate).
(a)(b)Figure5(a) shows that the x-direction rectangular shape of the SRV results in a slower wellbore pressure decline rate. Figure 5(b) is the schematic diagram of five different SRV shapes, and the SRV size is fixed at 25600 m3. Under constant production, increasing the ratio of length to width can effectively slow down formation energy loss. In the production process, the production wave spreads to the edges of the SRV region, and unstimulated regions participate in the flow. Region 2 is the first and dominant region. When the shape of the enhanced region is similar to the formation, Region 2 has a more extensive scope and sufficient supply, resulting in better performance.


#### 3.1.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=0.01,0.1,1,10,125) under the condition of constant flow rate is shown in Figure 6. At the initial stage of production, the bottom hole flow pressure changes slowly, and the reservoir of energy loss is small. At the middle stage, the wellbore flow pressure decreases rapidly. When the fracture conductivity increases, the rate of wellbore pressure decline decreases by about 32%. For fracture conductivities of 10 and above, pressures are the same. That is, for CfD≥10, there is a pronounced reduction in formation energy loss.Figure 6
The effect of fracture conductivity on the wellbore pressure (constant rate).The influence of mobility ratio between the enhanced area and matrix area on wellbore pressure is shown in Figure7. Here, the SRV permeability of k4 is varied while all other properties are kept constant. It can be seen that, as the value of mobility ratio increases, the rate of wellbore pressure decline slows down by 83.9%. Since the matrix fluid supply capacity is sufficient and the enhanced SRV permeability helps with fluid flow, the wellbore pressure drops slowly as a result. Thus, it is essential to improve the SRV permeability in the fracturing design to extend the constant rate period.Figure 7
Effect of mobility ratio on the wellbore pressure (constant rate).


### 3.2. Constant Pressure Working System

#### 3.2.1. Shape and Size of Stimulated Reservoir Volume

Under a constant wellbore pressure, the effects of SRV size and shape on the production are shown in Figures8 and 9. In Figure 8, the width of the SRV remains fixed and the length is varied (ym=20,40,100,180,240m). When the fracture length is relatively small (ym<180m), a larger SRV size leads to higher cumulative oil production. When the fracture length is increased, increasingly induced fractures in the SRV region become effectively connected with artificial hydraulic fractures, forming a larger supply area and improving the cumulative production. This, however, has little effect on the production when the width is greater than 180 m as a result of the reservoirs controlled by a single are very limited. Considering the cost and difficulty of fracturing factors, an appropriate fracturing length should be considered.Figure 8
The effect of SRV size on cumulative production (constant pressure).Figure 9
The effect of SRV shape on production (constant pressure).As shown in Figure9, the enhanced area’s volume is fixed at 256000 m3, while the length-width ratio of the stimulated area is varied (xm/ym=0.38, 0.64, 1.00, 2.56, and 7.11). The schematic diagram of SRV shapes is displayed in Figure 5(b). Results show that an SRV rectangular shape along the x-direction enables higher production during the middle stage. In the production process, the production wave spreads to the edges of the SRV region. The unstimulated regions participate in the flow in which Region 2 is the first and dominant region. When the shape of the enhanced region is similar to the formation, Region 2 has a more extensive scope and sufficient supply, leading to improved performance.


#### 3.2.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=10,50,100,300,500) under the condition of constant pressure is shown in Figure 10. As can be seen, an increase in fracture conductivity can increase productivity by up to 814%. The effect of hydraulic fractures is to create a high permeability channel near the wellbore formation, which is convenient for oil to flow from the far formation to the bottom hole. For CfD>300, the cumulative production curves coincide, which means that hydraulic fractures can be regarded as uniform flow fractures.Figure 10
The effect of fracture conductivity on cumulative production (constant pressure).Figure11 represents the impact of the mobility ratio between the SRV region and other regions on the flow rate curves. In these scenarios, the formation permeability, k, is kept constant while the value of k4 is varied to yield different mobility ratios. The larger the mobility ratio is, the more pronounced the stimulation effect. Due to the high permeability and strong fluid supply capacity of the formation in the enhanced area, the productivity is naturally higher.Figure 11
The effect of mobility ratio on cumulative production (constant pressure).


### 3.3. Compound Working System

Under the compound working system, an oil well’s production performance can be summarized in a chronological manner. In the first stage, production is begun at a constant rate, and well bottom pressure declines rapidly. When wellbore pressure declines to a limiting value, the producing rate begins to decline as the well produces against a constant pressure [32]. The compound working system developed in this work enables the prediction of long-term wells’ performance with a high level of accuracy. It is very consistent with the actual production situation of oil wells.


#### 3.3.1. Shape and Size of Stimulated Reservoir Volume

The effects of SRV size and shape on the production rate are shown in Figures12 and 13, which demonstrate their significant influences on both the flow rate and cumulative production. Comparing the three models shown in Figure 12, the larger the SRV size is, the longer the constant flow period and the higher cumulative production. The larger size effect is to extend the constant period by about 166% and increase the cumulative production by about 30%. However, the rate difference becomes minimal after 500 days. This is because the SRV region’s oil yield has been primarily produced, leaving other regions to be the dominant factors for production.Figure 12
The effect of SRV size on the flow rate and cumulative yield (compound working system).Figure 13
The effect of SRV shape on flow rate and cumulative yield (compound working system).Figure13 shows that when the SRV region’s length-width ratio increases, the flow rate is higher. In other words, an SRV shape similar to the formation can help increase productivity. The effect of the x-direction rectangular shape of the SRV increases the cumulative production by about 56%. In the initial period, flow only happens in the hydraulic fractures and the stimulated region near the well. Therefore, producing rate is not affected by the shape and size of the SRV. In the exploitation process, the production wave spreads to the SRV region’s edges, and reservoirs with a similar shape to the SRV have a higher harvest. Therefore, in the fracturing design, the stimulated reservoir volume size and shape should be seriously discussed.


#### 3.3.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=0.01,20,300) under the compound working system is shown in Figure 14. For a fracture conductivity of 0.01 (and other properties and parameters remaining constant), the production period at a constant rate lasts only 270 d. For fracture conductivities of 20 and 300, the period increases to 410 d, i.e., an increase of 51.9%. Moreover, the flow rate decreases rapidly in the constant pressure stage. These results demonstrate that the fracture conductivity helps extend the stable production period and achieve higher production within a short time.Figure 14
The effect of fracture conductivity on well performance (compound working system).The influence of the mobility ratio between the SRV region and matrix area on the dimensionless pressure is shown in Figure15. Here, the SRV permeability of k4 is varied while keeping the well and oil characteristics constant. It can be seen that with the value of mobility ratio increasing, the producing period at a constant rate extends by 150%, and the rate of pressure decline decreases. This is because improving the permeability of the SRV aids fluid flow.Figure 15
The effect of mobility ratio on well performance (compound working system).


#### 3.3.3. Initial Flow Rate and Limiting Wellbore Pressure

Plots of producing rate and wellbore pressure vs. time for different initial flow rates (qoi=15,25,35m3/day) are shown in Figure 16. The producing rate decreases from 35 m3/day to 15 m3/day while the stable production period increases from 220 d to 1140 d. The effect of the lower initial flow rate is to extend the producing period at a constant rate of about 418%. Pressure decline at the wellbore steadily decreases during this period. Within a short period of time, a higher cumulative yield can be obtained. However, the ultimate recovery is not significantly different.Figure 16
Effect of initial flow rate on the performance curve (compound working system).Plots of producing rate and wellbore pressure vs. time for limiting wellbore pressure (pwfsc=6,8,10MPa) are shown in Figure 17. The conversion pressure decreases from 10 MPa to 6 MPa, while the stable production stage increases by 46% from 280 d to 410 d. The pressure drop rate is consistent, but the final output increases by 23% from 2.6856×104m3 to 3.2668×104m3. This shows that the conversion pressure can effectively extend the first period and significantly increase the well’s ultimate production at the same initial flow rate.Figure 17
Effect of limiting wellbore pressure on well performance curve (compound working system).


## 3.1. Constant Production Working System

### 3.1.1. Shape and Size of Stimulated Reservoir Volume

Figures4 and 5 show the effect of the SRV region’s different sizes and shapes on the wellbore pressure. In Figure 4, the SRV length remains fixed at 150 m, and the SRV width is varied (100 m, 150 m, 180 m, 200 m, and 240 m). As the SRV area increases, the rate of pressure at the wellbore decline slows down by 54%. This is mainly because, in the early stages, flow only occurs in the hydraulic fractures and the stimulated region near the well, which is not affected by the size of SRV. In the exploitation process, the production wave spreads to the SRV region edges, where reservoirs with larger effective SRVs have a lower pressure decline rate. Therefore, an appropriate SRV size should be ensured in the fracturing design.Figure 4
The effect of SRV size on the curves (constant rate).Figure 5
(a) The effect of SRV shape on the flow rate and cumulative yield. (b) Schematic diagram of five different SRV shapes (constant rate).
(a)(b)Figure5(a) shows that the x-direction rectangular shape of the SRV results in a slower wellbore pressure decline rate. Figure 5(b) is the schematic diagram of five different SRV shapes, and the SRV size is fixed at 25600 m3. Under constant production, increasing the ratio of length to width can effectively slow down formation energy loss. In the production process, the production wave spreads to the edges of the SRV region, and unstimulated regions participate in the flow. Region 2 is the first and dominant region. When the shape of the enhanced region is similar to the formation, Region 2 has a more extensive scope and sufficient supply, resulting in better performance.


### 3.1.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=0.01,0.1,1,10,125) under the condition of constant flow rate is shown in Figure 6. At the initial stage of production, the bottom hole flow pressure changes slowly, and the reservoir of energy loss is small. At the middle stage, the wellbore flow pressure decreases rapidly. When the fracture conductivity increases, the rate of wellbore pressure decline decreases by about 32%. For fracture conductivities of 10 and above, pressures are the same. That is, for CfD≥10, there is a pronounced reduction in formation energy loss.Figure 6
The effect of fracture conductivity on the wellbore pressure (constant rate).The influence of mobility ratio between the enhanced area and matrix area on wellbore pressure is shown in Figure7. Here, the SRV permeability of k4 is varied while all other properties are kept constant. It can be seen that, as the value of mobility ratio increases, the rate of wellbore pressure decline slows down by 83.9%. Since the matrix fluid supply capacity is sufficient and the enhanced SRV permeability helps with fluid flow, the wellbore pressure drops slowly as a result. Thus, it is essential to improve the SRV permeability in the fracturing design to extend the constant rate period.Figure 7
Effect of mobility ratio on the wellbore pressure (constant rate).


## 3.1.1. Shape and Size of Stimulated Reservoir Volume

Figures4 and 5 show the effect of the SRV region’s different sizes and shapes on the wellbore pressure. In Figure 4, the SRV length remains fixed at 150 m, and the SRV width is varied (100 m, 150 m, 180 m, 200 m, and 240 m). As the SRV area increases, the rate of pressure at the wellbore decline slows down by 54%. This is mainly because, in the early stages, flow only occurs in the hydraulic fractures and the stimulated region near the well, which is not affected by the size of SRV. In the exploitation process, the production wave spreads to the SRV region edges, where reservoirs with larger effective SRVs have a lower pressure decline rate. Therefore, an appropriate SRV size should be ensured in the fracturing design.Figure 4
The effect of SRV size on the curves (constant rate).Figure 5
(a) The effect of SRV shape on the flow rate and cumulative yield. (b) Schematic diagram of five different SRV shapes (constant rate).
(a)(b)Figure5(a) shows that the x-direction rectangular shape of the SRV results in a slower wellbore pressure decline rate. Figure 5(b) is the schematic diagram of five different SRV shapes, and the SRV size is fixed at 25600 m3. Under constant production, increasing the ratio of length to width can effectively slow down formation energy loss. In the production process, the production wave spreads to the edges of the SRV region, and unstimulated regions participate in the flow. Region 2 is the first and dominant region. When the shape of the enhanced region is similar to the formation, Region 2 has a more extensive scope and sufficient supply, resulting in better performance.


## 3.1.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=0.01,0.1,1,10,125) under the condition of constant flow rate is shown in Figure 6. At the initial stage of production, the bottom hole flow pressure changes slowly, and the reservoir of energy loss is small. At the middle stage, the wellbore flow pressure decreases rapidly. When the fracture conductivity increases, the rate of wellbore pressure decline decreases by about 32%. For fracture conductivities of 10 and above, pressures are the same. That is, for CfD≥10, there is a pronounced reduction in formation energy loss.Figure 6
The effect of fracture conductivity on the wellbore pressure (constant rate).The influence of mobility ratio between the enhanced area and matrix area on wellbore pressure is shown in Figure7. Here, the SRV permeability of k4 is varied while all other properties are kept constant. It can be seen that, as the value of mobility ratio increases, the rate of wellbore pressure decline slows down by 83.9%. Since the matrix fluid supply capacity is sufficient and the enhanced SRV permeability helps with fluid flow, the wellbore pressure drops slowly as a result. Thus, it is essential to improve the SRV permeability in the fracturing design to extend the constant rate period.Figure 7
Effect of mobility ratio on the wellbore pressure (constant rate).


## 3.2. Constant Pressure Working System

### 3.2.1. Shape and Size of Stimulated Reservoir Volume

Under a constant wellbore pressure, the effects of SRV size and shape on the production are shown in Figures8 and 9. In Figure 8, the width of the SRV remains fixed and the length is varied (ym=20,40,100,180,240m). When the fracture length is relatively small (ym<180m), a larger SRV size leads to higher cumulative oil production. When the fracture length is increased, increasingly induced fractures in the SRV region become effectively connected with artificial hydraulic fractures, forming a larger supply area and improving the cumulative production. This, however, has little effect on the production when the width is greater than 180 m as a result of the reservoirs controlled by a single are very limited. Considering the cost and difficulty of fracturing factors, an appropriate fracturing length should be considered.Figure 8
The effect of SRV size on cumulative production (constant pressure).Figure 9
The effect of SRV shape on production (constant pressure).As shown in Figure9, the enhanced area’s volume is fixed at 256000 m3, while the length-width ratio of the stimulated area is varied (xm/ym=0.38, 0.64, 1.00, 2.56, and 7.11). The schematic diagram of SRV shapes is displayed in Figure 5(b). Results show that an SRV rectangular shape along the x-direction enables higher production during the middle stage. In the production process, the production wave spreads to the edges of the SRV region. The unstimulated regions participate in the flow in which Region 2 is the first and dominant region. When the shape of the enhanced region is similar to the formation, Region 2 has a more extensive scope and sufficient supply, leading to improved performance.


### 3.2.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=10,50,100,300,500) under the condition of constant pressure is shown in Figure 10. As can be seen, an increase in fracture conductivity can increase productivity by up to 814%. The effect of hydraulic fractures is to create a high permeability channel near the wellbore formation, which is convenient for oil to flow from the far formation to the bottom hole. For CfD>300, the cumulative production curves coincide, which means that hydraulic fractures can be regarded as uniform flow fractures.Figure 10
The effect of fracture conductivity on cumulative production (constant pressure).Figure11 represents the impact of the mobility ratio between the SRV region and other regions on the flow rate curves. In these scenarios, the formation permeability, k, is kept constant while the value of k4 is varied to yield different mobility ratios. The larger the mobility ratio is, the more pronounced the stimulation effect. Due to the high permeability and strong fluid supply capacity of the formation in the enhanced area, the productivity is naturally higher.Figure 11
The effect of mobility ratio on cumulative production (constant pressure).


## 3.2.1. Shape and Size of Stimulated Reservoir Volume

Under a constant wellbore pressure, the effects of SRV size and shape on the production are shown in Figures8 and 9. In Figure 8, the width of the SRV remains fixed and the length is varied (ym=20,40,100,180,240m). When the fracture length is relatively small (ym<180m), a larger SRV size leads to higher cumulative oil production. When the fracture length is increased, increasingly induced fractures in the SRV region become effectively connected with artificial hydraulic fractures, forming a larger supply area and improving the cumulative production. This, however, has little effect on the production when the width is greater than 180 m as a result of the reservoirs controlled by a single are very limited. Considering the cost and difficulty of fracturing factors, an appropriate fracturing length should be considered.Figure 8
The effect of SRV size on cumulative production (constant pressure).Figure 9
The effect of SRV shape on production (constant pressure).As shown in Figure9, the enhanced area’s volume is fixed at 256000 m3, while the length-width ratio of the stimulated area is varied (xm/ym=0.38, 0.64, 1.00, 2.56, and 7.11). The schematic diagram of SRV shapes is displayed in Figure 5(b). Results show that an SRV rectangular shape along the x-direction enables higher production during the middle stage. In the production process, the production wave spreads to the edges of the SRV region. The unstimulated regions participate in the flow in which Region 2 is the first and dominant region. When the shape of the enhanced region is similar to the formation, Region 2 has a more extensive scope and sufficient supply, leading to improved performance.


## 3.2.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=10,50,100,300,500) under the condition of constant pressure is shown in Figure 10. As can be seen, an increase in fracture conductivity can increase productivity by up to 814%. The effect of hydraulic fractures is to create a high permeability channel near the wellbore formation, which is convenient for oil to flow from the far formation to the bottom hole. For CfD>300, the cumulative production curves coincide, which means that hydraulic fractures can be regarded as uniform flow fractures.Figure 10
The effect of fracture conductivity on cumulative production (constant pressure).Figure11 represents the impact of the mobility ratio between the SRV region and other regions on the flow rate curves. In these scenarios, the formation permeability, k, is kept constant while the value of k4 is varied to yield different mobility ratios. The larger the mobility ratio is, the more pronounced the stimulation effect. Due to the high permeability and strong fluid supply capacity of the formation in the enhanced area, the productivity is naturally higher.Figure 11
The effect of mobility ratio on cumulative production (constant pressure).


## 3.3. Compound Working System

Under the compound working system, an oil well’s production performance can be summarized in a chronological manner. In the first stage, production is begun at a constant rate, and well bottom pressure declines rapidly. When wellbore pressure declines to a limiting value, the producing rate begins to decline as the well produces against a constant pressure [32]. The compound working system developed in this work enables the prediction of long-term wells’ performance with a high level of accuracy. It is very consistent with the actual production situation of oil wells.


### 3.3.1. Shape and Size of Stimulated Reservoir Volume

The effects of SRV size and shape on the production rate are shown in Figures12 and 13, which demonstrate their significant influences on both the flow rate and cumulative production. Comparing the three models shown in Figure 12, the larger the SRV size is, the longer the constant flow period and the higher cumulative production. The larger size effect is to extend the constant period by about 166% and increase the cumulative production by about 30%. However, the rate difference becomes minimal after 500 days. This is because the SRV region’s oil yield has been primarily produced, leaving other regions to be the dominant factors for production.Figure 12
The effect of SRV size on the flow rate and cumulative yield (compound working system).Figure 13
The effect of SRV shape on flow rate and cumulative yield (compound working system).Figure13 shows that when the SRV region’s length-width ratio increases, the flow rate is higher. In other words, an SRV shape similar to the formation can help increase productivity. The effect of the x-direction rectangular shape of the SRV increases the cumulative production by about 56%. In the initial period, flow only happens in the hydraulic fractures and the stimulated region near the well. Therefore, producing rate is not affected by the shape and size of the SRV. In the exploitation process, the production wave spreads to the SRV region’s edges, and reservoirs with a similar shape to the SRV have a higher harvest. Therefore, in the fracturing design, the stimulated reservoir volume size and shape should be seriously discussed.


### 3.3.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=0.01,20,300) under the compound working system is shown in Figure 14. For a fracture conductivity of 0.01 (and other properties and parameters remaining constant), the production period at a constant rate lasts only 270 d. For fracture conductivities of 20 and 300, the period increases to 410 d, i.e., an increase of 51.9%. Moreover, the flow rate decreases rapidly in the constant pressure stage. These results demonstrate that the fracture conductivity helps extend the stable production period and achieve higher production within a short time.Figure 14
The effect of fracture conductivity on well performance (compound working system).The influence of the mobility ratio between the SRV region and matrix area on the dimensionless pressure is shown in Figure15. Here, the SRV permeability of k4 is varied while keeping the well and oil characteristics constant. It can be seen that with the value of mobility ratio increasing, the producing period at a constant rate extends by 150%, and the rate of pressure decline decreases. This is because improving the permeability of the SRV aids fluid flow.Figure 15
The effect of mobility ratio on well performance (compound working system).


### 3.3.3. Initial Flow Rate and Limiting Wellbore Pressure

Plots of producing rate and wellbore pressure vs. time for different initial flow rates (qoi=15,25,35m3/day) are shown in Figure 16. The producing rate decreases from 35 m3/day to 15 m3/day while the stable production period increases from 220 d to 1140 d. The effect of the lower initial flow rate is to extend the producing period at a constant rate of about 418%. Pressure decline at the wellbore steadily decreases during this period. Within a short period of time, a higher cumulative yield can be obtained. However, the ultimate recovery is not significantly different.Figure 16
Effect of initial flow rate on the performance curve (compound working system).Plots of producing rate and wellbore pressure vs. time for limiting wellbore pressure (pwfsc=6,8,10MPa) are shown in Figure 17. The conversion pressure decreases from 10 MPa to 6 MPa, while the stable production stage increases by 46% from 280 d to 410 d. The pressure drop rate is consistent, but the final output increases by 23% from 2.6856×104m3 to 3.2668×104m3. This shows that the conversion pressure can effectively extend the first period and significantly increase the well’s ultimate production at the same initial flow rate.Figure 17
Effect of limiting wellbore pressure on well performance curve (compound working system).


## 3.3.1. Shape and Size of Stimulated Reservoir Volume

The effects of SRV size and shape on the production rate are shown in Figures12 and 13, which demonstrate their significant influences on both the flow rate and cumulative production. Comparing the three models shown in Figure 12, the larger the SRV size is, the longer the constant flow period and the higher cumulative production. The larger size effect is to extend the constant period by about 166% and increase the cumulative production by about 30%. However, the rate difference becomes minimal after 500 days. This is because the SRV region’s oil yield has been primarily produced, leaving other regions to be the dominant factors for production.Figure 12
The effect of SRV size on the flow rate and cumulative yield (compound working system).Figure 13
The effect of SRV shape on flow rate and cumulative yield (compound working system).Figure13 shows that when the SRV region’s length-width ratio increases, the flow rate is higher. In other words, an SRV shape similar to the formation can help increase productivity. The effect of the x-direction rectangular shape of the SRV increases the cumulative production by about 56%. In the initial period, flow only happens in the hydraulic fractures and the stimulated region near the well. Therefore, producing rate is not affected by the shape and size of the SRV. In the exploitation process, the production wave spreads to the SRV region’s edges, and reservoirs with a similar shape to the SRV have a higher harvest. Therefore, in the fracturing design, the stimulated reservoir volume size and shape should be seriously discussed.


## 3.3.2. Fracture Conductivity and Mobility Ratio

The effect of fracture conductivity (CfD=0.01,20,300) under the compound working system is shown in Figure 14. For a fracture conductivity of 0.01 (and other properties and parameters remaining constant), the production period at a constant rate lasts only 270 d. For fracture conductivities of 20 and 300, the period increases to 410 d, i.e., an increase of 51.9%. Moreover, the flow rate decreases rapidly in the constant pressure stage. These results demonstrate that the fracture conductivity helps extend the stable production period and achieve higher production within a short time.Figure 14
The effect of fracture conductivity on well performance (compound working system).The influence of the mobility ratio between the SRV region and matrix area on the dimensionless pressure is shown in Figure15. Here, the SRV permeability of k4 is varied while keeping the well and oil characteristics constant. It can be seen that with the value of mobility ratio increasing, the producing period at a constant rate extends by 150%, and the rate of pressure decline decreases. This is because improving the permeability of the SRV aids fluid flow.Figure 15
The effect of mobility ratio on well performance (compound working system).


## 3.3.3. Initial Flow Rate and Limiting Wellbore Pressure

Plots of producing rate and wellbore pressure vs. time for different initial flow rates (qoi=15,25,35m3/day) are shown in Figure 16. The producing rate decreases from 35 m3/day to 15 m3/day while the stable production period increases from 220 d to 1140 d. The effect of the lower initial flow rate is to extend the producing period at a constant rate of about 418%. Pressure decline at the wellbore steadily decreases during this period. Within a short period of time, a higher cumulative yield can be obtained. However, the ultimate recovery is not significantly different.Figure 16
Effect of initial flow rate on the performance curve (compound working system).Plots of producing rate and wellbore pressure vs. time for limiting wellbore pressure (pwfsc=6,8,10MPa) are shown in Figure 17. The conversion pressure decreases from 10 MPa to 6 MPa, while the stable production stage increases by 46% from 280 d to 410 d. The pressure drop rate is consistent, but the final output increases by 23% from 2.6856×104m3 to 3.2668×104m3. This shows that the conversion pressure can effectively extend the first period and significantly increase the well’s ultimate production at the same initial flow rate.Figure 17
Effect of limiting wellbore pressure on well performance curve (compound working system).


## 4. Conclusions

This paper presents a semianalytical model to simulate the flow through a vertical well with branch fractures surrounded by a nonstimulated zone. The specific conclusions are as follows:(1)
A practical semianalytical model of ultralow-permeability reservoir with SRV is established based on the five-linear flow model and the fracture conductivity influence function(2)
The effects of critical parameters including SRV size and shape, fracture conductivity, mobility ratio, initial flow rate, and limiting wellbore pressure under different working systems are analysed(3)
The prediction results of constant production and constant pressure working systems show that a large SRV size can effectively stabilize oil production and reduce formation energy losses. Specifically, (i) an SRV similar in shape to the formation can stimulate the ultimate recovery rate, and (ii) for an SRV with constant shape and size, the larger the effective SRV, the better the well performance(4)
The prediction of long-term oil well performance under the compound working system shows that (i) SRV size, permeability, and initial flow rate can effectively prolong the constant rate period. Among them, the initial flow rate can extend by 418%. (ii) The SRV size and shape and conversion pressure can significantly improve the ultimate recovery rate. The SRV shape similar to the reservoir shape increases it by 56%



---

*Source: 5593663-2021-04-22.xml*