Susceptibility assessment of geological hazards based on susceptibility quantitative factors: A case study in Qijiang District, Chongqing City
-
摘要: 綦江区是重庆市地质灾害高发区域,本研究以綦江区为研究区,根据区内地形地貌、地质构造并结合DEM数据,提取了斜坡单元内的斜坡结构、坡度、地灾点密度、工程地质岩组、距水系距离、岩层倾角和地质构造共7项主要致灾因子,进行量化分析并结合地质灾害野外现场核查和修正其各因子权重和分级赋值后,采用层次分析法进行斜坡单元地质灾害易发性评价,探索建立适合綦江区的地质灾害易发性评价体系。Abstract: Qijiang District is a high-risk area for geological hazards in Chongqing. In this study, Qijiang District is selected as the research area. Based on the topography, geomorphology, geological structure, and DEM data of the region, seven main factors of geological hazards within slope units were extracted, including slope structure, slope gradient, density of geological disaster points, engineering geological lithology, distance to the water system, inclination angle of rock strata, and geological structure. These factors were quantitatively analyzed and combined with field inspection of geological hazards to revise the weights and grading of each factor. The analytic hierarchy process was used to evaluate the susceptibility of geological hazards within slope units. Through this process, this study aims to establish a geological disaster susceptibility evaluation system suitable for Qijiang District.
-
0. 引言
近年来,随着国家对天然气管道工程建设的重视力度不断提高,我国的天然气管道工程建设步伐不断加快,目前国家干线管道系统已初具规模,并加速形成“主干互联、区域成网”的全国天然气基础网络[1]。
由于管道线路较长,通常跨越不同地区,沿线地质地貌条件复杂,因此在建设过程中难以避免会穿越不良地质地段,在不稳定因素的影响下极有可能诱发地质灾害,其中滑坡是山区油气管道工程中常见的地质灾害类型[2]。在滑坡作用下,管道容易产生变形甚至破损,导致燃气泄漏,严重威胁沿线民众的正常生活。因此,如何对管道滑坡进行有效的支护并保护管道安全值得关注。
花管微型桩是一种新型压力灌浆式抗滑桩,集中了钢花管桩与微型桩的优点,在20世纪80年代开始在国内外应用并很快在基坑支护、滑坡治理等领域被广泛采用。目前,国内外学者对花管微型桩加固滑坡的承载性能已经开展了部分研究工作,但针对管道—滑坡体系支护性能研究方面还未涉及。潘锋[3]提出了一种注浆钢花管桩加固滑坡的理论计算公式;陈强等[4]通过离心模型试验研究发现注浆钢管桩在控制滑坡土体位移方面的有益作用;Wang等[5]针对钢花管桩微型桩的双注浆技术展开了研究,并通过试验证明了此工艺可以有效提高微型桩承载性能。
螺纹微型桩来源于地基基础工程中使用的螺纹桩,因其利用“螺丝钉比钉子更牢固”的原理[6],使螺纹桩拥有较好的承载性能,但螺纹桩目前只在地基工程中应用较多,还未被作为支挡结构应用在滑坡治理工程中。国内外对其进行的研究多集中在竖向承载性能方面,如Krasiński[7]采用数值模拟方法对螺纹桩承载机理进行了研究,并通过现场载荷试验验证了方法的可靠性;方崇等[8]通过静载试验分析了螺杆桩的竖向承载力传递特性与受力特征;叶阳升等[9]对高速铁路中的螺杆桩复合地基进行了原位测试试验,研究获得了直杆端与螺纹段在不同受力状态下的桩侧摩阻力关系;Malik等[10]进行了螺纹桩与直杆桩承载性能的对比研究,结果证明螺纹桩轴径比为2~4.1的螺旋桩端承载力比类似桩轴径的直杆桩高2~12倍;孟振等[11]通过模型试验手段针对砂土地基中的螺杆桩承载特性展开了研究。
为探究花管微型桩与螺纹微型桩加固管道—滑坡体系下桩体变形特征、桩体两侧土压力的空间分布规律以及在两种不同微型桩加固保护下管道的变形特点及破坏模式,进一步分析对比两种不同微型桩的支护性能,本文以中贵天然气管道K558+700滑坡工程为依托,以千斤顶模拟滑坡水平推力[12],开展7种不同水平加载压力工况下的室内大型推桩模型试验,并对桩顶位移、桩体应变、桩体两侧土压力、天然气管道应变以及管道两侧的土压力进行测试。进一步分析桩体弯矩以及桩体两侧峰值土压力空间分布规律,确定不同微型桩桩体变形特征,并结合管道变形破坏情况,探究两种微型桩在管道滑坡下的承载特性,总结试验成果,为我国管道地段滑坡灾害治理提供理论参考。
1. 工程概况
中贵天然气管道K558+700滑坡治理工程位于甘肃省陇南市成县黄陈镇中湾村。2020年8月12—18日中湾村附近连降暴雨,导致中贵天然气管道K558+700处斜坡发生明显滑动变形,管道滑坡区域现场见图1。共产生2处滑坡及1处滑塌体。其中,中贵天然气管道位于H1滑坡中部,横坡敷设通过,管径1016 mm,管道走向186°,滑坡区管道埋深2.2~4.1 m。H1滑坡纵长约190 m,宽约186 m,滑坡平面面积约2.25×104 m2,滑体平均厚度约8.5 m,滑体总体积约19.13×104 m2,属于中型土质滑坡。该滑坡地表变形强烈,尤其是与管道斜交的乡村水泥道路已完全损毁,受滑坡推挤作用,威胁管道安全运营。
2. 试验设计
2.1 相似设计
本次试验设计以中贵天然气管道K558+700处H1滑坡典型断面为工程原型,参考试验模型箱尺寸,确定试验模型的几何相似比尺为:
$$ {C}_{L}=\frac{{L}_{P}}{{L}_{M}}=30 $$ (1) 式中:CL——几何相似常数;
LP——原型尺寸;
LM——模型尺寸。
根据模型试验中各因素对现象影响的大小,抓住其起主要作用的因素,略去其次要因素的原则。以微型桩几何尺寸(L)、密度(ρ)、重力加速度(g)为主要控制参数,其相似比分别为30∶1、1∶1、1∶1,其余参数相似比可根据Buckinghamπ定理导出,如表1所示。
表 1 相似比设计Table 1. Design of similarity constant物理量 相似比 物理量 相似比 几何尺寸 CL=30 变形模量 $ {{C}}_{{E}}=30 $ 质量密度 $ {{C}}_{{\rho }}=1 $ 摩擦角 ${{C} }_{{\varphi } }=1$ 重度 Cy=1 黏聚力 $ {{C}}_{{C}}=30 $ 应变 $ {{C}}_{{\epsilon }}=1 $ 时间 ${ {C} }_{ {t} }=\sqrt {30}$ 位移 $ {{C}}_{{u}}=30 $ 重力加速度 $ {{C}}_{{g}}=1 $ 2.2 相似材料及配比设计
制作与原型形状完全相同的螺纹桩、花管桩较为困难,因此对其桩体模型进行简化,采用阻氧双色PP-R管材通过注水泥浆模拟桩体,其中花管桩模型为在PP-R管材长度方向上每隔8 cm环绕管道螺旋打孔并采用压力灌注水泥浆模拟;螺纹桩模型为在PP-R管材外径螺旋缠绕塑料软管并灌注水泥浆进行模拟;承台采用硬质木板模拟;天然气管道采用直径20 cm的PVC管模拟。制作完成的模型桩见图2。
滑坡岩土体相似材料较为复杂,在配比过程中不易满足推导得到的相似关系,因此试验选取滑体、滑带及基岩材料时,以最易影响原型滑坡岩土体稳定性的若干参数 (容重、黏聚力、内摩擦角、弹性模量等)为基础,参考其他学者已得到的部分研究成果[13],通过正交配比设计,并经直剪试验和三轴试验检验参数取值的准确合理性,确定使用红粉土、石英砂、水泥、石膏、水的混合物模拟Ⅳ级基岩,具体质量配比为70∶30∶5∶3∶10;滑体采用红粉土、石英砂、水的混合物模拟,具体质量配比为70∶20∶10;滑动带采用石英砂、土、滑石粉、水的混合物模拟,具体质量配比为27∶52∶35∶15。得到的模型材料相关力学参数及其与原型材料的对比见表2。
表 2 模型材料与原型材料相关物理性质参数Table 2. Physical property parameters related to model material and prototype material物理力学参数 重度/(kN·m−3) 内摩擦角/(°) 黏聚力/kPa 弹性模量/MPa 滑体 原型 19.0 30.0 40.0 / 模型 19.0 30.5 1.2 / 滑带 原型 18.5 25.0 30.0 / 模型 18.0 25.0 1.0 / 基岩 原型 27.0 / / 5000 模型 26.7 / / 160 2.3 模型设计
根据相似设计,试验使用的模型箱尺寸为 4 m×2 m×2 m(长×宽×高),箱体两侧由钢板及有机玻璃组成,箱体内填筑土体,从上往下分三层:滑体、滑带及基岩层。箱体内共布置4组微型桩群模型,其中左侧2组为花管桩,右侧2组为螺纹桩;桩后滑体内埋置有天然气管道模型。箱体后侧设置有反力墙,反力墙上安装有液压千斤顶,千斤顶通过作用在承压板上从而给滑坡模型施加水平推力。模型整体示意图如图3。
微型桩群桩模型单桩布置形式为5行×3列,共布置4组微型群桩,以模型箱中线为分界线分为左右两侧,左侧2组为花管桩,右侧2组为螺纹桩;单桩桩体模型长160 cm,直径2.5 cm。微型桩群桩埋深160 cm,从桩顶到桩底依次穿过滑体、滑带与基岩层,其中滑体层厚75 cm,滑体坡面角40°,滑带厚度5 cm,基岩层厚100 cm;天然气管道布置在桩前滑体中,距离桩体20 cm,埋深40 cm。桩顶采用承台方式将群桩连接为整体,单桩中心间距为7 cm,承台长38 cm,宽24 cm,厚2 cm,如图4所示。
2.4 传感器布设位置及监测设备
试验数据采用DH3816N静态应变测试系统采集,该系统具有60个采集通道,可同时采集应变、土压力等测试数据。
2.4.1 百分表与应变片布设位置
选取模型箱体中间两组微型桩群桩作为研究对象,在桩顶承台设置百分表测试桩顶位移,并在群桩中选取不同位置的典型单桩作为应变测试桩粘贴应变片,百分表设置位置及测试单桩的位置如图5所示。
测试单桩为2#花管桩及3#螺纹桩群桩中的第二列单桩,并对其进行编号,从桩后到桩前依次为1、2、3号测试单桩,每一测试单桩在桩体两侧沿桩体深度间隔一定距离粘贴应变片,每根测试桩体上共粘贴2×6=12个应变片,所有测试桩体共布设12×6=72个应变片,具体布设位置如图6所示(以1号测试单桩为例,其余测试桩相同)。
天然气管道山侧与河侧对称粘贴3组应变片,一组位于花管桩侧中部,一组位于螺纹桩侧中部,以监测在两种不同微型桩群桩支护作用下天然气管道的变形并作为对比,进一步优化天然气管道滑坡的支护方案设计;另一组粘贴在天然气管道接口处,以监测管道接口处的变形,进一步研究在滑坡作用下天然气管道接口处的处理方法。天然气管道应变片布设位置见图7。
2.4.2 土压力布设位置
为了研究花管微型桩与螺纹微型桩支护下桩周土压力空间分布规律,在2#群桩与3#群桩前后沿桩体深度间隔一定距离对称布置土压力盒,布设位置如图8,共布设6×2×2=24个土压力盒,以监测在滑坡水平推力作用下两种微型桩群桩前后不同深度位置土压力的变化,以进一步分析两种不同桩体变形特征。
在天然气管道前后两侧对称布设土压力盒,布设位置分别位于花管桩支护侧、螺纹桩支护侧及管道接口附近,共布设3×2=6个土压力盒,以监测天然气管道附近土体的应力变化。具体布设位置如图9。
2.5 加载工况设计
试验采用油压千斤顶施加水平荷载模拟滑坡推力,加载方法为慢速维持加载法分级加载[14],初始加载压力为1 MPa,每级加载压力增加0.5 MPa,试验一共设置7组加载工况。每级加载后观察桩顶百分表读数稳定后进行下级加载,具体加载工况见表3。
表 3 各级加载工况Q值Table 3. Q values of loading conditions at all levels工况 加载压力/MPa 工况 加载压力/MPa 1 1.0 5 3.0 2 1.5 6 3.5 3 2.0 7 4.0 4 2.5 − − 3. 试验结果分析
3.1 桩顶水平位移
由图10荷载下桩顶水平位移数据分析可知,在各级水平荷载作用下,螺纹桩侧3#与4#群桩顶水平位移明显大于花管桩侧1#与2#群桩,在4 MPa加载压力作用下螺纹桩侧位移量为花管桩侧位移量的两倍之多,说明在相同水平推力荷载下花管桩抗滑承载力更强;在2 MPa以下水平荷载作用下桩顶位移变化微弱,2~3 MPa时桩顶位移开始有小幅度增长,当荷载达到3 MPa时,桩顶位移变化量出现转折,增长幅度明显加大。根据桩顶位移容许值为30~50 mm,以1.67 mm(相当于原型的50 mm)为作为试验模型桩顶位移容许值,螺纹桩侧桩顶位移量在水平荷载为2.5 MPa达到容许值,花管桩侧桩顶位移量在水平荷载为3 MPa达到容许值,表明在本次试验条件下花管微型桩的加固效果优于螺纹微型桩。
3.2 微型桩群桩两侧土压力分析
3.2.1 桩体两侧土压力变化时程特性
根据图11可以看出,除个别测点外,在加载时段微型桩山侧土压力变化呈阶梯状上升,反映了随着加压荷载的增大,微型桩靠山侧土体应力在不断增加,并且有明显分级现象。
对比两侧不同桩体类型的土压力变化可知,花管桩侧不同埋深的土压力值普遍高于螺纹桩侧,花管桩侧在TH2测点达到最大值为46.09 kPa,相比螺纹侧在TL2测点达到的最大值20.54 kPa,花管桩侧土体最大应力高出螺纹桩侧两倍以上,结合螺纹桩桩顶最大位移与花管桩桩顶位移比值为2.4,说明花管桩侧土体受挤压程度大,分析为土拱效应[15]导致,在花管微型桩支护作用下桩后土体应力得不到有效释放,进一步表明花管桩的抗弯刚度高于螺纹桩。
由图11不同测点深度的土压力变化曲线可以看出,螺纹微型桩与花管微型桩的中上部位置土压力值最大,桩顶次之,桩底一般较小。
由图12,河侧螺纹桩测点4因土压力盒出现故障导致采集数据失真,故不作为分析对象。由图12可知,桩体河侧土压力变化趋势与山侧几乎保持一致,在加载阶段都有明显的上升趋势且各级加载土压力有明显的上升幅度加大的现象,卸载阶段开始下降并最终稳定在零值附近。
两侧桩体都在测点3(滑带附近位置)的最大级加载阶段达到峰值土压力,且螺纹桩侧峰值土压力(38.81 kPa)大于花管桩侧峰值土压力(32.57 kPa),高出比值约120%,对比山侧峰值土压力,花管桩河侧峰值土压力比山侧减少41%,螺纹桩河侧峰值土压力增加89%,由于抗滑桩承担并抵消了一部分滑坡推力,桩体河侧的土压力明显小于山侧,但花管桩两侧土压力差值明显大于螺纹桩两侧土压力差值,表明花管桩在承担滑坡推力作用方面发挥了良好的效果。
3.2.2 桩体两侧峰值土压力分析
提取桩体两侧各测点在各级加压荷载下土压力的峰值土压力进行单独分析,如图13绘制沿桩身深度不同位置随荷载改变的峰值土压力变化曲线分布图,峰值土压力基本随加压荷载的增大而增大,花管桩侧在4 MPa荷载下0.4 m桩深位置达到最大值46.09 kPa,螺纹桩侧在4 MPa荷载下0.4 m桩深位置达到最大值20.54 kPa。
花管桩侧峰值土压力沿桩深分布图大体呈“S”曲线形,且在各级荷载下分布规律保持一致,桩体中上部(位于滑带与桩顶中部)土压力值最大,在滑带附近土压力值有所减小,又在桩体中下部(位于滑带与桩底中部)有小幅度增大,往桩底位置开始减小。螺纹桩侧峰值土压力沿桩深分布图大体呈双“S”曲线形,同样在桩体中上部位置土压力值最大,滑带附近减小,桩体中下部土压力值最小,桩底土压力值开始增大。总的来说,花管桩侧土压力分布规律为桩中>桩顶>桩底,螺纹桩侧土压力分布规律为桩中>桩底>桩顶。
滑面以上的土压力由推力荷载的传递引起,滑面以下的土压力主要由微型桩群的变形引起[16]。根据上述土压力分布规律,由此说明,花管桩在滑面以上桩中上部位置所受的推力荷载最大,滑面以下桩体中下部位置变形最大;螺纹桩在滑面以上桩中上部位置所受的推力荷载最大,滑面以下桩底位置变形最大。
根据图14可知,花管桩侧在4 MPa荷载下0.65 m桩深位置达到最大值32.57 kPa,螺纹桩侧在4 MPa荷载下0.65 m桩深位置达到最大值38.81 kPa。
花管桩河侧峰值土压力分布形式与山侧分布形式大体呈现相似,土压力大小分布规律仍为桩中>桩顶>桩底,但土压力最大值出现位置下移到滑面附近,说明在花管微型桩支护作用下推力荷载在土体的分布形式发生了改变,由于花管桩刚度大于桩周土体,土体应力传递到桩体上产生应力集中,桩后土体应力重分布,最大水平应力作用点下移。螺纹桩侧河侧峰值土压力相比山侧分布形式产生了较大差异,土压力大小分布规律变为桩中>桩顶>桩底,桩底的土压力相比山侧有很大幅度减小。与花管桩侧类似,由于桩前土体应力重分布,土压力最大值出现位置下移到滑面附近且比山侧土压力最大值高出89%,说明螺纹微型桩并没有有效承担推力荷载,变形较大,在推力荷载作用桩前土体应力仍较大。
3.3 桩身弯矩分析
选取花管微型桩1号测试单桩与螺纹微型桩1号测试单桩作为研究对象,根据桩身两侧应变数据可由以下公式计算测点弯矩值[17]:
$$ M = \frac{{EI({\varepsilon _1} - {\varepsilon _2})}}{d} $$ (2) 式中:M——测试截面的弯矩/(kN·m);
E——桩身材料的弹性模量/MPa;
I——测试界面的惯性矩/m4;
ε1——微型桩后的桩身应变;
ε2——微型桩前的桩身应变;
d——桩身直径/m。
图15—16为两种桩体各截面测点在分级加压荷载下的弯矩分布图,图中弯矩正值表示桩后侧受拉。
从图15可以看出,在1~3.5 MPa加压荷载下,弯矩分布特征几乎保持一致,花管桩桩身弯矩沿深度方向呈“M”形分布,桩顶与桩底弯矩值较小,桩身中上部与中下部出现较大的负弯矩,桩身离模拟滑面以上5 cm位置处产生最大正弯矩;在4 MPa加载压力作用下,弯矩分布曲线有所偏移,桩身中上部位置负弯矩极具增大,最大正弯矩位置下移到桩中位置(模拟滑面处)且桩身正弯矩长度范围增大到0.5 m,滑面以下负弯矩最大值减小,且在桩底出现正弯矩。说明在1~3.5 MPa推力荷载下,花管桩桩体变形特征一致,桩后桩身中上部与中下部位置受压变形,桩后桩身离模拟滑面0~10 cm段受拉变形;在4 MPa推力荷载下,桩身弯曲变形加剧,在桩深0.4 m 处出现最大负弯矩−71.96 kN·m,在桩深0.8 m处(模拟滑面位置)出现最大正弯矩16.47 kN·m。花管桩设计时对桩身变形较大处可开展进一步的优化工作。
由图16可知,螺纹桩桩身弯矩分布沿深度方向呈“S”形,桩体正负弯矩位置在模拟滑面附近大致呈旋转对称分布,滑面以上大部分区段为负弯矩,滑面以下为正弯矩;螺纹桩在桩深0.4处出现最大负弯矩−64.44 kN·m,在桩深0.9m处(模拟滑面以下10 cm)出现最大正弯矩65.85 kN·m,相比花管桩,最大负弯矩相差不大,最大正弯矩高出近4倍,表明在相同推力荷载工况下,螺纹微型桩变形程度大于花管微型桩。
3.4 管道受力变形特征分析
3.4.1 管道两侧土压力变化规律及分布特性
为了研究天然气管道在两种不同微型桩支护作用下的变形特征,通过埋设土压力盒监测管道两侧土体应力变化,图17为管道山侧与河侧各测点的土压力变化时程曲线。其中TG1/TG1'测点位于花管桩支护侧,TG3/TG3'测点位于螺纹桩支护侧,TG2/TG2'测点位于管道接口位置两侧。
根据图17可以看出,管道山侧与河侧各测点土压力在加载阶段都有明显的随时间增大的趋势。山侧各测点变化趋势大致相同,在初期加载土压力值增长表现平稳,而在后期加载阶段土压力出现跳跃式增长,并在最后加载阶段达到峰值,位于螺纹桩支护侧的TG3测点峰值土压力达22.99 kPa,为花管桩支护侧的1.73倍;结合图18各级荷载下管道山侧峰值土压力分布曲线可知,管道山侧峰值土压力随荷载的增大而增大,并且螺纹桩支护侧峰值土压力整体大于花管桩支护侧及管中接口位置。说明花管桩在推力荷载下承载力表现优于螺纹桩,减小了传递到管道的滑坡推力,而螺纹桩在较大推力荷载下由于抗弯性能不足,变形严重,破坏后不能有效承担并抵消滑坡推力,传递到桩前管道的应力较大,从而导致管道易产生变形破坏。
由图17,管道河侧土压力变化趋势与山侧存在较大差异,在加载阶段整体呈波动性增长,测点最大峰值土压力在螺纹桩支护侧为12.29 kPa,最小峰值土压力在花管桩支护侧为1.54 kPa,比值7.98∶1。结合图19管道河侧峰值土压力分布曲线,在分级推力荷载作用下,峰值土压力规律表现为:花管桩支护侧<管中接口位置<螺纹桩支护侧。
3.4.2 管道弯矩分析
由图20分级荷载管道弯矩分布图可以看出,在较小推力荷载作用下(1~3 MPa),管道接口附近弯矩值整体较小,花管桩支护侧呈现负弯矩,螺纹桩支护侧呈现正弯矩,绝对值大小相差不大,说明在较小推力荷载下,管道变形程度较小;在较大推力荷载作用下(3.5~4 MPa),管道接口处弯矩急剧增大,最大值达−317.75 kN·m,说明在较大推力荷载下管道接口位置变形严重甚至破坏,在设计滑坡区管道工程时应注意管道接口位置的加固;另一方面,对比花管桩支护侧与螺纹桩支护侧管道弯矩可知,弯矩分布曲线明显向螺纹桩侧偏移,并整体大于花管桩支护侧,说明在本试验条件下花管微型桩支护作用下天然气管道的变形程度低于螺纹微型桩支护,在滑坡推力作用下 ,花管微型桩对管道的保护效益更突出,花管桩更适用于作为管道滑坡区域的支挡结构。
4. 结论
(1)花管微型桩山侧及河侧峰值土压力沿桩深分布形式基本相似,大体呈“S”曲线形,桩后土体土拱效应明显,且在各级荷载下分布形式大致保持一致,总体来说花管桩侧土压力分布规律为桩中>桩顶>桩底,滑带附近的桩体周围土压力较大,在抗滑桩设计工作中应重点考虑优化。
(2)螺纹桩山侧峰值土压力沿桩深分布图大体呈双“S”曲线形,河侧峰值土压力相比山侧分布形式产生了较大差异,桩底的土压力相比山侧减小很大幅度;随外部荷载的增加桩周土压力增加幅度较大,表明螺纹微型桩在横向承载性能方面有所欠缺。
(3)花管桩桩身弯矩沿深度方向呈“M”形分布,桩身离模拟滑面以上5 cm位置处产生最大正弯矩;螺纹桩桩身弯矩分布沿深度方向呈“S”形,桩体正负弯矩位置在模拟滑面附近大致呈旋转对称分布,滑面以上大部分区段为负弯矩,滑面以下为正弯矩;在相同推力荷载工况下,螺纹微型桩变形程度大于花管微型桩。
(4)在滑坡作用下花管微型桩可以有效减小传递到管道的坡体应力,在一定程度上预防管道受力破坏;而螺纹桩在较大推力荷载下抗弯性能不足,变形严重,破坏后不能有效承担滑坡推力,传递到桩前管道的应力较大,从而导致管道变形程度更为强烈;在本试验条件下,花管微型桩对管道的保护效益突出,更适用于作为管道-滑坡区域的支挡结构。
-
表 1 斜坡结构因子易发性量化指标
Table 1 Quantitative index of susceptibility factor of slope structure
地灾
类型斜坡
结构分级面积比 地灾面积
占比易发性指标
滑坡
数据近水平 0.06 0.01 −1.80 顺向坡 0.17 0.56 1.20 斜向坡 0.31 0.17 −0.60 横向坡 0.30 0.20 −0.40 逆向坡 0.16 0.06 −1.00 崩塌
数据近水平 0.06 0.004 −2.78 顺向坡 0.17 0.105 −0.47 斜向坡 0.31 0.340 0.09 横向坡 0.30 0.202 −0.39 逆向坡 0.16 0.349 0.76 表 2 孕灾体点密度因子易发性量化指标
Table 2 Quantitative index of susceptibility factor of the point density of the disaster-pregnant body
地灾
类型斜坡
结构分级面积比 地灾面积
占比易发性指标
滑坡
数据<0.2 0.47 0.14 −1.33 [0.2, 0.4) 0.34 0.30 −0.24 [0.4, 0.6) 0.11 0.24 0.61 [0.6, 0.8) 0.06 0.24 1.35 [0.8, 1.0] 0.02 0.08 1.47 崩塌
数据<0.2 0.47 0.49 0.30 [0.2, 0.4) 0.34 0.96 0.59 [0.4, 0.6) 0.11 0.05 0.03 [0.6, 0.8) 0.06 0.09 0.05 [0.8, 1.0] 0.02 0.05 0.03 表 3 工程地质岩组因子易发性量化指标
Table 3 Quantitative index of susceptibility factor of engineering geological lithology
地灾
类型地质
岩组主要岩性 分级面
积比地灾面积
占比易发性
指标滑坡
数据Ⅰ1 灰岩、白云岩 0.02 0.014 −0.49 Ⅰ2 含泥质灰岩 0.01 0.004 −1.22 Ⅱ1 泥岩、砂岩 0.17 0.175 0.08 Ⅱ2 砂岩夹页岩 0.63 0.761 0.20 Ⅱ3 厚层状砂岩 0.03 0.013 −0.80 Ⅲ1 泥灰岩夹粉砂岩、页岩 0.02 0.030 0.52 Ⅲ2 灰岩夹页岩 0.10 0.001 −4.55 Ⅲ3 页岩与灰岩 0.02 0.002 −2.40 崩塌
数据Ⅰ1 灰岩、白云岩 0.02 0.187 2.07 Ⅰ2 含泥质灰岩 0.01 0.007 −0.69 Ⅱ1 泥岩、砂岩 0.17 0.040 −1.40 Ⅱ2 砂岩夹页岩 0.63 0.491 −0.24 Ⅱ3 厚层状砂岩 0.03 / / Ⅲ1 泥灰岩夹粉砂岩、页岩 0.02 0.002 −1.92 Ⅲ2 灰岩夹页岩 0.10 0.273 1.05 Ⅲ3 页岩与灰岩 0.02 / / 表 4 坡度因子易发性量化指标
Table 4 Quantitative index of susceptibility factor of slope gradient
地灾类型 坡度/(°) 分级面积比 地灾面积占比 易发性指标 滑坡
数据<10 0.2199 0.161 −0.31 [10, 20) 0.3969 0.482 0.20 [20, 30) 0.2530 0.298 0.16 [30, 40) 0.0949 0.051 −0.61 [40, 50) 0.0264 0.007 −1.38 [50, 60) 0.0067 0.001 −1.62 [60, 70) 0.0017 / / >70 0.0005 / / 崩塌
数据<10 0.2199 0.015 −2.65 [10, 20) 0.3969 0.090 −1.48 [20, 30) 0.2530 0.203 −0.22 [30, 40) 0.0949 0.273 1.07 [40, 50) 0.0264 0.231 2.19 [50, 60) 0.0067 0.129 3.02 [60, 70) 0.0017 0.057 3.78 >70 0.0005 0.002 3.11 表 5 距水系距离因子易发性量化指标
Table 5 Quantitative index of susceptibility factor of distance from water system
地灾
类型距水系
距离/km分级面积比
地灾面积占比
易发性
指标滑坡
数据<0.2 0.33 0.612 0.64 [0.2, 0.4) 0.26 0.176 −0.39 [0.4, 0.6) 0.18 0.124 −0.38 [0.6, 0.8) 0.12 0.068 −0.57 [0.8, 1) 0.06 0.019 −1.16 ≥1.0 0.05 0.001 −4.35 表 6 岩层倾角因子易发性量化指标
Table 6 Quantitative index of susceptibility factor of inclination angle of rock strata
地灾
类型岩层
倾角/(°)分区面
积比地灾面积
占比易发性
指标滑坡
数据<10 0.191 0.225 0.18 [10, 20) 0.455 0.409 −0.09 [20, 30) 0.216 0.343 0.48 [30, 40) 0.081 0.015 −1.70 [40, 50) 0.038 0.005 −2.06 [50, 60) 0.015 0.002 −1.86 [60, 70) 0.003 / / >70 0.001 0.001 −0.47 崩塌
数据<10 0.191 0.4214 0.81 [10, 20) 0.455 0.1483 −1.10 [20, 30) 0.216 0.3421 0.48 [30, 40) 0.081 0.0666 −0.18 [40, 50) 0.038 0.0210 −0.59 [50, 60) 0.015 / / [60, 70) 0.003 0.0006 −1.83 >70 0.001 / / 表 7 地质构造因子易发性量化指标
Table 7 Quantitative index of susceptibility factor of geological structure
地灾
类型距断层
距离/km分区面
积比地灾面积
占比易发性
指标滑坡
数据(0, 0.5] 0.084 0.139 0.51 (0.5, 1] 0.079 0.079 0.01 (1, 1.5] 0.076 0.041 0.63 (1.5, 2] 0.075 0.068 0.08 (2, 2.5] 0.073 0.056 0.26 (2.5, 3] 0.068 0.085 0.23 (3, 3.5] 0.063 0.023 −0.99 (3.5, 4] 0.054 0.015 −1.30 (4, 4.5] 0.045 0.017 −0.95 (4.5, 5] 0.038 0.018 −1.77 >5 0.345 0.459 −0.29 崩塌
数据(0, 0.5] 0.084 0.072 0.09 (0.5, 1] 0.079 0.133 0.76 (1, 1.5] 0.076 0.069 0.14 (1.5, 2] 0.075 0.061 0.04 (2, 2.5] 0.073 0.003 −3.00 (2.5, 3] 0.068 0.228 −1.86 (3, 3.5] 0.063 0.134 −1.99 (3.5, 4] 0.054 0.088 −1.71 (4, 4.5] 0.045 0.105 −1.09 (4.5, 5] 0.038 0.004 −1.90 >5 0.345 0.103 −0.97 表 8 滑坡影响因素分级层次分析法权重取值表
Table 8 The weighted value table of the hierarchical analysis for landslide influencing factors
滑坡评价
因子权重 因子
分级易发性指标
易发性指标归一化
斜坡
结构0.20 近水平 −1.80 0.00 顺向坡 1.20 0.47 斜向坡 −0.60 0.19 横向坡 −0.40 0.22 逆向坡 −1.00 0.13 地灾点密度 0.17 <0.2 −1.33 0.00 [0.2, 0.4) −0.24 0.11 [0.4, 0.6) 0.61 0.22 [0.6, 0.8) 1.35 0.33 [0.8, 1.0) 1.47 0.34 工程地质岩组 0.15 Ⅰ1 −0.49 0.15 Ⅰ2 −1.22 0.12 Ⅱ1 0.08 0.17 Ⅱ2 0.20 0.17 Ⅱ3 −0.80 0.14 Ⅲ1 0.52 0.18 Ⅲ2 −4.55 0.00 Ⅲ3 −2.40 0.08 坡度
/(°)0.20 (0, 10] −0.31 0.21 (10, 20] 0.20 0.30 (20, 30] 0.16 0.29 (30, 40] −0.61 0.16 (40, 50] −1.38 0.04 (50, 60] −1.62 0.00 (60, 70] / / >70 / / 距水系距离
/km0.15 <0.2 0.64 0.25 [0.2, 0.4) −0.39 0.20 [0.4. 0.6) −0.38 0.20 [0.6, 0.8) −0.57 0.19 [0.8, 1.0) −1.16 0.16 ≥1.0 −4.35 0.00 岩层
倾角
/(°)0.08 (0, 10] 0.18 0.25 (10, 20] −0.09 0.22 (20, 30] 0.48 0.29 (30, 40] −1.70 0.04 (40, 50] −2.06 0.00 (50, 60] −1.86 0.02 (60, 70] / / >70 −0.47 0.18 距断层距离
/km0.05 (0, 0.5] 0.51 0.17 (0.5, 1] 0.01 0.13 (1, 1.5] −0.63 0.06 (1.5, 2] −0.08 0.12 (2, 2.5] −0.26 0.10 (2.5, 3] 0.23 0.15 (3, 3.5] −0.99 0.03 (3.5, 4] −1.30 0.00 (4, 4.5] −0.95 0.03 (4.5, 5] −0.77 0.05 >5 0.29 0.15 表 9 崩塌影响因素分级层次分析法权重取值表
Table 9 The weighted value table of the hierarchical analysis for the collapse influencing factors
崩塌评价
因子权重
因子
分级易发性指标
易发性指标归一化
斜坡结构 0.20 近水平 −2.78 0.00 顺向坡 −0.47 0.21 斜向坡 0.09 0.26 横向坡 −0.39 0.22 逆向坡 0.76 0.32 地灾点密度 0.15 <0.2 0.30 0.23 [0.2, 0.4) 0.59 0.45 [0.4, 0.6) 0.03 0.02 [0.6, 0.8) 0.05 0.04 [0.8, 1.0] 0.03 0.02 >1.0 0.30 0.23 工程地质
岩组0.20 Ⅰ1 2.07 0.38 Ⅰ2 −0.69 0.12 Ⅱ1 −1.40 0.05 Ⅱ2 −0.24 0.16 Ⅱ3 − − Ⅲ1 −1.92 0.00 Ⅲ2 1.05 0.29 Ⅲ3 − − 坡度
/(°)0.30 (0, 10] −2.65 0.00 (10, 20] −1.48 0.04 (20, 30] −0.22 0.08 (30, 40] 1.07 0.12 (40, 50] 2.19 0.16 (50, 60] 3.02 0.19 (60, 70] 3.78 0.21 >70 3.11 0.19 岩层
倾角
/(°)0.05 (0, 10] 0.81 0.31 (10, 20] −1.10 0.09 (20, 30] 0.48 0.27 (30, 40] −0.18 0.19 (40, 50] −0.59 0.14 (50, 60] − − (60, 70] −1.83 0 >70 − − 距断层距离
/km0.10 (0, 0.5] 0.09 0.10 (0.5, 1] 0.76 0.12 (1, 1.5] 0.14 0.10 (1.5, 2] 0.04 0.10 (2, 2.5] −3.00 0.00 (2.5, 3] −0.86 0.07 (3, 3.5] 0.99 0.13 (3.5, 4] 0.71 0.12 (4, 4.5] 1.09 0.14 (4.5, 5] −1.90 0.04 >5 −0.97 0.07 表 10 地质灾害易发性分区统计表
Table 10 Statistical table of geological disaster susceptibility zones
易发性等级 高易发区 中易发区 低易发区 非易发区 斜坡面积/km2 519.99 1021.29 598.93 47.27 占比/% 23.77 46.69 27.38 2.16 表 11 近年来灾(险)情易发性统计表
Table 11 Statistical table of geological disasters (hazards) susceptibility in recent years
易发分区 地灾点数量/个 占比/% 高易发区 30 58.82 中易发区 17 33.33 低易发区 4 7.85 非易发区 0 0 -
[1] 张倬元, 王士天, 王兰生. 工程地质分析原理[M]. 2版. 北京: 地质出版社, 1994 ZHANG Zhuoyuan. Principles of engineering geological analysis[M]. 2nd ed. Beijing: Geological Publishing House, 1994. (in Chinese )
[2] 李滨,殷跃平,高杨,等. 西南岩溶山区大型崩滑灾害研究的关键问题[J]. 水文地质工程地质,2020,47(4):5 − 13. [LI Bin,YIN Yueping,GAO Yang,et al. Critical issues in rock avalanches in the karst mountain areas of southwest China[J]. Hydrogeology & Engineering Geology,2020,47(4):5 − 13. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.202003060 LI Bin, YIN Yueping, GAO Yang, et al. Critical issues in rock avalanches in the Karst Mountain areas of southwest China[J]. Hydrogeology & Engineering Geology, 2020, 47(4): 5-13. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.202003060
[3] YIN Yueping,WANG Fawu,SUN Ping. Landslide hazards triggered by the 2008 Wenchuan earthquake,Sichuan,China[J]. Landslides,2009,6(2):139 − 152. DOI: 10.1007/s10346-009-0148-5
[4] 刘乐,杨智,孙健,等. 安徽黄山市徽州区地质灾害危险性评价研究[J]. 中国地质灾害与防治学报,2021,32(2):110 − 116. [LIU Le,YANG Zhi,SUN Jian,et al. Study on risk assessment of geological hazards in Huizhou District,Huangshan City,Anhui Province[J]. The Chinese Journal of Geological Hazard and Control,2021,32(2):110 − 116. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.15 LIU (L /Y), YANG Zhi, SUN Jian, et al. Study on risk assessment of geological hazards in Huizhou district, Huangshan city, Anhui Province[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(2): 110-116. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.15
[5] 周粤,王运生,赵逊,等. 怒江支流迪麻洛河流域泥石流易发性评价[J]. 地质通报,2022,41(4):702 − 712. [ZHOU Yue,WANG Yunsheng,ZHAO Xun,et al. Susceptibility assessment of debris flow in Dimaluo River, branch of Nujiang River[J]. Geological Bulletin of China,2022,41(4):702 − 712. (in Chinese with English abstract) ZHOU Yue, WANG Yunsheng, ZHAO Xun, et al. Susceptibility assessment of debris flow in Dimaluo River, branch of Nujiang River[J]. Geological Bulletin of China, 2022, 41(4)702-712(in Chinese with English abstract)
[6] 殷跃平,李滨,张田田,等. 印度查莫利“2·7”冰岩山崩堵江溃决洪水灾害链研究[J]. 中国地质灾害与防治学报,2021,32(3):1 − 8. [YIN Yueping,LI Bin,ZHANG Tiantian,et al. Study on flood disaster chain of “2·7” ice rock landslide in Chamorri,India[J]. The Chinese Journal of Geological Hazard and Control,2021,32(3):1 − 8. (in Chinese with English abstract) YIN Yueping, LI Bin, ZHANG Tiantian, et al. Study on flood disaster chain of “2.7” ice rock landslide in Chamorri, India[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(3): 1-8. (in Chinese with English abstract)
[7] 许冲. 汶川地震滑坡分布规律与危险性评价[J]. 岩石力学与工程学报,2012,31(2):432. [XU Chong. Distribution law and risk assessment for Wenchuan earthquake-triggered landslides[J]. Chinese Journal of Rock Mechanics and Engineering,2012,31(2):432. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-6915.2012.02.024 XU Chong. Distribution law and risk assessment for Wenchuan earthquake-triggered landslides[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(2): 432. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-6915.2012.02.024
[8] 孟祥瑞,裴向军,刘清华,等. GIS支持下基于因子分析法的都汶路沿线地质灾害易发性评价[J]. 中国地质灾害与防治学报,2016,27(3):106 − 115. [MENG Xiangrui,PEI Xiangjun,LIU Qinghua,et al. GIS-Based susceptibility assessment of geological hazards along the road from Dujiangyan to Wenchuan by factor analysis[J]. The Chinese Journal of Geological Hazard and Control,2016,27(3):106 − 115. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2016.03.16 MENG Xiangrui, PEI Xiangjun, LIU Qinghua, et al. GIS-Based susceptibility assessment of geological hazards along the road from Dujiangyan to Wenchuan by factor analysis[J]. The Chinese Journal of Geological Hazard and Control, 2016, 27(3): 106-115. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2016.03.16
[9] 田春山,刘希林,汪佳. 基于CF和Logistic回归模型的广东省地质灾害易发性评价[J]. 水文地质工程地质,2016,43(6):154 − 161. [TIAN Chunshan,LIU Xilin,WANG Jia. Geohazard susceptibility assessment based on CF model and Logistic Regression models in Guangdong[J]. Hydrogeology & Engineering Geology,2016,43(6):154 − 161. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2016.06.24 TIAN Chunshan, LIU Xilin, WANG Jia. Geohazard susceptibility assessment based on CF model and Logistic Regression models in Guangdong[J]. Hydrogeology & Engineering Geology, 2016, 43(6): 154-161. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2016.06.24
[10] 张瑛. “5·12”汶川大地震震裂山体灾害勘查评价与治理设计方法研究[D]. 成都: 成都理工大学, 2009 ZHANG Ying. Study on exploration, evaluation and treatment design method of mountain disaster caused by “5·12” Wenchuan earthquake[D]. Chengdu: Chengdu University of Technology, 2009. (in Chinese with English abstract)
[11] NILSEN T, BRABB E. 18 slope-stability studies in the San francisco bay region, California[J]. Reviews in Engineering Geology, 1977: 233-244.
[12] 黄润秋. 中国西南岩石高边坡的主要特征及其演化[J]. 地球科学进展,2005,20(3):292 − 297. [HUANG Runqiu. Main characteristics of hign rock slopes in southwestern China and their dynamic evolution[J]. Advance in Earth Sciences,2005,20(3):292 − 297. (in Chinese with English abstract) DOI: 10.3321/j.issn:1001-8166.2005.03.005 HUANG Runqiu. Main characteristics of hign rock slopes in southwestern China and their dynamic evolution[J]. Advance in Earth Sciences, 2005, 20(3): 292-297. (in Chinese with English abstract) DOI: 10.3321/j.issn:1001-8166.2005.03.005
[13] 邓茂林, 许强, 郑光, 等. 基于离心模型试验的武隆鸡尾山滑坡形成机制研究[J]. 岩石力学与工程学报, 2016, 35(增刊1): 3024 − 3035 DENG Maolin, XU Qiang, ZHENG Guang, et al. Study on the formation mechanism of Jiweishan landslide in Wulong, Chongqing, China-based on centrifugal model test[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(Sup 1): 3024 − 3035. (in Chinese with English abstract)
[14] 陈国庆,黄润秋,周辉,等. 边坡渐进破坏的动态强度折减法研究[J]. 岩土力学,2013,34(4):1140 − 1146. [CHEN Guoqing,HUANG Runqiu,ZHOU Hui,et al. Research on progressive failure for slope using dynamic strength reduction method[J]. Rock and Soil Mechanics,2013,34(4):1140 − 1146. (in Chinese with English abstract) DOI: 10.16285/j.rsm.2013.04.040 CHEN Guoqing, HUANG Runqiu, ZHOU Hui, et al. Research on progressive failure for slope using dynamic strength reduction method[J]. Rock and Soil Mechanics, 2013, 34(4): 1140-1146. (in Chinese with English abstract) DOI: 10.16285/j.rsm.2013.04.040
[15] 韩用顺, 梁川, 崔鹏, 等. 地形条件对次生山地灾害易发性分析[J]. 四川大学学报(工程科学版), 2010, 42(增刊1): 15 − 21 HAN Yongshun, LIANG Chuan, CUI Peng, et al. Susceptibility of mountain hazards triggered by Wenchuan earthquake to topographic factors[J]. Journal of Sichuan University (Engineering Science Edition), 2010, 42(Sup 1): 15 − 21. (in Chinese with English abstract)
[16] 黄河清,赵其华. 汶川地震诱发文家沟巨型滑坡-碎屑流基本特征及成因机制初步分析[J]. 工程地质学报,2010,18(2):168 − 177. [HUANG Heqing,ZHAO Qihua. Basic characteristics and preliminary mechanism analysis of large scale rockslide-sturzstrom at Wenjiagou triggered by Wenchuan earthquake[J]. Journal of Engineering Geology,2010,18(2):168 − 177. (in Chinese with English abstract) DOI: 10.3969/j.issn.1004-9665.2010.02.003 HUANG Heqing, ZHAO Qihua. Basic characteristics and preliminary mechanism analysis of large scale rockslide-sturzstrom at Wenjiagou triggered by Wenchuan earthquake[J]. Journal of Engineering Geology, 2010, 18(2): 168-177. (in Chinese with English abstract) DOI: 10.3969/j.issn.1004-9665.2010.02.003
[17] 郑颖人,赵尚毅. 有限元强度折减法在土坡与岩坡中的应用[J]. 岩石力学与工程学报,2004,23(19):3381 − 3388. [ZHENG Yingren,ZHAO Shangyi. Application of strength reduction fem in soil and rock slope[J]. Chinese Journal of Rock Mechanics and Engineering,2004,23(19):3381 − 3388. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-6915.2004.19.029 ZHENG Yingren, ZHAO Shangyi. Application of strength reduction fem in soil and rock slope[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(19): 3381-3388. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-6915.2004.19.029
[18] 李文娟,邵海. 基于遥感影像多尺度分割与地质因子评价的滑坡易发性区划[J]. 中国地质灾害与防治学报,2021,32(2):94 − 99. [LI Wenjuan,SHAO Hai. Landslide susceptibility assessment based on multi-scale segmentation of remote sensing and geological factor evaluation[J]. The Chinese Journal of Geological Hazard and Control,2021,32(2):94 − 99. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.13 LI Wenjuan, SHAO Hai. Landslide susceptibility assessment based on multi-scale segmentation of remote sensing and geological factor evaluation[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(2): 94-99. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.13
[19] 郭学飞,王志一,焦润成,等. 基于层次分析法的北京市地质环境质量综合评价[J]. 中国地质灾害与防治学报,2021,32(1):70 − 76. [GUO Xuefei,WANG Zhiyi,JIAO Runcheng,et al. Comprehensive evaluation method of geological environment quality in Beijing based on AHP[J]. The Chinese Journal of Geological Hazard and Control,2021,32(1):70 − 76. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.01.10 GUO Xuefei, WANG Zhiyi, JIAO Runcheng, et al. Comprehensive evaluation method of geological environment quality in Beijing based on AHP[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(1): 70-76. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.01.10
-
期刊类型引用(1)
1. 张天祥,王艳霞,张雪珂,林钏,周汝良. 基于样本优化和机器学习的地质灾害气象风险预报模型研究——以云南省怒江州为例. 广西师范大学学报(自然科学版). 2025(01): 201-215 . 百度学术
其他类型引用(0)