Reservoir Engineering Formula Guide
A reservoir engineering formula guide for Darcy flow, straight-line productivity index, Vogel IPR, material balance, and original oil in place checks.
By PetroCalcHub Editorial Team | Updated 2026-07-31
linked formula references
verified source links
Engineering context
Reservoir formulas answer different questions at different scales. Darcy and radial-flow relationships connect pressure drop to rock, fluid, and geometry. Productivity index summarizes a well's rate per unit drawdown. Vogel provides an empirical curved inflow relation for solution-gas-drive oil behavior. Material balance reconciles reservoir-wide production and injection with fluid and rock expansion.
A common failure is to choose an equation because its inputs are available rather than because its assumptions describe the reservoir. The workflow below starts with the physical question, identifies pressure and rate conditions, and then uses a simple equation as a transparent baseline before escalating to transient analysis, nodal analysis, or simulation.
Practical workflow
- 1
Use Darcy or radial-flow equations when the flow geometry, fluid properties, and pressure boundary assumptions match the case.
- 2
Use productivity index for quick linear inflow checks where rate is proportional to drawdown.
- 3
Use Vogel IPR for oil-well two-phase inflow screening when the solution-gas-drive assumptions are reasonable.
- 4
Use material-balance calculators as consistency checks against production, pressure, PVT, and aquifer assumptions.
Distinguish flow regime from inflow correlation
Steady-state radial flow assumes a maintained outer pressure, while pseudosteady flow represents a bounded drainage volume after pressure disturbance reaches the boundaries. Transient flow uses time-dependent pressure propagation. These descriptions concern reservoir pressure behavior and should not be confused with choosing a linear or Vogel-shaped IPR.
Record the boundary model, drainage dimensions, net thickness, permeability, viscosity, formation volume factor, wellbore radius, and skin basis with an analytical rate. The logarithmic radius term means geometry errors may be less dominant than permeability or skin, but they are still part of the model definition.
Use PI as a measured diagnostic
Straight-line PI is rate divided by average reservoir pressure minus flowing bottomhole pressure. A PI calculated from one test is most useful when compared with later tests collected on a consistent pressure, rate, completion, and fluid basis. A decline can indicate pressure depletion, mobility change, skin growth, liquid loading, scale, or test inconsistency.
Extrapolating a constant PI to zero flowing pressure can overstate deliverability when pressure crosses bubble point or multiphase flow develops. Outflow performance, facilities, artificial lift, sand, water, gas handling, and drawdown limits also constrain the achievable operating rate.
Apply Vogel within its scope
Vogel's normalized equation relates rate to the ratio of flowing pressure and average reservoir pressure for solution-gas-drive oil wells. The original relationship was developed from simulation behavior and is not a universal law for every oil, completion, pressure state, or drive mechanism.
For an undersaturated reservoir with flowing pressure below bubble point, a composite IPR can use a linear segment above bubble point and a curved segment below it. Calibrating the curve to a representative test point is generally more defensible than assuming an arbitrary absolute open-flow rate.
Keep material balance at reservoir scale
Material balance requires a consistent pressure history, cumulative production and injection, PVT properties, rock and connate-water compressibility, and any aquifer model. It estimates connected in-place volumes under the chosen tank assumptions; it does not by itself locate bypassed oil or describe spatial pressure variation.
Plotting diagnostic terms and checking straight-line behavior can reveal inconsistent data or a missing drive mechanism. A numerical answer without a pressure history and residual review is much weaker than a transparent model that shows which expansion and influx terms dominate.
Worked linear-versus-Vogel comparison
Average reservoir pressure is 3,000 psi. A stabilized single-phase test indicates PI = 1.0 STB/day/psi. For comparison, consider a Vogel curve with qmax = 1,500 STB/day and evaluate both at 1,500 psi flowing pressure.
- 1. Pressure drawdownDrawdown = 3,000 - 1,500 = 1,500 psi
Both comparisons use the same average and flowing pressures.
- 2. Linear IPRq = J x drawdown = 1.0 x 1,500 = 1,500 STB/day
This assumes productivity index remains constant over the full drawdown.
- 3. Vogel rate fractionq/qmax = 1 - 0.2(0.5) - 0.8(0.5)^2 = 0.70
The Vogel normalized curve reduces rate as the pressure ratio changes.
- 4. Vogel rateq = 1,500 x 0.70 = 1,050 STB/day
This is the curved IPR result for the stated qmax, not a forecast of surface deliverability.
Result
Straight-line inflow gives 1,500 STB/day. Vogel gives 1,050 STB/day because the normalized pressure ratio is 0.5 and the Vogel rate fraction is 0.70.
Interpretation
The linear model predicts 1,500 STB/day at the test drawdown, while the selected Vogel curve predicts 1,050 STB/day for the same 1,500 psi flowing pressure. The difference is not an error; it reflects different assumed inflow physics and shows why model selection must be explicit.
Method selection guide
Match the calculation method to the physical question and the evidence available.
| Condition | Use | Why |
|---|---|---|
| Single-phase rate proportional to drawdown | Straight-line productivity index | A constant PI is transparent and can be calibrated directly from a stabilized rate and pressure pair. |
| Solution-gas-drive oil inflow below bubble point | Vogel IPR or an appropriate composite IPR | The curved relationship represents declining oil mobility as gas evolves, within the correlation's assumptions. |
| Changing average pressure and cumulative production | Material balance with a documented tank model | Material balance evaluates reservoir-scale withdrawal and expansion rather than one well's instantaneous inflow. |
Before using the result
- Depth-match and time-match average reservoir pressure, flowing pressure, and rate.
- State whether pressure is above, at, or below bubble point.
- Use rate, viscosity, and formation volume factor on a compatible stock-tank or reservoir basis.
- Separate wellbore skin and completion effects from reservoir pressure behavior where possible.
- Identify the assumed outer boundary and steady, pseudosteady, or transient flow regime.
- Calibrate analytical results against tests, surveillance, material balance, and simulation evidence.
Start with these calculators
These links keep the guide close to the working calculator flow.