[发明专利]通过防混淆、防泄漏傅立叶变换内插地震数据的方法

权利要求书

1.一种用于在地震数据中内插轨迹的方法，包括：

通过获得时间-空间域中的输入地震数据；

将时间傅立叶变换应用于输入地震数据，从而产生变换后的地震 数据；

确定变换后的地震数据中的哪些频率分量相应于非混淆频率以 及变换后的地震数据中的哪些频率分量相应于混淆频率；

将第一防泄漏傅立叶变换方法应用于变换后的地震数据的非混 淆频率中的每一个；

将第二防泄漏傅立叶变换方法应用于变换后的地震数据的混淆 频率的每一个来产生估计频率-波数谱，其中第二防泄漏傅立叶变换 方法应用从非混淆频率外推到混淆频率的绝对频率-波数谱来加权混 淆频率的频率-波数分量；以及

将逆时间和空间傅立叶变换应用于估计频率-波数谱，从而产生 该地震数据的轨迹内插。

2.如权利要求1所述的方法，其中该第一防泄漏傅立叶变换方 法为标准防泄漏傅立叶变换方法。

3.如权利要求1所述的方法，其中非混淆频率基本上为地震数 据中较低的频率而混淆频率基本上为地震数据中较高的频率。

4.如权利要求1所述的方法，其中时间傅立叶变换为快速傅立 叶变换。

5.如权利要求1所述的方法，进一步包括：

为第一防泄漏傅立叶变换方法确定迭代数目N₁。

6.如权利要求5所述的方法，其中确定迭代数目N₁包括：

为第一防泄漏傅立叶变换方法确定阈值ε₁。

7.如权利要求6所述的方法，进一步包括：

从所述非混淆频率选择非混淆频率分量；

针对非混淆频率将所述估计频率-波数谱初始化为零；

将空间傅立叶变换应用于所述非混淆频率分量以获得变换后的 非混淆频率分量；以及

针对N₁次迭代执行下列步骤：

在所述变换后的非混淆频率分量中选择最大波数分量；

将最大波数分量增加到估计频率-波数谱中；

将逆空间傅立叶变换应用于所述最大波数分量以获得逆变 换后的最大波数分量；以及

从所述非混淆频率分量减去所述逆变换后的最大波数分量， 从而产生纠正的非混淆频率分量。

8.如权利要求7所述的方法，其中空间傅立叶变换为离散傅立 叶变换，以及逆空间傅立叶变换为逆离散傅立叶变换。

9.如权利要求7所述的方法，其中空间傅立叶变换为非均匀快 速傅立叶变换，以及逆空间傅立叶变换为逆非均匀快速傅立叶变换。

10.如权利要求7所述的方法，进一步包括：

组合所述纠正的非混淆频率分量的估计频率-波数谱，从而产生 非混淆估计频率-波数谱；

取所述纠正的非混淆频率分量的绝对值，从而产生绝对频率-波 数谱；

将绝对频率-波数谱外推到混淆频率，从而产生外推绝对频率-波 数谱；以及

将第二防泄漏傅立叶变换方法应用于来自变换后的地震数据的 外推绝对频率-波数谱的每一个混淆频率分量。

11.如权利要求10所述的方法，进一步包括：

在将绝对频率-波数谱外推到混淆频率之前将绝对频率-波数谱 平滑。

12.如权利要求10所述的方法，进一步包括：

为第二防泄漏傅立叶变换方法确定迭代数目N₂。

13.如权利要求12所述的方法，其中确定迭代数目N₂包括：

为该第二防泄漏傅立叶变换方法确定阈值ε₂。