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

Equations: Petroleum system and traps

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

Geological chance of success of a prospect
Pg=∏i=1kPi=P1 P2⋯PkP_g = \prod_{i=1}^{k} P_i = P_1\,P_2\cdots P_k

dimensionless

Holds when: The elements are independent. If they depend on each other (for example, a reservoir and a seal from one depositional environment), we must not multiply them. We use conditional probabilities instead (the chance of one element given that another is present). Each PiP_i is a subjective estimate. We calibrate it against drilling statistics in similar plays. Geological success is not the same as economic success: what is found may not be enough to pay for development.

Petroleum system and the play concept →
Common play risk and local prospect risk
Pg=Pplay⋅PlocP_g = P_\text{play}\cdot P_\text{loc}

dimensionless

Holds when: The common and local elements are split so that each element is counted once. PlocP_\text{loc} is a conditional probability: it is estimated on the assumption that the play works.

Petroleum system and the play concept →
Entry pressure of a pore channel
pd=2σcos⁡θrp_d = \frac{2\sigma\cos\theta}{r}

SI

Holds when: The channel is modelled as a cylinder. For a real rock, rr is the effective radius of the largest channels that form a connected path through the seal. The rock is water-wet (θ<90∘\theta < 90^\circ). In an oil-wet rock (one whose surface prefers oil) there is no capillary barrier for oil (Wettability and interfacial phenomena). Take σ\sigma and θ\theta at reservoir pressure and temperature.

Reservoir and seal →
Excess pressure of rising hydrocarbons (buoyancy)
Δpb=(ρw−ρh) g h\Delta p_b = (\rho_w - \rho_h)\,g\,h

SI

Holds when: The water in the layer is at rest (hydrostatics). When formation waters move, the contact tilts (Traps and timing). The column of hydrocarbons is continuous and connected to the water-bearing part of the layer. Densities are constant over the height of the accumulation.

Reservoir and seal →
Height of accumulation that the seal will hold
Hmax=pd,s−pd,r(ρw−ρh) g=2σcos⁡θ(1rs−1rr)(ρw−ρh) gH_{max} = \frac{p_{d,s} - p_{d,r}}{(\rho_w - \rho_h)\,g} = \frac{2\sigma\cos\theta\left(\dfrac{1}{r_s} - \dfrac{1}{r_r}\right)}{(\rho_w - \rho_h)\,g}

SI

Holds when: The water is at rest. The seal is continuous, with no fractures or conducting faults. The height is measured from the contact. Measured from the free water level, the height held is larger by pd,r/[(ρw−ρh) g]p_{d,r}/[(\rho_w - \rho_h)\,g]. The effective radius is the radius of the largest connected channels of the seal. Leakage follows the easiest path, not the average pore size. All properties are taken at reservoir conditions.

Reservoir and seal →
Converting entry pressure from mercury to reservoir fluids
pd,res=pd,Hg (σcos⁡θ)res∣σcos⁡θ∣Hg,r=2 ∣σcos⁡θ∣Hgpd,Hgp_{d,\text{res}} = p_{d,\text{Hg}}\,\frac{(\sigma\cos\theta)_{\text{res}}}{|\sigma\cos\theta|_{\text{Hg}}},\qquad r = \frac{2\,|\sigma\cos\theta|_{\text{Hg}}}{p_{d,\text{Hg}}}

SI (or any consistent units of pressure)

Holds when: The channels are modelled as cylinders, and the geometry of the pores is the same for mercury and for reservoir fluids. The contact angle of mercury is measured through the mercury, so we take the modulus of the cosine. The sample is representative. Without the load of the surrounding rock, its channels may be wider than in the reservoir.

Reservoir and seal →
Height of accumulation limited by fracturing of the seal
Hfrac=pfrac−pw(ρw−ρh) gH_{frac} = \frac{p_{frac} - p_w}{(\rho_w - \rho_h)\,g}

SI

Holds when: The pressures are taken at the depth of the top of the accumulation. pwp_w is the water pressure continued along the gradient of the water-bearing part of the layer. The seal opens when the pressure in it exceeds the fracture pressure. We ignore its tensile strength (how well the rock resists being pulled apart). The estimate is rough: the fracture pressure of the seal is known with a considerable error.

Reservoir and seal →
Hydrogen and oxygen indices
HI=100 S2TOC,OI=100 S3TOCHI = \frac{100\,S_2}{TOC},\qquad OI = \frac{100\,S_3}{TOC}

as shown; TOC/100TOC/100 is grams of organic carbon per gram of rock, so 100 S2/TOC100\,S_2/TOC is milligrams of HC per gram of carbon

Holds when: The sample is not contaminated by oil-based drilling fluid, migrated oil or weathering. To find the original kerogen type, you need an immature or only slightly mature sample, because HI falls as the rock matures.

Source rocks, maturation and migration →
Production index
PI=S1S1+S2PI = \frac{S_1}{S_1 + S_2}

dimensionless

Holds when: Some of the hydrocarbons that formed have already left the rock, so PI underestimates the share of kerogen that has converted. Migrated oil and drilling-fluid additives make S1S_1 and PI too high.

Source rocks, maturation and migration →
The Arrhenius law for the rate of kerogen breakdown
k=A e−E/(RT)k = A\,\e^{-E/(RT)}

SI; temperature only in kelvin (T=t+273.15T = t + 273.15)

Holds when: Treating the breakdown as a first-order reaction with one activation energy is a simplification: real kerogen is a mixture of structures with a range of activation energies. Laboratory pyrolysis at several heating rates gives AA and EE, and we then apply them to geological time.

Source rocks, maturation and migration →
Vitrinite reflectance profile
lg⁡Ro=a+b z\lg R_o = a + b\,z

RoR_o in %, zz in m; lg⁡\lg is the base-10 logarithm

Holds when: An empirical relation (based on data, not on theory) for a section with continuous burial and an unchanged geothermal regime (the same way temperature rises with depth). At erosion surfaces the profile jumps in a step. Where the heat flow has changed, the slope changes. The measurements have been screened, because caved cuttings (rock that fell into the hole from higher up), reworked vitrinite (vitrinite redeposited from older rocks) and underestimated values in hydrogen-rich rocks produce outliers (points far from the trend).

Source rocks, maturation and migration →
Mass of hydrocarbons generated in the kitchen
M=A h ρr S2 0 TRM = A\,h\,\rho_r\,S_2^{\,0}\,\mathrm{TR}

AA in m², hh in m, ρr\rho_r in t/m³ — their product gives the mass of rock in tonnes; S2 0S_2^{\,0} in kg/t, so MM comes out in kilograms

Holds when: Thickness, TOC and kerogen type are averaged over the kitchen. In practice they are mapped. We estimate the transformation ratio from the fall in HI or from a heating model. The result is the mass generated, not the mass expelled, and certainly not the mass accumulated.

Source rocks, maturation and migration →
Shift of the window temperature at a different heating rate
ΔT≈R T2E ln⁡β2β1\Delta T \approx \frac{R\,T^2}{E}\,\ln\frac{\beta_2}{\beta_1}

SI; ln⁡\ln is the natural logarithm

Holds when: A constant heating rate and a single activation energy. ΔT\Delta T is small compared with TT (we ignore the change of the factor T2T^2).

Source rocks, maturation and migration →
Tilt of the oil–water contact under moving formation water
 dzc dx=ρwρw−ρo⋅ dh dx\frac{\dd z_c}{\dd x} = \frac{\rho_w}{\rho_w - \rho_o}\cdot\frac{\dd h}{\dd x}

any consistent system; slopes are conveniently expressed in metres per kilometre

Holds when: The water flow is steady (it does not change with time), and the oil in the trap does not move. We ignore capillary effects, and the densities are constant. For a gas–water contact, use the gas density ρg\rho_g instead of ρo\rho_o.

Traps and timing →
Layer thickness at a different porosity (conservation of grain volume)
h1 (1−ϕ1)=h2 (1−ϕ2)h_1\,(1 - \phi_1) = h_2\,(1 - \phi_2)

thicknesses in any identical units, porosity as a fraction

Holds when: The grain volume is conserved: cement is neither added nor dissolved. The porosity is an average over the layer. We take it from the porosity–depth relation for the rock concerned.

Traps and timing →
Geological chance of success through conditional probabilities
Pg=P1⋅P2∣1⋅P3∣1, 2⋯Pk∣1,…,k−1P_g = P_1\cdot P_{2|1}\cdot P_{3|1,\,2}\cdots P_{k|1,\ldots,k-1}

dimensionless

Holds when: The elements are defined so that success happens when all of them have worked, and only then. Each uncertainty is counted in exactly one factor. The order of the elements does not matter, but we assess each following probability given that the elements before it have worked. For independent elements P2∣1=P2P_{2|1} = P_2 and so on, and the formula becomes the simple product Pg=∏iPiP_g = \prod_i P_i.

Geological risk and chance of success →
Probability that all prospects are dry
P0ind=∏i=1n(1−Pg,i),P0play=(1−Pplay)+Pplay(1−Ploc)nP_0^{\text{ind}} = \prod_{i=1}^{n} \left(1 - P_{g,i}\right), \qquad P_0^{\text{play}} = \left(1 - P_\text{play}\right) + P_\text{play}\left(1 - P_\text{loc}\right)^{n}

dimensionless

Holds when: In the first formula the outcomes of the prospects are independent. In the second the prospects are linked only by the common play risk, and their local elements are independent. The probability of at least one discovery equals 1−P01 - P_0.

Geological risk and chance of success →
Commercial chance of success
Pc=Pg⋅P(V≥Vmin)P_c = P_g\cdot P(V \ge V_\text{min})

dimensionless; VV and VminV_\text{min} in the same units

Holds when: The volume distribution is built given geological success and starts at the smallest volume that counts as success. VminV_\text{min} depends on price, costs, taxes and infrastructure: for the same trap it differs between companies and between years. The chance of success and the volume distribution relate to the same accumulation under one definition of success.

Geological risk and chance of success →
Risked volume
Vˉrisked=Pg Vˉ\bar V_\text{risked} = P_g\,\bar V

volume units

Holds when: If the well fails, the volume is zero, so the mean over all outcomes is PgVˉ+(1−Pg)⋅0P_g\bar V + (1 - P_g)\cdot 0. We use the mean, not the P50, because the mean volumes of prospects add up, but the P50 values do not.

Geological risk and chance of success →
Expected number of discoveries and its scatter
kˉ=∑i=1nPg,i,σk=∑i=1nPg,i(1−Pg,i)\bar k = \sum_{i=1}^{n} P_{g,i}, \qquad \sigma_k = \sqrt{\sum_{i=1}^{n} P_{g,i}\left(1 - P_{g,i}\right)}

dimensionless (a number of wells)

Holds when: kˉ\bar k holds whatever the dependence between the prospects. The formula for σk\sigma_k is for independent prospects. A common play risk increases the scatter. The probabilities were recorded before drilling and were not corrected afterwards.

Geological risk and chance of success →
Probability that the play works after a dry hole
Pplay′=Pplay(1−Ploc)1−Pplay PlocP_\text{play}' = \frac{P_\text{play}\left(1 - P_\text{loc}\right)}{1 - P_\text{play}\,P_\text{loc}}

dimensionless

Holds when: The well was drilled on a prospect of the same segment, and the cause of the failure has not been established. The local elements of different prospects are independent. If the well data show that exactly a local element failed (for example, there is no closure), the common risk hardly changes. The formula does not suit that case.

Geological risk and chance of success →

Symbols

SymbolMeaning and context
ΔT\Delta T
Vˉ\bar V
kˉ\bar k
(σcos⁡θ)res(\sigma\cos\theta)_{\text{res}}
∣σcos⁡θ∣Hg|\sigma\cos\theta|_{\text{Hg}}
  • for mercury in an evacuated sample we usually take σ≈0.485\sigma \approx 0.485 N/m and θ≈140∘\theta \approx 140^\circ (measured through the mercury), that is about 0.372 N/m, N/mReservoir and seal
AA
a, ba,\ b
Δpb\Delta p_b
  • excess of the hydrocarbon pressure over the water pressure at the same depth, PaReservoir and seal
EE
gg
hh
h1, h2h_1,\ h_2
HfracH_{frac}
  • height above the free water level at which the pressure of the hydrocarbons at the top reaches the fracture pressure of the seal, mReservoir and seal
HmaxH_{max}
  • greatest height of the column of hydrocarbons from the contact (oil–water, OWC, or gas–water, GWC) to the top under the seal, mReservoir and seal
HIHI
kk
MM
nn
OIOI
P0indP_0^{\text{ind}}
P0playP_0^{\text{play}}
P1P_1
P2∣1P_{2|1}
PcP_c
pdp_d
  • entry pressure: the excess of the hydrocarbon pressure over the water pressure at which the boundary between the hydrocarbon and the water passes through the narrowest point of the channel, PaReservoir and seal
pd,Hgp_{d,\text{Hg}}
pd,resp_{d,\text{res}}
pd,s, pd,rp_{d,s},\ p_{d,r}
pfracp_{frac}
  • fracture pressure of the seal at the depth of the top, close to the minimum horizontal stress, PaReservoir and seal
PgP_g
Pg,iP_{g,i}
PiP_i
PlocP_\text{loc}
PplayP_\text{play}
Pplay, PlocP_\text{play},\ P_\text{loc}
Pplay′P_\text{play}'
pwp_w
P(V≥Vmin)P(V \ge V_\text{min})
PIPI
rr
RR
RoR_o
rs, rrr_s,\ r_r
  • effective radii of the pore channels of the seal and of the reservoir, mReservoir and seal
Vˉrisked\bar V_\text{risked}
S1S_1
S2S_2
S2 0S_2^{\,0}
S3S_3
TT
TOCTOC
TR\mathrm{TR}
VV
VminV_\text{min}
xx
zz
zcz_c
β1, β2\beta_1,\ \beta_2
θ\theta
  • contact angle (the angle at which the boundary between the fluids meets the rock surface), measured through the water; close to zero for a water-wet rock (one whose surface prefers water), °Reservoir and seal
ρr\rho_r
ρw, ρh\rho_w,\ \rho_h
  • densities of formation water and of the hydrocarbon (oil ρo\rho_o or gas ρg\rho_g) at reservoir conditions, kg/m³Reservoir and seal
ρw, ρh, g\rho_w,\ \rho_h,\ g
ρw, ρo\rho_w,\ \rho_o
σ\sigma
  • interfacial tension between the hydrocarbon and water (how strongly the boundary between the two fluids pulls on itself), N/mReservoir and seal
σ, θ\sigma,\ \theta
  • interfacial tension and contact angle, as in the formula for entry pressure, N/m, °Reservoir and seal
σk\sigma_k
ϕ1, ϕ2\phi_1,\ \phi_2