库区公路K243滑坡形成机制及治理措施

2023-08-01王有昌

科技和产业 2023年11期

王有昌

(国家林业和草原局昆明勘察设计院, 昆明 650031)

随着澜沧江上游河段梯级水电站陆续开发建设和正常蓄水运营,库区公路路堑开挖引发了较多滑坡、崩塌等地质灾害问题,库岸再造致使岸坡坍塌变形,影响公路路基稳定和运营安全。因此,研究库区公路滑坡的成因机制和变形机理非常关键,而成功有效的治理措施和经验对保障库区公路路基稳定和运营安全至关重要。

众多学者对滑坡的成因机制、变形机理、治理措施、数值模拟、监测与预警等开展了大量研究,取得了较多成果。殷跃平等[1]、成永刚等[2]对滑坡的识别与机理分析、滑坡勘察与稳定性评价、防治工程设计理论与实践等进行了深入研究。王浩等[3-4]对滑坡孕育过程进行数值仿真反演分析,揭示滑坡形成的复杂力学机制,开展施工动态模拟预测及过程优化控制。许泰等[5]、吴冠仲[6]、丁恒[7]对滑坡影响因素、形成原因、变形机理进行研究,评价滑坡其稳定性。申永江等[8]、肖世国和何洪[9]、叶四桥和潘新恩[10]和邹盛堂和戴自航[11]提出了不同的双排桩计算模型,通过三维数值分析软件对双排桩的滑坡加载方式和受力体系进行研究,分析桩身内力及变形。于洋等[12]、李松林等[13]和顾中华[14]对双排桩结构设计参数进行分析,建立无连梁和有连梁双排抗滑桩有限差分模型,研究不同锚固深度组合时对双排抗滑桩变形、内力和抗滑力的影响规律,讨论双排刚架圆截面抗滑桩的计算方法。李松林等[15]、江强强等[16]和朱权威等[17]、李厚芝等[18]研究了库区滑坡成因机制、变形机理,分析库岸滑坡在不同水位波动、降雨强度及其联合作用下的变形演化特征和稳定性变化规律,提出滑坡治理对策。

现依托澜沧江黄登水电站库区沿江公路K243滑坡实例工程,对滑坡成因机制及治理措施进行研究,运用GEO5软件建立K243滑坡计算分析模型,对抗滑桩变形、抗剪与抗弯承载力进行分析,评价滑坡治理整体稳定性,对库区公路滑坡分析与治理具有重要意义。

1 工程实例

1.1 滑坡概况

澜沧江沿江公路K243+200～K243+450段位于维西县维登乡以南约4 km,距黄登水电站坝址约75 km,公路沿库区左岸岸坡展布,黄登水电站正常蓄水水位为1 619.00 m,K243位于库尾,回水水位按1 622.00 m考虑。沿江公路从古滑坡前缘部位通过,因路堑边坡开挖形成连续、坡高大于20 m的临空面,2011年7月中旬,K243+240～K243+390段边坡发生变形坍塌,清方后路堑仍持续滑坍,边坡变形加剧,古滑坡中前部位发生整体滑移,形成滑坡。K243滑坡由古滑坡及新生工程滑坡组成,滑坡全貌如图1所示。

图1 K243滑坡全貌

古滑坡后缘壁呈“圈椅”状,后缘壁坎高为50～70 m,水平滑移距离为50～100 m;滑坡后缘高程约1 803 m,前缘剪出口高程约1 584 m,地面落差约219 m;滑坡体平面周长为1 102 m,平均宽度为335 m,平均厚度为21 m,滑体最大厚度为27.1 m,属深层滑坡;滑坡体积为189.65×104m3,属巨型滑坡。古滑坡两侧基本以基岩为界,中后部存在多处台地,古滑坡周界不清晰、变形痕迹不显著,近期无活动痕迹。

新生工程滑坡后缘壁呈“圈椅”状,后缘壁坎高为2～3 m,后缘裂缝宽度为30～40 cm,前缘面状渗水,坡体坍塌严重;滑坡后缘高程约1 698 m,前缘剪出口高程约1 621 m,地面高差约76 m;滑坡体平面周长为504 m,平均宽度为155 m,平均厚度为11 m,最大厚度为16.2 m,属中层滑坡;滑坡体积为16.67×104m3,属中型滑坡。新生工程滑坡两侧基本以变形土体为界,周界清晰,变形特征明显、痕迹显著,为近期活动滑坡,其现状变形特征如图2和图3所示。

图2 新生工程滑坡后缘特征

图3 新生工程滑坡前缘特征

1.2 滑坡地质环境特征

1.2.1 地形地貌

滑坡区位于云贵高原西部、澜沧江上游地带的横断山脉南段,地处滇西北纵谷山原地貌单元,属澜沧江高山峡谷深切割地貌,构造侵蚀强烈,地形起伏大,岸坡上缓下陡。

1.2.2 地质构造

滑坡区位于青藏、滇缅、印尼巨型“歹”字形构造体系东支中段偏北与三江南北向经向构造体系复合部位,以北西向构造带与南北向构造带为主体的区域构造基本骨架。滑坡区附近无明显断裂构造。

1.2.3 地层岩性

地层岩性以第四系覆盖层和侏罗系砂泥岩地层组成。滑体物质为褐灰色、黄灰色混碎石粉质黏土、碎石土、碎块石层,碎石成分以砂岩、砂泥岩为主;滑带为黄灰、紫灰色含砾粉质黏土。滑床主要为黄灰色、灰色的碎石土和紫红色、紫灰色的全-强风化砂泥岩组成。

1.2.4 地震及区域稳定性

滑坡区地震基本烈度为7度,地震动峰值加速度为0.10g(g为重力加速度),地震动反应谱特征周期为0.45 s。滑坡区域地壳稳定性分区为Ⅲ区,属次稳定区。

1.3 水文地质条件

滑坡区地下水属第四系松散孔隙水、基岩裂隙水。孔隙水分布于第四系冲洪积、坡残积层及滑坡堆积层中,主要接受大气降水、地表水补给,水量随季节变化大;基岩裂隙水主要赋存于侏罗系上统坝注路组砂泥岩地层中,以大气降水补给为主。滑坡区中前部有多处泉点出露,稳定水位埋深为9.0～27.0 m,水位高差极大。

2 滑坡成因机制及变形破坏模式

2.1 滑坡成因机制

滑坡形成机制包括地质环境因素和诱发因素两个方面,地质环境因素包括地形地貌、岩土性质、地质构造、水文地质等;诱发因素包括持续强降雨及水流作用、地震、人类活动等。

2.1.1 地质环境因素

滑坡地貌上属于澜沧江高山峡谷岸坡地貌,岸坡陡峻,斜坡整体平衡状态脆弱,澜沧江长期冲刷,创造了临空面和滑动空间,也提供了滑动势能。地层岩性以第四系覆盖层和侏罗系砂泥岩地层组成,覆盖层结构松散,基岩风化强烈、岩性破碎,倾向河谷拉张卸荷裂隙十分发育,为滑坡形成创造了物质和结构条件。受澜沧江深大断裂影响,加剧了基岩地层破碎程度,局部小断裂、揉皱构造发育,为地下水渗透和运动创造了条件。地下水长期向下径流、侵蚀,软化岩土体,使其力学性质降低,是滑坡形成的重要条件。

2.1.2 诱发因素

持续降雨或强降雨,使斜坡内动水压力激增,岩土体重量增加,抗滑力降低,是滑坡发生的主要诱发因素。滑坡区处于强震区,地震能显著增加斜坡岩土体的下滑力,是滑坡发生偶然性、突发性的诱发因素。沿江公路从古滑坡体上通过,路堑开挖极为不利;地表开荒、植被破坏,进村道路施工及矿区堆载,也是滑坡变形的诱发因素。

2.2 滑坡变形破坏模式

澜沧江峡谷岸坡的地形地貌、地层岩性和地质构造为古滑坡形成创造了势能、物质和结构条件。澜沧江河谷深切,江水冲刷,形成超高临空面,为坡体滑动提供了空间;地表水和地下水长期作用,沿风化壳下渗,导致岩土体软化,且重量增加,加之暴雨、地震等多方面的影响,逐步形成软弱带,致使古滑坡孕育演化形成,因此古滑坡体属于牵引推动复合型滑坡。滑坡钻探勘察揭示了古滑坡特征。

沿江公路通过古滑坡体前缘,古滑坡水平滑移距离大,堆积体扰动强烈、结构松散,虽然古滑坡已经趋于稳定,但表层堆积体结构松散,开挖极易滑塌。路堑边坡切脚开挖,坡脚支撑力严重削弱,在古滑坡体前部形成连续、高陡的临空面,改变了坡体应力状态和地下水渗透路径,坡面大面积渗水、坍塌;在路堑开挖施工期间,滑坡区历经30年一遇的降雨量,持续降雨是滑坡形成的重要影响因素;在局部清方后持续滑坍,使岩土体进一步软化,抗剪强度持续降低,导致边坡失稳,沿最薄弱面形成贯通滑带,最终形成新生工程滑坡,属于牵引式滑坡。新生工程滑坡是多种因素共同作用形成的。

古滑坡前缘处于水库回水区,水库蓄水后,滑坡体前缘将被淹没,库区水位调节、波浪侵蚀可能引起库岸再造,将进一步改变坡体应力状态,使前缘抗滑力大幅降低,必将引起滑坡复活。

3 滑坡稳定性与剩余下滑力分析

3.1 滑坡计算参数反演分析

3.1.1 滑坡参数反演分析模型

根据滑坡现状变形特征、古滑坡变形行迹,以及钻探揭露滑动面位置,结合滑坡主轴地质剖面,划分条块,通过分析和试算合理确定新生工程滑坡和古滑坡的计算模型,如图4所示。

注:1～14为条块N的编号图4 滑面参数反演计算模型

3.1.2 滑坡计算参数反演

依据土工试验、原位测试所确定的滑体重度、滑带抗剪强度指标,结合经验数据类比,滑动面按折线形考虑,采用传递系数法进行参数反演分析,综合确定滑体及滑带不同工况下滑体重度和滑带参数建议值,如表1所示。

表1 滑体及滑动面参数反演分析成果

3.2 滑坡稳定性分析

滑坡稳定状态根据滑坡稳定系数进行划分,滑坡稳定系数Fs<1.00,滑坡稳定状态为不稳定;1.00≤Fs<1.05,为欠稳定;1.05≤Fs<1.15,为基本稳定;Fs≥1.15,为稳定。滑坡稳定性分析考虑正常工况(天然)、非正常工况Ⅰ(饱和)、非正常工况Ⅱ(地震)等三种工况,并考虑低水位和高水位(按回水水位1 622 m考虑),对滑坡进行多工况组合条件下的稳定性分析评价。

由表2可知,新生工程滑坡位于蓄水位以上,低水位和高水位状态下时稳定性系数一致;天然状态下滑坡为基本稳定状态,天然+地震状态下为欠稳定状态,饱和状态下为不稳定状态。而古滑坡前缘位于正常蓄水位以下,高水位工况下滑坡的稳定性分析比低水位工况下要低;天然状态下,在低水位、高水位或地震作用时滑坡稳定系数为1.053～1.199,滑坡处于基本稳定至稳定状态;饱和状态下,低水位时为基本稳定,高水位或地震作用时为不稳定。

表2 滑坡主剖面稳定性计算成果

3.3 滑坡推力计算

3.3.1 滑坡剩余下滑力分析

根据《公路滑坡防治设计规范》,滑坡防治工程设计稳定安全系数Ks,新生工程滑坡在正常工况时Ks取1.20,非正常工况Ⅰ时Ks取1.15,非正常工况Ⅱ时Ks取1.10;古滑坡在正常工况时Ks取1.25,非正常工况Ⅰ时Ks取1.15,非正常工况Ⅱ时Ks取1.10。根据滑坡变形行迹和稳定性分析,综合黄登水电站正常运行时库区公路路基稳定,考虑最不利条件高水位工况,新生工程滑坡和古滑坡的滑坡推力按正常工况、非正常工况Ⅰ、非正常工况Ⅱ进行计算。

由图5、图6可知,新生工程滑坡的剩余下滑力在非正常工况Ⅱ(饱和+地震+高水位)条件下的剩余下滑力最大,正常工况条件下剩余下滑力最小,故考虑非正常工况Ⅱ为最不利设计工况,古滑坡的剩余下滑力也基本一致。新生工程滑坡的剩余下滑力曲线正常工况下低于其他工况,其中条块5在非正常工况Ⅱ下的最大剩余下滑力为2 393.62 kN,大于其他工况的2 335.91 kN、1 591.36 kN。古滑坡各工况条件下的剩余下滑力曲线基本一致,条块4的剩余下滑力最大,条块6以后逐渐降低趋于平缓。

图5 新生工程滑坡不同工况下的剩余下滑力

图6 古滑坡不同工况下的剩余下滑力

3.3.2 滑坡设计推力确定

根据滑坡防治工程设计稳定安全系数下的滑坡剩余下滑力分析,确定新生工程滑坡及古滑坡设置支挡工程部位的最大水平设计推力T。K243滑坡治理考虑多排桩分级支挡,其中第1排抗滑桩位于新生工程滑坡的第3条块,设计推力T1-1=1 761 kN,而位于古滑坡第8条块,设计推力T1-2=2 678 kN,桩前剩余抗力781 kN。第2排桩设计推力分段计算,并考虑填土反压填方荷载,第2排抗滑桩位于新生工程滑坡前缘外,填土后设计推力T2-1=510 kN,古滑坡第11条块设计推力T2-2=305 kN,填土主动土压力为1 132 kN。

4 滑坡治理设计与稳定性分析

4.1 滑坡治理原则

K243滑坡新生工程滑坡为中型滑坡,古滑坡为巨型滑坡,公路以路堑型式通过滑坡前部,新生工程滑坡危害程度为严重,古滑坡危害程度为特严重;沿江公路等级为三级公路,且位于黄登水电站库区,滑坡防治安全等级按Ⅱ级考虑。沿江公路从古滑坡前缘部位通过,不具备路线改移绕避条件,对滑坡进行综合治理。滑坡治理按永久性防护工程考虑,本着安全可靠、技术可行、经济合理、动态设计的原则,充分考虑黄登水电站库区水位变化对滑坡稳定性的影响,确保库区沿江公路路基整体稳定。

4.2 滑坡治理措施

K243滑坡综合考虑采取多排抗滑桩分级支挡,经多方案计算分析论证,采用双排刚架抗滑桩、锚索桩板墙、填土反压、混凝土挡墙、浆砌片石护坡、植物防护和截排水等措施进行综合治理,滑坡治理平面布置和主剖面设计如图7和图8所示。

图7 滑坡治理工程布置平面

图8 滑坡治理主剖面设计断面

4.2.1 路基右侧治理措施

在沿江公路路基右侧布设桩径D=1.8 m的双排圆形截面抗滑桩,桩间距为5～6 m,前后排桩间距5 m呈品字型布置,桩顶设冠梁纵横向连接,形成双排刚架抗滑桩,前排桩顶至护坡间设现浇挡土板;冠梁1 m处设6×15.24型预应力锚索,锁定荷载为600kN;坡脚设高度3.5 m的C25混凝土挡土墙,墙顶至抗滑桩间采用30 cm厚的M10浆砌片石护坡,预留泄水孔。小里程方向设置2.0×2.5 m锚索抗滑桩,大里程方向布设高4 m的C20片石砼抗滑挡墙。坡面上设置仰斜排水孔,长度为10～20 m,孔间距为4 m。对坡面采用乔灌草混合植物防护。对坡面裂缝用黏土进行封闭,沿通村路设多条排水沟截排地表水。

4.2.2 路基下方治理措施

在路基左侧下方老路内侧布设2.0×3.0 m锚索抗滑桩,桩间距为5～6 m,距桩顶1 m、4 m处设两排6×15.24型预应力锚索,设置桩后预制挡土板。板后填土反压形成路堤,第一级边坡坡率为1∶1.5,高度为9.1 m,第二级坡率为1∶2,高为8 m,边坡间设置2 m宽平台,1 622 m高程以下采用透水性材料填筑。

4.3 滑坡治理设计分析

结合K243滑坡治理主剖面设计,运用GEO5软件中抗滑桩设计模块进行计算分析,建立双排圆截面抗滑桩和锚索抗滑桩模型,根据其所受荷载,分析其变形、抗剪承载力、抗弯承载力和嵌固端岩土体的承载力,并进行配筋校核;锚索施加在抗滑桩上,推力通过作用在抗滑桩上锚索传递到稳定地层中,计算模型如图9所示。

图9 滑坡计算模型3D视图

4.3.1 双排刚架抗滑桩计算分析

新生工程滑坡双排桩按设计推力值进行验算,由验算结果图10和图11可知,前排桩剪力最大值为735.82 kN/m,弯矩最大值为2 494.73 kN·m/m,位移最大值为38.18 mm;后排桩剪力最大值为743.18 kN/m,弯矩最大值为2 575.22 kN·m/m,位移最大值为38.53 mm,滑面处结构位移为9.1 mm。古滑坡按设计推力值适当折减进行验算,前排桩剪力最大值为878.57 kN/m,弯矩最大值为4 360.15 kN·m/m;后排桩剪力最大值为519.48 kN/m,弯矩最大值为2 292.65 kN·m/m,滑面处结构位移为3.2 mm。

图10 前排桩内力

图11 后排桩内力

4.3.2 锚索抗滑桩计算分析

新生工程滑坡锚索抗滑桩按设计推力值进行验算,设计时考虑桩后主动土压力,而古滑坡分段设计推力值较小。由验算结果图12可知,桩身剪力最大值为1 204.98 kN/m,弯矩最大值为6 606.15 kN·m/m,位移最大值为27.7 mm,滑面处结构位移为8.9 mm。

图12 锚索抗滑桩内力

4.4 滑坡治理稳定性分析

K243滑坡治理稳定性分析考虑施工前、锚索桩板墙+填土反压、双排抗滑桩、桩前路堑开挖4 m、路堑全部开挖完成等5个不同施工过程,在高水位条件下按正常工况、非正常工况Ⅰ、非正常工况Ⅱ等三种工况,对滑坡治理施工前后整体稳定性评价。图13为新生工程滑坡稳定性计算简图。

图13 新生工程滑坡稳定性计算

由图14可知,新生工程滑坡非正常工况Ⅰ时,在锚索抗滑桩+填土反压施工完成后稳定系数为1.08,双排抗滑桩施工完成后为1.22,路堑开挖完成后达到1.28,治理完成后滑坡稳定系数提高约32.0%;滑坡施工完成后整体稳定性正常工况时为1.33,非正常工况Ⅱ时为1.19,滑坡稳定系数提高分别为30.8%、23.2%,均满足滑坡设计稳定系数。而古滑坡施工完成后在正常工况、非正常工况Ⅰ、非正常工况Ⅱ下稳定系数分别为1.29、1.20和1.12,滑坡稳定性提高比例分别为15.6%、15.8%和13.1%。综合分析认为,K243滑坡治理稳定性满足规范要求,新生工程滑坡在不同施工过程稳定性提高较快,而古滑坡相对平缓。

注:1为施工前;2为锚索桩板墙+填土反应;3为双排抗滑桩;4为桩前路堑开挖4 m;5为路堑全部开挖完成图14 新生工程滑坡不同施工过程的稳定性系数

5 滑坡治理效果评价

K243滑坡是澜沧江沿江公路规模和治理投资较大的典型滑坡,治理施工自2012年6月开始,至2014年12月底竣工,图15为滑坡治理竣工后全景效果。

图15 滑坡治理竣工后全景(2014年12月)

根据施工过程滑坡变形监测,在下部锚索桩板墙与回填反压实施后,有效地抑制了新生工程滑坡变形;路基右侧双排抗滑桩施工完成,坡体抗力逐渐提高,地表位移减小逐渐趋于稳定。2017年11月28日黄登水电站导流洞下闸蓄水,2018年8月17日水库蓄水至正常蓄水位1 619.00 m,双排刚架抗滑桩与锚索抗滑桩受力稳定正常,路面稳定;2022年1月,经现场排查滑坡整体无变形迹象,图16为库区蓄水运营后滑坡全景。因此,沿江公路K243滑坡治理是成功的,对新生工程滑坡重点治理的同时兼顾古滑坡治理,有效保障沿江公路运营安全,积累了较为可靠的成功经验。