BZ参数瞬变电磁法的优势及在沿海滩涂区的验证试验

随着我国一系列发展战略的实施,在沿海滩涂区等低阻环境下开展电磁法深部探测工作的需求越发迫切.传统瞬变电磁法(TEM)一般观测感应电动势(dBz/dt)参数,其响应函数复杂,数据解译难度大,且衰减极快,勘探深度较小,不能满足深部探测的需求.而瞬变电磁法的感应磁场(Bz)参数具有响应函数单调、衰减速度慢等优点,深部探测应用前景广阔.本文首先从理论上讨论了TEM的Bz参数相比dBz/dt参数在响应曲线形态、探测深度等方面的优越性,然后在某沿海滩涂区开展了野外验证工作.理论分析计算和野外试验结果表明,TEM的Bz参数相较于dBz/dt参数有诸多优势,基于Bz参数的瞬变电磁法在滩涂区的极低阻条件下仍具有较强的分层能力、较大的探测深度和较强的适用性.     

海底油气藏地质模型的冲激响应

海洋可控源电磁法(mCSEM)的时间域冲激响应特征可以反映海底油气高阻薄层.本文计算了水平电偶极子源均匀大地半空间,海洋均匀双半空间和海洋四层模型的阶跃响应和冲激响应,提出了瞬变冲激时刻的概念.分析了水平电偶源瞬变冲激时刻与介质电导率的指示关系.对于海底油气高阻薄层宜采用多偏移距同时测量方式,由于在低电导率介质中电磁能量传播得要快,在适当的收发距瞬变冲激时刻会提前到达,提出的瞬变冲激时刻道间变化量可以明确指示高阻薄层的存在及埋深.文中还分析了海水深度对瞬变冲激时刻的影响.由于“天波”干扰,瞬变冲激响应受到一定收发距观测的限制.消除 “天波”影响是时间域和频率域mCSEM数据处理的研究热点.     

电性源地-井瞬变电磁全域视电阻率定义

地-井瞬变电磁(TEM)方法是在地表发射,探头沿钻孔(井)测量瞬变响应的一种井中物探方法.由于接收探头沿钻孔深人地下,因此能够获得更加可靠的地下目标体信息,尤其当存在低阻覆盖层、浅部矿化等地质干扰,或者勘查深部规模不大之良导矿体时,地-井TEM方法的优势更加明显.相对于磁性源系统,电性源系统发射功率大,辐射面积广,更适合地形复杂地区之深部找矿.本文的目的是研究电性源地-井TEM的近似解释方法,首先给出了层状介质电偶极子在地下的TEM响应,进而通过电偶极子叠加的方式获得了电性源条件下的TEM响应.利用感应磁场与均匀半空间电阻率的单调关系通过反函数定理进行了全域视电阻率定义.理论模型的视电阻率计算结果显示,不同深度测点的视电阻率曲线首支不同,而尾支基本一致.这说明了全域视电阻率是测点周围有限范围内介质的综合反映,所以不同深度测点的视电阻率曲线首支所反映的范围是不同的.绘制了不同时间道视电阻率剖面曲线,以期显示地下电性分布规律.计算结果表明全域视电阻率定义能够基本可靠地反映地下信息,从而为该技术在矿区深部找矿的应用提供有力支持.     

海洋直流电阻率法电测深装置响应分析

海洋直流电阻率法是海洋电磁勘探的重要方法之一,已经在许多海底资源探测中得到实践.本文推导了供电电极位于海底时四层水平层状介质电位表达式和视电阻率计算公式.根据快速汉克尔滤波法,计算含有贝塞尔函数的积分,编写MATLAB程序计算水平层海洋地电模型电位和视电阻率,探究海水深度和沉积层厚度对于视电阻率曲线的影响.研究表明海水越深,海平面对视电阻率曲线影响越小;沉积层厚度越薄,视电阻率异常响应越明显.建立了各种地电模型,并用二极装置、三极装置、偶极装置进行探测,对不同装置获得的视电阻率曲线进行分析,得出偶极装置对于滨海低阻层和高阻层异常反应能力较强.分析了低阻环境下不同装置在不同供电电流下接收端的电位信号强度,研究发现供电电流强度仅在收发距较小时对接收端电位信号强度影响较大,且供电电流越大电位信号强度越高.本文的结果为深入研究海洋直流电阻率法装置类型应用于海底资源探测奠定了基础.     

基于全区视电阻率的半航空瞬变电磁法对低阻薄层的探测能力探讨

针对使用感应电动势定义全区视电阻率会出现"二值性"的问题,本文采用与电阻率一一对应的垂直磁场分量结合二分搜索算法计算全区视电阻率,这使得全区视电阻率具有唯一性和可靠性.将全区视电阻率的概念引入到半航空瞬变电磁法中,在此基础上,建立以中间层为低阻薄层的三层地电模型,分别绘制了不同低阻薄层埋深与层厚比、围岩与低阻薄层电阻率比的全区视电阻率异常曲线图,并分析其全区视电阻率异常曲线形态及变化规律,通过计算其相对误差,分析并归纳了半航空瞬变电磁法对低阻薄层的探测能力.结果 表明全区视电阻率对地下低阻薄层分辨率较好,能较好反映地电断面信息,这对半航空瞬变电磁法有效探测低阻薄层具有实际的指导意义.     

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

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

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

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

逆Laplace变换新算法及其在时间域电磁响应计算中的应用

时间域电磁响应的正演计算多是由频率域响应经逆Laplace变换而得到.逆Laplace变换的计算精度和效率是时间域电磁响应计算中方法选择的重要指标.论文分析了几种逆Laplace变换的算法机制,并优选出Talbot算法计算了水平电偶源层状模型的时间域电磁响应.逆Laplace变换常用的算法有折线法、数字滤波算法和Gaver-Stehfest算法(简称G-S算法).折线法需要精细地确定分割步长以提高精度,数字滤波算法系数很多,适应频率范围受计算问题所限,而G-S算法受计算机字长和问题对象的影响大.本文在64位计算平台中计算比较了G-S算法、Euler算法和Talbot算法的节点数对于精度的影响,发现Talbot算法受节点数影响小,计算精度高,适应频率范围宽.最后利用21点Talbot算法计算了水平电偶源轴向偶极装置均匀大地模型径向电场的阶跃响应和冲激响应,计算精度及响应时间范围均优于G-S算法.计算了水平电偶源赤道偶极装置均匀大地模型垂直磁场的阶跃响应和冲激响应,冲激响应峰值时刻对于电阻率的变化响应灵敏,与轴向偶极径向电场响应能力相当,但垂直磁场随收发距增大,衰减较快.根据层状模型阶跃响应晚期渐近值计算的视电阻率,水平电偶源轴向偶极径向电场有能力发现大埋深高阻或低阻薄层,收发距应大于中间目标层埋深的5~6倍方可完整探测,类似的,采用水平电偶源赤道偶极装置测量垂直磁场也能达到与之相当的探测能力.计算结果证实了21点Talbot算法适应不同地电模型、不同观测方式的时间域电磁响应计算.     

各向异性层状介质中视电阻率与磁场响应研究

针对任意各向异性地层,利用极向型和环向型标量位函数,导出相应的直流电视电阻率和磁电阻率的磁场响应关系.计算了各向异性地层的直流电视电阻率和磁电阻率响应特征,重点分析了电阻率测深方法与磁电阻率探测方法对地下各向异性介质的探测能力.文中采用状态矩阵的分析方法,首先采用极向型和环向型标量位构造了各向异性层状介质电场与磁场的通解,利用各层界面电场、磁场的连续性及地面激励源的耦合条件,推导了不同层之间电磁场的状态矩阵,建立了空间电场与磁场的递归计算关系.其次,针对递归计算中指数项数值计算的不稳定性,借用状态矩阵的性质,导出了将不稳定指数计算项转化为稳定的指数项的转换关系.针对横向各向同性(TI)介质中极向位与环向位解耦的特点,导出了电磁场的直接积分解.最后,采用解析解验证了算法的正确性,给出了多层各向异性地层模型的视电阻率和磁场响应曲线,分析了直流电法探测裂缝性地层、估计裂缝分布性状的可能性.     

等值反磁通原理瞬变电磁法探测薄层能力

等值反磁通原理瞬变电磁法(OCTEM)的应用日益普遍.为了研究其探测能力,论文基于正演层状模型的瞬变二次磁场时间导数响应,以包含薄层的响应和不包含薄层的背景响应的相对误差作为评价指标,对比分析了不同埋深、不同厚度、不同电阻率薄层以及不同电阻率围岩等情况下的瞬变电磁响应相对异常.结果表明:OCTEM对低阻薄层较对高阻薄层更敏感,薄层与围岩电阻率差异越大,薄层厚度与埋深比越大,越容易被探测到.以响应相对异常值15%作为临界标准探讨了OCTEM方法的探测能力(薄层厚度与可探测埋深的百分比):对于低阻薄层,当与围岩电阻率比值为0.1时,OCTEM探测能力约为1.5%;对于高阻薄层,当与围岩电阻率比值为10时,OCTEM探测能力约为20%.相同情况下常规单磁源瞬变电磁法的探测能力分别为2%和25%,略低于OCTEM对薄层的分辨能力.