大地电磁测深拟二维模拟退火约束反演

大地电磁测深(MT)资料解释的核心是反演问题,其反演方法可分为两大类:局部方法和全局方法。局部方法基于最小方差原理,依赖于目标函数的导数求解,因而其解依赖于初始模型,容易陷入局部极值,难以得到全局最优解。全局方法如模拟退火算法,它是近年发展起来的全局最优化算法。其主要优点是:不用求目标函数的偏导数及解大型矩阵方程组,即能找到一个全局最优解,而且易于加入约束条件,方法易于移植。目前模拟退火算法已开始用于解决非线性地球物理反演问题,但还存在计算效率低的缺陷。以快速模拟退火算法为基础,对大地电磁测深理论数据及实际资料进行快速模拟退火约束多参量的反演,效果较为理想。对H型地电断面的理论曲线反演和根据某盆地预探井实际资料对中生界厚度及基底埋深的反演作了举例说明。     

基于MPI大地电磁二维正则化反演并行计算研究

运用并行技术能提高运算效率,在大地电磁数据处理中,偏导数矩阵的计算占据了大量的计算时间,使得反演速度较慢.这里在目标函数中采用基于阻尼高斯—牛顿法的正则化反演方法,对偏导数计算进行并行处理,其余部份保留串行计算模式,针对不同频率间实现完全对等的并行策略,模拟在多进程下的并行计算速度.通过两个理论模型试算,验证了并行算法的正确性、高效性和稳定性.     

起伏地形大地电磁二维反演

为了适应地形起伏的实际地质情况,开展带地形的最小二乘二维反演研究。鉴于大地电磁（ MT）反演的不适定问题,引入Tikhonov的正则化方法,从而获得关于总目标函数的方程,利用光滑约束最小二乘法求解总目标函数方程。由于正则化因子值与反演精度以及稳定性相关,采用主动约束平衡方法获取最优化的正则化因子,以确保反演精度和稳定性都达到最佳。与此同时,利用电磁场互易定理以节省反演迭代过程求解雅可比矩阵的计算时间。构建了若干地质构造模型进行试算,分别讨论TE、TM模式以及二者联合模式的反演结果,并与前人研究工作对比以说明本文方法的反演效果。     

大地电磁测深资料的超球逼近反演

超球逼近方法是一种全局搜索的二次优化方法。将该优化方法应用到大地电磁测深资料的反演中取得较好的效果。在介绍了超球逼近方法的基本原理后,应用数值结果对该方法的可行性进行了检验。     

采用生物遗传算法的大地电磁测深资料反演

本文用遗传算法（ＧＡ）进行大地电磁测深（ＭＴ）一维模型的反演。本反演方法无需初值，不用求导数，受噪声影响小，能有效克服反演中出现的局部极值，并能给出一组解     

混合遗传算法在大地电磁一维反演中的应用

遗传算法是近些年来产生和发展的一种模拟生物进化过程的自适应启发式全局优化的搜索算法。它不完全依赖于初始猜测,且具有全局收敛的特点,可以被用来解决各种复杂的实际问题,如工程优化设计,人工智能和决策系统,以及地球物理反演等。尽管遗传算法是一种效率很高的全局优化算法,但许多仿真结果表明,它具有计算时间长,局部搜索能力弱的缺点。而共轭梯度法属于非启发式全局优化搜索方法,收敛速度快,但容易陷入局部极值,且严重依赖初始猜测。根据遗传算法和共轭梯度法的特点,这里提出了一种混合遗传算法,用来进行地球物理反演。该算法既具有遗传算法的全局收敛性,又有共轭梯度法的快速收敛性,经实际应用,取得了良好的效果。     

大地电磁测深数据子库

"大地电磁测深数据子库研制"是国土资源部"深部地球物理探测数据共享与对比研究"项目的课题之一.在此介绍课题的研究概况,详细讨论大地电磁测深数据库的需求分析、开发平台的选取、系统结构和功能设计等问题.     

二维大地电磁测深交互式正、反演系统的实现及其应用

目前，大地电磁测深的资料处理与正、反演软件都是FORTRAN语言编写的，由于FORTRAN语言本身的局限性，使得这个软件存在很多的缺陷，但是因为用FORTRAN语言编写的算法已经比较成熟，将它全部用其他语言移植也比较费时费力，为了避免资料的浪费，提高大地电磁测深的资料处理与交互式系统的开发效率，并且更加完善其功能，在开发大地电磁测深二维交互式正反演系统时，对大地电磁测深中的计算部分仍用现有的FORTRAN程序，对有关的界面与图形诉开发部分用VC＋＋语言来编写，以增强和体现其交互式的功能。     

二维大地电磁正则化共轭梯度法反演算法

针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。     

大地电磁测深三维反演系统及应用

大地电磁测深是基于电磁感应原理,在地表采集正交的电、磁场来研究地壳乃至上地幔电性结构的地球物理方法,后期衍生出采集高频天然场信号的音频大地电磁测深和采集人工源信号的可控源大地电磁法。此类方法在深部构造研究、矿产资源勘查、油气资源探测等多领域得到了广泛应用。近几年来,虽然其野外工作趋近三维勘探水平,但其采集数据的处理和三维反演研究仍然存在不足,限制了地质解释的准确度和效率。针对这些问题,我们开发了大地电磁测深的三维反演软件系统,系统主要分为静位移校正、三维反演的初始模型构建和三维反演3个模块,有效地解决了数据处理中的静位移问题和数据解释中的三维反演问题,并且其中涉及的很多算法改进了现有技术,在计算精度、计算效率和成本节约方面有很大优势,具有良好的应用前景。     

大地电磁法的1D无偏差贝叶斯反演

应用贝叶斯理论对一维(1D)大地电磁反演问题进行无偏差不确定度分析。在贝叶斯理论中,测量数据和先验信息包含在后验概率密度函数(PPD)中,它可以解释成模型的单点估计和不确定度等贝叶斯推断,这些信息的获取需要对反演问题进行优化求最优模型和在高维模型空间中对PPD进行采样积分。采样的完全、彻底和效率,对反演结果有着重要的影响。为了使采样更有效、更完全,数值积分采用主分量参数空间的Metropolis Hastings采样,并采用了不同的采样温度。在反演中,同时采用了欠参数化和超参数化方法,数据误差和正则化因子被当成随机变量。反演结果得到各参数的不确定度、参数间的相关关系和不同深度模型的不确定度分布。COPROD1数据的反演结果表明模型空间中存在双峰结构。非地电参数在反演中得到了约束,说明数据本身不仅包含地球物理模型信息(电导率等),还包含了这些非地电参数的信息。     

基于并行处理的CSAMT拟二维反演解释

为提高野外CSAMT测深数据的反演解释效率,提出了将并行计算技术引入CSAMT测深一维自动迭代反演中进行CSAMT拟二维反演解释的思路,基于多处理器多核工作站,结合多进程和多线程并行编程技术,形成相应并行算法和程序。采用该并行计算方法进行实测数据反演解释,结果表明,CSAMT拟二维反演解释并行计算技术能较大幅度地提高CSAMT测深数据的反演解释效率。     

一维层状介质大地电磁模拟退火反演法

大地电磁模拟退火反演法是一种最优化的非线性反演方法,与传统的线性反演方法相比该方法具有：（１）不依赖于初始模型的选择;（２）能寻找全局最小点而不陷入局部极小;（３）在反演过程中不用计算雅可比偏数矩阵等优点;通过对各种类型的大地电磁测深理论曲线试算,结果表明模拟退火法能准确地自动反演地电参数（地层电阻率,厚度）最后对实际资料进行了处理,取得了较好的效果。     

地震—大地电磁数据非线性联合反演方法研究

在地球物理数据广义线性联合反演方法的基础上，研究了地球物理数据非线性联合反演方法。对这两种方法的反演效果进行了比较、分析，得出了非线性联合反演方法比广义线性联合反演方法更优越的结论。还对解决迭代公式中数量级不一致的问题进行了论述     

测高数据不准时的直升机航空瞬变电磁一维反演方法理论研究

航空电磁勘探中,高度计显示的发射线圈高度数据往往是不准确的,这可能会影响反演结果,出现假异常现象.针对这一情况,提出了一种高度计数据不准确条件下的航空瞬变电磁资料一维反演方法,在重构地电参数的同时也对发射线圈高度值进行校正.由于航空电磁反演问题是典型的病态问题,使用正则化反演方法,正则化因子由数据目标函数与模型目标函数之间的关系以自适应方式给出,使得反演迭代能稳定进行.以直升机航空瞬变电磁资料反演为例,用理论模型响应数据加入噪声进行反演,结果表明,无论初始高度值是高于真值还是低于真值,该反演方法均能够较好地重构地下介质参数和发射线圈高度值;作为反演初值的发射线圈高度测量值越准确时,反演迭代收敛速度就越快,反演效果也越好.     

复杂地形条件下的可控源音频大地电磁测深数据二维反演技术及应用效果

复杂地形严重地影响可控源音频大地电磁测深(CSAMT)的解释结果,为此,笔者考虑了地形的影响,将地电断面模型经过精细剖分,对模型正演计算的数据进行CSAMT二维反演,得出的结果与正演模型的形态基本一致,获得的电性断面能很好地反映地下地质结构,与实际地质资料相吻合,证实了该方法的正确性和有效性.     

多功能海底大地电磁数据处理软件包

作者在本文中详细地介绍了“863计划”项目“海底大地电磁数据处理和解释系统研究”所取得的成果。本系统功能强大,既能对海底数据进行畸变改正、压制噪声、二维正反演、三维正演,又能对中国南海、东海典型地段进行正演模拟。     

用改进的佐迪反演方法进行二维电阻率图像重建

针对佐迪反演方法在二维拟断面应用中的缺陷进行了三个方面的改进：固定佐迪反演方法中用于转换电极距和深度关系的深度转换因子,可使绝大多数数字模型的反演过程稳定收敛;在调整视电阻率的公式中加入一个迭代系数,可使计算值快速趋向于观测值,减少了迭代次数;对于含有噪声的数据,加入滤波因子可以在一定程度上抑制噪声。通过数字模拟和野外实测数据对改进后的佐迪反演方法在二维电阻率反演应用中的可行性进行了验证。模拟研究包含了简单模型,也包含一些复合模型,二者均取得较理想的应用效果。     

高温超导TEM瞬变磁场拟二维反演及应用效果

随着高温超导技术的应用及完善,在TEM测量中直接观测瞬变磁场信号已成为现实,而瞬变磁场数据反演解释方法成为当前急需解决的问题.在借鉴感应电动势资料反演解释方法的基础上,全面阐述了基于瞬变磁场的视电阻率定义,研究了TEM瞬变磁场拟二维反演技术,并从理论模型及实测资料两个方面验证了方法技术的实用性,为促进高温超导TEM的发展提供了方法技术支撑.     

瞬变电磁法利用不完整数据进行一维反演的可行性分析

为充分利用已经采集获得的并不完整的瞬变电磁衰减曲线,以一维三层理论模型以及实测数据为例进行了计算,计算结果说明,在不能获得整条衰减曲线时,利用部份时段进行反演获取一定的的地质信息,是一种可行的补救措施。试验的前提是瞬变响应曲线是正常的符合衰减规律的曲线,而对于畸变的曲线是不适合的。