2016年新疆阿克陶MS6.7地震同震形变特征和断层滑动分布反演

地震 ›› 2018, Vol. 38 ›› Issue (3): 81-91.

2016年新疆阿克陶M_S6.7地震同震形变特征和断层滑动分布反演

陈树¹, 董彦芳¹, 洪顺英¹, 刘泰²

1.中国地震局地震预测研究所, 北京 100036;
2.中国地震局地球物理研究所, 北京 100081

收稿日期:2018-03-08 发布日期:2019-08-15
通讯作者: 董彦芳,副研究员。E-mail:dyfgis@163.com
作者简介:陈树(1991-),男,湖北仙桃人,在读硕士研究生,主要从事InSAR地震应用方面研究
基金资助:
国家重点研发计划资助(2017YFB0502700), 国家自然科学基金(41304002), 中国地震局地震预测研究所基本科研业务费专项(2018IEF0303)共同资助

Co-seismic Deformation Field and Inversion of Fault Slip Distribution of the 2016 Aketao M_S6.7 Earthquake

CHEN Shu¹, DONG Yan-fang¹, HONG Shun-ying¹, LIU Tai²

1.Institute of Earthquake Forecasting, CEA, Beijing 100036, China;
2.Institute of Geophysics, CEA, Beijing 100081, China

Received:2018-03-08 Published:2019-08-15

摘要/Abstract

摘要： 本文针对2016年11月25日在新疆阿克陶发生的地震, 用差分干涉测量技术(D-InSAR)对3种不同观测模式的升、降轨数据进行处理, 提取了多视线向的同震形变场; 根据不同模式的LOS向形变量, 构建形变分解模型, 将其分解为垂直向形变量和沿断层走向形变量; 结合同震形变场特征与震源机制解, 采用单断层模型, 利用梯度下降法(SDM), 以Multi-LOS向形变进行约束, 反演了阿克陶地震的同震滑动分布特征。研究结果表明, 升、降轨LOS向同震形变场在发震断层两侧具有不同的形变特征, 发震断层走向近EW向; LOS向形变量分解表明, 此次地震破裂以右旋走滑为主; 滑动分布反演的形变残差介于0~5 cm之间, 发震断层的滑动量主要位于2~16 km深部, 最大滑动量可达1.02 m, 位于断层面深部5.83 km处, 最大滑动量处的滑动角为185.24°; 平均滑动角为181.32°, 平均滑动量为0.12 m; 滑动分布反演也证明该地震为右旋走滑破裂事件, 与LOS向形变分解结果一致; 当剪切模量μ=3.2 GPa时, 反演得到的地震矩震级约M_W6.6。

关键词: D-InSAR, 同震形变场, 滑动分布反演, 2016年阿克陶M_S6.7地震

Abstract: Utilizing differential interferometric synthetic aperture radar technique, we retrieved multi-LOS co-seismic deformation field of the 2016 Aketao Ms6.7 earthquake through processing 3 kinds of different observation pattern data. According to the strike-slip characteristic of the fault, the fault-perpendicular deformation was neglected, and the up-down and fault-parallel deformation was extracted based on the combination of Sentinel-1A and Sentinel-1B observation. Besides, with the combination of deformation field features and focal mechanism solution, we inverted effectively the seismic slip distribution using the steepest descent method. The results demonstrate that the ascending and descending LOS co-seismic deformation fields have distinct different deformation characteristics on both sides of the seismogenic fault. According to the co-seismic deformation fields, the seismogenic fault is in the near east-west direction. According to the deformation fields, the earthquake is dextral strike-slip. The results of inversion demonstrate that the deformation residuals of the inversion are effectively controlled, which are basically between 0±5 cm. Fault slips distribute at the depth of 2 to 16 km. Results show that the maximum slip is 1.02 m at the depth of 5.83 km, and the rake angle is 185.24°. The mean slip is 0.12 m and the mean rake angle is 181.32 °.The rake angles at the two sliding distribution centers are close to 180 °. These demonstrate that the earthquake is a dextral strike-slip event, which is consistent with the result of LOS decomposition. When the shear modulus is set to 3.2 GPa, the retrieved seismic moment magnitude is M_W6.60.

Key words: Co-seismic deformation field, D-InSAR, Inversion, The 2016 Aketao M_S6.7 earthquake

中图分类号:

P315.7

陈树, 董彦芳, 洪顺英, 刘泰. 2016年新疆阿克陶M_S6.7地震同震形变特征和断层滑动分布反演[J]. 地震, 2018, 38(3): 81-91.

CHEN Shu, DONG Yan-fang, HONG Shun-ying, LIU Tai. Co-seismic Deformation Field and Inversion of Fault Slip Distribution of the 2016 Aketao M_S6.7 Earthquake[J]. EARTHQUAKE, 2018, 38(3): 81-91.

参考文献

[1] Massonnet D, Rossi M, Carmona C, et al. The displacement field of the Landers earthquake mapped by radar interferometry[J]. Nature, 1993, 364(6433): 138-142.
[2] Ozawa S, Murakami M, Fujiwara S, et al. Synthetic aperture radar interferogram of the 1995 Kobe Earthquake and its geodetic inversion[J]. Geophys Res Lett, 1997, 24(18): 2327-2330.
[3] Tobita M, Fujiwara S, Ozawa S, et al. Deformation of the 1995 North Sakhalin earthquake detected by JERS-1/SAR interferometry[J]. Earth, Planets and Space, 1998, 50(4): 313-325.
[4] Massonnet D, Feigl K L. Satellite radar interferometric map of the coseismic deformation field of the M=6.1 Eureka Valley, California Earthquake of May 17, 1993[J]. Geophys Res Lett, 1995, 22(12): 1541-1544.
[5] Massonnet D, Feigl K L, Vadon H, et al. Coseismic deformation field of the M=6.7 Northridge, California Earthquake of January 17, 1994 recorded by two radar satellites using interferometry[J]. Geophys Res Lett, 1996, 23(9): 969-972.
[6] 王超, 刘智, 张红, 等. 张北—尚义地震同震形变场雷达差分干涉测量[J]. 科学通报, 2000, 45(23): 2550-2554
[7] 单新建, 柳稼航, 马超. 2001年昆仑山口西8.1级地震同震形变场特征的初步分析[J]. 地震学报, 2004, 26(5): 474-480.
[8] 洪顺英, 申旭辉, 单新建, 等. 基于D-InSAR技术的西藏改则地震同震形变场特征分析[J]. 地震, 2009, 29(4): 23-31.
[9] 洪顺英, 申旭辉, 单新建, 等. 基于升降轨ASAR的于田M_S7.3级地震同震形变场信息提取与分析[J]. 国土资源遥感, 2010, 2010(4): 98-102.
[10] Barrientos S E, Ward S N. The 1960 Chile earthquake: inversion for slip distribution from surface deformation[J]. Geophysical Journal of the Royal Astronomical Society, 2010, 103(3): 589-598.
[11] Jónsson S, Zebker H, Segall P, et al. Fault Slip Distribution of the 1999 M_W7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements[J]. Bull Seism Soc Am, 2002, 92(4): 1377-1389.
[12] 冯万鹏, 许力生, 许忠淮, 等. 利用InSAR资料反演2008年西藏改则M_W6.4和M_W5.9地震的断层参数[J]. 地球物理学报, 2009, 52(4): 983-993.
[13] 张国宏. Slip distribution and source parameters inverted from co-seismic deformation derived by InSAR technology of Wenchuan M_W7.9 earthquake基于InSAR同震形变场反演汶川M_W7.9地震断层滑动分布[J]. Chinese Journal of Geophysics, 2010, 53(2): 269-279.
[14] 洪顺英, 刘智荣, 申旭辉, 等. 多视线向InSAR形变场约束的改则地震滑动分布反演[J]. 遥感学报, 2015, 19(2): 288-301.
[15] 陈杰, 李涛, 孙建宝, 等. 2016年11月25日新疆阿克陶M_W6.6地震发震构造与地表破裂[J]. 地震地质, 2016, 38(4): 1160-1174.
[16] 张琳琳, 聂晓红, 高朝军, 等. 2016年11月25日阿克陶M_S6.7地震前后库仑应力变化分析[J]. 内陆地震, 2017, 31(2): 115-121.
[17] 李成绕, 薛东剑, 张露. 基于D-InSAR技术的新疆阿克陶地震形变监测[J]. 福建质量管理, 2017(17): 271.
[18] Wang S, Xu C, Wen Y, et al. Slip Model for the 25 November 2016 M_W6.6 Aketao Earthquake, Western China, Revealed by Sentinel-1 and ALOS-2 Observations[J]. Remote Sensing, 2017, 9(4): 325.
[19] Wang R, Diao F, Hoechner A. SDM-A geodetic inversion code incorporating with layered crust structure and curved fault geometry[C]. EGU General Assembly Conference, EGU General Assembly Conference Abstracts, 2013.
[20] 杨魁, 杨建兵, 江冰茹. Sentinel-1卫星综述[J]. 城市勘测, 2015(2): 24-27.
[21] Sandwell D, Mellors R, Tong X, et al. GMTSAR: An InSAR Processing System Based on Generic Mapping Tools[M]. Office of Scientific & Technical Information Technical Reports, 2011.
[22] Okada Y. Surface deformation due to shear and tensile faults in a half-space[J]. Bull Seism Soc Am, 1985, 92(2): 1018-1040.