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

    新疆卡魯安礦區(qū)807號偉晶巖脈近脈圍巖蝕變中成礦元素擴散模型及其影響因素

    2016-06-23 01:14:54趙景宇
    地球化學(xué) 2016年3期
    關(guān)鍵詞:號脈阿爾泰偉晶巖

    陳?陽, 張?輝, 趙景宇

    ?

    新疆卡魯安礦區(qū)807號偉晶巖脈近脈圍巖蝕變中成礦元素擴散模型及其影響因素

    陳?陽1,2, 張?輝1*, 趙景宇1,2

    (1. 中國科學(xué)院 地球化學(xué)研究所 地球內(nèi)部物質(zhì)高溫高壓重點實驗室, 貴州 貴陽?550002; 2. 中國科學(xué)院大學(xué), 北京?100049)

    通過對新疆卡魯安礦區(qū)807號偉晶巖脈體蝕變圍巖樣品的研究, 發(fā)現(xiàn)主元素Fe2O3、Al2O3、K2O、TiO2, 微量元素Li、Rb、Cs、Be等的含量在距接觸帶0.5 m處出現(xiàn)極大值, 向圍巖方向元素含量逐漸減少, 主元素SiO2含量則表現(xiàn)出相反的變化情況; 同時蝕變圍巖中發(fā)育白云母化和電氣石化蝕變, 結(jié)合礦相學(xué)及探針數(shù)據(jù)分析, 認為圍巖中可能發(fā)生3NaAlSi3O8+2H++K+=KAl3Si3O10(OH)2+3Na++6SiO2和2KFe3(AlSi3O10)(OH)2+ SiO2+5Al2SiO5+2Na++6H3BO3= 2NaFe3Al6[Si6O18][BO3]3(OH)4+2K++7H2O等反應(yīng)。對蝕變圍巖中組分含量變化率和距接觸帶的距離的關(guān)系進行擬合, 結(jié)果表明蝕變圍巖中組分的遷移方式以擴散為主, 組分含量變化率是影響擴散距離的主要因素, 有效擴散系數(shù)和擴散時間起次要作用, 同時組分的擴散時間和擴散距離呈正比。從807號脈蝕變圍巖Isocon圖解可看出Li、Rb、Cs等元素具有較高的遷入率, 與前人研究的Li、Rb、Cs具有較高的活動性的結(jié)論相一致。

    稀有金屬; 擴散模型; 圍巖蝕變; 偉晶巖; 阿爾泰

    0?引?言

    新疆阿爾泰地區(qū)是我國乃至世界上最重要的偉晶巖產(chǎn)地之一, 在38個偉晶巖田中, 分布著十萬余條偉晶巖脈, 是我國稀有金屬Li、Be、Rb、Cs等的重要產(chǎn)地[1]。由于長期開采, 阿爾泰地區(qū)已有的稀有金屬礦床均面臨枯竭。近年來, 在哈龍-青河稀有金屬成礦帶上發(fā)現(xiàn)了潛在的大型-超大型卡魯安鋰礦床, 預(yù)示在阿爾泰地區(qū)依然存在著良好的稀有金屬礦床找礦前景。

    由于偉晶巖型稀有金屬礦床無磁性、無導(dǎo)電性以及與偉晶巖造巖礦物之間無明顯的密度差, 對電、磁、重力等地球物理方法無響應(yīng)[2], 高精度 1﹕50000遙感區(qū)測也無法揭露隱伏的偉晶巖盲脈。因此, 尋找隱伏的偉晶巖型稀有金屬礦體一直是難中之難, 在世界范圍內(nèi)均未取得實質(zhì)性突破。

    稀有金屬偉晶巖侵入到變火山巖或變沉積巖圍巖中, 由巖漿分異出的富含稀有元素的流體相與圍巖發(fā)生水-巖相互作用, 形成分散暈, 并以富集高度活動的堿質(zhì)元素(Li、Rb、Cs)、揮發(fā)組分(H2O、B、F)為特征。在稀有金屬偉晶巖周圍的蝕變圍巖中, 通常發(fā)育有鋰閃石化、富Li-Rb-Cs黑云母化、白云母化以及電氣石化蝕變, 其中鋰閃石化的出現(xiàn)可很好地指示稀有金屬偉晶巖隱伏深度不會超過10 m[3]。

    顯然, 利用交代蝕變圍巖的全巖化學(xué)組成特征是預(yù)測偉晶巖盲脈含礦性的有效工具。但截至目前, 尚不清楚與偉晶巖侵入相關(guān)的圍巖蝕變過程中稀有金屬元素遷移方式。在我國阿爾泰, 開展與稀有金屬找礦相關(guān)的圍巖蝕變研究工作非常薄弱, 稀有金屬元素擴散定量建立模型方面還是空白。為此, 本研究擬以新疆阿爾泰卡魯安鋰礦區(qū)稀有金屬礦化偉晶巖807號脈的圍巖蝕變?yōu)檠芯繉ο? 在查明圍巖蝕變的類型、礦物學(xué)及地球化學(xué)特征的基礎(chǔ)上, 建立巖漿流體相與圍巖發(fā)生水-巖相互作用中稀有金屬元素的擴散模型, 揭示擴散模型的影響因素。

    1?地質(zhì)背景

    1.1?區(qū)域地質(zhì)

    新疆阿爾泰造山帶位于西伯利亞板塊西南緣與準(zhǔn)噶爾-哈薩克斯坦板塊結(jié)合部位, 南以額爾齊斯大斷裂為界與準(zhǔn)噶爾板塊相接, 北為西伯利亞板塊, 向南東延伸至蒙古的戈壁阿爾泰, 向北延伸至哈薩克斯坦的礦區(qū)阿爾泰和俄羅斯的山區(qū)阿爾泰。該區(qū)域主要由變質(zhì)變形的沉積巖、火山巖和花崗巖侵入體組成, 根據(jù)新疆阿爾泰造山帶的地層、巖漿巖、變質(zhì)巖及構(gòu)造, 該區(qū)域被劃分成北西-南東走向的4個地塊, 自北向南依次為北阿爾泰、中阿爾泰、瓊庫爾和額爾齊斯構(gòu)造單元(圖1a)[4]。

    北阿爾泰構(gòu)造單元主要出露上泥盆統(tǒng)庫馬蘇組和下石炭統(tǒng)紅山嘴組碎屑沉積巖、角閃巖相至綠片巖相變質(zhì)巖以及少量的黑云母花崗巖深成體; 中阿爾泰構(gòu)造單元出露低變質(zhì)程度的奧陶系哈巴河群砂巖、砂泥巖、泥巖和志留系庫魯姆提群碎屑沉積巖、片巖等, 以及大面積出露的花崗巖。哈巴河群碎屑鋯石206Pb/238U年齡為540~460 Ma,Hf()值主要為正值, 指示了早古生代, 活動大陸邊緣的環(huán)境[5?6]; 瓊庫爾構(gòu)造單元主要由康布鐵堡和阿爾泰建造組成。前者主要含泥盆系弧火山巖[7?8], 后者則主要由濁積巖系列和少量枕狀玄武巖、硅質(zhì)火山巖組成。該地塊也出現(xiàn)了一些含不常見組分的火成巖如埃達克巖、鎂鐵質(zhì)雜巖、高鎂安山巖和玻古安山巖等[9?13], 指示了活動陸緣的環(huán)境[14,15]; 額爾齊斯構(gòu)造單元分布在瓊庫爾構(gòu)造單元和額爾齊斯斷裂之間, 發(fā)育高級片麻巖和片巖, 其上覆蓋泥盆系含化石沉積巖和石炭系火山碎屑巖[7,16], 花崗巖分布規(guī)模小且分散。

    1.2?礦區(qū)地質(zhì)

    福??h卡魯安稀有金屬偉晶巖礦區(qū)處于新疆阿爾泰山中段, 卡拉額爾齊斯河上游地區(qū), 大地構(gòu)造單元屬于阿爾泰褶皺系, 喀納斯-可可托海地槽褶皺帶, 可可托海復(fù)背斜的西北段 (圖1b)。研究區(qū)內(nèi)出露的地層主要為中上志留統(tǒng)庫魯姆提群, 含有兩個亞群: 下亞群主要以云母片巖、石英片巖和石英巖為主, 上亞群主要以變質(zhì)砂巖及少量粉砂巖、石英砂巖為主并含海綿化石。該群與下伏哈巴河群呈斷層接觸, 與上覆康布鐵堡組呈整合接觸關(guān)系。

    卡魯安807號脈位于哈龍河西, 在2010年 1﹕10000地質(zhì)草測時發(fā)現(xiàn)偉晶巖露頭, 長約10 m, 厚1.5~2.7 m, 并在山坡處發(fā)現(xiàn)了近東西向的石英-鈉長石-鋰輝石偉晶巖轉(zhuǎn)石, 經(jīng)探槽及鉆孔揭露偉晶巖脈體長約650 m, 目前已控制礦脈長480 m, 向下延伸較穩(wěn)定, 脈體呈波狀彎曲, 走向近東西, 傾角60°~80° (圖2)。自上而下, 807號脈由鋰云母-石英-鈉長石型偉晶巖過渡到石英-鈉長石-鋰輝石型偉晶巖, 過渡界線在2400 m標(biāo)高左右, 上部含鋰云母-石英-鈉長石偉晶巖呈中粗粒結(jié)構(gòu), 由石英、鈉長石、少量白云母和鋰云母組成, 鋰云母多數(shù)為玫瑰色、少數(shù)淺紅色, 呈鱗片狀, 偶見小巢狀; 下部石英-鈉長石-鋰輝石偉晶巖為中細粒結(jié)構(gòu), 塊狀構(gòu)造, 主要由石英、鈉長石、鋰輝石組成, 脈體中未見綠柱石。圍巖以黑云母-石英片巖為主, 部分圍巖中含紅柱石。與卡魯安807號脈侵入相關(guān)的圍巖普遍發(fā)育白云母化和電氣石化(圖3), 其中白云母化蝕變分布廣泛, 在距接觸帶2.0 m范圍內(nèi)均有分布, 在距接觸帶0.5 m附近發(fā)育較好, 隨距接觸帶距離增加發(fā)育程度逐漸降低; 電氣石化主要發(fā)育在距接觸帶約0.5 m范圍, 隨距接觸帶距離增加發(fā)育程度降低, 在距接觸帶2.0 m處幾乎不可見。

    圖1?新疆阿爾泰造山帶

    (a)卡魯安礦區(qū), (b)地質(zhì)簡圖。數(shù)字1?4分別代表北阿爾泰、中阿爾泰、瓊庫爾和額爾齊斯構(gòu)造單元

    (a) Kalu’an ore area; (b) the geological sketch map of the Kalu’an ore area. Numbers 1?4 represent the North Altai, Central Altai, Qiongkuer and Erqis structural units, respectively

    圖2 卡魯安807號脈體地質(zhì)簡圖

    圖3?與卡魯安807號脈相關(guān)的蝕變圍巖中的白云母化、電氣石化蝕變(距離偉晶巖-圍巖接觸帶0.5 m處)

    2?樣品采集及分析方法

    2.1?樣品采集

    本研究樣品采自卡魯安807號偉晶巖脈侵入的圍巖中上志留統(tǒng)庫魯姆提群蝕變圍巖和未蝕變圍巖,其中蝕變圍巖樣品分別取自近脈體0~2 cm處及距離脈體約0.5、1.0、2.0 m處, 而取自距脈30.0 m處的樣品代表未蝕變的圍巖樣品。

    2.2?分析方法

    蝕變巖石的主要化學(xué)組成分析在中國科學(xué)院地球化學(xué)研究所礦床地球化學(xué)國家重點實驗室完成, 所用儀器是X熒光光譜儀, 儀器型號是AXIOS PANalytical。采用熔片法, 取0.7 g樣品和7 g助熔劑 (Li2B4O7), 攪拌均勻倒入鉑金坩堝中, 再加少量LiBr, 在熔樣爐中加熱至1100 ℃熔融制片。測試結(jié)果的相對標(biāo)準(zhǔn)偏差(RSD)小于1%, 燒失量(LOI)是將1 g樣品加熱到1100 ℃持續(xù)1 h后稱量獲得。

    蝕變巖石的微量元素含量測試在中國科學(xué)院地球化學(xué)研究所礦床地球化學(xué)國家重點實驗室完成。所用儀器是PerkinElmer ELAN DRC-e型電感耦合等離子體質(zhì)譜儀(ICP-MS), 采用HF+HNO3密閉高壓溶樣, 詳細的樣品溶解流程和分析過程見Qi.[15], 分析相對標(biāo)準(zhǔn)偏差(RSD)小于10%。

    蝕變礦物主要化學(xué)組成測試在長安大學(xué)西部礦產(chǎn)資源與地質(zhì)工程教育部重點實驗室完成, 采用JXI-8100型電子探針分析, 加速電壓15 kV, 束電流1.0×10?8A, 束斑直徑1 μm, 分析相對標(biāo)準(zhǔn)偏差(RSD)小于2%。

    3?分析結(jié)果

    3.1?蝕變圍巖巖相學(xué)特征

    圖4顯示黑云母-石英片巖及蝕變巖石背散射(BSE)照片, 利用網(wǎng)格法統(tǒng)計獲得的黑云母-石英片巖及蝕變巖石中主要礦物相相對含量及其電子探針(EMPA)分析結(jié)果列于表1。未蝕變巖石黑云母-石英片巖(樣品807-E)主要由石英、斜長石和黑云母組成, 它們相對體積含量分別約為55%、25%和19%, 不含蝕變礦物白云母和電氣石, 可據(jù)此認為圍巖中發(fā)育的白云母和電氣石均由熱液蝕變形成(表1)。由表1可知, 近脈圍巖 (=0.0 m, 807-A)和弱蝕變圍巖(=2.0 m, 807-D)具有與黑云母-石英片巖相似的礦物組成及相對含量, 石英變化于51%~53%(體積, 下同), 斜長石變化于26%~30%, 黑云母變化于17%~20%之間, 指示由807號偉晶巖脈侵入形成的蝕變帶僅發(fā)育在小于2.0 m范圍, 樣品807-A具有和樣品807-E相似的礦物組成和相對含量說明圍巖具有較好的均一性, 同時我們認為蝕變在807-A樣品中不發(fā)育而在807-B、807-C等樣品中發(fā)育是由于熱液出溶時在圍巖中的遷移速率和熱液組分與圍巖組分間反應(yīng)速率共同作用的結(jié)果。=0.5 m和=1.0 m處蝕變巖石顯示與黑云母-石英片巖顯著不同的礦物組成, 前者(807-B)以含有顯著低的石英(37%, 體積, 下同)、斜長石(21%)、黑云母(15%)、高的白云母(15%)、黑電氣石(9%)以及少量鈦鐵礦(3%)為特征; 后者(807-C)以低的石英(41%)、斜長石(22%)、黑云母(12%)、黑電氣石(7%)、鈦鐵礦(1%)及高的白云母(17%)為特征(表1)。圍巖蝕變前后長石組成未發(fā)生顯著變化, 為奧長石 (根據(jù)計算確定Ab、Or、An變化范圍Ab73.1~78.2Or0.3~0.8An21.2~26.4), 未見鉀長石形成。由此, 圍巖蝕變中可能包含斜長石溶解與白云母的形成, 以及黑云母溶解與黑電氣石的形成(圖5), 其蝕變反應(yīng)可分別寫成:

    3NaAlSi3O8+2H++K+=KAl3Si3O10(OH)2+Na++6SiO2(1)

    2KFe3(AlSi3O10)(OH)2+SiO2+5Al2SiO5+2Na++

    6H3BO3=2NaFe3Al6[Si6O18][BO3]3(OH)4+2K++7H2O (2)

    根據(jù)Shearer.[17]的研究, 上述反應(yīng)中NaCl和H3BO3來自偉晶巖分異的熱液, Al2SiO5來自于長石分解, 黑云母中的Mg、Fe等參與形成電氣石。

    3.2?蝕變圍巖主要化學(xué)組成特征

    未蝕變圍巖及蝕變圍巖的主要化學(xué)組成分析結(jié)果列于表2。由表2可見, 未蝕變巖石黑云母石英片巖(807-E)含有74.71% SiO2、0.51% TiO2、10.17% Al2O3、3.21% Fe2O3、2.61% K2O。與之對照,=0.5 m處蝕變巖石(807-B)具有最低的SiO2含量(63.22%), 最高的TiO2、Al2O3、Fe2O3和K2O含量 (分別為0.85%、15.46%、6.77%和5.46%)。近脈蝕變圍巖(=0.0 m, 807-A)中SiO2、TiO2含量為73.17%、0.44%, 與背景值相差不大; Fe2O3的含量為1.97%, 略低于背景值; 而Al2O3、K2O含量分別為14.25%、4.63%, 明顯高于背景值(圖6)。之后, 隨著距離的逐漸增大 (=1.0 m?30.0 m), 蝕變巖石 (807-C、807-D)中SiO2含量顯示遞增, 而Fe2O3、Al2O3、K2O和TiO2含量顯示遞減的趨勢。

    圖4?卡魯安807號脈體蝕變圍巖背散射照片

    (a) 807-A, (b) 807-B, (c) 807-C, (d) 807-D, (e) 807-E。礦物縮寫見表1

    (a) 807-A, (b) 807-B, (c) 807-C, (d) 807-D, (e) 807-E. Mineral abbreviations are shown in Table 1

    表2?蝕變圍巖中不同位置處(x=0、0.5、1.0和2.0 m)主要化學(xué)組成分析結(jié)果(%)

    注:為蝕變圍巖采樣位置距偉晶巖接觸帶的距離。LOI代表燒失量

    值得注意的是, 上述蝕變圍巖中SiO2、Al2O3、TiO2、Fe2O3和K2O含量隨距離的變化與基于礦物化學(xué)組成、相對含量計算的蝕變巖石中SiO2、Al2O3、TiO2、FeO和K2O含量(表1)隨距離的變化基本上是吻合的, 表明背散射(BSE)下主要礦物相鑒定及相對含量統(tǒng)計是可靠的。

    圖6?卡魯安807號脈體蝕變圍巖中部分主元素含量隨采樣位置變化圖

    3.3?蝕變圍巖稀有金屬元素分布特征

    未蝕變圍巖及蝕變圍巖的主要化學(xué)組成分析結(jié)果列于表3。由表3可知, 未蝕變巖石黑云母石英片巖(807-E)含有Li 137 μg/g、Rb 105 μg/g、Cs 30.6 μg/g和Be 2.4 μg/g。與之對照,= 0.5 m處蝕變巖石(807-B)具有最高含量的Li、Rb、Cs、Be, 分別為798、981、536和61.7 μg/g。近脈蝕變圍巖(= 0.0 m, 807-A)中含有相對低的Li、Rb、Cs、Be含量, 分別為313、841、402和6.8 μg/g。隨著逐漸增大(= 1.0 m?30.0 m), 蝕變巖石 (807-C、807-D)中Li、Rb、Cs、Be元素含量顯示逐漸降低的趨勢, 這種變化與白云母化和電氣石化蝕變在距接觸帶0.5 m處最為發(fā)育, 向圍巖方向發(fā)育程度逐漸降低這一現(xiàn)象相一致, 說明Li、Rb、Cs和Be等元素含量變化明顯受蝕變類型及蝕變強度的影響; 值得注意的是Nb、Ta元素含量從接觸帶向圍巖方向持續(xù)降低, 未受到蝕變類型及強度的影響(圖7)。

    表3?蝕變圍巖中不同位置處(x=0、0.5、1.0和2.0 m)巖石中稀有、稀土元素含量分析結(jié)果(μg/g)

    3.4?蝕變圍巖REE組成特征

    由表3可知, 未蝕變巖石黑云母石英片巖(807-E)以含有?REE (Y) 163.85 μg/g, 弱的Eu異常(Eu=0.57)為特征, 其(La/Sm)N、(La/Gd)N、(La/Yb)N分別為3.75、5.84和9.62, 顯示輕-輕稀土之間、輕-中稀土之間和輕-重稀土之間存在較大的分異; 球粒隕石標(biāo)準(zhǔn)化分布模式顯示為LREE富集的右傾模式(圖8)。由圖8可見,=2.0 m處的弱蝕變圍巖(807-D)以及近脈圍巖(= 0.0 m, 807-A)顯示與未蝕變巖石相似的REE分布模式; 而與之對照,=0.5 m和1.0 m處蝕變巖石(807-B, 807-C)以顯著低的(La/Yb)N(分別為2.96, 3.58)為特征, 指示巖石經(jīng)熱液蝕變后導(dǎo)致HREE相對富集。

    圖7?卡魯安807號脈體蝕變圍巖中部分微量元素含量隨采樣位置變化圖

    圖8?卡魯安807號脈體蝕變圍巖球粒隕石標(biāo)準(zhǔn)化圖解

    4?討?論

    4.1?元素的擴散-滲濾模式建立的理論依據(jù)

    4.1.1?線性等溫下純擴散作用

    在穩(wěn)態(tài)擴散的條件下, 單位時間內(nèi)通過垂直于擴散方向的單位面積的擴散物質(zhì)量 (通稱擴散通量)與該截面處的濃度梯度成正比, 這是菲克第一定律。假設(shè)擴散沿軸方向進行, 其表達式為:

    為擴散通量(mol/m2·s),為擴散系數(shù)(m2/s), dd為濃度梯度(mol/m3·m), “?”表示擴散方向為濃度梯度的反方向, 即擴散由高濃度向低濃度區(qū)進行;

    在非穩(wěn)態(tài)擴散條件下, 在距離處, 濃度隨時間的變化率等于該處的擴散通量隨距離變化率的負值 (菲克第二定律), 其表達式為:

    根據(jù)熱力學(xué)理論, 認為擴散過程與其他物理化學(xué)過程一樣, 其發(fā)生的根本推動力應(yīng)該是化學(xué)位梯度 (如電場、磁場等)。由此, 在其作用下質(zhì)點熱運動導(dǎo)致的定向遷移稱為擴散; 當(dāng)化學(xué)位梯度為零時, 擴散達到平衡。

    若巖石的孔隙流體處于靜止?fàn)顟B(tài), 體系發(fā)生純擴散作用, 可認為多孔介質(zhì)中的流體是服從亨利定律的無限稀溶液。純擴散作用則可用一維擴散作用菲克第二定律的標(biāo)準(zhǔn)形式來表示, 如式(4)所示。

    4.1.2?有效擴散系數(shù)推導(dǎo)

    為氣體常數(shù),為絕對溫度,為離子電價,為法拉第常數(shù)。

    對于在熱液中以絡(luò)合物形式遷移的元素, 根據(jù)Kohlrausch 離子獨立移動定律, 其極限摩爾電導(dǎo)率可以認為是組成絡(luò)合物的離子極限摩爾電導(dǎo)率之和。對于偉晶巖體系而言, 晚期出溶的流體相中揮發(fā)分F、B、P對于稀有金屬的遷移起著重要的作用, 尤其是F[18], 因此本文選取F作為主要絡(luò)合劑, 對形成的LiF、RbF、CsF絡(luò)合物的遷移狀態(tài)進行研究。

    極限摩爾電導(dǎo)率大小受溫度變化影響十分明顯, 一般具有以下表達式[19]:

    =0.02K?1, 代入之后可以得到極限摩爾電導(dǎo)率隨溫度變化的計算表達式。

    根據(jù)Arrhenius 方程建立有效擴散系數(shù)、溫度和擴散激活能()之間的相互關(guān)系, 可得到離子或絡(luò)合物在熱液中擴散系數(shù)隨溫度變化的經(jīng)驗公式:

    據(jù)Morgan.[20]對Tanco偉晶巖、朱金初等[21]對可可托海3號偉晶巖脈的研究, 推測流體相出溶及其隨之發(fā)生的圍巖蝕變反應(yīng)在500~550 ℃溫度范圍, 即=773 ~ 823 K。

    在地質(zhì)過程中, 熱液在地質(zhì)體中有效擴散系數(shù)主要受到介質(zhì)孔隙度和曲折因子的影響[22], 可表述為:

    為巖石的孔隙度,為巖石的曲折因子。由于曲折因子無法直接測量, 擴散率(′)常被表示為孔隙度的經(jīng)驗函數(shù)[23]:

    為經(jīng)驗指數(shù), 取值在2.0~2.5之間, 考慮到研究區(qū)的巖性本文取=2.0, 整理后獲得有效擴散系數(shù)為:

    利用該式, 所計算的絡(luò)合物在巖石中的擴散系數(shù)列于表4。

    表4?計算的LiF、RbF和CsF在黑云母-石英片巖中的擴散系數(shù)

    4.1.3?線性等溫下純滲濾作用

    除了在化學(xué)位梯度作用下發(fā)生擴散外, 元素在壓力梯度下發(fā)生的滲濾作用, 滲濾流可用下式表示[26]:

    I=(11)

    I為滲濾流,為流體質(zhì)點的平均移動速度,為介質(zhì)巖石的孔隙度,為組分濃度。

    4.1.4?擴散和滲濾的共同作用

    在實際地質(zhì)作用過程中, 流體的運移往往是擴散和滲濾共同作用以及其他化學(xué)反應(yīng)重疊的結(jié)果, 若假設(shè)孔隙流體是理想溶液, 擴散和滲濾作用的數(shù)學(xué)表達式如下:

    與邊界、初始條件相結(jié)合, 可組成一維擴散作用的微分方程組:

    解這個方程組得到:

    0和分別為熱液中離子或絡(luò)合物在扣除背景值后的初始含量及在位置處的實際含量,為擴散距離,為擴散時間,eff為有效擴散系數(shù), 并假定其在擴散過程中不發(fā)生變化。

    在利用野外數(shù)據(jù)對模型進行驗證過程中我們發(fā)現(xiàn)式(15)并不滿足野外數(shù)據(jù), 考慮到卡魯安807號脈體圍巖的孔隙度約為0.093, 在孔隙度較小時我們可忽略流體的滲濾作用, 對式(15)進行簡化后得到流體擴散作用的表達式:

    4.2?與卡魯安807號脈侵入相關(guān)的圍巖蝕變過程中元素遷移及方式

    4.2.1?圍巖蝕變中元素遷移方式-擴散模型

    由式(16)可知, 隨擴散進行稀有金屬元素濃度()變化取決于其初始濃度(0)、有效擴散系數(shù)(eff)、擴散距離以及擴散時間()。假定eff在擴散過程中不發(fā)生變化, 余補誤差函數(shù)erfc在區(qū)間內(nèi)為減函數(shù), 擴散距離與為負相關(guān)關(guān)系, 若隨擴散距離增加,值減少較慢則擴散較遠, 即值變化率是影響擴散距離的最重要因素。對于本次研究的圍巖蝕變, 由于距離偉晶巖接觸帶= 0 m處Li、Rb、Cs、Be異常低的分布(圖7), 無法獲得其準(zhǔn)確的初始濃度值(0)。但由于0為定值, 因此我們可以獲得0是有效擴散系數(shù)(eff)、擴散距離()以及擴散時間()的函數(shù)關(guān)系。將距卡魯安807號脈與圍巖接觸帶= 0.5、0、2.0和30.0 m處實際測得的全巖中Li、Rb、Cs含量及估算的接觸帶處初始的Li、Rb、Cs含量代入式(16), 擬合出0與擴散距離()的函數(shù)關(guān)系(圖9)。

    如圖9所示, LiF、RbF、CsF在卡魯安807號脈周圍的蝕變圍巖中的濃度變化基本上符合擴散的特征。從圖9中也可看出LiF在蝕變圍巖中的遷移距離最遠, CsF和RbF次之, 對應(yīng)的遷移時間分別為(1.73 ± 0.31)×103a, (0.85 ± 0.15)×103a和(0.80±0.15)× 103a, 遷移時間和遷移距離呈正比同時較好的擬合曲線也說明圍巖僅經(jīng)歷了一次熱液事件。

    4.2.2?圍巖蝕變中組分(元素)的遷移

    圍巖蝕變過程中不僅包含組分帶入和帶出, 而且不可避免涉及蝕變前后體積的變化。Gresens[27]最初提出了巖石交代蝕變的成分-體積分析法, 利用穩(wěn)定元素的質(zhì)量計算其他元素質(zhì)量的遷移。在此基礎(chǔ)上, Grant[28]提出的Isocon圖解法(或等濃度線法), 可不考慮圍巖蝕變過程中巖石體積的變化, 可直觀看出元素遷入與遷出。

    通常認為Al、Ti、Zr、Y、Nb組分在水-巖相互作用中穩(wěn)定[29]。通過對807號脈周圍蝕變圍巖主要化學(xué)組成及微量元素組成分析, 發(fā)現(xiàn)Y在蝕變過程中變化很小, 選其為不活動元素, 并取距接觸帶30.0 m處樣品中數(shù)據(jù)作為橫坐標(biāo)(背景值), 距接觸帶0.5 m處樣品中數(shù)據(jù)作為縱坐標(biāo), 對數(shù)據(jù)進行一定程度的等比例縮放。作出= 0.5 m處的蝕變圍巖的Isocon圖解。由圖10可見, 圍巖遭受熱液蝕變過程中, Li、Rb、Cs、Be、Ta顯示顯著的遷入, Na、Zr顯示明顯的帶出。此外, Al、Fe、K、P、Nb、Ti顯示少量遷入, Si、Mg、Ca顯示少量帶入。

    0代表扣除背景值后元素含量與初始值的比值

    0represents the ratio of element content after deducting the background content over its primary content

    顯然, 蝕變圍巖(807-B、807-C)中Li、Rb、Cs、Be、Ta含量顯著增大, 暗示偉晶巖演化晚期出溶的流體相中含有較高的Li、Rb、Cs、Be、Ta含量。這與Shearer.[17]研究South Dakota的Black Hills附近3個偉晶巖體Etta、Bob Ingersoll No. 1、Peerless圍巖蝕變, 發(fā)現(xiàn)石英-云母片巖中發(fā)育Li、Rb、Cs、F、B、As、U等元素的擴散暈 (<30.0 m)以及顯示Li > Cs≥Rb活動性是一致的。

    圖10?卡魯安807號脈體蝕變圍巖等濃度線圖解法

    o為未蝕變樣品或蝕變程度最低樣品組分,A為距接觸帶0.5 m內(nèi)蝕變樣品組分。氧化物以質(zhì)量百分比(%), 元素以濃度(μg/g)為單位進行投點, 圖中直線代表由元素Y所確定的Isocon線, 氧化物及元素前的數(shù)字代表比例系數(shù)

    whereorepresents the compositions of fresh or least altered samples,Arepresents the compositions of the samples collected at 0.5 m away from the contact zone in the country rocks. Oxides are plotted in weight percent and elements are expressed in per million. The straight line indicates the unified Isocon defined by Y; the numbers before the oxide and element symbols represent the scaling coefficients

    5?結(jié)?論

    (1) 通過對卡魯安807號脈體蝕變圍巖中所采集樣品中礦物及元素含量的分析發(fā)現(xiàn), 距離接觸帶0.5 m處樣品中發(fā)育有較強的白云母化、電氣石化蝕變伴隨著較低的石英、斜長石、黑云母含量, 出現(xiàn)主元素Fe2O3、Al2O3、K2O、TiO2含量的極大值, SiO2含量的極小值及微量元素Li、Rb、Cs、Be含量的極大值, 向圍巖方向除SiO2外, 以上組分的含量持續(xù)降低, HREE在蝕變發(fā)育較強烈的樣品中相對富集。

    (2) 通過對卡魯安807號脈體蝕變圍巖中組分遷移模式的分析, 發(fā)現(xiàn)組分的遷移方式以擴散為主, 滲濾作用較小, 影響組分擴散距離的主要因素是組分含量變化率, 組分的遷移時間和遷移距離呈正比, 同時較好的擬合曲線也說明圍巖僅經(jīng)歷了一次熱液事件; 根據(jù)蝕變圍巖中距接觸帶0.5 m處樣品中Li、Rb、Cs含量及背景值, 從Isocon圖解中可看出Li、Rb、Cs在圍巖蝕變過程中具有較高的遷入率, 與前人的研究認為Li、Rb、Cs具有較高的活動性的結(jié)論相一致。

    [1] 鄒天人, 李慶昌. 中國新疆稀有及稀土金屬礦床[M]. 北京: 地質(zhì)出版社, 2006: 1–284. Zou Tian-ren, Li Qing-chang. Rare and Rare Earth Metal Deposits in Xinjiang, China[M]. Beijing: Geological Publishing House, 2006: 1–284 (in Chinese).

    [2] Milkereit B. Ore deposits and exploration technology[M]// Galeschuk C, Vanstone P. Exploration Techniques for Rare-Element Pegmatite in the Bird River Greenstone Belt, Southeastern Manitoba. Canada: Proceedings of Exploration 07, 2007: 823–839.

    [3] London D. Holmquistite as a guide to pegmatitic rare metal deposits[J]. Econ Geol, 1986, 81(3): 704–712.

    [4] Cai Keda, Sun Min, Yuan Chao, Xiao Wenjiao, Zhao Guoshun, Long Xiaoping, Wu Fuyuan. Carboniferous mantle-derived felsic intrusion in the Chinese Altai, NW China: Implications for geodynamic change of the accretionary orogenic belt[J]. Gondwana Res, 2012, 22(2): 681–698.

    [5] Long Xiaoping, Sun Min, Yuan Chao, Xiao Wenjiao, Lin Shoufa, Wu Fuyuan, Xia Xiaoping, Cai Keda. Detrital zircon age and Hf isotopic studies for metasedimentary rocks from theChinese Altai: Implications for the Early Paleozoic tectonic evolution of the Central Asian Orogenic Belt[J]. Tectonics, 2007, 26(5): 1–20.

    [6] Long Xiaoping, Yuan Chao, Sun Min, Xiao Wenjiao, Zhao Guochun, Wang Yujing, Cai Keda, Xia Xiaoping, Xie Liewen. Detrital zircon ages and Hf isotopes of the early Paleozoic flysch sequence in the Chinese Altai, NW China: New constrains on depositionalage, provenance and tectonic evolution[J]. Tectonophysics, 2010, 480(1–4): 213–231.

    [7] Windley B F, Kroner A, Guo Jinghui, Qu Guosheng, Li Yingyi, Zhang Chi. Neoproterozoic to Paleozoic geology of the Altay orogen, NW China: New zircon age data and tectonic Evolution[J]. Geology, 2002, 110(6): 719–737.

    [8] Chai Fengmei, Mao Jingwen, Dong Lianhui, Yang Fuquan, Liu Feng, Geng Xinxia, Zhang Zhixin. Geochronology of metarhyolites from the Kangbutiebao Formation in the Kelang basin, Altay Mountains, Xinjiang: Implications for the tectonic evolution and metallogeny[J]. Gondwana Res, 2009, 16(2): 189–200.

    [9] Niu Hecai, Sato H, Zhang Haixiang, Ito J, Yu Xueyuan, Nagao T, Terada K, Zhang Xueyuan. Juxtaposition of adakite, boninite, high-TiO2and low-TiO2basalts in theDevonian southern Altay, Xinjiang, NW China[J]. J Asian Earth Sci, 2006, 28(4–6): 439–456.

    [10] Niu Hecai, Xu Jifeng, Yu Xueyuan, Chen Fanrong, Zheng Zouping. Discovery of Mg-rich volcanic rock series in western Altay area, Xinjiang and its geologic significance[J]. Chinese Sci Bull, 1999, 44(18): 1685–1688.

    [11] Xu Jifeng, Castillo P R, Chen Fanrong, Niu Hecai, Yu Xueyuan, Zheng Zuoping. Geochemistry of late paleozoic mafic igneous rocks from the Kuerti area, Xinjiang, northwest China: Implications for backarc mantle evolution[J]. Chem Geol, 2003, 193(1/2): 137–154.

    [12] Cai Keda, Sun Min, Yuan Chao, Zhao Guochun, Xiao Wenjiao, Long Xiaoping, Wu Fuyuan. Geochronological and geochemical study of mafic dykes from the northwest Chinese Altai: Implications for petrogenesis and tectonic evolution[J]. Gondwana Res, 2010, 18(4): 638–652.

    [13] Wong Kenny, Sun Min, Zhao Guochun, Yuan Chao, Xiao Wenjiao. Geochemical and geochronological studies of the Alegedayi ophiolitic complex and its implication for the evolution of the Chinese Altai[J]. Gondwana Res, 2010, 18(2/3): 438–454.

    [14] Yuan Chao, Sun Min, Xiao Wenjiao, Li Xianhua, Chen Hanlin, Lin Shoufa, Xia Xiaoping, Long Xiaoping. Accretionary orogenesis of the Chinese Altai: Insights from Paleozoic Granitoids[J]. Chem Geol, 2007, 242(1/2): 22–39.

    [15] Qi Liang, Hu Jing, Gregoire D C. Determination of trace elements in granites by inductively coupled plasma mass spectrometry[J]. Talanta, 2000, 51(3): 507–513.

    [16] 曲國勝, 崇美英. 阿爾泰造山帶的鉛同位素地質(zhì)及其構(gòu)造意義[J]. 現(xiàn)代地質(zhì), 1991, 5(1): 100–110. Qu Guo-sheng, Chong Mei-ying. Lead isotope geology and its tectonic implication in Altaids, China[J]. Geoscience, 1991, 5(1): 100–110 (in Chinese with English abstract).

    [17] Shearer C K, Papike J J, Simon S B. Pegmatite/wallrock interactions, Black Hills, South Dakota: Progressive boron metasomatism adjacent to the Tip Top pegmatite[J]. Am Mineral, 1986, 48(12): 518–539.

    [18] Linnen R L, van Lichtervelde M, ?erny P. Granitic pegmatites as sources of strategic metals[J]. Elements, 2012, 8(4): 275–280.

    [19] Coury L. Conductance measurements Part 1: Theory[J]. Curr Sep, 1999, 18(3): 91–96.

    [20] Morgan G B VI, London D. Alteration of amphibolitic wallrocks around the Tanco rare-element pegmatite, Bernic Lake, Manitoba[J]. Am Mineral, 1987, 72(11/12): 1097–1121.

    [21] 朱金初, 吳長年, 劉昌實, 李福春, 黃小龍, 周東山. 新疆阿爾泰可可托海3號偉晶巖脈巖漿—熱液演化和成因[J]. 高校地質(zhì)學(xué)報, 2000, 6(1): 46–47. Zhu Jin-chu, Wu Chang-nian, Liu Chang-shi, Li Fu-chun, Huang Xiao-long, Zhou Dong-shan. Magmatic-Hydrothermal evolution and genesis of Koktokay No. 3 rare metal pegmatite dyke, Altai, China[J]. Geol J China Univ, 2000, 6(1): 46–47 (in Chinese with English abstract).

    [22] Liu Jianguo, Nie Yongfeng. Fractal scaling of effective diffusion coefficient of solute inporous Media[J]. J Environ Sci, 2001, 13(2): 170–172.

    [23] Boving T B, Grathwohl P. Trace diffusion coefficients in sedimentary rocks correlationto porosity and hydraulic conductivity[J]. J Contam Hydrol, 2001, 53(1): 85–100.

    [24] Ejima T, Sato Y, Yaegashi S, Kijima T, Takeuchi E, Tamai K. Visosity of molten alkali fluorides[J]. J Japan Inst Met, 1987, 51(4): 328–337.

    [25] Marciniszyn T. Unusual application of porous matrix made from quartzite schist[J]. Physicochem Probl Miner Process, 2013, 49(1): 323–328.

    [26] Fletcher RC, Hojmann A W. Simple models of diffusion and combined diffusion–infiltration metasomatism[M]//Hofmann A W, Giletti B J, Yoder H S Jr, Yund R A. Geochemical Transport and Kinetics. Washington: Carnegie Institution of Washington Publication, 1974, 634: 243–259.

    [27] Gresens R L. Composition-volume relationships of metasomatism[J]. Chem Geol, 1967, 2(1): 47–65.

    [28] Grant J A. The Isocon Diagram — A simple solution to Gresens' Equation for metasomatic Alteration[J]. Econ Geol, 1986, 81(8): 1976–1982.

    [29] Maclean W H, Kranidiotis P. Immobile elements as monitors of mass transfer in hydrothermal alteration: Phelps Dodge massive sulfide adeposit, Matagami, Quebec[J]. Econ Geol, 1987, 82(4): 951–962.

    Altered country rocks of No. 807 pegmatite vein in the Kalu’an ore area, Xinjiang: Ore-forming element diffusion model and its influencing factors

    CHEN Yang1,2, ZHANG Hui1*and ZHAO Jing-yu1,2

    1. Laboratory for High Temperature and High Pressure Study of the Earth’s Interior, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang?550002, China; 2. University of Chinese Academy of Sciences, Beijing?100049, China

    Based on our studies on altered country rock samples from No.807 pegmatite vein in Kalu’an ore area in Xinjiang, it is found that major element oxides such as Fe2O3, Al2O3, K2O, TiO2and trace elements such as Li, Rb, Cs, Be reach their maximum contents at approximately 0.5m away from the contact zone between the pegmatite vein and the country rocks, and tend to decrease in the direction from the contact zone to the country rocks, while the content of SiO2shows an opposite tendency of variation. On the basis of the available mineralogical and electron probe analytical data, the occurrence of muscovite and tourmaline in the altered rocks probably indicates the following reactions, e.g., 3NaAlSi3O8+2H++K+=KAl3Si3O10(OH)2+3Na++6SiO2and 2KFe3(AlSi3O10)(OH)2+ SiO2+5Al2SiO5+2Na++6H3BO3=2NaFe3Al6[Si6O18][BO3]3(OH)4+2K++7H2O. By employing the data fitting of variation rates of component contents and the corresponding distance away from the contact zone, it is suggested that the component migration patterns are controlled mainly by diffusion. Among the factors that influence the diffusion distance, the variation rate for component contents is the primary factor, while both effective diffusion coefficient and diffusion time are the secondary factors. Meanwhile, we found have that the diffusion time is proportional to the diffusion distance. The Isocon diagram for the altered country rocks of No.807 pegmatite vein shows that elements such as Li, Rb and Cs exhibit higher migration rates, which is consistent with the results of the previous study that Li, Rb and Cs exhibit much higher mobility.

    rare metal; diffusion model; country rock alteration; pegmatite; Altai

    P595

    A

    0379-1726(2016)03-0268-13

    2015-03-19;

    2015-06-10;

    2015-07-28

    國家自然科學(xué)基金(41372104); 新疆有色金屬工業(yè)集團地質(zhì)科研項目(YSKY2011-02)

    陳陽(1989–), 男, 碩士研究生, 巖石地球化學(xué)研究方向。E-mail: desesile@126.com

    ZHANG Hui, E-mail: zhanghui@vip.gyig.ac.cn; Tel: +86-851-85891454

    猜你喜歡
    號脈阿爾泰偉晶巖
    剛果(金)馬諾諾(Manono)東部某偉晶巖型鈮鉭礦地質(zhì)特征
    媽媽的吻
    湘東北幕阜山巖體南部稀有金屬偉晶巖分帶特征研究
    盧旺達Gatumba地區(qū)花崗偉晶巖的地質(zhì)、地球化學(xué)特征及其成因研究綜述
    湖南省平江縣萬古礦區(qū)十號脈成礦規(guī)律及找礦方向分析
    媽媽的吻
    做人與處世(2020年3期)2020-03-07 05:23:55
    山東玲瓏金礦床地質(zhì)特征及深部找礦潛力分析
    新疆阿爾泰鐵礦成礦規(guī)律淺析
    給那些愛 “病”的人搭個脈
    婦女之友(2015年5期)2015-09-08 07:27:43
    阿爾泰發(fā)現(xiàn)大規(guī)模巖畫群
    404 Not Found

    404 Not Found


    nginx
    99久久中文字幕三级久久日本| 亚洲国产成人一精品久久久| 国产成人精品婷婷| 免费女性裸体啪啪无遮挡网站| 乱码一卡2卡4卡精品| 日韩欧美精品免费久久| 亚洲欧洲国产日韩| 久久女婷五月综合色啪小说| 国产亚洲最大av| 国产成人精品在线电影| 一级爰片在线观看| 亚洲欧美清纯卡通| 欧美精品国产亚洲| 亚洲激情五月婷婷啪啪| 秋霞伦理黄片| 国产精品一二三区在线看| 校园人妻丝袜中文字幕| 免费播放大片免费观看视频在线观看| 免费大片黄手机在线观看| 亚洲av在线观看美女高潮| 99久久综合免费| 黑人高潮一二区| a级毛片在线看网站| 久久精品国产亚洲av涩爱| 熟妇人妻不卡中文字幕| 婷婷色综合www| 人人澡人人妻人| 国产综合精华液| 丝袜美足系列| 午夜福利视频精品| 国产成人a∨麻豆精品| 日日摸夜夜添夜夜爱| 日本wwww免费看| 全区人妻精品视频| 久久久欧美国产精品| 如何舔出高潮| 久久国内精品自在自线图片| 国产精品欧美亚洲77777| av网站免费在线观看视频| 亚洲在久久综合| 欧美成人午夜免费资源| 99久国产av精品国产电影| 中国国产av一级| 日韩不卡一区二区三区视频在线| 亚洲国产精品一区二区三区在线| 18+在线观看网站| 久久久久久久大尺度免费视频| 亚洲国产最新在线播放| 有码 亚洲区| 波多野结衣一区麻豆| 精品国产露脸久久av麻豆| 国产成人精品一,二区| 青春草国产在线视频| 一区二区三区四区激情视频| 一级毛片 在线播放| 国产永久视频网站| 中国美白少妇内射xxxbb| 久久精品国产亚洲av天美| 欧美国产精品va在线观看不卡| 亚洲欧洲精品一区二区精品久久久 | 黄色 视频免费看| 一区二区三区乱码不卡18| 国国产精品蜜臀av免费| 少妇被粗大猛烈的视频| 亚洲国产精品一区二区三区在线| 欧美激情国产日韩精品一区| 精品人妻一区二区三区麻豆| 少妇熟女欧美另类| 亚洲婷婷狠狠爱综合网| 国产视频首页在线观看| 精品一区二区免费观看| a级毛片在线看网站| 夫妻性生交免费视频一级片| 赤兔流量卡办理| 男女午夜视频在线观看 | 亚洲欧美一区二区三区黑人 | 久久这里有精品视频免费| 成年av动漫网址| 熟女人妻精品中文字幕| 国产1区2区3区精品| 日韩av免费高清视频| 欧美日韩精品成人综合77777| 日本黄大片高清| 日韩,欧美,国产一区二区三区| 男人添女人高潮全过程视频| 国产成人欧美| 九色亚洲精品在线播放| 亚洲成人手机| 欧美另类一区| 七月丁香在线播放| 亚洲美女黄色视频免费看| 免费人妻精品一区二区三区视频| 欧美激情 高清一区二区三区| 91成人精品电影| 国产精品蜜桃在线观看| 亚洲人与动物交配视频| 黑人欧美特级aaaaaa片| 亚洲成国产人片在线观看| 两性夫妻黄色片 | 久久99蜜桃精品久久| 日本黄大片高清| 国产成人精品在线电影| 伦精品一区二区三区| 波野结衣二区三区在线| 国产精品成人在线| 一级片'在线观看视频| 午夜老司机福利剧场| 国产精品国产三级国产专区5o| 一级毛片我不卡| 国国产精品蜜臀av免费| 色吧在线观看| 国产一区二区三区av在线| 国产福利在线免费观看视频| 免费av不卡在线播放| 免费在线观看完整版高清| 午夜福利乱码中文字幕| 久久久国产欧美日韩av| 精品99又大又爽又粗少妇毛片| 日韩 亚洲 欧美在线| 国产精品久久久久久av不卡| 搡老乐熟女国产| 777米奇影视久久| 日本wwww免费看| 美女主播在线视频| 国产黄色免费在线视频| 亚洲av日韩在线播放| 免费av不卡在线播放| 在线天堂中文资源库| 18+在线观看网站| 秋霞伦理黄片| 美女中出高潮动态图| 亚洲少妇的诱惑av| 晚上一个人看的免费电影| 秋霞在线观看毛片| 熟女人妻精品中文字幕| 免费日韩欧美在线观看| 精品久久久精品久久久| 99久久综合免费| 啦啦啦视频在线资源免费观看| 9色porny在线观看| 国产精品欧美亚洲77777| 大码成人一级视频| 亚洲国产精品专区欧美| 制服诱惑二区| 日本与韩国留学比较| 男女边摸边吃奶| 久久99一区二区三区| 亚洲第一区二区三区不卡| 69精品国产乱码久久久| 国产成人精品无人区| 国产av精品麻豆| 国产免费现黄频在线看| 欧美日本中文国产一区发布| 久久久久久久久久成人| 免费看不卡的av| 国产午夜精品一二区理论片| 成年动漫av网址| 精品一品国产午夜福利视频| 精品一区二区三卡| 国产高清三级在线| 午夜福利视频在线观看免费| 亚洲国产精品国产精品| 亚洲精品一区蜜桃| 三上悠亚av全集在线观看| 久久精品久久久久久久性| 三级国产精品片| 韩国av在线不卡| 日韩 亚洲 欧美在线| 熟女人妻精品中文字幕| 欧美日韩av久久| 亚洲国产色片| 色婷婷久久久亚洲欧美| 亚洲av免费高清在线观看| 午夜福利影视在线免费观看| 精品久久久精品久久久| av电影中文网址| 2021少妇久久久久久久久久久| 亚洲四区av| 曰老女人黄片| 国产黄色视频一区二区在线观看| 国产黄色视频一区二区在线观看| 日韩 亚洲 欧美在线| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区国产| 在线免费观看不下载黄p国产| av女优亚洲男人天堂| 老司机影院成人| 在线观看免费日韩欧美大片| 插逼视频在线观看| 午夜福利在线观看免费完整高清在| 寂寞人妻少妇视频99o| 国产女主播在线喷水免费视频网站| 精品福利永久在线观看| 国产精品国产三级国产av玫瑰| 国产成人av激情在线播放| 午夜精品国产一区二区电影| 五月开心婷婷网| 国产精品无大码| 一区在线观看完整版| 最新中文字幕久久久久| 日韩制服骚丝袜av| 成人国语在线视频| 国产欧美另类精品又又久久亚洲欧美| 777米奇影视久久| 黑人猛操日本美女一级片| 成人18禁高潮啪啪吃奶动态图| 国产成人av激情在线播放| 亚洲,一卡二卡三卡| 亚洲国产最新在线播放| 欧美国产精品va在线观看不卡| 人妻少妇偷人精品九色| 久久久久国产精品人妻一区二区| 欧美xxxx性猛交bbbb| 久久精品久久久久久久性| 欧美日韩亚洲高清精品| 纯流量卡能插随身wifi吗| 久久精品国产亚洲av天美| 又黄又爽又刺激的免费视频.| 亚洲三级黄色毛片| 蜜桃在线观看..| 看免费成人av毛片| 99国产综合亚洲精品| 欧美精品一区二区大全| 亚洲精品,欧美精品| 欧美bdsm另类| 菩萨蛮人人尽说江南好唐韦庄| 大码成人一级视频| 亚洲精品国产色婷婷电影| 亚洲精品久久午夜乱码| videossex国产| 国产精品嫩草影院av在线观看| 男人爽女人下面视频在线观看| 水蜜桃什么品种好| 人妻系列 视频| 18禁动态无遮挡网站| 午夜福利视频在线观看免费| 久久婷婷青草| 不卡视频在线观看欧美| 国产熟女欧美一区二区| 午夜福利网站1000一区二区三区| 亚洲精品乱码久久久久久按摩| 国产精品成人在线| 亚洲国产精品专区欧美| 一级a做视频免费观看| 午夜免费鲁丝| 亚洲综合精品二区| 男人爽女人下面视频在线观看| 亚洲av综合色区一区| 久久精品国产a三级三级三级| 成人综合一区亚洲| 日韩成人伦理影院| 久久国产精品大桥未久av| 美女大奶头黄色视频| 观看av在线不卡| 一级毛片 在线播放| 亚洲欧洲精品一区二区精品久久久 | 成人国产av品久久久| 日本vs欧美在线观看视频| 国产精品三级大全| 午夜免费鲁丝| 国产亚洲精品久久久com| 麻豆精品久久久久久蜜桃| 久久久久久人妻| 在线观看一区二区三区激情| 国产成人一区二区在线| 美女脱内裤让男人舔精品视频| 1024视频免费在线观看| 国产国语露脸激情在线看| 高清av免费在线| 亚洲精品av麻豆狂野| a级毛片黄视频| 中文字幕亚洲精品专区| tube8黄色片| 伦精品一区二区三区| 中文字幕人妻丝袜制服| 18在线观看网站| 国产成人免费观看mmmm| av天堂久久9| 精品国产国语对白av| 久久人人爽人人片av| 97在线人人人人妻| 欧美老熟妇乱子伦牲交| 插逼视频在线观看| 国产精品成人在线| av又黄又爽大尺度在线免费看| 午夜免费鲁丝| 啦啦啦啦在线视频资源| 日韩,欧美,国产一区二区三区| 欧美bdsm另类| 中文字幕av电影在线播放| 老司机亚洲免费影院| 超碰97精品在线观看| 国产xxxxx性猛交| 最近中文字幕2019免费版| 免费人妻精品一区二区三区视频| 日韩大片免费观看网站| 亚洲情色 制服丝袜| 亚洲欧美清纯卡通| 欧美3d第一页| 青春草视频在线免费观看| 看十八女毛片水多多多| 亚洲精品视频女| www.色视频.com| 在线观看美女被高潮喷水网站| 五月天丁香电影| 九九爱精品视频在线观看| 少妇被粗大的猛进出69影院 | 欧美日韩精品成人综合77777| 国产精品免费大片| 美女福利国产在线| 韩国高清视频一区二区三区| 国产色婷婷99| 国产男女超爽视频在线观看| 亚洲美女黄色视频免费看| 高清av免费在线| 欧美日韩国产mv在线观看视频| 免费看av在线观看网站| 日韩电影二区| 欧美日韩综合久久久久久| 一区二区三区精品91| 国产国拍精品亚洲av在线观看| 精品国产国语对白av| 久热久热在线精品观看| 午夜免费男女啪啪视频观看| 日本-黄色视频高清免费观看| 日韩精品有码人妻一区| 青春草亚洲视频在线观看| 国内精品宾馆在线| av线在线观看网站| 捣出白浆h1v1| 18在线观看网站| 高清av免费在线| 亚洲欧美清纯卡通| 久久人人97超碰香蕉20202| 宅男免费午夜| 久久久久久久久久久久大奶| 99热网站在线观看| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡 | 亚洲av成人精品一二三区| 久久人人爽av亚洲精品天堂| 97人妻天天添夜夜摸| 在线观看一区二区三区激情| 久久久久久久久久人人人人人人| 国产熟女午夜一区二区三区| 成人影院久久| 国产 精品1| 青春草视频在线免费观看| 熟女电影av网| 日韩伦理黄色片| 亚洲av电影在线观看一区二区三区| 久久久久精品性色| 久久国产精品男人的天堂亚洲 | 亚洲色图综合在线观看| 性色av一级| 两个人免费观看高清视频| 国产成人精品无人区| 各种免费的搞黄视频| 黄色毛片三级朝国网站| 欧美丝袜亚洲另类| 97超碰精品成人国产| a级毛片在线看网站| 有码 亚洲区| 国产精品久久久久久精品电影小说| 国产精品麻豆人妻色哟哟久久| 亚洲经典国产精华液单| 亚洲天堂av无毛| 精品第一国产精品| 麻豆精品久久久久久蜜桃| 18+在线观看网站| 午夜日本视频在线| 男女国产视频网站| 最近的中文字幕免费完整| 一区二区日韩欧美中文字幕 | 日本wwww免费看| 91精品伊人久久大香线蕉| 国产爽快片一区二区三区| 女人被躁到高潮嗷嗷叫费观| 日日爽夜夜爽网站| 婷婷色av中文字幕| 亚洲欧美成人精品一区二区| 婷婷色麻豆天堂久久| 国产激情久久老熟女| 国产亚洲精品久久久com| 91成人精品电影| 美女脱内裤让男人舔精品视频| 少妇人妻久久综合中文| 免费观看无遮挡的男女| 在线亚洲精品国产二区图片欧美| 国产欧美亚洲国产| 99热国产这里只有精品6| 日本爱情动作片www.在线观看| 久热久热在线精品观看| www.av在线官网国产| 久久久久精品久久久久真实原创| 国产免费视频播放在线视频| 国产爽快片一区二区三区| 国产免费现黄频在线看| 菩萨蛮人人尽说江南好唐韦庄| 午夜91福利影院| av在线老鸭窝| 亚洲欧美精品自产自拍| 久久人人爽人人爽人人片va| a 毛片基地| 免费看不卡的av| 亚洲美女黄色视频免费看| 日本av手机在线免费观看| 999精品在线视频| 成年人免费黄色播放视频| 久久精品久久久久久噜噜老黄| 美女福利国产在线| 男女午夜视频在线观看 | 日韩人妻精品一区2区三区| 男女国产视频网站| 建设人人有责人人尽责人人享有的| 亚洲精品一二三| 一区在线观看完整版| 美女福利国产在线| 男女午夜视频在线观看 | av国产精品久久久久影院| 国产一级毛片在线| 色94色欧美一区二区| 桃花免费在线播放| 中文乱码字字幕精品一区二区三区| 考比视频在线观看| 欧美日韩成人在线一区二区| 成人18禁高潮啪啪吃奶动态图| 亚洲成人一二三区av| 亚洲中文av在线| 一区二区av电影网| 午夜老司机福利剧场| 亚洲国产看品久久| 欧美精品人与动牲交sv欧美| 在线观看国产h片| 在线亚洲精品国产二区图片欧美| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频 | 观看美女的网站| 国产深夜福利视频在线观看| 狠狠婷婷综合久久久久久88av| 午夜激情久久久久久久| 少妇熟女欧美另类| 欧美激情极品国产一区二区三区 | 黄网站色视频无遮挡免费观看| 人妻一区二区av| 国产国语露脸激情在线看| 久久国内精品自在自线图片| 两个人免费观看高清视频| 两性夫妻黄色片 | 黄网站色视频无遮挡免费观看| 大香蕉久久网| 大码成人一级视频| 青春草国产在线视频| 久久久国产精品麻豆| xxx大片免费视频| 久久99热6这里只有精品| 精品酒店卫生间| 成人国语在线视频| 2021少妇久久久久久久久久久| 欧美丝袜亚洲另类| 国产成人精品福利久久| 一区二区三区四区激情视频| 日韩,欧美,国产一区二区三区| 午夜激情久久久久久久| 丝袜喷水一区| 国产精品偷伦视频观看了| 亚洲高清免费不卡视频| 精品视频人人做人人爽| 午夜av观看不卡| videossex国产| 午夜免费男女啪啪视频观看| 免费在线观看黄色视频的| 国产老妇伦熟女老妇高清| 国产精品不卡视频一区二区| 国产一区二区在线观看av| 九色亚洲精品在线播放| 99国产精品免费福利视频| 国产精品国产三级国产av玫瑰| 秋霞在线观看毛片| 人人妻人人澡人人爽人人夜夜| 有码 亚洲区| 男人舔女人的私密视频| 伊人亚洲综合成人网| 久久国产亚洲av麻豆专区| 国产色爽女视频免费观看| 91国产中文字幕| 中文字幕精品免费在线观看视频 | 精品久久蜜臀av无| 99热这里只有是精品在线观看| 精品国产一区二区三区久久久樱花| 欧美精品av麻豆av| 国产乱人偷精品视频| 18在线观看网站| 一区二区三区乱码不卡18| 国产伦理片在线播放av一区| 成人国产av品久久久| 亚洲欧美日韩卡通动漫| 国产精品久久久av美女十八| 亚洲色图 男人天堂 中文字幕 | 久久精品国产亚洲av天美| 91午夜精品亚洲一区二区三区| 视频中文字幕在线观看| 少妇精品久久久久久久| 久久99精品国语久久久| 色视频在线一区二区三区| 国产色婷婷99| 久久久久精品久久久久真实原创| 黄色视频在线播放观看不卡| 90打野战视频偷拍视频| 狂野欧美激情性xxxx在线观看| 国产日韩欧美亚洲二区| 黄片无遮挡物在线观看| 色婷婷久久久亚洲欧美| av有码第一页| 中文字幕最新亚洲高清| 老女人水多毛片| av免费观看日本| 一二三四在线观看免费中文在 | 久久精品国产自在天天线| 亚洲欧美一区二区三区国产| 最近中文字幕2019免费版| 综合色丁香网| 亚洲丝袜综合中文字幕| 国产成人91sexporn| 亚洲成人一二三区av| 久久久久久久久久久久大奶| 成人无遮挡网站| 亚洲性久久影院| 九色成人免费人妻av| av在线app专区| 欧美人与性动交α欧美精品济南到 | 90打野战视频偷拍视频| 9色porny在线观看| 免费黄频网站在线观看国产| 日本色播在线视频| 人妻少妇偷人精品九色| 人体艺术视频欧美日本| av国产精品久久久久影院| 97超碰精品成人国产| 欧美另类一区| 国产精品久久久久久精品电影小说| 伦精品一区二区三区| √禁漫天堂资源中文www| 久久精品熟女亚洲av麻豆精品| av免费在线看不卡| 丝瓜视频免费看黄片| 精品一区二区三区四区五区乱码 | 校园人妻丝袜中文字幕| 免费人成在线观看视频色| 亚洲国产av影院在线观看| 美女视频免费永久观看网站| 日韩av不卡免费在线播放| 在现免费观看毛片| 美女脱内裤让男人舔精品视频| 中文字幕av电影在线播放| 人妻一区二区av| 91久久精品国产一区二区三区| 国产xxxxx性猛交| 十分钟在线观看高清视频www| 内地一区二区视频在线| 国产精品偷伦视频观看了| 狂野欧美激情性bbbbbb| 国产福利在线免费观看视频| 亚洲综合精品二区| 免费av不卡在线播放| 高清不卡的av网站| 中国美白少妇内射xxxbb| 一级毛片电影观看| 男女无遮挡免费网站观看| 七月丁香在线播放| 伊人久久国产一区二区| 五月伊人婷婷丁香| 国产精品久久久久久精品古装| 在线天堂最新版资源| 日日摸夜夜添夜夜爱| 亚洲精品美女久久久久99蜜臀 | 国产精品一区二区在线不卡| 两个人免费观看高清视频| 九色成人免费人妻av| 国产免费又黄又爽又色| 亚洲精品色激情综合| 少妇精品久久久久久久| 99久久中文字幕三级久久日本| av在线老鸭窝| 国国产精品蜜臀av免费| 日韩,欧美,国产一区二区三区| 亚洲精品国产av成人精品| 天天操日日干夜夜撸| 欧美成人午夜精品| 国产一区二区在线观看av| 日韩电影二区| 国产精品免费大片| 亚洲国产精品999| 日日爽夜夜爽网站| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 欧美激情 高清一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美成人精品一区二区| 久久午夜综合久久蜜桃| 久久毛片免费看一区二区三区| 激情五月婷婷亚洲| 巨乳人妻的诱惑在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 青春草视频在线免费观看| 精品久久蜜臀av无| 日韩视频在线欧美| 免费久久久久久久精品成人欧美视频 | 国产精品蜜桃在线观看| 中文字幕人妻熟女乱码| 久久久久久久国产电影| 草草在线视频免费看| 欧美成人午夜免费资源| 精品久久久久久电影网| 91成人精品电影|