Mathematical model of multiphase wellbore flow predicts the temperature, pressure and flow speed distribution along the wellbore trajectory with account for:
(\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 disambiguation fo the properties in the above equation is brought in The list of dynamic flow properties and model parameters.
Equations ??? define the continuity of the fluid components flow or equivalently represent the mass conservation of each mass component during its transportation in space.
Equations ??? define the motion dynamics of each phase, represented as linear correlation between phase flow speed and partial pressure gradient of this phase .
Equation defines the heat flow continuity or equivalently represents heat conservation due to heat conduction and convection with account for adiabatic and Joule–Thomson throttling effect.
The term defines the speed of change of heat energy volumetric density due to the inflow from formation into the wellbore.
The term represents heat convection defined by the wellbore mass flow.
The term represents the heating/cooling effect of the multiphase flow through the porous media. This effect is the most significant with light oils and gases.
The term represents the heating/cooling effect of the fast adiabatic pressure change. This usually takes effect in and around the wellbore during the first minutes or hours after changing the well flow regime (as a consequence of choke/pump operation). This effect is absent in stationary flow and negligible during the quasi-stationary flow and usually not modeled in conventional monthly-based flow simulations.
Stationary wellbore flow is defined as the flow with constant pressure and temperature: and .
This happens during the long-term (usually hours & days & weeks) production/injection or long-term (usually hours & days & weeks) shut-in.
The temperature dynamic equation is going to be:
\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} |
and its discrete computational scheme will be:
\bigg( \sum_{a = \{w,o,g \}} \rho_\alpha^{k-1} \ c_{p \alpha}^{k-1} \ q_\alpha^{k-1} \bigg) T^{k-1} - \bigg( \sum_{a = \{w,o,g \}} \rho_\alpha^k \ c_{p \alpha}^k \ q_\alpha^k \bigg) T^k = \sum_{a = \{w,o,g \}} \rho_\alpha^k \ c_{p \alpha}^k \ (q_\alpha^{k-1} - q_\alpha^k) \, (T_r^k + \epsilon_\alpha^k \delta p^k ) |
where is drawdown, – formation pressure in th grid layer, – bottom-hole pressure across th grid layer, – remote reservoir temperature of th grid layer.
The axis is pointing downward along hole with th grid layer sitting above the th grid layer.
If the flowrate is not vanishing during the stationary lift () then can be calculated iteratively from previous values of the wellbore temperature as:
T^{k-1} = \frac{\bigg( \sum_{a = \{w,o,g \}} \rho_\alpha^k \ c_{p \alpha}^k \ q_\alpha^k \bigg) T^k + \sum_{a = \{w,o,g \}} \rho_\alpha^k \ c_{p \alpha}^k \ (q_\alpha^{k-1} - q_\alpha^k) \, (T_r^k + \epsilon_\alpha^k \delta p^k )}{\bigg( \sum_{a = \{w,o,g \}} \rho_\alpha^{k-1} \ c_{p \alpha}^{k-1} \ q_\alpha^{k-1} \bigg) } |
The wellbore fluid velocity can be expressed thorugh the volumetric flow profile and tubing/casing cross-section area as:
so that
|
Beggs, H. D. and Brill, J. P.: "A Study of Two-Phase Flow in Inclined Pipes," J. Pet. Tech., May (1973), 607-617
|
time and space corrdinates , -axis is orientated towards the Earth centre, define transversal plane to the -axis | |
position vector at which the flow equations are set | |
speed of water-component mass change in wellbore draining points | |
speed of oil-component mass change in wellbore draining points | |
speed of gas-component mass change in wellbore draining points | |
volumetric water-component flow rate in wellbore draining points recalculated to standard surface conditions | |
volumetric oil-component flow rate in wellbore draining points recalculated to standard surface conditions | |
volumetric gas-component flow rate in wellbore draining points recalculated to standard surface conditions | |
volumetric water-phase flow rate in wellbore draining points | |
volumetric oil-phase flow rate in wellbore draining points | |
volumetric gas-phase flow rate in wellbore draining points | |
total well volumetric water-component flow rate | |
total well volumetric oil-component flow rate | |
total well volumetric gas-component flow rate | |
total well volumetric liquid-component flow rate | |
water-phase pressure pressure distribution and dynamics | |
oil-phase pressure pressure distribution and dynamics | |
gas-phase pressure pressure distribution and dynamics | |
water-phase flow speed distribution and dynamics | |
oil-phase flow speed distribution and dynamics | |
gas-phase flow speed distribution and dynamics | |
capillary pressure at the oil-water phase contact as function of water saturation | |
capillary pressure at the oil-gas phase contact as function of gas saturation | |
relative formation permeability to water flow as function of water and gas saturation | |
relative formation permeability to oil flow as function of water and gas saturation | |
relative formation permeability to gas flow as function of water and gas saturation | |
porosity as function of formation pressure | |
absolute formation permeability to air | |
gravitational acceleration vector | |
gravitational acceleration constant | |
mass density of -phase fluid | |
viscosity of -phase fluid | |
effective thermal conductivity of the rocks with account for multiphase fluid saturation | |
rock matrix thermal conductivity | |
thermal conductivity of -phase fluid | |
rock matrix mass density | |
differential adiabatic coefficient of -phase fluid | |
specific isobaric heat capacity of the rock matrix | |
specific isobaric heat capacity of -phase fluid | |
differential Joule–Thomson coefficient of -phase fluid дифференциальный коэффициент Джоуля-Томсона фазы |