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

    水蒸氣對(duì)煅燒石灰石碳酸化反應(yīng)影響的實(shí)驗(yàn)與模型分析

    2015-10-13 12:43:15劉洋楊勇平
    化工學(xué)報(bào) 2015年3期
    關(guān)鍵詞:碳酸水蒸氣轉(zhuǎn)化率

    劉洋,楊勇平

    ?

    水蒸氣對(duì)煅燒石灰石碳酸化反應(yīng)影響的實(shí)驗(yàn)與模型分析

    劉洋,楊勇平

    (華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院,北京102206)

    利用熱重分析儀,對(duì)比了不同反應(yīng)溫度、不同水蒸氣濃度對(duì)煅燒石灰石碳酸化反應(yīng)的影響。碳酸化反應(yīng)溫度在500℃時(shí),反應(yīng)初期水蒸氣對(duì)碳酸化反應(yīng)的影響并不明顯,反應(yīng)10 min之后,在含有1.5%、10%和20%(體積分?jǐn)?shù))水蒸氣條件下碳酸化轉(zhuǎn)化率比無水蒸氣條件下轉(zhuǎn)化率分別提高了19.8%、27.2%和30.5%。水蒸氣的存在有助于提高碳酸化反應(yīng)轉(zhuǎn)化率,但隨著水蒸氣濃度的增加轉(zhuǎn)化率增加幅度減小。利用隨機(jī)孔隙模型,對(duì)產(chǎn)物層擴(kuò)散階段擴(kuò)散系數(shù)及反應(yīng)活化能進(jìn)行了計(jì)算。高溫條件下,氣氛中含有1.5%的水蒸氣使反應(yīng)活化能從237.7 kJ·mol-1降低到179.9 kJ·mol-1,提高水蒸氣濃度到10%和20%后,反應(yīng)活化能從156.6 kJ·mol-1降低到148.6 kJ·mol-1。不同水蒸氣濃度條件下,碳酸化反應(yīng)存在兩個(gè)明顯特征:一是大約在550℃處存在一個(gè)明顯擴(kuò)散系數(shù)的斜率變化,這一溫度與氣氛中是否存在水蒸氣無關(guān);另一特征是隨著反應(yīng)溫度的提高,水蒸氣的促進(jìn)作用減弱。依據(jù)實(shí)驗(yàn)和模型計(jì)算結(jié)果,推測(cè)了當(dāng)反應(yīng)處于產(chǎn)物層擴(kuò)散階段時(shí)水蒸氣對(duì)碳酸化反應(yīng)影響的作用機(jī)理。

    二氧化碳捕集;化學(xué)反應(yīng);催化;產(chǎn)物層擴(kuò)散系數(shù);活化能

    引 言

    煅燒后的石灰石與CO2的碳酸化反應(yīng)是煤氣化制氫、吸收增強(qiáng)式甲烷重整、吸收增強(qiáng)式水氣變換過程中的重要反應(yīng)步驟,也可以應(yīng)用于煙氣中CO2捕集等方面,因此受到了國內(nèi)外研究人員的廣泛關(guān)注。反應(yīng)式為

    CaO(s)+CO2(g)CaCO3(s) (1)

    近年來國內(nèi)外學(xué)者分別利用固定床或者流化床反應(yīng)器,對(duì)比了碳酸化反應(yīng)過程中水蒸氣對(duì)碳酸化反應(yīng)的影響,發(fā)現(xiàn)添加水蒸氣后可明顯提高吸收劑吸收CO2的能力[1-3]。這些實(shí)驗(yàn)結(jié)果表明,水蒸氣可以顯著提高鈣基吸收劑循環(huán)吸收CO2的能力。因此可以采用添加水蒸氣的方式提高吸收劑的吸收能力,無須對(duì)吸收劑本身進(jìn)行預(yù)處理或改性,降低了吸收CO2的成本。在吸收增強(qiáng)式甲烷重整、吸收增強(qiáng)式水氣變換中,由于存在一定含量的水蒸氣,而水蒸氣對(duì)碳酸化反應(yīng)的促進(jìn)作用增強(qiáng),有利于在更寬的范圍內(nèi)選擇反應(yīng)溫度。

    對(duì)于水蒸氣對(duì)碳酸化反應(yīng)的促進(jìn)作用存在不同的機(jī)理認(rèn)識(shí):最早在20世紀(jì)70年代Dobner等[4]研究了在水蒸氣存在條件下煅燒白云石的碳酸化反應(yīng),發(fā)現(xiàn)水蒸氣具有催化CaO與CO2反應(yīng)的能力,尤其在低溫條件下(~550℃)更為明顯。Wang等[5]認(rèn)為瞬間形成的Ca(OH)2是水蒸氣促進(jìn)碳酸化反應(yīng)的原因,而Yang等[6]通過比較Ca(OH)2與CO2和CaO與CO2的反應(yīng)速率認(rèn)為Ca(OH)2不是加快碳酸化反應(yīng)的原因。Manovic等[7-8]根據(jù)實(shí)驗(yàn)結(jié)果否認(rèn)水蒸氣的催化作用,而是認(rèn)為其作用在于提高擴(kuò)散階段反應(yīng)速率。從目前發(fā)表的相關(guān)文獻(xiàn)中可以看出,研究人員對(duì)水蒸氣對(duì)碳酸化反應(yīng)機(jī)理的認(rèn)識(shí)存在明顯的分歧。

    為確定水蒸氣對(duì)碳酸化反應(yīng)的作用規(guī)律,還需要選擇合適的反應(yīng)模型。目前水蒸氣對(duì)碳酸化反應(yīng)影響的研究大都為定性分析,缺乏動(dòng)力學(xué)參數(shù)對(duì)比,為準(zhǔn)確反映水蒸氣的作用規(guī)律,必須選擇合適的反應(yīng)模型以獲取動(dòng)力學(xué)參數(shù)。與未反應(yīng)縮核模型[9]和微粒模型[10]相比,孔隙模型可以準(zhǔn)確反映碳酸化反應(yīng)界面相互交叉,孔隙結(jié)構(gòu)同時(shí)發(fā)生變化這一過程。而由Bhatia等[11-12]和Gavalas[13]同時(shí)提出的隨機(jī)孔隙模型認(rèn)為吸收劑顆粒內(nèi)部由一系列分配特征的孔隙隨機(jī)交叉所組成,可以利用結(jié)構(gòu)參數(shù)計(jì)算化學(xué)反應(yīng)速率及產(chǎn)物層擴(kuò)散速率。本文采用隨機(jī)孔隙模型分析水蒸氣對(duì)碳酸化反應(yīng)的作用規(guī)律,重點(diǎn)計(jì)算產(chǎn)物層擴(kuò)散速率。在采用隨機(jī)孔隙模型時(shí),反應(yīng)體系中不應(yīng)包含除碳酸化以外的其他反應(yīng)。根據(jù)熱力學(xué)分析,當(dāng)水蒸氣濃度低于平衡濃度時(shí),并不會(huì)與CaO反應(yīng)[14]。此時(shí)雖然存在水蒸氣,然而在宏觀層面上只有CaO與CO2的碳酸化反應(yīng)發(fā)生。在此條件下,利用隨機(jī)孔隙模型研究CaO與包含CO2和水蒸氣的混合氣體反應(yīng),建立水蒸氣濃度與碳酸化產(chǎn)物層擴(kuò)散速率的定量關(guān)系才有其合理性。

    本文采用天然石灰石,利用熱重分析儀研究在不同溫度條件下,不同水蒸氣濃度對(duì)碳酸化反應(yīng)的影響規(guī)律。利用壓汞實(shí)驗(yàn),通過計(jì)算確定樣品碳酸化反應(yīng)初始條件下單位體積孔長度、比表面積和孔隙率。利用隨機(jī)孔隙模型,分析碳酸化反應(yīng)實(shí)驗(yàn)數(shù)據(jù)。通過比較不同水蒸氣濃度條件下產(chǎn)物層厚度、產(chǎn)物層擴(kuò)散系數(shù)以及活化能大小,確定水蒸氣對(duì)碳酸化反應(yīng)作用階段及作用規(guī)律,本文重點(diǎn)分析擴(kuò)散階段反應(yīng)速率。

    1 計(jì)算模型

    考慮球形煅燒后石灰石與氣體反應(yīng)物CO2發(fā)生氣固反應(yīng),假定固體反應(yīng)物內(nèi)部由一致的孔隙組成。反應(yīng)初始發(fā)生在孔邊界上,導(dǎo)致顆粒內(nèi)部徑向表面積的增長。反應(yīng)產(chǎn)物CaCO3聚集在反應(yīng)表面和孔邊界之間,阻礙了氣體反應(yīng)物進(jìn)一步向固體反應(yīng)物的擴(kuò)散。在反應(yīng)進(jìn)行中,臨近的表面積相互交叉,伴隨著反應(yīng)表面積逐漸減少。

    假定反應(yīng)表面是由一系列具有孔徑分配()的圓柱形孔的隨機(jī)交叉組成,()d是單位體積空間半徑在和d之間圓柱表面積的總長度。在任何半徑位置增長的圓柱形孔的孔徑分配存在如下關(guān)系

    由于碳酸化反應(yīng)是可逆反應(yīng),反應(yīng)速率是反應(yīng)氣體A濃度與平衡濃度之差的一階函數(shù)并正比于反應(yīng)表面積,則有

    并且

    式中,是氣體反應(yīng)物CO2的濃度,e是CO2平衡濃度,假定在顆粒內(nèi)部任意半徑位置的反應(yīng)表面上濃度都相同。-e代表了氣體A在整個(gè)反應(yīng)表面的有效濃度。

    基于上面的關(guān)系,圓柱形孔的總長度0、表面積0和體積0適用于確定反應(yīng)速率,可以得到

    結(jié)構(gòu)參數(shù)由式(6)定義

    其中

    定義反應(yīng)過程形成的產(chǎn)物體積與反應(yīng)物體積之比是,則有

    -=(-0) (9)

    結(jié)合式(4)和式(9)得到形成反應(yīng)物的速率為

    定義為擴(kuò)散產(chǎn)物層平均厚度,并且認(rèn)為遠(yuǎn)小于反應(yīng)表面尺寸

    結(jié)合式(10)和式(11)可以得到

    由式(5)和式(12)可以得到

    考慮到當(dāng)初始=0時(shí),=0,積分得到

    式(14)表達(dá)了在任意轉(zhuǎn)化率時(shí),有效產(chǎn)物層厚度??紤]到反應(yīng)氣體A的質(zhì)量平衡,有

    結(jié)合式(5)、式(14)和式(15),可以得到反應(yīng)速率表達(dá)式

    在產(chǎn)物層擴(kuò)散階段,由于→∞,因此式(16)可以積分為

    其中

    2 實(shí) 驗(yàn)

    2.1 熱重實(shí)驗(yàn)

    實(shí)驗(yàn)采用妙峰山石灰石(MFS),煅燒溫度為900℃,碳酸化溫度為400~700℃。每次實(shí)驗(yàn)時(shí),質(zhì)量約為10 mg的石灰石樣品平鋪在樣品盤內(nèi)。CO2氣體體積分?jǐn)?shù)為20%,其余氣體為氮?dú)狻2捎脤O2氣體通過洗氣瓶的方式添加水蒸氣,將洗氣瓶放入恒溫油浴鍋中,通過控制油浴溫度控制攜帶水蒸氣含量。設(shè)定油浴溫度分別為20、45、60℃,對(duì)應(yīng)的氣體中水蒸氣體積分?jǐn)?shù)為1.5%、10%和20%。為防止水蒸氣在管路中冷凝,從洗氣瓶出口到TGA進(jìn)口纏上伴熱帶,溫度控制在100℃。在進(jìn)行碳酸化反應(yīng)之前,提前30 min將CO2和水蒸氣混合,通過旁路排空以保證水蒸氣濃度準(zhǔn)確。碳酸化反應(yīng)結(jié)束后,將CO2氣體切入旁路,TGA通入氮?dú)馀趴辗磻?yīng)爐中殘存的水蒸氣以確保安全。

    TGA實(shí)驗(yàn)同時(shí)得到時(shí)間、質(zhì)量和溫度3組數(shù)據(jù)。在碳酸化反應(yīng)過程中,樣品與通入的氣體反應(yīng)導(dǎo)致質(zhì)量增加,從增重點(diǎn)開始進(jìn)行轉(zhuǎn)化率的計(jì)算

    式中,0表示樣品煅燒后的原始質(zhì)量,g;p為樣品純度,樣品中的雜質(zhì)不參與碳酸化反應(yīng);m為時(shí)刻樣品的質(zhì)量,g。

    2.2 鈣基吸收劑

    MFS石灰石在破碎機(jī)中破碎,然后用標(biāo)準(zhǔn)篩篩分成粒徑為100~200mm的顆粒,采用X射線熒光光譜儀(型號(hào)為Rigaku ZSX Primus Π)測(cè)量樣品成分。樣品經(jīng)過煅燒后,利用壓汞儀(型號(hào)為AutoPore IV9510)測(cè)量孔隙結(jié)構(gòu)。表1列出了樣品成分質(zhì)量分?jǐn)?shù)組成及孔隙結(jié)構(gòu)參數(shù)。

    表1 石灰石成分組成及煅燒石灰石孔結(jié)構(gòu)參數(shù)Table 1 Compositional characteristics of limestone and pore structure parameters of calcined limestone

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

    當(dāng)氣氛中存在CO2和水蒸氣時(shí),煅燒的石灰石會(huì)與高于平衡分壓的CO2和水蒸氣分別反應(yīng)生成CaCO3和Ca(OH)2。圖1顯示了利用HSC軟件計(jì)算的碳酸化反應(yīng)和活化反應(yīng)氣體平衡分壓結(jié)果。利用軟件計(jì)算發(fā)現(xiàn),在溫度650℃時(shí),使CaO反應(yīng)產(chǎn)生氫氧化鈣的水蒸氣壓力高于一個(gè)大氣壓。因此在碳酸化反應(yīng)通常所選取的這一溫度,無論水蒸氣濃度高低,在常壓下不會(huì)有氫氧化鈣的形成。而在400℃和450℃時(shí),利用軟件計(jì)算得到的水蒸氣平衡濃度約為8.2%和19.3%。從熱力學(xué)平衡分壓角度看,水蒸氣壓力高于平衡分壓時(shí)可以反應(yīng)生成Ca(OH)2,低于平衡分壓時(shí)不生成Ca(OH)2。在同時(shí)存在CO2和高于平衡分壓的水蒸氣的條件下,CaO可以與CO2反應(yīng)生成CaCO3,也可以先與水蒸氣反應(yīng)生成Ca(OH)2,再與CO2反應(yīng)生成CaCO3。由于碳酸化反應(yīng)路徑發(fā)生變化,在研究水蒸氣對(duì)碳酸化反應(yīng)時(shí)無法形成比較的基礎(chǔ)。在利用模型進(jìn)行計(jì)算時(shí),無法有效地對(duì)比碳酸化反應(yīng)速率,因此在研究水蒸氣對(duì)碳酸化反應(yīng)影響時(shí)需要剔除400℃時(shí)水蒸氣濃度10%和20%兩個(gè)實(shí)驗(yàn)點(diǎn)以及450℃時(shí)水蒸氣濃度20%的實(shí)驗(yàn)點(diǎn)。而對(duì)于含有1.5%的水蒸氣的反應(yīng)氣氛,在400~700℃范圍內(nèi)均不會(huì)形成Ca(OH)2。在選擇的溫度范圍內(nèi),20%的CO2壓力均高于碳酸化反應(yīng)平衡分壓。本研究的CO2及水蒸氣壓力范圍在圖1所示的兩條平衡分壓線之間,保證CaO在宏觀上僅與CO2反應(yīng)生成CaCO3。

    3.1 水蒸氣濃度對(duì)碳酸化反應(yīng)速率的影響

    圖2顯示了在500℃時(shí)水蒸氣濃度對(duì)碳酸化反應(yīng)影響。在反應(yīng)的初始階段,不同水蒸氣濃度條件下轉(zhuǎn)化率曲線相互重合,水蒸氣對(duì)碳酸化反應(yīng)的影響可以忽略。在反應(yīng)10 s之后,可以看到水蒸氣的存在明顯減緩了碳酸化反應(yīng)速率的下降。而且隨著水蒸氣濃度的增加轉(zhuǎn)化率有所增加,但增加幅度減弱。在反應(yīng)10 min之后,含有1.5%、10%和20%的水蒸氣條件下碳酸化轉(zhuǎn)化率比無水蒸氣條件下轉(zhuǎn)化率分別提高了19.8%、27.2%和30.5%。

    對(duì)于水蒸氣對(duì)碳酸化反應(yīng)的作用,Wang等[15]和Symonds等[16]認(rèn)為水蒸氣在碳酸化過程中起到催化作用,即水蒸氣會(huì)與CaO反應(yīng)瞬間形成Ca(OH)2,而Ca(OH)2比CaO具有更高的活性。也就是說,在被化學(xué)反應(yīng)動(dòng)力學(xué)所控制的快速反應(yīng)階段,反應(yīng)速率由于水蒸氣的存在而提高。然而這一推測(cè)受到一些學(xué)者的質(zhì)疑,Yang等[6]在研究低水蒸氣濃度對(duì)碳酸化反應(yīng)催化能力時(shí)觀察到,CaO與水蒸氣反應(yīng)生成Ca(OH)2,再與CO2反應(yīng)的速率明顯低于CaO和CO2反應(yīng)的速率。因此,由于Ca(OH)2的形成造成反應(yīng)速率提高的推測(cè)并不準(zhǔn)確。除此之外,在碳酸化反應(yīng)的初始階段,可能包含反應(yīng)氣體在樣品顆粒之間的外擴(kuò)散、化學(xué)反應(yīng)以及內(nèi)擴(kuò)散等步驟。表觀反應(yīng)速率由最慢的步驟決定,而化學(xué)反應(yīng)速率這一步驟比其他擴(kuò)散步驟快得多[17],化學(xué)反應(yīng)速率步驟反應(yīng)速率的提高并不一定促進(jìn)整體轉(zhuǎn)化率的提高。本實(shí)驗(yàn)在初始反應(yīng)階段的實(shí)驗(yàn)結(jié)果顯示,在快速反應(yīng)階段水蒸氣并沒有促進(jìn)碳酸化反應(yīng)發(fā)生。圖2顯示在反應(yīng)10 s之后,實(shí)驗(yàn)結(jié)果明顯體現(xiàn)出水蒸氣對(duì)碳酸化反應(yīng)的影響。在此反應(yīng)階段,可能存在反應(yīng)氣體在樣品顆??紫秲?nèi)部擴(kuò)散以及產(chǎn)物層擴(kuò)散。Bhatia等[18]在研究產(chǎn)物層對(duì)CaO與CO2的影響中發(fā)現(xiàn),對(duì)于粒徑為125~149mm的吸收劑顆粒具有與74~88mm幾乎重合的轉(zhuǎn)化率曲線,表明采用粒徑100~200mm的樣品孔隙內(nèi)擴(kuò)散的影響并不大,因此本文認(rèn)為水蒸氣對(duì)碳酸化反應(yīng)促進(jìn)作用主要發(fā)生在產(chǎn)物層擴(kuò)散階段。

    3.2 不同水蒸氣濃度產(chǎn)物層厚度的計(jì)算

    CaCO3產(chǎn)物層的厚度對(duì)碳酸化反應(yīng)是一個(gè)重要參數(shù),產(chǎn)物層厚度計(jì)算公式[19]如下

    圖3顯示了不同水蒸氣濃度條件下,在d/d下降到0.02時(shí)碳酸化反應(yīng)產(chǎn)物層厚度與反應(yīng)溫度關(guān)系。從圖中可以看出,在一定的水蒸氣濃度時(shí),產(chǎn)物層厚度隨溫度的提高而增加。在沒有水蒸氣存在條件下,400℃時(shí)產(chǎn)物層厚度只有6.4 nm,當(dāng)反應(yīng)溫度提高到700℃,產(chǎn)物層厚度增加到23.2 nm。在同一反應(yīng)溫度,提高水蒸氣濃度,產(chǎn)物層厚度相應(yīng)增加,在低溫時(shí)尤為明顯。400℃時(shí),即便存在微量的水蒸氣仍然使產(chǎn)物層厚度從6.4 nm增加到13.1 nm。但隨著水蒸氣濃度的不斷提高,產(chǎn)物層厚度的增加變得不明顯。隨著溫度的提高,不同水蒸氣濃度下產(chǎn)物層隨溫度變化曲線會(huì)相交于一點(diǎn),這表明即便沒有水蒸氣的存在,提高溫度也能起到類似的作用。關(guān)于反應(yīng)溫度對(duì)碳酸化產(chǎn)物層的影響將在后面進(jìn)行討論。

    3.3 不同水蒸氣濃度擴(kuò)散參數(shù)計(jì)算

    假定反應(yīng)速率較慢的第2個(gè)階段被產(chǎn)物層擴(kuò)散所控制,那么根據(jù)式(17),產(chǎn)物層擴(kuò)散速率可以由與時(shí)間的平方根的斜率計(jì)算得到,這一推測(cè)被圖4的數(shù)據(jù)關(guān)系所證實(shí)。圖4對(duì)比了在1.5%水蒸氣條件下,不同溫度所對(duì)應(yīng)的擴(kuò)散速率的大小,遵循線性關(guān)系的轉(zhuǎn)化率在45%~60%之間。表2顯示了不同水蒸氣濃度不同反應(yīng)溫度條件下,根據(jù)式(17)所計(jì)算的斜率。

    表2 碳酸化反應(yīng)有效擴(kuò)散系數(shù)及活化能Table 2 Effective product layer diffusivities and activation energy for carbonation

    對(duì)于某一反應(yīng)溫度,在反應(yīng)的初期和末期均存在非線性部分。這是由于在碳酸化反應(yīng)的初期存在碳酸化反應(yīng)的誘導(dǎo)期以及可能由氣體外擴(kuò)散和化學(xué)反應(yīng)速率步驟決定的初始反應(yīng)較快的反應(yīng)階段。而在碳酸化反應(yīng)的末期,由于吸收劑顆粒內(nèi)部孔隙的閉合導(dǎo)致擴(kuò)散面積的減小,因此最后階段斜率也產(chǎn)生了偏差。這一推測(cè)也可以從Manovic等[7]的實(shí)驗(yàn)結(jié)果得到證實(shí),他們比較了8種不同的石灰石樣品,發(fā)現(xiàn)無論石灰石在N2還是在CO2氣氛中煅燒,水蒸氣的存在均會(huì)提高碳酸化轉(zhuǎn)化率。然而石灰石在CO2中煅燒對(duì)碳酸化過程中水蒸氣的影響更為明顯,而在N2中煅燒雖然轉(zhuǎn)化率更高,但影響并不明顯。這是由于石灰石在CO2中煅燒時(shí)更容易發(fā)生燒結(jié)現(xiàn)象,產(chǎn)生更大的孔隙,使得在碳酸化過程中孔隙并不容易堵塞[20],在此條件下水蒸氣的作用顯得更為明顯。

    3.4 不同水蒸氣濃度所對(duì)應(yīng)反應(yīng)活化能的計(jì)算

    不同水蒸氣濃度條件下有效擴(kuò)散系數(shù)根據(jù)式(17)計(jì)算得到,產(chǎn)物層擴(kuò)散階段碳酸化反應(yīng)活化能也由此得到。不同水蒸氣濃度條件下,碳酸化反應(yīng)活化能見表2。圖5顯示出兩個(gè)關(guān)于擴(kuò)散系數(shù)的特征,首先大約在550℃處存在一個(gè)明顯的斜率變化,這一溫度與氣氛中是否存在水蒸氣無關(guān)。當(dāng)反應(yīng)溫度從400℃升高到700℃,在550℃處斜率增加,這表明在高于550℃后反應(yīng)機(jī)理發(fā)生了變化。高溫條件下,1.5%水蒸氣的存在使得無水蒸氣條件下碳酸化反應(yīng)活化能從237.7 kJ·mol-1降低到179.9 kJ·mol-1,繼續(xù)提高水蒸氣濃度到10%和20%,反應(yīng)活化能從156.6 kJ·mol-1降低到148.6 kJ·mol-1。

    對(duì)比活化能可以發(fā)現(xiàn),提高水蒸氣濃度降低了碳酸化反應(yīng)活化能,但繼續(xù)增加水蒸氣濃度,活化能的降低速度減弱。對(duì)比圖5的高溫區(qū)域擴(kuò)散系數(shù)的斜率,碳酸化反應(yīng)的另一個(gè)特征是隨著反應(yīng)溫度的提高,水蒸氣的促進(jìn)作用減弱,可以預(yù)計(jì)的是當(dāng)溫度足夠高,不同濃度水蒸氣條件下擴(kuò)散系數(shù)的直線必然相交一點(diǎn)。從圖5中觀察到的現(xiàn)象與臨界產(chǎn)物層厚度趨勢(shì)一致,不同水蒸氣濃度下臨界產(chǎn)物層隨溫度變化曲線會(huì)相交于一點(diǎn),反應(yīng)溫度可以起到與水蒸氣添加類似的作用。

    有證據(jù)表明產(chǎn)物層擴(kuò)散過程是一個(gè)固體狀態(tài)反應(yīng)過程:在產(chǎn)物層中較低的擴(kuò)散系數(shù)和較高的活化能并不可能是由氣體分子擴(kuò)散導(dǎo)致;活化能的改變發(fā)生在550℃處,與CaCO3的Tamman溫度(527℃)非常接近;實(shí)驗(yàn)結(jié)果與固體CaCO3中導(dǎo)電現(xiàn)象非常相似。Bhatia等[18]認(rèn)為CO32-及O2-構(gòu)成在CaCO3中相對(duì)運(yùn)動(dòng)的離子對(duì),而在反應(yīng)過程中保持產(chǎn)物層電中性。他們推測(cè)在CaO-CaCO3界面上發(fā)生如下反應(yīng)

    CO32-+CaOCaCO3+ O2-(20)

    而在孔隙表面為

    CO2(g)(CO2)ads(21)

    (CO2)ads+ O2-CO32-(22)

    在CO2分壓遠(yuǎn)高于平衡分壓條件下,根據(jù)這一機(jī)制,在反應(yīng)表面上吸收氣體的能力達(dá)到飽和,而在兩個(gè)反應(yīng)表面上達(dá)到了反應(yīng)平衡。隨著溫度的提高,O2-及CO32-的擴(kuò)散速率加快,導(dǎo)致了碳酸化反應(yīng)速率的提高,同時(shí)也促進(jìn)了產(chǎn)物層厚度的增加。

    根據(jù)不同水蒸氣濃度下活化能的計(jì)算,在產(chǎn)物層擴(kuò)散階段,水蒸氣的添加導(dǎo)致碳酸化反應(yīng)活化能下降。水蒸氣的質(zhì)量在反應(yīng)過程中并不減少,因此水蒸氣有可能起到催化作用。水蒸氣濃度增加,活化能減少,但減少的幅度逐漸下降,表明水蒸氣濃度的影響受到反應(yīng)條件的限制。因此本文提出,在水蒸氣存在條件下,除存在以上的離子擴(kuò)散之外,在孔隙表面還可能存在

    H2O+CO2=2H++ CO32-(23)

    2H++ O2-=H2O (24)

    在反應(yīng)表面發(fā)生

    CO32-+CaO= CaCO3+ O2-(25)

    在孔隙表面,水蒸氣和CO2發(fā)生弱電離,產(chǎn)生碳酸根離子。表面產(chǎn)物層的CaCO3吸附電離的碳酸根離子,原來形成CaCO3的碳酸根離子向內(nèi)部擴(kuò)散。在反應(yīng)表面碳酸根離子與CaO反應(yīng)產(chǎn)生氧離子,為保持電中性,氧離子向外擴(kuò)散直到孔隙表面,并與氫離子反應(yīng)重新生成水蒸氣,這樣在整個(gè)循環(huán)過程中,水蒸氣從宏觀角度并不發(fā)生變化,僅是促進(jìn)了CO2與CaO的反應(yīng)。在產(chǎn)物層擴(kuò)散階段碳酸化反應(yīng)在存在及不存在水蒸氣條件下的示意圖如圖6所示。

    Sun等[21]使用標(biāo)記實(shí)驗(yàn)證實(shí)了碳酸化反應(yīng)過程中離子外擴(kuò)散機(jī)理,即CaO中的O2-外擴(kuò)散,而CO32-內(nèi)擴(kuò)散。當(dāng)氣相中含有水蒸氣時(shí),除了這一擴(kuò)散之外,還有可能存在水蒸氣和CO2的電離過程。在孔隙表面電離出的碳酸根離子增加了物質(zhì)轉(zhuǎn)移速率,在反應(yīng)界面形成的O2-向外擴(kuò)散與H+重新生成H2O,使得水蒸氣本身數(shù)量并不發(fā)生變化,而是起到促進(jìn)作用。隨著水蒸氣濃度的增加,碳酸根離子增加,導(dǎo)致碳酸化反應(yīng)速率增加,然而繼續(xù)增加水蒸氣濃度則對(duì)反應(yīng)速率的增加作用減弱。當(dāng)反應(yīng)溫度提高后,水蒸氣與CO2的電離作用增強(qiáng),可以產(chǎn)生更多的CO32-參與反應(yīng),同樣也提高了碳酸化速率。

    對(duì)于在高于550℃后碳酸化反應(yīng)機(jī)理發(fā)生了變化,Bhatia等[18]根據(jù)反應(yīng)活化能接近CaCO3煅燒所需熱量推測(cè),在較高溫度下,CO2在產(chǎn)物層中以連續(xù)CaCO3分解的方式向內(nèi)部擴(kuò)散。即已形成的CaCO3中的CO32-在高溫條件下瞬間分解,形成CO2和O2-,CO2分子移動(dòng)到相鄰的空位上,而同樣生成的CO2分子占據(jù)原來的位置,重新生成CaCO3,這一CO2分子這樣通過不同位點(diǎn)在產(chǎn)物層內(nèi)擴(kuò)散,最后到達(dá)CaO-CaCO3反應(yīng)表面。本文認(rèn)為在較高溫度時(shí),即便CaCO3并不分解,可逆的CaCO3=CaO+CO2反應(yīng)也存在動(dòng)態(tài)平衡。碳酸化反應(yīng)根本仍然是CO32-在產(chǎn)物層內(nèi)的擴(kuò)散過程,由于高溫條件下部分CaCO3并不穩(wěn)定,因此加劇了CO32-在徑向產(chǎn)物層擴(kuò)散的難度,導(dǎo)致碳酸化反應(yīng)活化能提高。

    4 結(jié) 論

    天然石灰石經(jīng)過煅燒后,通過比較不同水蒸氣濃度對(duì)碳酸化反應(yīng)的影響,結(jié)合隨機(jī)孔隙模型的分析,得到如下結(jié)論。

    (1)碳酸化反應(yīng)溫度為500℃時(shí),10 min之后在水蒸氣含量在1.5%、10%和20%的條件下碳酸化轉(zhuǎn)化率比無水蒸氣條件下轉(zhuǎn)化率分別提高了19.8%、27.2%和30.5%。在初始的快速反應(yīng)階段,水蒸氣并不能明顯提高碳酸化反應(yīng)速率,明顯的促進(jìn)作用主要發(fā)生在產(chǎn)物層擴(kuò)散階段。

    (2)在相同的碳酸化反應(yīng)溫度條件下,提高水蒸氣濃度,碳酸化反應(yīng)產(chǎn)物層厚度相應(yīng)增加,在低溫時(shí)尤為明顯。400℃時(shí),即便存在微量的水蒸氣便使產(chǎn)物層厚度從6.4 nm增加到13.1 nm。隨著水蒸氣濃度的不斷提高,產(chǎn)物層厚度的增加速度變慢。不同水蒸氣濃度下產(chǎn)物層厚度隨溫度變化曲線會(huì)相交于一點(diǎn),表明隨著水蒸氣濃度的增加,對(duì)碳酸化反應(yīng)的促進(jìn)作用減弱。

    (3)在產(chǎn)物層擴(kuò)散階段,水蒸氣的存在導(dǎo)致碳酸化反應(yīng)活化能下降。但隨著水蒸氣濃度增加,下降的幅度減緩。當(dāng)反應(yīng)溫度從400℃升高到700℃過程中,活化能在550℃處發(fā)生突變,表明在高于550℃后反應(yīng)機(jī)理發(fā)生了變化。

    致謝:感謝中澳第二輪低排放煤炭技術(shù)研究基金的支持。

    符 號(hào) 說 明

    C——CO2濃度,kmol·m-3 Ce——CO2平衡濃度,kmol·m-3 Ci——在反應(yīng)界面CO2濃度,kmol·m-3 Cs——孔表面擴(kuò)散濃度,kmol·m-3 Dp——擴(kuò)散系數(shù),m2·s-1 Dpe——有效擴(kuò)散系數(shù),m2·s-1 kS——表面反應(yīng)速率系數(shù),m4·mol-1·s-1 L0——t=0時(shí)單位體積孔隙總長度,m M——CaO摩爾質(zhì)量,kg·mol-1 MCaCO3——CaCO3摩爾質(zhì)量,kg·mol-1 MCaO——CaO摩爾質(zhì)量,kg·mol-1 r——圓孔表面半徑,m S——單位體積反應(yīng)表面積,m2·m-3 S0——t=0時(shí)單位體積表面積,m2·m-3 V——單位空間反應(yīng)表面包括的孔體積,m3·m-3 V0——t=0時(shí)孔體積,m3·m-3 v0——孔體積分配函數(shù) X——碳酸化轉(zhuǎn)化率 Z——生成物與反應(yīng)物體積比 b——畢沃模量 D——產(chǎn)物層有效厚度,m e——孔隙率 e 0——t=0時(shí)孔隙率 y——結(jié)構(gòu)參數(shù) r——單位體積CaO質(zhì)量,kg·m-3 rb——固體堆密度,kg·m-3

    References

    [1] Donat F, Florin N H, Anthony E J, Fennell P S. Influence of high-temperature steam on the reactivity of CaO sorbent for CO2capture [J].,2012, 46(2): 1262-1269

    [2] Chen Huichao (陳惠超), Zhao Changsui (趙長遂), Shen Peng (沈鵬). Effect of steam in flue gas on CO2capture for calcium based sorbent [J].(化工學(xué)報(bào)),2013, 64(4): 1364-1372

    [3] Wang Chunbo (王春波), Zhou Xing (周興), Zheng Zhimin (鄭之民), Cui Caiyan (崔彩艷), Chen Liang (陳亮). Influence of steam on the limestone calcination/carbonation looping cycle for CO2capture [J].(中國電機(jī)工程學(xué)報(bào)),2014, 34(8): 1224-1230

    [4] Dobner S, Sterns L, Graff R A, Squires A M. Cyclic calcination and re-carbonation of calcined dolomite [J].,1977, 16(4): 479-486

    [5] Wang C B, Jia L F, Tan Y W, Anthony E J. Carbonation of fly ash in oxy-fuel CFB combustion [J].,2008, 87(7): 1108-1114

    [6] Yang S J, Xiao Y H. Steam catalysis in CaO carbonation under low steam partial pressure [J].,2008, 47(12): 4043-4048

    [7] Manovic V, Anthony E J. Carbonation of CaO-based sorbents enhanced by steam addition [J].,2010, 49(19): 9105-9110

    [8] Manovic V, Fennell P S, Al-Jeboori M J, Anthony E J. Steam-enhanced calcium looping cycles with calcium aluminate pellets doped with bromides [J].,2013, 52(23): 7677-7683

    [9] Homma S, Ogata S, Koga J, Matsumoto S. Gas-solid reaction model for a shrinking spherical particle with unreacted shrinking core [J].,2005, 60(18): 4971-4980

    [10] Khoshandam B, Kumar R V, Allahgholi L. Mathematical modeling of CO2removal using carbonation with CaO: the grain model [J].,2010, 27(3): 766-776

    [11] Bhatia S K, Perlmutter D D. A random pore model for fluid-solid reactions(Ⅰ): Isothermal, kinetic control [J].,1980, 26(3): 379-386

    [12] Bhatia S K, Perlmutter D D. A random pore model for fluid-solid reactions(Ⅱ): Diffusion and transport effects [J].,1981, 27(2): 247-254

    [13] Gavalas G R. A random capillary model with application to char gasification at chemically controlled rates [J].,1980, 26(4): 577-585

    [14] Fennell P S, Davidson J F, Dennis J S, Hayhurst A N. Regeneration of sintered limestone sorbents for the sequestration of CO2from combustion and other systems [J].,2007, 80(2): 116-119

    [15] Wang C B, Jia L F, Tan Y W, Anthony E J. The effect of water on the sulphation of limestone [J].,2010, 89(9): 2628-2632

    [16] Symonds R T, Lu D Y, Hughes R W, Anthony E J, Macchi A. CO2capture from simulated syngascyclic carbonation/calcination for a naturally occurring limestone: pilot-plant testing [J].,2009, 48(18): 8431-8440

    [17] Duo W, Sevill J P K, Kirkby N F, Clift R. Formation of product layers in solid-gas reactions for removal of acid gases [J].,1994, 49(24A): 4429-4442

    [18] Bhatia S K, Perlmutter D D. Effect of the product layer on the kinetics of the CO2-lime reaction [J].,1983, 29(1): 79-86

    [19] Barker R. Reversibility of reaction CaCO3reversible CaO+CO2[J].,1973, 23(10): 733-742

    [20] Stewart M C, Manovic V, Anthony E J, Macchi A. Enhancement of indirect sulphation of limestone by steam addition [J].,2010, 44(22): 8781-8786

    [21] Sun Z, Luo S, Qi P, Fan L S. Ionic diffusion through calcite (CaCO3) layer during the reaction of CaO and CO2[J].,2012, 81(22): 164-168

    Experimental and model study of carbonation of calcined limestone in presence of steam

    LIU Yang, YANG Yongping

    School of EnergyPower and Mechanical EngineeringNorth China Electric Power UniversityBeijingChina

    Steam is present in combustion flue gas, oxy-fuel combustor/calciner and fuel gas. Some previous works in this field have examined the reactivity of calcium oxide in the presence of steam. There is general agreement that the presence of steam increases the rate of carbonation even at low concentrations. However, there is no description about the effect of steam on the carbonation reaction of calcined limestone. The effects of temperature, concentration of steam on carbonation were investigated in a thermogravimetric analyzer. The understanding of the mechanisms participating in the carbonation reaction could be arrived at through the use of model interpretations of the rate controlling process. The experimental data were analyzed by means of the random pore model. During carbonation reaction, a very initial rapid reaction was followed by the second stage of the reaction occurring in the next slower regime. The experiment results showed that the effect of steam on carbonation could be neglected during the first stage of reaction. However, the conversion of carbonation in the presence of 1.5%(vol), 10%(vol) and 20%(vol) steam were 19.8%, 27.2% and 30.5% higher than the conversion without steam after 10 min during the second stage of reaction at 500℃. Steam was beneficial to enhancing the conversion of carbonation, but the extent of increase became not significant with increasing concentration of steam. The linearity of the second stage reaction data strongly suggested that this stage was controlled by a diffusion process occurring in a layer of calcium carbonate surrounding the calcium oxide in the pores of the solid. The parameters related to effective product layer diffusivities were plotted in Arrhenius form and the changes in activation energy at various steam concentrations were also shown. Comparison with the experimental effective diffusivities calculated for the carbonation reaction showed that the activation energies in the lower temperature range were in agreement with those obtained in the conductivity measurements. This suggested that the product layer diffusion process was proceeded by a mechanism similar to that of conduction in calcium carbonate. The activation energy of carbonation decreased from 179.9 kJ·mol-1to 237.7 kJ·mol-1when 1.5% steam was added according to the random pore model. The activation energy was 156.6 kJ·mol-1and 148.6 kJ·mol-1respectively for the atmospheres of 10% and 20% steam. There were two characteristics of product layer diffusion during carbonation. One was the slope of diffusional coefficient to increase from about 550℃, which was irrelevant to the presence of steam. The other was the effect of steam on carbonation to turn weak with increasing concentration of steam. At the higher temperature range, the possible mechanism for effective diffusion could be associated with sequential decompositions of carbonate ions in the calcium carbonate layer. A carbonate ion momentarily decomposed to generate carbon dioxide and an oxygen ion. The carbon dioxide molecule then moved to a neighboring similarly vacated site, while another carbon dioxide so generated elsewhere moved to take its place and reform the carbonate ion. In this way by a site to site random walk the carbon dioxide molecule diffused through the product layer, before reaction at the interface of calcium oxide and calcium carbonate. Such a mechanism appeared more prominent than the motion of carbonate ion at the higher temperature range. Both the above mechanisms for the diffusion stage were likely, and the true situation might involve one or a combination of both.

    CO2capture;chemical reaction;catalysis;product layer diffusion coefficient; activation energy

    2014- 09- 22.

    Prof. YANG Yongping, yypncepu@163.com

    10.11949/j.issn.0438-1157.20141419

    TK 09

    A

    0438—1157(2015)03—1088—09

    國家杰出青年科學(xué)基金項(xiàng)目(51025624);111 引智項(xiàng)目(B12034);中央高校基金項(xiàng)目(2014QN08)。

    2014-09-22收到初稿,2014-11-24收到修改稿。

    聯(lián)系人:楊勇平。第一作者:劉洋(1977—),男,博士后。

    supported by the National Natural Science Foundation of China(51025624), the 111 Project (B12034) and the Fundamental Research Funds for the Central Universities (2014QN08).

    猜你喜歡
    碳酸水蒸氣轉(zhuǎn)化率
    什么!碳酸飲料要斷供了?
    藏在空氣里的水蒸氣
    我國全產(chǎn)業(yè)領(lǐng)域平均國際標(biāo)準(zhǔn)轉(zhuǎn)化率已達(dá)75%
    冒泡的可樂
    水蒸氣變戲法
    “碳酸鈉與碳酸氫鈉”知識(shí)梳理
    曲料配比與米渣生醬油蛋白質(zhì)轉(zhuǎn)化率的相關(guān)性
    透視化學(xué)平衡中的轉(zhuǎn)化率
    迷宮篇
    鑭石型碳酸鐠釹向堿式碳酸鐠釹的相轉(zhuǎn)變反應(yīng)特征及其應(yīng)用
    午夜老司机福利片| 国产精品久久视频播放| 成人国产综合亚洲| 可以免费在线观看a视频的电影网站| 狠狠狠狠99中文字幕| 身体一侧抽搐| АⅤ资源中文在线天堂| 久久热在线av| 久久九九热精品免费| 国产午夜精品久久久久久| 欧美成人午夜精品| 我的亚洲天堂| 亚洲人成伊人成综合网2020| 成人手机av| 制服丝袜大香蕉在线| 精品人妻在线不人妻| 搞女人的毛片| 久久久久精品国产欧美久久久| 可以在线观看的亚洲视频| 欧美激情极品国产一区二区三区| av视频免费观看在线观看| 亚洲熟女毛片儿| 久久国产亚洲av麻豆专区| 成年女人毛片免费观看观看9| 亚洲精华国产精华精| 在线天堂中文资源库| 国产野战对白在线观看| 亚洲欧洲精品一区二区精品久久久| 可以在线观看的亚洲视频| 国产欧美日韩一区二区精品| 电影成人av| 国产精品久久视频播放| 久久中文字幕一级| 亚洲熟妇中文字幕五十中出| 亚洲第一青青草原| 丁香欧美五月| 午夜福利18| 亚洲少妇的诱惑av| 在线观看免费视频日本深夜| 九色亚洲精品在线播放| 天天一区二区日本电影三级 | 国产精品野战在线观看| 久久精品国产99精品国产亚洲性色 | 淫妇啪啪啪对白视频| 国产区一区二久久| 啪啪无遮挡十八禁网站| 午夜福利免费观看在线| 国产欧美日韩精品亚洲av| 最好的美女福利视频网| 黄色 视频免费看| 亚洲五月色婷婷综合| 久久久久国产精品人妻aⅴ院| av电影中文网址| 精品高清国产在线一区| 午夜免费鲁丝| 国产成人av激情在线播放| 久久久久久久久免费视频了| 老司机福利观看| 夜夜夜夜夜久久久久| av福利片在线| 国产亚洲av嫩草精品影院| 这个男人来自地球电影免费观看| 国产三级黄色录像| 免费在线观看完整版高清| 国产99白浆流出| 国产精品一区二区精品视频观看| 欧美日韩精品网址| 777久久人妻少妇嫩草av网站| 欧美绝顶高潮抽搐喷水| 中文字幕精品免费在线观看视频| 亚洲中文日韩欧美视频| 午夜日韩欧美国产| 窝窝影院91人妻| 最近最新中文字幕大全电影3 | 极品教师在线免费播放| 亚洲精品美女久久久久99蜜臀| avwww免费| 两个人免费观看高清视频| 美女免费视频网站| 欧美日本亚洲视频在线播放| 大型黄色视频在线免费观看| 久久午夜亚洲精品久久| 99国产综合亚洲精品| 成在线人永久免费视频| 免费看美女性在线毛片视频| 国产成人精品无人区| 久久天躁狠狠躁夜夜2o2o| 18禁裸乳无遮挡免费网站照片 | 亚洲精品在线观看二区| 欧美亚洲日本最大视频资源| 国产一级毛片七仙女欲春2 | 亚洲色图av天堂| 无限看片的www在线观看| 久久久久国产精品人妻aⅴ院| 无限看片的www在线观看| 欧美不卡视频在线免费观看 | 国产三级黄色录像| 成人亚洲精品一区在线观看| 十八禁网站免费在线| 一进一出好大好爽视频| 久久中文看片网| 丰满的人妻完整版| 黄色女人牲交| 亚洲熟妇熟女久久| 亚洲精品美女久久久久99蜜臀| av有码第一页| 久久婷婷人人爽人人干人人爱 | 久久久精品国产亚洲av高清涩受| 97超级碰碰碰精品色视频在线观看| 国产欧美日韩综合在线一区二区| 精品熟女少妇八av免费久了| 日本一区二区免费在线视频| 欧美色欧美亚洲另类二区 | 妹子高潮喷水视频| 亚洲精品国产区一区二| 嫩草影视91久久| 色在线成人网| 午夜a级毛片| 欧美黄色片欧美黄色片| 免费高清视频大片| 久久精品人人爽人人爽视色| 成人亚洲精品一区在线观看| 国产亚洲精品久久久久5区| 视频在线观看一区二区三区| 国产色视频综合| 免费不卡黄色视频| 成人三级做爰电影| 欧美大码av| 级片在线观看| 精品高清国产在线一区| 9热在线视频观看99| 老司机靠b影院| 波多野结衣av一区二区av| 国产一区在线观看成人免费| 久久天堂一区二区三区四区| 国产成年人精品一区二区| 精品一品国产午夜福利视频| 99国产精品免费福利视频| 国产一区二区在线av高清观看| 亚洲狠狠婷婷综合久久图片| 两个人视频免费观看高清| 91国产中文字幕| 久久草成人影院| 国产激情欧美一区二区| 久久久久久久久免费视频了| 超碰成人久久| 一进一出抽搐gif免费好疼| 国产又色又爽无遮挡免费看| 岛国视频午夜一区免费看| 91成年电影在线观看| 国产精品久久视频播放| 可以在线观看的亚洲视频| 欧美一级毛片孕妇| 啦啦啦观看免费观看视频高清 | 宅男免费午夜| 长腿黑丝高跟| av在线天堂中文字幕| 国产成人精品无人区| aaaaa片日本免费| 国产在线精品亚洲第一网站| 欧美国产精品va在线观看不卡| 国产精品1区2区在线观看.| 免费在线观看影片大全网站| 国产三级在线视频| 久久婷婷人人爽人人干人人爱 | 亚洲第一电影网av| 国产日韩一区二区三区精品不卡| 精品国产乱子伦一区二区三区| 精品高清国产在线一区| 超碰成人久久| 一区二区三区高清视频在线| 精品久久久久久,| 久热这里只有精品99| 久久精品人人爽人人爽视色| 满18在线观看网站| 亚洲一卡2卡3卡4卡5卡精品中文| 非洲黑人性xxxx精品又粗又长| 丰满人妻熟妇乱又伦精品不卡| 男男h啪啪无遮挡| 在线十欧美十亚洲十日本专区| 非洲黑人性xxxx精品又粗又长| 国产欧美日韩综合在线一区二区| 成人手机av| 精品日产1卡2卡| 国产成人免费无遮挡视频| 午夜激情av网站| 午夜免费成人在线视频| 少妇被粗大的猛进出69影院| 国产一区二区三区在线臀色熟女| 亚洲av第一区精品v没综合| 黄片大片在线免费观看| www国产在线视频色| 人人澡人人妻人| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三区视频在线观看免费| 午夜福利18| 久久香蕉国产精品| 69av精品久久久久久| av福利片在线| 日本三级黄在线观看| 欧美绝顶高潮抽搐喷水| 亚洲精品国产区一区二| 美女免费视频网站| 天堂动漫精品| 十八禁人妻一区二区| 久久久久久亚洲精品国产蜜桃av| 黄色毛片三级朝国网站| 久久香蕉精品热| 国产av在哪里看| 欧美日本中文国产一区发布| 丝袜人妻中文字幕| 亚洲人成伊人成综合网2020| 久久久久久人人人人人| 女人精品久久久久毛片| 国产成人av教育| 亚洲av电影不卡..在线观看| 精品久久久久久久久久免费视频| 久久精品亚洲熟妇少妇任你| 中文亚洲av片在线观看爽| 巨乳人妻的诱惑在线观看| 黄色视频,在线免费观看| 国产熟女午夜一区二区三区| www.熟女人妻精品国产| 午夜福利影视在线免费观看| 丝袜在线中文字幕| 国产精品久久久久久精品电影 | 亚洲av成人av| 亚洲精品久久国产高清桃花| 久久久久国产一级毛片高清牌| 亚洲中文日韩欧美视频| 一区二区三区高清视频在线| 日韩 欧美 亚洲 中文字幕| 色综合亚洲欧美另类图片| av视频在线观看入口| 老司机午夜福利在线观看视频| 精品久久久久久久人妻蜜臀av | 日韩欧美一区视频在线观看| 欧美黑人欧美精品刺激| 午夜福利在线观看吧| 亚洲第一欧美日韩一区二区三区| 中文字幕久久专区| 亚洲五月色婷婷综合| av天堂在线播放| 露出奶头的视频| 桃色一区二区三区在线观看| 国产伦一二天堂av在线观看| 欧美在线黄色| 一夜夜www| 亚洲成国产人片在线观看| 亚洲一区二区三区不卡视频| 天天躁夜夜躁狠狠躁躁| 1024香蕉在线观看| 这个男人来自地球电影免费观看| a在线观看视频网站| 啦啦啦观看免费观看视频高清 | 国产精品一区二区三区四区久久 | 十八禁网站免费在线| 免费搜索国产男女视频| 色av中文字幕| 黑人操中国人逼视频| 大型黄色视频在线免费观看| 国产精品香港三级国产av潘金莲| 久久久精品欧美日韩精品| 非洲黑人性xxxx精品又粗又长| 精品国产乱子伦一区二区三区| 色哟哟哟哟哟哟| 免费高清视频大片| 午夜免费激情av| 9色porny在线观看| 身体一侧抽搐| 天堂√8在线中文| 亚洲成av人片免费观看| 久久影院123| 亚洲人成网站在线播放欧美日韩| 日韩欧美一区视频在线观看| 日本 欧美在线| 亚洲男人的天堂狠狠| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品在线美女| 99精品欧美一区二区三区四区| 在线观看免费日韩欧美大片| 亚洲欧洲精品一区二区精品久久久| 国产免费男女视频| 亚洲第一av免费看| 国产成人精品无人区| 日本三级黄在线观看| 777久久人妻少妇嫩草av网站| 男人舔女人下体高潮全视频| av有码第一页| 99精品久久久久人妻精品| 看免费av毛片| 亚洲欧美激情在线| 91麻豆精品激情在线观看国产| 欧美在线一区亚洲| 午夜久久久在线观看| 免费无遮挡裸体视频| 亚洲国产精品sss在线观看| 亚洲九九香蕉| 女人精品久久久久毛片| 国产熟女xx| 午夜影院日韩av| 免费看美女性在线毛片视频| 亚洲欧美激情综合另类| 中文字幕高清在线视频| ponron亚洲| 欧美黄色片欧美黄色片| 久久久国产成人免费| 免费观看精品视频网站| 久久久久久人人人人人| 欧美激情极品国产一区二区三区| 欧美日韩一级在线毛片| 国产99久久九九免费精品| 一区二区三区高清视频在线| 天天添夜夜摸| 少妇的丰满在线观看| 激情在线观看视频在线高清| 日韩欧美一区视频在线观看| 在线观看午夜福利视频| 久久草成人影院| www日本在线高清视频| 涩涩av久久男人的天堂| 黑丝袜美女国产一区| 免费不卡黄色视频| 波多野结衣巨乳人妻| 免费看a级黄色片| 亚洲精品中文字幕一二三四区| 亚洲av电影不卡..在线观看| 女同久久另类99精品国产91| 国产成年人精品一区二区| 欧美色视频一区免费| 久久久久久久久免费视频了| 亚洲 欧美 日韩 在线 免费| 久久久水蜜桃国产精品网| 操出白浆在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 不卡av一区二区三区| 国产在线精品亚洲第一网站| 又黄又爽又免费观看的视频| 国产高清视频在线播放一区| 日韩三级视频一区二区三区| 欧美成狂野欧美在线观看| 中文字幕色久视频| www.熟女人妻精品国产| 国产激情久久老熟女| 亚洲av美国av| 99久久久亚洲精品蜜臀av| 久久青草综合色| 中文字幕另类日韩欧美亚洲嫩草| 别揉我奶头~嗯~啊~动态视频| 精品欧美一区二区三区在线| 99国产精品99久久久久| 欧美日韩瑟瑟在线播放| 97人妻精品一区二区三区麻豆 | 老熟妇仑乱视频hdxx| 午夜免费观看网址| 久久香蕉激情| 精品卡一卡二卡四卡免费| 成人特级黄色片久久久久久久| 亚洲熟妇中文字幕五十中出| 免费在线观看黄色视频的| 亚洲av成人一区二区三| 久久久久亚洲av毛片大全| 亚洲精品粉嫩美女一区| 午夜老司机福利片| 无限看片的www在线观看| 免费在线观看影片大全网站| 国产精品二区激情视频| 久久九九热精品免费| 一区二区三区国产精品乱码| 亚洲伊人色综图| 最近最新中文字幕大全免费视频| 成人18禁高潮啪啪吃奶动态图| 老汉色∧v一级毛片| 99精品欧美一区二区三区四区| 亚洲五月婷婷丁香| 国产1区2区3区精品| 精品久久久久久,| 又大又爽又粗| 99国产精品99久久久久| 日韩成人在线观看一区二区三区| 久久精品国产亚洲av高清一级| 级片在线观看| 99久久久亚洲精品蜜臀av| 中文字幕另类日韩欧美亚洲嫩草| 中出人妻视频一区二区| 一边摸一边做爽爽视频免费| 国产免费av片在线观看野外av| 欧美日韩一级在线毛片| 在线观看舔阴道视频| 男人舔女人下体高潮全视频| 高清毛片免费观看视频网站| 免费在线观看日本一区| 久久久久久人人人人人| 中文字幕另类日韩欧美亚洲嫩草| 老司机福利观看| 亚洲精品中文字幕在线视频| 天堂动漫精品| 很黄的视频免费| 精品福利观看| 天堂动漫精品| 久热爱精品视频在线9| 99精品欧美一区二区三区四区| bbb黄色大片| 午夜激情av网站| 99在线人妻在线中文字幕| 久久精品国产99精品国产亚洲性色 | 国产私拍福利视频在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲中文日韩欧美视频| 欧美一区二区精品小视频在线| 免费av毛片视频| 午夜精品久久久久久毛片777| 午夜成年电影在线免费观看| 亚洲自偷自拍图片 自拍| 男人舔女人下体高潮全视频| 欧美人与性动交α欧美精品济南到| 国产精品亚洲美女久久久| 国产一卡二卡三卡精品| 99久久国产精品久久久| 激情在线观看视频在线高清| 久久精品影院6| 老司机午夜福利在线观看视频| 精品国产超薄肉色丝袜足j| 国内毛片毛片毛片毛片毛片| 亚洲av电影不卡..在线观看| 久久婷婷成人综合色麻豆| 精品国产乱码久久久久久男人| 欧洲精品卡2卡3卡4卡5卡区| ponron亚洲| 国产精品久久久久久人妻精品电影| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美日韩无卡精品| 亚洲五月婷婷丁香| 老司机深夜福利视频在线观看| 99国产精品免费福利视频| 国产亚洲欧美98| 999久久久精品免费观看国产| av在线天堂中文字幕| 久久影院123| 欧美日本亚洲视频在线播放| 免费看a级黄色片| 欧美绝顶高潮抽搐喷水| 美女高潮喷水抽搐中文字幕| 十分钟在线观看高清视频www| 女同久久另类99精品国产91| 国产精品av久久久久免费| 多毛熟女@视频| 成人三级黄色视频| 黄片大片在线免费观看| 欧美乱码精品一区二区三区| 在线免费观看的www视频| 成人三级做爰电影| 在线观看www视频免费| 免费看美女性在线毛片视频| 我的亚洲天堂| 精品国产一区二区三区四区第35| 久久香蕉激情| 一级黄色大片毛片| 男女之事视频高清在线观看| 日日爽夜夜爽网站| 在线永久观看黄色视频| 97超级碰碰碰精品色视频在线观看| 精品不卡国产一区二区三区| 欧美日韩亚洲综合一区二区三区_| 久久影院123| 成人三级黄色视频| 三级毛片av免费| 黑人巨大精品欧美一区二区mp4| 免费观看人在逋| 欧美久久黑人一区二区| av在线播放免费不卡| 亚洲色图 男人天堂 中文字幕| 激情在线观看视频在线高清| 久久久久久久午夜电影| 亚洲精品国产一区二区精华液| 丝袜在线中文字幕| 色精品久久人妻99蜜桃| 欧美国产精品va在线观看不卡| 国产精品九九99| 国产成人av激情在线播放| 中文字幕人成人乱码亚洲影| 99久久国产精品久久久| 一二三四在线观看免费中文在| 日本a在线网址| 热re99久久国产66热| 亚洲成av片中文字幕在线观看| 人妻久久中文字幕网| 天天躁夜夜躁狠狠躁躁| 久久久久九九精品影院| 精品高清国产在线一区| 亚洲自拍偷在线| av网站免费在线观看视频| 色尼玛亚洲综合影院| 国产精品亚洲一级av第二区| 亚洲熟妇中文字幕五十中出| 久久久国产成人免费| 丝袜美足系列| 无遮挡黄片免费观看| 满18在线观看网站| 欧美一级毛片孕妇| 欧美不卡视频在线免费观看 | or卡值多少钱| 午夜福利视频1000在线观看 | 日韩成人在线观看一区二区三区| 免费搜索国产男女视频| 制服人妻中文乱码| 一区二区三区高清视频在线| 亚洲性夜色夜夜综合| 老司机午夜福利在线观看视频| 国产乱人伦免费视频| 国产成年人精品一区二区| 黄色丝袜av网址大全| 男女做爰动态图高潮gif福利片 | 免费在线观看视频国产中文字幕亚洲| av有码第一页| 怎么达到女性高潮| 99热只有精品国产| 欧美一级毛片孕妇| 成人三级黄色视频| 操美女的视频在线观看| 亚洲一区二区三区色噜噜| 日日干狠狠操夜夜爽| 色综合婷婷激情| 亚洲精品国产精品久久久不卡| 精品久久久久久成人av| 成人三级黄色视频| 桃色一区二区三区在线观看| 国产精品综合久久久久久久免费 | 高清在线国产一区| 午夜a级毛片| 免费人成视频x8x8入口观看| 国产免费男女视频| 99国产综合亚洲精品| 国产伦一二天堂av在线观看| av在线天堂中文字幕| 热re99久久国产66热| 琪琪午夜伦伦电影理论片6080| 亚洲中文日韩欧美视频| 女生性感内裤真人,穿戴方法视频| 高潮久久久久久久久久久不卡| 在线免费观看的www视频| 国产精品一区二区在线不卡| 国产一卡二卡三卡精品| 国产av一区在线观看免费| avwww免费| 多毛熟女@视频| 精品国产美女av久久久久小说| 国产成人啪精品午夜网站| 精品国产一区二区三区四区第35| 欧美成人免费av一区二区三区| 亚洲第一av免费看| 精品欧美一区二区三区在线| 国产成人系列免费观看| 久久精品国产99精品国产亚洲性色 | 最近最新免费中文字幕在线| 一边摸一边抽搐一进一出视频| 日本a在线网址| 亚洲一区二区三区不卡视频| 老汉色av国产亚洲站长工具| 日日摸夜夜添夜夜添小说| av有码第一页| 国产亚洲精品综合一区在线观看 | 亚洲国产中文字幕在线视频| 日韩 欧美 亚洲 中文字幕| 人妻久久中文字幕网| 男人舔女人下体高潮全视频| 精品久久久久久成人av| av在线播放免费不卡| 多毛熟女@视频| 一级a爱片免费观看的视频| 成人三级做爰电影| 国产精品香港三级国产av潘金莲| 自拍欧美九色日韩亚洲蝌蚪91| 欧美在线黄色| 欧洲精品卡2卡3卡4卡5卡区| 欧美午夜高清在线| 亚洲在线自拍视频| 国产一区二区在线av高清观看| 成人欧美大片| 色综合站精品国产| 99久久99久久久精品蜜桃| 日韩欧美免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 免费少妇av软件| 成人国语在线视频| 97人妻精品一区二区三区麻豆 | 在线观看午夜福利视频| 亚洲国产高清在线一区二区三 | 90打野战视频偷拍视频| 好看av亚洲va欧美ⅴa在| 狠狠狠狠99中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品久久国产高清桃花| 激情视频va一区二区三区| 一区在线观看完整版| 精品电影一区二区在线| 免费不卡黄色视频| 美女免费视频网站| 级片在线观看| 99在线人妻在线中文字幕| 国产亚洲精品av在线| 免费女性裸体啪啪无遮挡网站| 欧美日韩瑟瑟在线播放| 又大又爽又粗| 国产高清视频在线播放一区| 99久久精品国产亚洲精品| 国产成人系列免费观看| 国产精品永久免费网站| 日韩欧美三级三区| 国产成人免费无遮挡视频| 国产精品av久久久久免费| 亚洲午夜精品一区,二区,三区| 性色av乱码一区二区三区2| 国产精品国产高清国产av| 18禁裸乳无遮挡免费网站照片 |