MatNERApor—A Matlab Package for Numerical Modeling of Nonlinear Response of Porous Saturated Soil Deposits to P- and SH-Waves Propagation

• Published in: Applied Sciences
• Main Authors: Artem A. Krylov, Sergey A. Kovachev, Elena A. Radiuk, Konstantin A. Roginskiy, Mikhail A. Novikov, Olga S. Samylina, Leopold I. Lobkovsky, Igor P. Semiletov
• Format: Article in Journal/Newspaper
• Language: English
• Published: MDPI AG 2022
• Subjects: site-specific analysis; geohazards; seismic microzonation; numerical modeling; nonlinear earthquake response; porous saturated soils; Technology; T; Engineering (General). Civil engineering (General); TA1-2040; Biology (General); QH301-705.5; Physics; QC1-999; Chemistry; QD1-999; Laptev Sea; Ice; laptev; permafrost
• Online Access: https://doi.org/10.3390/app12094614; https://doaj.org/article/b958b8ba518b4d1da7e8a4d00cdcdc21

The paper is devoted to the problem of numerical modeling of earthquake response of porous saturated soil deposits to seismic waves propagation. Site-specific earthquake response analysis is a necessary and important component of seismic hazard assessment. Accounting for the complex structure of porous saturated soils, i.e., the content in them, in addition to the solid matrix, pore water, gas mixture and ice, is especially important for the water areas in the zones of continuous or sparse permafrost, as well as the massive release of bubble gas from bottom sediments. The purpose of this study is to introduce an algorithm and its Matlab implementation for numerical modeling of the nonlinear response of porous saturated soil deposits to vertical P- and SH-waves propagation. The presented MatNERApor package consists of a set of Matlab scripts and functions. The package was tested and verified using the records of vertical seismic arrays of the Kik-net network. In addition, the records of local earthquakes obtained by ocean bottom seismographs in the Laptev Sea in 2019–2020 were used to demonstrate the effect of the water layer above the seabed sites on the reduction of vertical motions spectra. The results of the calculations showed good agreement with the data obtained from real seismic records, which justifies the correctness of the theoretical basis of the presented algorithm and its software implementation.