降雨入渗影响下边坡中的非饱和渗流特性

为了分析降雨入渗影响下非饱和土坡渗流特性,利用自制降雨模拟系统和实时监测系统,对降雨诱发非饱和土坡失稳过程进行全方位、多参量的实时监测,研究不同降雨条件下,不同坡度、不同压实度边坡坡体不同位置雨水入渗率和湿润峰的实变规律.结果表明:降雨入渗条件下,陡坡和高压实度土体不利于雨水入渗,而缓坡和低密实度土体入渗率变化快;实际土体吸力和含水量实时变化规律不同步,提出试验湿润峰概念,含水率（吸力）湿润峰点可按含水率（吸力）实时曲线的过渡区和雨后残余含水率（吸力）的线性交叉点确定;考虑单向吸湿或脱湿路径下土体含水率和吸力具有唯一对应关系,含水率湿润峰点与吸力湿润峰点的绝对值时差即为形成湿润峰所需时间;对比湿润峰实测值与Lumb半经验值散点分布规律,基于Lumb湿润峰深度计算公式提出非线性修正表达式.      

Investigation of wetting front behavior due to rainfall and ponding head effects in arid region wadi alluvium

Intensive experimental as well as numerical investigations have been carried out in order to develop empirical relationships to predict wetting front travel time and infiltration rate in deep arid region alluvium due to the effect of variable rainfall intensity and ponding head, respectively. A rainfall simulator that is capable of producing variable rainfall intensity up to 160 mm/h and a large diameter infiltrometer were constructed in the study area and then used to test the effect of rainfall and ponding head up to a depth of 3 m. Results of these experiments were utilized to calibrate an unsaturated zone model which was then used to test different scenarios for wetting front movement at depths beyond the scope of these experiments. Empirical infiltration relationships for rainfall and ponding head, respectively, were derived for shallow and deep soil columns from these experiments.     

基于数值模拟的湿润锋前进法测量精度分析

非饱和土的渗透系数函数测量难度大、耗时长,湿润锋前法(WFM)可在短时间内测得渗透系数函数(HF),但依赖肉眼识别的湿润锋,存在一定的局限性,且该方法的测量精度尚不明确,有待验证。文中针对湿润锋前进法进行研究,探讨初始含水率、湿润锋阀值、降雨入渗速率、传感器位置等因素对湿润锋前进法的测量精度的影响。采用Seep/W软件模拟均质土柱的入渗过程,分析湿润锋前进法数据,计算土体的渗透系数,将其和输入的渗透系数（可认为是真实解）进行比较,对湿润锋前进法的计算精度进行评估,并讨论误差的来源。分析结果表明,湿润锋前进法能够获得比较精确的计算结果;使用湿润锋特征含水率计算湿润锋前进速率,突破了原始湿润锋前进法存在“肉眼观察”的局限性,大大拓展了该方法的适用性;传感器间距对湿润锋前进法渗透系数函数计算精度没有直接影响;初始含水率越低,降雨入渗速率越大,渗透系数函数范围跨度越大。基于文中的分析结果,对湿润锋前进法的试验设计时,建议采用30～50 cm的土柱进行试验,传感器的数量建议为3～4,可以采用任意初始含水率进行试验接近干燥更好,为避免表面积水,建议降雨强度小于饱和渗透系数。     

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

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

Rainfall-triggered landslides in an unsaturated soil: a laboratory flume study

Extreme and/or prolonged rainfall events frequently cause landslides in many parts of the world. In this study, infiltration of rainfall into an unsaturated soil slope and triggering of landslides is studied through laboratory model (flume) tests, with the goal of obtaining the triggering rainfall intensity–duration (I–D) threshold. Flume tests with fine sand at two different relative densities (34 and 48%) and at slope angle of 56.5° are prepared, and rainfall (intensity in the range of 18 to 64 mm/h) is applied via a mist sprinkler system to trigger landslides. Soil water characteristic curve and hydraulic conductivity function of the fine sand are also presented. In flume tests, suction in the soil is measured with tensiometers, the progress of wetting front with time and deformations in the soil are also measured. Some of the findings of this study are: for the fine sand used in this study (a) the failure mechanism is infinite-slope type (mostly translational), and the failure surface is generally coincident with the wetting front or is in its vicinity, (b) the deformations leading to a landslide occurred abruptly, (c) both relatively high-intensity–short-duration rainfalls and relatively low-intensity–long duration rainfalls triggered landslides, (d) the shape of the I–D threshold is demonstrated to be a bilinear relation in log intensity–log duration plot, (e) below a certain rainfall intensity landslides are not triggered, (f) the effect of relative density of the soil on the I–D threshold is demonstrated by physical laboratory tests (as the relative density of the soil increases, the triggering rainfall intensity–duration threshold moves to larger rainfall events). The results of this study could be useful for accurate numerical modeling of rainfall-triggered landslides.     

降雨入渗下非饱和坡残积土湿润锋运移试验研究

我国东南沿海地区丘陵山地发育，降雨入渗到土体中的水分是导致滑坡灾害频发的关键因素。以福建省泉州市德化县地质灾害点为主要研究对象，考察典型地质灾害点具有代表性的坡积土与花岗岩残积土的渗透特性。利用自行研制的土体入渗装置，分别在降雨强度为15、30、60 mm/h条件下，考虑降雨历时一致（180 min）与过程降雨量一致（90 mm）两种工况开展一维土柱入渗试验，得到相应的各个土柱含水率、湿润锋、入渗率随时间变化的响应规律。试验结果表明：（1）土体渗透系数越大、雨强越大，土体湿润蔓延距离越深、速度越快。（2）降雨入渗过程中，土体含水率由浅及深逐次对降雨进行响应。不同雨强对含水率的影响主要体现在第一次响应时间以及饱和速度上，雨强越大，响应时间越快，饱和速度也越快。（3）提出可表征在不同雨强作用下，德化县马坪滑坡与崩土岭滑坡的湿润锋入渗公式。该研究成果对台风暴雨型滑坡的孕灾机制分析以及精细化监测预警具有重要的理论及实际意义。     

考虑优势流作用的降雨入渗边坡可靠度分析

降雨诱发滑坡是世界上最普遍的地质灾害,而关于降雨入渗对边坡的可靠度分析,优势入渗的影响被长期忽视。使用Comsol Multiphysics对降雨条件下优势入渗做数值求解,利用无限边坡模型计算边坡安全系数;并应用改进Cholesky分解法生成空间相关随机场,使用蒙特卡洛方法分析降雨过程中边坡的可靠度。结合确定性与可靠度计算对比均质入渗与优势入渗在降雨过程中边坡安全性变化：（1）降雨强度较低时,优势入渗安全性较好,而降雨强度较高时,均质入渗更稳定;（2）均质入渗中参数空间变异性是边坡失稳破坏的关键因素,而优势入渗的边坡失稳则由湿润峰快速推进所导致;（3）针对优势入渗模型研究,发现基质域与优势域水力交换强度较大时边坡有更大概率失稳,而较小的水力交换强度可能影响边坡底部的失稳破坏。     

压实黄土水分入渗规律及渗透性试验研究

为了研究压实黄土中的水分垂直入渗规律和非饱和渗透系数函数,在实验室内利用一维土柱垂直入渗模型试验装置,对两组压实黄土土柱试样分别进行了常水头入渗和降雨入渗试验。得到主要结论如下:（1）常水头入渗试验中,累积入渗量和湿润锋前进距离都随入渗时间呈幂函数形式增长,累积入渗量和湿润锋前进距离之间存在线性关系。入渗率在入渗初期最大,之后随入渗时间而快速降低,并在土柱试样底部出水以后达到稳定,且与湿润锋前进距离呈反比关系。（2）降雨入渗试验中,得到两组试样入渗过程中土-水特征曲线数据,分别用van Genuchten模型和Fredlund-Xing模型对两组试样进行了特征曲线拟合。并利用瞬时剖面法处理了入渗过程中水分和水势传感器的监测数据,得到两组试样的非饱和渗透系数,并拟合得到非饱和渗透系数与体积含水率之间的指数函数关系式。同时,采用van Genuchten和Fredlund等渗透系数模型分别对两组试样的非饱和渗透系数进行预测,通过对比模型预测结果和瞬时剖面法实测值,发现van Genuchten渗透系数模型预测结果更接近实测值。     

强降雨作用下堆积碎石土渗流规律研究

强降雨致堆(残坡)积碎石土渗透特性改变是触发滑坡、泥石流等降雨诱发型灾害的关键因素。为探究降雨入渗诱发滑坡的动态机制, 通过自制一维渗透仪与二维边坡降雨入渗模型, 对水竹湾隧道旁具有级配代表性的碎石土进行强降雨入渗试验, 并通过Seep/w软件进一步分析堆积碎石土湿润锋迁移入渗规律。研究结果表明:一维入渗深度-时间曲线存在阶梯状趋势, 即湿润锋迁移速度首先随表层土含水率的突增迅速增加至最大值, 然后缓慢减小, 最后湿润锋迁移速度略微增大; 湿润锋平均迁移速率随雨强的增大趋近线性增长; 碎石土沉降变形在入渗过程中经历缓慢增长、急剧增大、逐渐停滞3个阶段; 坡顶、坡腰及坡脚处入渗速率与一维入渗速率分别满足指数关系、对数关系及线性关系; 雨水最终入渗深度由大到小依次为坡脚、坡顶、坡腰。     

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.     

斜坡降雨入渗的改进Mein-Larson模型

土壤初始含水率分布是影响湿润锋下移的重要因素之一.传统的降雨入渗模型仅假定初始含水率均匀分布,并没有考虑非均匀分布的情况.以Mein-Larson降雨入渗模型为基础,假定初始含水率随垂直于坡面方向深度呈线性分布,推导了一种新的滑坡降雨入渗函数,弥补了原模型只能用于初始含水率均匀分布情况下的不足.研究结果表明:初始含水率分布对湿润锋下移有较大影响,在初始含水率呈线性分布的情况下,坡体表面含水率越大,含水率随垂直于坡面方向深度的变化率k值越接近零,湿润锋下渗速度越快,当k=0时,新模型退化为Mein-Larson降雨入渗模型,证明Mein-Larson降雨入渗模型是新模型的一个特例.与有限元法得到的结果对比表明,提出的降雨入渗新模型计算的湿润锋下渗深度与数值解结果相接近,证明了该方法的可靠性.      

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.     

Probabilistic failure analysis of infinite slopes under random rainfall processes and spatially variable soil

Rainfall-induced landslides occur during or immediately after rainfall events in which the pore water pressure builds up, leading to shallow slope failure. Thereby, low permeability layers result in high gradients in pore water pressure. The spatial variability of the soil permeability influences the probability such low permeability layers, and hence the probability of slope failure. In this paper, we investigate the influence of the vertical variability of soil permeability on the slope reliability, accounting for the randomness of rainfall processes. We model the saturated hydraulic conductivity of the soil with a one-dimensional random field. The random rainfall events are characterised by their duration and intensity and are modelled through self-similar random processes. The transient infiltration process is represented by Richards equation, which is evaluated numerically. The reliability analysis of the infinite slope is based on the factor of safety concept for evaluating slope stability. To cope with the large number of random variables arising from the discretization of the random field and the rainfall process, we evaluate the slope reliability through Subset Simulation, which is an adaptive Monte Carlo method known to be especially efficient for reliability analysis of such high-dimensional problems. Numerical investigations show higher probability of slope failure with increased spatial variability of the saturated hydraulic conductivity and with more uniform rainfall patterns.     

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.     

强降雨条件下基岩型层状边坡入渗模型及稳定性研究

强降雨作用下基岩型层状边坡易发生失稳破坏,给人们生命财产造成重大损失。为探索基岩型层状边坡的降雨入渗过程,以Green-Ampt入渗模型为基础,考虑了边坡几何特征及饱和带渗流作用,得到了边坡不同入渗阶段入渗率和湿润峰深度的计算式,建立了适用于基岩型层状边坡的降雨入渗计算模型;在分析基岩型层状边坡稳定性时,考虑了饱和带的渗流作用,联合降雨入渗模型与极限平衡法,分析了基岩面与湿润峰面稳定性的变化规律,得到了计算基岩型层状边坡安全系数的解析表达式;结果表明,该模型的计算与试验结果具有一致性。利用该计算模型与传统计算方法对张家湾滑坡进行了强度为30 mm/h的降雨入渗与稳定性分析,得到了不同计算方法下张家湾滑坡湿润峰深度及安全系数随降雨历时增加而发生的变化规律,结果表明该计算模型在分析基岩型层状边坡稳定性时优于传统分析方法。     

求取现场黄土饱和渗透系数的双环入渗法

土体渗透性参数的测定是解决与水有关的岩土工程问题的关键所在。为使现场原状土体饱和渗透系数的测定更为准确,通过一种易于安装和固定的防蒸发型双环入渗仪,研究黄土的入渗规律与特点,试验结果表明:黄土的入渗分区可分为饱和区、湿润区与干土区,其中饱和区的入渗深度约占总入渗深度的1/2,而湿润区土体的饱和度在70%以上,这一特征与Green-Ampt入渗模型的假设较为接近,但在应用基于Green-Ampt模型的规范法（SL237-042-1999）求取黄土的饱和渗透系数时,发现该法会高估黄土的饱和渗透系数。因此,结合Green-Ampt入渗模型与土水特征曲线主要增湿路径的特点,提出了能合理测定现场黄土饱和渗透系数的双环入渗法,对Green-Ampt入渗模型参数加以修正,即直接采用干土区的初始基质吸力水头值,且该值由张力计实测或由主要增湿曲线求得;需采用入渗试验前期的平均入渗率;所对应的湿润锋发展深度需经由水分传感器实测而得。研究结果表明该法能合理估算现场黄土的饱和渗透系数。     

降雨条件下考虑饱和渗透系数变异性的边坡可靠度分析

土体饱和渗透系数表现为天然的变异性,为此基于Green-Ampt模型建立了考虑饱和渗透系数变异性的降雨入渗物理模型,并藉此模型确定了坡体湿润锋深度和含水率分布。然后结合无限长非饱和土边坡稳定模型得到解析形式的反映边坡稳定性的极限状态函数。采用Monte Carlo法对饱和渗透系数进行随机抽样并最终建立降雨条件下考虑饱和渗透系数变异性的边坡概率分析框架。针对一假想边坡,探讨了饱和渗透系数的变异系数、降雨持时和降雨强度对边坡破坏概率以及破坏发生时间概率分布的影响,结果表明：在降雨初期,边坡的破坏概率随饱和渗透系数变异性的增强而逐渐增加,但随着降雨的持续,破坏概率开始随变异性的增强而显著降低;滑坡最可能发生时间的大小并不受饱和渗透系数变异性的影响,而是直接取决于降雨强度;滑坡最可能发生时间所对应的概率却随变异性的增强而逐渐减小。     

探讨变雨强条件下的入渗过程及影响因素

以积水条件下Green-Ampt入渗模型为基础,拓展其在变雨强中入渗率和累积入渗量计算形式。结合上述理论基础,选用饱和-非饱和土壤水、热和溶质运移模拟软件Hydrus-1D,探讨了常降雨与变雨强的入渗差异和饱和导水率、土壤类型对变雨强入渗过程的影响。研究表明：相对于常降雨,变雨强更利于雨水入渗,累积入渗量显著增加;其他水力参数不变,饱和导水率越大,入渗率与变雨强随时间变化曲线重合率越高,雨强控制时间越长,累积入渗量相应也越大;其他水力参数作用效果不容忽视,其影响作用可与饱和导水率达到相似效果。正确了解变雨强下入渗过程及影响因素,有助于准确地把握实际降雨入渗情况,对指导水力工程有重要的实际意义。     

降雨条件下滑坡排水沟排水数值模拟

地表排水沟是滑坡治理常见措施之一。在降雨入渗与坡面径流有限元耦合模型中,通过调整径流方程中的坡度参数,实现地表排水沟排水模型。模型具有易于理解、实现简单的优点。以地表排水模型为基础,对5种典型滑坡在设与不设排水沟情况下坡体渗流场、坡面径流场、坡体入渗量和入渗过程进行了数值模拟。结果表明,当降雨强度相对土体渗透性较小时,排水沟设置在坡体渗透性由小变大的位置,才能有效地截断上游径流,减小入渗量;而当降雨强度相对土体渗透性较大时,排水沟排水效果并不明显。     

适用于斜坡降雨入渗分析的修正Green-Ampt模型

降雨入渗分析是预测降雨诱发滑坡的关键因素之一。Green-Ampt模型原理简单、使用方便,在浅层滑坡的降雨入渗分析中有很大的应用潜力,但该方法主要适用于初始含水率为均匀分布的情况。基于这一不足,推导了初始含水率为非均匀分布条件下降雨入渗深度和时间的关系,并给出了基于Runge-Kutta原理的数值解法。当初始含水率为均匀分布时,提出的方法可简化为文献中已有的Green-Ampt模型。当初始含水率为非均匀分布时, Richards方程预测所得的孔隙水压力分布图中土体饱和区和未受降雨影响的非饱和区之间存在一个较窄的过渡段,由新模型计算所得的湿润锋穿过这一过渡段,且靠近饱和部区的底部。总体而言,新方法计算所得的孔隙水压力分布与Richards方程求解结果类似。