SIPPI: A Matlab toolbox for sampling the solution to inverse problems with complex prior information Part 2-Application to crosshole GPR tomography
Research output: Contribution to journal › Journal article › Research › peer-review
Standard
SIPPI : A Matlab toolbox for sampling the solution to inverse problems with complex prior information Part 2-Application to crosshole GPR tomography. / Hansen, Thomas Mejer; Cordua, Knud Skou; Looms, Majken Caroline; Mosegaard, Klaus.
In: Computers & Geosciences, Vol. 52, 03.2013, p. 481-492.Research output: Contribution to journal › Journal article › Research › peer-review
Harvard
APA
Vancouver
Author
Bibtex
}
RIS
TY - JOUR
T1 - SIPPI
T2 - A Matlab toolbox for sampling the solution to inverse problems with complex prior information Part 2-Application to crosshole GPR tomography
AU - Hansen, Thomas Mejer
AU - Cordua, Knud Skou
AU - Looms, Majken Caroline
AU - Mosegaard, Klaus
PY - 2013/3
Y1 - 2013/3
N2 - We present an application of the SIPPI Matlab toolbox, to obtain a sample from the a posteriori probability density function for the classical tomographic inversion problem. We consider a number of different forward models, linear and non-linear, such as ray based forward models that rely on the high frequency approximation of the wave-equation and 'fat' ray based forward models relying on finite frequency theory. In order to sample the a posteriori probability density function we make use of both least squares based inversion, for linear Gaussian inverse problems, and the extended Metropolis sampler, for non-linear non-Gaussian inverse problems. To illustrate the applicability of the SIPPI toolbox to a tomographic field data set we use a cross-borehole traveltime data set from Arrens, Denmark. Both the computer code and the data are released in the public domain using open source and open data licenses. The code has been developed to facilitate inversion of 2D and 3D travel time tomographic data using a wide range of possible a priori models and choices of forward models. (C) 2012 Elsevier Ltd. All rights reserved.
AB - We present an application of the SIPPI Matlab toolbox, to obtain a sample from the a posteriori probability density function for the classical tomographic inversion problem. We consider a number of different forward models, linear and non-linear, such as ray based forward models that rely on the high frequency approximation of the wave-equation and 'fat' ray based forward models relying on finite frequency theory. In order to sample the a posteriori probability density function we make use of both least squares based inversion, for linear Gaussian inverse problems, and the extended Metropolis sampler, for non-linear non-Gaussian inverse problems. To illustrate the applicability of the SIPPI toolbox to a tomographic field data set we use a cross-borehole traveltime data set from Arrens, Denmark. Both the computer code and the data are released in the public domain using open source and open data licenses. The code has been developed to facilitate inversion of 2D and 3D travel time tomographic data using a wide range of possible a priori models and choices of forward models. (C) 2012 Elsevier Ltd. All rights reserved.
KW - Inversion
KW - Non-linear
KW - Tomography
KW - Sampling
KW - A priori
KW - A posteriori
KW - FREQUENCY TRAVEL-TIMES
KW - SENSITIVITY KERNELS
KW - FRESNEL VOLUME
KW - SIMULATION
U2 - 10.1016/j.cageo.2012.10.001
DO - 10.1016/j.cageo.2012.10.001
M3 - Journal article
VL - 52
SP - 481
EP - 492
JO - Computers & Geosciences
JF - Computers & Geosciences
SN - 0098-3004
ER -
ID: 335428514