观测数据和模型参数估算值之方差课件.ppt

《观测数据和模型参数估算值之方差课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、地球物理反演理论地球物理反演理论课题组武汉大学测绘学院内容1，线性反演问题的最小方差解2，纯欠定问题的解法3，混定问题的解法4，先验信息在模型构置中的应用5，观测数据和模型参数估算值之方差;6，线性规划——范数解;7，范数解;1，线性反演问题的最小方差解在线性反演问题中，如果观测数据的个数多于模型参数的个数，更准确地说，在Ｍ＞Ｎ＝ｒ的情况下，最简单、最常用的反演方法是最小方差法。这里ｒ是数据方程数据核Ｇ的秩。设ｅ为观测数据ｄ与理论计算值Ｇｍ之误差向量，则方差（即目标函数）为：ｅ＝展开上式得：最小方差解必须满足：所以：例:一组观测数据为：现欲用一个平面方程：拟合之。式中：为模型参数；和是第个数据

2、对应的坐标。2，纯欠定问题的解法所谓纯欠定问题，是指在中，未知参数的个数大于观测数据的个数，且矩阵的秩的情况。换言之，在个方程中，既无相关方程，也无矛盾方程存在。从线性代数理论可知，此时有无限多个解能满足线性方程组式，且其误差均为零。这是因为，虽然观测数据提供了一些确定模型参数的信息，但其数量不足以全部确定模型参数，或未提供确定模型参数的足够充分的信息。因此，解不是惟一的，甚至有无限多能拟合观测数据的解。为了求得反演问题的一个解，我们必须从无限多个能拟合观测数据的解中，挑选出一个我们所需要的特定解。因此，解方程式时，必须加上一些在观测数据中未包含的信息．这种附加给反演问题的信息叫“先验信息”（

3、prioriinformation）。。说明：加先验信息的目的是为了补充那些为确定模型参数所不足的信息。因此，为了使反演问题的解更切合实际情况，就应本着“缺什么信息，就补充什么信息”的原则。为了更好地运用“择缺补充”的原则。几类可能的先验信息：第一类先验信息是待求地球物理参数的物理性质和其可能的数值范围。如速度，密度和电阻率等的非负性，它们都不可能小于零。而且它们的数值，根据一般物理常识可以限制在一定范围以内。第二类先验信息来自于其他已知的地质、地球物理和钻井资料。比如反演地区基底的埋深，油层的厚度，或金属矿的属性等等。第三类，某此参数比其他参数对解决地球物理问题更重要，此时可以对模型参数进行

4、加权，在一定权系数约束下求解。第四类，也是纯欠定问题解法中常应用的先验信息，假定的地球物理模型“最简单”。这里所谓最简单是指在保留实际地球物理模型基本特征不变的情况下，对地球物理模型的一种简化。解的长度，比如说解的欧几里得长度为最小的模型，应该是一种最简单的模型。设式是一纯欠定问题，此时的目标函数，在该式约束之下有极小，即根据极值理论，必须引入拉格朗日算子“λ”将条件极值问题化为无条件极值问题。目标函数应为：求上述目标函数的极小值问题可以化为求：3，混定问题的解法——马夸特（Marquardt）法当线性反演问题：呈现min（M，Ｎ）＞ｒ的情况时，称为混定问题。解混定问题的方法，通常称马夸特法，

5、或脊回归法（ridgeregression），又称为阻尼最小二乘法（damping）。从M、N和r的关系看，由于矩阵的秩r意味着在方程（2.1）式中，只有ｒ个线性无关的方程，只能确定ｒ个非零的解。因此，是超定问题，而又是欠定问题。许多地球物理线性反演问题既不完全是超定问题，也不完全是欠定问题，常常表现为一种混定形式（即混定问题）。鉴于混定问题的特殊性，它既有超定问题，也有欠定问题的性质，因此不难设想其目标函数应兼有方差项和模型长度项两项内容，即：求Ｅ相对于m的偏导数，并设其为零，简化后得：上式中：为阻尼系数或加权因子，它决定预测误差项和模型范数长度项在极小化目标函数时各自之相对重要性。4，线性

6、先验信息在模型构制中的应用在线性反演问题中，如果观测数据的个数多于模型参数的个数，更准确地说，在Ｍ＞Ｎ＝ｒ的情况下，最简单、最常用的反演方法是最小方差法。这里ｒ是数据方程数据核Ｇ的秩。在解欠定问题时，需要对模型参数强加一些先验信息，以增加（N-M）个不足的信息。然而，这并不表示，在解超定问题时，就不应该、不可能对模型参数强加任何先验信息。恰恰相反，在解欠定问题时，可对模型参数强加先验信息，在解超定问题时，也可对观测数据强加已知的先验信息。先验信息在模型构制中的应用:1.对模型参数的限制;2.对观测数据的限制;3.等式限制条件的应用;4.不等式限制条件的应用.5，观测数据和模型参数估算值之方差不

7、管是超定问题、欠定问题还是混定问题，均未涉及观测数据的统计特征。实际地球物理资料的反演中，观测数据是有误差的。有误差就要遵守一定的统计特性。在欧几里得空间解地球物理反问题时，对观测数据的统计特性有何要求，这是本次课要讨论的问题。假定每一个观测数据都是随机变量，且服从高斯分布规律，即：对M个独立观测数据来说其联合分布满足：欲使随机变量d的概率最大，则必须使：因范数就意味着观测数据（或模型参数）必须遵