TY - JOUR
ID - 89240
TI - Markov Chain Monte Carlo Non-linear Geophysical Inversion with an Improved Proposal Distribution: Application to Geo-electrical Data
JO - Journal of the Earth and Space Physics
JA - JESPHYS
LA - en
SN - 2538-371X
AU - Tafaghod Khabaz, Zahra
AU - Ghanati, Reza
AD - Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: zahratafaghod73@gmail.com
AD - Corresponding Author, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: rghanati@ut.ac.ir
Y1 - 2023
PY - 2023
VL - 48
IS - 4
SP - 107
EP - 124
KW - Markov Chain Monte Carlo
KW - Non-linear inverse problem
KW - Perturbation models
KW - Principal Component Analysis (PCA)
KW - Proposal distribution
DO - 10.22059/jesphys.2022.339477.1007407
N2 - Geophysical inverse problems seek to provide quantitative information about geophysical characteristics of the Earthâ€™s subsurface for indirectly related data and measurements. It is generally formulated as an ill-posed non-linear optimization problem commonly solved through deterministic gradient-based approaches. Using these methods, despite fast convergence properties, may lead to local minima as well as impend accurate uncertainty analysis. On the contrary, formulating a geophysical inverse problem in a probabilistic framework and solving it by constructing the multi-dimensional posterior probability density (PPD) allow for complete sampling of the parameter space and the uncertainty quantification. The PPD is numerically characterized using Markov Chain Monte Carlo (MCMC) approaches. However, the convergence of the MCMC algorithm (i.e. sampling efficiency) toward the target stationary distribution highly depends upon the choice of the proposal distribution. In this paper, we develop an efficient proposal distribution based on perturbing the model parameters through an eigenvalue decomposition of the model covariance matrix in a principal component space. The covariance matrix is retrieved from an initial burn-in sampling, which is itself initiated using a linearized covariance estimate. The proposed strategy is first illustrated for inversion of hydrogeological parameters and then applied to synthetic and real geo-electrical data sets. The numerical experiments demonstrate that the presented proposal distribution takes advantage of the benefits from an accelerated convergence and mixing rate compared to the conventional Gaussian proposal distribution.
UR - https://jesphys.ut.ac.ir/article_89240.html
L1 - https://jesphys.ut.ac.ir/article_89240_079b4882d3c1064f5110c0d2bc9d1fdd.pdf
ER -