[发明专利]一种火山岩型铀矿床磁性界面反演方法

权利要求书

1.一种火山岩型铀矿床磁性界面反演方法，其特征在于：具体包括以下步骤：

步骤一、在勘查区开展高精度磁法测量，获得磁测数据；

步骤二、在勘查区开展RTK测量，获得勘查区的测点精确经纬度坐标和高程数据；

步骤三、收集勘查区地层、岩体磁化率数据并把勘查区地层、岩体磁化率数据转换成CGS国际单位制；

步骤四、对步骤一获得的磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正，得到测点异常△T值；

步骤五、对步骤二得到的测点高程数据进行网格化处理，对步骤四所得的磁异常△T值进行数据平滑处理，得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理；

步骤六、对步骤五所得到的磁异常△T网格数据进行化极处理，化极处理所需的磁倾角、磁偏角由步骤三所得出；

步骤七、开始创建反演初始模型，并把模型分为两层，根据步骤三中不同岩体磁化率数据对两层模型进行磁化率赋值，针对反演初始模型的两层模型的交界面创建反演目标界面常量网格；

步骤八、对目标界面进行反演，在反演过程中加入余弦滤波器，余弦滤波器选择低通滤波，且设置上限截止波长，下限截止波长，反演出磁性界面网格；

步骤九、修改余弦滤波器参数上限截止波长与下限截止波长；

步骤十、重复步骤八、步骤九，迭代次数为M次，直至拟合误差达到5nT，获得最终磁性界面网格。

2.根据权利要求1所述的一种火山岩型铀矿床磁性界面反演方法，其特征在于：所述的步骤一中，磁测总精度应高于2nT，勘查区选在地形平坦的区域。

3.根据权利要求2所述的一种火山岩型铀矿床磁性界面反演方法，其特征在于：所述的步骤五、步骤六中，所述的网格化处理中的差值方法为最小曲率法。

4.根据权利要求3所述的一种火山岩型铀矿床磁性界面反演方法，其特征在于：所述的步骤七中，反演目标界面常量网格中，网格常量为-1000m、所需的反演观测面选择步骤五得出的高程网格。

5.根据权利要求4所述的一种火山岩型铀矿床磁性界面反演方法，其特征在于：所述的步骤八中，对目标界面进行反演采用帕克反演方法采用下式进行计算：

式中为磁性界面起伏的傅里叶变换，ΔT为磁异常，J是上下界面磁化率差值，S为径向波数，n为迭代次数。

6.根据权利要求5所述的一种火山岩型铀矿床磁性界面反演方法，其特征在于：所述的步骤九中，对步骤八得到的磁性界面网格利用帕克正演公式进行正演计算，所述的帕克正演方法采用下式进行计算：

式中为磁性界面起伏的傅里叶变换，为反演出磁性界面起伏深度，J是上下界面磁化率差值，S为径向波数，n为迭代次数。

7.根据权利要求6所述的一种火山岩型铀矿床磁性界面反演方法，其特征在于：所述的步骤十中，M≥10。