2008-05-23 基于量子蒙特卡罗的地球物理反演方法.pdfVIP

书第５１卷 第５期 ２００８年９月 地　球　物　理　学　报 ＣＨＩＮＥＳＥ　ＪＯＵＲＮＡＬ　ＯＦ　ＧＥＯＰＨＹＳＩＣＳ Ｖｏｌ．５１，Ｎｏ．５ Ｓｅｐ．，２００８ 魏　超，李小凡，张美根．基于量子蒙特卡罗的地球物理反演方法．地球物理学报，２００８，５１（５）：１４９４～１５０２ ＷｅｉＣ，ＬｉＸＦ，Ｚｈａｎｇ Ｍ Ｇ．ＴｈｅｇｅｏｐｈｙｓｉｃａｌｉｎｖｅｒｓｅｍｅｔｈｏｄｂａｓｅｄｏｎＱｕａｎｔｕｍ ＭｏｎｔｅＣａｒｌｏ．犆犺犻狀犲狊犲犑．犌犲狅狆犺狔狊．（ｉｎ Ｃｈｉｎｅｓｅ），２００８，５１（５）：１４９４～１５０２ 基于量子蒙特卡罗的地球物理反演方法 魏　超，李小凡，张美根 中国科学院地质与地球物理研究所，北京　１０００２９ 摘　要　本文将量子蒙特卡罗全局优化方法引入地球物理反问题，进而发展了一类新的地球物理非线性反演方 法．量子蒙特卡罗方法是基于量子力学机制的随机方法，包括变分蒙特卡罗方法、格林函数蒙特卡罗方法、扩散蒙 特卡罗方法、路径积分蒙特卡罗方法．本文简要回顾了量子蒙特卡罗方法的发展，阐述了其方法理论；随后的数值 试验结果表明，量子蒙特卡罗方法应用于地球物理反问题的求解是成功的，它适合于非线性、多极值的地球物理反 演问题，在收敛速度和避免陷入局部极小等方面有着一定的优势，且该方法也适用于其他领域非线性最优化问题 的求解，其算法具有较强的通用性；最后就量子蒙特卡罗方法在地球物理反问题中的应用前景以及存在的问题做 了简要概述． 关键词　量子蒙特卡罗，路径积分蒙特卡罗，变分蒙特卡罗，格林函数蒙特卡罗，扩散蒙特卡罗，地球物理反演 文章编号　０００１５７３３（２００８）０５１４９４０９ 中图分类号　Ｐ６３１ 收稿日期２００７１１０６，２００８０６２７收修定稿 基金项目　国家自然科学基金项目４０３０４００４，２００７ＣＢ２０９６０３）资助． 作者简介　魏　超，男，１９８１年生，中国科学院地质与地球物理研究所博士研究生，主要从事地球物理反演方面的研究． Ｅｍａｉｌ：ｃｗｅｉ＠ｍａｉｌ．ｉｇｇｃａｓ．ａｃ．ｃｎ 犜犺犲犵犲狅狆犺狔狊犻犮犪犾犻狀狏犲狉狊犲犿犲狋犺狅犱犫犪狊犲犱狅狀犙狌犪狀狋狌犿 犕狅狀狋犲犆犪狉犾狅 ＷＥＩＣｈａｏ，ＬＩＸｉａｏＦａｎ，ＺＨＡＮＧ ＭｅｉＧｅｎ 犐狀狊狋犻狋狌狋犲狅犳犌犲狅犾狅犵狔犪狀犱犌犲狅狆犺狔狊犻犮狊，犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犛犮犻犲狀犮犲狊，犅犲犻犼犻狀犵１０００２９，犆犺犻狀犪 犃犫狊狋狉犪犮狋　Ｑｕａｎｔｕｍ ＭｏｎｔｅＣａｒｌｏｍｅｔｈｏｄ（ＱＭＣ），ｉｎｃｌｕｄｉｎｇｖａｒｉａｔｉｏｎａｌＭｏｎｔｅＣａｒｌｏ（ＶＭＣ）， Ｇｒｅｅｎ′ｓｆｕｎｃｔｉｏｎＭｏｎｔｅＣａｒｌｏ（ＧＦＭＣ），ｄｉｆｆｕｓｉｏｎＭｏｎｔｅＣａｒｌｏ（ＤＭＣ）ａｎｄｐａｔｈｉｎｔｅｇｒａｌＭｏｎｔｅ Ｃａｒｌｏ（ＰＩＭＣ），ｉｓａｓｔｏｃｈａｓｔｉｃｍｅｔｈｏｄｂａｓｅｄｏｎｔｈｅｑｕａｎｔｕｍ ｍｅｃｈａｎｉｃｓ．Ｔｈｉｓｐａｐｅｒ，ｆｉｒｓｔｌｙ， ｇｉｖｅｓａｂｒｉｅｆｓｕｍｍａｒｉｚａｔｉｏｎａｂｏｕｔｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｑｕａｎｔｕｍ Ｍｏｎｔｅ Ｃａｒｌｏ ｍｅｔｈｏｄａｎｄ ｅｘｐａｔｉａｔｅｓｉｔｓｐｒｉｎｃｉｐｌｅｓ；ｔｈｅｎｄｏｅｓｓｏｍｅｍｏｄｅｌｔｅｓｔｓ，ｗｈｉｃｈｓｈｏｗｔｈａｔｔｈｅｕｓｉｎｇｏｆｑｕａｎｔｕｍ ＭｏｎｔｅＣａｒｌｏｍｅｔｈｏｄｆｏｒｇｅｏｐｈｙｓｉｃａｌｉｎｖｅｒｓｅｐｒｏｂｌｅｍｓｉｓｓｕｃｃｅｓｓｆｕｌ，ａｎｄｉｔｓａｌｇｏｒｉｔｈｍｉｓａｌｓｏ ａｃｃｅｐｔａｂｌｅｆｏｒｓｏｍｅｏｔｈｅｒｎｏｎｌｉｎｅａｒｏｐｔｉｍｉｚａｔｉｏｎｐｒｏｂｌｅｍｓ；ｌａｓｔｌｙ，ｇｉｖｅｓａｐｅｒｓｐｅｃｔｉｖｅｏｎｔｈｅｕｓｉｎｇ ｏｆｑｕａｎｔｕｍＭｏｎｔｅＣａｒｌｏｍｅｔｈｏｄｆｏｒｇｅｏｐｈｙｓｉｃａｌｉｎｖｅｒｓｉｏｎ，ａｎｄｐｏｉｎｔｓｏｕｔｉｔｓｓｈｏｒｔｃｏｍｉｎｇｓ． 犓犲狔狑狅狉犱狊　ＱｕａｎｔｕｍＭｏｎｔｅＣａｒｌｏ，ＰａｔｈｉｎｔｅｇｒａｌＭｏｎｔｅＣａｒｌｏ，ＶａｒｉａｔｉｏｎａｌＭｏｎｔｅＣａｒｌｏ，Ｇｒｅｅｎ′ｓ ｆｕｎｃｔｉｏｎＭｏｎｔｅＣａｒｌｏ，ＤｉｆｆｕｓｉｏｎＭｏｎｔｅＣａｒｌｏ，Ｇｅｏｐｈｙｓｉｃａｌｉｎｖｅｒｓｉｏｎ １　引　言 人类对地球内部物理性质以及矿产资源分布的 了解，大多来自地表地质和地球物理、地球化学资料 的反演和解释［１］．在科学研究及工程技术中，凡涉及 观测数据反推模型或体系时，均无法回避反演问题． ２０世纪六七十年代，人们把地球物理反演问题近似 为线性问题，采用线性化近似方法求解，如广义反演 理论；而后发展的非线性反演方法则是将反演问题 　５期 魏　超等