一种磁异常化极方法

摘要

本发明公开了一种磁异常化极方法。包括以下步骤：对采集到的磁场数据进行校正，得到磁异常数据；将磁异常数据形成规则的磁异常数据t(x,y)，记为矩阵T(1:m,1:n)；获得研究区的地磁场方向，并整理成与矩阵T(1:m,1:n)等大的矩阵I(1:m,1:n)；对矩阵T(1:m,1:n)进行镶边处理，转变成新矩阵T’(1‑m’:m+m’,1‑n’:n+n’)；设计频率域化极因子，并做反傅里叶变换，得到化极结果；计算化极结果的垂直梯度和总梯度，并以最大互相关系数确定剩磁方向；根据地磁场方向和剩磁方向，并对矩阵T’和I在频率域内进行变倾角化极处理，计算研究区磁异常的变倾角化极结果。通过上述方式，本发明能够提高磁法勘探的准确度。

主权项

一种磁异常化极方法，其特征在于，包括以下步骤：采集研究区的磁场数据，对采集到的所述磁场数据进行校正，得到磁异常数据；将所述磁异常数据网格化，形成规则的磁异常数据t(x,y)，其中x,y分别表示磁异常数据的平面网格点坐标，记为矩阵T(1:m,1:n)，大小为m×n，m表示行数，n表示点数；获得所述研究区的地磁场方向，并将所述地磁场方向整理成与矩阵T(1:m,1:n)等大的矩阵，记为I(1:m,1:n)；对所述矩阵T(1:m,1:n)进行镶边处理，将矩阵T(1:m,1:n)转变成行数和列数均为2的最小整数幂的新矩阵T’(1‑m’:m+m’,1‑n’:n+n’)；设计频率域化极因子，并对所述频率域化极因子做反傅里叶变换，得到化极结果；计算化极结果的垂直梯度和总梯度，并计算其互相关系数，以最大互相关系数确定剩磁方向；根据所述地磁场方向和所述剩磁方向，并根据所述频率域化极因子对所述矩阵T’和I在频率域内进行变倾角化极处理，计算所述研究区磁异常的变倾角化极结果。