《地球物理反演理论》数学基础知识和线性问题反演.pptVIP

每个模型参数的有两个解，对于p个参数，有2p个解，但总有，置参数投影向量，将产生两个解，可解释为最小平方解的上解包络和下解包络(即置信度界限)。若数据误差是单变的和不相关的,则可认为（观测数据数目）。将代入**式可得最大平方解的上解包络和下解包络这种方法可用于检验模型解的存在范围。四、参数估算的误差边界四、参数估算的误差边界考虑到非约束最小平方解的不稳定性，(病态)，可加入一些约束条件，例如Meju(1994)建议加入平滑约束：若反演问题最佳最小平方解为，剩余值为，寻找一个满足给定的剩余门槛值的最平滑解(用度量)。即:作业:试证明式中线形反演误差分析回顾一、观测误差通过加权得到处理二、拟和优度的定义：实质：拟和误差与观测误差相比拟三、参数解析矩阵:实质是反演能否得到完美的解四、参数结算的误差边界（1）参数协方差矩阵：观测误差是如何传递到模型参数中去的。?（2）?最小平方解的上解包络和下解包络????第四章约束线性最小平方反演主要内容一、带先验信息的反演二、带有平滑度措施的反演三、带有平滑度措施反演的几何解释一、带先验信息的反演：Dm=约束方程

D为除对角线元素外都是零的对角阵，与m作用得到包含在向量h中的m的先验值，而希望向偏量(与无偏估算相对应)

求解过程：lim=对m求偏导数，并令其为零，则有：标准方程若D为单位阵，则：约束解为：上式为约束线性反演公式，该解称为：偏置线性估算技术。它的主要有点是对于存在观测误差的超定或欠定问题，均能从无限多个（似是而非）解中获得单一解。单但其合理性则取决于约束条件（当然合理并不意味着与实际模型一致）。例：地震数据折射数据的直线约束拟合：已知表示地震折射波至，试用直线（折射波时距曲线）来拟合旅行时与炮检距的关系，并假定拟合必须经过点（注意“拟合必须经过点”是一个已知条件而非约束条件。反演结果不一定满足约束条件，但必定满足已知条件）。解：所以：Dm-h=0其中：所以-Lagrange算子二、带有平滑度措施的反演对于有限不精确观测数据集，最有效的反演解法是强制要求解是平滑的。对于多解性问题，尤其是具有很小奇异值（解不稳定）的反演问题，处理的基本思想是：有怀疑就平滑?非唯一性问题的补救方法1、公式化：A：期望模型参数随位置缓慢变化，获得物理参数上平滑的模型（显然各模型参数应属于同一类物理属性），则选择相邻两个参数模型之间差最小；D为(l=p-1)阶，h为维。B：模型参数不随位置缓慢变化：这实际上是无已知条件的偏置估算形式，但可有效的衰减解的长度，从而得到稳定的反演过程。这样原反演问题变成了约束反演问题，相当于增加了下式的平方度量约束：（解的平方度）平方度量的加入，使得多解性反演问题变成在不完整、不一致、不充足的观测数据条件下，在所有解中，根据找出最平滑的解。从数学上讲，等效于在或更广义的下求的极小值，式中H=为最大允许剩余错配值。小结为：2、求解对求导有：标准方程A，最平滑解：(D为一阶差分算子)B，D=I时有：实际上就是所共知的阻尼最小平方解，这与下面的C在数学上是一致的。C，在矩阵的本征值上附加了一个正的常数，用于改善反演条件(Marquaudt(1970)解法)。回顾非约束最小平方解：即对角矩阵用置换，这样的结果与B的情况稍有不同。而阻尼最小平方法，可解的在本征值上加一个小偏差已知G=，对应的最小平方解为：又从而从而产生约束解为：注意上面的推导并没有在向量d中添加任何零值，因此，c与a，b稍有差别，在应用上，前面的算法，a，b更有适应性，因为c要求出。但在结果分析中c更清晰。3.带有平滑度措施反演的几何解释：用下图表示和，图中等值线为常数和的l维超平面，和的约束条件