The statistical approach to seismic inversion considers the fact that the data is imperfect. Random noise, unwanted signals, missing data, and evanescent waves that cannot be backpropagated make a real deterministic inversion impossible. The goal of iterative inversion is to minimize the difference between the registered data vector and the synthetic data computed as the product of the propagator matrix and the earth model step by step in a Gaussian sense. To stabilize the computation, a small term has to be added:

(1) |

Nonlinear inversion goes one step further because it iterates not only for the model vector but also for the propagator . The start propagator matrix is usually the transpose of the forward operator (Claerbout, J.F., 1989). The method of steepest descent or better conjugate gradients are in widespread use to find the minimum. To guarantee convergence a reasonable start model is a necessity.

To make sure that the true earth model can be found, the assumption that the
earth behaves Gaussian has to be made. A *real* statistical approach to
seismic inversion should drop this assumption and should use non-Gaussian
methods
like the Monte Carlo method (Tarantola, A., 1990). Only if we iterate without
any a priori information (which can be wrong) can we hope to find the real
inverse operator.
Unfortunately, this makes the solution space infinite so theoretically
we will need a computer with infinite speed. In practice it should
be sufficient to
combine Gaussian iterative nonlinear inversion with the Monte Carlo
method to improve today's results.

1/13/1998