[發(fā)明專利]自適應(yīng)低頻異常振幅壓制的方法

有效

• 申請?zhí)枺?/td>: 201410083083.5
• 申請日: 2014-03-07
• 公開（公告）號: CN103869361A
• 公開（公告）日: 2014-06-18
• 發(fā)明（設(shè)計）人: 曹中林;何光明;陳愛萍;呂文彪;周強;張華;李振
• 申請（專利權(quán)）人: 中國石油集團川慶鉆探工程有限公司地球物理勘探公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京銘碩知識產(chǎn)權(quán)代理有限公司 11286
• 代理人: 張川緒;金光軍

1.一種自適應(yīng)低頻異常振幅壓制的方法，所述方法包括以下步驟：

對單炮地震數(shù)據(jù)的每道進行傅立葉正變換，得到頻率域的地震數(shù)據(jù)；

從頻率域的地震數(shù)據(jù)分離出低頻段的地震數(shù)據(jù)；

將分離出的低頻段的地震數(shù)據(jù)進行傅立葉反變換，得到時間域的低頻地震數(shù)據(jù)；

從時間域的低頻地震數(shù)據(jù)提取異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)；

根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值；

根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)；

從單炮地震數(shù)據(jù)中減去低頻異常振幅分布區(qū)數(shù)據(jù)與低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)之差，得到低頻異常振幅壓制后的地震數(shù)據(jù)。

2.如權(quán)利要求1所述的方法，其中，所述低頻段的頻率范圍為0-20HZ。

3.如權(quán)利要求1所述的方法，其中，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟包括：

（a）對提取的無干擾分布區(qū)數(shù)據(jù)進行傅立葉正變換，得到頻率域的無干擾分布區(qū)數(shù)據(jù)；

（b）計算頻率域的無干擾分布區(qū)數(shù)據(jù)的低頻異常振幅閥值。

4.如權(quán)利要求3所述的方法，其中，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟包括：

將無干擾分布區(qū)數(shù)據(jù)沿著時間的方向劃分為多個時間窗口；

對每個時間窗口的數(shù)據(jù)進行傅立葉正變換，得到相應(yīng)的頻率域的地震數(shù)據(jù)；

將得到的每個時間窗口的頻率域的地震數(shù)據(jù)劃分為多個頻段；

計算每個時間窗口的每個頻段的低頻異常振幅閥值。

5.如權(quán)利要求1所述的方法，其中，根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)的步驟包括：

（a）對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進行二維分數(shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分數(shù)域的一個時頻平面；

（b）從旋轉(zhuǎn)變換得到的數(shù)據(jù)查找低頻段的振幅大于低頻異常振幅閥值的區(qū)域；

（c）在查找到的區(qū)域填充預(yù)定值，從而獲取時頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)；

（d）對獲取的時頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)進行二維分數(shù)階傅立葉反變換，得到時間域的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)。

6.如權(quán)利要求5所述的方法，其中，對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進行二維分數(shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分數(shù)域的一個時頻平面的步驟包括：

將低頻異常振幅分布區(qū)數(shù)據(jù)沿著時間的方向劃分為多個時間窗口；

對每個時間窗口的數(shù)據(jù)進行二維分數(shù)階傅立葉變換，得到相應(yīng)的分數(shù)域的地震數(shù)據(jù)。

7.如權(quán)利要求6所述的方法，其中，從旋轉(zhuǎn)變換得到的數(shù)據(jù)查找低頻段的振幅大于低頻異常振幅閥值的區(qū)域的步驟包括：

將每個時間窗口的分數(shù)域的地震數(shù)據(jù)劃分為多個頻段；

將每個時間窗口的每個頻段的振幅值與每個時間窗口的每個頻段的低頻異常振幅閥值進行比較，在每個時間窗口的每個頻段中查找振幅值大于對應(yīng)的低頻異常振幅閥值的區(qū)域。

8.如權(quán)利要求5所述的方法，其中，所述預(yù)定值為零。

9.如權(quán)利要求5所述的方法，其中，在二維分數(shù)階傅立葉變換中的旋轉(zhuǎn)角度為pπ/2，其中，p為二維分數(shù)階傅立葉變換的階數(shù)。

10.如權(quán)利要求1所述的方法，其中，低頻異常振幅分布區(qū)數(shù)據(jù)是時間域的低頻地震數(shù)據(jù)中低頻異常振幅集中分布的區(qū)域中的地震數(shù)據(jù)，無干擾分布區(qū)數(shù)據(jù)是時間域的低頻地震數(shù)據(jù)中不存在低頻異常振幅干擾的區(qū)域中的地震數(shù)據(jù)。