S. Panday a（1994）

A three-dimensional multiphase flow model for assessing NAPL contamination in porous and fractured media, 2. Porous medium simulation examples

S. Panday a, Y.S. Wu a, P.S. Huyakorn a, E.P. Springer b

**标题：用于评估多孔和裂隙介质中NAPL污染的三维多相流模型，2. 多孔介质模拟示例**

**S. Panday a, Y.S. Wu a, P.S. Huyakorn a, E.P. Springer b**

**摘要：**
本文讨论了由Huyakorn等人提出的三维（3-D）多相流模型（本刊第I部分）的验证和应用，用于评估由于地下释放非水相液体（NAPL）而导致的污染。注意力集中在涉及一维、二维和三维流经多孔介质的情况。对高度非线性的现场条件测试了模型公式和数值方案。评估了对某些模拟场景的各种简化的效用和准确性。为了演示目的，包含了五个模拟示例。第一个示例验证了模型在垂直流动方面的性能，并比较了完全三相和被动空气相形式的性能。空气相边界条件对模拟结果有相当大的影响。第二个示例验证了模型对涉及LNAPL和DNAPL迁移的横截面分析的适用性。比较了不同网格长宽比的有限差分（5点）和有限元（9点）空间逼近。除非纠正，否则发现负传导率条件对有限元解决方案产生不良影响。第三个示例对油藏的5点式水驱处理进行了模型验证。对于有限差分方案，网格方向的敏感性得到了注意。第四个示例展示了模型在表面源处表征LNAPL和DNAPL的三维迁移的实用性。最后一个示例呈现了对空气冲刷的建模研究。研究了影响空气冲刷系统性能的关键参数。总体而言，建模结果表明在保水土壤中冲刷更为有效，并且在某些各向异性条件下可以实现更大的冲刷影响半径。

**1. 引言：**
对非水相（NAPL）污染物的行为和地下迁移进行建模需要适当地表示重要的物理机制、高效和稳健的数值解法，以及了解可能针对许多实际情况进行的各种简化假设。本研究的第一部分（本刊1994年Huyakorn等）提供了多相流问题的数学公式和实际场模拟所需的数值解法技术。已讨论了对流问题的各种简化（例如，用于三相系统的被动空气相形式，或者用于空气冲刷或蒸发过程的被动液相形式）。在模型中已实施了多种空间逼近选项，包括7点有限差分、27点有限元和11点（混合）三维（3-D）格子，以及它们的二维（2-D）对应物。本文的目的是验证模型及其数值方案，提供代码验证的基准问题，检查各种简化假设的有效性，展示数值方案的准确性，并讨论模型在各种场景中的适用性。

对于地下NAPL污染的文献记录的数字模型和建模工作的调查表明，很少有出版物详细说明各种简化假设和数值技术如何影响模拟结果。这是因为大多数先前的模型，如果不是所有的话，都没有结合灵活的公式，允许不仅使用最一般的三相流动方程，还可以使用其减少形式，适用于简化假设合理的特殊情况。目前，对于描绘三相流动行为和饱和度分布的被动空气相形式的表现，仍然缺乏评估，尽管如果空气相动态重要，被动空气假设是不合理的。还需要研究各种空间逼近方案在场问题中的表现。例如，即使对于在石油文献中常描述的典型的5点井流场景，也没有对各种空间逼近对网格方向效应的敏感性进行检查（Yason~k和McCracken，1979；Thomas，1982；Willhite，1986）。对于多相流动情况，在某些网格长宽比的负传导率条件的影响也需要进行研究。对于各向同性矩形棱柱元素，网格长宽比为x/y：1（Forsyth，1988）可能导致不可接受的负传导率条件，严重影响收敛并产生物理上不正确的结果。如果施加了严格的网格限制以避免负传导率，实际的三维场模拟将变得不切实际。当前的研究研究了这些效应，以更好地了解各种常用的数值方案的准确性和可靠性。

本文包括五个模拟示例。第一个示例用于验证代码及其解决方案方案与先前

发表的结果相比的性能。将被动空气相形式与完全三相形式进行比较。