陈明锴, 许光煜, 王乐洋. InSAR同震地表三维形变反演:一种顾及形变梯度的联合解算方法[J]. 武汉大学学报 ( 信息科学版), 2023, 48(8): 1349-1358. DOI: 10.13203/j.whugis20220284
引用本文: 陈明锴, 许光煜, 王乐洋. InSAR同震地表三维形变反演:一种顾及形变梯度的联合解算方法[J]. 武汉大学学报 ( 信息科学版), 2023, 48(8): 1349-1358. DOI: 10.13203/j.whugis20220284
CHEN Mingkai, XU Guangyu, WANG Leyang. InSAR 3D Coseismic Surface Deformation Inversion: A Combination Method Considering Deformation Gradient[J]. Geomatics and Information Science of Wuhan University, 2023, 48(8): 1349-1358. DOI: 10.13203/j.whugis20220284
Citation: CHEN Mingkai, XU Guangyu, WANG Leyang. InSAR 3D Coseismic Surface Deformation Inversion: A Combination Method Considering Deformation Gradient[J]. Geomatics and Information Science of Wuhan University, 2023, 48(8): 1349-1358. DOI: 10.13203/j.whugis20220284

InSAR同震地表三维形变反演:一种顾及形变梯度的联合解算方法

InSAR 3D Coseismic Surface Deformation Inversion: A Combination Method Considering Deformation Gradient

  • 摘要: 同震地表三维形变场是对地震引起的地表形变在真实空间中最直观的描述,在地震形变研究中,同震地表三维形变场的恢复具有重要意义。针对目前地表三维形变解算方法存在的不足和缺陷,提出了一种顾及形变梯度的同震地表三维形变场联合解算方法。该方法利用形变梯度信息将合成孔径雷达干涉(interferometric synthetic aperture radar,InSAR)观测的形变区域分解成形变梯度近似相等的互不相交的子区域,然后根据子区域内观测值数量选择加权最小二乘法(weighted least square,WLS)或基于应力-应变模型(stress-strain model,SM)和方差分量估计(variance components estimation,VCE)的InSAR三维形变解算方法(SM-VCE)获取同震地表三维形变场。模拟实验和实际震例研究结果表明,相比于传统的WLS和SM-VCE方法,通过形变梯度信息能有效恢复地震破裂带的地表三维形变,得到更加完整和可靠的同震地表三维形变场。

     

    Abstract:
      Objectives  The coseismic surface three-dimensional (3D) deformation field is the most intuitive description of the surface deformation caused by earthquakes in real space. In the study of seismic deformation, the recovery of coseismic surface 3D deformation field is of great significance. In view of the shortcomings and defects of the current coseismic 3D deformation calculation method, a coseismic surface 3D deformation field joint calculation method considering deformation gradient is proposed.
      Methods  First, we use the deformation gradient information to decompose the deformation region obtained by interferometric synthetic aperture radar (InSAR) into disjoint subregions with approximately equal deformation gradient. Then, according to the number of observations in the sub-region, the weighted least square (WLS) method or InSAR 3D deformation calculation method based on stress-strain model (SM) and variance component estimation (VCE) are selected to obtain 3D deformation field of the coseismic surface.
      Results  Compared with the traditional WLS and SM-VCE methods, the results of simulation experiments and actual earthquake examples show that the deformation gradient information can effectively restore the surface 3D deformation of the earthquake rupture zone and obtain more complete and reliable 3D deformation field of the coseismic surface.
      Conclusions  This proposed method can avoid the defect that the traditional SM-VCE method cannot identify the differential deformation information on both sides of the fault in the cross-fault area, and make up for the lack of adjacent point information in the traditional WLS solution method.

     

/

返回文章
返回