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.

Markov chain Monte Carlo methods applied to the stochastic inversion of 1D viscoelastic parameters

A 3-D Numerical Study of the Random Wave-Induced Response in a Spatially Heterogenous Seabed

Markov chain Monte Carlo for petrophysical inversion

Estimation of Reservoir Fracture Properties from Seismic Data Using Markov Chain Monte Carlo Methods

Bayesian inference and Markov chain Monte Carlo based estimation of a geoscience model parameter

A Workflow to Integrate Numerical Simulation, Machine Learning Regression and Bayesian Inversion for Induced Seismicity Study: Principles and a Case Study

Markov chain Monte Carlo for seismic facies classification

Geological structure-guided hybrid MCMC and Bayesian linearized inversion methodology

Bayesian seismic inversion: a fast sampling Langevin dynamics Markov chain Monte Carlo method

Preconditioning Markov Chain Monte Carlo Method for Geomechanical Subsidence using multiscale method and machine learning technique

Bayesian Inversion of Seismic and Electromagnetic Data for Marine Gas Reservoir Characterization Using Multi-Chain Markov Chain Monte Carlo Sampling

Bayesian Variational Time-lapse Full-waveform Inversion

Assessing the Effectiveness of the Shear-Tensile-compressive Model in Earthquake Source Inversions: Synthetic Experiments and Field Application

A probabilistic full waveform inversion of surface waves

Multi-scale uncertainty quantification in geostatistical seismic inversion

Two-dimensional Time-Domain Finite-Difference Modeling for Viscoelastic Seismic Wave Propagation

A gradient‐based Markov chain Monte Carlo algorithm for elastic pre‐stack inversion with data and model space reduction

Data-driven numerical site response

Markov Chain Monte Carlo-based Bayesian method for nonlinear stochastic model updating

Evolution of the stress and strain field in the tyra field during the Post-Chalk Deposition and seismic inversion of fault zone using informed-proposal Monte Carlo

Adjoint-based inversion for stress and frictional parameters in earthquake modeling