SS研 科学技術計算分科会 2021年度会合
富岳スペシャル3.0～新時代の防災・減災～
2022/1/20

大規模数値シミュレーションによる地震発生から地震動・地盤増幅評価までの統合的予測システムの構築とその社会実装のための取組

プロジェクト代表 堀 高峰
国立研究開発法人海洋研究開発機構

海域地震火山部門・地震津波予測研究開発センター

海域観測データを最大活用した海溝型巨大地震・津波の発生過程の現状把握＆推移予測システムの実現

観測システム開発研究G

地震津波モニタリング研究G

地震予測研究G

先行プロジェクト：ポスト京重点検討課題（2014年度〜2019年度）

地震・津波による複合災害の統合的予測システムの構築

K-computer (RIKEN)

Supercomputer Fugaku (RIKEN)

3

先行プロジェクト：ポスト京重点検討課題（2014年度〜2019年度）
地震・津波による被害の予測

index

distanceindex

attenuation relation

damagefragility curve

estimation

earthquake

scientific rationality
estimation

tusnamidisaster

将来：
地震・津波の物理過程のシミュレーションと構造物
一棟一棟の応答シミュレーションから，災害・被害を予測する．
予測精度が向上されることが期待される．
計算科学の進歩を最大限利用した計算ベースの予測は高い科学的合理性を持つ．
→都市のさまざまな変化に追随可能

FUTURE
simulation

現状：
都市を数百m四方の領域に分け，過去のデータの統計解析を基に，各領域での地震動と被害を計算する．
計算は経験ベースであり，精度は高くない．
他に代替がなく，経験ベースの予測に頼らざるを得なかった．

EARTHQUAKE DISASTER ASSESMENT

CURRENTempirical

地震・津波による複合災害の統合的予測システムの構築

先行プロジェクト：ポスト京重点検討課題（2014年度〜2019年度）

地震・津波による複合災害の統合的予測システムの構築・不確実さを考慮した多数のシナリオにもとづく強震動・津波による災害被害予測のための数値解析コンポーネントの高度化・増強

地震発生シナリオ

地殻変動＆波動伝播

3D不均質粘弾性構造での地震発生
シナリオ計算
コード開発

津波発生
2次元

地盤増幅＆建物応答

津波伝播

津波遡上粒子系

二次被害予測
サブ課題
Ｂ

構造物被害
メッシュ系

5

プロジェクト立ち上げの背景
昨年度のスーパーコンピュータ「京」を用いた調査成果○関東地域の広範囲を対象に３次元有限要素法を用いて地震動計算を実施○大規模、稠密、かつ地形を加味した条件でも安定した解析ができることを確認• ポスト「京」重点課題3において、「富岳」での活用を目指した地震災害予測のためのプログラム群を開発○地盤の塑性化を考慮するために、３次元有限要素法を用いて地震動計算を実施。

• サイエンスと計算機・計算科学双方での実績○大正関東地震を想定し、その影響が大きいと考えられる関東地域を計算範囲とした。

○周期0.4秒までの強震動の計• 世界初の大規模な地震解析をいくつか実現（サイエンス）
算を実施するため、・最小格子サイズ 28m• 世界最高峰の国際会議（SC）にて，5年連続で高い評価・水平格子数 6000×5000
・全ノード数 約170億

・時間ステップ0.01秒で200秒間• ゴードンベル賞ファイナリスト（’14, ’15, ’18），ベストポスター賞→36864ノード、14時間半で計算（’16, ’17），WACCPDベストペーパー賞（’16, ’17） 完了○このような大規模かつ高周波• 「富岳」コンピュータにより，さらに信頼性が向上すると期待成分も含み、地形を加味した条• 国の被害想定におけるニーズとの合致

件でも安定した解析が行えることが確認できた• 国の被害想定を担う内閣府（防災担当）では、長周期地震動評価のための推計手法の高度化を必要としていた。そのニーズに、重点課
題３で開発した大規模地震解析手法が合致。

関東地域の広範囲を対象に３次元有 9限要素法を用いて地震動計算を実施（「相模トラフ沿いの巨大地震による長周期地震動による被害予測の高度化」委員会資料 H30.6.12）

計算結果例：関東南部の工学的基板面上での揺れの大きさ（スナップショット速度）「富岳」でしかできない規模の地震動・地盤増幅計算開始が待たれる状況

「『富岳』成果創出加速プログラム(FY2020-2022)」の概要• 「早ければ2021年の運用開始を目指して開発が進められている『富岳』において、京の最大100倍の実効性能を（中略）最大限活用し、2020年度から試行的利用を通して、シミュレーションを中心とする計算科学とAIやデータ科学を組み合わせた新たな科学的パラダイムを構築し、早期に成果を創出する」
https://www.mext.go.jp/component/b_menu/other/__icsFi
les/afieldfile/2019/09/06/1420767_003.pdf

「『富岳』成果創出加速プログラム(FY2020-2022)」の概要

全19件（2件は計算資源配分のみ）

本課題の位置づけ・最終目標
プ ロ ジェ ク ト 代表
堀 高峰
• 「大規模数値シミュレーションによる地震発生から地震動・国立研究開発法人海洋研究開発機構
海域地震火山部門 地震津波予測研究開発セ ン タ ー セ ン タ ー長

地盤増幅評価までの統合的予測システムの構築とその社会実装」
プ ロ ジェ ク ト 副代表• 応募領域：②国民の生命・財産を守る取組の強化

市村 強
国の被害想定のための地震・ 津波災害予測における
ア プ リ ケーシ ョ ン 活用
東京大学地震研究所
・ 重点課題から 開発し て き た大規模有限要素を 用いた計算コ ード が、 国の被害想定のも と にな る 長周期地震動計算や津計算地球科学研究セン タ ー セ ン タ ー長波初期水位計算で用いら れる よ う ア プ リ ケーショ ン の改良を 国と 連携し て実施する 。

・ 建築・ 土木系企業等が、 国の想定と 同等の計算ができ る 仕組みを 構築する 。

地震に関する 災害被害予測のための大規模ア プ リ ケーシ ョ ン の開発
・ 「 富岳」 の性能を 引き 出すよ う に、 計算科学・ 計算機科学の最先端技術を 駆使し て 、 地震に関する 災害被害予測のための大規模ア プ リ ケーショ ン を 改良する 。

地震災害の定量的評価を超大規模計算で実施するための統合的予測システムを構築して、国の被害想定における地震動・地盤増幅評価の過程に実装し、政策判断や実務等に活用される基礎を築く。

実施体制・関係機関
共通基盤：超大規模高速FEMソルバー◆HPC・AI技術を駆使し、「富岳」の性能を引き出すFEMソルバ―の開発

カーネルコード開発：東大地震研（市村・藤田）

テーマ１：長周期地震動

テーマ２：地殻変動

テーマ３：浅層地盤増幅◆地震波伝播、長周期地震動手法の開発と適用◆地震発生、地殻変動手法の開発と適用◆浅層地盤における非線形地盤
増幅手法の開発と適用企業における長周期地震動
計算：
大崎総研・
小堀鐸二研・大成建設

国の長周期地震動計算への適用：
内閣府防災担当・東大地震研
災害科学

各テーマ向けコード改変

大規模計算
実計算用機能・実行環境整備

企業との共同研究

企業における非線形地盤増幅
計算：東電設計・東京ガス

プロジェクトの統括

国の被害想定への適用：
内閣府防災担当

東大地震研（市村・藤田）海洋研究開発機構（堀・縣）

国・企業・サイエンスコミュニティとの連携

計算結果のサイエンスコミュニティでの活用
東大理固体地球・科研費新学術「スロー地震」 ・地震調査研
究推進本部・防災科研・文科省委託事業「防災南海P」

国の津波被害想定への適用：内閣府防災担当

計算結果の都市被害予測への活用
理研R-CCS総合防災・防災科研代表機関 協力機関 連携機関 関係事業・機関○代表機関：

海洋研究開発機構（プロジェクト代表：堀高峰）○協力機関：

東京大学地震研究所（プロジェクト副代表：市村強）○連携機関：

内閣府（防災担当）、東京大学大学院理学系研究科、理研R-CCS
株式会社 大崎総合研究所、株式会社小堀鐸二研究所、大成建設株式会社、東京ガス株式会社、東電設計株式会社

by the combination of simple components

開発するコード群# of iterations

murfacemm
Elevation of interface
between 1st and 2nd layer

Elevation of interface
between 2nd and 3rd layer# ofsamples

Training data

Test data

Check data Total

22360

1000

1000

880
(88.0%)

891(89.1%)

E-Wave FEM

E-Cycle FEM

Case A19929(89.1%)

24360

wavelength of the obs
elevation1(x+50 i, y+50
elevation3(x+50 i, y+50
used as input data.
• To enhance sensitivity
thickness of layers andvalue and divided by sSTRIKE
21700(89.1%)

Summary：Case B shows
lop a practical approach to enable uncertainty analysis
Case B
22219
992
994
24205
Case B in following parts.
(99.4%)(99.2%)(99.4%)(99.4%)
accuracy was saturated amotion analysis
断層運動による地震波動場
非線形地盤増幅計算用の粘弾性変形計算のためのThe number of correct answers and correct ratios. Judged correct iffound to be improved to s
|target value - estimated value| < 2
昨年度のスーパーコンピュータ「京」をいた調査成果
用を高分解能で計算するため
大規模非線形有限要素法
and enhancing the sensit
scale analysis 大規模有限要素法
by heroic computing using full1000 points are picked up from 23,460 points as test data and check data.
○関東地域の広範囲を対象に３次元有限要素法を用いて地震動計算を実施の大規模有限要素法
tem○大規模、稠密、かつ地形を加味した条件でも安定した解析ができることを確認
By Supercomputing
By○地盤の塑性化を考慮するため
Comparison of
ed data to construct an AI methodologyに、３次元有限要素法を用いて地震動計算を実施。
earthquake
hazardogy to conduct Monte Carlo simulations○大正関東地震を想定し、そのclassification影響が大きいと考えられる関東地域を計算範囲とした。

HPC enables us to step towards high-quality
computing○周期0.4秒までの強震動の計算を実施するため、・最小格子サイズ 28m
・水平格子数 6000×5000
・全ノード数 約170億
・時間ステップ0.01秒で200秒間→36864ノード、14時間半で計算完了

itself not sufficient for HQC without considering uncertainty
hat enables heroic computing is strengthened by○このような大規模かつ高周波
hat makes uncertainty analysis possible
成分も含み、地形を加味した条Distribution by AI match
adequately tosupercomputing results
for whole domain

件でも安定した解析が行えることが確認できた

9

計算結果例：関東南部の工学的基板面上での揺れの大きさ（スナップショット速度）

地震学＆測地学への貢献

国の被害想定への貢献

地盤応答評価への貢献

（１）国の被害想定のための地震災害予測での活用• 国の被害想定を行う内閣府（防災担当）から、相模トラフの巨大地震による長周期地震動を、重点課題から協力して開発してきた大規模有限要素コードを用いて計算した結果が、初めて公表される予定• 先行プロジェクトで開発してきた統合的予測システムの数値解析コンポーネントが、国の被害想定のための基盤技術として、社会実装されることになる• これは、地形や地下構造の影響をより正確に取り込むためなど、内閣府（防災担当）が進める「長周期地震動評価のための推計手法の高度化」の第一歩であり、非線形地盤増幅評価における新たな方向性が示されることになる
これと並行して海洋研究開発機構で進めるサイト地震動評価の仕組みを用いて、本申請課題終了後に、国の被害想定で活用された計算コードと共通の地下構造モデル、さらには機構の地球シミュレータ（ES4）を、企業が実務や基盤研究で使用できるようにする。
12

被害想定に必要な精度検証・機能追加：「富岳」を用いた国の被害想定に必要な規模の地震動計算• 関東地方における長周期地震動のシミュレーション• 詳細な地形を反映させた計算モデルの構築• 588億自由度、147億個の四面体要素からなる• 「富岳」の4096計算ノードを5時間使用• 計算結果の可視化→国の長周期地震動高度化計算に本計算プログラムが使用される準備ができている「政策対応枠」で計算が実施されている

Vs1500m/s層の下端(縦5倍)

被害想定に必要な精度検証・機能追加：「富岳」を用いた国の被害想定に必要な規模の地震動計算

/2

東京大学地震研究所
計算地球科学研究セン タ ー セ ン タ ー長（２）富岳での大規模高速化の取り組み（地震研）
地震に関する 災害被害予測のための大規模ア プ リ ケーシ ョ ン の開発
・ 「 富岳」 の性能を 引き 出すよ う に、 計算科学・ 計算機科学の最先端技術を 駆使し て 、 地震に関する 災害被害予測のための大規模ア プ リ ケーショ ン を 改良する 。

• 152K-computer-node parallel scalable implicit solver for dynamic
nonlinear earthquake simulationdigital twin city（Ichimura et al., 2021）• 地下深部での地震発生から地下＆地上構造物の揺れまでのシームレスな超大規模計算を「富岳」で実施• アーキテクチャ依存が低く、他の問題にも適用可能性の高いソルバーの提案• 「富岳」全系での超大規模並列有限要素計算で高い並列性能を実現• Fast scalable implicit solver with convergence of physics-based
simulation & data-driven learning: toward high-fidelity simulation
with digital twin city
• SC20 research posterにaccept（Ichimura, Fujita et al.）

• A Fast Scalable Iterative Implicit Solver with Green's function-based
Neural Networks• SC20 the 11th Workshop on Latest Advances in Scalable Algorithms forLarge-Scale Systemsにaccept（Ichimura, Fujita et al.）2021/03/02

（２）富岳での大規模高速化の取り組み

京から富岳で何が変わったのか？

(Ichimura, et al., SC15 Gordon Bell Finalist)

地盤増幅計算→地表の建物の揺れ
両者を独立に計算「京」全系を用いた計算(Ichimura, et al., SC20 Technical Poster)

地表構造物＋地下構造物＋地盤の統合解析
変形特性の大きく異なるものを同時計算
収束しにくい高負荷大規模計算「富岳」2,048ノードでの実施

大規模ターミナル駅の地上・地下構造物を模擬した地震応答解析（「富岳」成果創出加速プログラム）

110億自由度、26億要素の計算モデルと「富岳」の98,304 CPU コア(2,048ノード)を用いて実施。
(Ichimura, Fujita, et al., SC20 Technical Poster)（２）富岳での大規模高速化の取り組み

京から富岳で何が変わったのか？

(Ichimura, et al., SC15 Gordon Bell Finalist)

地盤増幅計算→地表の建物の揺れ
両者を独立に計算「京」全系を用いた計算(Ichimura, et al., SC20 Technical Poster)

地表構造物＋地下構造物＋地盤の統合解析
変形特性の大きく異なるものを同時計算
収束しにくい高負荷大規模計算「富岳」2,048ノードでの実施

HPC Asia 2022:
The International Conference on High Performance Computing
in Asia-Pacific Region
Presented toT. Ichimura and K. Fujita (The University of Tokyo, RIKEN); K. Koyama (Fujitsu Ltd.); R. Kusakabe andY. Kikuchi (The University of Tokyo); T. Hori and M. Hori (Japan Agency for Marine-Earth Science andTechnology); L. Maddegedara (The University of Tokyo); N. Ohi, T. Nishiki and H. Inoue (Fujitsu Ltd.);
and K. Minami, S. Nishizawa, M. Tsuji and N. Ueda (RIKEN)

F or their paper entitled152K-computer-node Parallel Scalable Implicit Solver
for Dynamic Nonlinear Earthquake Simulation
January 14th, 2022

Kengo Nakajima, General Chair

Kento Sato, Program Co-Chair

Osni Marques, Program Co-Chair

関東平野を含む広域モデルで、地下深部の地震発生から構造物の揺れまでのシームレスな計算を「富岳」で実施(Ichimura, Fujita, et al., HPCAsia 2021, Best Paper Award)関東平野を含む広域モデルで、地下深部の地震発生から構造物の揺れまでのシームレスな計算を「富岳」で実施

T. Ichimura et al.

HPC Asia2022, January 12–14, 2022, Virtual Event, Japan𝑡 = 15 s

2m

Whole domain

Underground passage
with mesh visualized

𝑡 = 20 s• 震源断層• 地震波伝播• 地盤増幅• 地下構造物を含む都市の揺れ

𝑡 = 30 s

Underground structure𝑡 = 25 s

すべてを統合した計算を実現「富岳」

Domain ofapplication
simulation
0.2 kmClose-up view of the city

Ichimura et al. (2021)
メッシュ:
0.125m to 64m 3,240億自由度
Fi gure 2: Three-di m ensi onal l ow -order unstructur ed nonl i near ni te el em ent anal ysi s of an i nl and earthquak e on a faul t19
a highl y integrated city near a term inal station with soil and crust of Kanto Pl ain in Japan. The
structure m odel com prising「富岳」全系（京コンピュータの約1,000倍の解析能力を実現）

（３）企業が国の想定と同じ計算ができる仕組み作り（開発したアプリケーションの維持・普及に関する取組）• 重要構造物（超高層ビルやライフライン）に関わる建築・土木の主要企業が、国の想定と同じ大規模有限要素コードや地下構造モデルを用いて富岳で大規模計算ができる仕組みを構築し、実務で有用となる問題設定を検討した上で、富岳上での計算を最終年度に実施。
• これと並行して海洋研究開発機構で進めるサイト地震動評価の仕組みを用いて、本申請課題終了後に、国の被害想定で活用された計算コードと共通の地下構造モデル、さらには機構の地球シミュレータ（ES4）を、企業が実務や基盤研究で使用できるようにする。
• これにより、超高層ビルを対象にした長周期地震動など、実務に必要な検討に、本申請課題で開発したアプリケーションが継続的に活用されるような仕組みを構築する。

地球シミュレータ（ES4）20

企業との連携：E-wave FEMのユーザインタ
フェース整備・チュートリアル実施• Linux環境上でE-wave FEMを比較的容易に取り扱えるような計算環境を構築• 海洋研究開発機構所有の計算サーバ（プリプロセス用）と東京大学Oakforest-PACS上にまず構築し、現在では「富岳」・ES4へも展開済み。
• ユーザ用マニュアルの整備• 本課題連携企業を対象に、昨年10－12月にチュートリアルを計5回、今年10月に１回実施済み

WEB会議システムによるチュートリアルの模様• 構築済みの計算環境を使用して複数の例題解くためのハンズオントレーニング• 企業からのフィードバックを計算環境に反映中• 現実的な不均質構造でのモデル化まで• 非線形地盤応答は地震研との共同研究を発展（東京ガス・東電設計）
整備中のユーザー用マニュアル

企業による地震動計算の活用に向けて• 昨年度は、企業のユーザに使いやすい計算環境の構築を構築した。ただし、シンプルな速度構造を扱うことを想定した設計にとどまっていた。
• 今年度は、J-SHISデータベースや全国一次地下構造モデルなどの現実的なデータを取り込みやすいように、計算環境をアップデートした• 10/5はアップデートした環境に関するチュートリアルを実施• 今年度〜来年度にかけて具体的な問題設定を行い、連携企業の方々が「富岳」での大規模計算を実施予定

チュートリアルでの例題で計算対象とした３次元不均質地下構造の有限要素モデル。
東西180km、南北150km、鉛直50kmの領域に33の層（層内で地震波速度一様）からなる現実的な３次元不均質地形・地盤・地下構造モデルを構築。震源ならびに地下構造情報は防災科研からの提供

成果の利活用を含む他のプロジェクトとの連携
基盤S「震災軽減のためのヘテロ解析による地殻イメージング手法の開発とその適用 」（H30-R4）
地殻構造推定・推移予測の新たな手法開発開発された手法を現実的な問題に適用できるよう大規模化

JAMSTSEC サイト地震動評価プロジェクト
第２階層スパコン（ES4）を活用し、国の想定と同等の計算ができるプラットフォームを企業に提供
第２階層上での社会実装
産業利用へ

uu

サイエンスでの活用

科研費学術変革
Slow to Fast 地震学に参画
ゆっくり地震から通常地震までの統一的なシミュレーション・解析の実現

uu

国難である南海
地震への適用

防災対策に資する南海トラフ地震調査研究
プロジェクト（R2-R6）
現実的な地下構造モデルでの地殻活動解析の実現

社会科学系との連携
都市丸ごとシミュレーション

R-CCS総合防災・減災研究チームとの連携「富岳」での避難・交通・経済等の大規模社会科学系
シミュレーション

JAMSTEC サイト地震動評価PJ• 第２階層スパコン（ES4）を活用し、国の想定と同等の計算ができるプラットフォームを企業等に提供• 富岳加速プログラム終了後も継続した取組として進める• 国の被害想定で活用される長周期地
震動計算コードとプリ・ポスト環境の整備（数値解析リポジトリ）• 国にオーソライズされた地下構造モデルを活用できる環境整備（四次元仮想地球）

その他のプロジェクトとの連携
防災対策に資する南海トラフ地震調査研究プロジェクト（R2-R6）
現実的な地下構造モデルを用いた巨大地震震源域での固着・すべりの現状把握と推移予測のための手法開発
Fast crustal deformation computation by GPUs

15

固着・すべり状態の時空間変化の推定

地震時すべり分布の即時決定

現状：地下構造の曖昧さやプレート境界以外、分岐断層等の海域断層および粘弾性の効果が考慮されていない

現状：地下構造のモデルが単純（実際は三次元的に不均質）(a) case 1(b) case 2(c) case 3

粘弾性構造モデルでの固着・すべりの推定
00.386とその推定誤差を定量的に評価する手法を開発

三次元不均質地下構造モデルを国土地理院と用いつつ、推定不確実性の定量評
密接に連携
Figure 11. Results of one Green’s function derived from (a) the original material properties, (b) one case of the 1000 cases, and (c) an FE model with a ﬂat 価も含めた地震時すべりの即時推
surface and homogeneous material properties. Each arrow represents the direction and degree (m) of slip. the colour scale is logarithmic. The area inside thedashed box is the location of the input fault slip.

(a) original case043.12(b) difference

06.59

Figure 12. Comparison of two cases with different material properties. (a) The slip distribution using original material properties (m). (b) The largest differencein 1000 cases from the original slip distribution (m).

定手法を開発地下構造モデル・現状把握結果と過去の履歴を踏まえた推移予測その他のプロジェクトとの連携
科研費学術変革（A)Slow to Fast 地震学（R2-R6）
スロー地震から普通の地震まで、地震という現象を幅広くとらえ、深く理解し予測するための研究計画

富岳でしか計算することのできないサイエンスターゲットを、学術変革「Slow to Fast地震学」との共同研究にもとづいて提案

富岳を用いた大規模ミクロ経済シミュレーションの検討

https://www.eri.u-tokyo.ac.jp/cshpc/より「マルチエージェントモデルに基づく、災害時の全国規模の大規模ミクロ経済シミュレーション」の基礎検討をすすめており、次期PJで提案したいと考えている

• 従来の経験則にもとづいて行われている被害予測を、シミュレーションにもとづく科学的に合理的な予測に置き換えるための都市丸ごとシミュレーション技術の開発が進んでいる• 物理的な被害だけでなく、マルチエージェントでの避難や経済の「富岳」での大規模シミュレーションも開発が進められている
飯塚敦教授：IES（神戸モデル）社会実装「技術研究組合」より

成果の利活用を含む他のプロジェクトとの連携基盤S「震災軽減のためのヘテロ解析による地殻イメージング手法の開発とその適用 」（H30-R4）
地殻構造推定・推移予測の新たな手法開発開発された手法を現実的な問題に適用できるよう大規模化

JAMSTSEC サイト地震動評価プロジェクト
第２階層スパコン（ES4）を活用し、国の想定と同等の計算ができるプラットフォームを企業に提供
第２階層上での社会実装
産業利用へ

uu

サイエンスでの活用

科研費学術変革
Slow to Fast 地震学に参画
ゆっくり地震から通常地震までの統一的なシミュレーション・解析の実現

uu

国難である南海
地震への適用

防災対策に資する南海トラフ地震調査研究
プロジェクト（R2-R6）
現実的な地下構造モデルでの地殻活動解析の実現

社会科学系との連携
都市丸ごとシミュレーション

R-CCS総合防災・減災研究チームとの連携「富岳」での避難・交通・経済等の大規模社会科学系
シミュレーションマルチスケール・デジタルツイン日本列島
数百億自由度以上（最大２兆自由度）の超大規模高速有限要素法解析の実現（SC14・15・18GordonBell Finalist: 市村教授＠東大地震研）→列島スケールの応力蓄積過程からM9クラスの震源過程、都市の揺れから地盤・構造物の連成までをシームレスに扱う技術的基盤。さらに、ソルバーの高速化へのAIの活用などHPC+AIの技術開発。
都市スケール（〜1km）

列島スケール（数1,000km）

都市まるごとスケール（〜10km）
560億自由度・180億要素の地殻モデルによる地表面応答解析結果(最小要素サイズ: 5m)
T = 15 s.
0 km
-7 km

200万人の避難状況をシミュレーション

池袋
上野

地震時

地震後

構造物スケール（〜100m）
新宿

M9震源域スケール（数100ｋｍ）

東京駅

新橋
渋谷

非線形構造物モデルを使って計算した328,056棟の構造物応答結果
133,609,306,335自由度・33,212,898,352要素の地盤モデルによる地表面応答解析結果 (最小要素サイズ: 1m)

次のステップ・フォワードシミュレーションからデータ同化へ
・モデルや計算の定量的な品質保証・曖昧さの定量化