[論文レビュー] Upscaling of unsaturated flow in fractured porous media

Florian List, Kundan Kumar+2 more|arXiv (Cornell University)|Jul 16, 2018

Advanced Mathematical Modeling in Engineering参考文献 48被引用数 5

ひとこと要約

この論文は、薄い亀裂（幅 ε → 0）の極限をとることで、破壊的多孔質媒体における不飽和流れの有効モデルを厳密に導出している。導出されたアップスケールモデルは、ε に伴う密度比と透過率比のスケーリングに依存する。解析により、指数 κ と λ に応じて異なる有効方程式（例えば、亀裂界面での1次元リチャーズ方程式）が得られ、ε ≤ 10⁻³ の場合に数値的に検証され、誤差は 10⁻⁵ 未満であった。

ABSTRACT

In this work, we consider a mathematical model for flow in a unsaturated porous medium containing a fracture. In all subdomains (the fracture and the adjacent matrix blocks) the flow is governed by Richards' equation. The submodels are coupled by physical transmission conditions expressing the continuity of the normal fluxes and of the pressures. We start by analyzing the case of a fracture having a fixed width-length ratio, called $\varepsilon > 0$. Then we take the limit $\varepsilon o 0$ and give a rigorous proof for the convergence towards effective models. This is done in different regimes, depending on how the ratio of porosities and permeabilities in the fracture, respectively matrix scale with respect to $\varepsilon$, and leads to a variety of effective models. Numerical simulations confirm the theoretical upscaling results.

研究の動機と目的

• 薄い亀裂幾何における下位次元の亀裂モデルの使用を数学的に正当化し、厳密にそれらを導出すること。
• 密度比と透過率比が亀裂厚さ ε にどのようにスケーリングするかに応じて、異なる有効モデルがどのように出現するかを特定すること。
• 正の幅を持つ亀裂問題（ε-問題）が ε → 0 の極限で有効モデルに収束することを確立すること。
• 地下流れシミュレーションで用いられる離散的亀裂モデル手法の物理的に一貫した基礎を提供すること。
• 理論的結果を数値シミュレーションにより検証し、ε が小さい場合に高い精度を示すことを示すこと。

提案手法

• マトリックスブロックと幅 ε の薄い亀裂の両方でリチャーズ方程式を定式化し、界面で圧力と法線方向流量が連続であることを仮定する。
• スケーリング仮定を導入：密度比 ∝ ε^κ および透過率比 ∝ ε^λ で、κ > -1 および λ < 1 を満たす。
• 厳密な漸近解析とホメオジゼーション技術を用いて、ε → 0 の極限における有効モデルを導出する。
• κ と λ の値に応じて複数の有効モデルを導出し、λ < 1 かつ κ > -1 の場合に亀裂界面で1次元リチャーズ方程式が成立することを示す。
• 弱形式とエネルギー推定を用いて、ε-問題の弱解の存在と一意性を証明し、有効モデルへの収束を示す。
• van Genuchten–Mualem パrameters を用いた数値シミュレーションを行い、ε-解と有効モデルを比較し、部分領域ごとの L² 誤差を計算する。

実験結果

リサーチクエスチョン

• RQ1 亀裂の幅が小さいが正の値（ε > 0）である場合、密度比と透過率比のスケーリングに応じて、ε → 0 の極限でどのような有効モデルが得られるか？
• RQ2 密度比と透過率比のスケーリングにおける指数 κ と λ が定める異なる領域が、アップスケール方程式の構造にどのように影響するか？
• RQ3 リチャーズ方程式で記述される非退化な不飽和流れにおいて、ε-問題が有効モデルに厳密に収束することを証明できるか？
• RQ4 特に ε が小さい場合に、有効モデルが元の ε-問題を圧力ヘッドと水分含有率の観点でどれほど正確に近似できるか？
• RQ5 有効モデルの数値解が、特に流れのダイナミクスと水分含有率分布の観点で、元の2次元亀裂問題の本質的特徴を捉えているか？

主な発見

• 有効モデルは、指数 κ と λ に強く依存する。λ < 1 かつ κ > -1 の場合、亀裂は1次元界面として、縮小されたリチャーズ方程式に従う。
• ε ≤ 10⁻³ の場合、亀裂およびマトリックスブロックにおける L² 誤差は 10⁻⁵ 未満であり、非線形ソルバーの許容誤差と一致しており、高い数値的精度を示している。
• ε = 10⁻² に対しても L² 誤差は 10⁻⁴ 未満であり、ε → 0 の極限におけるアップスケールモデルの収束性が強く裏付けられている。
• ε-解における亀裂幅に沿った x 方向平均圧力ヘッドが、有効モデルの解とよく一致しており、平均化仮定の妥当性が検証された。
• 界面で圧力が連続であるにもかかわらず、低い水力的透過度のため、亀裂はマトリックスブロックよりも水分含有率が低く保たれる。これはモデルの予測通りである。
• 理論的枠組みは、実用的なシミュレーションにおける離散的亀裂モデルの使用を支持し、その物理的整合性に数学的根拠を提供している。