📦 Resource excel

PHREEQC Input File Generator for Mine Drainage Systems

The PHREEQC Input File Generator for Mine Drainage Systems is an Excel-based computational tool designed to automate the creation of input files compatible with PHREEQC—a widely used geochemical modeling software—for simulating water-rock interactions in acid mine drainage (AMD) and neutral mine drainage (NMD) environments. It translates site-specific hydrochemical data, mineralogy, and reaction parameters into standardized PHREEQC syntax (SOLUTION, PHASES, REACTION, TRANSPORT blocks), ensuring reproducibility and reducing manual error. The tool bridges field/lab data collection with reactive transport modeling, supporting risk assessment and remediation design.

📖 Overview

Mine drainage systems pose significant environmental challenges due to the release of metals, sulfate, and acidity from sulfide oxidation in waste rock and tailings. Accurate prediction of long-term water quality requires robust geochemical modeling—PHREEQC excels in speciation, saturation index calculation, kinetic mineral dissolution/precipitation, and one-dimensional reactive transport simulations. However, constructing syntactically correct PHREEQC input files manually is time-intensive and error-prone, especially when handling multiple sampling points, time-series data, or scenario-based sensitivity analyses. The Excel-based generator addresses this by providing structured worksheets for entering analytical data (e.g., pH, major ions, trace metals), solid-phase composition (e.g., pyrite, gypsum, ferrihydrite), thermodynamic database selection (e.g., minteq.v4.dat, phreeqc.dat), and user-defined reactions (e.g., O2 diffusion, organic matter oxidation). It applies built-in validation checks (e.g., charge balance tolerance, unit consistency, redox pairing) and auto-generates commented, modular PHREEQC scripts with embedded documentation. Advanced versions support batch generation for spatial grids or temporal scenarios, enabling rapid comparative analysis of treatment options (e.g., limestone dosing vs. anoxic limestone drains) or climate-driven changes in recharge chemistry.

📑 Key Components

1 Hydrochemical Data Input Sheet
2 Mineralogical & Reactive Surface Database
3 PHREEQC Syntax Translator Engine

🎯 Applications

  • Predictive AMD plume evolution under varying hydrological conditions
  • Design and optimization of passive treatment systems (e.g., constructed wetlands, alkalinity-generating substrates)
  • Regulatory compliance reporting and long-term monitoring plan development

📐 Key Formulas

Charge Balance Error (CBE)

CBE (%) = 100 × |Σcations − Σanions| / max(Σcations, Σanions)

Quantifies analytical reliability of water sample; values >5% typically indicate measurement error or unanalyzed species (e.g., colloidal Fe, organic acids)

Saturation Index (SI)

SI = log₁₀(IAP / K_sp)

Measures disequilibrium of a mineral phase; SI ≈ 0 indicates saturation, SI > 0 supersaturation (potential precipitation), SI < 0 undersaturation (potential dissolution)

Effective Diffusion Coefficient (D_eff)

D_eff = τ × D₀ / ϕ

Accounts for tortuosity (τ) and porosity (ϕ) in estimating solute transport through porous media; D₀ is molecular diffusion coefficient in water

🔗 Related Concepts

Acid Mine Drainage (AMD) Reactive Transport Modeling Geochemical Speciation

📚 References

#geochemistry #mine reclamation #water quality modeling #PHREEQC #Excel automation