Journal Article

Nullspace shuttles

Michael M. Deal and Guust Nolet

in Geophysical Journal International

Volume 124, issue 2, pages 372-380
Published in print February 1996 | ISSN: 0956-540X
Published online February 1996 | e-ISSN: 1365-246X | DOI:
Nullspace shuttles

Show Summary Details



In seismic tomography the problem is generally underdetermined. The solution to the tomographic problem depends on the specific optimization condition used and is inherently distorted due to noise in the data and approximations in the theory. Smoothing is often applied to reduce inversion artefacts with short correlation lengths. However, a posteriori smoothing generally affects the data fit. For more sophisticated, non-linear filters this effect can be severe. We present a technique to conserve data fit for filters of arbitrary complexity. The difference between the ‘optimal’ solution and a filtered version is projected onto the nullspace of the model space in order to preserve the data fit. Thus, we only allow changes to the image that do not conflict with the data. We demonstrate the benefits of such conservative filters using several different non-linear filters to reduce noise, smooth the image, and highlight edges.

The method is exact in small-scale experiments, where we can use the method of singular value decomposition: eigenvectors with large eigenvalues are used to project the difference between the original model and the filtered version onto the nullspace. With large-scale tomographic problems, calculation of all of the large eigenvectors is unrealistic. We show how to use the iterative method of conjugate gradients to apply conservative filters to large-scale tomographic problems with minimum computational effort.

Keywords: inversion; seismic tomography

Journal Article.  0 words. 

Subjects: Geophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.