Page tree

Versions Compared

Key

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

...

Although the actual reservoir fluid flow may not have an axial symmetry around the well-reservoir contact or around reservoir inhomogeneities (like boundary and faults and composite areas) but still in many practical cases the reservoir flow tends to become radial after some time which makes a Radial Flow Pressure Diffusion @model (in its general form or in particular BVP solution) a popular diagnostic tool. 


Inputs & Outputs

...


InputsOutputs

LaTeX Math Inline
bodyq_t

total sandface rate

LaTeX Math Inline
bodyp(t,r)

reservoir pressure

LaTeX Math Inline
body{p_i}

initial formation pressure

LaTeX Math Inline
body{p_{wf}(t)}

well bottomhole pressure

LaTeX Math Inline
body\sigma

transmissibility,

LaTeX Math Inline
body\sigma = \frac{k \, h}{\mu}



LaTeX Math Inline
body\chi

pressure diffusivity,

LaTeX Math Inline
body\chi = \frac{k}{\mu} \, \frac{1}{\phi \, c_t}



LaTeX Math Inline
bodyS

skin-factor

LaTeX Math Inline
bodyr_w

wellbore radius

LaTeX Math Inline
bodyr_e

drainage radius (could be infinite)


Expand
titleDetailing


LaTeX Math Inline
bodyk

absolute permeability

LaTeX Math Inline
bodyc_t

total compressibility,

LaTeX Math Inline
bodyc_t = c_r + c

LaTeX Math Inline
bodyh

effective thickness

LaTeX Math Inline
body{c_r}

pore compressibility

LaTeX Math Inline
body\mu

dynamic fluid viscosity

LaTeX Math Inline
bodyc

fluid compressibility

LaTeX Math Inline
body{\phi}

porosity




Physical Model

...

Radial fluid flowHomogenous reservoirInfinite boundarySlightly compressible fluid flowConstant rateConstant skin

LaTeX Math Inline
bodyp(t, {\bf r}) \rightarrow p(t, r)

LaTeX Math Inline
body{\bf r} \in ℝ^2 = \{ x, y\}

LaTeX Math Inline
bodyM(r, p)=M =\rm const

LaTeX Math Inline
body\phi(r, p)=\phi =\rm const

LaTeX Math Inline
bodyh(r)=h =\rm const

LaTeX Math Inline
bodyc_r(r)=c_r =\rm const

LaTeX Math Inline
bodyr \rightarrow \infty

LaTeX Math Inline
bodyr_w = 0

LaTeX Math Inline
bodyc_t(r,p) = \rm const

LaTeX Math Inline
bodyq_t = \rm const

LaTeX Math Inline
bodyS = \rm const


Mathematical Model

...


Expand
titleDefinition



LaTeX Math Block
anchorP1
alignmentleft
r_{wf} < r \leq r_e



LaTeX Math Block
anchor52112
alignmentleft
\frac{\partial p}{\partial t}  = \chi \, \left( \frac{\partial^2 p}{\partial r^2} + \frac{1}{r} \frac{\partial p}{\partial r} \right)



LaTeX Math Block
anchor88AEG
alignmentleft
p(t = 0, {\bf r}) = p_i



LaTeX Math Block
anchor3MUX9
alignmentleft
p(t, r \rightarrow r_e ) = p_i)

or

LaTeX Math Block
anchor3MUX9
alignmentleft
\left[ \frac{\partial p}{\partial r} \right]_{r =r_e} = 0



LaTeX Math Block
anchorEM415
alignmentleft
\left[ r\frac{\partial p(t, r )}{\partial r} \right]_{r \rightarrow r_w} = \frac{q_t}{2 \pi \sigma}



LaTeX Math Block
anchorPE
alignmentleft
p_{wf}(t)= p(t,r_w) - S \cdot r_w \, \frac{\partial p}{\partial r} \Bigg|_{r=r_w}




Expand
titleSolution

There is no universal analytical solution to the above problem

LaTeX Math Block Reference
anchorP1
LaTeX Math Block Reference
anchorPE
but it can be always presented as below:


LaTeX Math Block
anchorp_F
alignmentleft
p(t,r) = p_i - \frac{q_t}{4 \pi \sigma} \,  F \bigg( - \frac{r^2}{4 \chi t} \bigg)



LaTeX Math Block
anchorpwf
alignmentleft
p_{wf}(t) = p_i - \frac{q_t}{4 \pi \sigma} \, \bigg[2S +   F \bigg( - \frac{r_w^2}{4 \chi t} \bigg) \bigg]


where

LaTeX Math Inline
bodyF(\xi)
 is a single-argument function describing the peculiarities of the diffusion model (well geometry, penetration geometry, formation inhomogeneities, hydraulic fractures, boundary conditions, etc.).

The fact that solution of equations

LaTeX Math Block Reference
anchorP1
LaTeX Math Block Reference
anchorPE
can be presented as
LaTeX Math Block Reference
anchorp_F
LaTeX Math Block Reference
anchorpwf
finds a lot of practical applications in Well Testing.


...

Expand
titlePTA


PTA – Pressure Transient Analysis



Pressure Drop


LaTeX Math Block
anchor1EWTY
alignmentleft
\delta p = p_i - p_{wf}(t) \sim  \ln t + {\rm const}


Image Modified


Log derivative


LaTeX Math Block
anchorIBA4M
alignmentleft
t \frac{d (\delta p)}{dt}  \sim \rm const







Fig. 2. PTA Diagnostic plot for radial fluid flow




Expand
titleProductivity Index Analysis


The instantaneous Total Sandface Productivity Index for low-compressibility fluid and low-compressibility rocks  does not depend on formation pressurebottomhole pressure and the flowrate and can be expressed as:

LaTeX Math Block
anchorJ
alignmentleft
J_t(t) = \frac{q_t}{p_i - p_{wf}(t)} =\frac{ 2 \pi \sigma }{ S - 0.5 \, F \left( - \frac{r_w^2}{4 \chi t} \right)  }


...