CN103852791B - 用于在多勘探船地震勘探系统中管理爆破的方法 - Google Patents
用于在多勘探船地震勘探系统中管理爆破的方法 Download PDFInfo
- Publication number
- CN103852791B CN103852791B CN201310627253.7A CN201310627253A CN103852791B CN 103852791 B CN103852791 B CN 103852791B CN 201310627253 A CN201310627253 A CN 201310627253A CN 103852791 B CN103852791 B CN 103852791B
- Authority
- CN
- China
- Prior art keywords
- explosion
- exploration ship
- subordinate
- exploration
- ship
- 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.)
- Expired - Fee Related
Links
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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Oceanography (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种管理爆破的方法,其包括针对每个从属爆破勘探船:作为调度爆破勘探船及与调度爆破勘探船相关联的爆破点的速度值的函数,计算一系列第一理论爆破时间,每个第一理论爆破时间与调度爆破勘探船的接下来一个爆破相关联;作为从属爆破勘探船及与从属爆破勘探船相关联的爆破点的速度值的函数,计算一系列第二理论爆破时间,每个第二理论爆破时间与从属爆破勘探船的接下来一个爆破相关联;作为所述系列的第一理论爆破时间的函数,计算一系列内插虚拟爆破时间,至少包括与就在从属爆破勘探船的接下来的爆破之前及之后的爆破相关联的内插虚拟爆破时间。
Description
1.技术领域
本发明的领域是地球物理学数据采集领域。
更准确地说,本发明涉及海上地震勘探采集,包括几个勘探船以实现海床地壳的成像。
具体来说,本发明涉及一种在包括几个爆破勘探船(一个调度爆破勘探船及至少一个从属爆破勘探船)及至少一个监听勘探船的多勘探船地震勘探系统中管理爆破的方法。每个爆破勘探船包括至少一个在确定的瞬间执行爆破的震源。每个监听勘探船具有一个记录系统并且拖拽着至少一个集成有地震传感器的拖缆。
实际上,多勘探船地震勘探系统常常包括几个爆破勘探船及几个监听勘探船。此外,一些勘探船可以既是爆破勘探船又是监听勘探船。
本发明可以特别应用于石油勘探行业,但是也可以应用于任何在海洋环境中使用地球物理学数据采集网络的领域。
2.背景技术
2.1海洋地震勘探采集的原理
为了在测量区域中执行海洋地震勘探采集,一般使用一些震源(射枪、振动源等)及地震传感器。传感器容纳在电缆中,称为拖缆(或线性天线)。几个拖缆一起使用以形成数千个传感器的阵列。震源被一个或几个爆破勘探船拖拽,并且拖缆被一个或几个监听勘探船拖拽。同一个勘探船可以既是爆破勘探船又是监听勘探船(即可以拖拽一个或几个拖缆以及一个或几个震源)。
为了在海洋环境中收集地球物理学数据,震源(被至少一个爆破勘探船拖拽)被激活以产生单个脉冲或连续扫描的能量。每个震源产生的信号行进穿过地壳层,并且拖缆(被至少一个监听勘探船拖拽)中的传感器(水听器)捕获所反射的信号。
在每次爆破之后,会创建两个文件,第一文件含有地震传感器(包括在地震勘探拖缆中)提供的地震勘探数据;第二文件称为RH文件(即“记录标头”),含有关于爆破勘探船上的爆破的信息(射枪标头(GH)、真实爆破时间(FTB,即,“现场起爆信号”)及爆破时间的震源位置(SP))。这些第一及第二文件接着组合形成完整的第三文件(也称为SEG-D文件)。使用对SEG-D文件中包含的地震勘探数据的解释来计算地壳的3D图像。
每个震源必须爆破的理论位置是一个爆破点位置(也称为“爆破点”),通过其地理坐标(纬度/经度及/或东坐标北坐标)加以定义。当震源到达这个爆破点时,射枪被激活并且产生爆炸。所有震源的这组爆破点被称为“预定表”。
采集过程通过导航系统(也称为INS,即“集成式导航系统”)得到控制及监测,这个导航系统在主勘探船(也称为“主速度勘探船”)上,并且其作用是计算传感器及震源的位置,根据预定表的几何形状沿勘探船的采集路径驱动勘探船,并且在期望位置(爆破点)激活震源以执行地震勘探采集。
导航系统根据各种系统组件的位置确定每个爆破点的起动时刻。这个时刻也称为“爆破时间”,通常标注为T0。
通过众所周知的测量装置(GPS、RGPS、音响装置、罗盘、深度传感器……),可以知道所有设备(水听器及射枪)的实际位置。
2.2多勘探船操作
为了进一步提高地震勘探成像的质量,现在在多勘探船操作中执行地震勘探测量,以便得到地壳的宽方位角照明(这解释了为什么在这种情况下将预定表称为“宽方位角预定表”或“WAZ预定表”)。多勘探船地震勘探系统常常包括几个爆破勘探船及几个监听勘探船。同一个勘探船可以既是爆破勘探船又是监听勘探船。
宽方位角预定表定义一连串爆破点,其中各种勘探船的爆破交错发生。
在宽方位角预定表中还定义震源的爆破次序并且因此定义勘探船的爆破次序。必须尽可能靠近宽方位角预定表中指定的爆破点的地理坐标来遵照及执行勘探船的爆破次序。为了编译爆破次序,必须使各个勘探船同步。
在图4的宽方位角预定表的简化实例中,存在三个爆破勘探船V1、V2及V3,每个勘探船分别拖拽着一个震源S1、S2及S3。我们假设爆破勘探船V1是主速度勘探船。我们还假设爆破顺序与对应爆破点的顺序相同(举例来说,第七爆破称为“爆破7”并且必须是在“爆破点7”进行的)。宽方位角预定表的这个实例可以如下继续:
●爆破勘探船V1负责爆破1、爆破4、爆破7、爆破10、爆破13等,这些爆破必须分别在爆破点1、爆破点4、爆破点7、爆破点10、爆破点13等执行;
●爆破勘探船V2负责爆破2、爆破5、爆破8、爆破11等,这些爆破必须分别在爆破点2、爆破点5、爆破点8、爆破点11等执行;以及
●爆破勘探船V3负责爆破3、爆破6、爆破9、爆破12等,这些爆破必须分别在爆破点3、爆破点6、爆破点9、爆破点12等执行。
在这个实例中,三个爆破勘探船V1、V2及V3应当对准,但对应于三个连续爆破(每个爆破由三个爆破勘探船中的不同一个勘探船执行)的三个爆破点不对准(例如,爆破点1、2及3不对准)。但是,并且如图4中所示,我们假设爆破勘探船V2在前,爆破勘探船V3在后。
现在我们介绍与一个爆破点相关联的理论爆破时间的计算。
针对每个震源,定义一条路径(也称为“航线”或“导航线”),这条路径穿过《沿途点》,包含与这个震源相关联的爆破点。
给定点X(例如,震源的参考点)的速度在这条路径上的投影称为“沿途速度”(并且标注为SA(X))。
投影在这条路径上的两个点X及Y之间的距离称为“沿途距离”(并且标注为DA(X,Y))。
用于计算理论爆破时间(T0)的点称为“预测点”(并且标注为PP)。这个点可以是位于震源上或位于拖拽这个震源的爆破勘探船上的参考点。
对于给定震源S及给定预测点PP,根据下式计算与给定爆破点SP相关联的理论爆破时间T0:T0=(DA(SP,PP)/SA(PP))+当前时间
举例来说,在图4中,针对震源S1及对应于震源S1的参考点的给定预测点PP,根据下式(也称为“距离模式下的计算”)计算与“爆破7”相关联的理论爆破时间T0:
T0(爆破7)=(DA(“爆破点7”,PP)/SA(PP))+当前时间
必须指出,DA(“爆破点7”,PP)在图4中标注为DA1。
2.3定义
*靶心(标注为BE):主勘探船是其他勘探船(从属勘探船)的参考。主勘探船(或与主勘探船相关联的任何设备,例如震源)的点用作参考点,计算其他勘探船(从属勘探船)的理想位置,即用于从属勘探船的空间同步。
从属勘探船的理想位置用称为“靶心”(BE)的圆形靶子指定,所述靶子具有:
●一个中心,其取决于主勘探船的参考点在从属勘探船的航线上的投影。在从属爆破勘探船应当与主勘探船(参见图4)对准的特定情况下,“靶心”的中心与主勘探船的参考点在从属勘探船的航线上的投影一致。在爆破勘探船不应当对准(参见图1,2A-2D及3A-3C)的特定情况下,在“靶心”的中心与主勘探船的参考点在从属勘探船的航线上的投影之间存在沿从属勘探船航线的预定偏移。举例来说,在图1中,这个偏移对于从属爆破勘探船V2的“靶心”来说等于18,75m,对于从属爆破勘探船V3的“靶心”来说等于37,5m,并且对于从属爆破勘探船V4的“靶心”来说等于56,25m;以及
●通过合同要求可以确定的公差半径(例如10m)。
事先定义的从属勘探船的参考点必须位于“靶心”中,以确保从属勘探船的恰当同步。在图4的实例中,爆破勘探船V1是主勘探船。通过圆形靶子(“靶心”,标注为BE2)指定从属爆破勘探船V2的理想位置。通过圆形靶子(“靶心”,标注为BE3)指定从属爆破勘探船V3的理想位置。
*靶心沿途距离(标注为BE DA):针对从属爆破勘探船,这是投影在这个从属爆破勘探船的路径(航线)上的靶心的中心与这个从属爆破勘探船的参考点之间的距离。在图4的实例中,从属爆破勘探船V2的“靶心沿途距离”被标注为BE DA2。从属爆破勘探船V3的“靶心沿途距离”被标注为BE DA3。“从属勘探船的参考点必须位于‘靶心’中以确保从属勘探船的恰当同步”这个前述条件,也可表达为“靶心沿途距离(BE DA)应小于或等于‘靶心’(BE)的公差半径”。
*理论爆破间距:它是两个连续爆破点之间的理论距离。基本上,两个爆破点之间的每个真实距离应接近这个间距。换句话说,理论爆破间距是沿主勘探船航线的两个理论爆破之间的距离。例如,如果主勘探船必须爆破爆破点N及爆破点N+4,并且从属勘探船必须爆破点N+1、N+2及N+3,那么我们就得出:理论爆破间距=距离(爆破点N,爆破点N+4)/(4-0)。
*爆破时间间隔(标注为STI):这是两个爆破之间的实时时间间隔。
*最小爆破时间间隔(标注为“Min STI”,并且也称为“最小爆破循环时间”):这是必须在两个连续爆破之间维持以避免任何干扰的最小时间间隔。如果不遵照这个值,则存在爆破重叠,并且这两个爆破不被视为有效。
2.4用于避免爆破重叠的已知方法
2.4.1图1展示如下定义的特定上下文中的理想场景:存在四个爆破勘探船V1、V2、V3及V4,每个勘探船分别拖拽一个震源S1、S2、S3及S4。我们假设爆破勘探船V1是主速度勘探船。我们还假设爆破顺序与对应爆破点的顺序相同(举例来说,第四爆破称为“爆破4”并且必须是在“爆破点4”进行的)。我们还假设理论爆破间距等于18,75m。
为简单起见,仅图解说明前四个爆破点(爆破点1至爆破点4):
●爆破勘探船V1负载爆破1,爆破1要在爆破点1执行。这个第一爆破的实现在图1中通过术语“砰”表示;
●爆破勘探船V2负载爆破2,爆破2要在爆破点2执行;
●爆破勘探船V3负载爆破3,爆破3要在爆破点3执行;
●爆破勘探船V4负载爆破4,爆破4要在爆破点4执行。
在这个实例中,四个爆破点(爆破1至爆破4)应当对准,但所述四个震源S1至S4(及因而四个爆破勘探船V1至V4)不对准。
如果我们实现一种情形,其中每个从属爆破勘探船相对于主爆破勘探船保持恒定速度(使得震源之间的同线距离保持不变),并且其中在从属爆破勘探船与主爆破勘探船之间不存在通信中断,则没有问题(爆破2、3及4实际上分别在爆破点2、3及4执行)。但是,这是不现实的。
2.4.2现在参看图2A至2D,我们介绍一种在包括几个拍勘探船(一个调度爆破勘探船及至少一个从属爆破勘探船)及至少一个监听勘探船的多勘探船地震勘探系统中管理爆破的第一已知方法。
这个第一已知方法使得爆破位置上的“沿途距离”(DA)误差(即,对于从属爆破勘探船进行的给定爆破,投影在从属勘探船的航线上的在理论爆破点与实际上进行爆破的位置之间的距离)减到最小。
在这个第一已知方法中,导航系统(INS)根据从属爆破勘探船V2、V3及V4的沿线进程起动震源。
如图2A中所示,我们假设主爆破勘探船V1通常在处于爆破点1上时爆破,从属爆破勘探船V2较晚并且在其“靶心”后面30m(震源S2与爆破点2之间有48,75m,而不是图1中的18,75m),并且从属爆破勘探船V3较晚并且在其“靶心”后面10m(震源S3与爆破点3之间有47,5m,而不是图1中的37,5m)。
如图2B中所示,当主爆破勘探船V1向前移动47,5m时,从属爆破勘探船V3(更确切地说是其震源S3)先到达爆破点3,然后从属爆破勘探船V2(更确切地说是其震源S2)到达爆破点2。理论上,从属爆破勘探船V3应起动其震源(因为震源S3位于爆破点3上),但是事实上,主爆破勘探船V1的导航系统(INS)假设从属爆破勘探船V2马上就要起动(离爆破点2有1,25m),并且决定在爆破点3从属爆破勘探船V3不进行爆破。
如图2C中所示,当主爆破勘探船V1向前移动了48,75m时,从属爆破勘探船V2(更确切地说是其震源S2)到达爆破点2。从属爆破勘探船V2接着最终起动,但比计划推后了一个完整爆破循环以上。必须指出,从属爆破勘探船V4将在7,5m内(即大约4秒)到达其爆破点4。
如图2D中所示,当主爆破勘探船V1向前移动了56,25m时,从属爆破勘探船V4(更确切地说是其震源S4)到达爆破点4。主爆破勘探船V1的导航系统(INS)可以被配置成起动从属爆破勘探船V4的爆破,或者抑制从属爆破勘探船V4的爆破(考虑到它离从属爆破勘探船V2的爆破记录太近(4秒))。
在从属爆破勘探船V2未就位的时候,将重复图2A至2D的整个问题。
这些图2A至2D展现了当试图恰好在预定表靶上起动每个震源(即,恰好在与这个震源相关联的爆破点上)时的问题的复杂性。
如上文详述,这个第一已知方法的缺点是:
●由于从属爆破勘探船未就位(例如在图2B中,由于从属爆破勘探船V2未就位,所以从属爆破勘探船V3未进行爆破),导致未能爆破;以及
●由于因为来自前一次爆破的记录的重要部分当中起动了爆破(例如在图2D中,从属爆破勘探船V4的爆破接近与从属爆破勘探船V2的爆破相关的记录),所以爆破出现重叠。
在这个第一已知方法的替代实施例中,导航系统(INS)还检查是否满足以下条件:“BE DA≤BE半径”。这确保了合适的爆破时间间隔(STI),同时如果BE半径根据震源的速度得到调适(即,如果BE半径与勘探船的速度相比不太长),则可以保证起动次序。使用这个替代实施例,(从属爆破勘探船V2的)震源S2将未进行其爆破(因为BE DA2=30m>BE半径=10m)并且(从属爆破勘探船V3的)震源S3将未进行其爆破(因为:BE DA3=10m<=BE半径=10m)。
2.4.3现在参看图3A至3C,我们介绍一种在包括几个爆破勘探船(一个调度爆破勘探船及至少一个从属爆破勘探船)及至少一个监听勘探船的多勘探船地震勘探系统中管理爆破的第二已知方法。
这个第二已知方法使得可以根据主爆破勘探船的沿线进程起动震源。不论爆破器的位置在哪里,并且无需考虑DA误差,进行每个爆破点的爆破。换句话说,这个第二已知方法提议改变主爆破勘探船V1的导航系统(INS)的操作模式。更确切地说,导航系统(INS)仅根据主爆破勘探船V1的沿线进程起动震源,而并不检查是否满足“BE DA≤BE半径”这个条件。从属爆破勘探船V2、V3及V4的真实爆破时间设定成主爆破勘探船V1的导航系统(INS)预测的理论爆破时间(不管从属爆破勘探船的实际位置如何)。
我们假设与图2A中相同的场景:主爆破勘探船V1一般在其位于爆破点1上时爆破,从属爆破勘探船V2较晚(在其“靶心”后面30m),并且从属爆破勘探船V3也较晚(在其“靶心”后面10m)。
如图3A中所示,当主爆破勘探船V1向前移动了18,75m时,从属爆破勘探船V2的震源S2起动,即使从属爆破勘探船V3在从属爆破勘探船V2前面。这意味着,在这种情况下,从属爆破勘探船V2“离正确位置”30m(离爆破点2)起动,但第一已知方法的操作方面的顾虑消除了(不会错过爆破)。
如图3B中所示,当主爆破勘探船V1又向前移动18,75m时,从属爆破勘探船V3的震源S3被起动,这一点不同于图2A至2D的第一场景中,从而只“离正确位置”10m(离爆破点3)。
如图3C中所示,当主爆破勘探船V1又向前移动了18,75m时,从属爆破勘探船V4的震源S4起动。
这个第二已知方法的缺点是从属爆破勘探船可能错过其爆破点,即其同线靶子(例如,在图3A中,从属爆破勘探船V2“离开位置”30m起动,并且在图3B中,从属爆破勘探船V3“离开位置”10m起动)。
2.4.4第一及第二已知方法共有的其他缺点
在前述的第一及第二已知方法中,爆破管理集中在主爆破勘探船的导航系统(INS)中。主爆破勘探船进行所有决策,并且将激活信号(爆破命令)发送给从属爆破勘探船,以便命令起动从属爆破勘探船执行的各次爆破。实际上,这些激活信号包含在信息流中,主爆破勘探船与从属爆破勘探船之间经由无线电链路(无线信道)实时交换这个信息流。
无线电链路的缺点是它并不是100%可靠的,原因在于衰减、勘探船之间的长距离、多路径及漂浮障碍等。换句话说,在地震勘探测量过程中,勘探船之间的无线电链路可能丢失或出现故障(遭到破坏)。如果当爆破命令被发射到从属勘探船的震源时发生这种情况,则将接收不到爆破命令,将不会进行爆破,并且勘探船将错过其中必需得到数据的点(爆破点)。
无线电链路的另一个缺点是它需要定期校准。通常“离线”执行校准,结果是校准之间可能发生时序误差,这些时序误差未被检测到。
3.本发明的目的
在至少一个实施例中,本发明的目的尤其是克服现有技术的这些不同缺点。
更确切地说,本发明的至少一个实施例的目的是提供一种用于在包括几个爆破勘探船的多勘探船地震勘探系统中管理爆破的技术,这种技术可以避免错过爆破,即使勘探船未与预定表良好对准。
本发明的至少一个实施例的另一目的是提供这种类型的一种技术,使得可以将爆破位置上的“沿途距离”(DA)误差(如上文所定义)减到最小。
本发明的至少一个实施例的另一个目的是提供这种类型的一种技术,可以避免爆破重叠。
本发明的至少一个实施例的另一个目的是提供这种类型的一种技术,即使在勘探船之间的无线电链路出现临时丢失的情况下也可以继续操作。
本发明的至少一个实施例的另一个目的是提供这种类型的一种技术,实施起来简单而且成本不高。
4.发明内容
本发明的特定实施例提出一种在包括一个调度爆破勘探船及至少一个从属爆破勘探船的多勘探船地震勘探系统中管理爆破的方法,其中这种方法包括,针对每个所述从属爆破勘探船:
a)作为调度爆破勘探船及与所述调度爆破勘探船相关联的爆破点的速度值的函数,计算一系列第一理论爆破时间,每个第一理论爆破时间与调度爆破勘探船的接下来一个爆破相关联;
b)作为从属爆破勘探船及与从属爆破勘探船相关联的爆破点的速度值的函数,计算一系列第二理论爆破时间,每个第二理论爆破时间与从属爆破勘探船的接下来一个爆破相关联;
c)作为所述系列的第一理论爆破时间的函数,计算一系列内插虚拟爆破时间,至少包括与就在从属爆破勘探船的接下来的爆破之前及之后的爆破相关联的内插虚拟爆破时间;
d)作为所述系列的第一理论爆破时间的函数,计算所述系列内插虚拟爆破时间及最小爆破时间间隔,每个爆破时间窗与从属爆破勘探船的接下来一个爆破相关联;
e)针对从属爆破勘探船的每个接下来的爆破:
*如果相关联的第二理论爆破时间是在相关联的爆破时间窗中,则选择相关联的第二理论爆破时间作为预测爆破时间;
*否则的话,选择最接近相关联的第二理论爆破时间的相关联的爆破时间窗的边界作为预测爆破时间。
因而,这个特定实施例依赖于一种完全新颖并且有创造性的方法,其中,针对每个爆破,我们计算一个爆破时间窗(尤其取决于第一理论爆破时间)及第二理论爆破时间,并且我们比较这个爆破时间窗及第二理论爆破时间以决定必须使用哪个预测爆破:
●如果第二理论爆破时间是在爆破时间窗中,那么这意味着从属爆破勘探船与预定表对准良好,并且从属爆破勘探船可以在第二理论爆破时间(即在“距离模式”中)爆破。“沿途距离”误差将接近于零,并且爆破点将得到最佳匹配;
●如果第二理论爆破时间不在爆破时间窗中,则这意味着从属爆破勘探船未与预定表良好对准,并且从属爆破勘探船只能在爆破时间窗的适当边界(即,在“时间模式”中)爆破,并且不是在第二理论爆破时间(即,不在“距离模式”中)爆破。“沿途距离”误差将尽可能最佳,并且爆破点将不被匹配。
因而,不会错过爆破,即使从属爆破勘探船未与预定表良好对准,并且爆破位置上的“沿途距离”(DA)误差减到最少。
根据一个特定特征,所述方法包括:
f)当已执行前一次爆破n-1,并且从属爆破勘探船负责接下来一次爆破n时:
*如果前一次爆破n-1的最终爆破时间与接下来的爆破n的预测爆破时间之间的时间间隔大于或等于最小爆破时间间隔,则选择接下来的爆破n的预测爆破时间作为从属爆破勘探船的接下来的爆破n的最终爆破时间;
*否则的话,选择前一次爆破n-1的最终爆破时间加最小爆破时间间隔作为从属爆破勘探船的接下来的爆破n的最终爆破时间。
因而,不存在爆破重叠。
根据一个特定特征,所述多勘探船地震勘探系统包括一个主速度勘探船,并且:
-如果所述主速度勘探船是爆破勘探船,那么调度爆破勘探船是主速度勘探船;
-否则的话,调度爆破勘探船是具有最小靶心沿途距离的爆破勘探船。
因而,调度爆破勘探船通常是主速度勘探船,这简化了所提议的技术的实施方案。
根据一个特定特征,如果确定条件得到验证,那么所述方法进一步包括:
-在步骤a)中计算代替真实爆破预测点使用的虚拟爆破预测点,以便计算第一系列理论爆破时间。
原则是在所述第一系列理论爆破时间中故意引入误差(使用虚拟爆破预测点而不是真实爆破预测点),以便使爆破时间窗移位。因而,在第二理论爆破时间不在爆破时间窗中的情况下,从属爆破勘探船的爆破位置上的“沿途距离”(DA)误差进一步减少。
根据一个特定特征,所述多勘探船地震勘探系统包括一个主速度勘探船,并且所述确定的条件是:
|BE DA平均值|>“重叠公差”,其中:
●|BE DA平均值|,“BE DA平均值”的绝对值,“BE DA平均值”被定义为:
○所述多勘探船地震勘探系统中包括的从属爆破勘探船的“靶心沿途距离”的平均值,如果主速度勘探船是爆破勘探船的话;或者
○从属爆破勘探船的“靶心沿途距离”的平均值减去调度爆破勘探船的“靶心沿途距离”,如果主速度勘探船不是爆破勘探船的话;
●“重叠公差”=X-(Y*Z),其中:
○X是理论爆破间距;
○Y是最小爆破时间间隔;
○Z是沿途速度,在调度爆破勘探船的真实爆破预测点确定。
因而容易决定在预测调度爆破勘探船时引入误差是否适当。
在第一特定实施方案中,调度爆破勘探船执行步骤a),并且接着执行另一步骤a'),即将第一系列理论爆破时间发送到所述至少一个从属爆破勘探船,并且从属爆破勘探船执行另一步骤b'),即接收第一系列理论爆破时间,并且接着执行所述步骤b)至e)及f),如果执行了f)的话。
这是非集中化实施方案。每个从属爆破勘探船自主进行爆破决策。此外,因为从属爆破勘探船接收到一系列第一理论爆破时间,所以它可以继续操作,即使在与调度爆破勘探船的无线电链路出现临时丢失的情况下也是如此。
调度爆破勘探船的操作减少。调度爆破勘探船仅计算及发送与其自身的爆破相关的第一系列理论爆破时间(而不是像在现有技术解决方案中,发送与所有从属爆破勘探船的爆破相关的爆破命令)。
根据一个特定特征,调度爆破勘探船反复地执行步骤a)及a'),从而引起经由无线电链路将及时连续更新的多个第一系列理论爆破时间连续发送到从属爆破勘探船。
这允许在与调度爆破勘探船的无线电链路出现临时丢失的情况下优化所述系统的操作。
根据一个特定特征,如果调度爆破勘探船检测到无线电链路无法使用,则调度爆破勘探船停止产生新的经更新的的第一系列理论爆破时间,并且根据最后产生的第一系列理论爆破时间执行接下来的爆破,直到无线电链路再次可以使用为止。
因而,可以补偿在无线电链路切断的时间流逝过程中可能发生的时移。
在第二特定实施方案中,调度爆破勘探船执行步骤a)至e),及f),如果执行f)的话;并且接着执行另一将步骤e)(或f),如果执行f)的话)的输出发送到至少一个从属爆破勘探船的步骤。从属爆破勘探船执行接收及使用所述输出的另一步骤。
这是集中化实施方案。
在另一个实施例中,本发明涉及一种计算机程序产品,包括用于当在计算机或处理器上执行所述程序时实施上述方法(在任何其不同实施例中)的程序代码指令。
在另一个实施例中,本发明涉及一种非暂时计算机可读载体媒体,其存储一个程序,所述程序当被计算机或处理器执行时,使得所述计算机或处理器执行上述方法(在其任何不同实施例中)。
在另一个实施例中,本发明涉及一种多勘探船地震勘探系统,其包括一个调度爆破勘探船及至少一个从属爆破勘探船,其中所述系统包括以下用于管理爆破的装置,这些装置集成在所述调度爆破勘探船及/或所述至少一个从属爆破勘探船中:
-用于作为调度爆破勘探船及与所述调度爆破勘探船相关联的爆破点的速度值的函数,针对每个所述从属爆破勘探船计算一系列第一理论爆破时间的装置,每个第一理论爆破时间与调度爆破勘探船的接下来一次爆破相关联;
-用于作为从属爆破勘探船及与从属爆破勘探船相关联的爆破点的速度值的函数,计算一系列第二理论爆破时间的装置,每个第二理论爆破时间与从属爆破勘探船的接下来一个爆破相关联;
-用于作为所述第一系列理论爆破时间的函数,计算一系列内插虚拟爆破时间的装置,所述一系列内插虚拟爆破时间至少包括与就在从属爆破勘探船的接下来的爆破之前及之后的爆破相关联的内插虚拟爆破时间;
-用于作为所述第一系列理论爆破时间的函数,计算所述系列内插虚拟爆破时间及最小爆破时间间隔的装置,每个爆破时间窗与从属爆破勘探船的接下来的一个爆破相关联;
-用于针对从属爆破勘探船的每个接下来的爆破进行选择的装置,所述装置如下操作:
*如果相关联的第二理论爆破时间是在相关联的爆破时间窗中,则用于选择相关联的第二理论爆破时间作为预测爆破时间的装置;
*否则的话,用于选择最接近相关联的第二理论爆破时间的相关联的爆破时间窗的边界作为预测爆破时间的装置。
有利的是,多勘探船地震勘探系统包括用于在任何其各种实施例中实施如上文所描述的方法的步骤的装置。
5.附图说明
通过以下描述将体现本发明的实施例的其他特征及优点,所述描述是借助于指示并且非穷举的实例并且通过附图提供的,附图中:
-已经相对于现有技术论述的图1展示了用于在包括几个爆破勘探船的多勘探船地震勘探系统中管理爆破的理想场景;
-已经相对于现有技术论述的图2A至2D图解说明避免爆破重叠的第一种已知方法;
-已经相对于现有技术论述的图3A至3C图解说明避免爆破重叠的第二种已知方法;
-已经相对于现有技术论述的图4展示了宽方位角预定表的简化实例;
-图5是用于定义调度爆破勘探船的算法的特定实施例的流程图;
-图6是根据本发明的特定实施例的爆破管理方法的流程图,左侧部分是通过调度爆破勘探船实施的,而右侧部分是通过每个从属爆破勘探船实施的;
-图7展示几个时序图,其使用图4的宽方位角预定表的实例图解说明图6的爆破管理方法的操作;
-图8展示使用图4的宽方位角预定表的实例执行图6的爆破管理方法的结果;
-图9是用于调度爆破勘探船实施的部分(图6的左侧部分的变化形式)的所述爆破管理方法的替代实施例的流程图;
-图10展示使用图4的宽方位角预定表的实例执行组合图9的流程图(图6的左侧部分的变化形式)与图6的右侧部分的流程图的爆破管理方法的结果;
-图11展示根据本发明的一个特定实施例的爆破勘探船(主或从属)的简化结构。
6.具体实施方式
本发明涉及一种用于在包括几个爆破勘探船(一个调度爆破勘探船及至少一个从属爆破勘探船)的多勘探船地震勘探系统中管理爆破的方法。
现在参看图5,我们介绍用于定义调度爆破勘探船的算法的特定实施例。每个勘探船执行这种算法(在步骤50中,从未知状态),以便使得每个勘探船能够知道它是不是调度爆破勘探船。
在测试步骤51中,勘探船想知道它是不是爆破勘探船。在测试步骤51中的回答为否定的情况下,勘探船在步骤55中判断它不是调度爆破勘探船。在测试步骤51中的回答为肯定的情况下,勘探船继续进行测试步骤52,其中它想知道它是不是主速度勘探船(即,用于靶心计算的主勘探船)。
在测试步骤52中的回答为肯定的情况下,勘探船在步骤56中判断它不是调度爆破勘探船。在测试步骤52中的回答为否定的情况下,勘探船继续进行测试步骤53,其中它想知道主速度勘探船(即,用于靶心计算的主勘探船)是不是爆破勘探船。
在测试步骤53中的回答为肯定的情况下,勘探船在步骤55中判断它不是调度爆破勘探船。在测试步骤53中的回答为否定的情况下,勘探船继续进行测试步骤54,其中它想知道它是不是具有最小“靶心沿途距离”(BE DA)。
在测试步骤54中的回答为否定的情况下,勘探船在步骤55中判断它不是调度爆破勘探船。在测试步骤54中的回答为肯定的情况下,勘探船在步骤56中判断它是调度爆破勘探船。
换句话说,如果所述主速度勘探船是爆破勘探船,那么它就是调度爆破勘探船。否则的话,调度爆破勘探船是具有最小“靶心沿途距离”(BE DA)的爆破勘探船。
现在参看图6及7,我们介绍根据本发明的特定实施例的爆破管理方法。
我们首先介绍左侧部分(步骤61-62),这个部分是由调度爆破勘探船V1实施的。
在步骤61中,调度爆破勘探船计算第一系列预测,包括与其接下来的爆破点相关联的第一理论爆破时间。这个系列中的第一理论爆破时间的数目是预先确定的,或者取决于调度爆破勘探船的速度。
在图7中的标记为71的时序图中,这些第一理论爆破时间被标注为:T0(爆破4),T0(爆破7),T0(爆破10),T0(爆破13)等。
如上文详述,根据下式(距离模式下的计算)计算与给定爆破点SP相关联的第一理论爆破时间T0:
T0=(DA(SP,PP)/SA(PP))+当前时间
其中:
●SA(PP)是调度爆破勘探船V1或与其相关联的任何设备(例如,震源S1的参考点)的给定预测点PP(即,参考点)的速度在调度爆破勘探船航线上的投影;
●DA(SP,PP)是给定爆破点SP与给定预测点PP之间的投影在调度爆破勘探船航线上的距离。
在步骤62中,调度爆破勘探船将这个第一系列预测(通过标记为63的箭头表示)发送到从属爆破勘探船,并且还发送到监听勘探船。
反复地执行步骤61及62,从而引起经由无线电链路连续发送多个第一系列理论爆破时间,每个系列与调度爆破勘探船的接下来的爆破点相关联。所述连续系列被及时连续更新。
在特定实施例中,如果调度爆破勘探船检测到无线电链路不可使用,那么其停止产生新的经更新的第一系列理论爆破时间,并且根据最后产生的第一系列理论爆破时间执行其接下来的爆破,直到无线电链路再次可以使用为止。
我们现在介绍图6的右侧部分(步骤64-611),这个部分由每个从属爆破勘探船实施。在以下描述中,我们考虑从属爆破勘探船V2的情况作为实例。
在步骤64中,从属爆破勘探船V2计算一系列第二预测,包括与其接下来的爆破点相关联的第二理论爆破时间。
在图7中的标记为72的时序图中,这些第二理论爆破时间被标注为:T0(爆破5),T0(爆破8),T0(爆破11)等。
这个系列中的第二理论爆破时间的数目是预先确定的,或者取决于从属爆破勘探船的速度。
在前述公式中(距离模式中的计算),预测点PP(用于计算DA及SA)现在是从属爆破勘探船V2或与其相关联的任何设备的参考点(例如,震源S2的参考点)。
还在步骤64中,从属爆破勘探船V2计算其“靶心”(BE)及其“靶心沿途距离”(BEDA)。实际上,每个从属爆破勘探船计算其自身的BE及BE DA(并非主速度勘探船,其中BE DA定义上是无效的)。每个从属爆破勘探船将其计算结果(第二系列理论爆破时间,BE及BEDA)发送到其他勘探船。
在步骤65中,从属爆破勘探船V2作为第一系列理论爆破时间(通过调度爆破勘探船发送)的函数计算(减少)系列的内插虚拟爆破时间,其至少包括与就在从属爆破勘探船V2的接下来的爆破之前及之后的爆破相关联的内插虚拟爆破时间。
在替代实施例中,从属爆破勘探船V2计算(完整)系列的内插虚拟爆破时间,包括与所有从属爆破勘探船V2及V3的所有接下来的爆破相关联的所有内插虚拟爆破时间。
在图7中标记为73的时序图中,完整系列的内插虚拟爆破时间包括:T0int(爆破5),T0int(爆破6),T0int(爆破8),T0int(爆破9),T0int(爆破11),T0int(爆破12),等。
在这个实例中,我们得出:
T0int(爆破5)=T0(爆破4)+((1/3)*(T0(爆破7)-T0(爆破4))
T0int(爆破6)=T0(爆破4)+((2/3)*(T0(爆破7)-T0(爆破4))
相同公式适用于在调度爆破勘探船V1的两个连续爆破(例如,爆破7及10)之间,从属爆破勘探船V2及V3的每一对爆破(例如,爆破8及9)。
针对从属爆破勘探船V2,减少系列的内插虚拟爆破时间将包括:T0int(爆破6),T0int(爆破9),T0int(爆破12)等。
根据理论爆破间距与调度爆破勘探船的“沿途速度”之间的比率,计算在线路开头及末尾的内插虚拟爆破时间。
在步骤66中,从属爆破勘探船V2作为第一系列理论爆破时间(通过调度爆破勘探船V1发送)的函数计算所述系列内插虚拟爆破时间(步骤65的结果)及最小爆破时间间隔(Min STI),爆破时间窗各自与从属爆破勘探船V2的接下来一个爆破相关联。
在图7中标记为74的时序图中,从属爆破勘探船V2的爆破时间窗标注为W5、W8、W11等。在图7中,最小爆破时间间隔(Min STI)标注为“M”。
举例来说,如下获得爆破时间窗W5:
●第一(左侧)边界表示为:T0(爆破4)+Min STI
●第二(右侧)边界表示为:T0int(爆破6)-Min STI
在测试步骤67中,针对其每一个接下来的爆破,从属爆破勘探船V2检查相关联的第二理论爆破时间是否在相关联的爆破时间窗中,及:
●在检查结果为肯定的情况下,从属爆破勘探船V2转至步骤69,将这个第二理论爆破时间用作预测爆破时间;
●在检查结果为否定的情况下,从属爆破勘探船V2转至步骤68,使用最接近相关联的第二理论爆破时间的相关联爆破时间窗的边界作为预测爆破时间。
举例来说,在图7中,第二理论爆破时间T0(爆破5)(参见时序图72)是在爆破时间窗W5中,因而T0(爆破5)用作从属爆破勘探船V2的爆破5的预测爆破时间,标注为T0predicted(爆破5)。
在图7中标记为75的时序图展示了从属爆破勘探船V2的爆破5、8及9的预测爆破时间:T0predicted(爆破5)、T0predicted(爆破8)及T0predicted(爆破11)。
在步骤68或69之后(根据测试步骤67的结果),从属爆破勘探船V2在负责接下来的爆破时转至步骤610及611。在这些步骤610及611中,其在必要时检查及改变接下来的这个爆破的预测爆破时间,以便满足相对于前一爆破(由另一爆破勘探船根据预定表执行)的最小爆破时间间隔(Min STI)。
换句话说,在前述的实例中,当调度爆破勘探船V1已执行了前一爆破n-1时,从属爆破勘探船V2负责接下来的爆破n及:
-如果前一爆破n-1的最终爆破时间T0final(爆破n-1)(调度爆破勘探船V1提供的信息)与接下来的爆破n的预测爆破时间T0predicted(爆破n)之间的时间间隔d大于或等于最小爆破时间间隔(Min STI),则从属爆破勘探船V2使用接下来的爆破n的预测爆破时间T0predicted(爆破n)作为接下来的爆破n的最终爆破时间T0final(爆破n);
-否则的话,从属爆破勘探船V2使用前一爆破n-1的最终爆破时间T0final(爆破n-1)加上最小爆破时间间隔(Min STI)作为接下来的爆破n的最终爆破时间T0final(爆破n)。
在图7中标记为76的时序图中,作为一实例,我们假设针对爆破4,在预测爆破时间T0predicted(爆破4)与最终爆破T0final(爆破4)之间存在偏移Δ。我们还假设T0final(爆破4)与T0predicted(爆破5)之间的距离d小于最小爆破时间间隔(Min STI)。因此,从属爆破勘探船V2使用前一爆破4的最终爆破时间T0final(爆破4)加上最小爆破时间间隔(Min STI)作为接下来的爆破5的最终爆破时间T0final(爆破5)。
图8展示使用图4的宽方位角预定表的实例执行图6的爆破管理方法的结果。
通过标记为“爆破n”的十字展示预定表的爆破点,其中n为爆破的序号(调度爆破勘探船V1的爆破1、4、7、10、13…,从属爆破勘探船V2的爆破2、5、8、11…,及从属爆破勘探船V3的爆破3、6、9、12…)。
真实爆破的位置通过标注为“爆破n'”的圆展示,其中n是爆破的序号。
可以看出,执行了所有爆破,确保极少的覆盖。调度爆破勘探船V1的真实爆破的位置与预定表的爆破点一致。在这个实例中,出于图解说明的目的,将从属爆破勘探船V2及V3的真实爆破的位置展示为特别远离预定表的爆破点。
图9是用于调度爆破勘探船实施的部分(图6的左侧部分的变化形式)的所述爆破管理方法的替代实施例的流程图。
如图8中所图解说明,从属爆破勘探船的真实爆破的位置可能远离预定表的爆破点(理论爆破点)。为了克服这个问题,在从属勘探船不太远离其“靶心”的情况下,提议一个替代实施例,其中在从属爆破勘探船与调度爆破勘探船之间共用爆破点上的位置误差。换句话说,基本原则是在调度爆破勘探船的预测(第一理论爆破时间)中引入误差,以便使从属爆破勘探船的爆破时间窗移位。
在用前述公式(距离模式中的计算)计算第一理论爆破时间时,通过用虚拟爆破预测点替换真实爆破预测点来引入误差:
T0=(DA(SP,PP)/SA(PP))+当前时间
其中PP是调度爆破勘探船或与其相关联的任何设备的预测点(即参考点)(例如,震源的参考点)。
首先,在测试步骤91中,调度爆破勘探船必须评估引入误差是否必要及可能(即,不会起相反作用)。我们假设,针对每艘爆破勘探船,爆破预测点是与靶心的参考点相同的点。调度爆破勘探船检查以下条件是否为真:|BE DA平均值|>“重叠公差”,其中:
●|BE DA平均值|,“BE DA平均值”的绝对值,“BE DA平均值”被定义为:
○所述多勘探船地震勘探系统中包括的从属爆破勘探船的“靶心沿途距离”的平均值,如果主速度勘探船是爆破勘探船的话;或者
○从属爆破勘探船的“靶心沿途距离”的平均值减去调度爆破勘探船的“靶心沿途距离”,如果主速度勘探船不是爆破勘探船的话;
●“重叠公差”=X-(Y*Z),其中:
○X是理论爆破间距;
○Y是最小爆破时间间隔;
○Z是沿途速度,在调度爆破勘探船的真实爆破预测点确定。
在测试步骤91中的回答为否定的情况下,调度爆破勘探船不需要用虚拟爆破预测点替换真实爆破预测点。接着,其转至步骤92(等同于图6的步骤61),并且接着转至步骤93(等同于图6的步骤62)。
在测试步骤91中的回答为肯定的情况下,调度爆破勘探船必须用虚拟爆破预测点替换真实爆破预测点。步骤94-910允许决定适用于真实爆破预测点的移位,以便获得虚拟爆破预测点。
在测试步骤94中,调度爆破勘探船检查“BE DA平均值”是否大于0。
在测试步骤94中的回答为否定的情况下,调度爆破勘探船根据下式计算移位,然后转至步骤97:
移位=(“BE DA平均值”+“重叠公差”)/2
在测试步骤94中的回答为肯定的情况下,调度爆破勘探船根据下式计算移位,然后转至步骤97:
移位=(“BE DA平均值”-“重叠公差”)/2
在测试步骤97中,调度爆破勘探船检查在步骤95或96中计算的移位是否大于“理论爆破间距”/2。
在测试步骤97中的回答为肯定的情况下,调度爆破勘探船转至步骤98,其中调度爆破勘探船将在步骤95或96中计算的移位改变成“理论爆破间距”/2,然后转至步骤911。
在测试步骤97中的回答为否定的情况下,调度爆破勘探船转至测试步骤99,其中调度爆破勘探船检查在步骤95或96中计算的移位是否低于“理论爆破间距”/2。
在测试步骤99中的回答为肯定的情况下,调度爆破勘探船转至步骤910,其中调度爆破勘探船将在步骤95或96中计算的移位改变成“理论爆破间距”/2,然后转至步骤911。
在测试步骤99中的回答为否定的情况下,调度爆破勘探船保持在步骤95或96中计算的移位,并且转至步骤911。
在步骤911中,调度爆破勘探船通过沿调度爆破勘探船的航线以在步骤94-910中确定的移位值(真实爆破预测点)移位来计算虚拟爆破预测点。
在步骤912中,调度爆破勘探船使用虚拟爆破预测点(而不是步骤92中使用真实爆破预测点)计算第一系列预测(包括与其接下来的爆破点相关联的第一理论爆破时间)。
图10展示使用图4的宽方位角预定表的实例执行组合图9的流程图(图6的左侧部分的变化形式)与图6的右侧部分的流程图的爆破管理方法的结果。
与图8中一样,通过标记为“爆破n”的十字展示预定表的爆破点,其中n为爆破的序号(调度爆破勘探船V1的爆破1、4、7、10、13…,从属爆破勘探船V2的爆破2、5、8、11…,及从属爆破勘探船V3的爆破3、6、9、12…)。
真实爆破的位置通过标注为“爆破n'”的圆展示,其中n是爆破的序号。
可以看出,执行了所有爆破,确保极少的覆盖。调度爆破勘探船V1的真实爆破的位置与预定表的爆破点不一致(原因在于真实爆破预测点与虚拟爆破预测点之间的前述移位)。在这个实例中,从属爆破勘探船V2及V3的真实爆破的位置靠近预定表的爆破点。
图11展示根据本发明的一个特定实施例的爆破勘探船V(主或从属)的简化结构。
爆破勘探船V包括:
-集成导航系统INS;
-无线电通信系统R,被设计成经由无线电链路与其他勘探船交换数据;
-时钟参考CR;
-震源G,包括例如一个或几个射枪;
-射枪控制器GC,能够命令震源G及产生爆破数据。射枪控制器GC命令经由集成导航系统INS发送的二进制电信号起动爆破。当起动了爆破时,使用射枪控制器GC的二进制输出来确定真实爆破时间(也称为FTB,即“现场起爆信号”)。在每次爆破之后,射枪控制器GC产生关于射枪的信息(压力,激活的射枪,深度),这些信息聚集在“射枪标头”中(也称为GH)。射枪控制器GC向集成导航系统INS提供与每个爆破相关的爆破数据,即爆破FTB的真实时间以及射枪标头GH。
集成导航系统INS包括只读存储器(ROM)110、随机存取存储器(RAM)111及处理器112。只读存储器110是非暂时计算机可读载体媒体。其存储可执行程序代码指令,这些指令由处理器112执行以便可以实施上文所描述的方法(用于在包括几个爆破勘探船的多勘探船地震勘探系统中管理爆破的方法),及例如:
●图5的步骤50-56及图6的步骤61-62(或图9的步骤91-912),如果勘探船V是调度爆破勘探船的话;或
●图5的步骤50-56及图6的步骤64-611,如果勘探船V是从属爆破勘探船的话。
当初始化时,将前述程序代码指令从只读存储器110传送到随机存取存储器111以便由处理器112执行。随机存取存储器111同样包含用于存储这次执行所需的变量及参数的寄存器。
可以同样好地实施上述爆破管理方法的所有步骤:
●通过执行例如PC类型设备、DSP(数字信号处理器)或微控制器等可重新编程的计算机所执行的一组程序代码指令。这个程序代码指令可以存储在可拆卸的非暂时计算机可读载体媒体(例如,软盘、CD-ROM或DVD-ROM)或不可拆卸的非暂时计算机可读载体媒体中;或
●通过专用机器或组件,例如FPGA(现场可编程门阵列)、ASIC(专用集成电路)或任何专用硬件组件。
Claims (12)
1.一种用于在包括一个调度爆破勘探船及至少一个从属爆破勘探船的多勘探船地震勘探系统中管理爆破的方法,所述方法包括,针对每个所述从属爆破勘探船:
a)作为所述调度爆破勘探船及与所述调度爆破勘探船相关联的爆破点的速度值的函数,计算一系列第一理论爆破时间,每个第一理论爆破时间与所述调度爆破勘探船的接下来一个爆破相关联;
b)作为所述从属爆破勘探船及与所述从属爆破勘探船相关联的爆破点的速度值的函数,计算一系列第二理论爆破时间,每个第二理论爆破时间与所述从属爆破勘探船的接下来一个爆破相关联;
c)作为所述系列的第一理论爆破时间的函数,计算一系列内插虚拟爆破时间,至少包括与就在所述从属爆破勘探船的接下来的爆破之前及之后的爆破相关联的内插虚拟爆破时间;
d)作为所述系列的第一理论爆破时间的函数,计算所述系列内插虚拟爆破时间及最小爆破时间间隔,每个爆破时间窗与所述从属爆破勘探船的接下来一个爆破相关联;
e)针对所述从属爆破勘探船的每个接下来的爆破:
如果所述相关联的第二理论爆破时间是在所述相关联的爆破时间窗中,则选择所述相关联的第二理论爆破时间作为预测爆破时间;
否则的话,选择最接近所述相关联的第二理论爆破时间的相关联的爆破时间窗的边界作为预测爆破时间。
2.根据权利要求1所述的方法,其包括:
f)当已执行前一次爆破n-1,并且所述从属爆破勘探船负责接下来一次爆破n时:
如果前一次爆破n-1的最终爆破时间与接下来的爆破n的预测爆破时间之间的时间间隔大于或等于最小爆破时间间隔,则选择接下来的爆破n的预测爆破时间作为所述从属爆破勘探船的接下来的爆破n的最终爆破时间;
否则的话,选择前一爆破n-1的最终爆破时间加最小爆破时间间隔作为所述从属爆破勘探船的接下来的爆破n的最终爆破时间。
3.根据权利要求1所述的方法,所述多勘探船地震勘探系统包括一个主速度勘探船,其中:
如果所述主速度勘探船是爆破勘探船,那么所述调度爆破勘探船是主速度勘探船;
否则的话,所述调度爆破勘探船是具有最小靶心沿途距离的爆破勘探船。
4.根据权利要求1所述的方法,其中,则所述方法进一步包括以下计算虚拟爆破预测点的步骤:
在步骤a)中计算代替真实爆破预测点使用的虚拟爆破预测点,以便计算所述系列的第一理论爆破时间。
5.根据权利要求4所述的方法,所述多勘探船地震勘探系统包括一个主速度勘探船,其中满足以下条件时所述计算虚拟爆破预测点的步骤被执行:
|BE DA平均值|>“重叠公差”,其中:
|BE DA平均值|是“BE DA平均值”的绝对值,“BE DA平均值”被定义为:
所述多勘探船地震勘探系统中包括的所述从属爆破勘探船的“靶心沿途距离”的平均值,如果所述主速度勘探船是爆破勘探船的话;或者
所述从属爆破勘探船的“靶心沿途距离”的平均值减去所述调度爆破勘探船的“靶心沿途距离”,如果所述主速度勘探船不是爆破勘探船的话;
“重叠公差”=X-(Y*Z),其中:
X是理论爆破间距;
Y是最小爆破时间间隔;
Z是沿途速度,在所述调度爆破勘探船的真实爆破预测点确定。
6.根据权利要求1所述的方法,其中所述调度爆破勘探船执行步骤a),并且接着执行另一步骤a'),即将所述系列的第一理论爆破时间发送到所述至少一个从属爆破勘探船,并且其中所述从属爆破勘探船执行另一步骤b'),即接收所述系列的第一理论爆破时间,并且接着执行所述步骤b)至e)。
7.根据权利要求2所述的方法,其中所述调度爆破勘探船执行步骤a),并且接着执行另一步骤a'),即将所述系列的第一理论爆破时间发送到所述至少一个从属爆破勘探船,并且其中所述从属爆破勘探船执行另一步骤b'),即接收所述系列的第一理论爆破时间,并且接着执行所述步骤b)至f)。
8.根据权利要求6或7所述的方法,其中所述调度爆破勘探船反复地执行步骤a)及a'),从而引起经由无线电链路将及时连续更新的多个系列的第一理论爆破时间连续发送到所述从属爆破勘探船。
9.根据权利要求8所述的方法,其中,如果所述调度爆破勘探船检测到无线电链路无法使用,则所述调度爆破勘探船停止产生新的经更新的系列的第一理论爆破时间,并且根据最后产生的系列的第一理论爆破时间执行接下来的爆破,直到所述无线电链路再次可以使用为止。
10.根据权利要求1所述的方法,其中所述调度爆破勘探船执行步骤a)至e),并且接着执行将步骤e)的输出发送到至少一个从属爆破勘探船的另一步骤,并且其中所述从属爆破勘探船执行接收及使用所述输出的另一步骤。
11.根据权利要求2所述的方法,其中所述调度爆破勘探船执行步骤a)至f),并且接着执行将步骤f)的输出发送到至少一个从属爆破勘探船的另一步骤,并且其中所述从属爆破勘探船执行接收及使用所述输出的另一步骤。
12.一种包括一个调度爆破勘探船及至少一个从属爆破勘探船的多勘探船地震勘探系统,其中所述系统包括以下用于管理爆破的装置,所述装置集成在所述调度爆破勘探船及/或所述至少一个从属爆破勘探船中:
用于作为所述调度爆破勘探船及与所述调度爆破勘探船相关联的爆破点的速度值的函数,针对每个所述从属爆破勘探船计算一系列第一理论爆破时间的装置,每个第一理论爆破时间与所述调度爆破勘探船的接下来的一次爆破相关联;
用于作为所述从属爆破勘探船及与所述从属爆破勘探船相关联的爆破点的速度值的函数,计算一系列第二理论爆破时间的装置,每个第二理论爆破时间与所述从属爆破勘探船的接下来一个爆破相关联;
用于作为所述系列的第一理论爆破时间的函数,计算一系列内插虚拟爆破时间的装置,所述一系列内插虚拟爆破时间至少包括与就在所述从属爆破勘探船的接下来的爆破之前及之后的爆破相关联的内插虚拟爆破时间;
用于作为所述系列的第一理论爆破时间的函数,计算所述系列内插虚拟爆破时间及最小爆破时间间隔的装置,每个爆破时间窗与所述从属爆破勘探船的接下来的一个爆破相关联;
用于针对所述从属爆破勘探船的每个接下来的爆破进行选择的装置,所述装置执行如下操作:
如果所述相关联的第二理论爆破时间是在所述相关联的爆破时间窗中,则用于选择所述相关联的第二理论爆破时间作为预测爆破时间的装置;
否则的话,用于选择最接近所述相关联的第二理论爆破时间的相关联的爆破时间窗的边界作为预测爆破时间的装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP12306473.5A EP2738575B1 (en) | 2012-11-28 | 2012-11-28 | Method for managing shots in a multi-vessel seismic system |
EP12306473.5 | 2012-11-28 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103852791A CN103852791A (zh) | 2014-06-11 |
CN103852791B true CN103852791B (zh) | 2017-06-27 |
Family
ID=47435841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310627253.7A Expired - Fee Related CN103852791B (zh) | 2012-11-28 | 2013-11-28 | 用于在多勘探船地震勘探系统中管理爆破的方法 |
Country Status (8)
Country | Link |
---|---|
US (1) | US9285494B2 (zh) |
EP (1) | EP2738575B1 (zh) |
CN (1) | CN103852791B (zh) |
BR (1) | BR102013030322A2 (zh) |
CA (1) | CA2835096C (zh) |
MX (1) | MX2013013952A (zh) |
NO (1) | NO2738575T3 (zh) |
RU (1) | RU2013152725A (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2657725B1 (en) * | 2012-04-27 | 2015-01-14 | Sercel | Method and system for managing a multi-vessel seismic system |
EP2937717A1 (en) | 2014-04-24 | 2015-10-28 | Sercel | Method for managing a master vessel change in a multi-vessel seismic system |
GB2528160B (en) * | 2014-05-14 | 2021-02-10 | Sercel Rech Const Elect | Method for calculating a seismic survey |
AU2015329670A1 (en) * | 2014-10-07 | 2017-04-20 | Cgg Services Sas | Method and device for boosting low-frequencies for a marine seismic survey |
US10605941B2 (en) * | 2014-12-18 | 2020-03-31 | Conocophillips Company | Methods for simultaneous source separation |
CA2999920A1 (en) | 2015-09-28 | 2017-04-06 | Conocophillips Company | 3d seismic acquisition |
US10809402B2 (en) | 2017-05-16 | 2020-10-20 | Conocophillips Company | Non-uniform optimal survey design principles |
CN108279623A (zh) * | 2018-01-22 | 2018-07-13 | 大昌建设集团有限公司 | 一种无人机爆破安全警戒调度指挥方法 |
US11481677B2 (en) | 2018-09-30 | 2022-10-25 | Shearwater Geoservices Software Inc. | Machine learning based signal recovery |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO176860C (no) * | 1992-06-30 | 1995-06-07 | Geco As | Fremgangsmåte til synkronisering av systemer for seismiske undersökelser, samt anvendelser av fremgangsmåten |
DE69702142T2 (de) | 1996-10-17 | 2001-02-01 | Whitaker Corp | Elektrischer verbinder mit einem gehäuse und einem elektrischen kontaktelement |
US5761152A (en) * | 1996-10-29 | 1998-06-02 | Pgs Exploration (Us), Inc. | Method and system for increasing fold to streamer length ratio |
US6629037B1 (en) * | 2000-06-26 | 2003-09-30 | Westerngeco, L.L.C. | Optimal paths for marine data collection |
US7391673B2 (en) * | 2005-12-12 | 2008-06-24 | Bp Corporation North America Inc. | Method of wide azimuth seismic acquisition |
US7400552B2 (en) * | 2006-01-19 | 2008-07-15 | Westerngeco L.L.C. | Methods and systems for efficiently acquiring towed streamer seismic surveys |
US8175765B2 (en) * | 2007-12-13 | 2012-05-08 | Westerngeco L.L.C. | Controlling movement of a vessel traveling through water during a seismic survey operation |
-
2012
- 2012-11-28 NO NO12306473A patent/NO2738575T3/no unknown
- 2012-11-28 EP EP12306473.5A patent/EP2738575B1/en active Active
-
2013
- 2013-11-22 CA CA2835096A patent/CA2835096C/en active Active
- 2013-11-26 BR BR102013030322-4A patent/BR102013030322A2/pt not_active Application Discontinuation
- 2013-11-27 RU RU2013152725/28A patent/RU2013152725A/ru not_active Application Discontinuation
- 2013-11-27 US US14/092,416 patent/US9285494B2/en active Active
- 2013-11-27 MX MX2013013952A patent/MX2013013952A/es active IP Right Grant
- 2013-11-28 CN CN201310627253.7A patent/CN103852791B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CA2835096A1 (en) | 2014-05-28 |
MX2013013952A (es) | 2014-08-28 |
BR102013030322A2 (pt) | 2014-10-14 |
CA2835096C (en) | 2020-12-15 |
US20140146638A1 (en) | 2014-05-29 |
EP2738575B1 (en) | 2017-11-08 |
NO2738575T3 (zh) | 2018-04-07 |
EP2738575A1 (en) | 2014-06-04 |
US9285494B2 (en) | 2016-03-15 |
CN103852791A (zh) | 2014-06-11 |
RU2013152725A (ru) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103852791B (zh) | 用于在多勘探船地震勘探系统中管理爆破的方法 | |
EP3324215B1 (en) | Method of first arrival picking of multi-channel seismic survey data | |
US9164185B2 (en) | Near-simultaneous acquisition for borehole seismic | |
US9121956B2 (en) | Method and apparatus for determining a location to acquire geophysical data | |
NO339072B1 (no) | Hurtig 3-D-overflatemultippel-prediksjon | |
MXPA06009724A (es) | Manejo de correcciones estaticas en la prediccion de multiplos. | |
CN105723249A (zh) | 使用增强现实装置的地震勘测 | |
RU2712793C2 (ru) | Заполнение в реальном времени данных при морских сейсмических съемках с использованием независимого сейсмического источника | |
US20160259076A1 (en) | Method for determining a sail path of at least one vessel of a fleet of vessels, corresponding device, computer program product and non-transitory computer-readable carrier medium | |
CN103376467B (zh) | 用于管理多船地震系统的方法 | |
US20170219731A1 (en) | Method of visualizing and interpreting wide azimuth profile (wap) | |
EP3164737B1 (en) | Offset footprint analysis for selecting candidate lines for seismic survey | |
CN105445782A (zh) | 一种无定位拖缆多道地震勘探观测系统生成方法 | |
JP2008014830A (ja) | ハイドレートの存在領域探査方法及び探査システム | |
Smit et al. | Seismic data acquisition using ocean bottom seismic nodes at the Deimos Field, Gulf of Mexico | |
EP1040368B1 (en) | Seismic data acquisition method | |
KR102652357B1 (ko) | 머신러닝 기반 3차원 탄성파 탐사자료 내삽방법 | |
KR102630288B1 (ko) | 인공지능을 이용한 3차원 탄성파 탐사자료 내삽방법 | |
EP3338114B1 (en) | Traveling ocean bottom seismic survey | |
US20150331125A1 (en) | Method for calculating a seismic survey | |
BR112018005726B1 (pt) | Método e sistema para aquisição de um levantamento sísmico marinho, e, meio legível por computador não transitório | |
GB2528160A (en) | Method for calculating a seismic survey | |
EP2864816A1 (en) | A method for determining an equipment constrained acquisition design | |
WO2013192096A1 (en) | A method for determining an equipment constrained acquisition design |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170627 Termination date: 20201128 |
|
CF01 | Termination of patent right due to non-payment of annual fee |