CN116774280B - 一种鬼波压制定量质控方法、装置、电子设备及存储介质 - Google Patents
一种鬼波压制定量质控方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN116774280B CN116774280B CN202310754376.0A CN202310754376A CN116774280B CN 116774280 B CN116774280 B CN 116774280B CN 202310754376 A CN202310754376 A CN 202310754376A CN 116774280 B CN116774280 B CN 116774280B
- Authority
- CN
- China
- Prior art keywords
- mean square
- root mean
- ghost
- amplitude
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003908 quality control method Methods 0.000 title claims abstract description 44
- 238000009472 formulation Methods 0.000 title claims description 6
- 239000000203 mixture Substances 0.000 title claims description 6
- 230000000694 effects Effects 0.000 claims abstract description 76
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000003631 expected effect Effects 0.000 claims abstract description 22
- 238000004590 computer program Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 8
- 230000001629 suppression Effects 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 8
- 238000007906 compression Methods 0.000 description 5
- 230000006835 compression Effects 0.000 description 5
- 238000013016 damping Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004031 devitrification Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种鬼波压制定量质控方法、装置、电子设备及存储介质。该方法包括：获取原始地震数据；对原始地震数据进行提取，得到第一海底反射波均方根振幅；确定鬼波压制参数，根据鬼波压制参数对原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；对目标地震数据进行提取，得到第二海底反射波均方根振幅；根据第一海底反射波均方根振幅和第二海底反射波均方根振幅，确定地震数据振幅衰减率；根据地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；在鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。本技术方案，实现对鬼波压制效果的定量质控，在鬼波压制参数选择上，降低人为因素干扰，提高参数选择精度。
Description
技术领域
本发明涉及油气勘探开发技术领域，尤其涉及一种鬼波压制定量质控方法、装置、电子设备及存储介质。
背景技术
近年来，随着海上油气勘探开发对地震资料的频带及保真度要求越来越高，海上地震资料的鬼波影响也逐渐受到重视。
由于海上地震采集时气枪、检波器一般放置于水面以下特定的深度，这一特殊的观测方式使得记录到的地震信号包含了源端及检端的鬼波波场，鬼波的存在造成地震资料有效频道窄、低频能量弱、波组特征不清楚的问题。现阶段针对鬼波的压制技术得到迅速发展，但是针对鬼波压制效果的质量控制，只能基于剖面的视觉效果凭借质控人员的经验来选择压制参数。
基于剖面的视觉效果凭借质控人员的经验来选择压制参数，较大的限制了去鬼波参数的优选效果。
发明内容
本发明提供了一种鬼波压制定量质控方法、装置、电子设备及存储介质，能够实现对鬼波压制效果的定量质控，在鬼波压制参数选择上，降低了人为因素干扰，提高了参数选择精度。
根据本发明的一方面，提供了一种鬼波压制定量质控方法，该方法包括：
获取原始地震数据；
对所述原始地震数据进行提取，得到第一海底反射波均方根振幅；
确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；
对所述目标地震数据进行提取，得到第二海底反射波均方根振幅；
根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率；
根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的；
在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。
根据本发明的另一方面，提供了一种鬼波压制定量质控装置，该装置包括：
原始地震数据获取模块，用于获取原始地震数据；
第一海底反射波均方根振幅得到模块，用于对所述原始地震数据进行提取，得到第一海底反射波均方根振幅；
目标地震数据得到模块，用于确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；
第二海底反射波均方根振幅得到模块，用于对所述目标地震数据进行提取，得到第二海底反射波均方根振幅；
地震数据振幅衰减率确定模块，用于根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率；
鬼波压制效果确定模块，用于根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的；
目标地震数据输出模块，用于在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。
根据本发明的另一方面，提供了一种电子设备，所述电子设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行本发明任一实施例所述的一种鬼波压制定量质控方法。
根据本发明的另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现本发明任一实施例所述的一种鬼波压制定量质控方法。
本发明实施例的技术方案，通过获取原始地震数据；对原始地震数据进行提取，得到第一海底反射波均方根振幅；确定鬼波压制参数，根据鬼波压制参数对原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；对目标地震数据进行提取，得到第二海底反射波均方根振幅；根据第一海底反射波均方根振幅和第二海底反射波均方根振幅，确定地震数据振幅衰减率；根据地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；在鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。本技术方案，能够实现对鬼波压制效果的定量质控，在鬼波压制参数选择上，降低了人为因素干扰，提高了参数选择精度，以及提高了鬼波压制效果。
应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是根据本发明实施例一提供的一种鬼波压制定量质控方法的流程图；
图2是本申请实施例一提供的鬼波压制定量质控的流程图；
图3是本申请实施例一提供的均方根振幅衰减率误差计算示意图；
图4为本发明实施例二提供的一种鬼波压制定量质控过程的示意图；
图5是本申请实施例二提供的均方根振幅统计结果；
图6是本申请实施例二提供的实际地震数据近道剖面对应的不同参数去缆鬼波效果；
图7为本发明实施例三提供的一种鬼波压制定量质控装置的结构示意图；
图8是实现本发明实施例的一种鬼波压制定量质控方法的电子设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例一
图1是根据本发明实施例一提供的一种鬼波压制定量质控方法的流程图，本实施例可适用于对鬼波压制效果进行定量质控的情况，该方法可以由一种鬼波压制定量质控装置来执行，该鬼波压制定量质控装置可以采用硬件和/或软件的形式实现，该鬼波压制定量质控装置可配置于电子设备中。如图1所示，该方法包括：
S110、获取原始地震数据。
在本方案中，原始地震数据可以是指去鬼波前的地震数据。可以基于地震勘探获取原始地震数据。
S120、对所述原始地震数据进行提取，得到第一海底反射波均方根振幅。
在本实施例中，可以基于原始地震数据，提取海底反射小时窗对应的均方根振幅，得到第一海底反射波均方根振幅A1。其中，小时窗的大小可以根据数据提取需求进行设置。
可选的，对所述原始地震数据进行提取，得到第一海底反射波均方根振幅，包括：
基于所述原始地震数据，统计预先确定的小时窗内的海底反射波均方根振幅，得到第一海底反射波均方根振幅。
在本方案中，在确定小时窗后，可以基于原始地震数据统计小时窗内的每个样点的振幅值，并计算出该小时窗内的均方根振幅。
通过确定第一海底反射波均方根振幅，能够基于第一海底反射波均方根振幅对鬼波压制效果进行定量质控。
S130、确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据。
在本实施例中，鬼波压制参数可以是指阻尼因子，不同参数对应了去鬼波程度的强弱。其中，阻尼因子的取值范围一般是0.5到10，采用不同的阻尼因子处理后得到不同的鬼波压制效果。
在本方案中，鬼波压制定量控制主要是对鬼波压制参数进行优选，从而能够使得鬼波压制效果达到最优。
S140、对所述目标地震数据进行提取，得到第二海底反射波均方根振幅。
在本方案中，可以使用鬼波压制后的目标地震数据，提取海底反射小时窗对应的均方根振幅，得到第二海底反射波均方根振幅A2。
可选的，对所述目标地震数据进行提取，得到第二海底反射波均方根振幅，包括：
基于所述目标地震数据，统计所述小时窗内的海底反射波均方根振幅，得到第二海底反射波均方根振幅。
其中，小时窗与原始地震数据进行提取时设置的小时窗大小相同。
在本方案中，在确定小时窗后，可以基于目标地震数据统计小时窗内的每个样点的振幅值，并计算出该小时窗内的均方根振幅。
通过确定第二海底反射波均方根振幅，能够基于第二海底反射波均方根振幅对鬼波压制效果进行定量质控。
S150、根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率。
在本方案中，可以使用第一海底反射波均方根振幅A1和第二海底反射波均方根振幅A2，计算地震数据去鬼波后的振幅衰减率ΔA。
可选的，根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率，包括：
将所述第二海底反射波均方根振幅与所述第一海底反射波均方根振幅进行相除，得到地震数据振幅衰减率。
具体的，可以采用如下公式计算地震数据振幅衰减率ΔA：
ΔA＝A2/A1。
通过确定地震数据振幅衰减率，能够基于地震数据振幅衰减率实现鬼波压制效果的定量质控，提高了鬼波压制效果。
S160、根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的。
其中，地震子波可以是模拟的远场子波，也可以是从原始地震数据中提取的子波。
在本方案中，可以根据地震子波的均方根振幅，以及去除鬼波能量后的地震子波的均方根振幅，计算地震子波去鬼波后的振幅衰减率，得到地震子波振幅衰减率。
具体的，可以将地震数据振幅衰减率和地震子波振幅衰减率进行比较，并根据地震数据振幅衰减率和地震子波振幅衰减率的误差量判断鬼波压制效果。
在本实施例中，图2是本申请实施例一提供的鬼波压制定量质控的流程图，如图2所示，模拟远场子波，即地震子波，统计去鬼波前后均方根振幅衰减率；以及，输入去鬼波前的地震数据，即原始地震数据，基于鬼波压制参数对原始地震数据进行鬼波压制处理，统计去鬼波前后均方根振幅衰减率；将均方根振幅衰减率进行比对，确定鬼波压制效果。
可选的，所述地震子波振幅衰减率的确定过程包括步骤B1-B4：
步骤B1、获取地震子波；
步骤B2、对所述地震子波的均方根振幅进行计算，得到第一地震子波均方根振幅；
具体的，使用地震子波W，计算含鬼波时窗对应的均方根振幅，得到第一地震子波均方根振幅Aw1。
步骤B3、对所述地震子波进行鬼波能量去除处理，得到鬼波能量去除后的地震子波，并对鬼波能量去除后的地震子波的均方根振幅进行计算，得到第二地震子波均方根振幅；
进一步的，使用地震子波W，采用褶积法去除相应的鬼波能量，统计和上述相同时窗对应的均方根振幅，得到第二地震子波均方根振幅Aw2。
步骤B4、根据所述第一地震子波均方根振幅和所述第二地震子波均方根振幅，确定地震子波振幅衰减率。
在本实施例中，可以使用第一地震子波均方根振幅Aw1和第二地震子波均方根振幅Aw2，计算地震子波去鬼波后的振幅衰减率ΔAw。
通过计算地震子波振幅衰减率，能够将地震子波振幅衰减率和地震数据振幅衰减率进行比较，确定鬼波压制效果，从而实现鬼波压制效果的定量质控。
可选的，根据所述第一地震子波均方根振幅和所述第二地震子波均方根振幅，确定地震子波振幅衰减率，包括：
将所述第二地震子波均方根振幅与所述第一地震子波均方根振幅进行相除，得到地震子波振幅衰减率。
具体的，可以采用如下公式计算地震子波振幅衰减率ΔAw：
ΔAw＝Aw2/Aw1。
通过计算地震子波振幅衰减率，能够将地震子波振幅衰减率和地震数据振幅衰减率进行比较，确定鬼波压制效果，从而实现鬼波压制效果的定量质控。
进一步的，图3是本申请实施例一提供的均方根振幅衰减率误差计算示意图，如图3所示，模拟远场子波，即地震子波，对地震子波进行鬼波压制，确定统计时窗，并分别统计鬼波压制前后地震子波的均方根振幅Aw1和Aw2，计算均方根振幅衰减率ΔAw＝Aw2/Aw1。以及，获取实际地震数据，即获取原始地震数据，对原始地震数据进行鬼波压制，确定统计时窗，并分别统计鬼波压制前后原始地震数据的均方根振幅A1和A2，计算均方根振幅衰减率ΔA＝A2/A1。从而得到均方根振幅衰减率误差量Δ＝|ΔAw-ΔA|，并基于均方根振幅衰减率误差量判断鬼波压制效果。
S170、在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。
在本方案中，预期效果可以根据鬼波压制需求进行设置。如图2所示，当鬼波压制效果满足预期效果时，则可以输出鬼波压制后的目标地震数据，以用于油气勘探开发指导。
本发明实施例的技术方案，通过获取原始地震数据；对原始地震数据进行提取，得到第一海底反射波均方根振幅；确定鬼波压制参数，根据鬼波压制参数对原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；对目标地震数据进行提取，得到第二海底反射波均方根振幅；根据第一海底反射波均方根振幅和第二海底反射波均方根振幅，确定地震数据振幅衰减率；根据地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；在鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。通过执行本技术方案，能够实现对鬼波压制效果的定量质控，在鬼波压制参数选择上，降低了人为因素干扰，提高了参数选择精度，以及提高了鬼波压制效果。
实施例二
图4为本发明实施例二提供的一种鬼波压制定量质控过程的示意图，本实施例与上述实施例之间的关系是对鬼波压制参数确定过程的详细描述。如图4所示，该方法包括：
S410、获取原始地震数据。
S420、对所述原始地震数据进行提取，得到第一海底反射波均方根振幅。
S430、确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据。
S440、对所述目标地震数据进行提取，得到第二海底反射波均方根振幅。
S450、根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率。
S460、根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的。
S470、在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。
S480、在所述鬼波压制效果不满足预先确定的预期效果的情况下，根据所述地震数据振幅衰减率和所述地震子波振幅衰减率的误差量修改鬼波压制参数，并基于修改后的鬼波压制参数对所述原始地震数据进行鬼波压制，直至鬼波压制效果满足预先确定的预期效果，输出目标地震数据。
在本方案中，如图2所示，在鬼波压制效果不满足预先确定的预期效果的情况下，重新修改鬼波压制参数。
具体的，可以根据地震数据振幅衰减率和地震子波振幅衰减率的误差量修改鬼波压制参数，并基于修改后的鬼波压制参数对原始地震数据进行鬼波压制，并判断修改后的鬼波压制参数对原始地震数据的鬼波压制效果，直至鬼波压制效果满足预先确定的预期效果，并输出鬼波压制后的目标地震数据。其中，可以根据地震数据振幅衰减率和地震子波振幅衰减率的误差量的大小选择不同的鬼波压制参数。
在本方案中，图5是本申请实施例二提供的均方根振幅统计结果。如图5所示，首先统计远场子波选定时窗(包含一次反射及鬼波)内的均方根振幅，其中含源缆鬼波的远场子波均方根振幅为13.1727，只含源鬼波的远场子波均方根振幅为9.39499，不含鬼波的远场子波对应的均方根振幅为7.29474，通过计算，去缆鬼波后的振幅衰减率为71.3％，去源缆鬼波的振幅衰减率为55.4％。
进一步的，以去缆鬼波处理为例，应用优选的稀疏Tau-P反演自适应鬼波压制技术采用不同的鬼波压制参数进行去缆鬼波处理，得到不同鬼波压制参数去鬼波后的数据。图6是本申请实施例二提供的实际地震数据近道剖面对应的不同参数去缆鬼波效果，图6中的(a)表示去鬼波前的去缆鬼波效果，图6中的(b)表示DF0.5的去缆鬼波效果，图6中的(c)表示DF1.0的去缆鬼波效果，图6中的(d)表示DF1.5的去缆鬼波效果，图6中的(e)表示DF2.0的去缆鬼波效果，图6中的(f)表示DF2.5的去缆鬼波效果，图6中的(g)表示DF3.0的去缆鬼波效果，图6中的(h)表示DF4.0的去缆鬼波效果，图6中的(i)表示DF6.0的去缆鬼波效果，其中，DF0.5、DF1.0、DF1.5、DF2.0、DF2.5、DF3.0、DF4、DF6分别表示鬼波压制参数。如图6所示，选择部分海底平坦位置的近道剖面，统计选定时窗内的均方根振幅以及振幅衰减率。选定时窗内的均方根振幅以及振幅衰减率如表1所示，从振幅衰减率来看，应用参数DF1.5去缆鬼波后振幅衰减率和远场子波去缆鬼波的振幅衰减率最接近，因此可以优选参数DF1.5作为最佳参数。
表1
本发明实施例的技术方案，通过获取原始地震数据；对原始地震数据进行提取，得到第一海底反射波均方根振幅；确定鬼波压制参数，根据鬼波压制参数对原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；对目标地震数据进行提取，得到第二海底反射波均方根振幅；根据第一海底反射波均方根振幅和第二海底反射波均方根振幅，确定地震数据振幅衰减率；根据地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；在鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。在鬼波压制效果不满足预先确定的预期效果的情况下，根据地震数据振幅衰减率和地震子波振幅衰减率的误差量修改鬼波压制参数，并基于修改后的鬼波压制参数对原始地震数据进行鬼波压制，直至鬼波压制效果满足预先确定的预期效果，输出目标地震数据。通过执行本技术方案，能够实现对鬼波压制效果的定量质控，在鬼波压制参数选择上，降低了人为因素干扰，提高了参数选择精度，以及提高了鬼波压制效果。
实施例三
图7为本发明实施例三提供的一种鬼波压制定量质控装置的结构示意图。如图7所示，该装置包括：
原始地震数据获取模块710，用于获取原始地震数据；
第一海底反射波均方根振幅得到模块720，用于对所述原始地震数据进行提取，得到第一海底反射波均方根振幅；
目标地震数据得到模块730，用于确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；
第二海底反射波均方根振幅得到模块740，用于对所述目标地震数据进行提取，得到第二海底反射波均方根振幅；
地震数据振幅衰减率确定模块750，用于根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率；
鬼波压制效果确定模块760，用于根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的；
目标地震数据输出模块770，用于在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据。
可选的，第一海底反射波均方根振幅得到模块720，具体用于：
基于所述原始地震数据，统计预先确定的小时窗内的海底反射波均方根振幅，得到第一海底反射波均方根振幅。
可选的，第二海底反射波均方根振幅得到模块740，具体用于：
基于所述目标地震数据，统计所述小时窗内的海底反射波均方根振幅，得到第二海底反射波均方根振幅。
可选的，地震数据振幅衰减率确定模块750，具体用于：
将所述第二海底反射波均方根振幅与所述第一海底反射波均方根振幅进行相除，得到地震数据振幅衰减率。
可选的，鬼波压制效果确定模块760，包括：
地震子波获取单元，用于获取地震子波；
第一地震子波均方根振幅得到单元，用于对所述地震子波的均方根振幅进行计算，得到第一地震子波均方根振幅；
第二地震子波均方根振幅得到单元，用于对所述地震子波进行鬼波能量去除处理，得到鬼波能量去除后的地震子波，并对鬼波能量去除后的地震子波的均方根振幅进行计算，得到第二地震子波均方根振幅；
地震子波振幅衰减率确定单元，用于根据所述第一地震子波均方根振幅和所述第二地震子波均方根振幅，确定地震子波振幅衰减率。
可选的，地震子波振幅衰减率确定单元，具体用于：
将所述第二地震子波均方根振幅与所述第一地震子波均方根振幅进行相除，得到地震子波振幅衰减率。
可选的，所述装置还包括：
鬼波压制效果判断模块，用于在所述鬼波压制效果不满足预先确定的预期效果的情况下，根据所述地震数据振幅衰减率和所述地震子波振幅衰减率的误差量修改鬼波压制参数，并基于修改后的鬼波压制参数对所述原始地震数据进行鬼波压制，直至鬼波压制效果满足预先确定的预期效果，输出目标地震数据。
本发明实施例所提供的一种鬼波压制定量质控装置可执行本发明任意实施例所提供的一种鬼波压制定量质控方法，具备执行方法相应的功能模块和有益效果。
实施例四
图8示出了可以用来实施本发明的实施例的电子设备10的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
如图8所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(ROM)12、随机访问存储器(RAM)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(ROM)12中的计算机程序或者从存储单元18加载到随机访问存储器(RAM)13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出(I/O)接口15也连接至总线14。
电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如一种鬼波压制定量质控方法。
在一些实施例中，一种鬼波压制定量质控方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的一种鬼波压制定量质控方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行一种鬼波压制定量质控方法。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，CRT(阴极射线管)或者LCD(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(LAN)、广域网(WAN)、区块链网络和互联网。
计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。
Claims (7)
1.一种鬼波压制定量质控方法，其特征在于，包括：
获取原始地震数据；
对所述原始地震数据进行提取，得到第一海底反射波均方根振幅；
确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；
对所述目标地震数据进行提取，得到第二海底反射波均方根振幅；
根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率；
根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的；
在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据；
其中，根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率，包括：
将所述第二海底反射波均方根振幅与所述第一海底反射波均方根振幅进行相除，得到地震数据振幅衰减率；
其中，所述地震子波振幅衰减率的确定过程包括：
获取地震子波；
对所述地震子波的均方根振幅进行计算，得到第一地震子波均方根振幅；
对所述地震子波进行鬼波能量去除处理，得到鬼波能量去除后的地震子波，并对鬼波能量去除后的地震子波的均方根振幅进行计算，得到第二地震子波均方根振幅；
根据所述第一地震子波均方根振幅和所述第二地震子波均方根振幅，确定地震子波振幅衰减率；
其中，根据所述第一地震子波均方根振幅和所述第二地震子波均方根振幅，确定地震子波振幅衰减率，包括：
将所述第二地震子波均方根振幅与所述第一地震子波均方根振幅进行相除，得到地震子波振幅衰减率。
2.根据权利要求1所述的方法，其特征在于，对所述原始地震数据进行提取，得到第一海底反射波均方根振幅，包括：
基于所述原始地震数据，统计预先确定的小时窗内的海底反射波均方根振幅，得到第一海底反射波均方根振幅。
3.根据权利要求2所述的方法，其特征在于，对所述目标地震数据进行提取，得到第二海底反射波均方根振幅，包括：
基于所述目标地震数据，统计所述小时窗内的海底反射波均方根振幅，得到第二海底反射波均方根振幅。
4.根据权利要求1所述的方法，其特征在于，在确定鬼波压制效果之后，所述方法还包括：
在所述鬼波压制效果不满足预先确定的预期效果的情况下，根据所述地震数据振幅衰减率和所述地震子波振幅衰减率的误差量修改鬼波压制参数， 并基于修改后的鬼波压制参数对所述原始地震数据进行鬼波压制，直至鬼波压制效果满足预先确定的预期效果，输出目标地震数据。
5.一种鬼波压制定量质控装置，其特征在于，包括：
原始地震数据获取模块，用于获取原始地震数据；
第一海底反射波均方根振幅得到模块，用于对所述原始地震数据进行提取，得到第一海底反射波均方根振幅；
目标地震数据得到模块，用于确定鬼波压制参数，根据所述鬼波压制参数对所述原始地震数据进行鬼波压制，得到鬼波压制后的目标地震数据；
第二海底反射波均方根振幅得到模块，用于对所述目标地震数据进行提取，得到第二海底反射波均方根振幅；
地震数据振幅衰减率确定模块，用于根据所述第一海底反射波均方根振幅和所述第二海底反射波均方根振幅，确定地震数据振幅衰减率；
鬼波压制效果确定模块，用于根据所述地震数据振幅衰减率和地震子波振幅衰减率，确定鬼波压制效果；其中，所述地震子波振幅衰减率是根据地震子波进行确定的；
目标地震数据输出模块，用于在所述鬼波压制效果满足预先确定的预期效果的情况下，输出目标地震数据；
其中，地震数据振幅衰减率确定模块，具体用于：
将所述第二海底反射波均方根振幅与所述第一海底反射波均方根振幅进行相除，得到地震数据振幅衰减率；
其中，鬼波压制效果确定模块，包括：
地震子波获取单元，用于获取地震子波；
第一地震子波均方根振幅得到单元，用于对所述地震子波的均方根振幅进行计算，得到第一地震子波均方根振幅；
第二地震子波均方根振幅得到单元，用于对所述地震子波进行鬼波能量去除处理，得到鬼波能量去除后的地震子波，并对鬼波能量去除后的地震子波的均方根振幅进行计算，得到第二地震子波均方根振幅；
地震子波振幅衰减率确定单元，用于根据所述第一地震子波均方根振幅和所述第二地震子波均方根振幅，确定地震子波振幅衰减率；
其中，地震子波振幅衰减率确定单元，具体用于：
将所述第二地震子波均方根振幅与所述第一地震子波均方根振幅进行相除，得到地震子波振幅衰减率。
6.一种电子设备，其特征在于，所述电子设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-4中任一项所述的一种鬼波压制定量质控方法。
7.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-4中任一项所述的一种鬼波压制定量质控方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310754376.0A CN116774280B (zh) | 2023-06-25 | 2023-06-25 | 一种鬼波压制定量质控方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310754376.0A CN116774280B (zh) | 2023-06-25 | 2023-06-25 | 一种鬼波压制定量质控方法、装置、电子设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116774280A CN116774280A (zh) | 2023-09-19 |
| CN116774280B true CN116774280B (zh) | 2024-05-24 |
Family
ID=87989131
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310754376.0A Active CN116774280B (zh) | 2023-06-25 | 2023-06-25 | 一种鬼波压制定量质控方法、装置、电子设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116774280B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103308944A (zh) * | 2013-05-20 | 2013-09-18 | 中国石油天然气股份有限公司 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
| CN107450103A (zh) * | 2017-09-07 | 2017-12-08 | 中海石油(中国)有限公司 | 一种基于边界积分反算子的鬼波压制方法 |
| CN109143372A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 海洋上下缆拖缆宽频采集的缆深组合选取方法及系统 |
| CN109633752A (zh) * | 2019-01-04 | 2019-04-16 | 吉林大学 | 基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法 |
| CN111948711A (zh) * | 2020-08-06 | 2020-11-17 | 中海石油(中国)有限公司 | 一种利用深度学习法提取地震数据低频部分的方法和系统 |
| CN112817047A (zh) * | 2020-12-31 | 2021-05-18 | 北京东方联创地球物理技术有限公司 | 海洋地震自适应去鬼波方法、装置、电子设备及介质 |
| CN112835103A (zh) * | 2020-12-31 | 2021-05-25 | 北京东方联创地球物理技术有限公司 | 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统 |
| CN115980845A (zh) * | 2023-02-21 | 2023-04-18 | 中海石油(中国)有限公司 | 基于激发相位先验约束的子波提取方法、系统、设备 |
| CN116148928A (zh) * | 2023-02-16 | 2023-05-23 | 中国石油大学(北京) | 一种地震数据鬼波压制方法、装置、设备及介质 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9234977B2 (en) * | 2012-01-12 | 2016-01-12 | Westergeco L.L.C. | Processing collected survey data |
-
2023
- 2023-06-25 CN CN202310754376.0A patent/CN116774280B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103308944A (zh) * | 2013-05-20 | 2013-09-18 | 中国石油天然气股份有限公司 | 一种海洋地震资料处理中鬼波衰减的方法及装置 |
| CN109143372A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 海洋上下缆拖缆宽频采集的缆深组合选取方法及系统 |
| CN107450103A (zh) * | 2017-09-07 | 2017-12-08 | 中海石油(中国)有限公司 | 一种基于边界积分反算子的鬼波压制方法 |
| CN109633752A (zh) * | 2019-01-04 | 2019-04-16 | 吉林大学 | 基于三维快速Radon变换的海上拖缆资料自适应鬼波压制方法 |
| CN111948711A (zh) * | 2020-08-06 | 2020-11-17 | 中海石油(中国)有限公司 | 一种利用深度学习法提取地震数据低频部分的方法和系统 |
| CN112817047A (zh) * | 2020-12-31 | 2021-05-18 | 北京东方联创地球物理技术有限公司 | 海洋地震自适应去鬼波方法、装置、电子设备及介质 |
| CN112835103A (zh) * | 2020-12-31 | 2021-05-25 | 北京东方联创地球物理技术有限公司 | 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统 |
| CN116148928A (zh) * | 2023-02-16 | 2023-05-23 | 中国石油大学(北京) | 一种地震数据鬼波压制方法、装置、设备及介质 |
| CN115980845A (zh) * | 2023-02-21 | 2023-04-18 | 中海石油(中国)有限公司 | 基于激发相位先验约束的子波提取方法、系统、设备 |
Non-Patent Citations (4)
| Title |
|---|
| Ghost imaging with different frequencies through non-degenerated four-wave mixing;Yu, Y (Yu, Ya) 等;OPTICS EXPRESS;20160808;第24卷(第15期);全文 * |
| 南海深水盆地改善中深层地震成像关键技术研究;邓盾;裴健翔;邓勇;赫建伟;李潇潇;任婷;彭海龙;黎孝璋;;地球物理学进展;20181102(第02期);全文 * |
| 基于HHT方法的地层Q值补偿方法研究;颜中辉;徐华宁;李攀峰;刘鸿;杨佳佳;陈珊珊;杨传胜;;海洋地质与第四纪地质;20200228(第01期);全文 * |
| 基于波动方程的上下缆地震数据鬼波压制方法研究;管西竹;陈宝书;符力耘;陶杰;李列;;地球物理学报;20151015(第10期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116774280A (zh) | 2023-09-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7414907B2 (ja) | 事前訓練済みモデルの決定方法及びその決定装置、電子機器並びに記憶媒体 | |
| CN114242098A (zh) | 一种语音增强方法、装置、设备以及存储介质 | |
| CN114360562A (zh) | 语音处理方法、装置、电子设备和存储介质 | |
| CN117214950B (zh) | 一种多次波压制方法、装置、设备及存储介质 | |
| CN115270076A (zh) | 一种品质因子确定方法、装置、电子设备及存储介质 | |
| KR20220078538A (ko) | 생체 감지 모델 트레이닝 방법, 장치, 전자 기기 및 저장 매체 | |
| CN116774280B (zh) | 一种鬼波压制定量质控方法、装置、电子设备及存储介质 | |
| CN113657468A (zh) | 预训练模型的生成方法、装置、电子设备和存储介质 | |
| CN117786500A (zh) | 一种盾构施工地表沉降预测方法、装置、设备及存储介质 | |
| CN116230001A (zh) | 一种混合语音分离方法、装置、设备及存储介质 | |
| EP4064282A2 (en) | Time delay estimation for echo cancellation | |
| CN116363444A (zh) | 模糊分类模型训练方法、识别模糊图像的方法及装置 | |
| CN115599687A (zh) | 一种软件测试场景的确定方法、装置、设备及介质 | |
| CN114463215A (zh) | 一种图像去噪方法、装置、设备及存储介质 | |
| CN116736374B (zh) | 多方位地震数据的融合方法、装置、设备和介质 | |
| CN115267908A (zh) | 余弦相似度地震属性计算方法、装置、电子设备及介质 | |
| CN116468927A (zh) | 一种叠层石储层识别方法、装置、设备及存储介质 | |
| CN116719085B (zh) | 一种地震记录高分辨率处理方法、装置、设备及存储介质 | |
| CN117972469A (zh) | 一种火成岩分类识别方法、装置、设备及存储介质 | |
| CN117990378A (zh) | 一种发动机声品质评价方法、装置、电子设备及存储介质 | |
| CN118295009A (zh) | 速度谱处理方法、装置、设备及存储介质 | |
| CN117993478A (zh) | 基于双向知识蒸馏和联邦学习的模型训练方法及装置 | |
| CN116819627A (zh) | 地震弱信号的增强方法、装置、设备和介质 | |
| CN117988714A (zh) | 一种水平井的钻井轨迹确定方法、装置、设备及介质 | |
| CN114418084A (zh) | 非结构化剪枝模型获取方法、装置、电子设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |