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

摘要

本发明公开了一种基于块匹配算法和奇异值分解的地震资料噪声压制方法，解决现有方法由于没有考虑到地震资料同相轴固有的局部和非局部相似性，造成的噪声压制不充分和有效信息衰减的问题，该方法包括相似性分组、三维变换奇异值滤波去噪、根据基础估计重新分组并计算维纳收缩系数和维纳滤波四个步骤。该方法在理论上采用块匹配算法，考虑到地震资料同相轴的局部和非局部相似性的固有特点，采用二维奇异值分解和一维奇异值分解结合，充分利用地震资料的结构信息，实现地震资料的稀疏表示，能够较大程度地减小有效信号衰减，并提高噪声压制效果。

主权项

一种基于块匹配算法和奇异值分解的地震资料噪声压制方法，其特征在于，包括以下步骤：步骤一、相似性分组：将整个地震资料切割成过完备的子块，遍历地选择参考子块，计算参考子块与一定邻域内的其他子块的欧氏距离，判断两子块是否相似，并把相似二维子块按照三维数据的方式存放起来成为相似组；步骤二、三维奇异值滤波去噪：对相似组中每个子块先进行二维离散余弦变换，选用频域二维奇异值分解滤波，再将相似块向量化，进行频域一维奇异值分解，并进行奇异值保留去噪；二维离散余弦反变换后，将过完备子块加权聚合成原始尺寸地震资料，得到基础估计；步骤三、计算维纳收缩系数：将基础估计重新分块分组，并且通过新的相似组计算变换域的维纳收缩系数；步骤四、维纳滤波：按照新的相似分组，对原始信号进行三维离散余弦变换，采用维纳滤波实现噪声压制，三维离散余弦反变换并聚合重构得到最终噪声压制后的地震资料。