CN112198550B - 一种时移地震数据可重复性度量方法及装置 - Google Patents
一种时移地震数据可重复性度量方法及装置 Download PDFInfo
- Publication number
- CN112198550B CN112198550B CN201910608751.4A CN201910608751A CN112198550B CN 112198550 B CN112198550 B CN 112198550B CN 201910608751 A CN201910608751 A CN 201910608751A CN 112198550 B CN112198550 B CN 112198550B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- information
- monitoring
- mean square
- root mean
- 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
- 238000000691 measurement method Methods 0.000 title abstract description 7
- 238000012544 monitoring process Methods 0.000 claims abstract description 103
- 238000005314 correlation function Methods 0.000 claims abstract description 74
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 47
- 238000005259 measurement Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000004069 differentiation Effects 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 17
- 238000012545 processing Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
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. for interpretation or for event detection
- G01V1/30—Analysis
-
- 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. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
技术领域
本发明涉及油气勘探技术领域,尤其涉及一种时移地震数据可重复性度量方法及装置。
背景技术
随着油气的开采与开发,地下储层中油气也会发生运移和流动,时移地震可以精确地描述储层流体的变化情况,而时移地震的可重复性高低是时移地震成功与否的关键。因此,提供一种时移地震数据可重复性度量方法是十分必要的。
时移地震至少涉及到两个3D数据:开始数据和监控数据。开始数据是第一次项目勘探采集的地震数据,监控数据是经过一段时间之后采集的地震数据,而差数据就是开始数据与监控数据之差。
现有技术通常采用归一化均方根振幅差和可预测性两个指标进行定量的度量时移地震数据的可重复性。而由于这两个指标具有局限性,因此,现有技术难以实现对时移地震数据可重复性的直接度量。
发明内容
本发明实施例提供一种时移地震数据可重复性度量方法,用以准确地实现对时移地震数据可重复性的直接度量,该方法包括:
获取开始地震数据和监控地震数据;
根据开始地震数据和监控地震数据获取开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息;
利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息;
利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息;
其中,开始地震数据和监控地震数据相关函数的差异化信息定义为:相关函数差异化特征因子;归一化均方根振幅差异信息定义为:归一化均方根振幅差异化特征因子;可重复性时移地震数据度量信息定义为:可重复性度量特征因子,为归一化均方根振幅差异化特征因子和相关函数差异化特征因子的一种组合。
可选的,所述方法还包括:将可重复性时移地震数据度量信息值与预设截断值进行比较,判断时移地震数据的可重复性。
可选的,在获取开始地震数据和监控地震数据后,对开始地震数据和监控地震数据进行预处理。
可选的,利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取归一化均方根振幅差异信息,包括:
利用开始地震数据和监控地震数据获取差数据;
根据开始地震数据、监控地震数据和差数据分别获取开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息;
利用开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息获取归一化均方根振幅差异信息。
本发明实施例还提供一种时移地震数据可重复性度量装置,用以准确地实现对时移地震数据可重复性的直接度量,该装置包括:
数据获取模块,用于获取开始地震数据和监控地震数据;
函数信息获取模块,用于根据开始地震数据和监控地震数据获取开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息;
差异信息获取模块,利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息;
度量信息获取模块,用于利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息;
其中,开始地震数据和监控地震数据相关函数的差异化信息定义为:相关函数差异化特征因子;归一化均方根振幅差异信息定义为:归一化均方根振幅差异化特征因子;可重复性时移地震数据度量信息定义为:可重复性度量特征因子,为归一化均方根振幅差异化特征因子和相关函数差异化特征因子的一种组合。
可选的,所述装置还包括:判断模块,用于将可重复性时移地震数据度量信息值与预设截断值进行比较,判断时移地震数据的可重复性。
可选的,数据获取模块进一步用于:在获取开始地震数据和监控地震数据后,对开始地震数据和监控地震数据进行预处理。
可选的,差异信息获取模块进一步用于:
利用开始地震数据和监控地震数据获取差数据;
根据开始地震数据、监控地震数据和差数据分别获取开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息;
利用开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息获取归一化均方根振幅差异信息。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明实施例提供的时移地震数据可重复性度量方法及装置,通过利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息,进而利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息。由于开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息对数据的极小变化都非常敏感,因此,利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息,可以准确地实现对时移地震数据可重复性的直接度量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中时移地震数据可重复性度量方法的流程示意图;
图2为本发明实施例中时移地震数据可重复性度量装置的结构示意图;
图3为本发明实施例中开始地震数据的具体事例图;
图4为本发明实施例中监控地震数据的具体事例图;
图5为本发明实施例中相关函数差异化信息的具体事例图;
图6为本发明实施例中归一化均方根振幅差异信息的具体事例图;
图7为本发明实施例中可重复性度量特征信息的具体事例图;
图8为本发明实施例中可重复性时移地震数据度量信息与截断值的第一交汇对比示例图;
图9为本发明实施例中可重复性时移地震数据度量信息与截断值的第二交汇对比示例图;
图10为本发明实施例中可重复性时移地震数据度量信息与截断值的第一对比示例图;
图11为本发明实施例中可重复性时移地震数据度量信息与截断值的第二对比示例图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
随着油气的开采与开发,地下储层中油气也会发生运移和流动,时移地震可以精确地描述储层流体的变化情况,而时移地震的可重复性高低是时移地震成功与否的关键。时移地震至少涉及到两个3D数据:开始数据和监控数据。开始数据是第一次项目勘探采集的地震数据,监控数据是经过一段时间之后采集的地震数据,而差数据就是开始数据与监控数据之差。
附图1是本发明实施例提供的一种时移地震数据可重复性度量方法的流程示意图,该方法包括:
步骤101、获取开始地震数据和监控地震数据。
在本发明实施例中,可以用地震人工震源激发和采集开始地震数据和监控地震数据。而为了保证后续作业的顺利进行,在获取开始地震数据和监控地震数据后,可以对开始地震数据和监控地震数据进行预处理。
该预处理包括:设置标签、定义观测系统、速度分析、叠加等。
步骤102、根据开始地震数据和监控地震数据获取开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息。
在本发明实施例中,开始地震数据的自相关函数信息指的是:开始地震数据的平均自相关函数ai(k),其按照以下公式确定:
监控地震数据的自相关函数信息指的是:监控地震数据的平均自相关函数ci(k),其按照以下公式确定:
开始地震数据和监控地震数据的互相关函数信息指的是:开始地震数据与监控地震数据的平均互相关函数bi(k),其按照以下公式确定:
上述公式中,i是时窗参数道的顺序号,i=1,2,…,II,II是时窗参数总道数;k是相关函数延迟顺序号,k=0,±1,±2,…,±KK,(2KK+1)是相关函数长度,相关函数包括自相关函数和互相关函数,j是时窗参数时间样点的顺序号,j=1,2,…,JJ,JJ是时窗参数时间样点数。
步骤103、利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息。
基于上述,在本发明实施例中,开始地震数据和监控地震数据相关函数的差异化信息指的是:相关函数差异化特征因子Predi,其可以表示为:
其中,BBi为互相关函数乘积和特征因子,其可以表示为:
ACi为自相关函数乘积和特征因子,其可以表示为:
使用百分比表示时,相关函数差异化特征因子数值位于0-100%。实际数值可能大于100%,从可知,互相关函数乘积和特征因子是平方和,和的每个元素都是正数,相加只能产生加强作用;而从方程/>可知,自相关函数乘积和特征因子是乘积和,和的一部分元素是正数,一部分元素是负数,相加可能产生抵消作用,这样可能使得互相关函数乘积和特征因子大于自相关函数乘积和特征因子,此时相关函数差异化特征因子数值会大于100%。相关函数差异化特征因子对相关时窗长度和相关函数延迟个数,非常敏感,因此绝对数值大小没有意义。相关函数差异化特征因子对平均静校正、相位或者振幅差异并不敏感,而对噪声和地下反射系数变化非常敏感。如果开始地震数据与监控地震数据不相关,相关函数差异化特征因子是零。如果开始地震数据与监控地震数据反相关,相关函数差异化特征因子是100%。如果一个数据是另外一个数据振幅的一半,相关函数差异化特征因子是100%。
在本发明实施例中,利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取归一化均方根振幅差异信息,包括:
利用开始地震数据和监控地震数据获取差数据;
根据开始地震数据、监控地震数据和差数据分别获取开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息;
利用开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息获取归一化均方根振幅差异信息。
基于上述,开始地震数据均方根振幅特征信息,即开始数据均方根振幅特征因子BRMSi,可以表示为
监控数据均方根振幅特征信息,即监控数据均方根振幅特征因子MRMSi,可以表示为
差数据均方根振幅特征信息,即差数据均方根振幅特征因子DRMSi,可以表示为
归一化均方根振幅差异化特征信息,即归一化均方根振幅差异化特征因子NRMSDi,可以表示为
归一化均方根振幅差异化特征因子数值不太直观,并且没有把范围限制在0-100%。理论上,其最大数值范围为0-200%。例如,如果开始地震数据和监控地震数据中包含随机噪声,其数值为141%如果开始地震数据和监控地震数据反相关(即180度相位差,或者一个数据全部值是0),其数值是200%,理论最大值。如果一个道是另外一个道振幅的一半,其数值是66.7%。归一化均方根振幅差异化特征因子对数据的极小变化都非常敏感。
步骤104、利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息。
在本发明实施例中,定义可重复性时移地震数据度量信息,即可重复性度量特征因子,作为定量的度量时移地震数据可重复性的一个指标。把可重复性度量特征因子定义为归一化均方根振幅差异化特征因子和相关函数差异化特征因子的一种组合,即,相关函数差异化特征因子与归一化均方根振幅差异化特征因子同相关函数差异化特征因子与归一化均方根振幅差异化特征因子之比的乘积,可重复性度量特征因子RepMtri可以表示为:
在本发明实施例中,该方法还包括:将可重复性时移地震数据度量信息值与预设截断值进行比较,判断时移地震数据的可重复性。
具体地,在判断过程中,对于给定的可重复性度量特征因子的截断门槛值,把大于或者等于截断门槛值的点,指示可重复性设置为1,否则,设置为0。即,可重复性度量指示特征因子IndRepMtri可以表示为
式中,RepMtrCut是可重复性度量特征因子的截断门槛值。
本发明实施例提供的时移地震数据可重复性度量方法,通过利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息,进而利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息。由于开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息对数据的极小变化都非常敏感,因此,利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息,可以准确地实现对时移地震数据可重复性的直接度量。
下面以一具体实验结果对本发明进行说明:
绘制相关函数差异化特征因子、归一化均方根振幅差异化特征因子、可重复性度量特征因子和可重复性度量指示特征因子数据曲线,并进行数据存储。数据曲线如下:
图3和图4是共炮点道集数据对比,其中,图3是开始地震数据,图4是监控地震数据。
图5是相关函数差异化特征因子示意图。
图6是归一化均方根振幅差异化特征因子示意图。
图7是可重复性度量特征因子示意图。
图8和图9是可重复性度量特征因子与可重复性度量截断门槛值交汇对比示意图,其中,图8是可重复性度量截断门槛值是第1道可重复性度量特征因子的1.1倍,图9是可重复性度量截断门槛值是第1道可重复性度量特征因子的1.2倍.
图10和图11是可重复性度量指示特征因子对比示意图,其中,图10是可重复性度量截断门槛值是第1道可重复性度量特征因子的1.1倍,图11是可重复性度量截断门槛值是第1道可重复性度量特征因子的1.2倍。
图10和图11中可重复性度量指示特征因子数值为1的地震道,其对应时移地震数据是可重复的,可重复性度量指示特征因子数值为0的地震道,其对应时移地震数据没有重复性。由图5、图6、图7、图8、图9、图10、图11可以看出,利用本发明可以有效地度量时移地震数据的可重复性。因此,本发明提供的时移地震数据可重复性度量技术,满足了实际数据处理的需要。
基于同一发明构思,本发明实施例中还提供了一种时移地震数据可重复性度量装置,如下面的实施例所述。由于时移地震数据可重复性度量装置解决问题的原理与时移地震数据可重复性度量方法相似,因此,时移地震数据可重复性度量装置的实施可以参见时移地震数据可重复性度量方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
本发明实施例提供了一种时移地震数据可重复性度量装置,如附图2所示,该装置包括:
数据获取模块201,用于获取开始地震数据和监控地震数据。
函数信息获取模块202,用于根据开始地震数据和监控地震数据获取开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息。
差异信息获取模块203,利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息。
度量信息获取模块204,用于利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息。
在本发明实施例中,如附图2所示,时移地震数据可重复性度量装置还包括:判断模块205,用于将可重复性时移地震数据度量信息值与预设截断值进行比较,判断时移地震数据的可重复性。
在本发明实施例中,数据获取模块201进一步用于:在获取开始地震数据和监控地震数据后,对开始地震数据和监控地震数据进行预处理。
在本发明实施例中,差异信息获取模块203进一步用于:
利用开始地震数据和监控地震数据获取差数据;
根据开始地震数据、监控地震数据和差数据分别获取开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息;
利用开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息获取归一化均方根振幅差异信息。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述方法。
本发明实施例还提供了一种计算机可读存储介质,计算机可读存储介质存储有执行上述方法的计算机程序。
综上,本发明通过对相关函数差异化特征因子、归一化均方根振幅差异化特征因子、可重复性度量特征因子和可重复性度量指示特征因子的计算,实现了对时移地震数据可重复性的直接度量,具有计算量小、计算速度快、稳定性好和计算精度高的特点。
此外,本发明的技术方案也适用于海面双检拖缆、海底电缆/海底节点等海洋采集的水中检波器与陆地检波器地震数据。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种时移地震数据可重复性度量方法,其特征在于,包括:
获取开始地震数据和监控地震数据;
根据开始地震数据和监控地震数据获取开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息;
利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息;
利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息;
其中,开始地震数据和监控地震数据相关函数的差异化信息定义为:相关函数差异化特征因子;归一化均方根振幅差异信息定义为:归一化均方根振幅差异化特征因子;可重复性时移地震数据度量信息定义为:可重复性度量特征因子,为归一化均方根振幅差异化特征因子和相关函数差异化特征因子的一种组合。
2.如权利要求1所述的方法,其特征在于,还包括:将可重复性时移地震数据度量信息值与预设截断值进行比较,判断时移地震数据的可重复性。
3.如权利要求1所述的方法,其特征在于,在获取开始地震数据和监控地震数据后,对开始地震数据和监控地震数据进行预处理。
4.如权利要求1所述的方法,其特征在于,利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取归一化均方根振幅差异信息,包括:
利用开始地震数据和监控地震数据获取差数据;
根据开始地震数据、监控地震数据和差数据分别获取开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息;
利用开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息获取归一化均方根振幅差异信息。
5.一种时移地震数据可重复性度量装置,其特征在于,包括:
数据获取模块,用于获取开始地震数据和监控地震数据;
函数信息获取模块,用于根据开始地震数据和监控地震数据获取开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息;
差异信息获取模块,利用开始地震数据的自相关函数信息、监控地震数据的自相关函数信息、以及开始地震数据和监控地震数据的互相关函数信息,获取开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息;
度量信息获取模块,用于利用开始地震数据和监控地震数据相关函数的差异化信息和归一化均方根振幅差异信息获取可重复性时移地震数据度量信息;
其中,开始地震数据和监控地震数据相关函数的差异化信息定义为:相关函数差异化特征因子;归一化均方根振幅差异信息定义为:归一化均方根振幅差异化特征因子;可重复性时移地震数据度量信息定义为:可重复性度量特征因子,为归一化均方根振幅差异化特征因子和相关函数差异化特征因子的一种组合。
6.如权利要求5所述的装置,其特征在于,还包括:判断模块,用于将可重复性时移地震数据度量信息值与预设截断值进行比较,判断时移地震数据的可重复性。
7.如权利要求5所述的装置,其特征在于,数据获取模块进一步用于:在获取开始地震数据和监控地震数据后,对开始地震数据和监控地震数据进行预处理。
8.如权利要求5所述的装置,其特征在于,差异信息获取模块进一步用于:
利用开始地震数据和监控地震数据获取差数据;
根据开始地震数据、监控地震数据和差数据分别获取开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息;
利用开始地震数据均方根振幅特征信息、监控地震数据均方根振幅特征信息和差数据均方根振幅特征信息获取归一化均方根振幅差异信息。
9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910608751.4A CN112198550B (zh) | 2019-07-08 | 2019-07-08 | 一种时移地震数据可重复性度量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910608751.4A CN112198550B (zh) | 2019-07-08 | 2019-07-08 | 一种时移地震数据可重复性度量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112198550A CN112198550A (zh) | 2021-01-08 |
CN112198550B true CN112198550B (zh) | 2023-09-26 |
Family
ID=74004739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910608751.4A Active CN112198550B (zh) | 2019-07-08 | 2019-07-08 | 一种时移地震数据可重复性度量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112198550B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301635A (zh) * | 2015-10-09 | 2016-02-03 | 中国石油天然气集团公司 | 一种海上拖缆时移地震数据的构建方法和装置 |
CN105549083A (zh) * | 2015-12-25 | 2016-05-04 | 中国石油天然气集团公司 | 一种确定地震观测系统重复性整体量度的方法和装置 |
CN105866837A (zh) * | 2016-04-06 | 2016-08-17 | 中国海洋石油总公司 | 一种海上非重复性时移地震多层系互均衡处理方法 |
CN106383362A (zh) * | 2016-08-19 | 2017-02-08 | 中国海洋石油总公司 | 一种提高薄储层时移地震差异识别能力的方法 |
CN108614296A (zh) * | 2018-06-06 | 2018-10-02 | 中国石油集团东方地球物理勘探有限责任公司 | 观测系统重复性确定方法及装置 |
CN109407152A (zh) * | 2018-12-18 | 2019-03-01 | 吉林大学 | 基于零均值归一化互相关目标函数的时间域全波形反演方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140003191A1 (en) * | 2012-06-28 | 2014-01-02 | Chevron U.S.A. Inc | System and method for evaluating a time-lapse seismic signal recording using shifted normalized root mean square metric |
US20160161620A1 (en) * | 2013-06-28 | 2016-06-09 | Cgg Services Sa | System and method for estimating repeatability using base data |
US10768325B2 (en) * | 2016-10-27 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to estimate 4D seismic acquisition repeatability specifications from high-resolution near-water-bottom seismic images |
-
2019
- 2019-07-08 CN CN201910608751.4A patent/CN112198550B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301635A (zh) * | 2015-10-09 | 2016-02-03 | 中国石油天然气集团公司 | 一种海上拖缆时移地震数据的构建方法和装置 |
CN105549083A (zh) * | 2015-12-25 | 2016-05-04 | 中国石油天然气集团公司 | 一种确定地震观测系统重复性整体量度的方法和装置 |
CN105866837A (zh) * | 2016-04-06 | 2016-08-17 | 中国海洋石油总公司 | 一种海上非重复性时移地震多层系互均衡处理方法 |
CN106383362A (zh) * | 2016-08-19 | 2017-02-08 | 中国海洋石油总公司 | 一种提高薄储层时移地震差异识别能力的方法 |
CN108614296A (zh) * | 2018-06-06 | 2018-10-02 | 中国石油集团东方地球物理勘探有限责任公司 | 观测系统重复性确定方法及装置 |
CN109407152A (zh) * | 2018-12-18 | 2019-03-01 | 吉林大学 | 基于零均值归一化互相关目标函数的时间域全波形反演方法 |
Non-Patent Citations (5)
Title |
---|
ED KRAGH and PHIL CHRISTIE,等.Seismic repeatability, normalized rms, and predictability.THE LEADING EDGE.2002,第640-646页. * |
JUAN CANTILLO, Total E&P.Throwing a new light on time-lapse technology, metrics and 4D repeatability with SDR.The Leading Edge.2012,第405-411页. * |
Peter Gagliardi and Don C. Lawton .Experimental comparison of repeatability metrics .CREWES Research Report.2011,第1-16页. * |
朱振宇等.海上时移地震关键技术研究与应用.中国海上油气.2018,第76-85页. * |
郭念民等.时移地震数据重复性度量及一致性分析方法.物探化探计算技术.2012,第186-193页. * |
Also Published As
Publication number | Publication date |
---|---|
CN112198550A (zh) | 2021-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10386262B2 (en) | Leak identification in a mass transport network | |
Cheng et al. | Experimental strain modal analysis for beam-like structure by using distributed fiber optics and its damage detection | |
CN112946783A (zh) | 一种水合物饱和度确定方法、装置及设备 | |
CN110470581B (zh) | 确定储层应力敏感程度的方法、装置及存储介质 | |
CN112198550B (zh) | 一种时移地震数据可重复性度量方法及装置 | |
CN112049623B (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
CN111241652B (zh) | 一种确定地层原油粘度的方法及装置 | |
CN106415322B (zh) | 识别反射信号的方法 | |
CN113552036A (zh) | 中低孔渗储层油水相对渗透率的确定方法及装置 | |
CN111077574A (zh) | 一种确定地层弹性参数的方法、装置及系统 | |
CN106837318B (zh) | 岩层稠油含量获取方法和装置 | |
CN103498658B (zh) | 一种伽马测井数据处理方法 | |
US10487643B2 (en) | Two-dimensional imaging with multi-stage processing | |
Berlin | Application of the random decrement technique for experimental determination of damping parameters of bearing structures | |
CN112268608A (zh) | 一种强干扰dvs的高精度解调算法 | |
CN105510986B (zh) | 一种确定海绿石砂岩中海绿石含量的方法和装置 | |
CN113050191A (zh) | 基于双参数的页岩油toc预测方法及装置 | |
CN104533405A (zh) | 泥页岩储层有机碳含量的检测方法和装置 | |
CN109597119A (zh) | 加卸载响应比计算方法和系统 | |
CN114019578A (zh) | 一种超导航磁梯度张量重复测线一致性检验方法和装置 | |
Agneni et al. | Modal parameters directly estimated from power spectral densities or correlation functions in output-only analysis | |
CN111856571B (zh) | 测井频段速度向地震频段速度转换方法、装置及设备 | |
CN112630853B (zh) | 基于无阻流量的含气饱和度校正方法和装置 | |
CN112147172B (zh) | 基于核磁共振t2谱评估含水饱和度的方法和装置 | |
CN111781643B (zh) | 地震层速度确定方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |