留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

余鹏飞 陈威 乔学军 赵斌 李刚 熊维

余鹏飞, 陈威, 乔学军, 赵斌, 李刚, 熊维. 基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
引用本文: 余鹏飞, 陈威, 乔学军, 赵斌, 李刚, 熊维. 基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
YU Pengfei, CHEN Wei, QIAO Xuejun, ZHAO Bin, LI Gang, XIONG Wei. Slip Model of the 2022 Menyuan Ms 6.9 Earthquake Constrained by Mulit-Source SAR Data[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
Citation: YU Pengfei, CHEN Wei, QIAO Xuejun, ZHAO Bin, LI Gang, XIONG Wei. Slip Model of the 2022 Menyuan Ms 6.9 Earthquake Constrained by Mulit-Source SAR Data[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

doi: 10.13203/j.whugis20220093
基金项目: 

国家自然科学基金 42130101

国家自然科学基金 41731071

国家自然科学基金 41874018

国家自然科学基金 42074015

中国地震局地震研究所基本科研业务费专项 IS201956148

详细信息
    作者简介:

    余鹏飞,硕士,助理研究员,主要从事大地测量与地球动力学研究。yupengfei@cgps.ac.cn

    通讯作者: 熊维,硕士,助理研究员。xiongwei_19881229@163.com
  • 中图分类号: P315

Slip Model of the 2022 Menyuan Ms 6.9 Earthquake Constrained by Mulit-Source SAR Data

Funds: 

The National Natural Science Foundation of China 42130101

The National Natural Science Foundation of China 41731071

The National Natural Science Foundation of China 41874018

The National Natural Science Foundation of China 42074015

the Scientific Research Fund of Institute of Seismology, China Earthquake Administration IS201956148

More Information
  • 摘要: 2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达(synthetic aperture radar,SAR)遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量(differential interferometry synthetic aperture radar,D-InSAR)技术获取了门源地震的同震形变场,视线(line of sight,LOS)向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。
  • 图  1  2022年门源地震构造背景

    Figure  1.  Tectonic Setting of the 2022 Menyuan Earthquake

    图  2  2022年门源地震同震形变场

    Figure  2.  Coseismic Deformation Field of the 2022 Menyuan Earthquake

    图  3  Sentinel-2影像东西向形变场

    Figure  3.  Coseismic Deformation Field in EW Direction from Sentinel-2 Images

    图  4  InSAR反演的2022年门源地震滑动分布

    Figure  4.  Slip Distribution of the 2022 Menyuan Earthquake Inverted by InSAR

    图  5  InSAR同震形变场的观测值、模拟值及残差分布

    Figure  5.  Observed and Modelled InSAR Coseismic Deformation and Their Residuals

    图  6  2022年门源地震引起的静态库伦应力变化

    Figure  6.  Static Coulomb Stress Changes Induced by the 2022 Menyuan Earthquake

    表  1  2022年门源地震震源参数

    Table  1.   Source Parameters for the 2022 Menyuan Earthquake

    机构 震中位置 深度/km 节面Ⅰ 节面Ⅱ 震级(Mw)
    走向/(°) 倾角/(°) 滑动角/(°) 走向/(°) 倾角/(°) 滑动角/(°)
    USGS 101.278°E,37.815°N 11.5 104 88 15 13 75 178 6.6
    GCMT 101.31°E,37.80°N 14.8 104 82 1 14 89 172 6.7
    下载: 导出CSV

    表  2  SAR影像数据信息

    Table  2.   Parameters of SAR Images

    轨道 飞行方向 获取时间 时间基线/d 垂直基线/m
    主影像 辅影像
    26 升轨 2021-12-29 2022-01-10 12 -110.34
    33 降轨 2021-12-29 2022-01-10 12 56.24
    128 升轨 2022-01-05 2022-01-17 12 38.34
    下载: 导出CSV
  • [1] 李智敏, 盖海龙, 李鑫, 等. 2022年青海门源Ms 6.9级地震发震构造和地表破裂初步调查[J]. 地质学报, 2022, 96(1): 330-335 https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202201023.htm

    Li Zhimin, Gai Hailong, Li Xin, et al. Seismogenic Fault and Coseismic Surface Deformation of the Menyuan Ms 6.9 Earthquake in Qinghai, China[J]. Acta Geologica Sinica, 2022, 96(1): 330-335 https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202201023.htm
    [2] 张培震, 邓起东, 张国民, 等. 中国大陆的强震活动与活动地块[J]. 中国科学(D辑: 地球科学), 2003, 33(S1): 12-20 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1001.htm

    Zhang Peizhen, Deng Qidong, Zhang Guomin, et al. Active Tectonic Blocks and Strong Earthquakes in the Continent of China[J]. Science in China(Ser. D), 2003, 33(S1): 12-20 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1001.htm
    [3] 徐纪人, 姚立珣, 汪进. 1986年8月26日门源6.4级地震及其强余震的震源机制解[J]. 西北地震学报, 1986, 8(4): 82-84 https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ198604012.htm

    Xu Jiren, Yao Lixun, Wang Jin. Earthquake Source Mechanisms of Menyuan Earthquake (Ms=6.4, on Aug. 26, 1986)and Its Strong Aftershocks[J]. Northwestern Seismological Journal, 1986, 8(4): 82-84 https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ198604012.htm
    [4] Zhang W, Jiao D, Zhang P, et al. Displacement Along the Haiyuan Fault Associated with the Great 1920 Haiyuan, China, Earthquake[J]. Bulletin of the Seismological Society of America, 1987, 77(1): 117-131
    [5] Gaudemer Y, Tapponnier P, Meyer B, et al. Partitioning of Crustal Slip Between Linked, Active Faults in the Eastern Qilian Shan, and Evidence for a Major Seismic Gap, the 'Tianzhu Gap', on the Western Haiyuan Fault, Gansu (China)[J]. Geophysical Journal International, 1995, 120(3): 599-645 doi:  10.1111/j.1365-246X.1995.tb01842.x
    [6] Xu X, Yeats R S, Yu G. Five Short Historical Earthquake Surface Ruptures near the Silk Road, Gansu Province, China[J]. Bulletin of the Seismological Society of America, 2010, 100(2): 541-561 doi:  10.1785/0120080282
    [7] Guo P, Han Z J, Mao Z B, et al. Paleoearthquakes and Rupture Behavior of the Lenglongling Fault: Implications for Seismic Hazards of the Northeastern Margin of the Tibetan Plateau[J]. Journal of Geophysical Research: Solid Earth, 2019, 124(2): 1520-1543 doi:  10.1029/2018JB016586
    [8] Lasserre C, Morel P H, Gaudemer Y, et al. Postglacial Left Slip Rate and Past Occurrence of M ≥ 8 Earthquakes on the Western Haiyuan Fault, Gansu, China[J]. Journal of Geophysical Research: Solid Earth, 1999, 104(B8): 17633-17651 doi:  10.1029/1998JB900082
    [9] Molnar P, Tapponnier P. Cenozoic Tectonics of Asia: Effects of a Continental Collision: Features of Recent Continental Tectonics in Asia can Be Interpreted as Results of the India-Eurasia Collision[J]. Science, 1975, 189(4201): 419-426 doi:  10.1126/science.189.4201.419
    [10] Clark M K. Continental Collision Slowing Due to Viscous Mantle Lithosphere rather than Topography[J]. Nature, 2012, 483(7387): 74-77 doi:  10.1038/nature10848
    [11] Tapponnier P, Xu Z Q, Roger F, et al. Oblique Stepwise Rise and Growth of the Tibet Plateau[J]. Science, 2001, 294(5547): 1671-1677 doi:  10.1126/science.105978
    [12] Zheng W J, Zhang P Z, He W G, et al. Transformation of Displacement Between Strike-Slip and Crustal Shortening in the Northern Margin of the Tibetan Plateau: Evidence from Decadal GPS Measurements and Late Quaternary Slip Rates on Faults[J]. Tectonophysics, 2013, 584: 267-280 doi:  10.1016/j.tecto.2012.01.006
    [13] 郭鹏. 北祁连山冷龙岭断裂大震复发行为与危险性研究[D]. 北京: 中国地震局地质研究所, 2019

    Guo Peng. Earthquake Recurrence Behavior and Seismic Hazards of the Lenglongling Fault, Northern Qilian Shan[D]. Beijing: Institute of Geology, China Earthquake Administration, 2019
    [14] Peltzer G, Tapponnier P. Formation and Evolution of Strike-Slip Faults, Rifts, and Basins During the India-Asia Collision: An Experimental Approach[J]. Journal of Geophysical Research: Solid Earth, 1988, 93(B12): 15085-15117 doi:  10.1029/JB093iB12p15085
    [15] 何文贵, 袁道阳, 葛伟鹏, 等. 祁连山活动断裂带中东段冷龙岭断裂滑动速率的精确厘定[J]. 地震, 2010, 30(1): 131-137 doi:  10.3969/j.issn.1000-3274.2010.01.015

    He Wengui, Yuan Daoyang, Ge Weipeng, et al. Determination of the Slip Rate of the Lenglongling Fault in the Middle and Eastern Segments of the Qilian Mountain Active Fault Zone[J]. Earthquake, 2010, 30(1): 131-137 doi:  10.3969/j.issn.1000-3274.2010.01.015
    [16] Jiang W L, Han Z J, Guo P, et al. Slip Rate and Recurrence Intervals of the East Lenglongling Fault Constrained by Morphotectonics: Tectonic Implications for the Northeastern Tibetan Plateau[J]. Lithosphere, 2017, 9(3): 417-430 doi:  10.1130/L597.1
    [17] 郭鹏, 韩竹军, 姜文亮, 等. 青藏高原东北缘冷龙岭断裂全新世左旋滑动速率[J]. 地震地质, 2017, 39(2): 323-341 doi:  10.3969/j.issn.0253-4967.2017.02.005

    Guo Peng, Han Zhujun, Jiang Wenliang, et al. Holocene Left-Lateral Slip Rate of the Lenglongling Fault, Northeastern Margin of the Tibetan Plateau[J]. Seismology and Geology, 2017, 39(2): 323-341 doi:  10.3969/j.issn.0253-4967.2017.02.005
    [18] Lasserre C, Gaudemer Y, Tapponnier P, et al. Fast Late Pleistocene Slip Rate on the Leng Long Ling Segment of the Haiyuan Fault, Qinghai, China[J]. Journal of Geophysical Research: Solid Earth, 2002, 107(B11): ETG4-1
    [19] Zhang P Z, Shen Z K, Wang M, et al. Continuous Deformation of the Tibetan Plateau from Global Positioning System Data[J]. Geology, 2004, 32(9): 809 doi:  10.1130/G20554.1
    [20] Daout S, Jolivet R, Lasserre C, et al. Along-Strike Variations of the Partitioning of Convergence Across the Haiyuan Fault System Detected by InSAR[J]. Geophysical Journal International, 2016, 205(1): 536-547 doi:  10.1093/gji/ggw028
    [21] Cavalié O, Lasserre C, Doin M P, et al. Measurement of Interseismic Strain Across the Haiyuan Fault (Gansu, China), by InSAR[J]. Earth and Planetary Science Letters, 2008, 275(3/4): 246-257
    [22] Jolivet R, Lasserre C, Doin M P, et al. Shallow Creep on the Haiyuan Fault (Gansu, China) Revealed by SAR Interferometry[J]. Journal of Geophysical Research: Solid Earth, 2012, 117: B06401
    [23] Zheng G, Wang H, Wright T J, et al. Crustal Deformation in the India-Eurasia Collision Zone from 25 Years of GPS Measurements[J]. Journal of Geophysical Research: Solid Earth, 2017, 122(11): 9290-9312 doi:  10.1002/2017JB014465
    [24] Hao M, Wang Q L, Shen Z K, et al. Present Day Crustal Vertical Movement Inferred from Precise Leveling Data in Eastern Margin of Tibetan Plateau[J]. Tectonophysics, 2014, 632: 281-292 doi:  10.1016/j.tecto.2014.06.016
    [25] Werner C, Wegmüller U, Strozzi T, et al. GAMMA SAR and Interferometric Processing Software[C]//ERS-ENVISAT Symposium, Gothenburg, Sweden, 2000
    [26] Goldstein R M, Werner C L. Radar Interferogram Filtering for Geophysical Applications[J]. Geophysical Research Letters, 1998, 25(21): 4035-4038 doi:  10.1029/1998GL900033
    [27] Pepe A, Lanari R. On the Extension of the Minimum Cost Flow Algorithm for Phase Unwrapping of Multitemporal Differential SAR Interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(9): 2374-2383 doi:  10.1109/TGRS.2006.873207
    [28] Jonsson S. Fault Slip Distribution of the 1999 Mw 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements[J]. Bulletin of the Seismological Society of America, 2002, 92(4): 1377-1389 doi:  10.1785/0120000922
    [29] 谭悦. 南宁市地表形变InSAR监测及时空分析[J]. 测绘地理信息, 2021, 46(6): 46-50 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202106010.htm

    Tan Yue. Surface Deformation Monitoring and Temporal Spatial Analysis in Nanning City with InSAR[J]. Journal of Geomatics, 2021, 46(6): 46-50 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202106010.htm
    [30] 王乐洋, 邹阿健, 许光煜. 利用光学影像相关获取2019年Ridgecrest地震序列同震形变[J]. 测绘工程, 2021, 30(4): 1-8 https://www.cnki.com.cn/Article/CJFDTOTAL-CHGC202104001.htm

    Wang Leyang, Zou Ajian, Xu Guangyu. Coseismic Deformation of 2019 Ridgecrest Earthquake Sequence Obtained by Optical Images Correlation[J]. Engineering of Surveying and Mapping, 2021, 30(4): 1-8 https://www.cnki.com.cn/Article/CJFDTOTAL-CHGC202104001.htm
    [31] Leprince S, Berthier E, Ayoub F, et al. Monitoring Earth Surface Dynamics with Optical Imagery[J]. Eos, Transactions American Geophysical Union, 2008, 89(1): 1-2 doi:  10.1029/2008EO010001
    [32] Leprince S, Barbot S, Ayoub F, et al. Automatic and Precise Orthorectification, Coregistration, and Subpixel Correlation of Satellite Images, Application to Ground Deformation Measurements[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(6): 1529-1558 doi:  10.1109/TGRS.2006.888937
    [33] Fan L, Li B, Liao S, et al. Precise Relocation of the Aftershock Sequences of the 2022 M 6.9 Menyuan Earthquake[J]. Earthquake Science, 2022, 35(3): Q20220008
    [34] Wang R, Diao F, Hoechner A. SDM-A Geodetic Inversion Code Incorporating with Layered Crust Structure and Curved Fault Geometry[C]//EGU General Assembly 2013, Vienna, Austria, 2013
    [35] Yang H F, Wang D, Guo R M, et al. Rapid Report of the 8 January 2022 ​Ms 6.9 Menyuan Earthquake, Qinghai, China[J]. Earthquake Research Advances, 2022, 2(1): 100113 doi:  10.1016/j.eqrea.2022.100113
    [36] 李振洪, 韩炳权, 刘振江, 等. InSAR数据约束下的2016年和2022年青海门源地震震源参数及其滑动分布[J]. 武汉大学学报·信息科学版, 2022, DOI: 10.13203/j.whugis20220037

    Li Zhenhong, Han Bingquan, Liu Zhenjiang, et al. Source Parameters and Slip Distributions of the 2016 and 2022 Menyuan, Qinghai Earthquakes Constrained by InSAR Observations[J]. Geomatics and Information Science of Wuhan University, 2022, DOI: 10.13203/j.whugis20220037
    [37] Wang R J, Lorenzo-Martín F, Roth F. PSGRN/PSCMP—A New Code for Calculating Co- and Post-Seismic Deformation, Geoid and Gravity Changes Based on the Viscoelastic-Gravitational Dislocation Theory[J]. Computers & Geosciences, 2006, 32(4): 527-541
    [38] 石富强, 邵志刚, 占伟, 等. 青藏高原东北缘活动断裂剪切模量及应力状态数值模拟[J]. 地球物理学报, 2018, 61(9): 3651-3663 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201809014.htm

    Shi Fuqiang, Shao Zhigang, Zhan Wei, et al. Numerical Modeling of the Shear Modulus and Stress State of Active Faults in the Northeastern Margin of the Tibetan Plateau[J]. Chinese Journal of Geophysics, 2018, 61(9): 3651-3663 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201809014.htm
  • [1] 杜玉玲, 闫世勇, 张豪磊, 赵峰, 仇春平.  . 武汉大学学报 ● 信息科学版, 2022, 0(0): 0-0. doi: 10.13203/j.whugis20210365
    [2] 柴海山, 陈克杰, 魏国光, 方荣新, 邹蓉, 祝会忠.  北斗三号与超高频GNSS同震形变监测:以2021年青海玛多Mw 7.4地震为例 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 946-954. doi: 10.13203/j.whugis20220140
    [3] 李振洪, 韩炳权, 刘振江, 张苗苗, 余琛, 陈博, 刘海辉, 杜静, 张双成, 朱武, 张勤, 彭建兵.  InSAR数据约束下2016年和2022年青海门源地震震源参数及其滑动分布 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 887-897. doi: 10.13203/j.whugis20220037
    [4] 李志才, 丁开华, 张鹏, 温扬茂, 赵利江, 陈建峰.  GNSS观测的2021年青海玛多地震(Mw 7.4)同震形变及其滑动分布 . 武汉大学学报 ● 信息科学版, 2021, 46(10): 1489-1497. doi: 10.13203/j.whugis20210301
    [5] 赵超英, 刘晓杰, 张勤, 彭建兵, 许强.  甘肃黑方台黄土滑坡InSAR识别、监测与失稳模式研究 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 996-1007. doi: 10.13203/j.whugis20190072
    [6] 王乐洋, 赵雄, 高华.  大地测量地震断层同震滑动分布反演的两步解法 . 武汉大学学报 ● 信息科学版, 2019, 44(9): 1265-1273, 1311. doi: 10.13203/j.whugis20170382
    [7] 吴文豪, 李陶, 龙四春, 周志伟.  实时轨道条件下Sentinel-1卫星影像干涉配准 . 武汉大学学报 ● 信息科学版, 2019, 44(5): 745-750. doi: 10.13203/j.whugis20170098
    [8] 王乐洋, 李海燕, 陈汉清.  2013年芦山Ms 7.0级地震断层参数模型反演 . 武汉大学学报 ● 信息科学版, 2019, 44(3): 347-354. doi: 10.13203/j.whugis20160485
    [9] 刘洋, 许才军, 温扬茂.  门源Mw5.9级地震形变InSAR观测及区域断裂带深部几何形态 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 1035-1042. doi: 10.13203/j.whugis20190069
    [10] 苏利娜, 张勇.  GPS时间序列中同震和震后形变的自动识别和估计 . 武汉大学学报 ● 信息科学版, 2018, 43(10): 1504-1510. doi: 10.13203/j.whugis20170016
    [11] 熊新, 靳国旺, 张红敏, 王岩.  基于史赖伯规则的机载InSAR区域网平差 . 武汉大学学报 ● 信息科学版, 2017, 42(9): 1292-1299, 1305. doi: 10.13203/j.whugis20150803
    [12] 许才军, 邓长勇, 周力璇.  利用方差分量估计的地震同震滑动分布反演 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 37-44. doi: 10.13203/j.whugis20150500
    [13] 李志才, 张鹏, 温扬茂, 廖瑛.  利用GPS和海底基准点观测形变反演日本大地震(Mw 9.0)同震断层滑动分布 . 武汉大学学报 ● 信息科学版, 2013, 38(1): 40-43.
    [14] 许才军, 何平, 温扬茂, 张磊.  日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束 . 武汉大学学报 ● 信息科学版, 2012, 37(12): 1387-1391.
    [15] 张红敏, 靳国旺, 徐青, 秦志远.  利用差分滤波进行多基线InSAR相位解缠 . 武汉大学学报 ● 信息科学版, 2011, 36(9): 1030-1034.
    [16] 许才军, 林敦灵, 温扬茂.  利用InSAR数据的汶川地震形变场提取及分析 . 武汉大学学报 ● 信息科学版, 2010, 35(10): 1138-1142.
    [17] 刁法启, 熊熊, 李军, 郑勇.  利用GPS观测反演Tokachi-Oki地震(日本)及强余震的同震滑移分布 . 武汉大学学报 ● 信息科学版, 2010, 35(3): 270-273.
    [18] 温扬茂, 许才军.  基于敏感度的迭代拟合法反演玛尼Ms7.9级地震滑动分布 . 武汉大学学报 ● 信息科学版, 2009, 34(6): 732-735.
    [19] 周春霞, 鄂栋臣, 廖明生.  InSAR用于南极测图的可行性研究 . 武汉大学学报 ● 信息科学版, 2004, 29(7): 619-623.
    [20] 韩炳权, 刘振江, 陈博, 李振洪, 余琛, 张勇, 彭建兵.  . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20220636
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  358
  • HTML全文浏览量:  131
  • PDF下载量:  82
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-04-29
  • 刊出日期:  2022-06-05

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

doi: 10.13203/j.whugis20220093
    基金项目:

    国家自然科学基金 42130101

    国家自然科学基金 41731071

    国家自然科学基金 41874018

    国家自然科学基金 42074015

    中国地震局地震研究所基本科研业务费专项 IS201956148

    作者简介:

    余鹏飞,硕士,助理研究员,主要从事大地测量与地球动力学研究。yupengfei@cgps.ac.cn

    通讯作者: 熊维,硕士,助理研究员。xiongwei_19881229@163.com
  • 中图分类号: P315

摘要: 2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达(synthetic aperture radar,SAR)遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量(differential interferometry synthetic aperture radar,D-InSAR)技术获取了门源地震的同震形变场,视线(line of sight,LOS)向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。

English Abstract

余鹏飞, 陈威, 乔学军, 赵斌, 李刚, 熊维. 基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
引用本文: 余鹏飞, 陈威, 乔学军, 赵斌, 李刚, 熊维. 基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
YU Pengfei, CHEN Wei, QIAO Xuejun, ZHAO Bin, LI Gang, XIONG Wei. Slip Model of the 2022 Menyuan Ms 6.9 Earthquake Constrained by Mulit-Source SAR Data[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
Citation: YU Pengfei, CHEN Wei, QIAO Xuejun, ZHAO Bin, LI Gang, XIONG Wei. Slip Model of the 2022 Menyuan Ms 6.9 Earthquake Constrained by Mulit-Source SAR Data[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 898-906. doi: 10.13203/j.whugis20220093
  • 据中国地震台网测定,北京时间2022-01-08T13:45,在青海省海北州门源县(101.26°E,37.77°N)发生Ms 6.9地震,震源深度10 km。此次地震损坏房屋4 000余间,5 800余人受灾,多处道路、桥梁以及铁路等受损。据野外地质考察,地震造成了长约22 km的地表破裂[1]。地震发生后,多家机构发布的震源机制解均表明此次门源地震是一次走滑型为主的地震事件。

    2022年门源地震震中位于青藏高原东北缘北祁连山冷龙岭地区(图 1),该地区经受着强烈的地壳缩短和左旋剪切,构造变形复杂[2]。此次地震东南方向仅20~35 km距离,曾先后于1986年和2016年发生两次Ms 6.4地震,美国地质调查局(United States Geological Survey,USGS)和全球矩心矩张量项目(global centroid-moment-tensor,GCMT)公布的震源机制解均以逆冲为主(具体参数见表 1),但也有学者认为1986年门源地震有拉张性质[3]。历史地震资料显示该区域曾发生过2次8级以上地震,分别为1920年海原8.5级地震和1927年古浪8.0级地震。1920年海原地震是青藏高原东北缘有历史记载以来最大的地震,地震造成长约220 km的地表破裂,最大同震位移达10 m,并造成22万人死亡[4]。1927年古浪地震震中位于冬青顶山附近,皇城-双塔断裂、古浪断裂及冷龙岭断裂上均有地表破裂产生,发震断裂及震源机制尚存争议[5-7]。海原地震和古浪地震之间存在一个由冷龙岭、金强河、毛毛山和老虎山等断裂组成的长约260 km的天祝地震空区,其全新世构造活动剧烈,8级以上地震复发周期为1 050±450 a,已平静了至少800年,是未来最有可能发生强震的地区[58],此次门源地震即发生在天祝地震空区西端。因此,深入研究此次地震的形变及震源特征,对理解青藏高原东北缘现今应变分配模式及区域地震危险性评估都有重要科学意义。

    图  1  2022年门源地震构造背景

    Figure 1.  Tectonic Setting of the 2022 Menyuan Earthquake

    表 1  2022年门源地震震源参数

    Table 1.  Source Parameters for the 2022 Menyuan Earthquake

    机构 震中位置 深度/km 节面Ⅰ 节面Ⅱ 震级(Mw)
    走向/(°) 倾角/(°) 滑动角/(°) 走向/(°) 倾角/(°) 滑动角/(°)
    USGS 101.278°E,37.815°N 11.5 104 88 15 13 75 178 6.6
    GCMT 101.31°E,37.80°N 14.8 104 82 1 14 89 172 6.7

    合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)技术具有全天时、全天候、覆盖广等优点,其获取的同震地表形变在震源机制等研究中已获得广泛的应用。本文利用震后及时获取的升降轨Sentinel-1观测,获取门源地震的同震形变场,并以此为约束,基于均匀弹性半无限位错模型,反演同震滑移分布和库伦应力变化,并对区域地震危险性进行初步探讨。

    • 新生代以来,印度板块对欧亚大陆的碰撞和持续汇聚作用[9]造就了全世界海拔最高的青藏高原。随着高原地壳的逐渐增厚及北部刚性地块的阻挡[10],高原地壳向东挤出[11],并在青藏高原东北缘发育了阿尔金断裂带、祁连-海原断裂带等巨型走滑断裂带与之相适应。

      图 1所示,青藏高原东北缘是高原向大陆内部扩展的前缘部位。晚第四纪以来,受北侧戈壁-阿拉善和东侧鄂尔多斯块体的阻挡,该地区构造变形十分强烈,主要分解为祁连-海原断裂带的左旋剪切和祁连山褶皱-逆冲带的缩短增厚[12-13]。祁连-海原断裂带由托莱山断裂、冷龙岭断裂、金强河断裂、毛毛山断裂、老虎山断裂和海原断裂组成,并在冷龙岭断裂东端分支出古浪-香山-天景山断裂[13],全长约1 200 km,是青藏高原东北缘主要的边界断裂,以左旋走滑为主,在控制高原东北缘构造格局和调节地壳物质东向挤出运动中起重要作用[14]。此次门源地震即发生在祁连-海原断裂带中部的冷龙岭断裂西段。冷龙岭断裂长约120 km,西端与托莱山断裂呈左阶斜列,东端与古浪断裂和金强河断裂相接,其东西两侧构造变形分别以左旋剪切和褶皱逆冲缩短为主,是青藏高原东北缘的构造转换枢纽部位[13]。冷龙岭断裂晚第四纪滑动速率存在较大争议,范围从4~7 mm/a[1215-17]到15~20 mm/a[518],但一般都认为是祁连-海原断裂带上滑动速率最大的断裂[17-18]

      青藏高原东北缘吸收和调节着印度板块向欧亚板块NNE向约8.7 mm/a的挤压汇聚[19],在此作用下,柴达木地块以12~14 mm/a的速度向NE 61.45°方向运动。祁连地块的速度为7~14 mm/a,优势运动方向为NE 70°~90°[2]。除文献[20]利用InSAR技术估算的沿祁连-海原断裂带的滑动速率相对较高(10 mm/a)外,其他学者[1221-23]利用全球导航卫星系统(global navigation satellite system,GNSS)或InSAR技术估算的滑动速率均在4~6 mm/a左右。1970—2012年间的精密水准测量结果表明,祁连-海原断裂以南区域最大出现5 mm/a左右的隆升,而在断裂以北则迅速衰减至0~1 mm/a左右[24],分界特征明显。

    • 2022年门源地震发生后,欧空局迅速发布了覆盖此次地震的Sentinel-1A(C波段)卫星升降轨影像数据(表 2)。合成孔径雷达(synthetic aperture radar,SAR)影像数据在IW(interferometric wide)观测模式下采用TOPS(tera operations per second)技术进行观测,能够有效降低宽视场成像中出现的扇贝效应,提高影像的干涉性能。干涉像对的时间基线均为12 d,能够尽量减少震后形变的影响。

      表 2  SAR影像数据信息

      Table 2.  Parameters of SAR Images

      轨道 飞行方向 获取时间 时间基线/d 垂直基线/m
      主影像 辅影像
      26 升轨 2021-12-29 2022-01-10 12 -110.34
      33 降轨 2021-12-29 2022-01-10 12 56.24
      128 升轨 2022-01-05 2022-01-17 12 38.34

      本文采用二轨法和GAMMA软件对SAR数据进行差分干涉处理[25],采用分辨率为1″的SRTM(shuttle radar topography mission)数字高程模型用于地形相位的消除和地理编码,使用加权功率谱技术对干涉图进行滤波处理以抑制噪声[26],采用最小费用流法对干涉图进行相位解缠[27],以多项式拟合消除残余轨道误差的影响,并将雷达坐标系下的相位图通过DEM(digital elevation model)地理编码至WGS-84坐标系,最终获取门源地震视线(line of sight,LOS)向同震形变场,如图 2所示。

      图  2  2022年门源地震同震形变场

      Figure 2.  Coseismic Deformation Field of the 2022 Menyuan Earthquake

      T26升轨SAR影像仅捕获到本次地震西侧的同震形变场,其最大形变为-0.67 m。T128升轨和T33降轨影像则完整捕捉了门源地震的同震形变场。升、降轨干涉图显示门源地震的同震形变场以对称的蝴蝶状分布为特征,对称轴线大致沿113°(或293°)方向,地震造成的地表形变范围约30 km×20 km。升降轨干涉图条纹总体连续光滑,但靠近发震断层的位置因严重失相关导致条纹不连续,说明有地表破裂发生,长约20 km;条纹十分密集,说明形变梯度较大。由于升降轨观测的视线向不一致,升轨和降轨图像上发震断层两侧的形变符号相反,量值接近,说明发震断层以水平运动为主。T128升轨最大形变为-0.6 m,T33降轨最大形变为0.75 m。

      InSAR形变场包含数以百万的观测值,将全部观测值用于反演计算不仅影响计算效率,还会使反演结果难以收敛。为此,利用四叉树降采样法[28]对InSAR结果进行降采样处理,在保留原始形变场主要特征的同时大幅降低参与反演的InSAR观测值的数量。

    • 虽然光学影像易受云雾遮挡、太阳照度等因素的影响,其在形变监测领域的应用受到限制,但在气候条件允许情况下,它能很好地克服InSAR技术失相干和最大形变梯度限制等因素的影响,并逐步应用于地震形变监测[29-30]。Sentinel-2卫星携带的多光谱成像仪(multi-spectral imager,MSI)可覆盖13个光谱波段,有4个可见光近红外波段影像(Band 2/3/4/8)空间分辨率高达10 m,且有研究表明Band 8的地表形变监测效果最好[29]

      本文基于光学影像配准和相关算法[31],选取Sentinel-2卫星Band 8波段2022-01-02—2022-01-12和2022-01-05—2022-01-15两组数据,利用亚像素的频率域相关匹配技术对其进行处理,获取门源地震的东西向同震形变场,如图 3所示。由图 3可以看出,东西向最大同震位移达2.5 m;东西向同震形变场有明显的分界,其位置与野外调查的北支地表破裂位置一致[1];发震断裂北侧向西运动,南侧向东运动,表明此次地震具有左旋走滑的运动特征。

      图  3  Sentinel-2影像东西向形变场

      Figure 3.  Coseismic Deformation Field in EW Direction from Sentinel-2 Images

      由于基于亚像素的相关性匹配算法获取同震形变的匹配精度一般为1/20像素[32],即Sentinel-2卫星同震形变监测精度约为0.5 m,远低于由D-InSAR(different-InSAR)技术获取的同震形变,故由Sentinel-2卫星影像获取的同震形变结果未参与断层滑动分布的反演。

    • 2022年门源地震野外考察显示,门源地震的地表破裂带由北支主破裂和南支次破裂组成[1],分别对应于冷龙岭断裂西段和托莱山断裂,余震精定位结果[33]勾勒的断层形态与野外考察的地表破裂相吻合。因此,在反演断层滑动分布时采用双断层模型,断层迹线由野外考察的地表破裂确定。参考余震深度范围将断层宽度设定为20 km。参考USGS和GCMT发布的震源机制解及余震精定位结果[33]将断层倾角设置在70°~90°范围,以1°为步长,通过格网搜索,确定北支断层和南支断层倾角分别为78°和90°时获得最佳拟合效果。

    • 基于匀质半空间弹性位错模型,本文采用基于约束条件的最小二乘及最速下降法进行门源地震同震滑动分布的反演[34]。通过将东西两段发震断层按2 km×2 km的格网离散为20×10和10×10个子断层,以获取发震断层滑动分布的细节。根据L曲线法确定的最佳光滑因子为0.08,获取基于均匀介质模型的门源地震断层滑动分布如图 4所示。结果显示,此次同震破裂主要发生在北支断层,破裂范围约26 km×10 km,集中在0~10 km深度范围,平均滑动量0.41 m,最大滑动量3.25 m,滑动角为10.44°,对应深度为4.89 km。南支断层破裂规模相对较小,破裂范围约12 km×8 km,破裂集中在0~8 km深度范围,最大滑动量1.42 m。由滑动分布导出的矩震量为1.07×1019 N·m,相应的矩震级为6.6,与USGS结果一致。

      图  4  InSAR反演的2022年门源地震滑动分布

      Figure 4.  Slip Distribution of the 2022 Menyuan Earthquake Inverted by InSAR

      图 5为InSAR同震形变场的观测值、模拟值及残差分布。由图 5可以看出,InSAR观测的形变场与模拟值具有较好的一致性,T26升轨、T128升轨和T33降轨的整体残差分别为1.8 cm、2.5 cm和3.2 cm。观测值与模拟值的相关系数为0.9719。

      图  5  InSAR同震形变场的观测值、模拟值及残差分布

      Figure 5.  Observed and Modelled InSAR Coseismic Deformation and Their Residuals

    • 滑动分布模型表明,门源地震同震滑动以左旋走滑为主,破裂主要发生在东南走向的北支断层,南支断层破裂规模相对较小。北支断层破裂范围约26 km×10 km,并造成了长约20 km的地表破裂,与野外考察的结果一致。南支断层破裂范围约12 km×8 km,并破裂至地表。本文反演的滑动分布与文献[35]基本一致,最大滑动量、破裂范围和破裂深度与文献[36]一致,但由于采用的断层模型不同,在破裂的细部特征上存在差异。

      冷龙岭断裂西端出现分叉现象,走向也发生明显变化,由北西西向转变为近东西向[13],断裂与西侧托莱山断裂左阶斜列分布。重定位的余震与冷龙岭断裂及托莱山断裂相吻合,北支主破裂带仅东南段与冷龙岭断裂相吻合,西北段则继续向西北方向延伸而没有随冷龙岭断裂走向变化向西偏转,且西北段基本没有余震发生(见图 1)。利用反投影方法对门源地震破裂过程的成像显示地震发生后破裂向西北方向传播[35]。基于以上分析,初步判定冷龙岭断裂为此次地震的发震断裂。门源地震震中位于冷龙岭断裂走向发生偏转的位置,走向的变化导致形成挤压阶区并积累弹性应变。当积累的弹性应变超过岩石强度时发生破裂,地震发生。同震破裂沿冷龙岭断裂的分支断裂继续向西北方向传播,应变能得到充分释放,余震较少。破裂传播过程中还引起了与冷龙岭断裂斜接的托莱山断裂发生同震错动。由于走向的变化阻碍了破裂的向西传播,托莱山断裂同震错动规模相对较小,积累的应变能更多以余震的形式释放。

    • 冷龙岭断裂被认为是天祝地震空区的重要组成部分[58]。第四纪滑动速率的研究和现代大地测量观测表明冷龙岭断裂滑动速率为4~7 mm/a,若冷龙岭断裂在1927年古浪地震中发生了全段破裂[7],则1927年以来积累的应变为0.38~0.66 m,与发震断层上的平均滑动量相当,表明在冷龙岭断裂西段积累的应变已得到相当程度的释放,再次发生强震的可能性较小。地震发生后,会在近场和远场引起静态同震应力变化,通过计算和分析同震库仑应力变化可以评价地震对周围断层的影响。基于反演的滑动模型,利用PSGRN/PSCMP软件[37]计算了门源地震引起的同震库伦应力变化(图 6)。结果表明,门源地震使冷龙岭断裂东段和托莱山断裂西段处于应力加载状态,未来仍有发生强震的可能性。此外,冷龙岭断裂东侧的金强河断裂、毛毛山断裂和老虎山断裂晚第四纪活动强烈,800年来未发生过大地震[5],现代GPS观测也表明该区段剪应力率水平较高,预示着这些断裂当前滑动受阻,速率亏损较大[38],是未来有可能发生大地震的地区。

      图  6  2022年门源地震引起的静态库伦应力变化

      Figure 6.  Static Coulomb Stress Changes Induced by the 2022 Menyuan Earthquake

    • 基于欧空局Sentinel-1卫星升降轨数据和D-InSAR技术,获取了2022年1月8日门源Ms 6.9地震的LOS向同震形变场,最大形变约0.75 m。利用Sentinel-2卫星数据及光学影像配准和相关技术获取了本次地震的东西向同震形变场,形变场具有左旋走滑的特征,最大同震位移达2.5 m。基于LOS向同震形变场反演的断层滑动模型显示,门源地震为一次左旋走滑事件,同震位错达26 km,并造成了长约20 km的地表破裂。最大滑动量为3.25 m,对应深度4.89 km。门源地震释放的矩震量为1.07×1019 N·m,对应矩震级6.6。门源地震的发震断层为冷龙岭断裂,并引起托莱山断裂发生同震错动。同震库伦应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力明显升高,未来存在发生强震的可能性。

参考文献 (38)

目录

    /

    返回文章
    返回