S-s方式瞬变电磁测深资料处理方法研究

根据半无限均匀介质空间中电磁场的表达式,本文分析了S-s方式瞬变电磁测深法的全区视电阻率定义的基本性质并讨论了求解方法.在此基础上,引入了基于感应电动势的"虚拟全区视电阻率"概念,从而解决了全区视电阻率的解的唯一性和存在性等问题.理论算例表明,"虚拟全区视电阻率"能较好的反映地下一维电性层电性垂向变化特征,故对瞬变电磁测深资料的定性解释有一定的实际应用价值,同时对一维反演(构制初始模型和确立目标函数)也有重要意义.     

瞬变电磁测深早期数据的修正

通过时间-频率转换关系,TEM数据可以转换成平面波场测深数据,从而可以对TEM资料进行拟平面波场处理解释.在对瞬变电磁视电阻率数据向平面波场测深视电阻率数据转换时,发现由于瞬变电磁使用晚期计算公式及装置问题,使测深曲线早期数据发生畸变.文中建立了视电阻率曲线进入晚期所满足的关系式,从理论上给出不同情况下瞬变电磁测深视电阻率曲线进入晚期的临界点.以瞬变电磁、大地电磁、CSAMT为例,对大量的模型进行正演计算,对计算结果进行对比分析,建立不同地表电性结构、不同时间延迟情况下,瞬变电磁早期数据误差的校正量板.     

频率域线源近区(过渡区)测深理论初步研究

对频率域无限长线源近区(过渡区)测深的理论基础进行了探索性研究。通过计算均匀半空间、层状介质等多个模型的线源响应,进而定义并计算全区视电阻率。结果发现,利用水平电场、垂直磁场分别计算的全区视电阻率均能较好地反映地下电性结构的变化。把水平磁场与水平电场联合通过二分搜索算法计算的全区视电阻率也能较好地反映地下电性结构的变化。在近区通过单分量定义的全区视电阻率与一维MT曲线吻合得很好,因此可以用成熟的MT反演技术来进行线源近区电磁资料的反演,从而把有源的问题转化到无源的问题。这初步表明频率域线源近区(过渡区)测深是可行的,但还有一些问题需要解决。     

长偏移距瞬变电磁测深甚晚期响应及视电阻率的数值计算

文中讨论了提高长偏移距瞬变电磁测深甚晚期响应计算精度及相应的全区视电阻率的计算方法。在频率域 ,利用直接数值积分结合连分式展开提高低频段响应的计算精度 ,从而提高时间域甚晚期响应的计算精度。同时 ,利用连分式渐进展开 ,将均匀半空间电阻率表示为该模型长偏移距瞬变电磁测深响应的反函数 ,利用迭代方法求出适合全时间段的全区视电阻率。数值结果表明 ,文中的方法可有效地提高长偏移距瞬变电磁测深 (甚 )晚期数据的模拟与解释精度     

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

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

层状全空间瞬变电磁响应理论研究

研究全空间效应的影响因素及变化规律,可以有效地丰富地下瞬变电磁法的正、反演理论体系.为此,首先建立层状全空间模型,推导了全空间的瞬变电磁响应公式,采用快速汉克尔变换和余弦变换求出了回线中心处的感应电动势,对比分析了叠加方法得到的全空间响应,并引入全区视电阻率来分析瞬变电磁场扩散规律和全空间效应影响因素,最后通过实测模型来验证模型的计算效果.结果表明,层状全空间模型计算出的响应值相对误差低于3％,精度较高;叠加方法计算的瞬变电磁响应不是真正的全空间,且全空间响应均大于半空间响应的值;在低阻层存在时,全空间响应早期受电阻率和厚度的影响较强烈,由此发现在全空间条件下减少巷道内低电阻率干扰体是减小盲区的重要手段;在超高阻层存在时,半空间响应仅仅是全空间条件下观测点以上的地层电阻率为无穷大的特殊情况;低阻目标层的视电阻率受目标层埋深与层厚差异的影响,若埋深与厚度相差小则分层效果越好;实测数据的建模结果显示计算方法有一定的正确性.本研究填补了矿井瞬变电磁一维全空间响应研究内容,为实际应用提供了验算基础.     

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

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

考虑关断时间的瞬变电磁视电阻率计算及不同波形浅层分辨特征分析

发射电流波形与关断时间会影响瞬变电磁早期响应,导致瞬变电磁浅部信息缺失形成探测盲区．为了分析发射波形对瞬变电磁探测盲区的影响,文中实现了考虑发射波形影响的瞬变电磁全域视电阻率定义方法,并基于该算法对常用的发射波形瞬变响应对浅层电性介质的分辨特征进行分析．首先,采用电流波形与阶跃响应的褶积实现了任意发射电流的全波形瞬变电磁响应计算,将地下电阻率作为全波形响应的函数,利用反函数原理实现了全波形响应的视电阻率计算．然后,以航空瞬变电磁装置为例,利用文中视电阻率计算方法分析了层状介质不同发射波形瞬变电磁视电阻率曲线特征．根据不同盖层厚度视电阻率曲线可知,当覆盖层较厚时,不同发射波形的视电阻率成像结果与模型电性特征基本相近;当盖层厚度较小时,半正弦波、三角波和梯形波等off-time较长的发射波形瞬变电磁响应不能反映浅层的电阻率信息,早期道视电阻率受到深部电性影响,形成了一定范围的探测盲区．最后,利用文中视电阻率计算方法对不同发射波形的复杂三维模型正演的数据进行分析,半正弦波、三角波和梯形波的瞬变电磁响应都能反映异常的地表位置,但会丢失浅部地层信息,深部异常出现时间提前且纵向拉伸．这是由于关断时...     

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

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

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

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

有限元直接迭代算法及其在线源频率域电磁响应计算中的应用

采用有限元直接迭代算法实现了线源频率域测深电磁响应的二维正演计算. 首先给出了线源正演问题的有限元直接迭代格式,然后由迭代法进行求解. 在处理奇异源问题上,采用向内递推的组合网格技巧,在源点附近可进行局部加密,并实现粗细网格的对接,从而较好地解决了奇异源附近的计算问题. 还提出一种迭代求取全区视电阻率的方法,避免了远近区的划分. 通过对均匀半空间、层状介质和二维模型电磁响应的计算,获得了与大地电磁测深相似的视电阻率曲线,验证了算法的正确性;通过对计算结果的分析,在理论上说明了线源频率域近区测深的可行性.     

电偶源频率电磁测深二维阻抗视电阻率计算的源效应校正法

研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上，提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法（大地电磁测深二维ＴＥ极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车）。模型试验表明，利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。     

CSAMT全区电阻率法数值模拟及应用探讨

常规的可控源电磁法理论在计算视电阻率公式上,多半采用其电磁场的渐近特征,难以直接反映全区视电阻率的值,及直观地显现地下介质的地质构造.文中采用水平偶极子激发的电磁场,提出了电场的全区精确表达式,直接计算出大地电阻率.利用汉克尔数值滤波算法和逆样条插值算法对水平层状电磁场进行正演计算,并与计算的卡尼亚视电阻率的对比和野外试验结果表明：该方法的结果在远区等价卡尼亚电阻率,在近区和过渡带则明显地改善了卡尼亚电阻率的非波场区场畸变,从而能更好地接近基底的真电阻率,更形象地反映了地下介质的垂向电性变化.     

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

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

从瞬变电磁测深数据到平面电磁波场数据的等效转换

中心回线源装置的瞬变电磁场是一种涡流场,在地下主要以扩散形式传播.平面电磁波在地下介质中传播时具有反射、折射等电磁学特性,为了借鉴平面电磁波场测深成熟的解释方法,文中开展了从瞬变电磁涡流场测深数据到平面波场数据的快速等效转换的研究.通过大量的理论模型的正演计算、曲线对比,误差统计,对瞬变电磁场、平面波场在地下传播机制特性的分析,建立了一种由瞬变电磁测深视电阻率数据向平面波场测深视电阻率数据转换的时间-频率等效对应关系.利用这一关系式,把扩散场数据快速转换成等效平面波场测深数据. 为进一步进行瞬变电磁拟平面波成像解释打下理论基础.     

瞬变电磁法在煤矿采空区及富水性探测中的应用

瞬变电磁法作为地球物理勘查的重要手段之一,利用对电阻的异常灵敏反应,将采集到的数据分析处理、绘制出视电阻率测深剖面、视电阻率解释顺层切片,可直观地对煤系地层的采空区展布、含水层及其富水性、断层的富、导水性进行评价,可有效避免矿井开采过程中带来的地质危害及资源浪费.     

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

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

线源时间域电测深全期全区视电阻率求解

由于线源一维的正演国内外研究尚少,文中主要利用水平电场分量求取全区视电阻率。通过计算发现水平电场与电阻率之间是一一对应的关系,这使得定义全区视电阻率具有惟一性和可靠性。文中以水平电场定义视电阻率,然后利用二分搜索算法进行全区视电阻率求解。模型计算表明,该方法定义的全区视电阻率能较准确地反映地下电阻率信息。进一步研究发现:对时间域线源而言,在近区算不好晚期的数据,因为电磁波的传播、衰减很快,同样在极远区算不好早期的数据,而且探测深度是随收发距的增加而增大。最后,加入5%的随机噪声后与理论结果进行对比,得到了收发距与最大探测深度的关系。     

EDA各向异性层状介质MT正演模拟

推导了ＥＤＡ型层状各向异性介质大地电磁测深正演计算公式,对多种模型从不同的频率,测量方向计算了地表视电阻率曲线。结果表明,一般的各向异性介质的地表视电阻率值不仅与测量频率有关,而且与测量方向相对电性主轴的交角有关。     

试论大地电磁测深法在地震前兆探索中的进一步应用

从大地电磁测深理论分析和实际资料数字试验出发,我们探讨了有关提高大地电磁测深法监测深部电性变化前兆效能的某些问题.研究结果表明,在所定义的几种大地电磁视电阻率中,尚未加以利用的阻抗实部视电阻率具有某些独特的优点,可能成为监测深部电性变化的一个较好的参数指标;应根据各记录道误差的分配,计算过程误差的传播和累加对不同大地电磁响应函数产生不同的影响,选择相对精度较高的参数作为主要监测参量;而识别深部电性变化的多参数综合判别方法应是基本的分析方法.实际观测到的反映地震前兆的某些大地电磁资料证实了上述论证.