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ˉ=N1i=1∑Nxi
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.
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.
s is in the units of the measured quantity, s2 is in their squares
Holds when: The measurements belong to one population. The denominator is N−1, not N. The deviations are measured from the mean of the same sample, so on average they are slightly smaller than from the true mean. For N above a few dozen, the difference does not matter.
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.
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).
x, μ and σ 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 μ is greater than about 3σ. The distribution is unimodal (it has one peak) and symmetrical.
Median, mean and mode of the log-normal distribution
x50=eμ,xˉ=eμ+σ2/2,xmode=eμ−σ2
x is in the units of the quantity (for example, mD); μ and σ refer to lnx and are dimensionless
Holds when: The logarithm of the quantity is normally distributed. We calculate μ and σ from natural logarithms (base e). For decimal logarithms, μln=2.303μlg and σln=2.303σlg.
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.
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.
Number of samples needed for a given precision of the mean
N0=(P100tCv)2
P 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. Cv is known at least approximately, from a trial series of measurements or by analogy with similar rocks.
Rule for the number of samples: precision of the mean ±20% at 95% probability
N0=(10Cv)2,Ps=Ns200Cv
Ps is in per cent of the mean
Holds when:t=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=2. For very heterogeneous rocks (Cv noticeably greater than 1, for example complex carbonates), a tolerance of ±20% is often unattainable. Then we accept ±50%, and N0=(4Cv)2.
Holds when: Permeability is log-normally distributed. For other distributions the coefficient can be misleading. We take k50 and k84.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.1 at the 16% mark (z=−1).
Standard error of the mean for dependent measurements
SE=Ns1+(N−1)ρˉ,Neff=1+(N−1)ρˉN
SE 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.
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 (h 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.
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>0, and the variogram itself is zero at zero. The exponential and Gaussian models are often written with the parameter a/3 or a/3. When we move parameters between programs, we check which range is meant.
in squared units of the property; ρ 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.
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.
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.
Holds when: The weights and μ 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.
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 ϕ and Sw refer to the same volume of reservoir as the NTG. Boi 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.
Holds when: The pores of the reservoir are filled with gas and connate water (a gas accumulation or a gas cap). Bgi 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⋅Rsi, where Rsi is the initial solution gas–oil ratio of the oil (how much gas is dissolved in each cubic metre of oil).
any consistent units: pressures in the same units, temperatures absolute (K or °R)
Holds when: The real-gas equation of state pV=ZnRT (it links the pressure, volume and temperature of the gas); at standard conditions the gas is taken as ideal (Z=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).
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.
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.
u and F are dimensionless; x is in the units of the parameter
Holds when:u is uniformly distributed: any value between 0 and 1 is equally likely. F(x) is continuous and increasing, so that each u corresponds to one value of x. For independent parameters, the random numbers are independent of one another.
x, a, m and b are in the units of the parameter; F is dimensionless
Holds when: The probability density rises linearly from zero at a to its largest value at m, and falls linearly to zero at b. Values smaller than a and larger than b are impossible: there F=0 and F=1. a<m<b. If m is in the middle, the distribution is symmetrical; otherwise it is skewed.
Error of the mean and of a probability from a simulation
SE=ns,SEp=np(1−p)
SE is in the units of the result (million m³); SEp and p are fractions
Holds when: The trials are independent of one another: a simple random draw. n 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.
Mean and coefficient of variation of a product of independent quantities
μN=j∏μj,1+Cv2(N)=j∏[1+Cv2(xj)]
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−Sw, not Sw, and 1/Boi, not Boi. The type of the distributions does not matter: the formulas are exact for any independent factors.
Holds when:u1 and u′ are independent and uniformly distributed between 0 and 1. ρ is the correlation of the normal quantities z1 and z2. 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 meaning of the categories of reserves
Pr(Q≥1P)≥0.9,Pr(Q≥2P)≥0.5,Pr(Q≥3P)≥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.
Holds when: For prospective resources both factors apply. For contingent resources the discovery has already taken place, and Pc=Pd. For reserves no separate chance is stated: the class needs high confidence that the project will be carried out. Pd 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.
Starting signature bonus for an area with approved reserves
Bstart=max(0.0004C+0.0001Cc2c3;10,000MCI)
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.
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.
range: for the spherical model, the distance at which it reaches the sill; for the exponential and Gaussian models, the distance at which it reaches 95% of the partial sill, mSpatial correlation and the variogram
starting amount of the signature bonus under a production contract, or an exploration and production contract, for hydrocarbons, tengeClassification of reserves and resources
value of the reserves of category C2 (approved, or taken into account from the operational estimate) and of the resources C3, tengeClassification of reserves and resources
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
monthly calculation index on the date of publishing the auction conditions or signing the protocol of direct negotiations; in 2026 it is 4,325 tenge, tengeClassification of reserves and resources
tolerance — the permissible deviation of the mean of the samples from the true mean, % of the meanDistributions and sample size
Pc
chance of commerciality: the probability that the project will be approved and will reach commercial production, fractionClassification of reserves and resources
the quantity that will actually be recovered by the project from the date of the estimate to the end of production, t, m³ or bblClassification of reserves and resources
Student's coefficient (a table value) for the chosen confidence probability and N−1 degrees of freedom (the number of independent values left after the mean has been found), —Distributions and sample size