Lumped Parameter Leakoff Model
Lumped Parameter Leakoff Model
Lumped Parameter Leakoff Model
SPE 15069
A Complete Integrated Model for Design and RealTime Analysis
of Hydraulic Fracturing Operations
by AR. Crockett and N.M. Okusu, Resources
Massachusetts
Inst. of Technology
SPE Members
~opyright 1986, Society of Petroleum Engineers
rhis paper was prepared for presentation at the 5Sth Califomis Regional Meeting of the Society of Petroleum Engineers held m Oakland, CA, April 24,
!9ss,
rhis pa~r was selected for presentation by an SPE Program Commmes following review of mlormation contained in an abstract submitted by the
suthor(a). Contents of the psper, as presented, have not bean reviewad by the Sosiefy of Petroleum Enginaare and are aubjact to correction by tha
guthor(a). The material, aa presented, doas not necessarily reflect any poaifion of the $ociaty of Petroleum Engineers, its offiiere, or members. Papera
praaentad at SPE meetings are subject to publication review by Editorial Committees of the Socialy of Petroleum Enginaera. Permission to copy is
restricted to an abstract of not more than 200 words, Illustrations may not be copied. The abstract should contain conspicuous acknowledgment of where
nndby whom the pa~r is presented. Write Publicatiina Manager, SPE, P.O. Sox 83SSSS,FUchardaon, TX 7S0S3-3S3S. Telex, 7309S9, SPEDAL.
Abstract
Introduction
A new comprehensive model of hydraulic fracturing is presented which has been developed for the Gas Research Institute
(GRI) mobile fracture monitoring and analysis facility. The
main purpose of the model is to simulate the hydraulic fracturing process in real-time, that is on+ite during the fracturing
operation, but the model can also be used for pre-fracture design and post-fracture analysis. Sensor data obtained during
the course of the job such as wellhead pressure, flow rates,
frac-fluid viscosity, and proppant staging can be received
directly by the model se input, superseding the prefrac job
design schedule, and making possible more accurate model estimatea of current fracturing conditions and predictions of final
fracture geometry, as the job proceeds.
A substantial amount of effort has been invested, especially over the past decade, in the development of models for
design and analysis of hydraulic fracturing. The resulting models have varied in at least three major aspects: realism and
generality of the assumptions made in formulating the modela; complexity of the resulting computer codes and machine
requirements; and flexibility of the input-output characterietics, especially in relation to real job conditions and operator
.
interfacing. Although some good progress haa been made by
many groups, there have not been any models which were satisfactory in all of the three important areas, and most models
h ave, at beet, been adequate in one aspect only.
Examples of previous work may be found in the considThe overall model haa four major components describing:
erable literature which has evolved on this subject. The sim flow of fluids and slurry in tubular goods
plest models, assuming 2D geometry with a constant specified
creation and propagation of the hydraulic fracture
height, were those of Christianovich, Geertsma, de Klerk and
transport of proppant, deposition, and fracture closure
Daneshy (CGD, Refs. 1, 2), having width dependent on length,
heat and fluid exchange between fracture and reservoir.
and those of Perkine, Kern and Nordgren (PKN, Refs. 3, 4),
This fully-integrated, numerically robust model of the hy- having width dependent on height: these two models are really
draulic frauturing process takes directly into account as much valid only when length greatly exceeds height (or vita-versa if
of the essential physics as possible, given the computational
the models are turned on end) so clearly neither are approlimitations of the real-time application. Additional informa- priate for the typical field condition where the fracture is not
tion, pertaining to very complex reservoir characteristics, can very @l contained.
be indirectly supplied to the model through data-based results
Since both the height and length may grow substantially
obtained prior to the job from other more comprehensive [e,g.,
3D, croes+wctional) fracture simulations. Having thereby em- during the job in contradiction to these 2D models an
a priori specification of height is hardly acceptable, even if
bedded more elaborate fracture analysis into a simple lumped
based
on (usually tenuow) deductions from well 10SS. To remodel, accurate pre~lctions are achieved at execution speeds
faster than real-time, allowing on-the-job analysis and even solve thm problem, Cleary (Ref. 5) developed a pseudo-three.
dlmeneional hydrafrac model (P3DH) of simultaneow length
real-time history matching for unknown reservoir parameters.
and height growth, having width dependent mainly on the
Sample results are presented from model simulations of a
lesser of height or length. This model wee adopted by varifracture treatment performed in the Travis Peak formation of
ous groups in industry (e.g., Refis. 6, 7) while others developed
East Texas. Actual sensor data from the job was used as input
specialii
variationa on the same theme (e.g., R@. 8, 9). Aland model predictions are compared with field measurements.
though these P3DH-type models, if fully developed, do indeed
References and illustrations at end of paper.
9
m
.A
. fY3MPl,F!TF!
. . . . . . ----
REAL-TIME
.. ..
1. Frac-fluid/slurry
flow in tubular ~oods and perforations, including effects of cross-linking, foaming and proppant
concentration in gelled fluid flow rheology. This incorporates
both laboratory and field data on pipe flow, the latter based
on simultaneous uphole and downhole pressure measurement
wherever such data has been collected: a major purpose is
to eliminate the need for downhole data collection, by being
confidently able to calculate true downhole fracturing pressure
from measured uphole pressures.
of Model
2fv2
p
=exp(~)~l+~
(1)
Pr
The foregoing 3D, P3DH and 2D efforts have substanally enhanced abilities to model many aspects of hydrafrac;
owever, for different reasons, none has produced an adequate
~odel for realtime analysis of fracture operations. Thus, the
urpose of the model described in this paper is not only to be
sed for comprehensive hydrafrsc design, but also for real-time
nalysis of the treatment. To achi~ve this, the models had to
atisfy the three major conditions ~iet!orih at the outset: they
re realistic and general in the I,hysics and geometry which
hey represent; they run effectively on all levels of computer
ystem, from mainframe to portable; and they operate on the
ctual data generated by the job itself, running in real-time as
ceded and providing an instantaneous view of the job status to
he operator including pre-frac design and post-frac analyis. By offering all of these essential capabilities, the rnodelling
f hydrafrac finally takes on a greater and more meaningful
ractical significance.
Overall Description
SPE 15069
OF HYDRAULIC FRACTURING
MODEL
its solution is
dp
~-C1p=CO;
Pj=
co =
(Pr9
Cl = (Prg +
Pj-1(1
cl~B)
cO~B,
(2b]
uj-1)2/Pr~)/c2,
zf(Pj-1
zf(Pj-lvj-1)2/Pr~)
(Z4
/~c2t
(2C]
and
(2d:
C2 = 1 (~j-lvj-1)2/prB
Capabilities
The overall modelling capability consists of four fully integrated modules, each representing a major aspect of the hydraulic fracturing process. In order to be run during the fracture treatment at rates faster than real-.me and :.ing aCtual sensor data as input, these modules weie .. . kematical ly
formulated to accept a wide range of input values, and were
efficiently implemented on the computer to minimize execution time. None of thw was accomplished, however, at t he
expense of accuracy or completeness; on the contrary, we ha ve
taken considerz- ; care to incorporate as much of the easential physics as possible dkectly in the modules an indir ect
method of data-baaing results provides a means to account for
whatever remains. The physics which have been incorporateed
into each of these modules is described in the following four
sect ions:
en
MB =
/0
pAwdzz
(38
yielding
(Pj-lcl
co)t~
+ (~+ pj_,/B)tB
2B
999
.-.
=o.45-
2,75
n
(3b
deification of Clapp
B
-=0.
prAw
4.-.
n
q(ReCfC
(n-.2)/2)
(4a
RPF.
a
-
R flRf)f!KRTT-
1*V
MIRO
.A. . . . .
-A
-----
-.
-------
..-
.-.
-----
E =
acbclOg
Re
R,,
= ;;:_:m
(4b)
n which p: is the apparent viscosity for gel and s is the volunetric sand concentration. The effect of foam quality on visosity is accounted for in a similar fashion; however, the precise
orm of this effect has not yet been finalised,
!. Fracture creation. extension. shutin and closure:
rhis component is a completely general integrated 3D model of
Lydraulic fracturing, incorporating explicitly to first crder all
~fthe essential physics, including the rock mechanics of deformation and fracture, the coupled fluid flow within the fracture,
md the fluid and heat exchange with the surrounding reserfoir. Any required degree of accuracy may be obtained by
lmbedding, in functions or data-bases, the results of comprehensive (e, g,, fully 3D or cross%ectional) fracture simulations
nto the integration parameters or gamma (-Y)factors of the
nodel, thus taking into account the effects of complex reser{oircharacteristics (i. c,, multiple strata, inhomogeneities, etc.)
m the spatial variation of fracture attributes (i. e., crack openng, excess pressure, fracture height, etc.). Default values for
~amma factors have also been determined, which are either
:onstants when variation in the factors does not greatly afkct fracture propagation, or functions of relevant parameters
when the parameter variation dominantly affects the fracture
configuration (i. e., confinement, frictional drag, etc.). Having
relegated the complex numerical analysis to the determination
~f functions or data-bases for the integration parameters, it
is possible to perform accurate, faster-than-real-time
analyses during the fracture treatment, as well as routine pre- and
post-frac analyses on PC-based systems.
A 2n+1
~pF
~Xi
i=l,2
(7a)
+2/n)*
r = K~n
(7b)
(8a)
where the mass flow rate into the fracture wing is divided by
the fracture height, and
(8b)
where W2equals the time rate of change of the vertical crosssection of the fracture plus the mass loss rate per unit length.
7x and 72s are the vertical crocxwwction and fluid loss shape
f~tors respectively. Specifying WI as the mass flow rate input
W, however, reduces the order of the governing set of equations
by one; this causes discontinuous jumps in estimates oi fracture
height and width when the flowrate W varies with time. We
use another closure condition similar to Eqn. (8b) for both
lateral and vertical directions
7i4
(5)
Wi =
PF
~is A hi +
7i6 wL/Lkv
i=l,2;
k=2,1
(8c)
-----
..- .-.
A _ ~ .
i
(Ioc)
[[71 -Y&D
2pF7u ~
~
nd
Bi = N_
~w2~w~)
i = 1, 2; k = 2, 1. (lOd)
Once the length and height are determined from Eqns. 10,
le other parametem, such as the width and excess pressure,
}Iloweasily from bad substitution into Eqns, 5-9. In actualy, there are three differential equations governing length and
oth upper and lower height extensions as shown in Fig. 1; our
lodel possesses this general feature, but, for simplicity, only
~e symmetric fracture formulation is presented.
C1early, the fluid 10SS2WL plays a dominant role in the
dutiona just derived and fluid 10SSwill be discussed in the
ext section. However, it is interesting to derive some analyt:al results which can be obtained in the special case where a
Vewtonian fluid is assumed (n = 1) and fluid 10SSis neglected.
rhe solution for a circular fracture (no confinement) is found
,0 be
=
A=
(ha)
[:-:[+13],
[[i%w2[a]
*d
E 4 T1T13 P lt_~
71 [9 qlzqlq ~ 1
(llC
1=
[[:-lwi]i~
[:-;(&)2+
(12a)
qis = Li7~8*
%)9
12
1* 29
(13b)
(13C)
(12b)
md
u=%~*:($2#
(llb)
md
0=-=
(12c)
---
-.
----
~een finalised; however, to consolidate ideas, several illustraive examples based on cross+ectional analyzes are described
~elow:
1) If the confining stress is greater in the adjacent strata
,han in the reservoir by an amount Auc, vertical height growth
s impeded. Such a step stress contrast affects the pressure
yadient 7i2 and the crack opening 71 coefficients when Lz >
H,; the corresponding S-factors take on the following forms:
S~l = ~ [~- AoC[l -e~p(~.i(l
S82= : u[[
AuC 1-
- ~))]],
(14a)
sin-(~,z~)
sin-l (0,2)
(14b)
1]
and
s, = S*2+
2 Au=
;~~,$b!n
l-t
(1-(LLE)2)*
~)2
Spi =
p{ew]
p(eR)
+
viscosity variation
kwlexp(-~ia
f)
+
.
[ 1- (1-(/3,
(14C)
$1
% = 1 + [~ -1]
L
exp(-~c(%
L1/Lz, Lz < LI
= { L2/Ll, L1 < L2
1)),
(lg)
~~
PR 6h 6L
PC
&
Since there is no equilibrium height, the modulus contras t relating the fracture pressure pp to the interface preesure PI
can only impede, not stop fkacture growth; consequently, in th e The viscosity of the fluid occupying each of the zones is mitabll
limit of a very large fracture the adjacent strata characterizti Cs averaged to take into account its thermal sensitivity
again predominate arid the fracture grows as if it seen only
1
(lsb;
~a (0) dzs
the modulus of the adjacent strata in thw knit both S~l an d
*=C6
/
sdz must have returned to a value of unity and $j ta E/EA
in
where a may stand for the filtercake, the leak~ff, or the ther
order to have replaced E with ~A everywhere in the govemin g
mal zone (C, L, and h subscript respectively).
equations. The exponential coefiiciente fid and ~di control th e
rate of thu replacement and are determined from the rnore
detailed &nalyaes.
Zza
A CX3MPLETE
REAL-TIME
----.-- -
.
PR
~ qL (r)
7~ (~, T) d?,
Yf(h ~) = *
[4n :_
~)
(15C)
eff-t3(Z3,
6L=A ~d
t
~ qL
(~)
7h(~3,
The
T)] ~h (z3,
dX d?,
Z; t, T)
z; ~, ~)
(a)
(lb)
(P C) FVF
+ AFhL
(eF
eR)]
- (@) UILAF(eF
(Pc)wd(%
eR)
(17C)
- eR)
(18a)
OB =Ug+u:+ug
where KL is the overall fluid loss coefficient. This !OSScoefficient is related to filtercake, kink-off, and reservoir zone properties, in addition to the preaure dMerence PF - pff, by the
following quadratic expression in KL:
(18b)
Ug = 2J#CCi.?dL/%t,
= 2fi~(1 - #C,)~(pF
0;
u:
+pR@iZ
~R
O; t, r)-
~L (k)
$[&(-@c)+
eff e (Z,
(15e)
(lsf)
-qL(~J
~ hL (~) ~~ (2s,
except there is a second integral representing the effect of thermal convection which is absent in the case of fiuid-diffusion.
The effect of fracture area creation is taken into account just
as in fluid loss:
Yote that since the depth of the frac-pressure diffusion into the
:eservoir is typically small relative to the fracture length during
;he treatment time, it is sufficient to aesume ID loss normal to
;he fracture, provided that the creation of new fracture area is
Lakeninto account; this is accomplished by dividing the volume
!bssrate which occured at time level tm by the current fracture
mea Ap(tN):
qL (tn)
A
~k
[-rexp[~~~:l
W)
d~.
I%[r&[
vhere If=lfor
6~-6rj>Oand
H=0 for6h-6LS0.
eak-off penetration can be found from the expression
t)
+ [@ce)+&cc],L
SPE 15069
bn
&(PI -
PR)[l
ep(-t-2)
-pi)/?rt(l
-t
- V),
$erfc
(~)]
(MC)
(18d)
@c]K~
KL(pF-pR)=O.
(16b )
---
where ~ = 2 @/t
is the dimeaionless diffwion time; 4 is the
lesser of the fracture height or length; q = (1 - ~/&)(l
2u)/2(1 - v) is the poroelaetic induced stress factor; and b =
0.56 is a numerical factor (Ref. 20). Thermally-induced strasws
do not greatly affect hydraulic fracture growth, except in steam
and water fioodhg, because the time scale of the thermal
convective-&ffuaion in reservoir rock is much larger thau the
fracture growth time.
vP.
t
=
Q-
/[
(19)
Qbn]~r
tpm
where tpm is the time when the rnth proppant stage entered
the fracture, and Q~~ is the volume leak~ff rate through the
hcture area adjacent to VP-. Once VP- is obtained, the extension into the fracture of each concentration band Lpm is
Foundfrom mess conservation yielding
v
(723 -
%)
(%)3
-723
(%)2
(Zoa)
(Pp-Pf)~
(h+l)d
108 n
l(J1.s2(l-f.)
72 K
(21a)
and
Vef:
h
(21b)
Sample
Implementation
of the Models
To illustrate the capabilities of the overall integrated motlel, we analyse a fracture treatment performed in the Travis
Peak formation of Eeet Texas on which the models were implemented in real-time, using the Gas Research Institute (GRI)
fracture treatment monitoring and analysis facility developed
by Resources Engineering Systems (RES) and described in
Ref. 22. The well was one of many monitored and analysed
under the auspices of the GRI Tight Gas Sand Program, in conjunction with a number of ~,ther GRI contractors (e.g., Ref. 23).
A COMPLETE
REAL-TIME
Conclusions
Nomenclature
AF
A
SPE 1506$
cc
c~
CW, CF, CL
E
E
erf
erfc
crack-opening moduli for upper, lower and adjacent strata (when symmetric)
friction factor for wellbore flow
f.
9
&
gravitational constant
half height to modulus barrier
hL
Ha
&J,
EL,
EA
>W
A .R
. .
I KnRo
LJ
. .
J.uvvcz
fiRfW!Kl?TT.
----
-----
_..
_..
!Lm
L
P
P.
PC,
pL>
#R
Poissons ratio
tc=
d, Sdi
p, pF
porosity
wellbore, fracture slurry densities
Pp,
inl
Pr
*.
La
Pw,
~,
Jpa
,)#j
,
pi
time
time of entry of the mth proppant stage into
the fracture
Pm
Pf
PL
temperature
uB
u:, Ug, Ug
{F
fracture volume
/Pm
/
u.
AOCU,AOCL
l,, vah
Au.
shear stress
*,
ifiL , 2WL
01, tiz
cl, Z2, z~
dc
~d, ~di
flev@.i
0Hi
Acknowledgements
The work presented in this paper has been performed as part
of the GRI Tight Gas Sands Project, a comprehensive effort to
imp we many aspects of hydraulic fracturing operations. We
are grateful to H, D, Vo and R. M. Willis for their assistance in
implementing the model on the computer and to Marie Dokoupil for helping to prepare this manuscript.
Figures
Yf
Yh
7P
7.
71
712, 722
?13,
723
714,
724
%S, 726
A
6C, lfL, dh
gamma factors for lateral, vertical stress gra- 4(a). Flow, rheology, and proppant data; time in minutes.
dient
gamma factors for lateral, vertical cross~ectio n 4(b). Uphole pressure data showing comparison of model w
actual pressurea; time in minutes.
area
gamma factors for lateral, vertical channel flow 4(c).
gamma factom for fluid 10ss
width of fracture at wellbore
4(d).
filtercake, leak-off and thermal penetration
deptha
~L
Zzv
REAL-TIME
HYDRAULIC
.-- OF
-
.A. COMP1.F!TR
-----. -.
-. -.-..- MODEL
References
. Geertsma, J. and F. de Klerk, A Rapid Method of Predicting Width and Extent of Hydraulically-Induced Fracturesn,
J. Rt. Tech., Dec. 1969, pp. 1571-1581.
FRACTURING
SPE 15069
nisrna and Procedures for Producing Favorable Shapes of 18. Crockett, A. R., Willis, R. M. and M. P. Cleary, lrnprovHydraulic Fractures ).
ing Hydraulic Fracture Predictions by Real-Time HietoryMatching on Observed Pressures, SPE Paper No. 15264,
1. Settari, A. and M. P, Cleary, Three-Dimensional Simulapresented at the SPE Unconventional G= Recovery Sym.
tion of Hydraulic Fracturing, J. of Rt. Tech., July 1984,
posium in Louisville, KY, May 1986.
pp. 1177-1190 (see also SPE Paper No. 10505, Development and Testing of a Pseudo-Three-Dimensional Model of 19. Settari, A., A General Model of Fluid Loss in Hydraulic
Hydraulic Fracture Geometry, presented at the SPE SymFracturingn, SPE/DOE Paper No. 11625, presented at
posium on Numerical Simulation, New Orleans, February
the SPE/DOE Joint Symposium on Low Permeabiliw Gas
1982).
Reservoirs, Denver, March 1983.
~. Meyer, B. R., Frac Model in 3D - Parts 1-4, Oil and Gas 20. Crockett, A. R., Vo, H., and M, P. Cleary, Studies of Fluid
Journal,,
June 17 (p, 87), July 1 (p. 62), July 22 (p. 83) and
Flow, Heat Transfer and Induced Stresses In and Around
JUIY 29 (p. 132), 1985.
Underground Fractures, Report No. REL-84-1, MIT Resource
Extraction Laboratory, March 1984.
J. Palmer, I. D., and C. T. Luis, A Model of the Hydraulic
Fracturing Processes for Elongated Vertical Fractures and
Comparison of Results with Other Models, SPE/DOE Paper No, 13864, presented at the SPE/DOE Low Permeability Gas Reservoir Symposium, Denver, May 1985.
). Advani, S. H., Finite Element Model Simulations Associated with Hydraulic Fracturing, SPE/DOE Paper No.
8941, presented at the SPE/DOE Unconventional Gae Recovery Symposium, Pittsburgh, May 1980,
1. Cleary, M. P., Kavvadaa, M., and K. Y. Lam, Development of a Fully Three-Dimensional Simulator for Analysis and Design of Hydraulic Fracturing, SPE Paper No.
11631, presented at the Symposium of Low Permeability
Gas Reservoirs, Denver, March 1983. (See also Report
No. REL-82-12, MIT Resource Extraction Laboratory $
Dec. 1982.)
1. Abou-Sayed, A. S., Sinha, K. P., and R. J. Cifton, Evalu ation of Hydraulic Fractures Using a 3-D Simulator: Part 1
Technical Approach, SPE/DOE/GRI Paper No, 12877,
and Abou-Sayed, A. S., Clifton, R. J., Dougherty, R. L., an d
R. H. Morales, Evaluation of Hydraulic Fractures Using a
3-D Simulator: Part 2 Caae Studies, SPE/DOE/GRI
Paper No. 12878, presented at the SPE/DOE/GRI Unconventional Gas Recovery Symposium, PM.sburgh, 1984.
i5u69
WE
_..
.-
L,,,
. .
3
--Q
,.
-.~~
,-hi.
\::
x.?
.
L
+%L
..O
PL
+
j-&:
RESERVOIR
1
F ,:_~yR::KE
i. ...
~.1
~
SC ~
LLAKOFF
pPENETRATIoN
TAERMAL
PEN ETRATION---
4.4
Ah
a
PR
1. L,
.
,,, ~
x,
SpE
siu
I@e,
90.00
SLURRYFLOU (SPII)
m.ue
70.00
SAND CONC
2.08
(PPG )
d2.JA+&~
@
f--l
30
Fis.
TIME
@-Flow, dwo~,
(MINUIES)
5800
UPHOLE
PRESSURE
(PSI)
4500
4s00
3500
3080
2500
LENSOR
Zeea
1sea
1000
sea
15069
~ 0[ X [Ill)
H19N31 .[13)
1H913H 3WlA3ti#d
tLL
1-1
nL
0
lx
L
(Ld)
(ISd)
3WISS3Ud SS32X3
U9N31
lwddoud