西藏多不雜銅礦床鍶、釹、鉿同位素特征及其地質意義

何陽陽,溫春齊,劉顯凡

(1.內江師范學院 地理與資源科學學院， 四川 內江 641112；2.成都理工大學 地球科學學院， 成都 610059)

多不雜銅礦床所處位置隸屬于西藏改則縣，據西藏地調院報告，該礦床目前提交的銅、金資源儲量均已達到超大型規模[1]。巖石或礦物的Rb-Sr同位素組成常用于判別樣品的源區類型，由于自然界存在87Rb衰變形成87Sr，因此，含銣巖石或礦物中精確的鍶同位素組成取決于那種巖石或礦物的年齡及wRb/wSr比值[2](w表示質量分數)。釤和釹都是稀土元素，地殼中含量很少，在礦床研究中Sm-Nd同位素體系常用來分析物質來源、確定基底巖石年齡和成礦時代[2]；鋯石中的wLu/wHf值一般很小，w(176Hf)/w(177Hf)比值常用來分析其Hf同位素組成，進而反演原巖的成因性質[3]。本文通過對多不雜銅礦床相關樣品中的銣鍶、釤釹、镥鉿同位素進行研究，以此來示蹤其成礦物質來源。

1　礦床地質概況

班公湖-怒江縫合帶是拉薩地塊和羌塘地塊之間的分界線，總長度約為2 000 km，呈近東西方向展布，主要出露SSZ型蛇綠混雜巖[4-6]。伴隨著一系列構造運動及巖漿活動，大量中酸性侵入體在班怒西段羌塘地塊的南緣增生形成構造-巖漿弧[7-8]，本文研究的多不雜銅礦床位于其中(圖1-A)。

礦區依次出露下侏羅統曲色組第二巖段(J1q2)、下白堊統美日切組(K1m)、新近系康托組(N1k)、第四系等地層。F2和F10為多不雜礦區內的主要斷裂構造，周邊次級裂隙較為發育。礦區出露巖漿巖主要有花崗閃長斑巖、石英閃長巖、輝綠巖等，呈巖脈、巖株產出。

礦體走向為東-西方向，傾向南，長度約1.6 km，厚度約0.5 km，延深約0.6 km，主要產于花崗閃長斑巖體及與其接觸的變質砂巖中(圖1-B)，產出部位主要為鉀化帶、石英絹云巖化帶及青磐巖化帶。主要礦石構造包括：脈狀及細脈浸染狀構造，次為塊狀及團塊狀構造等；主要礦石結構包括：交代結構和結晶結構。黃銅礦和黃鐵礦為主要的礦物成分，其次為輝鉬礦、斑銅礦、黝銅礦、閃鋅礦、磁鐵礦等[9]。

多不雜銅礦床圍巖蝕變十分強烈，主要包括鉀硅化帶、絹英巖化帶和青磐巖化帶3種類型。從平面上看，總體表現出近環帶狀；從剖面上看，總體表現為“中心式”的面型蝕變。

2　樣品特征

2.1　花崗閃長斑巖

花崗閃長斑巖(γδπ53)是多不雜銅礦區內的主要成礦斑巖體(圖2)，呈東西向展布，長約2 km，寬100～320 m，面積約0.5 km2。斑巖體總體上為向南傾斜的“板狀體”，北側界線處產狀較緩，約40°；南側界線處相對較陡。巖體與曲色組第二巖段的變質砂巖呈侵入接觸關系，接觸帶部位蝕變強烈[10]。

塊狀構造和斑狀結構為花崗閃長斑巖的主要礦石組構，顯微鏡下觀測見不等粒結構。斑晶和基質是花崗閃長斑巖的2個重要組成部分，其中，斑晶所占體積比例為35%～40%，基質占整個巖石體積的比例為60%～65%。斑晶的構成主要為斜長石(體積分數約60%)、石英(約30%)、鉀長石(約5%)、角閃石(約3%)、黑云母(約2%)?；|主要為花崗結構、粒晶結構等。

2.2　黃鐵礦

黃鐵礦為礦床中分布最廣、數量最多的金屬礦物之一(圖3)?？傮w上看，從花崗閃長斑巖到接觸帶再到變質砂巖，從鉀化帶依次到黃鐵絹英巖化帶、硅化泥化帶、青磐巖化帶，黃鐵礦表現為礦化程度逐漸增強。其中在黃鐵絹英巖化帶的外圈礦化強度最高，可見大量呈團塊狀、脈狀、浸染狀的黃鐵礦，其次是青磐巖化帶的內圈。隨著遠離礦體到圍巖，黃鐵礦的礦化強度逐漸減弱。黃鐵礦化分布區的中心地帶總體上與銅礦化帶相對應[10]。

圖1　多不雜銅礦床大地構造位置及礦體平面圖Fig.1　Tectonic location and plane graph of ore body in the Duobuza copper depositQ4.第四系； N1k.康托組； K1m1.美日切組第一段； J1q2.曲色組第二段； γδπ53.花崗閃長斑巖； ν.灰綠色輝長巖； αβ53.玄武質安山巖； βμ53.墨綠色輝綠巖。1.斷層； 2.礦體范圍； 3.勘探線； 4.鉆孔

圖2　多不雜銅礦床花崗閃長斑巖勘探剖面示意圖Fig.2　Sketch showing exploration sections for the granodiorite porphyry in the Duobuza copper deposit

在原生礦體中，黃鐵礦的粒徑多數為0.01～0.8 mm，最大粒徑1 mm左右，常呈脈狀、浸染狀分布。脈石礦物及相關硫化物常對呈他形粒狀的浸染狀黃鐵礦進行交代，可見黃銅礦等礦物分布其中。在表生帶中，黃鐵礦常被褐鐵礦交代。局部可見硬石膏與黃鐵礦伴生。

2.3　鋯石

河北省區調所利用重選和磁選的方法分選出相關鋯石，陰極發光圖像在西北大學完成。由圖4可見，鋯石呈長柱狀，長度在100～350 μm，長軸與短軸的比例為(2∶1)～(3.5∶1)，有震蕩環帶現象。

圖3　多不雜銅礦床黃鐵礦宏觀及微觀照片Fig.3　Macroscopic and microscopic photos showing pyrite from the Duobuza copper deposit(A)標本照片;(B)反射單偏光鏡下照片

圖4　DB048樣品代表性鋯石陰極發光圖像及測試位置Fig.4　DB048 CL images showing representative zircons and tested positions

3　鍶、釹、鉿同位素特征

3.1　鍶、釹同位素

核工業北京地質研究院參考相關規范和方法對本文相關樣品進行了鍶、釹同位素測試，結果如表1所示，其中銣、鍶的誤差以2σ計。

由表1可見，5件黃鐵礦樣品(包括1件圍巖)Rb的質量分數(wRb)為(0.355～9.19)×10-6，平均為2.279×10-6；Sr的質量分數(wSr)為(4.3～74.7)×10-6，平均為26.92×10-6；87Rb/86Sr為0.081 8～0.355 9，平均為0.176 2；87Sr/86Sr為0.708 429～0.709 167，平均為0.708 909。其中圍巖黃鐵礦樣品中的Rb、Sr含量與其他4個樣品相比明顯偏高。本文測試的全巖樣品中Rb的質量分數為(107～142)×10-6，平均為122.6×10-6；Sr的質量分數為(203～393)×10-6，平均為302.4×10-6；87Rb/86Sr比值為0.791 8～1.938 8，平均為1.252 1；87Sr/86Sr比值為0.707 419～0.709 968，平均為0.708 344。

周雄等[12]認為鍶初始值在花崗質巖石起源的認識上具有重要指示意義。多不雜銅礦床相關樣品的Sr初始值位于全球殼源硅鋁質巖和幔源鎂鐵質巖Sr同位素的平均值之間(其值分別為0.720±0.005和0.703 50)，且相對靠近幔源巖石[13]，說明樣品具有殼?；烊咎卣?。

在地質作用過程中，鍶同位素與成礦元素的遷移方式較為相似[14]。圍巖樣品鍶初始值較斑巖樣品鍶初始值高，表明成礦流體的遷移方式是由斑巖體→圍巖，在這一遷移過程中與圍巖物質發生充填交代作用從而導致Sr含量增高。

143Nd/144Nd比值在一定程度上可以示蹤地質體殼幔物質之間的相互作用?；◢忛W長斑巖中(87Sr/86Sr)i和143Nd/144Nd原始地?，F代值分別為0.704 5和0.512 638。將表1數據與其對比，顯示礦區樣品具有高(87Sr/86Sr)i、低143Nd/144Nd的特征。在圖5中，樣品投影點集中在地幔演化線周圍。富集地幔單元EMⅠ和EMⅡ通過俯沖帶混染了大陸殼物質，是形成上述特征的主要原因。這樣可以認為多不雜銅礦區樣品的特征是交代富集地幔的反映。

3.2　镥鉿同位素

委托西北大學大陸動力學國家重點實驗室進行鋯石原位Lu-Hf同位素微區測定，使用的方法和儀器設備參數等詳見Yuan等[16]，測試結果見表2。

將表2數據在εHf(t)-年齡圖解中投點(圖6)，可以看到，樣品投影點均處于下地殼和虧損地幔之間，靠近球粒隕石線的上面，進一步表明樣品具備虧損地幔的特征，同時有地殼物質加入混染。另外，從表2以及圖7中可以看到，鋯石εHf(t)值全部為正，分布范圍較廣，集中于4～5之間。正的εHf(t)值指示有幔源組分參與成巖成礦，且幔源成分占據主動[3]。該礦床的鋯石樣品具有上述特征，指示伴有多種組分參與寄主巖石的形成，進一步暗示地幔流體作用下有殼源流體參與混染。

*數據來自文獻[11];εNd(t)=[(143Nd/144Nd)t/(143Nd/144Nd)CHUR-1]×104,t=120 Ma

圖5　(87Sr/86Sr)i -εNd(t)圖解Fig.5　(87Sr/86Sr)i -εNd(t) diagram(作圖方法據文獻[15])BSE．原始未分異地幔單元；HIMU．高U/Pb地幔單元；EMⅠ．Ⅰ型富集地幔；EMⅡ．Ⅱ型富集地幔； PREMA．流行地幔單元； DM．虧損地幔單元

圖6　花崗閃長斑巖鋯石年齡與εHf(t)關系圖解Fig.6　Relationship diagram showing zircon age and εHf(t) of granodiorite porphyry(作圖方法據文獻[17])

4　討 論

前人在該區進行勘查和研究時發現，礦體同時賦存于花崗閃長斑巖和變質砂巖中[18-22]，且對該礦床進行過成因方面的研究，在成礦物質來源方面有3種觀點，分別認為其來源于斑巖巖漿、圍巖以及二者混合，但大多數數據都顯示幔源組分參與了成礦[23-26]。

表2　多不雜銅礦床花崗閃長斑巖的鋯石Hf同位素組成Table 2　Zircon Hf isotopic compositions of the granitic porphyry in the Duobuza copper deposit

續表2

樣號t/Ma176Yb/177Hf176Lu/177Hf176Hf/177HfεHf(0)εHf(t)tDM1/MatDM2/Maf Lu/HfDW2-814114.00.03440.00130.2828482.75.1577845-0.96DW2-815116.00.02070.00110.2828954.36.8508739-0.97DW2-816130.00.01590.00080.2828021.13.8635938-0.98DW2-817123.00.01410.00070.282742-1.11.67181077-0.98DW2-818131.00.00920.00050.2828402.45.2577850-0.99DW2-819123.00.01700.00090.2828272.04.6601886-0.97DW2-820233.00.01820.00090.282633-4.90.18741257-0.97Dbz-γδπ-01121.00.00980.00050.2828151.54.1612914-0.99Dbz-γδπ-02121.00.03350.00160.2828623.25.7562813-0.95Dbz-γδπ-03124.00.01790.00090.2828251.94.5604890-0.97Dbz-γδπ-04128.00.01260.00060.2828312.14.8592874-0.98Dbz-γδπ-05128.00.01290.00060.2828091.34.1622922-0.98Dbz-γδπ-06124.00.02210.00100.2828201.74.3614903-0.97Dbz-γδπ-07122.00.01710.00080.2828141.54.1618916-0.98Dbz-γδπ-08119.00.02620.00120.2828322.14.6599879-0.96Dbz-γδπ-09125.00.01790.00080.2828181.64.3613905-0.97Dbz-γδπ-10122.00.03260.00150.2828302.04.6607884-0.96Dbz-γδπ-11120.00.01680.00080.2828312.14.7594879-0.98Dbz-γδπ-12118.00.02250.00110.2828613.15.6556814-0.97Dbz-γδπ-13127.00.02120.00100.2828302.04.7599878-0.97Dbz-γδπ-14122.00.02640.00130.2828382.34.9592864-0.96Dbz-γδπ-15120.00.00970.00050.2828623.25.8546808-0.99Dbz-γδπ-16125.00.01720.00080.2828141.54.2618913-0.98Dbz-γδπ-17111.00.02470.00110.2828733.65.9538789-0.97Dbz-γδπ-18118.00.02100.00100.2828753.76.2535781-0.97Dbz-γδπ-19120.00.03280.00150.2828984.57.0509731-0.95Dbz-γδπ-20120.00.02020.00090.2828292.04.6599884-0.97

DB048系列18個樣品為本文數據，其余數據引自文獻[11]；t為鋯石的結晶年齡

本文鍶、釹、鉿同位素綜合研究亦顯示出花崗閃長斑巖具有殼?；烊咎卣?，幔源組分參與了成礦。有學者認為在花崗質巖石形成及成礦時，幔源組分不僅提供熱能，而且提供礦質[27]。多不雜礦區發現有較強的磁鐵礦化現象，指示巖漿-熱液體系有較強烈的氧化性，地幔流體加入地表流體循環，進而進行能量和礦質的轉換，這一地質事件記錄了地幔流體作用導致該礦區大規模成礦的事實。在這一過程中，不僅可以通過降低原流體中還原性礦物的溶解度而促使其沉淀結晶，而且可導致相對氧化的礦物和典型還原礦物的相互疊加共存，如黃鐵礦與硬石膏的伴生。

多不雜銅礦床符合地幔流體成礦作用的條件[28]，具有地幔流體成礦作用的特征[29]，其成礦動力學模型如圖8所示。羌塘地體和拉薩地體發生碰撞閉合以及俯沖，導致該區發育深大斷裂，誘發深部地質作用。在其作用下，原始地幔一部分轉變為交代富集地幔，另一部分發生熔融，形成含有礦質的地幔流體；同時，板塊的俯沖形成了一定的富銅鈣堿性熔漿。它們相互混染，誘發地幔底辟作用，導致深部地殼部分熔融，三者混染在一起沿斷裂通道上升結晶，形成含銅的花崗閃長斑巖。另有部分混合流體溢出巖漿及巖體，滲入到圍巖的裂隙之中，相互之間會發生能量和礦質的轉換，進而形成多不雜變質砂巖型礦體，即接觸帶成礦體系[30-31]。

5　結 論

a.多不雜銅礦床樣品中的Sr初始值接近幔源鎂鐵質巖石Sr同位素平均值，暗示礦區含礦斑巖的原始巖漿可能來源于地幔，在上升侵位時與硅鋁質地殼物質之間發生充填交代，使其具有殼?；烊咎卣?。

圖7　多不雜銅礦床花崗閃長斑巖的鋯石Hf同位素組成直方圖Fig.7　Zircon Hf isotopic compositions of the granitic porphyry in the Duobuza copper deposit

圖8　多不雜銅礦床成礦動力學模型Fig.8　Metallogenic dynamic model of the Duobuza copper deposit1.地幔流體； 2.板塊俯沖流體； 3.富硅成礦流體； 4.富堿巖漿； 5.銅金礦； 6.矽卡巖銅金礦； 7.矽卡巖鐵礦； 8.矽卡巖銅礦； 9.斑巖銅礦； 10.矽卡巖銅多金屬礦； 11.金礦； 12.城鎮。Ⅲ-1.南巴顏額拉-雅魯藏布江； Ⅲ-3.西金烏蘭-金沙江； Ⅲ-4.喬爾天山-甜水海； Ⅲ-5.北羌塘-昌都； Ⅲ-6.南羌塘-昌都； Ⅲ-7.班公湖-怒江；Ⅲ-8.北岡底斯；Ⅲ-9.中岡底斯； Ⅲ-10.南岡底斯-念青唐古拉； Ⅲ-11.雅魯藏布江； Ⅲ-12.喜馬拉雅

b.花崗閃長斑巖具有相對高(87Sr/86Sr)i以及低143Nd/144Nd的特征，鍶同位素的遷移行為與成礦元素相似，含礦流體由斑巖體向圍巖進行了遷移，并與圍巖之間發生充填交代，導致變質砂巖含礦。

c.鋯石εHf(t)值為正值，且變化范圍較大，是有較多幔源組分參與成巖的標志，指示幔源巖漿作用過程中伴有地殼流體的混染。綜合分析表明，在地幔流體作用下該礦床符合透巖漿流體成礦理論中的接觸帶成礦體系。