第五章地球物理反演理论-非线性反演问题.pps-武汉大学.ppt

多尺度反演（Multi-Scale Inversion） 1．尺度的概念 尺度是指当我们以离散方式描述某一空间（或时间）函数时，均匀离散点之间的距离。尺度是分辨率的倒数。分辨率被定义为单位距离内离散点的个数。 多尺度反演时把目标函数分解成不同尺度的分量，根据不同尺度上目标函数的特征逐步搜索全局极小。一般情况下，在大尺度（或低波数）上，目标函数极值点少，且分得很开，且显示大尺度上的极小。用通常的方法很容易直接搜索出大尺度（总体背景）上的“全局极小点”。 多尺度反演（Multi-Scale Inversion） 1．尺度的概念 在相对较小尺度（稍大波数）上，目标函数极值点较多，如无大尺度搜索结果之指导，直接寻找全局极值点，虽比作尺度分解之前容易，但仍然比较困难。但是，只要我们以大尺度搜索到的背景“全局极小点”为起点，在其附近继续搜索，也较容易的搜索到中等尺度上到的“全局极小点”。如此，不断的缩小尺度，提高分辨力，目标函数的尺度降至原始尺度（即最下尺度）时，对应搜索出的全局极小点，就是真正的总体极小点。 多尺度反演（Multi-Scale Inversion） 1．尺度的概念 优点是，在大尺度（或低波数）上，反演稳定，反演结果不受初始模型的影响，在一定程度上，能避免其后的反演受局部极小所困扰，使收敛速速加快。 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 Goupillaud，Grossmanne和Morlet（1984）引进了小波的概念。 根据定义，我们称满足条件： （5.31） 的函数 为小波函数或母小波，式中 是 的傅氏变换。 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 若 （5.32） 式中： 为由母小波 生成的依赖于参数a,b的连续小波；a称为尺度变量；b为位置变量。尺度变量a的改变，标志 着相对于母小波发生了伸缩改变；而位置变量b发生变化时，标志着 相对于母小波发生了平移。 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 若 （5.32） 式中： 是归一化因子，它使所有的小波具有相同的模。因此，所有的连续小波 都具有相同的能量。 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 对于任一 的函数而言，其小波变换定义为： （5.33） 式中： 为内积；而 与 是共轭。 其逆变换公式为： （5.34） 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 式中： （5.35） 在实际计算中，常用其离散形式，若令： ， 则有著名的二进制小波，即： （5.36） 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 二进制小波构成L2(R)的一个正交基。利用 可以将在无穷大处衰减得充分快的任意函数f(t)分解为： （5.37） （5.38） 设： 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 则（5.37）式可写为： （5.39） 基于（5.39）式的分析方法称之为尺度分析法。 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 如果基于尺度上对函数f(t)的光滑近似属于空间Vj，f(t)的细节部分属于的正交补空间Wj，则有： （5.40） 当j=1时， （5.41） 这种分解法可用图5-7表示： 图5-7 多尺度分解示意图 多尺度反演（Multi-Scale Inversion） 2．小波与多尺度分析 当我们把反演问题从Vj空间分解为时，由于Vj+1空间仅包含Vj空间中的低波数成分。因此，在Vj+1空间求得的解应是Vj空间中反问题在所有模型参数附近局部化的某种平均值解。又由于Vj+1在Vj中的正交空间Wj+1含有Vj空间中的高波数成分，因此在Vj+1空间中不能得到分辨的小尺度特征，却能通过在Vj空间中的求解得到分辨。这就是多尺度方法逐步提高分辨率的原理。 多尺度反演（Multi-Scale Inversion） 3．多尺度反演法 多尺度反演算法可以表示为三个基本算子的操作过程：第一个算子将地球物理反演问题从尺度零（小尺度）开始，依次分解为尺度1，2，…等大尺度的问题；第二个算子求取各尺度上反问题的解；第三个算