3D modeling of time-harmonic elastic waves in anisotropic media using a direct Helmholtz solver

Abstract

We consider the modeling of time-harmonic elastic waves in 3D anisotropic media, using a multi-component finite-difference method. Here, we are concerned with solving the resulting algebraic equations for a large number of seismic sources and on a large 3D domain in the context of modeling seismic wave propagation with applications in full waveform inversion (FWI) in mind. The key result of this paper is the development of a massively parallel structured direct Helmholtz solver accommodating anisotropy and elasticity in the low to mid-frequency range.