文章速递|青年编委李宗超等:2025年1月7日西藏定日MS6.8地震珠峰地震台加速度时程快速模拟

学术   2025-01-25 00:00   湖北  

点击蓝字 关注我们






202517日,西藏日喀则市定日县发生MS6.8地震,震中附近遭受了严重的地震破坏。地震影响区位于海拔4000米以上的高原地带,地震台站分布稀少,100公里内缺乏完整的基准台或基本台地震记录,而这些记录对于快速评估地震灾害程度和指导震后重建有至关重要的意义。本次地震影响区域内的高山多覆盖有永久冰川(如珠峰地区),同时存在大量潜在的滑坡体,研究地震对冰川及滑坡体稳定性的影响,是值的关注的重要科学问题。基于上述背景和研究需求,本文选择距离震中55公里的珠峰地震台记录到的一次MS3.9小震作为格林函数采用经验格林函数方法,考虑参数不确定性,合成了本次地震在珠峰地震台的加速度波形。研究结果表明:珠峰地震台在本次地震中水平向PGA范围为70-190 cm/s2,垂直向PGA范围为40-140 cm/s2,对应的仪器地震烈度为Ⅵ-Ⅶ度。这些结果反映了珠峰地震台位置的地震动强度,与第五代中国地震动参数区划图中不同概率水准的地震动反应谱进行对比,并结合当地房屋建筑的结构性能,分析认为该强度地震动可能对珠峰地震台附近的一般民用建筑造成较严重损伤,震中附近的地震破坏可能会更为严重。同时该地震可能对高海拔冰川和潜在滑坡体稳定性造成影响,建议对此问题给予特别关注。



1

引言


据中国地震台网测定:2025179516秒,西藏自治区日喀则市定日县发生MS6.8地震,震中位于北纬28.50˚、东经87.45˚,距离定日县城约36公里,震源深度9.3公里。截至20251140800分,共记录到余震3614次,其中2.0级以下3104次,2.0-2.9458次,3.0-3.945次,4.0-4.96次,5.0-5.91。此次地震余震分布呈现南北向分布特征(1)。

通过震后现场震害调查,并结合震区断裂构造、仪器烈度、余震分布、震源机制和遥感等数据,中国地震局发布了西藏定日6.8级地震烈度图(http://www.mem.gov.cn)。此次地震引发约26公里的地表破裂,最高烈度为Ⅸ度(9度),Ⅵ度(6度)区及以上面积约24万平方公里,等震线长轴呈近南北走向,长轴约191公里,短轴约152公里,共涉及西藏自治区1个市6个县45个乡镇。此次地震造成126名同胞不幸遇难,约6.2万人不同程度受灾。

图1 定日地震震中位置、宏观地震烈度、历史地震和余震目录。历史地震和余震目录来源于中国地震台网,最后一次访问时间为2025年1月16日14时。宏观地震烈度来源于http://www.mem.gov.cn。
震后,中国地震预警台网迅速公布了震区的台站分布情况https://www.iem.net.cn。在震中100公里范围内,仅有珠峰地震台(XZ-ZHF)和日喀则地震台2个地震台站,其中仅日喀则地震台(1)记录到了完整的强震波形,珠峰地震台(XZ-ZHF)因故未能记录到本次地震的完整波形。这意味着震中西部和南部等方位缺少完整的强震记录,导致本次地震近场强震特征存在较大不确定性。

近场大震记录极为珍贵,而中国在这一方面,尤其是西部和高原地区的近场强震记录相对缺乏。自2000年以来发生多次7-8级的强震事件未能获得完整的近场强震记录。例如,2021年的青海玛多MW7.4地震在170公里范围内无完整强震记录(李宗超等,2022);2001年昆仑山大地震同样缺少近场强震记录。鉴于以上情况,我们提出了合成定日地震在珠峰地震台主震加速度波形的研究设想。主要采用经验格林函数方法,选取珠峰地震台之前记录到的一次MS3.9小地震作为格林函数,并适当考虑了参数的不确定因素,从地震应急的角度初步合成了珠峰地震台位置处主震的加速度波形。

珠峰地震台位于西藏日喀则地区定日县扎西宗乡,距珠峰大本营约40公里,海拔4255米。作为中国地震预警台网的基准台,珠峰地震台配备了测震仪和强震仪。尽管西藏日喀则南部地区中强地震活动较为频繁,但此前一直未设立地震台。珠峰地震台的建成填补了日喀则南部地区的地震监测空白,并提升了我国中尼边界地区的地震监测能力。此外,它为地学工作者提供了宝贵的一手地震观测数据,显著增强了我国在青藏高原国际地学及防震减灾领域的科技竞争能力。自2010年运行以来,珠峰地震台已记录了大量小震波形数据,为多学科的发展提供了重要的数据支撑。

本文认为,合成珠峰地震台的加速度波形有以下两方面的重要科研应用:(1)服务珠峰地区冰川及滑坡体稳定性研究。根据我国第二次冰川编目和尼泊尔冰川编目(刘时银等,2020显示,珠峰地区南、北坡有2438条冰川,总面积达3271.4平方公里,其中,我国境内珠峰自然保护区分布有1476条冰川,面积2030.5平方公里。此次破坏性地震作用可能对珠峰地区的冰川稳定性产生影响。2017年米林MS6.9地震曾对色东普冰川稳定性造成了一定影响,导致其一年后发生冰崩灾害(Li et al., 2022; Hu, 2019)。通过获取珠峰地区较为准确的加速度波形数据,可为本区域冰川及滑坡体灾害风险预测提供科学依据。(2)为珠峰地区工程建筑抗震设计提供支持。珠峰地震台位于边境区域,区域内的相关工程建设需服务于国防安全,因而需要准确的地震动输入。地震工程学界普遍认识到,近场强地震动预测对于减轻未来大地震灾害及指导重大工程建设抗震设计至关重要,特别是高频率地震动的科学评估具有极高的应用价值。参考第五代地震动参数区划图的标准(GB 18306-2015),该地区设计地震动为0.2g,设防水准为度,表明其地震危险性较高。珠峰地区准确的地震波形数据将为工程建筑的抗震设计提供可靠的数据基础,有效提升防灾减灾能力。


2

珠峰波形模拟方法及相关参数


小地震记录能够捕捉真实的震源破裂过程、传播路径介质以及浅表场地响应的复杂性。日本地震学家IrikuraMiyake等人的研究表明 (Miyake et al., 2003)利用震源区附近的小地震记录是模拟强地震动的有效方法之一(2)。因此,基于小地震记录的合成地震动在一定程度上保留了这些复杂影响,这也是经验格林函数方法的优势之一。尤其是在小地震记录丰富的地区,通过考虑震源参数的不确定性,利用经验格林函数法可以较为精确地模拟破坏性地震的近场高频强地面运动特征

本文的研究目标是合成定日MS6.8地震在珠峰地震台位置处的加速度波形。我们选择了一个MS3.9的小地震作为格林函数(2)。该小震的震中与主震位置高度一致,且珠峰地震台完整记录到这一小地震的波形。通过使用珠峰地震台的实际记录数据作为格林函数,保证了大小地震传播路径和场地条件的一致性,这对于提高合成波形的准确性尤为关键。

本研究参考了北京大学张勇团队和美国地质调查局(USGS)提供的震源破裂过程相关参数(3),包括断层位错、震源上升时间、破裂模式及震源机制等。这些参数为模拟珠峰地震台位置处的强地面运动提供了科学依据。小震数据来源于李宗超等人建成的格林函数数据库和地震科学国际数据中心。在地震动合成过程中,我们还综合考虑了震后地震局系统和高校科研单位的地震应急产出,结合震源参数的不确定性,最终给出了合理取值范围内的927加速度波形。本次地震中,珠峰地震台与震中的距离为54.9公里,方位角为253°,所选用的格林函数震级为MS3.9,小震震中为(东经87.50°北纬28.5°),与主震几乎一致。

小震合成大震示意图及本次合成珠峰地震台主震波形所有的格林函数

图3 定日地震强震破裂模型. (a)张勇团队破裂模型;(b) USGS破裂模型


3

珠峰地震台加速度波形合成结果



经验格林函数方法的主要优势体现在大于1.0 Hz的频段。该方法在这一频段内能够较好的合成地震动,其高频加速度幅值能够有效代表宽频带地震动加速度的幅值(Ji et al., 2023)。通过该方法,本研究最终合成了9个三分量的加速度波形(4a)以及对应的频率>1.0 Hz的加速度反应谱(4b。结果表明:合成的水平向地震动的PGA取值范围为(70-190cm/s2,竖直向地震动的PGA范围为(40-140cm/s2。根据中国地震烈度(GB/T 17742-2020)的判定,合成地震动对应的地震烈度为-Ⅶ度,以度为主。模拟结果与中国地震局公布的定日地震烈度空间分布高度一致。根据野外地震调查结果,珠峰地震台区域的烈度位于度区,模拟结果与实际调查的烈度反映了相同水平的破坏能力。本研究通过合成的加速度波形和烈度分析验证了经验格林函数方法在珠峰地震台区域的适用性,为后续地震动特征研究和抗震设计提供了可靠的理论支持和参考数据。

(4a)

(4b)

珠峰地震台的模拟地震动.4a)黑色曲线表示EW分量,蓝色曲线表示NS分量,红色曲线表示UD分量;(4b) 珠峰地震台的模拟地震动的加速度反应谱.a)黑色曲线表示EW分量,(b)红色曲线表示NS分量,(c)蓝色曲线表示UD分量.

本文采用了珠峰地震台测震仪记录的宽频带小震波形作为格林函数,合成了珠峰地震台基岩水平向加速度波形。需要注意的是,自由地表或者Ⅱ类场地的地震动一般因为场地放大效应显著高于4中展示的基岩水平地震动强度。结合第五代地震动参数区划图,将合成的基岩水平地震动与珠峰地震台区域Ⅱ类场地不同概率水准的设计地震动进行对比(5)发现:(1)合成的基岩水平地震动强度已经超过基本地震动强度水平;(2)若考虑场地放大效应,自由地表的地震动强度可能会接近或者超过罕遇地震动的强度水平。本文认为该水平的地震动破坏能力对珠峰区域的冰川稳定性可能会造成不确定性的影响,应当被引起重视。同时参考本地区房屋建设水平,该强度的地震动强度会对本地区的一般民用建筑造成一定的损伤。

珠峰地震台的模拟基岩层面的地震动的加速度反应谱与Ⅱ类场地设计谱


4

讨论


本文选取珠峰地震台记录到的与主震震中位置高度一致的一次MS3.9小震作为格林函数,采用经验格林函数方法,并适当考虑参数不确定性,合成了本次地震在珠峰地震台(XZ.ZHF)位置处的加速度波形。通过将合成结果与本地区不同概率水准的设计反应谱进行对比,并结合当地房屋建设水准,发现该强度的地震动可能对珠峰地震台附近的一般民用建筑造成较严重损伤。此外,地震可能影响本地区高海拔冰川及潜在滑坡体的稳定性,应引起特别关注。为了清晰展示研究内容和结果,本文采用6所示的成果展示模式。

本研究主要基于定日地震后的应急产出反演结果,目前部分参数仍有进一步提精度的空间,随着相关研究的深入,定日地震的震源参数将更加精确,我们也进一步优化珠峰地震台加速度波形的合成精度。本研究揭示了小震波形积累的重要性及其在强震动地震学中的应用潜力。中国已经建成全球最大规模的预警台网,全国基准台和基本台数量已超过5000个,未来高质量地震台站的数量还将继续增加。本研究充分发挥现有小震记录应用潜力,不仅致力于模拟历史大地震在高质量台站处的加速度波形,还可用于未来可能发生的破坏性情景地震的模拟需要高精度地震动输入的各类科研工作提供支持

图6 最终成果展示

文章信息:

李宗超, 孙吉泽, 纪志伟, 高孟潭, 谢俊举, 赵延娜, 2025. 2025年1月7日西藏定日MS6.8地震珠峰地震台加速度时程快速模拟. 地球科学. doi: 10.3799/dqkx.2025.009

LI Zongchao, SUN Jize, JI Zhiwei, GAO Mengtan, XIE Junju, ZHAO Yanna, 2025. Rapid simulation of acceleration waveform at Everest Seismic Station of Xizang Tingri MS6.8 earthquake on January 7, 2025. Earth Science. doi: 10.3799/dqkx.2025.009

更多论文信息,请点击左下方(阅读原文


加入学术交流群

备注“姓名+单位+ES”



扫码添加小编微信,拉你进入《地球科学》编辑部学术交流群,群内不定期开展答疑会、分享学术研究、学术直播等高质量活动。


往期回顾

矿物学·岩石学·矿床学

   →近70年中国变质岩石学-变质地质学的研究进展

造山带橄榄岩岩石学与构造过程:以松树沟橄榄岩为例

片麻岩穹窿与伟晶岩型锂矿的成矿规律探讨

岩浆作用与地球深部过程

从岩浆岩看青藏高原地壳的生长演化

地壳生长及深部物质架构研究与问题:以中亚造山带(北疆地区)为例

→辽东白云-小佟家堡子矿集区控矿构造及成矿有利区预测

→晶粥储存、侵入体累积组装与花岗岩成因

→交代岩石圈地幔与金成矿作用

实验地球科学的前沿与发展战略

实验地球化学的发展历史和研究展望

实验岩石学发展现状与趋势

实验矿床学的发展现状和前景展望

实验流变学的发展现状与趋势

高温高压装置研制和技术创新的发展现状与趋势


构造地质学

纪念马杏垣先生

关于大陆构造研究的一些思考与讨论

秦岭造山带是印支碰撞造山带吗?

关于发展洋板块地质学的思考

大陆俯冲带壳幔相互作用

增生型造山带结构解析与时空制约

洋板块地质与川藏铁路工程地质关键区带

华北克拉通太古宙构造热事件时代及演化

中条山地区涑水杂岩新太古代烟庄正长花岗岩年龄及成因:对华北克拉通地壳演化的制约

→俯冲带结构演变解剖与研究展望


地学与健康

"同一健康"视角下医学地质学的创新发展

优化生态环境保障人民健康


古生物学

地史时期生物对冰室气候形成的作用


沉积和盆地动力学

下扬子地区早古生代晚期前陆盆地沉积特征与盆山过程

陆相断陷盆地陆源碎屑与碳酸盐混合沉积研究进展

渤海海域秦南凹陷东南缘沙一段混积岩优质储层成因及勘探意义


石油地质学

非常规源岩层系油气形成分布与前景展望

→”双碳”目标下我国油气产业发展的思考

→我国陆上未来油气勘探领域探讨与攻关方向


地球物理学

从高分辨率地震层析成像看青藏高原软流圈的物质运动


《地球科学》

青年学者论坛回放

→黄康俊|金属稳定同位素示踪地球增氧事件

→张   文|激光微区锆石Zr同位素分析测试技术及参考物质研发

→Marize Muniz da Silva|巴西São Francisco地块 Mineiro带后碰撞岩浆作用

→Dionatan Ferri Padilha|巴西南部Dom Feliciano带后碰撞玄粗质和超钾质岩浆作用

→王选策|地球系统中的挥发分和物质循环——对非俯冲带岩浆的新观点

→陈   印|中国北方砂岩型铀成矿作用及其与构造演化和油气藏的联系

→荣   辉|基性岩侵入对铀储层成岩-成矿作用的影响:以松辽盆地钱家店铀矿床为例

→兰中伍|南秦岭岩浆作用与浅表环境响应

→余文超|贵州遵义下石炭统九架炉组铝土矿物源:来自地球化学与锆石U-Pb同位素年代学的证据

→岳  汉|从2021年青海玛多地震到松潘-甘孜块体的长期形变模式

→王  墩|2021年Mw7.4青海玛多地震的震级、地震动场及破裂过程

→邓阳凡|华北东南缘综合地球物理特征及其对徐淮弧形构造带的启示

→谢国治|胶东半岛高Ba-Sr牙山侵入体的成因及其与华北克拉通破坏的关系

→白江昊|长江中下游地区A型花岗岩成因:放射性成因Ca-Nd同位素制约

→李振洪|川藏交通廊道重大灾害防控专题

→王   栋|川藏交通廊道冰雪灾害风险评价

→汤高明|青藏高原典型地表地质灾害机理及分类

→崔一飞|川藏交通廊道泥石流级联调控风险防控原理与方法

→刘金高|金伯利岩与碱性玄武岩的成因联系及深部地幔演化

→丁   兴|富氟流体活动及其成岩成矿效应 

柳保军,雷超,郑金云 | 专题十:南海陆缘盆地动力学及其能源效应

施小清,杨晓帆,戴恒 | 专题十一:地下水环境的刻画与模拟

高晓英,罗涛 |专题十二:微区分析在地球科学中的应用

梁琛岳,刘博 | 专题十中亚造山带南缘构造与成矿

Alexander Humphreys ,Johann Hohenegger | 专题十四:热带大型底栖有孔虫

 付伟 ,李登峰,李占轲,等 | 专题十五:战略性矿产资源研究与勘查进展

雷宏武,张晓博 | 论坛专题十六:地热能地球科学

李斌 ,张志遥 | 专题十七 :深层油气论坛

李博 ,吴琼,李严严 | 专题十八:强震区滑坡崩塌灾害防治技术方法研究

程晓敢,石许华,林秀斌 | 专题十九:盆山过程与动力学

黄发明,张东明,李旭 | 专题二十 :机器学习与灾害风险

郑义,虞鹏鹏 | 专题二十一:华南大地构造演化与关键金属成矿

黄发明,李为乐,谭道远 | 专题二十二:地质灾害风险智能评价

陈红汉,罗群,李文正,王海学,唐大卿 | 专题二十三:走滑断裂与油气

侯晓坤,贾智杰 | 专题二十四 :工程地质与防灾减灾

王全荣,杜尧 | 专题二十五:地下水系统参数获取与污染成因

董志文,王欢业,蒋宏忱 | 专题二十六:青藏高原环境-地质-微生物

喜欢就点个在看吧!

地球科学编辑部
《地球科学》(中文版)和Journal of Earth Science(英文版)是由教育部主管,中国地质大学主办的综合性地球科学学术刊物。中文版被EI 等数据库收录,英文版被SCI收录。
 最新文章