[发明专利]微震记录的去噪方法及装置

说明书

本发明公开了一种微震记录的去噪方法及装置。其中，该方法包括：采集地震道的微震记录，在微震记录中的微震数据中选择参照块，并在参照块的预定范围进行搜索，得到数据块，其中，参照块是按照预定规格从微震数据中选取的微震数据；将数据块和参照块进行整合，得到三维矩阵；确定三维矩阵对应的二维矩阵，并对二维矩阵进行奇异值分解，得到奇异值和分解矩阵，其中，奇异值为多个；对分解矩阵进行重构，并根据重构后的分解矩阵得到滤波因子；根据滤波因子对微震记录中的微震数据进行去噪。本发明解决了相关技术利用传统去噪方式进行微震资料去噪不理想，无法达到突出有效信号的技术问题。

技术领域

本发明涉及微地震监测领域，具体而言，涉及一种微震记录的去噪方法及装置。

背景技术

在微震资料处理中，微震资料去噪是整个微震资料处理流程中直观重要的步骤。但对于地面接收的微震资料而言，有效信号的能量较弱，完全淹没在噪声之中，微震剖面的信噪比极低。而传统的去噪方法往往要求资料具有较高的信噪比，致使通过传统的去噪方法很难得到理想的去噪效果。

针对上述相关技术利用传统去噪方式进行微震资料去噪不理想，无法达到突出有效信号的问题，目前尚未提出有效的解决方案。

发明内容

本发明实施例提供了一种微震记录的去噪方法及装置，以至少解决相关技术利用传统去噪方式进行微震资料去噪不理想，无法达到突出有效信号的技术问题。

根据本发明实施例的一个方面，提供了一种微震记录的去噪方法，包括：采集地震道的微震记录，在所述微震记录中的微震数据中选择参照块，并在所述参照块的预定范围进行搜索，得到数据块，其中，所述参照块是按照预定规格从所述微震数据中选取的微震数据；将所述数据块和所述参照块进行整合，得到三维矩阵；确定所述三维矩阵对应的二维矩阵，并对所述二维矩阵进行奇异值分解，得到奇异值和分解矩阵，其中，所述奇异值为多个；对所述分解矩阵进行重构，并根据重构后的分解矩阵得到滤波因子；根据所述滤波因子对微震记录中的微震数据进行去噪。

可选地，将所述数据块和所述参照块进行整合，得到三维矩阵包括：基于所述参照块确定所述数据块的选取范围；在所述选取范围在所述参照块的预定范围内进行搜索，得到差异度小于预定阈值的数据块；对所述数据块进行整合，同时对所述参照块进行整合，得到所述三维矩阵。

可选地，对所述二维矩阵进行奇异值分解，得到奇异值包括：通过第一公式对所述二维矩阵进行奇异值分解，其中，所述第一公式为：A表示所述二维矩阵，U表示第一正交矩阵，由所述二维矩阵A和所述二维矩阵的转置矩阵A^T的乘积得到的特征向量构成；V表示第二正交矩阵，由所述二维矩阵的转置矩阵A^T和所述二维矩阵A的乘积得到的特征向量构成；∑是由所述二维矩阵A和所述二维矩阵的转置矩阵A^T的乘积确定的奇异值，或所述二维矩阵的转置矩阵A^T和所述二维矩阵A的乘积确定奇异值σ_i按照递减顺序组成的对角矩阵，σ_i表示第i个奇异值，i为正整数。

可选地，所述对角矩阵∑为Σ＝diag(σ₁,σ₂,…σ_r)，其中，r表示所述二维矩阵A的秩，ρ表示奇异值。

可选地，所述二维矩阵A和所述二维矩阵的转置矩阵A^T的乘积确定的奇异值或所述二维矩阵的转置矩阵A^T和所述二维矩阵A的乘积确定奇异值σ_i与特征值λ_i的对应关系为：

可选地，对所述分解矩阵进行重构包括：选取分解得到的奇异值中位于中间位置的部分奇异值；利用所述部分奇异值对所述分解矩阵进行重构。

可选地，根据重构后的分解矩阵得到滤波因子包括：根据所述二维矩阵和重构得到的二维矩阵确定所述滤波因子。