2021年5月22日青海玛多MS7.4地震同震变形的数值模拟及其动力学启示

doi:10.7523/j.ucas.2022.046

摘要/Abstract

摘要： 2021年5月22日青海玛多M_S7.4地震是近年来唯一一次发生在巴颜喀拉块体内部的强震，其同震变形及应力调整受到很多学者的高度关注。采用分裂节点技术研发了一套研究大地震同震变形的三维并行弹性有限元程序，并使用弹性地震位错模型的解析解验证该程序的正确性和有效性。在此基础上，研究2021年5月22日青海玛多M_S7.4地震同震变形及研究区域主要断层上的同震库仑应力变化，结果表明：计算得到的地表同震变形与不同大地测量数据和有限断层模型反演得到的结果一致，显示该地震主要为左旋走滑断层，有部分正断层分量；9~11 km深度范围内除主震断层带两侧大于1.0级的余震基本位于10 km深度水平面内平行于主震的库仑应力变化大于零的区域；玛多地震引起不同断裂带上的库仑应力变化的分布是不均匀的，其中鄂拉山断裂南段、昆中断裂东段、东昆仑断裂带的大部分(除中部靠西的小部分外)、玛多—甘德断裂北西段、南东段及主震震中附近、达日断裂带北西段库仑应力增加明显，地震危险性提高，需要密切关注。

关键词: 玛多地震, 同震变形, 有限元, 库仑应力变化, 地震危险性

Abstract: The Maduo M_S7.4 earthquake in Qinghai on May 22, 2021 is the only strong earthquake that occurred in Bayan Har block in recent years. The coseismic deformation and stress adjustment of the Maduo earthquake have been highly concerned by many scholars. In this paper, a three-dimensional parallel elastic finite element program for coseismic deformation simulation of large earthquakes is developed using split node technique. The correctness and effectiveness of the program were verified by comparing with the analytical solutions of the elastic seismic dislocation models. We calculated the coseismic deformation of the Maduo M_S7.4 earthquake and the coseismic Coulomb stress changes on the major faults in the research area. The results show that the surface coseismic deformation is consistent with those by different geodetic observation data and inversion results of some finite fault models. It also shows the earthquake is mainly a sinistral strike-slip fault with some normal fault components. Most aftershocks (magnitude>1.0 and depth of 9-11 km) except those besides the mainshock fault located in the regions of positive Coulomb stress changes at the depth of 10 km. The distribution of the Coulomb stress changes on the different fault zones is not uniform. The Coulomb stress changes on the southern end of Elashan Fault zone, the eastern end of Kunzhong Fault zone, most segments except some central-western parts of East Kunlun Fault zone, the northwestern segment of Dari Fault zone and the northwest and southeast segments and regions near the mainshock of Maduo-Gande Fault zone increase obviously, which requires more attention.

Key words: Maduo earthquake, coseismic deformation, finite element model, Coulomb stress change, seismic hazard

中图分类号:

P315

孟思晨, 孟秋, 陈启志, 胡才博. 2021年5月22日青海玛多M_S7.4地震同震变形的数值模拟及其动力学启示[J]. 中国科学院大学学报, 2024, 41(1): 70-80.

MENG Sichen, MENG Qiu, CHEN Qizhi, HU Caibo. Coseismic deformation of Maduo M_S7.4 earthquake in Qinghai, on May 22, 2021: numerical simulation analysis and geodynamic enlightenment[J]. Journal of University of Chinese Academy of Sciences, 2024, 41(1): 70-80.

参考文献

[1] 邓文泽, 杨志高, 席楠, 等. 2021年5月22日青海玛多M7.4地震的快速测定与数据产品产出[J]. 中国地震, 2021, 37(2): 541-550.DOI:10.3969/j.issn.1001-4683.2021.02.025.
[2] 王未来, 房立华, 吴建平, 等. 2021年青海玛多M_S7.4地震序列精定位研究[J]. 中国科学：地球科学, 2021,51(7): 1193-1202. DOI: 10.1360/SSTe-2021-0149.
[3] 徐志国, 梁姗姗, 张广伟, 等. 2021年5月22日青海玛多M_S7.4地震发震构造分析[J]. 地球物理学报, 2021, 64(8): 2657-2670. DOI:10.6038/cjg2021P0390.
[4] 詹艳, 梁明剑, 孙翔宇, 等. 2021年5月22日青海玛多M_S7.4地震深部环境及发震构造模式[J]. 地球物理学报, 2021, 64(7): 2232-2252.DOI:10.6038/cjg2021O0521.
[5] 李志才, 丁开华, 张鹏, 等. GNSS观测的2021年青海玛多地震(M_W 7.4)同震形变及其滑动分布[J]. 武汉大学学报·信息科学版, 2021, 46(10): 1489-1497.DOI:10.13203/j.whugis20210301.
[6] 尹欣欣, 王维欢, 蔡润, 等. 2021年青海玛多M_S7.4地震精定位和发震构造初探[J]. 地震工程学报, 2021, 43(4): 834-839. DOI:10.3969/j.issn.1000-0844.2021.04.834.
[7] 潘家伟, 白明坤, 李超, 等. 2021年5月22日青海玛多M_S7.4地震地表破裂带及发震构造[J].地质学报, 2021, 95(6): 1655-1670.DOI:10.19762/j.cnki.dizhixuebao.2021166.
[8] 李智敏, 李文巧, 李涛, 等. 2021年5月22日青海玛多M_S7.4地震的发震构造和地表破裂初步调查[J]. 地震地质, 2021, 43(3): 722-737.DOI:10.3969/j.issn.0253-4967.2021.03.016.
[9] 周保, 李五福, 董福辰, 等. 青海省玛多县“5.22 M_S7.4 地震”地表破裂与次生灾害发育特征[J/OL]. 地质通报, (2021-07-27)[2021-08-12]. https://kns.cnki.net/kcms/detail/11.4648.P.20210726.1650.002.html.
[10] 盖海龙, 姚生海, 杨丽萍, 等. 青海玛多“5·22” M_S7.4 级地震的同震地表破裂特征、成因及意义[J]. 地质力学学报, 2021, 27(6): 899-912. DOI:10.12090/j.issn.1006-6616.2021.27.06.073.
[11] 华俊, 赵德政, 单新建, 等. 2021年青海玛多M_W7.3地震InSAR的同震形变场、断层滑动分布及其对周边区域的应力扰动[J]. 地震地质, 2021, 43(3): 677-691. DOI:10.3969/j.issn.0253-4967.2021.03.013.
[12] Liu J H, Hu J, Li Z W, et al. Complete three-dimensional coseismic displacements due to the 2021 Maduo earthquake in Qinghai Province, China from Sentinel-1 and ALOS-2 SAR images[J]. Science China Earth Sciences, 65. DOI:10.1007/s11430-021-9868-9.
[13] Liu C Y, Bai L, Hong S Y, et al. Coseismic deformation of the 2021 M_W7.4 Maduo earthquake from joint inversion of InSAR, GPS, and teleseismic data[J]. Earthquake Science, 2021, 34(5): 1-11. DOI:10.29382/eqs-2021-0050.
[14] Zhao D Z, Qu C Y, Chen H, et al. Tectonic and geometric control on fault kinematics of the 2021 M_W7.3 Maduo (China) earthquake inferred from interseismic, coseismic, and postseismic InSAR observations[J]. Geophysical Research Letters, 2021, 48(18):e2021GL095417. DOI:10.1029/2021GL095417.
[15] 杨君妍, 孙文科, 洪顺英, 等. 2021年青海玛多7.4级地震的同震变形分析[J].地球物理学报, 2021, 64(8): 2671-2683. DOI:10.6038/cjg2021P0416.
[16] 杨光远, 李一帆, 王斌, 等. 青海玛多7.4级地震静态库仑应力分析[J]. 四川地震, 2021, 3: 1-4.DOI:10.13716/j.cnki.1001-8115.2021.03.001.
[17] Toda S, Stein R S, Richards-Dinger K, et al. Forecasting the evolution of seismicity in southern California: animations built on earthquake stress transfer [J]. Journal of Geophysical Research: Solid Earth, 2005, 110(B5): B05S16.DOI:10.1029/2004JB003415.
[18] Lin J, Stein R S. Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults[J]. Journal of Geophysical Research: Solid Earth, 2004, 109(B2): B02303. DOI:10.1029/2003JB002607.
[19] 谈洪波, 王嘉沛, 杨光亮, 等. 2021年玛多M_S7.4地震的震后效应模拟[J]. 地震地质, 2021, 43(4): 936-957. DOI:10.3969/j.issn.0253-4967.2021.04.013.
[20] Li Y J, Huang L Y, Ding R, et al. Coulomb stress changes associated with the M7.3 Maduo earthquake and implications for seismic hazards[J]. Natural Hazards Research, 2021, 1(2): 95-101. DOI:10.1016/j.nhres.2021.06.003.
[21] Wang R J, Martín F L, Roth F. Computation of deformation induced by earthquakes in a multi-layered elastic crust： FORTRAN programs EDGRN/EDCMP[J]. Computers & Geosciences, 2003, 29(2): 195-207.DOI:10.1016/S0098-3004(02)00111-5.
[22] Melosh H J，Raefsky A. A simple and efficient method for introducing faults into finite element computations[J]. Bulletin of the Seismological Society of America, 1981, 71(5): 1391-1400.DOI:10.1785/bssa0710051391.
[23] Okada Y. Internal deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 1992, 82(2): 1018-1040.DOI:10.1785/bssa0820021018.
[24] King G C P, Stein R S, Lin J. Static stress changes and the triggering of earthquakes [J]. Bulletin of the Seismological Society of America, 1994, 84(3): 935-953. DOI:10.1785/bssa0840030935.
[25] 邓起东, 高翔, 陈桂华, 等. 青藏高原昆仑—汶川地震系列与巴颜喀喇断块的最新活动[J]. 地学前缘, 2010, 17(5): 163-178.
[26] 刘雷, 李玉江, 朱良玉, 等. 1947年达日M7.7地震对巴颜喀拉块体边界断裂应力影响的数值模拟[J]. 地球物理学报, 2021, 64(7): 2221-2231. DOI:10.6038/cjg2021P0194.