Hydraulic Fracturing
Hydraulic Fracturing
Hydraulic Fracturing
Fracturing Overview
Hydraulic fracturing involves the
injection of fluids into the formation
at a rate and pressure above the
fracture pressure of the reservoir in
order to create a fracture within the
rock itself. We then pack this space
with proppant which:
Improves productivity
Interconnects formation permeability
Improves EUR (Estimated Ultimate Recovery)
3
Rekahan me-bypass
damage
Rekahan menjalar
10
11
open fracture
1/2"
during job
(frac width = wf)
sand used to
prop the
frac open
Hydraulic Fracturing
Frac perm. = kf
lxf = Fracture half length
}
lxf
wf
kf wf
Fcd
k elxf
Where;
kf = fracture permeability
Fracture conductivity
wf = fracture width
------------------------------ke = formation permeability
Formation conductivity
lxf = fracture half length
rw = 0.5*lxf
15
rw = 0.5*xf
rw = 0.5*800 ft
rw = 400 ft
width = 800 ft!!!
For a half length of 800 feet,
800 ft.
Hydraulic Fracturing
Materials
Base Fluid Systems
Chemical Additives
Proppants
17
Slickwater Applications
Low Friction
Low Viscosity (<5cp)
Low Residue, less damaging
Low Proppant Transport
capabilities
Linear Gel Applications
Mild Friction Pressures
Adjustable Viscosity (10<x<60cp)
High Residue, more damaging
Crosslinked Applications
High Friction
High Viscosity (>100cp)
Excellent Proppant Transport
capabilities
High Residue, more damaging
Expensive
Complex Chemical Systems
pH and Temperature
dependent
18
Chemical Additives
Gelling Agents
Friction Reducers
Crosslinker Control
pH Adjusting Agents
Clay Control
Breakers
Scale Inhibitors
Corrosion Inhibitors
Bactericide
Oxygen Scavengers
Surfactants
Recovery Agents
Foaming Agents
Acids
Anti-Sludge Agents
Emulsifiers
Fluid Loss Agents
Resin
Activator
Multiple varieties of all additives
differing
in
concentrations, ionic charges, limitations, and
compatibilities
19
Proppant Permeability
20
Proppants
21
Horizontal fracture
Vertical fracture
23
Change streamlines
Radial flow disappears
Wellbore radius is not a
factor any more
Increased PI can be
utilized
p or q
q J post p
24
Fracture Morphology
25
Upwar
d
fractur
e
growth
Horizont
al
fractures
?
?
?
Out-ofzone
growth
?
Twisting
fractures
Pay
Pay
?
Perfect
fracture
Multiple fractures
dipping from vertical
T-shaped
fractures
Main questions
Which wellbore-fracture orientation is
favorable?
Which can be done?
How large should the treatment be?
What part of the proppant will reach the pay?
Width and length (optimum dimensions)?
How can it be realized?
27
Circular:
2kh
J D p
q
B
JD
1
re 3
s
ln
rw 4
Drawdown
Dimensionless
Productivity Index
28
Hawkins formula
k
rs
s
1 ln
rw
ks
ks
Damage
penetration
distance
rw
rs
29
Exercise 1
Calculate the skin factor due to radial damage if
rs
Wellbore radius
0.328 ft
Permeability impairment
k
5 folds
ks
0.5 ft
Damage penetration
Solution of Exercise 1
rs
s
1 ln
rw
ks
rs 0.828 ft
0.828
s 5 1 ln[
] 3.7
0.328
Exercise 2
Assume pseudo-steady state and drainage radius r e = 2980 ft in
Exercise 1. What portion of the pressure drawdown is lost in the
skin zone? What is the damage ratio? What is the flow efficiency?
Solution 2
The fraction of pressure drawdown in the skin zone is given by (Since we
deal only with ratios, we do not have to convert units.):
3.7
2980
ln[
] 0.75 3.7
0.328
0.31
Exercise 3
Assume that the well of Exercise 2 has been matrix acidized and the
original permeability has been restored in the skin zone.
What will be the folds of increase in the Productivity Index?
(What will be the folds of increase in production rate assuming the
pressure drawdown is the same before and after the treatment?)
Solution 3
We can assume that the skin after the acidizing
treatment becomes zero. Then the folds of
increase is:
Folds of Increase :
r
ln[ e ] 0.75 s
rw
FOI
r
ln[ e ] 0.75
rw
2980
3.7
0.328
1.44
2980
0.75 ln
0.328
0.75 ln
32
Exercise 4
Assume that the well of Exercise 2 has been fracture treated and a negative pseudo
skin factor has been created: sf = -5. What will be the folds of increase in the
Productivity Index with respect to the damaged well?
Solution 4
The ratio of Productivity Indices after and before the treatment is
2980
] 0.75 3.7
FOI 0.328
3.6
2980
ln[
] 0.75 5
0.328
ln[
2Vfp
2xf
34
Dimensionless fracture
conductivity
2 xf
Dimensionless fracture
conductivity
C fD
kf w
kx f
fracture conductivity
no name
35
q Jp
2kh
2kh
1
J D
J
B ln[ re ] 0.75 s
B
f
rw
JD is a function of what?
half-length,
dimensionless fracture
conductivity
Drainage radius, re sf is a function of what?
half-length,
dimensionless fracture
conductivity
wellbore radius, rw
36
2kh
re
B ln 0.472 s f
rw
or
2kh
re
B ln 0.472
r 'w
Prats
f (C fD )
2kh
J
B
xf
0.472re
ln
s f ln
xf
rw
2kh
0.472re
ln
f
xf
Cinco-Ley
37
Notation
rw
r'w
f s f ln
xf
rw
sf
rw
xf
Dimensionless Productivity
Index, sf and f and rw
JD
1
re
ln 0.472 s f
rw
or
JD
1
re
ln 0.472
r 'w
Prats
f (C fD )
1
1
JD
0.472re
xf
0.472re
ln
f
ln
s f ln
xf
xf
rw
Cinco-Ley
39
Penetration Ratio
Dimensionless Fracture Conductivity
Proppant Number
Ix
2x f
C fD
N prop
xe
y e = xe
kf w
2 xf
kx f
4k f V f,prop,1 wing
kVres
xe
2k f V f,prop,2 wing
kVres
(I x )2C fD
40
41
rw'
0.5
xf
'
w
r
xf
0.1
rw'
0.25 C fD
xf
0.01
0.1
1.0
C fD
kf w
kx f
10
100
42
1.65-0.328u+0.116u 2
f (C fD )
1+0.18u+0.064u 2+0.005u 3
where u ln C fD
1
use f = ln(2) for CfD > 1000
0
0.1
10
CfD
100
1000
43
44
N prop I C fD
2
x
N prop
4k f x f w
2
e
kx
4k f V1 wing , propped
2
e
kx h
2k f V2 wing , propped
kVreservoir
45
46
47
V fp hw p x f
C fD
k f wp
kx f
V fp k f
1/ 2
C fDV fp k
1/ 2
xf
C hk
fD
wp
hk f
48
0.5
X e=Y e
Ye
0.4
Ix =1
2Xf
Medium perm
0.1
Xe
0.06
0.03
0.3
0.01
0.006
High perm
Frac&Pack
0.003
0.001
0.2
0.0006
0.0003
N prop=0.0001
-4
10
-3
10
-2
-1
10
10
10
Dimensionless Fracture Conductivity, CfD
10
10
49
X e=Y e
Ye
Ix =1
100
2Xf
1.5
60
Xe
30
10
Low perm
Massive HF
1.0
3
1
0.6
0.3
0.5
Medium perm
N prop =0.1
0.1
10
100
1000
50
Ix=1
Ye
0.4
2Xf
0.1
Xe
0.06
0.03
0.3
0.01
0.006
0.003
0.001
0.2
0.0006
0.0003
N prop =0.0001
-4
10
-3
10
-2
-1
10
10
10
Dimensionless Fracture Conductivity, CfD
10
10
CfDopt=1.6
52
1.8
1.6
30
X e=Y e
Ye
10
2Xf
Xe
1.4
1.2
1.0
1
0.8
0.6
0.6
0.3
0.4
0.01
N prop =0.1
0.1
Penetration Rate, IX
53
xf
V fp k f
1/ 2
C fDoptV fp k
wp
hk
fDopt
1/ 2
hk f
V fp k f
x f
1.6hk
1/ 2
wp
1.6V fp k
1/ 2
hk f
54
Economic optimization
Production forecast
Transient regime
Stabilized
Economics: Converting additional production
into value
Time value of money
Discounted revenue
NPV
55
Treatment Sizing
N
Rev n
- Cost
n
n 1 (1 i)
NPV
57
Rock Properties
Youngs modulus, Poisson ratio,
Fracture toughness, poroelastic const
Stress State
Leakoff
Proppant and Other Fluid properties
Operational constraints
58
Rock Properties
Linear Elasticity
Poroelasticity
Fracture Mechanics
59
xx l
l
yy D
D
xx
E
xx
l
D
xx F
A
yy
xx
D/2
G,
E ,G
Shear modulus, G
E
21
Young's modulus, E
2G 1
Poisson ratio,
E 2G
2G
E
1 2
2G
1
4G 2
4G E
Required
Known
61
Formation Classification
Two types
Consolidated and tight
E = 106 + psi
Unconsolidated and soft E = 105 - psi
62
p
Total Stress = Effective Stress + [Pore Pressure]
63
Force
Biots constant
Pore Fluid
64
65
v v p
v p
h
1
'
v p p
h
1
1) Poisson ratio changes from layer to layer
2) Pore pressure changes in time
66
0
Ground Surface
Critical Depth
977 m
-1000
m
mu
-1500
Tr
u
Ho
riz
-2000
on
tal
s
es
St r
-2500
-3000
0
20x106
O
rig
i
Ve
rt
ic
al
-500
-1000
na
lV
er
tic
St
re
ss
-1500
al
40x106
60x106
Stress, Pa
St
re
ss
Current Depth , m
-500
ni
Mi
-2000
-2500
80x106
67
Stress Gradients
Overburden gradient gradient
Slope of the Vertical Stress line
1.1 psi/ft
Frac gradient
69
E
1 2
Half length c
pn(x)
Deformation (distribution)
net pressure (distribution)
70
Width
4 pn
w( x)
c2 x2
E'
Max Width
4c
w0
pn
E'
pn : net pressure
c : half length
characteristic dimension
c
w
linearity preserved
71
Questions:
Upper tip
Pinch point
Lower tip
Contained?
Breakthrough?
Run-away?
Up or Down?
Width?
Hydrostatic
pressure?
Height
control?
What can be
measured?
72
73
K I 2c
pn ( x )
c x
2
dx
1
cx
psi - in.1/2
Weighting function: the nearer
to tip, the more important the
pressure value
x
c
KI : proportionality const
74
Stability of Crack,
Propagation
Critical value of stress intensity
factor:
Fracture Toughness KIC
Propagation: when stress
intensity factor is larger than
fracture toughness
75
Application:
Fracture Height Prediction
Height containment: why is it critical?
Can it be controlled?
Passive: safety limit on injection
pressure
Active: proppant (light and heavy)
76
fluid pressure
p
r
o
f
i
l
e
77
Injection rate
Bottomhole pressure
Time
78
Bottomhole pressure
Injection rate
79
3 ISIP
Fall-off (minifrac)
4 Closure
5 Reopening
6 Forced closure
8 Rebound
3
2nD injection
cycle
Injection rate
Injection rate
1st
injection
cycle
Bottomhole pressure
shut-in
flow-back
Time
80
Proppant Placement
Concepts
81
No settling in perfect
transport fluid
Viscosity (rheology) and density
difference
(Foams: visc good, dens: bad)
82
Two approximations:
Perkins-Kern-(Nordgren)
Vertical plane strain
characteristic half-length ( c ) is half
height, h/2
elliptic cross section
Kristianovich-Zheltov (Gertsmaa-deKlerk)
Horizontal plane strain
characteristic half length ( c ) is xf
rectangular cross section
84
Vf = w(h f x f )
Perkins-Kern-Nordgren PKN
ww,0 = 3.27
qi x f
E'
1/ 4
w 0.628ww, 0
Kristianovich-Zheltov
Geertsma-De-Klerk KGD
qi x
ww = 3.22
E' h
f
w 0.785ww
1/ 4
2
f
85
1
2n 2
ww, 0
3.98
n
2n2
1 2.14n
n
2n2
1
2n2
qi h1f n x f
1
2n 2
E'
86
Vf = w(h f x f )
Vf = w A
2qi
V i = qi t e
xf
Vfe = Vi - Vlost
Average
w(xf)
hf
qi
A
Lost: spurt +leakoff
87
Given:
Mass of proppant, target length, frac height, inj
rate, rheology, elasticity modulus, leakoff coeff,
max-possible-proppant-added-conc
88
hf
hp
rp= hp /hf
89
Design Outcome
Constraints allow optimum placement
of the given amount of proppant
Some improvement is necessary
Consider higher quality proppant
Better fluid loss control
Better rheology
Larger allowable proppant concentration
90
Width to accept
proppant
At the end of pad stage the
created width has to be at least
2-3 times the proppant diameter
At the end of pumping the
proppant reaches only that part
which has a width at least 2-3
times the proppant diameter
Propped length less than
hydraulic length
92
Width ratio
criterion
Considering material coordinate,
Accounting for fluid loss
Calculate ratio of (Dry width) to
(Dynamic width)
Criterion: cannot exceed critical
value (about 0.5)
93
94
Hydraulic Horsepower
Energy: (Power) (Time)
Power = (Pumping Pressure) (Injection
rate)
(Pumping Pressure) =
Minimum Stress + Net Pressure +
Friction Losses - Hydrostatic Pressure
Friction Losses : in tubulars, through
perforations and possibly in near
wellbore tortuous flow path
95
Limited entry
Few perforations in small groups
High perforation friction loss
Uniform coverage
97
Monitoring
Calculate proppant concentration at
bottom (shift)
Calculate bottomhole injection
pressure, net pressure
Calculate proppant in formation,
proppant in well
Later: Add and synchronize gauge
pressure
98
Nolte-Smith plot
Log net
pressur
e
Normal
frac
propagatio
n
Tip
screenout
Wellbore
screenout
Unconfined
height
growth
Log injection
time
99
Sizing
Pad volume for generic design
More aggressive or defensive proppant
schedule
Proppant change (resin coated, high
strength etc.)
Fluid system modification (crosslinked,
foam)
Proppant carrying capacity
Leakoff
100
Or
500 ft
450 ft
Pay
zone
What we want
What we get
102