西藏改則縣木如地區(qū)巖屑地球化學(xué)分析

劉 洪 李光明 張智林 黃瀚霄 肖萬峰

(1.中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，四川 成都 610081； 2.西藏自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局第五地質(zhì)大隊(duì)，青海 格爾木816000；3.河南省地質(zhì)調(diào)查院，河南 鄭州450007)

西藏改則縣木如地區(qū)巖屑地球化學(xué)分析

劉 洪1李光明1張智林2黃瀚霄1肖萬峰3

(1.中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，四川 成都 610081； 2.西藏自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局第五地質(zhì)大隊(duì)，青海 格爾木816000；3.河南省地質(zhì)調(diào)查院，河南 鄭州450007)

利用木如地區(qū)1∶1萬巖屑地球化學(xué)測量數(shù)據(jù)，運(yùn)用因子分析等多種統(tǒng)計方法，研究該地區(qū)地球化學(xué)特征，共劃分出Sb-Pb-Ag、Cu-Zn-Bi-As、W-Au-Mo 3個元素組合，繪制因子得分等值線圖，并在此基礎(chǔ)上探討找礦方向，縮小了找礦靶區(qū)。結(jié)合該地區(qū)地質(zhì)特征認(rèn)為：與該地區(qū)金成礦有關(guān)的成礦熱液應(yīng)該為巖漿熱液，具有沿著北東向斷裂，由南西向北東運(yùn)移的趨勢；在木如地區(qū)除了尋找與金礦體外，還有尋找與巖漿熱液有關(guān)的Cu-Zn礦體的潛力。

木如地區(qū) 多元統(tǒng)計 地球化學(xué)分析

相對于傳統(tǒng)的數(shù)據(jù)處理方法，運(yùn)用多元統(tǒng)計分析方法能有效提高數(shù)據(jù)信息利用率，卿成實(shí)[1]等利用多元統(tǒng)計分析在扎西康礦床圍巖元素研究中取得了良好效果。本次利用西藏阿里改則縣洞措鄉(xiāng)境內(nèi)的木如地區(qū)地表6 221件1∶1萬巖屑地球化學(xué)測量數(shù)據(jù)進(jìn)行R型因子分析、R型聚類分析等多元統(tǒng)計分析，研究該地區(qū)的地球化學(xué)元素組合，探討找礦方向。

1 地質(zhì)概況

木如地區(qū)位于班公湖—怒江成礦帶的中段，近東西向改則—尼瑪斷裂(班公湖—怒江縫合帶南界)與北東向的日干配錯—扎西錯斷裂的交匯區(qū)域。按中國大地構(gòu)造劃分方案[2-3]，研究區(qū)屬于西藏—三江造山系(Ⅰ級)、班公湖—怒江—昌寧—孟連結(jié)合帶(Ⅱ級)、班公湖—怒江俯沖增生雜巖帶(Ⅲ級)內(nèi)。區(qū)域地層時代主要為中生代，其中，該地區(qū)廣泛出露的中下侏羅統(tǒng)木嘎崗日群(J1-2m)集中了較豐富的Au、Ag、Cu等元素[4]，已發(fā)現(xiàn)達(dá)查[5]、商旭[6-8]、拉嘎等多處巖金礦床(點(diǎn))。同時，木如地區(qū)具有大量砂金開采遺跡，具有一定的金找礦潛力。

木如地質(zhì)出露的木嘎崗日群地層呈近東西向分布。根據(jù)巖性組合特征，可以劃分為3個巖性段：第一巖性段(J1-2m1)為一套淺紅色中粒長石石英砂巖；第二巖性段(J1-2m2)為一套灰色至淺灰色含炭變砂巖與粉砂泥質(zhì)板巖互層；第三巖性段(J1-2m3)為一套灰白色中粒變石英砂巖。該地區(qū)構(gòu)造總體上屬于一個軸面向南傾的背斜，并在褶皺上發(fā)育一系列的脆性斷裂，脆性斷裂可以分為2期，早期斷裂為一系列北東向擠壓破碎帶(F101，F(xiàn)102)；晚期(F201)為切穿地層和東北向斷裂的北西向斷裂。該地區(qū)內(nèi)無巖漿巖發(fā)育，但附近的達(dá)查金礦發(fā)育小型中酸性巖脈[5]。

2 元素地球化學(xué)分析

2.1 采樣情況

本次地球化學(xué)測量方法選擇為巖屑地球化學(xué)測量，即采樣層位為C層(母質(zhì)層)，布置測線距為100 m，測點(diǎn)距為40 m。在木如地區(qū)共采集6 221個樣品，分析了包括Au、Ag、As、Sb、Cu、Pb、Zn、W、Mo、Bi共10種元素。

2.2 因子分析

化探數(shù)據(jù)能否用于因子分析，首先應(yīng)判斷其是否通過KMO檢驗(yàn)和Bartlett球度檢驗(yàn)。本次數(shù)據(jù)KMO值為0.728(KMO≥0.9者為非常適合因子分析，KMO 為0.8～0.9者為適合因子分析，KMO<0.6者為不適合因子分析[9])；Bartlett卡方值為29 097.8，在自由度為45的條件下，概率值為0.000(小于0.05)。故本次數(shù)據(jù)滿足因子分析的檢驗(yàn)要求，適合進(jìn)行因子分析。

通過對地球化學(xué)測量數(shù)據(jù)的因子分析，抽取得出了4個主因子，這4個主因子累計貢獻(xiàn)率達(dá)到了97.82%，即4個濃縮的變量共解釋了原有10個變量所有信息的97.82%，且旋轉(zhuǎn)前后總的信息量(即總的累計貢獻(xiàn)率)與旋轉(zhuǎn)前相比，沒有發(fā)生信息丟失，信息提取合理完整，見表1。根據(jù)旋轉(zhuǎn)因子載荷得到4個因子變量對應(yīng)的元素組合類型，F(xiàn)1代表Ag、Pb、Sb元素組合；F2代表Bi、Cu、Zn元素組合；F3代表As元素；F4代表Au、W、Mo元素組合，見表2。

表1 工作區(qū)R型因子分析特征值和累計方差貢獻(xiàn)率

2.3 聚類分析

對本次數(shù)據(jù)進(jìn)行R型聚類分析，結(jié)果顯示：Sb、Pb、Ag元素，Cu、Zn、Bi、As元素，W、Au、Mo元素分別聚為一類，與因子分析結(jié)果基本一致，不同的是，因子分析中，F(xiàn)2因子為Cu、Zn、Bi元素組合，而F3因子為As元素，見圖2。結(jié)合因子分析，把木如地區(qū)地球化學(xué)組合分為：Sb、Pb、Ag元素組合(F1因子)，Cu、Zn、Bi、As元素組合(F2+F3因子)和W、Au、Mo元素組合(F4因子)。

表2 木如地區(qū)R型因子分析正交旋轉(zhuǎn)因子載荷矩陣

圖1 木如地區(qū)元素(變量)聚類分析

2.4 元素組合特征

在因子分析中，除了得出因子變量對應(yīng)的元素組合外，還得到了每個采樣點(diǎn)對各個因子變量的得分，根據(jù)每個采樣點(diǎn)的因子得分，繪制了木如地區(qū)的各因子得分等值線圖。由該圖可知，Sb、Pb、Ag元素組合(F1因子)地球化學(xué)強(qiáng)度較弱，呈零星分布于該地區(qū)中西部，具有呈近北東向的線性特征，見圖2；Cu、Zn、Bi、As元素組合(F2+F3因子)地球化學(xué)強(qiáng)度較高，主要分布于該地區(qū)北東部、南西部及南東部，其中高異常區(qū)位于北東部，整體分布特征為2條近平行的近北東的線狀帶，見圖3；W、Au、Mo元素組合(F4因子)地球化學(xué)強(qiáng)度相對較強(qiáng)，主要分布于北東部和南西部，高異常值出現(xiàn)在南西部，整體分布具有北東向的線性特征，見圖4。各元素組合(公共因子)得分等值線大致沿木如地區(qū)北東向的F101和F102斷裂分布，顯示出成礦元素運(yùn)移受北東構(gòu)造控制的特征。

圖2 木如地區(qū)F1因子(Sb-Pb-Ag)得分等值線

圖3 木如地區(qū)F2+F3因子(Cu-Zn-Bi-As)得分等值線

圖4 木如地區(qū)F4因子(W-Au-Mo)得分等值線

3 討 論

Sb、Pb、Ag元素組合(F1因子)的方差貢獻(xiàn)率為31.75%，為低溫元素組合，因子得分等值線呈零星分布，在本區(qū)的地球化學(xué)強(qiáng)度較弱，其元素分布規(guī)律具近北東的線性特征，北東方向地球化學(xué)強(qiáng)度略高，南西方向較低，可代表金礦體前緣暈—近礦暈元素的特征[10-11]。

Cu、Zn、Bi、As元素組合(F2+F3因子)的方差貢獻(xiàn)率為51.05%，為本地區(qū)占最主要地位的公共因子，地球化學(xué)強(qiáng)度極高，該元素組合表現(xiàn)為中低溫元素組合，分布范圍較大，表現(xiàn)為2條近平行的近北東的線狀帶，北東地球化學(xué)強(qiáng)度最高，南西方向略低，可代表金礦體近礦暈元素的特征[10-11]。Cu、Zn元素均表現(xiàn)為親硫性，遷移能力強(qiáng)。該元素組合的分布可能與地表的褐鐵礦化、方鉛礦化蝕變作用有關(guān)。

W、Au、Mo元素組合(F4因子)的方差貢獻(xiàn)率為17.68%，地球化學(xué)強(qiáng)度較高，為中高溫元素組合，W、Mo地球化學(xué)性質(zhì)相近，易在酸性巖體中富集，該元素組合可代表與巖漿熱液有關(guān)的元素組合，因此，本地區(qū)與金有關(guān)的成礦熱液可能與巖漿作用有關(guān)，這與前人對木如地區(qū)附近的達(dá)查金礦羅布日俄么礦段和達(dá)查金礦屋素拉礦段的成礦熱液的研究結(jié)果一致。該元素組合與另外2個元素組合一樣，整體表現(xiàn)為北東向的線性分布特征。該組合可代表金礦體近礦暈—礦尾暈的特征[10-11]，大面積的異常主要分布在北東部，但是高異常均分布在南西部。該地球化學(xué)組合與金礦體直接相關(guān)，因子得分等值線分部與北東向的F101和F102斷裂展布一致，因此，根據(jù)現(xiàn)有研究成果，尋找金礦體應(yīng)該在F101和F102斷裂及附近。

各元素組合呈近北東向線性分布，與區(qū)內(nèi)的北東向斷裂(F101和F102)分部一致。分析表明，本地區(qū)與金有關(guān)的成礦熱液應(yīng)該與巖漿作用有關(guān)。而前緣暈—近礦暈(Sb、Pb、Ag元素)組合高地球化學(xué)強(qiáng)度區(qū)在北東部、，近礦暈(Cu、Zn、Bi、As元素)組合高地球化學(xué)強(qiáng)度區(qū)比前緣暈組合偏南西部位，近礦暈—礦尾暈(W、Au、Mo元素)組合高地球化學(xué)強(qiáng)度區(qū)在更南西的部位。因此，與本地區(qū)金成礦具有關(guān)的巖漿熱液具有沿著北東向斷裂(F101、F102)、由南西向北東運(yùn)移的趨勢。此外，Cu、Zn、Bi、As元素組合(F2+F3因子)在本研究區(qū)地球化學(xué)強(qiáng)度最高，因此，在木如地區(qū)除了尋找與金礦體外，還有尋找與巖漿熱液有關(guān)的Cu、Zn礦體的潛力。

4 結(jié) 論

根據(jù)現(xiàn)有研究成果，尋找金礦體應(yīng)該在F101和F102斷裂及附近；與本地區(qū)金成礦有關(guān)的成礦熱液應(yīng)該為巖漿熱液，具有沿著北東向斷裂，由南西向北東運(yùn)移的趨勢；在木如地區(qū)除了尋找與Au礦體外，還有尋找與巖漿熱液有關(guān)的Cu、Zn礦體的潛力。

[1] 卿成實(shí),丁 俊,李光明,等．西藏隆子縣扎西康礦床地層圍巖元素分布特征[J]．金屬礦山,2014(6):97-100． Qing Chengshi，Ding Jun，Li Guangming，et al.Element distribution of the surrounding rock of Zhaxikang deposit in Longzi country，Tibet[J]．Metal Mine,2014(6):97-100．

[2] 潘桂棠，李興振，王立全，等．青藏高原及鄰區(qū)大地構(gòu)造單元初步劃分[J]．地質(zhì)通報，2002,21(11)：701-707． Pan Guitang，Li Xingzhen，Wang Liquan，et al.Preliminary division of tectonic units of the Qinghai-Tibet Plateau and its adjacent regions[J]．Geological Bulletin of China,2002,21(11)：701-707．

[3] 潘桂棠,肖慶輝,陸松年,等．中國大地構(gòu)造單元劃分[J]．中國地質(zhì),2009,26(1):1-28． Pan Guitang，Xiao Qinghui，Lu Songnian，et al.Subdivision of tectonic units in China[J].Geology in China,2009,26(1):1-28．

[4] 李勝榮,肖 潤,周 肅,等．西藏改則地區(qū)金成礦作用[J].礦床地質(zhì),2005,24(1)：1-14. Li Shengrong，Xiao Run，Zhou Su，et al.Gold mineralization in Gaize area，Tibet[J]．Mineral Deposits,2005,24(1)：1-14.

[5] 肖 潤,李勝榮,傅璐珈,等．西藏尼瑪縣達(dá)查金礦床同位素地球化學(xué)研究[J]．現(xiàn)代地質(zhì),2002,16(2)：165-169． Xiao Run，Li Shengrong，F(xiàn)u Lujia，et al.Isotope geochemistry of Dacha gold deposits in Tibet[J]．Geoscience,2002，16(2):165-169．

[6] 黃瀚霄,李光明,劉 波,等．藏北商旭造山型金礦床的發(fā)現(xiàn)及意義[J].礦床地質(zhì),2014,33(3)：486-496． Huang Hanxiao，Li Guangming，Liu Bo，et al.Discovery of Shangxu orogenic type gold deposit in northern Tibet and its significance[J]．Mineral Deposits,2014,33(3)：486-496.

[7] 肖萬峰,李光明,黃瀚霄,等．藏北商旭金礦床控礦因素及其找礦標(biāo)志[J]．黃金,2013,43(10):17-21． Xiao Wanfeng，Li Guangming，Huang Hanxiao，et al.Prospecting indicator and ore-controlling factors of Shangxu gold deposit in north Tibet[J]．Gold,2013,43(10):17-21．

[8] 劉 洪,黃瀚霄,張智林,等．藏北商旭金礦床地質(zhì)特征、控礦因素及發(fā)現(xiàn)意義[J]．礦物學(xué)報,2013,33(S2):337-338． Liu Hong，Huang Hanxiao，Zhang Zhilin，et al.The geological characteristics，ore-controlling factors and discovery significance of Shangxu gold deposit in north Tibet[J]．Acta Mineralogica Sinica,2013,33(S2):337-338．

[9] 趙少卿,魏俊浩,高 翔,等．因子分析在地球化學(xué)分區(qū)中的應(yīng)用:以內(nèi)蒙古石板井地區(qū)1∶5萬巖屑地球化學(xué)測量數(shù)據(jù)為例[J]．地質(zhì)科技情報,2012,31(2):27-34． Zhao Shaoqing，Wei Junhao，Gao Xiang，et al.Factor analysis in the geochemical subdivisions：taking 1∶50 000 debris geochemical survey in the Shibanjing area of Inner Mongolia as an example[J]．Geological Science and Technology Information,2012,31(2):27-34．

[10] 李 惠,張文華,劉寶林,等．金礦床軸向地球化學(xué)參數(shù)疊加結(jié)構(gòu)的理想模式及其應(yīng)用準(zhǔn)則[J]．地質(zhì)與勘探,1999(6):41-44． Li Hui，Zhang Wenhua，Liu Baolin，et al.Ideal model and principle of staching structure of axial geochemical parameter of gold deposit[J]．Geology and Exploration,1999(6):41-44．

[11] 劉 洪,呂新彪,李春誠,等.河南羅山金城金礦床成礦條件與深部找礦前景分析[J]．地質(zhì)與勘探,2013(2):265-273． Liu Hong，Lu Xinbiao，Li Chuncheng，et al.Metallogenic conditions and ore-searching prospect at depth of the Jincheng gold ore deposit in Luoshan county，Henan province[J]．Geology and Exploration,2013(2):265-273．

(責(zé)任編輯 鄧永前)

Geochemical Analysis of Rock Debris in Muru Area，Gerze County，Tibet

Liu Hong1Li Guangming1Zhang Zhilin2Huang Hanxiao1Xiao Wanfeng3

(1.Chengdu Center，China Geological Survey，Chengdu 610081，China；2.No.5 Geological Party，Tibet Geological Exploration Bureau，Golmud 816000，China；3.Henan Institute of Geology Survey，Zhengzhou 450007，China)

According to the geochemical data of debris at 1∶10 000 in Muru area，the multivariate statistics methods，like the factor analysis method is applied to investigate the geochemical characteristics of this area.Three element groups are marked out：Sb-Pb-Ag，Cu-Zn-Bi-As and W-Au-Mo.Furthermore，the contour plot of the factor score is drawn.Based on this，the prospecting direction is discussed so as to narrow the range of prospecting target.The geological features of this area indicated that the ore-forming fluid which is related to the gold mineralization sources from magmatic hydrotherm，and migrated from southwest to northeast along the NE trending fault.Besides discovering the gold ore bodies，there is the prospecting potential for finding the magmatic-hydrothermal related Cu-Zn ore bodies in the Muru area.

Muru area，Multivariate statistics，Geochemical analysis

2014-07-24

國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃(973計劃)項(xiàng)目(編號：2011CB403105)，中國地質(zhì)調(diào)查項(xiàng)目(編號： 12120114041501、1212011086074)，中國地質(zhì)調(diào)查局成都地質(zhì)礦產(chǎn)研究所青年學(xué)科基金(編號：所控基[2014]-03)。

劉 洪(1987—)，男，碩士研究生。

P624

A

1001-1250(2014)-11-105-04