胡利芳，趙彥巧，曹芳芳，郭凱旋，張俊濤

(1.河南省有色金屬地質(zhì)礦產(chǎn)局第一地質(zhì)大隊， 河南 鄭州 450016； 2.河南省有色金屬礦產(chǎn)探測工程技術(shù)研究中心， 河南 鄭州 450016)

0 引言

崤山申家窯金礦床位于河南省三門峽市南150°方向，直距30 km，行政區(qū)劃隸屬于陜州區(qū)西張村鎮(zhèn)管轄。作為崤山地區(qū)主要的金礦床之一，申家窯金礦床經(jīng)過多年勘查，現(xiàn)已查明金金屬量 1100 kg，銀金屬量296 t，均達(dá)到中型礦床規(guī)模(王春永等，2021)。前人在申家窯金礦區(qū)持續(xù)開展了多年的勘查和科研工作，新增了資源量，在成礦規(guī)律(陳衍景和富士谷，1992；丁振舉和方金云，1996；王學(xué)貞，2004；王宏運等，2013)、成礦年代學(xué)(朱嘉偉等，1999；魏俊浩等，2003)、成礦構(gòu)造(石銓曾等，1993；朱嘉偉等，2001；徐書奎等，2015)、礦床成因(郭方方和王春永，2020)、元素地球化學(xué)(王振凱，2016)等方面取得了多項研究成果，但關(guān)于該礦床原生暈地球化學(xué)方面的研究還很少。多年來人們對原生暈地球化學(xué)的研究多集中于其分帶特征及對深部找礦預(yù)測的作用(李惠等，2010；賴亞等，2018；鄭廣明等，2019；林成貴等，2020)，但對原生暈組合特征與成礦階段劃分的關(guān)系聯(lián)系起來的研究還比較少，僅有王振凱(2016)根據(jù)黃鐵礦、黃銅礦等礦物特征初步劃分了3個成礦期，但暫未有效地指導(dǎo)地質(zhì)勘查工作。本文在前人研究的基礎(chǔ)上，從申家窯金礦床成礦地質(zhì)背景入手，對礦床地質(zhì)特征深入研究，利用SPSS軟件對礦床原生暈地球化學(xué)數(shù)據(jù)進(jìn)行統(tǒng)計處理，開展元素相關(guān)性分析、聚類分析、因子分析等組合特征研究，結(jié)合以往勘查工作成果，對成礦階段進(jìn)行總結(jié)探討，深化崤山地區(qū)成礦規(guī)律認(rèn)識，以期為申家窯金礦乃至崤山地區(qū)金礦地質(zhì)勘查提供支持。

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

豫西崤山地區(qū)的大地構(gòu)造位置處于華北陸塊南緣，經(jīng)歷了多階段、多期次運動。華北陸塊南緣構(gòu)造帶北部為三門峽—魯山斷裂為界，緊鄰華北陸塊，南部則是欒川斷裂，北秦嶺構(gòu)造帶與其相連(圖1)。結(jié)晶基底和蓋層的雙層結(jié)構(gòu)為其突出特點，二者之間為不整合接觸。蓋層巖系變形相對不復(fù)雜、變質(zhì)微弱，主要為淺層次的脆韌性—脆性變形。結(jié)晶基底主巖系為上太古界太華巖群，而中元古界熊耳群和官道口群為其蓋層。構(gòu)造走向以NWW向為主，同時還是其北側(cè)與三門峽—魯山斷裂、南側(cè)與黑溝—欒川斷裂的分界線，發(fā)育多期推覆構(gòu)造和伸展構(gòu)造。巖漿活動集中于新太古代，其次為中生代燕山期。

圖1 研究區(qū)大地構(gòu)造位置示意圖(據(jù)高帥等，2015)1—第四系；2—元古界熊耳群、官道口群；3—元古界寬坪群；4—古生界二郎坪群；5—元古界秦嶺群；6—華北陸塊基底太華巖群；7—南秦嶺造山帶；8—山麓逆沖斷層；9—欒川逆沖斷層；10—區(qū)域大斷裂；11—拆離斷層；12—地質(zhì)界線；13—研究區(qū)范圍

2 礦區(qū)地質(zhì)

2.1 礦區(qū)地質(zhì)概況

區(qū)內(nèi)地層主要為太古界太華巖群(Arth)和第四系(Q)。太華巖群巖性以斜長角閃巖、斜長角閃片麻巖和黑云斜長片麻巖為主。區(qū)內(nèi)構(gòu)造較為簡單，分為NNW、SN、NE向三組，主要斷裂構(gòu)造為NNW向，產(chǎn)狀230°～265°∠60°～75°，這些構(gòu)造破碎帶控制了區(qū)內(nèi)主要金礦脈的分布。侵入巖主要發(fā)育新太古代變質(zhì)花崗巖類，集中分布于崤山斷隆區(qū)結(jié)晶基底內(nèi)，其次為中生代燕山期侵入巖，巖性主要為石英閃長斑巖，次為輝綠巖脈(圖2)。

圖2 申家窯金礦床地質(zhì)簡圖(據(jù)王春永，2021①)1—第四系；2—太古界太華巖群； 3—中元古代輝綠巖脈；4—中生代白堊紀(jì)石英閃長斑巖；5—斷層及編號；6—礦體及編號；7—構(gòu)造蝕變帶及編號；8—勘探線及編號

2.2 礦體特征

申家窯金礦區(qū)發(fā)育崤09、崤01、S2、S3、S4等5條含礦構(gòu)造帶，走向整體呈NNW、近SN向，主要賦存在蝕變斷裂帶的頂?shù)装?，產(chǎn)狀嚴(yán)格受斷裂構(gòu)造控制。賦存礦脈長870～3000 m，厚10.67～25.95 m，產(chǎn)狀220°～285°∠15°～60°。其中崤01脈規(guī)模最大，其特征簡述如下：

崤01脈長約3500 m，產(chǎn)狀265°∠15°～36°，主要呈NNW向展布，局地近SN向，性質(zhì)為張扭性，局部有膨脹狹縮現(xiàn)象，沿傾斜方向上寬下窄。礦脈具分支復(fù)合現(xiàn)象，呈舒緩波狀平行產(chǎn)出，礦化嚴(yán)格受構(gòu)造帶控制，多在頂?shù)装甯患晒I(yè)礦體。圍巖蝕變主要有硅化、絹云母化、黃鐵礦化、碳酸鹽化，局部見高嶺土化。與金礦化關(guān)系密切的圍巖蝕變?yōu)楣杌杌陧數(shù)装灏l(fā)育，以底板最為強烈，次為黃鐵絹英巖化。圍巖蝕變嚴(yán)格在斷裂帶上下盤發(fā)育，礦化受斷裂構(gòu)造控制。

圖3 申家窯金礦崤01礦體剖面簡圖1—覆蓋層；2—白云石英片巖；3—綠泥片巖；4—斜長角閃巖；5—混合花崗巖；6—黑云斜長片麻巖；7—蝕變巖；8—破碎蝕變巖；9—二長石英片巖；10—石英脈；11—構(gòu)造蝕變帶；12—礦體；13—鉆孔位置及編號；14—產(chǎn)狀

通過勘查在崤01礦脈中發(fā)現(xiàn)兩條礦體，編號分別為崤01-Ⅰ和崤01-Ⅱ。崤01-Ⅰ金銀鉛鋅多金屬礦體是區(qū)內(nèi)規(guī)模最大的工業(yè)礦體，呈NNW—SSE向展布，長1275 m，埋深336.26～525.98 m，賦存標(biāo)高+810～+226 m(圖3)。礦體形態(tài)為似層狀、透鏡狀、脈狀，主要產(chǎn)于構(gòu)造破碎帶中，產(chǎn)狀260°∠15°～36°，總體北緩南陡。崤 01-Ⅱ金銀礦體呈近EW向展布，長約160 m，寬約1050 m，近平行產(chǎn)出于崤01-Ⅰ礦體的下部，上部為金礦體，深部為銀礦體，礦體形態(tài)為細(xì)脈狀，產(chǎn)狀260°∠22°～24°。

2.3 礦石特征

區(qū)內(nèi)礦石結(jié)構(gòu)類型較簡單，以自形、半自形粒狀結(jié)構(gòu)為主。構(gòu)造主要有脈狀、團(tuán)塊狀，以及浸染狀，次為角礫狀、條帶狀等。金礦石中主要金屬礦物以黃鐵礦為主，次為黃銅礦、赤鐵礦，少量方鉛礦、閃鋅礦、褐鐵礦等。非金屬礦物以石英為主，約占25%～49%，次為白云石、方解石，約占15%～20%。銀鉛鋅礦石中主要金屬礦物為方鉛礦、閃鋅礦、黃鐵礦，非金屬礦物主要為石英、長石等。

3 原生暈元素組合特征

原生暈元素組合特征分析包括相關(guān)性分析、聚類分析和因子分析，主要用以揭示元素的分組性質(zhì)，表明元素的聚合規(guī)律，尤其是與金成礦相關(guān)的元素。前人利用SPSS多元統(tǒng)計軟件開展原生暈組合特征研究也取得了一定的成果(蔡朝陽等，2010；戴金旺和樊小磊，2017；耿國帥，2020)。本文通過多元統(tǒng)計分析軟件SPSS對相關(guān)數(shù)據(jù)進(jìn)行了元素組合特征的分析，具體有相關(guān)性、聚類和因子分析，并得到較為簡單明了的數(shù)據(jù)，將其結(jié)果作為地質(zhì)成礦階段劃分的依據(jù)，以便指導(dǎo)礦床的勘查工作。

3.1 樣品整理與分析

樣品采集主要布設(shè)在申家窯金礦床崤01、崤09、S3、S2等4條礦體上的ZKS302、ZKS802、ZKS806、ZK2005 4個鉆孔中，采用鉆探巖心法，采樣點距一般為5.0 m，在礦化層、斷裂構(gòu)造帶部位采樣點距加密為0.5～1.0 m。在全孔中系統(tǒng)采集了307件樣品，所采樣品巖性有二長石英片巖、綠泥片巖、斜長角閃巖、蝕變巖、構(gòu)造蝕變巖等。

據(jù)原地礦部地球化學(xué)探礦工作手冊中確定的各類礦床的地球化學(xué)暈指示元素組合，結(jié)合申家窯金礦體礦化和蝕變特點，選擇了Au、Ag、Cu、Pb、Zn、As、Sb、Sn、W、Bi和Mo等11種元素作為分析測試對象。

3.2 相關(guān)性分析

本文利用SPSS軟件對申家窯礦區(qū)特高值處理后的233件原生暈樣品中的11種地球化學(xué)元素的分析結(jié)果進(jìn)行相關(guān)性分析，結(jié)果如表1所示。

表1 申家窯金礦床原生暈R型因子分析相關(guān)系數(shù)矩陣表

從表1中可以較為直觀地看出，與Au相關(guān)性較好的元素有Cu、Pb、Zn、Sb、Ag和Bi，這幾種元素同Au相關(guān)系數(shù)比較大，說明它們的分布規(guī)律與Au較為相似，也意味著它們能在一定程度上較好的指示Au的富集規(guī)律，可以作為尋找金礦的指示元素；與W、Mo的相關(guān)性也較好，但與Sn幾乎不相關(guān)。Pb和Zn的相關(guān)系數(shù)最大，為0.784，表明二者成礦關(guān)系緊密。Sb和Bi幾乎與各主成礦元素相關(guān)性都較好，指示了礦床成礦時的環(huán)境復(fù)雜。Ag與Pb、Zn的相關(guān)系數(shù)為均大于0.5，表明三者密切的相關(guān)關(guān)系，Ag以類質(zhì)同象形式包含于方鉛礦中，或者以輝鉍銀鉛礦的形式與方鉛礦構(gòu)成固溶體結(jié)構(gòu)，是本區(qū)銀的主要載體礦物(王振凱，2016)。除此之外，Sn與各元素多為負(fù)相關(guān)，或者相關(guān)性弱，則反映出Sn的成礦環(huán)境差異較大。

3.3 聚類分析

本文選擇利用R型聚類分析對元素含量進(jìn)行處理分析。由圖4可見，當(dāng)相似水平距離取1時，Pb和Zn劃為一類，表明其成礦較為緊密，內(nèi)在聯(lián)系較大；當(dāng)相似水平距離取10時，Au、Pb、Zn、Sb歸為一類，Cu、Bi劃為一類；當(dāng)相似水平取15時，Pb、Zn、Au、Sb、Cu、Bi、Ag歸為一類，表現(xiàn)為主成礦元素之間的關(guān)系，其余W、As、Mo、Sn各為一類，與其它元素聚為一類的相關(guān)性水平均較低。

圖4 申家窯金礦床原生暈各元素聚類分析譜系圖

故Pb、Zn、Au、Sb、Cu、Bi、Ag成礦關(guān)系緊密，為主要的成礦指示元素組合。W、As、Mo、Sn各為一類，也反映了礦體的形成是一個多階段成礦的過程?？傮w來說R聚類分析得出的結(jié)果與相關(guān)性分析基本一致。

3.4 因子分析

因子分析在礦床學(xué)、地球化學(xué)等地學(xué)領(lǐng)域得到了成功應(yīng)用(崔來運等，2005；伊燕平等，2012；祁軼宏等，2017)。對233個樣品各元素數(shù)據(jù)進(jìn)行R型因子分析，計算得KMO值為0.803，Bartlett’s球度檢驗給出的相伴概率為0.000，適合于因子分析。

對233個數(shù)據(jù)進(jìn)行方差解釋得出因子的特征值和方差貢獻(xiàn)率如表2。由表2可以看出，由于前兩個因子的特征值和方差貢獻(xiàn)均較大，且二者的累計方差貢獻(xiàn)達(dá)到57.943%，因此前兩個因子的影響是較為重要的。特征值是主成分影響力度大小的主要指標(biāo)之一，按照特征值大于1的原則，選取前3個因子為主因子；從累積方差貢獻(xiàn)率上來看，前3個因子的累積方差貢獻(xiàn)率達(dá)到了68.885%，因此認(rèn)為前3個因子可以較好地反映11種因子的大部分信息。

表2 申家窯金礦床元素方差解

在此基礎(chǔ)上，為了變量特征更明顯，對3個因子采用主成分分析法提取后得到了旋轉(zhuǎn)后的載荷矩陣，同時為了因子含義更為簡明與清楚，利用方差最大法得到了因子旋轉(zhuǎn)(表3)。

表3 申家窯金礦床元素因子荷載表

以相關(guān)系數(shù)0.5為界限，由表3可得出三個主因子，分別為第一主因子(F1)Ag、Zn、Pb、Bi、Au、Cu、Sb，第二主因子(F2)As、Sb、W、Mo，以及第三主因子(F3)Mo。

第一個因子(F1)的方差貢獻(xiàn)率為33.575%，占絕對主要位置，主要元素為Ag、Zn、Pb、Bi、Au、Cu，表明Ag、Au、Pb、Zn礦化密切相關(guān)，Sb和Bi分別代表低溫和高溫元素，但Bi荷載較大，推測為一次中高溫成礦階段。Ag、Zn、Pb、Bi、Au、Cu為親硫元素，是研究區(qū)主要的成礦元素組合，體現(xiàn)了主要的成礦作用，可以作為研究區(qū)主要成礦元素空間聚集的綜合指示。

第二主因子F2的貢獻(xiàn)率為23.279%，主要載荷為As、Sb、W、Mo，為低、高溫元素的融合。As元素上有最高載荷值為0.854，表明低溫環(huán)境持續(xù)時間較長，為主要成礦環(huán)境，形成了As的硫化物毒砂；另外As、Sb為前緣暈元素，易遷移，可推測深部可能的成礦區(qū)域。

第三主因子F3的貢獻(xiàn)率為12.031%，主要載荷為Mo。Mo為高溫指示元素，容易在酸性巖體中富集，為巖漿熱液指示元素。在礦區(qū)北西側(cè)出露面積有0.2 km2的后河花崗斑巖體，礦區(qū)北側(cè)出露龍臥溝花崗巖體，有研究表明兩處巖體形成時代相近，在130 Ma左右，屬燕山期(盧仁等，2013，2014)，該期燕山花崗巖的就位可能與礦床的形成密切相關(guān)，為礦床的形成提供熱源、成礦物質(zhì)等。

4 討論

4.1 成礦階段分析

由相關(guān)性分析可知，Pb和Zn成礦關(guān)系緊密，Cu、Pb、Zn、Sb、Ag、Bi在一定程度上能較好的指示Au的富集規(guī)律，可以作為尋找金礦的指示元素。由聚類分析得知，Pb、Zn、Au、Sb、Cu、Bi、Ag歸為一類，表現(xiàn)為主成礦元素之間的關(guān)系，其中Pb和Zn成礦關(guān)系最為緊密，也表明Au、Ag、Pb、Zn、Cu、Bi、Sb元素組合異常區(qū)是地球化學(xué)找礦的重要標(biāo)志。因子分析可得出成礦元素組合的分帶特征為(Ag、Pb、Zn、Bi-Au-Cu-Sb)-(As-Sb-W-Mo)-Mo。Ag、Pb、Zn、Bi-Au-Cu-Sb為近礦帶，形成了黃鐵礦、黃銅礦、方鉛礦、閃鋅礦、毒砂、輝鉍銀鉛礦等金屬礦物，礦石類型主要為含金多金屬硫化物—砷化物型礦石。同時也形成了含Sb、Bi礦物，Bi主要與方鉛礦伴生；As-Sb-W-Mo為前緣暈元素組合，顯示中高溫成礦的疊加，但以中低溫為主，主要形成了As、Sb的硫化物，例如毒砂和輝鉬礦；Mo為尾暈元素，指示一次成巖活動。

如前文所述，本區(qū)主要礦石礦物為黃鐵礦、黃銅礦、方鉛礦、閃鋅礦、毒砂，脈石礦物主要為石英和菱鐵礦，局部可見少量絹云母。金主要以銀金礦和自然金的形式出現(xiàn)，銀礦物通常和方鉛礦伴生，也見有獨立的輝銀礦產(chǎn)出在黃鐵礦裂隙中。因此，成礦階段大致可以分為三個階段，早期為高溫環(huán)境，形成了大量菱鐵礦，為石英—菱鐵礦階段；中期為石英—多金屬硫化物階段，這一階段成礦作用以自形晶石英、黃鐵礦及黃銅礦、方鉛礦、閃鋅礦等硫化物的大量產(chǎn)出為主要特征，這一階段也是金、銀的主要成礦階段，伴隨含礦熱液的活動，在早期礦物裂隙中形成自然金、銀金礦、輝銀礦、輝鉍銀鉛礦等；晚期為石英—黃鐵礦階段，以細(xì)粒半自形—他形石英、黃鐵礦及部分毒砂為特征。

4.2 深部找礦預(yù)測分析

由原生暈組合特征分析研究可知，申家窯金礦形成環(huán)境復(fù)雜，經(jīng)歷了多期成巖成礦活動。前緣暈指示元素Sb和尾暈指示元素Mo較為活躍(李惠等，1999)，出現(xiàn)在主成礦階段等多個成礦活動中，表現(xiàn)了前緣暈和尾暈疊加的特征，也指示了深部找礦的可能性。

同時利用C.B.格里戈良分帶指數(shù)法對申家窯金礦崤01、崤09兩個礦體開展了分帶特征研究。崤09礦體分帶序列為：As-Ag-Cu-Sn-Mo-Bi-Sb-Pb-Au-Zn-W，作為尾暈元素的Sn、Bi遷移至前暈，而前暈元素Sb則出現(xiàn)在近礦暈部位，表明了成礦的多期疊加作用。崤01礦體為Au-Cu-Zn-Sb-Mo-As-Pb-Sn-W-Ag-Bi，顯示了前緣暈和尾暈疊加。

綜上所述，申家窯金礦床原生暈元素組合特征和分帶特征都表明了成礦元素的疊加特征，顯示了良好的深部找礦前景。

5 結(jié)論

(1)相關(guān)性特征顯示與Au相關(guān)性較好的元素有Cu、Pb、Zn、Sb、Ag和Bi，它們的分布規(guī)律與Au較為相似，能在一定程度上較好的指示Au的富集規(guī)律，可以作為尋找金礦的指示元素。

(2)聚類分析研究把Pb、Zn、Au、Sb、Cu、Bi、Ag歸為一類，表現(xiàn)為主成礦元素之間的關(guān)系，推測為同一成礦期次，是地球化學(xué)找礦的重要標(biāo)志。

(3)因子分析得出成礦元素組合的分帶特征為(Ag、Pb、Zn、Bi-Au-Cu-Sb)-(As-Sb-W-Mo)-Mo，體現(xiàn)了成礦元素疊加的特征。

(4)成礦活動劃分為3個階段，分別為石英—菱鐵礦階段、石英—多金屬硫化物階段、石英—黃鐵礦階段。

(5)礦床原生暈元素組合特征和分帶特征都表明了成礦元素的“反向”—疊加特征，顯示了良好的深部找礦前景。

注 釋

① 王春永，白英建，張磊，馮昂.2021.河南省陜縣申家窯-葫蘆峪礦區(qū)深部金礦普查報告[R].鄭州：河南省有色金屬地質(zhì)礦產(chǎn)局第一地質(zhì)大隊.