CN105301650A - 一种海上拖缆采集的时移地震数据的质量监测方法和装置 - Google Patents
一种海上拖缆采集的时移地震数据的质量监测方法和装置 Download PDFInfo
- Publication number
- CN105301650A CN105301650A CN201510649411.8A CN201510649411A CN105301650A CN 105301650 A CN105301650 A CN 105301650A CN 201510649411 A CN201510649411 A CN 201510649411A CN 105301650 A CN105301650 A CN 105301650A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- lapse seismic
- data
- phase time
- shot point
- 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
Links
Abstract
本申请实施例提供一种海上拖缆采集的时移地震数据的质量监测方法和装置。该方法包括:获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,采用第一方式、第二方式、第三方式和第四方式中至少一种方式得到的两期地震数据的炮点位置的误差分布数据、羽角分布数据、羽角差分布数据、相对炮点误差分布数据、相对检波点误差分布数据以及相对炮检对误差分布数据的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。利用本申请实施例提供的技术方案可以对海上拖缆采集过程中的施工质量以及时移地震数据的采集一致性进行系统地质量监测,从而满足了对海上拖缆采集地震数据进行质量监测的需求。
Description
技术领域
本发明涉及时移地震技术领域,尤其涉及一种海上拖缆采集的时移地震数据的质量监测方法和装置。
背景技术
时移地震技术是针对同一勘探区域的多期地震响应的差异信息进行油藏动态特征研究的技术。当非储层岩石的地震响应特征不随时间变化时,通过时移地震技术可以分析出因开采因素所造成的储层内孔隙压力、孔隙度、流体成分、流体饱和度及温度等油藏特性的变化情况,从而了解储层中流体流动状况以及流体成分的变化情况,为制定或更新油气开发方案提供有效依据。
现有技术中,一般将重复性好的两期时移地震数据相减可以消除由非储层岩石的构造、岩性等形成的地震响应,得到差异地震响应。利用所述差异地震响应可以研究储层中流体流动状况以及流体成分的变化情况。但在时移地震勘探实施过程中,特别是海上拖缆采集时移地震数据时,由于受到潮汐、风浪、洋流等特定施工环境和采集条件的限制,船体在行进过程中激发的炮点位置和拖缆上各检波器的接收位置都难以精确控制。这就导致无法保证海上拖缆采集过程中的施工质量以及时移地震数据的采集一致性,从而使得海上拖缆采集的时移地震数据质量难以得到保障,无法保证两期的海上时移地震数据的重复性。
因此,现有技术中亟需一种海上拖缆采集的时移地震数据的质量监测方法,保障海上拖缆采集的时移地震数据的质量,保证采集的两期海上时移地震数据良好的重复性,能为制定或更新油气开发方案提供有效的依据。
发明内容
本申请的目的是提供一种海上拖缆采集的时移地震数据的质量监测方法和装置,以改善海上拖缆时移地震数据的质量,保证采集的两期海上时移地震数据良好的重复性,能为制定或更新油气开发方案提供有效的依据。
为了实现上述目的,本申请提供了一种海上拖缆采集的时移地震数据的质量监测方法,所述方法包括:
获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,采用下述至少一种方式得到质量监测数据:
第一方式:
分别计算所述两期时移地震数据中炮点与预设点之间的误差距离;
统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据;
第二方式:
分别计算两期时移地震数据的羽角值;
统计所述羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据;
第三方式:
获取所述第一期时移地震数据中第一炮点的羽角值;
确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值;
计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据;
统计所述羽角差数据的数值分布,得到所述两期时移地震数据的羽角差分布数据,将所述羽角差分布数据作为质量监测数据;
第四方式:
根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围;获取所述预设处理范围内的第二期时移地震数据;以及确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对;以及计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;以及统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据;
根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
在一个优选的实施例中,所述第四方式还包括:
获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量;
获取所述第一测线上炮检对的数量,将所述数量作为第二数量;
根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率;
统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
在一个优选的实施例中,所述分别计算所述两期时移地震数据中炮点与预设点之间的误差距离包括:
对两期施工重叠范围进行预设网格剖分,将所述两期施工重叠范围划分为预设数量的预设网格区间;
计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
在一个优选的实施例中,所述统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据包括:
统计所述第一期的误差距离的数值分布,得到所述第一期时移地震数据中炮点位置的误差分布数据;
统计所述第二期的误差距离的数值分布,得到所述第二期时移地震数据中炮点位置的误差分布数据。
在一个优选的实施例中,所述分别计算两期时移地震数据的羽角值包括:
获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值;
计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值;
计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
在一个优选的实施例中,所述统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据包括:
统计所述第一期测线羽角值的数值分布,得到所述第一期时移地震数据的羽角值的数值分布数据;
统计所述第二期测线羽角值的数值分布,得到所述第二期时移地震数据的羽角值的数值分布数据。
在一个优选的实施例中,所述确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点包括:
将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
在一个优选的实施例中,所述根据所述第一期时移地震数据中第一炮点的时移地震数据确定预设处理范围包括:
从所述第一期时移地震数据中获取第一炮点的位置;
将以所述第一炮点的位置为圆心,且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
在一个优选的实施例中,所述从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对包括:
判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配;
当判断结果为不匹配时,从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。
在一个优选的实施例中,所述方法还包括:
当判断结果为匹配时,将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。
在一个优选的实施例中,所述从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对包括:
计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
当所述炮点距离小于所述预先设置的炮点误差阈值,且所述检波点距离小于预先设置的检波点误差阈值,且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时,将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
一种海上拖缆采集的时移地震数据的质量监测装置,所述装置包括:预设模块、第一方式处理模块、第二方式处理模块、第三方式处理模块、第四方式处理模块以及质量监测模块;
所述预设模块用于获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,采用下述至少一种处理模块得到质量监测数据。
所述第一方式处理模块包括第一计算单元和第一统计单元;
所述第一计算单元用于分别计算所述两期时移地震数据中炮点与预设点之间的误差距离;
所述第一统计单元用于统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据;
所述第二方式处理模块包括第一获取单元和第二统计单元;
所述第一获取单元用于分别计算两期时移地震数据的羽角值;
所述第二统计单元用于统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据;
所述第三方式处理模块包括第二获取单元、第三获取单元、第二计算单元和第三统计单元;
所述第二获取单元用于获取所述第一期时移地震数据中第一炮点的羽角值;
所述第二获取单元用于确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值;
所述第二计算单元用于计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据;
所述第三统计单元用于统计所述羽角差数据的数值分布,得到所述两期时移地震数据的羽角差分布数据,将所述羽角差分布数据作为质量监测数据;
所述第四方式处理模块包括第四统计单元、第三获取单元、第四获取单元、第三计算单元和第五统计单元;
所述第四统计单元用于根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围;
所述第三获取单元用于获取所述预设处理范围内的第二期时移地震数据;
所述第四获取单元用于确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对;
第三计算单元用于计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
所述第五统计单元用于统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据;
所述质量监测模块用于根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
在一个优选的实施例中,所述第四方式处理模块还包括:
第一数量获取单元,用于获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量;
第二数量获取单元,用于获取所述第一测线上炮检对的数量,将所述数量作为第二数量;
第四计算单元,用于根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率;
第六统计单元,用于统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
在一个优选的实施例中,所述第一计算单元包括:
划分单元,用于对两期施工重叠范围进行预设网格剖分,将所述两期施工重叠范围划分为预设数量的预设网格区间;
距离计算单元,用于计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
在一个优选的实施例中,所述第一统计单元包括:
第一数据处理单元,用于统计所述第一期的误差距离的数值分布,得到所述第一期时移地震数据中炮点位置的误差分布数据;
第二数据处理单元,用于统计所述第二期的误差距离的数值分布,得到所述第二期时移地震数据中炮点位置的误差分布数据。
在一个优选的实施例中,所述第一获取单元包括:
第一羽角值获取单元,用于获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值;
第二羽角值获取单元,用于计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值;
第三羽角值获取单元,用于计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
在一个优选的实施例中,所述第二统计单元包括:
第三数据处理单元,用于统计所述第一期测线羽角值的数值分布,得到所述第一期时移地震数据的羽角值的数值分布数据;
第四数据处理单元,用于统计所述第二期测线羽角值的数值分布,得到所述第二期时移地震数据的羽角值的数值分布数据。
在一个优选的实施例中,所述从所述第二获取单元包括:
第一设置单元,用于将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
在一个优选的实施例中,所述第四统计单元包括:
位置获取单元,用于从所述第一期时移地震数据中获取第一炮点的位置;
第二设置单元,用于将以所述第一炮点的位置为圆心,且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
在一个优选的实施例中,所述第四获取单元包括:
判断单元,用于判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配;
第三设置单元,用于当所述判断单元的判断结果为不匹配时,从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。
在一个优选的实施例中,所述装置还包括:
第四设置单元,用于当所述判断单元的判断结果为匹配时,将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。
在一个优选的实施例中,所述第三设置单元包括:
第五数据处理单元,用于计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
第六数据处理单元,用于当所述炮点距离小于所述预先设置的炮点误差阈值,且所述检波点距离小于预先设置的检波点误差阈值,且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时,将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
由以上本申请实施例提供的技术方案可见,本申请实施例通过计算及统计出海上拖缆采集两期地震数据的炮点位置的误差分布数据、羽角分布数据、羽角差分布数据、相对炮点误差分布数据、相对检波点误差分布数据以及相对炮检对误差分布数据,可对海上拖缆采集过程中的施工质量以及时移地震数据的采集一致性进行系统地质量监测,从而满足了对海上拖缆采集地震数据进行质量监测的需求。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程示意图;
图2是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图;
图3是本申请实施例提供的一种所述两期时移地震数据中炮点位置的误差分布数据图;
图4是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图;
图5是本申请实施例提供的分别计算两期时移地震数据的羽角值过程的流程图;
图6是本申请实施例提供的一种两期时移地震数据的羽角分布图;
图7是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图;
图8是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图;
图9是本申请实施例提供的一种相对炮点误差分布数据、相对检波点误差分布数据以及相对炮检对误差分布数据示意图;
图10是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图;
图11是本申请一种海上拖缆采集的时移地震数据的质量监测装置的实施例的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
以下首先介绍本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例。图1是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程示意图,虽然下文描述流程包括以特定顺序出现的多个操作,但是应该清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行。结合附图1,该实施例包括:
获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,在一个具体的实施例中,所述第一期时移地震数据可以包括间距100米的4条拖缆采集的时移地震数据,且每条拖缆有间距12.5米的检波器368个。所述第二期时移地震数据可以包括间距100米的10条拖缆采集的时移地震数据,且每条拖缆有间距12.5米的检波器372个。采用下述至少一种方式得到质量监测数据:
图2是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图,结合附图2,第一方式包括:
S201:分别计算所述两期时移地震数据中炮点与预设点之间的误差距离。
在一些实施例中,可以分别计算所述两期时移地震数据中炮点与预设点之间的误差距离。具体的,可以包括:对两期施工重叠范围进行预设网格剖分,将所述两期施工重叠范围划分为预设数量的预设网格区间;计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
S202:统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据。
在一些实施例中,在步骤S201获得所述炮点与预设点之间的距离之后,可以统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据。具体的,可以统计所述第一期的误差距离的数值分布,得到所述第一期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据;可以统计所述第二期的误差距离的数值分布,得到所述第二期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据。
如图3所示的是本申请实施例提供的一种所述两期时移地震数据中炮点位置的误差分布数据图。图中包括所述第一期时移地震数据中炮点位置的误差分布数据和所述第二期时移地震数据中炮点位置的误差分布数据。
由此可见,通过本申请第一方式获得的质量监测数据可以很好的反映炮点位置偏差情况,从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
图4是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图,结合附图4,第二方式包括:
S203:分别计算两期时移地震数据的羽角值。
在一些实施例中,可以分别计算两期时移地震数据的羽角值。图5是本申请实施例提供的分别计算两期时移地震数据的羽角值过程的流程图;具体的,结合附图5,可以包括:
S510:获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值。
在一些实施例中,可以依次获取所述两期时移地震数据的每条测线上每个炮点所对应的检波点羽角值。所述检波点的羽角值可以包括航向矢量与由所述检波点坐标至其所在拖缆测线上第一个检波点坐标确定的检波矢量之间的夹角。
进一步的,所述航向矢量可以包括由采集观测系统设计时确定的当前拖缆测线上拖缆船的行进方向矢量。
进一步的,所述拖缆测线上第一个检波点可以包括所述拖缆上按距离拖船由近至远排列的第一个检波点。
进一步的,所述航向矢量与所述检波矢量间夹角范围为-180°至+180°。具体的,当由航向矢量在180°范围内按顺时针方向可旋转至所述检波矢量时,可以将所述航向矢量与所述检波矢量间夹角记作正;当由所述航向矢量在180°范围内按逆时针方向可旋转至所述检波矢量时,可以将所述航向矢量与所述检波矢量间夹角记作负。
S520:计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值。
在实际应用中,一般一个炮点对应多个检波点,可以计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值。
S530:计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
在实际应用中,一般一条测线上对应多个炮点,可以计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
S204:统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据。
在一些实施例中,在步骤S203之后,可以统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据。具体的,可以统计所述第一期测线羽角值的数值分布,得到所述第一期时移地震数据的羽角值的数值分布数据;可以统计所述第二期测线羽角值的数值分布,得到所述第二期时移地震数据的羽角值的数值分布数据。
如图6所示的是本申请实施例提供的一种两期时移地震数据的羽角分布图;图中包括第一期时移地震数据的羽角分布和第二期时移地震数据的羽角分布。由图中可见,第一期采集过程中羽角值主要分布于±10°以内,第二期期采集过程中羽角值主要分布于±5°以内。
由此可见,通过本申请第二方式获得的质量监测数据可以很好的反映检波点位置偏差情况,从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
图7是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图,结合附图7,第三方式包括:
S205:获取所述第一期时移地震数据中第一炮点的羽角值。
在一些实施例中,可以获取所述第一期时移地震数据中第一炮点的羽角值。具体的,所述第一炮点可以包括所述第一期时移地震数据中的任一炮点。
S206:确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值。
在一些实施例中,可以确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值。具体的,可以将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
S207:计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据。
在一些实施例中,在步骤S206之后,可以计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据。
S208:统计所述羽角差数据的数值分布,得到所述两期时移地震数据的羽角差分布数据,将所述羽角差分布数据作为质量监测数据。
在一些实施例中,在步骤S207之后,可以根据所述第一期时移地震数据的炮点羽角差确定所述两期时移地震数据的羽角差分布,所述羽角差分布作为质量监测数据。
由此可见,通过本申请第三方式获得的质量监测数据可以很好的反映两期检波点位置之间的偏差情况,从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
图8是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图,结合附图8,第四方式包括:
S209:根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围。
在一些实施例中,可以根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围。具体的,可以包括:从所述第一期时移地震数据中获取第一炮点的位置;将以所述第一炮点的位置为圆心,且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
进一步的,所述第一炮点可以包括所述第一期时移地震数据中的任一炮点。
S210:获取所述预设处理范围内的第二期时移地震数据。
在一些实施例中,在步骤S209之后,可以获取所述预设处理范围内的第二期时移地震数据。
S211:确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对。
在一些实施例中,在步骤S210之后,可以确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对。具体的,可以包括:判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配;当判断结果为不匹配时,从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。当判断结果为匹配时,将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。具体的,所述第一炮检对可以包括所述第一炮点中的任一炮检对。
进一步的,所述从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对可以包括:计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;当所述炮点距离小于所述预先设置的炮点误差阈值,且所述检波点距离小于预先设置的检波点误差阈值,且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时,将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
进一步的,所述预先设置的炮点误差阈值可以包括两期时移地震数据中炮点的最大位置间距误差值,且当两个炮点间的距离小于所述炮点误差阈值时,可以判断两个炮点的位置为重合。具体的,所述炮点误差阈值可以根据应用场景进行预先设定。在一个具体的实施例中,所述炮点误差阈值可以设置为70m,但本申请实施例并不以此为限。
进一步的,所述预先设置的检波点误差阈值可以包括所述两期时移地震数据中检波点的最大位置间距误差值,且当两个检波点间的距离小于所述炮点误差阈值时,可以判断两个检波点的位置为重合。具体的,所述检波点误差阈值可以根据应用场景进行预先设定。在一个具体的实施例中,所述检波点误差阈值可以设置为50m,但本申请实施例并不以此为限。
进一步的,所述预先设置的炮检对误差阈值可以包括两期时移地震数据中炮检对的最大位置间距误差阈值,所述最大位置间距误差阈值包括两期炮检对的炮点位置间距误差值与所述两期炮检对的检波点位置间距误差值之和。当所述两期炮检对的炮点位置间距误差值与所述两期炮检对的检波点位置间距误差值之和小于所述炮点误差阈值时,可以判断两个炮检对的位置为重合。具体的,所述炮检对误差阈值可以根据应用场景进行预先设定。在一个具体的实施例中,所述炮检对误差阈值可以设置为100m,但本申请实施例并不以此为限。
S212:计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和。
在一些实施例中,在步骤S210之后,可以计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和。
S213:统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据。
在一些实施例中,在步骤S212之后,可以统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据。
如图9所示的是本申请实施例提供的一种相对炮点误差分布数据、相对检波点误差分布数据以及相对炮检对误差分布数据示意图。图中包括相对炮点误差曲线ds、相对检波点误差曲线dr以及相对炮检对误差曲线ds+dr。
由此可见,通过本申请第四方式获得的质量监测数据可以很好的反映两期炮点间、检波点间以及炮检对间的位置偏差情况,从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
S214:根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
在一些实施例中,可以利用上述第一方式、第二方式、第三方式以及第四方式中的至少一种方式获得的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
由此可见,可以采用本发明实施例中的任意一种方式得到质量检测数据,然后根据质量检测数据进行质量监控,评估。也可以采用上述实施例中的一种或者多种方式的组合获取质量检测数据,这样,根据本发明提供的多种实施方式的至少两种方式获取的包括多维度信息的质量监测数据,可以更加准确的对时移地震数据进行质量分析,评估。
由此可见,本申请一种海上拖缆采集的时移地震数据的质量监测方法实施例提供的技术方案通过计算及统计出海上拖缆采集两期地震数据的炮点位置的误差分布数据、羽角分布数据、羽角差分布数据、相对炮点误差分布数据、相对检波点误差分布数据以及相对炮检对误差分布数据,可对海上拖缆采集过程中的施工质量以及时移地震数据的采集一致性进行系统地质量监测,从而满足了对海上拖缆采集地震数据进行质量监测的需求。
以下介绍本申请提供的一种海上拖缆采集的时移地震数据的质量监测方法的利用实施例。图10是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图,虽然下文描述流程包括以特定顺序出现的多个操作,但是应该清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行。结合附图10,所述方法包括:
S1010:根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围。
S1011:获取所述预设处理范围内的第二期时移地震数据。
S1012:确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对。
S1013:计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和。
S1014:统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据。
S1015:获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量。
在一些实施例中,可以获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量。
S1016:获取所述第一测线上炮检对的数量,将所述数量作为第二数量。
在一些实施例中,可以获取所述第一测线上炮检对的数量,将所述数量作为第二数量。
S1017:根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率。
在一些实施例中,可以根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率。具体的,可以将所述第一数量除以所述第二数量。
S1018:统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
在一些实施例中,可以统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
S1019:根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
在一些实施例中,可以根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。具体的,可以根据上述中至少一种质量监测数据对采集的时移地震数据的进行质量监测。
由此可见,本申请一种海上拖缆采集的时移地震数据的质量监测方法实施例提供的技术方案通过计算及统计出海上拖缆采集两期地震数据的炮点位置的误差分布数据、羽角分布数据、羽角差分布数据、相对炮点误差分布数据、相对检波点误差分布数据、相对炮检对误差分布数据以及重复率分布数据,可对海上拖缆采集过程中的施工质量以及时移地震数据的采集一致性进行系统地质量监测,从而满足了对海上拖缆采集地震数据进行质量监测的需求。
以下介绍本申请一种海上拖缆采集的时移地震数据的质量监测装置的实施例,图11是本申请一种海上拖缆采集的时移地震数据的质量监测装置的实施例的示意图,结合附图11,所述装置可以包括:预设模块1101、第一方式处理模块1102、第二方式处理模块1103、第三方式处理模块1104、第四方式处理模块1105以及质量监测模块1106;
所述预设模块1101可以用于获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,采用下述至少一种处理模块得到质量监测数据。
所述第一方式处理模块1102可以包括第一计算单元和第一统计单元;
所述第一计算单元可以用于分别计算所述两期时移地震数据中炮点与预设点之间的误差距离;
所述第一统计单元可以用于统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据;
所述第二方式处理模块1103可以包括第一获取单元和第二统计单元;
所述第一获取单元可以用于分别计算两期时移地震数据的羽角值;
所述第二统计单元可以用于统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据;
所述第三方式处理模块1104可以包括第二获取单元、第三获取单元、第二计算单元和第三统计单元;
所述第二获取单元用于获取所述第一期时移地震数据中第一炮点的羽角值;
所述第二获取单元用于可以确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值;
所述第二计算单元用于可以计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据;
所述第三统计单元用于可以统计所述羽角差数据的数值分布,得到所述两期时移地震数据的羽角差分布数据,将所述羽角差分布数据作为质量监测数据;
所述第四方式处理模块1105可以包括第四统计单元、第三获取单元、第四获取单元、第三计算单元和第五统计单元;
所述第四统计单元可以用于根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围;
所述第三获取单元可以用于获取所述预设处理范围内的第二期时移地震数据;
所述第四获取单元可以用于确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对;
第三计算单元可以用于计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
所述第五统计单元可以用于统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据;
质量监测模块可以用于根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
在一个优选的实施例中,所述第四方式处理模块1105还可以包括:
第一数量获取单元,可以用于获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量;
第二数量获取单元,可以用于获取所述第一测线上炮检对的数量,将所述数量作为第二数量;
第四计算单元,可以用于根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率;
第六统计单元,可以用于统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
在一个优选的实施例中,所述第一计算单元可以包括:
划分单元,可以用于对两期施工重叠范围进行预设网格剖分,将所述两期施工重叠范围划分为预设数量的预设网格区间;
距离计算单元,可以用于计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
在一个优选的实施例中,所述第一统计单元可以包括:
第一数据处理单元,用于统计所述第一期的误差距离的数值分布,得到所述第一期时移地震数据中炮点位置的误差分布数据;
第二数据处理单元,用于统计所述第二期的误差距离的数值分布,得到所述第二期时移地震数据中炮点位置的误差分布数据。
在一个优选的实施例中,所述第一获取单元可以包括:
第一羽角值获取单元,可以用于获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值;
第二羽角值获取单元,可以用于计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值;
第三羽角值获取单元,可以用于计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
在一个优选的实施例中,所述第二统计单元可以包括:
第三数据处理单元,用于统计所述第一期测线羽角值的数值分布,得到所述第一期时移地震数据的羽角值的数值分布数据;
第四数据处理单元,用于统计所述第二期测线羽角值的数值分布,得到所述第二期时移地震数据的羽角值的数值分布数据。
在一个优选的实施例中,所述从所述第二获取单元可以包括:
第一设置单元,可以用于将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
在一个优选的实施例中,所述第四统计单元可以包括:
位置获取单元,可以用于从所述第一期时移地震数据中获取第一炮点的位置;
第二设置单元,可以用于将以所述第一炮点的位置为圆心,且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
在一个优选的实施例中,所述第四获取单元可以包括:
判断单元,可以用于判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配;
第三设置单元,可以用于当所述判断单元的判断结果为不匹配时,从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。
在一个优选的实施例中,所述装置还可以包括:
第四设置单元,可以用于当所述判断单元的判断结果为匹配时,将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。
在一个优选的实施例中,所述第三设置单元可以包括:
第五数据处理单元,可以用于计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
第六数据处理单元,可以用于当所述炮点距离小于所述预先设置的炮点误差阈值,且所述检波点距离小于预先设置的检波点误差阈值,且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时,将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
由此可见,本申请一种海上拖缆采集的时移地震数据的质量监测方法和装置的实施例提供的技术方案通过计算及统计出海上拖缆采集两期地震数据的炮点位置的误差分布数据、羽角分布数据、羽角差分布数据、相对炮点误差分布数据、相对检波点误差分布数据、相对炮检对误差分布数据以及重复率分布数据,可对海上拖缆采集过程中的施工质量以及时移地震数据的采集一致性进行系统地质量监测,从而满足了对海上拖缆采集地震数据进行质量监测的需求。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (22)
1.一种海上拖缆采集的时移地震数据的质量监测方法,其特征在于,所述方法包括:获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,采用下述至少一种方式得到质量监测数据:
第一方式:
分别计算所述两期时移地震数据中炮点与预设点之间的误差距离;以及统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据;
第二方式:
分别计算两期时移地震数据的羽角值;以及统计所述羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据;
第三方式:
获取所述第一期时移地震数据中第一炮点的羽角值;以及确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值;以及计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据;以及统计所述羽角差数据的数值分布,得到所述两期时移地震数据的羽角差分布数据,将所述羽角差分布数据作为质量监测数据;
第四方式:
根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围;获取所述预设处理范围内的第二期时移地震数据;以及确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对;以及计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;以及统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据;
根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
2.根据权利要求1所述的方法,其特征在于,所述第四方式还包括:
获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量;
获取所述第一测线上炮检对的数量,将所述数量作为第二数量;
根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率;
统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
3.根据权利要求1或2所述的方法,其特征在于,所述分别计算所述两期时移地震数据中炮点与预设点之间的误差距离包括:
对两期施工重叠范围进行预设网格剖分,将所述两期施工重叠范围划分为预设数量的预设网格区间;
计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
4.根据权利要求1或2所述的方法,其特征在于,所述统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据包括:
统计所述第一期的误差距离的数值分布,得到所述第一期时移地震数据中炮点位置的误差分布数据;
统计所述第二期的误差距离的数值分布,得到所述第二期时移地震数据中炮点位置的误差分布数据。
5.根据权利要求1或2所述的方法,其特征在于,所述分别计算两期时移地震数据的羽角值包括:
获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值;
计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值;
计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
6.根据权利要求1或2所述的方法,其特征在于,所述统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据包括:
统计所述第一期测线羽角值的数值分布,得到所述第一期时移地震数据的羽角值的数值分布数据;
统计所述第二期测线羽角值的数值分布,得到所述第二期时移地震数据的羽角值的数值分布数据。
7.根据权利要求1或2所述的方法,其特征在于,所述确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点包括:
将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
8.根据权利要求1或2所述的方法,其特征在于,所述根据所述第一期时移地震数据中第一炮点的时移地震数据确定预设处理范围包括:
从所述第一期时移地震数据中获取第一炮点的位置;
将以所述第一炮点的位置为圆心,且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
9.根据权利要求1或2所述的方法,其特征在于,所述从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对包括:
判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配;
当判断结果为不匹配时,从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。
10.根据权利要求9所述的方法,其特征在于,所述方法还包括:
当判断结果为匹配时,将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。
11.根据权利要求9所述的方法,其特征在于,所述从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对包括:
计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
当所述炮点距离小于所述预先设置的炮点误差阈值,且所述检波点距离小于预先设置的检波点误差阈值,且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时,将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
12.一种海上拖缆采集的时移地震数据的质量监测装置,其特征在于,所述装置包括:预设模块、第一方式处理模块、第二方式处理模块、第三方式处理模块、第四方式处理模块以及质量监测模块;
所述预设模块用于获取包括第一期时移地震数据和第二期时移地震数据的两期时移地震数据,采用下述至少一种处理模块得到质量监测数据;
所述第一方式处理模块包括第一计算单元和第一统计单元;
所述第一计算单元用于分别计算所述两期时移地震数据中炮点与预设点之间的误差距离;
所述第一统计单元用于统计所述误差距离的数值分布,得到所述两期时移地震数据中炮点位置的误差分布数据,将所述炮点位置的误差分布数据作为质量监测数据;
所述第二方式处理模块包括第一获取单元和第二统计单元;
所述第一获取单元用于分别计算两期时移地震数据的羽角值;
所述第二统计单元用于统计所述测线羽角值的数值分布,得到所述两期时移地震数据的羽角值的数值分布数据,将所述羽角值的数值分布数据作为质量监测数据;
所述第三方式处理模块包括第二获取单元、第三获取单元、第二计算单元和第三统计单元;
所述第二获取单元用于获取所述第一期时移地震数据中第一炮点的羽角值;
所述第二获取单元用于确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点,获取所述第二炮点的羽角值;
所述第二计算单元用于计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值,得到所述两期时移地震数据的羽角差数据;
所述第三统计单元用于统计所述羽角差数据的数值分布,得到所述两期时移地震数据的羽角差分布数据,将所述羽角差分布数据作为质量监测数据;
所述第四方式处理模块包括第四统计单元、第三获取单元、第四获取单元、第三计算单元和第五统计单元;
所述第四统计单元用于根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围;
所述第三获取单元用于获取所述预设处理范围内的第二期时移地震数据;
所述第四获取单元用于确定所述第一炮点的第一炮检对,从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对;
第三计算单元用于计算所述第一炮检对的炮点与所述第二炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第二炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
所述第五统计单元用于统计所述炮点距离的数值分布,得到所述相对炮点误差分布数据,统计所述检波点距离的数值分布,得到所述相对检波点误差分布数据,统计所述炮点距离与所述检波点距离之和的数值分布,得到所述两期时移地震数据的相对炮检对误差分布数据,将所述相对炮点误差分布数据、所述相对检波点误差分布数据以及所述相对炮检对误差分布作为质量监测数据;
所述质量监测模块用于根据采用上述方式得到的质量监测数据对采集的时移地震数据的进行质量监测处理,获取质量监测结果。
13.根据权利要求12所述的装置,其特征在于,所述第四方式处理模块还包括:
第一数量获取单元,用于获取所述第二期时移地震数据中第一测线上与所述第一炮检对相对应的第二炮检对的数量,将所述数量作为第一数量;
第二数量获取单元,用于获取所述第一测线上炮检对的数量,将所述数量作为第二数量;
第四计算单元,用于根据所述第一数量和所述第二数量计算得到所述第一测线炮检对重复率;
第六统计单元,用于统计所述重复率的数值分布,得到所述两期时移地震数据重复率分布,将所述重复率分布作为质量监测数据。
14.根据权利要求12或13所述的装置,其特征在于,所述第一计算单元包括:
划分单元,用于对两期施工重叠范围进行预设网格剖分,将所述两期施工重叠范围划分为预设数量的预设网格区间;
距离计算单元,用于计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
15.根据权利要求12或13所述的装置,其特征在于,所述第一统计单元包括:
第一数据处理单元,用于统计所述第一期的误差距离的数值分布,得到所述第一期时移地震数据中炮点位置的误差分布数据;
第二数据处理单元,用于统计所述第二期的误差距离的数值分布,得到所述第二期时移地震数据中炮点位置的误差分布数据。
16.根据权利要求12或13所述的装置,其特征在于,所述第一获取单元包括:
第一羽角值获取单元,用于获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值;
第二羽角值获取单元,用于计算所述炮点所对应的检波点羽角值的均值,将所述检波点羽角值的均值作为所述炮点的羽角值;
第三羽角值获取单元,用于计算所述两期时移地震数据中测线上炮点的羽角值的均值,将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
17.根据权利要求12或13所述的装置,其特征在于,所述第二统计单元包括:
第三数据处理单元,用于统计所述第一期测线羽角值的数值分布,得到所述第一期时移地震数据的羽角值的数值分布数据;
第四数据处理单元,用于统计所述第二期测线羽角值的数值分布,得到所述第二期时移地震数据的羽角值的数值分布数据。
18.根据权利要求12或13所述的装置,其特征在于,所述从所述第二获取单元包括:
第一设置单元,用于将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
19.根据权利要求12或13所述的装置,其特征在于,所述第四统计单元包括:
位置获取单元,用于从所述第一期时移地震数据中获取第一炮点的位置;
第二设置单元,用于将以所述第一炮点的位置为圆心,且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
20.根据权利要求12或13所述的装置,其特征在于,所述第四获取单元包括:
判断单元,用于判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配;
第三设置单元,用于当所述判断单元的判断结果为不匹配时,从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。
21.根据权利要求20所述的装置,其特征在于,所述装置还包括:
第四设置单元,用于当所述判断单元的判断结果为匹配时,将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。
22.根据权利要求20所述的装置,其特征在于,所述第三设置单元包括:
第五数据处理单元,用于计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离,计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离,计算所述炮点距离与所述检波点距离之和;
第六数据处理单元,用于当所述炮点距离小于所述预先设置的炮点误差阈值,且所述检波点距离小于预先设置的检波点误差阈值,且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时,将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510649411.8A CN105301650A (zh) | 2015-10-09 | 2015-10-09 | 一种海上拖缆采集的时移地震数据的质量监测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510649411.8A CN105301650A (zh) | 2015-10-09 | 2015-10-09 | 一种海上拖缆采集的时移地震数据的质量监测方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105301650A true CN105301650A (zh) | 2016-02-03 |
Family
ID=55199126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510649411.8A Pending CN105301650A (zh) | 2015-10-09 | 2015-10-09 | 一种海上拖缆采集的时移地震数据的质量监测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105301650A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106527976A (zh) * | 2016-10-13 | 2017-03-22 | 中国船舶工业系统工程研究院 | 一种基于sata接口的海洋物探数据存储处理系统 |
CN108614296A (zh) * | 2018-06-06 | 2018-10-02 | 中国石油集团东方地球物理勘探有限责任公司 | 观测系统重复性确定方法及装置 |
CN113484914A (zh) * | 2021-07-13 | 2021-10-08 | 中海石油(中国)有限公司 | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
CN113568041A (zh) * | 2021-07-20 | 2021-10-29 | 中海石油(中国)有限公司 | 时移地震三维拖缆采集数据的可重复性分析方法及系统 |
CN113568042A (zh) * | 2021-07-21 | 2021-10-29 | 中海石油(中国)有限公司 | 一种提高时移地震数据可重复性的拖缆采集方法和系统 |
CN113568040A (zh) * | 2021-07-20 | 2021-10-29 | 中海石油(中国)有限公司 | 一种时移地震采集数据的可重复性分析方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030182064A1 (en) * | 2000-06-30 | 2003-09-25 | Andreas Laake | Quality control of data |
CN104487870A (zh) * | 2012-06-28 | 2015-04-01 | 雪佛龙美国公司 | 用于使用移位归一化均方根度量估计时移地震信号记录的系统和方法 |
-
2015
- 2015-10-09 CN CN201510649411.8A patent/CN105301650A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030182064A1 (en) * | 2000-06-30 | 2003-09-25 | Andreas Laake | Quality control of data |
CN104487870A (zh) * | 2012-06-28 | 2015-04-01 | 雪佛龙美国公司 | 用于使用移位归一化均方根度量估计时移地震信号记录的系统和方法 |
Non-Patent Citations (3)
Title |
---|
FRANS SMIT ET AL.: "Repeatability monitoring during marine 4D streamer acquisition", 《EAGE 67TH CONFERENCE & EXHIBITION— MADRID,SPAIN》 * |
ZHANG LIBIN ET AL.: "Repeatability analysis and reconstruction of towed streamer 4D seismic data", 《SEG DENVER 2014 ANNUAL MEETING》 * |
董凤树等: "拖缆四维地震数据采集方法实践与探讨", 《中国地球物理2013》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106527976A (zh) * | 2016-10-13 | 2017-03-22 | 中国船舶工业系统工程研究院 | 一种基于sata接口的海洋物探数据存储处理系统 |
CN106527976B (zh) * | 2016-10-13 | 2019-05-24 | 中国船舶工业系统工程研究院 | 一种基于sata接口的海洋物探数据存储处理系统 |
CN108614296A (zh) * | 2018-06-06 | 2018-10-02 | 中国石油集团东方地球物理勘探有限责任公司 | 观测系统重复性确定方法及装置 |
CN113484914A (zh) * | 2021-07-13 | 2021-10-08 | 中海石油(中国)有限公司 | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
CN113484914B (zh) * | 2021-07-13 | 2023-09-12 | 中海石油(中国)有限公司 | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
CN113568041A (zh) * | 2021-07-20 | 2021-10-29 | 中海石油(中国)有限公司 | 时移地震三维拖缆采集数据的可重复性分析方法及系统 |
CN113568040A (zh) * | 2021-07-20 | 2021-10-29 | 中海石油(中国)有限公司 | 一种时移地震采集数据的可重复性分析方法及系统 |
CN113568040B (zh) * | 2021-07-20 | 2024-01-26 | 中海石油(中国)有限公司 | 一种时移地震采集数据的可重复性分析方法及系统 |
CN113568041B (zh) * | 2021-07-20 | 2024-01-30 | 中海石油(中国)有限公司 | 时移地震三维拖缆采集数据的可重复性分析方法及系统 |
CN113568042A (zh) * | 2021-07-21 | 2021-10-29 | 中海石油(中国)有限公司 | 一种提高时移地震数据可重复性的拖缆采集方法和系统 |
CN113568042B (zh) * | 2021-07-21 | 2024-01-26 | 中海石油(中国)有限公司 | 一种提高时移地震数据可重复性的拖缆采集方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105301650A (zh) | 一种海上拖缆采集的时移地震数据的质量监测方法和装置 | |
CN107918156B (zh) | 检测海底节点采集地震数据极性的方法及装置 | |
NO342539B1 (no) | Konsistent fallvinkelestimering for seismisk avbildning | |
CN108254798B (zh) | 一种快速定位地下光缆的方法及装置 | |
CN112013286B (zh) | 管道泄漏点的定位方法及装置、存储介质、终端 | |
US11733417B2 (en) | Quality control and preconditioning of seismic data | |
CN105549083A (zh) | 一种确定地震观测系统重复性整体量度的方法和装置 | |
US9482774B2 (en) | Variable turn radius for marine vessels | |
CN105425298A (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
CN104297790A (zh) | 一种起伏海底地震观测系统中炮检距分布的定量评估方法 | |
GB2536157A (en) | Electromagnetic geophysical survey systems and methods employing electric potential mapping | |
CN109782355A (zh) | Obs检波点漂移的检测方法及装置 | |
CN107850687B (zh) | 资源推定系统以及资源推定方法 | |
CN105301635A (zh) | 一种海上拖缆时移地震数据的构建方法和装置 | |
CN109581484B (zh) | 地震纵波分量多次波总照明度指标分析方法及装置 | |
CN113568041B (zh) | 时移地震三维拖缆采集数据的可重复性分析方法及系统 | |
CN110132509A (zh) | 一种埋地管道破损点定位系统和方法 | |
CN104266646A (zh) | 一种确定终端运动方向的方法和装置 | |
CN114966855A (zh) | 探测煤岩体高应力区的方法、装置、设备和介质 | |
CN104570091B (zh) | 一种获得初至波射线的方法和装置 | |
CN112241020B (zh) | 稀疏地震数据采集中确定欠采样率的方法和装置 | |
US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
CN113568040B (zh) | 一种时移地震采集数据的可重复性分析方法及系统 | |
CN105204070B (zh) | 一种炮偏检测方法及装置 | |
CN111596357B (zh) | 一种分析海底采集节点工作状态的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160203 |
|
WD01 | Invention patent application deemed withdrawn after publication |