We start with
LaTeX Math Block Reference |
---|
anchor | rho_dif |
---|
page | Derivation of Single-phase pressure diffusion @model |
---|
|
outside with reservoir pressure diffusion 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)0 |
| LaTeX Math Block |
---|
| \int_{\Sigma_k} \, {\bf u} \, d {\bf A} = q_k(t) |
|
...
| |
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 Then use the following equality:
...
...
|
d(\rho \, \phi) = \rho \, d \phi + \phi \, d\rho = \rho \, \phi \, \left( \frac{d \phi }{\phi} + \frac{d \rho }{\rho} \right)
= \rho \, \phi \, \left( \frac{1}{\phi} \frac{d \phi}{dp} \, dp + \frac{1}{\rho} \frac{d \rho}{dp} \, dp \right)
= \rho \, \phi \, (c_{\phi} \, dp + c \, dp) = \rho \, \phi \, c_t \, dp |
to
page | Derivation of Single-phase 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 |
---|
| \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} \, p \cdot d {\bf A} = q_k(t) |
|
Let's express the density via Z-factor:
...
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 p\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} \, 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.
...
LaTeX Math Block |
---|
|
\tau(t) = \int_0^t \frac{dt}{\mu (p_{BHP}(t) ) \, c_t (p_{BHP}) } \, , \ \ p_{BHP} = 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
...