[发明专利]基于波动方程延拓的近地表噪音压制方法

说明书

本发明提供一种基于波动方程延拓的近地表噪音压制方法，该基于波动方程延拓的近地表噪音压制方法包括：步骤１，输入近地表速度厚度模型；步骤２，将数据进行选排，分选为共炮，共检道集；步骤３，共炮道集对检波点进行延拓，延拓到低速带底界面；步骤４，共炮道集对检波点进行延拓，向上延拓到地表；步骤５，将完成了炮点和检波点延拓的数据，进行重新选排，重新抽取回共炮道集，完成输出。该基于波动方程延拓的近地表噪音压制方法几乎不损失有效信号就可以将与近地表相关噪音进行压制，从原理上就可以保护有效波，去除噪音信号。

技术领域

本发明涉及油气勘探地震资料处理技术领域，特别是涉及到一种基于波动方程延拓的近地表噪音压制方法。

背景技术

地震资料中的噪音主要与近地表相关，如强面波干扰、强折射波、声波以及近地表绕射源产生的噪音，这些噪音的存在严重影响了地震资料的信噪比。针对不同的噪音特征，目前已经提出了许多种噪音压制的方法，如低F-K域，Radon变换，曲波变换，小波变换等不同域进行去噪的处理。这些方法都具有各自的特点并取得了明显的效果，但往往只是侧重考虑某一方面的特性或依据某种假设条件，都有其局限性。F-K滤波根据有效信号和噪音的视速度差异进行滤波，在很多实际资料中面波和有效信号视速度差异不够明显，造成去噪效果差；Radon变换和F-K滤波器相似，也是依据经过变换后在Randon域噪音和有效信号会形成不同的聚焦能量，但是实际噪音往往无法和有效信号彻底分离，造成去噪效果并不理想；曲波变换和小波变换方法的问题和前述方法类似，均是通过有效信息和噪音的某种差异进行信噪分离，在不能精确分离时，必然伤及有效信号。

由于地震资料信噪比对处理解释的影响巨大，而目前常用的去噪方法存在诸多问题。随着油气勘探工作的不断深化，地震资料处理中对任何有效信号的伤害必然对后续的处理解释带来无法弥补的损失。

地震记录中噪音种类繁多，强折射、面波、声波等干扰各具特点，由于这些干扰的存在严重影响了地震的信噪比、分辨率，如果叠前不对其进行有效地压制，最终必然会影响到整个资料处理的质量。在压制噪音的同时，如何避免对有效信号的损害，一直是地震资料处理的重要内容。为此我们发明了一种新的基于波动方程延拓的近地表噪音压制方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种最大限度保护有效信息，同时压制多种噪音干扰的基于波动方程延拓的近地表噪音压制方法。

本发明的目的可通过如下技术措施来实现：基于波动方程延拓的近地表噪音压制方法，该基于波动方程延拓的近地表噪音压制方法包括：步骤1，输入近地表速度厚度模型；步骤2，将数据进行选排，分选为共炮，共检道集；步骤3，共炮道集对检波点进行延拓，延拓到低速带底界面；步骤4，共炮道集对检波点进行延拓，向上延拓到地表；步骤5，将完成了炮点和检波点延拓的数据，进行重新选排，重新抽取回共炮道集，完成输出。

本发明的目的还可通过如下技术措施来实现：

在步骤1中，输入的模型为离散网格数据，每个网格点的数值为当前位置纵波速度。

在步骤2中，将数据分别选排为共炮道集以进行检波点延拓和共检波点道集以进行共炮点延拓，进行延拓时采用的基本运算公式为：

式中，u表示在某一深度位置记录波场的傅里叶变换，v表示纵波在介质中的传播速度，ω表示圆频率，k_x表示水平方向上的圆波数，ΔZ表示深度域延拓步长，Z表示当前深度位置，公式中指数位置取正为向上延拓，取负值为向下延拓。

在步骤3中，选择有限差分方法在炮集和检波点道集记录中进行波场记录的向下延拓，将检波点和炮点延拓的低速带底界面，实现波场向下延拓过程，在延拓过程，仅保留正时间的地震数据，即延拓过程中对当前延拓面之上的信息进行清除。