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

    內(nèi)蒙古東珺鉛鋅銀礦床質(zhì)量平衡計算與圍巖蝕變特征

    2018-09-07 07:54:18段天緒溫守欽周鵬朱恩靜崔顯德孫家全趙桐魏燦坤
    中南大學學報(自然科學版) 2018年8期
    關(guān)鍵詞:化帶云母熱液

    段天緒,溫守欽,周鵬,朱恩靜,崔顯德,孫家全,趙桐,魏燦坤

    ?

    內(nèi)蒙古東珺鉛鋅銀礦床質(zhì)量平衡計算與圍巖蝕變特征

    段天緒1,溫守欽1,周鵬2, 3,朱恩靜1,崔顯德1,孫家全1,趙桐1,魏燦坤1

    (1. 東北大學 深部金屬礦山安全開采教育部重點實驗室,遼寧 沈陽,110819; 2. 中國海洋大學 海洋地球科學學院,山東 青島,266100;3.國家海洋局第二海洋研究所 海底科學重點實驗室,浙江 杭州,310012)

    為研究內(nèi)蒙古東珺鉛鋅銀礦床圍巖蝕變類型與蝕變分帶地球化學特征,將近礦、遠礦蝕變圍巖與紫紅色安山巖進行對比,以Al2O3為惰性組分,利用標準化Isocon圖解法對圍巖蝕變過程中物質(zhì)組分帶入、帶出進行質(zhì)量平衡計算;分析微量元素的遷移機制。研究結(jié)果表明:蝕變的空間位置關(guān)系自內(nèi)向外可以劃分為3個蝕變帶即鉀硅化?絹云母化帶、絹英巖化帶、青磐巖化帶,其中與礦床成礦關(guān)系最為密切的是絹英巖化帶;巖漿熱液沿裂隙侵位交代圍巖,近礦蝕變圍巖中各組分質(zhì)量變化都非常明顯,MgO,P2O5,TFe2O3,CaO和Na2O組分都在蝕變時隨熱液大量遷出,而SiO2和K2O則在熱液對圍巖進行交代過程中保留下來;近礦蝕變圍巖多發(fā)生硅化、鉀化、絹云母化、黑云母化蝕變;近礦蝕變圍巖中明顯減少的組分TFe2O3,MgO和P2O5在遠礦蝕變圍巖中都略有增加,CaO質(zhì)量分數(shù)顯著上升,SiO2質(zhì)量分數(shù)減少;遠礦蝕變多發(fā)生綠簾石化、綠泥石化、碳酸鹽化;近礦蝕變圍巖輕重稀土分餾最強,遠礦蝕變圍巖輕重稀土分餾減弱;近礦蝕變圍巖主要富集U和Th等高場強元素,多強烈虧損Ba和Sr等大離子親石元素,其中La,Ce,Sr,Nd,Zr,Sm和Y元素從近礦圍巖中組分流失而在遠礦圍巖處相對富集,具備火山—次火山熱液礦床地球化學與圍巖蝕變特征。

    東珺鉛鋅銀礦床;質(zhì)量平衡;圍巖蝕變;標準化Isocon圖解法

    圍巖蝕變是指在熱液成礦過程中,熱液流體或溶液在巖石中流動循環(huán),導致巖石物理化學性質(zhì)(包括礦物組合、巖石結(jié)構(gòu)等)改變的復雜過程[1]。對礦區(qū)巖石蝕變過程中元素行為的研究可以查明蝕變過程中物質(zhì)的帶入和帶出[2],是判斷成礦物質(zhì)來源、元素遷移方式和沉淀機制等問題的重要途徑之一。內(nèi)蒙古東珺鉛鋅銀礦床近幾年找礦工作取得重大突破,其中以X號礦帶礦體規(guī)模最大,礦石品位最高。本文作者結(jié)合野外實際地質(zhì)現(xiàn)象與室內(nèi)巖礦鑒定,總結(jié)圍巖蝕變類型、空間蝕變分布規(guī)律;對礦床大面積出露紫紅色安山巖、近礦蝕變圍巖、遠礦蝕變圍巖進行元素測試分析;利用標準化Isocon圖解法[3]進行質(zhì)量平衡計算,以便揭示成礦過程中圍巖所遭受的物理化學變化,為研究該礦的成礦物質(zhì)來源以及為成礦過程提供依據(jù)。

    1 區(qū)域地質(zhì)概況

    內(nèi)蒙古東珺鉛鋅銀礦床位于內(nèi)蒙古自治區(qū)額爾古納市北東方向20 km處,研究區(qū)及其周緣大地構(gòu)造分區(qū)見圖1。地理坐標為E120°17′00″~120°23′00″;N50°21′30″~50°23′00″[4]。區(qū)域性NE向得爾布干深斷裂從本區(qū)的北西側(cè)通過[5?6],成礦區(qū)隸屬東烏旗—興安成礦帶,西部與得爾布干成礦帶毗鄰[7?8]。該礦床為海拉爾侏羅紀火山斷陷盆地北部的火山—次火山熱液脈型鉛鋅銀礦床。

    研究區(qū)內(nèi)大面積分布侏羅系中統(tǒng)塔木蘭溝組地層,以中基性火山巖為主,其巖性以安山巖、粗安巖、凝灰?guī)r與沉凝灰?guī)r為主[9]。同時出露小面積侏羅系上統(tǒng)滿克頭鄂博組地層,其主要巖性為角礫熔巖、安山質(zhì)晶屑巖屑角礫巖、晶屑凝灰?guī)r與集塊巖等。礦床周邊構(gòu)造以線性構(gòu)造為主,多與NEE向的根河大斷裂有關(guān)。礦區(qū)位于根河大斷裂的北側(cè)。

    研究區(qū)燕山期的構(gòu)造巖漿活動強烈,與之相關(guān)的內(nèi)生有色金屬礦產(chǎn)豐富,礦床常呈帶狀分布。

    2 礦床地質(zhì)特征

    東珺礦區(qū)先后發(fā)現(xiàn)10條礦化帶(編號分別為I,II,III,IV,V,VI,VII,VIII,IX和X),近幾年發(fā)現(xiàn)的X礦帶礦體規(guī)模最大,礦石品位最高。東珺鉛鋅銀礦床地質(zhì)簡圖見圖2。所發(fā)現(xiàn)礦體均賦存于侏羅統(tǒng)塔木蘭溝組火山巖中,其中與蝕變安山巖及沉凝灰?guī)r關(guān)系最為密切。除侏羅系塔木蘭溝組火山巖外,侵入巖不甚發(fā)育,已圈定礦體可達80多條,多被殘坡積物覆蓋,為隱伏—半隱伏礦體。

    礦區(qū)發(fā)育NW和NE向2組斷裂構(gòu)造,其中以NW向斷裂構(gòu)造為主,也是重要的控礦和容礦構(gòu)造[10]。礦體產(chǎn)狀多受NW向構(gòu)造影響,礦體總體呈NWW走向,局部走向為NE,傾向多為SW,傾角為60°~80°。

    礦體長度及厚度差別較大,礦體長度大多為50~500 m,厚度則幾米至幾十米不等。礦體多呈脈狀、網(wǎng)脈狀與透鏡狀,局部呈囊狀、似層狀。

    礦石中金屬元素多以硫化物形式存在,主要礦石礦物為黃鐵礦、方鉛礦、閃鋅礦,其次為黃銅礦、毒砂、磁黃鐵礦、銀銻黝銅礦、硫銻銅銀礦等。脈石礦物主要為石英、長石、絹云母、方解石等。礦石結(jié)構(gòu)有半自形?它形粒狀結(jié)晶結(jié)構(gòu)、填隙結(jié)構(gòu)、交代殘余結(jié)構(gòu)、交代溶蝕結(jié)構(gòu)等。礦石構(gòu)造以塊狀、浸染狀、細(網(wǎng))脈狀構(gòu)造為主。

    1—主要斷裂及編號;2—額爾古納地塊;3—鄂倫春晚華力西褶皺帶;4—東烏旗早華力西褶皺帶; 5—東烏旗南晚華力西褶皺帶;6—海拉爾—拉布大林中生代火山巖盆地;7—根河中生代火山巖盆地;8—城鎮(zhèn)。 主要斷裂帶:①—呼倫湖西—額爾古納河斷裂;②—得爾布干斷裂;③—鄂倫春—頭道橋斷裂; ④—大興安嶺主脊斷裂;⑤—嫩江斷裂;⑥—阿榮旗斷裂

    圖2 東珺鉛鋅銀礦床地質(zhì)簡圖

    3 圍巖蝕變特征

    3.1 圍巖蝕變類型

    礦區(qū)圍巖蝕變發(fā)育,可見硅化、絹云母化、碳酸鹽化、黏土化、綠泥石化、綠簾石化、鉀化等,東珺礦區(qū)圍巖蝕變類型鏡下顯微照片見圖3。

    1) 硅化。在工作區(qū)內(nèi)十分常見,與礦化關(guān)系十分密切。表現(xiàn)為細粒石英集合體呈脈狀或團粒狀集合體交代巖石,或沿巖石裂隙充填交代巖石。硅化后使巖石顏色變淺、硬度變大。硅化強烈時,可見幾乎全由石英構(gòu)成的硅化脈,經(jīng)常可見其與碳酸鹽化、絹云母化等蝕變共同出現(xiàn)(見圖3(a))。

    2) 絹云母化。蝕變作用普遍,主要表現(xiàn)為巖石中的長石蝕變?yōu)榻佋颇?,絹云母呈細小鱗片狀定向排列(見圖3(b))。經(jīng)常與硅化和黏土化相伴生。

    3) 碳酸鹽化。主要表現(xiàn)為方解石等碳酸鹽類礦物呈粒狀集合體或浸染狀交代安山巖等火山巖中的角閃石、輝石和斜長石礦物(見圖3(c));或呈脈狀交代巖石。蝕變強烈時方解石脈較寬,在礦帶內(nèi)經(jīng)??梢娞妓猁}類礦物組成的脈和硅化脈與礦石礦物伴生,碳酸鹽化與礦化關(guān)系十分密切。

    4) 黏土化。工作區(qū)內(nèi)普遍存在,主要表現(xiàn)為巖石中長石蝕變成黏土礦物組合,蝕變后巖石呈土狀,鏡下礦物表現(xiàn)較臟,其經(jīng)常與絹云母化伴生,黏土化之后難以辨別原巖礦物。

    5) 綠泥石化。工作區(qū)內(nèi)較為普遍,主要表現(xiàn)為巖石中的角閃石、輝石等暗色礦物蝕變成綠泥石(見圖3(d)),綠泥石呈脈狀、浸染狀分布,或與硅質(zhì)、碳酸鹽礦物等充填在火山巖氣孔中形成杏仁構(gòu)造。

    6) 綠簾石化。在工作區(qū)內(nèi)較為常見,表現(xiàn)為巖石中角閃石、輝石等蝕變成綠簾石(見圖3(e)),蝕變后巖石顏色變淺綠色,部分蝕變強烈的巖石,幾乎全由綠簾石組成,新鮮面呈嫩綠色。

    7) 鉀化。鉀化(見圖3(f))是礦區(qū)內(nèi)深部、近礦蝕變類型,處于含礦巖體的中心部位。鉀化巖石中常伴有較強的浸染狀及細脈狀硫化物礦化。鉀化主要表現(xiàn)為新生成的呈細斑狀集合體的正長石交代斜長石斑晶和基質(zhì)。與硅化關(guān)系最為密切,巖石中次生鉀長石和石英多呈不均勻團塊分布,或呈他形晶充填于石英顆粒間,或次生鉀長石和碳酸鹽等礦物一起呈脈狀分布。

    (a) 硅化(+);(b) 絹云母化(+);(c) 碳酸鹽化(+);(d) 綠泥石化(+);(e) 綠簾石化(+);(f) 鉀化(+)

    3.2 圍巖蝕變組合及分帶

    東珺鉛鋅銀礦床圍巖蝕變強烈,在空間分布上具有明顯的分帶性,見表1。東珺礦區(qū)X號帶地表及P剖面蝕變簡圖分別見圖4和圖5。通過對東珺礦區(qū)X號帶縱向(1 500 m深鉆孔)和橫向(地表)系統(tǒng)采樣和綜合研究發(fā)現(xiàn)(見圖3),從深部到淺部,從近礦到遠礦方向,圍巖蝕變依次可以分為鉀硅化—絹云母化帶、絹英巖化帶和青磐巖化帶(見圖4)。

    3.2.1 鉀硅化—絹云母化帶

    該蝕變帶中,其主要蝕變類型包括鉀化、硅化、絹云母化、黑云母化。該蝕變帶在地表不可見,多發(fā)育在礦床深部,為近礦蝕變。原巖主要為安山巖或安山質(zhì)的火山碎屑巖,金屬礦物主要有輝鉬礦、黃鐵礦、黃銅礦,同時還有少量金紅石、磁鐵礦、毒砂、鈦鐵礦等。蝕變巖一般呈青灰色,塊狀構(gòu)造,石英?黃鐵礦脈較發(fā)育,但該蝕變帶內(nèi)礦石品位低。

    鉀硅化—絹云母化帶在鉆孔中主要位于鉆孔 650 m以下位置。該蝕變帶的劃分主要以黑云母的出現(xiàn)及鉀化的發(fā)育為標志而區(qū)別于絹英巖化帶,通常伴隨有脈狀的硅化、鉀長石化以及比較大范圍的絹云母化。

    其蝕變帶形成主要為2個階段。第1階段巖漿熱液上涌,其中的臨界?超臨界流體在冷凝過程中必然與圍巖發(fā)生反應,此時主要發(fā)生鉀交代作用。

    (Ca,Na)Al2Si2O8+2K++4SiO2→2KAlSi3O8+(Ca,Na)2+;(鈣長石) (石英) (鉀長石) (1)

    隨熱液上涌地下水混入,圍巖斑晶、基質(zhì)中的斜長石普遍遭受絹云母化、硅化。該過程中仍舊是K進入圍巖當中,Ca被置換出來,形成鉀硅化?絹云母化帶。

    Ca2(Fe,Mg)4Al2Si7O22(OH)2+2H2O+2KCl+2(Fe2+,(角閃石)

    Mg2+)→2K(Fe2+,Mg2+)3AlSi3O10(OH)2+2HCl+2Ca2+;(黑云母) (2)

    交代作用結(jié)束之后,由于地下水的循環(huán)流動使研究區(qū)發(fā)生廣泛水解作用,在晚期階段黑云母、斜長石等被分解形成一系列含羥基的鱗片狀礦物如綠泥石以及各類云母等。

    2K(Fe2+,Mg2+)3AlSi3O10(OH)2+Ca2++4H++H2→(黑云母)

    (Fe,Mg,Ca)5Al2Si3O10(OH)8+2(Fe,Mg)2++2K+(綠泥石) (3)

    3.2.2 絹英巖化帶

    熱液經(jīng)過鉀化蝕變改造后,在壓力的作用下熱液繼續(xù)沿裂隙向上運移、擴散,密度的差異使上升熱液中更加富含硫和金屬氯化物絡合離子,當溫度、壓力降低到一定程度時,氣體(H2S等)、水分大量溢出,圍巖中的成礦物質(zhì)在地表水的淋濾作用下不斷在熱液中富集,熱液中金屬礦物的濃度、熱液pH都有所升高,絹英巖化發(fā)生,大量金屬離子以硫化物(PbS,ZnS,F(xiàn)eS和CuS等)的形式發(fā)生沉淀。

    表1 東珺礦床圍巖蝕變帶

    圖4 東珺礦區(qū)X號帶地表蝕變簡圖

    圖5 東珺礦區(qū)X號帶P1剖面蝕變帶簡圖

    該蝕變帶的形成與礦體關(guān)系最為密切,礦體往往產(chǎn)生于該蝕變帶內(nèi),在宏觀上主要表現(xiàn)為褪色蝕變,蝕變巖石呈白色、黃白色的淺色調(diào)。在空間上該蝕變帶位于鉀硅化?絹云母化帶外圍,在縱向上鉆孔顯示主要位于孔深298~650 m處。該蝕變帶的形成指示進入了氫交代階段。

    鉀交代基本完結(jié),原巖中的礦物成分如斜長石、角閃石、黑云母等均遭受了強烈的蝕變作用改造。絹云母化主要表現(xiàn)為呈細小粒狀集合體大面積的交代斑晶、基質(zhì)中的斜長石,硅化主要表現(xiàn)為巖石裂隙中與礦體共生的寬大的石英脈以及基質(zhì)中形成的一些細小的次生石英顆粒。在早期的絹英巖化過程中黃鐵礦、毒砂、方鉛礦、閃鋅礦、黝銅礦、黃銅礦等集中沉淀,構(gòu)成了礦區(qū)的主礦體。

    3.2.3 青磐巖化帶

    完成絹英巖化蝕變后,被改造的熱液繼續(xù)沿地表的巖石裂隙向外運移、滲透,產(chǎn)生了大范圍的青磐巖化帶,在蝕變組合上主要表現(xiàn)為綠簾石化、綠泥石化、硅化、碳酸鹽化。在空間上該蝕變帶位于最外側(cè),同時也是研究區(qū)內(nèi)分布最廣的蝕變帶,在鉆孔內(nèi)主要分布于0~?298 m范圍內(nèi)。

    該蝕變帶內(nèi)原巖主要為安山巖及粗安巖,富含Ca質(zhì)的熱液首先在靠近絹英巖化的部位造成了強烈的綠簾石化,將原巖改造為綠簾石巖,整體都為嫩綠色,或在圍巖裂隙中形成粗大的綠簾石脈體,部分杏仁構(gòu)造內(nèi)充填有晶體顆粒較大的綠簾石。

    2Ca2++3(Ca,Na)Al2Si2O8+2H2O→(斜長石)

    2Ca2(Fe,Al)Al2Si3O12(OH)+2H++2Na+(綠簾石)(5)

    熱液中的Na質(zhì)在綠簾石化的過程中再次富集,并進入外側(cè)圍巖的綠泥石化帶中。最后流體中的Mg質(zhì)參與到廣泛的碳酸鹽蝕變帶中,此時,熱液引起的圍巖蝕變結(jié)束,形成蝕變范圍最廣的青磐巖化帶。該蝕變帶能指示礦體大致范圍,對找礦具有一定實踐 意義。

    4 蝕變過程質(zhì)量平衡計算

    4.1 計算方法選擇

    研究圍巖蝕變最主要的方法是通過質(zhì)量平衡計算來推斷蝕變過程中不同組分的遷移規(guī)律[11]。為了計算出蝕變過程中巖石和流體之間物質(zhì)傳遞的信息,須校正巖體總質(zhì)量改變后的全巖地球化學數(shù)據(jù)。因此,質(zhì)量平衡計算的根本目的就是為了消除總質(zhì)量變化帶來的影響[12]。此前,人們主要采用成分-體積圖解法[13]和Isocon圖解法來[14]進行蝕變過程中的質(zhì)量平衡 計算。

    采用成分?體積圖解法計算時,需分別測定原巖和蝕變巖的質(zhì)量、體積、密度和成分,然后根據(jù)不活潑元素計算出體積因子,之后才能計算物質(zhì)遷移量,這大大增加了工作難度和工作量。

    Isocon圖解法放棄“體積因子”轉(zhuǎn)而通過擬合“質(zhì)量因子”對物質(zhì)的遷移量進行計算,在數(shù)據(jù)測量上比較簡單,但它只適合于2個樣品之間的物質(zhì)遷移計算,并不適用于實際中復雜得多的組分遷移計算。

    標準化Isocon圖解法[14]可以克服上述方法的缺陷。該方法可以將多個樣品的Isocon線整合到1個圖解中,在標準化過程中,可以在各獨立樣品之間建立1個可直接對比的共同參照[15]。已有研究表明,這種方法非常適合于熱液交代作用過程中所產(chǎn)生的圍巖蝕變[16]。通常認為圍巖蝕變過程中的惰性礦物主要有P2O5,TiO2,MnO和Al2O3等[17],其中常用的有Al2O3和TiO2[18?19]。本次測試中P2O5,TiO2和MnO的質(zhì)量分數(shù)較低,測試誤差對實驗結(jié)果影響較大??紤]到蝕變過程中蝕變圍巖總體質(zhì)量可能發(fā)生變化,為了精確評價蝕變過程中圍巖組分的帶入、帶出情況,選取Al2O3作為惰性礦物進行計算。

    4.2 計算結(jié)果分析

    元素分析測試工作在澳實礦物分析檢測(廣州)有限公司完成。主量元素采用硼酸鋰熔融,X線熒光光譜分析法測定,檢出限為0.01%。不同巖層主量元素質(zhì)量分數(shù)測試結(jié)果見表2。微量元素采用四酸消解,質(zhì)譜/光譜儀綜合分析法測定,檢出限為0.01 μg/g。不同巖層微量元素質(zhì)量分數(shù)測試結(jié)果見表3。稀土元素采用硼酸鋰熔融,采用等離子質(zhì)譜定量法測定,檢出限為0.01 μg/g。不同巖層稀土元素質(zhì)量分數(shù)測試結(jié)果見表4。

    為了研究東珺礦區(qū)安山巖與蝕變圍巖物質(zhì)轉(zhuǎn)移量,主要對研究區(qū)大面積出露塔木蘭溝組紫紅色安山巖、近礦圍巖與遠礦圍巖建立標準化Isocon圖解。為保證對比具有針對性,其中所選取近礦、遠礦蝕變圍巖均由安山巖蝕變而來。各巖層主量元素平均質(zhì)量分數(shù)見表5。

    惰性組分的Isocon線斜率計算公式如下:

    k=O/A=CA/CO

    式中:k為惰性組分的Isocon線斜率;O為原巖總質(zhì)量;A為蝕變巖總質(zhì)量;CA為惰性組分在蝕變巖中的質(zhì)量分數(shù);CO為惰性組分在原巖中的質(zhì)量分數(shù)。

    將遠礦蝕變圍巖和紫紅色安山巖的Isocon線按比例調(diào)整到與近礦蝕變圍巖和紫紅色安山巖的Isocon線重合,并按調(diào)整比例相應調(diào)整對應的活動組分,這樣就得到1條標準化Isocon線。

    近礦蝕變圍巖選Al2O3為不活潑組分,通過計算可得Al的Isocon線斜率為0.744 0。為了精確評價蝕變過程中圍巖組分的帶入、帶出情況,須先對選定的惰性組分Al2O3質(zhì)量分數(shù)進行校正。

    以近礦蝕變圍巖中Al2O3為基準,將遠礦蝕變圍巖中Al2O3質(zhì)量分數(shù)乘以近礦、遠礦蝕變圍巖中惰性元素Al的Isocon線斜率之比,得到0.754 4,從而得到與近礦蝕變圍巖一致的Al2O3質(zhì)量分數(shù)標準化值。同理,將紫紅色安山巖也轉(zhuǎn)化為與近礦蝕變圍巖一致的Al2O3質(zhì)量分數(shù)標準化值。不同巖層主量元素平均質(zhì)量分數(shù)標準化值見表6。

    表2 不同巖層主量元素質(zhì)量分數(shù)

    表3 不同巖層微量元素質(zhì)量分數(shù)

    表4 不同巖層稀土元素質(zhì)量分數(shù)

    注:Σ(REE)為稀土元素總質(zhì)量分數(shù);Σ(LREE)為輕稀土元素總質(zhì)量分數(shù);Σ(HREE)為重稀土元素總質(zhì)量分數(shù);(LaN)和(YbN)為鑭元素和釔元素球粒隕石[20]標準化后的比值;δ(·)為元素異常。

    表5 不同巖層主量元素平均質(zhì)量分數(shù)

    表6 不同巖層主量元素平均質(zhì)量分數(shù)標準化值

    4.2.1 主量元素

    利用一致的Al2O3質(zhì)量分數(shù)標準化值繪制標準化Isocon圖解,蝕變圍巖主量元素遷移的標準化Isocon圖解見圖6(其中,橫坐標O為原巖標準化后的主量元素平均質(zhì)量分數(shù),縱坐標A為蝕變巖中標準化后的主量元素平均質(zhì)量分數(shù))。

    蝕變圍巖主量元素遷移定量計算結(jié)果見表7。由表7可知:近礦蝕變圍巖中各組分質(zhì)量變化都非常明顯,MgO,P2O5,TFe2O3,CaO和Na2O組分都在蝕變時隨熱液大量遷出,而SiO2和K2O組分在熱液對圍巖進行交代過程中保留下來,其質(zhì)量分數(shù)有所增加,顯示流體富鉀和硅。結(jié)合礦脈、近礦褪色蝕變巖的礦物特征(礦脈中主要金屬礦物種類為黃鐵礦、閃鋅礦、方鉛礦、黃銅礦、毒砂等硫化物,非金屬礦物主要為石英、碳酸鹽以及發(fā)生了比較徹底的絹云母化的斜長石;圍巖中主要礦物為石英、碳酸鹽、發(fā)生了比較徹底的絹云母化的斜長石以及少量金屬硫化物)可知:熱液中的鐵、鎂、鉛、鋅等金屬元素的在礦脈中以金屬硫化物的形式沉淀,而在熱液交代圍巖過程中富含F(xiàn)e和Mg的角閃石、輝石,以及富Ca和Na的斜長石均發(fā)生強烈的蝕變,其中所富含的元素嚴重流失,P元素也由于磷灰石的蝕變作用進入流體中流失,而硅化、絹云母化的發(fā)生,則使巖石的SiO2和K質(zhì)量分數(shù)明顯上升。

    圖6 蝕變圍巖主量元素遷移的標準化Isocon圖解

    熱液繼續(xù)向外運移交代外部圍巖,產(chǎn)生了大面積的青磐巖化安山巖,遠礦蝕變圍巖發(fā)生的蝕變范圍廣泛,發(fā)生的蝕變強度相對較弱,在近礦蝕變圍巖中明顯減少的組分即TFe2O3,MgO,P2O5和TiO2在遠礦蝕變巖中都略有增加,CaO質(zhì)量分數(shù)顯著上升,SiO2質(zhì)量分數(shù)減少。由于遠礦蝕變過程中原巖中的斜長石普遍發(fā)生了較強的綠簾石化,前期蝕變過程中置換出的Ca質(zhì)、Fe質(zhì)等在該蝕變帶內(nèi)沉積,故Ca和Fe富集而Na和K元素虧損。可見在早期近礦圍巖蝕變過程中已經(jīng)充分析出的組分多在晚期蝕變帶富集,Mg隨熱液進入青磐巖化帶外圍巖石中可引起廣泛的綠泥石化及碳酸鹽化。

    4.2.2 稀土元素

    3類巖石的稀土元素數(shù)據(jù)顯示,稀土總量在遠礦蝕變圍巖最高,在近礦蝕變圍巖最低,近礦蝕變圍巖輕重稀土分餾最強,其次為遠礦蝕變圍巖,紫紅色安山巖稀土元素分餾相對較弱。Eu異常只發(fā)生在近礦蝕變圍巖中。

    對3類巖石各稀土元素質(zhì)量分數(shù)求平均值以探討稀土元素遷移結(jié)果,見表8。從表8可以看出:熱液交代過程給圍巖帶來了大量的稀土元素組分的遷入與遷出,總體表現(xiàn)為近礦蝕變圍巖輕稀土虧損、重稀土富集;遠礦蝕變圍巖輕稀土富集,重稀土 虧損。

    近礦蝕變圍巖輕稀土元素發(fā)生明顯流失。Eu元素易賦存在斜長石中,交代過程中斜長石遭到破壞,故其中的Eu元素發(fā)生流失;Eu元素流失說明了流體具有低pH、低氧逸度的特征。而Pr,Nd,Sm和Gd易賦存于褐簾石、磷灰石、榍石等副礦物中,這些元素的流失說明在交代過程中這些褐簾石、磷灰石、榍石等副礦物也遭受到了不同程度的破壞。近礦蝕變圍巖中其礦物以石英和絹云母為主,石英具有穩(wěn)定的晶體格架,微量元素不能以類質(zhì)同象的形式進入晶格中,只能殘留在石英包裹體內(nèi)部;絹云母的層間可以容納稀土元素,但其穩(wěn)定性遠不如褐簾石、磷灰石、榍石等副礦物,故造成輕稀土元素普遍流失。近礦蝕變圍巖中流失的這些元素多進入遠礦蝕變圍巖中富集并堆積下來,造成其質(zhì)量分數(shù)普遍升高。由此可見,成礦熱液交代近礦圍巖將輕稀土元素帶入而將重稀土元素帶出,同時推動輕稀土元素向外圍擴散,輕重稀土元素分異對找礦具有指示意義。

    4.2.3 微量元素

    蝕變圍巖微量元素遷移定量計算結(jié)果見表9。從表9可以看出:近礦蝕變圍巖主要富集Rb,U,Th,Ta和Nb元素,強烈虧損Ba,Ce,Sr,Nd和Zr元素;蝕變引起的差異反映了熱流體的微量元素特征。巖漿熱液既有殼源物質(zhì)又有幔源物質(zhì)。在巖漿演化過程中,U,Th,Ta和Nb高場強元素多在近礦圍巖富集,Ba元素在蝕變過程易溶于巖漿熱液,組分主要從圍巖中帶出,顯示強烈虧損。Rb與K地球化學性質(zhì)相近,因而參與到鉀礦物的晶格中,一般含在云母與長石中,在近礦圍巖富集。Sr元素與Eu性質(zhì)相似,易進入斜長石晶體中,近礦圍巖斜長石多發(fā)生絹云母化,Sr受巖漿熱液交代作用組分大量帶出并在遠礦蝕變圍巖富集,故近礦圍巖富集Rb,U和Th而貧Ba和Sr。另外,蝕變過程中絹英巖化帶的斜長石等礦物都遭到了較嚴重的破壞,其中La,Ce,Sr,Nd,Zr,Sm和Y元素組分從近礦圍巖中流失而在遠礦圍巖處相對 富集。

    通過微量元素分析可以發(fā)現(xiàn):隨著成礦熱液交代圍巖,高場強元素由于性質(zhì)穩(wěn)定多在近礦圍巖富集,而大離子親石元素由于性質(zhì)活潑多被帶入成礦熱液或運移至遠礦圍巖。微量元素組分遷移受成礦熱液活動影響,與成礦關(guān)系密切。

    表7 蝕變圍巖主量元素遷移定量計算結(jié)果

    注:正值表示主量元素質(zhì)量分數(shù)增加;負值表示主量元素質(zhì)量分數(shù)減少。

    表8 蝕變圍巖稀土元素遷移定量計算結(jié)果

    注:正值表示主量元素質(zhì)量分數(shù)增加;負值表示主量元素質(zhì)量分數(shù)減少。

    表9 蝕變圍巖微量元素遷移定量計算結(jié)果

    注:正值表示主量元素質(zhì)量分數(shù)增加;負值表示主量元素質(zhì)量分數(shù)減少。

    5 結(jié)論

    1) 按照蝕變的空間位置關(guān)系,自內(nèi)向外可以將礦區(qū)劃分為3個蝕變帶即鉀硅化—絹云母化帶、絹英巖化帶、青磐巖化帶。絹英巖化帶與礦床成礦關(guān)系最為密切;鉀硅化?絹云母化帶能見一定輝鉬礦化,但不成規(guī)模;青磐巖化帶沒有明顯礦化現(xiàn)象,但可以圈定礦體大致范圍。

    2) 巖漿沿裂隙侵位交代圍巖,以Al2O3為惰性組分進行質(zhì)量平衡計算發(fā)現(xiàn),近礦蝕變圍巖中各組分質(zhì)量的變化都非常明顯,MgO,P2O5,TFe2O3,CaO和Na2O組分都在蝕變時隨熱液大量遷出,而SiO2和K2O在熱液對圍巖進行交代過程中保留下來,所以,近礦圍巖多發(fā)生硅化、鉀化、絹云母化、黑云母化蝕變。近礦蝕變圍巖中明顯減少的組分TFe2O3,MgO和P2O5在遠礦蝕變圍巖中都略有增加,CaO質(zhì)量分數(shù)顯著上升,SiO2質(zhì)量分數(shù)減少,遠礦蝕變多為綠簾石化、綠泥石化、碳酸鹽化。

    3) 近礦蝕變圍巖輕重稀土分餾最強,其次為遠礦蝕變圍巖,紫紅色安山巖稀土元素分餾相對較弱;Eu異常只發(fā)生在近礦蝕變圍巖中。總體表現(xiàn)為近礦蝕變圍巖輕稀土虧損、重稀土富集。遠礦蝕變圍巖輕稀土富集、重稀土虧損。近礦蝕變圍巖主要富集Rb,U,Th,Ta和Nb元素,強烈虧損Ba,Ce,Sr,Nd和Zr元素,其中的La,Ce,Sr,Nd,Zr,Sm和Y元素從近礦圍巖中組分流失而在遠礦圍巖處相對富集。

    [1] PIRAJNO F. Hydrothermal processes and mineral systems[M]. Berlin, Germany: Springer-Verlag, 2009: 73?84.

    [2] 龔慶杰, 周連壯, 胡楊, 等. 膠東玲瓏金礦田煌斑巖蝕變過程元素遷移行為及其意義[J]. 現(xiàn)代地質(zhì), 2012, 26(5): 1065?1077. GONG Qingjie, ZHOU Lianzhuang, HU Yang, et al. Element transfer behaviors and its application during lamprophyre alteration in the Linglong gold deposit, Jiaodong Peninsula, China[J]. Geoscience, 2012, 26(5): 1065?1077.

    [3] GUO Shun, YE Kai, CHEN Yi, et al. A normalization solution to mass illustration of multiple progressively altered samples using the isocon diagram[J]. Economic Geology, 2009, 104(6): 881?886.

    [4] 張榮慶. 內(nèi)蒙古東珺鉛鋅銀礦綜合找礦模型[J]. 西部探礦工程, 2014(6): 79?82. ZHANG Rongqing. Integrated ore prospecting model for Donjun Pb-Zn-Ag deposit in Inner Mongolia[J].West-china Exploration Engineering, 2014(6): 79?82.

    [5] 佘宏全, 李紅紅, 李進文, 等. 內(nèi)蒙古大興安嶺中北段銅鉛鋅金銀多金屬礦床成礦規(guī)律與找礦方向[J]. 地質(zhì)學報, 2009, 83(10): 1456?1472. SHE Hongquan, LI Honghong, LI Jinwen, et al. The metallogenetical characteristics and prospecting direction of the copper-lead-zinc polymetal deposits in the northern-central Daxing’anling mountain, Inner Mongolia[J]. Acta Geologica Sinica, 2009, 83(10): 1456?1472.

    [6] 張斌. 內(nèi)蒙古東珺鉛鋅銀礦礦床地質(zhì)特征及其成因研究[D]. 北京: 中國地質(zhì)科學院礦產(chǎn)資源研究所, 2011: 1?6. ZHANG Bin. The geological features and genesis of the Dongjun lead-zinc-sliver deposit in Inner Mongolia[D]. Beijing: Chinese Academy of Geological Sciences. Institute of Mineral Resources, 2011: 1?6.

    [7] 張光亮. 額爾古納東珺地區(qū)火山巖巖石地球化學特征研究[D]. 沈陽: 東北大學資源與土木工程學院, 2013: 1?10. ZHANG Guangliang. Study on geochemical characteristics of volcanic rocks from Dongjun region in Ergun[D]. Shenyang: Northeastern University. School of Resources and Civil Engineering, 2013: 1?10.

    [8] 唐克東, 王瑩, 何國琦, 等. 中國東北及鄰區(qū)大陸邊緣構(gòu)造[J]. 地質(zhì)學報, 1995, 69(1): 16?30. TANG Kedong, WANG Ying, HE Guoqi, et al. Continental-margin structure of northeast China and its adjacent areas[J]. Acta Geologica Sinica, 1995, 69(1): 16?30.

    [9] 周鵬, 朱恩靜. 內(nèi)蒙古東珺地區(qū)火山巖噴發(fā)期次探討[J]. 礦物學報, 2015, 35(增刊): 190?191. ZHOU Peng, ZHU Enjing. Discussion on volcanic eruption period in Dongjun, Inner Mongolia[J]. Acta Mineralogica Sinica, 2015, 35(Suppl): 190?191.

    [10] 梁玉偉, 張斌, 田京, 等. 內(nèi)蒙古東珺鉛鋅銀礦床地質(zhì)特征及成因探討[J]. 現(xiàn)代地質(zhì), 2014, 28(6): 1112?1121. LIANG Yuwei, ZHANG Bin, TIAN Jing, et al. Geological characteristics and genesis of the Dongjun Pb-Zn-Ag deposit in Inner Mongolia[J]. Geoscience,2014, 28(6): 1112?1121.

    [11] BOHLKE J K. Comparison of metasomatic reactions between a common CO2-rich vein fluid and diverse wall rocks: intensive variables, mass transfers, and Au mineralization at Alleghany, California[J]. Economic Geology, 1989, 84(2): 291?327.

    [12] 張婷, 彭建堂. 湘西合仁坪鈉長石-石英脈型金礦床圍巖蝕變及質(zhì)量平衡[J]. 地球科學與環(huán)境學報, 2014, 36(4): 32?44. ZHANG Ting, PENG Jiantang. Wall-rock alteration and mass balance of Herenping albite-quartz vein-type gold deposit in the western Hunan[J]. Journal of Earth Sciences and Environment, 2014, 36(4): 32?44.

    [13] GRESENS R L. Composition-volume relationships of metasomatism[J]. Chemical Geology, 1967, 2(1): 47?65.

    [14] GRANT J A. The isocon diagram-a simple solution to Gresens’ equation for metasomatic alteration[J]. Economic Geology, 1986, 81(8): 1976?1982.

    [15] 郭順, 葉凱, 陳意, 等. 開放地質(zhì)體系中物質(zhì)遷移質(zhì)量平衡計算方法介紹[J]. 巖石學報, 2013, 29(5): 1486?1498. GUO Shun, YE Kai, CHEN Yi, et al. Introduction of mass-balance calculation method for component transfer during the opening of a geological system[J]. Acta Petrologica Sinica, 2013, 29(5): 1486?1498.

    [16] GUO Shun, YE Kai, CHEN Yi, et al. Fluid-rock interaction and element mobilization in UHP metabasalt: constraints from an omphacite-epidote vein and host eclogites in the Dabie Orogen[J]. Lithos, 2013, s136/139(4): 145?167.

    [17] 張婷. 湘西合仁坪鈉長石?石英脈型金礦床圍巖蝕變及其與成礦關(guān)系研究[D]. 湖南: 中南大學地球科學與信息物理學院, 2014: 36?37. ZHANG Ting. Wallrock alteration and its relation with gold mineralization in the Herenping albite-quartz lode gold deposit in western Hunan[D]. Hunan: Central South University. School of Earth Science and Information Physics, 2014: 36?37.

    [18] 陳隨海, 韓潤生, 申屠良義, 等. 滇東北礦集區(qū)昭通鉛鋅礦區(qū)蝕變巖分帶及元素遷移特征[J]. 吉林大學學報(地球科學版), 2016, 46(3): 711?721. CHEN Suihai, HAN Runsheng, SHENTU Liangyi, et al. Alteration zoning and geochemical element migration in alteration rock of Zhaotong lead-zinc deposit in northeastern Yunnan mineralization concentration area[J]. Journal of Jilin University (Earth Science Edition), 2016, 46(3): 711?721.

    [19] 周永章, 涂光熾, CHOWN E H, 等. 熱液圍巖蝕變過程中數(shù)學不變量的尋找及元素遷移的定量估計——以廣東河臺金礦田為例[J]. 科學通報, 1994, 39(11): 1026?1028.ZHOU Yongzhang, TU Guangchi, CHOWN E H, et al. The search for mathematical invariants and quantitative estimation of element migration in the process of hydrothermal wall rock alteration: taking Guangdong Hetai gold deposit as an example[J]. Chinese Science Bulletin,1994, 39(11): 1026?1028.

    [20] BOYNTON W V. Cosmochemistry of the rare earth elements: meteorite studies. In: P Henderson (ed.), Rare Earth Element Geochemistry[M]. Amsterdam, Netherlands: Elsevier, 1983: 63?114.

    (編輯 伍錦花)

    Wall rock alteration characteristics and mass balance calculation ofDongjun Pb-Zn-Ag deposit in Inner Mongolia

    DUAN Tianxu1, WEN Shouqin1, ZHOU Peng2, 3, ZHU Enjing1, CUI Xiande1, SUN Jiaquan1, ZHAO Tong1, WEI Cankun1

    (1. Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines, Northeastern University, Shenyang 110819, China; 2. School of Marine Geosciences, Ocean University of China, Qingdao 266100, China; 3. Key Laboratory of Submarine Geosciences, Second Institute of Oceanography, State Oceanic Administration, Hangzhou 310012, China)

    To study the types and geochemical characteristics of wallrock alternation in Dongjun Pb-Zn-Ag deposit, Inner Mongolia, altered near-ore or far-ore wallrock and purplish red andesite were compared. Al2O3was used as the inert component, and mass transfer in the process of wallrock alteration was calculated by normalized Isocon diagram method. In addition, migration mechanism of trace elements was analyzed. The results show that there exist three alteration zones from the inside to the outside, i.e., K-silication-sericitization zone, phyllic alteration zone and propylitization zone. The phyllic alternation zone displays the closest relationships with mineralization. The components in the altered near-ore wallrock fluctuate obviously. A great amount of MgO, P2O5, TFe2O3, CaO and Na2O move out along with hydrotherm during alternation, while SiO2and K2O are preserved in the process of metasomatism of magmatic hydrothermal along fracture emplacement. Potassic, silicified, sericite and biotite alterations are common in the altered near-ore wallrock. Contents of TFe2O3, MgO and P2O5increase slightly in altered far-ore wallrock, which all reduce dramatically in altered near-ore wallrock. There is an obvious increase in CaO mass fraction, but a decrease in SiO2mass fraction in far-ore alternation zone in which epidotization, chloritization and carbonatization are easy to be found. Rare earth elements(REE) fractionation of altered near-ore wallrock is the strongest, and REE fractionation of altered far-ore wallrock is weakened. High field strength elements (e.g. U, Th) are enriched in altered near-ore wall rock, while large ion lithophile elements (e.g. Ba, Sr) are depleted. Elements (e.g. La, Ce, Sr, Nd, Zr, Sm, Y) run from near-ore wallrock and relatively enrich in far-ore wallrock, which shows the geochemical and wallrock alternation characteristics of volcanic?subvolcanic hydrothermal deposit.

    Dongjun Pb-Zn-Ag deposit; mass balance; wallrock alteration; standardized Isocon diagram method

    P594.1

    A

    1672?7207(2018)08?1991?12

    2017?08?13;

    2017?09?29

    國家自然科學基金資助項目(41372098)(Project(41372098) supported by the National Natural Science Foundation of China)

    溫守欽,博士,副教授,從事礦床地球化學研究;E-mail:wensq1963@sina.com

    10.11817/j.issn.1672?7207.2018.08.020

    猜你喜歡
    化帶云母熱液
    基于局部化帶全局追蹤的露天煤礦邊坡臨界滑面確定方法研究
    中國煤炭(2023年11期)2023-12-06 13:36:50
    巖-煤組合體試樣變形場與能量演化特征試驗研究
    裂紋角度對紅砂巖試樣局部化帶演化影響的試驗研究
    等幅循環(huán)加載巖石變形局部化帶位移演化規(guī)律
    煤炭學報(2019年4期)2019-05-08 00:21:22
    納米Al2O3改性多膠粉云母主絕緣的研究
    WS2/TiO2/絹云母復合物的制備及性能表征
    安徽化工(2018年4期)2018-09-03 07:11:50
    塔東熱液地質(zhì)作用機制及對儲層的改造意義
    “氟晶云母特性及應用進展”
    層結(jié)背景下熱液柱演化的實驗模擬*
    海洋與湖沼(2017年6期)2017-03-31 06:18:19
    熱液循環(huán)助采洗井裝置的分析與應用
    石油知識(2016年2期)2016-02-28 16:20:19
    国产91av在线免费观看| videossex国产| 嫩草影院精品99| 国产av不卡久久| 亚洲欧美成人综合另类久久久 | av福利片在线观看| 亚洲成人av在线免费| 国产免费男女视频| 成人高潮视频无遮挡免费网站| av在线亚洲专区| av在线亚洲专区| 日本熟妇午夜| 看免费成人av毛片| 亚洲国产精品合色在线| 看十八女毛片水多多多| www.av在线官网国产| 亚洲,欧美,日韩| 国产精品不卡视频一区二区| 日韩成人伦理影院| 亚洲丝袜综合中文字幕| 免费av观看视频| 熟妇人妻久久中文字幕3abv| 精品少妇黑人巨大在线播放 | 两性午夜刺激爽爽歪歪视频在线观看| 一个人看视频在线观看www免费| 内射极品少妇av片p| 国产成人freesex在线| 国产精品一区二区三区四区免费观看| 看免费成人av毛片| av天堂中文字幕网| 中文亚洲av片在线观看爽| 午夜福利在线观看免费完整高清在 | 黄片无遮挡物在线观看| www日本黄色视频网| 色综合亚洲欧美另类图片| 在线观看一区二区三区| 国产亚洲5aaaaa淫片| 91av网一区二区| 久久精品国产亚洲av香蕉五月| 国产精品免费一区二区三区在线| 深夜精品福利| 日本五十路高清| 日韩精品青青久久久久久| 国产极品天堂在线| 国产 一区精品| 精品久久久久久久久亚洲| 日韩高清综合在线| 男插女下体视频免费在线播放| or卡值多少钱| 亚洲精品色激情综合| 国产成人福利小说| 青春草亚洲视频在线观看| 美女xxoo啪啪120秒动态图| 99国产精品一区二区蜜桃av| 97热精品久久久久久| 最后的刺客免费高清国语| 能在线免费看毛片的网站| 亚洲欧美清纯卡通| 国产人妻一区二区三区在| 色视频www国产| 精品日产1卡2卡| 高清在线视频一区二区三区 | 少妇熟女欧美另类| 免费看日本二区| 日本三级黄在线观看| 亚洲欧美日韩高清在线视频| 99久久精品一区二区三区| 欧美+日韩+精品| 夫妻性生交免费视频一级片| 国内少妇人妻偷人精品xxx网站| av天堂在线播放| 亚洲av男天堂| 国国产精品蜜臀av免费| 国产精品女同一区二区软件| 91久久精品国产一区二区成人| 午夜老司机福利剧场| 老司机影院成人| 国产精品一二三区在线看| 校园春色视频在线观看| 又黄又爽又刺激的免费视频.| 欧美日本视频| 国产精品1区2区在线观看.| 午夜福利在线观看免费完整高清在 | 亚洲国产欧美在线一区| 综合色av麻豆| 亚洲国产精品成人综合色| 一本一本综合久久| 久99久视频精品免费| 亚洲一区高清亚洲精品| 亚洲不卡免费看| 国产伦精品一区二区三区四那| 一个人看视频在线观看www免费| 国产精品一区二区在线观看99 | 欧美zozozo另类| 男人和女人高潮做爰伦理| 三级毛片av免费| 国内久久婷婷六月综合欲色啪| 又黄又爽又刺激的免费视频.| 日韩在线高清观看一区二区三区| 网址你懂的国产日韩在线| 我的老师免费观看完整版| 亚洲av中文av极速乱| 亚洲美女视频黄频| 国产成人福利小说| 欧美激情在线99| 亚洲最大成人手机在线| 99热6这里只有精品| 亚洲第一电影网av| 赤兔流量卡办理| 少妇高潮的动态图| 成人毛片60女人毛片免费| 亚洲国产精品sss在线观看| 爱豆传媒免费全集在线观看| 亚洲av一区综合| 免费观看的影片在线观看| 久久精品久久久久久久性| 亚洲人成网站在线播| 国产av麻豆久久久久久久| 99久久九九国产精品国产免费| 亚洲欧美精品自产自拍| 亚洲欧美精品专区久久| 免费不卡的大黄色大毛片视频在线观看 | 桃色一区二区三区在线观看| 蜜桃久久精品国产亚洲av| 97超碰精品成人国产| 国语自产精品视频在线第100页| ponron亚洲| 国产不卡一卡二| 国产精品三级大全| 久久午夜福利片| 亚洲成人久久爱视频| 亚洲av电影不卡..在线观看| 成人综合一区亚洲| 看十八女毛片水多多多| 91午夜精品亚洲一区二区三区| 亚洲欧美成人综合另类久久久 | 免费看美女性在线毛片视频| 18禁裸乳无遮挡免费网站照片| 欧美丝袜亚洲另类| 在线播放国产精品三级| 成人永久免费在线观看视频| 99久久久亚洲精品蜜臀av| 日韩精品有码人妻一区| 亚洲无线在线观看| 精品久久久久久久末码| 中文字幕精品亚洲无线码一区| 91精品国产九色| 免费人成在线观看视频色| 老熟妇乱子伦视频在线观看| 国产精品综合久久久久久久免费| 免费搜索国产男女视频| 国产精品美女特级片免费视频播放器| 国产av一区在线观看免费| 日韩国内少妇激情av| 男女啪啪激烈高潮av片| 色吧在线观看| а√天堂www在线а√下载| 能在线免费观看的黄片| 尤物成人国产欧美一区二区三区| 国产午夜精品论理片| 亚洲国产精品合色在线| 成人毛片a级毛片在线播放| 少妇被粗大猛烈的视频| 国产成人福利小说| 国产精品久久视频播放| 欧美丝袜亚洲另类| 免费观看精品视频网站| 色综合色国产| 男人舔奶头视频| 日日撸夜夜添| 丝袜美腿在线中文| a级毛色黄片| 99热精品在线国产| av.在线天堂| av国产免费在线观看| 亚洲欧洲日产国产| 一区二区三区免费毛片| 啦啦啦啦在线视频资源| 亚洲在线自拍视频| av福利片在线观看| 精品久久久久久久久亚洲| 亚洲最大成人av| 国产老妇伦熟女老妇高清| 亚洲国产精品久久男人天堂| 欧美+日韩+精品| 男人舔奶头视频| 99在线人妻在线中文字幕| 一个人观看的视频www高清免费观看| 少妇人妻精品综合一区二区 | 你懂的网址亚洲精品在线观看 | 国产蜜桃级精品一区二区三区| 我要搜黄色片| 日韩 亚洲 欧美在线| 亚洲av成人精品一区久久| 变态另类成人亚洲欧美熟女| 亚洲成人中文字幕在线播放| 舔av片在线| 国产美女午夜福利| 男女那种视频在线观看| 九草在线视频观看| 国产精品福利在线免费观看| 色哟哟·www| 亚洲精品国产av成人精品| 午夜精品国产一区二区电影 | 日日摸夜夜添夜夜添av毛片| 91久久精品电影网| 精品午夜福利在线看| 插阴视频在线观看视频| av黄色大香蕉| 国产高清不卡午夜福利| 亚洲美女视频黄频| 国产黄片视频在线免费观看| 一级黄片播放器| 亚洲人成网站在线观看播放| 欧美三级亚洲精品| 国产高清激情床上av| 美女xxoo啪啪120秒动态图| 国产精品一区www在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩东京热| 自拍偷自拍亚洲精品老妇| 国产一区二区在线av高清观看| 亚洲最大成人手机在线| 级片在线观看| 国产精品一区二区在线观看99 | 久久久国产成人免费| 老女人水多毛片| 国产亚洲精品久久久久久毛片| 欧美bdsm另类| а√天堂www在线а√下载| 人妻系列 视频| 小蜜桃在线观看免费完整版高清| 一进一出抽搐gif免费好疼| 黑人高潮一二区| 一区二区三区免费毛片| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 18禁黄网站禁片免费观看直播| 中国国产av一级| 少妇人妻精品综合一区二区 | 国产精品福利在线免费观看| 国国产精品蜜臀av免费| 乱系列少妇在线播放| 赤兔流量卡办理| 亚洲av中文av极速乱| 美女xxoo啪啪120秒动态图| 亚洲成a人片在线一区二区| 97超碰精品成人国产| 此物有八面人人有两片| 嫩草影院精品99| 夜夜夜夜夜久久久久| 国产亚洲精品av在线| 2021天堂中文幕一二区在线观| 我要看日韩黄色一级片| 最近2019中文字幕mv第一页| 久久这里有精品视频免费| 成人毛片60女人毛片免费| 久久久久久国产a免费观看| 日韩 亚洲 欧美在线| 可以在线观看毛片的网站| 一区二区三区高清视频在线| 亚洲内射少妇av| 精品人妻熟女av久视频| 看片在线看免费视频| 日本五十路高清| 日日摸夜夜添夜夜爱| 亚洲久久久久久中文字幕| 亚洲色图av天堂| 国产精品av视频在线免费观看| 日本免费一区二区三区高清不卡| 国产精品,欧美在线| 国产亚洲精品久久久com| 国产在线精品亚洲第一网站| 熟女电影av网| 丰满人妻一区二区三区视频av| 亚洲欧美清纯卡通| 国产精品电影一区二区三区| 欧美一区二区精品小视频在线| 99热只有精品国产| 欧美xxxx黑人xx丫x性爽| 最近最新中文字幕大全电影3| 亚洲在久久综合| 国产三级在线视频| 免费观看人在逋| 如何舔出高潮| 国产精品1区2区在线观看.| 久久久久久久久大av| 一本一本综合久久| 欧美潮喷喷水| 男人舔奶头视频| 12—13女人毛片做爰片一| 亚洲精品日韩av片在线观看| av在线播放精品| 亚洲精品乱码久久久v下载方式| 成人亚洲欧美一区二区av| 色吧在线观看| 99久国产av精品国产电影| 亚洲欧美日韩高清在线视频| 国产免费一级a男人的天堂| 麻豆国产97在线/欧美| 秋霞在线观看毛片| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 久久6这里有精品| 色尼玛亚洲综合影院| 美女脱内裤让男人舔精品视频 | 插阴视频在线观看视频| 精品久久久噜噜| av在线蜜桃| 狂野欧美激情性xxxx在线观看| 啦啦啦啦在线视频资源| 大香蕉久久网| 色吧在线观看| 99在线视频只有这里精品首页| 国产精品野战在线观看| 人人妻人人澡欧美一区二区| 91在线精品国自产拍蜜月| 在线播放国产精品三级| 麻豆久久精品国产亚洲av| 国产精品女同一区二区软件| 成人亚洲精品av一区二区| 变态另类成人亚洲欧美熟女| 亚洲国产欧美在线一区| 国产麻豆成人av免费视频| 成人永久免费在线观看视频| 少妇被粗大猛烈的视频| 亚洲va在线va天堂va国产| 成人一区二区视频在线观看| av天堂在线播放| 亚洲aⅴ乱码一区二区在线播放| 中文在线观看免费www的网站| 在线天堂最新版资源| 久久精品综合一区二区三区| 干丝袜人妻中文字幕| 一个人看的www免费观看视频| 97超视频在线观看视频| 亚洲人与动物交配视频| 国产精品99久久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 春色校园在线视频观看| 美女被艹到高潮喷水动态| 国产 一区精品| 亚洲va在线va天堂va国产| av专区在线播放| 国产真实伦视频高清在线观看| 亚洲高清免费不卡视频| 成人美女网站在线观看视频| 亚洲激情五月婷婷啪啪| 日韩成人av中文字幕在线观看| 99久久人妻综合| 天天躁夜夜躁狠狠久久av| 少妇裸体淫交视频免费看高清| 99久久人妻综合| 天天躁夜夜躁狠狠久久av| 天美传媒精品一区二区| 又爽又黄a免费视频| 国产爱豆传媒在线观看| 最近中文字幕高清免费大全6| 精品人妻一区二区三区麻豆| 悠悠久久av| 欧美zozozo另类| 精品久久久噜噜| 少妇人妻一区二区三区视频| 免费观看的影片在线观看| 久久久久免费精品人妻一区二区| 久久精品夜夜夜夜夜久久蜜豆| 一本一本综合久久| 久久亚洲国产成人精品v| 亚洲美女视频黄频| 国产真实乱freesex| 国产精品.久久久| 床上黄色一级片| 蜜桃亚洲精品一区二区三区| 夜夜看夜夜爽夜夜摸| 一个人免费在线观看电影| 精品久久久久久久久久久久久| 亚洲精品亚洲一区二区| 亚洲五月天丁香| 免费一级毛片在线播放高清视频| 免费看a级黄色片| 最后的刺客免费高清国语| 国产精品人妻久久久久久| www.色视频.com| 51国产日韩欧美| 黄色一级大片看看| 亚洲精品粉嫩美女一区| 中文字幕制服av| 日韩人妻高清精品专区| 18禁在线无遮挡免费观看视频| 美女被艹到高潮喷水动态| 黄片wwwwww| 国产精品一区www在线观看| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 午夜福利成人在线免费观看| 色哟哟哟哟哟哟| 欧美zozozo另类| kizo精华| 久久久久久伊人网av| 国产视频内射| 天堂中文最新版在线下载 | 成人漫画全彩无遮挡| 级片在线观看| 国产精品福利在线免费观看| 99久久精品国产国产毛片| 精品少妇黑人巨大在线播放 | 九色成人免费人妻av| 校园春色视频在线观看| 国产高潮美女av| 日韩欧美精品免费久久| 波多野结衣巨乳人妻| 久久久久久久久大av| 成年免费大片在线观看| 亚洲第一区二区三区不卡| 精品免费久久久久久久清纯| 国产色婷婷99| 成人国产麻豆网| 非洲黑人性xxxx精品又粗又长| 成人亚洲欧美一区二区av| 一级毛片久久久久久久久女| 最近最新中文字幕大全电影3| 黄片无遮挡物在线观看| 成人av在线播放网站| www.av在线官网国产| 成人特级av手机在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品乱码久久久久久按摩| 日韩高清综合在线| 国产免费男女视频| 少妇的逼好多水| 别揉我奶头 嗯啊视频| av在线观看视频网站免费| 国内精品宾馆在线| .国产精品久久| 综合色丁香网| 可以在线观看毛片的网站| 日韩制服骚丝袜av| 青春草国产在线视频 | 少妇的逼好多水| 在线播放国产精品三级| 日本与韩国留学比较| 99在线视频只有这里精品首页| 网址你懂的国产日韩在线| 免费大片18禁| 99精品在免费线老司机午夜| 久久久久久伊人网av| 免费看日本二区| 亚洲av中文av极速乱| 男人和女人高潮做爰伦理| 老熟妇乱子伦视频在线观看| 免费大片18禁| 亚洲丝袜综合中文字幕| 国产精品免费一区二区三区在线| 美女大奶头视频| 精品一区二区三区视频在线| 婷婷精品国产亚洲av| 亚洲av一区综合| 成人午夜精彩视频在线观看| 精品午夜福利在线看| 99久久人妻综合| 看非洲黑人一级黄片| 午夜激情欧美在线| 精品一区二区三区人妻视频| 午夜福利在线观看吧| 黄色配什么色好看| 偷拍熟女少妇极品色| 18+在线观看网站| 一边亲一边摸免费视频| 男女做爰动态图高潮gif福利片| 熟女人妻精品中文字幕| 久久精品91蜜桃| 天天躁日日操中文字幕| 日本熟妇午夜| 免费看美女性在线毛片视频| 成人午夜精彩视频在线观看| 亚洲18禁久久av| 国产大屁股一区二区在线视频| 18禁黄网站禁片免费观看直播| 久久精品国产亚洲网站| 在线观看av片永久免费下载| 国内精品宾馆在线| 插逼视频在线观看| 国产成人a∨麻豆精品| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品综合一区二区三区| 天堂网av新在线| 国产私拍福利视频在线观看| 久久久久久久久久久丰满| 婷婷色综合大香蕉| 国产一级毛片在线| 久久99精品国语久久久| 一进一出抽搐动态| 亚洲欧美精品自产自拍| 午夜老司机福利剧场| 夫妻性生交免费视频一级片| 成人午夜高清在线视频| 亚洲,欧美,日韩| 亚洲最大成人av| 麻豆一二三区av精品| av又黄又爽大尺度在线免费看 | 国产淫片久久久久久久久| 全区人妻精品视频| 春色校园在线视频观看| 国产精品一区二区在线观看99 | 久久久久久久久中文| 久久久久久国产a免费观看| 青春草视频在线免费观看| 亚州av有码| 亚洲人成网站高清观看| 国产成人91sexporn| 国产久久久一区二区三区| 亚洲欧美成人综合另类久久久 | 国产色爽女视频免费观看| 91麻豆精品激情在线观看国产| 亚洲美女搞黄在线观看| 热99在线观看视频| 校园春色视频在线观看| av免费观看日本| 晚上一个人看的免费电影| 欧美bdsm另类| 九九爱精品视频在线观看| 看非洲黑人一级黄片| 午夜精品国产一区二区电影 | 日本色播在线视频| 晚上一个人看的免费电影| 3wmmmm亚洲av在线观看| 亚洲欧美精品综合久久99| 亚洲精品成人久久久久久| 男人狂女人下面高潮的视频| 日韩欧美精品v在线| 国产亚洲欧美98| av在线观看视频网站免费| 桃色一区二区三区在线观看| av天堂在线播放| 青春草国产在线视频 | 综合色丁香网| 国产伦一二天堂av在线观看| 成人漫画全彩无遮挡| 国产精品99久久久久久久久| 超碰av人人做人人爽久久| 亚洲性久久影院| 偷拍熟女少妇极品色| 久久精品国产亚洲av涩爱 | 内射极品少妇av片p| 青青草视频在线视频观看| 久久精品国产亚洲av涩爱 | 色5月婷婷丁香| 免费看美女性在线毛片视频| 国产成人一区二区在线| 听说在线观看完整版免费高清| 日韩成人伦理影院| 免费看美女性在线毛片视频| 精品少妇黑人巨大在线播放 | 亚洲高清免费不卡视频| 网址你懂的国产日韩在线| 亚洲高清免费不卡视频| 在线a可以看的网站| 一本久久精品| 午夜福利在线观看免费完整高清在 | 又爽又黄a免费视频| 色吧在线观看| 亚洲欧美日韩高清在线视频| 亚洲国产欧洲综合997久久,| 午夜福利视频1000在线观看| 永久网站在线| 亚洲国产精品成人久久小说 | 给我免费播放毛片高清在线观看| 亚洲一区二区三区色噜噜| 亚洲欧美日韩高清专用| 久久久久久九九精品二区国产| 亚洲人与动物交配视频| 免费看a级黄色片| 蜜桃亚洲精品一区二区三区| 色噜噜av男人的天堂激情| 别揉我奶头 嗯啊视频| 99久久中文字幕三级久久日本| 日本一本二区三区精品| 九色成人免费人妻av| kizo精华| 亚洲欧美日韩卡通动漫| 亚洲三级黄色毛片| 91久久精品国产一区二区三区| 97超碰精品成人国产| 超碰av人人做人人爽久久| 黄色一级大片看看| 国产成人影院久久av| 国产精品一区二区性色av| 热99re8久久精品国产| 国产精品福利在线免费观看| 亚洲欧美精品综合久久99| 菩萨蛮人人尽说江南好唐韦庄 | 精品一区二区三区视频在线| 成人毛片a级毛片在线播放| 99热全是精品| 久久精品国产自在天天线| 天堂影院成人在线观看| 亚洲av成人av| 亚洲成av人片在线播放无| 日日干狠狠操夜夜爽| .国产精品久久| 精品欧美国产一区二区三| 日韩精品青青久久久久久| 性欧美人与动物交配| 搡女人真爽免费视频火全软件| 乱码一卡2卡4卡精品| 国产国拍精品亚洲av在线观看| 日韩欧美在线乱码| 午夜亚洲福利在线播放| 国语自产精品视频在线第100页| 一本久久精品| 精品人妻偷拍中文字幕| 欧美成人免费av一区二区三区| 国产精品,欧美在线| 国产伦精品一区二区三区四那| 一级毛片久久久久久久久女|