CN111624659A - 一种地震数据的时变带通滤波方法及装置 - Google Patents
一种地震数据的时变带通滤波方法及装置 Download PDFInfo
- Publication number
- CN111624659A CN111624659A CN202010502784.3A CN202010502784A CN111624659A CN 111624659 A CN111624659 A CN 111624659A CN 202010502784 A CN202010502784 A CN 202010502784A CN 111624659 A CN111624659 A CN 111624659A
- Authority
- CN
- China
- Prior art keywords
- time
- data
- leveling
- seismic data
- horizon
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000011218 segmentation Effects 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/21—Frequency-domain filtering, e.g. band pass
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
本发明公开了一种地震数据的时变带通滤波方法及装置,所述方法包括以下步骤:S1.输入多道地震数据;时变带通滤波是对地震数据的中深层开时窗进行的,在时窗上边界附近选取一套连续的层位,对于每道地震数据沿层拾取对应的层位时间;S2.确定步骤S1得到层位时间中的最大值,然后确定层拉平时间,计算层拉平时移量,对输入的地震数据沿层时移拉平;S3.对时移拉平后的数据做多道数据分段时变带通滤波;S4.利用层位时间、拉平时间,计算层反拉平时移量,对时变带通滤波后的数据做层反拉平。本发明采用沿层拾取、层拉平、多道数据分段时变带通滤波、层反拉平的方式,能快速处理多道数据、可变分段的时变带通滤波,便于地震资料后续处理。
Description
技术领域
本发明涉及地球物理勘探的地震数据处理方法,属于地震处理技术领域,特别是涉及一种地震数据的时变带通滤波方法及装置。
背景技术
地震数据反射波的主频随着旅行时增加变低、频宽随着旅行时增加变窄,即频率是时变的。因此,在地震数据处理中经常使用时变带通滤波。通常,时变带通滤波是分段进行的,在段与段之间的过渡带会出现明显的接缝;多数据、多段时变带通不易实现。
发明内容
本发明的目的在于克服现有技术的不足,提供一种地震数据的时变带通滤波方法及装置,根据地震数据频率分布的特点,采用沿层拾取、层拉平、多道数据分段时变带通滤波、层反拉平的方式,实现了多数据、可变分段的时变带通滤波。
本发明的目的是通过以下技术方案来实现的:一种地震数据的时变带通滤波方法,包括以下步骤:
其中,Ti layer是第i道数据的层位时间,i=1,2,3,...,N;
S2.确定步骤S1得到层位时间中的最大值,然后确定层拉平时间Tlp,计算层拉平时移量,对输入的地震数据沿层时移拉平;
进一步地,所述步骤S2包括以下子步骤:
S202.确定拉平时间Tlp,Tlp≥Tmax;
S205.计算每一道数据的层拉平时移量:
其中,Ti Shift是第i道数据的层拉平时移量;
S206.对于每一道地震数据,利用步骤S205得到的时移量时移步骤S204得到的地震数据。
进一步地,所述步骤S3包括以下子步骤:
S301.对时移拉平后的每一道地震数据做傅里叶变换,计算对应的振幅谱A、相位谱φ;
S302.确定时移拉平后各道数据分段时变带通滤波振幅谱At为:
S303.对于每一道数据,利用振幅谱At和对应的相位谱φ做反傅里叶变换。
进一步地,所述步骤S4包括以下子步骤:
S401.对于每一道时变带通滤波后的数据,计算对应的层反拉平时移量:
{Ti InvShift=Ti layer-Tlp};i=1,2,3,...,N;
其中,Ti InvShift是第i道数据的层反拉平时移量;
S402.对于每一道时变带通滤波后的数据,利用层反拉平时移量时移时变带通滤波后的地震数据。
一种地震数据的时变带通滤波装置,包括:
数据输入单元,用于输入多道地震数据;
层位时间拾取单元,用于在时窗上边界附近选取一套连续的层位,对于每道地震数据沿层拾取对应的层位时间,最终得到N道输入地震数据的层位时间集合;
时移拉平单元,用于确定层位时间中的最大值和层拉平时间,并计算层拉平时移量,对输入的地震数据沿层时移拉平;
时变带通滤波单元,用于对时移拉平后的数据做多道数据分段时变带通滤波;
层反拉平单元,用于根据层位时间和拉平时间,计算层反拉平时移量,对时变带通滤波后的数据做层反拉平。
本发明的有益效果是:本发明采用沿层拾取、层拉平、多道数据分段时变带通滤波、层反拉平的方式,能快速处理多道数据、可变分段的时变带通滤波,便于地震资料后续处理。
附图说明
图1为本发明的方法流程图;
图2为实施例中输入地震数据示意图;
图3为实施例中沿层拾取时间层位示意图;
图4为实施例中层拉平地震数据示意图;
图5为实施例中时变带通滤波的层拉平地震数据示意图;
图6为实施例中层反拉平地震数据示意图;
图7为本发明的装置原理示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,一种地震数据的时变带通滤波方法,包括以下步骤:
其中,Ti layer是第i道数据的层位时间,i=1,2,3,...,N;
在本申请的实施例中,输入的地震数据如图2所示,其中横坐标为数据号;纵坐标为时间(单位:ms);进过步骤S1处理后,沿层拾取的层位时间,如图3所示,图3中,横坐标为数据号;纵坐标为时间(单位:ms)。
S2.确定步骤S1得到层位时间中的最大值,然后确定层拉平时间Tlp,计算层拉平时移量,对输入的地震数据沿层时移拉平:
S202.确定拉平时间Tlp,Tlp≥Tmax;
S205.计算每一道数据的层拉平时移量:
{Ti Shift=Tlp-Ti layer},i=1,2,3,...,N;
其中,Ti Shift是第i道数据的层拉平时移量;
S206.对于每一道地震数据,利用步骤S205得到的时移量时移步骤S204得到的地震数据。
在本申请的实施例中,经过步骤S2处理,得到层拉平地震数据,如图4所示,横坐标为数据号;纵坐标为时间(单位:ms);
S301.对时移拉平后的每一道地震数据做傅里叶变换,计算对应的振幅谱A、相位谱φ;
S302.确定时移拉平后各道数据分段时变带通滤波振幅谱At为:
S303.对于每一道数据,利用振幅谱At和对应的相位谱φ做反傅里叶变换。
在本申请的实施例中,通过时变带通滤波得到的层拉平数据如图5所示,该图中横坐标为数据号;纵坐标为时间(单位:ms)。
S401.对于每一道时变带通滤波后的数据,计算对应的层反拉平时移量:
{Ti InvShift=Ti layer-Tlp};i=1,2,3,...,N;
其中,Ti InvShift是第i道数据的层反拉平时移量;
S402.对于每一道时变带通滤波后的数据,利用层反拉平时移量时移时变带通滤波后的地震数据。
在本申请的实施例中,层反拉平得到的地震数据如图6所示,图中,横坐标为数据号;纵坐标为时间(单位:ms)。
如图7所示,一种地震数据的时变带通滤波装置,包括:
数据输入单元,用于输入多道地震数据;
层位时间拾取单元,用于在时窗上边界附近选取一套连续的层位,对于每道地震数据沿层拾取对应的层位时间,最终得到N道输入地震数据的层位时间集合;
时移拉平单元,用于确定层位时间中的最大值和层拉平时间,并计算层拉平时移量,对输入的地震数据沿层时移拉平;
时变带通滤波单元,用于对时移拉平后的数据做多道数据分段时变带通滤波;
层反拉平单元,用于根据层位时间和拉平时间,计算层反拉平时移量,对时变带通滤波后的数据做层反拉平。
综上,本发明根据地震数据频率分布的特点,采用沿层拾取、层拉平、多道数据分段时变带通滤波、层反拉平的方式,实现了多数据、可变分段的时变带通滤波,便于地震资料后续处理。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其他不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (5)
4.根据权利要求1所述的一种地震数据的时变带通滤波方法,其特征在于:所述步骤S4包括以下子步骤:
S401.对于每一道时变带通滤波后的数据,计算对应的层反拉平时移量:
{Ti InvShift=Ti layer-Tlp};i=1,2,3,...,N;
其中,Ti InvShift是第i道数据的层反拉平时移量;
S402.对于每一道时变带通滤波后的数据,利用层反拉平时移量时移时变带通滤波后的地震数据。
5.一种地震数据的时变带通滤波装置,采用如权利要求1~4中任意一项所述的方法,其特征在于:包括:
数据输入单元,用于输入多道地震数据;
层位时间拾取单元,用于在时窗上边界附近选取一套连续的层位,对于每道地震数据沿层拾取对应的层位时间,最终得到N道输入地震数据的层位时间集合;
时移拉平单元,用于确定层位时间中的最大值和层拉平时间,并计算层拉平时移量,对输入的地震数据沿层时移拉平;
时变带通滤波单元,用于对时移拉平后的数据做多道数据分段时变带通滤波;
层反拉平单元,用于根据层位时间和拉平时间,计算层反拉平时移量,对时变带通滤波后的数据做层反拉平。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010502784.3A CN111624659B (zh) | 2020-06-05 | 2020-06-05 | 一种地震数据的时变带通滤波方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010502784.3A CN111624659B (zh) | 2020-06-05 | 2020-06-05 | 一种地震数据的时变带通滤波方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111624659A true CN111624659A (zh) | 2020-09-04 |
CN111624659B CN111624659B (zh) | 2022-07-01 |
Family
ID=72258211
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010502784.3A Active CN111624659B (zh) | 2020-06-05 | 2020-06-05 | 一种地震数据的时变带通滤波方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111624659B (zh) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3539982A (en) * | 1969-03-03 | 1970-11-10 | Texas Instruments Inc | Method and apparatus for consistent static correction of common depth point seismic signals |
EP1605279A2 (de) * | 2004-06-09 | 2005-12-14 | Ernst D. Rode | Verfahren zur Erfassung und Erkundung von unterirdischen Kohlenwasserstoffvorkommen durch Empfänger von akustischen Wellen in einem Frequenzbereich von 0,2 bis 30 Hz |
CN102109612A (zh) * | 2009-12-23 | 2011-06-29 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
US20110218737A1 (en) * | 2010-03-05 | 2011-09-08 | Vialogy Llc | Active Noise Injection Computations for Improved Predictability in Oil and Gas Reservoir Discovery and Characterization |
CN102778693A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 |
CN104267431A (zh) * | 2014-10-13 | 2015-01-07 | 成都北方石油勘探开发技术有限公司 | 一种地震构造解释的小断层获取方法 |
CN104280765A (zh) * | 2013-07-11 | 2015-01-14 | 中国石油化工股份有限公司 | 基于变子波反射系数反演的地震高分辨处理方法 |
US9105075B1 (en) * | 2013-02-06 | 2015-08-11 | Ihs Global Inc. | Enhancing seismic features using an optical filter array |
CN105093324A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种剩余静校正方法 |
CN105629300A (zh) * | 2015-12-30 | 2016-06-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高复杂构造偏移数据信噪比的方法 |
CN106709566A (zh) * | 2016-12-07 | 2017-05-24 | 广东工业大学 | 一种基于深度学习的数据缺失值回填方法 |
CN106780347A (zh) * | 2017-02-09 | 2017-05-31 | 浙江科技学院 | 一种基于oct图像处理的枇杷早期瘀伤鉴别方法 |
CN107607992A (zh) * | 2017-08-24 | 2018-01-19 | 电子科技大学 | 基于卷积神经网络的多波匹配方法 |
CN107942387A (zh) * | 2017-10-23 | 2018-04-20 | 中国石油天然气集团公司 | 一种续至波衰减处理方法及装置 |
CN108571000A (zh) * | 2018-05-15 | 2018-09-25 | 中铁第六勘察设计院集团有限公司 | 一种沿河地段道路下穿既有铁路时地下水路堑的施工方法 |
CN109829398A (zh) * | 2019-01-16 | 2019-05-31 | 北京航空航天大学 | 一种基于三维卷积网络的视频中的目标检测方法 |
CN110399798A (zh) * | 2019-06-25 | 2019-11-01 | 朱跃飞 | 一种基于深度学习的离散图片文件信息提取系统及方法 |
-
2020
- 2020-06-05 CN CN202010502784.3A patent/CN111624659B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3539982A (en) * | 1969-03-03 | 1970-11-10 | Texas Instruments Inc | Method and apparatus for consistent static correction of common depth point seismic signals |
EP1605279A2 (de) * | 2004-06-09 | 2005-12-14 | Ernst D. Rode | Verfahren zur Erfassung und Erkundung von unterirdischen Kohlenwasserstoffvorkommen durch Empfänger von akustischen Wellen in einem Frequenzbereich von 0,2 bis 30 Hz |
CN102109612A (zh) * | 2009-12-23 | 2011-06-29 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
US20110218737A1 (en) * | 2010-03-05 | 2011-09-08 | Vialogy Llc | Active Noise Injection Computations for Improved Predictability in Oil and Gas Reservoir Discovery and Characterization |
CN102778693A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 |
US9105075B1 (en) * | 2013-02-06 | 2015-08-11 | Ihs Global Inc. | Enhancing seismic features using an optical filter array |
CN104280765A (zh) * | 2013-07-11 | 2015-01-14 | 中国石油化工股份有限公司 | 基于变子波反射系数反演的地震高分辨处理方法 |
CN105093324A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种剩余静校正方法 |
CN104267431A (zh) * | 2014-10-13 | 2015-01-07 | 成都北方石油勘探开发技术有限公司 | 一种地震构造解释的小断层获取方法 |
CN105629300A (zh) * | 2015-12-30 | 2016-06-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高复杂构造偏移数据信噪比的方法 |
CN106709566A (zh) * | 2016-12-07 | 2017-05-24 | 广东工业大学 | 一种基于深度学习的数据缺失值回填方法 |
CN106780347A (zh) * | 2017-02-09 | 2017-05-31 | 浙江科技学院 | 一种基于oct图像处理的枇杷早期瘀伤鉴别方法 |
CN107607992A (zh) * | 2017-08-24 | 2018-01-19 | 电子科技大学 | 基于卷积神经网络的多波匹配方法 |
CN107942387A (zh) * | 2017-10-23 | 2018-04-20 | 中国石油天然气集团公司 | 一种续至波衰减处理方法及装置 |
CN108571000A (zh) * | 2018-05-15 | 2018-09-25 | 中铁第六勘察设计院集团有限公司 | 一种沿河地段道路下穿既有铁路时地下水路堑的施工方法 |
CN109829398A (zh) * | 2019-01-16 | 2019-05-31 | 北京航空航天大学 | 一种基于三维卷积网络的视频中的目标检测方法 |
CN110399798A (zh) * | 2019-06-25 | 2019-11-01 | 朱跃飞 | 一种基于深度学习的离散图片文件信息提取系统及方法 |
Non-Patent Citations (1)
Title |
---|
张孝珍: "地震勘探中的去噪技术新进展", 《勘探地球物理进展》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111624659B (zh) | 2022-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109765615A (zh) | 一种地层品质因子反演方法及装置 | |
Kapetanidis et al. | Estimation of arrival-times in intense seismic sequences using a Master-Events methodology based on waveform similarity | |
CN108427140B (zh) | 一种用于小尺度缝洞储集体地震识别的方法 | |
CN105093318B (zh) | 一种自适应波动方程波场延拓静校正方法 | |
CN110568490B (zh) | 一种高速层顶薄储层的识别方法 | |
CN111624659B (zh) | 一种地震数据的时变带通滤波方法及装置 | |
CN104656140A (zh) | 一种基于垂直地震衰减规律的中值滤波方法以及系统 | |
CN110618457B (zh) | 一种可控震源地震数据的联合去噪方法及其装置 | |
CN105093312A (zh) | 基于频率域多阶微分的地震相对波阻抗预测方法和装置 | |
CN102012520B (zh) | 一种海上单点采集地震数据室内动态组合的方法 | |
CN109884701B (zh) | 地质体散射角度导向深度成像方法 | |
CN110471105B (zh) | 一种层位约束下基于地震解析道的相干增强断层识别方法 | |
CN106772596A (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN112213776A (zh) | 叠前道集和vsp资料联合层控q模型建立方法 | |
Clark et al. | Eocene turbidite-population statistics from shelf edge to basin floor, Spitsbergen, Svalbard | |
CN112859172A (zh) | 纵横波微测井资料处理方法及装置 | |
CN111060966B (zh) | 一种用于静校正技术的初至拾取方法 | |
CN110673211A (zh) | 一种基于测井与地震数据的品质因子建模方法 | |
CN112241025A (zh) | 一种井震联合地层压力确定方法及系统 | |
CN111694055A (zh) | 一种强反射屏蔽下avo属性提取和分析的方法和装置 | |
CN111060970B (zh) | 一种储层约束反褶积的方法及应用 | |
CN113589384B (zh) | 基于信号随偏移距变化特征的叠前道集保幅去噪方法 | |
US11542816B2 (en) | System and method for hydrocarbon exploration based on imaging tunnel valleys | |
CN112083491B (zh) | 基于信号属性特征的井筒波压制方法及装置 | |
CN117492084A (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 |