[发明专利]基于块匹配算法和奇异值分解的地震资料噪声压制方法

权利要求书

1.一种基于块匹配算法和奇异值分解的地震资料噪声压制方法，其特征在于，包括以下步骤：

步骤一、相似性分组：将整个地震资料切割成过完备的子块，遍历地选择参考子块，计算参考子块与一定邻域内的其他子块的欧氏距离，判断两子块是否相似，并把相似二维子块按照三维数据的方式存放起来成为相似组；

步骤二、三维奇异值滤波去噪：对相似组中每个子块先进行二维离散余弦变换，选用频域二维奇异值分解滤波，再将相似块向量化，进行频域一维奇异值分解，并进行奇异值保留去噪；二维离散余弦反变换后，将过完备子块加权聚合成原始尺寸地震资料，得到基础估计；

步骤三、计算维纳收缩系数：将基础估计重新分块分组，并且通过新的相似组计算变换域的维纳收缩系数；

步骤四、维纳滤波：按照新的相似分组，对原始信号进行三维离散余弦变换，采用维纳滤波实现噪声压制，三维离散余弦反变换并聚合重构得到最终噪声压制后的地震资料。

2.如权利要求1所述的一种基于块匹配算法和奇异值分解的地震资料噪声压制方法，其特征在于，所述步骤一相似性分组包括以下具体步骤：

1)预处理：先将整个地震资料按照一定的步进长度过完备的分割成N₁×N₁大小的时-空域子块Y_x；然后对各个子块的二维频谱进行阈值r′处理，将大于阈值的频谱值保留下来，小于阈值的频谱值置0；

2)计算块间距离：将预处理后的子块频谱r′(T(Y_x))作为计算块间距离的元素，选定一个参考块Y_R，选定搜索窗大小为N_h×N_h，通过2-范数的计算参考块与其相应搜索窗内其他子块的块间距离：

d(YR,Yx)=||r′(T(YR))-r′(T(Yx))||22(N1)2]]>

3)堆叠成三维数组：当两个子块的距离d(Y_R,Y_x)小于阈值τ时，判定两块相似，然后将所有与同一参考块相似的二维子块，连同参考块一起存放成一个三维数组，即相似组Z_R。