Reliability analysis based on joint uncertainty: A case study of a rock slope in Tibet
-
摘要: 边坡稳定性一直是边坡安全的重点研究对象,针对边坡评价中常见的不确定性因素,可靠度分析是值得利用的方法。为评价某节理发育的岩质岸坡稳定性,通过有限元计算软件,结合现场勘探测绘数据,建立以边坡节理强度参数c、φ为输入变量,安全系数为输出变量的点估计(PEM)计算概率模型,计算结果表明:节理发育对该边坡变形具有明显控制作用;边坡整体可靠性较好,破坏概率极低。最后,通过蒙托卡罗法对可靠度结果进行验证,结果表明两种方法的计算结果不存在显著性差异。研究结果表明节理对岩质边坡稳定具有良好的敏感性,基于节理不确定性的点估计法分析边坡可靠度是一种有效的方法。Abstract: Slope stability has always been the key research object of slope safety. In view of the common uncertain factors in slope evaluation, reliability analysis is a valuable method. To evaluate the stability of a jointed rock slope, a point estimation (PEM) calculation probability model is established by using the finite element analysis software and the field survey data, The PEM method takes the joint strength parameters c, φ and the safety factor as input and output variables. The results suggest that, the joints have a significant control effect on the deformation of the slope; The overall reliability of slope is great, and the failure probability is extremely low. Finally, the reliability results are verified by Monte Carlo method, and the results show that there is no significant difference between the two methods. The results indicate that the joints have good sensitivity to the stability of rock slope, and the point estimation method based on joint uncertainty is an effective method to analyze the slope reliability.
-
Keywords:
- Jointed rock slope /
- reliability /
- baecher model /
- point estimate /
- Monte carlo
-
0. 引言
地质灾害链是指一系列在时间上有先后,空间上彼此相依,成因上相互关联、互为因果,呈连锁反应,依次出现的几种地质灾害组成的链式灾害体系[1]。在高山极高山区,暴雨-崩塌-滑坡-泥石流等灾害链较为常见,具有高位、高速、远程、气浪、容易链状成灾等特点,能够在高位滑坡启动以后转化成泥石流,再形成堰塞湖,这类地质灾害隐蔽性强、破坏力大,每年都会造成巨大的经济损失和人员伤亡[2]。
金沙江上游河道坡降大、地面切割强烈,山高谷深,为典型的高山峡谷区,V型河谷发育,谷坡地势陡峻、岩体支离破碎、分水岭狭窄,极易发生山体滑坡,该流域堵江滑坡受学者的广泛关注,研究成果主要集中于堵江滑坡的形成机制与稳定性评价[3-7]、灾害灾情应急监测与治理[8-10]、古滑坡堰塞湖溃坝堆积物粒度特征[11]等方面,多以典型堵江滑坡为研究对象,对具有成因联系的不同灾种之间的链式特征鲜有报道。在遥感方面,姚鑫等以金沙江及其支流定曲河谷为例,利用多源遥感影像的图斑特征配合实地考察,总结了该地区古滑坡、活动滑坡的识别特征[12],但该地区滑坡相对不发育,样本数量有待进一步补充。本文以光学遥感技术为研究手段,采用光学卫星、无人机等遥感数据,分析了金沙江流域直门达-石鼓段堵江滑坡地质灾害的分布规律,总结了堵江滑坡地质灾害的链式特征,以期为金沙江流域直门达-石鼓段堵江滑坡调查评价、风险评估、防灾减灾提供光学遥感基础资料。
1. 研究区概况
金沙江流域直门达—石鼓段位于青藏高原东部三江地区,地处冈瓦纳大陆与泛华夏大陆结合部位的金沙江结合带,地势总体西高东低,由横断山脉、云贵高原、四川盆地等三大地貌组成。研究区气候条件复杂多变,属季风气候区,夏秋季受印度洋和太平洋暖湿气流影响,降水较多,冬春季受亚欧大陆中心及蒙古高原干冷气团控制,降水稀少。整体上,从北西向南东,降水量逐渐增大,气温逐步升高。
研究区经历了漫长的地质变迁史、构造复杂,是古特提斯洋闭合形成的重要缝合带,地震活跃。构造行迹以近南北向的金沙江断裂带为代表,该断裂带由东、西边界断裂以及金沙江主断裂等3条主要断裂和其他近南北向次级断裂组成,宽约50~60 km,总体走向近南北向,断裂带内主要出露岩浆岩体、基性岩脉、古生代地层等地质体[13],空间上自西向东可分为陆缘火山弧带、构造混杂岩带、裂谷盆地带、蛇绿构造混杂岩带、前陆褶皱-冲断层带等5个地质构造单元。根据各地质构造单元内岩性特点、成因类型与岩土物理力学性质,可划分出层状碎屑岩类、层状碳酸盐岩类、结晶岩类三大岩体类型。
2. 数据来源
研究区选用高分一号、高分二号、WorldView-2、无人机航摄等多源卫星数据,数据分辨率包括中高分辨率2.0 m、0.8 m、0.5 m及0.2 m。其中,卫星影像云覆盖<1%,雪覆盖<5%,满足研究要求。
3. 研究方法
根据遥感影像上滑坡地质灾害隐患的色调、形态、纹理、阴影等特征,经过野外验证,建立遥感解译标志,开展金沙江上游堵江滑坡地质灾害隐患的遥感解译。研究区滑坡隐患点在影像上表现为绿色、灰色、灰白色、黄褐色、灰褐色等色调,具扇形、舌形、半圆形、梨形、倒梨形、簸箕形、矩形、复合型等平面形态,常见滑坡壁、剪切裂缝、拉裂缝、崩滑破坏、滑坡台阶、滑坡舌、滑坡鼓丘、后缘洼地等变形标志和微地貌。有时可见滑坡堰塞湖、河道弯曲等重要间接解译标志(图1b—e)。
4. 结果
4.1 堵江滑坡地质灾害隐患遥感调查结果分析
4.1.1 堵江滑坡地质灾害隐患基本特征
采用多源遥感数据,本次在金沙江上游共识别出87处堵江滑坡地质灾害隐患点(图1、表1)。其中已发生滑动造成堵江事件的9处,局部发生变形、整体未滑动的潜在滑坡78处。9处已发生的堵江滑坡中堵塞金沙江干流有6处,分别是特米滑坡(HP11)、岗达滑坡(HP26)、王大龙滑坡(HP27)、旺各滑坡(HP28)、白格滑坡(HP30)、昌波滑坡(HP71);堵塞金沙江支流有3处,分别是通错滑坡(HP16)、上卡岗滑坡(HP31)、苏洼龙滑坡(HP70)。根据光学遥感解译堵江滑坡地质灾害隐患点的平面面积,估算其体积,按照滑坡规模划分标准,金沙江上游堵江滑坡地质灾害隐患主要为大型、特大型,其中大型40处,特大型47处。以堵江滑坡地质灾害隐患失稳后可能形成的灾害链为划分依据,研究区潜在滑坡-堵江灾害链54处、崩塌-滑坡-堵江灾害链20处、滑坡-泥石流-堵江灾害链13处。
表 1 金沙江上游滑坡地质灾害隐患光学遥感解译一览表Table 1. Optical remote sensing interpretation list of landslide in upper reach of Jinsha River滑坡编号 滑坡名称 地理位置 前缘高程/m 后缘高程/m 相对高差/m 规模 金沙江
位置可能灾害链 HP01 藏曲河口滑坡 江达县波罗乡 2915 3179 265 特大型 右岸 滑坡-堵江 HP02 圭利滑坡 江达县波罗乡 2892 3285 393 特大型 右岸 滑坡-堵江 HP03 肖莫久滑坡 江达县波罗乡 2890 3643 754 特大型 右岸 滑坡-堵江 HP04 色拉滑坡 贡觉县沙东乡 2650 3460 810 特大型 右岸 滑坡-堵江 HP05 敏都滑坡 贡觉县敏都乡 2648 3023 375 大型 左岸 滑坡-堵江 HP06 麦绿席滑坡 贡觉县敏都乡 2637 3316 679 特大型 左岸 崩塌-滑坡-堵江 HP07 瓦堆村滑坡 贡觉县敏都乡 2640 2823 183 大型 右岸 滑坡-堵江 HP08 拍若滑坡 巴塘县拉哇乡 2558 2785 227 大型 左岸 崩塌-滑坡-堵江 HP09 通中拉卡滑坡 巴塘县拉哇乡 2568 3554 986 特大型 左岸 崩塌-滑坡-堵江 HP10 曲引郎滑坡 巴塘县拉哇乡 2577 3096 520 大型 右岸 崩塌-滑坡-堵江 HP11 特米滑坡 巴塘县拉哇乡 2541 3013 472 特大型 左岸 崩塌-滑坡-堵江 HP12 打呷龙滑坡 巴塘县拉哇乡 2424 2660 236 大型 左岸 崩塌-滑坡-堵江 HP13 萨荣白拉滑坡 巴塘县昌波乡 2371 2708 337 大型 右岸 滑坡-堵江 HP14 霍荣滑坡 巴塘县昌波乡 2314 2855 541 大型 右岸 崩塌-滑坡-堵江 HP15 尼吕滑坡 得荣县贡波乡 2234 2639 405 大型 右岸 崩塌-滑坡-堵江 HP16 通错滑坡 贡觉县雄松乡 3304 3930 626 特大型 左岸支沟 滑坡-堵江 HP17 岗苏村滑坡 贡觉县雄松乡 3206 3545 338 大型 右岸支沟 滑坡-泥石流-堵江 HP18 扎木塘隆滑坡 贡觉县雄松乡 2973 3590 617 特大型 右岸支沟 滑坡-泥石流-堵江 HP19 下缺所滑坡 贡觉县雄松乡 2737 2847 110 大型 右岸支沟 滑坡-泥石流-堵江 HP20 约德同滑坡 曲麻莱县巴干乡 3824 4117 293 大型 左岸 滑坡-堵江 HP21 沙东乡2#滑坡 贡觉县敏都乡 2920 3626 706 特大型 右岸支沟 滑坡-泥石流-堵江 HP22 优巴村滑坡 白玉县绒盖乡 3228 4167 939 大型 左岸支沟 滑坡-泥石流-堵江 HP23 宁巴村滑坡 白玉县波罗乡 3059 3265 206 大型 右岸支沟 滑坡-泥石流-堵江 HP24 尼增滑坡 巴塘县拉哇乡 2536 2986 450 特大型 右岸 滑坡-泥石流-堵江 HP25 丁巴滑坡 贡觉县罗麦乡 2771 2970 198 大型 左岸 滑坡-堵江 HP26 岗达滑坡 芒康县索多西乡 2436 3362 926 特大型 右岸 崩塌-滑坡-堵江 HP27 王大龙滑坡 巴塘县昌波乡 2370 2860 490 特大型 左岸 滑坡-堵江 HP28 旺各滑坡 巴塘县昌波乡 2481 3270 788 特大型 左岸 滑坡-堵江 HP29 波罗寺滑坡 江达县波罗乡 2940 3072 132 大型 右岸支沟 滑坡-堵江 HP30 白格滑坡 江达县波罗乡 2893 3704 811 特大型 右岸 滑坡-堵江 HP31 上卡岗滑坡 江达县波罗乡 3026 3488 461 大型 左岸 滑坡-堵江 HP32 卡岗村滑坡 江达县波罗乡 2925 3558 632 特大型 左岸 滑坡-堵江 HP33 染莫滑坡 江达县波罗乡 3431 3949 518 特大型 右岸 滑坡-堵江 HP34 沙丁麦2#滑坡 白玉县金沙乡 2962 3346 384 特大型 右岸 滑坡-堵江 HP35 扎永滑坡 白玉县金沙乡 2971 3230 260 大型 左岸 滑坡-堵江 HP36 金沙江乡滑坡 白玉县金沙乡 2965 3167 202 大型 左岸支沟 滑坡-堵江 HP37 若翁洞滑坡 白玉县金沙乡 2988 3383 395 特大型 左岸支沟 滑坡-堵江 HP38 探戈滑坡 白玉县金沙乡 3170 3950 780 特大型 右岸支沟 滑坡-泥石流-堵江 HP39 甘棒滑坡 白玉县金沙乡 2968 3797 829 特大型 左岸 滑坡-堵江 HP40 冷普东滑坡 江达县岩比乡 2978 3227 249 大型 右岸 滑坡-堵江 HP41 沃达滑坡 江达县岩比乡 2991 3653 662 特大型 右岸 滑坡-堵江 HP42 然贡寺滑坡 江达县岩比乡 3009 3549 541 特大型 右岸 滑坡-堵江 HP43 格纳滑坡 江达县岩比乡 3007 3378 372 特大型 右岸 滑坡-堵江 HP44 延若滑坡 江达县岗托镇 3019 3346 327 特大型 右岸 滑坡-堵江 HP45 底巴滑坡 江达县岗托镇 3055 3456 401 特大型 右岸 滑坡-堵江 HP46 来格村滑坡 江达县岗托镇 3080 3710 630 特大型 左岸 滑坡-堵江 HP47 波桥滑坡 江达县岗托镇 3012 3655 643 特大型 右岸 滑坡-堵江 HP48 燃灯寺滑坡 江达县汪布顶乡 3062 3510 448 大型 右岸 滑坡-堵江 HP49 热公下游滑坡 江达县汪布顶乡 3248 3482 234 大型 右岸支沟 滑坡-泥石流-堵江 HP50 热公滑坡 江达县汪布顶乡 3075 3687 612 特大型 右岸 崩塌-滑坡-堵江 HP51 林切滑坡 江达县汪布顶乡 3120 3542 422 特大型 左岸 滑坡-堵江 HP52 汪布顶滑坡 江达县汪布顶乡 3108 3909 801 特大型 右岸 崩塌-滑坡-堵江 HP53 索达滑坡 石渠县麻呷乡 3203 3520 317 大型 左岸 滑坡-堵江 HP54 俄波贡村滑坡 石渠县麻呷乡 3617 3749 132 大型 左岸支沟 滑坡-泥石流-堵江 HP55 东仲滑坡 石渠县麻呷乡 3565 3829 264 大型 右岸支沟 滑坡-泥石流-堵江 HP56 擦塘滑坡 石渠县麻呷乡 3264 3471 207 大型 右岸 滑坡-堵江 HP57 挡底滑坡 石渠县奔达乡 3359 3710 351 特大型 右岸 滑坡-堵江 HP58 莫罗宫滑坡 石渠县奔达乡 3976 4203 227 大型 右岸支沟 滑坡-泥石流-堵江 HP59 色考滑坡 石渠县奔达乡 3472 3875 403 大型 右岸支沟 滑坡-泥石流-堵江 HP60 色巫滑坡 石渠县奔达乡 3375 4200 825 特大型 右岸 滑坡-堵江 HP61 热科塘滑坡 石渠县麻呷乡 3285 3762 477 大型 右岸 滑坡-堵江 HP62 麦巴村滑坡 贡觉县敏都乡 2646 3437 791 特大型 右岸 滑坡-堵江 HP63 上缺所村对岸滑坡 贡觉县敏都乡 2653 3555 902 特大型 左岸 崩塌-滑坡-堵江 HP64 下松哇滑坡 巴塘县拉哇乡 2563 2923 360 大型 右岸 崩塌-滑坡-堵江 HP65 下松哇对岸滑坡 巴塘县拉哇乡 2566 2956 390 大型 左岸 崩塌-滑坡-堵江 HP66 刀出顶滑坡 芒康县竹巴龙乡 2470 2949 479 特大型 左岸 滑坡-堵江 HP67 琼冲阁对岸滑坡 芒康县竹巴龙乡 2426 2799 373 特大型 左岸 滑坡-堵江 HP68 老衣对岸滑坡 芒康县竹巴龙乡 2445 2810 365 特大型 左岸 滑坡-堵江 HP69 巴窝滑坡 曲麻莱县巴干乡 3878 4300 422 大型 左岸 崩塌-滑坡-堵江 HP70 苏洼龙滑坡 芒康县苏哇龙乡 2815 3139 324 大型 左岸支沟 滑坡-堵江 HP71 昌波滑坡 巴塘县昌波乡 2395 2750 355 特大型 左岸 崩塌-滑坡-堵江 HP72 明惹滑坡 石渠县真达乡 3522 3919 398 特大型 右岸 滑坡-堵江 HP73 羊拉登桥滑坡 巴塘县昌波乡 2310 2683 373 大型 左岸 滑坡-堵江 HP74 吉布卡堆滑坡 巴塘县地巫乡 2367 3236 869 特大型 左岸 滑坡-堵江 HP75 拿荣滑坡 得荣县徐龙乡 2191 2630 439 特大型 左岸 崩塌-滑坡-堵江 HP76 格亚顶滑坡 得荣县徐龙乡 2195 2903 708 特大型 右岸 崩塌-滑坡-堵江 HP77 莫丁滑坡 得荣县徐龙乡 2196 2754 558 特大型 左岸 崩塌-滑坡-堵江 HP78 雪拉滑坡 芒康县古龙乡 2041 2466 425 特大型 左岸 滑坡-堵江 HP79 维热纳滑坡 石渠县奔达乡 3378 3537 159 大型 右岸 滑坡-堵江 HP80 民则口滑坡 石渠县奔达乡 3390 3603 213 大型 右岸 滑坡-堵江 HP81 奔达滑坡 石渠县奔达乡 3411 3576 165 大型 左岸 滑坡-堵江 HP82 僧提滑坡 石渠县奔达乡 3422 3686 264 大型 右岸 滑坡-堵江 HP83 牙们儿滑坡 香格里拉市尼西乡 1982 2353 371 大型 左岸 崩塌-滑坡-堵江 HP84 供仁滑坡 德钦县拖顶傈僳族乡 1935 2681 746 特大型 右岸 滑坡-堵江 HP85 达叶滑坡 察雅县真达乡 3511 3784 273 大型 右岸 滑坡-堵江 HP86 嘎达滑坡 玉树市安冲乡 3706 3971 265 大型 右岸 滑坡-堵江 HP87 沙东滑坡 贡觉县沙东乡 3290 4170 880 特大型 右岸 滑坡-堵江 4.1.2 堵江滑坡地质灾害隐患空间分布规律
金沙江上游堵江滑坡地质灾害隐患分布在青海省、四川省、西藏自治区、云南省等4个省12个市(县),其中西藏自治区42处、四川省40处,少量分布在青海省(3处)和云南省(2处),具体可划分出6个集中分布区,自北向南依次为:石渠县奔达乡集中分布区(Ⅰ)、石渠县麻呷乡集中分布区(Ⅱ)、江达县汪布顶乡—白玉县绒盖乡集中分布区(Ⅲ)、贡觉县沙东乡—雄松乡集中分布区(Ⅳ)、巴塘县拉哇乡集中分布区(Ⅴ)、芒康县竹巴龙乡—巴塘县地巫乡集中分布区(Ⅵ)(图1),均位于河流深切割的高山峡谷区,岸坡陡峻,地质构造强烈。从金沙干流及其支流左右两岸分布特征来看,金沙江干流分布堵江滑坡70处,其中左岸分布31处,右岸分布39处;金沙江支流分布17处,其中左岸支流分布6处,右岸支流分布11处。
4.2 堵江滑坡地质灾害隐患链式特征
4.2.1 滑坡-堵江灾害链
此类灾害链是金沙江上游主要的潜在堵江滑坡灾害链之一,分布在金沙干流左右两侧,表现为坡度较陡的新生滑坡或老滑坡复活,滑坡滑动后堵塞河道,形成堰塞湖,堰塞湖满库后导致堰塞坝溃决,下游形成洪水灾害,淹没村庄、耕地、公路、桥梁等,造成人员和财产损失。部分滑坡后缘可再次发生滑动,进一步造成人员和财产损失。以白格滑坡(HP30)、色拉滑坡(HP04)为代表。
采用2020年6月无人机航测影像,分辨率为0.2 m,对色拉滑坡进行精细解译,滑坡平面面积约1.5×106 m2,按照平均厚度50 m估算,体积约75×106 m3,为特大型滑坡,滑坡前缘高程约2650 m,后缘高程约3460 m,高差810 m,根据滑坡的变形特征,将滑坡分为Ⅰ、Ⅱ、Ⅲ三个变形区(图2)。
Ⅰ区分布在滑坡体前部,该区域变形较强烈,坡体发生多次滑动。根据变形特征又将该区分为Ⅰ-1、Ⅰ-2两个亚区。Ⅰ-1亚区主要分布在滑坡前部上游侧,属于强变形区,前缘受江水冲刷发生溜滑,溜滑区呈灰白色带状,沿坡向长约40~90 m。Ⅰ-1亚区内发育多级圈椅状陡坎,形成HT-3、HT-4两个次级滑体。其中,HT-3次级滑体分布在Ⅰ-1亚区中前部,后缘陡坎呈驼峰状,可见多条拉裂缝,宽约2~5 m,两侧剪切裂缝明显,次级滑体长约490 m,宽约530 m,平面面积约25×104 m2,按照平均厚度30 m估算,体积约7.50×106 m3,为大型滑坡,该次级滑体坡整体破碎,多处发生滑塌,呈逐步解体破坏的趋势。HT-4次级滑体位于HT-3次级滑体前部,后缘滑坡陡坎明显,错落高约10~20 m,坡体呈散体状,发育多级拉裂缝,次级滑体长约220 m,宽约360 m,平面面积约7.5×104 m2,按照平均厚度30 m估算,体积约2.25×106 m3,为大型滑坡,该次级滑体整体滑动可能性大。Ⅰ-2亚区分布于滑坡前部下游侧,坡体破碎,裂缝发育,根据变形特征可以明显划分出HT-1、HT-2两个次级滑体。HT-1次级滑体分布在下游侧前缘,次级滑体长约260 m,宽约230 m,平面面积约5.5×104 m2,按照平均厚度30 m估算,体积约1.65×106 m3,该次级滑体后缘滑坡壁明显,宽约15 m,有整体滑动的趋势。HT-2次级滑体平面形态呈舌形,次级滑体长约560 m,宽约130 m,平面面积约8.2×104 m2,按平均厚度30 m估算,体积约2.46×106 m3,次级滑体后缘发育圈椅状拉裂陡坎,坡体发育多级拉裂缝,可能发生多级滑动。
Ⅱ区主要分布在滑坡体中部,该区域变形相对Ⅰ区较弱。该区坡体冲沟、陡坎发育,光学影像上未见明显的变形迹象。
Ⅲ区分布在滑坡后部,后缘坡体发生局部溜滑,形成明显的滑坡壁,该区斜坡整体较稳定。
通过光学遥感解译分析,色拉滑坡受河流下切作用的影响,在重力作用下局部表层发生滑塌变形,堆积体沿斜坡分布,堆积体前缘受金沙江水冲刷。滑坡所在斜坡为斜向坡,受金沙江断裂带影响,岩体较破碎,节理裂隙发育,滑坡前缘发生强烈的变形,发生多级滑动,形成多级拉裂陡坎,如进一步扩展,将有可能形成牵引式的大规模滑坡,破坏模式类似白格滑坡,形成滑坡-堵江灾害链。
4.2.2 崩塌-滑坡-堵江灾害链
此类灾害链分布在高位斜坡区,在高陡斜坡处常发育危岩带,受几组控制性结构面作用,危岩体发生高位崩塌,岩体解体沿斜坡、沟道形成松散堆积体,剖面呈陡-缓-陡形态。在自重力作用下,加之后部局部不断再次发生崩塌加载,前缘受金沙江河流侵蚀破坏,滑坡体整体缓慢蠕变。在强降雨、地震等诱发因素下,可能发生大规模崩塌,并以巨大能量冲击加载在下部斜坡上,形成崩塌-滑坡-堵江灾害链,以汪布顶滑坡(HP52)为代表。
采用2020年6月无人机航测影像,分辨率为0.2 m,对滑坡进行精细解译,滑坡平面面积约80×104 m2,按照平均厚度40 m估算,体积约3.2×106 m3,为特大型滑坡,滑坡前缘高程约3108 m,后缘高程约3909 m,高差801 m,根据滑坡的变形特征,将滑坡分为Ⅰ、Ⅱ、Ⅲ、Ⅳ四个变形区(图3a)。
Ⅰ区主要分布在滑坡前缘,该区域变形较强烈,受公路切坡、金沙江侧蚀作用的影响,坡体发生多次滑动,形成4处坡表溜滑体,滑塌体平面面积介于0.5×104~2.8×104 m2,总面积约5.7×104 m2,按平均厚度30 m估算,体积约1.71×106 m3。
Ⅱ区主要分布在滑坡中下部。根据地形和变形特征又将该区分为Ⅱ-1、Ⅱ-2两个亚区。Ⅱ-1亚区分布于滑坡上游侧,该区域发育1条拉裂缝、3处坡表溜滑体,滑塌体平面面积介于0.3×104~0.7×104 m2,总面积约1.7×104 m2,按平均厚度30 m估算,体积约51×104 m3,整体变形相对较小,可见由块石碎石组成的崩滑堆积体。Ⅱ-2亚区主要分布在滑坡下游侧,该区地形较陡,坡体变形强,以滑塌变形为主,发育1处坡表溜滑体,面积约1×104 m2,按平均厚度30 m估算,体积约30×104 m3。
Ⅲ区主要分布在滑坡中部。该区域受基岩锁固作用变形较弱,整体较稳定,坡体局部发生滑塌现象,可见2处坡表溜滑体,滑塌体平面总面积约0.4×104 m2,按平均厚度30 m估算,体积约12×104 m3。
Ⅳ区主要分布在滑坡中后部(图3b)。根据变形特征可细分Ⅳ-1、Ⅳ-2两个亚区。Ⅳ-1亚区分布在滑坡后部上游侧,该区后缘发育危岩带,高约30~50 m,前方有大量崩塌堆积体,变形区平面面积介于0.3×104~12×104 m2,总面积约13×104 m2,按平均厚度30 m估算,体积约3.9×106 m3,以块石、碎石为主,最大块径约20 m。Ⅳ-2亚区分布在滑坡中后部下游侧,该区变形强烈,发育危岩带和多条拉裂缝,可分为崩塌变形区和拉裂变形区。崩塌变形区变形最强烈,地形陡峻,基于高分辨率的无人机航空倾斜摄影三维模型,危岩带发育4组优势结构面(图3c),分别为:311°~5°∠48°~55°、130°~150°∠58°~70°、32°~50°∠65°~80°、200°~225°∠81°~87°,在其控制下危岩体失稳发生高位崩塌,产生大量崩塌堆积体,崩塌变形区平面面积介于0.6×104~5.1×104 m2,总面积约11×104 m2,按平均厚度30 m估算,体积约3.30×106 m3,危岩体下方堆积大量松散块石,最大块石块径约30 m。拉裂变形区地形相对较缓,后部发育多条NE向拉裂缝,裂缝长度不一,未完全贯通。根据高分辨率无人机影像初步估计,Ⅳ区的危岩带平面面积约18×104 m2,按平均厚度30 m估算,体积约5.4×106 m3,规模为特大型,位于“陡-缓-陡”斜坡形态的上部较陡处(图3b),该危岩体失稳后可将以巨大的冲积载荷直接加载于斜坡中部较缓处,预计直接影响面积约16×104 m2,涉及Ⅱ-1区后缘、Ⅲ区、Ⅳ-2区前缘。
根据本次解译分析,汪布顶滑坡所在斜坡为斜向坡,坡体前缘受金沙江侧蚀作用和公路坡体开挖等人类工程活动的影响,局部发生滑塌变形。滑坡后缘斜坡较陡,发育大量的危岩带,一旦发生较大规模的崩滑失稳,滑塌体将会铲刮中下部的堆积体,形成大规模的滑坡碎屑流,形成崩塌-滑坡-堵江灾害链。
4.2.3 滑坡-泥石流-堵江灾害链
此类灾害链分布在金沙江两岸的支沟中,高位滑坡在坡体中下部和沟道内堆积了大量的松散堆积物,在强降雨的影响下很可能转化为沟道泥石流冲出支沟,挤压、堵塞金沙江,以探戈滑坡(HP38)为代表。
采用WorldView-2卫星影像,分辨率为0.5 m,对滑坡进行精细解译,滑坡纵剖面呈“陡-缓-陡”阶梯状地形,平面面积约1.8×106 m2,按照平均厚度50 m估算,体积约90×106 m3,为特大型滑坡,滑坡前缘高程约3170 m,后缘高程约3950 m,高差780 m,结合滑坡各部位的变形特征,将滑坡分为Ⅰ、Ⅱ、Ⅲ三个变形区(图4)。
Ⅰ区分布在滑坡体前部,该区域地形较陡,可见多处坡体滑动。根据变形强弱程度又将该区分为Ⅰ-1、Ⅰ-2两个亚区。Ⅰ-1亚区主要分布在滑坡前部支沟上游侧,该区滑坡变形强烈,属于强变形区,前缘临空条件较好,在重力作用下发生多处溜滑。溜滑区呈灰色带状,沿坡向长约150~200 m,后缘发育拉裂缝,裂缝断续分布,未整体贯通;溜滑体直达沟底,在坡脚形成扇形、近扇形松散堆积体。根据拉裂缝、溜滑体的分布和变形特征,可识别出HT-1、HT-2、HT-3三个次级滑体。HT-1次级滑体分布在Ⅰ-1亚区左前部,后缘陡坎明显,早期滑动面交错呈树枝状分布,次级滑体长约280 m,宽约130 m,平面面积约3.5×104 m2,按照平均厚度30 m估算,体积约1.05×106 m3,为大型滑坡,该次级滑体坡整体破碎,呈逐步解体破坏的趋势。HT-2次级滑体分布在Ⅰ-1亚区中前部,后缘滑坡陡坎明显,错落高约50 m,坡体左侧呈散体状,发育多级拉裂缝,变形最为强烈,次级滑体长约500 m,宽约500 m,平面面积约20×104 m2,按照平均厚度30 m估算,体积约为6×106 m3,为大型滑坡。HT-3次级滑体分布在Ⅰ-1亚区右前部,后缘滑坡陡坎高约15~20 m,次级滑坡呈哑铃状,中后部整体破碎,溜滑现象明显,次级滑体长约550 m,宽约240 m,平面面积约14×104 m2,按照平均厚度30 m估算,体积约4.2×106 m3,为大型滑坡。Ⅰ-2亚区分布于滑坡前部下游侧,坡体整体较稳定,前缘地形变陡处可见多处小规模溜滑。
Ⅱ区主要分布在滑坡体中部,该区域地势较平缓、坡体无明显变形迹象,从光学影像上看斜坡整体较稳定。
Ⅲ区分布在滑坡后部,根据变形强度的不同可识别出HT-4次级滑坡体。HT-4次级滑体分布在Ⅲ区左前部,滑坡壁清晰可见,后缘滑坡陡坎错落高约80 m,坡体整体破碎,变形强烈,与HT-1、HT-2、HT-3次级滑坡体有逐步贯通,整体下滑的趋势。HT-4次级滑体长约770 m,宽约630 m,平面面积约40×104 m2,按平均厚度30 m估算,体积约12×106 m3,为特大型滑坡。
根据本次解译分析,探戈滑坡在重力作用下局部发生溜滑变形,堆积体沿斜坡分布,前缘直达沟底。滑坡所在斜坡为顺向坡,滑坡体前缘和后缘均发生强烈变形,表现为滑坡陡坎、拉裂缝、多级滑动。据遥感影像可知,坡体下部支沟中泥石流堆积物已十分丰富,物源主要以溜滑堆积体为主,面积约16×104 m2,在沟口可见新近泥石流扇汇入金沙江。在地震、强降雨等因素诱发下,潜在滑坡体可能失稳大规模下滑至坡前沟道中堵塞沟道,并形成泥石流物源,探戈滑坡所在沟道中流域面积约50.11×106 m2,主沟两侧发育大小支沟近22条(图5),探戈滑坡位于沟口附近,具备堵塞19条支沟汇水水量的条件,相当于可拦截该流域约86%的水流量,在冰雪融化、强降雨等充足水流量条件下,探戈滑坡可最终演化成滑坡-泥石流-堵江灾害链。
4.3 金沙江断裂带与堵江滑坡地质灾害链的关系
从地理位置上看,金沙江断裂带明显控制了金沙江干流直门达-石鼓段的平面展布,断裂带附近构造混杂岩发育、岩体结构复杂、斜坡完整性差、地形坡度大,具备发生堵江滑坡的有利条件。受印度板块向欧亚板块的NNE向持续挤压作用,金沙江断裂带新生代以来发生持续隆升和挤压构造变形,是现今地球表面地形地貌和地质构造演化最复杂、构造活动最强烈的地区之一[6],新构造运动在研究区主要表现为断裂复活、地震等,具有继承性、新生性和节奏性。中国地震台网公开资料显示,自2012年以来研究区发生2≤Ms<3地震11起,3≤Ms<4地震58起,4≤Ms<5地震19起,5≤Ms<6地震4起,地震活动会破坏山体的稳定性,使岩体变得更加松散、破碎,节理、裂隙更加发育,增加了地质灾害发生的风险。
金沙江断裂带大致以NE向的巴塘断裂和NW向的德钦—中甸—大具断裂为界,可分为北、中、南3段[13],各段活动周期、强度存在差异性。地形地貌特征显示,断裂带活动性及强度由东向西、由北向南呈递增趋势,具体表现在:①研究区具整体间歇性抬升,发育两级夷平面,同时以金沙江为界,同一级夷平面差异升降明显,海拔高度表现出西高东低、南高北低的特点;②金沙江“V”型河谷的相对高差由北往南分别为500~1000 m、1000~2000 m,增幅近两倍。GPS水平速度场显示,金沙江断裂带北段活动性不明显,中段及南段表现出明显的右旋走滑运动,现今滑动速率分别为4.9 mm/a和5.5 mm/a,中段及南段存在应变积累,具有地震危险性[14]。采用面波振幅谱能量衰减特征拟合研究区及邻区已发生地震的震源深度一般在5~10 km范围内[15],均属浅源地震,破坏性强。前人研究表明,在断层活动性相对较弱的金沙江断裂带北段,白格滑坡-碎屑流-堵江灾害链的形成仍与多次地震作用关系密切[5],那么在断层活动性较强的金沙江断裂带中段、南段,应变积累更快,地震作用可能相对频繁,在巴塘以南的金沙江两岸有利斜坡区可能更具发生大规模堵江滑坡的风险,应加强研究和监测。
5. 结论
(1)根据遥感影像上滑坡地质灾害隐患的色调、平面形态、变形标志和微地貌,建立遥感解译标志,在金沙江上游直门达—石鼓段共识别出堵江滑坡地质灾害隐患点87处,分析了其基本特征和空间分布规律。
(2)以堵江滑坡地质灾害隐患点失稳后可能形成的链式灾害为基础,划分出滑坡-堵江灾害链、崩塌-滑坡-堵江灾害链、滑坡-泥石流-堵江灾害链等3类灾害链,基于光学遥感技术分析了其变形特征和链式特征。
(3)研究区地处金沙江活动断裂带,地震等地质作用与研究区堵江滑坡地质灾害隐患的稳定性关系密切,活动性较强、应变积累更快的中段和南段可能更具备发生堵江滑坡灾害链的区域地质环境背景。
-
表 1 岩体参数取值
Table 1 Parameter values of rock mass
材料 弹性模量/MPa 泊松比 φ/(°) c/MPa 容重/(kN·m−3) 均值 标准差 均值 标准差 强风化花岗岩 10000 0.28 38 2 0.12 0.01 23.5 弱风化花岗岩 13000 0.26 42 2 0.15 0.02 25.0 微新花岗岩 20000 0.25 49 3 0.50 0.04 27.0 表 2 结构面几何特征值
Table 2 Geometric characteristc values of structural plane
边坡分层 节理组 产状 迹长 密度/(条·m−2) 平均倾向/(°) 平均倾角/(°) 分布类型 Fisher常数 均值/m 分布类型 强风化层 J1 235 67 Fisher 95 8 负指数分布 0.20 J2 125 45 45 5.5 0.25 J1 245 74 50 15 负指数分布 0.06 弱风化层 J2 125 58 Fisher 40 5 0.03 J3 105 30 90 5 0.05 微新基岩 J1 236 78 Fisher 30 6 负指数分布 0.02 J2 135 55 110 4 0.01 表 3 结构面强度特征值
Table 3 Strength characteristic values of structural plane
节理组号 φ/(°) c/MPa 法向刚度/(MPa·m−1) 切向刚度/(MPa·m−1) 抗拉强度/MPa 均值 标准差 均值 标准差 J1 30 2 0.04 0.002 8000 2500 0 J2 33 2 0.08 0.003 10000 3000 0 J3 38 3 0.10 0.005 13000 3800 0 表 4 敏感性计算结果
Table 4 Sensitivity calculation results
节理强度参数 岩体强度参数 稳定性系数 c/MP φ/(°) c/MPa φ/(°) 0.04 30 0.15 42 1.35 0.034 30 0.15 42 1.33 0.046 30 0.15 42 1.38 0.04 24 0.15 42 1.23 0.04 36 0.15 42 1.45 0.04 30 0.144 42 1.34 0.04 30 0.156 42 1.37 0.04 30 0.150 36 1.27 0.04 30 0.150 48 1.44 表 5 敏感性指标计算结果
Table 5 Sensitivity index calculation results
节理强度参数 岩体强度参数 c/MP φ/(°) c/MP φ/(°) 0.0370 0.1630 0.0222 0.1259 表 6 可靠度对比分析结果
Table 6 Results of reliability comparison analysis
稳定性参数 JEF-PEM点估计 MCS蒙托卡罗 平均值 1.344 1.590 标准差 0.0573 0.0854 可靠性指数 6.007 6.907 破坏概率 0.001 <0.001 -
[1] 陈祖煜, 陈立宏, 王玉杰, 等. 滑坡和建筑物抗滑稳定分析中的可靠度分析和分项系数设计方法[C]//陈祖煜. 水利水电工程风险分析及可靠度设计技术进展. 北京: 中国水利水电出版社, 2010: 27–39 CHEN Zuyu, CHEN Lihong, WANG Yujie, et al. Reliability analysis and partial factor design methods for slope stability[C]// CHEN Zuyu. Proceedings of the symposium on risk analysis and reliability based design for water resources and hydropower projects. Beijing: China Water Power Press, 2010: 27–39. (in Chinese)
[2] 刘超, 袁颖, 左朝晖, 等. 考虑土参数自相关距离影响的单桩可靠性设计方法[J]. 水文地质工程地质,2020,47(3):122 − 127. [LIU Chao, YUAN Ying, ZUO Zhaohui, et al. Pile reliability design considering the influence of soil auto-correlation distance[J]. Hydrogeology & Engineering Geology,2020,47(3):122 − 127. (in Chinese with English abstract) [3] 李涛, 刘国栋, 王聪. 基于可靠度理论的尾矿坝失稳概率及敏感性分析[J]. 中国地质灾害与防治学报,2019,30(3):81 − 86. [LI Tao, LIU Guodong, WANG Cong. Instability probability and sensitivity analysis of tailings dam based on reliability theory[J]. The Chinese Journal of Geological Hazard and Control,2019,30(3):81 − 86. (in Chinese with English abstract) [4] JIANG S H, LI D Q, ZHANG L M, et al. Slope reliability analysis considering spatially variable shear strength parameters using a non-intrusive stochastic finite element method[J]. Engineering Geology,2014,168:120 − 128. DOI: 10.1016/j.enggeo.2013.11.006
[5] LI Y J, HICKS M A, NUTTALL J D. Comparative analyses of slope reliability in 3D[J]. Engineering Geology,2015,196:12 − 23. DOI: 10.1016/j.enggeo.2015.06.012
[6] ZHAO L H, YU C H, CHENG X, et al. A method for seismic stability analysis of jointed rock slopes using Barton-Bandis failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences,2020,136:104487. DOI: 10.1016/j.ijrmms.2020.104487
[7] LU R L, WEI W, SHANG K W, et al. Stability analysis of jointed rock slope by strength reduction technique considering ubiquitous joint model[J]. Advances in Civil Engineering,2020,2020:8862243.
[8] 冯开帅, 姜谙男, 吴洪涛, 等. 顺层非贯通节理边坡稳定性及破坏模式研究[J]. 公路工程,2020,45(6):52 − 58. [FENG Kaishuai, JIANG Annan, WU Hongtao, et al. Study on failure mode and stability of slope based on intermittent joints[J]. Highway Engineering,2020,45(6):52 − 58. (in Chinese with English abstract) [9] CHANG Y, CHANG L S, REN F Q. Strength anisotropy of jointed rock slope considering mining damage: A numerical analysis method[J]. Geomatics, Natural Hazards and Risk,2020,11(1):2587 − 2614. DOI: 10.1080/19475705.2020.1856200
[10] 陈国良, 廖国华. 岩体节理网络的计算机模拟[J]. 有色金属(矿山部分),1989,41(6):23 − 28. [CHEN Guoliang, LIAO Guohua. Computer simulation of rock mass joint network[J]. Nonferrous Metals (Mining),1989,41(6):23 − 28. (in Chinese) [11] 王宏, 陶振宇. 边坡稳定分析的节理网络模拟原理及工程应用[J]. 水利学报,1993,24(10):20 − 27. [WANG Hong, TAO Zhenyu. The principle of joint network simulation in the stability analysis of rock slope and its engineering application[J]. Journal of Hydraulic Engineering,1993,24(10):20 − 27. (in Chinese with English abstract) DOI: 10.3321/j.issn:0559-9350.1993.10.003 [12] 严豪, 宋彦辉, 陈子玉. 基于Voronoi节理模型的碎裂岩边坡稳定性分析[J]. 中国地质灾害与防治学报,2018,29(1):34 − 39. [YAN Hao, SONG Yanhui, CHEN Ziyu. Stability analysis of broken rock slope based on Voronoi joint net model[J]. The Chinese Journal of Geological Hazard and Control,2018,29(1):34 − 39. (in Chinese with English abstract) [13] 李源亮, 任光明, 刘艳领, 等. 基于节理网络有限元的岩质边坡稳定性分析[J]. 长江科学院院报,2019,36(1):78 − 83. [LI Yuanliang, REN Guangming, LIU Yanling, et al. Stability analysis of rock slope based on finite elements with joint network[J]. Journal of Yangtze River Scientific Research Institute,2019,36(1):78 − 83. (in Chinese with English abstract) DOI: 10.11988/ckyyb.20170823 [14] 肖术, 吴顺川, 高永涛, 等. 基于PEM-JFEM方法的节理岩质边坡稳定性评价[J]. 工程科学学报,2015,37(7):844 − 850. [XIAO Shu, WU Shunchuan, GAO Yongtao, et al. Jointed rock slope stability evaluation based on PEM- JFEM method[J]. Chinese Journal of Engineering,2015,37(7):844 − 850. (in Chinese with English abstract) [15] 张宜杰, 任光明, 常文娟, 等. 节理岩质边坡稳定性概率分析[J]. 成都理工大学学报(自然科学版),2021,48(2):235 − 241. [ZHANG Yijie, REN Guangming, CHANG Wenjuan, et al. Probabilistic analysis of stability of jointed rock slope[J]. Journal of Chengdu University of Technology (Science & Technology Edition),2021,48(2):235 − 241. (in Chinese with English abstract) [16] PRZEWLOCKI J, ZABUSKI L, WINKELMANN K. Reliability analysis of sea cliff slope stability by point estimate method[J]. IOP Conference Series:Materials Science and Engineering,2019,471:042003. DOI: 10.1088/1757-899X/471/4/042003
[17] 许湘华, 曲广琇, 方理刚. 基于节理几何参数不确定性的边坡可靠度分析[J]. 中南大学学报(自然科学版),2010,41(3):1139 − 1145. [XU Xianghua, QU Guangxiu, FANG Ligang. Reliability analysis of rock slope based on uncertainty of joint geometric parameters[J]. Journal of Central South University (Science and Technology),2010,41(3):1139 − 1145. (in Chinese with English abstract) [18] 王双, 陈征宙, 吴强, 等. 基于节理产状不确定性的边坡稳定性及敏感度分析[J]. 岩土工程学报,2013,35(2):348 − 354. [WANG Shuang, CHEN Zhengzhou, WU Qiang, et al. Stability and sensitivity analysis of slopes based on uncertainty of joint orientations[J]. Chinese Journal of Geotechnical Engineering,2013,35(2):348 − 354. (in Chinese with English abstract) [19] 王双. 节理产状概率模型研究及其在产状分组和岩坡不确定分析中的应用[D]. 南京: 南京大学, 2013 WANG Shuang. Research on distribution model of joint orientations and its application to joint set clustering and rock slope uncertainty analysis[D]. Nanjing: Nanjing University, 2013. (in Chinese with English abstract)
[20] 申矫健. 基于参数不确定性的岩质边坡稳定性分析方法研究[D]. 赣州: 江西理工大学, 2014 SHEN Jiaojian. Study on stability analysis method of rock slope based on parameter uncertainty[D]. Ganzhou: Jiangxi University of Science and Technology, 2014. (in Chinese with English abstract)
[21] BAECHER G B, LANNEY N A, EINSTEIN H H. Statistical description of rock properties and sampling[J]. 18th U S Symposium on Rock Mechanics, USRMS 1977,1977:5 − 8.
[22] 王宇, 李晓, 刘帅, 等. 岩体离散裂隙网络稳定性计算的节理有限元法[J]. 岩石力学与工程学报,2013,32(增刊 2):3337 − 3345. [WANG Yu, LI Xiao, LIU Shuai, et al. Stability calculation for discrete fracture network of rock masses based on joint finite element method[J]. Chinese Journal of Rock Mechanics and Engineering,2013,32(Sup 2):3337 − 3345. (in Chinese with English abstract) [23] CACAS M C, LEDOUX E, DE MARSILY G, et al. Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model[J]. Water Resources Research,1990,26(3):479 − 489.
[24] DERSHOWITZ W S, EINSTEIN H H. Characterizing rock joint geometry with joint system models[J]. Rock Mechanics and Rock Engineering,1988,21(1):21 − 51. DOI: 10.1007/BF01019674
[25] ROSENBLUETH E. Point estimates for probability moments[J]. Proceedings of the National Academy of Sciences of the United States of America,1975,72(10):3812 − 3814. DOI: 10.1073/pnas.72.10.3812
[26] 李侃, 巨能攀. 基于蒙特卡洛方法的边坡可靠性评价[J]. 中国地质灾害与防治学报,2014,25(1):23 − 27. [LI Kan, JU Nengpan. Integrated application of Monte-Carlo simulation for landslide reliability analysis[J]. The Chinese Journal of Geological Hazard and Control,2014,25(1):23 − 27. (in Chinese with English abstract) -
期刊类型引用(3)
1. 张曦,骆建文,潘俊义,刘斌,杨楠. 陕西榆林毛乌素沙漠南缘风积沙的湿陷规律及其影响因素. 中国地质灾害与防治学报. 2024(04): 75-84 . 本站查看
2. 刘鑫,覃泽华,汪潇杰,兰恒星. 循环荷载下饱和黄土的剪切特性与破坏模式. 吉林大学学报(地球科学版). 2024(05): 1604-1614 . 百度学术
3. 任森,邓龙胜,范文,徐良鑫,刘少鹏. 陕西富平黄土震陷特性及震陷小区划. 中国地质灾害与防治学报. 2023(04): 134-143 . 本站查看
其他类型引用(6)