CN106353818A - 一种地震数据时差校正的方法及系统 - Google Patents
一种地震数据时差校正的方法及系统 Download PDFInfo
- Publication number
- CN106353818A CN106353818A CN201610663124.7A CN201610663124A CN106353818A CN 106353818 A CN106353818 A CN 106353818A CN 201610663124 A CN201610663124 A CN 201610663124A CN 106353818 A CN106353818 A CN 106353818A
- Authority
- CN
- China
- Prior art keywords
- time
- data
- residual quantity
- seismic
- geological data
- 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.)
- Granted
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/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
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
技术领域
本申请涉及石油地震勘探中时移地震数据处理技术领域,特别涉及一种地震数据时差校正的方法及系统。
背景技术
四维地震技勘探术,也称时移地震勘探技术,是利用针对同一勘探区域的两期或多期地震响应的差异信息进行油藏动态特征研究的技术。该技术通过分析由于开采因素所造成的储层内孔隙压力、孔隙度、流体成分、流体饱和度及温度等油藏特性的变化,从而了解储层中流体流动状况以及流体成分的变化情况,为制定或更新油气开发方案提供有效依据。如果两期时移地震数据之间的地震响应差异仅包括油藏变化产生的期望振幅差异,时移地震勘探技术即可通过两期时移地震数据相减以消除由非储层岩石的构造、岩性等形成的共同地震响应,进而从相减得到的期望振幅差异来研究储层内部流体性质的变化。在时移地震勘探技术实施过程中,由于受到施工环境变化、采集仪器和采集参数差异等因素的影响,所述地震响应差异往往既包括油藏变化产生的期望振幅差异,又包括由两期地震数据时差等非油藏变化产生的非期望振幅差异。因此,需要消除两期时移地震数据时差,以获得无时差影响的期望振幅差异。
针对上述两期时移地震数据时差,目前常用的一种时差校正方法是小时窗互相关方法,其方法是利用互相关方法求取基础数据和监测数据在相同小时窗(比如时窗的窗长为8毫秒至64毫秒的时窗)位置处的时间样点的互相关极值,将求取的互相关极值作为所述小时窗中心点对应的时间样点的时差量,将小时窗移动一个时间样点,计算移动后的小时窗中心点对应的时间样点的时差量,直至完成监测地震数据中所有时间样点的时差量的计算,利用计算的时差量对监测地震数据进行时差校正。
时移地震数据时差产生的原因包括采集环境和采集参数变化引起的采集相关时差,也包括上覆储层由于流体变化引入的局部时差。
发明人发现现有技术中至少存在如下问题:我国油田往往含多个薄油层,开发方式多采用多层联采,上覆储层流体变化将引起下层地震反射信号的到达时间发生变化,形成与储层相关的两期时移地震数据局部时差。目前这种上覆储层流体变化引起的局部时差问题并没有引起工业界重视,但这部分局部时差对薄油层的两期时移地震数据的振幅差异影响较大。针对上述上覆储层流体变化引起的局部时差问题,需要选取时窗的窗长较小的时窗计算时差量,但小时窗互相关方法中所采用的互相关方法随时窗的窗长减小而计算精度降低。因此,小时窗互相关方法对上覆储层流体变化引起的两期时移地震数据局部时差的计算精度较低。
发明内容
本申请实施例的目的是提供一种地震时差校正的方法及系统,以提高地震数据时差校正的精度。
为解决上述技术问题,本申请实施例提供一种地震数据时差校正的方法及系统是这样实现的:
一种地震数据时差校正的方法,包括:
获取目标区域第一时间对应的基础地震数据,和第二时间对应的监测地震数据;
分别获取基础地震数据和监测地震数据中相同位置对应的第一地震道数据和第二地震道数据,分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点;
预设多个第一时差量,计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值;
根据所述第一方差值确定最小方差值,将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗位置处的时间样点对应的目标时差量;
根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据。
优选方案中,所述计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值,包括:
获取第一地震道数据中第一时窗位置处的时间样点的第一采样时间对应的第一地震道数据的第一振幅值;
根据第二地震道数据中第一时窗位置处的时间样点的第一采样时间减去第一时差量,得到第二地震道数据中第一时窗位置处的第二采样时间,获取第二采样时间对应的第二地震道数据的第二振幅值;
根据第一振幅值和第二振幅值,计算第一时差量对应的第一时窗位置处的时间样点的第一方差值。
优选方案中,所述根据第一振幅值和第二振幅值,计算第一时差量对应的第一时窗位置处的时间样点的第一方差值,采用下述公式实现:
公式中,第一振幅值i表示第一地震道数据中第一时窗位置处的第i个时间样点的第一振幅值,第二振幅值i为第二地震道数据中第一时窗位置处的第i个时间样点的第二振幅值。
优选方案中,所述预设多个第一时差量,包括:
选取第一时差量的第一时差量范围和第一时差增量;
根据第一时差量范围和第一时差增量,得到多个第一时差量。
优选方案中,第一时差增量为所述监测地震数据的采样间隔的四分之一。
优选方案中,所述方法还包括:检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适。
优选方案中,所述检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适,包括:
获取所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量,设定所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量构成的曲线的平滑阈值;
当所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时量差构成的曲线的平滑度大于平滑阈值,第一时窗的窗长的选取不合适;
根据所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量对所述监测地震数据进行时差校正;
将所述基础地震数据与所述监测地震数据相减的地震数据在第一主测线或第一联络测线上的非储层变化引起的校正前地震响应,同所述基础地震数据与所述经过时差校正后的所述监测地震数据相减的地震数据在第一主测线或第一联络测线上的非储层变化引起的校正后地震响应对比,设定所述校正后地震响应对比所述校正前地震响应的地震响应减弱阈值;
当所述校正后地震响应对比所述校正前地震响应的地震响应减弱度小于所述地震响应减弱阈值,第一时差量范围中最大第一时差量、第一时差量范围中最小第一时差量或第一时差增量的选取不合适。
优选方案中,所述根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据,包括:
根据所述目标时差量对所述监测地震数据进行时差校正,得到校正地震数据;
对所述校正地震数据进行辛格函数插值,得到目标地震数据。
一种地震数据时差校正的系统,包括:地震数据获取单元、时间样点获取单元、方差值计算单元、目标时差量获取单元和目标地震数据获取单元;其中,
所述地震数据获取单元,用于获取目标区域第一时间对应的基础地震数据,和第二时间对应的监测地震数据;
所述时间样点获取单元,用于分别获取基础地震数据和监测地震数据中相同位置对应的第一地震道数据和第二地震道数据,分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点;
所述方差值计算单元,用于预设多个第一时差量,计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值;
所述目标时差量获取单元,用于根据所述第一方差值确定最小方差值,将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗位置处的时间样点对应的目标时差量;
所述目标地震数据获取单元,用于根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据。
优选方案中,所述方差值计算单元包括:第一振幅值获取单元、第二振幅值获取单元和第一方差值计算单元;其中,
所述第一振幅值获取单元,用于获取第一地震道数据中第一时窗位置处的时间样点的第一采样时间对应的第一地震道数据的第一振幅值;
所述第二振幅值获取单元,用于根据第二地震道数据中第一时窗位置处的时间样点的第一采样时间减去第一时差量,得到第二地震道数据中第一时窗位置处的第二采样时间,获取第二采样时间对应的第二地震道数据的第二振幅值;
所述第一方差值计算单元,用于根据第一振幅值和第二振幅值,计算第一时差量对应的第一时窗位置处的时间样点的第一方差值。
优选方案中,所述系统还包括:检验单元;其中,
所述检验单元,用于检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适。
优选方案中,所述目标地震数据获取单元包括:校正地震数据获取单元和函数插值单元;其中,
所述校正地震数据获取单元,用于根据所述目标时差量对所述监测地震数据进行时差校正,得到校正地震数据;
所述函数插值单元,用于对所述校正地震数据进行辛格函数插值,得到目标地震数据。
本申请提供了一种地震数据时差校正的方法及系统,所述方法基于小时窗方差的方式获取并校正两期地震时差,不仅能够校正各类与采集相关的较大的时差量,也可以校正上覆储层流体引起的局部较小的时差量,消除了非储层变化引起的非期望振幅差异,显著提高了两期地震期望振幅差异精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种地震数据时差校正的方法实施例的流程图;
图2是本申请实施例中Marmousi2模型基础地震数据的振幅剖面、监测地震数据的振幅剖面以及基础地震数据和监测地震数据的振幅差剖面的示意图;
图3是本申请实施例中Marmousi2模型基础地震数据与监测地震数据之间的目标时差量剖面以及基础地震数据与目标地震数据的振幅差剖面的示意图;
图4是本申请实施例中基础地震数据和监测地震数据的振幅差在目的储层顶的沿层切片显示的示意图;
图5是本申请实施例中基础地震数据和目标地震数据的振幅差在目的储层顶的沿层切片显示的示意图;
图6是本申请一种地震数据时差校正的系统实施例的组成结构图;
图7是本申请地震数据时差校正的系统实施例中方差值计算单元的组成结构图;
图8是本申请地震数据时差校正的系统实施例中目标地震数据获取单元的组成结构图。
具体实施方式
本申请实施例提供一种地震数据时差校正的方法及系统。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请一种地震数据时差校正的方法实施例的流程图。如图1所示,所述地震数据时差校正的方法,包括以下步骤。
步骤S101:获取目标区域第一时间对应的基础地震数据,和第二时间对应的监测地震数据。
具体地,可以采用Marmousi2模型模拟目标区域第一时间对应的基础地震数据的振幅剖面,和第二时间对应的监测地震数据的振幅剖面。其中Marmousi2模型是基于Marmousi模型增加了一些储层并改进后的模型;Marmousi模型最早是由法国石油研究院(IFP)所属协会根据西非海上油田制作的,这种模型及其声波有限差分合成数据被应用于许多地球物理科研项目,例如,石油地震勘探数据处理等科研项目;Marmousi模型处理地震数据的过程是地质人员根据实际地震数据画出地质模型,如泥岩地层,砂岩地层,断层等,再根据测井速度对地质模型填上速度,得到速度模型,对速度模型做正演模拟,获得了合成的地震数据。第一时间对应的基础地震数据和第二时间对应的监测地震数据分别为同一目标区域在不同时间的地震数据。这样可以为后续时差校正的步骤做好地震数据的准备。
例如,图2是本申请实施例中Marmousi2模型基础地震数据的振幅剖面、监测地震数据的振幅剖面以及基础地震数据和监测地震数据的振幅差剖面的示意图。图2中(a)图为基础地震数据的振幅剖面图,图2中(b)图为监测地震数据的振幅剖面图,图2中(c)图为基础地震数据和监测地震数据的振幅差剖面图。图2中(a)图和(b)图的箭头所示为实际流体发生变化的储层地震数据的振幅剖面。图2中(c)图中基础地震数据和监测地震数据的振幅差可以包括:采集参数和采集环境引起的采集相关时差导致的振幅差,以及上覆储层流体变化引起的局部时差导致的振幅差。可以在所述监测地震数据中设置采集参数和采集环境引起的采集相关时差为2毫秒(ms)。
步骤S102:分别获取基础地震数据和监测地震数据中相同位置对应的第一地震道数据和第二地震道数据,分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点。
具体地,根据所述基础地震数据和所述监测地震数据,可以分别获取所述基础地震数据和所述监测地震数据中相同位置对应的第一地震道数据和第二地震道数据。可以分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点。其中,第一地震道数据和第二地震道数据可以分别为所述基础地震数据和所述监测地震数据中任意相同位置对应的地震道数据。第一时窗位置可以为第一地震道数据和第二地震道数据中任意时窗位置,其中任意时窗位置对应的时窗的窗长可以为32ms。这里选取第一时窗位置处的时间样点为后续计算第一时窗位置处的时间样点的时差做准备,同时也体现了后续计算的第一时窗位置处的时间样点的时差是随时间、空间变化的。
步骤S103:预设多个第一时差量,计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值。
具体地,可以选取第一时差量的第一时差量范围为-4ms~4ms。可以选取第一时差量的第一时差增量为所述监测地震数据的采样间隔的四分之一,即为0.25ms;其中,所述监测地震数据的采样间隔时间可以为1ms。根据第一时差量范围和第一时差增量,可以得到多个第一时差量,即为-4ms、-3.75ms、-3.5ms、…0ms、…3.5ms、3.75ms、4ms。其中省略号部分为按照所述时差增量依次增加后的第一时差量。
根据所述选取的第一地震道数据和第二地震道数据中第一时窗位置处的时间样点,可以获取第一地震道数据中第一时窗位置处的时间样点的第一采样时间对应的第一地震道数据的第一振幅值。其中第一采样时间为第一地震道数据中第一时窗位置处的任意时间样点对应的采样时间。根据第二地震道数据中第一时窗位置处的时间样点的第一采样时间减去第一时差量,可以得到第二地震道数据中的第二采样时间。可以获取第二采样时间对应的第二地震道数据中的时间样点的第二振幅值。第一振幅值和第二振幅值分别为第一地震道数据中第一采样时间和第二地震道数据中第二采样时间对应的振幅值。
根据第一振幅值和第二振幅值,可以采用下述公式计算第一时差量对应的第一时窗位置处的时间样点的第一方差值:
公式中,第一振幅值i表示第一地震道数据中第一时窗位置处的第i个时间样点的第一振幅值,第二振幅值i为第二地震道数据中第一时窗位置处的第i个时间样点的第二振幅值。第一方差值可以为第一时差量范围-4ms~4ms内不同时差量对应的所述计算的多个方差值。这个步骤可以通过所述计算的方差值体现不同时差量的真实性,从而为后续步骤中确定准确的时差量做准备。
步骤S104:根据所述第一方差值确定最小方差值,将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗位置处的时间样点对应的目标时差量。
具体地,根据所述计算的每一个第一时差量对应的第一时窗位置处的时间样点的第一方差值,可以比较所述第一方差值的大小,可以确定最小方差值。最小方差值可以表明其对应的时差量最能真实体现第二地震道数据中第一时窗的中间位置处的时间样点的时差。可以将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗的中间位置处的时间样点对应的目标时差量。其中,选取第一时窗位置处的时间样点个数为奇数。第二时差量为最小方差值对应的时差量。
可以将第一时窗的位置移动一个时间样点。采用上述步骤可以计算移动后的第一时窗的中间位置处的时间样点对应的目标时差量,直至完成所述监测地震数据中第二地震道中每一个时间样点对应的目标时差量。从而可以完成所述监测地震数据中每一个地震道中每一个时间样点对应的目标时差量。不仅能够得到各类与采集相关的较大的时差量,也可以得到上覆储层流体引起的局部较小的时差量。
例如,图3中(a)图为图2中(a)图中的基础地震数据和图2中(b)图中的监测地震数据之间的目标时差量剖面图。如图3中(a)图所示,可以看出本发明不仅可以准确计算人为设置的2ms系统时差,也可以准确得到上覆储层流体变化引入的局部时差。
在另一个实施方式中,还可以检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适。
具体地,可以获取所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量。可以为所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量构成的曲线设定一平滑阈值。当所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时量差构成的曲线的平滑度大于平滑阈值时,第一时窗的窗长的选取可能不合适,可以重新选择所述监测地震数据中第二时窗位置的第二时窗的窗长。所述平滑阈值的取值范围可以为大于或等于2毫秒。第一主测线和第一联络测线可以分别为所述监测地震数据中任意主测线和任意联络测线。所述监测地震数据可以包括三个维度,分别为:主测线维度、联络测线维度和时间维度;其中,所述时间维度、主测线维度和所述联络测线维度两两相互垂直。
根据所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量可以对所述监测地震数据进行时差校正。将所述基础地震数据与所述监测地震数据相减的地震数据在第一主测线或第一联络测线上的非储层变化引起的校正前地震响应,可以同所述基础地震数据与所述经过时差校正后的监测地震数据相减的地震数据在第一主测线或第一联络测线上的非储层变化引起的校正后地震响应对比。可以为所述校正后地震响应对比所述校正前地震响应的减弱程度设定一地震响应减弱阈值。当所述校正后地震响应对比所述校正前地震响应的地震响应减弱度小于所述地震响应减弱阈值时,第一时差量范围中最大第一时差量、第一时差量范围中最小第一时差量或第一时差增量的选取可能不合适,可以重新选择所述监测地震数据中第二时差量范围中最大第二时差量、第二时差量范围中最小第二时差量或第二时差增量。所述地震响应减弱阈值的取值范围可以为大于或等于90%。
步骤S105:根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据。
具体地,根据所述监测地震数据中每一个地震道中每一个时间样点对应的目标时差量,可以对所述监测地震数据进行时差校正。可以得到校正地震数据。可以对所述校正地震数据进行辛格函数(Sinc)插值,可以得到目标地震数据。Sinc函数是一个简单的最小相位的子波序列;由于所述时差校正后的校正地震数据采样间隔小于所述监测地震数据采样间隔且采样时间序列可能不连续,可以采用Sinc函数的八点截断实现所述校正地震数据的垂向插值。所述Sinc函数可以采用下述公式表示:
公式中,t表示为所述校正地震数据中的采样时间,Sinc(t)表示为在所述校正地震数据中的采样时间t处插入的振幅值。
例如,图3中(b)图为图2中(a)图中的基础地震数据和2中(b)图中的目标地震数据的振幅差剖面图,图3中(b)图的箭头所示为实际储层变化的期望振幅差异剖面。可以与图2中(c)图中的基础地震数据和监测地震数据的振幅差比较。如图3中(b)图所示,可以看出由所述采集脚印相关时差和上覆储层流体引起的时差导致的非储层振幅差异可以得到消除,实际储层变化的期望振幅差异可以得到凸显。
图4是本申请实施例中基础地震数据和监测地震数据的振幅差在目的储层顶的沿层切片显示的示意图。图5是本申请实施例中基础地震数据和目标地震数据的振幅差在目的储层顶的沿层切片显示的示意图。图4和图5的椭圆形圈内区域部分包含了井眼处实际储层变化引起的振幅差。如图4和图5所示,可以看出由所述采集脚印相关时差和上覆储层流体引起的时差导致的非储层振幅差异可以得到有效消除。椭圆形圈内区域部分所示的井眼处实际储层变化引起的振幅差异更为突出。
所述地震数据时差校正的方法实施例,基于小时窗方差的方式获取并校正基础地震数据和监督地震数据之间时差量,不仅能够校正各类与采集相关的较大的时差量,也可以校正上覆储层流体引起的局部较小的时差量,消除了非储层变化引起的非期望振幅差异,显著提高了两期地震期望振幅差异精度。
图6是本申请一种地震数据时差校正的系统实施例的组成结构图。如图6所示,所述地震数据时差校正的系统可以包括:地震数据获取单元100、时间样点获取单元200、方差值计算单元300、目标时差量获取单元400和目标地震数据获取单元500。
所述地震数据获取单元100,可以用于获取目标区域第一时间对应的基础地震数据,和第二时间对应的监测地震数据。
所述时间样点获取单元200,可以用于分别获取基础地震数据和监测地震数据中相同位置对应的第一地震道数据和第二地震道数据。可以分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点。
所述方差值计算单元300,可以用于预设多个第一时差量,可以计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值。
所述目标时差量获取单元400,可以用于根据所述第一方差值确定最小方差值,将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗位置处的时间样点对应的目标时差量。
所述目标地震数据获取单元500,可以用于根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据。
在另一种实施方式中,所述地震数据时差校正的系统还可以包括:检验单元。
所述检验单元,可以用于检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适。
图7是本申请地震数据时差校正的系统实施例中方差值计算单元的组成结构图。如图7所示,图6中方差值计算单元300包括:第一振幅值获取单元310、第二振幅值获取单元320和第一方差值计算单元330。
所述第一振幅值获取单元310,可以用于获取第一地震道数据中第一时窗位置处的时间样点的第一采样时间对应的第一地震道数据的第一振幅值。
所述第二振幅值获取单元320,可以用于根据第二地震道数据中第一时窗位置处的时间样点的第一采样时间减去第一时差量,得到第二地震道数据中第一时窗位置处的第二采样时间。可以获取第二采样时间对应的第二地震道数据的第二振幅值。
所述第一方差值计算单元330,可以用于根据第一振幅值和第二振幅值,可以计算第一时差量对应的第一时窗位置处的时间样点的第一方差值。
图8是本申请地震数据时差校正的系统实施例中目标地震数据获取单元的组成结构图。如图8所示,图6中目标地震数据获取单元500包括:校正地震数据获取单元510和函数插值单元520。
所述校正地震数据获取单元510,可以用于根据所述目标时差量对所述监测地震数据进行时差校正,得到校正地震数据。
所述函数插值单元520,可以用于对所述校正地震数据进行辛格函数插值,得到目标地震数据。
所述地震数据时差校正的系统实施例与所述地震数据时差校正的方法实施例相对应,可以实现校正各类与采集相关的较大的时差量,也可以实现校正上覆储层流体引起的局部较小的时差量,消除了非储层变化引起的非期望振幅差异,显著提高了两期地震期望振幅差异精度。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(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)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(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),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (12)
1.一种地震数据时差校正的方法,其特征在于,包括:
获取目标区域第一时间对应的基础地震数据,和第二时间对应的监测地震数据;
分别获取基础地震数据和监测地震数据中相同位置对应的第一地震道数据和第二地震道数据,分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点;
预设多个第一时差量,计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值;
根据所述第一方差值确定最小方差值,将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗位置处的时间样点对应的目标时差量;
根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据。
2.所述权利要求1所述一种地震数据时差校正的方法,其特征在于,所述计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值,包括:
获取第一地震道数据中第一时窗位置处的时间样点的第一采样时间对应的第一地震道数据的第一振幅值;
根据第二地震道数据中第一时窗位置处的时间样点的第一采样时间减去第一时差量,得到第二地震道数据中第一时窗位置处的第二采样时间,获取第二采样时间对应的第二地震道数据的第二振幅值;
根据第一振幅值和第二振幅值,计算第一时差量对应的第一时窗位置处的时间样点的第一方差值。
3.根据权利要求2所述一种地震数据时差校正的方法,其特征在于,所述根据第一振幅值和第二振幅值,计算第一时差量对应的第一时窗位置处的时间样点的第一方差值,采用下述公式实现:
公式中,第一振幅值i表示第一地震道数据中第一时窗位置处的第i个时间样点的第一振幅值,第二振幅值i为第二地震道数据中第一时窗位置处的第i个时间样点的第二振幅值。
4.根据权利要求1所述一种地震数据时差校正的方法,其特征在于,所述预设多个第一时差量,包括:
选取第一时差量的第一时差量范围和第一时差增量;
根据第一时差量范围和第一时差增量,得到多个第一时差量。
5.根据权利要求4所述一种地震数据时差校正的方法,其特征在于,第一时差增量为所述监测地震数据的采样间隔的四分之一。
6.根据权利要求1所述一种地震数据时差校正的方法,其特征在于,所述方法还包括:检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适。
7.根据权利要求6所述一种地震数据时差校正的方法,其特征在于,所述检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适,包括:
获取所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量,设定所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量构成的曲线的平滑阈值;
当所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时量差构成的曲线的平滑度大于平滑阈值,第一时窗的窗长的选取不合适;
根据所述监测地震数据中第一主测线或第一联络测线上的时间样点的目标时差量对所述监测地震数据进行时差校正;
将所述基础地震数据与所述监测地震数据相减的地震数据在第一主测线或第一联络测线上的非储层变化引起的校正前地震响应,同所述基础地震数据与所述经过时差校正后的所述监测地震数据相减的地震数据在第一主测线或第一联络测线上的非储层变化引起的校正后地震响应对比,设定所述校正后地震响应对比所述校正前地震响应的地震响应减弱阈值;
当所述校正后地震响应对比所述校正前地震响应的地震响应减弱度小于所述地震响应减弱阈值,第一时差量范围中最大第一时差量、第一时差量范围中最小第一时差量或第一时差增量的选取不合适。
8.根据权利要求1所述一种地震数据时差校正的方法,其特征在于,所述根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据,包括:
根据所述目标时差量对所述监测地震数据进行时差校正,得到校正地震数据;
对所述校正地震数据进行辛格函数插值,得到目标地震数据。
9.一种地震数据时差校正的系统,其特征在于,包括:地震数据获取单元、时间样点获取单元、方差值计算单元、目标时差量获取单元和目标地震数据获取单元;其中,
所述地震数据获取单元,用于获取目标区域第一时间对应的基础地震数据,和第二时间对应的监测地震数据;
所述时间样点获取单元,用于分别获取基础地震数据和监测地震数据中相同位置对应的第一地震道数据和第二地震道数据,分别选取第一地震道数据和第二地震道数据中第一时窗位置处的时间样点;
所述方差值计算单元,用于预设多个第一时差量,计算与每个所述第一时差量分别对应的第一时窗位置处的时间样点的第一方差值;
所述目标时差量获取单元,用于根据所述第一方差值确定最小方差值,将所述最小方差值对应的第二时差量作为第二地震道数据中第一时窗位置处的时间样点对应的目标时差量;
所述目标地震数据获取单元,用于根据所述目标时差量对所述监测地震数据进行时差校正,得到目标地震数据。
10.根据权利要求9所述一种地震数据时差校正的系统,其特征在于,所述方差值计算单元包括:第一振幅值获取单元、第二振幅值获取单元和第一方差值计算单元;其中,
所述第一振幅值获取单元,用于获取第一地震道数据中第一时窗位置处的时间样点的第一采样时间对应的第一地震道数据的第一振幅值;
所述第二振幅值获取单元,用于根据第二地震道数据中第一时窗位置处的时间样点的第一采样时间减去第一时差量,得到第二地震道数据中第一时窗位置处的第二采样时间,获取第二采样时间对应的第二地震道数据的第二振幅值;
所述第一方差值计算单元,用于根据第一振幅值和第二振幅值,计算第一时差量对应的第一时窗位置处的时间样点的第一方差值。
11.根据权利要求9所述一种地震数据时差校正的系统,其特征在于,所述系统还包括:检验单元;其中,
所述检验单元,用于检验第一时窗的窗长、第一时差量范围中最大第一时差量和最小第一时差量,以及第一时差增量的选取是否合适。
12.根据权利要求9所述一种地震数据时差校正的系统,其特征在于,所述目标地震数据获取单元包括:校正地震数据获取单元和函数插值单元;其中,
所述校正地震数据获取单元,用于根据所述目标时差量对所述监测地震数据进行时差校正,得到校正地震数据;
所述函数插值单元,用于对所述校正地震数据进行辛格函数插值,得到目标地震数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610663124.7A CN106353818B (zh) | 2016-08-12 | 2016-08-12 | 一种地震数据时差校正的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610663124.7A CN106353818B (zh) | 2016-08-12 | 2016-08-12 | 一种地震数据时差校正的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106353818A true CN106353818A (zh) | 2017-01-25 |
CN106353818B CN106353818B (zh) | 2018-10-16 |
Family
ID=57844827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610663124.7A Active CN106353818B (zh) | 2016-08-12 | 2016-08-12 | 一种地震数据时差校正的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106353818B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107102358A (zh) * | 2017-05-08 | 2017-08-29 | 中国石油集团东方地球物理勘探有限责任公司 | 基于地震数据的压制气泡干扰的方法和装置 |
CN108375789A (zh) * | 2018-02-11 | 2018-08-07 | 中国石油化工股份有限公司 | 联合采集地震数据的同步匹配方法 |
CN108426689A (zh) * | 2017-12-02 | 2018-08-21 | 青岛理工大学 | 基于抗震规范的地震动选择方法 |
CN109655912A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 三维地震数据动态拉伸时差校正方法及系统 |
CN111983684A (zh) * | 2019-05-24 | 2020-11-24 | 中国石油天然气集团有限公司 | 四维地震数据的时间校准方法及系统 |
CN113740908A (zh) * | 2020-05-29 | 2021-12-03 | 中国石油化工股份有限公司 | 地震切片的二维变差分析方法、电子设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8483964B2 (en) * | 2009-06-02 | 2013-07-09 | Exxonmobil Upstream Research Company | Estimating reservoir properties from 4D seismic data |
CN103744116A (zh) * | 2014-01-23 | 2014-04-23 | 成都晶石石油科技有限公司 | 一种叠前道集的全时间域相位一致性校正方法 |
CN104570116A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 基于地质标志层的时差分析校正方法 |
CN104914471A (zh) * | 2015-05-25 | 2015-09-16 | 中国石油天然气股份有限公司 | 适于黄土塬非纵测线的地滚波压制方法 |
-
2016
- 2016-08-12 CN CN201610663124.7A patent/CN106353818B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8483964B2 (en) * | 2009-06-02 | 2013-07-09 | Exxonmobil Upstream Research Company | Estimating reservoir properties from 4D seismic data |
CN104570116A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 基于地质标志层的时差分析校正方法 |
CN103744116A (zh) * | 2014-01-23 | 2014-04-23 | 成都晶石石油科技有限公司 | 一种叠前道集的全时间域相位一致性校正方法 |
CN104914471A (zh) * | 2015-05-25 | 2015-09-16 | 中国石油天然气股份有限公司 | 适于黄土塬非纵测线的地滚波压制方法 |
Non-Patent Citations (2)
Title |
---|
B. PAZETTI EL.: ""The Impact of Time-shift Estimation and"", 《78TH EAGE CONFERENCE & EXHIBITION 2016》 * |
张汛汛 等: ""基于快速匹配追踪算法的地震道集剩余时差校正"", 《石油物探》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107102358A (zh) * | 2017-05-08 | 2017-08-29 | 中国石油集团东方地球物理勘探有限责任公司 | 基于地震数据的压制气泡干扰的方法和装置 |
CN107102358B (zh) * | 2017-05-08 | 2019-05-14 | 中国石油集团东方地球物理勘探有限责任公司 | 基于地震数据的压制气泡干扰的方法和装置 |
CN109655912A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 三维地震数据动态拉伸时差校正方法及系统 |
CN108426689A (zh) * | 2017-12-02 | 2018-08-21 | 青岛理工大学 | 基于抗震规范的地震动选择方法 |
CN108375789A (zh) * | 2018-02-11 | 2018-08-07 | 中国石油化工股份有限公司 | 联合采集地震数据的同步匹配方法 |
CN108375789B (zh) * | 2018-02-11 | 2020-10-09 | 中国石油化工股份有限公司 | 联合采集地震数据的同步匹配方法 |
CN111983684A (zh) * | 2019-05-24 | 2020-11-24 | 中国石油天然气集团有限公司 | 四维地震数据的时间校准方法及系统 |
CN111983684B (zh) * | 2019-05-24 | 2023-09-26 | 中国石油天然气集团有限公司 | 四维地震数据的时间校准方法及系统 |
CN113740908A (zh) * | 2020-05-29 | 2021-12-03 | 中国石油化工股份有限公司 | 地震切片的二维变差分析方法、电子设备及介质 |
CN113740908B (zh) * | 2020-05-29 | 2024-05-07 | 中国石油化工股份有限公司 | 地震切片的二维变差分析方法、电子设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106353818B (zh) | 2018-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353818A (zh) | 一种地震数据时差校正的方法及系统 | |
CN109100796A (zh) | 一种天然气水合物地震数据处理方法及装置 | |
CN106054251B (zh) | 一种初至波拾取方法及装置 | |
CN106646611B (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN106772589B (zh) | 一种叠前地震反演方法及装置 | |
CN108181649B (zh) | 一种介质分界面振幅补偿方法及装置 | |
CN106338763A (zh) | 一种确定地震剖面显示数据的方法及装置 | |
CN107092031A (zh) | 一种声波测井数据校正方法及装置 | |
CN109061734A (zh) | 一种地震数据的处理方法及装置 | |
CN107179546B (zh) | 一种确定速度场的方法及装置 | |
CN108897036A (zh) | 一种地震数据处理方法及装置 | |
CN108469633B (zh) | 一种地层品质因子的计算方法及装置 | |
CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
CN106291691A (zh) | 一种地震偏移成像方法及装置 | |
CN109100797A (zh) | 一种确定致密储层含气性的方法及装置 | |
CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
CN106908842B (zh) | 一种确定单井沉积微相曲线的方法及装置 | |
CN108181656B (zh) | 一种近偏移距转换波动校正方法及装置 | |
CN107329170A (zh) | 一种地震波形匹配方法及装置 | |
CN106772596A (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN108181643B (zh) | 一种地震勘探数据采集处理方法及装置 | |
CN107942387B (zh) | 一种续至波衰减处理方法及装置 | |
CN108107470B (zh) | 一种地震数据处理方法及装置 | |
CN106249298B (zh) | 一种微震数据噪声压制方法及系统 | |
CN107290783A (zh) | 一种油气检测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |