Downhole Ballistics Simulator

Performance modeling

Sign in
Back
Part 4 (Geomech/Sens/MC/Bench)
Continuous-time simulation

Job Timeline — Standard NaNO₃ / Flour

Lumped-parameter ignition → burn → vent simulation. Hold over the charts to read instantaneous values.

Charge & Confinement
5 kg
30 L
80 ms
400 ms
Well baseline
Depth
12,000 ft
BHP
7,676 psi
BHT
88 °F
Peak P
7,676psi
Δ over BHP
+0psi
Rise time
0.0ms
Peak @
0.0ms
Impulse
0psi·ms
Quench
1ms
Pressure (psi)
Temperature (°F)
Gas mass (kg) & burn fraction
Model: Lumped ignition (Arrhenius hot-spot) → pressure-dependent burn rate r₀·(P/P_ref)^n → ideal-gas pressure rise in fixed effective volume → exponential vent with time constant τ. First-order; not a CFD. Use for relative comparison between formulations and confinement geometries.

Multi-cluster stage

analytical stress shadow · Lorentzian decay

Sequentially fired perforation clusters along a lateral. Earlier clusters induce a local minimum-stress increase on later neighbors, derating their effective performance. The shadow kernel can be the simple Lorentzian heuristic or — when fracture geometry is solved below — a Sneddon (2-D plane-strain) elastic kernel with Olson (2004) finite-height 3-D decay.

Fracture Geometry

PKN regime
propagation statusPropagating (depletion-adjusted)
UndepletedGO
K_I/K_IC32.98
Depletion-adj.GO
K_I/K_IC32.98
Rock mechanics
4.00MMpsi
0.22
1000psi·√in
7500psi
100ft
Pump conditions
15.00bpm
50cP
45min
Fracture-plane depletion (Phase 11)
active

Pore-pressure ΔP and temperature ΔT at the fracture face shift σ_min via poro-/thermo-elastic Δσ_h. Width, net pressure and K_I/K_IC are re-solved against σ_min_eff.

-2000psi
-50°F
Δσ_poro
-1149psi
Δσ_thermo
-1282psi
Δσ_total
-2431psi
σ_min_eff
5069psi
K_I/K_IC eff
32.98
Δ(K_I/K_IC)
0.00
Phase 12 BC impact§12.2–§12.5
K_IC eff (§12.5)
1000psi·√in
BHP drive (§12.3)
-1069psi
  • · §12.3 BHP drive = BHP − σ_min_eff = 4000 − 5069 = -1069 psi (below closure).
Solved geometrydepletion-adjusted
Half-length L
2622ft
Height h
100ft
Width w_max
0.085in
Width w_avg
0.067in
Net pressure
148psi
K_I / K_IC
32.98
Depletion-adjusted pressures
p_net,eff (above σ_min,eff)
148psi
p_wb shift vs σ_min₀
-2282psi
p_wb − σ_min₀ = p_net,eff + Δσ_total = 148 + (-2431) = -2282 psi. Depleted rock — surface pressure to maintain this width drops by 2431 psi vs the undepleted reference.
ΔL vs base0.0ft
Δw_avg vs base0.0mil
Δp_net vs base0.0psi
Δ(K_I/K_IC)0.00
  • · PKN selected: L > h ⇒ height-contained planar growth.
  • · Toughness-dominated: net pressure exceeds elastic prediction.
  • · Depletion Δσ_h = -2431 psi (poro -1149 + thermo -1282); σ_min_eff = 5069 psi.
  • · §12.3 BHP drive = BHP − σ_min_eff = 4000 − 5069 = -1069 psi (below closure).

Depletion sensitivity sweep

w_avg and p_net,eff vs ΔP (pore pressure). Vertical line = current run value.

Stage geometry
6
20 ft
Shadow kernel
0.02 %/psi
Mean perf
101.4%
Cluster eff
100%
Worst
c6
Rec. spacing
45ft
Lateral · 6 clustersheel → toe (left → right)
HEELTOE0ft25ft50ft75ft100ft1061c11022c21013c31004c41005c51006c6
≥95% lead85–95%<85%number = fire order · badge = effective %
Stress Shadow Matrix · Δσ (psi)row = receiver · col = source · max 180 psi
c1c2c3c4c5c6c1c2c3c4c5c618062180296218017296218011172962180
Select a cluster

Click any cluster marker on the lateral to inspect its per-source stress contributions and effective performance.

Model: Each cluster fired before cluster i contributes Δσ via a Lorentzian kernel σ₀·a²/(d²+a²). Effective performance η = η_base − k·Σ Δσ. Cluster efficiency = fraction firing within 90% of the lead cluster. For pre-job spacing screening; calibrate k against measured production logs before field decisions.