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

有效

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

技術領域

本發明屬于石油地震資料勘探領域，更具體地講，涉及一種自動、有效地去除采集的地震數據中包含的低頻異常振幅干擾，同時提高地震數據處理效率的自適應低頻異常振幅壓制方法。

背景技術

低頻異常振幅是石油地震資料勘探中廣泛存在的一種干擾波，它主要集中在地震記錄的近偏移距附近，具有低頻、低速、衰減較慢的特點，強低頻異常振幅的存在，將嚴重地影響反褶積和疊加效果，對保幅處理尤為不利。在地震資料處理中低頻異常振幅壓制是提高地震資料信噪比和分辨率的重要步驟，因此在地震資料處理中如何去除或提取低頻異常振幅就成了目前地震勘探數據處理的主要問題之一。

目前在地震資料處理過程中大多采用高通濾波的方法來壓制低頻異常振幅，主要是利用低頻異常振幅和反射波的頻率差異及速度差異進行濾波，如一維濾波和頻率波數濾波，但因低頻異常振幅的頻帶往往和反射波中的低頻帶攪在一起，即便是頻率波數域濾波也不能只對與有效信號相同頻帶的噪聲進行壓制，其實在濾波過程的同時也對有效波的低頻成分造成了損害，而且斜坡長的濾波響應在反褶積或譜白化處理后，低頻異常振幅往往又被突出出來。因此，需要一種能夠有效的去除采集的地震數據中包含的低頻異常振幅干擾，同時提高地震數據處理效率的地震數據處理方法。

發明內容

本發明的目的在于提供一種能夠自動去除采集的地震數據中包含的低頻異常振幅干擾，同時提高地震數據處理效率的自適應低頻異常振幅壓制方法，從而滿足日漸興起的高分辨率勘探的精度要求。

本發明的一方面提供一種自適應低頻異常振幅壓制的方法，包括：對單炮地震數據的每道進行傅立葉正變換，得到頻率域的地震數據；從頻率域的地震數據分離出低頻段的地震數據；將分離出的低頻段的地震數據進行傅立葉反變換，得到時間域的低頻地震數據；從時間域的低頻地震數據提取低頻異常振幅分布區數據、淺層折射區數據以及無干擾分布區數據；根據提取的無干擾分布區數據獲取低頻異常振幅閥值；根據低頻異常振幅閥值對提取的低頻異常振幅分布區數據進行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區數據；從單炮地震數據中減去低頻異常振幅分布區數據與低頻異常振幅壓制后的低頻異常振幅分布區數據之差，得到低頻異常振幅壓制后的地震數據。

可選地，所述低頻段的頻率范圍為0-20HZ。

可選地，根據提取的無干擾分布區數據獲取低頻異常振幅閥值的步驟包括：對提取的無干擾分布區數據進行傅立葉正變換，得到頻率域的無干擾分布區數據；計算頻率域的無干擾分布區數據的低頻異常振幅閥值。

可選地，根據提取的無干擾分布區數據獲取低頻異常振幅閥值的步驟包括：將無干擾分布區數據沿著時間的方向劃分為多個時間窗口；對每個時間窗口的數據進行傅立葉正變換，得到相應的頻率域的地震數據；將得到的每個時間窗口的頻率域的地震數據劃分為多個頻段；計算每個時間窗口的每個頻段的低頻異常振幅閥值。

可選地，根據低頻異常振幅閥值對提取的低頻異常振幅分布區數據進行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區數據的步驟包括：對提取的低頻異常振幅分布區數據進行二維分數階傅立葉變換，將提取的低頻異常振幅分布區數據旋轉變換到分數域的一個時頻平面；從旋轉變換得到的數據查找低頻段的振幅大于低頻異常振幅閥值的區域；在查找到的區域填充預定值，從而獲取時頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區數據；對獲取的時頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區數據進行二維分數階傅立葉反變換，得到時間域的低頻異常振幅壓制后的低頻異常振幅分布區數據。

可選地，對提取的低頻異常振幅分布區數據進行二維分數階傅立葉變換，將提取的低頻異常振幅分布區數據旋轉變換到分數域的一個時頻平面的步驟包括：將低頻異常振幅分布區數據沿著時間的方向劃分為多個時間窗口；對每個時間窗口的數據進行二維分數階傅立葉變換，得到相應的分數域的地震數據。

可選地，從旋轉變換得到的數據查找低頻段的振幅大于低頻異常振幅閥值的區域的步驟包括：將每個時間窗口的分數域的地震數據劃分為多個頻段；將每個時間窗口的每個頻段的振幅值與每個時間窗口的每個頻段的低頻異常振幅閥值進行比較，在每個時間窗口的每個頻段中查找振幅值大于對應的低頻異常振幅閥值的區域。

可選地，所述預定值為零。

可選地，在二維分數階傅立葉變換中的旋轉角度為pπ/2，其中，p為二維分數階傅立葉變換的階數。