# An Analytical Method for Parameter Interpretation of Fracture Networks in Shale Gas Reservoirs considering Uneven Support of Fractures

Authors: Shijun Huang; Jiaojiao Zhang; Sidong Fang; Xifeng Wang
Journal: Geofluids
 (2021)

---

## Abstract

In shale gas reservoirs, the production data analysis method is widely used to invert reservoir and fracture parameter, and productivity prediction. Compared with numerical models and semianalytical models, which have high computational cost, the analytical model is mostly used in the production data analysis method to characterize the complex fracture network formed after fracturing. However, most of the current calculation models ignore the uneven support of fractures, and most of them use a single supported fracture model to describe the flow characteristics, which magnifies the role of supported fracture to a certain extent. Therefore, in this study, firstly, the fractures are divided into supported fractures and unsupported fractures. According to the near-well supported fractures and far-well unsupported fractures, the SRV zone is divided into outer SRV and inner SRV. The four areas are characterized by different seepage models, and the analytical solutions of the models are obtained by Laplace transform and inverse transform. Secondly, the material balance pseudotime is introduced to process the production data under the conditions of variable production and variable pressure. The double logarithmic curves of normalized production rate, rate integration, the derivative of the integration, and material balance pseudotime are established, and the parameters are interpreted by fitting the theoretical curve to the measured data. Then, the accuracy of the method is verified by comparison the parameter interpretation results with well test results, and the influence of parameters such as the half-length and permeability of supported and unsupported fractures on gas production is analyzed. Finally, the proposed method is applied to four field cases in southwest China. This paper mainly establishes an analytical method for parameter interpretation after hydraulic fracturing based on the production data analysis method considering the uneven support of fractures, which is of great significance for understanding the mechanism of fracturing stimulation, optimization of fracturing parameters, and gas production forecast.


---

## Body

## 1. Introduction

China is rich in shale gas resources and has a broad prospect of gas reserves and production. Shale formations are so tight that fracturing of horizontal wells is the key technique for their efficient development. And the parameter interpretation of stimulated reservoir and hydraulic fractures is the key to hydraulic fracturing evaluation [1].Parameter inversion methods can be divided into pressure dynamic inversion (well testing analysis method) and production dynamic inversion (production data analysis method). The shale formations are so tight that the well test requires a long time to shut in the well and measure the pressure, which seriously affects production [2–5]. The production data analysis method is similar to the well testing analysis method. It is based on seepage theory, introduces the material balance time and normalized rate to process production and pressure data, and inverts the reservoir and fracture parameters via theoretical curve and measured curve fitting [6–12]. The production data analysis method does not require closing the well. Through comprehensive analysis of the daily production data, fracture network parameters equivalent to the well test data can be obtained, which cannot only reduce production loss but can also make full use of the production data accumulated during the development process.After fracturing, parameters such as length, azimuth, and dip angle of each fracture are different in the stimulated reservoir volume (SRV) [13]; so, it is impossible to accurately describe each fracture. The first priority of parameter inversion is to describe the complex fracture network and establish a corresponding seepage model to analyze the production or pressure of fractured horizontal wells [14, 15]. According to the solution method, seepage models can be classified into three types of numerical solution model, semianalytical solution model, and analytical solution model. The stimulated reservoir area is divided into many small grids by the numerical model, and then the finite volume, finite element, and other numerical methods are used to solve the model. Discrete fracture model (DFM) [16–19] and embedded discrete fracture model (EDFM) [20–23] are the most widely used numerical models. The DFM deals with the fractures by the local mesh refinement approach. The large-scale fractures are extracted individually and dealt with discrete media when multiscale fractures are simulated in the EDFM. In most semianalytical models of fractured horizontal wells, source function models are used to treat fractures as linear sources, and the point source function is used to solve the model. Then, the superposition principle is applied to deal with multistage fractures [24–29]. However, when a large number of fractures are generated, the numerical models and semianalytical models require a large amount of computation and have low computational efficiency. At present, most parameter inversion methods based on production data analysis mostly adopt analytical models. In most analytical models, the SRV are divided into different areas to characterize the heterogeneity of the fractured reservoir and solved through coupling of the production and pressure. Bello and Wattenbarger [30] used the fracture matrix dual model to approximate the productivity of multistage fractured horizontal shale gas wells. Based on the dual media model, Alahmadi and Wattenbarger [31] further considered the secondary fractures in the reservoir and established a three-hole model. Brown et al. [32] considered the contribution of the unreconstructed area outside of the reformed area to the production and established a classic trilinear flow model. Su et al. [33] regarded the flow perpendicular to the hydraulic fracture as compound flow with internal double holes and homogeneous external areas, and they extended the three linear flow model to the four linear flow model. Stalgorova and Louis [34] further considered the unreconstructed area between the fractures, divided the model into five areas, and established the five linear flow model. Zeng et al. [35] extended the five-linear flow model and established the seven linear flow model. These analytical models take into account the heterogeneity of the reservoir after hydraulic fracturing, but ignore the uneven support of fractures. The length and conductivity of the fractures obtained based on these seepage models are the equivalent parameters of supported fractures and unsupported fractures. In fact, affected by factors such as ground stress, rock strength, and the degree of development of natural fractures, various fractures including proppant-filled fractures and unfilled fractures are formed after hydraulic fracturing [36, 37]. Previous studies on the conductivity of propped fractures by means of experiments and numerical simulations have shown that the conductivity of propped fractures plays a major role in the productivity of gas wells [38, 39]. However, an increasing number of studies have shown that unpropped fractures have conductivity even after closure, which also plays an important role in productivity [40]. And the more complex the fracture network is, the smaller the average proppant concentration in the fracture, and the majority of fractures are unpropped.In this study, the complex fracture network after fracturing is equivalent to the SRV area, and the SRV area is divided into the outer SRV and inner SRV based on the uneven support of the fractures. First, a seepage model considering the shale gas adsorption and gas physical character nonlinearity is established. The analytical solution of the model is obtained by Laplace transform and Stehfest numerical inversion. Then, a corresponding interpretation method of key parameters after fracturing based on production data analysis method is established. Finally, the proposed method was applied to four field examples which is provided to demonstrate the application of the proposed analytical method for parameter interpretation, and six parameters were obtained, including the half-lengths and permeabilities of the supported and unsupported fractures and the permeabilities of the outer and inner SRV, which can quantitatively distinguish between supported fractures and unsupported fractures.


## 2. Methodology

### 2.1. Physical Model

Due to the limitations of fracturing technology, it is difficult to evenly distribute the proppant in the fractures. Current studies believe that the near-well zone is the main fracture, which is filled with proppant, and most of the fractures in the far-well zone are not filled with proppant [41]. For shale gas wells in the early stage of production, the unreconstructed area outside of the stimulated area has a small contribution to gas production. In order to reduce the computation cost of parameter inversion and the multiplicity of solutions, the model in this study does not consider the unreconstructed region. According to the distribution of proppants in the hydraulic fracture, we divide the fractures into supported fracture, which is filled with proppant and unsupported fracture, which is unfilled with proppant. And the SRV area is divided into the outer SRV and inner SRV based on the near-well supported fractures and the far-well unsupported fractures (Figure 1), and we assume that the flow in each area is linear. The other assumptions are as follows. (1) The reservoir is closed vertically, and the horizontal well is located in the center of the reservoir. (2) The gas flows in a single phase and in an unsteady state within the stimulated reservoir and fractures. (3) The desorption of the adsorbed gas in the SRV area obeys Langmuir theory. (4) The fluid in the SRV area only flows linearly into the fractures, and then it flows into the wellbore from the fractures.Figure 1
Schematic diagram of the fractured horizontal well model considering the uneven support of fractures.


### 2.2. Diffusivity Equation in SRV and Fractures

#### 2.2.1. Outer SRV

In the outer SRV, the continuity equation is given below:(1)∂∂y⋅−ρvom+ρqa=∂ρϕom∂t.In Equation (1), qa is desorption rate of adsorbed gas per unit volume per unit time, vom is the gas velocity, ρ is gas density, and they are defined as
(2)qa=−Bg∂V∂t,(3)vom=−komμ∂pom∂y,(4)ρ=pomMzRT,Volume factor is defined as(5)Bg=psczTpzscTsc.The adsorption model obeys the Langmuir isotherm:(6)V=VLpompom+pL.Pseudopressure is defined as(7)ψom=2∫0pompomμzdpom.There is not only free gas but also much absorbed gas in the shale reservoirs.ctom is the modified compressibility defined by Gerami et al. [42] considering the desorption effect, and it is defined as follows:
(8)ctom=com+cg+cd.In Equation (7), com is the rock compressibility, cg is the gas compressibility, and cd is the desorption compressibility. Their definition is as follows:
(9)com=1ϕomdϕomdpom,cg=1pom−1zdzdpom,cd=psczTpomzscTscVLpLϕompL+pom2.Substituting Equations (2)–(11) into Equation (1) and adopt the engineering unit system, to arrive at
(10)∂2ψom∂y2=ϕomμgctom0.0864kom∂ψom∂t.The fluid in the outer SRV flows into the unsupported fracture, and the boundary conditions are(11)ψomy,tt=0=ψi,∂ψomy,t∂yy=ye=0,ψomy,ty=woF/2=ψoF.


#### 2.2.2. Inner SRV

Similarly, the diffusivity equation in the inner SRV can be deduced as follows:(12)∂2ψIm∂y2=ϕImμgctIm0.0864kIm∂ψIm∂t.The boundary conditions are(13)ψImy,tt=0=ψi,∂ψImy,t∂yy=ye=0,ψImy,ty=wIF/2=ψIF.


#### 2.2.3. Unsupported Fracture

Shale gas flow linearly from the outer SRV into the unsupported fracture and the continuity equation of the unsupported fracture are given below:(14)∂∂x⋅−ρvoF+Qmf=∂ρϕoF∂t.Velocity is defined as(15)voF=−koFμ∂poF∂x.Equation (14) is the mass of the gas flowing into the unsupported fracture from per unit volume of the outer SRV per unit time, and its definition is
(16)Qmf=−2ρkomωoFμ∂pom∂yy=ωoF2.Substituting Equations (4), (7), (15), and (16) into Equation (14), the diffusivity equation becomes
(17)∂2ψoF∂x2=10.0864ϕctμoFkoF∂ψoF∂t−2komwoFkoF∂ψom∂yx=ωoF2.The gas flow from unsupported fracture into supported fracture and the unsupported and supported fractures are coupled through the rate and pressure. Therefore, the boundary conditions are as follows:(18)ψoFx,tt=0=ψi,ψoFx=xIF=ψIFx=xIF,koFwoF∂ψoF∂xx=xIF=kIFwIF∂ψIF∂xx=xIF.


#### 2.2.4. Supported Fracture

Similarly, the diffusivity equation in the supported fracture can be obtained as follows:(19)∂2ψIF∂x2=10.0864ϕctμIFkIF∂ψIF∂t−2kImwIFkIF∂ψIm∂yx=ωIF2.The horizontal well is produced at constant pressure, and the boundary conditions are as follows:(20)ψIFx,tt=0=ψi,ψoFx=xIF=ψIFx=xIF,koFwoF∂ψoF∂xx=xIF=kIFwIF∂ψIF∂xx=xIF,ψIFx=0=ψwf.When the shale gas fractured horizontal well is produced at constant pressure, the productivity is(21)qg=0.0864TsczscpsckIFwIFTH∂ψIF∂xx=0.


### 2.3. Dimensionless Mathematical Models

The dimensionless mathematical models of the outer and inner SRV and unsupported and supported fractures are shown in Equations (24)–(27), and the dimensionless parameters are shown in Table1.
(22)∂2ψomD∂yD2=1ηomD∂ψomD∂tD,ψomDyD,tDtD=0=0,∂ψomDyD,tD∂xDyD=yeD=0,ψomDyD,tDyD=woFD/2=ψoFD,∂2ψImD∂yD2=1ηImD∂ψImD∂tD,ψImDyD,tDtD=0=0,∂ψImDyD,tD∂xDyD=yeD=0,ψImDyD,tDyD=wIFD/2=ψIFD,∂2ψoFD∂xD2=1ηoFD∂ψoFD∂tD−2CoFD∂ψomD∂yyD=woFD/2,ψoFDxD,tDtD=0=0,∂ψoFD∂xDxD=xoFD=0,ψoFDxD=xIFD=ψIFDxD=xIFD,∂2ψIFD∂xD2=1ηIFD∂ψIFD∂tD−2CIFD∂ψImD∂yyD=wIFD/2,ψIFDxD,tDtD=0=0,koFDwoFD∂ψoFD∂xDxD=xIFD=kIFDwIFD∂ψIFD∂xDxD=xIFD,ψIFDxD=0=1.Table 1
Definition of the dimensionless variables.
Dimensionless variablesDefinitionDimensionless rate at constant pressure1/qD=krHψi−ψwf/1.291×10−3qscTDimensionless pseudo pressureψD=ψi−ψξ/ψi−ψwfDimensionless timetD=ηr/Lr2tDimensionless width of an unsupported fractureωoFD=ωoF/LrDimensionless width of a supported fractureωIFD=ωIF/LrDimensionless half-width of the fractured sectionyeD=ye/LrDimensionless half-length of an unsupported fracturexoFD=xoF/LrDimensionless half-length of a supported fracturexIFD=xIF/LrDimensionless permeabilitykξD=kξ/krDimensionless pressure transmitting coefficientηξD=0.0864kξ/ϕξctξμgηrDimensionless unsupported fracture conductivityCoFD=koFDwoFD/komDDimensionless supported fracture conductivityCIFD=kIFDwIFD/kImDThe dimensionless rate of shale gas fractured horizontal well is(23)qgD=−kIFDwIFDπ∂ψIFD∂xDxD=0.The values ofξ are om (outer SRV),oF (unsupported fracture),Im (inner SRV), IF (supported fracture), and r (reference).


### 2.4. Mathematical Model Linearization

As in Equation (10), μgctm is a function of pressure, and it is nonlinear. Therefore, pseudotime [43–46] is used to linearize Equation (10) as follows:
(24)ta=∫0tμgictmiμgp¯mctmp¯mdt.By substituting Equation (24) into Equation (10), the diffusivity equation becomes
(25)∂2ψom∂y2=ϕomμgictomi0.0864kom∂ψom∂ta.The dimensionless pseudotime is defined as(26)tD=ηrLr2ta.Therefore, the linearized dimensionless outer SRV mathematical model becomes(27)∂2ψomD∂yD2=1ηomD∂ψomD∂taD,ψomDyD,taDtaD=0=0,∂ψomDyD,taD∂xDyD=yeD=0,ψomDyD,taDyD=woFD/2=ψoFD.The linearized dimensionless inner SRV model is(28)∂2ψImD∂yD2=1ηImD∂ψImD∂taD,ψImDyD,taDtaD=0=0,∂ψImDyD,taD∂xDyD=yeD=0,ψImDyD,taDyD=wIFD/2=ψIFD.


## 2.1. Physical Model

Due to the limitations of fracturing technology, it is difficult to evenly distribute the proppant in the fractures. Current studies believe that the near-well zone is the main fracture, which is filled with proppant, and most of the fractures in the far-well zone are not filled with proppant [41]. For shale gas wells in the early stage of production, the unreconstructed area outside of the stimulated area has a small contribution to gas production. In order to reduce the computation cost of parameter inversion and the multiplicity of solutions, the model in this study does not consider the unreconstructed region. According to the distribution of proppants in the hydraulic fracture, we divide the fractures into supported fracture, which is filled with proppant and unsupported fracture, which is unfilled with proppant. And the SRV area is divided into the outer SRV and inner SRV based on the near-well supported fractures and the far-well unsupported fractures (Figure 1), and we assume that the flow in each area is linear. The other assumptions are as follows. (1) The reservoir is closed vertically, and the horizontal well is located in the center of the reservoir. (2) The gas flows in a single phase and in an unsteady state within the stimulated reservoir and fractures. (3) The desorption of the adsorbed gas in the SRV area obeys Langmuir theory. (4) The fluid in the SRV area only flows linearly into the fractures, and then it flows into the wellbore from the fractures.Figure 1
Schematic diagram of the fractured horizontal well model considering the uneven support of fractures.


## 2.2. Diffusivity Equation in SRV and Fractures

### 2.2.1. Outer SRV

In the outer SRV, the continuity equation is given below:(1)∂∂y⋅−ρvom+ρqa=∂ρϕom∂t.In Equation (1), qa is desorption rate of adsorbed gas per unit volume per unit time, vom is the gas velocity, ρ is gas density, and they are defined as
(2)qa=−Bg∂V∂t,(3)vom=−komμ∂pom∂y,(4)ρ=pomMzRT,Volume factor is defined as(5)Bg=psczTpzscTsc.The adsorption model obeys the Langmuir isotherm:(6)V=VLpompom+pL.Pseudopressure is defined as(7)ψom=2∫0pompomμzdpom.There is not only free gas but also much absorbed gas in the shale reservoirs.ctom is the modified compressibility defined by Gerami et al. [42] considering the desorption effect, and it is defined as follows:
(8)ctom=com+cg+cd.In Equation (7), com is the rock compressibility, cg is the gas compressibility, and cd is the desorption compressibility. Their definition is as follows:
(9)com=1ϕomdϕomdpom,cg=1pom−1zdzdpom,cd=psczTpomzscTscVLpLϕompL+pom2.Substituting Equations (2)–(11) into Equation (1) and adopt the engineering unit system, to arrive at
(10)∂2ψom∂y2=ϕomμgctom0.0864kom∂ψom∂t.The fluid in the outer SRV flows into the unsupported fracture, and the boundary conditions are(11)ψomy,tt=0=ψi,∂ψomy,t∂yy=ye=0,ψomy,ty=woF/2=ψoF.


### 2.2.2. Inner SRV

Similarly, the diffusivity equation in the inner SRV can be deduced as follows:(12)∂2ψIm∂y2=ϕImμgctIm0.0864kIm∂ψIm∂t.The boundary conditions are(13)ψImy,tt=0=ψi,∂ψImy,t∂yy=ye=0,ψImy,ty=wIF/2=ψIF.


### 2.2.3. Unsupported Fracture

Shale gas flow linearly from the outer SRV into the unsupported fracture and the continuity equation of the unsupported fracture are given below:(14)∂∂x⋅−ρvoF+Qmf=∂ρϕoF∂t.Velocity is defined as(15)voF=−koFμ∂poF∂x.Equation (14) is the mass of the gas flowing into the unsupported fracture from per unit volume of the outer SRV per unit time, and its definition is
(16)Qmf=−2ρkomωoFμ∂pom∂yy=ωoF2.Substituting Equations (4), (7), (15), and (16) into Equation (14), the diffusivity equation becomes
(17)∂2ψoF∂x2=10.0864ϕctμoFkoF∂ψoF∂t−2komwoFkoF∂ψom∂yx=ωoF2.The gas flow from unsupported fracture into supported fracture and the unsupported and supported fractures are coupled through the rate and pressure. Therefore, the boundary conditions are as follows:(18)ψoFx,tt=0=ψi,ψoFx=xIF=ψIFx=xIF,koFwoF∂ψoF∂xx=xIF=kIFwIF∂ψIF∂xx=xIF.


### 2.2.4. Supported Fracture

Similarly, the diffusivity equation in the supported fracture can be obtained as follows:(19)∂2ψIF∂x2=10.0864ϕctμIFkIF∂ψIF∂t−2kImwIFkIF∂ψIm∂yx=ωIF2.The horizontal well is produced at constant pressure, and the boundary conditions are as follows:(20)ψIFx,tt=0=ψi,ψoFx=xIF=ψIFx=xIF,koFwoF∂ψoF∂xx=xIF=kIFwIF∂ψIF∂xx=xIF,ψIFx=0=ψwf.When the shale gas fractured horizontal well is produced at constant pressure, the productivity is(21)qg=0.0864TsczscpsckIFwIFTH∂ψIF∂xx=0.


## 2.2.1. Outer SRV

In the outer SRV, the continuity equation is given below:(1)∂∂y⋅−ρvom+ρqa=∂ρϕom∂t.In Equation (1), qa is desorption rate of adsorbed gas per unit volume per unit time, vom is the gas velocity, ρ is gas density, and they are defined as
(2)qa=−Bg∂V∂t,(3)vom=−komμ∂pom∂y,(4)ρ=pomMzRT,Volume factor is defined as(5)Bg=psczTpzscTsc.The adsorption model obeys the Langmuir isotherm:(6)V=VLpompom+pL.Pseudopressure is defined as(7)ψom=2∫0pompomμzdpom.There is not only free gas but also much absorbed gas in the shale reservoirs.ctom is the modified compressibility defined by Gerami et al. [42] considering the desorption effect, and it is defined as follows:
(8)ctom=com+cg+cd.In Equation (7), com is the rock compressibility, cg is the gas compressibility, and cd is the desorption compressibility. Their definition is as follows:
(9)com=1ϕomdϕomdpom,cg=1pom−1zdzdpom,cd=psczTpomzscTscVLpLϕompL+pom2.Substituting Equations (2)–(11) into Equation (1) and adopt the engineering unit system, to arrive at
(10)∂2ψom∂y2=ϕomμgctom0.0864kom∂ψom∂t.The fluid in the outer SRV flows into the unsupported fracture, and the boundary conditions are(11)ψomy,tt=0=ψi,∂ψomy,t∂yy=ye=0,ψomy,ty=woF/2=ψoF.


## 2.2.2. Inner SRV

Similarly, the diffusivity equation in the inner SRV can be deduced as follows:(12)∂2ψIm∂y2=ϕImμgctIm0.0864kIm∂ψIm∂t.The boundary conditions are(13)ψImy,tt=0=ψi,∂ψImy,t∂yy=ye=0,ψImy,ty=wIF/2=ψIF.


## 2.2.3. Unsupported Fracture

Shale gas flow linearly from the outer SRV into the unsupported fracture and the continuity equation of the unsupported fracture are given below:(14)∂∂x⋅−ρvoF+Qmf=∂ρϕoF∂t.Velocity is defined as(15)voF=−koFμ∂poF∂x.Equation (14) is the mass of the gas flowing into the unsupported fracture from per unit volume of the outer SRV per unit time, and its definition is
(16)Qmf=−2ρkomωoFμ∂pom∂yy=ωoF2.Substituting Equations (4), (7), (15), and (16) into Equation (14), the diffusivity equation becomes
(17)∂2ψoF∂x2=10.0864ϕctμoFkoF∂ψoF∂t−2komwoFkoF∂ψom∂yx=ωoF2.The gas flow from unsupported fracture into supported fracture and the unsupported and supported fractures are coupled through the rate and pressure. Therefore, the boundary conditions are as follows:(18)ψoFx,tt=0=ψi,ψoFx=xIF=ψIFx=xIF,koFwoF∂ψoF∂xx=xIF=kIFwIF∂ψIF∂xx=xIF.


## 2.2.4. Supported Fracture

Similarly, the diffusivity equation in the supported fracture can be obtained as follows:(19)∂2ψIF∂x2=10.0864ϕctμIFkIF∂ψIF∂t−2kImwIFkIF∂ψIm∂yx=ωIF2.The horizontal well is produced at constant pressure, and the boundary conditions are as follows:(20)ψIFx,tt=0=ψi,ψoFx=xIF=ψIFx=xIF,koFwoF∂ψoF∂xx=xIF=kIFwIF∂ψIF∂xx=xIF,ψIFx=0=ψwf.When the shale gas fractured horizontal well is produced at constant pressure, the productivity is(21)qg=0.0864TsczscpsckIFwIFTH∂ψIF∂xx=0.


## 2.3. Dimensionless Mathematical Models

The dimensionless mathematical models of the outer and inner SRV and unsupported and supported fractures are shown in Equations (24)–(27), and the dimensionless parameters are shown in Table1.
(22)∂2ψomD∂yD2=1ηomD∂ψomD∂tD,ψomDyD,tDtD=0=0,∂ψomDyD,tD∂xDyD=yeD=0,ψomDyD,tDyD=woFD/2=ψoFD,∂2ψImD∂yD2=1ηImD∂ψImD∂tD,ψImDyD,tDtD=0=0,∂ψImDyD,tD∂xDyD=yeD=0,ψImDyD,tDyD=wIFD/2=ψIFD,∂2ψoFD∂xD2=1ηoFD∂ψoFD∂tD−2CoFD∂ψomD∂yyD=woFD/2,ψoFDxD,tDtD=0=0,∂ψoFD∂xDxD=xoFD=0,ψoFDxD=xIFD=ψIFDxD=xIFD,∂2ψIFD∂xD2=1ηIFD∂ψIFD∂tD−2CIFD∂ψImD∂yyD=wIFD/2,ψIFDxD,tDtD=0=0,koFDwoFD∂ψoFD∂xDxD=xIFD=kIFDwIFD∂ψIFD∂xDxD=xIFD,ψIFDxD=0=1.Table 1
Definition of the dimensionless variables.
Dimensionless variablesDefinitionDimensionless rate at constant pressure1/qD=krHψi−ψwf/1.291×10−3qscTDimensionless pseudo pressureψD=ψi−ψξ/ψi−ψwfDimensionless timetD=ηr/Lr2tDimensionless width of an unsupported fractureωoFD=ωoF/LrDimensionless width of a supported fractureωIFD=ωIF/LrDimensionless half-width of the fractured sectionyeD=ye/LrDimensionless half-length of an unsupported fracturexoFD=xoF/LrDimensionless half-length of a supported fracturexIFD=xIF/LrDimensionless permeabilitykξD=kξ/krDimensionless pressure transmitting coefficientηξD=0.0864kξ/ϕξctξμgηrDimensionless unsupported fracture conductivityCoFD=koFDwoFD/komDDimensionless supported fracture conductivityCIFD=kIFDwIFD/kImDThe dimensionless rate of shale gas fractured horizontal well is(23)qgD=−kIFDwIFDπ∂ψIFD∂xDxD=0.The values ofξ are om (outer SRV),oF (unsupported fracture),Im (inner SRV), IF (supported fracture), and r (reference).


## 2.4. Mathematical Model Linearization

As in Equation (10), μgctm is a function of pressure, and it is nonlinear. Therefore, pseudotime [43–46] is used to linearize Equation (10) as follows:
(24)ta=∫0tμgictmiμgp¯mctmp¯mdt.By substituting Equation (24) into Equation (10), the diffusivity equation becomes
(25)∂2ψom∂y2=ϕomμgictomi0.0864kom∂ψom∂ta.The dimensionless pseudotime is defined as(26)tD=ηrLr2ta.Therefore, the linearized dimensionless outer SRV mathematical model becomes(27)∂2ψomD∂yD2=1ηomD∂ψomD∂taD,ψomDyD,taDtaD=0=0,∂ψomDyD,taD∂xDyD=yeD=0,ψomDyD,taDyD=woFD/2=ψoFD.The linearized dimensionless inner SRV model is(28)∂2ψImD∂yD2=1ηImD∂ψImD∂taD,ψImDyD,taDtaD=0=0,∂ψImDyD,taD∂xDyD=yeD=0,ψImDyD,taDyD=wIFD/2=ψIFD.


## 3. Inversion Method for the Fracture Network Parameters

### 3.1. Derivation of the Mathematical Model

In order to facilitate the derivation, the model was transformed using the Laplace transform, and the solution was obtained in the Laplace domain [47, 48].The solution of the dimensionless rate in the Laplace domain at constant bottom hole flow pressure is(29)qD¯=−kIFDwIFDπαIFs⋅tanhαIF⋅xIFD−fsfstanhαIF⋅xIFD−1.αom, αIm, αoF, andαIF are defined as follows:
(30)αom=sηomD,αIm=sηImD,αoF=sηoFD−2αomCoFDtanhαomwoFD2−yeD,αIF=sηIFD−2αImCIFDtanhαImwIFD2−yeD.Furthermore, the Stehfest numerical inversion method [49] was used to obtain the solution in the time domain.


### 3.2. Parameter Interpretation

The general procedure of the production data analysis method is as follows: firstly, a seepage model in the SRV area is established according to actual situation of reservoir, and then the production data is fitted based on the established seepage model to obtain parameters of stimulated reservoir and fractures. During the production of shale gas wells, the gas production rate and bottom hole pressure are constantly changing. However, the model proposed above was established under the condition of a constant flowing bottom hole pressure; so, the actual rate curve and the theoretical curve cannot be directly fitted to interpret the parameters [50–52]. Blasingame et al. [53] introduced the material balance pseudotime to establish the equivalent relationship between the constant and the variable production rate and found that the curve for the constant pressure case and the curve for the constant rate case overlap after the material balance time conversion. Therefore, the material balance pseudotime can be used to convert the variable rate and the variable pressure production curves into constant pressure curves for parameter interpretation. We substituted the material balance pseudotime for the actual time, substituted the normalized rate for the actual rate, and established data processing methods for the variable rate or pressure cases.The material pseudotime is defined as follows:(31)tca=μctiqsc∫0tqscμctdt=−GctiqscμiZi∗2piψi−ψ¯,In Equation (15), Z∗ is the modified gas compressibility factor considering the influence of desorption, which is defined as follows:
(32)Z∗=Z1+pscTVLZ/p+pLϕTsc.The pseudopressure is used to normalize the rate:(33)qD=qscψi−ψwf.In order to reduce the influence of the noise of the production data and to reduce the nonuniqueness of the solution, the dimensionless rate integration and the derivative of the integration are applied, which are defined as follows:(34)qDi=∫0tcaqsc/Δψdtcatca,qDid=−dqsc/Δψidlntca=−tcadqsc/Δψidtca.Establish the chart of the relationship between normalized production rate, rate integration, and the derivative of the integration and material balance pseudotime. The parameters are interpreted by curve fitting. The theoretical curve is drawn according to the proposed mathematical model. The production data are processed, and then the measured curve is drawn.


## 3.1. Derivation of the Mathematical Model

In order to facilitate the derivation, the model was transformed using the Laplace transform, and the solution was obtained in the Laplace domain [47, 48].The solution of the dimensionless rate in the Laplace domain at constant bottom hole flow pressure is(29)qD¯=−kIFDwIFDπαIFs⋅tanhαIF⋅xIFD−fsfstanhαIF⋅xIFD−1.αom, αIm, αoF, andαIF are defined as follows:
(30)αom=sηomD,αIm=sηImD,αoF=sηoFD−2αomCoFDtanhαomwoFD2−yeD,αIF=sηIFD−2αImCIFDtanhαImwIFD2−yeD.Furthermore, the Stehfest numerical inversion method [49] was used to obtain the solution in the time domain.


## 3.2. Parameter Interpretation

The general procedure of the production data analysis method is as follows: firstly, a seepage model in the SRV area is established according to actual situation of reservoir, and then the production data is fitted based on the established seepage model to obtain parameters of stimulated reservoir and fractures. During the production of shale gas wells, the gas production rate and bottom hole pressure are constantly changing. However, the model proposed above was established under the condition of a constant flowing bottom hole pressure; so, the actual rate curve and the theoretical curve cannot be directly fitted to interpret the parameters [50–52]. Blasingame et al. [53] introduced the material balance pseudotime to establish the equivalent relationship between the constant and the variable production rate and found that the curve for the constant pressure case and the curve for the constant rate case overlap after the material balance time conversion. Therefore, the material balance pseudotime can be used to convert the variable rate and the variable pressure production curves into constant pressure curves for parameter interpretation. We substituted the material balance pseudotime for the actual time, substituted the normalized rate for the actual rate, and established data processing methods for the variable rate or pressure cases.The material pseudotime is defined as follows:(31)tca=μctiqsc∫0tqscμctdt=−GctiqscμiZi∗2piψi−ψ¯,In Equation (15), Z∗ is the modified gas compressibility factor considering the influence of desorption, which is defined as follows:
(32)Z∗=Z1+pscTVLZ/p+pLϕTsc.The pseudopressure is used to normalize the rate:(33)qD=qscψi−ψwf.In order to reduce the influence of the noise of the production data and to reduce the nonuniqueness of the solution, the dimensionless rate integration and the derivative of the integration are applied, which are defined as follows:(34)qDi=∫0tcaqsc/Δψdtcatca,qDid=−dqsc/Δψidlntca=−tcadqsc/Δψidtca.Establish the chart of the relationship between normalized production rate, rate integration, and the derivative of the integration and material balance pseudotime. The parameters are interpreted by curve fitting. The theoretical curve is drawn according to the proposed mathematical model. The production data are processed, and then the measured curve is drawn.


## 4. Field Case Studies

### 4.1. Model Validation

A shale gas well YY-1 in southwestern China was taken as an example for verification. The well test analysis method assumes that the fractures are all supported fractures. In order to compare the parameters interpreted by the proposed method with the well test method, the uneven support of fractures is ignored, and the SRV model proposed in Figure1 is simplified (Figure 2). YY-1 horizontal well section is 1502 m long, with 63 clusters of 23 fracturing stages. According to the gas testing, the initial formation pressure is 75 MPa. The thickness of high-quality shale is 38 m. Water saturation is 56%, matrix porosity is 5.15%, and pore compressibility is 3×10−4 MPa-1. Langmuir pressure is 6.5 MPa, and Langmuir volume is 3 m3/t obtained by the adsorption isotherm experiment. The proposed parameter interpretation method is used to analyze the gas production rate and flowing bottom hole pressure of well YY-1(Figure 3), and three parameters were obtained through curve fitting, including the half-length and permeability of the fracture and the permeability of the SRV area.Figure 2
Schematic diagram of the simplified model.Figure 3
Production data for gas well YY-1.The fitting curve is shown in Figure4, in which the abscissa is the pseudotime of the material balance, and the ordinate is the normalized rate, rate integral, and integral derivative. By comparing the inversion parameters with the values obtained from well testing (Table 2), it can be seen that the half-length and permeability of the fracture and the permeability of the stimulated reservoir obtained through inversion are in good agreement with the values obtained from well testing, and the simplified model fits the production data well. Therefore, the parameter interpretation method proposed in this study can be applied to the quantitative characterization of parameters after hydraulic fracturing in shale gas reservoirs.Figure 4
Type curve fitting of the production data for well YY-1.Table 2
Inversion parameters of well YY-1.
ParameterInversion parametersWell testing resultsHalf-length of the fracture (m)66.0660Permeability of the SRV (mD)1.44×10−51×10−5Permeability of the fracture (mD)252.28234.72


### 4.2. Effect of the Uneven Support of Fractures

The half-length and conductivity of hydraulic fractures and permeability of the SRV area are the key parameters for evaluation of fracturing result, which also have a significant impact on shale gas productivity forecast. In this section, we analyzed the influence of six parameters, including the half-lengths and permeabilities of the supported and unsupported fractures and the permeabilities of the outer and inner SRV, on the average gas production rate in the first year and the cumulative gas production in twenty years.The variable values are shown in Table3. A single factor analysis method is used to analyze the influence of each parameter on the production. The default value of half-lengths and permeabilities of the supported and unsupported fractures and the permeabilities of the outer and inner SRV are 100 m, 30 m, 10 mD, 1000 mD, 2×10−5 mD, and 1×10−4mD, respectively.Table 3
Value of influencing factors.
FactorsParameter valuesHalf-length of the unsupported fracture (m)80, 100, 120, 130, 150Half-length of the supported fracture (m)20, 30, 40, 50, 60Permeability of the outer SRV (mD)2×10−5, 4×10−5, 6×10−5, 8×10−5, 1×10−4Permeability of the inner SRV (mD)2×10−5, 4×10−5, 6×10−5, 8×10−5, 1×10−4Permeability of the unsupported fracture (mD)0.1, 1, 10, 100, 500Permeability of the supported fracture (mD)100, 300, 500, 800, 1000The effects of six parameters on gas production are shown in Figures5~10. Fracture conductivity or permeability is a very important parameter in fracturing. When the permeability of unsupported fracture varies in a wide range, it has a great influence on productivity. As it increases, the gas production increases, but the increasing amplitude slows down (Figure 5). As permeability of supported fracture increases, gas production increases with small gradient. This is because a given minimum supported fracture permeability is sufficient for fluid transport (Figure 6). An important purpose of fracturing stimulation is to increase the extent and permeability of the SRV. It can be seen from Figures 7 and 8, with the increase of the permeability of the outer SRV or inner SRV, the gas production increased, but the increase rate tends to be reduced. The half-length of fracture determines the range size of the SRV area, which has a great influence on productivity. As shown in Figures 9 and 10, along with the increase of the half-length of supported or unsupported fracture, the gas production rate and cumulative production increase. In particular, gas production increases almost linearly as the half-length of the supported fracture increases.Figure 5
Effects of the permeability of unsupported fractures on gas production rate (a) and cumulative production (b).
(a)(b)Figure 6
Effects of the permeability of supported fractures on gas production rate (a) and cumulative production (b).
(a)(b)Figure 7
Effects of the permeability of outer SRV on gas production rate (a) and cumulative production (b).
(a)(b)Figure 8
Effects of the permeability of inner SRV on gas production rate (a) and cumulative production (b).
(a)(b)Figure 9
Effects of the half-length of unsupported fractures on gas production rate (a) and cumulative production (b).
(a)(b)Figure 10
Effects of the half-length of supported fractures on gas production rate (a) and cumulative production (b).
(a)(b)In order to further study the effect of the half-length of supported and unsupported fractures on gas production, the half-length of fracture, which is the sum of the half-length of supported and unsupported fractures, is introduced, and the logarithmic curve of normalized production rate, rate integration, and the derivative of the integration and material balance pseudotime is drawn. Assuming that the half-length of the fracture is 100 m, the half-length of supported fracture is set at 20%, 40%, 60%, 80% and 100% of the half-length of the fracture, respectively, in Figure11. Assuming that the half-length of the supported fracture accounts for 50% of the half-length of the fracture, the half-length of the fracture is set at 80 m, 100 m, 120 m, 130 m, and 150 m, respectively, in Figure 12. As can be seen from the logarithmic curve, the fracture linear flow, matrix fracture bilinear flow, matrix linear flow, transition flow, and system boundary control flow appear from left to right. The value of xIF/xF mainly affects the matrix fracture bilinear flow stage, the matrix linear flow stage, and the boundary control flow stage. It mainly affects the occurrence time of the boundary control flow stage and has a greater impact on the gas production. When xIF/xF is constant, the value of xF affects the matrix linear flow stage and the boundary control flow stage, mainly affecting the appearance time of the matrix linear flow stage, and also has a great influence on the production.Figure 11
Effects of thex_IF/x_F on normalized gas production.Figure 12
Effects of thex_F on normalized gas production.


### 4.3. Field Examples

The stimulated reservoir and fracture parameter inversion method proposed in this study, which considers the uneven support of fractures, are used to fit the production data of four fractured horizontal wells in Southwest China. In order to illustrate the difference in fitting results of different models to the same production data, this section compares the two models: considering and ignoring the uneven support of fractures. The gas production and bottom hole flow pressure data of the four wells are shown in Figure13, the production fitting results are shown in Figures 14~17, and the parameter inversion results are shown in Tables 4 and 5.Figure 13
Production data for gas well YY-2 ~ YY-5.Figure 14
Type curve fitting of YY-2 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Figure 15
Type curve fitting of YY-3 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Figure 16
Type curve fitting of YY-4 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Figure 17
Type curve fitting of YY-5 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Table 4
Parameter inversion results for wells YY-2 ~ YY-5 (considering the uneven support of fractures).
WellHalf-length of the fracture (m)Permeability of the fracture (mD)Permeability of the stimulated reservoir (mD)Supported fractureUnsupported fractureSupported fractureUnsupported fractureInner SRVOuter SRVYY-219.4395.5493.6141.981.12×10−41.06×10−4YY-315.4378.3460.093.521.18×10−41.15×10−4YY-46.5641.3140.871.971.39×10−41.06×10-4YY-530.01102.45100.2310.961.15×10−41.09×10−4Table 5
Parameter inversion results for wells YY-2 ~ YY-5 (ignoring the uneven support of fractures).
WellHalf-length of the fracture (m)Permeability of the fracture (mD)Permeability of the stimulated reservoir (mD)YY-2106.2980.441.08×10−4YY-385.2544.211.17×10−4YY-446.8333.541.22×10−4YY-5116.0476.991.13×10−4It should be noted that when the uneven support of fractures is ignored, the integral derivative data of production is too noisy; so, only the normalized production and integral derivative curves of production are fitted. As we can see from the fitting curve (Figures14~17), when the production time is more than 100 days, except for the well YY-2 without considering the uneven support of fractures, the two models can obtain better fitting results in the matrix linear flow and boundary control flow stages. The fitting effect of production data in the early stage of fracture linear flow is not good, which is mainly related to the flowback of fracturing fluid in the fracture at the initial stage of production. However, the model considering the uneven support of fracture has a better fitting effect. Comparing the parameter inversion results in Tables 4 and 5, it can be seen that the inverted fracture parameters are the equivalent parameters of supported fractures and unsupported fractures when the uneven support of fractures is not considered. Although different models can fit the same production data, the fitting effects are different, and the key parameters of the inversion are also different. In order to accurately invert reservoir and fracture parameters and further understand the mechanism of fracturing stimulation measures, it is necessary to consider the influence of the uneven support of fractures in the production analysis of shale gas fractured horizontal wells.


## 4.1. Model Validation

A shale gas well YY-1 in southwestern China was taken as an example for verification. The well test analysis method assumes that the fractures are all supported fractures. In order to compare the parameters interpreted by the proposed method with the well test method, the uneven support of fractures is ignored, and the SRV model proposed in Figure1 is simplified (Figure 2). YY-1 horizontal well section is 1502 m long, with 63 clusters of 23 fracturing stages. According to the gas testing, the initial formation pressure is 75 MPa. The thickness of high-quality shale is 38 m. Water saturation is 56%, matrix porosity is 5.15%, and pore compressibility is 3×10−4 MPa-1. Langmuir pressure is 6.5 MPa, and Langmuir volume is 3 m3/t obtained by the adsorption isotherm experiment. The proposed parameter interpretation method is used to analyze the gas production rate and flowing bottom hole pressure of well YY-1(Figure 3), and three parameters were obtained through curve fitting, including the half-length and permeability of the fracture and the permeability of the SRV area.Figure 2
Schematic diagram of the simplified model.Figure 3
Production data for gas well YY-1.The fitting curve is shown in Figure4, in which the abscissa is the pseudotime of the material balance, and the ordinate is the normalized rate, rate integral, and integral derivative. By comparing the inversion parameters with the values obtained from well testing (Table 2), it can be seen that the half-length and permeability of the fracture and the permeability of the stimulated reservoir obtained through inversion are in good agreement with the values obtained from well testing, and the simplified model fits the production data well. Therefore, the parameter interpretation method proposed in this study can be applied to the quantitative characterization of parameters after hydraulic fracturing in shale gas reservoirs.Figure 4
Type curve fitting of the production data for well YY-1.Table 2
Inversion parameters of well YY-1.
ParameterInversion parametersWell testing resultsHalf-length of the fracture (m)66.0660Permeability of the SRV (mD)1.44×10−51×10−5Permeability of the fracture (mD)252.28234.72


## 4.2. Effect of the Uneven Support of Fractures

The half-length and conductivity of hydraulic fractures and permeability of the SRV area are the key parameters for evaluation of fracturing result, which also have a significant impact on shale gas productivity forecast. In this section, we analyzed the influence of six parameters, including the half-lengths and permeabilities of the supported and unsupported fractures and the permeabilities of the outer and inner SRV, on the average gas production rate in the first year and the cumulative gas production in twenty years.The variable values are shown in Table3. A single factor analysis method is used to analyze the influence of each parameter on the production. The default value of half-lengths and permeabilities of the supported and unsupported fractures and the permeabilities of the outer and inner SRV are 100 m, 30 m, 10 mD, 1000 mD, 2×10−5 mD, and 1×10−4mD, respectively.Table 3
Value of influencing factors.
FactorsParameter valuesHalf-length of the unsupported fracture (m)80, 100, 120, 130, 150Half-length of the supported fracture (m)20, 30, 40, 50, 60Permeability of the outer SRV (mD)2×10−5, 4×10−5, 6×10−5, 8×10−5, 1×10−4Permeability of the inner SRV (mD)2×10−5, 4×10−5, 6×10−5, 8×10−5, 1×10−4Permeability of the unsupported fracture (mD)0.1, 1, 10, 100, 500Permeability of the supported fracture (mD)100, 300, 500, 800, 1000The effects of six parameters on gas production are shown in Figures5~10. Fracture conductivity or permeability is a very important parameter in fracturing. When the permeability of unsupported fracture varies in a wide range, it has a great influence on productivity. As it increases, the gas production increases, but the increasing amplitude slows down (Figure 5). As permeability of supported fracture increases, gas production increases with small gradient. This is because a given minimum supported fracture permeability is sufficient for fluid transport (Figure 6). An important purpose of fracturing stimulation is to increase the extent and permeability of the SRV. It can be seen from Figures 7 and 8, with the increase of the permeability of the outer SRV or inner SRV, the gas production increased, but the increase rate tends to be reduced. The half-length of fracture determines the range size of the SRV area, which has a great influence on productivity. As shown in Figures 9 and 10, along with the increase of the half-length of supported or unsupported fracture, the gas production rate and cumulative production increase. In particular, gas production increases almost linearly as the half-length of the supported fracture increases.Figure 5
Effects of the permeability of unsupported fractures on gas production rate (a) and cumulative production (b).
(a)(b)Figure 6
Effects of the permeability of supported fractures on gas production rate (a) and cumulative production (b).
(a)(b)Figure 7
Effects of the permeability of outer SRV on gas production rate (a) and cumulative production (b).
(a)(b)Figure 8
Effects of the permeability of inner SRV on gas production rate (a) and cumulative production (b).
(a)(b)Figure 9
Effects of the half-length of unsupported fractures on gas production rate (a) and cumulative production (b).
(a)(b)Figure 10
Effects of the half-length of supported fractures on gas production rate (a) and cumulative production (b).
(a)(b)In order to further study the effect of the half-length of supported and unsupported fractures on gas production, the half-length of fracture, which is the sum of the half-length of supported and unsupported fractures, is introduced, and the logarithmic curve of normalized production rate, rate integration, and the derivative of the integration and material balance pseudotime is drawn. Assuming that the half-length of the fracture is 100 m, the half-length of supported fracture is set at 20%, 40%, 60%, 80% and 100% of the half-length of the fracture, respectively, in Figure11. Assuming that the half-length of the supported fracture accounts for 50% of the half-length of the fracture, the half-length of the fracture is set at 80 m, 100 m, 120 m, 130 m, and 150 m, respectively, in Figure 12. As can be seen from the logarithmic curve, the fracture linear flow, matrix fracture bilinear flow, matrix linear flow, transition flow, and system boundary control flow appear from left to right. The value of xIF/xF mainly affects the matrix fracture bilinear flow stage, the matrix linear flow stage, and the boundary control flow stage. It mainly affects the occurrence time of the boundary control flow stage and has a greater impact on the gas production. When xIF/xF is constant, the value of xF affects the matrix linear flow stage and the boundary control flow stage, mainly affecting the appearance time of the matrix linear flow stage, and also has a great influence on the production.Figure 11
Effects of thex_IF/x_F on normalized gas production.Figure 12
Effects of thex_F on normalized gas production.


## 4.3. Field Examples

The stimulated reservoir and fracture parameter inversion method proposed in this study, which considers the uneven support of fractures, are used to fit the production data of four fractured horizontal wells in Southwest China. In order to illustrate the difference in fitting results of different models to the same production data, this section compares the two models: considering and ignoring the uneven support of fractures. The gas production and bottom hole flow pressure data of the four wells are shown in Figure13, the production fitting results are shown in Figures 14~17, and the parameter inversion results are shown in Tables 4 and 5.Figure 13
Production data for gas well YY-2 ~ YY-5.Figure 14
Type curve fitting of YY-2 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Figure 15
Type curve fitting of YY-3 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Figure 16
Type curve fitting of YY-4 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Figure 17
Type curve fitting of YY-5 production data considering (a) and ignoring (b) the uneven support of fractures.
(a)(b)Table 4
Parameter inversion results for wells YY-2 ~ YY-5 (considering the uneven support of fractures).
WellHalf-length of the fracture (m)Permeability of the fracture (mD)Permeability of the stimulated reservoir (mD)Supported fractureUnsupported fractureSupported fractureUnsupported fractureInner SRVOuter SRVYY-219.4395.5493.6141.981.12×10−41.06×10−4YY-315.4378.3460.093.521.18×10−41.15×10−4YY-46.5641.3140.871.971.39×10−41.06×10-4YY-530.01102.45100.2310.961.15×10−41.09×10−4Table 5
Parameter inversion results for wells YY-2 ~ YY-5 (ignoring the uneven support of fractures).
WellHalf-length of the fracture (m)Permeability of the fracture (mD)Permeability of the stimulated reservoir (mD)YY-2106.2980.441.08×10−4YY-385.2544.211.17×10−4YY-446.8333.541.22×10−4YY-5116.0476.991.13×10−4It should be noted that when the uneven support of fractures is ignored, the integral derivative data of production is too noisy; so, only the normalized production and integral derivative curves of production are fitted. As we can see from the fitting curve (Figures14~17), when the production time is more than 100 days, except for the well YY-2 without considering the uneven support of fractures, the two models can obtain better fitting results in the matrix linear flow and boundary control flow stages. The fitting effect of production data in the early stage of fracture linear flow is not good, which is mainly related to the flowback of fracturing fluid in the fracture at the initial stage of production. However, the model considering the uneven support of fracture has a better fitting effect. Comparing the parameter inversion results in Tables 4 and 5, it can be seen that the inverted fracture parameters are the equivalent parameters of supported fractures and unsupported fractures when the uneven support of fractures is not considered. Although different models can fit the same production data, the fitting effects are different, and the key parameters of the inversion are also different. In order to accurately invert reservoir and fracture parameters and further understand the mechanism of fracturing stimulation measures, it is necessary to consider the influence of the uneven support of fractures in the production analysis of shale gas fractured horizontal wells.


## 5. Conclusion

In this study, a simple and general inversion method of reservoir and fracture parameters is established based on the production data analysis method considering the uneven support of fractures. We can draw some conclusions as follows:(1)
Well test analysis and production data analysis are the two most commonly used method for fracture network parameter inversion. The production data analysis method can be better used for parameter interpretation because of its easy data collection and simple process in shale gas reservoirs(2)
Six parameters can be interpreted by curve fitting, including the half-lengths and permeabilities of the supported and unsupported fractures and the permeabilities of the outer and inner SRV, and thus the method proposed can quantitatively distinguish between supported fractures and unsupported fractures(3)
In the parameter inversion, the production data can be well fitted either with or without uneven support of fractures, but the model considering the uneven support of fracture has a better fitting effect. If the uneven support of fracture is not considered, the predicted output in the early stage of production will be much higher than the field data(4)
The pseudotime and pseudopressure can deal with the nonlinearity of gas physical characteristics in the model. And a method for processing production data with variable production and pressure can be established by introducing material balance pseudo time and normalized production rate



---

*Source: 3800525-2021-10-14.xml*