引言

地磁场可分成变化非常缓慢的稳定磁场和经常变化的变化磁场.变化磁场携带了地球空间环境的丰富信息，提供了大地电磁测深的源场，在震磁研究中能够反映震源孕震过程电磁性质的变化以及与地震孕育发生有关的内外环境因素的变化，因此变化磁场的研究具有重要的理论意义和实用价值(丁鉴海等, 1994).

许多与地震有关的电磁现象在不同频段范围均已被观测到, 特别是在超低频(ultra-low frequency, 简写为ULF)频段，许多大地震例如1988年Spitak地震(Molchanov et al, 1992)和1989年Loma Prieta地震(Fraser-Smith et al, 1990)前也记录到了岩石圈电磁辐射异常信号. ULF辐射异常引起了国内外研究人员的关注，Hayakawa等(1996)通过对GAM台的数据进行极化分析，发现1993年关岛M8.0地震前该台数据出现了频带范围为0.02—0.05 Hz、幅值约为0.1 nT的辐射异常，并认为极化值在地震辐射识别方面的效果较好. Molchanov和Hayakawa(1995)通过岩石破裂实验对震前电磁辐射异常的物理机制进行了讨论，从微破裂过程中会产生随机电流的角度出发，在对膨胀和黏滑不稳定性进行简单假设后，发现了ULF电磁辐射现象，证实了在岩石破裂过程中会产生ULF电磁异常信号. 杨少峰等(1998)对1996年新疆喀什地区地震活动前观测到的地磁脉动异常进行研究分析，并初步探讨了异常的产生机制，其结果表明喀什地区记录到的特定范围内的磁信号实际上是孕震过程中地球内部破裂运动所产生的电磁扰动. 杜爱民等(2004)对喀什地区的ULF电磁辐射异常进行了更深入的分析研究，发现该地区震前水平磁场分量的偏振主轴方向与地震活动方向近乎垂直，并给出一种电磁辐射激发机理. Hattori(2004)利用极化分析方法发现1997年鹿儿岛(Kagoshima)M6.5和M6.3地震前垂水台地磁数据极化值显著增大，岩手M6.1地震前松川台站的极化值在震前两周出现了3次增强现象；使用主成分分析法对伊豆半岛西部的地磁台阵数据进行分析，其结果表明第三个本征值与地震活动关系密切，并且通过推算得到的ULF磁场信号的来源方向与震群活动方向相对应. Thomas等(2009)对1993年关岛地震前后GAM台地磁台数据进行了重新检验，认为其极化值的长期趋势与全球地磁场活动相关. 冯志生等(2010)利用极化方法从喀什台数据中提取了震前ULF磁场异常信号，结果显示周期为5—100 s的极化值具有年变化规律，采用傅里叶拟合技术消除年变，能够起到去除外空场变化的作用. Li等(2011)利用天津静海台和新疆喀什台的数据计算地震发生前后的极化值，在此基础上提取了较为可靠的地震地磁异常信息；李琪等(2015)利用极化方法对云南地磁台阵数据的分析结果显示，云南宾川M_S5.0地震前也存在地磁异常.

地震前兆现象是复杂多样的，不同的地质构造条件和力学条件等均会导致孕震立体环境的差异，一种前兆异常提取方法的普适性需要更多震例的验证.本文拟将极化方法应用到成都台地磁数据中，对其进行极化分析，从中提取震前ULF电磁辐射异常，以期研究异常的统计有效性，并给出适合成都台的磁场极化值异常判定标准.

1.   台站构造背景及数据

本文以成都地磁台的秒数据为研究对象.成都台位于龙门山构造带东侧前缘的成都新生代凹陷盆地内，距龙门山断裂带约20 km，西距安县—灌县断裂20余km，东距龙泉山断裂约60 km，周围广泛分布第四纪冲积层.该台站附近区域中强地震时有发生，为震前ULF异常的提取、研究提供了充分的条件.

成都台FHDZ-M15地磁组合观测系统于2007年开始观测，在投入观测初期仪器工作正常，经过一年多的运行后，垂直Z分量观测数据的噪声逐步变大，致使观测数据基本无法使用，经台站排除干扰后数据恢复正常. FHDZ-M15系统是一种连续测量并记录地磁场偏角(D)、水平强度(H)、垂直强度(Z)和总强度(F)随时间变化的相对记录仪器，采样间隔为1 s，分辨力为0.1 nT，仪器峰-峰值噪声＜0.1 nT.

Rikitake(1997)对各个频段范围内的观测资料研究结果进行了统计，分析表明震级M越大，可观测到电磁信号的最大距离R_max越大，并给出

按照式(1) 的条件，本文对成都台周边M_L≥5.0地震进行提取，图 1给出了2010—2015年满足式(1) 的8次地震的震中分布图.

图  1  2010—2015年成都台周边满足式(1) 的M_L≥5.0地震震中分布图

F₁：灌县—安县断裂；F₂：映秀—北川断裂；F₃：汶川—茂汶断裂；F₄：龙泉山断裂

Figure  1.  Epicentral distribution of M_L≥5.0 earth-quakes satisfying equation (1) around the station Chengdu from 2010 to 2015

F₁: Guanxian-Anxian fault; F₂: Yingxiu-Beichuan fault; F₃: Wenchuan-Maowen fault; F₄: Longquanshan fault

下载: 全尺寸图片 幻灯片

2.   计算方法

地磁极化值Y_zh定义为某测点地磁垂直分量与总水平分量的谱值之比：

式中，Z(ω)为地磁垂直分量谱值，H(ω)为地磁总水平分量谱值，H_x(ω)为地磁水平分量南北向谱值，H_y(ω)为地磁水平分量东西向谱值.

已有观测和分析表明(Molchanov et al, 1992)当信号源自高空电离层/磁层时，ULF磁场极化值通常小于1，而来自孕震区的ULF磁场极化值Y_zh≥1，因此认为极化值是区别源于空间电离层/磁层的地磁脉动与源于地下岩石层的震磁辐射的关键.

在识别震磁信息时，排除变化磁场本身的背景场及外空场的影响等干扰因素是至关重要的.由于子夜时段的地磁数据包含了最少的空间环境干扰，其次可以避开日间的人为干扰和工业干扰，本研究选取子夜时段的数据进行分析，并对极化值结果进行傅里叶拟合处理，消除背景变化的影响.

3.   结果分析

对成都台2010—2015年秒采样资料进行计算，选取地方时子夜0—2时的数据，每15分钟一个数据段，每天包含8个数据段，计算每段数据的傅里叶谱幅值，并在频率域求平均，然后计算0.005—0.01 Hz频段内垂直分量Z和总水平分量H的谱值积分，剔除低比值后得到每天的极化值，并对结果进行傅里叶拟合来消除长期变化.

图 2a给出了0.005—0.01 Hz频段内的极化值5日滑动平均结果，图 2b为傅里叶拟合结果，图 2c为极化值与傅里叶拟合残差5日滑动的平均值.

图  2  成都台地磁极化值变化

(a)极化值5日滑动平均值；(b)极化值傅里叶拟合; (c)极化值与傅里叶拟合残差5日滑动平均值

Figure  2.  The variations of polarization value of the station Chengdu

(a) Five days moving average of polarization value; (b) The Fourier fitting values of polarization value; (c) Five days moving average of residuals between polarization value and Fourier fitting

下载: 全尺寸图片 幻灯片

变化磁场由外源场和内源场两部分组成，在讨论震磁效应时着重研究内源场及其变化，因此需要抑制外源场的影响，而外源场最明显的特征即为周期变化；本文使用傅里叶拟合方法，滤除长周期变化，从而达到消除外源场变化的目的.由图 2b可见，成都台0.005—0.01 Hz频段的极化值长期曲线具有年变化特征，因此在对数据进行分析的过程中，消除年变是非常必要的，可在抑制外源场影响的同时提高异常识别的准确性.

图 3给出了2010—2015年期间逐年的极化值傅里叶拟合残差5日滑动平均曲线与地震的关系，图中虚线为每年计算结果的3倍均方差线，初次提取异常时将其定为异常阈值线.可以清楚地看出：2010年5月13日曲线出现高值异常，持续时间为5天，异常结束后8天发生汶川M_L5.2地震，震中距为39 km(图 3a)；2011年11月1日青川发生M_L5.6地震，震前几天10月26—29日极化值的傅立叶拟合残差曲线出现高值异常，异常幅值为0.36(图 3b)；2013年4月20日四川芦山发生M_S7.0地震，震中位于龙门山断裂带南段，距离成都台约100 km，利用极化方法分析可知震前18天极化值的傅立叶拟合残差曲线曾出现高值异常变化，异常幅值为0.45，且异常持续4天.为了排除外空场影响，与表示地磁场扰动强度的A_p指数进行对比，由图 3d可看出，2013年4月1日—20日地磁场处于相对平静状态，表明芦山M_S7.0地震前极化值的高值变化与外源场无关，很可能是与地震相关的ULF电磁异常信号相关. 2014年康定M_S6.3地震前也出现了极化值增强的现象，异常持续5天，异常结束后2天发生了地震(图 3e)；2015年成都台附近区域没有满足条件的M_L≥5.0地震发生，极化值曲线表现出相对平稳的特征，两次高值异常持续时间均较短，小于3天(图 3f).

图  3  2010—2015年成都台极化值傅里叶拟合残差(上)和A_p指数(下)的逐年变化

(a) 2010年；(b) 2011年；(c) 2012年；(d) 2013年；(e) 2014年；(f) 2015年

Figure  3.  The Fourier fitting residual of polarization values (upper) and the A_p index (lower) from 2010 to 2015

(a) 2010; (b) 2011; (c) 2012; (d) 2013; (e) 2014; (f) 2015

下载: 全尺寸图片 幻灯片

对2010—2015年的极化值变化进行分析，结果显示几次显著地震前的异常幅值均高于0.3，异常持续时间大于3天.根据这两条标准重新进行异常提取，成都台极化值异常与M_L≥5.0地震的对应情况列于表 1.可以看出：符合判据标准的异常共8次，其中异常发生后2个月内对应地震6次，应报地震8次，漏报2次；异常幅值与震中距呈负相关，相关系数为-0.93，震中距越小，异常幅值越大，反之亦然，如图 4所示.

表  1  成都台磁场极化值高值异常与M_L≥5.0地震对应情况统计

Table  1.  Statistics of corresponding relation between the high value of polarization anomaly at the station Chengdu and M_L≥5.0 earthquakes

异常出现时间	异常持续时间/d	异常幅值	对应地震	发震时间间隔/d	震中距 /km	震源深度 /km	震源机制
2010-01-04	3	0.83	虚报
2010-05-13	5	0.60	汶川ML5.2	12	38.7	20.5	逆冲
2010-07-09	4	0.47	虚报
2011-10-26	4	0.36	青川ML5.6	6	230	21.5	逆冲
			彭州ML5.2	61	51	25.5	逆冲
2012-06-19	3	0.37	虚报
2013-01-02	3	0.39	三台ML5.1	49	141	29.5	走滑
2013-04-02	4	0.45	芦山MS7.0	18	100	17.5	逆冲
漏报			理县MS4.9		104	22.5	走滑
漏报			梓潼ML5.2		150	21.5	走滑
2014-11-16	5	0.33	康定MS6.3	6	211	20.5	走滑

下载: 导出CSV 

| 显示表格

图  4  极化异常幅值与震中距关系

Figure  4.  Scatter diagram for the polarization value anomaly with epicentral distance

下载: 全尺寸图片 幻灯片

4.   讨论与结论

通过提取成都台附近几次中强震前的ULF电磁异常信号可知，成都台附近区域发生的6次地震前均出现了极化值增强现象，异常幅值大于0.3，持续时间为3—5天，在异常结束后2个月内发震，异常幅值与震中距呈负相关.本文给出的震前异常幅值相比李琪等(2015)的研究结果偏低，这是因为本文对极化值进行了傅里叶拟合处理，消除周期变化后再对拟合残差结果进行异常分析的结果.本文结果显示成都台极化值长期曲线具有一定的年变化特征，因此在异常提取的过程中，消除年变是非常必要的，只有这样才能提高异常识别的准确性.

本文中与极化值异常对应的6次地震中，4次为逆冲型，2次为走滑型，逆冲型地震前的异常幅值高于走滑型地震，除与震中距有关外，一种可能的解释是完整岩石破裂所需要的应力水平明显高于岩石黏滑所需的应力水平，由此导致走滑型地震的前兆现象不如破裂型地震丰富(张国民等，2001)，这也可能是导致梓潼M_L5.2和理县M_S4.9地震漏报的原因所在. 6次地震的震源深度为17—30 km，而0.005—0.01 Hz频段内ULF电磁信号穿透深度最大可达70 km，因此本文中给出的震前极化值异常源于孕震区的可能性较大.同时极化值与A_p指数对比结果显示，震前极化值高值异常发生在地磁场平静时段，表明这种异常与外空场关系不大.

从本文得到的结果来看，地震前兆异常与地震的发生并非一一对应，电磁信号的传播途径、地质构造和发震机制的差异均会导致电磁异常现象的不同(Hattori，2004)，有限的震例中得出的认识还不能反映客观存在的复杂情况，本文得到的结果与认识仅仅是初步的，在以后的研究工作中仍然需要继续探索.