Specific case of general multi-phase pressure diffusion assuming the equivalent single-phase diffusion with constant dynamic parameters, thus resulting in linear partial differential equation:
\phi c_t \partial_t p - \nabla \big( M \cdot ( \nabla p - \rho \cdot \mathbf{g} ) \big) = q_t(\mathbf{r}) \delta(\mathbf{r}) |
where
q_t(\mathbf{r}) = q_w + q_o + q_g = B_w \, q_W + (B_o - R_v \, B_g) \, q_O + (B_g - R_s \, B_o) \, q_G |
| |
| |
| |
s(\mathbf{r}) = \{ s_w(\mathbf{r}), \ s_o(\mathbf{r}), \ s_g(\mathbf{r}) \} |
|
reservoir saturation as a function of location 
|
c_t = c_r + c_w s_w + c_o s_o + c_g s_g + s_o [ R_{sp} + (c_r + c_o) R_{sn} ] + s_g [ R_{vp} + R_{vn}(c_r + c_g) ] |
| |
| |
| |
M = M_w + M_o \big( 1 + R_{sn} \big) + M_g \big( 1 + R_{vn} \big) |
| |
| |
| |
| |
M_{rw} = \frac{k_{rw}(s)}{\mu_w} |
| |
M_{ro} = \frac{k_{ro}(s)}{\mu_o} |
| |
M_{rg} = \frac{k_{rg}(s)}{\mu_g} |
| |
|
absolute permeability as a function of location at reference pressure 
|
| |
R_{sn} = \frac{R_s B_g}{B_o} \ , \quad R_{vn} = \frac{R_v B_o}{B_g} |
| |
R_{sp} = \frac{\dot R_s B_g}{B_o} \ , \quad R_{vp} = \frac{\dot R_v B_o}{B_g} |
| |
\rho = \frac{ M_{rw} \rho_w + M_{ro} (1 + R_{sn}) \rho_o + M_{rg} (1+R_{vn}) \rho_g }{ M_{rw} + M_{ro} (1 + R_{sn}) + M_{rg} (1+R_{vn}) }
|
| |
g = 9.81 \ \textrm{m} / \textrm{s}^2 |
|
standard gravity |
\big ( \big)^{\LARGE \cdot} = \frac{d}{dp} |
|
differentiation with respect to the pressure
|
All the above dynamic properties are calculated at reference pressure and temperature 
thus making 
a linear partial differential equation.
The reference temperature is more or less in practical cases.
The choice of the reference pressure 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:
p_{\rm ref} = p_{wf}(t=0) |
For accurate modelling of late time pressure response (LTR) it is recommended to use current formation pressure:
Все вышеприведенные параметры рассчитываются при одном давлении (которое предполагается не сильно меняющимся во время исследования), как правило при пластовом давлении.
При небольших депрессиях (репрессиях) это условие можно считать разумным.
Вышеприведенные формулы 
– 
представляют собой обобщение оригинальной модели ( [1], [2] ) на случай летучей нефти.
Для случае черной нефти (
) некторые из вышеприведенных формул упрощаются:
q_t = B_w \, q_W + B_o \, q_O + B_g \, ( q_G - R_s \, q_O) |
| суммарный отбор воды, нефти и газа в пластовых условиях |
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 |
| эффективная сжимаемость пласта с мультифазным насыщением как функция насыщенности при опорном давлении  |
\alpha(s) = \Big \langle \frac{k} {\mu} \Big \rangle = \alpha_w(s) + \alpha_o(s) \big( 1 + R_{sn} \big) + \alpha_g(s) |
| эффективная проводимость пород как функция насыщенности при опорном давлении |
\rho_{\alpha} = \frac{ \alpha_{rw} \rho_w + \alpha_{ro} \rho_o (1 + R_{sn}) + \alpha_{rg} \rho_g }{ \alpha_{rw} + \alpha_{ro} (1 + R_{sn}) + \alpha_{rg} }
|
| гравитационная компонента потока при опорном давлении |
Для режима двухфазной фильтрации недонасыщенной нефти (а именно в этом случае модель Перрина дает наиболее удовлетворительное количественное приближение) вышеприведенные формулы упрощаются еще больше:
q_t = B_w \, q_W + B_o q_O |
| суммарный отбор воды, нефти и газа в пластовых условиях |
c_t(s) = c_r + c_w s_w + c_o s_o |
| эффективная сжимаемость пласта с мультифазным насыщением как функция насыщенности при опорном давлении |
\alpha(s) = \Big \langle \frac{k} {\mu} \Big \rangle = \alpha_w(s) + \alpha_o(s) |
| эффективная проводимость пород как функция насыщенности при опорном давлении |
\rho_{\alpha} = \frac{ \alpha_{rw} \rho_w + \alpha_{ro} \rho_o }{ \alpha_{rw} + \alpha_{ro}}
|
| гравитационная компонента потока при опорном давлении |
\sigma = \Big \langle \frac{k} {\mu} \Big \rangle \, h = k \, h\, \Bigg[ \frac{k_{rw}}{\mu_w} + \frac{k_{ro}}{\mu_o} \Bigg] |
|
гидропроводность пласта при опорном давлении |
Несмотря на то, что приближение Перрина на практике чаще всего приводит к недостаточно аккуратным количественным оценкам давления оно, тем не менее:
- дает правильное представление о характере отклика давления на мультифазные отборы (или закачку),
- вводит понятие мультифазной проводимости, сжимаемости, гидропроводности и пьезопроводности пласта,
- а также от чего именно и в какой степени они зависят
Все это делает модель Перрина полезным аналитическим и учебно-методическим инструментом.
При численных оценках надо помнить, что чем выше содержание легкой нефти и газа в пласте и чем больше перепад давления в течении исследуемого интервала времени, тем хуже будет оценка давления по модели Перрина.
Более аккуратные оценки давления при мультифазной фильтрации можно получить на основе:
See also
Linear Perrine multi-phase diffusion @model derivation
Reference
- Perrine, R.L. 1956. Analysis of Pressure Buildup Curves. Drill. and Prod. Prac., 482. Dallas, Texas: API.
- 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