(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210897117.9 (22)申请日 2022.07.28 (71)申请人 中南大学 地址 410000 湖南省长 沙市岳麓区麓山 南 路932号 (72)发明人 邓浩　谷晓堂　毛先成　 (74)专利代理 机构 长沙轩荣专利代理有限公司 43235 专利代理师 姚兴 (51)Int.Cl. G01N 33/24(2006.01) G06T 17/00(2006.01) (54)发明名称 基于勘查数据的岩浆岩型矿床岩体原始产 状定量推断方法 (57)摘要 本发明公开了一种基于勘查数据的岩浆岩 型矿床岩体原始产状定量推断方法， 具体包括以 下步骤： 基于原始勘查数据， 得到样品在岩浆岩 型矿床岩体中的空间定位； 将经过空间定位后的 数据采用分形含量 ‑总量法确定矿 元素品位值异 常界限； 将经过异常界限筛选的数据采用反距离 加权法进行估值后， 确定成矿元素在空间分布的 品位值梯度； 基于成矿元素的空间分布， 确定在 宏观上品位值变化最快的方向模拟成矿期重力 方向； 基于成矿期重力方向恢复并推断岩体原始 产状。 本发明利用 岩浆岩型矿床已有的大量的、 广泛分布的勘查数据， 从中挖掘出与含矿岩体原 始产状信息相关的成矿期重力方向， 为岩体原始 产状的研究提供了定量化的思路， 实现对含矿岩 体原始产状的推断。 权利要求书2页 说明书6页 附图2页 CN 115267132 A 2022.11.01 CN 115267132 A 1.一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 具体 包括以下步骤： 步骤一： 基于原 始勘查数据， 得到样品在岩浆岩型矿床 岩体中的空间定位； 步骤二： 将经过步骤一空间定位后的数据采用分形含量 ‑总量法确定成矿元素品位值 异常界限； 步骤三： 将经过步骤二异常界限筛选的数据采用反距离加权法进行估值后， 确定成矿 元素在空间分布的品位 值梯度； 步骤四： 基于成矿元素空间分布的品位值梯度， 确定在宏观上品位值变化最快的方向 模拟成矿期重力方向； 步骤五： 基于成矿期重力方向恢复 并推断岩体原 始产状。 2.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 所述步骤一中， 勘查数据包括钻孔孔号信息、 测斜信息、 采样深度、 采样样长、 样品平均品位 值。 3.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 所述步骤二具体为： 将空间定位后的数据按照成矿元素品位值进 行升序排列； 再按其指数规律将数据划分为若干级别， 并统计大于分组级别的数据总量及其与分组级别 的比值——相对组成个数； 最后以分组级别为横轴， 对应的相对组成个数为纵轴组成双对 数坐标曲线， 使用最小二乘法在曲线图中进行分段拟合， 确定成矿元 素品位值的异常界限。 4.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 所述步骤三中， 将经过步骤二异常界限筛选的空间中存在样品的某处按照如 下(1)式采用反距离加权法进行估值； 并按(2) ‑(4)式使用定义法求该处的梯度t＝(tx,ty, tz)： 其中， I(x,y,z)为三维空间中某处待估的成矿元素品位值， I(xi,yi,zi)为待估点邻域 内任一成矿元 素品位值， di采用距离平方。 5.根据权利要求4所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 所述待估点邻域内为以待估点 为球心， 勘探线间距为半径的球形 领域内。权　利　要　求　书 1/2 页 2 CN 115267132 A 26.根据权利要求4所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 所述步骤四具体为： 取空间中任一方向d计算方向导数， 当且仅当方向导数为 正时参与以下求和： 其中， m为经过异常界限筛选后的成矿元素品位数据总量， d为单位方向向量， α 为方向d 的水平投影所对应的方位角， β 为方向d倾伏角， α、 β 均取整数角度， tk＝(tx,ty,tz)k为第k个 数据点的梯度； 当e取最大值emax时， 其所对应的方向 为成矿元素的空间分布在宏观上品位值 变化最快的方向， 并以该 方向模拟成矿期重力方向g， 即 7.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法， 其特征在于， 所述步骤五具体为： 以重心 为旋转中心， 将岩 体线框模 型从成矿期 重力方向旋 转至现今竖直向下的重力方向， 即获得岩体原 始产状。权　利　要　求　书 2/2 页 3 CN 115267132 A 3