A hybrid parallel algorithm for simulating seismic wave propagation in complex 3D models containing intrusions

This paper presents an original algorithm for the simulation of seismic waves in models containing geological formations with complex properties such as anisotropy, attenuation, and small-scale inhomogeneities. Each of these structures requires special treatment: either small grid steps or computationally intense models and algorithms. These formations typically take as little as 25% of the model; thus computationally expensive approaches are only used locally, while more efficient algorithms are applied elsewhere. The designed hybrid algorithm is heterogeneous; the main focus of this research is the development of an efficient parallel implementation via domain decomposition.