| 研究課題/領域番号 |
21K03722
|
|---|
| 研究機関 | 京都大学 |
|---|
研究代表者 |
澁谷 拓郎 京都大学, 防災研究所, 教授 (70187417)
|
|---|
| 研究期間 (年度) |
2021-04-01 – 2025-03-31
|
| キーワード | レシーバ関数 / フルウェーブインバージョン / 3次元地震波速度構造 / 紀伊半島 / フィリピン海スラブ |
|---|
| 研究実績の概要 |
本研究では、紀伊半島において10年間にわたり行った臨時観測等で取得した地震波形データから得た多数のレシーバ関数を用いて、フルウェーブインバージョンを実行し、紀伊半島下に沈み込むフィリピン海プレートの上面やその上方に位置する大陸モホ面の形状とそれらを含む深さ100kmまでの領域の地震波速度の3次元構造を従来の解析手法の数分の1の分解能で推定することを試みる。それにより、南海トラフ巨大地震の規模や強震動の予測の高度化とその背景にあるフィリピン海プレートの形状やプレートから放出された流体の挙動や流体とスロー地震の関係などの基礎的な研究の深化に貢献することを目的とする。
令和5年度は、有限差分法に基づく3次元粘弾性体中の地震波伝播の数値シミュレータであるOpenSWPCを用いて、フィリピン海スラブを単純化した傾斜構造モデルに対してレシーバ関数を計算することを試みた。高周波近似の波線理論に基づいて同じ傾斜構造モデルに対して計算されたレシーバ関数と比較すると、主要なPS変換波については、走時はほぼ合っているが、振幅にはかなりの差異が認められた。計算ノイズと思われる小振幅の波が目立つ結果が散見された。OpenSWPCによる計算では、モデル領域の設定、速度構造のグリッド間隔、計算する地震波のサンプリング間隔などについての調整が必要であることがわかった。
フルウェーブインバージョンの初期モデルの作成に有用となる、レシーバ関数から1次元S波速度構造を推定する遺伝アルゴリズム（GA）の改良を行った。各世代の最も良いモデルを残しつつ、ランダムな大きさで発生する「絶滅」により多様性も確保するというアイデアを実装した。Computational Geosciencesに掲載される予定である。
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
4: 遅れている
理由
データの準備においては、フルウェーブインバージョンで用いるレシーバ関数の選別作業中にデータサーバが故障し、バックアップデータをフルウェーブインバージョン用の計算機にコピーし、再度、レシーバ関数の選別をし直したため。
理論波形の計算においては、当初予定していたプログラムより小規模の計算機でも実行可能であるOpenSWPCに乗り換えたため。
|
| 今後の研究の推進方策 |
令和6年度には、OpenSWPCによるの理論計算において、モデル領域の設定、速度構造のグリッド間隔、計算する地震波のサンプリング間隔などについての調整を行い、単純化されたモデルから始めて、フルウェーブインバージョンの初期モデルについても、レシーバ関数を計算できるようにする。
加えて、改良版のGAを用いて、各観測点の到来方向ごとのレシーバ関数から1次元S波速度構造を推定し、それらをマッピングして、フルウェーブインバージョンの初期モデルを作成する。
本研究のまとめとして、得られた速度構造モデルから、フィリピン海プレートの上面や大陸モホ面等の形状を求める。また、低速度異常域の分布からプレートから放出された流体の挙動を推察し、スロー地震との関係や南海トラフ巨大地震の震源断層であるプレート境界面の状態について議論する。
|
| 次年度使用額が生じた理由 |
令和5年度は研究の進捗が遅れて、学会発表や英文校閲ができなかったため、残額が生じ、次年度に繰り越すこととした。
令和6年度には、5月に千葉市で開催される日本地球惑星科学連合と、10月に新潟市で開催される日本地震学会に参加するための投稿料・参加費・旅費を支出予定である。また、論文の英文校閲料と掲載料も予定している。
|
