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

黃慶華，宋麗云，吳 銳，何 洪，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