CN107976716A - 一种远偏移距转换波动校正方法及装置 - Google Patents
一种远偏移距转换波动校正方法及装置 Download PDFInfo
- Publication number
- CN107976716A CN107976716A CN201711203741.XA CN201711203741A CN107976716A CN 107976716 A CN107976716 A CN 107976716A CN 201711203741 A CN201711203741 A CN 201711203741A CN 107976716 A CN107976716 A CN 107976716A
- Authority
- CN
- China
- Prior art keywords
- wave
- offset
- converted
- travelling
- zero
- 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
- 238000000034 method Methods 0.000 title claims abstract description 108
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 100
- 230000009471 action Effects 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 46
- 238000012545 processing Methods 0.000 claims description 26
- 238000012937 correction Methods 0.000 claims description 25
- 230000008569 process Effects 0.000 abstract description 34
- 238000004364 calculation method Methods 0.000 abstract description 14
- 230000008859 change Effects 0.000 abstract description 6
- 238000003860 storage Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 13
- 230000006872 improvement Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 238000013508 migration Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 239000000047 product Substances 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000000205 computational method Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 239000000463 material 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
- 230000000644 propagated effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000003068 static effect Effects 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
技术领域
本申请属于地震勘探处理技术领域，尤其涉及一种远偏移距转换波动校正方法及装置。
背景技术
多波地震勘探是利用纵波、横波、转换波等多种信息对油气藏进行精细勘察的一种勘探技术。这种技术可以识别薄互层、小断层、低幅度构造，准确预测储层边界和储层物性、监测流体，甚至还可以直接用于油气预测等。进行地震勘探时需要经过地震数据采集、数据处理和地震资料解释等过程。数据处理包括多种地震波的数据处理，例如：横波、纵波、转换波等。由于转换波的信噪比低，频率低，速度低，转换波的处理是地震数据处理的一项技术难题。
转换波的速度与单纯纵波或者横波不同，转换波的速度不符合双曲线规律。由于转换波的信噪比低，频率低，速度低，速度不符合双曲线规律，导致转换波的初始速度难于把握，同相轴难以较平，难于形成一个好的叠加剖面，给后续进行地震资料解释带来不便。进行转换波的地震数据处理时，需要进行转换波动校正。动校正是地震数据处理的重要组成部分，动校正是将非零炮检距上所接收到的反射波旅行时校正到零炮检距上的自激自收反射波旅行时。进行动校正的目的是消除炮检距对反射波旅行时的影响，将共深度点反射波时距曲线校成一条直线形成同相轴。现有技术中，进行转换波动校正的方法通常是根据转换波时距曲线理论公式，计算不同偏移距对应的旅行时差，根据旅行时差进行转换波动校正。但是，转换波时距曲线理论公式，比较复杂，进行转换波动校正的计算量比较大，过程比较复杂。并且，现有技术中，转换波动校正的方法主要针对的是全偏移距数据，不具有针对性，对于远偏移距转换波的校正没有针对性的方法。因此，如何提供一种方案，能够减少转换波动校正的计算量，简化转换波动校正的计算过程，尤其是远偏移距转换波的动校正，是亟待解决的技术问题。
发明内容
本申请目的在于提供一种远偏移距转换波动校正方法及装置，根据纵波数据和转换波数据分析获得转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时。通过对比转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时，获得不同偏移距的偏移距旅行时差，根据偏移距旅行时差完成转换波动校正。减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。
一方面本申请提供了一种远偏移距转换波动校正方法，包括：
根据采集到的地震数据，获取纵波数据和转换波数据；
根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时；
根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时；
根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差；
根据所述偏移距旅行时差对所述转换波数据进行动校正，以校平偏移距大于预设阈值的转换波数据。
进一步地，所述方法的另一个实施例中，所述根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时，包括：
根据所述纵波数据获取纵波速度和纵波零偏移距旅行时；
根据所述转换波数据获取转换波零偏移距旅行时；
根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
进一步地，所述方法的另一个实施例中，所述根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时，包括：
根据所述纵波速度、所述纵波零偏移距旅行时计算地层厚度；
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
进一步地，所述方法的另一个实施例中，所述根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时，包括：
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，利用计算所述转换波横波段零偏移距旅行时；
上式中，ts0表示所述转换波横波段零偏移距旅行时，tps表示所述转换波零偏移距旅行时，h表示所述地层厚度，vp表示所述纵波速度。
进一步地，所述方法的另一个实施例中，所述根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时，包括：
根据所述纵波速度、所述地层厚度、所述转换波横波段零偏移距旅行时，利用计算不同偏移距对应的所述转换波偏移距旅行时；
上式中，t表示偏移距x对应的所述转换波偏移距旅行时，vp表示所述纵波速度，h表示所述地层厚度，ts0表示所述转换波横波段零偏移距旅行时。
进一步地，所述方法的另一个实施例中，所述根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差，包括：
将所述转换波偏移距旅行时与所述转换波零偏移距旅行时的差作为所述偏移距旅行时差。
另一方面，本申请提供了一种远偏移距转换波动校正装置，包括：
数据获取模块，用于根据采集到的地震数据，获取纵波数据和转换波数据；
零偏移距旅行时计算模块，用于根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时；
偏移距旅行时计算模块，用于根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时；
旅行时差计算模块，用于根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差；
动校正模块，用于根据所述偏移距旅行时差对所述转换波数据进行动校正，以校平偏移距大于预设阈值的转换波数据。
进一步地，所述装置的另一实施例中，所述零偏移距旅行时计算模块包括：
纵波数据处理模块，用于根据所述纵波数据获取纵波速度和纵波零偏移距旅行时；
转换波数据处理模块，用于根据所述转换波数据获取转换波零偏移距旅行时；
旅行时计算模块，用于根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
进一步地，所述装置的另一个实施例中，所述旅行时计算模块具体用于：
根据所述纵波速度、所述纵波零偏移距旅行时计算地层厚度；
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
进一步地，所述装置的另一个实施例中，所述旅行时计算模块具体用于：
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，利用计算所述转换波横波段零偏移距旅行时；
上式中，ts0表示所述转换波横波段零偏移距旅行时，tps表示所述转换波零偏移距旅行时，h表示所述地层厚度，vp表示所述纵波速度。
进一步地，所述装置的另一个实施例中，所述旅行时计算模块具体用于：
根据所述纵波速度、所述地层厚度、所述转换波横波段零偏移距旅行时，利用计算不同偏移距对应的所述转换波偏移距旅行时；
上式中，t表示偏移距x对应的所述转换波偏移距旅行时，vp表示所述纵波速度，h表示所述地层厚度，ts0表示所述转换波横波段零偏移距旅行时。
进一步地，所述装置的另一个实施例中，所述旅行时差计算模块具体用于：
将所述转换波偏移距旅行时与所述转换波零偏移距旅行时的差作为所述偏移距旅行时差。
再一方面，本申请还提供了一种远偏移距转换波动校正装置，包括：处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述远偏移距转换波动校正方法。
本申请提供的远偏移距转换波动校正方法及装置，通过对采集到的地震数据进行分析获得纵波数据和转换波数据，进一步根据纵波数据和转换波数据获得纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时。再利用纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时，计算获得不同偏移距对应的转换波偏移距旅行时。将不同偏移距对应的转换波偏移距旅行时与转换波零偏移距旅行时进行对比，可以获得偏移距旅行时差，利用偏移距旅行时差完成转换波数据的动校正。纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时可以根据获得的纵波数据和转换波数据简单分析获得，不需要大量的计算过程。不需要入射角等复杂参数的计算，不需要复杂的公式计算，只需要简单的数据和计算过程。减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。同时，采用本申请的方法进行转换波的动校正，可以校平远偏移距转换波数据，而不影响近偏移距转换波数据，具有更好的针对性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的一种远偏移距转换波动校正方法一个实施例的方法流程示意图；
图2是现有技术中理论转换波时距曲线和转换波动校正效果示意图；
图3是本申请一个实施例中的理论转换波时距曲线和转换波动校正效果示意图；
图4是本申请提供的远偏移距转换波动校正装置一个实施例的模块结构示意图；
图5是本申请一个实施例中零偏移距旅行时计算模块的结构示意图；
图6是本申请提供的另一种远偏移距转换波动校正装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
在进行地震勘探时，可以预先在地表设置激发点和检波点。以人工方法激发地震波，在各个检波点处接收地震数据，通过对地震数据的处理和解释，完成地震勘探。在地表以人工方法激发地震波，地震波在向地下传播时，遇有介质性质不同的岩层分界面，收到的地震波信号与震源特性、检波点的位置、地震波经过的地下岩层的性质和结构有关。通过对地震波记录进行处理和解释，可以推断地下岩层的性质和形态，实现地震勘探。
各个检波点接收到的地震数据中，包括不同地震波的数据，例如：包括纵波数据、横波数据和转换波数据等。与入射波型相同的波称为同类波，波型改变的如则称为转换波。如：若入射波为纵波，经地下岩层进行传播，出射波变为横波，则波形发生转变，检波点接收到的地震波称为转换波。在地震勘探过程中，从激发点到地面各观测点即检波点的距离称为偏移距，又称炮检距。偏移距常常分解为两个分量：垂直偏移距，即以直角到排列线的距离；纵偏移距，从激发点在排列线的投影到第一个检波器组中心的距离。本申请实施例中，可以将纵横比大于1的转换波数据称为远偏移距数据，即可以将纵偏移距与垂直偏移距之间的比值大于1的转换波数据称为远偏移距转换波数据。自激自收称为零炮检距又称零偏移距。各接收点属于同一相位振动的连线即地震记录上各道振动相位相同的极值(俗称波峰或波谷)的连线，称为同相轴。
由于炮检距即偏移距的不同，对各个检波点接收到的反射波的旅行时和相位会有影响，进一步影响后续地震资料的处理和解释。在地震勘探时，可以将非零炮检距上所接收到的反射波旅行时校正到零炮检距上的自激自收反射波旅行时，这一过程称为动校正。转换波动校正可以包括将非零炮检距上所接收到的转换波旅行时校正到零炮检距上的自激自收转换波旅行时，将转换波时距曲线较平，以获得同相轴的转换波，方便后续数据叠加处理。
图1是本申请提供的一种远偏移距转换波动校正方法一个实施例的方法流程示意图，本申请提供的远偏移距转换波动校正方法包括：
S1、根据采集到的地震数据，获取纵波数据和转换波数据。
具体地，本申请一个实施例中可以通过在地表以人工方法激发地震波，通过预先设置的检波点接收地震数据。可以利用三分量地震采集方法采集地震数据，三分量地震勘探可以包括利用纵波激发，通过三分量检波器记录地震数据。三分量地震勘探通常包括纵波数据和转换波数据。可以将地面接收到的地震数据通过数学变换获取纵波数据和转换波数据。进一步对纵波数据和转换波数据进行分析处理，可以获得纵波零偏移距旅行时、纵波速度即纵波的传播速度以及转换波零偏移距旅行时。
纵波零偏移距旅行时即纵波的零偏移距旅行时，纵波零偏移距旅行时可以包括激发点处采集到的纵波的时间，可以通过分析激发点采集到的地震数据获得。转换波零偏移距旅行时可以包括激发点处接收到转换波的时间，可以通过分析激发点处接收到的地震数据获得。如：可以直接根据激发点接收到纵波数据和转换波数据的时间，进行纵波数据和转换波数据的同相轴对比，获得纵波零偏移距旅行时和转换波零偏移距旅行时。
S2、根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时。
具体地，当入射波是纵波，经过地下岩层传播后，以横波的形式出射，则检波器接收到的波形可以认为是转换波。转换波旅行时可以包括入射的纵波段旅行时和出射的横波段旅行时，转换波横波段零偏移距旅行时可以包括激发点接收到的转换波横波段的时间。可以将纵波和转换波进行道集同相轴对比，获取转换波横波段零偏移距旅行时。如：可以获取纵波零偏移距旅行时和同相轴的转换波零偏移距旅行时，通过对比分析纵波零偏移距旅行时和同相轴的转换波零偏移距旅行时，获得转换波横波段零偏移距旅行时。当然，还可以使用其他方法，获得转换波横波段零偏移距旅行时，本申请不作具体限定。
S3、根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时。
具体地，根据纵波数据和转换波横波段零偏移距旅行时，可以计算不同偏移距对应的转换波偏移距旅行时。具体方法可以包括：根据纵波速度计算出地层厚度，利用地层厚度和激发点到检波点之间的距离即偏移距，可以计算岩心到检波点处的距离。根据岩心到检波点处的距离和纵波速度，结合获得转换波横波段零偏移距旅行时，可以计算不同偏移距对应的转换波偏移距旅行时。当然，还可以采用其他方法计算不同偏移距对应的转换波偏移距旅行时，本申请不作具体限定。例如：可以根据纵波数据获得纵波零偏移距旅行时和纵波速度，根据纵波零偏移距旅行时、纵波速度和获得的转换波横波段零偏移距旅行时，可以估算转换波的传播速度。利用估算的转换波的传播速度和偏移距，可以计算不同偏移距对应的转换波偏移距旅行时。
S4、根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差。
具体地，计算获得不同偏移距对应的转换波偏移距旅行时后，可以根据转换波数据获得转换波零偏移距旅行时。对比转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时，获得同偏移距对应的偏移距旅行时差。具体可以将转换波偏移距旅行时与转换波零偏移距旅行时的差作为偏移距旅行时差，即不同偏移距对应一个偏移距旅行时差。
S5、根据所述偏移距旅行时差对所述转换波数据进行动校正，以校平偏移距大于预设阈值的转换波数据。
具体地，获得偏移距旅行时差后，对转换波数据进行动校正。具体可以将转换波数据中行不同偏移距处采集到的转换波数据的时间减去偏移距旅行时差，作为动校正的结果。例如：若某一偏移距处的检波点采集到的转换波数据的时间对应为2800毫秒，本申请一个实施例计算获得该偏移距对应的偏移距旅行时差为300毫秒，则可以将该偏移距该时刻的转换波数据校正为2500毫秒。本申请的转换波动校正方法，可以校平偏移距大于预设阈值的转换波数据，即可以校平远偏移距转换波数据。其中，预设阈值可以根据实际需要获得，本申请实施例中可以将预设阈值取为1，即可以校平纵横比大于1的转换波数据。
本申请提供的远偏移距转换波动校正方法，对采集到的地震数据进行分析，获得纵波数据和转换波数据，并根据纵波数据和转换波数据获得转换波横波段零偏移距旅行时，进一步获得不同偏移距对应的转换波偏移距旅行时。根据不同偏移距对应的转换波偏移距旅行时和转换波零偏移距旅行时，可以获得不同偏移距对应的偏移距旅行时差。根据不同偏移距对应的偏移距旅行时差，进行转换波动校正。不需要计算获取地震波的入射角等数据参数，也不需要复杂的计算过程。只需要纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时，就可以实现转换波的动校正。并且纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时可以通过地震数据的简单分析获得。减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。同时，采用本申请的方法进行转换波的动校正，可以校平远偏移距转换波数据，而不影响近偏移距转换波数据，具有更好的针对性。
在上述实施例的基础上，步骤S2即所述根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时，包括：
根据所述纵波数据获取纵波速度和纵波零偏移距旅行时；
根据所述转换波数据获取转换波零偏移距旅行时；
根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
具体地，在本申请一个实施例中，可以通过对纵波数据进行速度分析，获得纵波速度，同时可以对纵波数据进行处理，获得纵波零偏移距旅行时即激发点处采集到的纵波的时间，具体可以通过分析激发点采集到的地震数据获得纵波零偏移距旅行时。还可以通过分析转换波数据获取转换波零偏移距旅行时即激发点处接收到转换波的时间。因为转换波旅行时可以包括入射的纵波段旅行时和出射的横波段旅行时，因此，转换波横波段零偏移距旅行时与纵波零偏移距旅行时以及转换波零偏移距旅行时有关。可以根据获得的纵波速度、纵波零偏移距旅行时以及转换波零偏移距旅行时，计算转换波横波段零偏移距旅行时。具体可以通过分析激发点处接收到的地震数据，将接收到的纵波数据和转换波数据进行道集同相轴对比，获得转换波零偏移距横波段旅行时。
在本申请一个实施例中，所述根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时，包括：
根据所述纵波速度、所述纵波零偏移距旅行时计算地层厚度；
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
具体地，可以根据获得的纵波速度和纵波零偏移距旅行时计算地层厚度，具体可以利用如下公式(1)计算地层厚度：
h＝vp×tp0/2 (1)
上式中，h可以表示地层厚度，vp可以表示纵波速度，tp0可以表示纵波零偏移距旅行时。
计算获得地层厚度后，可以根据地层厚度、纵波速度以及转换波零偏移距旅行时，计算转换波横波段零偏移距旅行时。具体可以根据地层厚度和纵波速度计算入射的纵波段的零偏移距旅行时，将获得转换波零偏移距旅行时减去入射的纵波段的零偏移距旅行时，可以获得转换波横波段零偏移距旅行时。在本申请一个实施例中，可以采用如下公式(2)计算获得转换波横波段零偏移距旅行时：
上式中，ts0可以表示转换波横波段零偏移距旅行时，tps可以表示转换波零偏移距旅行时，h可以表示地层厚度，vp可以表示纵波速度。
由于转换波旅行时包括入射的纵波段旅行时和出射的横波段旅行时，在获取到纵波零偏移距旅行时和转换波零偏移距旅行时后，也可以直接利用转换波零偏移距旅行时和纵波零偏移距旅行时计算获得示转换波横波段零偏移距旅行时。当入射的纵波经过地下岩层反射转换成横波，激发点处获得的转换波横波段零偏移距旅行时可以认为包括：转换波零偏移距旅行时与入射纵波段零偏移距旅行时的差。此时，入射纵波段零偏移距旅行时可以为纵波零偏移距旅行时的一半。具体可以结合公式(1)和公式(2)获得如下公式(3)：
上式中，各参数的含义与公式(1)和(2)相同，此处不再赘述。
在本申请一个实施例中，获得转换波横波段零偏移距旅行时后，可以根据所述纵波速度、所述地层厚度、所述转换波横波段零偏移距旅行时，利用如下公式(4)计算获得不同偏移距对应的所述转换波偏移距旅行时：
上式中，t可以表示偏移距x对应的转换波偏移距旅行时，vp可以表示纵波速度，h可以表示地层厚度，ts0可以表示转换波横波段零偏移距旅行时。
获得不同偏移距对应的转换波偏移距旅行时后，对比不同偏移距对应的转换波偏移距旅行时和转换波零偏移距旅行时，获得偏移距旅行时差，具体可以参照如下公式(5)：
Δt＝t-tps (5)
上式中，Δt可以表示偏移距旅行时差，t可以表示不同偏移距对应的转换波偏移距旅行时，tps可以表示转换波零偏移距旅行时。
获得不同偏移距对应的偏移距旅行时差Δt后，可以对相应时刻的采样点的转换波数据进行校正。具体校正方法可以参照上述实施例。
需要说明的是，本申请中的各个公式只是一种实施例，还可以根据需要对公式进行变形或变换，本申请不作具体限定。
图2是现有技术中理论转换波时距曲线和转换波动校正效果示意图，图3是本申请一个实施例中的理论转换波时距曲线和转换波动校正效果示意图，图2和图3中水平轴是偏移距，从左到右增大，纵轴是时间，从上到下增大。图2中下倾的是理论转换波时距曲线，而水平的是理论上较平后的转换波时距曲线，是完全水平的。图3中，下倾的是理论的转换波时距曲线，几乎水平的是本申请远偏移距转换波动校正方法较平的转换波时距曲线。对比图2和图3可以看出，本申请的远偏移距转换波动校正方法能够在很大偏移距范围内校平转换波时距曲线，与理论上转换波动校正的效果几乎相同，效果明显。即采用本申请的方法进行转换波的动校正，可以校平远偏移距转换波数据，而不影响近偏移距转换波数据，具有更好的针对性。
本申请提供的远偏移距转换波动校正方法，通过对采集到的地震数据进行分析获得纵波数据和转换波数据，进一步根据纵波数据和转换波数据获得纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时。再利用纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时，计算获得不同偏移距对应的转换波偏移距旅行时。将不同偏移距对应的转换波偏移距旅行时与转换波零偏移距旅行时进行对比，可以获得偏移距旅行时差，利用偏移距旅行时差完成转换波数据的动校正。纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时，可以根据获得的纵波数据和转换波数据简单分析获得，不需要大量的计算过程。不需要复杂的公式计算，只需要简单的数据和计算过程，减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。
基于上述所述的远偏移距转换波动校正方法，本说明书一个或多个实施例还提供一种远偏移距转换波动校正装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
具体地，图4是本申请提供的远偏移距转换波动校正装置一个实施例的模块结构示意图，如图4所示，本申请中提供的远偏移距转换波动校正装置包括：数据获取模块41、零偏移距旅行时计算模块42、偏移距旅行时计算模块43、旅行时差计算模块44和动校正模块45。
数据获取模块41，可以用于根据采集到的地震数据，获取纵波数据和转换波数据；
零偏移距旅行时计算模块42，可以用于根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时；
偏移距旅行时计算模块43，可以用于根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时；
旅行时差计算模块44，可以用于根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差；
动校正模块45，可以用于根据所述偏移距旅行时差对所述转换波数据进行动校正，以校平偏移距大于预设阈值的转换波数据。
本申请提供的远偏移距转换波动校正装置，利用纵波数据和转换波数据，获得转换波零偏移距旅行时、转换波横波段零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时。通过对比转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时，获得转换波的偏移距旅行时差，根据偏移距旅行时差完成转换波数据的动校正。转换波零偏移距旅行时、转换波横波段零偏移距旅行时和不同偏移距对应的偏移距旅行时可以通过纵波数据和转换波数据简单的分析计算获得，不需要复杂的计算过程。也不需要计算地震波的入射角等数据参数，减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。
图5是本申请一个实施例中零偏移距旅行时计算模块的结构示意图，如图5所示，在上述实施例的基础上，所述零偏移距旅行时计算模块42包括：
纵波数据处理模块51，可以用于根据所述纵波数据获取纵波速度和纵波零偏移距旅行时；
转换波数据处理模块52，可以用于根据所述转换波数据获取转换波零偏移距旅行时；
旅行时计算模块53，可以用于根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
本申请提供的远偏移距转换波动校正装置，可以通过对纵波数据和转换波数据进行简单的数据分析获得纵波零偏移距旅行时、纵波速度和转换波零偏移距旅行时。根据纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时，可以计算获得转换波横波段零偏移距旅行时。不需要入射角等复杂参数的计算，参数的获取不需要复杂的计算过程，为后续转换波动校正提供了准确的数据基础。同时，减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。并且，采用本申请的方法进行转换波的动校正，可以校平远偏移距转换波数据，而不影响近偏移距转换波数据，具有更好的针对性。
在上述实施例的基础上，所述旅行时计算模块具体用于：
根据所述纵波速度、所述纵波零偏移距旅行时计算地层厚度；
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
本申请提供的远偏移距转换波动校正装置，提供了简单的计算转换波横波段零偏移距旅行时的计算方法，只需要根据纵波速度、地层厚度、转换波零偏移距旅行时可以完成计算。不需要入射角等复杂参数的计算，参数的获取不需要复杂的计算过程，为后续转换波动校正提供了准确的数据基础。同时，减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上，所述旅行时计算模块具体用于：
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，利用计算所述转换波横波段零偏移距旅行时；
上式中，ts0表示所述转换波横波段零偏移距旅行时，tps表示所述转换波零偏移距旅行时，h表示所述地层厚度，vp表示所述纵波速度。
本申请提供的远偏移距转换波动校正装置，提供了简单的计算转换波横波段零偏移距旅行时的计算方法，不需要复杂的开平方计算等，也不需要入射角等复杂参数的计算，参数的获取不需要复杂的计算过程，为后续转换波动校正提供了准确的数据基础。同时，减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上，所述旅行时计算模块具体用于：
根据所述纵波速度、所述地层厚度、所述转换波横波段零偏移距旅行时，利用计算不同偏移距对应的所述转换波偏移距旅行时；
上式中，t表示偏移距x对应的所述转换波偏移距旅行时，vp表示所述纵波速度，h表示所述地层厚度，ts0表示所述转换波横波段零偏移距旅行时。
本申请提供的远偏移距转换波动校正装置，提供了简单的计算不同偏移距对应的转换波旅行时的计算方法，不需要复杂的开平方计算等，也不需要入射角等复杂参数的计算，参数的获取不需要复杂的计算过程，为后续转换波动校正提供了准确的数据基础。同时，减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上，所述旅行时差计算模块具体用于：
将所述转换波偏移距旅行时与所述转换波零偏移距旅行时的差作为所述偏移距旅行时差。
本申请提供的远偏移距转换波动校正装置，在获得不同偏移距对应的转换波偏移距旅行时差后，通过对比转换波偏移距旅行时和转换波零偏移距旅行时，可以获得不同偏移距对应的偏移距旅行时差。根据获得的偏移距旅行时差，可以完成转换波的动校正，消除偏移距对转换波数据的影响。减少了转换波动校正的计算量，简化了转换波动校正的计算过程，进一步提高了转换波动校正的效率和准确性，为后续地震数据的处理和解释提供了更准确的数据基础。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书实施例提供的上述远偏移距转换波动校正方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在PC端实现、linux系统实现，或其他例如使用android、iOS系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种远偏移距转换波动校正装置的另一个实施例中，图6是本申请提供的另一种远偏移距转换波动校正装置实施例的模块结构示意图，如图6所示，本申请另一实施例提供的远偏移距转换波动校正装置可以包括处理器61以及用于存储处理器可执行指令的存储器62，
处理器61和存储器62通过总线63完成相互间的通信；
所述处理器61用于调用所述存储器62中的程序指令，以执行上述各远偏移距转换波动校正方法实施例所提供的方法，例如包括：根据采集到的地震数据，获取纵波数据和转换波数据；根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时；根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时；根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差；根据所述偏移距旅行时差对所述转换波数据进行动校正。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。
在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 (13)
1.一种远偏移距转换波动校正方法，其特征在于，所述方法包括：
根据采集到的地震数据，获取纵波数据和转换波数据；
根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时；
根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时；
根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差；
根据所述偏移距旅行时差对所述转换波数据进行动校正，以校平偏移距大于预设阈值的转换波数据。
2.如权利要求1所述的一种远偏移距转换波动校正方法，其特征在于，所述根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时，包括：
根据所述纵波数据获取纵波速度和纵波零偏移距旅行时；
根据所述转换波数据获取转换波零偏移距旅行时；
根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
3.如权利要求2所述的一种远偏移距转换波动校正方法，其特征在于，所述根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时，包括：
根据所述纵波速度、所述纵波零偏移距旅行时计算地层厚度；
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
4.如权利要求3所述的一种远偏移距转换波动校正方法，其特征在于，所述根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时，包括：
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，利用计算所述转换波横波段零偏移距旅行时；
上式中，ts0表示所述转换波横波段零偏移距旅行时，tps表示所述转换波零偏移距旅行时，h表示所述地层厚度，vp表示所述纵波速度。
5.如权利要求3所述的一种远偏移距转换波动校正方法，其特征在于，所述根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时，包括：
根据所述纵波速度、所述地层厚度、所述转换波横波段零偏移距旅行时，利用计算不同偏移距对应的所述转换波偏移距旅行时；
上式中，t表示偏移距x对应的所述转换波偏移距旅行时，vp表示所述纵波速度，h表示所述地层厚度，ts0表示所述转换波横波段零偏移距旅行时。
6.如权利要求2-5任一项所述的一种远偏移距转换波动校正方法，其特征在于，所述根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差，包括：
将所述转换波偏移距旅行时与所述转换波零偏移距旅行时的差作为所述偏移距旅行时差。
7.一种远偏移距转换波动校正装置，其特征在于，包括：
数据获取模块，用于根据采集到的地震数据，获取纵波数据和转换波数据；
零偏移距旅行时计算模块，用于根据所述纵波数据和所述转换波数据，计算转换波横波段零偏移距旅行时；
偏移距旅行时计算模块，用于根据所述纵波数据和所述转换波横波段零偏移距旅行时，计算不同偏移距对应的转换波偏移距旅行时；
旅行时差计算模块，用于根据所述转换波偏移距旅行时和所述转换波数据，计算不同偏移距对应的偏移距旅行时差；
动校正模块，用于根据所述偏移距旅行时差对所述转换波数据进行动校正，以校平偏移距大于预设阈值的转换波数据。
8.如权利要求7所述的一种远偏移距转换波动校正装置，其特征在于，所述零偏移距旅行时计算模块包括：
纵波数据处理模块，用于根据所述纵波数据获取纵波速度和纵波零偏移距旅行时；
转换波数据处理模块，用于根据所述转换波数据获取转换波零偏移距旅行时；
旅行时计算模块，用于根据所述纵波速度、所述纵波零偏移距旅行时、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
9.如权利要求8所述的一种远偏移距转换波动校正装置，其特征在于，所述旅行时计算模块具体用于：
根据所述纵波速度、所述纵波零偏移距旅行时计算地层厚度；
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，计算所述转换波横波段零偏移距旅行时。
10.如权利要求9所述的一种远偏移距转换波动校正装置，其特征在于，所述旅行时计算模块具体用于：
根据所述纵波速度、所述地层厚度、所述转换波零偏移距旅行时，利用计算所述转换波横波段零偏移距旅行时；
上式中，ts0表示所述转换波横波段零偏移距旅行时，tps表示所述转换波零偏移距旅行时，h表示所述地层厚度，vp表示所述纵波速度。
11.如权利要求9所述的一种远偏移距转换波动校正装置，其特征在于，所述旅行时计算模块具体用于：
根据所述纵波速度、所述地层厚度、所述转换波横波段零偏移距旅行时，利用计算不同偏移距对应的所述转换波偏移距旅行时；
上式中，t表示偏移距x对应的所述转换波偏移距旅行时，vp表示所述纵波速度，h表示所述地层厚度，ts0表示所述转换波横波段零偏移距旅行时。
12.如权利要求8-11任一项所述的一种远偏移距转换波动校正装置，其特征在于，所述旅行时差计算模块具体用于：
将所述转换波偏移距旅行时与所述转换波零偏移距旅行时的差作为所述偏移距旅行时差。
13.一种远偏移距转换波动校正装置，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现如权利要求1至6中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711203741.XA CN107976716B (zh) | 2017-11-27 | 2017-11-27 | 一种远偏移距转换波动校正方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711203741.XA CN107976716B (zh) | 2017-11-27 | 2017-11-27 | 一种远偏移距转换波动校正方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107976716A true CN107976716A (zh) | 2018-05-01 |
| CN107976716B CN107976716B (zh) | 2019-08-06 |
Family
ID=62011798
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711203741.XA Active CN107976716B (zh) | 2017-11-27 | 2017-11-27 | 一种远偏移距转换波动校正方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107976716B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4802146A (en) * | 1985-05-23 | 1989-01-31 | Mobil Oil Corporation | Method for moveout correction and stacking velocity estimation of offset VSP data |
| CN101251604A (zh) * | 2008-04-12 | 2008-08-27 | 中国石油集团西北地质研究所 | 二参数转换波速度分析及动校正方法 |
| CN103630934A (zh) * | 2012-08-23 | 2014-03-12 | 中国石油天然气集团公司 | 一种确定转换波检波点大的横波静校正量的方法 |
| CN104216016A (zh) * | 2014-08-12 | 2014-12-17 | 中国石油天然气集团公司 | 一种纵波约束扫描的转换波动校正方法以及系统 |
| CN104316967A (zh) * | 2014-08-27 | 2015-01-28 | 中国石油化工股份有限公司 | 大偏移距动校正速度时移误差补偿方法 |
| CN104678439A (zh) * | 2013-11-27 | 2015-06-03 | 中国石油化工股份有限公司 | 一种多波连续介质速度分析方法 |
-
2017
- 2017-11-27 CN CN201711203741.XA patent/CN107976716B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4802146A (en) * | 1985-05-23 | 1989-01-31 | Mobil Oil Corporation | Method for moveout correction and stacking velocity estimation of offset VSP data |
| CN101251604A (zh) * | 2008-04-12 | 2008-08-27 | 中国石油集团西北地质研究所 | 二参数转换波速度分析及动校正方法 |
| CN103630934A (zh) * | 2012-08-23 | 2014-03-12 | 中国石油天然气集团公司 | 一种确定转换波检波点大的横波静校正量的方法 |
| CN104678439A (zh) * | 2013-11-27 | 2015-06-03 | 中国石油化工股份有限公司 | 一种多波连续介质速度分析方法 |
| CN104216016A (zh) * | 2014-08-12 | 2014-12-17 | 中国石油天然气集团公司 | 一种纵波约束扫描的转换波动校正方法以及系统 |
| CN104316967A (zh) * | 2014-08-27 | 2015-01-28 | 中国石油化工股份有限公司 | 大偏移距动校正速度时移误差补偿方法 |
Non-Patent Citations (1)
| Title |
|---|
| 黄中玉 等: "三维多分量地震资料处理技术研究", 《石油物探》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107976716B (zh) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106324676B (zh) | 一种确定断层封闭性的方法及装置 | |
| CN108181649B (zh) | 一种介质分界面振幅补偿方法及装置 | |
| CN107420092B (zh) | 井下管柱泄漏点定位方法及装置 | |
| CN107797154A (zh) | 一种电成像测井图像刻度方法及装置 | |
| CN108614293A (zh) | 砂体预测方法及装置 | |
| CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
| CN107783186A (zh) | 一种近地表速度建模方法及装置 | |
| CN104199083A (zh) | 一种确定宽线地震观测系统的方法和装置 | |
| CN108089238A (zh) | 一种利用综合地球物理资料确定深层裂谷的方法及装置 | |
| CN108828668A (zh) | 一种叠前时间偏移数据处理方法及装置 | |
| CN107956465A (zh) | 基于关联井的全区多井测井曲线标准化方法及装置 | |
| CN106842304A (zh) | 一种叠前深度偏移方法及装置 | |
| CN107179546B (zh) | 一种确定速度场的方法及装置 | |
| CN108107471A (zh) | 一种分方位初至数据体的获取方法及装置 | |
| CN108469633A (zh) | 一种地层品质因子的计算方法及装置 | |
| CN108873085A (zh) | 一种沉积岩薄层电磁数据处理方法及装置 | |
| CN107329174A (zh) | 在地震剖面上显示测井信息的方法及装置 | |
| CN108693557A (zh) | 一种地震数据处理方法及装置 | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| CN104122584A (zh) | 根据地震数据确定方向性的方法及装置 | |
| CN108710149A (zh) | 一种薄储层的预测方法及装置 | |
| CN108181656B (zh) | 一种近偏移距转换波动校正方法及装置 | |
| CN116224439B (zh) | 强各向异性层状vti介质高效射线追踪方法、设备和介质 | |
| CN107976716B (zh) | 一种远偏移距转换波动校正方法及装置 | |
| CN108492014A (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 |