Physically-Consistent Fracture Capillary Pressure Type Curves and Relative Permeability: A Unified Aperture-Based Model without Empirical Matching

技術要約：物理的一貫性を備えた割れ目毛細管圧タイプカーブおよび相対浸透率

問題の定義
割れ目（フラクチャー）ネットワークの複雑な幾何学的形状と、マトリクス・割れ目間の複雑な相互作用により、割れ目系における流体流動の正確なモデリングは依然として大きな課題となっている。割れ目の毛細管圧（$P_c$）および相対浸透率（$k_r$）は、相間の干渉と置換効率を支配する基本的な構成関係であるが、従来の石油工学の実務では、過度に単純化された仮定に依存することが多い。具体的には、シミュレータはしばしば、割れ目における毛細管圧をゼロとし、相対浸透率を線形（Xカーブ）の挙動と仮定する。近年の実験およびモデリングの証拠は、これらの仮定に健全な物理的根拠が欠けており、性能予測において重大な誤差を招く可能性があることを示唆している。さらに、既存のモデルは、物理的な意味が不明確であったり、ケースごとに異なる複数の「フィッティングパラメータ」に依存したりすることが多く、異なる割れ目形状に対する予測能力を制限している。割れ目の微視的構造特性（特に開口幅の分布）と、これらの流動パラメータとの間の直接的な関連性に決定的な空白が存在する。

手法
本研究は、割れ目の開口幅分布と基礎的な流体力学を統合することにより、割れ目の$P_c$および$k_r$モデルのための統一された物理ベースのフレームワークを開発する。コアとなる手法は、以下のステップを含む：

1. 統計的特性評価： 割れ目システムは、個々の割れ目が空間的に変化する開口幅を持つ、完全に連結されたネットワークを含む単位体積として概念化される。本研究では、割れ目の開口幅、幅、および長さを特徴付けるために対数正規分布を用い、割れ目の配向には正規分布を用いる。開口幅の標準偏差（$\sigma_b$）が主要な支配パラメータとして特定される。
2. 相飽和度の導出： 毛細管圧の閾値によって定義される特定の開口幅範囲に基づき、相体積を計算することによって、濡れ性相および非濡れ性相の移動性および不動性飽和度を導出する。これらの範囲は、残留濡れ性相、残留非濡れ性相、および全相飽和度に対応する。
3. タイプカーブの開発：
   • 毛細管圧： 無次元の濡れ性相飽和度は、有効飽和度および無次元毛細管圧（$P_{cDf}$）の関数として定式化される。対数正規開口幅分布を組み込むことで、誤差関数（erf）を用いたドレナージ（排水）およびインビジション（吸水）の割れ目毛細管圧タイプカーブが導出される。これらの曲線は、単一の物理パラメータである$\sigma_b$と、無次元境界圧力（$P_{cDminf}$, $P_{cDmaxf}$）に依存する。
   • 相対浸透率： 本モデルは、導管内の流体流動に関するハーゲン・ポアズイユ方程式と、多孔質媒体におけるダルシーの法則を統合することで、割れ目の相対浸透率を導出する。この積分は、ドレナージおよびインビジションのプロセスにおける開口幅分布と相飽和度を考慮している。曲がりやすさ（Tortuosity）係数が組み込まれ、その指数は実験データとの検証を通じて決定される。
4. 検証およびマッチング： モデルは、幅広い理論的枠組み（例：Pruess and Tsang, Zhou and Wheater, Liu et al.）および実験データセット（例：Reitsma and Kueper, Huo et al., Chen, Diomampo）に対して検証される。
   • $\sigma_b$が既知の場合、モデルは直接比較可能な曲線を作成する。
   • $\sigma_b$が未知の場合、タイプカーブ・マッチング手順が用いられる。実験的な$P_c$対飽和度のデータを対数-対数スケール上にプロットし、開発された無次元タイプカーブにマッチさせることで、最適な$\sigma_b$およびその他のパラメータを抽出する。
5. 定量的指標： モデルの精度は、平均二乗誤差（MSE）および平均絶対相対誤差（AARE）を用いて評価される。

主な貢献

• 統一された開口幅ベースのモデル： 本研究は、任意の経験的なマッチング係数を排除し、割れ目の開口幅分布から割れ目の毛細管圧と相対浸透率の両方を直接導出する単一のフレームワークを提示する。
• 単一の物理パラメータ： 本モデルは、割れ目の開口幅の標準偏差（$\sigma_b$）を唯一の入力パラメータとして利用する。このパラメータは物理的に解釈可能であり、幾何学的な不均質性を表しており、複数のケース固有のフィッティングパラメータを必要とする従来のモデルとは対照的である。
• タイプカーブの手法： 割れ目の毛細管圧タイプカーブの導入により、以下の二重の適用が可能となる：(1)既知の開口幅分布から$P_c$および$k_r$を予測すること、および(2)実験的な$P_c$データからタイプカーブ・マッチングを通じて割れ目の開口幅分布を特性評価すること。
• 数学的な簡潔さと堅牢性： 導出された方程式は標準的な統計関数（誤差関数）を利用しており、商用の油層シミュレータへのシームレスな統合が可能となるよう設計されている。

結果
開発されたモデルは、多様な割れ目のタイプおよび流動条件下で、極めて高い堅牢性と有効性を示した：

• 毛細管圧： モデルは、理論モデル（例：Zhou and Wheater）をほぼ完璧に一致させて再現し（MSEは$10^{-9}$から$10^{-10}$の範囲）、実験データ（例：Huo et al., Reitsma and Kueper）とも低いAARE値（通常20%未満）で一致した。特筆すべきは、有効応力の増加に伴い最適化された$\sigma_b$が増加するという物理的傾向を捉えており、これは割れ目の閉塞と不均質性の増大を反映している。
• 相対浸透率： モデルは、ドレナージおよびインビジションの両方のプロセスについて、相対浸透率曲線を正確に予測した。また、理論的限界（例：$\lambda \to \infty$のBrooks and Corey）および実験データセット（例：Chen, Diomampo）との強い一致を示した。本モデルは、相間の干渉を二次的な実験変数から分離することに成功しており、従来の飽和度ベースのプロットよりも厳密な検証を提供した。
• パラメータの効率性： すべての検証ケースにおいて、本モデルは単一のパラメータ$\sigma_b$のみを用いながら、複雑な既存モデルと同等またはそれ以上の精度を達成した。

意義および主張
本論文は、提案されたフレームワークが、割れ目の毛細管圧および相対浸透率を決定するための、コストのかかる時間のかかる実験に代わる、信頼性が高く費用対効果の高い選択肢を提供すると主張している。モデルを経験的な曲線フィッティングではなく、割れ目の開口幅分布という物理的実態に根ざすことで、本研究は石油エンジニアにとって物理的一貫性のあるツールを提供している。

著者らは、モデルの数学的構造が商用の油層シミュレータへの直接実装を容易にし、割れ目系のより現実的な大規模シミュレーションを可能にすることを強調している。本研究は、割れ目の毛細管圧はゼロでも一定でもなく、割れ目の相対浸透率は線形の「Xカーブ」から大きく逸脱することを結論付けており、これらの現象は、開口幅ベースのモデルを用いて正確に捉えられなければならない。本研究は、複雑な解析解と実用的な石油工学アプリケーションとの間の溝を埋め、任意のフィッティングパラメータの曖昧さなしに、割れ目媒体を特性評価するための標準化された手法を提供している。