1、Mar.2023PLATEAU EARTHOUAKE RESEARCH2023年3月Vol.35No.第3 5 卷第1 期原地震高大连地磁台极低频观测天然源电磁场数据分析孙国涛，曲骞，于治信（大连地震台，辽宁大连116000)摘要：利用大连极低频电磁观测台的连续观测数据，采用天然源电磁法对观测的电磁场数据进行了处理和分析。结果表明天然源电磁场具有明显的周期性变化特征，表现为白天电磁场谱值比夜间的大，日变化和夏季谱值是冬季谱值的2 倍左右的季节性的特征。对地磁台站观测天然源电磁场数据的非地震影响因素主要为雷电、磁暴以及电磁环境（如高压直流输电）等。关键词：极低频电磁观测；天然源电磁场；时间变化特
2、征；影响因素中图分类号：P315.72*5文献标识码：A文章编号：1 0 0 5-5 8 6 X(2023)01-0034-090引言极低频电磁观测的布置和观测是参照大地电磁测深法（MT）开展工作方式，其天然源观测部分是以天然电磁场为场源，研究地球内部电性结构的一种重要的地球物理手段1 。依据不同频率的电磁波在导电介质中具有不同趋肤深度的原理，在地表测量由高频至低频的地球电磁场时间序列经过相关的数据处理，获得不同频率的大地电磁响应,对其响应分析和反演，得到由浅至深的地下电性结构。极低频电磁观测是一项新技术，其观测量不同于其它已有电磁法，它既可观测极低频发射台发射的人工源电磁信号，也可开展类似大
3、地电磁测深法来观测天然源的电磁场。为得到人工源和天然源电磁场随时间的变化，在大连地震台需要开展极低频电磁场的连续观测。极低频天然电磁场观测不仅能获得电磁场的时间序列（类似于地电场观测，还能获取地下电性结构及其变化2 。大连地磁台（以下简称大连台）是极低频探地工程地震预测分系统在首都圈布置的一个极低频电磁观测站，位于大连瓦房店市九龙街道办事处九龙村大连地磁台院内（1 2 1 0 8 E,3920N）。于2013年9 月1 0 日开始观测。主要开展的研究工作有：提高抗干扰能力的数据处理新技术的应用效果研究，获得极低频频段（30 0 0.1 Hz）天然电磁场谱密度分布和变化；研究非地震事件过程可能引
4、起的天然电磁场变化及特点；研究天然电磁场的背景变化特点；研究电磁场和地下电阻率等物性参数可能出现的异常特点及其成因，探索提高捕捉地震电磁异常的方法等。这不仅对大连台极低频电磁观测资料的分析研究有重要意义，同时对华北地区台网获得高水平的研究结果也具有重要的价值。1数据的获取与处理1.1观测装置观测仪器为德国Metronix生产的ADU-07e电磁观测仪，测量彼此正交的电磁场2 。取正北为x轴，正东为y轴，采用十字型布极方式，外接南、北、东、西四个方向的电极，观测南北向Ex和东西向Ey电场分量，南北、东西和垂直三个方向的感应式磁传感器，观测南北Hx、东西Hy和垂直Hz磁场分量。电极采用Pb-PbC
5、l,不极化电极和纯铅板两套电极，两套电极的埋深相距5 0 8 0 cm，为确保接地可靠稳定，防止冬季冰冻和地下水的影响，底部的铅板电极埋深为4m，且在布设时添加了膨润土做降阻剂降低接地电阻，安装之初使用不极化电极，其南北向接地电阻为6 3,东西向接地电阻为1 1 5 Q。目前使用铅板电极，其南北向接地电阻在40 0 Q左右，东西向接地电阻在300左右，南北和东西电极距都是8 0 m。为了收稿日期：2 0 2 2-0 5-30作者简介：孙国涛(1 9 7 4），男，辽宁大连人，工程师，主要从事台站管理及地磁日常观测、资料分析工作。E-mail:35孙国涛极低频观测天然源电磁场数据分析地磁期第减少
6、风噪、温度等因素的影响，磁探头采用深埋，埋深达1.7 m。1.2数据获取台站观测采用预设定的观测时间表，每天观测天然源数据和人工源数据3天然源观测采用三个采用率进行观测，其采样率分别为40 9 6 Hz（高频）2 5 6 Hz（中频）和1 6Hz（低频）,其中高、中频段每1 0 分钟一次记录。每十分钟40 9 6 Hz记录4s,256Hz记录6 4Hz，16Hz采样除每周一次的CPS时钟同步有间断以外，一直连续记录。每天自动记录的起始和结束时间表详见表1。记录的数据以天为单位形成一个单独的时间序列文件，每1 0 分钟的高中频记录处理出一个功率谱和阻抗edi数据文件，此外利用全天高、中和低频记录
7、处理出功率谱和阻抗edi数据文件4表1天然源数据记录名称采样率开始时间结束时间观测时长循环天然源NS-Cyclic4096 HzXX:XO:00XX:X0:044s循环天然源NS-Cyclic256HzXX:X1:30XX:X2:3464s状态文件NS-Status LogXX:X5:00XX:X5:00一连续天然源NS-Cont16Hz00:00:00InfiniteSplit at end of day1.3数据处理数据处理基本与采用大地电磁测深数据处理类似的方法和流程，首先对各观测分量的时间序列数据进行快速傅氏变换，再在频率域进行仪器响应进行校正，然后计算各分量的自功率谱和分量间的互功率
8、谱，最后计算视电阻率和阻抗相位。计算采用Metronix公司最近开发的资料处理软件Procmt，它有三种数据处理方式：全叠加（s t a c k a l l）、相关度阈值法（Coherency Threshold）和中值法（MedianProcessing）。处理结果以EDI(Electronic Data Interchange）形式存储 。本文只分析电磁场自功率谱的时间变化特征，不对视电阻率和阻抗相位等大地电磁响应进行分析和研究，因此，只使用EDI文件中的自功率谱结果1.3.1全叠加法（stackall)全叠加法就是在谱张量、阻抗张量等叠加时选取其平均值作为最后的估算值。大部分情况下，这种
9、方法作为另外一些方法的参考，因为它可以给出数据最整体的情况，让我们清晰的了解到记录到的所有数据质量，指示出数据好的频段范围。然而有时候全叠加方法在死频段或者高频部分给出令人惊奇的好结果。其原因是谱分布范围虽然非常广、非常混乱，但是其分布很平衡，这就可以使得大小相互抵消。但是这些谱估算方法（全叠加法，中值法，相关性阈值法）都会在估算前剔除一些肯定是错误的飞点，因为当你计算过程中，将相关性剔除这个功能打开时候，谱的全叠加也不是真正意义上的全叠加，而是选取了所有正确的数据进行的全叠加。其中均值的计算公式为：1Nx=(1)标准差的计算公式为：1NV1N-1N-1,)/2)2(2);N-=0=01.3.
10、2相干度阈值法(Coherency Threshold)相干度阈值方法，是谱叠加中最经典的一种方法，它计算Ex和Hy、Ey 和Hx的相干度，设置相干度阈值对数据进行挑选。在数据量比较大的时候，可以选取较大的阈值，例如0.8 或者更大，这样可以得到更高品质的数据，但是阈值也不能选取得太大，否则会因为选取条件太过苛刻，选不到相应的值。但是相干度阈值法也有很大的限制，当测点处近场干扰较强的时候，其相干度值也很大，但是得到的结果也很差。1.3.3中值法(MedianProcessing)中值法是通过计算阻抗张量的中值来估算阻抗张量Zxy和Zyx，即将阻抗张量进行排序然后选取第N/2个值。这是一种稳定估
11、算方法，可以有效剔除飞点值。但是如果叠加次数不够，中值和平均值的计算都是没有意义的，建议至少有2 0地原36高震第3 5 卷次叠加。但是计算过程中将上限和下限的值都选取为1则得不到中值结果，因为处理方法需要该中值对Zxy和Zyx二个值进行估计，他们不会有相同的中值,所有这样设置的条件下得到的edi文件为空，没有数据能够满足这种条件。但是若将视电阻率上下限设置特别严苛时，也会得不到相应的结果。使用三种谱估算法进行视电阻率和相位曲线计算,FFT长度和Parzen半径这两个参数对处理结果影响很大。短的窗口长度会使得目标频率分辨率下降，但相邻频段会变得光滑，对于低能量的谱线，短的窗口长度及大的Parz
12、en半径是有利的；存在一定干扰的时序曲线,增大窗口长度,降低Parzen半径，会减少尖峰对目标频率的影响，从而提高资料的信噪比。然而如果窗口长度过大，将降低时间序列的可能叠加次数，这样统计错误就会增加，而如果时序非常短，短的窗口长度也会导致统计错误。为此，在资料处理过程中需要找出参数的最佳组合以期得到理想的结果。实际数据处理是高中频天然源数据每1 0 分钟循环记录的数据，记录完即传输到台站服务器，进行处理，形成一个edi文件。且在全天（一天的记录从当日8：0 0 次日8：0 0）数据记录完后上传到服务器，进行处理，得到一天的edi文件。最终每日可以得到ct、m e d i a n 和stack
13、_all三种处理方式的EDI结果36 9 个：全天人工源（CompleteDayns，单个文件不能小于5 0 Kb)3个、全天上午人工源（CompleteDay_ca)3个、全天下午人工源（Com-pleteDay_cp)3个以及循环天然源（ns)360个。数据存放在/mtdata/CEAcopy/Site_*/年/月/日的文件夹下的各文件夹中。经数据统计及整理，自2 0 1 3年9 月8 日开始，截止到2 0 1 6 年8 月31 日，大连台极低频观测天然源数据缺失5 6 天，上午人工源数据缺失6 4天，下午人工源数据缺失8 0 天，无数据5 4天。在2015年之前，由于软件问题以及仪器运行
14、不稳定等原因造成，数据缺失多，在运行的47 7 天之中无数据的达到了40 天，而2 0 1 5 年以来运行至2 0 1 6年8 月8 日6 0 5 天之中无数据的仅为6 天，运行率大大提高5 经数据计算分析表明，中值法数据处理方式较其他两种数据处理方式得到的结果更为合理，故在以下数据分析中使用中值法数据处理方式的数据。2时间变化特征分析2.1电磁场谱日变特征通过对每十分钟一次采样获得数据处理结果，我们提取出干扰小、相对稳定的时间段作图，发现电磁场谱具有明显的日变特征。图1 和图2给出了频率9 8 Hz连续8 天的不同时期的十分钟电磁场原始谱图。电场谱单位：【mV/km*sqrt(Hz)】ser
15、Ex磁场谱单位：【nT/sqrt（H z)】serEy1.0E+0serHxExserHy1.0E-0serHz1.0E-1Ey1.0E-21.0E-6Hx1.0E-71.0E-81.0E-6Hy1.0E-71.0E-61.0E-7Hz1.0E-82015-01-012015-01-022015-01-032015-01-042015-01-052015-01-062015-01-072015-01-082015-01-0900:0000:0000:0000:0000:0000:0000:0000:0000:00日期/时间图12015年1 月1 日至8 日谱密度日变曲线（9 8 Hz）期37孙
16、国涛地磁台极低频观测天然源电磁场数据分析第电场谱单位：【mV/km*sqrt(Hz)】磁场谱单位：【nT/sqrt（H z）】净serExExserEy1.0E-1serHx净serHy1.0E+0serHz1.0E-1Ey1.0E-21.0E-51.0E-6Hx1.0E-71.0E-81.0E-6Hyy1.0E-71.0E-81.0E-7Hz 1.0E-82015-11-21 2015-11-222015-11-232015-11-24 2015-11-252015-11-26 2015-11-272015-11-28 2015-11-2900:0000:0000:0000:0000:0000:0000:0000:0000:00日期/时间图22015年1 1 月2 1 日至2 8 日谱密度日变曲线（9 8 Hz）从图中整体上看，电磁场功率谱的日变曲线呈现类似正弦波的曲线形态，白天的电磁场谱强（高值一般出现在世界时7 时左右），晚上的稍弱（低值一般出现在世界时1 8 时左右），这与太阳的活动相关。2.2电磁场谱季节变化特征我们以频率9 8 Hz为例，绘制了2 0 1 5 年冬季（2 0
