何其芬，尹維民，張秀文，孫波，邱愷毅

(山東省物化探勘查院，山東 濟(jì)南　250013)

?

R型因子分析和聚類分析在內(nèi)蒙古伊山林場(chǎng)土壤地球化學(xué)測(cè)量中的應(yīng)用

何其芬，尹維民，張秀文，孫波，邱愷毅

(山東省物化探勘查院，山東 濟(jì)南250013)

土壤地球化學(xué)測(cè)量分析元素眾多，所成圖件數(shù)量龐大，因此往往不易判斷各元素之間的相關(guān)關(guān)系。筆者利用R型因子分析和聚類分析方法，對(duì)內(nèi)蒙古伊山林場(chǎng)金礦普查區(qū)土壤樣品進(jìn)行分析，研究其元素間的共生組合關(guān)系。結(jié)果顯示，R型因子分析和聚類分析二者對(duì)應(yīng)性較好，有效地劃分了元素共生組合類型。

土壤地球化學(xué)測(cè)量；因子分析；聚類分析；伊山林場(chǎng)；內(nèi)蒙古

引文格式：何其芬，尹維民，張秀文，等.R型因子分析和聚類分析在內(nèi)蒙古伊山林場(chǎng)土壤地球化學(xué)測(cè)量中的應(yīng)用[J].山東國土資源，2015,31(8)：61-64.HE Qifen ，YIN Weimin，ZHANG Xiuwen,etc. Application of R Type Factor Analysis and Cluster Analysis in Soil Geochemical Survey at Inner Mongolia of Yishan Forest Farm[J].Shandong Land and Resources, 2015,31(8):61-64.

伊山林場(chǎng)金礦普查區(qū)位于內(nèi)蒙古自治區(qū)東北部原始森林區(qū)，由于普查區(qū)植被覆蓋嚴(yán)重，基巖出露極少，多年來開展了大量的找礦工作都沒有取得明顯突破，覆蓋區(qū)找礦已經(jīng)成為礦產(chǎn)勘查界公認(rèn)的難題。我國近年來一些找礦實(shí)踐表明[1-4]，土壤地球化學(xué)方法能有效縮小找礦靶區(qū)，圈定異?；拘螒B(tài)和規(guī)模，查明異常源，對(duì)異常的成礦遠(yuǎn)景作出初步評(píng)價(jià)，為后續(xù)礦產(chǎn)勘查工作提供可靠的地球化學(xué)依據(jù)。然而土壤測(cè)量分析元素多，所成圖件數(shù)量龐大，往往不易判斷各元素之間的相關(guān)關(guān)系，從而給后續(xù)的異常圈定造成困難。筆者通過R型因子分析和聚類分析方法，對(duì)普查區(qū)的元素共生組合類型進(jìn)行劃分，取得了良好的效果。

1　地質(zhì)概況與土壤測(cè)量方法

普查區(qū)地層以中生代早白堊世梅勒?qǐng)D組為主，分布面積占60%以上，其次為中生代晚侏羅世滿克頭鄂博組和新生代第四紀(jì)全新世地層；出露巖漿巖為中生代侏羅紀(jì)中粒二長花崗巖，僅在普查區(qū)南北兩端有少量分布。由于普查區(qū)地表覆蓋嚴(yán)重，未見明顯斷裂構(gòu)造，僅在普查區(qū)北端有兩條物探推斷斷裂(圖1)①。

在普查區(qū)20km2范圍內(nèi)，按100m×40m網(wǎng)度進(jìn)行采樣，共采集樣品4628件，樣品主要采自25cm以下的“B”層(坡積層和殘積層)，部分采自“C”層(母質(zhì)層)，樣品的組成成分由褐色或黃褐色粘土、巖石碎屑組成。樣品充分曬干后進(jìn)行加工過篩，篩取-5～+40目樣品，過篩后樣品重量大于160g。樣品分析Au，Ag，Cu，Pb，Zn，Mo，As，Sb，W，Bi 10個(gè)元素。

2　分析方法及結(jié)果

由于各元素間有著千絲萬縷的聯(lián)系，如何有效地劃分元素共生組合類型，是圈定異常的關(guān)鍵。該文利用R型因子分析和聚類分析來判別它們的相關(guān)關(guān)系，從而確定元素共生組合類型。

2.1R型因子分析

R型因子分析是研究元素共生組合的有效手段和方法[4-7],其中每一個(gè)因子所包含主要元素,不僅僅表示它們的一種組合關(guān)系,而且反映了一種內(nèi)在的成因聯(lián)系[8]。筆者采用R型因子分析對(duì)原有10個(gè)元素(變量)進(jìn)行濃縮，提取有代表性的公共因子，根據(jù)因子負(fù)載矩陣中所反映的不同元素組合來確定各元素的相關(guān)關(guān)系，進(jìn)而劃分元素共生組合類型[9]。

1—現(xiàn)代沖積層:砂、礫；2—梅勒?qǐng)D組；3—滿克頭鄂博組；4—侏羅紀(jì)變粒二長花崗巖；5—地質(zhì)界線；6—不整合接觸界線；7—推斷斷層；8—普查區(qū)范圍圖1　普查區(qū)地質(zhì)概況

2.1.1因子分析的前提條件

因子分析的主要任務(wù)之一是對(duì)原有變量進(jìn)行濃縮，即將工作區(qū)10個(gè)元素中的信息重疊部分進(jìn)行提取并綜合成因子，進(jìn)而達(dá)到減少變量個(gè)數(shù)的目的，因此要求原有變量之間應(yīng)存在較強(qiáng)的相關(guān)關(guān)系，否則，如果原有變量之間相互獨(dú)立，不存在信息重疊，也就無法從中綜合出能夠反映某些變量共同特性的幾個(gè)較少的公共因子[10]。該文利用巴特利特球度檢驗(yàn)和KMO檢驗(yàn)對(duì)所選數(shù)據(jù)進(jìn)行相關(guān)關(guān)系檢驗(yàn)，概率P小于給定的顯著性水平a時(shí)，認(rèn)為原有變量適合做因子分析；KMO值越接近于1，意味著變量間的相關(guān)性越強(qiáng)[11]。工作區(qū)KMO值為0.779，概率P值為0，適合做因子分析。

2.1.2元素組合類型的確定

用全區(qū)4628件樣品10元素的分析數(shù)據(jù)做R型因子分析，前5個(gè)因子特征根代表的方差已占總方差的80%以上(表1)，因此視這5個(gè)因子為主要因子。由于正交旋轉(zhuǎn)因子負(fù)載矩陣比初始因子負(fù)載矩陣所反映的元素組合更具合理性和可解釋性[12]，因此該文采用了正交旋轉(zhuǎn)因子載荷矩陣(表2)來劃分元素組合，可以認(rèn)為這5個(gè)因子分別代表了工作區(qū)的5種元素組合類型，以因子載荷絕值γ>0.4的元素為該因子主要載荷元素，得出因子結(jié)構(gòu)式(表3)，即確定了5種元素組合類型：F1-As，Sb，Mo；F2-Pb，Ag，Bi；F3-Zn，Cu；F4-W；F5-Au。

2.2R型聚類分析

因子分析提取了各元素對(duì)區(qū)內(nèi)地球化學(xué)變差的貢獻(xiàn)，但各元素間的親疏關(guān)系不明，利用原始數(shù)據(jù)做聚類分析反映各元素間的親疏關(guān)系，研究各因子分配情況(圖2)。R型聚類分析在γ=0.34的相似水平上，可分為5大聚類：Ⅰ簇-Cu，Zn；Ⅱ簇-Pb，Bi，Ag；Ⅲ簇-W；Ⅳ簇-Mo，As，Sb；Ⅴ簇-Au。

表1　因子解釋原有變量總方差情況

表2　正交因子載荷矩陣

表3　因子結(jié)構(gòu)式

圖2　R型聚類分析圖

2.3元素組合特征分析

聚類分析中的Ⅰ簇元素組合為Cu，Zn，對(duì)應(yīng)因子分析中的變量F3；Ⅱ簇元素組合為Pb，Bi，Ag，與F2相對(duì)應(yīng)，兩種元素組合為區(qū)內(nèi)主要多金屬礦化元素，說明普查區(qū)內(nèi)存在較強(qiáng)的熱液活動(dòng)。區(qū)域上已發(fā)現(xiàn)多處銅、鉛、鋅礦點(diǎn)，主要為中溫?zé)嵋盒偷V床。從礦床地質(zhì)特征看，該類型礦床主要分布在中酸性巖體與地層接觸帶附近。普查區(qū)發(fā)育侏羅紀(jì)變粒二長花崗巖，該巖體與晚白堊世梅勒?qǐng)D組、晚侏羅世滿克頭鄂博組的接觸帶和土壤地球化學(xué)異常是尋找該類型礦床的標(biāo)志。

Ⅲ簇元素組合僅為W元素，對(duì)應(yīng)變量F4。由表2看出，W元素在其他變量中的載荷要遠(yuǎn)遠(yuǎn)小于F4，充分說明W元素在該區(qū)的富集具有較大的獨(dú)立性。

Ⅳ簇元素組合為Mo，As，Sb，對(duì)應(yīng)變量F1，該組合主要反映區(qū)內(nèi)的巖漿作用，但不排除獨(dú)立成礦的可能性。

Ⅴ簇元素組合為Au元素，對(duì)應(yīng)變量F5，充分說明Au元素本身離散的特性，但區(qū)內(nèi)金多處高值點(diǎn)，也有可能在局部富集成礦。

3　應(yīng)用效果

筆者用樣品化驗(yàn)數(shù)據(jù)，繪制了F1-As，Sb，Mo；F2-Pb，Ag，Bi；F3-Zn，Cu 3張組合異常圖，結(jié)果顯示，異常明顯，元素套合程度高，有效的圈定了組合異常，為下一步的異常查證工作指明了方向，這也充分證明了R型因子分析和聚類分析劃分元素共生組合類型的有效性。

4　結(jié)論

(1)R型因子分析和聚類分析是劃分元素共生組合的有效方法，二者對(duì)應(yīng)性較強(qiáng)，通過對(duì)比，大大增強(qiáng)了劃分元素共生組合類型的準(zhǔn)確性。

(2)由R型因子分析和聚類分析得知，普查區(qū)10種元素可劃分為5種元素共生組合類型：As，Sb，Mo；Pb，Ag，Bi；Zn，Cu；W；Au，為圈定組合異常類型提供了依據(jù)。

(3)Cu，Zn，Pb，Bi，Ag為區(qū)內(nèi)主要的多金屬礦化元素，中酸性巖體與地層接觸帶附近為成礦有利部位，建議對(duì)接觸帶上的該類異常區(qū)進(jìn)行重點(diǎn)查證。

[1]付山嶺，胡斌，龔玉爽，等.白沙縣長塘嶺鎢多金屬礦區(qū)土壤地球化學(xué)特征及其指示意義[J]．物探與化探，2012，36(5)：723-727.

[2]呂軍，王建民，王洪波，等．土壤地球化學(xué)測(cè)量在三道彎子金礦床的應(yīng)用[J]．物探與化探，2005,29(6)：515-518.

[3]刁理品，韓潤生，方維萱．溝系土壤地球化學(xué)測(cè)量在貴州普晴銻金礦勘查區(qū)應(yīng)用于找礦效果[J]．地質(zhì)與勘探，2010，46(1)：120-127.

[4]于林松，邱成貴，劉偉．青海省北部綠草山地區(qū)化探元素組合分類信息探討[J]．山東國土資源，2014,30(2)：54-57.

[5]藏興運(yùn)，王建新，趙利剛，等．火山巖地區(qū)土壤地球化學(xué)測(cè)量數(shù)據(jù)的處理與找礦——以琿春柳樹河子金銅研究區(qū)為例[J]．黃金，2007，28(4)：10-13.

[6]吳錫生.化探數(shù)據(jù)處理方法[M].北京:地質(zhì)出版社，1993：38-39.

[7]趙鵬大.定量地質(zhì)學(xué)理論與方法[M].北京:地質(zhì)出版社，2004.178-180.

[8]董慶吉,陳建平,唐宇.R型因子分析在礦床成礦預(yù)測(cè)中的應(yīng)用——以山東黃埠嶺金礦為例[J]．地質(zhì)與勘探,2008，44(4)：64-68.

[9]董毅.因子分析在水系沉積物測(cè)量地球化學(xué)分區(qū)中的應(yīng)用探討[J]．礦產(chǎn)與地質(zhì)，2008，22(1)：78-82.

[10]董毅,范麗琨,段煥春,王寶勛.青海大坂山地區(qū)水系沉積物測(cè)量元素組合分區(qū)[J]．地質(zhì)與勘探，2009，45(1)：70-74.

[11]薛薇.SPSS統(tǒng)計(jì)分析方法及應(yīng)用[M].北京：電子工業(yè)出版社，2004：9-12.

[12]戴慧敏，鮑慶中，宮傳東，等.因子分析法對(duì)內(nèi)蒙古查巴奇地區(qū)水系沉積物地球化學(xué)分區(qū)的應(yīng)用研究[J]．現(xiàn)代地質(zhì)，2010，24(2)：245-251.

Application of R Type Factor Analysis and Cluster Analysis in Soil Geochemical Survey at Inner Mongolia of Yishan Forest Farm

HE Qifen，YIN Weimin,ZHANG Xiuwen,SUN Bo,QIU Kaiyi

(Shandong Geophysical and Geochemical Exploration Institute，Shandong Jinan 250013，China)

Soil geochemical surve analysis numerous elements and The resulting maps are huge, it is often difficult to determine the relationship between elements. The soil samples in the yi shan forest gold census area of Inner Mongolia were analyzed by using R Type Factor Analyses and Cluster Analysis,the writer tries to study the association relationship between the elements. Results show that the R Type Factor Analyses and Cluster Analysis has good correspondence between them,To divide the element assemblage.

Soil geochemical surve Factor Analyses Cluster Analysis

2014-08-11；

2014-08-30；編輯：曹麗麗

何其芬(1986—)，女(土家族)，貴州銅仁人，助理工程師，主要從事地質(zhì)礦產(chǎn)勘查工作；E-mail：heqifen1987@163.com

P208

B

①內(nèi)蒙古第六地質(zhì)礦產(chǎn)勘查開發(fā)有限責(zé)任公司，內(nèi)蒙古自治區(qū)鄂倫春自治旗伊山林場(chǎng)金礦普查，2012年。