Closure Twin Geochemical Boundary Condition Generator
The Closure Twin Geochemical Boundary Condition Generator is an Excel-based computational resource designed to dynamically define and propagate time-variant geochemical boundary conditions—such as porewater chemistry, redox potential, and mineral saturation states—into Mine Digital Twin frameworks for post-closure environmental performance assessment. It bridges site-specific geochemical monitoring data, thermodynamic models, and regulatory compliance timelines to generate auditable, scenario-aware boundary inputs for reactive transport simulations. Its core purpose is to ensure that digital twin representations of mine closure systems reflect realistic, evolving geochemical constraints over centuries-scale horizons.
📖 Overview
📑 Key Components
🎯 Applications
- ✓ Input generation for reactive transport models in mine closure certification
- ✓ Scenario testing of cover system longevity under climate change projections
- ✓ Audit-ready documentation of geochemical assumptions for third-party review and permitting
📐 Key Formulas
Redox Potential Conversion (Eh-pH)
Eh = E⁰ − (0.05916 / n) × log₁₀(Q) − 0.05916 × pH
Calculates aqueous redox potential (Eh, volts) at a given pH using Nernst equation for dominant redox couples (e.g., Fe³⁺/Fe²⁺, SO₄²⁻/HS⁻), where E⁰ is standard potential, n is electron transfer count, and Q is reaction quotient.
Mineral Saturation Index (SI)
SI = log₁₀(IAP / K_sp)
Quantifies departure from equilibrium for a mineral phase; SI = 0 indicates saturation, SI > 0 supersaturation (precipitation likely), SI < 0 undersaturation (dissolution likely), where IAP is ion activity product and K_sp is solubility product.
Infiltration-Driven Boundary Flux
J_i(t) = I(t) × C_i(t) × A_eff
Estimates time-dependent solute mass flux (kg/yr) across a boundary layer, where I(t) is infiltration rate (m/yr), C_i(t) is dissolved concentration of species i (mg/L), and A_eff is effective area (m²) accounting for cover heterogeneity.