信号反褶积方法、装置、计算机设备和存储介质与流程

1.本技术涉及地震资料处理技术领域，特别涉及一种信号反褶积方法、装置、计算机设备和存储介质。


背景技术：

2.在通过信号反褶积方法处理地震数据时，该方法能够有效处理地震数据中每个地震道的地震子波数据，得到尖脉冲信号和反射系数；再基于尖脉冲信号和反射系数可以得到高分辨率的地震数据。
3.相关技术中，信号反褶积方法主要采用一维信号反褶积方法来处理地震子波数据，一维信号反褶积方法仅处理了在垂直方向上的地震子波数据。由于不同的地震子波在激发方向上会有方位角和出射角上的差异，这样，通过一维信号反褶积方法处理多个不同方位角和出射角上的地震子波数据，会降低信号反褶积方法处理地震子波数据的准确度。


技术实现要素：

4.本技术实施例提供了一种信号反褶积方法、装置、计算机设备和存储介质，能够提高信号反褶积方法处理地震数据的准确度。所述技术方案如下：
5.一方面，提供了一种信号反褶积方法，所述方法包括：
6.确定地震道集初始的第一地震数据，所述地震道集包括多个地震道；
7.基于所述第一地震数据，确定所述地震道集的第二地震数据，以及从所述第一地震数据中确定所述多个地震道的方位角和出射角，所述第二地震数据包括每个地震道对应的射线参数；
8.对于每个地震道，基于所述地震道的方位角和出射角，确定所述地震道对应的角度组合，所述角度组合包括所述方位角和所述出射角；
9.对于多个预设角度组合中的每个预设角度组合，基于所述多个地震道对应的角度组合，确定目标地震道，所述目标地震道的角度组合属于所述预设角度组合；
10.对于每个预设角度组合，叠加所述预设角度组合对应的多个目标地震道的射线参数，得到所述预设角度组合对应的单元子波数据；
11.基于多个预设角度组合对应的单元子波数据，确定所述地震道集的远场子波数据；
12.基于所述远场子波数据，生成所述地震道集的反褶积算子；
13.通过所述反褶积算子，处理所述第一地震数据，得到所述地震道集的反褶积地震数据。
14.在一种可能的实现方式中，所述基于所述远场子波数据，生成所述地震道集的反褶积算子，包括：
15.对所述远场子波数据进行气泡消除处理、压制鬼波处理和零相位化处理，得到期望子波数据；
16.确定将所述远场子波数据转换为所述期望子波数据所使用的所述反褶积算子。
17.在一种可能的实现方式中，所述确定地震道集的第一地震数据，包括：
18.获取所述地震道集的原始地震数据；
19.对所述原始地震数据进行预处理，得到第三地震数据；
20.对所述第三地震数据进行τ-p域正变换处理，得到所述地震道集的第一地震数据。
21.在一种可能的实现方式中，所述对所述原始地震数据进行预处理，得到第三地震数据，包括：
22.对所述原始地震数据进行规则化处理和插值处理，得到第四地震数据；
23.从所述第四地震数据中确定目标地震数据，所述目标地震数据为波形发生预设波动的数据；
24.确定所述目标地震数据中的下行波数据；
25.对所述下行波数据中的多次波数据进行压制，得到所述第三地震数据。
26.在一种可能的实现方式中，所述确定所述目标节点数据中的下行波数据，包括：
27.对所述目标地震数据进行波场分离处理，得到所述目标地震数据中的下行波数据。
28.在一种可能的实现方式中，所述基于所述第一地震数据，确定所述地震道集的第二地震数据，包括：
29.对所述第一地震数据进行噪音压制处理，得到压制噪音后的第一地震数据；
30.对所述压制噪音后的第一地震数据进行同相轴拉平处理，得到同相轴拉平后的第一地震数据；
31.对所述同相轴拉平后的第一地震数据进行畸变切除处理，得到切除畸变后的第一地震数据；
32.对所述切除畸变后的第一地震数据进行边界外推处理，得到所述地震道集的第二地震数据。
33.在一种可能的实现方式中，所述通过所述反褶积算子，处理所述第一地震数据，得到所述地震道集的反褶积地震数据，包括：
34.通过所述反褶积算子，处理所述第一地震数据，得到第五地震数据；
35.对所述第五地震数据进行τ-p域反变换处理，得到所述反褶积地震数据。
36.另一方面，提供了一种信号反褶积装置，所述装置包括：
37.第一确定模块，用于确定地震道集初始的第一地震数据，所述地震道集包括多个地震道；
38.第二确定模块，用于基于所述第一地震数据，确定所述地震道集的第二地震数据，以及从所述第一地震数据中确定所述多个地震道的方位角和出射角，所述第二地震数据包括每个地震道对应的射线参数；
39.第三确定模块，用于对于每个地震道，基于所述地震道的方位角和出射角，确定所述地震道对应的角度组合，所述角度组合包括所述方位角和所述出射角；
40.第四确定模块，用于对于多个预设角度组合中的每个预设角度组合，基于所述多个地震道对应的角度组合，确定目标地震道，所述目标地震道的角度组合属于所述预设角度组合；
41.叠加模块，用于对于每个预设角度组合，叠加所述预设角度组合对应的多个目标地震道的射线参数，得到所述预设角度组合对应的单元子波数据；
42.第五确定模块，用于基于多个预设角度组合对应的单元子波数据，确定所述地震道集的远场子波数据；
43.生成模块，用于基于所述远场子波数据，生成所述地震道集的反褶积算子；
44.处理模块，用于通过所述反褶积算子，处理所述第一地震数据，得到所述地震道集的反褶积地震数据。
45.在一种可能的实现方式中，所述生成模块，用于：
46.对所述远场子波数据进行气泡消除处理、压制鬼波处理和零相位化处理，得到期望子波数据；
47.确定将所述远场子波数据转换为所述期望子波数据所使用的所述反褶积算子。
48.在一种可能的实现方式中，所述第一确定模块，包括：
49.获取单元，用于获取所述地震道集的原始地震数据；
50.第一处理单元，用于对所述原始地震数据进行预处理，得到第三地震数据；
51.第二处理单元，用于对所述第三地震数据进行τ-p域正变换处理，得到所述地震道集的第一地震数据。
52.在一种可能的实现方式中，所述第一处理单元，包括：
53.处理子单元，用于对所述原始地震数据进行规则化处理和插值处理，得到第四地震数据；
54.第一确定子单元，用于从所述第四地震数据中确定目标地震数据，所述目标地震数据为波形发生预设波动的数据；
55.第二确定子单元，用于确定所述目标地震数据中的下行波数据；
56.压制单元，用于对所述下行波数据中的多次波数据进行压制，得到所述第三地震数据。
57.在一种可能的实现方式中，所述第二确定子单元，用于：
58.对所述目标地震数据进行波场分离处理，得到所述目标地震数据中的下行波数据。
59.在一种可能的实现方式中，所述第二确定模块，用于：
60.对所述第一地震数据进行噪音压制处理，得到压制噪音后的第一地震数据；
61.对所述压制噪音后的第一地震数据进行同相轴拉平处理，得到同相轴拉平后的第一地震数据；
62.对所述同相轴拉平后的第一地震数据进行畸变切除处理，得到切除畸变后的第一地震数据；
63.对所述切除畸变后的第一地震数据进行边界外推处理，得到所述地震道集的第二地震数据。
64.在一种可能的实现方式中，所述处理模块，用于：
65.通过所述反褶积算子，处理所述第一地震数据，得到第五地震数据；
66.对所述第五地震数据进行τ-p域反变换处理，得到所述反褶积地震数据。
67.另一方面，提供了一种计算机设备，所述计算机设备包括一个或多个处理器和一
个或多个存储器，所述一个或多个存储器中存储有至少一条指令，所述至少一条指令由所述一个或多个处理器加载并执行以实现上述任一实现方式所述的信号反褶积方法所执行的操作。
68.另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储至少一条指令，所述至少一条指令由处理器加载并执行，以实现上述任一实现方式所述的信号反褶积方法所执行的操作。
69.另一方面，提供了一种计算机程序产品或计算机程序，所述计算机程序产品或所述计算机程序包括计算机程序代码，所述计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取所述计算机程序代码，处理器执行所述计算机程序代码，使得所述计算机设备执行上述的信号反褶积方法所执行的操作。
70.本技术实施例提供的技术方案的有益效果至少包括：
71.本技术实施例提供了一种信号反褶积方法，由于该方法基于地震道的方位角和出射角对地震道进行了划分，这样，再基于划分方位角和出射角的地震道来依次确定单元子波数据、远场子波数据，使得确定的远场子波数据充分考虑了方位角和出射角的因素；进而使得基于远场子波数据确定的反褶积算子与方位角和出射角相关，这样，基于反褶积算子处理第一地震数据时，能够有效补偿方位角和出射角带来的差异，从而使得基于反褶积算子确定的反褶积地震数据更加准确。
附图说明
72.为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
73.图1是本技术实施例提供的一种信号反褶积方法的流程图；
74.图2是本技术实施例提供的一种地震道集的示意图；
75.图3是本技术实施例提供的一种地震道集的示意图；
76.图4是本技术实施例提供的一种地震道集的示意图；
77.图5是本技术实施例提供的一种射线参数的划分图；
78.图6是本技术实施例提供的一种单元子波的示意图；
79.图7是本技术实施例提供的一种远场子波的示意图；
80.图8是本技术实施例提供的一种远场子波的频谱图；
81.图9是本技术实施例提供的一种期望子波的频谱图；
82.图10是本技术实施例提供的一种远场子波的相位谱；
83.图11是本技术实施例提供的一种期望子波的相位谱；
84.图12是本技术实施例提供的一种期望子波的示意图；
85.图13是本技术实施例提供的一种地震道集的示意图；
86.图14是本技术实施例提供的一种地震道集的示意图；
87.图15是本技术实施例提供的一种地震子波的示意图；
88.图16是本技术实施例提供的一种地震子波的叠加图；
89.图17是本技术实施例提供的一种地震数据的频谱图；
90.图18是本技术实施例提供的一种地震数据的剖面图；
91.图19是本技术实施例提供的一种信号反褶积装置的框图；
92.图20是本技术实施例提供的一种计算机设备的框图。
具体实施方式
93.为使本技术的目的、技术方案和优点更加清楚，下面将结合附图对本技术实施方式作进一步地详细描述。
94.本技术的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们的任意变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
95.本技术实施例提供一种信号反褶积方法，参见图1，方法包括：
96.步骤101：计算机设备确定地震道集初始的第一地震数据。
97.其中，该地震道集包括多个地震道，该地震道集为待研究的目标区域的任同一检波器接收的所有地震道形成的道集。目标区域为待处理地震数据的区域，可以为海洋区域。该目标区域包括多个检波器，即该目标区域的地震数据包括多个地震道集的第一地震数据。该目标区域的地震数据可以为气枪震源产生的地震波信号在目标区域传播后经检波器检收得到的数据。
98.该步骤可以通过以下步骤(1)-(3)实现：
99.(1)计算机设备获取该地震道集的原始地震数据。
100.其中，原始地震数据包括该地震道集的多个地震道的旅行时间和横向旅行路程。
101.(2)计算机设备对原始地震数据进行预处理，得到第三地震数据。
102.该步骤可以通过以下步骤a1-a4实现：
103.a1：计算机设备对原始地震数据进行规则化处理和插值处理，得到第四地震数据。
104.需要说明的是，原始地震数据采集过程中由于存在障碍物、禁区、海洋拖缆羽状漂移和经济成本等因素，使得原始地震数据沿空间方向通常是不规则或稀疏采样的，且由于道缺失和剔除坏道等也会造成地震数据不规则，不规则的地震数据处理过程中会引入空间假频；因此，需要对原始地震数据进行规则化处理，使能够得到一个规则且完整的地震道集；以及，需要对原始地震数据进行插值处理，以填补原始地震数据中缺失的数据和加大原始地震数据的密度，使能够得到完整的地震数据以及提高原始地震数据的精密度。这样，通过规则化处理和插值处理，能够避免原始地震数据进行τ-p域正变换后出现的空间假频和其它假象。
105.参见图2，图2的左边图为原始地震数据经τ-p域正变换后的地震道集，存在着空间假频现象。图2的右边图为对原始地震数据进行规则化处理和插值处理后，得到的第四地震数据经τ-p域正变换的地震道集；从图中可以看出，规则化处理和插值处理后，空间假频现象得到了有效消除。
106.a2：计算机设备从第四地震数据中确定目标地震数据，目标地震数据为波形发生
预设波动的数据。
107.其中，发生预设波动的目标地震数据为第四地震数据中噪音弱、气泡能量较强的品质好的数据。
108.a3：计算机设备确定目标地震数据中的下行波数据。
109.在该步骤中，计算机设备对目标地震数据进行波场分离处理，得到目标地震数据中的下行波数据。其中，计算机设备基于原始地震数据中的水检数据和陆检数据对目标地震数据进行交叉鬼波化处理，得到分离后的下行波数据。
110.在本技术实施例中，由于下行波数据的气泡效应相对更强，这样，通过确定下行波数据，能够使进行τ-p域正变换的第三地震数据更加合理。
111.a4：计算机设备对下行波数据中的多次波数据进行压制，得到第三地震数据。
112.其中，压制多次波数据的方法可以为共中心点叠加法、拉东变换法等。
113.需要说明的是，地震道的地震子波的长度对信号反褶积方法的处理效果影响较大，地震子波越长，其包含的气泡信息越丰富，进而使得信号反褶积方法压制的气泡信息的效果越好。然而多次波会制约地震子波的长度，因此，将多次波数据进行压制后，通过信号反褶积方法能够更好的压制第一地震数据中的气泡信息。
114.(3)计算机设备对第三地震数据进行τ-p域正变换处理，得到地震道集的第一地震数据。
115.参见图3，图3的左边图为第三地震数据的地震道集，纵坐标为旅行时间，横坐标为横向旅行路程；图3的右边图为进行τ-p域正变换得到的第一地震数据的地震道集，纵坐标为地震道的截距，横坐标为射线参数。
116.需要说明的是，由于地震道集为三维的地震道集，包括横向的地震道集和纵向的地震道集，即进行的τ-p域正变换为三维τ-p域正变换，即第一地震数据位于三维的τ-p
x-py域内。
117.在本技术实施例中，通过对第三地震数据进行三维τ-p域正变换处理，实现了将时空域的第三地震数据转换到了三维τ-p域，信号反褶积大部分工作在三维τ-p域进行，三维τ-p域变换是信号反褶积的基础工作，三维τ-p域变换是信号反褶积成败的关键，进而便于对三维τ-p域内的第一地震数据进行处理。
118.步骤102：计算机设备基于第一地震数据，确定地震道集的第二地震数据，以及从第一地震数据中确定多个地震道的方位角和出射角，
119.其中，第二地震数据包括每个地震道对应的射线参数，每个地震道的射线参数包括如图3的右边图所示的每个截距对应的射线参数。
120.计算机设备确定地震道集的第二地震数据通过以下步骤(1)-(4)实现：
121.(1)计算机设备对第一地震数据进行噪音压制处理，得到压制噪音后的第一地震数据。
122.需要说明的是，在地震勘探中，检波器在采集来自地下地层的有效波数据的同时也记录了来自地上和地下的各种各样的与地下地层无关的干扰波噪音；干扰波噪音混淆了有效波，使得难以提取到有效的地震数据；第一地震数据中包括线性干扰波噪音、异常振幅干扰波噪音等干扰波噪音。在本技术实施例中，可以采用三维f-k(频率-波数)滤波方法压制线性干扰波噪音，采用中值滤波法压制异常振幅噪音。在本技术实施例中，将第一地震数
据进行噪音压制处理，能够有效突出有效波数据，进而能够得到有效的第一地震数据。
123.(2)计算机设备对压制噪音后的第一地震数据进行同相轴拉平处理，得到同相轴拉平后的第一地震数据。
124.需要说明的是，地震道集的同相轴不平会影响后续对地震道集的处理，因此，需要将对第一地震数据进行同相轴拉平处理，将多个地震道的轴拉到同一时间线上，以得到同相轴拉平的地震道集。其中，若同时处理多个第一地震数据，则可以将多个第一地震数据对应的地震道集间进行同相轴拉平。
125.(3)计算机设备对同相轴拉平后的第一地震数据进行畸变切除处理，得到切除畸变后的第一地震数据。
126.需要说明的是，第一地震数据中会存在影响有效波数据的边界信息和深层信息等畸变数据，因此，需要对第一地震数据进行畸变切除处理，以得到有效的第一地震数据。
127.(4)计算机设备对切除畸变后的第一地震数据进行边界外推处理，得到地震道集的第二地震数据。
128.其中，计算机设备对切除畸变后的第一地震数据的地震道集沿着切除边界向两侧进行外推处理，进而实现了将第一地震数据中的畸变数据的有效去除。
129.参见图4，图4的左边图为τ-p域正变换得到的第一地震数据的地震道集，图4的右边图为第一地震数据经过上述处理后得到的第二地震数据的地震道集；可以看出，第二地震数据的地震道集规则且平滑，与气泡信息无关的噪音数据、畸变数据等信息被去除，得到了气泡信息丰富的地震道集。多个地震道的方位角和出射角可以从如图4的右边图所示的地震道集中每个地震道的角度数据中提取得到。
130.步骤103：计算机设备对于每个地震道，基于地震道的方位角和出射角，确定地震道对应的角度组合。
131.需要说明的是，地震道的地震子波在发射方向上会有方位角和出射角的差异，而不同的方位角和出射角对地震道的影响不同，因此，在处理地震数据时，需要充分考虑方位角和出射角的因素。
132.其中，角度组合包括地震道的方位角和出射角。方位角的角度范围为0
°‑
360
°
，出射角的角度范围为0
°‑
90
°
，角度组合为地震道的方位角和出射角组合得到的；例如，若地震道的方位角为35
°
，出射角为70
°
，则地震道的角度组合为35
°×
70
°
；若地震道的方位角为30.2
°
，出射角为25.9
°
，则该地震道的角度组合为30.2
°×
25.9
°
。
133.步骤104：计算机设备对于多个预设角度组合中的每个预设角度组合，基于多个地震道对应的角度组合，确定目标地震道。
134.其中，目标地震道的角度组合属于该预设角度组合。
135.多个预设角度组合包括360
°
范围内的方位角和90
°
范围内的出射角。多个预设角度组合包括0
°×0°
、0
°×1°
、0
°×2°
、1
°×0°
、2
°×0°
、1
°×1°…
等预设角度组合。每个预设角度组合包括一个角度范围内的多个方位角和一个角度范围内的多个出射角；例如，对于预设角度组合1
°×1°
，则该预设角度组合包括[1,2)区间内的方位角和[1,2)区间内的出射角，如1.25的方位角和1.36的出射角等；角度组合为1.1
°×
1.5
°
、1.38
°×
1.69
°
等的地震道为该预设角度组合的目标地震道。
[0136]
参见图5，图5的上半部分为依次排序的多个方位角对应的地震道的射线参数的划
分图，图5的下半部分为依次排序的多个出射角对应的地震道的射线参数的划分图。从图中可以看出，每个角度范围内的方位角和出射角均包括多个地震道的射线参数。
[0137]
在本技术实施例中，通过基于地震道的方位角和出射角对地震道进行划分，充分考虑了出射角和方位角对地震道带来的影响，进而能够提高后续对地震道处理的准确度。
[0138]
步骤105：计算机设备对于每个预设角度组合，叠加该预设角度组合对应的多个目标地震道的射线参数，得到该预设角度组合对应的单元子波数据。
[0139]
其中，每个预设角度组合对应一个单元子波；图6为多个预设角度组合对应的多个单元子波数据生成的多个单元子波的示意图，从图中可以看出，多个单元子波规则且有序排列。
[0140]
步骤106：计算机设备基于多个预设角度组合对应的单元子波数据，确定地震道集的远场子波数据。
[0141]
继续参见图6，计算机设备将如图6所示的多个单元子波的单元子波数据进行叠加，得到一个远场子波数据，该远场子波数据对应的远场子波如图7所示。
[0142]
步骤107：计算机设备基于远场子波数据，生成地震道集的反褶积算子。
[0143]
该步骤可以通过以下步骤(1)-(2)实现：
[0144]
(1)计算机设备对远场子波数据进行气泡消除处理、压制鬼波处理和零相位化处理，得到期望子波数据。
[0145]
其中，计算机设备可以采用预测反褶积方法进行气泡消除处理，采用模拟炮点虚反射方法来进行压制鬼波处理，以及进行相位调整来实现零相位化处理。
[0146]
图8为远场子波的频谱图，从图中可以看出，低频和高频部分均存在抖动，低频的抖动代表气泡的影响，高频的抖动代表鬼波的影响，说明远场子波存在气泡效应和鬼波影响。图9为期望子波的频谱图，从图中可以看出，气泡得到了有效的消除，鬼波得到了有效的压制。图10为远场子波的相位谱，从图中可以看出，远场子波的相位偏移，并没有归于零位。图11为期望子波的相位谱，从图中可以看出，期望子波的相位归于零位。图12为期望子波数据对应的期望子波，从图中可以看出，该期望子波的气泡得到了消除，鬼波得到了压制，相位为零相位，为待得到的理想的地震子波。
[0147]
(2)计算机设备确定将远场子波数据转换为期望子波数据所使用的反褶积算子。
[0148]
需要说明的是，远场子波数据和期望子波数据存在一定差距，通过对远场子波数据和期望子波数据进行相匹配处理，得到能够使远场子波数据到达期望子波数据的因子，即反褶积算子。下面以公式一为例进行说明。
[0149]
公式一：y＝a*x。
[0150]
其中，y代表期望子波数据，x代表远场子波数据，a代表反褶积算子。则反褶积算子可以通过期望子波数据与远场子波数据相除，得到反褶积算子，即实现了相匹配处理，得到了反褶积算子。
[0151]
在本技术实施例中，由于基于地震道的方位角和出射角对地震道进行了划分，这样，再基于划分方位角和出射角的地震道来依次确定单元子波数据、远场子波数据，使得确定的远场子波数据充分考虑了方位角和出射角的因素，这样，使得基于远场子波数据得到的反褶积算子与方位角和出射角有关，使该反褶积算子有效补偿了地震道的地震子波在方位角和出射角上的差异。并且，由于确定反褶积算子的期望子波为进行了气泡消除处理、压
制鬼波处理和零相位化处理的子波，这样，使得该反褶积算子处理地震数据时，能够有效消除地震数据中的气泡信息、压制鬼波数据、以及使得地震数据归于零相位。
[0152]
步骤108：计算机设备通过反褶积算子，处理第一地震数据，得到地震道集的反褶积地震数据。
[0153]
该步骤可以通过以下步骤(1)-(2)实现：
[0154]
(1)计算机设备通过反褶积算子，处理第一地震数据，得到第五地震数据。
[0155]
需要说明的是，由于反褶积算子为与方位角和出射角有关的算子，这样，基于反褶积算子处理第一地震数据，能够有效补偿地震道的地震子波在方位角和出射角上的差异，使得到的第五地震数据的整体的最终地震子波趋于一致，从而提高了第五地震数据的准确性，实现了信号反褶积方法从一维到三维的跨越，使得三维信号反褶积方法处理宽频领域的地震数据更加理想。
[0156]
(2)计算机设备对第五地震数据进行τ-p域反变换处理，得到反褶积地震数据。
[0157]
由于第五地震数据为τ-p域内的地震数据，因此，需要对第五地震数据进行τ-p域反变换处理，以得到原有的时空域的地震数据，即反褶积地震数据。
[0158]
需要说明的是，该地震道集所属的目标区域包括多个地震道集，将多个地震道集按照上述所述的反褶积方法进行处理，得到每个地震道集对应的反褶积地震数据。
[0159]
参见图13，图13从上至下依次为信号反褶积之前的地震数据的地震道集、一维信号反褶积方法处理后的地震数据的地震道集、本技术实施例提供的信号反褶积方法处理后的地震数据的地震道集。参见图14，图14从左至右依次为信号反褶积之前的地震数据的地震道集、一维信号反褶积方法处理后的地震数据的地震道集、本技术实施例提供的信号反褶积方法处理后的地震数据的地震道集的局部放大图。从图13和图14中可以看出，通过本技术实施例提供的信号反褶积方法处理的地震数据的地震道集的气泡得到了明显压制，压制鬼波的效果更佳理想。
[0160]
参见图15，图15为从上至下依次为一维信号反褶积方法处理后的地震数据对应的地震子波的示意图、本技术实施例提供的信号反褶积方法处理后的地震数据对应的地震子波的示意图。图16为一维信号反褶积方法处理后的地震数据对应的地震子波的和本技术实施例提供的信号反褶积方法处理后的地震数据对应的地震子波的叠加图。参见图17，图17中包括信号反褶积之前的地震数据的频谱、一维信号反褶积方法处理后的地震数据的频谱、本技术实施例提供的信号反褶积方法处理后的地震数据的频谱。从图15、16和17可以看出，通过本技术实施例提供的方法得到的地震子波的旁瓣更小，峰值更尖锐，频率信息更加丰富，低频信息更强。说明通过本技术实施例提供的方法处理地震数据，得到的地震子波效果更加理想。
[0161]
参见如18，图18从左至右依次为一维信号反褶积方法处理后的地震数据的地层剖面图、本技术实施例提供的信号反褶积方法处理后的地震数据的地层剖面图。对比可以看出，右图的地层剖面中有效消除了气泡和压制了鬼波，且每一相轴更精细，从而提高了地震数据对应的地层剖面图的分辨率。
[0162]
本技术实施例提供了一种信号反褶积方法，由于该方法基于地震道的方位角和出射角对地震道进行了划分，这样，再基于划分方位角和出射角的地震道来依次确定单元子波数据、远场子波数据，使得确定的远场子波数据充分考虑了方位角和出射角的因素；进而
使得基于远场子波数据确定的反褶积算子与方位角和出射角相关，这样，基于反褶积算子处理第一地震数据时，能够有效补偿方位角和出射角带来的差异，从而使得基于反褶积算子确定的反褶积地震数据更加准确。
[0163]
本技术实施例还提供了一种信号反褶积装置，参见图19，装置包括：
[0164]
第一确定模块1901，用于确定地震道集初始的第一地震数据，地震道集包括多个地震道；
[0165]
第二确定模块1902，用于基于第一地震数据，确定地震道集的第二地震数据，以及从第一地震数据中确定多个地震道的方位角和出射角，第二地震数据包括每个地震道对应的射线参数；
[0166]
第三确定模块1903，用于对于每个地震道，基于地震道的方位角和出射角，确定地震道对应的角度组合，角度组合包括方位角和出射角；
[0167]
第四确定模块1904，用于对于多个预设角度组合中的每个预设角度组合，基于多个地震道对应的角度组合，确定目标地震道，目标地震道的角度组合属于预设角度组合；
[0168]
叠加模块1905，用于对于每个预设角度组合，叠加预设角度组合对应的多个目标地震道的射线参数，得到预设角度组合对应的单元子波数据；
[0169]
第五确定模块1906，用于基于多个预设角度组合对应的单元子波数据，确定地震道集的远场子波数据；
[0170]
生成模块1907，用于基于远场子波数据，生成地震道集的反褶积算子；
[0171]
处理模块1908，用于通过反褶积算子，处理第一地震数据，得到地震道集的反褶积地震数据。
[0172]
在一种可能的实现方式中，生成模块1907，用于：
[0173]
对远场子波数据进行气泡消除处理、压制鬼波处理和零相位化处理，得到期望子波数据；
[0174]
确定将远场子波数据转换为期望子波数据所使用的反褶积算子。
[0175]
在一种可能的实现方式中，第一确定模块1901，包括：
[0176]
获取单元，用于获取地震道集的原始地震数据；
[0177]
第一处理单元，用于对原始地震数据进行预处理，得到第三地震数据；
[0178]
第二处理单元，用于对第三地震数据进行τ-p域正变换处理，得到地震道集的第一地震数据。
[0179]
在一种可能的实现方式中，第一处理单元，包括：
[0180]
处理子单元，用于对原始地震数据进行规则化处理和插值处理，得到第四地震数据；
[0181]
第一确定子单元，用于从第四地震数据中确定目标地震数据，目标地震数据为波形发生预设波动的数据；
[0182]
第二确定子单元，用于确定目标地震数据中的下行波数据；
[0183]
压制单元，用于对下行波数据中的多次波数据进行压制，得到第三地震数据。
[0184]
在一种可能的实现方式中，第二确定子单元，用于：
[0185]
对目标地震数据进行波场分离处理，得到目标地震数据中的下行波数据。
[0186]
在一种可能的实现方式中，第二确定模块1902，用于：
[0187]
对第一地震数据进行噪音压制处理，得到压制噪音后的第一地震数据；
[0188]
对压制噪音后的第一地震数据进行同相轴拉平处理，得到同相轴拉平后的第一地震数据；
[0189]
对同相轴拉平后的第一地震数据进行畸变切除处理，得到切除畸变后的第一地震数据；
[0190]
对切除畸变后的第一地震数据进行边界外推处理，得到地震道集的第二地震数据。
[0191]
在一种可能的实现方式中，处理模块1908，用于：
[0192]
通过反褶积算子，处理第一地震数据，得到第五地震数据；
[0193]
对第五地震数据进行τ-p域反变换处理，得到反褶积地震数据。
[0194]
图20示出了本技术一个示例性实施例提供的计算机设备2000的结构框图。该计算机设备2000可以是便携式移动计算机设备，比如：智能手机、平板电脑、mp3播放器(moving picture experts group audio layer iii，动态影像专家压缩标准音频层面3)、mp4(moving picture experts group audio layer iv，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。计算机设备2000还可能被称为用户设备、便携式计算机设备、膝上型计算机设备、台式计算机设备等其他名称。
[0195]
通常，计算机设备2000包括有：处理器2001和存储器2002。
[0196]
处理器2001可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器2001可以采用dsp(digital signal processing，数字信号处理)、fpga(field－programmable gate array，现场可编程门阵列)、pla(programmable logic array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器2001也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称cpu(central processing unit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器2001可以集成有gpu(graphics processing unit，图像处理器)，gpu用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器2001还可以包括ai(artificial intelligence，人工智能)处理器，该ai处理器用于处理有关机器学习的计算操作。
[0197]
存储器2002可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器2002还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器2002中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器2001所执行以实现本技术中方法实施例提供的信号反褶积方法。
[0198]
在一些实施例中，计算机设备2000还可选包括有：外围设备接口2003和至少一个外围设备。处理器2001、存储器2002和外围设备接口2003之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口2003相连。具体地，外围设备包括：射频电路2004、显示屏2005、摄像头组件2006、音频电路2007、定位组件2008和电源2009中的至少一种。
[0199]
外围设备接口2003可被用于将i/o(input/output，输入/输出)相关的至少一个外围设备连接到处理器2001和存储器2002。在一些实施例中，处理器2001、存储器2002和外围
设备接口2003被集成在同一芯片或电路板上；在一些其他实施例中，处理器2001、存储器2002和外围设备接口2003中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
[0200]
射频电路2004用于接收和发射rf(radio frequency，射频)信号，也称电磁信号。射频电路2004通过电磁信号与通信网络以及其他通信设备进行通信。射频电路2004将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路2004包括：天线系统、rf收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路2004可以通过至少一种无线通信协议来与其它计算机设备进行通信。该无线通信协议包括但不限于：万维网、城域网、内联网、各代移动通信网络(2g、3g、4g及5g)、无线局域网和/或wifi(wireless fidelity，无线保真)网络。在一些实施例中，射频电路2004还可以包括nfc(near field communication，近距离无线通信)有关的电路，本技术对此不加以限定。
[0201]
显示屏2005用于显示ui(user interface，用户界面)。该ui可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏2005是触摸显示屏时，显示屏2005还具有采集在显示屏2005的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器2001进行处理。此时，显示屏2005还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏2005可以为一个，设置在计算机设备2000的前面板；在另一些实施例中，显示屏2005可以为至少两个，分别设置在计算机设备2000的不同表面或呈折叠设计；在另一些实施例中，显示屏2005可以是柔性显示屏，设置在计算机设备2000的弯曲表面上或折叠面上。甚至，显示屏2005还可以设置成非矩形的不规则图形，也即异形屏。显示屏2005可以采用lcd(liquid crystal display，液晶显示屏)、oled(organic light-emitting diode,有机发光二极管)等材质制备。
[0202]
摄像头组件2006用于采集图像或视频。可选地，摄像头组件2006包括前置摄像头和后置摄像头。通常，前置摄像头设置在计算机设备的前面板，后置摄像头设置在计算机设备的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及vr(virtual reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件2006还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。
[0203]
音频电路2007可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器2001进行处理，或者输入至射频电路2004以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在计算机设备2000的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器2001或射频电路2004的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路2007还可以包括耳机插孔。
[0204]
定位组件2008用于定位计算机设备2000的当前地理位置，以实现导航或lbs
(location based service，基于位置的服务)。定位组件2008可以是基于美国的gps(global positioning system，全球定位系统)、中国的北斗系统或俄罗斯的伽利略系统的定位组件。
[0205]
电源2009用于为计算机设备2000中的各个组件进行供电。电源2009可以是交流电、直流电、一次性电池或可充电电池。当电源2009包括可充电电池时，该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池，无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。
[0206]
在一些实施例中，计算机设备2000还包括有一个或多个传感器2010。该一个或多个传感器2010包括但不限于：加速度传感器2011、陀螺仪传感器2012、压力传感器2013、指纹传感器2014、光学传感器2015以及接近传感器2016。
[0207]
加速度传感器2011可以检测以计算机设备2000建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器2011可以用于检测重力加速度在三个坐标轴上的分量。处理器2001可以根据加速度传感器2011采集的重力加速度信号，控制显示屏2005以横向视图或纵向视图进行用户界面的显示。加速度传感器2011还可以用于游戏或者用户的运动数据的采集。
[0208]
陀螺仪传感器2012可以检测计算机设备2000的机体方向及转动角度，陀螺仪传感器2012可以与加速度传感器2011协同采集用户对计算机设备2000的3d动作。处理器2001根据陀螺仪传感器2012采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变ui)、拍摄时的图像稳定、游戏控制以及惯性导航。
[0209]
压力传感器2013可以设置在计算机设备2000的侧边框和/或显示屏2005的下层。当压力传感器2013设置在计算机设备2000的侧边框时，可以检测用户对计算机设备2000的握持信号，由处理器2001根据压力传感器2013采集的握持信号进行左右手识别或快捷操作。当压力传感器2013设置在显示屏2005的下层时，由处理器2001根据用户对显示屏2005的压力操作，实现对ui界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
[0210]
指纹传感器2014用于采集用户的指纹，由处理器2001根据指纹传感器2014采集到的指纹识别用户的身份，或者，由指纹传感器2014根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器2001授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器2014可以被设置在计算机设备2000的正面、背面或侧面。当计算机设备2000上设置有物理按键或厂商logo时，指纹传感器2014可以与物理按键或厂商logo集成在一起。
[0211]
光学传感器2015用于采集环境光强度。在一个实施例中，处理器2001可以根据光学传感器2015采集的环境光强度，控制显示屏2005的显示亮度。具体地，当环境光强度较高时，调高显示屏2005的显示亮度；当环境光强度较低时，调低显示屏2005的显示亮度。在另一个实施例中，处理器2001还可以根据光学传感器2015采集的环境光强度，动态调整摄像头组件2006的拍摄参数。
[0212]
接近传感器2016，也称距离传感器，通常设置在计算机设备2000的前面板。接近传感器2016用于采集用户与计算机设备2000的正面之间的距离。在一个实施例中，当接近传感器2016检测到用户与计算机设备2000的正面之间的距离逐渐变小时，由处理器2001控制
显示屏2005从亮屏状态切换为息屏状态；当接近传感器2016检测到用户与计算机设备2000的正面之间的距离逐渐变大时，由处理器2001控制显示屏2005从息屏状态切换为亮屏状态。
[0213]
本领域技术人员可以理解，图20中示出的结构并不构成对计算机设备2000的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
[0214]
本技术实施例还提供了一种计算机可读存储介质，计算机可读存储介质中存储至少一条指令，至少一条指令由处理器加载并执行，以实现上述任一实现方式的信号反褶积方法所执行的操作。
[0215]
本技术实施例还提供了一种计算机程序产品或计算机程序，计算机程序产品或计算机程序包括计算机程序代码，计算机程序代码存储在计算机可读存储介质中。计算机设备的处理器从计算机可读存储介质读取计算机程序代码，处理器执行计算机程序代码，使得计算机设备执行上述的信号反褶积方法所执行的操作。
[0216]
在一些实施例中，本技术实施例所涉及的计算机程序可被部署在一个计算机设备上执行，或者在位于一个地点的多个计算机设备上执行，又或者，在分布在多个地点且通过通信网络互连的多个计算机设备上执行，分布在多个地点且通过通信网络互连的多个计算机设备可以组成区块链系统。
[0217]
本技术实施例提供了一种信号反褶积方法，由于该方法基于地震道的方位角和出射角对地震道进行了划分，这样，再基于划分方位角和出射角的地震道来依次确定单元子波数据、远场子波数据，使得确定的远场子波数据充分考虑了方位角和出射角的因素；进而使得基于远场子波数据确定的反褶积算子与方位角和出射角相关，这样，基于反褶积算子处理第一地震数据时，能够有效补偿方位角和出射角带来的差异，从而使得基于反褶积算子确定的反褶积地震数据更加准确。
[0218]
以上仅为本技术的可选实施例，并不用以限制本技术，凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。