模拟退火算法原理及应用.pptVIP

二 MT-重力联合反演研究 联合反演研究现状： 线性反演算法居多、非线性算法少，用模拟退火算法进行联合反演研究更少。 尽管目前开展的联合反演研究已有多种，但研究内容主要集中在地震与重力、地震与MT联合反演的研究。 有关电磁测深与重力的联合反演研究只查阅到一篇相关的论文，因此对这方面的研究基本上还是空白。 1.MT与重力联合可以使两方法相互弥补 MT纵向分辨率与重力横向的分辨率的互补 实际工作中MT方法的测点点距一般较稀，而野外重力数据的采集点较密。 MT-重力联合反演必要性: 2.充分利用野外资料：在生产实际中，非地震方法一般同时开展。 MT-重力联合反演可能性: 电性与密度同源界面是两种方法联合的前提 地下结构有可能以同源或部分同源的形式出现 示例: 电性与密度界面有几个是一致的： 上第三系-下第三系的物性界面 下第三系-白垩系的物性界面 侏罗系-三叠系物性界面等。 开展MT与重力两者的联合反演是可能的。 MT与重力联合反演技术难点： 如何处理电性与密度界面的不一致情况是首先要遇到的技术难点 如何构筑一个共同的目标函数—因为地下界面的变化使MT与重力两者的场值变化幅度是不同的 MT的正演可以写为以下变分问题： 经比较，本文选用有限元法作为二维MT正演方法，既保持较高的计算精度，又适应于复杂结构的地电模型。 电性与密度界面的GIS编码技术 若是共同的物性界面，MT与重力均需要作正演计算 若是单独的电性界面，只参加MT正演计算 若是单独的密度界面，只参加重力正演计算 对考虑三种情况： 由此产生的三个问题： 如何记录由MT与重力合成的模型，这个模型中有两者一致和不一致的物性界面 如何扰动由两种不同地球物理方法合成的模型 如何让扰动后模型分别参加MT与重力正演计算以求得有关的目标函数 解决办法：借鉴GIS中的属性编码思想 记录模型数据的MTG(i,j),按照MVFSA算法的具体要求进行全局或在当前模型周围扰动产生新模型 编码数组Code(k)（k=1,2…N)，对记录电性与密度界面情况的数组MTG(i,j)进行界面属性编码 各自的正演计算 具体的界面属性编码方法： 电性与密度界面一致：Code(k)＝0，MT与重力均作正演计算； 只是电性界面：Code(k)＝1，只有MT作正演计算； 只是密度界面：Code(k)＝－1，只有重力作正演计算； 目标函数及其加权联合： WMT+Wg=1，体现了联合后的目标函数与Φmtφg的尺度相同，以防止目标函数的发散。 WMT、Wg的大小决定了何者为主反演方法，相应的权系数越大，则该方法为主反演方法。 落实到WMT、Wg的具体取值则是一个非常复杂的问题。 背斜模型的反演结果 二层模型的反演结果 三层模型反演结果 电阻率拟合曲线 三层模型相位拟合曲线 三层模型目标函数与迭代次数关系曲线 结 论 GIS技术可为MT-重力模拟退火联合反演建立初始模型。 改进的MVFSA算法，计算速度比VFSA算法快20-30%；在同等条件下，MVFSA算法更稳健。 针对MT曲线与重力场共源与非共源的复杂情况，借鉴GIS中的属性编码技术，对电性与密度界面的一致性与否进行编码，简洁而有效地处理了两者的正演计算。 无论是重力模拟退火反演还是MT-重力模拟退火联合反演，在模型试验取得成功的基础上，对实际资料进行了处理，地震与钻井等有关资料验证了反演结果的合理性，表明本文开发的计算软件均可应用于实际资料的处理。 借鉴了GIS的可视化基本思想，应用Surfer软件成功实现了综合物探反演成果的剖面、平面及三维曲面的可视化。 建 议 SA算法的模型扰动方式及相应的目标函数计算，约束条件的加入等方面尚有进一步改进的余地 GIS与综合地球物理研究的结合有相当大的空间，需要进一步去开拓，从而深化地球物理的综合研究 MT与重力曲线的变化十分复杂，对如何更快速有效地进行MT与重力的联合应进行深入研究。 多谢指正！ * 模拟退火算法原理及应用研究 主讲: 陈 华 根 同济大学海洋与地球科学学院 一 模拟退火算法及VFSA算法 模拟退火算法在反演中的应用： 非线性组合优化算法：模型扰动，模拟退火，全局寻优。 能量函数—目标函数 模拟退火过程—反演迭代 传统模拟退火流程图 Yes No 随机选择初始模型m0 计算能量函数E(m0) 模型扰动产生新模型 m1=m0+△m0 计算能量函数E(m1) △E= E(m1) -E(m0) △E≤0? m0= m1 新模型按 Metropolis准则接受 缓慢降低温度 满足收敛条件为止 Metropolis接受准则： ⊿E0，新模型无条件被接受—接受能量值较小状态 否则， 产生随机数ξ∈[0，1] 若rξ，新模型被接受，否则被舍弃。 —接受能量值