CN108181656B - 一种近偏移距转换波动校正方法及装置 - Google Patents

一种近偏移距转换波动校正方法及装置 Download PDF

Info

Publication number
CN108181656B
CN108181656B CN201711210811.4A CN201711210811A CN108181656B CN 108181656 B CN108181656 B CN 108181656B CN 201711210811 A CN201711210811 A CN 201711210811A CN 108181656 B CN108181656 B CN 108181656B
Authority
CN
China
Prior art keywords
offset
wave
converted
travelling
offset distance
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
Application number
CN201711210811.4A
Other languages
English (en)
Other versions
CN108181656A (zh
Inventor
贺照权
曾友爱
曹孟起
张保庆
赵贻水
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201711210811.4A priority Critical patent/CN108181656B/zh
Publication of CN108181656A publication Critical patent/CN108181656A/zh
Application granted granted Critical
Publication of CN108181656B publication Critical patent/CN108181656B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy

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 Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本申请提供一种近偏移距转换波动校正方法及装置。所述方法包括:根据纵波数据、转换波数据计算地层厚度、转换波零偏移距旅行时;根据地层厚度和转换波零偏移距旅行时,计算转换波零偏移距速度;根据地层厚度、转换波零偏移距旅行时、转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;根据目标转换波偏移距旅行时和转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差,进一步对转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据。利用本申请中各个实施例,可以简化转换波动校正的计算过程,提高转换波动校正的效率和准确性。同时,可以针对性的校平近偏移距转换波数据。

Description

一种近偏移距转换波动校正方法及装置
技术领域
本申请属于地震勘探处理技术领域,尤其涉及一种近偏移距转换波动校正方法及装置。
背景技术
多波地震勘探是利用纵波、横波、转换波等多种信息对油气藏进行精细勘察的一种勘探技术。这种技术可以识别薄互层、小断层、低幅度构造,准确预测储层边界和储层物性、监测流体,甚至还可以直接用于油气预测等。进行地震勘探时需要经过地震数据采集、数据处理和地震资料解释等过程。数据处理包括多种地震波的数据处理,例如:横波、纵波、转换波等。由于转换波的信噪比低,频率低,速度低,转换波的处理是地震数据处理的一项技术难题。
转换波的速度与单纯纵波或者横波不同,转换波的速度不符合双曲线规律。由于转换波的信噪比低,频率低,速度低,速度不符合双曲线规律,导致转换波的初始速度难于把握,同相轴难以较平,难于形成一个好的叠加剖面,给后续进行地震资料解释带来不便。进行转换波的地震数据处理时,需要进行转换波动校正。动校正是地震数据处理的重要组成部分,动校正是将非零炮检距上所接收到的反射波旅行时校正到零炮检距上的自激自收反射波旅行时。进行动校正的目的是消除炮检距对反射波旅行时的影响,将共深度点反射波时距曲线校成一条直线形成同相轴。现有技术中,进行转换波动校正的方法通常是根据转换波时距曲线理论公式,计算不同偏移距对应的旅行时差,根据旅行时差进行转换波动校正。但是,转换波时距曲线理论公式,比较复杂,进行转换波动校正的计算量比较大,过程比较复杂。并且,现有技术中,转换波动校正的方法主要针对的是全偏移距数据,不具有针对性,对于近偏移距转换波的校正没有针对性的方法。因此,提供一种方案,能够减少转换波动校正的计算量,简化转换波动校正的计算过程,尤其是近偏移距转换波的动校正,是亟待解决的技术问题。
发明内容
本申请目的在于提供一种近偏移距转换波动校正方法及装置,根据纵波数据和转换波数据分析获得转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时。通过对比转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时,获得不同偏移距的偏移距旅行时差,根据偏移距旅行时差完成转换波动校正。减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
一方面本申请提供了一种近偏移距转换波动校正方法,包括:
根据采集到的地震数据,获取纵波数据和转换波数据;
根据所述纵波数据计算地层厚度;
根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时;
根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度;
根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;
根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差对所述转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据。
进一步地,所述方法的另一个实施例中,所述预先建立的转换波旅行时模型,包括:
上式中,t表示偏移距x对应的目标转换波偏移距旅行时,c表示比例系数,0<c<1,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
进一步地,所述方法的另一个实施例中,所述根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时,包括:
S201、设置所述比例系数的初始值;
S202、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时;
S203、若判断获知所述初始转换波偏移距旅行时不满足预设条件,则调整所述比例系数,将所述调整后的比例系数的值作为所述初始值;
S204、重复步骤S202和S203,直至所述初始转换波偏移距旅行时满足所述预设条件,将所述初始转换波偏移距旅行时作为所述目标转换波偏移距旅行时。
进一步地,所述调整所述比例系数,将所述调整后的比例系数的值作为所述初始值,包括:
预先设置比例系数的增量;
若判断获知所述初始转换波偏移距旅行时不满足预设条件,则将所述比例系数的初始值增加所述增量,将增加所述增量后的比例系数作为所述比例系数的初始值。
进一步地,所述方法的另一个实施例中,所述方法还包括:
若判断获知所述比例系数增加到1,所述初始转换波偏移距旅行时仍不满足所述预设条件,则调整所述增量的大小。
进一步地,所述方法的另一个实施例中,所述初始转换波偏移距旅行时不满足预设条件,包括:
将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时后,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据;
若所述校正后的转换波数据小于校平程度阈值,则所述初始转换波偏移距旅行时不满足所述预设条件。
进一步地,所述方法的另一个实施例中,所述根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度,包括:
根据所述地层厚度和所述转换波零偏移距旅行时,利用v=2h/t0计算所述转换波零偏移距速度;
上式中,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
另一方面,本申请提供了一种近偏移距转换波动校正装置,包括:
数据采集模块,用于根据采集到的地震数据,获取纵波数据和转换波数据;
地层厚度计算模块,用于根据所述纵波数据计算地层厚度;
零偏移距旅行时计算模块,用于根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时;
零偏移距速度计算模块,用于根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度;
偏移距旅行时计算模块,用于根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;
偏移距旅行时差计算模块,用于根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差;
校正模块,用于根据所述偏移距旅行时差对所述转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据。
进一步地,所述装置的另一实施例中,所述偏移距旅行时计算模块中预先建立的转换波旅行时模型,包括:
上式中,t表示偏移距x对应的转换波偏移距旅行时,c表示比例系数,0<c<1,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
进一步地,所述装置的另一个实施例中,所述偏移距旅行时计算模块具体用于:
S201、设置所述比例系数的初始值;
S202、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时;
S203、若判断获知所述初始转换波偏移距旅行时不满足预设条件,则调整所述比例系数,将所述调整后的比例系数的值作为所述初始值;
S204、重复步骤S202和S203,直至所述初始转换波偏移距旅行时满足所述预设条件,将所述初始转换波偏移距旅行时作为所述目标转换波偏移距旅行时。
进一步地,所述装置的另一个实施例中,所述偏移距旅行时计算模块具体用于:
预先设置比例系数的增量;
若判断获知所述初始转换波偏移距旅行时不满足预设条件,则将所述比例系数的初始值增加所述增量,将增加所述增量后的比例系数作为所述比例系数的初始值。
进一步地,所述装置的另一个实施例中,所述偏移距旅行时计算模块还用于:
若判断获知所述比例系数增加到1,所述初始转换波偏移距旅行时仍不满足所述预设条件,则调整所述增量的大小。
进一步地,所述装置的另一个实施例中,所述偏移距旅行时计算模块具体用于:
将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时后,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据;
若所述校正后的转换波数据小于校平程度阈值,则所述初始转换波偏移距旅行时不满足所述预设条件。
进一步地,所述装置的另一个实施例中,所述零偏移距速度计算模块具体用于:
根据所述地层厚度和所述转换波零偏移距旅行时,利用v=2h/t0计算所述转换波零偏移距速度;
上式中,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
再一方面,本申请还提供了一种近偏移距转换波动校正装置,包括:处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现上述近偏移距转换波动校正方法。
本申请提供的近偏移距转换波动校正方法及装置,对采集到的地震数据进行分析,获得纵波数据和转换波数据。并根据纵波数据和转换波数据获得纵波速度、纵波零偏移距旅行时、转换零偏移距旅行时等,进一步获得地层厚度和转换波零偏移距旅行速度。根据获得的地层厚度、转换波零偏移距旅行速度、转换波零偏移距旅行时,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时。根据不同偏移距对应的目标转换波偏移距旅行时和转换波零偏移距旅行时,可以获得不同偏移距对应的偏移距旅行时差。根据不同偏移距对应的偏移距旅行时差,进行转换波动校正。不需要计算获取地震波的入射角等数据参数,也不需要复杂的计算过程。只需要纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时,就可以实现转换波的动校正。并且纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时可以通过地震数据的简单分析获得。同时,还可以通过调整转换波旅行时模型的比例参数,获得满足校平程度阈值的转换波校正数据。减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。同时,采用本申请的方法进行转换波的动校正,可以校平近偏移距转换波数据,而不影响远偏移距转换波数据,具有更好的针对性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的一种近偏移距转换波动校正方法一个实施例的方法流程示意图;
图2是本申请一个实施中目标转换波旅行时计算流程示意图;
图3是本申请另一实施例中近偏移距转换波动校正方法的流程示意图;
图4是本申请另一实施例中近偏移距转换波动校正方法的流程示意图;
图5是现有技术中理论转换波时距曲线和转换波动校正效果示意图;
图6是本申请一个实施例中的理论转换波时距曲线和转换波动校正效果示意图;
图7是本申请提供的近偏移距转换波动校正装置一个实施例的模块结构示意图;
图8是本申请提供的另一种近偏移距转换波动校正装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
在进行地震勘探时,可以预先在地表设置激发点和检波点。以人工方法激发地震波,在各个检波点处接收地震数据,通过对地震数据的处理和解释,完成地震勘探。在地表以人工方法激发地震波,地震波在向地下传播时,遇有介质性质不同的岩层分界面,收到的地震波信号与震源特性、检波点的位置、地震波经过的地下岩层的性质和结构有关。通过对地震波记录进行处理和解释,可以推断地下岩层的性质和形态,实现地震勘探。
各个检波点接收到的地震数据中,包括不同地震波的数据,例如:包括纵波数据、横波数据和转换波数据等。与入射波型相同的波称为同类波,波型改变的如则称为转换波。如:若入射波为纵波,经地下岩层进行传播,出射波变为横波,则波形发生转变,检波点接收到的地震波称为转换波。在地震勘探过程中,从激发点到地面各观测点即检波点的距离称为偏移距,又称炮检距。偏移距常常分解为两个分量:垂直偏移距,即以直角到排列线的距离;纵偏移距,从激发点在排列线的投影到第一个检波器组中心的距离。本申请实施例中,可以将纵横比小于1的转换波数据称为近偏移距数据,即可以将纵偏移距与垂直偏移距之间的比值小于1的转换波数据称为近偏移距转换波数据。自激自收称为零炮检距又称零偏移距。各接收点属于同一相位振动的连线,即地震记录上各道振动相位相同的极值(俗称波峰或波谷)的连线,称为同相轴。
由于炮检距即偏移距的不同,对各个检波点接收到的反射波的旅行时和相位会有影响,进一步影响后续地震资料的处理和解释。在地震勘探时,可以将非零炮检距上所接收到的反射波旅行时校正到零炮检距上的自激自收反射波旅行时,这一过程称为动校正。转换波动校正可以包括将非零炮检距上所接收到的转换波旅行时校正到零炮检距上的自激自收转换波旅行时,将转换波时距曲线较平,以获得同相轴的转换波,方便后续数据叠加处理。
图1是本申请提供的一种近偏移距转换波动校正方法一个实施例的方法流程示意图,本申请提供的近偏移距转换波动校正方法包括:
S1、根据采集到的地震数据,获取纵波数据和转换波数据。
具体地,本申请一个实施例中可以通过在地表以人工方法激发地震波,通过预先设置的检波点接收地震数据。可以利用三分量地震采集方法采集地震数据,三分量地震勘探可以包括利用纵波激发,通过三分量检波器记录地震数据。三分量地震勘探通常包括纵波数据和转换波数据。可以将地面接收到的地震数据通过数学变换获取纵波数据和转换波数据。进一步对纵波数据进行分析处理,可以获得纵波零偏移距旅行时以及纵波速度即纵波的传播速度。
纵波零偏移距旅行时即纵波的零偏移距旅行时,纵波零偏移距旅行时可以包括激发点处采集到的纵波的时间,可以通过分析激发点采集到的地震数据获得。如:可以直接根据激发点接收到纵波数据和转换波数据的时间,获得纵波零偏移距旅行时。
S2、根据所述纵波数据计算地层厚度。
具体地,获得纵波数据后,可以对纵波数据进行速度分析获得纵波速度,并且可以根据纵波数据获得纵波零偏移距旅行时。根据纵波速度和纵波零偏移距旅行时可以计算获得地层厚度,具体可以采用如下公式(1)计算获得:
h=vp×tp0/2 (1)
上式中,h可以表示地层厚度,vp可以表示纵波速度,tp0可以表示纵波零偏移距旅行时。
S3、根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时。
具体地,可以通过分析激发点处接收到的地震数据,将获取到的纵波数据和转换波数据进行道集同相轴对比,获得转换波零偏移距旅行时。可以根据地震数据获取同相轴的纵波数据和转换波数据,通过对比分析激发点处接收到纵波的时间和接收到转换波的时间,获得转换波零偏移距旅行时。
S4、根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度。
具体地,转换波零偏移距旅行时可以认为是转换波在地下岩层内通过的时间,获得地层厚度和转换波零偏移距旅行时后,可以根据转换波零偏移距旅行时和地层厚度,就算获得转换波零偏移距速度。具体可以参考如下公式(2)计算转换波零偏移距速度:
v=2h/t0 (2)
上式中,v可以表示转换波零偏移距速度,h可以表示地层厚度,t0可以表示转换波零偏移距旅行时。
S5、根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时。
具体地,转换波旅行时模型可以包括不同偏移距时转换波对应的旅行时与其他参数之间的关系,其他参数可以包括地震厚度、偏移距的大小、转换波零偏移距速度、转换波零偏移距旅行时等。本申请一个实施例中可以预先通过实验或根据现有地震数据分析不同偏移距的转换波旅行时与地层厚度、转换波零偏移距旅行时、转换波零偏移距速度等之间的关系,构建转换波旅行时模型。具体可以通图形、表格、数值模拟、函数拟合等方式,构建转换波旅行时模型。当然,还可以通过其他方法构建转换波旅行时模型,本申请不作具体限定。在根据纵波数据、转换波数据获得地层厚度、转换波零偏移距旅行时、转换波零偏移距速度后,可以利用预先建立的转换波旅行时模型,计算获得不同偏移距对应的目标转换波偏移距旅行时。
S6、根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差。
具体地,计算获得不同偏移距对应的转换波偏移距旅行时后,对比转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时,获得同偏移距对应的偏移距旅行时差。具体可以将转换波偏移距旅行时与转换波零偏移距旅行时的差作为偏移距旅行时差,即不同偏移距对应一个偏移距旅行时差。具体可以参照如下公式(3):
Δt=t-tps (3)
上式中,Δt可以表示偏移距旅行时差,t可以表示不同偏移距对应的转换波偏移距旅行时,tps可以表示转换波零偏移距旅行时。
S7、根据所述偏移距旅行时差对所述转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据。
具体地,获得偏移距旅行时差后,对转换波数据进行动校正。具体可以将转换波数据中行不同偏移距处采集到的转换波数据的时间减去偏移距旅行时差,作为动校正的结果。例如:若某一偏移距处的检波点采集到的转换波数据的时间对应为2800毫秒,本申请一个实施例计算获得该偏移距对应的偏移距旅行时差为300毫秒,则可以将该偏移距该时刻的转换波数据校正为2500毫秒。本申请的转换波动校正方法,可以校平偏移距小于预设阈值的转换波数据,即可以校平近偏移距转换波数据。其中,预设阈值可以根据实际需要获得,本申请实施例中可以将预设阈值取为1,即可以校平纵横比小于1的转换波数据。
本申请提供的近偏移距转换波动校正方法,对采集到的地震数据进行分析,获得纵波数据和转换波数据。并根据纵波数据和转换波数据获得纵波速度、纵波零偏移距旅行时、转换零偏移距旅行时等,进一步获得地层厚度和转换波零偏移距旅行速度。根据获得的地层厚度、转换波零偏移距旅行速度、转换波零偏移距旅行时,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时。根据不同偏移距对应的目标转换波偏移距旅行时和转换波零偏移距旅行时,可以获得不同偏移距对应的偏移距旅行时差。根据不同偏移距对应的偏移距旅行时差,进行转换波动校正。不需要计算获取地震波的入射角等数据参数,也不需要复杂的计算过程。只需要纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时,就可以实现转换波的动校正。并且纵波速度、纵波零偏移距旅行时和转换波零偏移距旅行时可以通过地震数据的简单分析获得。减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。同时,采用本申请的方法进行转换波的动校正,可以校平近偏移距转换波数据,而不影响远偏移距转换波数据,具有更好的针对性。
在上述实施例的基础上,本申请一个实施例中预先建立的转换波旅行时模型可以包括如下公式(4):
上式中,t可以表示偏移距x对应的目标转换波偏移距旅行时,c可以表示比例系数,0<c<1,v可以表示转换波零偏移距速度,h可以表示所述地层厚度,t0可以表示转换波零偏移距旅行时。比例系数c可以根据实际需要进行设置,例如可以预先根据实验或现有的地震数据,估算出比例系数的取值。也可以,选取多个比例系数的值,获得不同的目标转换波偏移距旅行时,进行后续的转换波动校正,比较校正效果,选取校正效果好的作为校正后的转换波数据。
根据地震数据获得地层厚度、转换波零偏移距旅行时、转换波零偏移距速度后,可以将获得的数据代入上述公式(4),计算不同偏移距x对应的目标转换波偏移距旅行时。进一步根据不同偏移距对应的目标转换波偏移距旅行时和转换波零偏移距旅行时,获得偏移距旅行时差,根据偏移距旅行时差完成转换波动校正。可以看出,本申请在进行转换波动校正时不需要地震波入射角等复杂的参数,也不需要大量的开平方等计算过程。减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
图2是本申请一个实施中目标转换波旅行时计算流程示意图,如图2所示,在上述实施例的基础上,图1中的步骤S5即所述根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时,包括:
S201、设置所述比例系数的初始值。
具体地,可以预先设置比例系数c的初始值,具体可以根据实际需要进行设置,本申请不作具体限定,只需要满足比例系数0<c<1的取值范围。
S202、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时。
具体地,可以将计算获得的地层厚度、转换波零偏移距旅行时、转换波零偏移距速度代入转换波旅行时模型如公式(4),计算获得比例系数初始值时,不同偏移距的初始转换波旅行时。
S203、若判断获知所述初始转换波偏移距旅行时不满足预设条件,则调整所述比例系数,将所述调整后的比例系数的值作为所述初始值。
具体地,计算获得比例系数初始值对应的初始转换波偏移距旅行时后,判断计算获得的初始转换波旅行时是否满足预设条件。若不满足预设条件,则调整比例系数的大小,将调整后的比例系数的值作为比例系数的初始值。
S204、重复步骤S202和S203,直至所述初始转换波偏移距旅行时满足所述预设条件,将所述初始转换波偏移距旅行时作为所述目标转换波偏移距旅行时。
具体地,重复步骤S202,根据调整后的比例系数的初始值,计算不同偏移距对应的初始转换波偏移距旅行时。判断调整后的比例系数的初始值对应的初始转换波偏移距旅行时是否满足预设条件,若不满足,则继续调整比例系数的初始值。重复步骤S202和S203,直至初始转换波偏移距旅行时满足预设条件,将此时的初始转换波偏移距旅行时作为目标偏移距旅行时。
在本申请一个实施例中,所述初始转换波偏移距旅行时不满足预设条件,包括:
将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时后,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据;
若所述校正后的转换波数据小于校平程度阈值,则所述初始转换波偏移距旅行时不满足所述预设条件。
具体地,根据获得的初始转换波偏移距旅行时完成后续的转换波动校正,获得校正后的转换波数据。判断校正后的转换波数据是否小于校平程度阈值,若小于,则可以认为初始转换波偏移距旅行时满足预设条件;若不小于,则可以认为初始转换波偏移距旅行时不满足预设条件,继续调整比例系数的大小。其中,校平程度阈值的大小可以根据实际需要进行设置,本申请不作具体限定。实际可以根据校正后的转换波数据的较平程度,设置校平程度阈值,判断初始转换波偏移距旅行时是否满足预设条件。
如:可以根据校正后的转换波时距曲线与理论转换波时距曲线之间的夹角,或校正后的转换波数据的能量团的大小,判断校正后的转换波数据是否达到校平程度阈值,即转换波时距曲线是否较平,是否满足预设条件。若是根据校正后的转换波时距曲线与理论转换波时距曲线之间的夹角确定校平程度阈值,可以将较平阈值设置为根据校正后的转换波时距曲线与理论转换波时距曲线之间的夹角的负数。还可以根据校正后的转换波数据进行叠加处理,根据叠加数据获得地震坡面图,根据地震剖面图判断转换波数据的较平程度,进一步判断初始转换波偏移距旅行时是否满足预设条件。
具体地,图3是本申请另一实施例中近偏移距转换波动校正方法的流程示意图,如图3所示,本申请进行转换波动校正的方法包括:
S301、根据采集到的地震数据,获取纵波数据和转换波数据。
S302、根据所述纵波数据计算地层厚度。
S303、根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时。
S304、根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度。
S305、设置所述比例系数的初始值。
S306、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时。
S307、将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差。
S308、根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据。
S309、判断所述校正后的转换波数据是否小于校平程度阈值,若小于则执行步骤S310,若不小于,则执行步骤S311。
S310、调整所述比例系数,将所述调整后的比例系数的值作为所述初始值,返回步骤S306。
S311、输出步骤S308获得的校正后的转换波数据。
本申请提供的近偏移距转换波动校正方法,通过不断调整转换波旅行时模型的比例系数,进行转换波数据的动校正,直至校正后的转换波数据满足校平程度阈值的要求。提高了的转换波动校正的准确性,同时不需要计算获取地震波的入射角等复杂数据参数,也不需要复杂的计算过程减少了转换波动校正的计算量,简化了转换波动校正的计算过程。
在上述实施例的基础上,所述调整所述比例系数,将所述调整后的比例系数的值作为所述初始值,包括:
预先设置比例系数的增量;
若判断获知所述初始转换波偏移距旅行时不满足预设条件,则将所述比例系数的初始值增加所述增量,将增加所述增量后的比例系数作为所述比例系数的初始值。
具体地,本申请一个实施例中,可以通过设置比例系数的增量,进行比例系数的调整,可以认为比例系数在0-1区间内进行比例系数逐步增加增量的扫描,直至获得满足较平要求的转换波数据。图4是本申请另一实施例中近偏移距转换波动校正方法的流程示意图,如图4所示,本申请进行转换波动校正的方法包括:
S401、根据采集到的地震数据,获取纵波数据和转换波数据。
S402、根据所述纵波数据计算地层厚度。
S403、根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时。
S404、根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度。
S405、设置所述比例系数的初始值c0和比例系数的增量Δc。
S406、将所述比例系数的初始值c0、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时。
S407、将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差。
S408、根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据。
S409、判断所述校正后的转换波数据是否小于校平程度阈值,若小于则执行步骤S410,若不小于,则执行步骤S412。
S410、c0=c0+Δc。
S411、c0<1是否成立。若是,则返回步骤S406,若否,则执行步骤S412。
S412、输出步骤S408获得的校正后的转换波数据。
步骤S411中,若c0<1不成立,说明比例系数增加到1,超出比例系数的取值范围,仍没有获取到满足校平程度阈值的校正后的转换波数据。此时,可以直接输出上一次循环获得的校正后的转换波数据,也可以输出没有获取到目标结果的提示,提醒操作人员进行比例系数的再次调整。
在上述实施例的基础上,在以增量的方式调整比例系数的大小时,若判断获知所述比例系数增加到1,所述初始转换波偏移距旅行时仍不满足所述预设条件,则调整所述增量的大小。即比例系数从初始值c0逐步增加增量Δc,直至增加到1,超出比例系数c的取值范围时,仍没有获得满足预设条件的初始转换波偏移距旅行时,即校正后的转换波数据仍达不到校平程度阈值。则可以通过调整增量Δc的大小,调整比例系数增加的密度。例如:当比例系数的增量设置为Δc=0.1时,没有获得满足预设条件的初始转换波偏移距旅行时,则可以将增量设置为Δc=0.01,提高扫描密度,增大获取满足校平程度阈值的校正后的转换波数据的可能性。
需要说明的是,本申请中的各个公式只是一种实施例,还可以根据需要对公式进行变形或变换,本申请不作具体限定。
图5是现有技术中理论转换波时距曲线和转换波动校正效果示意图,图6是本申请一个实施例中的理论转换波时距曲线和转换波动校正效果示意图,图5和图6中水平轴是偏移距,从左到右增大,纵轴是时间,从上到下增大。图5中下倾的是理论转换波时距曲线,而水平的是理论上较平后的转换波时距曲线,是完全水平的。图6中,下倾的是理论的转换波时距曲线,几乎水平的是本申请近偏移距转换波动校正方法较平的转换波时距曲线。对比图5和图6可以看出,本申请的近偏移距转换波动校正方法能够在一定偏移距范围内校平转换波时距曲线,尤其是转换波偏移距较小时相比理论上转换波动校正的较平结果有明显的改善。并且图6中校正后的转换波时距曲线比图5中校正后的转换波时距曲线的上翘程度低,可以认为,本申请的近偏移距转换波动校正方法较平的范围比理论方法较平的范围有所增大。
本申请提供的近偏移距转换波动校正方法,通过对采集到的地震数据进行分析获得纵波速度、纵波零偏移距旅行时、地层厚度、转换波零偏移距旅行时、转换波零偏移距速度。再根据地层厚度、转换波零偏移距旅行时、转换波零偏移距速度,利用预先建立的转换波旅行时模型,通过调整比例系数的值,计算获得满足预设条件的不同偏移距对应的转换波偏移距旅行时。将不同偏移距对应的转换波偏移距旅行时与转换波零偏移距旅行时进行对比,可以获得偏移距旅行时差,利用偏移距旅行时差完成转换波数据的动校正。纵波速度、纵波零偏移距旅行时、地层厚度、转换波零偏移距旅行时、转换波零偏移距速度,可以根据获得的纵波数据和转换波数据简单分析获得,不需要大量的计算过程。不需要复杂的公式计算,只需要简单的数据和计算过程,减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。同时,采用本申请的方法进行转换波的动校正,可以校平近偏移距转换波数据,而不影响远偏移距转换波数据,具有更好的针对性。
基于上述所述的近偏移距转换波动校正方法,本说明书一个或多个实施例还提供一种近偏移距转换波动校正装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思,本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似,因此本说明书实施例具体的装置的实施可以参见前述方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
具体地,图7是本申请提供的近偏移距转换波动校正装置一个实施例的模块结构示意图,如图7所示,本申请中提供的近偏移距转换波动校正装置包括:数据获取模块71、地层厚度计算模块72、零偏移距旅行时计算模块73、零偏移距速度计算模块74、偏移距旅行时计算模块75、偏移距旅行时差计算模块76和校正模块77。
数据采集模块71,可以用于根据采集到的地震数据,获取纵波数据和转换波数据;
地层厚度计算模块72,可以用于根据所述纵波数据计算地层厚度;
零偏移距旅行时计算模块73,可以用于根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时;
零偏移距速度计算模块74,可以用于根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度;
偏移距旅行时计算模块75,可以用于根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;
偏移距旅行时差计算模块76,可以用于根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差;
校正模块77,可以用于根据所述偏移距旅行时差对所述转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据。
本申请提供的近偏移距转换波动校正装置,利用纵波数据和转换波数据,获得转换波零偏移距旅行时、转换波零偏移距速度、地层厚度,利用预先建立的转换波旅行时模型,可以计算不同偏移距对应的偏移距旅行时。通过对比转换波零偏移距旅行时和不同偏移距对应的转换波偏移距旅行时,获得转换波的偏移距旅行时差,根据偏移距旅行时差完成转换波数据的动校正。转换波零偏移距旅行时、转换波零偏移距速度、地层厚度等数据,可以通过纵波数据和转换波数据简单的分析计算获得,不需要复杂的计算过程。也不需要计算地震波的入射角等数据参数,减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。并且,采用本申请的方法进行转换波的动校正,可以校平近偏移距转换波数据,而不影响远偏移距转换波数据,具有更好的针对性。
在上述实施例的基础上,所述偏移距旅行时计算模块中预先建立的转换波旅行时模型,包括:
上式中,t表示偏移距x对应的转换波偏移距旅行时,c表示比例系数,0<c<1,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
转换波旅行时模型的建立方法可以同上述实施例一致,此处不再赘述。
本申请提供的近偏移距转换波动校正装置,预先建立的转换波旅行时模型,只需要转换波零偏移距速度、地层厚度、转换波零偏移距旅行时,不需要入射角等复杂参数。模型简单,参数的获取不需要复杂的计算过程,为后续转换波动校正提供了准确的数据基础。同时,减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上,所述偏移距旅行时计算模块具体用于:
S201、设置所述比例系数的初始值;
S202、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时;
S203、若判断获知所述初始转换波偏移距旅行时不满足预设条件,则调整所述比例系数,将所述调整后的比例系数的值作为所述初始值;
S204、重复步骤S202和S203,直至所述初始转换波偏移距旅行时满足所述预设条件,将所述初始转换波偏移距旅行时作为所述目标转换波偏移距旅行时。
本申请提供的近偏移距转换波动校正装置,可以通过简单的比例系数的调整,获得满足预设条件的初始转换波偏移距旅行时,进一步获得满足较平要求的转换波校正数据。提高了转换波动校正效率,减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上,所述偏移距旅行时计算模块具体用于:
预先设置比例系数的增量;
若判断获知所述初始转换波偏移距旅行时不满足预设条件,则将所述比例系数的初始值增加所述增量,将增加所述增量后的比例系数作为所述比例系数的初始值。
本申请提供的近偏移距转换波动校正装置,通过设置比例系数的增量,逐步调整比例系数的值,实现扫描式的计算,提高了转换波动校正效率,减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上,所述偏移距旅行时计算模块还用于:
若判断获知所述比例系数增加到1,所述初始转换波偏移距旅行时仍不满足所述预设条件,则调整所述增量的大小。
本申请提供的近偏移距转换波动校正装置,可以通过调整比例系数的增量,调整扫描的密度,增大获取达到较平要求的转换波校正数据。提高了转换波动校正效率和准确性,减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上,所述偏移距旅行时计算模块具体用于:
将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时后,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据;
若所述校正后的转换波数据小于校平程度阈值,则所述初始转换波偏移距旅行时不满足所述预设条件。
本申请提供的近偏移距转换波动校正装置,通过设置校平程度阈值,判断校正后的转换波数据是否达到较平要求,简单快速的完成转换波的动校正。减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
在上述实施例的基础上,所述零偏移距速度计算模块具体用于:
根据所述地层厚度和所述转换波零偏移距旅行时,利用v=2h/t0计算所述转换波零偏移距速度;
上式中,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
本申请提供的近偏移距转换波动校正装置,通过地层厚度和转换波零偏移距旅行时的简单计算,获得转换波零偏移距速度,为后续转换波动校正提供可准确的数据基础。减少了转换波动校正的计算量,简化了转换波动校正的计算过程,进一步提高了转换波动校正的效率和准确性。
需要说明的,上述所述的装置根据方法实施例的描述还可以包括其他的实施方式,具体的实现方式可以参照相关方法实施例的描述,在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本说明书实施例提供的上述近偏移距转换波动校正方法或装置可以在计算机中由处理器执行相应的程序指令来实现,如使用windows操作系统的c++语言在PC端实现、linux系统实现,或其他例如使用android、iOS系统程序设计语言在智能终端实现,以及基于量子计算机的处理逻辑实现等。本说明书提供的一种近偏移距转换波动校正装置的另一个实施例中,图8是本申请提供的另一种近偏移距转换波动校正装置实施例的模块结构示意图,如图8所示,本申请另一实施例提供的近偏移距转换波动校正装置可以包括处理器81以及用于存储处理器可执行指令的存储器82,
处理器81和存储器82通过总线83完成相互间的通信;
所述处理器81用于调用所述存储器82中的程序指令,以执行上述各近偏移距转换波动校正方法实施例所提供的方法,例如包括:根据采集到的地震数据,获取纵波数据和转换波数据;根据所述纵波数据计算地层厚度;根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时;根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度;根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差;根据所述偏移距旅行时差对所述转换波数据进行动校正。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式,具体的实现方式可以参照方法实施例的描述,在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例,仍然可以属于本说明书实施例的可选实施方案范围之内。
在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.一种近偏移距转换波动校正方法,其特征在于,所述方法包括:
根据采集到的地震数据,获取纵波数据和转换波数据;
根据所述纵波数据计算地层厚度;
根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时;
根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度;
根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;
根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差对所述转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据;
其中,所述预先建立的转换波旅行时模型,包括:
上式中,t表示偏移距x对应的目标转换波偏移距旅行时,c表示比例系数,0<c<1,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
2.如权利要求1所述的一种近偏移距转换波动校正方法,其特征在于,所述根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时,包括:
S201、设置所述比例系数的初始值;
S202、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时;
S203、若判断获知所述初始转换波偏移距旅行时不满足预设条件,则调整所述比例系数,将所述调整后的比例系数的值作为所述初始值;
S204、重复步骤S202和S203,直至所述初始转换波偏移距旅行时满足所述预设条件,将所述初始转换波偏移距旅行时作为所述目标转换波偏移距旅行时。
3.如权利要求2所述的一种近偏移距转换波动校正方法,其特征在于,所述调整所述比例系数,将所述调整后的比例系数的值作为所述初始值,包括:
预先设置比例系数的增量;
若判断获知所述初始转换波偏移距旅行时不满足预设条件,则将所述比例系数的初始值增加所述增量,将增加所述增量后的比例系数作为所述比例系数的初始值。
4.如权利要求3所述的一种近偏移距转换波动校正方法,其特征在于,所述方法还包括:
若判断获知所述比例系数增加到1,所述初始转换波偏移距旅行时仍不满足所述预设条件,则调整所述增量的大小。
5.如权利要求2-4任一项所述的一种近偏移距转换波动校正方法,其特征在于,所述初始转换波偏移距旅行时不满足预设条件,包括:
将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时后,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据;
若所述校正后的转换波数据小于校平程度阈值,则所述初始转换波偏移距旅行时不满足所述预设条件。
6.如权利要求1所述的一种近偏移距转换波动校正方法,其特征在于,所述根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度,包括:
根据所述地层厚度和所述转换波零偏移距旅行时,利用v=2h/t0计算所述转换波零偏移距速度;
上式中,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
7.一种近偏移距转换波动校正装置,其特征在于,包括:
数据采集模块,用于根据采集到的地震数据,获取纵波数据和转换波数据;
地层厚度计算模块,用于根据所述纵波数据计算地层厚度;
零偏移距旅行时计算模块,用于根据所述纵波数据和所述转换波数据计算转换波零偏移距旅行时;
零偏移距速度计算模块,用于根据所述地层厚度和所述转换波零偏移距旅行时,计算转换波零偏移距速度;
偏移距旅行时计算模块,用于根据所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度,利用预先建立的转换波旅行时模型,计算不同偏移距对应的目标转换波偏移距旅行时;
偏移距旅行时差计算模块,用于根据所述目标转换波偏移距旅行时和所述转换波零偏移距旅行时,计算不同偏移距对应的偏移距旅行时差;
校正模块,用于根据所述偏移距旅行时差对所述转换波数据进行动校正,以校平偏移距小于预设阈值的转换波数据;
其中,所述偏移距旅行时计算模块中预先建立的转换波旅行时模型,包括:
上式中,t表示偏移距x对应的转换波偏移距旅行时,c表示比例系数,0<c<1,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
8.如权利要求7所述的一种近偏移距转换波动校正装置,其特征在于,所述偏移距旅行时计算模块具体用于:
S201、设置所述比例系数的初始值;
S202、将所述比例系数的初始值、所述地层厚度、所述转换波零偏移距旅行时、所述转换波零偏移距速度代入所述转换波旅行时模型,计算获得所述比例系数初始值对应的不同偏移距的初始转换波偏移距旅行时;
S203、若判断获知所述初始转换波偏移距旅行时不满足预设条件,则调整所述比例系数,将所述调整后的比例系数的值作为所述初始值;
S204、重复步骤S202和S203,直至所述初始转换波偏移距旅行时满足所述预设条件,将所述初始转换波偏移距旅行时作为所述目标转换波偏移距旅行时。
9.如权利要求8所述的一种近偏移距转换波动校正装置,其特征在于,所述偏移距旅行时计算模块具体用于:
预先设置比例系数的增量;
若判断获知所述初始转换波偏移距旅行时不满足预设条件,则将所述比例系数的初始值增加所述增量,将增加所述增量后的比例系数作为所述比例系数的初始值。
10.如权利要求9所述的一种近偏移距转换波动校正装置,其特征在于,所述偏移距旅行时计算模块还用于:
若判断获知所述比例系数增加到1,所述初始转换波偏移距旅行时仍不满足所述预设条件,则调整所述增量的大小。
11.如权利要求8-10任一项所述的一种近偏移距转换波动校正装置,其特征在于,所述偏移距旅行时计算模块具体用于:
将所述初始转换波偏移距旅行时作为所述目标偏移距旅行时后,根据所述目标转换波偏移距旅行时和所述转换波数据,计算不同偏移距对应的偏移距旅行时差;
根据所述偏移距旅行时差进行转换波动校正,获得校正后的转换波数据;
若所述校正后的转换波数据小于校平程度阈值,则所述初始转换波偏移距旅行时不满足所述预设条件。
12.如权利要求7所述的一种近偏移距转换波动校正装置,其特征在于,所述零偏移距速度计算模块具体用于:
根据所述地层厚度和所述转换波零偏移距旅行时,利用v=2h/t0计算所述转换波零偏移距速度;
上式中,v表示所述转换波零偏移距速度,h表示所述地层厚度,t0表示所述转换波零偏移距旅行时。
13.一种近偏移距转换波动校正装置,其特征在于,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现如权利要求1至6中任意一项所述方法的步骤。
CN201711210811.4A 2017-11-28 2017-11-28 一种近偏移距转换波动校正方法及装置 Active CN108181656B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711210811.4A CN108181656B (zh) 2017-11-28 2017-11-28 一种近偏移距转换波动校正方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711210811.4A CN108181656B (zh) 2017-11-28 2017-11-28 一种近偏移距转换波动校正方法及装置

Publications (2)

Publication Number Publication Date
CN108181656A CN108181656A (zh) 2018-06-19
CN108181656B true CN108181656B (zh) 2019-08-06

Family

ID=62545390

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711210811.4A Active CN108181656B (zh) 2017-11-28 2017-11-28 一种近偏移距转换波动校正方法及装置

Country Status (1)

Country Link
CN (1) CN108181656B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109725354B (zh) * 2018-11-20 2020-07-10 中国石油天然气集团有限公司 各向异性速度建模方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4881209A (en) * 1988-10-12 1989-11-14 Mobil Oil Corporation Method of processing converted mode seismic data
CN101251604B (zh) * 2008-04-12 2011-02-16 中国石油集团西北地质研究所 二参数转换波速度分析及动校正方法
US9753166B2 (en) * 2014-05-20 2017-09-05 Cgg Services Sas P-wave and S-wave separation of seismic data in the presence of statics and irregular geometry
CN104216016B (zh) * 2014-08-12 2017-01-25 中国石油天然气集团公司 一种纵波约束扫描的转换波动校正方法以及系统
CN104155691B (zh) * 2014-08-27 2017-12-26 中国石油天然气集团公司 转换波各向异性速度分析方法及装置
CN106443777B (zh) * 2016-08-15 2018-02-23 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波各向异性参数谱的获取方法及转换波速度分析方法

Also Published As

Publication number Publication date
CN108181656A (zh) 2018-06-19

Similar Documents

Publication Publication Date Title
CN107765298B (zh) 一种确定近地表速度模型的方法及装置
CN108181649B (zh) 一种介质分界面振幅补偿方法及装置
US10288752B2 (en) Automatic focus identification method and system for Karst cave reservoir
CN102129084A (zh) 一种井控获取地震薄储层速度的方法及装置
CN107797154B (zh) 一种电成像测井图像刻度方法及装置
CN109143358B (zh) 一种获得深层、超深层碎屑岩地层压力结构的方法及装置
CN108828668A (zh) 一种叠前时间偏移数据处理方法及装置
CN107817525B (zh) 一种近地表常速填充方法及装置
CN108983285A (zh) 一种基于矩张量的多种震源波场模拟方法及装置
CN108287365B (zh) 一种基于波动方程的vsp三参数同步反演方法及装置
CN108107471B (zh) 一种分方位初至数据体的获取方法及装置
CN108089238A (zh) 一种利用综合地球物理资料确定深层裂谷的方法及装置
CN109725349A (zh) 一种恢复生烃关键期古地貌的方法及装置
CN104199103A (zh) 转换波静校正方法及装置
CN105549084A (zh) 一种三维高精度速度建模方法和系统
CN108710149A (zh) 一种薄储层的预测方法及装置
CN108181656B (zh) 一种近偏移距转换波动校正方法及装置
CN107179546A (zh) 一种确定速度场的方法及装置
CN108828669A (zh) 一种二维相交测线静校正处理方法、装置及系统
CN108469633A (zh) 一种地层品质因子的计算方法及装置
CN108508481B (zh) 一种纵波转换波地震数据时间匹配的方法、装置及系统
CN108693557A (zh) 一种地震数据处理方法及装置
CN106772596B (zh) 一种确定叠前时间偏移速度场的方法及装置
CN107976716B (zh) 一种远偏移距转换波动校正方法及装置
CN111208558B (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