[論文レビュー] Impact of background flow on dissolution trapping of carbon dioxide injected into saline aquifers

Saikiran Rapaka, Rajesh Pawar|arXiv (Cornell University)|Jan 13, 2011

CO2 Sequestration and Geologic Interactions参考文献 15被引用数 6

ひとこと要約

本研究は、背景地下水の流れが塩水層におけるCO2の溶解をどのように促進するかを調査しており、フリーコンベクション（フィンガリング）からフォースドコンベクションへの移行に注目している。理論的分析と高解像度の数値シミュレーションを用いて、システムがフォースドコンベクションに移行する臨界背景流速を同定し、特にフィンガリングが遅いか効果的でない低レイノルズ数の層では、溶解速度が顕著に向上することを示している。

ABSTRACT

While there has been a large interest in studying the role of dissolution-driven free convection in the context of geological sequestration, the contribution of forced convection has been largely ignored. This manuscript considers CO$_2$ sequestration in saline aquifers with natural background flow and uses theoretical arguments to compute the critical background velocity needed to establish the forced convective regime. The theoretical arguments are supported by two dimensional high-resolution numerical simulations which demonstrate the importance of forced convection in enhancing dissolution in aquifers characterised by low Rayleigh numbers.

研究の動機と目的

• 背景地下水の流れが塩水層におけるCO2の溶解効率に与える影響を定量化すること。
• 低レイノルズ数系において、フリーコンベクションからフォースドコンベクションに移行する臨界背景流速を特定すること。
• 低対流環境において、フィンガリングに起因する遅い溶解と比較して、フォースドコンベクションが溶解速度を顕著に向上させるかどうかを評価すること。
• 地質的CO2隔離における溶解のスケーリング法則に対する理論的および数値的妥当性を検証すること。
• 流れに起因する溶解ダイナミクスの理解を通じて、モニタリング井戸の配置やリザポアシミュレーションの設計を支援すること。

提案手法

• 対流速度、深さ、溶解度スケールを用いたダルシーの法則、質量保存則、拡散移流方程式の次元なし化。
• ブッシネスク近似を適用し、CO2溶解に起因するブリネの密度変化をモデル化。密度変化は浮力にのみ寄与する。
• 浮力力と粘性力・拡散力の比としてのレイノルズ数（Ra）を定義し、対流状態の遷移を支配する。
• スケーリング解析に基づき、フォースドコンベクションが支配的になる臨界背景流速（u_c）の理論的導出。
• 次元なし方程式を用いた高解像度2次元数値シミュレーションにより、ヌセルト数（Nu）を計算し、時間経過に伴う溶解ダイナミクスを追跡。
• フリーコンベクション状態でフィンガリングを維持するために、定期的にシミュレーションを強制的に駆動し、人工的な定常状態を回避。

実験結果

リサーチクエスチョン

• RQ1 塩水層におけるCO2の溶解がフリーコンベクションからフォースドコンベクションに移行する臨界背景流速は何か？
• RQ2 低レイノルズ数系において、フォースドコンベクションはフリーコンベクションと比較してどのように溶解速度を向上させるか？
• RQ3 理論的スケーリング法則が、数値シミュレーションにおける溶解速度をどの程度正確に予測できるか？
• RQ4 背景流は、溶解CO2を検出するためのモニタリング井戸の配置の信頼性にどのように影響するか？
• RQ5 境界層が解像されない大規模なリザポアシミュレータへのスケーリングにおいて、フォースドコンベクションの意味は何か？

主な発見

• 理論的に予測された臨界背景流速は、$ u_c = 2.5 $（$ Ra = 400 $、$ L = 10 $）であり、シミュレーションがこの値を非常に高い精度で確認した。
• $ u_0 > u_c $ の場合、システムはフォースドコンベクションの挙動を示し、安定した溶解速度を示すが、$ u_0 < u_c $ の場合、フィンガリングに起因する一時的な最小値が観察され、状態遷移が確認された。
• 次元なしのヌセルト数 $ Nu/RaL $ は理論的予測に従う：フォースド状態では $ Nu/RaL = \sqrt{u_0/RaL} $、フリーモードでは定数となる。
• 低レイノルズ数系（例：$ Ra = 400 $）では、フォースドコンベクションが、遅くフィンガリングに支配されるフリーコンベクションと比較して、溶解速度を顕著に向上させる。
• シミュレーションでは、フォースドコンベクションが安定的かつ予測可能な溶解速度をもたらすことが示され、大規模なリザポアモデルへのスケーリングにおいて有利である。
• 流動的拡散および水平方向透水率の異方性が、フォースドコンベクションの効率をさらに高め、臨界流速を低下させ、理想化されたモデルを上回る溶解速度を実現する。