CN107942387A - 一种续至波衰减处理方法及装置 - Google Patents
一种续至波衰减处理方法及装置 Download PDFInfo
- Publication number
- CN107942387A CN107942387A CN201710990343.0A CN201710990343A CN107942387A CN 107942387 A CN107942387 A CN 107942387A CN 201710990343 A CN201710990343 A CN 201710990343A CN 107942387 A CN107942387 A CN 107942387A
- Authority
- CN
- China
- Prior art keywords
- data
- noise
- later arrivals
- processing
- layer position
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 116
- 238000000034 method Methods 0.000 claims abstract description 80
- 230000008569 process Effects 0.000 claims abstract description 19
- 230000002441 reversible effect Effects 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 10
- 238000013519 translation Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 33
- 230000000694 effects Effects 0.000 description 15
- 238000003860 storage Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 14
- 238000011161 development Methods 0.000 description 11
- 230000006872 improvement Effects 0.000 description 11
- 239000011435 rock Substances 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000004590 computer program Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 239000000047 product Substances 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000004480 active ingredient Substances 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000003610 charcoal Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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
-
- 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
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供一种续至波衰减处理方法及装置。所述方法包括：根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；根据所述续至波数据拾取目标层位；根据所述目标层位和所述续至波数据，建立地层模型；根据所述地层模型，进行噪音预测处理，获得预测噪音数据；将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。利用本申请中各个实施例，衰减了地震数据中的续至波噪音，提高了地震数据处理和解释的准确性。
Description
技术领域
本申请属于石油勘探技术领域，尤其涉及一种续至波衰减处理方法及装置。
背景技术
现有技术中，可以利用测井技术和地震资料，进行油气的勘探。具体可以根据地震资料，获取地震波的剖面图，根据剖面图，分析地质结构，进行油气的勘探。通过对不同盆地地震资料进行处理解释分析发现，原始地震资料中间在强波阻抗界面下伏地层中，往往发育一套与强界面反射地层产状相同、频率较低、振幅较强的地震波。该地震波是由续至波产生，由续至波带来的地震波影响下伏地层的真实成像，而这些不整合界面附近往往是油气聚集的部位。地质人员需要准确识别每一套地震反射的地质含义，进而落实油藏的类型和规模。但是，续至波的影响，给地质含义的确定带来一定的难度。因此，地震资料的处理中，续至波的处理是一项关键的技术。
在现有技术中，由于处理人员缺乏对地区地质背景认识，难以采取针对性的措施。而且在现有技术的处理流程，如：提频处理，容易使这些低频续至波演变为高频噪音，与真实地层表现为串层现象，使处理解释人员更加难以判别。或者对这些低频噪音进行整体低频滤波，能达到一定的衰减效果，但是容易对深层地层低频有效成分造成损失，因此两种方法都不能达到衰减续至波噪音保留有效信号的目的。因此，业内亟需一种能够衰减待分析地震数据中的续至波，进一步提高地震数据的处理精度的实施方案。
发明内容
本申请目的在于提供一种续至波衰减处理方法及装置，通过层位解释建立地层模型，对频率特征和层位产状进行约束，有非常强的针对性效果，达到衰减强波阻抗界面产生的续至波噪音的效果。
一方面本申请提供了一种续至波衰减处理方法，包括：
根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；
根据所述续至波数据拾取目标层位；
根据所述目标层位和所述续至波数据，建立地层模型；
根据所述地层模型，进行噪音预测处理，获得预测噪音数据；
将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；
根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
进一步地，所述方法的另一个实施例中，所述根据所述目标层位和所述续至波数据，建立地层模型，包括：
拾取所述目标层位，将所述目标层位的数据进行拉平处理，获取所述目标层位的拉平层位数据；
给定空间时窗，根据所述拉平层位数据和所述续至波数据，获取所述空间时窗内的波形特征，建立所述地层模型，所述地层模型包括所述波形的参数信息。
进一步地，所述方法的另一个实施例中，所述根据所述地层模型，进行噪音预测处理，获得预测噪音数据，包括：
将所述地层模型中的数据由时间空间域转化为频率空间域；
根据所述地层模型，在频率域内利用最小二乘反演法，进行噪音预测滤波处理，获得所述预测噪音数据。
进一步地，所述方法的另一个实施例中，所述方法还包括：
获得噪音衰减地震数据后，将所述噪音衰减地震数据进行反拉平处理。
进一步地，所述方法的另一个实施例中，所述方法还包括：
判断所述噪音衰减地震数据是否满足预设的衰减要求，若不满足，则调整所述噪音预测数据，根据调整后的噪音预测数据进行续至波噪音的衰减处理。
另一方面，本申请提供了一种续至波衰减处理装置，包括：
续至波识别模块，用于根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；
目标层位拾取模块，用于根据所述续至波数据拾取目标层位；
地层模型建立模块，用于拾取目标层位，根据所述目标层位和所述续至波数据，建立地层模型；
噪音预测模块，用于根据所述地层模型，进行噪音预测处理，获得预测噪音数据；
噪音匹配模块，用于将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；
噪音衰减模块，用于根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
进一步地，所述装置的另一实施例中，所述地层模型建立模块包括：
数据拉平模块，用于将所述目标层位的数据进行拉平处理，获取所述目标层位的拉平层位数据；
数据处理模块，用于给定空间时窗，根据所述拉平层位数据和所述续至波数据，获取所述空间时窗内的波形特征，建立所述地层模型，所述地层模型包括所述波形的参数信息。
进一步地，所述装置的另一个实施例中，所述噪音预测模块包括：
频率转换模块，用于将所述地层模型中的数据由时间空间域转化为频率空间域；
反演预测模块，用于根据所述地层模型，在频率域内利用最小二乘反演法，进行噪音预测滤波处理，获得所述预测噪音数据。
进一步地，所述装置的另一个实施例中，所述续至波衰减处理装置还包括：
反拉平处理模块，用于在所述噪音衰减模块获得噪音衰减地震数据后，将所述噪音衰减地震数据进行反拉平处理。
进一步地，所述装置的另一个实施例中，所述噪音衰减模块还用于：
判断所述噪音衰减地震数据是否满足预设的衰减要求，若不满足，则调整所述噪音预测数据，根据调整后的噪音预测数据进行续至波噪音的衰减处理。
再一方面，本申请还提供了一种续至波衰减处理装置，包括：处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述续至波衰减处理方法。
本申请提供的续至波衰减处理方法及装置，根据地震资料、测井数据、地质分层数据等，识别出待分析地震数据中的续至波数据。并进行目标层位的拾取，在目标层位的基础上，结合获得的续至波数据的发育特征，进行噪音预测，获得预测噪音数据。根据预测噪音数据和待分析地震数据，可以进行续至波噪音的衰减处理。并在获得续至波噪音衰减后的地震数据后，进行反拉平操作，确保续至波噪音衰减后的地震数据保持原始的待分析地震数据的特征，仅仅只是少了续至波噪音数据。为后续地震数据的处理和解释提供了准确的数据基础。通过层位解释，对频率特征和层位产状进行约束，有非常强的针对性效果，达到衰减强波阻抗界面产生的续至波噪音的效果。衰减了续至波后的地震剖面，火成岩内幕反射结构更加清楚，地层接触关系、地层尖灭点更加清楚，石炭系火成岩体资料品质的改善。有利于火山岩体的边界刻画和油藏认识，衰减了地震假象给油气藏认识带来的影响，使得地震数据更加真实反映地下地质的变化，提高了地震数据处理的准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是现有技术中西部某盆地地区A中三叠系与石炭系接触的部位的剖面示意图；
图2是现有技术中西部某盆地地区B中二叠系与石炭系地层不整合界面剖面示意图；
图3是本申请提供的一种续至波衰减处理方法一个实施例的方法流程示意图；
图4是本申请一个实施例中通过测井标定后沿层拾取的目标层位示意图；
图5是本申请一个实施例中目标层位拉平处理后获得的拉平目标层位示意图；
图6是本申请一个实施例中续至波噪音衰减后的剖面示意图；
图7是本申请一个实施例中基于层位模拟的预测噪音干扰波示意图；
图8(a)-图8(b)是本申请一个实施例中续至波衰减前后剖面对比示意图；
图9是本申请一个实施例中测井标定结果图对比示意图；
图10(a)-图10(b)是本申请一个实施例中续至波衰减前后效果对比示意图；
图11(a)-图11(b)是本申请又一实施例中续至波衰减前后效果对比示意图；
图12是本申请提供的续至波衰减处理装置一个实施例的模块结构示意图；
图13是本申请一个实施例中地层模型建立模块的结构示意图；
图14是本申请一个实施例中噪音预测模块的结构示意图；
图15是本申请又一个实施例中的续至波衰减处理装置的结构示意图；
图16是本申请提供的另一种续至波衰减处理装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
现有技术中，通常可以利用测井技术和地震资料，进行油气的勘探。具体可以根据地震资料和测井资料，获取地震波的剖面图，根据剖面图，分析地质结构，进行油气的勘探。但是，由于地震资料是混合相位，本身不满足最小相位要求。在经过预测反褶积处理后，会带有一定延续时间的低频伴随相位，与不整合强波阻抗界面下伏地层反射波发生干涉。产生具有频率低、振幅强、与不整合界面地层产状一致等特点的地震波，即续至波。
续至波是一种噪音没具体的地质含义，在测井曲线上也没有地球物理响应，对于地质含义的解释带来一定的影响。例如：图1是现有技术中西部某盆地地区A中三叠系与石炭系接触的部位的剖面示意图，图2是现有技术中西部某盆地地区B中二叠系与石炭系地层不整合界面剖面示意图。图1中的T、P、C；分别表示三叠系、二叠系、石炭系三个不同地质时期，纵轴表示时间，横轴CDP(Common Depth Point)表示共深度点。图1的低频噪音表现特征较为突出，石炭系顶界面反射非常清楚，地质解释人员较容易识别地震假象。但是，图2中单凭地震资料认识和测井标定难以准确确定石炭系顶界与地震反射的对应关系，进一步难以进行油气的勘探。图2中纵轴LINE CMP(Common Middle Point)表示共中心点道集测线。
续至波的存在使得地震数据无法准确的反应地质的真实含义，本申请提供的一个或多个实施例，可以基于层位约束，通过层位解释建立地层模型，对频率特征和层位产状进行约束。达到衰减强波阻抗界面产生的续至波噪音的效果，为后续地质含义的解释，完成油气勘探提供准确的数据基础。
图3是本申请提供的一种续至波衰减处理方法一个实施例的方法流程示意图，本申请提供的续至波衰减处理方法包括：
S1、根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据。
具体地，续至波是由于上覆层强波阻抗界面经过反褶积后的一种低频伴随相位，与下伏层干涉后产生的具有频率低、振幅强、和不整合界面地层产状一致等特点的地震波。续至波是一种噪音，没有具体的地质含义，在测井曲线上也没有地球物理响应。本实施例中可以通过获得的测井数据如：测井曲线，进行声波和密度与地震子波做褶积得到合成记录，再将得到的合成记录与地震资料进行互相关匹配。再根据地质分层数据，可以通过标定识别出待分析地震数据中存在的噪音。同时，本实施例可以通过频谱扫描和时频分析确定有效信号和噪音的频带范围，如：有效信号的频率为6-60Hz，噪音的频率范围为6-12Hz。有效信号是指可以包括地质层位的反射界面同相轴。续至波是一种噪音，将噪音进行归类处理，可以识别出续至波，也可以获取续至波在待分析地震数据中位置。并可以获取续至波的发育特征，例如：续至波的振幅强弱、频率高低、产状等。
S2、根据所述续至波数据拾取目标层位。
具体地，续至波是由于上覆层强波阻抗界面经过反褶积后的一种低频伴随相位，与下伏层干涉后产生的地震波，通过地震资料的解释可以拾取目标层位。目标层位可以包括在噪音发育的上覆层有确切的地质代表含义且波阻抗界面清楚的特征的地质层，在待分析地震数据体上表现为强波峰或强波谷特征，且地层的产状与噪音的产状基本一致的地质层。本实施例中，所说的拾取目标层具体的可以包括将目标层位在空间上连续的一组最大波峰值或者波谷值记录下来，可以将这一组数值(x，y，t0)定义成目标层位。
S3、根据所述目标层位和所述续至波数据，建立地层模型。
具体地，在目标层位的基础上，可以根据续至波数据建立地层模型。地层模型可以包括：某一范围内的同相轴的振幅、频率等物理特征，具体可以参考续至波数据的发育特征获取目标层位中的波形特征，建立地层模型。地层模型可以是对地层实体的组成、结构、时空存在状况的简化表达和综合解释，具体可以是用一定垂向地层间隔的代表性基本层序表示该间隔地层组成、结构、沉积环境的剖面模型，或根据续至波数据建立的可以反映地质结构特征的数据图像模型。
S4、根据所述地层模型，进行噪音预测处理，获得预测噪音数据。
具体地，可以根据建立的地层模型，结合强波阻抗界面的产状、振幅、频率特征进行噪音预测，获得噪音预测数据。噪音预测数据可以包括与目标层位产状、振幅特征相似的同相轴。
S5、将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据。
具体地，本申请一个实施例可以将预测噪音数据和待分析地震数据进行匹配，具体可以利用最小二乘滤波将预测噪音数据和待分析地震数据进行匹配，获得匹配噪音数据。原始的待分析地震数据与预测噪音数据因为振幅、相位不匹配不能直接相减。本申请进行预测噪音的匹配，使得匹配后的匹配噪音数据能够与原始的待分析地震数据直接相减，实现续至波噪音的衰减。
S6、根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
具体地，续至波是一种噪音，续至波衰减处理又可以理解为地震数据的噪音衰减处理。获得匹配噪音数据后，可以将待分析地震数据和匹配噪音数据利用自适应法进行相减，使得匹配噪音数据从待分析地震数据中删除，获得噪音衰减地震数据，实现待分析地震数据的续至波衰减处理。
在本申请一个实施例中，在获得噪音衰减地震数据后，还可以判断获得的噪音衰减地震数据是否满足预设的要求。若不满足，可以通过分析噪音衰减地震数据中的噪音特征，通过调整预测噪音数据，重新进行续至波的噪音衰减操作，直至获得的噪音衰减地震数据满足预设的要求。预设的要求可以根据实际需要进行设置，如：可以根据获得的噪音衰减地震数据中噪音数据的频率是否满足预设的频率阈值等进行设置，本申请不做具体限定。
本申请提供的续至波衰减处理方法，根据地震资料、测井数据、地质分层数据等，识别出待分析地震数据中的续至波数据。并进行目标层位的拾取，在目标层位的基础上，结合获得的续至波数据的发育特征，进行噪音预测，获得预测噪音数据。根据预测噪音数据和待分析地震数据，可以进行续至波噪音的衰减处理。通过层位解释，对频率特征和层位产状进行约束，有非常强的针对性效果，达到衰减强波阻抗界面产生的续至波噪音的效果。衰减了续至波后的地震剖面，火成岩内幕反射结构更加清楚，地层接触关系、地层尖灭点更加清楚，石炭系火成岩体资料品质的改善。有利于火山岩体的边界刻画和油藏认识，衰减了地震假象给油气藏认识带来的影响，使得地震数据更加真实反映地下地质的变化，提高了地震数据处理的准确性。
在上述实施例的基础上，所述根据所述目标层位和所述续至波数据，建立地层模型，包括：
将所述目标层位的数据进行拉平处理，获取所述目标层位的拉平层位数据；
给定空间时窗，根据所述拉平层位数据和所述续至波数据，获取所述空间时窗内的波形特征，建立所述地层模型，所述地层模型包括所述波形的参数信息。
具体地，由于实际的层位在空间上的时间是变化的，后期处理比较复杂。本申请一个实施例中可以通过定义一个固定时间值作为参考面，参考时间与实际数据层位的时间值作为时间差。将目标层位的时间值通过时间差值调整到统一的参考面，即可实现目标层位的样点时间为一常数(也就是拉平层位)。即将目标层位的数据进行拉平处理，获取到目标层位的拉平层位数据。空间时窗可以包括对拉平后的数据体即拉平层位数据，在目标层上下定义时窗，时窗的顶、底尽量包全目标层位。给定空间时窗的作用是限定后面建立的地层模型的建立的空间大小，即只在这一空间时窗范围内进行特征分析，建立地层模型，缩小数据处理的范围，提高数据处理的速度。具体可以在给定的空间时窗内，根据拉平层位数据和续至波数据，获取空间时窗内与续至波数据的发育特征相似的波形的特征数据，构建地层模型。地层模型的具体形式可以包括：波形的空间坐标、时窗的时间顶底值、目标层的振幅值、频率值以及目标层横向连续一致性等。地层模型可以限定某一范围并在这一范围类统计同相轴的振幅、频率物理特征。
图4是本申请一个实施例中通过测井标定后沿层拾取的目标层位示意图，图5是本申请一个实施例中目标层位拉平处理后获得的拉平目标层位示意图。如图4和图5所示，目标层位进行数据拉平处理后，层位剖面数据在时间上更加整齐。
本申请提供的续至波衰减处理方法，可以通过给定空间时窗，在目标层位上结合获得续至波数据的发育特征，在空间时窗内进行特征统计，快速的构建地层模型。提高了续至波衰减处理的速度，同时为后续进行续至波噪音的衰减提供了准确的数据基础。
在上述实施例的基础上，所述根据所述地层模型，进行噪音预测处理，获得预测噪音数据，包括：
将所述地层模型中的数据由时间空间域转化为频率空间域；
根据所述地层模型，在频率域内利用最小二乘反演法，进行噪音预测滤波处理，获得所述预测噪音数据。
具体地，可以将获得地层模型中的数据由时间空间域转换成频率空间域，在频率域进行噪音预测滤波。可以利用最小二乘反演预测滤波，具体可以利用常规预测方法求解预测算子，用最小二乘法优化滤波目标函数反演求解信号，再迭代求出新的反演预测算子，利用新的预测算子求出反演信号，将反演信号与实际数据匹配相减得到预测噪音数据。
本申请提供的续至波衰减处理方法，可以采用最小二乘反演法进行噪音预测，获得预测噪音数据。为后续进行续至波噪音的衰减提供了准确的数据基础，提高了地震数据处理结果的准确性。
在上述实施例的基础上，所述方法还包括：
获得噪音衰减地震数据后，将所述噪音衰减地震数据进行反拉平处理。
通常的，由于在建立地层模型时，将目标层位的数据进行了拉平处理。本申请一个实施了中，在获得噪音衰减地震数据后，将进行续至波噪音衰减后的噪音衰减地震数据进行了反拉平处理。将目标层位数据拉平处理时得到的调整时间差值反应用到拉平数据体上，可以将去噪后的数据恢复成处理前的特征。保证获得噪音衰减地震数据与处理前的待分析地震数据的差别只在于是否包含续至波噪音，保证了待分析地震数据处理的准确性。
在本申请一个示例中，续至波衰减处理方法可以参考如下的实施过程：
1)加载数据。本申请中加载的数据可以包括常规偏移处理完的共反射点叠加数据。
2)识别续至波。井震结合，即测井资料和地震资料结合进行数据处理识别续至波的具体位置和待分析地震数据中续至波的发育特征。即可以利用测井曲线和地震资料标定地质层位，确定续至波发育时间位置和续至波的发育特征等。
3)拾取目标层位，拉平目标层位的待分析地震数据体，建立地层模型。具体可以通过拾取的层位时间t0信息，通过修改道头值方式，将拾取时间值赋予数据体(数据体是指整个流程中处理数据)的低频量。然后通过进行时间域校正，从而实现数据体的拉平，获得目标层位的拉平层位数据。指定空间时窗，在指定的空间时窗内，根据获取的续至波数据，拾取目标层位强波阻抗界面反射波的最大波峰或者波谷值，建立地层模型。
4)续至波的预测。根据地层模型，即强波阻抗界面的产状、振幅、频率特征进行噪音预测(预测的是与目标层位产状、振幅特征相似的同相轴)，获得噪音预测数据。噪音预测数据的主要参数可以包括横向选取道数NT、空间重叠道数TOL、最大最小预测频率Fmax、Fmin等。
5)预测噪音数据与待分析地震数据进行匹配，通过最小二乘滤波法匹配原始数据衰减噪音。原始的待分析地震数据与预测噪音数据相减，因为振幅、相位二者不匹配不能直接相减，因此必须将预测噪音数据做最小二乘法匹配后才能相减。
6)将处理后的数据反拉平及恢复原始数据特征。
7)对比处理前后效果，根据残余的噪音特征，可以通过修改预测噪音数据的最大最小预测频率Fmax、Fmin和空间选取道数NT重复上述过程，直到处理结果满意，实现待分析地震数据的续至波噪音衰减。
图6是本申请一个实施例中续至波噪音衰减后的剖面示意图，图7是本申请一个实施例中基于层位模拟的预测噪音干扰波示意图。图8(a)-图8(b)是本申请一个实施例中续至波衰减前后剖面对比示意图，图8(a)是本申请一个实施例中续至波衰减前剖面示意图，图8(b)是本申请一个实施例中续至波衰减后剖面示意图。图9是本申请一个实施例中测井标定结果图对比示意图，图9左边部分是本申请一个实施例中续至波衰减前标定剖面示意图，右边部分是本申请一个实施例中续至波衰减后标定剖面示意图。图10(a)-图10(b)是本申请一个实施例中续至波衰减前后效果对比示意图，图10(a)是续至波衰减前剖面示意图，图10(b)是续至波衰减后剖面示意图。图11(a)-图11(b)是本申请又一实施例中续至波衰减前后效果对比示意图，图11(a)是续至波衰减前剖面示意图，图11(b)是续至波衰减后剖面示意图。上述各图中LINECMP表示供中心点道集测线。通过上述各个示意图的对比可以发现，本申请的续至波衰减方法，强反射续至波得到很好压制且其他产状有效信号成分基本没有损伤。
本申请提供的续至波衰减处理方法，可以根据测井数据、地质分层数据和地震资料，分析获得待分析地震数据中的续至波数据。根据获得的续至波数据，进行目标层位当然拾取，并通过设定空间时窗，目标层位的拉平处理等操作，建立地层模型。根据地层模型，利用最小二乘反演法进行噪音预测，获得预测噪音数据。根据预测噪音数据进行续至波噪音的衰减处理。将待分析地震数据中的续至波进行衰减处理，且只衰减给定频率端并且与上覆层位产状一致的续至波噪音，其他信息保持不变。提高了地震数据的精度，避免因为地震数据存在续至波噪音，影响地质解释人员对地震数据的分析。
基于上述所述的续至波衰减处理方法，本说明书一个或多个实施例还提供一种续至波衰减处理装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
具体地，图12是本申请提供的续至波衰减处理装置一个实施例的模块结构示意图，如图12所示，本申请中提供的续至波衰减处理装置包括：续至波识别模块121、目标层位拾取模块122，地层模型建立模块123、噪音预测模块124、噪音匹配模块125、噪音衰减模块126。
续至波识别模块121，可以用于根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；
目标层位拾取模块122，可以用于根据所述续至波数据拾取目标层位；
地层模型建立模块123，可以用于拾取目标层位，根据所述目标层位和所述续至波数据，建立地层模型；
噪音预测模块124，可以用于根据所述地层模型，进行噪音预测处理，获得预测噪音数据；
噪音匹配模块125，可以用于将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；
噪音衰减模块126，可以用于根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
本申请提供的续至波衰减处理装置，根据地震资料、测井数据、地质分层数据等，识别出待分析地震数据中的续至波数据。并进行目标层位的拾取，在目标层位的基础上，结合获得的续至波数据的发育特征，进行噪音预测，获得预测噪音数据。根据预测噪音数据和待分析地震数据，可以进行续至波噪音的衰减处理。通过层位解释，对频率特征和层位产状进行约束，有非常强的针对性效果，达到衰减强波阻抗界面产生的续至波噪音的效果。衰减了续至波后的地震剖面，使得火成岩内幕反射结构更加清楚，地层接触关系、地层尖灭点更加清楚，石炭系火成岩体资料品质的改善。有利于火山岩体的边界刻画和油藏认识，衰减地震假象给油气藏认识带来的影响，使得待分析地震数据更加真实反映地下地质的变化，提高了待分析地震数据处理的准确性。
图13是本申请一个实施例中地层模型建立模块的结构示意图，如图13所示，在上述实施例的基础上，所述地层模型建立模块123包括：
数据拉平模块131，可以用于将所述目标层位的数据进行拉平处理，获取所述目标层位的拉平层位数据；
数据处理模块132，可以用于给定空间时窗，根据所述拉平层位数据和所述续至波数据，获取所述空间时窗内的波形特征，建立所述地层模型，所述地层模型包括所述波形的参数信息。
本申请提供的续至波衰减处理装置，通过拾取目标层位、拉平数据处理，在给定的空间时窗内对拉平后的数据进行数据的分析，构建地层模型。为后续续至波噪音的衰减提供了准确的数据基础，提高了地震数据处理的精度。
图14是本申请一个实施例中噪音预测模块的结构示意图，如图14所示，在上述实施例的基础上，所述噪音预测模块124包括：
频率转换模块141，可以用于将所述地层模型中的数据由时间空间域转化为频率空间域；
反演预测模块142，可以用于根据所述地层模型，在频率域内利用最小二乘反演法，进行噪音预测滤波处理，获得所述预测噪音数据。
本申请提供的续至波衰减处理装置，根据获得的地层模型和续至波数据的特征，通过最小二乘反演法，进行噪音预测，可以获得准确的预测噪音数据。为后续续至波噪音的衰减操作提供了准确的数据基础，提高了地震数据处理的精度。
图15是本申请又一个实施例中的续至波衰减处理装置的结构示意图，如图15所示，在上述实施例的基础上，所述续至波衰减处理装置还包括：
反拉平处理模块151，用于在所述噪音衰减模块126获得噪音衰减地震数据后，将所述噪音衰减地震数据进行反拉平处理。
本申请提供的续至波衰减处理装置，在获得噪音衰减地震数据后，通过反拉平数据处理，可以保证进行续至波噪音衰减处理后的数据，保持原数据的特征，与衰减前的数据相比，仅仅只是缺少了续至波噪音数据。保证了数据的准确性，为后续地震资料的解释等操作，提供了准确的数据基础。
在上述实施例的基础上，所述噪音衰减模块还用于：
判断所述噪音衰减地震数据是否满足预设的衰减要求，若不满足，则调整所述噪音预测数据，根据调整后的噪音预测数据进行续至波噪音的衰减处理。
本申请提供的续至波衰减处理装置，可以通过调整噪音预测数据的频率，重复上述实施例中的续至波衰减操作，直至波衰减操作后的噪音衰减地震数据满足预设的衰减要求。确保续至波衰减操作的精度，保证续至波衰减操作后的地震数据中的续至波噪音数据尽可能的被衰减，以保证后续地震数据的处理和解释的准确性。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书实施例提供的上述续至波衰减处理方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在PC端实现、linux系统实现，或其他例如使用android、iOS系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种续至波衰减处理装置的另一个实施例中，图16是本申请提供的另一种续至波衰减处理装置实施例的模块结构示意图，如图16所示，本申请另一实施例提供的续至波衰减处理装置可以包括处理器161以及用于存储处理器可执行指令的存储器162，
处理器161和存储器162通过总线163完成相互间的通信；
所述处理器161用于调用所述存储器162中的程序指令，以执行上述各续至波衰减处理方法实施例所提供的方法，例如包括：根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；根据所述续至波数据拾取目标层位；根据所述目标层位和所述续至波数据，建立地层模型；根据所述地层模型，进行噪音预测处理，获得预测噪音数据；将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device，PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书一个或多个实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。第一，第二等词语用来表示名称，而并不表示任何特定的顺序。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本发明是参照根据本发明实施例的方法、装置(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储、石墨烯存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书一个或多个实施例的实施例而已，并不用于限制本本说明书一个或多个实施例。对于本领域技术人员来说，本说明书一个或多个实施例可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在权利要求范围之内。
Claims (11)
1.一种续至波衰减处理方法，其特征在于，所述方法包括：
根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；
根据所述续至波数据拾取目标层位；
根据所述目标层位和所述续至波数据，建立地层模型；
根据所述地层模型，进行噪音预测处理，获得预测噪音数据；
将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；
根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
2.如权利要求1所述的一种续至波衰减处理方法，其特征在于，所述根据所述目标层位和所述续至波数据，建立地层模型，包括：
将所述目标层位的数据进行拉平处理，获取所述目标层位的拉平层位数据；
给定空间时窗，根据所述拉平层位数据和所述续至波数据，获取所述空间时窗内的波形特征，建立所述地层模型，所述地层模型包括所述波形的参数信息。
3.如权利要求1所述的一种续至波衰减处理方法，其特征在于，所述根据所述地层模型，进行噪音预测处理，获得预测噪音数据，包括：
将所述地层模型中的数据由时间空间域转化为频率空间域；
根据所述地层模型，在频率域内利用最小二乘反演法，进行噪音预测滤波处理，获得所述预测噪音数据。
4.如权利要求2所述的一种续至波衰减处理方法，其特征在于，所述方法还包括：
获得噪音衰减地震数据后，将所述噪音衰减地震数据进行反拉平处理。
5.如权利要求1-4任一项所述的一种续至波衰减处理方法，其特征在于，所述方法还包括：
判断所述噪音衰减地震数据是否满足预设的衰减要求，若不满足，则调整所述噪音预测数据，根据调整后的噪音预测数据进行续至波噪音的衰减处理。
6.一种续至波衰减处理装置，其特征在于，包括：
续至波识别模块，用于根据测井数据、地质分层数据、地震资料，识别出待分析地震数据中的续至波数据；
目标层位拾取模块，用于根据所述续至波数据拾取目标层位；
地层模型建立模块，用于根据所述目标层位和所述续至波数据，建立地层模型；
噪音预测模块，用于根据所述地层模型，进行噪音预测处理，获得预测噪音数据；
噪音匹配模块，用于将所述预测噪音数据和所述待分析地震数据进行匹配，获得匹配噪音数据；
噪音衰减模块，用于根据所述待分析地震数据和所述匹配噪音数据进行续至波噪音衰减处理，获得噪音衰减地震数据。
7.如权利要求6所述的一种续至波衰减处理装置，其特征在于，所述地层模型建立模块包括：
数据拉平模块，用于将所述目标层位的数据进行拉平处理，获取所述目标层位的拉平层位数据；
数据处理模块，用于给定空间时窗，根据所述拉平层位数据和所述续至波数据，获取所述空间时窗内的波形特征，建立所述地层模型，所述地层模型包括所述波形的参数信息。
8.如权利要求6所述的一种续至波衰减处理装置，其特征在于，所述噪音预测模块包括：
频率转换模块，用于将所述地层模型中的数据由时间空间域转化为频率空间域；
反演预测模块，用于根据所述地层模型，在频率域内利用最小二乘反演法，进行噪音预测滤波处理，获得所述预测噪音数据。
9.如权利要求7所述的一种续至波衰减处理装置，其特征在于，所述续至波衰减处理装置还包括：
反拉平处理模块，用于在所述噪音衰减模块获得噪音衰减地震数据后，将所述噪音衰减地震数据进行反拉平处理。
10.如权利要求6-9任一项所述的一种续至波衰减处理装置，其特征在于，所述噪音衰减模块还用于：
判断所述噪音衰减地震数据是否满足预设的衰减要求，若不满足，则调整所述噪音预测数据，根据调整后的噪音预测数据进行续至波噪音的衰减处理。
11.一种续至波衰减处理装置，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现如权利要求1至5中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710990343.0A CN107942387B (zh) | 2017-10-23 | 2017-10-23 | 一种续至波衰减处理方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710990343.0A CN107942387B (zh) | 2017-10-23 | 2017-10-23 | 一种续至波衰减处理方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107942387A true CN107942387A (zh) | 2018-04-20 |
| CN107942387B CN107942387B (zh) | 2019-05-07 |
Family
ID=61935551
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710990343.0A Active CN107942387B (zh) | 2017-10-23 | 2017-10-23 | 一种续至波衰减处理方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107942387B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111624659A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种地震数据的时变带通滤波方法及装置 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1338640A (zh) * | 2001-06-29 | 2002-03-06 | 中国石油天然气股份有限公司 | 一种模型约束动静校正方法 |
| CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
| CN101627322A (zh) * | 2007-03-09 | 2010-01-13 | 格库技术有限公司 | 可控震源采集方法 |
| CN103630937A (zh) * | 2013-12-09 | 2014-03-12 | 中国石油天然气集团公司 | 对地震记录进行横波静校正处理的方法和装置 |
| CN104459789A (zh) * | 2014-12-08 | 2015-03-25 | 翟明岳 | 一种地震波初至拾取的方法 |
| CN104730579A (zh) * | 2013-12-18 | 2015-06-24 | 中国石油化工股份有限公司 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
| CN104808245A (zh) * | 2015-05-19 | 2015-07-29 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 道集优化处理方法及其装置 |
| CN105068119A (zh) * | 2015-07-07 | 2015-11-18 | 中国石油天然气集团公司 | 低频地震数据中的面波的衰减方法及其装置 |
| RU2601232C2 (ru) * | 2012-08-20 | 2016-10-27 | Лэндмарк Графикс Корпорейшн | Способы и системы для включения мест выделенных точек псевдоповерхности в модели скорости сейсмических волн |
| CN106255902A (zh) * | 2014-05-08 | 2016-12-21 | 雪佛龙美国公司 | 使用地震数据分析地质特征的系统和方法 |
-
2017
- 2017-10-23 CN CN201710990343.0A patent/CN107942387B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1338640A (zh) * | 2001-06-29 | 2002-03-06 | 中国石油天然气股份有限公司 | 一种模型约束动静校正方法 |
| CN101627322A (zh) * | 2007-03-09 | 2010-01-13 | 格库技术有限公司 | 可控震源采集方法 |
| CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
| RU2601232C2 (ru) * | 2012-08-20 | 2016-10-27 | Лэндмарк Графикс Корпорейшн | Способы и системы для включения мест выделенных точек псевдоповерхности в модели скорости сейсмических волн |
| CN103630937A (zh) * | 2013-12-09 | 2014-03-12 | 中国石油天然气集团公司 | 对地震记录进行横波静校正处理的方法和装置 |
| CN104730579A (zh) * | 2013-12-18 | 2015-06-24 | 中国石油化工股份有限公司 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
| CN106255902A (zh) * | 2014-05-08 | 2016-12-21 | 雪佛龙美国公司 | 使用地震数据分析地质特征的系统和方法 |
| CN104459789A (zh) * | 2014-12-08 | 2015-03-25 | 翟明岳 | 一种地震波初至拾取的方法 |
| CN104808245A (zh) * | 2015-05-19 | 2015-07-29 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 道集优化处理方法及其装置 |
| CN105068119A (zh) * | 2015-07-07 | 2015-11-18 | 中国石油天然气集团公司 | 低频地震数据中的面波的衰减方法及其装置 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111624659A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种地震数据的时变带通滤波方法及装置 |
| CN111624659B (zh) * | 2020-06-05 | 2022-07-01 | 中油奥博(成都)科技有限公司 | 一种地震数据的时变带通滤波方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107942387B (zh) | 2019-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107976713B (zh) | 一种高维地震数据输入下去除沉积背景的方法及装置 | |
| CN104950334B (zh) | 一种预测储层分布的方法及装置 | |
| CN108181649B (zh) | 一种介质分界面振幅补偿方法及装置 | |
| CN109100796A (zh) | 一种天然气水合物地震数据处理方法及装置 | |
| CN107728210B (zh) | 多仪器采集数据中缺失道的确定方法和装置 | |
| CN107783185B (zh) | 一种层析静校正的处理方法及装置 | |
| CN106338763B (zh) | 一种确定地震剖面显示数据的方法及装置 | |
| CN108196304B (zh) | 一种多次波压制方法及装置 | |
| CN106353818B (zh) | 一种地震数据时差校正的方法及系统 | |
| CN108089238A (zh) | 一种利用综合地球物理资料确定深层裂谷的方法及装置 | |
| CN108445538A (zh) | 基于反射地震资料建立深度域层q模型的方法和系统 | |
| CN108897036A (zh) | 一种地震数据处理方法及装置 | |
| CN107797154A (zh) | 一种电成像测井图像刻度方法及装置 | |
| CN103675904B (zh) | 一种井震匹配解释性目标处理方法及装置 | |
| CN107179546B (zh) | 一种确定速度场的方法及装置 | |
| CN107656315B (zh) | 储层流体因子的确定方法和装置 | |
| CN107144878A (zh) | 一种断层识别方法及装置 | |
| CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
| CN107942387B (zh) | 一种续至波衰减处理方法及装置 | |
| CN108693557A (zh) | 一种地震数据处理方法及装置 | |
| CN108254785A (zh) | 储层确定方法和装置 | |
| CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
| CN108710149A (zh) | 一种薄储层的预测方法及装置 | |
| CN109100797A (zh) | 一种确定致密储层含气性的方法及装置 | |
| CN107991716A (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 |