基于污染羽形态对比的地下水污染源识别研究.pdf

水 利 学 报 2014年6月 SHUILI XUEBAO 第45卷 第6期 文章编号：0559-9350（2014）06-0735-07 基于污染羽形态对比的地下水污染源识别研究 江思珉，张亚力，周念清，赵 姗 （同济大学 水利工程系，上海 200092） 摘要：地下水污染源位置与强度的确定有助于提高地下水污染治理和修复的效果。本文在场地信息不确定的前提 下，基于卡尔曼滤波技术与单纯形法提出一种地下水污染源识别的新方法。该方法利用卡尔曼滤波技术连续作用 采样点来估计复合污染羽与误差协方差矩阵，进而利用误差协方差矩阵进行新采样点的选择；模糊集被用来进行 污染羽的表示，通过复合污染羽与单个污染羽的形态对比以确定污染源位置；在进行污染源位置反演时，嵌入单 纯形法以进行污染源强度的同步反演。算例研究表明，该方法在判断污染源位置时考虑了污染羽的整体形态，从 而降低了场地局部信息不确定对于污染源识别结果的影响，并能够通过合理的采样点的选择，正确地识别出污染 源位置与强度，反演结果具有较高的可靠性。 关键词：污染源识别；模糊对比；单纯形法；卡尔曼滤波 中图分类号：P641.2 文献标识码：A doi：10.13243/j.cnki.slxb.2014.06.013 1 研究背景 与地表水相比，地下水一旦遭受污染其治理和修复往往需要更长的时间（常常持续30年以上）， 其治理和修复费用往往极为昂贵。由于污染物随地下水流向多呈羽状扩散，污染羽状体一般被称为 污染羽。通常来说，要对地下水进行有效的治理和修复，需要大致弄清污染源信息（包括污染源分布 以及污染物排放信息等）以及含水层中污染羽的时空分布情况。这些都需要通过钻孔取样建立地下水 污染监测网，以获取地下水的物理特性、化学特性和生物特性在时间、空间上的变化特征数据，仅 ［1］ 仅通过钻探取样的方法很难直接确定污染源的位置以及污染物释放强度 。 在地下水水质的监测过程中，过量取样和取样不足是比较常见的问题。过度取样会造成不必要 的取样和分析费用，而取样不足则不能获得污染物随时空变化的准确信息，进而难以确立地下水污 ［2］ 染治理和修复的合理方案 。文献［3-7］通过对地下水污染监测网进行优化设计，可以在不降低监测 精度及不影响评价风险的前提下减少不必要的取样或增加必要的取样，从而以尽可能少的费用识别 出地下水污染源。 地下水系统是一个受多种因素影响的复杂系统，存在着各种不确定性。在实际问题中，因实测 ［8-10］ 数据有限，含水层参数（如渗透系数等）常具有不确定性 。场地信息的不确定性还包括由于抽（注） 水井、排水沟渠等调查数据的不完备所引起的不确定性。以上这些因素都会增加地下水污染源识别 问题的复杂性和困难度。本文在场地信息不确定的前提下，基于污染羽形态对比方法（以往的研究多 采用观测值偏差最小化［11-12］），采用卡尔曼滤波技术、模糊集合理论和数学优化方法提出一种地下水 污染源识别的新方法。 收稿日期：2013-06-04 基金项目：国家自然科学基金资助项目；中央高校基本科研业务费专项资金资助 作者简介：江思珉（1980-），男，安徽歙县人，博士，主要从事地下水数值模拟研究。E-mail：jiangsimin@tongji.edu.cn — 735 — 2 污染源反演方法 在进行地下水污染治理和修复过程中，通常要求以尽可能少的水质采样点达到准确识别地下水 污染源的目的。在考虑场地信息不确定性的情形下，一般选择使污染浓度场的不确定度减少最大的 点作为采样点［13-14］。本文基于污染羽形态对比方法，采用卡尔曼滤波技术、模糊集合理论和数学优化 方法提出一类地下水污染源识别的新方法，