Bayesian inversion of log-normal eikonal equations

DOI10.1088/1361-6420/ACC888zbMATH Open1516.35545arXiv2111.11087OpenAlexW3217430866MaRDI QIDQ6101034FDOQ6101034

Authors: Viet Ha Hoang 

Publication date: 31 May 2023

Published in: Inverse Problems (Search for Journal in Brave)

Abstract: We study the Bayesian inverse problem for inferring the log-normal slowness function of the eikonal equation given noisy observation data on its solution at a set of spatial points. We study approximation of the posterior probability measure by solving the truncated eikonal equation, which contains only a finite number of terms in the Karhunen-Loeve expansion of the slowness function, by the Fast Marching Method. The error of this approximation in the Hellinger metric is deduced in terms of the truncation level of the slowness and the grid size in the Fast Marching Method resolution. It is well known that the plain Markov Chain Monte Carlo procedure for sampling the posterior probability is highly expensive. We develop and justify the convergence of a Multilevel Markov Chain Monte Carlo method. Using the heap sort procedure in solving the forward eikonal equation by the Fast Marching Method, our Multilevel Markov Chain Monte Carlo method achieves a prescribed level of accuracy for approximating the posterior expectation of quantities of interest, requiring only an essentially optimal level of complexity. Numerical examples confirm the theoretical results.

Full work available at URL: https://arxiv.org/abs/2111.11087

Recommendations

• Analysis of a multilevel Markov chain Monte Carlo finite element method for Bayesian inversion of log-normal diffusions
• A stochastic Newton MCMC method for large-scale statistical inverse problems with application to seismic inversion
• Bayesian methods for a particular inverse problem: seismic topography
• Complexity analysis of accelerated MCMC methods for Bayesian inversion
• Sequential Monte Carlo methods for Bayesian elliptic inverse problems

zbMATH Keywords

Gaussian prioroptimal convergenceeikonal equationfast marching methodBayesian inverse problemsmultilevel Markov chain Monte Carlolognormal coefficients

Mathematics Subject Classification ID

Inverse problems for PDEs (35R30) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Boundary value problems for nonlinear first-order PDEs (35F30)

Cites Work

• Statistical and computational inverse problems.
• Title not available (Why is that?)
• Inverse problems: a Bayesian perspective
• A fast marching level set method for monotonically advancing fronts.
• Title not available (Why is that?)
• Spectral gaps for a Metropolis-Hastings algorithm in infinite dimensions
• Sparse tensor discretizations of high-dimensional parametric and stochastic PDEs
• Stochastic Galerkin discretization of the log-normal isotropic diffusion problem
• Uniqueness and error analysis for Hamilton-Jacobi equations with discontinuities
• Bayesian inverse problems for functions and applications to fluid mechanics
• Complexity analysis of accelerated MCMC methods for Bayesian inversion
• Multilevel Monte Carlo methods
• Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems
• A Hierarchical Multilevel Markov Chain Monte Carlo Algorithm with Applications to Uncertainty Quantification in Subsurface Flow
• Multilevel sequential Monte Carlo samplers
• Multilevel ensemble Kalman filtering
• Numerical analysis of an inverse problem for the eikonal equation
• Multilevel Markov chain Monte Carlo for Bayesian inversion of parabolic partial differential equations under Gaussian prior
• Improved Efficiency of a Multi-Index FEM for Computational Uncertainty Quantification
• Tikhonov regularization within ensemble Kalman inversion
• An introduction to multilevel Monte Carlo methods
• Bayesian inverse problems in measure spaces with application to Burgers and Hamilton-Jacobi equations with white noise forcing
• Analysis of a multilevel Markov chain Monte Carlo finite element method for Bayesian inversion of log-normal diffusions
• Binary recovery via phase field regularization for first-arrival traveltime tomography
• Reconciling Bayesian and perimeter regularization for binary inversion

Cited In (2)

• Title not available (Why is that?)
• Polynomial surrogates for Bayesian traveltime tomography

This page was built for publication: Bayesian inversion of log-normal eikonal equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6101034)