Pressure Buildup Analysis With Wellbore Phase Redistribution
Pressure Buildup Analysis With Wellbore Phase Redistribution
Pressure Buildup Analysis With Wellbore Phase Redistribution
fOR
ldt) lOCI) 1()()
OOOOOQ(3)
_____ L---"-"-i'-Cjl!O 10
"wo l8"" t."W lOOPS'
'I) 1]0"" .... ' 1HR
Fig. 9 - Pressure buildup data, Example 2.
with the buildup data. In gas-lifted oil wells, this
gradient must be measured below the point of gas-lift
gas entry. Any differences in the flowing and static
gradients which cannot be attributed to frictional
effects generally will give an indication of the flow of
free gas from the reservoir. The gradients used in the
following examples were measured in conjunction
with the pressure surveys.
Note also that the type curves presented are not
meant to replace semilog analysis methods or the use
of previously published type curves. If it is possible
to analyze well test data using semilog methods,
greater accuracy will be obtained in nearly all cases,
mainly due to the similarity of the shape of the type
curves which makes type-curve matching difficult.
When such simple analyses are not possible,
however, the type curves presented in this work may
permit approximate analysis which would not be
possible otherwise.
Example Analyses
Example 1 is an actual set of pressure buildup data
measured in a gas-lifted oil well in southeast
Louisiana. The basic data are shown in Table 7 and a
log-log plot of the pressure data is shown in Fig. 8.
From the data plot in Fig. 8, a point on the unit
slope straight line is estimated to be !::t.p = 153 psi
(lOSS kPa) at !::t.t = 0.1 hour. The wellbore storage
coefficient is calculated as in Eq. 30 and the apparent
storage coefficient as in Eq. 31. The gradient used in
Eq. 30 is calculated from flowing- and static-pressure
surveys measured in conjunction with the buildup
test.
Aw
C= - = 0.01173 bbllpsi
PI
(0.000 270 5 m
3
IkPa). . ............... (30)
268
TABLE 8 - PRESSURE BUILDUP DATA FOR FIELD EXAMPLE 2
q = 14BID(2.23m
3
/d)
Ii = 4 cp (5 x 10- 4 Pa s)
B = 1.05 RB/STB (1.05 res m
3
/stocktank m
3
)
h = 20 ft (6.10 m)
Aw = 0.00387 bbl/ft (0.002 02 m
2
)
<I> = 0.28
Ct = 150x10-
6
psi-
1
(2.176x10-
5
kPa-
1
)
r w 2 = 0.085 ft2 (0.007 90 m
2
)
PI = 0.420 psi/ft (9.50 kPa/m), measured
M
(hours)
o
0.25
0.50
0.75
1
2
3
4
6
8
10
12
14
16.5
C4Jw
psi (kPa)
102 (703)
190 (1310)
254 (1751)
278 (1917)
306 (2110)
292 (2013)
284 (1958)
273 (1882)
276 (1903)
276 (1903)
276 (1903)
278 (1917)
281 (1937)
qB!::t.t
C
a
= 2-- =0.00635 bbllpsi
4!::t.p
(0.000 1464 m
3
IkPa). . .............. (31)
Since C
a
< C, we can conclude that phase
redistribution effects are significant. These values
yield CD =752 and CaD =407.
The data then are matched to the type curves for
CaD =400 and CD =750 as indicated, with a match
point chosen as C<I>D = 10, S = 0, t D = 6,800 at !::t.t = 1
hour, and PwD = 1.02 at !::t.pw = 100 psi (689 kPa).
From the standard definitions of tD and PwD' the
permeability is calculated as follows.
FrompwD match: k= 134 md.
From tD match: k= 144 md.
Example 2 consists of data taken in a well in
southeast Louisiana producing at low rates and high
water cuts from a shaly sand. Pressure buildup data
is given in Table 8 and the log !::t.pw vs. log !::t.t plot is
shown in Fig. 9. From the static and flowing surveys
taken in conjunction with the pressure buildup, C is
calculated as shown and C
a
also is estimated from an
extrapolation of the short-time data.
Aw 3
C= - =0.00921 bbllpsi (0.000 212 m IkPa).
PI
qB!::t.t
C
a
= -- = 0.00130 bbllpsi
24!::t.p
(0.000 030 m
3
IkPa).
Since C
a
< C, phase redistribution effects are
believed to be significant, so CD and CaD are
calculated to be
SOCIETY OF PETROLEUM ENGINEERS JOURNAL
CD = 115 and CaD = 16.
The data are matched on. the type curves for
CD = 100, CaD = 20 and the best match is estimated
as shown in Fig. 9 to be Cq,D = 10, S = 5, PwD = 2.80
at Llp
w
= 100 psi (689 kPa), and tD = 320 at Llt= 1
hour. From the standard definitions of P wD and t D'
the permeability is calculated as follows.
FromPwD match: k= 1.45 md.
From tD match: k=2.16 md.
Fig. 9 also seems to indicate that the last two data
points may be close to the semi log straight line. Using
semilog analysis techniques, the permeability and
skin are estimated to be 1.46 md and 3.4, respec-
tively. Since the true straight line may not have been
reached and only two points are used to determine
the semi log straight line, these estimates are in
adequate agreement with the estimates obtained from
type-curve matching.
Summary and Conclusions
It has been recognized for some time that wellbore
phase redistribution can cause anomalous pressure
buildup behavior in oil and gas wells. General aspects
of the phenomenon have been presented
previously2,3; however, no technique for the analysis
of such tests has been available.
The work presented in this paper provides an
analysis of the well bore phase redistribution problem
and, with an assumed behavior based upon physical
arguments, provides a general method for the
analysis and description of such anomalous pressure
buildup tests. It has been shown that the well bore
phase redistribution problem is a complex well bore
storage phenomenon, and mathematical methods
previously applied to well bore storage problems have
been extended to solve this more general problem.
In the analysis of buildup surveys, I have found
that the observed storage constant often does not
agree with that calculated from the well completion
properties. One possible explanation for this ob-
servation lies in the apparent storage observed to be
associated with phase redistribution. Even though a
hump may not be observed, phase redistribution
effects may cause an inobvious distortion in the data
plot. Analysis of such data by other type curve
techniques may yield totally meaningless results. In
view of this, it is recommended that the true and
apparent storage coefficients always be calculated
and checked for consistency before proceeding with
detailed analysis of a buildup survey.
The main assumption of this work is the ex-
ponential form used in representing the phase
redistribution pressure function. I have found that
this form apparently represents phase redistribution
in a gas-lifted oil well very well; however, no
meaningful experimental data are available to
substantiate this completely. Such data would be
useful either in verifying this function or in
proposing a new function for the phase redistribution
pressure. This data could be collected by measuring
the pressure in a well shut in simultaneously at the
surface and at the bottom of the tubing string using
APRIL 1981
equipment described in Ref. 3. Laboratory ex-
periments also could measure this pressure. This data
and further analytical work is definitely needed to
determine the range of well conditions over which the
assumed form is applicable and to extend the basic
technique to other conditions.
In this work, only positive values of the skin effect
factor have been considered. It would pose no major
problem to calculate dimensionless well bore
pressures for negative skin factors by the technique
described by Agarwal, et at.
5
Note, however, that
such an approach places a great emphasis on the
accuracy of the various functions used at small times
and these functions are inherently more difficult to
evaluate with great preCISIOn. In the phase
redistribution problem, such an approach would
require the evaluation of the dimensionless pressures
at extremely small dimensionless times, dimen-
sionless storage coefficients, and dimensionless phase
redistribution time parameters.
Although the numerical work presented here is
based upon an infinite, homogeneous, radial
reservoir model, the basic concepts are much more
general. In particular, it is possible to apply the
techniques used in this study to other reservoir
models and thereby to obtain techniques for the
analysis of data in fractured systems as well as other
practical situations.
N omencJature
A =
w
B=
C=
cross-sectional area of the well bore,
bbl/ft (m
2
)
formation volume factors, RSB/STB
(res m
3
Istock-tank m
3
)
well bore storage coefficient, bbl/psi
(m
3
/kPa)
apparent storage coefficient, bbl/psi
(m
3
/kPa)
apparent dimensionless storage coef-
ficient,
5.6146C
a
C D= 2
a 27r(j>c thr w
= ( C q,D + _1 ) _ 1
aD CD
(CaD = p ( j > ~ ~ r w 2)
C
eD
= effective dimensionless storage coef-
ficient defined in Eq. 5
c ( compressibility, psi - 1 (kPa - 1 )
CD dimensionless well bore storage coef-
ficient,
(
CD = C )
2p(j>c(hr w 2
phase redistribution pressure parameter,
psi (kPa)
269
dimensionless phase redistribution
pressure parameter,
khC<I>
C<I>D = 141.2 qBJl
(
_ 7.27rx 1O-6khC<I
C<I>D - ---q-B-Jl--:r:...
reservoir thickness, ft (m)
reservoir permeability, md
modified Bessel function of second kind
of order n
(PD)
(PwD)
(P<I>D)
P
Pgej
Laplace transform of P D
Laplace transform of P wD
Laplace transform of P<I>D
pressure, psi (kPa)
flowing pressure at point of gas entry,
psi (kPa)
270
Pw
Pwhj
PwD
PD
P<I>
P<I>D
well bore pressure, psi (kPa)
flowing wellhead pressure, psi (kPa)
dimensionless well bore pressure
dimensionless pressure
phase redistribution pressure, psi (kPa)
dimensionless phase redistribution
pressure,
_ khP<I>
P<I>D - 141.2 qBJl
(
_ 7.27rx 1O-6khP<I
P<I>D-
qBJl
q = flow rate, BID (m
3
Id)
qsj = sand-face flow rate, BID (m
3
Id)
r = radius, ft (m)
r w well bore radius, ft (m)
r D = dimensionless radius, r D = r / r w
s = Laplace transform variable
S = skin factor
t time, hours
t D = dimensionless time,
0.OO0264kt
t
D
= 2
Jlctr w
(
_ 3.6x 1O-6kt)
tD - 2
Jlctr w
phase redistribution time parameter,
hours
dimensionless phase-redistribution time
parameter,
Jl
Pj
7 =
0.OO0264ka
aD= 2
Jlctr w
(
_ 3.6x 1O-6ka)
fXD- 2
Jlctr w
Eulers constant,
,,(=0.577 215 664901 53
fluid viscosity, cp (Pas)
fluid density, psilft (kPa/m)
dummy variable of integration
porosity (fraction)
pressure difference, psi (kPa)
shut-in time, hours
dimensionless shut-in time
Acknowledgment
The author thanks the management of Shell Oil Co.
for permission to publish this paper.
References
I. Earlougher, R.C. Jr.,: Advances in Well Test Analysis, SPE
Monograph Series, Dallas (1977) 5.
2. Stegemeier, G.L. and Matthews, C.S.: "A Study of
Anomalous Pressure Build-Up Behavior," Trans., AI ME
(1958) 213, 44-50.
3. Pitzer, S.C., Rice, J.D., and Thomas, C.E. : "A Comparison
of Theoretical Pressure Build-Up Curves with Field Curves
Obtained from Bottom-Hole Shut-In Tests," Trans., AIME,
216,416-419.
4. Van Everdingen, A.F. and Hurst, W.: "The Application of the
Laplace Transformation to Flow Problems in Reservoirs,"
Trans., AIME (1949) 186, 305-324.
5. Agarwal, R.G., Rafi AI-Hussainy, and Ramey, H.J. Jr.: "An
Investigation of Well Bore Storage and Skin Effect in Unsteady
Liquid Flow: l. Analytical Treatment," Soc. Pet. Eng. J. (Sept.
1970) 279-290; Trans., AIME, 249.
6. Stehfest, H.: "Algorithm 368 - Numerical Inversion of Laplace
Transforms," Communications oj the A eM (1970) 13, 47.
7. Handbook oj Mathematical Functions with Formulas, Graphs
and Mathematical Tables, Abramowitz, M. and Stegun, l.A.
(eds)., Nat!. Bureau of Standards Applied Mathematics Series
55, U.S. Dept. of Commerce (1972).
SI Metric Conversion Factor
psi x 6.894 757 E+OO kPa
SPEJ
Original manuscript received in the Society of Petroleum Engineers office
July 16,1979. Paper accepted for publication Feb. 21,1980. Revised manuscript
received Nov. 10, 1980. Paper (SPE 8206) first presented at the SPE 54th Annual
Technical Conference and Exhibition, held in Las Vegas, Sept. 23-26, 1979.
SOCIETY OF PETROLEUM ENGINEERS JOURNAL