CN101852866A - 一种叠后地震数据滤波方法 - Google Patents
一种叠后地震数据滤波方法 Download PDFInfo
- Publication number
- CN101852866A CN101852866A CN201010176516A CN201010176516A CN101852866A CN 101852866 A CN101852866 A CN 101852866A CN 201010176516 A CN201010176516 A CN 201010176516A CN 201010176516 A CN201010176516 A CN 201010176516A CN 101852866 A CN101852866 A CN 101852866A
- Authority
- CN
- China
- Prior art keywords
- data acquisition
- seismic data
- amplitude
- acquisition point
- domain
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001914 filtration Methods 0.000 title claims abstract description 24
- 238000013016 damping Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种叠后地震数据滤波方法,包括:获取叠后地震数据;将获取的叠后地震数据由时间域转换到时间切片域;将时间切片域中每一时间切片的地震数据从空间域转换到波数域;在波数域上获取满足下述条件的地震数据采集点:1)在以原点为中心,半径为R0的圆形区域外;2)地震数据采集点的振幅为局部振幅极大值;3)地震数据采集点的振幅大于全局振幅最大值的d倍;4)地震数据采集点的振幅大于以满足条件1)、2)、3)的地震数据采集点为中心,半径为R1的区域内的振幅平均值的f倍;根据获取的地震数据采集点构造滤波器;利用该滤波器对所述的地震数据采集点进行滤波。本发明能够消除采集脚印对地震数据的干扰,提高数据的信噪比。
Description
技术领域
本发明是关于石油勘探技术,特别是关于地震数据处理技术,具体的讲是关于一种叠后地震数据滤波方法。
背景技术
在地震数据的采集过程中所使用的观测系统会对采集到的地震数据产生干扰,可以形象地将这种干扰被喻为在地震数据中留下了脚印,即采集脚印干扰。这种脚印表现为在叠后地震数据体的时间切片上出现有规律的振幅变化,是调制在地层反射信号上的一种系统性噪声,它严重的影响了地震属性反演的精确度。
现有压制采集脚印干扰的技术一般是根据观测系统,设计出压制采集脚印干扰的滤波器,在地震数据叠后压制采集脚印干扰。但是在实际的运用过程中,由于斜交等不规则观测系统的使用,以及受观测系统改变,地表起伏(如地面遇到建筑物)等因素的影响,采集脚印干扰的周期性分布会受到扰乱,导致根据观测系统设计滤波器变得非常困难。
发明内容
本发明提供一种叠后地震数据滤波方法,以有效去除有数据采集、处理过程中引入的脚印干扰。
为了实现上述目的,本发明提供了一种叠后地震数据滤波方法,所述的方法包括:
获取叠后地震数据;
将获取的叠后地震数据由时间域转换到时间切片域;
将所述时间切片域中每一时间切片的地震数据从空间域转换到波数域;在所述波数域上获取满足下述条件的地震数据采集点:
1)在以原点为中心,半径为R0的圆形区域外;
2)所述地震数据采集点的振幅为局部振幅极大值;
3)所述地震数据采集点的振幅大于全局振幅最大值的d倍;
4)所述地震数据采集点的振幅大于以满足条件1)、2)、3)的地震数据采集点为中心,半径为R1的区域内的振幅平均值的f倍;
根据获取的地震数据采集点构造滤波器;
利用所述的滤波器对所述的地震数据采集点进行滤波。
本发明的有益技术效果:
本发明可以有效的去除由采集、处理过程中引入的采集脚印,提高了信噪比,消除了采集脚印对地震数据属性的影响,有利于开展地震资料的综合研究,提高了储层预测精度,能够更加准确地刻画断裂系统,精确预测砂体的分布范围和展布特征。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例一正演模型的叠后地震数据示意图;
图2为本发明实施例一切除第700ms处的时间切片的叠后地震数据示意图;
图3为本发明实施例一变换得到的二维波数谱示意图;
图4为本发明实施例一二维波数域滤波器示意图;
图5为本发明实施例一由二维波数域(Kx,Ky)上被滤波了的数据变换回时间切片上的二维空间域的结果示意图;
图6为本发明实施例的方法流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
实施例一
输入正演模型的叠后地震数据,总道数260×300,共中心点(Common MidPoint CMP)线数260,每条CMP线有CMP点数300,正演模型的叠后地震数据如图1所示,图1显示了距离地面不同深度的不同叠后地震数据。
将上述正演模型的叠后地震数据由时间域转换到时间切片域。根据上述叠后地震数据,可以获得时间点的时间切片,如图2所示,为第700ms处切除的时间切片。
利用二维傅立叶变换把该时间切片的地震数据由空间域(x,y)转换到波数域(Kx,Ky),图3为由空间域(x,y)转换到波数域(Kx,Ky)变换得到的二维波数谱示意图。
在二维波数谱中,按如下条件查找被采集脚印噪音干扰的点(地震数据采集点),用于构造二维波数域滤波器。
①在以原点(0,0)点为中心半径为R0(单位:周期数/道)的圆形保护区域外。其中,半径R0的取值范围为0.03~0.08(周期数/道),该取值范围仅为本发明的实施例,并非用于限定本发明。
②地震数据采集点的振幅为局部振幅极大值的点。
③地震数据采集点的振幅大于全局振幅最大值的d倍;其中,d的取值范围为0.0007~0.005,该取值范围仅为本发明的实施例,并非用于限定本发明。
④地震数据采集点的振幅大于以满足条件①、②、③的地震数据采集点为中心,半径为R1(单位:周期数/道)的区域内的振幅平均值的f倍;其中,半径R1的取值范围为0.02~0.07(周期数/道);f的取值范围为1.0~1.3,该取值范围仅为本发明的实施例,并非用于限定本发明。
根据满足①、②、③、④的这些地震数据采集点构造一个二维波数域滤波器,将以这些地震数据采集点为中心,半径为r(单位:周期数/道)的压制区域内振幅衰减系数设定为e。其中,半径r的取值范围为0.01~0.04(周期数/道);e的取值范围为0.01~0.05。构造的二维波数域滤波器如图4所示。
在二维波数域(Kx,Ky)中,应用上述构造的二维波数域滤波器,对波数域(Kx,Ky)的叠后地震数据进行滤波,压制采集脚印噪音干扰。
利用二维傅立叶反变换,把在二维波数域(Kx,Ky)上被滤波了的叠后地震数据变换回时间切片上的二维空间域(x,y),得到二维空间域的叠后地震数据如图5所示。
将图5中滤波后的叠后地震数据与图2中滤波前的叠后地震数据进行比较,可以发现利用本申请的叠后地震数据滤波方法消除了脚印干扰的图5明显比图2明显清晰。
完成叠后地震数据地滤波后,把地震数据由时间切片域(t,x,y)转置回时间(x,y,t)域,可以根据从时间切片域转换到时间域的地震数据生成地层剖面图像,得到压制了采集脚印干扰的地层剖面图像,便于对叠后地震数据进行准确的分析。
实施例二
上述实施例以正演模型的叠后地震数据对本发明的叠后地震数据的滤波方法进行了说明,本实施例以实际的叠后地震数据为例说明叠后地震数据的滤波方法。
输入实际的叠后地震数据,总道数420×600,其中CMP线数420,每条CMP线有CMP点数600。
将上述正演模型的叠后地震数据由时间域转换到时间切片域。根据上述叠后地震数据,可以获得时间点的时间切片,例如切除第600ms处的时间切片。
利用二维傅立叶变换把第600ms处的时间切片的地震数据由空间域(x,y)转换到波数域(Kx,Ky)。
在二维波数谱中,按如下条件查找被采集脚印噪音干扰的点(地震数据采集点),用于构造二维波数域滤波器。
①在以原点(0,0)点为中心半径为R0(单位:周期数/道)的圆形保护区域外。其中,半径R0的取值范围为0.03~0.08(周期数/道),该取值范围仅为本发明的实施例,并非用于限定本发明。
②地震数据采集点的振幅为局部振幅极大值的点。
③地震数据采集点的振幅大于全局振幅最大值的d倍;其中,d的取值范围为0.0007~0.005,该取值范围仅为本发明的实施例,并非用于限定本发明。
④地震数据采集点的振幅大于以满足条件①、②、③的地震数据采集点为中心,半径为R1(单位:周期数/道)的区域内的振幅平均值的f倍;其中,半径R1的取值范围为0.02~0.07(周期数/道);f的取值范围为1.0~1.3,该取值范围仅为本发明的实施例,并非用于限定本发明。
根据满足①、②、③、④的这些地震数据采集点构造一个二维波数域滤波器,将以这些地震数据采集点为中心,半径为r(单位:周期数/道)的压制区域内振幅衰减系数设定为e。其中,半径r的取值范围为0.01~0.04(周期数/道);e的取值范围为0.01~0.05。
在二维波数域(Kx,Ky)中,应用上述构造的二维波数域滤波器,对波数域(Kx,Ky)的叠后地震数据进行滤波,压制采集脚印噪音干扰。
利用二维傅立叶反变换,把在二维波数域(Kx,Ky)上被滤波了的叠后地震数据变换回时间切片上的二维空间域(x,y)。
滤波后的叠后地震数据与滤波前的叠后地震数据相比,已经滤除了采集脚印的干扰。
完成叠后地震数据地滤波后,把地震数据由时间切片域(t,x,y)转置回时间(x,y,t)域,可以根据从时间切片域转换到时间域的地震数据生成地层剖面图像,得到压制了采集脚印干扰的地层剖面图像,便于对叠后地震数据进行准确的分析。
图6为本发明实施例的方法流程图,如图6所示,所述的方法包括:
步骤S601:获取叠后地震数据;
步骤S602:将获取的叠后地震数据由时间域转换到时间切片域;
步骤S603:将所述时间切片域中每一时间切片的地震数据从空间域转换到波数域;
步骤S604:在所述波数域上获取满足下述条件的地震数据采集点:
1)在以原点为中心,半径为R0的圆形区域外。半径R0的取值范围为0.03~0.08(周期数/道)。
2)所述地震数据采集点的振幅为局部振幅极大值;
3)所述地震数据采集点的振幅大于全局振幅最大值的d倍,d的取值范围为0.0007~0.005。
4)所述地震数据采集点的振幅大于以满足条件1)、2)、3)的地震数据采集点为中心,半径为R1的区域内的振幅平均值的f倍,半径R1的取值范围为0.02~0.07(周期数/道);f的取值范围为1.0~1.3。
步骤S605:根据获取的地震数据采集点构造滤波器;将以这些地震数据采集点为中心,半径为r(单位:周期数/道)的压制区域内振幅衰减系数设定为e。其中,半径r的取值范围为0.01~0.04(周期数/道);e的取值范围为0.01~0.05。
步骤S606:利用所述的滤波器对所述的地震数据采集点进行滤波。
本发明的有益技术效果:本发明可以有效的去除由采集、处理过程中引入的采集脚印,提高了信噪比,消除了采集脚印对地震数据属性的影响,有利于开展地震资料的综合研究,提高了储层预测精度,能够更加准确地刻画断裂系统,精确预测砂体的分布范围和展布特征。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种叠后地震数据滤波方法,其特征在于,所述的方法包括:
获取叠后地震数据;
将获取的叠后地震数据由时间域转换到时间切片域;
将所述时间切片域中每一时间切片的地震数据从空间域转换到波数域;在所述波数域上获取满足下述条件的地震数据采集点:
1)在以原点为中心,半径为R0的圆形区域外;
2)所述地震数据采集点的振幅为局部振幅极大值;
3)所述地震数据采集点的振幅大于全局振幅最大值的d倍;
4)所述地震数据采集点的振幅大于以满足条件1)、2)、3)的地震数据采集点为中心,半径为R1的区域内的振幅平均值的f倍;
根据获取的地震数据采集点构造滤波器;
利用所述的滤波器对所述的地震数据采集点进行滤波。
2.如权利要求1所述的方法,其特征在于,在利用所述的滤波器对所述的地震数据采集点进行滤波之后,所述的方法还包括:
将滤波器滤波后的地震数据从波数域转换到空间域;
将滤波器滤波后的地震数据从时间切片域转换到时间域;
根据从时间切片域转换到时间域的地震数据生成地层剖面图像。
3.如权利要求1所述的方法,其特征在于,根据获取的地震数据采集点构造滤波器,包括:
将以获取的地震数据采集点为中心,半径为r的压制区域内的振幅衰减系数设定为e。
4.如权利要求1所述的方法,其特征在于,半径R0的取值范围为0.03~0.08(周期数/道)。
5.如权利要求1所述的方法,其特征在于,d的取值范围为0.0007~0.005。
6.如权利要求1所述的方法,其特征在于,半径R1的取值范围为0.02~0.07(周期数/道)。
7.如权利要求1所述的方法,其特征在于,f的取值范围为1.0~1.3。
8.如权利要求3所述的方法,其特征在于,半径r的取值范围为0.01~0.04(周期数/道)。
9.如权利要求3所述的方法,其特征在于,e的取值范围为0.01~0.05。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101765163A CN101852866B (zh) | 2010-05-13 | 2010-05-13 | 一种叠后地震数据滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101765163A CN101852866B (zh) | 2010-05-13 | 2010-05-13 | 一种叠后地震数据滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101852866A true CN101852866A (zh) | 2010-10-06 |
CN101852866B CN101852866B (zh) | 2012-01-04 |
Family
ID=42804438
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101765163A Active CN101852866B (zh) | 2010-05-13 | 2010-05-13 | 一种叠后地震数据滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101852866B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558636A (zh) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种从叠后地震数据采集脚印衰减的方法 |
CN104570120A (zh) * | 2014-12-19 | 2015-04-29 | 中国石油天然气集团公司 | 低信噪比地区微测井资料减噪方法及系统 |
CN104865603A (zh) * | 2014-02-26 | 2015-08-26 | 中国石油化工股份有限公司 | 一种针对倾斜层的svd滤波方法及装置 |
CN105425294A (zh) * | 2015-11-30 | 2016-03-23 | 中国石油天然气集团公司 | 采集脚印的确定方法和装置 |
CN105607125A (zh) * | 2016-01-15 | 2016-05-25 | 吉林大学 | 基于块匹配算法和奇异值分解的地震资料噪声压制方法 |
CN106353816A (zh) * | 2016-08-09 | 2017-01-25 | 中国石油天然气集团公司 | 一种地震采集脚印噪音压制方法及系统 |
CN106908836A (zh) * | 2015-12-23 | 2017-06-30 | 中国石油天然气股份有限公司 | 采集脚印压制方法及系统 |
CN112255684A (zh) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | 一种地震数据的采集脚印的处理方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1058131B1 (en) * | 1999-06-03 | 2004-09-01 | Jason Geosystems B.V. | Method of determining spatial changes in subsurface structure stratigraphy, lithology and fluid content and of reducing seismic noise |
CN101419292B (zh) * | 2007-10-25 | 2011-05-25 | 中国石油天然气集团公司 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
US20090192717A1 (en) * | 2008-01-24 | 2009-07-30 | Terje Iversen | Methods and apparatuses for dynamic filtering of geometric primitives in 3d space |
-
2010
- 2010-05-13 CN CN2010101765163A patent/CN101852866B/zh active Active
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558636B (zh) * | 2013-11-08 | 2016-05-25 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种从叠后地震数据采集脚印衰减的方法 |
CN103558636A (zh) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种从叠后地震数据采集脚印衰减的方法 |
CN104865603A (zh) * | 2014-02-26 | 2015-08-26 | 中国石油化工股份有限公司 | 一种针对倾斜层的svd滤波方法及装置 |
CN104570120A (zh) * | 2014-12-19 | 2015-04-29 | 中国石油天然气集团公司 | 低信噪比地区微测井资料减噪方法及系统 |
CN105425294B (zh) * | 2015-11-30 | 2018-01-05 | 中国石油天然气集团公司 | 采集脚印的确定方法和装置 |
CN105425294A (zh) * | 2015-11-30 | 2016-03-23 | 中国石油天然气集团公司 | 采集脚印的确定方法和装置 |
CN106908836B (zh) * | 2015-12-23 | 2019-03-15 | 中国石油天然气股份有限公司 | 采集脚印压制方法及系统 |
CN106908836A (zh) * | 2015-12-23 | 2017-06-30 | 中国石油天然气股份有限公司 | 采集脚印压制方法及系统 |
CN105607125B (zh) * | 2016-01-15 | 2018-04-13 | 吉林大学 | 基于块匹配算法和奇异值分解的地震资料噪声压制方法 |
CN105607125A (zh) * | 2016-01-15 | 2016-05-25 | 吉林大学 | 基于块匹配算法和奇异值分解的地震资料噪声压制方法 |
CN106353816A (zh) * | 2016-08-09 | 2017-01-25 | 中国石油天然气集团公司 | 一种地震采集脚印噪音压制方法及系统 |
CN112255684A (zh) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | 一种地震数据的采集脚印的处理方法及装置 |
CN112255684B (zh) * | 2020-10-26 | 2024-03-01 | 中国石油天然气集团有限公司 | 一种地震数据的采集脚印的处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101852866B (zh) | 2012-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101852866B (zh) | 一种叠后地震数据滤波方法 | |
Mousavi et al. | Automatic noise-removal/signal-removal based on general cross-validation thresholding in synchrosqueezed domain and its application on earthquake data | |
US11880011B2 (en) | Surface wave prediction and removal from seismic data | |
Roten et al. | Two-dimensional resonances in Alpine valleys identified from ambient vibration wavefields | |
CN104007469A (zh) | 一种基于曲波变换的弱地震信号重构方法 | |
Kirkpatrick et al. | Subduction megathrust heterogeneity characterized from 3D seismic data | |
Pang et al. | Automatic passive data selection in time domain for imaging near-surface surface waves | |
CN106896409B (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
CN102043165B (zh) | 基于基追踪算法的面波分离与压制方法 | |
CN108387933A (zh) | 一种确定地层品质因子的方法、装置及系统 | |
CN109669212A (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
Wu et al. | Microseismic source locations with deconvolution migration | |
CN107179551B (zh) | 一种利用微震记录对地下构造直接成像的方法 | |
Picozzi et al. | Deblurring of frequency–wavenumber images from small-scale seismic arrays | |
CN104345341A (zh) | 一种基于区域约束的分频段能量地震面波处理方法 | |
CN104597502A (zh) | 一种新的石油地震勘探数据去噪方法 | |
Atterholt et al. | A unified wavefield-partitioning approach for distributed acoustic sensing | |
Wang et al. | Ground roll attenuation using polarization analysis in the tfk domain | |
Cheng et al. | Phase-weighted slant stacking for surface wave dispersion measurement | |
CN105445801A (zh) | 一种消除二维地震资料随机噪音的处理方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
Martorana et al. | Assessment of a geological model by surface wave analyses | |
CN102841380B (zh) | 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 | |
CN106950600B (zh) | 一种近地表散射面波的去除方法 | |
Bai et al. | Least-squares Gaussian beam transform for seismic noise attenuation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |