Fast Bayesian optimal experimental design for seismic source inversion
From MaRDI portal
Abstract: We develop a fast method for optimally designing experiments in the context of statistical seismic source inversion. In particular, we efficiently compute the optimal number and locations of the receivers or seismographs. The seismic source is modeled by a point moment tensor multiplied by a time-dependent function. The parameters include the source location, moment tensor components, and start time and frequency in the time function. The forward problem is modeled by elastodynamic wave equations. We show that the Hessian of the cost functional, which is usually defined as the square of the weighted L2 norm of the difference between the experimental data and the simulated data, is proportional to the measurement time and the number of receivers. Consequently, the posterior distribution of the parameters, in a Bayesian setting, concentrates around the "true" parameters, and we can employ Laplace approximation and speed up the estimation of the expected Kullback-Leibler divergence (expected information gain), the optimality criterion in the experimental design procedure. Since the source parameters span several magnitudes, we use a scaling matrix for efficient control of the condition number of the original Hessian matrix. We use a second-order accurate finite difference method to compute the Hessian matrix and either sparse quadrature or Monte Carlo sampling to carry out numerical integration. We demonstrate the efficiency, accuracy, and applicability of our method on a two-dimensional seismic source inversion problem.
Recommendations
- Fast Bayesian experimental design: Laplace-based importance sampling for the expected information gain
- Guided Bayesian optimal experimental design
- Efficient D-optimal design of experiments for infinite-dimensional Bayesian linear inverse problems
- Fast estimation of expected information gains for Bayesian experimental designs based on Laplace approximations
- Optimal design of simultaneous source encoding
Cites work
- A Laplace method for under-determined Bayesian optimal experimental designs
- A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data
- A stochastic collocation method for the second-order wave equation with a discontinuous random speed
- An energy absorbing far-field boundary condition for the elastic wave equation
- Bayesian experimental design: A review
- Fast estimation of expected information gains for Bayesian experimental designs based on Laplace approximations
- High dimensional polynomial interpolation on sparse grids
- scientific article; zbMATH DE number 1126861 (Why is no real title available?)
- scientific article; zbMATH DE number 3321507 (Why is no real title available?)
- Inverse problems: a Bayesian perspective
- Multilevel Monte Carlo Path Simulation
- On Information and Sufficiency
- On the optimal polynomial approximation of stochastic PDEs by Galerkin and collocation methods
- Source estimation by full wave form inversion
- Stable Difference Approximations for the Elastic Wave Equation in Second Order Formulation
- Stochastic Spectral Galerkin and Collocation Methods for PDEs with Random Coefficients: A Numerical Comparison
- Updating the Inverse of a Matrix
Cited in
(25)- Optimal Bayesian experimental design for subsurface flow problems
- Direct imaging for the moment tensor point sources of elastic waves
- Multimodal information gain in Bayesian design of experiments
- Calibration experimental design considering field response and model uncertainty
- Fast Bayesian experimental design: Laplace-based importance sampling for the expected information gain
- Small-noise approximation for Bayesian optimal experimental design with nuisance uncertainty
- Optimal design of simultaneous source encoding
- Goal-oriented optimal design of experiments for large-scale Bayesian linear inverse problems
- Efficient D-optimal design of experiments for infinite-dimensional Bayesian linear inverse problems
- Sequentially optimized projections in x-ray imaging
- Hyper-differential sensitivity analysis for inverse problems constrained by partial differential equations
- Optimal experimental design for infinite-dimensional Bayesian inverse problems governed by PDEs: a review
- Edge-promoting adaptive Bayesian experimental design for X-ray imaging
- Calibration of elastoplastic constitutive model parameters from full‐field data with automatic differentiation‐based sensitivities
- Stability estimates for the expected utility in Bayesian optimal experimental design
- Recovering source location, polarization, and shape of obstacle from elastic scattering data
- A Fast and Scalable Computational Framework for Large-Scale High-Dimensional Bayesian Optimal Experimental Design
- Bayesian design of measurements for magnetorelaxometry imaging *
- Large-scale Bayesian optimal experimental design with derivative-informed projected neural network
- Optimal Bayesian experimental design for priors of compact support with application to shock-tube experiments for combustion kinetics
- Greedy selection of optimal location of sensors for uncertainty reduction in seismic moment tensor inversion
- Optimising seismic imaging design parameters via bilevel learning
- Sequential infinite-dimensional Bayesian optimal experimental design with derivative-informed latent attention neural operator
- Virtual sensing based on hierarchical Bayesian modeling framework using a Laplace-based Gibbs sampler
- Accurate, scalable, and efficient Bayesian optimal experimental design with derivative-informed neural operators
This page was built for publication: Fast Bayesian optimal experimental design for seismic source inversion
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1734492)
