深反射地震数据最优频带选取方法、电子设备及存储介质与流程

1.本技术涉及深反射地震数据处理领域，尤其涉及深反射地震数据最优频带选取方法、电子设备及存储介质。


背景技术：

2.海洋海盆区岩石圈壳幔结构成像是研究区域构造活动的重要依据，目前，海上地震数据莫霍面成像普遍存在成像难的问题，根据地震波传播原理，当地震波从软的沉积层传播到硬的火成岩基底时，由于两种介质巨大的速度与密度属性差异，会产生强烈的反射现象，仅有少数低频能量能透过界面，向地壳深处传播。
3.由于莫霍面的存在，地震波的主体能量均以反射波的形式，从该界面反射至检波点，同时，由于南海海盆区通常沉积层较薄，导致其地震波高频成分衰减不充分。
4.目前的深反射地震数据处理，还是采用全频段处理的方式，没有全面压制高频能量，无法凸显基底以下的低频信息。深反射地震数据处理通常无法确定目的层位置，而且基底以下存在大量噪声，信噪比远低于常规地震资料。
5.因此，目前亟需一种深反射地震数据最优频带选取方法。


技术实现要素：

6.为克服相关技术中存在的问题，本技术提供深反射地震数据最优频带选取方法、电子设备及存储介质，该方法能够精准获取深反射地震数据的最优频带，从而凸显深反射信息。
7.本技术第一方面提供深反射地震数据最优频带选取方法，包括：
8.将深反射地震数据重新排序，将共炮点道集数据转换为共检波点道集数据，形成单道剖面；
9.通过设置滤波器频带对所述单道剖面进行分频滤波，得到滤波剖面，根据所述滤波剖面的信噪比和对比度选取初始优势频带；
10.扩展所述初始优势频带的中间频率，得到扩展优势频带，所述扩展优势频带的范围为最低扩展频率到最高扩展频率；根据所述最低扩展频率和所述初始优势频带的最低频率得到低频优势频带，根据所述最高扩展频率和所述初始优势频带的最高频率得到高频优势频带；根据所述低频优势频带和所述高频优势频带求取最优频带。
11.通过设置滤波器频带对所述单道剖面进行分频滤波，包括：
12.设置6个所述滤波器频带，所述滤波器频带分别设置为1-6hz、6-12hz、12-24hz、24-36hz、36-60hz和60hz以上；
13.分别使用设置的6个所述滤波器频带对所述单道剖面进行分频滤波。
14.所述扩展所述初始优势频带的中间频率，包括：
15.计算所述初始优势频带的中间频率，再将所述中间频率向低频率和高频率的方向各自进行扩展；
16.所述中间频率为所述初始优势频带的最低频率和所述初始优势频带的最高频率的平均值。
17.所述再将所述中间频率向低频率和高频率的方向各自进行扩展，包括：
18.将所述中间频率减去修正频率，得到最低扩展频率；
19.将所述中间频率增加所述修正频率，得到最高扩展频率；
20.将所述最低扩展频率作为所述扩展优势频带的最低频率，将所述最高扩展频率作为所述扩展优势频带的最高频率；
21.所述修正频率为所述中间频率的一半。
22.所述根据所述最低扩展频率和所述初始优势频带的最低频率得到低频优势频带，包括：
23.将所述最低扩展频率作为所述低频优势频带的最高频率，将所述初始优势频带的最低频率作为所述低频优势频带的最低频率。
24.所述根据所述最高扩展频率和所述初始优势频带的最高频率得到高频优势频带，包括：
25.将所述最高扩展频率作为所述高频优势频带的最低频率，将所述初始优势频带的最高频率作为所述高频优势频带的最高频率。
26.所述根据所述低频优势频带和所述高频优势频带求取最优频带，包括：
27.将所述低频优势频带的中间频率作为所述最优频带的最低频率，将所述高频优势频带的中间频率作为所述最优频带的最高频率；
28.所述低频优势频带的中间频率为所述低频优势频带的最低频率和所述低频优势频带的最高频率的平均值；
29.所述高频优势频带的中间频率为所述高频优势频带的最低频率和所述高频优势频带的最高频率的平均值。
30.所述根据所述滤波剖面的信噪比和对比度选取初始优势频带，包括：
31.对所述单道剖面进行分频滤波，得到滤波剖面；
32.选取初始优势频带，所述初始优势频带为信噪比和对比度最高的所述滤波剖面对应的所述滤波器频带。
33.本技术第二方面提供一种电子设备，包括：
34.处理器；以及
35.存储器，其上存储有可执行代码，当所述可执行代码被所述处理器执行时，使所述处理器执行如上所述的方法。
36.本技术第三方面提供一种非暂时性机器可读存储介质，其上存储有可执行代码，当所述可执行代码被电子设备的处理器执行时，使所述处理器执行如上所述的方法。
37.本技术提供的技术方案可以包括以下有益效果：
38.本技术技术方案通过设置滤波器频带对单道剖面进行分频滤波，根据滤波剖面的信噪比和对比度选取初始优势频带，将滤波器频带的范围进行缩小。然后扩展初始优势频带的中间频率，得到扩展优势频带。根据最低扩展频率和初始优势频带的最低频率得到低频优势频带。低频优势频带包含较多的目标层深反射信息，但是图像信噪比较低。根据最高扩展频率和初始优势频带的最高频率得到高频优势频带，高频优势频带包含的目标层深反
射信息较少，但是图像信噪比较高。再根据低频优势频带和高频优势频带求取最优频带，结合低频优势频带包含较多的目标层深反射信息和高频优势频带图像信噪比较高的优点，兼顾深反射地震数据的可靠性和图像的可观测性。
39.应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本技术。
附图说明
40.通过结合附图对本技术示例性实施方式进行更详细的描述，本技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本技术示例性实施方式中，相同的参考标号通常代表相同部件。
41.图1是本技术实施例示出的深反射地震数据最优频带选取方法的流程示意图；
42.图2是本技术实施例示出的滤波剖面的示意图；
43.图3是本技术实施例示出的滤波剖面的另一示意图；
44.图4是本技术实施例示出的迭代滤波器频带后的单道剖面的示意图；
45.图5是本技术实施例示出的原始叠加剖面的示意图；
46.图6是本技术实施例示出的优势频带叠加剖面的示意图；
47.图7是本技术实施例示出的电子设备的结构示意图。
具体实施方式
48.下面将参照附图更详细地描述本技术的优选实施方式。虽然附图中显示了本技术的优选实施方式，然而应该理解，可以以各种形式实现本技术而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本技术更加透彻和完整，并且能够将本技术的范围完整地传达给本领域的技术人员。
49.在本技术使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本技术。在本技术和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
50.应当理解，尽管在本技术可能采用术语“第一”、“第二”、“第三”等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本技术范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本技术的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
51.实施例一
52.目前深反射地震数据处理还是采用全频段处理的方式，没有全面压制高频能量，无法凸显基底以下的低频信息。
53.针对上述问题，本技术实施例提供深反射地震数据最优频带选取方法，以下结合附图详细描述本技术实施例的技术方案。
54.图1是本技术实施例示出的深反射地震数据最优频带选取方法的流程示意图，请
参阅图1，深反射地震数据最优频带选取方法具体如下：
55.a1、数据分选，得到单道剖面。
56.将深反射地震数据重新排序，由共炮点道集数据转换为共检波点道集数据，形成单道剖面。
57.a2、设置滤波器频带对单道剖面进行分频滤波。
58.设置滤波器频带：将滤波器频带分为1-6hz、6-12hz、12-24hz、24-36hz、36-60hz及60hz以上6个不同的滤波器频带，使用6个不同的滤波器频带对单道剖面进行分频滤波，得到滤波剖面。
59.a3、根据滤波剖面的信噪比和对比度选取初始优势频带。
60.选择信噪比和对比度最高的滤波剖面对应的滤波器频带作为初始优势频带。
61.示例性，经过对比发现，使用1-6hz、6-12hz、12-24hz和24-36hz频带得到的滤波剖面的效果较好，基底以下层位隐约可见潜在的深反射信号，将基底以下层位作为目标层位，选取1-36hz为初始优势频带。
62.优选地，对多个滤波剖面做能量均衡调整，调节滤波剖面的色彩比重，突出目的层。
63.图2是本技术实施例示出的滤波剖面的示意图，包含了1-6hz、6-12hz和12-24hz分别对应的滤波剖面的示意图。图3是本技术实施例示出的滤波剖面的另一示意图，包含了24-36hz、36-60hz和60hz以上分别对应的滤波剖面的示意图。
64.a4、扩展初始优势频带的中间频率，得到扩展优势频带。
65.计算初始优势频带的最低频率和初始优势频带的最高频率的平均值，得到中间频率。将中间频率减去修正频率，得到最低扩展频率。将中间频率加上修正频率，得到最高扩展频率。所述修正频率为所述中间频率的一半。
66.将最低扩展频率作为扩展优势频带的最低频率，将最高扩展频率作为扩展优势频带的最高频率。
67.示例性，初始优势频带为1-36hz，中间频率为18.5hz，修正频率为9.25hz。最低扩展频率为18.5-9.25＝9.25hz，与最低扩展频率相对应的，最高扩展频率为18.5 9.25＝27.75hz。扩展优势频带为9.25-27.75hz。
68.a5、求取低频优势频带和高频优势频带。
69.所述低频优势频带为包含较多的目标层深反射信息，但是图像信噪比较低的频带。
70.所述高频优势频带为包含的目标层深反射信息较少，但是图像信噪比较高的频带。
71.将所述最低扩展频率作为所述低频优势频带的最高频率，将所述初始优势频带的最低频率作为所述低频优势频带的最低频率。
72.将所述最高扩展频率作为所述高频优势频带的最低频率，将所述初始优势频带的最高频率作为所述高频优势频带的最高频率。
73.示例性，低频优势频带的最高频率为9.25hz，低频优势频带的最低频率为1hz，低频优势频带为1-9.25hz，高频优势频带的最低频率为27.75hz，高频优势频带的最高频率为36hz，高频优势频带为27.75-36hz。
74.a6、根据低频优势频带和高频优势频带求取最优频带。
75.最优频带为包含较多的目标层深反射信息，并且信噪比较高的频带。最优频带同时具备低频优势频带和高频优势频带的优点。
76.求取低频优势频带的最低频率和低频优势频带的最高频率的平均值，将该平均值作为低频优势频带的中间频率。
77.求取高频优势频带的最低频率和高频优势频带的最高频率的平均值，将该平均值作为高频优势频带的中间频率。
78.将低频优势频带的中间频率作为所述最优频带的最低频率，将高频优势频带的中间频率作为所述最优频带的最高频率。
79.示例性，低频优势频带的中间频率为(1 9.25)/2＝5.125hz，高频优势频带的中间频率为(27.75 36)/2＝31.875hz，最优频带为5.125-31.875hz。
80.图4是本技术实施例示出的迭代滤波器频带后的单道剖面的示意图，请参阅图4，图4左为使用初始优势频带进行滤波得到的滤波剖面示意图，图4中为使用扩展优势频带进行滤波得到的滤波剖面示意图，图4右为使用最优频带进行滤波得到的滤波剖面示意图。
81.a7、将最优频带对应的滤波剖面转换为叠加剖面。
82.叠加剖面为将单道剖面的横向序列打散，单道剖面纵向上时间序列相对不变，然后将打散后的时间序列重组和叠加得到的剖面。
83.使用最优频带对单道剖面进行滤波，得到最优频带对应的滤波剖面，将滤波剖面的横向序列打散，滤波剖面纵向上时间序列相对不变，然后将打散后的时间序列重组和叠加，得到优势频带叠加剖面。
84.滤波剖面是根据滤波器频带进行滤波后的单道剖面。
85.图5是本技术实施例示出的原始叠加剖面的示意图，图6是本技术实施例示出的优势频带叠加剖面的示意图。
86.通过对比图5和图6可以看出，优势频带叠加剖面比原始叠加剖面具有更多的基底以下的低频信息。
87.本技术实施例通过设置滤波器频带对单道剖面进行分频滤波，根据滤波剖面的信噪比和对比度选取初始优势频带，将滤波器频带的范围进行缩小。然后扩展初始优势频带的中间频率，得到扩展优势频带。根据最低扩展频率和初始优势频带的最低频率得到低频优势频带。低频优势频带包含较多的目标层深反射信息，但是图像信噪比较低。根据最高扩展频率和初始优势频带的最高频率得到高频优势频带，高频优势频带包含的目标层深反射信息较少，但是图像信噪比较高。再根据低频优势频带和高频优势频带求取最优频带，结合低频优势频带包含较多的目标层深反射信息和高频优势频带图像信噪比较高的优点，兼顾深反射地震数据的可靠性和图像的可观测性。
88.实施例二
89.与前述应用功能实现方法实施例相对应，本技术还提供了一种电子设备及相应的实施例。
90.图7是本技术实施例示出的电子设备的结构示意图。
91.参见图7，电子设备1000包括存储器1010和处理器1020。
92.处理器1020可以是中央处理单元(centralprocessingunit，cpu)，还可以是其他
通用处理器、数字信号处理器(digitalsignalprocessor，dsp)、专用集成电路(applicationspecificintegratedcircuit，asic)、现场可编程门阵列(field-programmablegatearray，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
93.存储器1010可以包括各种类型的存储单元，例如系统内存、只读存储器(rom)，和永久存储装置。其中，rom可以存储处理器1020或者计算机的其他模块需要的静态数据或者指令。永久存储装置可以是可读写的存储装置。永久存储装置可以是即使计算机断电后也不会失去存储的指令和数据的非易失性存储设备。在一些实施方式中，永久性存储装置采用大容量存储装置(例如磁或光盘、闪存)作为永久存储装置。另外一些实施方式中，永久性存储装置可以是可移除的存储设备(例如软盘、光驱)。系统内存可以是可读写存储设备或者易失性可读写存储设备，例如动态随机访问内存。系统内存可以存储一些或者所有处理器在运行时需要的指令和数据。此外，存储器1010可以包括任意计算机可读存储媒介的组合，包括各种类型的半导体存储芯片(dram，sram，sdram，闪存，可编程只读存储器)，磁盘和/或光盘也可以采用。在一些实施方式中，存储器1010可以包括可读和/或写的可移除的存储设备，例如激光唱片(cd)、只读数字多功能光盘(例如dvd-rom，双层dvd-rom)、只读蓝光光盘、超密度光盘、闪存卡(例如sd卡、minsd卡、micro-sd卡等等)、磁性软盘等等。计算机可读存储媒介不包含载波和通过无线或有线传输的瞬间电子信号。
94.存储器1010上存储有可执行代码，当可执行代码被处理器1020处理时，可以使处理器1020执行上文述及的方法中的部分或全部。
95.上文中已经参考附图详细描述了本技术的方案。在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详细描述的部分，可以参见其他实施例的相关描述。本领域技术人员也应该知悉，说明书中所涉及的动作和模块并不一定是本技术所必须的。另外，可以理解，本技术实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减，本技术实施例装置中的模块可以根据实际需要进行合并、划分和删减。
96.此外，根据本技术的方法还可以实现为一种计算机程序或计算机程序产品，该计算机程序或计算机程序产品包括用于执行本技术的上述方法中部分或全部步骤的计算机程序代码指令。
97.或者，本技术还可以实施为一种非暂时性机器可读存储介质(或计算机可读存储介质、或机器可读存储介质)，其上存储有可执行代码(或计算机程序、或计算机指令代码)，当所述可执行代码(或计算机程序、或计算机指令代码)被电子设备(或电子设备、服务器等)的处理器执行时，使所述处理器执行根据本技术的上述方法的各个步骤的部分或全部。
98.本领域技术人员还将明白的是，结合这里的申请所描述的各种示例性逻辑块、模块、电路和算法步骤可以被实现为电子硬件、计算机软件或两者的组合。
99.附图中的流程图和框图显示了根据本技术的多个实施例的系统和方法的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标记的功能也可以以不同于附图中所标记的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流
程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
100.以上已经描述了本技术的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。