1. 引言
银川以北至宁蒙交界地区是宁夏社会和经济建设的重要区域，银川以北至宁蒙交界地区被列为2006~2020年全国地震重点监视防御区。该区域自
1991 年 9 月 14 日
内蒙古阿拉善左旗5.0级地震后平静了23年7个月，于
2015 年 4 月 15 日
发生了内蒙古阿拉善左旗5.8级地震和
2017 年 6 月 3 日
内蒙古阿拉善左旗5.0级地震。根据以往震例，该区域未来几年内5.0级以上地震存在连续发生的可能，甚至不排除发生6~7级地震的可能。
2015 年 4 月 15 日
内蒙古阿拉善左旗5.8级地震发生后，需要对该区域未来中强以上地震趋势进一步进行多角度分析。
2. 宁夏北部及邻区强震构造危险性分析
宁夏地处青藏地块、华北地块和阿拉善地块的交接部位，受多种动力学因素的制约，分布有多条控制现今强震活动的活动断层，以牛首山为界，南部属于祁连山断裂带的部分，以北则属于华北地块的鄂尔多斯地台西缘断裂带的部分。由于宁夏北部长期受青藏块体北东向的推挤作用和阿拉善地块之间强烈的相互作用影响，造成该区地质构造十分复杂，地震活动十分强烈。上溯1000年来，宁夏北部及邻区(37.5˚~40.7˚N, 103.5˚~107.7˚E)地区发生过8次6级以上地震，包括1次8级大地震。与历史地震相对照，该区现今地震活动具有继承性，无论频度和强度仍然显示出较高的活动水平，自1950年以来共发生Ms ≥ 5.0地震27次。
2015 年 4 月 15 日
内蒙古阿左旗5.8级地震发生在宁夏以北华北地块的鄂尔多斯西北缘，而
2017 年 6 月 3 日
内蒙古阿左旗5.0级地震发生在祁连山地震带。
宁夏北部及邻区7级以上地震都与几条主要活动断裂带息息相关。如贺兰山东麓断裂带和香山–天景山断裂带等均发生过7级以上地震，大多数6级地震也是沿几条主要的活动断裂带分布。另外，根据近几十年野外地质调查发现，宁夏境内的几条主要活动断裂带上揭露出多期7级以上古地震事件，进一步说明活动断裂与强地震的依存关系。第四纪盆地是易发生强地震的场所。银川盆地和吉兰泰盆地是本区两个活动强烈的拉张型断陷盆地。
银川盆地历史上强震非常活跃，曾发生过1739年银川–平罗8级大震。位于盆地南端的吴忠、灵武地区，自1921年以来发生过8次5级地震。在中卫和中宁的压陷型盆地，历史上曾分别发生过1709年中卫南7 1/2级和1561年中宁7 1/4级。
宁夏城市活动断层探测及地震危险性评价工作走在了全国前列，对研究区主要活动断裂第四纪以来的活动性研究成果为地震危险地点的判定提供了重要依据。根据银川市活动断层探测与地震危险性评价结果，未来200年(2006年以后)，贺兰山东麓断裂北段发生特征地震的危险程度很低，贺兰山东麓断裂南段复发特征地震的可能性为中等，银川隐伏断裂北段发生特征地震的可能性低，黄河断裂南段特征地震的复发危险程度高，存在发生7级地震的可能。
根据M7专项工作组(M7专项工作组，2012) [1] 对南北地震带北段中长期大地震危险性研究成果，即南北地震带北段历史与史前地震破裂区展布及大地震空区判定结果，未来可能发生强震的空段有6个，宁夏北部及邻近地区强震空段为磴口–五原地震空段、同心–灵武地震空段。
依据地震行业科研专项《2016~2025年中国大陆地震危险区预测》研究成果，未来7年宁夏北部及邻区存在3个强震危险区，分别为鄂尔多斯西北缘、香山天景山断裂东南段–黄河灵武断裂南段和天桥沟黄羊川断裂–老虎山断裂。
3. 宁夏北部及邻区1950年以来5级以上地震活动特征
从图1可以看出，宁夏北部及邻区1950年以来5级以上地震可以分为4个活跃时段，分别为：1952~1962年、1971~1976年、1984~1991年和2015~2017年，分别持续10、5、7和2年，最后一个活跃时段时间最短，根据前3个活跃时段的持续时间，第4个活跃时段有可能持续几年，但不应超过第1活跃时段5级以上余震持续时间(8年)；3个平静时段分别为：1962~1971、1976~1984和1991~2017年，持续时间分别为9、8和24年，最后1个平静时段最长。
从图2(a)可以看出，1952~1954年宁夏北部及邻区的5级以上地震主要分布在主震——
1954 年 7 月 31 日
7.0级地震余震区内(以震中为中心，南北距离
80 km
，东西
40 km
的矩形内)。1955~1957年宁夏北部及邻区没有5级以上地震发生(图1)，1958~1962年宁夏北部及邻区的5级以上地震依然主要分布在
1954 年 7 月 31 日
7.0级地震余震区内，同时在外围出现了4个响应区，宁夏北部及邻区有3个响应区(图2(b))。
从图3(a)和图3(b)可以看出，第2、3活跃时段宁夏吴忠–灵武响应区都发生了5级以上地震，内蒙古巴音木仁地区第2活跃时段也发生了5级以上地震，第3活跃时段在内蒙古阿拉善左旗发生了3次5级地震，位于1954年民勤7.0级地震响应区之间；1991年之后宁夏北部及邻区5级以上地震活动平静26年之后，
2015 年 4 月 15 日
内蒙古阿拉善左旗发生5.8级地震，之后在1954民勤7.0级地震余震区5级地震平静56年后2017年6月3日发生了5.0级地震，第4活跃时段的2次5级以上地震依然发生在1954年民勤7.0级地震的响应区内(图4)，因此可以判断第4活跃时段5级以上地震活动地点有可能维持目前的格局。
综合以上分析，可以推断未来几年宁夏北部及邻区5级以上地震活动有可能继续活跃，但强度不超过6级，地点有可能继续发生在1954年民勤7.0级地震的余震区和响应区内或者在内蒙古阿拉善左旗地区。
Figure 1. M-t map of earthquakes with MS ≥ 5.0 in northern Ningxia and adjacent areas from 1950 to 2018
图1. 1950~2018年宁夏北部及邻区5级以上地震M-t图
Figure 2. Distribution of earthquakes with MS ≥ 5.0 in northern Ningxia and adjacent areas from 1952 to 1962
图2. 1952~1962年宁夏北部及邻区5级以上地震分布图
Figure 3. Distribution of earthquakes with MS ≥ 5.0 in northern Ningxia and adjacent areas from1971 to 1991
图3. 1971~1991年宁夏北部及邻区5级以上地震分布图
Figure 4. Distribution of earthquakes with MS ≥ 5.0 in northern Ningxia and adjacent areas from 2015 to 2018
图4. 2015~2018年宁夏北部及邻区5级以上地震分布图
4. 近期宁夏北部及邻区存在的广义前震特征
通过分析2015年5月以来宁夏北部及邻区ML ≥ 3.0地震活动情况，发现2015年5月以来甘肃民勤至阿拉善左旗一带形成了北东向的较大范围ML ≥ 3.0地震围空区，长轴约240 km (图5)。该空区的存在表明未来在空区内部及边缘有发生6级左右地震的可能。
Figure 5. Distribution of earthquakes with ML ≥ 3.0 in northern Ningxia and adjacent areas since 2009
图5. 2009年以来宁夏北部及邻区的ML ≥ 3.0地震分布
5. 近期宁夏北部及邻区存在的地震活动性参数异常
地震活动性参数是指采用地震的时空强等信息根据统计模型回归得到的参数值，其与地质构造、岩石力学性质、地下介质物性以及震源机制等有密切的关系。研究地震活动性参数的时空变化对于了解研究区的构造环境、应力环境等有重要的科学意义。以下给出了一些重要的地震活动性参数异常分析结果。
1) 沿主要断裂带的b值空间扫描图像
特定岩体内构造应力的大小与破裂产生地震的震级—频度关系为lgN(M) = a − bM，式中N为震级M ≥ Mmin (最小完整性震级)的地震累积频度，a、b为经验常数。由于应力的高低与震级–频度关系中的b值成反比(Scholz, 1968; Urbancic et al., 1992) [2] [3] ，因此可以根据低b值的空间分布来揭示和推断活动断裂带不同段落现今应力积累的相对水平，描绘出活动断裂带上的高应力积累区或闭锁段落(Wiemer, 1997; Wyss et al., 2000) [4] [5] 。
利用全国正式目录，选取1970年以来宁夏及邻近地区(34˚~41.3˚N, 103.5˚~109˚E)记录的弱震资料，将研究区以0.1˚ × 0.1˚划分网格，固定每个单元的搜索地震数，并计算相应的最小完整性震级，然后利用震级—频度关系基于最大似然估计法计算每个单元格内的b值，得到宁夏及邻区沿主要断裂带的b值空间分布图像(图6)。结果显示：鄂尔多斯西北缘存在高应力集中区。
2) 地震加速矩释放时空扫描图像
蒋长胜等(2009) [6] 用“破裂时间分析”方程中的幂指数m作为描述震前加速矩释放程度的参量，结
Figure 6. b value space scan of Ningxia and adjacent areas (since 1970)
图6. 宁夏及邻区b值空间扫描图(1970年以来)
果表明地震前可观测到m值逐渐减小，即加速特征逐渐明显的过程，并且出现AMR现象的空间区域与震中位置似有较好的对应。
利用宁夏及邻区近5年来ML2.0以上地震资料，取空间扫描半径为60 km，得到该区域地震矩加速释放模型参数m值空间分布图像(图7)。结果显示，内蒙古阿拉善右旗和内蒙古磴口–杭锦后旗一带存在地震矩加速释放低m值异常区。
3) A值平面扫描图像
地震空间集中度C值反应了地震在空间上的集中程度，地震危险度D值反应了地震在时间上的丛集程度，地震强度因子Mf值反应了地震活动增强的程度，它们分别反应了地震在时、空、强方面的震兆。为综合考虑地震活动的时、空、强特征，王炜教授(1986) [7] 定义地震活动因子A值为C、D、Mf值的线性组合。A值异常场区与未来地震有较好的相关性，统计1971年以来的预报效能，R大于0.207，通过R检验。选取宁夏及邻区ML2.5以上地震，以时间尺度1年、步长0.5˚和窗长1˚进行A值空间滑动扫描，取大于0.40为异常值。从图8可以看出异常区为内蒙古乌海–磴口和宁夏银川周围一带。
4) 区域Pg波速度横向变化
利用宁夏测震台网记录的1970年以来的Pg波到时数据，得到了宁夏及邻区Pg波速度横向变化图像(图9)。其中，蓝色代表高速区，红色代表低速区，白色圆圈为1970年以来发生于宁夏及邻区的5.0级以上地震。速度成像结果显示：除黄河断裂南段–灵武断裂呈现为高速区外，其它断裂带上则以低速为主；1970年以来5级以上地震多发生于高、低速的转换地带。
Figure 7. Feature m-value spatial distribution of acceleration release of seismic moment in Ningxia and adjacent regions
图7. 宁夏及邻区地震矩加速释放特征m值空间分布
Figure 8. Abnormal distribution of A value (2017.10~2018.10)
图8. A值异常分布(2017.10~2018.10)
Figure 9. Horizontal variation diagram of Pg wave speed in Ningxia and adjacent regions
图9. 宁夏及邻区Pg波速度横向变化图
基于地震易发生于高、低速区的过渡带和坚固体孕震模式(梅世蓉，1995) [8] ——高速体(坚固体)的存在是高应力集中的重要条件，宁夏吴忠–灵武以北的高速区及其周围很可能是未来地震孕育区域。
5) 视应力空间分布
视应力常被用来研究震源区的应力水平，以及了解地震发生过程中的构造应力释放水平，即地震视应力越高，震源区的应力水平就越高。研究(王琼等，2005；秦嘉政等，2006) [9] [10] 发现，视应力与震级之间存在一定的依赖关系。因此为尽量消除视应力与震级之间的相关性，特引入差视应力概念，即指实际计算得到的视应力与理论视应力的差值。理论视应力可根据视应力与震级的线性关系得到。
选取2009年以来宁夏测震台网记录的ML2.5以上地震，依次截取波形，扣除仪器响应，采用Brune圆盘模型，计算每个中小地震的视应力值，拟合视应力与震级标度关系，视应力值扣除震级影响后即为差视应力值。最后以10 km为扫描半径、0.1˚为扫描步长进行空间扫描，得到差视应力空间分布图(图10)。结果显示，在宁夏石嘴山–银川、青铜峡–吴忠为未来可能发生中等地震的区域。
6) 地震分形维度D值
地震活动性的分形特性指的是地震在时间和空间分布的随机自相似结构，其可以使用分形维度D值来表示，根据D值的大小可确定地震在空间是随机分布还是聚集分布的(Ogata, 1988, Hirata, 1989) [11] [12] 。分形维度也可以用来定量描述某一断层系统地震活动的差异程度以及地震地表破裂的复杂程度，这时D主要受应力场分布、地质、力学及构造差异等因素的控制。分形维度值为判断一个地区是否处于稳定状态提供了十分重要的信息。分形维度值的任何变化都意味着地区系统状态的改变。因此研究余震区的地震分形维度对于认识地震地表破裂的复杂性以及判断余震活动趋势具有重要的意义。
Figure 10. Abnormal distribution of high apparent stress
图10. 视应力高值异常分布
Figure 11. D value space scan of Ningxia and adjacent areas (since 1970)
图11. 宁夏及邻区D值空间扫描图(1970年以来)
通过计算1970年以来宁夏及邻区(103.5˚~107.7˚E, 35.0˚~40.7˚E) ML ≥ 3.0地震的D值(图11)，从图中可以看出，1970年以来宁夏及邻区的D值偏高，各个区域的特征不同：宁夏中北部地区、甘肃景泰地区和内蒙古阿拉善左旗5级以上地震发生在D值1.9左右的区域；宁夏南部和内蒙古乌海一带的5级以上地震的5级以上地震发生在D值大于2.1的区域。以上特征对今后判定危险区提供了一定的参考。
6. 结论与讨论
通过前面的背景分析，得到了以下结论：可以推断未来几年宁夏北部及邻区5级以上地震活动有可能继续活跃，但强度不超过6级，地点有可能继续发生在1954年民勤7.0级地震的余震区和响应区内或者在内蒙古阿拉善左旗地区。2015年5月以来甘肃民勤至阿拉善左旗一带存在北东向的ML ≥ 3.0地震围空区(长轴约240 km，图7)，该空区的存在表明未来在空区内部及边缘有发生6级左右地震的可能。该空区内部和边缘包含了1954年民勤7.0级地震的余震区和2个响应区，因此该区域未来发生6级以上地震的可能性较大。再结合宁夏北部及邻区近几年的地震活动性参数异常，得到以下结论：
1) 内蒙古乌海至磴口一带为未来宁夏及邻区发生6级以上地震的危险区域，其异常指标最多，有ML ≥ 3.0地震空区、地震矩加速释放低m值异常区、A值平面扫描异常区、Pg波速度横向变化和分形维度D值异常，采用的异常分析绝大多数出现了异常，仅b值和视应力的空间分布图像未出现异常，因此需重点跟踪该区域的震情发展。
2)
1954 年7 月31 日
内蒙古腾格里沙漠北7.0级地震的余震区未来几年也有发生6级以上地震的危险，其异常指标有ML ≥ 3.0地震空区和分形维度D值异常，因此对该区域的震情发展需进一步跟踪。
基金项目
宁夏回族自治区重点研发计划(科技惠民)“内蒙古阿拉善左旗5.8级地震后宁夏北部及邻区地震危险性分析”(项目编号：2016KJHM135)资助。