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

Equations: Bounded drainage area

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

The pseudo-steady-state stretch of a drawdown test
pwf(t)=pint−m∗t,m∗=qBctVp=qBct ϕhAp_{wf}(t) = p_{int} - m^*t,\qquad m^* = \frac{qB}{c_tV_p} = \frac{qB}{c_t\,\phi hA}

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

Holds when: a closed area with no inflow through the boundaries a constant rate since start-up a single-phase slightly compressible fluid, and ctc_t is constant pseudo-steady-state flow has started (we check it on the derivative)

Reservoir limit test →
Intercept of the depletion line on the pressure axis
pint=pi−qBμ4πkh(ln⁡4AγCArw2+2S)p_{int} = p_i - \frac{qB\mu}{4\pi kh}\left(\ln\frac{4A}{\gamma C_A r_w^2} + 2S\right)

SI; ln⁡\ln is the natural logarithm

Holds when: the same as for the depletion line; a homogeneous formation, and the well penetrates the whole thickness.

Reservoir limit test →
Shape factor from two lines
CA=4π mem∗ t1 exp⁡pint−p1meC_A = \frac{4\pi\,m_e}{m^*\,t_1}\,\exp\frac{p_{int} - p_1}{m_e}

any consistent units: m∗m^* is in units of pressure per the unit of time in which t1t_1 is given

Holds when: both lines are drawn from one record with a constant rate a homogeneous formation; radial flow and pseudo-steady-state flow are really visible in the record

Reservoir limit test →
How the slope in ordinary coordinates approaches its limit
∣ dpwf/ dt∣m∗−1≈2cos⁡2πxwL  exp⁡(−π2η tL2)\frac{|\dd p_{wf}/\dd t|}{m^*} - 1 \approx 2\cos^2\frac{\pi x_w}{L}\;\exp\left(-\frac{\pi^2\eta\,t}{L^2}\right)

SI; in metric units η=3.553⋅10−3 k/(ϕμct)\eta = 3.553\cdot10^{-3}\,k/(\phi\mu c_t) m²/h for tt in hours

Holds when: a closed rectangular area, a homogeneous formation a late transition period: the faster “modes” have already decayed the well is not at the middle of the length: for a well at the middle of the length, on the long axis, this mode does not appear, and the equalisation goes four times faster

Reservoir limit test →
Estimate of the time to reach pseudo-steady-state flow
tpss≈L2π2η ln⁡(200cos⁡2πxwL)t_{pss} \approx \frac{L^2}{\pi^2\eta}\,\ln\left(200\cos^2\frac{\pi x_w}{L}\right)

metric: LL and xwx_w in m, η\eta in m²/h, tpsst_{pss} in hours

Holds when: a closed rectangular area; the well is clearly shifted from the middle of the length (cos⁡2(πxw/L)\cos^2(\pi x_w/L) is not small) for a well at the middle of the length, we use the table of tDAt_{DA} of the previous section.

Reservoir limit test →
Lower estimate of the connected pore volume
Vp≳πri2hϕ,ri=0.119k tϕμctV_p \gtrsim \pi r_i^2 h\phi,\qquad r_i = 0.119\sqrt{\frac{k\,t}{\phi\mu c_t}}

metric oilfield units: kk in mD, tt in h, μ\mu in mPa·s, ctc_t in 1/MPa; rir_i in m

Holds when: there are no signs of boundaries in the whole record: the derivative is on the plateau an estimate of the order of magnitude: different definitions of the radius of investigation differ by tens of per cent

Reservoir limit test →
Rate-normalised pressure and material balance time
pi−pwf(t)q(t)=b+BctVp tmb,tmb=Q(t)q(t)\frac{p_i - p_{wf}(t)}{q(t)} = b + \frac{B}{c_tV_p}\,t_{mb},\qquad t_{mb} = \frac{Q(t)}{q(t)}

metric: for qq in m³/d, QQ in m³ and ctc_t in 1/MPa, the time tmbt_{mb} is in days, and the slope is in MPa/(m³/d) per day

Holds when: a closed area, pseudo-steady-state flow the rate changes slowly compared with the time of pressure equalisation; we do not use the points right after a sharp change of the rate (tens of hours) a single-phase fluid, and BB and ctc_t are constant

Reservoir limit test →
A build-up in a closed area as the sum 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: a single-phase slightly compressible fluid, a homogeneous formation, a closed area a constant rate before shut-in (for a changing rate, the equivalent time Np/qN_p/q) we ignore wellbore storage (the well acting as a small tank): the well is shut in at the bottom, or the early points are not used

Build-up in a bounded reservoir and average reservoir pressure →
The pressure to which a build-up in a closed area levels off
pws(Δt→∞)=pi−2πtpDA qBμ2πkh=pi−qB tpctVp=pˉp_{ws}(\Delta t \to \infty) = p_i - 2\pi t_{pDA}\,\frac{qB\mu}{2\pi kh} = p_i - \frac{qB\,t_p}{c_t V_p} = \bar p

SI; in metric units the pressure drop equals qBtp/(24 ctVp)qBt_p/(24\,c_tV_p) for qq in m³/d, tpt_p in hours and ctc_t in 1/MPa

Holds when: ΔtDA=kΔt/(ϕμctA)\Delta t_{DA} = k\Delta t/(\phi\mu c_tA) is larger than the time to reach pseudo-steady-state flow for the given shape of the area the area is closed, the neighbouring wells do not work, and the fluid is single-phase

Build-up in a bounded reservoir and average reservoir pressure →
The dimensionless MBH correction
pDMBH(tpDA)=4πkhqBμ (p∗−pˉ)=2.303 (p∗−pˉ)mp_{DMBH}(t_{pDA}) = \frac{4\pi kh}{qB\mu}\,(p^* - \bar p) = \frac{2.303\,(p^* - \bar p)}{m}

SI; the second form is right in any units, if p∗p^*, pˉ\bar p and mm are in the same pressure units

Holds when: the build-up has a radial-flow line (we confirm it with the pressure derivative) a homogeneous formation, a single-phase fluid, a closed area, and one well in it (or the area is a part of a block, separated from the neighbouring wells by dividing lines of flow, the lines between the drainage areas of neighbouring wells) the rate was constant before shut-in; on the stretch of the line Δt≪tp\Delta t \ll t_p

Build-up in a bounded reservoir and average reservoir pressure →
The MBH function through the solution for a producing well
pDMBH(tpDA)=4πtpDA−2[pD(tpD)−12(ln⁡tpD+0.809)]p_{DMBH}(t_{pDA}) = 4\pi t_{pDA} - 2\Big[p_D(t_{pD}) - \tfrac12\big(\ln t_{pD} + 0.809\big)\Big]

dimensionless; ln⁡\ln is the natural logarithm

Holds when: the same as for the MBH function; pDp_D is calculated for a specific shape of the area and position of the well.

Build-up in a bounded reservoir and average reservoir pressure →
Average pressure after production in pseudo-steady-state flow
pˉ=p∗−m lg⁡(CA tpDA)\bar p = p^* - m\,\lg\left(C_A\,t_{pDA}\right)

any consistent units: p∗p^*, pˉ\bar p and mm in the same pressure units

Holds when: before shut-in the well had reached pseudo-steady-state flow: tpDAt_{pDA} is larger than the time from which the MBH function coincides with ln⁡(CAtpDA)\ln(C_At_{pDA}) (the table above) p∗p^* and tpDAt_{pDA} are calculated with the same producing time

Build-up in a bounded reservoir and average reservoir pressure →
Dietz time
Δtpˉ=ϕμctAk CA\Delta t_{\bar p} = \frac{\phi\mu c_t A}{k\,C_A}

SI

Holds when: before shut-in the well worked in pseudo-steady-state flow, and the producing time is much larger than Δtpˉ\Delta t_{\bar p} the MDH line is drawn through the points of radial flow

Build-up in a bounded reservoir and average reservoir pressure →
The local slope of a Horner plot after production in pseudo-steady-state flow
local slopem=(1+Δttp′)(1−4π ΔtDA)\frac{\text{local slope}}{m} = \left(1 + \frac{\Delta t}{t_p'}\right)\left(1 - 4\pi\,\Delta t_{DA}\right)

dimensionless; Δt\Delta t and tp′t_p' in the same units

Holds when: before shut-in the well worked in pseudo-steady-state flow after shut-in there is still radial flow around the well: the disturbance from the shut-in has not reached the boundaries

Build-up in a bounded reservoir and average reservoir pressure →
A late build-up in a closed area
pˉ−pws(Δt)≈a e−λ Δt,λ≈π2ηL2\bar p - p_{ws}(\Delta t) \approx a\,\e^{-\lambda\,\Delta t},\qquad \lambda \approx \frac{\pi^2\eta}{L^2}

SI; in metric units η=3.553⋅10−3 k/(ϕμct)\eta = 3.553\cdot10^{-3}\,k/(\phi\mu c_t) m²/h, and λ\lambda comes out in 1/h

Holds when: a late period: the other, faster “modes” have already decayed (for a rectangle this is ΔtDA\Delta t_{DA} of about 0.2–0.5 and more) the formula for λ\lambda is for a rectangular area and a well that is not at the middle of its length. At the middle, this mode does not appear, and the decay is faster. For a well at the centre of a circle, λ=14.7 η/re2\lambda = 14.7\,\eta/r_e^2 one well in a closed area, and the neighbours do not work

Build-up in a bounded reservoir and average reservoir pressure →
Average pressure of a closed area and cumulative production
pˉ=pi−BctVp Np\bar p = p_i - \frac{B}{c_t V_p}\,N_p

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

Holds when: a closed area with no water inflow; single-phase oil above the bubble point pressure (the pressure below which gas comes out of the oil); BB and ctc_t are constant there is one well in the area, or NpN_p and VpV_p refer to the whole closed block

Build-up in a bounded reservoir and average reservoir pressure →
A build-up after several flow periods in the Horner form
pws(Δt)=p∗−m XBU(Δt),XBU=∑j=1nqj−qj−1qn lg⁡θj+ΔtΔtp_{ws}(\Delta t) = p^* - m\,X_{BU}(\Delta t),\qquad X_{BU} = \sum_{j=1}^{n}\frac{q_j - q_{j-1}}{q_n}\,\lg\frac{\theta_j + \Delta t}{\Delta t}

metric oilfield units: pressure and mm in MPa, qq in m³/d, μ\mu in mPa·s, kk in mD, hh in m; time in any equal units

Holds when: radial flow in an infinite formation over the whole history and the build-up: the disturbance has not reached any boundaries a linear problem (the effects of separate events add up): a single-phase slightly compressible fluid, a constant skin wellbore storage (the well acting as a small tank) is over

Equivalent producing time and synthetic flowing time →
The superposition derivative
ΔpX′=− d Δp d (XBUln⁡10)=Δp′∑j=1nqj−qj−1qn θjθj+Δt\Delta p'_X = -\frac{\dd\,\Delta p}{\dd\,(X_{BU}\ln 10)} = \frac{\Delta p'}{\displaystyle\sum_{j=1}^{n}\frac{q_j - q_{j-1}}{q_n}\,\frac{\theta_j}{\theta_j + \Delta t}}

any consistent system; derivatives in units of pressure

Holds when: the conditions of the previous equation; the minus sign is there because XBUX_{BU} decreases as Δt\Delta t grows.

Equivalent producing time and synthetic flowing time →
The Horner equivalent producing time
tp=Npqnt_p = \frac{N_p}{q_n}

NpN_p in m³, qnq_n in m³/d, tpt_p in days; for the Horner ratio in hours we multiply by 24

Holds when: the last rate was held long enough compared with the duration of the build-up (an estimate below) the history went on in an infinitely acting formation; after a long production in a bounded area we use the synthetic time (the sections below)

Equivalent producing time and synthetic flowing time →
The error of the slope with one producing time
R(Δt)=tp+Δttp ∑j=1nqj−qj−1qn θjθj+ΔtR(\Delta t) = \frac{t_p + \Delta t}{t_p}\,\sum_{j=1}^{n}\frac{q_j - q_{j-1}}{q_n}\,\frac{\theta_j}{\theta_j + \Delta t}

dimensionless; times in the same units

Holds when: radial flow in an infinite formation; wellbore storage is over.

Equivalent producing time and synthetic flowing time →
Synthetic flowing time
tsDA=k tsϕμctA=1CA exp⁡pe−pˉmet_{sDA} = \frac{k\,t_s}{\phi\mu c_tA} = \frac{1}{C_A}\,\exp\frac{p_e - \bar p}{m_e}

SI; in metric oilfield units tsDA=3.553⋅10−3 kts/(ϕμctA)t_{sDA} = 3.553\cdot10^{-3}\,kt_s/(\phi\mu c_tA): mD, h, mPa·s, 1/MPa, m²

Holds when: before shut-in the well had reached pseudo-steady-state flow we analyse the build-up by the radial-flow line, while the shut-in has not yet been felt at the boundaries the ratio (pe−pˉ)/me(p_e - \bar p)/m_e depends only on the shape of the area, the position of the well and the chosen point of the boundary. We calculate it with a model of the area (the table below)

Equivalent producing time and synthetic flowing time →
A build-up as the sum of the trend and the response to the shut-in
pws(Δt)=ptr(Δt)+qBμ2πkh[pD(ΔtD)+S]p_{ws}(\Delta t) = p_{tr}(\Delta t) + \frac{qB\mu}{2\pi kh}\Big[p_D(\Delta t_D) + S\Big]

SI

Holds when: the problem is linear: a single-phase slightly compressible fluid, and the properties of the formation do not depend on the pressure the rates of the neighbours do not change during the build-up, and if they changed earlier, their history is in the trend we leave out wellbore storage (the well acting as a small tank): the well is shut in at the bottom, or we do not use the early points

Testing in a developed reservoir: Slider's method →
The pressure decline rate of a block in pseudo-steady-state flow
m∗=− dp dt=∑jqjBctVpm^* = -\frac{\dd p}{\dd t} = \frac{\sum_j q_jB}{c_tV_p}

SI

Holds when: the block is closed, the fluid is single-phase, and the rates have been constant so long that the transient effects in the whole block have died out there is no water influx (water flowing in from outside) and no gas cap (free gas above the oil)

Testing in a developed reservoir: Slider's method →
Slider's correction: the reconstructed drawdown test
ΔpDD(Δt)=pws(Δt)−pwf(Δt=0)+m∗Δt\Delta p_{DD}(\Delta t) = p_{ws}(\Delta t) - p_{wf}(\Delta t = 0) + m^*\Delta t

any consistent units: pressure in MPa, m∗m^* in MPa/h and Δt\Delta t in hours

Holds when: the trend is linear: the block is in pseudo-steady-state flow, and the rates of the neighbours did not change either before the build-up or during it m∗m^* is found from the record before shut-in or from the material balance of the block (the section below)

Testing in a developed reservoir: Slider's method →
The correction for the trend as a fraction of the plateau level
m∗Δtme=4π ΔtDAw,ΔtDAw=k ΔtϕμctAw,Aw=A q∑jqj\frac{m^*\Delta t}{m_e} = 4\pi\,\Delta t_{DA_w},\qquad \Delta t_{DA_w} = \frac{k\,\Delta t}{\phi\mu c_tA_w},\qquad A_w = A\,\frac{q}{\sum_j q_j}

dimensionless; in metric oilfield units ΔtDAw=3.553⋅10−3 k Δt/(ϕμctAw)\Delta t_{DA_w} = 3.553\cdot10^{-3}\,k\,\Delta t/(\phi\mu c_tA_w) for kk in mD, Δt\Delta t in hours, μ\mu in mPa·s, ctc_t in 1/MPa, AwA_w in m²

Holds when: the block is in pseudo-steady-state flow, and all the wells work at constant rates in one homogeneous formation.

Testing in a developed reservoir: Slider's method →
The derivative of a raw build-up in a developed accumulation
ΔpBU′=ΔpDD′−m∗Δt\Delta p_{BU}' = \Delta p_{DD}' - m^*\Delta t

any consistent units

Holds when: a linear trend.

Testing in a developed reservoir: Slider's method →

Symbols

SymbolMeaning and context
pˉ\bar p
Δtpˉ\Delta t_{\bar p}
Δp\Delta p
Δt\Delta t
Δp′\Delta p'
ΔpX′\Delta p'_X
12(ln⁡tpD+0.809)\tfrac12(\ln t_{pD} + 0.809)
aa
AA
AwA_w
A, k, ϕ, μ, ctA,\ k,\ \phi,\ \mu,\ c_t
A, q, ∑jqjA,\ q,\ \sum_j q_j
A, CAA,\ C_A
A, rw, SA,\ r_w,\ S
bb
  • drawdown per unit rate: 0.933 Bμkh(ln⁡4AγCArw2+2S)0.933\,\frac{B\mu}{kh}\left(\ln\frac{4A}{\gamma C_Ar_w^2} + 2S\right) for kk in mD, μ\mu in mPa·s, hh in m; 0.933=1/(4π⋅0.0853)0.933 = 1/(4\pi\cdot0.0853) is the constant of the slope per unit of natural logarithm (where the constants come from), MPa/(m³/d)Reservoir limit test
BB
B, ctB,\ c_t
B, ct, VpB,\ c_t,\ V_p
  • formation volume factor, total compressibility, pore volume, m³/m³, 1/MPa, m³Reservoir limit test
ΔpBU′\Delta p_{BU}'
CAC_A
ctc_t
ct, Vpc_t,\ V_p
ΔtDA\Delta t_{DA}
ΔtDAw\Delta t_{DA_w}
ΔpDD′\Delta p_{DD}'
ΔpDD(Δt)\Delta p_{DD}(\Delta t)
k, hk,\ h
k, ϕ, μ, ctk,\ \phi,\ \mu,\ c_t
k, ϕ, μ, ct, hk,\ \phi,\ \mu,\ c_t,\ h
  • permeability, porosity, viscosity, total compressibility, thickness, mD, —, mPa·s, 1/MPa, mReservoir limit test
LL
L, xw, ηL,\ x_w,\ \eta
mm
m∗m^*
m∗Δtm^*\Delta t
m∗, pintm^*,\ p_{int}
mem_e
NpN_p
p∗p^*
p∗, mp^*,\ m
p1p_1
pD(ΔtD)p_D(\Delta t_D)
pD(tpD)p_D(t_{pD})
pD(t)p_D(t)
pe−pˉp_e - \bar p
pip_i
pintp_{int}
ptr(Δt)p_{tr}(\Delta t)
  • the trend: the bottomhole pressure that the well would have had at the same moment if it had kept producing at the rate qq; at the moment of shut-in it equals pwf(Δt=0)p_{wf}(\Delta t = 0), PaTesting in a developed reservoir: Slider's method
pwf(Δt=0)p_{wf}(\Delta t = 0)
pws(Δt)p_{ws}(\Delta t)
qjq_j
qnq_n
q, Bq,\ B
  • rate at standard conditions and formation volume factor; qBqB is the rate at reservoir conditions, m³/s, m³/m³Reservoir limit test
q, B, μq,\ B,\ \mu
q, B, μ, k, hq,\ B,\ \mu,\ k,\ h
q(t)q(t)
Q(t)Q(t)
rir_i
R(Δt)R(\Delta t)
SS
tt
tmbt_{mb}
tpt_p
tp′t_p'
tpDt_{pD}
tpDAt_{pDA}
tpsst_{pss}
  • time from which the slope in ordinary coordinates differs from the limit by less than 1%, hReservoir limit test
tst_s
tsDAt_{sDA}
VpV_p
Vp=ϕhAV_p = \phi hA
XBUX_{BU}
xwx_w
γ\gamma
η=k/(ϕμct)\eta = k/(\phi\mu c_t)
θj\theta_j
θj, qj\theta_j,\ q_j
λ\lambda