Fast Bayesian Optimal Experimental Design for Seismic Source Inversion

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.

References in corpus (1)

• Simulation-based optimal Bayesian experimental design for nonlinear systems

Cited by in corpus (18)

• Fast Bayesian experimental design: Laplace-based importance sampling for the expected information gain
• Goal-Oriented Optimal Design of Experiments for Large-Scale Bayesian Linear Inverse Problems
• Variational Bayesian experimental design for geophysical applications: seismic source location, amplitude versus offset inversion, and estimating CO2 saturations in a subsurface reservoir
• Optimal Bayesian experimental design for subsurface flow problems
• Hyper-Differential Sensitivity Analysis for Inverse Problems Constrained by Partial Differential Equations
• Calibration of Elastoplastic Constitutive Model Parameters from Full-field Data with Automatic Differentiation-based Sensitivities
• Direct imaging for the moment tensor point sources of elastic waves
• A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design
• Simplified algorithms for adaptive experiment design in parameter estimation
• A layered multiple importance sampling scheme for focused optimal Bayesian experimental design
• Recovering source location, polarization, and shape of obstacle from elastic scattering data
• Sequentially optimized projections in X-ray imaging
• hIPPYlib: An Extensible Software Framework for Large-Scale Inverse Problems Governed by PDEs; Part I: Deterministic Inversion and Linearized Bayesian Inference
• Stability estimates for the expected utility in Bayesian optimal experimental design
• Small-noise approximation for Bayesian optimal experimental design with nuisance uncertainty
• Observation Site Selection for Physical Model Parameter Estimation toward Process-Driven Seismic Wavefield Reconstruction
• Single and Multi-Objective Optimization of Distributed Acoustic Sensing Cable Layouts for Geophysical Applications
• Multimodal Information Gain in Bayesian Design of Experiments