Production Engineering & Artificial Lift
Getting fluids from the reservoir to the surface. Learn well completions (casing, tubing, packers, perforations, sand control), the inflow performance relationship (productivity index, Vogel), vertical lift performance and multiphase flow, nodal analysis, artificial lift (rod pumps, ESP, gas lift, PCP, jet pumps), stimulation (acidising, hydraulic fracturing), surface facilities (separators, water treatment) and flow assurance (wax, hydrates, scale) — with Python for Vogel IPR and a nodal operating point.
Production engineers maximise the rate each well can sustainably deliver. Flow is controlled by the reservoir's ability to deliver fluid (inflow) and the well's ability to lift it to the surface (outflow); their intersection is the operating point.
As fields mature and pressure falls, wells need artificial lift, stimulation and careful management of water, gas, sand and deposits.
An analogy
A well is like drinking through a straw from a pressurised bottle: the bottle pressure pushes liquid into the straw (inflow), while gravity and friction in the straw resist lifting it (outflow). Artificial lift is like sucking harder on the straw.
Completions
After drilling, the well is completed: production casing and cement, perforations (shaped charges) through casing into the reservoir, production tubing with a packer, safety valves (SCSSV), wellhead and Christmas tree. Sand control (gravel packs, screens) is needed in weak sandstones; smart completions have downhole valves and gauges.
Inflow and outflow
Inflow: straight-line PI above bubble point; Vogel (or Fetkovich) curve below. Skin damage (drilling mud invasion, perforation damage) reduces inflow. Outflow: vertical lift performance from multiphase flow correlations (Beggs–Brill, Hagedorn–Brown) — pressure loss from gravity (dominant), friction and acceleration. Nodal analysis finds the rate and tests sensitivities (tubing size, wellhead pressure, skin, water cut).
Artificial lift
Choose by rate, depth, GOR, sand and power availability:
- Sucker-rod (beam) pumps — onshore, low rates, simple
- Electrical submersible pumps (ESP) — high rates (up to 30 000+ bbl/d)
- Gas lift — inject gas to lighten the fluid column; tolerant of sand and gas
- Progressive cavity pumps — viscous/sandy oil
- Jet and hydraulic pumps; plunger lift for gas wells
Stimulation, facilities and flow assurance
Matrix acidising removes near-wellbore damage (HCl for carbonates, mud acid for sandstones); hydraulic fracturing creates propped fractures (essential in shale and tight rock). Surface: three-phase separators, heater treaters, gas compression and dehydration, produced-water treatment and reinjection. Flow assurance manages hydrates (inhibitors, insulation), wax, asphaltenes, scale and corrosion; digital oilfield tools monitor wells in real time.
Formulas
- Productivity index:
q = J·(p̄_r − p_wf) (undersaturated) - Vogel IPR (saturated):
q/q_max = 1 − 0.2(p_wf/p̄_r) − 0.8(p_wf/p̄_r)² - Outflow (VLP):
p_wf = p_wh + ΔP_hydrostatic + ΔP_friction - Nodal analysis:
operating point where inflow p_wf = outflow p_wf - Hydraulic horsepower (ESP):
HHP = q(bbl/d) × ΔP(psi) × 1.7×10⁻⁵
Worked example: Productivity index
Given: q 800 STB/d, p̄_r 2500, p_wf 1700 psi
Find: J.
- 800/800
Answer: 1.0 STB/d/psi
Worked example: Vogel AOF
Given: 1200 STB/d at 2000 psi, p̄_r 3000
Find: q_max.
- 1200/(1 − 0.1333 − 0.3556)
Answer: 2348 STB/d
Worked example: ESP hydraulic power
Given: 2000 bbl/d, 2500 psi lift
Find: HHP.
- 2000 × 2500 × 1.7×10⁻⁵
Answer: 85 hp (motor ~ 150 hp at 55% efficiency)
Vogel IPR and nodal operating point
# Inflow performance: productivity index and Vogel IPR, then operating point with a tubing curve
Pr, Pb = 3000.0, 3000.0 # reservoir and bubble-point pressure (psi) -> saturated reservoir
q_test, pwf_test = 1200.0, 2000.0
qmax = q_test / (1 - 0.2 * (pwf_test / Pr) - 0.8 * (pwf_test / Pr) ** 2)
print(f"Vogel AOF qmax = {qmax:.0f} STB/d")
for pwf in (2500, 1500, 500):
q = qmax * (1 - 0.2 * pwf / Pr - 0.8 * (pwf / Pr) ** 2)
print(f"pwf {pwf} psi -> q {q:.0f} STB/d")
# Nodal: tubing requires pwf = 800 + 0.9*q (simple linear outflow)
lo, hi = 0.0, qmax
for _ in range(60):
mid = (lo + hi) / 2
pwf_in = Pr * (-0.2 + (0.04 + 3.2 * (1 - mid / qmax)) ** 0.5) / 1.6
if pwf_in > 800 + 0.9 * mid: lo = mid
else: hi = mid
print(f"operating point ~ {lo:.0f} STB/d at pwf {800+0.9*lo:.0f} psi")
# Vogel AOF qmax = 2348 STB/d
# pwf 2500 psi -> q 652 STB/d
# pwf 1500 psi -> q 1643 STB/d
# pwf 500 psi -> q 2217 STB/d
# operating point ~ 1262 STB/d at pwf 1936 psiOne well test (1200 STB/d at 2000 psi) calibrates the Vogel curve. Drawdown gives diminishing returns below the bubble point (gas blocks the pores). With the tubing outflow curve, the well naturally flows at ~1260 STB/d; lowering p_wf with artificial lift could move it toward 2000+ STB/d. (Outflow line is a simplified illustration; real VLP uses multiphase correlations.)
Common mistakes
- Wrong: Using a straight-line PI below the bubble point. — Right: Use Vogel/Fetkovich.
- Wrong: Oversizing tubing in low-rate gas wells. — Right: Liquid loading — check critical velocity.
- Wrong: Ignoring hydrate risk at shutdown. — Right: Inhibit and insulate subsea lines.
Applications
- Well optimisation: Choke, tubing and lift tuning.
- Mature fields: Artificial lift and water management.
- Shale wells: Multi-stage hydraulic fracturing.
- Subsea tie-backs: Flow assurance over long distances.
Key takeaways
- q = J(p̄_r − p_wf) or Vogel below P_b.
- Operating point = inflow ∩ outflow.
- Artificial lift lowers p_wf.
- Manage sand, water, hydrates, wax, scale.
Going further
- Multiphase flow modelling: Mechanistic models and OLGA transient simulation.
- Hydraulic fracture design: Fracture geometry and proppant transport.
- Production optimisation: Integrated asset models and gas-lift allocation.