The general form of the Linear single-phase pressure diffusion @model with the finite number of sources/sinks $\begin{array}{l}\Sigma_k\end{array}$ is given by:

(1)	$\begin{array}{l}\displaystyle \phi \cdot c_t \cdot \partial_t p + \nabla {\bf u} = \sum_k q_k(t) \cdot \delta({\bf r}-{\bf r}_k)\end{array}$

(2)	$\begin{array}{l}\displaystyle {\bf u} = - M \cdot ( \nabla p - \rho \, {\bf g})\end{array}$

(3)	$\begin{array}{l}\displaystyle \int_{\Gamma} \, {\bf u} \, d {\bf \Sigma} = q_\Gamma(t)\end{array}$

or

(4)	$\begin{array}{l}\displaystyle \phi \cdot c_t \cdot \partial_t p + \nabla {\bf u} = 0\end{array}$

(5)	$\begin{array}{l}\displaystyle {\bf u} = - M \cdot ( \nabla p - \rho \, {\bf g})\end{array}$

(6)	$\begin{array}{l}\displaystyle \int_{\Sigma_k} \, {\bf u} \, d {\bf \Sigma} = q_k(t)\end{array}$

(7)	$\begin{array}{l}\displaystyle \int_{\Gamma} \, {\bf u} \, d {\bf \Sigma} = q_\Gamma(t)\end{array}$

where

$\begin{array}{l}p(t, {\bf r})\end{array}$	reservoir pressure	$\begin{array}{l}t\end{array}$	time
$\begin{array}{l}\rho({\bf r})\end{array}$	fluid density	$\begin{array}{l}{\bf r }\end{array}$	position vector
$\begin{array}{l}\phi({\bf r})\end{array}$	effective porosity	$\begin{array}{l}{\bf r }_k\end{array}$	position vector of the $\begin{array}{l}k\end{array}$ -th source
$\begin{array}{l}c_t({\bf r})\end{array}$	total compressibility	$\begin{array}{l}\delta ( \bf r )\end{array}$	Dirac delta function
$\begin{array}{l}M({\bf r})\end{array}$	reservoir fluid mobility $\begin{array}{l}M({\bf r}) = \frac{k({\bf r})}{\mu}\end{array}$	$\begin{array}{l}\nabla\end{array}$	gradient operator
$\begin{array}{l}k({\bf r})\end{array}$	formation permeability to a given fluid	$\begin{array}{l}{ \bf g }\end{array}$	gravity vector
$\begin{array}{l}\mu\end{array}$	dynamic viscosity of a given fluid	$\begin{array}{l}{ \bf u }\end{array}$	fluid velocity under Darcy flow
$\begin{array}{l}q_k(t)\end{array}$	sandface flowrates of the $\begin{array}{l}k\end{array}$ -th source	$\begin{array}{l}\Gamma\end{array}$	reservoir boundary
$\begin{array}{l}q_\Gamma(t)\end{array}$	flow through the reservoir boundary $\begin{array}{l}\Gamma\end{array}$ , which is aquifer or gas cap

Derivation

Physical models of pressure diffusion can be split into two categories: Newtonian and Rheological (non-Newtonian) based on the fluid stress model.

Mathematical models of pressure diffusion can be split into three categories: Linear, Pseudo-Linear and Non-linear.

These models are built using Numerical, Analytical or Hybrid pressure diffusion solvers.

Many popular 1DR solutions can be approximated by Radial Flow Pressure Diffusion @model which has a big methodological value.

The simplest analytical solutions for pressure diffusion are given by 1DL Linear-Drive Solution (LDS) and 1DR Line Source Solution (LSS)

The table below shows a list of popular well and reservoir pressure diffusion models.

Wellbore storage model	Well model	Reservoir model	Boundary model
Constant	Skin-factor	Homogeneous	Infinite
Fair	Vertical well	Dual-porosity	Circle No Flow
Rate-dependant	Fractured vertical well	Dual-permeability	Circle Constant P_i
	Limited entry well	Anisotropic reservoir	Single fault
	Horizontal well	Multi-layer reservoir	Parallel faults
	Slanted well	Linear-composite	Intersecting Faults
	Multifrac horizontal well	Radial-composite

Page tree

See also