第一篇：寧蕪盆地寧藏中段次火山巖含礦性貝葉斯兩組判別分析報(bào)告

4.2寧蕪盆地寧藏中段次火山巖含礦性貝葉斯兩組判別分析

寧蕪盆地位于淮陽(yáng)山字型構(gòu)造的前弧東翼與新華夏構(gòu)造復(fù)合處，為地跨蘇皖的中生代繼承式陸相火山巖盆地。盆地明顯受北北東及北西向斷裂控制，呈北北向東延伸。邊緣主要由三迭—侏羅系地層組成，中部則廣泛分布火山巖系?；鹕綆r系及少數(shù)下伏沉積巖系中發(fā)育有次火山巖，它們晚于該火山旋回侵入，并與火山巖在成分上相似。以Fe、Cu、S、P為主的多種礦產(chǎn)與次火山巖有直接關(guān)系。

在寧蕪中段取得31件次火山巖樣品，其中1-15為含礦次火山巖及晚期脈巖，16-31為非礦次火山巖及侵入巖。每個(gè)樣品分析SiO2、TiO2、Al2O3、Fe2O3、FeO、MnO、MgO、CaO、Na2O、K2O、H2O+等11項(xiàng)指標(biāo)。原始數(shù)據(jù)見表2-1。

表2-1 原始數(shù)據(jù)

樣品/變量 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 經(jīng)過30次迭代，得到有關(guān)參數(shù)如表2-2。表2-2 由迭代法求出的各極小值

角度FIMIN=2583996

函數(shù)值RMIN=9.228292E-02 一類判別函數(shù)均值=3.041054

二類判函數(shù)均值=-4.515806 一類判別函數(shù)方差=3.145484

二類判別函數(shù)方差=17.67206 此時(shí)的線性判別函數(shù)的系數(shù)如表2-3。樣品回代判別結(jié)果如表2-4。

從回代結(jié)果看，90%以上的樣品歸類適當(dāng)，僅有含礦類第二號(hào)樣品（輝長(zhǎng)閃安山玢巖，采自南山504孔），不含礦類第17號(hào)樣品（角閃安山玢巖，采自鶯子山龍王山組），第24號(hào)樣品（鈉長(zhǎng)斑巖，采自李洼大王山組）三件樣品歸屬有誤。具體原因待查。當(dāng)然，并不排除統(tǒng)計(jì)分析方法本身的概率風(fēng)險(xiǎn)。

此外，從線性判別函數(shù)（表2-3）可以看到，對(duì)于判含礦與不含礦時(shí)的11項(xiàng)分析指標(biāo)中，MnO、Na2O、MgO、H2O+及TiO2起重要作用。對(duì)于含礦，前四項(xiàng)起正作用，最后一項(xiàng)起負(fù)作用。其中關(guān)于Na2O的正作用可以理解為鈉化蝕變是與鐵的富集緊密相關(guān)聯(lián)的。表2-3 判 別 函 數(shù) 系 數(shù)（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

樣品回代和判別結(jié)果 樣品號(hào) 原類別 判別函數(shù)值 判別類別 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 樣品號(hào) 原類別 判別函數(shù)值 判別類別 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

四、參考文獻(xiàn)

1.潘恩沛《判別分析》講座筆記（未刊），1986年，西南冶金地質(zhì)科學(xué)研究所。2.趕漢地質(zhì)學(xué)院巖石教研室 寧蕪中段某些安山質(zhì)火山巖、次火山巖巖石學(xué)特征及其與玢巖鐵礦的關(guān)系，1978年，地質(zhì)出版社《寧蕪火山巖鐵銅礦床會(huì)議選集》，第154-177頁(yè)。/ 3