Wellbore Genius
Physics Verification
Runs canonical test cases against the bedding-plane (He-Hutchinson) and TIV (Thomsen) libraries with expected ranges. Use this as a quick smoke check when the underlying libs change.
14 / 14 PASS
Bedding-plane criterion (He-Hutchinson 1989)
6 / 6
Kink-vs-pierce decisions at a weak laminated horizon. Uses ζ = Gc_iface / Gc_matrix and the He-Hutchinson energy bridge.
| Case | Expected | Actual | Status |
|---|---|---|---|
Weak bed, shallow dip (ζ=0.10) Gc_iface/Gc_matrix ≤ 0.25 ⇒ kink along bedding | kink | kink | PASS |
↳ energy ratio ζ Expected ≈ 0.10 | [0.090, 0.110] | 0.100 | PASS |
Tough bed, shallow dip (ζ=1.50) Gc_iface ≥ Gc_matrix ⇒ pierce through | pierce | pierce | PASS |
Mixed regime, sufficient G (ζ=0.50) G_avail ≥ Gc_iface in 0.25<ζ<1 band ⇒ kink | kink | kink | PASS |
Mixed regime, starved G (ζ=0.50) G_avail < Gc_iface in 0.25<ζ<1 band ⇒ pierce | pierce | pierce | PASS |
Steep NF (dip 60°) dip > 15° tol ⇒ defer to RP/GW | not-bedding | not-bedding | PASS |
TIV anisotropic elastic moduli (Thomsen 1986)
8 / 8
Five-constant TIV stiffness from (Vp0, Vs0, ρ, ε, γ, δ) and Mavko et al. (2009) engineering-moduli closed forms.
| Case | Expected | Actual | Status |
|---|---|---|---|
Isotropic collapse: E_h/E_v ε=γ=δ=0 ⇒ ratio ≡ 1 (within 1e-9) | [0.999999999, 1.000000001] | 1.000000000 | PASS |
Shale E_v Typical laminated shale: 25-55 GPa | [25.00, 55.00] GPa | 28.90 GPa | PASS |
Shale E_h Stiffer in bedding plane: 35-80 GPa | [35.00, 80.00] GPa | 42.71 GPa | PASS |
Shale anisotropy E_h/E_v Expected 1.2-2.0 for ε=0.20, γ=0.25 | [1.200, 2.000] | 1.478 | PASS |
E(θ=0°) vs E_v relative deviation Closed-form must agree to <1e-6 | [0.000000000, 1.000000000e-6] | 1.320043280e-16 | PASS |
E(θ=90°) vs E_h relative deviation Closed-form must agree to <1e-6 | [0.000000000, 1.000000000e-6] | 0.000000000 | PASS |
ν_v (vertical Poisson) Physical band 0.10-0.45 | [0.100, 0.450] | 0.236 | PASS |
ν_h (in-plane Poisson) Physical band 0.05-0.45 | [0.050, 0.450] | 0.176 | PASS |