Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.


Specific case of general multi-phase pressure diffusion with 3-phase Oil + Gas + Water fluid model which assumes that pressure diffusion is equivalent to single-phase diffusion with specifically averaged dynamic parameters, thus resulting in linear partial differential equation with constant coefficients: 

LaTeX Math Block
anchorDiffP
alignmentleft
\phi \, c_t \, \partial_t p - \nabla \big( M \cdot ( \nabla p - \rho \cdot \mathbf{g} ) \big)  = \sum_k \, q_k(t) \cdot \delta(\mathbf{r} -\mathbf{r}_k)

where

LaTeX Math Inline
bodyt

time

LaTeX Math Inline
body--uriencoded--%7B\rm r%7D = (x,y,z)

reservoir location

LaTeX Math Inline
body--uriencoded--\mathbf%7Br%7D_k

Well–reservoir contact for 

LaTeX Math Inline
bodyk
-th  well

Если в исследуемой системе давление и распределение насыщенности в пласте (как вблизи, так и вдали от скважин) не успели сильно измениться за время исследования, то можно принять, что коэффициенты уравнения 

LaTeX Math Block Reference
anchorPZ
pageMulti-phase diffusion model
 не зависят от давления (рассчитываются при одном значении опорного давления в пласте) и нелинейная мультифазная фильтрация может быть представлена как линейная фильтрация некого эффективного однофазного флюида, называемая моделью Перрина (Perrine):

LaTeX Math Block
anchorD3Z52
alignmentleft
\phi c_t \partial_t P - \nabla \big( \alpha ( \nabla P - \rho_{\alpha} \mathbf{g} ) \big)  = q(\mathbf{r})

где

LaTeX Math Block
anchor1
alignmentleft
q_t =
 B_w \, q_W
+ \frac{B_o - R_s B_g}{1 - R_v R_s} \, q_O
+ \frac{B_g - R_v B_o}{1 - R_v R_s} \, q_G =
 B_w \, q_W
+ B_o \, q_O \, \frac{1 - R_{sn}}{1 - R_{vn} R_{sn}} 
+ B_g  \, q_G \, \frac{1 - R_{vn}}{1 - R_{vn} R_{sn}}

суммарный отбор воды, нефти и газа в пластовых условиях


LaTeX Math Block
anchor1
alignmentleft
B_w
p = 
B_w(P_{\rm ref}), \ B_o = B_o(P_{\rm ref}), \ B_g = B_g(P_{\rm ref})
\frac{1}{3} \cdot \left( p_w + p_o + p_g \right)



3-phase average reservoir pressure


LaTeX Math Block
anchorqt

объемный фактор воды, нефти и газа при опорном давлении

LaTeX Math Inline
bodyP_{\bf ref}

LaTeX Math Blockanchor1
alignmentleft
R
q_t(\mathbf{
vn
r}) =
\frac{R_v B_o}{B_g} \bigg| _{P_{\bf ref}} \ , \quad R_{sn} = \frac{R_s B_g}{B_o} \bigg| _{P_{\bf ref}} \ , \quad R_{vp} = \frac{\dot R_v B_o}{B_g} \bigg| _{P_{\bf ref}} \ , \quad R_{sp} = \frac{\dot R_s B_g}{B_o} \bigg| _{P_{\bf ref}}нормированные обменные коэффициенты межфазного обмена при опорном давлении
 q_w + q_o + q_g = B_w \, q_W + (B_o - R_s \, B_g) \, q_O + (B_g - R_v \, B_o) \, q_G



total sandface flowrate
at reservoir location

LaTeX Math Inline
body\mathbf{r}


LaTeX Math Block
anchor1
alignmentleft
B_w, \ B_o, \ B_g 



formation volume factors
at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchor1
alignmentleft
\phi(\mathbf{r})



effective porosity
in reservoir location

LaTeX Math Inline
body\bf r
at reference pressure
LaTeX Math Inline
bodyp_{\rm

LaTeX Math InlinebodyP_{\bf

ref}


LaTeX Math Block
anchor0OM3S
alignmentleft
s(\mathbf{r}) = \{ s_w(\mathbf{r}), \ s_o(\mathbf{r}), \ s_g(\mathbf{r})  
\}распределение насыщенности водяной, нефтяной и газовой фаз в объеме пород
\}



reservoir saturation
as a function of location

LaTeX Math Inline
body\bf r



LaTeX Math Block
anchor0OM3S
alignmentleft
c_t
(s)
 = c_r
(1
 + c_w s_w +  
R
c_
{sn}
o s_o +  c_g s_g  + s_o [ R_{
vn
sp}
s_g)
 + (c_
w
r 
s_w
 + c_o) 
s_o
 
(1+
R_{sn}
)
 ] + 
c
s_g [ 
s
R_
g (1
{vp} + R_{vn
}) + R_{sp} s_o + R_{vp} s_gэффективная сжимаемость пласта с мультифазным насыщением как функция насыщенности при опорном давлении
}(c_r + c_g) ]



LaTeX Math Block
anchor1
alignmentleft
с_r



reservoir pore compressibility
at reference pressure

LaTeX Math Inline
body

P

p_{\

bf

rm ref}


LaTeX Math Block
anchor1
alignmentleft
с_
r
w, 
= - \frac{1}{\phi} \frac{d \phi}{dP} \bigg| _{P_{\bf ref}}, \quad с_w = \frac{1}{\rho_w} \frac{d \rho_w}{dP} \bigg| _{P_{\bf ref}}, \quad с_o = \frac{1}{\rho_o} \frac{d \rho_o}{dP} \bigg| _{P_{\bf ref}}, \quad с_g = \frac{1}{\rho_g} \frac{d \rho_g}{dP} \bigg| _{P_{\bf ref}} сжимаемости породы, воды, нефти и газа при опорном давлении
\ с_o, \ с_g



fluid compressibilities
at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchorRHTVX
alignmentleft
M = M_w + M_o \big( 1 + R_{sn} \big) + M_g \big( 1 + R_{vn} \big)



total fluid mobility
at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchor6OXKP
alignmentleft
M_w = k_a \cdot M_{rw}



water mobility
at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm

LaTeX Math InlinebodyP_{\bf

ref}


LaTeX Math Block
anchor
1
6OXKP
alignmentleft
\beta
M_o = 
\phi
k_a \
,
cdot 
c_t \bigg| _{P_{\bf ref}}упругоемкость мультифазного пласта
M_{ro}



oil mobility at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchor
RHTVX
6OXKP
alignmentleft
\alpha(s) = \Big \langle \frac{k} {\mu} \Big \rangle
M_g = 
\alpha_w(s) + \alpha_o(s) \big( 1 + R_{sn} \big) + \alpha_g(s) \big( 1 + R_{vn} \big)эффективная проводимость пород  как функция насыщенности при опорном давлении
k_a \cdot M_{rg}



gas mobility at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchorQBU02
alignmentleft
M_{rw} = \frac{k_{rw}(s)}{\mu_w}



relative water mobility
at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm

LaTeX Math InlinebodyP_{\bf

ref}


LaTeX Math Block
anchor
6OXKP
QBU02
alignmentleft
\alpha_w = k_a \alpha_{rw}(s), \quad \alpha_o
M_{ro} = \frac{k_
a \alpha_
{ro}(s)
, \quad \alpha_g = k_a \alpha_{rg}(s)фазовые проводимости пласта по каждой фазе как функции насыщенности при опорном давлении
}{\mu_o}



relative oil mobility at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
body

P

T_{\

bf

rm ref}


LaTeX Math Block
anchor
1
QBU02
alignmentleft
M_{rg} = \frac{k_
a(\mathbf
{
rраспределение воздушной проницаемости пласта в объеме пород
rg}(s)
\bigg| _{P_{\bf ref}}
}{\mu_g}



relative gas mobility at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchor
QBU02
1
alignmentleft
\alpha_{rw}(s) = \frac{k_{rw}(s)}{\mu_w}, \quad \alpha_{ro}(s) = \frac{k_{ro}(s)}{\mu_o}, \quad \alpha_{rg}(s) = \frac{k_{rg}(s)}{\mu_g}
относительные фазовые проводимости пласта по каждой фазе как функции насыщенности при опорном давлении
k_a(\mathbf{r})



absolute permeability
as a function of location

LaTeX Math Inline
body\bf r
at reference pressure
LaTeX Math Inline
bodyp_{\rm ref}



LaTeX Math Block
anchor1
alignmentleft
\mu_w, \ \mu_o, \ \mu_g



water
, oil, gas dynamic viscosity at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm

LaTeX Math InlinebodyP_{\bf

ref}



LaTeX Math Block
anchor1
alignmentleft
\mu_w = \mu_w(P_{\bf ref}), \quad \mu_o
R_{sn} = \
mu_o(P_{\bf ref})
frac{R_s B_g}{B_o} \ , \quad 
\mu_g
R_{vn} = \
mu_g(P
frac{R_v B_o}{B_g}



normalized cross-phase exchange ratios
at reference pressure

LaTeX Math Inline
bodyp_{\

bf

rm ref}

)вязкость воды, нефти и газа при опорном давлении
P

T_{\

bf

rm ref}


LaTeX Math Block
anchor
QL1DV
1
alignmentleft
\rho_{\alpha
R_{sp} = \frac{\dot R_s B_g}{B_o} \ , \quad R_{vp} = \frac{
\
alpha_{rw} \rho_w + \alpha_{ro} \rho_o (1 + R_{sn}) + \alpha_{rg} \rho_g (1+R_{vn}) }{ \alpha_{rw} + \alpha
dot R_v B_o}{B_g}



LaTeX Math Block
anchorQL1DV
alignmentleft
\rho = \frac{ M_{rw} \rho_w + M_{ro}  (1 + R_{sn}) \rho_o  + 
\alpha
M_{rg}  (1+R_{vn})
}

гравитационная компонента потока при опорном давлении

LaTeX Math Inline
bodyP_{\bf ref}

LaTeX Math Block
anchorQL1DV
alignmentleft
\rho_w =
 \rho_
w(P
g }{ M_{
\rm ref}), \ \rho_o = \rho_o(P_{\rm ref}), \ \rho_g = \rho_g(P_{\rm ref}) плотность воды, нефти а газа при опорном давлении
rw}  + M_{ro}  (1 + R_{sn})  + M_{rg}  (1+R_{vn}) }



mobility-weighted fluid density
at reference pressure

LaTeX Math Inline
body

P

p_{\

bf

rm ref}

mathblock
anchor1
alignmentleft
 g = 9.81 \ \textrm{m} / \textrm{s}^2
ускорение свободного падения (константа)

LaTeX Math Inline
bodyT_{\rm ref}


LaTeX Math Block
anchor
chialignmentleft
 \chi = \frac{\alpha}{\beta} = \Big \langle \frac{k} {\mu} \Big \rangle \frac{1}{\phi c_t} \bigg| _{P_{\bf ref}}
пьезопроводность пласта при опорном давлении
1
alignmentleft
g = 9.81 \ \textrm{m} / \textrm{s}^2



standard gravity
LaTeX Math InlinebodyP_{\bf ref}


LaTeX Math Block
anchor
1
der
alignmentleft
\sigma
 
=
\
alpha
big 
\
( 
h
 
=
 \
Big
big)^{\LARGE \
langle
cdot} = \frac{
k} {\mu} \Big \rangle \ h \bigg| _{P_{\bf ref}}гидропроводность пласта при опорном давлении
d}{dp}



differentiation
with respect to the pressure


All the above dynamic properties are calculated at reference pressure

LaTeX Math Inline
bodyp_{\rm ref}
and temperature
LaTeX Math Inline
bodyT_{\rm ref}
 thus making 
LaTeX Math Block Reference
anchorDiffP
 a linear partial differential equation.

The reference temperature 

LaTeX Math Inline
body

...

T_{\

...

rm ref}

...

LaTeX Math Block
anchor18
alignmentleft
h

  is more or less in practical cases.

The choice of the reference pressure 

LaTeX Math Inline
bodyp_{\rm ref}
 depends on the task.


For accurate modelling of early time pressure response (ETR) it is recommended to use initial bottom hole pressure at the moment of the test: 

LaTeX Math Block
anchor1
alignmentleft
p_{\rm ref} = p_{wf}(t=0)


For accurate modelling of late time pressure response (LTR) it is recommended to use current formation pressure: 

LaTeX Math Block
anchor2
alignmentleft
p_{\rm ref} = p_e


The above equations 

LaTeX Math Block Reference
anchorDiffP
 – 
LaTeX Math Block Reference
anchorder
 is generalization of original model ( [1][2])  to the case of Volatile Oil fluid model



In case of Black Oil fluid model (

LaTeX Math Inline
bodyR_v = R_{vn} = R_{vp} = 0
) some equations are simplified:


LaTeX Math Block
anchor1
alignmentleft
q_t

толщина пласта

Все вышеприведенные параметры рассчитываются при одном давлении (которое предполагается не сильно меняющимся во время исследования), как правило при пластовом давлении.

При небольших депрессиях (репрессиях) это условие можно считать разумным.

Вышеприведенные формулы 

LaTeX Math Block Reference
anchor1
 – 
LaTeX Math Block Reference
anchor18
 представляют собой обобщение оригинальной модели ( [1], [2] )  на случай летучей нефти.

Для случае черной нефти (

LaTeX Math Inline
bodyR_v = R_{vn} = R_{vp} = 0
) некторые из вышеприведенных формул упрощаются:

LaTeX Math Block
anchor1
alignmentleft
q_t = B_w \, q_W + (B_o - R_s B_g) \, q_O + B_g \, q_G
 =  B_w \, q_W + B_o \, q_O + B_g \, ( q_G - R_s \, q_O)
суммарный отбор воды, нефти и газа в пластовых условиях



LaTeX Math Block
anchor0OM3S
alignmentleft
c_t(s) = c_r (1 + R_{sn} s_o ) + c_w s_w + c_o s_o (1+R_{sn}) + c_g s_
g + R_{sp} s_o эффективная сжимаемость пласта с мультифазным насыщением как функция насыщенности при опорном давлении LaTeX Math InlinebodyP_{\bf ref}
g  + R_{sp} s_o 



LaTeX Math Block
anchorRHTVX
alignmentleft
\alpha
M(s) = \Big \langle \frac{k} {\mu} \Big \rangle = 
\alpha
M_w(s) + 
\alpha
M_o(s) \big( 1 + R_{sn} \big) + \alpha_g(s) 
эффективная проводимость пород  как функция насыщенности при опорном давлении LaTeX Math InlinebodyP_{\bf ref}



LaTeX Math Block
anchorQL1DV
alignmentleft
\rho
_{\alpha}
 = \frac{ 
\alpha
M_{rw} \rho_w + 
\alpha
M_{ro}
\rho_o
(1 + R_{sn}) \rho_o  + 
\alpha
M_{rg} \rho_g  }{ 
\alpha
M_{rw}  + 
\alpha
M_{ro}  (1 + R_{sn})  + 
\alpha
M_{rg}  }

гравитационная компонента потока при опорном давлении

LaTeX Math Inline
bodyP_{\bf ref}



For the 2-phase Oil + Water fluid model (where  Linear Perrine multi-phase diffusion mode is the most accurate ) the above equations are getting even simpler:Для режима двухфазной фильтрации недонасыщенной нефти (а именно в этом случае модель Перрина дает наиболее удовлетворительное количественное приближение) вышеприведенные формулы упрощаются еще больше: 


LaTeX Math Block
anchor1
alignmentleft
q_t = B_w \, q_W + B_o  q_
Oсуммарный отбор воды, нефти и газа в пластовых условиях
O



LaTeX Math Block
anchor0OM3S
alignmentleft
c_t(s) = c_r  + c_w s_w + c_o s_o

эффективная сжимаемость пласта с мультифазным насыщением как функция насыщенности при опорном давлении

LaTeX Math Inline
bodyP_{\bf ref}



LaTeX Math Block
anchor
RHTVX
1
alignmentleft
\alpha
M(s) = \Big \langle \frac{k} {\mu} \Big \rangle = 
\alpha
M_w
(s) + \alpha_o(s)эффективная проводимость пород  как функция насыщенности при опорном давлении LaTeX Math InlinebodyP_{\bf ref}
(s) + M_o(s)



LaTeX Math Block
anchorQL1DV
alignmentleft
\rho
_{\alpha}
 = \frac{ 
\alpha
M_{rw} \rho_w + 
\alpha
M_{ro} \rho_o   }{ 
\alpha
M_{rw}  + 
\alpha
M_{ro}}

гравитационная компонента потока при опорном давлении

LaTeX Math Inline
bodyP_{\bf ref}

Несмотря на то, что приближение Перрина на практике чаще всего приводит к недостаточно аккуратным количественным оценкам давления оно, тем не менее:
 

  • дает правильное представление о характере отклика давления на мультифазные отборы (или закачку),
     
  • вводит понятие мультифазной проводимости, сжимаемости, гидропроводности и пьезопроводности пласта, 
     
  • а также от чего именно и в какой степени они зависят

Все это делает модель Перрина полезным аналитическим и учебно-методическим инструментом. 

При численных оценках надо помнить, что чем выше содержание легкой нефти и газа в пласте и чем больше перепад давления в течении исследуемого интервала времени, тем хуже будет оценка давления по модели Перрина. 

Более аккуратные оценки давления при мультифазной фильтрации можно получить на основе:
 

...



Despite the fact that in many practical cases the Linear Perrine multi-phase diffusion model leads to inaccurate pressure predictions it still:

  • provides fair understanding on how pressure responds to multi-phase intakes to or offtakes from formation

  • introduce a concept of multiphase transmissibility, diffusivity, compressibility, mobility and total sandface flowrate which are very helpful in multiphase dynamic analysis

  • provides fair understanding on how the above properties depend on reservoir porosity, permeability, saturation, compressibility and PVT/SCAL model

This makes Linear Perrine multi-phase diffusion model a helpful analytical and  methodological tool for multiphase dynamic analysis.


One should remember that high content of light oil and gas in reservoir and high drawdowns deteriorate the accuracy of Linear Perrine multi-phase diffusion model.


More accurate pressure estimations are provided by the following pressure diffusion models:
 

See also

...

Petroleum Industry / Upstream / Subsurface E&P Disciplines / Well Testing

Pressure diffusion models ] [ Linear Perrine multi-phase diffusion @model derivation ]


Reference

...

...


  1. Anchor
    Perrine
    Perrine
  2. Perrine
  3.  Perrine, R.L. 1956. Analysis of Pressure Buildup Curves. Drill. and Prod. Prac., 482. Dallas, Texas: API.

  4. Anchor
  5.  
  6. Martin
    Martin
     SPE-1235-G, Martin, J.C.
    1959. , Simplified Equations of Flow in Gas Drive Reservoirs and the Theoretical Foundation of Multiphase Pressure Buildup Analyses. SPE-1235-G. Trans., AIME, 216: 321–323.
     Пример оформления ОФП.xls, 1959


Show If
grouparax


Panel
bgColorpapayawhip
titleARAX

Perrine.xls