Well Stimulation - Fracturing: Design of A Fracture Treatment Hydraulic Fracturing
Well Stimulation - Fracturing: Design of A Fracture Treatment Hydraulic Fracturing
Well Stimulation - Fracturing: Design of A Fracture Treatment Hydraulic Fracturing
Introduction
Design of a fracture treatment
Hydraulic fracturing
Fracturing pressure
Fundamental modes of fracture propagation
Rock stresses
Fracture geometry
Fracture length
Fracture area
Fracture width
Fracture fluid efficiency
Fracture fluids
Proppants
Operational Procedures
9.1 Introduction
There are three main reasons to fracture a well:
Heterogeneous Reservoir
Shale
1. Increased recovery from low permeability layers – most reservoirs have a non-
uniform distribution of permeability. Fracturing may increase the production from
layers of low permeability. Figure 9.2 illustrates this:
Low-Permeability Standstone
Fig. 9.2 – Fracturing in low permeability zone enhances production from that zone.
2. Longer field life – most wells are abandoned when the flow rate is no longer
economical. Fracturing may prolong the operating life of a well. This is illustrated
from an actual case in Fig. 9.3
INCREMENTAL PRODUCTION
PRODUCTION (BOPD)
DUE TO FRACTURING
PRE-FRAC
POST-FRAC PRODUCTION
PRODUCTION
TIME
• Fracture conductivity, kfw – this is a measure of how much reservoir fluid will
be carried through the fracture. Obviously, it is desirable to have a high
fracture conductivity, because it means more reservoir fluids can flow and a
higher rate of production results. (kf = fracture permeability; w = fracture
width)
There are other variables that influence design, such as well spacing, fracture azimuth,
fracture height, and transient flow conditions.
Because every reservoir is unique, so is the nature of the fracture job required for a
particular reservoir. It is not always the case that longer fracture lengths increase
productivity. For example, there is no point in fracturing beyond the limits of a
reservoir because we do not gain any additional production from beyond the reservoir.
Also it should not be taken that high permeability formations cannot gain from
fracturing – it may be the case that smaller, fatter fractures are required and the length
of fracture is less important.
kfw
FCD = (9.1)
kx f
We shall see that this variable is important in the design of a fracture treatment
process.
Actual wellbore
Vertical planer
fracture
A fracture = 2 ( 2 x f ) h
A eqborehole = 2π rw' h
A fracture = A eqborehole
∴ 2π rw' h = 4 x f h
i.e. , rw' = 0 . 64 x f
Thus for an infinite conductivity fracture, the effective wellbore radius is given as:
The above relation however, is a bit optimistic. In industry, the following relation is
used and is stated without proof:
From the chart we notice that as FCD approaches infinity, i.e. as the fracture becomes
infinitely conductive, the ratio of effective wellbore radius rw' to ½ fracture length xf
approaches 0.5. This is the same result we would expect for infinitely conductive
fractures, i.e. rw' = 0.5 xf.
Folds of Increase
The folds of increase (FOI) is a measure the improvement in production due to
fracturing. It is simply the ratio of the production after fracturing to the production
before fracturing, i.e.
r
ln e
FOI =
Q fractured
=
rw
(9.4)
Q unfractured
r
ln e'
rw
Obviously, the higher the FOI is, the more attractive a fracture job looks. But all this
comes at a cost. The final stage in the design process is to consider the economics
associated with well stimulation.
Cost/revenue information can then be calculated and a decision can be made based on
a suitable economic indicator. A commonly used economic indicator is the
Discounted Return On Investment (DROI):
FutureValue−Costs
DROI = (9.5)
Costs
The optimal design is based on the scenario that gives the best DROI.
2. Assume a value for kfw, xf (this can be obtained from a reservoir simulator).
4. Use value from charts/empirical data to determine rw' /xf ratio (eg., Cinco-Ley
chart).
5. Calculate rw' .
6. Calculate FOI.
8. Repeat steps 1 – 7 until optimal economic scenario is reached (i.e. the best
DROI).
The design process can easily be modelled in reservoir simulator. The design process
described above is very useful for steady-state and pseudo-steady-state reservoirs.
FOI must be used with extreme caution in the cases of transient flow reservoirs, as the
results may lead to erroneous conclusions. Fracture geometry and azimuth also must
be considered in the design process.
The pressure behavior of a fracture treatment is illustrated in the Fig. 9.7 below:
BREAKDOWN
PRESSURE
BOTTOMHOLE PRESSURE
FRACTURE PROPAGATION
PRESSURE
INSTANTANEOUS SHUT-IN
PRESSURE
PUMPING TIME
As fluids are injected into the reservoir with increasing pressure, three stages in the
pressure profile can be identified:
The pressure profile described above is highly idealistic. Different formations may
exhibit different pressure profiles for fracturing. The profiles will be governed by a
number of underlying factors. For example, if the formation was stimulated before,
there may little or no difference at all between the breakdown pressure and the
propagation pressure.
If the instantaneous shut in pressure (pISIP) is the pressure required to just hold open a
fracture at the surface, then the bottomhole shut-in pressure (pBISIP) can be calculated:
From the bottomhole shut-in pressure, we can find the fracture gradient (FG) which is
defined as:
p BISIP
FG = (9.7)
D
There are three fundamental modes in which a fracture will propagate. These are
illustrated in Fig. 9.8.
MODE III –
MODE II – In-planer Anti-planer shear
MODE I - Tensile shear dilation dilation
Fig. 9.8 – Fundamental modes of fracture propagation.
These are termed Mode I (tensile), Mode II (in-plane shear dilation), and Mode III
(anti-plane shear dilation). As will be discussed in the following section, a fracture
will propagate in the direction normal to the plane of least compressive stress.
Mode II- and Mode III-type fractures are at present impossible to initiate using
existing stimulation technology. However, Mode II and Mode III fracturing may
occur naturally in areas of high tectonic activity. We may, in future, develop the
technology to take advantage of these natural in-situ fractures. In subsequent
discussion, the focus of this course will be primarily on Mode I fractures.
σz
σx
σy z
Fig. 9.9 – Rock fracture shown perpendicular to the direction of least stress.
The nature of stress regimes in the subsurface is extremely complex and is governed
by depth and tectonic activity in the region. There are several theories used to
simplify and predict the nature of stress in the subsurface. One approach is the elastic
theory of materials:
v (9.8)
σ i − p = 2G (∈i + e)(i = x, y , z )
1 − 2v
E
G=
2(1 + v)
e = ∈x + ∈ y + ∈z
To really understand the stress regimes in the subsurface, one must get into a serious
study of the mechanical properties of rocks, tectonics and geology of the region of
interest. The fact that these areas are still subjects of active research (e.g., Mode II and
Mode III shear dilation) is testament to their complexity and is far beyond the scope
of this course. However, the parameters in the equations above can be found from
core samples and knowledge of tectonics and geology for a given region. Further
simplifications and assumptions can be made to find the compressive stresses in the
principal directions. Once the principal stresses have been identified, we can
determine the plane of least stress. Since a fracture will propagate normal to the plane
of least stress, we can predict the direction in which fracturing will occur.
In general, it is reasonable to assume that the stress regime will be governed primarily
by overburden stresses – and hence the depth of burial. This means that vertical
compressive stresses are usually greater than the horizontal stresses – unless there is
significant tectonic activity in the region which contributes to an increase in
horizontal compressive stresses. This implies that most fractures will be vertically
oriented and, in general, this is the case.
• Fracture length
• Fracture area
• Fracture width
The determination of the fracture geometry stems from the continuity equation:
Volume Pumped = qt p
Volume Lost = (3CH p L) t p + 2 S p H p L
Volume of Fracture = w HL
It can be shown that the length of the fracture can be expressed as:
qt p
L= = Tip- to -Tip Length (9.10)
3CH p t p + 2 S p H p + w H
Often, H = Hp and the height of a fracture govern the length. The height of a fracture
is extremely important to consider – especially when the fracture treatment is near oil-
water or oil/gas contacts. Under these conditions, if a fracture is allowed to grow
vertically, it may gain access to gas/water zones which may be undesirable from a
production point of view. Hence, if we know the height restrictions for a given
reservoir, we can establish the length which would give the desired height.
There are many equations that could be used to estimate the fracture area. These are
two that are often used in the industry:
qw x 2 2x
A= 2
e erfc( x) + − 1
4πC π (9.11)
πt
where x = 2C
w
qt p
A= (9.12)
3C t p + 2 S p + w
Like the equation for length, these equations are derived from the continuity
equation. It will be shown that the fracture area combined with the average
fracture width help to define the fracture fluid efficiency.
ELLIPTICAL CROSS
SECTION
AREA OF LARGEST
xf = L/2 FLOW RESISTANCE
xf
R = L/2
H
H
If a slit inside a rock formation is opened by internal pressure, a “slit” will open up to
an elliptical cross section, as shown in Fig. 9.11 below.
2d ( P − σ c )
w=
E'
d W
σC d = characteristic dimension (length or height)
P E' = plane strain modulus (≅ E, Young' s modulus)
P − σ c = net pressure
The characteristic dimension “d” in the mathematical model above is length, height,
or radius and is dependent on the fracture width model used. Table 9.1 gives a
general summary of the relations used in 2-D fracture geometry models.
Table 9.1 – Fracture width relations for 2-D fracture geometry models.
d (ellipse) H L R
qt p w w
Aw = = qt p
3C t p + 2S p + w 3C t p + 2S p + w
Aw = V FRACTURE , qt p = V INJECTED (9.13)
w
∴ V FRACTURE = V INJECTED
3C t p + 2S p + w
w
ef = (9.14)
3C t p + 2 S p + w
This is an important parameter in the design process as it can help determine the
volume of fluid required to complete a fracture treatment.
The key properties of a fracture fluid are its viscosity, fluid loss, and its ability to
transport and hold proppant.
Viscosity
Generally, viscosity plays two important roles in a fracturing process:
2. Provide sufficient fracture width to ensure proppant entry into the fracture.
Fluid Loss
Fluid loss is an extremely important variable in the design of a fracture treatment
process. It is measured by a fluid loss coefficient (C), which is a function of
formation permeability, reservoir temperature, and the “wall building
characteristics” of the fracturing fluid. Typically, values of C range between
0.0005 to 0.003 ft.min-½. Assuming that linear flow dominates, fluid loss can be
determined using the following equation:
CA
qloss = (9.15)
t −τ
Fluid additives must be used with extreme caution. Care must be taken so that the
use of additives has minimal effect on the conductivity of the fracture.
2. Stable viscosity – to minimise the need for using various fluids to achieve a
desired viscosity over the entire fracture length.
6. Short breaking time – to minimize the shut in time (after the fracture closes)
before the fluid viscosity “breaks” to a lower viscosity and the well can flow.
Of course not all fluids satisfy these desirable qualities. There are always
compromises made due to economic and design constraints.
9.8 Proppants
After a fracture has been created, it must be supported to hold it open when the
instantaneous shut-in pressure is removed – this is the role of proppant. The proppant
provides support to fractures by filling up the fracture volume, for this reason the
proppant must have high permeability so that the fluid conductivity remains high
through the fracture. To address these needs, there are industry set guidelines for
suitable proppants and selection. Some important physical properties are listed based
on API recommended practices:
1. Grain size and grain-size distribution – will have an effect on the fracture
conductivity. Proppants with larger grain sizes provide a more permeable pack.
However, their use must be evaluated in relation to the formation that os propped
and the increased difficulties encountered in proppant transport and placement.
Larger grain sizes can be more difficult to use in deeper wells because of greater
susceptibility to crushing. Figure 9.12 illustrates the effect of grain size at
increasing closure stresses.
3. Turbidity – is a measure of the presence of silt, clay and any other fine particle
contaminants.
f = (1 − e 2 ) + f (9.16)
p f c
The term ef in Eq. (9.16) is the fracture fluid efficiency as defined by Eq. (9.14).
Banking Fluids
Fluids with viscosity less than 40 cp are termed banking fluids. Banking fluids do not
have a high enough viscosity to suspend solids, therefore the proppant will settle as
the pad is being replenished. This practice is illustrated below:
PAD
Replenished PAD
The goal of this procedure is for the proppant to settle near the wellbore. Thus, the
final conductivity is a function of fracture width and proppant settling rate, and less
dependent on the final proppant concentration.
The two general types of fracture treatment schedules are illustrated in Fig. 9.15:
V1 and V2 are the pad volumes for a tip screenout and standard fracture treatment
respectively.
REFERENCES
1. Economides, M.J., and Nolte, K.G., 1989: “Reservoir Stimulation”, 2nd Edition.
Schlumberger Educational Services, Houston, TX.
2. “Hydraulic Fracturing” training manual by NSI Technologies. Inc. 2nd Edition, 1992.