Skip to content
Brent$117.6Gas TTF$1,192USD453.87 ₸
MunaiHub

Equations: Boundaries and late-time behaviour

The equations of the “Boundaries and late-time behaviour” lessons in “Well testing and pressure transient analysis”, with their conditions and units, and what each symbol means in them.

A well near a sealing fault: dimensionless form
pD=12(ln⁡tD+0.809)+S⏟well+12 E1 ⁣(LD2tD)⏟image,LD=Lrwp_D = \underbrace{\tfrac12\left(\ln t_D + 0.809\right) + S}_{\text{well}} + \underbrace{\tfrac12\,E_1\!\left(\frac{L_D^2}{t_D}\right)}_{\text{image}},\qquad L_D = \frac{L}{r_w}

no units; natural logarithm

Holds when: The fault is a straight sealing boundary through the whole thickness of a homogeneous reservoir, and the disturbance has not reached any other boundaries. For the well itself the logarithmic approximation holds (tD>25t_D > 25), and wellbore storage is over. The image is a line source at a distance 2L2L. The argument of E1E_1 is (2LD)2/(4tD)=LD2/tD(2L_D)^2/(4t_D) = L_D^2/t_D.

A sealing fault near the well: the method of images →
The first and second semilog straight lines
first:Δp=me[ln⁡tD+0.809+2S]second:Δp=2me[ln⁡tD+0.809+S−ln⁡2Lrw]\begin{aligned}\text{first:}\quad &\Delta p = m_e\left[\ln t_D + 0.809 + 2S\right]\\ \text{second:}\quad &\Delta p = 2m_e\left[\ln t_D + 0.809 + S - \ln\frac{2L}{r_w}\right]\end{aligned}

any consistent system of units; natural logarithm

Holds when: The first line holds while ηt/L2<0.2\eta t/L^2 < 0.2, once wellbore storage is over. The second line holds at ηt/L2≫1\eta t/L^2 \gg 1. Its slope differs from 2m2m by less than 5% only after ηt/L2≈10\eta t/L^2 \approx 10. The skin enters the second line once, because only the real well has a pressure drop across the skin.

A sealing fault near the well: the method of images →
Distance to the fault from the intersection time of the lines
tx=γ ϕμctL2k⇒L=k txγ ϕμct,γ=1.781t_x = \frac{\gamma\,\phi\mu c_tL^2}{k}\quad\Rightarrow\quad L = \sqrt{\frac{k\,t_x}{\gamma\,\phi\mu c_t}},\qquad \gamma = 1.781

SI

Holds when: Both lines are identified reliably, and the second is drawn through points where the slope has really doubled. ϕ\phi and ctc_t are known independently. Their error passes into LL as a square root.

A sealing fault near the well: the method of images →
The derivative step at a sealing fault, and the distance from a single point on it
Δp′=me[1+exp⁡(−L2ηt)]⇒L=η t ln⁡meΔp′−me\Delta p' = m_e\left[1 + \exp\left(-\frac{L^2}{\eta t}\right)\right]\quad\Rightarrow\quad L = \sqrt{\eta\,t\,\ln\frac{m_e}{\Delta p' - m_e}}

SI

Holds when: The reservoir is homogeneous and there is one straight sealing fault. Wellbore storage is over, and mem_e and kk are known from the first plateau. We take points in the middle of the step, where Δp′/me\Delta p'/m_e is between 1.2 and 1.7. At the start and at the end of the step, small errors in the derivative change LL greatly.

A sealing fault near the well: the method of images →
Hemiradial flow: flow capacity and skin for a known distance to the fault
kh=qBμ2π Δp2′,S=Δp(t)Δp2′−ln⁡ktϕμctrw2−0.809+ln⁡2Lrwkh = \frac{qB\mu}{2\pi\,\Delta p_2'},\qquad S = \frac{\Delta p(t)}{\Delta p_2'} - \ln\frac{kt}{\phi\mu c_tr_w^2} - 0.809 + \ln\frac{2L}{r_w}

SI; natural logarithm

Holds when: The first plateau is hidden by wellbore storage, and the doubled plateau has formed. The distance LL is known independently, because we cannot find it from such a record. The reservoir is homogeneous, and the fault is straight and sealing.

A sealing fault near the well: the method of images →
Build-up of a well near a sealing fault
pi−pws=me[ln⁡tp+ΔtΔt+E1 ⁣(L2η (tp+Δt))−E1 ⁣(L2η Δt)]p_i - p_{ws} = m_e\left[\ln\frac{t_p + \Delta t}{\Delta t} + E_1\!\left(\frac{L^2}{\eta\,(t_p + \Delta t)}\right) - E_1\!\left(\frac{L^2}{\eta\,\Delta t}\right)\right]

SI; natural logarithm

Holds when: Before start-up the pressure in the reservoir was pip_i. The rate stays constant during tpt_p. For the well itself the logarithmic approximation holds, and wellbore storage is over. There are no other boundaries.

A sealing fault near the well: the method of images →
Distance to the fault from the intersection of the Horner lines
L=η Δtexγ,Δtex=tp Δtxtp+ΔtxL = \sqrt{\frac{\eta\,\Delta t_{ex}}{\gamma}},\qquad \Delta t_{ex} = \frac{t_p\,\Delta t_x}{t_p + \Delta t_x}

SI

Holds when: The fault was felt long before shut-in: L2≪ηtpL^2 \ll \eta t_p. Then E1(L2/(ηtp))E_1\big(L^2/(\eta t_p)\big) can be replaced by a logarithm. Both Horner lines have fully formed. A second Horner line that has not fully formed gives too small a value of LL, as in a drawdown test.

A sealing fault near the well: the method of images →
Largest deviation of the build-up from the line, and the maximum distance at which a fault can be detected
δpmax⁡=qB tp4πe ϕcth L2⇒Lmax⁡=qB tp4πe ϕcth δp\delta p_{\max} = \frac{qB\,t_p}{4\pi\e\,\phi c_th\,L^2}\quad\Rightarrow\quad L_{\max} = \sqrt{\frac{qB\,t_p}{4\pi\e\,\phi c_th\,\delta p}}

SI

Holds when: The fault was not felt during production: ηtp≪L2\eta t_p \ll L^2. In practice the formula already works for L2>3ηtpL^2 > 3\eta t_p. The record lasts at least L2/ηL^2/\eta after shut-in. The reservoir is homogeneous. There are no other disturbances, such as neighbouring wells in operation or gauge drift (a slow creep of the gauge reading).

A sealing fault near the well: the method of images →
The pressure derivative of a well with several sealing boundaries
Δp′=me[1+∑jexp⁡(−dj24ηt)],me=qBμ4πkh,η=kϕμct\Delta p' = m_e\left[1 + \sum_{j} \exp\left(-\frac{d_j^2}{4\eta t}\right)\right],\qquad m_e = \frac{qB\mu}{4\pi kh},\quad \eta = \frac{k}{\phi\mu c_t}

SI (any consistent system)

Holds when: The boundaries are straight sealing planes through the whole thickness of a homogeneous reservoir. Images describe exactly a channel and a wedge with an angle 180∘/n180^\circ/n. Wellbore storage is over. The images are line sources (djd_j is much larger than the wellbore radius). The rate is constant. For a build-up we add up the same terms by superposition, and we allow for the producing time (Pressure build-up analysis).

Several boundaries: parallel faults and wedge-shaped blocks →
A well in a channel: linear flow
pD=2πtDW+S+Sch,tDW=ktϕμctW2⇔Δp=mLt+qBμ2πkh (S+Sch)p_D = 2\sqrt{\pi t_{DW}} + S + S_{ch},\qquad t_{DW} = \frac{kt}{\phi\mu c_tW^2}\quad\Leftrightarrow\quad \Delta p = m_L\sqrt t + \frac{qB\mu}{2\pi kh}\,(S + S_{ch})

SI

Holds when: The walls are parallel sealing boundaries through the whole thickness of a homogeneous reservoir; along the channel the reservoir is unbounded. The disturbance has reached both walls: the derivative lies on a line with slope ½ (by the table below, at ηt/W2\eta t/W^2 from 0.1 to 0.4, depending on the position of the well). A constant rate; wellbore storage is over.

Several boundaries: parallel faults and wedge-shaped blocks →
Geometric skin of a channel
Sch=ln⁡W2πrw−ln⁡sin⁡πL1WS_{ch} = \ln\frac{W}{2\pi r_w} - \ln\sin\frac{\pi L_1}{W}

any equal units of length

Holds when: Steady linear flow; the well is a line source in a homogeneous channel (a pseudo-skin of the geometry, and not damage of the reservoir).

Several boundaries: parallel faults and wedge-shaped blocks →
The channel width from the slope of the line Δp against √t
W=qBh mLμπkϕctW = \frac{qB}{h\,m_L}\sqrt{\frac{\mu}{\pi k\phi c_t}}

SI

Holds when: The line is drawn through the points where the derivative already has a slope of ½. The permeability was found independently: from the radial flow of the same record or from other data. ϕ\phi, ctc_t and hh are known: their error enters WW as 1/ϕct1/\sqrt{\phi c_t} and 1/h1/h.

Several boundaries: parallel faults and wedge-shaped blocks →
A build-up in a channel in the section of linear flow
pi−pws=mL(tp+Δt−Δt)p_i - p_{ws} = m_L\left(\sqrt{t_p + \Delta t} - \sqrt{\Delta t}\right)

SI (metric oilfield units: MPa, h, mLm_L in MPa/√h)

Holds when: A constant rate during tpt_p; before start-up the pressure in the reservoir was pip_i. Both tpt_p and Δt\Delta t are large enough for both terms of the superposition to be in the linear-flow section (ηΔt/W2\eta\Delta t/W^2 not less than 0.1–0.4 by the table above).

Several boundaries: parallel faults and wedge-shaped blocks →
The late level of the derivative in a wedge, and the angle from it
Δp∞′=360∘θ me⇒θ=360∘ meΔp∞′\Delta p'_{\infty} = \frac{360^\circ}{\theta}\,m_e\quad\Rightarrow\quad \theta = 360^\circ\,\frac{m_e}{\Delta p'_{\infty}}

any equal units of pressure

Holds when: Two straight sealing faults through the whole thickness of a homogeneous reservoir meet at one point (in plan view). The disturbance has reached no other boundaries. The level has really been reached: on the way to it, the derivative is still lower, and the angle from it comes out too large. The angle is less than 180°: the well is inside an “acute” (less than 90°) or “obtuse” (more than 90°) angle, not beyond a bend of a fault.

Several boundaries: parallel faults and wedge-shaped blocks →
The approach of the derivative to the late level of a wedge, and the angle from two late points
Δp′me≈R∞(1−at),a=d22η;a=r1−r2r1/t2−r2/t1,R∞=r21−a/t2\frac{\Delta p'}{m_e} \approx R_\infty\left(1 - \frac{a}{t}\right),\quad a = \frac{d^2}{2\eta};\qquad a = \frac{r_1 - r_2}{r_1/t_2 - r_2/t_1},\quad R_\infty = \frac{r_2}{1 - a/t_2}

SI; tt and aa in the same units of time

Holds when: The shortfall is small, no more than 15–20%. At a larger one, the next terms of the expansion matter, and the estimate of R∞R_\infty comes out too low and the angle too large. The derivative from measurements has no noticeable noise: the difference r2−r1r_2 - r_1 is small, and its error strongly affects the result. The disturbance has reached no other boundaries than the two faults.

Several boundaries: parallel faults and wedge-shaped blocks →
Pressure in a well with images, in dimensionless form
pwD=12(ln⁡tD+0.809)+S+12∑jE1 ⁣(βj24τ),τ=ηtL2=tDLD2,βj=djLp_{wD} = \tfrac12\left(\ln t_D + 0.809\right) + S + \tfrac12\sum_j E_1\!\left(\frac{\beta_j^2}{4\tau}\right),\qquad \tau = \frac{\eta t}{L^2} = \frac{t_D}{L_D^2},\quad \beta_j = \frac{d_j}{L}

dimensionless

Holds when: the boundaries are straight and sealing, through the whole thickness of a homogeneous reservoir; the images are line sources (dj≫rwd_j \gg r_w) the rate is constant, wellbore storage is over, tD>25t_D > 25 for images that go off to infinity (channels), we add terms until they become negligible

Semi-infinite systems: theory and test analysis →
Product solution: the derivative of a rectangular system
Δp′me=Fx(t) Fy(t),Fx=∑iexp⁡ ⁣(−ξi24ηt),Fy=∑lexp⁡ ⁣(−ζl24ηt)\frac{\Delta p'}{m_e} = F_x(t)\,F_y(t),\qquad F_x = \sum_i \exp\!\left(-\frac{\xi_i^2}{4\eta t}\right),\quad F_y = \sum_l \exp\!\left(-\frac{\zeta_l^2}{4\eta t}\right)

SI (any consistent system)

Holds when: all the boundaries are straight sealing planes, parallel to xx or yy, through the whole thickness of a homogeneous reservoir the well is a line source, wellbore storage is over, and the rate is constant

Semi-infinite systems: theory and test analysis →
A well in a channel closed at one end: late linear flow
Δp=2mLt+qBμ2πkh (S+Sch)−qBμkhW L3,Δp′=mLt\Delta p = 2m_L\sqrt t + \frac{qB\mu}{2\pi kh}\,(S + S_{ch}) - \frac{qB\mu}{khW}\,L_3,\qquad \Delta p' = m_L\sqrt t

SI

Holds when: the walls and the closed end are sealing planes through the whole thickness of a homogeneous reservoir; in the other direction the channel is unbounded the closed end is farther away than the width of the channel (L3>WL_3 > W), and ηt≫L32\eta t \gg L_3^2, so the closed end has long been “felt”; wellbore storage is over, and the rate is constant

Semi-infinite systems: theory and test analysis →
Distance to the closed end of the channel: from a point of the derivative and from the intersection of the lines of Δp against √t
L3=ηt ln⁡1r−1,r=Δp′(t)12mLt;L3=η txπL_3 = \sqrt{\eta t\,\ln\frac{1}{r - 1}},\quad r = \frac{\Delta p'(t)}{\tfrac12 m_L\sqrt t};\qquad L_3 = \sqrt{\frac{\eta\,t_x}{\pi}}

SI

Holds when: the record contains the linear flow of the open channel. We draw the line 12mLt\tfrac12m_L\sqrt t along it on the derivative, and the line with slope mLm_L on the plot of Δp\Delta p against t\sqrt t the closed end is farther away than the width of the channel; the permeability is known from the radial-flow plateau the second formula holds if the late line has really been reached. If we draw a line with slope 2mL2m_L through the last reading at time tendt_{end}, it gives Lx=L3−L32/πηtendL_x = L_3 - L_3^2/\sqrt{\pi\eta t_{end}} (the correction is in the example)

Semi-infinite systems: theory and test analysis →
A well at a sealing boundary as a well with a doubled rate
pwD=2 pD∗(tD)−ln⁡2δrw+Sp_{wD} = 2\,p_D^{*}(t_D) - \ln\frac{2\delta}{r_w} + S

dimensionless

Holds when: the boundary is straight and sealing; δ\delta is small compared with the other dimensions of the system the time is longer than a few δ2/η\delta^2/\eta: the image at distance 2δ2\delta has already merged with the well (the derivation is in the advanced section)

Semi-infinite systems: theory and test analysis →
Build-up of a well in a system with images
pi−pws=me[ln⁡tp+ΔtΔt+∑j(E1 ⁣(dj24η (tp+Δt))−E1 ⁣(dj24η Δt))]p_i - p_{ws} = m_e\left[\ln\frac{t_p + \Delta t}{\Delta t} + \sum_j\left(E_1\!\left(\frac{d_j^2}{4\eta\,(t_p + \Delta t)}\right) - E_1\!\left(\frac{d_j^2}{4\eta\,\Delta t}\right)\right)\right]

SI; natural logarithm

Holds when: before start-up the pressure in the formation was pip_i; the rate is constant during tpt_p; for the well itself the logarithmic approximation holds, and wellbore storage is over for channels the sum is infinite but converges: distant images give equal terms at tp+Δtt_p + \Delta t and at Δt\Delta t, and the difference vanishes

Semi-infinite systems: theory and test analysis →
Largest trace of the closed end of a channel on a build-up after a short producing time
δpmax⁡=e−1/222π⋅qB tpϕcth WL3=0.121 qB tpϕcth WL3,Δtmax⁡≈2L32η\delta p_{\max} = \frac{\e^{-1/2}}{2\sqrt{2\pi}}\cdot\frac{qB\,t_p}{\phi c_t h\,WL_3} = 0.121\,\frac{qB\,t_p}{\phi c_t h\,WL_3},\qquad \Delta t_{\max} \approx \frac{2L_3^2}{\eta}

SI

Holds when: during the producing time the closed end has hardly been felt: ηtp≪L32\eta t_p \ll L_3^2 (in practice it is enough that L32>3ηtpL_3^2 > 3\eta t_p) the closed end is farther away than the width of the channel, and by the time Δtmax⁡\Delta t_{\max} the channel has long been in linear flow a homogeneous formation; no gauge drift and no interference from neighbouring wells

Semi-infinite systems: theory and test analysis →
Flow through the fault zone
u=kfμ⋅p1−p2lfu = \frac{k_f}{\mu}\cdot\frac{p_1 - p_2}{l_f}

SI (any consistent system)

Holds when: the zone is thin compared with the distance to the well, and it stores no fluid single-phase flow; the zone is the same along the whole length and thickness of the fault on both sides there is a formation with the same properties and thickness

Partially communicating faults →
Dimensionless fault transmissibility
α=kf/lfk/L=kf Lk lf\alpha = \frac{k_f/l_f}{k/L} = \frac{k_f\,L}{k\,l_f}

dimensionless; kfk_f and kk, and lfl_f and LL, in the same units

Holds when: the thickness of the formation is the same on both sides of the fault (otherwise the numerator and the denominator contain their own khkh) the fluid has the same properties on both sides

Partially communicating faults →
A well near a partially communicating fault: dimensionless pressure
pD=12(ln⁡tD+0.809)+S+12E1 ⁣(1tDL)⏟image−α∫0∞e−2αξ E1 ⁣((2+ξ)24tDL) dξ⏟leakp_D = \tfrac12\left(\ln t_D + 0.809\right) + S + \underbrace{\tfrac12 E_1\!\left(\frac{1}{t_{DL}}\right)}_{\text{image}} - \underbrace{\alpha\int_0^\infty \e^{-2\alpha\xi}\,E_1\!\left(\frac{(2 + \xi)^2}{4t_{DL}}\right)\dd\xi}_{\text{leak}}

dimensionless; natural logarithm

Holds when: the fault is a straight semipermeable barrier through the whole thickness of a homogeneous formation that is the same on both sides; there are no other boundaries for the well itself the logarithmic approximation holds, and wellbore storage is over the image and the leak are line sources (LL is much larger than the well radius)

Partially communicating faults →
Pressure derivative near a partially communicating fault
Δp′me=1+e−1/tDL[1−2π αtDL  erfcx(2αtDL+1tDL)]\frac{\Delta p'}{m_e} = 1 + \e^{-1/t_{DL}}\left[1 - 2\sqrt{\pi}\,\alpha\sqrt{t_{DL}}\;\mathrm{erfcx}\left(2\alpha\sqrt{t_{DL}} + \frac{1}{\sqrt{t_{DL}}}\right)\right]

dimensionless

Holds when: the same as for the pressure; wellbore storage is over.

Partially communicating faults →
The late stretch: the rest of the derivative and the apparent fault skin
Δp′me−1≈1+4α8α2 tDL,Δp=me[ln⁡tD+0.809+2S+2Sf],Sf=e4αE1(4α)\frac{\Delta p'}{m_e} - 1 \approx \frac{1 + 4\alpha}{8\alpha^2\,t_{DL}},\qquad \Delta p = m_e\left[\ln t_D + 0.809 + 2S + 2S_f\right],\quad S_f = \e^{4\alpha}E_1(4\alpha)

the remainder and SfS_f are dimensionless; Δp\Delta p and mem_e are in the same pressure units; natural logarithm

Holds when: tDL≫1/α2t_{DL} \gg 1/\alpha^2 and tDL≫1t_{DL} \gg 1: the remainder of the derivative is small; the formula for the remainder is approximate (it overstates it by a few per cent) the disturbance has not reached any other boundaries

Partially communicating faults →
Transmissibility of a section of the fault
Tf=kflf A,qr=Tfμ (p1−p2)T_f = \frac{k_f}{l_f}\,A,\qquad q_r = \frac{T_f}{\mu}\,(p_1 - p_2)

SI

Holds when: single-phase flow; the zone is the same over the area AA the jump across the zone is smaller than the difference of the average pressures in the blocks: on the way to the fault and away from it, some pressure drop is also lost in the blocks themselves

Partially communicating faults →
Fault transmissibility multiplier
F=Δx/kΔx/k+lf/kf=α Δxα Δx+LF = \frac{\Delta x/k}{\Delta x/k + l_f/k_f} = \frac{\alpha\,\Delta x}{\alpha\,\Delta x + L}

dimensionless; Δx\Delta x and LL in the same units

Holds when: neighbouring cells have the same size and permeability; for different cells, we add the half-resistances Δxi/(2ki)\Delta x_i/(2k_i) separately the width of the zone is small compared with the cell

Partially communicating faults →

Symbols

SymbolMeaning and context
Δp\Delta p
δp\delta p
Δx\Delta x
Δp′\Delta p'
Δp∞′\Delta p'_{\infty}
Δp′(t)\Delta p'(t)
Δp(t)\Delta p(t)
Δp2′\Delta p_2'
aa
AA
  • area of the section: its length along the fault (the strike) times the thickness of the formation that is in contact across the fault, m²Partially communicating faults
dd
djd_j
dj, βjd_j,\ \beta_j
e\e
E1E_1
erfcx\mathrm{erfcx}
  • scaled complementary error function, ex2erfc(x)\e^{x^2}\mathrm{erfc}(x); for large xx it is close to 1/(xπ)1/(x\sqrt\pi), —Partially communicating faults
Δtex\Delta t_{ex}
FF
Fx, FyF_x,\ F_y
kk
kf, lfk_f,\ l_f
k, ϕ, μ, ctk,\ \phi,\ \mu,\ c_t
k, ϕ, μ, ct, rwk,\ \phi,\ \mu,\ c_t,\ r_w
LL
L1L_1
L, tL,\ t
L, LDL,\ L_D
L, rwL,\ r_w
mem_e
me, ηm_e,\ \eta
me, η, Lm_e,\ \eta,\ L
mLm_L
δpmax⁡\delta p_{\max}
Δtmax⁡\Delta t_{\max}
p1−p2p_1 - p_2
p1, p2p_1,\ p_2
pDp_D
pD∗(tD)p_D^{*}(t_D)
pip_i
pwDp_{wD}
pwsp_{ws}
qrq_r
q, B, h, μ, ϕ, ctq,\ B,\ h,\ \mu,\ \phi,\ c_t
q, B, μ, k, hq,\ B,\ \mu,\ k,\ h
q, B, μ, k, h, ϕ, ctq,\ B,\ \mu,\ k,\ h,\ \phi,\ c_t
q, B, tpq,\ B,\ t_p
qB tpqB\,t_p
rr
R∞R_\infty
r1, r2r_1,\ r_2
rwr_w
SS
SchS_{ch}
SfS_f
S, SchS,\ S_{ch}
tt
tDt_D
tD, St_D,\ S
tD, tDLt_D,\ t_{DL}
  • dimensionless time kt/(ϕμctrw2)kt/(\phi\mu c_tr_w^2), and the same time based on L2L^2: tDL=kt/(ϕμctL2)=ηt/L2t_{DL} = kt/(\phi\mu c_tL^2) = \eta t/L^2, —Partially communicating faults
tD, tDL, St_D,\ t_{DL},\ S
tDLt_{DL}
tDWt_{DW}
TfT_f
tp, Δtt_p,\ \Delta t
txt_x
uu
WW
W, L3W,\ L_3
Δtx\Delta t_x
α\alpha
α, L\alpha,\ L
γ\gamma
δ\delta
η\eta
η, γ\eta,\ \gamma
θ\theta
μ\mu
ξ\xi
ξi, ζl\xi_i,\ \zeta_l
τ\tau
ϕ, ct, h\phi,\ c_t,\ h