...
  • 主题
高级搜索  >
历史搜索 清空历史
首页> 外文期刊>Near surface geophysics >A geostatistical Markov chain Monte Carlo inversion algorithm for electrical resistivity tomography
【24h】

A geostatistical Markov chain Monte Carlo inversion algorithm for electrical resistivity tomography

机译：电阻率断层扫描的地统计马尔可夫链蒙特卡罗反演算法

获取原文
获取原文并翻译 | 示例
           
页面导航
  • 摘要
  • 著录项
  • 相似文献
  • 相关主题

摘要

Electrical resistivity tomography is an ill-posed and nonlinear inverse problem commonly solved through deterministic gradient-based methods. These methods guarantee a fast convergence towards the final solution, but the local linearization of the inverse operator impedes accurate uncertainty assessments. On the contrary, numerical Markov chain Monte Carlo algorithms allow for accurate uncertainty appraisals, but appropriate Markov chain Monte Carlo recipes are needed to reduce the computational effort and make these approaches suitable to be applied to field data. A key aspect of any probabilistic inversion is the definition of an appropriate prior distribution of the model parameters that can also incorporate spatial constraints to mitigate the ill conditioning of the inverse problem. Usually, Gaussian priors oversimplify the actual distribution of the model parameters that often exhibit multimodality due to the presence of multiple litho-fluid facies. In this work, we develop a novel probabilistic Markov chain Monte Carlo approach for inversion of electrical resistivity tomography data. This approach jointly estimates resistivity values, litho-fluid facies, along with the associated uncertainties from the measured apparent resistivity pseudosection. In our approach, the unknown parameters include the facies model as well as the continuous resistivity values. At each spatial location, the distribution of the resistivity value is assumed to be multimodal and non-parametric with as many modes as the number of facies. An advanced Markov chain Monte Carlo algorithm (the differential evolution Markov chain) is used to efficiently sample the posterior density in a high-dimensional parameter space. A Gaussian variogram model and a first-order Markov chain respectively account for the lateral continuity of the continuous and discrete model properties, thereby reducing the null-space of solutions. The direct sequential simulation geostatistical method allows the generation of sampled models that honour both the assumed marginal prior and spatial constraints. During the Markov chain Monte Carlo walk, we iteratively sample the facies, by moving from one mode to another, and the resistivity values, by sampling within the same mode. The proposed method is first validated by inverting the data calculated from synthetic models. Then, it is applied to field data and benchmarked against a standard local inversion algorithm. Our experiments prove that the proposed Markov chain Monte Carlo inversion retrieves reliable estimations and accurate uncertainty quantifications with a reasonable computational effort.
机译：None

著录项

相似文献

  • 外文文献
  • 中文文献
  • 专利
获取原文

客服邮箱：kefu@zhangqiaokeyan.com

京公网安备：11010802029741号 ICP备案号：京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有

  • 客服微信

  • 服务号