[發(fā)明專利]一種地震信號(hào)處理方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201611071443.5
• 申請(qǐng)日: 2016-11-29
• 公開(kāi)（公告）號(hào): CN108072894A
• 公開(kāi)（公告）日: 2018-05-25
• 發(fā)明（設(shè)計(jì)）人: 朱博華;向雪梅;楊江峰;楊勤林;李洋;陳哲;張薇;胡瑋;齊晴
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油物探技術(shù)研究院
• 主分類(lèi)號(hào): G01V1/28
• 分類(lèi)號(hào): G01V1/28
• 代理公司: 北京聿宏知識(shí)產(chǎn)權(quán)代理有限公司 11372
• 代理人: 吳大建

1.一種地震信號(hào)處理方法，所述方法以希爾伯特-黃變換為基礎(chǔ)，將原始的地震信號(hào)通過(guò)經(jīng)驗(yàn)?zāi)B(tài)分解得到本征模式函數(shù)分量，分別對(duì)每個(gè)本征模式函數(shù)分量進(jìn)行振幅增強(qiáng)處理，獲得時(shí)間域本征模式函數(shù)分量信號(hào)，對(duì)時(shí)間域本征模式函數(shù)分量信號(hào)進(jìn)行疊加，然后進(jìn)行多道均衡處理獲得最終的地震信號(hào)。

2.根據(jù)權(quán)利要求1所述的地震信號(hào)處理方法，其特征在于，所述方法包括如下步驟：

2.1設(shè)計(jì)白化濾波器，白化濾波器表達(dá)式如下

其中，為外包絡(luò)振幅值，v是外包絡(luò)的最大值即ε為白噪因子，ω表示頻率；

2.2對(duì)第j個(gè)本征模式函數(shù)分量的瞬時(shí)振幅進(jìn)行振幅增強(qiáng)，獲得白化后的瞬時(shí)振幅，所述白化后的瞬時(shí)振幅表達(dá)式如下

a ~ j ( t ) = a j ( t ) f ( ω j ( t ) ) ]]>

其中，aj(t)為第j個(gè)本征模式函數(shù)分量的瞬時(shí)振幅，ωj(t)為第j個(gè)本征模式函數(shù)分量的瞬時(shí)頻率，j＝1,2,3…n；

2.3獲得重構(gòu)后時(shí)間域本征模式函數(shù)分量，表達(dá)式如下

s ~ j ( t ) = Re [ a ~ j ( t ) e iθ j ( t ) ] ]]>

其中，e為自然常數(shù)2.71828，i為虛數(shù)符號(hào)，為第j個(gè)本征模式函數(shù)分量譜白化后的瞬時(shí)振幅，θj(t)是對(duì)應(yīng)的瞬時(shí)相位；

2.4將時(shí)間域本征模式函數(shù)分量疊加，獲得譜白化后的地震信號(hào)，表達(dá)式為

y ( t ) = Σ j = 1 n s ~ j ( t ) ]]>

2.5對(duì)譜白化后的單道地震信號(hào)進(jìn)行多道均衡處理，獲得最終的地震信號(hào)，多道均衡處理表達(dá)式為

z m ( t ) = ( Σ j = m - 1 m + 1 y j ( t ) ) / 3 , m = 2 , 3 , 4 ... n - 1 ]]>