R. K. Rorncu and B. Noctingcr （1995）

Calculation of internodal transmissivities in finite difference models of flow in heterogeneous porous media

**标题：有限差分模型中非均质多孔介质中节间传导率的计算**

**R. K. Rorncu 和 B. Noctingcr**

**摘要：**
我们对非均质多孔介质的等效导电率的有限差分计算进行了分析和数值研究。评估有限差分节间传导率的传统谐波方案在数值结果中产生系统性偏差，除非使用极细的网格。

为了量化这种影响，我们开发了一种解析方法，将等效数值导电率进行了导电度方差的幂级数展开。这导致将数值答案表达为网格块大小的函数。计算确认存在明显的偏差和非常缓慢的收敛。我们提出了一种简单的纠正方法，非常适合上溯。

对比使用更对比度较大的非均匀介质进行的数值实验显示了类似的结果。这使得可以使用更粗的网格，从而在数值模拟方法中获得可观的加速。

**引言：**
自然含水层或油藏表现出广泛的导电度变化，在流体流动计算中必须考虑这些变化。在这些非均质情况下的分析处理的复杂性只允许得出很少的确切分析结果，大多数研究和应用是通过数值方法进行的。

其中一个相关问题是在非均匀介质中放大水力导电度，因为由于计算原因，网格块的数量有限。与含水层的典型尺寸相比，即使在高分辨率模型中，一个基本网格块的尺寸也在10-100米的数量级。因此，它很可能是非均匀的，我们必须设计一个上溯过程，以使用实验室和地质数据为块分配导电度。这使我们能够计算非均质多孔介质的有效或等效导电度。我们在论文中进一步讨论了使用能量标准和周期性假设的等效导电度的定义。

除了在许多应用中具有重要实际意义之外，等效导电度还可以作为任何数值模拟器的准确性的代表性度量来使用，而不依赖于任何上溯应用。在非特定的应用中，等效导电度可能是最能总结结果的单一值，提供了所有计算出的局部值的一种规范，并具有明确定义的物理解释。它必须被正确计算，因为准确的能量平衡是检查任何数值程序的正确性的必要条件（但不是充分条件）。

我们预计一些传统方案在履行这一基本标准时系统性失败，这表明数值模拟中存在一些偏差。对于标准有限差分方案而言是真的，在这种方案中，通过使用块导电度的谐波平均值来计算节间传导率。正如我们将看到的那样，这导致了等效导电度的低估。通常使用精细网格来克服这个问题，但仍然没有实际的算法收敛性的实用估计器。对于特定的应用，有一个简单的公式，可以根据网格间距、对数传导度方差和相关长度的函数给出所选数值模型的准确性的估计器，将是有用的。这是本文的一个目标。

我们对非均质多孔介质的等效导电度的数值计算进行了分析和研究。更确切地说，我们关注了由给定算法产生的等效导电度结果随网格间距的收敛性。由于在最一般的情况下，该问题似乎难以解决，我们通过在传导度对比度或方差的幂级数中进行展开来进行攻击。尽管严格来说只对小对比度有效，但数值实验表明，当考虑到介质具有较强的导电度异质性时，我们的结果是相当稳健的。我们提出了一种允许在粗网格中正确评估等效导电度的最佳程序。我们只考虑中心有限差分算法，它作为数值流体流动模拟的传统参考，但相同的方法可以逐步用于任何其他离散化方案。

在有限差分模型中，评估节间传导率似乎非常明确但脆弱：在非均匀介质中，它涉及对介质最小尺度结构的一些隐藏假设，这可能导致结果中的重要偏差。大多数模拟器使用谐波平均规则，理论上仅适用于一维流动。已经提出并测试了其他规则，如算术平均、几何平均，甚至更一般的功率平均公式，但它们的理论依据仍然是经验性的。在这方面，必须区分一些关于传导度输入数据的确切含义以及小尺度上的流场模式的重要假设。本工作的一个重要目标是通过提供客观标

准和一个明确的框架来为这个讨论做出贡献。

另一个重要的概念是网格细化。在均质情况下没有问题，但在非均匀介质中这是一个相当微妙的问题。我们可以有比异质性测量支持小的网格块，因此我们必须对最小尺度的异质性模式做出假设。

本文的组织如下：首先，我们介绍了这里描述的非均匀介质的方式。然后，我们讨论文献中找到的在有限差分模型中计算节间传导率的不同方案。

我们展示了在标准情况下传统方案导致不正确的结果。在介绍我们的数学方法之前，先介绍一些定义。我们通过一些特定的非均质介质来说明一般的方法，即二维各向同性不相关和相关的对数正态介质（还简要提到了三维不相关的对数正态介质和方格图案的一些结果）。我们为小和大对数传导度方差进行了分析和数值研究。

收敛公式的明确推导在附录中找到。最后，我们总结了我们的主要结果，强调了它们的影响，并详细讨论了细化问题。