Acoustic Determination of Producing

Bottomhole Pressure
dames N. McCoy, SPE, Echometer Co.
Augusto L. Podlo, SPE, U. of Texas
Ken L. Huddleston, SPE, Echometer Co.

Summary. This paper discusses the acoustic determination of producing bottornhole pressure (BHP). Two different techniques
are presented for wells that have liquid above the formation and gas flowing upward through the gaseous liquid column. One
technique involves the acoustic measurement of the liquid level and the casing-pressure buildup rate when the casinghead valve is
closed. When these data are used along with an empirically derived correlation given here, the gradient of the gaseous liquid
column in the annulus can be obtained. This technique offers a reasonably accurate procedure for determining the producing BHP
of a well by acoustic means. The second method involves two acoustic measurements. A backpressure valve is used in the casing
head to depress and to stabilize the liquid level at two positions while the well is produced at a constant rate. The gradient of the
gaseous liquid column is then calculated and extrapolated to the formation depth.
This paper discusses results from the field testing of numerous wells where the actual gradients of gaseous liquid columns were
measured in a variety of casing/tubing sizes, oil gravities, gas flow rates, and pressures.

Introduction
The producing-rate efficiency of a well can be determined with the Liquid Level Above Formation With Casinghead Gas Flow (Case
curve of inflow performance relationship, 1 which requires C). This condition results in a gaseous annular liquid column. At
knowledge of the producing and static BHP's. Techniques for de- stabilized producing conditions, the oil in the casing annulus
termination of static BHP's by acoustic means have been presented becomes saturated with the gas that is continuously flowing to the
in Refs. 2 through 5. These techniques have prO'ved to be suffi- surface. Consequently, if gas is being vented at the surface at a
ciently accurate for most conditions. constant rate, free gas is being produced from the formation si-
The producing BHP is the sum of the surface casing pressure multaneously with the oil. Generally, most oil is produced through
plus the pressure from the column of fluids in the annulus. the pump while most free gas is produced up the casing annulus.
The fluid distribution in the annulus is a function of the producing BHP calculation is undertaken from a measurement of casinghead
conditions of the particular well. Three situations are generally found pressure, knowledge of oil and gas properties, and an estimate of
in the field: (1) the liquid level is at or near the formation and the oil fraction in the annular liquid. The fraction estimate is re-
casinghead gas mayor may not produced; (2) the liquid level is quired to obtain the gradient of the gas/liquid mixture. This problem
above the formation and casinghead gas is not produced; and (3) has received considerable attention by numerous authors. 7- 11 These
the liquid level is above the formation and casinghead gas is techniques involve the determination of the gas flow rate up the
produced. annulus and, in tum, the calculation of the amount of liquid present
Fig. I illustrates these three cases. For Cases A and B, the in the gaseous liquid column by use of such well conditions as
pressure distribution is well defined from a measurement of the casing/tubing sizes, liquid properties, and pressure. All these
pressure at the surface, a knowledge of the properties of the fluids, methods are based on a combination of theoretical and empirical
and the position of the liquid level. Case C, on the other hand, in- models and yield different results for a given set of conditions as
volves the uncertainty of the gaseous liquid column gradient as a shown by Kabir and Hasan. 10
result of the annular gas flow. Because of the disagreement of these techniques, a comprehensive
field study was performed to determine directly the gradient of
Liquid Level at Formation (Case A). The casinghead pressure gaseous liquid columns.
constitutes the major portion of the producing BHP in normal-<lepth The wells tested during this study included casing sizes from 4.5
wells because the pressure from the gas column is relatively small. to 7 in. [11.4 to 18 em] and oil gravities between 32 and 43 0 API
Even when gas is being vented, the frictional pressure losses are [0.86 and 0.81 g/cm 3 ]. Long gaseous liquid columns of more than
minimal. BHP calculation is undertaken from a measurement of 5,000 ft [1525 m] were studied in wells up to 9,000 ft [2745 m]
the casinghead pressure, the knowledge of the gas composition, and deep. Annular gas flow rates ranged from 13 to 120 Mcf/D [368
the temperature distribution as described in Ref. 2. Or the BHP to 3400 m3 /d] and oil fractions ranging from 20 to 77% were
calculation can be performed by a computer program given in Ref. measured. The wells were located in regions of normal temper-
6. This program also includes Cases B and C. ature gradients in the range of 0.9 to 1.2°F/100 ft [16.4 to 21.9
The liquid level will always be at the tubing perforations when mK/m].
a well is being produced with the casing valves closed and free gas When a gaseous liquid column exists in the annulus of a well
is flowing from the formation. producing at stabilized conditions. the pressure at any depth in the
gaseous column is independent of the surface pressure. This is il-
Liquid Level Above Formation Without Free Gas Flow From lustrated in Fig. 2, which is a schematic of a well producing at three
Reservoir (Case B). At stabilized producing conditions, the liquid different values of casinghead pressure. The producing BHP remains
above the tubing perforations is 100% oil. This producing BHP unaffected by the changes in surface pressure and liquid level as
is calculated from measurement of the surface casinghead pressure, long as the production rates through the tubing and the casing an-
measurement of the depth to the liquid level by an acoustic survey, nulus remain constant.
and a knowledge of the oil and gas properties. Details of the calcu- The annular pressure for three cases is plotted as a function of
lation are given in Ref. 2. depth in Fig. 3. The gradient of the gaseous column (gas/oil mixture)
can be obtained from the change in liquid level and in pressure at
Copyright 1908' Society 01 Pelrolaum Engineers the gas/liquid interface.

SPE Ponnation Evaluation. September \988 617

 - I-
=
Pc' 100 Pel =200

10 PSI [

80 PSI

•
390 PSI

•••
TYPE A
LIQUID LEVEL
TYPE B
LIOUID LEVEL
IlJ
TYPE C
LIQUID LEVEL
20 PSI [

AT PERFORATIONS ABOVE FORMATION ABOVE FORMATION

AND NOT PRODUCING AND PRODUCING Pwt =500 pwt=500 Pwt =500
ANNULAR GAS ANNULAR GAS
Fig. 2-Pressure relations In a pumping well with a gaseous
Fig. 1-Fluld distribution In the annulus of a pumping well. fluid column.

The gaseous oil column gradient is equal to the difference in gas to vent at its original rate. When casing pressure ceases to in-
pressures at the two fluid-level depths divided by the distance be- crease, liquid from the annulus is no longer forced into the pump
tween them: and the producing BHP returns to its original value. The well is
in a stabilized condition and a true gradient can be calculated.
[(Pef2 +Pgc2)-(Pefl +Pgel)]
'ggle= , (1) During earlier tests, BHP sensors, dynamometers, and portable
(D L 2- DLI) well testers were used to determine when stabilization occurred.
As experience was gained, the results showed that a stabilized
where D L2 and D LI are the depths to the liquid level at pressures flowing condition was accurately indicated by a stabilized liquid
P cf2 and P efl' respectively, and Pgel and Pge2 are the pressures level, obviating the need for this auxiliary equipment.
from the hydrostatic gas column. Figs. 4 and 5 show the most recent field data, obtained by a
Eq. 1 suggests the procedure for field determination of the gaseous microprocessor-controlled test system 12 that allows data to be
column gradient: the liquid level can be depressed by increasing taken as frequently as once a minute, thus yielding an accurate and
the casing pressure. Then, after waiting until stabilized flow con- detailed description of the behavior of the fluids in the well. The
ditions are re-established, the new liquid level is measured and the casing pressure is plotted as a function oftime in Fig. 4. Measure-
gaseous column gradient is calculated with Eq. 1. ments were made at three stabilized values of casinghead pressure:
It will generally take considerable time (more than 24 hours) for 80 psi [551 kPa] (Points 1 and 2), 180 psi [1.2 MPa] (Points 3 and
stabilized conditions to be established when an increase in the 4), and 250 psi [1.7 MPa] (Points 5 and 6). The corresponding fluid-
casinghead pressure occurs. During this transient period, the level behavior is presented in Fig. 5, which shows that for this par-
pressure buildup at the surface (caused by the accumulation of gas ticular well, fluid-level stabilization occurs in a relatively short time,
flowing into the gas column from the gaseous liquid column) will 5 to 10 hours.
result in a pressure increase in the annulus and at the formation. The measured gradient includes both the hydrostatic (mixture
This increase causes liquid to be depressed from the casing annulus density) and the dynamic (friction and kinetic energy) components.
into the pump, thus restricting fluid flow and possibly causing liquid The linear behavior with depth shows that the total gradient remains
backflow into the formation. During this time, the liquid level in practically constant, even though the pressure varies by a factor
the annulus will drop to a lower level. After gas begins to vent at of 3 or 4 as shown in Fig. 6. This behavior was observed in all
the increased stabilized surface pressure, the surface gas flow rate the field tests covered by this study. The measured total gradient
will stabilize at the original value and the liquid level will stabilize was expressed in terms of the gas-free oil gradient by introducing
as the producing BHP returns to its original value. In some instances, a term defined as the effective oil fraction:
this process may take 1 or more days.
lo=gmtlgo • ...................................... (2)
Testing Procedure
With the use of a backpressure vll;lve, the casing pressure was in- where
creased to a specific value and then stabilized by allowing annular 10 = effective oil fraction,

Pet, Pc12 Pc13 , 215.0

=100

constant ~,
~ ~=~ 15,90 I 2
g,gg 31 63 94 125
LO--1-0,l..O---2-0.L.O--3-0.L.O--4-0..L-O--5--0~O psi
ELAPSEO TIME, HOURS
PRESSURE )
Fig. 4-Fleld recording of casinghead pressure during fluid-
Fig. 3-Pressure/depth traverse In the annulus for the cases level stabilization procedure.
Illustrated In Fig. 2.

618 SPE Formation Evaluation, September 1988

 2l1llg 425g
).I~~2
1 '
1-2
..,.; 2563 3688
u. : GRADIENT =0.0962 PSI/FT
..,.;
......r ~ lJ-

...>
~
3125 I:, .
..
~.
x 3125
" I.\tt

c I I ", "'~I\I," \' '~. 4 l.!l

3 I, ~
:5 x "
e
::i 3688
1 3-4
\ 2563
I"
5....,"...IIa\I...... "·I' "" ....... ,6
4m 2ggB
g,gg 31 63 94 125 15 131 188 a44
ELAPSED TIME, HOURS PRESSURE, PSIG

Fig. 5-Fleld recording of annular fluid-level depth corre- Fig. 6-Gaseous liqUid column height as a function of pres-
sponding to Fig. 4. sure at the gaslllquid Interface.

8 ml = measured total gradient, and to the volume of gas above the fluid level, yielding an adjusted depth
80 = gradient of gas-free oil. to the liquid level:
This term was then correlated directly with the corresponding an- D a =DL +(l-fo)L, (3)
nular gas flow rate and annular area as shown in Fig. 7. Note that
implied in the correlation is the effect of temperature and gas com- which is then used to calculate the annular gas flow rate.
pressibility factor, because it is derived with field data from regions
with normal temperature gradients and wells producing casinghead 0,00068 x ApADa
gas with specific gravity corresponding to oils in the 32- to 43 0 API Q= , (4)
[0.86- to 0.81-g/cm 3 J gravity range. . At
Direct measurement of annular gas flow in the field is a tedious where
process that can be avoided by estimating the gas flow from a short . Ap = change in pressure,
casing-pressure-buildup test. The test is conducted by closing the At = elapsed time, and
casing valve while the well continues to pump. The rate at which
casing pressure increases is measured. During the early phase of
A = cross-sectional area of annulus.
this study, it was observed that annular gas flow rates calculated To facilitate use of the effective-oil-fraction correlation, it has been
by short casing-pressure-buildup tests generally were less than the plotted with its axis representing the terms corresponding to the
rates measured with a critical flow prover. This error was a result casing-pressure-buildup rate times the adjusted depth to the liquid
of considering only the volume of gas above the liquid level and level:
neglecting the gas bubbles present in the gaseous liquid column.
The procedure was corrected as follows: once the gradient present
was determined, an effective oil fraction was calculated and, in tum,
the volume of gas present in the column was estimated and added

'g. FRACTION Of' GAS IN

GASEOUS LIQUID COLUMN

100'OOO".~0~JI~~"~'7~1"~~·!l~~·4~'~5~'2~~§
10,000
4D,DDO
50,000 1-+-+-+-+-+-+-+-+-+---1
2D,DDO I-+~-+-+-+-+-+-+-+-J

10,000
a,ooo
~m~~!~D
4,000 '\.
5,000 I-+--+--i---"ot:~---Ir-+-+--i-+--;;

l'OOOmfRm
[PSI-FTl 2,000 f-+--+--i-f--''''\oolt::-+-+---+--+-J
[MIN J ""i'..
aoo
400 ,
:~~ \ !DO
\
aoo 1-+--+-1-+---1--+-+--+--\t--t 200
200 1-+--+-+-+-+-+-+--+-;-;
,~ 100 .01
..\
100 L-..l--l----I_.L---..l---l._.L..-..L....J...»J o .2 .4 .6 .8 1.0
o .1 .2 .5 .4 • a .• .7 .• .S I,D
I" EFFECnve OIL FRACTION
CORRECTION FACTOR
Fig. 7-Echometer gaseous liquid column gradient correction
curve. Fig. 8-Comparlson of Gilbert's gradlent-eorrectlon·factor
curve With data points from this experimental stUdy.

SPE Formation Evaluation, September 1988 619

 Using the backpressure-valve technique presented here will increase
TABLE i-SAMPLE CALCULATION the accuracy of the calculated BHP and afterflow. Also, if the
gaseous liquid column is stabilized a short distance above the pump,
Depth to formation, ft 7,000
Depth to gaseous liquid column, ft 3,000 then the influence of possible errors in the effective-oil-gradient
Oil gravity, GAPI 42 correlation is diminished.
Oil gradient, psi/ft 0.353
Casing-pressure buildup, psi 10 Conclusions
Buildup time, minutes 30
Casing pressure, psig 100 Calculation of producing BHP can be performed in 90 % of pumping
Gas column pressure, psi 10 wells without consideration for gaseous liquid columns. Most wells
tip/At, psi/min 0.333 are in the first two categories described, either having no liquid
0 8 , ft 3,000+(1-1)4,000=3,000 above the formation or having liquid above the formation without
(tip/At)D 8' psi.ftImin (assume oil fraction = 1) 1,000 annular gas flow.
011 fraction (from Fig. 7) 0.72 In wells where gaseous liquid columns are present, the producing
0 8 , ft 3000+(1-0.72)4,000=4,120 BHP can be obtained by use of the correlation presented here. For
Iterate until oil fraction converges as shown below additional accuracy, the difficult, expensive, and lengthy Walker's
-E.L. (tip/At)D 8 .!.s- method may be used to obtain the producing BHP. Accuracy from
4,120 1,371 0.63 Walker's method depends on depressing and properly stabilizing
4,480 1,492 0.62 the liquid level near the pump.
4,520 1,505 0.62
This study, based on actual gradients measured in the annulus,
Gaseous column offers increased accuracy over previously reported correlations
pressure, psig = 0.62(0.353)4,000 = 875 based on theory and surface models.
Pwr, psi 100+ 10+ 875= 985
Nomenclature
A = cross-sectional"area of annulus, in. 2 [cm 2 ]
To date, the most commonly used method to estimate the gra-
dient in a gaseous liquid column is Gilbert's curve. For this reason, Da = adjusted depth to liquid level, ft [m]
the data points were plotted on Gilbert's curve and are shown in Df = depth to formation, ft [m]
Fig. 8. Gilbert's curve predicts gradients greater than those D L = depth to top of liquid, ft [m]
measured in the field tests, sometimes by a factor of two. fo = effective oil fraction
8g1c = gradient of gaseous liquid column, psilft [kPa/m]
Determining Acoustic Producing BHP's 80 = gradient of gas-free oil coluJ.lll1, psi/ft [kPa/m]
The simplest technique to calculate a producing BHP is to determine L = length of gaseous liqUid column, ft [m]
acoustically the distance to the top of the gaseous liquid column Pc] == casing pressure at producing conditions, psig [kPa]
and measure the casing-pressure-buildup rate. Then, use Fig. 7 to Pgc = pressure exerted by gas column, psi [kPa]
determine the effective oil fraction. Pglc = pressure exerted by gaseous liquid column, psi [kPal
The producing BHP will then be Pwf = producing BHP, psig [kPa]
Q = annular gas flow rate, Mcf/O [m 3 /d]
Pwf=Pef+Pge +gofoL. (5) Ilr == time, minutes
1,2 = denotes separate stabilized conditions
The gas column pressure, Pge' can be calculated by various well-
known methods that include the effect of gas gravity, z factor, and Acknowledgments
temperature. 2 ,4 Because the adjusted depth, D a , is a function of
the oil fraction and vice versa, the effective oil fraction will have Special thanks go to Dave Warner with Marathon Oil Co. and the
to be obtained with a short iterative process starting with an oil companies that have provided assistance and allowed use of their
fraction ofone as shown in the sample calculation in Table I. Note wells in these tests. These companies include Shell Oil Co., Gulf
that care must be taken in measuring the casing-pressure-buildup Oil Corp., Sun Oil Co., Bass Enterprises, R.E. Smith, Burk Royalty
rate. At a buildup rate of 1 psi [6.9 kPa] in 10 minutes, which might Co., Texaco Inc., Quinoco Petroleum Inc. and others.
be overlooked, an effective oil fraction of 0.80 could exist. A
minimum of 10 psi [69 kPa] and/or 10 minutes should be used in References
each casing-pressure buildup for the best results. I. Vogel, J.V.: "Inflow Performance Relationships for Solution-Gas Drive
The second method, known as Walker's method, 13 involves the Wells," JPT (Jan. 1968) 83-87; Trans., AIME, 243.
determination of the liquid level at two different casing pressures. 2. McCoy, J.N. er al.: "Acoustic Static Bottomhole Pressures," paper
The use of a backpressure valve is required to stabilize the well SPE 13810 presented at the 1985 SPE Production Operations Sym-
at the higher casing pressure. Be sure that the casing pressure and posium, Oklahoma City, March 10-12.
the gaseous liquid level have stabilized at both conditions. 3. Podio, A.L., Weeks, S.G., and McCoy, J.N.: "Low-Cost Wellsite De-
The gaseous liquid column gradient is calculated with Eq. 1. Using tennination of Bottomhole Pressures from Acoustic Surveys in High-
this gaseous liquid gradient, the producing BHP may be calculated as Pressure Wells," paper SPE 13254 presented at the 1984 SPE Annual
Technical Conference and Exhibition, Houston, Sept. 16-19.
4. "Calculating Subsurface Pressure Via Fluid Level Recorders," Energy
Resources Conservation Board, Calgary (1978).
5. Weeks, S.G., Podio, A.L., and McCoy, J.N.: "Fluid-Level Detenni-
The accuracy of this method increases as the depth to the second nation Through Internal Flush Tubing Without Depth, Temperature,
or Pressure Limitations," paper SPE 12912 presented at the 1984 SPE
liquid level approaches the formation.
Rocky Mountain Regional Meeting, Casper, WY, May 21-23.
Note that care must be taken to leave a small amount of liquid 6. Huddleston, K.L., Podio, A.L., and McCoy, J.N.: "Application of
above the pump to be sure that additional gas is not being forced Portable MicroComputers to Field Interpretation of Acoustic Well
into the pump. Surveys," SPE, 1985 Computer Technology Symposium, Lubbock, TX.
7. Gipson, F.W. and Swaim, H.W.: "The Beam Pumping Design Chain,"
Application to Pressure-Buildup Work Southwestern Petroleum Short Course, Texas Tech, Lubbock, TX
(1980).
To obtain the most accurate pressure data at the earliest times, the 8. Podio, A.L., Tarrillion, M,J., and Roberts, E.T.: "Bottomhole Pressure
top of the gaseous liquid level should be stabilized near, but above, Methods and Corrections are Reviewed," Oil & GasJ. (Aug. 25, 1980).
the pump before the well is shut in. Otherwise, oscillations of the 9. Godbey, 1.K. and Dimon, C.A.: "The Automatic Liquid Level Monitor
gaseous liquid column prevent obtaining accurate early-time data. for Pumping Wells," JPT (Aug. 1977) 1019-24.

620
10. Kabir, C.S. and Hasan, A.R.: "Predicting Liquid Gradient in a Pumping 51 Metric Conversion Factors
Well Annulus," SPEPE (Feb. 1988) 113-20; Trans., AIME, 285.
II. McCoy, J.N.: "Detennining Producing Bottomhole Pressures in Wells
°API 141.5/(l31.5+oAPI) g/cm 3
Having Gaseous Columns," JpT(Jan. 1978) 117-19. ft x 3.048* E-Ol m
12. Podio, A.L., McCoy, J.N., and Huddleston, K.L.: "Automatic psi x 6.894 757 E+OO kPa
Pressure-Buildup Data Acquisition and Interpretation Using a
MicroComputer-Based Acoustic Liquid Level Instrument," paper SPE •Conversion factor Is exact. SPEFE
16228 presented at the 1987 SPE Production Operations Symposium,
Oklahoma City, March 8-10. Original SPE manuscript received for review Sept. 22, 1965. Paper accepted for pUblication
July 9, 1967. Revised manuscript received Nov. 2, 1967. Paper (SPE 14254) first presented
13. Nind, T.E.W.: Principles of011 Well Production, McGraw-Hili Book at the 1965 SPE Annual Technical Conference and Exhibition held In Las Vegas, Sept.
Co. Inc., New York City (1964) 290-91. 22-25.

SPE Formation Evaluation. Septemher 19RR