TY - JOUR
T1 - An efficient method to solve large linearizable inverse problems under Gaussian and separability assumptions
AU - Zunino, Andrea
AU - Mosegaard, Klaus
PY - 2019/1/1
Y1 - 2019/1/1
N2 - Inverse problems where relationships are linear arise in many fields of science and engineering and, consequently, algorithms for solving them are widespread. However, when the size of the problem increases, the computational challenge becomes huge, hence, unless some simplifying assumptions are yielded, it becomes impossible to solve the problem. Our study addresses these issues for large linear(izable) inverse problems when the uncertainties can be modeled as Gaussian and the forward relationship and covariance matrices can be expressed in terms of Kronecker products. Under these conditions, we illustrate an algorithm capable of addressing very large problems with very limited storage requirements and much faster than the traditional approach. The result is a complete characterization of the posterior distribution in a probabilistic sense in terms of mean and covariance. We extend this method also to nonlinear problems where a Gauss-Newton algorithm is employed. Applications to reflection seismology, magnetic anomaly inversion and image restoration are presented.
AB - Inverse problems where relationships are linear arise in many fields of science and engineering and, consequently, algorithms for solving them are widespread. However, when the size of the problem increases, the computational challenge becomes huge, hence, unless some simplifying assumptions are yielded, it becomes impossible to solve the problem. Our study addresses these issues for large linear(izable) inverse problems when the uncertainties can be modeled as Gaussian and the forward relationship and covariance matrices can be expressed in terms of Kronecker products. Under these conditions, we illustrate an algorithm capable of addressing very large problems with very limited storage requirements and much faster than the traditional approach. The result is a complete characterization of the posterior distribution in a probabilistic sense in terms of mean and covariance. We extend this method also to nonlinear problems where a Gauss-Newton algorithm is employed. Applications to reflection seismology, magnetic anomaly inversion and image restoration are presented.
U2 - 10.1016/j.cageo.2018.09.005
DO - 10.1016/j.cageo.2018.09.005
M3 - Journal article
SN - 0098-3004
VL - 122
SP - 77
EP - 86
JO - Computers & Geosciences
JF - Computers & Geosciences
ER -