降雨入渗对非饱和土边坡稳定性的影响

孔郁斐　宋二祥　杨军　张龙英　施洪刚　刘剑

摘要：研究降雨入渗对边坡稳定性的影响规律。采用有限元法进行非饱和土边坡的二维非稳态渗流计算，考虑基质吸作用利用极限平衡法进行非饱和土边坡稳定安全系数计算，进而通过算例计算，分析了降雨过程中及降雨之后，边坡内孔隙水压分布、潜在滑裂面位置以及边坡稳定安全系数的变化情况。着重分析了降雨强度和降雨持续时间的影响，并特别注意分析降雨结束后的边坡稳定性。算例表明某些情况下边坡安全系数最小值出现在降雨之后的数小时或数天，而非降雨的过程中或降雨刚刚结束之时。

关键词：非饱和土；降雨；降雨入渗；土坡稳定

中图分类号：TU 433文献标志码：A文章编号：16744764（2013）06001606

降雨导致的滑坡是一种常见的自然灾害，在广东、广西、云贵川等地，降雨引至的滑坡时有发生，社会影响重大。近年来我国频繁遭遇强降雨等极端天气，由此引发的滑坡灾害越来越为人们所关注。所以分析降雨影响下的边坡稳定性是很值得重视的课题。

通常地下水位以上的土是非饱和土，它具有异于饱和土的特殊性质——降雨时边坡中非饱和土的饱和度逐渐上升，与此同时土的强度降低、渗透性增加，边坡变得更加危险。从20世纪开始，国际上便已展开针对非饱和土及非饱和土边坡的分析研究[1]。已有研究涵盖了很多方面，如非饱和土性质、稳态渗流条件下边坡稳定性分析、降雨过程中坡体内基质吸力分布及稳定性的变化等。研究方法包括足尺模型试验，缩尺模型离心机试验、数值模拟等。实验方面， Moriwaki等[2]、林鸿州等[3]都在足尺实验中实现了边坡的降雨型破坏，实验结果与真实情况很接近；张建民等[4]、Ling等[5] 分别进行过含软弱层边坡和长边坡的离心机实验，利用人工降雨让边坡破坏，发现了特殊边坡的破坏模式和规律。数值计算方面，也有很多研究者做了很多工作，例如Ng等[6]研究了降雨量、降雨持时对边坡稳定性的影响，通过计算发现当总降雨量一定时，降雨持续时间会对降雨结束时的安全系数有显著影响，且存在一个临界持时，当持时等于临界值时土坡稳定安全系数最小。Tsaparas等[7]利用假想的边坡分析了降雨量、降雨持时、初始水位、渗透性对降雨诱发滑坡的影响，发现上述参数的取值对计算结果影响较大。进行这类问题的数值计算，可采用有限元法或简化方法进行非稳态渗流计算[89]，用极限平衡法或强度折减有限元法进行边坡稳分析[1011]。孔郁斐，等：降雨入渗对非饱和土边坡稳定性的影响

前人的研究和计算已经覆盖了非饱和土边坡稳定性分析的很多方面，也曾研究过降雨过程中边坡稳定性的变化[12]，不过鲜有人关注雨后的情况。本文不仅考虑了降雨过程中，还分析了降雨结束后边坡稳定性的变化，展示了一个全过程。通过一个典型边坡的计算发现：虽然降雨过程中边坡的安全系数逐渐减小，但是雨停的时刻并不一定是安全系数最低的时刻，雨后随着雨水继续下渗，边坡安全系数还有可能在较长时间内继续减小。1非饱和土理论和边坡稳定分析方法

1.1非饱和土基质吸力

土是三相体，土中水压和气压分别为uw和ua。非饱和土中，uw低于ua，ua-uw被称作基质吸力，记作Ψ（单位：kPa）。一般若空气相互连通，可认为ua=0，Ψ=-uw。此时孔隙水压为负，基质吸力为孔隙水压的绝对值。在饱和土中ua=uw≥0，Ψ=ua-uw=0。后文中为描述方便，在同时涉及饱和区和非饱和区时使用“孔隙水压”一词，非饱和区的基质吸力等于孔隙水压的绝对值；强调非饱和区的变化特征时使用“基质吸力”一词，对应的孔隙水压等于基质吸力的相反数。

基质吸力分布或含水量分布是计算的初始条件，天然土层的吸力分布存在较大的随机性[14]，它与土的性质、周边环境以及水文历史有关，可以通过实地测量或经验估计来获得[15]，也可根据长期的气象条件，利用软件进行模拟给出其估计。易知，若上表面降雨量与蒸发量均为0，地下水面以上基质吸力将呈直线分布；但自然状态下土的基质吸力一般有上限，霜露、植被等自然环境的影响可以使得在相当的深度范围内基质吸力为常数，故通常可认为基质吸力在水面以上一定范围呈直线分布，某点的基质吸力水头等于该点到水面的距离，到达上限后，基质吸力与高度无关、为常数。

计算中边界条件的选取与地质情况、降雨量等有关。如果雨水能够全部从上表面入渗，不产生地表径流，则入渗的速度等于降雨速率，计算中上表面使用流量边界条件（第2类边界条件）；如果因降雨量较大或土体渗透性较小，雨水不能全部入渗，则可能产生地表径流，流过的区域表面几乎处于饱和状态，使用水头边界条件（第1类边界条件），令表面压力水头为0。除上述2种情况外，降雨时还可能出现更为复杂的情况，可酌情适当偏于保守地简化处理，或者进行考虑地表径流的更复杂的计算。

本文的渗流分析使用数值计算软件GeoStudio 2007的地下水渗流分析模块SEEP/W，计算1个均质的边坡，不考虑地表径流，用有限元法计算渗流过程中的瞬态孔压，进而得出不同时刻的孔隙水压力分布或基质吸力分布。

1.4边坡稳定分析

基于渗流计算结果及土体强度参数可以进行极限平衡法或强度折减有限元法的边坡稳定性计算。极限平衡法将滑体划分为较多土条，通过静力分析判断边坡的稳定性，计算滑动力与抗滑力的比值、滑动力矩与抗滑力矩的比值得到安全系数。常用的极限平衡法包括：瑞典条分法、简化毕肖甫法、简布法、MorgensternPrice法（简称MP法）等等。本文计算中使用MP法，它是一种严格条分法，同时考虑了力的平衡和力矩的平衡；在建立极限平衡方程时，同时考虑土条间法向力和切向力、并建立反映二者关系的条间力函数X=Eλf（x）。力安全系数与力矩安全系数随着条间力函数权重λ的变化而变化，MP法的安全系数位于这两条曲线的交点，该方法适用于计算包括圆弧滑裂面在内的任意滑裂面。2边坡在降雨全过程中的状态变化

3影响规律研究

3.1不同深度的滑裂面，受降雨影响程度不同

在极限平衡法计算安全系数时，可以假设无数个滑裂面（圆弧的或者任意形状的），在一定限制条件下可以找到安全系数最小的滑裂面。由于土体的性质随时间改变，每个滑裂面的安全系数都在随时间变化。含水率的增大导致土体抗剪强度减小——对较深的滑裂面，若降雨尚未渗到滑裂面所处深度，其安全系数随降雨持时的变化不大；而对较浅的滑裂面，雨水很快可以渗入它通过的区域，其安全系数在降雨开始后很快就降低。这里所说的“深”和“浅”并非几何概念，而是与受降雨影响的程度有关。边坡土体不透水的情况下，即便降雨持续很久，坡面以下几十cm处含水率依旧保持不变；如果边坡渗透系数较大，坡内孔隙水压对降雨的响应就较快。

通过一些计算，可以总结出深层滑裂面和浅层滑裂面的安全系数变化规律。图11是两类滑裂面安全系数随时间变化的示意图。降雨开始后浅层滑裂面安全系数逐渐减小，如果降雨时间较短、降雨量小，则雨水可以很快下渗，土坡的安全系数在雨停后即可回升；如果降雨时间较长，安全系数在降雨过程中会趋于一个下限值，雨停后安全系数回升。深层滑裂面在降雨初期不会受影响，如果降雨时间较短，雨停之前该处土体还不受影响，水在雨停后继续下渗，安全系数才逐渐减小、继而回升；如果降雨时间较长，安全系数在雨停之前就开始减小。

3.2降雨持续时间的影响

降雨的强度和持续时间对边坡稳定性的变化趋势也有较大的影响。如果降雨速率不变，持续时间延长，则较长时间后边坡内的渗流会趋于稳态，安全系数也达到一个固定值：新的算例如图12所示，降雨时间由基本工况下的1 d延长为3 d，降雨速率不变，2 d后安全系数趋近于1.78，并保持恒定，这是11.12 mm/h的降雨速率下本算例边坡可能达到的最低安全系数。第72 h降雨停止后安全系数立刻反弹，不再有下降段。通过该算例和其他相关算例的计算，我们发现，在特定降雨强度下，可找到一个与之对应的临界持续时间，在此临界时间之前安全系数几乎呈线性变化，该时间点以后如果降雨速率不变，则安全系数不再继续减小。

4结论

研究了降雨对均质非饱和土边坡稳定的影响的相关理论和计算方法，结合算例计算了不同工况下降雨过程中及降雨结束后边坡稳定性变化，可以得到以下结论：

1）降雨强度和持续时间对边坡安全系数有显著影响。若降雨速率不变、降雨时间足够长，安全系数会在一段时间内持续下降，其后趋于恒定。

2）安全系数不一定在雨停的时刻达到最低点。雨停后随着雨水继续下渗，边坡安全系数还有可能在数小时至数天内继续减小，正因如此，很多滑坡灾害发生在雨后。

3）浅层滑裂面对降雨的响应较快，其安全系数在降雨开始后很快开始下降，降雨结束后也很快开始回升；深层滑裂面的安全系数在更长时间后才会对降雨情况的变化有所响应。

参考文献：

[1]Fredlund D G，Rahardjo H. Soil mechanics for unsaturated soil [M]. Beijing： China Architecture and Building Press， 1997.

[2]Moriwaki H， Inokuchi T， Hattanji T， et al. Failure processes in a fullscale landslide experiment using a rainfall simulator[J]. Landslides， 2004， 1（4）： 277288.

[3]林鸿州，于玉贞，李广信，等. 降雨特性对土质边坡失稳的影响[J]. 岩石力学与工程学报，2009， 28（1）： 198204.

[4]张建民，王睿，张嘎. 降雨条件下含软弱夹层土坡的离心模型试验研究[J]. 岩土工程学报，2010（10）： 15821587.

[5]Ling H， Ling H I. Centrifuge Model Simulations of RainfallInduced Slope Instability [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2012， 138（9）： 11511157.

[6]Ng C W W， Shi Q. Influence of rainfall intensity and duration on slope stability in unsaturated soils [J]. Quarterly Journal of Engineering Geology， 1998， 31（Part 2）： 105113.

[7]Tsaparas I， Rahardjo H， Toll D G， et al. Controlling parameters for rainfallinduced landslides [J]. Computers and Geotechnics， 2002， 29（1）： 127.

[8]Santoso A M， Phoon K K， Quek S T. Effects of soil spatial variability on rainfallinduced landslides [J]. Computers & Structures， 2011， 89（11/12）： 893900.

[9]孔郁斐. 降雨对非饱和土边坡稳定性影响研究[D]. 北京： 清华大学， 2012.

[10]于玉贞，林鸿州，李荣建，等. 非稳定渗流条件下非饱和土边坡稳定分析[J]. 岩土力学，2008， 29（11）： 28922898.

[11]宋二祥. 土工结构完全系数的有限元计算[J]. 岩土工程学报，1997， 19（2）： 17.

[12]Zhao H， Song E X. A method for predicting coseismic displacements of slopes for landslide hazard zonation [J]. Soil Dynamics and Earthquake Engineering， 2012， 40： 6277.

[13]Nuth M， Laloui L. Effective stress concept in unsaturated soils： Clarification and validation of a unified framework [J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2008， 32（7）： 771801.

[14]汤明高，许强，黄润秋，等. 滑坡体基质吸力的观测试验及变化特征分析[J]. 岩石力学与工程学报，2006， 25（2）： 355362.

[15]吴俊杰，王成华，李广信. 非饱和土基质吸力对边坡稳定的影响[J]. 岩土力学，2004， 25（5）： 732736， 744.

（编辑胡玲）