Abstract:Over the past decades, surface wave methods have been routinely employed to retrieve the physical characteristics of the first tens of meters of the subsurface, particularly the shear wave velocity profiles. Traditional methods rely on the application of the multichannel analysis of surface waves to invert the fundamental and higher modes of Rayleigh waves. However, the limitations affecting this approach, such as the 1D model assumption and the high degree of subjectivity when extracting the dispersion curve, motivate us to apply the elastic full‐waveform inversion, which, despite its higher computational cost, enables leveraging the complete information embedded in the recorded seismograms. Standard approaches solve the full‐waveform inversion using gradient‐based algorithms minimizing an error function, commonly measuring the misfit between observed and predicted waveforms. However, these deterministic approaches lack proper uncertainty quantification and are susceptible to get trapped in some local minima of the error function. An alternative lies in a probabilistic framework, but, in this case, we need to deal with the huge computational effort characterizing the Bayesian approach when applied to non‐linear problems associated with expensive forward modelling and large model spaces. In this work, we present a gradient‐based Markov chain Monte Carlo full‐waveform inversion where we accelerate the sampling of the posterior distribution by compressing data and model spaces through the discrete cosine transform. Additionally, a proposal is defined as a local, Gaussian approximation of the target density, constructed using the local Hessian and gradient information of the log posterior. We first validate our method through a synthetic test where the velocity model features lateral and vertical velocity variations. Then we invert a real dataset from the InterPACIFIC project. The obtained results prove the efficiency of our proposed algorithm, which demonstrates to be robust against cycle‐skipping issues and able to provide reasonable uncertainty evaluations with an affordable computational cost.

Iterative PDE-Constrained Optimization for Seismic Full-Waveform Inversion

Full Waveform Inversion Based on Time-Integral-damping Wavefield

Full Waveform Inversion of the Second-Order Time Integral Wavefield

Full waveform inversion using extended and simultaneous sources

Full-Wavefield Inversion: An Extreme-Scale PDE-Constrained Optimization Problem

Application of an RBF-FD solver for the Helmholtz equation to full-waveform inversion

Lift and Relax for PDE-constrained inverse problems in seismic imaging

Application of Optimal Transport and the Quadratic Wasserstein Metric to Full-Waveform Inversion

General Optimization Framework for Robust and Regularized 3D Full Waveform Inversion

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

A probabilistic full waveform inversion of surface waves

A finite-difference iterative solver of the Helmholtz equation for frequency-domain seismic wave modeling and full-waveform inversion

$ω$-FWI: Robust full-waveform inversion with Fourier-based metric

Reflection full waveform inversion with second‐order optimization using the adjoint‐state method

Multipliers waveform inversion

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

Isogeometric Multi-Resolution Full Waveform Inversion based on the Finite Cell Method

Physics-Consistent Data-driven Waveform Inversion with Adaptive Data Augmentation

Immersed boundary parametrizations for full waveform inversion

Full Waveform Inversion Using a High-Dimensional Local-Coherence Misfit Function

Multiscale Data-driven Seismic Full-waveform Inversion with Field Data Study