[發(fā)明專利]地震資料處理方法及其闕值和有效信號(hào)頻率的求取方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201710659366.3
• 申請(qǐng)日: 2017-08-04
• 公開（公告）號(hào): CN107526103B
• 公開（公告）日: 2019-08-02
• 發(fā)明（設(shè)計(jì)）人: 陳立;熊晶璇;唐虎;段鵬飛;李敏;王光銀;李敬益
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/32;G01V1/36
• 代理公司: 成都中璽知識(shí)產(chǎn)權(quán)代理有限公司 51233
• 代理人: 譚昌馳;張敏

本發(fā)明涉及一種地震資料處理方法及其闕值和有效信號(hào)頻率的求取方法。所述闕值求取方法包括獲取地震資料預(yù)定時(shí)窗內(nèi)的頻譜曲線；計(jì)算頻譜曲線上不同頻率的疊加頻譜，然后對(duì)疊加頻譜進(jìn)行求導(dǎo)，得到疊加頻譜變化率曲線；選取疊加頻譜變化率曲線上疊加頻譜變化率最大值所對(duì)應(yīng)的頻率作為閾值。所述有效信號(hào)頻率求取方法包括根據(jù)所述闕值確定出有效頻帶范圍。所述地震資料處理方法采用所述地震資料處理中有效信號(hào)頻率的求取方法求取有效信號(hào)頻率。本發(fā)明能夠定量的得到闕值大小，使得求取的闕值更加準(zhǔn)確有效；并且能夠有效地將噪聲與有效信號(hào)區(qū)分開來(lái)，求取的有效頻率更為準(zhǔn)確、可靠，從而保證了地震資料處理與解釋的準(zhǔn)確性。

技術(shù)領(lǐng)域

本發(fā)明屬于石油地震勘探地震資料處理領(lǐng)域，具體地說(shuō)，涉及一種地震資料處理中闕值的求取方法、地震資料處理中有效信號(hào)頻率的求取方法以及地震資料處理方法。

背景技術(shù)

當(dāng)前，在石油勘探領(lǐng)域，噪聲是影響地震資料品質(zhì)的一個(gè)重要因素，噪聲越多、越強(qiáng)，對(duì)有效信號(hào)的影響越大，疊加產(chǎn)生的地質(zhì)剖面也就越差。在整個(gè)地震資料處理的過(guò)程中，如果能有效地壓制噪聲，保護(hù)有效信號(hào)，處理結(jié)果就能較為真實(shí)的還原地質(zhì)形態(tài)面貌，這對(duì)于地震解釋分析是有利的。但是由于野外采集到的地震資料中，噪聲通常會(huì)與有效信號(hào)混雜在一起，因此需要選取一個(gè)闕值將噪聲與有效信號(hào)區(qū)分開來(lái)，這對(duì)于地震處理的過(guò)程質(zhì)量控制也是十分有必要的。

而在闕值的求取過(guò)程中，通常是根據(jù)經(jīng)驗(yàn)人為地進(jìn)行噪聲與有效信號(hào)的劃分，從而確定闕值大小。這種判定方法很模糊、很主觀，并且闕值的取值缺乏定量的判定和嚴(yán)密的計(jì)算，因此，在這種情況下，處理效果往往會(huì)大打折扣。

因此，需要一種地震資料處理中闕值和有效頻率的求取方法來(lái)對(duì)地質(zhì)勘探過(guò)程中產(chǎn)生的地震數(shù)據(jù)進(jìn)行有效地處理和準(zhǔn)確的解釋。

發(fā)明內(nèi)容

本發(fā)明的目的在于解決現(xiàn)有技術(shù)存在的上述不足中的至少一項(xiàng)。例如，本發(fā)明的目的之一是解決在地震資料處理過(guò)程中確定闕值大小時(shí)的隨機(jī)性與不準(zhǔn)確性，缺乏定量判定和嚴(yán)密計(jì)算的技術(shù)問(wèn)題。

總體來(lái)講，本發(fā)明通過(guò)對(duì)時(shí)窗內(nèi)數(shù)據(jù)進(jìn)行傅里葉變換得到頻譜曲線，頻譜曲線進(jìn)行疊加求和后，進(jìn)而計(jì)算疊加頻譜曲線的變化率的方式求取闕值。

為了實(shí)現(xiàn)上述目的，本發(fā)明一方面提供了一種地震資料處理中闕值的求取方法，所述方法可以包括以下步驟：獲取地震資料預(yù)定時(shí)窗內(nèi)的頻譜曲線；計(jì)算所述頻譜曲線上不同頻率的疊加頻譜，然后對(duì)所述疊加頻譜進(jìn)行求導(dǎo)，得到疊加頻譜變化率曲線；選取所述疊加頻譜變化率曲線上疊加頻譜變化率最大值所對(duì)應(yīng)的頻率作為閾值。

在本發(fā)明的地震資料處理中闕值的求取方法的一個(gè)示例性實(shí)施例中，所述獲取地震資料預(yù)定時(shí)窗內(nèi)的頻譜曲線的步驟可以包括：在地震資料上選取一個(gè)時(shí)窗；對(duì)所選時(shí)窗內(nèi)的地震資料進(jìn)行傅里葉變換，以將所選時(shí)窗內(nèi)的地震資料由時(shí)間域變換到頻率域，得到頻譜曲線。其中，可以通過(guò)等式(1)對(duì)所選時(shí)窗內(nèi)的地震資料進(jìn)行傅里葉變換：

在等式(1)中，x(t)為所選時(shí)窗的地震資料，t表示時(shí)間，單位為秒，w表示頻率，單位為Hz，

在本發(fā)明的地震資料處理中闕值的求取方法的一個(gè)示例性實(shí)施例中，可以通過(guò)等式(2)計(jì)算所述頻譜曲線上不同頻率的疊加頻譜：

E_w＝∑X(w)w＝w₁,w₂,…,w_max 等式(2)

在等式(2)中，E_w為頻譜曲線上不同頻率的疊加頻譜，w表示頻率，w₁,w₂,…,w_max表示頻譜曲線上的不同頻率，w_max表示頻譜曲線上頻率的最大值。

在本發(fā)明的地震資料處理中闕值的求取方法的一個(gè)示例性實(shí)施例中，可以通過(guò)等式(3)對(duì)所述疊加頻譜進(jìn)行求導(dǎo)，得到疊加頻譜變化率曲線：