[發明專利]微震記錄的去噪方法及裝置

在審

• 申請號: 201811545159.6
• 申請日: 2018-12-17
• 公開（公告）號: CN111324598A
• 公開（公告）日: 2020-06-23
• 發明（設計）人: 謝斌;潘勇;蘆志偉;段勝男;王寧博;潘樹林;李士建;侯萬勇;闕軍仁;陳浩
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G06F16/215
• 分類號: G06F16/215
• 代理公司: 北京康信知識產權代理有限責任公司 11240
• 代理人: 趙囡囡;董文倩

1.一種微震記錄的去噪方法，其特征在于，包括：

采集地震道的微震記錄，在所述微震記錄中的微震數據中選擇參照塊，并在所述參照塊的預定范圍進行搜索，得到數據塊，其中，所述參照塊是按照預定規格從所述微震數據中選取的微震數據；

將所述數據塊和所述參照塊進行整合，得到三維矩陣；

確定所述三維矩陣對應的二維矩陣，并對所述二維矩陣進行奇異值分解，得到奇異值和分解矩陣，其中，所述奇異值為多個；

對所述分解矩陣進行重構，并根據重構后的分解矩陣得到濾波因子；

根據所述濾波因子對微震記錄中的微震數據進行去噪。

2.根據權利要求1所述的方法，其特征在于，將所述數據塊和所述參照塊進行整合，得到三維矩陣包括：

基于所述參照塊確定所述數據塊的選取范圍；

在所述選取范圍在所述參照塊的預定范圍內進行搜索，得到差異度小于預定閾值的數據塊；

對所述數據塊進行整合，同時對所述參照塊進行整合，得到所述三維矩陣。

3.根據權利要求1所述的方法，其特征在于，對所述二維矩陣進行奇異值分解，得到奇異值包括：

通過第一公式對所述二維矩陣進行奇異值分解，其中，所述第一公式為：A表示所述二維矩陣，U表示第一正交矩陣，由所述二維矩陣A和所述二維矩陣的轉置矩陣A^T的乘積得到的特征向量構成；V表示第二正交矩陣，由所述二維矩陣的轉置矩陣A^T和所述二維矩陣A的乘積得到的特征向量構成；∑是由所述二維矩陣A和所述二維矩陣的轉置矩陣A^T的乘積確定的奇異值，或所述二維矩陣的轉置矩陣A^T和所述二維矩陣A的乘積確定奇異值σ_i按照遞減順序組成的對角矩陣，σ_i表示第i個奇異值，i為正整數。

4.根據權利要求3所述的方法，其特征在于，所述對角矩陣∑為Σ＝diag(σ₁,σ₂,…σ_r)，其中，r表示所述二維矩陣A的秩，ρ表示奇異值。

5.根據權利要求4所述的方法，其特征在于，所述二維矩陣A和所述二維矩陣的轉置矩陣A^T的乘積確定的奇異值或所述二維矩陣的轉置矩陣A^T和所述二維矩陣A的乘積確定奇異值σ_i與特征值λ_i的對應關系為：

6.根據權利要求5所述的方法，其特征在于，對所述分解矩陣進行重構包括：

選取分解得到的奇異值中位于中間位置的部分奇異值；

利用所述部分奇異值對所述分解矩陣進行重構。

7.根據權利要求6所述的方法，其特征在于，根據重構后的分解矩陣得到濾波因子包括：根據所述二維矩陣和重構得到的二維矩陣確定所述濾波因子。

8.一種微震記錄的去噪裝置，其特征在于，包括：

選擇單元，用于采集地震道的微震記錄，在所述微震記錄中的微震數據中選擇參照塊，并在所述參照塊的預定范圍進行搜索，得到數據塊，其中，所述參照塊是按照預定規格從所述微震數據中選取的微震數據；

整合單元，用于將所述數據塊和所述參照塊進行整合，得到三維矩陣；

分解單元，用于確定所述三維矩陣對應的二維矩陣，并對所述二維矩陣進行奇異值分解，得到奇異值和分解矩陣，其中，所述奇異值為多個；

重構單元，用于對所述分解矩陣進行重構，并根據重構后的分解矩陣得到濾波因子；

去噪單元，用于根據所述濾波因子對微震記錄中的微震數據進行去噪。

9.一種存儲介質，其特征在于，所述存儲介質包括存儲的程序，其中，所述程序執行權利要求1至7中任意一項所述的微震記錄的去噪方法。

10.一種處理器，其特征在于，所述處理器用于運行程序，其中，所述程序運行時執行權利要求1至7中任意一項所述的微震記錄的去噪方法。