[發明專利]一種地震記錄的方波壓制方法

有效

• 申請號: 202110102152.2
• 申請日: 2021-01-26
• 公開（公告）號: CN112925024B
• 公開（公告）日: 2023-10-20
• 發明（設計）人: 周開明;郭利榮;霍浩;陳鵬;彭連滔;黃倩
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司江漢油田分公司物探研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 湖北武漢永嘉專利代理有限公司 42102
• 代理人: 劉秋芳;張宇

本發明公開了一種地震記錄的方波壓制方法，屬于石油勘探開發領域，該方法包括以下步驟：在單炮地震數據體上逐道、逐樣點判斷方波是否存在；如果該樣點存在方波則進行方波預測與壓制，如果不存在方波，則該樣點值保持不變，接著處理下一個樣點。本發明利用方波能量強，且持續時間長，取值恒定的特點，采取差分結合能量的方式進行判斷，可以較為準確地判斷方波是否存在；在預測方波時，通過均值結合中值技術可以快速的預測方波值，由于是逐樣點處理，去噪比較徹底；該技術克服了常規方法的邊界效應及對有效信號的破壞，保真度較高。

技術領域

本發明屬于石油勘探開發領域，更具體地，涉及一種地震記錄中方波的壓制方法。

背景技術

在地震記錄中，方波通常是由于儀器或檢波器漏碼故障引起，在地震記錄上表現為一個非常大的常數干擾值疊加在有效信號上，形同廢道廢炮。不同道的干擾值相差較大，可能只是影響局部數據，也可能對整道數據都有影響。

常規處理方法往往有去道去炮，內切除，高通濾波，中值濾波等等。去道去炮和內切除都造成了有效信號的喪失和經濟上的損失。而高通濾波因為大值的影響會損害有效信號，且會造成嚴重的邊界效應。中值濾波往往耗時較大，當方波連續多道出現時，去噪結果往往存在較大誤差，且在沒有方波的地方對有效信號有較大影響。

發明內容

針對現有技術的以上缺陷或改進需求，本發明提出了一種地震記錄的方波壓制方法，用于解決去噪對有效信號損害較大的問題。

為實現上述目的，按照本發明的一個方面，提供了一種地震記錄的方波壓制方法，包括：

(1)在單炮地震數據體上逐道、逐樣點進行處理，以各樣點為中心，通過預設長度時窗內的數據進行方波存在與否判斷；

(2)如果該樣點存在方波則進行方波預測與壓制，如果不存在方波，則該樣點值保持不變。

在一些可選的實施方案中，所述通過預設長度時窗內的數據進行方波存在與否判斷，包括：

通過預設長度時窗內的數據，采用差分結合能量的方式進行方波存在與否判斷。

在一些可選的實施方案中，所述通過預設長度時窗內的數據，采用差分結合能量的方式進行方波存在與否判斷，包括：

設地震記錄為s_i，其中，下標i為樣點序號，設置滑動時窗長度參數為2*L+1，識別方波的門檻因子為C，令：L為自然數；

對于第i個樣點，在C_i＞C時，則第i個樣點存在方波，否則第i個樣點不存在方波。

在一些可選的實施方案中，所述進行方波預測與壓制，包括：

通過均值結合中值技術進行方波預測與壓制。

在一些可選的實施方案中，所述通過均值結合中值技術進行方波預測與壓制，包括：

取步驟(1)中的同一個時窗內樣點數據，令

將s_j，(j＝i-L，...，i+L)分為兩個序列A_i，(i＝1，...，M)和B_i，(i＝1，...，N)，M+N＝2*L+1，其中，A_i為序列s_j中大于的樣點值，B_i為序列s_j中小于的樣點值，M，N為自然數；

若M＞N，令IM＝M-L，則在序列A_i中找到第IM小的值為樣點i處方波值，記為V_i；若M≤N，令IN＝N-L，在序列B_i中找到第IN大的值為樣點i處方波值，記為V_i；