An Overview of Numerical Methods Suitable for Geophysical Imaging

Modelling methods are nowadays at the heart of any geophysical interpretation approach. These are heavily relied upon by imaging techniques in elastodynamics and electromagnetism, where they are crucial for the extraction of subsurface characteristics from ever larger and denser datasets. While highfrequency or one-way approximations are very powerful and efficient, they reach their limits when complex geological settings and solutions of full equations are required for high resolution imaging. A review of three important formulations will be carried out during this presentation: the spectral method, which is very efficient and accurate but generally restricted to simple earth structures, and often layered earth structures; the pseudo-spectral, finite-difference and finite-volume methods based on strong formulation of the partial differential equations, which are easy to implement and currently represent a good compromise between accuracy, efficiency and flexibility; and the continuous or discontinuous Galerkin finite-element methods that are based on the weak formulation, which lead to more accurate earth representations and therefore to more accurate solutions, although with higher computational costs and more complex use. The choice between these different approaches is still difficult and depends on the applications. Guidelines are given here through discussion of the requirements for imaging/ inversion.