PRC_MTMV和EMTF在长周期大地电磁数据处理中的应用

长周期大地电磁信号强度弱,频带宽,容易受到多种噪声的干扰,观测数据不满足加性、平稳和最小相位特性.因而,其数据处理仍有较多难点.常规的大地电磁数据处理方法是相似的,其主要差别是在时域或频域压制噪声、提高信噪比的数学方程和算法不同.为了提高数据处理的可靠性,拓宽方法的选择范围,本文介绍了广泛使用的PRC_MTMV和较少使用的EMTF两套软件系统的处理原理、参数配置以及用于处理长周期大地电磁资料的基本流程.利用两套系统对实测数据进行处理,对比发现,PRC_MTMV处理得到的视电阻率和相位曲线总是比EMTF的曲线平滑;EMTF的远参考处理对视电阻率曲线高频部分有较大的改善,说明EMTF能够有效克服磁场干扰;将同一测点长周期和宽频测深曲线拼接,二者的效果是相当的,表明EMTF能够满足需要并用于处理长周期大地电磁数据.

一种基于亥姆霍兹分解的大地电磁测深有限元正演预条件解法

本研究针对大地电磁测深法有限元数值模拟中,迭代法求解线性方程组效率较低的问题,利用亥姆霍兹分解原理,将电场矢量双旋度方程的预条件问题转化为基于矢量位的泊松问题和基于标量位的拉普拉斯问题,并在四面体非结构化棱边元离散的情况下,借助节点元辅助网格离散上述预条件问题,进一步利用代数多重网格方法(AMG)实施求解,最终实现预条件算法.利用经典的COMMEMI理论模型进行试算并与前人的积分方程解进行对比,验证了本文数值模拟程序与预条件方法的正确性和可靠性.此外,利用不同自由度规模的实验模型对这一预条件算法的效率进行了测试.结果表明,这一算法可以有效地提升大地电磁测深法棱边有限元数值模拟迭代法的收敛性,计算效率较通用的不完全LU分解预条件算法明显更高;在较大自由度网格(>1000万)数值模拟计算中,其算法效率及内存占用相对直接解法有较大优势,也使小型工作站上利用较大自由度的有限元网格进行大地电磁测深数值模拟计算成为可能.

基于MATLAB的SSMT-2000大地电磁资料批处理技术

资料处理是大地电磁测深工作中的重要环节。在野外测量数据量大,并且要求相互做远参考处理的情况下,由于SSMT-2000大地电磁资料处理软件在批处理上的限制,室内资料处理的工作量巨大。为了提高工作效率,节约成本,通过研究SSMT-2000的参数文件格式、文件操作及内部模块调用方式,以MATLAB脚本语言为基础编写接口,实现了时间序列傅氏变换,单点阻抗张量ROBUST估计,固定测点或固定参考点的多点远参考阻抗张量ROBUST估计以及多个测点相互远参考阻抗张量ROBUST估计的批量处理。采用批量处理,可为室内人员节省大量时间,更能为整个大地电磁测深工作节约资源。

长周期大地电磁测深数据采集及处理技术

研究岩石圈底部的上地幔盖层甚至软流圈的电性结构,需要采集周期长达几万秒的长周期大地电磁场信号,这对仪器设备、采集以及处理技术都提出了更高的要求.本文论述了用于采集长周期信号的仪器设备及其特殊性,讨论了长周期大地电磁测深测站的布置要求,并采用长周期大地电磁测深仪器试采了长周期大地电磁场信号.对试采信号的处理结果表明,采用长周期仪器设备和稳定性好的不极化电极,合理布置测站,能够采集到最长周期达几万秒的长周期大地电磁场信号.

冈底斯成矿带东段的电性结构特征研究

冈底斯成矿带为我国著名的成矿区域,其东段存在一系列规模较大的矿床.为了研究冈底斯成矿带东段的电性结构特征,对覆盖主要矿集区的大地电磁测深数据进行全面的处理分析,通过二维与三维反演的综合对比得到了可靠的电性结构模型.结合其他地质与地球物理资料,对电性结构模型进行分析得到冈底斯成矿带东段的矿床分布规律:矿床主要分布在地壳浅表的电性分界面附近;中地壳高导体可能通过上地壳隐伏的南-北向断裂控制着与地壳伸展作用有关的矿床;南-北向张性构造与东-西向逆冲-推覆构造的交汇部位可能是冈底斯成矿带东段重要的成矿区域;下地壳可能受到软流圈物质上涌的影响发生部分熔融,从而与中地壳的高导体共同影响地壳浅部的成矿作用.

EMTF在宽频大地电磁测深数据处理中的应用

如何有效地压制各种电磁噪声,提高信噪比,保证大地电磁测深数据质量,始终是大地电磁资料采集与处理的核心问题。为寻求更好的处理效果,使用EMTF系统(单点,远参考以及多站同步MMT)对华北多个宽频MT测点数据进行了处理,并编写Matlab接口来方便地改写处理参数。与SSMT-2000对比,结果表明,EMTF系统在多数情况下可以压制噪声,并获得较好的处理效果,提高数据处理的效率。