西藏知不拉和浪母家果矽卡巖銅礦地質(zhì)特征及驅(qū)龍－知不拉－浪母家果成礦系統(tǒng)的確定

王登芳，蔣宗洋，岳寧飛，路 文，蔣光武

(1.成都理工大學(xué)，四川成都 610059;2.西藏巨龍銅業(yè)有限公司，西藏拉薩 850000)

1 前言

南－北美安第斯斑巖Cu－Mo－Au成礦帶、西南太平洋斑巖Cu－Au成礦帶的找礦和研究工作都已發(fā)現(xiàn)，巨型的斑巖礦床的外圍往往分布著相關(guān)的矽卡巖型礦床，例如:秘魯Cerro de Pasco和Colquijirca斑巖礦床外圍的矽卡巖型Pb－Zn礦和淺成低溫Au－Ag礦(Baumgartner et al.，2009);亞利桑那 Bisbee斑巖Cu－Mo礦外圍有矽卡巖Cu－Au礦產(chǎn)出(Sillitoe，2010);猶他州 Tintic和 Bingham斑巖型 Cu－Mo－Au礦外圍分布一系列矽卡巖型Cu－Au或Cu－Pb－Zn礦(Babcock et al.，1995);西南太平洋島弧背景的印尼Ertsberg－Grasberg斑巖Cu－Au礦外圍的矽卡巖 Cu－Au礦(Meinert et al.，1997)。我國的長江中下游地區(qū)，矽卡巖Cu礦常與斑巖Cu礦同時產(chǎn)在一個礦床內(nèi)，如城門山、武山封三洞銅礦和銅山口Cu－Mo礦(秦克章等，1999;周濤發(fā)等，2008)。

斑巖－矽卡巖系統(tǒng)中相關(guān)聯(lián)的地質(zhì)要素常表現(xiàn)為由斑巖體侵位到碳酸鹽地層形成氧化或還原性矽卡巖(Sillitoe，2010)，如藏東玉龍斑巖銅礦中成礦斑巖體和灰?guī)r地層中發(fā)育廣泛的矽卡巖礦體(唐仁鯉等，1995)。這些成礦系統(tǒng)內(nèi)的矽卡巖型礦床不僅為礦山開采提供高品位的礦體，也為淺部矽卡巖礦床在深部或周圍找尋斑巖型礦床提供優(yōu)勢的研究對象。對這類礦床深入研究不僅有利于理解斑巖型－矽卡巖型成礦體系的形成機理，而且對于找礦勘查具有重要實際意義。

圖1 西藏驅(qū)龍中新世驅(qū)龍－知不拉－浪母家果斑巖－矽卡巖成礦系統(tǒng)地質(zhì)概況(改自西藏自治區(qū)墨竹工卡縣驅(qū)龍礦區(qū)銅多金屬礦勘探報告，2008)Fig.1 Simple geological map of the Miocene Qulong-ZHibula-Langmujiaguo porphyry-skarn Cu－Mo district(modified from geologic exploration report of the Qulong porphyry copper polymetallic deposit，2008)

本文以西藏岡底斯驅(qū)龍巨型斑巖銅鉬礦床南部2 km和東南4 km發(fā)育的知不拉和浪母家果矽卡巖型銅礦為研究對象(圖1)。知不拉礦床已開采多年且已經(jīng)發(fā)現(xiàn)的礦體范圍還在擴大，浪母家果礦床系新近發(fā)現(xiàn)。通過研究兩個礦床的地質(zhì)特征、侵入巖特征與巖石地球化學(xué)及硫化物硫同位素，并結(jié)合驅(qū)龍礦床已有的研究結(jié)果，確立驅(qū)龍－知不拉－浪母家果斑巖－矽卡巖成礦系統(tǒng);結(jié)合前人對斑巖－矽卡巖系統(tǒng)的研究現(xiàn)狀，總結(jié)形成該系統(tǒng)的地質(zhì)條件及勘查意義。

2 知不拉－朗母家果矽卡巖礦床地質(zhì)特征

2.1 礦區(qū)地層和構(gòu)造

知不拉－浪母家果矽卡巖礦床的含礦圍巖為葉巴組凝灰?guī)r－晶屑凝灰?guī)r(圖2)，由安山質(zhì)凝灰熔巖、流紋質(zhì)巖屑凝灰?guī)r、安山質(zhì)晶屑凝灰?guī)r、英安質(zhì)凝灰?guī)r，夾有變石英砂巖和灰?guī)r。矽卡巖體和矽卡巖型礦化都是沿凝灰?guī)r地層分布，受控于層間斷層或巖性(圖2)。而火山巖的蝕變主要為黑云母角巖、長英質(zhì)角巖和斑點角巖(圖2)。其中知不拉礦Cu 品位較富，礦脈產(chǎn)狀(0°～20°∠70°～80°)，已經(jīng)開采多年。受附近區(qū)域性大斷裂以及多期次巖漿活動的影響，次級斷裂構(gòu)造較復(fù)雜，部分地段形成寬5m～20m不等，長達(dá)1000m左右的斷裂破碎帶。知不拉礦段斷裂構(gòu)造總體可分為兩組，即NWW向?qū)娱g滑動斷裂和NW向平推(見圖2)右行走滑斷裂。近東西向?qū)娱g斷裂為成礦期斷裂，北西向斷裂為成礦期后斷裂，對礦體和巖性有短距離錯動(圖2)。

浪母家果地表出露明顯的矽卡巖礦化，產(chǎn)狀355°～25°∠60°～85°，圍巖地層是葉巴組凝灰?guī)r及灰?guī)r地層，灰?guī)r地層發(fā)生明顯的大理巖化(圖1)。主要斷裂分為北西向的層間斷裂和北東向的切穿礦體的成礦后斷裂。

2.2 礦區(qū)侵入巖

圖2 驅(qū)龍外圍知不拉矽卡巖型銅礦床礦區(qū)地質(zhì)簡圖(改自西藏自治區(qū)墨竹工卡縣知不拉礦區(qū)銅多金屬礦資源儲量核實報告，2012)Fig.2 Geological map of the Zhibula skarn copper deposit(modified from reserves verification report of the Zhibula skarn copper polymetallic deposit，2012)

知不拉礦區(qū)地表并未發(fā)現(xiàn)有侵入巖，而在深部主要為花崗閃長巖和花崗細(xì)晶巖兩種巖相，鉆孔中可見到花崗細(xì)晶巖侵位到花崗閃長巖中?；◢忛W長巖呈灰白色，中粗?；◢徑Y(jié)構(gòu)，局部為似斑狀，塊狀構(gòu)造，主要礦物有斜長石(40～45%)、鉀長石(15～20%)、石英(10～15%)、角閃石(7～10%)、黑云母(7～10%);副礦物主要有榍石、磷灰石、鋯石、磁鐵礦及不透明金屬礦物等。巖相學(xué)上與驅(qū)龍礦區(qū)中的花崗閃長巖基本一致，為榮木錯拉巖體的主要組成巖性?；◢徏?xì)晶巖呈白色－粉紅色，細(xì)晶結(jié)構(gòu)，塊狀構(gòu)造，主要礦物為鉀長石(40～45%)、斜長石(15～20%)、石英(20～25%)、黑云母(7～10%);副礦物為磷灰石和鋯石。

圖3 驅(qū)龍外圍知不拉矽卡巖型銅礦床8線勘探線剖面簡圖(改自西藏自治區(qū)墨竹工卡縣知不拉礦區(qū)銅多金屬礦資源儲量核實報告，2012)Fig.3 Geologic cross section 8 from the Zhibula skarn－Cu deposit(modified from reserves verification report of the Zhibula skarn copper polymetallic deposit，2012)

浪母家果礦區(qū)地表出露花崗閃長巖，局部有石英閃長巖(地表工程未完全揭示)，且在鉆孔深部亦發(fā)現(xiàn)有花崗閃長巖、二長花崗斑巖和細(xì)晶巖。目前揭示的巖性為一套復(fù)式侵入體，主要巖性為石英閃長巖－花崗閃長巖－二長花崗斑巖－細(xì)晶巖。石英閃長巖呈灰色－灰黑色，粗粒自形－半自形粒狀結(jié)構(gòu)，塊狀構(gòu)造，主要礦物有斜長石(45～50%)、鉀長石(3～5%)、石英(1～5%)、角閃石(15～20%)、黑云母(15～20%);副礦物主要有磷灰石、鋯石等?；◢忛W長巖呈灰白色，中粗?；◢徑Y(jié)構(gòu)，塊狀構(gòu)造，主要礦物有斜長石(40～45%)、鉀長石(15～20%)、石英(10～15%)、角閃石(7～10%)、黑云母(7～10%);副礦物主要有榍石、磷灰石、鋯石、磁鐵礦及不透明金屬礦物等。與知不拉的花崗閃長巖巖石學(xué)特征一致。二長花崗斑巖巖石呈淺灰－灰白色、淺肉紅色，斑狀結(jié)構(gòu)。斑晶含量約為35～45%，主要有:斜長石，自形－半自形粒狀，粒度(2～6mm)，含量約為15～20%;鉀長石，自形－半自形板柱狀，粒度(4～8mm)，含量約為15%;石英，半自形粒狀，粒度(1～8mm)，含量 10～15%;黑云母，半自形片狀，粒度(1～4mm)，含量5%?；|(zhì):主要為隱晶質(zhì)石英和長石類礦物。巖體中常含有石膏脈體。細(xì)晶巖呈灰白色－粉紅色，細(xì)晶結(jié)構(gòu)，塊狀構(gòu)造，主要礦物為鉀長石(40～45%)、斜長石(15～20%)、石英(20～25%)、黑云母(7～10%);副礦物為磷灰石和鋯石。

圖4 驅(qū)龍外圍浪母家果矽卡巖型銅礦床礦區(qū)地質(zhì)簡圖(改自西藏自治區(qū)墨竹工卡縣知不拉礦區(qū)銅多金屬礦資源儲量核實報告，2012)Fig.4 Geological map of the Langmujiaguo skarn copper deposit(modified from reserves verification report of the Zhibula skarn copper polymetallic deposit，2012)

圖5 驅(qū)龍外圍浪母家果矽卡巖型銅礦床0號線和7號線剖面簡圖(改自西藏自治區(qū)墨竹工卡縣知不拉礦區(qū)銅多金屬礦資源儲量核實報告，2012)Fig.5 Geologic cross section 0 and 7 from the Langmujiaguo skarn－Cu deposit(modified from reserves verification report of the Zhibula skarn copper polymetallic deposit，2012)

2.3 礦體特征

知不拉矽卡巖礦床中矽卡巖占主體，大理巖發(fā)育較少。主要見2條較厚矽卡巖礦體，礦體呈較規(guī)則的板狀(厚8～12m)，向北陡傾(傾角70°～80°)，矽卡巖及礦體傾斜方向指向驅(qū)龍礦床，暗示深部可能與驅(qū)龍相連。圍巖為厚層狀的中－粗粒侏羅紀(jì)葉巴組晶屑凝灰?guī)r，具明顯的流紋構(gòu)造。通過西藏巨龍礦業(yè)公司2012年儲量核查工作(表1)，結(jié)合前人資料及礦化情況確定，礦區(qū)內(nèi)共有四個銅礦體(八層銅礦)。其中Ⅰ、Ⅱ、Ⅲ、Ⅳ號礦體均呈近東西向展布于礦區(qū)中部，呈不規(guī)則似層狀，賦存于中侏羅世葉巴組一巖性段(J2y1)安山質(zhì)晶屑凝灰?guī)r及矽卡巖帶中、在南段Ⅳ號礦體穿越安山質(zhì)晶屑凝灰?guī)r賦存于淺灰色流紋質(zhì)巖屑凝灰熔巖內(nèi)。礦體受近東西向斷裂裂隙控制，產(chǎn)于F1、F2、F3、F6、F8、F9、F10 層間斷層破碎帶中，近東西向?qū)娱g斷層為礦體的儲礦構(gòu)造。礦體產(chǎn)狀嚴(yán)格受近東西向?qū)娱g斷層破碎帶控制。礦體走向一般為284°～104°，傾向 10°～15°，傾角 68°～80°(圖 3)。礦體厚度沿走向和傾向變化較穩(wěn)定，局部具分枝現(xiàn)象。礦體巖性為綠簾石矽卡巖，石榴石矽卡巖，圍巖為淺灰色大理巖化灰?guī)r、安山質(zhì)晶屑凝灰構(gòu)造蝕變巖，礦體與兩側(cè)圍巖界線為漸變關(guān)系。礦體Cu平均品位1.50wt%。各礦體詳細(xì)產(chǎn)狀及主要的礦石礦物特征見表1。在礦體深部的局部地段見有花崗閃長巖與黑云母角巖相接觸(20～48勘探線)，對礦體的完整性沒有明顯影響。知不拉矽卡巖中的金屬礦物主要有黃銅礦、黃鐵礦、斑銅礦及磁鐵礦，以及少量的閃鋅礦、方鉛礦和輝鉬礦(表1)。金屬硫化物與石榴石、石英關(guān)系密切，礦石以塊狀、條帶狀為主，少量的侵染狀。黃銅礦含量遠(yuǎn)大于黃鐵礦，斑銅礦較發(fā)育，磁鐵礦很發(fā)育是知不拉矽卡巖礦床的特點，局部見條帶狀分布的磁鐵礦+石榴石矽卡巖，很少見磁黃鐵礦。

浪母家果矽卡巖礦，位于驅(qū)龍斑巖礦床東南4 km、知不拉礦床以東約2 km處。圍巖為薄層狀的中－細(xì)粒葉巴組凝灰?guī)r和大理巖化灰?guī)r，凝灰?guī)r粒度很細(xì)且晶屑含量少。野外露頭矽卡巖呈順層條帶狀產(chǎn)出，礦區(qū)地表見2條矽卡巖礦體(圖4)。其中Ⅱ號礦體圍巖主要是薄層狀的細(xì)粒凝灰?guī)r，矽卡巖礦體與地層產(chǎn)狀一致，呈互層關(guān)系。地表露頭長(＞1000m)，厚度為20～25m，產(chǎn)狀穩(wěn)定(20°～35°∠65°～85°)。其中Ⅰ號礦體，規(guī)模較小，以磁鐵礦為主，為磁鐵礦－石榴石－透輝石含Cu矽卡巖礦體，寬約6～8m，產(chǎn)狀100°∠65°，產(chǎn)出在大理巖化的厚層狀灰?guī)r地層中(圖5)。該礦體具有明顯的礦化分帶(由西向東):磁鐵礦矽卡巖漸變?yōu)槭袷篙x石矽卡巖。探槽中常見普遍而強烈的孔雀石化，銅礦化主要是與磁鐵礦矽卡巖有關(guān)。浪母家果矽卡巖礦床中的金屬礦物有黃銅礦、黃鐵礦、斑銅礦、磁黃鐵礦及少量的黝銅礦、鏡鐵礦，磁鐵礦極少見。礦石以條帶狀、浸染狀及細(xì)脈狀為主，很少見塊狀礦化。在知不拉矽卡巖中及角巖化的凝灰?guī)r地層中，普遍發(fā)育磁黃鐵礦，呈團斑狀、細(xì)脈狀。

表1 知不拉礦床礦體主要特征表格Table 1 Characteristics of main ore bodies in the Zhibula skarn copper deposit

3 侵入巖巖石地球化學(xué)特征

在野外地表和系統(tǒng)的鉆孔觀察基礎(chǔ)上，選擇知不拉－浪母家果礦區(qū)內(nèi)出露的各種巖石類型進行全巖的地球化學(xué)分析。全巖主量元素數(shù)據(jù)采用順序式X射線熒光光譜儀(XRF－1500)測得;微量元素采用美國Finnian MAT公司生產(chǎn)的ICP－MS Element測試，測試結(jié)果見表2。以上實驗均在中國科學(xué)院地質(zhì)與地球物理研究所完成。

表2 知不拉和浪母家果礦區(qū)主要侵入巖巖石化學(xué)特征Table 2 Petrochemistry of intrusive rocks in the Zhibula and Langmujiaguo skarn copper depoits

續(xù)表2Continued Table 2

分析結(jié)果(表2)顯示兩個矽卡巖礦區(qū)內(nèi)各類侵入巖與驅(qū)龍巨型斑巖型銅－鉬礦中相同巖性的巖石有著一致的巖石地球化學(xué)組成(圖6)。在SiO2－K2O+Na2O圖解上(圖6)，主要落在閃長巖、花崗閃長巖、石英二長閃長巖和花崗巖區(qū)域;在SiO2－K2O圖上(圖6c)，與驅(qū)龍礦區(qū)內(nèi)中新世巖石一致，絕大多數(shù)屬于高鉀鈣堿性系列，少數(shù)為鈣堿性和鉀玄巖的區(qū)域;在ACNK圖解中(圖6d)，這些巖石絕大多數(shù)屬于準(zhǔn)鋁質(zhì)巖石。

中新世侵入巖的稀土配分模式呈典型的“鏟狀”形式，表現(xiàn)為強烈的LREE和HREE的分異，HREE明顯的相對虧損，具有很弱的Eu負(fù)異常，除知不拉的細(xì)晶巖有很強的Eu負(fù)異常。在MORB標(biāo)準(zhǔn)化spider圖解上(圖7)，驅(qū)龍中新世巖石明顯富集大離子親石元素(LILE)和高場強元素(HFSE)，K、Pb的正異常，Ta、Nb、Ti的負(fù)異常。而這些微量元素的特征與驅(qū)龍中新世除石英閃長玢巖的其余巖體特征相似。

4 硫化物S同位素特征

本次工作挑選知不拉礦床中的硫化物開展硫同位素分析測試，分析樣品特征描述見表3。硫化物的分選采用傳統(tǒng)的重液方法，經(jīng)檢查大多數(shù)硫化物的純度大于98%。硫同位素在中國科學(xué)院地質(zhì)與地球物理研究所穩(wěn)定同位素實驗室完成。所用儀器為Delta－S氣體同位素比值質(zhì)譜儀，可以完成C、N、S元素含量與同位素比值的在線分析。儀器靈敏度為 1500mol/ion;離子源真空＜3×10－8mba;分析室真空＜5×10－8mba;90°扇形磁場，Rm=180mm，屬于二級方向聚焦型氣體同位素質(zhì)譜儀。對于SO2測定的內(nèi)部測量精度(n=10)為50μg S時，δ34S的精度為0.20‰。在真空系統(tǒng)和高溫條件下把硫化物轉(zhuǎn)化為純凈的SO2氣體，以備IRMS分析測量，測定其34S與32S的比值。硫同位素數(shù)據(jù)以傳統(tǒng)的per mil的形式表示 δ34S‰，硫同位素的參考標(biāo)準(zhǔn)為?Diable Troilite(CDT)，根據(jù)實驗室的對同一樣品的重復(fù)測試，總的分析誤差小于0.20‰。

圖6 知不拉－浪母家果礦區(qū)侵入體主量元素地球化學(xué)特征(驅(qū)龍巖石化學(xué)數(shù)據(jù)來自楊志明，2008;肖波，2011)Fig.6 Major element plots of the magmatic rocks associated with different ore deposit types(Qulong data from Yang et al.，2008;Xiaobo，2011)

結(jié)果顯示，石榴石矽卡巖中浸染狀和與磁鐵礦共生的黃銅礦、斑銅礦的 δ34S值分別為－4.2～－3.4‰和－3.7～－2.1‰，較為集中;脈狀礦石中，早期含銅脈體中S同位素黃鐵礦和黃銅礦δ34S值相似，為－5.6‰和－5.5‰，晚期多金屬硫化物脈體中黃鐵礦和黃銅礦δ34S值相似，為－0.1‰和－1.6‰，方鉛礦為－0.1‰。由結(jié)果看出，早期有磁鐵礦礦化的矽卡巖(1607－110)和晚期多金屬礦化階段(2801－146)中的硫化物硫同位素并不平衡。前人對驅(qū)龍礦區(qū)中與花崗閃長巖接觸的矽卡巖中黃銅礦的硫同位素(－6.2‰)和晚期不平衡的脈狀礦石(黃銅礦－0.3‰和黃鐵礦－2.4%)相比較，也發(fā)現(xiàn)有類似特征。

同時結(jié)果顯示，知不拉矽卡巖礦床中的硫同位素與驅(qū)龍斑巖銅鉬礦的硫同位素有明顯的差異(圖8)。驅(qū)龍礦床硫同位素變化較小且基本平衡，處于巖漿硫的范圍內(nèi)(－3.0～+3.0‰)。而知不拉礦床的硫同位素值低于驅(qū)龍礦床硫化物與成礦斑巖體。這樣的差異可能是由于成礦熱液氧逸度的變化造成的(Ohmoto，1972)，即當(dāng)矽卡巖形成過程中，由于熱液系統(tǒng)氧逸度升高(磁鐵礦沉淀)，熱液硫以為主要種屬，此時沉淀出的硫化物的δ34S將大大低于系統(tǒng)的總硫值(近似以成礦斑巖值為總硫值)，繼而出現(xiàn)整體δ34S值降低的特征。而在晚期高δ34S值的天水加入可能是造成矽卡巖系統(tǒng)晚期多金屬硫化物階段同位素不平衡和δ34S值升高的重要因素。

5 討論

5.1 驅(qū)龍－知不拉－浪母家果斑巖－矽卡巖成礦系統(tǒng)的確定

斑巖型Cu±Mo±Au的礦化中心的外圍熱液蝕變范圍很大，且在其周圍、廣泛分布著矽卡巖型、淺成熱液型(高硫和低硫型)和脈狀礦床(Sillitoe，2010)。因此研究者往往需要通過厘定斑巖型礦床與外圍礦床之間是否存在時空要素間的耦合特征，查明各個礦床之間是否有著相互關(guān)聯(lián)的地質(zhì)因素，研究各礦床成礦物質(zhì)是否具有同源性和連續(xù)演化的特點，才得以確定各個礦床之間是否具有成因上的聯(lián)系，進而構(gòu)成斑巖成礦系統(tǒng)。

圖7 知不拉－浪母家果礦區(qū)侵入巖微量元素和稀土元素標(biāo)準(zhǔn)化圖解(標(biāo)準(zhǔn)化數(shù)據(jù)來自Sun et al.，1989)Fig.7 Rock/primitive mantle-normalized spider plots and rock/chondrite-normalized REE plot for the magmatic rocks associated with different ore deposit types(after Sun et al.，1989)

斑巖成礦系統(tǒng)中成礦物質(zhì)往往具有同源性和連續(xù)演化的特點。如在斑巖鉬礦和外圍脈狀鉛鋅銀礦所構(gòu)成的礦田中，硫化物Pb－S同位素特征表征系統(tǒng)內(nèi)各礦床成礦物質(zhì)為巖漿成因的特點(Warehem et al.，1998)。從巖體內(nèi)向外，成礦流體的溫度、鹽度和H－O同位素都指示連續(xù)降低的過程(葉會壽，2006)。若在礦床尺度范圍內(nèi)，可能出現(xiàn)有相似鹽度的流體由鉬礦化到鉛鋅礦化溫度降低的特點(Lawley et al.，2010)。

結(jié)合現(xiàn)階段上述對斑巖成礦系統(tǒng)的研究，我們確定驅(qū)龍－知不拉－浪母家果斑巖－矽卡巖成礦系統(tǒng)主要基于如下幾點:

1)時代相近，空間毗鄰:驅(qū)龍礦床的成巖時代為17.9 Ma～13.1 Ma，(孟祥金等，2003;芮宗瑤等，2003;Yang et al.，2009);輝鉬礦 Re－Os年齡為16.4±0.5 Ma(孟祥金等，2003)，而知不拉矽卡巖礦床的輝鉬礦年齡為 16.9±0.6 Ma(李光明等，2005)。兩者在誤差范圍內(nèi)一致，且兩礦床在空間上如此相近(與驅(qū)龍礦床相距約2～4 km)，說明兩者在時空上隸屬于同一成礦系統(tǒng)。

2)工程揭示驅(qū)龍斑巖銅礦花崗閃長巖接觸帶發(fā)育有明顯矽卡巖礦化:驅(qū)龍南部為花崗閃長巖與葉巴組凝灰?guī)r地層直接接觸(圖1)。雖然知不拉與浪母家果矽卡巖礦床地表與驅(qū)龍中新世侵入雜巖體有一定距離，但是在驅(qū)龍最南部鉆孔 ZK820和ZK442揭示花崗閃長巖與接觸帶附近都發(fā)育有典型的矽卡巖化，矽卡巖中有少量金屬硫化物。這說明該成礦系統(tǒng)的矽卡巖化是由驅(qū)龍中新世侵入巖的侵位造成的。

3)與礦化有關(guān)的侵入巖地球化學(xué)特征相似:知不拉和浪母家果礦床深部出現(xiàn)和驅(qū)龍礦床相似的花崗閃長巖、二長花崗斑巖和細(xì)晶巖(圖1)。通過本次研究并對比前人對驅(qū)龍礦床成礦巖石的地球化學(xué)研究，發(fā)現(xiàn)三個礦床中的巖漿巖地球化學(xué)特征相近，說明三者成礦巖體屬于同一中新世多期侵位的巖漿系統(tǒng)。

4)成礦物質(zhì)演化趨勢相似:硫同位素的研究揭示驅(qū)龍礦區(qū)中花崗閃長巖接觸的矽卡巖中硫化物的硫同位素值與變化趨勢，和知不拉矽卡巖礦床整體的硫同位素變化趨勢及晚期脈狀礦體硫化物不平衡的特征相一致。

表3 西藏知不拉矽卡巖礦床中的硫化物的硫同位素Table 3 S isotopic analysis result for the Zhibula skarn copper deposit

圖8 驅(qū)龍成礦系統(tǒng)硫化物硫同位素圖解Fig.8 S isotopic diagram for the Qulong porphyryskarn Cu－Mo mineralization system

5)矽卡巖連續(xù)演化的物理化學(xué)條件:驅(qū)龍成礦系統(tǒng)內(nèi)的知不拉矽卡巖中普遍出現(xiàn)磁鐵礦，而浪姆家果矽卡巖中則大量出現(xiàn)磁黃鐵礦。這正好體現(xiàn)出斑巖－矽卡巖成礦系統(tǒng)中，從接觸帶由內(nèi)向外，矽卡巖成礦體系由氧化狀態(tài)遞變?yōu)檫€原狀態(tài)(肖波等，2011)。

5.2 斑巖－矽卡巖成礦系統(tǒng)的特征與形成條件

斑巖－矽卡巖系統(tǒng)中兩類礦床的蝕變礦化特征有著明顯差別。斑巖型礦床以其典型的環(huán)狀蝕變?yōu)樘卣鳎詭r體中心往外可依次發(fā)育鉀化蝕變帶、絹云母蝕變帶、粘土化蝕變帶(這兩者為長石分解蝕變帶)和青磐巖化蝕變帶。許多大型－巨型斑巖型礦床蝕變帶更為復(fù)雜(如美國Henderson斑巖鉬礦中出現(xiàn)黃玉蝕變，Seedorff et al，2004;驅(qū)龍斑巖銅礦中的深部鈉鈣蝕變，肖波，2011)，但礦體主要集中在鉀化蝕變帶和絹云母蝕變帶中。而與斑巖有關(guān)的矽卡巖型礦床通常產(chǎn)于成礦巖體與碳酸鹽巖圍巖接觸帶的矽卡巖及其附近的交代巖中，形成矽卡巖型礦化，且礦化多產(chǎn)于退變矽卡巖中(Einaudi，1982)。這樣的斑巖－矽卡巖成礦系統(tǒng)常發(fā)育明顯的礦化分帶，表現(xiàn)為隨著逐漸遠(yuǎn)離巖體，近端矽卡巖常表現(xiàn)為Cu(Mo)礦，遠(yuǎn)端矽卡巖常表現(xiàn)為Cu－Pb－Zn礦，遠(yuǎn)端矽卡巖礦體最遠(yuǎn)距離可達(dá)5km(如美國的Bingham礦床，Atkinson et al.，1978)，在斑巖體外圍巖中常發(fā)育脈狀鉛鋅銀等礦化，并構(gòu)成Pb－Zn－Ag礦體。礦化特征上，距斑巖體由近至遠(yuǎn)，矽卡巖礦化中以黃銅礦占主導(dǎo)(常多于斑銅礦)變?yōu)辄S鐵礦占主導(dǎo)，較少斑銅礦產(chǎn)出，出現(xiàn)閃鋅礦和方鉛礦(Einaudi，1982)。

形成斑巖－矽卡巖型成礦系統(tǒng)的條件主要有:

1)巖漿巖侵位的地層中有碳酸鹽類或鈣質(zhì)地層，即圍巖特征控制。李靖輝(2008)對大別山北麓的斑巖型鉬礦進行研究后發(fā)現(xiàn)斑巖型鉬礦與矽卡巖型鉬銅礦、脈狀Pb－Zn礦具有密切的空間分布特征，在巖體或巖體與圍巖接觸帶，形成受巖體控制的斑巖型鉬礦，在遠(yuǎn)離巖體的非碳酸鹽圍巖地段形成脈狀(細(xì)脈浸染型、大脈型鉬礦)，或在遠(yuǎn)離巖體的碳酸鹽圍巖地段形成矽卡巖型鎢(鉬、銅)礦，組成了斑巖型鉬(銅)礦－矽卡巖型鎢(鉬、銅)礦－脈狀鉬(錸)礦－含銅、鉛鋅石英脈的成礦系列。在驅(qū)龍成礦系統(tǒng)中，圍巖中侏羅統(tǒng)葉巴組地層中的灰?guī)r、泥灰?guī)r夾層為該系統(tǒng)形成矽卡巖礦床提供了有力的圍巖條件。

2)巖漿多期次侵位且侵位較淺。多期次巖漿活動帶來大量成礦熱液。如中新世(18～13 Ma)驅(qū)龍－知不拉－浪母家果斑巖系統(tǒng)的中酸性復(fù)式雜巖體，主要包含:成礦早期的石英閃長巖、花崗閃長巖和黑云母二長花崗巖;主成礦期為二長花崗斑巖;成礦晚期的花崗閃長斑巖和細(xì)晶巖;成礦后的石英閃長玢巖。江西城門山同源兩期六次巖體侵位形成斑巖－矽卡巖礦化(趙一鳴等，1990)。Qin et al.(1998)根據(jù)矽卡巖礦床形成深度將矽卡巖分為頂部矽卡巖、淺部矽卡巖和深部矽卡巖，并認(rèn)為頂部矽卡巖向深部可能出現(xiàn)斑巖型礦體，淺部矽卡巖礦體往往與斑巖型礦體相伴，而深部矽卡巖礦床則不會與斑巖型礦床伴生。如長江中下游地區(qū)目前的勘探顯示，矽卡巖的深部或者周緣總會發(fā)育中酸性含礦斑巖體(Song et al.，2013)。這樣的頂部和淺部矽卡巖型礦化可以作為尋找深部斑巖型礦床的標(biāo)志之一。

3)區(qū)域內(nèi)構(gòu)造活動及礦區(qū)內(nèi)斷裂、裂隙發(fā)育，即構(gòu)造因素控制。在美國Bingham斑巖Cu－Mo－Au礦、矽卡巖Cu礦床及周邊脈狀或交代成因的Pb－Zn－Ag礦由北－北東向的一系列陡傾角斷裂控制，這些北東向的斷裂為成礦巖體的侵入和礦質(zhì)運移－就位提供空間(Redmond et al，2010)。Bingham斑巖礦床本身廣泛發(fā)育的裂隙為高品位礦體的形成提供就位空間(Gruen et al.，2010)。而在美國的Twin Buttes礦床中矽卡巖的蝕變－礦化和成礦巖體侵入主要與北西西向的斷層控制有關(guān)(Barter et al.，1982)。江西城門山銅礦中北西向和北東向構(gòu)造交匯部位限定礦體和成礦巖體的侵入(趙一鳴等，1990)。

印度－歐亞大陸碰撞強烈影響驅(qū)龍礦區(qū)范圍內(nèi)及外圍的侏羅紀(jì)葉巴組地層呈現(xiàn)明顯的北西西向褶皺構(gòu)造，以及礦區(qū)內(nèi)韌性斷裂(北西西和近東西向)與脆性斷層(近東西、北西西和北西向)。而驅(qū)龍礦區(qū)巖石中破裂裂隙范圍大、發(fā)育程度好(在50條/m以上)，并且破裂裂隙率與Cu、Mo元素的品位成正比(肖波等，2008)。這些構(gòu)造因素為驅(qū)龍斑巖－矽卡巖成礦系統(tǒng)的形成提供巖漿侵位和流體運移通道。

5.3 斑巖－矽卡巖成礦系統(tǒng)勘查意義

成礦系統(tǒng)的理念為礦集區(qū)的找礦突破提供有力支撐。對于斑巖－矽卡巖系統(tǒng)的找礦過程中，往往能夠利用矽卡巖礦床找尋到斑巖礦床。例如印度尼西亞的Ertsberg－Grasberg礦床，早期 Ertsberg有數(shù)個Cu－Au－Ag礦體，并按照矽卡巖型礦床開采，而后發(fā)現(xiàn)Grasberg斑巖礦化(Jowitt et al.，2013);澳大利亞的Cadia礦床中，利用早期的矽卡巖Cu－Au－Ag礦化找到深部的斑巖Au－Cu礦化(Foster et al.，2004)。而對與斑巖－矽卡巖成礦系統(tǒng)的確定能夠真實反映礦區(qū)的儲量和礦化類型，如在巴布亞新幾內(nèi)亞Ok Tedi Cu－Au－Ag礦床目前普遍認(rèn)為是斑巖銅礦，但實際上其大部分礦體在矽卡巖中(Bamford et al.，1972;England et al.，1991)。

在西藏岡底斯地區(qū)發(fā)育有大量奧陶－志留紀(jì)、侏羅紀(jì)－白堊紀(jì)的海相碳酸鹽地層(Pan et al.，2012)，中新生代巖漿活動強烈(紀(jì)偉強，2010)，加之歐亞大陸的碰撞作用形成大量的斷裂，這些都為形成斑巖－矽卡巖的復(fù)合成礦系統(tǒng)提供了有利的條件。近期的勘查和研究工作，如努日斑巖－矽卡巖銅鉬鎢礦床的確認(rèn)(陳雷等，2012)、甲瑪矽卡巖－斑巖銅多金屬成礦系統(tǒng)的研究(唐菊興等，2011)和驅(qū)龍－知不拉－浪母家果斑巖－矽卡巖成礦系統(tǒng)的確定，都已揭示這些都說明西藏岡底斯地區(qū)具有尋找斑巖－矽卡巖型礦床的成礦潛力。勘查研究工作應(yīng)著力于在圍巖為碳酸鹽類巖石，已發(fā)現(xiàn)有矽卡巖或斑巖型礦床的周邊是尋找斑巖型或矽卡巖型礦床有利地段。

6 結(jié)論

(1)知不拉與浪母家果矽卡巖礦床中的侵入巖地球化學(xué)特征與驅(qū)龍礦區(qū)的侵入巖地球化學(xué)特征相似，屬于準(zhǔn)鋁質(zhì)、鈣堿性－高鉀鈣堿性巖石;稀土元素表現(xiàn)為強烈的LREE和HREE的分異，除知不拉的細(xì)晶巖有很強的Eu負(fù)異常，其他具有很弱的Eu負(fù)異常;明顯富集大離子親石元素(LILE)和高場強元素(HFSE)，K、Pb 的正異常，Ta、Nb、Ti的負(fù)異常。

(2)知不拉硫化物硫同位素與驅(qū)龍斑巖銅鉬礦的硫同位素相比有明顯的差異，石榴石矽卡巖中浸染狀和與磁鐵礦共生的黃銅礦、斑銅礦的δ34S值分別為－4.2～－3.4‰和－3.7～－2.1‰，較為集中;脈狀礦石中，早期含銅脈體中S同位素黃鐵礦和黃銅礦δ34S值相似，為－5.6‰和－5.5‰，晚期多金屬硫化物脈體中黃鐵礦和黃銅礦δ34S值相似，為－0.1‰和－1.6‰，方鉛礦為－0.1‰。這樣的差異可能是由于成礦熱液氧逸度的變化和晚期高δ34S值的天水加入造成的。

(3)驅(qū)龍－知不拉－浪母家果斑巖－矽卡巖成礦系統(tǒng)的確立，基于如下五點:時代相近，空間毗鄰;工程揭示驅(qū)龍斑巖銅礦花崗閃長巖接觸帶發(fā)育有明顯矽卡巖礦化;與礦化有關(guān)的侵入巖地球化學(xué)特征相似;成礦物質(zhì)演化趨勢相似;矽卡巖連續(xù)演化的物理化學(xué)條件。

(4)鈣質(zhì)圍巖、多期次侵位且侵位較淺的巖漿、與區(qū)域內(nèi)構(gòu)造活動及礦區(qū)內(nèi)斷裂、裂隙發(fā)育，共同決定了斑巖－矽卡巖系統(tǒng)的產(chǎn)出。

致謝中國科學(xué)院地質(zhì)與地球物理研究院李光明副研究員、秦克章研究員、趙俊興博士、西藏巨龍銅業(yè)有限公司夏代祥教授級高級工程師參與了野外調(diào)查研究和樣品分析測試并給予悉心指導(dǎo)，馬進全高級工程師、鄭前明工程師參加了野外工作，在此一并致以衷心感謝!

Atkinson，W.，Jr.，Einaudi，M.T.1978.Skarn formation and mineralization the contact in aureole Carr Fork，Bingham，Utah［J］.Economic Geology，73:1326－1365

Babcock，R.C.，Jr.，Ballantyne，G.H.，Phillips，C.H.1995.Summary of the geology of the Bingham district［C］.Arizona Geological Society Digest，20:316－335

Bamford R.W..1972.The Mount Fubilan(Ok Tedi)porphyry copper deposit，Territory of Papua and New Guinea［J］.Economic Geology，67:1019－1033

Barter，C.F.，Kelly，J.L.1982.Geology of the Twin Buttes mineral deposit，Pima mining district，Pima County，Arizona，in Titley，S.R.，ed.，Advances in geology of the porphyry copper deposits，southwestern North America［M］.Tucson:University of Arizona Press:407－432

Baumgartner， R.， Fontboté， L.， Spikings， R.， Ovtcharova， M.，Schaltegger，U.，Schneider，J.，Page，L.，Gutjahr，M..2009.Bracketing the age of magmatic-hydrothermal activity at the Cerro de Pasco epithermal polymetallic deposit，central Peru:A U－Pb and40Ar/39Ar study［J］.Economic Geology，104:479－504

Chen Lei，Qin Ke-zhang，Li Guang-ming，Li Jin-xiang，Xiao Bo，Jiang Hua-zhai，Zhao Jun-xing，F(xiàn)an Xin，Jiang Shan-yuan.2012.Geological and skarn mineral characteristics of Nuri Cu－W－Mo deposit in southeast Gangdese，Tibet［J］.Mineral Deposits，31(3):417－437(in Chinese with English abstract)

Einaudi，M.T.1982.Description of skarns associated with porphyry copper plutons，in Titley，S.R.，ed.，Advances in Geology of the Porphyry Copper Deposits［M］.Tuscon:University of Arizona Press:139－183

England J.K.，Kilgour I.，Kanau J.L..1991.Processing copper-gold ore at Ok Tedi:Australasian Institute of Mining and Metallurgy，PNG Geology，Exploration and Mining Conference，Rabaul，Papua New Guinea［C］.Proceedings:183－190

Forster D.B.，Seccombe P.K.，Phillips D..2004.Controls on skarn mineralization and alteration at the Cadia deposits，New South Wales，Australia［J］.Economic Geology，99:761－788

Gruen，G.，Heinrich，C.A.，Schroeder，K.2010.The Bingham Canyon Porphyry Cu－Mo－Au Deposit.II.Vein Geometry and Ore Shell Formation by Pressure-Driven Rock Extension［J］.Economic Geology，105:69－90

Ji Wei-qiang.2010.Geochronology and Petrogenesis of Granitic Rocks from East Segment of the Gangdese Batholith，Southern Tibet［D］.Institute of Geology and Geophysics，Chinese Academy of Sciences.Ph.D.thesis:141

Jowitt，S.M.，Mudd，G.M.，Weng，Z..2013.Hidden mineral deposits in cu-dominated porphyry-skarn systems:How resource reporting can occlude important mineralization types within mining camps［J］.E-conomic Geology，108:1185－1193

Lawley，C.J.M.，Richards，J.P.，Anderson，R.G.，Creaser，R.A.，Heaman，L.M..2010.Geochronology and Geochemistry of the MAX Porphyry Mo Deposit and its Relationship to Pb－Zn－Ag Mineralization，Kootenay Arc，Southeastern British Columbia［J］.Canada.E-conomic Geology，105:1113－1142

Li JH.2008.Metallogenic Systems and Regularity of Porphyry-Type Molybdenum Deposit in Dabieshan Mountain［J］.Journal of East China Institute of Technology，31(1):25－30

Li Guang-ming，Rui Zong-yao，Wang Gao-ming，Lin Fang-cheng，Liu Bo，She Hong-quan，F(xiàn)eng Cheng-you，Qu Wen-jun.2005.Molybdenite Re－Os dating of Jiama and Zhibula polymetallic copper deposits in Gangdese metallogenic belt of Tibet and its significance［J］.Mineral Deposits，24:481－489(in Chinese with English abstract)

Meinert，L.D.，Hefton，K.K.，Mayes，D.，Tasiran，I.1997.Geology，zonation，and fluid evolution of the Big Gossan Cu－Au skarn deposit，Ertsberg district，Irian Jaya［J］.Economic Geology，92:509－534

Meng Xiang-jin，Hou Zeng-qian，Gao Yong-feng，Huang Wei，Qu Xiaoming，Qu Wen-jun.2003.Re－Os dating for molybdenite from Qulong porphyry copper deposit in Gangdese metallogenic belt，Xizang and its metallogenic significance［J］.Geological Review，49:660－666(in Chinese with English abstract)

Meng Xiang-jin，Hou Zeng-qian，Li Zhen-qing.2006.Sulfur and lead isotope compositions of the Qulong porphyry copper deposit，Tibet:implications for the sources of plutons and metals in the deposit［J］.Acta Geologica Sinica，80:554－560(in Chinese with English abstract)

Ohmoto H.1972.Systematics of sulfur and carbon isotopes in hydrothermal ore deposits［J］.Economic Geology，67:551－578

Pan G T，Wang L Q，Li R S，Yuan S H，Ji W H，Yin F G，Zhang W P，Wang B D.2012.Tectonic evolution of the Qinghai-Tibet Plateau［J］.Journal of Asian Earth Sciences，53:3－14

Qin K Z，Ishihara S..1998.On the possibility of porphyry copper mineralization Japanese Islands［J］.International Geology Review，40:539－551

Qin Ke-zhang，Wang Dong-po，Wang Zhi-tian，Sun Shu.1999.Types，geological background，metallogenic provinces and ore-forming systematics of major copper deposits in eastern China［J］.Mineral Deposits，18(4):359－371(in Chinese with English abstract)

Qin Ke-zhang，Xia Dai-xiang，Li Guang-ming，Xiao Bo，Duo Ji，Jiang Guang-wu，Zhao Jun-xing.2014.Qulong porphyry-skarn Cu－Mo deposit，Tibet［M］.Beijing:Science Press:1－316((in Chinese)

Redmond，P.B.，Einaudi，M.T.2010.The Bingham Canyon Porphyry Cu－Mo－Au Deposit.I.Sequence of Intrusions，Vein Formation，and Sulfide Deposition［J］.Economic Geology，105:43－68

Rui Zong-yao，Hou Zeng-qian，Qu Xiao-ming，Zhang Li-sheng，Wang Long-sheng，Liu Yu-lin.2003.Metallogenic epoch of the Gangdese porphyry copper belt and uplift of the Qinghai－Tibet plateau［J］.Mineral Deposits，22:224－232(in Chinese with English abstract)

Seedorff，E.，Einaudi，M.T.2004.Henderson porphyry molybdenum system，Colorado:I.Sequence and abundance of hydrothermal mineral assemblages，flow paths of evolving fluids，and evolutionary style［J］.Economic Geology and the Bulletin of the Society of Economic Geologists，99:3－37

She Hong-quan，F(xiàn)eng Cheng-you，Zhang De-quan，Pan Gui-tang，Li Guang-ming.2005.Characteristics and metallogenic potential of skarn copper-lead-zinc polymetallic deposits in central eastern Gangdese［J］.Mineral Deposits，25:508－520(in Chinese with English abstract)

Sillitoe R H.2010.Porphyry Copper Systems［J］.Economic Geology，105:3－41

Song GX，Qin KZ，Li GM，Li XH，Qu WJ.2013.Geochronology and ore-forming fluids of the Baizhangyan W－Mo deposit in the Chizhou Area，Middle-Lower Yangtze Valley，SE-China［J］.Resource Geology，63(1):57－71

Sun S S.，McDonough W F.1989.Chemical and isotopic systematics of ocean basalts:implications for mantle composition and processes.In:Saunders，AD，Norry MJ(Eds.)，Magmatism in the Ocean Basins.Special Publications［M］.Geological Society，London，42:313－345

Tang Ju-xing，Deng Shi-lin，Zheng Wen-bao，Ying Li-juan，Wang Xiong-wu，Zhong Kang-hui，Qin Shi-peng，Ding Feng，Li Feng-ji，Tang Xiao-qian，Zhong Yu-feng，Peng Hhui-juan.2011.An exploration model for Jiama copper polymetallic deposit in Maizhokunggar County，Tibet［J］.Mineral Deposits，32(2):179－196(in Chinese with English abstract)

Tang Ren-li，Luo Huai-song.1995.The geology of Yulong porphyry copper(molybdenum)ore belt，Xizang(Tibet)［M］.Beijing:Geological Publishing House:1－320(in Chinese with English abstract)

Tibet Julong Copper Industry Co.，LTD.2008.Geologic exploration report of the Qulong porphyry copper polymetallic deposit［R］

Tibet Julong Copper Industry Co.，LTD.2012.Reserves verification report of the Zhibula skarn copper polymetallic deposit［R］

Wareham，C.，Rice，C.，Boyce，A.and Rogers，G..1998.S，C，Sr，and Pb sources in the Pliocene Silver Creek porphyry，Mo system Rio，Colorado［J］.Economic Geology，93:32－46

Xiao Bo，Qin Ke-zhang，Li Guang-ming，Li Jin-xiang，Chen Lei，Zhao Jun-xing，F(xiàn)an Xin.2011.Distribution and characters of Zhibula-Langmujiaguo skarn Cu deposits environing the Qulong porphyry Cu－Mo deposit and their implications for ore-search towards the deep subsurface［J］.Geology and Exploration，47(1):43－53

Xiao Bo.2011.High oxidized magmatic-hydrothermal ore-forming process at the Qulong giant porphyry Cu－Mo deposit，Gangdese，Tibet［D］.Institute of Geology and Geophysics，Chinese Academy of Sciences.Ph.D.thesis:1－275

Xiao Bo，Li Guang-ming，Qin Ke-zhang，Li Jin-xiang，Zhao Jun-xing，Liu Xiao-bin，Xia Dai-xiang，Wu Xing-bing，Peng Zhen-wei.2008.Magmatic intrusion center and mineralization center of Qulong porphyry Cu－Mo deposit in Tibet:evidence from fissure veinlets and mineralization intensity［J］.Mineral Deposits，27:200－208(in Chinese with English abstract)

Yang Z M，Hou Z Q，White N C，Chang Z S，Li Z Q，Song Y..2009.Geology of the postcollisional porphyry copper-molybdenum deposit at Qulong，Tibet［J］.Ore Geology Reviews，36:133－159

Yang Zhi-ming 2008.The Qulong giant porphyry copper deposit in Tibet:magmatism and mineralization［D］.Beijing:Unpublished Ph.D.Thesis，Chinese Academy of Geological Sciences:1－145(in Chinese with English abstract)

Ye Hui-shou.2006.The Mesozoic Tectonic Evolution and Pb－Zn－Ag Metallogeny in the south margin of North China Craton［D］.Beijing:Unpublished Ph.D.Thesis，Chinese Academy of Geological Sciences:1－217(in Chinese with English abstract)

Zhao Yi-ming，Lin Wen-wei，Bi Cheng-si，Li Da-xin，Jiang Chong-jun.1990.Skarn deposits of China［M］.Beijing:Geological Publishing House:1－354(in Chinese)

Zhou Tao-fa，F(xiàn)an Yu，Yuan Feng.2008.Advances on petrogenesis and metallogeny study of the mineralization belt of the Middle and Lower Reachese of the Yangtze River area［J］.Acta Petrologica Sinica，24(8):1665－1678(in Chinese with English abstract)

［附中文參考文獻(xiàn)］

陳 雷，秦克章，李光明，李金祥，肖 波，江化寨，趙俊興，范 新，江善元.2012.西藏岡底斯南緣努日銅鎢鉬礦床地質(zhì)特征與矽卡巖礦物學(xué)研究［J］.礦床地質(zhì)，31(3):417～437

紀(jì)偉強.2010.藏南岡底斯巖基東段花崗巖時代與成因［D］.中國科學(xué)院地質(zhì)與地球物理研究所，博士學(xué)位論文:1－141

李光明，芮宗瑤，王高明，林方成，劉 波，佘宏全，豐成友，屈文俊.2005.西藏岡底斯成礦帶甲瑪和知不拉銅多金屬礦床的Re－Os同位素年齡及其意義［J］.礦床地質(zhì)，24(5):481－489

李靖輝，2008.大別山(北麓)斑巖型鉬礦床成礦系列及成礦規(guī)律［J］.東華理工大學(xué)學(xué)報(自然科學(xué)版)，31(1):25－30

孟祥金，侯增謙，李振清.2006.西藏驅(qū)龍斑巖銅礦S、Pb同位素組成:對含礦斑巖與成礦物質(zhì)來源的指示［J］.地質(zhì)學(xué)報，80(4):554－560

孟祥金，候增謙，高永豐，黃 衛(wèi)，曲曉明，屈文俊.2003.西藏岡底斯成礦帶驅(qū)龍銅礦Re－Os年齡及成礦學(xué)意義［J］.地質(zhì)評論，49(6):660－666

秦克章，汪東波，王之田，孫 樞.1999.中國東部銅礦床類型、成礦環(huán)境、成礦集中區(qū)與成礦系統(tǒng)［J］.礦床地質(zhì)，18(4):359－371

秦克章，夏代祥，李光明，肖 波，多 吉，蔣光武，趙俊興.2014.西藏驅(qū)龍斑巖－矽卡巖銅鉬礦床［M］.科學(xué)出版社:1－316

芮宗瑤，侯增謙，曲曉明，張立生，王龍生，劉玉琳.2003.岡底斯斑巖銅礦成礦時代及青藏高原隆升［J］.礦床地質(zhì)，22(3):217－225

佘宏全，豐成友，張德全，潘桂棠，李光明.2005.西藏岡底斯中東段矽卡巖銅－鉛－鋅多金屬礦床特征及成礦遠(yuǎn)景分析［J］.礦床地質(zhì)，24(5):508－520

唐菊興，鄧世林，鄭文寶，應(yīng)麗娟，汪雄武，鐘康惠，秦世鵬，丁 楓，黎楓佶，唐曉倩，鐘裕峰，彭慧娟.2011.西藏墨竹工卡縣甲瑪銅多金屬礦床勘察模型［J］.礦床地質(zhì)，32(2):179－196

唐仁鯉，羅懷松.1995.西藏玉龍斑巖銅(鉬)礦帶地質(zhì)［M］.北京:地質(zhì)出版社:1－320

西藏巨龍銅業(yè)有限公司.2008.西藏自治區(qū)墨竹工卡縣驅(qū)龍礦區(qū)銅多金屬礦勘探報告［R］.

西藏巨龍銅業(yè)有限公司.2012.西藏自治區(qū)墨竹工卡縣知不拉礦區(qū)銅多金屬礦資源儲量核實報告［R］

肖 波，李光明，秦克章，李金祥，趙俊興，劉小兵，夏代祥，吳興炳，彭震威.2008.岡底斯驅(qū)龍斑巖銅鉬礦床的巖漿侵位中心和礦化中心:破裂裂隙和礦化強度證據(jù)［J］.礦床地質(zhì)，27(2):200－208

肖 波，秦克章，李光明，李金祥，陳 雷，趙俊興，范 新.2011.岡底斯驅(qū)龍斑巖銅－鉬礦區(qū)外圍矽卡巖型銅礦的分布、特征及深部找礦意義［J］.地質(zhì)與勘探，47(1):43－53

肖 波.2011.岡底斯驅(qū)龍巨型斑巖銅－鉬礦床高氧化巖漿－熱液成礦過程［D］.中國科學(xué)院地質(zhì)與地球物理研究所，博士論文:1－275

楊志明.2008.西藏驅(qū)龍超大型銅礦床－巖漿作用與礦床成因［D］.中國地質(zhì)科學(xué)院，博士論文:1－145

葉會壽.2006.華北陸塊南緣中生代構(gòu)造演化與鉛鋅銀礦成礦作用［D］.中國地質(zhì)科學(xué)院博士研究生論文:1－270

趙一鳴，林文蔚，畢承思，李大新，蔣崇俊.1990.中國矽卡巖礦床［M］.北京:地質(zhì)出版社:1－354

周濤發(fā)，范 裕，袁 峰.2008.長江中下游成礦帶成巖成礦研究進展［J］.巖石學(xué)報，24(8):1665－1678