CN105425297A - 一种压制虚反射信号的方法及装置 - Google Patents

Abstract

本发明提供一种压制虚反射信号的方法及装置。所述方法包括：获取频率-波数域的海上地震数据；基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子；根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。利用所述本发明实施例提供的技术方案，计算得到比较准确的虚反射压制因子和稳定性较强的压制虚反射信号后的海上地震数据，可以压制虚反射信号，消除虚反射信号的干扰。

Description

一种压制虚反射信号的方法及装置

技术领域

本发明涉及海上地震资料处理技术领域，尤其涉及一种压制虚反射信号的方法及装置。

背景技术

地球物理勘探技术领域中，虚反射往往是影响海上地震资料品质的一个重要因素。在海上进行地震勘探时，一般将震源和检波器放置于海面以下一定深度处，由于海水与空气的波阻抗差异较大，往往在海水与空气之间形成较强的反射面，所述反射面构成自由界面。震源激发的地震波经过所述自由界面反射后可以形成虚反射。所述虚反射的波形、频率、视速度等参数一般和一次反射波相似，因此，所述虚反射可能会干扰一次反射波，造成地震记录的低频响应，降低地震剖面的辨率。所述虚反射有时甚至会产生虚假同相轴，给地震资料反演与解释造成较大的困难。因此，压制虚反射信号是海上地震资料处理最重要的步骤之一。

现有技术中常用反褶积方法压制虚反射信号，具体可以设计虚反射滤波器，在时间域利用预测反褶积来压制虚反射信号。所述反褶积方法假设地震记录是反射波与虚反射滤波器褶积以及所述虚反射滤波器为最小相位，实际上所述虚反射滤波器的振幅谱存在一系列的“零值点”，所以所述虚反射滤波器不是最小相位的，不满足反褶积的条件，因此时间域反褶积方法不能完全消除虚反射。

发明内容

本发明的目的在于提供一种压制虚反射信号的方法及装置，可以压制虚反射信号，消除虚反射信号的干扰。

本发明提供的一种压制虚反射信号的方法及装置具体是这样实现的：

一种压制虚反射信号的方法，所述方法包括：

获取频率-波数域的海上地震数据；

基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子；

根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。

可选的，在本发明的一个实施例中，所述基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子包括：

采用下式计算所述虚反射压制因子：

$G(F,K)=1+{\mathrm{Re}}^{-\frac{2iz}{c}\sqrt{F^2-c^2{K}^2}}$

式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地震数据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在海水中的深度，c为海水速度。

可选的，在本发明的一个实施例中，所述根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理包括：

采用下式计算虚反射压制后的海上地震数据：

$P(F,K)=\frac{D(F,K)}{G(G,K)+white}$

式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，G(F，K)为虚反射压制因子在频率-波数域的值，white为白噪声系数。

可选的，在本发明的一个实施例中，此后还包括：

对所述压制虚反射信号后的频率-波数域的海上地震数据进行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。

可选的，在本发明的一个实施例中，所述获取频率-波数域的海上地震数据，包括：

获取时间-空间域的海上地震数据；

对所述时间-空间域的海上地震数据进行二维傅里叶变换。

可选的，在本发明的一个实施例中，所述对所述时间-空间域的海上地震数据进行二维傅里叶变换包括：

采用下式对所述时间-空间域的海上地震数据进行二维傅里叶变换：

D(F,K)＝∫∫h(x,t)e^-j2π(Ft+Kx)dxdt

式中，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，h(x,t)为压制虚反射信号前的时间-空间域的海上地震数据，F和K分别为海上地震数据的频率和波数，t和x分别为海上地震数据的纵向采样时间和采样点深度。

一种压制虚反射信号的装置，所述装置包括：

地震数据获取单元，用于获取频率-波数域的海上地震数据；

压制因子计算单元，用于基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子；

压制后地震数据计算单元，用于根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。

可选的，在本发明的一个实施例中，所述压制因子计算单元包括：

所述虚反射压制因子的计算公式包括：

$G(F,K)=1+{\mathrm{Re}}^{-\frac{2iz}{c}\sqrt{F^2-c^2{K}^2}}$

式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地震数据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在海水中的深度，c为海水速度。

可选的，在本发明的一个实施例中，所述压制后地震数据计算单元包括：

所述虚反射压制因子的计算公式包括：

$P(F,K)=\frac{D(F,K)}{G(G,K)+white}$

式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，G(F，K)为虚反射压制因子在频率-波数域的值，white为白噪声系数。

可选的，在本发明的一个实施例中，所述装置还包括：

傅里叶反变换单元，用于对所述压制虚反射信号后的频率-波数域的海上地震数据进行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。

可选的，在本发明的一个实施例中，所述地震数据获取单元，包括：

时深域地震数据获取单元，用于获取时间-空间域的海上地震数据；

傅里叶变换单元，用于对所述时间-空间域的海上地震数据进行二维傅里叶变换。

可选的，在本发明的一个实施例中，所述傅里叶变换单元包括：

所述对所述时间-空间域的海上地震数据进行二维傅里叶变换的计算公式包括：

D(F,K)＝∫∫h(x,t)e^-j2π(Ft+Kx)dxdt

式中，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，h(x,t)为压制虚反射信号前的时间-空间域的海上地震数据，F和K分别为海上地震数据的频率和波数，t和x分别为海上地震数据的纵向采样时间和采样点深度。

由此可见，本发明的一种压制虚反射信号的方法和装置的技术方案通过对频率-波数域内的海上地震数据进行分析处理，考虑到反射系数、检波器深度、海水速度等各种因素的影响，计算得到比较准确的虚反射压制因子，再根据所述虚反射压制因子计算得到稳定性较强的压制虚反射信号后的海上地震数据，可以消除虚反射信号的干扰。消除虚反射信号干扰后，不仅可以消除虚反射的同相轴，使得有效信号的同相轴变得更加清晰、连续，还可以消除由于虚反射信号影响得到的陷波点，补偿陷波点处的能量，拓宽海上地震数据的有效带宽，提升海上地震数据的低频信息。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本发明提供的压制虚反射信号方法的一种实施例的方法流程图；

图2是本发明提供的采用理论模型合成的海上地震数据的同相轴曲线图；

图3是本发明提供的采用理论模型合成的海上地震数据的频谱图；

图4是本发明对海上地震数据采用本发明实施例方法处理后的同相轴曲线图；

图5是本发明对海上地震数据采用本发明实施例方法处理后的频谱图；

图6是本发明实施例提供的压制虚反射信号装置的一种实施例的模块结构示意图；

图7是本发明提供的海上地震数据获取单元的一种实施例的模块结构示意图。

具体实施方式

为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

在时间-空间域，有用信号和干扰信号一般不易被分离，同理，在时间-空间域进行虚反射压制时，不能轻易分离反射波和虚反射。如果将在时间-空间域中不易分离的反射波和虚反射转换至频率-波数域，将有利于反射波信号和虚反射信号的处理。本发明提供了一种压制虚反射信号的方法及装置。图1是本发明提供的压制虚反射信号方法的一种实施例的方法流程图。如图所示，所述方法包括：

步骤S1：获取频率-波数域的海上地震数据。

海上地震数据的采集和陆上采集相似，在震源处激发地震波，在接收器接收反射波，不同的地方在于，海上地震数据采集时所用的震源和接收器都是位于海面以下，并且海上地震数据中一般存在虚反射信号的干扰。在一次海上地震数据采集过程中，一般是多炮激发，每炮激发可以产生多个地震道，一个地震道中可以包括多个采样点。

在本发明的一个实施例中，所述获取频率-波数域的海上地震数据可以包括以下步骤：

步骤SS1：获取时间-空间域的海上地震数据；

海上地震数据的原始数据往往是时间-空间域的采集信号，可以用h(x,t)表示，其中x可以表示检测器到震源点的距离，t可以表示纵向采集的时间。在时间-空间域内，海上地震数据h(x,t)中的反射波信号和虚反射信号往往是相互穿插叠合的，不容易分离。但是，海上地震数据h(x,t)在时间和深度上都可以分解为正弦和余弦形式的波动组分，转换成二维的频谱。单个波动组分在时间上的额度，以每秒中的波动次数来计量，可以生成信号的频率，在深度上的额度，以每米中的波动次数来计量，可以生成信号的波数。

步骤SS2：对所述时间-空间域的海上地震数据进行二维傅里叶变换。

本实施例中可以将海上地震数据h(x,t)从时间-空间域转换至频率-波数域，具体可以对海上地震数据h(x,t)进行傅里叶变换，计算公式如下式(1)所示：

D(F,K)＝∫∫h(x,t)e^-j2π(Ft+Kx)dxdt(1)

式中，D(F,K)可以表示频率-波数域的海上地震数据。

运用二维傅里叶变换将海上地震数据从时间-空间域转换至频率-波数域，可以将所述海上地震数据的能量集中于频率-波数域中的一个矩形范围内。所述海上地震数据中视速度不同的信号，能量集中方向也不相同，同时，在所述地震信号中，反射波信号和虚反射信号的视速度一般是有较大差异的，因此，对海上地震数据从时间-空间域转换至频率-波数域可以为压制虚反射信号提供了可能。

步骤S2：基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子。

在本发明的一个实施例中，所述虚反射压制因子的计算公式如下式(2)所示：

$G(F,K)=1+{\mathrm{Re}}^{-\frac{2iz}{c}\sqrt{F^2-c^2{K}^2}}---\left(2\right)$

式中，G(F，K)可以为虚反射压制因子在频率-波数域的值，F和K可以分别为海上地震数据的频率和波数，R可以为空气和海面间自由界面的反射系数，i可以为虚数单位，z可以为检波器在海水中的深度，c可以为海水速度，海水速度一般取值为1500m/s。

在压制虚反射信号过程中，虚反射压制因子的计算精度可以影响虚反射的压制能力。在海上地震激发过程中，空气和海面间自由界面的反射系数和震源及检波器的放置深度、炮检距、海水面状况等因素有关，观察式(2)，可以发现，R值越大，虚反射压制因子值也越大。式(2)中，检波器在海水中的深度可以从所述海上地震数据中直接获取，检波器在海水中的深度一般可以记录在单炮中每一记录道记录数据的开始处。在本发明的一个实施例中，反射系数R的取值范围可以包括0.001到1.0，优选的，反射系数R的值可以包括0.97。

所述虚反射压制因子的计算考虑了反射系数、检波器深度、海水速度等因素的影响，可以提高虚反射压制因子的准确度，提高虚反射的压制能力。

步骤S3：根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。

在本发明的一个实施例中，所述压制虚反射信号后的频率-波数域的海上地震数据的计算公式如下式(3)所示：

$P(F,K)=\frac{D(F,K)}{G(G,K)+white}---\left(3\right)$

式中，P(F，K)可以为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)可以为压制虚反射信号前的频率-波数域的海上地震数据，G(F，K)可以为虚反射压制因子在频率-波数域的值，white可以为白噪声系数。所述白噪声系数可以增加计算的稳定性，所述白噪声系数的取值范围可以包括0.0001到0.1。

后续地，所述方法还包括：对所述压制虚反射信号后的海上地震数据进行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。所述傅里叶反变换的公式可以如下式(4)所示：

h'(x,t)＝∫∫P(F,K)e^j2π(Ft+Kx)dfdk(4)

式中，h’(x,t)可以为压制虚反射信号后的时间-空间域的海上地震数据。

下面通过结合具体的场景应用上述方法步骤，图2是本发明提供的采用理论模型合成的海上地震数据的同相轴曲线图，如图2所示，图中的同相轴曲线可以是所述海上地震数据中各记录道振动相位相同的极值的连线，由于虚反射的影响，从所述同相轴曲线来看，所述虚反射信号的同相轴和一次波(属于有用反射波)同相轴被混在一起，图2中深色曲线可以表示一次波同相轴曲线，图2中浅色曲线可以表示虚反射的同相轴。同时，图3是本发明提供的采用理论模型合成的海上地震数据的频谱图，如图3所示，由于虚反射的影响，在所述频谱图中的65Hz左右的位置出现明显的陷波点，所述陷波点可以限制有效频带的宽带，损失海上地震数据中一部分有意义的低频信息和高频信息，所述低频信息中可以包括一部分对于较深地质体成像有较大影响的低频信息。图4是本发明对海上地震数据采用本发明实施例方法处理后的同相轴曲线图。从图4可以发现，通过采用本发明实施例提供的将海上地震数据在频率-波数域中进行虚反射压制处理后，虚反射的同相轴已经被消除，所述一次波的同相轴曲线变得更加清晰、更加连续。可以推断已经消除虚反射对有用反射波的干扰。图5是本发明对海上地震数据采用本发明实施例方法处理后的频谱图。从图5可以发现，通过采用本发明实施例提供的将海上地震数据在频率-波数域中进行虚反射压制处理后，所述陷波点被消除，所述陷波点处的能量得到补偿，拓宽海上地震数据的有效带宽，提升海上地震数据的低频信息。

本发明另一方面还提供一种压制虚反射信号的装置，图6是本发明提供的压制虚反射信号装置的一种实施例的模块结构示意图，结合附图6，装置60可以包括：地震数据获取单元61、压制因子计算单元62、压制后地震数据计算单元63，其中，

地震数据获取单元61，用于获取频率-波数域的海上地震数据。

压制因子计算单元62，用于根据所述海上地震数据的分析处理结果，计算虚反射压制因子。

压制后地震数据计算单元63，用于根据所述海上地震数据以及所述虚反射压制因子计算压制虚反射信号后的海上地震数据。

图7是本发明提供的地震数据获取单元的一种实施例的模块结构示意图，如图7所示，所述地震数据获取单元61可以包括：时深域地震数据获取单元71、傅里叶变换单元72，其中，

时深域地震数据获取单元71，用于获取时间-空间域的海上地震数据。

傅里叶变换单元72，用于对所述时间-空间域的海上地震数据进行二维傅里叶变换。

所述时深域地震数据获取单元71和所述傅里叶变换单元72运用二维傅里叶变换将海上地震数据从时间-空间域转换至频率-波数域，可以将所述海上地震数据的能量集中于频率-波数域中的一个矩形范围内，对海上地震数据从时间-空间域转换至频率-波数域可以为压制虚反射信号提供了可能。

如图6所示，所述装置还可以包括傅里叶反变换单元64，所述傅里叶反变换单元64用于对所述压制虚反射信号后的频率-波数域的海上地震数据进行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。

由此可见，本发明的一种压制虚反射信号的方法和装置的技术方案通过对频率-波数域内的海上地震数据进行分析处理，考虑到反射系数、检波器深度、海水速度等各种因素的影响，计算得到比较准确的虚反射压制因子，再根据所述虚反射压制因子计算得到稳定性较强的压制虚反射信号后的海上地震数据，可以消除虚反射信号的干扰。消除虚反射信号干扰后，不仅可以消除虚反射的同相轴，使得有效信号的同相轴变得更加清晰、连续，还可以消除由于虚反射信号影响得到的陷波点，补偿陷波点处的能量，拓宽海上地震数据的有效带宽，提升海上地震数据的低频信息。

本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。尤其，对于***实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrativelogicalblock)，单元，和步骤可以通过电子硬件、电脑软件，或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability)，上述的各种说明性部件(illustrativecomponents)，单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个***的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。

本发明实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器，数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列(FPGA)或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。

本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。

在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (12)

1.一种压制虚反射信号的方法，其特征在于，所述方法包括：

获取频率-波数域的海上地震数据；

基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子；

根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。

2.根据权利要求1所述的一种压制虚反射信号的方法，其特征在于，所述基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子包括：

采用下式计算所述虚反射压制因子：

$G(F,K)=1+{\mathrm{Re}}^{-\frac{2iz}{c}\sqrt{F^2-c^2{K}^2}}$

式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地震数据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在海水中的深度，c为海水速度。

3.根据权利要求1所述的一种压制虚反射信号的方法，其特征在于，所述根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理包括：

采用下式计算虚反射压制后的海上地震数据：

$P(F,K)=\frac{D(F,K)}{G(G,K)+white}$

式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，G(F，K)为虚反射压制因子在频率-波数域的值，white为白噪声系数。

4.根据权利要求1至3中任意一项所述的一种压制虚反射信号的方法，其特征在于，所述方法还包括：

对所述压制虚反射信号后的频率-波数域的海上地震数据进行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。

5.根据权利要求1至3中任意一项所述的一种压制虚反射信号的方法，其特征在于，所述获取频率-波数域的海上地震数据，包括：

获取时间-空间域的海上地震数据；

对所述时间-空间域的海上地震数据进行二维傅里叶变换。

6.根据权利要求5所述的一种压制虚反射信号的方法，其特征在于，所述对所述时间-空间域的海上地震数据进行二维傅里叶变换包括：

采用下式对所述时间-空间域的海上地震数据进行二维傅里叶变换：

D(F,K)＝∫∫h(x,t)e^-j2π(Ft+Kx)dxdt

式中，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，h(x,t)为压制虚反射信号前的时间-空间域的海上地震数据，F和K分别为海上地震数据的频率和波数，t和x分别为海上地震数据的纵向采样时间和采样点深度。

7.一种压制虚反射信号的装置，其特征在于，所述装置包括：

地震数据获取单元，用于获取频率-波数域的海上地震数据；

压制因子计算单元，用于基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子；

压制后地震数据计算单元，用于根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。

8.根据权利要求7所述的一种压制虚反射信号的装置，其特征在于，所述压制因子计算单元包括：

所述虚反射压制因子的计算公式包括：

$G(F,K)=1+{\mathrm{Re}}^{-\frac{2iz}{c}\sqrt{F^2-c^2{K}^2}}$

式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地震数据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在海水中的深度，c为海水速度。

9.根据权利要求7所述的一种压制虚反射信号的装置，其特征在于，所述压制后地震数据计算单元包括：

所述虚反射压制因子的计算公式包括：

$P(F,K)=\frac{D(F,K)}{G(G,K)+white}$

式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，G(F，K)为虚反射压制因子在频率-波数域的值，white为白噪声系数。

10.根据权利要求7至9中任意一项所述的一种压制虚反射信号的装置，其特征在于，所述装置还包括：

傅里叶反变换单元，用于对所述压制虚反射信号后的频率-波数域的海上地震数据进行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。

11.根据权利要求7至9中任意一项所述的一种压制虚反射信号的装置，其特征在于，所述地震数据获取单元，包括：

时深域地震数据获取单元，用于获取时间-空间域的海上地震数据；

傅里叶变换单元，用于对所述时间-空间域的海上地震数据进行二维傅里叶变换。

12.根据权利要求11所述的一种压制虚反射信号的装置，其特征在于，所述傅里叶变换单元包括：

所述对所述时间-空间域的海上地震数据进行二维傅里叶变换的计算公式包括：

D(F,K)＝∫∫h(x,t)e^-j2π(Ft+Kx)dxdt

式中，D(F，K)为压制虚反射信号前的频率-波数域的海上地震数据，h(x,t)为压制虚反射信号前的时间-空间域的海上地震数据，F和K分别为海上地震数据的频率和波数，t和x分别为海上地震数据的纵向采样时间和采样点深度。