非饱和土渗流研究进展

高边坡山体饱和非饱和渗流场的初步分析_张家发

◇　科　研 高边坡山体饱和非饱和渗流场的初步分析 张家发　李思慎　叶自桐 摘　要　对一高边坡山体花岗岩全风化带土样和基岩裂隙测得了水份特征曲线,并进而得到了水力传导 率的解析模型。采用有限元方法模拟分析了类似多年平均降雨条件下高边坡山体中的稳定渗流状态,以及强降雨入渗条件下的非稳定入渗过程和渗流场变化趋势,并对设计的排水措施的效果进行了初步分析。模拟结果说明,在一定条件下现已布置的排水措施作用可能有限。针对这一情况,为下一步的观测和研究工作提出了相应的建议。 主题词　渗流观测　排水设施　有限元法　风化岩　试验研究 边坡体中的水压力和水流分布是影响边坡稳定和变形的因素之一。因此边坡排水设计及其方案论证是边坡设计与科研的重要内容。以往高边坡渗流场研究通常采用的是稳定流模型。对降雨入渗补给的作用,仅考虑多年平均降水量对应的入渗条件,且入渗边界假设在地下水面上。实际上在一些山地暴雨区,雨量丰沛且在时间分布上很集中,以阵发性暴雨为主。强降雨过程中高边坡山体接受入渗补给以及边坡体中的地下水瞬态运动和水压分布将会更加恶化边坡的排水条件和边坡体的稳定状态。另一方面,在多山地区,尤其是在工程开挖边坡造成的地形深切割条件下,旱季地下水位通常是很低的,形成了深厚的非饱和区。在继之而来的雨季强降雨过程中,非饱和区的水份运动将对降水入渗补给过程以及地下水压力分布发生影响。考虑这些影响后,在以往研究成果基础上设计的排水措施的效果如何?这是本文的研究重点。 本文在已往研究成果基础上,开展降雨入渗条件下高边坡岩体饱和非饱和渗流的研究。由于这是参数高度非线性问题,且非饱和参数的测定不是常规地质勘探工作的内容,通常资料很少,此次工作只能是初步的。文中首先进行了非饱和参数研究;然后在参数资料不足的情况下通过假设,用有限元数值模拟分析了类似多年平均降雨强度下边坡中的稳定渗流状态以及强降雨过程中的非稳定入渗补给过程及边坡中地下水运动和水压分布的变化趋势。 1　参数的研究 本文所研究的高边坡岩体为花岗岩体,其渗透性与岩体风化程度和裂隙发育状况关系密切,总体上可分为与全强风化带、弱风化带和微新岩体相对应的三个渗透性分区。每一渗透性分区的非饱和参数包括K (θ)和h (θ)这两组非线性函数,通常它们是分别通过实验测定的。本文首次介绍了对该高边坡岩体的饱和非饱和渗流参数试验研究,包括分别对新鲜岩体的裂隙和全强风化带的松散介质进行的室内实验。 新鲜岩样中有一个完整的贯穿裂隙面。通过饱和渗流试验 并根据立方定理推算出裂隙的水力开度为179μm 。假设通过岩样的基质孔隙的水流可以忽略不计,采用不互溶驱替法,用 非湿润流体(变压器油)驱替裂隙中的湿润流体(水),通过流量和水压力的观测得到裂隙的毛管压力～饱和度关系(如图1所示)。然后采用va n Ge nuch ten 模型[1]和Broo ks -Co rey 模型[2]对实验数据进行拟合,进而得到了裂隙水力传导率的解析模型(如图2所示) 。 图1　 裂隙负压水头与饱和度关系实验曲线 图2　裂隙的K r 曲线 将全强风化带土样按 1.88的容重装填,采用压力板出流 · 44·V o l .29　N o .1 人民长江　Y AN G T Z E RIV ER J a nuar y 1998DOI:10.16232/https://www.wendangku.net/doc/db16710798.html, k i .1001-4179.1998.01.016

非饱和渗流

非饱和渗流中渗透系数计算的推导 （1） 拟合 由实验测出测点的含水率和基质吸力的实验数据，所测得的含水率可算出其有效饱和度，即有效饱和度可由含水率表示出来，然后再用VG 模型拟合出土体的水分特征曲线SWCC 。 式中: 为有效饱和度，，为基质吸力。拟合出VG 模型中的三个参数，即可得到有效饱和度与基质吸力的关系SWCC 曲线。 用所得到的有效饱和度，再由VG 模型可得到相对渗透系数与有效饱和度的关系 而非饱和渗透系数与相对渗透系数的关系是： k w = k r w k s 由土常规物理实验可测得土体的饱和渗透系数，即可得到非饱和渗透系数与含水率的函数。 （2）达西定律直接计算 由法国水力学家 H.-P.-G.达西在1852～1855年通过大量实验得出，1856年总结得出渗透能量损失与渗流速度之间的相互关系即为达西定律。反映水在岩土孔隙中渗流规律的实验定律。这个定律说明水通过多孔介质的速度同水力梯度的大小及介质的渗透性能成正比。 达西定律是多孔介质中流体所应满足的运动方程。质量守恒是物质运动和变化普遍遵循的原理，将质量守恒原理具体应用在多孔介质中的流体流动即为连续方程。达西定律和连续方程相结合便导出了土体中水分运动的基本方程。 根据达西(Darcy)定律和质量守恒定律，对于二维问题非饱和土壤水运动的基本微分方程如下： ()()x y K K t x x y y θ??θθ?????????=+????????????? 式中θ为体积含水量；φ为总水势(总水头)，由基质势和重力势组成，φ= y+h ，y 为重力势(位置势)，h 为基质势；x K ，y K 为x ，y 方向的渗透系数，若土体为各向同性，则x K =y K =K （θ） 由于非饱和渗透系数是基质吸力或者含水率的函数，故此方程为一个二阶非线性的偏微分方程，除少量问题外，一般情况下对此方程的解析求解是困难的，很多的问题需要用数值法求解。 由于非饱和土的渗透系数K 可以是基质吸力(负压水头)的函数，因此方程(5.1)的左端可以改写为：

comsol在非饱和土渗流的应用

基于comsol的非饱和土渗流研究 /comsol在岩土工程渗流的应用 摘要：岩土工程的核心难点即解决地下水问题，一般岩土工程事故都是由于对地下水的影响重视不够而造成的，然而解决这一难点关键在于解决地下水渗流问题。目前对于非饱和土渗流研究的理论仍相对落后，本文结合非饱和土渗流场基本方程以及由水土特征曲线得到的相关渗流参数（渗透系数，体积含水量），阐明了如何解决渗透模型要求渗流场方程的连续性与现场实测数据的非连续性之间的矛盾，并利用comsol Multiphysics 软件对某工程中非饱和土渗流问题进行了模拟，并验证了Fredlund和xing（1994）土水特征曲线方程的正确性。这种解决非饱和土渗流问题的思想可供学者参考。 关键词：非饱和土；渗流场；渗流参数；连续性矛盾；Comsol Multiphysics Study on seepage of unsaturated soil seepage based on comsol Abstract:The core difficulty of geotechnical engineering is to solve groundwater problems, the general geotechnical engineering accidents are due to the impact of groundwater caused by insufficient attention, however, the key to solve this difficult problem is to deal with the groundwater flow. At present, for the study of unsaturated soil seepage theory is still relatively backward, this paper combines basic equation of unsaturated soil seepage with soil-water characteristics curve and obtains the relevant flow parameters (hydraulic conductivity, volumetric water content) from them, and illustrates how to solve the conflict between the seepage field penetration model requiring Equation of continuity and the measured data of non-continuity, and using the software comsol Multiphysics to simulate unsaturated soil seepage problems in one project and verified the right of Fredlund and xing (1994) soil-water characteristic curve equation. The idea of solving unsaturated soil seepage problems may be referred by similar projects. Key words: unsaturated soil; seepage field; seepage parameters; continuous conflict; Comsol Multiphysics 1引言 岩土工程设计与施工的难点在于解决地下水问题，一般岩土工程事故都是由于对地下水的影响重视不够而造成的，像2003年7月14日上海轨道交通４号线工程事故；2007年8月17日山东新汶煤矿透水事故；2008年11月15日杭州地铁工地塌陷事故以及2011年1月1日杭州余杭区-工地土方坍塌事故等等都是由于忽视地下水的影响而造成的。然而解决这一难题的关键在于解决地下水渗流问题。虽然众多国内外学者对土的渗流问题做了大量的研究，但是目前对于非饱和土渗流研究的理论以及实践应用仍相对落后。一般来说，解决非饱和土渗流设计的问题以及与其相关的工程实践问题，可以归结于就具体的非饱和土渗流工程概况而建立渗流场基本方程，然后解这一渗流场基本方程，从而得出相关的渗流流线（水位）分布、水流渗流力矢量分布、流速矢量分布和相关的趋势，最后以此来指导实践施工。然而在求解非饱和土渗流场基本方程时，首要要解决两个重要未知参数，即体积含水量θ和渗透系数k，这两个参数在实际工程中是通过实验得到的，试验得到的是一系列孤立的点，然而这与渗流场基本方程建立于连续性模型相悖，这就给求解渗流场基本方程带来了很大的困难，于是国内外很多学者对此进行了大量的研究。为了解决非饱和土的

土的渗透性和渗流问题

第四章 土的渗透性和渗流问题 第一节 概述 土是由固体相的颗粒、孔隙中的液体和气体三相组成的，而土中的孔隙具有连续的性质，当土作为水土建筑物的地基或直接把它用作水土建筑物的材料时，水就会在水头差作用下从水位较高的一侧透过土体的孔隙流向水位较低的一侧。 渗透：在水头差作用下，水透过土体孔隙的现象 渗透性：土允许水透过的性能称为土的渗透性。 水在土体中渗透，一方面会造成水量损失，影响工程效益；另一方面将引起土体内部应力状态的变化，从而改变水土建筑物或地基的稳定条件，甚者还会酿成破坏事故。 此外，土的渗透性的强弱，对土体的固结、强度以及工程施工都有非常重要的影响。 本章将主要讨论水在土体中的渗透性及渗透规律，以及渗透力渗透变形等问题。 第二节 土的渗透性 一、土的渗透规律——达西定律 (一)渗流中的总水头与水力坡降 液体流动的连续性原理：（方程式） dw v dw v w w ??=2 211 2211v w v w = 1 221w w v v = 表明：通过稳定总流任意过水断面的流量是相等的；或者说是稳定总流的过水断面的 平均流速与过水断面的面积成反比。 前提：流体是连续介质 流体是不可压缩的； 流体是稳定流，且流体不能通过流面流进或流出该元流。 理想重力的能量方程式（伯努利方程式1738年瑞士数学家应用动能定理推导出来的。） c g v r p Z =++22 饱和土体空隙中的渗透水流，也遵从伯努利方程，并用水头的概念来研究水体流动中 的位能和动能。 水头：实际上就是单位重量水体所具有的能量。 按照伯努利方程，液流中一点的总水头h ，可以用位置水头Z ，压力水头U/r w 和流速水

成层非饱和土渗流的耦合解析解

2011年8月 Rock and Soil Mechanics Aug. 2011 收稿日期：2010-06-07 基金项目：国家自然科学基金（No. 40902087）；香港Research Grants Council （No. 622207）；教育部科学技术研究重点项目（No. 110186）；教育部博士点新教师基金（No. 20095122120007）；四川省杰出青年学术技术带头人培育计划（No. 2010JQ0034）。 第一作者简介：吴礼舟，男，1975年生，博士，副教授，主要从事工程地质和岩土工程科研教学工作。E-mail: wulizhoucn@https://www.wendangku.net/doc/db16710798.html, 文章编号：1000－7598 (2011) 08－2391－06 成层非饱和土渗流的耦合解析解 吴礼舟1，张利民2，黄润秋1 （1.成都理工大学 地质灾害防治国家重点实验室，成都 610059；2.香港科技大学 土木系，香港 九龙） 摘 要：成层土在工程中很常见，研究降雨过程中成层非饱和土的渗流-变形耦合对非饱和土土力学的发展具有重要的意义。由流体质量守恒，Darcy 定律和Lloret 等的非饱和土本构模型可得成层非饱和土渗流-变形耦合的控制方程。采用Gardner 的非饱和土的渗透系数公式以及Boltzman 模型，基于Laplace 变换得到耦合方程的解析解。解析及其参数分析表明，渗流和变形耦合是具有时间效应的。与吸力变化相关的土的模量F ，对成层土的孔隙水压力分布有明显影响。两层土的F 差异越大，孔隙水压力消散得越慢，耦合效应越不显著。增大表层土的F 值有利于降低耦合效应。成层土饱和体积含水率变化对吸力变化产生有限的影响。 关 键 词：非饱和土；渗流和变形；耦合；成层土；降雨入渗 中图分类号：TU 46+2 文献标识码：A Analytic solution to coupled seepage in layered unsaturated soils WU Li-zhou 1 , ZHANG Li-min 2, HUANG Run-qiu 1 (1. State Key Laboratory of Geological Hazard Prevention and Geological Environment Protection, Chengdu University of Technology, Chengdu 610059, China; 2. Department of Civil and Environmental Engineering, HongKong University of Science and Technology, Kowloon, Hong Kong, China) Abstract: Layered soil such as landfill and cracked soils are often found in engineering. Its coupled infiltration and deformation during rainfall is significant for development of unsaturated soil mechanics. Based on fluid mass conservation, Darcy’s law, and the constitutive model proposed by Lloret et al., coupled governing equations for seepage and deformation in unsaturated soils are obtained. The unsaturated coefficient of permeability is expressed using Gardner’s model and the water retention characteristics are expressed using Boltzman’s model. The analytic solution to the coupled equation is developed by Laplace transformation. The analytic solution and parameter analysis results show that the effect of coupling between unsaturated seepage and deformation is related with time. The modulus related to suction changes F has a marked effect on the pore water pressure. The larger the F ratio values for two-layer soils are, the more slowly the suction dissipates. The results indicate that a large F for the top-layer soil can effectively reduce the coupling effect. V olumetric moisture content changes in two-layer soils play a limited role in the suction distribution. Key words: unsaturated soil; seepage and deformation; coupling; layered soil; rainfall infiltration 1 引 言 非饱和土在地球表面广泛分布，非饱和土由土骨架、孔隙水、孔隙气和水气膜组成[1]。降雨入渗过程中非饱和土体中的渗流场和位移场均是变化的，且相互影响。孔隙水压力变化导致应力变化及非饱和土体变形；应力变化及孔隙改变反过来又影响渗流状态[2]。因此，降雨过程中非饱和土的渗流-变形耦合问题是一个重要的课题，成层土常见于工程中，如垃圾、废料填埋土。干湿循环易诱发表层 土产生裂隙，因而表层裂隙的土层作为连续介质，其渗透系数增大，与下层未出现裂隙土层一起构成等效的成层土。研究成层非饱和土的渗流-变形耦合对推动非饱和土土力学的发展有着重要的意义。 关于渗流-变形耦合的数值解有不少研究[3 －6] 。 Kim [3]提出非饱和土地表加载引起水位波动和变形的耦合数值模型。Thomas [4]提出了可变形的非饱和土热、水和气转化的理论表达式。田东方等[5]提出坡面径流-非饱和渗流分析与应力场的耦合计算方法，并编制了相应的有限元程序。张玉军[6]从建立

ABAQUS在饱和-非饱和渗流分析中的应用

ABAQUS 在饱和－非饱和渗流分析中的应用 徐海奔 河海大学水工结构工程专业，南京 (210024) E-mail ：hohaixhb@https://www.wendangku.net/doc/db16710798.html, 摘 要：本文首先对大型通用有限元软件ABAQUS 在土石坝渗流分析中的应用进行分析，着重从多孔介质的饱和渗流，非饱和渗流及二者的混合问题（渗流自由面的计算）等方面论述。结合一个土石坝库水位下降时二维渗流计算实例，考虑流体重力作用下，采用非线性定律求解总孔隙压力及库水位下降过程渗流自由面变化过程。 关键词：非饱和；渗流；ABAQUS ；土石坝；自由面 1．引言 ABAQUS 大型通用有限元软件,在我国土木工程结构分析方面应用日益广泛。本文对它在土石坝渗流计算分析中的应用进行评述。 近年来，在国内外随着孔隙介质非饱和渗流和土体饱和渗流理论的发展，人们逐渐意识到堤坝稳定性与非饱和区渗流作用密切相关。在研究堤坝非饱和渗流问题时，主要采用数值模拟的方法。长期蓄水的土坝，当库水位以太快的速度下降时，坝体内孔隙水压力常常不能很快消散，因而坝体的浸润线高于上游库水水位。在这种情况下，渗流的动水压力或渗透力的作用对上游坝坡造成浮起及下滑的趋势，甚至酿成滑坡事故。因此在实际工程中必须防止因库水位下降速度太快而导致这类事故发生。为进行上游坝坡的稳定分析，需要确定库水位下降过程中各时段坝体浸润线的位置，也就是通常所说的进行土坝不稳定渗流计算。 坝体浸润线下降的速度，一般决定于库水位下降的速度V 、土坝坝体渗透系数k 以及土体的给水度u 等因素[1]，与坝体的结构形式特别是坝体及地基上游面的排水条件也有很大关系。 2．ABAQUS 在均质土坝饱和－非饱和渗流计算原理 在饱和土壤中，引起水分转移的力是重力和水的压力。在非饱和土中，支配着土壤水在液态下整体转移的是重力和水的表面张力。Richards 等曾在1931年就证明非饱和土中的渗流与饱和土一样符合达西定律和连续方程[2]。若将达西定律代入连续方程（忽略渗透过程中总应力的改变和土颗粒骨架的变形）并以总水头h 作为未知量，当渗透的主方向与坐标轴一致时，非饱和土渗流的二维微分方程就可表示为： t y h k x x h k x w y x ??=????????????+??????????θ （1） 式中，x k ，y k 分别为x ，y 方向的渗透系数；w θ为体积含水量；h 为总水头；t 为时间。 令y 为位置水头，则：y u h w w +=γ，若w m 为土水特征曲线斜率，则： ()y h m u m w w w w w ??=?=?γθ。式（1）就可以写为： ()t y h m y h k x x h k x w w y x ???=??? ?????????+??????????γ （2） 因为y 为常数，式（2）可简化为：

裂隙岩体非饱和渗流研究综述_胡云进

收稿日期:19990614 基金项目:国家自然科学基金资助项目(59879004);高等学校博士点学科专项基金资助项目(98029408);水利部水利技术开发基金资助项目(97472603).作者简介:胡云进(1974—),男,浙江东阳人,博士研究生,水力学及河流动力学专业,主要从事裂隙岩体非饱和渗流研究. 裂隙岩体非饱和渗流研究综述 胡云进,速宝玉,詹美礼 (河海大学水利水电工程学院,江苏南京210098) 摘要:综述了国内外裂隙岩体(主要是细、微裂隙岩体)非饱和渗流的研究情况.首先,评述了现有的 测定和确定单裂隙非饱和水力参数的几种方法的优缺点,为单裂隙非饱和水力参数的确定提供了 理论依据.其次,分析了目前用于求解裂隙岩体非饱和渗流的四种数学模型的优缺点,为选取合理 的数学模型用于求解具体的裂隙岩体非饱和渗流问题提供了参考依据.最后,在上述基础上提出了 一些需要进一步研究的问题. 关键词:裂隙岩体;非饱和渗流;非饱和水力参数;数学模型 中图分类号:TV139.1 文献标识码:A 文章编号:10001980(2000)01004007 由于构造、风化、卸荷等作用,天然岩体中存在着大量的裂隙.裂隙岩体中地下水位以上部分是未被水充满的非饱和带,降雨入渗和地面水体的下渗都是通过该带到达稳定地下水面的非饱和渗流过程.可以说,在自然界中裂隙岩体非饱和渗流是普遍、客观存在的.以往,由于裂隙岩体非饱和渗流的复杂性,在许多工程问题中都作了简化处理[1].近十多年来,随着西方工业发达国家核电等核工业的不断发展,大量核废料亟待深埋处理(一般埋在裂隙岩体深厚的非饱和带中),故评价核废料深埋对地下水环境的污染以及处理核废料的选址等都要求对裂隙岩体非饱和渗流作深入细致的研究[2,3,4].另外,研究降雨入渗对地面污染物的淋滤[5],石油二次开采[6]和地热能开发[7]等也涉及裂隙岩体非饱和渗流或两(多)相流问题.在国内,近年来随着孔隙介质非饱和渗流和裂隙岩体饱和渗流理论的发展,已越来越清楚地认识到雨季的岩坡滑坡、地下洞室巷道的塌方以及泄洪雾化雨导致岩质边坡的失稳等均与裂隙岩体非饱和渗流密切相关[1](即降雨入渗等会导致地下水位以上非饱和区孔隙水压力的升高,产生暂态的附加水荷载,同时降低岩体的力学强度指标). 由于上述工程应用领域的需要,国外已有不少学者相继从80年代中期开始对裂隙岩体非饱和渗流进行试验和理论研究.近年来,国内也有学者开始了这方面的研究工作.目前的工作主要有:(a )对单裂隙非饱和渗流进行试验和理论研究,主要集中在单裂隙非饱和水力参数(即毛细压力饱和度和相对渗透率(非饱和渗透率与饱和渗透率的比值)饱和度(或毛细压力)的函数关系)的测定和确定方面;(b )提出各种求解裂隙岩体非饱和渗流的数学模型并进行相应的数值分析. 1　单一裂隙非饱和渗流研究 在裂隙岩体非饱和渗流研究中,最关键的是单裂隙毛细压力饱和度和相对渗透率饱和度(或毛细压力)关系的建立.目前,建立上述关系主要有以下三种方法:(a )物模试验法,即直接通过单裂隙拟稳态驱替试验和非饱和渗流试验(准确地说是二相流试验),借用孔隙介质拟合模型拟合出经验关系式;(b )数值试验法,即通过建立单裂隙概化模型,利用数值模拟法和孔隙介质拟合模型拟合出经验关系式;(c )数学推导法,即在某些假设简化的前提下,根据裂隙开度分布推导出上述关系式. 1.1　物模试验法 天然裂隙壁面是凹凸不平的,两粗糙裂隙面间的空隙空间的开度是逐点变化的[8,9].天然裂隙可概化为第28卷第1期 2000年1月河海大学学报JOUR NAL OF HOHAI UNIVERSITY Vol .28No .1Jan .2000

非饱和土力学

石家庄铁道大学 研究生课程考试答题纸 培养单位土木工程学院 学科专业岩土工程 课程名称非饱和土力学 任课教师冯怀平 考试日期 2011-12-28 学生姓名李幸吉 学号 120110404 研究生学院 4 非饱和渗流场与应力场耦合分析

非饱和渗流场与应力场的耦合机制为：介质内渗流会产生渗透力，进而改变应力场及其分布；应力场的变化又引起介质孔隙率等的变化，从而改变其渗透系数及渗流场，这一点可以从考虑了应力的土-水特征曲线（即式（3））得到体现。从而两场是耦合的。目前，解决两场耦合问题常用的方法还是迭代解法。本文同样采用迭代的方法。该程序先假定一初始应力场，计算围压后，通过式（3）、式（7）计算渗透系数，即从考虑固结压力的土-水特征曲线入手，建立应力与渗透系数的关系，从而计算渗流场，实现两场耦合计算。其中，应力计算采用弹塑性有限元计算程序，非饱和径流-渗流耦合计算采用第3节中所述方法。笔者在现有程序基础上，自编两场耦合有限元分析程序CPSS，用于耦合分析。其计算流程图如图 1 所示。

目前在工程实践中渗流对土坡稳定的影响分析一般采用传统的极限平衡法(LEM) ,强度折减有限元方法也已成功用于分析干土坡的稳定性。Gri ffiths D V和Lane P A[1 ]分析了稳定渗流对土坝稳定性影响的算例,但文中假定土体内的浸润线为直线且没有采用非饱和渗流分析以及没有考虑非饱和区的影响。Cai F和Ugai K[2 ]采用该技术并结合非饱和非稳定渗流有限元程序用于分析降雨作用下的土坡稳定性。Lane P A和Gri ffiths D V[3 ]采用强度折减技术用于确定水位缓降以及水位刚骤降瞬时这两种极端情况的土坡稳定性,但其中没有考虑非饱和非稳定渗流的影响。黄茂松和贾苍琴[4 ,5 ]利用强度折减有限元方法分析了饱和- 非饱和非稳定渗流条件下土坝的安全系数,并与极限平衡法进行初步的比较。韦立德[6～8 ]先后研制了饱和稳定渗流下的强度折减有限元程序、饱和- 非饱和稳定渗流下的强度折减有限元程序、考虑泡水强度降低的饱和稳定渗流下的强度折减有限元程序和考虑饱和稳定渗流下饱和区膨胀变形的强度折减有限元程序等,并把程序计算结果和常规极限平衡法的计算结果进行对比,并得到了一些有价值的结论。 饱和- 非饱和渗流计算基本方程和耦合原理 基本方程 假设地下水渗流在微段压力梯度上遵从Darcy定律。对于稳定饱和- 非饱和渗流场(无内源时)的简单情况计算归结为求解满足边界条件的拟调和方程:

第二章 土的渗透性和渗流问题

第二章 土的渗透性和渗流问题 第一节 概 述 土是多孔介质，其孔隙在空间互相连通。当饱和土体中两点之间存在能量差时，水就通过土体的孔隙从能量高的位置向能量低的位置流动。水在土体孔隙中流动的现象称为渗流；土具有被水等液体透过的性质称为土的渗透性。 土的渗透性是土的重要力学性质之一。在水利工程中，许多问题都与土的渗透性有关。渗透问题的研究主要包括以下几个方面： 1．渗流量问题。 例如对土坝坝身、坝基及渠道的渗漏水量的估算（图2－la 、b ），基坑开挖时的渗水量及排水量计算（图2－1C ），以及水井的供水量估算（图2－1d ）等。渗流量的大小将直接关系到这些工程的经济效益。 2．渗透变形（或称渗透破坏）问题。 流经土体的水流会对土颗粒和土体施加作用力，这一作用力称为渗透力。当渗透力过大时就会引起土颗粒或土体的移动，从而造成土工建筑物及地基产生渗透变形。渗透变形问题直接关系到建筑物的安全，它是水工建筑物和地基发生破坏的重要原因之一。由于渗透破坏而导致土石坝失事的数量占总失事工程数量的25%~30%。 3．渗流控制问题。 当渗流量和渗透变形不满足设计要求时，要采用工程措施加以控制，这一工作称为渗流控制。 渗流会造成水量损失而降低工程效益；会引起土体渗透变形，从而直接影响土工建筑物和地基的稳定与安全。因此，研究土的渗透规律、对渗流进行有效的控制和利用，是水利工程及土木工程有关领域中的一个非常重要的课题。 第二节 土的渗透性 一、土的渗透定律—达西定律 （一）渗流中的总水头与水力坡降 液体流动除了要满足连续原理外，还必须要满足液流的能量方程，即伯努里方程。在饱和土体渗透水流的研究中，常采用水头的概念来定义水体流动中的位能和动能。水头是指单位重量水体所具有的能量。按照伯努里方程，液流中一点的总水头h ，可用位置水头Z 、压力水头w u γ和流速水头g v 22 之和表示，即 1)-(2 22 g v u z h w ++=γ 式(2—1)中各项的物理意义均代表单位重量液体所具有的各种机械能，其量纲为长度。 对于流经土体中A 、B 二点渗流（图2－2），按照式（2－1），A 、B 两点的总水头可分别表示为： g v u z h g v u z h B w B B A w A A 22222 1++=++=γγ 且 h h h ?+=21

基于饱和-非饱和理论的尾矿库渗流场模拟分析

90 工业安全与环保2013年第39卷第10期 I ndus t r i al Saf et y a nd Envi w nnm [1t al Pr o t ect i on O ct ober 2013 基于饱和一非饱和理论的尾矿库渗流场模拟分析 朱毅川1 徐克2程镇2杨楠2 (1．中国地质大学工程学院武汉430074；2．湖北省安全科学技术研究中心武汉430070) 摘要在稳定性分析过程中，尾矿库力学参数选取困难，尾矿库边界条件复杂，需要能够较为精确预测尾矿坝浸润线位置和稳定性计算的方法。在饱和一非饱和渗流理论研究的基础上，对尾矿库渗流场的分布进行了数值模拟，通过对比研究，饱和一非饱和模型比饱和渗流模型更加符合尾矿库的渗流规律。研究了初期 坝透水性与干滩长度对尾矿库浸润线的影响，对比分析了影响尾矿库渗流的一些关键因素，以期能更好地解 决尾矿库渗流稳定性评估问题。 关键词尾矿库稳定性分析数值模拟饱和一非饱和理论 T he S i m ul a t i on of T ai l i ngs Seepage Fi el d B as ed on t he Saturat ed —-嘲turat ed T heor y Z H U Y i c hua nl X U K e2 C H EN G Zben2 Y A N G N an2 (1．Eng ／neer ／ng Facul t y ，C hi na U ni ver si t y of C ．eos ci ences W uba n 430074) A bst r act I n t he s t abil i t y anal y si s of t he r ai l i ng s ，t he sel e ct i on of t he m ec ha nl e al para m et e r s of t he t ai li ngs is usua l l y i na c-cu _rl l t e a nd SO i t i s ba dl y i n need of a m or e ac cu r at e m et h od t o f or ecast t aat ,峥dam s ee pa ge l l ne l ocat i on a nd s t abi l it y ca l cu- l at i on becaus e of t he bou ndar y condi t i ons of t he t a i l i n 铲．Base d o n t he sa t ura t e d —unsa t ura t e d s ee pa ge t heor y ，t he di st r i bu-t i m 0f t he t ai li ng pond s ee pa ge f i el d is si m ul a t e d i n t his paper ．．I l II oL 曲com parat i ve gt udy ，i t is f ound t hat sat ur at e d —t i n sa t ura t ed m odel c an bet t er agr e e t o t he t aaa 啦seepage r ul e t han sat ur at ed s ee pa ge m odel ；t he i nf l uence s of t he i ni ti al da m l m m eabi l i t y a nd dr y l e ngt h o n t he t aj l i I 骖i nf i l t r at i ng l i ne ar e r es ea r ched ．The paper ha s al so 舡l a 】咖脚80nl e key f act or s af- f ect i ng the 咖seepage i n or de r t o bet t er addr es s t he r ai li ngs s ee pa ge st abi l i t y asse ssm en t ．K ey W or ds t ai l ings pond s t abi l it y mla ：Iys i 8 num er i c al si m ul a t i on sa t ura t e d —unsa t ura t e d t heor y 0引盲尾矿库浸润线的高低是影响尾矿堆积坝稳定性的主要因素之一。尾矿水对尾矿坝体不仅会产生动水压力，降低坝体的稳定性，同时还会产生管涌、流沙和坝面沼泽化等危害，给尾矿坝的安全带来严重威胁。目前尾矿坝渗流分析有理论求解、物理模拟和数值模拟3种方法。理论求解浸润线的方法对工程进行概化后采用相应的近似公式求解，它只能适用于均质的简单的工程，对于多介质尾矿场，该法受限于理论条件；物理模拟方法保持了原型介质具有的连续性，并能较精确地模拟原型的复杂边界和内部结构。不足之处是制作模型的工作量大，费用高，而且费时费力；数值模拟法在我国20世纪六七十年代被应用到渗流计算中，并且随着计算机的高速发展，这种方法越来越广泛地被科研人员利用。 渗流数值计算中常用的数值方法有差分法、有 限单元法、边界元法等。过去，在渗流稳定性分析中，仅侧重于考虑饱和地下水对边坡稳定的影响，随着非饱和渗流理论的发展，人们越来越清楚地认识到边坡滑坡、泥石流等地质灾害与土体的非饱和渗 流密切相关。1962年，M i l l er 提出非饱和介质的渗透系数是含水量或压力水头的函数，达西定律同样适用于非饱和渗流…1，只不过要将达西定律的形式进行转变，渗透系数不能假定为常数而是与含水量或者是基质吸力有关的函数。1973年，N eum an 首先提出用有限元的方法解决土坝饱和一非饱和渗流问题的数值方法[2】2，采用伽辽金法对R i char ds 方程进行空间区域离散，用C r ank —N i coson 有限差分格式对时间域进行离散，并在拓展一维渗流有限元法之后提出了解决二维渗流问题的有限元法【3“J 。1997年，A kai 提出三维饱和一非饱和渗流的有限元 法[5】，Lai n 和Fr edl und 对饱和一非饱和土渗流问题做了较完整的论述[6】6，认为±体中水的流动是受压力水头大梯度和位置水头梯度共同控制的，并将非饱和土壤水运动理论和非饱和土固结理论结合，得到符合岩土工程使用习惯的饱和一非饱和渗流控制方程。 本文结合前人有关理论工作，在选取适当的土水特征曲线和渗透函 曲线模型的基础上，对尾矿 库渗流场的分布进行了值模拟，并对比分析了影

GeoStudio非饱和渗流关键问题探讨

GeoStudio 软件使用体会交流
汇报题目：非饱和渗流关键问题探讨
罗红明 中国科学院武汉岩土力学研究所 2010年6月

汇报提纲 汇报提纲
一、非饱和土特性 二、非饱和渗流基本理论 三、非饱和数值求解关键问题 四、降雨数值分析中几个问题 五、应用实例分析

一、非饱和土特性
1、土水特征曲线 基质吸力直接影响着非饱和土的力学性质。对于非饱和土体，其基质 吸力的大小则是土体含水量的函数。非饱和土的基质吸力随着含水量 的变化而变化，含水量和基质吸力的关系曲线称为土水特性曲线。 土水特性曲线对于研究非饱和土的物理力学特性至关重要。根据土水 特性曲线可以确定非饱和土的强度和渗透系数，甚至可以确定地下水 面以上水分分布。因此，研究含水量对非饱和土力学性质的影响，就 是研究非饱和土力学性质与基质吸力以及土水特性曲线的相互关系。 土水特性曲线的研究，起源于土壤学和土壤物理学。当时主要着重于 天然状态下表层土壤吸力的变化、土壤的持水特性及水分运动特征的 研究，基质吸力值一般小于100kPa。 近年来，由于非饱和土力学理论在斜坡稳定性评价以及降雨型滑坡预 测等方面的广泛应用，对非饱和土的土水特性曲线进行了更加深入的 研究。

一、非饱和土特性
2、土水特征曲线数学模型-1 Fredlund等通过对土体孔径分布曲线的研究，用统计分析理论推导出 适用于全吸力范围的任何土类的土水特征曲线的表达式：
θ 1 = F (? ) = C (? ) n m θs {Ln[e + (? / α ) ]}
C (? ) = 1 ? Ln(1 + ? / ? r ) Ln(1 + 10 6 / ? r )
θ 公式中，体积含水量的取值范围为： ∈ [0,θ s ] ，基质吸力的取值范围 ? ? 为： ∈ [0, ? max ] ， max 为土体含水量 θ =0时，所能达到的最大基质吸力。由 此可见，公式适用于全吸力范围的任何土类。但公式形式较为复杂， 给实际应用者带来诸多不便。

SciFEA软件模拟非饱和土壤渗流

SciFEA 饱和—非饱和土壤水分运动计算 计算机是现代科学技术的重要支柱，它不仅能带来巨大的经济效益，而且能带来深远的社会效益，对科学技术的发展起着巨大的推动作用。其在农业领域求解难于得到严格解析解的非线性偏微分方程的广泛应用成为现代化中计算问题的重要标志。国内外计算机已在农业科学实验、农业政策制定、资源普查与监测、植物和环境保护、土地规划及合理利用、农业信息处理与预测、自动化温室控制、农业气候与作物种植、土壤墒情监测与农田灌溉、农业机械化检测以及农药化肥生产控制等领域广迂使用。以节水灌溉为目标的现代灌溉管理需要科学地进行用水管理和定量描述并预测灌溉—蒸发条件下的土壤水分在土壤耕作层内的运动规律。土壤水分运动同时也是制约土壤溶质迁移的主要因素，对于干旱、半干旱的次生盐碱化威胁较严重的地区，这是至关重要的。因此，对于土壤水分的模拟研究不仅有助于现代节水灌溉管理，而且是土壤溶质迁移规律研究的基础。 随着科学的进步与发展，土壤水分运动的研究已由定性描述的形态学观点逐步发展成为定量研究阶段。土壤水分运动基本方程是一个非线性偏微分方程，在许多定解条件下，很难求得解析解。近年来．国内外许多学者都对求解非饱和土壤水分运动方程的数值方法进行了研究。 1问题描述 土壤水分运动基本方程 以土壤水势为因变量的一维饱和—非饱和等温、均质、各向同性土壤水分运动基本方程(Richards，1931)如下： ()()()C K Z t z z ?Ψ????Ψ=ΨΨ???????? (1) 其中：Ψ是土壤水势(cm)；()C Ψ是容水度；()K Ψ是非饱和导水率(cm/min)； Z是垂向坐标(cm)，向下为正；t 是时间(min)。并假设Ψ和K 均为土壤含水率θ的 单值函数。

土的渗透性和渗流问题

第二篇 土力学 第四章 土的渗透性和渗流问题 第一节 概述 土是由固体相的颗粒、孔隙中的液体和气体三相组成的，而土中的孔隙具有连续的性质，当土作为水土建筑物的地基或直接把它用作水土建筑物的材料时，水就会在水头差作用下从水位较高的一侧透过土体的孔隙流向水位较低的一侧。 渗透：在水头差作用下，水透过土体孔隙的现象 渗透性：土允许水透过的性能称为土的渗透性。 水在土体中渗透，一方面会造成水量损失，影响工程效益；另一方面将引起土体内部应力状态的变化，从而改变水土建筑物或地基的稳定条件，甚者还会酿成破坏事故。 此外，土的渗透性的强弱，对土体的固结、强度以及工程施工都有非常重要的影响。 本章将主要讨论水在土体中的渗透性及渗透规律，以及渗透力渗透变形等问题。 第二节 土的渗透性 一、土的渗透规律——达西定律 (一)渗流中的总水头与水力坡降 液体流动的连续性原理：（方程式） dw v dw v w w ?? =2 21 1 2211v w v w = 1 2 21w w v v = 表明：通过稳定总流任意过水断面的流量是相等的；或者说是稳定总流的过水断面的 平均流速与过水断面的面积成反比。 前提：流体是连续介质 流体是不可压缩的； 流体是稳定流，且流体不能通过流面流进或流出该元流。 理想重力的能量方程式（伯努利方程式1738年瑞士数学家应用动能定理推导出来的。） c g v r p Z =++22 饱和土体空隙中的渗透水流，也遵从伯努利方程，并用水头的概念来研究水体流动中 的位能和动能。 水头：实际上就是单位重量水体所具有的能量。 按照伯努利方程，液流中一点的总水头h ，可以用位置水头Z ，压力水头U/r w 和流速水 头V 2/2g 之和表示，即 g v r u Z h w 22 ++= 4-1 此方程式中各项的物理意义均代表单位重量液体所具有的各种机械能，而其量纲都是 长度。

三维非稳定饱和-非饱和渗流有限元法及其改进

三维非稳定饱和－非饱和渗流有限元法及其改进 江沆，沈振中，邱乾勇 河海大学水利水电工程学院，江苏南京（210098） E-mail ：jw01021514@https://www.wendangku.net/doc/db16710798.html, 。 摘 要：阐述了非稳定饱和－非饱和渗流的基本方程、定解条件和有限元迭代计算过程。针对有限元计算过程中经常出现的数值弥散现象，研究了其出现的规律和处理方法，提出了容水度的改进算法，并编制了有限元程序。最后详细模拟复合土工膜缺陷渗漏量室内物理模型试验，给出了试验过程中饱和－非饱和渗流场的变化过程。计算表明改进后的计算程序能有效地消除数值弥散现象，计算成果与试验结果基本一致。 关键词：非稳定饱和－非饱和渗流；容水度；有限元法 饱和－非饱和渗流在计算方法上把饱和－非饱和区作为一个整个流场研究，不需要进行自由面的调整，自由面仅是负压区正压区的分界面，不再是像饱和渗流场计算时流场的一个边界，避免了饱和渗流计算中试找自由面和选取给水度补给自由面边界的麻烦。但是，通常 k h ?和h θ?的函数关系是由离散的点线性插值得来，h θ?的线性插值易造成 h θ ??的非连续性，在计算中产生振荡，造成非饱和区孔隙负压计算结果分布规律较乱、与实际物理现象不符的现象，即所谓的数值弥散现象[ 3] [ 4 ]。为消除这一现象，黄康乐[ 3]采用变坐标的特征有限元法，朱学愚等[ 4 ]采用在Galerkin 有限元法的权函数中加上一个扰动的SUPG 有限元法，这些方法较好地改善了数值弥散现象，但是均较复杂。吴梦喜[ 5 ]提出了类似式（10）函数来计算容水度C ，能较好的解决这一问题，但表达式不严密，没有分析表达式中分母等于零时的情况。这里在文献[ 5 ]中对计算容水度C 的处理方法加以改进。 1. 非稳定饱和－非饱和土渗流定解问题 对于服从达西定律的非均质各向异性饱和－非饱和土，在忽略水的压缩性和孔隙气体对水流运动的影响情况下渗流基本微分方程为[ 1]： ()()()3s s c c ij r c i r c c s i j h h k k h k k h Q C h S x x t β????? +?=+????????????? （1） 式中：h 为总水头；c h 为压力水头；s ij k 为饱和渗透系数张量；r k 为相对透水率，是非饱和土的渗透系数与同一种土饱和时的渗透系数的比值，在非饱和区0