物化探測(cè)量元素基于R型聚類分析的分類

童 宇

(廣東省核工業(yè)地質(zhì)局二九二大隊(duì)，廣東河源517001)

1 概述

聚類分析是通過(guò)不同地質(zhì)體或地質(zhì)現(xiàn)象的某些共同屬性（親疏關(guān)系）的相似程度來(lái)進(jìn)行分類的一種多元統(tǒng)計(jì)分析方法。該方法已演變處許多種，這些聚類方法與判別分析不同點(diǎn)在于后者必須遵循“從已知到未知”的原則，而聚類分析則直接對(duì)“樣品”或變量進(jìn)行分類，然后進(jìn)行地質(zhì)解釋。

在物化探工作中，聚類分析主要是解決：①研究原生暈或次生暈中異常的共生組合及其礦體的關(guān)系；②評(píng)價(jià)巖體含礦性，評(píng)價(jià)或區(qū)分物化探異常點(diǎn)、帶和含礦帶；③了解成礦元素與哪些地質(zhì)因素有關(guān)；④應(yīng)用常量元素或其共生組合作為找礦標(biāo)志，尋找稀有元素的礦產(chǎn)等。

2 聚類的一般原則和繪制分群圖的約定準(zhǔn)則

2.1 聚類的一般原則

（1）相關(guān)類統(tǒng)計(jì)量由大到小成對(duì)逐次合并。

（2）若選出的兩樣品在分好的類中均未出現(xiàn)過(guò)，則把它們合并為一個(gè)獨(dú)立的新類；若有一個(gè)已經(jīng)出現(xiàn)在分好的類中，則另一個(gè)樣品與前一個(gè)樣品所在的類合并；若兩樣品都分別出現(xiàn)在已經(jīng)分好的類中，則這兩類合并成一個(gè)大類。

（3）若選出的兩樣品都出現(xiàn)在同一類中，則這對(duì)樣品就不再分類了。

（4）給定相似程度臨界水平。合并過(guò)程中，一旦相似水平小于臨界水平時(shí)，則停止合并。在此水平上可得到若干類。若不給定臨界水平，則到全部分類樣品合并成一類為止。

2.2 繪制分群圖的約定準(zhǔn)則

（1）在劃出的統(tǒng)計(jì)量坐標(biāo)軸上，對(duì)相關(guān)類統(tǒng)計(jì)量由大到小分度。

（2）連接各類時(shí)，遵循“保留較小序號(hào)”準(zhǔn)則。

（3）兩個(gè)單一樣品合并時(shí)，約定較大序號(hào)的樣品統(tǒng)一放在分群圖的左側(cè)（或上面）；兩個(gè)不同的類合并時(shí)，相似水平高的一類放在分群圖左側(cè)（或上面）；一個(gè)樣品與一類樣品合并時(shí)，一類樣品放在左側(cè)（或上面）。

3 相似系數(shù)矩陣

以XX盆地1∶5萬(wàn)水系沉積物為例，該水系沉積物共分析1363個(gè)樣品，分析元素有14種（Be、Bi、Co、Cu、La、Mn、Mo、Nb、Ni、Pb、Sn、W、Y、Zn）。利用分析結(jié)果計(jì)算相似系數(shù)矩陣，其計(jì)算公式如下：

通過(guò)計(jì)算可得相似矩陣：

即其相似系數(shù)矩陣為表1。

表1 水系沉積物14種元素相似系數(shù)矩陣表

4 一次形成法分類

根據(jù)前期所得的相似系數(shù)矩陣表1，用一次形成法分類，其過(guò)程如下：

(1)記下矩陣R中的最大值RY,La=0.6063,劃去Y行和Y列；

(2)記下矩陣R中剩余元素的最大值RNi,Co=0.5810，劃去Ni行和Ni列；

(3)記下矩陣R中剩余元素的最大值RSn,Nb=0.5475，劃去Sn行和Sn列；

(4)記下矩陣R中剩余元素的最大值RMn,Co=0.4531，劃去Mn行和Mn列；

(5)記下矩陣R中剩余元素的最大值RW,Nb=0.3499，劃去W行和W列；

(6)記下矩陣R中剩余元素的最大值RZn,Nb=0.3427，劃去Zn行和Zn列；

(7)記下矩陣R中剩余元素的最大值RPb,La=0.2810，劃去Pb行和Pb列；

(8)記下矩陣R中剩余元素的最大值RNb,Mo=0.2749，劃去Nb行和Nb列；

(9)記下矩陣R中剩余元素的最大值RCo,Bi=0.2032，劃去Co行和Co列；

(10)記下矩陣R中剩余元素的最大值RLa,Be=0.1571，劃去La行和La列；

(11)記下矩陣R中剩余元素的最大值RMo,Bi=0.1172，劃去Mo行和Mo列；

(12)記下矩陣R中剩余元素的最大值RBi,Be=0.0899，劃去Bi行和Bi列；

(13)記下矩陣R中剩余元素的最大值RCu,Be=0.0431，劃去Cu行和Cu列；

總結(jié)一次形成法分類，其結(jié)果見(jiàn)表2。

5 相關(guān)系數(shù)分群圖繪制

從表2可以看出：

表2 相似系數(shù)一次分群結(jié)果表

（1）相關(guān)系數(shù)均大于0.0431。

（2）可分成六級(jí)相關(guān)。Y、La，Ni、Co和Sn、Nb關(guān)系最為密切，為一級(jí)相關(guān)；其次Y、La和Pb，Ni、Co和Mn，Sn、Nb和W為二級(jí)相關(guān)；Y、La、Pb和Be，Ni、Co、Mn和Bi，Sn、Nb、W和Zn為三級(jí)相關(guān)；Be，Ni、Co、Mn和Bi，Sn、Nb、W、Zn和Mo為四級(jí)相關(guān)，相關(guān)系數(shù)為0.1172；Be、Ni、Co、Mn、Bi、Sn、Nb、W、Zn、Mo和Y、La、Pb、Be為五級(jí)相關(guān)，相關(guān)系數(shù)為0.0899；最后為六級(jí)相關(guān)，Cu和Be、Ni、Co、Mn、Bi、Sn、Nb、W、Zn、Mo、Y、La、Pb、Be相關(guān)，相關(guān)系數(shù)為0.0431。

依照繪制分群圖的約定準(zhǔn)則，可以劃出圖1。

6 地質(zhì)解釋

由分群圖圖1可以看出，1363個(gè)水系沉積物樣品間14種元素關(guān)系中，Y、La，Ni、Co和Sn、Nb關(guān)系最為密切，Cu和其它元素較為疏遠(yuǎn)。通過(guò)R型聚類分析結(jié)果與該盆地水系沉積物測(cè)量綜合異常圖相比較，發(fā)現(xiàn)各種元素的聚類情況基本一致，特別是Co、Ni，La、Y，Nb、Sn異常暈完全吻合。其它元素或重疊或相嵌，說(shuō)明分類結(jié)果是可信的，因此在探索某一元素礦化（體）的同時(shí)，必須考慮聚類中其它元素礦化（體）的可能。

圖1 相關(guān)系數(shù)一次形成法分群圖