基于光学遥感和SBAS-InSAR的川渝输电网滑坡隐患早期识别
Early identification of potential landslides for the Sichuan-Chongqing power grid based on optical remote sensing and SBAS-InSAR
通讯作者: 汤明高(1978-),男,博士,教授,主要从事地质灾害机理、评价预测与防治研究。Email:tomyr2008@163.com。
责任编辑: 李瑜
收稿日期: 2022-08-10 修回日期: 2023-02-9
基金资助: |
|
Received: 2022-08-10 Revised: 2023-02-9
作者简介 About authors
赵华伟(1998-),男,硕士研究生,主要从事岩土工程专业相关研究。Email:
近年来,山区输电网工程遭遇滑坡而导致输电杆塔倒塌、电力中断的事故屡有发生,早期识别滑坡隐患、防患于未然,对于保障电力工程安全具有重要意义。为此,依托川渝输电网工程,采用光学遥感与小基线集(small baseline subset,SBAS)-合成孔径雷达差分干涉测量(interferometric synthetic aperture radar,InSAR)相结合的SBAS-InSAR方法,对川渝输电网沿线区域进行了滑坡隐患点的早期识别。通过对高分辨率光学遥感影像的解译,共识别出电网沿线杆塔附近的滑坡隐患点28处。在此基础上,采用SBAS-InSAR技术针对研究区进行地表形变探测,发现滑坡隐患点27处。上述2种方法共识别出40处滑坡隐患点,其中15处隐患点为2种方法共同识别。最后,通过现场复核、变形迹象及稳定性定性分析,认为7处隐患点对电网杆塔具有威胁而存在风险、其中2处风险较大。该成果对于川渝输电线路的滑坡地质灾害防治具有指导和参考价值。
关键词:
Power grid projects in mountainous regions have encountered numerous landslides in recent years, leading to collapsed transmission towers and power outages. Hence, early identification of potential landslides is crucial for ensuring the safety of power engineering. For this purpose, this study conducted early identification of potential landslides along the Sichuan-Chongqing power grid based on optical remote sensing and the small baseline subset (SBAS) - interferometric synthetic aperture radar (InSAR) technology. The interpretation of high-resolution optical remote sensing images revealed 28 potential landslide sites near the transmission towers along the power grid. Based on this, this study detected the study area’s surface deformation using the SBAS-InSAR technology, identifying 27 potential landslide sites. Except for 15 repeated results, the above two methods identified a total of 40 potential landslide sites. Finally, through field check and the qualitative analysis of deformation signs and stability, this study determined that seven potential landslide sites threaten the safety of transmission towers, with two of them presenting higher risks. These findings provide valuable guidance and references for the prevention and control of landslides along the Sichuan-Chongqing power grid.
Keywords:
本文引用格式
赵华伟, 周林, 谭明伦, 汤明高, 童庆刚, 秦佳俊, 彭宇辉.
ZHAO Huawei, ZHOU Lin, TAN Minglun, TANG Minggao, TONG Qinggang, QIN Jiajun, PENG Yuhui.
0 引言
输电杆塔是高压输电网中非常重要的基础设施之一,输电杆塔能否安全可靠运营关系着地方经济发展[1]。川渝输电网杆塔多位于我国西南山区的山坡后缘或者山脊处,而西南山区地质条件复杂,地质构造发育,尤其是以滑坡为代表的地质灾害频繁发生,对川渝输电网长期正常运营构成很大威胁。因此,开展输电网沿线滑坡隐患调查和早期识别具有重要的意义。
调查滑坡隐患点的传统方式包括实地踏勘和传统测量手段(水准仪、GPS、裂缝计、观测站等),但是传统测量方式只能对小范围区域的滑坡进行调查,存在一定的局限性; 且效率低、成本高、易受外界环境条件的制约,尤其对于高位远程、高隐蔽性的滑坡识别起来更为困难[2]。随着科学技术的不断发展,光学遥感技术和合成孔径雷达干涉测量技术(interferometric synthetic aperture radar,InSAR)等被更广泛地应用于滑坡隐患的识别中,具有识别效率较高、识别范围很广、识别成本较低等优势[3-4]。但两者也存在一定的局限[5-6]。光学遥感技术具有分辨率高、覆盖范围广、影像信息量丰富等特点,通过地形地貌特征分析可以快速识别出发生过整体变形的古(老)滑坡以及具有明显变形特征的滑坡; 但是,对于正在发生缓慢变形而又无明显变形迹象的滑坡隐患则很难识别出来,且西南地区的光学遥感影像质量受云雾遮蔽的影响较大,这对识别准确率也存在一定影响[7]。现如今InSAR技术已经较为成熟,可以利用雷达影像的相位信息来探测mm级的地表形变量,具有测量精度高、不受云层影响、可以全天运作、测量范围广等特点[8⇓-10],能够很好地弥补光学遥感技术的不足,获得的地表形变量对分析滑坡体形变的演化过程、发生机理及机制、早期识别及预警有重大意义[11]。但InSAR技术也存在一定的缺点: ①容易受植被、大气信号或地形效应的影响,导致探测结果不准确; ②监测范围有限,如果隐患点变形过快,则会造成失相干现象; ③由于技术本身的限制,生成的有效干涉点很难覆盖整个研究区域[12]。由于光学遥感技术与InSAR技术都存在一定的优势与局限性,为避免在滑坡隐患的早期识别过程中出现遗漏的现象,许强等[13-14]和葛大庆等[15⇓-17]提出,综合利用具有高分辨率的光学遥感影像与InSAR技术所获得的形变信息,充分发挥2种技术手段的优势,实现对地质灾害隐患的识别。
其中,小基线集(small baseline subset,SBAS)时序InSAR方法能够有效利用冗余干涉对进行干涉相位建模,可有效抑制形变恢复模型中的误差项,已在滑坡运动监测等领域中得到广泛应用[18]。刘沛源等[19]利用87景降轨的Sentinel-1A影像,采用SBAS-InSAR技术实现了对成汶高速公路成汶段滑坡隐患点的早期识别,共识别滑坡隐患点10处,经野外验证均为处于缓慢变形中的滑坡,为高速公路选线提供了依据。张家勇等[20]利用2017—2020年间获取的97景Sentinel-1A影像,基于SBAS-InSAR技术,对纳雍、黔西2县进行了滑坡隐患点的早期识别,经探测发现了60处隐患点,并通过现场调查验证了InSAR探测结果的有效性,为该地区的防灾减灾提供了依据。
川渝输电网位于西南地区,植被覆盖茂密,云雾遮挡严重,仅采用单一技术手段很难对该地区滑坡隐患进行准确、全面的识别。因此,本文参考前人研究的经验,利用多源遥感数据,包括Google Earth影像、高分2号光学遥感影像和Sentinel-1A数据SBAS-InSAR综合分析,开展川渝输电网黄万线路沿线的滑坡隐患识别,旨在为川渝输电网地质灾害的监测与预防提供参考依据。
1 研究区概况
川渝输电网“黄万线”线路位于E106°50'~108°16',N30°32'~30°50'区域,穿过的行政区包括四川省广安市、达州市以及重庆市,路线总体为SW—NE走向,西起广安市前锋区,经达州市大竹县、重庆市梁平区,止于重庆市万州区,全长141 km(图1)。研究区输电网黄万线位于四川盆地向川东平行岭谷过度地带,输电网穿越了5大山脉,其中包括华蓥山山脉、铜锣山山脉、明月山山脉、铁峰山山脉以及挖断山山脉,地势北高南低,西高东低。研究区处于亚热带湿润气候与亚热带季风湿润气候之间,每年降雨量为1 000~1 200 mm。地表植被覆盖面积较大,包括阔叶林、针叶林、竹林、灌丛、草丛和草甸,植被分布受海拔的影响明显。在区域构造上分属川东褶带、川中褶带、川东南陷褶带和大巴山弧形断褶带单元,出露地层主要为中生界侏罗系到三叠系的砂泥岩、灰岩、板岩,少部分上古生界二叠系地层,以及第四系崩坡积和残积层等。研究区属于我国地质灾害易发区,但是对于川渝输电网以前并未开展遥感调查与监测,在电网运行过程中调查发现存在滑坡变形迹象,因此开展本项研究对于保障输电线路的正常运营具有非常重要的意义。
图1
图1
研究区位置图及雷达影像覆盖图
Fig.1
Map of location of study area and Radar image coverage
2 多源遥感识别方法
首先利用Google Earth、高分二号2种光学遥感影像,通过目视解译对研究区滑坡隐患点进行早期识别。但光学遥感影像对于一些正在发生缓慢变形但又无明显变形迹象的滑坡无法准确识别,且识别出来的滑坡也可能是处于稳定状态的古(老)滑坡; 因此在光学遥感影像识别的基础上,利用多景Sentinel-1A雷达影像数据,基于SBAS-InSAR技术对研究区的正在变形的滑坡隐患进行早期的识别。
2.1 光学遥感影像目视解译
图2
正在发生缓慢变形的滑坡人工目视解译识别相对困难,对于有明显变形迹象的滑坡,可以通过对比多期次高分辨率光学遥感影像的色调和纹理的微变化差异,其解译标志为后缘拉裂缝,前缘小规模滑塌、崩塌等迹象(图3); 而对于没有明显变形迹象的滑坡,则可以采用InSAR技术进行识别。
图3
图3
H37滑坡隐患点光学遥感影像
Fig.3
Optical remote image of H37 landslide hidden danger point
2.2 SBAS-InSAR时序变形探测
2.2.1 采用的数据
表1 Sentinel-1A影像基本参数
Tab.1
参数 | 内容 |
---|---|
卫星 | Sentinel-1A |
模式 | IW |
波段 | C |
波长/cm | VV |
地面分辨率/m | 5×20 |
重访周期/d | 12 |
入射角/(°) | 33.896 |
轨道方向 | 升轨 |
极化方式 | VV |
2.2.2 数据处理流程
SBAS技术处理流程: 首先获取在tA和tB这2个时间段内N+1幅SAR影像,选取其中一幅影像作为主影像,然后将剩余的N幅影像与其配准。根据实际情况设定合理的时间阈值和空间阈值,最终生成M幅差分干涉图,且满足:
对所生成的M幅干涉图中的第j(j=1,2,…,M)幅干涉图,其任意位置的干涉相位可表示为:
式中:
在消除了大气延迟相位、残余地形相位和噪声相位后,最终差分干涉相位可简化为:
最后利用奇异值分解法估算形变速率,在时域积分得形变时间序列。
3 基于光学遥感的滑坡识别结果
通过对高分辨率的光学遥感影像进行人工目视解译,在杆塔附近共识别滑坡隐患点28处,因电网路线较长,将其分为4段进行展示(图4)。
图4
图4
基于光学遥感人工解译识别的滑坡隐患分布图
Fig.4
Distribution map of landslide hidden danger based on manual interpretation and identification of optical remote sensing
4 基于SBAS-InSAR的滑坡隐患探测结果
通过SBAS-InSAR技术对Sentinel-1A升轨道数据处理后,获得研究区在2020年2月—2022年3月期间地表年平均形变速率,最大形变速率为102.5 mm/a。由于研究区范围过大,将其分为4块区域展示(图5)。
图5
图5
哨兵1号A升轨年均形变速率图及识别滑坡隐患分布图
Fig.5
Maps of annual deformation rate of Sentinel-1A Ascent and distribution of landslide hidden danger dentification
图5中,形变量正负值在雷达视线上有不同意义,负值代表远离卫星、正值代表靠近卫星。形变集中区域可能存在滑坡隐患,也可能是由于地表自然沉降、风化作用以及人为建造改造活动导致的形变; 同时,可能受技术水平、气候变化、密集植被覆盖等因素影响,也会形成一些无关形变集中区。因此,结合人工目视解译以及野外调查,排除其他地表运动、各类因素造成的干扰形变项,再从中识别滑坡隐患,最终共识别滑坡隐患点27处,其中15处滑坡隐患点与光学遥感影像解译相符。
5 重点滑坡隐患分析
通过光学遥感影像人工目视解译和SBAS-InSAS技术对研究区滑坡隐患点的探测,共识别滑坡隐患点40处,对电网正常运营有影响的隐患点共7处(表2),其中2处对电网正常运营的风险较大。H08,H10,H18和H36为2种技术手段共同识别出; H16,H23和H24为光学遥感影像单独识别出。
表2 对电网正常运营有影响的滑坡隐患
Tab.2
编号 | 经度 | 纬度 | 长度/m | 宽度/m | 面积/m2 | 遥感解译信息 |
---|---|---|---|---|---|---|
H08 | 107°16'53.58″ | 30°40'58.33″ | 240 | 80 | 17 586.53 | 土质滑坡,前缘存在局部滑塌,公路裂缝明显,杆塔处于滑坡内部后缘,该滑坡一旦发生过大变形影响电网正常运营的风险较大 |
H10 | 107°17'50.83″ | 30°40'57.79″ | 937 | 950 | 668 049.35 | 古滑坡堆积体,由2个子滑坡组成,后缘陡峭,前缘鼓起,后缘有变变形迹象,杆塔位于后缘,影响电网正常运营的风险较大 |
H16 | 107°20'30.33″ | 30°41'49.56″ | 172 | 445 | 69 764.23 | 古滑坡堆积体,两侧存在冲沟,杆塔存在于滑坡内部,但变形迹象不明显,有影响电网正常运营的风险 |
H18 | 107°21'1.76″ | 30°41'57.51″ | 358 | 175 | 64 833.38 | 土质滑坡,局部存在滑塌现象,但整体无明显变形,杆塔处于滑坡后缘,有影响电网正常运营的风险 |
H23 | 107°34'18.04″ | 30°45'49.50″ | 800 | 381 | 218 886.36 | 土质滑坡,前缘存在局部小型滑塌,整体无明显变形,InSAR识别,杆塔位于滑坡中后部,影响电网正常运营的风险较大 |
H24 | 107°34'27.97″ | 30°46'2.87″ | 445 | 627 | 228 396.85 | 土质滑坡,两侧以冲沟为界,前缘存在局部滑塌,杆塔位于滑体中后部,有影响电网正常运营的风险 |
H37 | 107°56'31.03″ | 30°43'35.50″ | 123 | 232 | 24 831.82 | 土质滑坡,前缘位于河中,前缘有明显的变形迹象,杆塔位于中部,有影响电网正常运营的风险 |
5.1 H08滑坡
图6
图6
H08滑坡隐患点光学遥感影像
Fig.6
Optical satellite images of H08 landslide hidden danger point
图7
图8
通过Sentinel-1A升轨方向形变速率结果发现(图9),该滑坡在2020—2022年间LOS方向最大的形变速率为36.5 mm/a,形变主要集中区域为滑坡的前缘与中部。在该滑坡上选取3个特征点P1、P2和P3,获得其时间序列的形变曲线(图10),发现该滑坡堆积体整体处于一个缓慢变形的状态,其中处于滑坡前缘的P1处累积变形形变量为82.2 mm,处于滑坡中部的P2处累积变形形变量为65.4 mm,处于滑坡后缘的P3处累积形变量为28.6 mm。结合该滑坡的形变速率与3个特征点的累积形变量分析发现,该滑坡前缘变形要明显大于中部和后缘,且结合该滑坡前缘较陡、后缘较缓的地形,推测该滑坡堆积体为牵引式滑坡。
图9
图9
H08滑坡隐患点形变速率图
Fig.9
Map of deformation rate of H08 landslide hidden danger point
图10
图10
H08滑坡隐患特征点时序形变曲线
Fig.10
Time series deformation curve of characteristic points of H08 landslide hidden danger
5.2 H10滑坡
该滑坡位于达州市大竹县乌木镇大广子村,E107°17'36.79″,N30°40'54.89″。后缘呈圈椅状,前缘高程527 m,后缘高程829 m,相对高差302 m,横宽约950 m,纵长约937 m,面积约6.7×105 m2,主滑方向为254°,两侧以陡坎为界(图11)。综合InSAR结果与光学遥感影像,推测该滑坡为一古滑坡堆积体,由Ⅰ及Ⅱ2个子区组成,后缘较陡,坡度25°~30°,中部前缘较缓10°~15°。该滑坡在2020—2022年间LOS方向最大的形变速率为66.1 mm/a(图12)。在滑坡堆积体Ⅰ区上选取P1,P2,P3这3个特征点,获得其时间序列的形变曲线(图13),其中P1处累积形变量为58.6 mm,P2处累积形变量为32.8 mm,P3处累积形变量为19.8 mm,在2021年4月后形变速度加快,推测是受到降雨影响,使得该滑坡变形加剧。结合形速率结果以及3个特征点的累积形变量发现,滑坡堆区后缘变形明显大于前缘,且结合其前缘有鼓包,后缘陡峭的地形,推测Ⅰ区为推移式滑坡。在滑坡堆积体Ⅱ区上选取P4,P5,P6这3个特征点,获得其时间序列的形变曲线(图14),P4处累积形变量为60.2 mm,P5处累积形变量为36.8 mm,P6处累积形变量为21.8 mm,目前整体处于匀速变形阶段,发现Ⅱ区同Ⅰ区一样后缘的变形明显大于前缘,但是Ⅰ区后缘变形更加剧烈,位于Ⅰ区后缘的杆塔则可能由于后缘不断变形而发生失稳,所以应对该滑坡堆积体变形情况进行重点监测,以保证电网的正常运营。
图11
图11
H10滑坡隐患点光学遥感影像
Fig.11
Optical satellite images of H10 landslide hidden danger point
图12
图12
H10滑坡隐患点形变速率图
Fig.12
Deformation rate Map of H10 landslide hidden danger point
图13
图13
H10滑坡隐患Ⅰ区特征点时序形变曲线
Fig.13
Time series deformation curves of characteristic points of H10 landslide hidden danger point in zone Ⅰ area
图14
图14
H10滑坡隐患Ⅱ区特征点时序形变曲线
Fig.14
Time series deformation curves of characteristic points of H10 landslide hidden danger point in zone Ⅱ area
6 结论
1)本文采用光学遥感影像人工目视解译与SBAS-InSAR技术相结合的方法,辅以现场调查,实现了对川渝电网“黄万线”杆塔附近地区滑坡隐患的早期识别。利用光学遥感影像在杆塔附近共识别滑坡隐患点28处; 再基于SBAS-InSAR技术,利用Sentinel-1A影像数据在杆塔附近探测出滑坡隐患点27处,2种技术手段共识别滑坡隐患40处,其中15处滑坡隐患为2种技术手段共同识别,有7处隐患点威胁电网杆塔安全、2处风险较大。
2)通过光学遥感影像对研究区进行人工目视解译所识别的滑坡隐患点多为已经发生过滑动的老(古)滑坡,或者有明显形变迹象的滑坡,但对于没有明显形变迹象的滑坡隐患点很难识别出来; SBAS-InSAR对形变探测精度高,可以对缓慢变形的滑坡隐患进行有效探测,但是对于突发型滑坡,由于变形速度太快,则无法进行有效探测。因此,综合利用光学遥感影像目视解译法与SBAS-InSAR技术,可以更加准确、全面地对研究区进行滑坡隐患点进行识别,这对于滑坡隐患早期识别具有借鉴和指导意义。
3)针对植被覆盖率较高的西南地区,采用C波段的Sentinel-1A数据成像时间间隔过长会出现较为严重的失相干现象,所以为了保持较好相干性,应控制好影像的时间间隔。
参考文献
时序InSAR边坡滑移监测下输电杆塔失稳机理研究
[J].
Research on the instability mechanism of transmission towers under time-series InSAR slope slip monitoring
[J].
基于时序InSAR技术的滑坡隐患识别对比研究
[J].
A comparative study of the identification of hidden landslide hazards based on time series InSAR techniques
[J].
遥感技术在滑坡灾害解译中的应用——以陕西省延安市子长县为例
[J].
Application of remote sensing technology to interpretation of landslide disaster:A case study of Zichang County,Shanxi Province
[J].
无人机数字摄影测量技术在滑坡灾害调查中的应用研究
[J].
Research on the application of unmanned aerial vehicle digital photogrammetry in landslide disaster investigation
[J].
四川茂县叠溪镇新磨村滑坡特征与成因机制初步研究
[J].
The Xinmocun landslide on June 24,2017 in Maoxian,Sichuan:characteristics and failure mechanism
[J].
基于InSAR和光学遥感的贵州鬃岭采煤滑坡识别与危险性评价
[J].
Identification and risk assessment of coal mining-induced landslides in Guizhou Province by InSAR and optical remote sensing
[J].
Land Subsidence in Tianjin for 2015 to 2016 Revealed by the Analysis of Sentinel-1A with SBAS-InSAR
[J].DOI:10.1117/1.JRS.11.026024 URL [本文引用: 1]
光学遥感与InSAR结合的金沙江白格滑坡上下游滑坡隐患早期识别
[J].
Early detection of landslides in the upstream and downstream areas of the baige landslide,the Jinsha River based on optical remote sensing and InSAR technologies
[J].
基于时间序列InSAR分析的西部山区滑坡灾害隐患早期识别——以四川丹巴为例
[J].
Early detection of landslide hazards in mountainous areas of West China using time series SAR interferometry:A case study of Danba, Sichuan
[J].
基于天-空-地一体化的重大地质灾害隐患早期识别与监测预警
[J].
Integrated space-air-ground early detection, monitoring and warning system for potential catastrophic geohazards
[J].
对地质灾害隐患早期识别相关问题的认识与思考
[J].
Understanding and consideration of related issues in early identification of potential geohazards
[J].
地质灾害早期识别与监测预警中的综合遥感应用
[J].
Comprehensive application of remote sensing application in early identification,monitoring and early warning of geological disasters
[J].
重大地质灾害隐患早期识别中综合遥感应用的思考与建议
[J].
Early identification of serious geological hazards with integrated remote sensing technologies: Thoughts and recommendations
[J].
重大地质灾害隐患早期识别中综合遥感应用的思考
[J].
Thoughts on the application of comprehensive remote sensing in the early identification of hidden dangers of major geological disasters
[J].
基于SBAS-InSAR的新疆哈密砂墩子煤田开采沉陷监测与反演
[J].
SBAS-InSAR based monitoring and inversion of surface subsidence of the Shadunzi Coal Mine in Hami City,Xinjiang
[J].
基于SBAS-InSAR,技术的成汶高速汶川段滑坡易发区选线研究
[J].
Route selection of landslide prone area in Wenchuan section of Chengdu Wenchuan expressway based on SBAS-InSAR
[J].
利用时序InSAR进行毕节市潜在滑坡识别与形变监测
[J].
DOI:10.13474/j.cnki.11-2246.2022.0183.
[本文引用: 1]
贵州省是全国地质灾害多发、频发最为严重的省份之一,还具有很强的隐蔽性、突发性和不可预见性。为加强贵州纳雍、黔西两县级区域内潜在滑坡的监测与评估,本文利用2017年3月—2020年7月Sentinel-1 SAR数据基于SBAS-InSAR技术获取了纳雍、黔西两县级区域内的时序形变结果。结果表明,在研究区域内共发现60处疑似滑坡点,在监测时段内,滑坡最大累积形变量达23 cm。最终通过区域内的滑坡点进行实际踏勘,表明本文方法可为贵州省地质灾害防治和管理等工作提供参考。
Potential landslides deformation monitoring in Bijie City with InSAR time series
[J].
DOI:10.13474/j.cnki.11-2246.2022.0183.
[本文引用: 1]
Guizhou province is one of the provinces with frequent and most serious geological disasters in China, featured with strong concealment, abruptness and unpredictability. In order to strengthen the monitoring and evaluation of potential landslides in Nayong and Qianxi county, the surface deformation of Nayong and Qianxi county are obtained using Sentinel-1 SAR data from March 2017 to July 2020 and SBAS-InSAR technology. The results show that a total of 60 suspected landslide locations are found in the study area, and the maximum cumulative deformation variable of landslide is up to 23 cm in the monitoring time. Through the field survey of landslide locations in the study region, it shows that this method can provide reference for the prevention and management of geological disasters in Guizhou province.
基于SBAS-InSAR技术和Logistic模型的矿区沉降动态预测方法
[J].
Method for dynamic prediction of mining subsidence based on the SBAS-InSAR technology and the logistic model
[J].
A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms
[J].DOI:10.1109/TGRS.2002.803792 URL [本文引用: 1]
/
〈 |
|
〉 |
