大地电磁法的线性化和非线性误差估计的差异性分析.pdfVIP

第十届中国国际地球电磁学术讨论会2011 大地电磁法的线性化和非线性误差估计的差异性分析 柳建新，郭荣文‘，童孝忠，邓小康，刘文劫 有色资源与地质灾害湖南省重点实验室，中南大学地球科学与信息物理学院，湖南长沙410083． 擅要：本文针对lD大地电磁反演问题，直接对比了线性化和非线性不确定估计：为了达到这个目的，我们对两种 反演方法都采用了最新方法、技术。在反演中，由于模型的合理参数化一般情况下是未知的，因此考虑了两种常 用的参数化方法(欠参数化和超参数化)。对理论模型和大地电磁法的COPRODl实测数据进行了反演，并对比了 非线性和线性化反演结果。在有砦例子中，它们的反演结果相差比较大，对于低电导率薄层，线性化反演不敏感， 且线性化反演不能解决多峰PPD的反演问题。 关■词：贝叶斯反演，参数化，非线性，大地电磁法 Onlinearizedandnonlinear forMTinversion uncertainty LIU Jian—Xin，GUO Xiaokan，LIU Rong-Wen，TONGXiao-Zhong，DENGWenjie and South School Geomatics oflnfo-physics Engineering,CentralUniversity,Changsha,410038 Abstract：This inferenceto model paper theory uncertainties，and employsBayesian study parameterisation,parameter fortheone-dimensional inverse this nonlinearity magnetotelluric(MDproblem．In a糟considered．Inthose differencesareindicatedbetweenthenonlinearandlinearised approaches cases．important multi—modalPPDswhichcannotbeaddressedwithinalinearised approaches，including approximation method inversion，Nonlinearity,Parameterization,Magnetotelluric Keywords：Bayesian l 引言 一维(ID)大地电磁法反演问题是通过地表测量的电磁场，估计层状大地的地电模型，在地球物理界已经获 得了广泛的关注∞。大部分大地电磁反演方法是基于线性化参数估计∞。而很少有人关注非线性分析和不确定度分 析。贝叶斯推断理论为这些研究提供了天然的框架，可以有效地解决反演结果的不确定估计。 2贝叶斯理论 这节简单地描述贝叶斯反演基本理论和方法，简单介绍贝叶斯反演推断信息。m表示未知的膨维模型参数 空间，d表示JⅣ维的数据空间，在贝叶斯反演中都看成是随机变量。为了更具一般性，假设模型参数化0也是未 知随机变量。它们之间的关系通过贝叶斯框架o P(mId，p)：—P(dlm,O—)P(m,O) (1) 。P(d，口) 进行关联，其中P(mfd，回表示m对d和0的条件概率密度函数，P(dIm，研表示d对m和口的条件概率密度 函数，P(d，国和P(m，国分别表示数据和模型与0的联合先验概率密度函数。在实际应用中，如果数据d是预先 给定，参数