Skip to content
Brent$117.6Gas TTF$1,192USD453.87 ₸
MunaiHub

Equations: Statistics, volumes and reserves

The equations of the “Statistics, volumes and reserves” lessons in “Petroleum geoscience”, with their conditions and units, and what each symbol means in them.

Arithmetic mean of a sample
xˉ=1N∑i=1Nxi\bar x = \frac{1}{N}\sum_{i=1}^{N} x_i

the same as those of the measured quantity

Holds when: Each measurement represents the same share of the studied volume (the plugs are taken evenly). Otherwise we need weights, such as thicknesses or volumes. All the measurements belong to one population: one layer, one rock type.

Descriptive statistics in geology →
Geometric mean
kˉgeom=(∏i=1Nki)1/N=exp⁡(1N∑i=1Nln⁡ki)\bar k_{\text{geom}} = \left(\prod_{i=1}^{N} k_i\right)^{1/N} = \exp\left(\frac{1}{N}\sum_{i=1}^{N}\ln k_i\right)

the same as for kik_i; ln⁡\ln is the natural logarithm

Holds when: All the values are positive. A single zero turns the mean into zero, so instead of “zeros” we put in the measurement limit of the instrument. The plugs represent equal volumes of rock. Otherwise the logarithms are weighted by thickness.

Descriptive statistics in geology →
Harmonic mean
kˉharm=N(∑i=1N1ki)−1\bar k_{\text{harm}} = N\left(\sum_{i=1}^{N}\frac{1}{k_i}\right)^{-1}

the same as for kik_i

Holds when: All the values are positive. The plugs represent equal thicknesses. Otherwise kˉharm=∑hi/∑(hi/ki)\bar k_{\text{harm}} = \sum h_i / \sum (h_i/k_i).

Descriptive statistics in geology →
Relation between the three means
kˉharm≤kˉgeom≤kˉarith\bar k_{\text{harm}} \le \bar k_{\text{geom}} \le \bar k_{\text{arith}}

any units, the same for all three

Holds when: All the values are positive and are taken with equal weights. Equality holds only if all the values are equal.

Descriptive statistics in geology →
Variance and standard deviation of a sample
s2=1N−1∑i=1N(xi−xˉ)2,s=s2s^2 = \frac{1}{N-1}\sum_{i=1}^{N}\left(x_i - \bar x\right)^2, \qquad s = \sqrt{s^2}

ss is in the units of the measured quantity, s2s^2 is in their squares

Holds when: The measurements belong to one population. The denominator is N−1N - 1, not NN. The deviations are measured from the mean of the same sample, so on average they are slightly smaller than from the true mean. For NN above a few dozen, the difference does not matter.

Descriptive statistics in geology →
Coefficient of variation
Cv=sxˉC_v = \frac{s}{\bar x}

dimensionless; sometimes it is expressed in per cent

Holds when: The quantity is positive and is measured from a natural zero (permeability, porosity, thickness). For temperature in degrees Celsius or for elevations above sea level, the coefficient of variation has no meaning. The mean and the deviation are calculated from the same measurements on the ordinary scale, not the logarithmic one.

Descriptive statistics in geology →
Cumulative shares of flow capacity and storage capacity (Lorenz curve)
FJ=∑j=1Jkjhj∑i=1Ikihi,CJ=∑j=1Jϕjhj∑i=1IϕihiF_J = \frac{\sum_{j=1}^{J} k_j h_j}{\sum_{i=1}^{I} k_i h_i}, \qquad C_J = \frac{\sum_{j=1}^{J} \phi_j h_j}{\sum_{i=1}^{I} \phi_i h_i}

fractions; kk, ϕ\phi and hh are in any units that are the same for all the layers

Holds when: The layers are ordered by decreasing k/ϕk/\phi. Each plug represents its own layer. For equal thicknesses, hh cancels.

Descriptive statistics in geology →
Lorenz coefficient
Lc=2AL_c = 2A

dimensionless, from 0 to 1

Holds when: We find the area from the plotted curve, for example by the trapezoid rule (adding up the areas of the trapezoids between neighbouring points).

Descriptive statistics in geology →
Probability density of the normal distribution
f(x)=1σ2π exp⁡[−(x−μ)22σ2]f(x) = \frac{1}{\sigma\sqrt{2\pi}}\,\exp\left[-\frac{(x - \mu)^2}{2\sigma^2}\right]

xx, μ\mu and σ\sigma are in the units of the measured quantity

Holds when: The quantity can take any values, including negative ones. For porosity or thickness this is only an approximation, and it is good as long as μ\mu is greater than about 3σ3\sigma. The distribution is unimodal (it has one peak) and symmetrical.

Distributions and sample size →
Standardised value
z=x−μσz = \frac{x - \mu}{\sigma}

dimensionless

Holds when: The distribution is normal. Then zz follows the standard normal distribution, which has a mean of 0 and a standard deviation of 1.

Distributions and sample size →
Median, mean and mode of the log-normal distribution
x50=eμ,xˉ=eμ+σ2/2,xmode=eμ−σ2x_{50} = \e^{\mu}, \qquad \bar x = \e^{\mu + \sigma^2/2}, \qquad x_{\text{mode}} = \e^{\mu - \sigma^2}

xx is in the units of the quantity (for example, mD); μ\mu and σ\sigma refer to ln⁡x\ln x and are dimensionless

Holds when: The logarithm of the quantity is normally distributed. We calculate μ\mu and σ\sigma from natural logarithms (base e). For decimal logarithms, μln⁡=2.303 μlg⁡\mu_{\ln} = 2.303\,\mu_{\lg} and σln⁡=2.303 σlg⁡\sigma_{\ln} = 2.303\,\sigma_{\lg}.

Distributions and sample size →
Percentiles and coefficient of variation of the log-normal distribution
xp=eμ+zpσ,Cv=eσ2−1x_p = \e^{\mu + z_p\sigma}, \qquad C_v = \sqrt{\e^{\sigma^2} - 1}

xpx_p is in the units of the quantity; the rest is dimensionless

Holds when: The logarithm of the quantity is normally distributed.

Distributions and sample size →
Standard error of the mean
SE=sNSE = \frac{s}{\sqrt{N}}

the same as those of the measured quantity

Holds when: The measurements are independent: the value of one sample does not predict the value of another. The sample is representative: every part of the population had an equal chance of being included. For strongly skewed data, the normal approximation for the mean works only with tens of samples.

Distributions and sample size →
Confidence interval of the mean
xˉ−t sN  ≤  μ  ≤  xˉ+t sN\bar x - t\,\frac{s}{\sqrt{N}} \;\le\; \mu \;\le\; \bar x + t\,\frac{s}{\sqrt{N}}

the same as those of the measured quantity

Holds when: The same as for the standard error. The data are distributed approximately normally, or there are enough samples. A confidence probability of 95% means this: if we build a hundred such intervals from different samples, about 95 of them will cover the true mean.

Distributions and sample size →
Number of samples needed for a given precision of the mean
N0=(100 t CvP)2N_0 = \left(\frac{100\,t\,C_v}{P}\right)^2

PP is in per cent; the other quantities are dimensionless

Holds when: We are estimating the arithmetic mean. The measurements are independent and the sample is representative. There are enough samples for the mean to be distributed approximately normally. CvC_v is known at least approximately, from a trial series of measurements or by analogy with similar rocks.

Distributions and sample size →
Rule for the number of samples: precision of the mean ±20% at 95% probability
N0=(10 Cv)2,Ps=200 CvNsN_0 = (10\,C_v)^2, \qquad P_s = \frac{200\,C_v}{\sqrt{N_s}}

PsP_s is in per cent of the mean

Holds when: t=2t = 2: probability 95%, more than thirty samples. The first formula refers to a tolerance of 20%, the second to any number of samples at t=2t = 2. For very heterogeneous rocks (CvC_v noticeably greater than 1, for example complex carbonates), a tolerance of ±20% is often unattainable. Then we accept ±50%, and N0=(4 Cv)2N_0 = (4\,C_v)^2.

Distributions and sample size →
Dykstra–Parsons coefficient
VDP=k50−k84.1k50=1−e−σln⁡kV_{DP} = \frac{k_{50} - k_{84.1}}{k_{50}} = 1 - \e^{-\sigma_{\ln k}}

dimensionless, from 0 to 1

Holds when: Permeability is log-normally distributed. For other distributions the coefficient can be misleading. We take k50k_{50} and k84.1k_{84.1} from the line drawn through the points of the probability plot, not from individual samples. The percentage 84.1 is the proportion of samples with a greater permeability. It is the exceedance probability, not the cumulative proportion. On a plot with an axis of cumulative proportion, as in this lesson, we read k84.1k_{84.1} at the 16% mark (z=−1z = -1).

Distributions and sample size →
Confidence interval of the geometric mean
exp⁡(m−t sln⁡N)  ≤  x50  ≤  exp⁡(m+t sln⁡N)\exp\left(m - t\,\frac{s_{\ln}}{\sqrt N}\right) \;\le\; x_{50} \;\le\; \exp\left(m + t\,\frac{s_{\ln}}{\sqrt N}\right)

the limits are in the units of the quantity

Holds when: The quantity is log-normally distributed. The measurements are independent and representative.

Distributions and sample size →
Standard error of the mean for dependent measurements
SE=sN 1+(N−1) ρˉ,Neff=N1+(N−1) ρˉSE = \frac{s}{\sqrt{N}}\,\sqrt{1 + (N - 1)\,\bar\rho}, \qquad N_{\text{eff}} = \frac{N}{1 + (N - 1)\,\bar\rho}

SESE is in the units of the measured quantity

Holds when: We are estimating the mean of a population that is wider than the sampled stretch: of a reservoir from one well, or of an accumulation from several wells. All the measurements have the same variance (the same spread), and the correlation between them is not negative.

Distributions and sample size →
Experimental variogram
γ^(h)=12N(h)∑i=1N(h)[z(xi+h)−z(xi)]2\hat\gamma(h) = \frac{1}{2N(h)}\sum_{i=1}^{N(h)}\left[z(x_i + h) - z(x_i)\right]^2

values in their own units, distances in metres

Holds when: The property is stationary in the studied area: the mean does not change in a regular way, and the differences depend only on the lag. Irregularly placed data are grouped into classes of lag (hh plus or minus a tolerance) and of direction (plus or minus an angular tolerance). An estimate from lags with few pairs is unreliable. A practical rule is to use lags no longer than half the size of the area and to keep several tens of pairs per lag.

Spatial correlation and the variogram →
Spherical, exponential and Gaussian models
γsph(h)=C0+C[3h2a−12(ha)3] for h≤a;C0+C for h>aγexp(h)=C0+C[1−e−3h/a]γgauss(h)=C0+C[1−e−3h2/a2]\begin{aligned} \gamma_{\text{sph}}(h) &= C_0 + C\left[\frac{3h}{2a} - \frac{1}{2}\left(\frac{h}{a}\right)^{3}\right]\ \text{for}\ h \le a;\qquad C_0 + C\ \text{for}\ h > a\\ \gamma_{\text{exp}}(h) &= C_0 + C\left[1 - e^{-3h/a}\right]\\ \gamma_{\text{gauss}}(h) &= C_0 + C\left[1 - e^{-3h^2/a^2}\right] \end{aligned}

the variogram in squared units of the property, distances in metres

Holds when: A model is a sum of permitted structures with non-negative coefficients (nested structures). The nugget effect is a jump at any arbitrarily small h>0h > 0, and the variogram itself is zero at zero. The exponential and Gaussian models are often written with the parameter a/3a/3 or a/3a/\sqrt3. When we move parameters between programs, we check which range is meant.

Spatial correlation and the variogram →
Relation between covariance, correlogram and variogram
C(h)=C(0)−γ(h),ρ(h)=C(h)C(0)=1−γ(h)C(0),C(0)=C0+C=σ2C(h) = C(0) - \gamma(h),\qquad \rho(h) = \frac{C(h)}{C(0)} = 1 - \frac{\gamma(h)}{C(0)},\qquad C(0) = C_0 + C = \sigma^2

in squared units of the property; ρ\rho is dimensionless

Holds when: Second-order stationarity: a constant mean and a finite variance, which the variogram reaches at the sill. With a trend, or with a variogram that has no sill, the covariance is not defined but the variogram is, and we use the variogram.

Spatial correlation and the variogram →
Reduced distance and range in a direction for geometric anisotropy
h′=(h∥amax)2+(h⊥amin)2,a(α)=[cos⁡2(α−α0)amax2+sin⁡2(α−α0)amin2]−1/2h' = \sqrt{\left(\frac{h_{\parallel}}{a_{max}}\right)^2 + \left(\frac{h_{\perp}}{a_{min}}\right)^2},\qquad a(\alpha) = \left[\frac{\cos^2(\alpha - \alpha_0)}{a_{max}^2} + \frac{\sin^2(\alpha - \alpha_0)}{a_{min}^2}\right]^{-1/2}

distances in metres, azimuths in degrees

Holds when: Geometric anisotropy: the sill and the nugget effect are the same in all directions. The axes of the ellipse coincide with the principal directions, usually along the elongation of the bodies and across it. In three dimensions we add a vertical axis with a much smaller range. For zonal anisotropy, we add to the model a structure that acts in one direction only.

Spatial correlation and the variogram →
The ordinary kriging system
∑j=1nλj γ(xi−xj)+μ=γ(xi−x0),i=1,…,n;∑j=1nλj=1\sum_{j=1}^{n}\lambda_j\,\gamma(x_i - x_j) + \mu = \gamma(x_i - x_0),\quad i = 1, \dots, n;\qquad \sum_{j=1}^{n}\lambda_j = 1

variogram in squared units of the property

Holds when: The mean is unknown but constant in the neighbourhood of the estimate (ordinary kriging). When the mean is known, we use simple kriging: without the condition on the sum of the weights, and with the covariance instead of the variogram. The variogram is known and refers to the same quantity and the same support as the data. The estimate is a linear combination of the data. Measurement errors are included in the nugget effect.

Spatial correlation and the variogram →
The kriging estimate and its variance
z∗(x0)=∑i=1nλi z(xi),σK2=∑i=1nλi γ(xi−x0)+μz^*(x_0) = \sum_{i=1}^{n}\lambda_i\,z(x_i),\qquad \sigma_K^2 = \sum_{i=1}^{n}\lambda_i\,\gamma(x_i - x_0) + \mu

as for the property and the variogram

Holds when: The weights and μ\mu are the solution of the ordinary kriging system. The kriging variance is correct if the variogram is correct. The distribution of the errors need not be normal.

Spatial correlation and the variogram →
Oil initially in place (volumetric method)
N=GRV⋅NTG⋅ϕ⋅(1−Sw)BoiN = \frac{\mathrm{GRV}\cdot \mathrm{NTG}\cdot\phi\cdot(1 - S_w)}{B_{oi}}

metric: volumes in m³, the other quantities are dimensionless

Holds when: The pores of the reservoir contain only oil and connate water (the water that has stayed in the pores since the rock was formed). If there is a gas cap (a layer of free gas on top of the oil), the volume above the gas–oil contact is calculated separately, with the formula for gas. The mean ϕ\phi and SwS_w refer to the same volume of reservoir as the NTG. BoiB_{oi} corresponds to the initial reservoir conditions and to the adopted scheme of separating the oil at the surface. We use mean values for the whole accumulation. If the parameters change noticeably across the area, we calculate from maps or from the cells of the geological model and add up the results.

Volumetric estimation of hydrocarbons in place →
Free gas initially in place
G=GRV⋅NTG⋅ϕ⋅(1−Sw)BgiG = \frac{\mathrm{GRV}\cdot \mathrm{NTG}\cdot\phi\cdot(1 - S_w)}{B_{gi}}

metric: volumes in m³

Holds when: The pores of the reservoir are filled with gas and connate water (a gas accumulation or a gas cap). BgiB_{gi} is calculated for the same standard conditions to which the volumes in place are referred. Gas dissolved in the oil is not included here. We count it separately as N⋅RsiN\cdot R_{si}, where RsiR_{si} is the initial solution gas–oil ratio of the oil (how much gas is dissolved in each cubic metre of oil).

Volumetric estimation of hydrocarbons in place →
Gas formation volume factor
Bg=pscTsc⋅Z TpB_g = \frac{p_{sc}}{T_{sc}}\cdot\frac{Z\,T}{p}

any consistent units: pressures in the same units, temperatures absolute (K or °R)

Holds when: The real-gas equation of state pV=ZnRTpV = ZnRT (it links the pressure, volume and temperature of the gas); at standard conditions the gas is taken as ideal (Z=1Z = 1). The pressures are absolute (measured from zero), not gauge (measured above atmospheric pressure). The temperatures are in kelvins or in degrees Rankine (the absolute scale for degrees Fahrenheit).

Volumetric estimation of hydrocarbons in place →
Mean porosity over the net thickness
ϕˉ=∑kϕk hk∑khk\bar\phi = \frac{\sum_k \phi_k\,h_k}{\sum_k h_k}

any units, the same for all intervals

Holds when: We add up only the intervals of net thickness. All the intervals belong to one well (one point on the map).

Volumetric estimation of hydrocarbons in place →
Mean water saturation weighted by pore volume
Sˉw=∑kSw,k ϕk hk∑kϕk hk\bar S_w = \frac{\sum_k S_{w,k}\,\phi_k\,h_k}{\sum_k \phi_k\,h_k}

any units, the same for all intervals

Holds when: We add up the same intervals as for the mean porosity. Saturation is a share of the pore volume, so the weight of an interval is its pore volume, not its thickness.

Volumetric estimation of hydrocarbons in place →
Hydrocarbon pore thickness
H=∑khk ϕk (1−Sw,k)=heff ϕˉ (1−Sˉw)H = \sum_k h_k\,\phi_k\,(1 - S_{w,k}) = h_{\text{eff}}\,\bar\phi\,(1 - \bar S_w)

m (or ft)

Holds when: ϕˉ\bar\phi and Sˉw\bar S_w are averaged by the formulas above.

Volumetric estimation of hydrocarbons in place →
Recoverable volume
Qrec=N⋅RFQ_{\text{rec}} = N\cdot\text{RF}

the same as for NN

Holds when: The RF refers to a particular development system and particular economic conditions. When these change, the recoverable volume changes too. For gas we use the gas recovery factor in the same way.

Volumetric estimation of hydrocarbons in place →
Value of a parameter from a random number
x=F−1(u)x = F^{-1}(u)

uu and FF are dimensionless; xx is in the units of the parameter

Holds when: uu is uniformly distributed: any value between 0 and 1 is equally likely. F(x)F(x) is continuous and increasing, so that each uu corresponds to one value of xx. For independent parameters, the random numbers are independent of one another.

Probabilistic volumetrics by Monte Carlo simulation →
Distribution function of the triangular distribution
F(x)={(x−a)2(b−a)(m−a),a≤x≤m1−(b−x)2(b−a)(b−m),m<x≤bF(x) = \begin{cases} \dfrac{(x - a)^2}{(b - a)(m - a)}, & a \le x \le m \\[10pt] 1 - \dfrac{(b - x)^2}{(b - a)(b - m)}, & m < x \le b \end{cases}

xx, aa, mm and bb are in the units of the parameter; FF is dimensionless

Holds when: The probability density rises linearly from zero at aa to its largest value at mm, and falls linearly to zero at bb. Values smaller than aa and larger than bb are impossible: there F=0F = 0 and F=1F = 1. a<m<ba < m < b. If mm is in the middle, the distribution is symmetrical; otherwise it is skewed.

Probabilistic volumetrics by Monte Carlo simulation →
Drawing from the triangular distribution
x={a+u (b−a)(m−a),u≤F(m)b−(1−u)(b−a)(b−m),u>F(m)F(m)=m−ab−ax = \begin{cases} a + \sqrt{u\,(b - a)(m - a)}, & u \le F(m) \\[4pt] b - \sqrt{(1 - u)(b - a)(b - m)}, & u > F(m) \end{cases} \qquad F(m) = \frac{m - a}{b - a}

xx, aa, mm and bb are in the units of the parameter

Holds when: The same as for the distribution function above. uu is uniformly distributed between 0 and 1.

Probabilistic volumetrics by Monte Carlo simulation →
Error of the mean and of a probability from a simulation
SE=sn,SEp=p (1−p)nSE = \frac{s}{\sqrt{n}}, \qquad SE_p = \sqrt{\frac{p\,(1 - p)}{n}}

SESE is in the units of the result (million m³); SEpSE_p and pp are fractions

Holds when: The trials are independent of one another: a simple random draw. nn is large enough for the mean and the share to follow a distribution close to normal (hundreds of trials or more). The formulas describe the error of the calculation for the given distributions, not the uncertainty of the volume in place.

Probabilistic volumetrics by Monte Carlo simulation →
Mean and variance of the triangular distribution
μ=a+m+b3,σ2=a2+m2+b2−am−ab−mb18\mu = \frac{a + m + b}{3}, \qquad \sigma^2 = \frac{a^2 + m^2 + b^2 - am - ab - mb}{18}

in the units of the parameter; the coefficient of variation Cv=σ/μC_v = \sigma/\mu is dimensionless

Holds when: A triangular distribution with a<m<ba < m < b.

Probabilistic volumetrics by Monte Carlo simulation →
Mean and coefficient of variation of a product of independent quantities
μN=∏jμj,1+Cv2(N)=∏j[1+Cv2(xj)]\mu_N = \prod_j \mu_j, \qquad 1 + C_v^2(N) = \prod_j \left[1 + C_v^2(x_j)\right]

the means are in the units of the quantities; the coefficients of variation are dimensionless

Holds when: The factors are independent. With dependence, the mean of the product is no longer equal to the product of the means. The factor is exactly what enters the product: 1−Sw1 - S_w, not SwS_w, and 1/Boi1/B_{oi}, not BoiB_{oi}. The type of the distributions does not matter: the formulas are exact for any independent factors.

Probabilistic volumetrics by Monte Carlo simulation →
A pair of linked random numbers
z2=ρ z1+1−ρ2  z′,u2=Φ(z2)z_2 = \rho\,z_1 + \sqrt{1 - \rho^2}\;z', \qquad u_2 = \Phi(z_2)

all quantities are dimensionless

Holds when: u1u_1 and u′u' are independent and uniformly distributed between 0 and 1. ρ\rho is the correlation of the normal quantities z1z_1 and z2z_2. The correlation of the parameters themselves is close to it, but does not coincide with it exactly. The distributions of both parameters stay the same: only the way their values combine changes.

Probabilistic volumetrics by Monte Carlo simulation →
Probabilistic meaning of the categories of reserves
Pr⁡(Q≥1P)≥0.9,Pr⁡(Q≥2P)≥0.5,Pr⁡(Q≥3P)≥0.1\Pr(Q \ge \text{1P}) \ge 0.9,\quad \Pr(Q \ge \text{2P}) \ge 0.5,\quad \Pr(Q \ge \text{3P}) \ge 0.1

any units of quantity, the same for all the estimates

Holds when: This is how estimates made by a probabilistic method are checked: 1P, 2P and 3P coincide with the P90, P50 and P10 of the distribution of recoverable quantities (Probabilistic volumetrics by Monte Carlo simulation). In a deterministic estimate, proved reserves are quantities that will be recovered with “reasonable certainty”, that is, with a high degree of confidence. For 2P, the actual recovery is equally likely to be above or below the estimate. The same thresholds apply to contingent resources (1C, 2C, 3C) and prospective resources (1U, 2U, 3U). The probabilities of prospective resources refer to the volume on the condition that the accumulation is discovered.

Classification of reserves and resources →
Chance of commerciality
Pc=Pg⋅PdP_c = P_g\cdot P_d

probabilities, fractions

Holds when: For prospective resources both factors apply. For contingent resources the discovery has already taken place, and Pc=PdP_c = P_d. For reserves no separate chance is stated: the class needs high confidence that the project will be carried out. PdP_d includes the size of the discovery (whether it is enough to pay back; this is how the lesson Geological risk and chance of success estimates the commercial chance of success), and also the market and the approvals.

Classification of reserves and resources →
Starting signature bonus for an area with approved reserves
Bstart=max⁡(0.0004 C+0.0001 Cc2c3; 10,000 MCI)B_\text{start} = \max\left(0.0004\,C + 0.0001\,C_\text{c2c3};\ 10{,}000\ \text{MCI}\right)

tenge (KZT)

Holds when: The value of the reserves is found from the arithmetic mean of the price quotations on the day before the auction conditions are published (or before the protocol of direct negotiations is signed), and at the official exchange rate on the day before the bonus is paid (Article 761 of the 2026 Tax Code). If the area is divided into blocks, then for each block after the three-hundredth the starting amount is increased by 10 MCI. The final bonus is set by the result of the auction or negotiations, and it is not lower than the starting amount.

Classification of reserves and resources →
Reserves replacement ratio
RRR=ΔRrev+ΔRdisc+ΔRdealsQyr\mathrm{RRR} = \frac{\Delta R_\text{rev} + \Delta R_\text{disc} + \Delta R_\text{deals}}{Q_\text{yr}}

the same units for all quantities; reserves of one category (usually 1P or 2P)

Holds when: There is no single mandatory formula: companies state whether deals and revisions caused by prices are included. The ratio of reserves to annual production (R/P) shows for how many years the reserves will last at the present production, if they are not replaced.

Classification of reserves and resources →

Symbols

SymbolMeaning and context
xˉ\bar x
Pr⁡(⋅)\Pr(\cdot)
1P, 2P, 3P\text{1P},\ \text{2P},\ \text{3P}
aa
AA
amax, amina_{max},\ a_{min}
a, m, ba,\ m,\ b
a(α)a(\alpha)
kˉarith\bar k_{\text{arith}}
bb
BgiB_{gi}
BoiB_{oi}
BstartB_\text{start}
CC
C0C_0
C0+CC_0 + C
Cc2c3C_\text{c2c3}
CJC_J
CvC_v
Cv(N)C_v(N)
Cv(xj)C_v(x_j)
C(0)=σ2C(0) = \sigma^2
C(h)C(h)
ΔRdeals\Delta R_\text{deals}
ΔRdisc\Delta R_\text{disc}
F−1(u)F^{-1}(u)
FJF_J
F(m)F(m)
f(x)f(x)
F(x)F(x)
GG
kˉgeom\bar k_{\text{geom}}
GRV\mathrm{GRV}
GRV, NTG, ϕ, Sw\mathrm{GRV},\ \mathrm{NTG},\ \phi,\ S_w
hh
HH
  • the thickness of the layer that the hydrocarbons would form if we gathered them out of the pores into a solid layer; the volume of hydrocarbons at reservoir conditions per unit area, mVolumetric estimation of hydrocarbons in place
h∥, h⊥h_{\parallel},\ h_{\perp}
heffh_{\text{eff}}
hkh_k
h′h'
kˉharm\bar k_{\text{harm}}
II
JJ
k50k_{50}
k84.1k_{84.1}
kik_i
kj, ϕj, hjk_j,\ \phi_j,\ h_j
LcL_c
mm
MCI\text{MCI}
nn
NN
N0N_0
NeffN_{\text{eff}}
NsN_s
N(h)N(h)
NTG\mathrm{NTG}
pp
PP
PcP_c
PdP_d
PgP_g
PsP_s
psc, Tscp_{sc},\ T_{sc}
QQ
QrecQ_{\text{rec}}
QyrQ_\text{yr}
ΔRrev\Delta R_\text{rev}
RF\text{RF}
RRR\mathrm{RRR}
ss
s2s^2
sln⁡s_{\ln}
SwS_w
Sw,kS_{w,k}
s, Ns,\ N
SESE
SEpSE_p
tt
  • Student's coefficient (a table value) for the chosen confidence probability and N−1N - 1 degrees of freedom (the number of independent values left after the mean has been found), —Distributions and sample size
  • Student's coefficient; about 2 at a probability of 95% and more than thirty samples, —Distributions and sample size
TT
t, Nt,\ N
uu
u2u_2
VDPV_{DP}
xx
x50x_{50}
xix_i
xmodex_{\text{mode}}
xpx_p
zz
ZZ
z∗(x0)z^*(x_0)
z1, z′z_1,\ z'
z2z_2
zpz_p
z(xi), z(xi+h)z(x_i),\ z(x_i + h)
α, α0\alpha,\ \alpha_0
γ^(h)\hat\gamma(h)
γ(h)\gamma(h)
γ(xi−x0)\gamma(x_i - x_0)
γ(xi−xj)\gamma(x_i - x_j)
λj\lambda_j
μ\mu
μ, σ\mu,\ \sigma
μj\mu_j
μN\mu_N
ρˉ\bar\rho
  • mean correlation coefficient between pairs of measurements (0 means independent, 1 means that they repeat one another), —Distributions and sample size
ρ\rho
ρ(h)\rho(h)
σ\sigma
σln⁡k\sigma_{\ln k}
σ2\sigma^2
σK2\sigma_K^2
ϕ\phi
Φ\Phi
ϕk\phi_k
ϕk hk\phi_k\,h_k