How is reservoir simulation used in field development planning?

I. High-level purpose and where simulation fits in field development planning

Reservoir simulation is the decision engine of field development planning (FDP). It converts geology, fluids, and well/facility concepts into production and injection forecasts to select, optimize, and de-risk the development concept before capital is committed.

• I.1 Establishes how many wells to drill, where to place them, completion type, and how to control them (rates/BHP/smart valves) to maximize recovery and value.
• I.2 Tests depletion, waterflood, gas injection, or EOR options; quantifies incremental recovery and surface handling requirements.
• I.3 Sizes facilities and export capacity (oil, gas, water, compression, water treatment) under realistic reservoir deliverability and constraints.
• I.4 Frames uncertainty envelopes (P10–P50–P90) for volumes and rates, enabling risked economics and phased developments.
• I.5 Underpins reserves/resources classification and provides surveillance plans for closed-loop updates post-sanction.

II. Step-by-step process flow

• II.1 Define objectives, decision levers, and KPIs
  • II.1.1 Objectives: maximize NPV, accelerate cash flow, meet plateau targets, minimize emissions and water cut.
  • II.1.2 Decision levers: well count/spacing, lateral length, vertical targets, completion type (open hole/cased perf/fractures), injection scheme, facility capacities.
  • II.1.3 KPIs: NPV, recovery factor, plateau duration, water/gas handling loads, flaring, energy intensity.
• II.2 Integrate data and build the static model
  • II.2.1 Seismic-guided structural framework; facies and petrophysical models (porosity, permeability, net-to-gross, water saturation).
  • II.2.2 Upscale properties to the simulation grid preserving flow capacity and storage.
• II.3 Prepare fluids (PVT) and SCAL
  • II.3.1 Generate black-oil tables or tune EOS for compositional cases; include impurities (CO2, H2S) if material.
  • II.3.2 Derive relative permeability/capillary pressure (SCAL), including hysteresis and wettability trends.
• II.4 Build the dynamic model
  • II.4.1 Grid design (structured/unstructured, local refinement around wells/fractures).
  • II.4.2 Initialization: contacts, pressure, temperature, aquifer models and transmissibilities.
  • II.4.3 Well representations: vertical/horizontal/multilateral, completions, skin/friction, stimulation, inflow control.
• II.5 History match (if production data exist)
  • II.5.1 Calibrate to rates, pressures, WOR/GOR, RFT/PLT profiles, tracers, interference tests.
  • II.5.2 Use assisted history-matching and ensembles to maintain geologic realism and quantify non-uniqueness.
• II.6 Design and simulate development scenarios
  • II.6.1 Drill sequencing, well spacing/sweep patterns, injector/producer ratios, lift methods, artificial lift timing.
  • II.6.2 EOR screening: WAG, miscible/immiscible gas, polymer/surfactant; evaluate incremental recovery and facility impacts.
• II.7 Couple surface network and facilities constraints
  • II.7.1 Integrate tubing/lift performance and surface network to honor backpressure, compression, and water handling limits.
  • II.7.2 Constrain flaring and emissions; include gas reinjection/curtailment strategies.
• II.8 Optimize controls and layout
  • II.8.1 Apply automated optimization for well placement, rates/BHP, and injection allocation under constraints.
  • II.8.2 Use proxy models or streamline-based methods to accelerate screening.
• II.9 Quantify uncertainty and risk
  • II.9.1 Build ensembles varying structure, facies, petrophysics, PVT/SCAL, aquifer strength, faults.
  • II.9.2 Produce P10–P50–P90 forecasts; compute Expected Monetary Value (EMV) and downside protection.
• II.10 Select FDP concept and plan surveillance
  • II.10.1 Choose the concept that maximizes risk-adjusted value; define phased development triggers.
  • II.10.2 Design surveillance (PLT, PTA, 4D seismic) to reduce dominant uncertainties during execution.
• II.11 Closed-loop updates post-sanction
  • II.11.1 Calibrate with new data; update forecasts and adjust drilling and facility debottlenecking plans.

III. Major equipment/components and their functions

• III.1 Reservoir simulators
  • III.1.1 Black-oil and compositional engines for multiphase flow; optional thermal/chemical modules for EOR and heavy oil.
  • III.1.2 Dual-porosity/dual-permeability and discrete fracture options to represent fractured systems and hydraulically fractured wells.
• III.2 Pre-/post-processors
  • III.2.1 Grid builders, upscaling tools, well-path/completion editors, and visualization for diagnostics (streamlines, saturation fronts).
• III.3 Optimization and data assimilation
  • III.3.1 Assisted history match, ensemble Kalman filters/smoothers, global/local optimizers for control and placement.
• III.4 Compute infrastructure
  • III.4.1 Workstations and HPC clusters/cloud for parallel runs, ensemble management, and rapid turnaround.
• III.5 Data inputs
  • III.5.1 Seismic, well logs/cores, well tests, production surveillance (rates/pressures), PVT/SCAL lab results, tracer data.
• III.6 Surface network and facility models
  • III.6.1 Nodal analysis and network solvers to enforce tubing, manifold, compressor, separator, and water plant constraints.
• III.7 Geomechanics coupling (as needed)
  • III.7.1 Compaction, subsidence, fault reactivation risks, and fracture conductivity retention under drawdown.

IV. Key performance drivers (efficiency, cost, safety, emissions)

• IV.1 Model fidelity vs. runtime
  • IV.1.1 Grid resolution and timestep controls must capture key heterogeneities and physics without excessive runtime; use local refinement where it matters (near wells, contacts, fronts).
• IV.2 Subsurface data quality
  • IV.2.1 PVT/EOS and SCAL realism dominate forecast reliability; poor endpoints or wettability assumptions skew waterflood/EOR performance.
• IV.3 Constraints integration
  • IV.3.1 Honoring tubing, facility, and export limits avoids unrealistic plateaus and over-sizing; include downtime/maintenance factors where material.
• IV.4 Economic and environmental metrics
  • IV.4.1 Decision metrics: NPV, payout, unit technical cost; include carbon intensity (kg CO2e/boe) and flaring penalties in scenarios.
• IV.5 Operational safety
  • IV.5.1 Simulated drawdown management limits sanding, coning, and H2S breakthrough; injection pressure control reduces integrity risks.

V. Typical challenges and mitigation strategies

• V.1 Non-uniqueness in history match
  • V.1.1 Mitigation: multi-objective assisted matching with ensembles; constrain with PLT/RFT, tracers, and geological priors.
• V.2 Scale mismatch and heterogeneity
  • V.2.1 Mitigation: robust upscaling, local grid refinement around wells/fractures, transmissibility multipliers validated by diagnostics.
• V.3 Numerical stability and runtime
  • V.3.1 Mitigation: CFL-aware timestepping, well index tuning, solver/preconditioner selection, and proxy/streamline screening before full-physics.
• V.4 Complex fluids and EOR physics
  • V.4.1 Mitigation: EOS/PVT tuning, black-oil equivalents where acceptable, and stepwise physics activation with lab calibration.
• V.5 Fractures and compartmentalization
  • V.5.1 Mitigation: dual-porosity/dual-permeability or DFN hybrids; scenario bracketing of fault transmissibility and sealing behavior.
• V.6 Surface–subsurface coupling gaps
  • V.6.1 Mitigation: two-way coupling with network solvers; rate control strategies aligned with compressor and water plant curves.
• V.7 Uncertainty communication
  • V.7.1 Mitigation: P10–P50–P90 forecast ribbons, tornado charts for sensitivities, and risked economics to prevent over/under-sizing.

VI. Why this matters economically and operationally

• VI.1 Avoids mis-investment: right-sizes facilities and well count to realistic reservoir capacity, preventing stranded capex or bottlenecks.
• VI.2 Accelerates cash flow: optimized well placement and controls increase early-time rates and extend plateau.
• VI.3 Increases ultimate recovery: informed sweep patterns and EOR timing lift recovery factor by several percentage points, often worth hundreds of millions on mid-size fields.
• VI.4 Reduces operating costs and emissions: better water/gas management lowers lifting cost and flaring, improving carbon intensity.
• VI.5 Supports reserves booking and phased decisions: credible forecasts and uncertainty ranges underpin reserves classification and gating.

Key equations used in FDP-focused reservoir simulation

• Flow in porous media (Darcy’s law)

  Single-phase: \( q = - \dfrac{k A}{\mu} \dfrac{\mathrm{d}p}{\mathrm{d}x} \). Multiphase: \( q_\alpha = - k \, k_{r\alpha}(S_\alpha) \dfrac{A}{\mu_\alpha} \left(\dfrac{\mathrm{d}p_\alpha}{\mathrm{d}x} - \rho_\alpha g \dfrac{\mathrm{d}z}{\mathrm{d}x}\right) \).

• Material balance (control-volume form)

  For phase \(\alpha\): \( \dfrac{\partial}{\partial t} \left( \phi \rho_\alpha S_\alpha \right) + \nabla \cdot \left( \rho_\alpha \mathbf{v}_\alpha \right) = q_\alpha^{\text{well}} \), where \( \mathbf{v}_\alpha \) follows Darcy’s law.

• Fractional flow (waterflood diagnostics)

  \( f_w(S_w) = \dfrac{ \dfrac{k_{rw}(S_w)}{\mu_w} }{ \dfrac{k_{rw}(S_w)}{\mu_w} + \dfrac{k_{ro}(S_o)}{\mu_o} } \). Shock and breakthrough times guide injector/producer spacing and pattern efficiency.

• Economic objective (NPV)

  \( \mathrm{NPV} = \sum_{t=1}^{T} \dfrac{ \left[ p_o q_o(t) + p_g q_g^{\text{sales}}(t) - c_{\text{lift}}(t) - c_{\text{OPEX}}(t) - c_{\text{carbon}}(t) \right] - \mathrm{CAPEX}(t) }{ (1+r)^t } \).

• Risked value (Expected Monetary Value)

  \( \mathrm{EMV} = \sum_{i} p_i \, \mathrm{NPV}_i \), using P10–P50–P90 or scenario probabilities to inform concept selection.

How simulation evidence maps to FDP decisions