引言

地震孕育及发生过程中伴随着不同程度的电磁辐射异常（Molchanov et al，1992；Karakelian et al，2002；Matsushima et al，2002；汤吉等，2008，2010），天然地震发生前，也常常观测到电磁辐射异常信号（Han et al，2015；Wang et al，2018；Chang et al，2017；Nagao，2002；Uyeda et al，2002；Hattori，2004）。依据信号波长与传播距离的关系，周期介于几百秒到几百赫兹的信号不满足电磁波要求，地震研究人员将其称为地震电磁扰动（姚休义，冯志生，2018）。

在地震电磁扰动研究中，通过有效的解析方法能够从较强的电磁干扰背景中提取出相对较弱的地震电磁扰动信号。近年来，应用较为广泛的方法有极化法、主成分分析法、分形分析法和梯度定向法等。随着区域台站的增加和仪器采样精度的提高，各种解析方法均积累了很多成功的震例（杜爱民等，2004；Cianchini et al，2012；Hattori et al，2013；潘晖等，2014；李琪等，2015；袁桂平等，2015），其中能反映地震地磁扰动异常的产生、传播和分布特征的方法主要有地磁垂直强度极化法和水平椭圆极化法（Hayakawa et al，2000；Gladychev et al，2001；Molchanov et al，2003；Prattes et al，2008；冯志生等，2010；Li et al，2011）。

Hayakawa等（1996）利用震中距为65 km的日本关岛地磁台的磁通门磁力仪观测资料，分析了1993年8月8日关岛M8.0地震前后地磁场极化值Y_ZH的变化特征，结果显示震前两个月0.01—0.05 Hz （100—20 s）频段的Y_ZH逐渐增大，发震时达到最大，震后逐渐恢复。此后的众多震例研究验证了地震前数天至三个月内地磁场存在高极化值异常现象（Hayakawa et al，2000；Ismaguilov et al，2001；Molchanov et al，2003；Hattori，2004；Prattes et al，2008；李琪等，2015；何畅，冯志生，2017；姚休义，滕云田，2017）。冯志生等（2010）的研究表明，垂直强度极化存在季节变化，提取异常前需予以剔除。上述大部分研究仅涉及震中周围部分台站的观测数据。Hattori （2004）的研究涉及两个台站的数据，最大震中距为118 km；李琪等（2015）利用2009年9月至12月云南西部地磁台阵的连续观测数据，分析了距离2009年11月2日宾川M5.0地震震中143 km范围内的四个地磁台资料，由于台站少，空间特征研究仅限于异常幅值的线性拟合分析。许多研究甚至仅采用单台数据进行异常分析，无法给出异常的空间分布特征。

2015年前后，中国地震局国家地磁台网中心新增了大量磁通门观测仪器，不仅大幅度提升了南北地震带地磁场监测的密度，填补了青藏高原及新疆部分地区地磁监测的空白，而且将采样频率从1 min提高至1 s，使得整个中国大陆西部强震多发地区的地磁场观测水平得到了大幅提升。上述磁通门秒采样观测资料为开展不高于1 Hz的地磁场信号分析工作提供了数据支持。

中国大陆的强震主要分布于西部，对西部地区开展极化异常时空演化特征分析具有十分重要的现实意义。鉴于此，本文拟利用极化法对110°E以西区域的磁通门秒采样观测数据进行处理分析，获取震前地震磁扰动异常信号，并分析其时空演化特征，旨在为该地区后续震情跟踪分析提供参考。

1.   数据及分析方法

1.1   地磁秒采样观测数据

中国地震局国家地磁台网中心总计约150个地磁观测站，观测仪器主要有磁通门磁力仪、质子矢量磁力仪以及Overhauser磁力仪，大部分台站架设了两台以上的观测仪器，以保证数据记录的完整性。其中，磁通门磁力仪产出的秒采样数据可用于地磁超低频（ultra-low frequency，缩写为ULF）异常分析。2011年，国家地磁台网中心仅有磁通门磁力仪112套；2015年前后，随着大规模的台网建设，磁通门磁力仪增至191套。截至2019年，国家地磁台网中心配备磁通门磁力仪189套（图1）。

图  1  2011年至2019年国家地磁台网中心磁通门磁力仪数量

Figure  1.  The number of fluxgate magnetometers in National Geomagnetic Networks Center of China from 2011 to 2019

下载: 全尺寸图片 幻灯片

目前，中国地震局国家地磁台网中心在110°E以西区域的76个台点共安装120套磁通门磁力仪。为了研究大范围地磁场ULF信号特征，本文采用台站分布相对合理（图2）、观测资料较为丰富的2015—2019年数据进行分析。

图  2  2019年中国西部磁通门磁力仪（秒采样）分布

Figure  2.  Distribution of fluxgate magnetometers （second sampling） in western Chinese mainland in 2019

下载: 全尺寸图片 幻灯片

1.2   地震事件

根据中国地震台网地震目录，2015年至2020年1月，中国大陆西部共发生M6.0以上强震11次。除此之外，本文将中国大陆边界附近发生的4次M6.0以上强震也纳入分析范围，但不包括其强余震（表1）。

表  1  2015年至2020年1月中国大陆西部及周边M≥6.0强震信息

Table  1.  Information of strong earthquakes （M≥6.0） during the period from 2015 to January of 2020 in western Chinese mainland and its surrounding areas

发震时刻 　年-月-日　时:分:秒	M	北纬/°	东经/°	震源深度/km	参考位置
2020-01-19　21:27:55	6.4	39.83	77.21	16	新疆伽师县
2019-11-21　07:50:45	6.0	19.50	101.35	10	老挝
2019-06-17　22:55:43	6.0	28.34	104.9	16	四川长宁县
2019-04-24　04:15:48	6.3	28.40	94.61	10	西藏墨脱县
2017-11-18　06:34:19	6.9	29.75	95.02	10	西藏米林县
2017-08-09　07:27:52	6.6	44.27	82.89	11	新疆精河县
2017-08-08　21:19:46	7.0	33.2	103.82	20	四川九寨沟县
2016-12-08　13:15:03	6.2	43.83	86.35	6	新疆呼图壁县
2016-11-25　22:24:30	6.7	39.27	74.04	10	新疆阿克陶县
2016-10-17　15:14:49	6.2	32.81	94.93	9	青海杂多县
2016-06-26　19:17:11	6.7	39.43	73.4	10	吉尔吉斯斯坦
2016-01-21　01:13:13	6.4	37.68	101.62	10	青海门源县
2015-12-07　15:50:03	7.4	38.2	72.9	30	塔吉克斯坦
2015-07-03　09:07:46	6.5	37.6	78.2	10	新疆皮山县
2015-04-25　14:11:26	8.1	28.2	84.7	20	尼泊尔

下载: 导出CSV 

| 显示表格

1.3   极化分析

提取地震磁扰动信号的关键是突出源于震源区的异常信号，压制源于外源场的信号。Hayakawa等（1996）由此发展出了地磁垂直分量与水平分量幅值比分析方法：

式中，Z（ω）为地磁垂直分量的傅里叶谱值，H（ω）为地磁水平分量全矢量的傅里叶谱值，H_X（ω）为地磁水平分量南北向傅里叶谱值，H_Y（ω）为地磁水平分量东西向傅里叶谱值，地磁各分量单位为nT，ω为圆频率。由于该方法反映了磁场强度或能量在垂直方向的分配比例，将其简称为地磁垂直强度极化法，将地磁垂直分量与水平分量的幅值比Y_ZH称为地磁垂直强度极化值。

1.4   信号处理

根据趋肤效应，可以推算源于震源处的信号频率范围。趋肤深度定义为

式中δ为趋肤深度，ρ为电阻率，μ为磁导率，μ＝μ₀＝4π×10⁻⁷ H/m，μ₀为真空中的磁导率，ω为角频率。由于不能准确地测得震源处的电阻率，我们用沉积岩的电阻率ρ＝100 Ω·m代替。根据表2所示的2015年以来中国大陆西部Y_ZH高值异常信息，本研究的地震震源深度处于6—30 km范围内。将上述参数代入式（2），计算得到ω的范围为4.42—0.18 Hz。考虑到资料采样率的局限性以及电阻率的不确定性，并综合前人的研究成果（冯志生等，2010；李琪等，2015），我们将频带范围限定到0.01—0.2 Hz。数据处理的步骤如下：① 将磁通门三分量Z，H_X和H_Y每天的秒资料分为96段，每段15 min，分段计算垂直分量和水平分量的傅里叶谱值；② 计算每天5—100 s （0.01—0.2 Hz）周期内各频点的极化值均值，获得极化值的逐日变化序列；③ 对极化值逐日变化序列进行傅里叶拟合（数据一般不短于半年），获得周期大于半年的傅里叶拟合变化曲线，即为极化值的年变化曲线；④ 将5—100 s周期内各频点极化值减去年变化后，剔除低于二倍均方差的部分，仅保留高值部分；⑤ 对筛选出的极化高值的日均值逐日变化序列再次进行傅里叶拟合，如果仍然存在年变化成分，则再次扣除极化值高值周期大于半年的成分，所得残差为最终处理得到的极化高值序列，即Y_ZH逐日值。

表  2  2015年至2019年中国大陆西部$Y_{ZH2}$高值异常信息

Table  2.  Information of high Y_ZH2 value anomalies in western Chinese mainland from 2015 to 2019

序号	异常出现时间 　年-月-日	异常持续时间/d	异常台站最多日（异常日） 　年-月-日	异常台站数	总台数	YZH2
1	2015-01-16	6	2015-01-18	28	36	0.73
2	2015-02-09	9	2015-02-13	14	36	0.39
3	2015-07-17	5	2015-07-18	8	36	0.29
4	2015-09-26	5	2015-09-29	13	27	0.47
5	2015-11-18	11	2015-11-23	21	27	0.92
6	2016-01-27	8	2016-01-30	17	41	0.75
7	2016-08-14	9	2016-08-17	15	33	0.92
8	2016-09-11	9	2016-09-15	14	33	0.88
9	2016-12-01	6	2016-12-02	14	37	0.75
10	2017-01-13	5	2017-01-15	8	38	0.41
11	2017-02-13	5	2017-02-14	10	37	0.49
12	2017-03-11	9	2017-03-18	14	34	0.66
13	2017-08-24	6	2017-08-28	22	40	0.65
14	2017-10-29	6	2017-10-30	14	40	0.70
15	2018-04-16	3	2018-04-17	11	43	0.62
16	2018-06-10	6	2018-06-12	17	43	0.44
17	2019-09-22	3	2019-09-24	13	54	1.15
18	2019-11-12	5	2019-11-12	26	56	0.85

下载: 导出CSV 

| 显示表格

以往震例（姚休义，冯志生，2018）表明，震前的超低频异常往往持续数天，且不同台站出现时间可能略有差异。虽然这种差异并不会影响到时间演化过程，但在研究异常空间特征时，可能会因为前后出现时间的差异而导致异常的空间范围与实际不相符，因此对逐日值进行5日滑动平均处理，将高值异常进行时间上的延拓，得到Y_ZH1。

为便于异常空间特征的分析，利用各台Y_ZH1的二倍均方差对其进行归一化，对归一化后的残差曲线减1，得到Y_ZH2，该值被称为归一置零极化值。具体算法为：

Y_ZH2高于0，等价于Y_ZH1高于2倍方差。图3中展示了喀什台2017年Y_ZH1与Y_ZH2的变化曲线，可以看到两者具有相同的变化趋势。 利用现有台站的归一置零极化值，在空间上进行简单的数学插值（克里金法），即可得到大陆西部地区每日的归一置零极化值的空间分布。

图  3  喀什台2017年Y_ZH1与Y_ZH2的对比

Figure  3.  Comparison of Y_ZH1 with Y_ZH2 results of Kashi station in 2017

下载: 全尺寸图片 幻灯片

2.   结果分析

2.1   极化法计算结果

对2015年以来的76个台点120套观测仪器所产出的数据进行上述处理后，每套观测资料计算得到一个Y_ZH1序列。因观测环境变化或仪器故障等原因，每年能够正常产出数据的台站数出现小幅波动。

Y_ZH序列与Y_ZH1序列在中国大陆西部各台具有较好的同步性。以2017年为例（图4），由南向北分布的通海台（THA）、西昌台（XIC）、江油台（JOY）、天水台（TYS），以及由西向东分布的喀什台（KSH）、且末台（QMO）、都兰台（DUL）、英鸽台（YGE），其Y_ZH与Y_ZH1序列均具有相同的变化趋势，但变化幅度未呈现随经度或纬度规律变化的特征。

图  4  沿经线和纬线分布的部分台站及其Y_ZH和Y_ZH1 2017年结果

Figure  4.  Location of some stations and their Y_ZH and Y_ZH1 results in 2017

下载: 全尺寸图片 幻灯片

基于以上规律，对Y_ZH1序列进一步处理，并再一次筛选参与分析的台站。因原始数据存在一些干扰、突跳甚至错误，Y_ZH1时序曲线出现显著偏离背景的情况，会出现与周边台站出现显著的不一致现象，因此在确认原始数据出错的情况下，删除显著偏离背景时段的Y_ZH1结果。对于有两套以上仪器的台站，仅选择数据连续性更好、与周边台站同步性更高的一套参与后续的空间特征分析。

2.2   极化高值与地磁扰动

为了研究极化值与外源场扰动之间的关系，我们在东京地磁国际数据中心网站http://wdc.kugi.kyoto-u.ac.jp下载了2015年至2019年间的地磁K_p指数和D_st指数，并将两者与Y_ZH （逐日值）进行相关性分析。计算结果显示每天的K_p指数最大值与Y_ZH的相关系数为−0.449，表明K_p指数与Y_ZH存在弱的负相关，即K_p指数高值与Y_ZH低值可能存在弱的相关性，这与冯志生等（2010）的研究结果一致，即“当K指数为高值时，极化结果为低值”。这一结果同时也符合该方法的预期，即来自电离层外空磁场一次源的垂直分量幅值与水平分量幅值的比值较小。此外，D_st指数最大值与Y_ZH的相关系数为0.176，表明两者不相关。

从2017年K_p指数、D_st指数与喀什台Y_ZH1结果（图5）来看，极化高值往往出现在K_p指数小于3的时段和D_st指数处于±20之间的时段，即地磁扰动较弱的时段。

图  5  喀什台2017年Y_ZH1序列与K_p （a）和D_st （b）指数曲线的对比

Figure  5.  Comparison of Y_ZH1 series with K_p （a） and D_st （b） curves in Kashgar station in 2017

下载: 全尺寸图片 幻灯片

2.3   中国大陆西部超低频高值异常

Y_ZH1单台高值现象较多，7个及以下台站出现高值时，高值台站也往往孤立分布，其可靠性有待验证。因此本文仅讨论2015年至2019年中国大陆西部出现的8个以上台点同步Y_ZH1高值现象。

表2给出了2015年至2019年中国大陆西部共出现的18次8个以上台点同步Y_ZH2高值现象，图6对每月出现高值异常的频次进行了统计，结果表明：异常出现频次最高的月份为1月和9月，均为3次；其次是2月、8月、11月，分别出现2次；除5月外，其余各月均有1次高值异常。总体呈春夏低秋冬高的趋势，但并未呈明显的优势分布。

图  6  Y_ZH1高值异常在12个月中的频次分布

Figure  6.  The frequency distribution of Y_ZH1 high value anomalies in 12 months

下载: 全尺寸图片 幻灯片

图7给出了2015年至2019年期间Y_ZH1高值频次与地震频次的对比，可见在Y_ZH1高值频次较低的2018年，后续无M6.0以上强震发生。

图  7  Y_ZH1高值异常频次与地震频次图

Figure  7.  The frequency of Y_ZH1 high value anomalies and earthquake frequency

下载: 全尺寸图片 幻灯片

为便于后续描述，我们按时间顺序分别将18次高值现象编号为1—18，并取高值台站最多日作为研究日，对上述高值异常的空间分布特征进行研究。

图8给出了2015年至2019年期间部分台站的Y_ZH1时序图，图9和10给出了18次高值异常的空间分布图。从空间等值线图不难看出，每次高值现象出现时，高值区域的位置和面积均不相同：研究区域东部的高值区分布特征较为复杂，总体来说高值区域面积较小；西部地区由于台站稀疏，更容易出现大面积高值区，这种大面积高值可能与插值有关，但也从一定程度上反映了相关区域实际Y_ZH2偏高的情况。从时间变化过程来看，2015年1月、11月和2016年9月高值区面积较大（图9），2018年的两次高值区面积较小（图10），其余高值区的面积相差不明显。

图  8  2015年至2019年部分台站Y_ZH1时序曲线（图中数字序号与表2中的异常序号相对应）

Figure  8.  Y_ZH1 time series curves of some stations from 2015 to 2017 （The numbers correspond with the series numbers of the anomalies in Table 2）

下载: 全尺寸图片 幻灯片

图  9  2015年至2016年磁场高值日Y_ZH2的空间分布图

Figure  9.  Spatial distribution of Y_ZH2 on the magnetic field high value day from 2015 to 2016

下载: 全尺寸图片 幻灯片

图  10  2017年至2019年磁场高值日Y_ZH2的空间分布图

Figure  10.  Spatial distribution of Y_ZH2 on the magnetic field high value day from 2017 to 2019

下载: 全尺寸图片 幻灯片

为了进一步验证高值现象与强震相关的可信度，本文选择中国东部的静海台、高邮台与西部的喀什台进行对比。以2017年为例（图11），喀什台的Y_ZH1随时间变化的趋势与高邮台大体相同，但不同时段的相对变化幅值差异较大。图中红色虚线标示出Y_ZH1出现高值异常的时间，喀什台Y_ZH1序列高值时间与其基本吻合，而东部的高邮台Y_ZH1序列在上述时间并无明显的相对高值出现。计算上述两台站的Y_ZH1序列相关系数为0.26，表明两台站资料存在弱相关性；而静海台与喀什台的曲线差别较大，Y_ZH1序列相关系数为0.05，表明两者不相关。

图  11  喀什台、静海台和高邮台的空间分布（a）以及2017年极化Y_ZH1时序曲线对比图（b）

图中红色竖线标示中国大陆西部Y_ZH1出现高值异常的时间

Figure  11.  The locations of the stations Kashi，Jinghai and Gaoyou （a） and Y_ZH1 time series curves in 2017 （b）

The red vertical line indicates the time when the Y_ZH1 high value anomaly appeared in western Chinese mainland

下载: 全尺寸图片 幻灯片

2.4   强震与Y_ZH高值的时空关系

2.4.1   时空关系

2015—2019年中国大陆西部共出现18次8个以上台站Y_ZH1同步高值现象，共发生15次M6.0以上强震；2015—2017年内共出现14次多台同步高值现象，发生11次M6.0以上强震；2018—2019年共出现4次高值现象，较2015—2017年频次明显降低，强震频次也有所降低，共发生4次强震。Y_ZH高值频次与地震频次的相关系数为0.47，即二者存在弱相关。

对2015年1月至2020年1月中国大陆西部及周边15次强震发生前半年的Y_ZH1高值现象进行回溯性研究，结果显示，除2019年4月24日墨脱M6.3地震和2019年6月17日长宁M6.0地震以外，其余13次强震发生前半年内，Y_ZH1均出现多台同步高值。上述13次强震中，除2016年11月25日新疆阿克陶M6.7地震外，其余12次强震震中及周边均存在高值区。

值得注意的是，在2015年7月3日皮山M6.5地震后15天、2016年1月21日门源M6.4地震后9天、2016年12月8日呼图壁M6.2地震后37天、2017年8月9日精河M6.6地震后19天以及2017年8月8日九寨沟M7.0地震后20天内，出现了Y_ZH1多台同步高值现象且地震震中位于高值区（表3）。

表  3  单次高值异常及短期内多次高值异常叠加后发生强震的高值区面积

Table  3.  The high value areas of single high-value anomalies and the superposition of several high-value anomalies in a short period which followed with strong earthquakes

序号	出现高值日期 年-月-日	震后出现高值区的 地震（提前的天数）	半年内高值区强震	高值区面积① /（104 km2）	高值区总面积② /（104 km2）
1	2015-01-18	无	2015-07-03皮山M6.5	251.6	313
2	2015-02-13		2015-04-25尼泊尔M8.1
3	2015-07-18	2015-07-03皮山M6.5 （15）	2015-12-07塔吉克斯坦M7.4	13	20.6
4	2015-09-29	无	2016-01-21门源M6.4	20	144
5	2015-11-23		2015-12-07塔吉克斯坦M7.4	123.7
6	2016-01-30	2016-01-21门源M6.4 （9）	2016-06-26吉尔吉斯斯坦M6.7	17.5	84.7
7	2016-08-17	无	2016-12-08呼图壁M6.2	24.1	82.2
8	2016-09-15		2016-10-17杂多M6.2	58
9	2016-12-02	无	2016-12-08呼图壁M6.2	18.2	80.8
10	2017-01-13	2016-12-08呼图壁M6.2 （37）	无	47.2	47.2
11	2017-02-14	无	2017-08-08九寨沟M7.0	41.6	84.4
12	2017-03-18		2017-08-09精河M6.6	30.5
13	2017-08-28	2017-08-08九寨沟M7.0 （19） 2017-08-09精河M6.6 （20）	无		39.5
14	2017-10-31		无		16.1
15	2018-04-17	无	无		32.4
16	2018-06-12
17	2019-09-24	无	2019-11-21老挝M6.0	19.2	72.3
18	2019-11-12
注：① 指地震震中所在高值区的面积；② 指一次高值现象出现的多个高值区面积之和。

下载: 导出CSV 

| 显示表格

2.4.2   短期内反复出现高值的区域易发强震

上述研究还存在一个现象，Y_ZH2多台高值现象可能在短期内多次出现，且后续强震的震中位置往往为多次出现高值的区域。为了进一步分析这种现象，我们将短期（两个月）内出现且高值区有重叠的多次Y_ZH2高值进行简单叠加，具体做法为：对同一点的N次多台高值日Y_ZH2求平均，经过上述处理，多次均为高值的区域仍为大于0的高值区。但为了排除可能的震后效应产生的影响，强震发生后震中附近区域出现的高值不叠加。

经过筛选，短期内（两个月）出现且高值区有重叠的Y_ZH2多台同步高值现象共有6组，序号分别为①-②，④-⑤，⑦-⑧，⑩-⑪，⑮-⑯，⑰-⑱，叠加结果如图12所示。叠加过程使得多次高值的区域突出，高值区面积缩小。不难看出，后续强震震中大部分位于多次叠加后的高值区域及其边缘，2016年1月25日阿克陶M6.7地震和2019年11月21日伽师M6.0地震除外。其中阿克陶地震前出现了两次高值现象，但在震中附近均未出现高值区，而伽师地震前58天在震中出现了高值区。

图  12  2015—2019年6组高值异常叠加前、后的Y_ZH2与后续6个月内地震的空间关系

Figure  12.  Spatial relationship of six groups of high Y_ZH2 anomalies before and after superposition as well as the distribution of earthquakes in the following six months from 2015 to 2019

下载: 全尺寸图片 幻灯片

2.4.3   一次或一组多台高值现象中多个高值区发生强震

图12显示，一次高值现象中会出现多个Y_ZH2高值区域。多次高值现象累加后，也产生多个高值区域。从上述研究中我们不难发现，一次多台高值现象中的高值区或一组高值现象叠加后出现的多个高值区均有可能发生强震，一个高值区内也可能发生多次强震。这种现象在叠加后的结果中表现得更明显，例如：①与②的叠加结果中，尼泊尔地震和皮山地震位于同一个高值区内；④与⑤的叠加结果中，塔吉克斯坦地震和门源地震的震中分别位于两个高值区内；⑦与⑧的叠加结果中，呼图壁地震和杂多地震的震中分别位于两个高值区内；⑩与⑪的叠加结果中，九寨沟地震和精河地震的震中分别位于两个高值区内。

图  12  2015—2019年6组高值异常叠加前、后的Y_ZH2与后续6个月内地震的空间关系

Figure  12.  Spatial relationship of six groups of high Y_ZH2 anomalies before and after superposition as well as the distribution of earthquakes in the following six months from 2015 to 2019

下载: 全尺寸图片 幻灯片

2.4.4   Y_ZH2高值区面积与后续强震震级的关系

我们基于叠加结果，将18次高值异常重新分组，产生12组高值，每组均有一幅Y_ZH2空间分布图。统计单张Y_ZH2空间分布图中高值区的总面积及后续半年内发生强震的高值区域面积（表3），对上述两个参量与后续地震震级的关系进行了分析，结果如图13所示。可以看出，高值区面积有随震级增大而增大的趋势，相关系数为0.89。

图  13  叠加后的Y_ZH2高值区面积与后续6个月内地震震级的关系

Figure  13.  The relationship between Y_ZH2 high-value areas of the superposition of several high-value anomalies and the magnitude of strong earthquakes in six months

下载: 全尺寸图片 幻灯片

3.   讨论与结论

本文基于2015至2019年期间中国大陆西部磁通门秒采样观测资料，开展了5—100 s频段地磁垂直强度的极化分析。分析结果表明，中国大陆西部各台站的极化结果具有较好的一致性，沿经向和纬向均无明显的形态及幅值变化，往往出现多台同步高值现象。虽然出现同步高值的台站较多、空间面积较大，但对极化结果与K_p指数及D_st指数的互相关系数分析结果表明，极化高值与外源场扰动的关系不大。对Y_ZH2高值与后续强震的时空关系研究显示：Y_ZH2高值现象出现后的半年内，高值区域及周边容易发生强震；有部分地震后短时间内，震中周边也出现了极化高值；高值现象出现后，多个高值区均有可能发生强震；后续强震的震级与高值区面积呈正相关。值得注意的是，上述结果仅仅是基于5年的数据所得，因统计样本有限，可能不能完全反映客观规律。随着资料的进一步积累，将进一步加深对地磁垂直强度极化高值异常与地震关系的认识。

尽管地磁垂直强度极化法经过了很多震例的验证，但仍有研究人员质疑该方法提取震磁异常信息的可靠性（Thomas et al，2009；Masci，2011），主要存在二点疑问：第一，因磁层电流体系的影响无法排除，认为极化值升高很可能与外源场变化有关；第二，以往研究中震前异常通常仅存在于震中距最小的单一台站，在无法排除数据干扰的情况下，异常的可靠性也常被质疑（姚休义，冯志生，2018）。

关于第一点，本文中进行了简单的讨论，认为极化高值与外源场扰动的关系不大。冯志生等（2010）也曾就此问题展开研究，其结果表明地磁垂直强度极化高值与地磁活动性指数之间存在较好的负相关性。Currie和Waters （2014） 对各频段极化值与K_p指数、AE指数及SYM-H指数的相关性分析结果表明，各频段极化值与上述各种地磁活动性指数之间的相关系数均低于0.3，即两者并不相关。

关于第二点，本文开展了大区域多台站极化高值与地震相关性的研究，结果显示高值往往在较大范围内的多个台站同时出现，与台站附近的干扰无关。李琪等（2015）利用极化法在多个台站均提取到了震前极化异常，且异常特征符合地震电磁扰动信号的衰减特征（Huang，Ikeya，1998）。上述研究均表明极化值升高的结果真实可信。

但与以往研究结果不同的是，本研究显示震前出现高值的空间范围不局限于数十千米，时间也不局限于数日之内，有可能达到半年左右。梅世蓉（1996）在华北地区开展的研究表明，一次M6.0以上地震前常出现持续时间长达数十年的椭圆形活动增强区，其长度至少长达160 km，M8.0以上强震前甚至能够达到1 000 km以上。而活动增强区的微震破裂很可能产生电磁辐射并出现极化高值。

同时出现的多个异常区有可能发生多次强震。根据马宗晋（1980）的研究结果，震前区域应力会在多点集中，跨区域的多次强震之间很可能存在联系而出现强震连发的情况。如果这种联系是存在的，那么在某种条件的激发下，就有可能在多个区域同时出现极化高值。而同时出现的多个极化高值区后续是否发震，取决于每个孕震体的发育情况。震级与高值区面积呈正相关的结果也支持上述推论。

此外，震后一个月左右震中附近有可能出现极化值再次升高的情况。关于这一现象的机理尚需要进一步研究。

在数据下载及处理过程中，上海市地震局朱培育高级工程师提供了软件支持；文中插图主要使用GMT （Wessel，Smith，1995）绘制，作者在此表示衷心的感谢。