《一种用于预测非饱和土壤水力导度的封闭形式方程》

A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils1

M. TH. VAN GENUCHTEN2

摘要：本文描述了一种新的、相对简单的土壤-含水量-压力头曲线8(h)的方程。该方程的特定形式使得可以在N.T. Burdine或Y. Mualem的预测导度模型中代入，从而导出相对水力导度Kr的封闭形式解析表达式。所得到的Kr(h)的表达式包含三个独立参数，可以通过将提出的土壤-含水保留模型与实验数据进行拟合来获得。利用基于Mualem理论的封闭形式解析表达式得到的结果与具有广泛水力性质范围的五种土壤的观察到的水力导度数据进行了比较。在五个案例中，非饱和水力导度的预测在四个案例中得到了良好的结果。发现在低含水量下合理描述土壤-含水保留曲线对于准确预测非饱和水力导度至关重要。

附加索引词：土壤-水扩散性，土壤-含水保留曲线。

van Genuchten, M. Th. 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44:892-898.

数值模型用于模拟非饱和带中的流体流动和质量输运在最近几年变得越来越受欢迎。

近期的文献确实表明，人们在开发这类模型方面付出了很多努力（Reeves和Duguid，1975；Segol，1976；Vauclin等，1979）。

不幸的是，似乎充分表征模拟系统的能力没有跟上数值和建模专业知识的发展。导致非饱和流理论成功应用于实际领域问题的最重要因素之一可能是关于进入支配传递方程的参数的信息不足。由于在现场具有广泛变异性，并且测量该参数是耗时且昂贵的，因此可靠估计非饱和水力导度尤其难以获得。由于这些原因，一些研究人员使用模型从更容易测量的土壤-水保留曲线计算非饱和导度。其中，Millington-Quirk方法（Millington和Quirk，1961）在这些模型中非常受欢迎，其各种形式已在一些研究中取得了一定成功（参见Jackson等，1965；Jackson，1972；Green和Corey，1971；Bruce，1972）。不幸的是，该方法的缺点是产生表格结果，例如，当应用于多维非齐次流模型中的非饱和流时，使用起来相当繁琐。

还已经开发了预测非饱和水力导度的封闭形式解析表达式。例如，Brooks和Corey（1964）以及Jeppson（1974）分别使用了基于Burdine理论（Burdine，1953）的导度的分析表达式。Brooks和Corey（1964, 1966）使用其方程获得了相当准确的预测，尽管在一些负的压力头值（这点通常称为起泡压力）存在土壤-水保留曲线和非饱和水力导度曲线的斜率的不连续性。这种不连续有时会阻碍在数值饱和-非饱和流问题中的迅速收敛。此外，似乎基于Brooks和Corey方程的预测相对于（修改的）Millington-Quirk方法的各种形式来说略显不准确。

最近，Mualem（1976a）从土壤-水保留曲线和饱和时导度的知识中推导出了一个用于预测水力导度的新模型。Mualem的推导导致了非饱和水力导度的一个简单积分公式，该公式使得可以推导出封闭形式的解析表达式，只要有适当的土壤-水保留曲线的方程可用。本文的目的是使用土壤-水保留曲线的方程（该方程既连续又具有连续的斜率）推导出这样的表达式。由Mualem理论提出的封闭形式解析表达式通常包含三个独立参数，可以通过将提出的土壤-水保留曲线与实验数据进行匹配来获得。将基于Mualem理论的封闭形式方程得到的结果与具有广泛水力性质的几种土壤的观测数据进行比较。

理论

基于Mualem模型的方程

Mualem（1976a）导出了以下方程，用于根据土壤-水保留曲线的知识来预测相对水力导度（Kr）：

其中\[h_t\]是压力头，这里被给出为无量纲含

水量的函数：

在这个方程中，s和r分别表示土壤-水含量（θ）的饱和和残余值。为了解Eq. [1]，需要一个将无量纲含水量与压力头关联的表达式。在本研究中采用的一类吸引人的Θ(Ψ)-函数由以下一般方程给出：

其中α、n和m是尚未确定的参数。为了简化后续的符号表示，假设Eq. [3]中的\[h_t\]是正的。

方程\[3\]在m=1的情况下已经成功地在许多研究中用于描述土壤-水保留数据（Ahuja和Swartzendraber，1972；Endelman等，1974；Haverkamp等，1977）。

根据Eq. [2]和[3]得到的典型的ah曲线如图1所示。注意到得到了一个几乎对称的“S”形曲线，并且当\(\theta\)接近其饱和和残余值时，斜率（d\(\theta\)/dh）变为零。