大地电磁测深理念综合处理软件体系研究.docVIP

　　大地电磁测深理念综合处理软件体系研究 第 1 章 引言 大地电磁测深（Magotelluric, MT）作为一种重要的地球物理勘探方法，自提出以来（Tikhonov, 1950; Cagniard, 1953），已经在地球科学的众多领域取得了显著成果（魏文博等，2002, 2010; T 数据反演、电性结构模型分析及后续地球物理、地质解释工作的基础和保障，能否获得可靠的数据处理及分析结果成为 MT 工作成败的关键。 1.1选题背景 大地电磁测深（MT）是一种被动源的电磁探测技术，通过在地表测量正交电、磁场分量的扰动值，计算得到地下介质的电性结构信息。MT 采集由天然变化的电磁场而引起的地下介质的感应信号，这一信号具有很宽的频率范围（~103Hz 至~105s），因而根据电磁信号趋肤深度的定义，MT 探测深度的范围可以从地表到达上地幔；既可用于浅部矿产、油气、地热等资源的勘探，又可用于地球岩石圈结构探测、地壳与上地幔物质交换模式、地球动力学机制等重大地球科学问题的研究。一般认为不同频率的天然电磁场信号，其激发源是不同的，其中高频信号（大于 1Hz）主要由某些气象条件产生，如闪电等；而低频信号（小于 1Hz）的场源机制则较为复杂，通常认为来自于太阳风与地球磁层及电离层之间的相互作用而产生的电磁扰动。MT 之所以得到广泛应用，主要由于具有以下优势：1）信号频带宽，因而探测深度范围广，可作为探测地球深部电性结构的有效手段；2）不受高阻层的屏蔽作用影响；3）对高导异常体的反映灵敏；4）野外施工较为简便。近年来，随着电子、计算机、数字处理等技术的飞速发展，大地电磁测深（MT）在资料处理、正反演计算等方面都取得了长足的进步，特别是三维大地电磁测深反演方法的研究（Lin et al., 2011; Siripunvaraporn, 2012），进一步拓展了大地电磁测深实际应用和研究工作的发展前景。目前，在一般情况下大地电磁测深资料的处理技术已相对成熟；但对于复杂环境条件，往往存在由人文干扰所引起的电磁噪声影响，这大大削弱了大地电磁测深的应用效果。因此，人们需要进一步深入研究大地电磁测深数据处理的方法，以期获得高质量的张量阻抗、视电阻率和阻抗相位等参数的频率响应。此外，由于三维正、反演算法的复杂性和海量数据运算的制约，大地电磁测深三维反演尚未得到推广应用；因而，现今对大地电磁测深的反演仍主要采用二维反演方法，这就要求对地下电性结构的维性作较准确的分析，以保障所提取的区域电性主轴方向的可靠性。因此，也仍然必需对大地电磁测深数据分析方法作更深入的研究。显然，MT 数据处理与分析，对于提高 MT 探测精度具有非常重要的意义，是保障 MT 探测效果的关键环节。 自 2008 年以来，我国启动了深部探测技术与实验研究（SinoProbe）专项，专项定位为即将开展的《地壳探测工程》的培育性计划，包括《大陆电磁参数标准网实验研究（SinoProbe-01）》等 9 个项目，其中 SinoProbe-01 旨在通过研究大区域岩石圈尺度的综合物性成像方法，构建中国大陆岩石圈电性、密度等物性结构的三维标准模型。根据项目研究目标要求，需要进行大地电磁测深标准网观测的关键技术研究，并开展华北和青藏实验区阵列式大地电磁测深标准网构建的示范性研究工作。因此，完成Sinoprobe-01 项目的研究任务，需要在各种环境条件下采集大量的大地电磁测深数据，这使我们面临海量、质量参差不齐的 MT 数据精细处理与分析的难题；尤其当《地壳探测工程》启动后，我们将面临的 MT 数据处理与分析的工作压力更是现在所不可比拟的。为此，迫切需要有一套集成化程度高、功能强大、性能稳定、运算结果准确、数据处理和分析效率高的大地电磁测深数据处理和分析计算机软件系统。 第 2 章 大地电磁测深数据处理方法研究 大地电测测深（MT）野外采集得到磁场及电场分量的时间序列数据，这一原始数据最为直观的特点是其数据文件一般较大。简单地讲，MT 数据处理工作便在于从时间序列数据中提取信息，并最终得到依赖于频率的传递函数（Transfer function, TF），并以 TF或者其衍生参数（视电阻率、相位、倾子等数据）作为反演模型的输入参数。将时域时间序列转化为频率域信号，通常是 MT 数据处理的第一步工作，而这一过程的常规工具便是傅立叶变换。经傅氏变换后，对于某一频率，通常有多个时窗的傅立叶系数（Fouriercoefficients, FC）与之对应，这就要求需要对不同窗口的数据进行叠加处理；需要强调的是，实际上同一时窗内，相邻频点亦涉及到基于中心频点的数据叠加问题。在进行 MT 数据处理方法的详细阐述、研究之前，有必要对传递函数（TF）作进一步阐述。假设地球为一线性系统，以随时间变化的磁场分量作为输入信