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

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

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

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

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

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

多种频率域电磁法在冷水坑矿集区的应用效果对比

江西省冷水坑矿集区具有浅部斑岩型矿体与深部层控叠生型矿体的复杂组合.利用现代地球物理方法对矿集区地下深部结构进行了探测,对圈定地下容矿、示矿地质信息和认识地下深部成矿过程具有重要意义.本文以位于江西省贵溪市境内的冷水坑银铅锌矿田为试验研究区,开展了混场源音频大地电磁测深法(EH4)、音频大地电磁测深法(AMT)以及可控源音频大地电磁测深法(CSAMT)三种不同频率域的探测试验,研究和总结了频率域方法以及不同仪器电磁法数据采集特点和应用效果,为频率域电磁法勘探提供借鉴和参考.通过对冷水坑矿集区典型勘探线上三种电磁方法的数据处理,绘制了视电阻率和相位单点测深曲线、原始拟断面及二维反演断面图.发现单点曲线形态基本一样.由于AMT在低频段受矿区工业干扰,视电阻率曲线呈45°上升;视电阻率、相位拟断面对地下电性分层结构基本一样,总体可分为三层,且分层结构较清晰;二维反演断面能够较好反映地下高低阻异常响应,但响应效果有所差异且CSAMT法效果较好.因此在干扰严重地区,建议选用人工源发射的可控源音频大地电磁测深法(CSAMT).     

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

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

二维水平电偶极变频测深阻抗视电阻率的有限元正演计算

在二维地电断面情况下,电偶源变频测深阻抗视电阻率的有限元计算,存在着两个困难:一个是对源的处理;另一个是边界条件的确定。本文提出了解决这两个问题的方法,构成了含源波动方程的边值问题,用有限元法求出了相应的电磁场分量,从而计算出二维地电断面的阻抗视电阻率。对各种地电断面计算的结果表明,计算方法是正确的,可对实际地电断面进行有限元模拟,求得其电磁响应。     

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

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

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

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

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

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

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

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