Full-wave 3D earthquake simulation using the double-couple model and the grid-characteristic method

Bagaev R.A., Golubev V.I., Golubeva Y.A.
 pdf (600K)

One of the destroying natural processes is the initiation of the regional seismic activity. It leads to a large number of human deaths. Much effort has been made to develop precise and robust methods for the estimation of the seismic stability of buildings. One of the most common approaches is the natural frequency method. The obvious drawback of this approach is a low precision due to the model oversimplification. The other method is a detailed simulation of dynamic processes using the finite-element method. Unfortunately, the quality of simulations is not enough due to the difficulty of setting the correct free boundary condition. That is why the development of new numerical methods for seismic stability problems is a high priority nowadays.

The present work is devoted to the study of spatial dynamic processes occurring in geological medium during an earthquake. We describe a method for simulating seismic wave propagation from the hypocenter to the day surface. To describe physical processes, we use a system of partial differential equations for a linearly elastic body of the second order, which is solved numerically by a grid-characteristic method on parallelepiped meshes. The widely used geological hypocenter model, called the “double-couple” model, was incorporated into this numerical algorithm. In this case, any heterogeneities, such as geological layers with curvilinear boundaries, gas and fluid-filled cracks, fault planes, etc., may be explicitly taken into account.

In this paper, seismic waves emitted during the earthquake initiation process are numerically simulated. Two different models are used: the homogeneous half-space and the multilayered geological massif with the day surface. All of their parameters are set based on previously published scientific articles. The adequate coincidence of the simulation results is obtained. And discrepancies may be explained by differences in numerical methods used. The numerical approach described can be extended to more complex physical models of geological media.

Keywords: mathematical simulation, numerical method, grid-characteristic method, seismic waves, earthquake, seismic resistance
Citation in English: Bagaev R.A., Golubev V.I., Golubeva Y.A. Full-wave 3D earthquake simulation using the double-couple model and the grid-characteristic method // Computer Research and Modeling, 2019, vol. 11, no. 6, pp. 1061-1067
Citation in English: Bagaev R.A., Golubev V.I., Golubeva Y.A. Full-wave 3D earthquake simulation using the double-couple model and the grid-characteristic method // Computer Research and Modeling, 2019, vol. 11, no. 6, pp. 1061-1067
DOI: 10.20537/2076-7633-2019-11-6-1061-1067
URL: http://crm-en.ics.org.ru/journal/article/2864/
Creative Commons License This work is licensed under a Creative Commons Attribution-NoDerivs 3.0 Unported License.

Indexed in Scopus

Full-text version of the journal is also available on the web site of the scientific electronic library eLIBRARY.RU

The journal is included in the Russian Science Citation Index

The journal is included in the RSCI

International Interdisciplinary Conference "Mathematics. Computing. Education"

Copyright © 2009–2024 Institute of Computer Science