一种地球物理勘探中的矩形窗带通保幅滤波数据处理方法

摘要

本发明是物探技术中矩形窗带通保幅滤波数据处理方法，采用插值方法将重磁资料网格化，利用低精度的重磁资料补充镶边，快速傅氏正变换，在波数域计算两个初始带通滤波频率响应函数，对两个初始的滤波器进行补偿并相乘获得矩形窗带通滤波器，经过快速傅氏反变换，获得带通滤波后的重磁资料，反镶边处理提取原数据大小的带通滤波后的重磁数据。本发明滤波器在通带内信号保幅，在阻带内信号完全抑制，能够压制重磁异常低频背景和高频噪声，特别适合突出在波数域为矩形窗的中频有用信号。

主权项

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、m分别为u、v方向的补偿次数，为正整数；8)在波数域，计算矩形窗带通滤波频率响应函数：9)在波数域，将频谱F_b(u，v)与带通滤波频率响应函数相乘：10)经过快速傅氏反变换(IFFT)，获得带通滤波后的重磁资料f_c(x，y)，f_c(x，y)＝IFFT[F_c(u，v)]11)反镶边处理，即从f_c(x，y)中提取原数据大小的带通滤波后的重磁数据f_d(x，y)。