留言板

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

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

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

张庆云 李永生 罗毅 张景发 杨建钦

张庆云, 李永生, 罗毅, 张景发, 杨建钦. 2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演[J]. 武汉大学学报 ● 信息科学版, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
引用本文: 张庆云, 李永生, 罗毅, 张景发, 杨建钦. 2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演[J]. 武汉大学学报 ● 信息科学版, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
ZHANG Qingyun, LI Yongsheng, LUO Yi, ZHANG Jingfa, YANG Jianqin. Source Parameters for the 2016 Mw 6.2 Italy Amatrice Earthquake Revealed by InSAR Observations[J]. Geomatics and Information Science of Wuhan University, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
Citation: ZHANG Qingyun, LI Yongsheng, LUO Yi, ZHANG Jingfa, YANG Jianqin. Source Parameters for the 2016 Mw 6.2 Italy Amatrice Earthquake Revealed by InSAR Observations[J]. Geomatics and Information Science of Wuhan University, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

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

中央级公益性科研院所基本科研业务专项 ZDJ2017-29

中央级公益性科研院所基本科研业务专项 ZDJ2015-15

民用航天项目 D010102

中国地震局项目 Y201711

详细信息

Source Parameters for the 2016 Mw 6.2 Italy Amatrice Earthquake Revealed by InSAR Observations

Funds: 

The Special Fund for Basic Scientific Research of Institute of Crustal Dynamics ZDJ2017-29

The Special Fund for Basic Scientific Research of Institute of Crustal Dynamics ZDJ2015-15

Civil Aerospace Project D010102

the China Earthquake Administration Project Y201711

More Information
  • 摘要: 2016年8月24日,意大利中部阿马特里切(Amatrice)地区发生Mw 6.2地震。采用ALOS-2条带模式和SENTINEL-1A宽幅模式的合成孔径雷达(synthetic aperture radar,SAR)数据分别进行SAR差分干涉测量处理,获取了该地震的同震形变场。结果显示,本次地震造成意大利中部地区发生明显的地壳形变,在雷达视线向最大沉降量达19.6 cm。基于合成孔径雷达干涉测量(interferometry synthetic aperture radar,InSAR)和GPS同震形变场数据对此次地震的发震断层进行联合反演,通过改进倾角和平滑系数获取方法,得到了最优滑动分布模型。通过使用单断层模型和双断层模型进行反演可知,双断层模型反演结果优于单断层反演结果,两种模型下反演模型相关系数分别为0.85和0.89,发震断层走向分别为160°和158°,倾角分别为44°和46°,倾滑分布主要位于地下5~7 km,平均倾滑角为-80°,最大倾滑量0.9 m位于地壳深度5 km处,该发震断层是亚平宁冲断带的一部分,为NW-SE向延伸的正断层,断层长约20 km。综合使用地震同震形变场和GPS数据对震源机制进行反演、模拟和分析,获取了高精度的震源参数,可以为分析地震危险性和断层破裂参数等提供数据支持。
  • 图  1  2016年8月24日阿马特里切地震的构造背景场

    Figure  1.  Tectonic Background of the Italy Earthquake on August 24, 2016

    图  2  阿马特里切Mw 6.2地震InSAR同震形变场

    Figure  2.  InSAR Coseismic Deformation Field of the Amatrice Mw 6.2 Earthquake

    图  3  GPS站点得到的水平形变场及模拟得到的水平形变场

    Figure  3.  Horizontal Deformation Fields Obtained by GPS Site and Simulation

    图  4  基于InSAR反演的单断层和双断层空间滑动模型

    Figure  4.  Single and Double Fault Slip Models Determined with InSAR Datasets

    图  5  双断层模型InSAR资料约束下同震形变场、最优拟合模型预测位移场以及残差分布

    Figure  5.  Coseismic Deformation Field, Predicted Displacement Field and Residual Distribution for InSAR in Double Fault Model

    表  1  2016年阿马特里切地震干涉图详细信息

    Table  1.   Interferograms Used for the 2016 Italy Earthquake

    卫星类型 轨道号 飞行方向 成像模式 波长/cm 获取时间 时间基线/d 入射角/(°)
    主影像 从影像
    ALOS-2 197 升轨 条带 24.2 20160127 20160824 210 35
    SENTINEL-1A 44 升轨 宽幅 5.6 20160815 20160827 12 43
    下载: 导出CSV

    表  2  阿马特里切地震断层震源参数

    Table  2.   Fault Geometric Parameters of the Italy Earthquake

    来源 位置/(°) 深度/km 震源机制/(°) 断层大小/km 震级
    经度 纬度 走向 倾角 滑动角 长度 宽度
    USGS[1] 13.19 42.72 4.4 165 49 -78 6.2
    GCMT[28] 13.22 42.64 12.0 145 38 -101 6.2
    INGV[29] 13.22 42.71 5.0 155 41 -93 6.0
    Tinti等[30] 13.23 42.70 7.3 156 50 -120~-70 26 12 6.1
    Lavecchia等[11] 13.25/13.22 42.70/42.78 0.5/0.6 161/161 52/42 -85/-85 15/13 10/10 6.2
    Xu等[27] 13.26 42.73 2.7 164 43.6 -89.7 36 18 6.2
    本文 单断层 13.22 42.73 5 165 45 -80 25 10 6.2
    双断层 断层1 13.20 42.68 5 160 44 -80 10 10 6.2
    断层2 13.22 42.73 5 158 46 -85 10 10
    注:Lavecchia等[11]使用的双断层模型,因此有两组参数
    下载: 导出CSV
  • [1] USGS. Poster of the Central Italy Earthquake of 24 August 2016-Magnitude 6.2[OL]. http://earthquake.Usgs.gov/earthquakes/eqarchives/poster/2016/20160824.php, 2016
    [2] 谭凯, 王琪, 申重阳.用大地测量数据反演2001年昆仑山地震[J].大地测量与地球动力学, 2004, 24(3):47-50 http://d.old.wanfangdata.com.cn/Periodical/dkxbydz200403010

    Tan Kai, Wang Qi, Shen Chongyang. Using Geodetic Data to Inverse Co-seismic Dislocation of 2001 Kunlun Earthquake[J]. Journal of Geodesy and Geodynamics, 2004, 24(3):47-50 http://d.old.wanfangdata.com.cn/Periodical/dkxbydz200403010
    [3] 屠泓为, 汪荣江, 刁法启, 等.运用SDM方法研究2001年昆仑山口西Ms 8.1地震破裂分布:GPS和InSAR联合反演的结果[J].地球物理学报, 2016, 56(6):2103-2112 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201606016.htm

    Tu Hongwei, Wang Rongjiang, Diao Faqi, et al. Slip Model of the 2001 Kunlun Mountain Ms 8.1 Earthquake by SDM:Joint Inversion from GPS and InSAR Data[J]. Chinese Journal of Geophysics, 2016, 56(6):2103-2112 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201606016.htm
    [4] Li Zhenhong, Elliott J R, Feng Wanpeng, et al. The 2010 Mw 6.8 Yushu (Qinghai, China) Earthquake:Constraints Provided by InSAR and Body Wave Seismology[J]. Journal of Geophysical Research:Solid Earth, 2011, 116:B10302 doi:  10.1029/2011JB008358/full
    [5] 李永生, 申文豪, 温扬茂, 等. 2015年尼泊尔Mw 7.8地震震源机制InSAR反演及强地面运动模拟[J].地球物理学报, 2016, 59(4):1359-1370 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201604017.htm

    Li Yongsheng, Shen Wenhao, Wen Yangmao, et al. Source Parameters for the 2015 Nepal Earthquake Revealed by InSAR Observations and Strong Ground Motion Simulation[J]. Chinese Journal of Geophysics, 2016, 59(4):1359-1370 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201604017.htm
    [6] D'Agostino N, Mantenuto S, D'Anastasio E, et al. Contemporary Crustal Extension in the Umbria-Marche Apennines from Regional CGPS Networks and Comparison Between Geodetic and Seismic Deformation[J]. Tectonophysics, 2009, 476(1-2):3-12 http://www.sciencedirect.com/science/article/pii/S0040195108004575
    [7] Huang Monghan, Fielding E J, Liang Cunren, et al. Coseismic Deformation and Triggered Landslides of the 2016 Mw 6.2 Amatrice Earthquake in Italy[J]. Geophysical Research Letters, 2017, 44(3):1266-1274 doi:  10.1002/2016GL071687/abstract
    [8] Pucci S, de Martini P M, Civico R, et al. Coseismic Ruptures of the 24 August 2016, Mw 6.0 Amatrice Earthquake (Central Italy)[J]. Geophysical Research Letters, 2017, 44(5):2138-2147 doi:  10.1002/2016GL071859/pdf
    [9] Ekstrom G, Andrea M, Enzo B, et al. Moment Tensor Analysis of the Central Italy Earthquake Sequence of September-October 1997[J]. Geophysical Research Letters, 1998, 25(11):1971-1974 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=16ce55982bab3fbc7d93303a83001a48
    [10] Amato A, Azzara R, Chiarabba C, et al. The 1997 Umbria-Marche, Italy, Earthquake Sequence:A First Look at the Main Shocks and Aftershocks[J]. Geophysical Research Letters, 1998, 25(15):2861-2864 doi:  10.1029/98GL51842/full
    [11] Lavecchia G, Castaldo R, de Nardis R, et al. Ground Deformation and Source Geometry of the August 24, 2016 Amatrice Earthquake (Central Italy) Investigated Through Analytical and Numerical Modeling of DInSAR Measurements and Structural-Geological Data[J]. Geophysical Research Letters, 2016, 43:12389-12398 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=42d31df786e281366597e88879266387
    [12] Michetti A M. Surface Faulting During the August 24, 2016, Central Italy Earthquake[J]. Annals of Geophysics, 2016, 59:1-8
    [13] Yano T E, Shao Guangfu, Liu Qiming, et al. Finite Fault Kinematic Rupture Model of the 2009 Mw 6.3 L'Aquila Earthquake from Inversion of Strong Motion, GPS and InSAR Data[C]. American Geophysical Union Fall Meeting, San Francisco, California, USA, 2009
    [14] Temblor: Italy Earthquake Leaves Seismic Gaps that were Last Filled by Three Large Earthquakes in 1703[OL].http://temblor.net/earthquake-insights/gaps-persist-beyond-ends-m6-2-rieti-italy-earthquake-progressive-sequence-large-shocks-struck-1703-1216/, 2016
    [15] Cinti F R, Cucci L, Marra F, et al. The 1997 Umbria-Marche (Italy) Earthquake Sequence:Relationship Between Ground Deformation and Seismogenic Structure[J]. Geophysical Research Letters, 1999, 26(7):895-898 doi:  10.1029/1999GL900142/full
    [16] Wei Shengji, Barbot S, Graves R, et al. The 2014 Mw 6.1 South Napa Earthquake:A Unilateral Rupture with Shallow Asperity and Rapid Afterslip[J]. Seismological Research Letters, 2015, 86(2A):344-354
    [17] 李永生, 冯万鹏, 张景发, 等. 2014年美国加州纳帕Mw 6.1地震断层参数的SENTINEL-1A InSAR反演[J].地球物理学报, 2015, 58(7):2339-2349 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201507012.htm

    Li Yongsheng, Feng Wanpeng, Zhang Jingfa, et al. Coseismic Slip of the 2014 Mw 6.1 Napa, California Earthquake Revealed by SENTINEL-1A InSAR[J]. Chinese Journal of Geophysics, 2015, 58(7):2339-2349 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201507012.htm
    [18] RING-Rete Integrata Nazionale GPS[OL].http://ring.gm.ingv.it/?page_id=702, 2017
    [19] Cheloni D, Serpelloni E, Devoti R, et al. GPS Observations of Coseismic Deformation Following the 2016, August 24, Mw 6 Amatrice Earthquake (Central Italy):Data, Analysis and Preliminary Fault Model[J]. Annals of Geophysics, 2016, 59(5):1-8 https://www.annalsofgeophysics.eu/index.php/annals/article/view/7269
    [20] Cheloni D. Geodetic Model of the 2016 Central Italy Earthquake Sequence Inferred from InSAR and GPS Measurements[C]. European Geosciences Union General Assembly 2017, Vienna, Austria, 2017
    [21] Lohman R B, Simons M. Some Thoughts on the Use of InSAR Data to Constrain Models of Surface Deformation:Noise Structure and Data Downsampling[J]. Geochemistry Geophysics Geosystems, 2005, 6(1):359-361 doi:  10.1029-2004GC000841/
    [22] Feng Wanpeng, Li Zhenhong, Elliott J R, et al. The 2011 Mw 6.8 Burma Earthquake:Fault Constraints Provided by Multiple SAR Techniques[J]. Geophysical Journal International, 2013, 195(1):650-660 http://gji.oxfordjournals.org/content/195/1/650
    [23] 冯万鹏, 李振洪. InSAR资料约束下震源参数的PSO混合算法反演策略[J].地球物理学进展, 2010, 25(4):1189-1196 doi:  10.3969/j.issn.1004-2903.2010.04.007

    Feng Wanpeng, Li Zhenhong. A Novel Hybrid PSO/Simplex Algorithm for Determining Earthquake Source Parameters Using InSAR Data[J]. Progress in Geophysics, 2010, 25(4):1189-1196 doi:  10.3969/j.issn.1004-2903.2010.04.007
    [24] Li Yongsheng, Jiang Wenliang, Zhang Jingfa, et al. Space Geodetic Observations and Modeling of 2016 Mw 5.9 Menyuan Earthquake:Implications on Seismogenic Tectonic Motion[J]. Remote Sensing, 2016, 8(6):519-530 http://adsabs.harvard.edu/abs/2016RemS....8..519L
    [25] Okada Y. Surface Deformation due to Shear and Tensile Faults in a Half-Space[J]. Bulletin of the Seismological Society of America, 1985, 92(2):1018-1040 http://gji.oxfordjournals.org/cgi/ijlink?linkType=ABST&journalCode=ssabull&resid=75/4/1135
    [26] Harris R A, Segall P. Detection of a Locked Zone at Depth on the Parkfield, California, Segment of the San Andreas Fault[J]. Journal of Geophysical Research:Solid Earth, 1987, 92(B8):7945-7962 doi:  10.1029/JB092iB08p07945/full
    [27] Xu Guangyu, Xu Caijun, Wen Yangmao, et al. Source Parameters of the 2016-2017 Central Italy Earthquake Sequence from the Sentinel-1, ALOS-2 and GPS Data[J]. Remote Sensing, 2017, 9(11):1182-1203 http://adsabs.harvard.edu/abs/2017RemS....9.1182X
    [28] Global Centroid Moment Tensor[OL].http://www.globalcmt.org/CMTsearch.html, 2017
    [29] Institute of Geophysics and Volcanology[OL].http://info.terremoti.ingv.it/, 2017
    [30] Tinti E, Scognamiglio L, Michelini A, et al. Slip Heterogeneity and Directivity of the ML 6.0, 2016, Amatrice Earthquake Estimated with Rapid Finite-Fault Inversion[J]. Geophysical Research Letters, 2016, 43:10745-10752 doi:  10.1002/2016GL071263/pdf
  • [1] 李鹏, 高梦瑶, 李振洪, 王厚杰.  阿尔金断裂带宽幅InSAR对流层延迟估计方法评估 . 武汉大学学报 ● 信息科学版, 2020, 45(6): 879-887. doi: 10.13203/j.whugis20190236
    [2] 刘洋, 许才军, 温扬茂, 李志才.  杂多Mw 5.9级地震断层滑动的InSAR反演及边界元分析 . 武汉大学学报 ● 信息科学版, 2020, 45(11): 1678-1686. doi: 10.13203/j.whugis20190368
    [3] 刘洋, 许才军, 温扬茂.  门源Mw5.9级地震形变InSAR观测及区域断裂带深部几何形态 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 1035-1042. doi: 10.13203/j.whugis20190069
    [4] 温扬茂, 冯怡婷.  地震破裂模型约束的中国阿里地震三维形变场 . 武汉大学学报 ● 信息科学版, 2018, 43(9): 1369-1375. doi: 10.13203/j.whugis20160450
    [5] 何永红, 朱建军, 靳鹏伟.  一种使用剪切波变换的干涉图滤波算法 . 武汉大学学报 ● 信息科学版, 2018, 43(7): 1008-1014. doi: 10.13203/j.whugis20160236
    [6] 李水平, 王琪, 陈刚, 乔学军, 杨少敏, 何平, 陈超.  尼泊尔Mw 7.9级地震同震垂直位移与断层运动模型 . 武汉大学学报 ● 信息科学版, 2017, 42(10): 1489-1496. doi: 10.13203/j.whugis20160057
    [7] 许才军, 何平, 温扬茂, 张磊.  日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束 . 武汉大学学报 ● 信息科学版, 2012, 37(12): 1387-1391.
    [8] 温扬茂, 许才军, 刘洋, 何平.  利用断层自动剖分技术的2008年青海大柴旦Mw6.3级地震InSAR反演研究 . 武汉大学学报 ● 信息科学版, 2012, 37(4): 458-462.
    [9] 许才军, 江国焰, 王浩, 温扬茂.  基于GIS的InSAR结果分析方法及在汶川Mw7.9级地震同震解释中的应用 . 武汉大学学报 ● 信息科学版, 2011, 36(4): 379-383.
    [10] 吴宏安, 汤益先, 张红, 王超.  基于ALOS/PALSAR轨道参数的干涉平地效应消除研究 . 武汉大学学报 ● 信息科学版, 2010, 35(1): 92-96.
    [11] 许才军, 王乐洋.  大地测量和地震数据联合反演地震震源破裂过程研究进展 . 武汉大学学报 ● 信息科学版, 2010, 35(4): 457-462.
    [12] 温扬茂, 许才军.  基于敏感度的迭代拟合法反演玛尼Ms7.9级地震滑动分布 . 武汉大学学报 ● 信息科学版, 2009, 34(6): 732-735.
    [13] 鄢子平, 李振洪.  InSAR大气水汽改正模型的比较应用研究 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 723-726.
    [14] 许小兰, 李斐.  InSAR及其差分技术在仲巴地震中的应用 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 727-729.
    [15] 许才军, 温扬茂.  基于InSAR数据的西藏玛尼Ms7.9级地震的地壳不均匀性研究 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 846-849.
    [16] 王健, 向茂生, 李绍恩.  一种基于InSAR相干系数的SAR阴影提取方法 . 武汉大学学报 ● 信息科学版, 2005, 30(12): 1063-1066.
    [17] 李振洪, 刘经南, 许才军.  InSAR数据处理中的误差分析 . 武汉大学学报 ● 信息科学版, 2004, 29(1): 72-76.
    [18] 许才军, 王华.  InSAR相位解缠算法比较及误差分析 . 武汉大学学报 ● 信息科学版, 2004, 29(1): 67-71.
    [19] 许才军, 王华, 黄劲松.  GPS与INSAR数据融合研究展望 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 58-61,78.
    [20] 王乐洋, 孙龙翔, 许光煜.  利用GPS观测数据反演震源参数的单纯形组合加权距离灰狼新算法 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210114
  • 加载中
图(5) / 表(2)
计量
  • 文章访问数:  733
  • HTML全文浏览量:  56
  • PDF下载量:  228
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-12-06
  • 刊出日期:  2019-01-05

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

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

    中央级公益性科研院所基本科研业务专项 ZDJ2017-29

    中央级公益性科研院所基本科研业务专项 ZDJ2015-15

    民用航天项目 D010102

    中国地震局项目 Y201711

    作者简介:

    张庆云, 博士生, 主要从事InSAR数据处理研究。zhangqingyun1107@163.com

    通讯作者: 李永生, 副研究员。yongsheng2009@gmail.com
  • 中图分类号: P237

摘要: 2016年8月24日,意大利中部阿马特里切(Amatrice)地区发生Mw 6.2地震。采用ALOS-2条带模式和SENTINEL-1A宽幅模式的合成孔径雷达(synthetic aperture radar,SAR)数据分别进行SAR差分干涉测量处理,获取了该地震的同震形变场。结果显示,本次地震造成意大利中部地区发生明显的地壳形变,在雷达视线向最大沉降量达19.6 cm。基于合成孔径雷达干涉测量(interferometry synthetic aperture radar,InSAR)和GPS同震形变场数据对此次地震的发震断层进行联合反演,通过改进倾角和平滑系数获取方法,得到了最优滑动分布模型。通过使用单断层模型和双断层模型进行反演可知,双断层模型反演结果优于单断层反演结果,两种模型下反演模型相关系数分别为0.85和0.89,发震断层走向分别为160°和158°,倾角分别为44°和46°,倾滑分布主要位于地下5~7 km,平均倾滑角为-80°,最大倾滑量0.9 m位于地壳深度5 km处,该发震断层是亚平宁冲断带的一部分,为NW-SE向延伸的正断层,断层长约20 km。综合使用地震同震形变场和GPS数据对震源机制进行反演、模拟和分析,获取了高精度的震源参数,可以为分析地震危险性和断层破裂参数等提供数据支持。

English Abstract

张庆云, 李永生, 罗毅, 张景发, 杨建钦. 2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演[J]. 武汉大学学报 ● 信息科学版, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
引用本文: 张庆云, 李永生, 罗毅, 张景发, 杨建钦. 2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演[J]. 武汉大学学报 ● 信息科学版, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
ZHANG Qingyun, LI Yongsheng, LUO Yi, ZHANG Jingfa, YANG Jianqin. Source Parameters for the 2016 Mw 6.2 Italy Amatrice Earthquake Revealed by InSAR Observations[J]. Geomatics and Information Science of Wuhan University, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
Citation: ZHANG Qingyun, LI Yongsheng, LUO Yi, ZHANG Jingfa, YANG Jianqin. Source Parameters for the 2016 Mw 6.2 Italy Amatrice Earthquake Revealed by InSAR Observations[J]. Geomatics and Information Science of Wuhan University, 2019, 44(1): 118-124. doi: 10.13203/j.whugis20170252
  • 2016年8月24日3时,意大利阿马特里切(Amatrice)发生Mw 6.2级地震。美国地质调查局给出的震源机制解显示此次地震震中位置为42.72°N、13.18°E,震源深度约4.4 km[1]。意大利地球物理与火山学研究所(National Institute of Geophysics and Vulcanology,INGV)的分析结果显示此次地震是由西倾正断层破裂所致。由于震区建筑相对脆弱,所产生的破坏较大,共造成300人死亡,数百人受伤,大量建筑物损坏。地震震源机制反演有助于确定发震断层的几何形态,理解发震断层的运动学特征,从而有利于提高该地区地震危险性认知,规避地震风险。

    大量研究表明,以合成孔径雷达干涉测量(interferometry synthetic aperture radar,InSAR)同震形变场为约束,可以确定地震的滑动分布[2-5]。目前常用位错模型进行分析,通过调整反演和断层参数,确定最优的同震滑动分布。模拟退火法[2]和最速下降方法[3](steepest descent method,SDM)发展较早,但这两种方法算法相对复杂,反演中对算法自身控制参数选择的要求较高。本文选用随机搜索粒子群优化算法(particle swarm optimization,PSO),该反演方法基于InSAR资料进行震源参数反演,对断层破裂的几何参数和倾滑分布进行反演,被广泛应用到同震形变反演中,如中国青海省玉树地震、尼泊尔地震等[4-5]

    • 阿马特里切地震的震中位于意大利亚平宁山脉,位于欧亚板块和非洲板块的交接地带,地质构造非常复杂。地中海地区分布着若干个子块体,彼此间交互作用,使得亚平宁山脉所在的亚平宁半岛呈逆时针旋转,形成了意大利地区主要的地震活动带——亚平宁冲断带。亚平宁冲断带的形成主要与中新世-上新世时期亚得里亚板块向欧亚板块的持续俯冲有关[6],其中第勒尼安盆地的拉张运动速率大于非洲板块和欧亚板块之间北北西向的挤压运动速率,导致亚平宁冲断带发生东西向伸展,从而引发本次地震[7-8]

      其构造背景如图 1所示,黑白相间的沙滩球表征地震震源机制解;主要活动断层有古比奥断层F1、瓜多尔塔迪诺断层F2、科尔弗里奥托断层F3、博韦-维托尔断层F4、诺尔恰断层F5、高扎诺断层F6、拉奎拉断层F7。

      图  1  2016年8月24日阿马特里切地震的构造背景场

      Figure 1.  Tectonic Background of the Italy Earthquake on August 24, 2016

      意大利发震断层大多为正断层,属于意大利中部的正断层系[9]。本次阿马特里切地震是该地区自2009年拉奎拉Mw 6.3地震以来破坏性最强的地震。地震发生在高扎诺断层(F6)北段、博韦-维托尔断层(F4)的南段,阿马特里切地震的地表破裂带位于2009年拉奎拉地震以及1997年翁布里亚-马尔凯地震序列范围内(见图 1)。意大利中部的历史地震以及发震断层分段分布特征显示该地区的地震常以序列或者震群的形式呈现,这说明震级在5~6级之间的余震的危险性不容忽视,有可能会触发周边断层的较大震级的破裂性地震[10-15]。阿马特里切Mw 6.2级地震发生后1 h以内又发生了一次Mw 5.4级余震,震后两天内共记录到2级以上余震600余次。余震空间分布显示了发震断层是一个NW-SE走向的活断层。

    • 欧洲太空局的SENTINEL-1A卫星和日本宇宙航空研究开发机构的ALOS-2卫星都获取了大量阿马特里切震区SAR影像。本文利用开源软件JPL/Caltech ROI_PAC和商业软件Gamma对ALOS-2数据和SENTINEL-1A数据进行差分干涉处理,得到不同卫星平台获取的同震形变场(见图 2)。

      图  2  阿马特里切Mw 6.2地震InSAR同震形变场

      Figure 2.  InSAR Coseismic Deformation Field of the Amatrice Mw 6.2 Earthquake

      阿马特里切地震的余震主要集中在震后1~2 h内,而ALOS-2震后影像获取时间为8月24日23时,故较准确地捕获了同震形变场。选用20160127-20160824干涉对,时间基线为210 d,详细参数见表 1。由于ALOS-2卫星的SAR传感器为L波段载荷,因此此干涉图具有很高的相干性,且条带模式可以获取完整的地震形变区域(见图 2(a))。

      表 1  2016年阿马特里切地震干涉图详细信息

      Table 1.  Interferograms Used for the 2016 Italy Earthquake

      卫星类型 轨道号 飞行方向 成像模式 波长/cm 获取时间 时间基线/d 入射角/(°)
      主影像 从影像
      ALOS-2 197 升轨 条带 24.2 20160127 20160824 210 35
      SENTINEL-1A 44 升轨 宽幅 5.6 20160815 20160827 12 43

      SENTINEL-1A卫星获取的宽幅数据在地质灾害领域应用广泛[16-17]。本文使用SENTINEL-1A卫星20160815-20160827干涉对,时间基线为12 d,干涉对参数见表 1。数据可以完整覆盖地震形变场,且相干性高,形变场特征明显(见图 2(b))。

      图 2中,黑线表示断层,蓝色五角星为震中,红色负值区表示沿卫星飞行方向发生沉降,蓝色正值区表示沿卫星飞行方向发生抬升。两种卫星平台获得的形变场结果有较高的一致性,在整个极震区有两个形变中心,在主形变区域沿卫星视线向(line of sight,LOS)表现为沉降,阿库莫利地区的沉降量最大。ALOS-2升轨数据的LOS向最大形变量为19.6 cm,SENTINEL-1A升轨数据的LOS向最大形变量为18.9 cm。

    • 2016年8月24日意大利中部地震发生后,意大利地球物理与火山学研究所公布了震前6天至震后4天的GPS数据[18],本文选取研究区内的27个GPS站点,处理得到水平同震形变场(图 3中红色箭头),并基于InSAR和GPS结果联合约束反演得到GPS站点的同震位移(图 3中紫色箭头),图 3中蓝色五角星表示主震,彩色圆点表示余震位置。对比模拟得到的水平同震形变场与实测值,两者高度一致。测得本次地震位移最大的GPS台站为NRCI和LNSS,位于震中的西南侧,GPS测得的形变场形态与NW-SE走向正断层的拉张运动特征一致[19-20]

      图  3  GPS站点得到的水平形变场及模拟得到的水平形变场

      Figure 3.  Horizontal Deformation Fields Obtained by GPS Site and Simulation

    • 以ALOS-2和SENTINEL-1A卫星数据得到的同震形变场为反演约束,分析阿马特里切地震的断层参数及滑动机制,主要分为两步:①采用数据分辨率约束的四叉树方法对InSAR形变场降采样[21],降低计算量,提高反演速度。ALOS-2和SENTINEL-1A数据降采样后,数据量分别为813和688。②采用两步反演法对断层几何参数及倾滑分布进行反演[22],首先基于均一断层模型假设反演得到断层空间几何参数,然后通过分布式断层模型计算断层面的滑动量。整个反演过程采用基于粒子群优化算法的PSOKINV软件来进行[4-5, 23-24]

      使用均一断层模型进行反演主要是确定断层位置(经纬度)、顶部埋深、走向、倾角等参数。基于Okada弹性位错模型[25],使用PSO非线性优化算法比较模拟结果,确定反演最优参数,使适配函数达到最小解。适配函数定义为:

      $$ \sigma = \sqrt {\frac{{{{({\mathit{\boldsymbol{W}}}({\mathit{\boldsymbol{D}}} - {\mathit{\boldsymbol{GS}}}))}^2}}}{N}} $$ (1)

      式中,σ为适配函数;G为系数矩阵,表示均匀断层上倾滑量达到1 m时引起的地表运动响应;S表示倾滑矢量矩阵;W表示每个数据集的相对权重矩阵;D表示地表形变观测值矩阵;N表示形变观测值的个数[22]。以InSAR同震形变场为约束,用PSO方法在全参数域内找到使适配函数最小的解。由均一断层模型得到断层最优拟合参数:断层走向157°,倾角42°,震级Mw 6.2。

      使用分布式断层模型反演地震滑动空间分布。在考虑模型空间分辨率的基础上,使用较小单元尺寸对浅层子断层进行分割,对深部断层的分割尺寸随深度线性增长。分别对断层面在走向和倾向上进行扩展,选取大小合适的离散子断层进行反演。使用二次差分拉普拉斯算子约束滑动粗糙度[26],反演算法为:

      $$ \left[ \begin{array}{l} {\mathit{\boldsymbol{G}}}\\ {\alpha ^2}{\mathit{\boldsymbol{L}}} \end{array} \right]{\mathit{\boldsymbol{S}}} = \left[ \begin{array}{l} {\mathit{\boldsymbol{D}}}\\ {\mathit{\boldsymbol{0}}} \end{array} \right] $$ (2)

      式中,L表示二阶微分算子,用于评价滑动粗糙度;α2为平滑系数,用于平滑解析解[22]。均一模型可以确定最优拟合解,但是使用滑动分布模型时,最优倾角会有微小变动。

      依据模型的粗糙度、残差、倾角与平滑系数之间的关系,通过固定不同倾角值,改变不同的平滑系数,分析模型的粗糙度和残差变化趋势。使用均方根误差表示残差ξ,模型的粗糙度ψ定义为:

      $$ \psi = \sum\limits_{i = 1}^n {|{p_i}|} /(2n) $$ (3)

      式中,p= LSn表示子断层个数。使用log函数模型f(δ, α)=log(ψ+ξ)获得反演的最优粗糙度ψ和残差ξ,其中δ表示断层倾角。对模型归一化处理,获得残差和粗糙度折中曲线,从而确定最优平滑系数,并通过log函数获得最优倾角。

    • InSAR获取的形变场显示LOS向沉降区域有两个沉降中心,且本次地震处于高扎诺断层和博韦-维托尔断层之间。因此,考虑是否为双断层引起的地震。本文依据上述两步反演策略,联合InSAR和GPS数据,分别使用单断层和双断层进行震源机制反演。

      图 4表示横轴墨卡托投影坐标系(universal transverse Mercatol Projection, UTM)下断层位置。图 4(a)为单断层最优滑动分布模型。结果显示,本次地震有两个明显的滑动中心。断层的走向为167°,倾角为45°,断层倾滑以正断层为主,最大倾滑量为0.9 m,最优平滑系数为1.5。本结果与Xu等[27]的研究结果相符,详细参数见表 2

      图  4  基于InSAR反演的单断层和双断层空间滑动模型

      Figure 4.  Single and Double Fault Slip Models Determined with InSAR Datasets

      表 2  阿马特里切地震断层震源参数

      Table 2.  Fault Geometric Parameters of the Italy Earthquake

      来源 位置/(°) 深度/km 震源机制/(°) 断层大小/km 震级
      经度 纬度 走向 倾角 滑动角 长度 宽度
      USGS[1] 13.19 42.72 4.4 165 49 -78 6.2
      GCMT[28] 13.22 42.64 12.0 145 38 -101 6.2
      INGV[29] 13.22 42.71 5.0 155 41 -93 6.0
      Tinti等[30] 13.23 42.70 7.3 156 50 -120~-70 26 12 6.1
      Lavecchia等[11] 13.25/13.22 42.70/42.78 0.5/0.6 161/161 52/42 -85/-85 15/13 10/10 6.2
      Xu等[27] 13.26 42.73 2.7 164 43.6 -89.7 36 18 6.2
      本文 单断层 13.22 42.73 5 165 45 -80 25 10 6.2
      双断层 断层1 13.20 42.68 5 160 44 -80 10 10 6.2
      断层2 13.22 42.73 5 158 46 -85 10 10
      注:Lavecchia等[11]使用的双断层模型,因此有两组参数

      图 4(b)为双断层最优滑动分布模型。结果显示,两个滑动中心分布明显,发震的双断层走向分别为160°和158°,倾角分别为44°和46°,倾滑主要分布在地下5~7 km,断层倾滑以正断层为主,最大倾滑量为0.9 m,最优平滑系数为4.1。该结果与Lavecchia等[11]的研究成果相近,参数见表 2

      单断层和双断层模型下的反演模型相关系数分别为0.85和0.89,残差均方根误差分别为0.025 cm和0.021 cm。对比USGS[1]、GCMT[28]、INVG[29]等机构的结果(见表 2),本研究认为双断层模型更适合。

      图 5为双断层模型下最优拟合滑动模型给出的结果。图 5(a)图 5(d)分别表示ALOS-2和SENTINEL-1A卫星的InSAR观测值;图 5(b)图 5(e)分别表示基于ALOS-2和SENTINEL-1A形变场反演模拟得到的干涉形变场;图 5(c)图 5(f)分别表示图 5(b)图 5(e)两种情况下的残差。在双断层模型下,ALOS-2和SENTINEL-1A模拟结果的残余信号相对较小,可见双断层模型能更好地模拟出形变区域的双沉降中心。ALOS-2模拟结果中存在少量的残余信号,可能由部分残留的大气误差导致; SENTINEL-1A模拟的结果基本没有残余信号,吻合度较高。总体而言,使用两步反演策略基于双断层模型能很好地解释干涉形变场,证明本文确定的反演结果稳定、可靠。

      图  5  双断层模型InSAR资料约束下同震形变场、最优拟合模型预测位移场以及残差分布

      Figure 5.  Coseismic Deformation Field, Predicted Displacement Field and Residual Distribution for InSAR in Double Fault Model

    • 基于ALOS-2和SENTINEL-1A两种卫星传感器获取的地震同震形变场,采用两步反演策略进行地震破裂机制反演,并使用log函数获取最优倾角和平滑系数,获取完整的震源机制反演结果。对比单断层模型和双断层模型反演结果,双断层模型下残余误差的均方根误差小于单断层模型下的结果,且双断层模型反演相关系数为0.89,优于单断层的0.85。最终证明双断层模型获取的形变场和震源参数更加符合本次地震的破裂机制。本研究认为意大利阿马特里切Mw 6.2地震的发震断层为双断层,断层走向分别为160°和158°,倾角分别为44°和46°,断层倾滑角为-80°,最大倾滑量为0.9 m,深度5~7 km。从InSAR形变场、GPS形变场以及地震滑动分布反演结果可以看出,该地震造成的最大沉降量达19.6 cm,且发震断层为NW-SE向正断层,断层特征为拉张。本文的InSAR同震形变场反演结果与GCMT[28]、USGS[1]等机构的结果有很高的吻合度,说明该方法切实可行。本研究可以用于补充原始SAR形变场中的失相干区域,保证同震形变场的完整性,了解更详细的震源机制参数,为进一步研究该地区断层性质以及不同大地震之间的空间关系提供依据。

参考文献 (30)

目录

    /

    返回文章
    返回