本申请公开了一种基于多重K均值聚类分析的找矿预测方法及系统，用以解决现有的找矿方法存在误差且效率较低的技术问题，属于找矿预测技术领域。方法包括：在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围；基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区；收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。本申请提出的方法针对赋煤区展开研究，在中生代赋煤断陷盆地中，概选出存在于含煤地层及煤系矿产高变质岩中的1:50万航磁低缓异常区，快速圈定富铁矿远景区；通过对物探资料进行聚类分析，识别次级异常，提高找矿预测精度。
1.一种基于多重K均值聚类分析的找矿预测方法，其特征在于，所述方法包括：步骤S1：在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围；确定有利区的过程，具体为：步骤S11：基于赋煤盆地内煤系变质及航磁异常信息，确定远景区；
步骤S12：在所述远景区对重磁数据分别进行统计筛选，得到数据值异常区域；
步骤S13：将ΔG垂向一阶导数与ΔT化极异常分别归一化后进行相关性分析，确定重磁同高异常范围，从而得到所述有利区；其中，所述ΔG垂向一阶导数为重力异常沿垂直方向求得的一阶导数，所述ΔT化极异常是与岩石层中的磁性物质相关的异常磁场；
对所述磁场总量异常ΔT数据进行位场转换，以得到磁场参数；其中，所述磁场参数包括：ΔT 化极异常、ΔT 化极垂向一阶导数、ΔT 化极垂向二阶导数、ΔT化极总梯度模、ΔT 化极剩余异常、ΔT 化极标准差、ΔT 化极水平梯度模以及ΔT 化极下延不同深度数据；
基于K均值聚类方法依次对所述磁场参数进行识别，得到各磁场参数对于矿体的反映，继而筛选出第一地质单元；其中，所述第一地质单元为成矿概率超过预设值的单元；
步骤S2：基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区；
步骤S3：收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
2.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述方法还包括：采集重力值，并计算布格重力异常ΔG数据；
对所述布格重力异常ΔG数据进行位场转换，以得到重力参数；其中，所述重力参数包括：剩余重力异常、ΔG垂向一阶导数、ΔG总梯度模、ΔG标准差、ΔG水平梯度模、ΔG下延不同深度数据；
基于K均值聚类方法依次对所述重力参数进行识别，得到第二地质单元。
3.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述步骤S12，具体包括：开展高精度重力测量、磁法测量，两相邻测线距离为100 250m，两相邻测点距离为两相~邻测线距离的1/5 1/4。
4.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述步骤S3中，确定找矿靶区的过程，具体为：步骤S31：对所述重点区内地质、钻探资料以及物探结果进行多元综合分析；
步骤S33：结合矽卡岩型富铁矿成矿因素，预测重要找矿标志。
5.根据权利要求4所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述重要找矿标志包括：大理岩化碳酸盐岩、角岩化碎屑岩、蚀变中基性岩、磁铁矿化蚀变、矽卡岩化蚀变、钠长石化蚀变、蛇纹石化蚀变、构造交汇部位、层间滑脱带、脆性裂隙。
6.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述找矿靶区岩体全铁含量较周边低3%以上、镁含量较周边低2%以上，钾含量较周边低1%以上，钠含量较周边高1.5%以上；围岩地层钙含量较周边低1%以上，铁含量较周边低1%以上。
7.一种基于多重K均值聚类分析的找矿预测系统，应用如权利要求1‑6任意一项所述的方法，其特征在于，所述系统包括：有利区确定模块，用于在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围；确定有利区的过程，具体为：基于赋煤盆地内煤系变质及航磁异常信息，确定远景区；
在所述远景区对重磁数据分别进行统计筛选，得到数据值异常区域；
将ΔG垂向一阶导数与ΔT化极异常分别归一化后进行相关性分析，确定重磁同高异常范围，从而得到所述有利区；其中，所述ΔG垂向一阶导数为重力异常沿垂直方向求得的一阶导数，所述ΔT化极异常是与岩石层中的磁性物质相关的异常磁场；
对所述磁场总量异常ΔT数据进行位场转换，以得到磁场参数；其中，所述磁场参数包括：ΔT 化极异常、ΔT 化极垂向一阶导数、ΔT 化极垂向二阶导数、ΔT化极总梯度模、ΔT 化极剩余异常、ΔT 化极标准差、ΔT 化极水平梯度模以及ΔT 化极下延不同深度数据；
基于K均值聚类方法依次对所述磁场参数进行识别，得到各磁场参数对于矿体的反映，继而筛选出第一地质单元；其中，所述第一地质单元为成矿概率超过预设值的单元；
聚类模块，用于基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区；
预测模块，用于收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
一种基于多重K均值聚类分析的找矿预测方法及系统
技术领域
[0001] 本申请涉及找矿预测技术领域，尤其涉及一种基于多重K均值聚类分析的找矿预测方法及系统。
背景技术
[0002] 我国富铁矿以矽卡岩型铁矿为主，以往受理论和方法限制，以查证高磁异常和就矿找矿为主，一直未取得明显的找矿突破；随着找矿工作向深部推进，深部信息的获取难度越来越大，深部微弱异常信号难以获取和识别，加之地球物理探测的多解性，严重影响了找矿的效率和精度，亟需突破传统找矿思路，创新地球物理探测技术，为富铁矿深部找矿提供支撑。
[0003] 卡岩型铁矿一般赋存于侵入岩体与灰岩接触带内，我国大量中生代赋煤盆地内煤系之下发育稳定了灰岩地层，同时也发现局部侵入岩对煤层产生了一定的影响，灰岩与侵入岩体为矽卡岩型富铁矿成矿奠定了物质基础，在含煤盆地内分布的低缓磁异常成为寻找矽卡岩型富铁矿的重要新方向。但受上覆地层屏蔽影响，深部目标地质体异常信息微弱，相关研究中无系统性针对矽卡岩型富铁矿找矿预测方法；单一地球物理方法难以满足深部探测需求，同时使用多种参数进行识别异常时，存在人工拾取具体位置存在误差且效率较低等问题，亟需提高预测精度。
发明内容
[0004] 本申请实施例提供了一种基于多重K均值聚类分析的找矿预测方法及系统，用以解决现有的找矿方法存在误差且效率较低的技术问题。
[0005] 一方面，本申请实施例提供了一种基于多重K均值聚类分析的找矿预测方法，所述方法包括：步骤S1：在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围。步骤S2：基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区。 步骤S3：收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
[0006] 在本申请的一种实现方式中，确定有利区的过程，具体为：步骤S11：基于赋煤盆地内煤系变质及航磁异常信息，确定远景区。步骤S12：在所述远景区对重磁数据分别进行统计筛选，得到数据值异常区域。步骤S13：将ΔG垂向一阶导数与ΔT化极异常分别归一化后进行相关性分析，确定重磁同高异常范围，从而得到所述有利区；其中，所述ΔG垂向一阶导数为重力异常沿垂直方向求得的一阶导数，所述ΔT化极异常是与岩石层中的磁性物质相关的异常磁场。
[0007] 在本申请的一种实现方式中，所述步骤S13之后，所述方法还包括：采集磁场值，并计算磁场总量异常ΔT数据。对所述磁场总量异常ΔT数据进行位场转换，以得到磁场参数；其中，所述磁场参数包括：ΔT 化极异常、ΔT 化极垂向一阶导数、ΔT 化极垂向二阶导数、ΔT化极总梯度模、ΔT 化极剩余异常、ΔT 化极标准差、ΔT 化极水平梯度模以及ΔT 化极下延不同深度数据。基于K均值聚类方法依次对所述磁场参数进行识别，得到各磁场参数对于矿体的反映，继而筛选出第一地质单元；其中，所述第一地质单元为成矿概率超过预设值的单元。
[0008] 在本申请的一种实现方式中，所述方法还包括：采集重力值，并计算布格重力异常ΔG数据。对所述布格重力异常ΔG数据进行位场转换，以得到重力参数；其中，所述重力参数包括：剩余重力异常、ΔG垂向一阶导数、ΔG总梯度模、ΔG标准差、ΔG水平梯度模、ΔG下延不同深度数据。基于K均值聚类方法依次对所述重力参数进行识别，得到第二地质单元。
[0009] 在本申请的一种实现方式中，所述步骤S12，具体包括：开展高精度重力测量、磁法测量，两相邻测线距离为100 250m，两相邻测点距离为两相邻测线距离的1/5 1/4。~ ~
[0010] 在本申请的一种实现方式中，所述步骤S2中，确定重点区的过程，具体包括：步骤S21：确定K个聚类中心。
[0011] 步骤S22：对磁场参数及重力参数设定不同权重系数。
[0012] 步骤S23：随机选择某一数据点的位置作为单元中心点，计算所有数据点与该中心的加权距离和Li，同时确定与该中心点最远距离的数据点位置为第2个中心点，计算公式如下：
[0013] 。
[0014] 。
[0015] 。
[0016] 其中， 为重力或磁法数据点参数值， 为数据点参数值中的最小值，l(m ,n)为单个数据点与聚类中心的加权距离，Li为所有数据点与聚类中心的加权距离之和， 为数据点参数值中的最大值， 为数据点坐标， 为第i个聚类中
心的坐标，m为数据种类，n为数据量， 为参数权重系数， 为第i个聚类中心参数值，a、b分别对应实际计算中选择的磁法与重力的权重系数。
[0017] 步骤S24：重复步骤S23，直至K个加权距离和达到最小值，确定K个聚类中心，并将两组数据分配到其距离最近的聚类中心，使得重力与磁法的处理结果被分为K组，确定多个重点区。
[0018] 在本申请的一种实现方式中，所述步骤S3中，确定找矿靶区的过程，具体为：步骤S31：对所述重点区内地质、钻探资料以及物探结果进行多元综合分析。
[0019] 步骤S32：推断由浅至深地层及构造分布情况。
[0020] 步骤S33：结合矽卡岩型富铁矿成矿因素，预测重要找矿标志。
[0021] 在本申请的一种实现方式中，所述重要找矿标志包括：大理岩化碳酸盐岩、角岩化碎屑岩、蚀变中基性岩、磁铁矿化蚀变、矽卡岩化蚀变、钠长石化蚀变、蛇纹石化蚀变、构造交汇部位、层间滑脱带、脆性裂隙。
[0022] 在本申请的一种实现方式中，所述找矿靶区岩体全铁含量较周边低3%以上、镁含量较周边低2%以上，钾含量较周边低1%以上，钠含量较周边高1.5%以上；围岩地层钙含量较周边低1%以上，铁含量较周边低1%以上。
[0023] 另一方面，本申请还提供了一种基于多重K均值聚类分析的找矿预测系统，应用如前述的方法，所述系统包括：有利区确定模块，用于在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围；聚类模块，用于基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区；预测模块，用于收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
[0024] 本申请实施例提供的一种基于多重K均值聚类分析的找矿预测方法及系统，具有以下有益效果：针对赋煤区展开研究，在中生代赋煤断陷盆地中，概选出存在于含煤地层及煤系矿产高变质岩中的1:50万航磁低缓异常区，快速圈定富铁矿远景区；通过对物探资料进行聚类分析，识别次级异常，提高找矿预测精度。
附图说明
[0025] 此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
[0026] 图1为本申请实施例提供的一种基于多重K均值聚类分析的找矿预测方法流程图；
[0027] 图2为本申请实施例提供的一种基于多重K均值聚类分析的找矿预测系统组成图。
具体实施方式
[0028] 为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
[0029] 本申请实施例提供了一种基于多重K均值聚类分析的找矿预测方法及系统，用以解决现有的找矿方法存在误差且效率较低的技术问题。
[0030] 下面通过附图对本申请实施例提出的技术方案进行详细的说明。
[0031] 图1为本申请实施例提供的一种基于多重K均值聚类分析的找矿预测方法流程图。如图1所示，该方法主要包括以下步骤。
[0032] 步骤S1：在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围。
[0033] 步骤S2：基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区。
[0034] 步骤S3：收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
[0035] 在本申请实施例中，确定有利区的过程，具体为：步骤S11：基于赋煤盆地内煤系变质及航磁异常信息，确定远景区。
[0036] 首先，收集分析已有工作区地质、钻探及物探资料，重点分析已有航磁测量资料。在中生代赋煤断陷盆地中，圈定存在于含煤地层及煤系矿产高变质岩中的1:50万航磁低缓正异常区，即与背景值相比高100 300nT的区域，初步概选远景区。其中含煤地层的变质特~
征主要包括角岩化、红柱石化等，煤系矿产变质特征包括煤的煤质程度较周边高2个煤阶以上或直接变为天然焦。
[0037] 步骤S12：在所述远景区对重磁数据分别进行统计筛选，得到数据值异常区域。
[0038] 步骤S13：将ΔG垂向一阶导数与ΔT化极异常分别归一化后进行相关性分析，确定重磁同高异常范围，从而得到所述有利区；其中，所述ΔG垂向一阶导数为重力异常沿垂直方向求得的一阶导数，所述ΔT化极异常是与岩石层中的磁性物质相关的异常磁场。
[0039] 在本申请实施例中，所述步骤S13之后，所述方法还包括：采集磁场值，并计算磁场总量异常ΔT数据；对所述磁场总量异常ΔT数据进行位场转换，以得到磁场参数；其中，所述磁场参数包括：ΔT 化极异常、ΔT 化极垂向一阶导数、ΔT 化极垂向二阶导数、ΔT化极总梯度模、ΔT 化极剩余异常、ΔT 化极标准差、ΔT 化极水平梯度模以及ΔT 化极下延不同深度数据；基于K均值聚类方法依次对所述磁场参数进行识别，得到各磁场参数对于矿体的反映，继而筛选出第一地质单元；其中，所述第一地质单元为成矿概率超过预设值的单元。包括：ΔT 异常、ΔT 化极垂向一阶导数异常靠正磁异常中部的等值线密集区或梯级带部位;ΔT 化极异常、ΔT 化极垂向二阶导数异常的中间部位或正异常中心偏边部位置;ΔT 化极剩余磁异常的中心位置;ΔT 化极标准差的两个高异常值之间中心位置；ΔT 化极水平梯度的正异常中心偏边部位置；ΔT 化极下延不同深度异常的中间部位。
[0040] 在本申请实施例中，所述方法还包括：采集重力值，并计算布格重力异常ΔG数据；对所述布格重力异常ΔG数据进行位场转换，以得到重力参数；其中，所述重力参数包括：剩余重力异常、ΔG垂向一阶导数、ΔG总梯度模、ΔG标准差、ΔG水平梯度模、ΔG下延不同深度数据。
[0041] 基于K均值聚类方法依次对所述重力参数进行识别，得到第二地质单元。包括：剩余重力异常、ΔG垂向一阶导数等值线密集区或梯级带位置； ΔG总梯度模、ΔG标准差及ΔG水平梯度模的正异常中心偏边部位置； ΔG下延不同深度的正异常中心偏边部等值线密集区或梯级带部位。
[0042] 在本申请实施例中，所述步骤S12，具体包括：开展高精度重力测量、磁法测量，两相邻测线距离为100 250m，两相邻测点距离为两相邻测线距离的1/5 1/4。在实际探测中，~ ~测线与测点距离可以根据现场地形情况及探测需求进行调整；异常中心区域加密探测，测点数量为其他区域的2 3倍。
~
[0043] 在本申请实施例中，所述步骤S2中，确定重点区的过程，具体包括：步骤S21：确定K个聚类中心。
[0044] 步骤S22：对磁场参数及重力参数设定不同权重系数。
[0045]
[0046] 表1权重系数
[0047] 步骤S23：随机选择某一数据点的位置作为单元中心点，计算所有数据点与该中心的加权距离和Li，同时确定与该中心点最远距离的数据点位置为第2个中心点，计算公式如下。
[0048] 。
[0049] 。
[0050] 。
[0051] 其中， 为重力或磁法数据点参数值， 为数据点参数值中的最小值，l(m ,n)为单个数据点与聚类中心的加权距离，Li为所有数据点与聚类中心的加权距离之和， 为数据点参数值中的最大值， 为数据点坐标， 为第i个聚类中
心的坐标，m为数据种类，n为数据量， 为参数权重系数， 为第i个聚类中心参数值，a、b分别对应实际计算中选择的磁法与重力的权重系数。
[0052] 步骤S24：重复步骤S23，直至K个加权距离和达到最小值，确定K个聚类中心，并将两组数据分配到其距离最近的聚类中心，使得重力与磁法的处理结果被分为K组，确定多个重点区。
[0053] 在本申请实施例中，所述步骤S3中，确定找矿靶区的过程，具体为：步骤S31：对所述重点区内地质、钻探资料以及物探结果进行多元综合分析；步骤S32：推断由浅至深地层及构造分布情况；步骤S33：结合矽卡岩型富铁矿成矿因素，预测重要找矿标志。根据找矿靶区进行钻探验证，对预测结果进行修正。
[0054] 在本申请实施例中，所述重要找矿标志包括：大理岩化碳酸盐岩、角岩化碎屑岩、蚀变中基性岩、磁铁矿化蚀变、矽卡岩化蚀变、钠长石化蚀变、蛇纹石化蚀变、构造交汇部位、层间滑脱带、脆性裂隙。
[0055] 在本申请实施例中，所述找矿靶区岩体全铁含量较周边低3%以上、镁含量较周边低2%以上，钾含量较周边低1%以上，钠含量较周边高1.5%以上；围岩地层钙含量较周边低1%以上，铁含量较周边低1%以上。
[0056] 以上是本申请实施例提供的一种基于多重K均值聚类分析的找矿预测方法，基于同样的发明构思，本申请实施例还提供了一种基于多重K均值聚类分析的找矿预测系统，图2为本申请实施例提供的一种基于多重K均值聚类分析的找矿预测系统组成图，如图2所示，所述系统主要包括：有利区确定模块201，用于在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围；聚类模块202，用于基于K均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区；预测模块203，用于收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
[0057] 本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
[0058] 还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0059] 以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
