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

    西藏尼雄鐵礦成礦花崗巖成因及其對成礦構(gòu)造背景的啟示

    2015-01-19 03:40:59范淑芳曲曉明辛洪波
    大地構(gòu)造與成礦學(xué) 2015年2期
    關(guān)鍵詞:洋殼班公湖雅魯藏布江

    范淑芳, 曲曉明, 宋 揚(yáng), 辛洪波

    (中國地質(zhì)科學(xué)院 礦產(chǎn)資源研究所, 北京 100037)

    西藏尼雄鐵礦成礦花崗巖成因及其對成礦構(gòu)造背景的啟示

    范淑芳, 曲曉明, 宋 揚(yáng), 辛洪波

    (中國地質(zhì)科學(xué)院 礦產(chǎn)資源研究所, 北京 100037)

    本文分析了西藏尼雄鐵礦與成礦有關(guān)花崗閃長巖的主量元素、微量元素、稀土元素及Sr-Nd-Pb同位素特征, 并做了鋯石LA-ICPMS U-Pb年齡測定。巖石地球化學(xué)分析結(jié)果顯示樣品為亞鋁質(zhì)中鉀–高鉀鈣堿性巖系列, 屬I-型花崗巖。地球化學(xué)組成上其富集大離子親石元素(LILE)Rb、Ba、Sr、Th、U、K、Pb, 虧損高場強(qiáng)元素(HFSE)Nb、Ta、Ti, 具有典型的島弧巖漿作用的特征; 稀土元素球粒隕石標(biāo)準(zhǔn)化圖解表現(xiàn)為富集輕稀土的右傾型式(LREE/HREE=5.67~8.37), 無Eu異常, 顯示活動大陸邊緣巖漿巖的稀土配分特征。巖體ISr為0.707625~0.710997, εNd(t)為–6.6~–8.7,206Pb/204Pb,207Pb/204Pb,208Pb/204Pb比值分別為18.786~18.955、15.694~15.726、39.355~39.676, 顯示出富集地幔特征(EMⅡ), 表明地殼組分對巖漿生成有重要影響。所測巖體的鋯石206Pb/238U 加權(quán)平均年齡為112.09±0.54 Ma(MSWD=0.45), 表明尼雄花崗巖體形成于早白堊世晚期。綜合分析班公湖–怒江中特提斯洋和雅魯藏布江新特提斯洋的演化歷史, 作者認(rèn)為尼雄鐵礦是早白堊世雅魯藏布江洋殼板塊向北側(cè)拉薩地塊之下俯沖的構(gòu)造背景下的產(chǎn)物。而成礦巖體主要是俯沖組分(流體和熔融)對地幔楔交代改造的結(jié)果, 是地幔巖漿底侵引起下地殼物質(zhì)部分熔融, 兩種巖漿混合而形成的。

    巖石地球化學(xué); Sr-Nd-Pb同位素; 鋯石U-Pb年齡; 尼雄鐵礦; 西藏

    地處青藏高原南部的岡底斯帶(也稱為拉薩地塊)夾持于雅魯藏布江縫合帶和班公湖–怒江縫合帶之間, 區(qū)內(nèi)構(gòu)造巖漿活動強(qiáng)烈, 成礦地質(zhì)條件優(yōu)越,孕育了大量的銅、金、銀、鉬、富鐵、鉛鋅等礦產(chǎn)資源, 是青藏高原上最重要的成礦區(qū)帶。目前對岡底斯成礦帶的研究, 多集中于新生代斑巖型銅金成礦作用(曲曉明等, 2001; 侯增謙等, 2003)和矽卡巖型鉛鋅銀成礦作用(佘宏全等, 2005; 李光明等, 2011), 對鐵礦的研究明顯偏少。尼雄鐵礦是近年來江西省地質(zhì)調(diào)查研究院在1:25萬填圖過程中發(fā)現(xiàn)的一個超大型矽卡巖型富鐵礦床(曹圣華等, 2007;袁健芽等, 2008; 王樂等2012)。該礦床南距雅魯臧布江縫合帶約150 km, 北距班公湖–怒江縫合帶約130 km, 基本處于岡底斯帶的中帶位置。目前該礦床研究中爭議最大的問題是它的產(chǎn)出構(gòu)造背景, 即它是由雅魯藏布江洋殼向北俯沖形成的?還是由班公湖–怒江洋殼向南俯沖形成的?亦或是拉薩地塊板內(nèi)巖漿作用的產(chǎn)物?從目前的研究結(jié)果看, 多數(shù)人傾向于該鐵礦與班公湖–怒江洋殼向南俯沖有關(guān)(袁健芽等, 2008; 張曉倩等, 2010; 曹圣華, 2012;于玉帥等, 2012), 也有人認(rèn)為是雅魯藏布江洋殼向北俯沖的結(jié)果(吳旭玲和陳振華, 2005)。然而, 對該礦床形成構(gòu)造背景的正確認(rèn)識不能僅依靠礦床本身,還涉及到對班公湖–怒江中特提斯洋演化歷史的厘定。本文根據(jù)曲曉明等(2012)獲得的有關(guān)班公湖–怒江洋盆閉合的證據(jù), 通過對尼雄鐵礦成礦巖體的巖石地球化學(xué)分析和LA-ICPMS鋯石U-Pb定年, 闡述了尼雄鐵礦與班公湖–怒江和雅魯藏布江兩條縫合帶之間的關(guān)系, 為該礦床的形成環(huán)境提供了新的依據(jù)。并利用成礦巖體的Sr、Nd、Pb同位素特征, 進(jìn)一步探討了該礦床的巖漿源區(qū)組成。通過本文的研究, 揭示岡底斯帶鐵礦的成礦地質(zhì)條件和成礦規(guī)律,更好地為區(qū)域地質(zhì)找礦服務(wù)。

    1 地質(zhì)背景及巖石學(xué)特征

    1.1 成礦地質(zhì)條件

    岡底斯帶夾持于南側(cè)的印度河–雅魯藏布江縫合帶與北部的班公湖–怒江縫合帶之間, 是一條東西長約2500 km, 南北寬150~300 km, 面積達(dá)45萬平方公里的巨型構(gòu)造–巖漿巖帶(潘桂棠等, 2006)。中生代以來印度大陸向北側(cè)歐亞大陸的俯沖以及新生代的陸–陸碰撞孕育了豐富的銅金多金屬礦產(chǎn)資源。同時,由于其處于兩個大陸碰撞造山的關(guān)鍵部位, 岡底斯帶一直被看作是窺探青藏高原隆升機(jī)制的窗口和研究全球大陸動力學(xué)的天然實(shí)驗(yàn)室。

    尼雄鐵銅礦田處于隆格爾–工布江達(dá)斷裂帶南側(cè)、岡底斯陸緣火山–巖漿弧以北、措勤–納木錯初始弧間盆地南部(吳旭玲和陳振華, 2005)。礦田由北西至南東方向劃分為5個礦區(qū): 滾糾鐵礦區(qū)、尼雄鐵礦區(qū)、沙松南鐵礦區(qū)、毛加崍鐵礦區(qū)及日阿銅礦區(qū)(圖1)。尼雄鐵礦區(qū)內(nèi)出露的主要地層為石炭系永珠組, 二疊系拉嘎組、昂杰組、下拉組、敵布錯組及古近系日貢拉組。中二疊統(tǒng)下拉組和上二疊統(tǒng)敵布錯組是區(qū)內(nèi)磁鐵礦床形成的重要條件。下拉組為灰–淺灰色厚層狀(含)生物碎屑微晶灰?guī)r、中厚層狀含菱鐵礦灰?guī)r、含鮞狀赤鐵礦灰?guī)r; 敵布錯組為深灰–灰白色中厚層狀含礫長石砂巖、含鮞狀赤鐵礦長石石英砂巖、含炭巖屑砂巖與黑色薄層狀含炭泥質(zhì)粉砂巖、含炭粉砂質(zhì)泥巖互層, 下部夾深灰色薄層狀硅質(zhì)巖, 上部夾深灰色微晶灰?guī)r。

    圖1 西藏尼雄礦田地質(zhì)簡圖(據(jù)曹圣華等, 2010修改)Fig.1 Sketch geological map of the Nixiong ore field

    鐵礦體主要產(chǎn)于中細(xì)粒黑云母花崗閃長巖、中細(xì)粒黑云母二長花崗巖與二疊紀(jì)下拉組、敵布錯組地層的外接觸帶、下拉組和敵布錯組的層間破碎帶及二者的接觸面上。斷裂構(gòu)造主要表現(xiàn)為NWW向的推覆構(gòu)造帶。礦區(qū)燕山晚期巖漿活動強(qiáng)烈, 一系列中酸性侵入體沿NWW向呈串珠狀分布, 巖性主要為中細(xì)粒角閃黑云花崗閃長巖、中細(xì)粒二長花崗巖、細(xì)粒斑狀花崗巖和花崗斑巖脈。鐵礦體呈 NWW 向產(chǎn)于矽卡巖帶內(nèi)及敵布錯和下拉組地層中(圖1)。

    1.2 成礦巖體地質(zhì)及巖相學(xué)

    尼雄鐵銅礦田早白堊世晚期侵入的成礦巖漿巖出露面積約180 km2, 呈多個不規(guī)則橢圓形、圓形巖株產(chǎn)出, 沿NWW向呈串珠狀分布。本文研究的成礦巖體樣品采自尼雄鐵礦區(qū)(采樣位置見圖1), 為中細(xì)粒花崗閃長巖, 具中細(xì)?;◢徑Y(jié)構(gòu), 塊狀構(gòu)造。巖株存在巖相分帶, 內(nèi)部巖性為二長花崗巖。巖石主要礦物組成為斜長石(45%~50%)、鉀長石(15%~20%)、石英(20%~25%)、黑云母(<5%)、角閃石(5%), 副礦物為磷灰石、鋯石、磁鐵礦、榍石等(<1%)。巖體與二疊紀(jì)地層接觸帶附近蝕變強(qiáng)烈, 主要有矽卡巖化、大理巖化、黃鐵絹英巖化、綠泥石化、綠簾石化、磁鐵礦化等。

    2 分析測試方法

    筆者在三個采樣點(diǎn)共采集12件新鮮樣品做了全巖主量元素、微量元素和稀土元素含量分析(分析結(jié)果見表1)。樣品分析在北京核工業(yè)地質(zhì)分析研究中心完成。主量元素采用PHILLIPSPW-2404型X-熒光光譜儀分析完成, 精度好于1%; 微量元素和稀土元素采用ICP-MS測定, 儀器型號為ELEMENT-2質(zhì)譜儀, 分析精度優(yōu)于2%, 分析流程見Qu et al. (2004)。

    尼雄鐵礦成礦花崗巖12件樣品的Pb、Sr、Nd同位素分析在北京核工業(yè)地質(zhì)分析測試研究中心完成, 所用儀器是MAT-261固體同位素質(zhì)譜儀。Nd和Sr同位素分析, 樣品(200目粉末)先用HF+HNO3完全溶解, 蒸干后再用HCl溶解。樣品中的Nd、Sr采用常規(guī)的離子交換技術(shù)分離。Nd同位素測定用146Nd/144Nd=0.7219進(jìn)行質(zhì)量分餾校正。標(biāo)樣J.M.Nd2O3測定結(jié)果143Nd/144Nd=0.511125±8(2σ), GBW04419測定結(jié)果143Nd/144Nd=0.512725±10(2σ), Sm、Nd流程空白為5×10–11g。Sr同位素測定用88Sr/86Sr= 8.37521進(jìn)行質(zhì)量分餾校正。標(biāo)樣NBS987 SrCO3測定結(jié)果為87Sr/86Sr=0.71025±(2σ); Rb、Sr流程空白為10–9~10–11g。Pb同位素分析, 樣品用三酸(HF+H2SO4+ HCl)分解后用離子交換樹脂分離出鉛, 蒸干后進(jìn)行同位素測試。分析精度對1 μg鉛含量204Pb/206Pb低于0.05%,208Pb/206Pb一般不大于0.005%。對國際標(biāo)樣NBS981的測試結(jié)果為:208Pb/206Pb=2.162189, 誤差為0.0027%;207Pb/206Pb=0.913626, 誤差為0.0059%;204Pb/206Pb=0.059201, 誤差為0.0015%所用鋯石從樣品NX11-04中分選出, 樣品經(jīng)人工破碎后用重力和磁選的方法分選出鋯石, 之后在雙目鏡下去除雜質(zhì),使鋯石單礦物純度達(dá)到95%以上。將待測鋯石及標(biāo)準(zhǔn)鋯石用環(huán)氧樹脂固定于載玻片上, 磨至鋯石顆粒的2/3使其內(nèi)部結(jié)構(gòu)充分暴露, 用透射光、反射光和陰極發(fā)光(CL)照相, 以選取合適部位進(jìn)行U-Pb同位素定年。鋯石CL照相在國土資源部同位素地質(zhì)重點(diǎn)實(shí)驗(yàn)室完成, LA-ICPMS鋯石U-Pb年齡測定在國土資源部成礦作用與資源評價重點(diǎn)實(shí)驗(yàn)室完成, 采用的儀器分別為Agilent 7500型ICPMS和ComPex 102 ArF準(zhǔn)分子激光器, 工作物質(zhì)ArF波長193 nm,與GeoLas 200M光學(xué)系統(tǒng)聯(lián)機(jī)進(jìn)行。鋯石年齡采用國際標(biāo)準(zhǔn)鋯石91500作外標(biāo), 元素含量采用NIST SRM 610作外標(biāo),29Si作內(nèi)標(biāo), 數(shù)據(jù)用Ludwig SQUID 1.0程序處理, 應(yīng)用Isoplot程序完成年齡計(jì)算以及諧和圖的繪制。

    3 分析結(jié)果

    3.1 主量元素

    從表1分析結(jié)果可知, 尼雄礦區(qū)成礦巖體的SiO2含量在63.63%~67.37%之間, 平均為65.72%; K2O含量為2.19%~3.07%, Na2O含量為3.09%~4%,全堿含量5.76%~6.81%, Na2O/K2O>1, 樣品在TAS圖(圖2a)上落入花崗閃長巖區(qū), 在K2O-SiO2圖(圖2b)上顯示樣品為中鉀–高鉀鈣堿性巖。巖體的CaO和MgO含量分別為3.73%~4.72%和1.51%~2.11%, FeOT含量2.28%~5.27%, Al2O3含量15.37%~16.24%,平均為15.71%, 在A/NK-A/CNK圖(圖2c)上位于亞鋁質(zhì)巖區(qū)。在R2-R1圖(圖2d)上它們都落在碰撞前階段, 表現(xiàn)為板塊俯沖階段巖漿作用產(chǎn)物。在P2O5-SiO2圖和Na2O-K2O圖(圖2e, 2f)上樣品都顯示出I-型花崗巖特征, 表明其源巖為火成巖。

    表1 尼雄花崗閃長巖主量元素(%)、微量元素和稀土元素(μg/g)分析結(jié)果Table 1 Major (%) and trace (μg/g) element results of the granodiorites in the Nixiong iron deposit

    圖2 尼雄花崗閃長巖地球化學(xué)圖: TAS圖(a, 據(jù)Middlemost, 1994); K2O-SiO2圖(b, 據(jù)Peccerillo and Taylor, 1976); A/NK-A/CNK圖(c, 據(jù)Maniar and Piccoli, 1989); R2-R1圖(d, 據(jù)Batchelor and Bowden, 1985); P2O5-SiO2圖(e, 據(jù)李獻(xiàn)華等, 2007); Na2O-K2O圖(f, 據(jù)Collins et al., 1982)Fig.2 Geochemical plots for the granodiorites in the Nixiong iron deposit: TAS (a); K2O vs SiO2(b); A/NK vs A/CNK (c); R2vs R1(d); P2O5vs SiO2(e); Na2O vs K2O (f)

    3.2 微量元素

    在原始地幔標(biāo)準(zhǔn)化微量元素蛛網(wǎng)圖中(圖3a),尼雄花崗閃長巖表現(xiàn)為富集大離子親石元素Rb、Ba、Th、U、K、Pb、Sr, 虧損高場強(qiáng)元素Nb、Ta、 Ti, 顯示出島弧巖漿巖的基本特征(Wilson, 1989)。大離子不相容元素 Sr 呈弱富集, Ba 則相對虧損, 說明巖漿經(jīng)歷過分異演化。在構(gòu)造環(huán)境判別圖上(圖4),巖體均落在島弧型巖漿巖區(qū)。

    圖3 尼雄花崗閃長巖的微量元素原始地幔標(biāo)準(zhǔn)化蛛網(wǎng)圖(a)和稀土元素球粒隕石標(biāo)準(zhǔn)化分布曲線圖(b) (原始地幔和球粒隕石標(biāo)準(zhǔn)化值引用 Sun and McDonough, 1989)Fig.3 Primitive mantle normalized trace element spider diagram (a) and chondrite-normalized REE patterns (b) of the granodiorites in the Nixiong iron deposit

    圖4 尼雄花崗閃長巖Nb-Y (a)、Ta-Yb (b)、Rb-(Y+Nb) (c)和Rb-(Yb+Ta) (d)圖(據(jù)Pearce et al., 1984)Fig.4 Nb vs Y (a), Ta vs Yb (b), Rb vs (Y+Nb) (c), and Rb vs (Yb+Ta) (d) diagrams for the granodiorites in the Nixiong iron deposit

    3.3 稀土元素

    樣品稀土元素含量較高, ΣREE為84.79~159.92 μg/g,稀土元素球粒隕石標(biāo)準(zhǔn)化分布曲線表現(xiàn)為富集輕稀土元素的右傾型(圖3b), LREE/HREE=5.67~8.37, (La/Yb)N=4.8~9.06。巖體的中、重稀土元素沒有明顯分異, 符合活動大陸邊緣巖漿巖的稀土配分型式。巖體δEu接近于1, 無Eu異常, 說明含礦巖漿形成演化過程中沒有發(fā)生明顯的斜長石分離結(jié)晶作用或在巖漿源區(qū)的無明顯殘留。

    3.4 鋯石U-Pb年齡

    尼雄鐵礦成礦花崗閃長巖的鋯石顆粒自形程度較好, 呈長柱狀或短柱狀, 大小在100~250 μm之間,發(fā)育震蕩環(huán)帶(圖5a), 鋯石U、Th含量分別為55~ 653 μg/g和52~1107 μg/g, Th/U比值變化于0.94~ 1.72之間(表2), 遠(yuǎn)大于0.4, 說明其為巖漿成因鋯石(吳元保和鄭永飛, 2004), 而這些鋯石的結(jié)晶年齡可以代表巖體的侵位時間。17顆鋯石的206Pb/238U加權(quán)平均年齡為112.09±0.54 Ma, MSWD= 0.45(圖5b), 說明尼雄花崗閃長巖體形成于早白堊世晚期。

    圖5 尼雄花崗閃長巖(NX11-04)中的鋯石CL圖像(a)和U-Pb年齡諧和圖(b)Fig.5 CL images (a) and U-Pb concordant diagram (b) for zircon from the granodiorites (NX11-04) in the Nixiong iron deposit

    表2 尼雄花崗閃長巖鋯石U-Pb同位素分析結(jié)果Fig.2 U-Pb isotope results of zircon from the granodiorites in the Nixiong iron deposit

    3.5 Sr、Nd同位素

    尼雄花崗閃長巖87Sr/86Sr比值為0.709539~ 0.712390, ISr比值為0.707625~0.710997(表3), 較高的ISr值指示成礦巖漿中有地殼物質(zhì)的參與。DePaolo (1981)研究指出, 地球上不同構(gòu)造位置的近代火山巖的εNd是不同的, 源自地幔的火山巖其εNd>0, 若εNd<0, 則說明這些巖石形成于地殼中或受到了地殼物質(zhì)的強(qiáng)烈混染。尼雄花崗閃長巖143Nd/144Nd比值為0.512135~0.512237, εNd(t)為–6.6~–8.7(表3), 顯示出地殼組分對巖漿生成的重要影響。在εNd(t)-ISr圖上(圖6a), 尼雄鐵礦成礦花崗閃長巖處于地幔玄武巖與下地殼的混合線附近, 指示成礦巖漿可能是由殼幔巖漿混合而成。

    3.6 Pb同位素

    尼雄礦區(qū)成礦花崗閃長巖Pb同位素組成變化不大,206Pb/204Pb、207Pb/204Pb、208Pb/204Pb比值分別為18.786~18.955, 15.694~15.726, 39.355~39.676(表3)。Pb同位素比值206Pb/204Pb>18、207Pb/204Pb>15.5、208Pb/204Pb>38, 反映了高放射成因鉛特征(陳能松等, 2007)。在Zartman and Doe (1981)的鉛構(gòu)造模式圖上(圖7a、b) 尼雄鐵礦區(qū)成礦花崗閃長巖樣品均處于EMⅡ單元內(nèi), 說明尼雄鐵礦的成礦巖漿源區(qū)具EMⅡ型富集地幔特征。而EMⅡ地幔端元的形成與俯沖和再循環(huán)的地殼物質(zhì)密切相關(guān)(Hart, 1988), 這些由同位素得出的巖漿源區(qū)信息與巖石地球化學(xué)所反映的島弧環(huán)境是完全一致的。在Sr-Pb、Nd-Pb同位素圖上(圖8a、b)也同樣反映出俯沖組分對巖漿源區(qū)的影響。

    4 討 論

    4.1 成礦構(gòu)造背景

    尼雄鐵礦床產(chǎn)于岡底斯帶內(nèi)部, 南距雅魯藏布江縫合帶約150 km, 北距班公湖–怒江縫合帶約130 km。作為矽卡巖型鐵礦, 尼雄鐵礦的成因已沒有太多異義, 而目前爭論最大的是關(guān)于該礦床形成的構(gòu)造環(huán)境, 即尼雄鐵礦的形成是與雅魯藏布江洋殼的向北俯沖有關(guān)?還是與班公湖–怒江洋殼的向南俯沖有關(guān)?抑或是拉薩地塊內(nèi)部板內(nèi)巖漿作用的產(chǎn)物?吳旭玲和陳振華(2005)通過對尼雄巖體巖石地球化學(xué)的研究并結(jié)合其所處的大地構(gòu)造位置, 將尼雄鐵礦的形成歸因于雅魯藏布江大洋板塊向北俯沖消減;袁建芽等(2008)通過對岡底斯南帶的斑巖型礦床與本區(qū)的矽卡巖型鐵銅礦床的對比研究, 認(rèn)為該礦田的形成受班公湖–怒江中特提斯洋殼向南俯沖控制;張曉倩等(2010)通過區(qū)域宏觀構(gòu)造演化研究, 推斷尼雄巖體很可能侵位于一種與班公湖–怒江洋殼巖石圈板片俯沖角度的變陡進(jìn)而斷離有關(guān)的伸展背景。曹圣華(2012)以及于玉帥等(2012)通過對滾糾鐵礦床地球化學(xué)特征的研究, 認(rèn)為早白堊世晚期羌塘陸殼與岡底斯島弧陸殼對接碰撞, 使兩者會聚速率減緩, 班公湖–怒江中特提斯洋殼繼續(xù)受重力或/和俯沖慣性拖拽作用下向南俯沖, 導(dǎo)致俯沖的洋殼板片回轉(zhuǎn), 誘使軟流圈地幔上涌, 引發(fā)巖石圈地幔和上覆地殼發(fā)生部分熔融, 進(jìn)而形成了以殼源為主的殼?;旌铣傻V巖漿。

    表3 尼雄花崗閃長巖Pb、Sr、Nd同位素組成Table 3 Pb, Sr and Nd isotopic compositions of the granodiorites in the Nixiong iron deposit

    圖6 尼雄花崗閃長巖 εNd(t)-ISr(a, 據(jù)Depaolo and Wassethurg, 1979)和 (La/Yb)N-δEu圖(b, 據(jù)張宇等, 2013)Fig.6 εNd(t) vs ISr(a) and (La/Yb)Nvs δEu (b) diagrams for the granodiorites in the Nixiong iron deposit

    圖7 尼雄花崗閃長巖207Pb/204Pb-206Pb/204Pb (a)和208Pb/204Pb-206Pb/204Pb (b)圖(底圖據(jù)Zartman and Doe, 1981)Fig.7207Pb/204Pb vs206Pb/204Pb (a) and208Pb/204Pb vs206Pb/204Pb (b) diagrams for the granodiorites in the Nixiong iron deposit

    圖8 Sr-Pb (a)和 Nd-Pb (b)同位素圖解(據(jù)Zindler and Hart, 1986)Fig.8 Sr vs Pb (a) and Nd vs Pb (b) isotopic diagrams for the granodiorites in the Nixiong iron deposit

    本文研究表明, 尼雄鐵礦與成礦有關(guān)的花崗閃長巖富集大離子不相容元素Rb、Ba、Th、U、K、Sr、Pb, 虧損高場強(qiáng)元素Nb、Ta、Ti, 具有島弧型花崗巖的鮮明特征。R1-R2圖(圖2d)反映了造山帶演化所經(jīng)歷的一個完整的Wilson旋回上, 尼雄鐵礦的成礦花崗閃長巖在該圖上處于碰撞前階段, 指示其為板塊俯沖階段巖漿作用的產(chǎn)物。在構(gòu)造環(huán)境判別圖上(圖4), 這些花崗閃長巖也表現(xiàn)出火山弧的形成環(huán)境。在微量元素比值Ce/Pb-Ce和Nb/Th-Nb圖上(圖9a, b), 它們也顯示出島弧火山巖的特性。

    圖9 尼雄花閃長崗巖Ce/Pb-Ce (a)和Nb/Th-Nb (b)圖(底圖據(jù) Boztug et al. 2007; 原始地幔數(shù)據(jù)據(jù) Hofmann, 1988; 大陸地殼、MORB+OIB和島弧型火山巖區(qū)域據(jù)Schmidberger and Henger, 1999)Fig.9 Ce/Pb vs Ce (a) and Nb/Th vs Nb (b) diagrams for the granodiorites in the Nixiong iron deposit

    在尼雄鐵礦本身的形成構(gòu)造環(huán)境已經(jīng)明確的情況下, 要查明該鐵礦的成礦背景, 即是與南面的雅魯藏布江洋殼俯沖有關(guān)還是與北面的班公湖–怒江洋殼俯沖有關(guān), 還需要先確定兩條縫合帶的演化歷史。作為拉薩地塊北緣的班公湖–怒江縫合帶構(gòu)造演化歷史十分復(fù)雜, 人們對該縫合帶的地質(zhì)問題都還存在較大的爭論, 包括班公湖–怒江中特提斯洋盆的打開時間、閉合時間、俯沖方向、俯沖極性乃至洋盆屬性(即是一個統(tǒng)一的大洋還是多個局限性洋盆)。對尼雄鐵礦的成礦來講, 重要的是要查明班公湖–怒江中特提斯洋殼向南側(cè)拉薩地塊之下俯沖的時間或者說該洋盆的閉合時間。曹圣華等(2006)在綜合分析前人成果的基礎(chǔ)上, 根據(jù)班公湖–怒江縫合帶南、北兩側(cè)大陸邊緣古特提斯地層沉積特征的差異, 推斷中侏羅世–晚侏羅世早期為中特提斯洋雙向俯沖階段, 班公湖–怒江洋向南俯沖消減于拉薩地塊之下。Kapp et al. (2003)根據(jù)沉積地層的覆蓋關(guān)系, 推斷班公湖–怒江洋盆的閉合時間為侏羅紀(jì)末–白堊紀(jì)初, 即145 Ma左右。Zhu et al.(2009, 2011)通過對拉薩地塊北緣巖漿巖的研究,認(rèn)為 110 Ma 左右班公湖–怒江特提斯洋殼可能發(fā)生過板片斷離, 此時本地區(qū)已經(jīng)進(jìn)入板內(nèi)環(huán)境。曲曉明等(2012, 2013)通過對班公湖–怒江縫合帶中段新發(fā)現(xiàn)的A型花崗巖的研究, 也認(rèn)為這些110Ma±形成的A-型花崗巖是俯沖的洋殼板片斷離, 軟流圈地幔沿板片窗上涌的結(jié)果。并進(jìn)一步推斷班公湖–怒江特中提斯洋盆的閉合時間應(yīng)為早白堊世初(140~130 Ma之間)。Sui et al.(2013)通過拉薩地體北部不同類型火山巖巖石地球化學(xué)的研究, 認(rèn)為140~110 Ma該地區(qū)已處于同碰撞階段。由以上討論可知早白堊世晚期班公湖–怒江縫合帶已處于碰撞后階段, 與尼雄鐵礦形成時所處的巖漿弧的構(gòu)造環(huán)境并不相符。

    關(guān)于南面雅魯藏布江新特提斯樣的演化歷史人們的認(rèn)識還比較一致。Garzanti et al. (1999) 通過對吉隆縣南部吉隆溝二疊紀(jì)玄武巖的研究, 提出這些早二疊世火山巖是岡瓦納大陸北緣大陸裂谷作用的產(chǎn)物, 代表了雅魯藏布江新特提斯的初始開啟。關(guān)于雅魯藏布江新特提斯洋的俯沖時間, 早期的研究者提出晚白堊世到始新世為特提斯洋殼的俯沖階段(髙延林, 1985, 1988), 而近期的研究大多認(rèn)為特提斯洋殼的早期俯沖開始于晚侏羅世到早白堊世(劉志飛等, 2000;朱弟成等, 2009; 徐向珍等, 2011; 張萬平等, 2011; McDermid et al., 2002)。對于雅魯藏布江洋盆的閉合時間, 郝杰等(1999)通過對日喀則弧前盆地海相沉積巖的研究, 提出特提斯–喜瑪拉雅洋消亡和印度板塊與歐亞大陸碰撞應(yīng)發(fā)生在始新世末期或漸新世初期。尹安(2001)認(rèn)為印度板塊與亞洲板塊之間最初碰撞可能開始于白堊紀(jì)晚期(約70 Ma)。李國彪等(2004)認(rèn)為印度–亞洲板塊碰撞啟動的時間可能在白堊紀(jì)與古近紀(jì)之交。由以上研究結(jié)果我們不難得出, 早白堊世晚期尼雄鐵礦形成時雅魯藏布江新特提斯洋殼正處于向北側(cè)拉薩地塊之下俯沖的早期階段, 符合尼雄鐵礦的形成環(huán)境, 因此,筆者認(rèn)為尼雄鐵礦的形成是雅魯藏布江洋殼板片向北側(cè)拉薩地塊之下俯沖所致。

    4.2 成礦巖漿源區(qū)

    尼雄鐵礦成礦花崗閃長巖屬中鉀–高鉀鈣堿性巖系列, 富堿、富集大離子親石元素(Rb、Th、U、K、Pb) 和輕稀土元素, 虧損重稀土和高場強(qiáng)元素(Nb、Ta、Ti), 顯示出島弧花崗巖的典型特征。巖體高的(87Sr/86Sr)i(0.707625~0.710997)、低的εNd(t) (–6.6~–8.7)、富集放射成因 Pb 同位素特征, 說明有較多地殼組分參與了成礦巖漿的生成。成礦花崗閃長巖Nb/U比值平均為6.09, MORB和OIB的Nb/U比值為47±(Hofmann et al., 1986), 兩者相差甚遠(yuǎn),而全球平均俯沖沉積物的Nb/U比值≈5(Plank and Langmuir, 1998), 與該巖體Nb/U比值接近, 因此推測尼雄鐵礦成礦巖漿生成過程中俯沖沉積物在地幔深部脫水和熔融產(chǎn)生的流體和熔體是造成巖漿源區(qū)EMⅡ型富集特征的首要因素。也就是說中生代早白堊世晚期雅魯藏布江洋殼向拉薩地塊之下低角度俯沖過程中, 來自俯沖板片的流體和熔體對地幔楔進(jìn)行交代改造, 使之具有了EMⅡ的源區(qū)特征, 并在流體和熔體作用下發(fā)生部分熔融形成了幔源巖漿。從圖10微量元素比值的相關(guān)變化(Ba-Nb/Y 和Th/Yb-Sr/Nd圖)可以看出尼雄鐵礦成礦巖漿生成過程中由俯沖板片釋放的流體和熔體起到了同等重要的作用。毫無疑問, 在島弧巖漿作用階段, 隨洋殼板片一起俯沖到深部的洋底沉積物無疑是地殼物質(zhì)參入造山帶巖漿作用最主要的方式, 然而, 也不能排除殼幔巖漿混合和含礦巖漿在上升過程中地殼物質(zhì)混染的可能性。對于成礦巖漿在上升過程中是否受到地殼物質(zhì)混染的情況, Briqueu and Lancelot (1979)指出可以利用ISr對1/Sr變異關(guān)系來進(jìn)行辯別(圖11a)。在巖漿上升過程中如果有地殼物質(zhì)混入, 在ISr-1/Sr圖上樣品呈正相關(guān)分布。圖11a中尼雄鐵礦的成礦花崗閃長巖沒有表現(xiàn)出正消長關(guān)系, 表明含礦巖漿在上升途中沒有受到地殼物質(zhì)混染。在εNd(t)-SiO2圖(圖11b)上, 隨著SiO2的增加尼雄花崗閃長巖的εNd(t)值沒有明顯的變化, 也說明其未遭受明顯的地殼物質(zhì)混染。而對于殼幔巖漿混合, 在(La/Yb)N-δEu圖中(圖6b)除了一個樣品落在殼源以外, 其他樣品都顯示殼?;煸吹奶卣? 在εNd(t)-ISr圖(圖6a)中, 樣品落在地幔玄武巖與下地殼的混合線附近, 也顯示出一種殼幔巖漿混合的特征。朱弟成等(2012)認(rèn)為尼雄鐵礦區(qū)所在的拉薩地塊中部結(jié)晶基底發(fā)育古元古代甚至太古宙微陸塊, 而成礦巖體Nd同位素模式年齡(tDM)集中在1207~1453 Ma(中元古代)之間, 也說明本地區(qū)的下地殼經(jīng)歷過早白堊世晚期幔源巖漿的底侵作用發(fā)生了部分熔融, 形成的殼源巖漿與幔源巖漿混合最后形成尼雄鐵礦的成礦巖漿。

    圖10 尼雄花崗閃長巖Ba-Nb/Y圖(a, 底圖據(jù)Kepezhinskas et al., 1997)和Th/Yb-Sr/Nd圖(b, 底圖據(jù)Woodhead et al., 1998)Fig.10 Ba vs Nb/Y (a) and Th/Yb vs Sr/Nd (b) diagrams for the granodiorites in the Nixiong iron deposit

    圖11 尼雄花崗閃長巖ISr-1/Sr (a, 據(jù)Briqueu and Lancelot, 1979)和εNd(t)-SiO2圖(b, 據(jù)張連昌等, 2007)Fig.11 I vs 1/Sr (a) and ε(t) vs SiO (b) diagrams for the granodiorites in the Nixiong iron deposit

    5 結(jié) 論

    (1) 尼雄鐵礦成礦花崗閃長巖為中鉀–高鉀鈣堿性巖, 亞鋁質(zhì), 是板塊俯沖作用的產(chǎn)物。巖體富集大離子不相容元素Rb、Ba、Th、U、K、Sr、Pb, 虧損高場強(qiáng)元素Nb、Ta、Ti, 具有島弧型巖漿作用的基本特征。稀土元素球粒隕石標(biāo)準(zhǔn)化分布曲線表現(xiàn)為輕稀土富集的右傾型, 無Eu異常, 顯示出活動大陸邊緣島弧巖漿作用特征。

    (2) 尼雄鐵礦成礦花崗閃長巖LA-ICPMS法測得的鋯石U-Pb年齡為112.09±0.54 Ma, 表明其形成于早白堊世晚期。

    (3) 綜合分析班公湖–怒江中特提斯洋和雅魯藏布江新特提斯洋的演化歷史, 得出尼雄鐵礦成礦構(gòu)造背景是雅魯藏布江洋殼板塊向北側(cè)拉薩地塊之下俯沖。

    (4) 尼雄鐵礦成礦花崗閃長巖的巖漿源區(qū)主要是俯沖組分(流體和熔融)對地幔楔交待改造的結(jié)果,地幔巖漿底侵引起下地殼物質(zhì)部分熔融, 兩種巖漿混合形成尼雄鐵礦的成礦巖漿。

    致謝: 中國地質(zhì)大學(xué)(北京)趙志丹教授和中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心李光明研究員提出了建設(shè)性修改意見和建議, 筆者在此表示衷心的感謝!

    曹圣華. 2012. 西藏尼雄式鐵礦及岡底斯中部鐵銅礦區(qū)域成礦規(guī)律研究. 北京: 中國地質(zhì)大學(xué)博士論文.

    曹圣華, 陳毓川, 樓法生. 2010. 西藏措勤縣尼雄矽卡巖型鐵(銅)礦田特征與找礦遠(yuǎn)景分析. 礦床地質(zhì), 29(S1): 65–66.

    曹圣華, 鄧世權(quán), 肖志堅(jiān), 廖六根. 2006. 班公湖–怒江結(jié)合帶西段中特提斯多島弧構(gòu)造演化. 沉積與特提斯地質(zhì), 26 (4): 25–32.

    曹圣華, 李德威, 余忠珍, 袁建芽, 吳旭鈴, 胡為正. 2007.西藏岡底斯尼雄超大型富鐵礦的成礦地質(zhì)特征. 大地構(gòu)造與成礦學(xué), 31(3): 328–334.

    陳能松, 王新宇, 張宏飛, 孫敏, 李曉彥, 陳強(qiáng). 2007. 柴–歐微地塊花崗巖地球化學(xué)和Nd-Sr-Pb同位素組成: 基底性質(zhì)和構(gòu)造屬性啟示. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào), 32(1): 7–21.

    高延林. 1985. 西藏南部雅魯藏布江縫合帶的板塊構(gòu)造標(biāo)志與演化. 西北大學(xué)學(xué)報(bào): 自然科學(xué)版, 15(4): 92–110.

    高延林. 1988. 雅魯藏布江縫合帶中段的板塊構(gòu)造證據(jù)及演化特征. 中國區(qū)域地質(zhì), 7(4): 53–59.

    郝杰, 柴育成, 李繼亮. 1999. 雅魯藏布江蛇綠巖的形成與日喀則弧前盆地沉積演化. 地質(zhì)科學(xué), 34(1): 4–12.

    侯增謙, 曲曉明, 王淑賢, 高永豐, 杜安道, 黃衛(wèi). 2003. 西藏高原岡底斯斑巖銅礦帶輝鉬礦 Re-Os 年齡: 成礦作用時限與動力學(xué)背景應(yīng)用. 中國科學(xué)(D輯), 33(7): 609–618.

    李光明, 秦克章, 陳雷, 陳金標(biāo), 范新, 琚宜太. 2011. 岡底斯東段山南地區(qū)第三紀(jì)矽卡巖–斑巖Cu-Mo-W(Au)多金屬礦床勘查模型及深部找礦意義. 地質(zhì)與勘探, 47(1): 20–30.

    李國彪, 萬曉樵, 劉文燦, 梁定益, H.Yun. 2004. 雅魯藏布江縫合帶南側(cè)古近紀(jì)海相地層的發(fā)現(xiàn)及其構(gòu)造意義. 中國科學(xué)(D輯), 34(3): 228–240.

    李獻(xiàn)華, 李武顯, 李正祥. 2007. 再論南嶺燕山早期花崗巖的成因類型與構(gòu)造意義. 科學(xué)通報(bào), 52( 9): 981–991.

    劉志飛, 王成善, 李祥輝. 2000. 西藏南部雅魯藏布江縫合帶的沉積–構(gòu)造演化. 同濟(jì)大學(xué)學(xué)報(bào): 自然科學(xué)版, 28(5): 010.

    潘桂棠, 莫宣學(xué), 侯增謙, 朱弟成, 王立全, 李光明, 趙志丹, 耿全如, 廖忠禮. 2006. 岡底斯造山帶的時空結(jié)構(gòu)及演化.巖石學(xué)報(bào), 22(3): 521–533.

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

    曲曉明, 辛洪波, 杜德道, 陳華. 2012. 西藏班公湖–怒江縫合帶中段碰撞后A–型花崗巖的時代及其對洋盆閉合時間的約束. 地球化學(xué), 41(1): 1–14.

    曲曉明, 辛洪波, 杜德道, 陳華. 2013. 西藏班公湖–怒江縫合帶中段 A–型花崗巖的巖漿源區(qū)與板片斷離. 地質(zhì)學(xué)報(bào), 87(6): 759–772.

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

    王樂, 何政偉, 劉婷婷, 倪忠云, 高慧, 蔡柯柯, 張船紅,吳華. 2012. 西藏尼雄銅鐵多金屬成礦帶遙感找礦模式與找礦方向. 現(xiàn)代地質(zhì), 26(3): 489–497.

    吳旭鈴, 陳振華. 2005. 西藏尼雄巖體巖石地球化學(xué)特征及其成因探討. 中國地質(zhì), 32(1): 122–127.

    吳元保, 鄭永飛. 2004. 鋯石成因礦物學(xué)研究及其對 U-Pb年齡解釋的制約. 科學(xué)通報(bào), 49(16): 1589–1604.

    徐向珍, 楊經(jīng)綏, 郭國林, 李金陽. 2011. 雅魯藏布江縫合帶西段普蘭蛇綠巖中地幔橄欖巖的巖石學(xué)研究.巖石學(xué)報(bào), 27(11): 3179–3196.

    尹安. 2001. 喜馬拉雅–青藏高原造山帶地質(zhì)演化——顯生宙亞洲大陸生長. 地球?qū)W報(bào), 22(3): 193–230.

    于玉帥, 楊竹森, 劉英超, 田世洪, 紀(jì)現(xiàn)華, 高原, 趙燦,趙武強(qiáng), 劉阿睢. 2012. 西藏措勤尼雄礦田滾糾鐵礦金云母礦物學(xué)特征及40Ar-39Ar年代學(xué). 巖石礦物學(xué)雜志, 31(5): 681–690.

    袁健芽, 曹圣華, 羅小川, 胡為正. 2008. 西藏措勤縣尼雄矽卡巖型鐵銅礦田的發(fā)現(xiàn)及地質(zhì)特征與找礦意義.中國地質(zhì), 35(1): 88–94.

    張連昌, 陳志廣, 周新華, 英基豐, 王非, 張玉濤. 2007.大興安嶺根河地區(qū)早白堊世火山巖深部源區(qū)與構(gòu)造–巖漿演化: Sr-Nd-Pb-Hf同位素地球化學(xué)制約. 巖石學(xué)報(bào), 23(11): 2823–2835.

    張萬平, 袁四化, 劉偉. 2011. 青藏高原南部雅魯藏布江蛇綠巖帶的時空分布特征及地質(zhì)意義. 西北地質(zhì), 44(1): 1–9.

    張曉倩, 朱弟成, 趙志丹, 王立全, 黃建村, 莫宣學(xué). 2010.西藏措勤尼雄巖體的巖石成因及其對富Fe成礦作用的潛在意義. 巖石學(xué)報(bào), 26(6): 1793–1804.

    張宇, 邵擁軍, 劉忠法, 周鑫, 鄭明泓. 2013. 安徽銅陵新橋銅硫鐵礦床地球化學(xué)特征及成因分析. 中國有色金屬學(xué)報(bào), 23(6): 1666–1680.

    朱弟成, 莫宣學(xué), 趙志丹, 牛耀齡, 潘桂棠, 王立全, 廖忠禮. 2009. 西藏南部二疊紀(jì)和早白堊世構(gòu)造巖漿作用與特提斯演化: 新觀點(diǎn). 地學(xué)前緣, 16(2): 1–20.

    朱弟成, 趙志丹, 牛耀齡, 王青, Yildirim DILEK, 董國臣,莫宣學(xué). 2012. 拉薩地體的起源和古生代構(gòu)造演化.高校地質(zhì)學(xué)報(bào), 18(1): 1–15.

    Batchelor R A and Bowden P. 1985. Petrogenetic interpretation of granitoid rock series using multicationic parameters. Chemical Geology, 48(1): 43–55.

    Boztug D, Harlavan Y, Arehart G B, Sat?r M and Avc? N. 2007. K-Ar age, whole-rock and isotope geochemistry of A-type granitoids in the Divri?i–Sivas region, eastern-central Anatolia, Turkey. Lithos, 97(1): 193–218.

    Briqueu L and Lancelot J R. 1979. Rb-Sr systematics and crustal contamination models for calc-alkaline igneous rocks. Earth and Planetary Science Letters, 43(3): 385–396.

    Collins W J, Beams S D, White A J R and Chappell B W. 1982. Nature and origin of A-type granites with particular reference to southeastern Australia. Contributions to Mineralogy and Petrology, 80(2): 189–200.

    DePaolo D J. 1981. Nd isotopic studies: Some new perspectives on earth structure and evolution. Eos, Transactions American Geophysical Union, 62(14): 137.

    DePaolo D J and Wasserburg G J. 1979. Petrogenetic mixing models and Nd-Sr isotopic patterns. Geochimica et Cosmochimica Acta, 43(4): 615–627.

    Garzanti E, Le Fort P and Sciunnach D. 1999. First report of Lower Permian basalts in South Tibet: Tholeiitic magmatism during break-up and incipient opening of Neotethys. Journal of Asian Earth Sciences, 17(4): 533–546.

    Hart S R. 1988. Heterogeneous mantle domains: Signatures, genesis and mixing chronologies. Earth and Planetary Science Letters, 90(3): 273–296.

    Hofmann A W. 1988. Chemical differentiation of the Earth: The relationship between mantle, continental crust, and oceanic crust. Earth and Planetary Science Letters, 90(3): 297–314.

    Hofmann A W, Jochum K P, Seufert M and White W M. 1986. Nb and Pb in oceanic basalts: New constraints on mantle evolution. Earth and Planetary Science Letters, 79(1): 33–45.

    Kapp P, Murphy M A, Yin A, Harrison T M, Ding L, and Guo J. 2003. Mesozoic and Cenozoic tectonic evolution of the Shiquanhe area of western Tibet. Tectonics, doi: 10.1029/2001TC001332.

    Kepezhinskas P, Mcdermott F, Defant M J, Hochstaedter A, Drummond M S, Hawkesworth C J, Koloskov A, Maury R C and Bellon H. 1997. Trace element and Sr-Nd-Pb isotopicconstraints on a three-component model of Kamchatka arc petrogenesis. Geochimica et Cosmochim Acta, 61(3): 577–600.

    Maniar P D and Piccoli P M. 1989. Tectonic discrimination of granitoids. Geological society of America bulletin, 101(5): 635–643.

    McDermid I R C, Aitchison J C, Davis A M, Harrison T M and Grove M. 2002. The Zedong terrane: A Late Jurassic intra-oceanic magmatic arc within the Yarlung-Tsangpo suture zone, southeastern Tibet. Chemical Geology, 187(3): 267–277.

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

    Pearce J A, Harris N B W and Tindle A G. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. Journal of petrology, 25(4): 956–983.

    Peccerillo A and Taylor S R. 1976. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, northern Turkey. Contributions to Mineralogy and Petrology, 58(1): 63–81.

    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): 325–394.

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

    Schmidberger S S and Hegner E. 1999. Geochemistry and isotope systematics of calc-alkaline volcanic rocks from the Saar-Nahe basin (SW Germany)-implications for Late-Variscan orogenic development. Contributions to Mineralogy and Petrology, 135(4): 373–385.

    Sui Q L, Wang Q, Zhu D C, Zhao Z D, Chen Y, Santosh M and Mo X X. 2013. Compositional diversity of ca. 110 Ma magmatism in the northern Lhasa Terrane, Tibet: Implications for the magmatic origin and crustal growth in a continent-continent collision zone. Lithos, 168: 144–159.

    Sun S S and McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implication for mantle composition and processes.

    Wilson B M. 1989. Igneous petrogenesis a global tectonic approach. Springer.

    Woodhead J D, Eggins S M and Johnson R W. 1998. Magma genesis in the New Britain island arc: Further insights into meltingand mass transfer processes. Journal of Petrology, 39(9): 1641–1668.

    Zartman R E and Doe B R. 1981. Plumbotectonics—the model. Tectonophysics, 75(1): 135–162.

    Zhu D C, Mo X X, Niu Y, Zhao Z D, Wang L Q, Liu Y S and Wu F Y. 2009. Geochemical investigation of Early Cretaceous igneous rocks along an east–west traverse throughout the central Lhasa Terrane, Tibet. Chemical Geology, 268(3): 298–312.

    Zhu D C, Zhao Z D, Niu Y, Mo X X, Chung S L, Hou Z Q and Wu F Y. 2011. The Lhasa Terrane: Record of a microcontinent and its histories of drift and growth. Earth and Planetary Science Letters, 301(1): 241–255.

    Zindler A and Hart S. 1986. Chemical geodynamics. Annual Review of Earth and Planetary Sciences, 14: 493–571.

    Petrogenesis of the Ore-forming Granodiorite in the Nixiong Iron Deposit and its Implications for the Metallogenic Tectonic Background

    FAN Shufang, QU Xiaoming, SONG Yang and XIN Hongbo
    (Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 100037, China)

    Whole-rock major, trace and rare earth elements and Sr-Nd-Pb isotopic analyses, combined with zircon LA-ICPMS U-Pb dating, of granodiorite related to the Nixiong iron deposit, are carried out in this paper. The petrogeochemical analysis results indicate that the granodiorite related to Nixiong iron deposit belongs to subaluminous, medium-high K calc-alkaline I-type granite. Geochemically, the granodiorite is enriched in large ion lithophile elements(LILE) Rb, Ba, Sr, Th, U, K, Pb and depleted in high field strength elements(HFSE) Nb, Ta, Ti and shows basic characteristics of arc magmatic rocks. The rare earth element chondrite-normalized distribution patterns appear a notable enrichment of light rare earth elements relative to heavy rare earth elements without marked Eu anomalies and show signatures of magmatic rocks from an active continental margin. Zircon U-Pb LA-ICPMS dating yields a weighted206Pb/238U mean age of 112.09±0.54 Ma with MSWD=0.45, indicating that the granodiorite was formed in late stage of Early Cretaceous. These granites have high radiogenic Sr isotopic ratios (ISr=0.707625–0.710997) and low Nd isotopic ratios εNd(t)(–6.6 to –8.7) and primitive Pb isotopic compositions (206Pb/204Pb=18.786–18.955,207Pb/204Pb=15.694–15.726,208Pb/204Pb=39.355–39.676), reflecting an important influence of crustal components on the magmatic generation. Through a comprehensive analysis on the evolution histories of the Bangonghu-Nujiang Meso-Tethyan ocean and the Yaluzangbo Neo-Tethyan ocean, the authors of the paper reach a conclusion that the Nixiong iron deposit is genetically related to the northward subduction of the Yarlu-Zangbo oceanic crust. The enrichment component of fluid and melt released from sediments adhered to the subducted oceanic crust triggered the mantle wedge partially melt to generate the enriched magma. The mantle magma intruded and caused partical melting of lower crust and mixed with each other to form the metallogenic magma.

    petrogeochemistry; Sr-Nd-Pb isotopics; zircon U-Pb LA-ICPMS dating; Nixiong iron deposit; Xizang

    P613; P597

    A

    1001-1552(2015)02-0286-014

    2013-05-09; 改回日期: 2013-08-03

    項(xiàng)目資助: 中國地質(zhì)調(diào)查局青藏專項(xiàng)工作項(xiàng)目(編號: 12120113037300)和國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)(編號: 2011CB403103)聯(lián)合資助。第一作者簡介: 范淑芳(1988–), 女, 碩士研究生, 礦床地球化學(xué)專業(yè)。Email: fsf1005@126.com

    曲曉明(1960–), 男, 研究員, 地球化學(xué)專業(yè), 主要從事造山帶成礦學(xué)研究。Email: xiaominqu@163.com

    猜你喜歡
    洋殼班公湖雅魯藏布江
    洋殼活動斷層特征和動力作用
    給青藏高原的班公湖量體溫
    軍事文摘(2021年20期)2021-11-10 01:58:40
    印度洋、大西洋洋殼流運(yùn)動對地形地貌的影響
    巖漿與構(gòu)造作用對洋殼厚度的影響
    ——以西北印度洋為例
    在班公湖,與天堂偶遇
    太平洋洋殼流運(yùn)動對地形地貌的影響
    中國與孟加拉國在雅魯藏布江河流治理中的合作與問題探究
    商(2016年32期)2016-11-24 15:28:37
    西藏雅魯藏布江區(qū)橋梁水文計(jì)算
    西藏科技(2015年2期)2015-09-26 12:18:32
    雅魯藏布江特大橋的臨時鋼棧橋結(jié)構(gòu)計(jì)算
    西藏科技(2015年3期)2015-09-26 12:11:11
    大美尼洋河
    亚洲欧美精品综合久久99| 久久婷婷人人爽人人干人人爱| 久久精品久久久久久噜噜老黄 | 热99re8久久精品国产| 午夜精品在线福利| 成人无遮挡网站| 偷拍熟女少妇极品色| 日韩成人av中文字幕在线观看 | 男女做爰动态图高潮gif福利片| 日日啪夜夜撸| 国产真实伦视频高清在线观看| 99久久中文字幕三级久久日本| 91久久精品国产一区二区三区| 国产黄色小视频在线观看| 一本久久中文字幕| 国产中年淑女户外野战色| 蜜桃久久精品国产亚洲av| 欧美绝顶高潮抽搐喷水| 天堂影院成人在线观看| 麻豆一二三区av精品| 欧美一级a爱片免费观看看| 可以在线观看毛片的网站| 白带黄色成豆腐渣| 99久久九九国产精品国产免费| 老熟妇仑乱视频hdxx| 国产av在哪里看| av国产免费在线观看| av国产免费在线观看| 国产精品,欧美在线| 久久精品影院6| 哪里可以看免费的av片| 男女边吃奶边做爰视频| АⅤ资源中文在线天堂| 女人十人毛片免费观看3o分钟| 国产真实乱freesex| 直男gayav资源| 国产精品亚洲美女久久久| 国产乱人视频| 变态另类丝袜制服| 最近的中文字幕免费完整| 成人综合一区亚洲| 12—13女人毛片做爰片一| 久久久久久久久久黄片| 中国美女看黄片| 成年av动漫网址| 免费看光身美女| 国产午夜精品久久久久久一区二区三区 | 男人舔女人下体高潮全视频| 亚洲精品影视一区二区三区av| 成熟少妇高潮喷水视频| 国产视频一区二区在线看| 男女做爰动态图高潮gif福利片| 少妇猛男粗大的猛烈进出视频 | 国产高清不卡午夜福利| 99热只有精品国产| 亚洲人与动物交配视频| 色5月婷婷丁香| 亚洲真实伦在线观看| 欧美性猛交黑人性爽| 美女 人体艺术 gogo| 日韩一区二区视频免费看| 我的老师免费观看完整版| 深夜精品福利| 秋霞在线观看毛片| 亚洲成av人片在线播放无| 免费看a级黄色片| 久久久精品欧美日韩精品| 中文字幕精品亚洲无线码一区| a级一级毛片免费在线观看| 成人av一区二区三区在线看| 中文字幕熟女人妻在线| 91精品国产九色| 久久精品国产亚洲av涩爱 | 在线观看美女被高潮喷水网站| 久久久久性生活片| 黄色配什么色好看| 国产精品国产三级国产av玫瑰| 99久国产av精品国产电影| 亚洲人成网站在线播| 国产精品国产高清国产av| 精品久久久久久久久亚洲| 日韩av在线大香蕉| 午夜久久久久精精品| 婷婷精品国产亚洲av| 日韩av不卡免费在线播放| 一级毛片电影观看 | 免费大片18禁| 欧美色欧美亚洲另类二区| 免费一级毛片在线播放高清视频| 久久人人精品亚洲av| 亚洲最大成人av| 99热全是精品| 人人妻人人澡人人爽人人夜夜 | 精品午夜福利视频在线观看一区| 黄色一级大片看看| 搡女人真爽免费视频火全软件 | 亚洲国产精品合色在线| 最新中文字幕久久久久| 亚洲精品久久国产高清桃花| 国产在视频线在精品| 熟女电影av网| 97在线视频观看| 日韩 亚洲 欧美在线| 可以在线观看毛片的网站| 波多野结衣高清无吗| 国产精品精品国产色婷婷| 97超级碰碰碰精品色视频在线观看| 国产成人一区二区在线| 精品国产三级普通话版| 麻豆国产av国片精品| 在线播放国产精品三级| 欧美日韩国产亚洲二区| 婷婷六月久久综合丁香| 免费看美女性在线毛片视频| 性欧美人与动物交配| 久久精品国产亚洲av涩爱 | 久久久久国内视频| 看片在线看免费视频| 国产精品福利在线免费观看| 亚洲婷婷狠狠爱综合网| 欧美日韩精品成人综合77777| 深爱激情五月婷婷| 亚洲欧美精品综合久久99| 一级黄色大片毛片| 少妇丰满av| 欧美性感艳星| 国产美女午夜福利| 99热精品在线国产| 此物有八面人人有两片| a级毛片免费高清观看在线播放| 一本一本综合久久| 99久国产av精品| 啦啦啦观看免费观看视频高清| 日韩国内少妇激情av| 国内少妇人妻偷人精品xxx网站| 99热网站在线观看| 日韩成人av中文字幕在线观看 | 国内少妇人妻偷人精品xxx网站| 日韩欧美一区二区三区在线观看| 深爱激情五月婷婷| 亚洲欧美中文字幕日韩二区| 亚洲第一电影网av| 日本熟妇午夜| 日韩人妻高清精品专区| 国产精品三级大全| 日本黄色片子视频| 国产老妇女一区| 美女黄网站色视频| 99热这里只有是精品在线观看| 国产伦一二天堂av在线观看| 国产精品国产三级国产av玫瑰| 国产麻豆成人av免费视频| av.在线天堂| 日日撸夜夜添| 少妇熟女欧美另类| 一个人看的www免费观看视频| 免费av毛片视频| 国产美女午夜福利| 精品不卡国产一区二区三区| 偷拍熟女少妇极品色| 亚洲不卡免费看| 中文亚洲av片在线观看爽| 长腿黑丝高跟| 性欧美人与动物交配| 九九在线视频观看精品| 国产黄色小视频在线观看| 国产综合懂色| 人人妻人人澡欧美一区二区| 一本久久中文字幕| 91在线观看av| 热99在线观看视频| 51国产日韩欧美| 国产精品亚洲美女久久久| 免费av毛片视频| 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添av毛片| 亚洲av美国av| 亚洲一区高清亚洲精品| 高清午夜精品一区二区三区 | 能在线免费观看的黄片| 老熟妇乱子伦视频在线观看| 中国国产av一级| 尤物成人国产欧美一区二区三区| av天堂在线播放| 卡戴珊不雅视频在线播放| 国产aⅴ精品一区二区三区波| 此物有八面人人有两片| 国内精品宾馆在线| 久久亚洲国产成人精品v| 热99在线观看视频| 成人鲁丝片一二三区免费| .国产精品久久| 丰满乱子伦码专区| 国产又黄又爽又无遮挡在线| 亚洲精品日韩av片在线观看| 久久天躁狠狠躁夜夜2o2o| 中文字幕久久专区| 日本色播在线视频| 国产人妻一区二区三区在| 久久精品国产亚洲av涩爱 | 国产淫片久久久久久久久| 欧美日韩国产亚洲二区| 国产成人福利小说| 国产色爽女视频免费观看| 最好的美女福利视频网| 美女cb高潮喷水在线观看| 91在线观看av| 日韩三级伦理在线观看| 成年免费大片在线观看| 成人av一区二区三区在线看| 久久韩国三级中文字幕| 欧美高清成人免费视频www| 国产成人91sexporn| 三级毛片av免费| 成年版毛片免费区| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 日韩av不卡免费在线播放| 18禁在线无遮挡免费观看视频 | 男女之事视频高清在线观看| 美女被艹到高潮喷水动态| 3wmmmm亚洲av在线观看| 99热这里只有是精品在线观看| 99视频精品全部免费 在线| 男女边吃奶边做爰视频| 又黄又爽又刺激的免费视频.| 精品一区二区三区人妻视频| 69人妻影院| 乱码一卡2卡4卡精品| 国产一区二区三区av在线 | 国产又黄又爽又无遮挡在线| 一区福利在线观看| 可以在线观看的亚洲视频| 国产精品1区2区在线观看.| 女生性感内裤真人,穿戴方法视频| 亚洲精品456在线播放app| 欧美日韩一区二区视频在线观看视频在线 | 久久99热这里只有精品18| 精品人妻熟女av久视频| .国产精品久久| 黄色视频,在线免费观看| 乱码一卡2卡4卡精品| 亚洲精品日韩av片在线观看| av女优亚洲男人天堂| 日韩精品中文字幕看吧| 国产精品美女特级片免费视频播放器| h日本视频在线播放| 97在线视频观看| av黄色大香蕉| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 国产午夜福利久久久久久| 久久6这里有精品| av专区在线播放| 99久久精品热视频| 精品久久久久久成人av| 亚洲成人中文字幕在线播放| 1024手机看黄色片| 成人特级黄色片久久久久久久| 免费观看人在逋| 国产欧美日韩精品一区二区| 欧美精品国产亚洲| 国产精品永久免费网站| 亚洲欧美日韩无卡精品| 午夜激情福利司机影院| 国产高清不卡午夜福利| 99久久无色码亚洲精品果冻| 内射极品少妇av片p| 麻豆国产97在线/欧美| 久久精品国产亚洲av天美| 亚洲自拍偷在线| 国产片特级美女逼逼视频| 白带黄色成豆腐渣| 99国产精品一区二区蜜桃av| 日本熟妇午夜| 成人漫画全彩无遮挡| av黄色大香蕉| 国产精品av视频在线免费观看| 一级黄片播放器| 寂寞人妻少妇视频99o| 亚洲av熟女| 久久久欧美国产精品| 亚洲18禁久久av| 九九爱精品视频在线观看| 蜜臀久久99精品久久宅男| 禁无遮挡网站| 久久久久久久久久成人| 日韩强制内射视频| 可以在线观看毛片的网站| 亚洲av.av天堂| 午夜激情福利司机影院| 嫩草影视91久久| 99久久精品国产国产毛片| 一个人看的www免费观看视频| 观看美女的网站| 精品久久久噜噜| 亚洲综合色惰| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 给我免费播放毛片高清在线观看| 小说图片视频综合网站| 午夜福利在线观看吧| 联通29元200g的流量卡| 日本与韩国留学比较| 一个人免费在线观看电影| 久久精品国产亚洲av天美| 一进一出抽搐gif免费好疼| 国产精品一及| 久久久久精品国产欧美久久久| 人人妻,人人澡人人爽秒播| 成人性生交大片免费视频hd| 国产成人精品久久久久久| 最后的刺客免费高清国语| 看免费成人av毛片| 免费观看人在逋| 97超碰精品成人国产| 亚洲国产欧洲综合997久久,| 国产aⅴ精品一区二区三区波| 国产精品av视频在线免费观看| 国产v大片淫在线免费观看| 人人妻人人澡欧美一区二区| 日本黄色片子视频| 日本与韩国留学比较| 美女 人体艺术 gogo| 12—13女人毛片做爰片一| 毛片女人毛片| 麻豆久久精品国产亚洲av| 日本一二三区视频观看| 亚洲专区国产一区二区| 成人午夜高清在线视频| 白带黄色成豆腐渣| 久久天躁狠狠躁夜夜2o2o| 久久鲁丝午夜福利片| 少妇的逼好多水| 人妻丰满熟妇av一区二区三区| 日韩欧美三级三区| 国产 一区 欧美 日韩| 亚洲美女视频黄频| 亚洲精品乱码久久久v下载方式| 亚洲精品日韩av片在线观看| 成熟少妇高潮喷水视频| 日本 av在线| 欧美极品一区二区三区四区| 亚洲激情五月婷婷啪啪| 亚洲国产色片| 大型黄色视频在线免费观看| 国产精品美女特级片免费视频播放器| 男女啪啪激烈高潮av片| 国产麻豆成人av免费视频| 午夜福利高清视频| 亚洲最大成人手机在线| 一区二区三区高清视频在线| 亚洲国产色片| 亚洲精品亚洲一区二区| 97超级碰碰碰精品色视频在线观看| 久久欧美精品欧美久久欧美| 麻豆成人午夜福利视频| 精品人妻偷拍中文字幕| 最近的中文字幕免费完整| 日韩人妻高清精品专区| 久久久久性生活片| 国产真实伦视频高清在线观看| 亚洲一区高清亚洲精品| 精品久久久久久久久亚洲| 亚洲欧美精品综合久久99| 亚洲精品成人久久久久久| 极品教师在线视频| 精华霜和精华液先用哪个| 中出人妻视频一区二区| 午夜a级毛片| 成熟少妇高潮喷水视频| 国产蜜桃级精品一区二区三区| 亚洲国产精品国产精品| 深夜精品福利| 99久久精品一区二区三区| 变态另类丝袜制服| 欧美高清成人免费视频www| 舔av片在线| 91狼人影院| 嫩草影院精品99| 国产精品人妻久久久久久| 丰满乱子伦码专区| 国产精品亚洲美女久久久| 观看美女的网站| 天天一区二区日本电影三级| 六月丁香七月| 国产午夜福利久久久久久| 国产欧美日韩一区二区精品| 欧美+亚洲+日韩+国产| 少妇熟女欧美另类| 搡老岳熟女国产| 一卡2卡三卡四卡精品乱码亚洲| av国产免费在线观看| 成人亚洲精品av一区二区| 精品久久久久久久久亚洲| 免费黄网站久久成人精品| 国产蜜桃级精品一区二区三区| 国产在视频线在精品| 97超视频在线观看视频| 国产精品野战在线观看| 久久精品综合一区二区三区| 中国美白少妇内射xxxbb| 欧美性猛交╳xxx乱大交人| 日日摸夜夜添夜夜添av毛片| 欧美zozozo另类| 亚洲av五月六月丁香网| videossex国产| 网址你懂的国产日韩在线| 亚洲婷婷狠狠爱综合网| 寂寞人妻少妇视频99o| 人妻少妇偷人精品九色| 午夜影院日韩av| 女的被弄到高潮叫床怎么办| 99久久精品热视频| 国产色婷婷99| 亚洲精品国产成人久久av| 日韩欧美在线乱码| 一个人免费在线观看电影| 国产精品国产高清国产av| 18禁在线无遮挡免费观看视频 | 一级黄片播放器| 国产精品无大码| 亚洲欧美精品综合久久99| 国产伦精品一区二区三区视频9| 直男gayav资源| 日韩三级伦理在线观看| 少妇的逼水好多| 在线观看美女被高潮喷水网站| 婷婷亚洲欧美| 国内精品宾馆在线| 欧美又色又爽又黄视频| 午夜精品国产一区二区电影 | 亚洲精品国产成人久久av| 亚洲,欧美,日韩| 国产高清有码在线观看视频| 国产亚洲精品久久久久久毛片| 美女cb高潮喷水在线观看| 老熟妇仑乱视频hdxx| 国产视频内射| 少妇被粗大猛烈的视频| 国产高清三级在线| 淫秽高清视频在线观看| 久久久久国产网址| 一级毛片我不卡| 你懂的网址亚洲精品在线观看 | 亚洲婷婷狠狠爱综合网| 两个人的视频大全免费| 亚洲一区高清亚洲精品| 亚洲综合色惰| 真人做人爱边吃奶动态| 午夜亚洲福利在线播放| 国产精品亚洲一级av第二区| 欧美xxxx性猛交bbbb| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| 久久久久久久久大av| 成人一区二区视频在线观看| 综合色丁香网| 色综合色国产| 禁无遮挡网站| av天堂在线播放| 午夜精品国产一区二区电影 | 超碰av人人做人人爽久久| 又黄又爽又免费观看的视频| 给我免费播放毛片高清在线观看| 人人妻,人人澡人人爽秒播| 人妻久久中文字幕网| www.色视频.com| 国产乱人偷精品视频| 天堂动漫精品| 青春草视频在线免费观看| a级毛片免费高清观看在线播放| 免费看光身美女| 俺也久久电影网| 天天躁夜夜躁狠狠久久av| 久久精品国产亚洲av涩爱 | 精品久久久久久成人av| 波多野结衣高清作品| 国产精品不卡视频一区二区| 91久久精品国产一区二区成人| 国产老妇女一区| 亚洲成人久久爱视频| 91av网一区二区| 亚洲国产精品久久男人天堂| 国产男人的电影天堂91| 欧美绝顶高潮抽搐喷水| 国产亚洲av嫩草精品影院| 99国产极品粉嫩在线观看| 蜜桃久久精品国产亚洲av| 久久精品影院6| 欧美激情国产日韩精品一区| 男插女下体视频免费在线播放| 精品人妻偷拍中文字幕| 人妻久久中文字幕网| 亚洲18禁久久av| 国产成人一区二区在线| 国产成人影院久久av| 亚洲精品日韩在线中文字幕 | 三级经典国产精品| 好男人在线观看高清免费视频| 日日啪夜夜撸| 又爽又黄a免费视频| 欧美性猛交黑人性爽| 伦理电影大哥的女人| 国产淫片久久久久久久久| 少妇人妻一区二区三区视频| 直男gayav资源| 97在线视频观看| 99久久九九国产精品国产免费| 1024手机看黄色片| 97人妻精品一区二区三区麻豆| 精品人妻一区二区三区麻豆 | 又黄又爽又刺激的免费视频.| 嫩草影院精品99| 国产高清不卡午夜福利| 久久人人爽人人片av| avwww免费| 淫妇啪啪啪对白视频| 听说在线观看完整版免费高清| av福利片在线观看| 少妇熟女欧美另类| 91麻豆精品激情在线观看国产| 天堂av国产一区二区熟女人妻| 免费人成视频x8x8入口观看| 国内少妇人妻偷人精品xxx网站| 欧美成人免费av一区二区三区| 国产亚洲91精品色在线| 少妇熟女aⅴ在线视频| 禁无遮挡网站| 美女大奶头视频| 久久人妻av系列| 激情 狠狠 欧美| .国产精品久久| 99热6这里只有精品| 日本 av在线| 成人精品一区二区免费| videossex国产| 日本五十路高清| 国产aⅴ精品一区二区三区波| 欧美性猛交╳xxx乱大交人| 国产视频一区二区在线看| 99热这里只有是精品在线观看| 高清日韩中文字幕在线| 国产伦在线观看视频一区| 精品久久久久久久末码| 日本五十路高清| 国产又黄又爽又无遮挡在线| 免费观看人在逋| 国内精品美女久久久久久| 国产69精品久久久久777片| 国产私拍福利视频在线观看| 国产伦一二天堂av在线观看| 一本一本综合久久| 久久精品国产自在天天线| 青春草视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区精品小视频在线| 国产v大片淫在线免费观看| 国产国拍精品亚洲av在线观看| 国产精品国产高清国产av| 国产麻豆成人av免费视频| 丰满的人妻完整版| 看黄色毛片网站| 插阴视频在线观看视频| 亚洲欧美日韩无卡精品| 99精品在免费线老司机午夜| 国产精品女同一区二区软件| 亚洲电影在线观看av| 中出人妻视频一区二区| 国产精品三级大全| 一级黄色大片毛片| 亚洲一区高清亚洲精品| 国产国拍精品亚洲av在线观看| 亚洲色图av天堂| 麻豆乱淫一区二区| 亚洲18禁久久av| 伦理电影大哥的女人| 国产中年淑女户外野战色| 亚洲一级一片aⅴ在线观看| 搡老熟女国产l中国老女人| 久久午夜亚洲精品久久| 日韩精品有码人妻一区| 亚洲成av人片在线播放无| 国内精品宾馆在线| avwww免费| 国产亚洲精品av在线| 最近手机中文字幕大全| 国产真实伦视频高清在线观看| 欧美潮喷喷水| 嫩草影院入口| 亚洲精品在线观看二区| 老师上课跳d突然被开到最大视频| 嫩草影院新地址| 91久久精品电影网| 高清毛片免费看| 免费看美女性在线毛片视频| 麻豆久久精品国产亚洲av| 日本三级黄在线观看| 激情 狠狠 欧美| 三级国产精品欧美在线观看| 久久精品国产亚洲av天美| 日韩欧美免费精品| 人妻夜夜爽99麻豆av| 国产精品永久免费网站| 日韩成人av中文字幕在线观看 | 中文字幕免费在线视频6| 激情 狠狠 欧美| 我的老师免费观看完整版| 色尼玛亚洲综合影院| 中国国产av一级| 国产女主播在线喷水免费视频网站 | 日韩av不卡免费在线播放| 观看免费一级毛片| 免费人成视频x8x8入口观看| 2021天堂中文幕一二区在线观|