小波奇异性检测课件.pptVIP

* * 小波细节的微分特征及其在重力场断裂分析中的应用 断裂分析是地球物理资料解释中的一个重要环节。人们往往根据断裂在地球物理场上所产生的各种异常特征来进行定性和定量解释，还采用相关分析，水平一阶方向导数，垂向二阶导数，总梯度模，希尔伯特变换等各种方法来推断解释断裂位置、产状及性质[1－3]。为了分析深层的断裂或判断断裂向下延深的情况，往往通过向上延拓的方法来分析，但是，随着延拓高度增加，异常变得平坦光滑，即使再求水平一阶导数，其特征也变得模糊，不仅难以可靠地分析断裂空间位置，而且异常的信息也大大减少，这就是以往断裂分析中存在的问题。 小波分析是近年来发展起来的一种新的数学分析方法，广泛应用于信号与图象处理和石油地震勘探。利用小波分析人们成功地实现了地震勘探数据的压缩，地震道的奇性反演，分析石油储层的详细结构，信号重建与压制面波以及提高地震资料的信噪比与分辨率；利用小波多尺度分析方法对重力异常进行分解和多尺度反演，显示了小波分析方法在地球物理信号处理方面的广阔前景[3－6]。本文利用小波细节的微分特征对重力场进行断裂分析。 小波分析及其微分特征 ——小波分析原理及Mallat算法 对信号f (t)∈L2 (R)，其连续小波变换为 ， （1） 其中，ψ(t) ∈L2(R) 称为小波母函数，ψa,b(t)为由ψ(t)生成的依赖于参数ａ和ｂ的连续小波。在ψ(t)的傅立叶变换 满足条件 的前提下，有小波逆变换公式： ， （2） 取 ，则有小波变换离散形式 ， （3） 相应的小波逆变换为 ， （4） 在多分辨分析理论基础上，Mallat提出了所谓的塔式分解算法，简述如下： 设{Vj}是一多分辨分析，分别是相应的尺度函数和小波函数，对于整数J1J2∈Z，函数，有以下分解 ， （5） 其中 ， （6） ， （7） 而 ， （8） ， （9） 定义无穷矩阵和，则（8）和（9）式可写成矩阵形式 小波分析及其微分特征 ——小波分析原理及Mallat算法 ， （10） 对（10）式两端同时与函数作内积，则可得Mallat重构算法： ， （11） 上述一维空间的多分辨分析及Mallat算法可以很容易推广到二维空间[11]，只要将原来的{Vj}用{Vj2}取代就可以了。对于信号f(x,y)有类似一维的Mallat分解算法： ， （12） 和重构算法： ， （13） 小波分析及其微分特征 ——小波分析原理及Mallat算法 小波分析及其微分特征 ——小波细节的微分特征与信号奇异性检测 在信号处理中，信号中的奇异点及不规则的突变部分经常带有比较重要的信息，它是信号重要的特征之一。通常情况下，信号奇异性分两类：一类是信号在某一时刻内，其幅值发生突变；一类是信号幅值没有突变，但