Residual Saturation under Pressure-Controlled Drainage

技術要約：圧力制御排水下における残留飽和

問題提起
本論文は、多孔質媒体の物理学における根本的な未解決問題、すなわち無限サイズ極限における圧力制御排水下の残留飽和の挙動を取り上げている。標準的な浸透（Invasion Percolation: IP）モデルは、システムサイズが無限大に近づくにつれて侵入相の飽和が消失すること（突破クラスタのフラクタル構造によって有限の防御相が捕捉されるのみであること）を予測するが、キャピラリー閾値に対してインレットに接続されたすべての孔隙へのアクセスを許容する「圧力制御プロセス」において、熱力学的極限における防御相の有限な残留飽和が存在するかどうかは依然として不明である。この区別は、水攻めなどのフィールドスケールの現象を理解する上で極めて重要であり、侵入領域が貯留層の寸法に近づくにつれて残留油飽和が有限に留まるのかという問いを提起している。

手法
著者らは、圧力制御排水を、細孔ネットワークグラフ上での**トラッピングを伴うボンドパーコレーション（Bond Percolation with Trapping: BPT）**として定式化している。

• グラフ表現: 多孔質媒体は、ノードが孔隙領域を表し、エッジが割り当てられた入口半径（キャピラリー閾値）を持つスロートを表す無向グラフとしてモデル化されている。
• プロセスの定義: 標準的なIPとは異なり、BPTはグローバルな圧力制約の下で動作する。キャピラリー圧力が上昇するにつれ、「アクティブグラフ」（現在の圧力未満の入口閾値を持つスロートで構成されるグラフ）が拡大する。
• トラッピング機構: 浸入は、アクティブグラフ上でのインレットへの連結性によって定義される。トラッピングは、防御相がアウトレット集合への連結性を失ったときに発生する。アルゴリズムは以下の手順でシステムを更新する：(1) 圧力を上昇させて新しいスロートを活性化させる、(2) 侵入領域（インレットに接続されたアクティブな孔隙）を特定する、(3) 防御領域（アウトレットに接続された非侵入孔隙）を特定する、(4) 残りの非接続な孔隙をすべて「捕捉（trapped）」として分類する。
• シミュレーション: 本研究では、様々な2D格子（正方形、三角形、ボロノイ）および3D格子（ダイヤモンド、および配位数変化を伴う単純立方格子）を用いて、BPTを標準的な浸透（IP）と比較している。スロゲート半径には対数正規分布を用い、アンサンブル平均を用いて有限サイズ効果を分析している。

主な貢献と結果

1. 無限極限における有限の残留飽和:
   主要な知見は、圧力制御排水が空間充填型の侵入状態をもたらし、システムサイズ $L \to \infty$ において防御相の残留飽和が非ゼロの有限な極限値 ($0 < S_\infty < 1$) に収束することである。これは、IPにおいて突破時に侵入飽和が消失する（$S_{IP} \to 0$）こと（侵入クラスタがフラクタルであり、防御相の広大な領域をバイパスするため）と鮮明な対照をなしている。

2. 二次元における有限サイズスケーリング:
   2Dにおいて、残留飽和の熱力学的極限からの偏差は、以下の有限サイズスケーリング則に従う：
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   ここで、指数 $\delta \approx 0.25$ である。この指数は異なる格子タイプ（正方形、三角形、ボロノイ）間で普遍的であり、微視的な詳細には依存しない。この緩やかな収束（$\Delta S_r \sim L^{-0.25}$）は、2Dにおいて代表要素体積（REV）を達成するために、極めて大きなシステムサイズが必要であることを示唆している（例：システムサイズを倍にしても、有限サイズによるバイアスは $\approx 1/16$ の係数でしか減少しない）。

3. 三次元における次元性と配位数の影響:
   3Dでは、有限サイズ補正は指数 $\delta \approx 0.75$ で著しく速く減衰する。その結果、無限サイズ極限を近似するために必要な代表体積は、2Dよりも3Dの方がはるかに小さい。さらに、漸近的な残留飽和は3Dにおいても有限であるが、その値は格子の配位数に依存する。結合度が高いほど代替となる浸入経路が増え、捕捉される割合は減少するが、それが完全に消失することはない。

4. 浸透（Invasion Percolation）との区別:
   本論文は、BPTとIPの違いが単なる曲線の形状やバーストのダイナミクスの問題ではなく、漸近的な状態における根本的な違いであることを明らかにしている。IPは突破時にフラクタルな侵入クラスタを生じさせるが、BPTは有限の捕捉分を共存させたまま、空間充填型の侵入領域を生じさせる。

意義
本論文は、連結性とトラッピングが圧力制御排水下の残留飽和を決定する要因であることを示すことにより、パーコレーション理論と圧力・飽和関係の間の直接的な理論的関連性を確立している。有限の容積が無限システム極限においても捕捉されたまま残ることを示すことで、本研究は突破状態を超えて標準的な浸透（IP）の概念を拡張している。これは、代表要素体積の問題が次元に非常に敏感であり、2Dの方が3Dよりもはるかに制約が厳しいことを示唆している。これらの結果は、圧力制御シナリオにおいては、システムサイズに関わらず有限の割合の防御相が回収不能であることを意味しており、地下流体の再分配および回収プロセスに対する直接的な示唆を与えている。