2_非饱和水流运动基本方程精选.pptxVIP

非饱和土壤水流的基本方程§1非饱和土壤水流的达西定律§2土壤水分运动参数及其测定方法§3非饱和土壤水运动的基本方程§4土壤水分运动基本方程的定解条件1.非饱和土壤水流的达西定律非饱和土壤水分运动和饱和土壤水分运动一样，水分从水势高处向水势低处运动。一般认为，适用于饱和水流动的达西定律在很多情况下也同样适用于非饱和土壤水分流动。1931年，Richards最早将达西定律引入非饱和土壤水流动。非饱和土壤水分流动的达西定律:饱和土壤水分流动的达西定律:达西定律的推导流体中由动量和连续方程可以推导出渗流中与渗透率成反比速度v由流量/孔隙度替换忽略可压缩性 得的处理渗流边界上速度不为达西定律有速度上限，速度太高时需要考虑惯性和湍流效应。而对于多孔介质，骨架会阻止流体运动，这一项的存在就不合理了。去除项得对于稳态得达西定律ks，渗透系数=const.非饱和流饱和流水势组成：ψ=ψg+ψpψ=ψg+ψm流动准则：高低土水势总水头高低k(θ)为土壤含水率的函数ks≥k(θ)θ↗⇒k(θ)↗ψg:ψp:ψm=0相对参考平面之高度至地下水面的高度ψp=0ψm:ψm取决于土壤的干湿程度K:2.土壤水分运动参数及其测定方法非饱和/饱和水力传导度容水度土壤水分扩散度是反映土壤水分在水势梯度作用下流动的性能。一般在饱和土壤中导水率称为渗透系数，为常量。定义：在水势梯度作用下，单位断面面积上流过的水流通量。非饱和水力传导度的概念及特征在非饱和土壤中，导水率是负压或含水率的函数，随着含水率降低而减小。非饱和水力传导度及其测定K(h)K(θ)K(θ)\K(h)多种近似vanGenuchten-Mualem:BrooksandCorey(BC):非饱和水力传导度的测定在水平土柱两端有多孔板，分别由平水箱保持一定水位，使其负压为h1和h2，在梯度作用下，土柱中土壤水从l端向2端运移。土壤水通量q可由l端补给量或2端溢出量测得，两者相等时，水流处于稳定状态。非饱和土壤水力传导度可由达西定律求得。计算的k是平均的k试样中各点的是不相同的非饱和达西实验在不同的平均负压(吸力)值下，通量与负压梯度成正比，两者呈直线关系，但其斜率(即水力传导度)随平均负压而变。容水度(或比水容量)单位基膜势(负压值)变化所引起土壤含水率的变化，一般称为容水度或比水容量（C)，可以下式表示：表示在单位压力水头降低时自单位体积土壤中所释放出来的水的体积，它与饱和土壤的给水度相似。用测水分特征曲线的方法来测定土壤水分扩散度D土壤水分扩散度为单位含水率梯度下，通过单位面积的土壤水流量，其值为土壤含水率的函数，即土壤水分扩散度与土壤的关系可用以下经验公式表示直角坐标非饱和水分基本方程3.非饱和土壤水运动的基本方程理论基础：达西定律质量守恒定律（水流连续原理）直角坐标非饱和水分基本方程基本方程推导同理(y+y)、(z+z)面流速为：取微分单元体，体积：xyz设沿x、y和z方向流速分别为：vx，vy，vz，则(x+x)面流速为vxvzvy设六面体土壤含水量为，则t内六面体内土壤水质量变化量为：根据质量守恒原理有即——非饱和土壤水运动基本方程，可简写为：根据达西定律有:将上式代入假定土壤各向同性，则有：基本方程的不同形式用基质势h为变量的基本方程对于非饱和土壤水，总水头H由负压水头h和重力水头z组成：c（h）表示比水容量（也称容水度）令则有对上式求偏导,则有故用基质势h为变量的基本方程剖面二维：垂向一维：以基质势h为变量的基本方程，最突出的优点是适用于饱和-非饱和问题的求解，也可用于分层土壤的水分运动的计算，但非饱和土壤的导水率和容水度受滞后影响较大，计算中参数选取不当会造成较大误差。基本方程的不同形式用含水量θ为变量的基本方程D（θ）为土壤水的扩散率，令则有上式中代入上式有用含水量θ为变量的基本方程剖面二维：垂向一维：以含水量θ为变量的基本方程常用于求解均质土层或非饱和流问题，但不适宜层状土壤或求解饱和-非饱和问题。这里达西定律其实就是流体的动量定理的转变，当不考虑温度变化时与连续方程耦合。若要考虑能量方程，就得给出温度条件，再和温度耦合。4.土壤水分运动方程的定解条件初始条件（t＝０）边界条件以垂向一维流动为例：一类边界（变量已知边界）：在一维垂向土壤水分运动中，一类边界的情况发生在：地表形成积水时；地表含水率达到饱和含水率；当强烈蒸发时，表土达到风干土含水率。二类边界条件(边界上水流通量已知)在一维垂向土壤水分运动中，这种情况常发生在降雨、灌水入渗或蒸发强度已知的边界上。在降雨或灌水入渗时，(t)为负值，在蒸发时(t)为正值。在不透水边界和无蒸发入渗的边界，(t)=0，则三类边界条件：相当于水流通量随边界上的变量(含水率或压力)值而变化的情况在土壤蒸发强度为表土含水率