[發明專利]基于閾值振幅噪音壓制方法、系統、電子設備及存儲介質

在審

• 申請號: 202111221772.4
• 申請日: 2021-10-20
• 公開（公告）號: CN116010839A
• 公開（公告）日: 2023-04-25
• 發明（設計）人: 張林
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G06F18/24
• 分類號: G06F18/24;G06F18/10
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

本發明提供一種基于閾值振幅噪音壓制方法、系統、電子設備及存儲介質，包括從原始地震記錄中選擇某一道地震道，計算該地震道的中值絕對差；基于中值絕對差和振幅閾值因子，計算該地震道的閾值振幅；判斷該地震道的振幅值是否達到閾值振幅；根據閾值振幅對該地震道做衰減處理，并輸出衰減后的地震道。本發明的技術方案的有益效果在于，采用中值振幅和給定的振幅閾值因子定義閾值振幅，振幅閾值因子隨時間和空間而變化，根據閾值振幅對地震道的振幅值進行衰減處理以達到壓制異常噪音的目的，能夠保護有效信號，可以提高地震數據資料的信噪比，對探區野外采集的單炮中的異常噪音干擾具有非常好的應用效果。

技術領域

本發明屬于地震勘探資料處理方法技術領域，更具體地，涉及一種基于閾值振幅噪音壓制方法、系統、電子設備及存儲介質。

背景技術

在地震數據資料處理中，異常噪音是一種常見的噪音，它的形態表現為沒有一定的傳播方向和視速度，且能量極強，如果不在疊前壓制的話，對后續的速度建模和偏移成像都非常不利，因此異常噪音的壓制非常關鍵。

如何有效的去除異常噪音，一直是資料處理中的重點問題。異常噪音壓制的方法一般來說，都是將數據轉換到頻率域，然后應用空間中值濾波。給定特定閾值后，對中值振幅中提取的頻帶范圍，乘以一個特定的標定算子或者使用周圍道的插值頻帶替代原有振幅值。當采樣點的中值振幅超出頻帶內的中值振幅時，那么這一頻帶就會被認為是異常道。通常的做法都是采用閾值因子，閾值因子的類型可以是中值或者中值絕對差。計算中值振幅間的絕對差和計算中值寬度范圍內每一道的振幅來獲得中值絕對差。當數據空間變化很大時，中值絕對差也需要自動的調整；當數據不那么復雜時，中值絕對差也需要相應的調整。

因此期待一種基于閾值振幅噪音壓制方法，旨在對地震數據進行異常噪音壓制，消除異常噪音干擾。

發明內容

本發明提供一種基于閾值振幅噪音壓制方法、系統、電子設備及存儲介質，采用中值振幅和給定的振幅閾值因子定義閾值振幅，振幅閾值因子隨時間和空間而變化，根據閾值振幅對地震道的振幅值進行衰減處理以達到壓制異常噪音的目的，能夠保護有效信號，可以提高地震數據資料的信噪比。

為實現上述目的，本發明提供一種基于閾值振幅噪音壓制方法，包括：

步驟S1，從原始地震記錄中選擇某一道地震道，計算該地震道的中值絕對差；

步驟S2，基于所述中值絕對差和振幅閾值因子，計算該地震道的閾值振幅；

步驟S3，判斷該地震道的振幅值是否達到所述閾值振幅，若是，則轉入步驟S4進行處理，若否，則不做處理；

步驟S4，根據所述閾值振幅對該地震道做衰減處理，并輸出衰減后的地震道。

優選的，所述步驟S1包括：基于選擇的某一道地震道的時間t和頻帶f，先計算某一道地震道的中值振幅間的絕對差，然后根據中值振幅間的絕對差對中值寬度內的每一道地震道計算振幅，得到某一道地震道在時間t頻帶f的中值絕對差；

所述中值絕對差由下述公式(1)計算得到：

P_ftk＝med_k[|R_ftk-med_k(R_ftk)|]???(1)；

其中，P_ftk為某一道地震道在時間t頻帶f的中值絕對差，med_k為對某一道地震道在時間t頻帶f的中值寬度內的所有道做中值，R_ftk為某一道地震道在時間t頻帶f的振幅能量。

優選的，當某一道地震道的振幅值在閾值振幅范圍內時，所述振幅閾值因子由下述公式(2)計算得到：

其中，U_t為所述振幅閾值因子，x為滑動因子。