| 主权项 |
1.一种地球物理勘探中的矩形窗带通保幅滤波数据处理方法，特点是通过如下步骤实现：1)在工区利用重力仪或磁力仪在野外采集重力或磁力资料；2)对野外采集的重磁资料进行各种常规校正及改正；3)采用插值方法将重磁资料网格化，形成规则网重磁资料f<sub>a</sub>(x,y)，其中x、y分别为重磁资料的平面网格点坐标；4）利用低精度的重磁资料补充镶边：将规则网重磁资料f<sub>a</sub>（x,y）转变为其线数和点数扩大为2的整数幂重磁资料f<sub>b</sub>（x,y）；5）利用快速傅氏正变换（FFT）将镶边后的重磁资料f<sub>b</sub>(x,y）变换到波数域，获得f<sub>b</sub>(x,y）的频谱F<sub>b</sub>(u,v），F<sub>b</sub>(u,v）＝FFT［f<sub>b</sub>(x,y）］，其中u、v分别为x和y方向的波数；6）用下列公式计算两个初始的滤波器：<img file="FDA00002835632800011.GIF" wi="508" he="50" /><img file="FDA00002835632800012.GIF" wi="503" he="50" />其中：β为正有理数，u<sub>0</sub>、v<sub>0</sub>分别为x和y方向的中心波数；7）用下列公式分别对两个初始的滤波器进行补偿：<img file="FDA00002835632800013.GIF" wi="556" he="53" /><img file="FDA00002835632800014.GIF" wi="565" he="53" />其中：n、m分别为u、v方向的补偿次数，为正整数；8）在波数域，计算矩形窗带通滤波频率响应函数：<img file="FDA00002835632800015.GIF" wi="466" he="55" />9）在波数域，将频谱F<sub>b</sub>（u,v）与带通滤波频率响应函数<img file="FDA00002835632800021.GIF" wi="466" he="56" />相乘：<img file="FDA00002835632800022.GIF" wi="713" he="67" />10）经过快速傅氏反变换（IFFT），获得带通滤波后的重磁资料f<sub>c</sub>（x，y），f<sub>c</sub>(x，y）＝IFFT［F<sub>c</sub>(u,v)］11）反镶边处理，即从f<sub>c</sub>（x，y）中提取原数据大小的带通滤波后的重磁数据f<sub>d</sub>（x，y）。 |
