乡村展销大棚滑坡成因及稳定性评价
2024-02-29许德龙
许德龙, 殷 浩
(1.济南城建集团有限公司, 济南 250031; 2.山东省鲁岳资源勘查开发有限公司, 山东 泰安 271021)
近年来,各地都在施行美丽乡村建设。建设过程中涉及诸多工程建设活动,在丘陵、山区易出现边坡滑塌现象。国内外都有边坡滑塌所造成的生命财产损失的例子,因此对于边坡滑塌治理的重要性也就不言而喻。
近年来,国内外众多学者对滑坡的诱导因素及稳定性评价都做了大量研究[1-12]。宋国虎等[1]通过对五家坟滑坡进行现场调查分析,研究了滑坡的变形破坏特征和形成机制;张泽方等[2]通过对六盘水市水城区滑坡研究,提出了一种考虑最优影响因素组合的滑坡易发性评价模型;李厚芝等[5]通过对三峡库区滑坡成因机制研究,分析了强降雨条件下滑坡的稳定性;许泰等[6]通过对永登县苦水镇大路村滑坡得研究,定性、定量分析了滑坡的稳定性;吕文斌等[7]以西宁市张家湾滑坡为研究对象,结合极限平衡法和Midas/GTS有限元数值模拟方法,研究了滑坡的形成机理、发展过程;黄海韵[12]通过对尖山营滑坡的分析研究,得出了滑坡的稳定状态。现以莱芜区和庄镇马家峪村乡村旅游区建设所引起的山体滑坡为研究对象,分析了滑坡形成的原因,并对其稳定性进行评价,提出治理建议,以保证边坡安全。
1 工程概况
为实施乡村振兴战略,莱芜区和庄镇马家峪村在其后方山坡建设展销大棚,在建设区靠近斜坡边缘陆续填方堆载扩建。截至2022年8月位于斜坡顶部平台堆载完成并建设展销大棚设施,自此填方区上部展销大棚地表场外堆载平台斜坡开始出现缓慢开裂变形,至2022年年底变形加剧形成滑坡。滑坡变形主要表现为地表开裂和展销大棚场区内房屋破坏,第一排房屋结构性完全破坏,不能使用;场地靠近坡体一侧围墙全部垮塌,滑坡体下方道路内侧排水沟破坏,路面局部鼓胀开裂变形。此外,滑坡体前缘及东侧为乡村道路,因道路修建造成2~3 m切坡,形成临空面且无支护措施。
经过实地调查、测绘,目前整个滑坡体范围均存在不同程度的变形开裂,其中后缘变形最为严重,后缘呈现台阶式下错,最大高度约2.5 m,形成下陷区域,最宽约10 m,裂缝可见深2~3 m,滑坡下部公路内侧坡体见少许放射性裂缝,宽约10 cm,可见深约10 cm。
2 区域地质环境
2.1 地形地貌
治理区位于县城东南侧,区内地貌类型为山地,由于侵蚀作用强烈,形成了“脊尖”“坡陡”“谷窄”的地貌特征,尖棱状山脊与“V”形沟谷交替出现,山坡陡峻,常见崩塌和滑坡。
滑坡区整体地势南高北低,地形较为简单,海拔高度为498~550 m,相对高差52 m,山体坡坡度为15°~35°。滑坡区上部为展销大棚,下部为盘山乡道。区内植被茂盛。
2.2 地层岩性
场区岩土体类型主要为第四系残坡积、砂质板岩、硅质角岩。残坡积土及碎石直接覆盖于基岩之上,结构松散,物理力学性质差。砂质板岩、硅质角岩节理裂隙与孔隙较为发育,尤以风化节理最为发育,在强风化层之下,岩体物理力学性质较差。
场区内出露的岩层主要为二叠系上统包尔敖包组。岩性为紫色、灰绿色砂质板岩,与二云母角岩、硅质角岩呈渐变沉积接触或呈捕虏体分布其中。砂质板岩与硅质角岩不等厚互层,工程地质的特性为坚硬、半坚硬和软质岩类,按结构可分成块状到薄层状。
第四系土层分布广泛,其中斜坡中及斜坡中下部、斜坡凹部等部位较厚,岩性主要为亚黏土、亚砂土、砂砾石、碎块石土,表面堆积腐殖层,厚度一般为2~5 m。其不良工程地质现象主要表现为不均匀沉陷。在陡坡地段而局部松散堆积体较厚时,易出现滑坡。
中-厚层状砂质板岩与硅质角岩互层组主要分布于场区内大部分地区,易风化崩解,且遇水易软化,常形成较厚的风化层,沿其风化层,特别是在地下水、地表水渗透的作用下,易产生滑动。
2.3 水文地质条件
地下水主要为松散堆积层和基岩裂隙及构造破碎带潜水,风化裂隙发育深度在45 m左右,大气降水的垂直入滲直接补给地下水。山前洼地和沟谷地带主要分布孔隙潜水及局部的承压水,除接受大气降水补给外,还接受基岩裂隙水的补给。后沿沟谷向下游径流,在低洼地带受阻而形成沼泽湿地或以泉水的形式出露。
3 滑坡体特征
3.1 滑坡区地貌形态及边界特征
勘查区内地貌类型为山地,整体地势南高北低,展销大棚位于山体半山腰平台处,由于平台边缘填土及构筑物加载,导致岩土层滑动,挤压下部产生变形,形成滑坡。滑坡周界明显,滑坡前缘以下方盘山道路为界,滑坡后缘则以养殖场建筑边缘区域为界,两侧分别为西侧建筑构边缘延伸及东侧道路边缘。
整个滑坡地形呈上陡下缓,坡度为15°~23°。坡体植被以松树、灌木、杂草为主,原生植被覆盖率约60%。滑坡体前缘高程为512 m,后缘高程为541 m,最大高差29 m。滑坡体横向宽度约120 m,主滑方向为180°,顺主滑方向长约95 m,滑体厚度为1.5~6 m,体积约5.0×104m3,为小型推移式滑坡,如图1所示。
3.2 滑坡体变形特征
3.2.1 滑坡前缘
滑坡前缘为道路切坡,切坡坡度为83°,高1.5~2 m,表面裸露处风化严重,上部40 cm为耕植土,下部为卵石土、砂质板岩等,可见树根。局部切坡出现坍塌现象,导致部分水沟及道路挤压变形,上方土体表面也出现多条放射性裂隙,走向88°,宽约0.2 m,深0.7 m,长2~5 m,如图2所示。
图1 滑坡全貌
图2 滑坡前缘裂缝
3.2.2 滑坡后缘
滑坡体位移后出现明显分界面,滑坡壁整体呈圈椅状暴露,高1.2~1.5 m,同时滑坡体与滑坡壁被拉开形成反坡地形的沟槽状低洼封闭地形,导致局部水泥硬化路面及展销大棚下陷垮塌,出现多个台坎及一条沉陷带。在展销大棚外出现3条张拉裂缝,呈圆弧形,与滑坡壁大致平行,其中主裂缝LF1长约110 m,宽0.5~1 m,深约0.7 m,贯通整个后缘滑坡周界的裂缝;LF2长约50 m,宽约0.5 m,深1~2 m;LF3长约35 m,宽约0.3~0.5 m,深1~2 m,如图3所示。
3.3 滑坡体的物质结构特征
3.3.1 滑体特征
滑坡体主要为滑坡堆积层,滑体物质主要以人工填筑土以及岩石风化坡积、残积成因土质为主,成分以黏土及碎石等为主,次为砂质板岩、硅质角岩碎屑、粉粒,具有结构松散、孔隙率高、力学性质差等特性。堆积层在滑坡后缘最厚,主要集中在后缘堆积平台位置约7 m,然后向滑体前缘厚度逐渐减小,滑坡体前缘厚度1.5~3 m。
图3 滑坡后缘裂缝
3.3.2 滑带土特征
根据现场勘查,滑带位于上部松散堆积层底部与下伏强风化砂质板岩岩土接触面,滑带土物质主要为含卵石粉质黏土,卵石含量35%左右,粒径多小于3 cm,总体呈潮湿、近饱和软塑状[13]。
3.3.3 滑床特征
滑床为二叠系上统包尔敖包组砂质板岩夹硅质角岩。岩质较软,具有饱水软化、脱水开裂特征。岩层较为平缓,倾角13°左右,倾向南西180°左右,与滑坡方向180°基本一致,为顺层滑坡。根据风化程度可分为强风化带和弱风化带。①强风化带:岩体极为破碎,成碎石状,该层较为平缓,分布较均匀,平均厚度4 m;②弱风化带:该层风化程度较低,岩体较为完整,呈中厚层状。
4 滑坡成因及形成机制
4.1 地质环境条件
4.1.1 地形地貌
滑坡区域位于山体中部斜坡处,总体呈上陡下缓趋势,坡度为15°~35°。
4.1.2 地层岩性
滑坡体物质成分主要为人工填筑弱膨胀性黏性土和含卵石粉质黏土,结构松散,抗剪强度和抗风化能力较低,结构均匀性差,具饱水软化脱水开裂特征。根据现场勘查,滑体底部黏性土层相对富集,原生黏性土土质均匀,结构紧密,下伏滑床为风化砂质板岩。由于黏性土与砂质板岩透水性相对较差,滑体受强降雨影响,导致层内土体软化,抗剪强度大大降低,形成潜在软弱滑带,在外因的作用下滑体沿软弱带产生滑移[14-16]。
4.2 诱发因素
(1)为工程活动破坏坡体平衡。早期展销大棚建设在场地南侧边缘填方堆载渣石后未对坡体采取支护措施,且修建展销大棚基础设施对坡体加载;再加上乡村道路修建在斜坡下方切坡形成临空面,导致坡体失去支撑诱发滑坡。
(2)大气降水。滑坡区雨水充沛,大气降雨渗入山体斜坡内,导致斜坡土层遇水饱和,增大土体重度;且滑体浸泡软化易滑地层,使土层的抗剪强度大幅降低;同时水充满裂隙时形成静水压力,当出现水头差时形成动水压力从而诱发滑坡。
5 滑坡稳定性分析
5.1 稳定性计算方法
根据前文分析得知,该滑坡滑动面为土石接触界面,根据现场勘查及调查显示滑动面呈折线型,采用《滑坡防治工程勘查规范》(GB/T 32864—2016)所推荐的传递系数法进行稳定性计算,计算公式为
(1)
Rn={Wn[(1-ru)cosan-Asinan]-RDn}tanφn+cnLn
(2)
Tn=Wn(sinan+Acosan)+TDn
(3)
(4)
式中:Kf为滑坡稳定性系数;ψi为第i块段的剩余下滑力传递至第i+1块段时的传递系数(j=i),即φi=cos(ai-ai+1)-sin(ai-ai+1)tanφi+1;Ri为第i计算条块坡体抗滑力,kN/m;Ti为第i计算条块坡体下滑力,kN/m;ai为第i计算条块的滑面倾角(°);Ci为第i计算条块的内聚力,kPa;φi为第i计算条块内摩擦角,(°);Li为第i计算条块滑面长度,m;A为地震加速度;Wi为第i条滑块自重力;ru为孔隙压力比;RDi为第i条块渗透压力产生的垂直滑面分力;TDi为第i条块渗透压力产生的平行滑面分力;RDn为第n条块渗透压力产生的垂直滑面分力;TDn为第n条块渗透压力产生的平行滑面分力;Rn为第n计算条块坡体抗滑力,kN/m;Tn为第n计算条块坡体下滑力。
5.2 计算参数
滑动面抗剪强度根据反演分析及土工试验综合确定。结合滑坡变形特征,得知滑坡是在暴雨期间发生剧烈变形,滑坡在暴雨期间应处于极限稳定状态或临近极限稳定状态。因此,以能使得滑坡稳定性系数处于1.01左右的c、φ组合值作为分析参考对象,与室内土工试验结果分析对比,根据相关技术规范要求综合确定滑动面抗剪强度。反分析计算结果见表1。
表1 计算剖面滑坡稳定性反分析计算结果
室内土工试验统计的滑带土抗剪参数见表2。勘查得知滑带物质成分主要为含卵石粉质黏土,其粗颗粒含量33%左右,因此根据《滑坡防治工程勘查规范》(GB/T 32864—2016)第13.2.7条相关规定,须对试验得出的内摩擦角进行修正,增大系数取1.25,修正后的滑带土内摩擦角见表3。
根据表1可知,暴雨工况下,能使计算剖面滑坡稳定性系数处于1.01左右的c、φ组合值有3种。组合1:c=18 kPa,φ=6°,稳定性系数1.012;组合2:c=20 kPa,φ=5°,稳定性系数1.005;组合3:c=22 kPa,φ=4°,稳定性系数1.005。
天然状态:因《滑坡防治工程勘查规范》(GB/T 32864—2016)中规定了滑带土内摩擦角增大系数取值范围为1.15~1.25,然而实际勘查该滑坡滑带土粗颗粒含量在33%以上,本次修正取最大修正系数1.25,内摩擦角取修正后的标准值8°,黏聚力取范围值中的标准值,即c=37.3 kPa,φ=8°。
饱和状态:选择反演分析组合3,即c=22 kPa,φ=4°。因为该组合中φ=4°,与修正后的滑带土内摩擦角最大值4.1°最为接近。
最终选定的计算参数取值见表4。
5.3 计算工况及结果分析
5.3.1 计算工况
工况1:仅考虑滑坡的自重力,计算时采用滑带土天然状态下的c、φ值,采用坡体的天然容重。
工况2:根据水文地质资料,滑坡受地表水和地下水影响大,滑带土主要为松散卵石土,透水性强,在雨季,滑坡体大量充水,持续降雨可能使滑体全部饱水,计算时采用饱水容重。滑带土强度采用饱和状态下残余c、φ值。
工况3:同时考虑了暴雨和地震,计算时采用饱水容重和饱和状态下残余c、φ值,根据规范,设计基本地震加速度值取0.05g(g为重力加速度)。
5.3.2 计算结果
潜在滑坡在3种不同工况下的稳定性系数见表5。
表2 滑带土试验抗剪参数统计
表3 修正后的滑带土内摩擦角
表4 滑带土抗剪参数综合取值
表5 稳定性计算结果及评价
工况1条件下,稳定性系数大于1.2,滑坡处于稳定状态。
工况2条件下,稳定性系数小于1,说明滑坡整体处于不稳定-欠稳定状态,该工况是滑坡最为常见的工况,滑坡主要在该工况下发生变形。
工况3条件下,稳定性系数小于1,说明暴雨状态滑坡整体处于不稳定-欠稳定状态,在该工况下滑坡会发生整体变形。
6 结论
(1)滑坡体横向宽度约120 m,主滑方向180°,顺主滑方向长约95 m,滑体厚度1.5~6 m,体积约5.0×104m3,为小型推移式滑坡。
(2)滑坡形成的主要因素为人类工程活动破坏坡体平衡和大气降水。
(3)综合分析潜在剖面在不同工况下的稳定性系数,工况2、工况3条件下,潜在滑坡的稳定性系数位于0.93~0.96,说明潜在滑坡处于欠稳定状态,可能发生滑动变形。因此,大气降雨为展销大棚滑坡的主要诱发因素。