自适应全局混合反演3.pdf

中 国 科 学 (D 　辑) 第 28 卷 　第 2 期 SCIENCE IN CHINA ( Series D) 　 1998 年 4 月 自适应全局混合反演 艾印双　刘鹏程　郑天愉 ( 中国科学院地球物理研究所, 北京 100101) 摘要 　　地球物理学反演多是一个多参数、非线性、多极值、非连续函数的优化问题. 在模拟退火与单纯形、均匀设计中引入自适应的退火规则 ,形成自适应全局混合反演 方法. 数值检验与模型计算表明 ,该方法具有搜索速度快、效率高的特点. 关键词 　　模拟退火 　单纯形 　均匀设计 　自适应 　全局混合反演 地球物理学的反问题大多是一个多参数、非线性、多极值、非连续函数的优化问题. 对于 非线性反演问题 ,用局部方法反演的结果容易陷入局部极值解 ,而非全局最优解. 所以 ,用全 局反演方法来求此类问题的最优解就显得越来越重要. 在全局反演方法中 ,最典型的算法为 ( ) [1～5 ] [6 ] 模拟退火方法 SA . Liu 等 将全局方法搜索能力强、局部方法收敛速度快的特点结合 起来 ,形成全局混合反演方法 ; 纪晨、姚振兴[7 ] 将均匀设计方法引入到地球物理反演问题中. 这些方法极大改善了地球物理反演中的搜索能力和速度问题. 在大多数具有 D 维参数空间 的物理问题中 ,不同的参数具有不同的取值范围 ,具有不同的退火敏感度 , Ingber [1 ] 通过加大 不敏感参数的搜索范围 ,形成了 VFSR (Very fast simulated reannealing) 方法[1 ] . 本文在综合 模拟退火、单纯形和均匀设计的全局混合反演方法中 ,进一步引入自适应的退火规则 ,形成自 适应的全局混合反演方法 ,文中通过一些实例证明了该方法的有效性. 1 　自适应全局混合反演方法原理 刘鹏程等[8 ] 提出的综合反演方法将模拟退火、单纯形和均匀设计方法结合起来 ,新模型 的生成公式可以表示为 j j Δ j ( ) M = M d + M k , 1 j j j χ Δ kΦ ( ) M k = k , 2 j j χ π( ξ ) k - 0. 5 ( )