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

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

大定源装置下瞬变电磁法视电阻率定义

对于应用较广的大回线源瞬变电磁装置,虽然得到了非中心点电磁响应式,但由于响应的复杂性,一直没有推导出不同场点视电阻率定义式,通常借用中心回线方式的视电阻率公式进行解释,必然会在观测区的边缘处出现异常.为提高瞬变电磁法的解释精度,将非中心点垂直分量转化为人们熟悉的电阻率参数,本文提出并建立了大定源装置情况下视电阻率定义....     

地-空瞬变电磁法电阻率成像研究与应用

地-空瞬变电磁法（Semi-airborne transient electromagnetic,SATEM）凭借适应能力强、探测深度大、实时性强等特点,适用于湖泊、沼泽、山区等地形复杂地区观测.SATEM接收数据量大、精度要求高,传统成像方法基于地面视电阻率计算未考虑飞行高度,开展SATEM的高精度快速成像研究对实际运用有重要作用.针对SATEM中电性源发射装置感应电流分布,接收装置设置于空中以及两者不在相同平面的特点,本文考虑横向、纵向感应电流差异,飞行高度对视电阻率影响等因素,提出电性源SATEM电阻率成像方法.首先,基于电性源瞬变电磁场感应电流分布特征及各分量空间分布、扩散特性分析,推导出基于均匀半空间磁场与感应电压解析解,然后定义了电性源地-空瞬变电磁法早期、晚期、全期视电阻率;最后根据对感应电流分布分析定义了横向分量和纵向分量的成像深度,并研究了成像深度受飞行高度、发射磁矩、偏移距、偏移角度的影响;对煤炭采空区的实测结果表明,利用本文提出的电阻率成像方法可以更加精确定位目标体位置.

     

大回线源瞬变电磁全域视电阻率定义

为了进一步满足精细探测的需要,给出了一种大回线源瞬变电磁全域视电阻率数值计算方法,利用该方法可以快速求取大回线源非中心点处多分量的全域视电阻率.利用贝塞尔函数展开法计算大回线源非中心点垂直分量与水平分量的电磁响应,通过分析非中心点处的多分量电磁响应特征,得到基于磁感应强度的大回线源全域(时间上不分早晚,距离上不分远近)视电阻率计算方法.与现有全区视电阻率方法不同在于:对非中心点处多分量视电阻率进行定义,利用磁感应强度定义不存在假极值现象,无需迭代,计算速度快,精度高.     

矿井瞬变电磁法全空间视电阻率解释方法研究

从全空间中心回线源响应公式出发,推导了全空间晚期和全区视电阻率公式;通过与半空间公式比较和对实测资料的分析,研究了全空间视电阻率公式的应用效果.结果显示,全空间视电阻率为半空间视电阻率的(5/2)2/3倍,其晚期视电阻率与真值的相对误差小于0.56%,而相应的半空间视电阻率与真值误差均大于42%.与常规解释方法相比,采用全空间全区视电阻率换算方法能更好地圈定低阻异常范围,且电阻率差异大.因此,该公式更适用于矿井瞬变电磁法资料处理与解释,并具有更好的效果.     

航空瞬变电磁法的全时域视电阻率计算方法

鉴于ATEM数据量大质低的特点,视电阻率仍是最有效而直观的参数.本文给出了适用于ATEM多种观测条件下(包括各种装置类型,发送波形,关断电流等情况)的全时域视电阻率的计算方法,并对典型层状大地模型的理论瞬变电磁响应作了全时域视电阻率及其对应等效深度的计算.结果表明,全时域视电阻率随等效深度变化的曲线可定性反映地下介质电性特征,并可近似估计低阻或高阻层上界面深度.对于非单调下降的瞬变电磁曲线,该方法也能取得较好的计算结果.     

小回线源TEM法全域视电阻率成像在某金矿水害中的应用

在煤矿水害问题中,利用瞬变电磁法研究煤矿的底板水、顶板水、断层或陷落柱特征,已经有很多成功的案例;但是在金矿水害中的应用,尚未发现成功的案例.目前,在生产实践中,瞬变电磁法在煤矿水害资料解释还停留在利用晚期视电阻率近似处理、解释的水平.金矿水害数据采集位于深部采场,需要全空间瞬变电磁资料解释技术.基于以上背景,介绍利用反函数原理计算的全空间全域视电阻率成像,对金矿水害探测的瞬变电磁资料进行处理解释.研究发现:该解释方法不仅计算速度快,且精度满足要求,方便实测数据的快速处理.同时能够去除晚期视电阻率成像在早期、过渡期出现的假异常;通过钻孔验证和开挖,证明该成像解释方法对金矿水害防治,效果非常明显.     

多辐射场源地空瞬变电磁法快速成像方法研究

多辐射场源地空瞬变电磁法目前仍处于研究阶段,迫切需要完善该法的理论和解释体系,以便在生产中推广使用.本文利用地空瞬变电磁法数据量大的特点,结合微分电导与相关叠加合成算法,实现了多辐射场源地空瞬变电磁法快速成像解释.首先基于等效导电平面原理,从电偶极子源出发推导了较高精度的瞬变场快速近似正演计算公式,采用剖分叠加的思想,将其推广到多辐射场源地空瞬变电磁法空间任意一点瞬变场的快速近似正演计算,并用四种典型的层状模型验证了该法的计算精度.以此为基础,进一步推导了多辐射场源地空瞬变电磁法的视纵向电导和视深度计算公式.利用微分电导对电性界面敏感和相关叠加合成算法突出异常、压制干扰的特点,实现了多辐射场源地空瞬变电磁相关叠加快速成像.三维模型试验表明,本文提出的快速成像算法计算速度快,能有效压制噪声,且对地下界面有较好的分辨能力,有助于提高多辐射场源地空瞬变电磁法的解释精度.     

时间域瞬变电磁法全区视电阻率的平移算法

深入分析了均匀半空间中瞬变电磁场的解析表达式,发现瞬变响应曲线随地下电导率、发射回线边长与观测时间具有平移伸缩特性,并据此提出了一种新的直接计算全区视电阻率的方法.该法与现有计算全区视电阻率的方法的不同点在于：无需迭代和求解非线性方程,速度快,精度高,算法稳定,且只有惟一解.     

时间域瞬变电磁法中心方式全程视电阻率的数值计算

给出了一种时间域瞬变电磁法视电阻率的数值计算方法，利用该方法可以容易地求出中心方式的全程视电阻率. 根据中心方式磁场垂直分量时间变化率Bz/T的核函数Y′(Z)的表现特征，以参数Z把整个瞬变过程分为早期阶段（Z>1.6）、早期到晚期的转折点（Z=1.6）和晚期阶段（Z<1.6）. 首先分别得到早期视电阻率和晚期视电阻率的精确值，然后通过转折点构成一条完整的全程视电阻率曲线. 虽然磁场垂直分量Bz的核函数Y(Z)是参数Z的单值函数，但同样存在一个从早期到晚期的转折点Z=1.6，转换点两边仍然可以得到一条早期曲线和一条晚期曲线. 在数值计算中，当迭代步长ΔZ<0.005Z时，视电阻率的相对误差小于0.5髎. 理论模型和实际数据计算表明，与早期和晚期近似值比较，全程视电阻率具有更高的精度和分辨率.     

电性源瞬变电磁地空逆合成孔径成像

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.     

基于神经网络的视电阻率快速算法

本文从瞬变电磁均匀半空间二次磁场响应公式出发,提出了一种基于神经网络的视电阻率快速计算方法.以中心回线为例,根据瞬变响应公式的特点,简化网络结构,选用三层BP神经网络和误差训练算法,用均匀半空间样本数据进行训练,确定了收敛快、误差小的一步正割法和隐含单元数,得到基于不同采样时窗的一组网络参数.用本文方法与二分法、牛顿迭代法做模型计算比较,及最后的实验计算,说明算法的快速,准确.本文方法不依赖初始模型,避开了复杂的电磁场数值计算,实现了视电阻率的快速计算,对瞬变电磁法资料的快速解释有一定的参考价值.     

2.5维起伏地表条件下频率域地空电磁正演模拟

频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特...     

瞬变电磁横磁场响应特征与探测能力分析

瞬变电磁法是重建地下电阻率等电性结构的重要方法.传统磁性源和电性源瞬变电磁法主要观测横电极化场,横电极化场仅对良导目标敏感,对高阻目标的分辨能力有限.横磁极化场对高阻目标具有较强的分辨能力,但未得到有效利用.双线源瞬变电磁法可以增强观测电场中横磁极化场的占比,但对该方法的响应特征和分辨能力缺乏系统性的研究.为此,本文以双线源为例开展瞬变电磁横磁场响应特征与分辨能力分析.双线源瞬变电磁水平电场的响应强度要小于传统接地导线源,在发射源的中垂线上,层状大地模型的响应为零,观测的水平电场响应只能由地下的三维目标体产生.分别提取双线源和传统接地导线源激发电磁场中的横磁场和横电场,双线源瞬变电磁场中的横磁场占比要大于传统接地导线源,特别是在中晚期,横磁场远大于横电场,横磁场占比得到明显增强.通过均方根差和三维数值模拟的计算,双线源瞬变电磁水平电场显示出相较于传统接地导线源更强的高阻目标分辨能力,特别是在赤道向,水平电场对高阻目标分辨能力的增强效果更加明显.

     

Tensor time domain electromagnetic resistivity measurements at Ngatamariki geothermal field, New Zealand

Experimental measurements in the Ngatamariki geothermal field, North Island, New Zealand were made to test the applicability of the time domain electromagnetic method for detailed investigation of the resistivity structure within a geothermal field. Low-frequency square wave signals were transmitted through three grounded bipole current sources sited about 8 km from the measurement lines. Despite high levels of electrical noise, transient electric field vectors could be determined reliably for times between 0.02 and 3.3 s after each step in the source current. Instantaneous apparent resistivity tensors were then calculated. Apparent resistivity pseudosections along the two measurement lines show smooth variations of resistivity from site to site. Over most of the field the images consistently show a three-layer resistivity structure with a conductive middle layer (3–10 Ωm) representing the conductive upper part of the thermal reservoir. A deep-seated region of low resistivity in the northwest of the field may indicate a conductive structure at about 1 km associated with a deeper diorite intrusion. Measurements sited closer than about 100 m to drillholes appear to have been disturbed by metallic casing in the holes. A change in resistivity structure in the east of the field may indicate a major geological or hydrothermal boundary.     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.