基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法技术

本发明专利技术公开了一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，具体包括以下步骤：基于原始勘查数据，得到样品在岩浆岩型矿床岩体中的空间定位；将经过空间定位后的数据采用分形含量

【技术实现步骤摘要】
基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法


[0001]本专利技术涉及隐伏矿体三维预测研究
，具体涉及一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法。

技术介绍

[0002]成矿后构造运动对岩浆岩型矿床的扰动，往往会使得含矿岩体的产状发生改变，而含矿岩体的原始产状信息与矿床深边部及外围找矿息息相关。根据原始产状信息沿岩体深部延伸方向寻找隐伏含矿岩浆房是理想的深部找矿方案，同时原始产状信息的确定也有利于在岩体侧伏位置寻找到新的含矿岩体；现有关于含矿岩体原始产状的研究多从传统地质角度如断裂、岩相、矿石组构等进行，缺乏定量化的分析且较少关注岩体的成矿期重力方向。而基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法将有助于补充这一研究。
[0003]大部分岩浆岩型矿床经过长期勘探和基建，都拥有的大量的、广泛分布的勘查数据，这些数据不仅表达了矿化空间分布的三维结构信息，还蕴含了与原始产状信息相关的重要线索——成矿期重力信息。凡重力分异凸显的成矿作用，从宏观上讲，成矿元素会在成矿期重力方向上富集，其所形成的空间分布也会呈现相应的富集规律。通过勘查数据分析研究矿床的矿化空间分布，捕捉其中蕴含的成矿期重力信息，即能帮助我们推断含矿岩体的原始产状信息。

技术实现思路

[0004]鉴于目前存在的上述不足，本专利技术提供一种具体涉及一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，本专利技术能够利用岩浆岩型矿床已有的大量的、广泛分布的勘查数据，从中挖掘出与含矿岩体原始产状信息相关的重要线索——成矿期重力方向，为岩体原始产状的研究提供了定量化的思路，是勘查数据利用的创造性尝试，实现对含矿岩体原始产状的推断，进一步能为矿床深边部及外围找矿提供可靠建议。
[0005]为了达到上述目的，本专利技术提供一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，具体包括以下步骤：
[0006]步骤一：基于原始勘查数据，得到样品在岩浆岩型矿床岩体中的空间定位；
[0007]步骤二：将经过步骤一空间定位后的数据采用分形含量
‑
总量法确定矿元素品位值异常界限；
[0008]步骤三：将经过步骤二异常界限筛选的数据采用反距离加权法进行估值后，确定成矿元素在空间分布的品位值梯度；
[0009]步骤四：基于成矿元素空间分布的品位值梯度，确定在宏观上品位值变化最快的方向模拟成矿期重力方向；
[0010]步骤五：基于成矿期重力方向恢复并推断岩体原始产状。
[0011]依照本专利技术的一个方面，所述步骤一中，勘查数据包括钻孔孔号信息、测斜信息、采样深度、采样样长、样品平均品位值。
[0012]依照本专利技术的一个方面，所述步骤二具体为：将经过步骤一空间定位后的数据按照成矿元素品位值进行升序排列；再按其指数规律将数据划分为若干级别，并统计大于分组级别的数据总量及其与分组级别的比值——相对组成个数；最后以分组级别为横轴，对应的相对组成个数为纵轴组成双对数坐标曲线，使用最小二乘法在曲线图中进行分段拟合，确定成矿元素品位值的异常界限。
[0013]依照本专利技术的一个方面，所述步骤三中，将经过步骤二异常界限筛选的空间中存在样品的某处按照如下(1)式采用反距离加权法进行估值；并按(2)
‑
(4)式使用定义法求该处的梯度t＝(t
x
，t
y
，t
z
)：
[0014][0015][0016][0017][0018]其中，I(x，y，z)为三维空间中某处待估的成矿元素品位值，I(x
i
，y
i
，z
i
)为待估点邻域内任一成矿元素品位值，d
i
采用距离平方。
[0019]依照本专利技术的一个方面，所述待估点邻域内为以待估点为球心，勘探线间距为半径的球形领域内。
[0020]依照本专利技术的一个方面，所述步骤四具体为：取空间中任一方向d计算方向导数，当且仅当方向导数为正时参与以下求和：
[0021][0022]其中，m为经过异常界限筛选后的成矿元素品位数据总量，d为单位方向向量，α为方向d的水平投影所对应的方位角，β为方向d倾伏角，α、β均取整数角度，t
k
＝(t
x
，t
y
，t
z
)
k
为第k个数据点的梯度；
[0023]当e取最大值e
max
时，其所对应的方向为成矿元素的空间分布在宏观上品位值变化最快的方向，并以该方向模拟成矿期重力方向g，即
[0024]依照本专利技术的一个方面，所述步骤五具体为：以重心为旋转中心，将岩体线框模型从成矿期重力方向旋转至现今竖直向下的重力方向，即获得岩体原始产状。
[0025]本专利技术的有益效果：本申请面向成矿后构造扰动的岩浆岩型矿床中的含矿岩体，使用分形含量
‑
总量法确定含矿岩体中成矿元素品位值异常界限，筛选出受重力分异作用
影响显著的勘查数据，并通过计算其在三维空间中的方向导数，利用成矿元素的空间分布在宏观上品位值变化最快的方向模拟成矿期重力方向，并进一步结合含矿岩体现现今的产状信息对其原始产状进行推断。本专利技术能够利用岩浆岩型矿床已有的大量的、广泛分布的勘查数据，从中挖掘出与含矿岩体原始产状信息相关的重要线索——成矿期重力方向，为岩体原始产状的研究提供了定量化的思路，是勘查数据利用的创造性尝试，实现对含矿岩体原始产状的推断，进一步能为矿床深边部及外围找矿提供可靠建议。
附图说明
[0026]图1(a)为本申请实施例1中Ⅲ矿区岩体的铜采样数据的可视化效果图；图1(b)为图1(a)对应的Ⅲ矿区岩体线框模型图；
[0027]图2为本申请实施例1的分组级别(R)为横轴，对应的相对组成个数为纵轴的双对数坐标曲线图；
[0028]图3为本申请实施例1的采样数据的方向导数并求和绘制出的三维散点模型图；
[0029]图4(a)为本申请实施例1未恢复至原始产出状态的Ⅲ矿区岩体线框模型图；图4(b)为本申请实施例1恢复至原始产出状态的Ⅲ矿区岩体线框模型图。
具体实施方式
[0030]为使本专利技术更加容易理解，下面结合具体实施例，进一步阐述本专利技术。应理解，这些实施例仅用于说明本专利技术而不用于限制本专利技术的范围，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。除非另有定义，下文所用专业术语和本领域专业技术人员所理解的含义一致。
[0031]一种具体涉及一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法：
[0032](1)基于原始勘查数据，得到样品在岩浆岩型矿床岩体中的空间定位
[0033]收集某一岩浆岩型矿床中某一矿区或含矿岩体区域内已有钻孔柱状图资料并数字化，提取柱状图中本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，其特征在于，具体包括以下步骤：步骤一：基于原始勘查数据，得到样品在岩浆岩型矿床岩体中的空间定位；步骤二：将经过步骤一空间定位后的数据采用分形含量
‑
总量法确定成矿元素品位值异常界限；步骤三：将经过步骤二异常界限筛选的数据采用反距离加权法进行估值后，确定成矿元素在空间分布的品位值梯度；步骤四：基于成矿元素空间分布的品位值梯度，确定在宏观上品位值变化最快的方向模拟成矿期重力方向；步骤五：基于成矿期重力方向恢复并推断岩体原始产状。2.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，其特征在于，所述步骤一中，勘查数据包括钻孔孔号信息、测斜信息、采样深度、采样样长、样品平均品位值。3.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，其特征在于，所述步骤二具体为：将空间定位后的数据按照成矿元素品位值进行升序排列；再按其指数规律将数据划分为若干级别，并统计大于分组级别的数据总量及其与分组级别的比值——相对组成个数；最后以分组级别为横轴，对应的相对组成个数为纵轴组成双对数坐标曲线，使用最小二乘法在曲线图中进行分段拟合，确定成矿元素品位值的异常界限。4.根据权利要求1所述的基于勘查数据的岩浆岩型矿床岩体原始产状定量推断方法，其特征在于，所述步骤三中，将经过步骤二异常界限筛选的空间中存在样品的某处按照如下(1)式采用反距离加权法进行估值；并按(2)
‑
(4)式使用定义法求该处的梯度t＝(t
x
...

【专利技术属性】
技术研发人员：邓浩，谷晓堂，毛先成，
申请(专利权)人：中南大学，
类型：发明
国别省市：