《地球物理反演理论》基本概念.pptVIP

第二节基本概念若已知噪音的统计规律，并可得到一个观测数据样本集，则唯一可说的是真正的观测数据e存在于的某个区域内。如图1.12所示：但是阴影区（边界其实是模糊的）内的每个点都可能是真实的数据e。而且这个集合内的每个点都能映射到模型空间的一个点。模型空间的阴影表示所有可以接受的模型。因而需要考虑的问题变成：假如附加噪声较小，若，是否意味着？一般说来该问题的回答是否定的。令，，则大部分问题是：即使||很小，||也可能很大。也就是说反演问题本性是不稳定的。数学家则说反演问题是病态（ill-posed或ill-conditioned）的。产生病态的原因为地球物理实验本身的特性，核函数（公式1.1）一般比较宽，因而“平滑”了模型。观测数据在某种意义上可以说是模型的平滑。作为特殊的例子，考虑层速度与叠加速度的关系：第二节基本概念（1.4）式中tm是感兴趣的最大时间（或者说地震记录的最大时间），H为Heaviside阶跃函数，tj时刻的核函数为从原点到时刻tj的方波，因此V(tj)为该时间段上层速度平方加权平均值的平方根。下面给出一个实际的层速度和相应的叠加速度例子，显然V(t)比v(t)平滑很多。第二节基本概念图1.13层速度与均方根速度（ａ为层速度ｂ为均方根速度）第二节基本概念若正演映射是“平滑”，则反演映射必然是“粗化”。观测数据的微小起伏必然被放大以恢复模型。自然观测数据中的小误差能在反演模型中产生大的特征。任何用Dix公式对不精确数据反演RMS速度的人员都必须面对这个事实。另外一个考察细小误差在恢复的模型中产生巨大起伏的方法是利用公式（1.3）的解析解。由于要采用V(t)的导数，因此即使V(t)的误差比较小，也会产生较大的差值。更甚的是，在反演映射中的被用t相乘，因此随着t的变大，稳定性也将变差。这也正是使用DIX公式人员所熟知的现象。第二节基本概念许多地球物理问题的不稳定性远比RMS反演中严重。在RMS反演中，核函数对所有时间来说振幅值是常数。在大多数地球物理问题中核函数随深度（或时间）衰减，也就是说随深度增加，对地层的观测能力越差。这种核函数的衰减所引起的反演问题的不稳定性大于上面所举例子。第三节离散化和线性化地球物理模型一般用连续函数来描述，例如：地球质量与惯性动量和密度的关系：质量=惯性动量=R:地球半径；径向距离r处的密度。上面两式可用通式表示为:上式为积分，很容易用数字计算机计算。这涉及离散数学。令为m.，令，得到非常有用的公式：此时我们称理论问题被离散化。第三节离散化和线性化由于技术上的原因，现场或实验室的观测都是在有限区间上进行记录的（例如离散频率、固定带宽等），用以代替[0，∞]区间上的全部观测值。同时为了方便计算，将连续分布的物理特性p(z)用有限参数集来描述，例如用每一层具有特定密度和厚度的层状结构来描述地层，这种做法称为参数化。在反演问题的离散形式中，若观测值与模型值不成正比，则其关系称为非线性。第三节离散化和线性化例如旅行时与速度的关系：若定义模型参数，则上面问题改写为：从而具有了d=Gm形式，现在的关系便成了线性，这种变换称为线性化参数化。不是所有的非线性问题都能很容易地变成线性问题。有些问题借助于Taylor展开可以线性化，但一些高度非线性化问题则不能线性化，例如电阻率测深问题。第四节小结?????对于反演问题，怎么样描述与公式表达，主要考虑以下几个问题：1、应用哪种参数化方式——离散的还是连续的2、地球物理数据的性质是什么？观测中的误差是什么？3、问题能不能作为数学问题提出，如果能够，它是不是适定的？4、对问题有无物理约束？5、能获得什么类型的解，达到什么精度？要求得到近似解、解的范围、还是精确解？6、问题是线性的还是非线性的？??7、问题是欠定的、超定的还是适定的？????8、什么是问题的最优解法?9、解的置信界限是什么？能否用其他方法来评价？第四节小结对于反演问题来说，主要考察两个重点：非唯一性和不稳定性第四节小结一、为了控制唯一性，在地球物理反演过程中一般通过如下三个方面来研究特定反演问题的唯一性：1、通过采用不同的模来获得不同的模型，根据物理规律确定一个合理的解2、引入附