摘要：针对常规大地电磁（Magnetotelluric,MT）反演方法对电阻率异常体边界不太敏感的问题,本文尝试基于贝叶斯理论开展二维大地电磁电阻率尖锐边界反演研究.在反演中,模型参数由边界位置及内部电阻率组成,通过贝叶斯理论将模型参数与数据相联系,采用Markov Chain Monte Carlo（MCMC）的Metropolis-Hastings（MH）方法对后验概率密度函数（Posteriori Probability Density,PDD）进行采样.采样过程中无罚值函数约束,完全以数据自身所包含的信息对模型进行约束,同时与有限约束进行比较,并考虑不同起始采样点对结果的影响.以接受率为参考,用模型算例说明MH方法中建议分布函数选择的重要性.当模型参数间相关性较弱时,使用边缘概率分布对采样结果进行分析.该方法能给出模型参数的分布范围,并给出该模型参数范围对应的数据范围.通过与已知模型的对比及数据拟合情况分析检验了该反演方法的有效性.该方法有助于提高大地电磁尖锐边界反演的分辨能力.

关键词：贝叶斯 尖锐边界 大地电磁反演 接受率 

单位：北京大学地球与空间科学学院; 北京100871