Bounding the Null Space: Interval-Based Uncertainty Quantification for Non-Identifiable Groundwater Models

本論文は、区間演算とマコーミック緩和を用いて、非同定な地下水モデルに対してサンプリングフリーの保証された不確実性境界を提供し、同時に、特定の符号および非回転制約を通じて非物理的な回転流の課題に対処する、最適化ベースの境界タイトニング（OBBT）フレームワークを提案するものである。

要約

大きな問題：「失われたパズルのピース」

あなたが地下水の動きを解明しようとしているところを想像してみてください。手元にはいくつかの手がかりがあります。どこで水を汲み上げているか、どこに雨が降ったか、そしていくつかの観測井での水位の測定値です。

しかし、地面は広大であり、手元にある測定値はわずかです。これは「非識別性（non-identifiability）」と呼ばれる問題を引き起こします。それは、物体の角を3つだけ触って、その正確な形を推測しようとするようなものです。その3つの角に完璧にフィットする形状（岩石の種類、流速、水位の組み合わせ）は、何百万通りも存在します。

従来の方法（サンプリング）：
多くの科学者は、推測することでこの問題を解決しようとします。彼らは、地下の状態について少しずつ異なる推測に基づいたシミュレーションを何千回も実行します。そして結果を見て、「よし、水位はおそらく5メートルから10メートルの間だろう」と判断します。

• 欠点： もし1,000回しか推測を行わなかった場合、真の「極端なケース」を見逃してしまう可能性があります。水位は安全（5〜10m）だと考えていても、実際の現実は2〜15mであるかもしれません。十分に推測を行わなかったために、危険を過小評価してしまうのです。

新しい方法：「箱を絞り込む」（OBBT）

著者らは、**最適化に基づく境界絞り込み（OBBT: Optimization-based Bound Tightening）**と呼ばれる、全く異なるアプローチを提案しています。ランダムにシナリオを推測するのではなく、厳格なルールを持つ数学的なパズルとしてこの問題を取り扱います。

比喩：シュリンクラップの箱
考えられる答えが、巨大で透明な段ボール箱の中に浮いていると想像してください。

1. 初期の箱： 最初、私たちは何も知らないため、箱は巨大です。水位は0メートルから100メートルのどこにでもあり得ます。
2. ルールの追加： 次に、物理法則（水は低い方へ流れる）と、実際のデータ（ここで7メートルと測定された）に基づいた「ルール」を追加していきます。
3. 箱を縮める： ルールを追加するたびに、物理的に不可能な部分を箱から切り取っていくことができます。物理的に可能な答えだけが収まるように、箱が最もタイトにフィットするまで、箱を縮め続けます。
4. 結果： 得られるのは「推測のリスト」ではなく、保証された安全な範囲です。水位がこの最終的なタイトな箱の外側になることは決してないと、確信を持って言えるのです。

障害：「壊れたコンパス」

コンピュータ上でこの数学を機能させるために、著者らは複雑な地下水流動の法則を簡略化する必要がありました。彼らは**マコーミック緩和（McCormick relaxations）**と呼ばれる数学的なテクニックを使用しました。

比喩：
地下水の流れを、道路を走る車のように考えてみてください。車（水）は常に、道路の傾斜（下り坂）に従って進まなければなりません。

• 問題： 著者らが計算速度を上げるために数学を簡略化した際、彼らの「コンパス」は壊れてしまいました。その数学では、特定の非常に奇妙な速度と傾斜の組み合わせにおいて、車が「上り坂」を走ることができてしまったのです。
• 結果： 数学がこのような「不可能な」上り坂走行を許容してしまったため、コンピュータは効果的に箱を縮めることができませんでした。コンピュータが「ええと、ここでは水が上に向かって流れる可能性もあるので、何も除外できない」と考えてしまうため、箱は巨大なままになってしまったのです。

解決策：ルールを強制する

著者らは、コンピュータに対して「いいえ、水は上に向かって流れることはできません」と手動で教える必要があると気づきました。彼らは2つの具体的な修正策を追加しました。

1. 流向の符号（Flow Signs）： コンピュータに対し、早い段階で「水は北に向かって流れているのか、それとも南なのか？」を決定するように強制しました。一度この方向が固定されれば、「上り坂」というおかしな現象は消滅します。
2. 渦の禁止（No Swirls）： 水が理由もなく回転（渦巻きのように）することはない、というルールを追加しました。これにより、数学が流れの真の形状を理解できるようになります。

これらの修正を加えることで、ようやく「箱」がタイトに縮まり、信頼できる答えが得られるようになりました。

何をテストしたのか

チームはこの手法を3つの異なるシナリオでテストしました。

1. 1次元のストリップ： 単純なセルの列。完璧に機能し、従来の「推測」による手法よりもはるかに優れていました。
2. 2次元グリッド： 平面図。これによって、「渦の禁止」や「流向の符号」の修正がない場合、この手法が失敗することが示されました。修正を加えると、うまく機能しました。
3. タイムトラベル・グリッド： 時間とともに変化する2次元マップ（ビデオのようなもの）。水位が日ごとに変化していく様子を扱い、時間が経過するにつれて不確実性が減少していくことを示しました。

トレードオフ

良いニュース： この手法は保証された安全性を提供します。十分な回数の推測を行わなかったために、稀に起こる危険なシナリオを見逃す心配はありません。起こり得る事象の絶対的な限界を見つけ出すことができます。

悪いニュース： 計算コストがかかります。単に数千回の推測を実行するのに比べて、これらの数学的パズルを解くには長い時間がかかります。それは、ハサミを使って紙を切る代わりに、レーザーカッターを使って紙をトリミングするようなものです。遅いですが、結果は数学的に完璧です。

まとめ

本論文は、地下水モデルにおける不確実性を扱う新しい方法を提示しています。何千回も推測して最悪のシナリオを捕まえようと願うのではなく、厳格な数学的ルールを使用して、不可能な答えをすべて削ぎ落とし、水位の「保証された安全圏」を残す手法です。これを実現するためには、数学が不可能な水の流れを想像しないように追加のルールを課す必要があることを彼らは発見しましたが、一度そのルールを適用すれば、この手法は従来の推測法よりもはるかに信頼性の高いセーフティネットを提供します。

技術概要

技術要約：最適化に基づく境界タイトニングによる零空間の境界設定

1. 問題提起

地下水モデルは、観測データがまばらであるために、しばしば**非同定性（non-identifiable）**の問題に直面する。この等価性（equifinality）により、異なるパラメータの組み合わせ（例：透水係数、比貯留係数）、状態（水頭）、および境界条件が、同一の観測結果を生み出すことが存在する。

従来の不確実性定量化（UQ）手法は、有限のモデル実現例（例：モンテカルロ法、アンサンブル・データ同化、または正則化インバージョン）を探索することによってこれに対処している。著者らは、このパラダイムには構造的なジレンマがあると考えている。すなわち、システムの複雑さが増すにつれて未知数の数は増加する一方で、シミュレーションの計算コストが探索可能な実現例の数を制限してしまう。その結果、探索が不完全となり、許容される解の真の範囲を系統的に過小評価することにつながる。これは、リスクが稀な極端な事象（例：干ばつや洪水）から生じる地下水管理において、有限のサンプルではそれらを見逃す可能性があるため、特に危険である。

本論文は、保守的なUQを目指すのであれば、目標とすべきは「実現可能領域をサンプリングすること」ではなく、「不確確定な変数に対する保証された外側境界を通じて、実現可能領域の広がりそのものを表現すること」であると主張している。

2. 手法：最適化に基づく境界タイトニング (OBBT)

提案されたフレームワークは、サンプリングに代わって**最適化に基づく境界タイトニング（OBBT）**を採用する。実現例を列挙する代わりに、この手法は不確実性を区間として扱い、物理法則と観測をエンコードした制約系に基づいて変数を極値化することで、それらをタイトニング（絞り込み）する。

2.1 数学的定式化

地下水流動問題を以下の制約系のシステムとして定式化する：

• 変数: 水頭 ($h$)、フラックス ($q$)、涵養量 ($R$)、透水係数 ($T$)、および比貯留係数 ($S$)。
• 制約:
  • 質量保存: 有限体積法を用いて離散化（定常および非定常）。
  • ダルシーの法則: $q = -T \nabla h$。
  • 観測: 特定の場所における固定値または境界。
  • 事前境界: すべての変数に対する初期区間。

2.2 非線形性の処理：マコーミック緩和 (McCormick Relaxations)

支配方程式には双線形項（例：不確実な透水係数 $T$ と水頭勾配 $\nabla h$ の積）が含まれる。大規模なグリッドに対して、結果として生じる非凸最適化問題をグローバルに解くことは計算量的に困難である。

計算の実行可能性を維持しつつ保守的な境界を確保するために、著者らはマコーリック緩和を採用している。これは、双線形等式 $w = xy$ を、真の非線形実現可能領域の周囲に形成される凸包（エンベロープ）となる4つの線形不等式に置き換えるものである。

• 利点: 効率的な線形計画法（LP）ソルバーの使用を可能にする。
• 欠点: 緩和によって非物理的な解が許容される。具体的には、フラックスと水頭勾配の間の符号の結合が崩れる。これにより、「回転流」（勾配に逆らって水が循環するサイクル）が許容され、これがダルシー流の物理学に反し、効果的な境界タイトニングを妨げる。

2.3 アルゴリズム

OBBTアルゴリズム（アルゴリズム1）は、以下の手順で反復的に動作する：

1. 初期化: 変数を時空間グリッドグラフ上に定義し、初期境界を設定する。
2. 緩和の構築: 現在の変数境界に基づき、マコーリック制約を構築する。
3. 極値化: 各変数について、完全な制約系に制約された2つのLP（最小化および最大化）を解く。
4. タイトニング: 新しい極値を用いて変数境界を更新する。
5. 後処理: 区間演算を用いて、フラックスと勾配の間の符号の一貫性を強制する。
6. 反復: 収束（境界のタイトニングが停止）するか、最大反復回数に達するまで繰り返す。

このプロセスは、各反復内での個々の変数の極値化に関して、**タブラサス・パラレル（並列化が容易）**である。

3. 主要な貢献と救済策

本論文は、マコーリック緩和を用いた単純なOBBTの適用は、緩和が回転流を許容してしまうため、有益な境界を提供できないことを指摘している。著者らは、物理的一貫性を回復するために、2つの具体的な救済策を提案し評価している。

1. 無回転制約 (Irrotationality Constraints): グリッドグラフ内の任意の基本サイクルにおけるフラックスの総和がゼロであることを明示的に強制する。

   • 有効性: 均質媒体において非常に効果的であり、有益な境界を回復させる。
   • 限界: ゼロの純フラックスが厳密に無回転性から導かれない、不均質または異方性媒体には一般に適用できない。

2. 流向の符号指定 (Flow Sign Prescription): リファレンス・シミュレーションまたは事前知識に基づいて、各インターフェースにおける流向（方向）の符号を事前に指定する。

   • 有効性: マコーリック・エンベロープ内の非物理的な象限を除去し、境界を大幅にタイトニングする。
   • 汎用性: 無回転制約よりも不均質設定において堅牢であるが、流向に関する事前仮定を導入することになる。

4. 結果

本フレームワークは、3つの数値例でテストされた。

4.1 1次元定常モデル

• 設定: ソース/シンク・ウェルと2つの観測点を持つ10セル。
• 比較: OBBTの境界を、正確なモンテカルロ・サンプラー（真の零空間多様体から抽出）およびMCMCと比較した。
• 知見:
  • OBBTは零空間多様体を包含することに成功し、保証された外側境界を提供した。
  • MCMCおよび有限サンプル手法は、実現可能領域の広がりを著しく過小評価した（例：10万個のサンプルを用いても、周辺範囲の約45%しかカバーできていない）。
  • OBBTは2回の反復で収束し、この非同定的なシステムにおいて、境界を捉えるためにサンプリングは不要であることを示した。

4.2 2次元定常モデル

• 設定: ソース、シンク、および中央観測点を持つ5x5グリッド。
• 知見:
  • 基本OBBT (追加制約なし): 符号の曖昧さが回転流を許容するため、水頭や透水係数の境界タイトニングに失敗した。
  • 無回転性を伴うOBBT: 流向を解決し、境界をタイトニングすることに成功し、純粋なLPリファレンス解に近似した。
  • 流向指定を伴うOBBT: これも境界をタイトニングしたが、この均質設定においては無回転の場合よりもわずかに広い境界となった。これは、仮定の強さと境界のタイトネスの間のトレードオフを示している。

4.3 2次元非定常モデル

• 設定: 六角形グリッド、5タイムステップ、涵養ゾーン、抽出ウェル、および固定水頭境界。
• 知見:
  • OBTBは時変システムを正常に処理し、降下錐（drawdown cone）の発達と、時間を遡った制約の伝播を捉えた。
  • 水頭の境界は大幅にタイトニングされたが、透水係数の境界は（制約が不足しているパラメータ空間において予想通り）わずかにタイトニングされるにとどまった。
  • 計算コスト: 16コアで約71.7時間を要した。著者らは、これは相当なものだが、サンプリング手法における過小評価のリスクと比較して、「保守的なカバレッジのための代償」であると述べている。

5. 重要性と主張

本論文は、OBBTがアンサンブルベースのUQに代わる、保守的、決定的、かつ物理的根拠に基づいた代替案を提供すると主張している。その主な意義は、サンプリングを行うことなく、すべての不確実な変数に対して保証された外側境界を提供できる点にあり、これによりモンテカルロ法に固有の「探索不足による過小評価」という問題を回避できる。

著者らは、OBBTは現在のサンプリング手法よりも計算コストが高いものの、提案された抽出率が「あらゆる許容されるパラメータの組み合わせの下で」安全であることを規制当局が証明しなければならない意思決定支援シナリオにおいては、そのトレードオフは正当化されると強調している。この手法は零空間理論と自然に結びついており、点推定ではなく区間境界を通じて「零空間多様体」（等価な解の集合）を特徴付けている。

結論として、現在の実装は、ダルシー流の非凸性を扱うために特定の仮定（流向の符号指定など）に依存しているものの、このフレームワークは、より効率的な混合整数定式化の開発や、データ同化ワークフローへの統合といった、将来の研究に向けた有望な道筋を示すものであるとしている。