土壤水分运动参数的研究.docVIP

土壤水分运动参数的研究 　　摘要求解非饱和土壤水分运动方程进而预报非饱和土壤水分运动，必须首先获得土壤水分运动参数。参数的准确性决定于与这些参数相关的水分运动模型的可靠性。介绍了土壤水分入渗模型，概括了描述土壤水分运动的基本参数：土壤导水率（K）、土壤水分扩散率（D）、土壤比水容重（C）即水分特征曲线等。其中水分特征曲线被认为是土壤最基本的导水参数之一。 　　关键词土壤水分运动；基本参数；水分特征曲线 　　中图分类号 S152.7 文献标识码A文章编号 1007-5739（2009）04-0136-03 　　 　　1土壤水分入渗模型研究 　　1.1水分运动基本方程 　　Darcy（1856）通过饱和砂层的渗透试验，得出通量q 和水力梯度成正比，即达西定律：q=Ks ΔH/L，式中，L为渗流路径的直线长度，H为总水头，ΔH为渗流路径始末断面总水头差，ΔH/L是相应的水力梯度，Ks为饱和导水率。Richards（1931）将达西定律引入非饱和土壤水流动，表示为：q=-K（Ψm）?塄Ψ或q=-K（θ）?塄Ψ，式中，K（θ）为非饱和导水率，?塄Ψ为总水势梯度。它成为研究非饱和土壤水流动的基本定律。 　　达西定律是多孔介质中液体流动所应满足的运动方程，质量守恒是物质运动和变化普遍遵循的基本原理，将质量守恒原理具体应用在多孔介质中的液体流动即为连续方程。将土壤视为一种固相骨架不变形、各向同性的多由于滞后作用，基质势Ψm 和土壤含水量θ不是单值函数，土壤吸湿过程和脱湿过程不同，Richards 基本方程只用于吸湿和脱湿的单一过程。运用上述基本方程解决实际问题时，根据实际情况的不同及求解方便，基本方程可以有多种形式： 　　（1）以基质势Ψm为因变量的基本方程。非饱和土壤导水率K 和比水容量C 均可表示为土壤含水量θ的函数K（θ）（2）以土壤含水量θ为因变量的基本方程。非饱和土壤水分扩散率D（θ）定义为非饱和土壤导水率K（θ）和比水水分运动参数，用解析或数值方法对基本方程求解，就可得到土壤含水量θ或基质势Ψm 的空间分布及随时间的变化，即水分运动模型。 　　1.2Green-Ampt（1911）模型 　　Green-Ampt模型研究初始干燥土壤在薄层积水条件下入渗问题。基本假定是，入渗时存在明确的水平湿润锋面，同时具有固定不变的吸力Sf，土壤含水率θ的分布成阶梯状，湿润区为饱和含水量θs，湿润锋前为初始含水率θi。由达西定律得出地表处入渗率为： 　　式中，i为地表入渗率（cm/min），Sf为湿润锋处的土壤水吸力（cm），Zf为概化的湿润锋深度（cm），H 为积水衡原理，可得出累积入渗量I（t）为：I=（θs-θi）zf，对于入渗时间较长而Zf较大或H较小，公式②可简化为：i=Ks[1+(θs-θi）Sf/I］。 　　Green-Ampt 模型入渗公式简单，且有一定的物理模型基础，可应用于均质与非均质土壤或初始含水率不均匀的情况，均有较好结果。缺点是湿润锋处的土壤基质吸力Sf的确定较为困难，不能描述水分实际分布情况。 　　1.3Philip（1957）模型 　　Philip 模型中一维水分垂直运动基本方程取解的无穷级数形式为：I=St0.5+At+Bt0.5+……，式中，S 为吸湿率（cm/min），A为稳渗率（cm/min），各系数逐项递减且相邻2系数相差1~2个数量级，通常取前2项为常用的Philip 一维垂直入渗公式：I=St0.5+At，取第1项为忽略重力作用的Philip一维水平吸渗公式：I=St0.5，相应入渗率分别为：i= 1/2St0.5+A，i=1/2St0.5。 　　Philip 模型具有明确的物理意义，对于短历时的入渗情况较精确。缺点是只适用于均质土壤，在长历时的入渗情况下计算值与实际有较大偏差，对参数精度要求较高。 　　1.4Ghosh（1980）入渗公式 　　Ghosh 分析常用的Philip一维垂直入渗公式及Kostia-kov 经验公式得出，随着入渗时间的延长，两者预测的入渗过程与实际入渗过程逐渐偏离，对于长历时入渗，难以满足精度要求。因此，Ghosh 将两者综合考虑提出新的入渗公式为：I=atb=Kst，式中，a、b为常数，与土质及土壤含水量有关；Ks为饱和导水率。应用这一公式时，需要预先已知Ks值并根据土壤含水量确定出a、b值。Ghosh入渗公式能够很好地预测入渗过程，特别是对于长历时入渗，较其他公式具有更高的预测精度。但由于公式中含有3个参数，尤其当饱和导水率Ks 未知时，公式应用较繁琐且精度难以保证。 　　1.5其他入渗经验公式 　　在研究和解决实际问题计算过程中，不断提出的一些形式简单、方便实用的入渗经验公式得到了广泛应用。Kostiakov 入渗公式：i（t）=