电性源瞬变电磁B场全程视电阻率计算

在传统电性源瞬变电磁资料解释中,常常采用Bz/T求取视电阻率值,由此带来Bz/T计算全程视电阻率多解性.而Bz随电阻率变化是单值函数,由Bz计算全程视电阻率具有减少多解性的优势.为了研究接地源瞬变电磁B场全程视电阻率,文中首先研究了接地源瞬变电磁法Bz和Bz/T响应值随电阻率的变化特性,然后设计不同的地电模型,采用二分搜索算法对B场全程视电阻率进行计算,结果表明全程视电阻率对电性目标体分辨率高,分层能力强,能较好地反映地电断面的信息;另外,全程视电阻率仅与地层的电性结构有关,与收发距无关,便于野外进行接地源瞬变电磁近源测量.在实际应用中,通过把观测的二次感应电压转换成磁场微分数据,再对磁场微分数据进行积分,求取B场数据,并进而求取视电阻率,取得了较好的应用效果.     

SOTEM响应特性分析与最佳观测区域研究

电性源短偏移距瞬变电磁法(SOTEM)是目前研究和应用较为广泛的一种人工源时间域电磁法工作装置,对深部资源地球物理精细探测具有一定的实际意义.为了深入理解方法内涵并更好地进行推广应用,本文基于电性源瞬变电磁一维正演理论,研究了SOTEM地下感应电流扩散、多分量电磁响应平面分布、多偏移距衰减等特性,然后根据上述特性研究了SOTEM的最佳观测区域.研究结果表明:电性源在地下可以产生水平和垂直两个方向的感应电流.其中,水平感应电流又分为上部水平感应电流和下部水平感应电流(又称作返回电流),水平感应电流的极大值主要集中于发射源附近并垂直向下扩散;垂直感应电流极大值沿与地面呈45°角的方向向下、向外扩散,并且具有较低的振幅和较快的扩散速度.电性源激发的六个方向的电磁场分量都具有一定的探测能力,但是考虑到地面观测的方便性和各分量的传播、分布特点,大多数情况仅利用垂直磁场分量Hz(B/t)和水平电场分量Ex.其中,Hz仅对低阻目标体敏感,且敏感区域位于赤道向区域,并集中在发射源附近;Ex既对低阻体敏感也对高阻体敏感,对低阻体的敏感区域位于赤道向区域,而对高阻体的敏感区域位于轴向区域,并且敏感区域距发射源的距离与目标体埋深和围岩电性有关.     

接地长导线源半航空瞬变电磁三维有限体积正演

基于有限体积法开发了长导线源半航空瞬变电磁三维正演方法.首先将均匀半空间产生的背景场和异常体产生的异常场进行电磁场分离,可获得双旋度方程,后采用有限体积法在控制体积内进行积分,可得到控制方程,对该控制方程采用交错网格离散,时间离散采用后向欧拉离散,离散后形成半航空瞬变电磁三方向离散方程,系数矩阵根据Yee网格尺寸和控制体积尺寸进行计算,根据时频转换可得到接地长导线源半航空瞬变电磁背景场,选用PARDISO直接求解器对方程进行求解,其中迭代时间步长逐渐递增,以降低时间步长对正演效率的制约,任意点电磁场采用对角体积加权平均的方法进行插值.最后,选用均匀半空间模型和层状模型进行精度验证,并采用山峰、山谷、山峰山谷混合三类模型分析了接地长导线源半航空瞬变电磁的响应规律,结果显示开发的程序可用于半航空瞬变电磁勘探,可为半航空瞬变电磁数据解释的发展提供理论指导.     

电性源瞬变电磁发射源形变对观测结果影响分析

电性源瞬变电磁法的基本理论是建立在水平电偶极子的基础之上的,但是实际生产中发射源常为接地长导线,随着短偏移距观测方式的实现,一方面发射源的尺寸不能忽略,偶极子条件不再成立,需要新的方式计算电磁场响应;另一方面,受制于场地条件,发射源可能无法水平、笔直地铺设,这种源形态上的改变必然会对接收点处的响应带来影响,并最终影响到数据解释.本文从电性源一维正演理论出发,研究了发射源发生偏移和弯曲给电磁响应和视电阻率带来的影响,并以某金矿探测为实例做出了验证和说明.研究认为源形态的改变对电磁响应和视电阻率造成的影响随着偏移距的改变而改变,对于电磁响应,偏移距越大影响越小,而对于视电阻率,偏移距越大影响越大.因此在实际生产中应尽量保证源的正规铺设,详细记录实际源的坐标位置与形状,并尽量在偏移距较小的区域内观测.     

瞬变电磁法全空间三维伪谱法模拟

本文根据均匀介质中时间域麦克斯韦方程组,应用伪谱法模拟三维全空间瞬变电磁场的扩散.将电磁场的六分量整理成统一的方程,通过初值条件替代源项将源响应问题转化为初值问题.使用快速傅里叶变换近似空间偏导数,应用切比雪夫多项式近似演化算子来求解方程.通过与全空间均匀介质瞬态线源响应的解析解的比较,验证伪谱法的有效性.计算和分析不同时刻全空间均匀介质中高阻体和低阻体的响应,说明瞬变电磁法对高阻体灵敏度不高,但能较好的分辨低阻体的特点,并反映出时间域伪谱法的优点,即可以连续地计算出电磁场扩散及与低阻体的响应过程.最后讨论伪谱法未来可能应用于矿井瞬变电磁和海底瞬变电磁等领域.     

海底电性源频率域CSEM勘探建模及水深影响分析

为了探索我国海域油气和水合物等高阻目标体CSEM勘探的可行性和方法技术,本文研究了在海水中水平电性源激励下有限水深海洋地电模型的频率域电磁响应,为进一步的1D和3D仿真计算奠定了理论基础.在推导电磁响应公式时,首先给出了各层介质的Lorentz势,然后根据Coulomb势与Lorentz势的关系,得到了各层介质的Coulomb势.各层介质中的电磁场均可以由Lorentz势或者Coulomb势计算得到,但在有限元计算时Coulomb势具有优势.长导线源的电磁场和势函数可以由电偶源的电磁场和势函数沿导线长度积分得到.文中具体给出了海水中水平电偶源和长导线源在海水层的电磁场公式,并根据该公式计算了不同水深环境下海底表面的电磁场分布,分析了海水深度对海底油气储层电磁异常的影响.结果表明,随着水深减小,异常幅度和形态特征发生明显变化.当水深很浅时(如50 m),只有同线方向的E_x和E_z两个电场分量存在明显异常.最后,以两个已知海底油田为例,计算了不同水深环境下可观测到的电场异常,展示了电性源频率域CSEM在海底勘探中(包括浅海环境)的良好应用前景.对于该方法实用化过程中还需进一步解决的问题,文中结尾部分也进行了初步探讨.     

任意形状水平接地导线源瞬变电磁法一维正反演研究

本文发展了一种适应于任意形状水平接地导线源瞬变电磁法的一维正反演方法,可用于长偏移距、短偏移距、多通道测量接地导线源瞬变电磁数据.正演误差通常在1%以内.弧形导线源与直导线源的正演对比表明:导线源赤道方向电磁场容易受到导线源形态的影响,而轴向方向电场所受影响较小.采用修改的矩阵传播法求取水平电场和垂直磁感应场的敏感度.结果表明:相比扰动法,在层数越多时,矩阵传播法优势越明显,当层数达到64层时,效率提高近30倍.反演方法采用正则化Gauss-Newton迭代法实现.反演结果表明:单纯使用磁感应场时间导数进行反演对高导异常层恢复良好,而对低导异常层重构并不理想,且用反演模型生成的电场预测数据与原始电场数据差距较大.使用单一接收点电场进行反演结果也并不理想,而当使用多个接收点处电场值进行共同反演时,高导异常层与低导异常层均能得到较好的体现,相应电场及磁感应场的预测数据与原始数据均吻合的比较好.对有磁导率异常的模型进行反演表明,电场与磁感应场共同反演能够较好地体现地层电导率与磁导率的变化规律.     

矿井瞬变电磁法三维时域有限差分数值模拟

给出了矿井全空间瞬变电磁场的时域有限差分法(FDTD)算法,并推导了Mur吸收边界条件.利用Mur吸收边界条件,选用均匀全空间电偶极源作为初始激发源,模拟了均匀介质中巷道底板岩层内部和层状介质中三维低阻异常体的全空间响应特性,分析了瞬变电磁场在均匀介质中1.9μs和27.6μs两时刻的传播规律及在层状介质中1.2μs和0.023 ms两时刻的传播规律.结果表明:巷道对电场的影响在初期并不明显,在20μs后才表现出来;瞬变场能较好地分辨低阻层,对高阻层的穿透能力强,并且对低阻体反映灵敏;吸收边界条件只有在场域较大、网格节点较多时才能产生明显的效果.所研究成果为矿井瞬变电磁法资料的解释提供了理论依据.     

地-井瞬变电磁多分量响应数值分析（英文）

本文对地-井瞬变电磁法多分量响应进行计算分析。规则局部体瞬变场响应的计算方法与解释模型对于实际导电围岩模型的适用性较差,针对该问题,本文提出了一种基于地下瞬态电磁场数值模拟的计算分析方法。瞬变电磁场模拟方面,本文以时域有限差分法实现正演模拟,引入采用Gaver-Stehfest逆拉氏变换与Prony法的离散镜像法求解初始电磁场,应用透射边界条件保证迭代计算精度。通过均质半空间模型算例,证明该套方法可行。响应分析方面,设定含井旁目标体和导电围岩的地电模型,以上述方法对地下瞬态电场进行正演,以多分量观测装置为例换算感应电动势。通过对比各条件下瞬态电场与多分量响应,得出结论:地-井瞬变电磁多分量感应电动势响应反映了地下瞬态电场沿水平、垂直方向的梯度变化;响应特征取决于地层中瞬变场在不同条件下的"扩散、衰减、畸变"过程和观测位置的电磁场状态。本文的计算分析方法兼顾围岩背景场与局部体异常场,较之传统局部体瞬变场原理能够更全面的反映地质信息。     

ATTEM系统中电流关断期间瞬变电磁场响应求解的研究

在瞬变电磁法中,由于发射电流关断时间不为零、接收线圈的谐振频率有限,早期瞬变电磁信号发生畸变,只能舍弃,因此存在着探测盲区. 针对这一问题,研究了瞬变电磁方法中发射电流关断期间总磁场的形成过程,论证了一次场、二次场和总瞬变场的关系,分析了接收线圈的频率特性和关断时间对瞬变电磁场的影响,提出从总磁场中剔除一次磁场影响的方法,从而获得电流关断期间和电流关断后的早期瞬变电磁场. 采用吉林大学自主研制的瞬变电磁测量系统(ATTEM)在长春市伊通河活断层进行勘探,进一步验证了算法的有效性,缩短了瞬变电磁法的勘探盲区,实现了近地表4 m以下的勘探,可以清晰地分辨近地表的低阻异常,提高了浅层探测精度和分辨率.     

地—井瞬变电磁三维非结构有限元正演模拟研究（英文）

为分析复杂地电模型条件下地—井瞬变电磁响应特征,本文基于非结构网格矢量有限元法实现了时间域地—井瞬变电磁三维正演模拟。为求解双旋度电场扩散方程,我们使用非结构四面体网格对计算空间进行离散,同时选择无条件稳定的后退欧拉方法进行时间离散。数值实验中我们首先以长导线源为例分析电性源电磁场分量的扩散特征和探测能力;然后重点研究Ex、d By/dt零值带的产生机理及分布规律。目标在于指导实际地—井电磁探测中选择合理的测量区域。最后,我们通过模拟复杂模型的异常特征,分析地—井瞬变电磁法探测典型脉状矿体的可行性,为深部矿产资源探测提供参考。     

电性源地-井瞬变电磁法全分量响应特性与探测能力分析

结合地面电性源瞬变电磁法（TEM）探测深度大和地-井TEM探测精度高的优点,本文提出对电性源地-井TEM进行研究.基于一维正演理论,本文对电性源在地下激发的六个电磁场分量的扩散、分布特性和探测能力进行分析研究.研究结果表明地下瞬变场各分量的分布特征与地面情况类似,但是对地层的探测能力与地面情况存在较大差别.各分量的探测能力与它们的扩散特性有关,其中受积累电荷作用,垂直电场E_z和水平磁场Ḣ_x分量对目标层的反映最为明显,E_x和Ḣ_y分量在地下扩散时受"返回电流"的影响会出现变号现象,使得异常体带来的影响被掩盖,不易被分辨.本文研究结果为发展电性源地-井TEM的施工技术、数据处理与解释,建立了相应的理论基础,获得了对该方法的初步认识.     

多孔径瞬变电磁场物理模拟

为了提高瞬变电磁法对地探测精度,提出了一种新颖的瞬变电磁装置一多孔径发射装置,并进行了多孔径瞬变电磁场物理模拟观测.分别使用四孔径发射装置和常规回线发射装置,进行不同收一发距情况下的一次场观测,并对铜板在不同埋藏深度情况下的二次响应进行观测,通过比较两种不同装置下的响应曲线和计算观测数据误差,最终发现使用四孔径发射装置可使一次场强度提高34%,可使二次场强度提高31%.结论表明,多孔径发射源可以改变电磁场的空间分布,认为瞬变电磁场存在类相干现象.研究成果对提高瞬变电磁的探测精度具有很重要的意义.     

海底表面磁源瞬变响应建模及海水影响分析

根据电磁场理论,推导了磁偶源和接收点均位于海水中时层状海底模型的频域电磁场响应一般表达式,并通过此式,得到了海水为均匀半空间和有限海水深度两种情况下,垂直磁偶极装置、中心回线和重叠回线分别置于均匀半空间海底表面时的瞬变电磁响应（磁场和感应电压）表达式. 这些表达式将瞬变响应和海底的电导率等参数有机联系在一起,为海底瞬变电磁法的正演计算和反演解释提供了理论基础. 仿真计算表明,海水的存在不仅使得瞬变响应曲线形态发生变化,而且影响其对海底电导率的分辨能力.     

多通道瞬变电磁m序列全时正演模拟与反演

传统瞬变电磁方法主要用于金属矿勘查,无法满足油气资源高阻目标体的勘探需要.多通道瞬变电磁(MTEM)系统的出现解决了这一问题.该方法采用伪随机序列发射波形和拟地震观测方式,测量同线电场分量,记录全时发射电流及多道观测数据,实现对高阻薄层的高精度探测.鉴于国内对此方法的研究还处于理论探索阶段,尚未进行相应的仿真模拟和数据处理工作,本文针对m序列发射波形多通道瞬变电磁法的全时正演模拟和反演解释进行研究,为国内正在进行的MTEM仪器系统研发及数据解释提供理论指导.我们利用方波响应移位叠加和电流导数与阶跃响应褶积两种方法实现理论m序列和实际发射波形的全时正演模拟;再通过相关辨识技术,削弱噪声影响,计算脉冲响应;最后对积分得到的阶跃响应进行共中心点道集数据联合反演,获取地下电性分布信息.     

积水采空区瞬变电磁响应微分解释方法

瞬变电磁(TEM)视电阻率早、晚期计算公式在求取时进行了简化近似,造成部分地电信息的失真.瞬变电磁响应信号衰减过程与介质的电性特征有关,由此引入瞬变电磁响应微分衰减速度和二阶微分衰减速度两个解释参数,力图避开视电阻率计算,以衰减速度表示瞬变电磁探测结果的电性特征.首先研究了不同导电均匀半空间与低阻目标体地电模型瞬变电磁响应微分衰减速度在时间上的变化特征,在此基础上对瞬变电磁响应二阶微分衰减速度进行了分析,通过绘制TEM二阶微分衰减拟断面图实现了低阻目标体地电模型的电性层划分.研究结果表明:电性分界面与二阶微分衰减的"0"等值线相对应,可以用"0"值等值线作为电性分界面的划分依据,即瞬变电磁微分解释方法能够区分不同地电断面,可以用于瞬变电磁探测成果解释.     

接地源瞬变电磁短偏移深部探测技术

对于接地源时间域瞬变电磁法,当选取适当的激励波形后,可将辐射场与自有场分离开来, 实现频率域电磁法无法实现的近源深部勘探;水平分层大地的解析分析表明,随着偏移距的缩短,接地导线源的场对地层的反映变得更为灵敏;时间域瞬变电磁法的探测深度主要由观测时长决定.基于接地源近场测深的优越性,作者提出短偏移瞬变电磁探测技术并首次命名为SOTEM, 采用了1000 m的偏移距对埋深为1400 m的某盐矿溶腔进行探测, 在全期视电阻率-深度剖面上圈定的溶腔分布被钻孔所揭露, 验证了SOTEM方法的探测能力.该方法为大深度、高分辨探测地下矿产资源提供了新的技术手段.     

瞬变电磁场的直接时域数值分析

为了深入了解瞬变电磁场的勘探原理,直接在时间域对负阶跃脉冲激发的二维瞬态场进行了数值分析.采用的方法是从反映电磁场基本规律的麦克斯韦方程组出发,导出时域电场的齐次扩散方程,对所研究的空间区域作差分离散,源作为初始条件加入,利用准静态近似处理空中边界,然后进行时间的逐步递推,由此展现瞬变电磁场在地下扩散随时间发展的全过程.通过模拟计算不同时刻瞬态电场在地下的分布形态及地面上感生电动势相应的变化,揭示了低阻异常体对感应涡流的聚集作用,低阻覆盖层对瞬变场扩散的减速作用,及瞬变场的延时效应.因此,瞬变电磁法对低阻体是敏感的,有上覆低阻层时探测同样的深度需要较长的时间,而延时效应瞬变场的晚期时段可反映埋藏较浅的异常体.     

瞬变电磁快速成像解释系统研究

为了提高瞬变电磁测深法对地探测的垂向分辨能力,文中研究了磁偶源瞬变电磁法快速成像的方法.详细讨论了电磁波与地震波在介质中的传播规律,分析了磁偶源瞬变电磁场在远区的响应特性.根据磁偶源瞬变电磁场与反射地震波的传播规律的相似性,提出了TEM资料拟地震快速成像解释的思路.理论模型正、反演计算结果表明方法理论是可行的,野外实测资料处理结果进一步说明了该方法是有效的.     

层状介质中地下瞬变电磁场全空间效应

全空间效应是矿井瞬变电磁方法固有的特殊理论问题之一.建立了3种典型的层状地质模型,采用时域有限差分法求解Maxwell方程,并引入电磁波衰减因子EA(dB),来分析全空间效应的影响特征及影响因素.研究结果表明,地下瞬变电磁法工作装置中,分离回线装置对低阻目标层的探测效果优于重叠回线装置;相对于半空间低阻层瞬变响应,全空间低阻层瞬变响应的时间范围更宽,因而需要更长的观测时间;全空间效应影响因素为顶、底板内低阻层相对电阻率差异大小和目标层深度与厚度,当目标层电阻率大于顶板内低阻层电阻率时,所观测响应主要为低阻层响应,此时将会增大资料解释的难度.此外,研究还显示,瞬变电磁法对低阻层的纵向分辨率与低阻层厚度相关,厚度越小分辨率越高,当厚度较大时,其纵向分辨能力较差.研究成果为实际工作参数选择及资料处理与解释提供了参考.