基于同步大地电磁时间序列依赖关系的噪声处理

本文从信号与系统的角度讨论了同步大地电磁时间序列信号之间的依赖关系,选取高信噪比的时间序列信号作为先验数据,用最小二乘法估算依赖关系;结合参考道的数据,合成本地道含噪声时段的数据,最后用合成数据替代噪声段数据,组成新数据,从而在时域中去除大地电磁噪声.西藏地区高信噪比实测数据的试算结果表明,无论电场还是磁场,信号之间的依赖关系是相对稳定的,只与先验数据的长度有关,与时间无关;虽然不同参考点之间的依赖关系不同,但都可以精确合成本地点数据,与参考点地下电性结构和参考距离无关.仿真实验显示,去噪后的信号与原始信号基本一致.实测数据处理结果表明,该方法可以有效去除强噪声干扰,抑制中高频段的近场源效应,同时保留了微弱的有效信号,保证了处理结果的正确性.最后针对方差比方法无法识别的方波噪声,提出了一种简单的平移方法,成功去除了持续时间大于窗口长度的方波噪声;将该方法与远参考技术结合,可以有效抑制近场源噪声干扰,获得光滑连续并且可信的测深资料.     

ARM嵌入式系统在海底大地电磁信号采集中的应用研究初探

按照现代地电探测仪器的特点和技术要求,首次将ARM嵌入式系统引入地电信号探测领域．建立了基于ARM的海底大地电磁数据采集系统．面向系统开发过程。提出ARM嵌入式系统实现的硬件和软件设计方案。初步尝试实现了地电信号数据采集仪器的低功耗、小型化和智能化。     

海底MT探测仪器的结构可靠性设计

海底大地电磁探测仪器在严酷的海洋环境下工作．为保证海上作业的效率和海底探测的成功率,一方面应注意尽量小仪器的重量和体积,同时也应保证其结构的可靠．在高可靠度前提下求得最省的材料用量是本仪器系统的设计目标．根据可靠性设计理论,通过对仪器的主要结构部件进行力学验算,从抗压、抗拉、抗剪以及抗变形等技术角度确定出各结构部件的材料配置．海洋试验的结果表明,本仪器系统的结构可靠性设计方案基本合理。     

青藏高原地壳高导层的成因及动力学意义探讨——大地电磁探测提供的证据

大地电磁探测结果显示,青藏高原的中下地壳普遍存在大规模的高导层,其电阻率低于10 Ωm,远低于稳定构造区地壳的平均电阻率值.通过对可能形成地壳内大规模低阻异常地质原因的分析认为,青藏高原地壳中的高导层不可能是由金属矿、石墨层或者单纯的含盐水流体引起的,而很可能是由于地壳岩石的部分熔融或者是部分熔融与含盐水流体共同导致的.这些高导层的形成是与板块运动等动力学过程密切相关的.地壳内的高导层可能是印度板块和亚洲板块俯冲的电性痕迹,其成因与板块俯冲过程中由于摩擦生热导致的岩石部分熔融和脱水作用有关,也可能与岩石圈拆沉造成幔源物质上涌有关.沿高原内主要缝合带均存在东西向连续分布的大规模高导体,有可能是青藏高原下地壳物质向东"逃逸"的证据;其中班公-怒江缝合带可能是最重要的物质运移"通道".     

华北煤田瞬变电磁勘探深度研究

为研究瞬变电磁勘探在华北煤田的勘探深度,利用瞬变电磁感应电压晚期段表达式,通过指定仪器最小可分辨电压,推导出瞬变电磁勘探最晚可检测出有效信号时间及最大探测深度表达式,计算了特定情况下的极限探测深度和最晚有效信号检出时间:在发射磁矩为10002×20 A.m2,地层电阻率分别为10和100Ω.m时,瞬变电磁勘探最晚有效信号检出时间为209.6 ms和52.6 ms,最大探测深度为1823 m和2889 m;以实测典型华北煤系地层瞬变电磁曲线为例,其计算机反演深度达1400 m;结果表明:采用合适的仪器和正确的工作装置,瞬变电磁勘探可以满足华北煤田目前开采深度的勘探需要.     

基于站间天然电磁场单位脉冲响应的大地电磁时间序列去噪方法

针对天然大地电磁场信号在人文活动密集地区易受噪声干扰的问题,本文提出利用两个同步测点天然电磁场时间序列之间的单位脉冲响应,合成本地点受干扰时段的数据,从而去除大地电磁噪声.首先,选择高信噪比时段的数据,采用最小二乘法,估算本地点与参考点之间的单位脉冲响应,再根据卷积定律,结合参考磁场合成本地点的磁场和电场.最后用合成数据替换含噪声时段数据,实现时间域去噪.实测高信噪比数据和含噪数据的处理结果表明,该方法可以高精度合成本地点磁场与电场信号,有效去除本地点电场和磁场噪声,包括相关噪声,提高大地电磁数据质量.     

希尔伯特-黄变换在海底大地电磁测深数据处理中的应用

海水运动产生的电磁场,是海底大地电磁探测数据中的主要噪声之一,影响了数据处理与解释的精度.为提高海底MT探测的应用效果,本文在海底MT信号处理中引入HHT时频分析方法.结合海底MT噪声的特点与海底MT信号的特性,对南海实测数据进行分析,通过对比去噪前后信号的Hilbert时-频谱和边际谱,表明利用经验模态分解及其多尺度滤波特性,能够有效压制海水运动产生的电磁噪声.经HHT去噪后的测深曲线在一定程度上有所改善,仍存在"飞点"现象;再结合Robust估计,明显改善海底大地电磁探测数据质量,得到满意的观测结果.     

大地电磁三维数据空间反演并行算法研究

目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.     

藏北高原主要断裂带电性结构特征

对600线的部分测点及2100线的全部测点构成的五道梁-绿草山大地电磁深探测剖面进行了二维非线性共轭梯度反演,得到青藏高原中北部二维电性结构模型.根据该电性结构模型,结合研究区域重、磁及区域地质资料推断了青藏高原中北部主要断裂的位置、产状和切割深度等信息.研究结果表明,青藏高原中北部发育有F1-F16一系列深断裂.其中,F1(苟鲁山克错-囊谦断裂带)和F9(乌兰乌拉湖-玉树断裂带)共同构成金沙江缝合带的北界,是松潘-甘孜-可可西里地块与羌塘-唐古拉地块的分界线;F4、F10-F12共同构成昆仑断裂带,F4(阿尼玛卿断裂带)是南部松潘-甘孜-可可西里地体和北部北昆仑-柴达木地体的分界线;F6、F13-F16为柴北缘断裂带,由南倾的岩石圈深断裂F6和一系列产状相背、北倾的逆冲断裂F13-F16所构成;F7和F8可能反映岩石圈内产状平缓、隐伏的滑脱构造形迹.     

大地电磁测深二维反演方法求解复杂电性结构问题的适应性研究

为探讨二维反演方法在三维电性结构中的适应性问题,本研究中设计了一系列的二维/三维正演模型进行实验计算,分析了三维高阻/低阻异常体对模型响应的畸变作用,并从反演模式选择和数据旋转方向两个方面进行模型二维反演的对比分析,与三维反演的结果进行了比较,最后采用了实测数据进一步进行了二维和三维反演的比较实验.实验和研究结果表明,在剖面选择方面,在剖面方向与垂直主构造方向相差不大的情况下,截取剖面方向,将电性主轴旋转到垂直剖面方向的二维反演结果与垂直主构造方向的反演结果都可以较好地还原正演模型,在大的构造的反映上并无太大差异.在地下为二维或近三维条件时,正演模型的主要结构都可以较好地被二维和三维反演解析出来.二维的反演结果可能甚至会比三维的反演结果的边界更清晰,更精确.然而,对于具有较强的三维结构的模型而言,其二维反演结果与原始模型可能仍然存在较大差异,其中TM+TP或TM模式的二维反演结果相对更接近原始模型,而TE模式的结果往往会有较大误差,需要在解释时特别注意以免得出错误结论.