CN109557587A - 一种vsp地震资料井筒波频率域滤波方法及装置 - Google Patents

一种vsp地震资料井筒波频率域滤波方法及装置 Download PDF

Info

Publication number
CN109557587A
CN109557587A CN201811618332.0A CN201811618332A CN109557587A CN 109557587 A CN109557587 A CN 109557587A CN 201811618332 A CN201811618332 A CN 201811618332A CN 109557587 A CN109557587 A CN 109557587A
Authority
CN
China
Prior art keywords
amplitude spectrum
peak
data
channel record
statistical model
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
Application number
CN201811618332.0A
Other languages
English (en)
Other versions
CN109557587B (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.)
Yangtze University
Original Assignee
Yangtze University
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 Yangtze University filed Critical Yangtze University
Priority to CN201811618332.0A priority Critical patent/CN109557587B/zh
Publication of CN109557587A publication Critical patent/CN109557587A/zh
Application granted granted Critical
Publication of CN109557587B publication Critical patent/CN109557587B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/30—Noise handling
    • G01V2210/32—Noise reduction
    • G01V2210/324—Filtering

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

本发明提出一种VSP地震资料井筒波频率域滤波方法及装置，所述方法包括：采集VSP地震资料原始数据，使用初至时间对所述原始数据进行下行波拉平预处理；采用傅里叶变换将预处理后的数据变换到频率域，进行频谱分析得到每一道记录的振幅谱；对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；对处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；在时间域上对所述使用统计模型替换振幅谱的单道记录进行能量幅度修正；输出资料数据。本发明提供的方法操作简单、运算高效，可减少余波及假频等人为主观性的干扰，属于地震勘探资料处理技术领域。

Description

一种VSP地震资料井筒波频率域滤波方法及装置
技术领域
本发明属于地震勘探资料处理技术领域，主要涉及一种VSP地震资料中井筒波的频率域滤波方法及装置。
背景技术
地震勘探中的VSP(垂直地震剖面)资料由于特殊的采集方式和固井质量差异，数据采集时经常受到井筒波干扰，使得资料信噪比降低，对后续处理造成很大影响。
当前常用的井筒波压制方法有中值滤波和F-K滤波(频率-波数域滤波)。中值滤波具有较强的统计性，需要准确的井筒波初至时间，再按照初至对井筒波进行滤除，其缺点在于滤除井筒波的同时会改变原始数据的相位和能量关系，且会留下较强的残余井筒波。F-K滤波利用视速度的差异，清零井筒波视速度范围内的F-K谱以滤除井筒波，其缺点在于无法消除当F-K谱中的波数大于空间采样中的奈奎斯特波数时所产生的假频。
发明内容
鉴于此，本发明提供一种操作简单，运算高效的VSP地震资料中井筒波的频率域滤波方法及装置，以去除人为主观性的噪声干扰。
本发明提出一种VSP地震资料井筒波频率域滤波方法，所述方法包括：
S1、采集VSP地震资料原始数据，使用初至时间对所述原始数据进行下行波拉平预处理；
S2、采用傅里叶变换将所述预处理后的数据变换到频率域，进行频谱分析得到每一道记录的振幅谱；
S3、对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；
S4、对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；
S5、对处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；
S6、在时间域上根据能量随初至时间呈指数衰减的理论关系，对步骤S4中所述使用统计模型替换振幅谱的单道记录进行能量幅度修正；
S7、输出分离了井筒波的VSP资料数据。
可选的，步骤S4所述对单道记录振幅谱进行修正的具体过程为：
S41、确定所述统计模型主频范围、主频峰值MPV(Model Peak Value)及峰值位置MPP(Model Peak Position)；得到归一化处理后的单道记录振幅谱；
S42、确定所述统计模型主频范围内单道记录振幅谱的峰值SPV(Single PeakValue)及峰值位置SPP(Single Peak Position)；
S43、对SPP进行判断，如果SPP偏离MPP的距离超过预设阈值，则使用统计模型替换单道记录道振幅谱；
S44、如果SPP偏离MPP的距离在预设阈值范围内，计算SPV与MPV的比值Rmax，按照Rmax来修正单道记录振幅谱上其他异常峰值。
可选的，所述步骤S44中，所述按照Rmax来修正单道记录振幅谱上其他异常峰值的具体过程为：
S441、根据所述单道记录振幅谱的峰值SPV及峰值位置SPP寻找左右边界，确定SPV对应的峰值区间，保存所述峰值区间；
S442、除开所有保存过的峰值区间，在单道记录振幅谱其他范围内寻找最大峰值SMPV以及峰值位置SMPP，并确定峰值区间C，保存峰值区间C；
S443、确定统计模型在区间C上的峰值MMPV，并计算SMPV与MMPV的比值R；
S444、预设一个经验阈值DT(Difference Threshold)，计算比值R和Rmax的相对差异，如果相对差异大于所述经验阈值DT，则使用比值Rmax来修正单道记录在区间C内的振幅谱；如果相对差异不超过所述经验阈值DT，则不做任何处理；
S445、重复步骤S442至步骤S444，直到单道记录振幅谱上所有异常峰值修正完毕。
可选的，所述步骤S6中，所述对使用统计模型替换振幅谱的道进行能量幅度修正的具体过程为：
S61、在时域波形中选定包含前三个波峰的时间范围TS；
S62、对每道记录在TS范围内的数据求平方和，得到能量振幅向量A，并根据记录道初至时间得出初至时间向量t，A＝A0ta，其中A0为初至时间最小的记录道对应的能量振幅，a为衰减因子；
S63、分别对A和初至时间向量t进行对数运算，通过一次线性拟合得到每道记录能量理论大小的对数向量LA；
S64、确定通过统计模型替换振幅谱的单道记录n，找到A中对应位置的能量值，记为E；
S65、对LA中对应位置的值进行指数运算，记为LE，并计算LE与E的比值B；
S66、使用B对单道记录n进行能量修正；
S67、重复步骤S64至步骤S66，直到所有通过统计模型替换振幅谱的单道记录修正完毕。
本发明还提供一种VSP地震资料井筒波频率域滤波装置，所述装置包括：
采集处理模块：采集VSP地震资料原始数据，使用初至时间对所述原始数据进行下行波拉平预处理；采用傅里叶变换将所述预处理后的数据变换到频率域，频谱分析得到每一道记录的振幅谱；
统计建模模块：对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；
频域修正模块：对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；
时域变换模块：对经过频域修正模块处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；
时域修正模块：在时间域上根据能量随初至时间呈指数衰减的理论关系，对使用统计模型替换振幅谱的道进行能量幅度修正；输出分离了井筒波的VSP资料数据。
可选的，所述频域修正模块包括：
极值确定单元：确定所述统计模型主频范围、主频峰值MPV及峰值位置MPP；得到归一化处理后的单道记录振幅谱；确定所述统计模型主频范围内单道记录振幅谱的峰值SPV及峰值位置SPP；
位置判断单元：对SPP进行判断，判断SPP偏离MPP的距离是否超过预设阈值；
替换修正单元：如果SPP偏离MPP的距离超过预设阈值，则使用统计模型替换单道记录道振幅谱；
比值修正单元：
如果SPP偏离MPP的距离在预设阈值范围内，计算SPV与MPV的比值Rmax，确定峰值SPV对应的峰值区间，保存所述峰值区间，按照Rmax来修正单道记录振幅谱上的异常峰值；
所述按照Rmax来修正单道记录振幅谱上的异常峰值的过程为：除去所有保存过的峰值区间，在单道记录振幅谱的其他范围内寻找最大峰值SMPV以及峰值位置SMPP，并确定峰值区间C，保存峰值区间C；确定所述统计模型在所述峰值区间C上的峰值MMPV，并计算SMPV与MMPV的比值R；预设一个经验阈值DT，计算比值R和Rmax的相对差异，如果所述相对差异大于所述经验阈值DT，则使用比值Rmax来修正单道记录在区间C内的振幅谱；如果所述相对差异小于等于所述经验阈值DT，则不做任何处理；
重复所述按照Rmax来修正单道记录振幅谱上的异常峰值的过程，直到单道记录振幅谱上所有异常峰值修正完毕
可选的，所述时域修正模块包括：
计算单元：在时域波形中选定包含前三个波峰的时间范围TS；对每道记录在TS范围内的数据求平方和，得到能量振幅向量A，并根据记录道初至时间得出初至时间向量t，VSP地震资料总的振幅衰减公式为：A＝A0ta，其中A0为初至时间最小的记录道对应的能量振幅，a为衰减因子；分别对A和初至时间向量t进行对数运算，通过一次线性拟合得到每道记录能量理论大小的对数向量LA；
能量修正单元：确定通过统计模型替换振幅谱的单道记录n，找到A中对应位置的能量值，记为E；对LA中对应位置的值进行指数运算，记为LE，并计算LE与E的比值B，使用B对单道记录n进行能量修正；重复能量修正单元，直到所有通过统计模型替换振幅谱的单道记录修正完毕。
本发明的有益效果是：
本发明提供一种VSP地震资料中井筒波的频率域滤波方法及装置，与现有技术相比，本发明提供的方法操作简单、运算高效，可减少余波及假频等人为主观性的干扰。人为主观性指通过滤波器进行频域滤波时，随着不同记录道频谱上需要滤波的频段不同而设置不同的滤波器指标，就需要人为的去寻找这些不同的需要滤除的频段，但现有滤波技术往往会带来余波、假频等干扰导致信息畸变或失真，本发明不用滤波器，直接对频谱进行修正，通过频率域滤波来逐个修正异常峰值，因此可排除人为主观性的干扰。
附图说明
为了更清楚地说明本发明的技术方案，下面将对本发明技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的井筒波频率域滤波方法流程示意图；
图2为本发明提供的VSP原始下行波拉平数据；
图3(a)为本发明提供的振幅谱统计模型；
图3(b)为本发明提供的主频畸变的含井筒波单道记录振幅谱；
图3(c)为本发明提供的含井筒波的单道记录振幅谱谱；
图3(d)为本发明提供的含井筒波的单道记录振幅谱谱修正前后对比图；
图4为本发明提供的含井筒波单道记录异常峰值修正流程示意图；
图5为本发明提供的井筒波分离后的时域数据能量修正流程示意图；
图6为本发明提供的井筒波分离之后的数据；
图7为本发明提供的井筒波频率域滤波装置结构示意图。
具体实施方式
本发明提供了一种VSP地震资料井筒波频率域滤波方法及装置，用于对VSP地震资料井筒波进行频率域滤波，减少噪声干扰。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
请参见图1，本发明提供的VSP地震资料井筒波频率域滤波方法整体流程示意图，所述方法具体包括：
S1、采集VSP地震资料原始数据，使用初至时间对所述原始数据进行下行波拉平预处理；
具体的，请参见图2，图2是使用初至时间对VSP原始数据下行波拉平处理后的数据。
S2、采用傅里叶变换将所述拉平预处理后的数据变换到频率域，进行频谱分析得到每一道记录的振幅谱；
S3、对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；
具体的，请参见图3(a)，图3(a)为归一化均值处理后的振幅谱统计模型。
S4、对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；
步骤S4所述对单道记录振幅谱进行修正的具体过程为：
S41、确定所述统计模型主频范围、主频峰值MPV(Model Peak Value)及峰值位置MPP(Model Peak Position)；得到归一化处理后的单道记录振幅谱；
S42、确定所述统计模型主频范围内单道记录振幅谱的峰值SPV(Single PeakValue)及峰值位置SPP(Single Peak Position)；
S43、对SPP进行判断，如果SPP偏离MPP的距离超过预设阈值，则使用统计模型替换单道记录道振幅谱；具体的，如果由于井筒波过强导致的主频畸变使得SPP偏离MPP过远，如图3(b)所示，则使用统计模型替换单道记录道振幅谱，完成单道记录振幅谱修正。
S44、如果SPP偏离MPP的距离在预设阈值范围内，计算SPV与MPV的比值Rmax，按照Rmax来修正单道记录振幅谱上其他异常峰值。
上述步骤S44中，所述按照Rmax来修正单道记录振幅谱上其他异常峰值的流程示意图请参见图4，修正异常峰值的具体过程为：
S441、根据所述单道记录振幅谱的峰值SPV及峰值位置SPP寻找左右边界，确定SPV对应的峰值区间，保存所述峰值区间；
S442、除开所有保存过的峰值区间，在单道记录振幅谱其他范围内寻找最大峰值SMPV以及峰值位置SMPP，并确定峰值区间C，保存峰值区间C；
S443、确定所述统计模型在所述峰值区间C上的峰值MMPV，并计算SMPV与MMPV的比值R；
S444、预设一个经验阈值DT(Difference Threshold)，计算比值R和Rmax的相对差异，如果所述相对差异大于所述经验阈值DT，则使用比值Rmax来修正单道记录在区间C内的振幅谱；如果相对差异小于等于所述经验阈值DT，则不做任何处理；
S445、重复步骤S442至步骤S444，直到单道记录振幅谱上所有异常峰值修正完毕。具体的，图3(c)所示为含井筒波单道记录振幅谱，图3(d)为含井筒波的单道记录振幅谱修正前后对比，由图3(d)可知，本实施例中所述按照Rmax来修正单道记录振幅谱上异常峰值的方法对异常峰值具有极好的处理效果。
S5、对处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；
S6、在时间域上根据能量随初至时间呈指数衰减的理论关系，对使用统计模型替换振幅谱的道进行能量幅度修正；
步骤S6中所述对使用统计模型替换振幅谱的道进行能量幅度修正的具体过程为：
S61、在时域波形中选定包含前三个波峰的时间范围TS；
S62、对每道记录在TS范围内的数据求平方和，得到能量振幅向量A，并根据记录道初至时间得出初至时间向量t，A＝A0ta，其中A0为初至时间最小的记录道对应的能量振幅，a为衰减因子；
S63、分别对A和初至时间向量t进行对数运算，通过一次线性拟合得到每道记录能量理论大小的对数向量LA；
具体的，VSP地震资料总的振幅衰减公式为：A＝A0ta，对振幅衰减公式两边进行对数运算，得到公式：lnA＝alnt+ln A0，通过一次线性拟合结合初至时间向量t得到每道记录能量理论大小的对数向量LA；
S64、确定通过统计模型替换振幅谱的单道记录n，找到A中对应位置的能量值，记为E；
S65、对LA中对应位置的值进行指数运算，记为LE，并计算LE与E的比值B；
S66、使用B对单道记录n进行能量修正；
S67、重复步骤S64至步骤S66，直到所有通过统计模型替换振幅谱的单道记录修正完毕。
S7、输出分离了井筒波的VSP资料数据。
请参见图6，图6为采用本发明所述的井筒波频率域滤波方法得到的井筒波分离之后的数据，与图2的VSP原始下行波拉平数据比较可知，本发明所述的井筒波频率域滤波方法可较大程度地排除干扰，有较好的滤波效果。
请参见图7，本发明还提供一种VSP地震资料井筒波频率域滤波装置，所述装置包括：
采集处理模块710：采集VSP地震资料原始数据，使用初至时间对数据进行下行波拉平预处理；采用傅里叶变换将所述预处理后的数据变换到频率域，进行频谱分析得到每一道记录的振幅谱；
统计建模模块720：对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；
频域修正模块730：对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；
时域变换模块740：对处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；
时域修正模块750：在时间域上根据能量随初至时间呈指数衰减的理论关系，对使用统计模型替换振幅谱的道进行能量幅度修正；输出分离了井筒波的VSP资料数据。
所述频域修正模块730包括：
极值确定单元7310：确定所述统计模型主频范围、主频峰值MPV及峰值位置MPP；得到归一化处理后的单道记录振幅谱；确定所述统计模型主频范围内单道记录振幅谱的峰值SPV及峰值位置SPP；
位置判断单元7320：对SPP进行判断，判断SPP偏离MPP的距离是否超过预设阈值；
替换修正单元7330：如果SPP偏离MPP的距离超过预设阈值，则使用统计模型替换单道记录道振幅谱；
比值修正单元7340：如果SPP偏离MPP的距离在预设阈值范围内，计算SPV与MPV的比值Rmax，确定峰值SPV对应的峰值区间，保存所述峰值区间，按照Rmax来修正单道记录振幅谱上的异常峰值；
所述按照Rmax来修正单道记录振幅谱上的异常峰值的过程为：除去所有保存过峰值区间，在单道记录振幅谱的其他范围内寻找最大峰值SMPV以及峰值位置SMPP，并确定峰值区间C，保存峰值区间C；确定所述统计模型在所述峰值区间C上的峰值MMPV，并计算SMPV与MMPV的比值R；预设一个经验阈值DT，计算比值R和Rmax的相对差异，如果所述相对差异大于所述经验阈值DT，则使用比值Rmax来修正单道记录在区间C内的振幅谱；如果所述相对差异小于等于所述经验阈值DT，则不做任何处理；
重复所述按照Rmax来修正单道记录振幅谱上的异常峰值的过程，直到单道记录振幅谱上所有异常峰值修正完毕
所述时域修正模块750包括：
计算单元7510：在时域波形中选定包含前三个波峰的时间范围TS；对每道记录在TS范围内的数据求平方和，得到能量振幅向量A，并根据记录道初至时间得出初至时间向量t，A＝A0ta，其中A0为初至时间最小的记录道对应的能量振幅，a为衰减因子；分别对A和初至时间向量t进行对数运算，通过一次线性拟合得到每道记录能量理论大小的对数向量LA；
能量修正单元7520：确定通过统计模型替换振幅谱的单道记录n(n为道号)，找到A中对应位置的能量值，记为E；对LA中对应位置的值进行指数运算，记为LE，并计算LE与E的比值B，使用B对单道记录n进行能量修正；重复能量修正单元，直到所有通过统计模型替换振幅谱的单道记录修正完毕。
本发明提供一种操作简单、运算高效的VSP地震资料中井筒波的频率域滤波方法及装置，可减少余波及假频等人为主观性的干扰。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各实施例的模块、单元和/或方法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式。
以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种VSP地震资料井筒波频率域滤波方法，其特征在于，所述方法包括：
S1、采集VSP地震资料原始数据，使用初至时间对所述原始数据进行下行波拉平预处理；
S2、采用傅里叶变换将所述预处理后的数据变换到频率域，进行频谱分析得到每一道记录的振幅谱；
S3、对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；
S4、对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；
S5、对处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；
S6、在时间域上根据能量随初至时间呈指数衰减的理论关系对数据进行能量幅度修正；
S7、输出分离了井筒波的VSP资料数据。
2.根据权利要求1所述一种VSP地震资料井筒波频率域滤波方法，其特征在于，所述步骤S4中，所述对单道记录振幅谱进行修正的具体过程为：
S41、确定所述统计模型主频范围、主频峰值MPV及峰值位置MPP；得到归一化处理后的单道记录振幅谱；
S42、确定所述统计模型主频范围内单道记录振幅谱的峰值SPV及峰值位置SPP；
S43、对SPP进行判断，如果SPP偏离MPP的距离超过预设阈值，则使用统计模型替换所述单道记录道振幅谱；
S44、如果SPP偏离MPP的距离在预设阈值范围内，计算SPV与MPV的比值Rmax，按照Rmax来修正单道记录振幅谱上其他异常峰值。
3.根据权利要求2所述一种VSP地震资料井筒波频率域滤波方法，其特征在于，所述步骤S44中，所述按照Rmax来修正单道记录振幅谱上其他异常峰值的具体过程为：
S441、根据所述单道记录振幅谱的峰值SPV及峰值位置SPP寻找左右边界，确定SPV对应的峰值区间，保存所述峰值区间；
S442、除开所有保存过的峰值区间，在单道记录振幅谱其他范围内寻找最大峰值SMPV以及峰值位置SMPP，并确定对应的峰值区间C，保存峰值区间C；
S443、确定所述统计模型在所述峰值区间C上的峰值MMPV，并计算SMPV与MMPV的比值R；
S444、预设一个经验阈值DT，计算比值R和Rmax的相对差异，如果所述相对差异大于所述经验阈值DT，则使用比值Rmax来修正单道记录在区间C内的振幅谱；如果所述相对差异小于等于所述经验阈值DT，则不做任何处理；
S445、重复步骤S442至步骤S444，直到单道记录振幅谱上所有异常峰值修正完毕。
4.根据权利要求2所述一种VSP地震资料井筒波频率域滤波方法，其特征在于，所述步骤S6的具体过程为：
S61、在时域波形中选定包含前三个波峰的时间范围TS；
S62、对每道记录在TS范围内的数据求平方和，得到能量振幅向量A，并根据记录道初至时间得出初至时间向量t，VSP地震资料总的振幅衰减公式为A＝A0ta，其中A0为初至时间最小的记录道对应的能量振幅，a为衰减因子；
S63、分别对A和初至时间向量t进行对数运算，通过一次线性拟合得到每道记录能量理论大小的对数向量LA；
S64、确定通过统计模型替换振幅谱的单道记录n，找到A中对应位置的能量值，记为E；
S65、对LA中对应位置的值进行指数运算，记为LE，并计算LE与E的比值B；
S66、使用B对单道记录n进行能量修正；
S67、重复步骤S64至步骤S66，直到所有通过统计模型替换振幅谱的单道记录修正完毕。
5.一种VSP地震资料井筒波频率域滤波装置，其特征在于，所述装置包括：
采集处理模块：采集VSP地震资料原始数据，使用初至时间对所述原始数据进行下行波拉平预处理；采用傅里叶变换将所述预处理后的数据变换到频率域，频谱分析得到每一道记录的振幅谱；
统计建模模块：对所有记录道振幅谱进行归一化均值处理，建立振幅谱统计模型；
频域修正模块：对比所述振幅谱统计模型，对单道记录振幅谱进行修正，直到所有记录处理完毕；
时域变换模块：对经过频域修正模块处理完的频率域数据进行傅里叶反变换，得到分离井筒波的时间域数据；
时域修正模块：在时间域上根据能量随初至时间呈指数衰减的理论关系，对使用统计模型替换振幅谱的道进行能量幅度修正；输出分离了井筒波的VSP资料数据。
6.根据权利要求5所述一种VSP地震资料井筒波频率域滤波装置，其特征在于，所述频域修正模块包括：
极值确定单元：确定所述统计模型主频范围、主频峰值MPV及峰值位置MPP；得到归一化处理后的单道记录振幅谱；确定所述统计模型主频范围内单道记录振幅谱的峰值SPV及峰值位置SPP；
位置判断单元：对SPP进行判断，判断SPP偏离MPP的距离是否超过预设阈值；
替换修正单元：如果SPP偏离MPP的距离超过预设阈值，则使用统计模型替换单道记录道振幅谱；
比值修正单元：如果SPP偏离MPP的距离在预设阈值范围内，计算SPV与MPV的比值Rmax，确定峰值SPV对应的峰值区间，保存所述峰值区间，按照Rmax来修正单道记录振幅谱上的异常峰值；
所述按照Rmax来修正单道记录振幅谱上的异常峰值的过程为：除去所有保存过的峰值区间，在单道记录振幅谱的其他范围内寻找最大峰值SMPV以及峰值位置SMPP，并确定峰值区间C，保存峰值区间C；确定所述统计模型在所述峰值区间C上的峰值MMPV，并计算SMPV与MMPV的比值R；预设一个经验阈值DT，计算比值R和Rmax的相对差异，如果所述相对差异大于所述经验阈值DT，则使用比值Rmax来修正单道记录在区间C内的振幅谱；如果所述相对差异小于等于所述经验阈值DT，则不做任何处理；
重复所述按照Rmax来修正单道记录振幅谱上的异常峰值的过程，直到单道记录振幅谱上所有异常峰值修正完毕。
7.根据权利要求5所述一种VSP地震资料井筒波频率域滤波装置，其特征在于，所述时域修正模块包括：
计算单元：在时域波形中选定包含前三个波峰的时间范围TS；对每道记录在TS范围内的数据求平方和，得到能量振幅向量A，并根据记录道的初至时间得出初至时间向量t，VSP地震资料总的振幅衰减公式为：A＝A0ta，其中A0为初至时间最小的记录道对应的能量振幅，a为衰减因子；分别对A和初至时间向量t进行对数运算，通过一次线性拟合得到每道记录能量理论大小的对数向量LA；
能量修正单元：确定通过统计模型替换振幅谱的单道记录n，找到A中对应位置的能量值，记为E；对LA中对应位置的值进行指数运算，记为LE，并计算LE与E的比值B，使用B对单道记录n进行能量修正；重复能量修正单元，直到所有通过统计模型替换振幅谱的单道记录修正完毕。
CN201811618332.0A 2018-12-28 2018-12-28 一种vsp地震资料井筒波频率域滤波方法及装置 Active CN109557587B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811618332.0A CN109557587B (zh) 2018-12-28 2018-12-28 一种vsp地震资料井筒波频率域滤波方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811618332.0A CN109557587B (zh) 2018-12-28 2018-12-28 一种vsp地震资料井筒波频率域滤波方法及装置

Publications (2)

Publication Number Publication Date
CN109557587A true CN109557587A (zh) 2019-04-02
CN109557587B CN109557587B (zh) 2020-10-30

Family

ID=65871646

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811618332.0A Active CN109557587B (zh) 2018-12-28 2018-12-28 一种vsp地震资料井筒波频率域滤波方法及装置

Country Status (1)

Country Link
CN (1) CN109557587B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111965704A (zh) * 2019-05-20 2020-11-20 中国石油天然气集团有限公司 井中地震数据的光缆次生震荡噪声压制方法及系统
CN112083491A (zh) * 2019-06-12 2020-12-15 中国石油天然气集团有限公司 基于信号属性特征的井筒波压制方法及装置
CN112099088A (zh) * 2020-09-16 2020-12-18 中油奥博(成都)科技有限公司 一种基于高密度光纤地震数据的油气指示及表征方法
CN113219539A (zh) * 2020-02-05 2021-08-06 中国石油天然气集团有限公司 光纤声波传感地震数据中井筒波干扰的去除方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060203614A1 (en) * 2005-03-09 2006-09-14 Geo-X Systems, Ltd. Vertical seismic profiling method utilizing seismic communication and synchronization
CN101852864A (zh) * 2009-04-03 2010-10-06 中国石油集团东方地球物理勘探有限责任公司 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法
CN102692650A (zh) * 2011-03-23 2012-09-26 中国石油天然气集团公司 一种具有假频压制功能的井筒波分离方法
EP2548052A2 (en) * 2010-03-16 2013-01-23 BP Corporation North America Inc. System and method of 3d salt flank vsp imaging with transmitted waves
CN104090299A (zh) * 2014-07-16 2014-10-08 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于vsp初至波的地面地震数据振幅补偿方法
CN104459784A (zh) * 2014-12-11 2015-03-25 中国科学院地质与地球物理研究所 基于单台、双台和双事件数据二维Lg波Q值层析成像方法
US20160131780A1 (en) * 2004-02-06 2016-05-12 Sigma Cubed, Inc. Method for processing borehole seismic data
CN106291697A (zh) * 2015-06-26 2017-01-04 中国石油化工股份有限公司 一种确定地层品质因子q的值的方法及系统
CN107305223A (zh) * 2016-04-19 2017-10-31 天津大学 一种改进的相位差频率估计方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160131780A1 (en) * 2004-02-06 2016-05-12 Sigma Cubed, Inc. Method for processing borehole seismic data
US20060203614A1 (en) * 2005-03-09 2006-09-14 Geo-X Systems, Ltd. Vertical seismic profiling method utilizing seismic communication and synchronization
CN101852864A (zh) * 2009-04-03 2010-10-06 中国石油集团东方地球物理勘探有限责任公司 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法
EP2548052A2 (en) * 2010-03-16 2013-01-23 BP Corporation North America Inc. System and method of 3d salt flank vsp imaging with transmitted waves
CN102692650A (zh) * 2011-03-23 2012-09-26 中国石油天然气集团公司 一种具有假频压制功能的井筒波分离方法
CN104090299A (zh) * 2014-07-16 2014-10-08 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于vsp初至波的地面地震数据振幅补偿方法
CN104459784A (zh) * 2014-12-11 2015-03-25 中国科学院地质与地球物理研究所 基于单台、双台和双事件数据二维Lg波Q值层析成像方法
CN106291697A (zh) * 2015-06-26 2017-01-04 中国石油化工股份有限公司 一种确定地层品质因子q的值的方法及系统
CN107305223A (zh) * 2016-04-19 2017-10-31 天津大学 一种改进的相位差频率估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
凌云 等: "零井源距VSP上行波统计反褶积方法研究", 《石油物探》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111965704A (zh) * 2019-05-20 2020-11-20 中国石油天然气集团有限公司 井中地震数据的光缆次生震荡噪声压制方法及系统
CN111965704B (zh) * 2019-05-20 2023-09-26 中国石油天然气集团有限公司 井中地震数据的光缆次生震荡噪声压制方法及系统
CN112083491A (zh) * 2019-06-12 2020-12-15 中国石油天然气集团有限公司 基于信号属性特征的井筒波压制方法及装置
CN112083491B (zh) * 2019-06-12 2023-10-31 中国石油天然气集团有限公司 基于信号属性特征的井筒波压制方法及装置
CN113219539A (zh) * 2020-02-05 2021-08-06 中国石油天然气集团有限公司 光纤声波传感地震数据中井筒波干扰的去除方法及装置
WO2021155754A1 (zh) * 2020-02-05 2021-08-12 中国石油天然气集团有限公司 光纤声波传感地震数据中井筒波干扰的去除方法及装置
CN113219539B (zh) * 2020-02-05 2022-08-30 中国石油天然气集团有限公司 光纤声波传感地震数据中井筒波干扰的去除方法及装置
CN112099088A (zh) * 2020-09-16 2020-12-18 中油奥博(成都)科技有限公司 一种基于高密度光纤地震数据的油气指示及表征方法
CN112099088B (zh) * 2020-09-16 2022-04-12 中油奥博(成都)科技有限公司 一种基于高密度光纤地震数据的油气指示及表征方法

Also Published As

Publication number Publication date
CN109557587B (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
CN109557587A (zh) 一种vsp地震资料井筒波频率域滤波方法及装置
CN107272066B (zh) 一种含噪地震信号初至走时拾取方法及装置
CN108267784A (zh) 一种地震信号随机噪声压制处理方法
CN105425289B (zh) 确定低频波阻抗的方法和装置
CN105242309B (zh) 一种地震数据中规则干扰的压制方法及装置
CN102854533A (zh) 一种基于波场分离原理提高地震资料信噪比的去噪方法
CN104251934B (zh) 谐波分析方法和装置以及确定谐波间杂波的方法和装置
CN104502965A (zh) 一种振幅补偿因子的反演方法
CN104457643B (zh) 一种轨道几何检测数据的脉冲噪声滤波方法及装置
CN1839295A (zh) 流量计滤波系统和方法
CN104614767A (zh) 基于分段延拓的时变地震子波相位校正方法
CN103119848A (zh) 数据处理方法及装置
CN112183407B (zh) 一种基于时频域谱减法的隧道地震波数据去噪方法及系统
CN105068119B (zh) 低频地震数据中的面波的衰减方法及其装置
CN105093282A (zh) 基于频率约束的能量置换面波压制方法
Santillán-Guzmán et al. Real-time empirical mode decomposition for EEG signal enhancement
CN116299643A (zh) 处理核脉冲信号的数字尖峰-平顶双通道成形方法及系统
CN116136518A (zh) 色谱仪
EP3764542B1 (en) Signal enhancer and method for weak and blind signal denoising and extracting
CN107870356B (zh) 基于观测系统覆盖次数统计的真振幅补偿方法及系统
LU504061B1 (en) Frequency domain filtering method and device of wellbore wave in vsp seismic data
CN103869361B (zh) 自适应低频异常振幅压制的方法
CN109459788B (zh) 地层品质因子计算方法及系统
CN1924850A (zh) 音频快速搜索方法
CN116047604B (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