近地表Q场生成方法、装置、设备及介质与流程

本技术涉及地震资料处理，特别涉及一种近地表q场生成方法、装置、设备及介质。

背景技术：

1、在油气勘探领域，对勘探目标的精度要求越来越高，这就要求地震数据频带宽、分辨率高。然而地震波在地层中传播时，能量被地层吸收衰减，尤其是近地表的介质疏松，对地震波的吸收衰减会严重降低地震数据的频带宽和分辨率。目前，可以通过生成近地表q场，基于近地表q场对地震数据进行补偿来提高地震数据的频带宽和分辨率。

2、相关技术中，生成近地表q场的方式为：获取双井微测井数据，采用谱比法或频移法，计算井上多个单点的q值，并基于双井微测井数据解释得到每个单点的速度，拟合出q值与速度的关系曲线，基于q值与速度的关系曲线，将近地表速度模型转换成近地表q场。

3、由于拟合的q值与速度的关系曲线存在一定误差，因此，生成的近地表q场的准确性较差。

技术实现思路

1、本技术实施例提供了一种近地表q场生成方法、装置、设备及介质，提高了近地表q场的准确性。所述技术方案如下：

2、一方面，提供了一种近地表q场生成方法，所述方法包括：

3、从目标区域的地震数据中截取多个地震道的初至波，一个地震道对应一个炮点和一个检波点，所述地震道的初至波为所述地震道的炮点激发地震之后，所述地震道的检波点检测到的第一个地震波，包括多个采集点对应的振幅；

4、分别对每个地震道的初至波中的多个采集点对应的振幅进行统计处理，得到所述每个地震道对应的均方根振幅；

5、将所述每个地震道对应的均方根振幅进行分解处理，得到所述每个地震道对应的炮点振幅分量和检波点振幅分量，所述炮点振幅分量用于表示近地表对下行地震波的影响，所述检波点振幅分量用于表示近地表对上行反射波的影响；

6、对于所述每个地震道，基于所述地震道的炮点振幅分量，确定所述地震道对应的炮点的相对吸收衰减系数，基于所述地震道的检波点振幅分量，确定所述地震道对应的检波点的相对吸收衰减系数，所述相对吸收衰减系数用于表示近地表对地震波的吸收衰减情况；

7、基于所述目标区域的近地表速度模型，确定每个炮点对应的近地表垂直单程传播时间和每个检波点对应的近地表垂直单程传播时间，所述近地表垂直单程传播时间为地震波在近地表的上下表面之间垂直传播一次所需的时间；

8、基于所述每个炮点的相对吸收衰减系数、所述每个炮点对应的近地表垂直单程传播时间、所述每个检波点的相对吸收衰减系数以及所述每个检波点的对应的近地表垂直单程传播时间，生成所述目标区域的近地表q场；

9、基于所述目标区域的近地表q场进行处理。

10、在一种可能的实现方式中，所述基于所述每个炮点的相对吸收衰减系数、所述每个炮点对应的近地表垂直单程传播时间、所述每个检波点的相对吸收衰减系数以及所述每个检波点的对应的近地表垂直单程传播时间，生成所述目标区域的近地表q场，包括：

11、基于所述每个炮点的相对吸收衰减系数、所述每个炮点对应的近地表垂直单程传播时间和所述地震数据的主频，确定所述每个炮点对应的q值；

12、基于所述每个检波点的相对吸收衰减系数、所述每个检波点对应的近地表垂直单程传播时间和所述地震数据的主频，确定所述每个检波点对应的q值；

13、基于所述每个炮点对应的q值和所述每个检波点对应的q值，生成所述目标区域的近地表q场。

14、在一种可能的实现方式中，所述相对吸收衰减系数、所述近地表垂直单程传播时间、所述地震数据的主频以及所述q值之间的关系如下：

15、q＝-πft/ln(r)；

16、其中，q为炮点对应的q值或检波点对应的q值，为圆周率，f为地震数据的主频，t为炮点对应的近地表垂直单程传播时间或检波点对应的近地表垂直单程传播时间，ln为自然对数函数，r为炮点的相对吸收衰减系数或检波点的相对吸收衰减系数。

17、在一种可能的实现方式中，所述从目标区域的地震数据中截取多个地震道的初至波，包括：

18、从所述目标区域的地震数据中拾取多个地震道的初至时间，所述地震道的初至时间为所述地震道的炮点激发地震之后，所述地震道的检波点检测到第一个地震波的时间；

19、基于所述多个地震道的初至时间，从所述目标区域的地震数据中截取多个地震道的初至波。

20、在一种可能的实现方式中，所述基于所述初至时间，从所述目标区域的地震数据中截取多个地震道的初至波，包括：

21、显示初至波截取界面，所述初至波截取界面包括所述地震数据、所述多个地震道的初至时间以及第一时长输入框和第二时长输入框；

22、基于所述第一时长输入框，获取输入的第一时长，基于所述第二时长输入框，获取输入的第二时长；

23、对于任一地震道，将所述地震道的初至时间减去第一时长得到的第三时间作为所述地震道的初至波截取时窗的起始时间，将所述地震道的初至时间增加所述第二时长得到的第四时间作为所述地震道的初至波截取时窗的终止时间；

24、基于所述多个地震道的初至波截取时窗，从所述目标区域的地震数据中截取所述多个地震道的初至波。

25、在一种可能的实现方式中，所述基于所述地震道的炮点振幅分量，确定所述地震道对应的炮点的相对吸收衰减系数，基于所述地震道的检波点振幅分量，确定所述地震道对应的检波点的相对吸收衰减系数，包括：

26、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行归一化处理，得到所述地震道对应的炮点的相对吸收衰减系数和对应的检波点的相对吸收衰减系数。

27、在一种可能的实现方式中，所述将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行归一化处理，得到所述地震道对应的炮点的相对吸收衰减系数和对应的检波点的相对吸收衰减系数之前，所述方法还包括以下至少一项：

28、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行异常值剔除处理；

29、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行插值处理；

30、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行平滑处理。

31、另一方面，提供了一种近地表q场生成装置，所述装置包括：

32、截取模块，用于从目标区域的地震数据中截取多个地震道的初至波，一个地震道对应一个炮点和一个检波点，所述地震道的初至波为所述地震道的炮点激发地震之后，所述地震道的检波点检测到的第一个地震波，包括多个采集点对应的振幅；

33、统计模块，用于分别对每个地震道的初至波中的多个采集点对应的振幅进行统计处理，得到所述每个地震道对应的均方根振幅；

34、分解模块，用于将所述每个地震道对应的均方根振幅进行分解处理，得到所述每个地震道对应的炮点振幅分量和检波点振幅分量，所述炮点振幅分量用于表示近地表对下行地震波的影响，所述检波点振幅分量用于表示近地表对上行反射波的影响；

35、系数确定模块，用于对于所述每个地震道，基于所述地震道的炮点振幅分量，确定所述地震道对应的炮点的相对吸收衰减系数，基于所述地震道的检波点振幅分量，确定所述地震道对应的检波点的相对吸收衰减系数，所述相对吸收衰减系数用于表示近地表对地震波的吸收衰减情况；

36、时间确定模块，用于基于所述目标区域的近地表速度模型，确定每个炮点对应的近地表垂直单程传播时间和每个检波点对应的近地表垂直单程传播时间，所述近地表垂直单程传播时间为地震波在近地表的上下表面之间垂直传播一次所需的时间；

37、生成模块，用于基于所述每个炮点的相对吸收衰减系数、所述每个炮点对应的近地表垂直单程传播时间、所述每个检波点的相对吸收衰减系数以及所述每个检波点的对应的近地表垂直单程传播时间，生成所述目标区域的近地表q场；

38、处理模块，用于基于所述目标区域的近地表q场进行处理。

39、在一种可能的实现方式中，所述生成模块，用于基于所述每个炮点的相对吸收衰减系数、所述每个炮点对应的近地表垂直单程传播时间和所述地震数据的主频，确定所述每个炮点对应的q值；基于所述每个检波点的相对吸收衰减系数、所述每个检波点对应的近地表垂直单程传播时间和所述地震数据的主频，确定所述每个检波点对应的q值；基于所述每个炮点对应的q值和所述每个检波点对应的q值，生成所述目标区域的近地表q场。

40、在一种可能的实现方式中，所述相对吸收衰减系数、所述近地表垂直单程传播时间、所述地震数据的主频以及所述q值之间的关系如下：

41、q＝-πft/ln(r)；

42、其中，q为炮点对应的q值或检波点对应的q值，π为圆周率，f为地震数据的主频，t为炮点对应的近地表垂直单程传播时间或检波点对应的近地表垂直单程传播时间，ln为自然对数函数，r为炮点的相对吸收衰减系数或检波点的相对吸收衰减系数。

43、在一种可能的实现方式中，所述截取模块，包括：

44、时间确定单元，用于从所述目标区域的地震数据中拾取多个地震道的初至时间，所述地震道的初至时间为所述地震道的炮点激发地震之后，所述地震道的检波点检测到第一个地震波的时间；

45、截取单元，用于基于所述多个地震道的初至时间，从所述目标区域的地震数据中截取多个地震道的初至波。

46、在一种可能的实现方式中，所述截取单元，用于显示初至波截取界面，所述初至波截取界面包括所述地震数据、所述多个地震道的初至时间以及第一时长输入框和第二时长输入框；基于所述第一时长输入框，获取输入的第一时长，基于所述第二时长输入框，获取输入的第二时长；对于任一地震道，将所述地震道的初至时间减去第一时长得到的第三时间作为所述地震道的初至波截取时窗的起始时间，将所述地震道的初至时间增加所述第二时长得到的第四时间作为所述地震道的初至波截取时窗的终止时间；基于所述多个地震道的初至波截取时窗，从所述目标区域的地震数据中截取所述多个地震道的初至波。

47、在一种可能的实现方式中，所述系数确定模块，用于将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行归一化处理，得到所述地震道对应的炮点的相对吸收衰减系数和对应的检波点的相对吸收衰减系数。

48、在一种可能的实现方式中，所述装置还包括预处理模块，所述预处理模块，用于执行以下至少一项：

49、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行异常值剔除处理；

50、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行插值处理；

51、将所述地震道的炮点振幅分量和所述地震道的检波点振幅分量进行平滑处理。

52、另一方面，提供了一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述至少一条程序代码由所述处理器加载并执行，以实现如上述任一实现方式所述的近地表q场生成方法。

53、另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储至少一条程序代码，所述至少一条程序代码由处理器加载并执行，以实现如上述任一实现方式所述的近地表q场生成方法。

54、另一方面，提供了一种计算机程序产品，所述计算机程序产品包括至少一条程序代码，所述至少一条程序代码由处理器加载并执行，以实现如上述任一实现方式所述的近地表q场生成方法。

55、本技术实施例提供的技术方案的有益效果至少包括：

56、本技术实施例提供了一种近地表q场生成方法，通过分解初至波的均方根振幅，得到炮点和检波点的振幅分量，进而基于炮点和检波点的振幅分量确定炮点和检波点的振幅分量，确定炮点和检波点的相对吸收衰减系数，基于炮点和检波点的相对吸收衰减系数和地震波在近地表的垂直单程传播时间，生成近地表q场。由于初至波未经地下界面反射，是由炮点激发经近地表后到达检波点，因此，初至波主要受到近地表的影响，从而初至波包含丰富的近地表吸收衰减信息，保证了炮点和检波点的相对吸收衰减系数的准确性。由于地震波的近地表垂直单程传播时间也收到近地表地质影响，因此，综合炮点和检波点的相对吸收衰减系数以及地震波在近地表中的垂直单程传播时间，能够更加准确地得到近地表q场，也即是，提高了近地表q场的准确性。