Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

Motivation

...


Excerpt Include
Aquifer Drive
Aquifer Drive
nopaneltrue

Inputs & Outputs

...


InputsOutputs

LaTeX Math Inline
bodyp(t)

field-average formation pressure at time moment

LaTeX Math Inline
bodyt

LaTeX Math Inline
bodyQ^{\downarrow}_{AQ}(t)

Cumulative subsurface water influx from aquifer

LaTeX Math Inline
bodyp_i

initial formation pressure

LaTeX Math Inline
bodyq^{\downarrow}_{AQ}(t) = \frac{dQ^{\downarrow}_{AQ}}{dt}

Subsurface water flowrate from aquifer

LaTeX Math Inline
bodyB

water influx constant

LaTeX Math Inline
body\chi

aquifer diffusivity

LaTeX Math Inline
bodyA_e

net pay area
Expand
titleDetailing
Detailing Inputs

LaTeX Math Inline
bodyB = \frac{\theta}{2\pi} \cdot A_e \cdot h

...

\cdot \phi

...

\cdot c_t

water influx constant

LaTeX Math Inline
body\theta

central angle of net pay area ↔ aquifer contact

LaTeX Math Inline
bodyh

...

aquifer effective thickness

LaTeX Math Inline
body\phi

...

aquifer porosity

LaTeX Math Inline
bodyc_t=c_r +c_w

aquifer total compressibility

LaTeX Math Inline
bodyc_r

aquifer pore compressibility 

LaTeX Math Inline
bodyc_w

aquifer water compressibility

...

Physical Model

...


...









Fig. 1.Radial VEHaquifer drive schematic


Mathematical Model

...

LaTeX Math Block
anchor
1
VEH
alignmentleft
\frac{d
Q^{\downarrow}_{AQ}
}{dt} = q^{\downarrow}_{AQ}(t)
= B \cdot \int_0^t W_{eD} \left( \frac{(t-\tau)\chi}{r_e^2}\right) \dot p(\tau) d\tau
LaTeX Math Block
anchor
1
WeD
alignmentleft
q^{\downarrow}
W_{
AQ
eD}(t) = 
C_a \cdot
\int_0^{t} \frac{\partial p_
a(t,r)
1}{\partial r_D} \bigg|_{r
=r_e}
_D = 1} dt_D 
LaTeX Math Block
anchor
CT
1
alignmentleft
p_a
q^{\downarrow}_{AQ}(t
, r
)= 
p(0) +
\
int_0^t p_1 \left(\
frac{
(t-\tau) \cdot \chi}{r_e^2}, \frac{r}{r_e} \right) \dot p(\tau) d\tau
dQ^{\downarrow}_{AQ}}{dt}
LaTeX Math Block
anchor1
alignmentleft
p_1 = p_1(t_D, r_D)
LaTeX Math Block
anchor
RC
RC1
alignmentleft
\frac{\partial p_1}{\partial t_D} =  \frac{\partial^2 p_1}{\partial r_D^2} + \frac{1}{r_D}\cdot \frac{\partial p_1}{\partial r_D}
LaTeX Math Block
anchorCT
alignmentleft
p_1(t_D = 0, r_D)= 0



LaTeX Math Block
anchorCT
alignmentleft
p_1(t_D, r_D=1) = 1
LaTeX Math Block
anchor
1
gradp
alignmentleft
\frac{\partial p_1(t_D, r_D)}{\partial r_D} 
\
bigg
Bigg|_
{r_D=r
_a/r_e}
_{aD}} = 0

or

LaTeX Math Block
anchorgradp
alignmentleft
 p_1(t_D, r_D = \infty) = 0



Expand
titleDerivation
Derivation of Radial VEH Aquifer Drive @model


Computational Model

...

LaTeX Math Block
anchor

...

VEHD
alignmentleft
Q^{\downarrow}_{AQ}(t)= B \cdot \sum_\alpha W_{eD} 
\left( \frac{ (t-\tau_\alpha) \chi}{

...

r_e^2} \right)\Delta p_\alpha 


= B \cdot W_{eD} 
\left( \frac{ (t-\tau_1) \chi}{r_e^2}  \right)\Delta p_1 +
 B \cdot W_{eD} 
\left( \frac{ (t-\tau_2) \chi}{r_e^2}  \right)\Delta p_2
+ ... + B \cdot W_{eD} 
\left( \frac{ (t-\tau_N) \chi}{r_e^2} \right)\Delta p_N



This computational model is using a discrete convolution (also called superposition in some publications) with time-grid 

LaTeX Math Inline
body\{ \tau_1, \, \tau_2, \ ... \ , \tau_N \}
.

In practical exercises with manual or spreadsheet-assisted calculations the time-grid is usually uniform: 

LaTeX Math Inline
body\{ \tau_1 =\Delta \tau, \, \tau_2 = 2 \cdot \Delta \tau, \ ... \ , \tau_N = N \cdot \Delta \tau\}
 with the time step 
LaTeX Math Inline
body\Delta \tau
 of a month to ensure the formation pressure does not change much since the previous time step.

Moving to annual time step may accumulate a substantial mistake if formation pressure has varied substantially in some years.  


Expand
titlePolynomial approximations for WeD

Polynomial approximation of 

LaTeX Math Inline
bodyW_{eD}(t_D)
  are available (http://dx.doi.org/10.2118/15433-PA).


Table 1. Polynomial approximation of 

LaTeX Math Inline
bodyW_{eD}(t_D)
 for infinite aquifer

LaTeX Math Inline
bodyt_D< 0.01

LaTeX Math Inline
bodyW_{eD}=\sqrt{\frac{t_D}{\pi}}

LaTeX Math Inline
body0.01 < t_D<200

LaTeX Math Inline
body\displaystyle W_{eD}=\frac {1.2838 \cdot t_D^{1/2} + 1.19328 \cdot t_D +0.269872 \cdot t_D^{3/2} +0.00855294 \cdot t_D^2} {1+0.616599 \cdot t_D^{1/2}+0.0413008 \cdot t_D}

LaTeX Math Inline
bodyt_D > 200

LaTeX Math Inline
body\displaystyle W_{eD}=\frac{-4.29881+2.02566 \cdot t_D}{\ln t_D}


Scope of Applicability

...


The benefit of VEH approach is that net pay formation pressure history 

LaTeX Math Inline
bodyp(t)
is usually known so that water influx calculation based on aquifer properties 
LaTeX Math Inline
body\{ B, \, r_a, \, \chi \}
 is rather straightforward.

In the past the VEH approach was considered as tedious in calculating superposition during the manual exercises.

In modern computers the convolution is a fast fully-automated procedure and VEH model is considered as a reference in the range of analytical aquifer models.

Although the model is derived for linear and radial flow it also shows a good match for bottom-water drive depletions.

See Also

...

Petroleum Industry / Upstream / Subsurface E&P Disciplines / Field Study & Modelling / Aquifer Drive / Aquifer Drive Models

Reference

...

 1. van Everdingen, A.F. and Hurst, W. 1949. The Application of the Laplace Transformation to Flow Problems in Reservoirs. Trans., AIME 186, 305.

...