【研究領域課題番号】20K21785 (KAKENデータベースで見る）

【研究キーワード】

データ同化 / 4次元変分法 / マルコフ連鎖モンテカルロ法 / レプリカ交換モンテカルロ法 / 地震波伝播 / 不確実性評価 / 粒子フィルタ / Ginzburg-Landau方程式 / フェーズフィールドモデル

【研究成果の概要】

理工学分野において、物理・化学等の理論に基づくモデルと観測・実験に基づくデータを比較することは、対象の理解を深め、かつその将来予測をする上で基本的なアプローチであるが、近年急速に進んでいる数値モデルの大規模化と観測データの大容量化に応じた比較手法の高度化が急務となっている。大規模モデルと大容量データをベイズ統計学に基づいて融合する計算技術であるデータ同化は、数値モデルからあらかじめ生成した複数のシナリオを観測データが入力されるたびに再構成する「逐次型データ同化」と、数値モデルから生成した単一のシナリオを観測データに適合するように最適化する「非逐次型データ同化」とに大別される。大規模データ同化を実施する場合には、計算機資源や計算時間の制約の観点からも圧倒的に非逐次型データ同化の方が有効であるが、数値モデルのコードに大掛かりに手を加える必要があり、実装の手間の観点から採用が回避される傾向にある。一方、逐次型データ同化は実装の手間が比較的小さく、また数値モデルをブラックボックスのまま取り扱うことも可能であるが、大規模モデルと大容量データを融合するためには大量のシナリオ数を必要とする。
本研究の最終目標は、現実的な計算時間と計算機資源で実行可能な、逐次型と非逐次型を融合したハイブリッド型データ同化手法の開発に挑戦することである。2021年度は、前年度までに開発した地震波動伝播シミュレーションモデルの高速化を行なった上で、非逐次データ同化の代表的手法である4次元変分法を実装することにより、地下構造の物性パラメータの最適化ならびにその不確実性評価を実施可能なコードを開発した。現在、逐次データ同化の一種であるレプリカ交換モンテカルロ法を実装することにより、4次元変分法における反復計算の効率化を図っており、今後のハイブリッド型データ同化手法の確立と実モデルへの適用に向けた見通しが立った。

【研究種目】挑戦的研究(萌芽)

【研究期間】2020-07-30 - 2023-03-31

【配分額】6,500千円 (直接経費: 5,000千円、間接経費: 1,500千円)

【研究分野】統計科学

【研究領域課題番号】17H01703 (KAKENデータベースで見る）

【研究キーワード】

データ同化 / 4次元変分法 / 最適化 / 不確実性評価 / 地震 / 時系列解析 / モデリング / Ginzburg-Landau方程式 / マルコフ連鎖モンテカルロ法 / 地震波 / 摩擦

【研究成果の概要】

次世代型4次元変分法のプロトタイプ版を開発し、2次元空間内における地震波伝播シミュレーションコード等に実装した。双子実験によって地下構造に関する物性パラメータの最適値推定と不確実性評価が可能となることを示した。また、本手法を実装する新たな対象として、磁性体内における磁化パターンの時空間変化を記述する時間依存Ginzburg-Landau方程式の数値シミュレーションコードを新たに開発した。これらの成果については、国際誌ならびに国内外の学会において発表を行った。
本研究で開発した次世代型4次元変分法のソースコードについては、現在、公開に向けた準備を進めている。

【研究の社会的意義】

本研究で開発した次世代型4次元変分法により、従来型4次元変分法では不可能であった広域最適解探索と不確実性評価を一挙に解決することが可能となり、大規模数値モデルに基づくデータ同化の関連分野、特に気象予報業務や地震発生予測研究にとっては大きな転回点となるであろう。特に地震学においては、地震発生直後の迅速な地震波動場推定や、巨大地震発生メカニズムの理解に直結するプレート境界における断層すべりの時空間モデルの構築など、将来的に地震防災・減災に資する手法となる可能性が本研究によって拓かれた。

【研究種目】基盤研究(B)

【研究期間】2017-04-01 - 2020-03-31

【配分額】15,860千円 (直接経費: 12,200千円、間接経費: 3,660千円)

【研究分野】統計科学

【研究領域課題番号】26280006 (KAKENデータベースで見る）

【研究キーワード】

データ同化 / 固体地球科学 / 逐次ベイズフィルタ / スパースモデリング / データ駆動型モデリング / マルコフ連鎖モンテカルロ法 / レプリカ交換モンテカルロ法 / 地震 / 疎性モデリング / ベイズ統計学 / 火山 / プレート境界

【研究成果の概要】

本研究課題では、固体地球科学に資する次世代型データ同化手法の創出という所期の目的を、大きく進歩させることができた。特に，大自由度モデルにも適用可能なデータ同化手法である4次元変分法の高度化や、レプリカ交換モンテカルロ法との融合によるモデル／データ両駆動型データ同化の地震観測網データへの応用研究を実施した。また、国内および海外における多数の学会において成果発表を行うとともに，国際誌に論文を出版した。これらの成果をまとめた教科書を出版する予定である。さらに、日本地球惑星科学関連学会連合大会を始めとする学会において、本研究課題に関連するセッションを開催した。

【研究種目】基盤研究(B)

【研究期間】2014-04-01 - 2018-03-31

【配分額】16,380千円 (直接経費: 12,600千円、間接経費: 3,780千円)

【研究分野】感性情報学・ソフトコンピューティング

【研究領域課題番号】23700278 (KAKENデータベースで見る）

【研究キーワード】

データ同化 / クラウドコンピューティング / 地震音波 / 南海トラフ地震 / 粒子フィルタ / マルコフ連鎖モンテカルロ / 事後分布 / 微気圧 / ベイズ統計学 / 並列計算 / 細胞質流動 / ライブラリ / グラフィカルユーザインターフェース / ベイズ統計

【研究成果の概要】

本研究では、ベイズ統計学の枠組みで数値シミュレーションモデルと観測データを融合するための計算基盤技術であるデータ同化を、クラウドコンピューティングを利用して実行するためのデータ同化システムを開発した。大地震に伴う地震音波伝播を題材に双子実験を実施したところ、各震源パラメータに関する周辺化事後分布は正しく真の値の周辺に分布し、かつ事後分布の最大値を与えるパラメータから再現された理論波形と擬似観測波形は、明らかに一致することが確認された。これにより、本研究で開発したデータ同化システムが適切に動作していることが実証された。

【研究種目】若手研究(B)

【研究期間】2011-04-28 - 2015-03-31

【配分額】4,420千円 (直接経費: 3,400千円、間接経費: 1,020千円)