CN105093299B - 一种基于炮检距向量片技术优化观测系统的方法及装置 - Google Patents
一种基于炮检距向量片技术优化观测系统的方法及装置 Download PDFInfo
- Publication number
- CN105093299B CN105093299B CN201510440958.7A CN201510440958A CN105093299B CN 105093299 B CN105093299 B CN 105093299B CN 201510440958 A CN201510440958 A CN 201510440958A CN 105093299 B CN105093299 B CN 105093299B
- Authority
- CN
- China
- Prior art keywords
- observation system
- ovt
- mrow
- units
- target area
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种基于炮检距向量片技术优化观测系统的方法及装置。该方法包括获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度;将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统。通过本申请实施例所公开的技术方案,可以提供一个有利于后续OVT域数据处理和解释的观测系统。
Description
技术领域
本申请涉及地球物理勘探方法技术领域,特别涉及一种基于炮检距向量片技术优化观测系统的方法及装置。
背景技术
随着油气勘探的不断深入,油田开发对地震数据的精度要求越来越高,并逐步向宽频带、全波场方向发展。高密度宽方位三维地震勘探是现今提高地震资料品质的一种有效而现实的勘探方法。
炮检距向量片(Offset Vector Tiles,简称OVT)技术是用于对高密度宽方位三维地震勘探技术所采集的宽方位地震资料进行处理的一项新技术,该技术有利于提高地震成像精度,同时在处理过程中可以保留炮检距(即炮点和检波点之间的距离,也可以称为偏移距)和方位角的信息,在复杂构造地震成像和储层描述中可以发挥更大作用。
在进行高密度宽方位三维地震勘探时,往往需要对所设计的观测系统进行优化,以便于使所设计的观测系统更加可靠。而对观测系统进行优化主要体现在选择具有合适观测系统参数的观测系统。观测系统参数主要包括面元、炮线距、接收线距、最大偏移距及覆盖次数等。
目前,现有技术中对基于炮检距向量片技术的观测系统的优化主要是在OVT域通过分析确定反射系数密度较高、拟合曲线可能性较小的炮线距和接收线距等观测系统参数来实现对观测系统的优化。
在实现本申请过程中,发明人发现现有技术中至少存在如下问题:
现有技术中对基于炮检距向量片技术的观测系统的优化仅仅针对OVT域分析了炮线距和接收线距,无法保证利用所设计的观测系统采集到的地震数据是否可以满足后续OVT域的处理和解释。
发明内容
本申请实施例的目的是提供一种基于炮检距向量片技术优化观测系统的方法及装置,以实现提供一个有利于后续OVT域数据处理和解释的观测系统的目的。
为解决上述技术问题,本申请实施例提供一种基于炮检距向量片技术优化观测系统的方法及装置是这样实现的:
本申请实施例提供了一种基于炮检距向量片技术优化观测系统的方法,包括:
获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度;
将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统。
优选的,所述获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度包括:
依次对每个观测系统进行OVT单元划分;
依次获取每个观测系统中OVT单元所对应的OVT道集;
依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
优选的,所述对每个观测系统进行OVT单元划分包括:
对预设的观测系统进行观测系统参数检验;
根据检验结果判断是否需要对预设的观测系统进行初步优化;
在判断需要对预设的观测系统进行初步优化时,对预设的观测系统进行初步优化;
对初步优化的观测系统进行OVT单元划分。
优选的,所述获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度包括:
对所述多个观测系统中的第一观测系统进行OVT单元划分;
获取所述第一观测系统中OVT单元所对应的OVT道集;
按照获取第一观测系统所对应的OVT道集的方法,依次获取所述多个观测系统中剩余的观测系统各自所对应的OVT道集;
依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
优选的,所述依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度包括:
利用预先收集的测井和钻井资料,判断是否需要对目的层进行裂缝检测;
根据判断结果来计算所述目标区域中每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
进一步的,根据判断结果来计算所述目标区域中每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度包括:
在判断出不需要对目的层进行裂缝检测时,利用下述公式来计算所述目标区域中所述观测系统所对应的OVT单元中任意方位角所对应的炮检距均匀度:
上面两式中,N为目标区域中OVT单元的个数,P(Xi)为OVT单元中任意方位角所对应的炮检距个数,为OVT单元中任意方位角所对应的炮检距平均个数,Xi为第i个OVT单元的炮检距组,H为目的层埋深,δ为目的层埋深的变化量,为方位角,i为正整数。
优选的,根据判断结果来计算所述目标区域中每个观测系统OVT单元的炮检距与方位角的分布均匀度包括:
在判断出需要对目的层进行裂缝检测时,利用下述公式来计算所述目标区域里所述观测系统OVT单元中垂直于裂缝方向的方位角所对应的炮检距均匀度:
上面两式中,N为目标区域中OVT单元的个数,P(Xi)为OVT单元中任意方位角所对应的炮检距个数,为OVT单元中任意方位角所对应的炮检距平均个数,Xi为第i个OVT单元的炮检距组,H为目的层埋深,Δ为方位角的变化量,α为垂直于裂缝方向的方位角,i为正整数。
优选的,所述根据对比结果确定出一个适用于所述目标区域的观测系统包括:
在对比出相比剩余m-1个观测系统中OVT单元的炮检距与方位角的分布均匀度,第j个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值最小时,将第j个观测系统确定为适用于所述目标区域的观测系统,其中m为观测系统的总个数,m和j均为正整数。
优选的,所述根据对比结果确定出一个适用于所述目标区域的观测系统包括:
在对比出m个观测系统中存在有k个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值均相同时,计算所述k个观测系统各自所对应的勘探成本,其中k为正整数;
将所得到的k个勘探成本各自分别进行对比;
将勘探成本最小的观测系统确定为适用于所述目标区域的观测系统。
本申请实施例还提供了一种基于炮检距向量片技术优化观测系统的装置,包括:
获取单元,用于获取目标区域中目的层所对应的多个观测系统中炮检距向量片OVT单元的炮检距与方位角的分布均匀度;
对比单元,用于将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统。
由以上本申请实施例提供的技术方案可见,本申请实施例通过获取目标区域中目的层所对应的多个观测系统中炮检距向量片OVT单元的炮检距与方位角的分布均匀度;将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统,这实现了提供一个有利于后续OVT域数据处理和解释的观测系统的目的。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种基于炮检距向量片技术优化观测系统的方法的流程图;
图2是十字排列的炮线和检波线以及VOT单元的示意图。
图3是本申请实施例提供的一种基于炮检距向量片技术优化观测系统的装置的模块图。
具体实施方式
本申请实施例提供一种基于炮检距向量片技术优化观测系统的方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面结合附图对本申请实施例所提供的一种基于炮检距向量片技术优化观测系统的方法进行详细说明。虽然本申请提供了如下述实施例或流程图所述的方法操作步骤,但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤,操作步骤之间的执行顺序并没有限制。如图1所示,本申请实施例所提供的一种基于炮检距向量片技术优化观测系统的方法包括如下步骤:
S110:获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度。
所述目标区域可以是所述目标区域可以是整个勘探施工区域,也可以是勘探施工区域中的部分区域。
所述观测系统可以是指地震波的激发点和接收点之间的相互位置关系。所述观测系统可以是正交观测系统,优选的为具有最小满覆盖的宽方位(全方位)正交观测系统,以便于减少后续计算数据量,提高计算速度。
正交观测系统可以是指炮线和检波线垂直的观测系统。宽方位正交观测系统可以是指横纵比大于等于0.5的正交观测系统;全方位正交观测系统可以是指横纵比近似为1的正交观测系统。横纵比可以表示观测系统的横向宽度和纵向宽度之间的比值。
最小满覆盖的宽方位(全方位)正交观测系统一般可以是纵向滚动次数大于等于纵向覆盖次数、横向滚动次数大于等于横向覆盖次数的宽方位(全方位)正交观测系统。所述最小满覆盖是指使施工区域达到满覆盖的最小覆盖次数。所述满覆盖可以是指同一施工区域内的覆盖次数相同。
在项目施工前期,可以利用所收集的以往工区的地质资料、测井和钻井资料以及地震解释成果等资料,基于目标区域的勘探任务,按照常规观测系统设计方法针对目标区域预设多个观测系统,然后获取所述多个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
在一实施例中,所述多个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度可以是依次对每个观测系统进行OVT单元划分,并依次获取每个观测系统中OVT单元所对应的OVT道集;再依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元炮检距与方位角的分布均匀度。
在另一实施例中,所述获取所述多个观测系统所对应的OVT单元炮检距与方位角的分布均匀度可以是对所述多个观测系统中的第一观测系统进行OVT单元划分,并获取所述第一观测系统中OVT单元所对应的OVT道集;再按照获取第一观测系统所对应的OVT道集的方法,依次获取其他观测系统所对应的OVT道集;最后依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元炮检距与方位角的分布均匀度。
所述对每个观测系统进行OVT单元划分可以是直接对预设的观测系统进行OVT单元划分;也可以是对预设的观测系统进行观测系统参数检验,根据检验结果判断是否需要对预设的观测系统进行初步优化,在判断需要对预设的观测系统进行初步优化后,对预设的观测系统进行初步优化,然后对初步优化的观测系统进行OVT单元划分。所述对预设的观测系统进行初步优化可以是指调整所述观测系统的观测系统参数。
所述对预设的观测系统进行观测系统参数检验可以包括对面元、最大炮检距、接收线距以及覆盖次数(或覆盖密度)等参数进行检验。具体的,
(1)面元
可以检验预设的面元大小,例如检验预设的面元边长与目标地质体边长之间的大小关系,以确保预设的面元边长不得大于目标地质体边长的三分之一。
可以检验所设计的面元是否满足最高无混叠频率法则、横向分辨率的要求、考虑去噪无空间假频要求、陡倾角地层以及断层偏移成像要求等。
(2)最大炮检距
最大炮检距也可以称为最大偏移距。可以检验预设的最大炮检距是否满足目的层埋深、大于6%的速度分析精度、小于12.5%动校拉伸、反射系数稳定性对排列长度的限制以及偏移时可以归位95%的绕射波能量等需求;还可以在满足上述需求的基础上,结合波动方程正演、波场照明、实际地震资料道集、速度分析、叠加成像剖面等来检验最大炮检距。
(3)接收线距
可以检验预设的接收线距与菲涅尔带半径和最大非纵距之间的关系,以确保预设的接收线距大小小于等于一个菲涅尔带半径,同时还可能受到最大非纵距的限制。所述最大非纵距需要满足三维地震资料同一面元内不同方位角的反射同相叠加。
(4)覆盖次数
可以检验预设的覆盖次数是否是根据目标区域中构造的复杂程度、表层激发、接收条件、激发方式以及接收方式来设计的。具体的可以是检验预设的覆盖次数是否是通过分析以往不同覆盖次数的剖面信噪比或者借鉴相邻地区的覆盖次数来预设的;还可以检验预设的覆盖次数是否是根据目标地域的野外试验结果来预设的。
在检验出观测系统参数并不满足上述相应的要求后,可以判断需要对预设的观测系统进行初步优化,以确保观测系统参数中的每个参数均满足相应的要求。如有需要,还可以调整预设的观测系统参数,例如接收线距、炮线距、接收线数和/或面元等,以确保观测系统属性满足预设的对称性(例如炮点位于观测系统的观测中心)、均匀度(例如覆盖次数和/或炮检距的分布均匀度)以及采样充分性(即激发的地震波能尽可能的被接收)等要求。
所述对每个观测系统进行OVT单元划分可以是按照预设规则,将所述观测系统划分成多个小矩形,在所划分的矩形区域中标记出所述观测系统中每一炮检对所对应的共中心点,在每个矩形区域中抽取出相同或大致相同的炮检距和方位角所对应的共中心点,每一矩形区域中所抽取出的共中心点所围成的区域构成了该矩形区域中的OVT单元,所有矩形区域中的OVT单元构成了观测系统的OVT单元。
图2示出了十字排列的炮线和检波线以及对应的OVT单元。图2中,RLI表示接收线距;SLI表示炮线距。
所述对观测系统进行OVT单元划分可以包括对预设的观测系统或对优化后的观测系统进行OVT单元划分。
所述预设规则可以是按照炮线距和检波线距对观测系统进行矩形划分,例如所划分的一个矩形区域的面积为炮线距和检波线距的乘积,也可以为炮线距倍数和检波线距倍数之间的乘积。
每个炮检对均由一个炮点和一个检波点构成,每个炮检对里的炮点可以是炮线上第一预设间距内的炮点;每个炮检对里的检波点可以是检波线上第二预设间距内的检波点,所述炮线和检波线十字排列。所述第一预设间距可以是炮线距的倍数;所述第二预设间距可以是检波线距的倍数
所述OVT单元的大小可以由炮线距和检波线距来决定。优选的,所述OVT单元为炮线距和检波线距的2倍,以便于确保OVT单元属性均匀。此外,炮线距和检波线距的数值不宜过大,以免造成OVT单元内炮检距离散度和方位角离散度过大,不利于后续数据计算。
OVT单元的个数等于覆盖次数。每个OVT单元都是由炮线上的炮点和检波线上的检波点构成,由于炮线和检波线的长度均有所限制,因此OVT单元的大小取值也可能受到限制,也就是说OVT单元具有限定范围的炮检距和方位角。
在对观测系统进行OVT单元划分后,可以对所述OVT单元进行编号。例如可以按照每个OVT单元中心点的横坐标和纵坐标大小来所述OVT单元进行编号,也还可以按照每个OVT单元中心点的炮检距或方位角大小对所述OVT单元进行编号,但并不限于上述编号方式。
在对观测系统进行OVT单元划分后,也还可以根据所述观测系统来获取每个OVT单元所对应的OVT道集。所述OVT道集可以是一个5维空间数集,可以用OVT(X,Y,T,O,R)来表示。其中X可以表示每个OVT单元中心点的横向坐标,Y可以表示每个OVT单元中心点的纵向坐标,T为时间方向,O为每个OVT单元中心点的偏移距或每个OVT单元的偏移距组范围(即偏移距的数值范围),R为每个OVT单元中心点的方位角或每个OVT单元的方位角组范围(即方位角的数值范围)。
需要说明的是获取所述OVT道集与对所述OVT单元进行编号之间的执行顺序并没有限制。
在获取所有OVT单元所对应的OVT道集后,可以利用所获取的OVT道集计算所述目标区域(目标区域为满覆盖区域)中OVT单元的炮检距与方位角的分布均匀度。所述炮检距与方位角的分布均匀度可以是指在方位角分布范围内任意方位角所对应的炮检距均匀度。所述炮检距均匀度可以包括炮检对个数和/或炮检距大小均匀性,其具体的可以是指每个炮检对所对应的炮检距的间隔分布均匀性。
具体的,对于每个观测系统,在获取对应的OVT道集后,可以利用预先收集的测井和钻井资料,判断是否需要对目的层进行裂缝检测,根据判断结果来计算所述目标区域(目标区域为满覆盖区域)中OVT单元的炮检距与方位角的分布均匀度。对于每个观测系统,在判断出不需要对目的层进行裂缝检测时,可以利用下述公式来计算所述目标区域里所述观测系统的OVT单元中任意方位角所对应的炮检距均匀度:
其中
上面两式中,N为目标区域中OVT单元的个数,P(Xi)为OVT单元中任意方位角所对应的炮检距个数,为OVT单元中任意方位角所对应的炮检距平均个数,Xi为第i个OVT单元的炮检距组,H为目的层埋深,δ为目的层埋深的变化量,为方位角,i为正整数。
在判断出需要对目的层进行裂缝检测时,可以利用下述公式来计算所述目标区域里所述观测系统的OVT单元中垂直于主要裂缝方向的方位角所对应的炮检距均匀度:
其中
上面两式中,N为目标区域中OVT单元的个数,P(Xi)为OVT单元中任意方位角所对应的炮检距个数,为OVT单元中任意方位角所对应的炮检距平均个数,Xi为第i个OVT单元的炮检距组,H为目的层埋深,Δ为方位角的变化量,α为垂直于裂缝方向的方位角,i为正整数。
需要说明的是,上述计算是针对目标层为单层的情形进行的,对于目的层为多层的情形,而浅层目的层和深层目的层的计算结果差异不是很大时,可以以深层目地层的计算结果作为该观测系统中OVT单元的炮检距与方位角的分布均匀度。所述浅层和所述深层可以是根据目的层埋深大小来划分的。
S120:将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统。
在得到所有观测系统所对应的OVT单元的炮检距与方位角的分布均匀度后,可以依次将每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统。在确定出适用于所述目标区域的观测系统后,可以输出该观测系统的观测系统参数。所述预设阈值可以为1,但并不限于1。
所述依次将每个观测系统的炮检距与方位角的分布均匀度与预设阈值进行对比,根据对比结果确定出适用于所述目标区域的观测系统可以是在对比出相比其他观测系统(例如剩余的m-1个观测系统)OVT单元的炮检距与方位角的分布均匀度,某个观测系统(例如第j个观测系统,j为正整数)中OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值最小时,例如该观测系统的炮检距与方位角的分布均匀度近似为1时,可以将该观测系统确定为适用于所述目标区域的观测系统;在对比出所述多个(例如m个)观测系统中存在有若干个(例如k个,k小于等于m,m和k均为正整数)观测系统的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值的差值均相同时,计算所述k个观测系统各自所对应的勘探成本,然后将所得到的k个勘探成本进行对比,可以将勘探成本最小的观测系统确定为适用于所述目标区域的观测系统,以便于节约成本。需要说明的是,所确定的观测系统可以是预设的观测系统,也可以是对预设观测系统的观测系统参数进行调整后的观测系统。
计算每个观测系统所对应的勘探成本可以是首先计算对所述目标区域进行野外采集所需的道密度以及炮密度;然后利用所述道密度和所述炮密度,计算勘探成本。
所述道密度可以是所述目标区域的接收道数与满覆盖面积之比,即每平方公里道数,其可以表示如下:
所述炮密度可以是所述目标区域中的激发炮数与满覆盖面积之比,即每平方公里炮数,其可以表示如下:
所述勘探成本可以是满覆盖面积内炮成本与道成本总和;
M=Dr×mr+Ds×ms (7)
其中Dr为道密度,Ds为炮密度,RL为目标区域的接收线数,RLN为单个接收线的道数,SL目标区域的激发线数,SLN为单个激发线的炮数,S为满覆盖面积,mr为单个接收道数的勘探成本,ms为单个激发炮的勘探成本,M为满覆盖每平方公里的勘探成本。
通过上述步骤可以看出,本申请实施例通过对目标区域中多个观测系统进行OVT单元划分,并获取对应的OVT道集;然后利用所获取的OVT道集计算OVT单元炮检距与方位角的分布均匀度;再将每个观测系统所对应的分布均匀度与预设阈值进行对比,根据对比结果来确定出一个适用于所述目标区域的观测系统,这实现了提供一个有利于后续OVT域数据处理和解释的观测系统的目的。由于检验利用观测系统所采集到的地震数据是否满足OVT域的处理和解释,主要在于该观测系统中OVT单元的炮检距和方位角的分布情况,因此利用本申请实施例所提供的技术方案,可以实现保证利用所设计的观测系统采集到的地震数据满足后续OVT域的处理和解释的目的。此外,本申请实施例还通过计算观测系统所对应的勘探成本,并将勘探成本最小的观测系统确定为适用于所述目标区域的观测系统,因而利用该观测系统可以节约勘探成本。
本申请实施例还提供了一种基于炮检距向量片技术优化观测系统的装置,如图3所示。该装置包括获取单元310和对比单元320。其中,获取单元310用于获取目标区域中多个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度;对比单元320用于将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统。
在一实施例中,所述获取单元310可以包括(图中未示出):
划分子单元,用于依次对每个观测系统进行OVT单元划分;
获取子单元,用于依次获取每个观测系统中OVT单元所对应的OVT道集;
计算子单元,用于依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元炮检距与方位角的分布均匀度。
在一实施例中,所述划分子单元具体的可以用于对预设的观测系统进行观测系统参数检验;根据检验结果判断是否需要对预设的观测系统进行初步优化;在判断需要对预设的观测系统进行初步优化时,对预设的观测系统进行初步优化;对初步优化的观测系统进行OVT单元划分。
在另一实施例中,所述获取单元310可以包括(图中未示出):
划分子单元,用于对所述多个观测系统中的第一观测系统进行OVT单元划分;
第一获取子单元,用于获取所述第一观测系统中OVT单元所对应的OVT道集;
第二获取子单元,用于按照获取第一观测系统所对应的OVT道集的方法,依次获取所述多个观测系统中剩余的观测系统各自所对应的OVT道集;
计算子单元,用于依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元炮检距与方位角的分布均匀度。
在一实施例中,所述计算子单元具体的可以用于利用预先收集的测井和钻井资料,判断是否需要对目的层进行裂缝检测;根据判断结果来计算所述目标区域中每个观测系统OVT单元的炮检对和方位角的分布均匀度。
在一实施例中,所述对比单元320包括第一确定子单元,其用于在对比出相比剩余m-1个观测系统中OVT单元的炮检距与方位角的分布均匀度,第j个观测系统中OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值最小时,将第j个观测系统确定为适用于所述目标区域的观测系统。
在另一实施例中,所述对比单元320包括第二确定子单元,其在对比出m个观测系统中存在有k个观测系统的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值的差值均相同时,计算所述k个观测系统各自所对应的勘探成本,将所得到的k个勘探成本各自分别进行对比;将勘探成本最小的观测系统确定为适用于所述目标区域的观测系统。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (9)
1.一种基于炮检距向量片技术优化观测系统的方法,其特征在于,包括:
获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度;
将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统,
其中,所述根据对比结果确定出一个适用于所述目标区域的观测系统包括:
在对比出m个观测系统中存在有k个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值均相同时,计算所述k个观测系统各自所对应的勘探成本,其中,m为观测系统的总个数,k和m均为正整数,并且k小于或等于m;
将所得到的k个勘探成本各自分别进行对比;
将勘探成本最小的观测系统确定为适用于所述目标区域的观测系统。
2.根据权利要求1所述的方法,其特征在于,所述获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度包括:
依次对每个观测系统进行OVT单元划分;
依次获取每个观测系统中OVT单元所对应的OVT道集;
依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
3.根据权利要求2所述的方法,其特征在于,所述对每个观测系统进行OVT单元划分包括:
对预设的观测系统进行观测系统参数检验;
根据检验结果判断是否需要对预设的观测系统进行初步优化;
在判断需要对预设的观测系统进行初步优化时,对预设的观测系统进行初步优化;
对初步优化的观测系统进行OVT单元划分。
4.根据权利要求1所述的方法,其特征在于,所述获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度包括:
对所述多个观测系统中的第一观测系统进行OVT单元划分;
获取所述第一观测系统中OVT单元所对应的OVT道集;
按照获取第一观测系统所对应的OVT道集的方法,依次获取所述多个观测系统中剩余的观测系统各自所对应的OVT道集;
依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
5.根据权利要求2或4所述的方法,其特征在于,所述依次利用每个观测系统所对应的OVT道集,计算每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度包括:
利用预先收集的测井和钻井资料,判断是否需要对目的层进行裂缝检测;
根据判断结果来计算所述目标区域中每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度。
6.根据权利要求5所述的方法,其特征在于,根据判断结果来计算所述目标区域中每个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度包括:
在判断出不需要对目的层进行裂缝检测时,利用下述公式来计算所述目标区域中所述观测系统所对应的OVT单元中任意方位角所对应的炮检距均匀度:
上面两式中,N为目标区域中OVT单元的个数,P(Xi)为OVT单元中任意方位角所对应的炮检距个数,为OVT单元中任意方位角所对应的炮检距平均个数,Xi为第i个OVT单元的炮检距组,H为目的层埋深,δ为目的层埋深的变化量,为方位角,i为正整数。
7.根据权利要求5所述的方法,其特征在于,根据判断结果来计算所述目标区域中每个观测系统OVT单元的炮检距与方位角的分布均匀度包括:
在判断出需要对目的层进行裂缝检测时,利用下述公式来计算所述目标区域里所述观测系统OVT单元中垂直于裂缝方向的方位角所对应的炮检距均匀度:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>U</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>N</mi>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
上面两式中,N为目标区域中OVT单元的个数,P(Xi)为OVT单元中任意方位角所对应的炮检距个数,为OVT单元中任意方位角所对应的炮检距平均个数,Xi为第i个OVT单元的炮检距组,H为目的层埋深,Δ为方位角的变化量,α为垂直于裂缝方向的方位角,i为正整数。
8.根据权利要求1所述的方法,其特征在于,所述根据对比结果确定出一个适用于所述目标区域的观测系统包括:
在对比出相比剩余m-1个观测系统中OVT单元的炮检距与方位角的分布均匀度,第j个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值最小时,将第j个观测系统确定为适用于所述目标区域的观测系统,其中j为正整数。
9.一种基于炮检距向量片技术优化观测系统的装置,其特征在于,包括:
获取单元,用于获取目标区域中多个观测系统所对应的炮检距向量片OVT单元的炮检距与方位角的分布均匀度;
对比单元,用于将所述多个观测系统各自所对应的分布均匀度分别与预设阈值进行对比,根据对比结果确定出一个适用于所述目标区域的观测系统,
其中,所述对比单元还用于在对比出m个观测系统中存在有k个观测系统所对应的OVT单元的炮检距与方位角的分布均匀度与所述预设阈值之间的差值均相同时,计算所述k个观测系统各自所对应的勘探成本,其中,m为观测系统的总个数,k和m均为正整数,并且k小于或等于m;将所得到的k个勘探成本各自分别进行对比;将勘探成本最小的观测系统确定为适用于所述目标区域的观测系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510440958.7A CN105093299B (zh) | 2015-07-24 | 2015-07-24 | 一种基于炮检距向量片技术优化观测系统的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510440958.7A CN105093299B (zh) | 2015-07-24 | 2015-07-24 | 一种基于炮检距向量片技术优化观测系统的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093299A CN105093299A (zh) | 2015-11-25 |
CN105093299B true CN105093299B (zh) | 2017-11-10 |
Family
ID=54574153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510440958.7A Active CN105093299B (zh) | 2015-07-24 | 2015-07-24 | 一种基于炮检距向量片技术优化观测系统的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093299B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105467445B (zh) * | 2015-12-15 | 2018-02-02 | 中国石油天然气集团公司 | 一种三维地震观测系统的建立方法及装置 |
CN107797146B (zh) * | 2016-08-31 | 2019-07-23 | 中国石油化工股份有限公司 | 偏移距矢量片域振幅恢复方法及系统 |
CN106680875A (zh) * | 2017-01-05 | 2017-05-17 | 中国科学院地质与地球物理研究所 | 宽方位地震偏移速度的分析方法和装置 |
CN108181644B (zh) * | 2017-12-06 | 2020-01-07 | 中国石油天然气集团公司 | 一种采样方法及装置 |
CN110389392A (zh) * | 2018-01-31 | 2019-10-29 | 西南石油大学 | 一种基于炮检距矢量的观测系统评价方法 |
CN112415592B (zh) * | 2019-08-22 | 2022-07-08 | 中国石油化工股份有限公司 | 地震资料五维谱解析噪声压制方法、存储介质及计算设备 |
CN112444843B (zh) * | 2019-08-30 | 2024-06-18 | 中国石油化工股份有限公司 | 一种观测系统的调整方法 |
CN112578452B (zh) * | 2019-09-27 | 2024-06-18 | 中国石油化工股份有限公司 | 一种gpu加速的ovt域叠前深度偏移方法及系统 |
CN113589361B (zh) * | 2020-04-30 | 2024-03-15 | 中国石油化工股份有限公司 | 一种地震数据炮检距向量片快速划分方法 |
CN111538084B (zh) * | 2020-06-11 | 2022-02-11 | 中国石油大学(华东) | Ovt域数据转换成方位角度域成像道集的方法及系统 |
CN113009573B (zh) * | 2021-02-25 | 2023-07-18 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据采集方法、装置、计算机设备及存储介质 |
CN113960669A (zh) * | 2021-10-21 | 2022-01-21 | 中国石油化工股份有限公司 | 基于共成像点道集组合计算的反射信息增强方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6337636B1 (en) * | 1998-09-16 | 2002-01-08 | Cirrus Logic, Inc. | System and techniques for seismic data acquisition |
US6470276B1 (en) * | 1999-05-25 | 2002-10-22 | Westerngeco L.L.C. | Method of determining merit of offset and azimuth distributions in seismic data acquisition |
CN101598807A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 地震观测系统设计中确定不同面元属性差异的方法 |
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN103576213A (zh) * | 2012-07-27 | 2014-02-12 | 中国石油天然气集团公司 | 一种炮检距属性均匀性的确定方法 |
-
2015
- 2015-07-24 CN CN201510440958.7A patent/CN105093299B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6337636B1 (en) * | 1998-09-16 | 2002-01-08 | Cirrus Logic, Inc. | System and techniques for seismic data acquisition |
US6470276B1 (en) * | 1999-05-25 | 2002-10-22 | Westerngeco L.L.C. | Method of determining merit of offset and azimuth distributions in seismic data acquisition |
CN101598807A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 地震观测系统设计中确定不同面元属性差异的方法 |
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN103576213A (zh) * | 2012-07-27 | 2014-02-12 | 中国石油天然气集团公司 | 一种炮检距属性均匀性的确定方法 |
Non-Patent Citations (5)
Title |
---|
Reciprocal offset-vector tiles in various acquisition geometries;Gijs.O.Vermeer;《SEG Technicial Program Expanded Abstracts 2007》;20071231;第61-65页 * |
Wide-azimuth land processing:Fracture detection using offset vector tile technology;JAIME A.STEIN,et al.;《THE LEADING EDGE》;20101130;第1328-1337页 * |
三维观测系统属性均匀性的定量分析;吴永国等;《石油地球物理勘探》;20120630;第47卷(第3期);第361-365页 * |
基于属性评价分析的三维观测系统优化设计与应用效果;姚江;《石油物探》;20140731;第53卷(第4期);第385页第1栏第3行-第388页第1栏第15行 * |
黄土塬宽方位角地震资料OVT处理技术及效果;陈娟 等;《中国石油学会2015年物探技术研讨会论文集》;20150430;第142页第2栏第4行-第143页第1栏第27行 * |
Also Published As
Publication number | Publication date |
---|---|
CN105093299A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093299B (zh) | 一种基于炮检距向量片技术优化观测系统的方法及装置 | |
CN103076623B (zh) | 一种基于叠前相干的裂缝检测方法 | |
CN101551463B (zh) | 三维观测系统噪声压制估算方法 | |
CN104297789B (zh) | 一种三维倾角域稳相叠前时间偏移方法及系统 | |
CN103713323B (zh) | 一种全方位各向异性保幅成像与抽道集方法 | |
CN109100795B (zh) | 一种面元的炮检点布设方法、装置及系统 | |
CN107678057B (zh) | 三维地震中确定或优化检波器或设备投入量及排布的方法 | |
CN113341467A (zh) | 基于多插值方法的矿井瞬变电磁三维显示方法 | |
CN103954996B (zh) | 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置 | |
CN105467460A (zh) | 电磁勘探方法和装置 | |
CN102998698B (zh) | 一种束状海底电缆宽方位角地震数据采集方法 | |
CN105137479B (zh) | 一种面元覆盖次数的计算方法及装置 | |
CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
CN109163674B (zh) | 一种面结构光自动化三维测量中传感器测量视点规划方法 | |
CN101551464B (zh) | 一种复杂地表三维勘探初至波时间确定方法 | |
CN105467445B (zh) | 一种三维地震观测系统的建立方法及装置 | |
CN106842302A (zh) | 一种批量编辑初至的方法及装置 | |
CN102236101A (zh) | 多层裂缝预测方法和装置 | |
CN102914790B (zh) | 二维观测系统和三维观测系统一次采集的观测系统方法 | |
CN108375794B (zh) | 基于对称观测的vsp缝洞绕射成像技术方法 | |
CN105676278B (zh) | 一种处理地震波数据的方法和装置 | |
CN103217708B (zh) | 获得海上多缆拖缆地震采集全单缆单源地震剖面的方法 | |
CN109991663A (zh) | 工区地震速度体校正方法和装置 | |
CN107976714A (zh) | 一种复杂地表分级空间距离加权的超道集计算方法 | |
CN107942373A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |