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ˉ−Jq
any consistent units, for example m³/d and MPa
Holds when: Single-phase inflow (pwf above the bubble point pressure) and stabilised flow. The productivity index does not depend on the flow rate: there are no inertial losses.
pressures in MPa; the product ρgH comes out in pascals and is divided by 106
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). K 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.
Injection balance under reservoir pressure maintenance
j∑qw,jBw+qaq=qoBo
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.
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 at standard conditions.
SI (any consistent system of units); ln is the natural logarithm
Holds when: Steady state: a constant pressure pe 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).
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).
SI; absolute pressures only; 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 changes little between pw and pe, 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.
Holds when: Constant flow rate after start-up; a homogeneous reservoir. r is much smaller than the radius of investigation 4kt/(ϕμct) (how far into the reservoir the disturbance has travelled): the disturbance reached this point long ago.
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).
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.
Pressure distribution in pseudo-steady-state inflow
p(r)−pwf=2πkhqBμ(lnrwr−2re2r2+S)
SI; ln is the natural logarithm
Holds when: Pseudo-steady-state regime, circular area with the well at the centre. rw≪re: terms of the order of (rw/re)2 are left out. The skin is an additional pressure drop in a thin zone at the wellbore (Skin factor).
SI (any consistent system of units); 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 inflow in an area of arbitrary shape
pˉ−pwf=2πkhqBμ(21lnγCArw24A+S)
SI; A and rw2 in the same units of area; ln is the natural logarithm
Holds when: Pseudo-steady-state regime: the dimensionless time tDA 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.
Bottomhole pressure in the pseudo-steady-state regime in terms of the initial pressure
pi−pwf=2πkhqBμ(21lnγCArw24A+S+2πtDA)
SI; 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.
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.
dimensionless; L, W and rw in the same units; 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 W 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/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/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.
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 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 k and h that describe the inflow in the reservoir.
Skin of a zone of altered permeability (Hawkins formula)
S=(ksk−1)lnrwrs
any consistent units: only the ratios k/ks and rs/rw enter
Holds when: The altered zone is a cylinder around the wellbore with uniform permeability ks, and rs≪re. 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.
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′≪re).
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.
Holds when: Pseudo-steady-state flow; both productivities refer to the same drainage area. 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).
Holds when: Simple addition is only an approximation, because the components affect one another. If only part of the thickness hw is open, the damage skin, counted over the whole thickness h, grows by roughly a factor of h/hw. Shallow damage that the perforation tunnels pass right through hardly hinders the inflow.
Holds when: Piston-like displacement; the reservoir is uniform over its thickness. Usually rinv is much larger than rw, and the term rw2 can be left out.
dimensionless; lengths and permeabilities in identical units; ln is the natural logarithm
Holds when: A homogeneous reservoir with constant kh and kv between an impermeable top and bottom; a vertical well, one continuous interval. The formula approximates numerical solutions. For hD from 100 to tens of thousands and hw/h≥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 hD (a thin reservoir) it underestimates the pseudo-skin: at hD=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.
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 S, and partial penetration also magnifies it (Inflow at high rates: non-Darcy flow).
Pressure derivative in spherical and hemispherical flow
Δp′=8π3/2ks3/2tcqBμϕμct,ks=(kh2kv)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.
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 a1 is four times larger and a2 is about 20% smaller. Wellbore storage is not taken into account: it often hides early radial flow completely.
dimensionless; 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∘ and hD>40. In an isotropic reservoir θ′=θ. 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).
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θ′ equals cosθ (the ratio of the thickness to the length of the wellbore in the reservoir) multiplied by the ratio of kh to the effective permeability across the wellbore. In an anisotropic reservoir that permeability is smaller than kh, so cosθ′>cosθ.
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, and hc is the distance from the GOC to the top perforations.
Critical water-free rate from hydrostatic equilibrium
qc=μB(lnrwre+Spp)2πkhhΔρghc
SI; 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.
Critical rate of a cone at vertical equilibrium (Meyer–Garder)
qc=μBlnrwreπkhΔρg(h2−hw2)
SI; 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 kv 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 hw.
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. F grows as reD decreases, so a lower vertical permeability raises the critical rate.
Holds when: The groups and the exponent α 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 (Δρ in g/cm³, h in ft, k in mD, q in STB/d, t in days) the factors are 0.00307 (for Z) and 0.00137 (for tD). Some reference books write the dimensionless time with the same factor as Z (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.
Dimensionless breakthrough time (Bournazel–Jeanson correlation)
tD,BT=3−0.7ZZ,Z<0.73≈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) for Z<3.5. For Z 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.
Holds when: Both zones deliver liquid to the well side by side, in proportion to kkrh/μ. The pressure and the permeability are the same in both zones. As fluid is produced, the OWC rises. So haq grows, h 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).
Growth of water cut after breakthrough (Kuo–DesBrisay correlation)
fw,limfw=0.29+0.94lgtBTt,0.5≤tBTt≤5.7
dimensionless
Holds when: At t/tBT<0.5 the ratio is zero, and at t/tBT>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 tBT, and it is one more reminder of how accurate that estimate is. We allow separately for the rise of the limit fw,lim as the OWC rises, by recalculating h and haq from the cumulative production (the total produced so far).
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).
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). β is a property of the rock, like k. We find it by experiment (see below). As v→0, the second term vanishes and the equation becomes Darcy's law.
Gas test on core: the linear form of the Forchheimer equation
2ZRTLμGM(p12−p22)=k1+βμG
SI; pressures absolute
Holds when: The flow is isothermal and single-phase. Z and μ 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).
Correlation of the inertial resistance coefficient with permeability
β=k1.2017.64⋅1010
empirical: β in m⁻¹, k 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 β for gas increases, sometimes several times over.
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. β and ρ are constant.
Holds when:re≫rw, so the term 1/re is dropped. The well is an open hole over the whole thickness, and β at the wellbore is the same as in the reservoir. With damage and with partial penetration, D is larger (see the last section).
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.
SI; pressures absolute; 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ˉ) instead of the volume-averaged pseudo-pressure is the usual approximation. D is taken with the gas viscosity at the bottom of the well. The gas in the reservoir is single-phase.
metric oilfield (or any units consistent with a and b)
Holds when: The inflow is stabilised (the pressures no longer change with time) and the coefficients a and b are known. For b=0 the formula becomes q=Δm/a.
any; D 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 kh is the same for all of them. The skin and D do not change during the test (the well has already been cleaned up after completion).
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≫rs, and the hole is open over the whole thickness.
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.
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
1 for spherical flow (an interval in the middle of the reservoir), 2 for hemispherical flow (an interval at the top or bottom), —Partial penetration and well deviation
thickness of the water-saturated part of the reservoir below the oil within the drainage area, mWater and gas coning
hc
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−hw, mWater and gas coning
hD
dimensionless thickness of the symmetry element (the part of the reservoir that repeats by symmetry), with allowance for anisotropy, —Partial penetration and well deviation
distance from the GOC to the top perforations; h−hw−hg remains to the OWC, mWater and gas coning
hs
thickness of the symmetry element: h for an interval at the top or bottom, h/2 for an interval in the middle of the reservoir, mPartial penetration and well deviation
end-point mobility ratio (krw′/μw)/(kro′/μo): water over oil, each as end-point relative permeability divided by viscosity (Relative permeability), —Water and gas coning
water influx from the aquifer beyond the oil–water contour (the water-bearing rock outside the oil zone), at reservoir conditions, m³/dWell inflow: the well performance diagram
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
Sd
damage skin: a change in the permeability of the near-wellbore zone, —Skin factor
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
Sop
operating effects: inertial losses that grow with the flow rate, and two-phase flow near the wellbore, —Skin factor
the viscous and inertial pressure losses between rw and r, each as a share of the corresponding losses between rw and re, —Inflow at high rates: non-Darcy flow
Δρ,g
density difference between water and oil at reservoir conditions, and the acceleration due to gravity, kg/m³, m/s²Water and gas coning
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≈1), Pa·s, —Inflow at high rates: non-Darcy flow
μˉ,Zˉ
gas viscosity and gas deviation factor at the average pressure (pe+pw)/2, Pa·s, —Steady-state radial flow
μ(p),Z(p)
gas viscosity and gas deviation factor as functions of pressure at reservoir temperature, Pa·s, —Steady-state radial flow
gas density at standard conditions: 1.204γg at 20 °C (γ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