...
Expand |
---|
|
The Volatile Oil flow model simulates 3-component fluid : water, liquid hydrocarbon (called "oil") and gaseous hydrocarbons ( called "gas") that flow in 3 possible phases (water, gasified oil and free gas) and defined by the following set of equations: Section |
---|
Column |
---|
| LaTeX Math Block |
---|
| \partial_t \bigg [ \phi \ \rho_W \bigg ] + \nabla \bigg ( \rho_{Ww} \ \mathbf{u}_w \bigg ) = q_{mW}(\mathbf{r}) |
LaTeX Math Block |
---|
| \partial_t \bigg [ \phi \ \rho_O \bigg ] + \nabla \bigg ( \rho_{Oo} \ \mathbf{u}_o
+ \rho_{Og} \ \mathbf{u}_g \bigg ) = q_{mO}(\mathbf{r}) |
LaTeX Math Block |
---|
| \partial_t \bigg [ \phi \ \rho_G \bigg ] + \nabla \bigg ( \rho_{Go} \ \mathbf{u}_o
+ \rho_{Gg} \ \mathbf{u}_g \bigg ) = q_{mG}(\mathbf{r}) |
|
Column |
---|
| LaTeX Math Block |
---|
anchor | DarcyW1 |
---|
alignment | left |
---|
| \mathbf{u}_w = - k_a \ \frac{k_{rw}(s_w, s_g)}{\mu_w} \ ( \nabla P_w - \rho_w \mathbf{g} ) |
LaTeX Math Block |
---|
anchor | DarcyO1 |
---|
alignment | left |
---|
| \mathbf{u}_o = - k_a \ \frac{k_{ro}(s_w, s_g)}{\mu_o} \ ( \nabla P_o - \rho_o \mathbf{g} ) |
LaTeX Math Block |
---|
anchor | DarcyG1 |
---|
alignment | left |
---|
| \mathbf{u}_g = - k_a \ \frac{k_{rg}(s_w, s_g)}{\mu_g} \ ( \nabla P_g - \rho_g \mathbf{g} ) |
|
Column |
---|
| LaTeX Math Block |
---|
anchor | CapilarOW1 |
---|
alignment | left |
---|
| P_o - P_w = P_{cow}(s_w) |
LaTeX Math Block |
---|
anchor | CapilarOG1 |
---|
alignment | left |
---|
| P_o - P_g = P_{cog}(s_g) |
LaTeX Math Block |
---|
anchor | swsosg1 |
---|
alignment | left |
---|
| s_w + s_o + s_g = 1 |
|
|
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} |
Equations LaTeX Math Block Reference |
---|
| – LaTeX Math Block Reference |
---|
| define 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 in space. Equations LaTeX Math Block Reference |
---|
| – LaTeX Math Block Reference |
---|
| define the motion dynamics of each phase, represnted as linear correlation between phase flow speed and partial pressure gradient of this phase (which is also called Darci flow with account of the gravity and relative permeability).
Equations LaTeX Math Block Reference |
---|
| – LaTeX Math Block Reference |
---|
| define the hydrodynamic inter-facial balance between the phases with account of capillary pressure in porous formation . The key assumption is that capillary pressure at oil-water boundary is a function of water saturation alone LaTeX Math Inline |
---|
body | P_{cow} = P_{cow}(s_w) |
---|
| and capillary pressure at oil-gas boundary is a function of gas saturation alone LaTeX Math Inline |
---|
body | P_{cog} = P_{cog}(s_g) |
---|
|
In the absence of capillary pressure the inter-facial equilibrium simplifies and implies that all phases are at the same pressure at all times.
Equations LaTeX Math Block Reference |
---|
| implies that porous space is fully occupied by fluid at all times .
Equation LaTeX Math Block Reference |
---|
| 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 LaTeX Math Inline |
---|
body | \frac{\delta E_H}{ \delta V \delta t} |
---|
| defines the speed of change of heat energy volumetric density.In impermeable rocks ( LaTeX Math Inline |
---|
body | \phi =0, \; \bar u_\alpha = 0 |
---|
| ) heat flow is defined by heat conduction only: LaTeX Math Block |
---|
| \rho_r \, c_{pr} \frac{\partial T}{\partial t} - \nabla (\lambda_t \nabla T) = \frac{\delta E_H}{ \delta V \delta t} |
The effective specific heat capacity of formation with multiphase flow is a simple sum of its components: LaTeX Math Block |
---|
| (\rho \,c_{pt})_p = (1-\phi) \rho_r \, \ c_{pr} + \phi \ (s_w \rho_w \, c_{pw} + s_o \rho_o \, c_{po} + s_g \rho_g \, c_{pg} ) |
The effective thermal conductivity of formation with multiphase flow is assumed to be a sum of its components: LaTeX Math Block |
---|
| \lambda_{t} = (1-\phi) \ \lambda_r + \phi \ (s_w \lambda_w + s_o \lambda_o + s_g \lambda_g ) |
The term LaTeX Math Inline |
---|
body | \bigg( \sum_{a = \{w,o,g \}} \rho_\alpha \ c_{p \alpha} \ \mathbf{u}_\alpha \bigg) \ \bar \nabla T |
---|
| represents heat convection defined by the mass flow. The term LaTeX Math Inline |
---|
body | \bigg( \sum_{a = \{w,o,g \}} \rho_\alpha \ c_{p \alpha} \ \epsilon_\alpha \ \mathbf{u}_\alpha \bigg) \bar \nabla P |
---|
| represents the heating/cooling effect of the multiphase flow through the porous media. This effect is the most significant with light oil and gases.
The term LaTeX Math Inline |
---|
body | \ \phi \sum_{a = \{w,o,g \}} \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 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.
The set LaTeX Math Block Reference |
---|
| – LaTeX Math Block Reference |
---|
| represent the system of 16 scalar equations on 16 unknowns: LaTeX Math Inline |
---|
body | \{ T, \ P_w, \ P_o, \ P_g, \ s_w, \ s_o, \ s_g, \ u_w^x, \ u_w^y, \ u_w^z, \ u_o^x, \ u_o^y, \ u_o^z, \ u_g^x, \ u_g^y, \ u_g^z \} |
---|
| ,which are all functions of time and space coordinates LaTeX Math Inline |
---|
body | (t, \mathbf{r}) = (t,x,y,z) |
---|
| .
Expressing the molar densities with mass shares and phase density (see also "Volatile Oil Model") one gets:
Section |
---|
Column |
---|
| LaTeX Math Block |
---|
| \partial_t \bigg [ \phi \ \rho_W \bigg ] + \nabla \bigg ( \rho_w \ \mathbf{u}_w \bigg ) = q_{mW}(\mathbf{r}) |
LaTeX Math Block |
---|
| \partial_t \bigg [ \phi \ \rho_O \bigg ] + \nabla \bigg ( {\tilde m}_{Oo} \ \rho_o \ \mathbf{u}_o
+ {\tilde m}_{Og} \ \rho_{g} \ \mathbf{u}_g \bigg ) = q_{mO}(\mathbf{r}) |
LaTeX Math Block |
---|
| \partial_t \bigg [ \phi \ \rho_G \bigg ] + \nabla \bigg ( {\tilde m}_{Go} \ \rho_{o} \ \mathbf{u}_o
+ {\tilde m}_{Gg} \ \rho_g \ \mathbf{u}_g \bigg ) = q_{mG}(\mathbf{r}) |
|
Column |
---|
| LaTeX Math Block |
---|
anchor | DarcyW1 |
---|
alignment | left |
---|
| \mathbf{u}_w = - k_a \ \frac{k_{rw}(s_w, s_g)}{\mu_w} \ ( \nabla P_w - \rho_w \mathbf{g} ) |
LaTeX Math Block |
---|
anchor | DarcyO1 |
---|
alignment | left |
---|
| \mathbf{u}_o = - k_a \ \frac{k_{ro}(s_w, s_g)}{\mu_o} \ ( \nabla P_o - \rho_o \mathbf{g} ) |
LaTeX Math Block |
---|
anchor | DarcyG1 |
---|
alignment | left |
---|
| \mathbf{u}_g = - k_a \ \frac{k_{rg}(s_w, s_g)}{\mu_g} \ ( \nabla P_g - \rho_g \mathbf{g} ) |
|
Column |
---|
| LaTeX Math Block |
---|
anchor | CapilarOW1 |
---|
alignment | left |
---|
| P_o - P_w = P_{cow}(s_w) |
LaTeX Math Block |
---|
anchor | CapilarOG1 |
---|
alignment | left |
---|
| P_o - P_g = P_{cog}(s_g) |
LaTeX Math Block |
---|
anchor | swsosg1 |
---|
alignment | left |
---|
| s_w + s_o + s_g = 1 |
|
|
Substituting the values of mass densities and mass shares of fluid components (см. "Volatile Oil Model") and dividing each equation by density of corresponding component in standard conditions one gets the most popular form of Volatile Oil flow equations:
|
Initial Conditions
...
Начальное условие по температуре задается распределением температурного поляInitial temperature distribution is set as input:
LaTeX Math Block |
---|
|
T(0, \mathbf{r}) = T_0(\mathbf{r}) |
Начальное условие на давления, скоростей и насыщенности задается одним из двух вариантов.
Условие I – Стационарный старт
In case the simulation is performed over the undisturbed reservoir then initial temperature distribution is geothermal.
The initial condition on phase pressure, phase velocities and phase saturations is set by one of the following options: Equilibrium Start and Non-equilibrium Start.
Condition I – Equilibrium Start
Equilibrium Start means that phase pressure Стационарный старт означает, что до начального момента времени поле давлений
LaTeX Math Inline |
---|
body | \{ P_w, \ P_o, \ P_g \} |
---|
|
, скоростей and phase saturations LaTeX Math Inline |
---|
body | \{ \mathbf{u}s_w, \ \mathbf{u}s_o, \ \mathbf{u}s_g, \} |
---|
|
и насыщенностей were in stationary (not varying in time) conditions and phase velocities LaTeX Math Inline |
---|
body | \{ \mathbf{ su}_w, \ \ smathbf{u}_o, \ \ smathbf{u}_g , \} |
---|
|
находилось в стационарном (не меняющемся во времени) состоянии, соответствующем гидродинамическому равновесию were zero, corresponding to hydrodynamic equilibrium:
Section |
---|
Column |
---|
| LaTeX Math Block |
---|
| \nabla \cdot \bigg ( \frac{1}{B_w} \ \mathbf{u}_w \bigg )_{t=0} = 0 |
LaTeX Math Block |
---|
| \nabla \cdot \bigg ( \frac{1}{B_o} \ \mathbf{u}_o
+ \frac{R_v}{B_g} \ \mathbf{u}_g \bigg )_{t=0} = 0 |
LaTeX Math Block |
---|
| \nabla \cdot \bigg ( \frac{1}{B_g} \ \mathbf{u}_g
+ \frac{R_s}{B_o} \ \mathbf{u}_o \bigg )_{t=0} = 0 |
|
Column |
---|
| LaTeX Math Block |
---|
anchor | DarcyW |
---|
alignment | left |
---|
| \mathbf{u}_w(0, \mathbf{r}) = - k_a \ \frac{k_{rw}(s_w, s_g)}{\mu_w} \ (\nabla P_w(0, \mathbf{r}) - \rho_w \ \mathbf{g} ) = 0 |
LaTeX Math Block |
---|
anchor | DarcyO |
---|
alignment | left |
---|
| \mathbf{u}_o(0, \mathbf{r}) = - k_a \ \frac{k_{ro}(s_w, s_g)}{\mu_o} \ ( \nabla P_o(0, \mathbf{r}) - \rho_o \ \mathbf{g} ) = 0 |
LaTeX Math Block |
---|
anchor | DarcyG |
---|
alignment | left |
---|
| \mathbf{u}_g(0, \mathbf{r}) = - k_a \ \frac{k_{rg}(s_w, s_g)}{\mu_g} \ ( \nabla P_g(0, \mathbf{r}) - \rho_g \ \mathbf{g} ) = 0 |
|
Column |
---|
| LaTeX Math Block |
---|
anchor | CapilarOW |
---|
alignment | left |
---|
| P_o(0, \mathbf{r}) - P_w(0, \mathbf{r}) = P_{cow}(s_w) |
LaTeX Math Block |
---|
anchor | CapilarOG |
---|
alignment | left |
---|
| P_o(0, \mathbf{r}) - P_g(0, \mathbf{r}) = P_{cog}(s_g) |
LaTeX Math Block |
---|
anchor | swsosg |
---|
alignment | left |
---|
| s_w + s_o + s_g = 1 |
|
|
...
Condition II – Non-equilibrium Start
Equilibrium Start means that phase pressure
LaTeX Math Inline |
---|
body | \{ P_w, \ P_o, \ P_g \} |
---|
|
, phase velocities LaTeX Math Inline |
---|
body | \{ \mathbf{u}_w, \ \mathbf{u}_o, \ \mathbf{u}_g \} |
---|
|
and phase saturations LaTeX Math Inline |
---|
body | \{ s_w, \ s_o, \ s_g, \} |
---|
|
were in stationary (not varying in time) conditions, corresponding to hydrodynamic equilibrium:Нестационарный старт означает, что к начальному моменту времени поле насыщенностей
LaTeX Math Inline |
---|
body | \{ s_w, \ s_o, \ s_g, \} |
---|
|
является произвольным, с условием
...