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

Equations: Gas lift

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

Bottomhole pressure of a gas-lift well
pwf=pwh+Δpabove ⁣(q, GLRres+qgq)+Δpbelow(q)p_{wf} = p_{wh} + \Delta p_{\text{above}}\!\left(q,\ \text{GLR}_{\text{res}} + \frac{q_{\text{g}}}{q}\right) + \Delta p_{\text{below}}(q)

metric oilfield units; pressures are absolute, gas volumes are at standard conditions

Holds when: Steady state; all the injected gas enters the tubing through the operating valve. Gas from the reservoir is added to the injected gas wherever it is free. Below the injection point the pressure is usually above the bubble point (the pressure at which gas starts to come out of the oil), so the gas is dissolved in the oil. The pressure differences are calculated with correlations or models of multiphase flow (Multiphase flow in tubing). The examples in this lesson use a simple teaching drift-flux model (a model in which the gas drifts upward through the liquid at its own speed).

Gas lift: principles and choice of operating conditions →
Gas pressure in the annulus at a given depth
pg(H)=pg,wh exp⁡ ⁣(MgHZˉRTˉ)p_{\text{g}}(H) = p_{\text{g,wh}}\,\exp\!\left(\frac{M g H}{\bar Z R \bar T}\right)

SI; pressures are absolute

Holds when: The gas in the annulus is almost still. Its velocity is low, so friction is negligible. The temperature and ZZ are taken as averages over the column. ZZ depends on the average pressure, which is itself being found, so the calculation is repeated two or three times. The same gas fills the whole column. The formula is the same as for any column of gas (column of gas).

Gas lift: principles and choice of operating conditions →
Condition for injecting gas at a given depth
pg(Hinj)−Δpvalve≥pt(Hinj)p_{\text{g}}(H_{\text{inj}}) - \Delta p_{\text{valve}} \ge p_{\text{t}}(H_{\text{inj}})

any consistent units; the pressures are all absolute or all gauge

Holds when: Steady state; gas enters the tubing only through the operating valve. The pressure drop across the valve depends on the gas rate and on the size of the valve port (its opening). How it is calculated is explained in the lesson Unloading a gas lift well: valves and mandrels.

Gas lift: principles and choice of operating conditions →
Condition for the technical optimum
(∂Δpweight∂qg)q+(∂Δpfr∂qg)q=0\left(\frac{\partial \Delta p_{\text{weight}}}{\partial q_{\text{g}}}\right)_{q} + \left(\frac{\partial \Delta p_{\text{fr}}}{\partial q_{\text{g}}}\right)_{q} = 0

any consistent units

Holds when: The acceleration of the mixture is not taken into account. The column below the valve does not depend on the gas rate. The flow rate is at its maximum when extra gas stops lowering the pressure that is needed to lift this flow rate (the derivation is in the advanced section).

Gas lift: principles and choice of operating conditions →
Economic optimum of the gas rate
 dqo dqg=cgvo\frac{\dd q_o}{\dd q_{\text{g}}} = \frac{c_{\text{g}}}{v_o}

any consistent units: the ratio of the costs is in the same units as the slope of the curve

Holds when: The oil price and the cost of gas do not depend on the volumes. The gas supply of the system is not limited. If it is limited, the gas is shared between wells (next section). The limits of the surface facilities on water and gas are not taken into account. They move the optimum to the left.

Gas lift: principles and choice of operating conditions →
The best sharing of gas between wells
 dqo,1 dqg,1= dqo,2 dqg,2=⋯=β,∑iqg,i=Qg\frac{\dd q_{o,1}}{\dd q_{\text{g},1}} = \frac{\dd q_{o,2}}{\dd q_{\text{g},2}} = \dots = \beta, \qquad \sum_i q_{\text{g},i} = Q_{\text{g}}

any consistent units

Holds when: The curves of the wells bulge upward (they are concave): each next portion of gas gives less than the previous one. The wells do not affect each other: the pressure in the gathering system does not change with the way the gas is shared. A well whose marginal return is below β\beta even at a small gas rate gets no gas. If a well needs a minimum gas rate to start, the simple rule breaks down. We then compare the options with and without that well (below).

Gas lift: principles and choice of operating conditions →
Power of ideal isothermal compression of gas
Niso=psc qg ln⁡p2p1N_{\text{iso}} = p_{\text{sc}}\,q_{\text{g}}\,\ln\frac{p_2}{p_1}

SI; pressures are absolute

Holds when: The gas is ideal and the compression is isothermal at 20 °C (an infinite number of stages with cooling), with no losses. This is a lower limit for real compressors. If the gas is compressed at a temperature TT, the power is higher by a factor of T/(293 K)T/(293\ \text{K}). The compressibility of a real gas (Z<1Z < 1) reduces it by 5–10%. At 40–50 °C these two corrections almost cancel each other. A real multistage compressor with its drive uses more. In the example below it is assumed to use 1/0.651/0.65 times as much.

Gas lift: principles and choice of operating conditions →
Slope of the performance curve of a gas-lift well
 dq dqg=− ∂poutflow/∂qg∂poutflow/∂q+1/J\frac{\dd q}{\dd q_{\text{g}}} = -\,\frac{\partial p_{\text{outflow}}/\partial q_{\text{g}}}{\partial p_{\text{outflow}}/\partial q + 1/J}

any consistent units

Holds when: The IPR is a straight line, and the operating point is stable, so the denominator is positive (stability of the operating point).

Gas lift: principles and choice of operating conditions →
Depth to which the kick-off pressure pushes out the kill fluid
pkoexp⁡ ⁣(0.03418 γg HZˉ Tˉ)=pwh+ρkill g Hp_{\text{ko}}\exp\!\left(\frac{0.03418\,\gamma_{\text{g}}\,H}{\bar Z\,\bar T}\right) = p_{wh} + \rho_{\text{kill}}\,g\,H

metric oilfield units; pressures are absolute

Holds when: The left side is the pressure of the column of gas in the annulus, as in the lesson on gas lift. The right side is the pressure of the liquid in the tubing, which is full up to the wellhead. The balance is static. During displacement, flow friction adds to the pressure in the tubing, and the level stops a little higher still. There is a check valve (a one-way valve) at the bottom of the tubing, so the liquid does not go into the reservoir. Otherwise part of it would go into the reservoir, and the reservoir would first return that same liquid.

Unloading a gas lift well: valves and mandrels →
Gas pressure that opens the valve
pdomeAb=pg (Ab−Ap)+ptAp⇒pg,open=pdome−R pt1−R,R=ApAbp_{\text{dome}}A_{\text{b}} = p_{\text{g}}\,(A_{\text{b}} - A_{\text{p}}) + p_{\text{t}}A_{\text{p}}\quad\Rightarrow\quad p_{\text{g,open}} = \frac{p_{\text{dome}} - R\,p_{\text{t}}}{1 - R},\qquad R = \frac{A_{\text{p}}}{A_{\text{b}}}

any consistent units; the pressures are all absolute or all gauge, and the equation does not change, because the areas on each side add up to AbA_{\text{b}}

Holds when: The valve has no spring. The elasticity of the bellows itself is ignored. The pressures are uniform over the areas on which they act.

Unloading a gas lift well: valves and mandrels →
Closing pressure, and the difference between the opening and closing pressures
pg,close=pdome,Δpspread=pg,open−pg,close=kt (pdome−pt),kt=R1−Rp_{\text{g,close}} = p_{\text{dome}},\qquad \Delta p_{\text{spread}} = p_{\text{g,open}} - p_{\text{g,close}} = k_{\text{t}}\,(p_{\text{dome}} - p_{\text{t}}),\qquad k_{\text{t}} = \frac{R}{1 - R}

any consistent units

Holds when: When the ball has moved away from the seat, the gas pressure acts on the whole area of the bellows, and the valve stays open as long as pg>pdomep_{\text{g}} > p_{\text{dome}}. The valve has no spring. In reality the elasticity of the bellows shifts the closing pressure, as explained below.

Unloading a gas lift well: valves and mandrels →
Tubing pressure that opens a valve operated by the tubing pressure
pt,open=pdome−R pg1−Rp_{\text{t,open}} = \frac{p_{\text{dome}} - R\,p_{\text{g}}}{1 - R}

any consistent units; the pressures are all absolute or all gauge

Holds when: The valve has no spring. The tubing pressure acts on the bellows from outside, and the gas pressure acts on the seat.

Unloading a gas lift well: valves and mandrels →
Test-rack opening pressure
pTRO=pdomerack1−Rp_{\text{TRO}} = \frac{p_{\text{dome}}^{\text{rack}}}{1 - R}

gauge pressures, as shown by the gauges of the rack

Holds when: On the tubing side there is atmospheric pressure, that is, pt=0p_{\text{t}} = 0 on the gauge. The equation is a special case of the force balance.

Unloading a gas lift well: valves and mandrels →
Temperature correction for the nitrogen in the dome
CT=pdomerackpdome=11+0.00215 (T∘F−60)C_T = \frac{p_{\text{dome}}^{\text{rack}}}{p_{\text{dome}}} = \frac{1}{1 + 0.00215\,(T_{^{\circ}\text{F}} - 60)}

temperature in degrees Fahrenheit; pressures are gauge, as in the accepted field method

Holds when: It is a straight-line fit to measurements on nitrogen at the constant volume of the dome. The tables of factors for setting valves are made from it. It is accurate enough at a nitrogen pressure of about 4 MPa. At higher pressures real nitrogen raises its pressure more when heated, and the correction underestimates it (advanced section). If part of the dome is filled with silicone fluid, it also expands on heating and adds pressure.

Unloading a gas lift well: valves and mandrels →
Greatest depth of the next unloading valve
pt,des(Hn)+ρkill g (Hn+1−Hn)=pg(Hn+1; pg,wh,n+1)−Δpmargin,pg,wh,n=pko−(n−1) Δpstepp_{\text{t,des}}(H_n) + \rho_{\text{kill}}\,g\,(H_{n+1} - H_n) = p_{\text{g}}\big(H_{n+1};\ p_{\text{g,wh},n+1}\big) - \Delta p_{\text{margin}},\qquad p_{\text{g,wh},n} = p_{\text{ko}} - (n - 1)\,\Delta p_{\text{step}}

metric oilfield units; pressures are absolute

Holds when: The well is vertical. (In a deviated well, true vertical depths are converted into measured depths along the hole for the mandrels.) Below the working valve, the annulus and the tubing hold kill fluid, and the level in the tubing stays at this valve. The friction of the liquid while it is pushed out is not included, and the margin covers it. The tubing pressure at the working valve is taken from the design line, and enough gas passes through the valve for this mode.

Gas lift design →
Spacing to the next valve with constant gradients
Hn+1−Hn=pg,wh,n+1+ggHn−Δpmargin−pt,des(Hn)ρkill g−ggH_{n+1} - H_n = \frac{p_{\text{g,wh},n+1} + g_{\text{g}}H_n - \Delta p_{\text{margin}} - p_{\text{t,des}}(H_n)}{\rho_{\text{kill}}\,g - g_{\text{g}}}

metric oilfield units; pressures are absolute

Holds when: The gas gradient is constant. It is taken from the formula for a column of gas at the average pressure over the interval, so the error is small. For the first valve H0=0H_0 = 0 and pt,des(0)=pwhp_{\text{t,des}}(0) = p_{wh}.

Gas lift design →
The stock of gas in the annulus and the rate of change of the pressure in it
Vsc=Vann pˉpsc TscTˉ 1Zˉ, dpˉ dt≈pˉ qin−qoutVscV_{\text{sc}} = V_{\text{ann}}\,\frac{\bar p}{p_{\text{sc}}}\,\frac{T_{\text{sc}}}{\bar T}\,\frac{1}{\bar Z},\qquad \frac{\dd\bar p}{\dd t} \approx \bar p\,\frac{q_{\text{in}} - q_{\text{out}}}{V_{\text{sc}}}

any consistent units; pressures are absolute

Holds when: The temperature and Zˉ\bar Z hardly change during a swing. The pressure in the annulus changes almost equally over the whole height. The second equality follows from the first: with constant VannV_{\text{ann}}, Tˉ\bar T and Zˉ\bar Z, the pressure is proportional to the amount of gas.

Gas lift operation and optimisation →

Symbols

SymbolMeaning and context
Zˉ\bar Z
Tˉ\bar T
pˉ, Tˉ, Zˉ\bar p,\ \bar T,\ \bar Z
(⋅)q(\cdot)_{q}
0.034180.03418
Ab, ApA_{\text{b}},\ A_{\text{p}}
Δpabove\Delta p_{\text{above}}
Δpbelow\Delta p_{\text{below}}
cgc_{\text{g}}
CTC_T
Δpfr\Delta p_{\text{fr}}
gg
ggg_{\text{g}}
  • pressure gradient of the gas in the annulus, taken as constant, MPa/mGas lift design
GLRres\text{GLR}_{\text{res}}
HH
HinjH_{\text{inj}}
Hn, Hn+1H_n,\ H_{n+1}
  • true vertical depths of the working valve and of the next one, mGas lift design
JJ
ktk_{\text{t}}
MM
Δpmargin\Delta p_{\text{margin}}
  • margin: the amount by which the gas at the opening valve must exceed the pressure of the liquid, MPaGas lift design
NisoN_{\text{iso}}
∂poutflow/∂q\partial p_{\text{outflow}}/\partial q
∂poutflow/∂qg\partial p_{\text{outflow}}/\partial q_{\text{g}}
p1, p2p_1,\ p_2
pdomep_{\text{dome}}
pdome, pg, Rp_{\text{dome}},\ p_{\text{g}},\ R
pdomerackp_{\text{dome}}^{\text{rack}}
pgp_{\text{g}}
pg(Hinj)p_{\text{g}}(H_{\text{inj}})
pg(H; pg,wh)p_{\text{g}}(H;\ p_{\text{g,wh}})
pg(H)p_{\text{g}}(H)
pg,closep_{\text{g,close}}
pg,whp_{\text{g,wh}}
pg,wh,np_{\text{g,wh},n}
  • gas pressure at the annulus wellhead at which valve nn opens, MPaGas lift design
pg,wh,n+1, Δpmargin, pt,des(Hn)p_{\text{g,wh},n+1},\ \Delta p_{\text{margin}},\ p_{\text{t,des}}(H_n)
  • as in the equation for the greatest depth of the next valve, MPaGas lift design
pkop_{\text{ko}}
pscp_{\text{sc}}
psc, Tscp_{\text{sc}},\ T_{\text{sc}}
ptp_{\text{t}}
pt(Hinj)p_{\text{t}}(H_{\text{inj}})
pt,des(Hn)p_{\text{t,des}}(H_n)
  • tubing pressure at depth HnH_n on the design line; for the first valve, the tubing pressure at the wellhead pwhp_{wh} at H0=0H_0 = 0, MPaGas lift design
pt,openp_{\text{t,open}}
pTROp_{\text{TRO}}
pwfp_{wf}
pwhp_{wh}
qq
qgq_{\text{g}}
QgQ_{\text{g}}
qg,iq_{\text{g},i}
qin, qoutq_{\text{in}},\ q_{\text{out}}
qoq_o
qo,iq_{o,i}
RR
Δpspread\Delta p_{\text{spread}}
Δpstep\Delta p_{\text{step}}
  • pressure step: the amount by which each next valve opens at a lower pressure at the wellhead, MPaGas lift design
T∘FT_{^{\circ}\text{F}}
VannV_{\text{ann}}
vov_o
VscV_{\text{sc}}
Δpvalve\Delta p_{\text{valve}}
Δpweight\Delta p_{\text{weight}}
β\beta
γg, Zˉ, Tˉ\gamma_{\text{g}},\ \bar Z,\ \bar T
ρkill\rho_{\text{kill}}
ρkill g\rho_{\text{kill}}\,g