We start with reservoir pressure diffusion outside wellbore:
LaTeX Math Block |
---|
anchor | rho_dif |
---|
alignment | left |
---|
| \frac{\partial (\rho \phi)}{\partial t} + \nabla \, ( \rho \, {\bf u}) = 0 |
| LaTeX Math Block |
---|
| \int_{\Sigma_k} \, {\bf u} \, d {\bf A} = q_k(t) |
|
where
| |
LaTeX Math Inline |
---|
body | --uriencoded--d %7B\bf \Sigma%7D |
---|
|
| normal vector of differential area on the well-reservoir contact, pointing inside wellbore |
and use
LaTeX Math Block Reference |
---|
anchor | din_term |
---|
page | Derivation of Linear pressure diffusion @model |
---|
|
to arrive at:
LaTeX Math Block |
---|
| \rho \, \phi \, c_t \cdot \frac{\partial (p)}{\partial t} + \nabla \, ( \rho \, {\bf u}) = 0 |
| LaTeX Math Block |
---|
| \int_{\Sigma_k} \, {\bf u} \, d {\bf A} = q_k(t) |
|
where
Let's assume Darcy flow with constant permeability
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \frac%7Bdk%7D%7Bdp%7D = 0 |
---|
|
and ignore gravity forces:
LaTeX Math Block |
---|
|
{\bf u} = \frac{k}{\mu} \nabla \, p |
so that diffusion equation becomes:
LaTeX Math Block |
---|
| \rho \, \phi \, c_t \cdot \frac{\partial p}{\partial t} + \nabla \, ( k \cdot \frac{\rho}{\mu} \, \nabla \, p) = 0 |
| LaTeX Math Block |
---|
| \frac{k}{\mu} \cdot \int_{\Sigma_k} \, {\bf \nabla } \, p \cdot d {\bf A} = q_k(t) |
|
Let's express the density via Z-factor:
LaTeX Math Block |
---|
|
\rho = \frac{M}{RT} \, \frac{p}{Z(p)} |
where
and assuming the fluid temperature
does not change over time and space during the modelling period:
LaTeX Math Block |
---|
| \phi \, c_t \, \mu \cdot \frac{p}{\mu \, Z} \cdot \frac{\partial p}{\partial t} + \nabla \, ( k \cdot \frac{p}{\mu \, Z} \, \nabla \, p) = 0 |
| LaTeX Math Block |
---|
| \frac{k}{\mu} \cdot \int_{\Sigma_k} \, {\bf \nabla } \, p \cdot d {\bf A} = q_k(t) |
|
or
LaTeX Math Block |
---|
| \phi \, c_t \, \mu \cdot \frac{\partial \Psi}{\partial t} + \nabla \, ( k \cdot \nabla \, \Psi) = 0 |
| LaTeX Math Block |
---|
| \frac{k}{\mu} \cdot \int_{\Sigma_k} \, {\bf \nabla } \, p \cdot d {\bf A} = q_k(t) |
|
where
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \Psi(p) =2 \, \int_0%5ep \frac%7Bp \, dp%7D%7B\mu(p) \, Z(p)%7D |
---|
|
| Pseudo-Pressure |
In some practical cases the complex
can be considered as constant in time which makes
LaTeX Math Block Reference |
---|
|
a linear differential equation.
But during the early transition times the pressure drop is usually high and the complex
can not be considered as constant in time which leads to distortion of
pressure transient diagnostics at early times.
In this case one can use Pseudo-Time, calculated by means of the bottom-hole pressure
LaTeX Math Inline |
---|
body | --uriencoded--p_%7BBHP%7D(t) |
---|
|
:
LaTeX Math Block |
---|
|
\tau(t) = \int_0^t \frac{dt}{\mu(p_{BHP}(t)) \, c_t(p_{BHP}(t))}
|
to correct early-time transient behaviour which turn equation
LaTeX Math Block Reference |
---|
|
into:
LaTeX Math Block |
---|
|
\phi \cdot \frac{\partial \Psi}{\partial \tau} + \nabla \, ( k \cdot \nabla \, \Psi) = 0 |
See also
Physics / Mechanics / Continuum mechanics / Fluid Mechanics / Fluid Dynamics / Pressure Diffusion / Pressure Diffusion @model / Pseudo-linear pressure diffusion @model