Derivation of The Fluid Flow Equation
Derivation of The Fluid Flow Equation
Conservation of Mass
The law of conservation of mass is a material-balance equation written for a component in
a control volume of the system to be modeled. In petroleum reservoirs, the control volume is made
up of a porous medium containing one, two or three phases. This section concerns single-phase
fluid. The porous medium is treated as a continuum whose physical properties at any point are
those of a representative element of the medium. The control volume, whose shape depends on the
coordinate system used in the model.
Now, consider an element of a reservoir through which a single-phase is flowing in the x-
direction.
4
University of Baghdad Reservoir Simulation
Department of Petroleum Engineering 2020
Now the simultaneous limits over time and space that is, look at an instant of time as the
control volume shrinking to infinitesimal.
𝑣 ∆ 𝜌 ∆ 𝑣 𝜌 𝜌 ∆ 𝜌
lim lim ∅
∆ → ∆𝑥 ∆ → ∆𝑡
In the above equation, we can recognize the definition of the first-order partial derivative with
respect to time and space coordinates (x and t).
For the space variable:
𝜕𝑓 𝑓 𝑥 ∆𝑥 𝑓 𝑥
lim
𝜕𝑥 ∆ → ∆𝑥
5
University of Baghdad Reservoir Simulation
Department of Petroleum Engineering 2021
Rate Equation
Darcy’s law is an empirical relationship between fluid flow rate through a porous media
and potential gradient. For single-phase, one-dimensional flow, this law can be expressed in a
differential form as:
𝑞 𝑘 𝜕𝑝
𝑣
𝐴 𝜇 𝜕𝑥
Substitute the above equation into the continuity equation:
𝑘 𝜕𝑝
𝜕 𝜌 𝜕𝜌
𝜇 𝜕𝑥
∅
𝜕𝑥 𝜕𝑡
The flow equation can be written as follows by expanding the left-hand side:
𝜌 ∅ (B)
𝑘 𝜕 𝑝 𝑘 𝜕𝜌 𝜕𝑝 𝜕𝜌 𝜕𝑝
𝜌 ∅
𝜇 𝜕𝑥 𝜇 𝜕𝑝 𝜕𝑥 𝜕𝑝 𝜕𝑡
Neglecting the term since we are going to assume small pressure gradient, the above
equation becomes
𝜌 ∅ (C)
Equation of State
The equation of state is needed to express the density in terms of pressure. Most of oil
field systems are considered to be slightly compressible.
In this case, the equation of state is:
𝜌 𝜌 𝑒
6
University of Baghdad Reservoir Simulation
Department of Petroleum Engineering 2021
𝜕 𝑝 ∅ 𝜇 𝑐 𝜕𝜌
𝜕𝑥 𝑘 𝜕𝑡
𝑘 𝜕𝑝
𝜕 𝜕𝜌
𝜇 𝜕𝑥
∅𝑐
𝜕𝑥 𝜕𝑡
This equation is generally called the diffusivity equation.
The following equations are for other coordinate systems:
For radial flow,
𝜕 𝑝 1 𝜕𝑝 ∅ 𝜇 𝑐 𝜕𝜌
𝜕𝑟 𝑟 𝜕𝑟 𝑘 𝜕𝑡
For two- dimensional flow,
𝜕 𝑝 𝜕 𝑝 ∅ 𝜇 𝑐 𝜕𝑝
𝜕𝑥 𝜕𝑦 𝑘 𝜕𝑡
For three- dimensional flow,
𝜕 𝑝 𝜕 𝑝 𝜕 𝑝 ∅ 𝜇 𝑐 𝜕𝑝
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑘 𝜕𝑡