Processing math: 100%

基于SBAS-InSAR技术的金沙江流域沃达村巨型老滑坡形变分析

冯文凯, 顿佳伟, 易小宇, 张国强

冯文凯, 顿佳伟, 易小宇, 等. 2020. 基于SBAS-InSAR技术的金沙江流域沃达村巨型老滑坡形变分析[J]. 工程地质学报, 28(2): 384-393. doi: 10.13544/j.cnki.jeg.2019-411.
引用本文: 冯文凯, 顿佳伟, 易小宇, 等. 2020. 基于SBAS-InSAR技术的金沙江流域沃达村巨型老滑坡形变分析[J]. 工程地质学报, 28(2): 384-393. doi: 10.13544/j.cnki.jeg.2019-411.
Feng Wenkai, Dun Jiawei, Yi Xiaoyu, et al. 2020. Deformation analysis of Woda village old landslide in Jinsha river basin using SBAS-InSAR technology[J]. Journal of Engineering Geology, 28(2): 384-393. doi: 10.13544/j.cnki.jeg.2019-411.
Citation: Feng Wenkai, Dun Jiawei, Yi Xiaoyu, et al. 2020. Deformation analysis of Woda village old landslide in Jinsha river basin using SBAS-InSAR technology[J]. Journal of Engineering Geology, 28(2): 384-393. doi: 10.13544/j.cnki.jeg.2019-411.

基于SBAS-InSAR技术的金沙江流域沃达村巨型老滑坡形变分析

基金项目: 

国家自然科学基金 41977252

四川省青年科技创新研究团队专项计划项目 2017TD0018

地质灾害防治与地质环境保护国家重点实验室团队项目 SKLGP2016Z001

详细信息
    作者简介:

    冯文凯(1974-),男,博士,教授,博士生导师,主要从事区域稳定及与岩体稳定以及地质灾害评价与防治方面的教学与研究工作.E-mail:fengwenkai@cdut.cn

  • 中图分类号: P642.A

DEFORMATION ANALYSIS OF WODA VILLAGE OLD LANDSLIDE IN JINSHA RIVER BASIN USING SBAS-INSAR TECHNOLOGY

Funds: 

This study is supported by the National Natural Science Foundation of China 41977252

the Sichuan Provincial Youth Science and Technology Innovation Team Special Projects of China 2017TD0018

the Team Project of Independent Research of SKLGP SKLGP2016Z001

  • 摘要: 近年来突发性高位滑塌灾害日益频发,造成恶劣影响。这类地质灾害调查难度高、隐蔽性强,单靠群测群防和地质调查难以解决灾害的防治问题。随着雷达遥感卫星数据质量的不断提升,合成孔径干涉雷达测量(InSAR)中的SBAS-InSAR技术为特大型老滑坡灾前形变探测提供了新的技术途径。利用SBAS-InSAR技术对金沙江流域沃达村滑坡进行地表形变监测,获取了2017年3月30日至2019年9月28日内的形变结果,划定了强烈形变区(Ⅰ雷达)、均匀形变区(Ⅱ雷达),分析了滑坡复活区整体和局部滑塌地表形变速率、累积位移变化趋势和主裂缝形变情况。同时实地进行了工程地质调查和复核,发现老滑坡复活区变形迹象与SBAS-InSAR技术解译成果有着较好的一致性。表明SBAS-InSAR技术在复杂山区地质灾害监测预警领域有较为广阔的应用前景,为类似老滑坡监测预警提供了新的思路与借鉴。
    Abstract: In recent years, unexpected high-level landslides have been increasing, resulting in adverse effects. This kind of geological hazard investigation is very difficult and concealed. It is difficult to solve the problem of disaster prevention and control by group survey and geological investigation alone. With the continuous improvement of the quality of radar remote sensing satellite data, the SBAS-InSAR technology in synthetic aperture Interferometric Radar(InSAR)provides a new technical approach for deformation detection of large-scale old landslides before disaster. The surface deformation of Wodacun landslide in Jinsha River Basin is monitored by SBAS-InSAR technology. The deformation results from March 30, 2017 to September 28, 2019 are obtained. Strong deformation area(Ⅰradar) and uniform deformation area(Ⅱradar) are delineated. The surface deformation rate, cumulative displacement trend and main crack deformation of the whole and local landslide in the revived area are analyzed. At the same time, the engineering geological survey and review are carried out on the spot. It is found that the resurrection and deformation signs of the old landslide are in good agreement with the interpretation results of the SBAS-InSAR technology. It is clear that the SBAS-InSAR technology has broad application prospects in the field of early warning and monitoring of geological disasters in complex mountain areas, and provides new ideas and references for monitoring and early warning similar to old landslides.
  • 随着人类活动深度和广度的日益加深和拓展,我国西部地区大型突发地质灾害日益频发。这些灾害表现出隐蔽性强、危害性高,具有显著的突发性,严重影响了社会经济发展和人民安居乐业,给地质灾害防治与地质环境保护带来了巨大挑战。

    2017年6月24日清晨6时许,四川省茂县叠溪镇新磨村新村组后山约450×104m3山体突发顺层高位滑动,造成80余人失踪或遇难,引起国内外广泛关注(许强等,2017)。灾害发生后,意大利Tre Altamira (2017)通过InSAR(Interferometric Synthetic Aperture Radar,合成孔径雷达干涉)技术对2014年以来45景Sentinel-1雷达卫星影像进行了分析,证实新磨村滑坡的显著形变主要发生在滑前的数月之内。2018年10月11日凌晨4时许,金沙江流域西藏昌都市江达县波罗乡白格村发生特大库岸岩质滑坡,总方量约3165×104m3,滑坡体迅速阻断金沙江,回水造成多条道路中断,川藏滇等地紧急转移安置2.1万余人(冯文凯等,2019)。许强等(2018)通过分析白格滑坡灾前雷达卫星数据,指出该滑坡灾害发生前一年,滑源区斜坡最大位移量约25m。这些案例充分说明InSAR技术对处于高位、交通困难、人迹罕至的隐蔽性地质灾害的识别和监测研究具有显著优势。

    InSAR技术是近20a来兴起的一种新型对地观测技术,具有全天候工作、覆盖范围广、探测精度高、非接触和综合成本低等优点,已被广泛地应用于滑坡(李凌婧等,2017)、地面沉降(张静等,2018)、地裂缝(赵超英等,2011)等地质灾害的形变监测。InSAR技术中的小基线(Small Baseline Subsets,SBAS)技术具有获取微小形变信息和长时间序列缓慢地表变形的优势(莫玉娟等,2018),在复杂山区地质灾害预警和监测领域有较为广阔的应用前景。在基于SBAS-InSAR技术对滑坡早期识别与形变监测的研究领域中,很多学者取得了相应的实质性成果。余睿(2014)采用SBAS-InSAR技术对甘肃省西和县慢速型滑坡进行全局和局部的监测,揭示了影响西和县滑坡体移动的主导因素;刘筱怡等(2017)以鲜河水断裂带内蠕滑型滑坡为研究对象,采用SBAS-InSAR技术获取蠕滑型滑坡时间序列形变特征,分析其发展变化规律,对鲜水河断裂带以及类似构造活动地区的地质灾害防治起到了重要作用;聂兵其等(2018)基于SBAS-InSAR技术对丹巴县城进行滑坡形变探测及隐患识别,验证了InSAR监测结果与现场调查的一致性。可以说,SBAS-InSAR技术为地质灾害隐患点的普查和分析研究提供了重要支撑。在充分利用星载雷达对地观测技术的前提下,开展针对性的工程地质调查和复核工作,一方面有利于正确认识InSAR解译数据的规律性,另一方面也有利于灾害预警研判的准确性,从而大大提升我国科学防灾减灾的技术实力。

    沃达村巨型老滑坡是区域地质灾害调查隐患数据库的隐患点,位于江达县岩比乡沃达村、金沙江右岸,距离白格滑坡直线距离42km。本文在搜集沃达村老滑坡区域升轨哨兵一号(Sentinel-1A)2017年3月30日至2019年9月28日的39景数据的基础上,采用SBAS-InSAR技术分析了老滑坡堆积体的形变特征,并进一步针对性地开展了现场工程地质调查和复核工作,在此基础上分析了SBAS-InSRAR技术在复杂山区地质灾害预警和监测领域的适用性,为类似老滑坡监测预警提供了新的思路与借鉴。

    沃达村老滑坡位于西藏自治区昌都市江达县岩比乡沃达村,坐标东经E93°26′06″、北纬N31°26′06″(图 1)。前期区域地质灾害调查表明,该地段地质历史时期曾发生滑坡灾害,其滑坡堆积体在降水、地震和水位波动情况下易发生复活,应对其重点监测研究。

    图  1  滑坡地理位置图
    Figure  1.  Geographic location map of the landslide

    从综合光学遥感解译图(图 2)分析,老滑坡体纵长约1.9km,横宽约1.7km,主滑方向为NE30°;老滑坡堆积体局部复活,复活部分纵向长约0.88km,横宽约1.08km,面积约0.925km2。滑坡体平面形态为半圆状,前缘呈由缓变陡的陡崖形态,缓部与江面高差最大约240m,前端临空面陡倾,坡脚即为金沙江河谷,高程约2950m,中后部相对平缓。据前期调查显示,该老滑坡上部覆盖层为第四系滑坡堆积物,含碎石、砾石黏性土,厚约20~30m,下伏基岩为变质页岩,岩层产状约为240°∠80°。

    图  2  综合光学遥感解译图
    Figure  2.  Integrated optical remote sensing interpretation map

    本文使用的单视复数SAR影像数据来自欧空局(European Space Agency)2014年4月3日发射的对地观测卫星Sentinel 1A,数据覆盖区域如图 3所示。波段为C波段,所用影像为合成孔径雷达99轨道数据,时间跨度从2017年3月30日至2019年9月28日,24d一景,共39景,轨道方向为升轨,成像中心入射角为36.2°,成像模式为IW宽幅模式,地面分辨率为5m×20m,幅宽为250km×250km,极化方式为VV,方位向和距离向采样间距分别为2.329m、13.956m(表 1)。本研究采用了日本JAXA提供的30m分辨率ALOS World 3D(AW3D30)DEM作为外部数据来消除地形相位的影响。

    图  3  数据覆盖范围图
    Figure  3.  Data coverage chart
    表  1  沃达村滑坡Sentinel 1A数据列表
    Table  1.  Data list of sentinel 1 in Dawo Village landslide
    序号 成像时间 成像
    模式
    极化
    方式
    累积时
    间基线
    /d
    累积空
    /m
    1 20170330 IW VV 0 0
    2 20170423 IW VV 24 -46.119
    3 20170517 IW VV 48 -50.186
    4 20170610 IW VV 72 77.826
    5 20170704 IW VV 96 40.558
    6 20170809 IW VV 132 -10.291
    7 20170902 IW VV 156 -39.573
    8 20170926 IW VV 180 -50.931
    9 20171020 IW VV 204 114.923
    10 20171113 IW VV 228 -24.411
    11 20171207 IW VV 252 50.346
    12 20171231 IW VV 276 63.293
    13 20180124 IW VV 300 -165.698
    14 20180217 IW VV 324 27.433
    15 20180313 IW VV 348 30.264
    16 20180406 IW VV 372 -6.096
    17 20180430 IW VV 396 23.623
    18 20180524 IW VV 420 -147.052
    19 20180617 IW VV 444 116.508
    20 20180711 IW VV 468 43.829
    21 20180723 IW VV 492 -86.624
    22 20180816 IW VV 516 -15.263
    23 20180909 IW VV 540 64.119
    24 20181003 IW VV 564 -43.180
    25 20181027 IW VV 588 35.953
    26 20181120 IW VV 612 26.433
    27 20181214 IW VV 636 -18.007
    28 20190107 IW VV 660 31.303
    29 20190131 IW VV 684 -57.897
    30 20190224 IW VV 708 -62.838
    31 20190320 IW VV 732 61.137
    32 20190413 IW VV 756 -40.041
    33 20190507 IW VV 780 81.387
    34 20190531 IW VV 804 -93.348
    35 20190624 IW VV 828 27.419
    36 20190718 IW VV 852 18.360
    37 20190811 IW VV 876 -105.069
    38 20190904 IW VV 900 36.953
    39 20190928 IW VV 924 41.097
    下载: 导出CSV 
    | 显示表格

    为减少时间和空间去相干,Berardino et al.(2002)Casu(2000)研究的基础上,提出了短基线差分干涉法。该方法的基本原理是利用具有较短时-空基线的影像对产生干涉图提高相干性,再进行解缠干涉相位。依据相干像元相位和观测时间的关系,采用奇异值分解方法将多个小基线集合数据联合起来求解(许军强等,2019),有效地解决了各数据集之间空间基线过长造成的时间不连续问题,提高了监测的时间分辨率,从而求得影像序列间地表形变速率的最小范数最小二乘解(何秀凤等,2012)。

    其基本原理:假设有N+1景单视复数影像,其成像时间为(t0t1,….,tn),设某一差分解缠图为j,某像元(xr)对应的解缠相位可表示为式(1):

    δφj( x,r) =φ (tB,x,r) φ (tA,x,r)δφtopoj (x,r) +δφdispj (x,r) +δφatmj(x,r) +δφnoisej (x,r) (1)

    式中,δφjtopo(xr)为DEM误差所形成的相位;δφjdisp(xr)为坡体形变产生的相位;δφjatm(xr)为在tBtA时刻间大气影响产生的相位;δφjnoise(xr)为噪声导致的相位,式(1)右端前3项可表示为式(2):

    {δφtopoj (x,r)=4 π   BjΔZλ Rsinθδφdispj( x,r)=4 π λ  [d (tB,x,r)d( tA,x,r)]j=1,Mδφatmj(x,r)=φatm (tB,x,r)φatm (tA,x,r)  (2)

    式中λ为波长;θ为雷达视角;Δz为DEM误差;R为雷达距观测物的斜距;d(tBxr)和d(tAxr)分别表示相对于参考时间t0的视向累积变形量。M为干涉图数量。假定N代表全部SAR影像数,则有M2MN(N1)2,对M个解缠相位进行最小二乘或奇异值分解,即:

    由于在处理过程中产生M幅干涉图,根据式(2)可得M个方程,以矩阵形式表示为式(3):

    δφ(x,r)=Aφ(x,r) (3)

    其中,AM×N系数矩阵;φ(xr)为(xr)点在N个时刻对应的未知形变相位构成的矩阵,当MN时,通过最小二乘法得出式(4):

    φ(x,r)=(ATA)1ATδφ(x,r) (4)

    MN时,方程有无数解,使用SVD方法联合求解多个小基线,最终可得到不同时刻所对应的累积形变量。

    本文使用ENVI中SArscape模块,用SBAS算法对Sentinel 1A数据进行干涉处理,共生成109个干涉对,时间基线为75d,临界基线百分比为2%,各像对时间基线和空间基线的连接方式如图 4图 5所示,图中绿色的点代表 39景影像的成像时间,黄色点代表作为参考主影像的成像时间。对合成的干涉图做第1次反演、第2次反演、地理编码以及栅矢转换,生成地表沿雷达视向的平均形变速率(mm ·a-1),最终得到时间序列上的形变结果。具体处理流程如图 6所示。

    图  4  影像时间基线
    Figure  4.  Image time baseline
    图  5  影像空间基线
    Figure  5.  Image spatial baseline
    图  6  SBAS处理流程
    Figure  6.  SBAS-InSAR processing flow

    由于滑坡多沿斜坡面进行滑动,雷达视线方向的形变信息无法准确反映斜坡面的真实形变情况(戴可人等,2019),考虑雷达视线方向、坡度方向和垂直沉降方向的几何关系(图 7),假设运动沿单位矢量û指定的方向发生,采用式(5)、式(6)(Cascini et al., 2010)和式(8)(Colesanti et al., 2006)将视线方向的形变速率转化为坡度方向和垂直方向的形变速率。

    图  7  雷达视线方向和坡度方向几何示意图(Cascini et al., 2010)
    Figure  7.  Geometric sketches of radar line- of -sight direction and slope direction(Cascini et al., 2010)
    ˆu=[sinαcosφ cosαcosφsinφ]  (5)
    VSlope=VLos/cosβ (6)

    式(6)中,cosβ可用式(7)(Cascini et al., 2010)计算;

    cosβ=[(sinαcosφ)(sinθcosαs)+(cosαcosφ)(sinθsinαs)+ sinφcosθ] (7)
    Vu=VLos+VSlopesinθcos [δαs32 π ]cosθ (8)

    其中,VSlope代表沿坡度方向的形变速率;VLos代表沿雷达视线方向的形变速率;Vu代表垂直方向的形变速率;αs为方位向和正北方向的夹角; αs32 π 为方位向视线方向;δ为滑坡方位角;α为斜坡坡向;β为视坡夹角;θ为入射角;φ为斜坡坡度。

    避免在VLosVSlope转换过程中存在绝对值的极大异常现象,采用Herrera et al.(2013)提出cosβ=±0.3为固定阈值,即VSlope不能大于VLos的3.33倍,设- 0.3<cosβ<0时,cosβ=-0.3;0<cosβ<0.3时,cosβ=0.3。

    结合遥感影像分析,该滑坡中部和后缘部位表面植被覆盖密度不大,相干点较多,干涉效果较好,有较明显的速率积累,滑坡部分地方植被覆盖茂密,造成影像的失相干,速率点较为稀少。

    形变速率为负值表示形变点向远离卫星传感器方向运动,正值表示形变点向靠近卫星传感器方向运动。由InSAR视向形变分析结果可见,图 8中形变值均为负值,表明滑坡整体向着金沙江方向向下滑动,符合滑坡体的运动规律。从图 8中可以看出该滑坡形变较强部位位于滑坡中部和中前缘(图 8紫色虚线圈内),其中紫色虚线圈内的红色区域形变量最大,形变速率≥80mm ·a-1

    图  8  2017年3月~2019年9月形变速率图(LOS方向)
    Figure  8.  Deformation rate map from March 2017 to September 2019(LOS Direction)

    沃达村滑坡沿坡向年平均速率(VSlope)和垂直方向年平均速率(Vu)分别如图 9图 10所示,解译出的强变形区与LOS方向大致一致,但在形变速率量值上存在较大差异,沿坡向形变速率最大。

    图  9  2017年3月~2019年9月形变速率图(Slope方向)
    Figure  9.  Deformation rate map from March 2017 to September 2019(Slope Direction)
    图  10  2017年3月~2019年9月形变速率图(Vertical方向)
    Figure  10.  Deformation rate map from March 2017 to September 2019(Vertical Direction)

    为复核解译数据,开展了实地工程地质调查与复核工作。经实地踏勘,滑坡体前缘地形陡峭;中前部坡度较大,树木生长茂盛;中部为缓坡平台,主要为耕地,树木零星生长;后部为古滑坡滑源区(图 11a)。受冲沟下切侵蚀及前缘陡倾临空面综合作用,滑坡前缘发育多处向冲沟、临空面方向的崩滑变形体(B01~B04),崩滑变形体表部裂缝横张、错坎发育,裂缝张开2~10cm,最大错距1~2m,前端陡崖可见明显基覆界面(图 11b)。滑坡复活区中部发育多处小型滑动变形体(H01~H03),坡体发育横切梯级陡坎,高度几十厘米到几米不等,延伸几十米到数百米不等,陡坎上部多耕作,有灌木丛生,表明陡坎形成时间较长(图 11c);可见明显马刀树,结合树龄推测稳定时间至少10a以上(图 11d)。

    图  11  沃达村老滑坡照片
    Figure  11.  Photographs of the old landslide in Woda Village

    复活区主要位于老滑坡体中下部,强烈变形区(Ⅰ现场)后缘位于老滑坡后缘陡缓交界处,两侧边界依地形延伸至老滑坡边界(图 8),其内裂缝、下错台坎、局部滑塌现象明显。现场的变形迹象与SBAS-InSAR方法获取的形变速率图(图 8)特征有着较高的一致性。在滑坡视向形变速率图(图 8)形变量较大的红色区域取一点A作出该点历史累积位移形变曲线,与搜集的2017年1月~2019年8月月累计降雨曲线对比分析(图 12):随着2017-05~2017-06、2018-05~2018-08、2019-06~2019-07 3个不同时间段内降雨量急剧增加,A点形变量明显加速;降雨量较小的时间段内,A点形变量增加变缓,形变趋势与降雨具有较好的响应规律。据现场调查,滑坡所处金沙江水位波动最高高程为2930~2950m,而滑坡局部复活区位置位于3220~3580m,远离金沙江水位波动影响区。说明该滑坡中前部复活区主要受降雨影响,为老滑坡降雨复活型滑坡。

    图  12  InSAR时序形变与降雨响应
    Figure  12.  InSAR time series deformation and rainfall response

    从滑坡视向形变速率(图 8)分析,坡体中部(紫色虚线圈闭范围内)灰色(-40~-20mm ·a-1)、黄色(-60~-40mm ·a-1)、橙色(-80~-60mm ·a-1)、红色(-102~-80mm ·a-1)等颜色相互交错,分布杂乱,差异性很大,说明该区域变形强烈是滑坡体裂缝台坎发育的体现。滑坡区前缘、后缘雷达形变速率颜色为单一蓝色(-20~0mm ·a-1),故在雷达监测时间段,该区域形变较小,差异变形不明显。通过以上雷达形变速率及分区分析,将形变速率整体分布趋势(图 8)分成强烈形变区(Ⅰ雷达)和均匀形变区(Ⅱ雷达),分析2017年3月~2019年9月不同时段内沿滑坡主滑方向的视向形变速率(图 13a~ 图 15a),发现每年老滑坡前后缘形变量相对较小,中部的形变量较大,2017~2019年分别最大可达约- 125mm ·a-1、-100mm ·a-1、-55mm ·a-1。由于坡体前缘较陡,前缘坡度与坡体整体坡度相差太大,考虑地形因素引起的叠掩和阴影造成的前缘SAR影像失真(张毅,2018),老滑坡前缘形变速率结果存在不确定性。综合现场调查和形变速率分析结果,认为滑坡现阶段滑动变形以向冲沟、陡倾临空面方向为主,破坏形式以前缘散落、局部滑塌为主,整体表现老滑坡局部复活的牵引式滑动,而老滑坡中前部主要堆积体覆盖区则表现为中部向前部的强烈推挤作用。

    图  13  2017年3~12月形变速率及剖面形变图
    Figure  13.  Deformation rate and profile deformation map from March to December 2017
    图  14  2018年形变速率及剖面形变图
    Figure  14.  Deformation rate and profile deformation map of 2018
    图  15  2019年1~9月形变速率及剖面形变图
    Figure  15.  Deformation rate and profile deformation map from January to September 2019

    现场复核阶段对前缘临空面发育的4处崩滑变形体(B01~B04)和滑坡中、下部3处滑动变形体(H01~H03)进行了细致调查。4处崩滑变形体均为表层发生的溜滑,规模较小,树木歪斜的范围位于冲沟沟口,表明为集中降雨引起的变形。3处滑动变形体中,2处(H02和H03)因为前部临空面变形导致出现牵引式滑动变形,1处(H01)因降雨和人类工程活动因素发生牵引式滑动变形。下面主要分析滑动变形体H01的情况。

    H01滑动变形体位于坡体的中后部,沃达村公路上侧,滑动方向与坡向一致(图 16a)。据访,H01滑动变形体前缘为沃达村村道,村道挡土墙于2013年左右发生隆起破坏(图 16b);2018年中旬,H01滑动变形体发生了较为强烈的变形,导致挡墙上方较高处田地出现多处裂缝及下错台坎(图 16c)。调查发现,其下部紧邻的村民房屋完好无破坏,表明其仅仅是近年来表层岩土体在降雨等作用下发生的局部牵引式滑动变形。

    图  16  H01变形体SBAS-InSAR监测与野外现象对比
    Figure  16.  Comparison of SBAS-InSAR monitoring data of H01 landslide with field review

    为细致分析H01滑动变形体的滑动性质与雷达形变速率信号之间的相关性,沿H01滑动变形体主滑方向布置剖面线,并依次取1#~4#监测点,利用图 9绘制其3年内滑坡历史累积位移形变曲线(图 16d)。如图所示,4#监测点历史累积形变量较大,增加的速率也较大;1#~3#监测点历史累积形变量增加幅度基本一致。前缘的形变速率远大于中部和后部的形变速率,说明H01滑动变形体属于牵引式滑动变形。

    根据对滑塌体后缘下错台坎发育处沿坡向以及垂向形变速率点的分析:沿坡向的形变速率(图 9图 16a)颜色各异,形变速率差异性很大,形变速率140mm ·a-1、160mm ·a-1、180mm ·a-1均有分布,与裂缝形成特征吻合;沿垂向形变速率分析(图 10),在后缘位置黄色形变速率(-80~-60mm ·a-1)附近出现多级的橙色形变速率(-100~-80mm ·a-1),说明该区域发生多级下错沉降现象。两种现象印证了滑坡强变形区历史出现多级下错台坎,表现出蠕动变形,与地质调查的现象基本相符。

    (1) 通过对视向、坡向、垂向方向的形变速率分区以及形变点之间的差异性分析,初步将沃达村滑坡分为强烈形变区(Ⅰ雷达)和均匀形变区(Ⅱ雷达),与现场判定的强烈变形区(Ⅰ现场)和弱变形区(Ⅱ现场)基本吻合。

    (2) 根据2017年3月~2019年9月3a不同时段内沿滑坡主滑方向的剖面形变速率分析,后缘、前缘形变量较小,中部的形变量较大,分别最大可达约- 125mm ·a-1、-100mm ·a-1、-55mm ·a-1(LOS)。考虑地形因素对InSAR前缘形变速率的影响并结合现场调查分析认为滑坡现阶段变形以向冲沟、陡倾临空面方向为主,破坏形式以前缘散落、局部滑塌为主,整体表现为老滑坡局部复活的牵引式滑动,而老滑坡中前部主要堆积体覆盖区则表现为中部向前缘的强烈推挤作用。

    (3) 基于坡向形变速率图、垂向形变速率图中形变速率值的分布以及量值的差异性特征,初步判定了局部滑塌的性质以及多级台坎发育的位置,与野外地质调查吻合性较好。

    (4) 通过SBAS-InSAR技术分析以及现场调查复核,发现老滑坡复活变形迹象与SBAS-InSAR技术解译成果有着较好的一致性,在复杂山区地质灾害预警和监测领域有较为广阔的应用前景,为类似老滑坡的监测预警提供新的思路与借鉴。

    致谢: 本文在完成过程中得到了中国地震局地壳应力研究所、四川省国土空间生态修复与地质灾害防治研究院、江达县气象局等单位的帮助。欧洲太空局为本研究提供的Sentinel-1A数据,日本JAXA提供的30m分辨率“ALOS World 3D—30m” DEM。在此一并表示感谢。
  • 图  1   滑坡地理位置图

    Figure  1.   Geographic location map of the landslide

    图  2   综合光学遥感解译图

    Figure  2.   Integrated optical remote sensing interpretation map

    图  3   数据覆盖范围图

    Figure  3.   Data coverage chart

    图  4   影像时间基线

    Figure  4.   Image time baseline

    图  5   影像空间基线

    Figure  5.   Image spatial baseline

    图  6   SBAS处理流程

    Figure  6.   SBAS-InSAR processing flow

    图  7   雷达视线方向和坡度方向几何示意图(Cascini et al., 2010)

    Figure  7.   Geometric sketches of radar line- of -sight direction and slope direction(Cascini et al., 2010)

    图  8   2017年3月~2019年9月形变速率图(LOS方向)

    Figure  8.   Deformation rate map from March 2017 to September 2019(LOS Direction)

    图  9   2017年3月~2019年9月形变速率图(Slope方向)

    Figure  9.   Deformation rate map from March 2017 to September 2019(Slope Direction)

    图  10   2017年3月~2019年9月形变速率图(Vertical方向)

    Figure  10.   Deformation rate map from March 2017 to September 2019(Vertical Direction)

    图  11   沃达村老滑坡照片

    Figure  11.   Photographs of the old landslide in Woda Village

    图  12   InSAR时序形变与降雨响应

    Figure  12.   InSAR time series deformation and rainfall response

    图  13   2017年3~12月形变速率及剖面形变图

    Figure  13.   Deformation rate and profile deformation map from March to December 2017

    图  14   2018年形变速率及剖面形变图

    Figure  14.   Deformation rate and profile deformation map of 2018

    图  15   2019年1~9月形变速率及剖面形变图

    Figure  15.   Deformation rate and profile deformation map from January to September 2019

    图  16   H01变形体SBAS-InSAR监测与野外现象对比

    Figure  16.   Comparison of SBAS-InSAR monitoring data of H01 landslide with field review

    表  1   沃达村滑坡Sentinel 1A数据列表

    Table  1   Data list of sentinel 1 in Dawo Village landslide

    序号 成像时间 成像
    模式
    极化
    方式
    累积时
    间基线
    /d
    累积空
    /m
    1 20170330 IW VV 0 0
    2 20170423 IW VV 24 -46.119
    3 20170517 IW VV 48 -50.186
    4 20170610 IW VV 72 77.826
    5 20170704 IW VV 96 40.558
    6 20170809 IW VV 132 -10.291
    7 20170902 IW VV 156 -39.573
    8 20170926 IW VV 180 -50.931
    9 20171020 IW VV 204 114.923
    10 20171113 IW VV 228 -24.411
    11 20171207 IW VV 252 50.346
    12 20171231 IW VV 276 63.293
    13 20180124 IW VV 300 -165.698
    14 20180217 IW VV 324 27.433
    15 20180313 IW VV 348 30.264
    16 20180406 IW VV 372 -6.096
    17 20180430 IW VV 396 23.623
    18 20180524 IW VV 420 -147.052
    19 20180617 IW VV 444 116.508
    20 20180711 IW VV 468 43.829
    21 20180723 IW VV 492 -86.624
    22 20180816 IW VV 516 -15.263
    23 20180909 IW VV 540 64.119
    24 20181003 IW VV 564 -43.180
    25 20181027 IW VV 588 35.953
    26 20181120 IW VV 612 26.433
    27 20181214 IW VV 636 -18.007
    28 20190107 IW VV 660 31.303
    29 20190131 IW VV 684 -57.897
    30 20190224 IW VV 708 -62.838
    31 20190320 IW VV 732 61.137
    32 20190413 IW VV 756 -40.041
    33 20190507 IW VV 780 81.387
    34 20190531 IW VV 804 -93.348
    35 20190624 IW VV 828 27.419
    36 20190718 IW VV 852 18.360
    37 20190811 IW VV 876 -105.069
    38 20190904 IW VV 900 36.953
    39 20190928 IW VV 924 41.097
    下载: 导出CSV
  • 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 & Remote Sensing, 40(11): 2375-2383.

    Cascini L, Fornaro G, Peduto D. 2010. Advanced low-and full-resolution D-InSAR map generation for slow-moving landslide analysis at different scales[J].Engineering Geology, 112(1): 29-42.

    Casu F, Manzo M, Lanari R. 2000. A quantitative assessment of the SBAS algorithm performance for surface deformation retrieval from D-InSAR data[J].Remote Sensing of Environment, 102(3): 195-210.

    Colesanti C, Wasowski J. 2006. Investigating landslides with space-borne Synthetic Aperture Radar(SAR)interferometry[J].Engineering Geology, 88(3): 173-199.

    Dai K R, Zhuo G C, Xu Q, et al. 2019. Tracing the pre-failure two-dimensional surface displacements of Nanyu Landslide, Gansu province with radar interferometry[J].Geomatics and Information Science of Wuhan University, 44(12): 1778-1786.

    Feng W K, Zhang G Q, Bai H L, et al. 2019. A preliminary analysis of the formation mechanism and development tendency of the huge Baige landslide in Jinshan River on October 11, 2018[J]. Journal of Engineering Geology, 27(2): 415-425.

    He X F, He M. 2012. InSAR earth observation data processing method and comprehensive measurement[M].Beijing: Science Press.

    Herrera G, Gutierrez F, Garcia-Davalillo J C, et al. 2013. Multi-sensor advanced DinSAR monitoring of very slow landslides: The Tena Valley case study(Central Spanish Pyrenees)[J]. Remote Sensing of Environment, 128(none): 31-43.

    Kang Y. 2016. Application of InSAR technology in landslide detection and monitoring in southwest mountainous areas[D].Xi'an: Chang'an University.

    Li L J, Yao X, Zhou Z K, et al. 2017. The deformation characteristics of a large landslide before and afire impoundment in the Xiluodu reservoir area based on InSAR technology[J].Journal of Engineering Geology, 25 (S1): 458-462.

    Liu X Y, Yang Z H, Guo C B, et al. 2017. Study of slow-moving landslide characteristics based on the SBAR-InSAR in the Xianshuihe fault zone[J].Geoscience, 31(5): 965-977.

    Mo Y J, Wu Y, Liu X W. 2018. Monitoring the ground subsidence in Xiaojin County, Sichuan province based on small baseline subset technique[J].Engineering of Surveying and Mapping, 27(11): 46-50.

    Nie B Q. 2018. Landslide deformation detection and identification based on InSAR technology——A case of Danba County[D]. Chengdu: Chengdu University of Technology.

    Tre altamira. 2017. Data in focus: precursor of Maoxian landslide measured from space[EB/OL].http://tre-altamira.com/news/data-focus-precursor-maoxian-landslide-measured-space/.2017-06-29.

    Wang J. 2018. Long-term spaceborne InSAR technology landslide geological disaster monitoring research[D].Beijing: Beijing Jiaotong University.

    Xu J Q, Ma T, Lu Y K, et al. 2019. Land subsidence monitoring in North Henan plain based on SBAS-InSAR technology[J].Journal of Jilin University(Earth Science Edition), 49(4): 1182-1191.

    Xu Q, Li W L, Dong X J, et al. 2017. The Xinmocun landslide on June 24, 2017 in Maoxian, Sichuan: Characteristics and failure mechanism[J].Chinese Journal of Rock Mechanics and Engineering, 36(11): 2612-2628.

    Xu Q, Zheng G, Li W L, et al. 2018. Study on successive landslide damming events of Jinsha River in Baige Village on Octorber 11 and November 3, 2018[J]. Journal of Engineering Geology, 26(6): 1534-1551.

    Yu R. 2014. Xi'an County landslide monitoring research based on short baseline(SBAS)technology[D].Nanjing: Nanjing Normal University.

    Zhang J, Feng D X, Qi W, et al. 2018. Monitoring land subsidence in Panjin region with SBAS-InSAR method[J].Journal of Engineering Geology, 26(4): 999-1007.

    Zhang Y. 2018. Surface deformation monitoring based on InSAR technology and early identification of landslide[D].Lanzhou: Lanzhou University.

    Zhao C Y, Zhang Q, Zhang J. 2011. Deformation monitoring of ground fissure with SAR inter-ferometry in Qingxu, Shanxi province[J].Journal of Engineering Geology, 19(1): 70-75.

    戴可人, 卓冠晨, 许强, 等. 2019.雷达干涉测量对甘肃南峪乡滑坡灾前二维形变追溯[J].武汉大学学报(信息科学版):44(12): 1778-1786.
    冯文凯, 张国强, 白慧林, 等. 2019.金沙江"10 ·11"白格特大型滑坡形成机制及发展趋势初步分析[J].工程地质学报, 27(2): 415-425. DOI: 10.13544/j.cnki.jeg.2018-392
    何秀凤, 何敏. 2012. InSAR对地观测数据处理方法与综合测量[M].北京:科学出版社.
    康亚. 2016. InSAR技术在西南山区滑坡探测与监测的应用[D].西安: 长安大学.
    李凌婧, 姚鑫, 周振凯, 等. 2017.溪洛渡库区某大型滑坡蓄水前后变形特征InAR分析[J].工程地质学报, 25 (S1): 458-462. DOI: 10.13544/j.cnki.jeg.2017.s1.071
    刘筱怡, 杨志华, 郭长宝, 等. 2017.基于SBAS-InSAR的鲜水河断裂带蠕滑型滑坡特征研究[J].现代地质, 31(5): 965-977. DOI: 10.3969/j.issn.1000-8527.2017.05.008
    莫玉娟, 吴洋, 刘学武. 2018.基于SBAS技术的四川阿坝州小金县地表形变监测[J].测绘工程, 27(11): 46-50.
    聂兵其. 2018.基于InSAR的滑坡形变探测及隐患识别研究——以丹巴县城区为例[D].成都: 成都理工大学.
    许军强, 马涛, 卢意恺, 等. 2019.基于SBAS-InSAR技术的豫北平原地面沉降监测[J].吉林大学学报(地球科学版), 49(4): 1182-1191.
    许强, 李为乐, 董秀军, 等. 2017.四川茂县叠溪镇新磨村滑坡特征与成因机制初步研究[J].岩石力学与工程学报, 36(11): 2612-2628.
    许强, 郑光, 李为乐, 等. 2018.2018年10月和11月金沙江白格两次滑坡—堰塞堵江事件分析研究[J].工程地质学报, 26(6): 1534-1551. DOI: 10.13544/j.cnki.jeg.2018-406
    余睿. 2014.基于短基线(SBAS)技术的西和县滑坡监测研究[D].南京: 南京师范大学.
    张静, 冯东向, 綦巍, 等. 2018.基于SBAS-InSAR技术的盘锦地区地面沉降监测[J].工程地质学报, 26(4): 999-1007. DOI: 10.13544/j.cnki.jeg.2017-382
    张毅, 2018基于InSAR技术的地表变形监测与滑坡早期识别研究[D].兰州: 兰州大学.
    赵超英, 张勤, 张静. 2011.山西清徐地裂缝形变的InSAR监测分析[J].工程地质学报, 19(1): 70-75. DOI: 10.3969/j.issn.1004-9665.2011.01.011
图(16)  /  表(1)
计量
  • 文章访问数:  3292
  • HTML全文浏览量:  1397
  • PDF下载量:  381
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-07
  • 修回日期:  2019-12-26
  • 网络出版日期:  2021-04-21
  • 刊出日期:  2020-04-24

目录

/

返回文章
返回