...
We start with
LaTeX Math Block Reference |
---|
anchor | rho_dif |
---|
page | Derivation of Single-phase pressure diffusion @model |
---|
|
outside wellbore:
LaTeX Math Block |
---|
anchor | rho_dif |
---|
alignment | left |
---|
| \frac{\partial (\rho \phi)}{\partial t} + \nabla \, ( \rho \, {\bf u}) = \rho \cdot \sum_k q_k(t) \cdot \delta({\bf r}-{\bf r}_k) |
|
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 Single-phase pressure diffusion @model |
---|
|
to arrive at:
LaTeX Math Block |
---|
| \rho \, 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 and 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 \, 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 p} \, 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:
\rho \, 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 p} \, d {\bf A} = q_k(t) |
|
or
LaTeX Math Block |
---|
| 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 p} \, 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.
See also
...
Physics / Mechanics / Continuum mechanics / Fluid Mechanics / Fluid Dynamics / Pressure Diffusion / Pressure Diffusion @model / Pseudo-linear pressure diffusion @model
...