🎓 Lesson 1
D1
What Is a Physics-Informed Mining Digital Twin?
A physics-informed mining digital twin is a virtual copy of a real mine that uses real-world physics laws—like how rocks break or how explosives behave—to predict what will happen before it happens.
🎯 Learning Objectives
- ✓ Explain how conservation of momentum governs blast-induced rock motion in a digital twin
- ✓ Analyze the coupling between explosive energy deposition and rock fragmentation using physics-based PDEs
- ✓ Apply the 1D stress-wave reflection model to predict near-field crater formation in a twin simulation
- ✓ Validate twin-predicted muck pile geometry against survey-derived point clouds using RMS error metrics
📖 Why This Matters
In modern open-pit mines, unplanned overbreak, poor fragmentation, or excessive ground vibration cost millions annually in rehandling, delays, and regulatory penalties. A physics-informed digital twin doesn’t just mirror reality—it anticipates it: predicting how a specific blast design will fracture a unique rock mass *before* drilling begins. This shifts blasting from empirical trial-and-error to deterministic, auditable engineering—reducing risk, optimizing energy use, and enabling regulatory compliance through verifiable physics.
📘 Core Principles
A physics-informed digital twin rests on three foundational layers: (1) The geometric-digital layer—a registered 3D model of geology, infrastructure, and equipment; (2) The physics layer—governing partial differential equations (PDEs) for stress wave propagation (e.g., elastodynamic wave equation), fracture mechanics (Griffith criterion), and coupled thermo-chemical detonation modeling; and (3) The data assimilation layer—where real-time drill logs, seismic monitoring, and LiDAR surveys constrain and update the physics model via Kalman filtering or physics-informed neural networks. Crucially, the physics layer ensures causality and generalizability: unlike black-box ML models, it remains valid under unseen rock conditions because it encodes immutable material response laws.
📐 Stress-Wave Peak Particle Velocity (PPV) Prediction
The peak particle velocity at a given distance from a blast source is governed by the scaled-distance law rooted in elastic wave theory. While empirical regressions exist, the physics-informed form derives from spherical wave attenuation in a homogeneous half-space and incorporates rock modulus and density explicitly.
Physics-Based PPV Model
PPV = σ₀ / Z, where Z = ρ × VₚPredicts peak particle velocity at distance R from a blast based on rock impedance and estimated near-field stress amplitude.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| PPV | Peak Particle Velocity | m/s | Maximum ground vibration velocity at receiver location |
| σ₀ | Near-field Stress Amplitude | Pa | Estimated compressive stress wave amplitude generated by explosive energy deposition |
| Z | Rock Impedance | Pa·s/m | Product of rock density and P-wave velocity; governs wave transmission/reflection |
| ρ | Rock Density | kg/m³ | Mass per unit volume of intact rock |
| Vₚ | P-wave Velocity | m/s | Compressional wave speed in rock, directly related to stiffness |
Typical Ranges:
Hard sedimentary rock (coal measure): 3500 – 4800 m/s
Granite: 5000 – 6200 m/s
Weak shale: 1800 – 2800 m/s
💡 Worked Example
Problem: Given: ANFO charge weight = 150 kg, distance from blast face = 80 m, P-wave velocity = 4200 m/s, rock density = 2650 kg/m³, dynamic Young’s modulus = 45 GPa. Estimate PPV using the physics-based model.
1.
Step 1: Compute characteristic impedance Z = ρ × Vp = 2650 kg/m³ × 4200 m/s = 11.13 × 10⁶ Pa·s/m
2.
Step 2: Estimate far-field stress amplitude σ₀ ≈ (E × W^(1/3)) / (4πR²) where E ≈ 1.5×10⁹ J/kg for ANFO → σ₀ ≈ (1.5e9 × 150^(1/3)) / (4π × 80²) ≈ 128 kPa
3.
Step 3: Apply PPV = σ₀ / Z = 128,000 Pa / 11.13e6 Pa·s/m ≈ 0.0115 m/s (11.5 mm/s)
4.
Step 4: Compare to USBM limit (25 mm/s for residential structures): 11.5 mm/s < 25 mm/s → compliant
Answer:
The predicted PPV is 11.5 mm/s, which falls within the safe range of <25 mm/s per USBM criteria.
🏗️ Real-World Application
At BHP’s Mt. Arthur Coal Mine (Australia), a physics-informed digital twin integrated 3D geological modeling, explicit finite-element blast simulations (using LS-DYNA with calibrated rock damage models), and real-time microseismic monitoring. Before a critical highwall blast, the twin predicted localized tensile failure zones exceeding 12 MPa—prompting redesign of delay timing and burden. Post-blast LiDAR survey confirmed 92% accuracy in predicted backbreak geometry, reducing highwall remediation costs by AU$1.7M/year.
🔧 Interactive Calculator
🔧 Open Mine Digital Twin Implementation Calculator📋 Case Connection
📋 Canadian Iron Ore Mine: Blast Performance Twin for Fragmentation Optimization
Over-break damaging ore recovery infrastructure and under-break increasing crushing costs
📋 South African Coal Mine: Digital Twin for Methane Drainage & Ventilation Safety
Intermittent CH₄ spikes triggering false alarms and production halts; inability to distinguish between drainage ineffici...
📋 Norwegian Limestone Mine: Digital Twin for Sustainable Closure Planning
Regulatory requirement for 100-year water quality forecast post-closure; uncertainty in acid rock drainage (ARD) evoluti...