CN108828668B - 一种叠前时间偏移数据处理方法及装置 - Google Patents
一种叠前时间偏移数据处理方法及装置 Download PDFInfo
- Publication number
- CN108828668B CN108828668B CN201810256518.XA CN201810256518A CN108828668B CN 108828668 B CN108828668 B CN 108828668B CN 201810256518 A CN201810256518 A CN 201810256518A CN 108828668 B CN108828668 B CN 108828668B
- Authority
- CN
- China
- Prior art keywords
- ovt
- data
- processing
- information
- seismic data
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 242
- 238000013508 migration Methods 0.000 title claims abstract description 242
- 238000003672 processing method Methods 0.000 title claims description 15
- 238000012545 processing Methods 0.000 claims abstract description 418
- 238000012937 correction Methods 0.000 claims abstract description 154
- 230000003068 static effect Effects 0.000 claims abstract description 138
- 238000000034 method Methods 0.000 claims abstract description 120
- 238000004458 analytical method Methods 0.000 claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims description 61
- 230000015654 memory Effects 0.000 claims description 30
- 239000013598 vector Substances 0.000 claims description 12
- 230000004044 response Effects 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 25
- 230000008569 process Effects 0.000 description 23
- 238000003384 imaging method Methods 0.000 description 19
- 238000005516 engineering process Methods 0.000 description 17
- 238000003860 storage Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 12
- 230000006872 improvement Effects 0.000 description 10
- 238000004590 computer program Methods 0.000 description 8
- 230000010354 integration Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000005284 excitation Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 229910021389 graphene Inorganic materials 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002123 temporal effect Effects 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
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
本申请提供一种叠前时间偏移数据处理方法及装置。所述方法包括：根据采集到的宽方位地震数据，获得OVT域地震数据；将所述OVT域地震数据进行静校正处理，获得静校正处理数据；将所述静校正处理数据进行速度分析，获得均方根速度场；在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT集地震数据；根据所述OVT集地震数据和所述均方根速度场，进行叠前时间偏移处理。利用本申请中各实施例，提高了叠前时间偏移的处理效率，节省了叠前时间偏移数据处理的时间。
Description
技术领域
本申请属于石油物理勘探技术领域，尤其涉及一种叠前时间偏移数据处理方法及装置。
背景技术
地震数据处理的重要目的之一是通过对野外地震采集的数据做各种处理后提高反射波数据的信噪比、分辨率，以便于地震数据的后期解释。叠前时间偏移是复杂构造成像最有效的方法之一，适用于大倾角，横纵向速度变化较大的情况。随着地震勘探的目标复杂程度日益增加，对地震成像的精度要求也越来越高，使得叠前时间偏移已经成为地震资料常规处理中的一个重要环节。
现有技术中，叠前时间偏移方法中常用的是克希霍夫积分法偏移，该偏移方法主要是应用波动方程的克希霍夫积分解来解决反射层的偏移问题。现有技术中，叠前时间偏移数据处理方法通常都是以单个数据道为单位进行数据处理，随着地震勘探技术的不断发展，地震采集道数已从原来的几百道发展为上千道、上万道。随着地震数据采集量的日益增大，常规的叠前时间偏移方法在针对海量的数据资料时，运算时间过长，已经难以满足当前应用的需求。因此，业内亟需一种能够提高叠前时间偏移数据处理速度的实施方案。
发明内容
本申请目的在于提供一种叠前时间偏移数据处理方法及装置，提高了地震数据叠前时间偏移处理的处理速度，节省了数据处理时间。
一方面本申请提供了一种叠前时间偏移数据处理方法，包括：
根据采集到的宽方位地震数据，获得偏移距矢量片OVT域地震数据；
将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
将所述静校正处理数据进行速度分析，获得均方根速度场；
在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
进一步地，所述方法的另一个实施例中，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
对所述述静校正处理数据进行OVT计算，获取所述OVT信息，将所述OVT信息添加到所述述静校正处理后的OVT域地震数据中；
利用五维叠前面元数据规则化处理方法，对添加所述OVT信息的静校正处理数据进行规则化处理，获得所述OVT域偏移地震数据。
进一步地，所述方法的另一个实施例中，所述对添加所述OVT信息的静校正处理数据进行规则化处理，获得所述OVT域偏移地震数据，包括：
对添加所述OVT信息的静校正处理数据进行规则化处理后，对所述规则化处理后的地震数据进行OVT计算，将获得的OVT信息添加到所述规则化处理后的地震数据中，获得所述OVT域偏移地震数据。
进一步地，所述方法的另一个实施例中，所述方法还包括：
将所述静校正处理数据进行噪音处理，获得噪音处理数据；
相应地，所述将所述静校正处理数据进行速度分析，获得均方根速度场，包括：
将所述噪音处理数据进行速度分析，获得所述均方根速度场；
相应地，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
在所述噪音处理数据中添加所述OVT信息，对添加所述OVT信息的噪音处理数据进行规则化处理，获得OVT域偏移地震数据。
进一步地，所述方法的另一个实施例中，所述方法还包括：
将所述静校正处理数据进行振幅补偿，并将所述振幅补偿后的地震数据进行反褶积处理，获得补偿处理数据；
相应地，所述将所述静校正处理数据进行速度分析，获得均方根速度场，包括：
将所述补偿处理数据进行速度分析，获得所述均方根速度场；
相应地，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
在所述补偿处理数据中添加所述OVT信息，对添加所述OVT信息的补偿处理数据进行规则化处理，获得OVT域偏移地震数据。
另一方面，本申请提供了一种叠前时间偏移数据处理装置，包括：
宽方位数据处理模块，用于根据采集到的宽方位地震数据，获得偏移距矢量片OVT域地震数据；
静校正模块，用于将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
速度分析模块，用于将所述静校正处理数据进行速度分析，获得均方根速度场；
规则化处理模块，用于在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
叠前时间偏移模块，用于根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
进一步地，所述装置的另一个实施例中，所述规则化处理模块，包括：
OVT计算单元，用于对所述述静校正处理数据进行OVT计算，获取所述OVT信息，将所述OVT信息添加到所述静校正处理数据中；
规则化处理单元，用于利用五维叠前面元数据规则化处理方法，对添加所述OVT信息的静校正处理数据内进行规则化处理，获得所述OVT域偏移地震数据。
进一步地，所述装置的另一个实施例中，所述规则化处理模块还包括：
OVT二次计算，用于对添加所述OVT信息的静校正处理数据进行规则化处理后，对所述规则化处理后的地震数据进行OVT计算，将获得的OVT信息添加到所述规则化处理后的地震数据中，获得所述OVT域偏移地震数据。
又一方面，本申请还提供了一种叠前时间偏移数据处理装置，包括：包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
根据采集到的宽方位地震数据，获得偏移距矢量片OVT域地震数据；
将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
将所述静校正处理数据进行速度分析，获得均方根速度场；
在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
再一方面，本申请还提供了一种叠前时间偏移数据处理系统，包括：处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述叠前时间偏移数据处理方法。
本申请提供的叠前时间偏移数据处理方法及装置，利用宽方位采集技术获得的宽方位地震数据，获得OVT域地震数据，并将OVT信息应用到地震数据的规则化处理过程中，获得OVT域偏移地震数据。将OVT域偏移地震数据应用到地震数据的叠前时间偏移处理中，将现有技术中叠前时间偏移以单道处理模式改为OVT片处理模式，即将属于同一个OVT片的道数据分选出来，参与偏移计算。一个OVT片可以包括多个道集数据，可以大大提高叠前时间偏移计算的运算效率，提升了叠前时间偏移的数据处理速度，节省了运行时间。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的叠前时间偏移数据处理方法一个实施例的方法流程示意图；
图2是本申请又一个实施例中叠前时间偏移处理的流程示意图；
图3是本申请实施例中五维叠前面元数据规则化处理的流程示意图；
图4是本申请又一个实施例中叠前时间偏移的流程示意图；
图5是现有技术中叠前时间偏移方法生成的某工区的OVT域叠前时间偏移示意图；
图6是本申请一个实施例中叠前时间偏移方法生成的某工区的OVT域叠前时间偏移示意图；
图7是现有技术中叠前时间偏移方法生成的某工区的OVT域共偏移距道集；
图8是本申请一个实施例中叠前时间偏移方法生成某工区的OVT域共偏移距道集；
图9是本申请一个实施例中OVT域叠前时间偏移的均方根速度场的示意图；
图10是本申请实施例提供的叠前时间偏移数据处理装置的模块结构示意图；
图11是本申请一个实施例中规则化处理模块的结构示意图；
图12是本申请又一个实施例中规则化处理模块的结构示意图；
图13是本申请提供的一种叠前时间偏移数据处理系统实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
偏移有两种方式：叠加之前偏移和叠加之后偏移，前者就叫叠前，后者叫叠后。叠前时间偏移是复杂构造成像最有效的方法之一，能适应纵横向速度变化较大的情况，适用于大倾角的偏移成像。叠前时间偏移方法主要分为两类，即用于准确构造成像的叠前时间偏移和振幅保持叠前时间偏移，使用比较多的是采用克希霍夫积分法进行叠前时间偏移。
近些年来，通过宽方位高密度地震勘探方式获得了海量的高品质宽方位地震信息，通过对这些宽方位地震勘探获取的高品质数据进行处理分析，可以得到信噪比较高的叠前数据以及合理的偏移速度模型，为后续的地震资料处理提供了信噪比较高的规整数据。
本申请实施例中，利用OVT(Offset Vector Tile，偏移距矢量片)域宽方位地震采集数据，将OVT计算应用到数据规则化处理中，将数据划分为OVT域偏移地震数据(可以表示属于同一个OVT片的道数据)，并将OVT域规则化处理后的地震数据进行叠前时间偏移。改变了现有技术中叠前时间偏移对数据的处理模式，由单道处理模式改为OVT片处理模式，即将属于同一个OVT片的道数据分选出来，参与偏移计算，可以大大提高偏移计算的运算效率，节省运行时间。
具体地，图1是本申请提供的叠前时间偏移数据处理方法一个实施例的方法流程示意图，本申请提供的叠前时间偏移数据处理方法包括：
S1、根据采集到的宽方位地震数据，获得OVT域地震数据。
宽方位地震采集技术可以包括宽方位角三维地震采集方法，利用宽方位采集技术可以采集到宽方位地震数据。可以将野外采集得到的宽方位地震数据按照OVT域进行分选，得到OVT域宽方位高密度地震数据即OVT域地震数据。野外采集得到的宽方位地震数据所对应的观测系统参数可以包括：炮线距、接收线距、面元、最小覆盖次数以及偏移距等。对野外采集得到的宽方位地震数据按照OVT域进行分选可以包括：获取炮线和检波线相互垂直的所有地震道集的集合，即十字排列子集；然后对所述十字排列子集中每一个十字排列按炮线距和检波线距进行等距离划分，得到许多小矩形，每个小矩形可以称为OVT向量片，从所有的OVT向量片中获取偏移距和方位角大致相同的OVT向量片，所获取的OVT向量片可以称为OVT道集数据体即OVT域地震数据。
获得OVT域地震数据后，可以将OVT域地震数据加载到预先建立的工区测线下。通过读取SPS格式文件，应用野外采集提供的炮点文件信息、检波点文件信息，关系文件信息，进行SPS观测系统的定义，可以将野外的观测信息加载到处理系统中，使得观测信息与实际数据一一对应。即可以利用野外采集的SPS格式文件，将野外采集数据的道头信息应用到OVT域地震数据中，使得地震数据的道头信息更加丰富，便于后续的叠前时间偏移处理。
S2、将所述OVT域地震数据进行静校正处理，获得静校正处理数据。
在陆上地震数据处理过程中，通常要将地震数据校正到一个统一的基准面，这个基准面一般为水平面。地震勘探解释的理论都假定激发点与接收点是在一个水平面上，并且地层速度是均匀的。但实际上地面常常不平坦，各个激发点深度也可能不同，低速带中的波速与地层中的波速又相差悬殊，所以必将影响实测的时距曲线形状。对原始地震数据要进行地形校正、激发深度校正、低速带校正等，这些校正对同一观测点的不同地震界面都是不变的，可以称为静校正，通过对地震数据进行静校正处理可以消除上述影响。可以认为静校正量＝校正后时间-原始时间，可以通过对OVT域地震数据进行静校正处理，将计算得到的静校正量分离为低频分量(CMP(Common Middle Point，共中心点道集)参考面到浮动基准面的校正量)与高频分量(炮点和接收点到CMP参考面的校正量)。在静校正处理中可以应用高频分量，将高频分量数据校正到CMP参考面，实现静校正处理。静校正处理可以将二维的地震数据转化为三维的地震数据，为后续速度分析、叠前时间偏移提供数据基础。通过对OVT域地震数据进行静校正处理可以改进CMP同向叠加的叠加效果，提高地震剖面的信噪比、叠加剖面的垂向分辨率，从进一步而改进叠加速度分析的质量。
S3、将所述静校正处理数据进行速度分析，获得均方根速度场。
速度分析可以包括速度谱计算以及速度场的建立，可以对静校正处理后的地震数据即静校正处理数据，进行速度谱计算，获得速度谱，然后对产生的速度谱进行速度解释建速度场，得到精确的均方根速度场。速度场可以表示每一时刻、每一点上的速度矢量组成的物理场，均方根速度场可以表示把水平层状介质的反射波时距曲线近似的当作双曲线求出的波束。速度偏移是叠前时间偏移的关键因素，获得的均方根速度场可以用于后续的叠前时间偏移处理中。
S4、在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据。
获得静校正处理数据后，可以通过对静校正处理数据进行OVT计算，获得地震数据中的OVT信息，OVT信息可以包括面元信息如：偏移距和方位角信息。OVT计算主要是可以包括对三维地震数据进行OVT分组计算，适应笛卡尔坐标和极坐标两种坐标系统。在笛卡尔坐标系下，可以按照炮检距OffsetX、OffsetY进行分组；在极坐标系下，可以按照炮检距、方位角进行分组。通过OVT计算获得地震数据中的OVT信息，可以将获得的OVT信息添加到静校正处理数据中的道头信息中，对添加OVT信息的静校正处理数据进行规则化处理。数据规则化是地震数据处理中重要的处理技术，它可以改善面元属性、炮检距分组、提高叠加数据的信噪比和偏移成像质量等。例如：常用的地震数据规则化处理技术包括三维地震数据规则化处理。通过对地震数据进行OVT域的规则化处理，可以将地震数据分选成OVT集地震数据，将属于同一个OVT片的道数据分选出来作为一个OVT域地震数据，为后续叠前时间偏移提供数据基础。
在本申请一个实施例中，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，可以包括：
对所述述静校正处理数据进行OVT计算，获取所述OVT信息，将所述OVT信息添加到所述述静校正处理后的OVT域地震数据中；
利用五维叠前面元数据规则化处理方法，对添加所述OVT信息的静校正处理数据进行规则化处理，获得所述OVT域偏移地震数据。
具体地，本申请一个实施例中，将OVT计算应用到传统的五维叠前面元数据规则化的处理流程中，用以保护来自不同方位的地下反射信息。OVT计算的主要内容可以参考上述实施例的介绍，通过OVT计算，可以直接修改静校正处理数据中的道头信息，在道头中增加OVT信息，OVT信息的具体内容也可以参考上述实施例的介绍，此处不赘述。五维叠前面元数据规则化处理方法可以通过读取OVT计算后得到的OVT信息，利用非均匀傅立叶重构技术，在不同的坐标系统下同时进行四个空间方向的规则化处理。五叠前面元数据规则化处理技术能综合利用三维地震数据的“纵向、横向、时间、偏移距、方位角”五个维度的信息，通过将OVT计算应用到五维叠前面元数据规则化处理过程中，使得空间方向不均匀采样得到规则化重建，从而将炮检距、覆盖次数等属性的不均匀性得到改善，同时将缺失的地震道重新创建。规则化后的地震数据可以实现将地震数据分选成OVT域偏移地震数据，一方面提高了OVT面元的密度，同时可以改善后续叠前时间偏移的效果。
S5、根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
对宽方位地震数据进行速度分析和OVT域内的规则化处理后，可以将规则化处理后的OVT域偏移地震数据以及速度分析后获得的均方根速度场作为输入数据，进行叠前时间偏移处理。例如：可以使用积分法进行叠前时间偏移处理。图2是本申请又一个实施例中叠前时间偏移处理的流程示意图，如图2所示，可以将规则化处理后的地震数据和均方根速度场作为输入数据，进行叠前时间偏移处理，判断经过叠前时间偏移处理后的地震数据的CRP道集(共反射道集)是否被拉平，若是，则输出CRP道集，若否，则对叠前时间偏移后的地震数据进行反动校正以及速度分析，获得新的均方根速度场，继续进行叠前时间偏移处理。最后将叠前时间偏移处理的道集数据继续进行切除和叠加，获得叠前时间偏移剖面，完成地震数据的叠前时间偏移处理。切除指的是叠前时间偏移生成的CRP道集如果远偏移距有拉伸畸变，即不平整，我们需要用切除软件将CRP道集的远偏移距进行切除，切除后的数据用于叠加。叠加指的是将切除后的CRP道集数据叠加起来，以便达到提高信噪比，压制干扰波的目的。
本申请实施例提供的叠前时间偏移数据处理方法，利用宽方位采集技术获得的宽方位地震数据，获得OVT域地震数据，并将OVT信息应用到地震数据的规则化处理过程中，获得OVT域偏移地震数据。将OVT域偏移地震数据应用到地震数据的叠前时间偏移处理中，将现有技术中叠前时间偏移以单道处理模式改为OVT片处理模式，即将属于同一个OVT片的道数据分选出来，参与偏移计算。一个OVT片可以包括多个道集数据，可以大大提高叠前时间偏移计算的运算效率，提升了叠前时间偏移的数据处理速度，节省了运行时间。
在上述实施例的基础上，所述对添加所述OVT信息的静校正处理数据进行规则化处理，获得所述OVT域偏移地震数据，包括：
对添加所述OVT信息的静校正处理数据进行规则化处理后，对所述规则化处理后的地震数据进行OVT计算，将获得的OVT信息添加到所述规则化处理后的地震数据中，获得所述OVT域偏移地震数据。
图3是本申请实施例中五维叠前面元数据规则化处理的流程示意图，如图3所示，本申请一个实施例中，在对静校正处理后的地震数据进行规则化处理时，可以先对静校正处理数据进行动校正，以消除地震波到达各检波点的正常时差。动校正处理后，可以对地震数据进行OVT计算，将获得的OVT信息添加到地震数据的道头信息中。再通过五维叠前面元数据规则化处理模块读取地震数据中的OVT信息，对地震数据进行五维叠前面元数据规则化处理。五维叠前面元数据规则化处理后，重新计算道头信息，并利用OVT计算，将OVT信息添加到五维叠前面元数据规则化处理后的地震数据的道头信息中，获得OVT域偏移地震数据，确保地震数据规则化处理结果更加准确。将规则化处理后获得的OVT域偏移地震数据进行反动校正处理，规则化处理后的地震数据即OVT域偏移地震数据可以用于叠前时间偏移处理。
在上述实施例的基础上，所述方法还可以包括：
将所述静校正处理数据进行噪音处理，获得噪音处理数据；
相应地，所述将所述静校正处理数据进行速度分析，获得均方根速度场，包括：
将所述噪音处理数据进行速度分析，获得所述均方根速度场；
相应地，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
在所述噪音处理数据中添加所述OVT信息，对添加所述OVT信息的噪音处理数据进行规则化处理，获得OVT域偏移地震数据。
具体地，根据实际应用的需要，在对OVT域地震数据进行静校正处理后，可以对静校正处理后获得的静校正处理数据进行噪音处理，去除地震数据中的噪音数据，提高地震数据的信噪比，为后续叠前时间偏移提供准确的数据基础。相应地，噪音处理后，可以对噪音处理获得的噪音处理数据进行速度分析，获得均方根速度场。还可以对噪音处理数据进行规则化处理，即可以对噪音处理数据进行OVT计算，将OVT信息添加到噪音处理数据中，并对添加了OVT信息的噪音处理数据进行规则化处理，获得OVT域偏移地震数据。
其中，噪音处理的具体方法可以根据实际需要进行选择，如：滤波处理、小波变换等，本申请实施例不作具体限定。
在上述实施例的基础上，所述方法还包括：
将所述静校正处理数据进行振幅补偿，并将所述振幅补偿后的地震数据进行反褶积处理，获得补偿处理数据；
相应地，所述将所述静校正处理数据进行速度分析，获得均方根速度场，包括：
将所述补偿处理数据进行速度分析，获得所述均方根速度场；
相应地，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
在所述补偿处理数据中添加所述OVT信息，对添加所述OVT信息的补偿处理数据进行规则化处理，获得OVT域偏移地震数据。
具体地，根据实际应用的需要，在对OVT域地震数据进行静校正处理后，可以对静校正处理后获得的静校正处理数据进行振幅补偿处理，通过振幅补偿可以恢复地震数据的振幅的能量，使得振幅值基本达到一致水平，保证数据满足后续处理的需求。可以通过地表一致性振幅补偿方法消除地震记录在道间的振幅差异，通过调用地表一致性振幅补偿模块就可以实现振幅补偿。振幅补偿处理后，可以对振幅补偿后的地震数据进行反褶积处理，消除由于激发、接收等因素引起的地震记录在道间振幅、子波的差异。相应地，可以对振幅补偿以及反褶积处理获得的补偿处理数据进行速度分析，获得均方根速度场。还可以对补偿处理数据进行规则化处理，即可以对补偿处理数据进行OVT计算，将OVT信息添加到补偿处理数据中，并对添加了OVT信息的补偿处理数据进行规则化处理，获得OVT域偏移地震数据。
实际应用时，可以根据需要选择是否需要对地震数据进行噪音处理、振幅补偿处理、反褶积处理，可以选择其中一个或多个数据处理过程同时使用，本申请实施例不作具体限定。具体地，图4是本申请又一个实施例中叠前时间偏移的流程示意图，如图4所示，本申请一个示例中叠前时间偏移还可以包括以下过程：
1)将野外采集得到的宽方位地震数据按照OVT域进行分选，得到OVT域宽方位高密度地震数据(即OVT域地震数据)，将OVT域地震数据加载到工区测线下。通过读取野外采集到的地震数据的SPS格式文件，应用野外采集提供的炮点文件信息、检波点文件信息，关系文件信息，进行SPS观测系统的定义，将野外的观测信息加载到处理系统中，使得观测信息与实际数据一一对应。
2)根据实际需要，进行OVT域宽方位高密度采集地震数据的静校正处理，通常静校正量＝校正后时间-原始时间。通过应用静校正计算公式，计算得到的静校正量分离为低频分量(CMP参考面到浮动基准面的校正量)与高频分量(炮点和接收点到CMP参考面的校正量)，一般常规处理中可以应用高频分量，将数据校正到CMP参考面，用此方法改进CMP同向叠加的叠加效果，提高地震剖面的信噪比、叠加剖面的垂向分辨率，进一步改进叠加速度分析的质量。
3)根据实际需要，将静校正处理后的数据进行OVT域宽方位高密度采集叠前地震数据的噪音处理，确定合理的叠前去噪流程，从而达到提高信噪比的目的。
4)根据实际需要，进行OVT域宽方位高密度采集地震数据的振幅补偿，可以对噪音处理后的地震数据进行振幅补偿处理，恢复振幅的能量，使得振幅值基本达到一致水平，保证数据满足后续处理的需求。可以应用地表一致性振幅补偿方法消除地震记录在道间的振幅差异，通过调用地表一致性振幅补偿模块可以实现该步骤。
5)根据实际需要，进行OVT域宽方位高密度采集地震数据的反褶积处理，可以对振幅补偿后的地震数据进行反褶积处理，消除由于激发、接收等因素引起的地震记录在道间振幅、子波的差异。
6)根据实际需要，进行OVT域宽方位高密度采集地震数据的速度分析，即对上述静校正处理、噪音处理、振幅补偿处理、反褶积处理后的地震数据进行速度分析。可以应用速度分析的批量模型进行速度谱计算，然后用交互速度分析模块对产生的速度谱进行速度解释建速度场，得到精确的均方根速度场用于叠前时间偏移。
7)根据实际需要，将数据完成步骤1)至步骤6)的处理后，对地震数据进行OVT域宽方位高密度采集地震数据的五维叠前面元数据规则化处理，常规的规则化处理流程仅是应用五维叠前面元数据规则化模块对数据进行处理，得到规则化后的地震道。在本申请实施例中将OVT计算模块应用到传统的五维叠前面元数据规则化的处理流程中，用以保护来自不同方位的地下反射信息。OVT计算主要是对三维地震数据进行OVT分组计算，适应笛卡尔坐标和极坐标两种坐标系统。在笛卡尔坐标系下，按照炮检距OffsetX、OffsetY进行分组；在极坐标系下，按照炮检距、方位角进行分组。通过OVT计算，可以直接修改数据的道头信息，在道头中增加OVT信息。五维叠前面元数据规则化模块通过读取OVT模块计算后得到的OVT信息，利用非均匀傅立叶重构技术，在不同的坐标系统下同时进行四个空间方向的规则化处理，使得空间方向不均匀采样得到规则化重建，从而将炮检距、覆盖次数等属性的不均匀性得到改善，同时将缺失的地震道重新创建，规则化后的地震数据一方面提高了OVT面元的密度，同时改善了叠前时间偏移的效果。
需要说明的是，图4中叠前时间偏移处理过程中，叠前数据五维规则化处理即本申请实施例中的五维叠前面元数据规则化处理。图4中先进性五维叠前面元数据规则化处理再进行速度分析处理，实际应用中，也可以先进性速度分析处理再进行五维叠前面元数据规则化处理，还可以在对地震数据进行静校正处理、噪音处理、振幅补偿处理、反褶积处理后，同时进行五维叠前面元数据规则化处理和速度分析处理，本申请实施例不作具体限定。图中，速度分析后的动校正、叠加、建立均方根速度场、剩余静校正等迭代处理过程可以作为速度分析的过程。
8)根据实际需要，改变原有克希霍夫积分法时间偏移对数据的处理模式，首先将OVT计算模块应用到五维叠前面元数据规则化的处理流程中，将和OVT相关的道头信息加载到地震数据中。应用五维叠前面元数据规则化模块对数据在OVT域进行规则化，将数据分选成OVT域偏移地震数据，目的就是在频率域针对规则化后的数据由单道数据处理模式改为OVT片偏移。因为经过前面步骤1)-7)的处理后，已经将数据分选为OVT域数据，即将数据按照不同的OVT片进行分选，分选后的数据用于克希霍夫积分法叠前时间偏移，这样每个OVT片的地震道一次性偏移，会大大节省OVT域叠前时间偏移的运算时间。
本申请在实施过程中，用50个节点做过对比试验，同样的输入数据，原有的克希霍夫积分法叠前时间偏移全部运行结束需要90小时，而改进后的克希霍夫积分法叠前时间偏移全部运行结束仅仅需要要3小时，大大节省了运算时间，极大的提高了模块的运行速度。
图5是现有技术中叠前时间偏移方法生成的某工区的OVT域叠前时间偏移示意图，图6是本申请一个实施例中叠前时间偏移方法生成的某工区的OVT域叠前时间偏移示意图，如图5-图6所示，图5和图6的成像剖面结果基本完全一致，图5和图6的成像结果说明在同一数据，应用常规叠前时间偏移方法生成的偏移成像结果作业运行时间需要90小时，而本申请实施例中的叠前时间偏移方法生成的偏移成像结果作业运行时间仅需要3小时。因此本申请实施例提供的叠前时间偏移方法，在不改变成像结果的前提下，极大提高了叠前时间偏移模块的运行效率，从而也大大节省了作业的运行时间。
图7是现有技术中叠前时间偏移方法生成的某工区的OVT域共偏移距道集，图8是本申请一个实施例中叠前时间偏移方法生成某工区的OVT域共偏移距道集，如图7-图8所示，图中横坐标可以表示OVT域共偏移距道集的共中心点数CMP和道数Trace，纵坐标为共偏移距道集的深度值。图7和图8两图的OVT域共偏移距道集结果基本完全一致，应用图7和图8说明在OVT域共偏移距道集结果一致的前提下，本申请实施例提供的叠前时间偏移方法，在不改变成像结果的前提下，极大提高了叠前时间偏移模块的运行效率，从而也大大节省了作业的运行时间。
图9是本申请一个实施例中OVT域叠前时间偏移的均方根速度场的示意图，如图9所示，本申请实施例提供的OVT域叠前时间偏移方法，因叠前时间偏移模块运行效率的提升，导致均方根速度场的更新迭代速度也加快，因此从整体上说本申请实施例大大提高了整个叠前时间偏移流程的迭代更新速度，不仅大大缩减了整个生产项目的工期时间，节省了成本，也提高了经济效益。
本申请实施例，改进了传统的OVT域宽方位高密度采集地震数据叠前时间偏移的常规处理模式，通常的处理系统提供的叠前时间偏移方法都是基于克希霍夫积分法研发的，如果地震数据量是海量数据体，传统的单道处理模式会导致叠前时间偏移模块的计算时间过长，有时还会延迟生产项目的工期，给应用人员带来了极大的不方便。本申请实施例在进行OVT域叠前时间偏移处理时，改进了传统的单道时间偏移模式，采用OVT片偏移的多道处理方法进行OVT域叠前时间偏移，大大提高了OVT域叠前时间偏移的运算效率。本申请实施例通过对OVT域宽方位资料处理技术进行分析，解决了如何根据地表特征和数据资料特点，设计处理流程逐步提高了海量宽方位地震数据的品质，同时在保证时间域偏移成像质量理想的前提下，通过一种新的叠前时间偏移方法大大缩短了OVT域叠前时间偏移模块的计算时间，优化了OVT域叠前时间偏移模块的计算模式，扩大了OVT域叠前时间偏移模块的使用范围。
本申请提供的叠前时间偏移数据处理方法及装置是一种频率波数域偏移，从一个空间算子开始，得到一个固定相位的显式解。该算子包括空间和时间的相移，即：
其中：U为下行波场，ω是角频率，v是速度，z在时间偏移中对应的是时间，
为水平方向波数在方位角上的投影，且
而
上式中kx，ky是X和Y两个方向的波数。
偏移是一种成像的过程，成像结果以波场的值体现，在进行OVT域叠前时间偏移时，可以采用上述公式计算波场值U。
上述算子中的幂指数，可以分为三个部分，即三个相位的和。第一个部分的相位表示将能量沿着炮检方向进行移动，移动的量为一个倾角时差校正的距离b。第二个部分的相位等价于一个倾角时差校正，其中A是叠前偏移算子响应的椭圆的长半轴，B是共中心点与b的射线路径长度。第三个部分的相位等价于零炮检距的相位移偏移。
综上所述，对于每一个给定的参数集(z,v,ω,kx,ky)，都能够得到：
上式中：k＝2ω/v，h为半炮检距，A可以表示叠前偏移算子响应的椭圆的长半轴，B可以表示共中心点与b的射线路径长度，ka可以表示水平方向波数在方位角上的投影。
利用A、B、ka等参数，结合上述公式可以计算出波场U，即获得OVT域叠前时间偏移的波场值，进一步可以获得OVT域叠前时间偏移的成像结果，完成了偏移过程，实现了叠前时间偏移成像。
虽然上述算子是在常速条件下得到的，但是仍然能够应用于基于变速度的均方根速度模型的时间偏移，继而获得较高的计算效率。
本申请实施例提供的叠前时间偏移数据处理方法，提供的一种针对OVT域宽方位高密度采集数据的快速时间偏移成像技术，将OVT计算应用到数据规则化处理以及叠前时间偏移的过程中，与原有时间偏移技术在成像精度一致或更理想的情况下，具有计算量小，效率高，成本低的优点。
基于上述所述的叠前时间偏移数据处理方法，本说明书一个或多个实施例还提供一种叠前时间偏移数据处理装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
具体地，图10是本申请实施例提供的叠前时间偏移数据处理装置的模块结构示意图，如图10所示，本申请中提供的叠前时间偏移数据处理装置包括：宽方位数据处理模块101，静校正模块102，速度分析模块103，规则化处理模块104、叠前时间偏移模块105。
宽方位数据处理模块101，可以用于根据采集到的宽方位地震数据，获得OVT域地震数据；
静校正模块102，可以用于将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
速度分析模块103，可以用于将所述静校正处理数据进行速度分析，获得均方根速度场；
规则化处理模块104，可以用于在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
叠前时间偏移模块105，可以用于根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
本申请提供的叠前时间偏移数据处理装置，利用宽方位采集技术获得的宽方位地震数据，获得OVT域地震数据，并将OVT信息应用到地震数据的规则化处理过程中，获得OVT域偏移地震数据。将OVT域偏移地震数据应用到地震数据的叠前时间偏移处理中，将现有技术中叠前时间偏移以单道处理模式改为OVT片处理模式，即将属于同一个OVT片的道数据分选出来，参与偏移计算。一个OVT片可以包括多个道集数据，可以大大提高叠前时间偏移计算的运算效率，提升了叠前时间偏移的数据处理速度，节省了运行时间。
图11是本申请一个实施例中规则化处理模块的结构示意图，如图11所示，在上述实施例的基础上，本申请一个实施例中，所述规则化处理模块104可以包括：
OVT计算单元111，可以用于对所述述静校正处理数据进行OVT计算，获取所述OVT信息，将所述OVT信息添加到所述静校正处理数据中；
规则化处理单元112，可以用于利用五维叠前面元数据规则化处理方法，对添加所述OVT信息的静校正处理数据内进行规则化处理，获得所述OVT域偏移地震数据。
本申请提供的叠前时间偏移数据处理装置，将OVT计算应用到传统的五维叠前面元数据规则化的处理流程中，用以保护来自不同方位的地下反射信息。将地震数据进行了OVT域的分选，为后续叠前时间偏移提供了准确的数据基础。
图12是本申请又一个实施例中规则化处理模块的结构示意图，如图12所示，在上述实施例的基础上，本申请一个实施例中，所述规则化处理模块104可以包括：
OVT二次计算121，可以用于对添加所述OVT信息的静校正处理数据进行规则化处理后，对所述规则化处理后的地震数据进行OVT计算，将获得的OVT信息添加到所述规则化处理后的地震数据中，获得所述OVT域偏移地震数据。
本申请提供的叠前时间偏移数据处理装置，通过两次OVT计算，将OVT信息叠加到地震数据的道头信息中，为叠前时间偏移处理提供了准确的数据基础。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种叠前时间偏移数据处理装置，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
根据采集到的宽方位地震数据，获得OVT域地震数据；
将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
将所述静校正处理数据进行速度分析，获得均方根速度场；
在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书实施例提供的上述叠前时间偏移数据处理方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在PC端实现、linux系统实现，或其他例如使用android、iOS系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种叠前时间偏移数据处理系统的一个实施例中，图13是本申请提供的一种叠前时间偏移数据处理系统实施例的模块结构示意图，如图13所示，本申请另一实施例提供的叠前时间偏移数据处理系统可以包括处理器131以及用于存储处理器可执行指令的存储器132，
处理器131和存储器132通过总线133完成相互间的通信；
所述处理器131用于调用所述存储器132中的程序指令，以执行上述各叠前时间偏移数据处理方法实施例所提供的方法，例如包括：根据采集到的宽方位地震数据，获得OVT域地震数据；将所述OVT域地震数据进行静校正处理，获得静校正处理数据；将所述静校正处理数据进行速度分析，获得均方根速度场；在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。
在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 (10)
1.一种叠前时间偏移数据处理方法，其特征在于，包括：
根据采集到的宽方位地震数据，获得偏移距矢量片OVT域地震数据；
将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
将所述静校正处理数据进行速度分析，获得均方根速度场；
在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理；
其中，所述进行叠前时间偏移处理，包括采用下述公式计算波场值：
上式中，U表示下行波场值，ω表示角频率，v表示速度，z表示在时间偏移中对应的是时间，
表示水平方向波数在方位角上的投影，h表示半炮检距，且
而
kx，ky是X和Y两个方向的波数，b表示一个倾角时差校正的距离，A表示叠前偏移算子响应的椭圆的长半轴，B表示共中心点与b的射线路径长度。
2.如权利要求1所述的一种叠前时间偏移数据处理方法，其特征在于，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
对所述述静校正处理数据进行OVT计算，获取所述OVT信息，将所述OVT信息添加到所述述静校正处理后的OVT域地震数据中；
利用五维叠前面元数据规则化处理方法，对添加所述OVT信息的静校正处理数据进行规则化处理，获得所述OVT域偏移地震数据。
3.如权利要求2所述的一种叠前时间偏移数据处理方法，其特征在于，所述对添加所述OVT信息的静校正处理数据进行规则化处理，获得所述OVT域偏移地震数据，包括：
对添加所述OVT信息的静校正处理数据进行规则化处理后，对所述规则化处理后的地震数据进行OVT计算，将获得的OVT信息添加到所述规则化处理后的地震数据中，获得所述OVT域偏移地震数据。
4.如权利要求1所述的一种叠前时间偏移数据处理方法，其特征在于，所述方法还包括：
将所述静校正处理数据进行噪音处理，获得噪音处理数据；
相应地，所述将所述静校正处理数据进行速度分析，获得均方根速度场，包括：
将所述噪音处理数据进行速度分析，获得所述均方根速度场；
相应地，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
在所述噪音处理数据中添加所述OVT信息，对添加所述OVT信息的噪音处理数据进行规则化处理，获得OVT域偏移地震数据。
5.如权利要求1所述的一种叠前时间偏移数据处理方法，其特征在于，所述方法还包括：
将所述静校正处理数据进行振幅补偿，并将所述振幅补偿后的地震数据进行反褶积处理，获得补偿处理数据；
相应地，所述将所述静校正处理数据进行速度分析，获得均方根速度场，包括：
将所述补偿处理数据进行速度分析，获得所述均方根速度场；
相应地，所述在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据，包括：
在所述补偿处理数据中添加所述OVT信息，对添加所述OVT信息的补偿处理数据进行规则化处理，获得OVT域偏移地震数据。
6.一种叠前时间偏移数据处理装置，其特征在于，包括：
宽方位数据处理模块，用于根据采集到的宽方位地震数据，获得偏移距矢量片OVT域地震数据；
静校正模块，用于将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
速度分析模块，用于将所述静校正处理数据进行速度分析，获得均方根速度场；
规则化处理模块，用于在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
叠前时间偏移模块，用于根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理；
其中，所述进行叠前时间偏移处理，包括采用下述公式计算波场值：
上式中，U表示下行波场值，ω表示角频率，v表示速度，z表示在时间偏移中对应的是时间，
表示水平方向波数在方位角上的投影，h表示半炮检距，且
而
kx，ky是X和Y两个方向的波数，b表示一个倾角时差校正的距离，A表示叠前偏移算子响应的椭圆的长半轴，B表示共中心点与b的射线路径长度。
7.如权利要求6所述的一种叠前时间偏移数据处理装置，其特征在于，所述规则化处理模块，包括：
OVT计算单元，用于对所述述静校正处理数据进行OVT计算，获取所述OVT信息，将所述OVT信息添加到所述静校正处理数据中；
规则化处理单元，用于利用五维叠前面元数据规则化处理方法，对添加所述OVT信息的静校正处理数据内进行规则化处理，获得所述OVT域偏移地震数据。
8.如权利要求7所述的一种叠前时间偏移数据处理装置，其特征在于，所述规则化处理模块还包括：
OVT二次计算，用于对添加所述OVT信息的静校正处理数据进行规则化处理后，对所述规则化处理后的地震数据进行OVT计算，将获得的OVT信息添加到所述规则化处理后的地震数据中，获得所述OVT域偏移地震数据。
9.一种叠前时间偏移数据处理装置，其特征在于，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
根据采集到的宽方位地震数据，获得偏移距矢量片OVT域地震数据；
将所述OVT域地震数据进行静校正处理，获得静校正处理数据；
将所述静校正处理数据进行速度分析，获得均方根速度场；
在所述静校正处理数据中添加OVT信息，对添加所述OVT信息的静校正处理数据进行规则化处理，获得OVT域偏移地震数据；
根据所述OVT域偏移地震数据和所述均方根速度场，进行叠前时间偏移处理。
10.一种叠前时间偏移数据处理系统，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现如权利要求1至5中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810256518.XA CN108828668B (zh) | 2018-03-27 | 2018-03-27 | 一种叠前时间偏移数据处理方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810256518.XA CN108828668B (zh) | 2018-03-27 | 2018-03-27 | 一种叠前时间偏移数据处理方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108828668A CN108828668A (zh) | 2018-11-16 |
| CN108828668B true CN108828668B (zh) | 2020-04-10 |
Family
ID=64154216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810256518.XA Active CN108828668B (zh) | 2018-03-27 | 2018-03-27 | 一种叠前时间偏移数据处理方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108828668B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884692B (zh) * | 2019-04-10 | 2020-11-10 | 中煤科工集团西安研究院有限公司 | 一种煤田采空区探测中的四维地震数据处理方法及系统 |
| CN110187383B (zh) * | 2019-05-27 | 2020-09-29 | 中海石油(中国)有限公司 | 一种海上宽方位地震数据cov道集快速分选方法 |
| CN112578452B (zh) * | 2019-09-27 | 2024-06-18 | 中国石油化工股份有限公司 | 一种gpu加速的ovt域叠前深度偏移方法及系统 |
| CN111813551A (zh) * | 2020-07-08 | 2020-10-23 | 上海雪湖科技有限公司 | 用于地震勘探资料处理需求的叠前时间偏移调度系统 |
| CN115184994B (zh) * | 2022-06-08 | 2023-04-18 | 北京东方联创地球物理技术有限公司 | 一种多方位角采集的三维地震数据融合处理方法 |
| CN115963565B (zh) * | 2023-01-20 | 2023-06-16 | 中国地质大学（北京） | 基于双曲矢量中值滤波器的海底多次波衰减方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104280768A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种适用于逆时偏移的吸收边界条件方法 |
| CN106405641A (zh) * | 2016-10-26 | 2017-02-15 | 中国石油天然气集团公司 | 叠前时间偏移方法和装置 |
| CN109425896A (zh) * | 2017-08-25 | 2019-03-05 | 中国石油天然气股份有限公司 | 白云岩油气储层分布预测方法及装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101581792B (zh) * | 2008-05-15 | 2011-06-22 | 中国石油天然气集团公司 | 一种提高纵波勘探横向分辨率的方法 |
| US9733371B2 (en) * | 2013-09-05 | 2017-08-15 | Exxonmobil Upstream Research Company | Creating seismic images using expanded image gathers |
-
2018
- 2018-03-27 CN CN201810256518.XA patent/CN108828668B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104280768A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种适用于逆时偏移的吸收边界条件方法 |
| CN106405641A (zh) * | 2016-10-26 | 2017-02-15 | 中国石油天然气集团公司 | 叠前时间偏移方法和装置 |
| CN109425896A (zh) * | 2017-08-25 | 2019-03-05 | 中国石油天然气股份有限公司 | 白云岩油气储层分布预测方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 杨兴圣，等.高密度宽方位地震处理技术在王集地区的应用.《石油地质与工程》.2016, * |
| 高密度宽方位地震处理技术在王集地区的应用;杨兴圣，等;《石油地质与工程》;20160331;36-38 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108828668A (zh) | 2018-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108828668B (zh) | 一种叠前时间偏移数据处理方法及装置 | |
| US10386519B2 (en) | Automated near surface analysis by surface-consistent refraction methods | |
| US7768872B2 (en) | Offset-azimuth binning for migration and velocity analysis | |
| Liu et al. | The issues of prestack reverse time migration and solutions with graphic processing unit implementation | |
| Maupin | Upper-mantle structure in southern Norway from beamforming of Rayleigh wave data presenting multipathing | |
| US12007515B2 (en) | Optimal survey design | |
| WO2006054181A1 (en) | Method for processing at least two sets of seismic data | |
| CN107193045B (zh) | 一种地震数据处理方法及装置 | |
| US10795039B2 (en) | Generating pseudo pressure wavefields utilizing a warping attribute | |
| WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN112394414A (zh) | 两步法地震绕射波场叠前分离的工艺 | |
| Zhang et al. | Improving receiver function imaging with high-resolution Radon transform | |
| Perrone et al. | Linearized wave-equation migration velocity analysis by image warping | |
| Decker et al. | Quantifying and correcting residual azimuthal anisotropic moveout in image gathers using dynamic time warping | |
| US8380440B2 (en) | 3D residual binning and flatness error correction | |
| US20170176616A1 (en) | System and method for dip-guided seismic image stacking | |
| Smith et al. | Making seismic monitoring work in a complex desert environment—4D processing | |
| GB2425355A (en) | Residual move out processing of seismic data | |
| CN114002742B (zh) | 一种欧拉高斯束偏移成像方法及装置 | |
| CN109471173B (zh) | 一种剩余静校正方法及装置 | |
| CN113743193A (zh) | 一种叠前地震资料线性干扰压制方法及系统 | |
| CN114114421A (zh) | 基于深度学习的导向自学习地震数据去噪方法及装置 | |
| Zand et al. | Integrated algorithm for high‐resolution crustal‐scale imaging using complementary OBS and streamer data | |
| Li et al. | An efficient deep learning method for VSP wavefield separation: A DAS-VSP case | |
| CN103852786B (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 |