Massively parallel hybrid algorithm for simulation of seismic wave propagation in complex media

This paper presents a problem-oriented approach, designed for numerical simulation of seismic wave propagation in models containing geological formations with complex properties such as anisotropy, attenuation, and small-scale heterogeneities. Each of the named property requires special treatment which increases computational com-plexity of the algorithm in comparison with ideally elastic isotropic media. At the same time these formations are typically relatively small taking up to 25% of the model, thus the computationally expensive approaches can be used only locally to speed-up the simulations. In this paper we discuss both mathematical and numerical aspects of hybrid algorithm paying the main attention to its parallel implementation.