[發明專利]基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法

在審

• 申請號: 202211437162.2
• 申請日: 2022-11-17
• 公開（公告）號: CN115840249A
• 公開（公告）日: 2023-03-24
• 發明（設計）人: 段偉偉;張倩慧;彭仁艷;楊晨瑩;陳良;蘇燕;顏曉霞;趙思思;付婷
• 申請（專利權）人: 中石化石油工程技術服務有限公司;中石化石油工程地球物理有限公司;中石化石油工程地球物理有限公司科技研發中心
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30;G01V1/36
• 代理公司: 南京天華專利代理有限責任公司 32218
• 代理人: 許軻;徐冬濤

1.一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于包括如下步驟：

步驟一、根據“黑三角”強能量干擾分布特征對地震道給定一個隨機編號進行重排；

步驟二、采用時變、空變及多域處理等現有的處理方式，利用矢量中值濾波進行三角區強能量干擾衰減；

步驟三、利用高分辨率Radon變換，對干擾波和有效波進行分離，最終實現“黑三角”區的信噪分離。

2.根據權利要求1所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，所述步驟一和二中，首先將地震數據分為不同的頻段，在不同的頻段范圍內識別地震記錄中存在的強能量干擾，確定噪聲出現的空間位置，再定義相應的門檻值和衰減系數，然后采用時變、空變及多域的方式進行壓制，達到壓制高能干擾以及實現高保真去噪。

3.根據權利要求2所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，強能量干擾在振幅、頻率特征方面，呈隨機發育特征。

4.根據權利要求1所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，所述步驟三中，利用Radon變換進行波場分離時，通過在最小二乘拋物線拉冬變換的算子中引入非線性規則化對角矩陣W實現將拉冬域數據約束至一些離散的(τ，q)點，即：

(LL^H+W^HW)M＝LD (1)

其中，非線性對角矩陣W的物理意義理解為曲率參數q的權重向量，L為Radon變換算子，L^H表示L的共軛轉置算子，M是頻率域地震數據，D是Radon域數據和頻率空間域數據相互關系的矩陣形式。

5.根據權利要求4所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，采用廣義最小二乘原理，遞歸求取權重的方法，將Radon域數據的能量集中于其能量最強的頻帶內，實現對有效波與干擾波的高效分離。

6.根據權利要求5所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，具體的變換過程表示為：

M(ω_k)＝WL(L^HWL+εI)^-1D(ω_k) (2)

其中，M(ω_k)為頻率為k時的頻率域地震數據，I為單位矩陣，ε為阻尼因子，D(ω_k)是Radon域數據和頻率空間域數據在頻率為k時的相互關系的矩陣形式。

7.根據權利要求6所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，利用低頻帶求出的變換結果作為對高頻帶權重的約束因子。

8.根據權利要求7所述的一種基于隨機道矢量中值濾波及高分辨率Radon變換的黑三角噪音壓制方法，其特征在于，采用的權重求取方法表示為：

W_ii(ω_k)＝||M_i(ω_k-1)||，i＝1,…,N_q (3)

其中，W_ii(ω_k)為頻率為k時的權重。