基于和声搜索算法的地下水污染源与未知含水层参数的同步反演研究 Simultaneous identification of groundwater contaminant source and aquifer parameters by harmony search algorithm.pdfVIP

水 利 学 报 2012年 12月 SHUILI XUEBA0 第 43卷 第 12期 文章编号 ：0559—9350(2012)12-1470—08 基于和声搜索算法的地下水污染源与未知含水层 参数 的同步反演研究 江思珉，蔡 奕，王 敏2，周念清 (1．同济大学 水利工程系 ，上海 200092；2．山东科技大学 地质科学与工程学院，山东 青岛 266510) 摘要 ：地下水污染源反演 问题和含水层参数反演 问题都是典型的地下水逆 问题 。在未知含水层参数(渗透系数、 弥散度等)等先决信息的情况下进行地下水污染源反演计算时，需要根据已有 的监测数据 (水位和浓度等)对地下 水污染源和未知含水层参数进行 同步反演 。在 同步反演优化问题中，决策变量包括污染源位置 、强度 以及待求的 含水层参数 。本文对 同步反演模 型的框架组成 (包括污染物迁移模型和反演优化模 型)进行分析后 ，在对 已有 的各 种和声搜索改进算法进行研究的基础上结合 同步反演模型提出一种改进 的和声搜索算法 ，最后将 同步反演模型和 改进 的和声搜索算法应用于具体的算例研究。研究表明，改进的和声搜索算法具有算法稳定高效 、求解精度高等 特点 ，能够广泛应用于复杂的地下水污染源和含水层参数反演 问题 。 关键词 ：同步反演 ；污染源反演 ；含水层参数反演 ；和声搜索算法 ；地下水逆问题 中图分类号 ：P641．2 文献标识码 ：A 1 研究背景 目前 ，我国地下水污染已十分严重。与地表水体相比，地下水一旦遭受污染其治理和修复往往 需要更长的时间(常常持续30年以上)，其治理和修复费用也极为昂贵。通常来说 ，要对地下水进行 有效 的治理和修复 ，需要先大致弄清污染源的情况 ，包括污染源分布以及污染浓度排放情况等 ，该 过程被称为地下水污染源反演…。一般所说的地下水污染源反演过程 ，首先是假定含水层参数和边界 条件为已知 ，然后在此基础上根据有限且离散的实测数据进行污染源反演 。然而，在实际问题 中 含水层参数 (渗透系数、弥散度等)等先决信息往往是未知的，此时就需要根据已有的实测数据对地 下水污染源和含水层参数进行 同步反演 。 地下水污染源反演 问题和含水层参数反演 问题都属于地下水逆 问题 ，而最优化方法是求解地下 水逆问题 的有力工具 。传统的优化方法有Newton法 、单纯形法 、共轭梯度法 、Hooke—Jeeves方法 等 ， 述方法均属于局部邻域搜索算法 ，因而易陷入局部最优而难 以获得全局最优解 。近年来 ， 能搜索到全局最优解 的启发式搜索算法 (遗传算法 、模拟退火算法 、禁忌搜索算法 、粒子群算法等) 正受到人们越来越多的青睐 ，并得到了成功的应用…。。 本文采用的和声搜索算法 (HarmonySearchalgorithm，Hs)是一种全新的启发式全局搜索算法 。 Ayvaz 采用Hs算法结合地下水流模型MODFLOW对二维地下水管理模型进行 了求解 ，并与GA(Ge— neticAlgorithm，GA)、SA(SimulatedAnnealing，SA)、GASAPF(GeneticAlgorithmbasedSimulatedAn— nealingPenaltyFunction，GASAPF)以及SCE—uA(ShuffledComplexEvolution，SCE—UA)等算法进行 比 较 ，Hs算法体现 出更好 的优化性能 。Ayvaz还结合地下水污染物迁移模型MT3DMS，将 Hs算法推 广应用于地下水污染源反演问题 (已知含水层参数情形下)。骆乾坤等 采用改进的HS算法对二维承 收稿 日期 ：2012-02—10 基金项 目：国家 自然科学青年基金资助项 目；光华同济土木学院基金资助 作者简介 ：江思珉(1980-)，男 ，安徽歙县人，讲师，博士 ，主要从事地下水数值模拟研究 。E-mail：jiangsimin@tongji．edu．ca 一 1470 — 万方数据 压含水层进行 了含水层参数识别。上述研究均表明HS算法对初始参数值敏感度低 ，具有很强的全局 寻优