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

Equations: Facilities and processing of production

The equations of the “Facilities and processing of production” lessons in “Oil and gas production”, with their conditions and units, and what each symbol means in them.

Oil rate under the oil and liquid limits
qo=min⁡ ⁣[ Q, L (1−WC) ],qw=qo WC1−WCq_o = \min\!\left[\,Q,\ L\,(1 - \text{WC})\,\right],\qquad q_w = q_o\,\frac{\text{WC}}{1 - \text{WC}}

any consistent units; volumes at standard conditions

Holds when: The wells can give more liquid than the facilities can take, so the rate is limited by the facilities, not by the reservoir. The water cut depends only on the share of the reserves already produced, not on the rate of production. This is a simplification: in practice fast production can make the waterflood sweep (the part of the reservoir reached by the water) worse. In the late years the economic limit also restricts production: when the revenue no longer covers the costs, the field is shut down.

Field development concepts: surface facilities →
Fuel gas consumption of a power plant
Vfuel=Nη QLHVV_{\text{fuel}} = \frac{N}{\eta\,Q_{\text{LHV}}}

SI; to get thousand m³/d, multiply the result by 86.4

Holds when: The efficiency is constant. In gas turbines it falls at part load and in hot weather, and then the fuel consumption is higher. The heating value of associated gas depends on its composition: the more ethane, propane and butanes, the higher it is.

Field development concepts: surface facilities →
Plateau rate as a share of the reserves
d=365 QNr⋅100 %d = \frac{365\,Q}{N_r}\cdot 100\ \%

volumes at standard conditions

Holds when: The measure is handy for comparing fields of different size. A sensible share depends on the reservoir, the development method and the costs, and it is not a norm.

Field development concepts: surface facilities →
Time in which the water cut grows from upu_p to uu
t−tp=2a[(up−u)+ln⁡1−up1−u],a=365 LNr (1−Rbt)t - t_p = \frac{2}{a}\left[(u_p - u) + \ln\frac{1 - u_p}{1 - u}\right],\qquad a = \frac{365\,L}{N_r\,(1 - R_{\text{bt}})}

time in years; LL in m³/d, NrN_r in m³

Holds when: The water cut follows this lesson's curve WC=(R−Rbt)/(1−Rbt)\text{WC} = \sqrt{(R - R_{\text{bt}})/(1 - R_{\text{bt}})}. The facilities are full of liquid all the time.

Field development concepts: surface facilities →
Allowed gas velocity (Souders–Brown equation)
vmax⁡=KρL−ρgρgv_{\max} = K\sqrt{\frac{\rho_{\text{L}} - \rho_{\text{g}}}{\rho_{\text{g}}}}

any consistent units: KK and vmax⁡v_{\max} in the same units of velocity, the densities in the same units

Holds when: The form follows from the balance of the weight and the drag of a droplet (the derivation is in the advanced part of the lesson). Everything that depends on the droplet size and on the way the gas flows around it is collected in KK. Rough guides for an estimate: for a vertical vessel with a mesh mist extractor KK is about 0.05–0.1 m/s, and for a horizontal one about 0.12–0.15 m/s (0.40–0.50 feet per second). The manufacturers give exact values for their internal parts. At high pressure, with foaming oil and with a heavy liquid load on the mist extractor, we reduce KK. Without a mist extractor the allowed velocity is clearly lower: everything that has not settled leaves with the gas.

Oil treating: separators, dehydration, desalting and stabilisation →
Retention time of the liquid
t=Vqt = \frac{V}{q}

any consistent units

Holds when: The flow goes through the whole volume evenly. Dead zones (places where the liquid hardly moves) and short paths of the flow (where part of the liquid runs straight to the outlet) reduce the real time, which is why the inlet and distribution devices are so important. Guides from industry practice: in two-phase separators from 1 to 3–4 min; in three-phase separators 3–10 min for light oil and 10–30 min or more for heavy and cold oil. For a particular oil, a laboratory settling test of a sample gives a better value.

Oil treating: separators, dehydration, desalting and stabilisation →
Share of the cross-section of a horizontal vessel that the liquid occupies
ALA=θ−sin⁡θ2π,θ=2arccos⁡ ⁣(1−2hD)\frac{A_{\text{L}}}{A} = \frac{\theta - \sin\theta}{2\pi},\qquad \theta = 2\arccos\!\left(1 - \frac{2h}{D}\right)

geometry; angles in radians

Holds when: The cylindrical part of the vessel. The volume of the heads (the rounded ends of the vessel) is counted separately or included in the margin.

Oil treating: separators, dehydration, desalting and stabilisation →
Settling or rising velocity of a droplet (Stokes' law)
v=g d2 Δρ18 μcv = \frac{g\,d^2\,\Delta\rho}{18\,\mu_c}

SI

Holds when: Slow (creeping) flow around the droplet: the Reynolds number of the droplet, Re=ρcvd/μc\mathrm{Re} = \rho_c v d/\mu_c (a number that shows how turbulent the flow around it is), is less than about 1. At larger Re\mathrm{Re} the drag grows faster, and the formula overestimates the velocity. The droplet is a rigid sphere, and the droplets are far from each other (at a high concentration they hinder each other and settle more slowly). The flow in the vessel is calm. Turbulence and convection from uneven heating hinder settling.

Oil treating: separators, dehydration, desalting and stabilisation →
Content of chloride salts in oil
S=φ CwS = \varphi\,C_{\text{w}}

mg/dm³; with a small fraction of water, the volume of the mixture is almost equal to the volume of the oil

Holds when: The salts are only dissolved in the water. Salt crystals and salts on solid particles are not counted.

Oil treating: separators, dehydration, desalting and stabilisation →
Wash water for desalting
w=1E(φ1CwCallow−φ1),Callow=Sreqφ2w = \frac{1}{E}\left(\frac{\varphi_1 C_{\text{w}}}{C_{\text{allow}}} - \varphi_1\right),\qquad C_{\text{allow}} = \frac{S_{\text{req}}}{\varphi_2}

fractions and mg/dm³

Holds when: The wash water is fresh, and the droplets that have mixed merge and average their salinity. The residual water after washing has the salinity of the mixture. In practice, small droplets of produced water that have not mixed are harder to remove from the oil than the mixed ones. So a model with a single fraction EE is a simplification that underestimates the wash water needed. EE is chosen from the performance of a particular unit: the mixing valve, the pressure drop across it, the number of washing stages.

Oil treating: separators, dehydration, desalting and stabilisation →
Net mass of oil
mnet=mgross(1−Ww+Wcs+Wimp100),Wcs=0.1 φcsρom_{\text{net}} = m_{\text{gross}}\left(1 - \frac{W_{\text{w}} + W_{\text{cs}} + W_{\text{imp}}}{100}\right),\qquad W_{\text{cs}} = \frac{0.1\,\varphi_{\text{cs}}}{\rho_{\text{o}}}

t, %, mg/dm³, kg/m³

Holds when: The gross mass is measured by a mass meter, or worked out from the volume and the density at the same conditions. The factor 0.1 converts units: mg/dm³ divided by kg/m³ (that is, by g/dm³) gives thousandths of a fraction, which is tenths of a per cent.

Oil treating: separators, dehydration, desalting and stabilisation →
Settling velocity of a droplet for any flow around it
πd36 (ρL−ρg) g=CD πd24⋅ρgvt22⇒vt=4 g d3 CD⋅ρL−ρgρg\frac{\pi d^3}{6}\,(\rho_{\text{L}} - \rho_{\text{g}})\,g = C_D\,\frac{\pi d^2}{4}\cdot\frac{\rho_{\text{g}} v_t^2}{2}\quad\Rightarrow\quad v_t = \sqrt{\frac{4\,g\,d}{3\,C_D}}\cdot\sqrt{\frac{\rho_{\text{L}} - \rho_{\text{g}}}{\rho_{\text{g}}}}

SI

Holds when: A single spherical droplet; the gas is at rest on average. With creeping flow CD=24/ReC_D = 24/\mathrm{Re}, and the formula turns into Stokes' law.

Oil treating: separators, dehydration, desalting and stabilisation →
Water content of natural gas saturated with water (Bukacek correlation)
yw=ps(T)p+β(T),β=10 2.0179−3083.871.8 T−0.07,W=748.9 ywy_w = \frac{p_s(T)}{p} + \beta(T),\qquad \beta = 10^{\,2.0179 - \frac{3083.87}{1.8\,T - 0.07}},\qquad W = 748.9\,y_w

absolute pressure in MPa (in the ratio ps/pp_s/p any identical units), temperature in kelvins

Holds when: Sweet natural gas, without noticeable hydrogen sulfide and carbon dioxide (they increase the water content), in equilibrium with fresh water (salty water gives a little less). The tested range of the correlation is from 15.5 to 238 °C and from 0.1 to 70 MPa. Below 15.5 °C, and still more in the hydrate region, the result is an estimate. We take ps(T)p_s(T) from an accurate equation for water (in the models of the lesson, the international equation for the vapour pressure of water). Any error in it goes straight into the result. 748.9 g/m³ is the mass of a kilomole of water (18.015 kg) divided by the volume of a kilomole of gas at 20 °C and 0.101325 MPa (24.055 m³).

Gas processing and compression →
Lowering of the hydrate formation temperature by an inhibitor (Hammerschmidt formula)
ΔT=K WM (100−W)\Delta T = \frac{K\,W}{M\,(100 - W)}

ΔT\Delta T in °C, WW in per cent by mass

Holds when: A dilute (weak) solution: for methanol up to about 20–25% by mass. At higher concentrations the formula overestimates ΔT\Delta T, that is, it gives too little inhibitor (the advanced part of the lesson has a logarithmic form for methanol). For monoethylene glycol, reference books give its own constant, a few per cent smaller. With the same constant, about twice as much MEG by mass is needed as methanol, because its molar mass is twice as large. The hydrate formation temperature without inhibitor must be known in advance, from the gas composition, laboratory tests or an estimate from the specific gravity.

Gas processing and compression →
Partial pressure of a gas component
pi=yi pp_i = y_i\,p

any identical units of pressure

Holds when: An ideal mixture of gases. At high pressures it is more accurate to use fugacity (a corrected pressure that accounts for non-ideal behaviour), but the standards for materials use this simple measure.

Gas processing and compression →
Mole fraction of a component from its mass concentration
y=cM nsc,nsc=pscR Tsc=41.57 mol/m3y = \frac{c}{M\,n_{\text{sc}}},\qquad n_{\text{sc}} = \frac{p_{\text{sc}}}{R\,T_{\text{sc}}} = 41.57\ \text{mol/m}^3

g/m³, g/mol; standard conditions 20 °C and 101.325 kPa

Holds when: The gas at standard conditions is ideal.

Gas processing and compression →
Compression ratio of a stage when the ratios are equal
r=(poutpin)1/nr = \left(\frac{p_{\text{out}}}{p_{\text{in}}}\right)^{1/n}

absolute pressures in identical units

Holds when: The gas is cooled after each stage to the same temperature. The pressure losses in coolers and separators are small (in practice the ratios are raised slightly to cover them).

Gas processing and compression →
Discharge temperature and specific work of a stage
T2=T1[1+r(k−1)/k−1η],w=kk−1⋅Z R T1M⋅r(k−1)/k−1ηT_2 = T_1\left[1 + \frac{r^{(k-1)/k} - 1}{\eta}\right],\qquad w = \frac{k}{k-1}\cdot\frac{Z\,R\,T_1}{M}\cdot\frac{r^{(k-1)/k} - 1}{\eta}

SI; the power of a stage is ww multiplied by the mass flow rate of the gas

Holds when: An ideal gas with a constant kk and the correction ZZ at suction. At high pressures, and for rich gas, it is more accurate to calculate with an equation of state. The losses are covered by one efficiency, referred to adiabatic compression. Mechanical losses and losses in the driver are added separately. The gas is single-phase: the liquid is separated before the stage.

Gas processing and compression →
Lowering of the hydrate formation temperature (the logarithmic form of Nielsen and Bucklin)
ΔT=−72.0 ln⁡xw,xw=11+18.015 WM (100−W)\Delta T = -72.0\,\ln x_w,\qquad x_w = \frac{1}{1 + \dfrac{18.015\,W}{M\,(100 - W)}}

°C (in °F the constant is 129.6)

Holds when: It was proposed for methanol. Reference books recommend it for concentrations from about 25 to 50% by mass, where the Hammerschmidt formula is already inaccurate. The solution is treated as ideal. For glycols, and for mixed inhibitors with salts, a calculation with an equation of state and a water activity model is more accurate.

Gas processing and compression →
Rising velocity of an oil droplet in water (Stokes' law)
vt=g d2 (ρw−ρo)18 μwv_t = \frac{g\,d^2\,(\rho_w - \rho_o)}{18\,\mu_w}

SI

Holds when: The droplet is a sphere and moves like a solid ball. For oil, which is more viscous than water, the circulation inside the droplet speeds up the rise by only a few per cent. The flow around the droplet is laminar (smooth): the Reynolds number of the droplet, ρwvtd/μw\rho_w v_t d/\mu_w, is less than about 1. For oil droplets in water this holds up to sizes of the order of 150 µm. The droplets do not hinder each other and do not merge, and the water is at rest. In a vessel, flow and mixing hinder the rise, while the merging of droplets helps it.

Produced water: composition, treatment and disposal →
The smallest droplet that an ideal settling tank catches completely
dmin⁡=18 μw qg (ρw−ρo) Ad_{\min} = \sqrt{\frac{18\,\mu_w\,q}{g\,(\rho_w - \rho_o)\,A}}

SI

Holds when: An ideal settling tank: the water moves evenly over the whole area, the droplets enter evenly over the height, and they neither break up nor merge. With a horizontal flow, droplets smaller than dmin⁡d_{\min} are caught in part: those that entered close enough to the surface. Their share is (d/dmin⁡)2(d/d_{\min})^2. In a vertical tank, where the water moves down against the droplets, all the small droplets leave with the water (the derivation is in the section for advanced study). In a real vessel part of the flow moves faster than the average (“short-circuiting”), and thermal currents and eddies interfere, so it catches less than the ideal one. The calculation gives the best case: the real smallest droplet that is caught is larger than the calculated one.

Produced water: composition, treatment and disposal →
Partial pressure of a gas component
pi=yi pp_i = y_i\,p

any units of pressure, the same for pip_i and pp

Holds when: The gas is treated as a mixture of ideal gases (Dalton's law). At high pressures the standards also allow more accurate methods. The pressure is absolute: the gauge pressure from a pressure gauge gives too low a value, especially at low pressures (Systems of units: SI, oilfield and metric oilfield). For the water, we take the gas that is in equilibrium with it: in the separator or in the vapour space of the tank.

Produced water: composition, treatment and disposal →
Water from the source at a given voidage replacement
qsrc=VRR Bo qo+Bw qwBw−qwq_{\text{src}} = \mathrm{VRR}\,\frac{B_o\,q_o + B_w\,q_w}{B_w} - q_w

rates at standard conditions, m³/d

Holds when: The pressure is above the bubble point pressure (the pressure at which gas starts to come out of the oil), and there is no free gas in the reservoir. Otherwise the gas is also counted in the voidage (Why water is injected: pressure maintenance and oil displacement). All the produced water is returned to the reservoir. Losses in treatment (water that left with the oil, the reject of hydrocyclones, the washing of filters) increase the water from the source a little. The produced water is compatible with the reservoir into which it is injected.

Produced water: composition, treatment and disposal →
Pore throat diameter by the model of a bundle of straight capillaries
dth≈28kϕd_{\text{th}} \approx 2\sqrt{\frac{8k}{\phi}}

SI

Holds when: The rock is represented by a bundle of identical straight tubes. With tortuosity τ\tau (the winding of the flow paths) the radius is τ\tau times larger, so the formula gives a lower estimate, and the allowed particle sizes derived from it have a margin. This is an order-of-magnitude estimate. Real throats have a wide distribution of sizes, and the design is based on the large ones, which carry the main flow.

Water injection systems →
Equilibrium oxygen content of water (Henry's law)
C=Csat pO2pO2airC = C_{\text{sat}}\,\frac{p_{\mathrm{O_2}}}{p_{\mathrm{O_2}}^{\text{air}}}

concentrations in identical units, pressures in identical units

Holds when: Equilibrium: the water is in contact with the gas long enough. In a column it is approached gradually, stage by stage, so in practice more oxygen is left. The temperature and salinity are the same as for CsatC_{\text{sat}}. With heating, the solubility of oxygen falls.

Water injection systems →
Power lost in throttling
Nthr=∑iqi (phdr−pwh,i)86.4N_{\text{thr}} = \sum_i \frac{q_i\,(p_{\text{hdr}} - p_{\text{wh},i})}{86.4}

metric oilfield units: m³/d, MPa, kW

Holds when: 86.486.4 is the same constant as in the formula for the hydraulic power of a pump: the seconds in a day divided by a thousand. The losses in the water lines from the station to the block are the same for all wells and are not included here. The electric power is 1/η1/\eta times the hydraulic power, where η\eta is the efficiency of the pump unit.

Water injection systems →
Injectivity index
I=qpw−pˉI = \frac{q}{p_w - \bar p}

metric oilfield units

Holds when: The injection is close to steady state. Right after a change of the rate, the pressure has not settled yet. The reservoir pressure in the injection area is known, from a falloff test (pressure falloff after shut-in) or from neighbouring wells. The fracture has not opened. After it opens, the rate grows with almost no rise in pressure, and the injectivity stops being a constant.

Water injection systems →
Bottomhole pressure of an injection well
pw=pwh+ρgH−Δpfr,Δpfr=λ Hd ρv22p_w = p_{\text{wh}} + \rho g H - \Delta p_{\text{fr}},\qquad \Delta p_{\text{fr}} = \lambda\,\frac{H}{d}\,\frac{\rho v^2}{2}

SI; to get MPa, divide by 10610^6

Holds when: A single-phase flow of water. The viscosity and density are at the mean temperature in the tubing. Friction grows roughly as the square of the rate, so at a low rate it can be neglected, and at a high rate it cannot. Injection down the annulus (the space between the tubing and the casing) or into a deviated well needs its own lengths and diameters.

Water injection systems →
Slope of the Hall plot and injectivity
mH=Δ ⁣[∑(pw−pˉ) Δt]ΔW=1Im_H = \frac{\Delta\!\left[\sum (p_w - \bar p)\,\Delta t\right]}{\Delta W} = \frac{1}{I}

metric oilfield units

Holds when: The injection in each period is close to steady state, and the reservoir pressure in the injection area hardly changes. Its rise also increases the slope. We calculate the bottomhole pressure from the wellhead pressure, counting the hydrostatic pressure and the friction at the actual rate. The slope over a segment is an average. If the injectivity changed within the segment, the slope is the average value of 1/I1/I, weighted by the injection.

Water injection systems →

Symbols

SymbolMeaning and context
ΔT\Delta T
pˉ\bar p
ΔW\Delta W
Δt\Delta t
aa
AA
ALA_{\text{L}}
Bo, BwB_o,\ B_w
cc
CC
CallowC_{\text{allow}}
CDC_D
CsatC_{\text{sat}}
  • concentration in water saturated with air at atmospheric pressure, at the same temperature and salinity, mg/LWater injection systems
CwC_{\text{w}}
dd
DD
dmin⁡d_{\min}
dthd_{\text{th}}
d, vd,\ v
  • inner diameter of the tubing and the mean velocity of the water in it, v=4q/(πd2)v = 4q/(\pi d^2), m, m/sWater injection systems
EE
gg
hh
HH
  • vertical depth of the middle of the injection interval (for friction, the length of the tubing along the hole; in a vertical well they coincide), mWater injection systems
II
kk
  • adiabatic exponent of the gas, cp/cvc_p/c_v (the ratio of the heat capacities at constant pressure and at constant volume): about 1.2–1.3 for natural gas, —Gas processing and compression
  • permeability of the reservoir (1 mD = 0.987⋅10−150.987\cdot 10^{-15} m²), m²Water injection systems
KK
LL
ln⁡\ln
MM
mHm_H
mnet, mgrossm_{\text{net}},\ m_{\text{gross}}
nn
NN
NrN_r
nscn_{\text{sc}}
NthrN_{\text{thr}}
pp
phdrp_{\text{hdr}}
pip_i
pO2p_{\mathrm{O_2}}
  • partial pressure of oxygen in the gas over the water (the share of oxygen times the absolute pressure), kPaWater injection systems
pO2airp_{\mathrm{O_2}}^{\text{air}}
  • partial pressure of oxygen in air: 0.2095⋅101.325=21.230.2095\cdot 101.325 = 21.23 kPa, kPaWater injection systems
pout, pinp_{\text{out}},\ p_{\text{in}}
ps(T)p_s(T)
psc, Tscp_{\text{sc}},\ T_{\text{sc}}
pwp_w
  • bottomhole pressure during injection (for producing wells the bottomhole pressure is written pwfp_{wf}), MPaWater injection systems
pwhp_{\text{wh}}
pwh,ip_{\text{wh},i}
p, Tp,\ T
qq
QQ
qiq_i
QLHVQ_{\text{LHV}}
qo, qwq_o,\ q_w
QscQ_{\text{sc}}
qsrcq_{\text{src}}
rr
RR
SS
SreqS_{\text{req}}
tt
TT
T1, T2T_1,\ T_2
tpt_p
uu
upu_p
vv
VV
VfuelV_{\text{fuel}}
vmax⁡v_{\max}
vtv_t
VRR\mathrm{VRR}
ww
WW
  • water content: the mass of water in one cubic metre of gas at 20 °C and 0.101325 MPa, g/m³Gas processing and compression
  • concentration of the inhibitor in the water phase (water plus inhibitor). The letter is the one used in the Hammerschmidt formula; it is not the water content of the gas WW from the previous section, % by massGas processing and compression
WcsW_{\text{cs}}
Ww, WimpW_{\text{w}},\ W_{\text{imp}}
W, MW,\ M
WC\text{WC}
xwx_w
yy
yiy_i
ywy_w
ZZ
β(T)\beta(T)
Δρ\Delta\rho
η\eta
θ\theta
λ\lambda
μc\mu_c
μw\mu_w
μw, ρw, ρo, g\mu_w,\ \rho_w,\ \rho_o,\ g
ρ\rho
ρL, ρg\rho_{\text{L}},\ \rho_{\text{g}}
ρo\rho_{\text{o}}
ρw, ρo\rho_w,\ \rho_o
φ\varphi
ϕ\phi
φ1, φ2\varphi_1,\ \varphi_2
φcs\varphi_{\text{cs}}