[发明专利]一种地震波分解方法

摘要

一种地震波分解的方法，该方法包括选择以Gabor小波或Morlet小波作为子波形式；对获取的地震信号进行多子波的分解；在每一次分解，通过地震信号的复数属性计算小波的三个参数：时移u、中心频率ω、相移φ；根据这三个参数计算时间尺度参数σ；估计小波的振幅a；通过迭代确定最终参数及分解出不同主频的子波。本发明提供的地震子波分解方法，对传统的匹配追踪方法进行了改进，提高了计算效率；同时σ是可自适应可调节的。

主权项

一种地震波分解方法，该方法包括：1)选择子波形式，该子波形式为由公式(6)定义的Morlet小波形式： m ( t ) = exp [ - ( ln 2 π 2 ) ω m 2 ( t - u ) 2 σ 2 ] exp [ i ( ω m ( t - u ) + φ ) ] 公式(6)，其中m(t)是Morlet小波函数、t是时间变量、u是时移、σ是时间轴长度、ωm是平均角频率、φ是相移、i表示虚数单位；2)获取地震信号，并对获取的地震信号进行多次分解，直至满足分解完成条件，分解后的地震信号的形式为公式(7)所定义的形式： f ( t ) = Σ n = 1 N a n m n ( t ) + R ( N + 1 ) f 公式(7)，其中f(t)是地震信号的函数，n为分解次数，an是与第n次分解的小波mn对应的振幅值，N为正整数，R(N+1)f为第N次分解后的残差， R ( j ) f = a j g r j + R ( j + 1 ) f 且R(0)f＝f，标记j为1至N的正整数；根据所述地震信号的复数属性来计算第j次分解时的参数集 γ j = { u j , σ j , ω m j , φ j } ; 在计算出所述参数集 γ j = { u j , σ j , ω m j , φ j } 的基础上计算第j次分解时的小波的振幅值aj；3)将所述计算出的 γ j = { u j , σ j , ω m j , φ j } 代入到公式(6)以得到在第j次分解出的小波mj，将计算出的振幅值aj与该小波mj相乘以得到第j次分解出的子波。