CN101852866B - 一种叠后地震数据滤波方法 - Google Patents

一种叠后地震数据滤波方法 Download PDF

Info

Publication number
CN101852866B
CN101852866B CN2010101765163A CN201010176516A CN101852866B CN 101852866 B CN101852866 B CN 101852866B CN 2010101765163 A CN2010101765163 A CN 2010101765163A CN 201010176516 A CN201010176516 A CN 201010176516A CN 101852866 B CN101852866 B CN 101852866B
Authority
CN
China
Prior art keywords
seismic data
amplitude
post
domain
radius
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
Application number
CN2010101765163A
Other languages
English (en)
Other versions
CN101852866A (zh
Inventor
王彦仓
叶秋焱
张树森
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN2010101765163A priority Critical patent/CN101852866B/zh
Publication of CN101852866A publication Critical patent/CN101852866A/zh
Application granted granted Critical
Publication of CN101852866B publication Critical patent/CN101852866B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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 (7)

1.一种叠后地震数据滤波方法，其特征在于，所述的方法包括：
获取叠后地震数据；
将获取的叠后地震数据由时间域转换到时间切片域；
将所述时间切片域中每一时间切片的地震数据从空间域转换到波数域；在所述波数域上获取满足下述条件的地震数据采集点：
1)在以原点为中心，半径为R0的圆形区域外；
2)所述地震数据采集点的振幅为局部振幅极大值；
3)所述地震数据采集点的振幅大于全局振幅最大值的d倍；
4)所述地震数据采集点的振幅大于以满足条件1)、2)、3)的地震数据采集点为中心，半径为R1的区域内的振幅平均值的f倍；
将以获取的地震数据采集点为中心，半径为r的压制区域内的振幅衰减系数设定为e，构造滤波器；
利用所述的滤波器对所述的地震数据采集点进行滤波；
将滤波器滤波后的地震数据从波数域转换到空间域；
将滤波器滤波后的地震数据从时间切片域转换到时间域；
根据从时间切片域转换到时间域的地震数据生成地层剖面图像。
2.如权利要求1所述的方法，其特征在于，半径R0的取值范围为0.03～0.08周期数/道。
3.如权利要求1所述的方法，其特征在于，d的取值范围为0.0007～0.005。
4.如权利要求1所述的方法，其特征在于，半径R1的取值范围为0.02～0.07周期数/道。
5.如权利要求1所述的方法，其特征在于，f的取值范围为1.0～1.3。
6.如权利要求3所述的方法，其特征在于，半径r的取值范围为0.01～0.04周期数/道。 
7.如权利要求3所述的方法，其特征在于，e的取值范围为0.01～0.05。 
CN2010101765163A 2010-05-13 2010-05-13 一种叠后地震数据滤波方法 Active CN101852866B (zh)

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 CN101852866A (zh) 2010-10-06
CN101852866B true 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)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103558636B (zh) * 2013-11-08 2016-05-25 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种从叠后地震数据采集脚印衰减的方法
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 中国石油天然气集团公司 采集脚印的确定方法和装置
CN106908836B (zh) * 2015-12-23 2019-03-15 中国石油天然气股份有限公司 采集脚印压制方法及系统
CN105607125B (zh) * 2016-01-15 2018-04-13 吉林大学 基于块匹配算法和奇异值分解的地震资料噪声压制方法
CN106353816B (zh) * 2016-08-09 2018-10-16 中国石油天然气集团公司 一种地震采集脚印噪音压制方法及系统
CN112255684B (zh) * 2020-10-26 2024-03-01 中国石油天然气集团有限公司 一种地震数据的采集脚印的处理方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1058131A1 (en) * 1999-06-03 2000-12-06 Jason Geosystems B.V. Method and determining spatial changes in subsurface structure stratigraphy, lithology and fluid content and of reducing seismic noise
CN101419292A (zh) * 2007-10-25 2009-04-29 中国石油天然气集团公司 采用纵波源多分量地震数据生成横波地震剖面的方法
WO2009094067A2 (en) * 2008-01-24 2009-07-30 Schlumberger Canada Limited Methods and apparatuses for dynamic filtering of geometric primitives in 3d space

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1058131A1 (en) * 1999-06-03 2000-12-06 Jason Geosystems B.V. Method and determining spatial changes in subsurface structure stratigraphy, lithology and fluid content and of reducing seismic noise
CN101419292A (zh) * 2007-10-25 2009-04-29 中国石油天然气集团公司 采用纵波源多分量地震数据生成横波地震剖面的方法
WO2009094067A2 (en) * 2008-01-24 2009-07-30 Schlumberger Canada Limited Methods and apparatuses for dynamic filtering of geometric primitives in 3d space

Also Published As

Publication number Publication date
CN101852866A (zh) 2010-10-06

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
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) 一种基于波动方程边值反演的变深度缆鬼波压制方法
Williams et al. Scholte wave inversion and passive source imaging with ocean-bottom DAS
Wu et al. Microseismic source locations with deconvolution migration
CN109669212A (zh) 地震数据处理方法、地层品质因子估算方法与装置
CN102043165A (zh) 基于基追踪算法的面波分离与压制方法
Wang et al. Ground roll attenuation using polarization analysis in the tfk domain
Chambers Using DAS to investigate traffic patterns at Brady Hot Springs, Nevada, USA
CN104597502A (zh) 一种新的石油地震勘探数据去噪方法
Atterholt et al. A unified wavefield-partitioning approach for distributed acoustic sensing
Cheng et al. Phase-weighted slant stacking for surface wave dispersion measurement
Ikeda et al. Surface wave attenuation in the shallow subsurface from multichannel–multishot seismic data: a new approach for detecting fractures and lithological discontinuities
CN102841380B (zh) 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法
Zhang et al. Improving receiver function imaging with high-resolution Radon transform
Martorana et al. Assessment of a geological model by surface wave analyses
Taipodia et al. Influence of data acquisition and signal preprocessing parameters on the resolution of dispersion image from active MASW survey
Hu et al. Data-driven dispersive surface-wave prediction and mode separation using high-resolution dispersion estimation
Zhang et al. Separating the scattered wavefield from teleseismic P using curvelets on the long beach array data set
Ugalde et al. Temporal variations of seismic velocity at Paradox Valley, Colorado, using passive image interferometry
Ebadi Coherent and incoherent seismic noise attenuation using parabolic radon transform and its application in environmental geophysics

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