《Journal of Geochemical Exploration》：Long-term evolution of pH and metals (Al, Fe, Mn, U) in pyrite-bearing waste rocks using reactive transport modeling: Dealing with stock and process uncertainties
编辑推荐：
分析了一个前铀矿30年酸性矿山排水（AMD）数据；建立地球化学模型解释每种金属（U、Mn等）的源项；采用一维反应性运移（RT）模拟方法预测未来出水成分；整合不确定性评估出水处理潜在变化。研究人员利用有限的数据集和独特的长期地球化学数据集，确定控制pH和金属浓度
分析了一个前铀矿30年酸性矿山排水（AMD）数据；建立地球化学模型解释每种金属（U、Mn等）的源项；采用一维反应性运移（RT）模拟方法预测未来出水成分；整合不确定性评估出水处理潜在变化。研究人员利用有限的数据集和独特的长期地球化学数据集，确定控制pH和金属浓度时空演化的主要地球化学反应，并开发了简约的RT模型以确认黄铁矿氧化及其他矿物（尤其是U和Mn物种）的过程，同时评估岩石量与组成、热力学及水力参数不确定性的影响。
**论文解读文章**
**研究背景与问题**
矿山废石管理是全球性挑战，每年产生数十亿吨废石，其中含硫化物矿物（如黄铁矿（pyrite）、磁黄铁矿）的废石易生成酸性矿山排水（Acid Mine Drainage, AMD），导致水体酸化、硫酸盐及重金属（如Cu、Zn、As等）浓度升高，威胁地表水、地下水及生态系统。AMD的持续时间常被低估，传统经验模型和短期监测认为金属浓度呈线性下降，但未考虑控制浓度的反应性矿物（如黄铁矿）的长期缓冲作用，导致AMD持续数十年甚至更久。对于老矿山（如法国Bertholène铀矿），停产后废石堆的长期水文地球化学行为存在巨大不确定性，包括矿物学异质性、水文地质复杂性及有限监测数据，阻碍了可持续水管理策略的制定。因此，需要借助反应性运移（Reactive Transport, RT）模型预测长期演化，以优化处理成本与策略。
**研究内容与意义**
研究人员以法国Bertholène铀矿（1977–1994年运营）的尾矿储存设施（Tailings Storage Facility, TSF）为研究对象，该设施含2.8 Mt废石，黄铁矿导致AMD（pH~3，Fe高达0.5 g/L，Al、Mn、U为1–10 mg/L），废水经处理后排放。基于30年（1994–2023年）的强制环境监测数据（pH、Al、Fe、Mn、U浓度及流量），结合矿物学表征（X射线衍射（XRD）、ICP-AES/ICP-MS）和地球化学分析，构建了一维RT模型（HYTEC代码，PRODATA热力学数据库），模拟1992–2052年pH及金属浓度的长期演化，并评估初始黄铁矿储量等不确定性。研究确认了控制AMD的关键反应（黄铁矿氧化、次生矿物沉淀/溶解），并预测未来处理策略可从主动化学处理转向被动处理（如石灰石排水）。该成果发表于《Journal of Geochemical Exploration》。
**主要关键技术方法**
1. **长期环境监测**：1994–2023年，在TSF出口控制盆地（BC）进行每周采样，测pH、Al、Fe、Mn、U浓度及流量；额外在2019年和2023年对井（BE10–BE12、V1、V2）和地表点（DMCO、AVL）进行两次采样，水样经0.2 μm过滤后分析（离子色谱、ICP-SF-MS）。
2. **矿物学表征**：对废石、尾矿和泥浆进行XRD和ICP-AES/MS分析，确定氧化物组成及矿物相（如黄铁矿、针铁矿、石膏、自生铀矿等）。
3. **一维反应性运移模型**：采用HYTEC 4.7.5，将TSF简化为40 m高饱和多孔介质（孔隙度5%，渗透率10?? m/s），顶部施加恒定降雨（509 mm/yr，pH 5.6，O?饱和），底部加入地下水侧向流（P3井成分），模拟1992–2052年，校准矿物初始浓度、比表面积和动力学参数。
4. **不确定性分析**：手动调整黄铁矿初始浓度（0.2–23 mol/L孔隙水），评估对pH和金属浓度的影响，以检验模型鲁棒性。
**研究结果**
- **3.1 水化学时间演化**：通过30年监测数据，发现BC点Al、Fe、Mn、U浓度持续下降，Fe从1994年的~2×10?3 mol/L降至2012年的~2×10?? mol/L，U从1.5×10?? mol/L降至约4.5×10?? mol/L；pH从2007年的3.2升至2021年的3.5。金属总浓度（Al+Fe+Mn）接近法规上限15 mg/L。
- **3.2 化学表征与相平衡**
- **3.2.1 水端元识别**：通过硫酸盐与碱度关系，识别出两个端元——正交片麻岩地下水（低矿化、HCO?主导）和AMD水（高硫酸盐、酸性、零碱度）。TSF水（V1井）为最受影响的端元，BC点由TSF垂直入渗（~27%）与地下水侧向流（~73%）混合。
- **3.2.2 TSF水聚焦**：TSF水中Mg与SO?呈线性相关（指示铁镁矿物风化），Si浓度恒定（与玉髓平衡），Fe高度可变，U在低流期稳定（~1×10?? mol/L），高流期受稀释。
- **3.2.3 反应性矿物识别**：基于饱和指数（SI）计算，P3井地下水与玉髓、微斜长石、针铁矿、三水铝石、软锰矿及偏钙铀云母接近平衡；TSF出水（V1和BC）与玉髓、针铁矿平衡，略欠饱和于石膏和软锰矿，明显欠饱和于偏钙铀云母；与黄铁矿高度欠饱和（SI<-200），与羟铝矾（jurbanite）和黄钾铁矾（jarosite）接近平衡。
- **4. 地球化学模拟与讨论**
- **4.1 酸性排水地球化学**：
- 铝硅酸盐溶解动力学：以微斜长石代表云母溶解，控制K、Mg释放。
- 黄铁矿氧化动力学：采用Williamson和Rimstidt（1994）模型，仅考虑O?氧化（pH>3.4），反应速率由O?和H?催化。
- Fe和Al系统：黄铁矿氧化产生Fe3?，次生矿物包括针铁矿、黄钾铁矾；Al由铝硅酸盐释放，次生矿物为三水铝石和羟铝矾。
- **4.2 矿物学**：废石、尾矿和泥浆的化学分析（表5）显示，废石含大量Fe（8–10%，以黄铁矿为主），尾矿中黄铁矿已酸浸去除，泥浆含石膏、Fe和Al氧化物。初始矿物总量估算见表6。
- **4.3 储存建模**：一维RT模型成功再现了1992–2022年BC点pH缓慢上升（3.8→4.2）及Fe、Al、Mn、U浓度下降趋势（图8）。模型揭示：黄铁矿溶解提供酸度，针铁矿和微斜长石消耗酸度；石膏控制SO?和Ca浓度；偏钙铀云母溶解控制U；软锰矿动力学控制Mn。60年后，黄铁矿几乎耗尽，Fe主要存于黄钾铁矾（82%）。
- **4.4 讨论**
- **4.4.1 长期预测**：模型预测至2052年，pH仍低于法规限值（5.5），但金属总浓度降至15 mg/L以下，U浓度低于1.8 mg/L。未来处理可转向仅调节pH的被动系统（如石灰石排水）。Mn逐渐成为主要金属（2025年占27%），已启用浮石被动处理。
- **4.4.2 不确定性评估**：改变黄铁矿初始浓度（0.6–6 wt%），高含量导致pH过低（<2）、金属浓度过高，无法拟合观测数据。低含量（0.2 mol/L，即反应性部分）才能重现观测。这表明模型对黄铁矿反应性部分敏感，但整体支持处理策略聚焦酸度而非金属。
**结论总结**
Bertholène铀矿30年监测数据为AMD长期演化提供了宝贵约束。研究人员通过RT模型识别出关键过程：黄铁矿氧化驱动酸度，次生矿物（黄钾铁矾、羟铝矾）缓冲pH，石膏控制SO?，针铁矿、微斜长石、软锰矿和偏钙铀云母分别控制Fe、Al、Mn、U的源项。模型成功预测未来25年浓度下降趋势，并表明随Mn成为主要问题，被动处理（浮石、石灰石排水）可替代当前化学处理。本研究为数据有限的老矿山提供了简约建模方法，用于优化长期水管理策略。
