用于模拟可变饱和多孔介质中水流的三维有限元模型

A Three-Dimensional Finite-Element Model for Simulating Water Flow in Variably Saturated Porous Media

PETER S. HUYAKORN 
GeoTrans, Inc., Herndon, Virginia 
EVERETT P. SPRINGER 
Los Alamos National Laboratory, Los Alamos, New Mexico 
VARUT GUVANASEN AND TERRY D. WADSWORTH 
GeoTrans, Inc., Herndon, Virginia 

**PETER S. HUYAKORN**  
GeoTrans，Inc.，弗吉尼亚州赫恩登  
**EVERETT P. SPRINGER**  
洛斯阿拉莫斯国家实验室，新墨西哥州洛斯阿拉莫斯  
**VARUT GUVANASEN** 和 **TERRY D. WADSWORTH**  
GeoTrans，Inc.，弗吉尼亚州赫恩登

本文提出了一种用于模拟可变饱和多孔介质中水流的三维有限元模型。该模型的制定是通用的，并能够适应与渗漏面、土壤表面的渗透或蒸发相关的复杂边界条件。在这个制定中包括了一个改进的Picard算法，旨在处理严重非线性的土壤含水量关系。该算法适用于矩形和三角柱元素。使用“影响系数”技术评估元素矩阵，避免了昂贵的数值积分。

对三维区域的空间离散化是使用垂直分片方法执行的，该方法旨在适应具有不规则边界、分层和/或横向不连续性的复杂几何形状。使用切片逐次超松弛方案实现矩阵解，该方案允许在小型小型计算机上有效处理相当大数量的节点未知数（数量在几千数量级）。

提供了六个例子来验证和演示所提出的有限元模型的效用。前四个例子涉及用作样本问题以进行基准测试的一维和二维流动问题。其余的例子涉及三维问题。

这些问题用于说明所提算法在涉及渗漏面和各向异性土壤媒体的三维情况中的性能。

**引言**

在过去的二十多年里，关于可变饱和多孔介质中水流的问题进行了深入的研究。对于涉及高度非线性土壤含水特性和与渗漏面、渗透和蒸发相关的大气边界条件的情况，这些问题难以解决。文献中报告的大多数模拟仅限于一维和二维流动情况。

到目前为止，只有很少几个三维数值模型得到了开发。 Freeze [1971]首次提出了一个三维有限差分模型，并用它研究了非均质、各向异性含水层中的饱和-非饱和瞬时流动。他的模型采用了一个隐式的Picard迭代过程，其中使用了垂直线逐次超松弛矩阵解决方案方案。尽管此方案可以容纳大量的节点未知数，但在涉及具有渗漏面和/或蒸发和渗透的稳态流动的情况下，它在处理中有严重限制。

自Freeze [1971]发表了他的开创性工作以来，对三维可变饱和流动模拟及相关模型开发进行了有限数量的研究。其中包括Seqol [1977]，Frind和Verge [1978]，Reisenauer等人 [1982]和Davis和Seqol [1985]的研究。在这些研究中使用了两种建模方法。

第一种方法[Reisenauer等人，1982]采用了集成有限差分法（IFDM），并采用混合显式-隐式的Picard迭代过程和矩阵求逆或点逐次弛豫方案[Narasirnhan和Witherspoon，1976]。此方法允许使用不规则网格，以提供比传统的有限差分方法更有效的表示具有复杂几何形状的区域。但是，为了计算质量通量，需要节点坐标、节点连接和节点面积数据。坐标、连接和面积必须满足某些正交约束，以确保质量守恒。

这些约束使得难以，如果不是不可能的话，对一般的曲线形状进行建模。报告了在不使用时间步进的情况下获得非振荡稳态解的困难[Kincaid等，1984，第3-37页]。此外，在涉及高度非线性土壤特性和渗漏面的问题中，使用混合显式-隐式时间步进与点逐次弛豫解决方案方案可能存在潜在问题。

第二种方法（由Segol [1977]和Frind和Verge [1978]提出）基于伽辽金有限元法，结合全隐式的Picard迭代过程和直接矩阵求解器。与IFDM方法相比，这种方法更灵活，适应曲线形状。

然而，即使对于涉及几千个节点的问题，其直接矩阵解决方案程序所需的核心存储和中央处理单元（CPU）时间也过多，其潜在的实用性受到限制。

在本文中，我们提出了一种改进的三维有限元解法技术，旨在缓解计算限制的问题。目前的论文是三篇中的第三篇。第一篇论文[Huyakorn等，1984]描述了用于二维可变饱和流问题的改进有限元制定。在其中，提出了一种无需数值积分即可进行元素矩阵计算的影响系数技术。

随后，在第二篇论文[Huyakorn等人，1986]中，提出了一种有效且可靠的质量平衡计算方案的制定。在本文中，将呈现一种创新的三维有限元技术，用于模拟可变饱和多孔介质中的水流。第一篇论文中提出的Picard算法得到了增强和扩展，以适应具有渗漏面和与渗透和蒸发相关的大气边界条件的严重非线性三维问题。影响系数矩阵方法与切片逐次超松弛（SSOR）矩阵解决方案方案相结合，已适应三维问题。这种方法允许以经济的方式分析复杂的三维现场问题。提供了六个例子来验证所提出的有限元模型，并展示其多功能性和效率。

GOVERNING EQUATION AND BOUNDARY CONDITIONS

Adopting the assumptions from Huyakorn et al [1984], the governing equation for three-dimensional flow of water in a variably saturated porous medium can be written in the form 

\frac{\partial }{\partial x_i}[K_{i,j}(\frac{\partial \Psi}{\partial x_j}+e_j)]=\eta \frac{\partial  \Psi}{\partial t}-q\qquad(1)