Application of image-free control UAV aerial survey technology in emergency treatment of landslide-debris flow disaster in Lijie north hill，Zhouqu County
-
摘要: 近年来，无人机航测技术在地质灾害防治领域广泛应用，特别是在突发性地质灾害应急处置中优势显著。免像控无人机航测技术能够快速获取测绘数据，跟传统航测方法相比，极大的节约了时间和人力成本。文中利用大疆精灵Phantom 4 RTK无人机，对甘肃省舟曲县立节北山滑坡-泥石流灾害进行免像控航测数据获取和处理。根据航测成果快速获取了滑坡形态特征、泥石流沟道特征等相关参数, 分析了滑坡变形迹象，并应用实地测量数据进行免像控航测精度评价分析。研究结果表明：滑坡形态特征、泥石流沟道特征、变形迹象均与实地相符，检核点精度满足《低空数字航空摄影测量内业规范》（GH/Z 3003—2010）数据生产规范，成果数据可用于灾害体变形特征分析、几何参数快速提取，辅助分析评估险情、灾情。对突发性地质灾害应急处置中无人机测绘数据的快速获取工作具有一定的借鉴意义。Abstract: In recent years, UAV aerial survey technology has been widely used in the field of geological disaster prevention, especially in the emergency disposal of abrupt geological hazard. Image free UAV aerial survey technology can quickly obtain surveying and mapping data. Compared with traditional aerial survey methods, it greatly saves time cost and labor cost. In this paper, the DJI Phantom 4 RTK UAV is used to acquire and process the image free aerial survey data of landslide and debris flow disaster in Lijie North Hill landslide-debris flow disaster in Zhouqu County, Gansu Province. According to the aerial survey results, the relevant parameters such as landslide morphological characteristics and debris flow gully characteristics are quickly obtained, the signs of landslide deformation are analyzed, and the field survey data are used to evaluate and analyze the accuracy of image free aerial survey.The results show that the deformation signs, morphological characteristics and debris flow gully characteristics of landslide are consistent with the field, and the accuracy of verification points meets the data production specification of low altitude digital aerial photogrammetry (GH/Z 3003—2010). The result data can be used for the analysis of deformation characteristics of disaster body, rapid extraction of geometric parameters, and auxiliary analysis and evaluation of danger and disaster, It has a certain reference significance for the rapid acquisition of UAV mapping data in the emergency disposal of abrupt geological hazard.
-
0. 引言
一般来讲，滑坡的时空分布受限于内在控制因子和外在诱发因子[1]。对于三峡库区滑坡的内在控制因子，学界普遍认为滑坡受工程地质岩组影响较大，如侏罗系红层地区滑坡广泛发育[2-3]，也有学者分析了三峡库区滑坡与构造活动的关系[4-5]，以及斜坡结构类型与滑坡成因机制及空间分布的关系[6-7]。
外在诱发因子主要有降雨、地震、库水位升降以及人类工程活动等。大量案例已表明，随着人类活动对自然环境改造的加剧，滑坡发生频率日渐增高，人类活动已逐渐成为滑坡变形失稳的主要诱因[8]。例如，三峡库水位升降导致大量涉水老、古滑坡发生变形复活，代表性的有巴东县黄土坡滑坡[9]；大规模采矿导致的滑坡以2009年6月5日重庆武隆鸡尾山发生的大型滑坡-碎屑流为代表，体积7.00×106 m3，造成74人死亡[10]；还有部分滑坡是由人类工程建设切坡所致[11-12]，2020年7月13日发生在重庆武隆S529省道旁的白马镇牛儿湾滑坡，体积超过1.80×106 m3，造成4栋房屋倒塌，93人紧急撤离。
除了对单体滑坡的研究外，近年来大量学者从区域层面探讨了滑坡时空分布的驱动因子。Froude等[13]对2004—2016年全球灾难性滑坡的分析发现，在不断增强的城市建设、非法采矿和肆意削坡的扰动下，滑坡数量与日俱增，人类活动极有可能超过降雨成为滑坡的第一诱因；Li等[14]研究了2003—2014年全国28省市土地城市化与滑坡的关系，发现滑坡数量随城市建成区面积和路网密度的增大而增加，尤其以我国西部地区最为显著；王新胜等[15]针对三峡库区重庆主城段1703个历史滑坡点的研究发现，人口密度位居滑坡影响因子之首。这一观点在Li等[16]的研究中也得到证实，除库水位升降外，距道路的距离和土地利用类型等人类活动因子，对三峡库区滑坡分布的影响也已超过降雨因子。
不单是针对三峡库区，国内外其他研究也已表明，在快速城镇化的背景下，不断增强的人类活动是控制滑坡分布的关键要素[17]。这一认识在经济发达的城市区域是适用的，但是在城镇化率较低、人类活动强度较小的农村县城是否可行？为此，文中选择三峡库区腹心的武隆区作为研究区，该区的城镇化率较低（仅为全市的67.4%，2018年），人类活动强度位列三峡库区（重庆段）倒数第二[18]。以重庆武隆1991—2015年330处历史滑坡为样本点，运用重力模型和标准差椭圆模型描述滑坡的时空分布格局，利用地理探测器进行驱动力分析，试图解释该格局的演化机制，得出欠发达县域内的人类活动也是影响滑坡分布关键性因素的观点，其结论可为三峡库区滑坡灾害减灾防灾规划提供依据。
1. 研究区概况与数据来源
1.1 研究区概况
武隆区地处重庆市东南部乌江下游、三峡库区腹心（图1）。地势东北高、西南低，以山地为主，相对高差700～1000 m。气候属亚热带湿润季风气候，雨量丰沛、四季分明，多年平均降水量多在1000 mm以上，雨量在时间分配上不均匀，4—9月降雨占全年总量的70%以上。区内地质构造以NE-NNE走向的褶皱为主，断层不发育，属新华夏构造体系和南北经向构造体系——川黔南北构造带。出露地层全为沉积岩，岩性以碳酸盐岩、砂岩、页岩和泥岩为主。
![]()
图 1 研究区概况注：滑坡数据来自“重庆市武隆县地质灾害排查项目”；该图基于国家标准地图服务网的审图号为GS(2020)4630号的标准地图、重庆市标准地图服务网的审图号为渝S(2019)071号的标准地图制作，底图无修改，下同。Figure 1. Research area武隆境内地势起伏大、山高谷深、河流密布、降雨量大且时间集中，为滑坡提供了适宜的孕灾环境，武隆一直都是三峡库区地质灾害比较严重的地区，在重庆市地质灾害危险性分区中属于武隆—南川重灾带[19]。加之近年来，区内社会经济发展水平不断提高、城市建设规模逐步扩大，人类工程活动如道路开挖、矿山开采和水电站等建设活动日益频繁，人类对地表环境的干扰程度也不断加强，导致重大滑坡灾害时有发生。
长江干流库岸滑坡受库水位变动的影响大，这已是学界的共识[1,16]。武隆位于长江二级支流乌江的下游（图1），三峡蓄水前武隆站水位变化区间约为169～185 m，蓄水后水位变化区间缩小为169～177 m，且季节间变动幅度小（图2），区内滑坡受库水位升降影响较小，同时文中重在揭示城镇化背景下人类生活生产活动对滑坡分布的影响，因此可排除库水位升降的干扰。
![]()
图 2 武隆站与三峡大坝（坝前）的水位对比注：数据来自长江水文网。Figure 2. Water level comparison between Wulong Station and Three Gorges Dam (before the dam)1.2 影响因子选取
结合武隆区自然环境条件及社会经济发展状况，文中将影响滑坡分布的多个指标分为降雨、地质条件、地形地貌和人类活动等4个方面共14种因子[20-21]（表1）。
表 1 滑坡分布的影响影响因子分类及说明Table 1. Classification and description of impact factors affecting landslide distribution类别 名称 年份 来源 描述 降雨 年均降雨量 1990—2015 http://www.resdc.cn 1 km×1 km格网 地质条件 工程地质岩组 2002 重庆市工程地质图 栅格数据，1∶50万 距背向斜轴区距离 − 重庆市1∶50万地质图 栅格数据 距断层距离 − 重庆市1∶50万地质图 栅格数据 地形地貌 土壤类型 2000 http://www.resdc.cn 1 km×1 km格网 土壤侵蚀强度 2000 http://www.resdc.cn 1 km×1 km格网 地表起伏度 2003 由DEM提取 30 m×30 m格网 坡度 2003 由DEM提取 30 m×30 m格网 距河流距离 2014 http://www.geodata.cn 矢量数据 NDVI 2000—2015 http://www.resdc.cn 1 km×1 km格网 人类活动 人口密度 2000—2015 http://www.resdc.cn 1 km×1 km格网 GDP 2000—2015 http://www.resdc.cn 1 km×1 km格网 距道路距离 2000—2015 http://www.geodata.cn 武隆区1∶25万道路交通数据 距采矿点距离 − 全国矿产地数据库 矢量数据 （1）降雨：降雨往往是滑坡的直接诱发因素，由于强降雨，斜坡土体含水量增大，斜坡岩土体静、动水压力急剧增大，降低滑动面的抗剪强度[22]，增大了下滑力，从而诱发滑坡。
（2）地质条件：选取工程地质岩组、距背向斜轴区距离、距断层距离3个因子来衡量该区的地质条件。地层岩性影响深部岩土体的物理属性，斜坡上的岩体往往在软岩、软层的影响下容易出现塑性形变，不同岩性接触带容易出现片理风化作用，形成强烈的错动带和风化带[23]；背向斜轴区构造应力相对集中[24]，裂隙分布密集，层间常有构造滑移面，且在局部地段存在顺向坡；断层附近岩土体结构不稳定，在外力扰动作用下易发生滑坡。
（3）地形地貌：地表起伏度可反映地面相对高差，影响着地表物质的侵蚀、搬运、堆积等过程，在很大程度上决定了滑坡灾害的易发程度；斜坡的坡度越大，土体下滑力越大，斜坡稳定性越低[1]；河流、沟谷的两侧冲刷侵蚀作用强烈，斜坡前缘土体被带走，形成临空面，容易发生滑坡。土壤类型、土壤侵蚀强度和NDVI等6个因子可以综合衡量该区地表土体的抗蚀性和抗冲性[25]。
（4）人类活动
武隆区破坏地质环境的人类工程活动主要体现在不合理切坡、填土、盲目采矿和矿山弃渣不合理堆放以及滥垦滥伐，因此选取距道路距离、距采矿点距离、人口密度和GDP等4个因子作为衡量该区人类活动强度的指标[16,26-27]。
1.3 数据处理
研究区内共有330个历史滑坡灾害点，数据来自“重庆市武隆县地质灾害排查项目”调查成果，包括20世纪80年代以来武隆境内滑坡的类型、发生时间、位置、主要诱因、威胁人户及财产等信息。文中首先利用1991—2015年的330处历史滑坡，分析该时段滑坡的时空分布格局。其次，考虑到GDP、人口、道路等相关数据的匹配性，选取2001—2015年的287处历史滑坡分析该时间段滑坡空间分布的驱动因子。上述14种因子中，对于连续型因子进行离散化，并采用自然断点法分成9级[15,28]；对于离散型因子则保持其原有分级，部分因子见图3。
2. 研究方法
采用趋势分析、核密度分析、重力模型和标准差椭圆模型分析滑坡的时空分布格局，利用地理探测器对不同时段滑坡分布的驱动因子进行分析，揭示其演化机制。
2.1 滑坡时空分布格局分析
2.1.1 重力模型
重心变化是研究海量时空数据演变趋势的有效手段[29]，“重心”可表征地理要素的时空分布特征，重心移动方向指示了空间现象的“高密度”部位，而重心移动的距离反映了要素变化幅度的空间差异。定义武隆滑坡的重心坐标分别为：
$$ \overline{X}=\frac{{\displaystyle \sum _{i=1}^{n}{x}_{i}}}{n}，\overline{Y}=\frac{{\displaystyle \sum _{i=1}^{n}{y}_{i}}}{n} $$ (1) 式中：
$ \overline X $ 、$ \overline Y $ ——滑坡的重心坐标；xi、yi——第i个滑坡的空间坐标；
n——滑坡总数。
2.1.2 标准差椭圆模型
标准差椭圆可用于度量地理要素的中心趋势、离散和方向趋势等空间分布特征[30]。标准差椭圆的中心点表示所有要素的中心位置，长半轴表示地理要素分布的方向性，短半轴表示地理要素分布的范围。因此通过绘制滑坡空间分布的标准差椭圆，比较不同时间段的椭圆变化特征，可以揭示滑坡分布的年代际空间动态过程和总体趋势。计算公式如下：
$$ C = \frac{1}{n}\left( \begin{gathered} \sum\limits_{i = 1}^n {{{\widetilde {{x_i}}}^2}} \hfill \\ \sum\limits_{i = 1}^n {\widetilde {{x_i}}} \widetilde {{y_i}} \hfill \\ \end{gathered} \right.\left. \begin{gathered} \sum\limits_{i = 1}^n {\widetilde {{x_i}}} \widetilde {{y_i}} \hfill \\ \sum\limits_{i = 1}^n {{{\widetilde {{y_i}}}^2}} \hfill \\ \end{gathered} \right) \text{，} \left\{ \begin{gathered} \widetilde {{x_i}} = \left( {{x_i} - \overline {x'} } \right) \hfill \\ \widetilde {{y_i}} = \left( {{y_i} - \overline {y'} } \right) \hfill \\ \end{gathered} \right. $$ (2) 式中：C——标准差椭圆轴长；
n——滑坡总数；
$\widetilde {{x_i}}$ 、$\widetilde {{y_i}}$ ——平均中心和xi、yi坐标的差；xi、yi——第i起滑坡的坐标；
$\overline {x'} $ 、$\overline {y'} $ ——滑坡分布的平均中心。2.2 空间分布驱动因素探测
地理探测器是基于地理空间分异立论所提出的探测空间分异性、揭示其背后驱动因子的一种统计学方法[31]。其核心思想认为：若某个自变量对因变量有重要影响，那么二者在空间分布上应该具有相似性。地理探测器已广泛用于台风[32]、山洪[33]、滑坡[34]和泥石流[35]等自然灾害领域。地理探测器包括因子探测、风险探测、交互探测和生态探测。
文中采用因子探测和交互探测对滑坡空间分布的驱动因子进行分析，以降雨、地质条件、地形地貌和人类活动等4个方面共14种因子作为自变量X，以滑坡核密度值作为因变量Y，利用不同时间段的q值来反映驱动力的演化过程，通过探测各变量之间的空间分异程度和交互作用影响，从而找到影响滑坡空间分布的关键影响因子。
（1）因子探测
因子探测器通过因子解释力衡量自变量对因变量的贡献大小，进而判断某自变量是否为因变量空间分异性的驱动因子，以及所占权重。通过比较因变量在自变量各类别中的方差与全局方差的关系来计算因子解释力的大小，计算公式如下：
$$ q = 1 - \frac{{\displaystyle\sum\limits_{h = 1}^L {{N_h}\sigma _h^2} }}{{N{\sigma ^2}}}{\text{ = }}1 - \frac{{SSW}}{{SST}} $$ (3) 式中：q——解释力，为某影响因子在多大程度上解释了 滑坡的空间分布，q的值域为[0, 1]，q值越大， 表示影响因子对滑坡空间分布的解释力越强；
h——变量因子分类数，
$h= 1,\cdots, L $ ；Nh、N——类型h、全区的样本数；
$ \sigma_h^2$ 、$ \sigma^2$ ——类型h、全区的方差；SSW——类型方差之和；
SST——全区总方差。
（2）交互探测
交互探测器分别计算和比较各单因子q值及两因子叠加后的q值，以判断两因子是否存在交互作用以及交互作用的强度等[36]。X1∩X2是将X1和X2进行空间叠加后形成的新的空间因子，利用表2的判别依据，确定因子间的交互作用类型。
表 2 交互探测关系对应表Table 2. Interaction detection relationship correspondence table判据 交互作用 q(X1∩X2) < min(q(X1), q(X2)) 非线性减弱 Min(q(X1), q(X2)) < q(X1∩X2) < max(q(X1), q(X2)) 单因子非线性减弱 q(X1∩X2) > max(q(X1), q(X2)) 双因子增强 q(X1∩X2) = q(X1) + q(X2) 独立 q(X1∩X2) > q(X1) + q(X2) 非线性增强 3. 滑坡灾害的时空演变格局
3.1 时间变化特征
以武隆境内330个历史滑坡点数据为基础，分析1991—2015年滑坡逐年变化趋势及其演化过程（图4）。武隆滑坡发生频次按累积曲线变化特征可划分为1991—1997年、1998—2007年和2008—2015年3个阶段性变化过程。在1991—1997年处于低发阶段，滑坡发生频次为0～7起/年，平均每年发生2.28起；在1998—2007年进入快速增长阶段，发生滑坡277起，占总数的83.94%，滑坡发生频次为4～84起/年，平均每年发生27.7起；在2008—2015年处于缓慢增长阶段，滑坡发生频次为2～8起/年，平均每年发生4.63起。3个异常高值均出现在快速增长阶段，分别为1998年、2003年和2007年，对应的滑坡数量为18起、84起和53起，共计155起，占总数的46.97%。
![]()
图 4 滑坡累积量与降雨量的对应关系Figure 4. Correspondence between landslide accumulation and rainfall从滑坡累积量与年均降雨量的关系（图4）来看，1991—1997年的年均降雨量仅为973.38 mm，1998—2007年增至1325.04 mm，增幅超过30%，滑坡累积曲线斜率较上一阶段明显变大，与降雨量的变化趋势基本一致，2008—2015年间年均降雨量为1306.79 mm，与1998—2007年相比波动较小，然而滑坡累计曲线斜率却明显减缓。有研究表明，1300 mm是三峡库区地质灾害暴发的年降雨量临界值，在年均降雨量超过1300 mm后，地质灾害累积数与降雨量关系曲线斜率明显变缓[37]。因此，我们初步判断认为，降雨对于滑坡空间分布的影响在2007年以后存在一定程度的降低，对此还需结合滑坡发生机理和相关数据，进一步论证分析。
3.2 空间变化特征
3.2.1 滑坡分布密度分析
针对武隆境内1991—2015年的历史滑坡点，使用ArcGIS 计算核密度估计值以展示其聚集区域与聚集程度。据核密度估计结果（图5），滑坡的空间分布存在明显差异，呈现出西、中、东3个聚集区域：（1）西部：平桥镇—鸭江镇的带状聚集；（2）中部：羊角镇—白马镇—长坝镇的环形聚集；（3）东部：巷口镇—火炉镇的带状聚集。其中，白马镇西北部和巷口镇中部的灾害密度值>1.5，平桥镇、长坝镇、羊角镇、火炉镇西部和仙女山镇东南部等地的灾害密度值为1.0～1.5，鸭江镇、庙垭乡、白云乡、江口镇以及沧沟乡等地的灾害密度值为0.5～1，其余地区滑坡发生起数较少。
3.2.2 滑坡重心及其变化趋势
将1991—2015年的武隆滑坡数据以5年为间隔，采用重力模型分析滑坡灾害重心的移动轨迹和移动距离（图6），利用标准差椭圆分析滑坡方向及分布（表3），以揭示滑坡空间分布的变化特征。
表 3 1991—2015年武隆滑坡标准差椭圆参数变化Table 3. Variation of standard deviation ellipse parameters of Wulong landslide from 1991—2015年份 1991—1995 1996—2000 2001—2005 2006—2010 2011—2015 重心坐标 107°35′E
29°21′N107°32′E
29°22′N107°37′E
29°23′N107°43′E
29°21′N107°41′E
29°22′N移动方向 − 西北 东北 东南 西北 方向角度/（°） 101.4 97.5 88.3 94.4 77.7 长半轴/km 25.4 25.3 27.3 23.1 24.4 短半轴/km 6.5 11.3 11.4 13.4 11.9 偏心率 0.97 0.89 0.91 0.81 0.87 ![]()
图 6 历史滑坡灾害重心及方向分布特征Figure 6. Historical landslide hazard center of gravity and directional distribution characteristics1991—2015年武隆滑坡的重心在29°21′N～29°23′N、107°32′E～107°43′E之间移动。1996—2000年滑坡重心由白马镇向西北移动至羊角镇中部，2001—2005年向东北方向移动至羊角镇西部，2006—2010年由羊角镇向东南方向移动至巷口镇与仙女山镇交界处，2011—2015年由巷口镇向西北方向迁回至羊角镇南部。1991—2015年灾害重心移动7.6 km，重心整体向东迁移，其中1996—2000年、2001—2005年、2006—2010年和2011—2015年移动距离分别为6.2 km、8.0 km、11.4 km和4.8 km，共计移动30.4 km，平均每年移动1.2 km。
1991—2015年的标准差椭圆方向在77.7°～101.4°之间，武隆滑坡的空间分布呈现出由西北—东南方向转向东北—西南方向的变化过程。标准差椭圆的长半轴整体呈先增后减的趋势，2001—2005年最长（27.3 km），2006—2010年最短（23.1 km）；标准差椭圆的短半轴整体呈增长趋势，1991—1995年最短（6.5 km），2006—2010年最长（13.4 km）。结合偏心率分析，武隆滑坡空间分布的方向性呈减弱趋势，武隆滑坡的空间分布由聚集走向分散。
4. 滑坡时空分布的驱动力分析
4.1 单因子影响分析
考虑到人口密度、GDP等数据存在一定滞后性及其与灾害数据的匹配性，仅选取2001—2015年的287个历史滑坡点为样本，将其按5年为间隔划分为3个时间段，利用地理探测器定量评价因子对滑坡分布的影响，实现对武隆滑坡时空格局驱动力的演化分析，分析结果见图7。
![]()
图 7 滑坡分布影响因子q值的变化过程Figure 7. The process of changing q value of landslide distribution influence factors从2001—2015年来看，各影响因子对滑坡空间分布的解释力均值最高的前5位依次为：降雨（0.274）、工程地质岩组（0.220）、土壤类型（0.132）、河流（0.120）、GDP（0.116）；而背向斜轴区、断层、NDVI、采矿点、地表起伏度、坡度等因子的解释力均低于0.05，对滑坡空间分布的影响不显著。
从不同时间段来看，2001—2005年、2006—2010年和2011—2015年解释力最高的因子分别为降雨（0.382）、工程地质岩组（0.319）、工程地质岩组（0.204）。相比于2001—2010年，2011—2015年降雨、地质条件因子和地形地貌因子对滑坡空间分布的解释力均有不同程度降低，其中降雨和河流的降幅最为明显，而人类活动因子的解释力则逐渐增强。根据前文（图3、图6）不难看出，滑坡高密度区往往也是人类活动密集区，武隆区2001—2015年间各项社会经济指标均呈现不同幅度增长（表4），表明随着经济的不断发展，人类对地表环境的改造越来越强烈，人类不合理的工程建设活动如切坡、填土、采矿和矿山弃渣不合理堆放以及由于经济发展和人口增长，所带来的人均耕地面积的减少，植被破坏、水土流失等问题，使得人类活动对于滑坡空间分布的影响也越来越显著。
表 4 武隆区2005—2015年社会经济相关统计数据Table 4. Relevant socio-economic statistics of Wulong from 2005 to 2015年份 2005 2010 2015 GDP/亿元 29.88 72.42 131.40 人口数量/万人 40.10 41.27 41.43 公路里程/km 1153 2795 2857 4.2 因子交互影响分析
根据地理探测器交互探测结果，各时间段解释力最高的前10位如表5所示，各因子叠加后其解释力明显提高，表明滑坡的空间分布是多因素共同作用的结果[34]。2001—2015年间，主要以降雨和工程地质岩组与其它因子的交互作用最为显著，其中2001—2005年降雨与土壤侵蚀的交互结果以及2006—2010年工程地质岩组与降雨和河流的交互解释力均超过0.5，表明其对研究区50%以上的滑坡空间分布具有解释力。
表 5 交互探测后的q值Table 5. q-value after interaction detection2001—2005 2006—2010 2011—2015 降雨∩土壤侵蚀 0.505 工程地质岩组∩降雨 0.504 降雨∩工程地质岩组 0.449 降雨∩土壤类型 0.487 工程地质岩组∩河流 0.500 GDP∩工程地质岩组 0.390 降雨∩工程地质岩组 0.462 工程地质岩组∩道路 0.458 GDP∩降雨 0.318 降雨∩NDVI 0.451 工程地质岩组∩
人口密度0.422 人口密度∩
工程地质岩组0.317 降雨∩河流 0.445 工程地质岩组∩GDP 0.416 降雨∩土壤类型 0.310 降雨∩背向斜轴区 0.435 工程地质岩组∩
土壤类型0.402 GDP∩土壤类型 0.302 降雨∩人口密度 0.434 工程地质岩组∩
土壤侵蚀0.396 人口密度∩降雨 0.298 降雨∩GDP 0.431 降雨∩河流 0.392 GDP∩道路 0.293 降雨∩断层 0.403 降雨∩人口密度 0.388 工程地质岩组∩河流 0.288 降雨∩地表起伏度 0.401 工程地质岩组∩NDVI 0.385 降雨∩道路 0.272 该区岩土体均受水力侵蚀影响，2001—2005年有29%的滑坡分布于强度水力侵蚀区和剧烈水力侵蚀区，50%的滑坡分布于中度水力侵蚀区，在土壤侵蚀作用下，地表岩、土体等被剥蚀、破坏和搬运的过程加剧。随着降雨的入渗，土体发生软化形成软滑面，造成土体失稳进而引发滑坡。同时，从工程地质岩组分布来看，区内滑坡点多沿河谷呈带状分布，河谷地带岩性组合以灰岩与泥、页岩岩组为主，此类岩性组合往往存在软弱夹层，因差异风化使斜坡成陡—缓—陡形式，形成凹腔[38]，在后缘加载与强降雨触发双重影响下极易发生滑坡，因而滑坡灾害沿河流走向呈带状密集发育。
需要指出的是人类活动因子与其它因子的交互解释力在2001—2015年间比重逐步上升，其中2006—2010年，工程地质岩组与距道路距离的交互解释力高达0.458，同时，2011—2015年解释力较高的交互结果也多以人类活动因子为主，且所有因子交互后的q值都超过原因子之和，其解释力均呈非线性增长，进一步论证了上述人类活动对于滑坡空间分布的影响越来越明显的观点。
5. 讨论
尽管人们对影响滑坡发生的人为作用并不陌生，但大多属于感性认识，文中通过对武隆历史滑坡中人类活动的影响过程进行量化，有助于理解灾害系统中“人”的地位和作用，有利于形成合理的人地协调观。基于文中的认识，在武隆区未来的滑坡减灾防灾工作中，应当重视人类活动的影响，重点关注城镇建设、道路开挖等活动对自然坡体的扰动，同时在这些人类活动密集的区域内，人员、建筑物等承灾体也相对集中，因此，灾害造成的损失往往都很大。这也提醒我们在滑坡灾害的风险评价，尤其是危险性评价中，需要着重考虑人类活动的相关指标。
同时文中在以下两个方面尚存在不足之处：（1）本文降雨因子仅采用了年平均降雨数据，且由于无法获取研究区地层产状数据，仅以工程地质岩组、断层和背斜向斜轴区作为地质影响因子稍有欠缺，同时人类活动因子的选取也有一定的局限性。（2）仅以滑坡密度表征灾害强度，同时以网格作为评价单元，忽视了滑坡面的整体性和连续性，滑坡与各因子间的内在空间关系和影响机制也未得到良好体现。文中仅对武隆区滑坡驱动因子的演变格局和人类活动影响进行了初步的探讨，上述不足之处将在今后的研究中进一步完善。
6. 结论
（1）武隆1991—2015年的滑坡累计曲线在时间分布上呈现出“缓-陡-缓”的特征，以2008年为界，之前的滑坡发生速率随降雨量的增加而增长，2008—2015年，在降雨量保持稳定的情况下，滑坡发生速率明显减缓。
（2）在空间分布上，武隆历史滑坡呈现西、中、东3个聚集区域；滑坡分布的重心集中于中部的白马镇、羊角镇和东部的巷口镇等地；标准差椭圆分析表明：武隆滑坡的空间分布呈现出由西北—东南方向转向东北—西南方向的变化过程，同时表现出方向性减弱和离散化的趋势。
（3）针对上述滑坡的时空分布格局，利用地理探测器找寻了2001—2005、2006—2010和2011—2015三个时段滑坡分布的驱动因子，结果显示，降雨和“软-硬”互层结构岩组与其它因子的交互作用是影响武隆区滑坡分布的主要影响因素；同时降雨、地质因子和地形地貌等因子对滑坡空间分布的解释力均有不同程度降低，而人类活动因子的解释力不断增强，已逐渐成为影响滑坡分布的关键影响因素之一，是未来减灾防灾重点关注的对象。
-
![]()
图 3 研究区地理位置和高程分布图
Figure 3. Geographical location and elevation distribution map of the study area
-
[1] 陈文君, 陈锁忠, 张磊, 等. 突发性地质灾害应急处置系统的设计与实现[J]. 计算机应用与软件,2014,31(5):6 − 9. [CHEN Wenjun, CHEN Suozhong, ZHANG Lei, et al. Design and implementation of an emergency disposition system on sudden geological disasters[J]. Computer Applications and Software,2014,31(5):6 − 9. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-386x.2014.05.002 [2] 刘传正. 重大突发地质灾害应急处置的基本问题[J]. 自然灾害学报,2006,15(3):24 − 30. [LIU Chuanzheng. Basic problem on emergency disposition of abrupt heavy geological disaster[J]. Journal of Natural Disasters,2006,15(3):24 − 30. (in Chinese with English abstract) DOI: 10.3969/j.issn.1004-4574.2006.03.005 [3] 刘传正, 陈红旗, 韩冰, 等. 重大地质灾害应急响应技术支撑体系研究[J]. 地质通报,2010,29(1):147 − 156. [LIU Chuanzheng, CHEN Hongqi, HAN Bing, et al. Technical support system of emergency response for serious geo-hazards[J]. Geological Bulletin of China,2010,29(1):147 − 156. (in Chinese with English abstract) [4] 陈红旗, 徐永强, 庄茂国, 等. 地质灾害应急支撑体系建设基本问题分析[J]. 中国地质灾害与防治学报,2011,22(4):108 − 111. [CHEN Hongqi, XU Yongqiang, ZHUANG Maoguo, et al. Study on the construction of geological disaster emergency support system[J]. The Chinese Journal of Geological Hazard and Control,2011,22(4):108 − 111. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2011.04.023 [5] 李俊峰, 张小趁, 刘红岩, 等. 突发地质灾害应急数值模拟浅析[J]. 中国地质灾害与防治学报,2016,27(2):154 − 161. [LI Junfeng, ZHANG Xiao chen, LIU Hongyan, et al. Numerical simulation technology analysis of geological disaster emergency[J]. The Chinese Journal of Geological Hazard and Control,2016,27(2):154 − 161. (in Chinese with English abstract) [6] CHOU T Y, YE M L, CHEN Y C, et al. Disaster monitoring and management by the unmanned aerial vehicle technology[J]. International archives of the photogrammetry, remote sensing and spatial information sciences,2010,38(7B):137 − 142.
[7] DOMINICI, ALICANDRO, MASSIMI V. UAV photogrammetry in the post-earthquake scenario: Case studies in L'Aquila[J]. Geomatics Natural Hazards and Risk,2017,8(1):87 − 103.
[8] SALVINI R, MASTROROCCO G, SEDDAIU M, et al. The use of an unmanned aerial vehicle for fracture mapping within a marble quarry (Carrara, Italy): Photogrammetry and discrete fracture network modelling[J]. Geomatics Natural Hazards & Risk,2017,8(1):34 − 52.
[9] 曾涛, 杨武年, 简季. 无人机低空遥感影像处理在汶川地震地质灾害信息快速勘测中的应用[J]. 测绘科学,2009,34(增刊 2):64 − 65. [ZENG Tao, YANG Wunian, JIAN Ji. Applications of remote sensing images processing of unmanned air vehicle on the low attitude in quick acquirement of geological disasters information in Wenchuan earthquake[J]. Science of Surveying and Mapping,2009,34(Sup 2):64 − 65. (in Chinese with English abstract) [10] 高姣姣, 颜宇森, 盛新蒲, 等. 无人机遥感在西气东输管道地质灾害调查中的应用[J]. 水文地质工程地质,2010,37(6):126 − 129. [GAO Jiaojiao, YAN Yusen, SHENG Xinpu, et al. Application of UAV remote sensing in Geologic hazards survey along the project of west-east gas transmission[J]. Hydrogeology & Engineering Geology,2010,37(6):126 − 129. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-3665.2010.06.025 [11] 朱婵莲, 姜啸, 焦龙进. 无人机在地质灾害应急抢险调查中的应用[J]. 中国金属通报,2020(6):184 − 185. [ZHU Chanlian, JIANG Xiao, JIAO Longjin. Application of UAV in geological disaster emergency investigation[J]. China Metal Bulletin,2020(6):184 − 185. (in Chinese with English abstract) [12] 郭晨, 许强, 彭双麒, 等. 无人机摄影测量技术在金沙江白格滑坡应急抢险中的应用[J]. 灾害学,2020,35(1):203 − 210. [GUO Chen, XU Qiang, PENG Shuangqi, et al. Application research of UAV photogrammetry technology in the emergency rescue of Baige landslide[J]. Journal of Catastrophology,2020,35(1):203 − 210. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-811X.2020.01.038 [13] 李晓俊, 朱钱洪, 胡勇, 等. 多旋翼无人机山区城市地质灾害应急响应能力评估[J]. 地理空间信息,2019,17(3):12 − 15. [LI Xiaojun, ZHU Qianhong, HU Yong, et al. Geological disaster emergency response capability assessment of multirotor UAV in mountain city[J]. Geospatial information,2019,17(3):12 − 15. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-4623.2019.03.004 [14] 叶伟林, 宿星, 魏万鸿, 等. 无人机航测系统在滑坡应急中的应用[J]. 测绘通报,2017(9):70 − 74. [YE Weilin, SU Xing, WEI Wanhong, et al. Application of UAV aerial photograph system in emergency rescue and relief for landslide[J]. Bulletin of Surveying and Mapping,2017(9):70 − 74. (in Chinese with English abstract) [15] 董秀军, 王栋, 冯涛. 无人机数字摄影测量技术在滑坡灾害调查中的应用研究[J]. 地质灾害与环境保护,2019,30(3):78 − 85. [DONG Xiujun, WANG Dong, FENG Tao. Research on the application of unmanned aerial vehicle digital photogrammetry in landslide disaster investigation[J]. Journal of Geological Hazards and Environment Preservation,2019,30(3):78 − 85. (in Chinese with English abstract) [16] 吴振宇, 马彦山. 无人机遥感技术在地质灾害调查中的应用[J]. 宁夏工程技术,2012,11(2):133 − 136. [WU Zhenyu, MA Yanshan. Application of remote sensing technology of UAV in the investigation of the geological disaster[J]. Ningxia Engineering Technology,2012,11(2):133 − 136. (in Chinese with English abstract) DOI: 10.3969/j.issn.1671-7244.2012.02.010 [17] 崔红霞, 林宗坚, 孙杰. 无人机遥感监测系统研究[J]. 测绘通报,2005(5):11 − 14. [CUI Hongxia, LIN Zongjian, SUN Jie. Research on UAV remote sensing system[J]. Bulletin of Surveying and Mapping,2005(5):11 − 14. (in Chinese with English abstract) DOI: 10.3969/j.issn.0494-0911.2005.05.003 [18] 晏磊, 吕书强, 赵红颖, 等. 无人机航空遥感系统关键技术研究[J]. 武汉大学学报(工学版),2004,37(6):67 − 70. [YAN Lei, LYU Shuqiang, ZHAO Hongying, et al. Research on key techniques of aerial remote sensing system for unmanned aerial vehicles[J]. Engineering Journal of Wuhan University,2004,37(6):67 − 70. (in Chinese with English abstract) [19] 雷添杰, 李长春, 何孝莹. 无人机航空遥感系统在灾害应急救援中的应用[J]. 自然灾害学报,2011,20(1):178 − 183. [LEI Tianjie, LI Changchun, HE Xiaoying. Application of aerial remote sensing of pilotless aircraft to disaster emergency rescue[J]. Journal of Natural Disasters,2011,20(1):178 − 183. (in Chinese with English abstract) [20] 马泽忠, 王福海, 刘智华, 等. 低空无人飞行器遥感技术在重庆城口滑坡堰塞湖灾害监测中的应用研究[J]. 水土保持学报,2011,25(1):253 − 256. [MA Zezhong, WANG Fuhai, LIU Zhihua, et al. Research of unmanned aerial vehicle remote sensing technique applied to monitor landslip and barrier lake hazard in Chengkou County of Chongqing[J]. Journal of Soil and Water Conservation,2011,25(1):253 − 256. (in Chinese with English abstract) [21] 尹鹏飞, 尹球, 陈兴峰, 等. 无人机航空遥感技术在震后灾情调查中的应用[J]. 激光与光电子学进展,2010,47(11):130 − 134. [YIN Pengfei, YIN Qiu, CHEN Xingfeng, et al. Unmanned aerial vehicle aerial remote sensing techniques and its application on post-earthquake disaster investigation[J]. Laser & Optoelectronics Progress,2010,47(11):130 − 134. (in Chinese with English abstract) [22] 韩文权, 任幼蓉, 赵少华. 无人机遥感在应对地质灾害中的主要应用[J]. 地理空间信息,2011,9(5):6 − 8. [HAN Wenquan, REN Yourong, ZHAO Shaohua. Primary usages of UAV remote sensing in geological disaster monitoring and rescuing[J]. Geoscience Information,2011,9(5):6 − 8. (in Chinese with English abstract) [23] 刘刚, 徐宏健, 马海涛, 等. 无人机航测系统在应急服务保障中的应用与前景[J]. 测绘与空间地理信息,2011,34(4):177 − 179. [LIU Gang, XU Hongjian, MA Haitao, et al. Unmanned aerial system applications and prospects in the protection of emergency services[J]. Geomatics & Spatial Information Technology,2011,34(4):177 − 179. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-5867.2011.04.057 [24] 肖波, 朱兰艳, 黎剑, 等. 无人机低空摄影测量系统在地质灾害应急中的应用研究:以云南洱源特大山洪泥石流为例[J]. 价值工程,2013,32(4):281 − 282. [XIAO Bo, ZHU Lanyan, LI Jian, et al. Application of unmanned aerial vehicle photogram metric system in geological disaster: A case study in Er yuan devastating mudslides of Yunnan Province[J]. Value Engineering,2013,32(4):281 − 282. (in Chinese with English abstract) DOI: 10.3969/j.issn.1006-4311.2013.04.155 [25] 嵇沛, 胡一凡, 王晨. 无人机低空航拍遥感在应急测绘保障中的应用[J]. 科技资讯,2016,14(24):2 − 3. [JI Pei, HU Yifan, WANG Chen. Application of UAV low altitude aerial remote sensing in emergency surveying and mapping support[J]. Science & Technology Information,2016,14(24):2 − 3. (in Chinese with English abstract) [26] 田彦. 无人机遥感技术在地质环境灾害监测中的应用探讨[J]. 环境与发展,2019,31(4):102 − 104. [TIAN Yan. Application of UAV remote sensing technology in geological environmental disaster monitoring[J]. Environment and Development,2019,31(4):102 − 104. (in Chinese with English abstract) [27] 周嵩. 无人机技术在地质灾害防治中的应用[J]. 冶金管理,2020(21):87 − 88. [ZHOU Song. Application of UAV technology in geological disaster prevention and control[J]. Metallurgical Management,2020(21):87 − 88. (in Chinese with English abstract) [28] 郭晨, 许强, 董秀军, 等. 无人机在重大地质灾害应急调查中的应用[J]. 测绘通报,2020(10):6 − 11. [GUO Chen, XU Qiang, DONG Xiujun, et al. Application of UAV photogrammetry technology in the emergency rescue of catastrophic geohazards[J]. Bulletin of Surveying and Mapping,2020(10):6 − 11. (in Chinese with English abstract) [29] 张允涛. 免相控无人机航摄系统在棚户区改造1∶500地形图中的应用[J]. 测绘通报,2019(增刊 1):295 − 297. [ZHANG Yuntao. Application of phase-free UAV aerial photography system in 1∶500 topographic map of shantytown reconstruction[J]. Bulletin of Surveying and Mapping,2019(Sup 1):295 − 297. (in Chinese with English abstract) [30] 石鼎. 免像控无人机在地形图测绘中的应用[J]. 北京测绘,2020,34(2):163 − 166. [SHI Ding. Application of UAV image free control in topographic map surveying[J]. Beijing Surveying and Mapping,2020,34(2):163 − 166. (in Chinese with English abstract) [31] 朱晓康, 魏景帅. 1∶500免像控无人机航测技术应用研究[J]. 地理空间信息,2019,17(2):22 − 26. [ZHU Xiaokang, WEI Jingshuai. Key technologies and application research on 1∶500 unmanned aerial vehicle photography without image control point[J]. Geospatial Information,2019,17(2):22 − 26. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-4623.2019.02.006 [32] 范秀庆. 无人机免像控技术在地形图测量中的应用研究[J]. 测绘通报,2017(增刊 1):66 − 68. [FAN Xiuqing. Application of UAV image free control technique in topographic map surveying[J]. Bulletin Surveying and Mapping,2017(Sup 1):66 − 68. (in Chinese with English abstract) [33] 相涛, 栾元重, 许章平, 等. 基于Pix4Dmapper的无人机低空摄影测量数据处理[J]. 测绘与空间地理信息,2019,42(3):75 − 78. [XIANG Tao, LUAN Yuanzhong, XU Zhangping, et al. Unmanned aerial vehicle with low-altitude aerial photogrammetry data processing flow and precision analysis based on Pix4Dmapper[J]. Geomatics & Spatial Information Technology,2019,42(3):75 − 78. (in Chinese with English abstract) [34] 戴嵩, 魏冠军, 梁斌. 控制点布设方案对无人机精度测量的影响及其应用[J]. 中国地质灾害与防治学报,2021,32(5):113 − 120. [DAI Song, WEI Guanjun, LIANG Bin. Influence of control point number on UAV low-altitude photogrammetry and its application: A case study in subsidence monitoring of a tailing dam area in northwestern China[J]. The Chinese Journal of Geological Hazard and Control,2021,32(5):113 − 120. (in Chinese with English abstract) [35] 贾虎军, 王立娟, 范冬丽. 无人机载LiDAR和倾斜摄影技术在地质灾害隐患早期识别中的应用[J]. 中国地质灾害与防治学报,2021,32(2):60 − 65. [JIA Hujun, WANG Lijuan, FAN Dongli. The application of UAV LiDAR and tilt photography in the early identification of geo-hazards[J]. The Chinese Journal of Geological Hazard and Control,2021,32(2):60 − 65. (in Chinese with English abstract) [36] 张欢, 巨能攀, 陆渊, 等. 基于无人机的滑坡地形快速重建与稳定性分析[J]. 水文地质工程地质,2021,48(6):171 − 179. [ZHANG Huan, JU Nengpan, LU Yuan, et al. Rapid remodeling of three-dimensional terrain and stability analyses of landslide based on UAV[J]. Hydrogeology & Engineering Geology,2021,48(6):171 − 179. (in Chinese with English abstract) [37] 文宝萍, 王凡. 1965年烂泥沟滑坡前兆、高速远程运动及后期演化特征[J]. 水文地质工程地质,2021,48(6):72 − 80. [WEN Baoping, WANG Fan. Precursors and motion characteristics of the 1965 Lannigou rockslides and the subsequent evolution[J]. Hydrogeology & Engineering Geology,2021,48(6):72 − 80. (in Chinese with English abstract) -
期刊类型引用(13)
1. 张晓丽，曾祥宇，唐流东，刘亮，张阳. 基于信息量法的区域滑坡易发性评价——以宜宾市叙州区为例. 农业灾害研究. 2025(01): 265-267 . 
百度学术
2. 石文君，王宇栋，解晋航，李章杰，梁形形. 基于多种模型对比的寻甸县地质灾害易发性分析. 矿产勘查. 2024(06): 1092-1102 . 百度学术
3. 陆雄轩，朱永浩. 基于信息量模型开展的滑坡崩塌易发性评价. 科技资讯. 2024(12): 196-198 . 百度学术
4. 石文君，王宇栋，李章杰，解晋航. 基于I-LR模型耦合的斜坡类地质灾害易发性评价. 防灾减灾学报. 2024(03): 31-37 . 百度学术
5. 陈绰裕，黄强兵，解庆禹，李伦. 泾阳——渭南活动断裂带场地形变与地层应力特征数值分析. 水文地质工程地质. 2024(05): 136-149 . 百度学术
6. 余杰，丁艳伟. 者么箐泥石流沟发育特征及运动参数研究. 科技与创新. 2024(23): 76-79 . 百度学术
7. 丁思蒙，韩尚鹏，田仕雄，刘鑫. 基于GIS技术的乃东区斜坡类地质灾害易发性评价. 防灾减灾学报. 2024(04): 13-20 . 百度学术
8. 曾斌，吕权儒，寇磊，艾东，许汇源，袁晶晶. 基于Logistic回归和随机森林的清江流域长阳库岸段堆积层滑坡易发性评价. 中国地质灾害与防治学报. 2023(04): 105-113 . 本站查看
9. 胡杨，张紫昭，林世河. 基于证据权与逻辑回归耦合的新疆伊犁河谷地区滑坡易发性评价. 工程地质学报. 2023(04): 1350-1363 . 百度学术
10. 高波，张佳佳，陈龙，田尤，刘建康. 藏东红层地区断裂对泥石流物源的控制作用——以西藏贡觉县哇曲中游流域为例. 中国地质灾害与防治学报. 2023(05): 20-31 . 本站查看
11. 谭彬，刘国栋，余安亮，林平. 基于AHP-INF耦合模型的察雅县地质灾害易发性评价. 地理空间信息. 2023(11): 15-19 . 百度学术
12. 王志民，罗刚，王媛，胡卸文，陈仕阔. 切割斜坡断层的几何形态对斜坡地震响应影响研究. 水文地质工程地质. 2023(06): 147-157 . 百度学术
13. 铁永波，孙强，徐勇，张勇，魏云杰，杨秀元，张泰丽，谭建民. 南方山地丘陵区典型地质灾害成因机制与风险评价. 中国地质调查. 2022(04): 1-9 . 百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 414
- HTML全文浏览量: 195
- PDF下载量: 267
- 被引次数: 14
邮件订阅
RSS