【技术实现步骤摘要】
本申请涉及找矿预测，尤其涉及一种基于多重k均值聚类分析的找矿预测方法及系统。
技术介绍
1、我国富铁矿以矽卡岩型铁矿为主，以往受理论和方法限制，以查证高磁异常和就矿找矿为主，一直未取得明显的找矿突破；随着找矿工作向深部推进，深部信息的获取难度越来越大，深部微弱异常信号难以获取和识别，加之地球物理探测的多解性，严重影响了找矿的效率和精度，亟需突破传统找矿思路，创新地球物理探测技术，为富铁矿深部找矿提供支撑。
2、卡岩型铁矿一般赋存于侵入岩体与灰岩接触带内，我国大量中生代赋煤盆地内煤系之下发育稳定了灰岩地层，同时也发现局部侵入岩对煤层产生了一定的影响，灰岩与侵入岩体为矽卡岩型富铁矿成矿奠定了物质基础，在含煤盆地内分布的低缓磁异常成为寻找矽卡岩型富铁矿的重要新方向。但受上覆地层屏蔽影响，深部目标地质体异常信息微弱，相关研究中无系统性针对矽卡岩型富铁矿找矿预测方法；单一地球物理方法难以满足深部探测需求，同时使用多种参数进行识别异常时，存在人工拾取具体位置存在误差且效率较低等问题，亟需提高预测精度。
技术实现思路
1、本申请实施例提供了一种基于多重k均值聚类分析的找矿预测方法及系统，用以解决现有的找矿方法存在误差且效率较低的技术问题。
2、一方面，本申请实施例提供了一种基于多重k均值聚类分析的找矿预测方法，所述方法包括：步骤s1：在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围。步骤s2：基于k均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点
3、在本申请的一种实现方式中，确定有利区的过程，具体为：步骤s11：基于赋煤盆地内煤系变质及航磁异常信息，确定远景区。步骤s12：在所述远景区对重磁数据分别进行统计筛选，得到数据值异常区域。步骤s13：将δg垂向一阶导数与δt化极异常分别归一化后进行相关性分析，确定重磁同高异常范围，从而得到所述有利区；其中，所述δg垂向一阶导数为重力异常沿垂直方向求得的一阶导数，所述δt化极异常是与岩石层中的磁性物质相关的异常磁场。
4、在本申请的一种实现方式中，所述步骤s13之后，所述方法还包括：采集磁场值，并计算磁场总量异常δt数据。对所述磁场总量异常δt数据进行位场转换，以得到磁场参数；其中，所述磁场参数包括：δt 化极异常、δt 化极垂向一阶导数、δt 化极垂向二阶导数、δt化极总梯度模、δt 化极剩余异常、δt 化极标准差、δt 化极水平梯度模以及δt 化极下延不同深度数据。基于k均值聚类方法依次对所述磁场参数进行识别，得到各磁场参数对于矿体的反映，继而筛选出第一地质单元；其中，所述第一地质单元为成矿概率超过预设值的单元。
5、在本申请的一种实现方式中，所述方法还包括：采集重力值，并计算布格重力异常δg数据。对所述布格重力异常δg数据进行位场转换，以得到重力参数；其中，所述重力参数包括：剩余重力异常、δg垂向一阶导数、δg总梯度模、δg标准差、δg水平梯度模、δg下延不同深度数据。基于k均值聚类方法依次对所述重力参数进行识别，得到第二地质单元。
6、在本申请的一种实现方式中，所述步骤s12，具体包括：开展高精度重力测量、磁法测量，两相邻测线距离为100~250m，两相邻测点距离为两相邻测线距离的1/5~1/4。
7、在本申请的一种实现方式中，所述步骤s2中，确定重点区的过程，具体包括：步骤s21：确定k个聚类中心。
8、步骤s22：对磁场参数及重力参数设定不同权重系数。
9、步骤s23：随机选择某一数据点的位置作为单元中心点，计算所有数据点与该中心的加权距离和 li，同时确定与该中心点最远距离的数据点位置为第2个中心点，计算公式如下：
10、。
11、。
12、。
13、其中，为重力或磁法数据点参数值，为数据点参数值中的最小值， l(m ,n)为单个数据点与聚类中心的加权距离， li为所有数据点与聚类中心的加权距离之和，为数据点参数值中的最大值，为数据点坐标，为第i个聚类中心的坐标，m为数据种类，n为数据量，为参数权重系数，为第i个聚类中心参数值，a、b分别对应实际计算中选择的磁法与重力的权重系数。
14、步骤s24：重复步骤s23，直至k个加权距离和达到最小值，确定k个聚类中心，并将两组数据分配到其距离最近的聚类中心，使得重力与磁法的处理结果被分为k组，确定多个重点区。
15、在本申请的一种实现方式中，所述步骤s3中，确定找矿靶区的过程，具体为：步骤s31：对所述重点区内地质、钻探资料以及物探结果进行多元综合分析。
16、步骤s32：推断由浅至深地层及构造分布情况。
17、步骤s33：结合矽卡岩型富铁矿成矿因素，预测重要找矿标志。
18、在本申请的一种实现方式中，所述重要找矿标志包括：大理岩化碳酸盐岩、角岩化碎屑岩、蚀变中基性岩、磁铁矿化蚀变、矽卡岩化蚀变、钠长石化蚀变、蛇纹石化蚀变、构造交汇部位、层间滑脱带、脆性裂隙。
19、在本申请的一种实现方式中，所述找矿靶区岩体全铁含量较周边低3%以上、镁含量较周边低2%以上，钾含量较周边低1%以上，钠含量较周边高1.5%以上；围岩地层钙含量较周边低1%以上，铁含量较周边低1%以上。
20、另一方面，本申请还提供了一种基于多重k均值聚类分析的找矿预测系统，应用如前述的方法，所述系统包括：有利区确定模块，用于在赋煤盆地内确定有利区；其中，所述有利区用于划定找矿范围；聚类模块，用于基于k均值聚类分析对重磁多重位场处理结果进行识别，在所述有利区内确定重点区；预测模块，用于收集地质体矿化蚀变级元素含量变化信息，在所述重点区内确定找矿靶区。
21、本申请实施例提供的一种基于多重k均值聚类分析的找矿预测方法及系统，具有以下有益效果：针对赋煤区展开研究，在中生代赋煤断陷盆地中，概选出存在于含煤地层及煤系矿产高变质岩中的1:50万航磁低缓异常区，快速圈定富铁矿远景区；通过对物探资料进行聚类分析，识别次级异常，提高找矿预测精度。
本文档来自技高网...【技术保护点】
1.一种基于多重K均值聚类分析的找矿预测方法，其特征在于，所述方法包括：
2.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，确定有利区的过程，具体为：
3.根据权利要求2所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述步骤S13之后，所述方法还包括：
4.根据权利要求3所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述方法还包括：
5.根据权利要求2所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述步骤S12，具体包括：
6.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述步骤S2中，确定重点区的过程，具体包括：
7.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述步骤S3中，确定找矿靶区的过程，具体为：
8.根据权利要求7所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述重要找矿标志包括：大理岩化碳酸盐岩、角岩化碎屑岩、蚀变中基性岩、磁铁矿化蚀变、矽卡岩化蚀变、钠长石化蚀变、蛇纹
9.根据权利要求1所述的基于多重K均值聚类分析的找矿预测方法，其特征在于，所述找矿靶区岩体全铁含量较周边低3%以上、镁含量较周边低2%以上，钾含量较周边低1%以上，钠含量较周边高1.5%以上；围岩地层钙含量较周边低1%以上，铁含量较周边低1%以上。
10.一种基于多重K均值聚类分析的找矿预测系统，应用如权利要求1-9任意一项所述的方法，其特征在于，所述系统包括：
...【技术特征摘要】
1.一种基于多重k均值聚类分析的找矿预测方法，其特征在于，所述方法包括：
2.根据权利要求1所述的基于多重k均值聚类分析的找矿预测方法，其特征在于，确定有利区的过程，具体为：
3.根据权利要求2所述的基于多重k均值聚类分析的找矿预测方法，其特征在于，所述步骤s13之后，所述方法还包括：
4.根据权利要求3所述的基于多重k均值聚类分析的找矿预测方法，其特征在于，所述方法还包括：
5.根据权利要求2所述的基于多重k均值聚类分析的找矿预测方法，其特征在于，所述步骤s12，具体包括：
6.根据权利要求1所述的基于多重k均值聚类分析的找矿预测方法，其特征在于，所述步骤s2中，确定重点区的过程，具体包括：
7.根据权利要求1所述的基于多重k均值聚类...
【专利技术属性】
技术研发人员：朱裕振，闫冰，沈立军，张文艳，孙超，汝亮，张龙，单嘉琪，
申请(专利权)人：山东省煤田地质规划勘察研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。