Simulation of Seismic Processes with High-Order Grid-Characteristic Methods

This paper considers simulation of the seismic wave propagation in geological media with different rheological properties. The present work aims to construct a numerical scheme to model porous fluid-saturated medium, for the description of which the Dorovsky model was selected. We employed the grid-characteristic method, which includes choosing the appropriate operator splitting method for a 3D problem, deriving the transformation to the Riemann invariants analytically, and explicitly setting boundary and contact conditions. We simulated two scenarios. Firstly, we compared the wavefields generated by a point-source in the acoustic, linear elastic, and porous fluid-saturated approximations, noting the similarities in the longitudinal wave and differences in other wave types. Secondly, we simulated a part of the marine seismic survey process, including a source in the water layer, governed by the acoustic equations, a water-saturated layer described by the Dorovsky equations, and an explicit contact between these layers. To utilize the modern HPC multi-core and multi-processor systems, the hybrid MPI-OpenMP parallel algorithms were used.