"Numerical simulation of seismic wave propagation in media with viscoelastic intrusions"
Vishnevsky D.M., Lisitsa V.N., Reshetova G.V.

Seismic attenuation is a widespread property of the Earth's interior due to low consolidation and fluid saturation of the rocks. However, viscoelastic formations typically take as few as 10 to 20% of the model. Thus, the use of computationally intensive algorithms, oriented on viscoelasticity, in the entire computational domains makes the simulations inefficient. In this paper an original algorithm based on the coupling of the elastic and viscoelastic statements is presented. The computationally intensive approach is used only in subdomains embedding the viscoelastic formations, whereas the ideal elastic equations are solved elsewhere. As a result, the speed-up of the hybrid algorithm is 1.7 compared to a purely viscoelastic implementation.

Keywords: elasticity, seismic attenuation, finite-difference schemes, domain decomposition, parallel algorithms

Vishnevsky D.M., e-mail: vishnevskydm@ipgg.sbras.ru;   Lisitsa V.N., e-mail: lisitsavv@ipgg.sbras.ru – Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of Russian Academy of Sciences; prospect Koptyuga 3, Novosibirsk, 630090, Russia;
Reshetova G.V., e-mail: kgv@sscc.nsc.ru – Institute of Computational Mathematics and Mathematical Geophysics, Siberian Branch of Russian Academy of Sciences; prospekt Akademika Lavrentjeva 6, Novosibirsk, 630090, Russia