一种基于像素叠加的反射点拾取方法、装置、设备及介质与流程

本发明涉及地球物理勘探，尤其涉及一种基于像素叠加的反射点拾取方法、装置、电子设备及计算机可读存储介质。

背景技术：

1、本部分旨在为权利要求书中陈述的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。

2、目前，面向深度域速度建模的偏移速度分析依旧以层析反演技术为主，该技术包含三大要素：层析反射点拾取、层析核函数构建以及层析正则化技术。其中层析反射点拾取是层析反演的数据基础，也是层析反演技术能否应用成功的关键。层析核函数是层析反演核心算法，其计算精度很大程度上依托于层析反射点的拾取。层析正则化技术属于层析反演的实用化配套技术，主要是对计算精度的补充和完善。整体而言，从实际应用层面来看，层析反射点拾取精度及自动化程度是能否获得高精度速度模型的关键，拾取质量也是层析反演对于陆上数据工业化应用的瓶颈。目前层析反射点拾取方法有多种，主流的自动拾取方法均是基于振幅的同相叠加的思想，简单来讲，反射点与构造倾角的拾取为一个扫描线性叠加的过程，即在反射点处沿不同方向进行加窗叠加，叠加能量最大的方向即认为是同相轴的走向，对剖面中的每一点均实施这一操作，然后根据预先给定的门槛值筛选出潜在的反射点位置。这种方法的精度十分依赖于地震子波的形态与同相轴曲率，因此其精度取决于各类控制参数的选取，且容易受到剖面中的强振幅的影响，造成地层倾角估计的偏差。这样对于同一剖面采用同样的参数进行地层倾角估计，会出现部分地层倾角得到准确估计，同时存在很多地方的倾角则受到强振幅的影响，造成反射点拾取出现大的不稳定性，不适应低信噪比数据，对弱同相轴不敏感。

3、基于此，本专利通过发明一种局部像素叠加扫描的算法，通过对局部时窗进行扫描叠加，对强像素区域进行有效信号拾取，该技术利用扫描叠加的局部像素差异进行反射信号识别，能够适应弱信号或低信噪比区的反射点拾取，对于深层及山前带等低信噪比区的反射点拾取精度均有较大程度的提升。

技术实现思路

1、针对上述问题，本发明的实施例提供了一种基于像素叠加的反射点拾取方法、装置、电子设备及计算机可读存储介质。

2、第一方面，本发明实施例提供了一种基于像素叠加的反射点拾取方法，包括：

3、逐个选取预设的地震数据体中的地震数据组作为目标数据组，逐个选取所述目标数据组中的数据像素点位作为目标像素点位，利用预设的计算窗口对所述目标像素点位进行局部像素叠加，得到所述目标像素点位的数据像素叠加值；

4、利用所述计算窗口计算出所述目标像素点位的时窗像素值，根据所述时窗像素值和所述数据像素叠加值对所述目标像素点位进行数据像素扫描，并将所有的目标像素点位的扫描结果汇集成所述目标数据组的数据像素扫描结果；

5、利用预设的拾取窗口对所述数据像素扫描结果进行层析极值拾取，得到所述目标数据组的组数据拾取结果，并将所有的组数据拾取结果汇集成数据体拾取结果，完成层析反射点拾取。

6、根据本发明的实施例，所述利用预设的计算窗口对所述目标像素点位进行局部像素叠加，得到所述目标像素点位的数据像素叠加值，包括：

7、从所述计算窗口中分别提取出横轴半时窗长度、竖轴半时窗长度以及纵轴半时窗长度；

8、从所述目标数据组中选取出所述目标像素点位对应的地震数据值；

9、利用如下的像素叠加值公式根据所述地震数据值、所述横轴半时窗长度、所述竖轴半时窗长度以及所述纵轴半时窗长度计算出所述数据像素叠加值：

10、

11、其中，pixellocal(x,y,z)表示所述目标像素点位为(x,y,z)时的所述数据像素叠加值，x是指所述目标像素点位的横轴坐标，y是指所述目标像素点位的竖轴坐标，z是指所述目标像素点位的纵轴坐标，w1是指所述横轴半时窗长度，w2是指所述竖轴半时窗长度，w3是指所述纵轴半时窗长度，由所述横轴半时窗长度、所述竖轴半时窗长度以及所述纵轴半时窗长度生成的计算窗口大小为(2w1+1)(2w2+1)(2w3+1)，|·|是绝对值运算符，p(i,j,k)是指所述数据像素点位为(i,j,k)时对应的所述地震数据值，其中，i的取值范围在(x-w1,x+w1)之间，j的取值范围在(y-w2,y+w2)之间，k的取值范围在(z-w3,z+w3)之间。

12、根据本发明的实施例，所述从所述计算窗口中分别提取出横轴半时窗长度、竖轴半时窗长度以及纵轴半时窗长度，包括：

13、将所述计算窗口的横轴窗口长度减一，得到横轴叠加长度，并将所述横轴叠加长度除以二，得到所述横轴半时窗长度；

14、将所述计算窗口的竖轴窗口长度减一，得到竖轴叠加长度，并将所述竖轴叠加长度除以二，得到所述竖轴半时窗长度；

15、将所述计算窗口的纵轴窗口长度减一，得到纵轴叠加长度，并将所述纵轴叠加长度除以二，得到所述纵轴半时窗长度。

16、根据本发明的实施例，所述利用所述计算窗口计算出所述目标像素点位的时窗像素值，包括：

17、利用如下的时窗像素值算法根据所述目标像素点位的地震数据值、所述横轴半时窗长度、所述竖轴半时窗长度以及所述纵轴半时窗长度计算出所述时窗像素值：

18、

19、其中，是指当所述目标像素点位为(x,y,z)时的所述时窗像素值，x是指所述目标像素点位的横轴坐标，y是指所述目标像素点位的竖轴坐标，z是指所述目标像素点位的纵轴坐标，w1是指所述横轴半时窗长度，w2是指所述竖轴半时窗长度，w3是指所述纵轴半时窗长度，由所述横轴半时窗长度、所述竖轴半时窗长度以及所述纵轴半时窗长度生成的计算窗口大小为(2w1+1)(2w2+1)(2w3+1)，|·|是绝对值运算符，p(i,j,k)是指所述数据像素点位为(i,j,k)时对应的所述地震数据值，其中，i的取值范围在(x-w1,x+w1)之间，j的取值范围在(y-w2,y+w2)之间，k的取值范围在(z-w3,z+w3)之间。

20、根据本发明的实施例，所述根据所述时窗像素值和所述数据像素叠加值对所述目标像素点位进行数据像素扫描，包括：

21、利用如下的像素扫描公式根据所述时窗像素值和所述数据像素叠加值对所述目标像素点位进行数据像素扫描：

22、

23、其中，δpixel(x,y,z)是指所述目标像素点位的坐标为(x,y,z)时的所述扫描结果，x是指所述目标像素点位的横轴坐标，y是指所述目标像素点位的竖轴坐标，z是指所述目标像素点位的纵轴坐标，pixellocal(x,y,z)表示当所述目标像素点位为(x,y,z)时的所述数据像素叠加值，是指当所述目标像素点位为(x,y,z)时的所述时窗像素值。

24、根据本发明的实施例，所述利用预设的拾取窗口对所述数据像素扫描结果进行层析极值拾取，得到所述目标数据组的组数据拾取结果，包括：

25、从所述拾取窗口中分别提取出横轴半拾取窗长、竖轴半拾取窗长以及纵轴半拾取窗长；

26、逐个选取所述目标数据组中的数据像素点位作为拾取像素点位，从所述数据扫描结果中提取出所述拾取像素点位的扫描结果作为待提取扫描结果；

27、利用如下的层析拾取公式根据所述待提取扫描结果、所述横轴半拾取窗长、所述竖轴半拾取窗长以及所述纵轴半拾取窗长计算出所述拾取像素点位的拾取结果，并将所有的所述拾取像素点位的拾取结果汇集成所述目标数据组的组数据拾取结果：

28、

29、其中，ppick是指所述拾取像素点位为时的拾取结果，是指所述拾取像素点位的横轴坐标，是指所述拾取像素点位的竖轴坐标，是指所述拾取像素点位的纵轴坐标，max是取极大值符号，m1是指所述横轴半拾取窗长，m2是指所述竖轴半拾取窗长，m3是指所述纵轴半拾取窗长，是指所述拾取像素点位为所述拾取窗口大小为(2m1+1)(2m2+1)(2m3+1)时，所述拾取窗口内部的所述数据像素点位为时的所述扫描结果，且的取值范围在之间，的取值范围在之间，的取值范围在之间。

30、根据本发明的实施例，所述从所述拾取窗口中分别提取出横轴半拾取窗长、竖轴半拾取窗长以及纵轴半拾取窗长，包括：

31、将所述拾取窗口的横轴窗口长度减一，得到横轴拾取长度，并将所述横轴拾取长度除以二，得到所述横轴半拾取窗长；

32、将所述拾取窗口的竖轴窗口长度减一，得到竖轴拾取长度，并将所述竖轴拾取长度除以二，得到所述竖轴半拾取窗长；

33、将所述拾取窗口的纵轴窗口长度减一，得到纵轴拾取长度，并将所述纵轴拾取长度除以二，得到所述纵轴半拾取窗长。

34、第二方面，本发明实施例提供了一种基于像素叠加的反射点拾取装置，其特征在于，包括：

35、像素叠加模块，用于逐个选取预设的地震数据体中的地震数据组作为目标数据组，逐个选取所述目标数据组中的数据像素点位作为目标像素点位，利用预设的计算窗口对所述目标像素点位进行局部像素叠加，得到所述目标像素点位的数据像素叠加值；

36、时窗扫描模块，用于利用所述计算窗口计算出所述目标像素点位的时窗像素值，根据所述时窗像素值和所述数据像素叠加值对所述目标像素点位进行数据像素扫描，并将所有的目标像素点位的扫描结果汇集成所述目标数据组的数据像素扫描结果；

37、极值拾取模块，用于利用预设的拾取窗口对所述数据像素扫描结果进行层析极值拾取，得到所述目标数据组的组数据拾取结果，并将所有的组数据拾取结果汇集成数据体拾取结果，完成层析反射点拾取。

38、第三方面，本发明实施例提供了一种电子设备，其包括：

39、处理器；

40、用于存储所述处理器可执行指令的存储器；

41、其中，所述处理器被配置为执行所述指令，以实现如前第一方面所述的一种基于像素叠加的反射点拾取方法。

42、第四方面，本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时，实现如前第一方面所述的一种基于像素叠加的反射点拾取方法。

43、与现有技术相比，本发明的上述技术方案具有如下有益效果：

44、本发明的实施例通过利用预设的计算窗口对所述目标像素点位进行局部像素叠加，得到所述目标像素点位的数据像素叠加值，能够有效地从所述目标数据组中提取出信号特征，实现对弱信号或低信噪比地区的信号增强，从而提高了后续反射点拾取的拾取精度，通过利用所述像素扫描公式根据所述时窗像素值和所述数据像素叠加值对所述目标像素点位进行数据像素扫描，并将所有的目标像素点位的扫描结果汇集成所述目标数据组的数据像素扫描结果，能够计算出所述数据像素叠加值和所述时窗像素值之间的差值，从而更好的反应所述地震数据体中数据的波动程度，能够有效过滤所述地震数据体中的杂波数据，实现特征提取的效果，增强了地震数据体的纹理特征，从而提高反射点拾取的精度，通过利用预设的拾取窗口对所述数据像素扫描结果进行层析极值拾取，得到所述目标数据组的组数据拾取结果，能够对所述目标数据组对应的所述数据像素扫描结果进行极值扫描，从而增强所述地震数据中波形信号纹路的对比度，通过将所有的组数据拾取结果汇集成数据体拾取结果，完成层析反射点拾取，能够实现对整个地震数据体的局部像素增强，从而提升拾取精度，因此，本发明提出的像素叠加的反射点拾取方法、装置、电子设备及计算机可读存储介质，可以解决进行反射点拾取时精度较低的问题。