A-optimal encoding weights for nonlinear inverse problems, with applications to the Helmholtz inverse problem

Abstract

The computational cost of solving an inverse problem governed by PDEs, using multiple experiments, increases linearly with the number of experiments. A recently proposed method to decrease this cost uses only a small number of random linear combinations of all experiments for solving the inverse problem. This approach applies to inverse problems where the PDE solution depends linearly on the right-hand side function that models the experiment. As this method is stochastic in essence, the quality of the obtained reconstructions can vary, in particular when only a small number of combinations are used. We develop a Bayesian formulation for the definition and computation of encoding weights that lead to a parameter reconstruction with the least uncertainty. We call these weights A-optimal encoding weights. Our framework applies to inverse problems where the governing PDE is nonlinear with respect to the inversion parameter field. We formulate the problem in infinite dimensions and follow the optimize-then-discretize approach, devoting special attention to the discretization and the choice of numerical methods in order to achieve a computational cost that is independent of the parameter discretization. We elaborate our method for a Helmholtz inverse problem, and derive the adjoint-based expressions for the gradient of the objective function of the optimization problem for finding the A-optimal encoding weights. The proposed method is potentially attractive for real-time monitoring applications, where one can invest the effort to compute optimal weights offline, to later solve an inverse problem repeatedly, over time, at a fraction of the initial cost.

References in corpus (2)

• A computational framework for infinite-dimensional Bayesian inverse problems. Part I: The linearized case, with application to global seismic inversion
• Simultaneous Source for non-uniform data variance and missing data

Cited by in corpus (9)

• Optimal design of acoustic metamaterial cloaks under uncertainty
• Goal-Oriented Optimal Design of Experiments for Large-Scale Bayesian Linear Inverse Problems
• Taylor approximation and variance reduction for PDE-constrained optimal control under uncertainty
• Bayesian Poroelastic Aquifer Characterization from InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate
• Derivative-Informed Projected Neural Networks for High-Dimensional Parametric Maps Governed by PDEs
• A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design
• Taylor approximation for chance constrained optimization problems governed by partial differential equations with high-dimensional random parameters
• An Optimal Experimental Design Framework for Adaptive Inflation and Covariance Localization for Ensemble Filters
• Hessian-based sampling for high-dimensional model reduction