• ISSN 2097-1893
  • CN 10-1855/P

基于时序InSAR分析的高精度同震形变监测方法

吴雄骁 冯光财 贺礼家 卢昊

引用本文: 吴雄骁,冯光财,贺礼家,卢昊. 2023. 基于时序InSAR分析的高精度同震形变监测方法. 地球与行星物理论评(中英文),54(0):1-10
Wu X X, Feng G C, He L J, Lu H. 2023. High precision coseismic deformation monitoring method based on time-series InSAR analysis. Reviews of Geophysics and Planetary Physics, 54(0): 1-10 (in Chinese)

基于时序InSAR分析的高精度同震形变监测方法

doi: 10.19975/j.dqyxx.2022-023
基金项目: 国家自然科学基金资助项目(42174039,41574005)
详细信息
    作者简介:

    吴雄骁(1997-),男,硕士研究生,主要从事地震形变监测与参数反演工作.E-mail:wuxiongxiao0716@163.com

    通讯作者:

    冯光财(1981-),男,教授,主要从事高精度地表形变监测与参数反演方面的研究.E-mail:fredgps@csu.edu.cn

  • 中图分类号: P228,P315

High precision coseismic deformation monitoring method based on time-series InSAR analysis

Funds: Supported by the National Natural Science Foundation of China (Grant Nos. 42174039, 41574005)
  • 摘要: 合成孔径雷达干涉测量(InSAR)技术凭借全天时、全天候对地监测、高空间分辨率等特点,成为监测地表形变的重要手段,并广泛地应用到地震形变监测领域. 然而同震形变监测中最常用的D-InSAR技术在水域和植被覆盖严重等区域中容易受到时空失相关的影响,导致获取的同震形变场会受到严重的污染,此外还包含有大气延迟误差. 本文提出基于时序InSAR分析的高精度同震形变监测方法获取地震同震形变结果,主要是通过选择合适的干涉对和选择稳定点两步来提高形变场精度. 凭借充足的Sentinel-1A/B卫星SAR数据的支撑,利用大量震前和震后影像生成众多干涉图,按照一定标准挑选受误差影响较小的干涉图进行研究,减少大气延迟误差造成的影响;同时对震前影像的幅度图进行统计分析,从幅度值、相干性和幅度离散指数等方面设置阈值选择稳定点目标,削弱噪声干扰,提高形变场精度. 本文以2018年中国台湾花莲MW6.4地震为例,详细地介绍了高精度同震形变监测方法的数据处理流程,并与传统D-InSAR方法的结果进行了精度比较,结果表明本文方法能削弱形变场中的噪声误差,提高同震形变的信噪比. 应用本文方法获取了14个不同震级和位置的地震形变,结果表明通过选择稳定点的方式能提高形变场精度,且对于获取不同地震的同震形变场具有普遍适用性.

     

  • 图  1  高精度同震形变监测方法获得同震形变场的流程图. 整个过程包括输入、处理和输出三个部分,其中处理部分包括差分干涉图选择和稳定点选择两步

    Figure  1.  Flowchart of obtaining coseismic deformation by the proposed high-precision coseismic deformation monitoring method. The whole process includes three parts: input, processing and output. The processing part includes differential interferograms choosing and stable point selection

    图  2  地震案例研究位置分布图. 红色沙滩球是USGS给出的震源机制解,蓝色方框是每个地震对应的Sentinel-1A/B影像升降轨覆盖范围

    Figure  2.  The location distribution map of the selected earthquake. The red beach ball is the focal mechanism solution determined by USGS, and the blue box is coverage of the Sentinel-1A/B ascending and descending image for each earthquake

    图  3  2018年花莲地震升轨和降轨干涉图网,升轨和降轨各有33对干涉图,红色实线表示2018年地震发生的时刻

    Figure  3.  Interferogram network of the ascending and descending tracks of the 2018 Hualian earthquake. 33 pairs of interferograms for the ascending and descending tracks, and the red line indicates the occurrence time of the 2018 earthquake

    图  4  2018年花莲地震升轨33幅干涉图结果. 每个干涉图非变形区域的STD显示在子图的右下角,红色方框内是STD最小的干涉对

    Figure  4.  The 33 ascending interferogram pairs of the Hualien earthquake. The STD of the non-deformation area map of each interferogram is shown in the bottom right corner of the subfigure. The red box is the smallest interferogram pair of STD

    图  5  2018年花莲地震降轨33幅干涉图结果. 每个干涉图非变形区域的STD显示在子图的右下角,红色方框内是STD最小的干涉对

    Figure  5.  The 33 descending interferogram pairs of the Hualien earthquake. The STD of the non-deformation area map of each interferogram is shown in the bottom right corner of the subfigure. The red box is the smallest interferogram pair of STD

    图  6  Sentinel-1影像获取的2018年花莲地震的同震形变场. (a, b)用本文高精度同震形变监测方法得到的升轨和降轨形变结果;(e, f)传统D-InSAR获取的升轨和降轨结果;(c, d, g, h)分别对应图(a)中选取的A、B、C、D四个区域的光学图

    Figure  6.  The coseismic deformation fields of the 2018 Hualien earthquake derived from Sentinel-1 A/B images. (a, b) The ascending and descending interferograms obtained from high-precision coseismic deformation monitoring method, respectively; (e, f) The ascending and descending coseismic interferograms obtained by the classic D-InSAR, respectively; (c, d, g, h) Correspond to the optical images of the areas A, B, C and D selected in (a), respectively

    图  7  实验所用的14个地震案例的位置分布和同震形变结果

    Figure  7.  Location distribution and coseismic deformation results of 14 earthquake cases

    图  8  实验中14个地震案例升轨和降轨形变场的最大形变量及误差标准差

    Figure  8.  The maximum deformation and error standard deviations in the ascending and descending deformation fields of 14 seismic cases

    表  1  实验分析所用地震参数信息

    Table  1.   Seismic parameter information for experimental analysis

    编号日期参考位置经度/°E纬度/°N深度/km节面Ⅰ矩震级
    走向/(°)倾角/(°)滑动角/(°)
    1 2015-04-25 西藏定日县 87.3 28.4 20 - - - 5.9
    2 2015-07-03 新疆皮山县 78.2 37.6 10 317 69 102 6.5
    3 2016-01-21 青海门源县 101.62 37.68 10 337 41 103 6.4
    4 2016-11-25 新疆阿克陶县 74.04 39.27 10 199 84 14 6.7
    5 2018-02-06 台湾花莲县 121.71 24.13 11 209 73 22 6.4
    6 2018-12-24 西藏谢通门县 87.64 30.32 8 206 42 −66 5.8
    7 2020-01-19 新疆伽师县 77.21 39.83 16 221 20 72 6.4
    8 2020-03-20 西藏定日县 87.42 28.63 10 180 42 −77 5.9
    9 2020-06-26 新疆于田县 82.33 35.73 10 24 42 −108 6.3
    10 2020-07-23 西藏尼玛县 86.81 33.19 10 203 29 −88 6.6
    11 2021-03-19 西藏比如县 92.74 31.94 10 17 37 −113 6.1
    12 2021-03-30 西藏双湖县 87.68 34.38 10 173 55 −129 5.8
    13 2021-05-21 云南漾濞县 99.87 25.67 8 135 82 −165 6.4
    14 2021-05-22 青海玛多县 98.34 34.59 17 92 67 −40 7.4
    下载: 导出CSV

    表  2  2018年花莲地震Sentinel-1A/B卫星影像信息

    Table  2.   Sentinel-1A/B satellite image information of the 2018 Hualien earthquake

    卫星轨道震前影像时间震后影像时间干涉图数量
    Sentinel-
    1A/B
    升轨 2017-10-06 2017-10-18
    2017-10-30 2017-11-11
    2017-11-23 2017-12-05
    2017-12-17 2017-12-29
    2018-01-10 2018-01-22
    2018-02-03
    2018-02-09
    2018-02-15
    2018-02-27
    33
    降轨 2017-09-14 2017-09-26
    2017-10-08 2017-10-20
    2017-11-13 2017-11-25
    2017-12-07 2017-12-19
    2017-12-31 2018-01-12
    2018-02-05
    2018-02-11
    2018-02-17
    2018-03-01
    33
    下载: 导出CSV

    表  3  本文方法和传统D-InSAR方法所得形变场的精度比较

    Table  3.   The accuracy comparison of the deformation field obtained by the proposed method and the traditional D-InSAR method

    轨道方法STD/cm
    ABCD
    升轨 本文方法 0.91 0.83 0.89 0.52
    传统D-InSAR方法 1.05 0.84 1.01 0.56
    降轨 本文方法 0.86 0.72 1.08 0.80
    传统D-InSAR方法 1.34 0.75 1.2 0.88
    下载: 导出CSV
  • [1] Berardino P, Fornaro G, Lanari R, et al. 2002. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms[J]. IEEE Transactions on geoscience and remote sensing, 40(11): 2375-2383. doi: 10.1109/TGRS.2002.803792
    [2] Bie L, Ryder I, Nippress S E J, et al. 2014. Coseismic and post-seismic activity associated with the 2008 MW6.3 Damxung earthquake, Tibet, constrained by InSAR[J]. Geophysical Journal International, 196(2): 788-803. doi: 10.1093/gji/ggt444
    [3] Colesanti C, Ferretti A, Prati C, et al. 2003. Monitoring landslides and tectonic motions with the permanent scatterers technique[J]. Engineering geology, 68(1-2): 3-14. doi: 10.1016/S0013-7952(02)00195-3
    [4] Costantini M. 1998. A novel phase unwrapping method based on network programming[J]. IEEE Transactions on geoscience and remote sensing, 36(3): 813-821. doi: 10.1109/36.673674
    [5] Feng G C, Hetland E A, Ding X L, et al. 2010. Coseismic fault slip of the 2008 MW7.9 Wenchuan earthquake estimated from InSAR and GPS measurements[J]. Geophysical research letters, 37(1): L01302.
    [6] 冯光财, 许兵, 单新建, 等. 2015. 基于Landsat 8光学影像的巴基斯坦Awaran MW7.7地震形变监测及参数反演研究[J]. 地球物理学报, 58(5): 1634-1644 doi: 10.6038/cjg20150515

    Feng G C, Xu B, Shan X J, et al. 2015. Coseismic deformation and source parameters of the 24 September 2013 Awaran, Pakistan Mw7.7 earthquake derived from optical Landsat 8 satellite images[J]. Chinese Journal of Geophysics, 58(5): 1634-1644 (in Chinese). doi: 10.6038/cjg20150515
    [7] Ferretti A, Prati C, Rocca F. 2001. Permanent scatterers in SAR interferometry[J]. IEEE Transactions on geoscience and remote sensing, 39(1): 8-20. doi: 10.1109/36.898661
    [8] Fialko Y, Sandwell D, Simons M, et al. 2005. Three-dimensional deformation caused by the Bam, Iran, earthquake and the origin of shallow slip deficit[J]. Nature, 435(7040): 295-299. doi: 10.1038/nature03425
    [9] Gabriel A K, Goldstein R M, Zebker H A. 1989. Mapping small elevation changes over large areas: Differential radar interferometry[J]. Journal of Geophysical Research: Solid Earth, 94(B7): 9183-9191. doi: 10.1029/JB094iB07p09183
    [10] Goldstein R M, Zebker H A, Werner C L. 1988. Satellite radar interferometry: Two-dimensional phase unwrapping [J]. Radio science, 23(4): 713-720. doi: 10.1029/RS023i004p00713
    [11] Goldstein R M, Werner C L. 1998. Radar interferogram filtering for geophysical applications[J]. Geophysical research letters, 25(21): 4035-4038. doi: 10.1029/1998GL900033
    [12] Hanssen R F. 2001. Radar Interferometry: Data Interpretation and Error Analysis[M]. Springer Science & Business Media.
    [13] Hooper A, Segall P, Zebker H. 2007. Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos[J]. Journal of Geophysical Research: Solid Earth, 112: B07407.
    [14] Li S W, Xu W B, Li Z W. 2022. Review of the SBAS InSAR time-series algorithms, applications, and challenges[J]. Geodesy and Geodynamics, 13: 114-126. doi: 10.1016/j.geog.2021.09.007
    [15] Li Z W, Ding X L, Liu G X. 2004. Modeling atmospheric effects on InSAR with meteorological and continuous GPS observations: algorithms and some test results[J]. Journal of Atmospheric and Solar-Terrestrial Physics, 66(11): 907-917. doi: 10.1016/j.jastp.2004.02.006
    [16] Massonnet D, Rossi M, Carmona C, et al. 1993. The displacement field of the Landers earthquake mapped by radar interferometry[J]. Nature, 364(6433): 138-142. doi: 10.1038/364138a0
    [17] 欧阳伦曦, 李新情, 惠凤鸣, 等. 2017. 哨兵卫星Sentinel-1A数据特性及应用潜力分析[J]. 极地研究, 29(2): 286-295

    Ouyang L X, Li X Q, Hui F M, et al. 2017. Sentinel-1A data characteristics and applications potential analysis[J]. Chinese Journal of Polar Research, 29(2): 286-295 (in Chinese).
    [18] Rigo A, De Chabalier J B, Meyer B, et al. 2004. The 1995 Kozani-Grevena (northern Greece) earthquake revisited: An improved faulting model from synthetic aperture radar interferometry[J]. Geophysical Journal International, 157(2): 727-736. doi: 10.1111/j.1365-246X.2004.02220.x
    [19] Rosen P A, Hensley S, Joughin I R, et al. 2000. Synthetic aperture radar interferometry[J]. Proceedings of the IEEE, 88(3): 333-382. doi: 10.1109/5.838084
    [20] Sandwell D T, Price E J. 1998. Phase gradient approach to stacking interferograms[J]. Journal of Geophysical Research: Solid Earth, 103(B12): 30183-30204. doi: 10.1029/1998JB900008
    [21] 舒宁. 2003. 雷达影像干涉测量原理[M]. 武汉: 武汉大学出版社.

    Shu N. 2003. Principle of Radar Image Interferometry[M]. Wuhan: Wuhan University Press (in Chinese).
    [22] Simons M, Fialko Y, Rivera L. 2002. Coseismic deformation from the 1999 MW7.1 Hector Mine, California, earthquake as inferred from InSAR and GPS observations[J]. Bulletin of the Seismological Society of America, 92(4): 1390-1402. doi: 10.1785/0120000933
    [23] Wan Y G, Shen Z -K, Bürgmann R, et al. 2017. Fault geometry and slip distribution of the 2008 MW7.9 Wenchuan, China earthquake, inferred from GPS and InSAR measurements[J]. Geophysical Journal International, 208: 748-766. doi: 10.1093/gji/ggw421
    [24] Werner C, Wegmüller U, Strozzi T, et al. 2000. Gamma SAR and interferometric processing software[C]//Proceedings of the Ers-envisat Symposium, Gothenburg, Sweden. 1620: 1620.
    [25] Wright T J, Parsons B, England P C, et al. 2004. InSAR observations of low slip rates on the major faults of western Tibet[J]. Science, 305(5681): 236-239. doi: 10.1126/science.1096388
    [26] Xu W, Wu S, Materna K, et al. 2018. Interseismic ground deformation and fault slip rates in the greater San Francisco bay area from two decades of space geodetic data[J]. Journal of Geophysical Research: Solid Earth, 123(9): 8095-8109. doi: 10.1029/2018JB016004
    [27] 杨魁, 杨建兵, 江冰茹. 2015. Sentinel-1卫星综述[J]. 城市勘测, 24-27 doi: 10.3969/j.issn.1672-8262.2015.02.006

    Yang K, Yang J B, Jiang B R. 2015. Sentinel-1 Satellite Overview[J]. Urban Geotechnical Investigation & Surveying, 24-27 (in Chinese). doi: 10.3969/j.issn.1672-8262.2015.02.006
    [28] Zebker H A, Villasenor J. 1992. Decorrelation in interferometric radar echoes [J]. IEEE Transactions on Geoscience and Remote Sensing, 30(5): 950-959. doi: 10.1109/36.175330
  • 加载中
图(8) / 表(3)
计量
  • 文章访问数:  375
  • HTML全文浏览量:  192
  • PDF下载量:  99
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-03-19
  • 录用日期:  2022-04-12
  • 修回日期:  2022-04-12
  • 网络出版日期:  2022-04-24

目录

    /

    返回文章
    返回