2 FTP - 1-WT
2 FTP - 1-WT
2 FTP - 1-WT
Sh. Kord 1
Model
• Physical model
• Simplifying assumptions
• Mathematical model
• Choosing an appropriate element
• Governing equation
• Mass balance
• Momentum balance (Darcy’s law)
• Equation of state
• Initial and Boundary conditions
• Infinite acting
• Constant rate production
• Constant pressure production
• Finite acting
• Constant rate production
• Constant pressure production
• Solutions
• Laplace space solutions
• Time domain solutions
• Simplified solutions
• Applications (Drawdown (single rate & multi rate), Reservoir limit test, Build up,
Superposition (time & space), …),
ft psia cp 1/psia
1 ∂ ∂p ϕµ ct ∂p
r =
r ∂r ∂r 0.0002637 k ∂t
0.0002637 k
Hydraulic Diffusivity η= hr
φµ ct
md
ct = c f + co S o + cw S w
Well Testing: Flow through Porous Media 3 Sh. Kord
System of Units
For example, when the product is 0.1, the solution has about 5%
error.
r = rch rD
• Dimensionless Radius r = rw rD
t = tcht D φµct rw2
• Dimensionless Time t= tD
k
pi − p = pch pD qµB
• Dimensionless Pressure pi − p = pD
2πkh
Well Testing: Flow through Porous Media 11 Sh. Kord
Dimensionless Variables for Radial Flow
k (2πhrw ) ∂p ∂p D
Bq = rD = −1
µ ∂r rw ∂rD rD =1
The outer boundary condition with no flow boundary:
∂p ∂p D
= 0
=0
∂r r →re ∂rD reD
p (r , t = 0) = pi pD (rD , t D = 0) = 0
ft 2 cp/psi 0.0002637 k
1 hr = 0.0002637 tD = t
md rw φµct
2
1 md = 9.869233 × 10-16 m2
1 bbl = 0.1589873 m3
1 psi = 6894.757 Pa
1 ft = 0.3048 m
1 cp = 0.001 Pa s
1 day= 24×3600 s
1 hr = 3600 s
Well Testing: Flow through Porous Media 13 Sh. Kord
Initial And Boundary Conditions
Radial Flow In a Circular Reservoir
Initial Condition : p = pi , t = 0, r ≥ rw
Why
parallel?
∂2 p φµct ∂p ∂ 2 p D ∂pD
2 =
∂x ∂t
=
∂x 0.0002637 k ∂t
2
D D
x
xD =
• Dimensionless length A
0.0002637 k
• Dimensionless Time tD = t
φµct A
• Dimensionless Pressure pD =
k A
( pi − p )
2π (141.2) Bqµ
Well Testing: Flow through Porous Media 16 Sh. Kord
Dimensionless Variables for Linear Flow
∂p D
= −1
∂xD xD =1
The outer boundary condition with infinite acting
p D ( x D → ∞, t D ) = 0
pD ( xD , t D = 0) = 0
Well Testing: Flow through Porous Media 17 Sh. Kord
Solution of the Diffusivity Equation
Assumptions in the development of diffusivity equation:
• Homogeneous and isotropic reservoir
• Horizontal fluid flow with negligible effects of gravity
• Single phase fluid
• Isothermal and laminar flow
• Multiphase flow
• Horizontal wells
• Other cases
−α n t D ( )
2
qµB 2t D ∞
J12 α n reD
2 + n(reD ) − + 2∑ 2 2
3 e
pwf = pi − 141.2
kh reD
4 [
n =1 α n J 1 (α n reD ) − J 1 (α n )
2
]
Approximate Solutions
• The solution to differential equations treating the well as a vertical line through a
porous medium. The solution is nearly identical to the finite-wellbore solution. At
very early times, there is a notable difference in the solutions, but the differences
disappear soon after a typical well is opened to flow or shut in for a buildup test,
and in practice the differences are masked by wellbore storage.
∂p µ qBo ∂p µ qBo
=− =−
∂r r →0 2π rw hk ∂r r
w
2π rw hk
qµB 2t D 3
pwf = pi − 141.2 2 + n(reD ) −
kh reD 4
141.2qBµ 0.0005274 k re
pwf = pi − t + n − 0.75
φµ ct re
2
kh rw
Well Testing: Flow through Porous Media 22 Sh. Kord
Infinite Cylindrical Reservoir with Line-Source Well
Assumptions:
Well produces at constant rate
Well has zero radius
Reservoir is at uniform pressure
Well drains an infinite reservoir
1 φµ ct r 2
pD = − Ei − 948
2 k t
qBµ φµ c r 2
pwf = pi + 70.6 Ei − 948
t w
kh k t
The assumption of
zero wellbore
5
φµ c r 2
φµ c r 2
The reservoir is
limits the accuracy 3.79 ×10
t w
< t < 948 t e
no longer
k k
of the solution infinite acting
qBµ φµ ct rw2
pwf = pi + 70.6 Ei − 948
kh k t
∞
e −u ≈ n(1.781x ) for x < 0.02 (error ≈ 0.6%)
− Ei (− x ) = ∫ du =
x
u ≈ 0 for x > 10.9
qBµ 1.781φµ ct rw
2
pwf = pi + ln
4πkh 4k t
162.6qBo µ
pwf (t ) = pi − k
log t − 3.23
kh φµ c r 2
t w
162.6qBo µ
pi − pwf (t ) = log(t ) + log
k − 3.23
2
kh φµ c r
t w
Well Testing: Flow through Porous Media 24 Sh. Kord
Drawdown Test
162.6qBo µ k
pi − pwf (t ) = log(t ) + log − 3.23
kh φµ c r 2
t w
pwf (t ) = a + m log(t )
162.6qBo µ
m=−
kh
162.6qBo µ k
a = pi − log − 3.23
2
kh φµ ct rw
For: φµ ct re2
t > 948
k
qµB 2t D 3
pwf = pi − 141.2 2 + n(reD ) −
kh reD 4
141.2qBµ 0.0005274 k re 3
pwf = pi − t + n −
φµ ct re
2
kh rw 4
Well Testing: Flow through Porous Media 26 Sh. Kord
Pseudo Steady-state Flow Approximate Solution
∂pwf 0.0744qBo
=−
∂t φ ct re2
V p = π re2 hφ
∂pwf 0.234qBo
=−
∂t ctV p
During this time period, the rate of pressure decline is inversely proportional to
the liquid-filled pore volume, Vp. This result leads to a form of well testing
sometimes called reservoir-limits testing, which seeks to determine reservoir
size from the rate of pressure decline in a wellbore with time.
141.2qBµ 0.0005274 k re 3
pwf = pi − t + n −
φµ ct re
2
kh rw 4
141.2qBµ 0.0005274 k re 3
pwf = pi − t + n −
φµ ct re
2
kh rw 4
_
∆V 5.615 qB(t 24 ) 0.0744 qBt
pi − p = = =
ctV 2
(
ct π re hφ )
φ ct hre2
_
141.2qBµ re 3
pwf =p − n − + S
kh rw 4
Well Testing: Flow through Porous Media 30 Sh. Kord
Pseudo Steady-state Flow Approximate Solution
_
141.2qBµ re 3
pwf =p − n − + S
kh rw 4
141.2qBµ re
pwf = pi − n
kh rw
Well Testing: Flow through Porous Media 32 Sh. Kord
Transient linear flow, constant rate production from a hydraulically
fractured well
∂ 2 p D ∂pD
=
∂x ∂t
2
D D
p D ( x D → ∞, t D ) = 0
∂p D
The inner boundary condition: = −1
∂xD xD =1
Well Testing: Flow through Porous Media 33 Sh. Kord
Transient linear flow, constant rate production from a hydraulically
fractured well
xD2 xD
pD ( xD , t D ) = 2
tD
exp − − xD erfc −
π 4t D 2 t
D
For special case of xD=0:
qB µt
pwf = pi − 4.064
hL f kφct