土壤水动力学SWD土壤水运动的基本理论与方法要点.pptVIP

第二类边界条件（水流通量已知边界Neuman条件） q=-K(θ) ▽(ψm±z) 垂向一维： 举例： 通量为零的情况：如不透水边界、无蒸发入渗的边界 通量已知的情况：降雨、灌溉、蒸发强度已知 第三类边界条件（水流通量随边界上的变量变化而变化的情况 ） 举例： 当土壤蒸发强度为表土含水率（基质势）函数的情况 数值模拟中常遇到的其他边界条件，一般为以上边界条件的组合： 渗流边界 （seepage face） 如大坝表面，当非饱和时为二类（零通量）边界，当饱和时为一类（已知水头or含水率）边界 自由排水（重力排水）边界 （free drainage ） 假设垂向水势梯度为1，一般用于土壤深层，注意此时通量为含水率（基质势）函数，非常值。 3.4 土壤水分运动的通量法简介 直接利用Darcy定律和连续方程分析土壤水分运动特性 一维垂直运动： z*～z积分： 通量法：根据z*处通量及含水率变化估算其它深度通量 零通量面法：确定零通量面位置 表面通量法：估计地表蒸发/入渗通量 定位通量法：根据实测水势差、导水率估算某一位置通量 零通量面法： 零通量面--当水势梯度为0，该处通量为0，则该处为零通量面 分类： 单一聚合型零通量面 单一发散型零通量面 多个零通量面 零通量面法： 由于零通量面为已知通量（零）断面，若t1和t2时段内零通量位置不变，则根据两时刻的土壤含水率观测值，可计算出时段内任一断面处流过的土壤水通量。 ZFP ?(z,t) ? z ? ?(z,t1) ?(z,t2) ?s Qs= ？ Qg=? 0 Z0 H ZFP1 ?(z,t) ? z ? ?(z,t1) ?(z,t2) ?s Qs= ？ Qg=? 0 Z01 H ZFP2 Z02 表面通量法： 以地表处的入渗量（Qs0）或蒸发量（Qs0）作为已知条件，根据实测时段始末土壤含水率剖面数据，推求土壤其他断面处水流通量。 入渗、蒸发量需进行估算或实测 定位通量法： 在土壤剖面中选定一个合适的位置，上下安装两支负压计，监测两点的基质势。根据达西定律计算定位点处通量 需要确定非饱和导水率~基质势关系 适合于土壤含水率变化不大的情况（如靠近潜水面） 土壤水动力学 第3章 土壤水运动的基本理论与方法 土壤水动力学 (Soil Water Dynamics)第3章 土壤水运动的基本理论与方法 (Theory Analyzing Method of Soil Water Movement) 毛 晓 敏 中国农业大学水利与土木工程学院 Tel:Email: maoxiaomin@tsinghua.org.cn 第3章 土壤水运动的基本理论与方法 Darcy定律 土壤水运动基本方程 基本方程的定解条件 土壤水分运动的通量法简介 3.1 Darcy定律 地下水流动的Darcy定律 土壤水流动的Darcy定律 土壤导水率 3.1.1 地下水流动的Darcy定律 1856年：Darcy根据饱和砂土的渗透试验，得出了水流通量q与水力梯度成正比的结论，称之为Darcy定律： q=KsΔH/L L－渗流路径的直线长度 H－水头； ΔH －水头差 ΔH/L－水力梯度 Ks －饱和导水率（Hydraulic conductivity,也称渗透系数、水力传导度） ，表示多孔介质透水性能 Darcy定律的微分形式： 三维：q=-Ks grad H=-Ks▽H Hamilton (Nabla)算子： 梯度： Darcy定律的适用范围： 小Re数层流：与粘滞力相比，惯性力作用可以忽略不计。在紊流状态下，通量与水势梯度呈非线性关系 对于颗粒极细的土壤：克服一定的初始水头差才能发生流动 一般情况下，Darcy定律有效 通量 水势梯度 饱和导水率Ks ： 综合反映了多孔介质对流体流动的阻碍作用 多孔介质的基质特征：质地、结构… 流体物理性质：粘滞性、密度… 实验室测定： 现场测定： 双环入渗试验 Guelph渗透仪 抽水试验 1.5m 1.2m 土柱直径20cm， 测得水流通量为0.5L/h, 画出沿土柱水头变化曲线，并求土壤饱和导水率？ 1.5m 1.2m 若土柱分两层，上层Ks为6cm/h,厚1m；下层土壤Ks为0.6cm/h,厚0.2m. 求土壤水流通量，并画出沿土柱水头变化曲线 。 Ks=6 cm/h Ks=0.6 cm/h 1 m 3.1.2 土壤水流动的Darcy定律 1907年：E