[發明專利]自適應低頻異常振幅壓制的方法

有效

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

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

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

從頻率域的地震數據分離出低頻段的地震數據；

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

從時間域的低頻地震數據提取異常振幅分布區數據和無干擾分布區數據；

根據提取的無干擾分布區數據獲取低頻異常振幅閥值；

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

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

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

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

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

（b）計算頻率域的無干擾分布區數據的低頻異常振幅閥值。

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

將無干擾分布區數據沿著時間的方向劃分為多個時間窗口；

對每個時間窗口的數據進行傅立葉正變換，得到相應的頻率域的地震數據；

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

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

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

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

（b）從旋轉變換得到的數據查找低頻段的振幅大于低頻異常振幅閥值的區域；

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

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

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

將低頻異常振幅分布區數據沿著時間的方向劃分為多個時間窗口；

對每個時間窗口的數據進行二維分數階傅立葉變換，得到相應的分數域的地震數據。

7.如權利要求6所述的方法，其中，從旋轉變換得到的數據查找低頻段的振幅大于低頻異常振幅閥值的區域的步驟包括：

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

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

8.如權利要求5所述的方法，其中，所述預定值為零。

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

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