油藏数值模拟可以用来研究复杂的地下变化过程。经过长期的开采生产，储层受地下流体的冲刷或者流体饱和度的变化可能会对储层物性造成不同程度的影响。同时一些机理现象会对开发生产有影响，例如储层水锁。在没有考虑水锁效应时，模拟的产量要比实际的产量高，产量也不会发生突降。这个影响对于气藏的模拟更明显。
常规的数值模拟技术有一些功能可以模拟部分物性的变化，例如在不同的时间步改变储层网格的传导率，以此来近似物性发生的变化。但是复杂的情况很难考虑，一些数模工程师也尝试过用重启运算的方法来修改模型的物性，例如修改相渗曲线。这个过程比较繁琐，需要先在某一个时间步建立重启运算模型，然后修改重启运算模型中的相渗曲线，以此来模拟此时间步的岩性变化。如果需要多次修改参数来模拟，则需要多次地创建重启运算模型，修改和查看结果都不方便。那么在修改传导率和使用重启运算之外还有没有方法可以比较快捷的修改相渗曲线呢？这在tNavigator里可以较容易的实现。
在数值模拟中有很多可以表示不同物性的分区属性，例如相渗曲线分区、PVT分区等等。不同的分区编号代表了这个网格所使用的属性表的编号。如果我们可以在不同的时间步控制这些编号发生变化，就可以在不同时间步使用不同的物性参数，这里以相渗曲线分区为例，来模拟产生水锁时对产量的影响。
模型为油水两相黑油模型，数据体中有两组相渗曲线。初始状态时模型只使用了表1，随着生产的进行一些网格的物性会发生变化。两组油水相渗曲线如下图所示。表2为发生水锁时的油水相渗曲线。
模型的初始含水饱和度为0.22，当含水饱和度高于0.35时或者压力低于300Bar时，修改网格相渗曲线分区编号，使用油相渗透率更低的曲线。初始的相渗分区编号如下图所示（IK剖面）。
物性刚发生变化时（生产约8年后）的相渗分区编号如下所示。红色部分的网格开始使用第2组相渗曲线。
模型的产油量曲线和含水率曲线如下图所示：
图中产量降低和含水上升的时间点就是物性发生了变化。以上的过程只需要在数据体的Schedule部分添加如下关键字即可实现：
ARITHMETIC
SATNUM = IF( SWAT > 0.35 | PRESSURE < 300, 2, 1 ) /
除了相渗曲线分区的修改外，还可以修改其他的属性，例如前面提到的MULTX传导率参数。可以这样修改。
ARITHMETIC
MULTX = MULTX * SWAT + 0.2 * PRESSURE
/
这里是举了个例子来说明属性的计算方法非常灵活，tNavigator支持的算数方法有很多种。除了常规的加减乘除外，还有很多内置的公式如LOG10、SUM、MIN、SQRT等，经过组合可以实现复杂的计算。几乎所有的模型属性都可以用来组合公式计算新的物性参数（在等式右侧使用），可实现对生产过程中储层物性变化的模拟。本例简单地考虑了含水升高和压力降低孔隙变化的双重影响，更复杂的机理可以有更多的判断条件和不同条件下的相渗曲线组合来模拟。
| 
发表于 2016-6-3 17:02:27