突发水污染事故全周期生态损害评估框架：多模型耦合预测与三维定量评估——以尾矿库溃坝为例

《Toxics》：Full-Cycle Ecological Damage Assessment Framework for Sudden Water Pollution Accidents: Multi-Model Coupled Prediction and Three-Dimensional Quantitative Evaluation with a Case Study of Tailings Dam Breach

　　

突发尾矿库溃坝事故在耦合的地表水-地下水系统中引发大规模重金属复合污染，需要适用于不同污染类型的系统性全周期生态损害量化工具。本研究构建了一个用于突发水污染事故的集成全周期生态损害评估框架，整合了三个核心模块：多模型污染物迁移预测、多尺度水生生物损害诊断以及三维生态-经济损失核算。该框架采用模块化设计，可潜在地容纳重金属（镉Cd、铬Cr、砷As、铅Pb）和多环芳烃（PAHs）等有机污染物，并预留了标准化的分子、个体和种群水平生物终点以及相应的污染物剂量-响应模板作为参考计算模块。然而，其在本案例之外的适用性尚未验证，需要针对具体案例进行校准。为验证所提出集成系统的可操作性和准确性，研究人员选择了一个以六价铬（Cr(VI)）和砷（As）污染为主的典型尾矿库泄漏事件作为实际验证案例；该案例中所有现场监测、污染物模拟和最终经济损失量化均完全依赖于现场实测的Cr(VI)和As数据，而手稿中保留的Cd和PAH相关生物响应曲线及修复成本公式仅作为框架的通用模板组件，而非案例实测结果。针对该Cr(VI)/As污染案例，对流-扩散模型模拟显示，Cr(VI)污染羽在48小时内水平扩散250米，七天后延伸至560米，且厌氧地下水环境驱动了原生五价砷向毒性更高、迁移性更强的三价砷（As(III)）的转化。校准后的SWAT模型对溶解态Cr(VI)和颗粒态As的纳什-萨特克利夫效率（NSE）系数分别达到0.75和0.68。基于图论的快速预测模型将计算时间缩短至分钟级；经独立现场监测数据验证，其平均相对误差为14.2%，与SWAT模型的一致性相对偏差为10.5%，满足应急预警的精度要求。现场生物监测表明生态损害显著：鱼类组织中的金属硫蛋白（MT）表达显著升高（所报告的6.2倍诱导值来源于框架内标准Cd暴露模板测试，现场Cr(VI)/As共胁迫下也观察到类似的MT上调），且受污染河段底栖群落香农多样性下降了50%以上。该案例的标准化生态损害指数（EDI）计算为480.2，表明水生生态系统严重受损，综合生态-经济损失总计达1725万元人民币。本研究创新性地将高精度物理传输模型与快速应急预测算法相耦合，并建立了从分子生物标志物到群落完整性指标的完整多层次生物指标链；整合生态系统服务损害、修复支出和污染后恢复损失的三维损失核算系统实现了全周期损害评估的闭环。所提出的框架以Cr(VI)和As污染为例进行了验证，其模块化设计可通过调整模型参数潜在地扩展到Cd和PAHs等其他污染物，为水污染事故的应急处置、污染修复和生态补偿提供定量参考，但需要针对不同污染物和水文条件进行进一步验证。

**论文解读文章**

**研究背景与问题**

随着全球经济发展和工业化进程，人类活动对环境的干扰日益加剧，水生环境作为关键环境载体面临严重污染压力。重金属作为无机元素，在环境介质中无法被微生物或化学降解为无害终产物，具有持久性和环境持久性，易在土壤、沉积物和生物体中积累，并沿食物链传递和生物放大，对水生生态系统稳定性和人类健康构成长期风险。尾矿库溃坝等突发环境事件导致重金属以高释放率、大范围扩散方式进入水体，带来显著生态影响。然而，现有研究在污染物扩散预测、生态损害量化、修复评估和决策支持等方面仍存在关键问题：污染物扩散预测模型常受限于数据依赖性和时效性；生态损害评估多集中于单一或孤立的生物指标，缺乏分子、个体、种群、生态系统多层次的系统整合；修复技术多针对单一污染物或单一污染途径，缺乏系统性的快速响应评估与修复策略。

**本研究内容与结论**

针对上述问题，研究人员以典型尾矿库溃坝事故造成的六价铬（Cr(VI)）和砷（As）复合污染为研究对象，构建了“多模型耦合预测+多层次生物监测+三维综合评估”的集成研究框架。该框架包括：利用对流-扩散方程模拟地下水三维迁移；利用SWAT模型模拟复杂水文条件下地表水污染物动态扩散；构建基于图论的快速预测模型以满足应急响应需求；建立涵盖分子-个体-群体-生态系统的四层生物损害指标体系；采用“生态服务损害-修复成本-恢复期损失”三维评估框架计算污染总经济损失。研究得出以下结论：在砂质含水层中，水力弥散系数是控制污染物羽扩散的主导因素；五价砷在厌氧区被还原为高迁移性三价砷，扩散范围扩大了235%；SWAT模型对溶解态Cr(VI)和颗粒态As的纳什系数分别为0.75和0.68；基于图论的快速预测模型将计算时间从小时级降至分钟级，平均相对误差14.2%，与SWAT模型偏差10.5%；现场生物监测显示Cr(VI)/As共污染河段鱼类MT表达显著上调，底栖生物多样性下降超过50%；综合生态-经济损失达1725万元人民币，生态损害指数（EDI）为480.2，表明严重受损。该研究实现了水文学建模、生态毒理学和环境经济学的跨学科整合，填补了突发重金属污染全周期评估的空白，相关成果发表在《Toxics》期刊。

**关键技术方法**

研究人员主要采用了以下关键技术方法：（1）基于对流-扩散方程的地下水污染物迁移模型，利用SUFI-2算法进行参数校准和不确定性分析；（2）半分布式水文模型SWAT，用于地表水污染物迁移模拟，基于SRTM DEM、Landsat 8土地利用数据和HWSD土壤数据构建，设置3年预热期（2017-2019）；（3）基于图论的快速预测模型，将河流系统抽象为有向图，转化为图遍历问题，适用于应急响应前72小时快速预测；（4）多层次生物损害指标体系，涵盖分子/细胞水平（金属硫蛋白MT、氧化应激指标）、个体/生理水平（孵化抑制、回避行为）和种群/群落水平（香农多样性指数、生物完整性指数IBI）；（5）三维生态-经济损失核算框架，整合生态系统服务功能损害评估、环境修复成本核算（包括原位化学还原ISCR、可渗透反应墙PRB、分期疏浚）和生态系统功能恢复期损失评估。所有子方法均为已发表研究中的成熟工具，样本队列来源为长江中下游某工业园区的尾矿库泄漏事故现场监测数据。

**研究结果**

**3.1. 地下水污染物扩散预测结果**
基于对流-扩散方程模型，利用SUFI-2算法进行参数校准和不确定性分析，模拟了Cr(VI)和As在砂质含水层中的时空扩散特征。结果显示：泄漏后48小时，Cr(VI)污染羽水平扩散250米，浓度峰值降至28.7 mg/L；As扩散较慢，水平扩散140米，浓度峰值10.2 mg/L。7天后，Cr(VI)扩散至560米，垂直入渗深度8米，浓度峰值降至12.3 mg/L；初始五价砷因厌氧条件（Eh = -80 mV）被还原为高迁移性三价砷，扩散范围扩大至470米（净增加235%）。模型校准集NSE：Cr(VI)为0.82，As为0.76；验证集NSE分别为0.78和0.71。参数敏感性分析表明水力弥散系数贡献了42%的预测不确定性，95%预测不确定性区间包含89%的观测数据。

**3.2. 地表水污染物扩散预测结果**
基于SWAT模型和基于图论的快速预测模型，分析了不同水文条件下Cr(VI)和As的动态扩散过程。SWAT模型模拟显示：泄漏后24小时，Cr(VI)高浓度带沿河延伸8公里，支流分流使主河浓度峰值降至0.35 mg/L；48小时后峰值降至0.18 mg/L。暴雨事件后4.5小时，地表源Cr(VI)负荷峰值达1.1 mg/L，导致1-5公里河段普遍超标；暴雨后24小时，颗粒态As富集于河底沉积物，浓度达28.5 mg/kg。若降雨增至150 mm（20年一遇），Cr(VI)扩散范围扩大至12公里，底泥再悬浮导致As二次污染，水体As浓度升至0.3 mg/L。SWAT模型验证期NSE为0.78（流量），Cr(VI)和As模拟值与监测值吻合良好。基于图论的快速预测模型将河流系统抽象为有向图，泄漏后2小时Cr(VI)污染羽前锋到达下游1公里断面，浓度12.3 mg/L；8小时迁移至5公里断面，浓度0.42 mg/L，与SWAT模拟结果（0.38 mg/L）相对偏差10.5%，与现场实测浓度（0.39 mg/L）相对误差7.7%。模型验证显示平均相对误差14.2%，均方根误差0.055 mg/L，决定系数R²=0.87，应急响应窗口（72小时）内平均相对误差为13.1%。

**3.3. 水生生物损害指标分析**
以目标河流及周边水体为研究对象，选择铜绿微囊藻、霍甫水丝蚓和鲫鱼为指示生物，建立多层次生物损害指标体系。通过广义加性模型（GAM）和结构方程模型（SEM）耦合分析，揭示了环境因素对生物损害的综合作用机制。分子/细胞水平：框架标准Cd参考模板显示金属硫蛋白（MT）表达呈饱和增长趋势，低浓度区（Cd<0.1 mg/L）高度敏感，浓度>0.3 mg/L时趋于饱和（米氏动力学）；现场Cr(VI)/As共污染河段鱼类MT浓度与水中Cr(VI)浓度呈显著正指数相关。个体/生理水平：Cd 0.09 mg/L导致鲫鱼孵化抑制、畸形率增至28%，幼鱼体长缩短35%；As>20 mg/kg时霍甫水丝蚓回避率升至65%，掘穴深度从12 cm降至5 cm。种群/群落水平：污染组浮游藻类香农多样性指数H′=1.2（基线/对照组H′=2.8），均匀度指数J=0.45（基线/对照组J=0.82）；底栖动物H′=0.9（基线/对照组H′=2.5），耐污类群（如水丝蚓科）比例从15%增至65%，指示生态系统功能退化。基于毒代动力学-毒效动力学（TK-TD）框架计算，鲫鱼肌肉中Cr(VI)累积浓度30天后达0.85 mg/kg，超过食品安全国家标准限值（0.1 mg/kg）7.5倍。

**3.4. 综合评估**
基于“生态服务损害-修复成本-恢复期损失”三维评估框架，计算全周期污染损害。生态系统服务功能损害：以渔业损失为核心，采用营养传递效率η=0.12、底栖生物量损失、热力学损失系数等参数，计算一年内目标河流次级消费者经济损失为185万元；水体自净功能（氨氮降解）损害采用Cr(VI)特异性参数（最大抑制率V_max,Cr、半效应浓度EC_50,Cr等），年损失约677万元，累计贴现损失为862万元。环境修复成本：针对Cr(VI)污染采用原位化学还原（ISCR），计算还原剂（连二亚硫酸钠）投加量4.8吨，总成本52万元；针对As污染采用可渗透反应墙（PRB），成本86万元；沉积物分期疏浚成本185万元，总修复成本323万元。生态系统功能恢复期损失：以沉水植被恢复为核心，模拟恢复路径，计算恢复期损失540万元。综合生态-经济损失总计1725万元。生态损害指数（EDI）计算公式为EDI = w₁×L₁ + w₂×E₁，其中供给服务权重w₁=0.4，调节服务权重w₂=0.6，L₁=185万元（渔业损失），E₁=677万元（自净损失），计算得EDI=480.2，对应年基线生态系统服务价值（650万元/年）的74%，确认严重受损状态。

**讨论与结论**

研究人员将本研究结果与全球典型尾矿库溃坝事件（如巴西Fund?o溃坝、西班牙Aznalcóllar矿泄漏）进行了对比分析，证实本案例的污染物迁移距离、模型性能指标（SWAT NSE 0.68-0.75）以及生态损害程度（底栖多样性下降>50%等）均处于国际公认的合理范围内。研究结论指出：在砂质含水层中，水力弥散系数是控制污染羽扩散的主导因素，五价砷在厌氧区被还原为三价砷使扩散范围扩大235%；SWAT模型对Cr(VI)和As模拟效果良好，暴雨事件因沉积物再悬浮触发显著脉冲污染；基于图论的快速预测模型将计算时间从小时级降至分钟级，平均相对误差14.2%，与SWAT模型偏差10.5%，在应急响应场景中具有优异应用潜力。多尺度生物损害通过框架指标体系得到系统诊断，框架标准Cd参考模板显示MT表达升高6.2倍，现场Cr(VI)/As共污染河段确认MT显著上调、底栖生物多样性下降超过50%。综合生态-经济损失达1725万元，EDI为480.2，确认水生生态系统严重受损。该研究在理论上建立了重金属污染预测的耦合模型系统，阐明了参数敏感性的时空异质性，构建了多层次生物损害指标体系和三维评估框架，实现了水文学、生态毒理学和环境经济学的跨学科整合；在实践上，图论模型为应急响应提供快速预测工具，参数敏感性结果指导修复技术优化，综合评估方法为生态损害补偿、修复资金分配和司法仲裁提供定量依据。尽管所提出的耦合模型系统表现良好，但存在以下局限性：对流-扩散模型未考虑胶体促进重金属迁移，可能低估长距离迁移潜力；SWAT模型对极端暴雨事件模拟精度较低（NSE=0.62），源于气象数据时间分辨率有限；图论模型在复杂支流结构河流系统中性能下降。框架的模块化设计表明通过调整参数可潜在地适用于其他重金属和有机污染物，但该单案例研究仅展示潜在适用性而非普遍有效性，需要针对不同污染物、事故类型和水文条件进行进一步验证。