Every production decision—artificial lift selection, compression strategy, infill drilling timing, facility sizing—depends not just on today's well deliverability but on how it will change over 5, 10, or 20 years of reservoir depletion. This final topic builds the bridge from a single-point IPR to a life-of-field deliverability forecast.
Topics 4.1 through 4.5 developed the tools to characterise a well's inflow performance at a given moment: Vogel's curve, composite IPR, Standing's FE correction, pseudo-pressure for gas, and Fetkovich's empirical fit. All of these were static snapshots — the IPR at current reservoir pressure.
In reality, reservoirs deplete. Pressure falls. Gas saturation builds near the wellbore. Relative permeability to oil worsens. The same well that produced 2,000 stb/d on natural flow in Year 1 may struggle to produce 400 stb/d without aggressive artificial lift by Year 8. Production engineers must be able to predict this trajectory and design for it — not just react to it well by well as performance deteriorates.
This topic covers three quantitative methods for predicting future IPR curves:
Standing's J* scaling method — uses laboratory PVT data (kro, µo, Bo at future pressure) to scale the current PI to future conditions. Most rigorous when PVT data is available.
Fetkovich's depletion method — scales the deliverability coefficient J' linearly with (pe/pi). Simple, requires only current deliverability and pressure ratio. Covered conceptually in Topic 4.5, fully developed here.
Empirical methods — log PI vs. cumulative recovery (Kermit Brown method), pivot-point method (Uhri-Blount), and reservoir simulation output — practical tools for fields with production history.
▶
Lecture 4.6A: Why Future IPR Prediction Defines Field Development Value
21:00 · HD
Opens with a field case study from the Forties Field, North Sea: how an incorrect assumption of constant PI over field life led to a compression facility designed 40% undersized, requiring a costly topsides modification in Year 6. Explains the three principal mechanisms of IPR decline (pressure depletion, relative permeability worsening, viscosity increase), and frames the three prediction methods with their data requirements and uncertainty ranges. Sets up the Karama Field problem context for the final PBL deliverable.
LEARNING OBJECTIVES
After completing Topic 4.6, you will be able to:
1. Identify the three physical mechanisms driving IPR decline over reservoir life and explain how each affects the radial flow equation.
2. Apply Standing's J* scaling method using PVT data to compute the future PI ratio (J*f/J*p) from kro, µo, and Bo at future reservoir conditions.
3. Construct the complete future Vogel IPR from the scaled J* and projected reservoir pressure.
4. Apply Fetkovich's depletion scaling (J'f = J'i × pe,f/pi) to predict future deliverability without PVT data.
5. Use the empirical log-PI vs. cumulative recovery method to project PI decline from production history.
6. Compare all three methods quantitatively and select the most appropriate method for given data availability.
7. Produce a life-of-field IPR family (5–6 pressure stages) and use it to design artificial lift stages, compression timing, and infill well requirements.
8. Articulate the uncertainty range in future IPR predictions and communicate it appropriately to stakeholders.
PREREQUISITE
Required: All of Topics 4.1–4.5. Standing's J* concept (Topic 4.1), Vogel's equation (Topic 4.1), composite IPR (Topic 4.2), Fetkovich depletion (Topic 4.5). PVT concepts: kro, µo, Bo as functions of pressure from Course 02.
PBL CONNECTION — KARAMA FIELD FINAL DELIVERABLE
This topic delivers the final analytical component of the Karama Field problem set. Using the current IPR data developed across Topics 4.1–4.5 for wells KA-07 (oil) and KA-G2 (gas), you will: (a) apply Standing's method using the Karama PVT table to predict KA-07's IPR at p̄ = 3,200, 2,600, and 2,000 psia, (b) apply Fetkovich's method for the same stages and compare, (c) build a life-of-field IPR family for KA-07, (d) determine at which reservoir pressure artificial lift becomes essential to meet the 1,800 stb/d field target, and (e) recommend compression timing for the gas cap wells. This analysis forms the core of the Module 04 Final Integrated Report (graded deliverable due at end of module).
Scope
Three future IPR prediction methods. Life-of-field analysis. Artificial lift timing. Integration of all Module 04 IPR tools.
Integrates
Topics 4.1–4.5. The synthesis topic: every prior IPR method converges here in a life-of-field framework.
Time
~120 min: 40 min reading, 25 min simulation, 35 min worked examples, 20 min quiz.
Section 1
Why IPR Changes Over Field Life — The Three Mechanisms
Before building prediction equations, understanding the physical mechanisms driving deliverability decline provides the engineer with the intuition to judge whether predictions are physically plausible.
▶
Lecture 4.6B: Three Mechanisms of IPR Decline — Animation and Field Evidence
18:00 · HD
Three-part animated lecture. Part 1: pressure depletion and its direct effect on the IPR driving force. Part 2: gas saturation build-up reducing k_ro as solution gas liberates below bubble point — the S_g vs. k_ro relationship animated across depletion. Part 3: viscosity increase and B_o decrease in the depleting oil — PVT curves for Karama Field showing the trajectory. The lecture quantifies each mechanism's contribution to total PI decline for a typical North Sea solution-gas drive reservoir.
Mechanism 1: Reservoir Pressure Depletion
The most direct IPR change: as p̄ falls, the maximum available drawdown (p̄ − pwf,min) shrinks. Even if the PI remains constant, the maximum rate at any given BHFP decreases. For the Darcy segment (above bubble point), this appears as a parallel shift of the IPR line to lower rates. For the Vogel/two-phase segment, it compresses the curve toward lower qmax and lower p̄.
As reservoir pressure falls below the bubble point, solution gas liberates and builds up gas saturation Sg throughout the reservoir. Higher Sg means lower So, which means lower kro (from the relative permeability curves). The effective oil permeability ko = k × kro decreases, directly reducing the PI:
J* ∝ ko/(µoBo) = k × kro/(µoBo)
This is the dominant mechanism of PI decline in solution-gas drive reservoirs. The rate of kro decline with pressure depends critically on the reservoir's relative permeability curves — which vary by rock type, wettability, and pore structure. This is why Standing's method, which uses actual kro from laboratory PVT/core data, is preferred for accuracy.
Mechanism 3: Oil Viscosity Increase and Bo Decrease
As pressure falls below pb, gas evolves from the oil, leaving behind heavier, more viscous oil. µo increases (worsening deliverability). Simultaneously, Bo decreases (the oil shrinks as dissolved gas leaves), slightly partially compensating. The net µoBo product typically increases with depletion, further reducing J*.
Quantifying the Combined Effect
Standing's method captures all three mechanisms through the ratio:
This single ratio encapsulates all three mechanisms, the numerator and denominator are both evaluated from PVT and relative permeability data at the respective reservoir pressures. A typical North Sea solution-gas drive reservoir might show:
Reservoir Pressure
kro
µo (cp)
Bo (rb/stb)
kro/(µoBo)
J* Ratio vs. Present
2,250 psi (present)
0.815
3.11
1.173
0.2238
1.000
2,000 psi
0.740
3.28
1.158
0.1951
0.872
1,800 psi
0.685
3.59
1.150
0.1659
0.741
1,500 psi
0.600
3.95
1.135
0.1341
0.599
1,200 psi
0.480
4.38
1.115
0.0983
0.439
THE COMPOUNDING EFFECT
Note how the J* decline is faster than pressure decline alone. At 1,200 psi (53% of 2,250 psi initial), J* has fallen to only 44% of present value. The additional mechanisms (kro and µo changes) amplify the pure pressure effect. This means production at a given BHFP declines faster than reservoir pressure alone, a critical insight for artificial lift design.
In a waterflood or strong water drive reservoir, the mechanism is fundamentally different from solution-gas drive. Rather than gas saturation building up and reducing kro, water saturation increases Sw, which reduces kro through a different relative permeability path. The PI for oil production decreases, but the gross fluid PI may remain roughly constant (or even increase briefly) as water mobility contributes to total inflow.
The difference: (a) fractional flow where oil and water share the same pore paths (captured by the kro/krw ratio at current Sw) and (b) segregated flow where watered-out layers and oil-producing layers contribute separately. In waterflooded fields, reservoir simulation is typically needed to track the evolving Sw distribution and predict future kro, the simple analytical scaling methods of this topic have limited applicability.
Section 2
Standing's J* Scaling Method
The most rigorous analytical method for predicting future PI — using PVT laboratory data to scale the current productivity index to any future reservoir pressure.
▶
Lecture 4.6C: Standing's Scaling Method — Derivation, Application, and Guo Example 3.6
22:00 · HD
Full derivation of the J* ratio from the Darcy radial flow equation. Shows how k_ro/(µ_o B_o) is obtained from laboratory PVT reports and special core analysis, and how it changes with pressure. Walks through the complete Guo et al. Example 3.6 (p̄ declining from 2,250 to 1,800 psi) step by step. Discusses the uncertainty from using correlations vs. measured PVT data, and the critical importance of getting the right k_ro curve for the specific reservoir rock.
The Standing J* Ratio — Complete Derivation
From the Darcy pseudo-steady state radial flow equation, the productivity index J* is:
J* = 0.00708 k h / [µo Bo (ln(re/rw) − 0.75 + S)]
The terms k (absolute permeability), h (pay thickness), re/rw (drainage geometry), and S (skin) are all independent of reservoir pressure for most practical purposes over the production lifetime. Therefore, the ratio of future to present J* reduces to just the ratio of the fluid mobility terms:
Obtain PVT data at current and future reservoir pressuresFrom laboratory PVT reports: kro at current and future So (from relative permeability), µo (cp), Bo (rb/stb). Future conditions correspond to the projected reservoir pressure at the time of interest.
2
Compute [kro/(µoBo)] at each pressureCalculate the fluid mobility ratio at both present and future conditions. This single number encapsulates all three physical mechanisms of PI decline.
Compute qmax,futureqmax,f = J*f × p̄f / 1.8. This is the future Absolute Open Flow at zero BHFP.
5
Tabulate and plot the future Vogel IPRUse q = qmax,f × [1 − 0.2(pwf/p̄f) − 0.8(pwf/p̄f)²] for pwf from p̄f to 0.
6
Overlay present and future IPRs on one plotThe family of curves shows how the "IPR window" shrinks over time — critical input for artificial lift design and production forecasting.
What Data Does Standing's Method Require?
✓ Required Data
• Current J* or PI from a well test
• Laboratory PVT data: kro(So), µo(p), Bo(p) at current reservoir conditions
• Same PVT data at each future reservoir pressure
• Projected future reservoir pressures (from material balance or simulation)
⚠ Key Uncertainties
• kro at future So depends on tracking oil saturation depletion — requires material balance or simulation
• Laboratory kro curves may not represent in-situ wettability
• µo and Bo correlations (if not measured) introduce ±5–15% error
• Assumes skin S stays constant — may not hold for well with scale build-up or fines migration
Present mobility: kro/(µoBo) = 0.815/(3.11×1.173) = 0.815/3.648 = 0.2235
Future mobility: kro/(µoBo) = 0.685/(3.59×1.150) = 0.685/4.129 = 0.1659
J* ratio = 0.1659/0.2235 = 0.742
J*_future = 1.01 × 0.742 = 0.749 stb/d/psi
q_max,future = J*_f × p̄_f / 1.8 = 0.749 × 1,800 / 1.8 = 749 stb/d
Future IPR (Vogel at p̄ = 1,800 psi):
q = 749 × [1 - 0.2(pwf/1800) - 0.8(pwf/1800)²]
Selected points:
pwf(psi) | q(stb/d)
1,800 | 0
1,620 | 129
1,440 | 246
1,260 | 351
1,080 | 444
900 | 525
720 | 594
540 | 651
360 | 696
180 | 729
0 | 749 ← q_max (matches Guo Example 3.6 ✓)
INTERPRETATION:
Present q_max = 1.01 × 2,250/1.8 = 1,263 stb/d at p̄=2,250 psi
Future q_max = 749 stb/d at p̄=1,800 psi
Decline = (1263-749)/1263 = 40.7% in q_max for 20% pressure drop!
The PI ratio (0.742) does not equal the q_max ratio (749/1263=0.593)
because BOTH J* and p̄ declined — the compounding is severe.
Fetkovich (1973) and Eickmeier showed that if PVT data is unavailable, the J* ratio can be approximated as proportional to reservoir pressure ratio (the linear assumption of the Fetkovich depletion method):
The AOFP therefore falls as the square of the pressure ratio — an even faster decline than the simple pressure ratio alone. This approximation is the basis of Fetkovich's depletion method (Section 4.6.3).
Eickmeier further found empirically that in some fields, the exponent m in J*f/J*i = (p̄f/p̄i)m was not unity but ranged from 1 to 3, with m=2.5 consistent with Vogel's simulations and m=1.66 from Eickmeier's House Mountain Field data.
Section 3
Fetkovich's Depletion Method — Full Development
A complete development of Fetkovich's pressure-ratio scaling approach for future IPR — valid without PVT data, straightforward to apply, and useful for sensitivity analysis across a range of depletion scenarios.
The Fetkovich Depletion Framework — Theoretical Basis
From Guo et al., Fetkovich's depletion approach is based on integrating the reservoir inflow equation under the assumption that the pressure function f(p) = kro/(µoBo) is linear in pressure. This gives the initial deliverability coefficient:
FETKOVICH INITIAL DELIVERABILITY COEFFICIENT
J'i = [0.007082 k h / ln(re/rw)] × [kro/(µoBo)]i / (2pi)
The future J' at reservoir pressure pe scales as (Guo Eq. 3.61):
The AOFP decline formula reveals a powerful result: AOFP scales with the cube of the pressure ratio (pe/pi)³, because both J' and the driving force p²e decline with pressure. For a 30% pressure drop (pe/pi = 0.70):
AOFfuture/AOFinitial = (pe/pi)³ = 0.70³ = 0.343
A 30% pressure drop causes a 65.7% AOFP decline according to Fetkovich's method. This "cubic" decline is why deliverability of solution-gas drive reservoirs deteriorates so sharply in the later stages of field life.
Practical Fetkovich Depletion Table
For field planning, a table showing expected AOFP, q at target BHFP, and artificial lift requirement at each pressure stage is the standard deliverable. Using KA-07 current data as the reference:
Stage
p̄ (psi)
pe/pi
J'future
AOFP (stb/d)
q at pwf=800 psi
q vs. Target 1,800 stb/d
Action Required
Current
3,600
1.000
J'i
976
968
−46% below target
ESP immediately
Stage 2
3,000
0.833
0.833 J'i
565
559
−69% below target
ESP + review
Stage 3
2,400
0.667
0.667 J'i
289
286
−84% below target
Infill drilling?
Stage 4
1,800
0.500
0.500 J'i
122
121
−93% below target
Workovers essential
FETKOVICH VS STANDING — WHICH DECLINE IS "CORRECT"?
The table above uses Fetkovich's method, which gives the fastest decline (cubic). Standing's method (using actual PVT) typically gives a slower decline than Fetkovich for the same pressure drop, because the actual kro/(µoBo) function is usually not exactly linear in pressure — it declines more slowly initially and faster later. For decision-making: use Fetkovich as the pessimistic case and Standing as the base case when PVT data is available. The difference between them defines part of your uncertainty band.
Section 4
Empirical Methods & Reservoir Simulation
Three practical methods for predicting future IPR when PVT data is incomplete or when production history is available: the log-PI vs. cumulative method, the pivot-point method, and using simulator outputs.
Method A: Log-PI vs. Cumulative Recovery (Kermit Brown Method)
Standing (1970), citing Kermit Brown, observed that for many solution-gas drive reservoirs, a plot of log(J*) vs. cumulative oil recovery (Np) produces a straight line. This straight-line relationship can be extrapolated to predict future PI at any stage of depletion.
LOG-PI METHOD — Kermit Brown
log(J*) = log(J*initial) − m × Np
or equivalently:
J* = J*initial × 10−m × N_p
where:
Np = cumulative oil production (stb or MMSTB)
m = slope of log(J*) vs. Np plot — determined from historical PI measurements
J*initial = productivity index at the start of the period
Procedure: Plot measured J* values (from periodic well tests) on semi-log paper vs. cumulative production. Fit a straight line. Extrapolate to future Np values (from production forecasts) to predict future J*.
Why the Semi-Log Relationship Holds
The empirical observation of a straight-line log(J*) vs. Np relationship has its roots in the reservoir physics: as cumulative production increases (pressure decreases), the gas saturation builds up approximately proportionally to Np, and kro decreases roughly exponentially with increasing Sg. Since kro dominates the PI, the log-linear relationship emerges. It holds best in the middle phase of solution-gas drive, before aquifer influx or gas cap gas coning significantly alters the flow regime.
Method B: Pivot-Point Method (Uhri and Blount)
Uhri and Blount proposed a graphical pivot-point method that uses two well tests conducted at different reservoir pressures to establish the IPR decline trend, without requiring PI values or PVT data:
PIVOT-POINT METHOD
Given two test points at different reservoir pressures: (q1, pwf1, p̄1) and (q2, pwf2, p̄2) where p̄1 ≠ p̄2:
Step 1: Plot both IPR curves (e.g., Vogel) on the same pwf vs. q chart. Step 2: Find the "pivot point" — the rate at which both curves intersect, or equivalently, the BHFP where the IPR ratio is constant. Step 3: Use the slope of the two-point line to project future IPRs at other reservoir pressures.
This method is particularly useful when well tests are available at two distinct reservoir pressures but no PVT data exists. It captures the empirical decline trend without requiring any fluid property measurements.
Method C: Reservoir Simulation Output
For major development projects and fields with complex geology (layered reservoirs, faults, strong water drive, gas cap), reservoir simulation is the primary tool for future IPR prediction. The simulator tracks So(r,t), kro(So), and p(r,t) in three dimensions, providing a far more accurate prediction than any of the analytical methods.
Simulator-Based IPR
The simulator predicts p̄ and wellbore rate/pressure at each time step. The production engineer back-calculates J* and constructs the IPR at each stage. The simulator implicitly accounts for saturation changes, cross-flow, and multi-layer effects that analytical methods ignore.
Validating Simulator IPR
Well test data (PI measurements at different dates) provides the ground truth for validating simulator IPR predictions. If the simulator-predicted PI diverges from measured PI by more than 20%, the relative permeability curves or the saturation model need recalibration.
Limitation of Simulation
Reservoir simulation requires: a calibrated geological model, well-characterised relative permeability curves, history-matched pressure and production data, and significant computational resources. Not always available for smaller fields or early-stage appraisal.
The empirical log-PI vs. Np relationship breaks down when:
1. Flow regime changes: Transition from solution-gas drive to water influx (aquifer breakthrough) changes the krw/kro balance abruptly — the PI may stabilise or even temporarily increase at water breakthrough due to improved total mobility, before declining again as water cut increases further.
2. Workovers: Stimulation (acid jobs) or recompletions that change skin significantly will cause the PI to jump off the declining trend. Multiple trend lines may be needed for a well that has been worked over several times.
3. Pressure maintenance operations: Waterflooding or gas injection that maintains reservoir pressure will cause the log-PI vs. Np trend to be less steep than in natural depletion, or even flat. The method should only be used within a consistent drive mechanism regime.
4. Very early or very late depletion: In the early transient phase and in the very late depletion phase (when So approaches residual), the log-linear relationship may not hold well.
Section 5
Method Comparison & Selection Framework
Choosing the right future IPR prediction method based on data availability, field maturity, required accuracy, and the specific engineering decision to be made.
Comprehensive Method Comparison
Criterion
Standing's J* Scaling
Fetkovich Depletion
Log-PI vs. Np
Pivot-Point
Reservoir Simulation
Data required
PVT table + current J*
Current J' and p only
Historical J* measurements
Two well tests at different p̄
Full reservoir model
Accuracy (typical)
±10–20%
±20–35%
±15–30% (extrapolation)
±20–40%
±5–15% (calibrated)
Best for
Base case planning, appraisal
Pessimistic scenario, sensitivity
Mature fields with history
Limited data, quick assessment
Complex reservoirs, major FIDs
Handles water influx
Partial (kro from Sw)
No
Only if in PI history
Only if in test history
Yes (full tracking)
Speed
Hours
Minutes
Hours
Hours
Days–weeks
Recommended use
Phase 1 screening
P10 pessimistic case
Field-history validation
Data-poor quick scan
FEED and FID decisions
The P10/P50/P90 Framework for Future IPR Uncertainty
No single future IPR prediction is "correct", they all carry uncertainty from the underlying data and model assumptions. Best practice is to present a family of future IPRs spanning the uncertainty range:
P10 — Pessimistic Case
Fetkovich cubic decline (fastest). Use for stress-testing artificial lift designs and compression capacity. Ensures the installed system can still meet minimum targets even under pessimistic depletion.
P50 — Base Case
Standing's J* scaling with best-estimate PVT data. Use for field development planning, facilities sizing, and economic NPV calculation. The primary design basis.
P90 — Optimistic Case
Slower decline than Standing's (e.g., if reservoir has partial pressure support from an aquifer, or if kro decline is less steep than laboratory curves suggest). Used for upside production scenario and maximum production commitment to buyers.
DESIGN PRACTICE — Which Case to Size For?
Artificial lift equipment (ESP, gas lift mandrels): Size for P50 base case; verify P10 does not cause failure (e.g., ESP running at minimum allowable rate). The P10 pessimistic case defines the minimum rate the equipment must handle without damage.
Surface facilities (separators, flowlines, compression): Size for P90 optimistic on the high end (to avoid bottlenecking), and verify economics still hold at P10 on the low end. These are fixed infrastructure investments that must work across the full uncertainty range.
Economic NPV calculation: Run deterministic cases at P10, P50, P90 and probability-weight them appropriately for project sanction.
Section 6
Artificial Lift Re-Design Over Field Life
Using the future IPR family to optimally time and size artificial lift installations — the key engineering output of future IPR prediction work.
▶
Lecture 4.6D: Life-of-Field Lift Design — From IPR Family to ESP Roadmap
24:00 · HD
Full life-of-field artificial lift design case study for a North Sea oil well. Shows how the IPR family (Years 0, 3, 6, 10) intersects with the natural flow TPC and three successively more aggressive ESP curves. Demonstrates: (1) natural flow onset identification (when TPC no longer intersects IPR), (2) first-stage ESP selection, (3) ESP replacement timing as reservoir depletes, (4) the economic trigger point where well abandonment becomes preferable to further lift investment. Includes a Karama Field worked example with specific equipment recommendations.
The Life-of-Field Artificial Lift Framework
A life-of-field lift design uses the future IPR family to answer four critical questions:
Q1: When does natural flow fail?
The point at which the natural flow TPC (tubing performance curve at minimum wellhead pressure) no longer intersects the current IPR. This defines the "artificial lift onset date." Beyond this date, production rate falls rapidly without intervention.
Q2: What is the target BHFP trajectory?
Define the target BHFP at each depletion stage to maintain the production rate target. From the future IPR: pwf,target(p̄future) = f(qtarget, IPR). This defines the "BHFP roadmap" that artificial lift must achieve.
Q3: How many lift stages are needed?
As reservoir pressure declines, the target BHFP must be progressively reduced to maintain rate. Each "ESP upgrade" — higher head, larger capacity, or deeper set — represents a lift stage. The number of stages = number of times the BHFP target must be reduced.
Q4: When is the economic abandonment limit?
The economic abandonment rate (qecon) defines when lifting costs exceed revenue. From the future IPR: find the latest stage at which q(at BHFP=pwf,min) > qecon. Beyond this, the well is economically depleted.
The intersection of the natural flow TPC (tubing performance curve, fixed at wellhead pressure = 150 psia) with each future IPR gives the natural flow rate at that depletion stage. When the intersection point falls below qtarget = 1,800 stb/d, artificial lift is required.
Using KA-07 Fetkovich depletion (base case):
p̄ (psi)
Natural flow rate (stb/d)
vs. 1,800 target
Action
3,600
642
−1,158 stb/d
ESP immediately
3,000
~373
−1,427 stb/d
ESP essential
2,400
~191
−1,609 stb/d
Higher-head ESP
Note: KA-07 is already below target at current natural flow (642 stb/d vs. 1,800 target). The IPR analysis confirms that the ESP installation recommended in Topics 4.1–4.3 should have been installed at or before first oil.
For gas wells (KA-G2), the future IPR prediction feeds into compression timing decisions rather than artificial lift. As reservoir pressure declines, the wellhead flowing pressure (pwh) also falls. When pwh drops below the pipeline back-pressure or sales contract minimum pressure, compression is needed to boost gas to the required delivery pressure.
The future gas IPR (using Fetkovich depletion: Cfuture = Cinitial × pe,f/pi) predicts the wellhead pressure at the target delivery rate at each depletion stage. Compression is timed to when: wellhead pressure at target rate < pipeline inlet pressure requirement (typically 800–1,500 psia). This defines the "compression onset date" — critical for investment timing and facilities design.
For KA-G2 with current AOFP ≈ 14.8 MMscfd and p̄ = 3,800 psia: at p̄ = 2,800 psia, AOFP ≈ 14.8 × (2800/3800)³ = 14.8 × 0.400 = 5.9 MMscfd. If field target is 12 MMscfd from three wells, compression to lower backpressure becomes critical by Year 5–6 of production.
Section 7 — Interactive Tool
Life-of-Field IPR Simulator
Build and compare future IPR families using Standing's method and Fetkovich's depletion — enter current well data and PVT ratios to generate a full depletion forecast with lift timing analysis.
Future IPR — Standing vs. Fetkovich Comparison INTERACTIVE
Loading...
Scenarios:
• KA-07 base (Guo Ex 3.6): p̄_i=3600, J*=1.01 (approx), p_f=1800, PVT ratio=0.742
• Compare Standing vs Fetkovich — note Fetkovich always gives lower future q_max
• Set PVT ratio=1.0 (no PVT change) — which decline driver remains? (pressure only)
• What PVT ratio is needed for future q_max to equal current natural flow rate?
Life-of-Field Depletion Family — 5 Pressure Stages INTERACTIVE
Loading...
Design questions:
• At which p̄ does natural flow (BHFP~p̄×0.7) fall below target?
• What BHFP must the ESP achieve at each stage?
• When does even AOF fall below the target rate? (well abandonment signal)
Section 8
Worked Examples
Four comprehensive worked examples covering Standing's method, Fetkovich depletion, method comparison, and a complete life-of-field planning exercise for the Karama Field.
WORKED EXAMPLE 4.6-AStanding's J* Scaling — Guo et al. Example 3.6
Given: Current p̄ = 2,250 psi, J*p = 1.01 stb/d/psi. PVT: Present: kro=0.815, µo=3.11 cp, Bo=1.173. Future at p̄ = 1,800 psi: kro=0.685, µo=3.59 cp, Bo=1.150. Tasks: (a) Calculate J*future. (b) Build the future Vogel IPR at p̄=1,800 psi. (c) Compare AOFP present vs. future.
WORKED EXAMPLE 4.6-BFetkovich Depletion — Guo et al. Example 3.7 (Complete)
Given: pi = 2,000 psia, J'i = 5×10⁻⁴ stb/d/psia². Predict complete IPRs at pe = 2,000 (initial), 1,500, and 1,000 psia. Task: Build the full IPR table for all three stages and determine at which stage q at pwf=400 psi falls below 200 stb/d (economic abandonment rate).
WORKED EXAMPLE 4.6-CStanding vs. Fetkovich — Divergence Quantification
Given: A well has current J* = 0.85 stb/d/psi at p̄ = 3,400 psi. Future target: p̄ = 2,200 psi. PVT data gives: kro/(µoBo) ratio future/present = 0.650. Fetkovich assumes linear J' scaling with pressure. Tasks: (a) Compute future J* via Standing's method. (b) Compute future J' via Fetkovich. (c) Compute future qmax by both methods. (d) Quantify the divergence.
SOLUTION 4.6-C — Standing vs. Fetkovich Divergence
STANDING'S METHOD:
J*_future = J*_present × PVT_ratio
= 0.85 × 0.650 = 0.553 stb/d/psi
q_max,Standing = J*_f × p̄_f / 1.8
= 0.553 × 2,200 / 1.8
= 676 stb/d
FETKOVICH METHOD:
J'_future/J'_initial = pe_future/pi = 2200/3400 = 0.647
But Fetkovich's J' and J* are related differently:
J*_Fetkovich = J*_i × (pe_f/pi)
= 0.85 × (2200/3400)
= 0.85 × 0.647 = 0.550 stb/d/psi
q_max,Fetkovich = 0.550 × 2,200 / 1.8 = 673 stb/d
DIVERGENCE:
Standing q_max = 676 stb/d
Fetkovich q_max = 673 stb/d
Difference = 3 stb/d (0.4% — essentially identical!)
WHY THEY AGREE HERE:
The PVT ratio from Standing = 0.650
The pressure ratio from Fetkovich = 2200/3400 = 0.647
These are nearly equal → both methods give almost the same answer!
This means: for THIS particular well, the actual PVT data shows
that k_ro/(µoBo) declines approximately linearly with pressure
(Fetkovich's linear assumption is validated by the PVT data).
WHEN THEY DIVERGE:
If PVT ratio were 0.550 (not 0.650):
Standing q_max = 0.85 × 0.550 × 2200/1.8 = 572 stb/d
Fetkovich q_max = 673 stb/d
Difference = 101 stb/d (15%) — Fetkovich is OPTIMISTIC
If PVT ratio were 0.750:
Standing q_max = 0.85 × 0.750 × 2200/1.8 = 780 stb/d
Fetkovich q_max = 673 stb/d
Difference = 107 stb/d (14%) — Fetkovich is PESSIMISTIC
ENGINEERING LESSON:
Fetkovich assumes k_ro/(µoBo) ∝ p exactly.
If PVT ratio > pressure ratio: Standing > Fetkovich (optimistic PVT)
If PVT ratio < pressure ratio: Standing < Fetkovich (pessimistic PVT)
The divergence is typically 10-30% for moderate pressure drops.
ALWAYS use Standing when laboratory PVT data is available.
WORKED EXAMPLE 4.6-DKA-07 Complete Life-of-Field IPR Analysis (PBL Final Task)
Given: KA-07 current data: p̄ = 3,600 psi, J* = 0.72 stb/d/psi (after acid job to S=+1). Field target: 1,800 stb/d. ESP target BHFP = 2,000 psi currently, reducible to 800 psi minimum. Karama PVT table gives the following [kro/(µoBo)] values at each pressure: 3,600psi: 0.246; 3,000psi: 0.218; 2,400psi: 0.184; 1,800psi: 0.146. Tasks: (a) Compute J* at each future stage using Standing's method. (b) Compute qmax at each stage. (c) Determine rate at ESP BHFP=800 psi at each stage. (d) Identify when rate at BHFP=800 psi falls below target. (e) Make a recommendation.
SOLUTION 4.6-D — KA-07 Life-of-Field Analysis (PBL Final)
Current mobility: [kro/(µoBo)]_present = 0.246 (at p̄=3600)
Step 1 — J* at each stage (Standing):
Stage 1 (p̄=3600): J* = 0.72 stb/d/psi (current, reference)
Stage 2 (p̄=3000): J* = 0.72 × (0.218/0.246) = 0.72 × 0.886 = 0.638
Stage 3 (p̄=2400): J* = 0.72 × (0.184/0.246) = 0.72 × 0.748 = 0.539
Stage 4 (p̄=1800): J* = 0.72 × (0.146/0.246) = 0.72 × 0.594 = 0.428
Step 2 — q_max = J* × p̄ / 1.8 at each stage:
Stage 1: q_max = 0.720 × 3600/1.8 = 1,440 stb/d
Stage 2: q_max = 0.638 × 3000/1.8 = 1,063 stb/d
Stage 3: q_max = 0.539 × 2400/1.8 = 719 stb/d
Stage 4: q_max = 0.428 × 1800/1.8 = 428 stb/d
Step 3 — Rate at BHFP=800 psi using Vogel:
Stage 1 (p̄=3600): r=800/3600=0.222
Factor = 1-0.2(0.222)-0.8(0.222)² = 1-0.044-0.039 = 0.917
q = 1,440 × 0.917 = 1,320 stb/d
Stage 2 (p̄=3000): r=800/3000=0.267
Factor = 1-0.2(0.267)-0.8(0.267)² = 1-0.053-0.057 = 0.890
q = 1,063 × 0.890 = 946 stb/d
Stage 3 (p̄=2400): r=800/2400=0.333
Factor = 1-0.2(0.333)-0.8(0.333)² = 1-0.067-0.089 = 0.844
q = 719 × 0.844 = 607 stb/d
Stage 4 (p̄=1800): r=800/1800=0.444
Factor = 1-0.2(0.444)-0.8(0.444)² = 1-0.089-0.158 = 0.753
q = 428 × 0.753 = 322 stb/d
Step 4 — Compare to target of 1,800 stb/d:
Stage 1: 1,320 stb/d → 73% of target (shortfall 480 stb/d)
Stage 2: 946 stb/d → 53% of target (shortfall 854 stb/d)
Stage 3: 607 stb/d → 34% of target (shortfall 1,193 stb/d)
Stage 4: 322 stb/d → 18% of target (shortfall 1,478 stb/d)
EVEN AT MAXIMUM DRAWDOWN (BHFP=800 psi), KA-07 ALONE
CANNOT MEET THE 1,800 stb/d TARGET AT ANY STAGE.
Maximum achievable rate (at AOF) at each stage:
Stage 1: AOF = 1,440 stb/d → still 360 stb/d below target!
Step 5 — RECOMMENDATIONS:
SHORT TERM (p̄=3600):
→ Install ESP to pull BHFP from natural flow (~4200 psi) to 800 psi
→ Best achievable: 1,320 stb/d (73% of target)
→ CANNOT reach 1,800 stb/d even with maximum drawdown
MEDIUM TERM (p̄=3000):
→ ESP at 800 psi gives only 946 stb/d
→ Consider stimulation (acid job improved to S=−2? FE=1.4?)
→ Recalculate with improved J* after additional stimulation
LONG TERM (p̄<2400):
→ Well is economically marginal — evaluate infill drilling
→ Can KA-08 (proposed infill) supplement?
→ Consider waterflood to maintain pressure
CRITICAL FINDING:
The 1,800 stb/d target CANNOT be met by KA-07 alone at
any point in field life, even at maximum artificial lift.
Field development must rely on KA-07 + infill wells.
This is the central finding of the Module 04 PBL analysis.
MODULE 04 FINAL INTEGRATION — PBL DELIVERABLE
You now have all the tools for the Module 04 Final Integrated Report. Bring together your work from all six topics:
• Topic 4.1: Vogel IPR at current conditions
• Topic 4.2: Composite IPR (Darcy + Vogel)
• Topic 4.3: Standing's FE correction (pre- and post-acid)
• Topic 4.4: Gas well deliverability for KA-G2
• Topic 4.5: Fetkovich comparison for KA-07
• Topic 4.6: Life-of-field prediction (this topic)
Final Report structure:
Section 1: Current well deliverability (KA-07 oil + KA-G2 gas) using best method
Section 2: Stimulation economics (acid job recommendation with payback)
Section 3: Artificial lift requirement (ESP timing and sizing)
Section 4: Life-of-field IPR family (3–4 pressure stages)
Section 5: Gas compression timing (KA-G2 and KA-G3)
Section 6: Conclusions and recommendations to leadership forum
Submit before the Module 04 Leadership Forum presentation.
10 questions covering all three depletion methods, life-of-field design, and the Module 04 synthesis. Target 80%.
1. Standing's J* scaling method predicts future PI using the ratio of [kro/(µoBo)] at future vs. present conditions. Which physical mechanism does this ratio NOT capture?
Correct — D. The J* ratio captures the three PVT-based mechanisms (kro, µo, Bo changes) but the ratio alone does not account for the shrinking maximum drawdown as p̄ declines. The impact of lower p̄ on qmax is separately captured when computing qmax,future = J*future × p̄future/1.8. Both the J* decline AND the p̄ decline compound to reduce qmax — Standing's method correctly accounts for both steps, but the J* ratio alone only gives one of the two factors.
2. Using Standing's method, current [kro/(µoBo)]p = 0.220 and future [kro/(µoBo)]f = 0.165. Current J* = 0.85 stb/d/psi at p̄ = 3,200 psi. Future p̄ = 2,400 psi. Calculate future qmax.
Correct — B: 636 stb/d. J*f = 0.85 × (0.165/0.220) = 0.85 × 0.750 = 0.638 stb/d/psi. qmax,f = J*f × p̄f/1.8 = 0.638 × 2,400/1.8 = 851 stb/d... Wait: 0.638 × 2400 = 1,531 / 1.8 = 851 stb/d. Let me recalculate: 0.85 × (0.165/0.220) = 0.85 × 0.7500 = 0.6375. q_max = 0.6375 × 2400/1.8 = 1530/1.8 = 850 stb/d ≈ option A. Actually option B (636) would require J* = 636 × 1.8/2400 = 0.477, which needs ratio = 0.477/0.85 = 0.561 — doesn't match. The correct answer is approximately 850 stb/d (option A). This question demonstrates the two-step calculation: first get J*_f, then compute q_max = J*_f × p̄_f/1.8.
3. According to Fetkovich's depletion method, if reservoir pressure falls to 70% of initial (pe/pi = 0.70), the AOFP will decline to approximately:
Correct — C: 34%. AOFPfuture = J'i × (pe/pi) × p²e = J'i × p³e/pi. Therefore AOFPf/AOFPi = (pe/pi)³ = 0.70³ = 0.343 ≈ 34%. This cubic relationship means deliverability deteriorates much faster than reservoir pressure itself — a 30% pressure drop causes a 66% AOFP drop. This is the key message that drives the urgency of early artificial lift installation.
4. The log-PI vs. cumulative recovery (Kermit Brown) method requires:
Correct — C. The log-PI method is purely empirical — it plots measured J* values from periodic well tests (e.g., annual PI measurements from pressure buildups) against cumulative oil production Np. If the plot is linear on semi-log coordinates, extrapolation predicts future PI. The method requires historical production data and well tests but no PVT data, making it valuable for mature fields with production history but limited laboratory measurements.
5. When presenting future IPR uncertainty to a project team, the recommended approach is:
Correct — D. No single prediction is "correct" — all methods have uncertainty. The engineering-grade approach is to bracket the uncertainty with P10/P50/P90 cases: Fetkovich gives the fastest decline (pessimistic/P10), Standing with best-estimate PVT gives the base case (P50), and an optimistic case (P90, perhaps based on partial pressure support from an aquifer or optimistic kro) rounds out the uncertainty range. Different engineering decisions use different cases: equipment sizing uses P10, NPV calculation uses P50, maximum production commitment uses P90.
6. KA-07 analysis shows that even at maximum ESP drawdown (BHFP = 800 psia), the well cannot achieve the 1,800 stb/d field target at any depletion stage. What is the most appropriate engineering response?
Correct — C. KA-07's AOFP at current conditions is 1,440 stb/d (best case, at S=+1) — already below the 1,800 stb/d target. No amount of drawdown can exceed the AOFP. The correct responses are: (1) drill additional infill wells to supplement production, (2) implement pressure maintenance (waterflood) to sustain p̄ and extend the period of economic production, or (3) revise the development plan with realistic multi-well production forecasts. The IPR analysis has done its job — it has defined the physical limit and forced a development plan revision.
7. Standing's method requires the value of kro at future reservoir conditions. Where does this value come from?
Correct — B. kro is a function of oil saturation So, obtained from laboratory relative permeability curves (SCAL — Special Core Analysis). At future reservoir conditions, So will be lower (gas saturation Sg higher as solution gas has evolved), so kro will be lower. The future So is estimated from material balance (how much oil has been produced, accounting for gas saturation from liberated solution gas) or from reservoir simulation. This is why Standing's method requires close collaboration between the production engineer and the reservoir engineer who owns the PVT and SCAL data.
8. Fetkovich's depletion method assumes J' scales linearly with (pe/pi). This implies kro/(µoBo) varies:
Correct — C. Fetkovich's derivation assumes the pressure function f(p) = kro/(µoBo) is linearly proportional to pressure (the "straight-line assumption" — his key simplification). This is why J' ∝ (kro/(µoBo))/p scales as p/p = constant × pe/pi. The linear approximation is generally valid over moderate pressure ranges (within 30–40% of initial pressure) but deviates significantly for large depletion, which is why Standing's method with actual PVT data is preferred for accuracy.
9. For a gas well (KA-G2) with AOFP = 14.8 MMscfd at p̄ = 3,800 psia, Fetkovich's depletion predicts AOFP at p̄ = 2,660 psia (70% of initial pressure) as approximately:
Correct — B: ~5.1 MMscfd. Fetkovich AOFP ∝ p³e/pi. AOFPf/AOFPi = (pe,f/pi)³ = 0.70³ = 0.343. AOFPf = 14.8 × 0.343 = 5.08 ≈ 5.1 MMscfd. This confirms: for the Karama Field gas cap wells, a 30% reservoir pressure decline causes a 65.7% AOFP decline — compression from ~Year 3–5 is essential to maintain gas sales contract delivery when p̄ falls below ~3,200 psia.
10. Which sequence of Module 04 tools correctly represents the complete workflow for the Karama Field KA-07 engineering assessment?
Correct — D. The correct Module 04 workflow integrates all six topics: (1) Establish current IPR via Vogel or composite approach [Topics 4.1–4.2]; (2) Quantify skin impact and stimulation benefit via Standing's FE [Topic 4.3]; (3) Verify IPR accuracy with Fetkovich using multirate test data [Topic 4.5]; (4) For gas wells, apply pseudo-pressure LIT analysis [Topic 4.4]; (5) Predict future IPR using Standing's J* scaling and Fetkovich depletion [Topic 4.6]; (6) Design life-of-field artificial lift and compression schedule based on future IPR family [Topic 4.6]. This synthesis is the Module 04 Final Integrated Report deliverable.
MODULE 04 COMPLETE — CONGRATULATIONS
You have completed the full Module 04 IPR Framework. You now possess a professional-grade toolkit covering:
Topic 4.1: Vogel's two-phase IPR for solution-gas drive wells Topic 4.2: Composite IPR for wells above and below bubble point Topic 4.3: Standing's FE correction for damaged and stimulated wells Topic 4.4: Gas well deliverability — pseudo-pressure, non-Darcy, isochronal testing Topic 4.5: Fetkovich's empirical deliverability equation Topic 4.6: Future IPR prediction across field life
Final Deliverable: Complete and submit the Karama Field Module 04 Integrated Report including your IPR analysis, stimulation recommendation, ESP design basis, life-of-field deliverability forecast, and compression timing recommendation. This report goes to the Karama Field Leadership Forum — your professional communication and quantitative analysis skills are on the line.
Good luck, and remember: the best production engineers don't just calculate IPR curves — they use them to make better decisions, faster, with quantified uncertainty.