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

Equations: Inflow to a well

The equations of the “Inflow to a well” lessons in “Well testing and pressure transient analysis”, with their conditions and units, and what each symbol means in them.

Inflow line (single-phase inflow)
pwf=pˉ−qJp_{wf} = \bar p - \frac{q}{J}

any consistent units, for example m³/d and MPa

Holds when: Single-phase inflow (pwfp_{wf} above the bubble point pressure) and stabilised flow. The productivity index does not depend on the flow rate: there are no inertial losses.

Well inflow: the well performance diagram →
Lift curve for liquid without free gas
pwf=pwh+ρgH+Kq2p_{wf} = p_{wh} + \rho g H + K q^2

pressures in MPa; the product ρgH\rho g H comes out in pascals and is divided by 10610^6

Holds when: Vertical tubing down to the perforated interval; single-phase liquid without free gas, constant density. The wellhead pressure is held constant — by the choke or by the pressure in the gathering system (the network of pipes that collects the production). KK depends on the tubing diameter and on the density and viscosity of the liquid. The lesson Nodal analysis of the production system shows how to calculate it.

Well inflow: the well performance diagram →
Injection balance under reservoir pressure maintenance
∑jqw,jBw+qaq=qoBo\sum_{j} q_{w,j} B_w + q_{aq} = q_o B_o

any identical flow-rate units, for example m³/d

Holds when: We consider the period before water breaks through to the production wells. The pressure is above the bubble point pressure, so only oil with dissolved gas is produced. The reservoir pressure is held constant: the withdrawal is fully replaced.

Well inflow: the well performance diagram →
Underground withdrawal of oil, water and free gas
V=qoBo+qwBw+qo (Rp−Rs) BgV = q_o B_o + q_w B_w + q_o\,(R_p - R_s)\,B_g

m³/d; the producing and solution gas–oil ratios are in m³ of gas per m³ of oil at standard conditions

Holds when: The flow rates are at standard conditions; the formation volume factors and the solution gas–oil ratio are at the current reservoir pressure. Free gas is the gas in excess of the dissolved gas, (Rp−Rs) qo(R_p - R_s)\,q_o at standard conditions.

Well inflow: the well performance diagram →
Steady-state radial flow (Dupuit equation)
q=2πkh (pe−pw)μln⁡rerwq = \frac{2\pi kh\,(p_e - p_w)}{\mu\ln\dfrac{r_e}{r_w}}

SI (any consistent system of units); ln⁡\ln is the natural logarithm

Holds when: Steady state: a constant pressure pep_e is held at the outer boundary. The reservoir is homogeneous and isotropic (its permeability is the same in all directions) and has constant thickness. The well stands at the centre of a circular area and penetrates the whole thickness. The flow is single-phase. The fluid is incompressible or slightly compressible and obeys Darcy's law. The flow is horizontal. The near-wellbore zone is unchanged, so the skin factor is zero (the lesson shows below how to account for it).

Steady-state radial flow →
Pressure distribution in steady-state radial flow
p(r)=pw+qμ2πkhln⁡rrwp(r) = p_w + \frac{q\mu}{2\pi kh}\ln\frac{r}{r_w}

SI; ln⁡\ln is the natural logarithm

Holds when: The same as for the Dupuit equation.

Steady-state radial flow →
Productivity index in steady-state inflow
J=qspe−pw=2πkhBμln⁡rerwJ = \frac{q_s}{p_e - p_w} = \frac{2\pi kh}{B\mu\ln\dfrac{r_e}{r_w}}

SI; ln⁡\ln is the natural logarithm

Holds when: The conditions of the Dupuit equation hold, and the skin is zero.

Steady-state radial flow →
Steady-state inflow with a skin factor
qs=2πkh (pe−pwf)Bμ(ln⁡rerw+S)q_s = \frac{2\pi kh\,(p_e - p_{wf})}{B\mu\left(\ln\dfrac{r_e}{r_w} + S\right)}

SI; in metric and oilfield units the constants are the same as in the Dupuit equation

Holds when: The altered zone is small compared with the drainage area, and the flow in it can be treated as steady-state. The skin does not depend on the flow rate (there are no inertial losses).

Steady-state radial flow →
Average pressure in steady-state radial flow
pˉ−pw=qsBμ2πkh(ln⁡rerw−12),pe−pˉ=12⋅qsBμ2πkh\bar p - p_w = \frac{q_s B\mu}{2\pi kh}\left(\ln\frac{r_e}{r_w} - \frac12\right),\qquad p_e - \bar p = \frac12\cdot\frac{q_s B\mu}{2\pi kh}

SI; ln⁡\ln is the natural logarithm

Holds when: Steady-state inflow in a circular area; porosity and thickness are constant. rw≪rer_w \ll r_e: terms of the order of (rw/re)2(r_w/r_e)^2 are left out.

Steady-state radial flow →
Steady-state radial gas inflow
qsc=πkh Tsc (pe2−pw2)psc T μˉZˉln⁡rerwq_{sc} = \frac{\pi kh\,T_{sc}\,(p_e^2 - p_w^2)}{p_{sc}\,T\,\bar\mu\bar Z\ln\dfrac{r_e}{r_w}}

SI; absolute pressures only; ln⁡\ln is the natural logarithm

Holds when: Steady-state inflow at constant temperature (isothermal). The reservoir is homogeneous, and the well penetrates its whole thickness. μZ\mu Z changes little between pwp_w and pep_e, which holds roughly for pressures up to 14 MPa. At higher pressures we use the pseudo-pressure (advanced section). Darcy's law holds. In gas wells at high flow rates, inertial losses are added (Inflow at high rates: non-Darcy flow). The skin is added to the logarithm in the same way as for oil.

Steady-state radial flow →
Quasi-steady-state pressure profile near the well
p(r,t)=pwf(t)+qsBμ2πkh(ln⁡rrw+S)p(r, t) = p_{wf}(t) + \frac{q_sB\mu}{2\pi kh}\left(\ln\frac{r}{r_w} + S\right)

SI; ln⁡\ln is the natural logarithm

Holds when: Constant flow rate after start-up; a homogeneous reservoir. rr is much smaller than the radius of investigation 4kt/(ϕμct)\sqrt{4kt/(\phi\mu c_t)} (how far into the reservoir the disturbance has travelled): the disturbance reached this point long ago.

Steady-state radial flow →
Gas inflow in terms of pseudo-pressure
m(p)=2∫p0pp′  dp′μ(p′) Z(p′),qsc=πkh Tsc [m(pe)−m(pw)]psc T ln⁡rerwm(p) = 2\int_{p_0}^{p}\frac{p'\,\dd p'}{\mu(p')\,Z(p')},\qquad q_{sc} = \frac{\pi kh\,T_{sc}\,\big[m(p_e) - m(p_w)\big]}{p_{sc}\,T\,\ln\dfrac{r_e}{r_w}}

SI; absolute pressures

Holds when: Steady-state inflow at constant temperature, obeying Darcy's law. The gas properties are known from PVT data (laboratory measurements of how pressure, volume and temperature are linked).

Steady-state radial flow →
Rate of pressure decline in a closed area
 dpˉ dt=−qBctVp=−qBct ϕhA\frac{\dd \bar p}{\dd t} = -\frac{qB}{c_t V_p} = -\frac{qB}{c_t\,\phi h A}

any consistent units: for example, m³/d, 1/MPa and m³ give MPa/d

Holds when: A closed area: no fluid enters across the boundaries. Constant flow rate; single-phase fluid, constant compressibility. For the average pressure the equality holds at any time. For the pressure at each point it holds only after the pseudo-steady-state regime has set in.

Pseudo-steady-state flow and the shape factor →
Flow rate through a circle of radius r
q(r)=q(1−r2re2)q(r) = q\left(1 - \frac{r^2}{r_e^2}\right)

any consistent units

Holds when: Pseudo-steady-state regime, well at the centre of a circular area, homogeneous reservoir. rw≪rer_w \ll r_e.

Pseudo-steady-state flow and the shape factor →
Pressure distribution in pseudo-steady-state inflow
p(r)−pwf=qBμ2πkh(ln⁡rrw−r22re2+S)p(r) - p_{wf} = \frac{qB\mu}{2\pi kh}\left(\ln\frac{r}{r_w} - \frac{r^2}{2r_e^2} + S\right)

SI; ln⁡\ln is the natural logarithm

Holds when: Pseudo-steady-state regime, circular area with the well at the centre. rw≪rer_w \ll r_e: terms of the order of (rw/re)2(r_w/r_e)^2 are left out. The skin is an additional pressure drop in a thin zone at the wellbore (Skin factor).

Pseudo-steady-state flow and the shape factor →
Pseudo-steady-state inflow to a well at the centre of a circular area
pe−pwf=qBμ2πkh(ln⁡rerw−12+S)pˉ−pwf=qBμ2πkh(ln⁡rerw−34+S)\begin{aligned} p_e - p_{wf} &= \frac{qB\mu}{2\pi kh}\left(\ln\frac{r_e}{r_w} - \frac12 + S\right)\\ \bar p - p_{wf} &= \frac{qB\mu}{2\pi kh}\left(\ln\frac{r_e}{r_w} - \frac34 + S\right)\end{aligned}

SI (any consistent system of units); ln⁡\ln is the natural logarithm

Holds when: The same as for the pressure distribution. The fluid is single-phase and slightly compressible, Darcy's law holds, and the skin does not depend on the flow rate.

Pseudo-steady-state flow and the shape factor →
Pseudo-steady-state inflow in an area of arbitrary shape
pˉ−pwf=qBμ2πkh(12ln⁡4AγCArw2+S)\bar p - p_{wf} = \frac{qB\mu}{2\pi kh}\left(\frac12\ln\frac{4A}{\gamma C_A r_w^2} + S\right)

SI; AA and rw2r_w^2 in the same units of area; ln⁡\ln is the natural logarithm

Holds when: Pseudo-steady-state regime: the dimensionless time tDAt_{DA} is greater than the tabulated value for the given shape. A homogeneous isotropic reservoir, the well penetrates the whole thickness, and the boundary is impermeable along its whole perimeter.

Pseudo-steady-state flow and the shape factor →
Shape pseudo-skin
12ln⁡4AγCArw2=ln⁡rerw−34+SA,SA=12ln⁡4πγCA+34\frac12\ln\frac{4A}{\gamma C_A r_w^2} = \ln\frac{r_e}{r_w} - \frac34 + S_A,\qquad S_A = \frac12\ln\frac{4\pi}{\gamma C_A} + \frac34

dimensionless

Holds when: The same as for the equation with the shape factor.

Pseudo-steady-state flow and the shape factor →
Dimensionless time based on the area of the region
tDA=k tϕμctAt_{DA} = \frac{k\,t}{\phi\mu c_t A}

SI

Holds when: A homogeneous reservoir; time is counted from the start-up of the well or from the last change of flow rate.

Pseudo-steady-state flow and the shape factor →
Bottomhole pressure in the pseudo-steady-state regime in terms of the initial pressure
pi−pwf=qBμ2πkh(12ln⁡4AγCArw2+S+2πtDA)p_i - p_{wf} = \frac{qB\mu}{2\pi kh}\left(\frac12\ln\frac{4A}{\gamma C_A r_w^2} + S + 2\pi t_{DA}\right)

SI; ln⁡\ln is the natural logarithm

Holds when: Constant flow rate from the initial state; the pseudo-steady-state regime has already set in. A closed area, single-phase slightly compressible fluid.

Pseudo-steady-state flow and the shape factor →
Drainage volumes and the average pressure of the block
Vi=Vt qiqt,pˉblk=∑ipˉiVi∑iVi=∑ipˉiqi∑iqiV_i = V_t\,\frac{q_i}{q_t},\qquad \bar p_{\text{blk}} = \frac{\sum_i \bar p_i V_i}{\sum_i V_i} = \frac{\sum_i \bar p_i q_i}{\sum_i q_i}

any consistent units

Holds when: All wells of the block operate in the pseudo-steady-state regime at constant flow rates. Single-phase fluid; the compressibility and formation volume factor are the same throughout the block. Wells of different blocks separated by sealing faults are considered separately.

Pseudo-steady-state flow and the shape factor →
Approximate shape factor of a rectangular area
12ln⁡4AγCArw2≈2πLW(13−xD+xD2)+ln⁡W2πrwsin⁡(πyD)\frac12\ln\frac{4A}{\gamma C_A r_w^2} \approx \frac{2\pi L}{W}\left(\frac13 - x_D + x_D^2\right) + \ln\frac{W}{2\pi r_w\sin(\pi y_D)}

dimensionless; LL, WW and rwr_w in the same units; ln⁡\ln is the natural logarithm

Holds when: Pseudo-steady-state regime, homogeneous reservoir, the well penetrates the whole thickness. The well is at least the width WW away from both short sides. Then the error of the logarithmic term is below 0.01. For a well on the long axis (yD=1/2y_D = 1/2), half the width is enough. Closer to a short side the equation underestimates the resistance. For an off-axis well at a distance of W/2W/2 the error reaches 0.1–0.2. For a well in a corner of the area the equation is unsuitable, and tabulated values are needed.

Pseudo-steady-state flow and the shape factor →
Skin factor and the additional pressure drop
Δps=pw−pwf=qBμ2πkh S\Delta p_s = p_w - p_{wf} = \frac{qB\mu}{2\pi kh}\,S

SI (any consistent system of units)

Holds when: The zone of altered permeability is small compared with the drainage area, and the flow in it can be taken as steady-state at the current flow rate. No fluid builds up in the thin zone, so Δps\Delta p_s follows the flow rate immediately and disappears when the well is shut in. The flow is single-phase and obeys Darcy's law. If there are inertial losses, a term that grows with the flow rate is added to the skin. The skin is always referred to the same kk and hh that describe the inflow in the reservoir.

Skin factor →
Skin of a zone of altered permeability (Hawkins formula)
S=(kks−1)ln⁡rsrwS = \left(\frac{k}{k_s} - 1\right)\ln\frac{r_s}{r_w}

any consistent units: only the ratios k/ksk/k_s and rs/rwr_s/r_w enter

Holds when: The altered zone is a cylinder around the wellbore with uniform permeability ksk_s, and rs≪rer_s \ll r_e. The well penetrates the whole thickness of the reservoir, and the flow in the zone is radial. The flow is single-phase and obeys Darcy's law.

Skin factor →
Effective wellbore radius
rw′=rw e−Sr_w' = r_w\,\e^{-S}

any units of length (the same for rwr_w and rw′r_w')

Holds when: The skin does not change with time and does not depend on the flow rate. The zone that the effective radius replaces is small compared with the drainage area (rw′≪rer_w' \ll r_e).

Skin factor →
Radial inflow with skin
pe−pwf=qBμ2πkh(ln⁡rerw+S)pˉ−pwf=qBμ2πkh(ln⁡rerw−34+S)\begin{aligned} p_e - p_{wf} &= \frac{qB\mu}{2\pi kh}\left(\ln\frac{r_e}{r_w} + S\right)\\ \bar p - p_{wf} &= \frac{qB\mu}{2\pi kh}\left(\ln\frac{r_e}{r_w} - \frac34 + S\right)\end{aligned}

SI (any consistent system of units)

Holds when: A homogeneous reservoir, a circular drainage area with the well at the centre. The flow is single-phase and obeys Darcy's law. The well penetrates the whole thickness of the reservoir. Otherwise a component for imperfect penetration (the well opens only part of the layer) enters the skin. The second formula is for the pseudo-steady-state regime, when the pressure disturbance has already reached the boundaries of the closed area.

Skin factor →
Productivity index with skin (pseudo-steady-state flow)
J=qpˉ−pwf=2πkhBμ(ln⁡rerw−34+S)J = \frac{q}{\bar p - p_{wf}} = \frac{2\pi kh}{B\mu\left(\ln\dfrac{r_e}{r_w} - \dfrac34 + S\right)}

SI

Holds when: Pseudo-steady-state flow; in other regimes the logarithmic term changes, and the skin enters in the same way.

Skin factor →
Flow efficiency
FE=JJS=0=ln⁡rerw−34ln⁡rerw−34+S=pˉ−pwf−Δpspˉ−pwf\begin{aligned}\text{FE} &= \frac{J}{J_{S=0}} = \frac{\ln\dfrac{r_e}{r_w} - \dfrac34}{\ln\dfrac{r_e}{r_w} - \dfrac34 + S}\\ &= \frac{\bar p - p_{wf} - \Delta p_s}{\bar p - p_{wf}}\end{aligned}

dimensionless; pressures in any identical units

Holds when: Pseudo-steady-state flow; both productivities refer to the same drainage area. pˉ\bar p is the average pressure in the drainage area. For a new well we use instead the extrapolated pressure from the build-up curve (the pressure found by extending the build-up curve to an infinitely long shut-in).

Skin factor →
Components of the measured skin
S=Sd+Sgeom+Scompl+SopS = S_d + S_{\text{geom}} + S_{\text{compl}} + S_{\text{op}}

dimensionless quantities

Holds when: Simple addition is only an approximation, because the components affect one another. If only part of the thickness hwh_w is open, the damage skin, counted over the whole thickness hh, grows by roughly a factor of h/hwh/h_w. Shallow damage that the perforation tunnels pass right through hardly hinders the inflow.

Skin factor →
Skin pressure drop from test results
Δps=0.869 m S\Delta p_s = 0.869\,m\,S

any pressure units (the same as for mm)

Holds when: We get mm and SS from one test, on a section of radial flow. Δps\Delta p_s refers to the flow rate of that test.

Skin factor →
Radius of the filtrate invaded zone
rinv=rw2+Qlπ ϕ (1−Sor−Swc)r_{inv} = \sqrt{r_w^2 + \frac{Q_l}{\pi\,\phi\,(1 - S_{or} - S_{wc})}}

SI (m, m³/m)

Holds when: Piston-like displacement; the reservoir is uniform over its thickness. Usually rinvr_{inv} is much larger than rwr_w, and the term rw2r_w^2 can be left out.

Skin factor →
Pseudo-skin of partial penetration
Spp=(hhw−1)[ln⁡hD−G ⁣(hwh)],hD=hsrwkhkvG(x)=2.948−7.363 x+11.45 x2−4.675 x3\begin{gathered} S_{pp} = \left(\frac{h}{h_w} - 1\right)\left[\ln h_D - G\!\left(\frac{h_w}{h}\right)\right],\qquad h_D = \frac{h_s}{r_w}\sqrt{\frac{k_h}{k_v}}\\ G(x) = 2.948 - 7.363\,x + 11.45\,x^2 - 4.675\,x^3\end{gathered}

dimensionless; lengths and permeabilities in identical units; ln⁡\ln is the natural logarithm

Holds when: A homogeneous reservoir with constant khk_h and kvk_v between an impermeable top and bottom; a vertical well, one continuous interval. The formula approximates numerical solutions. For hDh_D from 100 to tens of thousands and hw/h≥0.1h_w/h \ge 0.1, it differs from the numerical solution with uniform pressure along the interval by no more than 0.6 (a few per cent), usually on the low side. For smaller hDh_D (a thin reservoir) it underestimates the pseudo-skin: at hD=20h_D = 20 and a small penetration ratio, by a third. Then a numerical calculation is needed. For an interval in an in-between position, the value lies between the values for an interval at the top and an interval in the middle. For greater accuracy, use a numerical model.

Partial penetration and well deviation →
Measured skin of a partially penetrating well
S=Sd hhw+Spp,Sd=(S−Spp) hwhS = S_d\,\frac{h}{h_w} + S_{pp},\qquad S_d = (S - S_{pp})\,\frac{h_w}{h}

dimensionless

Holds when: The damaged zone is thin: its radius is small compared with the interval and with the zone where the flow converges. The flow through the damaged zone is radial towards the interval. In a gas well an inertial skin is added to SS, and partial penetration also magnifies it (Inflow at high rates: non-Darcy flow).

Partial penetration and well deviation →
Pressure derivative in spherical and hemispherical flow
Δp′=c qBμϕμct8π3/2 ks3/2t,ks=(kh2 kv)1/3\Delta p' = \frac{c\,qB\mu\sqrt{\phi\mu c_t}}{8\pi^{3/2}\,k_s^{3/2}\sqrt{t}},\qquad k_s = \left(k_h^2\,k_v\right)^{1/3}

SI (any consistent system of units)

Holds when: A homogeneous anisotropic reservoir; the interval is short compared with the thickness. The time is later than the end of early radial flow and earlier than the start of late radial flow. Wellbore storage has ended.

Partial penetration and well deviation →
Duration of the regimes (estimate)
tearly≈a1 ϕμct hw2kv,tlate≈a2 ϕμct h2kvt_{\text{early}} \approx a_1\,\frac{\phi\mu c_t\,h_w^2}{k_v},\qquad t_{\text{late}} \approx a_2\,\frac{\phi\mu c_t\,h^2}{k_v}

metric oilfield units: tt in h, hh and hwh_w in m

Holds when: These are estimates from a model of uniform inflow along the interval, for a penetration ratio of up to 0.3. We take the boundaries of the regimes where the derivative departs from the plateau by 5%. If we allow 10% instead, the factor a1a_1 is four times larger and a2a_2 is about 20% smaller. Wellbore storage is not taken into account: it often hides early radial flow completely.

Partial penetration and well deviation →
Pseudo-skin of a deviated wellbore
Sθ=−(θ′41)2.06−(θ′56)1.865lg⁡hD100θ′=arctan⁡(kvkh tan⁡θ),hD=hrwkhkv\begin{gathered} S_\theta = -\left(\frac{\theta'}{41}\right)^{2.06} - \left(\frac{\theta'}{56}\right)^{1.865}\lg\frac{h_D}{100}\\ \theta' = \arctan\left(\sqrt{\frac{k_v}{k_h}}\,\tan\theta\right),\qquad h_D = \frac{h}{r_w}\sqrt{\frac{k_h}{k_v}}\end{gathered}

dimensionless; lg⁡\lg is the decimal logarithm, angles in degrees

Holds when: The wellbore opens and perforates the whole thickness of a homogeneous reservoir. The limits are θ′≤75∘\theta' \le 75^\circ and hD>40h_D > 40. In an isotropic reservoir θ′=θ\theta' = \theta. Allowing for anisotropy through the effective angle is a common approximation. The value refers to late, pseudo-radial flow. Other models describe horizontal wells (Horizontal wells).

Partial penetration and well deviation →
Measured skin of a deviated well
S=Sdcos⁡θ′+Sθ,Sd=S−Sθcos⁡θ′S = S_d\cos\theta' + S_\theta,\qquad S_d = \frac{S - S_\theta}{\cos\theta'}

dimensionless

Holds when: The damaged zone is thin and runs along the whole wellbore in the reservoir. The wellbore is perforated along its whole length in the reservoir. cos⁡θ′\cos\theta' equals cos⁡θ\cos\theta (the ratio of the thickness to the length of the wellbore in the reservoir) multiplied by the ratio of khk_h to the effective permeability across the wellbore. In an anisotropic reservoir that permeability is smaller than khk_h, so cos⁡θ′>cos⁡θ\cos\theta' > \cos\theta.

Partial penetration and well deviation →
Pseudo-skin for uniform inflow along the interval
Spp=2π2x2∑n=1∞[sin⁡(nπz2)−sin⁡(nπz1)]2n2 K0 ⁣(nπrwhkvkh)S_{pp} = \frac{2}{\pi^2 x^2}\sum_{n=1}^{\infty}\frac{\left[\sin(n\pi z_2) - \sin(n\pi z_1)\right]^2}{n^2}\,K_0\!\left(\frac{n\pi r_w}{h}\sqrt{\frac{k_v}{k_h}}\right)

dimensionless

Holds when: A homogeneous anisotropic reservoir. A line source with uniform inflow along the interval. The pressure is averaged over the interval.

Partial penetration and well deviation →
Pressure derivative of a partially penetrating well
pD′=12 e−1/(4tD) F(tD),F(tD)=1+∑n=1∞an22 e−n2π2tD/hD2p_D' = \tfrac12\,e^{-1/(4t_D)}\,F(t_D),\qquad F(t_D) = 1 + \sum_{n=1}^{\infty}\frac{a_n^2}{2}\,e^{-n^2\pi^2 t_D/h_D^2}

dimensionless

Holds when: A line source with uniform inflow along the interval; a reservoir unbounded in plan; no wellbore storage.

Partial penetration and well deviation →
Critical drawdown of a stationary cone
pe−pw=Δρ g hc,Δρ=ρw−ρop_e - p_w = \Delta\rho\,g\,h_c,\qquad \Delta\rho = \rho_w - \rho_o

SI (any consistent system of units)

Holds when: The cone is stationary, and the water beneath it is in hydrostatic equilibrium. The interface is sharp: we ignore the capillary transition zone (the zone where oil and water change over gradually). We take the pressure at the tip of the cone to be equal to the bottomhole pressure. In fact the tip is below the interval, and the pressure there is higher than the bottomhole pressure. So the estimate understates the permissible drawdown, which errs on the safe side. For a gas cone Δρ=ρo−ρg\Delta\rho = \rho_o - \rho_g, and hch_c is the distance from the GOC to the top perforations.

Water and gas coning →
Critical water-free rate from hydrostatic equilibrium
qc=2πkhh Δρ g hcμB(ln⁡rerw+Spp)q_c = \frac{2\pi k_h h\,\Delta\rho\,g\,h_c}{\mu B\left(\ln\dfrac{r_e}{r_w} + S_{pp}\right)}

SI; ln⁡\ln is the natural logarithm

Holds when: The inflow of single-phase oil to a partially penetrating well is steady-state, in a homogeneous reservoir. The pressure at the outer boundary is maintained. The water is stationary, and the cone does not distort the pressure field in the oil. The tip “feels” the whole drawdown at the bottomhole, including the resistance of the flow converging on the interval. This is why the result is too low. Damage skin does not enter the formula, because the pressure drop in the thin damaged zone at the wellbore does not reach the formation below the interval.

Water and gas coning →
Critical rate of a cone at vertical equilibrium (Meyer–Garder)
qc=πkh Δρ g (h2−hw2)μBln⁡rerwq_c = \frac{\pi k_h\,\Delta\rho\,g\,(h^2 - h_w^2)}{\mu B\ln\dfrac{r_e}{r_w}}

SI; ln⁡\ln is the natural logarithm

Holds when: Vertical equilibrium: the pressure in the oil is constant along the vertical, and the oil flows horizontally. Close to the well this assumption is rough. The reservoir is isotropic, so kvk_v does not enter the formula. The resistance of the flow converging on the interval is not taken into account. The interval is open from the top. At the critical rate the thickness of oil above the cone at the wellbore equals the interval hwh_w.

Water and gas coning →
Critical rate in an anisotropic reservoir (Chaperon's approximation)
qc=kh Δρ g (h−hw)2μB F,F=0.7311+1.943reD,reD=rehkvkhq_c = \frac{k_h\,\Delta\rho\,g\,(h - h_w)^2}{\mu B}\,F,\qquad F = 0.7311 + \frac{1.943}{r_{eD}},\qquad r_{eD} = \frac{r_e}{h}\sqrt{\frac{k_v}{k_h}}

SI (any consistent system of units)

Holds when: A homogeneous anisotropic reservoir. Steady-state inflow, with a constant pressure at the outer boundary. Stationary water. An interval open from the top. The formula approximates theoretical calculations of a stable cone for a well that is open at the top of the reservoir. FF grows as reDr_{eD} decreases, so a lower vertical permeability raises the critical rate.

Water and gas coning →
Dimensionless cone height and dimensionless time
Z=2πkh Δρ g h hcq μB,tD=kv Δρ g(1+M′ α)2 ϕ μ h tZ = \frac{2\pi k_h\,\Delta\rho\,g\,h\,h_c}{q\,\mu B},\qquad t_D = \frac{k_v\,\Delta\rho\,g\left(1 + M'^{\,\alpha}\right)}{2\,\phi\,\mu\,h}\,t

SI (any consistent system of units)

Holds when: The groups and the exponent α\alpha come from a correlation built on model experiments and calculations, for a homogeneous reservoir with bottom water and an interval open from the top. In oilfield units (Δρ\Delta\rho in g/cm³, hh in ft, kk in mD, qq in STB/d, tt in days) the factors are 0.00307 (for ZZ) and 0.00137 (for tDt_D). Some reference books write the dimensionless time with the same factor as ZZ (0.00307 instead of 0.00137). Then the breakthrough time comes out 2.2 times too small. The original definition gives 0.00137, which matches the form in consistent units above. If you use other units, work out the factor again from that form. A check of the units will not catch this error.

Water and gas coning →
Dimensionless breakthrough time (Bournazel–Jeanson correlation)
tD,BT=Z3−0.7Z,Z<30.7≈4.29t_{D,BT} = \frac{Z}{3 - 0.7Z},\qquad Z < \frac{3}{0.7} \approx 4.29

dimensionless

Holds when: This is an empirical correlation of model experiments. An earlier correlation with the same groups (Sobocinski–Cornelius) is tD,BT=(4Z+1.75Z2−0.75Z3)/(7−2Z)t_{D,BT} = (4Z + 1.75Z^2 - 0.75Z^3)/(7 - 2Z) for Z<3.5Z < 3.5. For ZZ above 4.29 (3.5 in the second correlation) the water does not reach the interval. This is one more, “dynamic”, estimate of the critical rate. The reservoir is homogeneous. Layers, fractures and vertical heterogeneity change the result several times over, so the accuracy is at best about a factor of two.

Water and gas coning →
Limiting water cut with a developed cone
fw,lim⁡=M′haqM′haq+hf_{w,\lim} = \frac{M' h_{aq}}{M' h_{aq} + h}

dimensionless; thicknesses in the same units

Holds when: Both zones deliver liquid to the well side by side, in proportion to k krh/μk\,k_r h/\mu. The pressure and the permeability are the same in both zones. As fluid is produced, the OWC rises. So haqh_{aq} grows, hh decreases, and the limit rises with time. The water fraction is at reservoir conditions. At standard conditions the water cut is somewhat higher, because when oil is brought to the surface its volume shrinks more than the volume of water (the oil formation volume factor is greater).

Water and gas coning →
Growth of water cut after breakthrough (Kuo–DesBrisay correlation)
fwfw,lim⁡=0.29+0.94 lg⁡ttBT,0.5≤ttBT≤5.7\frac{f_w}{f_{w,\lim}} = 0.29 + 0.94\,\lg\frac{t}{t_{BT}},\qquad 0.5 \le \frac{t}{t_{BT}} \le 5.7

dimensionless

Holds when: At t/tBT<0.5t/t_{BT} < 0.5 the ratio is zero, and at t/tBT>5.7t/t_{BT} > 5.7 it is 1. The reservoir is homogeneous with bottom water, and the liquid rate is constant. Formally the water appears already at half the calculated breakthrough time. This comes from how the correlation is scaled to the estimate of tBTt_{BT}, and it is one more reminder of how accurate that estimate is. We allow separately for the rise of the limit fw,lim⁡f_{w,\lim} as the OWC rises, by recalculating hh and haqh_{aq} from the cumulative production (the total produced so far).

Water and gas coning →
Critical rate and position of the interval in an oil rim
qc=πkh Δρeff g (h2−hw2)μBln⁡(re/rw),Δρeff=(ρw−ρo)(ρo−ρg)ρw−ρghg=(h−hw) ρw−ρoρw−ρg\begin{gathered} q_c = \frac{\pi k_h\,\Delta\rho_{\text{eff}}\,g\,(h^2 - h_w^2)}{\mu B\ln(r_e/r_w)},\qquad \Delta\rho_{\text{eff}} = \frac{(\rho_w - \rho_o)(\rho_o - \rho_g)}{\rho_w - \rho_g}\\ h_g = (h - h_w)\,\frac{\rho_w - \rho_o}{\rho_w - \rho_g}\end{gathered}

SI; ln⁡\ln is the natural logarithm

Holds when: Vertical equilibrium and an isotropic reservoir, as in the formula with one cone. Both phases, below and above the oil, are stationary. We find the position of the interval from the stability condition. We do not take into account the movement of the contacts during development (the cap expands, the water rises).

Water and gas coning →
Forchheimer equation
− dp dx=μk v+βρ v2-\frac{\dd p}{\dd x} = \frac{\mu}{k}\,v + \beta\rho\,v^2

SI (any consistent system of units)

Holds when: The flow is single-phase (one fluid) and isothermal (constant temperature). At every point the velocity points along the pressure gradient (the direction in which pressure changes fastest). β\beta is a property of the rock, like kk. We find it by experiment (see below). As v→0v \to 0, the second term vanishes and the equation becomes Darcy's law.

Inflow at high rates: non-Darcy flow →
Forchheimer number
Fo=βρv2μv/k=kβρvμ\text{Fo} = \frac{\beta\rho v^2}{\mu v/k} = \frac{k\beta\rho v}{\mu}

no units; the quantities are in consistent units

Holds when: This is an estimate at one point of the flow. Near a well, vv and Fo change with distance.

Inflow at high rates: non-Darcy flow →
Gas test on core: the linear form of the Forchheimer equation
M (p12−p22)2ZRT L μ G=1k+β Gμ\frac{M\,(p_1^2 - p_2^2)}{2ZRT\,L\,\mu\,G} = \frac{1}{k} + \beta\,\frac{G}{\mu}

SI; pressures absolute

Holds when: The flow is isothermal and single-phase. ZZ and μ\mu do not change over the range of pressures in the test. Gas slippage (the Klinkenberg correction: gas molecules slip along the pore walls, so the rock looks more permeable than it is) is small. To make sure of this, the test is run with back pressure (pressure held at the outlet), so that the mean pressure in the sample is about a megapascal. Otherwise slippage is allowed for separately (Rock permeability).

Inflow at high rates: non-Darcy flow →
Correlation of the inertial resistance coefficient with permeability
β=7.64⋅1010k1.201\beta = \frac{7.64\cdot10^{10}}{k^{1.201}}

empirical: β\beta in m⁻¹, kk in mD

Holds when: The rocks are consolidated and the gas flow is single-phase. Measurements scatter around the correlation by up to an order of magnitude (a factor of ten). Irreducible water (water held in the rock that cannot be produced) or condensate that has dropped out narrows the channels. Then β\beta for gas increases, sometimes several times over.

Inflow at high rates: non-Darcy flow →
Steady-state radial flow with inertial losses
pe−pw=qμ2πkhln⁡rerw+βρq2(2πh)2(1rw−1re)p_e - p_w = \frac{q\mu}{2\pi kh}\ln\frac{r_e}{r_w} + \frac{\beta\rho q^2}{(2\pi h)^2}\left(\frac{1}{r_w} - \frac{1}{r_e}\right)

SI; ln⁡\ln is the natural logarithm

Holds when: The inflow is steady-state and single-phase, and the liquid is slightly compressible (we treat gas below, through the mass flow rate). The reservoir is homogeneous (the same everywhere). The well is an open hole (not cased or perforated) over the whole thickness. β\beta and ρ\rho are constant.

Inflow at high rates: non-Darcy flow →
Shares of viscous and inertial losses within radius rr
δvisc(r)=ln⁡(r/rw)ln⁡(re/rw),δinert(r)=1/rw−1/r1/rw−1/re≈1−rwr\delta_{\text{visc}}(r) = \frac{\ln(r/r_w)}{\ln(r_e/r_w)},\qquad \delta_{\text{inert}}(r) = \frac{1/r_w - 1/r}{1/r_w - 1/r_e} \approx 1 - \frac{r_w}{r}

no units; radii in the same units

Holds when: Steady-state radial flow, as in the previous formula. The approximate equality holds when re≫rwr_e \gg r_w.

Inflow at high rates: non-Darcy flow →
Rate-dependent skin
pe−pwf=qμ2πkh(ln⁡rerw+S+Dq),D=βρk2πhμrwp_e - p_{wf} = \frac{q\mu}{2\pi kh}\left(\ln\frac{r_e}{r_w} + S + Dq\right),\qquad D = \frac{\beta\rho k}{2\pi h\mu r_w}

SI; ln⁡\ln is the natural logarithm

Holds when: re≫rwr_e \gg r_w, so the term 1/re1/r_e is dropped. The well is an open hole over the whole thickness, and β\beta at the wellbore is the same as in the reservoir. With damage and with partial penetration, DD is larger (see the last section).

Inflow at high rates: non-Darcy flow →
Non-Darcy coefficient of a gas well
D=β ρsc k2πh μw rwD = \frac{\beta\,\rho_{sc}\,k}{2\pi h\,\mu_w\,r_w}

SI

Holds when: The mass flow rate is the same at all radii, and the inflow is isothermal. The inertial losses are concentrated at the wellbore, so the viscosity is taken at the bottomhole pressure. The well is an open hole over the whole thickness, and the near-wellbore zone is undamaged.

Inflow at high rates: non-Darcy flow →
Two-term gas inflow equation in pseudo-pressures
m(pˉ)−m(pwf)=a q+b q2a=U(ln⁡rerw−34+S),b=UDU=psc Tπkh Tsc\begin{gathered} m(\bar p) - m(p_{wf}) = a\,q + b\,q^2\\ a = U\left(\ln\frac{r_e}{r_w} - \frac34 + S\right),\qquad b = U D\\ U = \frac{p_{sc}\,T}{\pi kh\,T_{sc}}\end{gathered}

SI; pressures absolute; ln⁡\ln is the natural logarithm

Holds when: The inflow is pseudo-steady-state, in a circular area with the well at the centre. For another shape of area we change the logarithmic term, as in the lesson on pseudo-steady-state flow. Using m(pˉ)m(\bar p) instead of the volume-averaged pseudo-pressure is the usual approximation. DD is taken with the gas viscosity at the bottom of the well. The gas in the reservoir is single-phase.

Inflow at high rates: non-Darcy flow →
Rate of a gas well from the two-term equation
q=−a+a2+4b Δm2b,Δm=m(pˉ)−m(pwf)q = \frac{-a + \sqrt{a^2 + 4b\,\Delta m}}{2b},\qquad \Delta m = m(\bar p) - m(p_{wf})

metric oilfield (or any units consistent with aa and bb)

Holds when: The inflow is stabilised (the pressures no longer change with time) and the coefficients aa and bb are known. For b=0b = 0 the formula becomes q=Δm/aq = \Delta m/a.

Inflow at high rates: non-Darcy flow →
Apparent skin of the flow periods
Sj′=S+D qjS'_j = S + D\,q_j

any; DD is in units that are the inverse of the rate units

Holds when: All flow periods are interpreted with one reservoir model, and the flow capacity khkh is the same for all of them. The skin and DD do not change during the test (the well has already been cleaned up after completion).

Inflow at high rates: non-Darcy flow →
Apparent skin with a damaged zone
S′=(kks−1)ln⁡rsrw+Dq,D=D0[βsβ(1−rwrs)+rwrs]S' = \left(\frac{k}{k_s} - 1\right)\ln\frac{r_s}{r_w} + Dq,\qquad D = D_0\left[\frac{\beta_s}{\beta}\left(1 - \frac{r_w}{r_s}\right) + \frac{r_w}{r_s}\right]

any consistent units; DD and D0D_0 in the same units

Holds when: The damaged zone is a cylinder around the wellbore, and the flow in it is steady-state. The first term is the Hawkins formula (Skin factor). re≫rsr_e \gg r_s, and the hole is open over the whole thickness.

Inflow at high rates: non-Darcy flow →
Apparent skin with partial penetration
S′=Sd hhw+Spp+D0(hhw)2qS' = S_d\,\frac{h}{h_w} + S_{pp} + D_0\left(\frac{h}{h_w}\right)^{2} q

no units; D0qD_0 q in consistent units

Holds when: The damaged zone is thin compared with the interval, and the flow near the wellbore is radial towards the interval. The inertial losses in the zone of flow convergence, farther from the wellbore, are small compared with the losses at the wall.

Inflow at high rates: non-Darcy flow →
Permeability and the β coefficient of a packed bed of grains
k=ϕ3dp2150 (1−ϕ)2,β=1.75 (1−ϕ)ϕ3dp=1.75 dp150 (1−ϕ) kk = \frac{\phi^3 d_p^2}{150\,(1 - \phi)^2},\qquad \beta = \frac{1.75\,(1 - \phi)}{\phi^3 d_p} = \frac{1.75\,d_p}{150\,(1 - \phi)\,k}

SI

Holds when: The pack is unconsolidated (loose), and the grains are of similar size. The constants 150 and 1.75 are empirical (found from experiments).

Inflow at high rates: non-Darcy flow →

Symbols

SymbolMeaning and context
pˉ\bar p
Δm\Delta m
pˉ−pwf\bar p - p_{wf}
pˉ, pwf\bar p,\ p_{wf}
pˉ, pwf, q, B, μ, k, h, S, rw\bar p,\ p_{wf},\ q,\ B,\ \mu,\ k,\ h,\ S,\ r_w
Δp′\Delta p'
AA
a1, a2a_1,\ a_2
ana_n
  • 2[sin⁡(nπz2)−sin⁡(nπz1)]/(nπx)2[\sin(n\pi z_2) - \sin(n\pi z_1)]/(n\pi x) — the amplitudes of the harmonics of the source, —Partial penetration and well deviation
a, ba,\ b
BB
  • fluid formation volume factor (reservoir volume per unit of volume at the surface); it appears in the oilfield forms, where the flow rate is at standard conditions, m³/m³Steady-state radial flow
  • fluid formation volume factor, m³/m³Steady-state radial flow
  • formation volume factor of the fluid, m³/m³Skin factor
BoB_o
Bo, Bw, BgB_o,\ B_w,\ B_g
BwB_w
B, μB,\ \mu
B, μ, k, h, re, rwB,\ \mu,\ k,\ h,\ r_e,\ r_w
pˉblk\bar p_{\text{blk}}
cc
CAC_A
CA, γC_A,\ \gamma
CA, γ, rwC_A,\ \gamma,\ r_w
ctc_t
DD
D0D_0
dpd_p
DqDq
FF
fwf_w
fw,lim⁡f_{w,\lim}
FE\text{FE}
Fo\text{Fo}
gg
GG
hh
HH
haqh_{aq}
  • thickness of the water-saturated part of the reservoir below the oil within the drainage area, mWater and gas coning
hch_c
  • distance from the bottom perforations of the interval to the original OWC: the critical height of the cone, mWater and gas coning
  • distance from the bottom perforations to the OWC; for an interval open from the top hc=h−hwh_c = h - h_w, mWater and gas coning
hDh_D
hgh_g
  • distance from the GOC to the top perforations; h−hw−hgh - h_w - h_g remains to the OWC, mWater and gas coning
hsh_s
hwh_w
hw/hh_w/h
h−hwh - h_w
h, hwh,\ h_w
h, rw, kh, kvh,\ r_w,\ k_h,\ k_v
pˉi\bar p_i
JJ
J, JS=0J,\ J_{S=0}
  • actual productivity and productivity of the same well without skin, m³/(d·MPa)Skin factor
kk
KK
K0K_0
khk_h
kh, g, μ, B, re, rwk_h,\ g,\ \mu,\ B,\ r_e,\ r_w
kh, Δρ, g, h, hc, μ, Bk_h,\ \Delta\rho,\ g,\ h,\ h_c,\ \mu,\ B
kh, Δρ, g, μ, B, rek_h,\ \Delta\rho,\ g,\ \mu,\ B,\ r_e
kh, Δρ, g, μ, B, re, rwk_h,\ \Delta\rho,\ g,\ \mu,\ B,\ r_e,\ r_w
kh, kvk_h,\ k_v
ksk_s
ks, βsk_s,\ \beta_s
kvk_v
k, hk,\ h
k, h, B, μ, re, rw, Sk,\ h,\ B,\ \mu,\ r_e,\ r_w,\ S
  • as in the inflow formula above, m², m, m³/m³, Pa·s, m, m, —Skin factor
k, h, B, μ, re, rw, pek,\ h,\ B,\ \mu,\ r_e,\ r_w,\ p_e
k, h, re, rwk,\ h,\ r_e,\ r_w
k, h, μk,\ h,\ \mu
k, h, μ, re, rw, pe, pwk,\ h,\ \mu,\ r_e,\ r_w,\ p_e,\ p_w
k, βk,\ \beta
k, β, ρ, v, μk,\ \beta,\ \rho,\ v,\ \mu
k, ϕ, μ, ctk,\ \phi,\ \mu,\ c_t
LL
L, WL,\ W
lg⁡\lg
mm
  • slope of the semilog straight line per decimal logarithmic cycle (absolute value), MPa/cycleSkin factor
MM
M′M'
  • end-point mobility ratio (krw′/μw)/(kro′/μo)(k_{rw}'/\mu_w)/(k_{ro}'/\mu_o): water over oil, each as end-point relative permeability divided by viscosity (Relative permeability), —Water and gas coning
  • end-point mobility ratio of water to oil, —Water and gas coning
m(p)m(p)
p0p_0
p1, p2p_1,\ p_2
pD′p_D'
  • dimensionless derivative, pD=2πkhh Δp/(qBμ)p_D = 2\pi k_h h\,\Delta p/(qB\mu) (over the whole thickness, without skin), —Partial penetration and well deviation
pep_e
pe, pwp_e,\ p_w
pip_i
pwp_w
pw, pep_w,\ p_e
pwfp_{wf}
pwf(t)p_{wf}(t)
pwhp_{wh}
p(r,t)p(r,t)
p(r)p(r)
qq
qaqq_{aq}
qcq_c
qi, qtq_i,\ q_t
qjq_j
QlQ_l
  • total volume of filtrate that has entered, per unit length of wellbore, m³/mSkin factor
qmaxq_{max}
  • absolute open flow: qq at Δmmax=m(pˉ)−m(patm)\Delta m_{max} = m(\bar p) - m(p_{\text{atm}}), thousand m³/dInflow at high rates: non-Darcy flow
qoq_o
qo, qwq_o,\ q_w
qsq_s
qs, B, μ, k, h, Sq_s,\ B,\ \mu,\ k,\ h,\ S
  • flow rate at standard conditions, formation volume factor, viscosity, permeability, thickness, skin, m³/s, m³/m³, Pa·s, m², m, —Steady-state radial flow
qscq_{sc}
qw,jq_{w,j}
q, B, μq,\ B,\ \mu
q, B, μ, k, h, Sq,\ B,\ \mu,\ k,\ h,\ S
q, B, μ, k, h, A, CA, rw, Sq,\ B,\ \mu,\ k,\ h,\ A,\ C_A,\ r_w,\ S
q, pˉ, pwfq,\ \bar p,\ p_{wf}
  • flow rate at standard conditions, average reservoir pressure and bottomhole flowing pressure, m³/s, PaSkin factor
q, μ, k, h, β, ρ, re, rwq,\ \mu,\ k,\ h,\ \beta,\ \rho,\ r_e,\ r_w
q, μ, k, h, pw, rwq,\ \mu,\ k,\ h,\ p_w,\ r_w
q(r)q(r)
rr
RR
rer_e
re, rwr_e,\ r_w
re, rw, Sr_e,\ r_w,\ S
  • radius of the drainage area, radius of the well, skin factor, m, m, —Skin factor
re=A/πr_e = \sqrt{A/\pi}
reDr_{eD}
rinvr_{inv}
RpR_p
rsr_s
RsR_s
rwr_w
rw′r_w'
  • effective radius: the radius of an ideal well without skin that has the same bottomhole pressure at the same flow rate, mSkin factor
r, rer,\ r_e
Δps\Delta p_s
  • additional pressure drop in the near-wellbore zone; positive for damage, negative for improvement, PaSkin factor
  • skin pressure drop at the current flow rate, MPaSkin factor
  • skin pressure drop at the flow rate of the test, MPaSkin factor
SS
SAS_A
ScomplS_{\text{compl}}
  • completion: the number, length and placement of the perforation tunnels, the crushed zone around them, a gravel pack (a filter of gravel placed around the well screen), —Skin factor
SdS_d
SgeomS_{\text{geom}}
  • geometry of penetration: only part of the thickness is open, a deviated wellbore, a well displaced towards the boundary of the drainage area, —Skin factor
SopS_{\text{op}}
  • operating effects: inertial losses that grow with the flow rate, and two-phase flow near the wellbore, —Skin factor
SorS_{or}
SppS_{pp}
SwcS_{wc}
SθS_\theta
S, DS,\ D
Sj′S'_j
tt
TT
tBTt_{BT}
tDt_D
tD, tt_D,\ t
tD,BTt_{D,BT}
  • dimensionless breakthrough time; we get the time with units from the definition of tDt_D, —Water and gas coning
tDAt_{DA}
tearlyt_{\text{early}}
tlatet_{\text{late}}
Tsc, pscT_{sc},\ p_{sc}
T, Tsc, pscT,\ T_{sc},\ p_{sc}
UU
vv
VV
ViV_i
Vp=ϕhAV_p = \phi hA
VtV_t
xx
xDx_D
x=z2−z1x = z_2 - z_1
yDy_D
ZZ
z1, z2z_1,\ z_2
α\alpha
  • 0.5 for M′≤1M' \le 1 and 0.6 for 1<M′≤101 < M' \le 10 (in the Bournazel–Jeanson correlation discussed below it is 0.7 for 0.14≤M′≤7.30.14 \le M' \le 7.3), —Water and gas coning
β\beta
β, k, h, rw\beta,\ k,\ h,\ r_w
β, ρ\beta,\ \rho
γ\gamma
δvisc, δinert\delta_{\text{visc}},\ \delta_{\text{inert}}
  • the viscous and inertial pressure losses between rwr_w and rr, each as a share of the corresponding losses between rwr_w and rer_e, —Inflow at high rates: non-Darcy flow
Δρ, g\Delta\rho,\ g
  • density difference between water and oil at reservoir conditions, and the acceleration due to gravity, kg/m³, m/s²Water and gas coning
Δρeff\Delta\rho_{\text{eff}}
θ\theta
θ′\theta'
μ\mu
μ, B\mu,\ B
μ, Z\mu,\ Z
  • gas viscosity and gas deviation factor (the correction that makes the ideal gas law fit a real gas; at the pressures of the test Z≈1Z \approx 1), Pa·s, —Inflow at high rates: non-Darcy flow
μˉ, Zˉ\bar\mu,\ \bar Z
  • gas viscosity and gas deviation factor at the average pressure (pe+pw)/2(p_e + p_w)/2, Pa·s, —Steady-state radial flow
μ(p), Z(p)\mu(p),\ Z(p)
  • gas viscosity and gas deviation factor as functions of pressure at reservoir temperature, Pa·s, —Steady-state radial flow
μw\mu_w
ρ\rho
ρsc\rho_{sc}
  • gas density at standard conditions: 1.204 γg1.204\,\gamma_g at 20 °C (γg\gamma_g is the gas specific gravity relative to air, that is, the gas density divided by the air density), kg/m³Inflow at high rates: non-Darcy flow
ρw, ρo\rho_w,\ \rho_o
ρw, ρo, ρg\rho_w,\ \rho_o,\ \rho_g
ϕ\phi
ϕ, h, A\phi,\ h,\ A
ϕ, μ, ct, kv\phi,\ \mu,\ c_t,\ k_v
ϕ, ct\phi,\ c_t