一种压制虚反射信号的方法及装置的制造方法
【技术领域】
[0001] 本发明涉及海上地震资料处理技术领域，尤其涉及一种压制虚反射信号的方法及 装置。
【背景技术】
[0002] 地球物理勘探技术领域中，虚反射往往是影响海上地震资料品质的一个重要因 素。在海上进行地震勘探时，一般将震源和检波器放置于海面以下一定深度处，由于海水与 空气的波阻抗差异较大，往往在海水与空气之间形成较强的反射面，所述反射面构成自由 界面。震源激发的地震波经过所述自由界面反射后可以形成虚反射。所述虚反射的波形、 频率、视速度等参数一般和一次反射波相似，因此，所述虚反射可能会干扰一次反射波，造 成地震记录的低频响应，降低地震剖面的辨率。所述虚反射有时甚至会产生虚假同相轴，给 地震资料反演与解释造成较大的困难。因此，压制虚反射信号是海上地震资料处理最重要 的步骤之一。
[0003] 现有技术中常用反褶积方法压制虚反射信号，具体可以设计虚反射滤波器，在时 间域利用预测反褶积来压制虚反射信号。所述反褶积方法假设地震记录是反射波与虚反射 滤波器褶积以及所述虚反射滤波器为最小相位，实际上所述虚反射滤波器的振幅谱存在一 系列的"零值点"，所以所述虚反射滤波器不是最小相位的，不满足反褶积的条件，因此时间 域反褶积方法不能完全消除虚反射。
【发明内容】
[0004] 本申请的目的在于提供一种压制虚反射信号的方法及装置，可以压制虚反射信 号，消除虚反射信号的干扰。
[0005] 本申请提供的一种压制虚反射信号的方法及装置是这样实现的：
[0006] -种压制虚反射信号的方法，所述方法包括：
[0007] 获取频率-波数域的海上地震数据；
[0008] 基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域 的虚反射压制因子；
[0009] 根据所述虚反射压制因子计算虚反射压制算子；
[0010] 根据所述虚反射压制算子对所述频率-波数域的海上地震数据进行虚反射信号 压制处理，得到虚反射信号压制后的海上地震数据。
[0011] 可选的，在本发明的一个实施例中，所述基于检波器在海水中的深度、所述海上地 震数据的频率、波数，计算频率-波数域的虚反射压制因子包括：
[0012] 采用下式计算所述虚反射压制因子：
[0014] 式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地 震数据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在 海水中的深度，C为海水速度。
[0015] 可选的，在本发明的一个实施例中，所述根据所述虚反射压制因子对所述频 率-波数域的海上地震数据进行虚反射信号压制处理包括：
[0016] 采用下式计算虚反射压制后的海上地震数据：
[0018] 式中，A(F，K)为虚反射压制算子在频率-波数域的值，G(F，K)为虚反射压制因子 在频率-波数域的值，G丨厂，幻为 G(F，κ)的共辄。
[0019] 可选的，在本发明的一个实施例中，所述根据所述虚反射压制算子对所述频 率-波数域的海上地震数据进行虚反射信号压制处理，包括：
[0020] 将所述频率-波数域的海上地震数据和所述虚反射压制算子做乘积运算。
[0021] 可选的，在本发明的一个实施例中，所述根据所述虚反射压制因子对所述频 率-波数域的海上地震数据进行虚反射信号压制处理包括：
[0022] 采用下式计算虚反射压制后的海上地震数据：
[0023] P (F, Κ) = D (F, k) [A (F, K) +white]
[0024] 式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压 制虚反射信号前的频率-波数域的海上地震数据，A (F，K)为虚反射压制算子在频率-波数 域的值，white为白噪声系数。
[0025] 可选的，在本发明的一个实施例中，所述白噪声系数的计算公式如下式所示：
[0027] 式中，white为白噪声系数，F为所述海上地震数据的频率。
[0028] 可选的，在本发明的一个实施例中，所述获取频率-波数域的海上地震数据，包 括：
[0029] 获取时间-空间域的海上地震数据；
[0030] 对所述时间-空间域的海上地震数据进行二维傅里叶变换。
[0031] -种压制虚反射信号的装置，所述装置包括：
[0032] 地震数据获取单元，用于获取频率-波数域的海上地震数据；
[0033] 压制因子计算单元，用于根据所述海上地震数据的分析处理结果，计算虚反射压 制因子；
[0034] 虚反射压制算子计算单元，用于根据所述虚反射压制因子计算虚反射压制算子；
[0035] 压制后地震数据计算单元，用于根据所述海上地震数据以及所述虚反射压制算子 计算压制虚反射信号后的海上地震数据。
[0036] 可选的，在本发明的一个实施例中，所述压制因子计算单元包括：
[0037] 所述虚反射压制因子的计算公式包括：
[0039] 式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为海上地震数 据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，ζ为检波器在海水 中的深度，c为海水速度。
[0040] 可选的，在本发明的一个实施例中，所述虚反射压制算子计算单元包括：
[0041] 所述虚反射压制算子的计算公式包括：
[0043] 式中，A(F，K)为虚反射压制算子在频率-波数域的值，G(F，K)为虚反射压制因子 在频率-波数域的值，Z 7. /〇为G (F，K)的共辄。
[0044] 可选的，在本发明的一个实施例中，所述压制后地震数据计算单元包括：
[0045] 计算子单元，用于将所述频率-波数域的海上地震数据和所述虚反射压制算子做 乘积运算。
[0046] 可选的，在本发明的一个实施例中，所述压制后地震数据计算单元包括：
[0047] 所述虚反射压制因子的计算公式包括：
[0048] P (F, K) = D (F, k) [A (F, K) +white]
[0049] 式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压 制虚反射信号前的频率-波数域的海上地震数据，A (F，K)为虚反射压制算子在频率-波数 域的值，white为白噪声系数。
[0050] 可选的，在本发明的一个实施例中，所述白噪声系数的计算公式如下式所示：
[0052] 式中，white为白噪声系数，F为所述海上地震数据的频率。
[0053] 可选的，在本发明的一个实施例中，所述地震数据获取单元，包括：
[0054] 时深域地震数据获取单元，用于获取时间-空间域的海上地震数据；
[0055] 傅里叶变换单元，用于对所述时间-空间域的海上地震数据进行二维傅里叶变 换。
[0056] 由此可见，本发明的一种压制虚反射信号的方法和装置的技术方案通过对频 率-波数域内的海