[发明专利]一种高维地震数据规则化方法

说明书

技术领域

本发明属于油气勘探开发中地震数据预处理领域，具体涉及一种高维地震数据规则化方法。

背景技术

在地震数据采集中，空间采样的不规则以及空间采样间隔过大(采集成本限制)导致地震数据的空间假频经常发生。由于偏移算子不能消除由空间采样引起的假频，以致这些假频能量以假象的形式出现在剖面中，这些假象会误导地震剖面的解释。因此为了实现真振幅成像，必须在偏移前通过地震数据插值等方法除去这些假频能量。在过去的十多年里，地球物理学家发展了很多解决地震数据非规则性和假频问题的方法，大体上可分为四类。描述如下：

第一类，基于褶积理论。这种方法的主要思想是，用地震数据的低频分量构造一个滤波器，并利用这个滤波器实现对地震数据的插值。属于这种类型的方法包括FX域预测滤波插值(Spitz，1991)，TX域预测误差滤波插值(Claerbout，1992)，FK域预测滤波插值(Gulunay，2003)。它们通常要求数据域的同相轴是线性的，只有满足线性假设条件，所构造的滤波器才有意义。这些方法的另一个要求是地震数据的空间采样方式必须是规则的，但是在数据采集过程中这个要求是很难实现的。

第二类，基于偏移算子和反偏移算子。这种方法先应用偏移算子将地震数据投影到成像域中，通过在成像域里完成地震数据的插值工作，然后将成像域的数据反偏移回数据域中。由Stolt(2002)和Trad(2003)介绍的方法都属于这种类型。这种方法的偏移算子种类繁多，包括正常时差校正(NMO)，倾角时差校正(DMO)，零偏移距偏移(MZO)，方位角时差校正(AMO)，叠前时间偏移(PSTM)，叠前深度偏移(PSDM)等。这类方法虽然可以很好的保持地震数据中同相轴的相位信息，但是振幅信息往往是失真的。

第三类，基于FK频谱估计。这种方法可以另外分为两小类：一种基于反演理论，另一种基于相关理论。基于反演理论的主流方法可参考A.J.W.Duijndam(1999)和Sacchi(2004，2007)的文章。这类方法的主要思想是，在反演理论的框架下，应用合适的正演算子和约束条件，估计地震数据的FK频谱。这类反演问题的模型空间和数据空间通常是FX域和FK域的地震数据，而正算子是傅立叶变换。这些方法的约束条件差异很大，不同的约束条件会产生不同的结果，因此约束条件是该方法的关键因素。总所周知，一个反问题经常会对应着多个解，不同的限制条件意味着在不同的解空间寻求结果，因此对于不同的约束条件而言解是完全不同的。另一类基于相关理论的方法包括抗泄露傅立叶变换(ALFT)(Sheng Xu，2005)，凸集投影(POCS)(Ray Abma，2006)，匹配追踪(MP)(Al i Ozbek，)以及其他一些与基于地震数据相关的方法。在相关过程中，这些方法将相关能量最大的基作为地震数据中能量最大的分量，然后选择这个最大能量作为FK谱估计的一个分量。实际上如果将相关理论看作解反演问题的一种方法，那么每次选择的基就是梯度，而相关能量就是这个反演问题的步长。如果数据中的假频不严重，那么这类方法能得到很好的结果。但是如果假频严重，最大相关能量的基可能就不再是地震数据有效信号的分量而是假频，则此方法得到的插值结果就很差。

最后一种，基于矩阵分解。如果将地震数据看作一个高维矩阵，那么在奇异值分解(SVD)中，数据插值可以通过插值奇异向量得到。由Sacchi(2011)提出的方法就属于这一类。如果数据域中的同相轴是线性的，那这种方法将得到更好的结果，因为如果同相轴是线性的，那么大奇异值的所对应的奇异向量的数量将大大减少。但是其缺陷也是显而易见的，因为地震数据的随机性很强，因此在奇异值分解时就会产生奇异值从而导致矩阵误解。

上述四类方法中，最后两类是当今地震数据插值中最为流行的，因为它们都具备高维地震数据插值的能力，并且插值效果也比较理想。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种高维地震数据规则化方法，解决地震数据采集中，空间采样的不规则以及空间采样间隔过大(采集成本限制)导致地震数据的空间假频的问题，这些假频能量以假象的形式出现在剖面中，这些假象会误导地震剖面的解释。本方法是在反演过程中以高维傅里叶变换为正演算子，数据线性拉东变换谱为约束的高维数据规则化方法。

本发明是通过以下技术方案实现的：

一种高维地震数据规则化方法，所述方法包括：

(1)输入原始三维炮道集；

(2)利用高维傅里叶变换将时间空间域的所述原始三维炮道集变换到频率空间域中，然后进行分选得到五维炮道集；