Decoherence recovery method of coseismic deformation field：Taking the 2009 L′Aquila MW6.3 earthquake as an example
-
摘要: 在合成孔径雷达干涉测量中，干涉相干是局部干涉条纹质量的重要评价标准，但是由于各种失相干源的影响，会出现失相干现象，而这可能会丢失重要的形变数据。本文首先利用ENVISAT ASAR的两景影像获得了2009年意大利拉奎拉MW6.3地震的带有失相干现象的同震形变场，对数据预处理后再使用移动窗口克里金法对拉奎拉地震的同震形变场进行插值处理，恢复了形变图的失相干区域，之后进一步通过正反演获得同震形变场，并与上一步的插值结果进行对比分析。结果显示，插值结果与反演结果在极震区吻合得很好，验证了移动窗口克里金插值法恢复失相干的可靠性，说明插值法可以成为实现地震同震形变场失相干恢复的一种途径。
-
关键词:
- InSAR技术 /
- 2009年拉奎拉MW6.3地震 /
- 移动窗口克里金法 /
- 正反演 /
- 失相干恢复
Abstract: In SAR interferometry, interference coherence is an important evaluation criterion for the quality of local interference fringes. However, due to the influence of various sources of decoherence, decoherent phenomena appear, which may lose important deformation data. In this paper, we first obtain the coseismic deformation field of the 2009 L′Aquila MW6.3 earthquake with the two scenes of ENVISAT ASAR. After preprocessing the data, we interpolate the coseismic deformation field of the L′Aquila earthquake using the moving window Kriging method, with a result that the decoherence area of the deformation map is restored, and then the co-seismic deformation field obtained by the forward and inversion is used to compare and analyze the interpolation results with the previous step. The results show that the interpolation results are in good agreement with the inversion results in the meizoseismal area, which verifies the reliability of the moving window Kriging interpolation method to restore the decoherence, and shows that the interpolation method can be used to realize the decoherence recovery of the coseismic deformation field. Therefore, it has important reference significance for the research on the recovery of remote sensing image data with decoherence phenomenon. -
-
![]()
图 3 拉奎拉地震同震形变场（a）及其部分放大区域（b）
Figure 3. Coseismic deformation field of the L′Aquila earthquake （a） and part of its magnified area （b）
![]()
图 5 经过裁剪、降采样处理所获得的结果形变图
Figure 5. Resulting deformation image obtained after the cropping and downsampling processing
表 1 ENVI ASAR雷达数据参数列表
Table 1 Data parameters of ENVI ASAR radar
轨道 升降轨模式 接收日期 时间基线/d 空间垂基线/m T079 降轨 2008-04-27
2009-04-12350 40.401 
下载: 导出CSV
表 2 拉奎拉地震断层模型几何参数列表
Table 2 List of geometric parameters of L′Aquila seismic fault model
采样点个数 断层参数 MW 误差/cm 东经/（°） 北纬/（°） 走向/° 倾角/° 断层长度/km 顶深/km 底深/km 滑动角/° 2 919 13.470 42.367 141.8° 55° 13.2 2.9 12.6 −118.7 6.3 1.2
下载: 导出CSV
-
冯万鹏,李振洪,李春来. 2010. 利用InSAR确定2009年4月6日MW6.3拉奎拉(Italy)地震最优震源模型[J]. 地球物理学进展,25(5):1550–1559. Feng W P,Li Z H,Li C L. 2010. Optimal source parameters of the 6 April 2009 MW6.3 L′Aquila,Italy earthquake from InSAR observations[J]. Progress in Geophysics,25(5):1550–1559 (in Chinese).
刘唐志,郭小宏. 2003. 均匀选点在全局优化中的应用[J]. 重庆交通学院学报,22(2):110–112. Liu T Z,Guo X H. 2003. The application of selecting oringal-points equally in global-optimization[J]. Journal of Chongqing Jiaotong University,22(2):110–112 (in Chinese).
刘晓萌,常占强,张景发,龚利霞. 2007. D-InSAR处理中失相干问题的研究[J]. 河北师范大学学报(自然科学版),31(2):260–263. Liu X M,Chang Z Q,Zhang J F,Gong L X. 2007. Investigation on decorrelation in processing of D-InSAR[J]. Journal of Hebei Normal University (Natural Science Edition)
,31(2):260–263 (in Chinese). 刘智荣,黄静宜,白相东,洪顺英,陈耀飞. 2012. 基于梯度下降法反演于田MW7.1地震同震滑动分布[J]. 大地测量与地球动力学,32(6):8–12. Liu Z R,Huang J Y,Bai X D,Hong S Y,Chen Y F. 2012. Coseismic slip distribution of Yutian MW7.1 earthquake based on steepest descend method[J]. Geodesy and Geodynamics,32(6):8–12 (in Chinese).
罗三明,杨国华,董运洪,李陶,袁油新,贾有良. 2012. 利用PSInSAR时间序列研究拉奎拉地震位移场变化特征[J]. 地球物理学进展,27(4):1307–1315. doi: 10.6038/j.issn.1004-2903.2012.04.004 Luo S M,Yang G H,Dong Y H,Li T,Yuan Y X,Jia Y L. 2012. Displacement field characteristics on L′Aquila earthquake from PSInSAR time series[J]. Progress in Geophysics,27(4):1307–1315 (in Chinese).
王超, 张红, 刘智. 2002. 星载合成孔径雷达干涉测量[M]. 北京: 科学出版社: 45−50. Wang C, Zhang H, Liu Z. 2002. Spaceborne Synthetic Aperture Radar Interferometry[M]. Beijing: Science Press: 45−50 （in Chinese）.
王乐洋,李海燕,温扬茂,许才军. 2017. 地震同震滑动分布反演的总体最小二乘方法[J]. 测绘学报,46(3):307–315. doi: 10.11947/j.AGCS.2017.20160212 Wang L Y,Li H Y,Wen Y M,Xu C J. 2017. Total least squares method inversion for coseismic slip distribution[J]. Acta Geodaetica et Cartographica Sinica,46(3):307–315 (in Chinese).
王秀萍. 2010. InSAR图像相位解缠的最小费用流法及其改进算法研究[J]. 测绘科学,35(4):129–131. Wang X P. 2010. Minimum cost flow and its ameliorative algorithm for phase unwrapping of InSAR image[J]. Science of Surveying and Mapping,35(4):129–131 (in Chinese).
王艳妮,谢金梅,郭祥. 2008. ArcGIS中的地统计克里格插值法及其应用[J]. 软件导刊,7(12):36–38. Wang Y N,Xie J M,Guo X. 2008. Application of geostatistical interpolation method in ArcGIS[J]. Software Guide,7(12):36–38 (in Chinese).
王永哲,李志伟,朱建军,胡俊. 2012. 融合多平台DInSAR数据解算拉奎拉地震三维同震形变场[J]. 武汉大学学报(信息科学版),37(7):859–863. Wang Y Z,Li Z W,Zhu J J,Hu J. 2012. Coseismic three-dimensional deformation of L′Aquila earthquake derived from multi-platform DInSAR data[J]. Geomatics and Information Science of Wuhan University,37(7):859–863 (in Chinese).
王子龙. 2011. 基于Surfer的GPS高程异常插值研究[J]. 矿山测量,(4):44–45. doi: 10.3969/j.issn.1001-358X.2011.04.015 Wang Z L. 2011. Research on GPS height anomaly interpolation based on Surfer[J]. Mine Survey,(4):44–45 (in Chinese).
杨久东,田辉,周春霞. 2010. 基于D-InSAR技术的意大利拉奎拉地震同震形变分析[J]. 大地测量与地球动力学,30(1):9–13. Yang J D,Tian H,Zhou C X. 2010. Coseismic deformation analysis of L′Aquila earthquake (Italy) by using D-InSAR[J]. Journal of Geodesy and Geodynamics,30(1):9–13 (in Chinese).
叶宝莹,沈利霞,郭芳芳,杨农. 2009. ArcGIS与Surfer图形文件的转换方法实现[J]. 测绘学科,34(5):212–213. Ye B Y,Shen L X,Guo F F,Yang N. 2009. The implementation of graphic transfer from Surfer to ArcGIS[J]. Science of Surveying and Mapping,34(5):212–213 (in Chinese).
张健, 彭军还. 2011. 基于D-InSAR技术的拉奎拉地震的同震形变初探[C]//2011成像雷达对地观测高级学术研讨会论文集. 北京: 国际数字地球学会中国国家委员会, 中国科学院: 378−383. Zhang J, Peng J H. 2011. Preliminary study on coseismic deformation of L′Aquila earthquake based on D-InSAR technology[C]//Proceedings of the 2011 Advanced Seminar on Imaging Radar Earth Observation. Beijing: Chinese National Committee of International Society for Digital Earth, Chinese Academy of Sciences: 378−383 （in Chinese）.
Goldstein R M,Werner C L. 1998. Radar interferogram filtering for geophysical applications[J]. Geophys Res Lett,25(21):4035–4038. doi: 10.1029/1998GL900033
Harris P,Charlton M,Fotheringham A S. 2010. Moving window Kriging with geographically weighted variograms[J]. Stoch Environ Res Risk Assess,24(8):1193–1209. doi: 10.1007/s00477-010-0391-2
Ryder I,Parsons B,Wright T J,Funning G J. 2010. Post-seismic motion following the 1997 Manyi (Tibet) earthquake:InSAR observations and modelling[J]. Geophys J Int,169(3):1009–1027.
Wang R J,Lorenzo-Martín F,Roth F. 2003. Computation of deformation induced by earthquakes in a multi-layered elastic crust:FORTRAN programs EDGRN/EDCMP[J]. Comput Geosci,29(2):195–207. doi: 10.1016/S0098-3004(02)00111-5
Wang R J,Lorenzo-Martín F,Roth F. 2006. PSGRN/PSCMP:A new code for calculating co-and post-seismic deformation,geoid and gravity changes based on the viscoelastic-gravitational dislocation theory[J]. Comput Geosci,32(4):527–541. doi: 10.1016/j.cageo.2005.08.006
Wang R J, Diao F Q, Hoechner A. 2013. SDM: A geodetic inversion code incorporating with layered crust structure and curved fault geometry[C]//EGU General Assembly Conference. Vienna, Austria: EGU.
Wang Z, Zhang J, Huang G, Chen T. 2010. Monitoring the coseismic deformation of L′Aquila earthquake using D-InSAR technique[C]//2010 2nd International Conference on Advanced Computer Control. Shenyang: IEEE: 521–525. doi: http://dx.doi.org/10.1109/ICACC.2010.5487119.
Zebker H A,Villasenor J. 1992. Decorrelation in interferometric radar echoes[J]. IEEE Trans Geosci Remote Sens,30(5):950–959. doi: 10.1109/36.175330
计量
- 文章访问数: 1217
- HTML全文浏览量: 447
- PDF下载量: 32
