井地大地电磁非线性共轭梯度二维反演研究

本文提出并验证了地下大地电磁测深数据的正、反演算法和程序,为井地大地电磁测深和井地可控源音频大地电磁测深这一前沿研究课题提供了地球物理解释方法.使其更加完善和有效.在Rodi(2001)发表的大地电磁场二维NLCG反演理论和算法的基础上,推导了井地大地电磁场二维NLCG反演算法,改编反演的代码,开发出二维NLCG反演程...     

利用一维正则化反演进行大地电磁测深数据拟二维反演解释

大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。     

大地电磁与地震反射波走时数据二维联合反演

由于地球物理反演存在固有的多解性问题,如何减小多解性,提高地球物理方法反演结果的可靠性和准确性是当今地球物理领域研究的热点。为了减小大地电磁与地震反射波走时资料反演的多解性,在Gallardo和Meju研究的交叉梯度理论的基础上,实现大地电磁与地震反射波走时数据的二维联合反演算法;设计理论地电模型和速度模型进行合成数据的单独反演和联合反演试算,并验证大地电磁和地震反射波走时数据二维联合反演算法的正确性和有效性。反演结果表明,利用交叉梯度方法耦合大地电磁电阻率参数和地震波速度参数能使联合反演过程中的这两个模型在更新的时候相互约束,联合反演结果比单一反演结果更接近真实模型,从而使反演结果更可靠,减少了反演解的多解性。     

层状介质大地电磁的自适应量子遗传反演法

将量子遗传算法引入到层状介质大地电磁数据的反演, 得到了层状介质的大地电磁量子遗传反演法.数值试验结果发现该算法仍然存在较严重的早熟收敛现象.为此, 将自适应思想引入到量子遗传算法中来, 通过动态调整量子遗传算法的模型搜索空间, 建立了一种新的改进型量子遗传算法——自适应量子遗传算法, 使算法在迭代过程中能自适应地寻找模型最优值.通过典型测试函数和层状介质大地电磁模型数值试验, 结果表明, 改进算法有效压制了常规量子遗传算法的早熟收敛性, 提高了算法的搜索效率和反演效果.采用该算法对实际的大地电磁资料进行了处理, 取得了较好的地质效果.      

基于遗传算法的地震面波和大地电磁数据一维联合反演算法研究

利用快速标量传递算法进行了一维层状模型的地震面波频散曲线正演计算,得到了基阶频散曲线;利用传统方法对一维层状模型进行了大地电磁正演计算,得到大地电磁的视电阻率曲线数据和相位数据。利用遗传算法对合成数据分别进行了单独反演和联合反演,得到拟合差分布图来分析两种数据的兼容性并研究联合反演的可行性。通过对比单独反演和联合反演的结果,可以看出地震面波频散曲线和大地电磁数据的联合反演结果较之单方法的反演结果要更加精确,收敛速度更快,两种数据的耦合情况较好,同时联合反演能够有效地限制多解性。     

航空大地电磁响应特征分析及反演研究

航空大地电磁法以其机动灵活、效率高等优点,适用于在地势崎岖的偏远环境中开展快速普查作业。在航空大地电磁探测中,起伏地形会对倾子响应产生影响,另外由于倾子是无量纲传输函数,只包含电性结构的空间变化信息,而与电阻率的数值水平无关,单独利用倾子数据进行反演时,反演多解性更大。为解决上述问题,首先采用四边形网格模拟地形起伏;其次在倾子数据的反演目标函数中加入少量地面大地电磁数据拟合项,克服只拟合倾子的多解性缺陷;最后利用二维地电模型,验证本文有限元算法的正演精度,并分析起伏地形对倾子响应数据的影响。针对倾子数据缺乏背景电阻率信息这一问题,开展了带地形的倾子反演和倾子与地面大地电磁数据反演研究。研究结果表明:向倾子数据中加入少量地面大地电磁数据进行反演,可有效降低单独倾子反演对初始模型的依赖。     

国外大地电磁研究现状

简要介绍了目前国外大地电磁的研究现状,重点对二维大地电磁反演研究现状进行了总结,并列举了当前国际流行的常用反演方法,讨论了大地电磁反演解释的发展方向。     

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

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

大地电磁交错采样有限差分二维正反演研究

交错采样网格能自动保证电磁场分布遵守能量守恒定律。本文基于交错采样网格,推导了大地电磁二维有限差分正演过程,实现了二维正演程序;通过与一维解析解对比,验证了算法的正确性且具有较高的计算精度。随后利用有限内存拟牛顿最优化算法,实现了交错采样有限差分二维反演;通过理论模型反演,验证了反演算法的稳定性,揭示有限内存拟牛顿反演计算效率明显优于非线性共轭梯度法。最后,通过对宕昌县官鹅沟大地电磁资料反演解释,查明了测区深部构造特征,表明算法具有较强实用性。     

基于电阻率测井曲线的大地电磁测深标定

大地电磁测深曲线畸变的校正是大地电磁勘探中的重要问题。通过在D盆地的电阻率测井曲线与大地电磁测深平移方式校正前后的反演测深曲线相比较,构造变换函数,实现了测井曲线对大地电磁测深曲线的标定。经过该技术改进的大地电磁反演结果与验证钻孔极为吻合,结合钻井、地震、重力、磁法资料进行综合分析,有效地揭示了D盆地的构造特征。     

基于差分进化的大地电磁反演

针对大地电磁线性迭代反演依赖初始模型,易陷入局部最优解的特点,设计一种基于差异进化的非线性反演方法。利用该方法对大地电磁一维层状介质G、H、HA型地电模型进行反演研究,在无噪声情况下反演结果和模型一致;在加入10%和20%噪声后,反演仍取得良好效果。数值实验结果表明,该反演方法不依赖于初始模型,具有较好的全局优化能力和抗噪声能力,能有效反演大地电磁数据。     

海洋CSEM和MT一维联合反演研究

传统频率域海洋可控源电磁(MCSEM)法利用极化方向电场实部与虚部信息进行反演,由于海底沉积层电阻率一般较低,在海底浅部能取得较为理想的反演效果,但难以探测到海底几千米下的深部基底。这里利用海洋可控源电磁和海洋大地电磁(MT)一维联合反演程序,采用奥克姆反演算法,尝试对海洋电阻率模型进行联合反演。通过反演试算发现,利用频率更低的海洋MT加入反演,不会影响浅部反演效果,同时能探测到深部基底,为海底油气勘探提供了更多技术手段。     

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

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

音频大地电磁法在泸沽湖地区地热勘查中的应用

利用音频大地电磁测深（AMT）在泸沽湖地区进行地热资源勘查,在分析地热形成条件的基础上,通过音频大地电磁测深和大地电磁测深（MT）数据联合反演,结合地热资源赋存有关的水文地质条件,对反演图件进行推断解释,查明了区内主要断裂构造位置、产状、规模及地下热储层深度,取得了较好的应用效果。     

基于边缘计算的大地电磁法三维反演分布式计算实现

随着大地电磁法三维勘探及处理技术的逐步推广,导致数据计算量成倍增长,需要具备高性能计算能力的系统才能满足计算需要。提出基于边缘计算的理念,采用容器技术构建具备大规模计算资源池的分布式计算环境,实现了大地电磁法三维反演的分布式并行计算。基于松耦合的软件架构设计研发了大地电磁法三维反演分布式计算软件,实现了人机交互操作和反演结果可视化。     

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

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

CSAMT拟二维反演并行计算程序设计

可控源音频大地电磁法（CSAMT）拟二维反演并行计算程序,基于多进程与多线程方式实现程序的并行化,并形成相应的算法和程序。采用该并行计算程序,在多处理器多核工作站上进行实测数据的计算测试,结果表明,CSAMT拟二维反演并行计算程序可以有效提高反演效率。     

可控源音频大地电磁法方位非各向同性层状介质Occam反演

目前,可控源音频大地电磁法测深资料反演主要基于各向同性介质,在复杂的地电体中,由于岩石裂隙或矿物形状等原因使得在各个方向的电阻率不同。本文主要探讨更为符合实际地质体情况的方位非各向同性层状介质在正交水平电偶源激发条件下的可控源音频大地电磁数据的反演。运用Occam反演方法,实现了可控源音频大地电磁法方位非各向同性层状介质反演。构建理论地电模型,通过理论模型合成数据的反演算例检验了所实现的可控源音频大地电磁法方位非各向同性层状介质反演算法的有效性和稳定性。同时,不同数据的反演结果表明:方位非各向同性介质的反演应同时使用视电阻率ρ_(sxy)、ρ_(syx)和相位φ_(xy)、φ_(yx)数据。     

EH-4大地电磁法探测地下暗河

通过一个实际工程勘查实例,介绍EH-4大地电磁法探测地下暗河的应用效果.给出电阻率二维反演断面与地质勘探剖面的对应情况,并指出EH-4大地电磁法反演工作中存在的主要问题.     

Transformer网络在大地电磁反演成像中的应用

传统大地电磁反演通常是基于确定性梯度的迭代求解,不仅需要大量时间计算雅可比矩阵,还依赖于初始模型的输入和正则化因子等参数的设置。近年来学者们不断引入机器学习方法以试图改善大地电磁反演,该方法不需要计算雅可比矩阵,不用输入初始模型,训练好的网络仅需几毫秒就可实现反演成像。这里利用Google团队提出的Transformer神经网络经典框架搭建大地电磁数据和模型之间的映射网络,以9240组正演数据为样本,对Transformer网络参数进行训练。采用南非开源大地电磁数据,实现了由视电阻率图像到电阻率模型的反演成像。研究表明:①经训练后的Transformer网络可以较准确的反映出异常体位置和大小;②网络实现了简单的矩阵并行化运算,大幅度提高训练的效率,且成像效率高于传统的反演。