2) 中国沈阳 110031 辽宁省地震局
2) Liaoning Earthquake Agency, Shenyang 110031, China
据中国地震台网正式测定,2023年12月18日23时59分(北京时间),甘肃临夏州积石山县(35.70°N,102.79°E)发生MS 6.2地震,震源深度为10 km。此次地震震中位于积石山县柳沟乡,距兰州市103 km。地震造成甘肃省117人遇难,781人受伤,155 393间房屋损毁;青海省伤亡212人,其中,死亡34人,受伤198人。根据甘肃省农业农村厅初步估算,此次地震造成甘肃农牧渔业直接经济损失5.32亿元。其中,种植业经济损失1.02亿元,畜牧业经济损失3.45亿元,渔业经济损失0.05亿元,农业系统办公用房0.8亿元(焦健等,2023)。根据青海省交通部门数据显示,积石山MS 6.2地震导致青海省3条高速公路、5条普通国省干线公路、58条农村公路共近300 km公路不同程度受损,造成损失7.37亿元(中国新闻网,2024)。
积石山MS 6.2地震发生后,中国地震台网中心立即启动应急响应,测震台网快速产出并陆续发布余震信息、震源区地震监测能力、序列精定位、震源机制、仪器烈度和库仑应力等数据产品。本文基于中国地震台网中心测震台网自主产出的地震数据产品,初步分析该地震的发震构造,为震源区地震危险性评估提供基础数据。
1 震源区地震背景积石山MS 6.2地震震源区附近历史上发生多次强度较大的地震,包括公元前193年临洮6—7级地震、公元前47年陇西6¾级地震、公元128年甘谷6½级地震、143年甘谷西7¼级地震、734年天水7级地震、1765年甘谷—武山6½级地震等(袁道阳等,2007;邵延秀等,2011)。据历史地震目录统计,自1900年以来,距此次积石山地震50 km范围内3.0级及以上地震共发生3次,震级范围3.0—3.9;100 km范围内3.0级及以上地震共发生33次,其中,6.0—6.9级地震1次,5.0—5.9级地震6次,4.0—4.9级地震1次,3.0—3.9级地震25次(图 1)。最大震级地震为1936年2月7日甘肃康乐西南6¾级地震,距此次地震震中约64 km(张波等,2015)。此次积石山MS 6.2地震的发生,打破了震源区100 km范围内继1936年2月7日甘肃康乐西南6¾级地震以来50多年无6级及以上地震发生的平静。
![]() |
图 1 甘肃积石山区域历史地震分布 Fig.1 Historical seismic map of the Jishishan region of Gansu |
在板块运动方面,积石山MS 6.2地震地震发生在青藏高原东北缘,该区受欧亚板块与印度板块的长期挤压作用,构造活动强烈,发育有西秦岭北缘断裂、倒淌河—临夏断裂、临潭—宕昌断裂、贵德断裂等多条深大断裂(邵延秀等,2011;王爱国等,2018)(图 2)。拉脊山断裂距此次地震震中最近,其位于祁连山褶皱带内,是西秦岭北缘断裂带的一部分(王爱国等,2018)。该断裂由拉脊山北缘断裂和拉脊山南缘断裂2条向NE凸出的弧形断裂所组成,分别长约230 km、220 km(袁道阳等,2005)。拉脊山北缘断裂西起日月山垭口山根村一带,向东沿西宁盆地南缘延伸,走向由NW60°渐变为近EW、NWW向,由多条不连续断裂带组成且活动性较高。自有地震记录以来至此次MS 6.2地震发生,沿拉脊山南北两侧发生20余次5级左右中小地震,无强烈地震发生。拉脊山南缘断裂由多条规模不等的次级断裂组成,其整体走向由近NWW向逐渐转为近EW向(涂德龙等,1998;袁道阳等,2005)。该断裂西起于日月山垭口克素尔村,向东终止于临夏大河家以南积石山前(李智敏等,2014)。拉脊山南缘断裂以千户村为界,分成东、西两段,东段多为直线状陡壁断崖,西段则多构成槽状负地形。
![]() |
图 2 2023年12月18日积石山MS 6.2地震背景构造 Fig.2 Background tectonic map of the Jishishan MS 6.2 earthquake on December 18, 2023 |
为了提升青藏高原地区地震监测能力,填补青藏交界地区监测能力空白,中国地震局自2019年起开始实施“青藏高原地震监测能力提升项目”,在我国西藏、青海、新疆3个地区建设地震台站,地震数据实时传输至中国地震台网中心(梁姗姗等,2021)。基于新建地震台站和原有中国测震台网波形数据,利用背景噪声和近震震级公式,对震源区地震监测能力进行计算,绘制积石山MS 6.2地震震源区地震监测能力图(图 3)。由图 3可见,震源区附近地震监测能力多在ML 1.5以上,与华北地区大部分区域地震监测能力一致。
![]() |
图 3 2023年12月18日积石山MS 6.2地震震源区监测能力 Fig.3 Earthquake monitoring capability in the source area with Jishishan MS 6.2 earthquake on December 18, 2023 |
2023年12月18日甘肃积石山MS 6.2地震余震频发,截至2023年12月29日23时59分,共记录余震749次(图 4),其中,3.0级及以上余震13次,分别是4.0—4.9级地震3次,3.0—3.9级地震10次;3.0级以下余震736次,分别是2.0—2.9级地震38次,1.0—1.9级地震149次,1.0级以下地震549次。目前,最大震级余震分别为2023年12月19日0时59分和12月21日4时2分发生的ML 4.1地震,距主震震中分别约14 km、16 km。此次甘肃积石山MS 6.2地震M—t图见图 4,余震序列震中分布见图 5。
![]() |
图 4 2023年12月18日积石山MS 6.2地震序列M—t图 Fig.4 M-T diagram of the sequence of the Jishishan MS 6.2 earthquake on December 18, 2023 |
![]() |
图 5 2023年12月18日积石山MS 6.2地震余震序列震中分布 Fig.5 Aftershocks sequence of the Jishishan MS 6.2 earthquake on December 18, 2023 |
选取波形信噪比相对较高且方位角覆盖分布相对较好的13个地震台站数据,采用gCAP方法进行震源机制反演(Zhu et al,2013)。所选台站震中距为73—210 km(图 6),在计算过程中,对这些台站记录的Pnl波和面波分别按30 s、60 s窗长进行截取,滤波范围为0.05—0.20 Hz,走向、倾角、滑动角的搜索间隔设为5°,震级搜索步长为0.1,深度步长为1 km,采用频率—波数法(F—K)计算格林函数,采样间隔设为0.1 s,采样点个数为1 024。地壳速度结构模型采用肖卓等(2017)对青藏高原东北缘及其邻区地壳速度结构的研究结果,波速比为1.70(表 1)。
![]() |
图 6 研究中震源机制计算和地震定位所用台站分布 Fig.6 Distribution of seismic source mechanisms and seismic positioning stations in this study |
![]() |
表 1 研究所用地壳速度结构模型 Table 1 Velocity model used in this research |
图 7给出积石山MS 6.2地震主震反演方差和震源机制解随不同深度取值的变化。由图 7可见,矩心深度为11 km时,震源机制解的均方根残差达最小值(RMS error = 46.49),相应的震源机制解为,节面Ⅰ:走向293°,倾角62°,滑动角52°;节面Ⅱ:走向172°,倾角46°,滑动角139°;矩震级为MW 5.96。
![]() |
图 7 积石山MS 6.2震源机制解均方根(RMS)残差随深度的变化 Fig.7 Variation of RMS error with focal depth during the focal mechanism in version for the MS 6.2 earthquake |
图 8示出最佳拟合深度处震源机制的理论波形与实际波形拟合情况。根据波形平均拟合相关系数对震源机制解结果质量进行分类(Zhu et al,2013),即平均拟合相关系数大于65%为A类,在60%—65%之间为B类,小于60%为C类。积石山MS 6.2地震震源机制最佳拟合深度处的平均波形拟合相关系数为73.3%,平均相关系数大于65%,属于A类,表明本研究震源机制解较可靠。用最小旋转角方法(万永革,2019)对比中国地震局地球物理研究所、美国地震调查局和全球矩心矩张量中心等国内外权威地震研究机构给出的积石山主震震源机制解参数结果(表 2),将得到的震源机制解为参考解,计算其与其他机构震源机制解的最小旋转角。由表 2可见,最小旋转角为13.79°—38.77°,表明本研究计算结果与其他研究机构较一致,进一步反映了本研究反演过程及结果的可靠性。
![]() |
图 8 2023年12月18日积石山MS 6.2地震主震双力偶机制解、矩张量解及理论波形(红色)与实际波形(黑色)拟合对比 波形下方数字表示理论波形相对实际波形的移动时间和二者间的相关系数,左侧大写字母为台站名,台站下方为震中距(km)和相对偏移时间(s) Fig.8 Double-coupled source mechanism solution and comparison between the theoretical (red) and observed (black) waveforms of the Jishishan MS 6.2 earthquake |
![]() |
表 2 不同机构反演的2023年12月18日积石山MS 6.2地震主震震源机制解 Table 2 The focal mechanism solutions of the 2023 Jishishan MS 6.2 earthquake from difference institutes |
采用中国地震台网提供的主震后27天内,即2023年12月18日至2024年1月14日的1 056个地震的观测报告进行地震序列精定位。选择震中距小于300 km的到时数据,且绝对走时残差为±2.0 s,经挑选,得到参与定位震相波形为:P波6 758条、S波4 463条。采用双差定位法,使用包含甘肃、青海、宁夏地震台网的24个地震台站观测数据(图 6)对此次地震进行重定位。
在地震定位过程中,S波震相因受到P波尾波干扰而不易识别,其读取精度要小于P波震相,因此,在双差重定位过程中,P、S波走时赋予先验权重分别设为1.0、0.5,事件对最大距离设为10 km。地壳速度结构模型与震源机制反演中所用模型相同。
采用双差定位方法,得到此次积石山MS 6.2地震序列623个地震事件的精定位结果(图 9)。由图 9可见,此次积石山MS 6.2地震主震及余震序列震源深度主要为5—17 km,均为浅源型地震;余震主要分布在主震西侧,呈近NNW方向展布特征;主震西北侧存在2个地震丛集,其地震深度浅于主震破裂初始深度;由剖面AA′可见,余震扩展范围约18 km;由剖面BB′可见,断层有略向东倾趋势。本研究结果与王世广等(2024)对该地震序列定位所得余震分布特征相似。此外,本研究所得震源机制节面Ⅰ走向与精定位地震序列空间分布的长轴取向一致。因此认为,节面Ⅰ(走向293°)为此次积石山MS 6.2地震的发震断层面。
![]() |
图 9 2023年12月18日积石山MS 6.2地震序列重定位后地震分布(a)及剖面图(b)、(c) 五角星表示积石山主震;圆圈表示余震,圆的大小对应震级大小;不同颜色代表余震相对主震的发生时间;AA′为沿余震走向的剖面;BB′为垂直余震走向的剖面 Fig.9 Top view(a) and cross sections(b and c) maps of the relocated Jishishan MS 6.2 earthquake sequence |
为深入理解余震扩展模式,将主震后2 h内发生的地震在水平空间展布出来(图 10),可见在该时间段内,余震主要沿NNW向展布,为单侧扩展模式,并在余震序列的NNW方向出现2簇小范围余震,其优势展布方向不明显,认为可能是主震触发所致。
![]() |
图 10 2023年12月18日积石山MS 6.2地震主震发生后2 h内余震震中空间分布 Fig.10 Earthquakes in map view within 2 hours after the mainshock |
Wang等(2003a)提出利用正交归一法计算地震应力场格林函数的方法,并在此基础上建立了粘弹松弛分层模型下的地震同震及震后形变的模型,并发展了相应的数值模拟方法。该方法同样可用于成层介质弹性半空间下的静态应力触发。基于上文所得余震精定位结果,余震深度范围为5—17 km,跨度较大,而库仑应力分布在不同深度上具有一定差别,因此分别计算5 km、10 km、15 km深度的库仑应力分布。
基于Crust1.0得到此次积石山地震所在地区地壳速度结构模型(表 3),结合623个地震的精定位结果,北京大学张勇研究组的徐晨雨等给出积石山地震的破裂模型(https://pku-geophysics-source.group/htmls/models.html),将节面Ⅰ的走向293°、倾角62°、滑动角52°设为接收断层参数,选取摩擦系数为0.4,使用PSGRN/PSCMP软件包(Wang,1999;Wang et al,2003a,2003b),计算得到此次地震在不同深度的库仑应力分布(图 11)。
![]() |
表 3 积石山MS 6.2地震所在地区成层地壳模型 Table 3 The layered crust model of the Jishishan MS 6.2 earthquake |
![]() |
图 11 积石山MS 6.2地震在5 km(a)、10 km(b)、15 km(c)深度的库仑应力分布 Fig.11 Coulomb stress distribution at depths of 5 km, 10 km and 15 km of the Jishishan MS 6.2 earthquake |
由图 11可见,甘肃积石山MS 6.2地震震源所在位置周边库仑应力变化均为负,说明地震的发生使震中周围地区长期积累的应力得到释放,库仑应力在5 km、10 km、15 km深度的最大卸载分别达0.24 MPa、0.86 MPa、1.70 MPa,同时,在沿发震断层NNW走向的两端有库仑应力增加的区域,在5 km、10 km、15 km深度上最大加载分别为0.40 MPa、0.30 MPa、0.22 MPa。从应力触发角度分析,大部分余震发生在库仑应力变化为负值的区域。这主要是文中只计算了主震的同震库仑应力变化,而主震的发生会使其周围地区长期积累的应力得到释放,从而使库仑应力变化为负。文中未计算震后库仑应力变化,而余震分布是主震后至2024年1月14日的精定位结果,有一定时间跨度。这些可能对应力触发结果产生影响。
4.4 地震烈度图随着国家地震烈度速报与预警工程的实施和运行,加速度计和烈度计记录数据已逐步得到应用。针对此次积石山MS 6.2地震,选取距震中300 km范围内站点波形资料,得到甘肃、青海、宁夏、四川共计1 306个台站(基本站140个、基准站197个、一般站969个)记录。通过对数据进行基线校正、记录转换、数字滤波处理,计算峰值加速度和峰值速度,得到积石山MS 6.2地震仪器烈度分布(图 12)。其中,GS/N0028台站地震仪器烈度最大,震中距为14 km。由图 12可见,震中NNW方向烈度值最大,与本研究地震精定位结果吻合度较好。需要说明的是,中国地震局正式发布的地震烈度图显示最大烈度为Ⅷ度(图 13)(https://www.cea.gov.cn/),而本研究计算的仪器烈度值为Ⅸ度,产生差异的原因是:①本研究计算的是某一点的烈度值,而调查烈度则为某一区域的综合烈度值,结果会有平均倾向;②本研究计算烈度时仅考虑土层影响,未考虑具体结构类型的影响,即烈度高值点若为抗震强度高的建筑,则实际调查烈度将偏低。
![]() |
图 12 2023年12月18日积石山MS 6.2地震仪器烈度 Fig.12 Distribution of instrumental intensity of the Jishishan MS 6.2 earthquake |
![]() |
图 13 2023年12月18日积石山MS 6.2地震烈度 Fig.13 Intensity map of the Jishishan MS 6.2 earthquake |
采用双差定位方法对2023年积石山MS 6.2地震序列进行精定位,结果显示,余震优势展布方向为NNW,长度约18 km,震源深度优势分布为5—17 km,深度剖面显示发震断层面倾向E。采用gCAP方法,获得此次地震震源机制解。结果表明,此次地震是一次逆冲型为主、兼有左旋分量的事件。初步推测节面Ⅰ为可能的实际发震断层面,发震断层为NNW走向(走向293°)的逆冲型断层。不同深度上的库仑应力分布结果显示,主震的发生使其周边区域出现较大范围的应力卸载。仪器烈度显示最大烈度值为Ⅸ,且烈度区域较小,集中分布在主震NNW方向。结合周围地质构造,认为此次地震的发震断层为拉脊山南缘断裂。
自2016年青海门源MS 6.4地震(梁姗姗等,2017)发生以来,在青藏高原东北缘相继发生2021年玛多MS 7.4(徐志国等,2021)、2022年德令哈MS 6.0、2023年门源MS 6.9(Liang等,2023)、2023年积石山MS 6.2地震,这些中强震的频繁发生,再次表明印度板块与欧亚板块碰撞导致青藏高原的快速隆升,强烈构造活动持续,致使青藏高原东北缘地区地震活动强烈。
本文撰写过程中,得到应急管理部国家自然灾害防治研究院张广伟副研究员,防灾科技学院靳志同教授,中国地震台网中心邹立晔高级工程师、刘艳琼高级工程师的指导和帮助,在此表示感谢。
丁国瑜, 卢演俦. 中国岩石圈动力学概论[M]. 北京: 地震出版社, 1991: 142-153.
|
冯建刚, 蒋长胜, 韩立波, 等. 甘肃测震台网监测能力及地震目录完整性分析[J]. 地震学报, 2012, 34(5): 646-658. DOI:10.3969/j.issn.0253-3782.2012.05.006 |
焦健, 王鹏. 积石山地震已造成甘肃农牧渔业直接经济损失5.32亿元[EB/OL]. (2023-12-23)[2024-01-15]. https://china.qianlong.com/2023/1223/8170379.shtml.
|
康来迅. 西秦岭北缘断裂带活动构造地貌特征及其形成机理[J]. 华东师范大学学报(自然科学版), 1994(2): 67-75. |
李智敏, 李延京, 田勤俭, 等. 拉脊山断裂古地震与喇家遗址灾变事件关系研究[J]. 地震研究, 2014, 37(Z1): 109-115. |
梁姗姗, 雷建设, 徐志国, 等. 2016年1月21日青海门源MS 6.4余震序列重定位和主震震源机制解[J]. 地球物理学报, 2017, 60(6): 2 091-2 103. |
梁姗姗, 邹立晔, 赵博, 等. 中国测震台网地震监测能力初步分析[J]. 地震地磁观测与研究, 2021, 42(6): 68-75. |
刘巧霞, 杨卓欣, 莘海亮, 等. 玉树MS 7.1级地震部分余震重新定位及发震构造分析[J]. 地球物理学报, 2012, 55(1): 146-154. |
邵延秀, 袁道阳, 王爱国, 等. 西秦岭北缘断裂破裂分段与地震危险性评估[J]. 地震地质, 2011, 33(1): 79-90. |
涂德龙, 王赞军, 曾包红, 等. 青海省湟水盆地全新世活动断裂分布及其活动特征研究[J]. 西北地震学报, 1998(4): 83-90. |
万永革. 同一地震多个震源机制中心解的确定[J]. 地球物理学报, 2019, 62(12): 4 718-4 728. |
王爱国, 袁道阳, 张波, 等. 西秦岭北缘断裂带地震危险性与危害性[J]. 城市与减灾, 2018(4): 48-54. |
肖卓, 高原. 利用双差成像方法反演青藏高原东北缘及其邻区地壳速度结构[J]. 地球物理学报, 2017, 60(6): 2 213-2 225. |
徐志国, 梁姗姗, 张广伟, 等. 2021年5月22日青海玛多MS 7.4地震发震构造分析[J]. 地球物理学报, 2021, 64(8): 2 657-2 670. |
许英才, 郭祥云. 2022年6月1日四川芦山MS 6.1地震的震源参数及其构造启示[J]. 地球物理学报, 2023, 66(8): 3 202-3 217. |
于晓. 积石山6.2级地震最后两名失联人员已找到, 青海海东34人遇难[EB/OL]. (2023-12-31)[2024-01-15]. https://www.chinanews.com/sh/2023/12-31/10138307.shtml.
|
袁道阳, 雷中生, 葛伟鹏, 等. 对143年甘谷西7级地震史料的新见解[J]. 西北地震学报, 2007, 29(1): 58-63. |
袁道阳, 张培震, 雷中生, 等. 青海拉脊山断裂带新活动特征的初步研究[J]. 中国地震, 2005, 21(1): 93-102. |
张波, 何文贵, 方良好, 等. 1936年甘肃康乐6(3/4)级地震地表破裂带调查[J]. 地震研究, 2015, 38(2): 262-271. |
张波, 何文贵, 袁道阳, 等. 西秦岭北缘断裂带西端晚第四纪活动特征及其西延问题[J]. 地震, 2012, 32(1): 136-143. |
张广伟, 雷建设, 梁姗姗, 等. 2014年8月3日云南鲁甸MS 6.5级地震序列重定位与震源机制研究[J]. 地球物理学报, 2014, 57(9): 3 018-3 027. |
中国新闻网. 积石山6.2级地震甘肃遇难人数升至117人现场排查搜救结束[EB/OL]. (2023-12-22)[2024-01-15]. https://www.chinanews.com.cn/sh/2023/12-22/10133726.shtml.
|
中国新闻网. 积石山6.2级地震致青海近300公里公路不同程度受损[EB/OL]. (2024-01-11)[2024-01-15]. https://www.sohu.com/a/751245438_121123826.
|
邹立晔, 黄志斌, 周静, 等. 2019年6月17日四川长宁M 6.0地震的快速测定与数据产品产出[J]. 中国地震, 2019, 35(3): 573-583. |
Liang S S, Xu Z G, Zhang G W, et al. The Menyuan, Qinghai MS 6.9 earthquake on January 8, 2022: a strike-slip cascading rupture event[J]. Pure and Applied Geophysics, 2023, 180(3): 829-850. |
Wang R J. A simple orthonormalization method for stable and efficient computation of Green's functions[J]. Bulletin of the Seismological Society of America, 1999, 89(3): 733-741. |
Wang R J, Küumpel H J. Poroelasticity: efficient modeling of strongly coupled, slow deformation processes in multilayered half-space[J]. Geophysics, 2003a, 68(2): 705-717. |
Wang R J, Lorenzo Martín F, Roth F. Computation of deformation induced by earthquakes in a multi-layered elastic crust-FORTRAN programs EDGRN/EDCMP[J]. Computers & Geosciences, 2003b, 29(2): 195-207. |
Zhu L P, Ben-Zion Y. Parametrization of general seismic potency and moment tensors for source inversion of seis-mic waveform data[J]. Geophysical Journal International, 2013, 194(2): 839-843. |