• <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)
    九九久久精品国产亚洲av麻豆| 综合色丁香网| 国产男女超爽视频在线观看| 久久这里有精品视频免费| 老女人水多毛片| 99九九在线精品视频 | 亚洲精品色激情综合| 制服丝袜香蕉在线| 亚洲第一av免费看| 狠狠精品人妻久久久久久综合| 在线观看一区二区三区激情| 国产精品一区二区在线观看99| 国产一区二区三区av在线| 久久人人爽人人爽人人片va| 18禁在线无遮挡免费观看视频| 久久人人爽av亚洲精品天堂| 一本色道久久久久久精品综合| 人人妻人人爽人人添夜夜欢视频 | 国产日韩欧美在线精品| 啦啦啦啦在线视频资源| 乱系列少妇在线播放| 中文字幕免费在线视频6| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 亚州av有码| 国产亚洲欧美精品永久| 亚洲美女视频黄频| 91久久精品电影网| 久久鲁丝午夜福利片| 在线观看免费视频网站a站| 赤兔流量卡办理| 国产成人精品无人区| 黑人高潮一二区| 综合色丁香网| av国产久精品久网站免费入址| 亚洲精品国产成人久久av| av免费观看日本| 高清午夜精品一区二区三区| 人人妻人人看人人澡| 国产午夜精品一二区理论片| 伊人亚洲综合成人网| 97超视频在线观看视频| 午夜免费鲁丝| 好男人视频免费观看在线| 青春草亚洲视频在线观看| 欧美bdsm另类| 久久久久久久久久久久大奶| 熟女电影av网| 国产精品一区二区在线观看99| 精品熟女少妇av免费看| 九色成人免费人妻av| 在现免费观看毛片| h视频一区二区三区| 国产精品一区二区性色av| 日韩三级伦理在线观看| 欧美日韩视频高清一区二区三区二| 精品午夜福利在线看| 一级毛片 在线播放| 男女啪啪激烈高潮av片| 三级经典国产精品| 尾随美女入室| 人人妻人人爽人人添夜夜欢视频 | 亚洲久久久国产精品| 精品人妻偷拍中文字幕| 18禁在线播放成人免费| 欧美日韩亚洲高清精品| 波野结衣二区三区在线| 午夜激情久久久久久久| 麻豆乱淫一区二区| 欧美人与善性xxx| 高清黄色对白视频在线免费看 | 亚洲综合精品二区| 久久久精品94久久精品| a级毛片免费高清观看在线播放| 久久久国产一区二区| 国国产精品蜜臀av免费| 大码成人一级视频| 亚洲国产精品一区三区| 在线亚洲精品国产二区图片欧美 | 日韩一本色道免费dvd| 日韩不卡一区二区三区视频在线| 3wmmmm亚洲av在线观看| 中文天堂在线官网| av天堂中文字幕网| 大香蕉久久网| 国产探花极品一区二区| 亚洲av二区三区四区| 3wmmmm亚洲av在线观看| 青春草亚洲视频在线观看| 精品卡一卡二卡四卡免费| 国产亚洲午夜精品一区二区久久| 亚洲伊人久久精品综合| av在线观看视频网站免费| 69精品国产乱码久久久| 日本午夜av视频| 日本wwww免费看| 性色avwww在线观看| 在线观看一区二区三区激情| 免费观看av网站的网址| 国产精品成人在线| 亚洲美女黄色视频免费看| 亚洲电影在线观看av| 国产免费福利视频在线观看| 亚洲欧美精品专区久久| 深夜a级毛片| 最黄视频免费看| 丰满人妻一区二区三区视频av| 亚洲av电影在线观看一区二区三区| 国产片特级美女逼逼视频| 97在线人人人人妻| 在线播放无遮挡| 国产av一区二区精品久久| 精品亚洲成国产av| 免费看日本二区| 中文字幕人妻熟人妻熟丝袜美| 一个人免费看片子| av视频免费观看在线观看| 欧美日本中文国产一区发布| 两个人免费观看高清视频 | 国产精品人妻久久久久久| 久久这里有精品视频免费| 另类亚洲欧美激情| 日本黄大片高清| 91aial.com中文字幕在线观看| 国产 一区精品| 日韩欧美 国产精品| 我的老师免费观看完整版| 午夜日本视频在线| 免费av不卡在线播放| av福利片在线| 欧美少妇被猛烈插入视频| 一级毛片久久久久久久久女| 国产亚洲91精品色在线| 色哟哟·www| av天堂久久9| 91在线精品国自产拍蜜月| 国语对白做爰xxxⅹ性视频网站| 日韩 亚洲 欧美在线| 伊人久久国产一区二区| 99热这里只有是精品在线观看| 免费看av在线观看网站| 欧美成人午夜免费资源| 久久国产精品男人的天堂亚洲 | 一级毛片我不卡| 久久人人爽人人片av| 天堂俺去俺来也www色官网| 在线观看www视频免费| 水蜜桃什么品种好| 中文乱码字字幕精品一区二区三区| 乱系列少妇在线播放| 丝袜在线中文字幕| 精品久久久噜噜| 蜜臀久久99精品久久宅男| 一区二区三区四区激情视频| 国产91av在线免费观看| 少妇的逼水好多| 国产在线免费精品| 99国产精品免费福利视频| 欧美激情极品国产一区二区三区 | 能在线免费看毛片的网站| 中文字幕精品免费在线观看视频 | 观看av在线不卡| 亚洲综合色惰| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡动漫免费视频| 午夜激情久久久久久久| 青青草视频在线视频观看| 成人无遮挡网站| 男人狂女人下面高潮的视频| 熟女av电影| 你懂的网址亚洲精品在线观看| 国产av精品麻豆| 成人免费观看视频高清| av天堂中文字幕网| 成人黄色视频免费在线看| 99久久精品热视频| 亚洲欧洲日产国产| 三上悠亚av全集在线观看 | 亚洲精品日本国产第一区| 亚洲精华国产精华液的使用体验| 精品久久久久久久久亚洲| 久久青草综合色| 国产精品国产三级国产专区5o| 亚洲av综合色区一区| 国产老妇伦熟女老妇高清| 丝袜脚勾引网站| 伊人久久国产一区二区| 国产又色又爽无遮挡免| 亚洲国产精品一区二区三区在线| 免费av中文字幕在线| 久久久久久久大尺度免费视频| 国产有黄有色有爽视频| 久久久久视频综合| 大片免费播放器 马上看| 国产美女午夜福利| a 毛片基地| 亚洲av综合色区一区| 日本欧美视频一区| 亚洲国产日韩一区二区| 制服丝袜香蕉在线| 我的女老师完整版在线观看| 免费观看的影片在线观看| 欧美老熟妇乱子伦牲交| 啦啦啦在线观看免费高清www| 久久99热6这里只有精品| 亚洲经典国产精华液单| 偷拍熟女少妇极品色| 伊人久久国产一区二区| 亚洲真实伦在线观看| 天堂8中文在线网| 在线免费观看不下载黄p国产| 蜜桃久久精品国产亚洲av| 国产在线男女| 国产成人一区二区在线| 人妻制服诱惑在线中文字幕| 亚洲欧美一区二区三区黑人 | 老司机影院毛片| 久久久国产一区二区| 极品教师在线视频| 精品少妇黑人巨大在线播放| 交换朋友夫妻互换小说| 国产黄色免费在线视频| 久久精品熟女亚洲av麻豆精品| 久久免费观看电影| 最新的欧美精品一区二区| 国产欧美日韩综合在线一区二区 | 久久久久久伊人网av| av天堂久久9| 日韩欧美精品免费久久| 99热这里只有精品一区| 国产女主播在线喷水免费视频网站| 国产高清有码在线观看视频| 亚洲人成网站在线观看播放| 国产成人精品婷婷| 国产深夜福利视频在线观看| 日本免费在线观看一区| 日本vs欧美在线观看视频 | 亚洲精品乱码久久久v下载方式| 国产深夜福利视频在线观看| 97超碰精品成人国产| 色网站视频免费| 深夜a级毛片| 久久人人爽av亚洲精品天堂| 水蜜桃什么品种好| kizo精华| 97在线人人人人妻| 久久毛片免费看一区二区三区| 久久热精品热| 亚洲av.av天堂| 2018国产大陆天天弄谢| 91久久精品国产一区二区成人| 丝袜脚勾引网站| 免费黄频网站在线观看国产| 欧美人与善性xxx| 美女大奶头黄色视频| 久久久久视频综合| 亚洲精品456在线播放app| 女人久久www免费人成看片| 久久精品国产自在天天线| 一级a做视频免费观看| 午夜免费鲁丝| 最近的中文字幕免费完整| 色视频www国产| 国产成人午夜福利电影在线观看| 欧美最新免费一区二区三区| 午夜av观看不卡| 一区二区三区免费毛片| 一个人看视频在线观看www免费| 18禁裸乳无遮挡动漫免费视频| 最近手机中文字幕大全| 丁香六月天网| 国产探花极品一区二区| 国产免费福利视频在线观看| 亚洲综合色惰| 免费观看无遮挡的男女| 亚洲人与动物交配视频| 久久精品熟女亚洲av麻豆精品| 午夜福利网站1000一区二区三区| 精品久久久久久久久av| 美女内射精品一级片tv| 国产精品国产三级国产专区5o| 丝瓜视频免费看黄片| 久久久精品免费免费高清| 日本黄大片高清| 日韩大片免费观看网站| 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 极品人妻少妇av视频| 在线观看www视频免费| av福利片在线观看| 五月玫瑰六月丁香| 夜夜骑夜夜射夜夜干| 精品99又大又爽又粗少妇毛片| 国产成人精品一,二区| 伦精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 在线观看三级黄色| 久久久久国产精品人妻一区二区| 午夜免费观看性视频| 日本与韩国留学比较| 一边亲一边摸免费视频| 亚洲高清免费不卡视频| 黑人巨大精品欧美一区二区蜜桃 | 日产精品乱码卡一卡2卡三| 日韩中文字幕视频在线看片| 热re99久久国产66热| 亚洲精品一区蜜桃| 色网站视频免费| 久久青草综合色| 97超碰精品成人国产| 五月玫瑰六月丁香| 少妇 在线观看| 国产精品国产av在线观看| 人妻夜夜爽99麻豆av| 最新中文字幕久久久久| 久久久国产一区二区| 王馨瑶露胸无遮挡在线观看| 自线自在国产av| 黑丝袜美女国产一区| 菩萨蛮人人尽说江南好唐韦庄| 99九九在线精品视频 | 天天操日日干夜夜撸| 免费人妻精品一区二区三区视频| 日本vs欧美在线观看视频 | 午夜视频国产福利| 天天躁夜夜躁狠狠久久av| 少妇人妻 视频| 国产精品久久久久成人av| 亚洲欧美成人精品一区二区| 我的女老师完整版在线观看| 一级黄片播放器| 内射极品少妇av片p| 只有这里有精品99| av又黄又爽大尺度在线免费看| 国产免费视频播放在线视频| 三上悠亚av全集在线观看 | 免费av中文字幕在线| 久久久久久久大尺度免费视频| av网站免费在线观看视频| 人人澡人人妻人| av黄色大香蕉| 天堂中文最新版在线下载| 亚洲欧洲国产日韩| 国产综合精华液| 国产精品99久久99久久久不卡 | 18禁在线无遮挡免费观看视频| 成人美女网站在线观看视频| 日本与韩国留学比较| 久久 成人 亚洲| 女的被弄到高潮叫床怎么办| 80岁老熟妇乱子伦牲交| 欧美日韩一区二区视频在线观看视频在线| xxx大片免费视频| 在现免费观看毛片| 国产精品麻豆人妻色哟哟久久| 久久久亚洲精品成人影院| 亚洲高清免费不卡视频| 大又大粗又爽又黄少妇毛片口| 亚洲av日韩在线播放| 五月玫瑰六月丁香| 免费观看的影片在线观看| 老熟女久久久| 精品熟女少妇av免费看| 各种免费的搞黄视频| 乱系列少妇在线播放| 欧美高清成人免费视频www| 久久久久久久久大av| 久久久久精品性色| 午夜免费男女啪啪视频观看| 夜夜爽夜夜爽视频| 国产在线视频一区二区| 成人国产麻豆网| 男男h啪啪无遮挡| 在线观看一区二区三区激情| av网站免费在线观看视频| 一级毛片 在线播放| 欧美日韩av久久| 内射极品少妇av片p| 国产日韩一区二区三区精品不卡 | 日韩在线高清观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 人妻少妇偷人精品九色| 另类亚洲欧美激情| 久久免费观看电影| 黑人猛操日本美女一级片| 久久人人爽人人片av| 人妻一区二区av| 最黄视频免费看| 在线观看免费高清a一片| 在线观看一区二区三区激情| av福利片在线| 在线 av 中文字幕| 国产伦在线观看视频一区| 国产精品无大码| 国产男人的电影天堂91| 亚洲婷婷狠狠爱综合网| 亚洲av在线观看美女高潮| 色网站视频免费| 激情五月婷婷亚洲| 精品亚洲成国产av| 精品国产一区二区久久| 国产高清三级在线| 少妇的逼好多水| 80岁老熟妇乱子伦牲交| 久久综合国产亚洲精品| 伦理电影大哥的女人| 黄色毛片三级朝国网站 | 亚洲av.av天堂| 日韩免费高清中文字幕av| 久久久久久人妻| 看免费成人av毛片| 日韩av不卡免费在线播放| 亚洲精品国产成人久久av| 99热这里只有是精品在线观看| 精品酒店卫生间| 国产视频首页在线观看| 内射极品少妇av片p| 国精品久久久久久国模美| 国产探花极品一区二区| 精品一区在线观看国产| 精品人妻一区二区三区麻豆| 交换朋友夫妻互换小说| 欧美国产精品一级二级三级 | 性高湖久久久久久久久免费观看| 国产伦精品一区二区三区视频9| 亚洲va在线va天堂va国产| 人妻 亚洲 视频| 日本wwww免费看| 久久久久久久国产电影| 丰满人妻一区二区三区视频av| 久久久久久久国产电影| 国产精品成人在线| 午夜老司机福利剧场| 另类亚洲欧美激情| 免费观看a级毛片全部| 国产在视频线精品| 国产精品国产av在线观看| 精品久久久久久久久亚洲| 国产高清有码在线观看视频| 国产色爽女视频免费观看| 国产精品久久久久久久电影| 久久午夜综合久久蜜桃| 欧美成人午夜免费资源| 久久热精品热| 一本一本综合久久| h日本视频在线播放| 一区二区三区四区激情视频| 日韩成人伦理影院| 国产男女内射视频| 18禁动态无遮挡网站| 寂寞人妻少妇视频99o| 日韩三级伦理在线观看| 欧美日韩av久久| 亚洲,欧美,日韩| 十分钟在线观看高清视频www | videos熟女内射| 91aial.com中文字幕在线观看| 狠狠精品人妻久久久久久综合| 视频区图区小说| 特大巨黑吊av在线直播| 国产伦在线观看视频一区| 午夜免费男女啪啪视频观看| 看非洲黑人一级黄片| 老司机亚洲免费影院| 少妇的逼水好多| xxx大片免费视频| 国产一区有黄有色的免费视频| 精品一区二区免费观看| 丝袜在线中文字幕| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 美女中出高潮动态图| 看非洲黑人一级黄片| 亚洲国产色片| 亚洲av二区三区四区| 欧美日韩视频高清一区二区三区二| 国产91av在线免费观看| 国产永久视频网站| 一本一本综合久久| 制服丝袜香蕉在线| 麻豆乱淫一区二区| 最新中文字幕久久久久| 日韩电影二区| 一级毛片电影观看| 精品一品国产午夜福利视频| 中文精品一卡2卡3卡4更新| 美女视频免费永久观看网站| 成人毛片60女人毛片免费| 久久午夜综合久久蜜桃| 欧美一级a爱片免费观看看| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美清纯卡通| 蜜桃在线观看..| 爱豆传媒免费全集在线观看| av在线app专区| 视频区图区小说| 欧美97在线视频| 高清av免费在线| 久久婷婷青草| 国产亚洲91精品色在线| 狠狠精品人妻久久久久久综合| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品一区三区| 又黄又爽又刺激的免费视频.| 午夜av观看不卡| 丝袜脚勾引网站| 熟女电影av网| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 精品国产国语对白av| 免费黄网站久久成人精品| av网站免费在线观看视频| 国产美女午夜福利| 婷婷色麻豆天堂久久| 日本黄色片子视频| 五月伊人婷婷丁香| 在线亚洲精品国产二区图片欧美 | 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 精品国产一区二区久久| 中文在线观看免费www的网站| 在线观看三级黄色| 国产av一区二区精品久久| 国产高清国产精品国产三级| 日日啪夜夜爽| 国产精品秋霞免费鲁丝片| 国产色婷婷99| 亚洲欧洲日产国产| 欧美老熟妇乱子伦牲交| 亚洲精品久久午夜乱码| 大陆偷拍与自拍| 老司机亚洲免费影院| 成年美女黄网站色视频大全免费 | 免费黄网站久久成人精品| 精品一区二区免费观看| av在线老鸭窝| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 亚洲欧美日韩卡通动漫| 插逼视频在线观看| 黄色视频在线播放观看不卡| 国产午夜精品久久久久久一区二区三区| 亚洲成人手机| 丝瓜视频免费看黄片| 一本色道久久久久久精品综合| 亚洲三级黄色毛片| 狂野欧美激情性xxxx在线观看| 国产伦理片在线播放av一区| 久久99热这里只频精品6学生| 日韩一区二区视频免费看| 国产有黄有色有爽视频| 亚洲精品国产av成人精品| 欧美三级亚洲精品| 丝袜脚勾引网站| 午夜日本视频在线| 嘟嘟电影网在线观看| 日韩精品有码人妻一区| 男人狂女人下面高潮的视频| 午夜av观看不卡| 欧美3d第一页| 亚洲熟女精品中文字幕| 中文字幕亚洲精品专区| 亚洲精品,欧美精品| 欧美xxxx性猛交bbbb| 这个男人来自地球电影免费观看 | 99热这里只有是精品50| 爱豆传媒免费全集在线观看| 在线观看一区二区三区激情| 国产国拍精品亚洲av在线观看| 成人午夜精彩视频在线观看| 日韩成人av中文字幕在线观看| 99国产精品免费福利视频| 国产精品久久久久久av不卡| 欧美区成人在线视频| 亚洲婷婷狠狠爱综合网| 精品国产一区二区三区久久久樱花| 国产在视频线精品| 欧美三级亚洲精品| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 亚洲人与动物交配视频| 国内少妇人妻偷人精品xxx网站| 亚洲精品第二区| 午夜日本视频在线| 日韩 亚洲 欧美在线| 亚洲精品自拍成人| 天堂俺去俺来也www色官网| 国产精品成人在线| 欧美成人精品欧美一级黄| 看非洲黑人一级黄片| 久久久久久久国产电影| 丰满人妻一区二区三区视频av| 午夜视频国产福利| 日韩免费高清中文字幕av| 成人黄色视频免费在线看| 秋霞伦理黄片| 最后的刺客免费高清国语| 亚洲激情五月婷婷啪啪| 美女视频免费永久观看网站| 国产一区二区在线观看日韩| 少妇精品久久久久久久| 激情五月婷婷亚洲| 免费大片黄手机在线观看| 日韩免费高清中文字幕av| 亚洲精品亚洲一区二区| 亚洲av不卡在线观看| 亚洲精品,欧美精品| 蜜桃久久精品国产亚洲av| 久久久国产一区二区| 精品人妻熟女毛片av久久网站| av国产久精品久网站免费入址| 精品人妻偷拍中文字幕| 中文字幕亚洲精品专区|