长江中下游南京（宁）-芜湖（芜）段深部壳幔电性结构——宽频大地电磁测深研究

为了较全面、客观地认识南京(宁)-芜胡(芜)地区岩石圈深部结构,探讨岩石圈热结构和壳、幔物质状态等重要科学问题,为矿集区成矿作用和成矿规律的研究提供依据,我们完成了6条宽频大地电磁测深剖面。通过分析各剖面电性成像结果,讨论了研究区地壳-上地幔导电性的"拟三维"结构,发现测区普遍存在连续的下地壳-上地幔低阻层,仅在巢湖冲褶体断开。结合区域热流结果,我们认为下地壳-上地幔有较好的导电性,可能存在局部"熔融体"或"含水剪切断裂",其物质状态很可能是热的、软弱的。此外,我们通过物质状态和电性界面推断了上地幔隆起的位置和长江深断裂带的分布范围,认为长江断裂带不但存在,而且由多条北深南浅的断裂组成,较为复杂。     

大地电磁面积性资料和稀疏测线资料的三维反演解释

对理论模型合成的稀疏测线数据和面积性数据以及实测单一测线数据进行了三维反演。通过对比分析合成稀疏测线数据和面积性数据的三维反演结果发现,使用4个张量阻抗元素作为反演数据的稀疏测线数据的三维反演结果可以较好地反映三维目标体,虽然效果可能不如面积性数据的三维反演结果。实测单一测线数据的三维反演结果也较好地揭示测线周围地下基岩界面的三维起伏情况。合成数据和实测资料的三维反演结果表明:可以采用三维反演解释的方法对大地电磁稀疏测线资料进行反演解释。同时,从合成数据的三维反演结果分析可以看出,对大地电磁稀疏测线资料进行三维反演,应尽可能反演所有张量数据;对大地电磁面积性资料进行三维反演,可以只使用两个非对角元素数据,但使用所有张量数据的效果更好。     

基于差分进化的大地电磁反演

针对大地电磁线性迭代反演依赖初始模型,易陷入局部最优解的特点,设计一种基于差异进化的非线性反演方法。利用该方法对大地电磁一维层状介质G、H、HA型地电模型进行反演研究,在无噪声情况下反演结果和模型一致;在加入10%和20%噪声后,反演仍取得良好效果。数值实验结果表明,该反演方法不依赖于初始模型,具有较好的全局优化能力和抗噪声能力,能有效反演大地电磁数据。     

基于经验模态分解法与小波变换的长周期大地电磁信号去噪方法

针对大地电磁信号具有非线性、非平稳和非最小相位的特点,提出了一种基于经验模态分解法结合小波变换的联合信号去噪方式,将时间序列信号通过经验模态分解,利用连续均方误差准则确定原始信号能量转折点,进而再使用小波阈值去噪法对剩余固有模态函数分量进行去噪,最后重构出消噪信号。通过对实测信号处理前后结果的对比,表明了本方法能够有效地应用于信号时域去噪。     

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

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

隐伏地质构造的大地电磁有限单元法正演模拟

对常见的隐伏地质构造体进行大地电磁场正演模拟,并研究和总结了构造体对大地电磁场的响应特征.由麦克斯韦方程组导出大地电磁场的边界条件,应用有限单元法求解变分问题并得出刚度矩阵,根据线性方程组计算出相应的场值,再由改进的视电阻率定义式求出视电阻率.模拟结果显示:改进的视电阻率定义式相对于传统的Cagniard定义式能有效地提高视电阻率曲线的收敛速度和逼近程度;隐伏地质构造体对大地电磁场的响应具有一定的规律性,通过分析两种极化模式下的正演结果,可以有效地识别地质构造体的存在和分布情况.通过实例分析,得出隐伏地质构造正演模拟与实际观测资料的互相验证关系,为实践应用中的反演工作提供了重要的指导意义.     

广义S变换及其在大地电磁测深数据处理中的应用

广义S变换是一种优于短时窗傅立叶变换和小波变换的时频分析方法,利用广义S变换能够准确定位大地电磁资料中存在的噪声,通过定义时频窗对噪声进行滤除,从而明显提高阻抗视电阻率与相位的估算质量.本文基于S变换和大地电磁测深资料处理的基本原理,研究了基于广义S变换的大地电磁测深资料的处理流程和方法.对理论模拟信号及实测大地电磁场时间序列数据的处理,证实了方法的有效性.     

广义S变换时频域滤波在MT数据处理中的应用

短时傅里叶变换是建立在稳态信号基础之上,它仅能提供信号的频域信息,对信号的时间分辨能力差。这影响了它在大地电磁测深数据处理中的应用效果。S变换是一种优于短时傅立叶变换的时频分析方法,能够提供信号时-频域信息。利用S变换对大地电磁测深数据进行时频分析,有助于实现大地电磁测深数据噪声的时频-域滤波,从而提高大地电、磁分量数据的频谱分析精度。从广义S变换理论出发,分析了各类波形噪声的时-频域特征及其对大地电磁测深数据的影响。针对大地电磁测深数据处理特点,利用广义S变换得到时频谱,采用时频比值和门槛值方法,研究适合压制电磁噪声的时频滤波器和滤波方法。对实际大地电磁测深数据的处理结果表明这个方法提高了阻抗张量的估算质量,验证了该方法的有效性。     

大地电磁测深在山东齐河法王南部地区地热普查中的应用

地热资源的勘查开发,对于合理利用自然资源,缓解能源紧张状况,促进经济可持续发展意义重大。在充分收集前人对区域地质、水文地质、地热地质条件研究成果的基础上,利用大地电磁测深(MT)查明齐河法王南部地区地层结构、主要断裂构造的产状、空间展布及断裂带两侧地层分布情况、岩溶裂隙分布、热储与盖层的分布规律、埋深、厚度、岩性等地热地质条件,取得了良好效果,为地热验证孔位的确定提供了可靠的依据。     

大地电磁场三维地形影响的矢量有限元数值模拟

地形对大地电磁场的影响规律研究对于实际大地电磁测深反演解释具有重要的意义。基于矢量有限元法进行了大地电磁三维地形响应数值模拟。首先与二维地形模拟结果比较,验证了算法的正确性;接着,分别对1个山峰地形和1个山谷地形进行了模拟;最后,比较了二维模拟与三维模拟结果。结果表明：三维阻抗张量的Z_xy模式和Z_yx模式都会受到地形影响;用二维地形模拟结果去近似三维地形响应,二维TE模式与真实三维响应结果差别很大,TM模式接近真实三维响应结果;如采用二维反演处理带地形的三维MT数据,TM极化模式反演更为合理。