[发明专利]一种压制叠前地震资料线性干扰方法、存储介质及设备

说明书

本发明公开了一种压制叠前地震资料线性干扰方法、存储介质及设备，对原始地震炮集数据进行处理得到新数据；对新数据使用离散傅里叶变换后排成Hankel矩阵，对Hankel矩阵进行奇异值分解降秩重构，对重构数据使用一维逆离散傅里叶变换得到线性干扰波的估计，将估计得到的线性干扰波数据反拉平后从原始炮集数据中减去，即完成该视倾角方向线性干扰的压制，重复以上步骤直到原始炮集数据中所有线性干扰波均被处理。本发明能够在有效信号与线性干扰信号空间混叠时，有效压制线性干扰且避免损伤有效信号，使原本被线性干扰波覆盖的信号得到更为清晰直观的显示，压制线性干扰后的炮集数据可用于后期速度分析、叠加处理、成像、地质结构分析。

技术领域

本发明属于地震勘探数据处理技术领域，具体涉及一种压制叠前地震资料线性干扰方法、存储介质及设备。

背景技术

在人工地震勘探的地震资料采集和处理过程中会遇到一些规则干扰波，其中在地震剖面上呈明显线性结构的规则干扰波称为线性干扰波，其一直是影响地震资料信噪比的关键因素之一。在线性干扰波发育的地区采集的原始地震资料上有效地震反射波信息受到较强的干涉作用，甚至被完全淹没，同时激发地震波能量的下传也被屏蔽，造成深部目的层有效信号能量弱，使后续叠加处理、偏移成像、地质解释极其困难。线性干扰波在从我国地表覆盖较薄的西部地区以及山地和沙丘区域采集的地震资料中普遍存在。例如在滕县煤田微山湖某矿区地震勘探项目中，采集的单炮记录上存在较强的多次反射折射波，对主要煤层反射波有较强的干涉作用，并导致煤层反射波能量变弱，信噪比低，以及陆上地震勘探采集的地震资料中常见沙丘鸣震线性干扰波，其往往淹没反射波场，需要采取一定的技术方法进行有效压制。

传统压制线性干扰的方法有频率波数域滤波、τ-p变换、奇异值分解滤波法和速度滤波，但这些方法均取得了一定的压制效果，但当线性干扰波和有效信号在空间上有混叠或数据采集间距不规则时，效果有限。仍需要进一步改进现有的技术方法在这些情况下取得更好的效果。

现有技术：

SVD滤波方法是基于特征值对应的特征向量进行不同的信号重建达到分离波场目的的一种方法，利用有效信号的相关性，识别噪声对应的特征值；然后利用特征值对应的特征向量重建地震信号

现有技术的缺点：

1、需要手动定位需要滤波区域，用以重构噪声的奇异值个数需要反复试验获取，不准确会造成去噪效果差或者损伤有效信号。

2、当有效信号与线性干扰波在空间上混叠时，易损伤有效信号。

发明内容

本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种压制叠前地震资料线性干扰方法、存储介质及设备。本发明将原始地震炮集数据主要看作两部分：线性干扰波和其它信号(其它信号中主要成分为有效信号时，近似为有效信号)。首先读取原始地震炮集资料，然后通过扫描分析获取原始地震炮集数据中任一线性干扰波视倾角方向并沿此方向拉平线性干扰波，通过随机抽样或者随机重排拉平后的数据得到新数据。对新数据中的每道数据沿时间方向使用一维离散傅里叶变换(DFT)，将每个空间方向上的频点排成Hankel矩阵，对Hankel矩阵进行奇异值分解，保留前三个范围内的奇异值重构，对重构的数据每个空间位置做一维逆离散傅里叶变换(IDFT)，得到线性干扰波的估计，将估计得到的线性干扰波数据反拉平后从原始炮集数据中减去，即完成该视倾角方向线性干扰的压制。重复以上除读取原始地震炮集资料外步骤直到原始炮集数据中所有需要压制的不同倾角方向的线性干扰波均被处理，得到压制线性干扰后的炮集数据用于后期速度分析、叠加处理、成像、地质结构分析。

本发明采用以下技术方案：

一种压制叠前地震资料线性干扰方法，包括以下步骤：

S1、读取原始地震炮集资料；