Mathematical modelling of elastic perturbations propagating from the earthquake hypocenter

The goal of this article is the development of the research software for the mathematical modelling of elastic waves propagation originated in the earthquake hypocenter through heterogeneous media. As a source of the perturbation the geophysical focal mechanism model based on the slipping along the fault is used. For the description of dynamic behavior of media the hyperbolic system of equations of elastic media is used with the explicit allocation of contact borders of heterogeneities. It is solved using the grid-characteristic method on curvilinear structured 3D meshes. One feature of used numerical algorithm is its high scalability per core. The usage of curvilinear meshes allows describing a wide variety of geometries with high precision. Mathematical formulation of problem, development of the research software and a set of numerical experiments were done by authors. The results of modelling of propagation of seismic perturbation through geological multilayered medium and assessment of seismic resistivity of ground facility are described in this article. The estimation of software scalability per core was carried out.