The Numerical Scheme of Approximate Inverse for 2D Linear Seismic Imaging

The inverse problem of seismic imaging entails reconstructing material structures, e.g. of the Earth's subsurface, from measured reflected wave fields.

In this thesis, we consider a linearized version of this problem about an a priori known wave speed.
We address the implementation of an approximate inversion scheme in two dimensions and cover the basics of the theoretical background, which is microlocal analysis.
We thoroughly test our implementation in a series of different numerical experiments, including an experiment on a field data set.