Residual Saturation under Pressure-Controlled Drainage

技术摘要：压力控制排水下的残余饱和度

问题陈述
本文探讨了多孔介质物理学中一个基本且尚未解决的问题：在无限尺寸极限下，压力控制排水过程中的残余饱和度行为。虽然标准的入侵渗透（Invasion Percolation, IP）模型预测，随着系统尺寸趋于无穷大，被入侵相的饱和度会趋于零（仅因突破簇的分形特性而留下有限比例的被防御流体），但目前尚不清楚压力控制过程——即允许入侵相访问所有在给定毛细管阈值下与入口相连的孔隙——是否会在热力学极限下导致被防御相具有有限的残余饱和度。这一区别对于理解诸如水驱等场尺度现象至关重要，因为这类现象涉及：当入侵区域接近储层维度时，残余油饱和度是否保持有限。

研究方法
作者将压力控制排水建模为孔隙网络图上的带捕获键渗透（Bond Percolation with Trapping, BPT）。

• 图表示法： 多孔介质被建模为一个无向图，其中节点代表孔隙区域，边代表具有分配半径（毛细管阈值）的喉道。
• 过程定义： 与通过局部最易接近路径逐个喉道推进的标准 IP 不同，BPT 在全局压力约束下运行。随着毛细管压力的增加，“活性图”（由阈值低于当前压力的喉道组成）不断扩张。
• 捕获机制： 入侵被定义为在活性图上与入口的连通性。当部分被防御相失去与出口集的连通性时，便发生了捕获。算法通过以下步骤更新系统：(1) 增加压力以激活新喉道；(2) 识别入侵区域（与入口连通的活性孔隙）；(3) 识别被防御区域（与出口连通的非入侵孔隙）；以及 (4) 将所有剩余的非连通孔隙归类为“被捕获”状态。
• 模拟： 本研究在各种 2D 格点（正方形、三角形、Voronoi）和 3D 格点（钻石型、具有不同配位数 的简单立方格点）上将 BPT 与标准入侵渗透（IP）进行了对比。研究利用对数正态分布的喉道半径和系综平均来分析有限尺寸效应。

主要贡献与结果

1. 无限极限下的有限残余饱和度：
   核心发现是，压力控制排水会导致一个充满空间的入侵态，即随着系统尺寸 $L \to \infty$，被防御相的残余饱和度趋于一个非零的有限极限 ($0 < S_\infty < 1$)。这与标准 IP 形成了鲜明对比，在 IP 中，由于入侵簇是分形的并绕过了大量被防御流体，入侵饱和度在突破时趋于零 ($S_{IP} \to 0$)。

2. 二维空间中的有限尺寸定标：
   在 2D 中，残余饱和度偏离其热力学极限的规律遵循有限尺寸定标律：
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   其中指数 $\delta \approx 0.25$。该指数在不同类型的格点（正方形、三角形、Voronoy）上是普适的，且独立于微观细节。这种缓慢的收敛性（$\Delta S_r \sim L^{-0.25}$）意味着，在 2D 中要获得代表性单元体积需要极其庞大的系统尺寸（例如，将系统尺寸增加一倍，其有限尺寸偏差仅减少约 $1/16$）。

3. 三维中的维度与配位数效应：
   在 3D 中，有限尺寸修正衰减得显著更快，其指数 $\delta \approx 0.75$。因此，在 3D 中，近似无限尺寸极限所需的代表性体积比 2D 小得多。此外，虽然在 3D 中渐近残余饱和度仍为有限值，但其数值取决于格点的配位数；较高的连通性提供了更多的替代入侵路径，从而降低了被捕获部分的比例，但并未将其完全消除。

4. 与入侵渗透的区别：
   本文阐明，BPT 与 IP 的区别不仅在于曲线形状或爆发动力学的差异，更是渐近状态的本质区别。IP 在突破时产生分形入侵簇，而 BPT 则产生一个充满空间的入侵区域，并在无限极限下与有限比例的被捕获部分共存。

意义
本文通过证明在压力控制排水下，连通性与捕获是决定残余饱和度的关键因素，建立了渗透理论与压力-饱和度关系之间的直接理论联系。通过展示即使在无限系统极限下仍存在有限比例的被捕获孔隙空间，这项工作将标准的入侵渗透图景扩展到了突破状态之后。这为解释宏观本构方程提供了严谨的框架，并表明“代表性单元体积”问题对维度高度敏感，在 2D 中的限制远比 3D 严苛。结果表明，在压力控制场景下，无论系统尺寸如何，都会有有限比例的被防御流体无法回收，这一发现对地下流体重新分布及采收过程具有直接的影响。