| Project/Area Number |
20K21785
|
|---|
| Research Category |
Grant-in-Aid for Challenging Research (Exploratory)
|
| Allocation Type | Multi-year Fund |
|---|
| Review Section |
Medium-sized Section 60:Information science, computer engineering, and related fields
|
|---|
| Research Institution | The University of Tokyo |
|---|
Principal Investigator |
長尾 大道 東京大学, 地震研究所, 准教授 (80435833)
|
|---|
| Project Period (FY) |
2020-07-30 – 2025-03-31
|
| Project Status |
Granted (Fiscal Year 2023)
|
| Budget Amount *help |
¥6,500,000 (Direct Cost: ¥5,000,000、Indirect Cost: ¥1,500,000)
Fiscal Year 2021: ¥3,250,000 (Direct Cost: ¥2,500,000、Indirect Cost: ¥750,000)
Fiscal Year 2020: ¥3,250,000 (Direct Cost: ¥2,500,000、Indirect Cost: ¥750,000)
|
|---|
| Keywords | データ同化 / 4次元変分法 / マルコフ連鎖モンテカルロ法 / レプリカ交換モンテカルロ法 / 地震波伝播 / 不確実性評価 / 粒子フィルタ / Ginzburg-Landau方程式 / フェーズフィールドモデル |
|---|
| Outline of Research at the Start |
理論に基づくモデルと、観測・実験に基づくデータを比較することが、対象の特性の把握や予測をする上で極めて重要であることは論をまたない。今日の大規模モデル・大容量データ時代において、両者の融合という必然的な要請に応える計算技術が「データ同化」である。データ同化は逐次型と非逐次型とに大別され、前者は大規模モデルに適しているが実装の手間が大きく、後者は実装は容易だが中・小規模モデルにしか適さない。そこで本研究では、現実的な計算時間と計算機資源で実行可能な、逐次・非逐次融合型データ同化手法の開発に挑戦する。また、結晶成長モデルおよび地震波動伝播モデルへ実装することにより、その性能を評価する。
|
| Outline of Annual Research Achievements |
本研究の最終目標は、現実的な計算時間と計算機資源で実行可能な、逐次型と非逐次型を融合したハイブリッド型データ同化手法の開発に挑戦することである。逐次型データ同化とは、数値モデルからあらかじめ生成した複数のシナリオを観測データが入力されるたびに再構成する手段であり、非逐次型データ同化とは、数値モデルから生成した単一のシナリオを観測データに適合するように最適化する手段である。2022年度までに、ハイブリッド型データ同化手法として4次元変分法に乱択アルゴリズムによる前処理を加えたものを提案し、大規模地震波動場データ同化における不確実性評価において、前処理がない既存方法よりも計算時間を10%以下に抑えることが可能であることを示した。
2023年度は、波動方程式特有の数理構造を利用することで、フォワードモデルとアジョイントモデルの計算量および必要メモリ量を大幅に削減する手法を開発した。具体的には、対象の波動方程式に対応するアジョイントモデルが"局所的にほぼ"線形系であることに着目し、元来非線形系では計算量的に利用困難な陰的解法を線形系と同程度の計算量で計算可能なアルゴリズムを開発し、適用することで計算不安定性を回避する。これにより時間刻みの縛りから解放され、時間刻みを大きくとれるようになったことにより、大幅な計算量およびメモリ量の削減を実現した。本提案手法を、現実の地下構造を考慮し地下方向へ完全吸収境界条件を課した１次元モデルへ適用したところ、既存の手法に比べ数十倍程度の高速化が実現された。
|
| Current Status of Research Progress |
Current Status of Research Progress
3: Progress in research has been slightly delayed.
Reason
本研究課題の目標であるハイブリッド型データ同化手法のアルゴリズムは、ほぼ完成しつつある。ただし、地震波伝播シミュレーションコードへの開発手法の実装を担当する研究員が2021年度上半期に退職した影響で、予定していたスケジュールからは1年以上遅れた状態となっている。その後、新たな研究協力者として研究員を雇用し、また高速化作業の一部を役務として外注することにより、本研究計画の継続的な推進に務めている。
|
| Strategy for Future Research Activity |
本研究課題の最終年度となる2024年度は、これまでに開発したハイブリッド型データ同化手法、および1次元の地震波伝播シミュレーションコードへの適用例を示した内容について、国際誌に投稿する論文としてまとめる。現実的な地震波伝播シミュレーションコードへの実装については、時間を要することが予想されるため、別プロジェクトにおいて継続実施していく。
|
