Ramey PDF
Ramey PDF
Ramey PDF
Transmission Revisited
Jacques Hagoort, Hagoort & Assocs. B.V.
Summary
In this work we assess Rameys classic method for the calculation
of temperatures in injection and production wells. We show that
Rameys method is an excellent approximation, except for an early
transient period in which the calculated temperatures are significantly overestimated. We present a simple graphical correlation to
estimate the length of this early transient period. We further demonstrate that Rameys procedure for the estimation of overall heattransfer coefficients holds good only for large values of the Fourier
dimensionless time. Finally, we illustrate the results of this work
by an example calculation of wellhead temperature in a flowing
oil well.
Introduction
Estimation of temperatures in a wellbore during injection or production is a recurring problem in petroleum engineering. Examples
are the prediction of bottomhole temperatures of injection fluids
(water, gas, and steam) and of wellhead temperatures in oil and gas
wells. Almost all practical methods for the calculation of temperatures in wellbores go back to the classic paper by Ramey1 on
wellbore heat transmission published in the early 1960s. In that
paper, Ramey presented a simple analytical equation for wellbore
temperatures based on a grossly simplified heat balance. Apart
from this analytical temperature equation, Ramey also proposed a
simple procedure to estimate an overall heat-transfer coefficient
for wellbore heat losses comprising both transient heat resistance
in the formation and near-wellbore heat resistance.2 Despite its
simplifications, in practice Rameys approach seems to work remarkably well.
In 1990, Wu and Pruess3 presented an analytical solution for
wellbore heat transmission in a layered formation with different
thermal properties without introducing the simplifying assumptions of Ramey. From their example calculations, they observed
that the Ramey method is valid at long times but can generate large
errors at early times. However, quantification of the conditions
under which Rameys method could be applied was beyond the
scope of their paper.
The objective of this work is to assess the Ramey method and
to establish criteria for its applicability. To this end we have first
performed an inspection analysis of the basic wellbore heat transmission equations without the simplifying Ramey assumptions.
This entails formulation of the equations in dimensionless form
and identification of the governing dimensionless numbers. Subsequently, we have developed a rigorous solution of these dimensionless equations and have compared this solution with the
Ramey solution for various ranges of the dimensionless numbers.
We have also checked Rameys procedure for estimating overall
heat-transfer coefficients. To illustrate the results of the paper, we
conclude with an example calculation of the estimation of wellhead temperature in a flowing oil well.
Our focus is primarily on heat transmission between the wellbore tubing and the formation. We have therefore kept the heat
transmission within the tubing as simple as possible by assuming
a constant fluid density and heat capacity, no heat conduction in
the flow direction, and no frictional heating. In addition, we have
NRA
TD
quDNGztD ;UD d. . . . . . . . . . . . . . . . . . . (1)
TfD
TfD
1
rD
=
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4)
rD rD
rD
tDFo
The initial condition is
tDFo = 0, rD 1:
TfD = 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (5)
TfD
rD
and
tDFo 0, rD = :
Tf D = 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (7)
TfD
rD
rD=1
. . . . . . . . . . . . . . . . . . . . . . . . . . . . (8)
= UD1 TfDrD = 1
The unit heat-loss function is a function of Fourier time and includes the dimensionless near-wellbore heat-transfer coefficient as
a parameter.
Inspection of the previous equations shows that the temperature
distribution in the tubing as a function of time is governed by three
dimensionless groups: NGz, NRa, and UD. The number NGz is
known as the Graetz number and arises from the two different time
scales that are at work in this problem: the Fourier time scale for
heat conduction in the formation and the time scale for convective
heat transport in the tubing. It is defined as the ratio of the dimensionless Fourier conduction time and the dimensionless time for
convection in the tubing:
NGz =
tDFo L
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (9)
=
tD rcf2 u
The number NRa expresses the ratio of the heat lost by conduction
in the formation to the heat transported up the tubing by convection and is defined by
2f L
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (10)
NRa =
cw
For small values of NRa (high production rates), the heat transport
is dominated by convection and little heat will be lost to the formation. By the same token, large values of NRa point to significant
cooling of the fluid in the tubing. In honor of Ramey, we will call
the number NRa the Ramey number.
466
Urcf
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (11)
f
1 expNRaquDxD
. . . . . . . . . . . . . . . . . . . . . . (13)
NRaquD
NRaquD 2 NRaquD2 3
xD +
xD . . . . . . . . . . . . . . . . . . . . (14)
2!
3!
tion calls for the solution of the radial Fourier equation, subject to
the radiation boundary condition. We have solved this equation by
means of Laplace transformation in combination with Stehfests
numerical Laplace inverter.5 (See Appendix B.) Rameys approximation of the unit heat-loss function will be discussed in a
later section.
Fig. 1 illustrates the general features of the Ramey solution.
Here we have plotted the temperature distribution in the tubing for
NRa1 and 1/UD0 (zero thermal skin) at three successive Fourier times: 1, 10, and 100. The profile is concave on the left,
reflecting the increase in heat losses at shallower depths. Because
of continuous heating of the formation by the hot flowing fluid in
the tubing, heat losses become progressively smaller with time.
Therefore, the temperature profile moves to the right, closer to the
inlet temperature of unity.
Fig. 2 shows the effect of the Ramey number on the time
profile of the wellhead temperature. It displays the profile for
NRa0.25, 1, and 4 and 1/UD0. As expected, the larger the
value of NRa, the stronger are the heat losses and thus the lower is
the wellhead temperature. We further see that the wellhead temperature increases gradually, a consequence of the decreasing heat
losses. After an infinite time, the wellhead temperature in all cases
must be equal to the inlet temperature of unity.
Fig. 3 illustrates the effect of skin factor on wellhead temperature. Here we have plotted wellhead temperature as a function of
time for NRa1 and for 1/UD0, 1, and 4. From Fig. 3 we see that
the larger the skin factor, the lower the heat losses are and thus the
closer the wellhead temperature is to the inlet temperature. In
addition, we observe that the time profile of the wellhead temperature is flatter for larger skins. This is a direct consequence of the
flatter time profile of the unit heat-loss function for larger skins. In
the extreme case of heat losses that are completely dominated by
near-wellbore heat resistance, the unit heat-loss function is constant and so is the wellhead temperature.
Rigorous Solution
We have rigorously solved the tubing heat balance (Eq. 1) and the
associated Fourier equation (Eq. 4) by using Laplace transformation. (See Appendix C.) In Laplace space, both equations can be
solved analytically and can be combined into a single analytical
equation for the Laplace-transformed temperature as a function of
distance. We have inverted this analytical solution to time by
means of the numerical inversion algorithm of Stehfest.5 Before
discussing the general character of the rigorous solution, we first
present its behavior at early and late times, and for large thermal
skins. Under these conditions, the Laplace-transformed temperature can be approximated to a form that can be analytically converted to real time. These solutions provide helpful insight into the
temperature distribution under limiting conditions. They also serve
as benchmarks for the Stehfest numerical-inversion algorithm.
The early-time solution for the temperature ahead of the fluid
front (xD>tD) is given by (see Appendix C)
TD = 1 xD + tD. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (15)
Early means a dimensionless time much smaller than unity and
thus long before the fluid front has reached the wellhead. Eq. 15
reflects the displacement of the initial fluid with the linear temperature profile from the tubing; the profile ahead of the front
retains its linear shape but moves up linearly proportional to time.
The late-time approximation of the temperature profile reads
(see Appendix C):
TD = 1
NRaquD 2
xD. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (16)
2!
Fig. 4Wellhead temperature vs. time for constant unit heatloss function.
of a constant temperature all the way from time zero. Fig. 4 depicts
the analytical and the approximate Ramey solution for NGz0.5
and NRaquD0.8.
Fig. 5 illustrates the general characteristics of the rigorous solution vis--vis the Ramey solution. Here we have plotted the
rigorously calculated wellhead temperature vs. Fourier time for
NRa1, 1/UD0, and NGz1, along with the Ramey solution
(Eq. 13) and the early- and late-time approximations of the rigorous solution (Eqs. 15 and 16, respectively). The rigorous solution
consists of three distinct time regimes: an early transient period, a
late-time period, and an intermediate period in between. In the
early transient period, the Ramey method underestimates heat
losses and thus overestimates the wellhead temperature. This we
have already observed in Fig. 4. The arrival of the front at the
wellhead (tDFoNGz1) is not attended by a clear kink in the
wellhead temperature as we saw in Fig. 4; however, it remains a
characteristic point. The late-time regime is the period with temperatures in excess of 0.9, in which there is always excellent agreement between the Ramey and the rigorous solution. In this example, the late-time regime begins at a Fourier time of approximately 1,500. The point at which the Ramey and the rigorous
solution converge marks the beginning of the intermediate-time
regime. This occurs, for the conditions of this case, at approximately a Fourier time of 10, one log cycle after the production of
one tubing volume.
Fig. 6 shows the effect of the Graetz number on wellhead
temperature. It depicts the wellhead temperature as a function of
time for NRa1, 1/UD0, and for NGz0.1, 1, 10 along with the
approximate Ramey solution, which is the same for all values of
NGz. As we can see, the length of the initial transient period de-
pends on NGz; the larger the NGz, the longer it takes for the two
solutions to converge and thus the longer is the early transient
period. This is in line with the observation in Fig. 5 that the
intermediate regime must begin after the arrival of the fluid front
at the wellhead, which occurs at tDFoNGz.
Fig. 7 shows how Ramey number affects wellhead temperature.
Here we have plotted wellhead temperature as a function of time
for NGz1, 1/UD0, and for NRa0.25, 1, 2, and 4, again compared with the approximate Ramey solution for each value of NRa.
Smaller Ramey numbers mean lower heat losses to the formation
and thus higher wellhead temperatures. We further observe that the
smaller the Ramey number, the shorter the early transient period.
In the case of the smallest Ramey number of 0.25, the intermediate-time regime is all but absent and the late-time regime directly
succeeds the early transient period.
Fig. 8 illustrates the effect of thermal skin factor on wellhead
temperature. It shows the time profile of the wellhead temperature
for NGz1, NRa1, and for 1/UD0, 1, and 4, along with the
approximate Ramey wellhead temperature for each value of 1/UD.
The trend is as expected: larger skins mean lower heat losses and
thus higher wellhead temperatures. The larger the thermal skin, the
shorter is the early transient period. For the largest skin of 4, the
early transient period is just above the minimum value of a Fourier
time of unity, and the intermediate-time regime is absent again.
This agrees with the behavior of the wellhead temperature for very
large skins as depicted in Fig. 4, where the early transient period
ends at tDFoNGz1.
In summary, comparison of the Ramey solution with the rigorous solution demonstrates that the Ramey solution is an excellent approximation, except at early times. In the early period, the
two Ramey assumptions are clearly not satisfied and cause significant errors in the calculated temperatures.
468
The vertical axis represents, on a log scale, the ratio of the length
of the initial period tDFoi to NGz, the arrival time of the fluid front.
The horizontal axis shows the wellhead temperature calculated by
the Ramey formula (Eq. 13) at a Fourier time of NGz. We have
determined the length of the initial period for three thermal skin
factors 1/UD0, 1, and 4, and for 0.1<NGz<10 and 0.25<NRa<4.
As a convergence criterion, we have used a difference between the
Ramey solution and the rigorous solution of 1%.
As we can see, the numerical data points in Fig. 9 cluster into
three groups, depending on the thermal skin. All clusters converge
to the ratio of unity at wellhead temperatures close to the inlet
temperature of unity. At these high temperatures, the duration of
the early transient period must be close to NGz. We further observe
that the larger the thermal skin, the shorter the initial transient
period. This stands to reason because in the limit of a very large
skin, the heat-transfer coefficient becomes independent of time,
and as a consequence, the transient period is equal to a Fourier
time of NGz. (See also Fig. 4).
Unit Heat-Loss Function
Application of the Ramey equation for the tubing temperatures
requires as input the unit heat-loss function at a certain time for a
given thermal skin factor. In the previous example calculations, we
have determined this function by rigorously solving the radial
Fourier heat-conduction equation. As a shortcut, Ramey1 proposed
to estimate the heat-loss function (or overall heat-transfer coefficient) by simply adding up the transient heat resistance for conduction in the formation for zero thermal skin and the total thermal
skin. In our dimensionless notation, the unit heat-loss function is
the reciprocal of the total dimensionless heat resistance. In formula
form, Rameys approximation thus reads
1
1
1
= 0
+
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . (17)
quDtDFo quD
tDFo UD
where the superscript zero denotes zero thermal skin. The total
thermal skin 1/UD is equal to the summation of the reciprocal of
the dimensionless heat-transfer coefficients of the individual resistance zones7:
1
=
UD
U
j=1
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (18)
jD
tDFo 105 :
1
0
quD
tDFo
= 0.5ln
4tDFo
. . . . . . . . . . . . . . . . . . . . . . (24)
dTwhD
quD =
dquD
1 TwhDquDNRa + 1
quD
. . . . . . . . . . . . . . . . . . . . . . . . . . . (21)
quD
0.79526
4
tDFo
+ 0.5
tDFo
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (22)
1 tDFo 105 :
1
0
quD
tDFo
= 0.0012ln tDFo3
T1
u
U
Ui
w
x
almost 0.4%, in agreement with Eq. 21. For the Fourier time of this
example, Eq. 23 is more accurate than Eq. 24.
Conclusions
1. Inspection analysis of the Ramey wellbore heat transmission
equations shows that temperature distribution in a wellbore is
governed by three dimensionless numbers: Graetz number,
Ramey number, and thermal skin factor.
2. Rameys approximate analytical solution for the temperature
distribution depends only on Ramey number and thermal skin
factor.
3. Comparison of Rameys analytical solution with a rigorous solution shows that the Ramey solution is an excellent approximation except for an early transient period. In this early period,
the Ramey solution significantly overestimates the temperatures.
4. The duration of the early transient period can be estimated from
a simple graphical correlation based on the calculated Ramey
wellhead temperature at the moment the tubing fluid has been
completely displaced.
5. Rameys method for estimating the overall heat transfer coefficient holds good only for large values of the Fourier time.
Nomenclature
a constant
b constant
c heat capacity, fluid
cf heat capacity, formation
h heat-transfer coefficient
j index of near-wellbore heat resistance zone
L tubing length
n number of near-wellbore heat resistance zones
NGz Graetz number
NRa Ramey number
q heat loss rate per unit of tubing length
qu unit heat-loss function
r radial distance from center of tubing
rcf radius of formation inner boundary (at interface of
cement zone and formation)
ri inner radius, circular cylinder
ro outer radius, circular cylinder
rt inner radius, tubing
s Laplace variable
t time
tD dimensionless time for convection
tDFo dimensionless time for conduction (Fourier time)
tDFoi duration of initial transient period (Fourier time)
T temperature in tubing
Tbot temperature at bottom of tubing
Tf temperature in formation
Ttop initial temperature at top of tubing
Twh wellhead temperature
T0 initial temperature of formation
December 2004 SPE Journal
Subscripts
D dimensionless quantity
F formation
J index of near-wellbore heat resistance zone
Superscript
* Laplace transform
0 zero thermal skin
References
1. Ramey Jr., H.J.: Wellbore Heat Transmission, JPT (April 1962) 427;
Trans., AIME, 225.
2. Willhite, G.P.: Overall Heat Transfer Coefficients in Steam and Hot
Water Injection Wells, JPT (May 1967) 607.
3. Wu, Y.-S. and Pruess, K.: An Analytical Solution for Wellbore Heat
Transmission in Layered Formations, SPERE (November 1990) 531.
4. Carlslaw, H.S. and Jaeger, J.C.: Conduction of Heat in Solids, second
edition, Oxford U. Press (1959) 184.
5. Stehfest, H.: Numerical Inversion of Laplace Transforms, Communications of the ACM (January 1970) 13, No. 1, 47.
6. Beek, W.J. and Muttzall, K.M.K.: Transport Phenomena, John Wiley
& Sons Ltd., Chichester, U.K. (1975) 176.
7. Earlougher Jr., R.C.: Advances in Well Test Analysis, SPE Monograph
Series, Vol. 5, Dallas, Texas (1977) 38.
8. Brill, J.P. and Mukherjee, H.: Multiphase Flow in Wells, SPE Monograph Series, Vol. 17, Richardson, Texas (1999).
9. Lauwerier, H.A.: The Transport of Heat in an Oil Layer Caused by the
Injection of Hot Fluid, Appl. Sci. Res. (1955) 5, Section A, 145.
10. Marx, J.W. and Langenheim, R.H.: Reservoir Heating by Hot Fluid
Injection, Trans., AIME (1959) 216, 312.
qT,x,t
T
T
+ cu
=
. . . . . . . . . . . . . . . . . . . . . . . . . . (A1)
t
x
rt2
Here, q is the heat-loss rate from the tubing to the formation per
unit length of tubing. In general, q depends on temperature, position, and time, and is determined by heat conduction in the formation and by near-wellbore heat resistance. The latter includes
heat transfer at the tubing wall, heat transfer by conduction through
the tubing wall, heat transfer in the annular space between tubing
and casing, heat transfer by conduction through the casing wall,
and finally, heat transfer through the cement zone.2
In the case of a constant and uniform temperature in the tubing,
the heat-loss rate is described by a unit heat-loss function, which
is defined as the heat loss per unit length and per unit drop in
temperature. This unit function follows from a solution of the
formation heat balance, which for the physical model adopted
reduces to the radial Fourier heat-conduction equation. Let us, for
the time being, assume that the temperature distribution for a unit
change in the temperature of the tubing at time zero is known. The
unit heat-loss function is then given by
qu =
q
2rcf f dT
=
T1 T0
T1 T0 dr
, . . . . . . . . . . . . . . . . . . . . (A-2)
r=rcf
471
q=
T
q t d . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-3)
u
T
T
+ cu
=
t
x
rt2
Likewise, we can write for the Fourier heat conduction (Eq. A-7),
the initial condition (Eq. A-8), and the boundary conditions (Eqs.
A-9 and A-10), respectively:
TfD
TfD
1
r
=
, . . . . . . . . . . . . . . . . . . . . . . . . . . (A-15)
rD rD D rD
tDFo
tDFo = 0, rD 1:
tDFo 0, rD = 1:
. . . . . . . . . . . . . . . . . . . . . . . (A-4)
x
. . . . . . . . . . . . . . . . . . . . . . . . (A-5)
L
quDtDFo;UD =
t 0, x = 0 : T = Tbot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-6)
The unit function in Eq. A-4 can be derived from the solution of
the radial Fourier heat-conduction equation for the formation surrounding the wellbore:
1
Tf
1 Tf
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-7)
r
=
r r
r
t
where f /(f cf).
Initially, the temperature is constant throughout the formation,
thus:
t = 0 : Tf = T0 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-8)
For the inner boundary condition at the interface between cement
zone and formation, we have
Tf
= UT1 Tf, . . . . . . . . . . . . . . . . . (A-9)
r
T Ttop
Tf T0
x
r
,T =
,x = ,r = ,
Tbot Ttop f D T1 T0 D L D rcf
u
1
t, tDFo = 2 t, and quDtDFo =
q .
L
2f u
rcf
. . . . . . . . . . . . . . . . . . . . . . . .(A-11)
The tubing heat balance (Eq. A-4) and the associated initial (Eq.
A-5) and boundary (Eq. A-6) conditions then become, respectively:
tD = 0 :
TD = 1 xD , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-13)
tD 0, xD = 0:
472
TfD
rD
rD=1
NRa =
NGz =
2L
curt2
2L
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-20)
cw
tDFo L
=
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-21)
tD urcf2
and
UD =
Urcf
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-22)
f
1 d
dT*D
rD
= sT*D , . . . . . . . . . . . . . . . . . . . . . . . . . . . . (B-1)
rD drD
drD
rD = 1 :
1
dT*D
= UD T*D , . . . . . . . . . . . . . . . . . . . . . . . (B-2)
drD
s
rD = : T*D = 0 , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (B-3)
where s is the Laplace variable and the asterisk denotes Laplace
transform.
Eq. B-1 is known as the modified Bessel equation of order zero.
Its general solution reads:
T*D = aI0rDs + bK0rDs , . . . . . . . . . . . . . . . . . . . . . . (B-4)
where I0 and K0 are the modified zero-order Bessel functions of
the first and second kind, respectively, and a and b are constants
that follow from the boundary conditions. The boundary condition
in Eq. B-3 requires that a must be zero. The constant b follows
from the radiation boundary condition (Eq. B-2) and is given by
UD
sK1ss + UDK0s
. . . . . . . . . . . . . . . . . . . (B-5)
TD()
quDNGztD d , . . . . . . . . . . . . . . . . . . (A-12)
NRa
b=
TD TD
+
=
tD xD
tD
= UD1 Tf D, . . . . . . . . . . . . . . . (A-17)
t 0, r = :
rD
TfD
tDFo 0, rD = : Tf D = 0 . . . . . . . . . . . . . . . . . . . . . . . . . (A-18)
T
q t d
u
t 0, r = rcf :
TfD = 0 , . . . . . . . . . . . . . . . . . . . . . . . . . (A-16)
TD = 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-14)
q*uD =
T*D
rD
rD=1
UDK1s
sK1ss + UDK0s
. . . . . . . . . . . . . . . . . . . . . (B-6)
Back transformation of Eq. B-6 with the help of the Stehfest inverter yields the unit heat loss as a function of dimensionless
Fourier time. We have used an implementation of the Stehfest
algorithm with 12 Stehfest terms. The algorithm very closely reproduces the analytical early- and long-time approximations of the
constant terminal temperature solution (zero thermal skin).
Appendix CRigorous Solution of the Tubing
and Formation Heat-Balance Equations
We have solved the tubing and formation heat-balance equations
by means of Laplace transformation. The subsidiary equations of
the tubing heat-balance equation (Eq. A-12) and its boundary condition (Eq. A-14) are
dT*D
sT*D 1 xD +
=
dxD
quDc
* =
UD
1
. . . . . . . . . . . . . . . . . . . . . . . (C-9)
NGz s NGz + UDs NGz
From Eq. C-9, we see that the Laplace transform of the unit heatloss function becomes negligibly small for large s. The Laplace
transformed temperature (Eq. C-6) then simplifies to
T*D =
1
2
1 exp sxD +
NRasT*D 1 xDquDc
* , . . . . . . . . . . . . . . . . . . . . . . . (C-1)
1
xD = 0 : T*D = , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (C-2)
s
T*DxD =
1
2
s NRaquDc
*
* sxD +
1 exp NRaquDc
where
* = LquDNGztD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (C-3)
quDc
Rearranging Eq. C-1 leads to the ordinary differential equation
dT*D
= NRaquDc
* + 1sT*D 1 xD . . . . . . . . . . . . . . . . . (C-4)
dxD
Integration of Eq. C-4 subject to the boundary condition (Eq. C-2)
yields
lnsNRaquDc
* + 1sT*D 1 + xD + 1
= sNRaquDc
* + 1xD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (C-5)
1
s NRaquDc
* + 1
* + 1sxD +
1 exp NRaquDc
1 xD
. . . . . . . . . . . . . . . . . (C-6)
s
Thus,
quDc
* =
T*DxD =
1 1
= NRaquDc
* xD2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (C-14)
s 2!
Back transformation yields
TDxD,tDFo = 1
1
N q t x2 . . . . . . . . . . . . . . . . . . . (C-15)
2! Ra uD DFo D
Eqs. C-6 and C-8 constitute the solution of the temperature distribution in Laplace space. Back transformation by the Stehfest inverter yields the temperature distribution in real space. We have
used the same implementation of the Stehfest algorithm as mentioned in Appendix B. We have checked this algorithm against the
analytical solution for the special case of a constant unit heat-loss
function. The algorithm does fairly well in calculating the linearly
increasing wellhead temperature during the transient period and
very well in the period thereafter when the wellhead temperature
remains constant. The difference between the constant temperature
and the calculated temperature is less than 1% for Fourier times
larger than 1.3 times the Graetz number.
For small times and thus large s, the ratio K1/K0 approaches
unity and the Laplace transformed unit function then becomes
December 2004 SPE Journal
1
N q * sx 2 + ..
2! Ra uD D
TD TD
+
NRaquDTD TD0 , . . . . . . . . . . . . . . . . . . . (D-1)
tD xD
1
1
q*
=
NGz uD ss NGz NGz s NGz
UDK1s NGz
1
1 xD
+ 2
s
s NRaquD
*
. . . . . . . . . . . . . . . . . . . . . . . .(C-13)
T*DxD =
1 xD
. . . . . . . . . . . . . . . . . . . . . (C-12)
s
For small s, the product of s and the Laplace transform of the unit
heat-loss function is also small. Expansion of the exponent in Eq.
C-12 gives
1 1 NRaquD
* sxD +
1 xD
. . . . . . . . . . . . . . . . . . . . (C-10)
s
TD0 = 1 xD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (D-2)
Eq. D-1 can be conveniently solved by Laplace transformation.
The subsidiary equation of Eq. D-1 is
sT*D TD0 +
dT*D
=
dxD
1
NRaqDT*D + NRaqDTD0 . . . . . . . . . . . . . . . . . . . . . . . (D-3)
s
The boundary condition (Eq. A-14) in Laplace space becomes
1
xD = 0 : T*D = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (D-4)
s
Rearranging Eq. D-3 leads to the ordinary differential equation
473
dT*D
= NRaqD + sT*D +
dxD
TD0 1 + NRaqD
TD = 1 xD +
1
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (D-5)
s
1 exp NRaquDtD
, . . . . . . . . . . . . . . . . . . . (D-9)
NRaquD
1 exp NRaquDxD
. . . . . . . . . . . . . . . . . . . (D-10)
NRaquD
The profile (Eq. D-10) is the same as the Ramey profile (Eq. 13).
NRaq*
uD + s xD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (D-6)
Solving Eq. D-6 for the Laplace transformed temperature, we obtain
1
T*D =
sNRaquD + s
1 xD
. . . . . . . . . . . . . . . . . . (D-7)
s
474
in. 2.54*
E+00
lbm 4.535 924
E01
m3
kJ
m
C
cm
kg