带约束的多辐射场源半航空瞬变电磁一维自适应正则化反演方法

基于多辐射场源半航空瞬变电磁法,研究了一种带约束的适用于垂直分量的一维自适应正则化反演算法.该算法在Occam反演的基础上,采用CMD自适应调节方案改进了拉格朗日乘子的处理方式,利用自然边界条件和模型修正量可行下降方向法对反演过程进行约束,在提高反演过程计算效率的同时,保证反演结果的稳定性和可靠性.层状模型试算结果表明...     

航空瞬变电磁数据一维Occam反演

分析了奥克姆(Occam)反演法的原理和计算方法,并将其应用于时间域航空电磁一维反演中,给出了一种简便的正则化参数搜索方法,对迭代过程进行了改进,大大提高了反演速度,5次迭代约20s。在反演初始模型中,电阻率数据由原始感应电动势数据通过自动迭代法转换而来,深度数据由初始电阻率计算得出。经算例表明,收敛过程中满足拟合差最小的拉格朗日乘子μ值随迭代次数增加而递减;反演过程非常稳定,一般5次迭代至10次迭代就能达到收敛,对低阻层的反演深度和电阻率值都接近真实值,但对高阻层的反演深度和电阻率值误差较大。     

考虑激发极化效应的瞬变电磁一维Occam反演

大地的感应激发极化效应有时会在瞬变电磁晚期响应上产生负值变号现象。常规实电阻率瞬变电磁反演由于没有考虑激发极化效应,对于观测数据负值部分的解释一直存在问题,这影响了反演解释的真实性和可信度。对此,本文首先将Cole-Cole复电阻率模型加入正演计算,并验证了计算的准确性。然后在阻尼最小二乘法的基础上加入Occam平滑约束来构建反演方程,能够同时反演出零频电阻率、充电率、时间常数以及频率相关系数,实现了一维瞬变电磁复电阻率反演算法。最后,建立具有不同程度激发极化效应的典型层状理论模型进行反演试算与结果分析,结果表明,在激发极化效果明显的低阻高极化地层中,复电阻率反演效果更好。与实电阻率反演结果的对比说明,瞬变电磁复电阻率反演既可以达到实电阻率的常规反演效果,也能解决实电阻率无法实现的负值拟合问题。     

多辐射场源半航空瞬变电磁法薄层探测能力

多辐射场源半航空瞬变电磁法具有探测深度大、地形适应性强和工作效率高的优点,已逐渐成为当前的研究热点。为了研究其探测能力,本文基于一维正演理论,从多辐射场源布设方式入手,选定典型测点与辐射源组合,通过定义瞬变响应相对异常,分析对比不同模型薄层电阻率、厚度、埋深等参数对薄层探测的影响。结果表明：根据测点位置合理设置发射电流方向可以有针对性地增大采集信号强度,获得形态相对简单的B和dB场响应曲线;与单辐射源相比,多辐射源激励时瞬变响应不仅幅值更大,衰减速度也更快,对大地电阻率的变化更灵敏;B场更适合用于相对异常计算,探讨多辐射场源半航空瞬变电磁法的薄层探测能力,其中x分量的薄层探测能力最强;薄层电阻率与围岩差异越大、厚度越大、埋深越浅,相对异常曲线的幅值越大,薄层越易于分辨;薄层分辨能力与测点位置有关,当测点位于多辐射源内侧时B场各分量相对异常曲线幅值更大,对薄层分辨能力更强。     

多辐射场源半航空瞬变电磁法多分量响应特征分析

多辐射场源半航空瞬变电磁法是对电性源半航空瞬变电磁法的发展和补充,具有探测深度大、地形适应性强和工作效率高的优点.基于一维正演理论,本文对平行辐射场源、带角度多辐射场源和单个源模型在空中产生的三分量瞬变响应特征进行分析,结果表明:整体来说z分量瞬变响应曲线形态最为简单,多辐射场源情况下响应增强的特点最明显,最利于数据处理和解释;x分量瞬变响应形态最为复杂,在全场域经常表现出极性反转的变号现象,但晚期信号幅值总是增大;y分量瞬变响应受电流方向和接收点位置的影响最大,曲线形态和响应幅值变化均较为复杂.结合偏移距变化下响应衰减幅度,建议工作区域偏移距控制在辐射源长度的5~6倍内为宜.本文的研究结果可为实际勘查工作的野外施工布置及定性分析提供参考,同时获得了对该方法的初步认识.     

半航空瞬变电磁一维聚焦反演研究

半航空瞬变电磁法(SATEM)是一种采用地面发射—空中接收的新兴地球物理勘探方法,具有灵活、高效的特点。目前应用于SATEM的反演方法往往基于最大平滑准则,使反演结果较平滑,无法较好地识别具体的层界面信息。本文将聚焦反演理论引入到半航空瞬变电磁一维反演中,通过选择合适的聚焦因子和正则化因子来确定聚焦反演稳定器,然后对包含聚焦反演稳定器反演目标函数进行求解,使得反演结果更能清晰识别层状地层突变的界面。设置多个层状地电模型验证聚焦反演的可靠性,同时与Occam反演结果进行对比,突出了聚焦反演识别界面的优势。对某地区地下水探测的实际数据进行聚焦反演计算,与水文地质资料、测井信息进行综合分析,查明了该区域地下含水层的位置及其空间展布情况,同时验证了半航空瞬变电磁对地下水探测的可行性。     

大回线源瞬变电磁一维自适应反演方法及应用

提出了一种适用于大回线源瞬变电磁法框内任意测点垂直Z分量的一维反演方法。该算法采用CMD正则化因子优化方案、光滑度函数和自然边界条件对反演过程进行了约束,减小了初值的依赖度,保证了反演过程的自主性和反演结果的稳定性;在反演过程中,在保证模型残差沿梯度下降的同时,对模型修正量的最大值进行了限制,使反演模型不易陷入局部极小,保证了反演结果全局的收敛性;同时反演的灵敏度矩阵采用牛顿法进行更新,减少了整个反演计算的工作量。最后通过五层的双低阻层模型、双高阻层模型和实测资料证明了该算法可靠。     

考虑关断时间效应的瞬变电磁一维反演

为克服因关断时间校正而带来的误差,尝试直接对包含关断时间效应的瞬变电磁资料进行反演。正演采用将斜阶跃波离散为多个阶跃波之和的方法,反演采用经典的马奎特法。对模型试算的结果表明,采用该方法直接对斜阶跃响应进行反演,结果比较理想,计算速度也可以接受。     

基于最小构造模型的地空瞬变电磁一维正反演技术研究

地空瞬变电磁法(Ground-airborne transient electromagnetic method)是一种新兴的物探方法,其探测范围大,工作效率高,在地形复杂的区域有着传统物探方法无法比拟的优势.这里通过对地空瞬变电磁的一维正演问题的研究,采用Hankel变换、Gauss积分和频-时转换实现了地空瞬变电磁...     

Occam反演及其在瞬变电磁测深中的应用

Occam反演技术被应用中心回线瞬变电磁测深数据的反演中.理论模型和实际数据的反演结果表明,利用Occam反演技术可以较为准确地获得地电断面的电阻率分布.一般5次迭代目标函数就可以收敛到5×10-2.     

地空瞬变电磁法研究进展

地–空瞬变电磁法在煤炭采空区勘探等领域受到了越来越多的关注, 有必要研究发射电流波形参数对地–空瞬变电磁三分量响应特征影响, 为地–空瞬变电磁数据处理与解释提供理论依据。以梯形波为例, 首先, 研究不同发射电流波形的频谱分布情况; 然后, 基于三维时域有限差分正演研究发射波形的上升沿时间、脉宽和关断时间对地–空瞬变电磁三分量磁场响应的影响。结果表明：上升沿时间对三分量二次场响应基本不产生影响; 关断时间对三分量二次场响应的影响主要集中在0.2 ms之前, 且关断时间越长对纯异常响应影响越大; 脉宽对三分量二次场响应的影响主要集中在0.1 ms之后, 且脉宽越短对纯异常响应影响越大。三维采空区模型结果表明：关断时间、脉宽对三分量异常场和背景场响应的影响特征基本一致; 通过三分量纯异常场响应多测道图和时间道图可以判断异常体的分布范围和深度。研究成果可为地–空瞬变电磁激励源波形的参数选取提供有价值的理论借鉴。     

回线源瞬变电磁测深快速反演方法

回线源瞬变电磁测深法(简称瞬变电磁法或TEM法)是以不接地回线为场源,在回线内、外测量电磁场的一种电磁测深方法,已经广泛用于固体矿产、工程地质和水文地质调查工作之中[1,2]。目前,回线源瞬变电磁法主要采用中心回线和重叠回线两种测量装置,测量回线内垂直磁场产生的感应电动势。“烟圈”反演方法是瞬变电磁反演的主要计算方法,“烟圈”反演利用地表瞬变电磁响应与某时刻电流源镜像等效原理,计算勘探深度和对应电阻率。这里给出了回线源瞬变电磁测深一维“烟圈”反演解释方法的基本原理、计算方法和应用实例。结果表明,一维“烟圈”反演解释结果纵向分辨率较高,能反映出地下电性变化情况,可以用于瞬变电磁测深资料解释工作。     

回线源瞬变电磁法一维反演算法

由于传统的阻尼最小二乘法只适合于模型较少的简单模型,因此当介质的层数较多时,反演就会受到多解性的影响,有时甚至出现不收敛的情况,并且反演十分耗时。为此,使用正则化思想引入模型约束进行反演,且正则化因子通过计算每次迭代的数据目标函数和模型目标函数自适应得到,使反演能够稳定地进行;引入拟牛顿法来更新雅可比矩阵,大大缩短反演所需要的时间,通过典型的3层与多层理论模型的反演试算,证明了拟牛顿法自适应正则化反演算法对初始模型的要求不高,拟合效果好,收敛速度快,适应性强,体现了良好的稳定性和可靠性。     

地面回线源瞬变电磁法一维反演系统及其应用

瞬变电磁法现阶段的资料处理解释仍以一维反演为主,因此开发一套功能齐全、高效稳定的一维反演系统对进一步提高国内瞬变电磁法的应用水平仍具有重要意义。本研究开发的地面回线源瞬变电磁法一维反演系统包括了基于高斯牛顿法的最小构造反演和Occam反演,也包括了基于阻尼最小二乘法的横向约束反演和空间约束反演。以内蒙古那仁宝力格煤田玄武岩岩体形态探测为例,将该反演系统的最小构造反演和Occam反演结果与商业软件IX1D进行了对比验证,发现不同反演方法获取的电阻率二维断面图中,玄武岩岩体形态相似、电阻率范围一致;结合钻孔资料,这些一维反演结果清晰反映了玄武岩岩体除岩浆上涌通道区域外的分布形态。将横向约束反演和空间约束反演应用于该实例,结果表明:相邻测点间玄武岩岩体电阻率差异缩小,玄武岩与沉积岩界面的连续性得到增强。     

基于在线惯序极限学习机的瞬变电磁非线性反演

基于梯度下降法的传统人工神经网络瞬变电磁反演方法计算效率低,不能保证全局收敛。为了解决上述问题,提出一种在线惯序极限学习机(online sequential extreme learning machine, OSELM)的瞬变电磁反演方法。该方法针对瞬变电磁法所获取的高维勘探数据进行建模反演,首先,通过随机设定隐层参数(输入权值和偏差)来简化模型的学习过程;然后,将测试得到的预测样本加入训练样本中,作为下一次的更新信息,建立在线贯序极限学习机预测模型,从而最大限度提高反演精度;最后,设计了两个经典的瞬变电磁层状地电模型并进行了拟二维地电模型的反演。反演结果表明,该方法能够较好地解决瞬变电磁法高维数据非线性建模的反演问题,同时相较极限学习机(extreme learning machine, ELM),非线性反演方法具有更加准确的反演结果、更好的泛化能力以及更高的计算效率,为神经网络在地球物理反演中的应用提供了新思路。     

井下瞬变电磁法中全期视电阻率的牛顿求解法

全期视电阻率的求解是瞬变电磁法进行地质反演的关键环节。针对井下瞬变电磁法中全期视电阻率的隐函数特点,引进非线性函数求根的牛顿法,给出了井下全期视电阻率的牛顿求解法及其算法步骤。同时还与二分法做了比较,说明牛顿法与二分法一样,可以满足求解全期视电阻率的任意精度要求,并且牛顿法比二分法具有更快的收敛速度。     

基于非结构三角网格的海洋CSEM和MT二维联合反演研究

采用非结构三角形网格结合快速Occam算法对海洋可控源电磁(CSEM)与大地电磁(MT)数据进行二维联合反演试验研究.非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,对反演目标区域采用精细网格剖分,其他区域采用粗网格剖分,在满足精度的前提下减少了不必要的计算量.将CSEM与MT数据加入到同一反演数据集中,通过对联合反演数据权重公式进行推导,构建CSEM及MT数据相关权重因子,控制不同数据的拟合权重来实现联合反演.最后对不同模型进行了反演计算,结果表明联合反演比单一反演对海底构造与异常体的还原度更高,验证了联合反演算法的可靠性.