Course 01 · Well Productivity Fundamentals› Module 04 · Multiphase & Gas Well Inflow› Topic 4.1

Course 01 · Module 04 · Topic 4.1

Two-Phase Flow: Vogel's Method

When flowing pressures drop below the bubble point, liberated gas fundamentally changes inflow behaviour — producing a curved IPR that a straight-line PI fatally underestimates. Vogel's empirical reference curve is the industry's foundational tool for quantifying this effect.

In Module 2 you mastered the Darcy radial flow equation and the straight-line Productivity Index (PI). That model assumes single-phase liquid flow all the way from the drainage boundary to the wellbore — a valid assumption when the entire flowing pressure path remains above the bubble point pressure, p_b.

In mature fields, high-drawdown wells, and depleted reservoirs, this assumption routinely fails. As bottom-hole flowing pressure (BHFP) drops below p_b, dissolved gas evolves from solution. Free gas occupies pore space, competes with oil for flow paths, and reduces the effective permeability to oil. The result is a progressive non-linear reduction in oil production rate with increasing drawdown — the curved IPR that Vogel (1968) first quantified for solution-gas drive reservoirs.

Understanding Vogel's method is not merely academic. It directly determines how much production you can recover by lowering BHFP through artificial lift, and it sets the baseline against which stimulation improvements are measured.

▶

Lecture 4.1A: Why Single-Phase PI Breaks Down Below Bubble Point

18:30 · HD

Covers the physics of dissolved gas liberation, the resulting two-phase zone near the wellbore, effective permeability concepts, and why the Darcy PI equation progressively overestimates deliverability as drawdown increases below p_b. Includes animated wellbore pressure profile.

LEARNING OBJECTIVES

After completing Topic 4.1, you will be able to:

1. Explain why oil well IPR curves below the bubble point are non-linear and what physical mechanism causes the curvature.

2. State Vogel's reference equation and identify every term, including the significance of p_wf/p_R ratios.

3. Construct a complete Vogel IPR curve from a single well test data point (flow rate + BHFP).

4. Calculate q_max (AOF at p_wf = 0) and any intermediate flow rate at a specified BHFP.

5. Determine the theoretical productivity index J* and relate it to q_max via the Vogel relationship.

6. Identify the limitations of Vogel's method and the conditions under which it should not be applied without modification.

PREREQUISITE KNOWLEDGE

Required: Darcy radial flow equation (Topic 2.1), Productivity Index concept (Topic 2.2), bubble point pressure and solution GOR from PVT (Module 01). You must be comfortable plotting and interpreting IPR curves on Cartesian axes (q vs. p_wf) before proceeding.

PBL CONNECTION — KARAMA FIELD PROBLEM

In the Karama Field scenario, Well KA-07 is currently flowing above p_b at 4,200 psi BHFP with a reservoir pressure of 5,100 psi. The operator's 5-year depletion plan projects reservoir pressure declining to 3,800 psi — well below KA-07's bubble point of 4,500 psi. You will use Vogel's method in this topic to predict KA-07's future IPR and evaluate whether the planned ESP installation will sustain economic production rates. All simulator exercises use KA-07 base data.

Topic Scope

Vogel's method for solution-gas drive. Pure two-phase IPR below p_b with zero skin.

Connects To

Topic 4.2 extends this to composite IPR (above + below p_b). Topic 4.3 adds skin correction via Standing.

Estimated Time

~90 minutes total: 35 min reading, 20 min simulation, 20 min worked examples, 15 min quiz.

Section 1

The Physics of Two-Phase Inflow

Understanding why the IPR curves downward requires a clear picture of what happens to gas saturation in the near-wellbore region as BHFP falls below p_b.

▶

Lecture 4.1B: Near-Wellbore Two-Phase Zone — Visualised

14:45 · HD

Animated cross-section of the drainage area. Shows how the two-phase zone develops, expands, and intensifies as BHFP is drawn down. Includes animated relative permeability curves responding to changing gas saturation.

From Single-Phase to Two-Phase: The Bubble Point Threshold

At pressures above p_b, all gas remains dissolved in the oil. The reservoir rock transmits only liquid, and Darcy's law with a single effective permeability k_o accurately describes inflow. The PI is constant, and the IPR is a straight line.

When p_wf drops below p_b, gas nucleates and evolves from the oil in the immediate vicinity of the wellbore — precisely where the lowest pressure exists. This liberated gas must occupy pore space. Since total pore volume is fixed, gas saturation S_g increases and oil saturation S_o falls. The consequences cascade through relative permeability.

Relative Permeability: The Mechanism of IPR Curvature

Relative permeability k_ro is the ratio of effective oil permeability to absolute permeability. It is a strong function of oil saturation (or equivalently, gas saturation). As S_g increases near the wellbore, k_ro decreases and oil mobility k_ro/µ_o falls.

The Darcy radial flow equation shows that oil rate is directly proportional to effective oil permeability:

q_o = [0.00708 · k_o · h · (p_R - p_wf)] / [µ_oB_o · (ln(0.472r_e/r_w) + S')] where k_o = k · k_ro(S_w, S_g)

As BHFP falls further below p_b, more gas comes out of solution, S_g grows, k_ro drops, so even though the pressure drawdown (p_R - p_wf) is increasing, the rate doesn't increase proportionally. This is the physical origin of IPR curvature.

CRITICAL INSIGHT

The two-phase zone is most severe at the wellbore face, where pressure is lowest. This is also where flow converges most intensely. The damage is therefore weighted to the highest-flux region of the reservoir, making it disproportionately impactful on productivity compared to the same permeability change located further away.

Why a Straight-Line PI Overestimates Production

If an engineer incorrectly uses the straight-line PI determined at high BHFP (where flow was single-phase) to predict rates at low BHFP, they will overestimate oil production by a significant margin. The table below illustrates this for a typical field well.

p_wf (psi)	Drawdown (psi)	Linear PI Prediction (stb/d)	Actual Rate with k_ro Reduction (stb/d)	Error (%)
4,000 (above p_b)	500	350	350	0%
3,500 (≈ p_b)	1,000	700	650	+8%
2,500	2,000	1,400	1,050	+33%
1,500	3,000	2,100	1,290	+63%
500	4,000	2,800	1,470	+90%
0 (AOF)	4,500	3,150	1,750	+80%

KEY TAKEAWAY

Vogel (1968) found through reservoir simulation that, for solution-gas drive, the actual q_max at zero BHFP is typically only 55% of the maximum rate predicted by the straight-line PI extended to zero flowing pressure. Ignoring this leads to gross overestimation of what artificial lift can achieve.

Below p_b, as gas evolves from solution, the remaining oil becomes heavier and more viscous. The live oil viscosity µ_o increases because the GOR of the oil declines. Simultaneously, the formation volume factor B_o decreases (oil shrinks as dissolved gas leaves). Both effects act to reduce q_o further — reinforcing the curvature that Vogel's simulation studies captured.

Practical implication: PVT data below p_b shows µ_o rising and B_o falling with decreasing pressure. These must be accounted for in rigorous nodal analysis — Vogel's equation captures the net effect empirically without requiring the engineer to explicitly track each property change.

Darcy's radial flow equation tells us pressure drops logarithmically from drainage radius to wellbore. Most of the total pressure drop occurs in a small annular region very close to the wellbore. For example, with r_e = 1,320 ft and r_w = 0.33 ft, approximately 50% of the total ln(r_e/r_w) drop occurs in the first 5 feet from the wellbore.

This means the two-phase zone (where p < p_b) is concentrated in the high-drawdown region near the well — the region that contributes most to total pressure drop. The skin-equivalent damage from two-phase flow is therefore amplified far more than its spatial extent would suggest.

Section 4.1.2

Vogel's IPR Reference Equation

The single most important equation in two-phase inflow performance — a dimensionless curve derived from simulation that anchors all subsequent IPR construction methods.

▶

Lecture 4.1C: Derivation and Significance of Vogel's Equation

21:00 · HD

Step-by-step walkthrough of Vogel's 1968 simulation study, how 21 reservoir simulations collapsed onto one dimensionless curve, and why the coefficients 0.2 and 0.8 appear in the equation. Includes a live comparison of Vogel vs. field-observed IPRs from three North Sea wells.

Vogel's Reference Curve — The Equation

In 1968, Kermit Vogel published the results of a comprehensive reservoir simulation study of solution-gas drive reservoirs. He ran 21 simulations varying reservoir fluid properties, relative permeability characteristics, and well conditions. Despite this variation, the dimensionless IPR data collapsed remarkably well onto a single reference curve, now universally known as Vogel's equation:

VOGEL'S EQUATION

q_o / q_o,max = 1 − 0.2(p_wf/p_R) − 0.8(p_wf/p_R)²

Terms:
  q_o = oil production rate at flowing pressure p_wf (stb/d)
  q_o,max = maximum oil rate at p_wf = 0 (stb/d) — the Absolute Open Flow
  p_wf = bottom-hole flowing pressure (psi)
  p_R = average reservoir pressure (psi)

Valid conditions: Solution-gas drive reservoir, p_wf ≤ p_b = p_R, skin ≈ 0 (FE = 1.0)

Understanding the Coefficients: 0.2 and 0.8

The two coefficients in Vogel's equation are not arbitrary — they emerge from fitting the simulation results and reflect fundamental characteristics of two-phase flow in solution-gas drive reservoirs.

The 0.2 coefficient (linear term)

Accounts for the laminar (Darcy) component of inflow. At zero drawdown, the curve passes through q/q_max = 1.0 (as expected). The slope at zero drawdown is −0.2/p_R — steeper than the Darcy PI slope because the gas saturation effect adds to the pressure-rate relationship from the start.

The 0.8 coefficient (quadratic term)

Captures the accelerating reduction in k_ro as p_wf decreases and gas saturation builds up. At very low p_wf/p_R, this term dominates, causing the curve to flatten dramatically — meaning large additional drawdown yields diminishing incremental oil rate.

Boundary Conditions — Check the Equation Makes Sense

It's always good practice to verify a new equation at its boundary conditions:

Condition	p_wf/p_R	Vogel Calculation	q_o/q_max	Interpretation
Well shut in (no flow)	1.0	1 − 0.2(1) − 0.8(1)² = 1 − 0.2 − 0.8	0.0	Correct: zero production at reservoir pressure
Moderate drawdown	0.6	1 − 0.2(0.6) − 0.8(0.36) = 1 − 0.12 − 0.288	0.592	59.2% of maximum rate
High drawdown	0.2	1 − 0.2(0.2) − 0.8(0.04) = 1 − 0.04 − 0.032	0.928	93% of AOF — diminishing returns region
AOF (absolute open flow)	0.0	1 − 0 − 0 = 1	1.0	Correct: q_max at zero BHFP

The Theoretical Productivity Index J*

Standing (1970) defined J* as the productivity index that would be measured at very small drawdowns (approaching zero), derived from the slope of the Vogel curve at p_wf = p_R:

KEY RELATIONSHIP — J* and q_max

Differentiating Vogel's equation at p_wf = p_R and rearranging:

J* = −d(q_o)/d(p_wf) |_{p_wf=pR} = q_max × (0.2 + 1.6×1) / p_R = 1.8 × q_max / p_R Therefore: q_max = J* × p_R / 1.8 (stb/d)

J* has units of stb/d/psi and is the theoretically correct PI for a solution-gas drive well with no skin. It is measured or estimated at low drawdown (above or just below p_b).

Starting with Vogel's equation: q_o = q_max [1 − 0.2(p_wf/p_R) − 0.8(p_wf/p_R)²]

Differentiating with respect to p_wf:
dq_o/dp_wf = q_max [−0.2/p_R − 1.6 p_wf/p_R²]

At p_wf = p_R:
(dq_o/dp_wf)|_pR = q_max [−0.2/p_R − 1.6/p_R] = −1.8 q_max/p_R

The productivity index is J* = −dq_o/dp_wf, so: J* = 1.8 q_max/p_R

Rearranging: q_max = J* p_R / 1.8

This is a key result — it allows you to obtain q_max from a measured PI value, without needing to flow the well to very low BHFP.

Vogel ran 21 reservoir simulations with the US Bureau of Mines reservoir simulator in 1968. The simulations varied: oil gravity (16–52° API), solution GOR (0–2000 scf/stb), bubble point pressure (1,230–6,997 psi), reservoir pressure (1,935–7,134 psi), and relative permeability characteristics. Despite this wide range, all 21 normalised IPR curves fell within a narrow band around his reference equation. This universality gave the equation its power — a single dimensionless curve representing all solution-gas drive reservoirs.

The original paper: Vogel, J.V.: "Inflow Performance Relationships for Solution-Gas Drive Wells," JPT (January 1968) 83–92.

Section 3

Constructing the Vogel IPR Curve from a Well Test

In practice you will rarely know q_max directly. This section teaches the step-by-step method to construct a full IPR curve from a single production test data point.

▶

Lecture 4.1D: Constructing the Vogel Curve — Step-by-Step Method

19:30 · HD

Screen-recording walkthrough of the construction procedure using a spreadsheet. Shows the tabular method, graphing technique, and common student errors (wrong axis orientation, forgetting to check p_wf vs p_b). Uses the Karama Field KA-07 data as a live example.

The Three-Step Construction Procedure

Given: reservoir pressure p_R, a measured test rate q_test at measured BHFP p_wf,test (where p_wf,test < p_b), construct the complete Vogel IPR.

STEP-BY-STEP METHOD

Step 1 — Calculate q_max from the test point:
From Vogel's equation rearranged for q_max:

q_max = q_test / [1 − 0.2(p_wf,test/p_R) − 0.8(p_wf,test/p_R)²]

Step 2 — Generate the IPR table:
For a range of p_wf values from p_R down to 0, calculate corresponding q_o using:

q_o = q_max × [1 − 0.2(p_wf/p_R) − 0.8(p_wf/p_R)²]

Step 3 — Plot and read off operating points:
Plot p_wf (y-axis) vs. q_o (x-axis). The intersection with the Tubing Performance Curve (TPC) gives the operating flow rate.

Practical Tabular Construction — Worked Template

Below is a standard tabular format for IPR construction. In practice, a spreadsheet automates this. Note that p_wf steps of p_R/10 give 11 data points — sufficient for a smooth curve.

p_wf (psi)	p_wf/p_R	(p_wf/p_R)²	0.2×(p_wf/p_R)	0.8×(p_wf/p_R)²	[1 − 0.2r − 0.8r²]	q_o = q_max × [...] (stb/d)
p_R	1.00	1.0000	0.200	0.800	0.000	0
0.9 p_R	0.90	0.8100	0.180	0.648	0.172	0.172 q_max
0.8 p_R	0.80	0.6400	0.160	0.512	0.328	0.328 q_max
0.7 p_R	0.70	0.4900	0.140	0.392	0.468	0.468 q_max
0.6 p_R	0.60	0.3600	0.120	0.288	0.592	0.592 q_max
0.5 p_R	0.50	0.2500	0.100	0.200	0.700	0.700 q_max
0.4 p_R	0.40	0.1600	0.080	0.128	0.792	0.792 q_max
0.3 p_R	0.30	0.0900	0.060	0.072	0.868	0.868 q_max
0.2 p_R	0.20	0.0400	0.040	0.032	0.928	0.928 q_max
0.1 p_R	0.10	0.0100	0.020	0.008	0.972	0.972 q_max
0	0.00	0.0000	0.000	0.000	1.000	q_max

COMMON ERRORS TO AVOID

1. Using p_R above p_b: If reservoir pressure exceeds bubble point, do NOT use Vogel alone — you need the composite IPR (Topic 4.2). Using Vogel when p_R > p_b will incorrectly curve the above-bubble-point portion.

2. Wrong axis orientation: Always plot p_wf on the y-axis and q_o on the x-axis. Reversing axes makes it impossible to find the TPC intersection correctly.

3. Not checking the test point: If p_wf,test > p_b, the two-phase flow equation doesn't apply to that test — the skin-corrected Standing method or composite IPR must be used instead.

4. Confusing q_max with q_test: q_max is the hypothetical AOF at zero BHFP — never a measured value from a routine well test.

When the reservoir pressure p_R significantly exceeds the bubble point p_b, the IPR has two sections: a straight-line (Darcy) segment from p_wf = p_R down to p_b, then a curved (Vogel) segment below p_b. The two segments meet at the bubble point — creating a visible kink. This composite IPR is covered in Topic 4.2.

If p_R ≈ p_b (the reservoir is currently at or just below bubble point — common in moderately depleted fields), the pure Vogel equation applies throughout the entire flowing pressure range down to zero.

Section 4

q_max, AOF and Practical Rate Predictions

Calculating the maximum deliverable rate and using it to predict production at any chosen operating BHFP — including for artificial lift design.

Three Ways to Determine q_max

In practice, q_max is never measured directly — flowing a well at zero BHFP is physically impossible. It is always a calculated or extrapolated quantity. There are three standard approaches:

Method A — From a Test Point

Rearrange Vogel: q_max = q_test / [1 − 0.2r − 0.8r²] where r = p_wf,test/p_R. Best if a reliable stabilised test below p_b is available.

Method B — From J*

q_max = J* × p_R / 1.8. Use J* from a Darcy PI test conducted at high BHFP (above p_b) or from the Darcy radial flow equation. Useful when only linear PI data is available.

Method C — From Reservoir Simulation

Simulator outputs predicted IPR directly. q_max is read from the curve at p_wf = 0. Vogel's equation provides a check. Most rigorous — preferred for field development planning.

Predicting Rate at Any BHFP — The Forward Calculation

Once q_max is known, predicting the oil production rate at any specific operating BHFP is a simple substitution into Vogel's equation:

q_o = q_max × [1 − 0.2(p_wf/p_R) − 0.8(p_wf/p_R)²]

EXAMPLE CALCULATIONRate prediction at target BHFP

Given: p_R = 3,600 psi, well tested at q = 820 stb/d with p_wf = 1,800 psi.
Required: Predicted rate at artificial lift target p_wf = 400 psi.

Step 1 — Calculate r_test = p_wf,test / p_R = 1800 / 3600 = 0.500 Step 2 — Calculate Vogel bracket at test: [1 - 0.2(0.5) - 0.8(0.5)²] = 1 - 0.10 - 0.200 = 0.700 Step 3 — Calculate q_max: q_max = 820 / 0.700 = 1,171 stb/d Step 4 — Calculate r at target: r = 400 / 3600 = 0.1111 Step 5 — Calculate q at target BHFP: Vogel bracket = 1 - 0.2(0.1111) - 0.8(0.1111)² = 1 - 0.0222 - 0.00988 = 0.968 q_o = 1,171 × 0.968 = 1,134 stb/d Answer: Lowering BHFP from 1,800 to 400 psi via artificial lift increases production from 820 → 1,134 stb/d (+314 stb/d, +38%)

Note how the large drawdown increase (1,400 psi extra) yields only a 38% rate uplift — this is the diminishing returns effect of the curved Vogel IPR at low BHFP.

Sensitivity: How Much Does Rate Increase With Lower BHFP?

This is arguably the most important practical use of the Vogel IPR — answering "what rate can we achieve if we install an ESP and pull BHFP down to X psi?" The table below shows the incremental gain for a reference well:

BHFP p_wf (psi)	p_wf/p_R	Vogel Factor	q_o (stb/d)	Incremental vs. 1,800 psi BHFP	Notes
1,800	0.500	0.700	820	—	Natural flow (test point)
1,400	0.389	0.789	924	+104 stb/d	Moderate ESP uplift
1,000	0.278	0.862	1,010	+190 stb/d	—
600	0.167	0.922	1,080	+260 stb/d	—
400	0.111	0.968	1,134	+314 stb/d	Target BHFP
200	0.056	0.987	1,156	+336 stb/d	Minimal gain vs 400 psi
0	0.000	1.000	1,171	+351 stb/d (AOF)	Theoretical maximum

DESIGN INSIGHT — Diminishing Returns

Going from natural flow BHFP (1,800 psi) to 400 psi delivers 314 stb/d additional production (89% of the maximum possible gain). Pulling BHFP all the way to 200 psi adds only 22 stb/d more. The Vogel curve tells you exactly where diminishing returns set in — critical for optimising ESP sizing (cost vs. benefit of extra drawdown capability).

The conventional measured PI from a well test (J = q/(p_R − p_wf)) will be higher than J* if the test BHFP is above p_b (single-phase, no curvature effect). It will be lower than J* if the test BHFP is below p_b and significant curvature is present, because the straight-line PI underestimates the true low-drawdown slope of the IPR.

For this reason, Vogel's J* (defined as the limit of dq/dp_wf as drawdown approaches zero) is the theoretically correct measure of a well's intrinsic deliverability, independent of where on the IPR the well happened to be tested.

Section 6

Worked Examples

Three fully worked problems of increasing complexity, matching PBL problem-set style. Attempt each calculation before revealing the solution.

WORKED EXAMPLE 4.1-ABasic Vogel IPR Construction

Scenario: A well in the Dunbar Field (solution-gas drive, depleted below p_b) has been tested. Given data:
  • Average reservoir pressure p_R = 2,850 psi (from latest pressure build-up test)
  • Bubble point pressure p_b = 3,200 psi (p_R is currently below p_b)
  • Production test: q_test = 640 stb/d at p_wf = 1,140 psi

Tasks: (a) Calculate q_max. (b) Tabulate IPR at 10 BHFP steps. (c) What rate is predicted at p_wf = 500 psi?

SOLUTION 4.1-A Given: pR = 2,850 psi, q_test = 640 stb/d, p_wf,test = 1,140 psi (a) Calculate q_max: r_test = 1,140 / 2,850 = 0.4000 Vogel bracket = 1 - 0.2(0.4000) - 0.8(0.4000)² = 1 - 0.0800 - 0.1280 = 0.7920 q_max = q_test / 0.7920 = 640 / 0.7920 = 808 stb/d (b) IPR table (selected points): p_wf r Vogel factor q_o (stb/d) 2,850 1.000 0.000 0 2,565 0.900 0.172 139 2,280 0.800 0.328 265 1,995 0.700 0.468 378 1,710 0.600 0.592 478 1,425 0.500 0.700 566 1,140 0.400 0.792 640 ← test point ✓ 855 0.300 0.868 701 570 0.200 0.928 750 285 0.100 0.972 786 0 0.000 1.000 808 = q_max (c) Rate at p_wf = 500 psi: r = 500 / 2,850 = 0.1754 Vogel bracket = 1 - 0.2(0.1754) - 0.8(0.1754)² = 1 - 0.0351 - 0.02461 = 0.9403 q_o = 808 × 0.9403 = 760 stb/d ANSWER: q_max = 808 stb/d; rate at 500 psi BHFP = 760 stb/d

WORKED EXAMPLE 4.1-BAOF from J* and Artificial Lift Assessment

Scenario: Karama Field Well KA-07. A drill stem test conducted during completion gave a PI of J = 0.72 stb/d/psi measured at BHFP = 4,600 psi (above p_b = 4,500 psi). Reservoir pressure is p_R = 5,100 psi.

The field engineer wants to evaluate whether installing an ESP to pull BHFP down to 2,000 psi will achieve the target rate of 2,800 stb/d during the depletion phase when p_R has declined to 3,800 psi (below p_b).

Note: At p_R = 3,800 psi, assume J* ≈ 0.68 stb/d/psi (corrected for depletion per Topic 4.6 methods). Use this J* to answer the question.

SOLUTION 4.1-B Future conditions: pR = 3,800 psi, J* = 0.68 stb/d/psi Step 1 — Calculate q_max at future conditions: q_max = J* × pR / 1.8 = 0.68 × 3,800 / 1.8 = 2,584 / 1.8 = 1,436 stb/d Step 2 — Calculate rate at target BHFP of 2,000 psi: r = 2,000 / 3,800 = 0.5263 Vogel bracket = 1 - 0.2(0.5263) - 0.8(0.5263)² = 1 - 0.10526 - 0.22159 = 0.6731 q_o = 1,436 × 0.6731 = 967 stb/d Step 3 — Even at AOF (pR = 3,800, pwf = 0): q_max = 1,436 stb/d CONCLUSION: Even at zero BHFP (AOF), the maximum deliverable rate at p_R = 3,800 psi is only 1,436 stb/d — well below the 2,800 stb/d target. At p_wf = 2,000 psi (ESP target), the rate is only 967 stb/d. RECOMMENDATION: The well cannot meet the 2,800 stb/d target from inflow alone during depletion. Options to investigate: 1. Stimulation to increase J* (acidize or frac, see Topic 4.3) 2. Revise rate target based on realistic IPR 3. Consider infill drilling This is exactly the kind of analysis the PBL problem set requires.

WORKED EXAMPLE 4.1-CBack-Calculating Reservoir Pressure from Two Test Points

Scenario: A well was tested twice at different drawdowns, both below the bubble point. Test 1: q₁ = 1,050 stb/d at p_wf1 = 2,400 psi. Test 2: q₂ = 1,420 stb/d at p_wf2 = 1,200 psi. Estimate the reservoir pressure p_R and q_max.

Note: This scenario arises when you don't have a recent pressure build-up test. The Vogel curve shape fixes the relationship between two test points and p_R — enabling a simultaneous solve.

SOLUTION 4.1-C From Vogel's equation, for any two test points: q₁/q_max = 1 - 0.2(p_wf1/pR) - 0.8(p_wf1/pR)² ...[1] q₂/q_max = 1 - 0.2(p_wf2/pR) - 0.8(p_wf2/pR)² ...[2] Dividing [1] by [2] eliminates q_max: q₁/q₂ = [1 - 0.2r₁ - 0.8r₁²] / [1 - 0.2r₂ - 0.8r₂²] where r₁ = 2400/pR and r₂ = 1200/pR 1050/1420 = 0.7394 0.7394 × [1 - 0.2(1200/pR) - 0.8(1200/pR)²] = [1 - 0.2(2400/pR) - 0.8(2400/pR)²] Let x = 1200/pR, so 2400/pR = 2x: 0.7394(1 - 0.2x - 0.8x²) = 1 - 0.4x - 3.2x² 0.7394 - 0.1479x - 0.5915x² = 1 - 0.4x - 3.2x² Rearranging: 2.6085x² + 0.2521x - 0.2606 = 0 Quadratic formula: x = [-0.2521 + √(0.0635 + 2.718)] / 5.217 = [-0.2521 + √2.7815] / 5.217 = [-0.2521 + 1.6678] / 5.217 = 1.4157 / 5.217 = 0.2714 Therefore: 1200/pR = 0.2714 pR = 1200 / 0.2714 = 4,420 psi Calculate q_max using Test 1: r₁ = 2400/4420 = 0.5430 Vogel bracket = 1 - 0.2(0.5430) - 0.8(0.5430)² = 1 - 0.1086 - 0.2357 = 0.6557 q_max = 1050 / 0.6557 = 1,601 stb/d Verify with Test 2: r₂ = 1200/4420 = 0.2715 Vogel bracket = 1 - 0.2(0.2715) - 0.8(0.2715)² = 1 - 0.0543 - 0.0590 = 0.8867 q₂ = 1601 × 0.8867 = 1,420 stb/d ✓ (matches test 2) ANSWERS: pR ≈ 4,420 psi; q_max ≈ 1,601 stb/d

SELF-STUDY CHALLENGE

Practice Problem 4.1-D: A well produces 780 stb/d at BHFP = 1,500 psi with p_R = 3,000 psi. The production engineer proposes installing a gas-lift system capable of pulling BHFP to 800 psi. Calculate: (a) q_max, (b) current J*, (c) rate increase achievable with gas lift, (d) percentage of AOF already being achieved at natural flow. Compare your answer with the simulator output before proceeding to the quiz.

Section 7

Field Applications, Limitations & Engineering Judgement

Vogel's method is powerful but bounded. Understanding where it holds, where it breaks down, and how engineers adapt it in practice is essential for professional-level competency.

▶

Lecture 4.1E: Field Case Studies — IPR Uncertainty and Engineering Decisions

23:00 · HD

Three real field case studies: (1) A North Sea well where Vogel significantly overestimated production due to high skin — leading to artificial lift undersizing. (2) A Middle East field where using Vogel correctly identified ESP potential, increasing production 60%. (3) A case where below-bubble-point testing during DST was mistakenly treated as steady-state, and the corrected Vogel analysis changed the development well count. Key engineering lessons from each case.

Where Vogel Applies: Valid Conditions

✓ Vogel Applies When:

• Drive mechanism is primarily solution-gas drive (dissolved gas liberating from oil)
• Reservoir pressure at or below bubble point (p_R ≤ p_b)
• Well skin is near zero (S ≈ 0, Flow Efficiency ≈ 1.0)
• Single-pay, reasonably homogeneous reservoir
• Pseudo-steady state flow has been established
• No significant water production changing relative permeability

✗ Vogel Requires Modification When:

• Skin is significantly positive or negative (use Standing's extension — Topic 4.3)
• p_R is above p_b (use composite IPR — Topic 4.2)
• Water drive or active waterflood (affects relative permeability differently)
• Gas cap drive as primary mechanism
• Naturally fractured reservoirs (k_ro behaviour differs)
• Highly heterogeneous multilayer completions

FIELD EXAMPLEForties Field, North Sea — Vogel IPR in Practice

The Forties Field provides a textbook example of solution-gas drive IPR management. As reservoir pressure declined from initial conditions (above p_b) into the two-phase regime, production engineers observed that wells on natural flow increasingly underperformed relative to straight-line PI predictions. Vogel-based IPR analysis revealed that actual q_max was approximately 58–62% of the linear PI extrapolation — consistent with Vogel's theoretical 55% figure. This insight drove the early introduction of ESPs on key producers and the revision of field production forecasts downward, avoiding significant capital misallocation.

Common Engineering Pitfalls and How to Avoid Them

Vogel's equation assumes pseudo-steady state (PSS) flow — the pressure disturbance has reached the drainage boundary. Many production tests, particularly short DSTs, are conducted in the transient flow regime. If q_test was measured during transient flow, it will be higher than the equivalent PSS rate, causing q_max to be overestimated.

Fix: Check the stabilisation time using: t_s = 948 φ µ c_t r_e² / k. Only use test data where flow has stabilised, or apply transient corrections. This issue is more acute in low-permeability reservoirs.

Vogel's original simulations assumed zero skin (undamaged, normally completed wells). Applying the raw Vogel equation to a well with S = +15 (heavy formation damage) will underpredict the improvement achievable from remediation. The Standing extension (Topic 4.3) corrects for this using a Flow Efficiency (FE) factor, which captures the ratio of the actual flowing pressure to the ideal undamaged flowing pressure.

Rule of thumb: Vogel is reliable for −3 ≤ S ≤ +3. Outside this range, Standing's extension is mandatory.

The reservoir pressure p_R in Vogel's equation must be referenced to the same datum depth as the BHFP p_wf. If p_R was measured by an RFT tool at a different depth, or if the pressure build-up datum differs from the perforation midpoint, a depth correction must be applied: Δp = SG × 0.433 psi/ft × (depth difference in ft vertical).

Failure to apply this correction introduces systematic error in all subsequent IPR calculations — particularly significant in deviated or highly dipping wells.

If the producing GOR significantly exceeds the solution GOR at current reservoir pressure, the well may be producing free gas from a gas cap or fracture network — not just liberated solution gas. In this case, the Vogel equation (which models only solution-gas drive) is no longer appropriate without modification. A separate free gas IPR must be constructed and combined with the oil IPR.

Quantifying IPR Uncertainty — Why a Range of Curves Matters

In field practice, q_max and J* are never perfectly known. Uncertainty arises from measurement accuracy of flowing pressures and rates, whether the test achieved pseudo-steady state, variations in skin over time, and representativeness of relative permeability data. A rigorous engineering analysis always presents a range of IPR curves — P10 (pessimistic), P50 (base), and P90 (optimistic) — rather than a single line. This uncertainty quantification is required for any artificial lift design, ESP selection, or production forecast.

ENGINEERING BEST PRACTICE

When constructing a Vogel IPR for a well that will be artificial-lifted, always:

1. Verify p_wf,test < p_b before applying Vogel (not just near p_b)
2. Use at least two test points to cross-check q_max
3. Bound the uncertainty with P10/P50/P90 curves
4. Check that the test was conducted at stabilised PSS conditions
5. Confirm skin is within Vogel's valid range (S ≈ 0); if not, use Standing's extension
6. For multilayer completions, consider zone-by-zone IPR summation before applying Vogel to the composite

IPR Method	Drive Mechanism	p_R vs p_b	Skin	Use Case
Darcy Linear PI	Any	p_R >> p_b	Any	Single-phase inflow; above bubble point
Vogel (pure)	Solution-gas drive	p_R ≤ p_b	≈ 0	Two-phase throughout; undamaged well
Composite (4.2)	Solution-gas drive	p_R > p_b	≈ 0	Above and below bubble point sections
Standing Extension (4.3)	Solution-gas drive	p_R ≤ p_b	Non-zero	Damaged or stimulated wells
Fetkovich (4.5)	Any (empirical)	Any	Any	Multirate test available; empirical fit

Assessment · Topic 4.1

Knowledge Check — Two-Phase Flow: Vogel's Method

10 questions covering all sub-topics. Immediate feedback with explanations. Aim for 80% before proceeding to Topic 4.2.

1. When does a well's IPR begin to deviate from a straight-line (Darcy) relationship?

Correct. IPR curvature arises because, below p_b, gas evolves from solution, increases gas saturation near the wellbore, and reduces relative permeability to oil (k_ro). This causes rate to increase more slowly than drawdown. The Darcy linear PI assumes constant single-phase permeability, which is violated once gas liberates.

2. In Vogel's reference equation: q_o/q_max = 1 − 0.2(p_wf/p_R) − 0.8(p_wf/p_R)², what does q_max represent physically?

Correct. q_max is the Absolute Open Flow Potential (AOF) — the hypothetical rate at p_wf = 0. It is never measured directly but is calculated by back-calculation from a well test or from J*. It sets the upper bound on deliverability for a given reservoir pressure.

3. A well has p_R = 3,200 psi and tested at q = 950 stb/d with p_wf = 1,600 psi. Calculate q_max (to nearest 10 stb/d).

Correct. r = 1,600/3,200 = 0.500. Vogel bracket = 1 − 0.2(0.5) − 0.8(0.25) = 1 − 0.10 − 0.20 = 0.700. q_max = 950/0.700 = 1,357 ≈ 1,360 stb/d.

4. The relationship between J* (theoretical PI) and q_max in Vogel's framework is:

Correct. Differentiating Vogel's equation at p_wf = p_R gives the slope dq/dp_wf = −1.8 q_max/p_R, which equals −J*. Therefore q_max = J* × p_R / 1.8.

5. At what fraction of p_R does the Vogel IPR curve produce 70% of q_max? (Select closest answer)

Correct. From the Vogel table: at p_wf/p_R = 0.50, the Vogel factor = 1 − 0.2(0.5) − 0.8(0.25) = 0.700, meaning q_o/q_max = 0.70. This is a key reference point worth memorising for rapid field estimates.

6. A completion engineer installs an ESP capable of pulling BHFP from 2,000 psi to 500 psi. With p_R = 4,000 psi and q_max = 2,200 stb/d, how much production is gained? (Calculate then select)

Correct. At p_wf=2,000: r=0.500, factor=0.700, q=2200×0.700=1,540 stb/d. At p_wf=500: r=0.125, factor=1−0.025−0.0125=0.9625, q=2200×0.9625=2,118 stb/d. Gain = 2,118 − 1,540 = 578 ≈ 580 stb/d. Option B (430) is closest among options; however 580 is the accurate answer — this illustrates why calculator accuracy matters.

7. Vogel's original simulation study found that actual q_max for solution-gas drive wells is approximately what fraction of the maximum rate predicted by extending the straight-line PI to zero BHFP?

Correct — 55%. This is the most practically important number from Vogel's work. At p_wf = 0, the Vogel curve gives q_max, while the straight-line PI (slope determined at low drawdown) extrapolated to zero gives q_Darcy,max = J* × p_R. Since q_max = J* × p_R / 1.8, the ratio is 1/1.8 = 0.556 ≈ 55%.

8. Which scenario INVALIDATES the direct application of Vogel's pure equation (without modification)?

Correct — Option D. When p_R > p_b, the IPR has two distinct sections: a straight-line Darcy region from p_R down to p_b, and a curved Vogel section below p_b. Applying pure Vogel across the entire pressure range would incorrectly curve the above-bubble-point section. This requires the composite IPR approach (Topic 4.2). Options A, B, and C are all within Vogel's applicable range with minor modifications.

9. A well's IPR is described by Vogel's equation with q_max = 1,800 stb/d at p_R = 3,600 psi. What is J* for this well?

Correct. J* = 1.8 × q_max / p_R = 1.8 × 1,800 / 3,600 = 3,240 / 3,600 = 0.90 stb/d/psi.

10. In the Karama Field PBL scenario, Well KA-07 has a J* of 0.68 stb/d/psi at a future p_R of 3,800 psi. The project team's target rate is 2,800 stb/d. What does Vogel's method tell you about this target?

Correct — Option C. q_max = J* × p_R / 1.8 = 0.68 × 3,800 / 1.8 = 1,436 stb/d. This is the absolute maximum the well can produce regardless of how low BHFP is pulled. Since 2,800 stb/d >> 1,436 stb/d (AOF), the inflow alone cannot achieve the target. Either stimulation (to increase J*) or a revised rate expectation is needed — this is the central decision point in the PBL problem set.

TOPIC 4.1 COMPLETE

Well done on completing Topic 4.1. You can now construct a Vogel IPR from a single well test, calculate q_max and AOF, predict rates at target BHFP, and identify when Vogel requires modification.

Next: Topic 4.2 — Composite IPR: Combining Darcy's linear PI (above bubble point) with Vogel's curve (below bubble point). This is the general case applicable to most producing wells where p_R > p_b.

PBL Task: Using the KA-07 data from the simulator, complete the IPR analysis worksheet in the Module 04 Problem Set, Tasks 1–3, before proceeding to Topic 4.2.