Open access peer-reviewed chapter

Well Test Analysis for Hydraulically-Fractured Wells

By Freddy Humberto Escobar

Submitted: May 31st 2018Reviewed: August 19th 2018Published: November 5th 2018

DOI: 10.5772/intechopen.80996

Downloaded: 509

Abstract

This chapter focuses on the application of Tiab’s direct synthesis (TDS) technique for practical and accurate interpretation of pressure tests on vertical wells in conventional reservoirs, so bilinear, linear, and elliptical flow regimes can be used for fracture characterization. Most fractured well interpretation tests are conducted using nonlinear regression analysis if the pressure model is available. This method has some drawbacks associated with the nonuniqueness of the solution. Also, the conventional straight-line method requires one plot for each individual flow regime observed in the pressure tests, and the estimated parameters cannot be verified. Tiab’s direct synthesis (TDS) methodology, which uses specific lines and intersection points found on the pressure and pressure derivative plot, is used in some direct equations which are obtained from the solution of the diffusivity equation for a given flow regime. It has been proven to provide accurate results, and its power allows verification of most results which is not possible from any other technique. The methodology has been successfully explained and tested by its application in two examples, although there exists more than a hundred articles that provide many useful applications.

Keywords

  • bilinear flow
  • linear flow
  • elliptical flow
  • half-length fracture
  • fracture conductivity
  • hydraulic fracturing

1. Introduction

Throughout their history, well test analyses for fractured wells have received many contributions. For practical purposes, let us name the most important ones for this chapter. A good place to start is by mentioning the work in [1], which described the pressure behavior for infinite-conductivity and uniform-flux fractured wells, so people started conducting interpretation tests on such wells by using type-curve matching. Later, [2] introduced the concept of finite-conductivity fractures and established the onset value of dimensionless conductivity as 300. Values lower than that are considered finite-conductivity values, and those above 300 are classified as infinite conductivity. In [2], a fine semi-analytical solution was introduced for describing the well-pressure behavior in hydraulically fractured wells. This solution was then applied in [3] to provide a well interpretation method using type-curve matching. Since then, other mathematical solutions have been presented for finite-conductivity fractures. Among them, the work in [4] using fractal theory is worth mentioning.

The way of conducting well test interpretation was changed by the introduction of Tiab’s direct synthesis (TDS) technique by [5]. This revolutionary and modern technique focuses on the different flow regimes seen on the pressure derivative curve. Defined lines are drawn through each individual flow regime, and the intersection points found among them are read and used for reservoir characterization. Additionally, reading arbitrary points on the pressure and pressure derivative of each flow regime also serve for reservoir parameter determination. A great number of applications of the TDS technique are given in [6]. The second work [7], by the same author of [5], presented TDS technique for infinite-conductivity and uniform-flux fractures in vertical wells. In [7], the elliptical or biradial flow regime was introduced and characterized. This elliptical flow is also seen in horizontal wells and was characterized in [8, 9, 10]. Because of the similarity between the mathematical models of hydraulic fractures and horizontal wells, this concept was applied by [11] to determine the average reservoir pressure in formations drained by horizontal wells using the TDS technique. The infinite-conductivity model in [7] also included the late-time pseudosteady-state period as well as some equations involved in the drainage area (conventional analysis for this case was included in [12]). This may be disadvantageous for inexperienced users of TDS technique when interpreting pressure tests without reaching reservoir boundaries because the equations involved the use of the unknown reservoir drainage area, although it can be still applied by using the intersection points. To overcome this drawback, [13] presented a new mathematical model excluding the late-time pseudosteady-state period.

TDS technique for finite-conductivity fractured wells is given in [14], with practical field applications to demonstrate the usefulness of the technique. The fracture parameters can be readily obtained by using an arbitrary point on the flow regimes. TDS technique plays an important role when analyzing short pressure tests because a user can “make up” nonexisting flow regimes since, for instance, the radial flow horizontal line can be obtained from the reservoir permeability even though radial flow regime is absent. [15, 16] extended the works of finite- and infinite-conductivity fractures in naturally fractured reservoirs. The equations provided by these works can also be applied to either homogeneous or naturally fractured formations since they involve a dummy variable that takes the value of one for the homogeneous case or the value of the dimensionless storativity coefficient for the case of a naturally fractured formation.

TDS technique has also been extended to several scenarios related to hydraulically fractured wells. For instance, when a finite-conductivity fracture intersects with a fault, the pressure trace changes; then, the equations developed in [17] apply for this case. There are cases where a threshold pressure is required to start the flow. The work in [18] includes this concept in uniform-fractured vertical wells, and the work in [19] includes the concept for horizontal wells. Also, when the fractured face is damaged, a pseudolinear flow regime develops along the fracture. [1] included TDS technique to characterize such systems. [16] presented TDS technique for fractured wells in gas composite reservoirs. TDS technique can also be usefully applied to transient-rate analysis, as seen in [20]. Application of TDS technique to horizontally isolated fractured wells was presented and characterized in [21] and in conventional analysis in [22]. The works in [23, 24] use TDS technique for shale reservoirs. Other applications of TDS technique to these systems are given by [25] under transient-rate analysis and [26] for pressure-transient analysis conditions. Other important applications of TDS Technique to fractured wells are given by [29, 30].

This chapter is devoted to the application of TDS technique to hydraulically fractured wells in either homogeneous or naturally fractured formations. Without given detailed derivations, the expressions for characterizing the hydraulic fracture parameters are presented along with the way they should be used. Important relationships and practical exercises are included.

2. TDS basis

The pioneer publication on the TDS technique, [5], explains in detail the derivation of the equations. The Laplace space solution of the arithmetic pressure derivative for a homogeneous and infinite reservoir with skin and wellbore storage is also presented in [5] and given by

PD=4π20eu2tDuuCDJ0u1CDsu2J1u2+uCDY0u1CDsu2Y1u2du.E1

However, we know that the pressure derivative is a horizontal line during radial flow regime. The dimensionless pressure derivative during radial line is easier represented by

tDPD=0.5.E2

Then, to obtain practical equations, dimensionless parameters must be used. The dimensionless time, based upon half-fracture length and reservoir drainage area, is given below:

tDxf=0.000263ktϕμctxf2E3

and

tDA=0.000263ktϕμctA.E4

The dimensionless pressure and pressure derivative parameters for oil reservoirs are given by

PD=khΔP141.2qμBE5

and

tDPD=khtΔP141.2qμB.E6

Finally, the dimensionless fracture conductivity introduced in [3] is defined as

CfD=kfwfkxf.E7

It is observed from Eq. (5) that the two key parameters of a hydraulic fracture are the half-fracture length, xf, and the fracture conductivity, kf wf. The total length of the fracture is given by 2 xf.

The easiest application of TDS technique is given by replacing the dimensionless pressure derivative defined by Eqs. (6) and (2), to provide an expression to readily determine formation permeability:

k=70.6qμBhtΔPR,E8

where (t*∆P′)R is the pressure derivative value during radial flow regime. The equations for the TDS technique are derived in the same manner Eq. (8) was obtained.

3. Biradial flow regime

Biradial or elliptical flow normally results in a hydraulically fractured well when areal anisotropy is present. This is recognized on the pressure derivative versus time log-log plot by a straight line with a slope of 0.36. In hydraulic fractures, the flow from the formation to the fracture is described by parallel flow lines resulting in a linear flow geometry better known as linear flow regime and characterized by a slope of 1/2 on the pressure derivative versus time log-log plot.

Both linear flow and biradial/elliptical flow regimes are seen on the plot of dimensionless pressure and pressure derivative versus dimensionless time based on half-fracture length for a naturally fractured formation. New expressions for the elliptical flow regime introduced in [13] excluding reservoir drainage area are given by.

PD=259πtDxf26ξ0.36E9

and

tDPD=πtDxf26ξ0.36,E10

being ξ a dummy variable that defines either a homogeneous or naturally fractured formation. When ξ = 1, a homogeneous reservoir is considered. For the case of naturally fractured formations, ξ = ω, the dimensionless storativity coefficient.

Once dimensionless parameters given by Eqs. (3), (5), and (6) are replaced into Eqs. (9) and (10), respectively, and solve for the half-fracture length, which yields

xf=22.5632qBhΔPBR1.3889tBRξϕctμk1.778E11

and

xf=5.4595qBhtΔPBR1.3889tBRξϕctμk1.778.E12

TDS technique is based on drawing a straight line throughout a given flow regime; then, the user is expected to read the pressure, ΔPBR, and pressure derivative, (tP’)BR, at a given time, tBR. A better way to reduce noise effects consists of extrapolating the mentioned straight line (biradial for this case) to the time of 1 h and read the pressure derivative value, (tP’)BR1, at 1 h. For this case, the pressure and pressure derivative set in Eqs. (11) and (12) is changed to ΔPBR1 and (tP’)BR1, respectively.

When bilinear flow is unseen, fracture conductivity can be found with an expression presented in [27]

kfwf=3.31739kesrw1.92173xf.E13

[5] also provided an equation for the determination of the skin factor using an arbitrary point read during radial flow regime:

s=0.5ΔPRtΔPRlnktRϕμctrw2+7.43.E14

The pseudosteady-state regime governing the pressure derivative equation is given by

tDAPDP=2πtDAP.E15

[7] used the point of intersection, tRPi, of Eqs. (2) and (15) to derive an equation for the estimation of the drainage area:

A=ktRPi301.77ϕμct.E16

The derivation of Eq. (16) follows a similar idea as that presented later in Section 4 for the use of the points of intersection.

4. Bilinear and linear flow regimes

Bilinear flow regime takes place when two linear flows, normal one flowing into the other, take place simultaneously. This situation occurs in low conductivity fractures where linear flow along the fracture and linear flow from the formation to the fracture are observed. Bilinear flow is recognized in the pressure derivative curve by a slope of 0.25. However, this is not shown in Figure 1 since bilinear flow is absent. The governing expressions for early bilinear and linear flow regimes for vertical fractures in naturally fractured systems were, respectively, presented in [16]

PD=2.45CfDtDxfξ1/4,E17
tDPD=0.6125CfDtDxfξ1/4,E18
PD=πtDxfξ1/2,E19

and

tDPD=12πtDxfξ1/2.E20

Figure 1.

Dimensionless pressure and pressure derivative behavior for an infinite-conductivity fractured vertical well in a naturally fractured bounded reservoir, λ = 1 × 10−9 and ω = 0.1 (taken from [13]).

Linear flow regime can be used to find the half-fracture length, and bilinear flow regime allows finding the fracture conductivity. Once the dimensionless quantities of Eqs. (1) and (3)(5) are replaced in Eqs. (16)(19), the fracture conductivity is solved for then

kfwf=1947.46ξϕμctkqμBhΔPBL12,E21
kfwf=121.74ξϕμctkqμBhtΔPBL12,E22

Once the fracture conductivity is found, Eq. (7) applies to find the dimensionless fracture conductivity if reservoir permeability and the half-fracture length are known. When bilinear flow is absent, the fracture conductivity may be found from Eq. (13), or the dimensionless fracture conductivity can be read from Figure 2:

xf=4.064qBhΔPL1μξϕctkE23

and

xf=2.032qBhtΔPL1μξϕctk.E24

Figure 2.

Effect of skin factor on fracture conductivity (taken from [28]).

5. Points of intersection

If bilinear flow also takes place, then the point of intersection between the pressure derivatives of the bilinear and biradial flow lines, tBLBRi, given by Eqs. (10) and (18), respectively, allows the development of an equation to find the half-fracture as follows:

πtDxf26ξ0.36=0.6125CfDtDxfξ1/4.E25

Simplifying,

tDxfξ0.11=0.2862CfD.E26

Replacing the dimensionless quantities, Eqs. (3) and (7) in Eq. (26) lead to

0.000263ktϕμctxf2ξ0.11=0.2862kxfkfwf.E27

Solving for the half-fracture from Eq. (27), we readily obtain

kfwf=10.5422ξϕμctk3.5454xf6.5454tBRBLi0.22.E28

By the same token, the intercept of Eq. (20) with Eq. (18), tLBRi, provides another expression to find the half-fracture length:

xf=ktLBRi39.044ωϕμct.E29

Bilinear flow regime is absent in the plot of Figure 1. Linear, biradial, and radial flow regimes along with the late pseudosteady-state period are seen. The interception points formed by the possible combinations of such periods can be represented schematically in this plot.

Another way to find the half-fracture length comes from the intersection of Eqs. (2) and (10), tRBRi, and Eq. (10) with Eq. (15), so that

xf=14584.16ktRBRiξϕμctE30

and

xf=41.0554A1.3889ξϕμctktBRPi0.8889.E31

The intercept point resulting between linear flow and bilinear flow lines given by the governing pressure derivative solutions, Eqs. (18) and (19), can be used to find either half-fracture length or permeability:

k=kfwfxf2216tBLLi13910ξϕμct.E32

tBLRi is the intersection of the bilinear pressure derivative line given by Eq. (18) with the radial flow regime line (Eq. (2)). This intersection point serves as the estimation of either permeability or fracture conductivity:

tBLRi=1677ξϕμctk3kfwf2.E33

6. Other estimations

The expressions for determination of the naturally fractured reservoir parameters cannot be included in this chapter for space reasons. However, they can be found in [15, 16], which also used intersection points and maximum and minimum data read from the pressure and pressure derivative curve.

Radial flow regime may be absent in short tests run in fractured wells with the sole purpose of determining fractured parameters. For these cases, the skin factor can be estimated from any of the two empirical correlations presented by [27]

s=lnrw1.92173xf3.31739kfwfE34

and

s=lnrwxf+1.650.32u+0.11u21+0.18u+0.064u2+0.005u3,E35

where

u=lnCfD.E36

Additionally, fracture conductivity can be read from the plot given in Figure 2.

Finally, space reasons prevent including TDS technique for fractured wells in unconventional shale formations. The reader is referred to [23, 24, 25, 26].

7. Examples

7.1 Field example

[14] presented a field example of a fractured well test. Pressure and pressure derivative data are given in Table 1 and Figure 3. Other relevant data are provided below:

q=101STB/Dϕ=0.08μ=0.45cp
ct=17.7×106psia1B=1.507bbl/STBh=42ft
rw=0.28fttp=2000hPi=2200psia
ξ=1
t, hΔP, psiatP’, psiat, hΔP, psiatP’, psia
0.2310226.315390117
0.391153020423112
0.613035.825446120
114540.830471141
1.818357.235493136.5
2.41956740510132
3.826083.345526135
4.126569.250540150
4.9628096.955556137.5
6.2308102.360565144
8.5320103.365580121.1
1034514971583

Table 1.

Pressure data for field example (taken from [14]).

Figure 3.

Pressure and pressure derivative against time log–log plot for field example (taken from [14]).

Using a commercial well test software, the following parameters were estimated by nonlinear regression analysis:

k=0.8md
xf=82.2ft
kfwf=300mdcp

The objective is to compute the hydraulic fracture parameters using the TDS technique and compare results obtained from the regression analysis.

7.1.1 Solution

7.1.1.1 Step 1: Obtain the characteristic points

Once the pressure and pressure derivative versus time log-log plot is built and reported in Figure 3, the characteristic points are read from such plot as follows:

tR=30hΔPR=471psiatΔPR=150psia
tΔPBL1=160psiaΔPBL1=40psiaΔPL1=120psia
tLRi=8.2htBLRi=195h

7.1.1.2 Step 2: Estimate permeability and skin factor

Permeability and skin factor are found in Eqs. (8) and (14) to be 0.76 md and −4.68, respectively.

7.1.1.3 Step 3: Estimate fracture conductivity

Fracture conductivity is estimated using Eqs. (21) and (22):

kfwf=121.740.080.4517.7×1060.761010.451.50742402=290.77mdft
kfwf=1947.460.080.4517.7×1060.761010.451.507421602=290.7mdft.

From Figure 3, tBLRi = 200 hr. A very close value is obtained from Eq. (33):

tBLRi=16770.080.4517.7×1060.763290.72=205.71hr,

which indicates that the calculation of the fracture conductivity is accurate. Notice that instead of estimating tBLRi the fracture conductivity can be found instead to obtain another value of fracture conductivity; then, Eq. (33) can also be expressed as

kfwf=k3tBLRi1677ξϕμct=0.763205167710.080.4517.7×106=290.2mdft.

7.1.1.4 Step 4: Half-fractured length and dimensionless fracture conductivity estimation

Find half-fracture length with Eqs. (23) and (24):

xf=4.0641011.507421200.450.0817.7×1060.76=79ft,
xf=ktLri1207ξϕμct=0.761012070.080.7617.7×106=76.5ft.

Solve for half-fracture length from Eq. (13) and find this:

xf=1.92173esrw3.31739kwfkf=1.92173e4.68440.283.317390.76290.7=79ft.

Find the dimensionless fracture conductivity using Eq. (5):

CfD=wfkfxfk=290.7790.76=4.8.

The above value confirms that the fracture has finite conductivity.

7.2 Synthetic example

[13] presented a synthetic example of a pressure test run in a bounded homogeneous reservoir with the information given below:

Bo=1.25bbl/STBq=300STB/D
h=30ftμ=5cp
rw=0.3ftct=3×106psi1
Pi=4000psiϕ=10%
k=33.334mdxf=200ft
A=592Acres

Estimate the half-fracture length by the TDS technique, and compare the answer with the value used for generating the test.

7.2.1 Solution

7.2.1.1 Step 1: Obtain the characteristic points.

A pressure and pressure derivative versus time log–log plot is presented in Figure 4, from which the following characteristic points are read:

tBR=1.01htΔPBR=64.63psitBRPi=3300h

Figure 4.

Pressure and pressure derivative vs. time for synthetic example (taken from [13]).

7.2.1.2 Step 2: Half-fractured length estimation.

The half-fracture length is estimated with Eq. (12) and confirmed with Eq. (31), as follows:

xf=5.45953001.253064.631.38891.0110.13×106533.3341.778=199ft,
xf=41.0554592×435601.388910.153×10633.33433000.8889=201.6ft.

8. Comments on the results

The two given examples show three aspects of the TDS technique: (1) practical use, (2) accuracy, and (3) self-confirmation.

As shown in the exercises, the process includes defining flow regimes, drawing a few lines, and finally computing the necessary parameters. Contrary to the conventional straight-line method, which requires a plot for each flow regime, TDS technique uses only the pressure and pressure derivative versus time log-log plot. Computations are straight forward.

Table 2 summarizes the main parameters obtained in the two worked examples. The results show a good agreement between the calculated results by TDS technique and the results obtained from commercial software packages, for the field case. The results of the half-fracture length for the synthetic case using TDS technique are even better compared to the input value use to simulate the test. This demonstrates that TDS technique is an accurate methodology which has been also presented in many publications, not only in the list reference but also in others not mentioned here.

Field example
Obtained from
ParameterCommercial softwareEq. (21)Eq. (22)Eq. (33)Eq. (23)Eq. (24)Eq. (13)
xf, ft82.27976.579
kf wf, md-ft300290.77290.7290.2
Synthetic example
Obtained from
ParameterCommercial softwareEq. (12)Eq. (31)
xf, ft200199201.6

Table 2.

Summary of results.

The last aspect dealt with is self-confirmation. In the field example, three values of half-fracture length and three values of fractured conductivity were found, and for the synthetic example, two values of half-fracture length were estimated from different equations. All the estimations match with the reference values.

9. Conclusion

It has been shown that TDS technique is a powerful, practical, and accurate tool for well test interpretation because manipulations are easy to do and parameters can be confirmed from different sources from the same pressure test. Compared to reference values, the worked examples provided accurate results of both half-fracture length and hydraulic fracture conductivity. Besides being accurate, TDS technique has the great advantage of being able to estimate a given parameter, such as half-fracture length or fracture conductivity, from more than one source or equation. This provides a means of verifying that the estimated parameter is in a good range.

Acknowledgments

The author wishes to express his gratitude to Universidad Surcolombiana for providing him the time to write this chapter.

Nomenclature

ADraining area (ft2)
BOil volume factor (rb/STB)
CfDDimensionless fracture conductivity
ctCompressibility (1/psi)
hFormation thickness (ft)
kFormation permeability (md)
kfwfFracture conductivity (md-ft)
PPressure (psi)
PwfWell-flowing pressure (psi)
qOil flow rate (STB/D)
qgGas flow rate (MSCF/D)
rwWellbore radius (ft)
xfHalf-fracture length (ft)
sSkin factor
tTest time (h)
tpProduction time (h)
t*∆P′Pressure derivative (psi)
tD*PD’Dimensionless pressure derivative
Greek symbols
Change
ϕPorosity (fraction)
λInterporosity flow parameter
μViscosity (cp)
ξVariable to identify homogeneous (ξ = 1) or heterogeneous (ξ = ω) reservoirs
ωDimensionless storativity coefficient
Suffixes
BLBilinear
BL1Bilinear at 1 h
BLLBilinear-linear intersection
BRBirradial
BR1Birradial at 1 h
BRBLiBirradial-bilinear intersection
BRPiBirradial-pseudosteady intersection
DDimensionless
DADimensionless based on area
DxfDimensionless based on half-fractured length
DLBRiDual linear-birradial intersection
LBRiLinear-birradial intersection
RRadial
RBRiRadial-birradial intersection
RPiIntersect of radial-pseudosteady-state lines
wWell
tTime
PPseudosteady state

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Freddy Humberto Escobar (November 5th 2018). Well Test Analysis for Hydraulically-Fractured Wells, Exploitation of Unconventional Oil and Gas Resources - Hydraulic Fracturing and Other Recovery and Assessment Techniques, Kenneth Imo-Imo Eshiet, IntechOpen, DOI: 10.5772/intechopen.80996. Available from:

chapter statistics

509total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Surface Drilling Data for Constrained Hydraulic Fracturing and Fast Reservoir Simulation of Unconventional Wells

By Ahmed Ouenes, Mohit Paryani, Yamina Aimene, Chad Hammerquist and Aissa Bachir

Related Book

First chapter

The Influence of Engine Fuel Manufacturing Processes on Their Performance Properties in Operating Conditions

By Krzysztof Biernat

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us