利用Sentinel-1A数据反演2016年青海门源MW5.9地震的同震形变场及断层参数

地震 ›› 2017, Vol. 37 ›› Issue (3): 12-21.

利用Sentinel-1A数据反演2016年青海门源M_W5.9地震的同震形变场及断层参数

白泽朝¹, 徐青¹, 靳国旺¹, 张红敏¹, 刘辉¹, 董彦芳²

1.信息工程大学地理空间信息学院, 河南郑州 450001;
2.中国地震局地震预测研究所, 北京 100036

收稿日期:2016-10-09 发布日期:2019-08-09
作者简介:白泽朝(1991-), 男, 河北唐山人, 在读硕士研究生, 从事InSAR地表形变监测与研究。
基金资助:
国家自然科学基金(41071296,41474010,41304002)

Co-seismic Deformation Field and Fault Parameter Inversion of the 2016 Menyuan M_W5.9 Earthquake from Sentinel-1A Data

BAI Ze-chao,XU Qing,JIN Guo-wang,ZHANG Hong-min,LIU Hui,DONG Yan-fang

1.Geospatial Information Institute, Information Engineering University, Zhengzhou 450001, China 2.Institute of Earthquake Science, CEA, Beijing 100036, China

Received:2016-10-09 Published:2019-08-09

摘要/Abstract

摘要： 利用Sentinel-1A卫星IW模式的升、降轨雷达复数影像, 结合ESA提供的精密轨道数据进行干涉处理, 根据30 m格网间距ASTER GDEM V2数据去除地形相位获得升轨和降轨差分干涉图, 采用带有掩膜的最小费用流算法解缠差分干涉图来反演2016年1月21日发生在青海门源M_W5.9地震的降轨同震形变场, 并基于弹性半空间矩形位错模型反演该地震的断层滑动分布。结果表明: 地震造成的同震形变场长约30 km, 宽约14 km, 震中呈现抬升, 为形变极大值区, 升轨雷达视线向最大形变量约为5 cm, 降轨雷达视线向最大形变量约为7 cm; 该地震断层滑动为左旋走滑, 为走向133.6°、倾角46.7°的东北倾断层, 主要滑动量集中在1~7 km深度, 最大滑动量约2.4 m, 位于4 km的深度; 反演得出的地震矩为1.19×10¹⁸ Nm, 震级为M_W6.0。

关键词: 同震形变, Sentinel-1A, InSAR, 断层滑动分布, 2016年青海门源地震

Abstract: Based on ascending and descending tracks of the Sentinel-1A satellite IW model, combined with the precision orbit data provided by ESA, the differential interferogram is obtained by removing the terrain phase from the ASTEM GDEM V2 data with 30 m grid spacing. The differential interferogram is unwrapped by minimum cost flow algorithm in order to derive the precise coseismic surface displacements of the Menyuan M_W5.9 earthquake occurred on January 21, 2016. Then elastic half-space dislocation model are used to estimate the fault geometry. The results show that the deformation field is about 30 km long and 14 km wide. The maximum LOS uplifting value is about 5 cm in the ascending coseismic deformation field and about 7 cm in the descending coseismic deformation field.The earthquake fault is a normal fault with a strike of 133.6°, a dip of 46.7° and dipping to northeast.The slip mainly occurs within 1~7 km along the orientation, with a maximum slip of 2.4 m at the depth of 4 km. The determined geodetic moment is 1.19×10¹⁸ Nm.

Key words: Co-seismic deformation, Fault slip distribution Sentinel-1A, InSAR, The 2016 Menyuan M_W5.9 earthquake

中图分类号:

P315.7

白泽朝, 徐青, 靳国旺, 张红敏, 刘辉, 董彦芳. 利用Sentinel-1A数据反演2016年青海门源M_W5.9地震的同震形变场及断层参数[J]. 地震, 2017, 37(3): 12-21.

BAI Ze-chao,XU Qing,JIN Guo-wang,ZHANG Hong-min,LIU Hui,DONG Yan-fang. Co-seismic Deformation Field and Fault Parameter Inversion of the 2016 Menyuan M_W5.9 Earthquake from Sentinel-1A Data[J]. EARTHQUAKE, 2017, 37(3): 12-21.

参考文献

[1] 靳国旺, 徐青, 张红敏. 合成孔径雷达干涉测量[M]. 国防工业出版社, 2014.
[2] 靳国旺, 张红敏, 徐青. 雷达摄影测量[M]. 测绘出版社, 2015.
[3] 靳国旺, 张红敏, 徐青. 雷达遥感测绘技术基础[M]. 解放军出版社, 2014.
[4] GABRIEL K, GOLDSTEIN R M, ZEBKER H A. Mapping small elevation changes over large areas: Differential radar interferometry[J]. Journal of Geophysical Research, 1989, 94 : 9183-9191.
[5] Massonnet, Didier, et al. The displacement field of the Landers earthquake mapped by radar interferometry[J]. Nature, 1993, 364(6433): 138-142.
[6] Fialko Y, Sandwell D, Simons M, et al. Three-dimensional deformation caused by the Bam, Iran, earthquake and the origin of shallow slip deficit.[J]. Nature, 2005, 435(7040): 295-9.
[7] 孙建宝, 梁芳, 徐锡伟, 等. 升降轨道ASAR雷达干涉揭示的巴姆地震(M_W6.5)3D同震形变场[J]. 遥感学报, 2006, 10(4): 489-496.
[8] 洪顺英, 申旭辉, 单新建, 等. 基于升降轨ASAR的于田M_S7.3级地震同震形变场信息提取与分析[J]. 国土资源遥感, 2010, 04: 98-102.
[9] 温扬茂, 许才军, 刘洋, 等. 升降轨InSAR数据约束下的2007年阿里地震反演分析[J]. 测绘学报, 2015, 06: 649-654.
[10] 周辉, 冯光财, 李志伟, 等. 利用InSAR资料反演缅甸M_W6.8地震断层滑动分布[J]. 地球物理学报, 2013, 09: 3011-3021.
[11] 李永生, 冯万鹏, 张景发, 等. 2014年美国加州纳帕M_W6.1地震断层参数的Sentinel-1A InSAR反演[J]. 地球物理学报, 2015, 07: 2339-2349.
[12] 李永生, 申文豪, 温扬茂, 等. 2015年尼泊尔M_W7.8地震震源机制InSAR反演及强地面运动模拟[J]. 地球物理学报, 2016, 04: 1359-1370.
[13] 李媛, 刘峡, 刘希康, 等. 2016年门源M_S6.4地震前祁连山周边断层活动分析[J]. 大地测量与地球动力学, 2016, 04: 288-293.
[14] 马保起, 李德文. 祁连山中段门源盆地新构造运动的阶段划分[J]. 地质力学学报, 2008, 03: 201-211.
[15] 胡朝忠, 杨攀新, 李智敏, 等. 2016年1月21日青海门源6.4级地震的发震机制探讨[J]. 地球物理学报, 2016, 05: 1637-1646.
[16] Lanari R, Berardino P, Bonano M, et al. Sentinel-1 results: SBAS-DInSAR processing chain developments and land subsidence analysis[C]∥Geoscience and Remote Sensing Symposium. IEEE, 2015: 2836-2839.
[17] Scheiber R, Moreira A. Coregistration of interferometric SAR images using spectral diversity. IEEE Trans. Geosci. Remote Sens., 2000, 38: 2179–2191.
[18] 郭树松, 祝意青, 赵云峰. 青海门源M6.4地震前重力场变化特征及其动力机制分析[J]. 大地测量与地球动力学, 2016, 04: 294-297+307.
[19] 郭安宁, 李鑫, 白雪见, 等. 2016年1月21日青海门源6.4级地震及相关参数[J]. 地震工程学报, 2016, 01: 150-158.
[20] OKADA Y. Surface Deformation due to Shear and Tensile Faults in a Half-space[J]. Bulletin of the Seismological Society of America, 1985, 75 (4): 1135-1154.
[21] Marquardt D W. An algorithm for least-squares estimation of nonlinear parameters[J]. Journal of the society for Industrial and Applied Mathematics, 1963, 11(2): 431-441.
[22] Wright T J, Lu Z, Wicks C. Source Model for the M_W6.7, 23 October 2002, Nenana Mountain Earthquake (Alaska) from InSAR [J]. Geophysical Research Letters, 2003, 30 (18): 1974
[23] Funning G J, Parsons B, Wright T J. Fault slip in the 1997 Manyi, Tibet earthquake from linear elastic modelling of InSAR displacements[J]. Geophysical Journal International, 2007, 169(3): 988-1008.