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

    長江中下游成礦帶龍角山矽卡巖鎢礦床成礦作用過程
    ——來自白鎢礦和石榴子石主微量元素的證據(jù)*

    2022-03-17 00:35:38聶利青周濤發(fā)蔡國軍孫孝峰宋玉龍王鳳云
    礦床地質(zhì) 2022年1期
    關(guān)鍵詞:子石白鎢礦矽卡巖

    聶利青,周濤發(fā),蔡國軍,孫孝峰,宋玉龍,蔡 毅,殷 帥,王鳳云

    (1安徽建筑大學(xué)土木工程學(xué)院智能地下探測技術(shù)重點(diǎn)試驗(yàn)室,安徽合肥 230601;2合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院,安徽合肥230009;3湖北省地質(zhì)局第一地質(zhì)大隊(duì),湖北大冶 435000)

    白鎢礦和石榴子石作為矽卡巖型礦床中普遍存在的金屬礦物和硅酸鹽礦物(Meinert et al.,2005;Goldmann et al.,2013;Poulin et al.,2016;Guo et al.,2016;Xiong et al.,2017;Sciuba et al.,2019;Liu et al.,2020;Zhang et al.,2020),隨著礦物原位分析(in-situ analyses)的不斷進(jìn)步和激光剝蝕電感耦合等離子質(zhì)譜儀(LA-ICP-MS)技術(shù)的推廣(Crowe et al.,2001;周 濤 發(fā) 等,2010;Reich et al.,2013;Huang et al.,2015;汪方躍等,2017),白鎢礦和石榴子石的原位微量元素特征被廣泛應(yīng)用于解釋成礦流體物理化學(xué)條件和組分,進(jìn)而反映成礦流體演化和水巖反應(yīng)程度(Ghaderi et al.,1999;Brugger et al.,2000;Song et al.,2014;Sun et al.,2017;Zhao et al.,2018;Li et al.,2018;Zhang et al.,2018;Harlaux et al.,2018;Wu et al.,2019;Han et al.,2020)。除巖相學(xué)特征外,陰極發(fā)光(CL)也是揭示礦物環(huán)帶特征的強(qiáng)有力技術(shù)手段,并且已經(jīng)成功應(yīng)用于多種礦物(例如石英、磷灰石、碳酸鹽礦物;McManuset al.,1992;Campbell et al.,1997;Gotze et al.,2002;Redmond et al.,2004;Rusk et al.,2008;Müller et al.,2010;Poulin et al.,2016),并用于指示成礦作用。特別是CL圖像直接與激光原位面掃描圖像和微量元素結(jié)合闡述礦物復(fù)雜的振蕩環(huán)帶,記錄矽卡巖成礦系統(tǒng)中脈石礦物和礦石礦物結(jié)晶時(shí)熱液流體的組成、性質(zhì)和演化(Cot‐trant,1981;Gaft et al.,1999;Brugger et al.,2000;Mair,2006;Roberts et al.,2006;MacRae et al.,2009;Zhaiet al.,2014;Poulin et al.,2016),并反演熱液礦床中流體特征和為成礦動(dòng)力學(xué)背景提供重要信息(Sverjensky,1984;Giere,1996;Uspensky et al.,1998;Ghaderi et al.,1999;Brugger et al.,2000;Hazarika et al.,2016;Qi et al.,2020;Su et al.,2021)。

    長江中下游成礦帶是中國礦床學(xué)研究的熱點(diǎn)地區(qū)之一(周濤發(fā)等,2012),也是“層控矽卡巖礦床”(常印佛等,1991)和“斑巖-矽卡巖復(fù)合成礦理論”(翟裕生等,1992)的發(fā)祥地。近年來在長江中下游成礦帶新發(fā)現(xiàn)的鎢礦床為該成礦帶的研究提出了新的課題(周濤發(fā)等,2019)。龍角山礦床是成礦帶鄂東南礦集區(qū)內(nèi)新發(fā)現(xiàn)的大型矽卡巖型鎢多金屬礦床。鄂東南礦集區(qū)從南東到北西,礦化類型呈鎢—鎢銅—銅鐵—銅—鐵扇形過渡,是研究長江中下游成礦帶鎢與銅、鐵成礦作用關(guān)系的理想對象。前人對龍角山礦床的研究主要集中于成巖成礦年代和巖體地球化學(xué)特征(丁麗雪等,2014)、穩(wěn)定同位素和流體包裹體(Lei et al.,2018)以及礦物電子探針主量元素分析(紀(jì)云昊等,2019)等方面,但礦床成礦作用的環(huán)境和礦質(zhì)沉淀機(jī)制的研究則相對薄弱。本文在礦床地質(zhì)特征和前人研究工作基礎(chǔ)上,對龍角山礦床的白鎢礦和石榴子石開展了系統(tǒng)的礦物學(xué)和原位微區(qū)地球化學(xué)特征研究,試圖厘定該礦床的成礦流體演化特征,分析鎢的沉淀機(jī)制,提高成礦帶鎢礦床成礦作用的研究程度并為找礦勘探提供理論依據(jù)。

    1 礦床地質(zhì)特征

    1.1 區(qū)域地質(zhì)背景

    長江中下游成礦帶內(nèi)金屬礦床(點(diǎn))達(dá)200余處,主要集中在八大礦集區(qū)(常印佛等,1991;周濤發(fā)等,2017;圖1a、b),成礦帶內(nèi)的礦床主要有4種類型:①矽卡巖-斑巖型銅-金-鎢礦床,其成礦作用與143~140 Ma的花崗閃長巖類侵入巖體關(guān)系密切,主要分布在銅陵礦集區(qū)、鄂東南礦集區(qū)、九瑞礦集區(qū)、安慶-貴池礦集區(qū)和宣城礦集區(qū)(Mao et al.,2006;Xie et al.,2007;王世偉等,2012;謝桂青等,2013;徐曉春等,2014;Xiao et al.,2021);②矽卡巖型鎢鉬礦床,其成礦作用與135 Ma的花崗閃長巖類侵入體關(guān)系密切,主要分布在皖南地區(qū)(宋國學(xué)等,2010;丁寧,2012);③玢巖型鐵(鎢)礦床(化)(磁鐵礦-磷灰石型礦床),成礦作用與130 Ma左右的輝長閃長玢巖關(guān)系密切,主要分布在寧蕪礦集區(qū)、廬樅礦集區(qū)和鄂東南礦集區(qū)(張樂駿等,2011;范裕等,2011;Xie et al.,2015;Nie et al.,2017);④與酸性巖相關(guān)的金(鈾)-鎢礦床,成礦時(shí)代為108~97 Ma,主要分布在寧鎮(zhèn)礦集區(qū)和廬樅礦集區(qū)(孫洋等,2014;聶利青等,2016a;張贊贊等,2018;Zhang et al.,2021)。前人對上述4類礦床開展了系統(tǒng)的研究工作,建立了多層樓矽卡巖模式(常印佛等,1991)、玢巖鐵礦模式(寧蕪玢巖鐵礦編寫組,1978)等,并指出燕山期的陸內(nèi)俯沖是造成長江中下游成礦帶大規(guī)模成巖和成礦作用的主導(dǎo)機(jī)制(Lüet al.,2015;2021)。近年來,一系列鎢礦床(化)在長江中下游成礦帶內(nèi)被相繼發(fā)現(xiàn)(周濤發(fā)等,2019),“南鎢北擴(kuò)”界限逐步向北擴(kuò)展,鎢礦床可獨(dú)立產(chǎn)出,也可與成礦帶內(nèi)各類金屬礦化共、伴生,成礦時(shí)代跨度大,成礦類型均為矽卡巖型(宋國學(xué)等,2014;顏代蓉等,2012;聶利青等,2016b;2018;陳雪鋒等,2017),如湖北省大冶市銅山口銅礦床、安徽省廬江縣龍橋鐵礦床、安徽省銅陵市姚家?guī)X鋅金礦床均發(fā)現(xiàn)共、伴生鎢礦化(朱喬喬等,2014;張維,2017;鐘國雄等,2014)。

    圖1 長中下游成礦帶礦床分布示意圖(據(jù)常印佛等,1991;周濤發(fā)等,2017)a.中國地質(zhì)略圖;b.長江中下游成礦帶地質(zhì)圖XGF—襄樊-廣濟(jì)斷裂;TLF—郯廬斷裂;HPF—黃-破斷裂;SMF—商麻斷裂;CCF—崇陽-常州斷裂;CHF—滁河斷裂;JNF—江南斷裂Fig.1 Geologic map of magmatic rocks and deposits in the Middle-Lower Yangtze River Metallogenic Belt(MLYB)(after Chang et al.,1991;Zhou et al.,2017)a.Geology sketch of China;b.Geologic map of MLYB XGF—Xiangfan-Guangji rupture;TLF—Tancheng-Lujiang rupture;HPF—Huang-Po rupture;SMF—ShangMa rupture;CCF—Chongyang-Changzhou rupture;CHF—ChuHe rupture;JNF—JiangNan rupture

    鄂東南礦集區(qū)發(fā)育有著名的“大冶式”大型富鐵矽卡巖型鐵礦床(如鐵山、程潮、金山店)、斑巖-矽卡巖復(fù)合型銅鉬礦床(銅山口和豐山洞)、矽卡巖型銅金礦床(雞冠嘴和雞籠山)、矽卡巖型銅鎢礦床(龍角山阮家灣)和中國最大的矽卡巖型銅鐵礦床(銅綠山),這些礦床均屬于與燕山期中酸性侵入巖有關(guān)的成礦系列(謝桂青等,2008)。根據(jù)地球化學(xué)元素分布特點(diǎn),鄂東南礦集區(qū)分為4個(gè)成礦地球化學(xué)區(qū),呈扇形展布,即鄂城-金山店-靈鄉(xiāng)親鐵元素地球化學(xué)區(qū)、黃石-大冶親鐵元素和親銅元素地球化學(xué)區(qū)、陽新親銅元素地球化學(xué)區(qū)、殷祖-豐山親銅元素和鎢鉬族元素地球化學(xué)區(qū)(舒全安等,1992)。

    1.2 礦床地質(zhì)特征

    龍角山鎢礦床位于鄂東礦集區(qū)陽新巖體西側(cè)(圖2),是鄂東南礦集區(qū)內(nèi)新發(fā)現(xiàn)的大型鎢礦床。龍角山礦區(qū)于1956~1972年進(jìn)行補(bǔ)充勘探和深部地質(zhì)找礦工作,2006年、2010~2012年湖北省鄂東南地質(zhì)大隊(duì)對龍角山礦區(qū)的港溝山礦段和港背山礦段進(jìn)行銅鉬鎢礦普查地質(zhì)工作,2018~2019年湖北省鄂東南地質(zhì)大隊(duì)對龍角山礦區(qū)進(jìn)行外圍找礦工作,現(xiàn)已探明礦床鎢資源量5.84萬噸,銅資源5.5萬噸(湖北省第一地質(zhì)大隊(duì),2018),目前該礦床鎢資源量已達(dá)到大型規(guī)模,進(jìn)一步的勘探工作正在進(jìn)行,有望擴(kuò)大資源量。

    圖2 鄂東南礦集區(qū)地質(zhì)簡圖(據(jù)謝桂青等,2008修改)Fig.2 Geological sketch map of the Edong orefield(modified after Xie et al.,2008)

    龍角山礦區(qū)內(nèi)地層主要為志留系中統(tǒng)墳頭組、石炭系中統(tǒng)黃龍組和大埔組、二疊系下統(tǒng)棲霞組和茅口組、二疊系上統(tǒng)龍?zhí)督M、下窯組和大隆組、三疊系下統(tǒng)大冶組。地層整體呈北東東向,傾向北北西。礦區(qū)出露巖體有花崗閃長斑巖、閃長玢巖、煌斑巖,鉆孔揭露成礦巖體為花崗閃長斑巖(圖3),該巖體成巖年齡為(144±1)Ma,成礦作用發(fā)生于(144.7±2.9)Ma(丁麗雪等,2014)。礦體主要賦存在付家山向斜的北西翼和港溝山背斜的南東翼。

    圖3 龍角山礦床地質(zhì)圖(據(jù)紀(jì)云昊等,2019修改)Fig.3 Geologic map of the Longjiaoshan deposit(modified after Ji et al.,2019)

    龍角山鎢礦床中鎢礦體產(chǎn)于花崗閃長斑巖與和二疊系下統(tǒng)茅口組(圖4a)、棲霞組(圖4b)接觸帶,由多個(gè)不連續(xù)的透鏡狀礦體組成;鎢鉬礦體產(chǎn)于花崗閃長斑巖與二疊系下統(tǒng)棲霞組、茅口組接觸帶(圖4a),礦體總體走向北東,傾向南西,礦體沿走向和傾向延伸數(shù)百米。金屬礦物主要為白鎢礦、輝鉬礦、黃鐵礦,少量黃銅礦、磁鐵礦、方鉛礦、閃鋅礦;非金屬礦物主要為鈣鋁榴石、鈣鐵榴石、透輝石、硅灰石、陽起石、綠簾石、綠泥石、磷灰石、方解石、硬石膏和石英等礦物。

    根據(jù)礦物的結(jié)構(gòu)構(gòu)造和共生組合,可以將龍角山礦床的熱液期分為矽卡巖階段、退蝕變階段、石英-硫化物階段和碳酸鹽階段(圖5)。矽卡巖階段礦物組合為石榴子石-輝石-硅灰石,由巖體向大理巖延伸,矽卡巖種類依次為紅棕色石榴子石內(nèi)矽卡巖(圖6a、d)、黃綠色石榴子石-透輝石外矽卡巖(圖6b、e)、脈狀石榴子石-硅灰石外矽卡巖(圖6c、f)。退蝕變階段礦物組合為白鎢礦、陽起石、綠簾石、綠泥石、磷灰石、磁鐵礦,白鎢礦交代先形成的石榴子石等矽卡巖礦物(圖6g),呈團(tuán)塊狀包裹在石榴子石粒間或環(huán)帶(圖6h、i)。石英-硫化物階段礦物組合為石英、黃銅礦、黃鐵礦、白鎢礦(圖6j、k)、輝鉬礦、方鉛礦、閃鋅礦、磁黃鐵礦、硬石膏,白鎢礦呈半自形與黃鐵礦等硫化物共生(圖6l)。碳酸鹽階段礦物組合為方解石、白云石、石英等礦物。

    圖5 龍角山礦床熱液蝕變期次Fig.5 Hydrothermal alteration periods of the Longjiaosan deposit

    圖6 龍角山礦床手標(biāo)本和鏡下照片a.紅棕色石榴子石內(nèi)矽卡巖;b.黃綠色石榴子石-透輝石外矽卡巖;c.脈狀石榴子石-硅灰石外矽卡巖;d.紅棕色石榴子石內(nèi)矽卡巖(正交偏光);e.黃綠色石榴子石-透輝石外矽卡巖(單偏光);f.脈狀石榴子石-硅灰石外矽卡巖(單偏光);g.白鎢礦化矽卡巖;h.白鎢礦呈團(tuán)塊狀包裹在石榴子石粒間(正交偏光);i.白鎢礦呈團(tuán)塊狀包裹在石榴子石環(huán)帶(正交偏光);j.白鎢礦與硫化物共生;k.白鎢礦與硫化物共生(熒光);l.白鎢礦呈半自形與黃鐵礦共生(正交偏光)Fig.6 Hand specimen and photomicrographs of the Longjiaoshan deposit a.Red-brown garnet from endoskarn;b.Yellow-green garnet from exoskarn;c.Veined wollastonite-garnet from exoskarn;d.Photomicrographs of red-brown garnet(orthogonal);e.Photomicrographs of yellow-green garnet(orthogonal);f.Photomicrographs of veined wollastonite-garnet(polarized);g.Hand specimen of scheelite mineralization of skarn;h.Photomicrographs of scheelite metasomatic garnet grain(orthogonal);i.Photomicrographs of scheelite encased in garnet ring(orthogonal);j.Hand specimen of scheelite intergrowth with sulfide;k.Fluorescent image of scheelite intergrowth with sulfide(fluorescence);l.Photomicrographs of subhedral scheelite intergrowth with pyrite(orthogonal)

    1.3 采樣特征

    本次所測石榴子石和白鎢礦樣品分別采自龍角山礦床的不同鉆孔(圖4a、b)。

    圖4 龍角山礦床剖面圖(據(jù)湖北省第一地質(zhì)大隊(duì),2018修改)a.66線剖面圖;b.54線剖面圖Fig.4 Cross section map of the Longjiaoshan deposit(modified after No.1 Geological Party of Hubei Geological Bureau,2018)a.Cross section map of 66 Line;b.Cross section map of 54 Line

    紅棕色石榴子石內(nèi)矽卡巖樣品號為5402-551、5403-487、6601-878、6602-681,黃綠色石榴子石-輝石外矽卡巖樣品號為5401-607、5402-456、5403-339、6601-705、6602-569,脈狀石榴子石-硅灰石外矽卡巖樣號分別為5403-283、6602-531;退蝕變階段白鎢礦樣品號為5403-377、5403-398、6601-780、6602-603;石英-硫化物階段白鎢礦樣品號為5401-645、5402-491、5403-417、5403-450,其中,樣號前4位為鉆孔號,后3位為鉆孔取樣深度。

    2 測試分析方法

    樣品微觀形貌分析在合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院礦床成因與勘查技術(shù)研究中心(OEDC)礦物微區(qū)分析實(shí)驗(yàn)室熱場發(fā)射掃描電鏡Tescan MI‐RA3(設(shè)備配置布魯克60mm2 EDX能譜儀,陰極發(fā)光儀)上完成。

    主量元素電子探針(EPMA)分析在合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院完成儀器型號為JEOL JXA 8230。實(shí)驗(yàn)條件為:加速電壓15 kV,束斑尺寸3μm,探針電流20 nA。所有元素的信號采集時(shí)間均為15 s,背景時(shí)間均為5 s,修正方法ZAF,檢測限優(yōu)于100×10-6。

    微量元素和面掃描激光剝蝕電感耦合等離子質(zhì)譜(LA-ICP MS)分析在合肥工業(yè)大學(xué)礦物原位分析實(shí)驗(yàn)室完成。電感耦合等離子體質(zhì)譜儀由美國安捷倫公司制造,型號為Agilent 7900。激光剝蝕系統(tǒng)為萊伯泰科公司制造的Analyte HE。ArF準(zhǔn)分子激光發(fā)生器產(chǎn)生193 nm深紫外光束,經(jīng)勻化光路聚焦于礦物表面。激光束斑直徑為30μm,頻率為8 Hz,剝蝕時(shí)間40 s,以高純He氣為載氣,與Ar氣和少量N2氣混合后進(jìn)入質(zhì)譜儀。測試過程中以NIST SRM 610作為信號漂移矯正,以NIST610作為外標(biāo),用無內(nèi)標(biāo)法測定主量和微量元素含量。測試數(shù)據(jù)利用ICPMS Data Cal 9.9離線處理(Liu et al.,2010)。

    3 分析測試結(jié)果

    3.1 陰極發(fā)光和面掃描特征

    龍角山礦床退蝕變階段(Sch-1)和石英-硫化物階段(Sch-2)的白鎢礦陰極圖像有明顯的環(huán)帶特征,Sch-1呈現(xiàn)明顯的環(huán)帶特征,且顆粒邊部CL顏色變淺,據(jù)此細(xì)分為Sch-1a、Sch-1b和Sch-1c(圖7b、e);Sch-2邊部環(huán)帶特征顯著,據(jù)此分為Sch-2a和Sch-2b(圖7h)。

    龍角山礦床退蝕變階段(Sch-1)和石英-硫化物階段(Sch-2)的白鎢礦的面掃描特征顯示了不同的元素環(huán)帶,與陰極發(fā)光圖像對應(yīng),特別是Mo元素,面掃描圖像中Mo元素越富集,陰極發(fā)光圖像中對應(yīng)的區(qū)域灰度越深(圖7a~i)。

    圖7 龍角山礦床白鎢礦陰極發(fā)光和面掃描特征a.退蝕變階段白鎢礦(Sch-1)交代石榴子石(Grt()正交偏光);b.退蝕變階段白鎢礦(Sch-1)陰極發(fā)光圖像;c.退蝕變階段白鎢礦(Sch-1)Mo元素面掃描圖像;d.退蝕變階段白鎢礦(Sch-1)充填石榴子石粒間(正交偏光);e.退蝕變階段白鎢礦(Sch-1)陰極發(fā)光圖像;f.退蝕變階段白鎢礦(Sch-1)Mo元素面掃描圖像;g.自形-半自形石英-硫化物階段白鎢礦(Sch-2);h.石英-硫化物階段白鎢礦(Sch-2)陰極發(fā)光圖像;i.石英-硫化物階段白鎢礦(Sch-2)Mo元素面掃描圖像Fig.7 Cathodoluminescence and mapping images of scheelite from the Longjiaoshan deposit a.Retrograde scheelite(Sch-1)metasomasis garnet(Grt)(perpendicular polarized);b.Cathodoluminescence of retrograde scheelite(Sch-1);c.Map‐ping image of retrograde scheelite(Sch-1)of Mo element;d.Retrograde scheelite(Sch-1)metasomatic garnet(perpendicular polarized);e.Cathodo‐luminescence of retrograde scheelite(Sch-1);f.Mapping image of retrograde scheelite(Sch-1)of Mo element;g.Idiomorphic-semiidiomorphic quartz-sulfide stage scheelite(Sch-2);h.Cathodoluminescence of idiomorphic-semiidiomorphic quartz-sulfide stage scheelite(Sch-2);i.Mapping image of quartz-sulfide scheelite(Sch-2)of Mo element

    3.2 主量元素

    龍角山礦床石榴子石主量元素特征顯示其為鈣鋁-鈣鐵榴石系列。矽卡巖階段的紅棕色石榴子石(Grt-1)矽卡巖、黃綠色石榴子石(Grt-2)-輝石矽卡巖和脈狀石榴子石-硅灰石矽卡巖中石榴子石(Grt-3)的成分分別為Adr30.6-84.1Gro13.9-50.7Pyr1.5-30.3、Adr38.3-100Gro0.0-39.4Pyr0.0-22.4和Adr75.3-100.0Gro0.0-13.9Pyr0.0-12.7(表1,圖8),石榴子石的環(huán)帶特征與主量元素特征具有很好的耦合型(圖8)。前人的研究表明早階段形成的石榴子石富Al,晚階段形成的石榴子石富Fe(Einaudi et al.,1981;Nakano et al.,1989;Meinert,1997;Meinert et al.,2005),與本次實(shí)驗(yàn)結(jié)果吻合。

    圖8 龍角山礦床石榴子石三角圖Fig.8 Section of the ternary diagram for garnets from the Longjiaoshan deposit

    表1 龍角山礦床石榴子石電子探針分析結(jié)果Table1 Results of EPAManalysis of garnet from the Longjiaoshan deposit

    龍角山礦床白鎢礦主量元素特征顯示W(wǎng)O3和MoO3呈明顯的負(fù)相關(guān)關(guān)系(表2,圖9a),且退蝕變階段白鎢礦的w(MoO3)(平均值1.99%)明顯高于石英硫化物階段白鎢礦的w(MoO3)(平均值為0.99%)。

    表2 龍角山礦床白鎢礦電子探針分析結(jié)果Table 2 Results of EPAM analysis of scheelite from the Longjiaoshan deposit

    3.3 微量元素

    龍角山礦床石榴子石微量元素(表3)顯示,富集大離子親石元素(ΣLILE=17184.40×10-6),虧損高場強(qiáng)元素(ΣHFSE=12.13×10-6)。上述3類石榴子石的稀土元素呈現(xiàn)明顯差異性,Grt-1稀土元素含量最高(ΣREEs=17.60×10-6),Grt-2含 量 較 高(ΣREEs=13.21×10-6),Grt-3稀土元素含量最低(ΣREEs=9.41×10-6)。上述3類石榴子石的稀土元素配分模型分布呈現(xiàn)中稀土元素富集、銪正異常的特征(圖10a~c)。在Grt-2中,金屬元素w(W)(平均值分別22.6×10-6)明顯高于其他類型的石榴子石。

    表3 龍角山礦床石榴子石微量元素分析結(jié)果Table 2 Trace elements of the garnet from the Longjiaoshan deposit

    龍角山礦床白鎢礦微量元素虧損大離子親石元素,如金屬元素w(Rb、Ba、U、Th)低于1×10-6,但w(Sr)(213.26×10-6)明顯高于其他大離子親石元素,且虧損高場強(qiáng)元素,如Zr、Hf、Ta等低于1×10-6,但w(Nb)(2.60×10-6)明顯高于其他高場強(qiáng)元素(表4)。上述2類白鎢礦的稀土元素配分模型呈現(xiàn)過渡特征,從退蝕變階段到石英硫化物階段白鎢礦的LREE含量,逐漸降低且正銪異常逐漸加強(qiáng)(2類白鎢礦的δEu平均值分別為0.526和2.423,圖10d~f)。

    4 討論

    龍角山礦床中石榴子石和白鎢礦含有多種微量元素,如REEs、LILEs(Sr)、HFSEs(Nb)、Mo、Sn、U等,可為巖漿熱液體系來源和演化提供證據(jù),如氧逸度、pH值、水巖反應(yīng)程度等(劉善寶等,2007;Song et al.,2014;Guo et al.,2016;Park et al.,2017;Xiao et al.,2018),因此,本文從上述幾個(gè)方面分別討論龍角山礦床的成礦流體特征。

    4.1 成礦流體特征

    4.1.1 氧逸度

    U是一種對氧化還原敏感的元素,因?yàn)閁4+比U6+更容易替代石榴子石中的鈣元素,降低f(O2)可以降低U在流體中的溶解度,從而增加石榴石中的w(U)(Smith et al.,2004;Gaspar et al.,2008)。由于從巖漿中分離出來的流體脈沖不同,形成的石榴子石的生成量不同(Shu et al.,2017),熱液流體中初始w(U)相似,石榴子石的w(U)受流體氧化還原狀態(tài)控制。因此,石榴子石中低w(U)反映了石榴子石形成于氧化環(huán)境,而高w(U)則反映了石榴子石形成于還原環(huán)境。龍角山Grt-1的w(U)最高,說明它是在相對還原的環(huán)境中形成的(圖9b)。Grt-3的w(U)最低,表明其形成于氧化環(huán)境,這與石榴子石的鈣鐵榴石的成分變化吻合,表明在矽卡巖階段成礦流體氧逸度逐漸升高。

    圖9 龍角山礦床白鎢礦和石榴子石主、微量元素圖解a.白鎢礦WO3-MoO3圖解;b.石榴子石ΣREE-U圖解Fig.9 Diagrams of major and trace elements of scheelite and garnet from the Longjiaoshan deposit a.Diagram of WO3-MoO3 of scheelite;b.Diagram ofΣREE-U of garnet

    Mo也是一種對氧化還原敏感的元素,在氧化條件下以Mo6+的形式遷移并替代W6+進(jìn)入白鎢礦。當(dāng)f(O2)降低時(shí),Mo6+還原為Mo4+,沉淀為輝鉬礦(MoS2)(Rempel et al.,2009;Song et al.,2014)。因次,顯示了成礦過程復(fù)雜的脈沖式的氧逸度變化特征。

    綜上所述,龍角山礦床矽卡巖階段到退蝕變階段再到石英-硫化物階段,成礦流體的氧逸度經(jīng)歷了2次先增加后降低的過程,顯示了成礦過程復(fù)雜的氧逸度變化特征。

    4.1.2 pH值

    除氧逸度外,成礦流體pH值的變化也控制了研究區(qū)成礦作用的形成,并且顯著影響矽卡巖熱液中稀土元素的分餾(Bau,1991)。通常來說,在中性條件下,石榴子石的稀土元素模式為富集重稀土元素和虧損輕稀土元素,且Eu呈負(fù)異?;驘o異常;在中等酸性pH值條件下,稀土元素模式更多地受Cl-的控制,Cl-的存在可以增強(qiáng)除REE3+外的可溶Eu2+離子的穩(wěn)定性,并導(dǎo)致明顯的Eu正異常,且富集輕稀土元素、虧損重稀土元素(Bau,1991;Gaspar et al.,2008;Zhang et al.,2017)。在龍角山礦床中,Grt-1、Grt-2和Grt-3均表現(xiàn)出輕稀土元素富集、重稀土元素虧損、銪正異常,且稀土元素含量逐漸減少的特征(圖10),說明在矽卡巖階段,成礦流體熱液pH值較穩(wěn)定,為酸性條件。低pH值水溶液可以攜帶高濃度的鎢(Wood et al.,2000;Wang et al.,2019),即pH值的任何增加都可能導(dǎo)致流體中鎢的飽和,從而導(dǎo)致白鎢礦的沉淀。pH值的改變的成礦系統(tǒng)在退蝕變階段發(fā)生,此階段成礦流體與碳酸鹽巖類圍巖反應(yīng)程度達(dá)到峰值,這一過程將增加成礦流體pH值(Legros et al.,2020)。

    圖10 龍角山礦床白鎢礦和石榴子石稀土元素配分圖解a.Grt-1稀土元素配分圖解;b.Grt-2稀土元素配分圖解;c.Grt-3稀土元素配分圖解;d.Sch-1稀土元素配分圖解;e.Sch-2稀土元素配分圖解Fig.10 Normalized REE pattern of scheelite and garnet from the Longjiaoshan deposit a.Normalized REE pattern of Grt-1;b.Normalized REE pattern of Grt-2;c.Normalized REE pattern of Grt-3;d.Normalized REE pattern of Sch-1;e.Normalized REE pattern of Sch-2

    白鎢礦的正銪異常特征為成礦流體與富鈣質(zhì)巖石或礦物反應(yīng)提供了證據(jù),這是由于富鈣質(zhì)巖石(如灰?guī)r)或礦物(如斜長石)中Eu以+2價(jià)離子形式存在,且Eu2+較Eu3+更易進(jìn)入白鎢礦晶格(Shannon,1976;Cottrant,1981;Raimbault et al.,1993;Ghaderi et al.,1999;Brugger et al.,2008),因此,正銪異常顯示了熱液流體與圍巖發(fā)生充分反應(yīng)(Sun et al.,2017;Wang et al.,2017;Zhang et al.,2018;Wu et al.,2019)。龍角山礦床2類白鎢礦顯示銪負(fù)異常減弱銪正異常增加的特征,這與上文前人提出的成礦階段熱液pH增加相吻合。綜上所述,龍角山礦床矽卡巖階段到退蝕變階段再到石英-硫化物階段,成礦流體的pH值逐漸增加,而非在退蝕變階段(成礦階段)達(dá)到峰值,因此,pH值不是影響白鎢礦沉淀的主要因素。

    4.1.3 水巖反應(yīng)程度

    Nb/La 0.114 0.180 0.189 0.087 0.234 0.060 0.161 0.064 0.060 1.882 0.181 0.351 0.138 0.152 0.091 Y 3.510 0.894 1.382 0.397 0.077 0.595 1.348 0.183 0.125 0.645 0.713 0.487 1.342 0.318 0.564 Rb/Sr 0 0 0 0 0 0 0 0 0 0.001 0 0.001 0.003 0.002 0 Lu 0.004 0.008 0 0 0.005 0.002 0.002 0.001 0 0.005 0 0 0 0.001 0.002 Zr/Hf ----8.160 2.715 3.342-2.421 14.098 3.463 13.309-0-Yb 0.009 0.036 0.019 0.028 0.053 0.017 0.008 0.008 0 0.019 0.016 0.012 0 0.072 0.028 δEu 0.800 0.495 0.175 0.499 0.730 0.701 0.303 0.361 0.551 3.944 1.015 1.892 4.230 2.333 2.360 Tm 0.005 0.015 0 0.001 0.011 0.005 0.001 0 0.001 0.006 0.005 0.003 0 0.013 0.007 Nb+Ta 2.807 5.546 2.076 1.093 2.120 0.987 3.875 1.233 1.069 3.321 2.363 3.431 2.260 2.634 1.655 0.063 Er 0.117 0.042 0.061 0.083 0.058 0.032 0.029 0.013 0.101 0.098 0.005 0 0.143 0.037 U 0.025 0.007 0.010 0.029 0.017 0.652 0.018 0.018 0.034 0.093 0.026 0.053 0.076 0.004 0.002 0.053 Ho 0.013 0.046 0.026 0.015 0.034 0.017 0.010 0.015 0.064 0.068 0.004 0 0.073 0.020 Th 0.040 0.028 0.054 0.011 0.006 0.012 0.003 0.009 0.015 0.003 0.005 0.024 0.019 0.018 0.010果結(jié)析Dy 0.169 0.355 0.055 0.230 0.092 0.263 0.071 0.081 0.057 0.377 0.312 0.043 0 0.506 0.201 W分620670 620800 621128 602261 611403 602539 630887 632591 622737 616078 612467 622659 634348 633808 634870素元w(B)/10-6量Tb 0.048 0.080 0.020 0.056 0.026 0.057 0.019 0.018 0.005 0.102 0.091 0 0 0.072 0.039 Ta 0.048 0.102 0.052 0.021 0.045 0.009 0.013 0.012 0.013 0.096 0.077 0.090 0.016 0.011 0.012微礦鎢白Gd 0.373 0.965 0.329 0.674 0.133 0.654 0.179 0.262 0.174 1.175 0.917 0.132 0.107 0.860 0.398 Hf 0 0 0 0 0.003 0.005 0.003 0 0.006 0.010 0.005 0.005 0 0.002 0床礦山w(B)/10-6 Cs角Eu 0.129 0.212 0.025 0.137 0.053 0.176 0.022 0.026 0.024 1.986 0.414 0.143 0.162 0.836 0.317 0.015 0 0 0 0 0 0.021 0 0.024 0.055 0.085 0.079 0.021 0.065 0龍4表TraceelementsofthescheelitefromtheLongjiaoshandeposit Sm 0.572 1.543 0.523 0.952 0.282 0.836 0.243 0.180 0.098 1.788 1.474 0.300 0.123 1.254 0.412 Mo 9590.008 9543.613 9445.386 19678.654 15948.468 19747.815 4417.656 4214.091 7692.262 12828.729 15254.013 7813.445 3667.033 3897.073 3610.831 Table4Nd 7.917 18.629 4.885 9.308 2.463 11.166 2.282 1.987 1.296 6.339 14.872 4.761 13.178 4.612 2.264 Nb 2.759 5.444 2.024 1.072 2.075 0.978 3.861 1.220 1.056 3.225 2.286 3.340 2.244 2.623 1.643 Pr 3.387 6.862 1.250 3.645 0.764 4.251 1.155 0.906 0.784 1.250 5.177 1.659 4.862 0.954 0.703 Zr 0 0.012 0.004 0.037 0.028 0.013 0.012 0.011 0.016 0.143 0.018 0.072 0.005 0 0.008 Ce 43.611 78.350 7.317 37.142 4.470 46.830 26.004 18.027 16.754 5.816 45.216 15.972 43.457 15.078 14.761 Sr 205.554 186.397 187.094 246.226 192.846 240.857 137.084 197.078 177.385 298.483 458.288 250.737 158.760 164.614 147.001 La 24.253 30.160 10.700 12.362 8.863 16.431 23.922 19.207 17.489 1.714 12.621 9.514 16.314 17.243 18.088 Rb 0.059 0.004 0 0.079 0 0.068 0.047 0 0.007 0.241 0.206 0.261 0.460 0.405 0.025號號樣5403-398-a1 5403-398-a2 5403-398-a3 5403-398-b1 5403-398-b10 5403-398-b2 5403-398-c1 5403-398-c2 5403-398-c3 5401-645-a1 5401-645-a10 5401-645-a2 5401-645-b1 5401-645-b2 5401-645-b3樣5403-398-a1 5403-398-a2 5403-398-a3 5403-398-b1 5403-398-b10 5403-398-b2 5403-398-c1 5403-398-c2 5403-398-c3 5401-645-a1 5401-645-a10 5401-645-a2 5401-645-b1 5401-645-b2 5401-645-b3 1。為位單值品品樣Sch-1a Sch-1a Sch-1a Sch-1b Sch-1b Sch-1b Sch-1c Sch-1c Sch-1c Sch-2a Sch-2a Sch-2a Sch-2b Sch-2b Sch-2b樣Sch-1a Sch-1a Sch-1a Sch-1b Sch-1b Sch-1b Sch-1c Sch-1c Sch-1c Sch-2a Sch-2a Sch-2a Sch-2b Sch-2b Sch-2b:比注

    矽卡巖的形成是一個(gè)動(dòng)態(tài)過程,包括不同階段的熱液活動(dòng)和巖漿熱液流體的持續(xù)演化(Meinert,1997;Meinert et al.,2005)。Y作為一種特殊的稀土元素,在平衡體系中,Y與REE存在正相關(guān)關(guān)系;在非平衡體系中,Y與REEs無正相關(guān)性且富集輕稀土元素(Gaspar et al.,2008;Yardley et al.,1991)。本文測的Grt-1稀土元素總量與Y存在正相關(guān)關(guān)系,而Grt-2、Grt-3稀土元素總量與Y不存在正相關(guān)關(guān)系(圖11a),反映了隨著石榴子石的結(jié)晶由平衡狀態(tài)向非平衡狀態(tài)轉(zhuǎn)變,即水巖反應(yīng)程度逐漸增加,由封閉的熱液系統(tǒng)逐步向開放的熱液系統(tǒng)過渡(Park et al.,2017),矽卡巖從反應(yīng)矽卡巖(擴(kuò)散)演化為交代矽卡巖(滲透)(Meinert et al.,2005),而開放體系石榴子石以富Fe為特征(Gaspar et al.,2008),這與本文所觀察到的Grt-3鈣鐵榴石成分高(Adr75.3-100.0Gro0.0-13.9Pyr0.0-12.7,圖8)相吻合。這一過程增加了矽卡巖蝕變灰?guī)r的孔隙度和滲透率,允許進(jìn)一步的熱液流動(dòng),從而增加了體系的水巖比,且對應(yīng)形成的富Fe石榴子石稀土元素總含量逐漸減少,這一現(xiàn)象在龍角山礦床石榴子石的稀土元素特征中得到印證(表3),且Sch-1和Sch-2的w(Sr)呈現(xiàn)增加的趨勢(圖11b),同時(shí),由于熱液礦物生長時(shí)物理和化學(xué)條件的波動(dòng),晶體內(nèi)部會產(chǎn)生微小的缺陷,而CL或BSE圖像可以精準(zhǔn)的反映礦物的微觀形貌學(xué)特征(Rusk et al.,2002;2008;Putnis,2009;Han et al.,2020),單個(gè)白鎢礦顆粒中不同灰度的區(qū)域代表不同世代白鎢礦在不同物理和化學(xué)條件下的流體平衡狀態(tài)。在Sch-1a和Sch-2a顆粒的周邊均表現(xiàn)不規(guī)則生長特征(即Sch-1b、Sch-1c和Sch-2b)(圖7a~i),且主成礦階段Sch-1的不規(guī)則生長環(huán)帶比Sch-2更發(fā)育,這也反映了主成礦階段(退蝕變階段)水巖反應(yīng)程度達(dá)到峰值。

    圖11 龍角山礦床石榴子石和白鎢礦圖解a.石榴子石ΣREE和Y元素圖解;b.白鎢礦Sr元素圖解Fig.11 Diagrams of scheelite and garnet from the Longjiaoshan deposit a.Diagram ofΣREE and Y of garnet;b.Diagram of Sr of scheelite

    4.2 龍角山礦床鎢的富集沉淀機(jī)制

    前人研究表明,鎢在熱液中的富集和遷移機(jī)制受控于流體物理化學(xué)條件(溫度、壓力、氧逸度、pH值)的改變(Bai et al.,1999;Wood et al.,2000;Zajacz et al.,2008);控制鎢的沉淀機(jī)制主要有:溫壓降低(Ni et al.,2015;Chen et al.,2018)、水巖反應(yīng)(Lecum‐berri et al.,2017)、流體混合(Wei et al.,2012;Legros et al.,2019;Pan et al.,2019)、流體不混溶(Korges et al.,2018;王國光等,2020),如大湖塘鎢銅多金屬礦床成礦流體發(fā)生強(qiáng)烈的水巖反應(yīng)導(dǎo)致礦質(zhì)沉淀(Sun et al.,2017)。柿竹園鎢錫多金屬礦床大氣降水參與成礦體系是白鎢礦沉淀的主要機(jī)制(祝新友等,2015)。南泥湖鎢鉬礦床經(jīng)歷了流體混合和流體不混溶作用(Yang et al.,2012;蔣少涌等,2020)。淘錫坑鎢礦和巖前鎢礦成礦流體的不混溶作用是礦質(zhì)沉淀的主要原因(魯麟等,2018;劉暢等,2018)。

    在上述4種鎢的沉淀機(jī)制中,降溫減壓過程可能不是白鎢礦礦床形成的有效因素(Foster,1977;Wood et al.,2000),因?yàn)?00~500℃范圍的白鎢礦溶解度隨著溫度降低而增加。由前文(4.1.1)可知,龍角山礦床成礦流體氧逸度變化復(fù)雜,存在多次升高降低的過程。而流體混合很可能伴著的f(O2)和pH值增加,溫度降低(Linnen et al.,1994;Singoyi et al.,2001;Wei et al.,2012),因此,龍角山白鎢礦氧逸度的多次升高很可能反映了有新的流體混入,前人研究也表明龍角山礦床成礦流體存在大氣水和有機(jī)質(zhì)混入(Lei et al.,2018),因此,流體混合是龍角山鎢礦床礦質(zhì)沉淀的重要原因之一。

    由前文(4.1.2)可知,pH值不是影響龍角山礦床白鎢礦沉淀的主要因素。而流體不混溶則導(dǎo)致壓力降低和pH值升高(Lu et al.,2003;Korges et al.,2017;Orhan,2017;Soloviev et al.,2017)。因此,龍角山礦床的流體不混溶作用可忽略。由前文(4.1.3)可知,水巖反應(yīng)在龍角山礦床主成礦階段達(dá)到峰值,有利于形成白鎢礦。水巖反應(yīng)主要伴隨著熱液非極性揮發(fā)分的加入、Ca離子富集和pH值增加(Gibert et al.,1992;O’Reilly et al.,1997),通常被認(rèn)為是形成鎢礦床的關(guān)鍵機(jī)制(Lecumberri et al.,2017)。

    綜上所述,流體混合和水巖反應(yīng)是控制龍角山礦床鎢沉淀的主要機(jī)制,二者協(xié)同控制白鎢礦沉淀成礦。流體混合作用被認(rèn)為是形成具有異常高品位鎢礦床的主要沉淀機(jī)制(Wei et al.,2012;Korges et al.,2017)。水巖反應(yīng)也是形成大型、超大型鎢多金屬礦床的重要過程(如大湖塘鎢銅多金屬礦床,Peng et al.,2018),通過本次工作可知,龍角山鎢礦床的礦質(zhì)沉淀機(jī)制兼具流體混合和水巖作用,該礦區(qū)具有成大礦、富礦的潛質(zhì)。

    5 結(jié)論

    (1)龍角山礦床矽卡巖階段的紅棕色石榴子石(Grt-1)矽卡巖、黃綠色石榴子石(Grt-2)-輝石矽卡巖和脈狀石榴子石-硅灰石矽卡巖中石榴子石(Grt-3)的成分分別為Adr30.6-84.1Gro13.9-50.7Pyr1.5-30.3、Adr38.3-100Gro0.0-39.4Pyr0.0-22.4和Adr75.3-100.0Gro0.0-13.9Pyr0.0-12.7,且富集大離子親石元素,虧損高場強(qiáng)元素,具有富集輕稀土元素、虧損重稀土元素、銪正異常的特征。

    (2)龍角山礦床矽卡巖階段石榴子石U含量逐漸降低,且與鈣鐵榴石的成分變化吻合,表明在矽卡巖階段成礦流體氧逸度逐漸升高;退蝕變階段到石英-硫化物階段白鎢礦w(Mo)先升高、后降低,對應(yīng)的氧逸度先增加后降低,顯示了成礦過程復(fù)雜的脈沖式的氧逸度變化特征。

    (3)龍角山礦床退蝕變階段和石英-硫化物階段對應(yīng)形成的白鎢礦顯示銪負(fù)異常減弱、銪正異常增加的變化特征,表明成礦流體的pH值逐漸增加,而非在退蝕變階段(主成礦階段)達(dá)到峰值,因此,pH值不是影響白鎢礦沉淀的主要因素。

    (4)龍角山礦床矽卡巖階段Grt-1稀土元素總量與Y存在正相關(guān)關(guān)系,Grt-2和Grt-3稀土元素總量與Y不存在正相關(guān)關(guān)系,且Grt-3富Fe,表明隨著石榴子石的結(jié)晶由熱液平衡狀態(tài)向非平衡條件轉(zhuǎn)變,即水巖反應(yīng)程度逐漸增加;退蝕變階段(主成礦階段)Sch-1的不規(guī)則生長環(huán)帶比石英硫化物階段Sch-2更發(fā)育,表明主成礦階段水巖反應(yīng)程度達(dá)到峰值,水巖反應(yīng)是控制該礦床礦質(zhì)沉淀的重要因素。

    (5)流體混合和水巖反應(yīng)是控制龍角山礦床鎢沉淀的主要機(jī)制,二者協(xié)同控制白鎢礦沉淀成礦,龍角山礦區(qū)具有成大礦、富礦的潛質(zhì)。

    致 謝野外地質(zhì)工作期間得到了湖北省地質(zhì)局第一地質(zhì)大隊(duì)楊偉衛(wèi)院長等多位工程師的大力支持和幫助。合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院汪方躍副教授在電鏡分析過程中提供了熱情的指導(dǎo)。審稿專家為本文提出了寶貴的修改意見。在此一并對上述人員志以謝忱!

    猜你喜歡
    子石白鎢礦矽卡巖
    激電聯(lián)合剖面在判斷矽卡巖型礦床礦體產(chǎn)狀中的應(yīng)用
    Fe-BHA用于白鎢礦的浮選分離
    礦冶(2020年4期)2020-08-21 08:16:24
    高鈣白鎢礦分解工業(yè)試驗(yàn)研究
    硫磷混酸分解白鎢礦過程中伴生稀土的行為
    Geochemical and geochronological studies of the Aketas granite from Fuyun County,Xinjiang:the implications of the petrogenesis and tectonic setting
    Acta Geochimica(2015年4期)2015-10-25 02:03:16
    打括號的警察之墓
    牡丹(2014年10期)2014-10-17 09:31:59
    石榴子石Lu-Hf年代學(xué)研究進(jìn)展
    WorldView—Ⅱ在青海東昆侖求勉雷克塔格地區(qū)矽卡巖型礦床找礦中的應(yīng)用
    河北遙感(2014年3期)2014-07-10 13:16:47
    選礦尾礦中的石榴子石可選性試驗(yàn)研究
    具有類白鎢礦結(jié)構(gòu)的KGd(MoO4)2的晶體結(jié)構(gòu)和能帶結(jié)構(gòu)
    欧美日韩综合久久久久久| 亚洲欧美一区二区三区黑人 | 日韩一本色道免费dvd| 美女中出高潮动态图| 久久精品国产a三级三级三级| 亚洲人与动物交配视频| 在线 av 中文字幕| 婷婷色综合www| 少妇人妻精品综合一区二区| av在线老鸭窝| 日韩欧美一区视频在线观看| 日韩一本色道免费dvd| 久久人人爽人人爽人人片va| 9热在线视频观看99| tube8黄色片| 成人免费观看视频高清| 国产一区二区在线观看av| 看免费av毛片| 亚洲精品国产av蜜桃| 亚洲国产精品999| 天天躁夜夜躁狠狠久久av| av又黄又爽大尺度在线免费看| 18禁在线无遮挡免费观看视频| 另类精品久久| 亚洲欧美成人精品一区二区| 成人毛片60女人毛片免费| 成人国语在线视频| 99久国产av精品国产电影| 国产乱人偷精品视频| 欧美+日韩+精品| 最新的欧美精品一区二区| 欧美少妇被猛烈插入视频| 中文乱码字字幕精品一区二区三区| 国产亚洲av片在线观看秒播厂| 日韩熟女老妇一区二区性免费视频| 日本与韩国留学比较| 多毛熟女@视频| www日本在线高清视频| 精品亚洲成国产av| 精品国产一区二区三区久久久樱花| 久久精品久久久久久噜噜老黄| 视频中文字幕在线观看| 欧美日韩亚洲高清精品| 夫妻性生交免费视频一级片| 免费播放大片免费观看视频在线观看| 一二三四中文在线观看免费高清| 熟妇人妻不卡中文字幕| 国产av精品麻豆| 久久这里只有精品19| 国产成人免费无遮挡视频| 如何舔出高潮| 久久女婷五月综合色啪小说| 99热网站在线观看| 欧美国产精品va在线观看不卡| 丁香六月天网| 22中文网久久字幕| 日韩制服骚丝袜av| 免费观看性生交大片5| 日韩一区二区三区影片| 亚洲国产日韩一区二区| 精品久久久久久电影网| 中文字幕精品免费在线观看视频 | 制服丝袜香蕉在线| 久久久精品免费免费高清| 成年人午夜在线观看视频| www日本在线高清视频| 亚洲成人一二三区av| 在线观看www视频免费| 青春草国产在线视频| 天天躁夜夜躁狠狠躁躁| tube8黄色片| a级片在线免费高清观看视频| 国产在线免费精品| 欧美精品人与动牲交sv欧美| 制服人妻中文乱码| 久久精品久久久久久久性| 韩国av在线不卡| 寂寞人妻少妇视频99o| 人人妻人人爽人人添夜夜欢视频| 国产精品国产av在线观看| 国产成人精品一,二区| 久久精品国产a三级三级三级| 亚洲第一av免费看| 考比视频在线观看| www.av在线官网国产| 久久99蜜桃精品久久| xxxhd国产人妻xxx| 国产精品99久久99久久久不卡 | 国产精品久久久久久精品古装| 高清欧美精品videossex| 国产精品久久久久久av不卡| 99久久综合免费| 国产精品久久久久久精品电影小说| 久久精品夜色国产| 五月天丁香电影| 最近2019中文字幕mv第一页| 欧美成人精品欧美一级黄| 免费少妇av软件| 国产不卡av网站在线观看| 欧美最新免费一区二区三区| 一级片免费观看大全| 观看av在线不卡| 日韩av免费高清视频| 亚洲成人手机| 一边摸一边做爽爽视频免费| 乱码一卡2卡4卡精品| 精品少妇黑人巨大在线播放| 在线看a的网站| 国产黄色免费在线视频| 99香蕉大伊视频| www.色视频.com| 99热国产这里只有精品6| 美女国产视频在线观看| 国产一区亚洲一区在线观看| 亚洲精品久久午夜乱码| 9191精品国产免费久久| 黑人欧美特级aaaaaa片| 亚洲精品久久午夜乱码| 免费观看性生交大片5| 精品久久国产蜜桃| 欧美最新免费一区二区三区| 欧美最新免费一区二区三区| 成人手机av| 日韩大片免费观看网站| 午夜激情av网站| kizo精华| 成人国产av品久久久| 最近中文字幕2019免费版| 欧美人与性动交α欧美精品济南到 | 亚洲激情五月婷婷啪啪| 免费大片18禁| 丝袜喷水一区| 22中文网久久字幕| 亚洲高清免费不卡视频| 青春草视频在线免费观看| 最新的欧美精品一区二区| 免费观看性生交大片5| 国产乱人偷精品视频| 在线看a的网站| 制服人妻中文乱码| 一个人免费看片子| av在线观看视频网站免费| 日本爱情动作片www.在线观看| 久久久久久久大尺度免费视频| 成人亚洲精品一区在线观看| 伦理电影免费视频| 乱人伦中国视频| 成年美女黄网站色视频大全免费| 97在线视频观看| 欧美丝袜亚洲另类| 日本黄大片高清| 91久久精品国产一区二区三区| 超碰97精品在线观看| 国产精品一区二区在线观看99| 免费不卡的大黄色大毛片视频在线观看| 五月天丁香电影| 人妻一区二区av| 亚洲国产av影院在线观看| www日本在线高清视频| 欧美老熟妇乱子伦牲交| 99热国产这里只有精品6| 我要看黄色一级片免费的| 欧美变态另类bdsm刘玥| 国产一区二区在线观看av| av女优亚洲男人天堂| 国产精品人妻久久久久久| 少妇被粗大猛烈的视频| 51国产日韩欧美| 中国美白少妇内射xxxbb| 69精品国产乱码久久久| av女优亚洲男人天堂| 欧美激情国产日韩精品一区| 精品午夜福利在线看| 天堂俺去俺来也www色官网| 欧美最新免费一区二区三区| 国产精品女同一区二区软件| 精品国产一区二区三区四区第35| 日产精品乱码卡一卡2卡三| 热re99久久国产66热| 成年女人在线观看亚洲视频| 国产极品天堂在线| 男女边吃奶边做爰视频| 久久久久人妻精品一区果冻| 国产成人av激情在线播放| 2018国产大陆天天弄谢| 中文字幕av电影在线播放| 国产亚洲精品久久久com| 51国产日韩欧美| h视频一区二区三区| 亚洲欧美精品自产自拍| 精品久久久久久电影网| 丝袜美足系列| 熟妇人妻不卡中文字幕| freevideosex欧美| 一本大道久久a久久精品| 免费看不卡的av| 又大又黄又爽视频免费| 人妻少妇偷人精品九色| 免费av不卡在线播放| 日本猛色少妇xxxxx猛交久久| 久久久久久人人人人人| 久久免费观看电影| 国产精品久久久久成人av| 熟女av电影| 日韩制服骚丝袜av| 美女国产高潮福利片在线看| 精品亚洲成国产av| 国产欧美日韩综合在线一区二区| 欧美日韩av久久| 大香蕉久久网| 搡老乐熟女国产| 成年av动漫网址| 日韩av不卡免费在线播放| 免费观看性生交大片5| 巨乳人妻的诱惑在线观看| 男人操女人黄网站| 亚洲精品第二区| 久久久久久久久久久免费av| 激情五月婷婷亚洲| 91成人精品电影| 国产黄频视频在线观看| 欧美 亚洲 国产 日韩一| 亚洲国产精品999| 亚洲一级一片aⅴ在线观看| 热re99久久精品国产66热6| 免费女性裸体啪啪无遮挡网站| 黑人欧美特级aaaaaa片| 午夜影院在线不卡| 亚洲精品视频女| 国产精品三级大全| 国产亚洲精品第一综合不卡 | 青春草亚洲视频在线观看| 日日啪夜夜爽| 哪个播放器可以免费观看大片| 免费日韩欧美在线观看| 波野结衣二区三区在线| 观看美女的网站| 精品国产一区二区三区四区第35| av片东京热男人的天堂| 欧美日韩亚洲高清精品| 一区二区三区精品91| 亚洲综合精品二区| 亚洲性久久影院| 亚洲精品乱久久久久久| 男人舔女人的私密视频| 男人操女人黄网站| 一级黄片播放器| 99精国产麻豆久久婷婷| 亚洲,一卡二卡三卡| 精品视频人人做人人爽| 日韩三级伦理在线观看| 亚洲图色成人| 亚洲色图 男人天堂 中文字幕 | 国产精品一区二区在线观看99| 国产视频首页在线观看| 国产xxxxx性猛交| 国产 一区精品| 下体分泌物呈黄色| 亚洲av电影在线进入| 久久午夜综合久久蜜桃| 两个人免费观看高清视频| 草草在线视频免费看| 日本与韩国留学比较| 女性被躁到高潮视频| 久久精品人人爽人人爽视色| 亚洲,一卡二卡三卡| 热99久久久久精品小说推荐| 有码 亚洲区| 国产深夜福利视频在线观看| 一区二区三区精品91| 国产片特级美女逼逼视频| 日韩熟女老妇一区二区性免费视频| 欧美精品国产亚洲| 精品亚洲成国产av| 精品一区二区免费观看| 纯流量卡能插随身wifi吗| 亚洲精品乱码久久久久久按摩| 99九九在线精品视频| 国产伦理片在线播放av一区| 国产视频首页在线观看| 80岁老熟妇乱子伦牲交| 中文精品一卡2卡3卡4更新| 女人被躁到高潮嗷嗷叫费观| 午夜av观看不卡| 曰老女人黄片| 天堂中文最新版在线下载| 免费观看无遮挡的男女| 欧美激情国产日韩精品一区| 欧美性感艳星| 日韩成人伦理影院| 久久影院123| 国产又爽黄色视频| 成人国语在线视频| 乱码一卡2卡4卡精品| 少妇被粗大猛烈的视频| 欧美 日韩 精品 国产| 91午夜精品亚洲一区二区三区| 国语对白做爰xxxⅹ性视频网站| 看免费成人av毛片| 91在线精品国自产拍蜜月| 日韩在线高清观看一区二区三区| 夫妻午夜视频| 国产淫语在线视频| 少妇的逼好多水| 男女午夜视频在线观看 | 亚洲av在线观看美女高潮| 久久99蜜桃精品久久| 久久久久网色| 精品国产露脸久久av麻豆| av国产精品久久久久影院| 日韩成人伦理影院| 久久热在线av| 亚洲欧美成人精品一区二区| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 少妇被粗大的猛进出69影院 | 中文字幕制服av| 汤姆久久久久久久影院中文字幕| 一二三四在线观看免费中文在 | 9191精品国产免费久久| 成人亚洲欧美一区二区av| 免费观看无遮挡的男女| 五月开心婷婷网| 精品少妇久久久久久888优播| 欧美老熟妇乱子伦牲交| 成人漫画全彩无遮挡| h视频一区二区三区| 黑人高潮一二区| 韩国精品一区二区三区 | 久久国产精品男人的天堂亚洲 | 18+在线观看网站| 国产男女内射视频| 免费人成在线观看视频色| 桃花免费在线播放| 国产日韩欧美亚洲二区| 最新的欧美精品一区二区| 亚洲欧美日韩卡通动漫| 妹子高潮喷水视频| 九草在线视频观看| 汤姆久久久久久久影院中文字幕| 亚洲 欧美一区二区三区| 亚洲av免费高清在线观看| 欧美精品高潮呻吟av久久| 久久精品国产自在天天线| 黄网站色视频无遮挡免费观看| 国产永久视频网站| 九色亚洲精品在线播放| 天天操日日干夜夜撸| 国产精品熟女久久久久浪| 色网站视频免费| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 热re99久久国产66热| 伦理电影大哥的女人| 亚洲,一卡二卡三卡| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 欧美+日韩+精品| 久久人人爽av亚洲精品天堂| 久久久久久久久久成人| 校园人妻丝袜中文字幕| 九草在线视频观看| 亚洲高清免费不卡视频| 香蕉丝袜av| 日韩免费高清中文字幕av| 老司机影院毛片| 久久韩国三级中文字幕| 在线观看免费视频网站a站| av播播在线观看一区| 亚洲欧洲精品一区二区精品久久久 | 欧美人与性动交α欧美软件 | 中国国产av一级| 在线观看人妻少妇| 中文欧美无线码| 高清欧美精品videossex| 捣出白浆h1v1| 少妇的丰满在线观看| 国产精品久久久久久久久免| 亚洲精品aⅴ在线观看| 免费观看无遮挡的男女| 久久精品aⅴ一区二区三区四区 | 国产女主播在线喷水免费视频网站| 日日摸夜夜添夜夜爱| 日韩电影二区| 国产精品秋霞免费鲁丝片| 亚洲精品成人av观看孕妇| a级毛片在线看网站| 少妇 在线观看| freevideosex欧美| 999精品在线视频| 少妇被粗大猛烈的视频| 97在线人人人人妻| 少妇的丰满在线观看| 夜夜爽夜夜爽视频| 91久久精品国产一区二区三区| 久久久欧美国产精品| 亚洲av电影在线观看一区二区三区| 97超碰精品成人国产| 中国国产av一级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 综合色丁香网| 国产精品.久久久| 好男人视频免费观看在线| 九九在线视频观看精品| 美女国产高潮福利片在线看| 18在线观看网站| 亚洲精品第二区| 免费观看在线日韩| 久久人人爽av亚洲精品天堂| 99久久人妻综合| 男女免费视频国产| 青青草视频在线视频观看| 国产爽快片一区二区三区| 亚洲精品视频女| 国产精品嫩草影院av在线观看| 亚洲性久久影院| 91aial.com中文字幕在线观看| 黄色毛片三级朝国网站| 亚洲图色成人| 蜜臀久久99精品久久宅男| 亚洲欧美清纯卡通| 国产欧美亚洲国产| 国产精品一区二区在线不卡| 少妇精品久久久久久久| 成人手机av| 黄色视频在线播放观看不卡| 久久久久久久久久久免费av| 两个人看的免费小视频| 久久99蜜桃精品久久| 亚洲av免费高清在线观看| 五月开心婷婷网| 色视频在线一区二区三区| h视频一区二区三区| 亚洲精品色激情综合| 国产色婷婷99| 亚洲国产精品一区三区| 成人免费观看视频高清| 日本欧美国产在线视频| 最新的欧美精品一区二区| 婷婷成人精品国产| 日本wwww免费看| 天天影视国产精品| 一级片'在线观看视频| 免费黄频网站在线观看国产| 亚洲欧洲国产日韩| 少妇 在线观看| 成年人免费黄色播放视频| 亚洲美女视频黄频| 美女国产高潮福利片在线看| 欧美日韩综合久久久久久| 视频中文字幕在线观看| 国产精品国产av在线观看| 我的女老师完整版在线观看| 欧美 日韩 精品 国产| 一区二区三区精品91| 水蜜桃什么品种好| 免费女性裸体啪啪无遮挡网站| 91精品伊人久久大香线蕉| 99re6热这里在线精品视频| 亚洲欧洲精品一区二区精品久久久 | 久久久精品免费免费高清| 毛片一级片免费看久久久久| 少妇高潮的动态图| 亚洲欧美清纯卡通| freevideosex欧美| 亚洲美女搞黄在线观看| 国产精品人妻久久久影院| 国产不卡av网站在线观看| 校园人妻丝袜中文字幕| 大陆偷拍与自拍| 午夜福利视频精品| 纵有疾风起免费观看全集完整版| 亚洲成人一二三区av| 亚洲伊人久久精品综合| 精品视频人人做人人爽| xxxhd国产人妻xxx| 久久人妻熟女aⅴ| 一级黄片播放器| 少妇被粗大猛烈的视频| 国产又爽黄色视频| 欧美日韩视频高清一区二区三区二| 五月开心婷婷网| 日本欧美国产在线视频| 成人手机av| 日日啪夜夜爽| 精品一区在线观看国产| www.av在线官网国产| 国产亚洲精品久久久com| 最新的欧美精品一区二区| 狠狠精品人妻久久久久久综合| 国产欧美另类精品又又久久亚洲欧美| 久久久欧美国产精品| 一二三四中文在线观看免费高清| 大话2 男鬼变身卡| 亚洲精品中文字幕在线视频| 免费不卡的大黄色大毛片视频在线观看| 热99久久久久精品小说推荐| 丝袜美足系列| 男女边摸边吃奶| 精品亚洲成国产av| 国产成人免费观看mmmm| kizo精华| 丰满少妇做爰视频| 最新的欧美精品一区二区| 免费黄网站久久成人精品| 国产成人精品婷婷| 欧美人与善性xxx| 国语对白做爰xxxⅹ性视频网站| 天天操日日干夜夜撸| 黄色 视频免费看| 亚洲精品乱久久久久久| 亚洲性久久影院| 人妻人人澡人人爽人人| av在线观看视频网站免费| 亚洲国产欧美日韩在线播放| 欧美日韩视频精品一区| 丝袜在线中文字幕| freevideosex欧美| 狂野欧美激情性bbbbbb| 欧美日韩综合久久久久久| 久久97久久精品| 久久久欧美国产精品| 国产男女内射视频| 久久久精品免费免费高清| 久久影院123| 欧美日韩综合久久久久久| 亚洲精品成人av观看孕妇| 亚洲国产av新网站| 久久久欧美国产精品| 好男人视频免费观看在线| 1024视频免费在线观看| 久久久久久人妻| 男女边吃奶边做爰视频| 久久精品人人爽人人爽视色| 欧美xxxx性猛交bbbb| 亚洲国产av新网站| 亚洲精华国产精华液的使用体验| 91精品伊人久久大香线蕉| 日韩人妻精品一区2区三区| 久久久久久人妻| 免费av中文字幕在线| 成人亚洲精品一区在线观看| 亚洲精品第二区| av片东京热男人的天堂| 人妻一区二区av| 国产亚洲精品久久久com| 亚洲色图 男人天堂 中文字幕 | 日本色播在线视频| 在线观看人妻少妇| 日韩成人av中文字幕在线观看| 永久网站在线| 女的被弄到高潮叫床怎么办| 国产精品三级大全| 久久女婷五月综合色啪小说| 日韩一区二区三区影片| 亚洲精品乱久久久久久| 国产在线视频一区二区| a 毛片基地| 日本wwww免费看| 亚洲精品一二三| 午夜老司机福利剧场| 91国产中文字幕| 天天影视国产精品| 啦啦啦视频在线资源免费观看| 亚洲精品aⅴ在线观看| 久久人人97超碰香蕉20202| 亚洲欧洲日产国产| 国产高清不卡午夜福利| 久久精品久久久久久久性| xxx大片免费视频| 欧美激情极品国产一区二区三区 | 亚洲色图 男人天堂 中文字幕 | 亚洲av.av天堂| 熟女av电影| 亚洲av免费高清在线观看| 少妇被粗大猛烈的视频| 国产午夜精品一二区理论片| 一本久久精品| av免费观看日本| 日韩欧美一区视频在线观看| 精品国产国语对白av| 国产精品成人在线| 美女中出高潮动态图| 99久久综合免费| 91国产中文字幕| 久久久久久久亚洲中文字幕| 亚洲精品av麻豆狂野| 日韩av在线免费看完整版不卡| 涩涩av久久男人的天堂| xxxhd国产人妻xxx| a级毛片在线看网站| 乱码一卡2卡4卡精品| 日本爱情动作片www.在线观看| 久久久久久久久久久久大奶| 老熟女久久久| 黄片播放在线免费| 交换朋友夫妻互换小说| 国产乱人偷精品视频| 人人妻人人添人人爽欧美一区卜| 成人毛片a级毛片在线播放| 高清在线视频一区二区三区| 人人妻人人添人人爽欧美一区卜| 22中文网久久字幕| 日本91视频免费播放| 精品少妇久久久久久888优播| 国产永久视频网站| 国产亚洲av片在线观看秒播厂| 啦啦啦中文免费视频观看日本| 1024视频免费在线观看| 国产精品久久久av美女十八| 一区二区三区四区激情视频| 久久久亚洲精品成人影院| 少妇精品久久久久久久| 国产1区2区3区精品|