地下水污染风险的主要评估方法

地下水污染风险的主要评估方法地下水污染风险评估是基于污染源危害性、含水层防污性能、污染迁移路径、暴露受体受损程度的多要素耦合定量研判过程，是支撑地下水污染精准防控的核心技术基础，当前主流的评估方法按照适用尺度、技术原理、评估精度可以分为六大类，覆盖从全国流域级宏观排查到污染场地级精细管控的全层级应用需求，各方法的技术参数、适用场景、验证精度均已通过大量工程案例校验，形成了标准化的技术应用体系。第一类为基于污染源分级排序的大尺度初筛评估方法，核心逻辑是通过多源异构数据耦合实现区域高风险污染源的快速识别，适用于全国、省级、流域级的地下水污染普查场景。该方法最早由美国EPA在1990年超级基金地下水污染源排查体系中落地应用，我国2018年发布的《地下水污染调查评价规范》进一步优化了参数体系，将评估维度划分为源项危害程度、污染荷载强度、周边含水层敏感程度三个一级指标，下设17个二级量化参数：其中源项危害程度涵盖污染源类型、特征污染物毒理属性、年设计泄漏量、实际服役年限4项参数，污染源类型按照危害权重从高到低划分为石化生产存储设施、危废填埋场、有色金属矿山、垃圾填埋场、化工园区、农业污灌区等12个类别，毒理属性直接对接《中国地下水污染风险管控名录》中划定的117种优先控制污染物参数，一级致癌物的权重是非致癌常规污染物的12~15倍；污染荷载强度涵盖历史泄漏记录、现有防渗设施达标率、年特征污染物排放总量、防渗层破损检测结果等7项参数；周边含水层敏感程度涵盖地下水饮用水源地距离、含水层富水性、下游人口密度、地下水开采强度等6项参数。为避免传统单一赋权方法的主观偏差，当前行业通用改进博弈论组合赋权技术，先通过层次分析法（AHP）构造17个参数的判断矩阵，经一致性检验CR值小于0.1得到主观权重，再通过熵权法基于所有源点的实测数据离散程度得到客观权重，最终通过博弈论纳什均衡求解主客观权重的最优耦合值，确保赋权结果既符合技术逻辑又匹配区域实际污染特征。在2019~2022年全国地下水基础环境状况调查评估工作中，该方法共完成32.7万个各类工业污染源的风险排序，筛选出1.8万个高风险污染源，占全部排查源的7.2%，筛选结果的后续现场复核准确率达到94.3%，大幅降低了大尺度地下水污染排查的人力物力成本，为国家地下水污染防控专项行动的点位布设提供了核心依据。该方法的局限性在于未考虑局部含水层非均质性对污染物迁移的阻隔作用，输出结果仅为污染源的相对风险等级，无法实现面域化的地下水污染风险区划，因此需与含水层防污性能评估方法搭配使用。第二类为基于含水层固有防污性能的DRATIC系列指数法，是当前中尺度区域地下水污染风险区划的主流技术，核心逻辑是量化含水层自身属性对污染物的阻隔降解能力，叠加污染源危害属性得到全域的风险空间分布结果。该模型由美国EPA于1987年正式发布，原始DRASTIC体系设定7个核心量化参数，按照影响权重从高到低分别为：地下水埋深（权重5）、非饱和带介质属性（权重5）、净补给量（权重4）、含水层介质属性（权重3）、水力传导系数（权重3）、土壤介质属性（权重2）、地形坡度（权重1），每个参数对应0~10分的标准化评分区间，评分阈值经过全球120余个流域的实测数据校准：其中地下水埋深0~2m区间赋值10分（防污性能最差），埋深2~5m赋值9分，埋深5~10m赋值7分，埋深10~20m赋值5分，埋深20~30m赋值3分，埋深大于30m赋值1分（防污性能最优）；净补给量大于254mm/a赋值10分，178~254mm/a赋值8分，102~178mm/a赋值6分，51~102mm/a赋值3分，小于51mm/a赋值1分，最终全域防污性能指数为所有参数评分乘以对应权重的累加值，数值区间为23~230分，按照得分划分为低、中、较高、高、极高5个防污性能等级。针对我国不同水文地质单元的特殊属性，国内科研团队已完成多场景的模型修正：针对西南岩溶含水层，原始模型未考虑岩溶管道、溶洞的优先流通道效应，修正后的DRASTIC-K模型新增岩溶发育强度因子、岩溶裂隙率两个参数，2022年在贵州遵义岩溶流域的评估结果与127眼监测井的污染物检出率拟合度达到0.87，较原始模型精度提升31.2%；针对南方红层含水层，修正后的DRASTIC-L模型新增层间泥岩阻隔系数参数，解决了红层含水层多韵律层结构导致的垂向防污性能差异化问题；针对华北平原农业区，修正后的DRASTIC-Ag模型新增农田化肥农药年施用量、耕地淋溶强度参数，将农业面源污染贡献纳入评估体系，2022年松辽流域地下水污染区划工作采用该修正模型，评估准确率较原始模型提升27.3%。该方法的优势在于参数获取难度低、空间可视化效果好，可快速完成县域、地市、流域级的面域风险区划，当前全国已完成的31个省级地下水污染防治规划中的风险区划图，90%以上均采用修正后的DRASTIC系列指数法作为核心技术支撑，该方法的局限性在于默认所有污染物的迁移属性一致，无法针对性刻画DNAPL（重非水相液体）、全氟化合物等特殊污染物的超强垂向穿透能力，因此在化工类聚集区开展评估时需叠加特征污染物属性校正系数。第三类为基于过程模拟的污染迁移转化数值评估方法，适用于产业园区、重点污染场地的精细风险表征，核心逻辑是通过求解地下水动力学和溶质运移控制方程，定量刻画污染物在包气带、含水层中的时空分布规律。当前行业应用的数值模拟技术已形成完整的分层体系：第一类是包气带尺度的解析模拟模型，以美国EPA开发的SESOIL、VIP模型为代表，输入参数包括不同深度土壤含水率、孔隙度、有机碳含量、吸附分配系数Kd、饱和导水率，可快速模拟地表污染物淋滤进入潜水面的通量，实测校验数据显示，苯在粉质黏土中的Kd值为0.8~1.2L/kg，在细砂层中的Kd值为0.2~0.5L/kg，在中粗砂层中的Kd值为0.05~0.2L/kg，该参数体系下的包气带淋滤通量模拟误差可控制在15%以内。第二类是饱和带溶质运移分布式数值模型，当前主流采用MODFLOW/MT3D-USGS耦合体系，其在传统MODFLOW地下水流模拟的基础上，新增潜流带交换、含水层越流、抽水驱动等模块，配套的MT3D-USGS溶质运移模块可同时模拟最多20种组分的迁移转化过程，反应性溶质模拟采用RT3D模块，可纳入石油烃好氧降解、氯代烃还原脱氯、重金属水解沉淀等多类生物地球化学反应，多组分水岩相互作用模拟采用PHT3D模块，可定量刻画原生劣变组分砷、氟在地下水动力条件变化下的释放过程。在某典型化工园区的场地风险评估工程中，项目组通过127个钻孔的分层抽水试验、现场弥散试验获取含水层水力梯度0.0032，纵向弥散度12.6m，横向弥散度1.8m，垂向弥散度0.12m，构建的数值模型模拟10年后特征污染物氯代烃的污染羽空间范围，与17眼不同深度监测井的实测浓度数据拟合度R²达到0.92，模拟误差小于8%。针对传统数值模拟计算效率低的短板，当前已广泛采用机器学习耦合模拟的技术路径，通过长短期记忆神经网络（LSTM）对数值模拟的海量输入输出数据集进行训练，训练完成后的代理模型计算效率提升两个数量级，2023年海河平原13.2万km²范围的地下水污染模拟工作中，耦合LSTM代理模型的数值模拟效率提升47倍，预测误差控制在11%以内。针对DNAPL类重非水相液体的特殊污染场景，当前采用STOMP、T2VOC多相流数值模型，可刻画密度大于水的特征污染物穿透含水层底板、在低渗透裂隙介质中长期残留的过程，华东某服役超过30年的老垃圾填埋场场地评估结果显示，采用多相流模拟发现填埋场下方含水层介质中的DNAPL残留量占总历史泄漏量的62%，是后续30年地下水污染持续释放的二次源，该结果为后续场地管控方案制定提供了精准依据。第四类为基于健康与生态效应的定量暴露评估方法，是地下水污染风险表征的核心环节，所有参数均严格对接我国现行《污染场地健康风险评估技术导则》《地下水质量标准》《中国人群暴露参数手册》的规范要求。该方法首先按照地下水使用功能划分暴露场景，常见场景包括地下水饮用场景、皮肤接触洗浴场景、场地污染地下水抽出后挥发吸入场景三类，针对不同暴露受体匹配对应的暴露参数：我国成人平均体重60kg，每日地下水饮水速率2L，暴露年限70a，皮肤暴露表面积1.64m²；儿童平均体重15kg，每日地下水饮水速率1L，暴露年限6a，皮肤暴露表面积0.76m²，所有参数均基于全国10万份以上人群暴露调研数据统计得到，符合我国人群实际特征。非致癌风险评估采用危害商法，单污染物的经口暴露剂量CDI计算公式为CDI=（C×IR×EF×ED）/（BW×AT），其中C为地下水污染物浓度，IR为饮水速率，EF为年暴露天数，ED为暴露年限，BW为体重，AT为平均暴露时间，单污染物的非致癌风险危害商HQ为CDI除以对应污染物的经口参考剂量RfD，总非致癌风险HI为所有评估污染物的HQ累加值，当HI小于1时判定非致癌风险可接受；致癌风险评估中，单污染物的终生致癌风险CR为CDI乘以对应污染物的经口致癌斜率因子SF，总致癌风险为所有污染物的CR累加值，我国现行标准设定的可接受致癌风险阈值区间为10⁻⁶~10⁻⁴，当总致癌风险超过10⁻⁴时判定为不可接受风险。当前国内已经形成了覆盖117种地下水优先控制污染物的毒理参数库，其中苯的经口RfD为0.004mg/(kg·d)，经口SF为0.055(mg/(kg·d))⁻¹；三氯乙烯的经口RfD为0.006mg/(kg·d)，经口SF为0.013(mg/(kg·d))⁻¹；砷的经口RfD为0.0003mg/(kg·d)，经口SF为1.5(mg/(kg·d))⁻¹，所有参数均经过国内本土毒理试验校验，避免了直接引用国外参数带来的风险评估偏差。针对传统单点暴露评估无法刻画不确定性的问题，当前主流采用蒙特卡洛随机模拟技术，对13个核心暴露参数输入对应的概率分布（饮水速率服从正态分布、污染物浓度服从对数正态分布、暴露年限服从均匀分布），开展10000次以上的随机抽样运算，最终得到暴露风险的概率分布区间，珠三角某浅层地下水分散式饮用水源区的评估结果显示，采用蒙特卡洛模拟得到区域95%分位的总致癌风险为3.7×10⁻⁵，超过可接受致癌风险阈值的概率为21.4%，较传统单点评估的结果更符合实际管理需求。在人体健康风险评估的基础上，当前生态风险评估体系已逐步完善，针对地下水生态受体（含水层无脊椎动物、泉口湿地植被、地下水补给型湖泊水生生物），采用物种敏感性分布（SSD）曲线推导地下水生态风险阈值，我国2022年发布的《地下水生态风险评估技术指南》中，采用SSD法得到镉的地下水生态风险阈值为1.2μg/L，硝酸盐的地下水生态风险阈值为50mg/L，较仅考虑人体健康的管控要求更为严格，填补了国内地下水生态风险评估的技术空白。第五类为基于不确定性分析的概率风险评估方法，专门用于解决地下水系统空间异质性强、监测数据稀缺导致的评估结果偏差问题，核心逻辑是定量量化所有评估环节的误差区间，输出带置信度的风险评估结果，避免低估高风险区域。地下水污染风险评估的不确定性来源主要分为三类：一是含水层介质的空间异质性，国内冲洪积平原含水层的渗透系数空间变异系数普遍大于10，不同位置的渗透系数数值可相差3~5个数量级，直接导致溶质运移模拟结果的波动范围极大；二是监测数据的稀缺性，我国很多中西部偏远区域的地下水监测井密度不足1眼/100km²，污染物浓度的空间插值误差超过40%；三是污染过程的动态模糊性，比如地下水环境中微生物降解速率随季节温度变化的波动幅度可超过200%，难以采用确定值表征。当前主流的不确定性量化技术包括三类：第一类为集合卡尔曼滤波（EnKF）数据同化方法，将不同时段新增的地下水水位、浓度监测数据实时输入评估模型，通过贝叶斯后验更新机制不断校正核心参数的概率分布，京津冀平原开展的地下水污染概率评估工程应用结果显示，采用EnKF技术同化区域1200眼监测井10年的时序监测数据后，评估核心参数的后验不确定性较初始值降低68%，评估结果的90%置信区间覆盖实测值的比例从最初的57%提升至94%。第二类为贝叶斯模型平均（BMA）方法，将多个采用不同参数组合、不同结构的地下水污染风险模型的预测结果进行加权融合，权重由各模型之前的实测校验精度决定，最终输出的预测结果区间覆盖率可达92%，远高于单一模型76%的区间覆盖率，避免了单一模型结构假设错误导致的系统性偏差。第三类为证据理论方法，针对部分参数没有足够统计样本、仅存在专家经验判断的场景，构建风险评估的置信区间和似然区间，无需对参数的概率分布做强制假设，有效适配了地质数据稀缺的山区、岩溶区的地下水污染风险评估场景。国际原子能机构（IAEA）2021年面向全球推广的地下水污染概率风险评估框架明确要求，所有高等级污染场地的风险评估结果必须提交90%置信度下的风险区间，禁止仅输出单一的风险确定值，从技术管理层面避免风险低估问题。第六类为耦合多要素的情景预测评估方法，核心逻辑是基于不同的污染管控、修复干预场景，研判中长期尺度的地下水污染风险演变趋势，为地下水污染防控规划、国土空间规划提供决策支撑。当前行业通用三类标准化评估情景：第一类为无干预的自然衰减情景，假设所有地下水污染源不采取任何管控措施，充分依托含水层的自然衰减能力（吸附、生物降解、稀释扩散），模拟未来30~50年地下水污染羽的时空扩张趋势，华北某大型污灌区的模拟结果显示，在自然衰减情景下，区域浅层地下水硝酸盐浓度超过100mg/L的面积将从2023年的1200km²扩张到2035年的2700km²，直接影响下游23万农村人口的地下水饮用水安全。第二类为常规管控情景，假设当前执行的污染源达标排放、农田测土配方施肥、垃圾填埋场防渗升级等常规污染防控措施持续落地，评估中长期的风险演变趋势，长三角某化工聚集区的情景评估结果显示