转换波中水层多次波的获取方法及装置与流程

本发明属于海洋地震数据处理，具体涉及一种转换波中水层多次波的获取方法及装置。

背景技术：

1、海洋地震数据处理中多次波是比较突出的问题，它会覆盖有效信号，在偏移成像中形成构造假象，因此多次波压制是地震资料处理非常关键的一步，通过压制多次波可以有效提高速度分析精度、较少成像假象，大幅改善海洋地震资料成像质量。多次波技术经过几十年的发展，纵波多次波衰减技术已经比较成熟。其中一些技术因计算速度快，计算效率高、压制效果好等优点已经广泛应用于实际成产，比如预测反褶积技术、srmr技术和mwd技术等。

2、转换波可以简单的理解为纵波入射，横波出射，其射线路径区别于纵波波场，是几何不对称的，因此横波和转换波在地层间的传播规律远比纵波复杂的多。我们知道海洋地震纵波波场中包含了大量的水层多次波和自由表面多次波，转换波中同样包含了大量的水层多次波和自由表面多次波，由于海面的反射系数显著强于包括海底在内的任一波阻抗界面，因此转换波中的多次波主要是以水层多次波为主。

3、转换波中水层多次波的存在会严重影响地震资料信噪比，且会对转换波成像产生影响，造成成像假象。因此发明一种转换波波场中水层多次波的压制方法对于提高转换波的成像质量非常必要。

技术实现思路

1、鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的一种转换波中水层多次波的获取方法及装置、计算设备及计算机存储介质。

2、根据本发明的第一个方面，提供了一种转换波中水层多次波的获取方法，包括：

3、根据海水深度建立海水深度模型，以及采集地震数据，将所述地震数据按照检波点、炮点的分选方式抽取成共检波点道集；

4、对于任意一个目标炮点，以所述目标炮点和检波点作为一个炮检对，根据所述深度模型和地震数据计算出所述炮检对之间任意一个所述中间炮点和所述目标炮点之间波场延拓的延拓值；

5、对于所述炮检对之间的任意一个所述中间炮点，根据该中间炮点与目标炮点之间的延拓值，以及所述共检波点道集中该中间炮点的地震数据，计算得到经过该中间炮点反射之后的水层多次波；

6、将所述炮检对之间的每一个所述中间炮点对应的所述水层多次波进行叠加，得到所述炮检对对应的水层多次波。

7、进一步的，所述对于任意一个目标炮点，以所述目标炮点和检波点作为一个炮检对，根据所述深度模型和地震数据计算出所述炮检对之间任意一个所述中间炮点和所述目标炮点之间波场延拓的延拓值进一步包括：

8、根据纯纵波情况下的基尔霍夫积分定理，推导得到波场延拓的延拓算子；

9、根据所述深度模型、所述地震数据以及所述延拓算子计算出所述炮检对之间任意一个所述中间炮点和所述目标炮点之间波场延拓的延拓值。

10、进一步的，所述根据纯纵波情况下的基尔霍夫积分定理，推导得到波场延拓的延拓算子进一步包括：

11、根据纯纵波情况下的基尔霍夫积分定理，将闭合曲面内任意一个目标点的压力场表示为两项之和，其中第一项与所述闭合曲面上的压力有关，第二项与所述闭合曲面上的质点运动速度有关；

12、通过中间函数使第二项等于0，则可以得到所述闭合曲面内任意一点的压力场与所述闭合曲面上的压力的第一表达式；

13、通过利用平面和无穷远界面构成所述闭合曲面，在无穷远界面的扰动尚未到达所述目标点的前提下，将所述第一表达式转化为第二表达式，在所述第二表达式中，所述闭合曲面内任意一点的压力场与所述平面和所述闭合曲面上的压力有关；

14、根据所述第二表达式，得到波场延拓的延拓算子。

15、进一步的，所述对于所述炮检对之间的任意一个所述中间炮点，根据该中间炮点与目标炮点之间的延拓值，以及所述共检波点道集中该中间炮点的地震数据，计算得到经过该中间炮点反射之后的水层多次波具体为：

16、对于所述炮检对之间的任意一个所述中间炮点，根据经过该中间炮点反射之后的水层多次波等于该中间炮点与目标炮点之间的延拓值与所述共检波点道集中该中间炮点的地震数据的乘积，计算得到经过该中间炮点反射之后的水层多次波。

17、进一步的，在采集地震数据之后，所述方法还包括：对所述地震数据进行插值处理，使相邻两个数据点之间的间隔小于等于预设的最大采样间隔，所述预设的最大采样间隔为地震波最小波长的二分之一；

18、所述将所述地震数据按照检波点、炮点的分选方式抽取成共检波点道集具体为：将插值处理后的地震数据按照检波点、炮点的分选方式抽取成共检波点道集。

19、进一步的，在所述对于所述炮检对之间的任意一个所述中间炮点，根据该中间炮点与目标炮点之间的延拓值，以及所述共检波点道集中该中间炮点的地震数据，计算得到经过该中间炮点反射之后的水层多次波之后，所述方法还包括：对所述水层多次波进行滤波处理；

20、所述将所述炮检对之间的每一个所述中间炮点对应的所述水层多次波进行叠加，得到所述炮检对对应的水层多次波具体为：将滤波处理后得到的水层多次波进行叠加，得到所述炮检对对应的水层多次。

21、根据本发明的第二发明，提供了一种转换波中水层多次波的获取装置，包括：

22、数据采集及海水深度模型建立模块，用于根据海水深度建立海水深度模型，以及采集地震数据，将所述地震数据按照检波点、炮点的分选方式抽取成共检波点道集；

23、第一计算模块，对于任意一个目标炮点，以所述目标炮点和检波点作为一个炮检对，用于根据所述深度模型和地震数据计算出所述炮检对之间任意一个所述中间炮点和所述目标炮点之间波场延拓的延拓值；

24、第二计算模块，对于所述炮检对之间的任意一个所述中间炮点，用于根据该中间炮点与目标炮点之间的延拓值，以及所述共检波点道集中该中间炮点的地震数据，计算得到经过该中间炮点反射之后的水层多次波；

25、水层多次波获取模块，用于将所述炮检对之间的每一个所述中间炮点对应的所述水层多次波进行叠加，得到所述炮检对对应的水层多次波。

26、根据本发明的第三方面，提供了一种计算设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；

27、所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行以上任一项所述的转换波中水层多次波的获取方法对应的操作。

28、根据本发明的第四方面，提供了一种计算机存储介质，所述存储介质中存储有至少一可执行指令，所述可执行指令使处理器执行以上任一项所述的转换波中水层多次波的获取方法对应的操作。

29、由上述技术方案可知，本发明提供的一种转换波中水层多次波的获取方法及装置，具有如下有益效果：

30、本发明转换波中水层多次波的获取方法利用波场延拓能够得到炮点端的水层多次波，从而可以用于压制检波点接收的转换波中的水层多次波，提高转换波成像的质量。

31、上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。