Inversion of fault slip parameters based on DInSAR and automated fault model discretization method
CHEN Danlei1, LIU Guoxiang2, WANG Xiaowen2, WANG Lei1, PU Huilong1
1. The Third Engineering Institution of Surveying and Mapping of Sichuan, Chengdu 610500, China;
2. Faculty of Geosciences and Environmental Engineering, Sourthwest Jiaotong University, Chengdu 610031, China
The inversion of geometry and motion parameters of the seismogenic fault based on elastic dislocation model and constrained by coseismic deformation data obtained by DInSAR has great significance for understanding the activity characteristics of the fault and assessing the risk of potential disasters. The existing inversion methods of fault slip parameters are mostly based on the uniform discretization with rectangular dislocation units or artificially setting discretizaiton, which will lead to "pseudo-slip" or the phenomenon that the slip distribution is too smooth to reflect the slip details on the fault plane. Therefore, the automated fault model discetization method is introduced in this paper, which takes into account the principle of the minimized residuals and the optimal smoothing scales of the model, so the reliable solution can be obtained under the effective constraints. The inversion of geometry and slip parameters of fault in Bam earthquake is taken as an example. The experimental results show that using the automated fault model discretization method to invert the motion parameters of the single fault can generate reliable results.
陈丹蕾, 刘国祥, 王晓文, 王蕾, 蒲慧龙. 基于DInSAR技术与断层自动剖分方法反演断层滑动参数[J]. 国土资源遥感, 2016, 28(3): 25-30.
CHEN Danlei, LIU Guoxiang, WANG Xiaowen, WANG Lei, PU Huilong. Inversion of fault slip parameters based on DInSAR and automated fault model discretization method. REMOTE SENSING FOR LAND & RESOURCES, 2016, 28(3): 25-30.
[1] 洪顺英,申旭辉,单新建,等.基于升降轨ASAR的于田Ms 7.3级地震同震形变场信息提取与分析[J].国土资源遥感,2010,22(4):98-102.doi:10.6046/gtzyyg.2010.04.20. Hong S Y,Shen X H,Shan X J,et al.The calculation and analysis of the co-seismic deformation field of Yutian Ms 7.3 earthquake basing on the ascending and descending orbit ASAR data[J].Remote Sensing for Land and Resource,2010,22(4):98-102.doi:10.6046/gtzyyg.2010.04.20.
[2] 罗旭巍,孙建宝,沈正康,等.基于InSAR同震形变观测反演2010年新西兰南岛Mw7.1 Darfield地震同震破裂分布[J].地球物理学报,2013,56(8):2613-2624. Luo X W,Sun J B,Shen Z K,et al.Co-seismic slip distribution of 2010 Darfield, New Zealand Mw7.1 earthquake inverted using InSAR measurements[J].Chinese Journal of Geophysics,2013,56(8):2613-2624.
[3] 薛莲,孙建宝,沈正康.2010年1月12日海地Mw7.0级地震InSAR同震形变观测及同震滑动分布反演[J].地震地质,2011,33(1):157-174. Xue L,Sun J B,Shen Z K.InSAR coseismic deformation observation of the Jan 12th, 2010 HaiTi earthquake and its coseismic slip distribution inversion[J].Seismology and Geology,2011,33(1):157-174.
[4] 周辉,冯光财,李志伟,等.利用InSAR资料反演缅甸Mw6.8地震断层滑动分布[J].地球物理学报,2013,56(9):3011-3021. Zhou H,Feng G C,Li Z W,et al.The fault slip distribution of the Myanmar Mw6.8 earthquake inferred from InSAR measurements[J].Chinese Journal of Geophysics,2013,56(9):3011-3021.
[5] 温扬茂,何平,许才军,等.联合Envisat和ALOS卫星影像确定L'Aquila地震震源机制[J].地球物理学报,2012,55(1):53-65. Wen Y M,He P,Xu C J,et al.Source parameters of the 2009 L'Aquila earthquake, Italy from Envisat and ALOS satellite SAR images[J].Chinese Journal of Geophysics,2012,55(1):53-65.
[6] Resor P G,Pollard D D,Wright T J,et al.Integrating high-precision aftershock locations and geodetic observations to model coseismic deformation associated with the 1995 Kozani-Grevena earthquake,Greece[J].Journal of Geophysical Research:Solid earth,2005,110(B9):B09402.
[7] Walters R J,Elliott J R,D'Agostino N,et al.The 2009 L'Aquila earthquake(central Italy):A source mechanism and implications for seismic hazard[J].Geophysical Research Letters,2009,36(17):L17312.
[8] 张国宏,屈春燕,汪驰升,等.基于GPS和InSAR反演汶川Mw7.9地震断层滑动分布[J].大地测量与地球动力学,2010,30(4):19-24. Zhang G H,Qu C Y,Wang C S,et al.Inversion of slip distribution of 2008 Wenchuan Mw7.9 earthquake constrained jointly by InSAR and GPS measurements[J].Journal of Geodesy and Geodynamics,2010,30(4):19-24.
[9] Barnhart W D,Lohman R B.Automated fault model discretization for inversions for coseismic slip distributions[J].Journal of Geophysical Research,2010,115(B10):B10419.
[10] Meade B J.Algorithms for the calculation of exact displacements,strains, and stresses for triangular dislocation elements in a uniform elastic half space[J].Computer and Geosciences,2007,33(8):1064-1075.
[11] Golub G H,Hansen P C,O'Leary D P.Tikhonov regularization and total least squares[J].SIAM Journal on Matrix Analysis and Applications,1999,21(1):185-194.
[12] Du Y J,Aydin A,Segall P.Comparison of various inversion techniques as applied to the determination of a geophysical deformation model for the 1983 Borah Peak earthquake[J].Bulletin of the Seismological Society of America,1992,82(4):1840-1866.
[13] 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.
[14] 凌勇,曾祺明,罗扬,等.巴姆地震变形场和应力场:I.用差分干涉雷达和Okada方法求解[J]. 岩石学报,2006,22(9):2367-2374. Ling Y,Zeng Q M,Luo Y,et al.Deformation and stress field of Bam earthquake:Ⅰ.InSAR and Okada calculation[J].Acta Petrologica Sinica,2006,22(9):2367-2374.
2016, Vol. 28 
