吴忠—灵武地区MS≥5.0地震前波速比变化特征分析

李文君 曾宪伟 任雪梅 许英才 朱鹏涛 李自芮 马小军

（中国银川750001 宁夏回族自治区地震局）

0 引言

地震波在传播过程中携带了地球内部的大量信息，可为研究地下介质变化提供基础数据。在地震孕育过程中，尤其是在临震阶段较短时间内，震源及其附近地区的地下介质由于区域应力场增强、孕震区应力增加和应变积累，会出现诸如微破裂、塑性硬化、裂隙密度及饱水层变化等现象，当地震波通过该介质及周边区域时，地震波的传播速度会发生变化。早在1956 年，日本地震学家宫本贞夫就发现了强震前地震波速度变化的异常现象（冯德益，1981）。李艳娥等（2014）研究发现，汶川MS8.0 地震前波速比显著变化台站的分布范围与孕震区范围基本一致，为确定此次地震前震源区是否存在地壳介质的显著时空变化提供了有利证据；王林瑛等（2014）对芦山MS7.0 地震前波速比变化进行回溯性分析，发现部分台站出现波速比低值异常随时间变化的现象，该研究结果对于判定未来强震孕育可能的展布空间具有指示意义。

近年来，随着数字地震台站逐渐增多，吴忠—灵武地区积累了大量中小地震数字化观测资料。值得关注的是，该区2020 年接连发生5 次ML3.0 以上地震（分别是3 月24 日吴忠ML3.0 地震，4 月5 日灵武ML3.8 地震，5 月8 日灵武ML3.2 地震，6 月12 日青铜峡ML4.0 和ML3.4 地震），2021 年7 月11 日出现小震群活跃现象，11 月18 日发生ML4.5 地震，震情形势严峻。鉴于此，有必要对吴忠、灵武地区开展波速比和地震危险性研究工作。

为此，笔者利用宁夏地区数字化地震观测台网记录的中小地震波形资料，采用单台多震和达法，分析1980—2008 年吴忠—灵武地区MS≥5.0 地震前波速比变化，以期为当前震情趋势研判提供参考。

1 计算方法与约束条件

1.1 单台多震和达法

单台多震和达法由日本地震学家和达清夫1928 年提出。在假定震源区到地表介质为理想均匀弹性的条件下，纵波速度vP和横波速度vS与介质泊松比σ、弹性模量E、介质密度ρ之间的关系（李善邦，1981）为

由式（3）可见，vP/vS为介质泊松比σ的函数，主要反映地壳中上层介质泊松比的变化。P 波和S 波的速度则与介质的泊松比σ、弹性模量E 和介质密度ρ密切相关。根据P 波走时与P 波、S 波到时差的线性关系，可得vP/vS波速比值计算公式，即

线性相关系数R为

计算误差γ为

式中，tPi为P 波走时，tSi为S 波走时，n为每次地震到时数据个数，冯德益，1981）。该方法示意图见图1。

图1 单台多震和达法示意图（据李艳娥等，2014）Fig.1 Schematic diagram of the Wadati method of single station multi earthquakes(according to Li et al,2014)

1.2 约束条件

资料质量对利用和达法求多台波速比的影响众所周知。单台波速比异常与地震台站方位和孕震区范围明显相关，且接收穿透孕震区较好射线的地震台短期异常较明显（李 丽 等，2016）。其可能影响因素为：①每组参与计算的地震空间位置可能不稳定；②单次地震发生时刻误差不同。因此，在地震时空分布相对均匀且稳定、每次参与计算的数据量基本相同情况下，其计算结果的可靠性和稳定性较高。

在研究过程中，由于不同研究者在计算波速比时选取的最大震中距、相关系统数、计算误差及最小完整性震级不同，为了使得计算结果更加合理，文中统计分析多个文献中关于参数的设定及适用范围（表1），并据前人经验，取震中距200 km 范围内的记录资料进行分析，以确保震相数据的科学性与合理性。

表1 不同文献中波速比计算的参数设置及适用范围Table1 Parameter settings and application scopes for the calculation of wave velocity ratio used by different researchers

2 资料选取

吴忠—灵武地区位于银川盆地南部，区内主要发育黄河断裂灵武段、银川隐伏断裂南段、新华桥断裂、牛首山断裂等（图2），1970 年以来中小地震频发，据宁夏地震台网地震目录，共发生ML≥0 地震1 213 次，其中：ML0—0.9 地震71 次，ML1.0—1.9 地震277 次，ML2.0—2.9 地震671 次，ML3.0—3.9 地震166 次，ML4.0—4.9 地震23 次，MS≥5.0 地震5 次。部分地震台站及1980 年以来4 级以上地震分布见图2。

利用宁夏地震台网1980 年1 月—2008 年12 月记录的中小地震数据，采用单台多震和达法，研究该区中强地震前波速比变化特征。考虑到ML＜2.0 的小地震一般背景噪声干扰大、信号弱，读取的P 波和S 波到时误差大，故震级下限取为ML2.0。1980—2008年研究区共发生4 次MS≥5.0 地震，其中3 次地震发生在1987—1988 年，时间间隔小，故作为一次地震丛事件来进行回溯分析。

由于天然地震发生位置的不稳定性和成丛性，会对单台多震和达法计算结果产生较大影响，考虑到地震分布的不稳定性，以及和达法自身固有的局限性，本研究在tS-P确定条件下，取tS-P≤20，计算样本组地震数取50，误差估计γ≤0.05，其中考虑了孕震体介质非均匀性也可导致波速比和达曲线的误差增大和相关性降低的因素。误差分析表明，计算得到的波速比是可靠的。

3 波速比时序变化分析

以1984 年11 月23 日灵武MS5.2 地震和1987—1988 年MS≥5.0 地震丛事件为例，利用发震前后10 年内台站周边200 km 范围内多次地震波形到时数据，回溯分析2 次震例发生前震中周边200 km 范围内地震台站记录波速比变化。扫描结果使用不同颜色的三角形表示：红色三角形示意震前存在波速比异常；橘色三角形示意震前波速比虽有下降，但未降至均值线以下；绿色三角形表示震前波速比无异常变化；灰色三角形表示研究时段该台站无震相数据资料。

3.1 1984 年11 月23 日灵武MS 5.2 地震

计 算1984 年11 月23 日 灵 武MS5.2 地震前波速比变化，结果显示：灵武台（LWT）、中卫台（ZW1）和寺口子台（SKZ）在此次地震发生前出现波速比异常现象，3 个台站呈三角形展布，与震中位置分别相距约13 km、120 km、200 km，震中位于三角形边缘（图3）。

图3 灵武MS 5.2 地震前波速比异常台站分布Fig.3 The distribution of stations with wave velocity ratio anomalies before Lingwu MS 5.2 earthquake

选取此次灵武MS5.2 地震发生前后10 年内，3 个台站周边中小地震（图4）长周期和短周期记录，分析波速比vP/vS随时间的变化，结果见图5。为了消除单点突跳，对波速比进行5 点平均，图5 中同时给出5 点平均1 点滑动平均线，3 条直线中实线为均值，上下2 条虚线为1.5 倍均方差。

图4 LWT、ZW1、SKZ 台站参与波速比计算的地震空间分布Fig.4 Spatial distribution maps of earthquakes that are used in the calculations at LWT,ZW1 and SKZ stations

图5 1984 年11 月23 日灵武MS 5.2 地震前后波速比异常时序Fig.5 Time sequence diagram of wave velocity ratio anomaly before and after Lingwu MS 5.2 earthquake on November 23,1984

由图5 可见：寺口子台（SKZ）波速比1982 年4 月开始下降，1984 年1 月达最低点，于恢复过程中发生地震，即震前21 个月出现波速比下降，震前10 个月达波速比值最低点，后缓慢恢复；中卫台（ZW1）波速比1982 年3 月开始下降，1984 年1 月达最低点，后出现转折，并于恢复过程中发生地震。即震前22 个月出现波速比下降，震前10 个月达波速比值最低点，后缓慢恢复；灵武台（LWT）波速比1982年10月开始下降，1983年7月达最低点，后缓慢恢复，并于恢复过程中发生地震。

3.2 1987—1988 年地震丛

1987—1988 年期间，吴忠—灵武地区共发生8 次4.3 级以上地震（5 级以上地震3 次），分别是1987 年8 月10 日5.4 级地震、1987 年8 月10 日4.3 级地震、1987 年8 月10 日2 次4.4级地震、1988 年1 月4 日5.4 级地震、1988 年1 月8 日4.3 级地震、1988 年1 月10 日5.0级地震和1988 年4 月12 日4.3 级地震（表2）。由于发震间隔较小，作为一次地震丛事件进行波速比回溯分析。该区多震型地震多发，此次地震丛即为一次多震型地震序列。

表2 1978—1988 年MS 4.0 以上地震三要素统计Table2 Statistics of MS ≥4.0 earthquakes from 1978 to 1988

此次地震丛半径约7.7 km，空间分布较为集中。从地震台站波速比异常扫描结果可知，震前未出现明显波速比异常变化，其中银川台（YCH）、灵武台（LWT）、中卫台（ZW2）和同心台（TXT）震前波速比稍有下降，且围绕地震丛呈四边形均匀展布[图6(a)]。从波形结果来看，银川台（YCH）、灵武台（LWT）波速比下降幅度较大。

将发震时段定义为1987 年8 月10 日至1988 年1 月10 日，选取银川台（YCH）、灵武台（LWT）、中卫台（ZW2）和同心台（TXT）发震前后记录的中小地震[图6(b)]，分析该地震丛波速比异常时序变化，结果见图7。由图7 可见，4 个台站波速比分别于1985 年9 月、1986 年4 月、1985 年3 月、1985 年9 月出现下降，其中银川台（YCH）、中卫台（ZW2）和同心台（TXT）波速比在下降过程中出现反弹后继续下降。波速比计算结果显示，震区北部的银川地区和灵武一带波速比异常较为明显。

图6 1987—1988 年地震丛发生前波速比异常台站分布（a）及参与波速比计算的地震分布（b）Fig.6 Distribution of stations with wave velocity ratio anomalies before the occurrence of earthquake clusters in 1987-1988(a) and distribution of earthquakes involved in the wave velocity ratio calculation (b)

图7 1987—1988 年地震丛波速比异常时序Fig.7 Time sequence diagram of wave velocity ratio anomaly before earthquakes from 1987 to 1988

3.3 近期波速比扫描结果

选取2015 年1 月—2021 年12 月吴忠—灵武地区发生的中小地震，对该区波速比进行扫描，结果发现，仅海原台（HYU）出现波速比异常（图8）。分析该时段海原台波速比时序变化，结果见图9，可见：2019 年7 月开始，该台波速比急速下降，至2021 年1 月达低值点，后恢复转平，且转平过程中，台站周边200 km 范围内至今无地震发生。

图8 2015—2021 年研究区地震台波速比异常分布（a）及参与波速比计算的地震分布（b）Fig.8 Distribution of seismic stations with wave velocity ratio anomalies in the study area from 2015 to 2021 (a) and distribution of earthquakes involved in the wave velocity ratio calculation (b)

图9 近年来海原台波速比时序Fig.9 Time sequence diagram of wave velocity ratio in recent years at HYU Seismic Station

从波速比空间扫描结果来看，未来短期内，需要重点关注宁夏南部及邻区震情形势变化。从动力学角度来看，此次波速比值低值点后转折上升的变化，可能与区域应力调整有关。追溯历史地震发现，1990 年4 月26 日青海共和MS7.0 地震发生后一年内，宁夏固原地区受区域应力调整的影响，小震群活动较为活跃。与此类似，2021 年5 月22 日青海玛多地震发生后，6 月11 日至7 月1 日宁夏固原市原州区中河乡发生小震丛集活动，共记录可定位地震348 次，震群分布范围约610 km2，这是否意味着海原台波速比低值异常对应此次青海玛多地震仍有待研究。

4 结论

利用单台多震和达法，对吴忠—灵武地区波速比进行分析，结合已有研究成果，得到如下结论：

（1）文中仅对研究区扫描了2 组MS≥5.0 地震事件，其中1 组为3 次MS≥5.0 地震丛事件，但从波速比计算结果来看，2 组地震前波速比异常具有共性特点，即震中周围分布的多个台站出现不同程度的波速比值下降变化。震中位置分布于波速比下降或异常台站多边形连线边缘或内部。当然，由于该区域震例较少，其结论不具有普适性，有待积累资料予以论证。

（2）吴忠—灵武地区MS≥5.0 地震前波速比存在正常→降低→上升（恢复）→发震的异常变化过程，与国内诸多学者所得结论一致（蔡杏辉，2013；翁钊强，2015；李艳娥，2016）。该研究结果对于判断和圈定未来强震孕育的可能空间范围具有一定指示意义。然而，该方法依然受到地震定位精度及震相到时判读精度的影响，波速及波速比变化的可靠性仍需波形相关等技术的进一步深入研究予以确认。

（3）根据对研究区近期中小地震扫描结果来看，海原台（HYU）波速比异常较显著，需密切关注宁夏南部及邻区地球物理异常变化。

文中分析使用了中国地震局地球物理研究所李艳娥提供的计算程序，且审稿专家提出诸多中肯意见，在此一并表示衷心感谢。