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

Equations: Transient testing: fundamentals

The equations of the “Transient testing: fundamentals” lessons in “Well testing and pressure transient analysis”, with their conditions and units, and what each symbol means in them.

Rate of pressure decline in a closed area
 dp dt=−qBctVp\frac{\dd p}{\dd t} = -\frac{qB}{c_t V_p}

any consistent units (for example, m³/d, MPa/d, 1/MPa, m³)

Holds when: The area is closed (nothing flows in across the boundaries), and the rate is constant. Pseudo-steady state has already been reached: the disturbance has reached all the boundaries. The fluid is single-phase and slightly compressible, with constant compressibility.

Well testing: purposes and types of test →
Radius of investigation
ri=4ktϕμctr_i = \sqrt{\frac{4kt}{\phi\mu c_t}}

SI

Holds when: The reservoir is homogeneous, and the disturbance has not yet reached the boundaries. The criterion kt/(ϕμctri2)=1/4kt/(\phi\mu c_t r_i^2) = 1/4 is a convention. If we choose a different level for what counts as a “noticeable” disturbance, the constant changes by tens of per cent. In practice, whether we can detect a disturbance also depends on the resolution of the gauge (the smallest change it can show).

Well testing: purposes and types of test →
Pressure drawdown in an infinite reservoir (logarithmic approximation)
pi−pwf=qBμ4πkh[ln⁡ktϕμctrw2+0.809+2S]p_i - p_{wf} = \frac{qB\mu}{4\pi kh}\left[\ln\frac{kt}{\phi\mu c_t r_w^2} + 0.809 + 2S\right]

SI; the logarithm is natural

Holds when: The rate is constant from start-up, and wellbore storage is ignored. The disturbance has not reached the boundaries, and the dimensionless time kt/(ϕμctrw2)kt/(\phi\mu c_t r_w^2) is greater than about 25. The reservoir is homogeneous, and the fluid is single-phase and slightly compressible.

Well testing: purposes and types of test →
Bottomhole pressure in pseudo-steady-state flow
pi−pwf=qB tctVp+qJ,J=2πkhBμ(ln⁡rerw−34+S)p_i - p_{wf} = \frac{qB\,t}{c_t V_p} + \frac{q}{J},\qquad J = \frac{2\pi kh}{B\mu\left(\ln\dfrac{r_e}{r_w} - \dfrac34 + S\right)}

SI

Holds when: The area is circular and closed, with the well at the centre. The time is longer than the time needed to reach pseudo-steady state.

Well testing: purposes and types of test →
Conversion of a dimensionless solution to pressure drop and time
Δp=qBμ2πkh pD,t=ϕμctrw2k tD\Delta p = \frac{qB\mu}{2\pi kh}\,p_D,\qquad t = \frac{\phi\mu c_t r_w^2}{k}\,t_D

SI

Holds when: The conditions of the diffusivity equation hold: the fluid is single-phase and slightly compressible, and the properties of the reservoir and the fluid are constant (Diffusivity equation). Δp\Delta p is positive for production. For an injection well the sign is opposite.

Dimensionless well solutions: skin, radius of investigation and units →
Flow rate through a circle in the line-source solution
q(r,t)q=exp⁡ ⁣(−r24ηt)=exp⁡ ⁣(−rD24tD)\frac{q(r, t)}{q} = \exp\!\left(-\frac{r^2}{4\eta t}\right) = \exp\!\left(-\frac{r_D^2}{4t_D}\right)

no units (a ratio of rates)

Holds when: This is the line-source solution: the reservoir is infinite and homogeneous, and the rate is constant from start-up (the derivation is in the advanced section).

Dimensionless well solutions: skin, radius of investigation and units →
Two ways of accounting for skin in the solution for a well
pwD(tD)=pD(tD)+SorpwD(tD)=pD ⁣(tD e2S)p_{wD}(t_D) = p_D(t_D) + S\qquad\text{or}\qquad p_{wD}(t_D) = p_D\!\left(t_D\,\e^{2S}\right)

no units

Holds when: The first form is the skin of an infinitely thin zone that follows the rate instantly. The second form is a well without skin, of radius rw′=rwe−Sr_w' = r_w\e^{-S}, in the same reservoir. The pressure is still referred to the kk and hh of the reservoir. Wellbore storage is not included. It can be added to either form (Wellbore storage).

Dimensionless well solutions: skin, radius of investigation and units →
Rate of pressure decline at a point in the reservoir
∂ Δp∂t=qBμ4πkh⋅1t exp⁡ ⁣(−r24ηt)\frac{\partial\,\Delta p}{\partial t} = \frac{qB\mu}{4\pi kh}\cdot\frac1t\,\exp\!\left(-\frac{r^2}{4\eta t}\right)

SI

Holds when: This is the line-source solution: the reservoir is infinite, and the rate is constant.

Dimensionless well solutions: skin, radius of investigation and units →
Share of the withdrawn volume supplied by the reservoir inside radius r
V(r)V=1−e−X+X E1(X),X=r24ηt=(rri)2\frac{V(r)}{V} = 1 - \e^{-X} + X\,E_1(X),\qquad X = \frac{r^2}{4\eta t} = \left(\frac{r}{r_i}\right)^2

no units

Holds when: This is the line-source solution: the reservoir is infinite, and the rate is constant from start-up (the derivation is in the advanced section).

Dimensionless well solutions: skin, radius of investigation and units →
The radius of drainage in transient flow
rd=4ηtγ=0.749 ri,pi−pwf=qBμ2πkh(ln⁡rdrw+S)r_d = \sqrt{\frac{4\eta t}{\gamma}} = 0.749\,r_i,\qquad p_i - p_{wf} = \frac{qB\mu}{2\pi kh}\left(\ln\frac{r_d}{r_w} + S\right)

SI; the logarithm is natural

Holds when: The rate is constant, and the reservoir is infinite and homogeneous. tD>25t_D > 25, so the logarithmic approximation is valid at the wall of the well.

Dimensionless well solutions: skin, radius of investigation and units →
Transient productivity index
Jt=qpi−pwf(t)=2πkhBμ[ln⁡(rd/rw)+S]J_t = \frac{q}{p_i - p_{wf}(t)} = \frac{2\pi kh}{B\mu\left[\ln(r_d/r_w) + S\right]}

SI

Holds when: The flow is transient: the disturbance has not reached the boundaries. We measure the drawdown from the initial pressure pip_i, because the average pressure in an infinite reservoir is not defined. If the drawdown is held constant, the rate falls by the same formula, with a small correction (the rate at constant pressure).

Dimensionless well solutions: skin, radius of investigation and units →
Bottomhole pressure at a constant rate in an infinite reservoir
pwf=pi−qBμ4πkh[ln⁡ktϕμctrw2+0.809+2S]p_{wf} = p_i - \frac{qB\mu}{4\pi kh}\left[\ln\frac{kt}{\phi\mu c_t r_w^2} + 0.809 + 2S\right]

SI; natural logarithm

Holds when: The rate is constant from the moment of start-up; before start-up the pressure throughout the reservoir equals pip_i. The reservoir behaves as an infinite one: the disturbance has not reached the boundaries. Wellbore storage is over. The condition kt/(ϕμctrw2)>25kt/(\phi\mu c_tr_w^2) > 25 is met within seconds. The reservoir is homogeneous, there is a single slightly compressible fluid, and the skin does not change with time.

Pressure drawdown analysis →
Flow capacity from the slope of the semilog straight line
kh=qBμ ln⁡104πmkh = \frac{qB\mu\,\ln 10}{4\pi m}

SI

Holds when: The slope is measured on the section of infinite-acting radial flow. The rate is constant. If the thickness of the producing part of the reservoir is not known, we report the result as khkh.

Pressure drawdown analysis →
Skin factor from a drawdown test
S=1.151[pi−p1hm−lg⁡kϕ μ ct rw2+2.098]S = 1.151\left[\frac{p_i - p_{1\text{h}}}{m} - \lg\frac{k}{\phi\,\mu\,c_t\,r_w^2} + 2.098\right]

metric oilfield units: pressure in MPa, kk in mD, μ\mu in mPa·s, ctc_t in 1/MPa, rwr_w in m, time in hours

Holds when: We read p1hp_{1\text{h}} on the line and do not take it from the measurements, because wellbore storage may not be over an hour after start-up. The well stood idle long enough before start-up for the pressure in the reservoir to equalise. Otherwise we need superposition with the previous history (adding up the effects of earlier rate changes) in place of pip_i. mm and kk are found from the same straight section.

Pressure drawdown analysis →
Flow capacity and skin from a plot of pressure against the natural logarithm of time
kh=0.933 qBμme,S=12[pi−p1hme−ln⁡kϕμctrw2+4.831]kh = 0.933\,\frac{qB\mu}{m_e},\qquad S = \frac12\left[\frac{p_i - p_{1\text{h}}}{m_e} - \ln\frac{k}{\phi\mu c_t r_w^2} + 4.831\right]

metric oilfield units; time in hours

Holds when: The same as for the formulas with the base-10 logarithm.

Pressure drawdown analysis →
Pressure build-up as the superposition of two rates
pi−pws(Δt)=qBμ2πkh[pD(tp+Δt)−pD(Δt)]p_i - p_{ws}(\Delta t) = \frac{qB\mu}{2\pi kh}\Big[p_D(t_p + \Delta t) - p_D(\Delta t)\Big]

SI (any consistent system of units)

Holds when: The problem is linear: one slightly compressible fluid, and the properties of the reservoir and the fluid do not depend on the pressure. The rate was constant before shut-in. The inflow from the reservoir stops at once at shut-in (we ignore wellbore storage).

Pressure build-up analysis →
The Horner equation
pws=pi−qBμ4πkh ln⁡tp+ΔtΔtp_{ws} = p_i - \frac{qB\mu}{4\pi kh}\,\ln\frac{t_p + \Delta t}{\Delta t}

SI (any consistent system of units)

Holds when: The reservoir behaves as an infinite one both while the well produces and during the shut-in: the disturbance has not reached the boundaries. The rate before shut-in is constant. The reservoir is homogeneous, and there is one slightly compressible fluid. Wellbore storage is over: the inflow from the reservoir into the wellbore has stopped after shut-in.

Pressure build-up analysis →
Flow capacity from the slope of the Horner line
kh=qBμ ln⁡104πmkh = \frac{qB\mu\,\ln 10}{4\pi m}

SI

Holds when: The slope is measured on the section of radial flow. hh is the thickness of the producing part of the reservoir. If it is not known, the result stays khkh.

Pressure build-up analysis →
Skin factor from a build-up
S=1.151[p1h−pwf(Δt=0)m−lg⁡kϕ μ ct rw2+2.098]S = 1.151\left[\frac{p_{1\text{h}} - p_{wf}(\Delta t = 0)}{m} - \lg\frac{k}{\phi\,\mu\,c_t\,r_w^2} + 2.098\right]

metric oilfield units: pressure in MPa, kk in mD, μ\mu in mPa·s, ctc_t in 1/MPa, rwr_w in m, time in hours

Holds when: The producing time tpt_p is much longer than 1 h. In the general case we add the term lg⁡[(tp+1)/tp]\lg[(t_p + 1)/t_p] to the brackets. Both the production and the shut-in take place as radial flow in an infinite reservoir. We read p1hp_{1\text{h}} on the line, even if the measurement itself at 1 h is still distorted by wellbore storage.

Pressure build-up analysis →
The MDH line
pws=p1h+m lg⁡Δtp_{ws} = p_{1\text{h}} + m\,\lg\Delta t

metric oilfield units (pressure in MPa, time in hours)

Holds when: Δt≪tp\Delta t \ll t_p over the whole section we use. If the shut-in time is no more than about a tenth of the producing time, the slope is distorted by no more than a few per cent. The other conditions are the same as for the Horner method.

Pressure build-up analysis →
Radius of the “probe” in a build-up
rp=4k Δteγ ϕμct,Δte=tp Δttp+Δtr_p = \sqrt{\frac{4k\,\Delta t_e}{\gamma\,\phi\mu c_t}},\qquad \Delta t_e = \frac{t_p\,\Delta t}{t_p + \Delta t}

SI

Holds when: The reservoir is infinite. At the point rpr_p the logarithmic approximation holds: the disturbance passed it long ago. Wellbore storage is over. The skin changes the pressure only in the well itself, and it does not enter the formula.

Pressure build-up analysis →

Symbols

SymbolMeaning and context
Δp\Delta p
Δt\Delta t
(tp+Δt)/Δt(t_p + \Delta t)/\Delta t
 dp/ dt\dd p/\dd t
BB
ctc_t
Δte\Delta t_e
E1(X)E_1(X)
e2S\e^{2S}
JJ
JtJ_t
kk
k, hk,\ h
k, h, B, μ, rw, Sk,\ h,\ B,\ \mu,\ r_w,\ S
k, h, μk,\ h,\ \mu
k, h, ϕ, ct, rwk,\ h,\ \phi,\ c_t,\ r_w
k, ϕ, μ, ctk,\ \phi,\ \mu,\ c_t
khkh
mm
mem_e
  • pressure drop on the line when ln⁡t\ln t increases by one, that is, when time grows by a factor of 2.718, MPaPressure drawdown analysis
m, p1hm,\ p_{1\text{h}}
  • in the oilfield forms: the pressure drop per decimal cycle of time and the pressure on the line at tt = 1 h, MPa or psiPressure drawdown analysis
p1hp_{1\text{h}}
pD, tDp_D,\ t_D
pD(tD)p_D(t_D)
pD(t)p_D(t)
pip_i
pi, q, B, μp_i,\ q,\ B,\ \mu
pi, pwfp_i,\ p_{wf}
pwDp_{wD}
pwfp_{wf}
pwf(Δt=0)p_{wf}(\Delta t = 0)
pwf(t)p_{wf}(t)
pwsp_{ws}
pws(Δt)p_{ws}(\Delta t)
qq
q, Bq,\ B
q, B, μq,\ B,\ \mu
q, B, μ, k, hq,\ B,\ \mu,\ k,\ h
q, B, μ, k, h, Sq,\ B,\ \mu,\ k,\ h,\ S
q(r,t)q(r, t)
rdr_d
rD, tDr_D,\ t_D
rer_e
rir_i
rpr_p
  • distance from the well at which the pressure at the moment of shut-in was equal to the reading pws(Δt)p_{ws}(\Delta t), mPressure build-up analysis
rwr_w
SS
tt
tp, Δtt_p,\ \Delta t
VV
VpV_p
V(r)V(r)
γ\gamma
η\eta
μ\mu
ϕ\phi
ϕ, ct, rw\phi,\ c_t,\ r_w