The Basic Differential Equation For Radial Flow in A Porous Medium
The Basic Differential Equation For Radial Flow in A Porous Medium
The Basic Differential Equation For Radial Flow in A Porous Medium
POROUS MEDIUM
The basic equation will be derived in radial form to simulate the flow of fluids in the
vicinity of a well (Figure 1).
Figure 1. Diagram of a circular reservoir with a well at the center (Ezekwe, 2010)
Simplifying assumptions
Consider the flow of fluid through a volume element of thickness “dr” situated at a
distance “r” from the center of the wellbore radius “rw” (Figure 2).
1
Figure 2. Radial flow of a single-phase fluid in the vicinity of a producing well (Dake, 1978)
Mass flow rate IN – Mass flow rate OUT = Rate of change of mass in the volume element
𝜕𝜌
𝜌𝑞]𝑟+𝑑𝑟 − 𝜌𝑞]𝑟 = 2𝜋𝑟ℎ∅𝑑𝑟 (Eq. 1)
𝜕𝑡
𝜕(𝜌𝑞) 𝜕𝜌
𝜌𝑞]𝑟 + 𝑑𝑟 − 𝜌𝑞]𝑟 = 2𝜋𝑟ℎ∅𝑑𝑟 (Eq. 2)
𝜕𝑟 𝜕𝑡
𝜕(𝜌𝑞) 𝜕𝜌
= 2𝜋𝑟ℎ∅ (Eq. 3)
𝜕𝑟 𝜕𝑡
2𝜋𝑟𝑘ℎ 𝜕𝑃
𝑞= (Eq. 4)
𝜇 𝜕𝑟
2
𝜕 2𝜋𝑟𝑘ℎ 𝜕𝑃 𝜕𝜌
( 𝜌 ) = 2𝜋𝑟ℎ∅ 𝜕𝑡 (Eq. 5)
𝜕𝑟 𝜇 𝜕𝑟
or
1 𝜕 𝑘𝜌 𝜕𝑃 𝜕𝜌
( 𝑟 𝜕𝑟 ) = ∅ 𝜕𝑡
𝑟 𝜕𝑟 𝜇
(Eq. 6)
We have the time derivative of the density at the RHS of equation, which can be
expressed in terms of time derivative of pressure by using isothermal coefficient of
compressibility.
1 𝜕𝑉 1 𝜕𝜌
𝐶=− ] = ] (Eq. 7)
𝑉 𝜕𝑃 𝑇 𝜌 𝜕𝑃 𝑇
This is the basic partial differential equation for the radial flow of any single-phase
fluid in a porous medium. It is a non-linear equation since the coefficients of the
𝑘𝜌
both sides ( and ∅𝑐𝜌) include pressure dependent terms C, µ and . In that respect,
𝜇
3
LINEARIZATION
Expand the left hand side of Equation 9 by using the chain rule of differentiation:
1 𝜕 𝑘 𝜕𝑃 𝑘 𝜕𝜌 𝜕𝑃 𝑘𝜌 𝜕𝑃 𝑘𝜌 𝜕2 𝑃 𝜕𝑃
⌈ ( ) 𝜌𝑟 𝜕𝑟 + 𝜇 𝜕𝑟 𝑟 𝜕𝑟 + + 𝑟 ⌉ = ∅𝑐𝜌 (Eq. 10)
𝑟 𝜕𝑟 𝜇 𝜇 𝜕𝑟 𝜇 𝜕𝑟 2 𝜕𝑡
1 𝜕 𝑘 𝜕𝑃 𝑘 𝜕𝑃 2 𝑘𝜌 𝜕𝑃 𝑘𝜌 𝜕2 𝑃 𝜕𝑃
⌈ ( ) 𝜌𝑟 𝜕𝑟 + 𝜇 𝑐𝜌𝑟 ( 𝜕𝑟 ) +
𝑟 𝜕𝑟 𝜇 𝜇 𝜕𝑟
+
𝜇
𝑟
𝜕𝑟 2
⌉ = ∅𝑐𝜌 𝜕𝑡 (Eq. 12)
These two assumptions eliminate the first two terms in the left hand side of Equation
12.
𝜕2 𝑃 1 𝜕𝑃 ∅𝜇𝑐 𝜕𝑃
+ = (Eq. 13)
𝜕𝑟 2 𝑟 𝜕𝑟 𝑘 𝜕𝑡
or
1 𝜕 𝜕𝑃 ∅𝜇𝑐 𝜕𝑃
𝑟 𝜕𝑟
(𝑟 𝜕𝑟 ) = 𝑘 𝜕𝑡
(Eq. 14)
∅𝜇𝑐
= 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝑘
4
For the flow of liquids, the assumptions that were made are quite reasonable and
have been frequently applied. However, it was shown in the literature that those
assumptions are valid for the conditions when the CP product is much less than 1
(CP<<1).
CONDITIONS OF SOLUTION
a) Transient condition:
This condition is only applicable for a relatively short period after some pressure
disturbance has been created in the reservoir. In the time for which the transient
condition is applicable it is assumed that the pressure response in the reservoir is not
affected by the presence of outer boundary (Figure 3), thus the reservoir appears
INFINITE in extent. The condition is mainly applied to the analysis of well tests in
which the production rate of the wellbore is deliberately changed and the resulting
pressure response in the wellbore is measured and analyzed during a short period of
time (a few hours). Then, unless the reservoir is extremely small, the boundary
effects will not be felt and the reservoir is, mathematically infinite.
In this solution, both the pressure and pressure derivative, with respect to time, are
functions of both position and time.
5
𝜕𝑃
𝑃 = 𝑓(𝑟, 𝑡) 𝑎𝑛𝑑 = 𝑓(𝑟, 𝑡) (Eq. 15)
𝜕𝑡
𝜕𝑃
=0 𝑎𝑡 𝑟 = 𝑟𝑒 (Eq. 16)
𝜕𝑟
If the well is producing at a constant rate then the cell pressure will decline in a such
way that
𝜕𝑃
≈ 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 𝑓𝑜𝑟 𝑎𝑙𝑙 𝑟 𝑎𝑛𝑑 𝑡 (Eq. 17)
𝜕𝑡
6
𝑑𝑃 𝑑𝑉
𝐶𝑉 =− = −𝑞 (Eq. 18)
𝑑𝑡 𝑑𝑡
Figure 4. Pseudo-steady state flow diagram showing pressure distribution (Dake, 1978)
The steady state condition applies, after the transient period, to a well draining a cell
has a completely open outer boundary. It is assumed that, for a constant rate of
production, fluid withdrawal from the cell will be exactly balanced by fluid entry
across the outer boundary and therefore,
and
𝜕𝑃
=0 𝑓𝑜𝑟 𝑎𝑙𝑙 𝑟 𝑎𝑛𝑑 𝑡 (Eq. 21)
𝜕𝑡
7
This condition is appropriate when pressure is being maintained in the reservoir due
to either natural water influx or the injection of some displacing fluid.
Pe = constant
P
Pwf Fluid
Influx
rw r re
REFERENCES