2022年1月8日青海门源6.9级地震短临异常跟踪分析

冯丽丽 李霞 赵玉红 张朋涛

1)中国地震局青海格尔木青藏高原内部地球动力学野外科学观测研究站，西宁 810001 2)青海省地震局，西宁 810001

极化法(Polarization)由Hayakawa等(1996)提出，因其能够反映地震地磁扰动异常的产生、传播和分布特征，被广泛应用于国内外地震电磁扰动异常提取(Hayakawa et al，2000；
Molchanov et al，2003；
Prattes et al，2008；

冯志生等，2010；

李琪等，2015；

廖晓峰等，2019；

艾萨·伊斯马伊力等，2020；

樊文杰等，2021)。但以往研究主要以单台资料分析为主，对震前电磁辐射异常的空间特征研究较少。2019年，基于2015—2018年中国地磁台网磁通门磁力仪产出的秒采样观测数据，在对分析结果进行了适当的处理后，总结了地磁垂直强度极化高值异常与后续地震的时空关系，建立了完整的地震预报指标体系(中国地震局监测预报司，2020)，使得该方法得到了较大发展。贺曼秋等(2019)在2017年九寨沟7.0级和精河6.6级地震前，在2个震中附近同时提取到了地磁垂直强度极化异常；

李霞等(2021)在2017年西藏米林6.9级地震前10天提取到了青藏高原出现的大范围地磁垂直强度极化高值异常；

管贻亮等(2021)对2015年尼泊尔地震前的地磁垂直强度极化异常进行了分析并从构造动力演化角度对地磁极化异常的孕育和发展进行了解释；

冯丽丽等(2021)总结了2015—2020年中国大陆西部地磁垂直强度极化异常和地震的关系，认为地磁垂直强度极化异常出现后半年内，空间分布图中0.2阈值线附近发生6级以上强震的可能性较大，且高值异常区面积与震级存在正相关关系。

地震前也常观测到地下流体的异常变化。在2021年青海玛多7.4级(钟骏等，2021)、2017年新疆精河6.6级(钟骏等，2018)、2010年青海玉树7.1级地震(李滔等，2012)等中国大陆西部强震前，均在震中周边多个台站观测到了地下流体化学组分的浓度变化、温度变化或水位变化等。

Hayakawa等(1996)提出的地磁垂直强度极化法能够突出源于震源区的异常信号，抑制来源于外源场的信号。具体数学表述如下

(1)

式中，Z(ω)为地磁垂直分量的谱幅度值，H(ω)为地磁水平分量的谱幅度值；
HX(ω)为地磁南北分量谱幅度值，HY(ω)为地磁东西分量谱幅度值。地磁分量单位为nT，ω为圆频率。

1.1 数据处理

以往震例表明，震前的极化异常在不同台站出现时间可能略有差异，在研究异常空间特征时，可能因为前后出现时间的差异导致异常的空间范围与实际不相符的情况，因此对逐日值进行了5日滑动平均处理，得到YZH1。

1.2 归一化

因观测环境及仪器的差异，各台YZH1波动范围不尽相同。在提取YZH1高值进行空间特征研究时，可能出现因台站差异导致的失真现象。为了消除台站影响，需要对各台YZH1进行归一化。为了使结果更直观，对归一化后的结果进行了置零处理。具体算法为

(2)

我们称YZH2为归一置零极化值，当其高于0时，等同于YZH1高于2倍方差。在空间上对现有台站的归一置零极化值进行克里金法数学插值，即可得到大陆西部地区每日的地磁垂直强度极化空间分布。

2.1 地磁垂直强度极化异常

2021年10月27日前后，全国范围内77个参与计算台站中共有17个台站出现高值异常。从青海地区的YZH1时序曲线来看(图1)，从2021年10月20日起青海地区多个台站的YZH1出现持续升高，至10月27日前后达到极大值，后持续回落。这是2021年6月(玛多7.4级地震)以来青海地区最显著、幅度最大的一次地磁垂直强度极化高值异常。

图1 2021年1月以来青海地区各台地磁垂直强度极化时序曲线

图2 2021年10月27日地磁垂直强度极化YZH2异常空间分布

经过后续处理和分析，发现在异常幅度最大、异常台站最多的时间，即2021年10月27日，在青海西北部形成1个面积大于6×104km2的高值异常区(图2)，异常面积为6.6×104km2。大于阈值0.2的范围内有3个异常台站，分别为德令哈台、贵德台和金银滩台，详细信息见表1。其中归一置零极化值最高的为金银滩台，达到0.9。

表1 2021年10月27日青海地区垂直强度极化信息

2.2 地下流体准同步短期异常

2021年7月至11月，青海地区陆续出现多项流体短临异常。其中佐署台动水位的下降异常及浅层水温的上升异常出现时间最早，2个测项出现异常的时间一致，均为7月初。而平安静水位破年变异常随后在7月底至8月初出现。门源静水位突降异常、乐都逸出气氡突升异常以及西宁逸出气氡浓度突降异常均出现在2021年9月。上述异常大多仅持续较短的时间，而西宁逸出气氡浓度突降异常一直持续到11月初，几乎与此同时，共和的2套中层水温也同时出现了突降异常，并于12月初再次出现同步变化。德令哈水温在12月上、中旬出现突升—转平—突降的“倒U”形变化，整个变化过程持续约半个月。至此，青海地区地下流体异常数量达到8项(共和2套水温观测仪器相同且装置于同一井孔中，因此计为1项)，空间分布见图3，时序曲线见图4。

图3 门源地震前青海地区地下流体异常台站分布

图4 2021年6月以来青海多个台站地下流体准同步异常

从震中距及异常出现时间来看，最远的德令哈台和相对稍远的共和台出现异常的时间最晚。而最近的门源台于震前4个月出现异常，在此次准同步变化过程中处于中段。西宁台、佐署台、平安台、乐都台及共和台震中距相差不大，但异常出现的时间最大相差了5个月，无明显的时空演化特征。

从异常类别来看，共3个台出现水位变化，3个台的4套仪器观测显示水温变化，2个台出现逸出气氡浓度变化。其中平安台、佐署台和门源台出现了水位变化，平安台打破了原本的年变化上升趋势，水位转平下降；

门源台出现了水位上升变化(平安台与门源台为静水位观测，数值下降表明水位上升)；

而佐署台则呈现了持续性的下降变化。水温变化在德令哈台和佐署台主要表现为升高变化，而共和台2套仪器的观测结果均显示阶降变化。而从逸出气氡浓度来看，乐都台出现稍早，主要为上升变化，在乐都台开始下降阶段西宁台出现了上升变化。

青海地区地下流体测项在2021年5月22日玛多7.4级地震前1～3个月，曾提出多项准同步异常变化，包括玉树水温(苏维刚等，2021)、佐署动水位、共和逸出气氡及共和中层水温等(刘磊等，2021)。2016年门源6.4级地震前1～3个月，也曾出现过3项准同步变化，包括玉树水温“V”型变化、乐都气氡快速升高异常(刘磊等，2020)、平安静水位破年变变化等(杨晓霞等，2016)，部分异常台站和测项与此次异常特征相似。且对应地震的时间均在3个月内。

3.1 时间判定

表2 列出了青海地区地磁垂直强度极化异常信息。自2015年以来，青海地区共出现8次地磁垂直强度极化异常，其中6次在异常出现后8个月内异常区内及附近发生了5.6级以上地震，对应率为75%，具有较高的信度。而其中4次地震发生在异常出现后的60天内，表明该异常具有一定的短期指示意义。从异常出现时间推测优势发震时间在2021年12月27日前。但考虑到目前震例较少，且部分震例对应时间较长，因此在发震时间判定中同时考虑了地下流体准同步变化的时间。2021年7月以来，青海地区陆续出现多项流体异常。2021年11月再次出现多项准同步变化，这种现象对于后续3个月内青海地区的强震有较好的指示意义。综合上述分析，判断发震时间为2021年12月至2022年1月。

表2 青海及周边地区极化异常信息

图5 青海地区地磁垂直强度极化YZH2异常空间分布

3.2 地点判定

图5给出了青海及周边地区垂直强度极化异常空间分布与后续地震的震中位置，震中全部位于异常区内或阈值线(归一置零极化值为0.2的等值线)附近。因此，预测发震地点在异常区内或阈值线附近。

3.3 震级判定

由表2 可知，青海地区震例震级在5.6～7.4级之间。根据以往对中国大陆西部的震例研究(冯丽丽等，2021)，后续地震震级与高值异常区面积有近似正相关关系，并存在拟合关系。此次异常区面积与夏河5.7级地震前异常区面积相当。综合分析认为，目标地震可能在6级左右，因此最终给出预测震级为5.6～6.4级。

3.4 综合预测结果

基于青海及周边地区地磁垂直强度极化异常，并综合2021年7—11月以来青海地区地下流体准同步异常变化，基于以往震例分析认为2021年11月23日至2022年1月23日，青海西北部地磁垂直强度极化高值区内可能发生5.6～6.4级左右地震，并填报地震短临预测卡片。此次门源6.9级地震发生于异常后的73天，震中位于预测区边界，发震时间和震中位置均符合预期，但震级比预期高0.5级。

(1)此次门源6.9级地震前73天在青海西北部出现了地磁垂直强度极化异常，高值区面积为6.6×104km2，空间上位于此次地震震中西南。震前，通过对该异常在青海地区的预测指标进一步分析认为异常出现后60天至稍长时间，高值区附近可能发生6级左右地震。

(2)此次门源6.9级地震前1～6个月，青海地区出现了多项地下流体异常，在时间上表现为准同步变化。根据异常指标，认为发震时间可能在2021年12月至2022年1月。

(3)此次门源6.9级地震的发生表明震前的分析和判断基本正确。地磁垂直强度极化异常空间阈值线是地震地点预测的良好指标。该异常在青海地区可作为短临异常，指示未来2～3个月内的中强以上地震，地下流体的准同步异常具有较强的时间指示意义。

优质的时间、空间预测指标是本次地震短临跟踪分析成功的关键。但同时应指出，此次预测的震级仍存在0.5级的偏差，主要原因是目前在青海地区尚未梳理出高质量的震级预测指标，需要在下一步工作中继续加强此方面工作。