基于小波分析的瞬变电磁法资料去噪技术.pdfVIP

增刊 解海军等：基于小波分析的瞬变电磁法资料去噪技术 I哕(菇)I (1)J 2dx∞，即平方可积，记为缈(并)∈L2(R)； (2)f”{警#幽∞，奴山)为妒(髫)的F0u五er变换． J一∞ l 6口f 小波变换是将信号以t)∈r(R)用小波函数集分解的运算．小波变换定义为 咖6)甜‰兆一肋z)以t号-b，)dt， 其逆变换为 州=长亡咖，b)g,o“t)争” 其中 c妒=仁哗华吵 实际应用中，常对尺度因子口进行二进离散，离散化后的小波变换为 加：以戈)=鸬t(菇)=方仁kt)qx2-。t，ldt． E1，2】 1．2 Mallat ．设巧为一个给定的多尺度分析，妒和哕分别为相应的尺度函数和小波函数． 设～个函数(信号)，∈ 巧．(，t为一确定整数)，数学上已证明 以戈)=Aj∥搿)=Ajl+lA菇)+DJl+l八石)， 其中 =∑C小l，。，9小l川 tit=一∞ +∞ D¨l以菇) =∑D』l+。，。，奶1+l，。， 而 C』l+l，。2 ∑hk_2mC”k， 2 D，J+l，m ∑gk-2mC，l，k． 写成简洁形式即为 (2) C^+l=ncj，DJl+l=Gq(，=J．，L1+l，…，以)． 式(2)便是Mallat分解算法． 由分解算法知，e可称为在型分辨率下的连续(离散)逼近，而Di可称为在2j分辨率下的连续(离 散)细节．由Mallat分解算法可见，函数。(信号)可被分解成不同频率通道的成分． 式(1)经两边同时与某个特定函数作内积后，再经数学处理，得到Mallat重构算法，即 cf=日+c『十l+G+2+l(J=以一1，…，．，l+1，J。)， 式中，日+和G+为日和G的共轭转置矩阵． ． 2基于小波分析的瞬变电磁法资料处理 2．1奇异信号在小波变换中的表象‘卜3】 奇异信号在瞬变电磁法深度剖面上表现为突变异常或者整体的抬升．那么如何判断哪些突变点是静态 煤 炭 学 报 2006年第31卷 a为正时，对应的异常为地质构造；当a为负时，对应的异常为干扰所引起的异常．下式是为了求a所构 造的目标函数． ， 2]2， k-jotlg E(k，a)=∑[19w2Jp。(∞)|-lg 其中应满足的条件为 2， lg k+jodg l w2ip。(∞)l≤lg 这里的训印。(o．1)为突变点的二进小波变换值． 2．2小波分析压制奇异信号的实际效果‘4，5】 本文选用的是对函数的奇异性检测很适合的具有一阶消失矩的小波函数缈