地球物理资料非线性反演方法讲座(三) 模拟退火法Simulated Annealing Method
师学明;王家映;
摘要(Abstract):
模拟退火法源于统计热力物理学,它模拟熔融状态下物体缓慢冷却达到结晶状态的物理过程。模拟退火反演算法的基本思想是:生成一系列参数向量模拟粒子的热运动,通过缓慢地减小一个模拟温度的控制参数,使模拟的系统最终冷却结晶达到系统能量最小值的过程。模拟退火反演算法实质是利用了地球物理反演问题求解过程与熔化固体退火过程的相似性,开辟了地球物理反演的新途径,是非线性反演算法中一种最常用的算法。本讲座概要地介绍了模拟退火法的基本原理,模型搜索及解的接受准则,模拟退火法的分类及实现方法,并给出了模拟退火法在地震资料、电法资料及其他地球物理资料反演中的实例,最后总结和归纳了拟退火法的特点以及该方法的局限性。
关键词(KeyWords): 模拟退火;非线性;反演;地球物理
基金项目(Foundation): 国家自然科学基金(编号:40274039,49674227,40204007)资助
作者(Authors): 师学明;王家映;
参考文献(References):
- [1]王家映.地球物理反演理论[M].北京:高教出版社,2002.
- [2]杨文采.线性地球物理反演方法回顾与展望[J].地球物理学进展,2002,17(2):255~261.
- [3]Metropolis N,Rosenbluth A,Rosenbluth M,et.al.Equation of state calculations by fast computing ma-chines[J].J.Chem.Phys.,1953,21:1087~1092.
- [4]Kirkpatrick S,Celatt C D,Vecchi MP.Opti mization by si mulated annealing[J].Sciences,1983,220:671~680.
- [5]Rothman D H.Nonlinear inversion,statistical me-chanics,and residual statics esti mation[J].Geophys-ics,1985,50(12):2784~2796.
- [6]Rothman D H.Automatic esti mation of large residual statics correction[J].Geophysics,1986,51(2):337~346.
- [7]Sen MK,and Stoffa P L.Nonlinear one-di mensional seismic waveforminversion using si mulated annealing[J].Geophysics,1991,56(10):1624~1638.
- [8]Send M K,Bhattacharya B B,Stoffa P L.Nonlinear inversion of resistivity sounding data[J].Geophys-ics,1993,58(4):496~507.
- [9]师学明,王家映.一维层状介质大地电磁模拟退火反演法[J].地球科学—中国地质大学学报,1998,23(5):543~545.
- [10]于鹏,王家林,吴健生,等.重力与地震资料的模拟退火约束联合反演[J].地球物理学报,2007,50(2):529~538.
- [11]Ingber L.Very fast si mulated annealing[J].Mathe-matical and Computer Modeling,1989,12(8):967~973.
- [12]王家映.地球物理反演问题概述[J].工程地球物理学报,2007,4(1):1~3.
- [13]王家映.蒙特卡洛法[J].工程地球物理学报,2007,4(2):81~85.