专利名称：一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法
技术领域：
本发明涉及地球物理勘探技术领域，具体是一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法。
背景技术：
地球物理勘探包括重力、磁力、电法、地震勘探等，重力、磁力是发展最早、理论成熟、应用广泛的两种地球物理方法。重力、磁力方法覆盖面积广、经济快速，广泛应用于油气、固体矿产资源勘探、区域与深部构造研究、环境与工程勘察等各个领域。重力、磁力勘探包括数据采集、处理与解释三部分。数字滤波是在数字信号情况下对信号作滤波处理，在信号数字处理中，低通、高通、带通滤波是最常用的最基本的滤波，可以在频率域和时间域设计这些滤波器。频域滤波是在频域上选定某一范围，使这一范围内的信号通过或去除，可通过FFT 变换及FFT反变换对信号进行处理。时域滤波是在时域上根据所需要的信号用褶积滤波或递归的方法对信号进行数字处理。理想的带通滤波器使某一频率范围内的所有信号无损通过，而让该频率范围以外的信号被完全阻止。理想的带通滤波器有以下两个特征其一，带宽在通带以内的信号可以无失真地通过滤波器；其二，滤波器的通带和阻带是截然分开的。其实，这种理想滤波器在实际中并不存在，它们只是一种分析问题的模型，实际应用中的滤波器都只是理想滤波器的某种近似。理想的频率响应是不可能实现的，由于它的波数响应是突变的，因此，滤波器产生的吉布斯现象不能有效地压制干扰，突出有用信号， 而且还可能使有效信号的频谱产生畸变。实际应用中，选用一个连续函数来模拟带通滤波器可以有效消除吉布斯现象。一般情况下，连续函数在通带和阻带之间有一部分频率衰减的区域，这部分区域的信号只是受到了部分衰减。理想的带通滤波器设计的要求就是这个连续函数应该有较陡的边界，从而使衰减区域尽可能的�。员愀寐瞬ㄆ髂茏畲笙薅冉咏昝劳ù纳杓�。根据局部异常与区域异常的波谱差异，Dean (1958), Byerly (1965), Fraser 等(1966)根据重磁异常的波谱特征来设计数字滤波器。Darby和Davies (1967)， Fuller (1967)，以及^irfIueh(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，ν)，Fb(u，ν) = FFT[fb(x，y)]，其中U、ν分别为χ和y方向的波数；6)用下列公式计算初始的滤波器
約)(m，v) = &χρ[~2πβ (u-ua)2 +(V-V0)2)1'2]其中β为正有理数，Uc^vtl分别为χ和y方向的中心波数；7)用下列公式对初始的滤波器进行补偿
Ψη ( , ν) = [2 - φη_λ (u, ν)]φη_λ (u, ν)其中η为补偿次数，为正整数；步骤7)中所述的补偿次数n，1 < η < 30。8)在波数域，将频谱Fb (U，ν)与带通滤波频率响应函数代(Μ, V) = [2-PdO, V)]队—(Μ, V)相乘，即
Fo (u, v) =Fb(u, ν) X φη(μ,ν)9)经过快速傅氏反变换(IFFT)，获得带通滤波后的重磁资料f。(x，y)，f。(x，y)= IFFT[Fc(u，ν)]10)反镶边处理，即从f。(x，y)中提取原数据大小的带通滤波后的重磁数据fd(x， y)。本发明滤波器在通带内信号保幅，在阻带内信号完全抑制，该滤波器的通带与阻带之间的衰减区域非常�。芄谎怪浦卮乓斐５推当尘昂透咂翟肷怀鲋衅涤杏眯藕牛乇鹗屎贤怀鲈诓ㄊ蛭残未暗闹衅涤杏眯藕�。本发明滤波器通带内的信号具有保幅特征，由于滤波器是个连续函数，因而，不存在由于滤波器的突变产生的吉布斯现象。本发明不仅可以用于重磁资料的数据处理，也可用于电法、地震等地球物理资料的处理。
图1是ID圆形窗波谱特征。随着补偿次数η的增加，滤波器的频带逐渐变宽。图2是2D圆形窗波谱特征。图3是使用本发明某工区滤波后剩余重力异常图。图4是使用本发明某工区3D地震Tg构造图。
具体实施例方式以下结合附图详细说明本发明。利用LCR型重力仪在某工区采集1/5万高精度重力资料。按照发明内容所述的步骤对重力资料进行处理，即1)在工区利用重力仪在野外采集重力资料；2)对野外采集的重力资料进行各种常规校正及改正；3)利用插值方法将重力资料网格化，形成规则网重力资料fa(x，y)，其中x、y分别为重力资料的平面网格点坐标；网格化数据大小为275行X327列。4)利用低精度的重力资料补充镶边；将规则网重力资料fa(x，y)转变为其线数和点数扩大为2的整数幂重力资料fb(x，y)；为了达到最好的镶边效果，利用该区其它重力资料补充镶边，将规则网重力资料fa(x，y)转变为其线数和点数扩大为2的整数幂重力资料 fb(x,y)；镶边后网格化数据大小为512行X512列；5)利用快速傅氏正变换(FFT)将镶边后的重力资料fb(x，y)变换到波数域，获得 fb(x, y)的频谱Fb (u，ν)，Fb(u，ν) = FFT[fb(x，y)]，其中U、ν分别为χ和y方向的波数；6)用下列公式计算初始的滤波器
炉0Ο,ν) = βχρ[-2πβ^(α - U0)2 + (ν - νη)2 )1/2]其中β为正有理数，Uc^vtl分别为χ和y方向的中心波数；7)用下列公式对初始的滤波器进行补偿
Ψη (u,v) = [2- φη_λ (u, ν)]φη_, (u, ν)其中η为补偿次数，为正整数；步骤7)中所述的补偿次数n，1 < η < 30。图1是一维圆形窗波谱特征。随着补偿次数的增加，滤波器的通带逐渐变宽。图2是二维圆形窗带通保幅滤波器的波谱特征，从图中可以看到，滤波器的滤波区域为一个圆形，滤波器是一个带通滤波器，通带与阻带之间的衰减区域非常小，因而滤波器具有保幅特征。
8)在波数域，将频谱Fb (U，ν)与带通滤波频率响应函数队(KV) = P-^v1O,V^^VJM，V)相乘，即
Fc(u, V) =Fb(11，V) X 队(Μ,V)9)经过快速傅氏反变换(IFFT)，获得带通滤波后的重力资料f。(x，y)，f。(x，y)= IFFT[Fc(u，ν)]10)反镶边处理，即从f。(x，y)中提取原数据大小的带通滤波后的重力数据fd(x， y)。fd(x, y)数据大小为275行X 327列；图3为某工区利用本方法获得的剩余重力异常图，图4为某工区三维地震Tg构造图。图3中发现的局部重力高被后来的地震资料所证实，验证了重力资料的成果。上述离散数据网格化、镶边处理、反镶边处理的过程，本专业技术领域分析人员能熟练完成。本发明所涉及的技术也可用于其它地球物理资料的数据处理以及信号数字处理领域。
权利要求
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分别为χ和y方向的波数；6)用下列公式计算初始的滤波器
2.根据权利要求1所述的方法，特点是步骤7)中所述的补偿次数n，1< η < 30。
全文摘要
本发明是地球物理勘探中的圆形窗带通保幅滤波数据处理方法，采用插值方法将重磁资料网格化，利用低精度的重磁资料补充镶边，快速傅氏正变换，在波数域计算初始的圆形带窗滤波频率响应函数并进行补偿，经过快速傅氏反变换，获得带通滤波后的重磁资料，反镶边处理提取原数据大小的带通滤波后的重磁数据。本发明滤波器在通带内信号保幅，在阻带内信号完全抑制，能够压制重磁异常低频背景和高频噪声，特别适合突出在波数域为圆形窗的中频有用信号。
文档编号G01V3/38GK102338890SQ20101051631
公开日2012年2月1日 申请日期2010年10月22日 优先权日2010年10月22日
发明者张研, 文百红, 杨辉, 王玲 申请人:中国石油天然气股份有限公司
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
