CN111596358B - 一种多次波的压制方法、设备以及系统 - Google Patents
一种多次波的压制方法、设备以及系统 Download PDFInfo
- Publication number
- CN111596358B CN111596358B CN201910128499.7A CN201910128499A CN111596358B CN 111596358 B CN111596358 B CN 111596358B CN 201910128499 A CN201910128499 A CN 201910128499A CN 111596358 B CN111596358 B CN 111596358B
- Authority
- CN
- China
- Prior art keywords
- common
- gather
- module
- model
- determining
- 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
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000001629 suppression Effects 0.000 title claims description 10
- 238000012937 correction Methods 0.000 claims description 46
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 abstract description 9
- 230000010354 integration Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 21
- 238000003825 pressing Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 230000006872 improvement Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000013016 damping Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种多次波的压制方法、设备以及系统，涉及能源勘探技术领域。所述方法包括：获取地震数据；将所述地震数据分选到共中心点域，形成共中心点道集；根据所述共中心点道集确定共中心点道集的多次波最终模型；将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集。本发明采用地震地质一体化理念，对多次波进行多数据域多方法的组合去除，从而克服单一方法去噪不完善的局限性，实现保真去除多次波提高成像分辨率的目的。
Description
技术领域
本发明关于能源勘探技术领域，特别是关于地震资料的处理技术，具体的讲是一种多次波的压制方法、系统、计算机设备以及计算机可读存储介质。
背景技术
多次波是地震资料中普遍存在的一种干扰波，其类型多、波场复杂，并且与有效波有非常相似的性质。因此，多次波的压制一直是地震资料处理中的难题。
现有技术中，去除多次波时首先需对多次波进行识别。在地震资料处理工作中，可以从道集、速度谱、叠加剖面以及偏移剖面中进行识别。多次波识别的越清晰，其压制效果往往也越好。对于多次波的压制方法经过近几十年的研究得到了长足的发展，基本可分为如下三类：
①基于多次波的周期性来压制多次波，如共中心点叠加技术、拉东变换技术。
②基于多次波与有效波的速度差异压制多次波，如τ-p域预测反褶积技术。
③多次波模型减去法，如FK倾角滤波技术。
由于多次波的周期、速度都有时变、空变的特征，且大部分压制方法都是基于数据或模型驱动。上述单一方法的应用无法完整的压制多次波。
因此，如何研究和开发出一种试验方案，其能够完整的压制多次波是本领域亟待解决的技术难题。
发明内容
本发明提供了一种多次波的压制方法、设备以及系统，采用地震地质一体化理念，对多次波进行多数据域多方法的组合去除，从而克服单一方法去噪不完善的局限性，实现保真去除多次波提高成像分辨率的目的。
本发明的目的是，提供一种多次波的压制方法，所述方法包括：
获取地震数据；
将所述地震数据分选到共中心点域，形成共中心点道集；
根据所述共中心点道集确定共中心点道集的多次波最终模型；
将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集；
其中，根据所述共中心点道集确定共中心点道集的多次波最终模型包括：
根据所述共中心点道集确定共中心点道集的多次波初始模型；
根据所述共中心点道集的多次波初始模型确定多次波速度；
利用所述多次波速度对所述共中心点道集进行多次波动校正；
利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测，得到共中心点道集的多次波最终模型；
其中，根据所述共中心点道集确定共中心点道集的多次波初始模型包括：
利用一次波速度对所述共中心点道集进行动校正；
将动校正后的共中心点道集进行分选，形成共偏移距道集；
在所述共偏移距道集中，分析产生多次波的目标层；
根据所述共偏移距道集结合所述目标层确定目标区域；
利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多次波初始模型；
将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模型；
其中，根据所述共中心点道集的多次波初始模型确定多次波速度包括：
利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正；
根据反动校正后的共中心点道集的多次波初始模型制作速度谱；
根据所述速度谱进行多次波速度解释，得到多次波速度。
在本发明的优选实施方式中，根据所述共中心点道集确定共中心点道集的多次波初始模型还包括当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。
在本发明的优选实施方式中，根据所述共中心点道集确定共中心点道集的多次波最终模型还包括：对所述共中心点道集的多次波最终模型进行多次波反动校正。
本发明的目的是，提供一种多次波的压制系统，所述系统包括：
地震数据获取模块，用于获取地震数据；
地震数据分选模块，用于将所述地震数据分选到共中心点域，形成共中心点道集；
最终模型确定模块，用于根据所述共中心点道集确定共中心点道集的多次波最终模型；
多次波压制模块，用于将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集；
其中，所述最终模型确定模块包括：
初始模型确定模块，用于根据所述共中心点道集确定共中心点道集的多次波初始模型；
多次波速度确定模块，用于根据所述共中心点道集的多次波初始模型确定多次波速度；
多次波动校正模块，用于利用所述多次波速度对所述共中心点道集进行多次波动校正；
多次波预测模块，用于利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测，得到共中心点道集的多次波最终模型；
其中，所述初始模型确定模块包括：
动校正模块，用于利用一次波速度对所述共中心点道集进行动校正；
道集分选模块，用于将动校正后的共中心点道集进行分选，形成共偏移距道集；
目标层分析模块，用于在所述共偏移距道集中，分析产生多次波的目标层；
目标区域确定模块，用于根据所述共偏移距道集结合所述目标层确定目标区域；
模型确定模块，用于利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多次波初始模型；
初始模型分选模块，用于将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模型；
其中，所述多次波速度确定模块包括：
初始模型反动校正模块，用于利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正；
速度谱制作模块，用于根据反动校正后的共中心点道集的多次波初始模型制作速度谱；
速度解释模块，用于根据所述速度谱进行多次波速度解释，得到多次波速度。
在本发明的优选实施方式中，所述初始模型确定模块还包括：
拉平处理模块，用于当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。
在本发明的优选实施方式中，所述最终模型确定模块还包括：
最终模型反动校正模块，用于对所述共中心点道集的多次波最终模型进行多次波反动校正。
本发明的目的之一是，提供一种计算机设备，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种多次波的压制方法。
本发明的目的之一是，提供一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于执行一种多次波的压制方法。
本发明的有益效果在于，提供了一种多次波的压制方法、设备以及系统，采用地震地质一体化理念，对多次波进行多数据域多方法的组合去除，从而克服单一方法去噪不完善的局限性，实现保真去除多次波提高成像分辨率的目的，具有保真压制多次波的目的能力，在有效去除多次波干扰的前提下不损伤有效信号，大幅提高成像分辨率，对地质构造的刻画更加准确。
为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种多次波的压制方法的流程图；
图2为图1中的步骤S103的具体流程图；
图3为图2中的步骤S201的具体流程图；
图4为图2中的步骤S202的具体流程图；
图5为本发明实施例提供的一种多次波的压制系统的结构框图；
图6为本发明实施例提供的一种多次波的压制系统中的最终模型确定模块的结构框图；
图7为本发明实施例提供的一种多次波的压制系统中的初始模型确定模块的结构框图；
图8为本发明实施例提供的一种多次波的压制系统中的多次波速度确定模块的结构框图；
图9为本发明提供的具体实施例中地震数据分选形成CMP道集的示意图；
图10为本发明提供的具体实施例中一次波动校正后的CMP道集的示意图；
图11为本发明提供的具体实施例中地震数据分选形成共偏移距道集的示意图；
图12、图13为本发明提供的具体实施例中共偏移距道集多次波预测示意图；
图14为本发明提供的具体实施例中多次波速度解释的示意图；
图15为本发明提供的具体实施例中CMP域最终多次波模型的示意图；
图16为本发明提供的具体实施例中多次波压制前的CMP道集示意图；
图17为本发明提供的具体实施例中多次波压制后的CMP道集示意图；
图18为本发明提供的具体实施例中多次波压制前偏移剖面的示意图；
图19为本发明提供的具体实施例中多次波压制后偏移剖面的示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本领域技术技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
本发明提供了一种多次波的压制方法、设备以及系统。下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。以下所使用的术语“模块”和“单元”，可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
由于单一方法的应用无法完整的压制多次波，多方法联合使用已成为目前压制多次波的趋势，本发明提供的一种多次波的压制方法，请参阅图1，所述方法包括：
S101：获取地震数据。
S102：将所述地震数据分选到共中心点域，形成共中心点道集。
S103：根据所述共中心点道集确定共中心点道集的多次波最终模型；
S104：将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集。
图2为步骤S103的具体流程图，请参阅图2，该步骤包括：
S201：根据所述共中心点道集确定共中心点道集的多次波初始模型。
图3为该步骤的具体流程图，请参阅图3，根据所述共中心点道集确定共中心点道集的多次波初始模型在本发明的一种实施方式中包括：
S301：利用一次波速度对所述共中心点道集进行动校正；
S302：将动校正后的共中心点道集进行分选，形成共偏移距道集；
S303：在所述共偏移距道集中，分析产生多次波的目标层；
S304：根据所述共偏移距道集结合所述目标层确定目标区域；
S305：利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多次波初始模型；
S306：将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模型。
在本发明的一种实施方式中，步骤S201在步骤S303之后还包括当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。
请参阅图2，该步骤S103还包括：
S202：根据所述共中心点道集的多次波初始模型确定多次波速度。图4为步骤S202的具体流程图，请参阅图4，该步骤包括：
S401：利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正；
S402：根据反动校正后的共中心点道集的多次波初始模型制作速度谱；
S403：根据所述速度谱进行多次波速度解释，得到多次波速度。
请参阅图2，该步骤S103还包括：
S203：利用所述多次波速度对所述共中心点道集进行多次波动校正；
S204：利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测，得到共中心点道集的多次波最终模型。
在本发明的一种实施方式中，步骤S204在对多次波进行预测之后还包括对所述共中心点道集的多次波最终模型进行多次波反动校正。
在得到共中心点道集的多次波最终模型后，步骤S104将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集，从而达到压制多次波，提高地震数据成像精度的目的。
如上所述，即为本发明提供的一种多次波的压制方法，具有保真压制多次波的目的能力，在有效去除多次波干扰的前提下不损伤有效信号，大幅提高成像分辨率，对地质构造的刻画更加准确。
本发明还提供了一种计算机设备，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种多次波的压制方法。
本发明还提供了一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于执行一种多次波的压制方法。
应当注意，尽管在附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
在介绍了本发明提供的多次波的压制方法、计算机设备以及计算机可读存储介质的示例性实施方式之后，接下来，参考附图对本发明示例性实施方式的系统进行介绍。该系统的实施可以参见上述整体的实施，重复之处不再赘述。
本发明解决其技术问题所采用的方案是：首先在共偏移距面上进行第一次多次波预测，然后将预测的模型分选到共中心点域，然后进行第二次预测，得到最终多次波模型，最后对数据进行组合处理，去除多次波。
由于单一方法的应用无法完整的压制多次波，多方法联合使用已成为目前压制多次波的趋势，本发明提供的一种多次波的压制系统，请参阅图5，所述系统包括：
地震数据获取模块100，用于获取地震数据。
地震数据分选模块200，用于将所述地震数据分选到共中心点域，形成共中心点道集。
最终模型确定模块300，用于根据所述共中心点道集确定共中心点道集的多次波最终模型；
多次波压制模块400，用于将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集。
图6为本发明实施例提供的一种多次波的压制系统中的最终模型确定模块300的结构框图，请参阅图6，最终模型确定模块300包括：
初始模型确定模块301，用于根据所述共中心点道集确定共中心点道集的多次波初始模型。
图7为本发明实施例提供的一种多次波的压制系统中的初始模型确定模块301的结构框图，请参阅图3，初始模型确定模块301在本发明的一种实施方式中包括：
动校正模块3011，用于利用一次波速度对所述共中心点道集进行动校正；
道集分选模块3012，用于将动校正后的共中心点道集进行分选，形成共偏移距道集；
目标层分析模块3013，用于在所述共偏移距道集中，分析产生多次波的目标层；
目标区域确定模块3014，用于根据所述共偏移距道集结合所述目标层确定目标区域；
模型确定模块3015，用于利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多次波初始模型；
初始模型分选模块3016，用于将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模型。
在本发明的一种实施方式中，初始模型确定模块还包括拉平处理模块，用于当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。
请参阅图6，最终模型确定模块300还包括：
多次波速度确定模块302，用于根据所述共中心点道集的多次波初始模型确定多次波速度。图8为本发明实施例提供的一种多次波的压制系统中的多次波速度确定模块的结构框图，请参阅图8，多次波速度确定模块302包括：
初始模型反动校正模块3021，用于利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正；
速度谱制作模块3022，用于根据反动校正后的共中心点道集的多次波初始模型制作速度谱；
速度解释模块3023，用于根据所述速度谱进行多次波速度解释，得到多次波速度。
请参阅图6，最终模型确定模块还包括：
多次波动校正模块303，用于利用所述多次波速度对所述共中心点道集进行多次波动校正；
多次波预测模块304，用于利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测，得到共中心点道集的多次波最终模型。
在本发明的一种实施方式中，最终模型确定模块在对多次波进行预测之后还包括对所述共中心点道集的多次波最终模型进行多次波反动校正。
在得到共中心点道集的多次波最终模型后，多次波压制模块将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集，从而达到压制多次波，提高地震数据成像精度的目的。
如上所述，即为本发明提供的一种多次波的压制系统，具有保真压制多次波的目的能力，在有效去除多次波干扰的前提下不损伤有效信号，大幅提高成像分辨率，对地质构造的刻画更加准确。
此外，尽管在上文详细描述中提及了系统的若干单元模块，但是这种划分仅仅并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样，上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。以上所使用的术语“模块”和“单元”，可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
下面结合具体的实施例，详细介绍本发明的技术方案。
在该实施例中，采集的地震资料是在新疆油田车排子凸起-四棵树凹陷地区采集的。图9为该具体实施例中对地震数据分选后形成CMP道集的示意图，图10为该具体实施例中一次波动校正后的CMP道集的示意图，图11为该具体实施例中地震数据分选形成共偏移距道集的示意图，图12、图13为该具体实施例中共偏移距道集多次波预测示意图，图14为该具体实施例中多次波速度解释的示意图，图15为该具体实施例中CMP域最终多次波模型的示意图，图16为该具体实施例中多次波压制前的CMP道集示意图，图17为该具体实施例中多次波压制后的CMP道集示意图，图18为该具体实施例中多次波压制前偏移剖面的示意图，图19为该具体实施例中多次波压制后偏移剖面的示意图。图12、图13中的白色的虚线即为多次波。
参阅图16至图19可知，在新疆油田车排子凸起-四棵树凹陷地区的地震资料处理的过程中，成功的保真去除了石炭系标志层下的多次波，有效改善了该地区地震数据的成像质量，辅助地震地质解释专家更好的认识地质构造。图16、图17是应用本发明前后的道集，图18、图19是应用本发明前后的偏移剖面，可以看出石炭系内幕多次波得到较好的压制，石炭系内幕一次波成像清楚。
对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(ProgrammableLogic Device,PLD)(例如现场可编程门阵列(Field Programmable Gate Array，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(HardwareDescription Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(AdvancedBoolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(JavaHardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby HardwareDescription Language)等，目前最普遍使用的是VHDL(Very-High-Speed IntegratedCircuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (8)
1.一种多次波的压制方法，其特征在于，所述方法包括：
获取地震数据；
将所述地震数据分选到共中心点域，形成共中心点道集；
根据所述共中心点道集确定共中心点道集的多次波最终模型；
将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集；
其中，根据所述共中心点道集确定共中心点道集的多次波最终模型，包括：
根据所述共中心点道集确定共中心点道集的多次波初始模型；
根据所述共中心点道集的多次波初始模型确定多次波速度；
利用所述多次波速度对所述共中心点道集进行多次波动校正；
利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测，得到共中心点道集的多次波最终模型；
其中，根据所述共中心点道集确定共中心点道集的多次波初始模型包括：
利用一次波速度对所述共中心点道集进行动校正；
将动校正后的共中心点道集进行分选，形成共偏移距道集；
在所述共偏移距道集中，分析产生多次波的目标层；
根据所述共偏移距道集结合所述目标层确定目标区域；
利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多次波初始模型；
将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模型；
其中，根据所述共中心点道集的多次波初始模型确定多次波速度包括：
利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正；
根据反动校正后的共中心点道集的多次波初始模型制作速度谱；
根据所述速度谱进行多次波速度解释，得到多次波速度。
2.根据权利要求1所述的方法，其特征在于，根据所述共中心点道集确定共中心点道集的多次波初始模型还包括：
当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。
3.根据权利要求1所述的方法，其特征在于，根据所述共中心点道集确定共中心点道集的多次波最终模型还包括：对所述共中心点道集的多次波最终模型进行多次波反动校正。
4.一种多次波的压制系统，其特征在于，所述系统包括：
地震数据获取模块，用于获取地震数据；
地震数据分选模块，用于将所述地震数据分选到共中心点域，形成共中心点道集；
最终模型确定模块，用于根据所述共中心点道集确定共中心点道集的多次波最终模型；
多次波压制模块，用于将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除多次波后的共中心点道集；
其中，所述最终模型确定模块包括：
初始模型确定模块，用于根据所述共中心点道集确定共中心点道集的多次波初始模型；
多次波速度确定模块，用于根据所述共中心点道集的多次波初始模型确定多次波速度；
多次波动校正模块，用于利用所述多次波速度对所述共中心点道集进行多次波动校正；
多次波预测模块，用于利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测，得到共中心点道集的多次波最终模型；
其中，所述初始模型确定模块包括：
动校正模块，用于利用一次波速度对所述共中心点道集进行动校正；
道集分选模块，用于将动校正后的共中心点道集进行分选，形成共偏移距道集；
目标层分析模块，用于在所述共偏移距道集中，分析产生多次波的目标层；
目标区域确定模块，用于根据所述共偏移距道集结合所述目标层确定目标区域；
模型确定模块，用于利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多次波初始模型；
初始模型分选模块，用于将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模型；
其中，所述多次波速度确定模块包括：
初始模型反动校正模块，用于利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正；
速度谱制作模块，用于根据反动校正后的共中心点道集的多次波初始模型制作速度谱；
速度解释模块，用于根据所述速度谱进行多次波速度解释，得到多次波速度。
5.根据权利要求4所述的系统，其特征在于，所述初始模型确定模块还包括：
拉平处理模块，用于当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。
6.根据权利要求4所述的系统，其特征在于，所述最终模型确定模块还包括：最终模型反动校正模块，用于对所述共中心点道集的多次波最终模型进行多次波反动校正。
7.一种计算机设备，其特征在于，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行如权利要求1至3任意一项所述的一种多次波的压制方法。
8.一种计算机可读存储介质，其特征在于，存储有计算机程序，所述计算机程序用于执行如权利要求1至3任意一项所述的一种多次波的压制方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910128499.7A CN111596358B (zh) | 2019-02-21 | 2019-02-21 | 一种多次波的压制方法、设备以及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910128499.7A CN111596358B (zh) | 2019-02-21 | 2019-02-21 | 一种多次波的压制方法、设备以及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111596358A CN111596358A (zh) | 2020-08-28 |
| CN111596358B true CN111596358B (zh) | 2023-04-07 |
Family
ID=72188561
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910128499.7A Active CN111596358B (zh) | 2019-02-21 | 2019-02-21 | 一种多次波的压制方法、设备以及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111596358B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112327362B (zh) * | 2020-10-30 | 2021-11-23 | 中国海洋大学 | 速度域的海底多次波预测与追踪衰减方法 |
| CN113514888B (zh) * | 2021-06-18 | 2024-04-16 | 中国石油化工股份有限公司 | 一种基于非剩余时差的多次波压制方法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4209854A (en) * | 1978-05-04 | 1980-06-24 | Mobil Oil Corporation | Method and system for enhancing primary seismic reflection signals in common depth point sets |
| US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
| US8300498B2 (en) * | 2009-06-30 | 2012-10-30 | Pgs Geophysical As | Method for dynamic aperture determination for three-dimensional surface-related multiple elimination |
| US9170345B2 (en) * | 2011-10-19 | 2015-10-27 | Exxonmobil Upstream Research Company | Near-offset extrapolation for free-surface multiple elimination in shallow marine environment |
| CN102636811B (zh) * | 2012-04-10 | 2014-01-29 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种海上二维地震资料中多次波的消除方法 |
| GB2501934B (en) * | 2012-05-11 | 2014-07-09 | Statoil Petroleum As | Method and apparatus for processing seismic data |
| CN104977615B (zh) * | 2015-07-01 | 2017-06-23 | 中国海洋石油总公司 | 一种基于模型统计拾取的深水obc资料多次波压制方法 |
| CN106353817A (zh) * | 2016-08-11 | 2017-01-25 | 北京软岛科技有限公司 | 一种地震勘探数据多次波自动搜索和压制方法 |
| CN107656313B (zh) * | 2017-10-12 | 2019-08-06 | 中国石油天然气集团公司 | 一种拉东域多次波压制方法及装置 |
-
2019
- 2019-02-21 CN CN201910128499.7A patent/CN111596358B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111596358A (zh) | 2020-08-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bell et al. | Three-dimensional adaptive mesh refinement for hyperbolic conservation laws | |
| CN108107474B (zh) | 一种基于稀疏反演的混叠数据分离方法及装置 | |
| CN109471171B (zh) | 一种混叠地震数据分离的方法、装置及系统 | |
| CN111596358B (zh) | 一种多次波的压制方法、设备以及系统 | |
| CN107563297B (zh) | 一种波形匹配方法及装置 | |
| CN110082823B (zh) | 一种地震数据插值方法及装置 | |
| JP2022552077A (ja) | 強化された同時並行性ガーベッジ収集スタック走査 | |
| CN107063306A (zh) | 一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 | |
| CN103557856A (zh) | 一种光纤陀螺随机漂移实时滤波方法 | |
| CN108181644B (zh) | 一种采样方法及装置 | |
| CN109061724A (zh) | 一种基于自适应变分模态分解的地震数据降噪方法 | |
| CN104391325B (zh) | 不连续非均质地质体检测方法和装置 | |
| CN115267912A (zh) | 基于压缩感知的avo特征保持去噪方法及系统 | |
| CN110068860B (zh) | 一种时移地震数据处理方法、装置及系统 | |
| CN107561588B (zh) | 一种地震数据噪声压制方法及装置 | |
| US20190331721A1 (en) | Noise spectrum analysis for electronic device | |
| CN112434478B (zh) | 仿真逻辑系统设计的虚拟接口的方法及相关设备 | |
| CN108035708A (zh) | 一种去除地层界面反射波的方法及装置 | |
| CN110618463A (zh) | 一种核磁共振数据反演方法、装置、存储介质及设备 | |
| CN111239832A (zh) | 3d3c-vsp成像处理方法、装置及设备 | |
| Zhang et al. | Delay‐free tracking differentiator design based on variational mode decomposition: Application on MEMS gyroscope denoising | |
| CN109917459A (zh) | 一种压制地震噪声的方法、装置及系统 | |
| CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
| CN109471173B (zh) | 一种剩余静校正方法及装置 | |
| CN106094026A (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 |