CN111624659B - 一种地震数据的时变带通滤波方法及装置 - Google Patents

一种地震数据的时变带通滤波方法及装置 Download PDF

Info

Publication number
CN111624659B
CN111624659B CN202010502784.3A CN202010502784A CN111624659B CN 111624659 B CN111624659 B CN 111624659B CN 202010502784 A CN202010502784 A CN 202010502784A CN 111624659 B CN111624659 B CN 111624659B
Authority
CN
China
Prior art keywords
time
data
leveling
seismic data
layer
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
CN202010502784.3A
Other languages
English (en)
Other versions
CN111624659A (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.)
Optical Science and Technology Chengdu Ltd of CNPC
Original Assignee
Optical Science and Technology Chengdu Ltd of CNPC
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 Optical Science and Technology Chengdu Ltd of CNPC filed Critical Optical Science and Technology Chengdu Ltd of CNPC
Priority to CN202010502784.3A priority Critical patent/CN111624659B/zh
Publication of CN111624659A publication Critical patent/CN111624659A/zh
Application granted granted Critical
Publication of CN111624659B publication Critical patent/CN111624659B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

一种地震数据的时变带通滤波方法及装置
技术领域
本发明涉及地球物理勘探的地震数据处理方法，属于地震处理技术领域，特别是涉及一种地震数据的时变带通滤波方法及装置。
背景技术
地震数据反射波的主频随着旅行时增加变低、频宽随着旅行时增加变窄，即频率是时变的。因此，在地震数据处理中经常使用时变带通滤波。通常，时变带通滤波是分段进行的，在段与段之间的过渡带会出现明显的接缝；多数据、多段时变带通不易实现。
发明内容
本发明的目的在于克服现有技术的不足，提供一种地震数据的时变带通滤波方法及装置，根据地震数据频率分布的特点，采用沿层拾取、层拉平、多道数据分段时变带通滤波、层反拉平的方式，实现了多数据、可变分段的时变带通滤波。
本发明的目的是通过以下技术方案来实现的：一种地震数据的时变带通滤波方法，包括以下步骤：
S1.输入N道地震数据；时变带通滤波是对地震数据的中深层开时窗进行的，在时窗上边界附近选取一套连续的层位，对于每道地震数据沿层拾取对应的层位时间，最终得到N道输入地震数据的层位时间集合
Figure BDA0002525369200000011
其中，Ti layer是第i道数据的层位时间，i＝1,2,3,...,N；
S2.确定步骤S1得到层位时间中的最大值，然后确定层拉平时间Tlp，计算层拉平时移量，对输入的地震数据沿层时移拉平；
S3.对时移拉平后的数据做时窗为
Figure BDA0002525369200000012
的多道数据分段时变带通滤波；
S4.利用层位时间
Figure BDA0002525369200000013
拉平时间Tlp，计算层反拉平时移量，对时变带通滤波后的数据做层反拉平。
进一步地，所述步骤S2包括以下子步骤：
S201.确定时间层位
Figure BDA0002525369200000014
的最大值Tmax；
S202.确定拉平时间Tlp，Tlp≥Tmax；
S203.确定每一道地震数据的最大时间
Figure BDA0002525369200000015
S204.在每一道地震数据的最大时间后补零至
Figure BDA0002525369200000016
S205.计算每一道数据的层拉平时移量：
Figure BDA0002525369200000021
其中，Ti Shift是第i道数据的层拉平时移量；
S206.对于每一道地震数据，利用步骤S205得到的时移量时移步骤S204得到的地震数据。
进一步地，所述步骤S3包括以下子步骤：
S301.对时移拉平后的每一道地震数据做傅里叶变换，计算对应的振幅谱A、相位谱φ；
S302.确定时移拉平后各道数据分段时变带通滤波振幅谱At为：
Figure BDA0002525369200000022
其中，fl是低频截止频率，fh是高频截止频率，t是时间，f是频率，
Figure BDA0002525369200000023
是第i道数据的时变带通滤波振幅谱。
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所示，一种地震数据的时变带通滤波方法，包括以下步骤：
S1.输入N道地震数据；时变带通滤波是对地震数据的中深层开时窗进行的，在时窗上边界附近选取一套连续的层位，对于每道地震数据沿层拾取对应的层位时间，最终得到N道输入地震数据的层位时间集合
Figure BDA0002525369200000031
其中，Ti layer是第i道数据的层位时间，i＝1,2,3,...,N；
在本申请的实施例中，输入的地震数据如图2所示，其中横坐标为数据号；纵坐标为时间(单位：ms)；进过步骤S1处理后，沿层拾取的层位时间，如图3所示，图3中，横坐标为数据号；纵坐标为时间(单位：ms)。
S2.确定步骤S1得到层位时间中的最大值，然后确定层拉平时间Tlp，计算层拉平时移量，对输入的地震数据沿层时移拉平：
S201.确定时间层位
Figure BDA0002525369200000032
的最大值Tmax；
S202.确定拉平时间Tlp，Tlp≥Tmax；
S203.确定每一道地震数据的最大时间
Figure BDA0002525369200000033
每一道输入的地震数据是有长度的，如6秒、8秒等，这里的地震数据最大时间，可以理解为地震数据的长度最大值；
S204.在每一道地震数据的最大时间后补零至
Figure BDA0002525369200000034
S205.计算每一道数据的层拉平时移量：
{Ti Shift＝Tlp-Ti layer}，i＝1,2,3,...,N；
其中，Ti Shift是第i道数据的层拉平时移量；
S206.对于每一道地震数据，利用步骤S205得到的时移量时移步骤S204得到的地震数据。
在本申请的实施例中，经过步骤S2处理，得到层拉平地震数据，如图4所示，横坐标为数据号；纵坐标为时间(单位：ms)；
S3.对时移拉平后的数据做时窗为
Figure BDA0002525369200000041
的多道数据分段时变带通滤波：
S301.对时移拉平后的每一道地震数据做傅里叶变换，计算对应的振幅谱A、相位谱φ；
S302.确定时移拉平后各道数据分段时变带通滤波振幅谱At为：
Figure BDA0002525369200000042
其中，fl是低频截止频率，fh是高频截止频率，t是时间，f是频率，
Figure BDA0002525369200000043
是第i道数据的时变带通滤波振幅谱。
S303.对于每一道数据，利用振幅谱At和对应的相位谱φ做反傅里叶变换。
在本申请的实施例中，通过时变带通滤波得到的层拉平数据如图5所示，该图中横坐标为数据号；纵坐标为时间(单位：ms)。
S4.利用层位时间
Figure BDA0002525369200000044
拉平时间Tlp，计算层反拉平时移量，对时变带通滤波后的数据做层反拉平。
S401.对于每一道时变带通滤波后的数据，计算对应的层反拉平时移量：
{Ti InvShift＝Ti layer-Tlp}；i＝1,2,3,...,N；
其中，Ti InvShift是第i道数据的层反拉平时移量；
S402.对于每一道时变带通滤波后的数据，利用层反拉平时移量时移时变带通滤波后的地震数据。
在本申请的实施例中，层反拉平得到的地震数据如图6所示，图中，横坐标为数据号；纵坐标为时间(单位：ms)。
如图7所示，一种地震数据的时变带通滤波装置，包括：
数据输入单元，用于输入多道地震数据；
层位时间拾取单元，用于在时窗上边界附近选取一套连续的层位，对于每道地震数据沿层拾取对应的层位时间，最终得到N道输入地震数据的层位时间集合；
时移拉平单元，用于确定层位时间中的最大值和层拉平时间，并计算层拉平时移量，对输入的地震数据沿层时移拉平；
时变带通滤波单元，用于对时移拉平后的数据做多道数据分段时变带通滤波；
层反拉平单元，用于根据层位时间和拉平时间，计算层反拉平时移量，对时变带通滤波后的数据做层反拉平。
综上，本发明根据地震数据频率分布的特点，采用沿层拾取、层拉平、多道数据分段时变带通滤波、层反拉平的方式，实现了多数据、可变分段的时变带通滤波，便于地震资料后续处理。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其他不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (4)

1.一种地震数据的时变带通滤波方法，其特征在于：包括以下步骤：
S1.输入N道地震数据；时变带通滤波是对地震数据的中深层开时窗进行的，在时窗上边界附近选取一套连续的层位，对于每道地震数据沿层拾取对应的层位时间，最终得到N道输入地震数据的层位时间集合
Figure FDA0003646617810000011
其中，Ti layer是第i道数据的层位时间，i＝1,2,3,...,N；
S2.确定步骤S1得到层位时间中的最大值，然后确定层拉平时间Tlp，计算层拉平时移量，
对输入的地震数据沿层时移拉平；
所述步骤S2包括以下子步骤：
S201.确定时间层位
Figure FDA0003646617810000012
的最大值Tmax；
S202.确定拉平时间Tlp，Tlp≥Tmax；
S203.确定每一道地震数据的最大时间
Figure FDA0003646617810000013
S204.在每一道地震数据的最大时间后补零至
Figure FDA0003646617810000014
S205.计算每一道数据的层拉平时移量：
{Ti Shift＝Tlp-Ti layer}，i＝1,2,3,...,N；
其中，Ti Shift是第i道数据的层拉平时移量；
S206.对于每一道地震数据，利用步骤S205得到的时移量时移步骤S204得到的地震数据；
S3.对时移拉平后的数据做时窗为
Figure FDA0003646617810000015
的多道数据分段时变带通滤波；
S4.利用层位时间
Figure FDA0003646617810000016
拉平时间Tlp，计算层反拉平时移量，对时变带通滤波后的数据做层反拉平。
2.根据权利要求1所述的一种地震数据的时变带通滤波方法，其特征在于：所述步骤S3包括以下子步骤：
S301.对时移拉平后的每一道地震数据做傅里叶变换，计算对应的振幅谱A、相位谱φ；
S302.确定时移拉平后各道数据分段时变带通滤波振幅谱At为：
Figure FDA0003646617810000017
其中，fl是低频截止频率，fh是高频截止频率，t是时间，f是频率，Ait是第i道数据的时变带通滤波振幅谱；
S303.对于每一道数据，利用振幅谱At和对应的相位谱φ做反傅里叶变换。
3.根据权利要求1所述的一种地震数据的时变带通滤波方法，其特征在于：所述步骤S4包括以下子步骤：
S401.对于每一道时变带通滤波后的数据，计算对应的层反拉平时移量：
{Ti InvShift＝Ti layer-Tlp}；i＝1,2,3,...,N；
其中，Ti InvShift是第i道数据的层反拉平时移量；
S402.对于每一道时变带通滤波后的数据，利用层反拉平时移量时移时变带通滤波后的地震数据。
4.一种地震数据的时变带通滤波装置，采用如权利要求1～3中任意一项所述的方法，其特征在于：包括：
数据输入单元，用于输入多道地震数据；
层位时间拾取单元，用于在时窗上边界附近选取一套连续的层位，对于每道地震数据沿层拾取对应的层位时间，最终得到N道输入地震数据的层位时间集合；
时移拉平单元，用于确定层位时间中的最大值和层拉平时间，并计算层拉平时移量，对输入的地震数据沿层时移拉平；
时变带通滤波单元，用于对时移拉平后的数据做多道数据分段时变带通滤波；
层反拉平单元，用于根据层位时间和拉平时间，计算层反拉平时移量，对时变带通滤波后的数据做层反拉平。
CN202010502784.3A 2020-06-05 2020-06-05 一种地震数据的时变带通滤波方法及装置 Active CN111624659B (zh)

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 CN111624659A (zh) 2020-09-04
CN111624659B true 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 (6)

* Cited by examiner, † Cited by third party
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
CN102778693A (zh) * 2011-05-13 2012-11-14 中国石油化工股份有限公司 一种基于反射波层拉平提取并消除的绕射波分离处理方法
CN104267431A (zh) * 2014-10-13 2015-01-07 成都北方石油勘探开发技术有限公司 一种地震构造解释的小断层获取方法
CN104280765A (zh) * 2013-07-11 2015-01-14 中国石油化工股份有限公司 基于变子波反射系数反演的地震高分辨处理方法
CN105629300A (zh) * 2015-12-30 2016-06-01 中国石油集团川庆钻探工程有限公司地球物理勘探公司 提高复杂构造偏移数据信噪比的方法
CN107942387A (zh) * 2017-10-23 2018-04-20 中国石油天然气集团公司 一种续至波衰减处理方法及装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004028034B4 (de) * 2004-06-09 2006-11-02 Ernst D. Rode Verfahren zur Bestimmung der Tiefe und der Mächtigkeit von unterirdischen Kohlenwasserstoffvorkommen mit einem oder mehreren Reservoirhorizonten durch an der Oberfläche angeordnete Empfänger für akustische Wellen in einem Frequenzbereich von 0,2 bis 30 Hz
CN102109612B (zh) * 2009-12-23 2012-10-17 中国石油天然气集团公司 一种地震波吸收衰减补偿方法
US8612156B2 (en) * 2010-03-05 2013-12-17 Vialogy Llc Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization
US9105075B1 (en) * 2013-02-06 2015-08-11 Ihs Global Inc. Enhancing seismic features using an optical filter array
CN105093324B (zh) * 2014-05-21 2017-06-16 中国石油化工股份有限公司 一种剩余静校正方法
CN106709566A (zh) * 2016-12-07 2017-05-24 广东工业大学 一种基于深度学习的数据缺失值回填方法
CN106780347B (zh) * 2017-02-09 2020-03-03 浙江科技学院 一种基于oct图像处理的枇杷早期瘀伤鉴别方法
CN107607992B (zh) * 2017-08-24 2020-08-18 电子科技大学 基于卷积神经网络的多波匹配方法
CN108571000A (zh) * 2018-05-15 2018-09-25 中铁第六勘察设计院集团有限公司 一种沿河地段道路下穿既有铁路时地下水路堑的施工方法
CN109829398B (zh) * 2019-01-16 2020-03-31 北京航空航天大学 一种基于三维卷积网络的视频中的目标检测方法
CN110399798B (zh) * 2019-06-25 2021-07-20 朱跃飞 一种基于深度学习的离散图片文件信息提取系统及方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
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
CN102778693A (zh) * 2011-05-13 2012-11-14 中国石油化工股份有限公司 一种基于反射波层拉平提取并消除的绕射波分离处理方法
CN104280765A (zh) * 2013-07-11 2015-01-14 中国石油化工股份有限公司 基于变子波反射系数反演的地震高分辨处理方法
CN104267431A (zh) * 2014-10-13 2015-01-07 成都北方石油勘探开发技术有限公司 一种地震构造解释的小断层获取方法
CN105629300A (zh) * 2015-12-30 2016-06-01 中国石油集团川庆钻探工程有限公司地球物理勘探公司 提高复杂构造偏移数据信噪比的方法
CN107942387A (zh) * 2017-10-23 2018-04-20 中国石油天然气集团公司 一种续至波衰减处理方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
地震勘探中的去噪技术新进展;张孝珍;《勘探地球物理进展》;20090625;全文 *

Also Published As

Publication number Publication date
CN111624659A (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
CN109061764B (zh) 一种分频融合波阻抗反演方法
CN103728659B (zh) 一种提高地下岩溶探测精度的方法
EA020070B1 (ru) Способ индексирования подземного участка с целью получения геологической информации
CN105116448B (zh) 一种转换波方位各向异性校正方法及装置
CN108845357B (zh) 一种基于同步挤压小波变换估计地层等效品质因子的方法
CN105116443B (zh) 一种低频信号的能量补偿方法及装置
CN103645507A (zh) 地震记录的处理方法
CN109765615A (zh) 一种地层品质因子反演方法及装置
EP1914568A2 (en) Method for determining impedence coefficients of a seismic trace
CN108427140B (zh) 一种用于小尺度缝洞储集体地震识别的方法
CN111624659B (zh) 一种地震数据的时变带通滤波方法及装置
CN105093318B (zh) 一种自适应波动方程波场延拓静校正方法
CN110568490B (zh) 一种高速层顶薄储层的识别方法
CN110618457B (zh) 一种可控震源地震数据的联合去噪方法及其装置
CN104062681B (zh) 一种基于分数阶导数的地震层位追踪预处理方法
CN112213776B (zh) 叠前道集和vsp资料联合层控q模型建立方法
CN102012520B (zh) 一种海上单点采集地震数据室内动态组合的方法
CN110471105B (zh) 一种层位约束下基于地震解析道的相干增强断层识别方法
CN112630832B (zh) 基于衰减因子随入射角变化的含气性预测方法和装置
CN109655885A (zh) 一种自动初至拾取方法及系统
CN111045083A (zh) 一种储层含气性检测方法
CN111060966B (zh) 一种用于静校正技术的初至拾取方法
CN111694055A (zh) 一种强反射屏蔽下avo属性提取和分析的方法和装置
CN111060970B (zh) 一种储层约束反褶积的方法及应用
CN112083491B (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