地下水渗流模型参数识别的模拟退火算法  

PARAMETER IDENTIFICATION PROCEDURE FOR GROUNDWATER FLOW MODEL WITH SIMULATED ANNEALING

作　　者：李守巨[1] 上官子昌[2,3] 刘迎曦[1] 于贺[1] 

机构地区：[1]大连理工大学工业装备结构分析国家重点实验室,辽宁大连116024 [2]大连理工大学土木水利学院,辽宁大连116024 [3]大连水产学院土木工程学院,辽宁大连116023

出　　处：《岩石力学与工程学报》2005年第A01期5031-5036,共6页Chinese Journal of Rock Mechanics and Engineering

基　　金：国家自然科学基金资助项目（10072014）

摘　　要：反问题的求解常常需要转化为非线性优化问题，其目标函数定义为观测数据与模型数据之间的残差平方和。地下水模型参数识别最常用的优化方法都是基于梯度搜索，其缺陷在于对模型参数初始估计比较敏感和局部极小问题。与传统的基于梯度搜索的优化方法相比，模拟退火算法具有良好的全局收敛特性。把含水层参数识别反问题转化为组合优化问题，提出模拟退火算法识别二维、非稳态地下水流动模型的渗透系数和储水系数的策略。反问题的不适定性由解的不唯一性和不稳定性来表征，模拟退火算法具有解决这一问题的能力。通过与梯度搜索算法相对比，数值模拟计算结果显示所提出反演方法的有效性和适用性。The solution of inverse problem usually requires nonlinear optimization of an objective function describing the difference between measured and simulated data. Most optimization algorithms used for parameter estimation in groundwater hydrology are gradient-type methods that have the disadvantages of being very sensitive to the initial guesses of parameters and being prone to converge to local minima. Compared with traditional optimization algorithms, simulated annealing algorithm is recognized to have better capability to find the global optimal solution. The inverse problem of identifying aquifer parameters is treated as a combinational optimization problem. The simulated annealing is presented to identify the transmissivity and storage coefficient for a two-dimensional unsteady state groundwater flow model. The ill-posedness of the inverse problem as characterized by instability and non-uniqueness is overcome by using simulated annealing algorithm. The effectiveness and flexibility of presented inversion technique are evaluated and compared with descent search methods.

关 键 词：水力学 参数识别 全局优化 模拟退火算法 导水系数 地下水模型 

分 类 号：TU452[建筑科学—岩土工程] TV13[建筑科学—土工工程]