[发明专利]基于矿物地球化学的斑岩系统含矿性评价方法

权利要求书

1.一种基于矿物地球化学的斑岩系统含矿性评价方法，其特征在于，包括以下步骤：

1)圈定找矿靶区；

2)采集靶区内所有中酸性岩体，选择无蚀变或弱蚀变的样品，进行化学分析，得到全岩主微量元素SiO₂、Sr、Y、V、Sc含量，记为c(SiO₂)、c(Sr)、c(Y)、c(V)、c(Sc)；锆石微量元素Ce、Nd、Y、Eu、Sm、Gd含量，记为c(Ce)、c(Nd)、c(Y＇)、c(Eu)、c(Sm)、c(Gd)；角闪石温度T和氧逸度△FMQ；

3)以F1-F4为判别因子，进行成矿岩体与贫矿岩体判别区分，具体为：

判别因子F1：F1＝-0.925*c(SiO₂)+113.75(公式1)；

判别因子F2：F2＝-0.355*c(SiO₂)+34.15(公式2)；

判别因子F3：F3＝-41.7*(c(Ce)/c(Nd))/c(Y＇)+4.707(公式3)；

判别因子F4：F4＝-0.0025*c(T)+3.8(公式4)；

将步骤2)获得的全岩主微量元素SiO₂含量带入上述公式1，计算判别因子F1，当c(Sr)/c(Y)F1时，则判别为成矿岩体，反之则为贫矿岩体；

将步骤2)获得的全岩主微量元素SiO₂含量带入上述公式2，计算判别因子F2，当c(V)/c(Sc)F2时，则判别为成矿岩体，反之则为贫矿岩体；

将步骤2)获得的锆石(c(Ce)/c(Nd))/c(Y＇)比值带入上述公式3，计算判别因子F3，当10000*(c(Eu)/c(Eu*))/c(Y＇)F3时，则判别为成矿岩体，反之则为贫矿岩体，其中

将步骤2)获得的角闪石温度T带入上述公式4，计算判别因子F4，当△FMQF4时，则判别为成矿岩体，反之则为贫矿岩体；

当四个判别因子均判别得到成矿岩体，则为成矿岩体；当有至少一个判别因子判别得到贫矿岩体，则为贫矿岩体；其中：

获得判别因子F1-F4具体步骤为：

(1)分别选择无蚀变或弱蚀变的成矿岩体样品大于60件、贫矿岩体样品大于60件；

(2)对步骤(1)选择的样品进行化学分析，得到全岩主微量元素SiO₂、Sr、Y、V、Sc含量，记为c(SiO₂)、c(Sr)、c(Y)、c(V)、c(Sc)；锆石微量元素Ce、Nd、Y、Eu、Sm、Gd含量，记为c(Ce)、c(Nd)、c(Y＇)、c(Eu)、c(Sm)、c(Gd)；角闪石温度T和氧逸度△FMQ；

(3)基于步骤(2)所得大量矿区发育的成矿和贫矿岩体的数据，绘制散点图，依据成矿岩体和贫矿岩体分界点附近多个数据点拟合为一条直线，求得该直线的方程即为判别因子，具体为：

通过拟合并计算得到判别因子F1-F4，具体为：

①计算判别因子F1

对步骤(2)获得的全岩主微量元素，以c(SiO₂)为横坐标c(Sr)/c(Y)为纵坐标进行投图，依据投图范围获得成矿岩体与贫矿岩体的分界线并进行拟合，利用如下公式计算得到判别因子F1：F1＝-0.925*c(SiO₂)+113.75(公式1)；

②计算判别因子F2

对步骤(2)获得的全岩主微量元素，以c(SiO₂)为横坐标，c(V)/c(Sc)为纵坐标进行投图，依据投图范围获得成矿岩体与贫矿岩体的分界线并进行拟合，利用如下公式计算得到判别因子F2：F2＝-0.355*c(SiO₂)+34.15(公式2)；

③计算判别因子F3

对步骤(2)获得的锆石微量元素，以(c(Ce)/c(Nd))/c(Y＇)为横坐标，10000*(c(Eu)/c(Eu*))/c(Y＇)为纵坐标进行投图，其中依据投图范围获得成矿岩体与贫矿岩体的分界线并进行拟合，利用如下公式计算得到判别因子F3：F3＝-41.7*(c(Ce)/c(Nd))/c(Y＇)+4.707(公式3)；

④计算判别因子F4

对步骤(2)获得的角闪石成分，计算其形成温度及氧逸度，以T(温度)为横坐标，△FMQ(氧逸度)为纵坐标进行投图，依据投图范围获得成矿岩体与贫矿岩体的分界线并进行拟合，利用如下公式计算得到判别因子F4：F4＝-0.0025*T+3.8(公式4)。