CN111751881A - 一种海上采集地震数据旅行时的校正方法、装置及*** - Google Patents

一种海上采集地震数据旅行时的校正方法、装置及*** Download PDF

Info

Publication number
CN111751881A
CN111751881A CN201910247385.4A CN201910247385A CN111751881A CN 111751881 A CN111751881 A CN 111751881A CN 201910247385 A CN201910247385 A CN 201910247385A CN 111751881 A CN111751881 A CN 111751881A
Authority
CN
China
Prior art keywords
seismic data
travel time
travel
correction
motion
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.)
Pending
Application number
CN201910247385.4A
Other languages
English (en)
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 CN201910247385.4A priority Critical patent/CN111751881A/zh
Publication of CN111751881A publication Critical patent/CN111751881A/zh
Pending legal-status Critical Current

Links

Images

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
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/52Move-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

一种海上采集地震数据旅行时的校正方法、装置及***
技术领域
本发明涉及地球物理勘探技术领域,特别涉及一种海上采集地震数据旅行时的校正方法、装置及***。
背景技术
在地震数据常规处理中,无论是叠加还是速度分析,都要求检波器位置固定。但是海上利用拖缆采集数据时,当震源激发后,在数据采集过程中,检波器不是固定的,而是随着拖缆一直在运动,这种运动就会引起地震记录旅行时的误差,在处理中需要把这种误差校正掉。然而,由于检波器一直在运动,造成的时间误差会随着采集数据旅行时的变化而变化,使得检波器运动校正是个很复杂的过程;此外,由于采集时拖缆的运动速度相对地震波在水中和地下固体介质中的传播速度小得多,使得这种时间误差相对较小,因此在以往的拖缆数据处理中未考虑这个问题。近年来,随着海上数据宽频处理技术的发展,对数据处理和解释的精度要求越来越高,迫切需要一种方法来消除海上拖缆地震数据的检波器运动效应。
因此,业内亟需一种可以消除海上拖缆采集地震数据时检波器运动造成旅行时误差的解决方案。
发明内容
本说明书实施例目的在于提供一种海上采集地震数据旅行时的校正方法、装置及***,可以消除海上拖缆采集地震数据时检波器运动造成的旅行时误差,为后续地震处理提供更加符合要求的数据,从而提高成像质量。
一方面本申请提供了一种海上采集地震数据旅行时的校正方法,包括:
获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
基于所述地震数据,计算所述地震数据中采样点的旅行时;
根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
本说明书提供的所述方法的另一个实施例中,还包括:
根据所述校正后地震数据的旅行时,获得所述校正后地震数据的旅行时所对应采样点的样点值。
本说明书提供的所述方法的另一个实施例中,按照下面公式计算所述地震数据中采样点的旅行时t0
t0=(n-1)×dt
其中,t0表示地震数据中采样点的旅行时,n表示第n个采样点,dt表示采样间隔。
本说明书提供的所述方法的另一个实施例中,所述地震数据对应的共中心点是检波器到震源的中点位置,所述动校正速度是通过速度分析法处理得到。
本说明书提供的所述方法的另一个实施例中,按照下面公式计算所述检波器运动校正时差Δt:
Figure BDA0002011426990000021
其中,Δt表示检波器运动校正时差,vB表示拖缆移动速度,x表示炮检距,v表示动校正速度。
本说明书提供的所述方法的另一个实施例中,按照下面公式获得所述校正后地震数据的旅行时t:
t=t0+Δt
其中,t表示校正后地震数据的旅行时,t0表示地震数据中采样点的旅行时,Δt表示检波器运动校正时差。
另一方面,本说明书实施例还提供一种海上采集地震数据旅行时的校正装置,所述装置包括:
数据信息获取模块,用于获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
旅行时计算模块,用于基于所述地震数据,计算所述地震数据中采样点的旅行时;
动校正速度获取模块,用于根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
校正时差计算模块,用于从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
校正后旅行时获得模块,用于利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
本说明书提供的所述装置的另一个实施例中,所述装置还包括:
校正后样点值获得模块,用于根据所述校正后地震数据的旅行时,获得所述校正后地震数据的旅行时所对应采样点的样点值。
另一方面,本说明书实施例还提供一种海上采集地震数据旅行时的校正设备,包括处理器及用于存储处理器可执行指令的存储器,所述指令被所述处理器执行时实现包括以下步骤:
获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
基于所述地震数据,计算所述地震数据中采样点的旅行时;
根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
另一方面,本说明书实施例还提供一种海上采集地震数据旅行时的校正***,包括至少一个处理器以及存储计算机可执行指令的存储器,所述处理器执行所述指令时实现上述任意一个实施例所述方法的步骤。
本说明书实施例提供的一种海上采集地震数据旅行时的校正方法、装置及***,可以通过从海上地震勘探地震数据和海上拖缆数据采集时的班报信息中获取拖缆移动速度、炮检距以及采样点的旅行时,利用地震数据对应的共中心点以及速度分析法获得旅行时对应的动校正速度,来计算由于检波器运动造成的旅行时误差,进而实现对海上拖缆采集地震数据时检波器运动造成的旅行时误差的校正,为后续地震处理提供更加符合要求的数据,从而提高成像质量。
附图说明
为了更清楚地说明本说明书实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本说明书提供的海上拖缆采集地震数据时一个实施例中检波器运动效应示意图;
图2是本说明书提供的一种海上采集地震数据旅行时的校正方法一个实施例的流程示意图;
图3是本说明书提供的一个实施例中检波器运动校正前CMP道集示意图;
图4是本说明书提供的一个实施例中检波器运动校正后CMP道集示意图;
图5是本说明书提供的一个实施例中检波器运动校正前CMP道集中部分放大显示的示意图;
图6是本说明书提供的一个实施例中检波器运动校正后CMP道集中部分放大显示的示意图;
图7是本说明书提供的一种海上采集地震数据旅行时的校正装置一个实施例的模块结构示意图;
图8是本说明书提供的一种海上采集地震数据旅行时的校正***一个实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本说明书中的一部分实施例,而不是全部的实施例。基于本说明书中的一个或多个实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本说明书实施例保护的范围。
目前,海上地震数据采集普遍采用拖缆的方式进行。这种海上地震数据采集方式,一般由采集船拖曳着震源和多条采集电缆组成采集***,震源相隔一定时间放炮,震源激发的同时接收电缆开始记录地下反射信号,并由导航定位***记录放炮时间和对应的炮点、接收点的位置信息等。采集过程中为了获得理想的观测***参数,采集船必须保持一定的航速使接收电缆处于拉直的状态。因此,与陆上地震数据采集时检波器放置在固定位置不同,海上拖缆数据采集时,接收电缆是随采集船移动的,即检波器不是固定的。然而,现有的勘探地球物理地震资料处理技术都是从陆上地震数据处理技术发展起来的,认为在地震资料采集的过程中检波器固定不动,这样,当处理海上拖缆采集数据时,就会引起地震记录旅行时的误差,随着检波器的运动,造成的时间误差会随着采集数据旅行时的变化而变化,使得检波器运动校正是个很复杂的过程,而且随着海上数据宽频处理技术的发展,对数据处理和解释的精度要求越来越高,迫切需要一种方法来消除海上拖缆地震数据的检波器运动效应,为后续的地震资料处理、解释提供高品质的基础数据,提高成像质量以及油气藏变化监测的精度。
为此,本申请实施例提供了一种海上采集地震数据旅行时的校正方法,根据拖缆移动速度、炮检距以及与采样点旅行时对应的动校正速度,求取检波器运动造成的旅行时误差,进而消除检波器运动效应,为后续地震处理提供更加符合地震波实际传播路径的数据,从而提高成像质量。具体地,如图1所示,图1是本说明书提供的海上拖缆采集地震数据时一个实施例中检波器运动效应示意图。其中,震源放炮时,检波器位于实心圆圈处,由于拖缆的运动,在数据采集的某一时刻,检波器移动到空心圆圈处,相应波场的射线路径从实线变成了虚线,检波器运动校正就是通过时差校正,把射线路经从虚线校正到实线,从而消除检波器运动的影响。
下面以一个具体的应用场景为例对本说明书实施方案进行说明。具体的,图2是本说明书提供的一种海上采集地震数据旅行时的校正方法一个实施例的流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构,但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中,这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时,可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
当然,下述实施例的描述并不对基于本说明书的其他可扩展到的技术方案构成限制。
具体的一种实施例如图2所示,本说明书提供的一种海上采集地震数据旅行时的校正方法一个实施例中,所述方法可以包括:
S1、获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度。
班报信息是指海上拖缆采集时,震源(如:***)组记录的相关施工信息,主要有船速,船的位置(如:经、纬度)、水深、测线的起始和终止经纬度等。通常情况,拖缆移动速度、检波器的速度、船速是相等的。所以,在一些实施例中,获得了船速就可以获得拖缆移动速度。
本说明书一个或多个实施例中,通过采集海上地震勘探数据,取得叠前单炮地震数据。所述采集海上地震勘探数据主要是通过采集船来完成的,具体的,采集海上地震勘探数据是把地震仪器安装在采集船上,使用海上专用的电缆和检波器,在采集船航行中连续地对震源发出的地震波进行激发和接收。所述叠前单炮地震数据是利用一个炮点激发,多个检波器进行采集所获得的多道地震数据,例如,利用一个炮点激发,三个检波器对地震数据进行采集,所获得的三道数据可以认为是叠前单炮地震数据。
地震数据有多种保存格式,比较常用的是SEG-Y格式,它是由SEG(Society ofExploration Geophysicists,勘探地球物理学家协会)提出的标准磁带数据格式之一,可以由三部分组成:400字节的二进制文件头、每个地震道的240字节道头和每个地震道的地震数据,它的数据都可以按地震道的顺序一个个存放。
本说明书实施例中,从海上拖缆数据采集时的班报信息中,读取船速就相当于读取拖缆的移动速度。
S2、基于所述地震数据,计算所述地震数据中采样点的旅行时。
采样点是一定时间(称采样周期)采集的数据点。采样点的旅行时是指炮点发射的地震波到被检波点检测到所用的时长。假设对一个模拟信号x(t)每隔Δt时间采样一次,则t=0、Δt、2Δt、3Δt……,时间间隔Δt被称为采样间隔或采样周期,x(t)的数值就被称为采样值,即采样点对应的数值,所有x(0)、xΔt)、x(2Δt)都是采样值。其中,获得的地震数据中可以包括多道地震数据,每道数据包括多个采样点。
本说明书一个或多个实施例中,在获得地震数据后,可以根据地震数据计算所述地震数据中采样点的旅行时。具体地,可以从所述地震数据中每个地震道的240字节的道头中提取所述地震道的采样点数、采样间隔(以微秒表示),按照下面公式计算所述地震数据中采样点的旅行时t0
t0=(n-1)×dt (1)
其中,t0可以表示地震数据中采样点的旅行时,n可以表示第n个采样点,dt可以表示采样间隔。
需要说明的是,采样点的旅行时还可以通过其他方式获得,例如由于仪器或者其他因素导致地震数据开始记录的时候存在延迟时,这时计算采样点的旅行时要去除延迟时的影响等等,本说明书对此不做限定。
进一步地,由于地震数据中可以包括多道地震数据,每道数据包括多个采样点,所以可以根据上述公式获得一道地震数据以及多道地震数据中所有采样点的旅行时,也就是说利用公式(1)可以获得不同地震道数据中不同采样点的旅行时,具体可以参见上述实施例的方法,对此不再赘述。
S3、根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度。
所述地震数据对应的共中心点(Common Middle Point,CMP)是检波器到震源的中点位置,一道地震数据对应的CMP是相同的。动校正是地震数据处理的重要组成部分,动校正是将非零炮检距上所接收到的反射波旅行时校正到零炮检距上的自激自收反射波旅行时,进行动校正的目的是消除炮检距对反射波旅行时的影响,其中,炮检距指检波器到震源的距离。对每一个CMP来说,动校正速度是按照时间排列的一系列速度值。同一CMP,动校正速度与时间是对应关系,因此,在获得CMP和对应的旅行时时,就可以从速度分析法处理得到的动校正速度文件中读取对应的动校正速度,即所述动校正速度是通过速度分析法处理得到。
需要说明的是,在一般情况下,同一CMP,随着时间的增大,动校正速度的数值也会随着增大,即动校正速度是变化的,也就是说,不同旅行时对应的动校正速度一般是不同的,但也可能存在相同的情况。
本说明书一个或多个实施例中,通过计算一道地震数据中一个采样点的旅行时,然后根据该道地震数据对应的CMP获取所述采样点旅行时对应的动校正速度,相应地,可以通过类似的方法获取另一采样点旅行时对应的动校正速度以及地震数据中所有采样点对应的动校正速度,对此不再赘述。
S4、从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差。
炮检距指检波器到震源的距离,其可以从每个地震道的240字节的道头中获得。运动校正时差也叫正常时差校正,其是用来消除接收点与激发点不在同一位置而造成的地震波到达接收点时,旅行时间与法线反射时间之间的差值。
具体地,由于海上利用拖缆采集数据时,检波器不是固定的,而是随着拖缆一直在运动,所以拖缆采集数据的旅行时与检波器固定时采集数据的旅行时之间是存在差异的,所以,检波器运动校正时差就是假设检波器固定不动时接收到的地震波旅行时与检波器随着拖缆移动后接收到地震波旅行时之间的差值。
本说明书一个或多个实施例中,可以按照下面公式计算所述检波器运动校正时差Δt:
Figure BDA0002011426990000071
其中,Δt可以表示检波器运动校正时差,vB可以表示拖缆移动速度,x可以表示炮检距,v可以表示动校正速度。
进一步地,由于同一道地震数据对应的炮检距、拖缆移动速度相同,不同采样点旅行时对应的动校正速度不同,因此,同一道地震数据中,不同采样点的检波器运动校正时差不相同,即在同一地震道中,下一采样点的检波器运动校正时差与上一采样点的检波器运动校正时差不同。此外,由于不同地震道数据对应的拖缆移动速度相同,炮检距和采样点旅行时对应的动校正速度不同,因此,不同地震道地震数据中,采样点的检波器运动校正时差也不相同。所以,地震数据中需要根据拖缆移动速度、炮检距、动校正速度,计算每一个采样点的检波器运动校正时差。
本说明书一个实施例中,按照公式(2)可以获得同一道地震数据中不同采样点的检波器运动校正时差,相应地,可以按照相同方法获得所有地震数据采样点对应的检波器运动校正时差,对此不再赘述。
S5、利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
校正后采样点的旅行时是指在震源放炮时检波器所在位置接收到地震波的时间,也就是常规地震数据处理中检波器位置固定不变时接收到地震波的时间。
本说明书一个或多个实施例中,可以按照下面公式获得所述校正后地震数据的旅行时t:
t=t0+Δt (3)
其中,t可以表示校正后地震数据的旅行时,t0可以表示地震数据中采样点的旅行时,Δt可以表示检波器运动校正时差。
本说明书另一个实施例中,在获得校正后地震数据的旅行时后,还包括根据所述校正后地震数据的旅行时,获得所述校正后地震数据的旅行时所对应采样点的样点值。具体地,可以根据所述校正后地震数据的旅行时,采用插值方法获得所述校正后地震数据的旅行时所对应采样点的样点值。其中,插值方法可以是sinc插值法或者拉格朗日插值方法,也可以是本领域人员所了解的其他插值方法,本说明书对此不做限定。
进一步地,按照相同的方法对同一道地震数据中不同采样点进行处理,可以获得一道校正后地震数据的旅行时所对应采样点的样点值以及所有地震数据的旅行时所对应采样点的样点值,为后续地震处理提供更加符合要求的数据,从而提高成像质量。
本申请实施例提供的一种海上采集地震数据旅行时的校正方法,通过从海上地震勘探地震数据和海上拖缆数据采集时的班报信息中获取拖缆移动速度、炮检距以及采样点的旅行时,并利用地震数据对应的共中心点以及速度分析法获得旅行时对应的动校正速度,来计算由于检波器运动造成的旅行时误差,进而实现对海上拖缆采集地震数据时检波器运动造成的旅行时误差的校正,为后续地震处理提供更加符合要求的数据,从而提高成像质量。
为了进一步来验证本方法的实用性和可行性,本说明书还提供了应用上述方案的一个具体实例,具体过程和步骤如下:
(1)采集海上地震勘探数据,取得叠前单炮地震数据。
(2)从海上拖缆数据采集时的班报信息中,读取拖缆的移动速度vB
(3)从所述地震数据的道头中,读取一道数据的炮检距x。
(4)按照公式t0=(n-1)×dt,计算所述一道地震数据中一个采样点的旅行时t0
(5)根据所述一道地震数据对应的CMP点,读取所述一个采样点的旅行时t0对应的动校正速度v。
(6)按照公式
Figure BDA0002011426990000091
计算检波器运动校正时差Δt。
(7)按照公式t=t0+Δt,计算检波器运动校正后的地震数据旅行时t。
(8)根据所述检波器运动校正后的地震数据旅行时t,采用sinc插值方法计算所述检波器运动校正后的地震数据旅行时t对应采样点的样点值。
(9)对一道地震数据上的所有采样点重复执行步骤(2)到步骤(8),完成该道数据的检波器运动校正。
(10)对地震数据上的所有道数据,都重复执行步骤(2)到步骤(9),完成所有数据的检波器运动校正。
需要说明的是,上述各实施例中的公式仅仅是一种示意性的表达公式,可以根据实际需要对各个公式进行调整和变形,本申请实施例不作具体的限定。
结果表明,利用本申请实施例提供的方案,可以有效可以消除海上拖缆采集地震数据时检波器运动造成的旅行时误差。如图3、图4、图5、图6所示,图3是本说明书提供的一个实施例中检波器运动校正前CMP道集示意图,图4是本说明书提供的一个实施例中检波器运动校正后CMP道集示意图,图5是本说明书提供的一个实施例中检波器运动校正前CMP道集中部分放大显示的示意图,图6是本说明书提供的一个实施例中检波器运动校正后CMP道集中部分放大显示的示意图。其中,纵坐标表示时间,单位为毫秒,横坐标表示跑检距,单位为米。从检波器运动校正前后CMP道集中部分放大显示图可以得出,检波器运动校正后的旅行时增加,而且检波器运动造成的旅行时误差被弥补,所以利用本申请实施例提供的方案可以有效消除检波器运动效应导致的误差。
基于上述所述的一种海上采集地震数据旅行时的校正方法,本说明书一个或多个实施例还提供一种海上采集地震数据旅行时的校正装置。所述的装置可以包括使用了本说明书实施例所述方法的***(包括分布式***)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思,本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似,因此本说明书实施例具体的装置的实施可以参见前述方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
具体地,图7是本说明书提供的一种海上采集地震数据旅行时的校正装置一个实施例的模块结构示意图,如图7所示,本说明书提供的海上拖缆采集地震数据旅行时的校正装置包括:数据信息获取模块121,旅行时计算模块122,动校正速度获取模块123,校正时差计算模块124、校正后旅行时获得模块125。其中,数据信息获取模块121,可以用于获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
旅行时计算模块122,可以用于基于所述地震数据,计算所述地震数据中采样点的旅行时;
动校正速度获取模块123,可以用于根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
校正时差计算模块124,可以用于从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
校正后旅行时获得模块125,可以用于利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
其中,所述一种海上采集地震数据旅行时的校正装置,还可以包括:
校正后样点值获得模块,可以用于根据所述校正后地震数据的旅行时,获得所述校正后地震数据的旅行时所对应采样点的样点值。
本说明书提供的一种海上采集地震数据旅行时的校正装置,通过从海上地震勘探地震数据和海上拖缆数据采集时的班报信息中获取拖缆移动速度、炮检距以及采样点的旅行时,并利用地震数据对应的共中心点以及速度分析法获得旅行时对应的动校正速度,来计算由于检波器运动造成的旅行时误差,进而实现对海上拖缆采集地震数据时检波器运动造成的旅行时误差的校正,为后续地震处理提供更加符合要求的数据,从而提高成像质量。
需要说明的,上述所述的装置根据方法实施例的描述还可以包括其他的实施方式,具体的实现方式可以参照相关方法实施例的描述,在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上,所述的存储介质可以计算机读取并执行,实现本说明书实施例所描述方案的效果。因此,本说明书还提供一种海上采集地震数据旅行时的校正设备,包括处理器及用于存储处理器可执行指令的存储器,所述指令被所述处理器执行时实现包括以下步骤:
获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
基于所述地震数据,计算所述地震数据中采样点的旅行时;
根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
所述存储介质可以包括用于存储信息的物理装置,通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括:利用电能方式存储信息的装置如,各式存储器,如RAM、ROM等;利用磁能方式存储信息的装置如,硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘;利用光学方式存储信息的装置如,CD或DVD。当然,还有其他方式的可读存储介质,例如量子存储器、石墨烯存储器等等。
需要说明的,上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述,在此不作一一赘述。
本说明书实施例提供的上述海上拖缆采集地震数据旅行时的校正方法、装置、设备可以在计算机中由处理器执行相应的程序指令来实现,如使用windows操作***的c++语言在PC端实现、linux***实现,或其他例如使用android、iOS***程序设计语言在智能终端实现,以及基于量子计算机的处理逻辑实现等。本说明书提供的一种海上采集地震数据旅行时的校正***的一个实施例中,图8是本说明书提供的一种海上采集地震数据旅行时的校正***一个实施例的模块结构示意图,如图8所示,本说明书提供的海上拖缆采集地震数据旅行时的校正***可以包括处理器131以及用于存储处理器可执行指令的存储器132,处理器131和存储器132通过总线133完成相互间的通信;
所述处理器131用于调用所述存储器132中的程序指令,以执行上述海上拖缆采集地震数据旅行时的校正方法实施例所提供的方法,例如包括:获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;基于所述地震数据,计算所述地震数据中采样点的旅行时;根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式,具体的实现方式可以参照方法实施例的描述,在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例,仍然可以属于本说明书实施例的可选实施方案范围之内。
在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.一种海上采集地震数据旅行时的校正方法,其特征在于,包括:
获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
基于所述地震数据,计算所述地震数据中采样点的旅行时;
根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
2.如权利要求1所述的一种海上采集地震数据旅行时的校正方法,其特征在于,还包括:
根据所述校正后地震数据的旅行时,获得所述校正后地震数据的旅行时所对应采样点的样点值。
3.如权利要求1所述的一种海上采集地震数据旅行时的校正方法,其特征在于,按照下面公式计算所述地震数据中采样点的旅行时t0
t0=(n-1)×dt
其中,t0表示地震数据中采样点的旅行时,n表示第n个采样点,dt表示采样间隔。
4.如权利要求1所述的一种海上采集地震数据旅行时的校正方法,其特征在于,所述地震数据对应的共中心点是检波器到震源的中点位置,所述动校正速度是通过速度分析法处理得到。
5.如权利要求1所述的一种海上采集地震数据旅行时的校正方法,其特征在于,按照下面公式计算所述检波器运动校正时差Δt:
Figure FDA0002011426980000011
其中,Δt表示检波器运动校正时差,vB表示拖缆移动速度,x表示炮检距,v表示动校正速度。
6.如权利要求1所述的一种海上采集地震数据旅行时的校正方法,其特征在于,按照下面公式获得所述校正后地震数据的旅行时t:
t=t0+Δt
其中,t表示校正后地震数据的旅行时,t0表示地震数据中采样点的旅行时,Δt表示检波器运动校正时差。
7.一种海上采集地震数据旅行时的校正装置,其特征在于,包括:
数据信息获取模块,用于获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
旅行时计算模块,用于基于所述地震数据,计算所述地震数据中采样点的旅行时;
动校正速度获取模块,用于根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
校正时差计算模块,用于从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
校正后旅行时获得模块,用于利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
8.如权利要求7所述的一种海上采集地震数据旅行时的校正装置,其特征在于,还包括:
校正后样点值获得模块,用于根据所述校正后地震数据的旅行时,获得所述校正后地震数据的旅行时所对应采样点的样点值。
9.一种海上采集地震数据旅行时的校正设备,其特征在于,包括处理器及用于存储处理器可执行指令的存储器,所述指令被所述处理器执行时实现包括以下步骤:
获取地震数据和班报信息,所述班报信息中至少包括拖缆移动速度;
基于所述地震数据,计算所述地震数据中采样点的旅行时;
根据所述地震数据对应的共中心点,获取所述旅行时对应的动校正速度;
从所述地震数据的道头中获取所述地震数据的炮检距,根据所述拖缆移动速度、所述炮检距、所述动校正速度,计算检波器运动校正时差;
利用所述旅行时与所述检波器运动校正时差,获得校正后地震数据的旅行时。
10.一种海上采集地震数据旅行时的校正***,其特征在于,包括至少一个处理器以及存储计算机可执行指令的存储器,所述处理器执行所述指令时实现权利要求1-6中任意一项所述方法的步骤。
CN201910247385.4A 2019-03-29 2019-03-29 一种海上采集地震数据旅行时的校正方法、装置及*** Pending CN111751881A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910247385.4A CN111751881A (zh) 2019-03-29 2019-03-29 一种海上采集地震数据旅行时的校正方法、装置及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910247385.4A CN111751881A (zh) 2019-03-29 2019-03-29 一种海上采集地震数据旅行时的校正方法、装置及***

Publications (1)

Publication Number Publication Date
CN111751881A true CN111751881A (zh) 2020-10-09

Family

ID=72671689

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910247385.4A Pending CN111751881A (zh) 2019-03-29 2019-03-29 一种海上采集地震数据旅行时的校正方法、装置及***

Country Status (1)

Country Link
CN (1) CN111751881A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114460646A (zh) * 2022-04-13 2022-05-10 山东省科学院海洋仪器仪表研究所 一种基于波场激发近似的反射波旅行时反演方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020126576A1 (en) * 2001-03-07 2002-09-12 Westerngeco L.L.C. Seismic receiver motion compensation
CN101630019A (zh) * 2008-07-16 2010-01-20 中国石油集团东方地球物理勘探有限责任公司 一种校正预测的海底多次波大时差的方法
CN105259570A (zh) * 2015-11-12 2016-01-20 中国海洋石油总公司 一种海上时移地震接收点位移校正方法
CN106855640A (zh) * 2017-02-03 2017-06-16 中国石油天然气集团公司 一种海底电缆地震数据初至计算方法及装置
CN107576986A (zh) * 2017-08-30 2018-01-12 中国石油集团东方地球物理勘探有限责任公司 一种确定海底地震反射波旅行时的方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020126576A1 (en) * 2001-03-07 2002-09-12 Westerngeco L.L.C. Seismic receiver motion compensation
CN101630019A (zh) * 2008-07-16 2010-01-20 中国石油集团东方地球物理勘探有限责任公司 一种校正预测的海底多次波大时差的方法
CN105259570A (zh) * 2015-11-12 2016-01-20 中国海洋石油总公司 一种海上时移地震接收点位移校正方法
CN106855640A (zh) * 2017-02-03 2017-06-16 中国石油天然气集团公司 一种海底电缆地震数据初至计算方法及装置
CN107576986A (zh) * 2017-08-30 2018-01-12 中国石油集团东方地球物理勘探有限责任公司 一种确定海底地震反射波旅行时的方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
O.R.伯格,D.G.伍尔弗顿: "数字信号处理与VXI自动化测试技术", 中国石油大学出版社, pages: 173 - 174 *
施剑等: "南黄海盆地地震勘探中的非双曲线时差分析", 《石油地球物理勘探》, no. 02, 15 April 2013 (2013-04-15), pages 192 - 199 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114460646A (zh) * 2022-04-13 2022-05-10 山东省科学院海洋仪器仪表研究所 一种基于波场激发近似的反射波旅行时反演方法
CN114460646B (zh) * 2022-04-13 2022-06-28 山东省科学院海洋仪器仪表研究所 一种基于波场激发近似的反射波旅行时反演方法

Similar Documents

Publication Publication Date Title
EP3120167B1 (en) Reconstructing impulsive source seismic data from time distributed firing airgun array data
Baker et al. Real-time earthquake location using Kirchhoff reconstruction
US9030910B2 (en) Jointly interpolating and deghosting seismic data
CN108828668B (zh) 一种叠前时间偏移数据处理方法及装置
US20100211320A1 (en) Reconstructing a seismic wavefield
US9405028B2 (en) Method and apparatus for multi-component datuming
CN107193045B (zh) 一种地震数据处理方法及装置
CN111239808B (zh) 一种地震反演谱分解方法、装置、电子设备及存储介质
AU2009212892B2 (en) Interbed multiple prediction
CN104199083A (zh) 一种确定宽线地震观测***的方法和装置
Özdemir et al. On data-independent multicomponent interpolators and the use of priors for optimal reconstruction and 3D up/down separation of pressure wavefields
CN108897036B (zh) 一种地震数据处理方法及装置
CN108181644B (zh) 一种采样方法及装置
Edwards et al. Effect of interarray elevation differences on infrasound beamforming
CN109696705A (zh) 深海obs上下行波场成像数据校正方法及装置
CN111751881A (zh) 一种海上采集地震数据旅行时的校正方法、装置及***
Trabattoni et al. From strain to displacement: using deformation to enhance distributed acoustic sensing applications
CN104459774A (zh) 一种地质岩性差异识别方法及***
CN111435171B (zh) 海底节点地震数据的时钟漂移校正方法、装置及存储介质
CN112764102A (zh) 地震资料采集方法及装置
CN103617332A (zh) 一种测量卫星导航算法性能的仿真平台和方法
GB2425355A (en) Residual move out processing of seismic data
AU2016370548B2 (en) Method for removing interference caused by time overlapping seismic recordings and seismic survey acquisition method associated therewith
CN109471173B (zh) 一种剩余静校正方法及装置
US20170176615A1 (en) Method for processing seismic images

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20201009