CN102338890B - 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 - Google Patents
一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 Download PDFInfo
- Publication number
- CN102338890B CN102338890B CN 201010516311 CN201010516311A CN102338890B CN 102338890 B CN102338890 B CN 102338890B CN 201010516311 CN201010516311 CN 201010516311 CN 201010516311 A CN201010516311 A CN 201010516311A CN 102338890 B CN102338890 B CN 102338890B
- Authority
- CN
- China
- Prior art keywords
- magnetic data
- data
- gravity
- frequency
- filtering
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是地球物理勘探中的圆形窗带通保幅滤波数据处理方法,采用插值方法将重磁资料网格化,利用低精度的重磁资料补充镶边,快速傅氏正变换,在波数域计算初始的圆形带窗滤波频率响应函数并进行补偿,经过快速傅氏反变换,获得带通滤波后的重磁资料,反镶边处理提取原数据大小的带通滤波后的重磁数据。本发明滤波器在通带内信号保幅,在阻带内信号完全抑制,能够压制重磁异常低频背景和高频噪声,特别适合突出在波数域为圆形窗的中频有用信号。
Description
技术领域
本发明涉及地球物理勘探技术领域,具体是一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法。
背景技术
地球物理勘探包括重力、磁力、电法、地震勘探等,重力、磁力是发展最早、理论成熟、应用广泛的两种地球物理方法。重力、磁力方法覆盖面积广、经济快速,广泛应用于油气、固体矿产资源勘探、区域与深部构造研究、环境与工程勘察等各个领域。重力、磁力勘探包括数据采集、处理与解释三部分。
数字滤波是在数字信号情况下对信号作滤波处理,在信号数字处理中,低通、高通、带通滤波是最常用的最基本的滤波,可以在频率域和时间域设计这些滤波器。
频域滤波是在频域上选定某一范围,使这一范围内的信号通过或去除,可通过FFT变换及FFT反变换对信号进行处理。时域滤波是在时域上根据所需要的信号用褶积滤波或递归的方法对信号进行数字处理。
理想的带通滤波器使某一频率范围内的所有信号无损通过,而让该频率范围以外的信号被完全阻止。理想的带通滤波器有以下两个特征:其一,带宽在通带以内的信号可以无失真地通过滤波器;其二,滤波器的通带和阻带是截然分开的。
其实,这种理想滤波器在实际中并不存在,它们只是一种分析问题的模型,实际应用中的滤波器都只是理想滤波器的某种近似。理想的频率响应是不可能实现的,由于它的波数响应是突变的,因此,滤波器产生的吉布斯现象不能有效地压制干扰,突出有用信号,而且还可能使有效信号的频谱产生畸变。
实际应用中,选用一个连续函数来模拟带通滤波器可以有效消除吉布斯现象。一般情况下,连续函数在通带和阻带之间有一部分频率衰减的区域,这部分区域的信号只是受到了部分衰减。理想的带通滤波器设计的要求就是这个连续函数应该有较陡的边界,从而使衰减区域尽可能的小,以便该滤波器能最大限度接近完美通带的设计。
根据局部异常与区域异常的波谱差异,Dean(1958),Byerly(1965),Fraser等(1966)根据重磁异常的波谱特征来设计数字滤波器。Darby和Davies(1967),Fuller(1967),以及Zurflueh(1967)设计二维滤波器来对网格化的重磁数据提取剩余重磁异常。
在波数域中,对于压制信号中的高频成分来说,现有的充分光滑的函数,虽然能起很大的压制作用,但同时也会在一定程度上压制中频和低频成分,这种对中频和低频成分的压制,却是我们所不希望发生的。为了减小这种对中、低频成分的压制,使有用信号的畸变尽可能降到最小,侯重初(1979,1981)(侯重初,一种压制干扰的频率滤波方法,物探与化探,1979,NO.5;侯重初,补偿圆滑滤波方法,石油物探,1981,NO.2。)提出了一种称为补偿圆滑的低通滤波方法,该方法是针对解释重力和磁力异常时常常需要压制高频干扰而研究的。
补偿圆滑滤波是针对压制重磁高频干扰而提出的一个低通滤波器(功能有限),而数字处理中更重要是突出某一频段的信号,亦即需要一个有不同滤波窗口的带通滤波器。除了补偿圆滑滤波之外的其它滤波器的主要问题有三个,一是滤波器不是一个连续函数,也即通带与阻带之间突然间断,从而造成吉布斯现象;二是不适合突出在波数域为圆形窗的中频有用信号;三是滤波器的通带与阻带之间的衰减区域比较大,通带与阻带没有截然分开,滤波器在通带内信号不保幅,在阻带内信号没有完全抑制,不仅影响后续的定性分析,更为严重的是影响定量分析。目前,地球物理资料保幅的处理已越来越重要。
发明内容
本发明目的是提供一种在通带内信号保幅,在阻带内信号完全抑制的地球物理勘探中的圆形窗带通保幅滤波数据处理方法。
本发明通过如下步骤实现:
1)在工区利用重力仪或磁力仪在野外采集重力或磁力资料;
2)对野外采集的重磁资料进行各种常规校正及改正;
3)采用插值方法将重磁资料网格化,形成规则网重磁资料fa(x,y),其中x、y分别为重磁资料的平面网格点坐标;
4)利用低精度的重磁资料补充镶边;将规则网重磁资料fa(x,y)转变为其线数和点数扩大为2的整数幂重磁资料fb(x,y);
5)利用快速傅氏正变换(FFT)将镶边后的重磁资料fb(x,y)变换到波数域,获得fb(x,y)的频谱Fb(u,v),Fb(u,v)=FFT[fb(x,y)],其中u、v分别为x和y方向的波数;
6)用下列公式计算初始的滤波器:
其中:β为正有理数,u0、v0分别为x和y方向的中心波数;
7)用下列公式对初始的滤波器进行补偿:
其中:n为补偿次数,为正整数;
步骤7)中所述的补偿次数n,1≤n≤30。
9)经过快速傅氏反变换(IFFT),获得带通滤波后的重磁资料fc(x,y),fc(x,y)=IFFT[Fc(u,v)]
10)反镶边处理,即从fc(x,y)中提取原数据大小的带通滤波后的重磁数据fd(x,y)。
本发明滤波器在通带内信号保幅,在阻带内信号完全抑制,该滤波器的通带与阻带之间的衰减区域非常小,能够压制重磁异常低频背景和高频噪声,突出中频有用信号,特别适合突出在波数域为圆形窗的中频有用信号。
本发明滤波器通带内的信号具有保幅特征,由于滤波器是个连续函数,因而,不存在由于滤波器的突变产生的吉布斯现象。
本发明不仅可以用于重磁资料的数据处理,也可用于电法、地震等地球物理资料的处理。
附图说明
图1是1D圆形窗波谱特征。随着补偿次数n的增加,滤波器的频带逐渐变宽。
图2是2D圆形窗波谱特征。
图3是使用本发明某工区滤波后剩余重力异常图。
图4是使用本发明某工区3D地震Tg构造图。
具体实施方式
以下结合附图详细说明本发明。
利用LCR型重力仪在某工区采集1/5万高精度重力资料。按照发明内容所述的步骤对重力资料进行处理,即:
1)在工区利用重力仪在野外采集重力资料;
2)对野外采集的重力资料进行各种常规校正及改正;
3)利用插值方法将重力资料网格化,形成规则网重力资料fa(x,y),其中x、y分别为重力资料的平面网格点坐标;网格化数据大小为275行×327列。
4)利用低精度的重力资料补充镶边;将规则网重力资料fa(x,y)转变为其线数和点数扩大为2的整数幂重力资料fb(x,y);为了达到最好的镶边效果,利用该区其它重力资料补充镶边,将规则网重力资料fa(x,y)转变为其线数和点数扩大为2的整数幂重力资料fb(x,y);镶边后网格化数据大小为512行×512列;
5)利用快速傅氏正变换(FFT)将镶边后的重力资料fb(x,y)变换到波数域,获得fb(x,y)的频谱Fb(u,v),Fb(u,v)=FFT[fb(x,y)],其中u、v分别为x和y方向的波数;
6)用下列公式计算初始的滤波器:
其中:β为正有理数,u0、v0分别为x和y方向的中心波数;
7)用下列公式对初始的滤波器进行补偿:
其中:n为补偿次数,为正整数;
步骤7)中所述的补偿次数n,1≤n≤30。
图1是一维圆形窗波谱特征。随着补偿次数的增加,滤波器的通带逐渐变宽。
图2是二维圆形窗带通保幅滤波器的波谱特征,从图中可以看到,滤波器的滤波区域为一个圆形,滤波器是一个带通滤波器,通带与阻带之间的衰减区域非常小,因而滤波器具有保幅特征。
9)经过快速傅氏反变换(IFFT),获得带通滤波后的重力资料fc(x,y),fc(x,y)=IFFT[Fc(u,v)]
10)反镶边处理,即从fc(x,y)中提取原数据大小的带通滤波后的重力数据fd(x,y)。fd(x,y)数据大小为275行×327列;
图3为某工区利用本方法获得的剩余重力异常图,图4为某工区三维地震Tg构造图。图3中发现的局部重力高被后来的地震资料所证实,验证了重力资料的成果。
上述离散数据网格化、镶边处理、反镶边处理的过程,本专业技术领域分析人员能熟练完成。
本发明所涉及的技术也可用于其它地球物理资料的数据处理以及信号数字处理领域。
Claims (1)
1.一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法,特点是通过如下步骤实现:
1)在工区利用重力仪或磁力仪在野外采集重磁资料;
2)对野外采集的重磁资料进行各种常规校正及改正;
3)采用插值方法将重磁资料网格化,形成规则网重磁资料fa(x,y),其中x、y分别为重磁资料的平面网格点坐标;
4)利用低精度的重磁资料补充镶边;将规则网重磁资料fa(x,y)转变为其线数和点数扩大为2的整数幂重磁资料fb(x,y);
5)利用快速傅氏正变换(FFT)将镶边后的重磁资料fb(x,y)变换到波数域,获得fb(x,y)的频谱Fb(u,v),Fb(u,v)=FFT[fb(x,y)],其中u、v分别为x和y方向的波数;
6)用下列公式计算初始的滤波器:
其中:β为正有理数,u0、v0分别为x和y方向的中心波数;
7)用下列公式对初始的滤波器进行补偿:
其中:n为补偿次数,为正整数;
步骤7)中所述的补偿次数n,1≤n≤30;
9)经过快速傅氏反变换(IFFT),获得带通滤波后的重磁资料fc(x,y),fc(x,y)=IFFT[Fc(u,v)]
10)反镶边处理,即从fc(x,y)中提取原数据大小的带通滤波后的重磁数据fd(x,y)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010516311 CN102338890B (zh) | 2010-10-22 | 2010-10-22 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010516311 CN102338890B (zh) | 2010-10-22 | 2010-10-22 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102338890A CN102338890A (zh) | 2012-02-01 |
CN102338890B true CN102338890B (zh) | 2013-04-24 |
Family
ID=45514717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010516311 Active CN102338890B (zh) | 2010-10-22 | 2010-10-22 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102338890B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103513288B (zh) * | 2012-06-21 | 2016-06-08 | 中国石油天然气集团公司 | 一种二维网格数据的补偿方向滤波方法 |
CN104618033B (zh) * | 2015-01-08 | 2015-09-09 | 东南大学 | 一种多层自适应形态滤波重力信号噪声抑制方法 |
CN104932010B (zh) * | 2015-06-09 | 2018-08-07 | 中国海洋石油集团有限公司 | 一种基于近道镶边稀疏Radon变换的绕射波分离方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797041A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种采用深度域滤波消除线性与非线性干扰波的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100399052C (zh) * | 2001-12-18 | 2008-07-02 | Bhp比利顿创新公司 | 处理海洋磁场梯度数据的方法和使用该数据的勘探方法 |
US7054750B2 (en) * | 2004-03-04 | 2006-05-30 | Halliburton Energy Services, Inc. | Method and system to model, measure, recalibrate, and optimize control of the drilling of a borehole |
-
2010
- 2010-10-22 CN CN 201010516311 patent/CN102338890B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797041A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种采用深度域滤波消除线性与非线性干扰波的方法 |
Non-Patent Citations (4)
Title |
---|
应用计算机"推导"位场圆滑滤波及二次导数公式的方法;杨辉;《物探化探计算技术》;19981231;第20卷(第04期);347-35页 * |
张凤旭等.波数域重力归一化总梯度法中的圆滑滤波因子.《物探与化探》.2005,第29卷(第03期),248-252页. |
杨辉.应用计算机"推导"位场圆滑滤波及二次导数公式的方法.《物探化探计算技术》.1998,第20卷(第04期),347-35页. |
波数域重力归一化总梯度法中的圆滑滤波因子;张凤旭等;《物探与化探》;20050630;第29卷(第03期);248-252页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102338890A (zh) | 2012-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105242309B (zh) | 一种地震数据中规则干扰的压制方法及装置 | |
CN105652322B (zh) | 多分量地震数据的t-f-k域极化滤波方法 | |
CN102692650A (zh) | 一种具有假频压制功能的井筒波分离方法 | |
CN110208856B (zh) | 一种基于流形分区2d-vmd的沙漠复杂噪声压制方法 | |
CN104345341A (zh) | 一种基于区域约束的分频段能量地震面波处理方法 | |
CN102338890B (zh) | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 | |
CN104597502A (zh) | 一种新的石油地震勘探数据去噪方法 | |
CN104635264B (zh) | 叠前地震数据的处理方法及设备 | |
CN103630936A (zh) | 基于波束定向原理的地震单炮记录随机噪声压制方法 | |
CN105445801A (zh) | 一种消除二维地震资料随机噪音的处理方法 | |
CN102338884B (zh) | 物探中的椭圆窗方向带通保幅滤波数据处理方法 | |
CN111708087A (zh) | 一种基于DnCNN神经网络对地震数据噪声压制的方法 | |
CN102323618B (zh) | 基于分数阶傅里叶变换的相干噪声抑制方法 | |
CN105093282B (zh) | 基于频率约束的能量置换面波压制方法 | |
CN102103215B (zh) | 三维高密度地震勘探记录叠前面波压制方法 | |
CN102368094B (zh) | 一种地球物理勘探中的矩形窗带通保幅滤波数据处理方法 | |
GUO et al. | The design and experiment of IIR lowpass digital filters for airborne gravity data | |
Simaan et al. | Optimum suppression of coherent signals with linear moveout in seismic data | |
CN105676292A (zh) | 一种基于二维曲波变换的三维地震数据去噪方法 | |
CN111257938A (zh) | 基于小波互相关时移地震虚拟震源波场重构方法和系统 | |
CN103513288B (zh) | 一种二维网格数据的补偿方向滤波方法 | |
CN107561578B (zh) | 一种针对宽方位地震资料的面波压制方法 | |
CN110716172B (zh) | 基于频率挑选的矢量水听器包络谱估计方法 | |
Zhao et al. | Two-dimensional nonstationary convolutional filtering and an adaptive linear interference suppression method | |
CN113093282A (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 |