Reproducing dynamic earthquake rupture processes by large scale simulations
Project/Area Number |
19K04031
|
Research Category |
Grant-in-Aid for Scientific Research (C)
|
Allocation Type | Multi-year Fund |
Section | 一般 |
Review Section |
Basic Section 17040:Solid earth sciences-related
|
Research Institution | The University of Tokyo |
Principal Investigator |
安藤 亮輔 東京大学, 大学院理学系研究科(理学部), 准教授 (10455256)
|
Co-Investigator(Kenkyū-buntansha) |
伊田 明弘 国立研究開発法人海洋研究開発機構, 付加価値情報創生部門(地球情報基盤センター), 副主任研究員 (80742121)
|
Project Period (FY) |
2019-04-01 – 2024-03-31
|
Project Status |
Granted (Fiscal Year 2022)
|
Budget Amount *help |
¥4,290,000 (Direct Cost: ¥3,300,000、Indirect Cost: ¥990,000)
Fiscal Year 2023: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2022: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2021: ¥650,000 (Direct Cost: ¥500,000、Indirect Cost: ¥150,000)
Fiscal Year 2020: ¥1,170,000 (Direct Cost: ¥900,000、Indirect Cost: ¥270,000)
Fiscal Year 2019: ¥650,000 (Direct Cost: ¥500,000、Indirect Cost: ¥150,000)
|
Keywords | 境界要素法 / 動的破壊シミュレーション / 地震 / 非平面断層形状 / 階層行列法 / 高速計算 / 並列計算 / 断層 / カイコウラ地震 / 地震発生 / H行列 / 地震サイクル / 高性能計算 / 高速領域分割法 / モデル / シミュレーション / 動的破壊 / 地震シーケンス / FDP / H-行列法 / 境界積分方程式法 / 北海道胆振東部地震 / Ridgecrest earthquake / 破壊 / スーパーコンピュータ |
Outline of Research at the Start |
本研究は,地震の発生,すなわち動的破壊過程は,物理モデルを用いた超大規模シミュレーションでどこまで再現可能かという問いに答えることを目指す.そのために,数値シミュレーション手法の困難から無視されてきていた3次元断層形状を,忠実に考慮した物理モデルを構築する.3次元断層上の動的破壊過程をシミュレーションで扱うために,超高効率な数値計算を実現するFDP=H-行列法を3次元の大規模な動的破壊シミュレーションを超高効率で実行するという,世界初の試みを行う.実際に発生し,高精度な観測データが得られている大地震(2016年熊本地震など)を対象に物理モデルの構築とシミュレーションによる再現性の検証を行う.
|
Outline of Annual Research Achievements |
合成開口レーダーや稠密地震計観測網の整備など近年の観測技術の発達により高分解能で地震時地殻変動や断層の三次元形状を観測的に捉えることが可能となってきている.このような高分解能データと比較可能な現実的シミュレーションを実現するため, 本研究では,高効率な数値解析手法の開発と,系統的な観測との比較によるモデルの実証的研究を行っている. 境界積分方程式法を用いた完全動的破壊シミュレーションの高効率手法の開発として,これまで開発を進めているFDP=格子階層行列法について,波動伝播の走時の近似精度を3次元問題に対して詳細に調べ,震源-観測点間距離の階層近似に対して良い収束性が得られることを確認した.さらに,FDP法に用いるI領域での畳み込み演算についてタイムステップMに対して従来のO(M)からO(1)へと縮減できる新たなアルゴリズムを構築した.従来のCPU演算のためのコードを,GPUとベクトル加速機上でも使用可能となるようにコードを改造した. 3次元断層形状モデルの構築法に関して,平均滑り速度を用いた新たな断層傾斜角の推定手法を開発した.可観測な広域応力場に対して可観測な断層滑り速度ベクトルが断層傾斜角の関数となることを利用して,逆問題として定式化したものである.本手法を中国四川地域に存在するWenchuan-Maoxien断層帯に適用して断層形状モデルを構築した.さらに,動的破壊シミュレーションを複数の条件設定において実行することで,本断層帯でのシナリオ地震を作成するとともに,地域での地震動の分布を評価した.本成果は,国際誌に投稿中である. 2016年熊本地震の3次元断層形状を詳細化して,合成開口レーダー観測で推定されている多数の分岐断層を考慮するモデル化を行った.動的破壊シミュレーションを実施し,合成開口レーダーの詳細観測と整合的な分岐断層の滑り方向の分布が得られることを確認した.
|
Current Status of Research Progress |
Current Status of Research Progress
2: Research has progressed on the whole more than it was originally planned.
Reason
コロナ禍による渡航制限を受けて,海外学会での発表が当初の予定通りには進んでいないが,研究成果については着実に上がっている.
|
Strategy for Future Research Activity |
前年度までに行えなかった研究成果の発表を精力的に行う.
|
Report
(4 results)
Research Products
(22 results)