Analysis on the effect of freeze-thaw on landslide prevention projects in seasonal frozen soil area
-
摘要: 川藏铁路规划线路穿越青藏高原东缘高寒带岛状山地冻土区，冻融环境相对复杂。文章通过对川西地区及藏东昌都地区G317、G318沿线滑坡防治工程25处典型案例进行调查，总结冻融条件下滑坡防治工程的破坏形式及受力机理，探讨了水平冻胀力、法向冻胀力的计算方法及其对挡墙、桩间板、排水沟等的影响，计算结果表明水平冻胀力远大于主动土压力，认为开放式冻胀较封闭式冻胀对防治工程的影响更大；阐明了冻胀融沉对锚索工程的影响，提出使用基于蓄能消能系统的锚索结构；探讨了季节性冻土区滑坡防治材料抗冻等级选择方法及减缓冻融效应的措施及建议，为川藏铁路沿线滑坡治理工程及拟建施工便道的边坡支护提供参考。Abstract: The planned route of the Sichuan-Tibet railway passes through the frozen soil on the eastern edge of the Qinghai-Tibet Plateau. The freezing-thawing environment is relatively complicated. More than 25 typical cases of landslide control projects along G317 and G318 in western Sichuan and Qamdo areas were investigated. The failure modes and force mechanisms of the projects under freeze-thaw conditions were summarized. The calculation method of horizontal frost heave force was discussed. The effects of horizontal frost heave retaining walls and drainage ditches were explored. The calculation results show that the horizontal frost heave force is much larger than the active earth pressure. The effects of frost heave and thawing on anchor cable engineering are clarified. It is considered that the effect of open frost heave on control engineering is greater than that of closed frost heave. Methods of material antifreeze grade and measures to reduce freeze-thaw effect were proposed. It can be used as a reference for landslide management projects along the Sichuan-Tibet railway.
-
Keywords:
- freeze-thaw /
- Sichuan-Tibet railway /
- seasonal frozen region /
- landslide /
- anchor
-
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三个时段滑坡分布的驱动因子，结果显示，降雨和“软-硬”互层结构岩组与其它因子的交互作用是影响武隆区滑坡分布的主要影响因素；同时降雨、地质因子和地形地貌等因子对滑坡空间分布的解释力均有不同程度降低，而人类活动因子的解释力不断增强，已逐渐成为影响滑坡分布的关键影响因素之一，是未来减灾防灾重点关注的对象。
-
![]()
图 2 川藏铁路沿线季节性冻土分布图
Figure 2. Seasonal frozen soil distribution along the Sichuan-Tibet railway
![]()
图 3 滑坡防治工程冻融破坏方式及影响因素
Figure 3. Freeze damage and influencing factors of landslide control projects
![]()
图 9 双向调节预应力锚索结构图
1—碟簧耗能系统；2—注浆管；3—保护套筒；4—退位卡套；5—安全卡套；6—封底板；7—钢绞线； 8—锚具；9—架线环；10—导向帽；11—承力转换板；12—自由段； 13—锚固段
Figure 9. Structure of bi-directional adjusting prestressed anchor cable
-
[1] 薛翊国, 孔凡猛, 杨为民, 等. 川藏铁路沿线主要不良地质条件与工程地质问题[J]. 岩石力学与工程学报,2020,39(3):445 − 468. [XUE Yiguo, KONG Fanmeng, YANG Weimin, et al. Main unfavorable geological conditions and engineering geological problems along Sichuan-Tibet railway[J]. Chinese Journal of Rock Mechanics and Engineering,2020,39(3):445 − 468. (in Chinese with English abstract) [2] 梁波, 王家东, 严松宏, 等. 多年冻土地区L型挡土墙土压力(冻胀力)的分析与试验[J]. 冰川冻土,2002,24(5):628 − 633. [LIANG Bo, WANG Jiadong, YAN Songhong, et al. Experiment and analysis of the earth pressure (frost heaving forces) on L-type retaining wall in permafrost regions[J]. Journal of Glaciolgy and Geocryology,2002,24(5):628 − 633. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-0240.2002.05.027 [3] 梁波, 王家东, 葛建军, 等. 青藏铁路L型挡土墙的土压力实测与分析[J]. 岩土工程学报,2004,26(5):627 − 631. [LIANG Bo, WANG Jiadong, GE Jianjun, et al. Testing and analysis of earth pressure about L type retaining wall in Qinghai-Tibet railway[J]. Chinese Journal of Geotechnical Engineering,2004,26(5):627 − 631. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-4548.2004.05.011 [4] 汪恩良, 钟华, 孙景路, 等. 加筋土挡墙冻融试验研究[J]. 岩土工程学报,2010,32(2):265 − 270. [WANG Enliang, ZHONG Hua, SUN Jinglu, et al. Experimental study on reinforced retaining walls suffering freeze-thaw cycling[J]. Chinese Journal of Geotechnical Engineering,2010,32(2):265 − 270. (in Chinese with English abstract) [5] 赵坚, 那文杰, 曹顺星. 季节冻土区挡土墙抗冻结构设计方法[J]. 黑龙江交通科技,2001,24(5):14 − 16. [ZHAO Jian, NA Wenjie, CAO Shunxing. Design of antifreezing structure for retaining walls in seasonal frost regions[J]. Communicatios Science and Technology Heilongjiang,2001,24(5):14 − 16. (in Chinese with English abstract) DOI: 10.3969/j.issn.1008-3383.2001.05.008 [6] ZHANG Y, DONG J H, DONG X G, et al. Analysis of freezing and thawing of slope improved by soil nailing structure in seasonal frozen soil region[J]. Rock and Soil Mechanics,2017,38(2):574 − 582.
[7] 赵延凯, 彭然, 周加伟, 等. 季节性冻土深基坑支护技术研究[J]. 建筑结构,2011,41(7):119 − 122. [ZHAO Yankai, PENG Ran, ZHOU Jiawei, et al. Investigation on technology of retaining and protection for seasonal frozen soil deep foundation pit[J]. Building Structure,2011,41(7):119 − 122. (in Chinese with English abstract) [8] 郭红仙, 宋二祥, 陈肇元. 季节冻土对土钉支护的影响[J]. 工程勘察,2006,34(2):1 − 6. [GUO Hongxian, SONG Erxiang, CHEN Zhaoyuan. Soil nailing under the influence of the seasonal frozen soil[J]. Journal of Geotechnical Investigation & Surveying,2006,34(2):1 − 6. (in Chinese with English abstract) [9] STOCKER M F, RIEDINGER G. The bearing behavior of nailed retaining structures[C]// Design and Performance of Earth Retaining Structures: Proceedings of a Conference Sponsored by the Geotechnical Engineering Division of the American Society of Civil Engineers. New York, June 18-21, 1990, edited by Philip C. Lambe and Lawrence A. Hansen: 612-628.
[10] KINGSBURY D W, SANDFORD T C, HUMPHREY D N. Soil nail forces caused by frost[J]. Transportation Research Record:Journal of the Transportation Research Board,2002,1808(1):38 − 46. DOI: 10.3141/1808-05
[11] LONG N T, T LIVET M, BOUTONNCT, et al. Repair of a reinforced earth wall[C]//Proc. Interactional Conference on Case Histories in Geotechnical Engineering, University of Missouri-Rolla, 1984: 335-346.
[12] 李宗龙. 川西高原川藏交通廊道季节性冻土分布规律及处置措施[D]. 成都: 西南交通大学, 2013. LI Zonglong. Distribution regularity and disposal measures of seasonal frozen soil of Sichuan-Tibet transportation corridor in west Sichuan plateau[D]. Chengdu: Southwest Jiaotong University, 2013. (in Chinese with English abstract)
[13] 中华人民共和国交通运输部. 季节性冻土地区公路设计与施工技术规范: JTG/T D31-06—2017[S]. 北京: 人民交通出版社, 2017. Ministry of Transport of the People's Republic of China. Technical Specifications for Design and Construction of Highway in seasonal Frozen Soil Region: JTG_T D31-06-2017[S]. Beijing: People's Communications Publishing Co. Ltd., 2017. (in Chinese)
[14] 张誉,蒋利学,等. 混凝土结构耐久性概论[M].上海:上海科学技术出版社, 2003:105-117. ZHANG Yu,JIANG Lixue, et al. Durability of concrete structure [M]. Shanghai:Shanghai Science and Technology Press, 2003:105−117. (in Chinese)
[15] 董旭光, 董建华, 何天虎, 等. 多年冻土边坡的水热力耦合分析及软件开发[J]. 中国地质灾害与防治学报,2020,31(4):77 − 84. [DONG Xuguang,DONG Jianhua,HE Tianhu,et al. Coupling analysis and software development of moisture-heat-stressfields of permafrost slope[J]. The Chinese Journal of Geological Hazard and Control,2020,31(4):77 − 84. (in Chinese with English abstract) [16] 张遂, 匡航, 靳占英, 等. 高含水量冻粉黏土应力-应变曲线特性的试验研究[J]. 水文地质工程地质,2020,47(5):116 − 124. [ZHANG Sui, KUANG Hang, JIN Zhanying, et al. An experimental study of the stress-strain characteristics of frozen silty clay with high moisture content[J]. Hydrogeology & Engineering Geology,2020,47(5):116 − 124. (in Chinese with English abstract) [17] 齐子萱, 周金龙, 季彦桢, 等. 西北内陆盆地降水入渗补给季节性变化—以新疆昌吉地下水均衡试验场为例[J]. 水文地质工程地质,2020,47(5):12 − 20. [QI Zixuan, ZHOU Jinlong, JI Yanzhen,et al. Seasonal variationin recharge of infiltration from precipitation for the inland basins of northwestern China: taking the Changji groundwater balance test station in Xinjiang as an example[J]. Hydrogeology & Engineering Geology,2020,47(5):12 − 20. (in Chinese with English abstract) [18] 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局. 水工建筑物抗冰冻设计规范: GB/T 50662—2011[S]. 北京: 中国计划出版社, 2012. Ministry of Housing and Urban-Rural Development of the People's Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. Code for design of hydraulic structures against ice and freezing action: GB/T 50662—2011[S]. Beijing: China Planning Press, 2012. (in Chinese)
[19] 何伟. 西藏自治区国道318线田妥镇滑坡机理及防治技术研究[D]. 重庆: 重庆交通大学, 2012. HE Wei. Research of the Tibet autonomous region national Highway 318 Line of the landslide mechanism and the field prevention and control technology[D]. Chongqing: Chongqing Jiaotong University, 2012. (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)
下载:
计量
- 文章访问数: 565
- HTML全文浏览量: 356
- PDF下载量: 269
- 被引次数: 14
邮件订阅
RSS