[发明专利]三维傅里叶变换保幅叠前道内插方法

说明书

技术领域

本发明涉及地震资料处理，特别是涉及到一种三维傅里叶变换保幅叠前道内插方法。

背景技术

就目前来讲，常用的地震道插值方法主要有Sinc插值方法、                                                域插值方法、域预测误差滤波插值方法。其中Sinc地震道插值方法虽然在进行地震道插值时速度快，易于实现，但是无法正确内插具有空间假频的地震道。Spitz的域地震道插值方法和Claerbout的域预测误差滤波插值方法虽然可以正确内插具有空间假频的地震道，但是这些方法计算量特别大，使其在实际数据中的实用性大大降低。基于这些问题，Porsani(1999)提出了半步长滤波技术实现道内插，虽然效果较高，但其反射信号同相轴必须是线性的且各道振幅不变的假设条件非常苛刻；相比而言，域地震道插值方法得到了更为广泛研究，Gulunay和Chambers（1996）提出了抗假频域地震道插值方法（UFKI），但抗假频域地震道插值方法（UFKI）的插值系数L不能为任意整数，即两相邻地震道之间不能被插入任意整数的插值道，广义域地震道插值方法GFKI是UFKI差值法的改进型，虽然其插值系数可以为任意整数，但同样要求插值的地震道数据必须是等间隔采样的，因此这两种抗假频域地震道插值方法均不能很好应用在三维地震数据中。

随着油气勘探开发的不断深入，对地下构造研究的精度要求越来越高。在地震勘探中，由于工作量的限制及施工条件的影响或是地震资料采集的孔径效应以及检波器可能出现的故障等因素的影响，在采集到的数据体中会出现某些测线方向上地震记录道间距过大、野值道或道缺失情况，从而造成了空间采样率的不足，严重影响了成像质量。对一些要求作精细地质解释的地区，只需对现有的地震资料作保幅地震道插值，然后进行偏移处理，而不必重新进行施工采样，这样不但可以节约大量人力、物力、财力，而且还可以大大缩短生产周期。保幅法地震道插值技术能保持地震道的相对能量关系，不产生附加的空间假频，使地球物理信息更加真实地反映地下地质体的地球物理特征，有利于进行构造解释和地震地层学的研究。因此，对保幅内插的研究极具意义。为此我们发明了一种新的三维傅里叶变换保幅叠前道内插方法。

发明内容

本发明的目的是提供一种地球物理意义明显，效率高的三维傅里叶变换保幅叠前道内插方法。    

本发明的目的可通过如下技术措施来实现：三维傅里叶变换保幅叠前道内插方法，该三维傅里叶变换保幅叠前道内插方法包括：对三维地震数据进行道头及数据分析，判断内插类型；对该三维地震数据进行三维傅里叶变换，得到F-Kx-Ky三维频率波数谱；在Kx与Ky方向对该三维频率波数谱充零拓展；对拓展后的该三维频率波数谱进行三维傅里叶逆变换，得到t-x-y域数据；对得到的该t-x-y域数据进行振幅增益处理；以及输出内插后的三维地震数据。益处理；以及输出内插后的三维地震数据。

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

该三维傅里叶变换保幅叠前道内插方法还包括在判断内插类型的步骤前，输入该三维地震数据。

将正演数据作为输入的该地震三维数据。

将采集的野外实际数据作为输入的该三维地震数据。

该三维地震数据为SEG-Y格式的地震数据。

该三维地震数据包括数据名称、道数、采样点数和采样率。

在进行振幅增益处理的步骤中，对得到的该t-x-y域数据按照内插前的振幅均值进行增益处理。

本发明中的三维傅里叶变换保幅叠前道内插方法立足三维地震数据体，通过三维傅里叶变换方法（3DFT），在时间t方向和空间x方向、y方向进行三维傅里叶变换，形成三维的F-Kx-Ky谱。通过对F-Kx-Ky谱拓展处理和三维逆变换，完成三维数据的内插。本发明中的三维傅里叶变换保幅叠前道内插方法可实现内插炮点、内插缺失道，内插后的数据保幅性好。本发明实现了保幅内插的需求，地球物理意义明显，效率高。

附图说明

图1为本发明的实施例中三维傅里叶变换保幅叠前道内插方法的流程图；

图2为本发明的实施例1中存在缺失道的正演炮集的示意图；

图3是本发明的实施例1中三维傅里叶变换后的频率波数谱的示意图；

图4为本发明的实施例1中充零后的频率波数谱的示意图；

图5为本发明的实施例1中内插后的炮集的示意图；

图6为本发明的实施例2中原始正演炮集记录的示意图；

图7为本发明的实施例2中内插出的炮集记录的示意图；