相同相似的部分互相参见即可。尤其，对于系统实 施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例 的部分说明即可。
[0094] 本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块 (illustrativelogicalblock)，单元，和步骤可以通过电子硬件、电脑软件，或两者的结 合进行实现。为清楚展示硬件和软件的可替换性（interchangeability)，上述的各种说明 性部件（illustrativecomponents)，单元和步骤已经通用地描述了它们的功能。这样的功 能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员 可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为 超出本发明实施例保护的范围。
[0095] 本发明实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器， 数字信号处理器，专用集成电路（ASIC)，现场可编程门阵列（FPGA)或其它可编程逻辑装 置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功 能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制 器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和 微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类 似的配置来实现。
[0096] 本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软 件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存 储器、EEPR0M存储器、寄存器、硬盘、可移动磁盘、⑶-ROM或本领域中其它任意形式的存储 媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信 息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储 媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以 设置于用户终端中的不同的部件中。
[0097] 在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软 件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的 媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以 是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但 不限于RAM、ROM、EEPR0M、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它 任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处 理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例 如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双 绞线、数字用户线（DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义 的电脑可读媒介中。所述的碟片（disk)和磁盘（disc)包括压缩磁盘、镭射盘、光盘、DVD、 软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的 组合也可以包含在电脑可读媒介中。
[0098] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种压制虚反射信号的方法，其特征在于，所述方法包括： 获取频率-波数域的海上地震数据； 基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚 反射压制因子； 根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制 处理，得到虚反射信号压制后的海上地震数据。2. 根据权利要求1所述的一种压制虚反射信号的方法，其特征在于，所述基于检波器 在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子包 括： 采用下式计算所述虚反射压制因子：式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地震数 据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在海水 中的深度，c为海水速度。3. 根据权利要求1所述的一种压制虚反射信号的方法，其特征在于，所述根据所述虚 反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理包括： 采用下式计算虚反射压制后的海上地震数据：式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压制虚 反射信号前的频率-波数域的海上地震数据，G (F，K)为虚反射压制因子在频率-波数域的 值，white为白噪声系数。4. 根据权利要求1至3中任意一项所述的一种压制虚反射信号的方法，其特征在于，所 述方法还包括： 对所述压制虚反射信号后的频率-波数域的海上地震数据进行二维傅里叶反变换，生 成压制虚反射信号后的时间-空间域的海上地震数据。5. 根据权利要求1至3中任意一项所述的一种压制虚反射信号的方法，其特征在于，所 述获取频率-波数域的海上地震数据，包括： 获取时间-空间域的海上地震数据； 对所述时间-空间域的海上地震数据进行二维傅里叶变换。6. 根据权利要求5所述的一种压制虚反射信号的方法，其特征在于，所述对所述时 间-空间域的海上地震数据进行二维傅里叶变换包括： 采用下式对所述时间-空间域的海上地震数据进行二维傅里叶变换：式中，D (F，K)为压制虚反射信号前的频率-波数域的海上地震数据，h(x，t)为压制虚 反射信号前的时间-空间域的海上地震数据，F和K分别为海上地震数据的频率和波数，t 和X分别为海上地震数据的纵向采样时间和采样点深度。7. -种压制虚反射信号的装置，其特征在于，所述装置包括： 地震数据获取单元，用于获取频率-波数域的海上地震数据； 压制因子计算单元，用于基于检波器在海水中的深度、所述海上地震数据的频率、波 数，计算频率-波数域的虚反射压制因子； 压制后地震数据计算单元，用于根据所述虚反射压制因子对所述频率-波数域的海上 地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。8. 根据权利要求7所述的一种压制虚反射信号的装置，其特征在于，所述压制因子计 算单元包括： 所述虚反射压制因子的计算公式包括：式中，G(F，K)为虚反射压制因子在频率-波数域的值，F和K分别为所述海上地震数 据的频率和波数，R为空气和海面间自由界面的反射系数，i为虚数单位，z为检波器在海水 中的深度，c为海水速度。9. 根据权利要求7所述的一种压制虚反射信号的装置，其特征在于，所述压制后地震 数据计算单元包括： 所述虚反射压制因子的计算公式包括：式中，P(F，K)为压制虚反射信号后的频率-波数域的海上地震数据，D(F，K)为压制虚 反射信号前的频率-波数域的海上地震数据，G (F，K)为虚反射压制因子在频率-波数域的 值，white为白噪声系数。10. 根据权利要求7至9中任意一项所述的一种压制虚反射信号的装置，其特征在于， 所述装置还包括： 傅里叶反变换单元，用于对所述压制虚反射信号后的频率-波数域的海上地震数据进 行二维傅里叶反变换，生成压制虚反射信号后的时间-空间域的海上地震数据。11. 根据权利要求7至9中任意一项所述的一种压制虚反射信号的装置，其特征在于， 所述地震数据获取单元，包括： 时深域地震数据获取单元，用于获取时间-空间域的海上地震数据； 傅里叶变换单元，用于对所述时间-空间域的海上地震数据进行二维傅里叶变换。12. 根据权利要求11所述的一种压制虚反射信号的装置，其特征在于，所述傅里叶变 换单元包括： 所述对所述时间-空间域的海上地震数据进行二维傅里叶变换的计算公式包括： D(F,K) =I f h(x, t)e j2ll(Ft+Kx)dxdt 式中，D (F，K)为压制虚反射信号前的频率-波数域的海上地震数据，h(x，t)为压制虚 反射信号前的时间-空间域的海上地震数据，F和K分别为海上地震数据的频率和波数，t 和X分别为海上地震数据的纵向采样时间和采样点深度。
【专利摘要】本发明提供一种压制虚反射信号的方法及装置。所述方法包括：获取频率-波数域的海上地震数据；基于检波器在海水中的深度、所述海上地震数据的频率、波数，计算频率-波数域的虚反射压制因子；根据所述虚反射压制因子对所述频率-波数域的海上地震数据进行虚反射信号压制处理，得到虚反射信号压制后的海上地震数据。利用所述本发明实施例提供的技术方案，计算得到比较准确的虚反射压制因子和稳定性较强的压制虚反射信号后的海上地震数据，可以压制虚反射信号，消除虚反射信号的干扰。
【IPC分类】G01V1/36
【公开号】CN105425297
【申请号】CN201510716957
【发明人】马光凯, 周铮铮, 赵波, 钱忠平, 任晓乔, 李合群
【申请人】中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司
【公开日】2016年3月23日
【申请日】2015年10月29日