Strong ground motion simulation and intensity distribution estimation for the MS6.9 Menyuan，Qinghai，earthquake on 8 January 2022
-
摘要: 2022年1月8日青海省海北藏族自治州门源回族自治县发生MS6.9地震。门源地震序列的重定位结果认为门源地区还存在一定的应力积累，未来该地区具有发生强震的可能。本文结合震源区地形数据、三维速度结构，根据门源地震震源破裂过程的初步结果，采用曲线网格有限差分方法模拟了门源地震的波场传播过程，得到烈度分布。结果表明：沿平行断层走向方向的地震动衰减明显小于垂直断层走向方向；门源地震的最大烈度为Ⅷ度，位于震源破裂起始点附近区域，理论烈度与野外调查的地震烈度分布基本一致；受强地面运动方向性效应和起伏地表的影响，地震灾害主要沿发震断层的WNW方向和ESE方向集中分布。
-
关键词:
- 门源地震 /
- 曲线网格有限差分方法 /
- 强地面运动模拟 /
- 地震烈度
Abstract: A MS6.9 earthquake occurred at Menyuan County, Haibei Prefecture, Qinghai Province on 8 January 2022. Relocation results of the Menyuan MS6.9 earthquake sequence suggest that there is still a possibility for another strong earthquake in Menyuan area. Based on topography data and three-dimensional velocity structure in the source zone, together with the preliminary result about the source rupture process of Menyuan earthquake, we simulated seismic wave propagation of Menyuan MS6.9 earthquake by curvilinear-grid finite difference method, and obtained seismic intensity distribution accordingly. The simulated results show that the ground motion attenuation along the direction of fault strike is significantly smaller than that along the direction of fault dip. The maximum seismic intensity caused by Menyuan MS6.9 earthquake is about Ⅷ degree, and located in the vicinity of the rupture initial point. Theoretical intensity distribution is roughly consistent with the seismic intensity distribution from field surveys. Affected by the directional effects of strong ground motion and undulated topography, earthquake disasters mainly concentrated in the WNW and ESE directions of the seismic fault. -
引言
2022年1月8日01时45分27秒青海省海北藏族自治州门源回族自治县发生MS6.9地震，震中位置为（37.77°N，101.26°E），震源深度为10 km。本次地震造成了17 069人受灾，由于震中距离人口稠密地区较远，尚无人员伤亡的相关报道（青海日报，2022），兰新高铁浩门至山丹军马场区间隧道群由于此次地震而发生局部塌方（央视网，2022）。门源MS6.9地震发生在青藏高原东北缘冷龙岭断裂、托莱山断裂和肃南—祁连断裂的阶区部位。现场工作队在冷龙岭断裂带西段探测到长约22 km的地表破裂带（青海省地震局，2022）。托莱山断裂和冷龙岭断裂现为左旋兼挤压断裂（李强等，2013），历史上无M7以上地震记载（姜文亮，2018），最近一次强震为2016年1月21日门源MS6.4地震，该地震的发生反映了青藏高原地块向NE向不断推挤生长的过程（胡朝忠等，2016）。门源MS6.9地震发生后，许英才等（2022）对门源地震早期序列（2022年1月8日至12日）进行了重定位和震源机制研究，其结果表明目前门源地区还存在一定的应力积累且应力尚未得到充分释放，该地区仍有发生强震的危险。由于地震灾害主要是地震在地表产生的强地面运动造成的，因此为了获得地震波场传播过程及其引起的地表响应，对门源MS6.9地震开展强地面运动初步模拟及烈度估计具有重要意义。
地震波场的正演模拟需要考虑复杂地表。曲线网格有限差分方法（Zhang，Chen，2006；Zhang et al，2012）适用于含起伏地形的地震波场模拟，该方法中地表的形状用任意曲线网格近似，并采用牵引力镜像法处理自由表面条件。曲线网格有限差分的强地面运动模拟方法在地震后的灾害评估中得到了非常广泛的应用。Zhang等（2008）结合三维介质模型、震源破裂模型和地表地形数据模拟了2008年汶川MS8.0地震的强地面运动，研究表明断层破裂方式和盆地构造控制了地表峰值速度（peak ground velocity，缩写为PGV）的分布，地表起伏剧烈地区往往对应较大的PGV数值，应关注其震害问题。张振国等（2014a，b）对2014年2月12日新疆于田MS7.3地震和2014年8月3日云南鲁甸MS6.5地震引起的强地面运动作了初步模拟和烈度预测，研究显示地震动在山峰、山脊处具有较大幅值，该结果可指导震后的灾区重建工作。赵宏阳和陈晓非（2017）利用1975年2月4日辽宁海城MS7.3地震的地震地质资料，模拟计算了海城地震的波场传播过程，分析了强地面运动的方向性效应、盆地效应和近断层效应，得出理论烈度分布同震后调查烈度分布基本一致，验证了研究所用的震源模型和速度结构的合理性。
本文拟根据张勇①提供的门源地震震源破裂过程的初步结果，利用曲线网格有限差分方法模拟门源MS6.9地震的强地面运动，再结合地表峰值速度和烈度间的关系，计算地震烈度，并在此基础上评估地震灾害分布特征，以期为门源地区的防震减灾提供科学依据。
1. 研究区概况及计算方法和模型
研究区范围如图1所示，可以看出，该区域的地貌特征变化大，地势西南高东北低，高程介于1.3—5.0 km之间。门源MS6.9地震发生后至2022年1月20日，该地区已发生MS≥4.0余震23次，其中MS≥5.0余震2次。
![]()
图 1 研究区及周边的地质概况和2022年1月8—20日震源区MS≥3.0 地震分布Figure 1. Geological overview of the study area and its surrounding regions，spatial distribution of MS≥3.0 earthquakes in source zone during the period of 8 to 20 January，2022采用曲线网格有限差分方法（Zhang，Chen，2006；Zhang et al，2012）对门源MS6.9地震的强地面运动进行模拟。计算中需要设置描述地形起伏的网格模型、反映地下物质属性的介质模型，以及表示地震破裂过程的震源模型。
地形选取GTOPO30地形数据，其水平分辨率大约为1 km。整个计算区域尺度为350 km×220 km，深度为60 km。将研究区域离散成700×440×120个网格，垂直方向采用等间距排列，单位网格为边长500 m的立方体。
考虑到面波对纵波速度的灵敏度和起伏地形的影响，Han等（2022）提出一种改进的体波和面波数据联合反演方法获得中国大陆地壳和上地幔水平分辨率为0.5°的速度结构（USTClitho2.0模型）。该模型对于青藏高原地块内体波射线覆盖相对稀疏的区域，加入面波频散数据可更好地提高纵波速度和横波速度的准确度，因此本文以USTClitho2.0模型为基础建立研究区域的网格化横波速度和纵波速度模型。每个网格的东向、北向和垂直方向的长度分别为43.75 ，55 和5 km，同时依据Brocher （2005）提出的密度和纵波速度的经验关系确定各网格的密度。模拟采用的横波速度结构如图2所示，可以看出，2—12 km深度存在高速层，沿北东方向地壳厚度逐渐减薄。虽然青海湖处于研究区域内（图1），但是由于其位于模拟计算区域的边缘，且水深较浅，因此在建立介质模型时忽略其影响。
张勇①反演了该地震的震源破裂过程，结果显示断层滑动出露地表，走向为103°，倾角为88°，滑动角为−6°，共有31×11＝341个子断层，子断层的空间分辨率为2 km×2 km，采样点为80个，每间隔0.25 s给出一个滑动速率，破裂过程持续时间约20 s，最大滑动量为1.8 m。基于这一结果，本文设置了相应的震源模型，如图3所示。
![]()
图 3 本文所用的门源MS6.9地震震源模型Figure 3. Source model of Menyuan MS6.9 earthquake used in this study2. 强地面运动模拟结果
强地面运动模拟计算的时间步长为0.01 s，总步数为1万步，模拟时长共100 s，使用640个计算核心进行计算。门源地震x分量的模拟速度场快照如图4所示，可见：地震的主要能量由位于震中附近下方的位错产生，且以水平走向错动为主，因此在断层垂向上产生了能量较强的S波；大约第2.36 s时，地震波到达地表，在初 始破裂时刻，速度的最大值集中在断层破裂的前锋上；随着远离发震断层，地震波场能量逐渐减小。
![]()
图 4 门源MS6.9地震x分量粒子速度波场快照Figure 4. Wavefield snapshots of the x component particle velocity of Menyuan MS6.9 earthquake通过强地面运动模拟得到各网格点x，y，z三个方向不同时刻的速度，然后对三个方向的速度分量求矢量和获得各网格点的运动速度随时间的变化，根据各网格点的速度最大值确定研究区域的地表峰值速度（PGV）分布，并结合国家市场监督管理总局和国家标准化管理委员会（2021）给出的PGV与烈度之间的关系得到对应的地震烈度，如图5所示，结果表明：沿平行断层走向方向的地震动衰减明显小于垂直断层走向方向；地震的最大烈度为Ⅷ度（PGV＝37.63 cm/s），位于震源破裂起始点附近区域；Ⅷ度区主要涉及门源县、祁连县和肃南县的部分区域；Ⅶ度区主要涉及大通县、永昌县、民乐县等部分区域。中国地震局（2022）野外调查的烈度分布显示等震线长轴呈WNW走向，同理论模拟显示的高烈度区主要沿WNW方向延伸的结果一致。理论烈度分布与中国地震局工程力学研究所强震动观测组（2022）给出的仪器烈度分布（图6）较为接近，即高烈度区主要沿断层走向展布，且断层南侧的地震烈度大于北侧。受模型分辨率和计算成本的限制，本文模拟的地震波场的最高有效频率为0.4 Hz （2.5 s），缺少高频成分，因此获得的烈度较强震仪和烈度仪观测的结果偏低，但相较于使用低分辨率介质模型的模拟结果（徐剑侠等，2015），地震动模拟的最高频率得到了一定的提升。
![]()
图 6 门源MS6.9地震仪器观测的烈度分布（引自中国地震局工程力学研究所强震动观测组，2022）Figure 6. Instrumental seismic intensity distribution of Menyuan MS6.9 earthquake （after Strong Motion Observation Group，Institute of Engineering Mechanics，China Earthquake Administration，2022）地震灾害区域集中在发震断层附近，并向WNW方向和ESE方向延伸，这主要由两个原因造成：一方面，该地震为高倾角的走滑型地震，地震波能量主要沿走向方向传播，表现出强地面运动的方向性效应；另一方面，起伏地表对地震波传播具有重要影响（Zhang et al，2008），发震断层的WNW方向和ESE方向分布有托莱山、大通山、达坂山和冷龙岭等山脉（图1），属于山地地貌，因此该区域的地震动高值可能与山脊地区地震波多次反射有关。
3. 讨论与结论
强地面运动模拟结果由地形数据、介质模型和震源模型共同决定，因此数据选择对结果的可靠性有较大影响，为使结果更为可信，本文参考美国加州综合地震破裂预测模型（Field et al，2014）中提出的“使用可获得最优解”原则对介质模型进行选择。常用的介质模型包括CRUST1.0模型（Laske et al，2012）和CRUST2.0模型（Bassin et al，2000），二者均基于水平层状介质的假设，水平分辨率分别为1.0°和2.0°，而Han等（2022）提出的USTClitho2.0模型的水平分辨率更高，因此研究中选择USTClitho2.0模型作为介质模型。本文基于门源地震震源区附近的地形数据、介质模型和震源模型，使用曲线网格有限差分方法计算了该地震的近场地震波传播过程，得到了理论的速度场快照和地震烈度分布。研究结果表明：总体上垂直断层走向方向的地震动衰减大于平行断层走向方向的地震动衰减，震中最大烈度为Ⅷ度；模拟得到的理论烈度同野外调查的地震烈度分布基本一致。受强地面运动方向性效应和起伏地表的影响，地震灾害主要沿断层的WNW方向和ESE方向分布。由于地震发生区域主要以山地地貌为主，该地震的发生造成了局部边坡崩塌、滚石和冻土开裂以及兰新高铁大桥桥面受损等次生灾害（颉满斌，2022）。考虑到研究区域仍有发生强震的危险（许英才等，2022），今后的防震减灾工作中有必要加强地表起伏剧烈区域特别是山脊地貌的震害防御工作。
北京大学张勇教授为本文提供了震源破裂过程的初步结果，中国科学技术大学张海江教授为本文提供了三维速度结构，中国地震局工程力学研究所马强研究员为本文提供了仪器观测的地震烈度分布图，中国地震局地震预测研究所徐岳仁研究员与作者就青藏高原东北缘的地质概况进行了讨论，审稿专家为本文提出了宝贵意见，作者在此一并表示感谢。
张勇2022年1月10日与作者的个人交流张勇2022年1月10日与作者的个人交流 -
![]()
图 1 研究区及周边的地质概况和2022年1月8—20日震源区MS≥3.0 地震分布
Figure 1. Geological overview of the study area and its surrounding regions，spatial distribution of MS≥3.0 earthquakes in source zone during the period of 8 to 20 January，2022
![]()
图 3 本文所用的门源MS6.9地震震源模型
Figure 3. Source model of Menyuan MS6.9 earthquake used in this study
![]()
图 4 门源MS6.9地震x分量粒子速度波场快照
Figure 4. Wavefield snapshots of the x component particle velocity of Menyuan MS6.9 earthquake
![]()
图 6 门源MS6.9地震仪器观测的烈度分布（引自中国地震局工程力学研究所强震动观测组，2022）
Figure 6. Instrumental seismic intensity distribution of Menyuan MS6.9 earthquake （after Strong Motion Observation Group，Institute of Engineering Mechanics，China Earthquake Administration，2022）
-
国家市场监督管理总局, 国家标准化管理委员会. 2021. GB/T 17742—2020中国地震烈度表[S]. 北京: 中国标准出版社: 1–11. State Administration for Market Regulatory, Standardization Administration. 2021. GB/T 17742−2020 China Seismic Intensity Scale[S]. Beijing: China Standard Press: 1–11 （in Chinese）.
胡朝忠,杨攀新,李智敏,黄帅堂,赵妍,陈丹,熊仁伟,陈庆宇. 2016. 2016年1月21日青海门源6.4级地震的发震机制探讨[J]. 地球物理学报,59(5):1637–1646. doi: 10.6038/cjg20160509 Hu C Z,Yang P X,Li Z M,Huang S T,Zhao Y,Chen D,Xiong R W,Chen Q Y. 2016. Seismogenic mechanism of the 21 January 2016 Menyuan,Qinghai MS6.4 earthquake[J]. Chinese Journal of Geophysics,59(5):1637–1646 (in Chinese).
姜文亮. 2018. 冷龙岭断裂带全新世破裂模式、大震复发特征研究及其区域构造意义[D]. 北京: 中国地震局地质研究所: 1–168. Jiang W L. 2018. Holocene Rupture Pattern, Seismic Recurrence Feature of the Lenglongling Fault Zone and Its Tectonic Implication for the Northeast Tibetan Plateau[D]. Beijing: Institute of Geology, China Earthquake Administration: 1–168 （in Chinese）.
颉满斌. 2022. 青海门源6.9级地震现场考察结果发布[EB/OL]. [2022-01-18]. http://www.stdaily.com/index/kejixinwen/202201/8b2a970831014fb4b6e9573276661278.shtml. Jie M B. 2022. Qinghai Menyuan MS6.9 earthquake site inspection results released[EB/OL]. [2022-01-18] .http://www.stdaily.com/index/kejixinwen/202201/8b2a970831014fb4b6e9573276661278.shtml （in Chinese）.
李强,江在森,武艳强,赵静,魏文薪,刘晓霞. 2013. 海原—六盘山断裂带现今构造变形特征[J]. 大地测量与地球动力学,33(2):18–22. Li Q,Jiang Z S,Wu Y Q,Zhao J,Wei W X,Liu X X. 2013. Present-day tectonic deformation characteristics of Haiyuan-Liupanshan fault zone[J]. Journal of Geodesy and Geodynamics,33(2):18–22 (in Chinese).
青海日报. 2022. 门源县抗震救灾工作基本结束[N/OL]. [2022-01-13]. https://epaper.tibet3.com/qhrb/html/202201/13/content_80039.html. Qinghai Daily. 2022. Earthquake relief work in Menyuan country is basically completed[N/OL]. [2022-01-13]. https://epaper.tibet3.com/qhrb/html/202201/13/content_80039.html （in Chinese）.
青海省地震局. 2022. 青海门源“1·8”6.9级地震发现22公里地表破裂[EB/OL]. [2022-01-08]. https://www.qhdzj.gov.cn/Item/2/20969.aspx. Qinghai Earthquake Agency. 2022. A 22 km surface rupture zone was discovered after “1·8” Qinghai Menyuan MS6.9 earthquake[EB/OL]. [2022-01-08]. https://www.qhdzj.gov.cn/Item/2/20969.aspx （in Chinese）.
徐剑侠,张振国,戴文杰,张伟,Akram N,文健,陈晓非. 2015. 2015年4月25日尼泊尔地震波场传播及烈度初步模拟分析[J]. 地球物理学报,58(5):1812–1817. doi: 10.6038/cjg20150531 Xu J X,Zhang Z G,Dai W J,Zhang W,Akram N,Wen J,Chen X F. 2015. Preliminary simulation of seismic wave propagation and the intensity map for the 25 April 2015 Nepal earthquake[J]. Chinese Journal of Geophysics,58(5):1812–1817 (in Chinese).
许英才,郭祥云,冯丽丽. 2022. 2022年1月8日青海门源MS6.9地震序列重定位和震源机制研究[J]. 地震学报,44(2):195–210. doi: 10.11939/jass.20210139 Xu Y C,Guo X Y,Feng L L. 2022. Relocation and focal mechanism solutions of the MS6.9 Menyuan earthquake sequence on January 8,2022 in Qinghai Province[J]. Acta Seismologica Sinica,44(2):195–210 (in Chinese).
央视网. 2022. 受青海门源地震影响兰新高铁列车停运[EB/OL]. [2022-01-08]. http://news.cctv.com/2022/01/08/ARTI3ZvEZmlKHWl62G7sSGqm220108.shtml. CCTV. 2022. Affected by Qinghai Menyuan earthquake, Lanxin high-speed trains are suspended[EB/OL]. [2022-01-08]. http://news.cctv.com/2022/01/08/ARTI3ZvEZmlKHWl62G7sSGqm220108.shtml （in Chinese）.
赵宏阳,陈晓非. 2017. 1975年海城MS7.3地震强地面运动模拟[J]. 地球物理学报,60(7):2707–2715. doi: 10.6038/cjg20170717 Zhao H Y,Chen X F. 2017. Simulation of strong ground motion by the 1975 Haicheng MS7.3 earthquake[J]. Chinese Journal of Geophysics,60(7):2707–2715 (in Chinese).
张振国,张伟,孙耀充,朱耿尚,文健,陈晓非. 2014a. 2014年2月12日新疆于田地震强地面运动初步模拟及烈度预测[J]. 地球物理学报,57(2):685–689. Zhang Z G,Zhang W,Sun Y C,Zhu G S,Wen J,Chen X F. 2014a. Preliminary simulation of strong ground motion for Yutian,Xinjiang earthquake of 12 February 2014,and hazard implication[J]. Chinese Journal of Geophysics,57(2):685–689 (in Chinese).
张振国,孙耀充,徐建宽,张伟,陈晓非. 2014b. 2014年8月3日云南鲁甸地震强地面运动初步模拟及烈度预测[J]. 地球物理学报,57(9):3038–3041. Zhang Z G,Sun Y C,Xu J K,Zhang W,Chen X. 2014b. Preliminary simulation of strong ground motion for Ludian,Yunnan earthquake of 3 August 2014,and hazard implication[J]. Chinese Journal of Geophysics,57(9):3038–3041 (in Chinese).
中国地震局. 2022. 中国地震局发布青海门源6.9级地震烈度图[EB/OL]. [2022-01-11].https://www.cea.gov.cn/cea/xwzx/fzjzyw/5646200/index.html. China Earthquake Administration. 2022. China Earthquake Administration released the intensity map of Qinghai Menyuan MS6.9 earthquake[EB/OL]. [2022-01-11]. https://www.cea.gov.cn/cea/xwzx/fzjzyw/5646200/index.html （in Chinese）.
中国地震局工程力学研究所强震动观测组.2022.2022年01月08日青海门源6.9级地震仪器烈度分布图[EB/OL]. [2022-01-20]. https://mp.weixin.qq.com/s/IGp1dw7KfFY3PiwWT2WY6A. Strong Motion Observation Group, Institute of Engineering Mechanics, China Earthquake Administration. 2022. Instrument intensity distribution map of the Menyuan, Qinghai MS6.9 earthquake on January 8th, 2022[EB/OL]. [2022-01-20]. https://mp.weixin.qq.com/s/IGp1dw7KfFY3PiwWT2WY6A （in Chinese）.
Bassin C,Laske G,Masters G. 2000. The current limits of resolution for surface wave tomography in North America[J]. EOS Trans AGU,81:F897.
Brocher T M. 2005. Empirical relations between elastic wavespeeds and density in the Earth’s crust[J]. Bull Seismol Soc Am,95(6):2081–2092. doi: 10.1785/0120050077
Field E H,Arrowsmith R J,Biasi G P,Bird P,Dawson T E,Felzer K P,Jackson D D,Johnson K M,Jordan T H,Madden C,Michael A J,Milner K R,Page M T,Parsons T,Powers P M,Shaw B E,Thatcher W R,Weldon R J,Zeng Y H. 2014. Uniform California earthquake rupture forecast,version 3 (UCERF3):The time-independent model[J]. Bull Seismol Soc Am,104(3):1122–1180. doi: 10.1785/0120130164
Han S C,Zhang H J,Xin H L,Shen W S,Yao H J. 2022. USTClitho2.0:Updated unified seismic tomography models for continental China lithosphere from joint inversion of body-wave arrival times and surface-wave dispersion data[J]. Seismol Res Lett,93(1):201–215. doi: 10.1785/0220210122
Laske G, Masters G, Ma Z, Pasyanos M E, 2012. CRUST1.0: An Updated Global Model of Earth’s Crust[C]//EGU General Assembly Conference Abstracts. Vienna: EGUGA: 3743.
Zhang W,Chen X F. 2006. Traction image method for irregular free surface boundaries in finite difference seismic wave simulation[J]. Geophys J Int,167(1):337–353. doi: 10.1111/j.1365-246X.2006.03113.x
Zhang W,Shen Y,Chen X F. 2008. Numerical simulation of strong ground motion for the MS8.0 Wenchuan earthquake of 12 May 2008[J]. Science in China:Series D,51(12):1673–1682. doi: 10.1007/s11430-008-0130-4
Zhang W,Zhang Z G,Chen X F. 2012. Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated‐grid finite‐difference method on curvilinear grids[J]. Geophys J Int,190(1):358–378. doi: 10.1111/j.1365-246X.2012.05472.x
-
期刊类型引用(4)
1. 谢张迪，于湘伟，章文波. 2022年青海门源M_W6.6地震震源动力学破裂过程. 地震学报. 2025(01): 21-36 . 
本站查看
2. 罗超，曹晓雨，高阳，徐飞，冯怀平，王昊. GP断层破裂模型在青海门源M_S6.9地震强地面运动模拟的适用性研究. 地球物理学报. 2024(02): 534-547 . 百度学术
3. 王子博，刘瑞丰，李赞，胡岩松，孔韩东. 2022年1月8日青海门源M_S6.9地震的静态和动态震源参数测定. 地球物理学报. 2023(06): 2420-2430 . 百度学术
4. 张海，高俊焱，王岱，刘中宪，李程程，张聪. 青海门源M_S6.9地震强地面运动精细化模拟. 震灾防御技术. 2023(03): 462-470 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 1091
- HTML全文浏览量: 291
- PDF下载量: 257
- 被引次数: 7
