基于和声搜索算法地下水污染源与未知含水层参数同步反演的研究.pdfVIP

水 利 学 报 2012年12月 SHUILI XUEBAO 第43卷 第12期 文章编号：0559-9350（2012）12-1470-08 基于和声搜索算法的地下水污染源与未知含水层 参数的同步反演研究 江思珉 ，蔡 奕 ，王 敏 ，周念清1 1 2 1 （1. 同济大学 水利工程系，上海 200092；2. 山东科技大学 地质科学与工程学院，山东 青岛 266510） 摘要：地下水污染源反演问题和含水层参数反演问题都是典型的地下水逆问题。在未知含水层参数（渗透系数、 弥散度等）等先决信息的情况下进行地下水污染源反演计算时，需要根据已有的监测数据（水位和浓度等）对地下 水污染源和未知含水层参数进行同步反演。在同步反演优化问题中，决策变量包括污染源位置、强度以及待求的 含水层参数。本文对同步反演模型的框架组成（包括污染物迁移模型和反演优化模型）进行分析后，在对已有的各 种和声搜索改进算法进行研究的基础上结合同步反演模型提出一种改进的和声搜索算法，最后将同步反演模型和 改进的和声搜索算法应用于具体的算例研究。研究表明，改进的和声搜索算法具有算法稳定高效、求解精度高等 特点，能够广泛应用于复杂的地下水污染源和含水层参数反演问题。 关键词：同步反演；污染源反演；含水层参数反演；和声搜索算法；地下水逆问题 中图分类号：P641.2 文献标识码：A 1 研究背景 目前，我国地下水污染已十分严重。与地表水体相比，地下水一旦遭受污染其治理和修复往往 需要更长的时间（常常持续30年以上），其治理和修复费用也极为昂贵。通常来说，要对地下水进行 有效的治理和修复，需要先大致弄清污染源的情况，包括污染源分布以及污染浓度排放情况等，该 ［1］ 过程被称为地下水污染源反演 。一般所说的地下水污染源反演过程，首先是假定含水层参数和边界 ［2-6］ 条件为已知，然后在此基础上根据有限且离散的实测数据进行污染源反演 。然而，在实际问题中 含水层参数（渗透系数、弥散度等）等先决信息往往是未知的，此时就需要根据已有的实测数据对地 ［7-8］ 下水污染源和含水层参数进行同步反演 。 地下水污染源反演问题和含水层参数反演问题都属于地下水逆问题，而最优化方法是求解地下 ［9］ 水逆问题的有力工具 。传统的优化方法有Newton法、单纯形法、共轭梯度法、Hooke-Jeeves方法 ［10］ 等，上述方法均属于局部邻域搜索算法，因而易陷入局部最优而难以获得全局最优解 。近年来， 能搜索到全局最优解的启发式搜索算法（遗传算法、模拟退火算法、禁忌搜索算法、粒子群算法等） 正受到人们越来越多的青睐，并得到了成功的应用［11-14］。 ［15］ 本文采用的和声搜索算法（HarmonySearchalgorithm，HS）是一种全新的启发式全局搜索算法 。 ［16］ Ayvaz 采用HS算法结合地下水流模型MODFLOW对二维地下水管理模型进行了求解，并与GA（Ge⁃ neticAlgorithm，GA）、SA（SimulatedAnnealing，SA）、GASAPF（GeneticAlgorithmbasedSimulatedAn⁃ nealingPenaltyFunction，GASAPF）以及SCE-UA（ShuffledComplexEvolution，SCE-UA）等算法进行比 ［5］ 较，H