Immiscible two-phase flow in porous media: a statistical mechanics approach

この論文は、シャノン情報エントロピーに基づく統計力学の枠組みを用いることで、多孔質媒体内の不混和二相流を、細孔スケールの物理と整合性を持ちながら管理可能な複雑さで巨視的に記述する新たなアプローチを概説し、従来の相対浸透率理論の限界を克服する可能性を示しています。

要約

1. 従来の考え方：「スポンジの迷路」の難しさ

まず、スポンジ（多孔質媒体）の中を水と油が流れる様子を想像してください。
スポンジには無数の小さな穴（ポア）があります。水と油は混ざらないので、互いに押し合い、穴の壁に引っかかったり、細い道を行き交ったりします。

• これまでの常識（相対浸透率理論）：
  これまでの科学者は、「スポンジ全体を一つの大きな管」として扱い、水と油の「混ざり具合（飽和度）」だけで、どれくらい流れやすいかを表す数式を作ってきました。
  しかし、これは**「スポンジの内部の複雑な迷路を無視して、外側から見た結果だけを当てはめる」**ようなもので、厳密には不自然な部分がありました。特に、水と油がぶつかり合う「界面」の動きを無視しすぎていたのです。

2. 新しいアプローチ：「統計力学」という新しい眼鏡

この論文の著者たちは、**「統計力学（Statistical Mechanics）」**という、通常は「熱」や「分子の動き」を扱う物理学の手法を、この「スポンジ内の流れ」に応用しました。

• ジャインズのアイデア：
  1950 年代、物理学者のジャインズは、「熱がなくても、情報の欠如（わからないこと）を数値化して、統計的な法則を見つけられる」と提案しました。
  彼らはこれをスポンジに応用しました。「スポンジの穴一つ一つで何が起きているか（ミクロ）」をすべて追うのは不可能ですが、**「全体としてどう振る舞っているか（マクロ）」**を、情報の統計から導き出そうというのです。

3. 発見された「新しい温度」と「新しい圧力」

この新しい視点で見ると、スポンジの中の流れには、驚くべき「熱力学のような法則」が働いていることがわかりました。

• 「アジチュア（Agiture）」＝ 揺さぶりの温度
  通常の熱力学には「温度（分子がどれくらい激しく揺れているか）」がありますが、ここでは**「アジチュア」**という新しい概念が登場しました。

  • 例え： スポンジに水を流すとき、圧力を強くすればするほど、水と油の境界線（界面）が激しく揺さぶられます。この「揺さぶりの強さ」が、この世界における「温度」の役割を果たしているのです。
  • 圧力勾配（流れを押し出す力）が大きいほど、この「アジチュア」は高くなり、流体は活発に動き回ります。

• 「共移動速度（Co-moving velocity）」＝ 二人組のダンス
  これが最も重要な発見です。水と油は、それぞれ独立して流れているように見えますが、実は**「互いに影響し合い、一緒に動く」**部分があることがわかりました。

  • 例え： 混雑した駅のホームで、水と油が流れていると想像してください。水が油を押し、油が水を引っ張る。彼らは単に並走しているだけでなく、**「互いの足並みを合わせて、少しずれて動く」**ような関係にあります。
  • この「ずれた動き」を表すのが**「共移動速度」**です。従来の理論では見逃されていたこの要素を考慮することで、実験結果と理論が完璧に一致するようになりました。

4. 「ガラス状態」という不思議な現象

さらに面白いことに、流れの速さや混ざり具合によって、スポンジ内の液体の状態が**「ガラス状態」と「普通の液体状態」**に分かれることがわかりました。

• ガラス状態（Phase Ib）：
  流れがゆっくりな時、水と油の境界線は「凍りついた」ように動かないように見えますが、実は微細に揺れています。これは、**「ガラスが固まっているが、分子レベルではまだ動いている」**状態に似ています。
  • この状態では、過去の履歴（どちらを先に流したか）によって結果が変わる「ヒステリシス」という現象が起きやすくなります。

5. まとめ：なぜこれが重要なのか？

この論文は、**「スポンジの中の複雑な流れも、実は『熱力学』のような美しい法則で説明できる」**と示しました。

• 従来の考え方： 「スポンジは汚くて複雑だから、エンジニアが経験則で適当に計算すればいい」というもの。
• 新しい考え方： 「スポンジの中にも、温度や圧力のような『新しい物理量』が存在し、それらを使って正確に予測できる」というもの。

最終的なメッセージ：
スポンジの中の水と油の動きは、単なる「泥臭い現象」ではなく、**「情報の統計と、新しい温度（アジチュア）が織りなす、高度なダンス」**だったのです。この新しい理論を使えば、石油の採掘や地下水の管理、さらには新しい材料の開発など、実社会の問題をより正確に解決できるようになるでしょう。

---

一言で言うと：
「スポンジの中の水と油の動きを、**『揺さぶり（アジチュア）』という新しい温度と、『二人組のダンス（共移動速度）』**という概念を使って、物理学の最高峰である『熱力学』と同じレベルで美しく説明することに成功した！」という論文です。

技術概要

論文要約：多孔質媒体中の不混和二相流に対する統計力学アプローチ

1. 背景と問題提起

多孔質媒体における不混和二相流（例：油と水）の物理において、中心的な課題は「スケールアップ問題」、すなわち、ポアスケール（微視的）の物理を、媒体を連続体とみなせるマクロなスケール（ダルシースケール）にどのように記述するかという点にあります。

従来のアプローチは、相対浸透率理論（Relative Permeability Theory）に依存しています。これは、1930 年代に確立された現象論的な方程式群であり、浸透圧（キャピラリ圧力）や相対浸透率が飽和度のみの関数であるという強い仮定に基づいています。しかし、この理論には以下の弱点があります。

• 現象論的であり、ポアスケールの物理と直接的な結びつきが弱い。
• 複雑な挙動（ヒステリシスや非線形性）を扱う際に限界がある。
• これを超えようとする試みは、計算複雑性の爆発により実用的な応用に至っていない。

本論文は、この問題に対し、エドウィン・T・ジェインズ（Edwin T. Jaynes）が提唱した情報理論的エントロピーに基づく統計力学の一般化を適用し、連続体スケールでの熱力学的な形式を構築することを目的としています。

2. 手法：非熱的統計力学の適用

著者らは、ポアスケールの流体力学をダルシースケールにスケールアップするために、ジェインズの最大エントロピー原理を流体力学問題に適用しました。

• 代表要素面積（REA）: 流れ方向に垂直な断面内に、統計的代表性を持つ「代表要素面積（REA）」を定義します。
• 変数の定義:
  • 湿潤相と非湿潤相の体積流量（$Q_w, Q_n$）およびそれらが占める断面積（$A_w, A_n$）。
  • 流体配置 $X$ に対する確率分布 $p(X)$。
• エントロピーの最大化: シュノン情報エントロピー $S = -\int p(X) \ln p(X) dX$ を、流量や面積の平均値に関する制約条件の下で最大化します。
• 分配関数の導出: これにより、確率分布 $p(X)$ が指数関数の形（ギブス分布に類似）で得られ、新しいラグランジュ乗数（$\lambda_u, \lambda_w, \lambda_p$）が導入されます。

このプロセスにより、ポアスケールの微視的記述から、熱力学に類似したマクロな形式が自然に導き出されます。

3. 主要な成果と発見

A. ダルシースケールにおける「非熱的熱力学」の構築

統計力学の枠組みから、以下の新しい変数と関係式が導かれました。これらは通常の熱力学の構造を模倣していますが、温度や化学ポテンシャルとは異なる物理的意味を持ちます。

1. アジチュア（Agiture, $\theta$）: 温度に相当する変数。圧力勾配の大きさ $|\nabla p|$ に比例すると解釈されます（$\theta \propto |\nabla p|$）。
2. 流動導数（Flow Derivative, $\mu$）: 化学ポテンシャルに相当する変数。湿潤相の断面積（または飽和度）に対する共役変数です。
3. 流動圧（Flow Pressure, $\pi$）: ポア面積に対する共役変数。通常の熱力学には対応する概念がありません。

これらを用いることで、平均流速 $v_p$ と圧力勾配、飽和度の間の構成方程式が、熱力学的な状態方程式として記述可能になります。

B. 共移動速度（Co-moving Velocity, $v_m$）の発見と重要性

本アプローチの最も重要な成果の一つは、共移動速度 $v_m$ の導入と、それが流体力学における中心的な役割を果たすことを示したことです。

• 定義: 各相の実際の浸透速度（$v_w, v_n$）と、熱力学的に定義された速度（$\hat{v}_w, \hat{v}_n$）の間には、以下の関係があります。
  $$v_w = \hat{v}_w - S_n v_m$$
  $$v_n = \hat{v}_n + S_w v_m$$
• 意味: $v_m$ は、二相が互いに「共動」する速度成分を表します。これは、従来の相対浸透率理論では無視され、あるいは暗黙的にゼロと仮定されてきた項です。
• 実験的検証: 相対浸透率データ、動的ポアネットワークモデル、および格子ボルツマンシミュレーションの分析により、$v_m$ は流動導数 $\mu$ に対して線形（または単純な関数形）で記述できることが示されました（図 3 参照）。
• 熱力学との類似性: 著者らは、この $v_m$ が、通常の熱力学における流体混合物の「共モル体積（co-molar volume, $\nu_m$）」と数学的に同型であることを示しました（図 4 参照）。これは、不混和二相流の理論が通常の熱力学の概念を拡張・補完する可能性を示唆しています。

C. 相図とガラス転移

ダルシースケールでの定常流の相図（図 1）が提示されました。

• ガラス状相（Glassy Phase）: 低キャピラリ数領域において、界面の運動が断続的（intermittent）であり、 Edwards-Anderson 秩序パラメータがゼロでない状態。この領域では、流速と圧力勾配の間にべき乗則（$v \propto |\nabla p|^\beta$）が観測されます。
• 非ガラス状相（Non-glassy Phase）: 高キャピラリ数領域。界面がすべて動き、流速と圧力勾配の間に線形関係（ダルシーの法則）が回復します。
• 臨界線: これらの相を分ける境界は、通常の熱力学的相転移に相当する「臨界線」として扱え、ガラス転移の特性（ヒステリシスや大きな揺らぎ）を説明できます。

4. 結論と意義

この論文は、多孔質媒体中の不混和二相流を記述する新たなパラダイムを提供しました。

1. 理論的統合: ポアスケールの複雑な流体力学を、統計力学の枠組みを通じて、熱力学的な形式（状態方程式、エントロピー、共役変数）へと統一的に記述することに成功しました。
2. 実用的な変数の発見: 「共移動速度（$v_m$）」という新しい変数が、従来の相対浸透率理論の欠陥を補い、実験データと理論を結びつける鍵であることを示しました。
3. 非平衡・非熱的現象の理解: 温度という概念が存在しない流体力学問題において、熱力学の数学的構造が有効に機能することを証明し、粉末や粒状物質などの他の非熱的システムへの応用可能性も示唆しています。
4. 将来展望: 定常流の枠組みを超え、非平衡熱力学を適用することで、キャピラリ圧力関数（$p_c$）のような問題のある概念を不要にし、より一般的な非定常流の記述が可能になると期待されています。

総じて、この研究は「多孔質媒体の流体力学は単なる工学的な問題ではなく、深い物理的・統計力学的構造を持つ分野である」ことを示し、基礎物理学と応用科学の架け橋となる重要な貢献です。