三维渗流问题的可变饱和有限元模型

**作者：Kalf, Frans Robert Peter 发表日期：1988**

**摘要**

通过使用四面体有限元制定一个可变饱和地下水流模型，可以实现对三维渗流问题的模拟。这样的模型克服了基于移动有限元网格或与传统的三维分层模型相关的经验性脱饱和方案的自由表面模型所面临的问题。一种称为TETFLOV的模型已经得到开发，能够模拟可变饱和流条件下的渗流面流动、渗透和蒸发、水源和汇以及非达西流动。该模型经过与现有的稳态和瞬态流动解的比较进行了验证。该模型被应用于通过和从堤坝的渗漏；土柱渗透；从溪流的瞬时渗流和非达西径向流向坑中的验证通过了罐实验和独立的数值解。该模型还可以应用于复杂的露天矿高墙的水平钻孔降水。该模型使用迭代和直接求解算法。

实施多项式预处理共轭梯度求解器提供了诸如快速收敛、速度、较小的计算机内存需求和节点编号的灵活性等优势。通过相同的元素还实现了可溶性饱和流的子问题的溶质输送子例程的制定。在溶质运输模拟器上进行的数值实验表明，对于支配性分散问题，可以获得良好的准确性，但在支配性的仿真中，数值分散和超调现象是明显的。使用截断消除过程（TCP）的高阶逼近在应用于三维可变速度流场时产生了不佳的结果。可以设计一种经验算法，模拟一些报告的顺风方案抑制超调的效果，但以牺牲数值分散为代价。溶质运输子模型已通过与文献中报告的一些数值解的比较进行了验证。TETFLOV模型已得到开发，并且可以在升级的PC兼容微型计算机系统上运行。

**第1章 引言**

**1.1 问题描述**

使用数值方法模拟地下水流已经成为近年来的常态。重点是使用有限差分和有限元方法解决二维和准三维问题。

解决完全三维问题更为艰巨，不仅因为增加了复杂性，而且因为大多数计算机的存储限制和速度。这些限制到目前为止尚未完全克服；然而，最近微型计算机系统的速度和可直接寻址的内存有了显著提高，而且价格相对较低。这些发展现在使得能够以低成本在这些标量计算机上解决实际的三维问题成为可能。

许多工程和矿业地下水渗流问题本质上是三维的。虽然通常有人争辩说缺乏数据不足以证明需要进行三维分析，但经验表明，对问题维度的限制往往是成功建模的更大障碍，而不是缺乏合适数据。

传统上，使用分析技术或使用二维数值程序来近似三维流问题。对于非受限流问题，在正确表示是一个饱和和非饱和多孔介质的水位面之间的移动边界时出现了困难。通常已经使用了三种数值方法来试图克服这个问题。在

横截面问题中，一种方法是使用“移动”有限元网格，其中网格被扩展或压缩以模拟水位的上升或下降。

这种方法在计算上是有效的，并且已成功地应用于许多地下水流问题（Neuman和Vitherspoon，1970a，b；France，1974；Chueng和Skjolingstad，1974；Taylor和Huyakorn，1978；Taylor，1980；Issacs，1979；Desai等，1983；Dudgeon，1985）。

在三维自由表面模型中，常规方法是在各个层中定义一个根据水位水平变化的传导率项。当发生脱饱和即水位下降到每个层以下时，这种假设就会出现特殊问题（lcDonald和Harbaugh 1984）。总的来说，复杂的渗流面不能被准确地建模，而且当水位在上层内时，这种方法是最令人满意的。

第三种方法是将问题制定为饱和/非饱和流系统，其中水位在固定的网格系统中简单地成为零压力线。虽然计算工作量增加了，但增加的灵活性使其能够应用于更复杂的流系统。这个解决方案过程既可以使用有限差分方法也可以使用有限元方法（Freeze，1971；Neuman等，1972，1974，1973，1977；Frind和Verge，1978；Cooley，1971，1983；Blanford，1984；Lam和Fredlund，1984）。

由于由于其提供非线性数量的元素平均，有限元方法或集成有限差分方法通常被用于这个高度非线性的流问题，它通常提供更好的整体稳定性（Cooley 1983）。

非受限的渗透问题提出了一个更大的困难，因为它们通常要求模拟渗透面或渗透面。由于这些面的高度事先是未知的，因此需要某种类型的迭代解决方案过程。

此外，对于许多应用程序，还需要渗透和蒸发边界。

对于孔隙速度很高的流动体系，速度和水力梯度之间的关系变得非线性，因此流动是非达西流动的。

此外，所描述的流动系统通常需要解决污染物流问题。

作者不知道有任何已发表的模型，该模型使用可变饱和方法，允许同时在三维流动体系中考虑上述各种流动条件。

除了上述之外，还进行了工作，以包括一种高效的迭代解决算法，作为文本中所述的直接解决方法的替代方法。直接技术可能需要大量的计算机存储器，因为对于具有大量节点的问题，获得了大的带宽。此外，随着向量和并行处理的出现，目标是提供一种可以最终利用这些技术的解决方案过程。结果，已包括多项式预处理共轭梯度解算算法。

综上所述，本论文涵盖了以下新的研究领域：1. 使用四面体有限元制定一个完全的3-D可变饱和渗流面模型。

2. 通过使用特定四面体元素方向创建的解决偏差的演示。

3. 通过可变饱和流和问题规模创建的缩放效应的演示。

4. 对一个已知的渗漏算法进行扩展和修改，使其在三维中处理水平钻孔渗漏问题。

5. 将该模型应用于模拟复杂地质渗漏面中的多个水平钻孔排水。

6. 为三维可变饱和流场开发非达西流动代码。

7. 对应用于可变饱和1 3-D流的标量多项式预处理共轭梯度解算算法的开发和性能分析。

8. 使用四面体有限元为TETFLOV模型开发饱和的溶质运输例程。