[發明專利]混疊噪聲壓制方法及裝置

在審

• 申請號: 202010170533.X
• 申請日: 2020-03-12
• 公開（公告）號: CN113391354A
• 公開（公告）日: 2021-09-14
• 發明（設計）人: 董烈乾;張慕剛;汪長輝;駱飛;周大同
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 周曉飛;許曼

1.一種混疊噪聲壓制方法，其特征在于，包括：

將混疊的共炮點域數據轉換為共偏移距域混疊數據；

對共偏移距域混疊數據進行混疊噪聲壓制；

若壓制混疊噪聲后的共偏移距域混疊數據的信噪比小于預設值，則對所述壓制混疊噪聲后的共偏移距域混疊數據進行混疊計算，得到混疊數據，將所述混疊的共炮點域數據更新為所述混疊數據；

重復執行以上步驟，直到壓制混疊噪聲后的共偏移距域混疊數據的信噪比大于等于預設值；

將信噪比大于等于預設值的壓制混疊噪聲后的共偏移距域混疊數據，轉換為壓制混疊噪聲后的共炮點域數據。

2.如權利要求1所述的方法，其特征在于，對共偏移距域混疊數據進行混疊噪聲壓制，包括：

對共偏移距域混疊數據進行奇異值分解，得到共偏移距域混疊數據的奇異值向量；

從共偏移距域混疊數據中選取一僅包括混疊噪聲的矩形窗，對矩形窗中的數據進行奇異值分解，得到混疊噪聲的奇異值向量；

根據共偏移距域混疊數據的奇異值向量和混疊噪聲的奇異值向量，確定壓制混疊噪聲后的共偏移距域混疊數據。

3.如權利要求2所述的方法，其特征在于，按如下公式，對共偏移距域混疊數據進行奇異值分解，得到共偏移距域混疊數據的奇異值向量：

d_offset＝U_bl V_blW^T_bl

其中，V_bl為共偏移距域混疊數據的奇異值向量；d_offset為共偏移距域混疊數據；U_bl為共偏移距域混疊數據的左特征向量；W^T_bl為共偏移距域混疊數據的右特征向量。

4.如權利要求2所述的方法，其特征在于，從共偏移距域混疊數據中選取一僅包括混疊噪聲的矩形窗，包括：

從共偏移距域混疊數據中，選取距離激發點最近的共偏移距域混疊數據，確定有效波最先到達的時間；

根據有效波最先到達的時間，從距離激發點最近的共偏移距域混疊數據中，選取一僅包括混疊噪聲的矩形窗。

5.如權利要求2所述的方法，其特征在于，按如下公式，對選取的矩形窗中的混疊數據進行奇異值分解，得到混疊噪聲的奇異值向量：

d_noise＝U_noiseV_noiseW^T_noise

其中，V_noise為混疊噪聲的奇異值向量；d_noise為選取的僅含有混疊噪聲的數據；U_noise為混疊噪聲的左特征向量；W^T_noise為混疊噪聲的右特征向量。

6.如權利要求2所述的方法，其特征在于，根據共偏移距域混疊數據的奇異值向量和混疊噪聲的奇異值向量，確定壓制混疊噪聲后的共偏移距域混疊數據，包括：

取混疊噪聲的奇異值向量的最大值為混疊噪聲的約束值；

采用硬閾值函數，根據混疊噪聲的約束值，對共偏移距域混疊數據的奇異值向量進行更新計算，得到壓制混疊噪聲后的共偏移距域混疊數據的奇異值向量；

根據壓制混疊噪聲后的共偏移距域混疊數據的奇異值向量，確定壓制混疊噪聲后的共偏移距域混疊數據。

7.如權利要求6所述的方法，其特征在于，按如下公式，取混疊噪聲的奇異值向量的最大值為混疊噪聲的約束值：

λ＝max(V_noise)

其中，λ為混疊噪聲的約束值；V_noise為混疊噪聲的奇異值向量。