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

Authors

DOI:

https://doi.org/10.14529/jsfi240401

Keywords:

Dorovsky model, grid-characteristic method, contact conditions

Abstract

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.

References

Aki, K., Richards, P.: Quantitative Seismology, Theory and Methods. University Science Books, second edn. (2002)

Auzinger, W., Koch, O., Thalhammer, M.: Defect-based local error estimators for high-order splitting methods involving three linear operators. Numerical Algorithms 70, 61–91 (2015). https://doi.org/10.1007/s11075-014-9935-8

Biot, M.: General theory of three dimensional consolidation. J. Appl. Phy. 12, 155–169 (1941). https://doi.org/10.1007/s11075-014-9935-8

Blokhin, A., Dorovsky, V.: Mathematical modelling in the theory of multivelocity continuum. Nova Science Publishers (1995)

Chelnokov, F.: Numerical modelling of deformational dynamic processes in media with complex structure. Ph.D. thesis, Moscow Institute of Physics and Technology, Moscow, Russia (2005)

Favorskaya, A., Zhdanov, M.S., Khokhlov, N.I., Petrov, I.B.: Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid-characteristic method. Geophysical Prospecting 66(8), 1485–1502 (2018). https://doi.org/10.1111/1365-2478.12639

Golubev, V., Shevchenko, A., Petrov, I.: Taking into account fluid saturation of bottom sediments in marine seismic survey. In: Doklady Mathematics. vol. 100, pp. 488–490. Springer (2019). https://doi.org/10.1134/S1064562419050107

Golubev, V., Shevchenko, A., Petrov, I.: Simulation of seismic wave propagation in a multicomponent oil deposit model. International Journal of Applied Mechanics 12(08), 2050084 (2020). https://doi.org/10.1142/S1758825120500842

Golubev, V., Shevchenko, A., Petrov, I.: Raising convergence order of grid-characteristic schemes for 2D linear elasticity problems using operator splitting. Computer Research and Modeling 14(4), 899–910 (2022). https://doi.org/10.20537/2076-7633-2022-14-4-899-910

LeVeque, R.: Finite Volume Methods for Hyperbolic Problems. Cambridge Texts in Applied Mathematics, Cambridge University Press (2002). https://doi.org/10.1017/CBO9780511791253

Sorokin, K., Imomnazarov, K.: Numerical solving of the liner two-dimensional dynamic problem in liquid-filled porous media. Journal of Siberian Federal University. Mathematics & Physics 3(2), 256–261 (2010). https://doi.org/10.1016/j.jcp.2012.11.001

Wolfram Research, I.: Mathematica, Version 14.1, https://www.wolfram.com/mathematica, Champaign, IL, 2024

Downloads

Published

2025-02-04

How to Cite

Golubev, V. I., & Shevchenko, A. V. (2025). Simulation of Seismic Processes with High-Order Grid-Characteristic Methods. Supercomputing Frontiers and Innovations, 11(4), 4–14. https://doi.org/10.14529/jsfi240401