一种重磁数据的空间外插方法

摘要

本发明涉及一种重磁数据的空间外插方法，依次包括：1.采集位场数据，得到分维值；2.确定空间外插的边长，并生成该范围内的高斯白噪声数据；3.对高斯白噪声数据进行二维付立叶变换，得到一组复系数；4.得到一组新的复系数；5.得到插值范围内的数据PB；6.采用最小曲率法对数据PA范围内的PB数据进行外插，得到数据PC；7.将数据PC从数据PB中减去，得到数据PA范围内为零值的数据PD；8.采用最小曲率法对位场数据PA进行外插，并与数据PD相加，就得到了空间外插的结果数据PE。本发明能够反映实际重磁数据的频谱分形特征，与实际数据之间连续且光滑，可以有效减少频率域处理中的Gibbs效应。

主权项

1.一种重磁数据的空间外插方法，依次包括以下步骤：（1）采集位场数据PA，包括重磁数据的平均功率谱密度S、波数f；得到分维值D；所述分维值D由式（1）获得Ln(S)＝c+(2D-8)Ln(f)    式（1）式（1）中，c为常数；（2）确定空间外插的边长K，并生成该范围内的高斯白噪声数据；（3）对高斯白噪声数据进行二维付立叶变换，得到一组复系数H_mn；$H_{\mathrm{mn}}={\left(\frac{1}{K}\right)}^2\underset{j=0}{\overset{K-1}{\Sigma }}\underset{k=0}{\overset{K-1}{\Sigma }}{h}_{\mathrm{jk}}\exp [-\frac{2\mathrm{\pi i}}{K}(\mathrm{mj}+\mathrm{nk})];$（4）得到一组新的复系数$H_{\mathrm{mn}}^{*}=H_{\mathrm{mn}}/f^{(7-2D)/2};$（5）对由复系数组成的矩阵进行二维付立叶反变换，得到插值范围内的数据PB；（6）采用最小曲率法对数据PA范围内的PB数据进行外插，得到数据PC；（7）将数据PC从数据PB中减去，得到数据PA范围内为零值的数据PD；（8）采用最小曲率法对位场数据PA进行外插，并与数据PD相加，就得到了空间外插的结果数据PE。