# Characteristic Value Method of Well Test Analysis for Horizontal Gas Well

Authors: Xiao-Ping Li; Ning-Ping Yan; Xiao-Hua Tan
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

This paper presents a study of characteristic value method of well test analysis for horizontal gas well. Owing to the complicated seepage flow mechanism in horizontal gas well and the difficulty in the analysis of transient pressure test data, this paper establishes the mathematical models of well test analysis for horizontal gas well with different inner and outer boundary conditions. On the basis of obtaining the solutions of the mathematical models, several type curves are plotted with Stehfest inversion algorithm. For gas reservoir with closed outer boundary in vertical direction and infinite outer boundary in horizontal direction, while considering the effect of wellbore storage and skin effect, the pseudopressure behavior of the horizontal gas well can manifest four characteristic periods: pure wellbore storage period, early vertical radial flow period, early linear flow period, and late horizontal pseudoradial flow period. For gas reservoir with closed outer boundary both in vertical and horizontal directions, the pseudopressure behavior of the horizontal gas well adds the pseudosteady state flow period which appears after the boundary response. For gas reservoir with closed outer boundary in vertical direction and constant pressure outer boundary in horizontal direction, the pseudopressure behavior of the horizontal gas well adds the steady state flow period which appears after the boundary response. According to the characteristic lines which are manifested by pseudopressure derivative curve of each flow period, formulas are developed to obtain horizontal permeability, vertical permeability, skin factor, reservoir pressure, and pore volume of the gas reservoir, and thus the characteristic value method of well test analysis for horizontal gas well is established. Finally, the example study verifies that the new method is reliable. Characteristic value method of well test analysis for horizontal gas well makes the well test analysis process more simple and the results more accurate.


---

## Body

## 1. Introduction

Recent years have seen the ever-growing application of horizontal wells technology, which aroused considerable interest in the exploration of horizontal well test analysis [1–4]. In order to surmount the challenges in estimating horizontal well productivity and parameters, analytical solutions for interpreting transient pressure behavior of horizontal wells have attracted great attention.Numerous studies on the pressure transient analysis of horizontal wells have been documented extensively in the literature. Combined with Newman’s product method, Gringarten and Ramey [5] found an access to solve the unsteady-flow problems in reservoirs by means of the use of source and Green’s function. Clonts and Ramey [6] presented an analytical solution for interpreting the transient pressure behavior of horizontal drain holes located in the heterogeneous reservoir. On the basis of finite Fourier transforms, Goode and Thambynayagam [7] addressed a solution for horizontal wells with infinite-conductivity in the semi-infinite reservoir. Ozkan and Rajagopal [8] demonstrated a derivative approach to analyze the pressure-transient behavior of horizontal wells, which revealed the relationship between the dimensionless well length and the horizontal-well pressure responses. Odeh and Babu [9] indicated that four significant flow periods could be observed during the process of horizontal well transient pressure behavior, which was further consolidated by the buildup and drawdown equations. Thompson and Temeng [10] introduced the automatic type curve matching method in analyzing multirate horizontal well pressure transient data through nonlinear regression analysis techniques. Rahavan et al. [11] employed a mathematical model to identify the features of pressure responses of a horizontal well with multiple fractures. Equipped with Laplace transformation and boundary element method, Zerzar and Bettam [12] addressed an analytical model for horizontal wells with finite conductivity vertical fractures. By extrapolating the transient pressure data, a simplified approach to predict well production was presented by Whittle et al. [13].Owing to the imperfection of common well test analysis methods including the semilog data plotting analysis technique [14, 15], type curve matching analysis method [16, 17], and automatic fitting analysis method [18], it is inconvenient to apply those methods during the process of analyzing and determining reservoir parameters. Therefore, this paper presents the characteristic value method of well test analysis for horizontal gas well for the sake of overcoming conventional limitations. This method involves two steps. The first step is to develop formulas to calculate gas reservoir fluid flow parameters according to the characteristic lines manifested by pseudopressure derivative curves of each flowing period. The next step is to utilize these formulas to complete the well test analysis for horizontal gas well by means of combining the measured pressure with the pseudopressure derivative curve. The characteristic value method of well test analysis for horizontal gas well enriches and develops the well test analysis theory and method.


## 2. Mathematical Models and Solutions of Well Test Analysis for Horizontal Gas Well

The hypothesis: the formation thickness ish, the initial formation pressure of gas reservoir is pi and equal everywhere, the gas reservoir is anisotropic, the horizontal permeability is Kh, the vertical permeability is Kv, horizontal section length is 2L, and the position of horizontal section in the gas reservoir which is parallel to the closed top and bottom boundary is zw. The surface flow rate of horizontal gas well is qsc and assumed to be constant. Single-phase compressible gas flow obeys Darcy law and the effect of gravity and capillary pressure is ignored. The physical model of horizontal gas well seepage is illustrated in Figure 1.Figure 1
Physical model of horizontal gas well seepage.Considering the complexity of the seepage flow mechanism of horizontal gas well and in order to make the mathematical model’s solving and calculation more simple, the establishment of mathematical models are divided into two parts: one is to ignore the effect of wellbore storage and skin effect; the other is to consider the effect of wellbore storage and skin effect [19, 20].


### 2.1. The Mathematical Models without Considering the Effect of Wellbore Storage and Skin

The diffusivity equation is expressed by Ozkan and Raghavan [21]:
(1)1rD∂∂rD(rD∂mD∂rD)+LD2∂2mD∂zD2=(hDLD)2∂mD∂tD.Initial condition is(2)mD(rD,0)=0.Inner boundary condition is(3)lim⁡ε→0[lim⁡rD→0∫zwD-ε/2zwD+ε/2rD∂mD∂rDdzwD]={0,zD>(zwD+ε2)-12,(zwD+ε2)≥zD≥(zwD-ε2)0,zD<(zwD-ε2),
where ε is a tiny variable.Infinite outer boundary condition in horizontal direction is(4)lim⁡rD→∞mD(rD,tD)=0.Closed outer boundary condition in horizontal direction is(5)∂mD∂rD|rD=reD=0.Constant pressure outer boundary condition in horizontal direction is(6)mD|rD=reD=0.Closed outer boundary conditions in vertical direction are(7)∂mD∂zD|zD=1=0,∂mD∂zD|zD=0=0.The dimensionless variables are defined as follows:(8)mD=78.489KhhqscT(mi-m),tD=3.6Khtϕμctrw2,LD=LhKvKh,hD=hrwKhKv,zD=zh,zrD=zwD+rwDLD,zwD=zwh,rD=rL,reD=reL,rwD=rwL.The defined gas pseudopressure is(9)m(p)=2∫prefppμ(p)Z(p)dp.


### 2.2. The Mathematical Model with Considering the Effect of Wellbore Storage and Skin

According to Duhamel’s principle [22] and the superposition principle, while using the definition of dimensionless variables, the mathematical model of horizontal gas well with considering the effect of wellbore storage and skin is derived as follows:
(10)mwD=mD+∫0tDCDdmwDdτDdmwD(tD-τD)dτDdτD+(1-CDdmwDdτD)hDS,
where
(11)CD=C2πϕCthL2.


### 2.3. The Solutions of the Mathematical Models

The solutions of the mathematical models [23, 24] at various outer boundary conditions can be obtained by applying source function and integral transform and taking the Laplace transform to s¯ with respect to tD.For gas reservoir with closed outer boundary in vertical direction and infinite outer boundary in horizontal direction, according to (1), (2), (3), (4), and (7), the dimensionless bottomhole pseudopressure of horizontal gas well in the Laplace space can be obtained. This results in
(12)m¯D=12s¯{∫-11K0((xD-α)2ε0)dα+2∑n=1∞∫-11K0((xD-α)2εn)cos⁡(βnzrD)×cos⁡(βnzwD)dα∑n=1∞∫-11K0((xD-α)2εn)cos(βnzrD)},
where
(13)βn=nπ,εn=s¯(hDLD)2+βnLD2.For gas reservoir with closed outer boundary both in vertical and horizontal direction, according to (1), (2), (3), (5), and (7), the dimensionless bottomhole pseudopressure of horizontal gas well in the Laplace space can be obtained. This results in
(14)m¯D=12s¯{∫-11K0((xD-α)2ε0)dα+K1(reDε0)I1(reDε0)∫-11I0((xD-α)2ε0)dα+2∑n=1∞[∫-11K0((xD-α)2εn)dα+K1(reDεn)I1(reDεn)∫-11I0((xD-α)2εn)dα·cos⁡(βnzrD)cos⁡(βnzwD)∫-11K0((xD-α)2εn)dα]}.For gas reservoir with closed outer boundary in vertical direction and constant pressure outer boundary in horizontal direction, according to (1), (2), (3), (6), and (7), the dimensionless bottomhole pseudopressure of horizontal gas well in Laplace space can be obtained. This results in
(15)m¯D=12s¯{∫-11K0((xD-α)2ε0)dα-K0(reDε0)I0(reDε0)∫-11I0((xD-α)2ε0)dα+2∑n=1∞[∫-11K0((xD-α)2εn)dα-K0(reDεn)I0(reDεn)∫-11I0((xD-α)2εn)dα·cos⁡(βnzrD)cos⁡(βnzwD)∫-11K0((xD-α)2εn)dα]}.Making the Laplace transform tos¯ with respect to tD/CD, (10) can be solved for the dimensionless bottomhole pseudopressure of horizontal gas well considering the effect of wellbore storage and skin in the Laplace space. This results in
(16)m¯wD=s¯m¯D+hDSs¯+s¯2(s¯m¯D+hDS)=1s¯(s¯+1/(s¯m¯D+hDS)).


## 2.1. The Mathematical Models without Considering the Effect of Wellbore Storage and Skin

The diffusivity equation is expressed by Ozkan and Raghavan [21]:
(1)1rD∂∂rD(rD∂mD∂rD)+LD2∂2mD∂zD2=(hDLD)2∂mD∂tD.Initial condition is(2)mD(rD,0)=0.Inner boundary condition is(3)lim⁡ε→0[lim⁡rD→0∫zwD-ε/2zwD+ε/2rD∂mD∂rDdzwD]={0,zD>(zwD+ε2)-12,(zwD+ε2)≥zD≥(zwD-ε2)0,zD<(zwD-ε2),
where ε is a tiny variable.Infinite outer boundary condition in horizontal direction is(4)lim⁡rD→∞mD(rD,tD)=0.Closed outer boundary condition in horizontal direction is(5)∂mD∂rD|rD=reD=0.Constant pressure outer boundary condition in horizontal direction is(6)mD|rD=reD=0.Closed outer boundary conditions in vertical direction are(7)∂mD∂zD|zD=1=0,∂mD∂zD|zD=0=0.The dimensionless variables are defined as follows:(8)mD=78.489KhhqscT(mi-m),tD=3.6Khtϕμctrw2,LD=LhKvKh,hD=hrwKhKv,zD=zh,zrD=zwD+rwDLD,zwD=zwh,rD=rL,reD=reL,rwD=rwL.The defined gas pseudopressure is(9)m(p)=2∫prefppμ(p)Z(p)dp.


## 2.2. The Mathematical Model with Considering the Effect of Wellbore Storage and Skin

According to Duhamel’s principle [22] and the superposition principle, while using the definition of dimensionless variables, the mathematical model of horizontal gas well with considering the effect of wellbore storage and skin is derived as follows:
(10)mwD=mD+∫0tDCDdmwDdτDdmwD(tD-τD)dτDdτD+(1-CDdmwDdτD)hDS,
where
(11)CD=C2πϕCthL2.


## 2.3. The Solutions of the Mathematical Models

The solutions of the mathematical models [23, 24] at various outer boundary conditions can be obtained by applying source function and integral transform and taking the Laplace transform to s¯ with respect to tD.For gas reservoir with closed outer boundary in vertical direction and infinite outer boundary in horizontal direction, according to (1), (2), (3), (4), and (7), the dimensionless bottomhole pseudopressure of horizontal gas well in the Laplace space can be obtained. This results in
(12)m¯D=12s¯{∫-11K0((xD-α)2ε0)dα+2∑n=1∞∫-11K0((xD-α)2εn)cos⁡(βnzrD)×cos⁡(βnzwD)dα∑n=1∞∫-11K0((xD-α)2εn)cos(βnzrD)},
where
(13)βn=nπ,εn=s¯(hDLD)2+βnLD2.For gas reservoir with closed outer boundary both in vertical and horizontal direction, according to (1), (2), (3), (5), and (7), the dimensionless bottomhole pseudopressure of horizontal gas well in the Laplace space can be obtained. This results in
(14)m¯D=12s¯{∫-11K0((xD-α)2ε0)dα+K1(reDε0)I1(reDε0)∫-11I0((xD-α)2ε0)dα+2∑n=1∞[∫-11K0((xD-α)2εn)dα+K1(reDεn)I1(reDεn)∫-11I0((xD-α)2εn)dα·cos⁡(βnzrD)cos⁡(βnzwD)∫-11K0((xD-α)2εn)dα]}.For gas reservoir with closed outer boundary in vertical direction and constant pressure outer boundary in horizontal direction, according to (1), (2), (3), (6), and (7), the dimensionless bottomhole pseudopressure of horizontal gas well in Laplace space can be obtained. This results in
(15)m¯D=12s¯{∫-11K0((xD-α)2ε0)dα-K0(reDε0)I0(reDε0)∫-11I0((xD-α)2ε0)dα+2∑n=1∞[∫-11K0((xD-α)2εn)dα-K0(reDεn)I0(reDεn)∫-11I0((xD-α)2εn)dα·cos⁡(βnzrD)cos⁡(βnzwD)∫-11K0((xD-α)2εn)dα]}.Making the Laplace transform tos¯ with respect to tD/CD, (10) can be solved for the dimensionless bottomhole pseudopressure of horizontal gas well considering the effect of wellbore storage and skin in the Laplace space. This results in
(16)m¯wD=s¯m¯D+hDSs¯+s¯2(s¯m¯D+hDS)=1s¯(s¯+1/(s¯m¯D+hDS)).


## 3. Type Curves of Well Test Analysis for Horizontal Gas Well

### 3.1. Gas Reservoir with Infinite Outer Boundary in Horizontal Direction

For gas reservoir with infinite outer boundary in horizontal direction, according to (12) which is the solution of horizontal well seepage mathematical model, while combining with (16), the type curve of well test analysis for horizontal gas well can be plotted with Stehfest inversion algorithm, as shown in Figure 2.Figure 2
Well test analysis type curve of horizontal well in gas reservoir with infinite outer boundary.As seen from Figure2, for gas reservoir with infinite outer boundary in horizontal direction, the pseudopressure behavior of horizontal gas well can manifest four characteristic periods: pure wellbore storage period (I), early vertical radial flow period (II), early linear flow period (III), and late horizontal pseudoradial flow period (IV).


#### 3.1.1. Pure Wellbore Storage Period

The characteristic of pure wellbore storage period of horizontal well is the same as vertical well, which is manifested as a 45° straight line segment on the log-log plot ofmwD, mwD′ versus tD/CD, and the duration of this period is affected by wellbore storage and skin effect.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(17)mwD=tDCD,mwD′=dmwDdln⁡(tD/CD)=tDCD.


#### 3.1.2. Early Vertical Radial Flow Period

The early vertical radial flow period appears after the effect of wellbore storage; the characteristic of this period is manifested as a horizontal straight line segment on the log-log plot ofmwD′ versus tD/CD. The pseudopressure behavior of this period is affected by formation thickness, horizontal section length, and the position of horizontal section in the gas reservoir. The flow regime of this period is shown in Figure 3.Figure 3
The schematic diagram of early vertical radial flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(18)mwD=14LD[ln⁡2.25(tDCD)+ln⁡CDe2S],mwD′=dmwDdln⁡(tD/CD)=14LD.


#### 3.1.3. Early Linear Flow Period

The early linear flow period appears after the early vertical radial flow period. The characteristic of this period is manifested as a straight line segment with a slope of 0.5 on the log-log plot ofmwD′ versus tD/CD. This characteristic describes the linear flow of fluid from formation to horizontal section. The pseudopressure behavior of this period is affected by dimensionless formation thickness hD, dimensionless horizontal section length LD, and the position of horizontal section in the gas reservoir zwD. The flow regime of this period is shown in Figure 4.Figure 4
The schematic diagram of early linear flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(19)mwD=2rwDπtD+S,mwD′=dmwDdln⁡(tD/CD)=rwDπtD.


#### 3.1.4. Late Horizontal Pseudoradial Flow Period

The late horizontal pseudoradial flow period appears after the early linear flow period. The characteristic of this period is manifested as a horizontal straight line segment with the value of 0.5 on the log-log plot ofmwD′ versus tD/CD. This characteristic describes the horizontal pseudoradial flow of fluid from formation horizontal plane in the distance to horizontal section. The flow regime of this period is shown in Figure 5.Figure 5
The schematic diagram of late horizontal pseudoradial flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(20)mwD=12[ln⁡(rwDtDCD)+ln⁡CDe2S],mwD′=dmwDdln⁡(tD/CD)=12.


### 3.2. Gas Reservoir with Closed Outer Boundary in Horizontal Direction

For gas reservoir with closed outer boundary in horizontal direction, according to (14) which is the solution of horizontal well seepage mathematical model, while combining with (16), the type curve of well test analysis for horizontal gas well can be plotted with Stehfest inversion algorithm, as shown in Figure 6.Figure 6
Well test analysis type curve of horizontal well in gas reservoir with closed outer boundary.As seen from Figure6, for gas reservoir with closed outer boundary in horizontal direction, the pseudopressure behavior of horizontal gas well can manifest five characteristic periods. The previous four periods of gas reservoir with closed outer boundary are exactly the same as gas reservoir with infinite outer boundary, but the pseudopressure behavior of the horizontal gas well adds the pseudosteady state flow period (V) which appears after the boundary response. The characteristic of the pseudosteady state flow period is manifested as a straight line segment with a slope of 1 on the log-log plot of mwD, mwD′ versus tD/CD. The greater the distance of the outer boundary, the later the appearance of the pseudosteady state flow period. The smaller the distance of the outer boundary, the sooner the appearance of the pseudosteady state flow period.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during the pseudosteady state flow period can be obtained. This results in(21)mwD=2π(rwDreD)2tD+S,mwD′=dmwDdln⁡(tD/CD)=2π(rwDreD)2tD.


### 3.3. Gas Reservoir with Constant Pressure Outer Boundary in Horizontal Direction

For gas reservoir with constant pressure outer boundary in horizontal direction, according to (15) which is the solution of horizontal well seepage mathematical model, while combining with (16), the type curve of well test analysis for horizontal gas well can be plotted with Stehfest inversion algorithm, as shown in Figure 7.Figure 7
Well test analysis type curve of horizontal well in gas reservoir with constant pressure outer boundary.As seen from Figure7, for gas reservoir with constant pressure outer boundary in horizontal direction, the pseudopressure behavior of horizontal gas well can manifest five characteristic periods; the previous four periods of gas reservoir with constant pressure outer boundary are exactly the same as gas reservoir with infinite outer boundary, but the pseudopressure behavior of the horizontal gas well adds the steady state flow period (V) which appears after the boundary response. The occurrence time of the steady state flow period is affected by the outer boundary distance in horizontal direction. The smaller the distance of the outer boundary, the sooner the appearance of the steady state flow period. The greater the distance of the outer boundary, the later the appearance of the steady state flow period.


### 3.4. Characteristic Value Method of Well Test Analysis for Horizontal Gas Well

The characteristic value method of well test analysis for horizontal gas well can determine the gas reservoir fluid flow parameters according to the characteristic lines which are manifested by pseudopressure derivative curve of each flow period on the log-log plot.


#### 3.4.1. Pure Wellbore Storage Period

The characteristic of pure wellbore storage period of horizontal well is manifested as a straight line segment with a slope of 1 on the log-log plot ofmwD, mwD′ versus tD/CD. The expression of dimensionless bottomhole pseudopressure during this period is
(22)mwD=tDCD.Equation (22) can be converted to dimensional form, and then according to the time and pressure data during pure wellbore storage period, the method to determine the wellbore storage coefficient can be obtained. By plotting the log-log plot of Δm, Δm′ versus t, the wellbore storage coefficient can be determined by the straight line segment with a slope of 1 on the log-log plot.The following can be obtained from the definitions of dimensionless variables:(23)tDCD=3.6Kht/ϕμctrw2C/2πϕcthL2=7.2πKhhL2μrw2tC.According to (22), (23), and the definition of dimensionless pseudopressure, the wellbore storage coefficient can be obtained. This results in
(24)C=0.288qscTμL2rw2tΔm,
where t/Δm represents the actual value on the log-log plot of Δm, Δm′ versus t during the pure wellbore storage period.


#### 3.4.2. Early Vertical Radial Flow Period

The Determination of Geometric Mean Permeability. The dimensionless bottomhole pseudopressure derivative curve is manifested as a horizontal straight line segment with the value of 1/(4LD) during the early vertical radial flow period. The expression of dimensionless bottomhole pseudopressure derivative during this period is
(25)mwD′=dmwDdln⁡(tD/CD)=14LD.According to the definitions of dimensionless variables, the dimensional form of (25) can be obtained. This results in
(26)78.489KhhqscT(tΔm′)er=14(L/h)Kv/Kh.The geometric mean permeability of gas reservoir can be determined by (26). This results in
(27)KhKv=3.185×10-3qscTL(tΔm′)er,
where (tΔm′)er represents the actual value on the log-log plot of Δm′ versus t during the early vertical radial flow period.The Determination of Skin Factor and Initial Reservoir Pressure. The expression of dimensionless bottomhole pseudopressure during the early vertical radial flow period is
(28)mwD=14LD[ln⁡2.25(tDCD)+ln⁡CDe2S].According to the definitions of dimensionless variables and (28), the skin factor can be obtained. This results in
(29)S=0.5[Δmer(tΔm′)er-ln⁡KhterϕμCtrw2-0.80907],
where Δmer and ter represent the pseudopressure difference and time corresponding to the (tΔm′)er, respectively.For pressure buildup analysis, whenΔt→∞, the ln⁡Δt/(Δt+tp)→0, (30) can be obtained through the use of the definitions of dimensionless variables and pseudopressure difference during the early vertical radial flow period:
(30)mi-mwftΔm′=ln⁡tpD+0.80907+2S.The initial reservoir pseudopressure can be determined by (30). This results in
(31)mi=mwf+(tΔm′)erb(ln⁡tpD+0.80907+2S),
where (tΔm′)erb represents the actual value on the pressure buildup log-log plot of Δm′ versus t during the early vertical radial flow period.


#### 3.4.3. Early Linear Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a straight line segment with the slope of 0.5 during the early linear flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, the following can be obtained:(32)78.489KhhqscT(tΔm′)l=rwL3.6πKhtϕμCtrw2.The horizontal permeability of gas reservoir can be determined by (32). This results in
(33)Kh=4.28×10-2qscTLhϕμCt[t(tΔm′)]l,
where (tΔm′)l represents the actual value on the log-log plot of Δm′ versus t during the early linear flow period.Combining (27) with (33), the vertical permeability can be obtained. This results in
(34)Kv=7.44×10-2hϕμCt[(tΔm′)/t]l(tΔm′)er.


#### 3.4.4. Late Horizontal Pseudoradial Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a horizontal straight line segment with the value of 0.5 during the late horizontal pseudoradial flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, (35) can be obtained:
(35)78.489KhhqscT(tΔm′)lr=0.5.The horizontal permeability of gas reservoir can be determined by (35). This results in
(36)Kh=6.37×10-3qscTh(tΔm′)lr,
where (tΔm′)lr represents the actual value on the log-log plot of Δm′ versus t during the late horizontal pseudoradial flow period.For pressure buildup analysis, whenΔt→∞, the ln⁡Δt/(Δt+tp)→0, (37) can be obtained through the use of the definitions of dimensionless variables and pseudopressure difference during the late horizontal pseudoradial flow period:
(37)mi-mwf(tΔm′)lrb=ln⁡rwDtpD+0.80907+2S.The initial reservoir pseudopressure can be determined by (37). This results in
(38)mi=mwf+(tΔm′)lrb(ln⁡rwDtpD+0.80907+2S),
where (tΔm′)lrb represents the actual value on the pressure buildup log-log plot of Δm′ versus t during the late horizontal pseudoradial flow period.


#### 3.4.5. Pseudosteady Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a straight line segment with the slope of 1 during the pseudosteady flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, (39) can be obtained:
(39)78.489KhhqscT(tΔm′)pp=7.2πKhtppre2hϕμCt.The pore volume of gas reservoir can be determined by (39). This results in
(40)πre2hϕ=0.905qscThμCttpp(tΔm′)pp,
where (tΔm′)pp represents the actual value on the log-log plot of Δm′ versus t during the pseudosteady flow period.


## 3.1. Gas Reservoir with Infinite Outer Boundary in Horizontal Direction

For gas reservoir with infinite outer boundary in horizontal direction, according to (12) which is the solution of horizontal well seepage mathematical model, while combining with (16), the type curve of well test analysis for horizontal gas well can be plotted with Stehfest inversion algorithm, as shown in Figure 2.Figure 2
Well test analysis type curve of horizontal well in gas reservoir with infinite outer boundary.As seen from Figure2, for gas reservoir with infinite outer boundary in horizontal direction, the pseudopressure behavior of horizontal gas well can manifest four characteristic periods: pure wellbore storage period (I), early vertical radial flow period (II), early linear flow period (III), and late horizontal pseudoradial flow period (IV).


### 3.1.1. Pure Wellbore Storage Period

The characteristic of pure wellbore storage period of horizontal well is the same as vertical well, which is manifested as a 45° straight line segment on the log-log plot ofmwD, mwD′ versus tD/CD, and the duration of this period is affected by wellbore storage and skin effect.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(17)mwD=tDCD,mwD′=dmwDdln⁡(tD/CD)=tDCD.


### 3.1.2. Early Vertical Radial Flow Period

The early vertical radial flow period appears after the effect of wellbore storage; the characteristic of this period is manifested as a horizontal straight line segment on the log-log plot ofmwD′ versus tD/CD. The pseudopressure behavior of this period is affected by formation thickness, horizontal section length, and the position of horizontal section in the gas reservoir. The flow regime of this period is shown in Figure 3.Figure 3
The schematic diagram of early vertical radial flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(18)mwD=14LD[ln⁡2.25(tDCD)+ln⁡CDe2S],mwD′=dmwDdln⁡(tD/CD)=14LD.


### 3.1.3. Early Linear Flow Period

The early linear flow period appears after the early vertical radial flow period. The characteristic of this period is manifested as a straight line segment with a slope of 0.5 on the log-log plot ofmwD′ versus tD/CD. This characteristic describes the linear flow of fluid from formation to horizontal section. The pseudopressure behavior of this period is affected by dimensionless formation thickness hD, dimensionless horizontal section length LD, and the position of horizontal section in the gas reservoir zwD. The flow regime of this period is shown in Figure 4.Figure 4
The schematic diagram of early linear flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(19)mwD=2rwDπtD+S,mwD′=dmwDdln⁡(tD/CD)=rwDπtD.


### 3.1.4. Late Horizontal Pseudoradial Flow Period

The late horizontal pseudoradial flow period appears after the early linear flow period. The characteristic of this period is manifested as a horizontal straight line segment with the value of 0.5 on the log-log plot ofmwD′ versus tD/CD. This characteristic describes the horizontal pseudoradial flow of fluid from formation horizontal plane in the distance to horizontal section. The flow regime of this period is shown in Figure 5.Figure 5
The schematic diagram of late horizontal pseudoradial flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(20)mwD=12[ln⁡(rwDtDCD)+ln⁡CDe2S],mwD′=dmwDdln⁡(tD/CD)=12.


## 3.1.1. Pure Wellbore Storage Period

The characteristic of pure wellbore storage period of horizontal well is the same as vertical well, which is manifested as a 45° straight line segment on the log-log plot ofmwD, mwD′ versus tD/CD, and the duration of this period is affected by wellbore storage and skin effect.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(17)mwD=tDCD,mwD′=dmwDdln⁡(tD/CD)=tDCD.


## 3.1.2. Early Vertical Radial Flow Period

The early vertical radial flow period appears after the effect of wellbore storage; the characteristic of this period is manifested as a horizontal straight line segment on the log-log plot ofmwD′ versus tD/CD. The pseudopressure behavior of this period is affected by formation thickness, horizontal section length, and the position of horizontal section in the gas reservoir. The flow regime of this period is shown in Figure 3.Figure 3
The schematic diagram of early vertical radial flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(18)mwD=14LD[ln⁡2.25(tDCD)+ln⁡CDe2S],mwD′=dmwDdln⁡(tD/CD)=14LD.


## 3.1.3. Early Linear Flow Period

The early linear flow period appears after the early vertical radial flow period. The characteristic of this period is manifested as a straight line segment with a slope of 0.5 on the log-log plot ofmwD′ versus tD/CD. This characteristic describes the linear flow of fluid from formation to horizontal section. The pseudopressure behavior of this period is affected by dimensionless formation thickness hD, dimensionless horizontal section length LD, and the position of horizontal section in the gas reservoir zwD. The flow regime of this period is shown in Figure 4.Figure 4
The schematic diagram of early linear flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(19)mwD=2rwDπtD+S,mwD′=dmwDdln⁡(tD/CD)=rwDπtD.


## 3.1.4. Late Horizontal Pseudoradial Flow Period

The late horizontal pseudoradial flow period appears after the early linear flow period. The characteristic of this period is manifested as a horizontal straight line segment with the value of 0.5 on the log-log plot ofmwD′ versus tD/CD. This characteristic describes the horizontal pseudoradial flow of fluid from formation horizontal plane in the distance to horizontal section. The flow regime of this period is shown in Figure 5.Figure 5
The schematic diagram of late horizontal pseudoradial flow.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during this period can be obtained. This results in(20)mwD=12[ln⁡(rwDtDCD)+ln⁡CDe2S],mwD′=dmwDdln⁡(tD/CD)=12.


## 3.2. Gas Reservoir with Closed Outer Boundary in Horizontal Direction

For gas reservoir with closed outer boundary in horizontal direction, according to (14) which is the solution of horizontal well seepage mathematical model, while combining with (16), the type curve of well test analysis for horizontal gas well can be plotted with Stehfest inversion algorithm, as shown in Figure 6.Figure 6
Well test analysis type curve of horizontal well in gas reservoir with closed outer boundary.As seen from Figure6, for gas reservoir with closed outer boundary in horizontal direction, the pseudopressure behavior of horizontal gas well can manifest five characteristic periods. The previous four periods of gas reservoir with closed outer boundary are exactly the same as gas reservoir with infinite outer boundary, but the pseudopressure behavior of the horizontal gas well adds the pseudosteady state flow period (V) which appears after the boundary response. The characteristic of the pseudosteady state flow period is manifested as a straight line segment with a slope of 1 on the log-log plot of mwD, mwD′ versus tD/CD. The greater the distance of the outer boundary, the later the appearance of the pseudosteady state flow period. The smaller the distance of the outer boundary, the sooner the appearance of the pseudosteady state flow period.Expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during the pseudosteady state flow period can be obtained. This results in(21)mwD=2π(rwDreD)2tD+S,mwD′=dmwDdln⁡(tD/CD)=2π(rwDreD)2tD.


## 3.3. Gas Reservoir with Constant Pressure Outer Boundary in Horizontal Direction

For gas reservoir with constant pressure outer boundary in horizontal direction, according to (15) which is the solution of horizontal well seepage mathematical model, while combining with (16), the type curve of well test analysis for horizontal gas well can be plotted with Stehfest inversion algorithm, as shown in Figure 7.Figure 7
Well test analysis type curve of horizontal well in gas reservoir with constant pressure outer boundary.As seen from Figure7, for gas reservoir with constant pressure outer boundary in horizontal direction, the pseudopressure behavior of horizontal gas well can manifest five characteristic periods; the previous four periods of gas reservoir with constant pressure outer boundary are exactly the same as gas reservoir with infinite outer boundary, but the pseudopressure behavior of the horizontal gas well adds the steady state flow period (V) which appears after the boundary response. The occurrence time of the steady state flow period is affected by the outer boundary distance in horizontal direction. The smaller the distance of the outer boundary, the sooner the appearance of the steady state flow period. The greater the distance of the outer boundary, the later the appearance of the steady state flow period.


## 3.4. Characteristic Value Method of Well Test Analysis for Horizontal Gas Well

The characteristic value method of well test analysis for horizontal gas well can determine the gas reservoir fluid flow parameters according to the characteristic lines which are manifested by pseudopressure derivative curve of each flow period on the log-log plot.


### 3.4.1. Pure Wellbore Storage Period

The characteristic of pure wellbore storage period of horizontal well is manifested as a straight line segment with a slope of 1 on the log-log plot ofmwD, mwD′ versus tD/CD. The expression of dimensionless bottomhole pseudopressure during this period is
(22)mwD=tDCD.Equation (22) can be converted to dimensional form, and then according to the time and pressure data during pure wellbore storage period, the method to determine the wellbore storage coefficient can be obtained. By plotting the log-log plot of Δm, Δm′ versus t, the wellbore storage coefficient can be determined by the straight line segment with a slope of 1 on the log-log plot.The following can be obtained from the definitions of dimensionless variables:(23)tDCD=3.6Kht/ϕμctrw2C/2πϕcthL2=7.2πKhhL2μrw2tC.According to (22), (23), and the definition of dimensionless pseudopressure, the wellbore storage coefficient can be obtained. This results in
(24)C=0.288qscTμL2rw2tΔm,
where t/Δm represents the actual value on the log-log plot of Δm, Δm′ versus t during the pure wellbore storage period.


### 3.4.2. Early Vertical Radial Flow Period

The Determination of Geometric Mean Permeability. The dimensionless bottomhole pseudopressure derivative curve is manifested as a horizontal straight line segment with the value of 1/(4LD) during the early vertical radial flow period. The expression of dimensionless bottomhole pseudopressure derivative during this period is
(25)mwD′=dmwDdln⁡(tD/CD)=14LD.According to the definitions of dimensionless variables, the dimensional form of (25) can be obtained. This results in
(26)78.489KhhqscT(tΔm′)er=14(L/h)Kv/Kh.The geometric mean permeability of gas reservoir can be determined by (26). This results in
(27)KhKv=3.185×10-3qscTL(tΔm′)er,
where (tΔm′)er represents the actual value on the log-log plot of Δm′ versus t during the early vertical radial flow period.The Determination of Skin Factor and Initial Reservoir Pressure. The expression of dimensionless bottomhole pseudopressure during the early vertical radial flow period is
(28)mwD=14LD[ln⁡2.25(tDCD)+ln⁡CDe2S].According to the definitions of dimensionless variables and (28), the skin factor can be obtained. This results in
(29)S=0.5[Δmer(tΔm′)er-ln⁡KhterϕμCtrw2-0.80907],
where Δmer and ter represent the pseudopressure difference and time corresponding to the (tΔm′)er, respectively.For pressure buildup analysis, whenΔt→∞, the ln⁡Δt/(Δt+tp)→0, (30) can be obtained through the use of the definitions of dimensionless variables and pseudopressure difference during the early vertical radial flow period:
(30)mi-mwftΔm′=ln⁡tpD+0.80907+2S.The initial reservoir pseudopressure can be determined by (30). This results in
(31)mi=mwf+(tΔm′)erb(ln⁡tpD+0.80907+2S),
where (tΔm′)erb represents the actual value on the pressure buildup log-log plot of Δm′ versus t during the early vertical radial flow period.


### 3.4.3. Early Linear Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a straight line segment with the slope of 0.5 during the early linear flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, the following can be obtained:(32)78.489KhhqscT(tΔm′)l=rwL3.6πKhtϕμCtrw2.The horizontal permeability of gas reservoir can be determined by (32). This results in
(33)Kh=4.28×10-2qscTLhϕμCt[t(tΔm′)]l,
where (tΔm′)l represents the actual value on the log-log plot of Δm′ versus t during the early linear flow period.Combining (27) with (33), the vertical permeability can be obtained. This results in
(34)Kv=7.44×10-2hϕμCt[(tΔm′)/t]l(tΔm′)er.


### 3.4.4. Late Horizontal Pseudoradial Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a horizontal straight line segment with the value of 0.5 during the late horizontal pseudoradial flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, (35) can be obtained:
(35)78.489KhhqscT(tΔm′)lr=0.5.The horizontal permeability of gas reservoir can be determined by (35). This results in
(36)Kh=6.37×10-3qscTh(tΔm′)lr,
where (tΔm′)lr represents the actual value on the log-log plot of Δm′ versus t during the late horizontal pseudoradial flow period.For pressure buildup analysis, whenΔt→∞, the ln⁡Δt/(Δt+tp)→0, (37) can be obtained through the use of the definitions of dimensionless variables and pseudopressure difference during the late horizontal pseudoradial flow period:
(37)mi-mwf(tΔm′)lrb=ln⁡rwDtpD+0.80907+2S.The initial reservoir pseudopressure can be determined by (37). This results in
(38)mi=mwf+(tΔm′)lrb(ln⁡rwDtpD+0.80907+2S),
where (tΔm′)lrb represents the actual value on the pressure buildup log-log plot of Δm′ versus t during the late horizontal pseudoradial flow period.


### 3.4.5. Pseudosteady Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a straight line segment with the slope of 1 during the pseudosteady flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, (39) can be obtained:
(39)78.489KhhqscT(tΔm′)pp=7.2πKhtppre2hϕμCt.The pore volume of gas reservoir can be determined by (39). This results in
(40)πre2hϕ=0.905qscThμCttpp(tΔm′)pp,
where (tΔm′)pp represents the actual value on the log-log plot of Δm′ versus t during the pseudosteady flow period.


## 3.4.1. Pure Wellbore Storage Period

The characteristic of pure wellbore storage period of horizontal well is manifested as a straight line segment with a slope of 1 on the log-log plot ofmwD, mwD′ versus tD/CD. The expression of dimensionless bottomhole pseudopressure during this period is
(22)mwD=tDCD.Equation (22) can be converted to dimensional form, and then according to the time and pressure data during pure wellbore storage period, the method to determine the wellbore storage coefficient can be obtained. By plotting the log-log plot of Δm, Δm′ versus t, the wellbore storage coefficient can be determined by the straight line segment with a slope of 1 on the log-log plot.The following can be obtained from the definitions of dimensionless variables:(23)tDCD=3.6Kht/ϕμctrw2C/2πϕcthL2=7.2πKhhL2μrw2tC.According to (22), (23), and the definition of dimensionless pseudopressure, the wellbore storage coefficient can be obtained. This results in
(24)C=0.288qscTμL2rw2tΔm,
where t/Δm represents the actual value on the log-log plot of Δm, Δm′ versus t during the pure wellbore storage period.


## 3.4.2. Early Vertical Radial Flow Period

The Determination of Geometric Mean Permeability. The dimensionless bottomhole pseudopressure derivative curve is manifested as a horizontal straight line segment with the value of 1/(4LD) during the early vertical radial flow period. The expression of dimensionless bottomhole pseudopressure derivative during this period is
(25)mwD′=dmwDdln⁡(tD/CD)=14LD.According to the definitions of dimensionless variables, the dimensional form of (25) can be obtained. This results in
(26)78.489KhhqscT(tΔm′)er=14(L/h)Kv/Kh.The geometric mean permeability of gas reservoir can be determined by (26). This results in
(27)KhKv=3.185×10-3qscTL(tΔm′)er,
where (tΔm′)er represents the actual value on the log-log plot of Δm′ versus t during the early vertical radial flow period.The Determination of Skin Factor and Initial Reservoir Pressure. The expression of dimensionless bottomhole pseudopressure during the early vertical radial flow period is
(28)mwD=14LD[ln⁡2.25(tDCD)+ln⁡CDe2S].According to the definitions of dimensionless variables and (28), the skin factor can be obtained. This results in
(29)S=0.5[Δmer(tΔm′)er-ln⁡KhterϕμCtrw2-0.80907],
where Δmer and ter represent the pseudopressure difference and time corresponding to the (tΔm′)er, respectively.For pressure buildup analysis, whenΔt→∞, the ln⁡Δt/(Δt+tp)→0, (30) can be obtained through the use of the definitions of dimensionless variables and pseudopressure difference during the early vertical radial flow period:
(30)mi-mwftΔm′=ln⁡tpD+0.80907+2S.The initial reservoir pseudopressure can be determined by (30). This results in
(31)mi=mwf+(tΔm′)erb(ln⁡tpD+0.80907+2S),
where (tΔm′)erb represents the actual value on the pressure buildup log-log plot of Δm′ versus t during the early vertical radial flow period.


## 3.4.3. Early Linear Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a straight line segment with the slope of 0.5 during the early linear flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, the following can be obtained:(32)78.489KhhqscT(tΔm′)l=rwL3.6πKhtϕμCtrw2.The horizontal permeability of gas reservoir can be determined by (32). This results in
(33)Kh=4.28×10-2qscTLhϕμCt[t(tΔm′)]l,
where (tΔm′)l represents the actual value on the log-log plot of Δm′ versus t during the early linear flow period.Combining (27) with (33), the vertical permeability can be obtained. This results in
(34)Kv=7.44×10-2hϕμCt[(tΔm′)/t]l(tΔm′)er.


## 3.4.4. Late Horizontal Pseudoradial Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a horizontal straight line segment with the value of 0.5 during the late horizontal pseudoradial flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, (35) can be obtained:
(35)78.489KhhqscT(tΔm′)lr=0.5.The horizontal permeability of gas reservoir can be determined by (35). This results in
(36)Kh=6.37×10-3qscTh(tΔm′)lr,
where (tΔm′)lr represents the actual value on the log-log plot of Δm′ versus t during the late horizontal pseudoradial flow period.For pressure buildup analysis, whenΔt→∞, the ln⁡Δt/(Δt+tp)→0, (37) can be obtained through the use of the definitions of dimensionless variables and pseudopressure difference during the late horizontal pseudoradial flow period:
(37)mi-mwf(tΔm′)lrb=ln⁡rwDtpD+0.80907+2S.The initial reservoir pseudopressure can be determined by (37). This results in
(38)mi=mwf+(tΔm′)lrb(ln⁡rwDtpD+0.80907+2S),
where (tΔm′)lrb represents the actual value on the pressure buildup log-log plot of Δm′ versus t during the late horizontal pseudoradial flow period.


## 3.4.5. Pseudosteady Flow Period

The dimensionless bottomhole pseudopressure derivative curve is manifested as a straight line segment with the slope of 1 during the pseudosteady flow period. According to the expression of dimensionless bottomhole pseudopressure derivative during this period and the definitions of dimensionless variables, (39) can be obtained:
(39)78.489KhhqscT(tΔm′)pp=7.2πKhtppre2hϕμCt.The pore volume of gas reservoir can be determined by (39). This results in
(40)πre2hϕ=0.905qscThμCttpp(tΔm′)pp,
where (tΔm′)pp represents the actual value on the log-log plot of Δm′ versus t during the pseudosteady flow period.


## 4. Example Analysis

The Longping 1 well is a horizontal development well in JingBian gas field, the well total depth is 4672 m, the drilled formation name is Majiagou group, the mid-depth of reservoir is 3425.63 m, and the well completion system is screen completion. According to the deliverability test during 26–29 December, 2006, the calculated absolute open flow was 94.26 × 104 m3/d. The commissioning data of Longping 1 well was in 12 May, 2007, the initial formation pressure was 29.39 MPa, before production, and the surface tubing pressure and casing pressure were both 23.90 MPa. The production performance curves of Longping 1 well are shown in Figures 8 and 9, respectively.Figure 8
The gas flow rate and water flow rate curve of Longping 1 well.Figure 9
The tubing head pressure and casing head pressure curve of Longping 1 well.Longping 1 well has been conducted pressure buildup test during 14 August, 2007, and 23 October, 2007. The gas flow rate of Longping 1 well was 40 × 104 m3/d before the shut-in. The bottomhole pressure recovered from 22.38 MPa to 27.83 MPa during the pressure buildup test. Physical parameters of fluid and reservoir are shown in Table 1. The pressure buildup log-log plot of Longping 1 well is shown in Figure 10.Table 1
Physical parameters of fluid and reservoir.
Parameter
Value
Initial formation pressurepi (MPa)
29.39
Formation temperatureT (°C)
95.80
Formation thicknessh (m)
6.31
Porosityϕ (%)
7.77
Initial water saturationSwi (%)
13.60
Well radiusrw (m)
0.0797
Gas gravityγg
0.608
Gas deviation factorZ
0.9738
Gas viscosityμg (mPa·s)
0.0222Figure 10
The pressure buildup log-log plot of Longping 1 well.As seen from the contrast between Figure10 and well test analysis type curves of horizontal well, the pseudopressure behavior of Longping 1 well manifests four characteristic periods during the pressure buildup test: pure wellbore storage period (I), early vertical radial flow period (II), early linear flow period (III), and late horizontal pseudoradial flow period (IV).Using the above characteristic value method of well test analysis for horizontal gas well, well test analysis results of Longping 1 well are shown in Table2. The pressure buildup semilog plot and pressure history matching plot of Longping 1 well are shown in Figures 11 and 12, respectively.Table 2
Well test analysis results of Longping 1 well.
Parameter
Parameter values
Wellbore storage coefficientC (m3/MPa)
1.229
Horizontal permeabilityKh (mD)
7.742
Vertical permeabilityKv (mD)
0.039
Flow capacityKhh (mD·m)
48.857
Skin factorS
-
2.49
Effective horizontal section lengthL (m)
198.37
Reservoir pressurepR (MPa)
28.385Figure 11
The pressure buildup semilog plot of Longping 1 well.Figure 12
The pressure history matching plot of Longping 1 well.


## 5. Summary and Conclusions

The four main conclusions and summary of this study are as follows.(1)
On the basis of establishing the mathematical models of well test analysis for horizontal gas well and obtaining the solutions of the mathematical models, several type curves which can be used to identify flow regime have been plotted and the seepage characteristic of horizontal gas well has been analyzed.(2)
The expressions of dimensionless bottomhole pseudopressure and pseudopressure derivative during each characteristic period of horizontal gas well have been obtained; formulas have been developed to calculate gas reservoir fluid flow parameters.(3)
The example study verifies that the characteristic value method of well test analysis for horizontal gas well is reliable and practical.(4)
The characteristic value method of well test analysis, which has been included in the well test analysis software at present, has been widely used in vertical well. As long as the characteristic straight line segments which are manifested by pressure derivative curve appear, the reservoir fluid flow parameters can be calculated by the characteristic value method of well test analysis for vertical well. The proposed characteristic value method of well test analysis for horizontal gas well enriches and develops the well test analysis theory and method.



---

*Source: 472728-2014-09-25.xml*