This Excel spreadsheet uses the Ei(-x) function (or its approx) to calculate p(r,t) values f

1. Given the argument x, calculate Ei(-x): Enter x argument here:

2
x= 0.006
−948 φμ c t r
[ ] [
p( r ,t )= pi +70 . 6
qB μ
kh
⋅Ei
kt ]
2
948 φμ ct r
x= [ kt ] Note:
ln(1.781x) is a good approx to Ei(x) when x ≤

2. Given reservoir and well info, calculate x, Ei(x), Dp, p(r,t), tmin, tmax, and check t.

3. Given reservoir and well info, calculate x, Ei(x), Dp, p(r,t) at different times.

Data: Value: Units:

pi 2500 psia r [ft]

q 150 STB/d 0.3

B 1.47 RB/STB 3
m
qB μ =
0.6 cp 70 .6 9.7296 60
k 16 md kh 300
h 60 ft 900
1500
f 0.2 fraction
ct 7.00E-06 1/psi
calculate p(r,t) values for a well 2 #VALUE!

Calculations: Note: The Ei(-x) result is negative and should be used as negative
Ei(-x) = #VALUE! in the equation. Sometimes it is written as positive (like in our table), but b
In the p(r,t) equation below, for a positive q, you can see that the entire
ln(1.781x) = -4.5388 negative to be subtracted from pi to get the expected p(r,t) for a flowing we
2
−948 φμ c t r
Error: #VALUE! p( r ,t )= pi +70 . 6
qB μ
kh [ ] [
⋅Ei
kt ]
od approx to Ei(x) when x ≤ 0.01.
−948 φμ c t r 2  3.975 105  ct rw2   948 c

max
, and check t.
Δp=70 . 6
qB μ
kh
⋅Ei
kt[ ] tmin  
 k


tmax  
 k

p(r ,t )= pi + Δp [ft] rw = [ft] re =

2
2 −948 φμ c t r
erent times. x=
948 φμ ct r
[ kt ] Δp=70 . 6
qB μ
kh
⋅Ei
kt[ ]
[hrs] Time = 10 [hrs] Time =
pi x Ei(-x) Dp p(r,t) x Ei(-x)

2500 4.48E-07 #VALUE! #VALUE! #VALUE! 2.24E-07 #VALUE!

2500 4.48E-05 #VALUE! #VALUE! #VALUE! 2.24E-05 #VALUE!
2500 1.79E-02 #VALUE! #VALUE! #VALUE! 8.96E-03 #VALUE!
2500 4.48E-01 #VALUE! #VALUE! #VALUE! 2.24E-01 #VALUE!
2500 4.03E+00 #VALUE! #VALUE! #VALUE! 2.02E+00 #VALUE!
2500 1.12E+01 #VALUE! #VALUE! #VALUE! 5.60E+00 #VALUE!

3000.000

2500.000

2000.000

p(r,t) 1500.000 pi
Column
N
Column
1000.000 R
Column
V
500.000

0.000
0 200 400 600 800 1000 1200 1400
Radius, r [ft]
p(r,t) 1500.000
Column
N
Column
1000.000 R
Column
V
500.000

0.000
0 200 400 600 800 1000 1200 1400
Radius, r [ft]
used as negative
ve (like in our table), but be sure to use it as negative.
can see that the entire Dp term must be
cted p(r,t) for a flowing well.
∞ n
x
Ei (−x )=−ln x−0 . 5772+ ∑ (−1 )n
n=1 n⋅n !
2 3
x x
Ei (−x )=−ln (1 . 781 x ) + x− + +. ..
 948 ct re2  4 18
tmax   
 k 

20 [hrs] Time = 60
Dp p(r,t) x Ei(-x) Dp p(r,t)

#VALUE! #VALUE! 7.47E-08 #VALUE! #VALUE! #VALUE!

#VALUE! #VALUE! 7.47E-06 #VALUE! #VALUE! #VALUE!
#VALUE! #VALUE! 2.99E-03 #VALUE! #VALUE! #VALUE!
#VALUE! #VALUE! 7.47E-02 #VALUE! #VALUE! #VALUE!
#VALUE! #VALUE! 6.72E-01 #VALUE! #VALUE! #VALUE!
#VALUE! #VALUE! 1.87E+00 #VALUE! #VALUE! #VALUE!

mn

mn

mn

1200 1400 1600

mn

mn

mn

1200 1400 1600