3.9.1 計算地震工学分野での大規模数値解析手法の開発に関する研究

断層-構造系システムとは,対象とする断層と構造物から成る地殻と構造物のモデルである.断層から生成される強震動と,その強震動に対する構造物の地震応答を計算するために使われる.開発されてきた独自のマルチスケール解析手法を改良し,大規模化・高速化を実現し,断層-構造系システムの解析を行っている.なお,大規模化・高速化の結果,従来の手法を凌駕する時間・空間分解能で,断層から伝播する地震動に対する構造物の地震応答を計算することに成功した.

断層-構造系システムの根幹である地震波動の計算では,地盤・地殻構造の幾何形状を詳細にモデル化することが重要であり,このためには有限要素法を用いる必要がある.しかし,有限要素法は差分法に比べ,計算コストが膨大となる.数理的な観点から分析し,計算コストを低減させる効率的なアルゴリズムを考案し,マルチスケール解析手法の計算コードに実装した.実装に際して並列化性能を上げることにも成功した.断層-構造系システムの大規模数値解析手法の開発では,このように基礎的な数理研究と計算科学研究にも重点が置かれている.

首都直下地震を対象として,山手線内の30万を超える構造物の地震動応答解析を行えるだけの解析技術が整いつつある.10Hzまでの精度保証可能な1000億自由度級の有限要素法モデルを用いて,断層から地表までの地震動解析,地表近傍の堆積層による地盤震動解析を行う.これらの解析技術は上記の基礎的な数理研究と計算科学研究に立脚する成果であり,ハイパフォーマンスコンピューティング分野における世界的な賞のひとつであるゴードンベル賞の最終選考論文5編に2014年2015年の2年連続で選ばれた.また,ポスト「京」重点課題アプリケーションのひとつとして選定され,ポスト「京」計算機上へ向けたチューニングが実施されている.

断層-構造系システムの具体的な対象として,大規模地下トンネルや原子力発電所といった実際の大規模構造物も挙げられる.実構造物に忠実な大自由度の解析モデルを構築し,改良されたマルチスケール解析手法を適用し,地震応答を計算している.構造物の特性を理解するためには,民間企業等の協力が必須であり,共同研究を介することで実構造物のより現実的な地震応答解析手法の構築をすすめている.地殻構造の幾何形状が地殻変動の弾性・粘弾性挙動に大きな影響を及ぼすことが指摘されていることから,構築中の技術のこれらの解析への展開も進められている.2016年には,日本列島全てを含む広領域において高詳細な地殻モデルから構築した100億自由度以上の有限要素モデルを用いた弾性・粘弾性地殻変動解析等が行われた.また,2兆自由度を超える有限要素モデル構築技術及びこれを用いた地殻変動解析技術を開発し,プレート境界の応力分布推定のための超高分解能有限要素解析が可能であることを示した.これらの成果は,ハイパフォーマンスコンピューティング分野における世界的な国際会議のひとつであるSCにおいて受賞するなど計算科学の分野においても高い評価を受けている.また,2017年には上述の山手線内の1000億自由度級の有限要素法モデルを用いた解析と人工知能を組み合わせた地震の揺れの推定高度化に関するする成果がSCにおいて受賞するなど,新たな研究の展開が進むと同時にその内容も高い評価を受けている.さらに,2018年には人工知能により高性能計算を高速化するというあらたな「人工知能と高性能計算の融合の在り方」を試みた超並列ソルバーを開発し,2018年時点で世界最速のスーパーコンピュータである米国Summitにおいて従来を凌駕する高性能を達成し,上述のゴードンベル賞の最終選考論文に三度選ばれるなど,新しい分野を開拓するとともに,継続的に高い国際的評価を受けている.