日本地球惑星科学連合2023年大会

講演情報

[J] 口頭発表

セッション記号 S (固体地球科学) » S-CG 固体地球科学複合領域・一般

[S-CG55] 機械学習による固体地球科学の牽引

2023年5月21日(日) 10:45 〜 12:15 302 (幕張メッセ国際会議場)

コンビーナ:久保 久彦(国立研究開発法人防災科学技術研究所)、小寺 祐貴(気象庁気象研究所)、直井 誠(京都大学)、矢野 恵佑(統計数理研究所)、座長:雨澤 勇太(国立研究開発法人産業技術総合研究所)、石塚 師也(京都大学大学院 工学研究科)

10:45 〜 11:00

[SCG55-01] 深層学習によるマルチチャンネル反射法地震探査データの速度モデル構築の高度化

*于 凡1ジャマリホンドリ エッサン2朴 進午1 (1.東京大学大気海洋研究所、2.株式会社ジオサイエンス)

キーワード:マルチチャンネル地震探査、深層学習、畳み込みニューラルネットワーク

At the Nankai Trough margin, the great megathrust earthquakes have generated strong motions and large tsunamis. The plate-boundary fault (decollement) is a source fault of the devastating Nankai megathrust earthquakes. There have been several seismic reflection studies conducted across the Nankai Trough margin to reveal the seismic structures of the Nankai Trough.
Although it is extremely important to construct a reliable velocity model for the seismic survey line and to calculate the physical properties of the subducting sediments, only a few studies on pore-fluid pressure estimation were conducted. One of the hurdles to accomplish this goal has been determining a P-wave velocity model from multi-channel seismic (MCS) reflection data, which takes a huge amount of time. Therefore, a more efficient approach is required to handle the large volume of MCS data along the Nankai Trough subduction zone. Based on these reasons, we are developing a method that can estimate the seismic velocity model, and eventually pore-fluid pressure, along the plate-boundary faults faster and more accurately by using deep learning.
In this research, we aim to incorporate convolutional neural networks (CNN) into MCS data processing and use the well-trained model to build P-wave velocity models faster. In the first step of our work, we used the Marmousi2 model (Fig. 1a) to generate the dataset for training our neural network. We first performed forward modeling using the Marmousi2 model to generate synthetic seismic dataset and then applied pre-stack depth migration (PSDM) on the synthetic data using a constant water velocity of 1500 (m/s). The common image gathers of the migrated data (with the sub-optimal water velocity) are used to calculate the vertical semblance panels for residual velocity analysis (Fig. 1b). The semblance panel (as the input) and vertical profile of the true interval velocity (as the response) for each CDP are used as pairs of training dataset. The whole dataset is randomly shuffled and split to training, validation, and test datasets. In the neural network training step, we randomly initialize the weights of the network and train the model with the shuffled dataset. After 200 training epochs, we examine the neural network with test dataset, which has not been included in the training set before. The network can predict the velocity profiles from Marmousi2 model within 100 m/s error (Fig. 1c).
Here, the neural network model was trained with synthetic data, however more diverse synthetic datasets, and eventually real data will be used to improve the model in the future. In the final step of our work, we aim to calculate the physical properties and structure of the decollement over a wide area of the Nankai Trough.