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

    沉淀-陳化-返溶作用和壓力對(duì)熱液中氟鈦絡(luò)合物高溫水解的影響及地質(zhì)意義*

    2015-07-21 08:50:36何俊杰丁興王玉榮孫衛(wèi)東
    巖石學(xué)報(bào) 2015年7期
    關(guān)鍵詞:沉淀物陳化淬火

    何俊杰 丁興 王玉榮 孫衛(wèi)東,4

    1.中國(guó)科學(xué)院廣州地球化學(xué)研究所礦物學(xué)與成礦學(xué)重點(diǎn)實(shí)驗(yàn)室,廣州 510640

    2.中國(guó)科學(xué)院大學(xué),北京 100049

    3.中國(guó)科學(xué)院廣州地球化學(xué)研究所同位素地球化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣州 510640

    4.中國(guó)科學(xué)院青藏高原地球科學(xué)卓越創(chuàng)新中心,北京 100101

    1 引言

    元素在地質(zhì)過(guò)程中以什么形式遷移以及如何遷移,一直以來(lái)就是地球化學(xué)界最基礎(chǔ)的科學(xué)問(wèn)題。在水溶液中,一些陽(yáng)離子、陰離子或電中性分子能絡(luò)合在一起形成一個(gè)穩(wěn)定的團(tuán)簇分子,稱之為絡(luò)合物(Ahrland et al.,1958;Baes and Mesmer,1976,1981);在地質(zhì)過(guò)程的高溫高壓條件下,由于組分復(fù)雜、揮發(fā)性組分濃度高、密度大等特點(diǎn),非常有利于絡(luò)合物的形成(王玉榮,1978)。因此,通常認(rèn)為,流體主導(dǎo)的地質(zhì)過(guò)程中,元素主要以配位型絡(luò)合物和多酸型絡(luò)合物存在并進(jìn)行遷移(王玉榮,1978;Taylor and Wall,1993;Bayot and Devillers,2006; Liu et al., 2007; Stefánsson, 2007;Antignano and Manning,2008;Williams-Jones et al.,2009;Yardley and Bodnar,2014)。自從20 世紀(jì)60 年代絡(luò)合物的概念被引入地球化學(xué)領(lǐng)域,元素絡(luò)合物在高溫高壓狀態(tài)下的賦存狀態(tài)以及穩(wěn)定性研究,對(duì)于理解地質(zhì)過(guò)程中元素的遷移以及豐富、完善成礦理論具有積極的意義。

    高場(chǎng)強(qiáng)元素(HFSE)通常被認(rèn)為是流體中最不活潑的金屬元素(Pearce and Peate,1995)。然而,在遭受后期流體作用的花崗巖中的黑云母會(huì)顯示不同程度的Ti 丟失(丁興等,2012);另一方面在自然界的高壓脈體中也經(jīng)??梢钥吹浇鸺t石(Gao et al.,2007;Guo et al.,2012;Huang et al.,2012),而金紅石是深部地質(zhì)過(guò)程中Ti、Nb、Ta 等HFSE 的主要載體(Xiao et al.,2006;Liang et al.,2009;Ding et al.,2009,2013;Meinhold,2010),這些觀測(cè)都說(shuō)明不管在低溫低壓環(huán)境還是高溫高壓環(huán)境,HFSE 在流體中都顯示一定的活動(dòng)性。而HFSE 在這些地質(zhì)過(guò)程中的活動(dòng)性可能事關(guān)大陸地殼和/或俯沖帶弧殼的形成和Nb、Ta、Ti 等稀有金屬礦床的成因(Ryerson and Watson,1987;Ayers and Watson,1991,1993;Brenan et al.,1994;Tropper and Manning,2005;Xiao et al.,2006;Gaetani et al.,2008;Liang et al.,2009;Ding et al.,2009,2013;Rapp et al.,2010;Linnen et al.,2013)。金紅石溶解度的高溫高壓實(shí)驗(yàn)發(fā)現(xiàn),金紅石在含F(xiàn) 流體中的溶解度要遠(yuǎn)遠(yuǎn)高于其他成分流體中的溶解度(Barsukova et al.,1979;Bright and Readey,1987;Rapp et al.,2010),這表明HFSE 易于在富氟流體中溶解并進(jìn)入流體,因而被認(rèn)為它們能在俯沖帶環(huán)境進(jìn)行大規(guī)模遷移的證據(jù)(Rapp et al.,2010)。然而,礦物溶解度的高低只能衡量源區(qū)中可以溶出元素最大的量(何俊杰等,2015),流體中元素的遷移量及遷移尺度需由穩(wěn)定的遷移形式——絡(luò)合物的穩(wěn)定性來(lái)決定(Baes and Mesmer,1976;Van Baalen,1993),而對(duì)于HFSE 絡(luò)合物在高溫高溫高壓流體中的穩(wěn)定性仍缺乏相關(guān)的基礎(chǔ)實(shí)驗(yàn)數(shù)據(jù)(何俊杰等,2015)。Ahrland et al.(1958)根據(jù)常溫常壓下的化學(xué)實(shí)驗(yàn)指出,HFSE 通常與C、O、F 形成穩(wěn)定的絡(luò)合物。Wang et al.(1993)使用快速淬火高壓釜通過(guò)高溫水解法研究了K2TaF7和K2NbOF5·5H2O 在不同溫度壓力下的水解行為,結(jié)果發(fā)現(xiàn)它們隨著溫度的升高發(fā)生強(qiáng)烈的水解,而壓力的升高也稍微增強(qiáng)水解。這暗示在自然地質(zhì)過(guò)程中,大多數(shù)HFSE 含氟絡(luò)合物可能是不穩(wěn)定的,高溫高壓的條件并不利于HFSE 的大規(guī)模遷移,這與目前自然樣品的研究結(jié)論相悖(Xiao et al.,2006;Gao et al.,2007;Liang et al.,2009;Ding et al.,2009,2013;Huang et al.,2012)。何俊杰等(2015)運(yùn)用類似方法研究了不同濃度的K2TiF6在不同溫度下的穩(wěn)定性,結(jié)果顯示溫度越高、初始濃度越低,氟鈦絡(luò)合物的水解程度越劇烈,該規(guī)律與Nb、Ta類似。何俊杰等(2015)認(rèn)為高溫高壓條件下HFSE 的遷移通常是受限制的,HFSE 的大規(guī)模遷移需要出現(xiàn)富氟流體才能實(shí)現(xiàn)。然而,何俊杰等(2015)并沒(méi)有考慮壓力的改變對(duì)K2TiF6穩(wěn)定性的影響,且絡(luò)合物水解的細(xì)節(jié)過(guò)程也并沒(méi)有關(guān)注及討論。從以上可以看出,由于自然條件和自然流體組分的復(fù)雜性,HFSE 絡(luò)合物在高溫高壓流體中的穩(wěn)定性仍需要進(jìn)一步深入研究。

    本文在以上研究的基礎(chǔ)上,進(jìn)一步開展對(duì)氟鈦絡(luò)合物的穩(wěn)定性及水解行為研究,著重探討沉淀、陳化和返溶作用對(duì)氟鈦絡(luò)合物水解行為的影響,以及壓力對(duì)氟鈦絡(luò)合物穩(wěn)定性的影響,試圖刻畫氟鈦絡(luò)合物更清晰的水解行為和流體中礦物生長(zhǎng)的機(jī)制,從而為元素絡(luò)合物的高溫水解以及HFSE 的活動(dòng)性研究提供新的信息和啟示。

    2 實(shí)驗(yàn)方法和分析方法

    2.1 實(shí)驗(yàn)方法

    為了研究自然地質(zhì)過(guò)程中Ti 的最大遷移行為,本研究采用K2TiF6作為實(shí)驗(yàn)的初始反應(yīng)物。K2TiF6具有六個(gè)氟配位結(jié)構(gòu),是目前已知最穩(wěn)定且具有最高氟配位的氟鈦無(wú)機(jī)絡(luò)合物(Lee et al.,2009;Banerjee,2011)。其中的氟酸基絡(luò)合陰離子[TiF6]2-在水中會(huì)發(fā)生逐級(jí)電離和水解,過(guò)程如下:

    以上過(guò)程會(huì)造成多種氟酸基和/或羥基氟鈦絡(luò)合物等種型在溶液中共存,同時(shí)水解的最終理想產(chǎn)物為TiO2沉淀。由于溫度、壓力、pH 等條件的不同,溶液中共存的絡(luò)合物種型以及TiO2的產(chǎn)出量也相應(yīng)不同,這就為利用高溫水解法研究絡(luò)合物的穩(wěn)定性提供了理論依據(jù)。

    本文所有高溫高壓實(shí)驗(yàn)和分析測(cè)試均在中國(guó)科學(xué)院廣州地球化學(xué)研究所同位素地球化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室水熱實(shí)驗(yàn)平臺(tái)完成。

    2.1.1 絡(luò)合物的高溫水解實(shí)驗(yàn)

    為了研究氟鈦絡(luò)合物在不同壓力條件下的穩(wěn)定性,我們使用了相同濃度的K2TiF6溶液在不同壓力下分別進(jìn)行高溫高壓實(shí)驗(yàn),根據(jù)以上反應(yīng)(1)~(6),可獲得綜合的理想反應(yīng)方程式:

    在溫度升高的情況下,K2TiF6會(huì)發(fā)生顯著水解,產(chǎn)生大量TiO2的沉淀物,顯示氟鈦絡(luò)合物在高溫高壓下穩(wěn)定性較常溫常壓下要降低(何俊杰等,2015)。我們通過(guò)測(cè)量特定溫壓條件下反應(yīng)(7)的反應(yīng)程度來(lái)了解氟鈦絡(luò)合物在目標(biāo)條件下的穩(wěn)定性,其中表征反應(yīng)程度的參數(shù)有水解率和累積水解常數(shù),其中:

    累積水解常數(shù)即反應(yīng)(7)的反應(yīng)平衡常數(shù),如方程(8)所示。由于水解反應(yīng)中各物質(zhì)濃度都比較低且缺乏高溫下的活度系數(shù)數(shù)據(jù),而在低濃度下溶液中離子的活度系數(shù)近似為1,因此方程(8)中使用了物質(zhì)的量濃度來(lái)代替活度(何俊杰等,2015)。

    本研究所用的K2TiF6粉晶均為即時(shí)合成:將分析純度的TiO2粉末和KHF2按1∶4.5 的比例在鉑金坩堝中充分混合并在350℃下保溫熔融15min,冷卻后使用5%氫氟酸將混合物溶解,過(guò)濾重結(jié)晶后制得K2TiF6粉晶,反應(yīng)方程式如方程(9)所示。

    合成產(chǎn)物經(jīng)粉晶X 射線衍射分析(XRD)后獲得的譜圖與參考譜圖一致,且成分較為純凈(何俊杰等,2015)。將48mg K2TiF6粉晶溶解于去100mL 離子水中可獲得鈦質(zhì)量濃度為960 ×10-6的K2TiF6溶液,經(jīng)分光光度計(jì)分析后顯示其實(shí)際濃度為(947 ±19)×10-6。

    反應(yīng)初始溶液在裝進(jìn)Ф4mm 金管并焊封后放進(jìn)水熱高壓釜中進(jìn)行高溫高壓實(shí)驗(yàn),在實(shí)驗(yàn)過(guò)程中的每一步操作后都對(duì)金管進(jìn)行稱重,以檢查金管有無(wú)泄露。實(shí)驗(yàn)中我們把反應(yīng)溫度統(tǒng)一控制為300℃,反應(yīng)壓力則依次設(shè)定為50MPa、100MPa、150MPa、200MPa、300MPa、400MPa。其 中50~200MPa 實(shí)驗(yàn)在快速淬火高壓釜反應(yīng)系統(tǒng)中完成,300~400MPa 實(shí)驗(yàn)在立式的冷封高壓釜反應(yīng)系統(tǒng)中完成。

    快速淬火高壓釜反應(yīng)系統(tǒng)采用氬氣作為壓力源,其反應(yīng)壓力極限為200MPa。反應(yīng)結(jié)束后樣品采用快速淬火處理:冷封端冰敷后倒置使得樣品迅速?gòu)母邷貐^(qū)掉落到低溫區(qū),樣品可在3s 內(nèi)迅速降溫至室溫以下。立式的Tuttle 型冷封高壓釜系統(tǒng)采用純水作為壓力介質(zhì),反應(yīng)壓力極限可達(dá)500MPa。反應(yīng)結(jié)束后直接移開加溫爐,使用冰桶對(duì)高壓釜熱端進(jìn)行快速淬火,可使樣品在5s 內(nèi)降溫至100℃以下,在半分鐘內(nèi)降溫至室溫以下。所有樣品的反應(yīng)時(shí)間均為12h,反應(yīng)結(jié)束后迅速打開金管,使用一次性針管提取上層清液,以4000 轉(zhuǎn)/分鐘轉(zhuǎn)速離心15min 后把溶液稀釋待測(cè)。

    為了了解并分析實(shí)驗(yàn)過(guò)程中產(chǎn)生的實(shí)驗(yàn)誤差,我們對(duì)每一個(gè)壓力下的實(shí)驗(yàn)設(shè)置了重復(fù)實(shí)驗(yàn):50~200MPa 的實(shí)驗(yàn)分別進(jìn)行四次重復(fù)實(shí)驗(yàn),300~400MPa 的實(shí)驗(yàn)分別進(jìn)行三次重復(fù)實(shí)驗(yàn),所有重復(fù)實(shí)驗(yàn)在實(shí)驗(yàn)條件設(shè)置及實(shí)驗(yàn)設(shè)備上均為一致。

    2.1.2 絡(luò)合物的返溶實(shí)驗(yàn)

    從反應(yīng)(7)中可以看到,實(shí)驗(yàn)產(chǎn)物中同時(shí)存在TiO2沉淀物和HF,在反應(yīng)條件下兩者可以平衡共存,而在反應(yīng)結(jié)束降溫的過(guò)程中由于發(fā)生平衡的移動(dòng)使得TiO2沉淀可能會(huì)重新溶解于HF 中而使Ti 重新進(jìn)入流體,從而使得實(shí)驗(yàn)測(cè)定結(jié)果偏高。為了了解TiO2在反應(yīng)產(chǎn)物流體中的返溶情況,我們?cè)O(shè)計(jì)了兩組實(shí)驗(yàn)分別測(cè)定了TiO2沉淀在淬火過(guò)程中高壓釜內(nèi)的返溶情況和淬火后沉淀物在常溫常壓下的返溶情況。

    所有的返溶實(shí)驗(yàn)均在100MPa 壓力條件下進(jìn)行,其初始溶液組成與前文反應(yīng)初始物一致,鈦質(zhì)量濃度均為(947 ±19)×10-6。高溫高壓下的返溶實(shí)驗(yàn)在立式的冷封高壓釜系統(tǒng)中進(jìn)行。我們先把樣品放進(jìn)高壓釜中,在300℃溫度下反應(yīng)12h,然后把體系溫度緩慢降低至200℃并保溫2h,隨后將樣品淬火并立即提取溶液待測(cè)。將測(cè)試結(jié)果與K2TiF6在200℃和300℃下的結(jié)果進(jìn)行對(duì)比即可了解緩慢降溫過(guò)程中TiO2的返溶情況。常溫常壓下的返溶實(shí)驗(yàn)使用的是快速淬火高壓釜,我們先把初始溶液裝進(jìn)Ф6mm 金管中在300℃溫度條件下進(jìn)行反應(yīng)12h,待淬火后依次在不同的時(shí)間分批提取溶液,溶液提取時(shí)間分別為0、2、5、9、12、23h。在取樣過(guò)程中,我們使用鎢針把金管戳出一個(gè)小孔,然后使用針管提取金管內(nèi)的溶液,每次提取溶液后均使用透明膠把小孔密封以排除金管內(nèi)溶液揮發(fā)造成的誤差。最后將不同時(shí)間提取溶液的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比即可了解在常溫常壓下TiO2的返溶情況。

    2.2 分析方法

    反應(yīng)后提取的溶液中的鈦濃度分析均在上海菁華7600-1CRT 型雙光束紫外可見光分光光度計(jì)(UV-Vis Spectrophotometer)上進(jìn)行。該設(shè)備波長(zhǎng)范圍為190~1100nm,準(zhǔn)確度為0.3nm,重復(fù)度為0.1nm,光譜帶寬1nm。分析使用二安替比林甲烷(DAPM)顯色法對(duì)溶液中的鈦濃度進(jìn)行測(cè)定,測(cè)試原理為DAPM 會(huì)取代TiF62-中的氟與鈦形成穩(wěn)定的黃色絡(luò)合物,在波長(zhǎng)為385 nm 處出現(xiàn)最大吸收,通過(guò)標(biāo)準(zhǔn)曲線法可準(zhǔn)確獲得溶液中的鈦濃度。分析過(guò)程中所用到的鈦標(biāo)準(zhǔn)溶液為國(guó)家標(biāo)準(zhǔn)樣品GSB 04-1757-2004,每一次標(biāo)準(zhǔn)曲線的建立其相關(guān)系數(shù)R2均大于0.9995。對(duì)每一個(gè)待測(cè)樣品我們分別進(jìn)行6 次平行測(cè)定以校正在濃度測(cè)定過(guò)程中出現(xiàn)的分析誤差。

    為了驗(yàn)證測(cè)試方法的分析準(zhǔn)確性,我們使用分光光度法對(duì)反應(yīng)初始溶液進(jìn)行了平衡測(cè)定。平衡測(cè)定次數(shù)為24 次,分三批進(jìn)行以避免不同測(cè)量條件下的分析誤差,測(cè)量結(jié)果中鈦的質(zhì)量濃度為(947 ±19)×10-6,相對(duì)標(biāo)準(zhǔn)誤差為1.3%;而且,Ti 濃度的分光光度計(jì)法與使用Vista Pro 型全譜直讀電感耦合等離子體原子發(fā)射光譜儀(ICP-AES)分析獲得的結(jié)果相比較,兩者在誤差范圍內(nèi)基本一致(何俊杰等,2015)。這些均說(shuō)明使用分光光度法測(cè)量反應(yīng)后提取的溶液中的鈦濃度是可靠的。

    3 實(shí)驗(yàn)結(jié)果與討論

    3.1 氟鈦絡(luò)合物水解過(guò)程中的沉淀、陳化和返溶作用

    隨著時(shí)間的增長(zhǎng),流體中的無(wú)定型沉淀會(huì)自發(fā)地進(jìn)行溶解-重結(jié)晶,使沉淀物轉(zhuǎn)化為當(dāng)前物理化學(xué)條件下最穩(wěn)定的晶型,即體系的自由能降至最低,該過(guò)程稱為沉淀的陳化(aging)(Jalava et al.,1998;Namin et al.,2008;Schott et al.,2009;Sileri et al.,2011;Fatmi et al.,2013;Palacios et al.,2013)。正如上述反應(yīng)式(6)所示,氟鈦絡(luò)合物水解反應(yīng)的結(jié)果首先是生產(chǎn)了Ti(OH)4含水膠狀沉淀(Knauss et al.,2001;Ryzhenko et al.,2006;Hayden and Manning,2011;Mysen,2012),隨后發(fā)生陳化脫水,繼而生成TiO2納米顆粒(Jalava et al.,1998;Pottier et al.,2003;Namin et al.,2008;Lee et al.,2009;Banerjee,2011),繼續(xù)的陳化作用會(huì)促進(jìn)晶核的形成(Schott et al.,2009;Palacios et al.,2013),直至在不同溫壓條件下轉(zhuǎn)化為銳鈦礦、板鈦礦或金紅石(Jalava et al.,1998;Pottier et al.,2003;Namin et al.,2008)。雖然在這過(guò)程中一直伴隨著Ti(OH)4、TiO2納米顆粒等的返溶,但隨著固體沉淀粒徑增大,沉淀的反應(yīng)活性也降低,因此,通常認(rèn)為,對(duì)于高價(jià)元素而言,其沉淀引起的陳化作用將更明顯;正因?yàn)槌恋砦锏年惢^(guò)程,沉淀物的返溶往往不易于進(jìn)行(王玉榮等,1995;Knauss et al.,2001)。

    表1 為480 ×10-6K2TiF6溶液在300℃和100MPa 條件下進(jìn)行12h 水解反應(yīng)后的產(chǎn)物在常溫常壓下開展的返溶實(shí)驗(yàn)結(jié)果。從圖1 中可以看到隨著返溶時(shí)間的增長(zhǎng),溶液中的鈦含量并沒(méi)有隨之升高,水解率變化范圍從90.0% 到90.4%,基本上可以認(rèn)為在誤差范圍內(nèi)一致。這一方面證實(shí)了淬火后TiO2沉淀物經(jīng)歷了顯著的陳化作用;另一方面也表明采用本文高溫水解法淬火后的溶液在待測(cè)、提取及離心等過(guò)程中并不會(huì)因?yàn)門iO2沉淀物的明顯返溶而造成溶液成分的變化。

    圖1 常溫常壓下TiO2 沉淀返溶實(shí)驗(yàn)結(jié)果Fig.1 Results of TiO2 precipitate re-dissolving experiment at ambient temperature and pressure condition

    表1 常溫常壓下TiO2 沉淀返溶實(shí)驗(yàn)結(jié)果Table 1 Results of TiO2 precipitate re-dissolving experiment at ambient temperature and pressure condition

    根據(jù)之前的研究發(fā)現(xiàn),Nb、Ta 和Ti 等HFSE 的富氟絡(luò)合物的水解程度與溫度成正相關(guān)性(Wang et al.,1993;何俊杰等,2015)。這表明,如果樣品在淬火過(guò)程中淬火速率過(guò)慢或者淬火操作失誤導(dǎo)致淬火延遲,則樣品溫度會(huì)在某一時(shí)間段內(nèi)處在反應(yīng)溫度和室溫之間,從而造成TiO2沉淀會(huì)發(fā)生返溶并促使反應(yīng)平衡向低溫方向移動(dòng)。圖2 為TiO2沉淀在高溫高壓條件下的返溶實(shí)驗(yàn)結(jié)果,樣品在300℃下反應(yīng)12h 后緩慢降溫至200℃,保溫2h 后淬火。實(shí)驗(yàn)產(chǎn)物溶液中Ti 的濃度為228.7 ×10-6,K2TiF6水解率為75.9%(圖2 中實(shí)線所示),而之前的研究結(jié)果顯示,在相同初始濃度和壓力條件下,K2TiF6在300℃和200℃的水解率分別是83.1%和61.5%,溶液中的Ti 濃度分別為162 ×10-6和370 ×10-6(何俊杰等,2015),由此可知本返溶實(shí)驗(yàn)中K2TiF6水解率介于兩者之間。這一實(shí)驗(yàn)現(xiàn)象一方面證明在高溫高壓條件下緩慢的降溫過(guò)程會(huì)使得TiO2沉淀發(fā)生返溶作用,即部分TiO2沉淀重新溶解于富含HF 和F-的高溫流體中;另一方面暗示由于沉淀物的陳化作用,返溶的速率并不快,歷經(jīng)2 個(gè)小時(shí),返溶還未達(dá)到平衡。按照Ti 的返溶量來(lái)估算,2 個(gè)小時(shí)的返溶速率為0.56 ×10-6min-1;如果按照這種返溶速率來(lái)計(jì)算,則從300℃降溫至200℃時(shí),TiO2沉淀返溶達(dá)到平衡則至少需要6 個(gè)小時(shí)。因此,盡管在高溫高壓條件下的降溫過(guò)程會(huì)造成返溶作用的發(fā)生,但是,由于沉淀物的陳化作用,其返溶速率并不快。因而采用快速淬火設(shè)備和/或快速淬火技術(shù)可以有效降低返溶作用的發(fā)生,使得淬火后流體中的成分能更接近代表高溫高壓下流體的情況。

    圖2 高溫高壓條件下TiO2 返溶實(shí)驗(yàn)結(jié)果Fig.2 Results of TiO2 re-dissolving experiment at HTHP conditions

    對(duì)于大多數(shù)高溫高壓實(shí)驗(yàn)而言,實(shí)驗(yàn)中往往容易出現(xiàn)較為顯著的誤差,其實(shí)驗(yàn)數(shù)據(jù)的誤差往往來(lái)源于:1)高溫高壓實(shí)驗(yàn)誤差;2)分析誤差。前者往往由于高溫高壓實(shí)驗(yàn)過(guò)程中的溫度及壓力控制存在偏差(Rudert et al.,1976;Charles and Vidale,1982;Dalton and Gibb,1996),或者是在樣品淬火過(guò)程中樣品成分發(fā)生改變導(dǎo)致實(shí)驗(yàn)結(jié)果偏離高溫高壓的真實(shí)情況(Jacobs and Kerrick,1980;Holloway et al.,1992)。而后者則主要取決于分析方法的誤差,在大多數(shù)情況下可以通過(guò)對(duì)測(cè)試設(shè)備的調(diào)試校正或測(cè)試方法的改進(jìn)進(jìn)行控制。盡管在采用高溫水解法研究絡(luò)合物的穩(wěn)定性過(guò)程中,使用快速淬火設(shè)備和/或快速淬火技術(shù)可以理論上有效規(guī)避沉淀物返溶作用造成的高溫高壓實(shí)驗(yàn)誤差,但是在本次研究中我們發(fā)現(xiàn)使用傳統(tǒng)的快速淬火設(shè)備也會(huì)造成較大的實(shí)驗(yàn)誤差。從表2 中我們可以看到,50~200MPa 的每一組平行實(shí)驗(yàn)測(cè)定值中都出現(xiàn)了明顯的數(shù)據(jù)偏差現(xiàn)象,如在50MPa 的實(shí)驗(yàn)結(jié)果中,結(jié)果最大值高達(dá)(130.6 ±1.0)×10-6,而結(jié)果最小值只有(36.1 ±0.9)×10-6,這明顯已經(jīng)超出了實(shí)驗(yàn)所允許的誤差范圍。這種現(xiàn)象在水熱實(shí)驗(yàn)中較為普遍,通常被認(rèn)為是快速淬火效率不同所致(Keppler and Audétat,2005)。因?yàn)榇慊鹚俾实牟煌?,?huì)造成流體中的溶解物質(zhì)發(fā)生不同程度的沉淀,或者沉淀物質(zhì)發(fā)生不同程度的返溶,從而出現(xiàn)實(shí)驗(yàn)結(jié)果的顯著偏離。在我們的實(shí)驗(yàn)中,50~200MPa 的實(shí)驗(yàn)是在快速淬火高壓釜上完成的,采用包裹著碎冰的毛巾冰敷在高壓釜的冷封端使冷封端溫度降低至室溫以下,然后將高壓釜倒置使金管自然掉落至低溫區(qū)從而實(shí)現(xiàn)淬火。這種淬火速率是目前已知最快的方式;而且,如上述,由于TiO2沉淀物的陳化作用,在降溫過(guò)程中的返溶速率并不快,因此不可能造成如此大的實(shí)驗(yàn)誤差。有意思的是,300~400MPa 的平行實(shí)驗(yàn)結(jié)果卻在很小范圍內(nèi)波動(dòng)(表2、圖3)。該實(shí)驗(yàn)是在立式的冷封高壓釜系統(tǒng)上完成的,采用冰桶直接冷卻高壓釜熱端的快速淬火方式。與上述實(shí)驗(yàn)方法相比,雖然該方法有稍微慢的快速淬火速率(王玉榮和周義明,1986),但是在淬火過(guò)程中,高壓釜體及樣品管一直處于靜止?fàn)顟B(tài)??紤]到礦物的溶解或水-巖反應(yīng)過(guò)程中,體系的振蕩可以加速反應(yīng)速率(Keppler and Audétat,2005;Schott et al.,2009),因此,我們認(rèn)為,50~200MPa 平行實(shí)驗(yàn)測(cè)量值中出現(xiàn)的明顯的變化是由于淬火過(guò)程中樣品管從高溫區(qū)掉落至低溫區(qū)引起的振蕩強(qiáng)化了返溶作用造成的。因此我們建議,為了降低這種機(jī)械性返溶作用造成的影響,采用快速淬火高壓釜系統(tǒng)進(jìn)行絡(luò)合物的高溫水解實(shí)驗(yàn),需要進(jìn)行多次平行實(shí)驗(yàn),并以最低值為基準(zhǔn)選取與最低值接近的平行測(cè)量值為有效數(shù)據(jù)點(diǎn),繼而可以有效降低高溫高壓實(shí)驗(yàn)誤差。

    表2 壓力對(duì)K2TiF6 高溫高壓下水解影響平行實(shí)驗(yàn)結(jié)果Table 2 Parallel experimental result of the effects of pressure on K2TiF6 hydrolysis at HTHP

    圖3 在300℃溫度條件下壓力對(duì)K2TiF6 水解影響的平行實(shí)驗(yàn)結(jié)果Fig.3 Parallel experimental result of the effects of pressure on K2TiF6 hydrolysis at HTHP

    綜上所述,我們的實(shí)驗(yàn)結(jié)果證實(shí),在氟鈦絡(luò)合物水解的過(guò)程中,由于TiO2沉淀物的陳化作用,在穩(wěn)定的溫度和壓力條件下,沉淀物的返溶很難發(fā)生;但在緩慢降溫和機(jī)械性振蕩過(guò)程中,沉淀物的返溶會(huì)明顯發(fā)生。為了降低返溶作用對(duì)實(shí)驗(yàn)結(jié)果的影響,可以采取快速淬火技術(shù)和多次平行實(shí)驗(yàn)取低值的方法。

    表3 不同壓力下K2TiF6 溶液在高溫高壓下水解情況Table 3 Hydrolysis results of K2TiF6 solution in HTHP under the varied pressures

    3.2 壓力對(duì)氟鈦絡(luò)合物水解的影響

    正如上述,為了降低在實(shí)驗(yàn)淬火過(guò)程中TiO2沉淀物的返溶作用對(duì)實(shí)驗(yàn)結(jié)果的影響,本文對(duì)50~200MPa 實(shí)驗(yàn)的結(jié)果采用了低值原則進(jìn)行處理,即排除明顯偏高的數(shù)據(jù),將相對(duì)較低的數(shù)據(jù)進(jìn)行平均計(jì)算。表3 為K2TiF6在不同壓力下的水解率和累積水解常數(shù)的結(jié)果。

    從表3 和圖4 中可以看出,K2TiF6的水解率并沒(méi)有跟壓力成明顯相關(guān)關(guān)系,其變化范圍為94.6%~96.4%,壓力對(duì)其的影響明顯小于溫度和濃度(何俊杰等,2015)。這說(shuō)明壓力的改變對(duì)K2TiF6在高溫高壓下的水解并沒(méi)有顯著的影響。Wang et al.(1993)曾對(duì)鈮和鉭的富氟絡(luò)合物水解情況受壓力的影響進(jìn)行研究,其高溫高壓實(shí)驗(yàn)結(jié)果顯示隨著壓力增大,絡(luò)合物的水解程度稍微增強(qiáng)。但值得注意的是,該實(shí)驗(yàn)結(jié)果顯示水解增強(qiáng)的幅度并不大,如果考慮實(shí)驗(yàn)誤差的話,基本可認(rèn)為在一個(gè)范圍內(nèi)(圖4)。因此,不管是Ti 還是Nb、Ta,其富氟絡(luò)合物在高溫高壓下的穩(wěn)定性都主要受溫度控制,而壓力對(duì)其影響并不顯著。這從另一側(cè)面可能暗示,HFSE 絡(luò)合物發(fā)生水解產(chǎn)生沉淀后,壓力的增高會(huì)促進(jìn)沉淀物的陳化,進(jìn)而抑制沉淀物的返溶。也就是說(shuō),壓力施加的效果更多體現(xiàn)在維持水解反應(yīng)平衡,而很可能并不會(huì)顯著改變水解反應(yīng)的進(jìn)行程度。這一特點(diǎn)支持了何俊杰等(2015)提出的K2TiF6高溫高壓下的累積水解常數(shù)與溫度之間的擬合關(guān)系。而這也暗示了熱液流體中的氟鈦絡(luò)合物可能不管是在俯沖帶這種高壓環(huán)境中還是巖漿熱液體系這種低壓環(huán)境中均遵循著近似的遷移規(guī)律。

    圖4 壓力改變對(duì)K2TiF6 水解率及累積水解常數(shù)的影響Fig.4 Effect of varied pressures on K2TiF6 hydrolysis ratios and apparent hydrolysis constants

    3.3 地質(zhì)意義

    金屬絡(luò)合物的形成-水解-沉淀/返溶-陳化是一個(gè)相互影響、相互制約的復(fù)雜過(guò)程,可能在水-巖作用或熱液成礦過(guò)程中扮演了關(guān)鍵的作用,是高場(chǎng)強(qiáng)元素活化、遷移和成礦的主要機(jī)制。對(duì)于HFSE 來(lái)說(shuō),其絡(luò)合物的形成,首先取決于金紅石、鈦鐵礦等賦存礦物的溶解。而之前的高溫高壓實(shí)驗(yàn)證實(shí),HFSE 在富氟流體中具有最大的溶解度(Barsukova et al.,1979;Bright and Readey,1987;Rapp et al.,2010)。因此富氟流體可以有效地活化全巖中的HFSE,從而形成含氟絡(luò)合物。期間,表現(xiàn)為HFSE 賦存礦物在流體里表面發(fā)生羥基化作用,HFSE 離子與水分子形成內(nèi)層絡(luò)合物和/或外層絡(luò)合物(Schott et al.,2009),隨后流體中的F-逐級(jí)置換羥基絡(luò)合物中的OH-,從而形成不同含氟絡(luò)合物(Ryzhenko et al.,2006;Schott et al.,2009;何俊杰等,2015)。在羥基絡(luò)合物與含氟絡(luò)合物形成的過(guò)程中,伴隨著水解的發(fā)生,部分HFSE 原位或半原位的沉淀下來(lái),并進(jìn)入次生的礦物晶格或生成獨(dú)立的金紅石(Keppler and Audétat,2005)。這一過(guò)程實(shí)質(zhì)上就是礦物的溶解/沉淀過(guò)程,代表了水-巖作用或變質(zhì)過(guò)程中流體對(duì)礦物的交代和次生礦物的原位或半原位的生長(zhǎng)。水解后殘留的HFSE 羥基絡(luò)合物全部轉(zhuǎn)化為含氟絡(luò)合物,并穩(wěn)定于流體中,可以進(jìn)行有效地遷移。何俊杰等(2015)通過(guò)計(jì)算認(rèn)為,在俯沖帶條件下,當(dāng)金紅石的溶解度大于1000 ×10-6時(shí),在流體中穩(wěn)定并可發(fā)生遷移的Ti 才多于1 ×10-6,而且隨著金紅石溶解度增大,遷移的量也會(huì)增大。這暗示從金紅石溶解出來(lái)的大部分Ti 在原位或半原位重新進(jìn)入了次生礦物,而高壓-超高壓巖石中靠近脈體且明顯受流體影響的榴輝巖往往也富集金紅石的現(xiàn)象也證實(shí)了這一點(diǎn)(Gao et al.,2007;Guo et al.,2012)。這些包含穩(wěn)定的HFSE 含氟絡(luò)合物的流體滲濾出巖石,并匯聚成脈,從而為HFSE 的遷移創(chuàng)造條件。根據(jù)本文研究的結(jié)果可知,壓力的變化并不會(huì)明顯影響HFSE 含氟絡(luò)合物的水解,因此,當(dāng)流體脈的溫度或pH 保持穩(wěn)定時(shí),HFSE含氟絡(luò)合物不會(huì)進(jìn)一步水解,HFSE 可以在流體中進(jìn)行遷移;當(dāng)流體脈的溫度降低時(shí),比如流體進(jìn)入淺部或者俯沖板塊折返,HFSE 含氟絡(luò)合物也不會(huì)進(jìn)一步水解,反而會(huì)發(fā)生返溶(何俊杰等,2015),即從原沉淀或圍巖中繼續(xù)溶解HFSE,從而增加流體中HFSE 的含量。然而,當(dāng)溫度升高或者pH 變化時(shí),導(dǎo)致反應(yīng)式(7)向右進(jìn)行,HFSE 含氟絡(luò)合物就會(huì)發(fā)生進(jìn)一步水解,形成沉淀。之后沉淀發(fā)生陳化作用形成晶核(Schott et al.,2009;Palacios et al.,2013),直至在不同溫壓條件下轉(zhuǎn)化為銳鈦礦、板鈦礦或金紅石(Jalava et al.,1998;Pottier et al.,2003;Namin et al.,2008;Meinhold,2010)。因此,在俯沖帶超高壓環(huán)境,高壓脈中的金紅石易于在增溫的進(jìn)變質(zhì)過(guò)程中形成;而在退變質(zhì)過(guò)程中,流體具有更高的HFSE 溶解和遷移能力,環(huán)境溫度的降低并不利于脈體礦物的結(jié)晶,因而部分退變質(zhì)脈的形成只可能歸因于流體成分的改變,比如與pH 高的流體混合,或者外來(lái)的Ca 吸收了流體中的F-,破壞了水解平衡,使得沉淀產(chǎn)生,繼而陳化成核,自然界高壓脈體中常見金紅石顆粒與富氟磷灰石共生(Gao et al.,2007)的現(xiàn)象很可能指示了這種機(jī)制。對(duì)于淺部的熱液成礦過(guò)程而言,流體從深部高溫環(huán)境進(jìn)入淺部低溫環(huán)境,溫度和壓力的降低,不會(huì)導(dǎo)致HFSE 的沉淀乃至成礦,反而會(huì)造成流體溶解和遷移HFSE 的能力增強(qiáng)。因此,要想使得HFSE 沉淀下來(lái)成礦,同樣需要外來(lái)組分的加入,比如與地表水等高pH 值的流體混合,或者外來(lái)的Ca 吸收了流體中的F。一個(gè)可能的典型例子是白云鄂博稀有金屬礦床的形成:含輕稀土、Nb 等組分的富氟流體進(jìn)入沉積碳酸鹽巖(Ling et al.,2013),造成酸性的富氟流體與堿性的碳酸巖發(fā)生中和反應(yīng),使得流體中的稀有金屬沉淀下來(lái)成礦,并生成大量伴生的螢石。

    4 結(jié)論

    本文對(duì)濃度為960 ×10-6的氟鈦絡(luò)合物(K2TiF6)溶液在300℃溫度和50~400MPa 壓力下的穩(wěn)定性進(jìn)行了研究,探討了絡(luò)合物水解過(guò)程中沉淀-陳化-返溶作用和壓力對(duì)氟鈦絡(luò)合物水解行為的影響,得到以下主要結(jié)論:

    (1)氟鈦絡(luò)合物水解的過(guò)程中,由于TiO2沉淀物的陳化作用,在穩(wěn)定的溫度和壓力條件下,沉淀物的返溶很難發(fā)生;但在緩慢降溫和機(jī)械性振蕩過(guò)程中,沉淀物的返溶會(huì)明顯發(fā)生。為了規(guī)避返溶作用對(duì)實(shí)驗(yàn)結(jié)果的影響,建議采取快速淬火技術(shù)和多次平行實(shí)驗(yàn)取低值的方法。

    (2)壓力的改變對(duì)K2TiF6在高溫高壓下的水解并沒(méi)有顯著的影響。結(jié)合前人的實(shí)驗(yàn)結(jié)果,認(rèn)為不管是Ti 還是Nb、Ta,其富氟絡(luò)合物在高溫高壓下的穩(wěn)定性都主要受溫度控制,而壓力對(duì)其影響并不顯著。而這也暗示了熱液中的氟鈦絡(luò)合物可能不管是在俯沖帶這種高壓環(huán)境中還是巖漿熱液體系這種低壓環(huán)境中均遵循著近似的遷移規(guī)律。

    (3)本文提出了金屬絡(luò)合物的形成-水解-沉淀/返溶-陳化是一個(gè)相互影響、相互制約的復(fù)雜過(guò)程,可能在水-巖作用或熱液成礦過(guò)程中扮演了關(guān)鍵的作用,也是高場(chǎng)強(qiáng)元素活化、遷移和成礦的主要機(jī)制。

    Ahrland S,Chatt J and Davies NR.1958.The relative affinities of ligand atoms for acceptor molecules and ions.Quarterly Review of the Chemical Society,12(3):265-276

    Antignano A and Manning CE.2008.Rutile solubility in H2O,H2OSiO2,and H2O-NaAlSi3O8fluids at 0.7~2.0GPa and 700~1000 degrees C:Implications for mobility of nominally insoluble elements.Chemical Geology,255(1-2):283-293

    Ayers JC and Watson EB.1991.Solubility of apatite,monazite,zircon,and rutile in supercritical aqueous fluids with implications for subduction zone geochemistry.Philosophical Transactions of the Royal Society of London Series a-Mathematical Physical and Engineering Sciences,335(1638):365-375

    Ayers JC and Watson EB.1993.Rutile solubility and mobility in supercritical aqueous fluids.Contributions to Mineralogy and Petrology,114(3):321-330

    Baes CF and Mesmer RE.1976.The Hydrolysis of Cations.New York:Wiley

    Baes CF and Mesmer RE.1981.The thermodynamics of cation hydrolysis.American Journal of Science,281(7):935-962

    Banerjee AN.2011.The design,fabrication,and photocatalytic utility of nanostructured semiconductors: Focus on TiO2-based nanostructures.Nanotechnology,Science and Applications,4:35-65

    Barsukova ML,Kuznetsov VA,Dorofeeva VA and Khodokovsky IL.1979.Experimental-study of solubility of Rutile-TiO2in fluoride solutions under High-temperature.Geokhimiya,(7):1017-1027(in Russian)

    Bayot D and Devillers M.2006.Peroxo complexes of niobium (V)and tantalum (V).Coordination Chemistry Reviews,250(19- 20):2610-2626

    Brenan JM,Shaw HF,Phinney DL and Ryerson FJ.1994.Rutileaqueous fluid partitioning of Nb,Ta,Hf,Zr,U and Th:Implications for High-field strength element depletions in Island-arc basalts.Earth and Planetary Science Letters,128(3- 4):327-339

    Bright E and Readey DW.1987.Dissolution kinetics of TiO2in HF-HCl solutions.Journal of the American Ceramic Society,70(12):900-906

    Charles RW and Vidale R.1982.Temperature calibration of a new rapid quench vessel.American Mineralogist,67(1-2):175-179

    Dalton JA and Gibb FGF.1996.Temperature gradients in large cold-seal pressure vessels.Mineralogical Magazine,60(399):337-345

    Ding X,Lundstrom C,Huang F,Li J,Zhang ZM,Sun XM,Liang JL and Sun WD.2009.Nature and experimental constraints on formation of the continental crust based on niobium-tantalum fractionation.International Geology Review,51(6):473-501

    Ding X,Sun WD,Wang FY,Chen LL,Li QL and Chen FK.2012.Single-grain mica Rb-Sr isochron ages and mineral chemistry for the Weishan pluton in Hunan Province and implications on petrogenesis and mineralization of Mesozoic composite granite in South China.Acta Petrologica Sinica,28(12):3823-3840 (in Chinese with English abstract)

    Ding X,Hu YH,Zhang H,Li CY,Ling MX and Sun WD.2013.Major Nb/Ta fractionation recorded in garnet amphibolite facies metagabbro.Journal of Geology,121(3):255-274

    Fatmi M,Chihi T,Ghebouli MA,Ghebouhi B,Derafa A,Sahnoune F,Belhouchet H and Gamoura A.2013.Study of microstructure,microhardness and phases of precipitation during ageing of Cu-Sn and Al-Ag alloys.American Journal of Mechanics and Applications,1(1):1-4

    Gaetani GA,Asimow PD and Stolper EM.2008.A model for rutile saturation in silicate melts with applications to eclogite partial melting in subduction zones and mantle plumes.Earth and Planetary Science Letters,272(3-4):720-729

    Gao J,John T,Klemd R and Xiong XM.2007.Mobilization of Ti-Nb-Ta during subduction: Evidence from rutile-bearing dehydration segregations and veins hosted in eclogite,Tianshan,NW China.Geochimica et Cosmochimica Acta,71(20):4974-4996

    Guo S,Ye K,Chen Y,Liu JB,Mao Q and Ma YG.2012.Fluid-rock interaction and element mobilization in UHP metabasalt:Constraints from an omphacite-epidote vein and host eclogites in the Dabie orogen.Lithos,136-139:145-167

    Hayden LA and Manning CE.2011.Rutile solubility in supercritical NaAlSi3O8-H2O fluids.Chemical Geology,284(1-2):74-81

    He JJ,Ding X,Wang YR,Sun WD and Fu B.2015.The effect of temperature and concentration on hydrolysis of fluorine-rich titanium complexes in hydrothermal fluids:Constraints on titanium mobility in deep geological processes.Acta Petrologica Sinica,31(3):802-810 (in Chinese with English abstract)

    Holloway JR,Dixon JE and Pawley AR.1992.An internally heated,rapid-quench,high-pressure vessel.American Mineralogist,77(5-6):643-646

    Huang J,Xiao YL,Gao Y,Hou Z and Wu W.2012.Nb-Ta fractionation induced by fluid-rock interaction in subduction-zones:Constraints from UHP eclogite-and vein-hosted rutile from the Dabie orogen,Central-Eastern China.Journal of Metamorphic Geology,30(8):821-842

    Jacobs GK and Kerrick DM.1980.A simple rapid-quench design for cold-seal pressure-vessels.American Mineralogist,65(9- 10):1053-1056

    Jalava JP,Heikkil? L,Hovi O,Laiho R,Hiltunen E,Hakanen A and Harma H.1998.Structural investigation of hydrous TiO2precipitates and their aging products by X-ray diffraction, atomic force microscopy,and transmission electron microscopy.Industrial &Engineering Chemistry Research,37(4):1317-1323

    Keppler H and Audétat A.2005.Fluid-mineral interaction at high pressure.European Mineralogical Union Notes in Mineralogy,7:225-251

    Knauss KG,Dibley MJ,Bourcier WL and Shaw HF.2001.Ti(IV)hydrolysis constants derived from rutile solubility measurements made from 100 to 300 degrees C.Applied Geochemistry,16(9-10):1115-1128

    Lee BG,Choi JW,Lee SE,Jeong YS,Oh HJ and Chi CS.2009.Formation behavior of anodic TiO2nanotubes in fluoride containing electrolytes.Transactions of Nonferrous Metals Society of China,19(4):842-845

    Liang JL,Ding X,Sun XM,Zhang ZM,Zhang H and Sun WD.2009.Nb/Ta fractionation observed in eclogites from the Chinese Continental Scientific Drilling Project.Chemical Geology,268(1-2):27-40

    Ling MX,Liu YL,Williams IS,Teng FZ,Yang XY,Ding X,Wei GJ,Xie LH,Deng WF and Sun WD.2013.Formation of the world’s largest REE deposit through protracted fluxing of carbonatite by subduction-derived fluids.Scientific Reports,3:doi:10.1038/srep01776

    Linnen RL,Samson IM,Williams-Jones AE and Chakhmouradian AR.2013.Geochemistry of the rare-earth element,Nb,Ta,Hf,and Zr deposits.In:Treatise on Geochemistry.2ndEdition.Elsevier Science,543-568

    Liu WH,Etschmann B,F(xiàn)oran G,Shelley M and Brugger J.2007.Deriving formation constants for aqueous metal complexes from XANES spectra:Zn2+and Fe2+chloride complexes in hypersaline solutions.American Mineralogist,92(5-6):761-770

    Meinhold G.2010.Rutile and its applications in earth sciences.Earth-Science Review,102(1-2):1-28

    Mysen B.2012.High-pressure and high-temperature titanium solution mechanisms in silicate-saturated aqueous fluids and hydrous silicate melts.American Mineralogist,97(7):1241-1251

    Namin HE,Hashemipour H and Ranjbar M.2008.Effect of aging and calcination on morphology and properties of synthesized nanocrystalline TiO2.International Journal of Modern Physics B,22(18-19):3210-3215

    Palacios E,Leret P,de la Mata MJ,F(xiàn)ernandez JF,de Aza AH and Rodriquez MA.2013.Influence of the pH and ageing time on the acid aluminum phosphate synthesized by precipitation.Cryst.Eng.Comm.,15(17):3359-3365

    Pearce JA and Peate DW.1995.Tectonic implications of the composition of volcanic arc magmas.Annual Review of Earth and Planetary Sciences,23:251-285

    Pottier AS,Cassaignon S,Chanéac C,Villain F,Tronc E and Jolivet JP.2003.Size tailoring of TiO2anatase nanoparticles in aqueous medium and synthesis of nanocomposites characterization by Raman spectroscopy.Journal of Materials Chemistry,13(4):877-882

    Rapp JF,Klemme S,Butler IB,and Harley SL.2010.Extremely high solubility of rutile in chloride and fluoride-bearing metamorphic fluids:An experimental investigation.Geology,38(4):323-326

    Ryerson FJ and Watson EB.1987.Rutile saturation in magmas:Implications for Ti-Nb-Ta depletion in island-arc basalts.Earth and Planetary Science Letters,86(2-4):225-239

    Ryzhenko BN,Kovalenko NI and Prisyagina NI.2006.Titanium complexation in hydrothermal systems.Geochemistry International,44(9):879-895

    Rudert V,Chou IM and Eugster HP.1976.Temperature-gradients in rapid-quench cold-seal pressure-vessels.American Mineralogist,61(9-10):1012-1015

    Schott J,Pokrovsky OS and Oelkers EH.2009.The link between mineral dissolution/precipitation kinetics and solution chemistry.Reviews in Mineralogy & Geochemistry,70(1):207-258

    Sileri D,Sahu KC and Matar OK.2011.Two-fluid pressure-driven channel flow with wall deposition and ageing effects.Journal of Engineering Mathematics,71(1):109-130

    Stefánsson,A.2007.Iron (III)hydrolysis and solubility at 25℃.Environ.Sci.Technol.,41(17):6117-6123

    Taylor JR and Wall VJ.1993.Cassiterite solubility,tin speciation,and transport in a magmatic aqueous phase.Economic Geology,88(2):437-460

    Tropper P and Manning CE.2005.Very low solubility of rutile in H2O at high pressure and temperature,and its implications for Ti mobility in subduction zones.American Mineralogist,90(2-3):502-505

    VanBaalen MR.1993.Titanium mobility in metamorphic systems:A review.Chemical Geology,110(1-3):233-249

    Wang YR.1978.Geochemical significance and experimental technique of complexes research at HTHP conditions.Chinese Science Bulletin,(11):682-686 (in Chinese)

    Wang YR.1995.The reversible and irreversible hydrolysis reaction of complexes under hydrothermal conditions.Bulletin of Mineralogy,Petrology and Geochemistry,(3):158-159 (in Chinese)

    Wang YR,Chen WH and Zhang HX.1995.Experimental study of irreversible hydrolysis reactions of Sn-,Nb-and Ta-fluorine complex at high temperature and high pressure.Geochimica,24(Suppl.):183-190 (in Chinese with English abstract)

    Wang YR and Chou IM.1986.Characteristics of hydrolysis of the complex Na2SnF6in hydrothermal solutions:An experimental study.Geochimica,(3):202-210 (in Chinese with English abstract)

    Wang YR,Gu F and Yuan ZQ.1993.Partitioning and hydrolysis of Nb and Ta and their implications with regard to mineralization.Chinese Journal of Geochemistry,12(1):84-91

    Williams-Jones AE,Bowell RJ and Migdisov AA.2009.Gold in solution.Elements,5(5):281-287

    Xiao YL,Sun WD,Hoefs J,Simon K,Zhang ZM,Li SG and Hofmann AW.2006.Making continental crust through slab melting:Constraints from niobium-tantalum fractionation in UHP metamorphic rutile.Geochimica et Cosmochimica Acta,70(18):4770-4782

    Yardley BWD and Bodnar RJ.2014.Fluids in the continental crust.Geochemical Perspectives,3(1):1-127

    附中文參考文獻(xiàn)

    丁興,孫衛(wèi)東,汪方躍,陳林麗,李秋立,陳福坤.2012.湖南溈山巖體多期云母的Rb-Sr 同位素年齡和礦物化學(xué)組成及其成巖成礦指示意義.巖石學(xué)報(bào),28(12):3823-3840

    何俊杰,丁興,王玉榮,孫衛(wèi)東,傅斌.2015.溫度、濃度對(duì)流體中氟鈦絡(luò)合物水解的影響:對(duì)深部地質(zhì)過(guò)程中鈦元素活動(dòng)的制約.巖石學(xué)報(bào),31(3):802-810

    王玉榮.1978.高溫高壓下絡(luò)合物研究的地球化學(xué)意義與實(shí)驗(yàn)方法.科學(xué)通報(bào),(11):682-686

    王玉榮.1995.水熱條件下配合物水解及可逆與不可逆反應(yīng).礦物巖石地球化學(xué)通報(bào),(3):158-159

    王玉榮,陳文華,張海祥.1995.Sn、Nb、Ta 氟配合物高溫水解不可逆反應(yīng)實(shí)驗(yàn)研究.地球化學(xué),24(增刊):183-190

    王玉榮,周義明.1986.Na2SnF6在水熱溶液中水解特征的實(shí)驗(yàn)研究.地球化學(xué),(3):202-210

    猜你喜歡
    沉淀物陳化淬火
    淬火礪劍
    ON THE EXISTENCE WITH EXPONENTIAL DECAY AND THE BLOW-UP OF SOLUTIONS FOR COUPLED SYSTEMS OF SEMI-LINEAR CORNER-DEGENERATE PARABOLIC EQUATIONS WITH SINGULAR POTENTIALS?
    六堡茶陳化工藝優(yōu)化研究
    高濃度硫酸體系中鈾與鐵的HDEHP萃取分離
    濕法冶金(2018年6期)2018-12-13 09:19:58
    胸腔積液沉淀物在惡性胸腔積液診斷中的臨床應(yīng)用價(jià)值
    淬火的目的
    淬火
    沉淀物對(duì)安塞油田原油計(jì)量的影響及應(yīng)對(duì)措施
    新型高頻淬火及表面強(qiáng)化技術(shù)與應(yīng)用
    應(yīng)用16S rDNA克隆文庫(kù)技術(shù)分析陳化煙葉細(xì)菌多樣性
    亚洲欧美成人精品一区二区| 99久国产av精品国产电影| av一本久久久久| 久久久成人免费电影| 国产高清不卡午夜福利| 亚洲精品一二三| 国产精品国产av在线观看| 高清毛片免费看| 99热全是精品| 美女内射精品一级片tv| 狂野欧美激情性xxxx在线观看| 国产成人91sexporn| 两个人的视频大全免费| 18+在线观看网站| av专区在线播放| 涩涩av久久男人的天堂| 久久国产精品男人的天堂亚洲 | 午夜精品国产一区二区电影| 国产精品熟女久久久久浪| 人人妻人人爽人人添夜夜欢视频 | 夫妻性生交免费视频一级片| 久久久亚洲精品成人影院| 亚洲国产色片| 男人舔奶头视频| 亚洲色图av天堂| 嘟嘟电影网在线观看| 国产成人a区在线观看| 久热这里只有精品99| 国产精品国产av在线观看| 涩涩av久久男人的天堂| 亚洲人成网站高清观看| 乱码一卡2卡4卡精品| 国产精品麻豆人妻色哟哟久久| 青春草视频在线免费观看| 亚洲精品第二区| 免费在线观看成人毛片| 中文乱码字字幕精品一区二区三区| 女性被躁到高潮视频| 男人爽女人下面视频在线观看| 能在线免费看毛片的网站| 少妇熟女欧美另类| 国产深夜福利视频在线观看| 精品一品国产午夜福利视频| 日本黄色日本黄色录像| 亚洲国产日韩一区二区| 在线天堂最新版资源| 永久免费av网站大全| 精品人妻视频免费看| 男女免费视频国产| 欧美三级亚洲精品| 国产成人精品久久久久久| 大话2 男鬼变身卡| 少妇的逼水好多| 中文资源天堂在线| 伊人久久精品亚洲午夜| 国产免费又黄又爽又色| 熟妇人妻不卡中文字幕| 国产91av在线免费观看| 高清毛片免费看| 国产高清三级在线| 亚洲国产色片| 欧美zozozo另类| 丝袜脚勾引网站| 欧美丝袜亚洲另类| 国产精品人妻久久久影院| 精品熟女少妇av免费看| 国产亚洲av片在线观看秒播厂| 高清在线视频一区二区三区| 久久精品国产亚洲av天美| 人妻少妇偷人精品九色| 女性被躁到高潮视频| av卡一久久| 日韩成人av中文字幕在线观看| 国产精品一及| 高清黄色对白视频在线免费看 | 久久久午夜欧美精品| av不卡在线播放| 国产av国产精品国产| 黄色怎么调成土黄色| 久久韩国三级中文字幕| 在线免费观看不下载黄p国产| 婷婷色av中文字幕| 97在线人人人人妻| 久久久久久久亚洲中文字幕| 老女人水多毛片| 亚洲欧美成人精品一区二区| 亚洲欧美成人精品一区二区| 成人高潮视频无遮挡免费网站| 身体一侧抽搐| 国产黄频视频在线观看| 久久99蜜桃精品久久| 欧美日韩在线观看h| 99热网站在线观看| 爱豆传媒免费全集在线观看| 亚洲国产欧美在线一区| 久久久成人免费电影| 交换朋友夫妻互换小说| 日产精品乱码卡一卡2卡三| 97超视频在线观看视频| 日本猛色少妇xxxxx猛交久久| 黑人猛操日本美女一级片| 十分钟在线观看高清视频www | 亚洲欧美精品专区久久| 日韩不卡一区二区三区视频在线| 一级毛片 在线播放| 精品一区二区免费观看| 婷婷色综合大香蕉| 黑丝袜美女国产一区| 在线播放无遮挡| 成年av动漫网址| 最近手机中文字幕大全| 国产男女内射视频| 亚洲av福利一区| 97在线视频观看| 国产精品一区二区在线不卡| 成人影院久久| 男女边吃奶边做爰视频| 国产精品一区二区在线不卡| 精品少妇久久久久久888优播| 舔av片在线| 18禁在线无遮挡免费观看视频| 国产精品熟女久久久久浪| 成人美女网站在线观看视频| 中文天堂在线官网| 免费大片18禁| 精品亚洲成国产av| 免费高清在线观看视频在线观看| 嫩草影院新地址| 直男gayav资源| av天堂中文字幕网| 免费高清在线观看视频在线观看| 日韩不卡一区二区三区视频在线| 国产精品精品国产色婷婷| 大片免费播放器 马上看| 国精品久久久久久国模美| 国产乱来视频区| 国产欧美日韩一区二区三区在线 | 免费av中文字幕在线| 舔av片在线| 免费观看的影片在线观看| 蜜桃亚洲精品一区二区三区| 久久婷婷青草| 精品酒店卫生间| 日本黄大片高清| 亚洲欧美日韩无卡精品| 又粗又硬又长又爽又黄的视频| 97超视频在线观看视频| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜爱| 一本久久精品| 日本黄色片子视频| 亚洲丝袜综合中文字幕| 国产一区二区三区综合在线观看 | 午夜福利在线观看免费完整高清在| 你懂的网址亚洲精品在线观看| 一区二区三区精品91| 国产精品三级大全| 中国国产av一级| 一级毛片久久久久久久久女| 美女福利国产在线 | 伊人久久精品亚洲午夜| 一级毛片 在线播放| 久久精品夜色国产| 国产熟女欧美一区二区| 深爱激情五月婷婷| 日本欧美视频一区| 精品久久久久久久久亚洲| 春色校园在线视频观看| 国产精品国产av在线观看| 26uuu在线亚洲综合色| 日韩免费高清中文字幕av| 少妇的逼水好多| 少妇丰满av| 国产免费一区二区三区四区乱码| 最近2019中文字幕mv第一页| 国产精品久久久久久久久免| 国产高清三级在线| 亚洲无线观看免费| 日本-黄色视频高清免费观看| 日韩av不卡免费在线播放| 五月玫瑰六月丁香| 午夜免费男女啪啪视频观看| 日日摸夜夜添夜夜添av毛片| 永久网站在线| 国产精品国产三级国产av玫瑰| 国产成人a∨麻豆精品| 国产精品蜜桃在线观看| 日韩 亚洲 欧美在线| 亚洲欧美中文字幕日韩二区| 又黄又爽又刺激的免费视频.| 欧美日本视频| 色5月婷婷丁香| 小蜜桃在线观看免费完整版高清| 日韩欧美 国产精品| 欧美最新免费一区二区三区| 赤兔流量卡办理| 狠狠精品人妻久久久久久综合| 国产色爽女视频免费观看| 免费黄频网站在线观看国产| 18禁动态无遮挡网站| 五月天丁香电影| 老司机影院成人| 18禁裸乳无遮挡免费网站照片| 久久国内精品自在自线图片| 久久人人爽人人片av| 97在线人人人人妻| 91精品一卡2卡3卡4卡| av福利片在线观看| 日本vs欧美在线观看视频 | 偷拍熟女少妇极品色| 国产视频内射| 日韩一区二区视频免费看| 精品熟女少妇av免费看| 国产精品国产三级国产av玫瑰| 欧美zozozo另类| 成年免费大片在线观看| 老女人水多毛片| 女人久久www免费人成看片| 在线观看免费日韩欧美大片 | 日日摸夜夜添夜夜添av毛片| a级毛色黄片| 精品视频人人做人人爽| 欧美最新免费一区二区三区| av又黄又爽大尺度在线免费看| 少妇熟女欧美另类| 久久女婷五月综合色啪小说| 国产精品人妻久久久影院| 久久韩国三级中文字幕| av在线老鸭窝| 身体一侧抽搐| 七月丁香在线播放| 国产精品女同一区二区软件| 欧美xxxx性猛交bbbb| 久久这里有精品视频免费| 午夜免费男女啪啪视频观看| 内地一区二区视频在线| 午夜免费男女啪啪视频观看| 亚洲精品视频女| 七月丁香在线播放| 国产精品一区二区在线不卡| 国产美女午夜福利| 男女无遮挡免费网站观看| 国产亚洲5aaaaa淫片| 99久久精品国产国产毛片| 久久鲁丝午夜福利片| 日本猛色少妇xxxxx猛交久久| 丝瓜视频免费看黄片| 男女啪啪激烈高潮av片| 日本猛色少妇xxxxx猛交久久| 在现免费观看毛片| 国产成人a区在线观看| 男女免费视频国产| 多毛熟女@视频| a级毛色黄片| 久久鲁丝午夜福利片| 国产一区二区在线观看日韩| 一级二级三级毛片免费看| 一级黄片播放器| 日日摸夜夜添夜夜添av毛片| 婷婷色综合大香蕉| 国产精品三级大全| 18禁在线无遮挡免费观看视频| 一级av片app| 国产成人一区二区在线| 久久国产乱子免费精品| 亚洲自偷自拍三级| 亚洲综合色惰| 国产成人aa在线观看| 亚洲,欧美,日韩| 菩萨蛮人人尽说江南好唐韦庄| 久久久久网色| 国产av一区二区精品久久 | av一本久久久久| 中文字幕制服av| 一级毛片黄色毛片免费观看视频| 国产精品欧美亚洲77777| 最近的中文字幕免费完整| 日本黄色日本黄色录像| 18禁裸乳无遮挡动漫免费视频| 一级黄片播放器| 亚洲,欧美,日韩| 高清黄色对白视频在线免费看 | 九色成人免费人妻av| 国产精品一区www在线观看| 成人一区二区视频在线观看| 欧美三级亚洲精品| 久久精品久久精品一区二区三区| 国产精品99久久久久久久久| 一级片'在线观看视频| 美女高潮的动态| 婷婷色综合www| 一本一本综合久久| 老女人水多毛片| 青春草视频在线免费观看| 国产成人午夜福利电影在线观看| 免费看光身美女| 欧美最新免费一区二区三区| 国产av码专区亚洲av| 国产高清不卡午夜福利| 五月天丁香电影| 少妇 在线观看| 亚洲av综合色区一区| 亚洲美女搞黄在线观看| 欧美高清性xxxxhd video| 亚洲av在线观看美女高潮| 欧美激情国产日韩精品一区| 1000部很黄的大片| 亚洲欧美日韩无卡精品| 午夜免费鲁丝| 亚洲国产成人一精品久久久| 欧美最新免费一区二区三区| 色婷婷av一区二区三区视频| 久久久色成人| 国产成人午夜福利电影在线观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美成人精品一区二区| 寂寞人妻少妇视频99o| 麻豆国产97在线/欧美| 亚洲人与动物交配视频| 高清在线视频一区二区三区| 精品亚洲乱码少妇综合久久| 午夜福利在线观看免费完整高清在| 国产男女内射视频| 欧美一区二区亚洲| 欧美日韩在线观看h| 久久青草综合色| 99久久人妻综合| 亚洲欧美精品自产自拍| 久久久久人妻精品一区果冻| 丝袜喷水一区| 精品久久久久久久久亚洲| 在线免费十八禁| 大陆偷拍与自拍| 搡老乐熟女国产| 少妇精品久久久久久久| 婷婷色麻豆天堂久久| 一本久久精品| 亚洲国产精品999| 久久久久久伊人网av| 欧美xxxx性猛交bbbb| 伦理电影大哥的女人| 麻豆国产97在线/欧美| 亚洲图色成人| 亚洲精品亚洲一区二区| 我要看黄色一级片免费的| 18+在线观看网站| 亚洲精品日韩av片在线观看| 久久久色成人| 人人妻人人爽人人添夜夜欢视频 | 色综合色国产| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 久久久久久九九精品二区国产| 观看美女的网站| 中文在线观看免费www的网站| 汤姆久久久久久久影院中文字幕| 亚洲aⅴ乱码一区二区在线播放| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线| 99国产精品免费福利视频| 国产亚洲91精品色在线| av线在线观看网站| 国产精品成人在线| 亚洲欧美一区二区三区国产| 国产一区二区三区综合在线观看 | 久久精品国产亚洲网站| 看免费成人av毛片| av一本久久久久| 新久久久久国产一级毛片| 日本av免费视频播放| 水蜜桃什么品种好| 国产av码专区亚洲av| 亚洲色图综合在线观看| 夜夜爽夜夜爽视频| 欧美日韩在线观看h| 国内精品宾馆在线| 国产在线男女| 亚洲美女黄色视频免费看| av国产精品久久久久影院| 美女中出高潮动态图| 国产一区二区在线观看日韩| 久久精品夜色国产| 国产精品.久久久| 精品一区在线观看国产| 亚洲精品乱码久久久v下载方式| 全区人妻精品视频| 天美传媒精品一区二区| 精品一品国产午夜福利视频| 少妇人妻一区二区三区视频| 一区二区三区乱码不卡18| 国产在线男女| 看非洲黑人一级黄片| 乱系列少妇在线播放| 岛国毛片在线播放| 日本猛色少妇xxxxx猛交久久| 国产爽快片一区二区三区| 亚洲一区二区三区欧美精品| 免费观看的影片在线观看| 夜夜看夜夜爽夜夜摸| 色婷婷久久久亚洲欧美| 亚洲av国产av综合av卡| 超碰97精品在线观看| 国产 一区精品| 欧美日韩亚洲高清精品| 亚洲伊人久久精品综合| 亚洲三级黄色毛片| 黑人猛操日本美女一级片| 中文资源天堂在线| av国产精品久久久久影院| 网址你懂的国产日韩在线| 一级毛片我不卡| 日本免费在线观看一区| 在线免费十八禁| freevideosex欧美| 亚洲国产成人一精品久久久| 色吧在线观看| 国产av一区二区精品久久 | 一个人免费看片子| 精品久久久久久久久av| 欧美+日韩+精品| 91久久精品国产一区二区成人| 成人亚洲欧美一区二区av| 久久国内精品自在自线图片| 日韩av在线免费看完整版不卡| 日本黄色片子视频| 亚洲av成人精品一二三区| 老女人水多毛片| 免费久久久久久久精品成人欧美视频 | 极品教师在线视频| 亚洲国产精品专区欧美| 大又大粗又爽又黄少妇毛片口| 男女边吃奶边做爰视频| 精品少妇黑人巨大在线播放| 一区二区三区免费毛片| 国产精品一区二区性色av| 国产成人aa在线观看| 国产精品成人在线| 久久精品国产鲁丝片午夜精品| 大码成人一级视频| 久久亚洲国产成人精品v| 校园人妻丝袜中文字幕| 成人亚洲欧美一区二区av| 在现免费观看毛片| 国产精品偷伦视频观看了| 特大巨黑吊av在线直播| 97超视频在线观看视频| 午夜免费男女啪啪视频观看| 高清不卡的av网站| 日本黄大片高清| 久久精品夜色国产| 26uuu在线亚洲综合色| 老女人水多毛片| 欧美国产精品一级二级三级 | 婷婷色综合大香蕉| 久久久a久久爽久久v久久| 久久久久精品久久久久真实原创| 久久精品国产亚洲av天美| 亚洲自偷自拍三级| 日日摸夜夜添夜夜添av毛片| 99久久综合免费| 天堂中文最新版在线下载| 一级毛片我不卡| 22中文网久久字幕| 亚洲av不卡在线观看| 色吧在线观看| 91精品一卡2卡3卡4卡| 国产免费视频播放在线视频| 男女免费视频国产| 久久精品久久久久久久性| 99久久精品热视频| 秋霞伦理黄片| 国产男女内射视频| 大又大粗又爽又黄少妇毛片口| 婷婷色综合大香蕉| 亚洲av成人精品一二三区| 国产黄片美女视频| 啦啦啦视频在线资源免费观看| 成人免费观看视频高清| 色婷婷久久久亚洲欧美| 深夜a级毛片| 乱码一卡2卡4卡精品| 丰满乱子伦码专区| 国产在线免费精品| 免费观看的影片在线观看| 这个男人来自地球电影免费观看 | 18禁在线无遮挡免费观看视频| 国产女主播在线喷水免费视频网站| 欧美丝袜亚洲另类| 黄片无遮挡物在线观看| 精品人妻一区二区三区麻豆| 亚洲精品国产成人久久av| 免费大片黄手机在线观看| 美女内射精品一级片tv| 黑丝袜美女国产一区| 一区二区三区乱码不卡18| 蜜桃久久精品国产亚洲av| 九九在线视频观看精品| 国模一区二区三区四区视频| 九九在线视频观看精品| 大香蕉97超碰在线| 亚洲精品aⅴ在线观看| 中文乱码字字幕精品一区二区三区| 午夜福利网站1000一区二区三区| 乱系列少妇在线播放| 免费高清在线观看视频在线观看| 啦啦啦啦在线视频资源| 在线观看免费日韩欧美大片 | 五月伊人婷婷丁香| 亚洲精品456在线播放app| 国产精品成人在线| 狂野欧美激情性bbbbbb| 91午夜精品亚洲一区二区三区| 水蜜桃什么品种好| 交换朋友夫妻互换小说| 能在线免费看毛片的网站| 亚洲成人av在线免费| 国产黄片视频在线免费观看| 精品久久久噜噜| 九九久久精品国产亚洲av麻豆| 91精品国产九色| 亚洲欧美精品专区久久| 亚洲av综合色区一区| 中文欧美无线码| 国产精品国产三级国产专区5o| 交换朋友夫妻互换小说| 久久午夜福利片| 五月开心婷婷网| 校园人妻丝袜中文字幕| 日韩成人伦理影院| av国产精品久久久久影院| 国产一区二区三区综合在线观看 | 天堂俺去俺来也www色官网| freevideosex欧美| 最近的中文字幕免费完整| 久久精品人妻少妇| 韩国av在线不卡| 自拍偷自拍亚洲精品老妇| 亚洲精品456在线播放app| 国产伦精品一区二区三区视频9| 又爽又黄a免费视频| 日日撸夜夜添| 国产日韩欧美在线精品| 激情五月婷婷亚洲| 欧美激情国产日韩精品一区| 国产成人免费无遮挡视频| 免费av不卡在线播放| 亚洲va在线va天堂va国产| 午夜福利在线观看免费完整高清在| 久久久久精品性色| 边亲边吃奶的免费视频| 天堂俺去俺来也www色官网| 97超视频在线观看视频| 国产成人freesex在线| 在线免费十八禁| 搡女人真爽免费视频火全软件| 午夜免费男女啪啪视频观看| 精品久久久精品久久久| 日本与韩国留学比较| 国产成人精品一,二区| 国产av国产精品国产| 天天躁夜夜躁狠狠久久av| 久久精品国产鲁丝片午夜精品| 亚洲av电影在线观看一区二区三区| 熟女电影av网| 国产精品伦人一区二区| 有码 亚洲区| 看十八女毛片水多多多| 久久精品国产亚洲网站| 国产精品人妻久久久影院| 国产在线免费精品| 亚洲精品成人av观看孕妇| 日日撸夜夜添| 十八禁网站网址无遮挡 | 久久久精品免费免费高清| 国产免费一区二区三区四区乱码| 国产精品福利在线免费观看| 欧美97在线视频| 国产精品久久久久久av不卡| 只有这里有精品99| 日韩视频在线欧美| 黄色怎么调成土黄色| 国产永久视频网站| 又爽又黄a免费视频| 色视频在线一区二区三区| 中文天堂在线官网| 高清在线视频一区二区三区| 人人妻人人添人人爽欧美一区卜 | 我要看黄色一级片免费的| 美女主播在线视频| 人人妻人人看人人澡| 国产亚洲精品久久久com| 26uuu在线亚洲综合色| 蜜桃在线观看..| 男人和女人高潮做爰伦理| 免费播放大片免费观看视频在线观看| 日韩在线高清观看一区二区三区| 亚洲四区av| 视频区图区小说| 欧美日本视频| 麻豆成人午夜福利视频| 久久久久久伊人网av| 国产成人freesex在线| 久久精品夜色国产| 日韩中文字幕视频在线看片 | 乱系列少妇在线播放| 99久久中文字幕三级久久日本| 国产欧美亚洲国产| 日韩一区二区三区影片| 大片电影免费在线观看免费| 黄色日韩在线| 国产69精品久久久久777片| 国产精品免费大片| 国产熟女欧美一区二区| 精品99又大又爽又粗少妇毛片|