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

Equations: The Earth and rocks

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

Velocities of compressional and shear waves
Vp=K+43Gρ,Vs=GρV_p = \sqrt{\frac{K + \tfrac{4}{3}G}{\rho}}, \qquad V_s = \sqrt{\frac{G}{\rho}}

SI; if the density is taken in g/cm³ and the velocity in km/s, ρV2\rho V^2 comes out directly in gigapascals

Holds when: A homogeneous, isotropic, elastic medium (the same at every point and in every direction, and springing back after a load). The moduli are dynamic: they apply to the small, fast strains that the wave produces.

Structure and composition of the Earth →
Mountain root under local compensation (the Airy model)
r=ρcρm−ρc hr = \frac{\rho_c}{\rho_m - \rho_c}\,h

any consistent units: rr and hh in the same unit of length, the densities in the same unit

Holds when: The density of the crust is the same everywhere and only its thickness varies (the Airy model). In the Pratt model it is the other way round: the density varies and the depth of the base stays the same. Compensation is local: each column is balanced on its own, and the strength of the lithosphere is ignored. The mantle is the same under the mountains and under the plain.

Structure and composition of the Earth →
Thickness of the sediments that fill a basin (local compensation)
S=ρm−ρwρm−ρs wS = \frac{\rho_m - \rho_w}{\rho_m - \rho_s}\,w

any consistent units

Holds when: The basin is filled up to the former water level, and sea level has not changed. Compensation is local; the lithosphere beneath the basin does not stretch or cool while the basin is filling. The mean density of the sediments is constant (in reality it increases as they compact).

Structure and composition of the Earth →
Temperature in a layer that releases radiogenic heat
T(z)=T0+q0λ z−A2λ z2T(z) = T_0 + \frac{q_0}{\lambda}\,z - \frac{A}{2\lambda}\,z^2

SI: q0q_0 in W/m², AA in W/m³, λ\lambda in W/(m·K), zz in m; temperature in °C

Holds when: Steady one-dimensional transfer of heat by conduction, with no movement of groundwater. λ\lambda and AA are constant in a layer of thickness DD; below it, where A≈0A \approx 0, the gradient is constant and equal to (q0−AD)/λ(q_0 - A D)/\lambda. There has been no recent heating or cooling: rifting, rapid sedimentation and erosion disturb the steady-state profile.

Structure and composition of the Earth →
Distance at which the head wave overtakes the direct wave
xc=2h V2+V1V2−V1x_c = 2h\,\sqrt{\frac{V_2 + V_1}{V_2 - V_1}}

any consistent units

Holds when: The boundary is horizontal and the layers are homogeneous. The source and receivers are at the surface.

Structure and composition of the Earth →
Heat production from the uranium, thorium and potassium content
A=10−5 ρ (9.52 CU+2.56 CTh+3.48 CK)A = 10^{-5}\,\rho\,\left(9.52\,C_U + 2.56\,C_{Th} + 3.48\,C_K\right)

AA in µW/m³, ρ\rho in kg/m³, CUC_U and CThC_{Th} in ppm, CKC_K in % by mass

Holds when: The coefficients are the heat released by the decay of 1 kg of natural uranium (about 95 µW), thorium (about 26 µW) and potassium (about 0.0035 µW), converted to the units of content used. Different tables of decay constants differ by a few per cent. The uranium and thorium decay series are in radioactive equilibrium.

Structure and composition of the Earth →
Oxygen-to-silicon ratio in a silicate structure
NONSi=4−s2\frac{N_{\text{O}}}{N_{\text{Si}}} = 4 - \frac{s}{2}

dimensionless

Holds when: All the tetrahedra are built around silicon. If aluminium replaces part of the silicon, the denominator is the sum of the silicon and aluminium atoms in the tetrahedra.

Minerals and rocks →
Density of a fragment by hydrostatic weighing
ρ=ρl WaWa−Wl\rho = \rho_l\,\frac{W_a}{W_a - W_l}

scale readings in the same unit, densities in the same unit

Holds when: The fragment is a single mineral, without pores and cracks that the liquid enters (otherwise the density includes the voids). The liquid does not dissolve or react with the fragment, and air bubbles have been removed. The buoyancy of air is neglected (the correction is less than 0.1%).

Minerals and rocks →
Calcite and dolomite content from the volume of carbon dioxide
wcc=V1k m,wdol=V21.086 k mw_{cc} = \frac{V_1}{k\,m}, \qquad w_{dol} = \frac{V_2}{1.086\,k\,m}

mL and g; find kk by calibration under the same conditions

Holds when: There is an excess of hydrochloric acid. Calcite reacts almost completely in the first minute and dolomite over tens of minutes, and the time boundary between them is arbitrary. Siderite and other slowly reacting carbonates are counted as “dolomite”. Part of the gas dissolves in the acid; calibration on pure calcite under the same conditions allows for this.

Minerals and rocks →
Spreading rate from a dated magnetic anomaly
v=xt,V=v1+v2v = \frac{x}{t}, \qquad V = v_1 + v_2

kilometres and millions of years: 1 km/million years = 1 mm/yr

Holds when: The sea floor forms at a narrow ridge axis, and its magnetisation does not change after cooling. The speed is the mean over the time tt; if it has changed, different anomalies give different values. The distance is measured along the direction of opening, parallel to the transform faults. An oblique profile across the stripes distorts the speed. If the ridge axis is perpendicular to the opening, the oblique profile makes the speed too high.

Plate tectonics and sedimentary basins →
Initial (syn-rift) subsidence when the lithosphere is stretched
Si=a[(ρm−ρc) tca(1−αVT1tc2a)−ρmαVT12]ρm (1−αVT1)−ρw(1−1β)S_i = \frac{a\left[(\rho_m - \rho_c)\,\frac{t_c}{a}\left(1 - \frac{\alpha_V T_1 t_c}{2a}\right) - \frac{\rho_m \alpha_V T_1}{2}\right]}{\rho_m\,(1 - \alpha_V T_1) - \rho_w}\left(1 - \frac{1}{\beta}\right)

SI; with aa and tct_c in kilometres, SiS_i comes out in kilometres

Holds when: The stretching is instant and uniform. The crust and the mantle part of the lithosphere thin by the same factor, and temperature in the lithosphere still rises linearly with depth. Isostasy is local; the strength of the lithosphere is ignored. Rock density falls with temperature as ρ(T)=ρ0(1−αVT)\rho(T) = \rho_0(1 - \alpha_V T). The depression is filled with water, and the initial surface is at sea level.

Plate tectonics and sedimentary basins →
Thermal (post-rift) subsidence
St(t)=E0 βπsin⁡πβ(1−e−t/τ),E0=4 a ρmαVT1π2(ρm−ρw),τ=a2π2κS_t(t) = E_0\,\frac{\beta}{\pi}\sin\frac{\pi}{\beta}\left(1 - e^{-t/\tau}\right),\qquad E_0 = \frac{4\,a\,\rho_m \alpha_V T_1}{\pi^2(\rho_m - \rho_w)},\qquad \tau = \frac{a^2}{\pi^2 \kappa}

SI (τ\tau is converted from seconds to millions of years: 1 million years ≈ 3.156⋅10133.156\cdot10^{13} s)

Holds when: The lithosphere cools by vertical conduction; the basin is wider than about 100 km, and no heat escapes sideways. The series that describes the cooling is cut to its first, slowest term. For β\beta from 1.3 to 3, the limiting thermal subsidence differs from the sum of the full series by a few per cent. For β\beta near 1.2 it is about 7% too low. In the first 20 million years after the stretching the difference is larger, up to 10–25%, but the thermal subsidence itself is then still small. Radiogenic heat is not included, nor is the fact that sediments of low thermal conductivity “insulate” the lithosphere and slow its cooling.

Plate tectonics and sedimentary basins →
Flexural rigidity and flexural parameter of the lithosphere
D=E Te312 (1−ν2),α=[4D(ρm−ρi) g]1/4D = \frac{E\,T_e^3}{12\,(1 - \nu^2)}, \qquad \alpha = \left[\frac{4D}{(\rho_m - \rho_i)\,g}\right]^{1/4}

SI

Holds when: A thin elastic plate of constant thickness on a fluid base (over the time it takes a basin to form, the asthenosphere behaves as a fluid). The whole basin is filled with material of a single density ρi\rho_i; the plate is pushed up by the difference ρm−ρi\rho_m - \rho_i.

Plate tectonics and sedimentary basins →
Deflection of a broken elastic plate under a load at its end
w(x)=w0 e−x/αcos⁡xα,x0=π2 α,xb=3π4 α,wb≈−0.067 w0w(x) = w_0\,e^{-x/\alpha}\cos\frac{x}{\alpha}, \qquad x_0 = \frac{\pi}{2}\,\alpha, \qquad x_b = \frac{3\pi}{4}\,\alpha, \qquad w_b \approx -0.067\,w_0

any consistent units: xx and α\alpha in the same unit of length

Holds when: The plate is broken beneath the thrust front, and the load is a line load (it acts along a line) at the end of the plate. This is the usual way to model foreland basins in front of mountain belts. If the plate is taken as continuous, the edge of the basin is farther (x0=3πα/4x_0 = 3\pi\alpha/4), the crest is at xb=παx_b = \pi\alpha, and the forebulge is lower — 4.3% of w0w_0.

Plate tectonics and sedimentary basins →
Isostatic balance of the original and the stretched column
M0=ρwSi+M0β+ρm(1−αVT1)(a−aβ−Si)M_0 = \rho_w S_i + \frac{M_0}{\beta} + \rho_m(1 - \alpha_V T_1)\left(a - \frac{a}{\beta} - S_i\right)

SI

Holds when: Local isostasy: at depth aa the pressure beneath the original and the stretched columns is the same.

Plate tectonics and sedimentary basins →
Speed of a point of a plate rotating about a pole
v=ω R sin⁡Δv = \omega\,R\,\sin\Delta

km and millions of years

Holds when: The plates are rigid; the Earth is a sphere.

Plate tectonics and sedimentary basins →
The law of radioactive decay, and age from the accumulated daughter isotope
P=P0 e−λt,t1/2=ln⁡2λ,t=1λln⁡(1+D∗P)P = P_0\,e^{-\lambda t}, \qquad t_{1/2} = \frac{\ln 2}{\lambda}, \qquad t = \frac{1}{\lambda}\ln\left(1 + \frac{D^*}{P}\right)

λ\lambda and tt in reciprocal units (yr⁻¹ and years); D∗/PD^*/P is a ratio of numbers of atoms; ln is the natural logarithm

Holds when: The system is closed: after the mineral formed, neither parent nor daughter isotope entered or left it. There was no daughter isotope at the start, or its initial amount is known (otherwise the isochron method is used — see the advanced section). Every decay gives an atom of exactly this daughter isotope; for ⁴⁰K only about 10% of decays turn into ⁴⁰Ar (the rest give ⁴⁰Ca), so for potassium–argon ages D∗/PD^*/P is divided by this fraction. The decay constant is known exactly; the age is counted from the moment when the mineral cooled enough for the daughter isotope to stop escaping from it.

Stratigraphy and geologic time →
Rate at which space for sediment is created
A˙=S˙+E˙\dot A = \dot S + \dot E

m/Myr (metres per million years, the same as mm per thousand years)

Holds when: Compaction of previously deposited sediment is included in the subsidence. The bending of the lithosphere under the weight of the sediment and the water is not separated from the tectonic subsidence (isostasy).

Stratigraphy and geologic time →
The isochron of the rubidium–strontium system
87Sr86Sr=(87Sr86Sr)0+(eλt−1)87Rb86Sr\frac{{}^{87}\mathrm{Sr}}{{}^{86}\mathrm{Sr}} = \left(\frac{{}^{87}\mathrm{Sr}}{{}^{86}\mathrm{Sr}}\right)_0 + \left(e^{\lambda t} - 1\right)\frac{{}^{87}\mathrm{Rb}}{{}^{86}\mathrm{Sr}}

ratios of numbers of atoms, dimensionless

Holds when: The samples are contemporaneous and had the same initial strontium isotope ratio. The system remained closed after it formed.

Stratigraphy and geologic time →

Symbols

SymbolMeaning and context
A˙\dot A
  • rate of relative sea-level rise, that is, of the creation of space for sediment; a negative value means a relative fall of sea level, m/MyrStratigraphy and geologic time
S˙\dot S
E˙\dot E
(87Sr86Sr)0\left(\frac{{}^{87}\mathrm{Sr}}{{}^{86}\mathrm{Sr}}\right)_0
87Rb86Sr\frac{{}^{87}\mathrm{Rb}}{{}^{86}\mathrm{Sr}}
87Sr86Sr\frac{{}^{87}\mathrm{Sr}}{{}^{86}\mathrm{Sr}}
aa
AA
a, ρm, ρw, αV, T1, βa,\ \rho_m,\ \rho_w,\ \alpha_V,\ T_1,\ \beta
CKC_K
CU, CThC_U,\ C_{Th}
DD
D∗D^*
  • number of atoms of the radiogenic daughter isotope accumulated over the time tt: D∗=P0−PD^* = P_0 - P, atomsStratigraphy and geologic time
EE
E0E_0
gg
GG
hh
kk
  • volume of gas from 1 g of pure calcite, found by calibrating the apparatus; about 240 mL/g at 20 °C and atmospheric pressure, mL/gMinerals and rocks
KK
mm
M0M_0
NO/NSiN_{\text{O}}/N_{\text{Si}}
P0, PP_0,\ P
q0q_0
rr
RR
ss
  • how many corners, on average, a tetrahedron shares with neighbouring tetrahedra (from 0 to 4), —Minerals and rocks
SS
SiS_i
Si, ρw, ρm, αV, T1, a, βS_i,\ \rho_w,\ \rho_m,\ \alpha_V,\ T_1,\ a,\ \beta
St(t)S_t(t)
tt
T0T_0
T1T_1
t1/2t_{1/2}
tct_c
TeT_e
T(z)T(z)
vv
VV
V1V_1
V1, V2V_1,\ V_2
V2V_2
Vp, VsV_p,\ V_s
ww
w0w_0
WaW_a
wcc, wdolw_{cc},\ w_{dol}
WlW_l
w(x)w(x)
xx
x0x_0
xb, wbx_b,\ w_b
xcx_c
zz
α\alpha
αV\alpha_V
β\beta
Δ\Delta
κ\kappa
λ\lambda
ν\nu
ρ\rho
ρc\rho_c
ρi\rho_i
ρl\rho_l
ρm\rho_m
ρm, ρc\rho_m,\ \rho_c
ρm, ρw, ρs\rho_m,\ \rho_w,\ \rho_s
ρw\rho_w
τ\tau
ω\omega