100% Automated Vitest Coverage • 18 Benchmark Cases
Mathematical Benchmark Verification
Transparent comparison of computed values against published literature, SPE reference datasets, and textbook worked examples.
| Test Case | Inputs | Expected | Computed | Standard Reference | Status |
|---|---|---|---|---|---|
| API Gravity to Specific Gravity | 35.0° API at 60°F | SG = 0.84985, Density = 53.0 lb/ft³ | SG = 0.84985, Density = 53.0 lb/ft³ | ASTM D1250 Table 5A/6A | Passed (0.00% err) |
| Sutton Pseudo-Critical Properties | Gas Gravity γ_g = 0.65 | Tpc = 365.11 °R, Ppc = 670.13 psia | Tpc = 365.11 °R, Ppc = 670.13 psia | Sutton (1985) SPE 14265 | Passed (0.00% err) |
| Dranchuk-Abou-Kassem (DAK) Z-Factor | Ppr = 2.0, Tpr = 1.30 | Z ≈ 0.684 ± 0.005 | Z = 0.6826 (converged in 4 iters) | Standing-Katz Chart / DAK (1975) | Passed (<0.2% err) |
| Hall-Yarborough Z-Factor | Ppr = 2.0, Tpr = 1.30 | Z ≈ 0.685 ± 0.005 | Z = 0.6849 (converged in 4 iters) | Hall & Yarborough (1973) | Passed (<0.1% err) |
| Colebrook-White Friction Factor | Re = 100,000, ε/D = 0.001 | f ≈ 0.0222 | f = 0.02222 | Moody Chart / Haaland (1983) | Passed (<0.1% err) |
| Fann Viscometer Bingham Rheology | θ600 = 50, θ300 = 30 | PV = 20 cP, YP = 10 lb/100ft² | PV = 20 cP, YP = 10 lb/100ft² | API RP 13D / Bourgoyne | Passed (0.00% err) |
| Directional Survey Minimum Curvature | MD 1000→1100 ft, Inc 10→15°, Az 45→50° | ΔTVD = 97.6 ft, DLS > 0 | ΔTVD = 97.64 ft, DLS = 5.1°/100ft | Bourgoyne Applied Drilling Eng. | Passed (0.00% err) |
| Deterministic STOIIP Volumetrics | 640 acres, 30 ft, φ=0.20, Sw=0.25, Boi=1.25 | 17,874,432 STB | 17,874,432 STB | Craft & Hawkins Ch. 1 | Passed (0.00% err) |
| Petroleum Monte Carlo Exceedance | 5,000 iterations Triangular distribution | P90 < P50 < P10 | P90 = 18.2 < P50 = 24.8 < P10 = 35.6 | Rose (2001) Risk Analysis | Passed (100% compliant) |