三维饱和非饱和降雨入渗渗流场分析.pdfVIP

水 利 学 报 2003 年12 月 SHUILI XUEBAO 第12 期 文章编号：0559-9350 (2003) 12-0066-06 三维饱和-非饱和降雨入渗渗流场分析 １ ２ １ 朱岳明 ，龚道勇 ，罗平平 (1.河海大学 水利水电工程学院，江苏 南京 210098； 2.长江水利委员会 规划设计研究院，湖北 武汉 430010) 摘要：用饱和-非饱和达西渗流理论对降雨入渗的复杂渗流场问题用有限单元法求解，在理论上进行了深入的研究。 文中特别提出了降雨入渗边界条件及非饱和逸出面边界条件(包括降雨中间停雨时期边界条件)在渗流场求解中的正 确处理方法，提出和定义了描述土体渗流特性的极限含水率和蒸逸系数的两个新概念。两个算例的数值试验也进一 步表明，文中所用的降雨入渗渗流场分析方法结果合理，实用性强。 关键词：饱和-非饱和渗流；降雨入渗；非饱和逸出面；极限含水率；环境湿度；蒸逸系数 中图分类号：TV139.14 文献标识码：A 边坡中的渗透水流及其渗透荷载是影响边坡变形和稳定的重要因素之一，而降雨，特别是历时长、强 度大的暴雨更是直接造成边坡失稳破坏的主要因素。如 1988 年龙羊峡水电站底孔首次放水，水流消能形 成雾雨，使下游右岸山坡发生滑坡，险些阻塞河道。随着人们对降雨过程中边坡失稳现象的重视，对降雨 入渗问题的研究越来越深入，取得了许多有益的成果，但因问题的复杂性，这一课题仍然是目前学术界研 究的重点课题之一。 降雨入渗是一个涉及两相流的渗透过程，即水在下渗过程中驱替空气的过程。目前，降雨入渗补给地 下水量的计算方法，国内通常多采用降雨入渗系数法，或采用统计的方法，求得降雨量或降雨补给量与地 下水位上升之间的关系，或通过成因分析探讨降雨补给与各因素之间的相互关系。本文根据饱和-非饱和 地下水运动微分方程和定解条件，模拟降雨入渗问题，并阐述降雨入渗边界条件的正确处理方法，特别是 对饱和-非饱和渗流入渗或逸出面的数值处理方法具有新的见解，使三维饱和-非饱和降雨入渗的复杂渗流 场问题得到了更加严密的有限单元法的数值求解。 1 饱和-非饱和降雨入渗渗流场分析的有限单元法 1.1 数学模型 认为降雨入渗补给地下水的过程即为非饱和土壤水连续运移的过程。式(1)为非稳定饱和 ［1～3］ -非饱和地下水运动的控制微分方程，式(2)～式(7)为它的初始条件和边值条件 。 ∂  ∂h  ∂h s ( ) c s ( ) [ ( ) β ] c k k h +k 3 k h +S C h + S (1) ∂xi  ij r c ∂xj i r c  c s ∂t 式中：h 为压力水头；t 为时间；x 为坐标；i=1，2，3；S 为源汇项；C 为容水度，C= θ/h , θ为体积含水 c i c 收稿日期：2002-08-06 基金项目：教育部年轻骨干教师资助项目(教技司65-10)；南京市人事局基金资助项目(宁人字-2000-4)；科技部中奥合作 资助项目(2001-Ⅵ-B-1)；江苏省科技厅资助项目(BZ2003046) 作者简介：朱岳明(1958-)，男，浙江长兴人，博士，教授，从事水工结构工程和岩土工程的渗控、温控、稳控和反问题等 研究。 66