Abstract:Summary Full waveform inversion has emerged as the state-of-the art high resolution seismic imaging technique, both in seismology for global and regional scale imaging and in the industry for exploration purposes. While gaining in popularity, full waveform inversion, at an operational level, remains a heavy computational process involving the repeated solution of large-scale 3D wave propagation problems. For this reason it is a common practice to focus the interpretation of the results on the final estimated model. This is forgetting full waveform inversion is an ill-posed inverse problem in a high dimensional space for which the solution is intrinsically non-unique. This is the reason why being able to qualify and quantify the uncertainty attached to a model estimated by full waveform inversion is key. To this end, we propose to extend at an operational level the concepts introduced in a previous study related to the coupling between ensemble Kalman filters and full waveform inversion. These concepts had been developed for 2D frequency-domain full waveform inversion. We extend it here to the case of 3D time-domain full waveform inversion, relying on a source subsampling strategy to assimilate progressively the data within the Kalman filter. We apply our strategy to an ocean bottom cable field dataset from the North Sea to illustrate its feasibility. We explore the convergence of the filter in terms of number of elements, and extract variance and covariance information showing which part of the model are well constrained and which are not. Analyzing the variance helps to gain insight on how well the final estimated model is constrained by the whole full waveform inversion workflow. The variance maps appears as the superposition of a smooth trend related to the geometrical spreading and a high resolution trend related to reflectors. Mapping lines of the covariance (or correlation matrix) to the model space helps to gain insight on the local resolution. Through a wave propagation analysis, we are also able to relate variance peaks in the model space to variance peaks in the data space. Compared to other posterior-covariance approximation scheme, our combination between Ensemble Kalman filter and full waveform inversion is intrinsically scalable, making it a good candidate for exploiting the recent exascale high performance computing machines.

Uncertainty Estimation for Full Waveform Inversion with a Prior Information from Depth Migration

A Bayesian Approach to Estimate Uncertainty for Full-Waveform Inversion Using A Priori Information from Depth Migration

Stochastic full waveform inversion with deep generative prior for uncertainty quantification

Local uncertainty quantification for 3D time-domain full waveform inversion with ensemble Kalman filters: application to a North sea OBC dataset

Efficient 3D Bayesian Full Waveform Inversion and Analysis of Prior Hypotheses

A probabilistic full waveform inversion of surface waves

Bayesian-based Full Waveform Inversion

Annealed Stein Variational Gradient Descent for Improved Uncertainty Estimation in Full-Waveform Inversion

Physically Structured Variational Inference for Bayesian Full Waveform Inversion

Bayesian Variational Time-lapse Full-waveform Inversion

Physics reliable frugal uncertainty analysis for full waveform inversion

WISE: full-Waveform variational Inference via Subsurface Extensions

Bayesian Inversion, Uncertainty Analysis and Interrogation Using Boosting Variational Inference

Linear Waveform Tomography Inversion Using Machine Learning Algorithms

WISER: multimodal variational inference for full-waveform inversion without dimensionality reduction

Full-waveform earthquake source inversion using simulation-based inference

Full Waveform Inversion with Combined Misfit Functions and Application in Land Seismic Data

A depth-variant seismic wavelet extraction method for basis pursuit inversion with impedance trend constraint

ADMM-based full-waveform inversion for microseismic imaging

Inexact Augmented Lagrangian Method-Based Full-waveform Inversion with Randomised Singular Value Decomposition