Целью данной работы является разработка программного комплекса для проведения математического моделирования распространения упругих возмущений, зарождающихся в гипоцентре землетрясения, в гетерогенных средах. В качестве источника возмущения используется геофизическая модель очага землетрясения в виде подвижки по уже существующему разлому. Для описания динамического поведения среды используется система уравнений (гиперболического типа) механики твердого тела с явным выделением контактных границ неоднородностей. Она решается с помощью сеточно-характеристического метода на криволинейных структурных трехмерных сетках. Одной из особенностей используемого вычислительного алгоритма является его высокая масштабируемость по количеству вычислителей. Использование криволинейных сеток позволяет с высокой точностью описать большой класс расчетных геометрий. Авторами была проведена математическая постановка задачи, разработан программный комплекс и выполнена серия численных расчетов. В статье приведены результаты расчета распространения сейсмического возмущения в многослойной геологической среде и оценки сейсмостойкости наземного сооружения. Проведено исследование ускорения процесса расчета при увеличении количества вычислительных узлов. 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. Голубев Василий Иванович, аспирант, Московский физико-технический институт, w.golubev@mail.ru Хохлов Николай Игоревич, инженер, Московский физико-технический институт, k_h@inbox.ru. V.I. Golubev, Moscow Institute of Physics and Technology (Dolgoprudny, Russian Federation), N.I. Khokhlov, Moscow Institute of Physics and Technology (Dolgoprudny, Russian Federation)