定源回线瞬变电磁x分量视电阻率计算方法

为获得定源回线瞬变电磁x分量的视电阻率,将定源回线等效为多个垂直磁偶极子,从而计算得到定源回线x分量瞬变响应结果.通过均匀半空间的瞬变响应值来拟合实测瞬变响应每一时间道的值,将均匀半空间的电阻率值定义为该时间道的视电阻率值.计算结果表明该方法有效,计算精度能够满足要求,可为定源回线瞬变电磁法多个分量综合解释提供必要的基础.     

大定源回线TEM地空系统全域视电阻率定义

介绍了一种快速有效的瞬变电磁地空系统全域视电阻率定义方法。利用磁场垂直分量与水平分量的时间域响应定义地空系统全域视电阻率,实现全空间、全时域视电阻率定义。通过深入分析均匀半空间中地空系统的电磁场表达式,发现地空系统的时间域磁场响应可以用与地面类似的多项式表示,利用视电阻率平移算法,实现全域视电阻率定义。理论模型计算表明,该算法无需迭代,速度快且精度高。最后,对比分析了磁场强度法与感应电动势法定义视电阻率对于异常体的分辨能力。     

中心回线瞬变电磁测深法快速电阻率成像方法及应用

利用导电全空间与均匀半空间中心回线源和磁偶极子在阶跃电流激发下磁场公式和扩散速度的定义,导出了不同条件下瞬变场的扩散速度公式。在此基础之上,引入瞬变电磁测深全区视电阻率定义数值计算方法,给出了电阻率成像的一阶与二阶近似公式,从而建立了一套中心回线瞬变电磁测深快速电阻率成像方法。模型检验结果、实测剖面电阻率成像以及对大地电磁测深视电阻率曲线进行静偏移校正效果均良好,说明该方法是可行的、有效的。      

一种地下瞬变电磁法场值曲线的解释方法

均匀大地表面上大回线源在地下形成的瞬变电磁场解析计算式很复杂,无法用常规方法定义其视电阻率,而场值曲线本身分辨率低,进行分层解释困难。对大回线源在地下形成的瞬变电磁场场值,用由第一层介质充满半无限大地时的场值进行归一化的方法进行解释,经理论计算表明,与地面装置条件下的电阻率曲线进行解释有异曲同工之妙     

斜阶跃波层状介质航空瞬变电磁响应数值计算

目前,航空瞬变电磁方法数据解释主要采用阶跃波均匀半空间模型计算视电导率值,而实际航空电磁系统发射波形的下降沿多为斜阶跃,导致解释时计算的视电导率值存在较大偏差。为此,笔者研究了航空瞬变电磁系统发射电流为斜阶跃时的电磁响应数值计算,将发射波形进行拉氏变换,利用G-S逆拉氏变换与241点汉克尔变换相结合的方法,实现斜阶跃波关断后的层状大地模型航空瞬变电磁响应计算;并对均匀半空间和层状大地模型下,不同关断时间和不同飞行高度对电磁响应的影响进行分析。得出结论:不同关断时间,关断后取样延时2 ms时,均匀半空间电磁响应的平均偏差为27.78%,三层模型的平均偏差为32.16%;当飞行高度从20 m增加到60 m时,均匀半空间和三层模型的感应电动势分别减小了43.6%和83.2%。     

等值反磁通瞬变电磁在隧道超前地质预报中的解释方法研究与应用

瞬变电磁法具有对低阻体反应灵敏的特点,在探测隧道掌子面前方含水或充泥断层、破碎带等低阻不良地质体上具有独特的优势。这里首先简单介绍了等值反磁通瞬变电磁法的原理及特点,然后通过对隧道全空间域高低阻模型一维正演曲线响应特征、视电阻率及视纵向电导曲线的解释特点和效果进行较为全面地分析研究,提出并分析了视电阻率等值线图与视纵向电导断面图联合解释方法。最后通过两个应用案例,总结出等值反磁通瞬变电磁在隧道超前地质预报成果解释中的一些特点和规律。结论表明,等值反磁通瞬变电磁方法适用于隧道超前质预报工作,能够有效提高预报成果的精度及准确度。     

大定源瞬变电磁三分量全域视电阻率定义与三分量联合反演

多数瞬变电磁仪器可以同时采集三个分量的响应,但目前瞬变电磁的解释方法却主要以中心回线垂直分量为主。为了充分利用瞬变电磁三分量的信息,介绍一种大定源瞬变电磁三分量约束反演方法。基本思想是利用瞬变电磁三个分量的场分别求取全区视电阻率,并以此为基础建立初始模型,利用阻尼最小二乘法进行三分量反演。首先,将回线源看成是多个电偶极子产生的场的叠加,利用偶极子叠加的原理求取回线的三分量电磁响应,并推导视电阻率的全区与全期定义(即全域视电阻率定义);然后,利用烟圈理论求取视深度,并以此建立反演的初始地电模型。最后,利用阻尼最小二乘法实现三分量约束反演。利用文中视电阻率定义方法对已知模型数据进行处理,视电阻率曲线与模型的电性参数变化一致;以视电阻率结果为依据建立初始模型,利用三分量约束反演方法对层状模型数据进行反演,结果与模型吻合较好。     

轴向各向异性地层瞬变电磁三分量响应特征

多数地层和裂隙积水区的电导率具有各向异性特性,地层的各向异性对瞬变电磁观测会产生较大的影响,尤其是对水平分量特征更加明显。为研究轴向各向异性地层的瞬变电磁三分量响应特征,基于时域有限差分算法,引入电导率张量构建控制方程,实现电导率轴向各向异性三维瞬变电磁三分量正演。通过与各向同性半空间模型和各向异性半空间模型的一维解析结果对比,验证了算法的准确性;建立各向异性半空间模型、层状模型、含水体模型并计算其回线源瞬变电磁三分量磁场响应。结果表明：水平方向电导率各向异性对三分量磁场响应均产生较大影响,其中x轴各向异性对?B_y/?t分量响应的影响大于?B_x/?t,y轴各向异性对?B_x/?t分量响应的影响大于?B_y/?t,z轴各向异性对三分量磁场响应几乎没有影响;浅部地层的电导率各向异性对三分量磁场响应的影响占主导成分;当采集点在x方向距异常体中心小于y方向时,x轴各向异性对三分量响应影响程度大于y轴各向异性,反之亦然。基于本文研究成果,可为各向异性地层条件下瞬变电磁三分量处理解释提供一些有价值的理论借鉴。     

三分量定源瞬变电磁解释技术及其在金属矿区的实验

目前的瞬变电磁仪器已经为三分量瞬变电磁解释提供了硬件基础,但瞬变电磁的实用解释技术仍局限于对垂直分量的解释或一些简单近似技术。为充分利用三分量数据,进一步提高定源瞬变电磁方法的分辨率,针对定源瞬变电磁装置给出三分量全域视电阻率定义和三分量联合一维反演的方法。在进行全域视电阻率定义时,利用迭代求取反函数值的方法替代直接求取反函数的过程;引入三分量瞬变电磁数据的权重以构建联合三分量数据的目标函数,利用可行域法实现三分量联合一维瞬变电磁反演;最后将三分量解释算法分别应用在模型计算和金属矿的实际勘查资料上,以检验其有效性。模型试验和实际资料处理结果表明,三分量解释方法能够利用更多的瞬变电磁数据,获得更高分辨率的处理结果,且具有较快的计算速度。     

轴向各向异性巷?孔瞬变电磁三分量响应特征研究

巷?孔瞬变电磁法在存在明显电导率各向异性的勘探区会产生较大的解释误差。基于时域有限差分算法,通过引入电导率各向异性张量构建控制方程、将矩形回线源电流密度加入Maxwell方程安培环路定理实现任意电流源的加载,以差分代替微分对控制方程进行离散,实现巷?孔瞬变电磁三维模型正演计算。在与解析解对比验证算法计算精度的基础上,构建全空间模型、层状模型和三维块状模型进行正演并分析电导率各向异性对巷?孔瞬变电磁三分量响应的影响程度与方式。结果表明:垂直轴电导率对巷?孔瞬变电磁三分量响应基本没有影响,水平轴电导率对巷?孔瞬变电磁三分量响应影响较大,其中,?B_x/?t响应主要受y轴电导率影响,?B_y/?t响应主要受x轴电导率影响;通过三分量响应的形态和幅值的相互关系可以辨别各向异性介质所在方位及主轴电导率方向;当异常体所处全空间介质为电导率各向异性时,异常体产生的异常响应会被全空间介质电导率各向异性产生的异常响应所淹没。电导率轴向各向异性特征在巷?孔瞬变电磁法解释过程中不可忽略,研究成果为巷?孔瞬变电磁法各向异性解释提供指导,也为各向异性反演提供参考。      

地一井瞬变电磁三维交错网格有限差分正演及响应特性

为了研究地-井瞬变电磁响应特征,获得不同井位低阻薄板异常体的响应规律:首先采用交错网格有限差分技术离散二次场满足的频率域赫姆霍兹方程;然后结合虚框叠加等效和虚拟界面法将发射源和接收位置扩展到任意层位,解决了全空间背景格林函数的计算问题;之后利用MUMPS求解器求解频率域二次场,再经过余弦变换,获得井中任意位置的瞬变电磁响应;最后采用三维模型对本文算法的正确性进行验证,设计了均匀半空间和嵌入低阻薄板的三维模型,获得垂直磁场的三维分布,分析均匀导电半空间中低阻薄板对地-井瞬变电磁三维响应的影响特点.结果表明:本文三维地-井瞬变电磁计算方案的计算精度与前人基本相同;水平导电薄板的存在主要影响板体附近及穿过薄板的井中瞬变响应,最大的特点是在异常体位置附近的中期时间道响应出现变号现象.本文的研究为定性解释地-井瞬变电磁法异常提供一个技术手段,也为地-井瞬变电磁三维反演奠定基础.     

对瞬变电磁测深几个问题的思考(三)——瞬变电磁场关断效应及全期视电阻率的普适算法

首先给出瞬变电磁（TEM）场卷积计算式,并将关断电流变化率函数看成冲击函数,垂直阶跃脉冲激发的场分量看成输入信号,卷积结果就是任意形状关断电流的TEM响应。把电流激发脉冲看成是许多垂直阶跃函数的叠加,通过叠加这些阶跃函数的响应求得任意电流脉冲激发的瞬变电磁响应。当已知TEM场的实测值和关断电流形状,就可按编程的方法计算全期视电阻率。文中给出了2个斜坡关断电流激发的电偶极源TEM垂直磁场算例,一个为均匀半空间的响应及其全期视电阻率,另一个为三层K型地电断面的响应及其全期视电阻率。      

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

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

基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟和响应特征

为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。     

全空间瞬变电磁法资料解释方法中的平移算法

全空间视电阻率公式及全区求解方法是地下瞬变电磁理论的重要组成部分。为探讨全空间全区视电阻率求解方法,笔者将半空间全区视电阻率的平移算法引入至全空间并对其进行了研究。同时,分别采用理论数据、正演模拟数据及实测数据对该方法在全空间资料处理与解释中的有效性和精确性进行了分析。结果显示,平移算法与常规的全区视电阻率解法具有相同的理论基础,但前者避免了传统全区计算方法中的多解性;此外,采用由截距参数绘制的时间剖面图和视电阻率断面图进行联合解释产生了较好的应用效果。     

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

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

高阻空洞的瞬变电磁响应特征及其探测效果

通过球体电磁响应的计算,对低阻围岩中高阻球体模型电磁响应测量结果进行了分析,获得了高阻球体的异常特征。对干旱地区深埋土洞的工程探测结果表明:高阻土洞的瞬变电磁响应特征与模拟实验结果类似,对高阻目标体具有较好的探测效果;高阻体引起的瞬变电磁场异常范围大于目标体实际规模。      

Effect of inclined conductivity anisotropy on frequency induction and TEM data

The harmonic electromagnetic field of a vertical magnetic dipole above an anisotropic half-space has been simulated using a forward algorithm for layered conductive media with inclined anisotropy. Inclined anisotropy has been found out to change the typical behavior of frequency and transient responses. Qualitative interpretation of FD loop–loop responses of a conducting earth with inclined anisotropy requires taking into account the receiver azimuth dependence of apparent resistivities. In the case of time-domain measurements, this dependence is absent but the apparent resistivities are higher at late times.