基于模拟退火的全局混合反演方法及其应用.PDF

第２８卷 第３期 桂 林 工 学 院 学 报 Ｖｏｌ２８Ｎｏ３ ２００８年８月 ＪｏｕｒｎａｌｏｆＧｕｉｌｉｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ Ａｕｇ．　２００８  文章编号：１００６－５４４Ｘ（２００８）０３－０３１２－０７ 基于模拟退火的全局混合反演方法及其应用 １ ２ １ １，２ 刘海飞 ，阮百尧 ，张赛民 ，吕玉增 （１中南大学 信息物理工程学院，长沙 ４１００８３；２桂林工学院资源与环境工程系，广西 桂林 ５４１００４） 摘　要：为了增强模拟退火算法的局部搜索性能，将局部搜索能力较强的单纯形法和鲍威尔法引入模 拟退火算法当中，形成了基于模拟退火的全局混合反演方法。给出了模拟退火算法与单纯形法的镶嵌 式结合方式以及模拟退火算法与鲍威尔法的串行式结合方式，并对其结合步骤进行了详细的阐述，以 便于编程实现。分别对模拟和实测的直流激电测深曲线进行了反演试算：对于不具有等值性的激电测 深曲线，反演效果很好，模型参数的相对误差基本控制在 １％以内；对于等值性较严重或含有观测误 差的激电测深曲线，全局混合反演方法也不能完全克服多解性，但可以通过改变模型参数空间和算法 参数进行多次反演，最终可给出合理的解估计。 关键词：模拟退火；单纯形法；鲍威尔法；全局混合反演 中图分类号：Ｐ６３１　　　　　　　　　　　　　　　文献标志码：Ａ 地球物理中的反问题大多都具有多参数、多 部直接优化方法具有收敛速度快的特点，并且它 极值和非线性的特点，线性或拟线性反演方法往 们都不需要计算偏导数矩阵，那么，就可将两者 往只能得到局部最优解，并且反演结果常常与初 有机的结合起来，起到取长补短的作用，笔者将 始模型有关，如果初始模型选择的合适，它可以 其称为全局混合反演方法。那么作为混合反演方 得到较好的结果，反之，可能会造成反演的假象。 法，如何根据不同算法的特点将它们有机的结合 另外，广义线性反演通常要用到模拟数据对模型 起来，以改善算法的整体结构和搜索机制，最终 参数的导数值，对于许多实际问题来说，计算导 达到提高算法的整体优化性能，这是本文开展混 数值有些困难，而以随机过程为基础的模拟退火 合反演方法研究的目的。本文将局部搜索能力较 算法不受确定规则的限制，自由的在模型空间内 强的单纯形法和鲍威尔法引入到模拟退火算法当 进行随机搜索，只要冷却进度表参数设置合理， 中，形成基于模拟退火的全局混合反演方法。最 可获得全局最优解或次最优解。近几年随着模拟 后，将其应用到直流激电测深反演当中，以检验 退火算法研究的不断深入，它已经越来越广泛的 算法的有效性。 应用于自然科学的各个领域，如计算机设计、图 ［１－３］ １　混合反演方法的结合方式 像处理、信号处理、地球物理等方面 。但在 实际应用当中，人们发现模拟退火算法往往只能 谈及多种方法的结合就涉及到结合方式问题， 得到全局次最优解而不是最优解，这主要与冷却 在这里给出了混合算法的串行和镶嵌结合方式： 进度表的参数设置有关。为减少模拟退火算法对 串行结构是混合算法中最简单的一种结合方式， 冷却进度表要求的苛刻性以及增强其局部搜索性 如图１ａ。它吸取不同算法的优点，用一种算法的 能，考虑将模拟退火与其他局部直接优化方法结 搜索结果作为另一种算法的初始解依次对问题进 其目的主要是在保证次全局最优解具有 合起来。由于模拟退火法全局搜索能力强，而局 行优化， 　收稿日期：２００７－０８－３０ 　基金项目：广西科技开发项目 （桂科攻０４２８００６－１７）；广西教育厅科研项目 （２００５０７０４４） 　作者简介：刘海飞 （１９７５－），男，博士，讲师，主要从事电磁法的资料处理及其正反演研究，Ｅｍａｉｌ：ｌｉｕｈａｉｆｅｉ＠１２６ｃｏｍ。 第３期　　　　　　　　　　　　刘海飞等：基于模拟退火的全局混合反演方法及其应用