柴达木—阿尔金地震带b 值统计计算1

谢江丽 阿布都瓦里斯·阿布都瓦衣提 黄帅堂

（新疆维吾尔自治区地震局, 乌鲁木齐 830011）

引言

b值是震级和频度关系式的比例系数，由实际地震资料统计得到，其与地震资料的可靠性、地震取样的时空分布特征、地震样本的起始震级等因素有关。根据研究区域特点，选择合适的分析方法确定b值，对于开展地震危险性分析具有重要意义。

1 柴达木—阿尔金地震带地震资料统计分析

柴达木-阿尔金地震带在行政区划图上北部涉及新疆维吾尔自治区和甘肃省，南部进入青海省和西藏自治区，在地震构造上由阿尔金山脉和柴达木盆地组成（姚远等，2014）。柴达木-阿尔金地震带相对于相邻地震带（西昆仑-帕米尔地震带、鲜水河-滇东地震带）地震活动强度小，发生频度低，但近10 年发生了2 次7 级地震，值得关注。本文通过对柴达木-阿尔金地震带历史地震资料的收集整理，结合柴达木-阿尔金地震带地震活动时间、空间分布特点，使用合适的统计方法对柴达木-阿尔金地震带b值进行统计计算，其中，现代小震震级ML不再转换为MS（汪素云等，2009），直接表示为M。

1.1 地震资料收集

本文采用的地震资料来源如下：《中国历史强震目录（公元前23 世纪-公元1911 年）》《中国近代地震目录（公元1912-1990 年，MS≥4.7）》《新疆维吾尔自治区地震目录（1970-1979）》《新疆维吾尔自治区地震目录（1980-1984）》《中国地震年报》《中国地震详目》《中国地震台网观测报告》《中国震例》。

柴达木-阿尔金地震带包括阿尔金山脉与柴达木盆地，东北起点为新疆与甘肃交界处，西南至琼木孜塔格一带，总体走向北西西（图1），地震带东、西段地震活动强于中部地区。大震主要发生在地震带东、西段，中部地区基本无大震发生，且中强地震活动性较低。有历史记录以来，柴达木-阿尔金地震带共记录到5 次7 级以上地震，地震带西段发生了4 次，包括2008、2014 年2 次于田7.3 级地震；地震带东段仅记录到1 次7 级地震。

图1 柴达木-阿尔金地震带范围示意图Fig. 1 Distribution diagram of Qaidam-Altun Seismic belt

柴达木-阿尔金地震带最早1 次地震记载为公元1832 年8 月昌马5½级地震，最大震级地震为2008 年3 月21 日和2014 年2 月12 日于田7.3 级地震。由于柴达木-阿尔金地震带人烟稀少，该地区地震记载时间较短，该地震带1920 年以前地震资料严重缺失，仅有4 次地震记录。1920 年以后5.0 级以上地震记录基本完整（王海涛等，2006；单新建等，1996）。至2019 年12 月，该地震带共记录到M≥4.0 地震885 次，其中，7.0～7.9 级地震5 次，6.0～6.9 级地震23 次，5.0～5.9 级地震167 次，4.0～4.9 级地震690 次。

由于不同的历史时期地震监测能力差异较大，为保证资料连续、可靠、充分，根据任雪梅等（2011）的研究，将柴达木-阿尔金地震带内1920 年以后M≥5.0 历史地震和1970 年以后区域台网记录的M≥4.0 地震作为本次统计基本资料。柴达木-阿尔金地震带M≥6.0 地震目录如表1 所示。

表1 柴达木—阿尔金地震带1920—2019 年M≥6.0 地震目录Table 1 Earthquake catalog M≥6.0 of Qaidam—Altun Seismic belt from 1920 to 2019

1.2 余震删除

由于余震对地震活动性统计特征的影响较大，因此对地震目录进行统计分析前，首先须删除余震。地震学家已发展了多种删除余震的方法（Keilis-Borok 等，1980），本文依据断层长度R与震级M的经验公式（M=3.3+2.1logR）确定删除余震的空间半径r与主震震级M的关系为（郭增建等，1979）：

考虑到地震实际定位精度，取r=5 km 作为空间窗半径最小值。删除余震时间窗，采用余震时间窗（陈凌等，1998），不同震级余震时间窗如表2 所示。

表2 余震时间窗Table 2 Aftershock time window

以2008 年3 月21 日于田7.3 级地震为例，本次统计基本资料取2008 年3 月21 日至2008 年12 月31日时间段记录的地震目录，剔除不在此区域范围的地震，绘制此次地震空间扫描图（图2），根据余震空间半径r与主震震级M的关系式得出本次地震余震空间半径为85 km，由图2 可知，余震基本在此范围内，证明余震分布符合此关系式，根据余震时间窗删除余震，删除前、后地震数目对比如表3 所示。

表3 2008 年于田7.3 级地震余震删除前、后地震数Table 3 Earthquakes before and after the deletion of Yutian M7.3 earthquake in 2008

图2 2008 年于田7.3 级地震空间扫描图Fig. 2 Spatial scan of Yutian M7.3 earthquake in 2008

选取柴达木-阿尔金地震带1920-2019 年M≥4.0 地震目录为基础资料，根据余震空间窗和时间窗进行扫描，手动删除余震，进行删除前、后地震数目对比（表4），本研究选取删除余震后地震目录作为原始资料。

表4 地震带地震目录余震删除前、后地震数Table 4 Earthquakes before and after the deletion of aftershock

1.3 地震目录完整性及地震带活动趋势分析

地震活动性研究中，地震目录的不完整会使地震活动性参数估计形成偏差（任雪梅等，2011），研究区域地震目录可靠性尤其重要。已有学者对新疆及相邻地区地震目录进行多次完整性分析（王海涛等，2006；徐伟进等，2014），根据本次研究区域特点，本文采用徐伟进等（2014）对新疆不同起始震级目录完整起始时间，即M为3.0、4.0、5.0、6.0、7.0 级时，地震完整起始时间分别为1980、1975、1945、1935、1840 年。

利用柴达木-阿尔金地震带地震目录绘制地震时序图，如图3 所示。由图3 可知，1980 年后地震分布较密集，这可能与地震监测手段逐步完善、地震记录不遗漏有关，说明历史地震遗漏较多。另外，地震在时间上分布不均匀，存在成丛性特点，地震活动性有强弱之分，有平静期也有活跃期。从地震发生趋势上看，本研究认为现阶段该地震带处于地震活跃期，故本次计算b值考虑未来一段时间该地震带处于地震活跃期。

图3 1920-2019 年M≥4.0 地震时序图Fig. 3 Sequence diagram of M≥4.0 earthquakes from 1920 to 2019

2 柴达木—阿尔金地震带地震活动性参数计算

在概率地震危险性方法（PSHA）中，以潜在震源区为基本单元建立地震活动性模型，每个潜在震源区具有独立的地震活动性参数，主要包括震级上限、b值和4.0 级以上地震的年平均发生率V4.0，合理确定活动性参数是地震危险性概率分析的基础工作（潘华等，2006，2013）。根据对国际上30 多个国家采用概率方法的地震区划工作的调研，目前较常见的做法是活动地区b值取0.8，稳定地区b值取0.95，这种方法存在的问题是未考虑地震活动的时空不均匀性。我国CPSHA 方法提出分层级的潜在震源区划分，同时兼顾了地震活动性参数统计与地震活动空间分布不均匀性表达的需要《工程场地地震安全性评价》（GB 17741-2005）（中华人民共和国国家质量监督检验检疫总局等，2006）。

2.1 地震带b 值

在概率地震危险性分析方法中，地震带活动性参数包括地震带震级上限Muz、起算震级M0、震级频度关系式中的b值（《工程场地地震安全性评价》（GB 17741-2005）；陈鲲，2014；《中国地震动参数区划图》（GB 18306-2015）（中华人民共和国国家质量监督检验检疫总局等，2016）），其统计表达式为：

b值作为地震震级频度分布特征和地震活动水平的重要参数。我国目前历史地震资料分析中地震带b值为0.5～0.6，不少学者利用现代仪器记录地震资料的b值为0.8～1.0（刘杰等，1997；陈培善等，2003；王海涛等，2006）。而在目前区划中，b值统计中存在许多影响因素，如b值统计的时空范围选择、震级上限和震级不确定性等。计算方法选择不同结果不同（郑确等，2018），常用方法为最大似然法和最小二乘法，其中最大似然法对样本数量要求较高（吴果等，2019），最小二乘法合理估计b值已相当成熟，且计算结果较合理（杨春柳等，2010；常金龙，2017）。本次研究区域历史记录时间较短，地震记录较少，综合实际资料情况，采用最小二乘法拟合b值较准确。

根据徐伟进等（2014）对新疆不同起始震级目录完整起始时间，本研究选取1975 年后4 级地震目录、1945 年后5 级地震目录（震级档选为0.5），对应的b值分别为0.739 3（图4（a））和0.771 1（图4（b））。本研究确定b值的原则是不低估已认识到的地震带地震危险性，同时合理考虑未来地震危险性，b值取值偏高往往会造成高估中强地震发震能力、低估大地震发震能力，b值取值偏低可能造成高估大地震发震能力、低估中强地震发震能力（《中国地震动参数区划图》（GB 18306-2015）），故本次地震带b值取平均值0.76。

图4 b 值拟合曲线Fig. 4 Fitted curve of b-value

柴达木-阿尔金地震带有地震记录以来，历史最大地震震级为7.3 级，2014 年于田7.3 级地震作为此地震带特征地震，对高震级段的地震发生率和地震带活动水平均有重大影响。此地震使地震带7.0～7.5 级震级段年发生率由每年0.04 上升至0.05，b值由0.720 5 下降为0.688 6（图5），由此可见，大地震的发生会降低地震带b值。故2014 年于田7.3 级地震的发生对该地震带活动强度和大震重现期等地震活动特征和活动水平的认识具有重要意义。

图5 2014 年于田7.3 级地震发生前、后b 值拟合曲线Fig. 5 Fitted curve of b-value before and after the Yutian M7.3 earthquake in 2014

2.2 不同地震构造区b 值

在现今地球动力学环境下，将新构造活动特点、构造应力场、深部地球物理场等构造环境与发震构造模型具有一致性或相似性特点的区域称为地震构造区。每个地震带内均由不同地震构造区组成，柴达木-阿尔金地震带包括青海湖地震构造区、柴达木地震构造区和阿尔金地震构造区（《工程场地地震安全性评价》（GB 17741-2005）），如图6 所示，根据不同构造区范围提取地震目录（表5），按照震级完整时段计算各自b值，如表6、图7 所示。

图6 柴达木-阿尔金地震带地震构造区分布Fig. 6 Distribution map of seismotectonic area in Qaidam-Alkin seismic belt

图7 不同地震构造区b 值拟合曲线Fig. 7 Curves of b-value fitting in different seismic tectonic regions

表5 不同构造区地震目录Table 5 Earthquake catalog of different structural areas

表6 不同构造区不同起算震级计算b 值Table 6 Calculation of b-values with different starting magnitudes in different structural areas

3 结论与讨论

（1）本文利用现有地震资料，对柴达木-阿尔金地震带b值进行统计研究，给出了该地震带b值统计评估结果。本研究确定b值的原则是不低估已认识到的地震带地震危险性，同时合理考虑其未来地震危险性，故本文地震带b值取平均值0.76。由于该地震带地震活动频度低、强度大，且地震记录时间短，故地震带上特征地震的发生频度及中、小地震活动频度对该地震带地震活动性参数的确定具有显著影响，2014 年于田7.3 级地震及近10 年中、小地震的发生，使b值取值减小，2014 年于田7.3 级地震的发生对该地震带活动强度和大震重现期等地震活动特征和活动水平的认识具有重要意义。

（2）本文根据柴达木-阿尔金地震不同构造区提取地震目录，按照震级完整时段计算各自b值，所得b值范围为0.536 5～0.801 9，最大差值为0.272 6。

（3）目前工程地震中使用的b值地震资料至2010 年，距今有10 多年的地震资料缺失，而新一代区划图尚未编制。本文统计计算b值补充了近年来地震资料，在下一代区划图编制实施前，可为科学建立重大建设工程地震危险性计算模型提供基础依据，对该地震带重大建设工程地震安全性评价与地震小区划工程均有重要现实意义。