本发明属于花岗岩型铀矿地质勘查领域，具体涉及一种基于物化探异常的花岗岩型铀矿靶区优选方法。
背景技术：
重力、航磁和元素地球化学方法是铀矿勘查有效的、直接的三种方法。重力勘探方法在解决大地构造单元划分、基底起伏、追索构造带和了解侵入岩空间分布等地质问题中取得肯定的效果；航磁法在解决圈定侵入岩、划分岩体边界和断裂构造等方面发挥重要作用，特别是近年来寻找铀及多金属矿产方面也取得显著的效果；元素地球化学方法是通过成矿元素和伴生元素的含量异常对铀矿成矿潜力评价的一种直接方法。南方花岗岩型铀矿是我国四大类型(花岗岩型、火山岩型、砂岩型和碳硅泥岩型)铀矿中最主要类型之一，在重力、航磁和元素地球化学异常评价工作中，以往并没有进行综合整体的考虑，没有明确的基于物化探异常优选铀成矿靶区的方法流程，直接影响其评价效果，因此亟需提供一种新型的花岗岩型铀矿靶区综合优选方法。
技术实现要素：
本发明要解决的技术问题是提供一种基于物化探异常的花岗岩型铀矿靶区优选方法，该方法覆盖面广、适用性强、准确性好。
为解决上述技术问题，本发明一种基于物化探异常的花岗岩型铀矿靶区优选方法，依次包括以下步骤：
步骤一、在铀资源勘查区域内，收集重力数据、航磁数据和元素地球化学数据；
步骤二、对重力数据进行向上延拓、求取不同方向水平一阶导数处理，编制布格重力异常不同方向一阶导数图、布格重力异常向上延拓图；
步骤三、对航磁数据进行化极、向上延拓、求取一阶导数处理，编制航磁异常不同方向一阶导数图、航磁异常向上延拓图；
步骤四、对花岗岩型铀矿找矿目标区全部地球化学元素编制单元素地球化学等值图；
步骤五、对步骤二中重力数据处理图件进行断裂构造推断和划分，并绘制重力综合推断断裂构造图；
步骤六、对步骤三中航磁数据处理图件进行断裂构造推断和划分，并绘制航磁综合推断断裂构造图；
步骤七、选择主成矿元素铀、钍、钾作为主要指示元素，辅助指示元素通过与铀、钍、钾元素相关性分析进行筛选，分别计算每种元素与主要指示元素相关度系数值，选取相关度值最大的若干元素作为辅助指示元素；
步骤八、编制主要和辅助指示元素地球化学等值图，将累频值>85％标定为异常区并得到各铀、钍、钾主要元素和辅助元素的地球化学异常等值图，将所有指示元素地球化学异常等值图叠加绘制指示元素地球化学综合异常图；
步骤九、对指示元素地球化学综合异常形态进行判断，处理步骤八中的指示元素的地球化学综合异常图中，保留符合特征异常形态；
步骤十、将步骤五的重力综合推断断裂构造图、步骤六航磁综合推断断裂构造图和步骤九所得的地球化学综合异常图均与该区域的数字地质底图投影参数进行投影与叠加，最终得到找矿目的区物化探异常综合推断图；
步骤十一、在步骤十的区物化探异常综合推断图中，进行定义。
所述步骤二中向上延拓高度依次选取为5km、10km、15km和20km，水平 一阶导数的方向依次为与正北方夹角0°、45°、90°和135°。
所述步骤三中向上延拓高度依次选取为0.5km、1km、2km和25km，水平一阶导数的方向依次选取为与正北方夹角0°、45°、90°和135°。
所述步骤五中，利用布格重力异常线性梯级带、不同重力场特征的分界线、水平方向导数的线性极值带和重力异常呈线性的同向弯曲、截断、局部膨大或收敛断裂构造识别依据对步骤二中重力数据处理图件进行断裂构造推断和划分。
所述步骤七中，相关系数求取公式为：
上式中，Xi为任一测点各主要指示元素含量，为任一指示元素所有测点含量平均值，Yi为任一测点各辅助指示元素含量，为任一辅助元素所有测点含量平均值。
所述步骤九中，保留符合以下特征的异常状态：
①异常形态呈规则或不规则块状；
②同一异常带包含主要、辅助指示元素圈定范围，并且异常带内有一个或多个元素浓集中心。
所述步骤十一中，将有重力推断两组以上断裂交汇部位、航磁推断两组以上断裂交汇部位和铀、钍、钾三种主要指示元素及四种辅助指示元素中的两种或两种以上的异常带重合区且异常形态呈块状的区域判断为一级花岗岩型铀成矿远景靶区；将有重力和航磁推断一组断裂通过和铀、钍、钾三种主元素及一种辅助元素异常带重合区且异常形态呈团块状的区域判断为二级花岗岩型铀成 矿远景靶区；其他情况判断为三级花岗岩型铀成矿远景靶区。
本发明的有益技术效果在于：
一、本发明适用于全国范围内花岗岩型铀矿预测区，覆盖面广、可操作性强、有效性高。
二、本发明利用物化探资料对铀矿远景靶区进行优选分级，评价方法流程操作简明，具有高效性，保证研究结果的真实客观。
具体实施方式
下面结合实施例对本发明作进一步详细说明。
本发明基于物化探异常的花岗岩型铀矿靶区优选方法，依次包括以下步骤：
步骤一、在铀资源勘查区域内，收集1：5万～1：20万的重力数据、航磁数据和元素地球化学数据；
步骤二、对重力数据进行向上延拓、求取不同方向水平一阶导数处理，编制布格重力异常不同方向一阶导数图、布格重力异常向上延拓图；
所述步骤二中向上延拓高度依次选取为5km、10km、15km和20km，水平一阶导数的方向依次为与正北方夹角0°、45°、90°和135°；
步骤三、对航磁数据进行化极、向上延拓、求取一阶导数处理，编制航磁异常不同方向一阶导数图、航磁异常向上延拓图；
所述步骤三中化极处理需查找地球基本磁场参数，通过利用http://www.ngdc.noaa.gov/seg/geomag/jsp/IGRF.jsp网站提供的球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角；向上延拓高度依次选取为0.5km、1km、2km和25km，水平一阶导数的方向依次选取为与正北方夹角0°、45°、90°和135°；
步骤四、对花岗岩型铀矿找矿目标区全部地球化学元素编制单元素地球化 学等值图；
步骤五、利用布格重力异常线性梯级带、不同重力场特征的分界线、水平方向导数的线性极值带和重力异常呈线性的同向弯曲、截断、局部膨大或收敛断裂构造识别依据对步骤二中重力数据处理图件进行断裂构造推断和划分，并绘制重力综合推断断裂构造图，推测断裂用红色实线绘制；
步骤六、利用不同磁场区的分界线、磁异常梯度带、磁异常突变带、错动带、放射状异常带组等断裂构造识别依据对步骤三中航磁数据处理图件进行断裂构造推断和划分，并绘制航磁综合推断断裂构造图，推测断裂用蓝色虚线表示；
步骤七、选择主成矿元素铀、钍、钾作为主要指示元素，辅助指示元素通过与铀、钍、钾元素相关性分析进行筛选，分别计算每种元素与主要指示元素相关度系数值，选取相关度值最大的四种元素作为辅助指示元素；
所述步骤七中，相关系数求取公式为：
上式中，Xi为任一测点各主要指示元素含量，为任一指示元素所有测点含量平均值，Yi为任一测点各辅助指示元素含量，为任一辅助元素所有测点含量平均值；
步骤八、编制主要和辅助指示元素地球化学等值图，将累频值>85％标定为异常区并得到各铀、钍、钾主要元素和辅助元素的地球化学异常等值图，将所有指示元素地球化学异常等值图叠加绘制指示元素地球化学综合异常图；
所述步骤八中，各元素地球化学等值图异常边界分别用不同颜色绘制，指 示元素地球化学综合异常图则综合显示各元素不同颜色异常带；
步骤九、对指示元素地球化学综合异常形态进行判断，处理步骤八中的指示元素的地球化学综合异常图中，保留符合以下特征的异常：
①异常形态呈规则或不规则块状；
②同一异常带包含主要、辅助指示元素圈定范围，并且异常带内有一个或多个元素浓集中心；
步骤十、将步骤五的重力综合推断断裂构造图、步骤六航磁综合推断断裂构造图和步骤九所得的地球化学综合异常图均与该区域的数字地质底图投影参数进行投影与叠加，最终得到找矿目的区物化探异常综合推断图；
步骤十一、在步骤十的区物化探异常综合推断图中，将有重力推断两组以上断裂交汇部位、航磁推断两组以上断裂交汇部位和铀、钍、钾三种主要指示元素及四种辅助指示元素中的两种或两种以上的异常带重合区且异常形态呈块状的区域判断为一级花岗岩型铀成矿远景靶区；将有重力和航磁推断一组断裂通过和铀、钍、钾三种主元素及一种辅助元素异常带重合区且异常形态呈团块状的区域判断为二级花岗岩型铀成矿远景靶区；其他情况判断为三级花岗岩型铀成矿远景靶区；
一级远景靶区是最有利的花岗岩型铀成矿靶区，二级远景靶区是有利的花岗岩型铀成矿靶区，三级远景靶区是较有利的花岗岩型铀成矿靶区，可按一、二、三级远景靶区依次开展工作。