• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    西藏玉龍成礦帶斑巖Cu-Mo礦床地質(zhì)地球化學(xué)特征及成礦機制探討
    ——玉龍和多霞松多對比研究

    2013-08-11 14:00:03吳偉中張玉泉董冰華夏中曦
    大地構(gòu)造與成礦學(xué) 2013年3期
    關(guān)鍵詞:玉泉含礦玉龍

    吳偉中, 夏 斌, 張玉泉 董冰華, 夏中曦

    (1.中國科學(xué)院 廣州地球化學(xué)研究所, 廣東 廣州 510640; 2.中國科學(xué)院大學(xué), 北京 100049; 3.中山大學(xué)海洋學(xué)院, 廣東 廣州 510275)

    西藏玉龍成礦帶斑巖Cu-Mo礦床地質(zhì)地球化學(xué)特征及成礦機制探討
    ——玉龍和多霞松多對比研究

    吳偉中1,2, 夏 斌1,3, 張玉泉1, 董冰華1,2, 夏中曦3

    (1.中國科學(xué)院 廣州地球化學(xué)研究所, 廣東 廣州 510640; 2.中國科學(xué)院大學(xué), 北京 100049; 3.中山大學(xué)海洋學(xué)院, 廣東 廣州 510275)

    玉龍斑巖銅礦帶位于特提斯-喜馬拉雅構(gòu)造域東緣三江構(gòu)造帶的中部, 是研究碰撞后喜山期富堿斑巖成巖成礦作用響應(yīng)的最佳地區(qū)。本文選擇玉龍和多霞松多兩個礦床, 分析了兩地含礦斑巖的成巖成礦年齡及巖石地球化學(xué)特征,并結(jié)合Sr-Nd-Pb、S-H-O同位素進一步探討了與岡底斯含礦斑巖巖漿源區(qū)組成、構(gòu)造控礦模式及深部動力學(xué)機制上的差異。賦礦巖體均具有富堿、高K、高Al等鉀質(zhì)堿性巖系列特征, 玉龍早期為石英二長斑巖, 晚期為花崗閃長斑巖和石英二長斑巖, 多霞松多早期為二長花崗斑巖和花崗閃長斑巖, 晚期為堿長花崗斑巖, 兩者明顯富集輕稀土和大離子親石元素Rb、Th、U、Pb, 相對虧損高場強元素Nb、Ta、Yb, 無明顯Eu、Ce負異常。略高于原始地幔現(xiàn)代值的87Sr/86Sr值、較低的εNd(t)、相對均一的206Pb/204Pb值, 指示斑巖巖漿可能源自接近于EMⅡ型富集地幔交代成因的石榴石角閃巖和/或角閃榴輝巖的部分熔融, 并發(fā)生低程度結(jié)晶分異, 在上升侵位途中受到地殼物質(zhì)的混染, 其形成于古近紀中始新世喜馬拉雅期的37~41 Ma左右, 可能受控于印度板塊與歐亞板塊碰撞誘發(fā)的紅河-哀牢山走滑斷裂系統(tǒng),尤其是妥壩-芒康斷裂(TBF)左行走滑產(chǎn)生的局部俯沖作用。

    富堿斑巖;地球化學(xué)特征;成礦機制;玉龍銅礦帶;藏東

    西藏東部玉龍斑巖銅礦帶是在陸-陸碰撞造山環(huán)境下形成的世界級大型銅礦帶, 帶內(nèi)現(xiàn)已發(fā)現(xiàn)1個超大型斑巖銅礦床(玉龍)、2個大型斑巖銅礦床(馬拉松多、多霞松多)和2個中小型斑巖銅礦床及20多個礦化巖體(芮宗瑤等, 1984; 馬鴻文, 1989, 1990;唐仁鯉等, 1995), 是哀牢山-金沙江鉀質(zhì)堿性巖帶的重要組成部分。自20世紀80 年代初以來, 眾多學(xué)者對玉龍銅礦帶的礦床地質(zhì)特征(芮宗瑤等, 1984;唐仁鯉等, 1995)、花崗巖類與成礦的關(guān)系(馬鴻文, 1990; 張玉泉等, 1998)、蝕變與成礦的關(guān)系(周宜吉, 1980, 1985)、流體包裹體特征(李蔭清, 1984, 1985)、區(qū)域構(gòu)造與成礦(唐仁鯉等, 1995; Hou et al., 2003)、斑巖型銅礦的成礦構(gòu)造背景(芮宗瑤等, 1984; 姜耀輝等, 2006; 陳建平等, 2009)、成巖成礦年齡(Liang et al., 2005; 郭利果等, 2006; 曾普勝等, 2006; 梁華英等, 2008; 伍靜等, 2011)等開展了大量的研究工作。目前大部分的研究工作主要集中在玉龍斑巖銅礦床上, 而對多霞松多斑巖銅礦床則研究相對較少, 對含礦斑巖產(chǎn)出的構(gòu)造背景、巖體屬性和成巖成礦時代等問題還存在爭議, 尤其是對玉龍富堿斑巖帶的巖漿起源演化與成因機制尚未取得一致的認識(唐仁鯉等, 1995; 張玉泉等, 1998; Hou et al., 2003; Liang et al., 2006a, b)。本文在以往的工作基礎(chǔ)上, 重點選擇了己知礦床(點)分布相對集中的玉龍和多霞松多兩子區(qū)作為研究區(qū), 通過研究兩地含礦斑巖體的巖石學(xué)、年代學(xué)及其地球化學(xué)特征及與岡底斯含礦斑巖的對比分析, 并進一步討論這三者物質(zhì)源區(qū)組成、構(gòu)造控礦模式與深部過程以及動力學(xué)機制上的區(qū)別和差異,從而為更好的認識青藏高原東部玉龍成礦帶含礦斑巖的起源及其構(gòu)造演化提供新的地質(zhì)依據(jù)。

    1 礦床地質(zhì)特征和形成時代

    圖1 玉龍斑巖銅礦帶區(qū)域構(gòu)造位置圖(據(jù)Wang et al., 2001; Hou et al., 2003; 唐菊興等, 2006修改)Fig. 1 Tectonic location of the Yulong porphyry copper belt (after Wang et al., 2001; Hou et al., 2003; Tang et al., 2006)

    表1 玉龍和多霞松多斑巖型礦床基本地質(zhì)特征對比Table 1 Comparison of the characteristics of the porphyry ore deposits in the Yulong and Duoxiasongduo areas

    玉龍斑巖銅礦床是藏東喜山期玉龍成礦帶中最大的斑巖銅礦床, 位于金沙江-紅河斷裂帶北段西側(cè), 恒星措-甘龍拉短軸背斜的軸部向南傾斜部位(圖1), 主要出露由黑云母二長花崗斑巖、石英二長斑巖、花崗閃長斑巖組成的平面上似梨形巖株產(chǎn)出的玉龍復(fù)式侵入巖體, 侵入時間介于41 Ma和38 Ma之間(唐仁鯉等, 1995; Hou et al., 2003; 梁華英等, 2008), 面積為0.64 km2。金屬銅儲量在628萬噸, 伴生金約100噸, 銅品位0.99%, 金品位 0.35 g/t, 具有世界級規(guī)模(芮宗瑤等, 1984)。地層主要為中-上三疊統(tǒng)波里拉組淺海相灰?guī)r夾砂巖、甲丕拉組陸相紅碎屑巖和阿堵拉組海相砂頁巖。多霞松多斑巖銅礦床處在玉龍銅礦帶偏南部及莽總背斜中段近軸部的南西翼(圖1), 含礦斑巖體主要侵入在上三疊統(tǒng)甲丕拉組下段砂泥巖層中, 是由角閃石、黑云母二長花崗斑巖和堿長花崗斑巖組成的復(fù)式斑巖體(張玉泉和鐘孫霖, 1997)。區(qū)內(nèi)出露的地層主要為上奧陶統(tǒng)(O1)板巖夾變質(zhì)灰?guī)r; 石炭系-二疊系(C-P)灰白色大理巖夾板巖、砂巖及基性火山巖; 上三疊統(tǒng)(T3)甲丕拉組: 中段泥質(zhì)粉砂巖、泥灰?guī)r、灰?guī)r和火山巖, 下段紫紅色砂礫巖、石英砂巖及頁巖夾凝灰?guī)r。礦床平均銅含量在0.378%左右, 含金 0.054 μg/g,含銀 3.0 μg/g, 含鉬 0.021%(唐仁鯉等, 1995), 斑巖成礦年齡介于37~38 Ma(梁華英, 2002; Liang et al., 2006b)。據(jù)眾多學(xué)者研究(芮宗瑤等, 1984; 唐仁鯉等, 1995; 張玉泉等, 1998; Hou et al., 2003), 產(chǎn)出于同一構(gòu)造背景和成礦環(huán)境上的玉龍和多霞松多兩礦床,其基本地質(zhì)特征具有相似性(表1): (1)礦區(qū)侵入巖體積小(直徑<2 km)、侵位淺(0.5~3 km)且多期次侵位、多巖相分帶; (2)含礦斑巖為富堿鉀玄巖系列, 以二長花崗斑巖、堿長花崗斑巖、石英二長斑巖和正長花崗斑巖等為主要礦物, 具斑狀結(jié)構(gòu); (3)斑巖侵入體及圍巖廣泛發(fā)育裂隙構(gòu)造控制的蝕變和礦化。蝕變具有明顯的典型中心式對稱型蝕變特征, 即以巖體為中心, 由內(nèi)向外分帶: 鉀硅酸鹽化→石英-絹云母化、黏土化→青磐巖化; (4)礦體主要產(chǎn)出于斑巖與圍巖接觸帶, 礦體形態(tài)總體呈筒狀、空心筒狀, 礦石以細脈浸染狀為主, 網(wǎng)脈狀次之; (5)金屬硫化物以黃鐵礦、黃銅礦、輝鉬礦、斑銅礦為主, 氧化物以赤鐵礦和磁鐵礦為主, 早期為黃鐵礦+黃銅礦+輝鉬礦+磁鐵礦組合, 晚期為黃鐵礦+斑銅礦+赤鐵礦組合; (6)成礦流體早期以高溫(>400 ), ℃高鹽度(含 34%~53% NaCl)巖漿水為主, 晚期混入大氣水, 并與地下水發(fā)生垂向的氧化淋濾作用, 形成氧化淋濾帶-硫化物富集帶-氧化物富集帶; (7)成巖、成礦年齡主要集中在37~40 Ma, 介于早喜馬拉雅運動高峰期(43.04 Ma)與中喜馬拉雅階段高峰期(28.44 Ma)之間, 成礦作用大約是在1~2 Ma內(nèi)完成的, 與賦礦巖體的年齡基本吻合。

    2 巖石地球化學(xué)特征

    本次獲得的樣品主量元素分析在中國科學(xué)院廣州地球化學(xué)研究所同位素地球化學(xué)國家重點實驗室利用X射線熒光光譜法(XRF)測定, 所有氧化物分析精度優(yōu)于1%; 微量(含稀土)元素采用酸溶法分解樣品, 在中國科學(xué)院廣州地球化學(xué)研究所超凈化室用Perkin-Elmer Sciex Elan 6000型電感耦合等離子體質(zhì)譜儀(ICP-MS)測定完成, 微量元素中含量>10 μg/g的樣品分析精度優(yōu)于5%(2σ), <10 μg/g的分析精度優(yōu)于10%(2σ), 所有稀土元素的分析精度都優(yōu)于5%(2σ), 具體分析流程參閱劉穎等(1996)。

    2.1 主量元素

    為了便于更好地討論, 我們將玉龍和多霞松多含礦斑巖的巖石化學(xué)組成進行對比分析, 從表2可以看出玉龍和多霞松多含礦巖體都具有高Al(Al2O3: 14.18%~18.00%、12.89%~15.24%)和富堿(Na2O+K2O: 7.84%~10.22%、6.81%~8.34%)、高K(K2O/Na2O>1), TiO2含量較低, 多小于0.4%的特征。Mg#值比較高(46.3~87.2、35.7~83.1), 大多數(shù)都大于44。其中玉龍斑巖體的A12O3、CaO、Na2O的含量明顯的高于多霞松多, 尤其是玉龍早期含礦斑巖體與多霞松多晚期相比, 而SiO2的含量略低于多霞松多(圖3)。

    表2 多霞松多含礦斑巖主量元素分析結(jié)果(%)Table 2 Major element contents (%) of the mineralized porphyries in the Duoxiasongduo area

    對玉龍和多霞松多含礦斑巖體類型, 目前有不同的看法, 馬鴻文(1990)和唐仁鯉等(1995)認為玉龍銅礦帶含礦斑巖屬鈣堿性系列; 張玉泉等(1998)和Hou et al.(2003)認為玉龍銅礦帶含礦斑巖均屬于鉀玄巖系列, 玉龍早期為石英二長斑巖, 晚期為正長花崗斑巖, 多霞松多早期為二長花崗斑巖, 晚期為堿長花崗斑巖; 姜耀輝等(2006)認為玉龍含礦斑巖體主要為二長花崗斑巖、堿長花崗斑巖和正長花崗斑巖, 少量的石英二長斑巖, 多霞松多為堿長花崗斑巖和正長花崗斑巖。在SiO2-K2O圖解中(圖2a),大部分巖體樣品均位于鉀玄巖范圍, 只有少量的玉龍晚期樣品落入高鉀鈣堿性巖系列內(nèi), 與鄰區(qū)的莽總(伍靜等, 2011)、馬拉松多(梁華英等, 2009)斑巖體和囊謙盆地火山巖(鄧萬明等, 2001; Sun et al., 2001)及日通-馬牧普富堿侵入巖(張玉泉和鐘孫霖, 1997;張玉泉等, 1998)相類似, 均顯示為鉀玄巖系列, 與岡底斯含礦斑巖分布的高鉀鈣堿性系列(Hou et al., 2004; Qu et al., 2004; Guo et al., 2007)有一定的區(qū)別,可能反映了兩類環(huán)境斑巖巖漿具有不同的起源和演化特性。含礦斑巖又具有埃達克巖的親合性, 屬于拉薩地塊埃達克巖(Middlemost, 1994; 王保弟等, 2010)的范疇, 并且具有很高的K2O含量和較高的K2O/Na2O比值(圖2a)。在侵入巖全堿-SiO2分類圖中(圖2b)(Le Bas et al., 1986), 玉龍早期含礦斑巖體主要落在石英二長巖內(nèi), 晚期主要落在石英二長巖和花崗閃長巖內(nèi), 少量位于花崗巖區(qū), 而多霞松多早期主要落在花崗巖區(qū)和花崗閃長巖區(qū)及兩者的分界線附近, 晚期主要落在花崗巖區(qū)。綜合考慮兩地斑巖體中鉀長石斑晶的含量分布及早期巖體相對貧SiO2、富TiO2、A12O3、MgO、CaO、Na2O, 晚期巖體相對富SiO2、K2O(表2), 可以認為玉龍早期為石英二長斑巖, 晚期主要為花崗閃長斑巖和石英二長斑巖, 多霞松多早期主要為二長花崗斑巖和花崗閃長斑巖, 晚期為堿長花崗斑巖。巖體的SiO2與Al2O3、MgO、CaO、Na2O、TiO2、P2O5均呈良好的負相關(guān)關(guān)系(圖3), 和K2O有較弱的正相關(guān)關(guān)系(圖2a), 而Fe2O3的變化趨勢不明顯(圖3), 表明在巖漿的演化過程中, 可能有角閃石、單斜輝石和鈦鐵礦等相對富鐵鎂礦物的結(jié)晶分異。

    圖2 玉龍和多霞松多含礦斑巖SiO2-K2O(a)和(Na2O+ K2O)-SiO2(b) 圖(據(jù)Le Bas et al., 1986)Fig.2 The diagrams of SiO2vs. K2O(a) and SiO2vs. Na2O+ K2O(b) for the mineralized porphyries from the Yulong and Duoxiasongduo areas in the eastern Tibet(after Le Bas et al., 1986)

    2.2 稀土元素和微量元素

    稀土元素的地球化學(xué)行為與巖石成因關(guān)系十分密切, 玉龍和多霞松多含礦斑巖體的稀土總量變化范圍∑REE分別為156.12~413.25 μg/g、160.89~304.89 μg/g, LREE/HREE為15.48~24.26、8.40~22.97,兩值均變化較大, 表明稀土元素具有強烈分異且LREE明顯富集, 尤其是玉龍早期含礦斑巖體。Sm/Nd值為0.14~0.24, Eu/Sm 值為0.12~0.27, 該比值范圍位于地殼和上地幔之間(王治華等, 2010)。δEu值分別為0.72~1.03、0.57~0.95, 具弱負或無明顯Eu負異常(表3), δCe值為0.89~1.07, 基本上無Ce異常(虧損), 在稀土元素配分模式圖(圖4a)中表現(xiàn)為較平緩的右傾型曲線, 與岡底斯帶含礦斑巖弱正或弱負Eu異常變化特征(Hou et al., 2004; Qu et al., 2004; Guo et al., 2007)相似。(La/Yb)N、(Ce/Yb)N的比值均遠遠大于1, 其中玉龍早期含礦斑巖體的該比值最大分別為46.55、32.40, 可能說明巖漿形成演化過程中結(jié)晶分異作用不占主導(dǎo)地位(馬鴻文, 1990)或源區(qū)缺少富含Eu的礦物。

    DWD無線隨鉆測量儀器由地面部分、井下部分及輔助工具、設(shè)備組成,其工作原理是利用儀器內(nèi)的渦輪發(fā)電機將部分鉆井液的能量轉(zhuǎn)化成電能為探管供電,同時鉆井液脈沖發(fā)生器將探管探測到的數(shù)據(jù)通過鉆井液傳遞到地面,地面上采用鉆井液壓力傳感器檢測來自井下儀器的鉆井液脈沖信息,并傳輸?shù)降孛娼獯a箱進行處理,最終井下儀器所測量的井斜角、方位角和工具面數(shù)據(jù)等信息可以顯示在計算機或DDU司鉆閱讀器上。

    在原始地幔標(biāo)準(zhǔn)化蛛網(wǎng)圖(圖4b)中, 玉龍和多霞松多含礦斑巖明顯富集大離子親石元素(LILE)Rb、Th、U、Pb, Hf略微富集, 相對虧損高場強元素(HFSE)Nb、Ta、Yb, Sr含量變化分別為457~995 μg/g、58~776 μg/g,除玉龍晚期樣品顯示為弱正異常外, 其他的具有負異常特征, 多霞松多早期含礦斑巖體則更明顯, 這與岡底斯含礦斑巖弱正Sr異常(Hou et al., 2004; Qu et al., 2004)有一定的區(qū)別, 且具有較明顯的Rb、Th、U、Pb峰與Ba、Nb谷。Y含量變化于9.33~19.69 μg/g、11.11~24.85 μg/g, Sm/Yb比值變化為4.34~7.15、3.45~6.78,均值為5.32, 與Kay and Mpodozis(2001)提出的成礦期斑巖的Sm/Yb比值(5~7)范圍相吻合, Zr/Sm比值變化為15.96~50.02、13.41~49.15, 均值為27.36, 證實玉龍和多霞松多含礦斑巖以較高Y>11 μg/g、中等Sm/Yb和Zr/Sm比值為特征。

    圖3 玉龍和多霞松多含礦斑巖Harker圖(SiO2-MgO, SiO2-Mg#圖據(jù)王強等, 2004, 2008; 侯增謙等, 2007)Fig.3 Harker diagrams for the mineralized porphyries in the Yulong and Duoxiasongduo areas (diagrams of SiO2vs. MgO and SiO2vs. Mg#based on Wang et al., 2004, 2008; Hou et al., 2007)

    圖4 含礦斑巖體稀土元素配分型式圖(a)和微量元素原始地幔蛛網(wǎng)圖(b)(標(biāo)準(zhǔn)化值據(jù)Sun and McDonough, 1989)Fig.4 Chondrite-normalized REE patterns (a) and primitive mantle-normalized trace element spider diagrams (b) of the mineralized porphyries (after Sun and McDonough, 1989)

    3 討 論

    3.1 含礦斑巖的埃達克(質(zhì))巖特征

    埃達克(質(zhì))巖作為一種俯沖洋殼板片MORB在榴輝巖相條件下部分熔融的產(chǎn)物, 具有獨特的地球化學(xué)特征, 如SiO2≥56%、A12O3≥15%、MgO< 3%(很少>6%)、較低的Y含量(≤18 μg/g)和HREE(如Yb≤1.9 μg/g)、相對虧損HFSE(Nb、Ta、Ti、P)以及較高的Sr(>400 μg/g)(Defant and Drummond, 1990)。最近的研究顯示, 埃達克質(zhì)巖與斑巖型Cu-Au-Mo礦床有著一定的關(guān)系(Defant and Kepezhinskas, 2001; Oyarzun et al., 2001; 張旗等, 2004; 王強等, 2008)。玉龍和多霞松多含礦斑巖較高的SiO2(>63%)和A12O3(>15%, 在SiO2=65%條件下), 較高的Sr/Y(大多數(shù)>35)和La/Yb比值, 虧損HREE(Yb≤1.6 μg/g)和Y(大多數(shù)<18 μg/g)以及LREE與HREE強烈分餾和無明顯Eu負異常諸特征(表4), 顯示其具有埃達克巖巖漿親合性。與典型的埃達克巖相比(Defant and Drummond, 1990), 兩地的含礦斑巖則相對高K(K2O/Na2O>1)、較高Mg(MgO=0.51%~1.49%)和富集LILE(Rb、Th), 但總體處于埃達克巖區(qū)范圍內(nèi), 與岡底斯帶斑巖明顯偏離(曲曉明等, 2001; Hou et al., 2004; 孟祥金, 2004)(圖5a, b), 暗示其巖漿源區(qū)MORB質(zhì)角閃榴輝巖或榴輝巖可能發(fā)生了較大程度的部分熔融(Defant and Drummond, 1990), 在地球化學(xué)上顯示出正常島弧火山巖與埃達克巖的過渡特征。

    表4 玉龍、多霞松多含礦斑巖與典型adakite地球化學(xué)特征對比Table 4 Comparison of geochemical characteristics of the typical adakite and the ore-bearing porphyries in the Yulong and Duoxiasongduo areas

    圖5 含礦斑巖Sr/Y-Y圖((a), 據(jù)Drummond and Defant, 1990; Martin, 1986)和(La/Yb)N-(Yb)N圖((b), 據(jù)Defant and Drummond, 1990; Petford and Atherton, 1996)Fig.5 Sr/Y vs. Y(a) (after Drummond and Defant, 1990; Martin, 1986) and (La/Yb)Nvs. (Yb)N(b) (after Defant and Drummond, 1990; Petford and Atherton, 1996) diagrams for the ore-bearing porphyries in the Yulong and Duoxiasongduo areas

    3.2 巖漿源區(qū)與深部過程

    造山型鉀質(zhì)巖通常認為來自于與俯沖板片流體(熔體)有關(guān)的富集地幔源區(qū)(Foley, 1992; Nelson, 1992), 而最近的一些研究表明青藏高原碰撞型含礦斑巖的成因與俯沖增生而導(dǎo)致的加厚下地殼的部分熔融作用有關(guān)(Hou et al., 2004; Qu et al., 2004; Guo et al., 2007; Gao et al., 2007, 2010)。目前對富堿斑巖的成因已提出多種模式, 如: 交代富集地幔部分熔融模式(張玉泉等, 2000; 姜耀輝等, 2006)、殼/幔過渡帶部分熔融模式(鄧萬明等, 1998; 鐘大賚等, 2000)、大陸板片沿紅河斷裂帶向東俯沖模式(Wang et al., 2001)和大規(guī)模走滑誘發(fā)“島弧型”地幔部分熔融模式(Hou et al., 2003; 王建等, 2003)。含礦斑巖的Al2O3、MgO、CaO、Na2O、TiO2、P2O5與SiO2之間存在明顯的負相關(guān)關(guān)系(圖3), 表明在巖漿演化過程中, 可能有角閃石、單斜輝石和鈦鐵礦等相對富鐵鎂礦物的結(jié)晶分異。但對部分熔融作用和結(jié)晶分離作用形成的巖漿巖的La與La/Sm的相關(guān)關(guān)系研究(圖7), 表明其形成的深度至少在40 km以上(陳建林等, 2011), 其形成機制主要以部分熔融為主, 和岡底斯斑巖帶相類似(Hou et al., 2004; Guo et al., 2007)。同時無明顯Eu負異常特性也說明巖漿形成演化過程中結(jié)晶分異作用不占主導(dǎo)地位(馬鴻文, 1990)。該區(qū)晚期具有較高的K2O含量, 并且隨著SiO2含量的增加而有較明顯的增加趨勢, 表明源區(qū)存在富鉀礦物相如金云母、富鉀角閃石和鉀長石的可能性。巖石中Ti的虧損可能說明源區(qū)有在含水條件下穩(wěn)定的金紅石礦物相的殘留, Nb、Ta的虧損與Hf的富集表明源區(qū)殘留相中可能有角閃石的存在。Yb的虧損則意味著巖漿源區(qū)有石榴石殘留, K、Rb的富集與Ba的虧損表明可能與源區(qū)金云母的熔融有關(guān)。由以上特征推斷, 玉龍和多霞松多斑巖源區(qū)含金云母, 殘留相中富集石榴石和角閃石。HREE極度虧損以及與LREE的強烈分異表明斑巖源區(qū)巖漿殘留相中存在有石榴石或/和普通角閃石(孟祥金, 2004), 這與微量元素所反映出的信息一致。據(jù)Y、Yb、Sm和Zr元素地球化學(xué)特征及與含礦斑巖共生的正長斑巖和二長斑巖中廣泛發(fā)現(xiàn)的角閃巖、石榴石角閃巖和榴輝巖等下地殼包體(蔡新平, 1992; 趙欣等, 2004), 表明具有較高Y(>11 μg/g)含量的玉龍和多霞松多含礦斑巖可能起源于含石榴石較少的石榴石角閃巖和角閃榴輝巖(圖5)。

    含礦斑巖具有略高于原始地?,F(xiàn)代值(0.7045)的87Sr/86Sr比值(0.7055~0.7075)、較低的εNd(t)(-2.0~-4.2)(表5), 不同于典型的大洋板片熔融(Kay, 1978; Kay et al., 1993; Stern and Kilian, 1996)和中國東部地殼熔融(朱訓(xùn)等, 1983; 常印佛等, 1991)形成的埃達克巖, 介于虧損地幔(MORB)、古老下地殼與富集地幔(EMⅡ)形成的三角形區(qū)域內(nèi), 與岡底斯含礦斑巖相比(Hou et al., 2004; Guo et al., 2007), 玉龍和多霞松多兩地則更接近于富集地幔源區(qū)EMⅡ(圖6a)。較窄變化范圍的206Pb/204Pb比值(18.713~18.890)(表5)和由MORB外圍向EMⅡ內(nèi)直線延伸的207Pb/204Pb比值變化(圖6b), 有著與日通-馬牧普富堿侵入巖(張玉泉和鐘孫霖, 1997; 張玉泉等, 1998)相類似的分布區(qū)域, 同樣暗示埃達克質(zhì)熔體曾與富集地幔物質(zhì)發(fā)生了物質(zhì)交換, 特別是EMⅡ型富集地幔部分熔融產(chǎn)生的鉀質(zhì)-超鉀質(zhì)玄武質(zhì)熔體(Turner et al., 1993), 這種混合作用導(dǎo)致了含礦斑巖的Mg#增高、K2O含量增加及87Sr/86Sr比值增高。斑巖體和圍巖的δ34S值變化范圍為-0.5‰~5.97‰(芮宗瑤等, 1984; 陳建平等, 2009), 與超基性巖、基性巖及石隕石的硫同位素組成相似(從柏林等, 1979), 反映出該區(qū)礦體的硫源主要來自上地?;蛳碌貧? 以地幔物質(zhì)為主。流體包裹體的氫氧同位素的分析測定(δD:-75.1‰~-116.7‰, δ18O: 7.1‰~23.6‰)(陳建平等, 2009), 表明含礦斑巖巖漿在上升侵位途中受到了地殼物質(zhì)的混染作用。

    3.3 成巖動力學(xué)機制

    印度板塊與歐亞板塊在65~41 Ma發(fā)生大規(guī)模碰撞之后, 青藏高原即進入陸內(nèi)碰撞造山階段, 之后相繼經(jīng)歷了晚碰撞(40~26 Ma)和后碰撞(25~0 Ma)(侯增謙等, 2006a)。始于40 Ma的晚碰撞造山作用, 發(fā)生于印度-亞洲大陸的持續(xù)會聚和南北擠壓背景之下, 以大陸內(nèi)部地體(陸塊)間的相對運動,即陸內(nèi)俯沖和逆沖-推覆-走滑活動為特征(侯增謙等, 2006b)。伴隨著第一幕大規(guī)模的應(yīng)力釋放(芮宗瑤等, 2004), 在青藏高原東南緣形成長達數(shù)千公里的金沙江-紅河走滑拉分帶(張玉泉和鐘孫霖, 1997),發(fā)育一系列派生性的含新生代富堿侵入巖和鉀質(zhì)火山巖的拉張盆地(如囊謙、貢覺、劍川)。玉龍斑巖銅礦帶位于紅河-哀牢山巨型走滑斷裂帶北延斷裂系統(tǒng)中, 西側(cè)為妥壩-芒康斷裂, 成礦巖體年齡從北西往南東逐漸降低, 微量元素及Sr-Nd-Pb同位素組成特征(張玉泉等, 2000; Jiang et al., 2006)和沿紅河-哀牢山左行走滑斷裂分布的鉀質(zhì)堿性巖相似(Wang et al., 2001), 都具交代富集地幔特征(張玉泉等, 2000; Hou et al., 2003; 姜耀輝等, 2006)。因此, 可以認為玉龍銅礦帶是由于受印度板塊與歐亞板塊碰撞誘發(fā)的紅河-哀牢山走滑斷裂系統(tǒng)尤其是妥壩-芒康斷裂左行走滑斷裂活動產(chǎn)生的局部俯沖作用, 把下地殼物質(zhì)帶入巖石圈地幔誘發(fā)其部分熔融并經(jīng)歷低程度的分離結(jié)晶作用演化(圖7), 之后大規(guī)模走滑斷裂和走滑拉分盆地(貢覺右行走滑拉分盆地、囊謙左行走滑拉分盆地)在富堿巖漿匯集上升過程中控制了巖漿熔體的淺成侵位和噴發(fā)(Richards, 2003; 梁華英等, 2006)(圖8)。

    表5 含礦斑巖Sr-Nd-Pb同位素組成Table 5 Sr-Nd-Pb isotopic compositions of the ore-bearing porphyries

    圖6 含礦斑巖的87Sr/86Sr-143Nd/144Nd圖(a)和207Pb/204Pb-206Pb/204Pb圖(b)Fig.6 Sr-Nd-Pb isotopic characteristics of the mineralized porphyries in the Yulong and Duoxiasongduo areas

    圖7 含礦斑巖La-La/Sm圖解Fig. 7 La vs. La/Sm diagram for the mineralized porphyries in the Yulong and Duoxiasongduo areas

    4 結(jié) 論

    (1) 玉龍和多霞松多含礦斑巖均屬于鉀玄質(zhì)侵入巖, 具有富堿、高K、高Al等鉀質(zhì)堿性巖系列的成分特征, 但是它們的主量元素含量在早、晚期有著較小的差異: 玉龍早期為石英二長斑巖, 多霞松多主要為二長花崗斑巖和花崗閃長斑巖, 富TiO2、A12O3、MgO、CaO、Na2O; 玉龍晚期主要為花崗閃長斑巖和石英二長斑巖, 多霞松多為堿長花崗斑巖,相對富SiO2、Fe2O3和K2O, 表明在巖漿的演化過程中, 可能有角閃石、單斜輝石和鈦鐵礦等富鐵鎂礦物的結(jié)晶分異。

    圖8 玉龍銅礦帶富堿斑巖巖漿發(fā)育構(gòu)造模式示意圖(據(jù)Richards, 2003; 梁華英等, 2006修改)Fig.8 Schematic cross section of the alkali-rich porphyry magmatic ascent and emplacement in the Yulong porphyry copper ore belt (modified from Richards, 2003; Liang et al., 2006)

    (2) 含礦斑巖稀土元素強烈分異且LREE明顯富集, 弱負或無明顯Eu、Ce負異常; 富集LILE(Rb、 Th、U、Pb、Hf), 相對虧損HFSE(Nb、Ta、Yb), Sr負異常, 具有較明顯的Rb、Th、U、Pb峰與Ba、Nb谷, 具有似埃達克巖地球化學(xué)特征, 略高于原始地?,F(xiàn)代值的87Sr/86Sr比值、較低的εNd(t)、較窄變化范圍的206Pb/204Pb比值及S-H-O同位素組成特征, 指示斑巖巖漿源區(qū)可能是接近于EMⅡ型富集地幔交代成因的石榴石角閃巖和/或角閃榴輝巖的部分熔融, 且在上升侵位途中受到了地殼物質(zhì)的混染作用。

    (3) 玉龍銅礦帶的形成可能與中新世受印度板塊與歐亞板塊碰撞誘發(fā)的紅河-哀牢山走滑斷裂系統(tǒng)尤其是妥壩-芒康斷裂(TBF)左行走滑斷裂活動產(chǎn)生的局部俯沖作用, 把下地殼物質(zhì)帶入巖石圈地幔誘發(fā)其部分熔融并經(jīng)歷低程度的分離結(jié)晶作用演化有關(guān), 在古近紀中始新世喜馬拉雅期的37~41 Ma左右, 巖漿侵入形成大型斑巖銅礦床。

    致謝: 野外和實驗工作中中國科學(xué)院廣州地球化學(xué)研究所和中山大學(xué)海洋學(xué)院課題組的師兄弟們給予大力協(xié)助, 特此致謝!張旗研究員及曲曉明研究員對本文提出了寶貴的修改意見, 特此向他們表示衷心感謝!

    蔡新平. 1992. 揚子地臺西緣新生代富堿斑巖中的深源包體及其意義. 地質(zhì)科學(xué), (2): 183-189.

    常印佛, 劉湘培, 吳言昌. 1991. 長江中下游成礦帶. 北京: 地質(zhì)出版社: 1-379.

    陳建林, 許繼峰, 任江波. 2011. 俯沖型和碰撞型含礦斑巖地球化學(xué)組成的差異. 巖石學(xué)報, 27(9): 2733-2742.

    陳建平, 唐菊興, 叢源, 董慶吉, 郝金華. 2009. 藏東玉龍斑巖銅礦地質(zhì)特征及成礦模型. 地質(zhì)學(xué)報, 83(12): 1887-1900.

    從柏林, 張雯華, 葉大年. 1992. 華北斷塊新生代玄武巖的研究. 地質(zhì)學(xué)報, 66(2): 112-123.

    鄧萬明, 黃萱, 鐘大賚. 1998. 滇西新生代富堿斑巖的巖石特征與成因. 地質(zhì)科學(xué), 33(2): 412-425.

    鄧萬明, 孫宏娟, 張玉泉. 2001. 囊謙盆地新生代鉀質(zhì)火山巖成因巖石學(xué)研究. 地質(zhì)科學(xué), 36(3): 304-318.

    董冰華, 夏斌, 張玉泉. 2013. 藏東玉龍鉀質(zhì)堿性花崗巖地球化學(xué)、黒云母40Ar-39Ar年齡及其地質(zhì)意義. 未刊.

    付建明, 馬昌前, 謝才富, 張業(yè)明, 彭松柏. 2005. 湖南金雞嶺鋁質(zhì)A型花崗巖的厘定及構(gòu)造環(huán)境分析. 地球化學(xué), 34(3): 215-226.

    郭利果, 劉玉平, 徐偉, 張興春, 秦克章, 李鐵勝, 石玉若. 2006. SHRIMP鋯石年代學(xué)對西藏玉龍斑巖銅礦成礦年齡的制約. 巖石學(xué)報, 22: 1009-1016.

    侯增謙, 莫宣學(xué), 楊志明, 王安建, 潘桂棠, 曲曉明, 聶鳳軍. 2006a. 青藏高原碰撞造山帶成礦作用: 構(gòu)造背景、時空分布和主要類型. 中國地質(zhì), 33(2): 348-359.

    侯增謙, 潘桂棠, 王安建, 莫宣學(xué), 田世洪, 孫曉明, 丁林, 王二七, 高永豐, 謝玉玲, 曾普勝, 秦克章, 許繼峰, 曲曉明, 楊志明, 楊竹森, 費紅彩, 孟祥金,李振清. 2006b. 青藏高原碰撞造山帶: II. 晚碰撞轉(zhuǎn)換成礦作用. 礦床地質(zhì), 25 (5): 521- 543.

    侯增謙, 潘小菲, 楊志明, 曲曉明. 2007. 初論大陸環(huán)境斑巖銅礦. 現(xiàn)代地質(zhì), 21(2): 332-351.

    姜耀輝, 蔣少涌, 凌洪飛, 戴寶章. 2006. 陸-陸碰撞造山環(huán)境下含銅斑巖巖石成因: 以藏東玉龍斑巖銅礦帶為例. 巖石學(xué)報, 22: 697-706.

    李蔭清. 1984. 玉龍某些中酸性侵入巖和火山雜巖中的熔融和流體包裹體. 中國地質(zhì)科學(xué)院院報, 85-106.

    李蔭清. 1985. 我國幾個重要(主要)斑巖銅(鉬)礦床的包裹體演化機制及成礦pH條件. 礦床地質(zhì), 4(2): 51-60.

    梁華英. 2002. 青藏高原東南緣斑巖銅礦成巖成礦研究取得新進展. 礦床地質(zhì), 21(4): 365.

    梁華英, 莫濟海, 孫衛(wèi)東, 喻亨祥, 張玉泉, Charllote M A.2008. 藏東玉龍超大型斑巖銅礦床成巖成礦系統(tǒng)時間跨度分析. 巖石學(xué)報, 24(10): 2352-2358.

    梁華英, 莫濟海, 孫衛(wèi)東, 張玉泉, 曾提, 胡光黔, Charllote M A. 2009. 玉龍銅礦帶馬拉松多斑巖體巖石學(xué)及成巖成礦系統(tǒng)年代學(xué)分析. 巖石學(xué)報, 25(2): 385-392.

    梁華英, 孫衛(wèi)東, 喻亨祥, 謝應(yīng)雯, 莫濟海, 張玉泉. 2006.西藏東緣玉龍斑巖銅礦帶含礦巖體時代及斑巖銅金礦床形成研究. 礦床地質(zhì), 25: 415-418.

    劉穎, 劉海臣, 李獻華. 1996. 用ICP-MS準(zhǔn)確測定巖石樣品40余種微量元素. 地球化學(xué), 25(6): 552-558.

    馬鴻文. 1989. 論藏東玉龍斑巖銅礦帶巖漿侵入時代. 地球化學(xué), (3): 210-216.

    馬鴻文. 1990. 西藏玉龍斑巖銅礦帶花崗巖類與成礦. 武漢: 中國地質(zhì)大學(xué)出版社: 1-158.

    孟祥金. 2004. 西藏碰撞造山帶岡底斯中新世斑巖銅礦成礦作用研究. 北京: 中國地質(zhì)科學(xué)院博士論文: 1-104.

    曲曉明, 侯增謙, 黃衛(wèi). 2001. 岡底斯斑巖銅成礦帶: 西藏第二條“玉龍”斑巖銅礦帶? 礦床地質(zhì), 20(4): 355-366.

    芮宗瑤, 黃崇軻, 齊國明. 1984. 中國斑巖銅(鉬)礦. 北京:地質(zhì)出版社: 1-300.

    芮宗瑤, 張立生, 陳振宇, 王龍生, 劉玉琳, 王義天. 2004.斑巖銅礦的源巖或源區(qū)探討. 巖石學(xué)報, 20(2): 229-238.

    唐菊興, 張麗, 李志軍, 陳建平, 黃衛(wèi), 王乾. 2006. 西藏玉龍銅礦床——鼻狀構(gòu)造圈閉控制的特大型礦床.礦床地質(zhì), 25(6): 654-662.

    唐仁鯉, 羅懷松, 李蔭清. 1995. 西藏玉龍銅(鉬)礦帶地質(zhì).北京: 地質(zhì)出版社: 76-96.

    王保弟, 陳陵康, 許繼峰, 劉鴻飛, 陳建林, 康志強. 2010.拉薩地塊麻江地區(qū)具有“超鉀質(zhì)”成分的鉀質(zhì)火山巖的識別及成因. 巖石學(xué)報, 27(6): 1662-1674.

    王建, 李建平, 王江海. 2003. 滇西大理-劍川地區(qū)鉀玄質(zhì)巖漿作用: 后碰撞走滑拉伸環(huán)境島弧型巖漿作用的地球化學(xué)研究. 巖石學(xué)報, 19(1): 61-69.

    王強, 唐功建, 賈曉輝, 資鋒, 姜子琦, 許繼峰, 趙振華. 2008. 埃達克質(zhì)巖的金屬成礦作用. 高校地質(zhì)學(xué)報, 14: 350-364.

    王強, 趙振華, 許繼峰, 白正華, 王建新, 劉成新. 2004.鄂東南銅山口、殷祖埃達克質(zhì)(adakite)的侵入巖地球化學(xué)特征對比: (拆沉)下地殼熔融與斑巖銅礦的成因.巖石學(xué)報, 20(2): 351-360.

    王增, 申屠保涌, 丁朝建. 1995. 藏東花崗巖類及其成礦作用. 成都: 西南交通大學(xué)出版社: 150-151.

    王治華, 郭曉東, 陳祥, 葛良勝, 鄒依林. 2010. 云南祥云馬廠箐富堿斑巖體的地球化學(xué)特征及其形成的構(gòu)造環(huán)境. 地質(zhì)評論, 56(1): 125-135.

    伍靜, 梁華英, 莫濟海, 張玉泉, 胡光黔. 2011. 玉龍斑巖銅礦帶莽總含礦斑巖體巖石學(xué)特征及鋯石U-Pb年齡研究. 大地構(gòu)造與成礦學(xué), 35(2): 300-306.

    曾普勝, 侯增謙, 高永峰, 杜安道. 2006. 印度-亞洲碰撞帶東段喜馬拉雅期銅-鉬-金礦床Re-Os年齡及成礦作用. 地質(zhì)論評, 52(1): 72-84.

    張旗, 秦克章, 王元龍, 張福勤, 劉紅濤, 王焰. 2004. 加強埃達克巖研究, 開創(chuàng)中國Cu, Au等找礦工作的新局面. 巖石學(xué)報, 20(2): 195-204.

    張玉泉, 謝應(yīng)雯, 李獻華, 邱華寧, 趙振華, 梁華英, 鐘孫霖. 2000. 青藏高原東部鉀玄巖漿巖同位素特征:巖石成因及構(gòu)造意義. 中國科學(xué)(D輯), 30(5): 493-498.

    張玉泉, 謝應(yīng)雯, 梁華英, 邱華寧, 李獻華, 鐘孫霖. 1998.藏東玉龍銅礦帶含礦斑巖及成巖系列. 地球化學(xué), 27(3): 236-243.

    張玉泉, 鐘孫霖. 1997. 藏東玉龍銅礦帶含礦斑巖演化與成礦關(guān)系. 西藏地質(zhì), 18: 73-86.

    趙欣, 喻學(xué)惠, 莫宣學(xué), 張瑾, 呂伯西. 2004. 滇西新生代富堿斑巖及其深源包體的巖石學(xué)和地球化學(xué)特征.現(xiàn)代地質(zhì), 18(2): 217-228.

    鐘大賚, 丁林, 劉福田, 劉建華, 張進江, 季建清, 陳輝. 2000. 造山帶巖石層多向?qū)蛹軜?gòu)造及其對新生代巖漿活動制約——以三江及鄰區(qū)為例. 中國科學(xué)(D輯), 30(1): 1-8.

    周宜吉. 1980. 玉龍礦帶中斑巖銅礦的控礦因素及其成因探討. 地質(zhì)論評, 26(4): 357-361.

    周宜吉. 1985. 試論玉龍斑巖銅礦帶內(nèi)礦床的蝕變和礦化分帶. 礦床地質(zhì), 4(2): 23-30.

    朱訓(xùn), 黃崇輸, 芮宗瑤. 1983. 德興斑巖銅礦. 北京: 地質(zhì)出版社: 1-314.

    Defant M J and Drummond M S. 1990. Derivation of some modern arcmagmas by melting of young subducted lithosphere. Nature, 347: 662-665.

    Defant M J and Drummond M S. 1993. Mount St. Helens: Potential example of the partial melting of the subducted lithosphere in a volcanic arc. Geology, 21(6): 547-550.

    Defant M J and Kepezhinskas P. 2001. Adakites: A review of slab melting over the past decade and the case for a slab-melt component in arcs. EOS Transactions, 82: 65, 68-69.

    Drummond M S and Defant M J. 1990. A model for trondhjenite-tonalite-dactite genesis and crustal growth via slab melting: Archean to modern composition. Journal of Geophysical Research, 95: 503-521.

    Foley S. 1992. Vein-plus-wall-rock melting mechanisms in the lithosphere and the origin of potassic alkaline magmas. Lithos, 28(3-6): 435-453.

    Gao Y F, Hou Z Q, Kamber B S, Wei R H, Meng X J and Zhao R S. 2007. Adakite-like porphyries from the southern Tibetan continental collision zones: Evidence for slab melt metasomatism. Contributions to Mineralogy and Petrology, 153(1): 105-120.

    Gao Y F, Yang Z S, Santosh M, Hou Z Q, Wei R H and Tian S L. 2010. Adakitic rocks from slab melt-modified mantle sources in the continental collision zone of southern Tibet. Lithos, 119(3-4): 651-663.

    Gill J B. 1981. Orogenic andesites and plate tectonic. New York: Springer-Verlag: 390.

    Guo Z F, Wilson M and Liu J Q. 2007. Post-collisional adakites in south Tibet: Products of partial melting of subduction-modified lower crust. Lithos, 96(1-2): 205-224.

    Hou Z Q, Gao Y F, Qu X M, Rui Z Y and Mo X X. 2004. Origin of adakitic intrusives generated during mid-Miocene east-west extension in southern Tibet. Earth and Planetary Science Letters, 220: 139-155.

    Hou Z Q, Ma H W, Zaw K, Zhang Y Q, Wang M J, Wang Z, Pan G T and Tang R L. 2003. The Himalayan Yulong porphyry copper belt: Product of large-scale strike-slip faulting in eastern Tibet. Economic Geology, 98(1): 125-145.

    Jiang Y H, Jiang S Y, Ling H F and Dai B Z. 2006. Low-degree melting of a metasomatized lithospheric mantle for the origin of Cenozoic Yulong monzogranite-porphyry, east Tibet: Geochemical and Sr-Nd-Pb-Hf isotopic constraints. Earth and Planetary Science Letters, 241: 617-633.

    Kay R W. 1978. Aleutian magnesium andesites: Melts from subducted Pacific oceanic crust. Journal of Volcanology and Geothermal Research, 4(1-2): 117-132.

    Kay S M and Mpodozis C. 2001. Central Andean ore deposits linked to evolved shallow seduction systems and thickening crust. GSA Today, 11: 4-9.

    Kay S M, Ramos V A and Marquez M. 1993. Evidence in Cerro Pampa volcanic rocks for slab-melting prior to ridge-trench collision in southern South Ameriea. The Journal of Geology, 101(6): 703-714.

    Le Bas M J, Le Maitre R W, Streckeisen A, Zanettin B and IUGS Subcommission on the Systematics of Igneous Rocks. 1986. A chemical classification of volcanic rocks based on the total alkali-silica diagram. Journal of Petrology, 27(3): 745-750.

    Liang H Y, Campbell I H, Allen C, Sun W D, Liu C Q, Yu H X, Xie Y W and Zhang Y Q. 2006a. Zircon Ce4+/Ce3+ratios and ages for Yulong ore-bearing porphyries in eastern Tibet. Mineralium Deposita, 41(2): 152-159.

    Liang H Y, Yu H X, Mo C H, Zhang Y Q and Xie Y W. 2006b. Zircon LA-ICP-MS U-Pb age Ce4+/Ce3+ratios and the geochemical features of the Machangqing complex associated the copper deposit. Chinese Journal of Geochemistry, 25(3): 223-229.

    Liang H Y, Zhang Y Q, Xie Y W, Lin W, Campbell I H and Yu H X. 2005. Geochronological and geochemical study on the Yulong porphyry copper ore belt in eastern Tibet, China // Mineral Deposit Research: Meeting the Global Challenge, 1: 1235-1237.

    Martin H. 1986. Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone magmas. Geology, 14(9): 753-756.

    Middlemost E A K. 1994. Naming materials in the magma/igneous rock system. Earth-Science Reviews, 37(3-4): 215-224.

    Miller C, Schuster R, Kl?tzli U, Frank W and Purtscheller F. 1999. Post-collisional potassic and ultrapotassic magmatism in SW Tibet: Geochemical and Sr-Nd-Pb-O isotopic constraints for mantle source characteristics and petrogenesis. Journal of Petrology, 40(9): 1399-1424.

    Nelson D R. 1992. Isotopic characteristics of potassic rocks: Evidence for the involvement of subducted sediments in magma genesis. Lithos, 28(3-6): 403-420.

    Oyarzun R, Márquez A, Lillo J, López I and Rivera S. 2001. Giant versus small porphyry copper deposits of Cenozoic age in northern Chile: Adakitic versus normal calc-alkaline magmatism. Mineralium Deposita, 36: 794-798.

    Petford N and Atherton M. 1996. Na-rich partial melts from newly underplated basaltic crust: The Cordillera Blanca Batholith, Peru. Journal of Petrology, 37(6): 1491-1521.

    Plank T and Langmuir C H. 1998. The chemical composition of subducting sediment and its consequences for the crust and mantle. Chemical Geology, 145(3-4): 325-394.

    Qu X M, Hou Z Q and Li Y G. 2004. Melt components derived from a subducted slab in late orogenic ore-bearing porphyries in the Gangdese copper belt,southern Tibetan plateau. Lithos, 74: 131-148.

    Richards J P. 2003. Tectono-magmatic precursors for porphyry Cu-(Mo-Au) deposit formation. Economic Geology, 98: 1515-1533.

    Stern C R and Kilian R. 1996. Role of the subducted slab, mantle wedge and continental crust in the generation of adakites from the Andean Austral Volcanic Zone. Contributions to Mineralogy and Petrology, 123(3): 263-281.

    Sun H J, Deng W M and Zhang Y Q. 2001. Petrogenesis of Cenozoic potassic volcanic rocks in the Nangqen Basin. Acta Geologica Sinca, 75(1): 27-40.

    Sun S S and McDonough W F. 1989. Chemical and isotopic systematic of oceanic basalts: Implication for mantle composition and processes // Saunders A D and Norry M J. Magmatism in Oceanic Basins. Geological Society, London, Special Publications, 42(1): 313-345.

    Turner S, Hawkesworth C, Liu J Q, Rogers N, Kelley S and Calsteren P. 1993. Timing of Tibetan uplift constrained by analysis of volcanic rocks. Nature, 364: 50-54.

    Wang J H, Yin A, Harrison T M, Grove M, Zhang Y Q and Xie G H. 2001. A tectonic model for Cenozoic igneous activities in the eastern Indo-Asian collision zone. Earth and Planetary Science Letters, 188: 123-133.

    Zindle A and Hart S R. 1986. Chemical geodynamics. Annual Review of Earth and Planetary Sciences, 14: 493-573.

    Geochemical Characteristics and Metallogenic Mechanism of the Porphyry Cu-Mo Deposits in the Yulong Ore Belt, Eastern Tibet: A Case Study of the Yulong and Duoxiasongduo Porphyries

    WU Weizhong1,2, XIA Bin1,3, ZHANG Yuquan1, DONG Binghua1,2and XIA Zhongxi3
    (1. Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, Guangdong, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. School of Marine Sciences, Sun Yat-Sen University, Guangzhou 510275, Guangdong, China)

    The Yulong porphyry copper belt, located in the middle of the Sanjiang tectonic zone of the eastern Tethys-Himalayan tectonic domain, is an ideal area to study the petrogenesis and mineralization of the Himalayan alkali-rich porphyries. This paper takes the Yulong and Duoxiasongduo as the study areas, presenting detailed dating on the porphyries and associated Cu-Mo mineralization as well as their geological and geochemical characteristics. Combined with Sr-Nd-Pb and S-H-O isotopic analysis, we discussed the difference of ore-bearing magma source, tectonic controls on mineralization, and dynamic mechanism between the Yulong and Gangdese ore-bearing porphyries. The results show that the ore-bearing porphyries which are of potassic alkali-granite and characterized by rich alkali, high potassium and high aluminum. In the Yulong area, the rocks mainly consist of quartz monzonite porphyry in the early stage, the granodiorite and quartz monzonite porphyry in the late stage. While in the Duoxiasongduo area, the rocks mainly consist of monzogranite and granodiorite porphyries, and less K-feldspar granite porphyry. Both of them are rich in light rare earth elements (LREE), large ion lithophile elements (LILE, e.g., Rb, Th, U and Pb), and depleted relatively in HFSE (e.g., Nb, Ta and Yb) without significant Eu and Ce negative anomalies. The slightly higher than modern primitive mantle initial87Sr/86Sr ratios, lower εNd(t) values, relatively uniform206Pb/204Pb ratios of the porphyries in the studying area, indicate that they were probably generated from partial melting of the garnet amphibolite and/or amphibolite eclogite, followed by a low degree of crystal fractionation and crustal contamination during the magma ascending and intrusion. We suggest that the formation of the Yulong large porphyry copper deposit belt is related to the Toba-Mangkang sinistral strike-slip fault (TBF) in the Red River-Ailaoshan fracture system caused by the India-Asia continental collision in middle Eocene of early Tertiary (37~38 Ma).

    alkali-rich porphyry; geochemical characteristics; metallogenic mechanism; Yulong copper belt; eastern Tibet

    P611; P595

    A

    1001-1552(2013)03-0440-015

    2013-03-13; 改回日期: 2013-04-10

    項目資助: 國家自然科學(xué)基金(批準(zhǔn)號: 4053019)、廣州市科技計劃項目(編號: 2012420004205008)和廣東省高校海洋石油勘探開發(fā)重點實驗基金(42000-3211101)聯(lián)合資助。

    吳偉中(1984-), 男, 博士研究生, 構(gòu)造地質(zhì)學(xué)專業(yè)。Email: 1604222@163.com

    猜你喜歡
    玉泉含礦玉龍
    從此蜀道不再難
    國企管理(2024年6期)2024-01-01 00:00:00
    紅山玉龍
    ——中華第一玉龍
    家教世界(2023年28期)2023-11-14 10:08:18
    紅山玉龍
    生日快樂,我的玉泉
    軍事文摘(2023年14期)2023-07-28 08:39:28
    喜報!玉泉橋牌選手獲獎啦
    軍事文摘(2021年22期)2022-01-18 06:22:58
    紅船“開進”了玉泉小學(xué)
    軍事文摘(2021年16期)2021-11-05 08:49:40
    斑巖型礦床含礦斑巖與非含礦斑巖鑒定特征綜述
    巖型礦床含礦斑巖與非含礦斑巖鑒定特征綜述
    玉泉豐收采摘節(jié)
    軍事文摘(2018年24期)2018-12-26 00:58:28
    新疆卡特巴阿蘇金礦床含礦巖石及圍巖地球化學(xué)特征與構(gòu)造環(huán)境簡析
    成人高潮视频无遮挡免费网站| 99热精品在线国产| 国产av麻豆久久久久久久| 精品福利观看| 国产高清激情床上av| 偷拍熟女少妇极品色| 在线免费十八禁| 女的被弄到高潮叫床怎么办| 小蜜桃在线观看免费完整版高清| 亚洲欧美精品自产自拍| 麻豆久久精品国产亚洲av| 亚洲成av人片在线播放无| 亚洲中文字幕日韩| 久久久成人免费电影| 亚洲av第一区精品v没综合| 狂野欧美白嫩少妇大欣赏| 真实男女啪啪啪动态图| 国产高清不卡午夜福利| 国产蜜桃级精品一区二区三区| 精品久久久久久久人妻蜜臀av| 精品不卡国产一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 99久久久亚洲精品蜜臀av| 亚洲av.av天堂| 69人妻影院| 久久精品国产亚洲网站| 看黄色毛片网站| 国产一级毛片七仙女欲春2| 热99在线观看视频| 69av精品久久久久久| 国产成人影院久久av| 看十八女毛片水多多多| 噜噜噜噜噜久久久久久91| 亚洲一级一片aⅴ在线观看| 午夜影院日韩av| 成人毛片a级毛片在线播放| 精品久久久久久久久久免费视频| 三级毛片av免费| 熟女人妻精品中文字幕| 别揉我奶头~嗯~啊~动态视频| 精品乱码久久久久久99久播| 丰满人妻一区二区三区视频av| 亚洲av免费在线观看| 又黄又爽又免费观看的视频| 一卡2卡三卡四卡精品乱码亚洲| 成人特级av手机在线观看| 乱人视频在线观看| 亚洲国产日韩欧美精品在线观看| 直男gayav资源| 少妇的逼水好多| 久久久久国产网址| 91久久精品电影网| 国产精品av视频在线免费观看| 99热只有精品国产| 少妇裸体淫交视频免费看高清| 97热精品久久久久久| 男人和女人高潮做爰伦理| 国产成人影院久久av| 丰满人妻一区二区三区视频av| 日日撸夜夜添| 亚洲18禁久久av| 亚洲最大成人av| 综合色丁香网| 亚洲国产色片| 欧美3d第一页| 一级黄片播放器| 真实男女啪啪啪动态图| 国产单亲对白刺激| 久久精品国产99精品国产亚洲性色| 亚洲美女搞黄在线观看 | 久久久久久久久中文| 欧美日韩综合久久久久久| 精品一区二区三区人妻视频| 日本免费一区二区三区高清不卡| www日本黄色视频网| 色综合亚洲欧美另类图片| 成人三级黄色视频| 97超视频在线观看视频| 激情 狠狠 欧美| av视频在线观看入口| 一级毛片我不卡| 小说图片视频综合网站| 国产高清不卡午夜福利| 精品人妻视频免费看| 国产av不卡久久| 国产精品三级大全| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品电影| 草草在线视频免费看| 亚洲美女视频黄频| 亚州av有码| 免费电影在线观看免费观看| 国产精品人妻久久久影院| 国产一区二区激情短视频| 日韩欧美免费精品| 一个人观看的视频www高清免费观看| 亚洲欧美清纯卡通| 国产亚洲欧美98| 中文资源天堂在线| 国产午夜精品久久久久久一区二区三区 | 国产伦精品一区二区三区视频9| 色播亚洲综合网| 国产乱人偷精品视频| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产亚洲av香蕉五月| 男女视频在线观看网站免费| 联通29元200g的流量卡| 人人妻人人澡欧美一区二区| 国产成人a∨麻豆精品| 欧美高清性xxxxhd video| 婷婷精品国产亚洲av| 亚洲综合色惰| 非洲黑人性xxxx精品又粗又长| 小蜜桃在线观看免费完整版高清| 亚洲真实伦在线观看| 成人漫画全彩无遮挡| 成人鲁丝片一二三区免费| 黄色一级大片看看| 精品午夜福利在线看| 精品欧美国产一区二区三| 久久久欧美国产精品| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 精品久久久久久成人av| 成人av一区二区三区在线看| 少妇裸体淫交视频免费看高清| 亚洲人成网站在线播| 欧美最新免费一区二区三区| АⅤ资源中文在线天堂| 最近的中文字幕免费完整| 欧美高清成人免费视频www| 日本色播在线视频| 亚洲av免费高清在线观看| 啦啦啦韩国在线观看视频| 国产高潮美女av| 91久久精品国产一区二区三区| 亚洲成a人片在线一区二区| 日韩国内少妇激情av| 日韩精品有码人妻一区| 麻豆成人午夜福利视频| 欧美绝顶高潮抽搐喷水| 高清日韩中文字幕在线| 少妇熟女欧美另类| 国产午夜精品论理片| av国产免费在线观看| 久久韩国三级中文字幕| 一级毛片aaaaaa免费看小| 亚洲av中文av极速乱| av在线天堂中文字幕| 国产精品久久久久久久久免| 亚洲国产高清在线一区二区三| 麻豆国产97在线/欧美| 国产日本99.免费观看| 日本熟妇午夜| 波多野结衣高清作品| 国产麻豆成人av免费视频| 欧美日本视频| 久久久久久伊人网av| 亚洲av免费高清在线观看| 国产精品人妻久久久久久| 亚洲乱码一区二区免费版| 久久6这里有精品| 国产乱人偷精品视频| 成人av一区二区三区在线看| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 男人舔女人下体高潮全视频| 久久久色成人| 欧美不卡视频在线免费观看| 色哟哟·www| 国产69精品久久久久777片| 伦理电影大哥的女人| 国产精品免费一区二区三区在线| 国产视频内射| 亚洲第一区二区三区不卡| 嫩草影院精品99| 色av中文字幕| 悠悠久久av| 黄色视频,在线免费观看| 亚洲精品乱码久久久v下载方式| 亚洲精品一卡2卡三卡4卡5卡| 亚洲内射少妇av| 免费看日本二区| 成人亚洲欧美一区二区av| 亚洲av成人av| 99久国产av精品| 免费观看在线日韩| 日本爱情动作片www.在线观看 | 久久人人爽人人爽人人片va| 晚上一个人看的免费电影| av.在线天堂| 尤物成人国产欧美一区二区三区| 两个人的视频大全免费| 激情 狠狠 欧美| 亚洲国产精品久久男人天堂| 日韩精品青青久久久久久| 国产高清有码在线观看视频| 国产又黄又爽又无遮挡在线| 欧美又色又爽又黄视频| 久久九九热精品免费| 精品国内亚洲2022精品成人| 看非洲黑人一级黄片| 高清毛片免费观看视频网站| 欧美日本视频| 国产av不卡久久| 99久国产av精品| 色哟哟哟哟哟哟| 少妇裸体淫交视频免费看高清| 国产aⅴ精品一区二区三区波| 22中文网久久字幕| 岛国在线免费视频观看| 国产av麻豆久久久久久久| 在线观看一区二区三区| 乱系列少妇在线播放| 不卡视频在线观看欧美| 亚洲18禁久久av| 极品教师在线视频| 一区二区三区高清视频在线| 性插视频无遮挡在线免费观看| 99热精品在线国产| 高清日韩中文字幕在线| 欧美极品一区二区三区四区| 亚洲真实伦在线观看| 99riav亚洲国产免费| 午夜日韩欧美国产| 国产精品人妻久久久影院| 老司机午夜福利在线观看视频| 亚洲国产高清在线一区二区三| 亚洲18禁久久av| 国产真实乱freesex| 精品久久久噜噜| 日本-黄色视频高清免费观看| 免费看a级黄色片| 亚洲熟妇熟女久久| 18禁裸乳无遮挡免费网站照片| 国产av一区在线观看免费| 国产精品国产高清国产av| 男女做爰动态图高潮gif福利片| 欧美性猛交黑人性爽| 国产精品久久久久久av不卡| 色视频www国产| 久久久欧美国产精品| 亚洲av中文av极速乱| 高清日韩中文字幕在线| 69av精品久久久久久| 在线看三级毛片| 中文字幕av在线有码专区| 亚洲国产高清在线一区二区三| 在线免费观看的www视频| 舔av片在线| 免费看日本二区| 看片在线看免费视频| 久久精品91蜜桃| 国产一区二区在线观看日韩| 成年女人永久免费观看视频| 国产中年淑女户外野战色| 婷婷精品国产亚洲av| 嫩草影院精品99| 男女那种视频在线观看| 国产高清激情床上av| aaaaa片日本免费| 国产精品久久久久久久电影| 男女之事视频高清在线观看| 日韩一区二区视频免费看| 午夜a级毛片| 禁无遮挡网站| 真实男女啪啪啪动态图| 全区人妻精品视频| 国国产精品蜜臀av免费| 欧美中文日本在线观看视频| 人人妻人人澡人人爽人人夜夜 | 日本免费一区二区三区高清不卡| 男女做爰动态图高潮gif福利片| 国产精品久久视频播放| a级毛色黄片| 老司机福利观看| 亚洲欧美日韩高清在线视频| 一个人看视频在线观看www免费| 国产女主播在线喷水免费视频网站 | 1000部很黄的大片| av卡一久久| 午夜精品在线福利| 国产午夜精品论理片| 天堂√8在线中文| 精品久久久久久成人av| 91av网一区二区| 久久精品夜色国产| 日本成人三级电影网站| 日日摸夜夜添夜夜爱| 亚洲人成网站在线观看播放| 精品午夜福利在线看| 高清毛片免费看| 久久人人精品亚洲av| 国产精品一区www在线观看| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频 | 久久精品国产亚洲av香蕉五月| 国产日本99.免费观看| 身体一侧抽搐| 一本久久中文字幕| 精品少妇黑人巨大在线播放 | 亚洲激情五月婷婷啪啪| 成人性生交大片免费视频hd| 女同久久另类99精品国产91| 欧美在线一区亚洲| 国产亚洲精品综合一区在线观看| 国产麻豆成人av免费视频| 黄色配什么色好看| 99九九线精品视频在线观看视频| 18+在线观看网站| 久久精品影院6| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 最新在线观看一区二区三区| 此物有八面人人有两片| 热99re8久久精品国产| 欧美高清成人免费视频www| 日本精品一区二区三区蜜桃| av在线老鸭窝| 99热只有精品国产| av在线老鸭窝| 黑人高潮一二区| 99久久成人亚洲精品观看| 精品久久久久久久久亚洲| 国产视频内射| 久久人人精品亚洲av| 最近的中文字幕免费完整| 人妻夜夜爽99麻豆av| 免费av不卡在线播放| 免费看美女性在线毛片视频| 亚洲无线观看免费| 嫩草影院精品99| a级一级毛片免费在线观看| 精品午夜福利在线看| 中文字幕人妻熟人妻熟丝袜美| 最新中文字幕久久久久| 精品午夜福利视频在线观看一区| 欧美日韩在线观看h| 99热6这里只有精品| ponron亚洲| 国产精品嫩草影院av在线观看| 在线观看免费视频日本深夜| 亚洲va在线va天堂va国产| 久久鲁丝午夜福利片| 日韩欧美免费精品| 伊人久久精品亚洲午夜| 最近视频中文字幕2019在线8| 天堂av国产一区二区熟女人妻| 免费高清视频大片| 亚洲欧美日韩高清在线视频| 久久人人爽人人爽人人片va| 成人欧美大片| 变态另类丝袜制服| 青春草视频在线免费观看| 国产精品免费一区二区三区在线| 国产精品一区二区免费欧美| 久久久久久伊人网av| 久久精品国产鲁丝片午夜精品| 婷婷色综合大香蕉| 精品久久久久久成人av| 久久婷婷人人爽人人干人人爱| 极品教师在线视频| 久久精品久久久久久噜噜老黄 | 最近2019中文字幕mv第一页| 国产成人精品久久久久久| 日韩欧美精品免费久久| 黄色日韩在线| 欧美成人一区二区免费高清观看| 在线观看66精品国产| 91久久精品国产一区二区成人| 国产午夜精品论理片| av在线天堂中文字幕| 亚洲精华国产精华液的使用体验 | 午夜激情福利司机影院| 久久久午夜欧美精品| 精品人妻熟女av久视频| 中文字幕av成人在线电影| 国产精品一区www在线观看| 日本黄大片高清| 日韩亚洲欧美综合| 成人av一区二区三区在线看| 51国产日韩欧美| 人妻夜夜爽99麻豆av| 1024手机看黄色片| 22中文网久久字幕| 欧美色欧美亚洲另类二区| 久久欧美精品欧美久久欧美| 欧美色欧美亚洲另类二区| 亚洲精品一区av在线观看| 国产精品嫩草影院av在线观看| 国产成人精品久久久久久| 国产精品嫩草影院av在线观看| 国产精品一区二区性色av| 好男人在线观看高清免费视频| 日本-黄色视频高清免费观看| 99久久精品国产国产毛片| 午夜福利在线观看免费完整高清在 | 日本欧美国产在线视频| 色哟哟哟哟哟哟| 国产淫片久久久久久久久| 国产精品一区二区性色av| 成人午夜高清在线视频| 丝袜喷水一区| 国产午夜精品久久久久久一区二区三区 | 久久精品综合一区二区三区| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品电影| 国产在视频线在精品| 丰满的人妻完整版| 黄色欧美视频在线观看| 99在线人妻在线中文字幕| 99热精品在线国产| 久久久精品欧美日韩精品| 国产伦在线观看视频一区| 中国国产av一级| 欧美激情国产日韩精品一区| 色播亚洲综合网| 欧美最新免费一区二区三区| 99久久精品国产国产毛片| 麻豆久久精品国产亚洲av| 亚洲成人久久爱视频| 人人妻,人人澡人人爽秒播| 日本黄色片子视频| 婷婷色综合大香蕉| 国产在视频线在精品| 日韩三级伦理在线观看| 99精品在免费线老司机午夜| 成人特级av手机在线观看| 国产精品人妻久久久久久| 免费看美女性在线毛片视频| 婷婷亚洲欧美| 美女cb高潮喷水在线观看| 亚洲中文字幕一区二区三区有码在线看| 午夜精品在线福利| 久久久午夜欧美精品| 欧美中文日本在线观看视频| 久久婷婷人人爽人人干人人爱| 老女人水多毛片| 深夜精品福利| av专区在线播放| 精品久久久久久久久av| 国内精品宾馆在线| 日本成人三级电影网站| 在线免费十八禁| 日日啪夜夜撸| 女人被狂操c到高潮| 搡老岳熟女国产| 一卡2卡三卡四卡精品乱码亚洲| 色综合亚洲欧美另类图片| 看非洲黑人一级黄片| 99热这里只有精品一区| 99久久中文字幕三级久久日本| 中文字幕久久专区| 综合色av麻豆| 又黄又爽又免费观看的视频| 97热精品久久久久久| 日本黄大片高清| 午夜激情欧美在线| 国产一区二区在线观看日韩| 又黄又爽又刺激的免费视频.| 九色成人免费人妻av| 亚洲中文日韩欧美视频| 成人午夜高清在线视频| 免费观看在线日韩| 99热全是精品| 美女免费视频网站| 丝袜喷水一区| 亚洲最大成人av| 国产毛片a区久久久久| 18禁在线无遮挡免费观看视频 | 插阴视频在线观看视频| 悠悠久久av| 国产精品一二三区在线看| 亚洲真实伦在线观看| 亚洲av电影不卡..在线观看| 久久中文看片网| 啦啦啦韩国在线观看视频| 欧美性感艳星| 亚洲人与动物交配视频| 国产 一区精品| 高清毛片免费看| 波多野结衣高清作品| 在线播放无遮挡| av在线播放精品| 亚洲中文字幕日韩| 一级黄片播放器| 一a级毛片在线观看| 亚洲精品乱码久久久v下载方式| 亚洲av美国av| 午夜日韩欧美国产| av卡一久久| 两个人的视频大全免费| 男人狂女人下面高潮的视频| 国产精品,欧美在线| 波多野结衣巨乳人妻| 成人午夜高清在线视频| 日本黄大片高清| 最新在线观看一区二区三区| 久久久久久久久大av| 美女被艹到高潮喷水动态| 天美传媒精品一区二区| 日本黄色视频三级网站网址| 亚洲成a人片在线一区二区| 热99在线观看视频| 免费一级毛片在线播放高清视频| 综合色丁香网| 国语自产精品视频在线第100页| 给我免费播放毛片高清在线观看| 亚洲精品一区av在线观看| 欧美色欧美亚洲另类二区| 亚洲精品一区av在线观看| 干丝袜人妻中文字幕| 色吧在线观看| 搡老熟女国产l中国老女人| 日韩av在线大香蕉| 国产高清不卡午夜福利| 美女 人体艺术 gogo| 精品人妻视频免费看| 十八禁国产超污无遮挡网站| 91久久精品国产一区二区三区| 99久久九九国产精品国产免费| 最近在线观看免费完整版| 久久国产乱子免费精品| 最近最新中文字幕大全电影3| 噜噜噜噜噜久久久久久91| .国产精品久久| 一进一出抽搐动态| 亚洲美女视频黄频| 赤兔流量卡办理| 夜夜爽天天搞| 国产精品综合久久久久久久免费| 久久草成人影院| 波多野结衣高清作品| 不卡视频在线观看欧美| 老司机影院成人| 亚洲乱码一区二区免费版| 国产综合懂色| 免费观看人在逋| 在线天堂最新版资源| 夜夜夜夜夜久久久久| 美女黄网站色视频| 国产精品99久久久久久久久| 日本欧美国产在线视频| 亚洲成人中文字幕在线播放| 精品久久久久久久久久免费视频| 日韩成人av中文字幕在线观看 | 国产69精品久久久久777片| 美女黄网站色视频| 熟妇人妻久久中文字幕3abv| 男人和女人高潮做爰伦理| 日本色播在线视频| 免费观看精品视频网站| 欧美zozozo另类| 天天躁夜夜躁狠狠久久av| 午夜精品国产一区二区电影 | 免费大片18禁| 高清午夜精品一区二区三区 | 国产成人一区二区在线| 给我免费播放毛片高清在线观看| 搡老岳熟女国产| 成人一区二区视频在线观看| av卡一久久| 国产精品久久久久久亚洲av鲁大| 午夜久久久久精精品| 嫩草影院精品99| 成人特级黄色片久久久久久久| 欧美中文日本在线观看视频| 久久精品影院6| 亚洲av五月六月丁香网| 久99久视频精品免费| 国产蜜桃级精品一区二区三区| 在线播放无遮挡| 嫩草影院新地址| 亚洲av五月六月丁香网| 亚洲va在线va天堂va国产| 寂寞人妻少妇视频99o| 国产淫片久久久久久久久| 国产女主播在线喷水免费视频网站 | 91狼人影院| 国产精华一区二区三区| 亚洲成人中文字幕在线播放| 国产伦精品一区二区三区四那| 男女视频在线观看网站免费| av天堂在线播放| 精品久久久久久久久久免费视频| 亚洲av熟女| 国内少妇人妻偷人精品xxx网站| 成年av动漫网址| 嫩草影院精品99| 午夜福利在线观看免费完整高清在 | 尤物成人国产欧美一区二区三区| 俄罗斯特黄特色一大片| 免费观看在线日韩| 午夜福利成人在线免费观看| 91久久精品国产一区二区三区| 国产精品久久久久久久电影| 极品教师在线视频| 色哟哟·www| 国产视频一区二区在线看| 亚洲精品456在线播放app| 色哟哟·www| 插逼视频在线观看| 晚上一个人看的免费电影| 91久久精品国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 老司机影院成人| 国产精品av视频在线免费观看| 久久久久久久久中文| 国产亚洲精品久久久com| 一a级毛片在线观看| 国产精品不卡视频一区二区| 老司机影院成人| 久久久成人免费电影| 亚洲性久久影院| 最近在线观看免费完整版|