Three dimensional finite-differential simulation of the wave field taking sphericity of the Earth into account

Authors

  • A.O. Verpakhovskaya Subbotin Institute of Geophysics of the National Academy of Sciences of Ukraine, Ukraine
  • V.N. Pilipenko Subbotin Institute of Geophysics of the National Academy of Sciences of Ukraine, Ukraine

DOI:

https://doi.org/10.24028/gzh.0203-3100.v42i6.2020.222293

Keywords:

seismology, modeling of the wave field, finite-differential method, keeping sphericity of the Earth into account

Abstract

Three dimensional modeling of the wave field makes possible to obtain a concept on distribution of seismic waves in geological medium with different complexity of its structure. The location of the source of fluctuations may be defined in this case at any depth. Therefore, three dimensional modeling of the wave field may be used while solving the problems of seismology related to the nature of strains in the Earth’s interior, while studying the dangerous phenomena of endogenous origin as well as in the case of studies of the deep structure of the Earth, the boundaries of its division and its tectonic structures.

All existing methods of numerical simulation of the wave field are based on the solving the equation of elastodynamics with the help of different methods defined by interfacial and initial conditions and are subdivided into axial-symmetrical and pseudo-spectral. The first ones are based on approximation of structural model as a rotary-symmetrical along the vertical axis of the source of fluctuations and solving of the equation of elastodynamics in sylindric or spheric coordinates and in the second ones while solving the equation of elastodynamics Fourier transformation is used.

At the Institute of Geophysics named after S.I. Subbotin NAS of Ukraine a method of finite-differential simulation of the wave field has been elaborated both in two-dimensional and three-dimensional variants. The choice of finite-differential method for proceeding of the wave field is explained by its high stability and accuracy as well as with possibility of its application to the wide class of heterogenous models of the medium. Finite-differential modeling is based on the solving of the scalar equation with application of the spacio-temporal net. While solving the problems of seismology in the case of three-dimensional modeling of the wave field we are to take sphericity of the Earth into account. For this purpose in the differential wave equation a transition from Cartesian to spherical coordinate system has been provided. The theory and algorithm of the elaborated three-dimensional modeling of the wave field with finite differential method taking into account sphericity of the Earth has been considered in the paper with demonstration on the model example.

References

Brekhovskikh, L.M., & Goncharov, V.V. (1982). Introduction to continuum mechanics (applied to the theory of waves). Moscow: Nauka, 336 p. (in Russian).

Verpakhovskaya, A.O., Sidorenko, G.D., Pilipenko, V.N., & Pilipenko, E.V. (2013). A finite-deference method of pre-stack migration procedures and modeling of the wave field with parallelizing of calculation on cluster. Geoinformatika, (3), 1-12 (in Russian).

Verpakhovskaya, A.O., Pilipenko, V.N., & Budkevich, V.B. (2015). 3D finite-difference migration of the field of refracted waves. Geofizicheskiy zhurnal, 37(3), 50-65. https://doi.org/10.24028/gzh.0203-3100.v37i3.2015.111102 (in Russian).

Verpakhovskaya, A.O., Pilipenko, V.N., & Pilipenko, E.V. (2016). 3D finite-difference migration with parallelization of computing process. Geoinformatika, (4), 29-38 (in Russian).

Godunov, S.K., & Ryabenkiy, V.S. (1977). Difference schemes. Moscow: Nauka, 440 p. (in Russian).

Pilipenko, V.N., & Verpakhovskaya, A.O. (2008). Study of the features of migration transformations of the refracted wave field using 2D and 3D finite-difference modeling of seismograms. Geofizicheskiy zhurnal, 30(1), 84-96 (in Russian).

Pilipenko, V.N., Verpakhovskaya, A.O., Pilipenko, E.V., & Sidorenko, G.D. (2013). Cluster computations for the implementation of procedures for the migration of source seismograms and wavefield modeling. Conference Proceedings, 12th EAGE International Conference on Geoinformatics - Theoretical and Applied Aspects, May 2013. https://doi.org/10.3997/2214-4609.20142419 (in Russian).

Pilipenko, V.N., & Verpakhovskaya, A.O. (2019). Three-dimensional modeling of temporal field by radial and finite-differential methods for solving the problems of seismology. Geofizicheskiy zhurnal, 41(5), 3-26. https://doi.org/10.24028/gzh.0203-3100.v41i5.2019.183636 (in Russian).

Samarskiy, A.A. (1983). The theory of difference schemes. Moscow: Nauka, 616 p. (in Russian).

Samarskiy, A.A., & Gulin, A.V. (1973). Stability of difference schemes. Moscow: Nauka, 416 p. (in Russian).

Alterman, Z.S., Aboudi, J., & Karal, F.C. (1970). Pulse propagation in a laterally heterogeneous solid elastic sphere. Geophysical Journal International, 21(3), 243-260. https://doi.org/10.1111/j.1365-246X.1970.tb01790.x.

Igel, H., & Gudmundsson, O. (1997). Frequency-dependent effects on travel times and waveforms of long-period S and SS waves. Physics of the Earth and Planetary Interiors, 104(1-3), 229-246. https://doi.org/10.1016/S0031-9201(97)00051-4.

Igel, H., & Weber, M. (1996). P-SV wave propagation in the Earth’s mantle using finite-differences: Application to heterogeneous lowermost mantle structure. Geophysical Research Letters, 23(5), 415-418. https://doi.org/10.1029/96GL00422.

Pylypenko, V.N., Verpakhovska, О.O., Starostenko, V.I., & Pavlenkova, N.I. (2011). Wave images of the crustal structure from refraction and wide-angle reflection migrations along the DOBRE profile (Dnieper-Donets paleorift). Tectonophysics, 508(1-4), 96-105. https://doi.org/10.1016/j.tecto.2010.11.009.

Takenaka, H., Furumura, T., & Fujiwara, H. (1998). Recent developments in numerical methods for ground motion simulation. In K. Irikura, K. Kudo, H. Okada, & T. Sasatani (Eds.), The Effects of Surface Geology on Seismic Motion (pp. 91-101). Rotterdam: Balkema.

Toyokuni, G., Takenaka, H., & Kanao, M. (2012). Quasi-axisymmetric finite-difference method for realistic modeling of regional and global seismic wavefield - review and application. In M. Kanao, (Ed.), Seismic Waves: Research and Analysis (pp. 85-112). IntechOpen. https://doi.org/10.5772/32422.

Verpakhovska, A., Pylypenko, V., Yegorova, T., & Murovskaya, A. (2018). Seismic image of the crust on the PANCAKE profile across the UKRAINIAN CARPATHIANS from the migration method. Journal of Geodynamics, 121, 76-87. https://doi.org/10.1016/j.jog.2018.07.006.

Wang, Y., & Takenaka, H. (2010). A scheme to treat the singularity in global seismic wavefield simulation using pseudospectral method with staggered grids. Earthquake Science, 23(2), 121-127. https://doi.org/10.1007/s11589-010-0001-x.

Published

2020-12-24

How to Cite

Verpakhovskaya, A., & Pilipenko, V. (2020). Three dimensional finite-differential simulation of the wave field taking sphericity of the Earth into account. Geofizicheskiy Zhurnal, 42(6), 176–191. https://doi.org/10.24028/gzh.0203-3100.v42i6.2020.222293

Issue

Section

Articles