[发明专利]一种重磁数据的空间外插方法

说明书

技术领域

本发明涉及一种重磁数据的空间外插方法，特别是涉及一种能够反映实际重磁数据的频谱分形特征的重磁数据的空间外插方法。

背景技术

在频率域中利用快速付立叶变换方法（FFT）进行对面积性重磁数据进行处理和转换时必须要求数据为网格数据，且大小为2的整数次方。空间域中的大部分重磁数据处理方法往往会造成实测数据的边部在处理后会损失。由于实测数据范围有限，为了尽量减少空间域内数据处理的边部损失和满足频率域中数据处理的要求，更充分地利用实测资料进行地质解释，往往需对实测数据进行空间外插。

常见的外插方法有补零外插方法、对折外插方法、区域场外插方法和余弦外插方法。上述方法均使得扩边结果与实际数据之间产生一定的畸变，从而产生一定的Gibbs效应，并导致重磁数据处理和转换结果的精度降低。因此，必须寻找既能反映重磁数据基本特征，又能有效减少边界效应的方法。

因此亟需提供一种新型的重磁数据的空间外插方法。

发明内容

本发明要解决的技术问题是提供一种能够反映实际重磁数据的频谱分形特征的重磁数据的空间外插方法。

为解决上述技术问题，本发明一种重磁数据的空间外插方法，依次包括以下步骤：

（1）采集位场数据PA，包括重磁数据的平均功率谱密度S、波数f；得到分维值D；

所述分维值D由式（1）获得

Ln(S)＝c+(2D-8)Ln(f)     式（1）

式（1）中，c为常数；

（2）确定空间外插的边长K，并生成该范围内的高斯白噪声数据；

（3）对高斯白噪声数据进行二维付立叶变换，得到一组复系数H_mn；

Hmn=(1K)2Σj=0K-1Σk=0K-1hjkexp[-2πiK(mj+nk)];]]>

（4）得到一组新的复系数

Hmn*=Hmn/f(7-2D)/2;]]>

（5）对由复系数组成的矩阵进行二维付立叶反变换，得到插值范围内的数据PB；

（6）采用最小曲率法对数据PA范围内的PB数据进行外插，得到数据PC；

（7）将数据PC从数据PB中减去，得到数据PA范围内为零值的数据PD；

（8）采用最小曲率法对位场数据PA进行外插，并与数据PD相加，就得到了空间外插的结果数据PE。