Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Pre-conditioned BFGS-based uncertainty quantification in elastic full-waveform inversion

Journal Article · · Geophysical Journal International
DOI:https://doi.org/10.1093/gji/ggab375· OSTI ID:1825948

SUMMARY

Full-waveform inversion has become an essential technique for mapping geophysical subsurface structures. However, proper uncertainty quantification is often lacking in current applications. In theory, uncertainty quantification is related to the inverse Hessian (or the posterior covariance matrix). Even for common geophysical inverse problems its calculation is beyond the computational and storage capacities of the largest high-performance computing systems. In this study, we amend the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm to perform uncertainty quantification for large-scale applications. For seismic inverse problems, the limited-memory BFGS (L-BFGS) method prevails as the most efficient quasi-Newton method. We aim to augment it further to obtain an approximate inverse Hessian for uncertainty quantification in FWI. To facilitate retrieval of the inverse Hessian, we combine BFGS (essentially a full-history L-BFGS) with randomized singular value decomposition to determine a low-rank approximation of the inverse Hessian. Setting the rank number equal to the number of iterations makes this solution efficient and memory-affordable even for large-scale problems. Furthermore, based on the Gauss–Newton method, we formulate different initial, diagonal Hessian matrices as pre-conditioners for the inverse scheme and compare their performances in elastic FWI applications. We highlight our approach with the elastic Marmousi benchmark model, demonstrating the applicability of pre-conditioned BFGS for large-scale FWI and uncertainty quantification.

Sponsoring Organization:
USDOE
Grant/Contract Number:
AC05-00OR22725
OSTI ID:
1825948
Alternate ID(s):
OSTI ID: 1982495
Journal Information:
Geophysical Journal International, Journal Name: Geophysical Journal International Journal Issue: 2 Vol. 228; ISSN 0956-540X
Publisher:
Oxford University PressCopyright Statement
Country of Publication:
United Kingdom
Language:
English

References (68)

Resolution analysis by random probing: RESOLUTION ANALYSIS BY RANDOM PROBING journal August 2015
Hamiltonian Monte Carlo Inversion of Seismic Sources in Complex Media journal April 2018
On the limited memory BFGS method for large scale optimization journal August 1989
Hybrid Monte Carlo journal September 1987
Seismic Tomography and the Assessment of Uncertainty book January 2014
Matrix probing: A randomized preconditioner for the wave-equation Hessian journal March 2012
Automated Large-Scale Full Seismic Waveform Inversion for North America and the North Atlantic journal July 2018
The Collaborative Seismic Earth Model: Generation 1 journal May 2018
Square‐Root Variable Metric‐Based Nullspace Shuttle: A Characterization of the Nonuniqueness in Elastic Full‐Waveform Inversion journal February 2020
Generalized nonlinear inverse problems solved using the least squares criterion journal January 1982
Structure of the European upper mantle revealed by adjoint tomography journal June 2012
Full-waveform inversion imaging of the human brain journal March 2020
Seismic wavefield imaging of Earth’s interior across scales journal December 2019
Introduction to the spectral element method for three-dimensional seismic wave propagation journal December 1999
Improved amplitude preservation for prestack depth migration by inverse scattering theory journal September 2001
Randomized algorithms for the low-rank approximation of matrices journal December 2007
Acoustic full-waveform inversion and its uncertainty estimation based on a vector-version square-root variable metric method journal August 2019
A family of variable-metric methods derived by variational means journal January 1970
Conditioning of quasi-Newton methods for function minimization journal September 1970
A new approach to variable metric algorithms journal March 1970
A Rapidly Convergent Descent Method for Minimization journal August 1963
Probing depth and lateral variations of upper-mantle seismic anisotropy from full-waveform inversion of teleseismic body-waves journal April 2020
Global adjoint tomography—model GLAD-M25 journal May 2020
Resolution tests revisited: the power of random numbers journal December 2012
Whole-mantle radially anisotropic shear velocity structure from spectral-element waveform tomography journal September 2014
Measuring the misfit between seismograms using an optimal transport distance: application to full waveform inversion journal February 2016
Seismic waveform inversion best practices: regional, global and exploration test cases journal June 2016
Global adjoint tomography: first-generation model journal September 2016
Synthetic inversions for density using seismic and gravity data journal March 2017
Square-root variable metric based elastic full-waveform inversion—Part 2: uncertainty estimation journal May 2019
Square-root variable metric based elastic full-waveform inversion – Part 1: theory and validation journal May 2019
Source encoding for adjoint tomography journal June 2019
The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations journal January 1970
Estimation of resolution and covariance for large matrix inversions journal May 1995
Seismic tomography, adjoint methods, time reversal and banana-doughnut kernels: Seismic tomography, adjoint methods, time reversal and banana-doughnut kernels journal December 2004
A review of the adjoint-state method for computing the gradient of a functional with geophysical applications journal November 2006
Seismic tomography of the southern California crust based on spectral-element and adjoint methods journal January 2010
Hessian kernels of seismic data functionals based upon adjoint techniques: Hessian kernels journal February 2011
Forward and adjoint simulations of seismic wave propagation on fully unstructured hexahedral meshes: SPECFEM3D Version 2.0 ‘Sesame’ journal May 2011
Resolution analysis in full waveform inversion: Resolution in full waveform inversion journal October 2011
Source encoding for viscoacoustic ultrasound computed tomography journal May 2020
Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions journal January 2011
Inverse Problem Theory and Methods for Model Parameter Estimation book January 2005
Convergence Conditions for Ascent Methods journal April 1969
Full Waveform Inversion and the Truncated Newton Method journal January 2013
A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems Part I: The Linearized Case, with Application to Global Seismic Inversion journal January 2013
A Time-Domain Preconditioned Truncated Newton Approach to Visco-acoustic Multiparameter Full Waveform Inversion journal January 2018
Automatic Preconditioning by Limited Memory Quasi-Newton Updating journal January 2000
Toward a Unified Theory of Reflector Mapping journal June 1971
Reverse time migration journal November 1983
Inversion of seismic reflection data in the acoustic approximation journal August 1984
Multiscale seismic waveform inversion journal September 1995
Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model journal May 1999
Multidimensional recursive filter preconditioning in geophysical estimation problems journal March 2003
Illumination‐based normalization for wave‐equation depth migration journal July 2003
Marmousi2: An elastic upgrade for Marmousi journal February 2006
An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation journal September 2007
Seismic wave modeling for seismic imaging journal May 2009
Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion journal November 2009
An overview of full-waveform inversion in exploration geophysics journal November 2009
Elastic full-waveform inversion application using multicomponent measurements of seismic data collection journal March 2014
Adaptive waveform inversion: Theory journal November 2016
A Bayesian approach to estimate uncertainty for full-waveform inversion using a priori information from depth migration journal September 2016
Elastic least-squares reverse time migration via linearized elastic full-waveform inversion with pseudo-Hessian preconditioning journal September 2017
2D Full-Waveform Inversion and Uncertainty Estimation using the Reversible Jump Hamiltonian Monte Carlo conference August 2017
Adjustment of an Inverse Matrix Corresponding to a Change in One Element of a Given Matrix journal March 1950
Improved transparent boundary formulations for the elastic-wave equation journal December 1988
Broadband, short-period or geophone nodes? Quality assessment of Passive Seismic signals acquired during the Maupasacq experiment journal April 2018

Similar Records

Compact representations of structured BFGS matrices
Journal Article · Fri Jul 09 00:00:00 EDT 2021 · Computational Optimization and Applications · OSTI ID:1824771

Related Subjects