中小地震三维形变场重构方法研究与同震滑动分布反演以2016年5月22日定日MW5.3地震为例

摘要/Abstract

摘要： 地震三维形变场对于研究地震发震机制等具有重要意义。已有的InSAR三维形变场重构研究中, 只有中强地震的实例。由于形变量级小、 InSAR方位向形变的误差较大, 中小地震的三维形变场重构易受噪声等影响。本文以2016年5月22日西藏定日M_W5.3地震为例, 开展中小地震三维形变场重构的尝试。首先基于InSAR技术获取了Sentinel-1升轨和降轨观测模式下的同震形变场, 再结合同震形变场的特点、区域构造特征等, 添加限定方程(走滑运动为0), 重构了同震三维形变场。结果显示, 震中附近以下降为主, 幅度达7 cm, 南北方向形变较小(约2 mm), 此区域还伴有2 cm的西向水平运动; 形变中心区域东西两侧部分区域均出现少许东向运动(1.5 cm)。由同震形变场特征判断此次地震以正断破裂为主。本文提出了基于连续性分层采样选取样本点方法, 以适应本地震形变场的实际情况。对所得的LOS向位移场和重构的三维形变场进行降采样, 反演得到了断层面上的滑动分布, 两种数据得到的结果相似, 最优发震断层的走向约181°, 倾角约45°, 断层错动平均滑动角约-87.1°, 平均滑动量约为3.6 cm, 最大滑动量位于深度6.5 km处, 相当于一次M_W5.4的地震。

关键词: 三维形变场重构, InSAR, 连续性分层采样, 反演, 中小地震

Abstract: Three-dimensional seismic deformation field is of great significance to the study of earthquake mechanism. It has been successfully reconstructed for only moderate or strong earthquakes in previous research. For small and medium earthquakes, the reconstruction of three-dimensional deformation field is vulnerable to noise and other effects, due to the small magnitude of displacement, and the large uncertainty in InSAR’s azimuthal deformation. Taking the May 22, 2016 Dingri M_W5.3 earthquake as an example, we attempt to reconstruct its three-dimensional deformation field. Firstly, the coseismic deformation field of the earthquake is obtained by using InSAR technology under the ascending and descending observation mode of Sentinel-1. According to the characteristics of coseismic deformation field, regional structure and so on, the three-dimensional deformation field of co-seismic is reconstructed by adding a limited equation (strike-slip motion is 0) in it. The results show that: the subsidence mainly distributes near the epicenter with the maximum amplitude reaching 7 cm, and the deformation in the north-south direction is small (about 2 mm), while accompanied by a horizontal westward motion of 2 cm; small eastward motions (1.5 cm) occur in the west and east parts of the deformation center region. From the deformation characteristics, the earthquake is illustrated as normal fault movement. In view of the characteristics of the seismic deformation field, a more suitable resampling point selection method based on continuous stratified sampling is proposed in this paper. Taking down-sample from the LOS displacement field and the reconstructed three-dimensional deformation field, the slip distribution on the fault plane is obtained, and the results from two data are similar. The best strike of the fault is about 181°, the optimal dip angle of the seismogenic fault about 45°, the average slip direction of the fault dislocation being about -87.1°, and the average slip is 3.6 cm. The maximum slip is located at a depth of 6.5 km. The observed fault slip is equivalent to an earthquake with a moment magnitude of M_W5.4.

Key words: Three-dimensional Deformation Field Reconstruction, InSAR, Continuous stratified sampling, Inversion, Small and Medium Earthquakes

中图分类号:

P315.7

纪润池, 申旭辉, 张景发, 田云锋. 中小地震三维形变场重构方法研究与同震滑动分布反演以2016年5月22日定日M_W5.3地震为例[J]. 地震, 2019, 39(4): 84-97.

JI Run-chi, SHEN Xu-hui, ZHANG Jing-fa, TIAN Yun-feng. Research on Reconstruction Method of Three-Dimensional Deformation Field for Small and Medium Earthquakes and Inversion of Co-seismic Slip DistributionA case study of May 22 2016 M_W5.3 Dingri earthquake[J]. EARTHQUAKE, 2019, 39(4): 84-97.

参考文献

[1] Fialko Y, Simons M, Agnew D. The complete (3-D) surface displacement field in the epicentral area of the 1999 M_W7.1 Hector Mine Earthquake, California, from space geodetic observations[J]. Geophysical Research Letters, 2001, 28(16): 3063-3066.
[2] Gudmundsson S, Sigmundsson F, Carstensen J M. Three-dimensional surface motion maps estimated from combined interferometric synthetic aperture radar and GPS data[J]. J Geophys Res, 2002, 107(B10): 2250.
[3] Wright T J, Elliott J R, Wang H, et al. Earthquake cycle deformation and the Moho: Implications for the rheology of continental lithosphere[J]. Tectonophysics, 2013, 609: 504-523.
[4] Bechor N B D, Zebker H A . Measuring two-dimensional movements using a single InSAR pair[J]. Geophys Res Lett, 2006, 33(16): L16311.
[5] 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 (London), 2005, 435(7040): 295-299.
[6] Wang H, Ge L, Xu C, et al. 3-D coseismic displacement field of the 2005 Kashmir earthquake inferred from satellite radar imagery[J]. Earth Planets & Space, 2007, 59(5): 343-349.
[7] Hu J, Wang Q J, Li Z W, et al. Retrieving three-dimensional coseismic displacements of the 2008 Gaize, Tibet earthquake from multi-path interferometric phase analysis[J]. Natural Hazards, 2014, 73(3): 1311-1322.
[8] Hu J, Li Z W, Ding X L, et al. 3D coseismic Displacement of 2010 Darfield, New Zealand earthquake estimated from multi-aperture InSAR and D-InSAR measurements[J]. Journal of Geodesy, 2012, 86(11): 1029-1041.
[9] Elliott J, Jolivet R, González P J, et al. Himalayan megathrust geometry and relation to topography revealed by the Gorkha earthquake[J]. Nature Geoscience, 2016.
[10] 罗海滨, 何秀凤, 刘焱雄. 利用DInSAR和GPS综合方法估计地表三维形变速率[J]. 测绘学报, 2008, 37(2): 168-171.
[11] 班保松, 伍吉仓, 陈永奇, 等. 联合GPS和InSAR观测结果计算汶川地震三维地表形变[J]. 大地测量与地球动力学, 2010, 30(4): 25-28.
[12] Guglielmino F, Bignami C, Bonforte A, et al. Analysis of satellite and in situ ground deformation data integrated by the SISTEM approach: The April 3, 2010 earthquake along the Pernicana fault (Mt. Etna—Italy) case study[J]. Earth & Planetary Science Letters, 2011, 312(3-4): 327-336.
[13] 王永哲. 基于InSAR的地表同震形变获取及震源参数反演研究[D]. 中南大学, 2012.
[14] 舒宁. 雷达影像干涉测量原理[M]. 武汉: 武汉大学出版社, 2003.
[15] Ito Y, Hosokawa M. Adegree-of-damage estimation model of earthquake damage using interferometric SAR data[J]. Electrical Engineering in Japan, 2003, 143(3): 49-57.
[16] 王超, 张红, 刘智. 星载合成孔径雷达干涉测量[M]. 北京: 科学出版社, 2002.
[17] 王晓文. 基于InSAR和MAI的电离层误差校正及同震三维形变场计算与断层滑动反演[D]. 西南交通大学, 2017.
[18] Hu J, Li Z W, Ding X L, et al. Resolving three-dimensional surface displacements from InSAR measurements: A review[J]. Earth-Science Reviews, 2014, 133(2): 1-17.
[19] Fialko Y, Simons M, Agnew D. The complete (3-D) surface displacement field in the epicentral area of the 1999 M_W7.1 Hector Mine Earthquake, California, from space geodetic observations[J]. Geophys Res Lett, 2001, 28(16): 3063-3066.
[20] Jonsson, S. Fault Slip Distribution of the 1999 M_W7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements[J]. Bull Seismol Soc Am, 2002, 92(4): 1377-1389.
[21] Wang C, Ding X, Shan X, et al. Slip distribution of the 2011 Tohoku earthquake derived from joint inversion of GPS, InSAR and seafloor GPS/acoustic measurements[J]. Journal of Asian Earth Sciences, 2012, 57(none).
[22] 汪致恒, 张瑞, 王晓文, 等. 2016年门源地震的InSAR同震形变监测与断层反演[J]. 遥感信息, 2018, 33(6).