CN105425294B - 采集脚印的确定方法和装置 - Google Patents
采集脚印的确定方法和装置 Download PDFInfo
- Publication number
- CN105425294B CN105425294B CN201510852245.1A CN201510852245A CN105425294B CN 105425294 B CN105425294 B CN 105425294B CN 201510852245 A CN201510852245 A CN 201510852245A CN 105425294 B CN105425294 B CN 105425294B
- Authority
- CN
- China
- Prior art keywords
- amplitude curve
- auto
- main profile
- cross
- correlation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 239000000284 extract Substances 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 35
- 238000000605 extraction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000011002 quantification Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000004804 winding 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
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
技术领域
本发明涉及地震勘探技术领域,特别涉及一种采集脚印的确定方法和装置。
背景技术
在地震勘探中,尤其是在三维地震勘探中的采集脚印问题,一直是学者和专家关注的问题,所谓采集脚印指的是在地震数据的采集、处理过程中在数据体上留下的痕迹,采集脚印在时间切片和深度切片上一般表现为振幅和相位的变化。因此,由于采集脚印的存在,从而使得在地震成像结果中会出现周期性的振幅假象,使得地震成像结果出现偏差。
目前,判断采集脚印的方法主要是基于从业人员对资料的感性认识和从业经验,是一种定性化的分析方法,没有一种统一的方法来准确衡量不同地区、不同方法资料中所存在的采集脚印的强弱程度。
针对如何定量确定是否存在采集脚印,以及采集脚印的强弱程度,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种采集脚印的确定方法,以达到准确定量化确定采集脚印的分布情况的目的,该方法包括:
在三维数据体上提取目标层位的等时切片;
在所述等时切片上分别抽取主测线方向的振幅曲线和联络测线方向的振幅曲线;
对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理;
根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布。
在一个实施方式中,对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理,包括:
分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理;
对预处理后的主测线方向的振幅曲线进行自相关处理得到主测线方向的自相关振幅曲线,对预处理后的联络测线方向的振幅曲线进行自相关处理得到联络测线方向的自相关振幅曲线。
在一个实施方式中,根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布包括:
根据主测线方向的自相关振幅曲线和联络测线方向的自相关振幅曲线、所述目标层位的炮线距和接收线距,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,以及存在的采集脚印的强度。
在一个实施方式中,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,包括:
将主测线方向的自相关振幅曲线的周期与目标层位的炮线距进行比较,如果变化趋势一致,则确定在主测线方向存在采集脚印;
将联络测线方向的自相关振幅曲线的周期与目标层位的接收线距进行比较,如果变化趋势一致,则确定在联络线方向存在采集脚印。
在一个实施方式中,根据主测线方向的自相关振幅曲线和联络测线方向的自相关振幅曲线、所述目标层位的炮线距和接收线距,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,以及存在的采集脚印的强度包括:
对主测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将主测线方向的自相关振幅曲线的横坐标转换为距离;
对联络测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将联络测线方向的自相关振幅曲线的横坐标转换为距离。
在一个实施方式中,确定存在的采集脚印的强度,包括:
对于存在采集脚印的三维数据体,根据归一化处理后的主测线方向的自相关振幅曲线的第二个波峰的峰值大小确定主测线方向存在的采集脚印的强弱程度,根据归一化处理后的联络测线方向的自相关振幅曲线的第二个波峰的峰值大小确定联络测线方向存在的采集脚印的强弱程度,峰值越大,采集脚印的强度越大,峰值越小,采集脚印的强度越小。
在一个实施方式中,所述峰值的最大值为1。
在一个实施方式中,分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理,包括:
分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理,使得预处理后的主测线方向的振幅曲线的振幅平均值为0,预处理后的联络测线方向的振幅曲线的振幅平均值为0。
本发明实施例还提供了一种采集脚印的确定装置,以达到准确定量化确定采集脚印的分布情况的目的,该装置包括:
提取模块,用于在三维数据体上提取目标层位的等时切片;
抽取模块,用于在所述等时切片上分别抽取主测线方向的振幅曲线和联络测线方向的振幅曲线;
处理模块,用于对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理;
确定模块,用于根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布。
在一个实施方式中,所述处理模块,包括:
预处理单元,用于分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理;
自相关处理单元,用于对预处理后的主测线方向的振幅曲线进行自相关处理得到主测线方向的自相关振幅曲线,对预处理后的联络测线方向的振幅曲线进行自相关处理得到联络测线方向的自相关振幅曲线。
在本发明实施例中,通过从等时切片上抽取的主测线方向的振幅曲线和联络测线方向的振幅曲线,以及目标层位的炮线距和接收线距确定目标层位采集脚印的分布。因为采集脚印的分布会对振幅曲线产生影响,而这种影响与炮线距和接收线距存在一定的一致性关系,因此,通过振幅曲线、炮线距和接收线距可以实现定量化确定采集脚印分布的目的,基于准确确定的采集脚印的分布可以在实际生产中有效消除因为采集脚印的存在对地下信息准确性产生的影响。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是根据本发明实施例的采集脚印的确定方法的方法流程图;
图2是根据本发明实施例的转换后得到的Inline方向的自相关振幅曲线示意图;
图3是根据本发明实施例的转换后得到的Crossline方向的自相关振幅曲线示意图;
图4是根据本发明实施例的采集脚印的确定装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本例中提供了一种采集脚印的确定方法,如图1所示,可以包括以下步骤:
步骤101:在三维数据体上提取目标层位的等时切片;
自从进行三维地震勘探以来,切片数据的应用就成为三维地震数据解释的重要手段,并且从水平时间切片不断发展为沿层切片和地层切片。目前,地层切片技术较为流行,该技术要求在地震剖面上找出地质上等时的或平行等时面的同相轴,在找出的同相轴之间内插做切片,对于这些切片一般都认为是等时切片,且这些等时切片是相应的地质信息的响应。
获取等时切片的作法可以是:先在地震剖面上找出与频率无关的同相轴,即用不同的频带对地震数据进行滤波,那些不随滤波参数变化的同相轴被认为是地质上等时的或平行等时面的同相轴,将这些同相轴选定为参考同相轴。然后,在参考同相轴之间内插做切片,认为这样的切片也都是等时的,即地层切片,并认为其中任何一张时间切片都代表这一地层时间模型中相对地质时间界面的地震响应。
基于上述的等时切片,在本例中,在获取了三维数据体后,可以在对应的目标层位上提取等时切片中,等时切片会存在主测线(Inline)方向的信息和联络测线(Crossline)方向的信息。
步骤102:在所述等时切片上分别抽取主测线方向的振幅曲线和联络测线方向的振幅曲线;
步骤103:对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理;
对主测线方向的振幅曲线和联络测线方向的振幅曲线进行处理可以包括:进行预处理、自相关处理、和归一化处理。
具体地,所谓的预处理可以是:分别对主测线方向的振幅曲线和联络测线方向的振幅曲线进行预处理,使得预处理后的主测线方向的振幅曲线的振幅平均值为0,预处理后的联络测线方向的振幅曲线的振幅平均值为0,即以振幅曲线的振幅平均值为0作为预处理的原则对主测线方向的振幅曲线和联络测线方向的振幅曲线进行处理。
上述预处理的目的主要是减少振幅变化对结果的影响,在预处理后,就可以对预处理后的主测线方向的振幅曲线和联络测线方向的振幅曲线进行自相关分析处理。
在实现的时候,分别对主测线方向的振幅曲线进行自相关分析和对联络测线方向的振幅曲线进行自相关分析,在自相关分析的过程中,所应用的自相关函数就是一般的自相关函数。
为了使得后续比较基础一致,在进行自相关处理后,可以对自相关分析处理得到的自相关振幅曲线进行归一化处理,进一步的,考虑到后续比较时候横坐标需要是距离的参数,因此,可以将归一化后的自相关振幅曲线的横坐标转换为距离。具体地,对主测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将主测线方向的自相关振幅曲线的横坐标转换为距离;对联络测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将联络测线方向的自相关振幅曲线的横坐标转换为距离。
步骤104:根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布。
在步骤103中对主测线方向的振幅曲线和联络测线方向的振幅曲线进行不同的处理,得到的结果自然也不同,那么作为比较依据的数据也是不同的。当然最好的最准确的比较数据就是对原始的主测线方向的振幅曲线和联络测线方向的振幅曲线进行预处理、再进行自相关处理,然后再进行归一化处理,通过最终得到的数据作为最后判断采集脚印的分布的判断依据。
在判断的时候,主要是将得到的主测线方向的自相关振幅曲线的周期和联络测线方向的自相关振幅曲线的周期分别与炮线距和接收线距进行对比,以确定是否存在采集脚印。这种比较方式所依据的原理是:采集脚印的存在会对振幅产生影响,因此,如果主测线方向的振幅曲线的变化趋势与炮线距的变化趋势相同,或者是联络测线方向的振幅曲线与接收线距的变化趋势相同,则说明是存在脚印的,如果变化趋势不一致,也就说明当前的层位未受到原始检测的影响,也就说明不存在采集脚印。
在一个实施方式中,可以将主测线方向的自相关振幅曲线的周期与目标层位的炮线距进行比较,如果变化趋势一致,则确定在主测线方向存在采集脚印;将联络测线方向的自相关振幅曲线的周期与目标层位的接收线距进行比较,如果变化趋势一致,则确定在联络线方向存在采集脚印。其中,所谓的周期在此处是一个距离的概念。
例如:炮线距为:500m,接收线距为:400m,如果如图2所示,主测线方向的振幅曲线的周期也为500m或接近500m(例如:误差可以控制在10%以内)时,那么可以说明在主测线方向是存在采集脚印的;如果如图3所示,联络线方向的振幅曲线的周期也为N*400m(N为大于0的正整数,一般可以取值为1-3)或接近N*400m(例如,误差也可以控制在10%以内)时,那么可以说明在联络线方向是存在采集脚印的。
在确定了是否存在采集脚印之后,还可以对存在的采集脚印的强弱程度进行判断,具体地,可以通过处理后的自相关振幅曲线的第二波峰的峰值大小来确定脚印强弱,例如,对于存在采集脚印的三维数据体,可以根据自相关处理和归一化处理后的振幅曲线的第二个波峰的峰值大小确定采集脚印的强弱程度,峰值越大,采集脚印的强度越大,峰值越小,采集脚印的强度越小。因为进行了归一化处理,因此峰值的最大值为1。
然而值得注意的是,上述仅是说明采集脚印的确定是根据振幅曲线和接收线距以及炮线距进行的,在实际实现的过程中,还需要结合目的层的信噪比等进行更为精准的判断,上述根据振幅曲线和接收线距以及炮线距仅是判断过程的一个重要依据,但不是全部,为了更精确地确定采集脚印的分布,还需要结合实际生产中得到的其它参数信息。
下面结合一个具体的实施例对上述采集脚印的确定方法进行具体说明,然而值得注意的是,该具体实施例仅是为了更好地说明本发明,并不构成对本发明的不当限定。
本例主要是考虑到当前地震勘探的目标一般都是岩性油气藏和隐蔽油气藏,这些对振幅、频率等地震属性的准确度提出了越来越高的要求,这也就要求要尽量保护地下信息的真实性,因为采集脚印的存在会使得地下信息的真实性受到影响,因此,需要有一种可以对采集脚印进行定量研究的方法,以便准确确定采集脚印的分布情况,在清楚确定采集脚印分布情况的前提下,才能更好地还原地下信息的真实性。
具体地,在本例中,针对采集脚印难以准确识别和难以量化判断强弱的问题,提出了一种在三维数据体等时切片上提取Inline振幅曲线、Crossline振幅曲线,并通过一系列的转换,最终达到定量化判断是否存在采集脚印及其强弱的目的。
在实现的时候,可以包括以下几个步骤:
S1:首先在三维数据体上提取目标层位的等时切片;
S2:在等时切片上分别抽取Inline方向振幅曲线、Crossline方向振幅曲线;
S3:分别对Inline方向振幅曲线和Crossline方向振幅曲线进行预处理,其中,预处理的要求是:预处理后的Inline方向振幅曲线和Crossline方向振幅曲线的振幅平均值都为零;
S4:对预处理后的振幅曲线进行自相关分析;
S5:对自相关分析后的结果进行归一化处理,并在归一化处理的过程中,将振幅曲线的横坐标转换为距离其中,距离可以等于点数*面元;
S6:根据归一化处理后的Inline方向振幅曲线的周期与炮线距的关系判断在Inline方向是否具有采集脚印,根据归一化处理后的Crossline方向振幅曲线的周期与接收线距的关系判断在Crossline方向是否具有采集脚印,具体地:如果周期与炮线距和接收线距一致或相近,则可以说明存在采集脚印;如果周期与炮线距和接收线距不相近,则说明不存在采集脚印或采集脚印很弱;
S7:对于存在采集脚印的数据体可根据自相关曲线第二个波峰值的大小判断采集脚印的强弱程度,峰值越大,采集脚印越严重,其最大值为1。
以台盆区为例,经过多年的勘探,该区域已经取得了很大的发展,勘探程度相对较高,勘探重点已经转向岩性勘探,岩性油气藏和隐蔽性油气藏勘探对地震资料的细节刻画要求更加精确,采集脚印的存在影响了地震资料振幅、相位属性,使得后期的储层预测的结果存在一定的误差,因为采集脚印的存在模糊了储层边界和储层构造特征,出现一些储层预测的假象,对油藏评价和储量预测有很大的影响,不利于油田的进一步勘探开发。
中古8高密度三维是目前塔里木盆地台盆区观测系统最强化的三维,通过对叠前时间偏移数据体不同时间的切片分析,发现中古8三维在2.5s处存在采集脚印现象,而深层无采集脚印现象,因此通过对中古8不同观测系统参数(如表1所示)的叠前时间偏移数据体2.5m处的切片进行分析,发现利用已有的方法很难判断那种观测系统采集脚印更强,而采用本例中所提供的方法可以很容易判断出自相关振幅曲线的周期与炮线距相关,观测系统A的采集脚印强。
表1
中古8观测方案A | 中古8观测方案B | |
接收线数 | 22 | 12 |
面元尺寸 | 15×15 | 15×15 |
覆盖次数 | 121 | 132 |
接收线距 | 480 | 240 |
炮线距 | 480 | 240 |
宽度系数 | 1 | 0.27 |
炮道密度 | 53.8 | 58.7 |
基于同一发明构思,本发明实施例中还提供了一种采集脚印的确定装置,如下面的实施例所述。由于采集脚印的确定装置解决问题的原理与采集脚印的确定方法相似,因此采集脚印的确定装置的实施可以参见采集脚印的确定方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图3是本发明实施例的采集脚印的确定装置的一种结构框图,如图4所示,可以包括:提取模块401、抽取模块402、处理模块403和确定模块404,下面对该结构进行说明。
提取模块401,用于在三维数据体上提取目标层位的等时切片;
抽取模块402,用于在所述等时切片上分别抽取主测线方向的振幅曲线和联络测线方向的振幅曲线;
处理模块403,用于对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理;
确定模块404,用于根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布。
在一个实施方式中,处理模块403可以包括:预处理单元,用于分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理;自相关处理单元,用于对预处理后的主测线方向的振幅曲线进行自相关处理得到主测线方向的自相关振幅曲线,对预处理后的联络测线方向的振幅曲线进行自相关处理得到联络测线方向的自相关振幅曲线。
在一个实施方式中,确定模块404可以根据主测线方向的自相关振幅曲线和联络测线方向的自相关振幅曲线、所述目标层位的炮线距和接收线距,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,以及存在的采集脚印的强度。
在一个实施方式中,确定模块404具体可以将主测线方向的自相关振幅曲线的周期与目标层位的炮线距进行比较,如果变化趋势一致,则确定在主测线方向存在采集脚印;将联络测线方向的自相关振幅曲线的周期与目标层位的接收线距进行比较,如果变化趋势一致,则确定在联络线方向存在采集脚印。
在一个实施方式中,处理模块403具体可以对主测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将主测线方向的自相关振幅曲线的横坐标转换为距离;对联络测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将联络测线方向的自相关振幅曲线的横坐标转换为距离。
在一个实施方式中,确定模块404具体可以对存在采集脚印的三维数据体,根据归一化处理后的主测线方向的自相关振幅曲线的第二个波峰的峰值大小确定主测线方向存在的采集脚印的强弱程度,根据归一化处理后的联络测线方向的自相关振幅曲线的第二个波峰的峰值大小确定联络测线方向存在的采集脚印的强弱程度,峰值越大,采集脚印的强度越大,峰值越小,采集脚印的强度越小。
在一个实施方式中,峰值的最大值为1。
在一个实施方式中,预处理单元具体可以分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理,使得预处理后的主测线方向的振幅曲线的振幅平均值为0,预处理后的联络测线方向的振幅曲线的振幅平均值为0。
在另外一个实施例中,还提供了一种软件,该软件用于执行上述实施例及优选实施方式中描述的技术方案。
在另外一个实施例中,还提供了一种存储介质,该存储介质中存储有上述软件,该存储介质包括但不限于:光盘、软盘、硬盘、可擦写存储器等。
从以上的描述中,可以看出,本发明实施例实现了如下技术效果:通过从等时切片上抽取的主测线方向的振幅曲线和联络测线方向的振幅曲线,以及目标层位的炮线距和接收线距确定目标层位采集脚印的分布。因为采集脚印的分布会对振幅曲线产生影响,而这种影响与炮线距和接收线距存在一定的一致性关系,因此,通过振幅曲线和炮线距和接收线距可以实现定量化确定采集脚印分布的目的,基于准确确定的采集脚印的分布可以在实际生产中有效消除因为采集脚印的存在对地下信息准确性产生的影响。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种采集脚印的确定方法,其特征在于,包括:
在三维数据体上提取目标层位的等时切片;
在所述等时切片上分别抽取主测线方向的振幅曲线和联络测线方向的振幅曲线;
对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理;
根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布;
其中,对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理,包括:
分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理;
对预处理后的主测线方向的振幅曲线进行自相关处理得到主测线方向的自相关振幅曲线,对预处理后的联络测线方向的振幅曲线进行自相关处理得到联络测线方向的自相关振幅曲线。
2.如权利要求1所述的方法,其特征在于,根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布包括:
根据主测线方向的自相关振幅曲线和联络测线方向的自相关振幅曲线、所述目标层位的炮线距和接收线距,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,以及存在的采集脚印的强度。
3.如权利要求2所述的方法,其特征在于,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,包括:
将主测线方向的自相关振幅曲线的周期与目标层位的炮线距进行比较,如果变化趋势一致,则确定在主测线方向存在采集脚印;
将联络测线方向的自相关振幅曲线的周期与目标层位的接收线距进行比较,如果变化趋势一致,则确定在联络线方向存在采集脚印。
4.如权利要求2所述的方法,其特征在于,根据主测线方向的自相关振幅曲线和联络测线方向的自相关振幅曲线、所述目标层位的炮线距和接收线距,确定所述目标层位的主测线方向和联络测线方向是否存在采集脚印,以及存在的采集脚印的强度包括:
对主测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将主测线方向的自相关振幅曲线的横坐标转换为距离;
对联络测线方向的自相关振幅曲线进行归一化处理,且在归一化处理的过程中,将联络测线方向的自相关振幅曲线的横坐标转换为距离。
5.如权利要求4所述的方法,其特征在于,确定存在的采集脚印的强度,包括:
对于存在采集脚印的三维数据体,根据归一化处理后的主测线方向的自相关振幅曲线的第二个波峰的峰值大小确定主测线方向存在的采集脚印的强弱程度,根据归一化处理后的联络测线方向的自相关振幅曲线的第二个波峰的峰值大小确定联络测线方向存在的采集脚印的强弱程度,峰值越大,采集脚印的强度越大,峰值越小,采集脚印的强度越小。
6.如权利要求5所述的方法,其特征在于,所述峰值的最大值为1。
7.如权利要求1至6中任一项所述的方法,其特征在于,分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理,包括:
分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理,使得预处理后的主测线方向的振幅曲线的振幅平均值为0,预处理后的联络测线方向的振幅曲线的振幅平均值为0。
8.一种采集脚印的确定装置,其特征在于,包括:
提取模块,用于在三维数据体上提取目标层位的等时切片;
抽取模块,用于在所述等时切片上分别抽取主测线方向的振幅曲线和联络测线方向的振幅曲线;
处理模块,用于对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行处理;
确定模块,用于根据所述目标层位的炮线距和接收线距、以及处理结果,确定所述目标层位采集脚印的分布;
其中,所述处理模块,包括:
预处理单元,用于分别对所述主测线方向的振幅曲线和所述联络测线方向的振幅曲线进行预处理;
自相关处理单元,用于对预处理后的主测线方向的振幅曲线进行自相关处理得到主测线方向的自相关振幅曲线,对预处理后的联络测线方向的振幅曲线进行自相关处理得到联络测线方向的自相关振幅曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510852245.1A CN105425294B (zh) | 2015-11-30 | 2015-11-30 | 采集脚印的确定方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510852245.1A CN105425294B (zh) | 2015-11-30 | 2015-11-30 | 采集脚印的确定方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105425294A CN105425294A (zh) | 2016-03-23 |
CN105425294B true CN105425294B (zh) | 2018-01-05 |
Family
ID=55503618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510852245.1A Active CN105425294B (zh) | 2015-11-30 | 2015-11-30 | 采集脚印的确定方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105425294B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353816B (zh) * | 2016-08-09 | 2018-10-16 | 中国石油天然气集团公司 | 一种地震采集脚印噪音压制方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173989A (zh) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | 地震采集三维观测系统定量分析方法 |
CN101852866A (zh) * | 2010-05-13 | 2010-10-06 | 中国石油天然气股份有限公司 | 一种叠后地震数据滤波方法 |
CN102023312A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统采集脚印定量分析方法 |
CN102645670A (zh) * | 2011-02-22 | 2012-08-22 | 中国石油天然气集团公司 | 一种基于叠加响应分析的观测系统优化设计方法 |
CN103558636A (zh) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种从叠后地震数据采集脚印衰减的方法 |
CN105044799A (zh) * | 2015-07-10 | 2015-11-11 | 中国石油大学(北京) | 确定三维地震观测系统面元属性均匀程度及均匀化的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9645268B2 (en) * | 2012-06-25 | 2017-05-09 | Schlumberger Technology Corporation | Seismic orthogonal decomposition attribute |
-
2015
- 2015-11-30 CN CN201510852245.1A patent/CN105425294B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173989A (zh) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | 地震采集三维观测系统定量分析方法 |
CN102023312A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统采集脚印定量分析方法 |
CN101852866A (zh) * | 2010-05-13 | 2010-10-06 | 中国石油天然气股份有限公司 | 一种叠后地震数据滤波方法 |
CN102645670A (zh) * | 2011-02-22 | 2012-08-22 | 中国石油天然气集团公司 | 一种基于叠加响应分析的观测系统优化设计方法 |
CN103558636A (zh) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种从叠后地震数据采集脚印衰减的方法 |
CN105044799A (zh) * | 2015-07-10 | 2015-11-11 | 中国石油大学(北京) | 确定三维地震观测系统面元属性均匀程度及均匀化的方法 |
Non-Patent Citations (5)
Title |
---|
"Suppression of the acquisition footprint for seismic sequence attribute mapping";Kurt J. Marfurt 等;《Geophysics》;19980630;第63卷(第3期);第1024-1035页 * |
"三维地震观测系统采集脚印定量分析";骆宗强 等;《地球物理学进展》;20120130;第27卷(第2期);第548-554页 * |
"三维观测系统采集脚印定量分析技术";碗学俭 等;《石油地球物理勘探》;20110630;第46卷(第3期);第357-363页 * |
"压制三维地震数据采集脚印的方法研究";董世泰 等;《石油地球物理勘探》;20070228;第42卷(第1期);第8页右栏第2段、第9页右栏第1段及附图4、6-7 * |
"基于地震物理模拟的采集脚印分析";熊金良 等;《石油地球物理勘探》;20061031;第41卷(第5期);第493-497页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105425294A (zh) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110178056B (zh) | 利用小波互相关进行虚拟源去噪 | |
CN106873033B (zh) | 一种基于神经网络的多属性地震信息融合裂缝预测的方法 | |
CN102176054B (zh) | 近地表综合信息处理解释方法 | |
EP3488267B1 (en) | Seismic spectral balancing | |
CN109001813B (zh) | 一种压制多次波的方法、装置及系统 | |
CN106405637B (zh) | 一种用于砂泥岩薄互层砂体厚度预测的方法及系统 | |
Takougang et al. | Characterization of small faults and fractures in a carbonate reservoir using waveform inversion, reverse time migration, and seismic attributes | |
CN102736107A (zh) | 一种能量约束非均质储层厚度识别系统 | |
Wang et al. | Separation and imaging of seismic diffractions using a localized rank-reduction method with adaptively selected ranks | |
CN109343115B (zh) | 一种基于测井约束的含气储层刻画方法 | |
CN107356965B (zh) | 基于加权叠加噪音压制策略的反射系数反演储层预测方法 | |
CN109164485A (zh) | 一种影响低序级断层识别精度的定量分析方法 | |
CN110799857A (zh) | 抑制地震数据中的噪声 | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN105474048B (zh) | 利用波束分解预测地震数据中的层间多次波 | |
Wu et al. | Incremental correlation of multiple well logs following geologically optimal neighbors | |
CN106772593A (zh) | 绕射波的成像方法及装置 | |
Martorana et al. | Assessment of a geological model by surface wave analyses | |
Alcalde et al. | 3-D reflection seismic imaging of the Hontomín structure in the Basque–Cantabrian Basin (Spain) | |
CN105425294B (zh) | 采集脚印的确定方法和装置 | |
CN108254785A (zh) | 储层确定方法和装置 | |
CN106353807A (zh) | 裂缝识别方法和装置 | |
CN104133244A (zh) | 一种震源横波勘探的近地表静校正处理方法 | |
CN109471162A (zh) | 层间多次波处理方法、系统、电子设备及可读介质 | |
CN108345032A (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 |