19、单相流体流动方程相关知识解析

单相流体流动方程相关知识解析

1. 网格块压力计算与物质平衡检查

在网格块 5 的计算中，通过一系列运算得到方程 (41.6388914p_{n + 1}^4 - 46.6943060p_{n + 1}^5 = -19290.5407)。求解方程((7.150))、((7.151))、((7.152))、((7.153))和((7.154))后，得到未知压力值：(p_{n+1}^1 = 3252.93) psia，(p_{n+1}^2 = 3300.96) psia，(p_{n+1}^3 = 3348.99) psia，(p_{n+1}^4 = 3396.89) psia，(p_{n+1}^5 = 3442.25) psia。

使用方程((7.85))对第二个时间步进行物质平衡检查：
[
IMB = \frac{\sum_{n = 1}^{5}V_{bn}\left[\frac{\phi^{\circ} n}{c + c {\phi}}\right]\frac{\alpha_cB^{\circ}\Delta t}{p_{n + 1}^n - p_n^n}}{\sum_{n = 1}^{5}q_{n + 1}^{scn}+\sum_{n = 1}^{5}\sum_{l\in\xi_n}q_{n + 1}^{scl,n}}
]
计算结果为(1.000006)，表明物质平衡得到满足。

2. 可压缩流体流动方程

可压缩流体在储层温度下的密度、体积系数（FVF）和粘度是压力的函数。在流动方程((7.12))的左侧，FVF、粘度和密度可假设为常数，但在每个时间步开始时至少更新一次。为了保持物质平衡，积累项用