CN102692650B - 一种具有假频压制功能的井筒波分离方法 - Google Patents

一种具有假频压制功能的井筒波分离方法 Download PDF

Info

Publication number
CN102692650B
CN102692650B CN201110071041.6A CN201110071041A CN102692650B CN 102692650 B CN102692650 B CN 102692650B CN 201110071041 A CN201110071041 A CN 201110071041A CN 102692650 B CN102692650 B CN 102692650B
Authority
CN
China
Prior art keywords
wave
wave field
field
tube
frequency
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
CN201110071041.6A
Other languages
English (en)
Other versions
CN102692650A (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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201110071041.6A priority Critical patent/CN102692650B/zh
Publication of CN102692650A publication Critical patent/CN102692650A/zh
Application granted granted Critical
Publication of CN102692650B publication Critical patent/CN102692650B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明是地震勘探数据处理利用地震波视速度差异结合频率-波数滤波的具有假频压制功能的井筒波分离方法。在垂直地震剖面法资料中，利用地震波的视速度差异，将波场每道记录的每个样点时间都减去该道井筒波初至时间，变换到频率域，得到振幅谱，在道号方向变换得到频率-波数谱，反傅里叶变换得其振幅谱，再进行反傅里叶变换，得到时间域的波场，然后波场相加去除井筒波。本发明可以很好的滤除VSP资料中的井筒波，并且对有效信号伤害较小。同时使用方便简洁，运算效率高。

Description

一种具有假频压制功能的井筒波分离方法
技术领域
本发明涉及地震勘探数据处理技术，是一种利用地震波视速度差异结合频率-波数(F-K)滤波的具有假频压制功能的井筒波分离方法。
背景技术
垂直地震剖面法(VSP)资料由于特殊的采集方式和固井质量不好造成记录经常产生井筒波干扰，会大大降低资料的信噪比。
当前对井筒波的压制主要都是基于原始波场的中值滤波和频率-波数(F-K)滤波。中值滤波是一种统计滤波，其需要拾取比较准确的井筒波的初至，然后按照井筒波的初至进行中值滤波，中值滤波滤除的波场即为井筒波，在滤除井筒波的同时，改变了原始数据中的能量关系以及相位，同时滤波后的记录中仍然会有很强的残余井筒波；基于原始波场的频率-波数滤波，选择井筒波的视速度范围，该视速度范围内的波场的频率-波数谱清零，其无法对假频进行压制；因此都很难取得良好的处理效果，影响后期的处理和解释。
发明内容
本发明目的是提供一种对有效信号伤害较小，同时使用方便简洁，运算效率高的具有假频压制功能的井筒波分离方法。
本发明具体步骤包括：
1)采集记录地震原始数据；
2)在原始数据上拾取每道记录中井筒波的起跳点作为初至时间；
3)将原始数据上每道记录中井筒波初至后的信号清零，得到第一波场；
4)将原始数据上每道记录中井筒波初至前的信号清零，得到第二波场；
5)将第二波场中每道记录的每个样点时间都减去该道井筒波初至时间，得到第三波场；
6)对第三波场先进行时间方向的傅里叶变换，从而变换到频率域，得到其振幅谱；然后对振幅谱在道号方向进行傅里叶变换，从而变换到波数域，得到其频率-波数谱；
7)在步骤6)得到的频率-波数(F-K)谱中，将井筒波对应的频率-波数谱清零；然后进行波数方向的反傅里叶变换，得到其振幅谱；再对得到的振幅谱进行频率方向的反傅里叶变换，得到时间域的波场；
8)将步骤7)中得到的波场，每道记录的每个样点时间都加上该道井筒波初至时间，得到第四波场；
9)将步骤8)中得到的第四波场与步骤2)得到的第二波场相加，得到去除井筒波后的波场，完成井筒波分离。
本发明可以很好的滤除VSP资料中的井筒波，并且对有效信号伤害较小。同时使用方便简洁，运算效率高。
附图说明
图1实际VSP资料；
图2实际VSP资料的F-K谱；
图3中值滤波后的波场；
图4基于原始波场的F-K滤波后的波场；
图5基于原始波场的F-K滤波后的F-K谱；
图6本发明井筒波分离后的波场；
图7本发明井筒波分离后的波场的F-K谱。
具体实施方式
本发明是在垂直地震剖面法(VSP)资料中，利用地震波的视速度差异，对波场进行一系列的数学运算，然后利用频率-波数(F-K)变换得到其频率-波数谱，在频率-波数域进行井筒波压制，从而提高资料的信噪比。
以下结合附图详细说明本发明。
本发明具体实施方法如下：
1)采集记录地震原始数据，图1的实际资料中有明显的井筒波；
2)在原始数据上拾取每道记录中井筒波的起跳点作为初至时间；
3)将原始数据上每道记录中井筒波初至后的信号清零，得到第一波场；
4)将原始数据上每道记录中井筒波初至前的信号清零，得到第二波场；
5)将第二波场中每道记录的每个样点时间都减去该道井筒波初至时间，得到第三波场；
6)对第三波场先进行时间方向的傅里叶变换，从而变换到频率域，得到其振幅谱；然后对振幅谱在道号方向进行傅里叶变换，从而变换到波数域，得到其频率-波数谱；
7)在步骤6)得到的频率-波数(F-K)谱中，将井筒波对应的频率-波数谱清零；然后进行波数方向的反傅里叶变换，得到其振幅谱；再对得到的振幅谱进行频率方向的反傅里叶变换，得到时间域的波场；
8)将步骤7)中得到的波场，每道记录的每个样点时间都加上该道井筒波初至时间，得到第四波场；
9)将步骤8)中得到的第四波场与步骤2)得到的第二波场相加，得到去除井筒波后的波场，完成井筒波分离。
如图2中，明显可以看到由于空间采样不足引起的井筒波的空间假频。
如图3中，利用中值滤波波场分离后，仍然有部分残留的井筒波。
如图4和图5，基于原始波场的F-K滤波，无法滤除空间假频对应的井筒波。
如图6和图7，基于本发明的F-K滤波，很好的滤除了井筒波。

Claims (1)

1.一具有假频压制功能的井筒波分离方法，特点是具体步骤包括：
1)采集记录地震原始数据；
2)在原始数据上拾取每道记录中井筒波的起跳点作为初至时间；
3)将原始数据上每道记录中井筒波初至后的信号清零，得到第一波场；
4)将原始数据上每道记录中井筒波初至前的信号清零，得到第二波场；
5)将第二波场中每道记录的每个样点时间都减去该道井筒波初至时间，得到第三波场；
6)对第三波场先进行时间方向的傅里叶变换，从而变换到频率域，得到其振幅谱；然后对振幅谱在道号方向进行傅里叶变换，从而变换到波数域，得到其频率-波数谱；
7)在步骤6)得到的频率-波数(F-K)谱中，将井筒波对应的频率-波数谱清零；然后进行波数方向的反傅里叶变换，得到其振幅谱；再对得到的振幅谱进行频率方向的反傅里叶变换，得到时间域的波场；
8）将步骤7)中得到的波场，每道记录的每个样点时间都加上该道井筒波初至时间，得到第四波场；
9）将步骤8）中得到的第四波场与步骤3）得到的第一波场相加，得到去除井筒波后的波场，完成井筒波分离。
CN201110071041.6A 2011-03-23 2011-03-23 一种具有假频压制功能的井筒波分离方法 Active CN102692650B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110071041.6A CN102692650B (zh) 2011-03-23 2011-03-23 一种具有假频压制功能的井筒波分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110071041.6A CN102692650B (zh) 2011-03-23 2011-03-23 一种具有假频压制功能的井筒波分离方法

Publications (2)

Publication Number Publication Date
CN102692650A CN102692650A (zh) 2012-09-26
CN102692650B true CN102692650B (zh) 2014-10-29

Family

ID=46858228

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110071041.6A Active CN102692650B (zh) 2011-03-23 2011-03-23 一种具有假频压制功能的井筒波分离方法

Country Status (1)

Country Link
CN (1) CN102692650B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412328B (zh) * 2013-08-01 2016-04-20 中国石油天然气集团公司 基于交错网格有限差分算法的波数域保幅波场分离方法
CN104459770B (zh) * 2013-09-24 2017-06-16 中国石油化工股份有限公司 一种高维地震数据规则化方法
CN104133247B (zh) * 2014-07-31 2017-01-25 中国石油天然气集团公司 垂直地震剖面数据中套管波的压制方法及装置
CN104345344B (zh) * 2014-10-30 2017-10-17 中国石油天然气集团公司 一种微震监测数据谐波干扰压制方法及装置
CN109407144B (zh) * 2018-12-05 2020-02-18 中国矿业大学 一种基于多波的单孔孤石立体探测方法
CN109557587B (zh) * 2018-12-28 2020-10-30 长江大学 一种vsp地震资料井筒波频率域滤波方法及装置
CN111965704B (zh) * 2019-05-20 2023-09-26 中国石油天然气集团有限公司 井中地震数据的光缆次生震荡噪声压制方法及系统
CN112083491B (zh) * 2019-06-12 2023-10-31 中国石油天然气集团有限公司 基于信号属性特征的井筒波压制方法及装置
CN113219539B (zh) * 2020-02-05 2022-08-30 中国石油天然气集团有限公司 光纤声波传感地震数据中井筒波干扰的去除方法及装置
CN112099086B (zh) * 2020-09-16 2022-03-29 中油奥博(成都)科技有限公司 一种高分辨率光纤井中地震数据深频分析方法
CN114200522B (zh) * 2020-09-17 2024-04-09 中国石油化工股份有限公司 深度域地震子波提取方法和装置、存储介质及电子设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
f-k域滤波假频的消除方法;李桂元;《石油地球物理勘探》;19940630;86-92 *
李桂元.f-k域滤波假频的消除方法.《石油地球物理勘探》.1994,86-92.
柴北缘地区干扰波特征分析以及压制方法;胡杰,夏斌等;《天然气地球科学》;20101031;第21卷(第5期);851-854 *
胡杰,夏斌等.柴北缘地区干扰波特征分析以及压制方法.《天然气地球科学》.2010,第21卷(第5期),851-854.

Also Published As

Publication number Publication date
CN102692650A (zh) 2012-09-26

Similar Documents

Publication Publication Date Title
CN102692650B (zh) 一种具有假频压制功能的井筒波分离方法
Shen et al. An improved method to extract very‐broadband empirical Green’s functions from ambient seismic noise
CN102707314B (zh) 一种多路径双谱域混合相位子波反褶积方法
CN102928873B (zh) 基于四维能量聚焦的地面微地震定位方法
CN113935467B (zh) 基于迭代多尺度注意力网络的das井中勘探数据噪声抑制方法
CN103376464A (zh) 一种地层品质因子反演方法
CN109031415B (zh) 一种基于深度卷积神经网络的可控震源数据振铃压制方法
CN104216010A (zh) 利用可控震源谐波提高地震数据质量的方法
CN102053276A (zh) 一种地震数字信号的复数道集二维滤波方法
CN101930081B (zh) 地震记录单频干扰的自动识别与压制方法
CN104597502A (zh) 一种新的石油地震勘探数据去噪方法
CN104345341A (zh) 一种基于区域约束的分频段能量地震面波处理方法
CN107605470A (zh) 一种纵横波径向速度变化成像方法
CN109581481B (zh) 一种便携式高频可控震源地震信号谐波干扰消除方法
CN100349007C (zh) 消除地震记录信号中单频干扰的方法
CN105093282A (zh) 基于频率约束的能量置换面波压制方法
CN107941511B (zh) 一种基于信号时频分解的频率—峭度图的实现方法
CN102323618A (zh) 基于分数阶傅里叶变换的相干噪声抑制方法
CN102998699A (zh) 一种含有套管谐波的垂直地震剖面资料层速度反演方法
CN103675917B (zh) 一种识别微地震监测信号谐振干扰并进行有效压制的方法
CN102323616B (zh) 提高灰岩出露区地震数据分辨率的分频匹配方法
CN104570118B (zh) 一种基于双因素的自动识别与去除工业干扰的方法
CN102998701A (zh) 一种套管谐波分离方法
CN102338884B (zh) 物探中的椭圆窗方向带通保幅滤波数据处理方法
CN103344988B (zh) 基于k-l分解的可控震源信号相位检测方法

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