持续降雨致我县公路多处山体滑坡过往司机请慢行
摘 要：滑坡作为最常见的地质灾害之一，其造成自然环境和社会环境的损失是毋庸置疑的，在渝东北地区，滑坡灾害影响了人们的正常生活，制约着当地的经济发展。强降雨作用通常是触发大型滑坡的主要诱因，为更好地研究持续降雨条件下滑坡的变形破坏机理，通过ABAQUS软件模拟该地区典型滑坡在连续8天降雨条件下的变形破坏，探究该地区降雨诱发滑坡的致灾机理，得出在渝东北地区长期持续的强降雨入渗是滑坡诱发的主要因素，其斜坡变形模式属于孔隙水压力诱发的平推式地质模式。研究成果对推动该地区滑坡预警预报，滑坡治理具有积极意义。
关键词：滑坡；降雨；ABAQUS；数值模拟
引言
大量事实和研究表明，降雨充沛且多暴雨的气候是渝东北地区滑坡灾害多发的主要触发因素。重庆自古就有十滑九雨的论断，国内外学者通过大量研究表明，滑坡的活动强度与降雨量呈正比关系。渝东北地区历年发生的极端降雨，均诱发了大范围区域性滑坡灾害。近年来，国内外学者对降雨诱发滑坡进行了大量的研究与实践，本文通过ABAQUS软件模拟该地区典型滑坡在持续降雨条件下的变形破坏，对探究该地区降雨诱发滑坡的致灾机理，推动该地区滑坡预警预报，滑坡治理具有积极意义。
1 工程实例
巫山县龙溪镇汪家坪滑坡位于龙溪镇双河村二组，整体形态成“扇”型，滑坡横宽约760.0m，纵长约620.0m，面积约24.7×104m2，厚度17.7m～44.8m，平均厚度35.0m，体积约864.5×104m3，属大型深层土质滑坡，主滑方向145°～157°。滑坡后缘以坡顶平台外侧变形裂缝为界，滑坡前缘剪出口位于长溪河河床左岸，滑坡东北侧以凸起山脊基岩出露为界，东南侧以地形转折处为界，西北侧以地形转折为界，西南侧以地形出现转折，且基岩出露部位为界。滑坡体局部区域变形较强烈，滑坡后缘发育一条长约169.6m的下挫裂缝，滑坡中前部发育多条贯通拉张下挫裂缝，并出现多处局部滑塌现象。
假设汪家坪滑坡滑体、滑床属于均值各相异性。坡体的有效黏聚力c=45KPa、有效内摩擦角φ=28°、重度γ=1.85g/cm3、弹性模量E=30MPa、泊松比υ=0.3、初始空隙比e=0.8、饱和渗透系数Kws=0.019（m/h）
一般来说，水力渗透系数在饱和土中是常量，而在非饱和土中则与基质吸力或体积含水量有关。在本文分析中，土中水的渗透系数引用了文献[1]的计算公式，定义了渗透系数和基质吸力之间的关系：
同样需要定义土-水特征曲线，采用文献[2]的计算公式，定义饱和度与基质吸力之间的关系：
为研究降雨强度及渗透系数对边坡稳定性的影响，根据气象部门搜集到的巫山县龙溪镇双河村气象资料，整理可以得到降雨工况，连续8天降雨量分别为15mm、46mm、49mm、18mm、5mm、74mm、251mm、18mm。
本文的计算是在ABAQUS软件平台上进行二次开发来进行的，将土体视为孔隙材料，假设孔隙材料是具有固体架构中均布孔隙的多孔介质，同时将节点位移和孔隙水压力作为节点自由度。通过在ABAQUS中编写用户子程序USDFLD，引入修正的Mohr-Coulomb破坏准则，进行降雨入渗下非饱和土坡渗流场和应力场的耦合分析。在ABAQUS中把控制剖面进行单元格划分，共21718个有限元单元，同时在滑坡剖面上设置4个监控单元A，B，C，D，有限元计算模型如图1所示。
其设定的边坡初始饱和度状况为坡顶8%左右，坡面介于8%～100%之间；初始地应力采用自重作用下的地应力场，侧压力系数采用经验公式K0=1-sinφ；初始孔隙比为1.0。
2 模拟结果分析
边坡位移的變化情况降雨入渗使得滑坡表层的非饱和区含水量增高，基质吸力下降，边坡表层形成一个暂态饱和区，导致边坡表层产生与基质吸力变化有关的应变，最终边坡产生位移。降雨入渗引起边坡岩土体的位移变化包括两部分：（1）边坡的孔隙水压力上升，基质吸力减小，边坡岩土体强度下降，引起边坡体的位移变化；（2）边坡的孔隙水压力上升，导致边坡的非饱和岩土体产生应变。
经过持续8天的强降雨，该滑坡体的孔隙水压力等值线云图经过调整达到相对稳定状态，与此同时，滑坡体应变状态随着降雨的持续，调整并达到平衡。
数值模拟结果表明，历时持续8天的强降雨，最大日降雨量251mm/h，通过提取A、B、C、D四点的特征点位移曲线，滑坡体表面发生10cm-30cm不等的水平位移量。根据重庆市地质灾害防治工程勘查设计院编制的重庆市巫山县龙溪镇汪家坪滑坡应急抢险勘查报告，滑坡变形破坏特征如下，裂缝L1发育滑坡后缘，裂缝延伸长度169.6m，宽1.0～8.3cm，局部下挫20cm，延伸方向约70°～75°，局部段分为两条平行裂缝。裂缝发育L2滑坡中前部，裂缝延伸长度317.0m，宽1.0～11.5cm，延伸方向约50°～60°，局部段分为多平行裂缝。裂缝L3发育滑坡前缘，裂缝延伸长度280.6m，宽1.0～10.4cm，局部下挫45cm，延伸方向53°～90°。
这一计算结果与那一阶段实际发生的滑坡现象相当吻合，从而很好验证了本文提出的降雨条件下边坡稳定性分析模型和相应的有限元计算程序的可靠性，可以用于降雨条件下边坡稳定性评价分析。
3 结束语
进一步对变形机理分析，得出如下结论：
（1）受持续强降雨及特大暴雨影响，滑坡存在变形破坏的可能性。
（2）持续强降雨及特大暴雨直接诱发汪家坪滑坡再次出现强变形迹象，滑坡后缘出现一组贯通长度达169.6m的拉张裂缝和局部滑塌体，滑坡前缘出现多处密集贯通拉张、下挫变形裂缝。滑坡前缘变形裂缝越靠前张开度越宽，下挫高度越高，变形迹象越明显，变形程度越强烈。
（3）汪家坪滑坡为一老滑坡发生整体滑动后形成的滑坡堆积体，滑坡体结构较为松散，抗剪力学强度较低，自身稳定性安全储备不够。且由于滑坡区地下水位较为丰富，形成的巨大动水压力进一步降低滑坡体的稳定性。受持续强降雨及特大暴雨影响，大量大气降雨入渗增加滑坡体的自重，降低滑体抗剪强度，抬升地下水位增加其动水压力，从而降低滑坡体整体或局部稳定性，直接诱发滑坡后缘、滑坡中部左侧及滑坡前缘长溪河河谷左岸三级较陡斜坡区域出现裂缝、局部滑塌等变形迹象，为老滑坡的局部失稳。
（4）长期持续的强降雨入渗是其诱发的主要因素，其斜坡变形模式属于孔隙水压力诱发的平推式地质模式。
参考文献
[1]Cho S E，Lee S R. Instability of unsaturated soil slopes due to infiltration[J].Computers and Geotechnics，2001，28（3）：185-208.
[2]谭万沛，罗晓梅，王成华.暴雨泥石流预报程式[J].自然灾害学报，2000（3）：106-111.
作者简介：罗晓龙（1992-），男，重庆南川人，硕士，主要从事滑坡、边坡方面的研究工作。
