CN102692650A - 一种具有假频压制功能的井筒波分离方法 - Google Patents
一种具有假频压制功能的井筒波分离方法 Download PDFInfo
- Publication number
- CN102692650A CN102692650A CN2011100710416A CN201110071041A CN102692650A CN 102692650 A CN102692650 A CN 102692650A CN 2011100710416 A CN2011100710416 A CN 2011100710416A CN 201110071041 A CN201110071041 A CN 201110071041A CN 102692650 A CN102692650 A CN 102692650A
- Authority
- CN
- China
- Prior art keywords
- wave
- wave field
- obtains
- tube
- field
- 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
Images
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)中得到的第四波场与步骤2)得到的第二波场相加,得到去除井筒波后的波场,完成井筒波分离。
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 true CN102692650A (zh) | 2012-09-26 |
CN102692650B 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) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412328A (zh) * | 2013-08-01 | 2013-11-27 | 中国石油天然气集团公司 | 基于交错网格有限差分算法的波数域保幅波场分离方法 |
CN104133247A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 垂直地震剖面数据中套管波的压制方法及装置 |
CN104345344A (zh) * | 2014-10-30 | 2015-02-11 | 中国石油天然气集团公司 | 一种微震监测数据谐波干扰压制方法及装置 |
CN104459770B (zh) * | 2013-09-24 | 2017-06-16 | 中国石油化工股份有限公司 | 一种高维地震数据规则化方法 |
CN109407144A (zh) * | 2018-12-05 | 2019-03-01 | 中国矿业大学 | 一种基于多波的单孔孤石立体探测方法 |
CN109557587A (zh) * | 2018-12-28 | 2019-04-02 | 长江大学 | 一种vsp地震资料井筒波频率域滤波方法及装置 |
CN111965704A (zh) * | 2019-05-20 | 2020-11-20 | 中国石油天然气集团有限公司 | 井中地震数据的光缆次生震荡噪声压制方法及*** |
CN112083491A (zh) * | 2019-06-12 | 2020-12-15 | 中国石油天然气集团有限公司 | 基于信号属性特征的井筒波压制方法及装置 |
CN112099086A (zh) * | 2020-09-16 | 2020-12-18 | 中油奥博(成都)科技有限公司 | 一种高分辨率光纤井中地震数据深频分析方法 |
CN113219539A (zh) * | 2020-02-05 | 2021-08-06 | 中国石油天然气集团有限公司 | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 |
CN114200522A (zh) * | 2020-09-17 | 2022-03-18 | 中国石油化工股份有限公司 | 深度域地震子波提取方法和装置、存储介质及电子设备 |
-
2011
- 2011-03-23 CN CN201110071041.6A patent/CN102692650B/zh active Active
Non-Patent Citations (2)
Title |
---|
李桂元: "f-k域滤波假频的消除方法", 《石油地球物理勘探》 * |
胡杰,夏斌等: "柴北缘地区干扰波特征分析以及压制方法", 《天然气地球科学》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412328B (zh) * | 2013-08-01 | 2016-04-20 | 中国石油天然气集团公司 | 基于交错网格有限差分算法的波数域保幅波场分离方法 |
CN103412328A (zh) * | 2013-08-01 | 2013-11-27 | 中国石油天然气集团公司 | 基于交错网格有限差分算法的波数域保幅波场分离方法 |
CN104459770B (zh) * | 2013-09-24 | 2017-06-16 | 中国石油化工股份有限公司 | 一种高维地震数据规则化方法 |
CN104133247A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 垂直地震剖面数据中套管波的压制方法及装置 |
CN104345344A (zh) * | 2014-10-30 | 2015-02-11 | 中国石油天然气集团公司 | 一种微震监测数据谐波干扰压制方法及装置 |
CN109407144A (zh) * | 2018-12-05 | 2019-03-01 | 中国矿业大学 | 一种基于多波的单孔孤石立体探测方法 |
CN109557587A (zh) * | 2018-12-28 | 2019-04-02 | 长江大学 | 一种vsp地震资料井筒波频率域滤波方法及装置 |
CN111965704B (zh) * | 2019-05-20 | 2023-09-26 | 中国石油天然气集团有限公司 | 井中地震数据的光缆次生震荡噪声压制方法及*** |
CN111965704A (zh) * | 2019-05-20 | 2020-11-20 | 中国石油天然气集团有限公司 | 井中地震数据的光缆次生震荡噪声压制方法及*** |
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 | 中国石油天然气集团有限公司 | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 |
CN112099086B (zh) * | 2020-09-16 | 2022-03-29 | 中油奥博(成都)科技有限公司 | 一种高分辨率光纤井中地震数据深频分析方法 |
CN112099086A (zh) * | 2020-09-16 | 2020-12-18 | 中油奥博(成都)科技有限公司 | 一种高分辨率光纤井中地震数据深频分析方法 |
CN114200522A (zh) * | 2020-09-17 | 2022-03-18 | 中国石油化工股份有限公司 | 深度域地震子波提取方法和装置、存储介质及电子设备 |
CN114200522B (zh) * | 2020-09-17 | 2024-04-09 | 中国石油化工股份有限公司 | 深度域地震子波提取方法和装置、存储介质及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102692650B (zh) | 2014-10-29 |
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 | |
CN103376464A (zh) | 一种地层品质因子反演方法 | |
CN101598812B (zh) | 去除数字检波器单点接收地震记录中的异常噪声方法 | |
CN113935467B (zh) | 基于迭代多尺度注意力网络的das井中勘探数据噪声抑制方法 | |
CN109031415B (zh) | 一种基于深度卷积神经网络的可控震源数据振铃压制方法 | |
CN104216010A (zh) | 利用可控震源谐波提高地震数据质量的方法 | |
CN102193108A (zh) | 一种提高石油勘探资料处理信噪比的方法 | |
CN104345341A (zh) | 一种基于区域约束的分频段能量地震面波处理方法 | |
CN105242309A (zh) | 一种地震数据中规则干扰的压制方法及装置 | |
CN107605470A (zh) | 一种纵横波径向速度变化成像方法 | |
CN101930081A (zh) | 地震记录单频干扰的自动识别与压制方法 | |
CN109581481B (zh) | 一种便携式高频可控震源地震信号谐波干扰消除方法 | |
CN100349007C (zh) | 消除地震记录信号中单频干扰的方法 | |
CN105093282A (zh) | 基于频率约束的能量置换面波压制方法 | |
CN102323618A (zh) | 基于分数阶傅里叶变换的相干噪声抑制方法 | |
CN103675917B (zh) | 一种识别微地震监测信号谐振干扰并进行有效压制的方法 | |
CN102998699A (zh) | 一种含有套管谐波的垂直地震剖面资料层速度反演方法 | |
CN102338884B (zh) | 物探中的椭圆窗方向带通保幅滤波数据处理方法 | |
CN102998701A (zh) | 一种套管谐波分离方法 | |
CN110737022B (zh) | 一种可控震源激发地震资料黑三角区噪音的压制方法 | |
CN104570118B (zh) | 一种基于双因素的自动识别与去除工业干扰的方法 | |
CN102338890B (zh) | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 | |
CN101907726B (zh) | 一种自动识别和消除地震勘探工业电干扰的方法 | |
CN112099086B (zh) | 一种高分辨率光纤井中地震数据深频分析方法 |
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 |