|
本年度は津波波源インバージョンの事前準備として必要なグリーン関数データベースの計算効率を相反原理に基づいて効率的に行う方法の背景となる理論を導出し，相反原理の成立条件を精査した．その結果，相反原理が成立しない方程式系でもアジョイントモデルの開発によりグリーン関数を逆向きに計算することでデータベースの計算効率を飛躍的に高められることがわかった．この方法により，高解像度のデータベースが容易に得られるようになり，これまで困難であった即時的な矩形断層モデルの推定が可能になった．マグニチュード9クラスの地震津波では，断層面が広く，面内の滑り量分布の多様性が大きいことから，様々な形状の津波波源を表現可能な自由度の高い津波波源推定手法が最適である．しかしながら，マグニチュード7クラスの地震津波では波源域が相対的に狭いために，多くの観測点で早期に波形を得ることが難しく，自由度の高い波源推定手法では推定結果が不安定になったり，解像度が極端に低くなったりする問題があった．一方で，マグニチュード7クラスの津波では一様な滑り分布を持つ矩形断層モデルでも津波の波形予測に十分な精度を有することが報告されており，自由度が少ない分，安定した推定が行えることが期待できる．このような観点から本年度は均一滑りの矩形断層モデルパラメータを津波波形から推定するプログラム開発を行った．断層パラメータの推定問題は複数の局所解が存在し，全体最適化が困難な問題であるが，レプリカ交換マルコフ連鎖モンテカルロ法を用いることで安定的に正解値が得られることを数値的な双子実験により確認した．
|
