宁芜盆地宁藏中段次火山岩含矿性贝叶斯两组判别分析报告

第一篇：宁芜盆地宁藏中段次火山岩含矿性贝叶斯两组判别分析报告

4.2宁芜盆地宁藏中段次火山岩含矿性贝叶斯两组判别分析

宁芜盆地位于淮阳山字型构造的前弧东翼与新华夏构造复合处，为地跨苏皖的中生代继承式陆相火山岩盆地。盆地明显受北北东及北西向断裂控制，呈北北向东延伸。边缘主要由三迭—侏罗系地层组成，中部则广泛分布火山岩系。火山岩系及少数下伏沉积岩系中发育有次火山岩，它们晚于该火山旋回侵入，并与火山岩在成分上相似。以Fe、Cu、S、P为主的多种矿产与次火山岩有直接关系。

在宁芜中段取得31件次火山岩样品，其中1-15为含矿次火山岩及晚期脉岩，16-31为非矿次火山岩及侵入岩。每个样品分析SiO2、TiO2、Al2O3、Fe2O3、FeO、MnO、MgO、CaO、Na2O、K2O、H2O+等11项指标。原始数据见表2-1。

表2-1 原始数据

样品/变量 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 54.14 52.48 56.42 55.48 57.12 45.86 43.84 52.28 54.78 60.66 57.98 51.24 64.22 65.25 58.68 52.90 55.94 56.69 54.34 53.07 58.52 58.66 54.72 54.02 62.34 73.46 72.25 66.20 67.36 60.26 2 0.89 0.88 0.82 0.80 0.75 0.88 0.60 0.66 0.89 0.64 0.81 0.85 0.65 0.55 0.65 0.92 0.86 1.06 0.82 0.80 1.01 1.00 0.85 0.68 0.55 0.02 0.25 0.23 0.36 0.65 17.86 17.30 19.59 16.65 17.98 14.24 16.18 15.35 20.12 15.78 16.65 16.76 17.17 17.15 15.50 18.80 16.84 15.24 18.54 17.45 20.15 19.43 18.06 18.70 15.28 13.08 13.94 13.49 14.59 15.59 3.85 5.20 3.63 4.54 4.25 13.80 2.85 7.04 1.55 4.86 4.60 9.72 1.82 1.93 3.53 6.35 3.19 11.55 9.61 8.55 1.12 1.76 8.25 6.70 2.39 1.62 1.23 1.18 1.62 2.66 4.39 6.03 2.25 5.23 3.84 7.83 3.88 1.74 3.13 3.12 4.64 7.61 1.33 2.05 6.22 3.58 5.19 1.53 0.55 3.61 1.73 1.73 3.46 3.18 5.71 1.11 1.82 1.86 2.34 4.33 0.15 0.17 0.10 0.09 0.09 0.08 0.08 0.04 0.04 0.08 0.06 0.05 0.03 0.03 0.20 0.42 0.20 0.01 0.02 0.08 0.04 0.11 0.08 0.07 0.08 0.03 0.05 0.08 0.06 0.06 3.46 2.57 1.54 3.03 2.04 2.56 2.22 2.83 2.84 2.14 3.14 2.90 1.69 1.66 2.09 2.09 2.94 1.23 0.76 1.90 2.64 2.55 1.10 1.62 2.14 0.17 0.52 1.51 1.65 2.22 6.20 5.04 4.48 4.08 3.44 2.92 1.72 5.06 3.84 0.85 0.94 0.55 1.98 1.92 1.34 6.54 4.27 0.51 0.30 4.72 4.94 4.26 2.07 1.34 1.95 0.11 1.09 3.56 2.61 3.27 5.36 5.72 7.82 6.00 6.25 6.44 5.00 4.00 6.14 4.70 5.31 8.00 8.25 8.04 3.90 2.89 4.12 3.34 0.65 6.14 7.00 6.71 7.08 6.70 3.25 0.48 0.77 0.26 3.90 4.12 1.24 1.64 0.78 2.24 2.83 1.36 3.48 2.18 2.19 3.20 1.09 0.25 0.98 0.55 4.10 2.41 3.00 4.79 11.54 0.83 0.77 0.63 1.88 1.46 4.00 7.59 6.32 4.80 4.39 4.78 1.80 0.57 1.78 1.79 1.30 0.97 2.04 2.22 2.04 2.28 3.56 1.94 1.72 0.68 2.54 1.19 2.12 2.91 1.40 1.60 1.32 0.68 1.94 0.76 2.26 2.34 1.54 2.59 0.46 0.80 / 3 31 48.12 1.20 17.28 4.93 6.16.0.18 4.86 9.20 3.80 2.01 0.86 经过30次迭代，得到有关参数如表2-2。表2-2 由迭代法求出的各极小值

角度FIMIN=2583996

函数值RMIN=9.228292E-02 一类判别函数均值=3.041054

二类判函数均值=-4.515806 一类判别函数方差=3.145484

二类判别函数方差=17.67206 此时的线性判别函数的系数如表2-3。样品回代判别结果如表2-4。

从回代结果看，90%以上的样品归类适当，仅有含矿类第二号样品（辉长闪安山玢岩，采自南山504孔），不含矿类第17号样品（角闪安山玢岩，采自莺子山龙王山组），第24号样品（钠长斑岩，采自李洼大王山组）三件样品归属有误。具体原因待查。当然，并不排除统计分析方法本身的概率风险。

此外，从线性判别函数（表2-3）可以看到，对于判含矿与不含矿时的11项分析指标中，MnO、Na2O、MgO、H2O+及TiO2起重要作用。对于含矿，前四项起正作用，最后一项起负作用。其中关于Na2O的正作用可以理解为钠化蚀变是与铁的富集紧密相关联的。表2-3 判 别 函 数 系 数（b 0）=64.3643（b 1）=-0.9825（b 2）=-18.3450（b 3）=-1.5013（b 4）=-0.8754（b 5）=-0.8883（b 6）=26.9320（b 7）=2.8232（b 8）=-0.3699（b 9）=4.0971（b10）=1.9838（b11）=2.5511 表2-4

样品回代和判别结果 样品号 原类别 判别函数值 判别类别 1 2 3 4 5 6 7 8 9 10 11 12 13 1 1 1 1 1 1 1 1 1 1 1 1 1 1.6365-0.8504 2.9691 4.7761 2.1364 2.8662 2.9704 3.9947 2.5380 1.8276 2.8434 5.8412 5.9430 2 1 1 1 1 1 1 1 1 1 1 1 / 3 14 15 1 1 1.6276 4.4959 1 1 样品号 原类别 判别函数值 判别类别 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2-6.8664 1.5847-10.8404-9.0242-2.4931-3.8847-4.5430-0.1080 0.4634-2.3713-9.2601-12.7419-5.3852-2.4846-1.0441-3.2543 1 2 2 2 2 2 2 1 2 2 2 2 2 2 2

四、参考文献

1.潘恩沛《判别分析》讲座笔记（未刊），1986年，西南冶金地质科学研究所。2.赶汉地质学院岩石教研室 宁芜中段某些安山质火山岩、次火山岩岩石学特征及其与玢岩铁矿的关系，1978年，地质出版社《宁芜火山岩铁铜矿床会议选集》，第154-177页。/ 3