[論文レビュー] Modelling and Numerical Simulation of Gas Migration in a Nuclear Waste Repository

Alain Bourgeat, Mladen Jurak+1 more|arXiv (Cornell University)|Jun 15, 2010

Hydraulic Fracturing and Reservoir Analysis参考文献 31被引用数 5

ひとこと要約

本稿では、深地中核燃料廃棄物処理施設における水素移動を対象として、全水素質量濃度と液相圧力を主変数として用いる、特異性のない構成的二相流モデルを提案する。この定式化は、飽和および不飽和状態の両方において相の出現・消失を的確に再現でき、数値的テストにより鋭い飽和フロントや動的平衡遷移の取り扱いにおいて高い安定性を示している。

ABSTRACT

We present a compositional compressible two-phase, liquid and gas, flow model for numerical simulations of hydrogen migration in deep geological radioactive waste repository. This model includes capillary effects and the gas diffusivity. The choice of the main variables in this model, Total or Dissolved Hydrogen Mass Concentration and Liquid Pressure, leads to a unique and consistent formulation of the gas phase appearance and disappearance. After introducing this model, we show computational evidences of its adequacy to simulate gas phase appearance and disappearance in different situations typical of underground radioactive waste repository.

研究の動機と目的

• 核燃料廃棄物処理施設のシミュレーションにおいて、多孔質媒体内での気相反応の出現および消滅をモデル化する課題に取り組む。
• 相飽和方程式における数値的特異性の問題を解消するため、非特異的変数を用いた定式化によるシステム再構築を実施する。
• 飽和および不飽和な流れ領域をまたがる水素移動を一貫的かつ安定的にシミュレートすることを保証する。
• 毛管効果、気体拡散、圧力不連続性といった複雑な物理現象の取り扱い能力をモデルが有するかを検証する。
• 地下放射性廃棄物処理施設に一般的に見られる、現実的で困難なシナリオへのモデル適用可能性を示す。

提案手法

• 水と水素（H₂）を成分とする液体相および気相反応の、構成的かつ圧縮性の二相流モデルを定式化する。
• 飽和方程式における特異性を回避するため、全水素質量濃度と液相圧力を主変数として採用する。
• 毛管圧力関係および理想気体法則を気相反応に適用し、液体相における溶解水素の拡散を組み込む。
• Darcy-Muskat法則を用いて相の流れ速度を定式化し、相対透過度関数を相飽和度に依存させる。
• 飽和フロントにおける不連続性や急勾配に対処できるよう、支配方程式の弱形式を採用する。
• 動的状態遷移を解像するため、時間ステップと空間メッシュの適応的制御を用いた数値的シミュレーションを実施する。

実験結果

リサーチクエスチョン

• RQ1 気相反応の出現および消滅に伴う数値的安定性を保つ二相流モデルは、どのように定式化できるか？
• RQ2 飽和および不飽和状態をまたがる飽和に基づく方程式において、一貫性を保ち特異性を回避するにはどのような変数選択が有効か？
• RQ3 毛管効果および溶解水素拡散は、気相反応度および圧力プロファイルの時間的変化にどのように寄与するか？
• RQ4 移動する飽和フロントを伴う飽和状態から不飽和状態への遷移を、モデルが正確に再現できるか？
• RQ5 対流的および拡散的フラックスは、飽和フロント付近における水素濃度勾配の形状に、どのように寄与するか？

主な発見

• モデルは、1.92 × 10⁵ s の間に x = 0.5 m から左境界（Γ_in）へ向かって気相反応フロントが移動することを的確に再現しており、フロント高さは約 0.16 から約 0.03 へと減少した。
• 液相圧力は、飽和領域で急激に上昇（t = 10³ s 時点で約 1 MPa から約 1.6 MPa へ）し、その後ゆっくり減少し、上層の気相反応による動的圧縮を示している。
• 溶解水素濃度に強い勾配が飽和フロント付近に発生し、t = 1.92 × 10⁵ s まで持続した。これは、左から右への対流的輸送と反対方向の拡散フラックスとの競合によるものである。
• t = 1.92 × 10⁵ s を過ぎると、全領域が不飽和状態（S_g > 0）となり、圧力、水素密度、飽和度といったすべての変数が空間的に一様な定常状態へ単調に収束した。
• 系は、閉じた系内での質量保存則と整合的な連続的な圧力および濃度プロファイルを介して平衡状態に達した。
• モデルは、圧力連続性の即時の回復と、濃度勾配の徐々な平滑化を的確に捉えており、相転移ダイナミクスの取り扱いにおけるモデルの強固さが裏付けられた。