Nonlinear regularization techniques for seismic tomography
- Department of Geosciences, Princeton University, Guyot Hall, Washington Road, Princeton, NJ 08544 (United States)
- Geosciences Azur, Universite de Nice-Sophia Antipolis, CNRS/IRD, 250 Rue Albert Einstein, Sophia Antipolis 06560 (France)
- Program in Applied and Computational Mathematics, Princeton University, Fine Hall, Washington Road, Princeton, NJ 08544 (United States)
- BP America Inc., 501 Westlake Park Blvd., Houston, TX 77079 (United States)
The effects of several nonlinear regularization techniques are discussed in the framework of 3D seismic tomography. Traditional, linear, l{sub 2} penalties are compared to so-called sparsity promoting l{sub 1} and l{sub 0} 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 l{sub 2} 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 l{sub 1} damping of Daubechies-4 wavelet coefficients. It has the additional advantage of yielding a noiseless reconstruction, contrary to simple l{sub 2} minimization ('Tikhonov regularization') which should be avoided. In some of our examples, the l{sub 0} method produced notable artifacts. In addition we show how nonlinear l{sub 1} methods for finding sparse models can be competitive in speed with the widely used l{sub 2} methods, certainly under noisy conditions, so that there is no need to shun l{sub 1} penalizations.
- OSTI ID:
- 21333927
- Journal Information:
- Journal of Computational Physics, Vol. 229, Issue 3; Other Information: DOI: 10.1016/j.jcp.2009.10.020; PII: S0021-9991(09)00572-5; Copyright (c) 2009 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Sparsity-based acoustic inversion in cross-sectional multiscale optoacoustic imaging
Image Reconstruction from Under sampled Fourier Data Using the Polynomial Annihilation Transform