Nonlinear regularization techniques for seismic tomography

الأصل · EN

The effects of several nonlinear regularization techniques are discussed in the framework of 3D seismic tomography. Traditional, linear, ℓ₂ penalties are compared to so-called sparsity promoting ℓ₁ and ℓ₀ penalties, and a total variation penalty. Which of these algorithms is judged optimal depends on the specific requirements of the scientific experiment. If the correct reproduction of model amplitudes is important, classical damping towards a smooth model using an ℓ₂ norm works almost as well as minimizing the total variation but is much more efficient. If gradients (edges of anomalies) should be resolved with a minimum of distortion, we prefer ℓ₁ damping of Daubechies-4 wavelet coefficients. It has the additional advantage of yielding a noiseless reconstruction, contrary to simple ℓ₂ minimization (`Tikhonov regularization') which should be avoided. In some of our examples, the ℓ₀ method produced notable artifacts. In addition we show how nonlinear ℓ₁ methods for finding sparse models can be competitive in speed with the widely used ℓ₂ methods, certainly under noisy conditions, so that there is no need to shun ℓ₁ penalizations.