基于SWCC试验数据的坝体非饱和非稳态渗流与稳定性研究

饱和-非饱和状态是土存在于自然界的真实状态,描述和解释这种状态的非饱和土力学理论在土石坝渗流、污染物传输、冻土渗流相变、边坡和路基稳定性分析等领域有着广泛的应用。非饱和土的土水特征曲线(SWCC)是非饱和土力学研究的基本内容之一,对非饱和土体渗流和稳定性分析至关重要。采用Python语言开发了非饱和渗流与稳定性分析软件包USSA和非饱和土土水特征曲线试验数据处理和模型拟合界面。对非饱和土的土水特征曲线数据进行处理,以模型拟合参数作为基本输入进行非饱和渗流场的模拟,再到非饱和土的稳定性分析,详细呈现了非饱和土渗流与稳定性分析及软件开发的全过程。基于SWCC模型的拟合参数对坝体渗流稳定性进行了分析,分析结果表明水位变化过程中坝体两侧斜坡具有明显不同的稳定性演化规律。当采用有限元强度折减法进行非饱和土斜坡稳定性分析时,最终安全系数为所有边坡安全系数演化曲线的最低包络线。     

渠坡非饱和膨胀土含水率与强度关系试验研究

南水北调中线渠首段渠坡土主要为膨胀土,渠坡膨胀土含水率的变化会影响坡体稳定性。为了探究含水率变化对渠坡稳定性影响的具体特征,对采集自南水北调中线渠首段的渠坡膨胀土进行了滴定直剪试验,获得了渠坡膨胀土抗剪强度指标与含水率的关系曲线,试验结果表明：抗剪强度随含水率增加而明显衰减,衰减过程先快后慢;试样初始含水率越低,黏聚力下降越慢,内摩擦角下降越快,不同初始含水率试样的抗剪强度均在增湿至30%左右时产生拐点,此时衰减速率降低并趋于稳定。设计进行了土-水特征曲线法和瞬时剖面法导水率测试试验,获得了非饱和渠坡膨胀土导水率随含水率变化的曲线,试验结果表明：含水率越低,导水率越小,导水率变化速率越快;反之,导水率变化速率越慢,土体含水率趋于稳定。研究成果可用于膨胀土渠坡稳定性与坡体地下水位关系的定量分析,应用在实际工程中可以更有效地获取渠坡非饱和土体中含水率和抗剪强度的分布特征,并获得含水率和抗剪强度随时间的变化规律,为引入土体空间动态抗剪强度分析,建立更准确的膨胀土渠坡稳定性评价模型奠定基础。     

一种非饱和土表层蒸发过程的预测方法

探索了一种由土壤基本物性指标、气象数据预测非饱和土表层蒸发过程的方法。通过土壤基本物理性质（颗粒分布、土粒相对密度、干密度）预测出土-水特征曲线（SWCC）,进而得到土壤气相对湿度与含水率的关系,采用Penman-Wilson模型预测出了非饱和土表面的蒸发曲线。通过该方法,只需实地采取土样,获取其基本物理性质,由任意时刻土壤的含水率及气象数据就能预测出该时刻土表的蒸发速率。采用自制的自动蒸发测量系统进行浅层土蒸发试验,得到了蒸发曲线,与预测结果进行比较发现预测蒸发过程与实测结果一样皆由临界含水率和风干含水率将其分为3个阶段：稳定阶段、减速阶段和残余阶段,且蒸发量吻合,证明了所提出方法的准确性和实用性,对工程界估计土表蒸发量,确定水流量边界有重要意义。     

Design Charts for Selecting Minimum Setback Distance from Side Slope to Horizontal Trench System in Bioreactor Landfills

The primary objective of bioreactor landfill is to achieve adequate and rapid distribution of moisture in landfilled municipal solid waste (MSW) to accelerate the anaerobic biodegradation of the organic fraction within MSW. A horizontal trench system (HT) is commonly adopted for leachate distribution in MSW under pressurized conditions. However, this approach should be implemented carefully due to the potential instability of landfill slopes that comes from the generation and distribution of excessive pore fluid pressures. In this study, HT design charts are presented that determine the optimal location of horizontal trench systems from the side slope (i.e., minimum lateral setback distance) under continuous leachate addition with maximum applied injection pressures, for which the landfill slopes remain stable [factor of safety (FOS) where FOS ≥ 1.5]. Use of any higher injection pressure and/or shorter lateral setback distance of HT than the one presented in the design charts would result in an unacceptable design of the bioreactor side slope (FOS < 1.5). The design chart was developed based on a parametric study that used a numerical two-phase flow model that involved different slope configurations and landfill waste depths. MSW heterogeneity and anisotropy, as well as unsaturated hydraulic properties, were taken into consideration in these simulations. Transient changes in pore water and gas pressures due to leachate recirculation were accounted for dually in the slope stability computations. The importance of these design charts is illustrated using a practical example. Site-specific conditions and the expertise and prior experience of a designer or operator must also be adequately considered and utilized with the design charts presented here for the safe design of a horizontal trench system in a bioreactor landfill.     

非饱和成都黏土瞬时剖面法渗透性试验研究

非饱和渗透系数是土体渗流分析的基础,成都黏土作为一种典型的非饱和膨胀土,具有吸水膨胀、失水收缩的特性,在受侧限的浸水过程中,土颗粒的膨胀致使孔隙体积减小,渗透性降低,使得直接对其进行非饱和渗透试验十分困难。根据瞬时剖面法的原理,利用EC-5土壤水分传感器测含水率、MPS-2电介质水势传感器直接同步测量同一位置的基质吸力,通过水平渗透试验研究了非饱和成都黏土在侧限条件下的渗透性。含水率和基质吸力的同步测量,保证了其测试条件的一致性,避免了采用其他土水特征曲线的影响。试验表明,试样的非饱和渗透系数为（1.33×10-11~3.14×10-9）m·s-1,非饱和渗透系数与基质吸力并非单调线性关系。基质吸力较高时,受膨胀土颗粒吸水膨胀的影响,渗透系数未出现明显变化,基质吸力降低到一定程度后,渗透系数快速增大。试验结束时土体已接近饱和,土中气体排出较慢,过水断面增加缓慢,促使渗透系数仍然持续增大。采用VG模型拟合k-s曲线,拟合参数α=0.048 kPa-1,n=1.79,m=0.48,试验结果可以用于成都黏土地区的渗流分析。     

Influence of Unsaturated Soil Properties Uncertainty on Moisture Flow Modeling

The quality of a numerical modeling solution of moisture flow through unsaturated soil, in part, depends on properly described unsaturated soil properties. The variability of the Soil Water Characteristic Curve, SWCC, is attributed to hysteresis and reproducibility of measurement. Because the unsaturated conductivity function is rarely directly measured, the variability of the unsaturated soil hydraulic conductivity function is attributed to the uncertainty associated with the estimation of this parameter with currently available fitting functions, and hence a range of reasonable variation was considered. One-dimensional modeling of expansive soil under dry initial conditions (suction of 1,500 kPa) was performed; both potential evaporation and infiltration boundary conditions were considered. It was found that small variations in the unsaturated soil hydraulic conductivity function result in significantly different modeling outputs, as expected, while substential variation in SWCC alone (assuming the same unsaturated soil hydraulic conductivity for all SWCCs) produced almost identical soil response in terms of soil suction when the slope of the SWCC is similiar. Thus, proper characterization of the slope of the SWCC is important to proper suction profile determination.     

Numerical Simulation of Moisture Movement in Unsaturated Expansive Soil Slope Suffering Permeation

INTRODUCTION Theinfiltrationandevaporationofwaterinasoil slopeareofparamountimportanceindetermining slopestability.Previousengineeringcasesandstud ieshaveshownthatrainwaterisoftenamajorfactor intheslopefailureofexpansivesoils.Expansivesoils intropicalandsubtropicalzonesareoftenunsaturat ed,andsoilslopesareinastablestateinnormalcli mateconditionsbecauseofthehighsuctionandshear strengthofthesoilmass.However,oncepermeation happensrainwaterwillinfiltrateintothesoilmasswhichleadstoanincreasein…     

Study on the selection of unsaturated flow model for the different types of soil and soft rock

Precise estimation of unsaturated hydraulic properties of porous media is indispensable in various study areas, such as analyzing the moisture flow, the drying process occurring from the surface, and the pollutant migration beneath the ground surface. Although many empirical/theoretical models describing the unsaturated hydraulic properties have been proposed by several previous researchers, the best model for the different types of soil/rock may not be identical. Thus, the model selection process and the estimation technique of the parameters included in the models should be developed. In the present study, the inverse technique based on the transient evaporation change was investigated to select the model and estimate the model parameters. The experimental work was based on a relatively low permeable soft rock and a relatively high permeable sandy soil (Toyoura standard sand). Experimental equipment was developed to precisely measure the evaporation rate for the high permeable sandy soil. The Genetic Algorithm (GA) was adopted in the inverse technique as an optimization tool. In order to simplify the problem, only the drying process from the saturated condition was considered. It was established that the information concerning the transient evaporation change could be used for the model selection and parameter estimation. Further, the saturation distribution could be used for the selection of the models. The present study provides important information for the development of the model selection process.     

基于孔隙分形维数的土壤大孔隙流水力特征参数研究

根据数字图像制备分析技术得到大孔隙数目、大小和形状等特征数据,用分形理论计算孔隙分形维数,并对两类不同质地土壤的孔隙分形维数进行分析比较,应用基于孔隙分形维数的大孔隙分形模型,建立了土壤水分特征曲线和非饱和导水率的预测模型。结果表明,孔隙分形维数D_v可定量描述不同质地土壤结构,其模型预测精度优于颗粒大小分布估计的分形模型,可用于实际问题的研究。     

Design Charts for the Stability Analysis of Unsaturated Soil Slopes

Simple limit equilibrium analyses can be performed to determine the Factor of Safety (FOS) against slope failure of unsaturated soil slopes. However, many of the input parameters needed for these analyses are highly variable, and the FOS value obtained is critically dependent on assumptions made by the designer. This paper describes a suite of reliability analyses on unsaturated soil slopes performed using an invariant reliability model. The results are presented in design charts from which a designer can choose the FOS value required to ensure a given target reliability index for a slope. The approach ensures that despite the variability of input parameters the slope will have a probability of failure of 2.23% or less.     

基于粒径分布曲线的非饱和砂土土水特征曲线概率预测模型

土水特征曲线定义了非饱和土的基质吸力和含水量之间的关系,与非饱和土的渗流和强度等特征有密切关系.本文通过对100组砂土的粒径分布曲线和土水特征曲线进行分析,结合常用的VG模型提出了基于粒径分布曲线的非饱和砂土土水特征曲线概率预测方法,并基于另外30组数据对提出的模型进行了验证.研究表明,基于粒径分布曲线无法唯一确定土体...     

东南沿海残积土地区降雨诱发型滑坡预报雨强-历时曲线的影响因素分析

东南沿海山地丘陵地区每年雨季期间有大量土质滑坡发生,如何对由降雨诱发的滑坡进行预报一直都是一个热点研究课题,雨强-历时曲线（简称I-D曲线）是目前国内外常用的降雨型滑坡预报的降雨量临界值曲线。针对东南沿海地区的浅层残积土滑坡,根据相关勘察数据及资料,概化得到了该类型边坡的地质剖面及岩土层性质,然后应用Geo-Studio软件分析了边坡初始湿润条件、土体抗剪强度、饱和渗透系数、边坡坡角、残坡积土层厚度及雨型等参数对I-D曲线的影响规律。分析结果表明：残坡积土抗剪强度参数、饱和渗透系数、边坡坡角、雨型等因素对I-D曲线的影响显著,边坡安全系数降至临界值所需降雨历时随抗剪强度参数值降低、表层残积土渗透系数增加或坡角增大而减少;当雨强较小时,初始湿润条件对I-D曲线的影响显著;当雨强大于残坡积土层饱和渗透系数时,入渗量主要由渗透系数控制,边坡安全系数降至临界值所需降雨历时不随雨强增大而变化。该研究结果为I-D曲线在东南沿海残积土地区降雨诱发滑坡预警预报中的应用奠定了基础。     

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

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

降雨入渗条件下非饱和土边坡稳定分析

针对降雨入渗土坡的稳定问题,建立一个考虑水力渗透系数特征曲线、土-水特征曲线以及修正的Mohr-Coulomb破坏准则的非饱和土流固耦合有限元计算模型,进行雨水入渗下非饱和土边坡渗流场和应力场耦合的数值模拟,得到非饱和土边坡变形与应力的若干重要规律。研究成果为降雨入渗条件下非饱和土边坡的稳定分析提供了基础。     

均质非饱和土一维稳态流水头垂直分布的解析解

渗流场水头分布计算是进行渗流量和渗流水力坡降计算的基础,准确、有效地求取渗流场水头分布是渗流计算的关键环节。对均质非饱和土体一维稳态流的流动方程进行分析,考虑到渗透系数是与基质吸力相关的函数,通过数学变换,给出了稳定渗流场的解析通式,并基于渗透性函数中的Gardner模型,给出了非饱和土一维稳态流水头垂直分布的解析解。该解析通式表明,均质非饱和土一维稳态流水头垂直分布主要受地表水头、深度和流动率3个因素控制。分别计算了一维稳态蒸发条件下粉土和黏土两种典型土类水头沿垂直方向的分布。计算结果表明：稳态蒸发条件下粉土层和黏土层内的水头分布表现出相似的变化规律,即自地表至地下水位处随着土层深度的增加,水头分布呈现出加速递减的趋势;在相同的蒸发条件下,对于相同深度处的黏土和粉土而言,黏土层内水头更高些;对同一种土类而言,在较大的蒸发状态下同一深度处土层内水头更高。反之,则较低。     

基于分层假设的Green-Ampt模型改进

传统Green-Ampt入渗模型由于假设过于简化而限制了计算精度,针对这一问题,在分层假设的基础上对模型进行了改进:将均质入渗土壤分为饱和层、过渡层和干土层,采用Richards模型分析了不同入渗条件下各层的厚度以及各层内含水率和导水率的变化规律,结果表明,入渗条件一定时,过渡层占湿润层(包括饱和层和过渡层)比例会随湿润层厚度的增加线性减小,积水深度和土壤初始含水率分别会影响这一线性关系的纵轴截距和斜率,而土壤饱和导水率则对其没有明显影响;过渡层内土壤含水率分布为椭圆形曲线,而导水率随深度增加线性减小。在此基础上,通过引入湿润层的等效导水率并相应给出锋面吸力的计算方法,对传统Green-Ampt模型进行了改进。结果显示,改进后模型计算精度显著提高。     

基于COMSOL Multiphysics的非饱和裂隙土降雨入渗特性研究

基于COMSOL Multiphysics软件对非饱和裂隙土降雨入渗特性进行数值模拟研究。通过将裂隙和基质分别离散成有限单元,建立了能充分模拟土中裂隙流、基质流以及裂隙-基质流量交换的离散裂隙-孔隙介质模型。结合"空气单元"的概念,对裂隙土的上边界进行模拟。该方法不仅能描述降雨初期雨水沿裂隙优先入渗的现象,还能描述当降雨量大于裂隙土入渗量时雨水沿地表流走的现象。通过对地表以下2 m深度内低渗含裂隙土体进行模拟,分析了裂隙的几何特征、基质的水力特性、前期水分条件以及降雨强度对非饱和裂隙土降雨入渗过程的影响。结果表明,在非饱和裂隙土中,存在两个主要的渗流过程:一是水沿裂隙优先流动;二是水不断从裂隙吸入基质中,基质吸收水的作用抑制了裂隙中优势流的发展。与裂隙的几何特征相比,基质的水力特性对非饱和裂隙土渗流的影响较大。增大基质的饱和渗透系数可能使由裂隙流主导的渗流过程转变为由基质流主导的渗流过程,而基质的非饱和特性与裂隙土的初始含水率改变了土体的储水能力,从而加速或延缓了降雨入渗至某一深度的时间。降雨强度对土体入渗速率和入渗量均有影响,当超过裂隙土的入渗能力时,多余积水沿地表流走,断面入渗率随...     

水库蓄水过程中堆积体边坡瞬态稳定性分析

由现场试验结合Fredlund和Xing方法确定了堆积体的非饱和渗透性函数分布曲线。根据概化的地质模型,建立了非饱和渗流分析的有限元模型,模拟得到蓄水过程中不同时刻坡体中孔隙水压力的分布。把暂态孔隙水压力的分布和非饱和土强度理论应用到普遍极限平衡法（GLE）中,进行了库水上升过程中边坡瞬态稳定性分析。结果表明,库水上升过程中,坡体整体稳定性是下降的,而在受库水变动影响的局部区域,其稳定性系数随着库水位的上升表现出先减小后增大的趋势,即存在一个危险水位,在一定条件下,可能会诱发坡体下部局部区域失稳。     

考虑参数空间变异性的非饱和土坡可靠度分析

在考虑多个土体参数空间变异性的基础上,提出了基于拉丁超立方抽样的非饱和土坡稳定可靠度分析的非侵入式随机有限元法。利用Hermite随机多项式展开拟合边坡安全系数与输入参数间的隐式函数关系,采用拉丁超立方抽样技术产生输入参数样本点,通过Karhunen-Loève展开方法离散土体渗透系数、有效黏聚力和内摩擦角随机场,并编写了计算程序NISFEM-KL-LHS。研究了该方法在稳定渗流条件下非饱和土坡可靠度分析中的应用。结果表明：非侵入式随机有限元法为考虑多个土体参数空间变异性的非饱和土坡可靠度问题提供了一种有效的分析工具。土体渗透系数空间变异性和坡面降雨强度对边坡地下水位和最危险滑动面位置均有明显的影响。当降雨强度与饱和渗透系数的比值大于0.01时,边坡失效概率急剧增加。当土体参数变异性或者参数间负相关性较大时,忽略土体参数空间变异性会明显高估边坡失效概率。     

邯郸地区土壤水入渗规律实验及数值模拟方法研究

在对土壤入渗理论进行总结和分析的基础上,依据邯郸市平原区的水文水情的基本特点,选定2~3处地点进行实验,完成了静态数据的收集和动态数据的观测,确定了实验区土壤质地类型、土壤容重等物理特性,结合实验测定的饱和导水率K、土壤水分特征曲线θ（ h）推算非饱和导水率K（θ）、扩散率 D（θ）等土壤水分运动参数。运用上述数据对Van Genuchten模型加以改进,并运用改进后的模型对实验中的水分变化过程进行数值模拟。结合实测成果和数值模拟成果,对邯郸市平原区已选定实验点土壤水入渗规律进行分析。