虚反射的压制方法和装置制造方法及图纸

本申请实施方式提供了一种虚反射的压制方法和装置，其中，该方法包括：获取目标区域的海洋地震数据，海洋地震数据包括多个单炮地震数据；根据单炮地震数据，通过波场反向延拓，确定反向延拓数据；根据反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱；对虚反射压制算子的相位谱进行希尔伯特变换，得到虚反射压制算子的振幅谱；确定虚反射压制算子；利用虚反射压制算子，对海洋地震数据进行虚反射压制处理。由于该方案通过相位驱动确定虚反射压制算子，再利用虚反射压制算子对海洋地震数据进行去虚反射处理，因而解决了现有的虚反射的压制方法中存在的不能完全消除虚反射、压制效果差、应用范围有限的技术问题。

【技术实现步骤摘要】
虚反射的压制方法和装置
本申请涉及地球物理勘探
，特别涉及一种虚反射的压制方法和装置。
技术介绍
在进行海底石油地震勘探时，常常需要获取海洋地震数据，根据海洋地震数据进行叠加成像、地震解释等，以进行具体的海底石油地震勘探。具体获取海底地震数据时，往往是在海上通过拖缆地震获取相关的数据。具体的，震源和检波器一般会设置于海水面以下一定的深度处。具体实施时，由于海水与空气波阻抗差异较大，两者之间的接触面可以视为一个比较良好的反射面。因此，震源信号到达该界面的一次波会被再次反射回地下，进而接收器除了接收到通常的一次反射波外，还将接收到与海平面反射作用相关的伴随波，这种伴随波称为虚反射，也称为鬼波。由于虚反射与一次波极性相反，并且虚反射与一次波的时差很小，导致在获取海底地震数据的过程中，虚反射会叠加在一次波的后部，会影响所采集的海洋地震数据的准确度。进而会导致，应用这种带有虚反射的海洋地震数据进行相关施工时，会引起地震记录的低频响应，降低地震剖面的分辨率，给地震资料反演与解释造成很大的困难。因此，如何对采集的海洋地震数据中的虚反射进行压制，一直是人们关注的一个问题。目前，为了压制海洋地震数据中的虚反射，通常采用时间域的反褶积方法进行虚反射压制处理。上述方法是基于地震记录是反射波与虚反射传播算子褶积的假设，在时间域上利用预测反褶积来压制虚反射。但是，带有虚反射的地震数据往往并不是最小相位，不满足反褶积的条件，因此，上述通过时间域反褶积进行虚反射压制的方法并不能完全消除虚反射。此外，现有的虚反射压制方法往往只能对通过水平缆采集的海洋地震数据进行相应的虚反射压制。对于通过其他方式采集的海洋地震数据，例如，通过斜缆采集的海洋地震数据，并不完全适用。综上可知，现有的虚反射的压制方法，具体实施时，往往存在不能完全消除虚反射、压制效果差、应用范围有限的技术问题。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本申请实施方式提供了一种虚反射的压制方法和装置，以解决现有的虚反射的压制方法中存在的无法完全消除虚反射、压制效果差、应用范围有限的技术问题。本申请实施方式提供了一种虚反射的压制方法，包括：获取目标区域的海洋地震数据，其中，所述海洋地震数据包括多个单炮地震数据；根据所述单炮地震数据，通过波场反向延拓，确定反向延拓数据；根据所述反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱；对所述虚反射压制算子的相位谱进行希尔伯特变换，得到虚反射压制算子的振幅谱；根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，确定虚反射压制算子；利用所述虚反射压制算子，对所述海洋地震数据进行虚反射压制处理，得到压制虚反射的海洋地震数据。在一个实施方式中，所述获取目标区域的海洋地震数据，包括：通过水平缆获取所述海洋地震数据，和/或，通过斜缆获取所述海洋地震数据。在一个实施方式中，根据所述单炮地震数据，通过波场反向延拓，确定反向延拓数据，包括：对所述单炮地震数据进行时间参量和空间参量的二维傅里叶变换，得到频率域和波数域的单炮地震数据；从所述海洋地震数据中获取道头信息，并根据所述道头信息，确定检波器的深度；根据所述检波器深度和所述频率域和波数域的单炮地震数据，进行波场反向延拓，得到所述反向延拓数据。在一个实施方式中，根据所述检波器深度和所述频率域和波数域的单炮地震数据，进行波场反向延拓，得到所述反向延拓数据，包括：根据所述检波器深度和所述频率域和波数域的单炮地震数据，按照以下公式，确定所述反向延拓数据：上式中，D1(ω,kx)为反向延拓数据，D(ω,kx)为频率域和波数域的单炮地震数据，z是检波器的深度，c为地震波在海水中的传播速度，ω为角频率，kx为x方向的波数，x为炮检距。在一个实施方式中，根据所述反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱，包括：对所述反向延拓数据进行二维傅里叶反变换，得到延拓后的时间域和空间域的单炮地震数据；按照预设时窗长度，将所述单炮地震数据划分为多个单时窗的单炮单道地震数据，将所述延拓后的时间域和空间域的单炮地震数据划分为多个单时窗的延拓后的时间域和空间域单炮单道地震数据；将所述单时窗的单炮单道地震数据进行一维傅里叶变换，得到单时窗的频率域单炮单道地震数据，将所述单时窗的延拓后的时间域和空间域的单炮单道地震数据进行一维傅里叶变换，得到单时窗的延拓后的频率域单炮单道地震数据；根据所述单时窗的频率域的单炮单道地震数据和所述单时窗的延拓后的频率域单炮单道地震数据，确定所述虚反射压制算子的相位谱。在一个实施方式中，根据所述单时窗的频率域的单炮单道地震数据和所述单时窗的延拓后的频率域的单炮单道地震数据，确定所述虚反射压制算子的相位谱，包括：根据所述单时窗的频率域的单炮单道地震数据和所述单时窗的延拓后的频率域的单炮单道地震数据，按照以下公式，确定所述反射压制算子的相位谱：上式中，为虚反射压制算子的相位谱，d1(f,x)为单时窗的频率域的单炮单道地震数据，d(f,x)为单时窗的延拓后的频率域的单炮单道地震数据，f为频率。在一个实施方式中，根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，确定虚反射压制算子，包括：根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，按照以下公式，确定所述虚反射压制算子：上式中，H(f,x)为虚反射压制算子，A(f,x)为虚反射压制算子的振幅谱，为虚反射压制算子的相位谱。在一个实施方式中，利用所述虚反射压制算子，对所述海洋地震数据进行虚反射压制处理，得到压制虚反射的海洋地震数据，包括：按照以下方式对所述海洋地震数据中的多个单炮地震数据中的同一炮的地震数据分别进行虚反射压制处理，得到压制虚反射的时间域的单炮地震数据：利用所述虚反射压制算子，对所述单炮地震数据的单时窗的各个频率域的单炮单道地震数据分别进行虚反射压制处理，得到压制虚反射的单时窗的频率域的单炮单道地震数据；对所述压制虚反射的单时窗的频率域的单炮单道地震数据进行一维傅里叶反变换，得到压制虚反射的单时窗的时间域单炮单道地震数据。在一个实施方式中，利用所述虚反射压制算子，对所述单炮地震数据的单时窗的频率域的单炮单道地震数据分别进行虚反射压制处理，得到压制虚反射的单时窗的频率域的单炮单道地震数据，包括：利用所述虚反射压制算子，按照以下公式，确定所述压制虚反射的单时窗的频率域的单炮单道地震数据：上式中，P(f,x)为压制虚反射的单时窗的频率域的单炮单道地震数据，d(f,x)为单时窗的频率域的单炮单道地震数据，H(f,x)为虚反射压制算子。在一个实施方式中，在得到压制虚反射的海洋地震数据后，所述方法还包括：根据所述压制虚反射的海洋地震数据，对所述目标区域进行海底石油地震勘探。本申请还提供了一种虚反射的压制装置，包括：获取模块，用于获取目标区域的海洋地震数据，其中，所述海洋地震数据包括多个单炮地震数据；反向延拓模块，用于根据所述单炮地震数据，通过波场反向延拓，确定反向延拓数据；第一确定模块，用于根据所述反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱；第二确定模块，用于对所述虚反射压制算子的相位谱进行希尔伯特变换，确定虚反射压制算子的振幅谱；第三确定模块，用于根据所述虚反射压制算子本文档来自技高网...

【技术保护点】
一种虚反射的压制方法，其特征在于，包括：获取目标区域的海洋地震数据，其中，所述海洋地震数据包括多个单炮地震数据；根据所述单炮地震数据，通过波场反向延拓，确定反向延拓数据；根据所述反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱；对所述虚反射压制算子的相位谱进行希尔伯特变换，得到虚反射压制算子的振幅谱；根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，确定虚反射压制算子；利用所述虚反射压制算子，对所述海洋地震数据进行虚反射压制处理，得到压制虚反射的海洋地震数据。

【技术特征摘要】
1.一种虚反射的压制方法，其特征在于，包括：获取目标区域的海洋地震数据，其中，所述海洋地震数据包括多个单炮地震数据；根据所述单炮地震数据，通过波场反向延拓，确定反向延拓数据；根据所述反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱；对所述虚反射压制算子的相位谱进行希尔伯特变换，得到虚反射压制算子的振幅谱；根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，确定虚反射压制算子；利用所述虚反射压制算子，对所述海洋地震数据进行虚反射压制处理，得到压制虚反射的海洋地震数据。2.根据权利要求1所述的方法，其特征在于，所述获取目标区域的海洋地震数据，包括：通过水平缆获取所述海洋地震数据，和/或，通过斜缆获取所述海洋地震数据。3.根据权利要求1所述的方法，其特征在于，根据所述单炮地震数据，通过波场反向延拓，确定反向延拓数据，包括：对所述单炮地震数据进行时间参量和空间参量的二维傅里叶变换，得到频率域和波数域的单炮地震数据；从所述海洋地震数据中获取道头信息，并根据所述道头信息，确定检波器的深度；根据所述检波器深度和所述频率域和波数域的单炮地震数据，进行波场反向延拓，得到所述反向延拓数据。4.根据权利要求3所述的方法，其特征在于，根据所述检波器深度和所述频率域和波数域的单炮地震数据，进行波场反向延拓，得到所述反向延拓数据，包括：根据所述检波器深度和所述频率域和波数域的单炮地震数据，按照以下公式，确定所述反向延拓数据：上式中，D1(ω,kx)为反向延拓数据，D(ω,kx)为频率域和波数域的单炮地震数据，z是检波器的深度，c为地震波在海水中的传播速度，ω为角频率，kx为x方向的波数，x为炮检距。5.根据权利要求1所述的方法，其特征在于，根据所述反向延拓数据和所述单炮地震数据，确定虚反射压制算子的相位谱，包括：对所述反向延拓数据进行二维傅里叶反变换，得到延拓后的时间域和空间域的单炮地震数据；按照预设时窗长度，将所述单炮地震数据划分为多个单时窗的单炮单道地震数据，将所述延拓后的时间域和空间域的单炮地震数据划分为多个单时窗的延拓后的时间域和空间域的单炮单道地震数据；将所述单时窗的单炮单道地震数据进行一维傅里叶变换，得到单时窗的频率域的单炮单道地震数据，将所述单时窗的延拓后的时间域和空间域的单炮单道地震数据进行一维傅里叶变换，得到单时窗的延拓后的频率域的单炮单道地震数据；根据所述单时窗的频率域的单炮单道地震数据和所述单时窗的延拓后的频率域的单炮单道地震数据，确定所述虚反射压制算子的相位谱。6.根据权利要求5所述的方法，其特征在于，根据所述单时窗的频率域的单炮单道地震数据和所述单时窗的延拓后的频率域的单炮单道地震数据，确定所述虚反射压制算子的相位谱，包括：根据所述单时窗的频率域的单炮单道地震数据和所述单时窗的延拓后的频率域的单炮单道地震数据，按照以下公式，确定所述反射压制算子的相位谱：上式中，为虚反射压制算子的相位谱，d1(f,x)为单时窗的频率域的单炮单道地震数据，d(f,x)为单时窗的延拓后的频率域的单炮单道地震数据，f为频率。7.根据权利要求1所述的方法，其特征在于，根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，确定虚反射压制算子，包括：根据所述虚反射压制算子的相位谱和所述虚反射压制算子的振幅谱，按照以下公式，确定所述虚反射压制算子：上式中，...

【专利技术属性】
技术研发人员：马光凯，周铮铮，钱忠平，任晓乔，方云峰，尹天奎，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11