[发明专利]基于频率空间域波场延拓的地震资料缺口补偿方法

说明书

技术领域

本发明涉及勘探地球物理地震数据处理技术领域，是一种利用时间-空间域内的波场延拓关系实现地震资料缺口补偿的方法。

背景技术

在地震数据处理中，地震道缺失和空间采样不足是很常见的，其表现为死道或者由于含有强烈的噪声而在预处理过程中被剔除的道记录。这样，基于多道处理算法的处理过程将受到缺失道的影响，其中受影响最严重的包括波动方程偏移，基于波动方程的抑制表面多次波的方法，与地表有关的多次波的消除和谱估计等处理过程。

地震资料补偿方法可以在不重新进行地震资料采集的情况下，补偿缺失的地震数据，保障后续处理的顺利进行。目前常用的地震资料补偿方法主要有Sinc插值方法、f-x域插值方法、t-x域预测误差滤波插值方法和基于Radon变换法的地震道插值方法。其中Sinc插值方法虽然在进行地震道插值时速度快，易于实现，但是无法正确内插具有空间假频的地震道。Spitz的f-x域地震道插值方法和Claerbout的t-x域预测误差滤波插值方法虽然可以正确内插具有空间假频的地震道，但是这些方法计算量特别大，使其在实际数据中的实用性大大降低。基于Radon变换法的地震道插值方法可以用于地震资料的补偿，但该方法本身不是保幅的，地震数据经过Radon正变换和Radon反变换之后，数据会有微小的改变。

此外多次波是地下反射层的多次反射，也蕴含了地层结构信息，在极其复杂的地质条件下，甚至包含了一次波不具有的地下信息，用适当的方法加以利用可以获得更加丰富的地质构造信息，但这种方法通常只适用于海上资料，大部分陆上地震资料的多次波并不发育。

发明内容

本发明的目的是针对现有技术的缺陷，提出一种为提高地震成像质量，解决三维地震勘探中由于仪器故障或地表条件影响造成地震道缺失，导致后续处理困难的基于频率空间域波场延拓的地震资料缺口补偿方法。

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

(1)对三维地震数据进行三维傅里叶变换，将(x,y,t)域地震数据转换到k_x-k_y-f域；

(2)利用地震信号的空间相干性，根据时间-空间域内的波场延拓关系，对非缺失道地震数据的波场进行延拓，得到缺失道的波场信息；

(3)进行三维傅里叶反变换，得到(x,y,t)域的补偿数据。

本发明的技术方案进一步优化为：

(1)对输入的三维地震数据进行三维傅里叶变换，将(x,y,t)域地震数据转换到k_x-k_y-f域，变换公式为；