线性反演理论3.pptVIP

线性反演理论 中国石油大学 2012年12月 概述 地球物理反演理论的目的： 根据地球物理观测数据求取相应的地球物理模型（例）。 地球物理反演问题： 1 反问题解的存在性； 2 地球模型的构制； 3 地球物理反问题解的非唯一性； 4 非唯一性解得评价。 概述 反问题解的非唯一性原因： 1 地球物理场的等效性（例） ； 2 观测数据的局限性（例） ； 3 数据测量和计算的误差（例） ； 反问题解的非唯一性解决方案： 1 扩大观测范围获取更多场的信息（高密度） ； 2 在反演过程中施加约束（测井约束） 概述 反问题解的不稳定性： 例： 1.11 x1 + x2 = 3.11 x1+ 0.90 x2 = 2.80 求解该方程组得到： x1 = 1 x2 = 2 概述 反问题解的不稳定性（续）： 例： 1.11 x1 + x2 = 3.10 x1+ 0.90 x2 = 2.81 求解该方程组得到： x1 = 20 x2 = -19.1 很小的数据误差引起解的极大波动。 概述 反问题解的不稳定性对策： — 正则化方法，使问题在一定先验条件下成为稳定的。 例： G m = dc 求解该方程组得到： m = G-1 dc 概述 反问题解的不稳定性对策： — 超定问题。 例： G m = dc 求解该方程组得到： m = ( GT G) -1 GT dc 先验条件： 观测数据dc与理论数据d方差达到极小。 概述 反问题解的不稳定性对策： — 欠定问题。 例： G m = dc 求解该方程组得到： m = GT(G GT) -1 dc 先验条件： 反演模型与先验模型方差达到极小。 概述 反问题解的不稳定性对策： — 混定问题。 例： G m = dc 转化为 (GGT + a I )m = dc 求解该方程组得到： m = (GTG + a I ) -1GT dc 阻尼最小二乘法。 概述 反问题解的求解： — 充分利用先验信息。 最下平方反演 Q = || d – G m ||2 —〉min （定量反演） 最大似然反演： Q = P ( m | d ) —〉max （随机反演） 本质是一样的，求解手段方法不同而已。 线性反演理论 线性反演理论 在连续介质条件下： 且有： 在观测数据有限的情况下 式中 表示内积。 线性反演理论 为讨论线性反演问题，对核函数和观测数据作如下假设： 利用线性无关的核函数构成正交函数组 式中 线性反演理论 以 ak j 为系数，对观测数据dj作一线性的组合： Ek是模型 m (z) 在正交坐标基 轴上的投影。 线性反演理论 而在区间 上，任何函数m(z)可以表示为级数： 为正交坐标基，可以分解为两项。 因此有： 线性反演理论 可以证明： 证明过程： 线性反演理论 因此可以得到： 从上式可以看出： 线性反演理论 （1）对于给定观测数据 dj ，只要满足： 就可找到如下式所示的解满足观测数据： 解是存在的 线性反演理论 （2）根据观测数据所构制模型m (z)，由两部分构成： 第一部分与观测数据有关，第二部分与观测数据无关，模型构制是对线性无关核函数实行正交变换，求得正交变换基以及模型在该正交变换基上的投影过程。 线性反演理论 （3）根据观测数据所构制模型m (z)，由两部分构成： 第一部分与观测数据有关，第二部分与观测数据无关，在特殊解基础上加任何零化子向量得到的模型都可拟合观测数据，因此反问题的解是非唯一的。 线性反演理论 （4）根据观测数据所构制模型m (z)，由两部分构成： 得到： 线性反演理论 （4）欲使 最小， ，最小模型能拟合观测数据，又无零化空间影响。在所有模型中L2范数为最小的模型，最小模型核函数的线性组合。 线性反演问题求解 参数化模型： 是指能用有限个参数进行表征的模型，包括连续模型和离散模型。 最小长度： 在处理超定问题和欠定问题时