• <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ì)算
    美女脱内裤让男人舔精品视频| 搡老岳熟女国产| 国产亚洲av片在线观看秒播厂| 国产亚洲一区二区精品| 一区二区三区乱码不卡18| 久久热在线av| 国产亚洲欧美在线一区二区| 婷婷成人精品国产| 一级毛片电影观看| 国产精品秋霞免费鲁丝片| 国产黄色免费在线视频| 国产男女超爽视频在线观看| 少妇精品久久久久久久| 久久久久久免费高清国产稀缺| 飞空精品影院首页| 久久精品熟女亚洲av麻豆精品| 性高湖久久久久久久久免费观看| 丝袜人妻中文字幕| 又大又爽又粗| 亚洲九九香蕉| 18禁国产床啪视频网站| 婷婷成人精品国产| 亚洲专区国产一区二区| av在线老鸭窝| 少妇裸体淫交视频免费看高清 | 国产男女超爽视频在线观看| 亚洲熟女精品中文字幕| 91大片在线观看| 天天躁日日躁夜夜躁夜夜| 99热全是精品| 国产精品久久久人人做人人爽| 99精国产麻豆久久婷婷| 亚洲伊人久久精品综合| netflix在线观看网站| 999久久久精品免费观看国产| 久久国产精品男人的天堂亚洲| 国产精品九九99| 日本a在线网址| 免费在线观看黄色视频的| 伦理电影免费视频| 精品福利观看| 免费黄频网站在线观看国产| 19禁男女啪啪无遮挡网站| 老鸭窝网址在线观看| 性少妇av在线| 午夜免费鲁丝| 又紧又爽又黄一区二区| 精品久久久久久电影网| 亚洲成国产人片在线观看| 老司机午夜十八禁免费视频| 国产一区二区在线观看av| av福利片在线| 亚洲精品乱久久久久久| 法律面前人人平等表现在哪些方面 | 纵有疾风起免费观看全集完整版| 国产免费福利视频在线观看| 成年美女黄网站色视频大全免费| 丝袜美腿诱惑在线| 国产免费福利视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲激情五月婷婷啪啪| 免费观看a级毛片全部| 国产欧美日韩一区二区精品| 亚洲一卡2卡3卡4卡5卡精品中文| 99热国产这里只有精品6| 日韩有码中文字幕| 色婷婷av一区二区三区视频| 在线看a的网站| 嫩草影视91久久| 无限看片的www在线观看| 亚洲欧洲日产国产| 高清在线国产一区| 国产男女内射视频| 一本色道久久久久久精品综合| 国产成人精品在线电影| 80岁老熟妇乱子伦牲交| 咕卡用的链子| 色视频在线一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品乱久久久久久| 在线看a的网站| 国产欧美日韩一区二区三 | 王馨瑶露胸无遮挡在线观看| 国产精品1区2区在线观看. | 午夜福利影视在线免费观看| 国产精品影院久久| 91av网站免费观看| 亚洲国产欧美在线一区| 日韩中文字幕视频在线看片| 亚洲avbb在线观看| 黑人猛操日本美女一级片| 亚洲第一欧美日韩一区二区三区 | 黄色视频不卡| 成年人免费黄色播放视频| 国产精品久久久人人做人人爽| 搡老岳熟女国产| 精品人妻1区二区| 午夜精品国产一区二区电影| 亚洲伊人色综图| 韩国精品一区二区三区| 亚洲精品国产av成人精品| 大香蕉久久成人网| 精品福利观看| 视频区图区小说| 男女无遮挡免费网站观看| 国产精品偷伦视频观看了| 亚洲五月色婷婷综合| 国产欧美日韩一区二区精品| 九色亚洲精品在线播放| 青草久久国产| 午夜福利乱码中文字幕| 亚洲精品乱久久久久久| 久久这里只有精品19| 首页视频小说图片口味搜索| 国产成人啪精品午夜网站| 欧美国产精品一级二级三级| 午夜免费鲁丝| 少妇粗大呻吟视频| 日韩视频一区二区在线观看| 国产成人欧美在线观看 | 高清在线国产一区| 飞空精品影院首页| 90打野战视频偷拍视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美激情久久久久久爽电影 | 9热在线视频观看99| 狠狠精品人妻久久久久久综合| 男女午夜视频在线观看| 亚洲第一av免费看| 久久精品国产a三级三级三级| 国产成人系列免费观看| 午夜福利在线观看吧| 正在播放国产对白刺激| 超碰97精品在线观看| 十八禁网站免费在线| 欧美日韩视频精品一区| 精品国产乱子伦一区二区三区 | 久久亚洲国产成人精品v| 国产一区二区三区综合在线观看| 成年人免费黄色播放视频| 久久毛片免费看一区二区三区| 人妻一区二区av| videosex国产| 欧美激情高清一区二区三区| 满18在线观看网站| 亚洲第一欧美日韩一区二区三区 | 午夜激情久久久久久久| 十分钟在线观看高清视频www| 黄色片一级片一级黄色片| 欧美久久黑人一区二区| 日韩制服丝袜自拍偷拍| 欧美在线一区亚洲| 久久久久久久精品精品| 精品亚洲成国产av| 一区二区三区精品91| 色视频在线一区二区三区| 亚洲精品国产av蜜桃| 国产精品av久久久久免费| 国产精品亚洲av一区麻豆| 国产精品熟女久久久久浪| h视频一区二区三区| 亚洲国产精品999| 黄片大片在线免费观看| 丁香六月欧美| 精品国产一区二区三区四区第35| 亚洲精品中文字幕一二三四区 | 欧美国产精品va在线观看不卡| 美女中出高潮动态图| 一级,二级,三级黄色视频| 精品第一国产精品| 91字幕亚洲| 少妇裸体淫交视频免费看高清 | 一级,二级,三级黄色视频| 成在线人永久免费视频| 女人精品久久久久毛片| 午夜激情av网站| 欧美老熟妇乱子伦牲交| 在线观看免费午夜福利视频| 成人影院久久| 成年人午夜在线观看视频| 亚洲专区中文字幕在线| 欧美日本中文国产一区发布| 免费看十八禁软件| 亚洲精华国产精华精| 99re6热这里在线精品视频| 免费在线观看完整版高清| 女人被躁到高潮嗷嗷叫费观| 欧美精品啪啪一区二区三区 | 可以免费在线观看a视频的电影网站| 国产黄频视频在线观看| 亚洲色图 男人天堂 中文字幕| 后天国语完整版免费观看| 高清视频免费观看一区二区| 久热爱精品视频在线9| av天堂在线播放| 国产亚洲av片在线观看秒播厂| 国产精品国产三级国产专区5o| 一边摸一边做爽爽视频免费| 国产一区有黄有色的免费视频| 国产免费现黄频在线看| 香蕉国产在线看| 高清在线国产一区| 美女国产高潮福利片在线看| 777久久人妻少妇嫩草av网站| 在线 av 中文字幕| 又黄又粗又硬又大视频| 在线观看舔阴道视频| 久久精品国产综合久久久| 亚洲va日本ⅴa欧美va伊人久久 | 巨乳人妻的诱惑在线观看| 亚洲五月色婷婷综合| 美女福利国产在线| 欧美一级毛片孕妇| 国产成人影院久久av| 在线天堂中文资源库| tube8黄色片| 一二三四在线观看免费中文在| 我的亚洲天堂| 免费不卡黄色视频| 男女免费视频国产| 久久亚洲国产成人精品v| 999久久久国产精品视频| 欧美 亚洲 国产 日韩一| 国产亚洲午夜精品一区二区久久| 国产精品偷伦视频观看了| 人人妻人人澡人人爽人人夜夜| 国产在线观看jvid| 99热全是精品| 亚洲成人手机| 国产男女内射视频| 久久香蕉激情| 51午夜福利影视在线观看| av视频免费观看在线观看| 12—13女人毛片做爰片一| 又黄又粗又硬又大视频| 欧美大码av| 国产一区二区在线观看av| av网站免费在线观看视频| 精品少妇黑人巨大在线播放| 日本91视频免费播放| 亚洲精品成人av观看孕妇| 久久国产亚洲av麻豆专区| 成人三级做爰电影| 两性夫妻黄色片| 99国产综合亚洲精品| 久久毛片免费看一区二区三区| 999久久久国产精品视频| 色老头精品视频在线观看| 国产成人av激情在线播放| 国产精品亚洲av一区麻豆| 青草久久国产| 久久人妻熟女aⅴ| 人妻一区二区av| 久久国产精品大桥未久av| 亚洲精品粉嫩美女一区| a级片在线免费高清观看视频| www.自偷自拍.com| 丝袜脚勾引网站| 美女福利国产在线| 男女国产视频网站| 日本vs欧美在线观看视频| 免费高清在线观看视频在线观看| 天天添夜夜摸| 国产又色又爽无遮挡免| 日本五十路高清| 一区二区三区激情视频| 免费在线观看视频国产中文字幕亚洲 | 99热网站在线观看| 午夜日韩欧美国产| 伊人久久大香线蕉亚洲五| 免费高清在线观看日韩| 亚洲伊人久久精品综合| 淫妇啪啪啪对白视频 | 超碰97精品在线观看| 欧美日韩国产mv在线观看视频| 一级毛片电影观看| 日韩一卡2卡3卡4卡2021年| 视频在线观看一区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲精品美女久久av网站| 久热这里只有精品99| 精品人妻一区二区三区麻豆| 超碰97精品在线观看| 日韩 亚洲 欧美在线| 咕卡用的链子| 脱女人内裤的视频| 成人影院久久| 亚洲avbb在线观看| 99国产精品一区二区三区| 人妻 亚洲 视频| 在线观看一区二区三区激情| 日韩制服丝袜自拍偷拍| 黄色 视频免费看| 一级黄色大片毛片| 国产老妇伦熟女老妇高清| 咕卡用的链子| 中文字幕高清在线视频| 免费在线观看日本一区| 成人18禁高潮啪啪吃奶动态图| 亚洲国产欧美在线一区| 国产麻豆69| tocl精华| 中文欧美无线码| 麻豆av在线久日| 精品福利永久在线观看| 欧美日韩亚洲综合一区二区三区_| 大香蕉久久网| 国产在线视频一区二区| 亚洲第一青青草原| 久久久精品免费免费高清| tube8黄色片| 精品亚洲成a人片在线观看| 老熟女久久久| 国产在线一区二区三区精| 成人国产一区最新在线观看| 久久久久久人人人人人| 啦啦啦啦在线视频资源| 亚洲五月色婷婷综合| 亚洲专区国产一区二区| 欧美日本中文国产一区发布| 亚洲国产欧美在线一区| 在线 av 中文字幕| 欧美+亚洲+日韩+国产| 亚洲av美国av| 欧美精品人与动牲交sv欧美| 亚洲欧美成人综合另类久久久| 久久热在线av| 悠悠久久av| 国产在线免费精品| 黄色视频,在线免费观看| 最近中文字幕2019免费版| 欧美日韩亚洲综合一区二区三区_| 精品国产乱码久久久久久男人| 久久人妻熟女aⅴ| cao死你这个sao货| 亚洲av欧美aⅴ国产| 久久久久久久久久久久大奶| av网站免费在线观看视频| 亚洲一区二区三区欧美精品| 久久精品成人免费网站| 日韩大码丰满熟妇| 如日韩欧美国产精品一区二区三区| 一进一出抽搐动态| 丰满少妇做爰视频| 热99久久久久精品小说推荐| 69精品国产乱码久久久| 热re99久久国产66热| 免费在线观看影片大全网站| 久久亚洲国产成人精品v| 国产精品 国内视频| 人妻久久中文字幕网| 黑人操中国人逼视频| 91字幕亚洲| 国产成人精品在线电影| 亚洲av日韩在线播放| 国产91精品成人一区二区三区 | 成人国产一区最新在线观看| 国产日韩欧美亚洲二区| 亚洲国产看品久久| 美女中出高潮动态图| 18禁国产床啪视频网站| 19禁男女啪啪无遮挡网站| 久久中文字幕一级| 永久免费av网站大全| 18在线观看网站| 免费在线观看黄色视频的| 国产97色在线日韩免费| 国产精品免费大片| 午夜免费鲁丝| 天天躁日日躁夜夜躁夜夜| a级毛片黄视频| 999久久久精品免费观看国产| 国产深夜福利视频在线观看| 国产精品久久久av美女十八| 欧美变态另类bdsm刘玥| 亚洲国产日韩一区二区| 国产av精品麻豆| 亚洲av欧美aⅴ国产| 波多野结衣一区麻豆| 国产区一区二久久| 欧美少妇被猛烈插入视频| 国精品久久久久久国模美| 中文字幕精品免费在线观看视频| 久热爱精品视频在线9| 18禁观看日本| 美女中出高潮动态图| 中文字幕人妻丝袜一区二区| 大陆偷拍与自拍| 久久精品成人免费网站| 热99re8久久精品国产| 国产精品免费大片| 日日夜夜操网爽| 欧美激情久久久久久爽电影 | 两个人免费观看高清视频| 亚洲久久久国产精品| 狠狠婷婷综合久久久久久88av| 免费高清在线观看视频在线观看| 热99re8久久精品国产| 午夜久久久在线观看| 真人做人爱边吃奶动态| 色老头精品视频在线观看| 悠悠久久av| 国产精品香港三级国产av潘金莲| 国产一区二区三区综合在线观看| 丝袜美腿诱惑在线| 免费高清在线观看视频在线观看| 性少妇av在线| 9色porny在线观看| 超碰成人久久| 久久久久久久大尺度免费视频| 少妇精品久久久久久久| 在线精品无人区一区二区三| 欧美xxⅹ黑人| 后天国语完整版免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久99热这里只频精品6学生| 精品人妻一区二区三区麻豆| 欧美亚洲 丝袜 人妻 在线| 日韩免费高清中文字幕av| 不卡一级毛片| 一区二区日韩欧美中文字幕| 亚洲中文av在线| 久久久久国内视频| 国产成人欧美在线观看 | 国产精品久久久久久精品古装| 国产精品久久久av美女十八| 黄色视频不卡| a在线观看视频网站| 久久久久视频综合| 国产免费视频播放在线视频| 久久国产精品影院| 一个人免费看片子| 国产在线免费精品| 汤姆久久久久久久影院中文字幕| 精品久久蜜臀av无| 美女福利国产在线| 激情视频va一区二区三区| 国产淫语在线视频| 国产一区二区三区在线臀色熟女 | 1024香蕉在线观看| 免费女性裸体啪啪无遮挡网站| 操出白浆在线播放| 青春草亚洲视频在线观看| 99国产精品一区二区三区| 男人舔女人的私密视频| 欧美黑人欧美精品刺激| 一级片免费观看大全| 精品久久蜜臀av无| 99九九在线精品视频| 在线亚洲精品国产二区图片欧美| 人人妻,人人澡人人爽秒播| 欧美黑人欧美精品刺激| 精品一区二区三区四区五区乱码| 超碰成人久久| 18禁国产床啪视频网站| 麻豆av在线久日| 婷婷丁香在线五月| 精品卡一卡二卡四卡免费| 丝袜美足系列| 男男h啪啪无遮挡| 少妇猛男粗大的猛烈进出视频| 一区二区三区四区激情视频| 午夜福利视频精品| 欧美乱码精品一区二区三区| 老司机深夜福利视频在线观看 | av天堂在线播放| av在线老鸭窝| 考比视频在线观看| 十八禁网站网址无遮挡| 性高湖久久久久久久久免费观看| 午夜福利在线免费观看网站| 中文字幕另类日韩欧美亚洲嫩草| 欧美在线一区亚洲| 久久九九热精品免费| 黑人猛操日本美女一级片| 国产在线观看jvid| 操美女的视频在线观看| 亚洲精品在线美女| 黑丝袜美女国产一区| 成人影院久久| 精品亚洲成国产av| 国产精品成人在线| 精品久久久精品久久久| 国产伦理片在线播放av一区| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看 | 午夜福利视频精品| 91成人精品电影| 天堂中文最新版在线下载| 国产男人的电影天堂91| 水蜜桃什么品种好| 欧美午夜高清在线| 亚洲天堂av无毛| 精品少妇内射三级| 丝袜喷水一区| 啪啪无遮挡十八禁网站| 中国国产av一级| 老司机靠b影院| 精品久久久久久久毛片微露脸 | 制服人妻中文乱码| 国产成人欧美| 一个人免费在线观看的高清视频 | 搡老乐熟女国产| 国产成人欧美| 日韩视频一区二区在线观看| 黑人操中国人逼视频| 五月开心婷婷网| 国产成人免费无遮挡视频| 人人妻人人爽人人添夜夜欢视频| 丰满迷人的少妇在线观看| 黄片小视频在线播放| 母亲3免费完整高清在线观看| 久久国产精品人妻蜜桃| 深夜精品福利| 午夜免费鲁丝| 人人妻人人爽人人添夜夜欢视频| 男女高潮啪啪啪动态图| 80岁老熟妇乱子伦牲交| 亚洲欧洲精品一区二区精品久久久| 蜜桃在线观看..| www.精华液| 国产免费av片在线观看野外av| 多毛熟女@视频| 99国产极品粉嫩在线观看| 国产精品久久久av美女十八| 国产精品久久久久久人妻精品电影 | 黄色 视频免费看| 叶爱在线成人免费视频播放| 国产在视频线精品| 国产av精品麻豆| 国产在视频线精品| 欧美xxⅹ黑人| 淫妇啪啪啪对白视频 | 欧美午夜高清在线| bbb黄色大片| 下体分泌物呈黄色| 亚洲精品美女久久av网站| 亚洲成国产人片在线观看| 亚洲欧美日韩另类电影网站| 女性生殖器流出的白浆| 久久久国产一区二区| 另类精品久久| 肉色欧美久久久久久久蜜桃| 制服人妻中文乱码| 精品人妻一区二区三区麻豆| 午夜福利在线观看吧| 精品少妇内射三级| 久久久精品免费免费高清| 精品熟女少妇八av免费久了| 女性被躁到高潮视频| 中文字幕人妻丝袜一区二区| 最新在线观看一区二区三区| 一区二区三区激情视频| 久久亚洲精品不卡| 欧美日韩福利视频一区二区| 热99国产精品久久久久久7| 亚洲av美国av| 老司机午夜福利在线观看视频 | 国产麻豆69| 麻豆乱淫一区二区| 香蕉丝袜av| 国产成人av激情在线播放| 日本欧美视频一区| 99香蕉大伊视频| 手机成人av网站| 一个人免费在线观看的高清视频 | 男男h啪啪无遮挡| 国产欧美日韩一区二区三 | 久久午夜综合久久蜜桃| 精品亚洲成国产av| 亚洲天堂av无毛| 国产精品一区二区精品视频观看| 视频区欧美日本亚洲| 中国国产av一级| 免费少妇av软件| 视频区图区小说| 亚洲欧美精品综合一区二区三区| 亚洲av电影在线进入| 精品第一国产精品| 美女午夜性视频免费| 亚洲人成电影免费在线| 九色亚洲精品在线播放| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美日韩高清在线视频 | 不卡一级毛片| 久久天躁狠狠躁夜夜2o2o| 一本大道久久a久久精品| 美女高潮到喷水免费观看| 免费av中文字幕在线| 午夜福利视频在线观看免费| 欧美人与性动交α欧美精品济南到| 一级a爱视频在线免费观看| 免费女性裸体啪啪无遮挡网站| 国产1区2区3区精品| 欧美黑人欧美精品刺激| 国产精品 欧美亚洲| 国产成人系列免费观看| 国产一区二区在线观看av| 蜜桃国产av成人99| 男女边摸边吃奶| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 亚洲精品av麻豆狂野| 美女高潮喷水抽搐中文字幕| 法律面前人人平等表现在哪些方面 | 日韩视频在线欧美| 成人国产av品久久久| 亚洲欧洲日产国产| 久久国产精品男人的天堂亚洲| 少妇人妻久久综合中文| 欧美日韩亚洲高清精品| 嫁个100分男人电影在线观看| 欧美精品亚洲一区二区| 国产福利在线免费观看视频| 一本综合久久免费|