Physics reliable frugal uncertainty analysis for full waveform inversion

TL;DR

This paper introduces a computationally efficient method for estimating uncertainty in full waveform inversion using a small number of particles with Stein Variational Gradient Descent, enabling practical industrial-scale applications.

Contribution

The work proposes a frugal uncertainty estimation approach in FWI using SVGD with limited particles, improving efficiency while maintaining qualitative reliability.

Findings

Uncertainty maps qualitatively align with physics of wave propagation.

Method achieves reasonable uncertainty estimation at lower computational cost.

Underestimation of uncertainties suggests cautious application in downstream tasks.

Abstract

Full waveform inversion (FWI) enables us to obtain high-resolution velocity models of the subsurface. However, estimating the associated uncertainties in the process is not trivial. Commonly, uncertainty estimation is performed within the Bayesian framework through sampling algorithms to estimate the posterior distribution and identify the associated uncertainty. Nevertheless, such an approach has to deal with complex posterior structures (e.g., multimodality), high-dimensional model parameters, and large-scale datasets, which lead to high computational demands and time-consuming procedures. As a result, uncertainty analysis is rarely performed, especially at the industrial scale, and thus, it drives practitioners away from utilizing it for decision-making. This work proposes a frugal approach to estimate uncertainty in FWI through the Stein Variational Gradient Descent (SVGD) algorithm by utilizing a relatively small number of velocity model particles. We warm-start the SVGD algorithm by perturbing the optimized velocity model obtained from a deterministic FWI procedure with random field-based perturbations. Such perturbations cover the scattering (i.e., high wavenumber) and the transmission (i.e., low wavenumber) components of FWI and, thus, represent the uncertainty of the FWI holistically. We demonstrate the proposed approach on the Marmousi model; we have learned that by utilizing a relatively small number of particles, the uncertainty map presents qualitatively reliable information that honours the physics of wave propagation at a reasonable cost, allowing for the potential for industrial-scale applications. Nevertheless, given that uncertainties are underestimated, we must be careful when incorporating them into downstream tasks of seismic-driven geological and reservoir modelling.