[发明专利]一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法

权利要求书

1.一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法，特点是通过如下步骤实现：

1)在工区利用重力仪或磁力仪在野外采集重力或磁力资料；

2)对野外采集的重磁资料进行各种常规校正及改正；

3)采用插值方法将重磁资料网格化，形成规则网重磁资料f_a(x，y)，其中x、y分别为重磁资料的平面网格点坐标；

4)利用低精度的重磁资料补充镶边；将规则网重磁资料f_a(x，y)转变为其线数和点数扩大为2的整数幂重磁资料f_b(x，y)；

5)利用快速傅氏正变换(FFT)将镶边后的重磁资料f_b(x，y)变换到波数域，获得f_b(x，y)的频谱F_b(u，v)，F_b(u，v)＝FFT[f_b(x，y)]，其中u、v分别为x和y方向的波数；

6)用下列公式计算初始的滤波器：

其中：β为正有理数，u₀、v₀分别为x和y方向的中心波数；

7)用下列公式对初始的滤波器进行补偿：

其中：n为补偿次数，为正整数；

步骤7)中所述的补偿次数n，1≤n≤30。

8)在波数域，将频谱F_b(u，v)与带通滤波频率响应函数相乘，即：

9)经过快速傅氏反变换(IFFT)，获得带通滤波后的重磁资料f_c(x，y)，f_c(x，y)＝IFFT[F_c(u，v)]

10)反镶边处理，即从f_c(x，y)中提取原数据大小的带通滤波后的重磁数据f_d(x，y)。

2.根据权利要求1所述的方法，特点是步骤7)中所述的补偿次数n，1≤n≤30。