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

权利要求书

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中任意一项所述的微震记录的去噪方法。