雨水入渗对非饱和土坡稳定性影响的参数研究

很多国家和地区的斜坡失稳与雨水入渗有密切关系。通过参数分析研究可以深化对这种关系的认识和理解,因而对滑坡灾害的预测和预防有重要意义。针对香港地区一种典型非饱和土斜坡,用有限元法模拟雨水入渗引起的暂态渗流场,然后将计算得到的暂态孔隙水压力分布用于斜坡的极限平衡分析。计算中采用延伸的摩尔－库伦破坏准则以便考虑基质吸力对抗剪强度的贡献,研究了降雨特征、水文地质条件及坡面防渗处理等因素对暂态渗流场和斜坡安全因数的影响。数值模拟结果表明：降雨强度、降雨历时和雨型对暂态渗流场及斜坡稳定性有明显的影响;土体的渗透系数,尤其是渗透系数各向异性的影响特别显著;斜坡中相对隔水层的存在以及斜坡防渗护面的效果等因素的影响均不容忽视。     

降雨过程中边坡临界滑动场

降雨入渗过程中孔隙水压力的升高与基质吸力的降低引起边坡稳定性的下降,是导致边坡滑塌的主要诱导因素。利用饱和-非饱和渗流有限元计算得到的孔隙水压力场,基于Fredlund提出的非饱和土抗剪强度理论,对边坡临界滑动场进行改进,提出可以考虑降雨过程的边坡临界滑动场数值模拟方法,能够方便、快速地计算出边坡局部、整体安全系数和相对应的临界滑动面在降雨过程中的变化历程。将该法用于一个典型均质边坡和一个非均质边坡在降雨过程中的稳定性计算,分析降雨持续时间、降雨强度和非饱和强度参数取值等因素对边坡稳定性的影响,并将计算结果与其他方法进行比较,结果表明临界滑动场方法能搜索任意形状最危险滑面,计算的安全系数合理。     

降雨条件下土坡饱和-非饱和渗流及稳定性分析

基于非饱和土力学理论,利用有限元法对降雨条件下边坡的饱和-非饱和渗流及稳定性进行了探讨。数值分析表明,雨水入渗使边坡非饱和区土体的基质吸力减小,是导致边坡稳定性降低的主要因素。进而考察了基质吸力的不同取值及裂隙的不同位置对边坡稳定性的影响,结果表明：在降雨初期基质吸力取值越大,得到的边坡安全系数越高,但经过长时间降雨即边坡土体基本达到饱和后,基质吸力的取值不同对边坡安全系数基本上没有影响;裂隙的存在为雨水入渗提供了良好的通道,因而导致边坡稳定性相对于无裂隙状态来说有所降低。降雨初期下游裂隙使边坡稳定性降低较快,但随着降雨历时的增加,即下游土体逐渐变为饱和状态后,边坡的稳定性变化则主要由上游裂隙控制。     

基于Green-Ampt模型的多层结构边坡降雨入渗改进计算方法及稳定性影响研究

降雨作用下,边坡土体的饱和度及含水率升高,基质吸力减小。随着降雨历时的增长,雨水入渗深度对坡体的稳定性产生影响。然而传统多层结构边坡的入渗计算方法并未考虑随入渗深度不断变化的基质吸力与层间积水点的形成,且忽略饱和层内沿坡体层面流动的部分雨水对入渗过程的影响,亦未考虑潜在滑动面位置随降雨历时的变化。将入渗过程分解为若干个子过程,并基于Green-Ampt（G-A）入渗模型对传统多层结构边坡的入渗计算方法进行改进,以对每个子过程进行求解,最后将其合并为整体入渗过程的解。在此基础上对层间积水点的形成时刻进行计算,进而分析雨水入渗深度与时间的关系,并研究降雨强度与雨水入渗深度对边坡不同位置处（湿润锋、饱和层）稳定系数和滑动面位置的影响。研究表明:（1）基于G-A模型的改进计算方法所得结果比传统多层结构边坡入渗计算方法所得结果更接近于数值模拟结果。（2）对于多层结构土质边坡,其安全系数随着雨水入渗深度的增加不断降低,并且在层间积水点形成时产生突变现象。（3）随着降雨历时和降雨强度的增大,边坡中潜在滑动面位置会产生变化,前期潜在滑动面位置出现在湿润锋处,后期则出现在饱和层交界面处。该方法提高了多层结构边坡传统降雨入渗计算方法的精度,更加全面的对多层结构边坡的稳定性进行评价,其工程应用范围亦得到进一步扩大     

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

Soil water balance model for precipitation-induced shallow landslides

Precipitation infiltration is one of the most significant triggering factors for slope failure occurrence in many places around the world. Knowledge of the mechanisms leading to precipitation-induced slope failures is of great importance to the management of landslide hazard. In this study, a soil water balance model is developed to estimate soil water flux during the process of infiltration from rainfall data, with consideration of storm periods and non-storm periods. Two important assumptions in this study are given: (1) instantaneous uniform distribution of the degree of saturation and (2) a linear relationship between evapotranspiration and the related degree of saturation. For storm periods, the Brooks and Corey model estimates both the soil water retention curve and soil water parameters. The infiltration partition is employed by an infinite-series solution of Philip in conjunction with the time compression approximation. For none-storm periods, evapotranspiration can be derived for the moisture depletion of soil water. This study presents a procedure for calculating the safety factor for an unsaturated slope suffering from precipitation infiltration. The process of infiltration into a slope due to rainfall and its effect on soil slope behavior are examined using modified Mohr–Coulomb failure criterion in conjunction with a soil water balance model. The results indicate that the matric suction, which is closely related to slope stability, is affected by the degree of saturation controlled by rainfall events.     

滑坡防治前后滑带土基质吸力特征研究

降雨使坡内地下水位上升,坡体的基质吸力和暂态孔隙水压力会随着降雨过程和时间呈现不同的变化趋势。采用抗滑桩、挡墙等工程防治后,针对滑坡体滑带土的性质变化以及抗滑结构的设置是否会影响到滑坡体的自然排水通道,导致坡体地下水位的抬升,从而影响滑坡防治效果等问题进行了研究。以鹰厦线K290滑坡为主线,通过开展滑坡饱和-非饱和渗流模型试验与数值分析,探讨了降雨入渗及地下水位变化下滑坡体及滑带土体积含水率与基质吸力的变化规律,以及其对滑坡防治前后坡体渗流场的影响。在此基础上,探讨了滑带土基质吸力对抗剪强度的贡献。研究表明：防治前后滑坡体及滑带土基质吸力受降雨强度等条件影响明显,不同深度处滑带土基质吸力变化呈现不同的变化规律,在土质滑坡的防治中应考虑抗滑桩的布置对滑带土性质的影响。     

降雨对均质各向异性土质边坡稳定性的影响

采用有限元法分析了非饱和均质各向异性土坡中的降雨渗透过程和土水相互作用的变化过程,使用修改的Mohr-Coulomb破坏准则考虑非饱和土的抗剪强度,根据基质吸力变化引起抗破坏强度变化的原理,寻找最危险的滑移面位置,计算最小稳定安全系数,研究了土坡遭受降雨渗透时的安全稳定问题。实例分析结果表明,该法在评判降雨时非饱和土坡的稳定性方面是合理可行的。同时分析数据表明裂隙的存在对土质边坡稳定性的影响是显著的,在对雨水入渗的土边坡稳定性分析时,必须考虑裂隙的影响。     

三峡库区黄土坡滑坡降雨入渗模型研究

传统的入渗模型未考虑坡角和降雨强度对滑坡入渗过程的影响,为了更好地描述黄土坡滑坡降雨入渗过程,在Green- Ampt入渗模型的基础上推导了考虑坡体倾角和小降雨强度影响的降雨入渗模型。为了获取改进的入渗模型参数,在黄土坡滑坡1#崩滑体上进行了双环渗透试验与降雨、土的含水率和基质吸力现场监测。结果表明,黄土坡滑坡1#崩滑体饱和渗透系数为4.81×10-5 m/s;降雨时体积含水率增加,降雨停止后体积含水率降低,深部表现出一定的滞后特性;基质吸力变化趋势与体积含水率相反,降雨使其减小,降雨停止后逐渐增大。通过双环渗透试验与现场监测,获取了黄土坡滑坡降雨入渗模型参数值。将入渗模型计算值与现场监测数据进行对比,该模型计算值与现场监测数据吻合,说明该降雨入渗模型可用于黄土坡滑坡降雨入渗分析。     

Rain-triggered slope failure of the railway embankment at Malda,India

The common slope stability analysis is incapable of accurately forecasting shallow slides where suction pressures play a critical role. This realization is used for elaborate stability analyses which include soil suction to better predict rainfall-induced slides at railway embankment at Malda where three known cases of slope failures and train derailments occurred after heavy rainfall. The relationship between the soil–water content and the matric suction is established for the embankment soil. It is then used in the coupled analyses of seepage and slope stability to estimate performances of the embankment at different intensity and duration of rainfall. The numerical simulations are performed with the FE code Geo-Studio. The numerical results show significant reduction in the factor of safety of the railway embankment with the increase in the intensity and duration of rainfall. The effectiveness of the proposed mitigation measures including placement of 2 m-wide free draining rockfill across the slopes and drilling 5-m-long sheet pile wall at the toe of the embankment is studied numerically. The study confirms that the proposed mitigation measures effectively increase the factor of safety of the embankment and stabilizing it even in case of a heavy rainfall of 25 mm/h over 12 h.     

Rainfall-induced landslides by deficit field matric suction in unsaturated soil slopes

Landslides are mainly triggered by decrease in the matric suction with deepening the wetting band by rainfall infiltrations. This paper reports rainfall-induced landslides in partially saturated soil slopes through a field study. A comprehensive analysis on Umyeonsan (Mt.) landslides in 2011 was highlighted. The incident involves the collapse of unsaturated soil slopes under extreme-rainfall event. Fundamental studies on the mechanism and the cause of landslides were carried out. A number of technical findings are of interest, including the failure mechanism of a depth of soil and effect of groundwater flow, the downward movement of wetting band and the increase of groundwater level. Based on this, an integrated analysis methodology for a rainfall-induced landslide is proposed in this paper that incorporates the field matric suction for obtaining hydraulic parameters of unsaturated soil. The field matric suction is shown to govern the rate of change in the water infiltration for the landslide analysis with respect to an antecedent rainfall. Special attention was given to a one-dimensional infiltration model to determine the wetting band depth in the absence of the field matric suction. The results indicate that landslide activities were primarily dependent on rainfall infiltration, soil properties, slope geometries, vegetation, and groundwater table positions. The proposed methodology has clearly demonstrated both shallow and deep-seated landslides and shows good agreement with the results of landslide investigations.     

不同渗流边界条件下粉砂边坡失稳模型试验研究

研究不同渗流边界条件下粉砂边坡的稳定特性对滑坡安全评价及防灾减灾具有重要意义。开发研制了室内边坡模型试验系统,分别进行了在上部边界入渗、侧向边界入渗和底部边界入渗的边坡失稳模型试验。分析了各试验中边坡内部不同位置体积含水率、孔隙水压力和基质吸力的变化规律以及坡体破坏过程。结果表明,不同入渗边界条件下边坡土体局部由非饱和变为饱和状态,基质吸力消失,最后在坡面形成不稳定区域塌落破坏。上部边界入渗、侧向边界入渗和底部边界入渗所引起的破坏模式根据其特点可分为局部浅层破坏、多级后退式破坏和坡面滑动破坏。通过监测边坡不同位置的体积含水率,可以为渗流引发的边坡失稳预警机制提供参考依据。     

生态护坡对边坡岩土体基质吸力的影响研究

对植被蒸腾作用以及阻滞地表径流作用进行了研究,提出生态护坡对边坡岩土体基质吸力的影响主要有3个方面：(1) 植被的蒸腾作用消耗边坡岩土体中的水分,使边坡体基质吸力升高,对边坡稳定有利;(2) 植被阻滞地表径流,增加降雨入渗量,使边坡体含水率增高,基质吸力降低,对边坡稳定不利;(3) 植被的存在究竟使边坡体基质吸力如何变化以及对边坡稳定是否有利,最终决定于降雨量和雨型。在此基础上,更进一步提出在降雨丰沛、暴雨较多的地区进行生态护坡时,应对边坡体采取防排水处理措施。     

An analytical solution for rainfall infiltration into an unsaturated infinite slope and its application to slope stability analysis

Surficial slope failures in residual soils are common in tropical and subtropical regions as a result of rainfall infiltration. This study develops an analytical solution for simulating rainfall infiltration into an infinite unsaturated soil slope. The analytical solution is based on the general partial differential equation for water flow through unsaturated soils. It can accept soil–water characteristic curve and unsaturated permeability function of the exponential form into account. Numerical simulations are conducted to verify the assumptions of the analytical solution and demonstrate that the proposed analytical solution is acceptable for the coarse soils with low air entry values. The pore‐water pressure (pwp) distributions obtained from the analytical solution can be incorporated into a limit equilibrium method to do infinite slope stability analysis for a rain‐induced shallow slip. The analysis method takes into account the influence of the water content change on unit weight and hence on factor of safety. A series of analytical parametric analyses have been performed using the developed model. The analyses indicate that when the residual soil slope, consisting of a completely decomposed granite layer underlain by a less permeable layer, is subjected to a continuous heavy rainfall, the loss of negative pwp and the reduction in factor of safety were found to be most significant for the shallow soil layer and during the first 12 h. The antecedent and subsequent rainfall intensity, depth of a less permeable layer and slope angle all have a significant influence on the pwp response and hence the slope stability. Copyright © 2012 John Wiley & Sons, Ltd.     

库水位下降过程中土坡稳定强度折减有限元分析

采用强度折减有限元法,考虑边坡土体的非饱和-非稳定渗流过程,从坡体内浸润线位置、土体的渗透系数、水位下降速率、下降比及基质吸力等5个方面研究了水位下降过程中岸坡的整体稳定性。结果表明,水位下降速率对高渗透性土坡内孔压的影响显著,对低渗透性土坡内孔压影响较小,对边坡安全系数的影响程度可达15%。     

降雨入渗条件下多级黄土高边坡稳定性分析

在西北地区进行基础工程建设,遇到了大量的多级黄土高边坡工程。边坡工程属于永久性工程,难免遇到降雨状况。降雨入渗会引起边坡土体饱和度与土体重度的变化,进而导致边坡体基质吸力与有效应力场的变化,使边坡的稳定性受到较大的影响。本文通过考虑边坡体不同深度处基质吸力随雨水入渗的变化情况,建立了渗流场与应力场的耦合计算模型,对降雨入渗条件下多级黄土高边坡的稳定性进行分析。最后,依托实际边坡工程,利用PLAXIS 3D岩土有限元软件,建立边坡稳定性计算模型,通过设置降雨量随时间的变化函数与渗流边界条件,进行降雨入渗条件下多级黄土高边坡的稳定性数值计算,并与理论分析结果进行了对比分析。根据数值计算结果可得到多级黄土高边坡在降雨入渗条件下的变形、基质吸力、有效应力、潜在滑移面以及安全系数的变化情况,从而对降雨入渗影响下多级黄土高边坡的稳定性做出分析。研究结果可为降雨入渗条件下多级黄土高边坡的稳定性的研究提供一定的指导作用。     

基于电导性能的原状土柱降雨入渗规律研究

我国东南沿海地质环境复杂脆弱,由台风暴雨引发的滑坡频频发生.雨水入渗到土体内可能导致坡体变形,进而产生滑坡、泥石流等灾害.以福建三明岩兜滑坡为研究对象,通过自行研制的人工降雨土柱入渗试验装置,对研究区滑坡的残积土柱进行了不同降雨强度下(20 mm·h-1、60 mm·h-1)的累次循环降雨实验,考虑了不同降雨量、降雨历...     

The role of infiltration processes in steep slope stability of pyroclastic granular soils: laboratory and numerical investigation

Rainfall-induced landslides on steep slopes are a common feature in much of Italy’s mountain areas covered by shallow-pyroclastic deposits. Generally, these deposits are unsaturated and have a slope angle higher than 40°–50°; hence their stability is due to the positive effect of matric suction on soil shear strength. During rainfall, rainwater infiltration causes a decrease in suction, which in turn causes changes in soil mechanical and hydraulic properties, leading towards an instability process. However, the response of pyroclastic soil slopes to rainwater infiltration is not fully understood. The aim of this study is to link slope instability to the infiltration process on the basis of advanced geotechnical characterization, in situ monitoring and numerical analysis calibrated through a back-analysis of well-instrumented flume tests.     

上方来水对坡面降雨入渗及土壤水分再分布的影响

在防止土壤侵蚀和雨后抑制蒸发的条件下,利用室内模拟降雨试验,研究了上方来水对坡面降雨入渗、湿润锋运移以及土壤水分再分布的影响。结果表明:对于初始含水量很低的土壤,与上方无来水相比,上方来水时降雨入渗过程中入渗率有一个上升的阶段,但平均入渗率反而降低;在降雨入渗初期,由于上方来水的沿程入渗,上方来水对坡面湿润锋运移的影响较大,但随后几乎没有影响,湿润锋的运移主要与基质势梯度有关;土壤水分沿坡面呈"波浪形"分布是坡面径流的波动性、上方来水(径流)的沿程入渗以及侧向沿坡向下流等综合作用的结果。     

Response parameters for characterization of infiltration

Rainwater infiltration is of particular interest with respect to slope stability hazard management. The process of rainfall infiltration into unsaturated soil is complex due to the number of soil parameters involved and the random nature of moisture flux boundary conditions. In this paper, the effect of rainfall intensity on infiltration is investigated through the use of finite element seepage modeling for a 5 m high soil column subjected to various rainfall intensities. The numerical modeling study identified three response parameters that can be used to describe one-dimensional infiltration into an unsaturated soil. The response parameters are the depth of wetting front, the matric suction reduction depth, and the sectional infiltration rate. The practical application of the rainfall response parameters in slope stability analysis is illustrated.