A Finite Analytic Method for Two-Dimensional Power-Law Fluid Flow in Heterogeneous Porous Media Using Unstructured Grids

地下深く、石油、ガス、水が蓄えられている岩石の複雑なネットワークの中では、流体は必ずしも単純で予測可能な直線を描いて動くわけではありません。岩石そのものが均一であることは稀であり、物質の性質が大きく異なるパッチワークのような構造をしています。ある層は開放された高速道路のようですが、別の層は密な壁のようです。エンジニアがこれらの流体の流れをモデル化しようとする際、岩石の特性が急激に変化することを考慮しなければなりません。さらに、流体自体も扱いにくい場合があります。水は容易に流れますが、製紙製造で使用されるスラリーや石油井で使用される掘削泥水など、多くの産業用流体は異なる挙動を示します。これらは「非ニュートン流体」であり、その粘性や流動への抵抗が、どれほど強く押されるかによって変化することを意味します。これらの複雑な流体が、性質が急変する岩石に遭遇すると、その動きを記述する数学は非常に困難になり、多くの場合、岩石の特性が変化するまさにその地点で破綻してしまいます。

数十年にわたり、科学者たちはこれらの流れを計算するための標準的な手法に頼ってきました。この手法は、本質的に2点間の岩石の特性を平均化するものです。このアプローチは、岩石が比較的均一であれば十分に機能しますが、岩石が高度に多様な場合には劇的に失敗します。そのような場合、標準的な手法は流体が実際にどれだけ移動できるかを一貫して過小評価し、膨大な誤差を生じさせます。これを解決するために、中国科学技術大学の研究者たちは、新しい計算手法を開発しました。この新手法は、単純な平均値に頼るのではなく、異なる2種類の岩石が接するエッジ（境界）において、流体がどのように振る舞うかという精密な数学的記述を使用します。これらの重要な境界における流体の真の、鋭い挙動を捉えることで、新手法は、岩石が極めて不均一で流体が複雑な場合であっても、地下流動のより明確で正確な姿を提供します。

研究者たちは、「べき乗則流体（power-law fluid）」として知られる特定のタイプの流体に焦点を当てました。これらは、流体を押し出す力と移動速度の関係が、特定の非線形な規則に従う流体です。現実の世界では、これは厚いポリマー溶液から特定の種類の原油に至るまで、幅広い物質をカバーしています。課題は、これらの流体が多孔質媒体（小さな穴の多い岩石）の中を流れる際に発生します。そこでは、浸透率（流れやすさ）が場所によって突然変化します。このような環境では、流体の圧力は滑らかに変化するのではなく、流れやすい岩石と流れにくい岩石の境界のすぐところで、鋭いスパイク、すなわち「特異点」を作り出します。従来のコンピュータモデルは、岩石を小さな三角形の格子（グリッド）に分割しますが、これらはこれらのスパイクに苦戦します。これらのモデルは、点の間で圧力が直線的に変化すると仮定していますが、これは境界においては物理的に不正確です。その結果、モデルは流速を低く見積もりすぎてしまい、まともな近似値を得るためだけに、膨大な数の極めて細かい格子点が必要になってしまいます。

これを修正するために、チームは「有限解析法（Finite Analytic Method）」を開発しました。核心となるアイデアは、境界における流体の振る舞いを推測するのをやめ、代わりにその振る舞いに関する既知の正確な数学的記述を使用することです。研究者たちは、岩石の性質が変わる点の近くでは、流体圧力が特定の曲線パターンに従って上昇または下降することに気づきました。彼らは、最も重要な流れの部分に対して、その曲線の正確な形状、すなわち「固有関数（eigenfunction）」を算出しました。この計算には、岩石の性質が異なる楔形の領域における流体の動きを記述する、複雑な一連の方程式を解くことが含まれました。これらの方程式の最小かつ最も重要な解を見つけることで、彼らは、問題の発生箇所における支配的な流体の振る舞いを捉えることができました。

この局所的な挙動が理解されると、研究者たちはその周囲に新しい数値スキームを構築しました。格子点同士を接続するために古い欠陥のある平均化手法を用いる代わりに、彼らは新しい解析解を使用して、グリッド内の各三角形セルの境界を横切る流れを計算しました。これにより、コンピュータは、高浸透率ゾーンと低浸透率ゾーンの間の狭い首の部分を流体がどのように通り抜けるのか、その挙動を正確に「知る」ことができるようになりました。その結果、近似値で問題を覆い隠すのではなく、最も困難な点における物理学を尊重したシミュレーションが可能となりました。

チームは、この新手法がどの程度うまく機能するかを確認するために、いくつかの困難なシナリオに対してテストを行いました。まず、高浸透率の岩石と低浸透率の岩石が交互に配置された、古典的な「チェッカーボード」パターンから開始しました。このセットアップでは、岩石の特性が100倍変化します。彼らが新手法を従来の手法と比較したとき、その差は歴然でした。従来の手法は、粗い格子において、真の値から時には90%近くも外れ、流動を大幅に過小評価していました。対照的に、新手法は、格子が細かくなくても、驚くほど正解に近い値を維持しました。格子を細かくしていくにつれて、新手法は非常に速く正解へと収束しましたが、従来の手法は極めて細かい格子を用いても追いつくのに苦労しました。

彼らはまた、規則的な六角形の領域、および実際の地下貯留層の混沌とした性質を模した、ランダムな対数正規分布を持つ浸透率を持つ領域においても、この手法をテストしました。あらゆるケースにおいて、新手法は優れていることが証明されました。それは岩石特性の極端な変動を容易に処理し、シミュレートされるべき乗則流体の具体的な種類に関わらず、正確性を維持しました。最も顕著な発見の一つは、新手法の精度が、流体や岩石の特定の特性に大きく依存しないことでした。流体が非常に粘り気が強いか薄いか、あるいは岩石のコントラストが緩やかか極端かに関わらず、手法は一貫した結果をもたらしました。研究者たちは、高い精度を得るためには、中程度の格子細分化で十分であることを発見しました。最悪のシナリオにおいても、誤差は8%未満に抑えられました。

この研究は、複雑な地下環境における流体の流れをシミュレートする方法における、大きな前進を意味します。欠陥のある平均化技術を、境界における真の物理学を尊重する解に置き換えることで、研究者たちは、より正確で効率的なツールを作り上げました。エネルギー生産から環境管理に至るまで、地球を通じて流体がどのように動くかを理解する必要がある業界にとって、これはシミュレーションがより信頼でき、より高速に実行できることを意味します。この手法は、まともな答えを得るために以前必要だった膨大な計算能力を必要とせず、複雑で不均質な系の高精度なモデリングをより身近なものにします。これらの知見は、世界が変化する特異な点を注意深く観察することで、全体に対するより優れた理解を築くことができることを示唆しています。