联合InSAR和质量守恒法估计西藏贡觉地区雄巴滑坡厚度  

Thickness Estimation of Xiongba Landslide in Gongjue Area of Tibet,China by Combining InSAR and Mass Conservation Method

作　　者：朱武[1,2,3] 杨璐瑶 张金敏 李振洪 许炫宇[1] ZHU Wu;YANG Lu-yao;ZHANG Jin-min;LI Zhen-hong;XU Xuan-yu(School of Geological Engineering and Geomatics,Chang'an University,Xi'an 710054,Shaanxi,China;Key Laboratory of Ecological Geology and Disaster Prevention of Ministry of Natural Resources,Chang'an University,Xi'an 710054,Shaanxi,China;Big Data Center for Geosciences and Satellites,Chang'an University,Xi'an 710054,Shaanxi,China)

机构地区：[1]长安大学地质工程与测绘学院,陕西西安710054 [2]长安大学自然资源部生态地质与灾害防控重点实验室,陕西西安710054 [3]长安大学地学与卫星大数据研究中心,陕西西安710054

出　　处：《地球科学与环境学报》2023年第3期535-547,共13页Journal of Earth Sciences and Environment

基　　金：国家自然科学基金项目(41941019,42074040);陕西省自然科学基础研究计划项目(2023-JC-JQ-24);陕西省科技创新团队项目(2021TD-51);陕西省地学大数据与地质灾害防治创新团队项目(2022)。

摘　　要：滑坡厚度是灾害风险评估和防控的关键参数。传统的滑坡厚度估计方法只能确定滑坡体局部稀疏点的深度,无法反映整个坡体深度分布情况,且存在探测成本高、效率低等问题。基于此,提出了一种联合合成孔径雷达干涉测量(InSAR)和质量守恒法估计平移式滑坡厚度分布的方法,并将其应用于西藏贡觉地区雄巴滑坡。使用时序InSAR技术分别计算Sentinel-1A升轨、Sentinel-1A降轨和ALOS-2升轨影像的一维雷达视线向(LOS)形变速度场,联合SAR成像参数与滑坡空间几何关系构建坡面坐标系,解算得到滑坡沿坡面的三维形变场;在此基础上,联合InSAR三维形变场和质量守恒准则估计滑坡厚度。实验共收集了覆盖雄巴滑坡的140景Sentinel-1A升轨、138景Sentinel-1A降轨和12景ALOS-2升轨数据。结果表明:雄巴滑坡总运动面积约为5.33 km^(2),当流变参数取值0.5时,滑坡厚度为0~109.57 m,集中分布在9~73 m,滑坡体积约3.12×10^(8)m^(3),估算结果与现场测量结果一致。本文提出的方法可以获取滑坡连续的厚度分布情况,可为灾害风险评估、分析与预防提供关键性数据支撑。The landslide thickness is a key parameter for disaster risk assessment,prevention and control.The traditional landslide thickness estimation method can only determine the depth of the local sparse points of the slope,and can not reflect the depth distribution of the whole landslide,and has the problems of high detection cost and low efficiency.Based on this,a method to estimate the thickness distribution of translational landslide was proposed by combining InSAR and mass conservation method,and was applied to Xiongba landslide in Gongjue area of Tibet.The time-series InSAR technology was used to calculate the one-dimensional LOS deformation velocity fields of SAR images of Sentinel-1A ascending orbit and descending orbit,and ALOS-2 ascending orbit,and the slope coordinate system was constructed by combining the SAR imaging parameters and the spatial geometric relationship of the landslide,and the three-dimensional deformation field of the landslide along the slope was calculated.On this basis,the landslide thickness was estimated by combining InSAR three-dimensional deformation field and mass conservation criterion.The experiment collected 140 views of Sentinel-1A ascending orbit,138 views of Sentinel-1A descending orbit and 12 views of ALOS-2 ascending orbit covering Xiongba landslide.The results show that the total movement area of Xiongba landslide is about 5.33 km^(2);when the rheological parameter is 0.5,the landslide thickness ranges from 0 to 109.57 m,distributes in 9 to 73 m,and the landslide volume is about 3.12×10^(8)m^(3),the estimated results are consistent with the field measurement results.The method proposed can obtain the continuous thickness distribution of landslides and provide key data support for disaster risk assessment,analysis and prevention.

关 键 词：INSAR技术 厚度反演 质量守恒法 三维形变场 坡面坐标系 空间几何参数 流变参数 雄巴滑坡 

分 类 号：P237[天文地球—摄影测量与遥感] P642.22[天文地球—测绘科学与技术]