[論文レビュー] Modern theory of hydraulic fracture modeling with using explicit and implicit schemes

A. M. Linkov|arXiv (Cornell University)|May 8, 2019

Hydraulic Fracturing and Reservoir Analysis参考文献 51被引用数 6

ひとこと要約

本稿は、水力破砲（HF）シミュレーションにおける速度方程式の中心的役割に焦点を当て、近端領域におけるCFL制約が緩やかであるため、明示的時間積分法が陰解法を上回る可能性があることを示している。波動的挙動が近端領域に見られるため、普遍的漸近的被覆（universal asymptotic umbrella）の概念を導入し、上流型スキームの効率的適用が可能となる。これにより、特に形状変化を伴う3次元問題において、明示的スキームが競争力あるか、あるいは優位な性能を示すことが明らかになった。

ABSTRACT

The paper presents novel results, obtained on the basis of the modified theory of hydraulic fractures (HF). The theory underlines significance of the speed equation. When applied to numerical simulation of HF, the theory reveals three distinct issues: (i) modeling the central part of a HF; (ii) modeling the near-front zone; and (iii) tracing changes in the shape of a fracture contour. Modeling the central part leads to a stiff system of ODE in time, what strongly complicates its integration. For explicit schemes, it requires small time steps to meet the CFL condition. For implicit schemes, it requires proper preconditioners. The gains and flaws of the two strategies are discussed. It is noted that a rough spatial mesh may be used in the central part. Modeling the near-front zone reveals the vital role of the speed equation for HF modeling by any method. Its asymptotic analysis has resulted in the fundamental concept of the universal asymptotic umbrella. For the near-front zone, it is also established that a notable part of the zone adjacent to the front propagates virtually as a simple wave. This implies that the CFL condition of stability for this zone is much less restrictive than for the central part of the fracture. Of essence is also that the wave-like propagation of the near front zone makes preferable upwind schemes of time stepping. On whole, the analysis implies that explicit time stepping may be complementing, competitive and even superior over implicit integration. Tracing changes in the front shape appears merely in 3D problems when the contour changes its form. The problem, being essentially geometrical, it may be solved separately by various methods, including fast marching, level set and the simplest string/marker methods. In 1D cases, it does not arise at all. Quantitative estimations and numerical examples illustrate the theoretical conclusions.

研究の動機と目的

• 水力破砲のシミュレーションにおける計算的課題、特に中央領域の剛性と前線領域の複雑な挙動に起因する問題を解決すること。
• HFモデリングにおける明示的および陰解法の時間積分スキームの性能と安定性を評価すること。
• 形状変化を伴う3次元問題において、破砲輪郭の進化を追跡するための堅牢なフレームワークを構築すること。
• 漸近解析に基づき、近端領域における上流型スキームの理論的・数値的基盤を確立すること。
• 従来の仮定とは対照的に、明示的スキームが特定の破砲領域において陰解法と同等または優れた性能を示す可能性を示すこと。

提案手法

• 動的挙動のモデリングに中心的な役割を果たす速度方程式を強調する修正された水力破砲理論を用いる。
• 近端領域における漸近解析を適用し、普遍的漸近的被覆の概念を導出する。
• 近端領域が単純な波動として振る舞うことが判明し、これによりCFL条件が著しく緩和されることを特定する。
• 波動的挙動の特徴に起因し、近端領域に上流型時間ステッピングスキームを提案する。
• 3次元問題における輪郭追跡には、ファストマーチィング、レベルセット、マーカー／ストリング法などの別個の幾何的手法を用いる。
• 数値的例題と定量的推定を用いて、異なるメッシュ解像度およびスキームにおける理論的結論の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 明示的および陰解法の時間積分スキームは、水力破砲モデリングにおいて、安定性および効率性の観点でどのように比較されるか？
• RQ2 速度方程式は、水力破砲の前線領域挙動を決定づける上で、果たす役割は何か？
• RQ3 安定性制約があるにもかかわらず、明示的スキームは水力破砲モデリングにおいて、陰解法と同等または優れた性能を示せるか？
• RQ4 前線領域の波動的挙動は、数値スキームの選択および時間刻みの制限にどのように影響を与えるか？
• RQ5 3次元水力破砲シミュレーションにおいて、破砲輪郭の進化を追跡するための最も効果的な手法は何か？

主な発見

• 近端領域では波動的挙動が見られ、中央部と比較してCFL条件が著しく緩和されている。
• 速度方程式が近端領域の動的挙動を支配しており、その漸近解析から「普遍的漸近的被覆」の概念が導かれ、異なるモデル間の挙動を統一的に記述できる。
• 近端領域における安定性制約の緩和により、明示的時間積分法が陰解法と同等または優れた性能を示す可能性がある。
• 波動的挙動の特徴に起因し、上流型スキームは近端領域に特に適しており、数値的安定性と効率性が向上する。
• 3次元シミュレーションでは破砲輪郭が時間とともに変化するため、ファストマーチィングやレベルセット法などの別個の幾何的追跡手法が必要となるが、1次元問題ではこのような追跡は不要である。
• 中央部では粗いメッシュで十分であり、精度を損なわず計算コストを削減できる。