Seismic Waveform Inversion by Stochastic Optimization
Reads0
Chats0
TLDR
It is found 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, which may lead to an order of magnitude speedup for waveform inversion in practice.Abstract:
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 the corresponding stochastic methods requires an order of magnitude fewer PDE solves than a comparable deterministic method applied to the full problem, which may lead to an order of magnitude speedup for waveform inversion in practice.read more
Citations
More filters
Journal ArticleDOI
An Effective Method for Parameter Estimation with PDE Constraints with Multiple Right-Hand Sides
TL;DR: Often, parameter estimation problems of parameter-dependent PDEs involve multiple right-hand sides and the computational cost and memory requirements increase linearly with the number of right- hand sides.
Journal ArticleDOI
Improved Bounds on Sample Size for Implicit Matrix Trace Estimators
TL;DR: New sufficient bounds are proved for the Hutchinson, Gaussian and unit vector estimators, as well as a necessary bound for the Gaussian estimator, which depend more specifically on properties of matrix “A” whose information is only available through matrix-vector products.
Journal ArticleDOI
Efficient least-squares imaging with sparsity promotion and compressive sensing
Felix J. Herrmann,Xiang Li +1 more
TL;DR: This work proposes a combination of randomized dimensionality‐reduction and divide‐and‐conquer techniques that take advantage of sophisticated sparsity‐promoting solvers that work on a series of smaller sub problems each involving a small randomized subset of data.
Journal ArticleDOI
Fast waveform inversion without source‐encoding
TL;DR: It is shown that conventional optimization strategies are bound to outperform stochastic methods in the long run, and an optimization strategy that combines the benefits of both conventional and Stochastic optimization is reviewed.
Journal ArticleDOI
Robust inversion, dimensionality reduction, and randomized sampling
TL;DR: A class of inverse problems in which the forward model is the solution operator to linear ODEs or PDEs is considered, which admits several dimensionality-reduction techniques based on data averaging or sampling, which are especially useful for large-scale problems.
References
More filters
Book
Compressed sensing
TL;DR: It is possible to design n=O(Nlog(m)) nonadaptive measurements allowing reconstruction with accuracy comparable to that attainable with direct knowledge of the N most important coefficients, and a good approximation to those N important coefficients is extracted from the n measurements by solving a linear program-Basis Pursuit in signal processing.
Book
Numerical Optimization
Jorge Nocedal,Stephen J. Wright +1 more
TL;DR: Numerical Optimization presents a comprehensive and up-to-date description of the most effective methods in continuous optimization, responding to the growing interest in optimization in engineering, science, and business by focusing on the methods that are best suited to practical problems.
Book
Iterative Methods for Sparse Linear Systems
TL;DR: This chapter discusses methods related to the normal equations of linear algebra, and some of the techniques used in this chapter were derived from previous chapters of this book.
Journal ArticleDOI
A Stochastic Approximation Method
Herbert Robbins,Sutton Monro +1 more
TL;DR: In this article, a method for making successive experiments at levels x1, x2, ··· in such a way that xn will tend to θ in probability is presented.
Journal ArticleDOI
Stable signal recovery from incomplete and inaccurate measurements
TL;DR: In this paper, the authors considered the problem of recovering a vector x ∈ R^m from incomplete and contaminated observations y = Ax ∈ e + e, where e is an error term.
Related Papers (5)
An overview of full-waveform inversion in exploration geophysics
Jean Virieux,Stéphane Operto +1 more