[發(fā)明專(zhuān)利]一種消除偽頻譜的頻譜分析方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210418514.X
• 申請(qǐng)日: 2012-10-26
• 公開(kāi)（公告）號(hào): CN103792577A
• 公開(kāi)（公告）日: 2014-05-14
• 發(fā)明（設(shè)計(jì)）人: 劉志成;謝金娥;方伍寶;賈春梅;宋林
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油物探技術(shù)研究院
• 主分類(lèi)號(hào): G01V1/30
• 分類(lèi)號(hào): G01V1/30
• 代理公司: 北京思創(chuàng)畢升專(zhuān)利事務(wù)所 11218
• 代理人: 劉明華

1.一種消除偽頻譜的頻譜分析方法，其特征在于，所述頻譜分析方法包括剔除偽頻譜步驟，所述偽頻譜為數(shù)字信號(hào)處理中正常頻譜以外存在的與正常頻譜呈周期關(guān)系的頻譜；所述剔除偽頻譜包括對(duì)數(shù)字信號(hào)進(jìn)行希爾伯特變換，然后進(jìn)行傅里葉變換，且將經(jīng)過(guò)希爾伯特變換的結(jié)果數(shù)據(jù)設(shè)置為傅里葉變換時(shí)輸入數(shù)據(jù)的虛部，經(jīng)過(guò)傅里葉變換后輸出頻譜，所述偽頻譜數(shù)字信號(hào)被剔除。

2.根據(jù)權(quán)利要求1所述的一種消除偽頻譜的頻譜分析方法，其特征在于，對(duì)所述數(shù)字信號(hào)x(t)進(jìn)行希爾伯特變換得到h(t)，表達(dá)式為：

h(t)=1π∫-∞+∞x(τ)t-τdτ---(1)]]>

其中，t為時(shí)間，τ為時(shí)間延遲，+∞為正無(wú)窮，-∞為負(fù)無(wú)窮；

所述傅里葉變換采用公式：

F(ω)=∫-∞+∞f(t)eiωtdt---(2)]]>

其中t為時(shí)間，ω為圓頻率，i為虛數(shù)，+∞為正無(wú)窮，-∞為負(fù)無(wú)窮；

采用歐拉公式e^±ix＝cosx±isinx簡(jiǎn)化（2）式，故（2）式可寫(xiě)為：

F(ω)=∫-∞+∞f(t)(cos(ωt)-isin(ωt))dt---(3)]]>

（3）式虛部為sin(ωt)，令sin(ωt)＝h(t)。

3.應(yīng)用權(quán)利要求1-2的方法進(jìn)行地震信號(hào)頻譜分析的方法，其特征在于，所述方法包括：

步驟（1）輸入含噪地震數(shù)據(jù)信號(hào)道集；

（2）對(duì)所述含噪地震數(shù)據(jù)信號(hào)道集進(jìn)行希爾伯特變換:

含噪地震數(shù)字信號(hào)x(t)進(jìn)行希爾伯特變換得到h(t)，

h(t)=1π∫-∞+∞x(τ)t-τdτ---(1)]]>

其中，t為時(shí)間，τ為時(shí)間延遲，+∞為正無(wú)窮，-∞為負(fù)無(wú)窮；

（3）進(jìn)行傅里葉變換：

完成步驟（2）后，進(jìn)行傅里葉變換，且設(shè)h(t)為傅里葉變換中的虛部；所述傅里葉變換的公式為：所述傅里葉變換采用公式：

F(ω)=∫-∞+∞f(t)eiωtdt---(2)]]>

其中t為時(shí)間，ω為圓頻率，i為虛數(shù)，+∞為正無(wú)窮，-∞為負(fù)無(wú)窮；

（4）輸出頻譜：將經(jīng)過(guò)傅里葉變換后的含噪地震數(shù)據(jù)信號(hào)頻譜圖輸出；

（5）頻譜分析：分析輸出的頻率-振幅譜，橫軸是頻率，縱軸是振幅。