張明明，龍瑾瀟，周國(guó)玉，焦峻嵚，方紅東

（1 合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院，安徽合肥 230009；2 合肥工業(yè)大學(xué)，礦床成因與勘查技術(shù)研究中心，安徽合肥 230009；3 合肥工業(yè)大學(xué)安徽省礦產(chǎn)資源與礦山環(huán)境工程技術(shù)研究中心，安徽合肥 230009；4 合肥工業(yè)大學(xué)空間信息集成與綜合分析平臺(tái)，安徽合肥 230009；5 安徽省地質(zhì)測(cè)繪技術(shù)院，安徽合肥 230022；6 宣城市華通礦業(yè)有限公司，安徽宣城 242000）

在當(dāng)前大數(shù)據(jù)的時(shí)代背景下，從數(shù)據(jù)出發(fā)，以數(shù)據(jù)驅(qū)動(dòng)模式去分析、挖掘地學(xué)研究積累的大量數(shù)據(jù)，可以從另一側(cè)面觀察、研究地質(zhì)問題，為地質(zhì)研究開拓了新的思路(Armstrong et al.,2014;Chen et al.,2014;邵燕等,2016;張旗等,2017;翟明國(guó)等,2018)。隨著礦產(chǎn)勘查方法與儀器分析測(cè)試技術(shù)的不斷進(jìn)步，地學(xué)研究已經(jīng)從傳統(tǒng)的定性研究向定量化研究快速發(fā)展(趙鵬大,2011;2014;2015)。地質(zhì)數(shù)據(jù)一般難以達(dá)到海量級(jí)數(shù)據(jù)的標(biāo)準(zhǔn)，但地學(xué)研究可以從數(shù)據(jù)出發(fā)，進(jìn)行數(shù)據(jù)之間各類關(guān)系的分析，挖掘地質(zhì)數(shù)據(jù)背后隱藏的信息，同樣是屬于大數(shù)據(jù)技術(shù)范疇(Li et al.,2009;Gantz et al.,2012;Carranza et al.,2015;Evans et al.,2016;Aryafar et al.,2017;Condie et al.,2017;Ernst et al.,2017;張旗等,2018a;2018b;2018c;Huang et al.,2021)。

通常，元素的遷移和富集規(guī)律研究、組合特征提取以及地化異常準(zhǔn)確圈定對(duì)于找礦勘查有著重要的意義(Zhao et al.,2012;Li et al.,2015;Wang et al.,2016)，如何快速分析巖石樣品中的元素成分及其含量，以判斷研究區(qū)的成礦潛力是一個(gè)非常關(guān)鍵的問題。實(shí)驗(yàn)室儀器分析方法由于測(cè)試的高精度使其在巖芯劈樣元素分析、礦石地化異常分析等分析中占據(jù)著主導(dǎo)作用，但實(shí)驗(yàn)室分析方法一般需要經(jīng)過(guò)“采集樣品-運(yùn)輸樣品-樣品加工-化學(xué)分析”等過(guò)程之后才能得到測(cè)試結(jié)果，過(guò)程耗時(shí)長(zhǎng)、費(fèi)用昂貴且需要破壞樣品，無(wú)法滿足對(duì)研究區(qū)進(jìn)行全數(shù)據(jù)采集(林延暢等,2002;Zhao et al.,2018;Cui et al.,2021)。便攜式X 射線熒光分析儀(pXRF)可以在野外快速、準(zhǔn)確地獲取原位巖芯數(shù)據(jù)，為野外礦床全鉆孔的元素采集提供了技術(shù)支持(Kalnicky et al.,2001;Hall et al.,2014;McAlister,2019)。

蕎麥山銅多金屬礦床位于長(zhǎng)江中下游成礦帶宣城礦集區(qū)內(nèi)，是區(qū)內(nèi)典型的矽卡巖型礦床(Mao et al.,2011;Pirajno et al.,2015;Nie et al.,2017;Liu et al.,2018)。前人對(duì)蕎麥山銅多金屬礦床的研究多從傳統(tǒng)地質(zhì)實(shí)驗(yàn)分析方面入手，采集礦床典型樣品進(jìn)行成巖成礦分析，成礦元素富集的特征及規(guī)律研究受到數(shù)據(jù)量及分析方法的限制。本次研究旨在從數(shù)據(jù)出發(fā)，運(yùn)用pXRF對(duì)蕎麥山銅礦床典型剖面進(jìn)行全鉆孔元素測(cè)試，以各元素之間的協(xié)同關(guān)系為基礎(chǔ)，研究各元素在深度空間上的分布趨勢(shì)和富集、虧損規(guī)律，深入分析元素之間相關(guān)性對(duì)礦物的指示，圈定地化異常，探索相關(guān)元素對(duì)成礦的指示作用。研究結(jié)果可以揭示蕎麥山銅多金屬礦床的元素分布特征，促進(jìn)對(duì)成礦過(guò)程中元素遷移、富集規(guī)律的認(rèn)識(shí)，明確矽卡巖型礦床的成礦模式和成礦條件，對(duì)進(jìn)一步的找礦勘查工作具有積極的指導(dǎo)意義(Gazley et al.,2014;Young et al.,2016;Bruno,2018;Perring et al.,2018;Gray et al.,2020)。

1 研究區(qū)地質(zhì)背景

長(zhǎng)江中下游成礦帶是中國(guó)最重要的銅鐵多金屬成礦帶之一，位于華北克拉通和揚(yáng)子克拉通的匯聚地帶，成礦帶內(nèi)有許多大規(guī)模、多類型的銅鐵多金屬礦床。蕎麥山銅多金屬礦床位于成礦帶內(nèi)宣城礦集區(qū)內(nèi)(圖1)。

長(zhǎng)江中下游成礦帶主要受特提斯構(gòu)造域、深部殼幔作用構(gòu)造過(guò)程和古太平洋構(gòu)造域的共同制約，經(jīng)歷了復(fù)雜的大陸動(dòng)力學(xué)過(guò)程和多階段的巖漿活動(dòng)。諸多因素控制了本區(qū)的構(gòu)造-巖漿-成礦作用，形成了許多大規(guī)模、多類型的銅鐵多金屬礦床(常印佛等,1991;翟裕生,1992;周濤發(fā)等,2008)。成礦帶自西向東主要包括鄂東南、九瑞、安慶-貴池、銅陵、廬樅、寧蕪和寧鎮(zhèn)7 個(gè)傳統(tǒng)的礦集區(qū)，以及近年來(lái)新增的宣城礦集區(qū)(圖1)(Pirajno et al.,2015;Nie et al.,2017;Liu et al.,2018)。

長(zhǎng)江中下游地區(qū)出露有晚太古代和古、中、新元古代變質(zhì)基底，前震旦系沉積變質(zhì)巖出露極少，震旦系—早三疊統(tǒng)穩(wěn)定沉積蓋層不整合覆蓋于揚(yáng)子克拉通不同基底之上，構(gòu)成了該區(qū)“一蓋多底”的地殼結(jié)構(gòu)(常印佛等,1996;陳志洪等,2017;朱安冬等,2017)。

雙基底層位是長(zhǎng)江中下游地區(qū)獨(dú)特的結(jié)構(gòu)，表現(xiàn)為南部是以新元古代復(fù)理石建造為主的江南式基底，北部是以太古宙—古元古代TTG 巖系為特征的崆嶺-董嶺式基底，長(zhǎng)江中下游成礦帶處于上述2種類型基底的邊界帶上(董樹文等,2011)。長(zhǎng)江中下游成礦帶從震旦系—白堊系都有各種類型的礦床產(chǎn)出，其中以銅、鐵、金、鎢等多金屬礦床為主，銅、金礦床主要產(chǎn)出于3 個(gè)層位，分別為上石炭統(tǒng)黃龍組、二疊系和中、下三疊統(tǒng)(包括東馬鞍山組、大冶組)，礦體一般處于碎屑巖和碳酸鹽巖的接觸部位(常印佛等，1991)。

宣城礦集區(qū)位于長(zhǎng)江中下游成礦帶東南側(cè)(圖1)。近年來(lái)，隨著區(qū)域找礦的進(jìn)展，在該礦集區(qū)發(fā)現(xiàn)了許多礦床，如茶亭斑巖銅金礦床、蕎麥山銅多金屬礦床、劉家山硫鐵礦床和麻姑山銅鉬礦床等。這些礦床主要為斑巖-矽卡巖型礦床(Mao et al.,2006;Zhou et al.,2011;李躍等,2019;徐曉春等,2019)，處于長(zhǎng)興-廣德凹斷褶束和寧蕪凹斷褶束2個(gè)Ⅳ級(jí)構(gòu)造的結(jié)合部位，區(qū)內(nèi)構(gòu)造主要表現(xiàn)為褶皺構(gòu)造、推覆體構(gòu)造和斷裂構(gòu)造(錢仕龍等,2017;徐曉春等,2019)，區(qū)內(nèi)中部疊置有北東向的九連山-貍橋、敬亭山、麻姑山推覆體構(gòu)造，亦為復(fù)背斜構(gòu)造。其中，九連山-貍橋推覆體往北東延伸至江蘇境內(nèi)為茅山推覆體(韓克從等,1985;朱志澄,1991;徐學(xué)思等,2001;黃潤(rùn)生等,2010)，斷續(xù)綿延超過(guò)100 km，目前已將宣城礦集區(qū)列為長(zhǎng)江中下游成礦帶的第八大礦集區(qū)(周濤發(fā)等，2017)。

圖1 長(zhǎng)江中下游成礦帶礦集區(qū)分布圖(據(jù)周濤發(fā)等，2017修改)Fig.1 Distribution map of ore concentration areas in the metallogenic belt in the middle and lower reaches of the Yangtze River(modified after Zhou et al.,2017)

蕎麥山銅多金屬礦床是區(qū)內(nèi)典型的矽卡巖型礦床，位于長(zhǎng)山復(fù)背斜的次一級(jí)連續(xù)褶皺構(gòu)造蕎麥山茶山倒轉(zhuǎn)背斜的南西傾伏端。礦床銅金屬量約5.9×104t、硫鐵礦礦石量約5×106t、鎢金屬量8.3×103t，為中型矽卡巖礦床(劉曉明等,2015)。礦區(qū)內(nèi)地層(圖2)由古至新分別為泥盆系上統(tǒng)五通組、石炭系中上統(tǒng)黃龍船山組、二疊系下統(tǒng)棲霞組、二疊系下統(tǒng)孤峰組、二疊系上統(tǒng)龍?zhí)督M及第四系全新統(tǒng)。礦體主要賦存于石炭系黃龍船山組與侵入體接觸帶內(nèi)。礦區(qū)侵入巖巖體平行于褶皺軸沿?cái)嗔褞ИM長(zhǎng)狀分布，呈巖枝、巖株或巖脈狀產(chǎn)出，屬于中-淺成相中酸性侵入巖，主要侵入蕎麥山倒轉(zhuǎn)背斜核部和北西翼部，與圍巖接觸面較陡，傾角50°～80°，有的呈脈狀穿插在石英砂巖和石灰?guī)r裂隙中，巖體為花崗閃長(zhǎng)斑巖，多位于蕎麥山南坡超覆于石英砂巖之上。礦區(qū)內(nèi)脈巖主要為石英斑巖，大致沿孤峰組和棲霞組之間分布，具有明顯的定向性，走向北東70°，厚約20 m，分布于礦區(qū)南部的茨山的東部的頂門山等地。

圖2 蕎麥山礦床地質(zhì)圖(據(jù)劉曉明等，2014修改)1—第四系；2—二疊系上統(tǒng)龍?zhí)督M；3—二疊系下統(tǒng)棲霞組；4—石炭系中上統(tǒng)黃龍組、船山組；5—泥盆系上統(tǒng)五通組；6—花崗閃長(zhǎng)斑巖；7—矽卡巖；8—砂巖；9—花崗斑巖；10—鐵帽；11—勘探線及編號(hào)；12—鉆孔及編號(hào)Fig.2 Geological map of the Qiaomaishan deposit(modified after Liu et al.,2014)1—Quaternary;2—Longtan Formation of Upper Permian;3—Qixia Formation of Lower Permian;4—Huanglong Formation and Chuanshan Formation of Middle and Upper Carboniferous;5—Wutong Formation of Upper Devonian;6—Granodiorite porphyry;7—Skarn;8—Sandstone;9—Granite porphyry;10—Gossan;11—Exploration line and number;12—Drill hole and number

2 數(shù)據(jù)采集

pXRF是最常見的野外巖芯原位元素測(cè)試儀器，其優(yōu)點(diǎn)在于其輕巧、便攜、測(cè)試速度快、精度高，在野外可直接獲取目標(biāo)物的元素含量，當(dāng)前在環(huán)境監(jiān)測(cè)、巖芯測(cè)試、分析沉積物序列等方面有著廣泛的應(yīng)用(Peinado et al.,2010;Yuan et al.,2014;劉琦,2014;Mauriohooho et al.,2016;Al-Musawi et al.,2020)。在合理設(shè)置參數(shù)及模式的基礎(chǔ)上，pXRF技術(shù)被認(rèn)為是一種滿足全數(shù)據(jù)采集要求、便捷有效的測(cè)試方法(Radu et al.,2009;Piercey et al.,2014;Perring et al.,2018;Zhou et al.,2018;Zhou et al.,2020)，如夏慶霖等運(yùn)用pXRF 技術(shù)對(duì)安徽泥河鐵礦進(jìn)行了巖芯礦化蝕變研究，獲得了元素組合異常，較好地指示了礦化(夏慶霖等,2011)；McNulty 等（2018）利 用pXRF 對(duì)加拿大溫哥華島上的邁拉瀑布火山巖礦床進(jìn)行了巖性判別研究，分析礦床Ti/Zr 趨勢(shì)確定火山巖的類型。

基于pXRF 技術(shù)的巖芯測(cè)試可以最大限度保持樣本的完整性及連續(xù)性，可用于野外鉆孔現(xiàn)場(chǎng)原位分析，具有效率高、投入小的特點(diǎn)，大大提高了找礦勘查的效率。本文運(yùn)用pXRF 技術(shù)在野外直接獲取鉆孔巖芯測(cè)試分析數(shù)據(jù)，并基于數(shù)據(jù)驅(qū)動(dòng)思維，利用數(shù)據(jù)挖掘方法分析該礦床的元素空間分布規(guī)律。本研究所選擇的儀器是日本奧林巴斯(OLYMPUS)科技公司生產(chǎn)的VANTA 系列便攜式X 射線熒光分析儀，采取地球化學(xué)模式進(jìn)行測(cè)量，在該模式下可測(cè)量元素包括P、Cl、Mg、Al、Si、Cr、K、Ca、S、Ti、Cu、V、Mn、Fe、As、Co、Ni、W、Bi、Zn、Th、Hg、Sr、Y、Pb、Se、Mo、U、Rb、Cd、Zr、Nb、Ag、Sn、Sb。

本文選取了蕎麥山銅多金屬礦床的代表性剖面16+線上的以銅硫礦石為主的3 個(gè)典型鉆孔，分別為ZK16+02、ZK16+03、ZK16+04(圖3)。在保持巖芯清潔、干燥的情況下，經(jīng)過(guò)等間距單點(diǎn)測(cè)試、校準(zhǔn)方程建立等步驟，獲得了測(cè)量有效的30 種元素Mg、Al、Si、P、S、Cl、Ca、Ti、V、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Rb、Sr、Y、Zr、Nb、Mo、Ag、W、Hg、Pb、Th、U 并進(jìn)行數(shù)據(jù)分析，測(cè)點(diǎn)間隔為1 m 左右，不同鉆孔不同深度位置間隔略有變化，每個(gè)點(diǎn)測(cè)試時(shí)間為60 s，累計(jì)測(cè)點(diǎn)433個(gè)(Bourke et al.,2016;周國(guó)玉,2021)。

圖3 16+線勘探剖面簡(jiǎn)圖(據(jù)劉曉明等，2014)1—石炭系中上統(tǒng)黃龍組、船山組；2—泥盆系上統(tǒng)五通組；3—花崗閃長(zhǎng)斑巖巖體；4—石英砂巖；5—大理巖；6—花崗閃長(zhǎng)斑巖；7—鉆孔及孔號(hào)；8—銅硫礦石；9—銅礦石；10—鎢礦石Fig.3 Simplified cross section of exploration line 16+through the Qiaomaishan Cu deposit(after Liu et al.,2014)1—Huanglong Formation and Chuanshan Formation of Middle and Upper Carboniferous;2—Wutong Formation of Upper Devonian;3—Granodiorite porphyry;4—Quartz sandstone;5—Marble;6—Granodiorite porphyry;7—Drill hole and number;8—Copper-sulfur ore;9—Copper ore;10—Tungsten ore

3 蕎麥山銅多金屬礦床元素富集規(guī)律分析方法研究

3.1 主成分分析

主成分分析法(PCA)對(duì)元素?cái)?shù)據(jù)進(jìn)行降維可以除去一些重復(fù)變量因素的影響，將信息進(jìn)行壓縮，建立盡可能少的新變量，使相關(guān)元素之間的關(guān)系更加明了(成秋明等,2007;韓小孩等,2012)。本文利用SPSS統(tǒng)計(jì)軟件對(duì)蕎麥山銅多金屬礦床16+線所測(cè)的30 種有效元素進(jìn)行全鉆孔主成分分析，最終將30 個(gè)元素壓縮得到8 個(gè)主成分信息，累計(jì)方差貢獻(xiàn)率達(dá)到70.576%，保留了鉆孔中元素含量的主要信息，其中第一主成分累計(jì)方差貢獻(xiàn)率達(dá)到27.079%，所得成分中包含了蕎麥山銅多金屬礦床成礦元素Cu、Fe、S、W 在內(nèi)的14 個(gè)正載荷元素（Cl 元素由于得分過(guò)低，暫不考慮）。主成分分析結(jié)果見表1。

在此基礎(chǔ)上利用全數(shù)據(jù)主成分分析得到的結(jié)果分別計(jì)算ZK16+02、ZK16+03、ZK16+04 深部空間的成分得分值，成分得分結(jié)果的可視化顯示能夠直觀反映各成分所包含的信息與鉆孔巖芯成分之間的對(duì)應(yīng)關(guān)系（圖4）。

蕎麥山銅多金屬礦床16+線主成分分析反映了如下特征（表1，圖4）：

圖4 蕎麥山銅多金屬礦床16+線主成分得分圖Fig.4 Principal component score of exploration line 16+through the Qiaomaishan Cu polymetallic deposit

表1 蕎麥山銅多金屬礦床16+線鉆孔主成分方差及成分載荷表Table 1 Principal component variance and component load of boreholes on exploration line 16+through the Qiaomaishan Cu polymetallic deposit

（1）綜合主成分得分（PCA）包含了全數(shù)據(jù)70.576%的信息在內(nèi)，高值區(qū)主要出現(xiàn)在包括銅硫礦石和鎢礦石在內(nèi)的礦化區(qū)，是對(duì)鉆孔主要信息的集中反映，計(jì)算公式如下：

（2）第一主成分(PC1)的高值富集區(qū)域與綜合主成分高度相似，與鉆孔的礦化區(qū)域明顯對(duì)應(yīng)，顯然PC1代表了富含元素Cu、Fe、S、W在內(nèi)的礦化信息，且伴隨有非成礦元素Mg、Ca、Mn、Co、Zn、As、Se、Ag、Hg、U。

（3）第二主成分(PC2)具有除Mg、Si、Rb、Sr在內(nèi)的較多的正載荷元素，結(jié)合主成分得分圖可以看出第二主成分信息并不具備有類似于第一主成分具有明顯礦化特征的元素組合，但高值區(qū)也多靠近礦石周圍，在偏向石英砂巖的區(qū)域有著較強(qiáng)的反映，總體上具有更為復(fù)雜的表現(xiàn)。

（4）第三主成分（PC3）具有正載荷元素Mg、S、Ti、V、Cr、Mn、Zn、Rb、Sr、Y、Nb、Mo、W、Hg，主成分得分高值區(qū)呈現(xiàn)出以銅硫礦石和鎢礦石為中心向兩側(cè)逐漸降低輻散的趨勢(shì)，結(jié)果表明，PC2 在一定程度上能夠反映銅硫礦石和鎢礦石附近復(fù)雜多樣的圍巖蝕變現(xiàn)象。

（5）第四主成分（PC4）所表現(xiàn)的不具備明顯特征，但總體上表現(xiàn)出主成分得分值由礦石位置向非礦石位置由高到低（ZK16+02、ZK16+03）或由低到高(ZK16+04)的現(xiàn)象。

（6）第五主成分（PC5）表現(xiàn)出得分值由石英砂巖段到銅硫礦石段、鎢礦石段再到花崗閃長(zhǎng)巖段逐漸降低變化的現(xiàn)象，反映出自淺向深得分值由高向低的變化。

（7）第六主成分（PC6）在銅硫礦石所對(duì)應(yīng)的位置表現(xiàn)出低值特征，而在鎢礦石所對(duì)應(yīng)的位置表現(xiàn)出高值特征。

（8）第七主成分（PC7）和PC6特征一致。

（9）第八主成分（PC8）所表現(xiàn)的特征與PC6、PC7 的特征相反，在銅硫礦石所對(duì)應(yīng)的位置表現(xiàn)為高值特征，在鎢礦石所對(duì)應(yīng)的位置表現(xiàn)出低值特征。

3.2 多元逐步線性回歸分析

主成分分析能夠保留諸多變量的主要信息，將30 種有效元素壓縮為8 個(gè)變量，能夠定性的反映元素富集規(guī)律與鉆孔巖芯之間的對(duì)應(yīng)關(guān)系。而多元逐步線性回歸分析方法能夠反映多個(gè)變量對(duì)一個(gè)因素的影響，能夠從大量可供選擇的變量中選取最重要的變量，以保證最后保留在模型中的解釋變量多重共線性降到最低，從而建立“最優(yōu)”的多元線性回歸方程，得到回歸分析的預(yù)測(cè)模型（Closs et al.,1975;付勇等,2009;Nazarpour et al.,2016;Ali et al,.2017），從而定量的反映元素的富集關(guān)系，量化礦化信息。為探究蕎麥山銅多金屬礦床中成礦元素Cu、Fe、S、W 與非成礦元素之間的富集伴隨關(guān)系，筆者利用SPSS 軟件將蕎麥山銅多金屬礦床16+線鉆孔的成礦元素與26 種非成礦元素進(jìn)行多元逐步線性回歸分析，并建立預(yù)測(cè)模型，以達(dá)到進(jìn)一步提取主要相關(guān)元素的目的。

為了避免不相關(guān)元素的干擾，在多元逐步回歸之前對(duì)所測(cè)得的30 種有效元素進(jìn)行相關(guān)性分析，結(jié)果見表2，相關(guān)性分析結(jié)果可知，與蕎麥山銅多金屬礦床與成礦元素Cu、Fe、S、W 具有正相關(guān)性元素主要有Mg、Ca、Mn、Co、As、Se、Ag、Hg、U。而由上述主成分分析表明，第一主成分集中反映了鉆孔礦化信息，并且從主成分載荷表中可以看出元素非成礦元素Mg、Ca、Mn、Co、Zn、As、Se、Ag、Hg、U均具有較高的正載荷系數(shù)，與相關(guān)性分析相比僅多出元素Zn，且載荷系數(shù)最小，因此，相關(guān)性分析與主成分分析所得出的結(jié)論具有高度一致性。

表2 蕎麥山銅多金屬礦床16+線鉆孔成礦元素相關(guān)性分析表Table 2 Correlation analysis of ore-forming elements in boreholes on exploration line 16+through the Qiaomaishan Cu polymetallic deposit

結(jié)合主成分分析和相關(guān)性分析的結(jié)果，可以初步反映出與成礦元素Cu、Fe、S、W 富集伴隨的正相關(guān)元素?；诖耍盟鶞y(cè)433 個(gè)數(shù)據(jù)點(diǎn)對(duì)成礦元素進(jìn)行多元逐步線性回歸，無(wú)差別的將全部鉆孔數(shù)據(jù)步入模型，以分析元素Mg、Ca、Mn、Co、As、Se、Ag、Hg、U 與成礦元素的量化關(guān)系，建立擬合模型，以探究諸多元素的富集成礦規(guī)律，分析結(jié)果見表3。

表3 多元逐步線性回歸分析系數(shù)表Table 3 Multiple stepwise linear regression analysis coefficient table

雖然Zn 元素在主成分分析中反映為低正載荷，但在初步的相關(guān)性分析中（表2），該元素與成礦元素并非均反映為顯著相關(guān)性，因此并未作為變量步入回歸模型。盡管在相關(guān)性分析中，元素Mg、Ca、Mn、Co、As、Se、Ag、Hg、U 對(duì)成礦元素Cu、Fe、S、W 均各自呈現(xiàn)出正相關(guān)關(guān)系，但在多元素共同作用的前提下，元素自身之間也存在相互作用即多重共線性作用，因此逐步回歸模型在回歸過(guò)程中對(duì)變量進(jìn)行邊引入邊剔除，直到既無(wú)新變量的引入也無(wú)舊變量的剔除，以建立“最優(yōu)”的多元線性回歸模型。從分析結(jié)果中可以看出，對(duì)于不同的成礦元素所保留的元素及系數(shù)也有所不同，模型R2系數(shù)為分別為0.589、0.888、0.624、0.340，且各元素的系數(shù)均在滿足顯著性檢驗(yàn)條件，因此，可以認(rèn)為模型具有較好的擬合效果。此外，各模型VIF 均小于5，即擬合模型中不具有多重共線性影響，說(shuō)明對(duì)于各成礦元素所保留下來(lái)的元素之間不存在明顯相互作用，且其共同的富集對(duì)銅多金屬礦床的形成產(chǎn)生著正向作用。逐步線性回歸分析反映出Cu 元素的富集伴隨著元素Mg、Co、As、Se、Ag、Hg、U；Fe 元素的富集伴隨著元素Ca、Mn、Ag、Se；S 元素的富集伴隨著元素Mg、Ca、Mn、As、Se、Ag、Hg；W 元素的富集伴隨著元素Mn、Hg。

此外，通過(guò)上述逐步線性回歸分析系數(shù)表建立多元逐步線性回歸擬合模型，對(duì)蕎麥山銅多金屬礦床ZK16+02、ZK16+03、ZK16+04 鉆孔中成礦元素進(jìn)行擬合，以Cu為例，建立擬合模型如下：

Cu=35.736×Co+23.277×As+61.993×Se+0.012×Mg+681.668×Ag+22.317×U-7.788×Hg-0.0045

各成礦元素模型建立同理，利用擬合結(jié)果對(duì)比鉆孔柱狀圖以及元素實(shí)測(cè)值，以此驗(yàn)證模型的準(zhǔn)確性（圖5）。從對(duì)比結(jié)果中可以看出，成礦元素的預(yù)測(cè)高值區(qū)與銅硫礦石和鎢礦石的位置均具有較好的對(duì)應(yīng)關(guān)系，并且以非成礦元素Mg、Ca、Mn、Co、As、Se、Ag、Hg、U 所建立的預(yù)測(cè)模型與實(shí)測(cè)值也和銅硫礦石和鎢礦石具有高度的吻合現(xiàn)象。由此可見，相關(guān)性分析及多元逐步線性回歸分析將26 種元素提取出了最為關(guān)鍵的相關(guān)元素以量化了元素富集信息，通過(guò)對(duì)成礦元素正相關(guān)元素進(jìn)行多元逐步線性回歸能夠進(jìn)一步提取對(duì)蕎麥山銅多金屬礦床具有正向影響的元素組合，并建立量化模型，定量反映了元素富集規(guī)律以及礦化信息。

圖5 蕎麥銅多金屬礦床16+線成礦元素?cái)M合對(duì)比圖Fig.5 Fitting comparison of metallogenic elements in boreholes on exploration line 16+through the Qiaomaishan Cu polymetallic deposit

4 討論

4.1 區(qū)別于實(shí)驗(yàn)室方法的全數(shù)據(jù)分析

當(dāng)前便攜式X 射線熒光分析技術(shù)的發(fā)展已經(jīng)取得了很大的進(jìn)步，測(cè)試檢出限也提高了很多，但是pXRF 數(shù)據(jù)精度與傳統(tǒng)實(shí)驗(yàn)室分析精度之間仍然有一定的差距，pXRF 分析總體來(lái)看屬于半定量分析。本研究并非脫離于實(shí)驗(yàn)室方法完全基于pXRF 的數(shù)據(jù)采集，為了提高pXRF 數(shù)據(jù)的可信度，以更好地反映研究區(qū)的元素富集分布情況，本研究選擇了以實(shí)驗(yàn)室數(shù)據(jù)為標(biāo)準(zhǔn)來(lái)校準(zhǔn)pXRF 數(shù)據(jù)。在數(shù)據(jù)采集的預(yù)處理階段，選取了蕎麥山銅多金屬礦床采集樣品中13 塊具有代表性的巖芯樣品進(jìn)行實(shí)驗(yàn)室全巖主微量元素測(cè)試，樣品分別位于礦床的不同深度（深度從10～140 m），并且樣品巖性從淺到深包含了蕎麥山銅多金屬礦床4 種主要巖性在內(nèi)的石英砂巖、矽卡巖型礦體、大理巖和花崗山長(zhǎng)斑巖。研究發(fā)現(xiàn)pXRF 分析結(jié)果與實(shí)驗(yàn)室分析結(jié)果之間具有較好的正相關(guān)關(guān)系，對(duì)于不同樣品中元素含量的高值、低值反映結(jié)果是一致的（周國(guó)玉,2021）。在巖芯數(shù)據(jù)采集中，對(duì)蕎麥山銅多金屬礦床16+線主要以銅硫礦石為主的3 個(gè)典型鉆孔進(jìn)行等間距連續(xù)單點(diǎn)測(cè)試，累計(jì)測(cè)點(diǎn)達(dá)到433 個(gè)，各點(diǎn)采集有效元素30 種，在獲取豐富數(shù)據(jù)的前提下最大限度地保持了巖芯的完整性以及測(cè)試樣本的連續(xù)性，這是實(shí)驗(yàn)室測(cè)樣難以做到的。

雖然，pXRF 化學(xué)分析是一種半定量、半定性的分析過(guò)程，但是其具有連續(xù)性測(cè)試、分析密度大、測(cè)驗(yàn)樣品多的特性。并且，pXRF可以在無(wú)需樣品制備的情況下很容易地在現(xiàn)場(chǎng)使用（Pinto,2018），其提供快速、實(shí)時(shí)分析的能力有助于收集適合地統(tǒng)計(jì)分析的有針對(duì)性的高數(shù)據(jù)密度讀數(shù)（Lemiere,2018），盡管不具備實(shí)驗(yàn)室方法的高精度特征，與地球化學(xué)分析值存在本質(zhì)區(qū)別，但能一定程度上能夠反映出某區(qū)域樣品的連續(xù)性特征與相關(guān)性特征，可以獲取更全面的巖性信息，突出相對(duì)變化，相比實(shí)驗(yàn)室測(cè)樣能夠更廣泛便捷地分析巖芯巖性變化趨勢(shì)特征，有助于全面反映礦化特征和蝕變現(xiàn)象，因此并不會(huì)降低用于勘探過(guò)程中元素分布、遷移、富集規(guī)律研究的數(shù)據(jù)價(jià)值。

4.2 元素富集規(guī)律與礦化、蝕變

作者通過(guò)利用pXRF 對(duì)蕎麥山銅多金屬礦床16+線3個(gè)典型鉆孔（ZK16+02、ZK16+03、ZK16+04）進(jìn)行連續(xù)、高密度的主微量元素全數(shù)據(jù)采集，運(yùn)用多元數(shù)據(jù)主成分分析法、多元逐步線性回歸法對(duì)各主微量元素的空間分布規(guī)律、元素異常組合情況和成礦指示元素等進(jìn)行了詳細(xì)分析。

本文對(duì)上述3 個(gè)鉆孔元素進(jìn)行主成分分析（圖4），研究結(jié)果表明，在蕎麥山銅多金屬礦床16+線鉆孔的主成分分析中，將30 種元素信息壓縮為8個(gè)主成分，每個(gè)主成分均代表了鉆孔的不同信息：①以正載荷元素Cu、Fe、S、W、Mg、Ca、Mn、Co、Zn、As、Se、Ag、Hg、U 為主的PC1 元素組合代表著蕎麥山銅多金屬礦床的礦化信息，研究結(jié)果從元素分布的角度細(xì)化了前人的研究；②PC2、PC3 反映了蕎麥山銅多金屬礦床中花崗閃長(zhǎng)斑巖、石英砂巖、大理巖相互接觸帶之間復(fù)雜多樣的蝕變信息及變化趨勢(shì)（劉曉明等,2015;2016）；③PC6、PC7、PC8 主成分信息一定程度上蘊(yùn)含了鎢礦石、銅硫礦石的礦化信息。

因此，通過(guò)主成分分析16+線鉆孔數(shù)據(jù)，將30 種元素壓縮為包括礦化信息、蝕變信息、銅硫礦石、鎢礦石、巖石信息等在內(nèi)的8 個(gè)主成分信息，并且綜合得分結(jié)果能夠包含并解釋70.576%的鉆孔信息，通過(guò)主成分分析確定的，以Cu、Fe、S、W、Mg、Ca、Mn、Co、Zn、As、Se、Ag、Hg、U 為正相關(guān)關(guān)系變化的元素組合，可輔助確定蕎麥山銅多金屬礦床的多樣性礦化。分析表明，礦床復(fù)雜多樣的礦化、蝕變是造成主成分分析中正載荷元素集中出現(xiàn)的主要原因，反過(guò)來(lái)多種類元素不同主成分中出現(xiàn)的正載荷現(xiàn)象也能夠一定程度上反映巖體多樣的礦化、蝕變特征，兩者之間相互佐證。

主成分分析側(cè)重于將諸多變量進(jìn)行壓縮合并、整合歸類，最大程度上保留原有信息，而多元逐步線性回歸側(cè)重于將諸多變量進(jìn)行對(duì)比、剔除、保留，減少變量的同時(shí)，建立最優(yōu)擬合方程，以提取相關(guān)信息。

在對(duì)蕎麥山銅多金屬礦床成礦元素Cu、Fe、S、W 進(jìn)行相關(guān)性分析發(fā)現(xiàn)，對(duì)成礦元素起正相關(guān)性影響的元素包括Mg、Ca、Mn、Co、As、Se、Ag、Hg、U，這一結(jié)果與主成分分析得出的結(jié)論大同小異，因此，可以進(jìn)一步判斷這些元素與成礦元素具有明顯的富集特征，可作為礦化指示元素?；谙嚓P(guān)性分析和主成分分析結(jié)果進(jìn)行多元逐步線性回歸分析并建立擬合模型，對(duì)比驗(yàn)證表明，盡管擬合R2系數(shù)大小不一，但都對(duì)礦化位置具有較好的反映效果。ZK16+02由于具有更為復(fù)雜的地質(zhì)環(huán)境，無(wú)論是實(shí)測(cè)元素值還是擬合值，都綜合反映了銅硫礦石和鎢礦石的整體礦化特征；ZK16+03、ZK16+04 中4 種成礦元素的擬合效果與實(shí)測(cè)值具有高度吻合的特征，并且對(duì)銅硫礦石的反映也十分顯著，甚至Cu 元素?cái)M合值對(duì)銅硫礦石位置的指示效果比Cu 元素的實(shí)測(cè)值更好，盡管W 元素的R2系數(shù)只有0.34，但整體上對(duì)于鎢礦床也具有大致的反映?？梢?，多元回歸分析能夠從諸多元素中提取與成礦元素具有最優(yōu)相關(guān)性的成礦指示性元素，與主成分分析對(duì)比研究共同指示出與成礦元素Cu、Fe、S、W 伴隨富集的元素組合為Mg、Ca、Mn、Co、As、Se、Ag、Hg、U。

筆者主要選擇了研究區(qū)代表性剖面的3 個(gè)鉆孔進(jìn)行研究，雖然數(shù)據(jù)涵蓋了礦床的主要巖層，但鉆孔數(shù)量相對(duì)于整個(gè)礦床來(lái)說(shuō)還是偏少，還有一些礦石類型如銅硫鎢組合礦石、銅鎢礦石、銅礦石未能全部覆蓋。在分析方法上主要選擇了主成分分析、相關(guān)性分析、多元逐步線性回歸分析，如若加入一些數(shù)據(jù)算法分析、機(jī)器學(xué)習(xí)分析，將提升數(shù)據(jù)分析的科學(xué)性及有效性。

5 結(jié)論

（1）蕎麥山銅多金屬鉆孔主成分分析結(jié)果顯示，綜合得分PCA、PA1 集中反映了元素礦化信息；PC3反映了礦石與接觸帶圍巖的蝕變信息；PC5一定程度上反映了巖芯中石英砂巖的成分；PC6、PC7 對(duì)鎢礦石具有一定的指示作用；PC8 對(duì)銅硫礦石具有一定的指示作用。

（2）多元逐步線性回歸分析提取了具有成礦指示作用的相關(guān)富集元素，結(jié)果表明，蕎麥山銅多金屬礦床成礦元素Cu 伴隨著元素Mg、Co、As、Se、Ag、Hg、U 的富集；成礦元素Fe 伴隨著元素Ca、Mn、Ag、Se 的富集；成礦元素S 伴隨著元素Mg、Ca、Mn、As、Se、Ag、Hg 的富集；成礦元素W 伴隨著元素Mn、Hg的富集。

（3）主成分分析和多元逐步線性回歸分析共同反映蕎麥山銅多金屬礦床中元素Mg、Ca、Mn、Co、As、Se、Ag、Hg、U 和成礦元素Fe、S、Cu、W 存在正相關(guān)關(guān)系，并且伴隨成礦元素呈現(xiàn)富集特征，對(duì)蕎麥山銅多金屬礦床銅硫礦石具有礦化指示作用。