3D Elastic Wavefield Tomography In the Time Domain

Wavefield inversion tries to minimise the difference between real and modelled seismic data, which means that its success relies in the accuracy of the wave equation solver used in the process. Then, the wave equation solver should take elastic phenomena into account in order to compare its results with real field data. The purpose of our project was to develop a 3D wavefield tomography algorithm with an elastic wave equation solver in order to recover elastic properties. First formulations of the elastic inverse problem were published more than 20 years ago (Tarantola (1986); Mora (1987)), but only recently the computational resources needed to apply it to 3D have become available.