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
- scientific article; zbMATH DE number 1126861 (Why is no real title available?)
- scientific article; zbMATH DE number 3321507 (Why is no real title available?)
- 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
- 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
(22)- Efficient D-optimal design of experiments for infinite-dimensional Bayesian linear inverse problems
- A Fast and Scalable Computational Framework for Large-Scale High-Dimensional Bayesian Optimal Experimental Design
- Sequentially optimized projections in x-ray imaging
- Optimal Bayesian experimental design for priors of compact support with application to shock-tube experiments for combustion kinetics
- Calibration of elastoplastic constitutive model parameters from full‐field data with automatic differentiation‐based sensitivities
- Bayesian design of measurements for magnetorelaxometry imaging *
- Optimal design of simultaneous source encoding
- Optimal experimental design for infinite-dimensional Bayesian inverse problems governed by PDEs: a review
- Direct imaging for the moment tensor point sources of elastic waves
- Multimodal information gain in Bayesian design of experiments
- Fast Bayesian experimental design: Laplace-based importance sampling for the expected information gain
- Calibration experimental design considering field response and model uncertainty
- Optimal Bayesian experimental design for subsurface flow problems
- Hyper-differential sensitivity analysis for inverse problems constrained by partial differential equations
- Goal-oriented optimal design of experiments for large-scale Bayesian linear inverse problems
- Edge-promoting adaptive Bayesian experimental design for X-ray imaging
- Stability estimates for the expected utility in Bayesian optimal experimental design
- Greedy selection of optimal location of sensors for uncertainty reduction in seismic moment tensor inversion
- Optimising seismic imaging design parameters via bilevel learning
- Small-noise approximation for Bayesian optimal experimental design with nuisance uncertainty
- Recovering source location, polarization, and shape of obstacle from elastic scattering data
- Large-scale Bayesian optimal experimental design with derivative-informed projected neural network
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)