The problem of numerical modeling of the process of initiation of seismic activity on Arctic shelf and its impact on engineering structures is considered. To describe the dynamic behavior of a geological array, the determining systems of equations of the theory of elasticity and acoustics are used with explicit identification of all geological layers. A distinctive feature of the developed approach is the full wave calculation of the propagation of seismic waves from the source of the earthquake to the day surface. Grid-characteristic method is used for numerical calculation on hexahedral and tetrahedral computational grids.Keywords: mathematical modeling, grid-characteristic method, earthquake, seismic stability, parallel algorithms Introduction. Currently, oil and natural gas are the main sources of energy. A significant part of their deposits is located on the shelf, which causes difficulties with their exploration and development, including building complex engineering structures on the shelf. An important issue is also the transportation of raw materials from the wellhead to the site of its accumulation, for which many kilometers of pipelines are laid along the seabed. One of the significant threats is deformation and destruction of the structures in question due to seismic activity (earthquake) in the region.Today, an innovative technological solution is the usage of composite materials for complex engineering structures. At the same time, a lot of scientific collectives around the world are engaged in studying the process of destruction of composite materials. The main results of a series of international projects WWFE (World Wide Failure Exercise), entirely devoted to the development and comparison of the criteria for volumetric destruction of polymer composite materials (PCM) with the reinforcement of long fibers, are given, for example, in [1, 2]. To date, many different destruction criteria have been developed, describing in-depth experiments with varying degrees of reliability. These include the criteria of Hashin, Pak, Tsai-Hill, Tsaya-Vu and Drucker-Prager.In this paper, we consider the problem of propagation of seismic waves from the source of an earthquake in the shelf region and the destruction of an engineering structure under their influence (composite pipeline is taken as an example structure). The mathematical formulation of the problem includes a system of equations of the theory of acoustics and linear elasticity, taking into account the anisotropy of the material [3]. A calculation algorithm is based on a sequential solution of two