《(完整版)地球物理学中的反演问题》由会员分享，可在线阅读，更多相关《(完整版)地球物理学中的反演问题（27页珍藏版）》请在金锄头文库上搜索。
1、地球物理学中的反演问题1、介绍物理科学的一个重要的方面是根据数据对物理参数做出推断。通常，物理定律提 供了计算给定模型的数据值的方法，这就被称为“正演问题”，见图-1。在反演 问题中，我们的目标是根据一组测量值重建物理模型。在理想情况下，存在一个 确定的理论规定了这些数据应该怎样转换从而重现该模型。从选择的一些例子来 看，这样一个存在的理论假定了（我们）所需要的无限的、无噪声的数据是可以 获得的。在一个空间维度中，当所有能量的反射系数已知时，量子力学势能可以 被重建Marchenko,1955; Brurridge,1980。这种手法可以推广到三维空间 Newton,1989,但是在那样的情形
2、下要求有多余数据组，其中的原因并不是很 理解。在一条一维的线上的质量密度可以通过对它的所有本征频率的测量来构建 Borg,1946,但是因为这个问题的对称性，因而只有偶数部分的质量密度可以 被确定。如果（地下的）地震波速只和深度有关，那么根据地震波的距离，运用 阿贝尔变换，这个速度可以通过测定震波的抵达时间来精确构建 Herglotz,1907;Wiechert,1907。从数学上看，这个问题和构建三维空间中的 球对称量子力学势是相同的Keller et al.,1956。然而，当波速随着深度单调 增 加 时，Herglotz-Wiechert 的构 建法只 能给 出唯一解Gerver and
3、 Markushevitch,1966。这种情况和量子力学是相似的，在量子力学中，当电势 没有局部最小值时，径向对称势只能被唯一建立Sabatier,1973（量子力学 相关概念不熟悉，翻译起来有点坑）尽管精确非线性反演法在数学表达上是美妙的，但它们的适用性是有限的。 原因有很多。第一，精确的反演法通常只在理想状态下适用，这在实际中可能无 法保持。比如，Herglotz-Wiechert反演假定了地下的波速只依赖于深度并且随 着深度单调增加。地震层析成像显示这两点要求在地幔层都不满足 Nolet et a l .,1994 。第二，精确反演方法常常很不稳定。 Do r r en et a l
4、1994已经清楚 地展示了 Marchenko方程解中这种不稳定性的存在。然而，第三个原因是最根本 的。在很多反演问题中，我们要确定的模型是空间变量的一个连续函数。这意味 着该模型有无穷多的自由度。然而，在实际实验中，能够用来确定模型的数据数 量通常都是有限的。通过变量的简单计算表明这些数据不能承担足够的信息来唯 一确定模型。在线性反演问题的背景下，Backus和Gilbert1967,1968提出了 这一观点，之后Parker1994也提出来这点。这个问题对于非线性反演问题同 样相关。在实际实验中有限多的数据可以用来重建具有无穷多自由度的模型这样的 事实必然表明反演问题不是唯一的，在这个意义
5、上讲，有很多模型同样可以很好 地解释这些数据。因此，从数据反演中得到的模型不一定等于我们想要的真实模 型。这意味着图1中展示的反演问题的观点太简单了。对于现实问题，反演实际 上包含两步。用m表示真实模型，d表示数据。由数据d我们得到一个估计的模 型m，这一步称为估计问题(estimationproblem)，看图2。除了估计一个和 数据一致的模型m，我们也需要探究估计模型m和真实模型m具有什么关系。 在评价问题中，我们会确定估计模型获得了真实模型的哪些性质以及附带了哪些 误差。这部分讨论的实质就是反演=估计+评价。当我们作出一个物理解释却不承认模型中存在误差的事实以及有限的精度，这是没有多少意
6、义的 Trampert,1998。因是反演问题的非唯一性，这使得一些（通常是无穷多的）模型满足这些数据。 从技术上来讲，这个模型因为模型空间的不充分取样所以零空间存在。第二个原 因是实际数据（以及物理理论比我们想要的更频繁）总是受到误差的污染，所以 估计模型也受到这些误差的污染。所以模型评价有两个方面，非唯一性和误差传 播。模型估计和模型评价对于具有有限自由度的离散模型和具有无穷多自由度 的连续模型在根本上是不同的。而且，模型评价的问题只有在线性反演问题上得 到很好的解决。因此，离散模型和连续模型的反演是分开处理的。线性反演和非 线性反演的情况也是分开处理的。在第 2节将讨论有限数量模型参数的
7、线性反 演。在第3节中将推广为处理带有无穷多自由度的连续模型的线性反演问题。实 际上，很多反演问题都不完全是线性的，但是这些问题常常可以通过做一些适当 的近似来线性化。在第4节中将推导出单次散射近似。这种方法形成了运用于反 射地震学中的成像工具的基础。Rayleigh原理将在第5节介绍，它是关于线性 化的，构成了使用正则模态频率对地球结构进行反演的基础。地震波传播时间层 析的线性化方法是基于Fermat原理的，这将在第6节介绍。非线性反演问题要 明显难于线性反演问题。第7节将会说明非线性可能是不适定性的一个来源。目 前，对于非线性反演问题的评价问题还没有令人满意的理论。在第8节将会介绍 三种可
8、用于非线性评价问题的方法。然而，这些方法没有一个是非常令人满意的， 表明非线性反演理论是一个有重要研究挑战的领域。2、解有限的线性方程组 在前面的章节中讨论过，反演问题将有限的数据映射到一个模型上。在地球物理 学大多数实际应用中，该模型是空间坐标的一个连续函数，因此具有无穷多的自 由度。我们暂时忽略这点并假定该模型的特征可以由有限个参数确定。我们将回 到这些模型的重要情形，在第 3 节中这些模型会是无限维的。2.1 线性模型估计对于一个有限维的模型，模型参数可以规定为向量m,类似地，数据可以规定为 向量d。矩阵A通过乘积Am将数据关联到模型上。这个矩阵常常被称为理论算 子。确实，在给定的问题上
9、，它包含了我们选择给模型的所有物理和数学信息。 实际上，这些数据包含了误差e，因此记录的数据和该模型的关系应该是：d = Am + e(1)有一点需要经常注意的是，我们对于包含在模型向量m中的模型参数的选择有某 种武断性。例如，若想要描述地球的密度，我们可以选择一个模型，在该模型中， 地幔和地核具有均匀密度，在这种情况下存在两个模型参数。或者，我们可以把 大量定义在球体上的特征方程中的地球密度展开，比如描述横向变化的球谐函数 以及描述深度方向变化的多项式，这种情况会有更多的模型参数。在同一个模型 上的这两种不同参数化方法对应于不同的模型参数m和不同的矩阵A。这个例子 表明模型m不一定是真实的模
10、型，但是对模型参数的选择通常包含了对于所能构 建的模型的等级的限制。以下我们将把m认为是真模型，虽然对于它的定义存在 很多困难。由记录的数据我们得到模型的一个估计。因为这个估计实际上跟真模型是不 同的，我们用m来表示估计模型。有很多方法来设计一个逆运算将数据映射到估 计模型上e.g. Menke,1984;Tarantola,1987;Parker,1994。无论选择什么估计 量，从数据到估计模型之间最一般的线性映射可以写做算子A-g称为矩阵A的广义逆。一般来说，数据的数量不等于模型参数的数量。 因此，A通常是一个非方阵矩阵，所以它的正常逆矩阵是不存在的。随后我们将 说明广义逆矩阵a-g如何来
11、选择,但目前a-g并不需要作详细说明。被估计模型m 与真模型m之间的关系遵循如下表达式（将等式（1）代入等式（2）m = A-gAm + A-ge（3）矩阵A-gA称为精度矩阵（resolution kernel），这个算子被定义为4)R 三 A-gA表达式（3）可以写成下列形式来进行解释有限分辨率 误差部分5)在理想情况下，估计模型等于真模型向量：m = m表示我们选择的参数（列在向 量m中）可以被相互独立估计。等式（5）中最后两项分别解释了估计模型中的模糊度（blurring）和伪差（artifacts）。（A-gA-1）m描述了估计模型向量的元素是真模型向量不同元素的线性组合。我们只能取
12、得模型估计中的参数平均值 和模糊度，因为我们无法映射出最完美的细节。在理想情况下，这一项是为零的， 此时A-gA等于单位矩阵。由（4）可知，对于完美解决的模型参数，精度矩阵为 单位矩阵，即R= I（6）如前所述，通常定义向量m的模型参数的定义存在某种歧义。精度算子告诉了我 们在估计过程中我们可以独立获得的模型参数的程度。但是，精度矩阵并没有完 全告诉我们估计模型和真实的潜在物理模型之间的关系是什么，因为它没有考虑 模型参数的选择对于在估计过程中能够得到的模型的限制程度。表达式（5）中的最后一项描述了误差e是怎样映射到估计模型上去的。这些 误差并不确知，否则它们就能从数据中减去。因为数据中存在误
13、差，所以需要一 个统计分析来描述估计模型中的这些误差。当数据d不相关且有标准差b ，则jdj根据数据误差传播，模型估计m .中的标准差b表达为imic2 =X （A一gG ）2（7）m ij di j j 理想上来看，我们希望同时获得：一个完美的精度，以及不存在误差的估计模型。不幸的是，实际上这是不可能实现的。比如，使用广义逆阵A-g = 0完全 抑制了误差传播。这导致（荒谬的）估计模型m = 0，这样确实不受误差的影响。 但是，这个特殊的广义逆阵对应的精度矩阵是R = 0，显然这和理想的精度矩阵 R = I相去甚远。因此，实际上我们需要在误差传播和精度限制之间找到一个可 接受的平衡点。2.2
14、 最小二乘估计 我们现在来考虑这样的情况：独立数据的数量多于未知数的数量。在这种情况下， 等式d = Am不总是对任意给定的模型m都满足，因为数据向量中包含的可能误 差使得方程左右矛盾。例如，我们来考虑下面的问题。我们有两个物块质量分别 是m和m。第一个物块的称重得出1千克质量。某人测量第二个物块，结果得12出 2 千克质量。接下来，某人把两个物块放在一起称重，结果发现总质量是 2 千克。这个问题中测量的结果可以用下列方程组表示m = d = 111m = d = 2（8）22m + m = d = 2123相应的矩阵A表示为1 0、A= 0 1（9）J1丿显然，这个方程组是不能满足的。不可能
15、第一个物块质量是m = 1，第二个物块1的质量是m = 2，而它们的质量之和m + m = 2。显然测量中存在误差，但是没2 1 2 理由舍弃三个方程中的一个而去支持另外两个。图 3（略）生动地阐述了这个问 题。在（m ,m ）平面中，三个方程对应三条实线。三条线不相交于同一点表示这12个线性方程组存在矛盾。所以，采用合理的方法调和这些方程是确定两个物块质 量的反演问题的一部分。通常估计模型的一种方法是寻找一个能够最佳拟合数据的模型m,在这个意 义上，数据向量d和再估算数据Am之间由L-范数(即欧几里德距离，延伸阅 读泛函分析)计量的差值要尽可能小。这表示由模型给出的最小二乘解最小化下 列目标函数 1 11I11J1-S -20 目 H 3ml图-3S = |d - Am|2(10)Strang 1988详细说明了这个量由下列模型估计最小化m = CatA)i ATd(11)在图3的例子中，最小二乘解是(m ,m )平面中到三条实线距离最短的点，12这个点用一个黑色方块来表示。使用矩阵(9)，我们很容易得出问题(8)的最小二乘估计量由下式给出m = Cat A Ard = $3 J1-1 1)d2 1丿12)代入数据向量，分别得到估计模型m 二 2/31m = 5/322.3 最小范数估计在一些问题中，未知量的数量少于参
