[发明专利]非线性离散推测模型找矿预测方法

摘要

本发明公开了一种非线性离散推测模型找矿预测方法，通过确定控矿要素的最佳组合数与最佳收缩距离，解决了满足较多控矿要素组合的区域范围小、较少控矿要素组合的区域范围大，收缩距离值过小导致最终圈定的成矿有利区个数多且分布离散、收缩距离过大导致所圈定的成矿有利区范围广而失去找矿预测研究意义等问题，从而在最佳控矿要素组合数与最佳收缩距离共同作用下克服矿区信息不对称以及连续插值模型的局限性，实现控矿要素最佳组合与最佳收缩距离双重约束下成矿有利区的圈定，达到找矿预测的目的。

主权项

1.非线性离散推测模型找矿预测方法，其特征在于包括以下步骤：步骤一：确定控矿要素最佳组合数(1)、块体计算：对研究区的地质背景以及成矿特征分析的基础上，总结研究区的找矿模型，确定控矿要素的类型，根据研究区现有的地质资料，构建局部三维地质实体模型，同时在网格划分的基础上构建三维块体模型。选取块体模型中的任一块体单元作为目标块体，并以该目标块体的中心位置作为圆心构建搜索椭球，根据研究区的范围大小，初始化设置搜索椭球的半径值；设置半径后，在搜索椭球范围内，分别统计任一不同控矿要素组合数下的块体单元(特征块体)到目标块体的距离值，并计算其最小距离值；(2)、距离(D)‑体积(V)分形拟合：将上述计算出的不同目标块体到特征块体的最小距离值，以控矿变量组合数为自变量，对同一控矿变量组合数下的所有块体单元按照其对应的最小距离值进行排序；对按照距离排序后的结果，统计某一距离值外所包含的的块体单元的个数(N)；由于在对研究区进行网格划分时，划分出的最小块体单元的尺寸(体积)是固定值(V0)，根据块体单元的个数与划分的尺寸计算出大于某一范围的块体的体积(V)；根据不同控矿要素组合的情况，按照公式1、公式2来分析距离和大于该距离的所有块体的体积的关系来拟合各种要素组合内和要素之间的自相似性V＝k×Dd   (公式1)其中：D为距离，V为大于D距离的块体体积，d为分形维数，k为一个比例常数对公式1取对数后得到公式2：ln V＝ln K+d ln D   (公式2)把数据投影到双log坐标系中，根据组间方差最小拟合出函数y＝a*x+b，所以分段拟合的分维数为d＝a；(3)、控矿要素组合分析：根据D‑V分形的结果，在满足不同个数的控矿要素组合下D‑V均满足分形分布特征的前提下，根据研究区内特征块体的分布范围，设置最大搜索距离(L)为搜索球的半径，在此数值设置上n等分最大搜索距离值L，分别构建以半径r＝L/n，2L/n，3L/n，4L/n，…的球体，其中最小半径长度的球体称为基础球体，其体积为V0，在此基础上构建等步长(L/n)变化的球体称为搜索球体，其体积计为Vi，分别统计相邻半径之间含有特定数目控矿变量组合的的特征块体的个数N，以及相邻搜索球体的体积增量，以单位体积内分布的特定数目的控矿要素特征块体的数目作为该数目控矿要素组合下的控矿密度(公式3)，以不同搜索椭球的控矿密度与基础椭球的控矿密度(公式4)之比为相对控矿密度(公式5)根据不同搜索尺度下相对控矿密度计算的结果，绘制不同控矿变量组合数下，相对控矿密度和和搜索尺度的统计曲线；步骤二：确定最佳收缩距离经上述研究确定出控矿要素的最佳组合数后，满足控矿变量最佳组合数的块体单元在找矿预测研究中都被称作是成矿有利单元。针对控矿要素最佳组合数约束后的成矿有利单元，非线性离散推测模型找矿预测方法中通过不同尺度的滑动窗体来实现成矿有利单元最佳收缩距离的确定；不同尺度的窗体对整个研究区进行滑动搜索的同时，滑动窗体中包含的控矿要素的组合情况也相应会发生变化，在不同尺度的窗体搜索下，满足最佳控矿要素组合数要求的块体单元称作成矿有利窗体，统计不同滑动搜索窗体的尺度与最终的成矿有利窗体个数之间的关系；在此基础上，选取成矿有利窗体个数峰值所对应的滑动搜索尺度，作为最佳控矿要素组合数下成矿有利单元收缩的最佳距离值；在确定控矿变量最佳组合数和最佳收缩距离的基础上，得到基于非线性离散推测模型找矿预测方法圈定出的成矿有利区。