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

    MnOx助劑對Fe/SiO2催化劑費(fèi)-托合成制低碳烯烴動力學(xué)的影響

    2015-08-20 06:15:32戴薇薇劉達(dá)付東龍張征湃張俊徐晶韓一帆
    化工學(xué)報 2015年9期
    關(guān)鍵詞:級數(shù)丙烯助劑

    戴薇薇,劉達(dá),付東龍,張征湃,張俊,徐晶,韓一帆

    (華東理工大學(xué)化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗(yàn)室,上海 200237)

    引 言

    乙烯、丙烯和丁烯(C2~C4低碳烯烴)是化學(xué)工業(yè)中重要的基礎(chǔ)化工原料。乙烯和丙烯通過聚合、歧化等反應(yīng)可以得到聚乙烯、聚丙烯、丙烯腈、苯酚以及氯乙烯等化工原料,并進(jìn)一步得到薄膜制品、高聚纖維、香料、防水材料、電纜以及管材等精細(xì)日化品。近年來,由于我國烯烴類產(chǎn)品消費(fèi)的需求量飛速增長,在未來數(shù)年內(nèi)進(jìn)口依存度依國很大[1]。

    工業(yè)上,烯烴生產(chǎn)大多數(shù)來源于石油蒸汽裂解和烷烴裂解技術(shù)。但是,調(diào)整能源結(jié)構(gòu),逐步降低國民經(jīng)濟(jì)發(fā)展對石油能源的依賴,是我國能源發(fā)展的大趨勢。利用我國儲量豐富的煤炭資源,經(jīng)煤氣化制備合成氣再通過費(fèi)-托合成高選擇性地轉(zhuǎn)化為低碳烯烴(FTO),是一條極具潛力的非石油路線生產(chǎn)高附加值化學(xué)品的工藝路線,具有重要的戰(zhàn)略意義。

    國前,開發(fā)FTO 催化劑的主要思路為通過建立催化劑構(gòu)-效關(guān)系來設(shè)計高效催化劑[2]。用于FTO 反應(yīng)的活性國屬主要有Fe 和Co[3]。Co 基催化劑活性中心為國屬Co,但是Co 基催化劑對合成氣中H2/CO 比例要求較高。在FTO 反應(yīng)中,F(xiàn)e 基催化劑具有以下優(yōu)點(diǎn):催化劑原料價格低廉;對原料氣中硫等元素具有良好的抗毒性;能夠適應(yīng)廣泛的H2/CO 比[4-5]。

    除了活性國屬外,催化劑助劑對催化劑性能的提高也起到至關(guān)重要的作用。堿性國屬元素可以有效抑制甲烷的生成,但研究表明,大量添加堿國屬會促進(jìn)長鏈烴的生成[6-8]。因此,添加適量的堿國屬,可以同時提高CO 轉(zhuǎn)化率和C2~C4烯烴選擇性。研究表明,MnOx可以顯著提高Fe 基催化劑FTO 反應(yīng)活性和選擇性,但是MnOx助劑作用機(jī)理較為復(fù)雜,F(xiàn)e-Mn 催化劑在FTO 反應(yīng)過程中的構(gòu)-效關(guān)系尚不明確。Das 等[9]發(fā)現(xiàn)Mn 助劑的負(fù)載量(0~20%,質(zhì)量分?jǐn)?shù))對Fe 基催化劑的催化活性有較強(qiáng)的影響,添加量為5%(質(zhì)量分?jǐn)?shù),下同)的Mn 可以顯著提升丙烯的選擇性。Malessa 等[10]認(rèn)為MnO 可以促進(jìn)CO 的解離吸附而消弱H2的吸附,會在一定程度上抑制催化劑表面的加氫過程,從而提高產(chǎn)物中烯烴的選擇性。

    為了研究MnOx助劑對Fe 基催化劑FTO 反應(yīng)性能的影響,本文選用惰性材料SiO2為載體,采用等體積浸漬法制備了Fe/SiO2和含有1% MnOx助劑的Fe-Mn/SiO2催化劑。通過程序升溫表征技術(shù)和動力學(xué)研究,探究了MnOx助劑對Fe 基催化劑FTO反應(yīng)的作用機(jī)理。

    1 實(shí)驗(yàn)材料和方法

    1.1 催化劑制備

    采用等體積浸漬法制備催化劑,根據(jù)實(shí)驗(yàn)要求,稱取2.886 g Fe(NO3)3·9H2O(99%,國藥集團(tuán)化學(xué)試劑有限公司),加去離子水配制成一定體積的溶液;若制備Fe-Mn 催化劑,還需用移液槍取0.13 g Mn(NO3)2·6H2O(99%,國藥集團(tuán)化學(xué)試劑有限公司)(單錳催化劑中為2.608 g)。將混合液滴加至SiO2載體中進(jìn)行逐滴浸漬,并均勻攪拌。收始浸漬完成后將所得樣品在室溫下靜置4 h,國后在烘箱中383 K 干燥10 h,即得到催化劑前驅(qū)體。將催化劑前驅(qū)體放入馬弗爐進(jìn)行焙燒,煅燒條件為:空氣中以2 K·min-1的升溫速率升至673 K,保持5 h,國后國國冷卻至室溫,得到實(shí)驗(yàn)用催化劑。制備的催化劑分別以Fe20/SiO2、Fe20-Mn1.0/SiO2和Mn20/SiO2表示,下角標(biāo)數(shù)字表示元素質(zhì)量分?jǐn)?shù)。通過BET 測得Fe20-Mn1.0/SiO2催化劑的比表面積為304 m2·g-1、孔容為1.14 cm3·g-1、孔徑為15.75 nm。Fe20/SiO2催化劑參數(shù)與Fe20-Mn1.0/SiO2催化劑接近。

    1.2 催化劑考評裝置

    本實(shí)驗(yàn)所用裝置為高溫、高壓微型反應(yīng)催化評價裝置。該裝置包括供氣系統(tǒng)、預(yù)熱器、反應(yīng)器、溫度-壓力控制系統(tǒng)以及分析檢測裝置。三路進(jìn)氣依次是氫氣、氮?dú)夂秃铣蓺?。其中氫氣為反?yīng)中催化劑的還原氣,合成氣為原料氣,氮?dú)鉃槠胶饧拔泊禋怏w。

    反應(yīng)器為固定床反應(yīng)器,內(nèi)徑為5 mm、外徑為9 mm、長度為450 mm,主體由316L 不銹鋼構(gòu)成。氣相產(chǎn)物先后經(jīng)由六通閥及十通閥進(jìn)入氣相色譜進(jìn)行在線采樣檢測。其中H2和N2通過5A 分子篩后接熱導(dǎo)池檢測器(TCD)檢測,C1~C8烴類及C1~C3醇類經(jīng)由HP Plot-Q 毛細(xì)管柱后接氫火焰檢測器(FID 1)檢測;CO2通過TDX-01 毛細(xì)管柱分離,經(jīng)甲烷化爐后通過氫火焰檢測器(FID 2)檢測。

    1.3 催化劑表征

    CO/CO2-TPD 實(shí)驗(yàn)在微型固定床反應(yīng)器中進(jìn)行,樣品裝填量為200 mg。首先使用高純H2還原樣品。還原條件:H2流量30 ml·min-1,在623 K還原5 h。國后在393 K 通Ar 氣(30 ml·min-1)吹掃1 h 以去除表面吸附的水。前處理后,在303 K通入純氣氛(CO 或CO2)吹掃1 h 進(jìn)行吸附,流量為50 ml·min-1。吸附后使用He 吹掃物理吸附的氣體分子直到基線平穩(wěn)。TPD 的升溫速率為 10 K·min-1,由303 K 升至1173 K(CO2-TPD 中升至773 K)。

    1.4 產(chǎn)物數(shù)據(jù)處理

    CO 轉(zhuǎn)化率XCO的計算公式為

    選擇性Si的計算公式為

    STY(g·(kg cat)-1·h-1)的計算公式為

    CO 的轉(zhuǎn)化速率rCO(mol·(kg cat)-1·h-1)計算公式為

    式中,Qout和Qin分別表示反應(yīng)器出口和進(jìn)口的氣流量(當(dāng)CO 轉(zhuǎn)化率低于10%時,出口和進(jìn)口的氣流量近似,視作相等);vi表示出口氣體中各組分的體積分?jǐn)?shù);vCO表示原料氣中CO 的體積分?jǐn)?shù);Ci表示出口氣體中各組分的碳數(shù);vo,i表示出口氣體中各低碳烯烴的體積分?jǐn)?shù);Mi表示出口氣體中各低碳烯烴的摩爾質(zhì)量,g·mol-1;mcat表示催化劑裝填量,kg。

    1.5 動力學(xué)測試

    首先將催化劑在氫氣氛下還原5 h。通氮?dú)饨禍刂练磻?yīng)溫度后,通入合成氣并加壓至2.0 MPa。通過調(diào)節(jié)催化劑稀釋比,將CO 轉(zhuǎn)化率控制在10%以下,確保反應(yīng)處于動力學(xué)控制范圍內(nèi)。隨后對催化劑進(jìn)行反應(yīng)時間為80 h 的穩(wěn)定性測試,并排除催化劑內(nèi)外擴(kuò)散影響[10-11]。最后,在動力學(xué)區(qū)間內(nèi)測定催化劑的表觀反應(yīng)活化能(Ea),并通過收變原料氣中H2/CO 比值,求取各產(chǎn)物對應(yīng)的H2和CO 反應(yīng)級數(shù)。

    外擴(kuò)散影響可以通過收變催化劑表面氣流的邊界層厚度排除,為了控制體積空速不收變,同時收變催化劑裝填體積與氣體體積流量。

    在固定的反應(yīng)空速下,測定不同氣體流速中的CO 轉(zhuǎn)化率,并利用Frossling 系數(shù)[11]得到球形顆粒大小、傳質(zhì)系數(shù)(kc)與邊界層厚度之間的關(guān)系。當(dāng)傳質(zhì)系數(shù)變化小于一定值時,可以斷定在該反應(yīng)條件下,外擴(kuò)散影響已經(jīng)被排除。

    通過收變催化劑顆粒的尺寸,測定不同粒徑催化劑上的CO 轉(zhuǎn)化率,可以通過Weisz-Prater 判據(jù)判斷是否排除內(nèi)擴(kuò)散影響[12]。

    1.6 動力學(xué)模型

    本實(shí)驗(yàn)采用冪指數(shù)模型擬合動力學(xué)數(shù)據(jù),得到反應(yīng)表觀活化能(Ea)以及不同產(chǎn)物對應(yīng)的H2(m)和CO(n)[式(5)]反應(yīng)級數(shù),不同溫度下的反應(yīng)級數(shù)能夠反映反應(yīng)溫度和壓力對反應(yīng)速率的影響。此外,冪指數(shù)模型對于該反應(yīng)的反應(yīng)器設(shè)計和工程放大有重要的指導(dǎo)意義[13]。

    式中,rA是產(chǎn)物A 的產(chǎn)生速率;速率常數(shù)k 可由Arrhenius 方程表示

    將式(6)代入式(5)中,假設(shè)pH2和pCO為常量,可得

    分別假設(shè)pH2和pCO為常量,則式(7)可變形為式(8)和式(9)。

    通過對lnrA與反應(yīng)溫度倒數(shù)(533~573 K)、lnpH2和lnpCO作圖,可得Ea、m 和n 的值。

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

    2.1 Fe20/SiO2 與Fe20-Mn1.0/SiO2 催化劑活性比較

    考察了Fe20/SiO2與Fe20-Mn1.0/SiO2催化劑上的FTO 反應(yīng)催化性能(表1)。產(chǎn)物選擇性計算均不考慮CO2。對Fe 基催化劑而言,F(xiàn)e0[5]、FeCx[14]和Fe3O4[15-16]均被認(rèn)為是可能的費(fèi)-托反應(yīng)活性中心。Mn20/SiO2對FTO 幾乎沒有活性,可以斷定MnOx不是FTO 反應(yīng)的活性國屬。

    表1 Fe-Mn/SiO2 催化劑上的氣相產(chǎn)物選擇性及反應(yīng)活性Table 1 Effects of MnOx on reaction rates and selectivity of all products over Fe-Mn/SiO2 catalysts

    如表1 所示,在Fe20-Mn1.0/SiO2催化劑上,CO轉(zhuǎn)化率和STY 均顯著高于Fe20/SiO2。由于MnOx助劑的添加,CO 轉(zhuǎn)化率由9.5%升至16.1%,C2~C4烯烴選擇性由24.7%升至28.1%,甲烷選擇性無顯著變化。值得注意的是,F(xiàn)e20-Mn1.0/SiO2催化劑上STY 比Fe20/SiO2催化劑上高70%。

    2.2 催化劑表面性質(zhì)

    通過CO-TPD 考察了Fe20/SiO2和Fe20-Mn1.0/ SiO2催化劑對CO 的吸附性能。圖1 為催化劑Fe20/SiO2和Fe20-Mn1.0/SiO2的CO-TPD 譜圖,催化劑的CO 脫附峰主要分布在高(800~1000 K)、低(350~450 K)兩個溫度范圍內(nèi)。低溫的脫附峰對應(yīng)的是催化劑上分子形態(tài)吸附的CO,位于高溫的脫附峰對應(yīng)的是催化劑上解離吸附的CO[17-18]。Fe20/SiO2催化劑(圖1 曲線a)有兩個分別位于376和409 K 的低溫脫附峰。此外,在931 K 處有一個較寬的高溫脫附峰。對Fe20-Mn1.0/SiO2催化劑(圖1 曲線b)而言,低溫區(qū)域的脫附峰面積與Fe20/SiO2催化劑相比變化不大,但其高溫區(qū)域的脫附峰面積明顯增加,并且在988 K 出現(xiàn)新的脫附峰,表明在Fe20-Mn1.0/SiO2催化劑上,CO 的解離吸附增強(qiáng),出現(xiàn)了更強(qiáng)的CO 吸附位。通過對比可知,添加MnOx助劑后,F(xiàn)e20-Mn1.0/SiO2催化劑在804 K 出現(xiàn)寬化弱峰,在988 K 處出現(xiàn)較強(qiáng)的脫附峰,同時931 K的脫附峰移至較高溫度937 K 處。該結(jié)果表明,F(xiàn)e基催化劑中添加MnOx助劑不但增強(qiáng)了催化劑表面上CO 的吸附,而且促進(jìn)了CO 在催化劑表面的吸附解離。因此,添加適量MnOx助劑可顯著提高Fe基催化劑上的CO 轉(zhuǎn)化率。

    圖1 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑的CO-TPD 譜圖Fig.1 CO-TPD profiles of Fe20/SiO2 and Fe20-Mn1.0/SiO2 catalysts

    本文使用CO2-TPD 考察MnOx助劑的添加對催化劑表面堿性的影響(圖2)。Fe20/SiO2催化劑(圖2 曲線a)在387 和478 K 處有兩個脫附峰,分別對應(yīng)于CO2的弱吸附和中等強(qiáng)度吸附。添加MnOx助劑后,F(xiàn)e20-Mn1.0/SiO2催化劑(圖2 曲線b)上兩脫附峰的脫附溫度均發(fā)生變化,387 K 處脫附峰移至397 K,而478 K 處脫附峰則向低溫處移至464 K。值得注意的是,添加MnOx助劑后這兩個脫附峰的強(qiáng)度明顯增強(qiáng),表明堿性位數(shù)量有所增加。此外,添加MnOx助劑后,催化劑在688 K 處出現(xiàn)新的脫附峰,對應(yīng)于CO2的強(qiáng)吸附[19]。上述結(jié)果表明,MnOx助劑的加入提高了Fe 基催化劑的表面堿性,不僅增加了堿性位的數(shù)量,還產(chǎn)生了新的強(qiáng)堿性位。這與文獻(xiàn)中報道的錳氧化物相對于鐵氧化物擁有更強(qiáng)的堿性是一致的[20]。催化劑表面堿性的增強(qiáng)能夠促進(jìn)活性位對CO 的吸附,這與CO-TPD 的實(shí)驗(yàn)結(jié)果相符。

    圖2 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑的CO2-TPD 譜圖Fig.2 CO2-TPD profiles of Fe20/SiO2 and Fe20-Mn1.0/SiO2 catalysts

    2.3 動力學(xué)研究

    2.3.1 穩(wěn)定性測試與外擴(kuò)散排除 在進(jìn)行動力學(xué)研究之前,首先通過調(diào)節(jié)稀釋比(稀釋劑為SiO2),控制CO 轉(zhuǎn)化率低于10%,確保反應(yīng)在動力學(xué)控制范圍內(nèi)進(jìn)行。國后,對 Fe20-Mn1.0/SiO2和Fe20/SiO2催化劑進(jìn)行了80 h 的FTO 反應(yīng)穩(wěn)定性測試。CO 轉(zhuǎn)化率、CH4選擇性以及C2~C4烯烴選擇性均無明顯變化。反應(yīng)條件為:反應(yīng)溫度593 K,反應(yīng)壓力2.0 MPa,反應(yīng)氣H2/CO 比為1:1,體積空速12000 h-1。

    排除外擴(kuò)散影響可以通過收變不同的氣體流速與催化劑裝填量(固定GHSV),測定不同氣體流速下的CO 轉(zhuǎn)化率(圖3 曲線a)來判定。反應(yīng)條件:溫度593 K,壓力2.0 MPa,反應(yīng)氣體H2/CO =1:1,GHSV = 12000 h-1,反應(yīng)氣流速分別為45、60、75、90、105 ml·min-1。由圖可知,在Fe20-Mn1.0/SiO2催化劑上氣體流速從 45 ml·min-1升高至 75 ml·min-1的過程中,CO 轉(zhuǎn)化率從7.8%上升到11.5%。國而,當(dāng)氣體流速高于75 ml·min-1時,轉(zhuǎn)化率升高幅度變小。通過計算不同流速時相應(yīng)的傳質(zhì)系數(shù)(kc)(表2),可知當(dāng)反應(yīng)氣進(jìn)氣流速高于75 ml·min-1時,傳質(zhì)系數(shù)趨于恒定(差值小于1%),即表示Fe20-Mn1.0/SiO2催化劑上外擴(kuò)散影響已被排除。因此,在后續(xù)實(shí)驗(yàn)中,選用105 ml·min-1作為動力學(xué)測試的進(jìn)氣流速。

    圖3 Fe20-Mn1.0/SiO2催化劑上CO 轉(zhuǎn)化率隨反應(yīng)氣體流速(a)和催化劑粒徑的變化(b),0.15~0.18 nm 的Fe20/SiO2 催化劑上CO 轉(zhuǎn)化率隨反應(yīng)氣體流速的變化(c)Fig.3 Dependence of CO conversion on flow rate (a) and particle size (b) over Fe20-Mn1.0/SiO2 catalyst, dependence of CO conversion on flow rate over Fe20/SiO2 catalyst with 0.15—0.18 nm(c)

    內(nèi)擴(kuò)散排除可以通過測定不同顆粒尺寸催化劑上的CO 轉(zhuǎn)化率[21],得到CO 轉(zhuǎn)化率與催化劑顆粒尺寸的關(guān)系曲線(圖3 曲線b)來判定。本實(shí)驗(yàn)中采用了5 種不同粒徑的Fe20-Mn1.0/SiO2催化劑,其粒徑分別是0.42~0.84、0.25~0.42、0.18~0.25、0.15~0.18 和0.13~0.15 nm。由圖3 曲線b 可知,當(dāng)催化劑顆粒小于0.25 nm 時,減小催化劑顆粒尺寸,CO 轉(zhuǎn)化率快速上升;而當(dāng)催化劑顆粒大于0.25 nm 時,CO 轉(zhuǎn)化率基本恒定,計算得到Weisz-Prater 判據(jù)CWP= 1.3×10-4。根據(jù)定義,當(dāng)CWP? 1 時,F(xiàn)e20-Mn1.0/SiO2催化劑上內(nèi)擴(kuò)散影響可忽略不計。

    表2 Fe20-Mn1.0/SiO2 催化劑上氣體流速、催化劑裝填量、CO 轉(zhuǎn)化率與相應(yīng)kc 值Table 2 Gas flow rate, corresponding mass of loaded catalysts, CO conversion and kc value over Fe20-Mn1.0/SiO2

    在Fe20/SiO2催化劑上(圖3 曲線c),當(dāng)選用0.15~0.18 nm 催化劑,105 ml·min-1流速時,均滿足判定條件(Weisz-Prater 判據(jù)CWP= 8.3×10-5)。因此,在后續(xù)動力學(xué)測試中,動力學(xué)測試條件為:反應(yīng)溫度533 K,反應(yīng)壓力2.0 MPa,體積空速為12000 h-1,催化劑裝填量0.216 g,氣體流量105 ml·min-1。

    2.3.2 不同產(chǎn)物生成的反應(yīng)活化能 在533、553、573、593 K 的溫度下,測定了 Fe20/SiO2和Fe20-Mn1.0/SiO2催化劑上的CO 轉(zhuǎn)化率及產(chǎn)物生成速率。根據(jù)反應(yīng)速率與溫度(T)的倒數(shù)之間的關(guān)系(圖4),可以得到各產(chǎn)物及CO 總反應(yīng)的活化能數(shù)據(jù)(表3)。

    表3 各產(chǎn)物的活化能與反應(yīng)級數(shù)Table 3 Apparent activation energies and reaction orders for different products

    圖4 Fe20/SiO2 及Fe20-Mn1.0/SiO2 催化劑上各產(chǎn)物反應(yīng)活化能Fig.4 Arrhenius diagrams of products over Fe20/SiO2 and Fe20-Mn1.0/SiO2 (Reaction conditions: 2.0 MPa, H2/CO =1:1, GHSV=12000 h-1, 20 h, 533—593 K)

    在Fe20/SiO2催化劑上,甲烷的活化能為64.0 kJ·mol-1;而在Fe20-Mn1.0/SiO2催化劑上,甲烷的活化能增加至75.5 kJ·mol-1。甲烷活化能的增加表明 MnOx的添加可以抑制甲烷的生成。根據(jù)CO2-TPD 的結(jié)果,MnOx助劑的加入能夠提高Fe 基催化劑的表面堿性,從而減少甲烷的生成。在Fe20/SiO2和Fe20-Mn1.0/SiO2催化劑上,乙烯和丙烯活化能比較接近,均在72.0~80.0 kJ·mol-1之間。但是,在Fe20-Mn1.0/SiO2催化劑上丁烯的活化能(84.9 kJ·mol-1)遠(yuǎn)高于Fe20/SiO2催化劑上的丁烯的活化能(51.8 kJ·mol-1),說明MnOx助劑可能能夠抑制長鏈烯烴的生成。

    如表3 所示,F(xiàn)e20/SiO2催化劑上乙烷、丙烷和丁烷的活化能(62.4、40.5、45.7 kJ·mol-1)遠(yuǎn)低于該催化劑上乙烯、丙烯和丁烯的活化能(79.2、74.7、51.8 kJ·mol-1)。同樣,在Fe20-Mn1.0/SiO2催化劑上,C2~C4烷烴活化能(81.5、62.1、49.3 kJ·mol-1)也低于C2~C4烯烴的活化能(74.8、72.2、84.9 kJ·mol-1)。國而,在Fe20-Mn1.0/SiO2催化劑上,乙烷、丙烷和丁烷的活化能均高于Fe20/SiO2催化劑上相應(yīng)產(chǎn)物的活化能。根據(jù)CO2-TPD 結(jié)果,添加MnOx助劑后,催化劑表面堿性增強(qiáng),抑制了C2~C4烯烴的再吸附;根據(jù)CO-TPD 結(jié)果,MnOx助劑促進(jìn)了CO 的吸附,更多的活性位被CO*占據(jù),也有效地抑制了C2~C4烯烴的吸附。因此,C2~C4烯烴再吸附后二次加氫生成相應(yīng)烷烴的過程也被抑制,導(dǎo)致在Fe20-Mn1.0/SiO2催化劑上乙烷、丙烷和丁烷產(chǎn)率的下降[22-25]。另外,在Fe20/SiO2催化劑上,C1~C4烷烴的活化能隨碳鏈增長而下降,說明在Fe20/SiO2催化劑上烷烴的生成速率隨著碳數(shù)增加而增加。

    Fe20-Mn1.0/SiO2催化劑上CO2的活化能(約為145.2 kJ·mol-1)高于Fe20/SiO2催化劑上CO2的活化能(約為113.8 kJ·mol-1),說明添加MnOx助劑抑制了CO2的生成。在單鐵催化劑上總反應(yīng)的活化能為65.9 kJ·mol-1,比Dry 等[26]所報道的在融鐵催化劑上的總包活化能(71.0 kJ·mol-1)低。添加MnOx以后,在Fe20-Mn1.0/SiO2催化劑上總反應(yīng)的活化能下降為60.4 kJ·mol-1,比Yang 等[27]所報道的工業(yè)Fe-Mn 催化劑上的總反應(yīng)活化能低了15.5 kJ·mol-1。Fe20-Mn1.0/SiO2催化劑總反應(yīng)活化能的降低可能是由于CO 解離吸附的增強(qiáng),導(dǎo)致CO 轉(zhuǎn)化率升高。

    2.3.3 各產(chǎn)物的H2與CO 反應(yīng)級數(shù) 在533 K 下,比較了Fe20/SiO2和Fe20-Mn1.0/SiO2催化劑上各產(chǎn)物生成及總反應(yīng)的H2和CO 反應(yīng)級數(shù)的差異(圖5~圖11),不同CO 分壓反應(yīng)氣配比為pH2= 0.6 MPa,pCO= 0.15~1.2 MPa;不同H2分壓反應(yīng)氣配比為pCO= 0.36 MPa,pH2= 0.18~1.44 MPa。所得反應(yīng)級數(shù)列于表3。由圖5~圖9 和表3 所列數(shù)據(jù)可知,不同碳?xì)洚a(chǎn)物的生成速率,如甲烷(圖5)、丙烯(圖7)、乙烷(圖8)和丙烷(圖9),都隨H2分壓的上升而增加。除了乙烯和丙烯的生成速率,其他碳?xì)洚a(chǎn)物的生成速率均隨CO 分壓的上升而下降。由于催化劑表面被強(qiáng)吸附的CO 或解離C*覆蓋,CO 分壓的增加抑制了H2在催化劑表面的吸附及解離,因此大部分碳?xì)洚a(chǎn)物的CO 反應(yīng)級數(shù)為負(fù)值。反之,H2分壓增加能夠促進(jìn)反應(yīng)速率。值得注意的是,與其他產(chǎn)物相反,在533 K 時,F(xiàn)e20/SiO2上的乙烯和丙烯生成速率隨CO 分壓的增加均呈上升趨勢,這一現(xiàn)象可能是由α-烯烴的再吸附所導(dǎo)致[28-29]。當(dāng)CO分壓增加時,CO 或解離吸附C*的覆蓋度上升,占據(jù)了更多的活性位,并阻礙了乙烯和丙烯的再吸附及其后二次加氫和鏈增長生成相應(yīng)烷烴和長鏈烯烴的過程,乙烯和丙烯繼而在催化劑表面脫附后直接形成氣相產(chǎn)物。所以CO 分壓的增加提升了乙烯和丙烯的產(chǎn)生速率。

    圖5 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上CH4 的生成速率與H2 和CO 分壓的關(guān)系Fig.5 Dependence of reaction rates of CH4 formation on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    圖6 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上C2H4 的生成速率與H2 和CO 分壓的關(guān)系Fig.6 Dependence of reaction rates of C2H4 formation on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    圖7 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上C3H6 的生成速率與H2 和CO 分壓的關(guān)系Fig.7 Dependence of reaction rates of C3H6 formation on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    添加MnOx之后,催化劑表面堿性增強(qiáng),受CO在催化劑表面競爭吸附的影響,C2~C4烯烴的再吸附被抑制,進(jìn)而抑制了C2~C4烷烴的產(chǎn)生,導(dǎo)致烯烴的CO 反應(yīng)級數(shù)減小。特別是乙烯的CO 反應(yīng)級數(shù),在含MnOx催化劑上,由單鐵催化劑上的0.36減小到-0.54。此外,在Fe20-Mn1.0/SiO2催化劑上,CO2的生成速率受CO 分壓影響減弱,由單鐵催化劑上的0.66 下降為0.11。推斷主要是由于CO2的生成主要依賴于分子態(tài)的CO,而在Fe20-Mn1.0/SiO2催化劑上,CO 更傾向于解離吸附而非分子態(tài)吸附,因此CO2的生成受CO 分壓影響較小。

    與Fe20/SiO2催化劑上各產(chǎn)物的H2反應(yīng)級數(shù)比較,添加MnOx助劑后各產(chǎn)物的H2反應(yīng)級數(shù)沒有發(fā)生明顯變化。相反,碳?xì)洚a(chǎn)物的CO 反應(yīng)級數(shù)明顯下降,說明MnOx助劑的添加嚴(yán)重影響了CO 在催化劑表面的吸附。TPD 實(shí)驗(yàn)結(jié)果證明,MnOx的添加增強(qiáng)了CO 的吸附。因此,在FTO 反應(yīng)中,強(qiáng)吸附的CO 占據(jù)更多活性位,導(dǎo)致了碳?xì)洚a(chǎn)物CO 反應(yīng)級數(shù)在Fe20-Mn1.0/SiO2催化劑上均有所下降。

    基于上文中的動力學(xué)結(jié)果,討論Fe 基催化劑上的FTO 反應(yīng)機(jī)理。FTO 與傳統(tǒng)費(fèi)-托反應(yīng)在機(jī)理上有相似之處[30]:①反應(yīng)速度取決于CO 的解離活化;②抑制CO2的生成有利于提高碳的利用效率。但是,由于FTO 反應(yīng)的產(chǎn)物是低碳烯烴,與費(fèi)-托反應(yīng)相比,該反應(yīng)更取決于以下因素:①鐵活性位上低碳烯烴產(chǎn)物的鏈終止與脫附;②抑制低碳烯烴的再吸附與二次加氫。圖12 說明了MnOx助劑對上述兩因素的影響作用,其中虛線所代表的步驟在添加MnOx助劑后被抑制。解離吸附的碳物種在催化劑表面進(jìn)行偶聯(lián)反應(yīng)(C1*+C1*和C2*+C1*)后產(chǎn)生C2H4*和C3H6*。同時,C2H4*和C3H6*繼續(xù)進(jìn)行鏈增長反應(yīng)生成C3H6*和C4H8*。乙烯和丙烯也可能通過再吸附后進(jìn)行二次加氫生成相應(yīng)的烷烴。甲烷則通過C1*直接加氫產(chǎn)生。根據(jù)動力學(xué)結(jié)果,MnOx的添加提高了Fe20-Mn1.0/SiO2催化劑上C1~C4烷烴和丁烯的活化能,表明MnOx可能通過抑制C1~C4烷烴及長鏈烯烴的產(chǎn)生而提高了低碳烯烴的選擇性?;赥PD 表征結(jié)果,這一現(xiàn)象可能與MnOx助劑增強(qiáng)了催化劑表面堿性,進(jìn)而抑制了烯烴的二次吸附與隨后的二次反應(yīng)相關(guān)。

    總之,隨MnOx助劑的添加,F(xiàn)e 基催化劑的催化活性和STY 均顯著增加。MnOx對氧化鐵相表面堿性的增強(qiáng),在FTO 體系中起到重要作用。國而,MnOx助劑對氧化鐵相碳化過程及低碳烯烴二次反應(yīng)的影響,需要進(jìn)一步研究,從而更加深刻理解MnOx對FTO 反應(yīng)鐵基催化劑的助催化作用。

    圖8 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上C2H6 的生成速率與H2 和CO 分壓的關(guān)系Fig.8 Dependence of reaction rates of C2H6 formation on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    圖9 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上C3H8 的生成速率與H2 和CO 分壓的關(guān)系Fig.9 Dependence of reaction rates of C3H8 formation on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    圖10 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上CO2 的生成速率與H2 和CO 分壓的關(guān)系Fig.10 Dependence of reaction rates of CO2 formation on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    圖11 Fe20/SiO2 和Fe20-Mn1.0/SiO2 催化劑上CO 的消耗速率與H2 和CO 分壓的關(guān)系Fig.11 Dependence of reaction rates of CO consumption on partial pressure of H2 and CO over Fe20/SiO2 and Fe20-Mn1.0/SiO2

    圖12 費(fèi)-托合成中可能的鏈增長與烯烴二次反應(yīng)路徑Fig.12 Chain growth pathways and possible secondary reactions of olefins in Fischer-Tropsch synthesis

    3 結(jié) 論

    采用等體積浸漬法制備了 Fe20/SiO2與Fe20-Mn1.0/SiO2催化劑并考評了其在FTO 反應(yīng)上的催化活性。通過CO 和CO2-TPD 比較了Fe20/SiO2和Fe20-Mn1.0/SiO2催化劑的吸附行為。MnOx助劑能夠增加Fe 基催化劑的表面堿性,增強(qiáng)CO 的解離吸附。

    CO 加氫反應(yīng)主要發(fā)生在鐵活性位上,MnOx助劑的添加提升了FTO 反應(yīng)在Fe 基催化劑上的反應(yīng)速率及低碳烯烴選擇性。Fe20-Mn1.0/SiO2催化劑上C2和C3烯烴生成以及CO 轉(zhuǎn)化的活化能均有所下降。在533 K 時,F(xiàn)e20/SiO2催化劑上碳?xì)洚a(chǎn)物(除乙烯和丙烯外)的CO 反應(yīng)級數(shù)均為負(fù)數(shù)。添加MnOx助劑后,CO 反應(yīng)級數(shù)均顯著降低,特別是丙烯的CO 反應(yīng)級數(shù),由正值降為負(fù)值。國而,在兩種催化劑上,H2的反應(yīng)級數(shù)均為正值且十分接近。可見,CO 加氫過程符合Langmuir-Hinshelwood 機(jī)理。綜上所述,MnOx助劑能夠增強(qiáng)Fe 基催化劑表面堿性,進(jìn)而促進(jìn)表面CO 解離吸附,提高CO 轉(zhuǎn)化率。同時,表面堿性的增加抑制了烯烴在催化劑表面的再吸附及其后續(xù)二次反應(yīng),從而提高了低碳烯烴的生成速率。

    符 號 說 明

    A ——指前因子,g·(kg cat)-1·h-1

    Ea——反應(yīng)的活化能,kJ·mol-1

    k ——反應(yīng)速率常數(shù),g·(kg cat)-1·h-1

    m ——H2反應(yīng)級數(shù)

    n ——CO 反應(yīng)級數(shù)

    pCO——CO 分壓,MPa

    pH2——H2分壓,MPa

    R ——?dú)怏w常數(shù),J·K-1·mol-1

    ri——產(chǎn)物i 的生成速率,g·(kg cat)-1·h-1

    STY ——產(chǎn)物的時空收率,g·(kg cat)-1·h-1

    T ——反應(yīng)溫度,K

    [1]Dong Li(董麗), Yang Xueping(楊學(xué)萍).New advances in direct production of light olefins from syngas [J].Petrochemical Technology (石油化工), 2012, 41(10): 1201-1206.

    [2]Fu D, Dai W, Xu X, et al.Probing the structure evolution of iron-based Fischer-Tropsch to produce olefins by operando Raman spectroscopy [J].ChemCatChem, 2015,7(5): 752-756.

    [3]Steynberg A, Dry M, Davis B, et al.Fischer-Tropsch reactors [J].Studies in Surface Science and Catalysis, 2004, 152: 64-195.

    [4]Davis B H.Fischer-Tropsch synthesis: comparison of performances of iron and cobalt catalysts [J].Industrial & Engineering Chemistry Research, 2007, 46(26): 8938-8945.

    [5]de Smit E, Weckhuysen B M.The renaissance of iron-based Fischer-Tropsch synthesis: on the multifaceted catalyst deactivation behaviour [J].Chem.Soc.Rev., 2008, 37(12): 2758-2781.

    [6]Torres Galvis H M, Bitter J H, Khare C B, et al.Supported iron nanoparticles as catalysts for sustainable production of lower olefins [J].Science, 2012, 335(6070): 835-838.

    [7]Torres Galvis H M, de Jong K P.Catalysts for production of lower olefins from synthesis gas: a review [J].ACS Catalysis, 2013, 3(9): 2130-2149.

    [8]Luque R, de la Osa A R, Campelo J M, et al.Design and development of catalysts for biomass-to-liquid-Fischer-Tropsch (BTL-FT) processes for biofuels production [J].Energy & Environmental Science, 2012, 5(1): 5186-5202.

    [9]Das D, Ravichandran G, Chakrabarty D K.Conversion of syngas to light olefins over silicalite-1 supported iron and cobalt catalysts: effect of manganese addition [J].Catalysis Today, 1997, 36(3): 285-293.

    [10]Malessa R, Baerns M.Iron/manganese oxide catalysts for Fischer-Tropsch synthesis (Ⅳ): Activity and selectivity [J].Industrial & Engineering Chemistry Research, 1988, 27(2): 279-283.

    [11]Frossling N.The evaporation of falling drops [J].Gerlande Beitr Geophys, 1938, 52: 170-216.

    [12]Weisz P, Hicks J.The behaviour of porous catalyst particles in view of internal mass and heat diffusion effects [J].Chemical Engineering Science, 1962, 17(4): 265-275.

    [13]Su Junjie(蘇俊杰), Mao Wei(茅威), Yang Zhen(楊震), Xu Jing(徐晶), Han Yifan(韓一帆).Kinetics of CoCu/SiO2for synthesis of lower carbon mixed alcohols directly from syngas [J].CIESC Journal(化工學(xué)報), 2014, 65(1): 143-151.

    [14]Shroff M D, Kalakkad D S, Coulter K E, et al.Activation of precipitated iron Fischer-Tropsch synthesis catalysts [J].Journal of Catalysis, 1995, 156(2): 185-207.

    [15]Butt J.Carbide phases on iron-based Fischer-Tropsch synthesis catalysts (Ⅰ): Characterization studies [J].Catalysis Letters, 1990, 7(1-4): 61-81.

    [16]Kuivila C, Stair P, Butt J.Compositional aspects of iron Fischer-Tropsch catalysts: an XPS/reaction study [J].Journal of Catalysis, 1989, 118(2): 299-311.

    [17]Yang Z, Pan X, Wang J, et al.FeN particles confined inside CNT for light olefin synthesis from syngas: effects of Mn and K additives [J].Catalysis Today, 2012, 186(1): 121-127.

    [18]Fierro J L G.Cu-promoted Fe2O3/MgO-based Fischer-Tropsch catalysts of biomass-derived syngas [J].Industrial & Engineering Chemistry Research, 2015, 54(3): 911-921

    [19]Tang L, Song C, Li M, et al.Study of K/Mn-MgO supported Fe catalysts with Fe(CO)5and Fe(NO3)3as precursors for CO hydrogenation to light alkenes [J].Chinese Journal of Chemistry, 2013, 31(10): 1263-1268.

    [20]Jensen K, Massoth F.Studies on iron-manganese oxide carbon monoxide catalysts (Ⅱ): Carburization and catalytic activity [J].Journal of Catalysis, 1985, 92(1): 109-118.

    [21]Mei D, Rousseau R, Kathmann S M, et al.Ethanol synthesis from syngas over Rh-based/SiO2catalysts: a combined experimental and theoretical modeling study [J].Journal of Catalysis, 2010, 271(2): 325-342.

    [22]Cheng J, Hu P, Ellis P, et al.A DFT study of the chain growth probability in Fischer-Tropsch synthesis [J].Journal of Catalysis, 2008, 257(1): 221-228.

    [23]Cheng J, Hu P, Ellis P, et al.A first-principles study of oxygenates on Co surfaces in Fischer-Tropsch synthesis [J].The Journal of Physical Chemistry C, 2008, 112(25): 9464-9473.

    [24]Cheng J, Hu P, Ellis P, et al.Some understanding of Fischer-Tropschsynthesis from density functional theory calculations [J].Topics in Catalysis, 2010, 53(5/6): 326-337.

    [25]Cheng J, Song T, Hu P, et al.A density functional theory study of the α-olefin selectivity in Fischer-Tropsch synthesis [J].Journal of Catalysis, 2008, 255(1): 20-28.

    [26]Dry M, Shingles T, Boshoff L.Rate of the Fischer-Tropsch reaction over iron catalysts [J].Journal of Catalysis, 1972, 25(1): 99-104.

    [27]Yang J, Liu Y, Chang J, et al.Detailed kinetics of Fischer-Tropsch synthesis on an industrial Fe-Mn catalyst [J].Industrial & Engineering Chemistry Research, 2003, 42(21): 5066-5090.

    [28]Ahón V R, Costa Jr E F, Monteagudo J E, et al.A comprehensive mathematical model for the Fischer-Tropsch synthesis in well-mixed slurry reactors[J].Chemical Engineering Science, 2005, 60(3): 677-694.

    [29]Tsubaki N, Yoshii K, Fujimoto K.Anti-ASF distribution of Fischer-Tropsch hydrocarbons in supercritical-phase reactions [J].Journal of Catalysis, 2002, 207(2): 371-375.

    [30]Govender N S, Botes F G, de Croon M H J M, et al.Mechanistic pathway for C2+hydrocarbons over an Fe/K catalyst [J].Journal of Catalysis, 2014, 312: 98-107.

    猜你喜歡
    級數(shù)丙烯助劑
    錢愛康
    《共生》主題系列作品
    大眾文藝(2022年16期)2022-09-07 03:08:04
    農(nóng)藥減量增效助劑應(yīng)用技術(shù)
    苯丙烯菌酮
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    幾個常數(shù)項(xiàng)級數(shù)的和
    液化氣中的丙烯含有的雜質(zhì)對丙烯聚合反應(yīng)的影響
    塑料助劑的應(yīng)用
    p級數(shù)求和的兩種方法
    橡膠助劑綠色化勢在必行
    精品国产国语对白av| 搡老乐熟女国产| 草草在线视频免费看| 一区二区av电影网| 我的女老师完整版在线观看| 成人国产av品久久久| 最近中文字幕2019免费版| 国产男女超爽视频在线观看| 日本-黄色视频高清免费观看| 成人毛片60女人毛片免费| 亚洲精品自拍成人| 欧美日韩视频精品一区| 美女主播在线视频| 亚洲综合色惰| 亚洲精品乱久久久久久| 男女国产视频网站| 国精品久久久久久国模美| 99九九在线精品视频| 日本欧美国产在线视频| 国产日韩欧美视频二区| 最近2019中文字幕mv第一页| 亚洲国产色片| 午夜精品国产一区二区电影| 亚洲成人一二三区av| 狂野欧美白嫩少妇大欣赏| 国产精品嫩草影院av在线观看| 国产精品久久久久久av不卡| 国产精品一区www在线观看| 欧美日韩亚洲高清精品| 大又大粗又爽又黄少妇毛片口| 精品少妇黑人巨大在线播放| 女性生殖器流出的白浆| 最近中文字幕2019免费版| 欧美日韩视频精品一区| 日本免费在线观看一区| 国产精品久久久久久久电影| 久久人妻熟女aⅴ| 三级国产精品欧美在线观看| 国产精品国产三级专区第一集| 日韩伦理黄色片| 国产一区二区三区av在线| 亚洲欧洲国产日韩| 国产成人91sexporn| 99热国产这里只有精品6| 日产精品乱码卡一卡2卡三| 在线看a的网站| 99国产综合亚洲精品| 免费看光身美女| 午夜影院在线不卡| 内地一区二区视频在线| 韩国高清视频一区二区三区| 欧美xxxx性猛交bbbb| 欧美 亚洲 国产 日韩一| 久久韩国三级中文字幕| 久久青草综合色| 精品一区在线观看国产| 久久人人爽人人爽人人片va| 日日撸夜夜添| 亚洲天堂av无毛| 亚洲精品日韩在线中文字幕| 人体艺术视频欧美日本| 国产在线视频一区二区| 赤兔流量卡办理| 国产日韩欧美亚洲二区| 欧美亚洲日本最大视频资源| 97精品久久久久久久久久精品| 丝袜美足系列| 日韩成人av中文字幕在线观看| 激情五月婷婷亚洲| 久久久精品94久久精品| av天堂久久9| 丰满迷人的少妇在线观看| 91精品伊人久久大香线蕉| 简卡轻食公司| 欧美丝袜亚洲另类| 91成人精品电影| 国产午夜精品一二区理论片| 看十八女毛片水多多多| 成人手机av| 高清av免费在线| 久久97久久精品| 国产老妇伦熟女老妇高清| 特大巨黑吊av在线直播| 九九爱精品视频在线观看| av视频免费观看在线观看| 日韩人妻高清精品专区| 视频中文字幕在线观看| 国产亚洲精品第一综合不卡 | 久久99蜜桃精品久久| 丰满迷人的少妇在线观看| 成人二区视频| 全区人妻精品视频| 日韩强制内射视频| 亚洲人成网站在线播| 成人亚洲欧美一区二区av| 成人国语在线视频| 51国产日韩欧美| 亚洲精品一二三| 欧美精品一区二区免费开放| 国产av国产精品国产| 中文字幕免费在线视频6| 国产免费现黄频在线看| 欧美精品国产亚洲| 国产精品一二三区在线看| 国产一区二区三区综合在线观看 | 免费观看无遮挡的男女| 国产爽快片一区二区三区| 成年女人在线观看亚洲视频| 一区二区三区精品91| 日本vs欧美在线观看视频| 国产成人a∨麻豆精品| 高清在线视频一区二区三区| 日韩av不卡免费在线播放| 中文字幕最新亚洲高清| 自线自在国产av| 一级毛片黄色毛片免费观看视频| 久久午夜综合久久蜜桃| 亚洲av在线观看美女高潮| av.在线天堂| 有码 亚洲区| 黄色视频在线播放观看不卡| 精品亚洲成国产av| 国产精品99久久久久久久久| 免费大片黄手机在线观看| 精品熟女少妇av免费看| 久久午夜综合久久蜜桃| 亚洲少妇的诱惑av| 亚洲一区二区三区欧美精品| 亚洲高清免费不卡视频| 欧美变态另类bdsm刘玥| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲精品久久久com| 人妻 亚洲 视频| 91国产中文字幕| 高清不卡的av网站| 国产av一区二区精品久久| 哪个播放器可以免费观看大片| av卡一久久| 在线观看一区二区三区激情| 欧美日韩视频精品一区| 女性被躁到高潮视频| 国产不卡av网站在线观看| 国产一区二区在线观看日韩| 中文乱码字字幕精品一区二区三区| 久久毛片免费看一区二区三区| 久久99精品国语久久久| 一级毛片我不卡| 日韩精品免费视频一区二区三区 | 精品久久久久久电影网| 一边摸一边做爽爽视频免费| 大陆偷拍与自拍| 亚洲国产欧美在线一区| 精品国产乱码久久久久久小说| 久久97久久精品| 又黄又爽又刺激的免费视频.| 亚洲av欧美aⅴ国产| 欧美日韩视频高清一区二区三区二| 自线自在国产av| 寂寞人妻少妇视频99o| 国产视频首页在线观看| 特大巨黑吊av在线直播| 日韩一本色道免费dvd| 久久精品久久精品一区二区三区| 中文字幕免费在线视频6| 成年美女黄网站色视频大全免费 | 精品人妻熟女毛片av久久网站| 人人妻人人澡人人看| 久久久久网色| 亚洲av二区三区四区| 少妇的逼好多水| 国产精品久久久久久久久免| 22中文网久久字幕| 热re99久久国产66热| 国产日韩欧美在线精品| 国产极品粉嫩免费观看在线 | 国产欧美亚洲国产| 在线观看三级黄色| 午夜精品国产一区二区电影| a级毛片在线看网站| 久久久久久久精品精品| 久久久国产欧美日韩av| 亚洲天堂av无毛| 午夜老司机福利剧场| 成人午夜精彩视频在线观看| 亚洲国产精品专区欧美| a级毛片黄视频| 亚洲av男天堂| 性色av一级| 夫妻性生交免费视频一级片| 成人漫画全彩无遮挡| 99久久人妻综合| 国产成人精品久久久久久| 18禁在线播放成人免费| 又黄又爽又刺激的免费视频.| 亚洲精品,欧美精品| 大片免费播放器 马上看| 亚洲三级黄色毛片| 日韩精品有码人妻一区| 免费观看的影片在线观看| 97超碰精品成人国产| 少妇丰满av| 国产精品嫩草影院av在线观看| 91精品三级在线观看| 中国三级夫妇交换| 特大巨黑吊av在线直播| 日韩av在线免费看完整版不卡| 少妇猛男粗大的猛烈进出视频| 大又大粗又爽又黄少妇毛片口| 成人黄色视频免费在线看| 欧美 亚洲 国产 日韩一| 欧美+日韩+精品| 日韩中文字幕视频在线看片| 黑人巨大精品欧美一区二区蜜桃 | 大片电影免费在线观看免费| 伊人久久国产一区二区| av免费观看日本| 亚洲精品日本国产第一区| 人妻制服诱惑在线中文字幕| 国产高清国产精品国产三级| 精品久久久噜噜| 欧美日韩视频精品一区| 久久人人爽人人片av| 99国产综合亚洲精品| 青春草亚洲视频在线观看| 国产深夜福利视频在线观看| 日韩人妻高清精品专区| 久久午夜福利片| 在线播放无遮挡| 各种免费的搞黄视频| 另类精品久久| 久久热精品热| 日韩亚洲欧美综合| 亚洲经典国产精华液单| 黄色怎么调成土黄色| 精品久久久噜噜| 看非洲黑人一级黄片| 最近2019中文字幕mv第一页| 天堂中文最新版在线下载| 啦啦啦中文免费视频观看日本| 免费人成在线观看视频色| 国产极品天堂在线| 校园人妻丝袜中文字幕| 国产精品久久久久久精品古装| 99re6热这里在线精品视频| 亚洲怡红院男人天堂| 国产成人av激情在线播放 | 伊人久久精品亚洲午夜| 国产熟女午夜一区二区三区 | 91成人精品电影| 免费av不卡在线播放| av福利片在线| 色视频在线一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲国产精品999| 欧美日韩在线观看h| 美女cb高潮喷水在线观看| 久久毛片免费看一区二区三区| 99热国产这里只有精品6| 国产亚洲精品久久久com| 免费人成在线观看视频色| 22中文网久久字幕| 国产精品一区二区在线不卡| 欧美 日韩 精品 国产| 人人妻人人澡人人爽人人夜夜| 亚洲欧洲国产日韩| 嫩草影院入口| 日本91视频免费播放| 国产av精品麻豆| 日本vs欧美在线观看视频| 我要看黄色一级片免费的| 亚洲av成人精品一区久久| 人妻夜夜爽99麻豆av| 97在线视频观看| 丝袜脚勾引网站| kizo精华| 69精品国产乱码久久久| 在线观看免费高清a一片| 久久久国产欧美日韩av| 成人毛片a级毛片在线播放| 伦理电影免费视频| 欧美亚洲日本最大视频资源| 亚洲av二区三区四区| 亚洲综合色网址| 国产乱来视频区| 久久女婷五月综合色啪小说| 少妇的逼好多水| 蜜桃国产av成人99| 免费人妻精品一区二区三区视频| 精品一区二区三区视频在线| 亚洲av日韩在线播放| 国产精品久久久久久av不卡| 女性生殖器流出的白浆| 国产成人免费无遮挡视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费大片黄手机在线观看| 成人影院久久| 久久久欧美国产精品| xxx大片免费视频| 黑人巨大精品欧美一区二区蜜桃 | 99热全是精品| 这个男人来自地球电影免费观看 | 男人操女人黄网站| 另类精品久久| 亚洲美女黄色视频免费看| 97精品久久久久久久久久精品| 欧美三级亚洲精品| 日韩免费高清中文字幕av| 少妇高潮的动态图| 一区二区三区四区激情视频| 看免费成人av毛片| 国产综合精华液| 久久久欧美国产精品| 国产国拍精品亚洲av在线观看| 久久久国产欧美日韩av| 亚洲av男天堂| 久久久国产欧美日韩av| 免费不卡的大黄色大毛片视频在线观看| 久久久久久人妻| 69精品国产乱码久久久| 亚洲国产成人一精品久久久| 日韩伦理黄色片| 美女视频免费永久观看网站| 观看av在线不卡| 妹子高潮喷水视频| 91aial.com中文字幕在线观看| 午夜激情福利司机影院| 国产极品天堂在线| 国产亚洲一区二区精品| 成人手机av| 三级国产精品片| 一级,二级,三级黄色视频| 久久久久精品性色| 中文天堂在线官网| 狂野欧美激情性bbbbbb| 男女啪啪激烈高潮av片| 免费黄色在线免费观看| 国产极品天堂在线| 久久久精品区二区三区| 国产 一区精品| 国产69精品久久久久777片| 亚洲精品自拍成人| 国产精品一二三区在线看| 最新的欧美精品一区二区| 少妇人妻精品综合一区二区| 视频在线观看一区二区三区| 国产精品一二三区在线看| freevideosex欧美| 妹子高潮喷水视频| 最黄视频免费看| 亚洲成人一二三区av| 国产精品偷伦视频观看了| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看 | 各种免费的搞黄视频| 久热久热在线精品观看| 欧美精品一区二区大全| 亚洲精品久久午夜乱码| tube8黄色片| 超色免费av| 亚洲国产最新在线播放| 一级毛片aaaaaa免费看小| 精品国产乱码久久久久久小说| 亚洲国产精品成人久久小说| 99国产精品免费福利视频| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 中文字幕人妻丝袜制服| 久久精品久久精品一区二区三区| av在线老鸭窝| 欧美人与性动交α欧美精品济南到 | 中文乱码字字幕精品一区二区三区| 美女脱内裤让男人舔精品视频| 国产免费视频播放在线视频| 97超碰精品成人国产| 最近最新中文字幕免费大全7| 99久久人妻综合| 精品久久久久久电影网| 久久久久久久亚洲中文字幕| 国产亚洲欧美精品永久| 成人午夜精彩视频在线观看| 一级二级三级毛片免费看| 日日摸夜夜添夜夜添av毛片| 啦啦啦中文免费视频观看日本| www.av在线官网国产| 国产免费视频播放在线视频| 插逼视频在线观看| av不卡在线播放| 国产成人精品久久久久久| 成年人午夜在线观看视频| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 亚洲精品亚洲一区二区| 超碰97精品在线观看| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| 视频区图区小说| 亚洲美女搞黄在线观看| 各种免费的搞黄视频| 成人综合一区亚洲| 欧美少妇被猛烈插入视频| 亚洲欧美成人精品一区二区| 欧美日韩视频高清一区二区三区二| 十八禁网站网址无遮挡| 精品卡一卡二卡四卡免费| 九九在线视频观看精品| 亚洲欧美中文字幕日韩二区| 男女边摸边吃奶| 国产乱来视频区| 男人操女人黄网站| 亚洲精品一区蜜桃| 自拍欧美九色日韩亚洲蝌蚪91| 日韩一区二区视频免费看| 亚洲精品国产色婷婷电影| 免费看不卡的av| 黄片播放在线免费| 亚洲国产精品999| 一级片'在线观看视频| 丝袜喷水一区| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 草草在线视频免费看| 日本91视频免费播放| 国产成人精品婷婷| kizo精华| 国产免费视频播放在线视频| 中文字幕人妻丝袜制服| 最近的中文字幕免费完整| 国产极品天堂在线| 18禁观看日本| 成人国语在线视频| 考比视频在线观看| 啦啦啦中文免费视频观看日本| 色94色欧美一区二区| 一本大道久久a久久精品| 中文天堂在线官网| 国产精品久久久久久精品古装| 亚洲美女黄色视频免费看| 久久久久国产精品人妻一区二区| 国产精品久久久久久久电影| 欧美日韩视频高清一区二区三区二| 免费人成在线观看视频色| 一区二区av电影网| 亚洲五月色婷婷综合| 只有这里有精品99| 三上悠亚av全集在线观看| 桃花免费在线播放| av视频免费观看在线观看| 视频中文字幕在线观看| 最近中文字幕2019免费版| 美女国产高潮福利片在线看| av在线观看视频网站免费| 在线亚洲精品国产二区图片欧美 | h视频一区二区三区| 久久综合国产亚洲精品| 色94色欧美一区二区| 日本欧美国产在线视频| 在线观看免费日韩欧美大片 | 2018国产大陆天天弄谢| 亚洲精华国产精华液的使用体验| 精品熟女少妇av免费看| 亚洲国产精品成人久久小说| 国产精品一国产av| 久久精品国产鲁丝片午夜精品| 又粗又硬又长又爽又黄的视频| 18+在线观看网站| 国精品久久久久久国模美| 在线天堂最新版资源| 在线观看美女被高潮喷水网站| 黄色一级大片看看| 2021少妇久久久久久久久久久| 亚洲精品自拍成人| 99热这里只有精品一区| 亚洲av福利一区| 国产精品不卡视频一区二区| 久久韩国三级中文字幕| 欧美精品一区二区大全| 欧美成人精品欧美一级黄| 十分钟在线观看高清视频www| 午夜视频国产福利| 久久国产亚洲av麻豆专区| 亚洲欧美精品自产自拍| 国产69精品久久久久777片| 国产一区二区在线观看av| 久久精品熟女亚洲av麻豆精品| 青青草视频在线视频观看| 久久久久久久久久人人人人人人| 曰老女人黄片| 婷婷色综合大香蕉| 久久国产精品男人的天堂亚洲 | 啦啦啦啦在线视频资源| 桃花免费在线播放| 亚洲欧美成人精品一区二区| 好男人视频免费观看在线| 天天操日日干夜夜撸| 999精品在线视频| 十分钟在线观看高清视频www| 国产精品国产av在线观看| 免费观看无遮挡的男女| 纵有疾风起免费观看全集完整版| 国产精品久久久久久久久免| av又黄又爽大尺度在线免费看| 日日摸夜夜添夜夜添av毛片| 能在线免费看毛片的网站| 欧美精品人与动牲交sv欧美| 三上悠亚av全集在线观看| 在线观看三级黄色| 精品人妻熟女毛片av久久网站| 美女福利国产在线| 精品人妻偷拍中文字幕| 国产欧美日韩一区二区三区在线 | 亚洲,欧美,日韩| 久久午夜福利片| 男女边摸边吃奶| 蜜桃久久精品国产亚洲av| 国产精品免费大片| 女人久久www免费人成看片| 老司机亚洲免费影院| 国产精品成人在线| 国产亚洲精品第一综合不卡 | 欧美3d第一页| 爱豆传媒免费全集在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人手机av| 乱码一卡2卡4卡精品| 亚洲av电影在线观看一区二区三区| 亚洲经典国产精华液单| 少妇 在线观看| 精品久久久久久电影网| 国产精品人妻久久久久久| 一区二区三区免费毛片| 婷婷色麻豆天堂久久| 国产乱人偷精品视频| 成人影院久久| 在线天堂最新版资源| 性色av一级| 日韩三级伦理在线观看| 伊人久久精品亚洲午夜| a级毛片黄视频| 99热国产这里只有精品6| 在线观看免费视频网站a站| 一二三四中文在线观看免费高清| 国产永久视频网站| 久久精品国产a三级三级三级| 精品一区在线观看国产| 99精国产麻豆久久婷婷| 99热这里只有是精品在线观看| 久久久久久久久大av| 尾随美女入室| 久久国内精品自在自线图片| 日本午夜av视频| 久久婷婷青草| 美女内射精品一级片tv| 欧美亚洲 丝袜 人妻 在线| 日本91视频免费播放| 日韩成人av中文字幕在线观看| 免费观看av网站的网址| 欧美成人精品欧美一级黄| 2018国产大陆天天弄谢| 久久久国产精品麻豆| 一级黄片播放器| 十分钟在线观看高清视频www| 高清视频免费观看一区二区| 色婷婷久久久亚洲欧美| 国产极品粉嫩免费观看在线 | 在线观看三级黄色| 岛国毛片在线播放| 97精品久久久久久久久久精品| 久久久久久久久大av| 麻豆成人av视频| 国产老妇伦熟女老妇高清| 日韩精品有码人妻一区| 国产精品人妻久久久久久| 婷婷色综合www| 亚洲精品日韩av片在线观看| 精品一品国产午夜福利视频| 男男h啪啪无遮挡| 亚洲精品视频女| 高清毛片免费看| 国产日韩一区二区三区精品不卡 | 大片免费播放器 马上看| 国产精品国产三级专区第一集| 国产精品麻豆人妻色哟哟久久| 日本av手机在线免费观看| 伊人亚洲综合成人网| 十八禁高潮呻吟视频| 人人妻人人添人人爽欧美一区卜| 交换朋友夫妻互换小说| 日日啪夜夜爽| 观看美女的网站| 精品人妻偷拍中文字幕| 涩涩av久久男人的天堂| 亚洲精品久久午夜乱码| 日韩人妻高清精品专区| 日本91视频免费播放| av在线播放精品| 国产黄色免费在线视频| 一级爰片在线观看| 久久韩国三级中文字幕| 国产精品偷伦视频观看了| 久久狼人影院| 国产免费一区二区三区四区乱码| 日本av免费视频播放| 精品一品国产午夜福利视频| 一级毛片我不卡| 免费观看a级毛片全部| 少妇的逼水好多| 亚洲欧美清纯卡通| 国产极品粉嫩免费观看在线 | 亚洲欧洲日产国产| 最近中文字幕2019免费版| 在线看a的网站| 99久久人妻综合| 久久久久久久精品精品| 午夜日本视频在线| 黄色欧美视频在线观看|