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

    SCR脫硝催化反應(yīng)模型的建立與驗(yàn)證

    2016-10-18 08:50:54黃慶華宋麗云
    關(guān)鍵詞:模擬計(jì)算反應(yīng)器煙氣

    黃慶華,宋麗云,吳 銳,何 洪,3

    (1.北京工業(yè)大學(xué)區(qū)域大氣復(fù)合污染防治北京市重點(diǎn)實(shí)驗(yàn)室,北京 100124;2.北京工業(yè)大學(xué)綠色催化與分離北京市重點(diǎn)實(shí)驗(yàn)室,北京 100124;3.北京電動(dòng)車輛協(xié)同創(chuàng)新中心,北京 100081)

    SCR脫硝催化反應(yīng)模型的建立與驗(yàn)證

    黃慶華1,2,宋麗云1,2,吳 銳1,2,何 洪1,2,3

    (1.北京工業(yè)大學(xué)區(qū)域大氣復(fù)合污染防治北京市重點(diǎn)實(shí)驗(yàn)室,北京 100124;2.北京工業(yè)大學(xué)綠色催化與分離北京市重點(diǎn)實(shí)驗(yàn)室,北京 100124;3.北京電動(dòng)車輛協(xié)同創(chuàng)新中心,北京 100081)

    為研究硫酸氫銨在低溫選擇性催化還原(selective catalytic reduction,SCR)催化反應(yīng)器中的生成條件、沉積及富集規(guī)律和催化劑再生技術(shù),設(shè)計(jì)和搭建參照實(shí)際低溫脫硝工程的中試規(guī)模的SCR反應(yīng)器系統(tǒng).在實(shí)驗(yàn)過(guò)程中,不能準(zhǔn)確測(cè)量出沿催化劑軸向或在反應(yīng)截面上NH3或NOx質(zhì)量濃度的分布情況,而該情況對(duì)實(shí)驗(yàn)研究是比較重要的.為此,采用化學(xué)反應(yīng)動(dòng)力學(xué)和計(jì)算流體力學(xué)建立了SCR脫硝催化反應(yīng)模擬模型,通過(guò)在相同反應(yīng)條件下分別設(shè)定不同氨氮比、溫度,比對(duì)實(shí)驗(yàn)測(cè)量數(shù)據(jù)間的變化趨勢(shì)與模擬計(jì)算結(jié)果的一致性,驗(yàn)證所建催化反應(yīng)模擬模型的正確性.該模型可為分析研究硫酸氫銨在反應(yīng)截面上或沿催化劑軸向的生成條件、沉積及富集規(guī)律提供數(shù)據(jù)支撐,同時(shí)也能用于實(shí)際SCR脫硝工程脫硝催化反應(yīng)過(guò)程的模擬.

    SCR脫硝;脫硝工程;氨逃逸;催化反應(yīng);反應(yīng)模擬

    在城市大氣污染中,NOx是主要污染物之一.據(jù)相關(guān)文獻(xiàn)報(bào)道,我國(guó)城市NOx排放主要來(lái)源于燃油(氣)交通運(yùn)輸工具、火力發(fā)電鍋爐和工業(yè)鍋(窯)爐[1].脫硝(DeNOx)主流技術(shù)NOx選擇性催化還原(selective catalytic reduction,SCR)技術(shù)包括中溫SCR脫硝工藝(催化劑工作溫度在300~400℃)和低溫SCR脫硝工藝(催化劑工作溫度低于300℃)[2].中溫SCR脫硝工藝在火力發(fā)電鍋爐脫硝領(lǐng)域已得到廣泛應(yīng)用.我國(guó)許多工業(yè)鍋(窯)爐(例如工業(yè)鍋爐、玻璃陶瓷爐窯、水泥爐窯、冶金燒結(jié)爐、煉焦和石化系統(tǒng)的裂解設(shè)備等)的脫硝,因其煙氣排放溫度較低(大多在200℃左右),比較適合低溫SCR脫硝工藝.研究和開(kāi)發(fā)低溫SCR催化劑已成為國(guó)內(nèi)外學(xué)術(shù)和工業(yè)界的熱點(diǎn)問(wèn)題.目前,已研究出工作溫度低至160℃的高效低溫SCR釩鈦基催化劑,但在工程應(yīng)用時(shí)遇到了必須克服的障礙,即在低溫條件下催化劑的SO2中毒和因硫酸氫銨(ammonium bisulfate,ABS)在催化劑的冷凝附著而造成的催化劑失活[3].

    為研究ABS在低溫SCR催化反應(yīng)器中的生成條件、沉積及富集規(guī)律和催化劑再生技術(shù),設(shè)計(jì)和搭建參照實(shí)際低溫脫硝工程的中試規(guī)模的SCR反應(yīng)器系統(tǒng).在實(shí)際實(shí)驗(yàn)過(guò)程中,由于反應(yīng)截面較小等客觀原因,不能準(zhǔn)確測(cè)量出在反應(yīng)截面上或沿催化劑軸向NH3和NOx的質(zhì)量濃度分布情況,而該情況對(duì)實(shí)驗(yàn)研究是比較重要的,為此作者建立了SCR脫硝催化反應(yīng)數(shù)值模擬模型.

    數(shù)值模擬在煙氣脫硝中的應(yīng)用,不但包括對(duì)脫硝工藝整體系統(tǒng)布置和流動(dòng)特性的模擬研究[4-6],還包括對(duì)SCR催化反應(yīng)過(guò)程的數(shù)值模擬研究.在SCR催化反應(yīng)過(guò)程中,各種組分在催化劑表面上的反應(yīng)是核心,采用數(shù)學(xué)模型可用于指導(dǎo)SCR催化劑的優(yōu)化設(shè)計(jì).

    Beeckman等[7]建立了SCR催化劑單孔道的一維模型,分析研究催化劑孔結(jié)構(gòu)對(duì)反應(yīng)活性的影響.沈伯雄等[8]建立了SCR催化劑單孔道的一維模型,模擬SCR催化劑孔道內(nèi)的催化反應(yīng)進(jìn)程.

    在SCR催化劑一維模型研究的基礎(chǔ)上,Dhanushkodi等[9]建立了SCR催化劑的二維模型,將催化劑孔道假設(shè)為圓柱狀,使模型具有便于簡(jiǎn)化的二維幾何旋轉(zhuǎn)對(duì)稱性,模擬計(jì)算結(jié)果與其實(shí)驗(yàn)數(shù)據(jù)一致.

    在以上研究基礎(chǔ)上,作者建立了SCR脫硝系統(tǒng)三維數(shù)學(xué)模型,該模型可為分析研究ABS在反應(yīng)截面上或沿催化劑軸向的生成條件、沉積及富集規(guī)律提供數(shù)據(jù)支撐,可用于指導(dǎo)SCR催化劑的優(yōu)化設(shè)計(jì),同時(shí)也能用于實(shí)際SCR脫硝工程脫硝催化反應(yīng)過(guò)程的模擬.

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

    1.1中試規(guī)模的SCR反應(yīng)器系統(tǒng)

    該系統(tǒng)的建立參照了實(shí)際低溫脫硝工程,包括低溫SCR反應(yīng)器、燃燒器、鼓風(fēng)機(jī)、NO反應(yīng)器、配氣系統(tǒng)和控制系統(tǒng),相當(dāng)于化學(xué)反應(yīng)工程中的單管反應(yīng)器.NO由NH3在空氣中氧化得到,通過(guò)調(diào)整NH3的量來(lái)保證NO在煙氣中的質(zhì)量濃度.

    該系統(tǒng)使用2塊150 mm×150 mm×800 mm催化劑模塊,如圖1所示.催化劑上下順序安裝,形成雙床層催化劑單體的反應(yīng)器結(jié)構(gòu).分別命名SCR脫硝反應(yīng)器系統(tǒng)中測(cè)溫點(diǎn)、取樣點(diǎn)(上)為SITE1,測(cè)溫點(diǎn)、取樣點(diǎn)(中)為SITE2,測(cè)溫點(diǎn)、取樣點(diǎn)(下)為SITE3.若以反應(yīng)器軸線為橫坐標(biāo)L,SITE1為橫坐標(biāo)原點(diǎn),正方向?yàn)闊煔饬鲃?dòng)方向,則SITE2處橫坐標(biāo)為1.2 m,SITE3處橫坐標(biāo)為2.4 m.

    1.2研究方法

    采用化學(xué)反應(yīng)動(dòng)力學(xué)、計(jì)算流體動(dòng)力學(xué)和系統(tǒng)實(shí)驗(yàn)的方法,建立SCR脫硝催化反應(yīng)模型.通過(guò)模擬計(jì)算得到在反應(yīng)截面中心點(diǎn)處的NOx質(zhì)量濃度,然后比對(duì)反應(yīng)截面中心點(diǎn)(即SCR脫硝反應(yīng)器系統(tǒng)中3個(gè)取樣點(diǎn))NOx質(zhì)量濃度模擬數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)的變化趨勢(shì)是否一致來(lái)驗(yàn)證模擬計(jì)算結(jié)果是否準(zhǔn)確,進(jìn)而驗(yàn)證SCR脫硝催化反應(yīng)模型的正確性,確定模型參數(shù),再通過(guò)模擬計(jì)算得到在反應(yīng)截面上NOx質(zhì)量濃度的分布情況.

    1.3催化反應(yīng)模擬計(jì)算模型的建立

    1.3.1理論推導(dǎo)

    1)催化反應(yīng)機(jī)理

    Topsoe等[10-11]、Maki等[12]采用光譜技術(shù)研究表面催化反應(yīng),提出了以下關(guān)于SCR催化反應(yīng)的結(jié)論:①對(duì)氨氣吸附的量在催化反應(yīng)條件下最多;②SCR催化反應(yīng)的活性位位于Bronsted酸位;③SCR催化反應(yīng)中的氨氣分子在Bronsted酸位上的吸附與脫附存在一個(gè)平衡過(guò)程;④被吸附的氨氣分子只與距其最近被吸附在表面釩元素的NO分子反應(yīng);⑤在催化反應(yīng)條件下不會(huì)發(fā)生催化劑表面對(duì)NO分子的大量吸附;⑥選擇性催化反應(yīng)是NO分子與活化后的氨氣分子發(fā)生氣態(tài)反應(yīng)生成氮?dú)夂退瑫r(shí)使催化劑部分被還原的反應(yīng)過(guò)程[13-16].他們認(rèn)為在SCR催化劑表面主要發(fā)生的反應(yīng)過(guò)程見(jiàn)圖2.

    基于以上催化反應(yīng)機(jī)理,建立的模擬計(jì)算模型忽略氣相反應(yīng),僅考慮在催化劑表面上進(jìn)行的催化反應(yīng),主要反應(yīng)公式[17]如下:

    2)催化反應(yīng)速率方程

    揭示反應(yīng)速率常數(shù)對(duì)溫度依賴關(guān)系的Arrhenius定理有3種數(shù)學(xué)表達(dá)式[18],如下:

    式中:k為T(mén)(K)時(shí)的反應(yīng)速率常數(shù);A為指數(shù)前因子;E為活化能.

    因Arrhenius定理由于假設(shè)指數(shù)因子A與溫度T無(wú)關(guān)是不精確的,實(shí)際上A與反應(yīng)溫度T的某次方成比例,即A=A0Tm.其中:A0是與溫度無(wú)關(guān)的常數(shù);|m|是小于等于4的整數(shù)或半整數(shù);A=A0Tm.

    所以,催化反應(yīng)速率方程為

    式中:A0為指前因子;m為溫度指數(shù);E為活化能.

    3)熱力學(xué)方程

    應(yīng)用量子力學(xué)原理,根據(jù)光譜分析的數(shù)據(jù),可以得到理想氣體比熱容與溫度的關(guān)系式.通常比定壓熱容、定壓焓和定壓熵可表示為溫度的多項(xiàng)式[19]

    4)流體動(dòng)力學(xué)方程

    煙氣流動(dòng)遵循物理守恒定律[20],這些守恒定律包括質(zhì)量守恒定律、動(dòng)量守恒定律、能量守恒定律.這些守恒定律通過(guò)相應(yīng)的控制方程來(lái)進(jìn)行數(shù)學(xué)描述,這些控制的通用形式為

    1.3.2模擬計(jì)算的假設(shè)條件

    為進(jìn)行催化反應(yīng)動(dòng)力學(xué)和計(jì)算流體動(dòng)力學(xué)模擬計(jì)算,建立催化反應(yīng)模擬計(jì)算模型時(shí),首先對(duì)系統(tǒng)進(jìn)行一些簡(jiǎn)化或假設(shè):1)忽略灰塵的影響;2)假設(shè)煙氣為理想氣體;3)假設(shè)煙氣的組成只有N2、CO2、O2和H2O;4)不考慮氣相化學(xué)反應(yīng)及影響;5)活性位均勻分布于催化劑表面;6)流動(dòng)是定常的.

    1.3.3三維建模

    采用ANSYS 14.0 ICEM軟件建模工具進(jìn)行三維建模.催化劑層按多孔介質(zhì)模擬.煙道、SCR反應(yīng)器都按照實(shí)際尺寸和圖紙建模.系統(tǒng)初次計(jì)算時(shí)網(wǎng)格數(shù)為31萬(wàn),驗(yàn)證網(wǎng)格數(shù)影響時(shí)網(wǎng)格數(shù)為56萬(wàn).

    1.3.4實(shí)際實(shí)驗(yàn)和模擬計(jì)算的條件

    模擬計(jì)算和實(shí)際實(shí)驗(yàn)采用相同的煙氣成分、流量、NO質(zhì)量濃度和反應(yīng)溫度.

    煙氣成分中,O2的體積分?jǐn)?shù)為18%,CO2的為3%,H2O的為8%,N2的為71%;純NH3;NO質(zhì)量濃度為1 g/m3;煙氣流量為140 m3/h;空塔速度為4 000/h.

    1.3.5模擬計(jì)算時(shí)的參數(shù)選擇

    采用ANSYS 14.0 fluent軟件進(jìn)行計(jì)算.三維湍流數(shù)值模擬方法采用標(biāo)準(zhǔn)κ-ε模型,模型常數(shù)分別取C1ε=1.44,C2ε=1.92,Cμ=0.09,σκ=1.0,σε= 1.3.采用SIMPLE壓力-速度耦合算法.收斂判斷標(biāo)準(zhǔn)為速度殘值絕對(duì)標(biāo)準(zhǔn)0.001,能量計(jì)算殘值絕對(duì)標(biāo)準(zhǔn)為10-6,NO、O2、N2、H2O、CO2和NH3的濃度計(jì)算殘值絕對(duì)標(biāo)準(zhǔn)都為0.000 1.

    1.3.6模擬計(jì)算時(shí)的操作方法

    通過(guò)以下步驟或方法實(shí)現(xiàn)SCR脫硝催化反應(yīng)過(guò)程的數(shù)值模擬:首先建立三維幾何數(shù)字模型、反應(yīng)機(jī)理文件和反應(yīng)組分的物性參數(shù)文件,然后在ANSYS 14.0 fluent中分別導(dǎo)入以上文件,接著進(jìn)行計(jì)算參數(shù)和邊界條件設(shè)置,進(jìn)行初步計(jì)算,最后選取初步計(jì)算的重要參數(shù)值與實(shí)際實(shí)驗(yàn)相應(yīng)參數(shù)值的差別進(jìn)行分析,不斷優(yōu)化反應(yīng)機(jī)理文件,使選取的計(jì)算參數(shù)值與實(shí)際實(shí)驗(yàn)相應(yīng)參數(shù)值變化趨勢(shì)一致.其中,建立反應(yīng)機(jī)理文件和反應(yīng)組分的熱力學(xué)參數(shù)文件為操作方法的關(guān)鍵.

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

    2.1催化反應(yīng)模擬模型的驗(yàn)證

    以反應(yīng)器軸線為橫坐標(biāo),正方向?yàn)闊煔饬鲃?dòng)方向,縱坐標(biāo)為煙氣中NO質(zhì)量濃度,SITE1代表第1層催化劑上中點(diǎn)位置,橫坐標(biāo)為0 mm;SITE2代表第2層催化劑中間平面中點(diǎn)位置,橫坐標(biāo)為1.2 m;SITE3代表第2層催化劑下中點(diǎn)位置,橫坐標(biāo)為2.4 m.F表示實(shí)際實(shí)驗(yàn)數(shù)據(jù);M表示模擬數(shù)據(jù).煙氣從SITE1通過(guò)2層催化劑向SITE3流動(dòng)過(guò)程中,通過(guò)比對(duì)NO質(zhì)量濃度分別在實(shí)際實(shí)驗(yàn)和模擬計(jì)算過(guò)程中的變化規(guī)律是否一致,驗(yàn)證SCR催化反應(yīng)模擬模型的正確性.

    2.1.1不同氨氮比時(shí)模擬結(jié)果與實(shí)驗(yàn)結(jié)果的比對(duì)

    在反應(yīng)溫度為170℃,氨氮比分別為0.5、0.7、0.9、1.1、1.3時(shí),SITE1、SITE2、SITE3三位置處NO質(zhì)量濃度的實(shí)際實(shí)驗(yàn)值和模擬計(jì)算值的關(guān)系見(jiàn)圖3~7.

    從圖3~7可以看出,在反應(yīng)溫度為170℃,氨氮比分別為0.5、0.7、0.9、1.1、1.3時(shí),模擬結(jié)果與實(shí)驗(yàn)結(jié)果的變化趨勢(shì)均一致.

    2.1.2不同溫度時(shí)模擬結(jié)果與實(shí)驗(yàn)結(jié)果的比對(duì)

    在氨氮比為1∶1.1,反應(yīng)溫度分別為170、180、200、220、250℃時(shí),SITE1、SITE2、SITE3三位置處NO質(zhì)量濃度的實(shí)際實(shí)驗(yàn)值和模擬計(jì)算值的關(guān)系見(jiàn)圖8~12.

    從圖8~12可以看出,在氨氮比為1∶1.1,反應(yīng)溫度分別為170、180、200、220、250℃時(shí),模擬結(jié)果與實(shí)驗(yàn)結(jié)果的變化趨勢(shì)均一致.

    2.1.3驗(yàn)證結(jié)果

    在不同氨氮比、溫度的情況下,實(shí)際實(shí)驗(yàn)和模擬計(jì)算的結(jié)果變化趨勢(shì)均一致,從不同角度說(shuō)明,模擬計(jì)算基本可以反映實(shí)驗(yàn)結(jié)果的變化趨勢(shì),驗(yàn)證了本文建立的SCR催化反應(yīng)的模擬計(jì)算模型是正確的.

    2.2NO質(zhì)量濃度和NH3質(zhì)量濃度分布情況

    NO質(zhì)量濃度分布和NH3質(zhì)量濃度分布情況(如圖13~16所示)的實(shí)驗(yàn)條件如下:

    煙氣成分中,O2的體積分?jǐn)?shù)為18%,CO2的為3%,H2O的為8%,N2的為71%;純NH3;NO質(zhì)量濃度為1 g/m3;氨氮比為1∶1;反應(yīng)溫度為170℃;煙氣流量為140 m3/h;空塔速度為4 000/h.

    2.3討論

    2.3.1計(jì)算結(jié)果網(wǎng)格關(guān)聯(lián)性驗(yàn)證

    因網(wǎng)格數(shù)的設(shè)置可能會(huì)影響計(jì)算結(jié)果,但在一定誤差范圍內(nèi),當(dāng)設(shè)置的網(wǎng)格數(shù)大于某個(gè)值后,網(wǎng)格數(shù)的變化對(duì)計(jì)算結(jié)果就不存在影響了.所以,需要對(duì)初次模擬計(jì)算時(shí)設(shè)置的網(wǎng)格數(shù)是否會(huì)對(duì)計(jì)算結(jié)果產(chǎn)生影響進(jìn)行驗(yàn)證.以上計(jì)算結(jié)果采用的均是驗(yàn)證后的計(jì)算結(jié)果.

    2.3.2實(shí)際實(shí)驗(yàn)引起誤差的因素

    在實(shí)際實(shí)驗(yàn)過(guò)程中,整個(gè)實(shí)驗(yàn)系統(tǒng)中對(duì)數(shù)據(jù)的測(cè)量可能存在較大誤差,主要原因?yàn)?1)煙氣管道和反應(yīng)器中存在湍流;2)燃燒不穩(wěn)定,煙氣流量變化比較大;3)NO是通過(guò)NH3與空氣中的O2在NO生成器中反應(yīng)產(chǎn)生的,NO質(zhì)量濃度不穩(wěn)定.

    2.3.3模擬模型正確性的驗(yàn)證

    由于實(shí)際實(shí)驗(yàn)系統(tǒng)中數(shù)據(jù)測(cè)量可能存在較大誤差,若采用模擬計(jì)算數(shù)值與實(shí)際測(cè)量數(shù)值進(jìn)行比對(duì),結(jié)果可能不準(zhǔn)確,但實(shí)際實(shí)驗(yàn)過(guò)程中的數(shù)據(jù)變化趨勢(shì)與模擬過(guò)程中的數(shù)據(jù)變化趨勢(shì)的比對(duì)可以弱化實(shí)際實(shí)驗(yàn)中數(shù)據(jù)測(cè)量不準(zhǔn)確的影響,若變化趨勢(shì)一致就能確定模擬計(jì)算模型的正確性.

    例如,在反應(yīng)溫度為170℃,氨氮比分別為0.5、0.7、0.9、1.1、1.3時(shí),SITE1、SITE2、SITE3三位置處NO質(zhì)量濃度的實(shí)際實(shí)驗(yàn)值和模擬計(jì)算值的關(guān)系見(jiàn)圖3~7,總的趨勢(shì)是SITE1和SITE3處模擬計(jì)算值略低于實(shí)際實(shí)驗(yàn)值,SITE2處模擬計(jì)算值略高于實(shí)際實(shí)驗(yàn)值,但圖3例外,圖中SITE2處的情況可能因?qū)嶋H實(shí)驗(yàn)測(cè)量不準(zhǔn)確引起的.但因?yàn)閷?shí)際實(shí)驗(yàn)值和模擬計(jì)算值總的變化趨勢(shì)是一致的,可以確定模擬計(jì)算模型在不同氨氮比的條件下是正確的,是可以適用的.

    3 結(jié)論

    1)實(shí)際實(shí)驗(yàn)各測(cè)量位間數(shù)值的變化趨勢(shì)與模擬計(jì)算的數(shù)值變化趨勢(shì)一致,說(shuō)明SCR脫硝催化反應(yīng)模擬模型是正確的,可以應(yīng)用于實(shí)際工程.

    2)模擬計(jì)算所得到的軸向剖面和各反應(yīng)截面NO質(zhì)量濃度分布情況可以為分析中試實(shí)驗(yàn)結(jié)果提供數(shù)據(jù)支撐.

    3)煙氣成分中,O2的體積分?jǐn)?shù)為18%,CO2的為3%,H2O的為8%,N2的為71%.在反應(yīng)溫度為170℃,NO質(zhì)量濃度為1 g/m3,純NH3,氨氮比為1∶1,煙氣流量為140 m3/h,空塔速度為4 000/h的條件下,SITE1處NO質(zhì)量濃度實(shí)驗(yàn)數(shù)據(jù)為1 g/m3,模擬數(shù)據(jù)為1 g/m3;SITE2處NO質(zhì)量濃度實(shí)驗(yàn)數(shù)據(jù)為210 mg/m3,模擬數(shù)據(jù)為218 mg/m3,基本一致;SITE3處NO質(zhì)量濃度實(shí)驗(yàn)數(shù)據(jù)為50 mg/m3,模擬數(shù)據(jù)為62 mg/m3,基本一致.所以在該條件下的NO質(zhì)量濃度分布模擬圖和軸向剖面NH3質(zhì)量濃度分布模擬圖可以用于分析中試實(shí)驗(yàn)結(jié)果.

    [1]郝吉明,馬大廣,王書(shū)肖.大氣污染控制工程[M].3版.北京:高等教育出版社,2011:378-379.

    [2]張曉鵬.基于Mn/Ce-ZrO2催化劑的低溫NH3-SCR脫硝性能研究[D].天津:南開(kāi)大學(xué),2013:6-12.ZHANG X P.Research of Mn/Ce-ZrO2for selective catalytic reduction of NOxwith NH3at low temperature[D].Tianjing:Nankai University,2013:6-12.(in Chinese)

    [3]何洪,王洪兵,李堅(jiān),等.水泥爐窯選擇性(非)催化還原脫硝技術(shù)[J].四川水泥,2012,46(1):46-52.HE H,WANG H B,LI J,et al.Cement kiln(non)selective catalytic reduction denitration technology[J].Sichuan Cement,2012,46(1):46-52.(in Chinese)

    [4]CHO J M,CHOI J W,HONG S H,et al.Application of computationalfluiddynamicanalysisforimproving performanceofcommercialscaleselectivecatalytic reduction[J].Korean Journal of Chemical Engineering,2005,23(1):43-56.

    [5]李芳芳,吳少華,曹慶喜,等.湍流混合限制下選擇性非催化還原過(guò)程的數(shù)值模擬[J].中國(guó)電機(jī)工程學(xué)報(bào),2008,28(23):46-52.LI F F,WU S H,CAO Q X,et al.Numerical simulation of selective non-catalytic reduction process with turbulence mixing influence[J].Proceedings of the CSEE,2008,28(23):46-52.(in Chinese)

    [6]雷達(dá),金保升.電站選擇性催化還原系統(tǒng)速度場(chǎng)測(cè)量方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2009,29(35):89-95.LEI D,JIN B S.Velocity distribution measuring method used in a power station selective catalytic reduction system[J].Proceedings of the CSEE,2009,29(35):89-95.(in Chinese)

    [7]BEECKMAN J W,HEGEDUS L L.Design of monolith catalysts for power plant NOxemission control[J].Industrial&Engineering Chemistry Research,1991,30(5):969-978.

    [8]沈伯雄,趙寧,劉亭.煙氣脫硝選擇性催化還原催化劑反應(yīng)模擬研究[J].中國(guó)電機(jī)工程學(xué)報(bào),2011,31(8): 31-37.SHEN B X,ZHAO N,LIU T.Mathematical simulation of flue fas denitration based on selective catalytic reduction catalyst[J].Proceedings of the CSEE,2011,31(8):31-37.(in Chinese)

    [9]DHANUSHKODI S R,MAHINPEY N,WILSON M.Kinetic and 2D reactor modeling for simulation of the catalytic reduction of NOxin the monolith honeycomb reactor[J].Process Safety and Environmental Protection,2008,86(4):303-309.

    [10]TOPSOE N Y,DUMESIC J A,TOPSOE H.Vanadia/titania catalysts for selective catalytic reduction of nitric oxide by ammonia[J].Journal of Catalysis,1995,151(1):241-252.

    [11]TOPSOE N Y,TOPSOE H,DUMESIC J A.Vanadia/titania catalysts for selective catalytic reduction(SCR)of nitric oxide by ammonia[J].Journal of Catalysis,1995,151(1):226-240.

    [12]MAKI T,MITSUYAKI S,TAKAHAHU O,et al.The mechanism of the reaction between NOx,and NH3,on V2O5in the presence of oxygen[J].Journal of Catalysis, 1977,50(3):441-446.

    [13]JOHN M.Comparison and analysis of intrinsic kinetics and fectiveness factors for the catalytic reduction of NO with ammonia in the presence of oxygen[J].Industrial& Engineering Chemistry Research,1992,31(4):987-994.

    [14]GIANGUIDO R,LI Yi,GUIDO B,et al.Adsorption,activation,and oxidation of ammonia over SCR catalysts[J].Journal of Catalysis,1995,157(2):523-535.

    [15]LINTZ H G,TUREK T.Intrinsic kinetics of nitric oxide reduction by ammonia on a vanadia-titania catalyst[J].Applied Catalysis A:General,1992,85(1):13-25.

    [16]MONICA C,BASMA M,CHRISTIAN M.Modeling catalytic reduction of NO by ammonia over V2O5[J].Surface Science Reports,2004,55(6):169-236.

    [17]DUMESIC J A,TOPSOE N Y,TOPSOE H,et al.Kinetics of selective catalytic reduction of nitric oxide by ammonia over vanadia/titania[J].Journal of Catalysis,1996,163(163):409-417.

    [18]趙學(xué)莊.化學(xué)反應(yīng)動(dòng)力學(xué)原理(上冊(cè))[M].北京:高等教育出版社,1984:19-29.

    [19]廉樂(lè)明,李力能,吳家正,等.工程熱力學(xué)[M].4版.北京:中國(guó)建筑工業(yè)出版社,1999:22-27.

    [20]王福軍.計(jì)算流體動(dòng)力學(xué)分析[M].北京:清華大學(xué)出版社,2011:7-13.

    (責(zé)任編輯 梁 潔)

    Establishment and Validation of the Model for SCR DeNOxCatalytic Reaction

    HUANG Qinghua1,2,SONG Liyun1,2,WU Rui1,2,HE Hong1,2,3
    (1.Beijing Key Laboratory of Regional Air Pollution Control,Beijing University of Technology,Beijing 100124,China;2.Beijing Key Laboratory of Green Catalysis and Separation,Beijing University of Technology,Beijing 100124,China;3.Collaborative Innovation Center of Electric Vehicles in Beijing,Beijing 100081,China)

    To study the conditions generated of ammonium bisulfate in the reactor of selective catalytic reduction(SCR)catalytic at low temperature,rules of its deposition and enrichment,and its renewable technologies,a pilot-scale system of SCR reactor was designed and constructed based on the actual SCR project at low temperature.During the experiment,the distribution of NH3or NOxconcentration could not be accurately measured which was axially along the catalyst or on the cross section of reactor,and which was more important to the experimental study.So a simulative model was created to follow kinetics of chemical reactions and computational fluid dynamics,whose correctness was verified by comparing consistency of the changing trends from experimental data and simulations,by different ammonia-nitrogen ratio,different temperatures and different rate of SCR set under the same conditions of reaction.The model can provide support for data analysis and research on the conditions generated of ammonium bisulfate in the reactor of SCR catalytic at low temperature,rules of its deposition and enrichment,while itcould be used to simulate the process of catalytic reaction in the actual SCR project.

    SCR DeNOx;DeNOxengineering;ammonia slip;catalyzed reaction;reaction simulation

    U 461;TP 308

    A

    0254-0037(2016)10-1532-08

    10.11936/bjutxb2016030007

    2016-03-03

    國(guó)家自然科學(xué)基金資助項(xiàng)目(21277009);北京市教育委員會(huì)資助項(xiàng)目(TJSHG201310005003)

    黃慶華(1972—),男,博士研究生,主要從事大氣污染控制、脫硝工藝、SCR脫硝工程、裝置流場(chǎng)設(shè)計(jì)、CFD模擬計(jì)算及物理模型實(shí)驗(yàn)方面的研究,E-mail:tsinghua.h@emails.bjut.edu.cn

    猜你喜歡
    模擬計(jì)算反應(yīng)器煙氣
    固體吸附劑脫除煙氣中SOx/NOx的研究進(jìn)展
    化工管理(2022年13期)2022-12-02 09:21:52
    R1234ze PVTx熱物性模擬計(jì)算
    能源工程(2022年1期)2022-03-29 01:06:26
    聚丙烯環(huán)管反應(yīng)器升密操作與控制
    云南化工(2021年11期)2022-01-12 06:06:36
    基于參數(shù)自整定模糊PID的SCR煙氣脫硝控制
    EGSB反應(yīng)器處理阿維菌素廢水
    上旋流厭氧反應(yīng)器在造紙廢水處理中的應(yīng)用
    擠出發(fā)泡片材褶皺分析及模擬計(jì)算
    基于非分散紫外吸收法的便攜式煙氣分析儀在煙氣二氧化硫監(jiān)測(cè)中的應(yīng)用研究
    費(fèi)托合成微反應(yīng)器研究進(jìn)展
    實(shí)際發(fā)射工況下底排藥柱結(jié)構(gòu)完整性的模擬計(jì)算
    女人精品久久久久毛片| 日本精品一区二区三区蜜桃| 久热爱精品视频在线9| 欧洲精品卡2卡3卡4卡5卡区| 中国美女看黄片| 在线免费观看的www视频| 黄网站色视频无遮挡免费观看| 国产高清videossex| 视频在线观看一区二区三区| 免费在线观看完整版高清| 午夜视频精品福利| 免费一级毛片在线播放高清视频 | 少妇 在线观看| 精品一区二区三区av网在线观看| 一区在线观看完整版| 欧美成人性av电影在线观看| 亚洲国产精品sss在线观看| www日本在线高清视频| 色尼玛亚洲综合影院| 色精品久久人妻99蜜桃| 一二三四社区在线视频社区8| 亚洲伊人色综图| 欧美国产精品va在线观看不卡| 成人18禁在线播放| 精品欧美国产一区二区三| 精品久久久久久久毛片微露脸| 成人三级黄色视频| 国产精品98久久久久久宅男小说| 9色porny在线观看| 国产精品秋霞免费鲁丝片| 搡老妇女老女人老熟妇| 黄色女人牲交| 日韩中文字幕欧美一区二区| 黄色毛片三级朝国网站| 在线av久久热| 色播亚洲综合网| 国产欧美日韩精品亚洲av| 天天躁狠狠躁夜夜躁狠狠躁| 午夜激情av网站| 精品国产乱子伦一区二区三区| 性少妇av在线| 免费看十八禁软件| 成人免费观看视频高清| avwww免费| 老汉色∧v一级毛片| 搞女人的毛片| 法律面前人人平等表现在哪些方面| 给我免费播放毛片高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美一区二区三区黑人| 18禁国产床啪视频网站| 伊人久久大香线蕉亚洲五| 欧美性长视频在线观看| 日韩欧美国产在线观看| 国产极品粉嫩免费观看在线| 18禁美女被吸乳视频| 欧美黑人欧美精品刺激| 欧美av亚洲av综合av国产av| 国产精品久久久av美女十八| 欧美av亚洲av综合av国产av| 婷婷六月久久综合丁香| 国产精品精品国产色婷婷| 一本久久中文字幕| 两人在一起打扑克的视频| 我的亚洲天堂| 国产成人精品久久二区二区免费| 在线十欧美十亚洲十日本专区| 精品高清国产在线一区| 夜夜看夜夜爽夜夜摸| 精品一品国产午夜福利视频| 国产三级黄色录像| 国产精品乱码一区二三区的特点 | 精品国产乱码久久久久久男人| 亚洲精品中文字幕一二三四区| 在线免费观看的www视频| 国产不卡一卡二| 长腿黑丝高跟| 亚洲久久久国产精品| 一夜夜www| 在线观看午夜福利视频| aaaaa片日本免费| 一卡2卡三卡四卡精品乱码亚洲| 人妻丰满熟妇av一区二区三区| 欧美乱码精品一区二区三区| 老熟妇乱子伦视频在线观看| 天天躁夜夜躁狠狠躁躁| 在线观看免费视频网站a站| 成年版毛片免费区| 自线自在国产av| 黑人巨大精品欧美一区二区蜜桃| 在线观看一区二区三区| 1024香蕉在线观看| 成在线人永久免费视频| 又大又爽又粗| av片东京热男人的天堂| 欧美日韩亚洲国产一区二区在线观看| 欧美中文综合在线视频| 亚洲在线自拍视频| 少妇 在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 午夜a级毛片| 91精品三级在线观看| 视频在线观看一区二区三区| 色尼玛亚洲综合影院| 日日爽夜夜爽网站| 免费在线观看黄色视频的| 午夜福利在线观看吧| 免费观看人在逋| 国产午夜精品久久久久久| 国产亚洲欧美在线一区二区| 1024视频免费在线观看| 久久人妻熟女aⅴ| 变态另类丝袜制服| 1024香蕉在线观看| 中文字幕人成人乱码亚洲影| 日本欧美视频一区| 国产精品,欧美在线| 一a级毛片在线观看| 99国产极品粉嫩在线观看| 精品国内亚洲2022精品成人| 亚洲一区高清亚洲精品| 国产精华一区二区三区| 在线观看免费午夜福利视频| 久久久久九九精品影院| a级毛片在线看网站| 中文字幕另类日韩欧美亚洲嫩草| 视频区欧美日本亚洲| 禁无遮挡网站| 国内精品久久久久久久电影| 一本久久中文字幕| 视频在线观看一区二区三区| 欧美在线一区亚洲| 亚洲九九香蕉| 国产成人精品无人区| 国产午夜精品久久久久久| 亚洲av美国av| 国产av一区二区精品久久| 久久国产精品男人的天堂亚洲| 久久亚洲精品不卡| 精品国产一区二区三区四区第35| 在线视频色国产色| www.自偷自拍.com| 一进一出抽搐gif免费好疼| 欧美日韩瑟瑟在线播放| 色综合婷婷激情| 99国产精品99久久久久| 91老司机精品| 人人澡人人妻人| 国产精品一区二区免费欧美| 亚洲欧美精品综合一区二区三区| 免费在线观看视频国产中文字幕亚洲| 午夜日韩欧美国产| xxx96com| 级片在线观看| 精品一区二区三区av网在线观看| 亚洲精品国产色婷婷电影| 国产成人啪精品午夜网站| 精品欧美一区二区三区在线| 国产精品亚洲av一区麻豆| 精品一区二区三区四区五区乱码| 99久久综合精品五月天人人| 成人国语在线视频| 性欧美人与动物交配| 亚洲成国产人片在线观看| 成人国语在线视频| www国产在线视频色| 最新在线观看一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 一进一出抽搐gif免费好疼| 一区二区三区高清视频在线| 色哟哟哟哟哟哟| 999精品在线视频| 久久婷婷成人综合色麻豆| 不卡一级毛片| 国产欧美日韩一区二区三区在线| 久久婷婷成人综合色麻豆| 亚洲成人久久性| 成人av一区二区三区在线看| 老司机靠b影院| 亚洲全国av大片| 亚洲精品美女久久久久99蜜臀| 99国产精品免费福利视频| 欧美黑人欧美精品刺激| 老汉色av国产亚洲站长工具| 啦啦啦 在线观看视频| 久久精品人人爽人人爽视色| 久久中文字幕一级| 成人免费观看视频高清| or卡值多少钱| 国内毛片毛片毛片毛片毛片| 欧美日韩乱码在线| 精品无人区乱码1区二区| 精品午夜福利视频在线观看一区| 可以免费在线观看a视频的电影网站| 色综合欧美亚洲国产小说| www.精华液| 最近最新中文字幕大全电影3 | 久久午夜亚洲精品久久| 精品熟女少妇八av免费久了| 黄色片一级片一级黄色片| 桃色一区二区三区在线观看| 国产单亲对白刺激| 日日夜夜操网爽| netflix在线观看网站| 久久精品亚洲精品国产色婷小说| 亚洲自偷自拍图片 自拍| 亚洲色图av天堂| 成人特级黄色片久久久久久久| 亚洲精品在线观看二区| 国产亚洲精品第一综合不卡| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲精品久久久久久毛片| 亚洲成人久久性| 日韩欧美国产一区二区入口| 国产av在哪里看| 久久狼人影院| 国产欧美日韩一区二区三| 人成视频在线观看免费观看| 最近最新免费中文字幕在线| 1024视频免费在线观看| 后天国语完整版免费观看| 极品人妻少妇av视频| 欧美成人免费av一区二区三区| 99国产精品一区二区三区| 久久国产精品影院| 国产高清视频在线播放一区| 国产三级在线视频| 国产精品香港三级国产av潘金莲| 国产精品久久久久久精品电影 | 亚洲中文av在线| 久久久久久大精品| 免费女性裸体啪啪无遮挡网站| 免费人成视频x8x8入口观看| 麻豆av在线久日| 久久热在线av| 丝袜美足系列| 亚洲国产精品合色在线| 女同久久另类99精品国产91| 老熟妇仑乱视频hdxx| 亚洲中文字幕一区二区三区有码在线看 | 国产日韩一区二区三区精品不卡| 久久国产亚洲av麻豆专区| 国产亚洲精品av在线| 亚洲熟女毛片儿| 国产精品久久久av美女十八| 欧美老熟妇乱子伦牲交| 女生性感内裤真人,穿戴方法视频| 免费看美女性在线毛片视频| cao死你这个sao货| 国产欧美日韩综合在线一区二区| 午夜久久久久精精品| 国产精品野战在线观看| 亚洲自偷自拍图片 自拍| 亚洲电影在线观看av| 国产av在哪里看| 精品国产美女av久久久久小说| 欧美成狂野欧美在线观看| 亚洲成a人片在线一区二区| 久久久久精品国产欧美久久久| 脱女人内裤的视频| 色婷婷久久久亚洲欧美| 悠悠久久av| 国产蜜桃级精品一区二区三区| 长腿黑丝高跟| 高清黄色对白视频在线免费看| 给我免费播放毛片高清在线观看| 两人在一起打扑克的视频| 久久久久精品国产欧美久久久| 丝袜美足系列| 露出奶头的视频| 亚洲aⅴ乱码一区二区在线播放 | 99久久久亚洲精品蜜臀av| 亚洲精品在线观看二区| 久久狼人影院| 欧美激情 高清一区二区三区| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 国产亚洲欧美98| 午夜成年电影在线免费观看| 久久影院123| 国产人伦9x9x在线观看| 极品人妻少妇av视频| 日日夜夜操网爽| 在线播放国产精品三级| 免费在线观看日本一区| 色播在线永久视频| 国产精品1区2区在线观看.| 亚洲国产欧美日韩在线播放| 1024香蕉在线观看| videosex国产| 久久久久久久精品吃奶| 一区在线观看完整版| www国产在线视频色| 黑人欧美特级aaaaaa片| 欧美乱码精品一区二区三区| 在线观看免费视频日本深夜| 男人舔女人的私密视频| 男女下面进入的视频免费午夜 | 熟女少妇亚洲综合色aaa.| www日本在线高清视频| 在线观看www视频免费| 亚洲欧美日韩无卡精品| 黄片小视频在线播放| 成人免费观看视频高清| 老鸭窝网址在线观看| 国产片内射在线| 在线观看舔阴道视频| 亚洲欧美激情在线| 91字幕亚洲| 久久精品aⅴ一区二区三区四区| 国产精品九九99| 精品国产一区二区久久| 亚洲第一电影网av| 欧美最黄视频在线播放免费| 老司机深夜福利视频在线观看| 黄网站色视频无遮挡免费观看| 成人国产一区最新在线观看| 男女下面进入的视频免费午夜 | 少妇粗大呻吟视频| 免费av毛片视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品美女久久久久99蜜臀| 中亚洲国语对白在线视频| 久久精品国产99精品国产亚洲性色 | 中出人妻视频一区二区| 色播在线永久视频| 国产男靠女视频免费网站| 欧美日韩中文字幕国产精品一区二区三区 | 满18在线观看网站| 操美女的视频在线观看| 久久青草综合色| 亚洲中文字幕日韩| 极品教师在线免费播放| 精品国产一区二区久久| 亚洲 欧美一区二区三区| 动漫黄色视频在线观看| 男女做爰动态图高潮gif福利片 | 日韩一卡2卡3卡4卡2021年| 又黄又爽又免费观看的视频| 日本欧美视频一区| 免费看美女性在线毛片视频| 久99久视频精品免费| 精品一区二区三区av网在线观看| 在线观看免费视频日本深夜| 男女下面进入的视频免费午夜 | 欧美在线黄色| 亚洲人成77777在线视频| 午夜精品久久久久久毛片777| 久久影院123| 国产精品亚洲一级av第二区| 精品一品国产午夜福利视频| 一进一出抽搐动态| 日日爽夜夜爽网站| 国产高清视频在线播放一区| 亚洲精品美女久久久久99蜜臀| 国产精华一区二区三区| 中文字幕最新亚洲高清| 国产乱人伦免费视频| 亚洲精品美女久久久久99蜜臀| 老司机午夜福利在线观看视频| 国产精品亚洲美女久久久| 免费看a级黄色片| 黑丝袜美女国产一区| 999久久久精品免费观看国产| 村上凉子中文字幕在线| 极品教师在线免费播放| 淫妇啪啪啪对白视频| 夜夜看夜夜爽夜夜摸| 欧美激情 高清一区二区三区| 又大又爽又粗| 午夜久久久久精精品| 亚洲成人免费电影在线观看| 国产av一区二区精品久久| 麻豆国产av国片精品| 成人欧美大片| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品成人免费网站| 欧美成人免费av一区二区三区| 亚洲中文字幕日韩| 国产一区二区三区综合在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 可以在线观看毛片的网站| 天天躁夜夜躁狠狠躁躁| 国产亚洲欧美在线一区二区| 热99re8久久精品国产| 99久久99久久久精品蜜桃| 黄片小视频在线播放| 成人手机av| 一级毛片女人18水好多| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久av网站| 免费在线观看黄色视频的| 久久久久久久久中文| 涩涩av久久男人的天堂| 露出奶头的视频| 桃红色精品国产亚洲av| 国产精品 欧美亚洲| 国产激情欧美一区二区| 男人舔女人的私密视频| 亚洲性夜色夜夜综合| www国产在线视频色| 一级a爱视频在线免费观看| 变态另类丝袜制服| 日本一区二区免费在线视频| 9191精品国产免费久久| 国产成人系列免费观看| 亚洲国产精品成人综合色| 日韩欧美国产在线观看| 变态另类成人亚洲欧美熟女 | 精品国产乱子伦一区二区三区| 一进一出抽搐gif免费好疼| 露出奶头的视频| 日本一区二区免费在线视频| 国产亚洲精品av在线| 村上凉子中文字幕在线| 欧美日韩乱码在线| 久久国产亚洲av麻豆专区| 老熟妇仑乱视频hdxx| 久久久精品欧美日韩精品| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看66精品国产| 亚洲av电影不卡..在线观看| 此物有八面人人有两片| 午夜精品在线福利| 欧美日本视频| 国产欧美日韩一区二区三| 操出白浆在线播放| 国产亚洲精品一区二区www| 婷婷六月久久综合丁香| a在线观看视频网站| 中出人妻视频一区二区| 激情视频va一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 国产欧美日韩综合在线一区二区| 熟妇人妻久久中文字幕3abv| 91麻豆精品激情在线观看国产| 热99re8久久精品国产| 国语自产精品视频在线第100页| 国产黄a三级三级三级人| 一卡2卡三卡四卡精品乱码亚洲| 丝袜在线中文字幕| av视频在线观看入口| 脱女人内裤的视频| 亚洲国产看品久久| av福利片在线| 欧美不卡视频在线免费观看 | 看黄色毛片网站| 久久婷婷人人爽人人干人人爱 | 欧美亚洲日本最大视频资源| 国产伦人伦偷精品视频| 欧美一级a爱片免费观看看 | 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区| 亚洲aⅴ乱码一区二区在线播放 | 黄色视频不卡| 法律面前人人平等表现在哪些方面| 日韩免费av在线播放| 一级片免费观看大全| 色播在线永久视频| 国产私拍福利视频在线观看| 99久久精品国产亚洲精品| 一区二区三区激情视频| 亚洲人成77777在线视频| av天堂久久9| 欧美av亚洲av综合av国产av| 国产av又大| 精品福利观看| 黑丝袜美女国产一区| 久久国产精品人妻蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩一区二区精品| 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看 | 搡老岳熟女国产| 亚洲人成电影免费在线| 久热爱精品视频在线9| 免费久久久久久久精品成人欧美视频| 女同久久另类99精品国产91| 91成年电影在线观看| 9色porny在线观看| 国产在线观看jvid| 别揉我奶头~嗯~啊~动态视频| 97人妻天天添夜夜摸| 国产av精品麻豆| 一级毛片女人18水好多| 久久婷婷人人爽人人干人人爱 | 亚洲 国产 在线| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 老汉色av国产亚洲站长工具| 法律面前人人平等表现在哪些方面| 91字幕亚洲| 大香蕉久久成人网| 757午夜福利合集在线观看| av有码第一页| 高清在线国产一区| 国产成人系列免费观看| 日本一区二区免费在线视频| 亚洲人成网站在线播放欧美日韩| 搡老妇女老女人老熟妇| 中文字幕高清在线视频| 深夜精品福利| 亚洲精品国产色婷婷电影| 国产精品久久久av美女十八| 18禁黄网站禁片午夜丰满| 国产亚洲av高清不卡| 国产伦人伦偷精品视频| 麻豆一二三区av精品| 亚洲男人天堂网一区| 国产成人啪精品午夜网站| 午夜影院日韩av| 亚洲人成电影免费在线| 亚洲成人免费电影在线观看| cao死你这个sao货| 一本大道久久a久久精品| 在线天堂中文资源库| 国产人伦9x9x在线观看| 美女午夜性视频免费| 一级片免费观看大全| 国产精品电影一区二区三区| 精品国产美女av久久久久小说| 男人操女人黄网站| 黑人操中国人逼视频| 亚洲熟妇熟女久久| 黑人欧美特级aaaaaa片| 国产免费av片在线观看野外av| 久久久国产成人精品二区| 亚洲av成人不卡在线观看播放网| 制服诱惑二区| 亚洲国产高清在线一区二区三 | 午夜福利一区二区在线看| 精品欧美一区二区三区在线| av中文乱码字幕在线| 久久久精品国产亚洲av高清涩受| 亚洲中文字幕一区二区三区有码在线看 | 欧美在线黄色| e午夜精品久久久久久久| 午夜福利高清视频| 欧美人与性动交α欧美精品济南到| 中文亚洲av片在线观看爽| 亚洲国产欧美一区二区综合| 色老头精品视频在线观看| 一夜夜www| 亚洲国产欧美网| 视频区欧美日本亚洲| 免费一级毛片在线播放高清视频 | 在线播放国产精品三级| av在线播放免费不卡| 老司机午夜福利在线观看视频| 国产亚洲精品第一综合不卡| 精品一区二区三区四区五区乱码| 天天躁夜夜躁狠狠躁躁| 99精品久久久久人妻精品| 丝袜在线中文字幕| 精品欧美国产一区二区三| 激情视频va一区二区三区| 桃红色精品国产亚洲av| 妹子高潮喷水视频| 久久久久九九精品影院| 久久香蕉精品热| 黄片大片在线免费观看| 黄色 视频免费看| 女警被强在线播放| 免费看a级黄色片| 多毛熟女@视频| 国产精品久久久久久精品电影 | 国产伦一二天堂av在线观看| 亚洲aⅴ乱码一区二区在线播放 | 大香蕉久久成人网| www.自偷自拍.com| 成人亚洲精品av一区二区| 深夜精品福利| 国产成人免费无遮挡视频| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 满18在线观看网站| 精品国产超薄肉色丝袜足j| bbb黄色大片| 亚洲自拍偷在线| 十分钟在线观看高清视频www| 男女做爰动态图高潮gif福利片 | 亚洲精品在线美女| 免费在线观看黄色视频的| 亚洲少妇的诱惑av| 少妇的丰满在线观看| 91麻豆精品激情在线观看国产| 一级毛片女人18水好多| 久久香蕉国产精品| 黑人操中国人逼视频| 精品国产乱码久久久久久男人| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产中文字幕在线视频| 久久亚洲精品不卡| 日本三级黄在线观看| 国产欧美日韩一区二区三| 两个人看的免费小视频| 国产一区二区三区综合在线观看| 黄色毛片三级朝国网站| 亚洲狠狠婷婷综合久久图片| 少妇熟女aⅴ在线视频| 午夜福利视频1000在线观看 | 在线av久久热| 国产不卡一卡二| 最新美女视频免费是黄的| videosex国产| 757午夜福利合集在线观看| 18禁裸乳无遮挡免费网站照片 | 9热在线视频观看99| 国产精华一区二区三区| 国产亚洲精品第一综合不卡| 国产精品 欧美亚洲| 欧美日本中文国产一区发布| 久久精品亚洲精品国产色婷小说| 一区二区三区精品91|