前言

地震是应力持续加载并最终导致断层失稳的结果，由于应力应变的持续积累过程伴随着压电、磁效应和动电效应等（Huang，2002；郝锦绮等，2003；Ren et al，2012；张建国，2017），岩土介质的物性参数会发生变化，孕震晚期会产生电磁辐射并引起电磁场变化（解滔等，2018）。一种有效的提取方法不仅能够顺利提取显著异常信号，而且能够深入认识震磁异常成因机理（姚休义，冯志生，2018）。

已有研究表明岩石破裂时会出现电磁波超低频（ultra-low frequency，缩写为ULF）信号异常，且临近破裂时岩石附近磁场强度会出现短周期前兆变化（郝锦绮等，2003）。地磁垂直强度极化方法是目前地震磁扰动定量分析中物理含义明确、分析过程完善且获取异常信息能力较强的方法之一（姚休义，冯志生，2018），该方法基于频谱分析，利用地磁场垂直分量Z和水平分量（H或G）的频谱振幅比计算得到。数值模拟结果（Molchanov，Hayakawa，1995）证明，对于一次源自地壳内频率约为1 Hz的磁信号，当其垂直分量幅值大于水平分量幅值时，则两者比值大于1，故可通过该比值突出岩石圈信号异常，同时抑制外源场的电磁信号。

自从在1988年斯皮塔克地震和1989年洛马普列塔地震前监测出超低频段电磁异常信号（Fraser-Smith et al，1990；Bernardi et al，1991；Molchanov et al，1992；Kopytenko et al，1993）以来，震前超低频异常信号引起了国内外众多科学家的关注。Hayakawa等（1996）利用距离震中65 km的日本关岛地磁台磁通门磁力仪的观测资料，对1993年8月8日关岛M_S8.0地震前后该地区地磁垂直强度极化值的变化特征进行分析，结果显示震前两个月0.01—0.05 Hz频段的极化值逐渐增大，发震时达到最大，震后又逐渐恢复到日常水平。Molchanov和Hayakawa （1995）在探讨震前电磁辐射异常信号的物理机理过程中，通过岩石破裂实验来研究微破裂过程中的微电流等问题，监测到超低频电磁辐射信号，证明了岩石在破裂过程中会伴随超低频电磁异常信号的产生。此后对众多震例的研究表明地震前三个月至震前数天在10—2 Hz电磁信号均呈现出高值极化异常特征（Hayakawa et al，2000；Molchanov et al，2003；Hattori，2004；Prattes et al，2008）。Hattori等（2013）分析了伊豆半岛和博索半岛近十年来地磁观测资料中0.01 Hz左右的超低频地磁信号，认为M6地震前存在地磁垂直分量能量增强的现象，且伴随电磁波超低频异常，该震例统计结果显示，地震基本发生在异常出现后一个月内。国内利用地磁垂直强度极化方法分析震前异常也已有较多成果，例如2004—2007年喀什一系列中强地震（冯志生等，2010）、2009年宾川M5.0地震（李琪等，2015）、2017年九寨沟M_S7.0和精河M_S6.6 （贺曼秋等，2019）、2017年阿拉善M_S5.0 （廖晓峰等，2019）等地震前均提取到了显著的电磁异常信号。

上述研究仅分析地磁垂直强度极化异常与地震的关系，但针对地磁垂直强度极化异常的时空演化过程研究不多。随着近几年我国地震地磁观测台网的大力建设，地磁秒采样观测仪器的密度有显著提高，积累了丰富的观测资料。本文拟利用目前机理明确、应用效果佳的地磁垂直强度极化方法提取震前磁异常信号，并探讨异常时空变化特征与地震震中的关系，以期能够更深入认识震前电磁异常信号的演化特征及展布规律。

1.   数据选取与计算方法

1.1   数据选取

自我国 “十五”以来，中国地震局在全国范围内建设了大量磁通门秒采样观测台站，收集到了丰富的电磁观测资料。本文通过对观测资料进行筛选、逐一对比并进行预处理分析，筛选处理基本原则如下：① 记录长期稳定，且断记、突跳情况较少；② 同台多套仪器同步观测时，选取其中同步性较好的一套资料。按照此原则最终筛选出2015—2020年中国大陆65套仪器的秒值观测资料进行垂直强度极化异常分析，并基于异常判据从全国层面上探讨异常的时空演化特征。

本文的65套观测数据来自磁通门GM4系列磁力仪及其改造仪器（59套）、FHD-M15地磁总场与分量组合观测系统（6套），均为秒值，观测频段为超低频段，台站分布如图1所示。

图  1  全国地磁计算台站分布图

Figure  1.  Distribution of national geomagnetic stations in China

下载: 全尺寸图片 幻灯片

1.2   处理方法

地磁垂直强度极化值为

式中，$Z{\text{（}}\omega {\text{）}}$为地磁场垂直分量的傅里叶谱值，$H{\text{（}}\omega {\text{）}}$为地磁场水平分量全矢量的傅里叶谱值，${H}_{x}{\text{（}}\omega {\text{）}}$为地磁场水平分量南北向傅里叶谱值，${H}_{y}{\text{（}}\omega {\text{）}} $为地磁场水平分量东西向傅里叶谱值，$ \mathrm{\omega } $为圆频率（何畅，冯志生，2017）。

以青海都兰台为例，具体计算步骤如下：① 将磁通门三分量Z，H_x和H_y每天的秒值观测资料划分为96段，每段15 min，分段计算垂直分量和水平分量的傅里叶谱值；② 计算每天5—100 s周期内各频点的极化值均值，获得极化值的逐日变化序列（图2a）；③ 对极化值逐日变化序列进行傅里叶拟合（数据不短于半年），获得周期大于半年的傅里叶拟合变化曲线，此即为极化值的年变化曲线（图2b）；④ 将5—100 s周期内各频点极化值减去年变化后，剔除低于二倍均方差的部分，仅保留高值部分（图2c）；⑤ 对筛选出的极化高值的日均值逐日变化序列再次进行傅里叶拟合，如果仍然存在年变化成分，则再次扣除，所得残差为最终处理得到的极化高值序列（图2d）（冯丽丽等，2021）。

图  2  地磁垂直强度极化方法计算过程（以青海都兰台为例）

（a） 极化值；（b） 极化值的傅里叶拟合；（c） 极化值与其傅里叶拟合的残差值；（d） 残差的5日滑动平均值

Figure  2.  The calculation process of geomagnetic vertical intensity polarization method （Data from Dulan observation station）

（a） Original polarization values；（b） Fourier fitting results of Fig. （a）；（c） The difference between Fig. （a） and Fig. （b）；（d） Five-day moving average of Fig. （c）

下载: 全尺寸图片 幻灯片

2.   结果分析

据中国地震台网测定，北京时间2017年11月18日6时34分19秒，西藏林芝市米林县发生M_S6.9地震，震中位置为（29.75°N，95.02°E），震源深度约10 km。通过整理2017年地磁垂直强度极化计算结果，发现西藏米林M_S6.9地震前全国地磁资料呈现出大范围同步极化异常现象，从2017年10月30日起，异常台站的时序计算曲线呈现单峰或双峰的高值形态，异常持续3天，在平静两天后随即又出现高值异常情况，此次异常持续4天（图3）。通过异常时间演化进程（图4）可以看出：全国地磁垂直强度极化超阈值异常首先出现在新疆天山以北、西藏东南和川滇交界以及黑龙江与内蒙交界北部地区，极化异常面积于2017年10月31日达到最大值，之后开始收缩并集中在青海大武台附近；随后异常范围再一次逐步增大，于11月6日达到异常面积次大值，主要集中于甘青川交界、新疆天山以北以及内蒙古东部和黑龙江地区，本次异常于11月7日结束。此外，本次地磁垂直强度极化高值异常区内还发生了2017年11月23日重庆武隆M_S5.0地震和12月15日青海泽库M_S4.9地震。

图  3  2017年10—12月部分地磁垂直强度极化异常台站时序曲线

Figure  3.  Time series of some stations with geomagnetic vertical strength polarization anomaly from October 1，2017 to December 31，2017

下载: 全尺寸图片 幻灯片

图  4  2017年10月28日至11月8日全国地磁垂直强度极化异常空间演化过程

Figure  4.  Spatial evolution process of geomagnetic vertical intensity polarization anomaly from October 28 to November 8，2017

下载: 全尺寸图片 幻灯片

图  4  2017年10月28日至11月8日全国地磁垂直强度极化异常空间演化过程

Figure  4.  Spatial evolution process of geomagnetic vertical intensity polarization anomaly from October 28 to November 8，2017

下载: 全尺寸图片 幻灯片

总之，此次异常的分布范围呈增大—收缩—消失—增大—消失的演化过程，高值异常多次出现在中国大陆西部，尤其青藏高原巴颜喀拉地块与羌塘地块交界处。此次异常地块共持续9天，自异常出现19天后西藏米林发生M_S6.9地震，震中位于2017年10月31日阈值零线5 km处。本次异常结果排除计算数据与高空电离层影响，表明其客观存在。而此次异常的分布面积呈增大—收缩—消失—增大—消失的特征，时间演化过程共持续9天（2017年10月30日至11月7日），具体参数列于表1，异常台站数量自2017年10月28日逐渐增加，于10月31日地磁极化超阈值异常面积达到最大，异常台站数量所占比例为48%，随后异常面积和异常台站数量减少，至11月6日又出现异常，异常台站数量所占比例为35%，11月7日极化异常消失（图5）。总体上，本次地磁极化异常分布范围从无到有、从小到大，达到异常台站数量、极化异常幅值的最大和次大值后异常消失，反映了区域应力的增强过程。

表  1  异常参数统计表

Table  1.  Statistical table of abnormal parameters

异常 日期	异常面积 /（104 km2）	异常台 个数	计算台 个数	异常 占比	异常最大 幅值	震中距 /km
10月28日	0	1	64	1.56%	0.02	1540
10月29日	0	0	64	0	-	-
10月30日	165.2	12	64	18.75%	0.40	380
10月31日	378.9	31	64	48.44%	0.87	1980
11月1日	70.9	8	65	12.31%	1.06	390
11月2日	0	0	65	0	-	-
11月3日	0	0	65	0	-	-
11月4日	2.3	1	63	1.59%	0.86	2300
11月5日	51.1	10	64	15.63%	1.26	970
11月6日	285.3	22	63	34.92%	1.22	2060
11月7日	0	1	64	1.56%	0.07	1540

下载: 导出CSV 

| 显示表格

图  5  异常参数统计

（a） 异常台站占比；（b） 异常面积；（c） 异常幅值最大值；（d） 震中距

Figure  5.  Statistical graph of abnormal parameters

（a） The proportion of abnormal stations；（b） Anomalous area；（c） Maximum anomalous amplitude；（d） The epicentral distance

下载: 全尺寸图片 幻灯片

地磁扰动会引起地磁观测资料波动变化。为了证明本次所研究的地磁垂直强度极化异常非磁暴所致，根据中国科学院国家空间科学中心环境预报中心（2020）提供的磁暴环电流指数D_st和磁情指数K_p变化情况（图6），异常出现期间地磁场活动较为平静（K_p＜5且D_st＞−30），故地磁扰动对此次异常的性质判定不产生影响。

图  6  极化高值异常期间D_st指数（a）和K_p指数（b）的变化情况

Figure  6.  Variation of D_st and K_p indices during polarization high value abnormal period

下载: 全尺寸图片 幻灯片

3.   讨论与结论

本文利用地磁垂直强度极化方法整理2017年计算结果，自2017年10月30日起中国大陆地磁观测台站出现大范围超阈值异常，超阈值异常分布范围呈现增大—收缩—消失—增大—消失的演化过程，且高值异常多次出现地区位于中国大陆西部，尤其青藏高原巴颜喀拉地块与羌塘地块交界处。异常变化共持续9天，自异常开始的19天后西藏米林发生M_S6.9地震，震中附近反复出现高值异常，且接近阈值零线。最近震中位置距2017年10月31日阈值零线仅5 km。本次异常分布范围从增大到减小，及至消失后又呈现出增大而后消失的过程，表明地震前的电磁现象是由多个孤立的电磁信号强激发源同时或不同时发射，并非单一激发源产生的电磁波（杜爱民等，2004），因而电磁扰动异常呈现出围绕在震源附近且片状分布的动态时空演化特征。众所周知，孕震过程中断层间的相互磨擦、挤压、形变和断裂等力学因素与构造活动息息相关，所以孕震期地壳的力学形变过程不可能仅发生在某一点上（杜爱民等，2004；苏维刚等，2020），而是点线面的组合形态。岩石力学的压力试验结果证明，当岩石接近破裂时，其电导率有所增加。根据电磁感应原理，地下电导率的变化必然会引起地磁感应场的变化（廖晓峰等，2019），孕震过程中断层间因磨擦、挤压、形变和断裂等力学因素会在孕震区附近产生临时感应电流，由于感应磁场的复杂性会导致异常的表现形式多样化，地磁垂直强度极化高值异常形态常在孕震区表现为时而间断、时而连续的片状分布，说明孕震过程中断层间因磨擦、挤压、形变和断裂等力学作用所产生的电磁异常信号是随机且不连续的，这可能是震前震源区附近观测台站不能一一接收到电磁异常信号的主要原因之一。

综合分析认为，地磁垂直强度极化高值异常与外源场无关，反映的是地下电磁信号。当地质块体因应力环境改变而发生或即将发生位移时，附近会产生临时感应电流，但受到地下介质的横向非均质性影响，感应磁场的表现形式多样化，故阈值零线附近会表现出不同相位、不同极化幅值的变化形态，地震也往往会发生在阈值零线附近。但是由于地震孕育过程的复杂性，地磁垂直强度极化异常与地震的发生并非一一对应的；并且，受限于不同的地磁台站布设间距，异常出现的概率不同，故单一震例无法代表该方法提取震磁异常的共性特征，但随着时间和资料的积累，判据指标及预测标准将得到进一步完善，期望能在未来实际地震预测工作中提供一定的参考价值。

审稿专家对本文提出了修改意见，作者在此表示感谢！