seismic waveform inversion by stochastic optimization地震波形反演的随机优化.pdfVIP

Hindawi Publishing Corporation International Journal of Geophysics Volume 2011, Article ID 689041, 18 pages doi:10.1155/2011/689041 Research Article Seismic Waveform Inversion by Stochastic Optimization Tristan van Leeuwen, Aleksandr Y. Aravkin, and Felix J. Herrmann Department of Earth and Ocean Sciences, University of British Columbia, Vancouver, BC, Canada V6T 1Z4 Correspondence should be addressed to Tristan van Leeuwen, tleeuwen@eos.ubc.ca Received 15 December 2010; Accepted 31 March 2011 Academic Editor: Sergey Fomel Copyright © 2011 Tristan van Leeuwen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. We explore the use of stochastic optimization methods for seismic waveform inversion. The basic principle of such methods is to randomly draw a batch of realizations of a given misfit function and goes back to the 1950s. The ultimate goal of such an approach is to dramatically reduce the computational cost involved in evaluating the misfit. Following earlier work, we introduce the stochasticity in waveform inversion problem in a rigorous way via a technique called randomized trace estimation. We then review theoretical results that underlie recent developments in the use of stochastic methods for waveform inversion. We present numerical experiments to illustrate the behavior of different types of stochastic optimization methods and investigate the sensitivity to the batch size and the noise level in the data. We find that it is possible to reproduce results that are qualitatively similar to the solution of the full problem with modest batch sizes, even on noisy data. Each iteration of