專利名稱：強(qiáng)能量干擾抑制方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種地震數(shù)據(jù)處理方法，更具體地講，涉及一種在處理地震數(shù)據(jù)時(shí)使用的強(qiáng)能量干擾抑制方法。
背景技術(shù)：
近年來(lái)，隨著油氣勘探形勢(shì)的發(fā)展，地震勘探的范圍也越來(lái)越廣泛，越來(lái)越深入， 地震勘探工作也逐漸從東部轉(zhuǎn)向西部，由平原轉(zhuǎn)向沙漠、丘陵、山地等復(fù)雜地表區(qū)，在這些復(fù)雜的地貌地區(qū)采集的地震資料中，普遍存在著相當(dāng)嚴(yán)重的強(qiáng)能量干擾，例如聲波干擾、 交流電干擾、脈沖噪聲及在記錄深部出現(xiàn)的高頻干擾，一般情況下這類噪聲嚴(yán)重干擾了有效反射信號(hào)，極大地降低了地震資料的信噪比。對(duì)于這些強(qiáng)能量干擾，目前常用的處理方法是濾波或道剔除，但無(wú)論是交互剔除廢道還是人工編輯，都需要花費(fèi)大量的時(shí)間，而且淹沒(méi)在能量干擾中的有效信號(hào)也會(huì)損失， 因此目前的干擾處理方法已不能滿足日漸興起的高分辨率勘探的精度要求，要解決這個(gè)問(wèn)題，只有尋求自動(dòng)識(shí)別和剔除這些異常能量的新方法。因此需要一種能夠有效的去處采集的地震數(shù)據(jù)中包含的強(qiáng)能量干擾，同時(shí)提高大量地震數(shù)據(jù)處理效率的地震數(shù)據(jù)處理方法。

發(fā)明內(nèi)容
根據(jù)本發(fā)明的一方面，提供了一種強(qiáng)能量干擾抑制方法，所述方法包括(a)將單炮地震數(shù)據(jù)沿地震道與時(shí)間的方向劃分為多個(gè)窗口 ；(b)分別對(duì)每個(gè)窗口的地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉變換，在分?jǐn)?shù)頻域中獲得與每個(gè)窗口的地震數(shù)據(jù)對(duì)應(yīng)的頻域地震數(shù)據(jù)；(C) 基于獲得的各個(gè)窗口的頻域地震數(shù)據(jù)，為每個(gè)窗口建立強(qiáng)能量衰減函數(shù)；(d)使用為每個(gè)窗口建立的強(qiáng)能量衰減函數(shù)，分別對(duì)每個(gè)窗口的頻域地震數(shù)據(jù)進(jìn)行衰減；(e)分別對(duì)衰減后的每個(gè)窗口的頻域地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉反變換，以獲得去除了強(qiáng)能量干擾的地震數(shù)據(jù)。此外，可通過(guò)以下等式對(duì)每個(gè)窗口的地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉變換
權(quán)利要求
1.一種強(qiáng)能量干擾抑制方法，包括(a)將單炮地震數(shù)據(jù)沿地震道與時(shí)間的方向劃分為多個(gè)窗口；(b)分別對(duì)每個(gè)窗口的地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉變換，在分?jǐn)?shù)頻域中獲得與每個(gè)窗口的地震數(shù)據(jù)對(duì)應(yīng)的頻域地震數(shù)據(jù)；(C)基于獲得的各個(gè)窗口的頻域地震數(shù)據(jù)，為每個(gè)窗口建立強(qiáng)能量衰減函數(shù)；(d)使用為每個(gè)窗口建立的強(qiáng)能量衰減函數(shù)，分別對(duì)每個(gè)窗口的頻域地震數(shù)據(jù)進(jìn)行衰減；(e)分別對(duì)衰減后的每個(gè)窗口的頻域地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉反變換，以獲得去除了強(qiáng)能量干擾的地震數(shù)據(jù)。
2.如權(quán)利要求I所述的方法，其中，通過(guò)以下等式對(duì)每個(gè)窗口的地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉變換
3.如權(quán)利要求2所述的方法，其中，對(duì)于所述多個(gè)窗口中的第k個(gè)窗口的頻域地震數(shù)據(jù)，建立強(qiáng)能量衰減函數(shù)的步驟包括將第k個(gè)窗口的頻域地震數(shù)據(jù)fk(x，u)劃分為多個(gè)頻段的頻域地震數(shù)據(jù)fk(Xl，u).....fk(xi7 u)、. . .、fk(xm, u),其中，m是劃分的頻段的數(shù)量，fk(Xi, u)表示第k個(gè)窗口的第i個(gè)頻段的頻域地震數(shù)據(jù)，并且KiSm;基于以下等式計(jì)算第k個(gè)窗口的每個(gè)頻段的頻域地震數(shù)據(jù)的振幅值，Ak(Xi, u) = I (fk(xi； u)) I,其中，Ak(Xi，u)表示第k個(gè)窗口的第i個(gè)頻段的頻域地震數(shù)據(jù)fk(Xi，u)的振幅值；基于計(jì)算的第k個(gè)窗口的多個(gè)頻段的頻域地震數(shù)據(jù)的振幅值，計(jì)算第k個(gè)窗口的多個(gè)頻段的頻域地震數(shù)據(jù)的振幅值的加權(quán)中值Mk(U)；根據(jù)計(jì)算的第k個(gè)窗口的頻域地震數(shù)據(jù)的振幅值的加權(quán)中值Mk (U)，計(jì)算第k個(gè)窗口的加權(quán)因子Sk(u) = ThrXMk(U)，其中，Thr是預(yù)設(shè)的閾值；基于以下等式為第k個(gè)窗口的每個(gè)頻段的頻域地震數(shù)據(jù)建立強(qiáng)能量衰減函數(shù)，
4.如權(quán)利要求3所述的方法，其中，對(duì)于所述多個(gè)窗口中的第k個(gè)窗口的頻域地震數(shù)據(jù)，進(jìn)行衰減的步驟包括根據(jù)以下等式對(duì)第k個(gè)窗口的各頻段的頻域地震數(shù)據(jù)進(jìn)行衰減，
5.如權(quán)利要求4所述的方法，其中，步驟(e)包括對(duì)各窗口的衰減后的每個(gè)頻段的頻域地震數(shù)據(jù)進(jìn)行組合，以獲得衰減后的每個(gè)窗口的頻域地震數(shù)據(jù)；分別對(duì)衰減后的各窗口的頻域地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉反變換，以獲得去除了強(qiáng)能量干擾的地震數(shù)據(jù)。
6.如權(quán)利要求3所述的方法，其中，計(jì)算第k個(gè)窗口的多個(gè)頻段的頻域地震數(shù)據(jù)的振幅值的加權(quán)中值的步驟包括為每個(gè)頻段的頻域地震數(shù)據(jù)的振幅值分別設(shè)置加權(quán)系數(shù)W、W2.....Wi.....wm，其中，Wi是為第i個(gè)頻段的頻域地震數(shù)據(jù)的振幅值A(chǔ)k(Xi，u)設(shè)置的加權(quán)系數(shù)，并且加權(quán)系數(shù)Wp w2、. . .、Wi、. . .、Wm 均為正整數(shù)；根據(jù)設(shè)置的加權(quán)系數(shù)，對(duì)振幅值序列
7.如權(quán)利要求I所述的方法，其中，在可編程圖形處理單元平臺(tái)下通過(guò)多線程執(zhí)行步驟(a)到(e)中的一個(gè)或多個(gè)步驟。
全文摘要
提供了一種強(qiáng)能量干擾抑制方法，所述方法包括將單炮地震數(shù)據(jù)沿地震道與時(shí)間的方向劃分為多個(gè)窗口；分別對(duì)每個(gè)窗口的地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉變換，在分?jǐn)?shù)頻域中獲得與每個(gè)窗口的地震數(shù)據(jù)對(duì)應(yīng)的頻域地震數(shù)據(jù)；基于獲得的各個(gè)窗口的頻域地震數(shù)據(jù)，為每個(gè)窗口建立強(qiáng)能量衰減函數(shù)；使用為每個(gè)窗口建立的強(qiáng)能量衰減函數(shù)，分別對(duì)每個(gè)窗口的頻域地震數(shù)據(jù)進(jìn)行衰減；分別對(duì)衰減后的每個(gè)窗口的頻域地震數(shù)據(jù)進(jìn)行分?jǐn)?shù)階傅里葉反變換，以獲得去除了強(qiáng)能量干擾的地震數(shù)據(jù)。
文檔編號(hào)G01V1/36GK102608658SQ20111042498
公開(kāi)日2012年7月25日 申請(qǐng)日期2011年12月16日 優(yōu)先權(quán)日2011年12月16日
發(fā)明者何光明, 劉鴻, 呂文彪, 張華 , 曹中林, 李振, 羅紅明, 陳愛(ài)萍 申請(qǐng)人:中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司