一种地球物理勘探中的重磁数据处理方法

摘要

一种地球物理勘探中的重磁数据处理方法，其包括步骤：在波数域将重磁资料频谱分别与任意方向的n阶导数或积分频率响应函数、延拓回返响应函数及差分法下延算子频率响应函数相乘，得到下延后的频谱；再经过快速傅氏反变换，获得向下延拓后的三维重磁数据体。利用向下延拓之后的三维重磁数据体结合地震、地质、钻井资料进行对比、分析并提交下一步勘探的有利区带及目标。该方法不仅有更大的向下延拓深度并且有较好的抗噪能力，从而获得稳定的下半空间三维重磁资料分布图像。

主权项

1、一种地球物理勘探中的重磁数据处理方法，该方法包括以下步骤：(1)首先利用重力仪或磁力仪在野外采集重磁资料，或者在室内对原有的重磁异常图进行数字化而得到重磁资料，然后对野外采集或室内数字化的重磁资料进行各种常规校正及改正；(2)利用Surfer中的Kriging插值方法将重磁资料网格化，形成规则网重磁资料f_a(x，y)，其中x、y分别为重磁资料的平面网格点坐标；(3)进行镶边处理：将实际重磁资料f_a(x，y)的线数和点数扩大为2的整数幂的重磁资料f_b(x，y)；利用低精度的重磁资料进行补充镶边；(4)利用快速傅氏正变换(FFT)将f_b(x，y)变换到波数域，获得实测重磁资料的频谱F_b(u，v)，F_b(u，v)＝FFT[f_b(x，y)]，其中u、v分别为x和y方向的波数；(5)在波数域，将镶边后的重磁资料频谱F_b(u，v)与重磁场任意方向的n阶导数或积分频率响应函数$q^n(u,v)={[2\mathrm{\pi i}(\mathrm{\alpha u}+\mathrm{\beta v})+\gamma \sqrt{u^2+v^2}]}^n$相乘，增强局部或区域重磁异常，即：F_c(u，v)＝F_b(u，v)×qⁿ(u，v)其中n为有理数，i为虚数；α、β、γ分别为所求导数或积分方向的三个方向余弦；n为导数或积分的阶数，即分辨率参数；(6)在波数域，根据向下延拓深度h，将增强后的局部或区域重磁异常F_c(u，v)与延拓回返频域响应函数相乘，抑制高频成分，使其满足向下延拓的条件，即：其中h为下延深度，p为延拓回返次数；h，p为正有理数；(7)在波数域，在抑制高频成分的基础上，与差分法向下延拓频域响应函数$d(u,v)=(6-\exp (-2\mathrm{\pi h}\sqrt{u^2+v^2}-2\cos \left(2\mathrm{\pi hu}\right)-2\cos \left(2\mathrm{\pi hv}\right)))$相乘，得到下延后的频谱F_e(u，v)F_e(u，v)＝F_d(u，v)×d(u，v)；(8)再经过快速傅氏反变换(IFFT)，获得向下延拓后的重磁资料f_e(x，y)，f_e(x，y)＝IFFT[F_e(u，v)]；(9)进行反镶边处理，即从f_e(x，y)中提取原数据大小的向下延拓深度h后的重磁数据f_f(x，y)；(10)重复步骤(6)，从而获得下半空间三维重磁数据体f_g(x，y，z)，其中，x，y，z为重磁异常的空间网格点坐标；(11)利用向下延拓之后的三维重磁数据体结合地震、地质、钻井资料进行对比、分析并提交下一步勘探的有利区带及目标。