一、正反演问题的概述.pptVIP

对于连续模型，两个连续函数d(x)和g(x)的关系可用积分方程表示 是积分函数（又称Green函数）。 实际模型的选择应考虑资料精度、分辨尺度、模型的物理意义、算法的经济指标以及折衷原则等因素。 （2）反演问题的数学适定性问题 （a）适定性，对于容许的数据d的每一集合，问题的解m存在且唯一，且连续依赖于d，则线性偏微分方程问题称为适定问题，或称其解是适定的。否则为不适定的。 不适定的原因：观测数据存在误差，计算机位数有限带来的误差等。因此应引入附加条件，使其成为适定问题（条件适定—Conditionally well posed），这时解变得不唯一。 （c）稳定性 首先考察问题的性态，设数据在定义域D上的数值是x，需要计算的数学问题为f（x），实际上只能知道数据的近似值x*，既只能计算f（x*），当x*接近x，且f（x*）也与f（x）相接近，那么问题是良态的即稳定的。否则，问题对数据的摄动很敏感，则问题是病态的、不稳定的。这时不论用什么算法，其解总是不精确的。 （c）稳定性 其次讨论算法的数值稳定性，反演算法相当于定义一个新的函数f*，使对已知的数据x产生一个近似解f*（x），当x*接近x，该算法使得f（x*）接近于f*（x），那么对于良态问题表示该算法能使计算解接近精确解—算法是数值稳定的。 （3）反演问题的解及其评价 对于反演问题的解（答），人们往往习惯于希望象解正问题一样得到解的确切数值。实际上得到的是：在下面两类信息之间作出的各种折衷，一类是人们确实希望得到的信息，另一类是从已知资料（数据）实际能得到的信息。不同折衷方案给出反问题的不同解的形式。 （a）模型参数的估计值 （b）模型参数的约束值 （c）模型参数的加权平均值 模型参数的估计值 最简单类型的解是模型参数的估计值。估算值本身不能给出有关解的质量信息，若给出解的适定性评价之后，此估算值是最有用的一类。 模型参数的约束值 给出解的约束值，来确定估值的可靠性。约束可以是确定型的，也可能是概率型的。 模型参数的加权平均值 在许多情况下，有可能辨认模型参数的组合或平均值，在某种意义上，这样比确定模型参数本身更好。 ?解的质量评价与求解是同等重要的问题。 ?解估计值对资料的拟合程度、对真实解的逼近程以及对误差的放大程度都是反演问题中解的质量的定量评价标准值。 （4）反演问题的局限性 反演问题中，始终认为正演问题是已知的。 仅限于研究那些可以用数值来表示的观测资料和问题。目前自然界中不是所有物理问题都已弄清它们的机理而能提出明确的数学模型（如地球起源、地震成因等），因此其反问题还无法研究。 地球物理反演概论 I 讲授提纲 一、正反演问题的概述 二、不适定问题及其稳定性 三、 Backus-Gilbert线性反演理论 四、解线性反问题的基本数值方法 五、直接非线性反演方法 六、地球物理中的正反演应用实例 课程特点 ?是《地球物理学》专业的必修课程； ?是学生在学习计算数学、高等物理等基础课程之后的专业基础课； ?理论性和实践性均较强。 ?了解并认识地球物理反演理论和方法； ?了解地球物理问题和数学物理问题的内在本质联系； ?把地球物理问题变为数学问题，反过来又学会把解数学问题与地球物理实际问题有机联系在一起。 课程重点 （1）地球物理反演问题的一般理论，模型空间，数据空间，零空间，解的存在性、惟一性、稳定性，适定与不适定问题。 （2）线性反演理论及方法，长度的概念，模型构制，解的非惟一性，稳定性和结果评价，参数化模型的最小长度解。 （3）广义线性反演，广义逆矩阵，广义逆的求解方法。 课程主要内容： （4）Backus-Gilbert反演理论，连续介质反演理论，BG线性评价。 （5）非线性反演方法，非线性问题的线性化，最优化法，阻尼最小二乘法，人工神经网络法，遗传算法，多尺度反演。 （6）典型地球物理反演实例。 课程主要内容： 主要参考文献 [1]《地球物理反演基本理论与应用方法》，中国地质大学“211”工程建设重点资助系列教材，姚姚主编，中国地质大学出版社，2002。 [2]《地球物理反演问题》，傅淑芳等编著，地震出版社，1998。 [3]《数字信号分析与处理》，朱光明，李庆春等编著，陕西人民教育出版社，2003。 [4]《地球物理反演理论》，王家映编著，高等教育出版社，2002。 [5]《反演理论及其应用》，徐果明编著，地震出版社，2003。 [6]《地球物理反演的理论和方法》，杨文采编著，地质出版社，1997。 [7]《地球物理数据分析—离散反演理论》，王光明等译，地质出版社，1988。 [8]Inverse Problem Theory—methods for data fitting a