Two-dimensional sharp boundary magnetotelluric inversion using Bayesian theory

Zhou SiJie,Huang QingHua
DOI: https://doi.org/10.6038/cjg2018M0241
2018-01-01
Abstract:As the conventional magnetotelluric (MT) inversion is not sensitive to anomalous body sharp boundary, we try to develop a two-dimensional MT inversion method based on Bayesian theory, aiming at better inversion of resistivity sharp boundary. To invert for the boundary of anomalous body, we deal with the boundary position and internal resistivity as model variable, and adopt Bayesian theory to link model variables with data. Then we sample the posteriori probability density (PDD) using the Metropolis-Hastings (MH) method of Markov Chain Monte Carlo method. In the sampling process, there is no penalty function constraint, i. e., the model is completely constrained by all the information of the data. We compare the results with those from the finite-constraint and investigate the influence of different initial sampling points on the results. Taking the acceptance rate as a reference, the importance of the selection of proposal distribution function in the MH method is illustrated by model examples. When cross-correlation of model variables is weak, marginal probability distribution is available to analyze sampling result. This method can give the distribution of the model parameters. The corresponding data distribution can be obtained from the distribution of the model parameters. The inversion method is validated by comparing with the known models and the evaluation of data fitting. This method is helpful to improve the resolution of sharp boundary MT inversion.
What problem does this paper attempt to address?