Tarantola (1987) formalizes the geophysical inverse problem by providing a theoretical approach to compensate for experimental deficiencies (e.g., acquisition geometry, complex overburden), while being consistent with the acquired data. His approach can be summarized as follows: given a linear modeling operator , compute synthetic data d using where m is a reflectivity model. Given the recorded data , a quadratic cost function,
(1) |
(2) |
The main difficulty with this approach is the explicit calculation of the inverse Hessian. In practice, it is more feasible to compute the least-squares inverse image as the solution of the linear system,
(3) |
Equation 3 states that if we convolve the Hessian matrix with a perfect model we should obtain the migration result ("Hessian impulse response"). In the next section we will study the approximations involved in the computation of the Hessian matrix, and will try to understand how far the "Hessian impulse response" computed using the approximated Hessian (from equation 3) is from the real migration result.