一种热液型铀矿床地球物理场源边界识别方法技术

本发明专利技术属于热液型铀矿床磁力数据处理领域，具体涉及一种热液型铀矿床地球物理场源边界识别方法。本发明专利技术包括以下步骤：步骤1、地表磁力异常网格化处理、RTK测量所有测点精确坐标、高程信息；步骤2、针对磁力异常网格进行化极处理；步骤3、对磁异常化极网格求取总水平方向导数；步骤4、利用汉宁平滑滤波器对步骤3得出的总水平方向导数网格进行滤波；步骤5、将步骤4网格中所有节点与其周边值进行对比；步骤6、将步骤5中数据库中的所有的峰值投影到布骤五所得到的磁异常△T化极网格图上，场源边界则显示在图中。本发明专利技术能够有效、真实地探测地下铀源岩体或地球物理场源体边界位置形态和范围。

【技术实现步骤摘要】
一种热液型铀矿床地球物理场源边界识别方法
本专利技术属于热液型铀矿床磁力数据处理领域，具体涉及一种热液型铀矿床地球物理场源边界识别方法。
技术介绍
热液型铀矿床是指含矿热液在一定的物理化学条件下，在各种有利的构造和岩体中，由充填、交代及沉积等成矿方式形成的有用矿物堆积体。包括花岗岩型和火山岩型铀矿床，它是我国主要工业类型。铀矿产出部位一般在岩体的内部或外接触带上。因此，铀源岩体或地球物理场源体边界位置形态、范围大小、平面展布特征等信息显得尤为重要，岩体边界及外接触带或重大断裂对岩体内部、外围成矿作用的控制具有重要的找矿意义。由于受到斜磁化、地下场源埋深、规模、磁异常叠加以及噪声等多种因素的影响，基于地表磁异常直接进行场源边缘的检测和识别是目前工作的难点。而且在热液型铀矿床地球物理勘查评价工作中，以往并没有明确的方法或没有明确的基于磁力异常进行地球物理场源边界识别的方法流程，直接影响热液型铀矿床调查评价效果，因此亟需提供一种新型的识别地球物理场源边界的方法。
技术实现思路
本专利技术要解决的技术问题：本专利技术提供一种热液型铀矿床地球物理场源边界识别方法，能够有效、真实地探测地下铀源岩体或地球物理场源体边界位置形态和范围。本专利技术采用的技术方案：一种热液型铀矿床地球物理场源边界识别方法，依次包括以下步骤：步骤1、地表磁力异常网格化处理、RTK测量所有测点精确坐标、高程信息；步骤2、针对磁力异常网格进行化极处理；步骤3、对磁异常化极网格求取总水平方向导数；步骤4、利用汉宁平滑滤波器对步骤3得出的总水平方向导数网格进行滤波，压制低振幅的浅部信息，汉宁滤波器的级数设定为0；步骤5、将步骤4网格中所有节点与其周边值进行对比，峰值检测级别选定为3级，选取出峰值，再将网格中所有峰值提取出来放入数据库；步骤6、将步骤5中数据库中的所有的峰值投影到布骤五所得到的磁异常△T化极网格图上，场源边界则显示在图中。所述步骤1包括如下步骤：步骤1.1、在勘查区开展地面高精度磁法测量，获得磁测数据，为了达到更好的效果，磁测总精度应高于2nT，测区选在地形平坦的区域；步骤1.2、在勘查区开展RTK测量，获得测点精确经纬度坐标和高程数据；步骤1.3、对步骤1.1获得的磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正，得到测点异常△T值，在对磁异常数据进行日变校正时，需查找地球基本磁场参数，通过球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角；步骤1.4、对步骤1.3所得的磁异常△T值进行数据平滑处理，消除高频干扰，得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理，同样形成401×402的磁异常△T网格文件，网格化差值方法优选最小曲率法。所述步骤2的具体步骤为：对布骤1.4所得到的磁异常△T网格数据进行化极处理，化极处理所需的磁倾角、磁偏角由步骤1.3所得出，形成401×402的磁异常化极网格文件，网格化差值方法优选最小曲率法。所述步骤3的具体步骤为：对布骤2所得到的磁异常△T化极网格数据求取总水平方向导数，求取总水平方向导数方法采用式(1)进行计算：式中THD表示计算点的磁异常总水平方向导数，ΔT为磁异常。所述步骤5中，周边值指的是：X方向，Y方向以及两条对角值；3级指的是：节点值比周围3个方向值高。所述步骤6中，绘制地球物理场源边界，给峰值标上色标，绘制的符号通过长轴方向将边缘走向显示出来，这些符号将追踪地质接触带边界。当将峰值投影到步骤3得出的总水平导数图中时，这些符号将追踪局部变化。本专利技术的有益技术效果在于：利用本专利技术技术方案后，对热液型铀矿床地球物理场源边界进行了特征识别，有效地评价了地下地球物理场源体的分布规律；从而达到快速、精确地评价热液型铀矿资源，对地质找矿工作具有重要的实际意义。附图说明图1为本专利技术所提供的一种热液型铀矿床地球物理场源边界识别方法流程图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。如图1所示，本专利技术提供的一种热液型铀矿床地球物理场源边界识别方法，依次包括以下步骤：步骤1、地表磁力异常网格化处理、RTK测量所有测点精确坐标、高程信息；步骤1.1、在勘查区开展地面高精度磁法测量，获得磁测数据，为了达到更好的效果，磁测总精度应高于2nT，测区选在地形平坦的区域；步骤1.2、在勘查区开展RTK测量(实时动态定位技术)，获得测点精确经纬度坐标和高程数据；步骤1.3、对步骤1.1获得的磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正，得到测点异常△T值，在对磁异常数据进行日变校正时，需查找地球基本磁场参数，通过球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角；步骤1.4、对步骤1.3所得的磁异常△T值进行数据平滑处理，消除高频干扰，得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理，同样形成401×402的磁异常△T网格文件，网格化差值方法优选最小曲率法；步骤2、针对磁力异常网格进行化极处理；对布骤1.4所得到的磁异常△T网格数据进行化极处理，化极处理所需的磁倾角、磁偏角由步骤1.3所得出，形成401×402的磁异常化极网格文件，网格化差值方法优选最小曲率法；步骤3、对磁异常化极网格求取总水平方向导数；对布骤2所得到的磁异常△T化极网格数据求取总水平方向导数，求取总水平方向导数方法采用式(1)进行计算：式中THD表示计算点的磁异常总水平方向导数，ΔT为磁异常。步骤4、利用汉宁平滑滤波器对步骤3得出的总水平方向导数网格进行滤波，压制低振幅的浅部信息，汉宁滤波器的级数设定为0；步骤5、将步骤4网格中所有节点与其周边值进行对比(X方向，Y方向以及两条对角值)，峰值检测级别选定为3级(即节点值比周围3个方向值高)，选取出峰值，再将网格中所有峰值提取出来放入数据库；步骤6、将步骤5中数据库中的所有的峰值投影到布骤五所得到的磁异常△T化极网格图上，场源边界则显示在图中。绘制地球物理场源边界，给峰值标上色标，绘制的符号通过长轴方向将边缘(地质接触关系)走向显示出来，这些符号将追踪地质接触带边界。当将峰值投影到步骤3得出的总水平导数图中时，这些符号将追踪局部变化。本文档来自技高网...

【技术保护点】
1.一种热液型铀矿床地球物理场源边界识别方法，其特征在于：依次包括以下步骤：步骤(1)、地表磁力异常网格化处理、RTK测量所有测点精确坐标、高程信息；步骤(2)、针对磁力异常网格进行化极处理；步骤(3)、对磁异常化极网格求取总水平方向导数；步骤(4)、利用汉宁平滑滤波器对步骤(3)得出的总水平方向导数网格进行滤波，压制低振幅的浅部信息，汉宁滤波器的级数设定为0；步骤(5)、将步骤(4)网格中所有节点与其周边值进行对比，峰值检测级别选定为3级，选取出峰值，再将网格中所有峰值提取出来放入数据库；步骤(6)、将步骤(5)中数据库中的所有的峰值投影到布骤五所得到的磁异常△T化极网格图上，场源边界则显示在图中。

【技术特征摘要】
1.一种热液型铀矿床地球物理场源边界识别方法，其特征在于：依次包括以下步骤：步骤(1)、地表磁力异常网格化处理、RTK测量所有测点精确坐标、高程信息；步骤(2)、针对磁力异常网格进行化极处理；步骤(3)、对磁异常化极网格求取总水平方向导数；步骤(4)、利用汉宁平滑滤波器对步骤(3)得出的总水平方向导数网格进行滤波，压制低振幅的浅部信息，汉宁滤波器的级数设定为0；步骤(5)、将步骤(4)网格中所有节点与其周边值进行对比，峰值检测级别选定为3级，选取出峰值，再将网格中所有峰值提取出来放入数据库；步骤(6)、将步骤(5)中数据库中的所有的峰值投影到布骤五所得到的磁异常△T化极网格图上，场源边界则显示在图中。2.根据权利要求1所述的热液型铀矿床地球物理场源边界识别方法，其特征在于：所述步骤(1)包括如下步骤：步骤(1.1)、在勘查区开展地面高精度磁法测量，获得磁测数据，为了达到更好的效果，磁测总精度应高于2nT，测区选在地形平坦的区域；步骤(1.2)、在勘查区开展RTK测量，获得测点精确经纬度坐标和高程数据；步骤(1.3)、对步骤(1.1)获得的磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正，得到测点异常△T值，在对磁异常数据进行日变校正时，需查找地球基本磁场参数，通过球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角；步骤(1.4)、对步骤(1.3)所得的磁异常△T值进行数据平滑处理，消除高...

【专利技术属性】
技术研发人员：喻翔，杨龙泉，赵丹，吴国东，陈聪，周俊杰，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11