异常体叠后成像构造方法、装置、电子设备及介质与流程

1.本发明涉及地震理论与技术领域，更具体地，涉及一种异常体叠后成像构造方法、装置、电子设备及介质。


背景技术：

2.塔里木盆地碳酸盐岩油藏埋藏深、非均质性强、裂缝及孔洞发育。这些特点给该类储层的精细预测带来很大困难，西部碳酸盐岩油藏主要分布于裂缝孔洞型(简称串珠)的储层中。对该类储层多年的研究发现，碳酸盐岩地层的塌陷构造会形成大规模的长条状分布，在长条状的塌陷带会发育有圆形和扁形的断层。当溶洞碎屑与碳酸盐岩储层围岩具有较大的波阻抗差时，这些塌陷的溶洞及其碎屑在地震剖面上会表现为串珠状反射，因而碳酸盐岩储层的串珠状反射通常也被解释为塌陷的古溶洞。另外，由于埋藏深度增加导致差异压实作用会在塌陷的溶洞上方产生裂缝，从而使得这些溶洞连片成为规模较大的古溶洞体系，进而形成良好的油藏储集空间，通常称这类溶洞体系为缝洞连通体。钻井证实了这类塌陷溶洞体系正是好的油气产层。
3.对于串珠状反射特征的表征，目前已经有多种技术方案，在实际工作中常采用叠后振幅剖面透明化显示的方式来表征或者通过振幅属性提取技术实现。串珠溶洞在地下属于异常体，对于这类异常体的成像，目前刚刚有一些保幅的叠前叠后的叠加成像、高分辨、分方位、提频等方法技术，其中叠后的突出表征技术虽然效率高，但效果不佳。
4.串珠类异常体在塔里木盆地，作为重要的油气储层，在地震成像中经常存在成像不突出，成像不准，串珠构造不清晰，被反射层掩盖，弱串珠若有若无而无法确认的现象，从而对油气钻探的井位目标造成一定的不确定性风险。
5.因此，有必要开发一种异常体叠后成像构造方法、装置、电子设备及介质。
6.公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。


技术实现要素：

7.本发明提出了一种异常体叠后成像构造方法、装置、电子设备及介质，其能够提取异常体反射能量，突出串珠类地质构造的地震成像特征，改善串珠类储层成像信噪比。
8.第一方面，本公开实施例提供了一种异常体叠后成像构造方法，包括：
9.根据初始成像体中串珠异常体的构造特征，将所述构造特征外的反射构造进行去除与压制，获得地震成像体；
10.提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体；
11.根据所述属性体提取串珠异常体；
12.将所述初始成像体或所述地震成像体与串珠异常体进行叠加，获得成像数据体。
13.优选地，提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体包
括：
14.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得所述属性体。
15.优选地，根据所述属性体提取串珠异常体包括：
16.针对所述属性体设置门槛值，剔除所述属性体中门槛值以外的属性，获得串珠异常体。
17.优选地，针对所述属性体设置门槛值，剔除所述属性体中门槛值以外的属性，获得串珠异常体包括：
18.将所述属性体中小于所述门槛值的属性值置为0，将其余属性值置为1；
19.提取属性值为1的属性体，即为所述串珠异常体。
20.优选地，所述门槛值的大小与串珠的表征尺度呈正相关。
21.优选地，通过公式(1)将所述初始成像体与串珠异常体进行加权叠加，获得成像数据体：
22.c＝m
×
b2
×
a+a
ꢀꢀꢀ
(1)
23.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a为初始成像体。
24.优选地，通过公式(2)将所述地震成像体与串珠异常体进行加权叠加，获得成像数据体：
25.c＝m
×
b2
×
a’+a
’ꢀꢀꢀ
(2)
26.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a’为地震成像体。
27.作为本公开实施例的一种具体实现方式，
28.第二方面，本公开实施例还提供了一种异常体叠后成像构造装置，包括：
29.去噪模块，根据初始成像体中串珠异常体的构造特征，将所述构造特征外的反射构造进行去除与压制，获得地震成像体；
30.属性体确定模块，提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体；
31.串珠异常体确定模块，根据所述属性体提取串珠异常体；
32.叠加模块，将所述初始成像体或所述地震成像体与串珠异常体进行叠加，获得成像数据体。
33.优选地，提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体包括：
34.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得所述属性体。
35.优选地，根据所述属性体提取串珠异常体包括：
36.针对所述属性体设置门槛值，剔除所述属性体中门槛值以外的属性，获得串珠异常体。
37.优选地，针对所述属性体设置门槛值，剔除所述属性体中门槛值以外的属性，获得串珠异常体包括：
38.将所述属性体中小于所述门槛值的属性值置为0，将其余属性值置为1；
39.提取属性值为1的属性体，即为所述串珠异常体。
40.优选地，所述门槛值的大小与串珠的表征尺度呈正相关。
41.优选地，通过公式(1)将所述初始成像体与串珠异常体进行加权叠加，获得成像数据体：
42.c＝m
×
b2
×
a+a
ꢀꢀꢀ
(1)
43.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a为初始成像体。
44.优选地，通过公式(2)将所述地震成像体与串珠异常体进行加权叠加，获得成像数据体：
45.c＝m
×
b2
×
a’+a
’ꢀꢀꢀ
(2)
46.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a’为地震成像体。
47.第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
48.存储器，存储有可执行指令；
49.处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的异常体叠后成像构造方法。
50.第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的异常体叠后成像构造方法。
51.其有益效果在于：
52.本发明能够加强串珠类异常体的成像反射，凸显串珠的构造特征，从而对串珠位置确认、弱串珠是否存在、凸显串珠目标等方面，作出成像改进。
53.本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
54.通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
55.图1示出了根据本发明的一个实施例的异常体叠后成像构造方法的步骤的流程图。
56.图2a和图2b分别示出了根据本发明的一个实施例的初始成像体与地震成像体的对比示意图。
57.图3a和图3b分别示出了根据本发明的一个实施例的地震成像体与属性体的对比示意图。
58.图4a、图4b、图4c分别示出了根据本发明的一个实施例的属性体、串珠异常体与成像数据体的对比示意图。
59.图5示出了根据本发明的一个实施例的一种异常体叠后成像构造装置的框图。
60.附图标记说明：
61.201、去噪模块；202、属性体确定模块；203、串珠异常体确定模块；204、叠加模块。
具体实施方式
62.下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
63.本发明提供一种异常体叠后成像构造方法，包括：
64.根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，获得地震成像体；
65.提取地震成像体中串珠异常体的构造特征属性，获得属性体；
66.根据属性体提取串珠异常体；
67.将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。
68.在一个示例中，提取地震成像体中串珠异常体的构造特征属性，获得属性体包括：
69.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得属性体。
70.在一个示例中，根据属性体提取串珠异常体包括：
71.针对属性体设置门槛值，剔除属性体中门槛值以外的属性，获得串珠异常体。
72.在一个示例中，针对属性体设置门槛值，剔除属性体中门槛值以外的属性，获得串珠异常体包括：
73.将属性体中小于门槛值的属性值置为0，将其余属性值置为1；
74.提取属性值为1的属性体，即为串珠异常体。
75.在一个示例中，门槛值的大小与串珠的表征尺度呈正相关。
76.在一个示例中，通过公式(1)将初始成像体与串珠异常体进行加权叠加，获得成像数据体：
77.c＝m
×
b2
×
a+a
ꢀꢀꢀ
(1)
78.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a为初始成像体。
79.在一个示例中，通过公式(2)将地震成像体与串珠异常体进行加权叠加，获得成像数据体：
80.c＝m
×
b2
×
a’+a
’ꢀꢀꢀ
(2)
81.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a’为地震成像体。
82.具体地，根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，串珠类异常体具有球形构造特征，那么可以对具有线性特征的层状反射构造进行信号压制，将此类信号作为噪声，进行去除，获得地震成像体。
83.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得属性体；对于串珠类异常体，可以采用诸如信号包络类属性技术，该类技术能够对不同构造特征的成像体计算出不同的属性值；信号包络属性计算是可以选用的一种计算方法，无论采用什么计算方法，只要能将该类异常体的属性特征凸显出来，其属性值的数值大小与异常体类以外信号属性值大小差异较大即可。
84.针对属性体设置门槛值，将属性体中小于门槛值的属性值置为0，将其余属性值置为1，剔除属性体中门槛值以外的属性，提取属性值为1的属性体，即为串珠异常体；只要能够将门槛值内外的数值，通过新的数值设置，明显的区别开就可以。
85.其中，门槛值的大小与串珠的表征尺度呈正相关：需要表征尺度较大的串珠，门槛
值较大，需要兼顾较小的中弱串珠，则门槛值较小。
86.将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。
87.本发明还提供一种异常体叠后成像构造装置，包括：
88.去噪模块，根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，获得地震成像体；
89.属性体确定模块，提取地震成像体中串珠异常体的构造特征属性，获得属性体；
90.串珠异常体确定模块，根据属性体提取串珠异常体；
91.叠加模块，将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。
92.在一个示例中，提取地震成像体中串珠异常体的构造特征属性，获得属性体包括：
93.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得属性体。
94.在一个示例中，根据属性体提取串珠异常体包括：
95.针对属性体设置门槛值，剔除属性体中门槛值以外的属性，获得串珠异常体。
96.在一个示例中，针对属性体设置门槛值，剔除属性体中门槛值以外的属性，获得串珠异常体包括：
97.将属性体中小于门槛值的属性值置为0，将其余属性值置为1；
98.提取属性值为1的属性体，即为串珠异常体。
99.在一个示例中，门槛值的大小与串珠的表征尺度呈正相关。
100.在一个示例中，通过公式(1)将初始成像体与串珠异常体进行加权叠加，获得成像数据体：
101.c＝m
×
b2
×
a+a
ꢀꢀꢀ
(1)
102.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a为初始成像体。
103.在一个示例中，通过公式(2)将地震成像体与串珠异常体进行加权叠加，获得成像数据体：
104.c＝m
×
b2
×
a’+a
’ꢀꢀꢀ
(2)
105.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a’为地震成像体。
106.具体地，根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，串珠类异常体具有球形构造特征，那么可以对具有线性特征的层状反射构造进行信号压制，将此类信号作为噪声，进行去除，获得地震成像体。
107.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得属性体；对于串珠类异常体，可以采用诸如信号包络类属性技术，该类技术能够对不同构造特征的成像体计算出不同的属性值；信号包络属性计算是可以选用的一种计算方法，无论采用什么计算方法，只要能将该类异常体的属性特征凸显出来，其属性值的数值大小与异常体类以外信号属性值大小差异较大即可。
108.针对属性体设置门槛值，将属性体中小于门槛值的属性值置为0，将其余属性值置为1，剔除属性体中门槛值以外的属性，提取属性值为1的属性体，即为串珠异常体；只要能够将门槛值内外的数值，通过新的数值设置，明显的区别开就可以。
109.其中，门槛值的大小与串珠的表征尺度呈正相关：需要表征尺度较大的串珠，门槛值较大，需要兼顾较小的中弱串珠，则门槛值较小。
110.将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。
111.本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的异常体叠后成像构造方法。
112.本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的异常体叠后成像构造方法。
113.为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
114.实施例1
115.图1示出了根据本发明的一个实施例的异常体叠后成像构造方法的步骤的流程图。
116.如图1所示，该异常体叠后成像构造方法包括：步骤101，根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，获得地震成像体；步骤102，提取地震成像体中串珠异常体的构造特征属性，获得属性体；步骤103，根据属性体提取串珠异常体；步骤104，将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。
117.以某实际油气勘探区块的一条二维数据线为例：
118.图2a和图2b分别示出了根据本发明的一个实施例的初始成像体与地震成像体的对比示意图。
119.根据图2a所示的初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，获得地震成像体，如图2b所示。
120.图3a和图3b分别示出了根据本发明的一个实施例的地震成像体与属性体的对比示意图。
121.针对地震成像体，提取其中异常体的构造特征属性，获得属性体，如图3b所示。
122.图4a、图4b、图4c分别示出了根据本发明的一个实施例的属性体、串珠异常体与成像数据体的对比示意图。
123.对异常体的构造特征属性体设置门槛值，将属性体中门槛值以外的属性消除。从而提取出串珠类异常体。这里对属性体设置门槛值为420，大于等于420以上的数据点位置，数据值大小设为1.9；属性体中小于420的数据点位置，数据值大小设为0，从而得到串珠异常体，如图4b所示。
124.将串珠异常体与初始成像体进行加权叠加，因为有弱串珠要进行突出成像表征，则设m＝3，从而得到成像数据体，如图4c所示，4.6s附近串珠得到成像加强凸显。
125.针对串珠类地质异常体的地震成像，本发明通过在叠后成像体中提取其异常能量特征，再对这类特征进行针对性能量加强，进而在地震成像体中能够凸显此类异常油气储层的构造，从而为此类井位目标落实，提供信噪比较高的构造成像数据体。
126.实施例2
127.图5示出了根据本发明的一个实施例的一种异常体叠后成像构造装置的框图。
128.如图5所示，该异常体叠后成像构造装置，包括：
129.去噪模块201，根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，获得地震成像体；
130.属性体确定模块202，提取地震成像体中串珠异常体的构造特征属性，获得属性体；
131.串珠异常体确定模块203，根据属性体提取串珠异常体；
132.叠加模块204，将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。
133.作为可选方案，提取地震成像体中串珠异常体的构造特征属性，获得属性体包括：
134.针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得属性体。
135.作为可选方案，根据属性体提取串珠异常体包括：
136.针对属性体设置门槛值，剔除属性体中门槛值以外的属性，获得串珠异常体。
137.作为可选方案，针对属性体设置门槛值，剔除属性体中门槛值以外的属性，获得串珠异常体包括：
138.将属性体中小于门槛值的属性值置为0，将其余属性值置为1；
139.提取属性值为1的属性体，即为串珠异常体。
140.作为可选方案，门槛值的大小与串珠的表征尺度呈正相关。
141.作为可选方案，通过公式(1)将初始成像体与串珠异常体进行加权叠加，获得成像数据体：
142.c＝m
×
b2
×
a+a
ꢀꢀꢀ
(1)
143.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a为初始成像体。
144.作为可选方案，通过公式(2)将地震成像体与串珠异常体进行加权叠加，获得成像数据体：
145.c＝m
×
b2
×
a’+a
’ꢀꢀꢀ
(2)
146.其中，c为成像数据体，m为加权系数，b2为串珠异常体，a’为地震成像体。
147.实施例3
148.本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述异常体叠后成像构造方法。
149.根据本公开实施例的电子设备包括存储器和处理器。
150.该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
151.该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
152.本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
153.有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
154.实施例4
155.本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的异常体叠后成像构造方法。
156.根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
157.上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
158.本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
159.以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

技术特征：
1.一种异常体叠后成像构造方法，其特征在于，包括：根据初始成像体中串珠异常体的构造特征，将所述构造特征外的反射构造进行去除与压制，获得地震成像体；提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体；根据所述属性体提取串珠异常体；将所述初始成像体或所述地震成像体与串珠异常体进行叠加，获得成像数据体。2.根据权利要求1所述的异常体叠后成像构造方法，其中，提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体包括：针对不同构造特征的成像体计算得到不同的属性值，进而确定串珠异常体对应的属性值，获得所述属性体。3.根据权利要求1所述的异常体叠后成像构造方法，其中，根据所述属性体提取串珠异常体包括：针对所述属性体设置门槛值，剔除所述属性体中门槛值以外的属性，获得串珠异常体。4.根据权利要求1所述的异常体叠后成像构造方法，其中，针对所述属性体设置门槛值，剔除所述属性体中门槛值以外的属性，获得串珠异常体包括：将所述属性体中小于所述门槛值的属性值置为0，将其余属性值置为1；提取属性值为1的属性体，即为所述串珠异常体。5.根据权利要求4所述的异常体叠后成像构造方法，其中，所述门槛值的大小与串珠的表征尺度呈正相关。6.根据权利要求1所述的异常体叠后成像构造方法，其中，通过公式(1)将所述初始成像体与串珠异常体进行加权叠加，获得成像数据体：c＝m
×
b2
×
a+a
ꢀꢀ
(1)其中，c为成像数据体，m为加权系数，b2为串珠异常体，a为初始成像体。7.根据权利要求1所述的异常体叠后成像构造方法，其中，通过公式(2)将所述地震成像体与串珠异常体进行加权叠加，获得成像数据体：c＝m
×
b2
×
a’+a
’ꢀꢀꢀ
(2)其中，c为成像数据体，m为加权系数，b2为串珠异常体，a’为地震成像体。8.一种异常体叠后成像构造装置，其特征在于，包括：去噪模块，根据初始成像体中串珠异常体的构造特征，将所述构造特征外的反射构造进行去除与压制，获得地震成像体；属性体确定模块，提取所述地震成像体中所述串珠异常体的构造特征属性，获得属性体；串珠异常体确定模块，根据所述属性体提取串珠异常体；叠加模块，将所述初始成像体或所述地震成像体与串珠异常体进行叠加，获得成像数据体。9.一种电子设备，其特征在于，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一项所述的异常体叠后成像构造方法。
10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的异常体叠后成像构造方法。

技术总结
本申请公开了一种异常体叠后成像构造方法、装置、电子设备及介质。该方法可以包括：根据初始成像体中串珠异常体的构造特征，将构造特征外的反射构造进行去除与压制，获得地震成像体；提取地震成像体中串珠异常体的构造特征属性，获得属性体；根据属性体提取串珠异常体；将初始成像体或地震成像体与串珠异常体进行叠加，获得成像数据体。本发明提取异常体反射能量，突出串珠类地质构造的地震成像特征，改善串珠类储层成像信噪比。善串珠类储层成像信噪比。善串珠类储层成像信噪比。


技术研发人员：姜大建 刘志远 王震 谭未一 莫延钢 肖彦君
受保护的技术使用者：中国石油化工股份有限公司石油勘探开发研究院
技术研发日：2022.03.02
技术公布日：2023/9/12