留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

尼泊尔地震触发滑坡识别和雪崩形变分析

丁超 冯光财 周玉杉 王会强 杜亚男 陈晨月

丁超, 冯光财, 周玉杉, 王会强, 杜亚男, 陈晨月. 尼泊尔地震触发滑坡识别和雪崩形变分析[J]. 武汉大学学报 ( 信息科学版), 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
引用本文: 丁超, 冯光财, 周玉杉, 王会强, 杜亚男, 陈晨月. 尼泊尔地震触发滑坡识别和雪崩形变分析[J]. 武汉大学学报 ( 信息科学版), 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
DING Chao, FENG Guangcai, ZHOU Yushan, WANG Huiqiang, DU Ya'nan, CHEN Chenyue. Nepal Earthquake Triggered Landslides Recognition and Deformation Analysis of Avalanches' Region[J]. Geomatics and Information Science of Wuhan University, 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
Citation: DING Chao, FENG Guangcai, ZHOU Yushan, WANG Huiqiang, DU Ya'nan, CHEN Chenyue. Nepal Earthquake Triggered Landslides Recognition and Deformation Analysis of Avalanches' Region[J]. Geomatics and Information Science of Wuhan University, 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031

尼泊尔地震触发滑坡识别和雪崩形变分析

doi: 10.13203/j.whugis20160031
基金项目: 

地震动力学国家重点实验室开放基金 LED2014B02

详细信息

Nepal Earthquake Triggered Landslides Recognition and Deformation Analysis of Avalanches' Region

Funds: 

the Open Foundation of the State Key Laboratory of Earthquake Dynamics LED2014B02

More Information
    Author Bio:

    DING Chao, postgraduate, specializes in the earthquake deformation monitoring. E-mail: dingchao_csu@csu.edu.cn

    Corresponding author: FENG Guangcai, PhD, associate professor. E-mail: fredgps@csu.edu.cn
  • 摘要: 2015-04-25尼泊尔Gorkha附近发生了Mw 7.8级强震,引发尼泊尔及中国西藏地区大量的山体滑坡和雪崩灾害,造成了严重的人员伤亡与经济损失。基于Landsat 8影像和亚像素相关性匹配技术识别形变失相关区域,并结合目视解译提取震害滑坡与雪崩沉积区域,统计分析了所提取的滑坡数目、面积等与地形因子的关系,同时对珠峰雪崩高发区的冰川形变进行了定量分析。实验结果表明,提取形变场中的形变失相关区域并结合目视解译验证是一种识别震害滑坡与雪崩沉积区域的有效方法;统计结果显示滑坡体随坡度分布呈现高斯特性,多发生在坡度大于30°的地形,易发于河谷两侧,且高坡度更容易触发单体面积大的滑坡体。该研究为利用Landsat 8影像数据进行震害滑坡与雪崩沉积区域提取、雪崩高发区的冰川流速监测预警工作提供了新视角。
  • 图  1  震区构造背景与Landsat 8影像覆盖图

    Figure  1.  Background Setting and Coverage of Landsat 8 for the 2015 Nepal Earthquake

    图  2  吉隆形变场和滑坡体提取图

    Figure  2.  Extraction of Deformation Fields and Landslides of Jilong

    图  3  吉隆滑坡分布图

    Figure  3.  Landslides Distribution of Jilong

    图  4  聂拉木滑坡分布图

    Figure  4.  Landslides Distribution of Nyalam

    图  5  吉隆、聂拉木地区不同坡度水平滑坡体数目、总面积与平均面积分布直方图

    Figure  5.  Landslides' Quantity, Gross and Average Area Distribution Histogram of Jilong and Nyalam in Different Slope Levels

    图  6  郎塘形变场和雪崩沉积区提取图

    Figure  6.  Extraction of Deformation Fields and Avalanches' Debris Areas in Langtang

    图  7  珠峰研究区域及形变场

    Figure  7.  Everest Mount Region and Its Deformation Fields

    图  8  昆布冰川区影像与形变场

    Figure  8.  Image and Deformation Fields of the Khumbu Glacier

    图  9  东绒布冰川区影像与形变场

    Figure  9.  Image and Deformation Fields of the East Rongbuk Glacier

    表  1  本文实验所采用的Landsat 8影像信息

    Table  1.   Information of Landsat 8 Images Used in This Study

    影像对 时段 轨道号 分幅号 获取时间 云量/%
    影像对1 震前 141 41 2014-03-26 2.4
    震后 141 41 2015-04-30 30.2
    影像对2 震前 141 40 2014-01-05 4.4
    震后 141 40 2015-04-30 30.2
    影像对3 震前 140 41 2015-03-22 11.8
    震后 140 41 2015-05-25 25.6
    影像对4 震前 140 41 2014-03-19 17.4
    震后 140 41 2015-05-22 13.1
    下载: 导出CSV

    表  2  吉隆山体滑坡编目表

    Table  2.   Landslides Catalog of Jilong

    序号 滑坡位置 滑坡体面积/m2 技术监测手段 高程/m 坡向/(°) 坡度/(°)
    东经/(°) 北纬/(°) SPC 目视判读
    J01 85.363 28.360 163 646 2 615 289.6 52.8
    J02 85.359 28.354 29 052 2 370 268.7 43.5
    J03 85.337 28.376 37 580 2 490 80.4 30.2
    J04 85.324 28.389 49 045 2 737 312.1 45.9
    J05 85.305 28.401 5 599 2 837 359.1 37.3
    J06 85.304 28.406 33 891 2 809 28.9 40.3
    J07 85.350 28.364 2 603 2 543 10.6 47.5
    J08 85.345 28.369 11 552 2 415 165.9 38.9
    J09 85.337 28.374 5 967 2 450 102.1 33.3
    J10 85.340 28.381 7 431 线状滑坡× 2 839 208.5 53.9
    J11 85.337 28.374 1 295 小滑坡× 2 603 24.7 27.1
    J12 85.355 28.364 680 小滑坡× 2 472 120.0 45.7
    J13 85.370 28.376 6 789 2 925 212.4 44.8
    J14 85.358 28.370 42 480 2 439 196.8 52.6
    J15 85.355 28.372 1 057 小滑坡× 2 550 172.0 41.4
    J16 85.308 28.404 2 263 小滑坡× 2 804 191.8 33.3
    J17 85.303 28.408 1 728 小滑坡× 2 805 49.1 42.9
    J18 85.312 28.399 3 628 小滑坡× 2 741 172.1 25.6
    J19 85.358 28.411 351 小滑坡× 2 984 135.1 37.0
    J20 85.364 28.382 21 275 3 152 256.7 48.2
    J21 85.365 28.382 17 310 3 152 256.7 42.7
    注:√表示可以监测到对应滑坡体,×表示无法监测到对应滑坡体。
    下载: 导出CSV

    表  3  聂拉木山体滑坡编目表

    Table  3.   Landslides Catalog of Nyalam

    序号 滑坡位置 滑坡体面积/m2 高程/m 坡向/(°) 坡度/(°)
    东经/(°) 北纬/(°)
    N01 85.997 28.090 46 565 3 457 123.3 37.6
    N02 85.998 28.087 7 153 3 264 177.7 31.7
    N03 85.990 28.084 13 099 3 769 166.3 60.8
    N04 85.997 28.085 6 360 3 238 270 26.9
    N05 85.991 28.046 38 639 3 265 274.9 57.1
    N06 85.994 28.039 53 518 3 049 171.8 32.8
    N07 85.981 28.032 51 028 3 157 100.9 68.8
    N08 85.980 28.021 12 191 2 602 106.2 45.4
    N09 85.963 27.983 309 600 2 375 157.2 60.8
    N10 86.039 27.975 381 165 3 462 106.5 42.2
    下载: 导出CSV
  • [1] Wikipedia. April 2015 Nepal Earthquake[OL]. https://en.wikipedia.org/wiki/April 2015 Nepal earthquake, 2016
    [2] Feng G C, Li Z W, Shan X J, et al. Geodetic Model of the 2015 April 25 Mw 7.8 Gorkha Nepal Earthquake and Mw 7.3 Aftershock Estimated from InSAR and GPS Data[J]. Geophysical Journal International, 2015, 203(2):896-900 doi:  10.1093/gji/ggv335
    [3] 陶舒. 汶川地震滑坡遥感信息提取及灾害危险性评价研究[D]. 北京: 首都师范大学, 2009

    Tao Shu. Information Extraction of Landslides Induced by the Wenchuan Earthquake and Susceptibi-lity Assessment Using Remotely Sensed Data[D]. Beijing: Capital Normal University, 2009
    [4] 许冲, 徐锡伟, 沈玲玲, 等. 2014年鲁甸MS 6.5地震触发滑坡编录及其对一些地震参数的指示[J].地震地质, 2014, 36(4):1186-1203 https://www.researchgate.net/profile/Chong_Xu/publication/287556048_Inventory_of_landslides_triggered_by_the_2014_MS65_Ludian_earthquake_and_its_implications_on_several_earthquake_parameters/links/580ea47a08ae7525273d29d3.pdf?origin=publication_list

    Xu Chong, Xu Xiwei, Shen Lingling, et al. Inventory of Landslides Triggered by the 2014 MS 6.5 Ludian Earthquake and Its Implications on Several Earthquake Parameters[J].Seismology and Geology, 2014, 36(4):1186-1203 https://www.researchgate.net/profile/Chong_Xu/publication/287556048_Inventory_of_landslides_triggered_by_the_2014_MS65_Ludian_earthquake_and_its_implications_on_several_earthquake_parameters/links/580ea47a08ae7525273d29d3.pdf?origin=publication_list
    [5] 武雪玲, 沈少青, 牛瑞卿. GIS支持下应用PSO-SVM模型预测滑坡易发性[J].武汉大学学报·信息科学版, 2016, 41(5):665-671 http://ch.whu.edu.cn/CN/abstract/abstract5447.shtml

    Wu Xueling, Shen Shaoqing, Niu Ruiqing. Landslide Susceptibility Prediction Using GIS and PSO-SVM[J]. Geomatics and Information Science of Wuhan University, 2016, 41(5):665-671 http://ch.whu.edu.cn/CN/abstract/abstract5447.shtml
    [6] Lacroix P, Zavala B, Berthier E, et al. Supervised Method of Landslide Inventory Using Panchromatic SPOT5 Images and Application to the Earthquake-Triggered Landslides of Pisco (Peru, 2007, Mw8.0)[J]. Remote Sensing, 2013, 5(6):2590-2616 doi:  10.3390/rs5062590
    [7] Leprince S, Barbot S, Ayoub F, et al. Automatic and Precise Orthorectification, Co-registration and Sub-pixel Correlation of Satellite Images, Application to Ground Deformation Measurements[J]. IEEE Transactions on Geoscience and Remote Sen-sing, 2007, 45(6):1529-1558 doi:  10.1109/TGRS.2006.888937
    [8] 董晓燕, 丁晓利, 李志伟, 等.一种新的SAR像素偏移量估计流程及其在同震形变监测中的应用[J].武汉大学学报·信息科学版, 2011, 36(7):789-792 http://ch.whu.edu.cn/CN/abstract/abstract592.shtml

    Dong Xiaoyan, Ding Xiaoli, Li Zhiwei, et al. A New Chain for SAR Pixel Offset Eestimation and Its Application in Co-seismic Deformation Measurement[J].Geomatics and Information Science of Wuhan University, 2011, 36(7):789-792 http://ch.whu.edu.cn/CN/abstract/abstract592.shtml
    [9] Ding C, Feng G C, Li Z W, et al. Spatio-Temporal Error Sources Analysis and Accuracy Improvement in Landsat 8 Image Ground Displacement Measurements[J].Remote Sensing, 2016, 8(11):937-958 doi:  10.3390/rs8110937
    [10] Binet R, Bollinger L. Horizontal Coseismic Deformation of the 2003 Bam (Iran) Earthquake Measured from SPOT-5 THR Satellite Imagery[J]. Geophysi-cal Research Letters, 2005, 320(2):287-294 doi:  10.1029/2004GL021897
    [11] Leung T K P, Admas J. The 2005, Mw 7.6 Kashmir Earthquake, Rupture Kinematics from Sub-pixel Correlation of ASTER Images and Seismic Waveforms Analysis[J]. Earth and Planetary Science Letters, 2006, 249(3):514-528 https://thesis.library.caltech.edu/2257/
    [12] 冯光财, 许兵, 单新建, 等.基于Landsat 8光学影像的巴基斯坦Awaran Mw7.7地震形变监测及参数反演研究[J].地球物理学报, 2015, 58(5):1634-1644 doi:  10.6038/cjg20150515

    Feng Guangcai, Xu Bing, Shan Xinjian, et al. 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 Geophysics, 2015, 58(5):1634-1644 doi:  10.6038/cjg20150515
    [13] Avouac J P, Meng L, Wei S, et al. Lower Edge of Locked Main Himalayan Thrust Unzipped by the 2015 Gorkha Earthquake[J].Nature Geoscience, 2015, 8(9):1-6 http://palgrave.nature.com/ngeo/journal/v8/n9/abs/ngeo2518.html
    [14] Kargel J S, Leonard G J, Shugar D H, et al. Geomorphic and Geologic Controls of Geohazards Induced by Nepal's 2015 Gorkha Earthquake[J]. Science, 2016, 351(6269):1-12 http://ecite.utas.edu.au/107060
    [15] Jane Q. Killer Landslides: The Lasting Legacy of Nepal's Quake[OL]. http://www.nature.com/news/killer-landslides-the-lasting-legacy-of-nepal-s-quake-1.19803, 2016
    [16] Scherler D, Leprince S, Strecker M R. Glacier-surface Velocities in Alpine Terrain from Optical Satellite Imagery-Accuracy Improvement and Quality Assessment[J]. Remote Sensing of Environment, 2008, 112(10):3806-3819 doi:  10.1016/j.rse.2008.05.018
  • [1] 王晨辉, 郭伟, 孟庆佳, 刘炎炎, 毕逢东.  基于虚拟参考站的GNSS滑坡变形监测方法及性能分析 . 武汉大学学报 ( 信息科学版), 2022, 47(6): 990-996. doi: 10.13203/j.whugis20220102
    [2] 何朝阳, 许强, 巨能攀, 解明礼.  滑坡实时监测预警模型调度算法优化研究 . 武汉大学学报 ( 信息科学版), 2021, 46(7): 970-982. doi: 10.13203/j.whugis20200314
    [3] 武雪玲, 杨经宇, 牛瑞卿.  一种结合SMOTE和卷积神经网络的滑坡易发性评价方法 . 武汉大学学报 ( 信息科学版), 2020, 45(8): 1223-1232. doi: 10.13203/j.whugis20200127
    [4] 秦宏楠, 马海涛, 于正兴.  地基SAR技术支持下的滑坡预警预报分析方法 . 武汉大学学报 ( 信息科学版), 2020, 45(11): 1697-1706. doi: 10.13203/j.whugis20200268
    [5] 陈涛, 钟子颖, 牛瑞卿, 刘桐, 陈胜云.  利用深度信念网络进行滑坡易发性评价 . 武汉大学学报 ( 信息科学版), 2020, 45(11): 1809-1817. doi: 10.13203/j.whugis20190144
    [6] 冯非凡, 武雪玲, 牛瑞卿, 许石罗.  一种V/S和LSTM结合的滑坡变形分析方法 . 武汉大学学报 ( 信息科学版), 2019, 44(5): 784-790. doi: 10.13203/j.whugis20170218
    [7] 刘国祥, 张波, 张瑞, 蔡嘉伦, 符茵, 刘巧, 于冰, 李志林.  联合卫星SAR和地基SAR的海螺沟冰川动态变化及次生滑坡灾害监测 . 武汉大学学报 ( 信息科学版), 2019, 44(7): 980-995. doi: 10.13203/j.whugis20190077
    [8] 张路, 廖明生, 董杰, 许强, 龚健雅.  基于时间序列InSAR分析的西部山区滑坡灾害隐患早期识别——以四川丹巴为例 . 武汉大学学报 ( 信息科学版), 2018, 43(12): 2039-2049. doi: 10.13203/j.whugis20180181
    [9] 刘渊博, 牛瑞卿, 于宪煜, 张凯翔.  旋转森林模型在滑坡易发性评价中的应用研究 . 武汉大学学报 ( 信息科学版), 2018, 43(6): 959-964. doi: 10.13203/j.whugis20160132
    [10] 刘坚, 李树林, 陈涛.  基于优化随机森林模型的滑坡易发性评价 . 武汉大学学报 ( 信息科学版), 2018, 43(7): 1085-1091. doi: 10.13203/j.whugis20160515
    [11] 段功豪, 牛瑞卿, 彭令, 付杰.  诱发因素影响下的滑坡参数优化预测模型研究 . 武汉大学学报 ( 信息科学版), 2017, 42(4): 531-536. doi: 10.13203/j.whugis20140913
    [12] 段功豪, 牛瑞卿, 赵艳南, 张凯翔, 咬登魁.  基于动态指数平滑模型的降雨诱发型滑坡预测 . 武汉大学学报 ( 信息科学版), 2016, 41(7): 958-962. doi: 10.13203/j.whugis20140276
    [13] 武雪玲, 沈少青, 牛瑞卿.  GIS支持下应用PSO-SVM模型预测滑坡易发性 . 武汉大学学报 ( 信息科学版), 2016, 41(5): 665-671. doi: 10.13203/j.whugis20130566
    [14] 彭令, 牛瑞卿, 赵艳南, 邓清禄.  基于核主成分分析和粒子群优化支持向量机的滑坡位移预测 . 武汉大学学报 ( 信息科学版), 2013, 38(2): 148-152,161.
    [15] 武雪玲, 任福, 牛瑞卿.  多源数据支持下的三峡库区滑坡灾害空间智能预测 . 武汉大学学报 ( 信息科学版), 2013, 38(8): 963-968.
    [16] 王利, 张勤, 管建安, 孔令杰.  基于GPS技术的滑坡动态变形监测试验结果与分析 . 武汉大学学报 ( 信息科学版), 2011, 36(4): 422-426.
    [17] 张正禄, 汪宏晨, 邓勇, 谢年生.  滑坡变形分析与预报的新方法 . 武汉大学学报 ( 信息科学版), 2009, 34(12): 1387-1389.
    [18] 张昆, 张松林, 刘祖强, 杨红.  滑坡变形的三维可视化研究 . 武汉大学学报 ( 信息科学版), 2006, 31(9): 795-798.
    [19] 曾旭平.  GPS滑坡高程监测的数据处理问题 . 武汉大学学报 ( 信息科学版), 2004, 29(3): 201-204.
    [20] 秦宏楠, 马海涛, 于正兴, 刘玉溪.  地基雷达干涉测量动态高频次数据用于滑坡早期预警方法研究 . 武汉大学学报 ( 信息科学版), 0, 0(0): -. doi: 10.13203/j.whugis20220152
  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  1061
  • HTML全文浏览量:  85
  • PDF下载量:  304
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-01-28
  • 刊出日期:  2018-06-05

尼泊尔地震触发滑坡识别和雪崩形变分析

doi: 10.13203/j.whugis20160031
    基金项目:

    地震动力学国家重点实验室开放基金 LED2014B02

    作者简介:

    丁超, 硕士生, 主要从事地震形变监测研究。dingchao_csu@csu.edu.cn

    通讯作者: 冯光财, 博士, 副教授。fredgps@csu.edu.cn
  • 中图分类号: P315

摘要: 2015-04-25尼泊尔Gorkha附近发生了Mw 7.8级强震,引发尼泊尔及中国西藏地区大量的山体滑坡和雪崩灾害,造成了严重的人员伤亡与经济损失。基于Landsat 8影像和亚像素相关性匹配技术识别形变失相关区域,并结合目视解译提取震害滑坡与雪崩沉积区域,统计分析了所提取的滑坡数目、面积等与地形因子的关系,同时对珠峰雪崩高发区的冰川形变进行了定量分析。实验结果表明,提取形变场中的形变失相关区域并结合目视解译验证是一种识别震害滑坡与雪崩沉积区域的有效方法;统计结果显示滑坡体随坡度分布呈现高斯特性,多发生在坡度大于30°的地形,易发于河谷两侧,且高坡度更容易触发单体面积大的滑坡体。该研究为利用Landsat 8影像数据进行震害滑坡与雪崩沉积区域提取、雪崩高发区的冰川流速监测预警工作提供了新视角。

English Abstract

丁超, 冯光财, 周玉杉, 王会强, 杜亚男, 陈晨月. 尼泊尔地震触发滑坡识别和雪崩形变分析[J]. 武汉大学学报 ( 信息科学版), 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
引用本文: 丁超, 冯光财, 周玉杉, 王会强, 杜亚男, 陈晨月. 尼泊尔地震触发滑坡识别和雪崩形变分析[J]. 武汉大学学报 ( 信息科学版), 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
DING Chao, FENG Guangcai, ZHOU Yushan, WANG Huiqiang, DU Ya'nan, CHEN Chenyue. Nepal Earthquake Triggered Landslides Recognition and Deformation Analysis of Avalanches' Region[J]. Geomatics and Information Science of Wuhan University, 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
Citation: DING Chao, FENG Guangcai, ZHOU Yushan, WANG Huiqiang, DU Ya'nan, CHEN Chenyue. Nepal Earthquake Triggered Landslides Recognition and Deformation Analysis of Avalanches' Region[J]. Geomatics and Information Science of Wuhan University, 2018, 43(6): 847-853, 950. doi: 10.13203/j.whugis20160031
  • 2015-04-25T14:11:00(北京时间)尼泊尔Gorkha附近发生了Mw 7.8级强震。截止2015-05-28,尼泊尔及周边国家已有8 964人丧生,21 952人受伤,300余万人口无家可归[1]。该地震所触发的大量山体滑坡、雪崩等次生灾害已造成了严重的人员伤亡与经济损失[1-2]

    快速实现对震后山体滑坡、雪崩分布的监测对震后救灾和次生灾害的预防非常关键,遥感技术是滑坡等灾害监测的重要手段[3-6]。目前,高分辨率光学影像已被广泛应用于震后滑坡的监测研究中,主要采取目视解译方法[3]。该方法可以准确提取滑坡信息,为灾害救援工作提供详实的技术支持,已成功应用于多个地震案例[3-4],但需要判读者具备丰富的遥感影像判读经验。除了目视解译的方法,遥感领域也发展了基于多波段分类与差值等自动化判读方法[3, 6]。近年来发展起来的影像亚像素相关性匹配技术[7-9](sub-pixel correlation,SPC)已广泛应用于地震形变场的提取,如2003年Bam地震[10]、2005年Kashmir地震[11]和2013年的巴基斯坦地震[9, 12]等。该技术的研究主要侧重形变信息的获取,失相关信息一般只作为噪声进行掩膜清除。由于尼泊尔震区特殊的地质和地势特点,本次震害滑坡、雪崩造成震害滑坡分布广泛,形变场存在较大比例的失相关噪声。基于研究现状与震区实际情况,本文以Landsat 8影像为主,综合遥感影像目视解译与SPC技术探讨从形变场失相关信息中提取震害滑坡与雪崩沉积区域的流程方法。

    • Landsat 8由美国国家航空航天局(NASA)在2013-02-11发射升空,全名为地球资源卫星数据连续性任务(Landsat Data Continuity Mission,LDCM)。与之前的系列卫星相比,Landsat 8的全色波段(Band 8)的空间分辨率更高(15 m),能更好地监测地表地物变化,在气象条件允许的情况下非常适合大范围形变监测[9, 12]。此外,Landsat 8影像具有丰富的多波谱信息,通过合成多波段影像,可凸显出相应的地表变化,有利于震害滑坡的识别。本文收集了覆盖地震区域的震前和震后的Landsat 8数据4对,轨道号为141、140;分幅号为40、41(如表 1图 1所示)。由于震后影像普遍受云层影响严重,本文选择了含云量较少的震前影像,以保证震前、震后数据的相关性。USGS提供的Landsat 8数据是已经做过正射校正的L1T产品,校正精度可满足SPC技术要求。

      表 1  本文实验所采用的Landsat 8影像信息

      Table 1.  Information of Landsat 8 Images Used in This Study

      影像对 时段 轨道号 分幅号 获取时间 云量/%
      影像对1 震前 141 41 2014-03-26 2.4
      震后 141 41 2015-04-30 30.2
      影像对2 震前 141 40 2014-01-05 4.4
      震后 141 40 2015-04-30 30.2
      影像对3 震前 140 41 2015-03-22 11.8
      震后 140 41 2015-05-25 25.6
      影像对4 震前 140 41 2014-03-19 17.4
      震后 140 41 2015-05-22 13.1

      图  1  震区构造背景与Landsat 8影像覆盖图

      Figure 1.  Background Setting and Coverage of Landsat 8 for the 2015 Nepal Earthquake

      图 1所示,右上方小图浅蓝色方框为Landsat 8影像覆盖区域,红色方框为大图所显示的区域,下方大图中洋红色方框为本文所选取的吉隆、郎塘、聂拉木和珠峰4个实验区,红色五角星分别显示的是2015-04-26 Mw 6.7级余震与2015-05-12 Mw 7.3级余震的震中位置,绿色点表示地震一个月内的余震分布。

    • 本文采用的SPC技术是基于窗口的影像纹理匹配,任何因素导致的纹理破坏、模糊、缺失都会造成匹配失败,进而产生失相关噪声。对该地震而言,本文将失相关噪声分为3类:一是时相失相关,例如云层、植被生长变化、人工设施拆建以及太阳角度变化在地形复杂区域所引起的地表辐射量变化等;二是辐射失相关,例如白雪、沙地、严重的地形阴影与大面积水域等辐射量过高或过低或为常数,造成匹配算法无法识别影像纹理;三是形变失相关,如震害滑坡、雪崩引起地表剧烈变化,造成影像纹理缺失或变化,进而无法匹配引入失相关噪声。根据相关研究[13-14],尼泊尔强震在震源中心附近造成2 m左右的水平同震形变。故本文假定,SPC技术获取的形变量绝对值小于等于2 m,认为是可靠测量值;否则认为是离群值。本文将离群值分为趋势性离群值与随机性离群值。趋势性离群值排除形变误差区域,可获取反映真实地表形变趋势的信息,可用于冰川形变提取;随机性离群值由SPC算法在匹配失败后对计算窗口的随机赋值得到,即失相关噪声,多对应低信噪比(signal noise ratio,SNR),本文从随机性离群值与低SNR中提取滑坡与雪崩沉积区域。利用震前、震后影像进行亚像素相关性匹配,提取震害滑坡体、雪崩沉积区域与雪崩区形变的具体流程如下。

      1) 影像预处理。选取震前、震后覆盖相同区域的影像做几何校正,并圈定云、阴影区域,以分析形变场中的失相关区域。

      2) 计算影像的偏移量。利用频域相关性算法对震前、震后的数据进行SPC计算,获取偏移量;窗口大小为32像素×32像素,步长为4像素,稳健迭代2次,频率掩膜阈值为0.9。

      3) 后处理。经过相关性计算后的结果中包含着许多误差,利用一阶多项式去除形变图中轨道趋势误差,通过均值滤波进一步降低形变图中其他噪声水平。

      4) 失相关区域。圈定同震形变场中的失相关区域,即低SNR与随机性离群值区域。

      5) 形变失相关区域。本文将失相关区域减去已确定的噪声失相关区域,初步确定形变失相关区域,再利用多源遥感影像目视解译方法对这些区域进行验证,进一步确定震害滑坡与雪崩沉积区域。

      6) 雪崩区形变。从同震形变场中的趋势性离群值中提取雪崩区形变场并进行分析。

    • 本文选取中国的吉隆和聂拉木两个典型滑坡受灾区,以及尼泊尔受雪崩危害严重的郎塘,利用§1.2方法提取出实验区内震害山体滑坡与雪崩沉积区。另外统计滑坡中心点位置、高程、坡度、坡向和对应的滑坡面积等信息,以此来分析震害滑坡与地形因子之间的关联性。

    • 吉隆是中尼边境的口岸城镇,位于尼泊尔强震中心的东北方向,距离约66 km。在本次强震与随后的余震中,该镇发生了严重的山体滑坡灾害。图 2中黑色矩形方框表示滑坡区域,椭圆圈出部分表示植被覆盖、人工建筑、阴影等因素导致的失相干噪声区域。图 2(c)中,黑色表示失相关区域。图 2(d)~2(g)分别对应图 2(c)相应绿色矩形方框的震前震后多波段(SWIR2/SWIR1/NIR)合成影像对比图,滑坡体在震后影像图中呈现亮色区域(白线圈出)。

      图  2  吉隆形变场和滑坡体提取图

      Figure 2.  Extraction of Deformation Fields and Landslides of Jilong

      通过分析低SNR与随机性离群值区域快速确定形变失相干,即山体滑坡区域,并通过高分影像目视解译的方法对以上结果进行验证。Landsat 8影像的空间分辨率15 m,实验获取的形变图空间分辨率为60 m×60 m,难以实现对小滑坡区域的监测,导致SPC技术无法监测出微小滑坡体。值得注意的是,震前震后影像的太阳高度角、方位角偏差也会造成趋势性离群值产生。此外,人工建筑设施的拆建、植被变化、云层覆盖等会造成时相失相关。这些误差的识别与去除也会一定程度影响本文方法提取滑坡体的精度。

      利用目视解译方法在吉隆镇识别了21处滑坡体(见图 3)。其中J01为该区域的最大滑坡体,位于冲色村附近河谷,该滑坡堵塞河道造成上游堰塞湖的生成,专家分析结果[14]对该滑坡体进行了确认。失相关区域能够较好地涵盖该区域的大滑坡体(大于3 600 m2),涵盖率为85.7%。该技术对大滑坡体的监测涵盖率未达到100%,主要是因为一部分滑坡体呈线状分布,直径幅宽低于本文设计的测量精度60 m,无法识别。除此之外,本文设定的窗口可能过大,漏检了部分滑坡区域。表 2列出了该区域21个滑坡体的参数与地形信息。

      图  3  吉隆滑坡分布图

      Figure 3.  Landslides Distribution of Jilong

      表 2  吉隆山体滑坡编目表

      Table 2.  Landslides Catalog of Jilong

      序号 滑坡位置 滑坡体面积/m2 技术监测手段 高程/m 坡向/(°) 坡度/(°)
      东经/(°) 北纬/(°) SPC 目视判读
      J01 85.363 28.360 163 646 2 615 289.6 52.8
      J02 85.359 28.354 29 052 2 370 268.7 43.5
      J03 85.337 28.376 37 580 2 490 80.4 30.2
      J04 85.324 28.389 49 045 2 737 312.1 45.9
      J05 85.305 28.401 5 599 2 837 359.1 37.3
      J06 85.304 28.406 33 891 2 809 28.9 40.3
      J07 85.350 28.364 2 603 2 543 10.6 47.5
      J08 85.345 28.369 11 552 2 415 165.9 38.9
      J09 85.337 28.374 5 967 2 450 102.1 33.3
      J10 85.340 28.381 7 431 线状滑坡× 2 839 208.5 53.9
      J11 85.337 28.374 1 295 小滑坡× 2 603 24.7 27.1
      J12 85.355 28.364 680 小滑坡× 2 472 120.0 45.7
      J13 85.370 28.376 6 789 2 925 212.4 44.8
      J14 85.358 28.370 42 480 2 439 196.8 52.6
      J15 85.355 28.372 1 057 小滑坡× 2 550 172.0 41.4
      J16 85.308 28.404 2 263 小滑坡× 2 804 191.8 33.3
      J17 85.303 28.408 1 728 小滑坡× 2 805 49.1 42.9
      J18 85.312 28.399 3 628 小滑坡× 2 741 172.1 25.6
      J19 85.358 28.411 351 小滑坡× 2 984 135.1 37.0
      J20 85.364 28.382 21 275 3 152 256.7 48.2
      J21 85.365 28.382 17 310 3 152 256.7 42.7
      注:√表示可以监测到对应滑坡体,×表示无法监测到对应滑坡体。

      聂拉木区域大部分滑坡体在影像中被云层覆盖,形变失相关区域不够明显,为增加参与后面滑坡因子分析的数据量,本文只采用高分辨率影像目视解译方法获取该区域的10处山体滑坡(如图 4表 3所示)。

      图  4  聂拉木滑坡分布图

      Figure 4.  Landslides Distribution of Nyalam

      表 3  聂拉木山体滑坡编目表

      Table 3.  Landslides Catalog of Nyalam

      序号 滑坡位置 滑坡体面积/m2 高程/m 坡向/(°) 坡度/(°)
      东经/(°) 北纬/(°)
      N01 85.997 28.090 46 565 3 457 123.3 37.6
      N02 85.998 28.087 7 153 3 264 177.7 31.7
      N03 85.990 28.084 13 099 3 769 166.3 60.8
      N04 85.997 28.085 6 360 3 238 270 26.9
      N05 85.991 28.046 38 639 3 265 274.9 57.1
      N06 85.994 28.039 53 518 3 049 171.8 32.8
      N07 85.981 28.032 51 028 3 157 100.9 68.8
      N08 85.980 28.021 12 191 2 602 106.2 45.4
      N09 85.963 27.983 309 600 2 375 157.2 60.8
      N10 86.039 27.975 381 165 3 462 106.5 42.2
    • 本文对吉隆与聂拉木滑坡体数目、滑坡面积、滑坡平均面积(共31处)与地形坡度的关系做简单的统计分析,结果如图 5所示。即每5°作为一个区间范围,通过计算范围内的滑坡个数与滑坡总个数的商,计算其所占百分比统计结果,同时统计区间范围内的滑坡总面积与平均面积。发现滑坡体分布呈现近似高斯正态分布,且主要分布在坡度大于30°(90%以上)的地形,坡度40°~ 45°为中型滑坡体(面积大于半个窗口)多发的坡度范围,且总体滑坡面积最大。同时发现单体(平均)滑坡面积随坡度的增大呈上升趋势,这说明坡度越高越容易触发单体面积大的滑坡,这类震后滑坡危害更为严重,需要重点监测。此外,滑坡体多沿河谷两侧分布,与坡向无明显关系。汶川地震等案例中的震害滑坡分布也得出了类似的统计特性[3, 6]

      图  5  吉隆、聂拉木地区不同坡度水平滑坡体数目、总面积与平均面积分布直方图

      Figure 5.  Landslides' Quantity, Gross and Average Area Distribution Histogram of Jilong and Nyalam in Different Slope Levels

    • 尼泊尔郎塘在此次地震中受灾严重[14-15], 该区域海拔高度在3 000 m以上,地势起伏很大,周边冰川广泛发育。与滑坡体不同,雪崩实际发育区域多为辐射亮度较高的冰雪区域,纹理特征变化不明显,难以识别。而雪崩碎屑沉积区多为海拔较低区域,当地的地面纹理特征一般会遭到严重破坏,因此用本文§2.1中滑坡体识别方法来识别雪崩沉积区域。

      图 6(a)~6(c)中黑线圈出的ABCD为排除云噪声后识别出的失相关区域。图 6(e)中红色虚线圈出部分为云噪声,白线圈出部分为目视解译出的雪崩碎屑物质覆盖区域, 分别对应6(a)、6(b)、6(c)中的A (Chyamki)、B (Langtang Village)、C (Mundu)、D (Singdum)所圈出的失相关区域,本文所获得的结论与Kargel等[14]通过0.5 m分辨率WorldView2目视解译的结果基本吻合。

      图  6  郎塘形变场和雪崩沉积区提取图

      Figure 6.  Extraction of Deformation Fields and Avalanches' Debris Areas in Langtang

    • 随着全球气候变暖,冰川消退速度加快,解冻岩层与冰川融水等因素本身就削弱了局部地质的稳定性。这时如果再发生此类大型地震,引发雪崩滑坡灾害的机率与风险会更大。本文对珠峰附近20 km×20 km区域进行形变监测(表 1中的影像对3),其中昆布冰川5 300 m海拔高度冰瀑区发生大型雪崩,本文重点对其进行形变监测,同时分析了东绒布冰川的疑似雪崩区。如图 7所示,对应冰川区有较为明显的形变量。冰川运动一般由自身重力引起,但是地震导致的大范围雪崩会对冰川局部运动产生影响。本文分别对昆布冰川与东绒布冰川选取了EF两个区域,提取并分析了EF两个典型冰川的水平形变场。

      图  7  珠峰研究区域及形变场

      Figure 7.  Everest Mount Region and Its Deformation Fields

      图 8(a)中红色箭头表示E区域的冰川运动方向与大小,图 8(b)分别表示东西向与南北向形变,黑线圈定区域为昆布冰川的雪崩高发区。对应区域冰川运动了约22~77 m(平均运动39 m,时间基线64 d)。而Scherler等[16]研究表明昆布冰川相应区域的水平运动速率约为20~130 m/a,换算到对应的时间段约为3~21 m,很有可能是由地震触发的雪崩造成的。

      图  8  昆布冰川区影像与形变场

      Figure 8.  Image and Deformation Fields of the Khumbu Glacier

      图 9(a)对应图 7中的F实验区,图 9(b)显示图 9(a)中红色虚线方框区域的震前(2014-03-19-2014-05-22,64 d)、震后(2015-03-22-2015-05-25,64 d)近似时间间隔的东西向与南北向形变场。鉴于两个影像对成像时间基本一致,认为其误差水平相同,而相同区域在地震发生前表现为5~20 m的形变量,地震发生后表现为失相关,说明该区域地表发生了剧烈变化,本文将该区域圈定为疑似雪崩区域。

      图  9  东绒布冰川区影像与形变场

      Figure 9.  Image and Deformation Fields of the East Rongbuk Glacier

    • 震后滑坡灾害的识别与提取对地震触发滑坡机理的探讨和防灾减灾具有重要的理论意义和应用价值。本文从Landsat 8影像所获得形变场中的失相干噪声中提取震害滑坡体与雪崩沉积区。实验结果表明本文方法可行,结果可靠。同时滑坡体与地形因子的统计分析结果表明, 实验区内震后山体滑坡分布与坡度呈近似高斯分布,多发生于坡度大于30°、海拔为2 500~3 500 m的区域,集中沿河谷两侧分布,且高坡度更容易触发单体面积大的山体滑坡,统计结果与其他研究成果[14]较为吻合。另外,本文对昆布冰川与东绒布冰川雪崩区形变进行定量化分析,为该技术应用于雪崩监测预警提供参考。因遥感数据源有限及成像条件不佳,本文只能对某些离散区域进行滑坡、冰川监测,同时影像分辨率与匹配窗口的设置也限制了本文方法对一些小型滑坡体的提取。如果条件允许,SPC技术可以获取整个震区的同震形变与震害滑坡分布,将有助于分析震后滑坡分布与地震机制的关系。此外,综合多源遥感影像提取震害滑坡、雪崩沉积区也值得做进一步研究。

参考文献 (16)

目录

    /

    返回文章
    返回