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

降雨入渗对非饱和土边坡稳定性影响 　　摘要：研究降雨入渗对边坡稳定性的影响规律。采用有限元法进行非饱和土边坡的二维非稳态渗流计算，考虑基质吸作用利用极限平衡法进行非饱和土边坡稳定安全系数计算，进而通过算例计算，分析了降雨过程中及降雨之后，边坡内孔隙水压分布、潜在滑裂面位置以及边坡稳定安全系数的变化情况。着重分析了降雨强度和降雨持续时间的影响，并特别注意分析降雨结束后的边坡稳定性。算例表明某些情况下边坡安全系数最小值出现在降雨之后的数小时或数天，而非降雨的过程中或降雨刚刚结束之时。 　　关键词：非饱和土；降雨；降雨入渗；土坡稳定 　　中图分类号：TU 433文献标志码：A文章编号2013 　　降雨导致的滑坡是一种常见的自然灾害，在广东、广西、云贵川等地，降雨引至的滑坡时有发生，社会影响重大。近年来我国频繁遭遇强降雨等极端天气，由此引发的滑坡灾害越来越为人们所关注。所以分析降雨影响下的边坡稳定性是很值得重视的课题。 　　通常地下水位以上的土是非饱和土，它具有异于饱和土的特殊性质——降雨时边坡中非饱和土的饱和度逐渐上升，与此同时土的强度降低、渗透性增加，边坡变得更加危险。从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种情况外，降雨时还可能出现更为复杂的情况，可酌情适当偏于保守地简化处理，或者进行考虑地表径流的更复杂的计算。 　　本文的渗