谢江丽 阿布都瓦里斯?阿布都瓦衣提 张文秀

摘要：利用1920—2020年地震资料，对柴达木—阿尔金地震带采用多方案（不同时间段、震级段）统计计算并综合分析，确定b值为0.78，较《中国地震动参数区划图》中确定的b值0.84有所降低，对地震带重大建设工程地震安全性评价与地震小区划的工程设防有重要意义。本文确定b值作为《中国地震动参数区划图》的补充，可为科学建立重大建设工程的地震危险性计算模型提供基础依据。

关键词：柴达木—阿尔金地震带；b值；多方案；地震动参数

doi：10.16256/j.issn.1001-8956.2023.03.010

b值的统计结果受时空范围、震级上限和震级的不确定性等多种因素影响^［1，2］，选择不同的计算方法结果也不同^［3］。常用的方法为最大似然法和最小二乘法，其中最大似然法对样本数量要求比较高^［4］，利用最小二乘法合理估计b值已经相当成熟，且计算结果也比较合理^［5］。本文研究区域历史记录时间较短，地震记录较少，综合实际资料情况，采用最小二乘法拟合b值比较准确。利用柴达木—阿尔金地震带1920—2020年地震资料，根据不同的统计时间段和震级段进行计算，并用地震年发生率约束计算结果，综合分析确定b值^［6］。此方法充分考虑了地震资料不完备性和认识不确定性，同时利用历史地震资料和现代仪器记录地震资料，充分发掘有价值的信息，克服样本资料不完备可能导致的认识偏差，更好的体现柴达木—阿尔金地震带的地震活动水平。通过补充近年来地震资料重新统计计算b值，在下一代区划图编制实施前，可为科学建立重大建设工程地震危险性计算模型提供基础依据，对该地震带重大建设工程地震安全性评价与地震小区划工程均有重要意义。

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

柴达木—阿尔金地震带西南从琼木孜塔格一带起，东北至甘肃与新疆交界，由柴达木盆地和阿尔金两个地震构造区组成^［7］，总体走向NWW。该地震带中部地区地震活动水平较低，大震主要集中在地震带的东段和西段。

通过对柴达木—阿尔金地震带历史地震资料收集整理，结合其地震活动时空分布特点^［8］，使用多方案的统计方法对柴达木—阿尔金地震帶的b值进行统计计算，汪素云等人对震级标度转换进行了专门的研究，认为M_L≈M_S^［9］，因此，本文在编制地震目录时不再进行面波震级与地方震级的转换^［10］，且统一表示成M。

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

柴达木—阿尔金地震带最早1次地震记录为公元1832年8月昌马5?级地震，历史上最大震级为7.3，分别为2008年3月21日和2014年2月12日于田M7.3地震。由于柴达木—阿尔金地震带历史地震资料较少、时间记载较短，资料缺失严重，1920年以前仅记录到4次地震。1920年以后M≥5.0地震记录基本完整^{［11、12］}。截至2020年12月，该地震带共发生901次M≥4.0地震，其中5次7级地震，24次6级地震，169次5级地震和703次4级地震。5次7级地震中4次发生在地震带西段，包括时间最近的2014、2008年于田的两次7.3级地震；1次7级地震发生在地震带东段，地震带西段强震发生率高于东段。柴达木—阿尔金地震带M≥6.0 地震震源参数见表1。

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

在研究地震活动性时，地震目录完整性尤其重要，其对地震活动性参数计算结果产生较大影响。地震监测能力在不同历史时期的差异性较大，不少学者多次对新疆及相邻地区地震目录进行完整性分析^［13］。为保证地震资料的可靠性、连续性和充分性^［14］，结合本次研究区域特点，使用统计基本资料为1920年以来柴达木—阿尔金地震带内M≥5.0历史地震和1970年以来区域台网记录M≥4.0地震。

利用柴达木—阿尔金地震带地震目录绘制地震时序图（图1），随着地震监测手段逐步完善，能记录到的地震最小完整性震级降低，地震分布越来越密集，间接说明历史地震遗漏较多。另外地震活动平静与活跃交替出现，存在成丛性特点，在时间上分布很不均匀。从地震发生趋势上看，本文认为该地震带正处于地震活跃期，在此基础上进行b值计算。

3 地震带余震删除

地震活动性统计特征与余震关系密切，余震对地震独立性有较大影响。假设震源之间具有相互独立性，余震被认为是和主震相关的事件，与主震之间不是独立关系，如果余震去除的不完整，地震的不独立性会使其不符合泊松分布。未删除余震会造成地震活动性统计分析不准确^［15］，故需先删除余震再对地震目录进行统计分析。目前删除余震的方法已有许多种。利用M（主震震级）的经验公式（M=3.3+2.1logR）与R（km）断层长度确定M（主震震级）与r（km）余震删除空间半径的关系为^［16］

logr=0.48M-1.57 .（1）

式中：r是余震删除空间半径，M是主震震级。根据地震实际定位精度，取r=5 km作为空间窗半径的最小值。采用陈凌等人研究的不同震级余震时间范围^［17］，确定4.0≤M＜4.5余震时间为42天，4.5≤M＜5.0为83天，5.0≤M＜5.5为155天，5.5≤M＜6.0为290天，6.0≤M＜6.5为510天，6.5≤M＜7.0为790天。

从本次统计基本资料选取2008年3月21日至2009年3月21日内记录的地震，剔除不在此区域范围内的地震，绘制于田7.3级地震余震分布图^［18］（图2），根据余震的空间半径r与主震震级M的关系式得出本次地震余震空间半径r=85 km，由图2可以看出，余震基本都在此范围内，证明余震分布符合此关系式。

本文选取了柴达木—阿尔金地震带1920—2020年12月M≥4.0地震目录为基础资料，余震空间窗和时间窗进行扫描，手动删除余震，同时进行删除前后地震数目对比（表2），本次研究选取的数据资料为删除余震后的地震目录。

4 地震带b值计算方法与结果

在概率地震危险性分析方法中，地震带活动性参数包括地震带震级上限M_uz、起算震级M₀、地震年发生率V₀及震级频度关系式中的b值^{［19，20］}，其统计表达式为

lgN=a-bM .（2）

式中：N为不同震级档地震发生的次数；M为地震震级；a、b为常数。

V₀是指在一定范围内（如地震带）平均每年发生震级等于或大于M₀（通常用M=4）的地震数，代表了统计区内的地震活动水平。V₀在本文中作为确定b值的一个约束参数。

为考虑地震资料的完备性，根据起算震级和地震完整性可信时段的不同选取多种统计方案，同时考虑每个起算震级的地震完整可信时段样本量是最充分的，也为避免忽略复发周期长的大地震对b值的影响，选取了1950—2020年M≥5.5、M≥5.0地震和1970—2020年M≥4.0、 M≥4.5等代表整个地震带可信时间段的统计样本。利用b值统计关系式和最小二乘法得出b值拟合曲线，图3为部分样本选取示意图，图4为部分b值拟合曲线图。

根据上述选取样本计算b值和其对应的地震年发生率，忽略因地震不完备性造成的b值过高或过低的样本，可信样本为33个（表3），同时对b值和地震年发生率结果进行散点绘图（图5）。

如图5（a）所示，b值结果处于0.47～0.81范围内，基本符合震级越高b值越低的规律。如图5（b）所示，不同震级的地震年发生率均在一定范围内变化，震级越高，其地震年发生率越低。75%～85%范围内的V₄所对应b值都在0.78附近，综合考虑b值结果和V₄，b值取0.78与之最接近V₄所对应b值都在0.78附近，故b值取0.78。该结果低于《中国地震动参数区划图》，可为判定该地区地震危险性提供参考。

5 结束语

利用现有地震资料，对柴达木—阿尔金地震带b值进行统计研究，采用多方案统计计算得到b值0.78，低于《中国地震动参数区划图》b值0.84，b值的降低，影响该地震带重大建设工程地震安全性评价与地震小区划工程设防参数的调整。针对地震记录时间短，历史地震资料中、小地震不全的地震带，本文所用多方案统计方法克服了样本资料不完备性可能導致的认识偏差，b值取值相对合理，历史地震资料情况类似的阿尔泰山地震带等可以用本方案计算b值。

参考文献：

［1］ 许洪泰，闵伟.利用地震矩和实时概率模型对中长期地震危险性进行评价——以东昆仑和阿尔金断裂为例［J］.地震研究，2007，30（2）：120-126.

［2］ 谢卓娟，吕悦军，兰景岩，等.b值和V₄的统计分析及其不确定性对地震危险性分析结果的影响研究［J］.地震研究，2013，36（1）：86-92.

［3］ 郑确，刘财，田有，等.地震活动性中震级-频度关系研究进展与再认识［J］.地球物理学进展，2018，33（5）：1 879-1 889.

［4］ 吴果，周庆，冉洪流.震级-频度关系中b值的极大似然法估计及其影响因素分析［J］.地震地质，2019，41（1）：21-43.

［5］ 常金龙.鹤岗及周边地区b值计算和地震危险性分析［J］.防灾减灾学报，2017，33（1）：32-35.

［6］ 潘华，高孟潭，谢富仁.新版地震区划图地震活动性模型与参数确定［J］.震灾防御技术，2013，8（1）：11-23.

［7］ 姚远，伊力亚尔·阿不力孜，吴传勇，等.2014年2月12日新疆于田M_S7.3地震发震构造初步研究［J］.内陆地震，2014，28（2）：121-126.

［8］ 谢江丽，阿布都瓦里斯·阿布都瓦衣提，黄帅堂.柴达木—阿尔金地震带b值统计计算［J］.震灾防御技术，2022，17（1）：95-103.

［9］ 汪素云，俞言祥.震级转换关系及其对地震活动性参数的影响研究［J］.震灾防御技术，2009，4（2）：141-149.

［10］金花，闫新义，任林，等.新疆地区近震震级与面波震级经验关系初步研究［J］.内陆地震，2021，35（2）：150-156.

［11］单新建，韩京，许静.新疆测震台网历史监测能力及现状［J］.内陆地震，1996，10（1）：61-67.

［12］王海涛，李莹甄，屠泓为.新疆历史地震目录完整性分析［J］.内陆地震，2006，20（1）：10-17.

［13］徐伟进，高孟潭.中国大陆及周缘地震目录完整性统计分析［J］.地球物理学报，2014，57（9）：2 802-2 812.

［14］任雪梅，高孟潭，冯静.地震目录的完整性对b值计算的影响［J］.震灾防御技术，2011，6（3）：257-268.

［15］冉慧敏，上官文明，赵晓成.2020年新疆于田M_S6.4地震及余震序列定位研究［J］.内陆地震，2020，34（3）：249-257.

［16］Keilis-Borok V I，Knopoff L，Rotvain I M.Bursts of aftershocks，long-term precursors of strong earthquakes［J］.Nature，1980，283：259-263.

［17］陈凌，陈颙，刘杰，等.地震活动性的统计分析：由过去推测将来的可能性研究［J］.地球物理学报，1998，41（1）：61-70.

［18］杨成荣，唐明帅.新疆于田7.3级地震及部分强余震震源参数的重新测定［J］.内陆地震，2010，24（1）：18-24.

［19］中华人民共和国国家质量监督检验检疫总局，中国国家标准化管理委员会.GB 17741-2005 工程场地地震安全性评价［S］.北京：中國标准出版社，2006.

［20］中华人民共和国国家质量监督检验检疫总局，中国国家标准化管理委员会.GB 18306-2015 中国地震动参数区划图［S］.北京：中国标准出版社，2016.

ANALYSIS OF B VALUE IN QAIDAM ALTUN SEISMIC

BELT WITH MULTIPLE SCHEMES

XIE Jiang-li，  Abuduwalisi Abuduwayiti，  ZHANG Wen-xiu

（Earthquake Agency of Xinjiang Uygur Autonomous Region，Urumqi 830011，Xinjiang，China）

Abstract： Based on seismic data from 1920 to 2020， b value of the Qaidam Altun seismic belt was calculated and analyzed with multiple schemes （different time periods and magnitude periods）， the b value was determined to be 0.78.It is lower than the value of 0.84 in engineering earthquakes comes from “Seismic ground motion parameters zonation map of China”， which is of great significance for seismic safety evaluation of major construction projects in seismic zones and engineering fortification of seismic zoning.As a supplement to “Seismic ground motion parameters zonation map of China”， b value determined in this paper can provide a scientific basis for the establishment of seismic hazard calculation model of major infra-structure.

Key words： Qaidam Basin-altun seismic belt；b value；Multiple schemes；Statistical calculations