We start with reservoir pressure diffusion outside wellborethe reservoir flow continuity equation:
LaTeX Math Block |
---|
anchor | rho_dif |
---|
alignment | left |
---|
| \frac{\partial (\rho \phi)}{\partial t} + \nabla \, ( \rho \, {\bf u}) = \sum_k \dot m_k(t) \cdot \delta({\bf r}-{\bf r}_k) |
|
percolation model:
| qk | \int_ | Sigma_k | \,= - M \cdot ( \nabla p - \rho |
| (p) | u} \, d |
and the reservoir boundary flow condition:
LaTeX Math Block |
---|
anchor | qGamma |
---|
alignment | left |
---|
| {\rm F}_{\Gamma}(p, {\bf |
| A | m_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 |
| mass rate flowrate at -th well LaTeX Math Inline |
---|
body | \dot m_k(t) = \rho(p) \cdot q_k(t) |
---|
|
|
| |
| |
...
LaTeX Math Block |
---|
anchor | rhophi |
---|
alignment | left |
---|
|
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 |
where
to arrive at:
S8TNB | \rho \, \phi \, c_t \cdot \frac{\partial p}{\partial t} + \nabla \, ( \rho \, {\bf u}) = \sum_k \rho \, |
|
\sum_k q_k(t) \cdot \delta({\bf r}-{\bf r}_k) |
| |
qkint_{\Sigma_k} \rm F}_{\Gamma}(p, {\bf u}) = 0 |
|
The left-hand side of equation
LaTeX Math Block Reference |
---|
|
can be transformed in the following way: LaTeX Math Block |
---|
anchor | nabla_rho |
---|
alignment | left |
---|
|
\nabla \, ( \rho \, {\bf u}) = \rho \, \nabla \, {\bf u} + (\nabla \rho, \, {\bf u}) = \rho \, \nabla \, {\bf u} + \frac{d\rho}{dp} \cdot (\nabla p, \, {\bf |
...
...
\rho \, \nabla \, {\bf u} + \rho \, c \cdot (\nabla p, \, {\bf u}) |
...
where
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle c |
---|
|
...
(p) = \frac%7B1%7D%7B\rho%7D \frac%7Bd\rho%7D%7Bdp%7D |
|
is fluid compressibility.By using the Dirac delta function property:
LaTeX Math Inline |
---|
body | f(x) \cdot \delta(x-x_0) = f(x_0) \cdot \delta(x-x_0) |
---|
|
the right-hand side of equation LaTeX Math Block Reference |
---|
|
can be transformed in the following wayThen use the following equality:
LaTeX Math Block |
---|
anchor | 1right_hand_side |
---|
alignment | left |
---|
|
\nablasum_k \rho(p(t, {\bf r})) \cdot q_k(t) \cdot \delta({\bf r}-{\bf r}_k)
= \sum_k \rho(p(t, {\bf r}_k)) \cdot q_k(t) \cdot \delta({\bf r}-{\bf r}_k)
= \sum_k \rho(p(t, {\bf u}) = ( \rho \, \nabla \,r})) \cdot q_k(t) \cdot \delta({\bf r}-{\bf r}_k)
= \rho(p) \cdot \sum_k q_k(t) \cdot \delta({\bf r}-{\bf r}_k) |
Substituting
LaTeX Math Block Reference |
---|
|
and LaTeX Math Block Reference |
---|
|
into LaTeX Math Block Reference |
---|
|
and reducing the density one arrives to: LaTeX Math Block |
---|
anchor | pre_filnal |
---|
alignment | left |
---|
| \phi \, c_t \cdot \frac{\partial p}{\partial t} + \nabla {\bf u}
+ c \cdot ( {\bf u} |
|
...
\, \nabla p) = \sum_k q_k(t) \cdot \delta({\bf r}-{\bf r}_k) |
| LaTeX Math Block |
---|
anchor | qGamma |
---|
alignment | left |
---|
| {\rm F}_{\Gamma}(p, {\bf u}) = 0 |
|
See also
...
Physics / Mechanics / Continuum mechanics / Fluid Mechanics / Fluid Dynamics / Pressure Diffusion / Pressure Diffusion @model / Single-phase pressure diffusion @model
...