• <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)境簡析
    久久久久久久精品吃奶| 伦理电影免费视频| 成人三级黄色视频| 很黄的视频免费| 大型黄色视频在线免费观看| 悠悠久久av| 亚洲精品在线美女| 搡老妇女老女人老熟妇| 国产伦在线观看视频一区| 欧美性猛交╳xxx乱大交人| 香蕉久久夜色| 高清毛片免费观看视频网站| 国产精华一区二区三区| 国产精品二区激情视频| 亚洲欧洲精品一区二区精品久久久| 中文字幕人成人乱码亚洲影| 国产激情久久老熟女| 午夜福利一区二区在线看| www.熟女人妻精品国产| 精品乱码久久久久久99久播| 精品久久久久久久久久免费视频| 99国产极品粉嫩在线观看| 国产成+人综合+亚洲专区| 亚洲欧美一区二区三区黑人| 手机成人av网站| 一个人观看的视频www高清免费观看 | 动漫黄色视频在线观看| 久久九九热精品免费| 日韩中文字幕欧美一区二区| 成年人黄色毛片网站| 国内少妇人妻偷人精品xxx网站 | 一级毛片高清免费大全| 亚洲人成77777在线视频| 国产伦在线观看视频一区| 日韩av在线大香蕉| 国产视频一区二区在线看| 国产av又大| 91字幕亚洲| 97超级碰碰碰精品色视频在线观看| 搡老岳熟女国产| 一卡2卡三卡四卡精品乱码亚洲| 国产免费av片在线观看野外av| 午夜免费鲁丝| 精品久久久久久久毛片微露脸| 久久精品亚洲精品国产色婷小说| 国产精品久久久久久人妻精品电影| 麻豆成人av在线观看| 久久久国产成人免费| 午夜影院日韩av| √禁漫天堂资源中文www| 18禁国产床啪视频网站| 成人欧美大片| 日本免费a在线| 国产乱人伦免费视频| 香蕉国产在线看| 人人妻人人澡欧美一区二区| 国产亚洲欧美精品永久| 18禁国产床啪视频网站| 成人精品一区二区免费| 久久精品国产亚洲av高清一级| 国产精品亚洲一级av第二区| 欧美日韩福利视频一区二区| 97人妻精品一区二区三区麻豆 | 99国产精品99久久久久| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 99久久国产精品久久久| 美女午夜性视频免费| 黑人欧美特级aaaaaa片| 精品久久久久久久毛片微露脸| 曰老女人黄片| 麻豆成人av在线观看| 999久久久国产精品视频| 国产免费av片在线观看野外av| 国产精品爽爽va在线观看网站 | 黑人欧美特级aaaaaa片| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文字幕日韩| 夜夜夜夜夜久久久久| 69av精品久久久久久| 成人手机av| 国产成人一区二区三区免费视频网站| 巨乳人妻的诱惑在线观看| 91老司机精品| 性色av乱码一区二区三区2| 久久中文字幕一级| 国产亚洲精品久久久久5区| 亚洲精品国产精品久久久不卡| 侵犯人妻中文字幕一二三四区| 看黄色毛片网站| 国产片内射在线| 精品午夜福利视频在线观看一区| 悠悠久久av| x7x7x7水蜜桃| 亚洲久久久国产精品| 久久久久国产精品人妻aⅴ院| 人人妻人人澡人人看| 欧美激情 高清一区二区三区| 又黄又粗又硬又大视频| 国产精华一区二区三区| 午夜福利欧美成人| 91九色精品人成在线观看| 久久久久国产精品人妻aⅴ院| 国产亚洲精品久久久久久毛片| 欧美精品啪啪一区二区三区| 夜夜躁狠狠躁天天躁| 精品日产1卡2卡| 很黄的视频免费| 久久精品国产99精品国产亚洲性色| 亚洲精品国产精品久久久不卡| 日日干狠狠操夜夜爽| 免费在线观看亚洲国产| 日韩大尺度精品在线看网址| 丝袜人妻中文字幕| 午夜a级毛片| 日韩国内少妇激情av| 日本撒尿小便嘘嘘汇集6| 男女床上黄色一级片免费看| 免费无遮挡裸体视频| 国产成人欧美在线观看| 色综合站精品国产| 亚洲精品在线美女| 国产精品永久免费网站| 久久精品国产亚洲av高清一级| 久久久久亚洲av毛片大全| 成人一区二区视频在线观看| 岛国视频午夜一区免费看| 啦啦啦观看免费观看视频高清| 91成年电影在线观看| 国产高清视频在线播放一区| 国产av一区在线观看免费| cao死你这个sao货| 波多野结衣av一区二区av| 一级毛片高清免费大全| 麻豆成人午夜福利视频| 久久亚洲精品不卡| 欧美绝顶高潮抽搐喷水| 桃色一区二区三区在线观看| 女性生殖器流出的白浆| 悠悠久久av| 91老司机精品| 男人的好看免费观看在线视频 | 又紧又爽又黄一区二区| 熟女少妇亚洲综合色aaa.| 男女之事视频高清在线观看| 男女午夜视频在线观看| 中文字幕高清在线视频| 一级作爱视频免费观看| 亚洲成av片中文字幕在线观看| 久久国产精品影院| 美女大奶头视频| 亚洲精品粉嫩美女一区| 国产精品 国内视频| 国产又色又爽无遮挡免费看| 亚洲va日本ⅴa欧美va伊人久久| 首页视频小说图片口味搜索| 俄罗斯特黄特色一大片| 亚洲美女黄片视频| 又黄又粗又硬又大视频| 首页视频小说图片口味搜索| 欧美激情高清一区二区三区| 波多野结衣av一区二区av| 亚洲自偷自拍图片 自拍| 欧美丝袜亚洲另类 | 日本黄色视频三级网站网址| 一本精品99久久精品77| 嫁个100分男人电影在线观看| 欧美午夜高清在线| 亚洲中文字幕日韩| 高潮久久久久久久久久久不卡| 亚洲精品在线观看二区| 精品国产乱子伦一区二区三区| tocl精华| 日韩三级视频一区二区三区| 久久久久久九九精品二区国产 | 女人被狂操c到高潮| 久久精品国产亚洲av高清一级| 国产不卡一卡二| 变态另类成人亚洲欧美熟女| 一本精品99久久精品77| 久久九九热精品免费| 亚洲熟妇熟女久久| 天天躁夜夜躁狠狠躁躁| 99久久99久久久精品蜜桃| 亚洲成a人片在线一区二区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产毛片av蜜桃av| 欧美不卡视频在线免费观看 | 久久香蕉激情| e午夜精品久久久久久久| 亚洲中文日韩欧美视频| 亚洲一区二区三区色噜噜| 九色国产91popny在线| 别揉我奶头~嗯~啊~动态视频| 久久久久久大精品| 亚洲欧美日韩无卡精品| 欧美日韩瑟瑟在线播放| 日本在线视频免费播放| 好看av亚洲va欧美ⅴa在| 精品福利观看| 黄色成人免费大全| 精品国产亚洲在线| 一区二区三区高清视频在线| 国产av在哪里看| 最好的美女福利视频网| 亚洲精品一卡2卡三卡4卡5卡| 欧美最黄视频在线播放免费| 国产区一区二久久| 麻豆av在线久日| 亚洲精品国产精品久久久不卡| 91成人精品电影| 变态另类成人亚洲欧美熟女| 免费观看精品视频网站| 十分钟在线观看高清视频www| 啪啪无遮挡十八禁网站| 久久精品影院6| 亚洲狠狠婷婷综合久久图片| 久久九九热精品免费| 精品一区二区三区四区五区乱码| 成年免费大片在线观看| 特大巨黑吊av在线直播 | 亚洲成人精品中文字幕电影| 国产私拍福利视频在线观看| 99热只有精品国产| 搡老岳熟女国产| 久久香蕉精品热| 国产精品亚洲美女久久久| 伦理电影免费视频| 亚洲成a人片在线一区二区| 高清毛片免费观看视频网站| 国产精品精品国产色婷婷| 国产精品久久久久久精品电影 | 免费人成视频x8x8入口观看| 99精品欧美一区二区三区四区| 亚洲成人精品中文字幕电影| 精品一区二区三区视频在线观看免费| 91在线观看av| 中文在线观看免费www的网站 | 欧美zozozo另类| 免费高清视频大片| 亚洲精品久久成人aⅴ小说| 亚洲全国av大片| 欧美精品亚洲一区二区| 热re99久久国产66热| 伦理电影免费视频| 欧美乱码精品一区二区三区| 伊人久久大香线蕉亚洲五| 19禁男女啪啪无遮挡网站| 日韩欧美一区视频在线观看| 日本免费一区二区三区高清不卡| 国产精品国产高清国产av| 成人欧美大片| 国产av一区在线观看免费| 桃色一区二区三区在线观看| 久久性视频一级片| 99国产精品一区二区蜜桃av| 久久香蕉国产精品| 国产高清有码在线观看视频 | 国产精品一区二区免费欧美| 欧美丝袜亚洲另类 | 欧美午夜高清在线| 在线视频色国产色| 国产爱豆传媒在线观看 | 国产激情欧美一区二区| 国产精品久久电影中文字幕| 很黄的视频免费| 好男人在线观看高清免费视频 | 久久国产亚洲av麻豆专区| 一本精品99久久精品77| 精品午夜福利视频在线观看一区| 日韩有码中文字幕| 亚洲精品久久国产高清桃花| 91成年电影在线观看| 深夜精品福利| 亚洲美女黄片视频| 久久久久久九九精品二区国产 | 一级毛片精品| 中文字幕最新亚洲高清| 精品午夜福利视频在线观看一区| 免费在线观看黄色视频的| 十八禁人妻一区二区| 国产高清有码在线观看视频 | 香蕉久久夜色| 成人欧美大片| 国产成人av激情在线播放| 午夜久久久在线观看| 亚洲一区高清亚洲精品| 麻豆成人午夜福利视频| 亚洲精品久久国产高清桃花| 亚洲欧美一区二区三区黑人| 淫秽高清视频在线观看| 波多野结衣高清作品| 一区二区日韩欧美中文字幕| 亚洲成人精品中文字幕电影| 免费看美女性在线毛片视频| 亚洲成a人片在线一区二区| 亚洲午夜精品一区,二区,三区| 国产精品一区二区精品视频观看| 好看av亚洲va欧美ⅴa在| 窝窝影院91人妻| 人人澡人人妻人| 国产又黄又爽又无遮挡在线| 一卡2卡三卡四卡精品乱码亚洲| 狂野欧美激情性xxxx| bbb黄色大片| 可以在线观看的亚洲视频| 男人操女人黄网站| 一本久久中文字幕| 最近最新免费中文字幕在线| 欧美在线一区亚洲| 99国产精品一区二区蜜桃av| 亚洲七黄色美女视频| 亚洲欧美一区二区三区黑人| 一夜夜www| 亚洲七黄色美女视频| 欧美日韩亚洲综合一区二区三区_| 日本精品一区二区三区蜜桃| 91麻豆精品激情在线观看国产| 国产av一区二区精品久久| 日韩欧美一区视频在线观看| 一级片免费观看大全| 黄色丝袜av网址大全| 免费在线观看视频国产中文字幕亚洲| 国产乱人伦免费视频| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久精品电影 | 午夜免费鲁丝| 免费观看精品视频网站| 一区二区三区国产精品乱码| 9191精品国产免费久久| 欧美绝顶高潮抽搐喷水| 国产精品国产高清国产av| 久久国产精品影院| 妹子高潮喷水视频| 欧美一区二区精品小视频在线| 精品电影一区二区在线| 99国产极品粉嫩在线观看| 亚洲成人精品中文字幕电影| 久久国产精品男人的天堂亚洲| avwww免费| 国产一区二区三区在线臀色熟女| 欧美成人免费av一区二区三区| 又黄又爽又免费观看的视频| 亚洲精品在线观看二区| 黄色视频不卡| 99在线视频只有这里精品首页| 级片在线观看| 欧美三级亚洲精品| 免费av毛片视频| 亚洲精品美女久久久久99蜜臀| 午夜福利免费观看在线| 日本a在线网址| 嫁个100分男人电影在线观看| 1024手机看黄色片| 免费在线观看成人毛片| 首页视频小说图片口味搜索| svipshipincom国产片| 国内精品久久久久精免费| 一区二区三区国产精品乱码| 不卡av一区二区三区| 一级a爱视频在线免费观看| 亚洲三区欧美一区| 麻豆久久精品国产亚洲av| 国产精品综合久久久久久久免费| 欧美日本亚洲视频在线播放| 一区二区日韩欧美中文字幕| 欧美黑人欧美精品刺激| 国产亚洲欧美在线一区二区| 免费看美女性在线毛片视频| 国产精品 国内视频| 久久香蕉激情| 精品欧美一区二区三区在线| 欧美久久黑人一区二区| 脱女人内裤的视频| 久久久久免费精品人妻一区二区 | 999久久久精品免费观看国产| 欧美又色又爽又黄视频| avwww免费| 在线天堂中文资源库| 午夜亚洲福利在线播放| 1024视频免费在线观看| 黄色片一级片一级黄色片| 欧美一级毛片孕妇| 两个人看的免费小视频| 在线国产一区二区在线| 日本熟妇午夜| 免费在线观看影片大全网站| 亚洲国产精品999在线| 日韩欧美三级三区| 欧美日韩亚洲综合一区二区三区_| 国产精品久久视频播放| 欧美日韩一级在线毛片| 精品日产1卡2卡| 日日爽夜夜爽网站| 又黄又爽又免费观看的视频| 久久婷婷成人综合色麻豆| 波多野结衣av一区二区av| 国产一卡二卡三卡精品| 亚洲性夜色夜夜综合| 视频在线观看一区二区三区| 国产麻豆成人av免费视频| 久久久久国产一级毛片高清牌| 久久人妻av系列| 超碰成人久久| 日韩欧美 国产精品| 欧美av亚洲av综合av国产av| 19禁男女啪啪无遮挡网站| 亚洲欧洲精品一区二区精品久久久| 欧美一区二区精品小视频在线| 国产精品久久视频播放| 国内精品久久久久久久电影| 在线播放国产精品三级| 亚洲中文日韩欧美视频| 久久香蕉精品热| 91成人精品电影| АⅤ资源中文在线天堂| cao死你这个sao货| 很黄的视频免费| 国产aⅴ精品一区二区三区波| 中文资源天堂在线| 欧美成人午夜精品| 波多野结衣av一区二区av| 99热只有精品国产| 91麻豆av在线| 亚洲熟女毛片儿| 一夜夜www| 免费看日本二区| 亚洲av中文字字幕乱码综合 | 黄色视频不卡| 俺也久久电影网| 91麻豆av在线| 99久久99久久久精品蜜桃| 国产午夜福利久久久久久| videosex国产| 90打野战视频偷拍视频| 欧美乱妇无乱码| 精品久久久久久久久久久久久 | 欧美不卡视频在线免费观看 | 老司机深夜福利视频在线观看| 午夜福利高清视频| 岛国视频午夜一区免费看| 日本 欧美在线| 99久久综合精品五月天人人| 亚洲午夜理论影院| 亚洲精品在线美女| 国产精品久久久人人做人人爽| 精品免费久久久久久久清纯| 精品一区二区三区四区五区乱码| 欧美最黄视频在线播放免费| 亚洲三区欧美一区| 国产高清激情床上av| 黄频高清免费视频| xxxwww97欧美| 色综合婷婷激情| 少妇粗大呻吟视频| 亚洲国产欧洲综合997久久, | 精品国产国语对白av| 老熟妇仑乱视频hdxx| 亚洲精品色激情综合| 久久久久国内视频| svipshipincom国产片| 老熟妇乱子伦视频在线观看| a级毛片a级免费在线| 免费一级毛片在线播放高清视频| 男女做爰动态图高潮gif福利片| 在线观看日韩欧美| 天天一区二区日本电影三级| 国产精品二区激情视频| 欧美丝袜亚洲另类 | 免费在线观看日本一区| 久久婷婷成人综合色麻豆| 淫秽高清视频在线观看| 日日爽夜夜爽网站| 国产av不卡久久| 国内精品久久久久精免费| 男人舔奶头视频| 久久中文字幕人妻熟女| 亚洲免费av在线视频| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品一区二区www| av天堂在线播放| 免费观看人在逋| 可以免费在线观看a视频的电影网站| 搞女人的毛片| 男人舔女人下体高潮全视频| 亚洲五月天丁香| 欧美国产精品va在线观看不卡| 精品久久久久久久毛片微露脸| 国产激情偷乱视频一区二区| 精品国产美女av久久久久小说| 三级毛片av免费| 侵犯人妻中文字幕一二三四区| 嫩草影院精品99| 中文在线观看免费www的网站 | 嫩草影视91久久| 亚洲色图 男人天堂 中文字幕| 亚洲色图av天堂| 亚洲av美国av| 在线观看免费午夜福利视频| 91在线观看av| 精品国产一区二区三区四区第35| 色婷婷久久久亚洲欧美| 搡老岳熟女国产| 人成视频在线观看免费观看| 国产片内射在线| 黄色 视频免费看| 视频区欧美日本亚洲| 在线av久久热| 国产精品98久久久久久宅男小说| 18禁国产床啪视频网站| 好男人在线观看高清免费视频 | 国产三级在线视频| 国产高清有码在线观看视频 | 变态另类丝袜制服| xxxwww97欧美| 欧美不卡视频在线免费观看 | 国产av又大| 国产一级毛片七仙女欲春2 | 亚洲精品在线美女| 视频区欧美日本亚洲| 久久久久国产一级毛片高清牌| 一级a爱视频在线免费观看| 亚洲国产欧洲综合997久久, | 女同久久另类99精品国产91| 国产av不卡久久| 在线av久久热| 国产精品亚洲一级av第二区| 亚洲精品国产精品久久久不卡| 在线av久久热| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 久久久久久亚洲精品国产蜜桃av| 777久久人妻少妇嫩草av网站| 国产精品一区二区三区四区久久 | 丰满人妻熟妇乱又伦精品不卡| 黄色 视频免费看| 午夜亚洲福利在线播放| 九色国产91popny在线| 人人妻,人人澡人人爽秒播| 夜夜爽天天搞| 日日爽夜夜爽网站| 久久狼人影院| 99精品久久久久人妻精品| 欧美性长视频在线观看| 国产91精品成人一区二区三区| 中文字幕精品亚洲无线码一区 | 日韩免费av在线播放| 久久久久久久精品吃奶| 91老司机精品| 在线视频色国产色| 香蕉国产在线看| 国产精品乱码一区二三区的特点| 88av欧美| 婷婷丁香在线五月| 国产高清视频在线播放一区| 满18在线观看网站| 可以在线观看毛片的网站| 搡老岳熟女国产| 久久久国产成人免费| 成人国产综合亚洲| 国产精品九九99| svipshipincom国产片| 日韩精品青青久久久久久| 夜夜爽天天搞| 亚洲午夜理论影院| 久久久久久久久久黄片| 中出人妻视频一区二区| 老熟妇仑乱视频hdxx| 国产成人精品久久二区二区91| 高潮久久久久久久久久久不卡| 亚洲第一av免费看| 国产麻豆成人av免费视频| 日本成人三级电影网站| 制服丝袜大香蕉在线| 成人一区二区视频在线观看| 婷婷六月久久综合丁香| e午夜精品久久久久久久| 欧美日韩黄片免| 欧美在线一区亚洲| 国产亚洲精品av在线| 亚洲男人天堂网一区| 婷婷丁香在线五月| 亚洲av成人不卡在线观看播放网| 欧美成人一区二区免费高清观看 | 日韩 欧美 亚洲 中文字幕| 一进一出抽搐动态| 亚洲中文字幕日韩| 一个人观看的视频www高清免费观看 | 国产精品爽爽va在线观看网站 | 中文资源天堂在线| 久久久水蜜桃国产精品网| 国产精品 欧美亚洲| 欧美av亚洲av综合av国产av| 日韩有码中文字幕| 国产99白浆流出| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| a级毛片a级免费在线| 一级毛片女人18水好多| 91麻豆精品激情在线观看国产| 精品午夜福利视频在线观看一区| 国产精品影院久久| 在线观看免费午夜福利视频| 精品国产亚洲在线| 午夜福利在线观看吧| 亚洲国产精品合色在线| 欧美日韩亚洲综合一区二区三区_| 亚洲精品av麻豆狂野| 男女床上黄色一级片免费看| 国产av一区在线观看免费| 男女下面进入的视频免费午夜 | 成人18禁高潮啪啪吃奶动态图| 午夜福利欧美成人| 国产v大片淫在线免费观看| 国产男靠女视频免费网站|