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

Equations: Structural geology

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

Stresses on an inclined plane (two-dimensional case)
σn=σ1+σ32+σ1−σ32cos⁡2θ,τ=σ1−σ32sin⁡2θ\sigma_n = \frac{\sigma_1 + \sigma_3}{2} + \frac{\sigma_1 - \sigma_3}{2}\cos 2\theta, \qquad \tau = \frac{\sigma_1 - \sigma_3}{2}\sin 2\theta

any identical units of stress; the angle in degrees or radians

Holds when: The plane contains the direction of the intermediate principal stress σ2\sigma_2, so it does not enter the formulas. The stresses are uniform near the plane. The formulas follow from the balance of forces on a triangular prism (advanced section). For effective stresses, subtract the pore pressure from σ1\sigma_1, σ3\sigma_3 and σn\sigma_n. The shear stress τ\tau does not depend on it.

Stress, strain and structural analysis →
Extension, stretch and shear strain
e=l−l0l0,S=ll0=1+e,γ=tan ψe = \frac{l - l_0}{l_0}, \qquad S = \frac{l}{l_0} = 1 + e, \qquad \gamma = \mathrm{tan}\,\psi

dimensionless quantities; lengths in the same units

Holds when: The quantities describe the finite (total, accumulated) strain, not the path by which it developed. In homogeneous strain (the same everywhere) they are the same throughout the volume. In inhomogeneous strain they refer to a small region.

Stress, strain and structural analysis →
Apparent dip on a vertical section
tan δ′=tan δ⋅sin⁡β=tan δ⋅cos⁡ω\mathrm{tan}\,\delta' = \mathrm{tan}\,\delta\cdot\sin\beta = \mathrm{tan}\,\delta\cdot\cos\omega

angles in degrees; tan⁡\tan is the tangent (written tg in Russian-language notation)

Holds when: The bed is flat (planar) within the area. For a curved surface, we apply the formula to the tangent plane (the flat plane that just touches the surface). The section is vertical. For a tilted section, for example along the hole of a deviated well, we need a vector calculation (advanced section).

Stress, strain and structural analysis →
True thickness of a bed from the vertical thickness
h=hvcos⁡δh = h_\text{v}\cos\delta

any identical units of length

Holds when: The top and base are parallel, and the dip is constant over the thickness of the bed. The well is vertical. In a deviated well, the length of the hole in the bed is related to the true thickness through the angle between the hole and the normal to the bed (advanced section).

Stress, strain and structural analysis →
Dip angle from structure contours
tan δ=Δzw\mathrm{tan}\,\delta = \frac{\Delta z}{w}

Δz\Delta z and ww in the same units; ww is the distance on the ground, that is, the distance on the map multiplied by the second number of the scale (25,000 in 1:25,000)

Holds when: The surface between neighbouring contours is taken to be planar. We take the distance at right angles to the contours. In any other direction, we would get an apparent dip.

Stress, strain and structural analysis →
Angle on a section with vertical exaggeration
tan δV=V⋅tan δ′\mathrm{tan}\,\delta_V = V\cdot\mathrm{tan}\,\delta'

angles in degrees

Holds when: The exaggeration is the same throughout the section. On a seismic section in time (where the vertical axis is time, not depth), the vertical scale also changes with the wave velocity, so angles are distorted unevenly (Seismic interpretation).

Stress, strain and structural analysis →
Angle between the hole and the normal to the bed, and the length of hole in the bed
cos⁡γ=∣sin⁡δ sin⁡θ cos⁡(a−φ)−cos⁡δ cos⁡θ∣,h=Δl cos⁡γ\cos\gamma = \left|\sin\delta\,\sin\theta\,\cos(a - \varphi) - \cos\delta\,\cos\theta\right|, \qquad h = \Delta l\,\cos\gamma

angles in degrees; lengths in the same units

Holds when: The bed is flat, and its top and base are parallel. The hole is straight over the interval of the bed. The azimuths of the bed and of the hole are measured from the same north.

Stress, strain and structural analysis →
Equal-area projection of a line onto the lower hemisphere
r=R2 sin⁡90°−ρ2,x=rsin⁡φ,y=rcos⁡φr = R\sqrt{2}\,\sin\frac{90° - \rho}{2}, \qquad x = r\sin\varphi,\quad y = r\cos\varphi

rr, RR, xx, yy in the same units of length; angles in degrees

Holds when: We use the lower-hemisphere projection, which is the standard one in structural geology. We show a plane either by an arc (the projections of all the lines lying in it) or by a single point, its pole. The pole is the projection of the normal, which plunges at an angle 90°−δ90° - \delta in the direction opposite to the dip.

Stress, strain and structural analysis →
Components of the slip on a fault
Dd=Dsin⁡r,Ds=Dcos⁡r,t=Ddsin⁡δ,w=Ddcos⁡δD_d = D\sin r,\quad D_s = D\cos r,\quad t = D_d\sin\delta,\quad w = D_d\cos\delta

any equal units of length; angles in degrees

Holds when: The fault plane is flat within the area. The values refer to the given point of the fault, since the slip changes along it. tt and ww are parts of the slip along the dip. The component DsD_s is horizontal and runs along the fault.

Faults →
Vertical and stratigraphic separation of an inclined bed
s=t±w tan δb,sn=scos⁡δb=Ddsin⁡(δ±δb)s = t \pm w\,\mathrm{tan}\,\delta_b,\qquad s_n = s\cos\delta_b = D_d\sin(\delta \pm \delta_b)

equal units of length; angles in degrees

Holds when: The sign “+” is for beds that dip towards the fault plane, and “−” for beds that dip the same way. The slip is along the dip. The beds are planar (not curved) and run parallel to the strike of the fault. Otherwise we take the apparent dip of the beds (the dip seen in a section that is not at right angles to the strike) in the section across the fault. The strike-slip part of the slip also changes ss. The well is vertical. In an inclined well, we convert the length of the interval along the hole to thickness at right angles to the beds, using the angle between the hole and the normal to the beds (inclined wells). The difference of vertical depths equals ss only for horizontal beds or for a hole along the strike. The stratigraphic separation at the crossing point does not depend on the well path.

Faults →
Power-law distribution of faults by throw
N(t)=N1(t1t)c,c=ln⁡(N1/N2)ln⁡(t2/t1)N(t) = N_1\left(\frac{t_1}{t}\right)^{c},\qquad c = \frac{\ln(N_1/N_2)}{\ln(t_2/t_1)}

throws in metres; the number of faults per area or per unit area; natural logarithms (the ratio of logarithms does not depend on the base)

Holds when: One population of faults: one setting and one stage of deformation. The power law holds over the whole range of extrapolation. Points below the resolution limit are not used for the fit, because the number of faults there is too low.

Faults →
Shale gouge ratio: the share of clay material in the fault zone
SGR=∑iVsh,i Δzit⋅100 %SGR = \frac{\sum_i V_{sh,i}\,\Delta z_i}{t}\cdot 100\,\%

thicknesses and throw in the same units; VshV_{sh} as a fraction

Holds when: The material of the interval that has slid past is fully mixed in the fault zone. VshV_{sh} is calculated in the same way as in the calibration data. SGR thresholds found with a different way of calculating shale volume must not be used. We use the throw at the given point. Along the fault, SGR changes together with it.

Faults →
Pressure difference that a fault holds, and the height of column above the point
Δpf=10 SGR/27−C bar,H=Δpf(ρw−ρh) g\Delta p_f = 10^{\,SGR/27 - C}\ \text{bar},\qquad H = \frac{\Delta p_f}{(\rho_w - \rho_h)\,g}

Δpf\Delta p_f in bar (to calculate HH in SI, multiply by 10510^5 Pa/bar); HH in m; densities in kg/m³

Holds when: The envelope is an empirical upper bound from fields of one region. A real fault may hold less. The growth of the pressure difference held with SGR is confirmed by data up to an SGR of about 40%. Above this value, buoyancy measurements show no further growth, and we do not extrapolate the envelope there. Across the fault there is a water-bearing layer of the same hydrodynamic system, and the pressure difference equals the buoyancy of the column of hydrocarbons (buoyancy). The rocks are water-wet. Calibrations for oil and for gas differ (because of the larger interfacial tension, gas needs a larger pressure difference to break through the same fault rock). Outside the region of the calibration, we check the envelope against local data.

Faults →
Extension and compression of a layer in bending
e=yR,emax=h2R,1R≈ d2z dx2≈z1−2z0+z2a2e = \frac{y}{R},\qquad e_{max} = \frac{h}{2R},\qquad \frac{1}{R} \approx \frac{\dd^2 z}{\dd x^2} \approx \frac{z_1 - 2z_0 + z_2}{a^2}

lengths in the same units; ee is dimensionless (multiply by 100 to get per cent)

Holds when: The neutral surface lies in the middle of the layer, strains are small and the layer does not slip within itself. The curvature from the second difference (the formula with z1−2z0+z2z_1 - 2z_0 + z_2) is valid when the slope at the middle point is small (at the crest  dz/ dx≈0\dd z/\dd x \approx 0). Otherwise use 1/R=z′′/(1+z′2)3/21/R = z''/(1 + z'^2)^{3/2}. We take the profile at right angles to the hinge. In any other direction the curvature comes out smaller.

Folds, salt diapirs and their relationship to faults →
Pressure difference in the salt beneath a minibasin and beneath a dome
Δp=σv(z1)−σv(z2)−ρs g (z1−z2)=g∫z2z1(ρb(z)−ρs)  dz\Delta p = \sigma_v(z_1) - \sigma_v(z_2) - \rho_s\,g\,(z_1 - z_2) = g\int_{z_2}^{z_1}\bigl(\rho_b(z) - \rho_s\bigr)\,\dd z

SI; in metric oilfield units, replace gg by 0.00981 (densities in g/cm³, pressure in MPa)

Holds when: Salt behaves as a fluid: the pressure in it is the same in all directions and increases with depth according to its density. The minibasin and the dome are wide enough for the weight of the column above a point to set the pressure. We ignore the strength of the roof and the friction against the walls, although they slow the flow. We compare the same depth z1z_1. Below it both columns consist of the same salt, so the difference does not change.

Folds, salt diapirs and their relationship to faults →
Dominant wavelength of buckling
Ld=2πh ηl6 ηm3L_d = 2\pi h\,\sqrt[3]{\frac{\eta_l}{6\,\eta_m}}

LdL_d and hh in the same units; only the ratio of the viscosities matters

Holds when: A single viscous layer in a uniform viscous medium; small amplitudes at the start of fold growth. The viscosity contrast is large (tens or more); with a small contrast the layer thickens instead of bending. Multilayer sequences, slip along layers and brittle rocks give other wavelengths. The formula shows a tendency, not an exact value.

Folds, salt diapirs and their relationship to faults →
Griffith envelope and the condition for a joint to form
τ2=4T (σn′+T);joint:  σ3′=−T,σ1−σ3≤4T\tau^2 = 4T\,(\sigma_n' + T);\qquad \text{joint:}\ \ \sigma_3' = -T,\quad \sigma_1 - \sigma_3 \le 4T

MPa; effective stresses use the full pore pressure, as in failure criteria

Holds when: The rock is intact, with no weak planes. On existing fractures and faults, slip starts earlier, by the law of friction (slip on a fracture). The parabola describes the tensile and weakly compressive part of the envelope. At high compressive stress, the real strength is closer to the Coulomb line, and the boundary between hybrid and shear fractures moves a little. The differential stress σ1−σ3\sigma_1 - \sigma_3 does not depend on pore pressure.

Fractures and naturally fractured reservoirs →
True linear intensity of a set from observations in a well (Terzaghi correction)
P10=Nl cos⁡γ,s=1P10,cos⁡γ=∣n⋅w∣P_{10} = \frac{N}{l\,\cos\gamma},\qquad s = \frac{1}{P_{10}},\qquad \cos\gamma = |\mathbf n\cdot\mathbf w|

metres and inverse metres; angles in degrees

Holds when: The fractures of the set are parallel and much longer than the diameter of the well. If the orientations scatter, we calculate the angle γ\gamma for each fracture separately and add up the weights 1/cos⁡γ1/\cos\gamma. When the well is almost parallel to the fractures (γ\gamma close to 90°), the correction is unreliable: the fractures are crossed rarely, and one chance fracture gives a huge contribution. So the weight is limited, usually when the angle between the well and the plane of the fractures is less than 10–15°. The correction cannot restore fractures that the well did not cross: a set parallel to the well may not be noticed at all.

Fractures and naturally fractured reservoirs →
The cubic law: flow rate through a fracture and permeability of a set
q=w3h12 μ⋅ΔpLf,k=∑iwi312 L,ϕf=∑iwiLq = \frac{w^3 h}{12\,\mu}\cdot\frac{\Delta p}{L_f},\qquad k = \frac{\sum_i w_i^3}{12\,L},\qquad \phi_f = \frac{\sum_i w_i}{L}

SI; to convert permeability to millidarcies, divide by 0.9869⋅10−150.9869\cdot10^{-15} m²

Holds when: Laminar flow (smooth and orderly) between smooth parallel walls. For rough fractures with contact points, we use the hydraulic aperture, which is smaller than the average mechanical one. The fractures of the set are parallel, and the flow is along them. Across the set the fractures hardly help, and the matrix decides the permeability in that direction. The matrix adds its own permeability to that along the fractures: kalong≈km+kk_{\text{along}} \approx k_m + k.

Fractures and naturally fractured reservoirs →
Closure of a fracture under load (hyperbolic law)
ΔV=σn′VmKniVm+σn′\Delta V = \frac{\sigma_n' V_m}{K_{ni} V_m + \sigma_n'}

stresses in MPa, closure and VmV_m in µm, KniK_{ni} in MPa/µm

Holds when: An empirical relation. VmV_m and KniK_{ni} are fitted to tests on samples with a fracture, and they depend on the roughness and strength of the walls. Closure is reversible only in part: under repeated load the fracture is stiffer, while slip along it can, on the contrary, open it through the bumps of the walls. We take the change of the hydraulic aperture to equal the change of the mechanical one. For rough fractures, the hydraulic aperture can fall faster as the contacts grow.

Fractures and naturally fractured reservoirs →
Stretching from the rotation of rigid blocks (domino model)
β=LL0=sin⁡δ0sin⁡δ,δ0=δ+δb\beta = \frac{L}{L_0} = \frac{\sin\delta_0}{\sin\delta},\qquad \delta_0 = \delta + \delta_b

angles in degrees; LL and L0L_0 in the same units of length

Holds when: The faults are planar (straight, not curved) and parallel, and the blocks rotate as rigid bodies about axes parallel to the strike. Before extension the beds were horizontal. We measure the tilt of the beds and the dip of the faults in a section across the strike of the faults, otherwise the angles are apparent (apparent dip). The formula describes the extension of only the upper brittle crust cut by these faults, not of the whole lithosphere.

Structural styles and tectonic history →
Directions of structures in a strike-slip zone
αR=φ2,αR′=90°−φ2,αP=−φ2,αT=45°,αF=−45°\alpha_R = \frac{\varphi}{2},\quad \alpha_{R'} = 90° - \frac{\varphi}{2},\quad \alpha_P = -\frac{\varphi}{2},\quad \alpha_T = 45°,\quad \alpha_F = -45°

angles in degrees

Holds when: Uniform simple shear in a wide zone. The greatest compression is horizontal and at 45° to the zone, and shears form by the Coulomb criterion (a rule for when rock fails in shear) at an angle of 45°−φ/245° - \varphi/2 to it. The angles refer to the beginning of the deformation. As the displacement grows, the structures rotate towards the strike of the strike-slip fault, and the shears merge into a through-going fault. If the blocks also move towards each other (convergence), folds and thrusts lie closer to the strike of the strike-slip fault. If they move apart (divergence), normal faults lie closer to it.

Structural styles and tectonic history →
Net throw of an inverted normal fault and the null point
tnet(y)=text(y)−trev,text(y)=hsyn−y,y0=hsyn−trevt_{net}(y) = t_{ext}(y) - t_{rev},\qquad t_{ext}(y) = h_{syn} - y,\qquad y_0 = h_{syn} - t_{rev}

metres

Holds when: There are no syn-rift deposits in the footwall, and the depression in the hanging wall was filled with sediment right to the top during extension. Then the throw of a horizon equals the thickness of the syn-rift layers above it. In compression the fault slides as a whole, and all horizons receive the same reverse throw. If part of the compression went into new cross-cutting reverse faults or folds, the reverse throw on the old fault is less than the total shortening.

Structural styles and tectonic history →
Expansion index and the increase in throw of a growth fault
EI=hhwhfw,Δt=hhw−hfw,tk=∑i>kΔtiEI = \frac{h_{hw}}{h_{fw}},\qquad \Delta t = h_{hw} - h_{fw},\qquad t_k = \sum_{i > k}\Delta t_i

thicknesses and throws in metres

Holds when: Each time, sediment filled the relief created by the fault right to the top, and the footwall was not eroded. The thicknesses in both walls are reduced to one degree of compaction. In the hanging wall the unit is buried deeper and compacted more, so we restore its thickness to the original where possible (restoring thicknesses). The youngest unit is not displaced by the fault. Otherwise its own throw is added to the sum.

Structural styles and tectonic history →
What fraction of the curvature of a fold the second difference shows
z1−2z0+z2a2=∑kAk(2πλk)2sk,sk=(sin⁡(πa/λk)πa/λk)2\frac{z_1 - 2z_0 + z_2}{a^2} = \sum_k A_k\left(\frac{2\pi}{\lambda_k}\right)^{2} s_k,\qquad s_k = \left(\frac{\sin(\pi a/\lambda_k)}{\pi a/\lambda_k}\right)^{2}

lengths in metres, curvature in m⁻¹

Holds when: We write the profile of the top as a sum of cosines. The formula is for a point where all the waves have a crest. At other points each wave enters the sum with its own factor, but the fraction sks_k is the same. The slopes are small, as in the formula for the crest in the lesson on folds. The depths have no errors. A random error in depth is divided by a2a^2, so at a very small spacing the noise distorts the curvature more than the loss of short waves does.

Practical structural geology →
Pressure in an isolated dipping bed and the excess at its top
pr(z)=psh(zc)+Gw (z−zc),Δptop=pr(ztop)−psh(ztop)=(Gsh−Gw) (zc−ztop)p_r(z) = p_{sh}(z_c) + G_w\,(z - z_c),\qquad \Delta p_{top} = p_r(z_{top}) - p_{sh}(z_{top}) = (G_{sh} - G_w)\,(z_c - z_{top})

pressures in MPa, gradients in MPa/m, depths in metres

Holds when: The body is connected and filled with water. Hydrocarbons at the top add their buoyancy (the extra pressure that comes from being lighter than water) to the pressure. The flows through the shales are steady, and the shales let water through equally well over the whole surface of the body. Then the centroid is the mean depth of the surface of the body, weighted by area (the derivation is in the advanced section). The pressure gradient of the shales is constant over the depth interval of the body. If the shales deeper down are more compacted and let less water through, the centroid moves up, and the excess at the crest is less than the calculated value.

Practical structural geology →
Normal stress on a vertical fracture
σn=σHsin⁡2θ+σhcos⁡2θ,σn′=σn−pp\sigma_n = \sigma_H\sin^2\theta + \sigma_h\cos^2\theta,\qquad \sigma_n' = \sigma_n - p_p

stresses and pressure in MPa, angle in degrees

Holds when: The fracture is vertical and one of the principal stresses (the three stresses that act at right angles to one another) is vertical. For an inclined fracture, the contribution of σv\sigma_v is added (advanced section). We measure the angle θ\theta from the strike of the fracture. In stresses on a plane the angle is measured from the normal to the fracture (the line at right angles to it), and from that line σH\sigma_H lies at 90°−θ90° - \theta. The effective stress on the fracture is the total stress minus the pore pressure, as in the lesson on fractures. The stresses are present-day, at a given point; near faults and salt they may differ from the regional ones.

Practical structural geology →
Shortening from the bed length
S=l0−W=∑ili (1−cos⁡δi),ε=Sl0S = l_0 - W = \sum_i l_i\,(1 - \cos\delta_i),\qquad \varepsilon = \frac{S}{l_0}

lengths in metres, angles in degrees

Holds when: Flexural slip: the length and thickness of the bed do not change. Before the deformation the bed lay horizontally; the section is drawn in depth, at true scale, along the direction of transport. We account for displacements on faults within the section; they are added to the bed length.

Practical structural geology →
Excess area and the depth of the detachment
A=S hd,hd=AS,A(z)=S (zd−z)A = S\,h_d,\qquad h_d = \frac{A}{S},\qquad A(z) = S\,(z_d - z)

lengths in metres, areas in m²

Holds when: The area of the unit between the horizon and the detachment does not change: plane strain, no pressure solution or compaction. The detachment is horizontal and there is no deformation beneath it. All horizons are shortened equally. Then the points of the area–depth graph lie on a straight line.

Practical structural geology →

Symbols

SymbolMeaning and context
Δz\Delta z
Δl\Delta l
Δp\Delta p
ΔV\Delta V
Δt\Delta t
Δp, Lf\Delta p,\ L_f
n, w\mathbf n,\ \mathbf w
aa
AA
AkA_k
Ak(2π/λk)2A_k(2\pi/\lambda_k)^2
cc
CC
  • constant that depends on the greatest burial depth: 0.5 for less than 3 km; 0.25 for 3 to 3.5 km; 0 for more than 3.5 km, —Faults
DD
  • slip on the fault plane (net slip), mFaults
DdD_d
  • slip along the dip of the fault plane, mFaults
Dd, DsD_d,\ D_s
  • components of the slip along the dip and along the strike of the fault plane, mFaults
ee
EIEI
Δpf\Delta p_f
  • greatest pressure difference across the fault that it holds (1 bar = 0.1 MPa), barFaults
FF
gg
GshG_{sh}
  • pore pressure gradient of the shales over the depth interval of the body; with undercompaction at a constant effective stress, it equals the overburden pressure gradient, MPa/mPractical structural geology
GwG_w
hh
HH
  • how far below the point the free water level of the accumulation can fall before the fault at this point starts to let hydrocarbons through, mFaults
hdh_d
hhw, hfwh_{hw},\ h_{fw}
hsynh_{syn}
hvh_\text{v}
Δzi\Delta z_i
  • thickness of layer ii within the interval of height tt that has slid past the point: in the section of the footwall, from the point upwards by tt, mFaults
kk
KniK_{ni}
ll
LL
l0l_0
l0, ll_0,\ l
LdL_d
li, δil_i,\ \delta_i
L, L0L,\ L_0
NN
N1, N2N_1,\ N_2
  • numbers of faults with throw not less than t1t_1 and t2t_2; both throws are at or above the resolution limit, —Faults
N(t)N(t)
  • number of faults whose greatest throw is not less than tt in the area, —Faults
pp
PP
P10P_{10}
ppp_p
pr(z)p_r(z)
psh(z)p_{sh}(z)
qq
rr
RR
R, R′R,\ R'
ss
  • vertical separation of the bed: the vertical distance between the parts of one bed in the two walls. In a vertical well it is the thickness of the missing interval (for a normal fault) or the repeated interval (for a reverse fault), mFaults
  • average spacing of the set, mFractures and naturally fractured reservoirs
SS
sks_k
sns_n
  • stratigraphic separation: the same thickness, measured at right angles to the layers, mFaults
SGRSGR
  • share of clay material in the fault zone at a given point of the fault plane, %Faults
  • SGR at this point, %Faults
tt
  • throw (vertical offset), mFaults
  • greatest throw of a fault, mFaults
  • throw of the fault at this point, mFaults
TT
text(y)t_{ext}(y)
tkt_k
tnet(y)t_{net}(y)
trevt_{rev}
t, wt,\ w
Δptop\Delta p_{top}
VV
VmV_m
Vsh,iV_{sh,i}
ww
WW
w, wiw,\ w_i
x, yx,\ y
yy
y0y_0
zz
z1z_1
z1, z0, z2z_1,\ z_0,\ z_2
z2z_2
zcz_c
  • depth of the centroid, where the pressures of the body and the shales are equal; for a body whose area is spread evenly over its depth range, the middle of its height, mPractical structural geology
zdz_d
ztopz_{top}
α\alpha
  • angle between the strike of the structure and the strike of the strike-slip fault; for a right-lateral fault it is measured clockwise, for a left-lateral one anticlockwise, °Structural styles and tectonic history
β\beta
γ\gamma
δ\delta
δ, a\delta,\ a
δ′\delta'
δ0\delta_0
δb\delta_b
  • dip of the beds in a section across the strike of the fault, °Faults
  • tilt of the beds in the blocks — the angle through which the blocks have rotated (the notation of the lesson on faults), °Structural styles and tectonic history
δV\delta_V
ε\varepsilon
ηl, ηm\eta_l,\ \eta_m
θ\theta
θ, φ\theta,\ \varphi
λk\lambda_k
μ\mu
ρ\rho
ρb(z)\rho_b(z)
ρs\rho_s
ρw, ρh\rho_w,\ \rho_h
  • densities of formation water and of the hydrocarbon in the reservoir, kg/m³Faults
σ1, σ3\sigma_1,\ \sigma_3
σ1′, σ3′\sigma_1',\ \sigma_3'
σH, σh\sigma_H,\ \sigma_h
σn\sigma_n
σn′\sigma_n'
σv(z)\sigma_v(z)
τ\tau
φ\varphi
ϕf\phi_f
ψ\psi
ω\omega