Reservoir engineering and stimulation
Volumetric reserves, Darcy inflow, reservoir oil and gas properties, fracturing and acidizing design, voidage replacement by injection.
35 calculators
Reserves and recovery
Stock tank volume of oil initially in place by the volumetric method, in million barrels.
N = C × A × h × φ × (1 − Swc) / Boi / 10⁶Free gas initially in place by the volumetric method, in million standard cubic feet.
G = C × A × h × φ × (1 − Swc) / Bgi / 10⁶Gas dissolved in oil initially in place by the volumetric method, in million standard cubic feet.
Gs = C × A × h × φ × (1 − Swc) × Rs / Bgi / 10⁶Initial oil in place by the volumetric formula of Russian and Kazakh reserves reports (GKZ): oil-bearing area, net oil pay, porosity and oil saturation, shrinkage factor and oil density. The result is in thousand tonnes and thousand cubic metres of stock-tank oil, with the dissolved gas in place.
Q = F · h · K_п · K_н · θ · ρThe recovery factor as the product of the displacement, sweep and flooding efficiencies; the displacement efficiency comes from the initial and residual oil saturations of core floods. Recoverable reserves are that share of the oil in place.
K_выт = (K_н − K_но) / K_нInflow and flow in porous media
Liquid flow rate in linear flow by Darcy's law, in barrels per day.
q = C × k × A × (Pe − Pw) / (μ × L)Well flow rate for radial liquid flow by Darcy's law, in barrels per day.
q = C × π × k × h × (Pe − Pw) / (μ × ln(Re / Rw))Well flow rate for radial flow by Darcy's law, corrected by the oil formation volume factor, in cubic metres per day.
q = C × k × h × (Pe − Pw) / (μ × Boi × ln(Re / Rw))Theoretical productivity index of a vertical well in a bounded drainage area (circle, square, rectangle) from permeability, thickness, fluid properties and skin, and the same without skin. The drainage radius comes from the area per well.
J = 2π · k · h / (μ · B · (½ · ln(4A / (γ · C_A · r_w²)) + s)), γ = 1.781Steady-state productivity index of a horizontal well by Joshi's equation with permeability anisotropy, compared with a vertical well in the same drainage area.
J_h = 2π · k_h · h / (μ · B · (ln((a + √(a² − (L/2)²)) / (L/2)) + (I · h / L) · (ln(I · h / (r_w · (I + 1))) + s)))Reservoir fluid properties (PVT)
Viscosity of degassed oil at standard conditions from its density, for oil with a density of 0.845 to 0.924 g/cm³, in mPa·s.
μ = (C₁ × ρo² / (C₂ − ρo²))²Viscosity of degassed oil at standard conditions from its density, for oil with a density of 0.78 to 0.845 g/cm³, in mPa·s.
μ = (C₃ × ρo² / (C₄ − ρo²))²Viscosity of degassed oil at a given temperature from its laboratory viscosity and the coefficient a, in mPa·s. The coefficient c is chosen from the viscosity: 10 for μ > 1000, 100 for 10 ≤ μ ≤ 1000, 1000 for μ < 10.
μₜ = (1 / c) × (c × μₜ₁)^aCoefficient a for converting the viscosity of degassed oil to another temperature. The coefficients b and c are chosen from the viscosity: b = 0.00252 and c = 10 for μ > 1000; b = 0.00144 and c = 100 for 10 ≤ μ ≤ 1000; b = 0.00076 and c = 1000 for μ < 10.
a = 1 / (1 + b × (t − t₁) × log₁₀(c × μₜ₁))Viscosity of gas-saturated oil at reservoir conditions from the degassed oil viscosity and the correlation coefficients A and B, in mPa·s.
μ = A × μₜ₁^BCoefficients A and B for the viscosity of gas-saturated oil, from the gas-oil ratio at 15 °C and atmospheric pressure.
A = exp(−0.008724 × Gf₁₅ + 0.0000129 × Gf₁₅²); B = exp(−0.004711 × Gf₁₅ + 0.0000083 × Gf₁₅²)Converts the gas content of oil from standard conditions to 15 °C and atmospheric pressure, in m³/m³.
Gf₁₅ = 0.983 × (1 + 5 × αo) × GfThermal expansion coefficient of degassed oil from its relative density.
αo = 10⁻³ × C₁ × (C₂ − ρo)Oil formation volume factor from the gas-oil ratio, for gas-oil ratios up to 400 m³/m³.
Boi = 1 + 0.00305 × GfOil formation volume factor from the gas-oil ratio, for gas-oil ratios above 400 m³/m³.
Boi = 1 + 0.00363 × (Gf − 58)Molecular weight of gas from its composition (in %), and the gas density at normal (0 °C) and standard (20 °C) conditions.
M = (Y₁ × 16.043 + Y₂ × 30.07 + Y₃ × 44.097 + Y₄ × 58.024 + Y₅ × 72.151 + Y₆ × 44.010 + Y₇ × 34.080 + Y₈ × 28.014) / 100; ρ₀ = M / 22.414; ρ₂₀ = M / 24.05Converts the saturation (bubble-point) pressure of oil from reservoir temperature to another temperature, in MPa.
Ps = Pi + (t − ti) / (9.157 + 701.8 / (Gf × (Y₁ − 0.8 × Y₂)))Converts the gas content of reservoir oil from m³/m³ to m³ per tonne at normal conditions (0 °C).
Gf₀ = 1000 × Gf / (ρo × T₂₀ / T₀), T₂₀ = 293.15 K, T₀ = 273 KDensity of gas-saturated oil from the oil and gas densities, the gas-oil ratio and the oil formation volume factor, in kg/m³.
ρ = (1 / Boi) × (ρo + ρg × Gf)Total compressibility of the rock and its fluids from the saturations and the oil, water, gas and pore compressibilities, and the storativity φ·c_t: inputs to well test interpretation, radius of investigation and diffusivity.
c_t = S_o · c_o + S_w · c_w + S_g · c_g + c_fFracturing and acidizing
Formation fracture pressure by Eaton's method: the minimum horizontal stress from the overburden and pore pressures through Poisson's ratio. The gradient is given in kPa/m or as an equivalent density.
P_frac = ν / (1 − ν) · (σ_v − P_p) + P_pExpected surface pressure while pumping: the bottomhole treating pressure minus the hydrostatic head plus the friction in the pipe and through the perforations, and the hydraulic power of the pump fleet.
P_s = G_frac · H − ρ · g · H + ΔP_fr + ΔP_pfPressure drop across the perforations: for limited-entry fracturing and for estimating the number of open perforations from a step-down test.
ΔP_pf = 8 · ρ · q² / (π² · C_d² · N² · d⁴)Proppant mass and slurry volume and density for a frac stage from the clean fluid volume and the proppant concentration in kilograms per cubic metre of fluid (in ppa, pounds added per gallon).
M = c · V_fDimensionless fracture conductivity F_CD, the equivalent skin and effective wellbore radius after fracturing from the Cinco-Ley and Samaniego chart, and the fold increase in productivity in pseudo-steady-state flow.
F_CD = k_f · w / (k · x_f)How much calcite or dolomite a cubic metre of acid dissolves: the gravimetric (β) and volumetric (X) dissolving power from the reaction stoichiometry, the mass of rock dissolved and the CO₂ released.
β₁₀₀ = ν_m · M_m / (ν_a · M_a)Volume of acid to fill the pores and dissolve the carbonate minerals in a ring around the well out to a given radius: in total and per metre of pay. The typical case is the hydrochloric acid preflush ahead of a mud acid treatment of a sandstone.
V = π · (r_s² − r_w²) · h · (φ + (1 − φ) · x_c / X)Waterflooding and pressure maintenance
Current voidage replacement at reservoir conditions: the injected water against the produced oil, water and free gas, and the injection needed for 100 % replacement.
VRR = q_i · B_w / (q_o · B_o + q_w · B_w + q_o · (R − R_s) · B_g)The water-oil mobility ratio from the end-point relative permeabilities, and the water fractional flow at a given water saturation with Corey curves: the stability of the displacement front and the water cut.
M = (k_rw° / μ_w) / (k_ro° / μ_o)Other areas
Results are engineering estimates from standard formulas; for design decisions check them against the codes, project documents and specialists' calculations. The formulas carried over from the original set are unchanged, and their errors are described in the notes.