[論文レビュー] Multilevel Monte Carlo methods for highly heterogeneous media

Aretha L. Teckentrup|arXiv (Cornell University)|Jun 7, 2012

Probabilistic and Robust Engineering Design参考文献 17被引用数 5

ひとこと要約

この論文は、地下水流モデリングで一般的な対象となる、非常に不均質で非一様に強制的なランダム係数（例：対数正規分布の透水係数）を有する楕円型PDEに、多層モンテカルロ（MLMC）法を拡張する。ランダム場の正則性に関する最小限の仮定の下で、テンソル値係数および圧力・ダルシー流速の点評価に対する収束性と計算コストの上限を証明する。

ABSTRACT

We discuss the application of multilevel Monte Carlo methods to elliptic partial differential equations with random coefficients. Such problems arise, for example, in uncertainty quantification in subsurface flow modeling. We give a brief review of recent advances in the numerical analysis of the multilevel algorithm under minimal assumptions on the random coefficient, and extend the analysis to cover also tensor--valued coefficients, as well as point evaluations. Our analysis includes as an example log--normal random coefficients, which are frequently used in applications.

研究の動機と目的

• 非常に不均質な透水係数を有する地下水流における不確実性定量化において、標準モンテカルロ法の高い計算コストを解決すること。
• ランダム係数が一様に強制的でない、あるいは有界でない場合（例：対数正規分布）に、多層モンテカルロ（MLMC）の理論的分析を拡張すること。
• 地下水流における異方性媒体のモデリングに適した、テンソル値係数を扱えるようにMLMCフレームワークを一般化すること。
• 圧力またはダルシー流速の点値評価を行う場合のMLMC推定子の収束性と計算コストの上限を確立すること。
• 解およびその有限要素近似のパスワイズ正則性とモーメントバウンドを最小限の仮定の下で厳密に誤差解析すること。

提案手法

• 地下水流における透水係数をモデル化するランダムテンソル値係数 $\mathbf{A}$ を有する楕円型PDE $-\mathrm{div}(\mathbf{A}\nabla u) = f$ に多層モンテカルロアルゴリズムを適用する。
• メッシュサイズ $h_\ell$ のネストされたメッシュの階層上で有限要素法を用いて解 $u$ を近似し、汎関数 $M$ に対して推定子 $Q_{h_\ell} = M(u_{h_\ell})$ を定義する。
• 期待値 $\mathbb{E}[Q]$ をレベル間の差分 $\mathbb{E}[Q_{h_\ell} - Q_{h_{\ell-1}}]$ を用いたテレスコピック和に分解し、分散を低減する。
• 解 $u(\omega,\cdot)$ のSobolev埋没およびSchauder型正則性推定を用いて、バイアス $\mathbb{E}[|Q - Q_{h_\ell}|]$ の収束速度を確立する。
• ランダム場のホルダー正則性および可積分性に関する仮定を用いて、解およびその有限要素誤差のモーメントバウンドを導出する。
• 最小限の仮定の下で、ルート平均二乗誤差 $\varepsilon$ を達成するための計算コストが $\mathcal{O}(\varepsilon^{-2})$ で上限づけられることを証明する。これは標準モンテカルロの最適レートを達成するが、定数が著しく小さい。

実験結果

リサーチクエスチョン

• RQ1 一様に強制的でない、あるいは有界でないランダム係数を有する楕円型PDEに対して、多層モンテカルロ法を厳密に解析可能か？
• RQ2 異方性媒体に用いられるテンソル値ランダム透水係数に対して、MLMC法はその最適な複雑度を維持するか？
• RQ3 グローバル汎関数ではなく、圧力またはダルシー流速の点評価に対して、収束性および計算コストの上限を確立できるか？
• RQ4 解が粗いランダム係数のため空間正則性が低い場合、MLMC推定子の収束速度はどのようになるか？
• RQ5 パスワイズ正則性推定（例：Schauder理論）は、MLMCフレームワークにおけるバイアスおよび分散の減衰にどのように寄与するか？

主な発見

• ランダム係数に関する最小限の仮定の下で、多層モンテカルロ推定子はルート平均二乗誤差 $\varepsilon$ を $\mathcal{O}(\varepsilon^{-2})$ の計算コストで達成でき、標準モンテカルロの最適レートを達成するが、定数が著しく小さい。
• 圧力またはダルシー流速の点評価に対しては、バイアスが $\alpha = 1 + s - d/2$（$s < 1$）またはSchauder理論からの $\alpha = 1 + t$ のレートで減少し、分散は $\beta = 2\alpha$ のレートで減少する。
• 一様に強制的でない対数正規ランダム場（コバリアンスがリプシッツ連続に限る）に対しても解析が成立し、従来の均一有界性を要請する結果を拡張する。
• 有限要素誤差の収束レートは $\|u - u_h\|_{L^p(\Omega, W^{1,\infty}(D))} \lesssim h^{s - d/2}$（$s < 1$）であり、Schauder正則性の下ではより良いパスワイズレート $h^t$ も達成可能である。
• テンソル値係数（例：異方性透水係数）に対しても、一様強制性や有界性を要件とせず、適用可能である。
• 理論的枠組みは、空間的変動が大きく、正則性が低い透水係数を有する実用的な地下水流シミュレーションにおけるMLMCの使用を支援する。