| 研究課題/領域番号 |
21K03686
|
|---|
| 研究種目 |
基盤研究(C)
|
| 配分区分 | 基金 |
|---|
| 応募区分 | 一般 |
|---|
| 審査区分 |
小区分17030:地球人間圏科学関連
|
|---|
| 研究機関 | 国立研究開発法人防災科学技術研究所 |
|---|
研究代表者 |
澤崎 郁 国立研究開発法人防災科学技術研究所, 地震津波火山ネットワークセンター, 主任専門研究員 (30707170)
|
|---|
| 研究期間 (年度) |
2021-04-01 – 2024-03-31
|
| 研究課題ステータス |
交付 (2022年度)
|
| 配分額 *注記 |
4,160千円 (直接経費: 3,200千円、間接経費: 960千円)
2023年度: 1,430千円 (直接経費: 1,100千円、間接経費: 330千円)
2022年度: 1,300千円 (直接経費: 1,000千円、間接経費: 300千円)
2021年度: 1,430千円 (直接経費: 1,100千円、間接経費: 330千円)
|
|---|
| キーワード | 極値統計解析 / 連続地震波形記録 / 最大振幅 / 超過確率 / マルコフ連鎖モンテカルロ法 / 区間最大振幅 / Hi-net / MeSO-net / 地震活動 / Frechet分布 |
|---|
| 研究開始時の研究の概要 |
本研究では、これまでに蓄積されてきた大量の地震計記録に極値統計解析を適用し、地震動記録の統計的性質のみに基づき、各地震観測点での地震動の活動度を定量的に評価する。その評価指標を外挿して2050年頃までの最大地震動の超過確率を算出し、実際の観測記録の傾向と従来の地震動予測法に基づく予測結果と照らし合わせ、地震動予測の高度化に資する情報を得ることを目的とする。
|
| 研究実績の概要 |
昨年度構築した極値統計解析のアルゴリズムにおいては、パラメータ推定値の分布関数として正規分布を仮定していた。しかしその仮定は厳密には成り立たないため、本年度は、パラメータの分布関数をマルコフ連鎖モンテカルロ（MCMC）法に基づき直接推定し、その各々の推定値から最大振幅の超過確率予測を行うアルゴリズムを開発した。また、任意の値以上の最大振幅が任意の回数以上発生する確率の計算もアルゴリズムに組み込んだ。このアルゴリズムを、200点余りのMeSO-net観測点において、10年以上の期間に記録された区間最大振幅に対し適用した。
その結果、30年以内に1000Gal以上の揺れを観測する確率は、99%以上の点から0.01%以下の点まで、観測点によりかなりばらついた。観測点ごとのサイト特性や地震活動の地域性を考慮しても、このばらつきは地震学的に不自然である。大きなばらつきの原因としては、2011年東北地方太平洋沖地震後の活発な地震活動や設置点周囲のノイズ環境の変化などにより、区間最大振幅がしたがう分布の定常性が保証できないことや、MCMC法のアルゴリズムの不備により、局所解周辺を探索している可能性があることなどが挙げられる。使用するデータセットの厳選と局所解に陥らないようなアルゴリズムの改良が今後必要である。
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
極値統計解析のアルゴリズム構築においては、MCMC法の導入により、パラメータ推定値の分布を考慮した予測が可能となった。また、任意の値以上の揺れが発生する回数も計算できるようになった。これらの点から、アルゴリズム開発においては大きな進展があったと考える。一方で、実データへの適用において、データの蓄積期間が長くなるほど区間最大振幅がしたがう分布の定常性が損なわれやすく、このことが予測結果に大きく影響する可能性を無視できないことが明らかとなった。そのため、定常性を表す客観的な指標を導入し、予測そのものの信頼度を示す必要が生じたため、予測結果の解釈には至っていない。
以上のことから、全体としてはおおむね順調に進展していると判断した。
|
| 今後の研究の推進方策 |
今後は、区間最大振幅データセットの定常性を客観的に見積もる手法を考案し、その程度に応じて、予測結果の信頼性をランク付けしたいと考えている。その上で、30年以内に任意の値以上の揺れが生じる確率を可視化し、その結果を確率論的地震動予測地図などの従来法による予測結果と比較し、その共通点や相違点などについて考察を進めたい。また、最大振幅に加えて計測震度の区間最大値のデータベースもすでに準備済みであるため、計測震度の超過確率を計算するアルゴリズムも構築したい。
|
