The multiphase wellbore flow in hydrodynamic and thermodynamic equilibrium is defined by the following set of equations:
LaTeX Math Block |
---|
anchor | MatBal |
---|
alignment | left |
---|
|
\frac{\partial (\rho_m A)}{\partial t} + {\nabla } \bigg( A \, \sum_\alpha \rho_\alpha \, { \bf u}_\alpha \bigg) = 0 |
LaTeX Math Block |
---|
|
\rho_\alpha \bigg[ \frac{\partial {\bf u}_\alpha}{\partial t} + ({\bf u}_\alpha \cdot \nabla) \ {\bf u}_\alpha - \nu_\alpha \Delta {\bf u}_\alpha\bigg] = - \nabla p + \rho_\alpha \, {\bf g} - \frac{ f_\alpha \, \rho_\alpha \, u_\alpha}{2 d} \ {\bf u}_\alpha |
LaTeX Math Block |
---|
|
(\rho \,c_{pt})_p \frac{\partial T}{\partial t}
- \sum_\alpha \rho_\alpha \ c_{p \alpha} \ \eta_{s \alpha} \ \frac{\partial P_\alpha}{\partial t}
+ \sum_\alpha \rho_\alpha \ c_{p \alpha} \ ({\bf u}_\alpha \cdot \nabla) \ T
\ = \ \frac{\delta E_H}{ \delta V \delta t} |
Equation
LaTeX Math Block Reference |
---|
|
defines the mass flow continuity or equivalently represents mass conservation during the mass transportation.
Equations
LaTeX Math Block Reference |
---|
|
are Navie-Stokes equations and define the momentum conservation during the mass transportation.
Equation
LaTeX Math Block Reference |
---|
|
defines the heat flow continuity or equivalently represents heat conservation due to heat
conduction exchange with surrounding rocks and
mass convection
from reservoir inflow with account for
wellbore adiabatic
and Joule–Thomson effects and Joule–Thomson reservoir throttling effect.
The term
LaTeX Math Inline |
---|
body | \frac{\delta E_H}{ \delta V \delta t} |
---|
|
defines the speed of change of heat energy
volumetric density due to
heat exchange with surrounding rocks and the inflow from formation into the wellbore.
Projecting the above equations to the well trajectory
with down-pointing along-hole coordinate
:
LaTeX Math Block |
---|
anchor | MatBal |
---|
alignment | left |
---|
|
\frac{\partial (\rho_m A)}{\partial t} + \frac{\partial}{\partial l} \bigg( A \, \sum_\alpha \rho_\alpha \, u_\alpha \bigg) = 0 |
LaTeX Math Block |
---|
|
\sum_\alpha \rho_\alpha \bigg[ \frac{\partial u_\alpha}{\partial t} + u_\alpha \frac{\partial u_\alpha}{\partial l} - \nu_\alpha \Delta u_\alpha\bigg] = - \frac{dp}{dl} + \rho_m \, g \, \sin \theta - \frac{ f_m \, \rho_m \, u_m^2 \, }{2 d} |
LaTeX Math Block |
---|
|
(\rho \,c_p)_m \frac{\partial T}{\partial t}
- \bigg( \sum_\alpha \rho_\alpha \ c_{p \alpha} \ \eta_{s \alpha}\bigg) \ \frac{\partial p}{\partial t}
+ \bigg( \sum_\alpha \rho_\alpha \ c_{p \alpha} \ u_\alpha \bigg) \frac{\partial T}{\partial l}
\ = \ \frac{1}{A} \ \sum_\alpha \rho_\alpha \ c_{p \alpha} T_\alpha \frac{\partial q_\alpha}{\partial l} |
Equations
LaTeX Math Block Reference |
---|
|
–
LaTeX Math Block Reference |
---|
|
define a closed set of 3 scalar equations on 3 unknowns: pressure
, temperature
and mixture-average fluid velocity
.
The disambiguation of the properties in the above equation is brought in The list of dynamic flow properties and model parameters.
Equation
LaTeX Math Block Reference |
---|
|
defines the continuity of the fluid components flow or equivalently represent the mass conservation of each mass component
LaTeX Math Inline |
---|
body | \{ m_W, \ m_O, \ m_G \} |
---|
|
during its transportation along wellbore.
Equation
LaTeX Math Block Reference |
---|
|
defines the motion dynamics of each phase (called Navier–Stokes equation), represented as linear correlation between phase flow speed
and pressure profile of mutliphase fluid
.
The term
LaTeX Math Inline |
---|
body | \sum_\alpha \rho_\alpha \ c_{p \alpha} \ u_\alpha \frac{\partial T}{\partial l} |
---|
|
represents heat convection defined by the wellbore mass flow.
The term
LaTeX Math Inline |
---|
body | \sum_\alpha \rho_\alpha \ c_{p \alpha} \ \eta_{s \alpha} \ \frac{\partial P_\alpha}{\partial t} |
---|
|
represents the heating/cooling effect of the fast adiabatic pressure change.
This usually takes effect in the wellbore during the first minutes or hours after changing the well flow regime (as a consequence of choke/pump operation).
The term
defines mass-specific heat capacity of the multiphase mixture and defined by exact formula:
LaTeX Math Block |
---|
anchor | rho_cp |
---|
alignment | left |
---|
|
(\rho \,c_p)_m = \sum_\alpha \rho_\alpha c_\alpha s_\alpha |
Stationary wellbore flow is defined as the flow with constant pressure and temperature:
LaTeX Math Inline |
---|
body | \frac{\partial T}{\partial t} = 0 |
---|
|
and
LaTeX Math Inline |
---|
body | \frac{\partial P_\alpha}{\partial t} = 0 |
---|
|
.
This happens during the long-term (usually hours & days & weeks) production/injection or long-term (usually hours & days & weeks) shut-in.
LaTeX Math Block |
---|
|
(\rho \,c_{pt})_p \frac{\partial T}{\partial t}
- \ \phi \sum_{a = \{w,o,g \}} \rho_\alpha \ c_{p \alpha} \ \eta_{s \alpha} \ \frac{\partial P_\alpha}{\partial t}
+ \bigg( \sum_{a = \{w,o,g \}} \rho_\alpha \ c_{p \alpha} \ \epsilon_\alpha \ \mathbf{u}_\alpha \bigg) \nabla P
+ \bigg( \sum_{a = \{w,o,g \}} \rho_\alpha \ c_{p \alpha} \ \mathbf{u}_\alpha \bigg) \ \nabla T
- \nabla (\lambda_t \nabla T) = \frac{\delta E_H}{ \delta V \delta t} |
The wellbore fluid velocity
can be expressed thorugh the volumetric flow profile
and tubing/casing cross-section area
as:
LaTeX Math Block |
---|
|
u_\alpha = \frac{q_\alpha}{\pi r_f^2} |
so that
LaTeX Math Block |
---|
|
\bigg( \sum_{a = \{w,o,g \}} \rho_\alpha \ c_{p \alpha} \ \mathbf{u}_\alpha \bigg) \ \nabla T
= \frac{\delta E_H}{ \delta V \delta t} |