You must log in to edit PetroWiki. Help with editing

Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information


Fluid flow in hydraulically fractured wells

PetroWiki
Revision as of 20:04, 4 June 2015 by Denise Watts (Denisewatts) (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Many wells, particularly gas wells in low-permeability formations require hydraulic fracturing to be commercially viable. Interpretation of pressure-transient data in hydraulically fractured wells is important for evaluating the success of fracture treatments and predicting the future performance of fractured wells. This page includes graphical techniques for analyzing post-fracture pressure transient tests after identifying several flow patterns that are characteristic of hydraulically fractured wells. Often, identification of specific flow patterns can aid in well test analysis.

Flow patterns in hydraulically fractured wells

Five distinct flow patterns (Fig. 1) occur in the fracture and formation around a hydraulically fractured well.[1] Successive flow patterns, which often are separated by transition periods, include fracture linear, bilinear, formation linear, elliptical, and pseudoradial flow. Fracture linear flow (Fig. 1a) is very short-lived and may be masked by wellbore-storage effects. During this flow period, most of the fluid entering the wellbore comes from fluid expansion in the fracture, and the flow pattern is essentially linear.

Because of its extremely short duration, the fracture linear flow period often is of no practical use in well test analysis. The duration of the fracture linear flow period is estimated by[1]

RTENOTITLE....................(1)

where tLfD is dimensionless time in terms of fracture half-length,

RTENOTITLE....................(2)

The dimensionless fracture conductivity, Cr, is

RTENOTITLE....................(3)

and ηfD is dimensionless hydraulic diffusivity defined by

RTENOTITLE....................(4)

Bilinear flow (Fig. 1b) evolves only in finite-conductivity fractures as fluid in the surrounding formation flows linearly into the fracture and before fracture tip effects begin to influence well behavior. Fractures are considered to be finite conductivity when Cr < 100. Most of the fluid entering the wellbore during this flow period comes from the formation. During the bilinear flow period, BHP, pwf, is a linear function of t1/4 on Cartesian coordinates.

A log-log plot of (pipwf) as a function of time exhibits a slope of 1/4 unless the fracture is damaged. The pressure derivative also has a slope of 1/4 during this same time period. The duration of bilinear flow depends on dimensionless fracture conductivity and is given by Eqs. 5a through 5c[1] for a range of dimensionless times and fracture conductivities:

RTENOTITLE....................(5a)

RTENOTITLE....................(5b)

and RTENOTITLE....................(5c)

Formation linear flow (Fig. 1c) occurs only in high-conductivity (Cr ≥ 100) fractures. This period continues to a dimensionless time of tLfD ≅ 0.016. The transition from fracture linear flow to formation linear flow is complete by a time of tLfD = 10–4 . On Cartesian coordinates, p wf is a linear function of t 1/2 , and a log-log plot of (pipwf) has a slope of 1/2 unless the fracture is damaged. The pressure derivative plot exhibits a slope of 1/2. Elliptical flow (Fig. 1d) is a transitional flow period that occurs between a linear or near-linear flow pattern at early times and a radial or near—radial flow pattern at late times.

Pseudoradial flow (Fig. 1e) occurs with fractures of all conductivities. After a sufficiently long flow period, the fracture appears to the reservoir as an expanded wellbore (consistent with the effective wellbore radius concept suggested by Prats et al.[2]). At this time, the drainage pattern can be considered as a circle for practical purposes. (The larger the fracture conductivity, the later the development of an essentially radial drainage pattern.) If the fracture length is large relative to the drainage area, then boundary effects distort or entirely mask the pseudoradial flow regime. Pseudoradial flow begins at tLfD ≅ 3 for high-conductivity fractures (Cr ≥ 100) and at slightly smaller values of tLfD for lower values of Cr.

These flow patterns also appear in pressure-buildup tests and occur at approximately the same dimensionless times as in flow tests. The physical interpretation is that the pressure has built up to an essentially uniform value throughout a particular region at a given time during a buildup test. For example, at a given time during bilinear or formation linear flow, pressure has built up to a uniform level throughout an approximately rectangular region around the fracture. At a later time during elliptical flow, pressure has built up to a uniform level throughout an approximately elliptical region centered at the wellbore. At a given time during pseudoradial flow, pressure has built up to a uniform level throughout an approximately circular region centered at the wellbore. The area of the region and the pressure level within that area increase with increasing shut-in time. Example 1 illustrates how to estimate the duration of flow periods for hydraulically fractured wells.

Example 1: Estimating duration of flow periods in a hydraulically fractured well

For each case, estimate the end of the linear flow period and the time at which pseudoradial flow period begins. Assume that pseudoradial flow begins when tLfD = 3. Table 1 gives the data for each case.

Solution. The end of the linear flow regime occurs at a dimensionless time of tLfD ≅ 0.016 or, using Eq. 2,

RTENOTITLE

Similarly, the time to reach pseudoradial flow is tLfD ≅ 3, or

RTENOTITLE

Table 1 summarizes the results.

Flow geometry and depth of investigation of a vertically fractured well

Fluid flow in a vertically fractured well has been described using elliptical geometry. [3] The equation for an ellipse with its major axis along the x-axis and minor axis along the y-axis is

RTENOTITLE....................(6)

where the endpoints of the major and minor axes are (±af, 0) and (0, ±bf), respectively. The foci of the ellipse are ±cf where cf2 = af2bf2. In terms of a well with a single vertical fracture with two wings of equal length, Lf, the relation becomes Lf2 = af2bf2, where Lf is the focal length of the ellipse. Fig. 8.76 shows the elliptical geometry of a vertically fractured well.

Hale and Evers[3] defined a depth of investigation for a vertically fractured well. Their definition is based on a definition of dimensionless time at a distance bf, the length of the minor axis:

RTENOTITLE....................(7)

Solving for the length of the minor axis,

RTENOTITLE....................(8)

Assuming that pseudosteady-state flow exists out to distance, bf, at dimensionless time tbD = 1/π as in linear systems, Eq. 8 becomes

RTENOTITLE....................(9)

which represents the depth of investigation in a direction perpendicular to the fracture at time, t, for a vertically fractured well. In gas wells, the terms μ and ct should be RTENOTITLE and RTENOTITLE, evaluated at average drainage-area pressure, RTENOTITLE.

The elliptical pattern of the propagating pressure transient can be fully described in terms of the lengths of the major axis, af, the minor axis, bf, and the focus, Lf. Using the estimate of bf from Eq. 9 and an estimate of Lf obtained by one of the methods described in sections that follow, the length of the major axis can be estimated from

RTENOTITLE....................(10)

Given values of af and bf, the depth of investigation at a particular time, t, in any direction from the fracture can be calculated using Eq. 6. Furthermore, the area, A, enclosed by the ellipse at time, t (the area of the reservoir sampled by the pressure transient), is given by

RTENOTITLE....................(11)

The coefficient 0.0002878 in Eq. 9 is strictly correct only for highly conductive fractures (Cr ≥ 100). As Cr becomes smaller, the ratio af/bf also becomes smaller. The lower bound of af/bf is 1 (a circle) as Cr approaches 0.

Fracture damage

Two major types of fracture damage are frequent: choked fracture damage and fracture-face damage. The choked-fracture damage means that the fracture has a reduced permeability in the immediate vicinity of the wellbore (Fig. 3). In this case, kf is used for the permeability in the propped portion of the fracture farther along the wellbore, and kfs for reduced permeability near the wellbore, out to a length, Ls, in the fracture.

The choked-fracture skin factor, sf, is[4]

RTENOTITLE....................(12)

Fracture face damage in a hydraulically fractured well (Fig. 4) is a permeability reduction around the edges of the fracture, usually caused by invasion of the fracture fluid into the formation or an adverse reaction with the fracturing fluid. The equation for fracture face skin is[4]

RTENOTITLE....................(13)

Specialized methods for post-fracture well-test analysis

Generally, the objectives of post-fracture pressure-transient test analysis are to assess the success of the fracture treatment and to estimate:

  • The fracture half-length
  • Fracture conductivity
  • Formation permeability

Three specialized methods of analyzing these post-fracture transient tests are included here:

  • Pseudoradial flow
  • Bilinear flow
  • Linear flow

Bilinear flow method

The bilinear flow method[5] applies to test data obtained during the bilinear flow regime in wells with finite-conductivity vertical fractures. Bilinear flow is indicated by a quarter-slope line on a log-log graph of pressure derivative vs. t or Δte.

During bilinear flow,

RTENOTITLE....................(14)

and RTENOTITLE....................(15)

The following procedure is recommended for analyzing test data obtained in the bilinear flow regime (that is, data in the time range with quarter slope on the diagnostic plot). In Step 1, note the use of "bilinear equivalent time," ΔtBe. Radial equivalent time is rigorously correct as a plotting function only for infinite-acting radial flow.

  1. For a constant-rate flow test, plot pwf vs. t1/4 on Cartesian coordinates. For a buildup test, plot pws vs. ΔtBe1/4, where
    RTENOTITLE....................(16)
  2. Determine the slope, mB, of the straight line region of the plot.
  3. Determine the pressure extrapolated to time zero, po, and the fracture skin, sf, from
    RTENOTITLE....................(17)
    for drawdown and buildup tests, respectively.
  4. From independent knowledge of k (for example, from a prefracture well test), estimate the fracture conductivity, wfkf, using mB and the relationship
    RTENOTITLE....................(18)
    where RTENOTITLE and RTENOTITLE, evaluated at RTENOTITLE, are used for a gas well test.

Fig. 5 is an example of bilinear flow analysis. The bilinear flow analysis method has the following important limitations.

  • No estimate of fracture half-length, Lf.
  • In wells with low-conductivity fractures, wellbore storage frequently distorts early test data for a sufficient length of time so that the quarter-slope line characteristic of bilinear flow may not appear on a log-log plot of test data.
  • An independent estimate of k is required. This suggests that prefracture well tests should be conducted before fracturing the well, thus obtaining independent estimates of formation properties.

Linear flow method

The linear flow method[5] applies to test data obtained during formation linear flow in wells with high-conductivity fractures (Cr ≥ 100). After wellbore storage effects have ended, formation linear flow occurs up to a dimensionless time of tLfD = 0.016, which means that a log-log plot of pressure derivative against time will have a slope of one-half. The plot of pressure change vs. time, however, will have a half-slope only if the fracture skin is zero. The pressure and pressure derivative are

RTENOTITLE....................(19)

and RTENOTITLE....................(20)

so that

RTENOTITLE....................(21)

which indicates that a log-log plot of the derivative against time will have a slope of one-half. Radial equivalent time applies rigorously only for radial flow in an infinite-acting reservoir. When linear flow is the flow pattern occurring at both times (tp + Δt) and Δt, a more useful equivalent time function is the linear equivalent time, ΔteL.

RTENOTITLE....................(22)

Test conditions in which linear flow occurs at both (tp + Δt) and Δt are rare, and, consequently, Eq. 22 is not necessarily rigorously correct for well-test analysis. Fortunately, when tp >> Δtmax, ΔteL ≈ Δt. Fig. 6 is an example of a plot used in linear flow analysis.

The linear flow analysis method also has limitations.

  • The method applies only for fractures with high conductivities. Strictly speaking, linear flow occurs for the condition of uniform flux into a fracture (same flow rate from the formation per unit cross-sectional area of the fracture at all points along the fracture) rather than for infinite fracture conductivity. Therefore, only very early test data (tLf D ≤ 0.016) exhibit linear flow in a high-conductivity fracture.
  • Some or all of these early data may be distorted by wellbore storage, further limiting the amount of linear-flow data available for analysis.
  • Estimating fracture half-length requires an independent estimate of permeability, k, which suggests the need for a prefracture well test.

Pseudoradial flow method

The pseudoradial flow method applies when a short, highly conductive fracture is created in a high-permeability formation, so that pseudoradial flow develops in a short time. The time required to achieve pseudoradial flow for an infinitely conductive fracture (Cr ≥ 100) in either a flow test or a pressure buildup test is estimated by

RTENOTITLE....................(23)

The beginning of pseudoradial flow is characterized by the flattening of the pressure derivative on a log-log plot and by the start of a straight line on a semilog plot. Hence, when the pseudoradial flow regime is reached, conventional semilog analysis can be used to calculate permeability and skin factor. For a highly conductive fracture, skin factor is related to fracture half-length by[2]

RTENOTITLE....................(24)

Fig. 7 shows an example.

A recommended procedure for analyzing test data from the pseudoradial flow regime is as follows.

  • For a drawdown test, plot pwf vs. log t. For a buildup test, plot pws vs. the Horner time ratio (HTR).
  • Determine the position and slope, m, of the semilog straight line and the intercept, p1hr on the line.
  • Using m, calculate values of k and s (or s′ for a gas well).
  • Calculate the fracture half-length, Lf, using Eq. 24.

The pseudoradial flow method has the following limitations that seldom make it applicable in practice. [5]

  • The conditions that are most favorable for the occurrence of pseudoradial flow are short, highly conductive fractures in high-permeability formations. These formations, however, are rarely fractured. The most common application of hydraulic fractures—wells with long fractures in low-permeability formations—require impractically long test times to reach pseudoradial flow.
  • For gas wells, the apparent skin factor, s′, calculated from test data is often affected by non-Darcy flow.
  • The pseudoradial method applies only to highly conductive (Cr ≥ 100) fractures. For lower conductivity fractures, fracture lengths calculated using the skin factor (Eq. 24) will be too low.

Using type curves for hydraulically fractured wells

Type curves are the most common method of analyzing hydraulically fractured wells. The independent variable for most type curves for analyzing hydraulically fractured wells is the dimensionless time based on hydraulic fracture half-length, tLf D. The dependent variable is usually the dimensionless pressure, pD.

For type curves used for manual type-curve matching, most vary only one parameter. The Cinco type curve[1] is obtained for zero CLf D and sf ; the only parameter is dimensionless fracture conductivity, Cr or FcD (where FcD = πcr). The choked-fracture skin is analyzed by assuming CLf D and infinite Cr with single parameter sf. The wellbore-storage type curve[6] sets sf to 0 and Cr (FcD) to infinity and varies the coefficient CLf D.

When using type curves in commercial software, the computer can set any two of the three parameters to fixed values (other than their limiting values) and vary the third parameter to obtain the matching stems.

Procedures for analyzing fractured wells with type curves

The following steps outline the procedure for analyzing fractured wells with type curves.

  • Graph field data pressure change and pressure derivatives.
  • Match field data to the appropriate type curve.
  • Find the match point and matching stem.
  • Calculate the formation permeability from the pressure match point.
  • Calculate Lf from the time match point.
  • Interpret the matching stem value appropriate for a given type curve. For one type curve, this can be wfkf, which will provide an estimate of fracture conductivity. For another, it can be sf, the choked-fracture skin, or, for a third, it can be C, the wellbore-storage coefficient.

To interpret the match points for a test with unknown permeability, use Eqs. 25 and 26. The formation permeability, k, is determined from the pressure match point; that is, the relationship between the pressure derivative and pressure change found at a match point given by

RTENOTITLE....................(25)

From the time match point, calculate the fracture half-length:

RTENOTITLE....................(26)

Matching can be ambiguous for hydraulically fractured wells; the data can appear to match equally well in several different positions. The ambiguity can be reduced or eliminated if a prefracture permeability is determined, and the post-fracture test data forced to match the permeability.

Type curves used for analysis in fractured wells

The Cinco type curve (Fig. 8), [1] assumes that CLf D = 0 and sf = 0. The type-curve stems on this curve are obtained by varying values of Cr or FcD. With the Cinco type curve, the fracture conductivity, wfkf, can be determined from the matching parameter:

RTENOTITLE....................(27)

Choked-fracture type curve

Fig. 9 shows the choked-fracture type curve. [4] The choked-fracture type curve is generated with wellbore-storage coefficient, CLf D, of zero and infinite fracture conductivity, Cr. On this type curve, the stems represent different values of the fracture skin, sf. The fracture skin, sf, can be used to find the additional pressure drop from

RTENOTITLE....................(28)

Wellbore-storage type curve

The wellbore-storage type curve (Fig. 10) takes into account the possibility of wellbore storage. The wellbore-storage type assumes sf = 0 and Cr = ∞. To interpret a best-fitting stem for this type curve, use the following:

RTENOTITLE....................(29)

Limitations of type-curve analysis

Although it is the most common methodology for analyzing hydraulically fractured well, type-curve analysis still has some limitations.

First, type-curves for analysis of hydraulically fractured wells are usually based on solutions for constant-rate drawdown tests. For buildup tests, shut-in time itself may possibly be used as a plotting function in those cases in which producing time is much greater than the shut-in time. Equivalent time can be used in some cases, but equivalent time has different definitions depending on the flow regime: radial, linear, and bilinear flow. Another possibility is to use a "superposition" type curve, which depends on the specific durations of flow and buildup periods. Superposition type curves can be readily generated with computer software.

Another problem with type curves is that they may ignore important behavior. The type curve that takes into account wellbore storage does not consider a variable wellbore storage coefficient. This can be caused by phase redistribution in the wellbore, for example. The widely available type curves that have been discussed do not include boundary effects. With gas wells, the probability of non-Darcy flow is high, but available type curves don’t take this into account.

An independent estimate of permeability may also be needed. A number of different type curves or a variety of stems on a given type curve may seem to match test data equally well. To remove this ambiguity, the best solution is to have an independent estimate of permeability.

Nomenclature

a = RTENOTITLE, stabilized deliverability coefficient, psia2-cp/MMscf-D
a = total length of reservoir perpendicular to wellbore, ft
ah = length of reservoir perpendicular to horizontal well, ft
af = RTENOTITLE, depth of investigation along major axis in fractured well, ft
at = RTENOTITLE, transient deliverability coefficient, psia2-cp/MMscf-D
A = drainage area, sq ft
A = πafbf , area of investigation in fractured well, ft2
Af = cross-sectional area perpendicular to flow, sq ft
Awb = wellbore area, sq ft
b = RTENOTITLE (gas flow equation)
bf = RTENOTITLE, depth of investigation of along minor axis in fractured well, ft
B = formation volume factor, res vol/surface vol
cf = formation compressibility, psi–1
ct = Soco + Swcw + Sgcg + cf = total compressibility, psi–1
CLfD = 0.8936 RTENOTITLE, dimensionless wellbore storage coefficient in fractured well
Cr = wfkf/πkLf, fracture conductivity, dimensionless
D = non-Darcy flow constant, D/Mscf
Ef = flow efficiency, dimensionless
hD = (h/rw)(kh/kv)1/2, dimensionless
hf = fracture height, ft
hp = perforated interval thickness, ft
hpD = hp/ht
ht = total formation thickness, ft
k = matrix permeability, md
RTENOTITLE = average permeability, md
kf = permeability of the proppant in the fracture, md
L = distance from well to no-flow boundary, ft
Lf = fracture half length, ft
Ls = length of damaged zone in fracture, ft
pf = formation pressure, psi
pi = original reservoir pressure, psi
ps = stabilized shut-in BHP measured just before start of a deliverability test, psia
pt = surface pressure in tubing, psi
pwf = flowing BHP, psi
pws = shut-in BHP, psi
pD = 0.00708 kh(pip)/qBμ, dimensionless pressure as defined for constant-rate production
q = flow rate at surface, STB/D
qw = water flow rate, STB/D
r = distance from the center of wellbore, ft
rp = radius of perforation tunnel, ft
rs = outer radius of the altered zone, ft
rsp = radius of source or inner boundary of spherical flow pattern, ft
rw = wellbore radius, ft
rD = r/rw, dimensionless radius
sf = skin of hydraulically fractured well, dimensionless
tbD = dimensionless time in linear flow, hours
tLfD = 0.0002637 kt/ϕμetLf2, dimensionless time for fractured wells
tp = constant-rate production period, t, hours
T = reservoir temperature, °R
u = dummy variable
V = volume, bbl
Vf = fraction of bulk volume occupied by fractures
wf = fracture width, ft
wkf = fracture conductivity, md-ft
ws = width of damaged zone around fracture face, ft
Δp = pressure change since start of transient test, psi
Δt = time elapsed since start of test, hours
Δta = RTENOTITLE, normalized or adjusted pseudotime, hours
Δtap = RTENOTITLE, pseudotime, hr-psia/cp
ΔtBe = bilinear equivalent time, hours
Δte = radial equivalent time, hours
ΔtLe = linear equivalent time, hours
η = 0.0002637 k/ϕμct, hydraulic diffusivity, ft2/hr
ηfD = hydraulic diffusivity, dimensionless
λ = interporosity flow coefficient
λt = RTENOTITLE, total mobility, md/cp
α = exponent in deliverability equation
α = parameter characteristic of system geometry in dual-porosity system
μ = viscosity, cp
μw = water viscosity, cp
RTENOTITLE = gas viscosity evaluated at average pressure, cp
ϕf = fraction of fracture volume occupied by pore space, ≅ 1
ϕ = porosity, dimensionless

References

  1. 1.0 1.1 1.2 1.3 1.4 1.5 Cinco-Ley, H. and Samaniego-V., F. 1981. Transient Pressure Analysis for Fractured Wells. J Pet Technol 33 (9): 1749-1766. SPE-7490-PA. http://dx.doi.org/10.2118/7490-PA
  2. 2.0 2.1 Prats, M., Hazebroek, P., and Strickler, W.R. 1962. Effect of Vertical Fractures on Reservoir Behavior--Compressible-Fluid Case. SPE J. 2 (2): 87-94. http://dx.doi.org/10.2118/98-PA
  3. 3.0 3.1 Hale, B.W. and Evers, J.F. 1981. Elliptical Flow Equations for Vertically Fractured Gas Wells. J Pet Technol 33 (12): 2489–2497. SPE-8943-PA. http://dx.doi.org/10.2118/8943-PA
  4. 4.0 4.1 4.2 Cinco-Ley, H. and Samaniego-V., F. 1981. Transient Pressure Analysis: Finite Conductivity Fracture Versus Damaged Fracture Case. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 4–7 October. SPE-10179-MS. http://dx.doi.org/10.2118/10179-MS
  5. 5.0 5.1 5.2 Lee, W.J. 1989. Postfracture Formation Evaluation. In Recent Advances in Hydraulic Fracturing, J.L. Gidley, S.A. Holditch, D.E. Nierode, and R.W. Veatch Jr. eds., Vol. 12. Richardson, Texas: Monograph Series, SPE.
  6. Ramey, H.J. Jr. and Gringarten, A.C. 1975. Effect of High Volume Vertical Fractures on Geothermal Steam Well Behavior. Proc., Second United Nations Symposium on the Use and Development of Geothermal Energy, San Francisco.

Noteworthy papers in OnePetro

Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read

External links

Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro

See also

Fluid flow through permeable media

Type curves

Hydraulic fracturing

Hydraulic fracturing in tight gas reservoirs

PEH:Fluid_Flow_Through_Permeable_Media