2D frequency domain finite-difference elastic-wave modeling

We present a Finite Differences implementation of the 2D elastic wave equation. Our methodology makes use of both a staggered and a rotated grid, combining the stability and low phase dispersion of the two methods. The numerical simulation is done in the frequency domain, which allows us to easily implement an efficient boundary condition (PML) and enables the simulation of many sources simultaneously. We apply a compressed sparse matrix solver (LAPACK) and a vertical/horizontal displacement formulation in order to decrease the matrix demand for the implicit method. Results show the modeled snapshots in the time and frequency domains, which can be used as input for any seismic imaging process such as migration or Full Waveform Inversion.