基于響應曲面法的柴油機SCR性能預測

聶學選，畢玉華，申立中，王 鵬，嚴 杰，彭益源

基于響應曲面法的柴油機SCR性能預測

聶學選1，畢玉華1※，申立中1，王 鵬1，嚴 杰1，彭益源2

（1. 昆明理工大學，云南省內燃機重點實驗室，昆明 650500；2. 昆明云內動力股份有限公司，昆明 650500）

針對柴油機選擇性催化還原（Selective Catalytic Reduction，SCR）系統在不同工況下運行時性能差異較大，搭建了帶SCR系統的柴油機測試臺架，在對SCR系統性能測試的基礎上，利用GT-POWER建立SCR系統模型，分析不同排氣溫度、不同排氣流量、不同氨氮比對SCR性能的影響，并基于Box-Behnken設計與響應面法對柴油機SCR系統進行了研究，以排氣溫度、排氣流量、氨氮比為變量因子，以NOX轉化效率與NH3逃逸率為優化目標進行響應曲面優化。結果表明：排氣溫度對SCR性能影響較大，250～450 ℃為SCR最佳轉化效率區間，NOX轉化效率均在80%以上，NH3逃逸率均在5%以內；排氣流量增加使NOX轉化效率下降，NH3逃逸率上升，排氣流量在200 kg/h以上尤為明顯，排氣流量每增加50 kg/h，NOX轉化效率平均下降3%，NH3逃逸率平均增加4%；氨氮比增加使得NOX轉化效率提升，同時NH3逃逸率增加，氨氮比在0.9以上能使NOX轉化效率保持較高水平，氨氮比在0.9以下能保證較低的NH3逃逸率，氨氮比的選擇尤為重要。根據響應曲面結果得出：不同的排氣溫度與排氣流量配合不同氨氮比可提高NOX轉化效率，降低NH3逃逸率，當排氣溫度為350 ℃，排氣流量為200 kg/h，氨氮比為1.0時，SCR性能最佳，NOX轉化效率達到96.4%，NH3逃逸率僅0.5%。該研究為SCR系統在柴油機不同工況下運行時的尿素控制提供有效的指導依據。

柴油機；響應曲面；SCR系統；NOX轉化效率；氨逃逸率

0 引 言

選擇性催化還原反應技術（SCR，Selective Catalytic Reduction）專門用于減少柴油機NOX排放的機外凈化手段[1-3]，技術成熟且效率高，在目前的柴油機國六排放控制技術中已成為不可缺少的氮氧化物（NOX）后處理技術[4]。然而SCR尿素噴射過多會造成NH3反應過剩，從而造成NH3泄漏，形成二次污染。NH3是大氣中的堿性氣體，與二氧化硫、氮氧化物迅速形成細顆粒物（PM2.5），對霧霾的形成和大氣污染有著重要影響[5-7]。為了保證尿素噴射后產生的NH3與NOX完全反應，既降低NOX排放，也不產生NH3泄漏，需要對發動機不同工況下的尿素噴射提供精準控制，提高其SCR系統的轉化效率，降低NH3逃逸率。而柴油機與SCR系統的協調控制卻極為復雜，柴油機在不同工況下運行時排氣狀態參數不同，這成為柴油機SCR系統在開發過程中必須考慮的因素[8-10]。

國內外學者對SCR性能的優化方法也頗有研究。Capetillo等[11]利用全因子試驗設計對SCR性能進行優化，以尿素噴射速度、尿素噴射角度、液滴直徑為變量因子，NH3均勻性與尿素液滴的壁膜厚度為優化目標，成功地預測了SCR系統的性能。Vedagir等[12]利用正交試驗方法和CFD仿真方法對SCR設計參數進行了研究，以排氣溫度、排氣流量、尿素噴射位置、尿素噴射角度以變量因子，以SCR入口NH3均勻性為優化目標進行優化，最終給出了最佳的NH3均勻性的組合參數。王國仰等[13-14]利用多目標遺傳算法對SCR系統氨覆蓋率進行優化，分析了催化器溫度、空速和SCR催化器入口NOX濃度對目標氨存儲的影響，最終優化了SCR系統NOX排放和NH3泄漏之間的此消彼長（tread-of）的關系。Wei等[15]基于狀態反饋控制提出了一種非線性模型預測控制算法，提高了滿足NH3排放小于10×10-6時的NOX轉化效率。仇滔等[16]在柴油發動機臺架上開展SCR系統催化箱的進出口溫度特性測試，并基于Mtalab/simulink模塊對SCR催化箱開展仿真計算，提出了催化箱出口溫度的計算模型。李軍等[17]基于ESC工況進行了動態SCR轉化效率溫度窗口試驗研究研究。Liu等[18]提出了一種基于支持向量機和遺傳算法的組合方法，建立了上下游NOX預測模型，優化了NOX排放和NH3逃逸檢測。Pio等[19]發明了一種增強型NH3-SCR反應，即通過向尿素水溶液中添加一定量的硝酸銨來促進低溫下的標準SCR反應，此方法可以在200～350℃溫度范圍內顯著提高NOX轉化率，同時有效抑制氨泄漏。Colombo等[20]通過仿真模型研究了NH3的存儲和泄漏過程。綜上所述，國內外學者針對SCR性能的多因素優化主要偏向于SCR結構、尿素噴射、氨存儲模型等優化方面，而針對柴油機排氣狀態參數與尿素之間的綜合性能影響規律的研究極少。排氣溫度、排氣流量的變化，會對SCR催化劑活性、化學反應過程造成影響，使得NOX轉化效率產生明顯差異[21-23]。因此探究發動機排氣狀態參數、不同尿素噴射量對NOX排放與NH3排放的影響規律有著非常重要的意義。

以滿足國五排放標準的高壓共軌柴油機為研究對象，搭建了帶SCR系統的柴油機測試臺架，在對SCR系統性能測試的基礎上，利用GT-POWER軟件分析了不同排氣溫度、排氣流量、氨氮比（ANR，Ammonia to Nitrogen oxides Ratio）對SCR性能的影響。利用響應曲面法以NOX轉化效率與NH3逃逸率為優化目標進行響應曲面優化，提出排氣溫度、排氣流量、氨氮比三因素交互作用的影響規律，通過對三因素的優化，提高了NOX轉化效率，降低了NH3逃逸率。該研究為SCR系統在柴油機不同工況下運行時的尿素控制提供有效的指導依據。

1 數值模型理論

SCR催化劑載體是圓柱形蜂窩狀多孔道結構，材料通常有陶瓷、堇青石及金屬等，載體用于承載催化劑，催化劑均勻涂敷在催化劑涂層，與流經的NH3與排氣混合氣體反應，從而減少NOX排放，載體的物性參數對SCR性能有重要影響，合適的載體選擇不僅能使SCR所占空間小，使用壽命長，還能保證SCR具有優異的性能并減小對發動機動力性和經濟性的影響。GT-POWER中催化劑模型的物理結構如圖1所示，催化劑模型中進行著一系列的物理和化學反應，在載體通道中進行氣相物質的對流、擴散和傳導，催化劑涂覆層進行化學反應、能量交換，并與基體進行熱量傳遞，NOX與NH3的混合氣體流經孔道后先擴散至催化劑外表面，進一步擴散至催化劑內表面后，NH3吸附在催化劑活性位上與NOX發生催化還原反應，生成H2O與N2，生成產物擴散至孔道隨排氣流出催化劑，NOX排放經過這一整個過程便被有效控制。

1.1 流體力學基礎

計算流體力學運用數值方法和計算機求解流體控制方程，并得到空間和時間離散位置處的數值解，是對流場離散的定量描述，流體動力學遵循三大基本方程：連續性方程、運動方程、能量方程，方程式見（1）、（2）、（3）。

式（3）中，為溫度，℃；為流體的傳熱系數；grad為梯度；S為流體黏性耗散項，表示當流體在流動過程中受到黏性摩擦力而使流體的機械能轉化成熱能。能量方程是能量守恒定律在流體力學中的運用，包含了熱力學定律，即流體的總能量等于動能與內能之和。

1.2 化學反應動力性模型

氣體在流經載體孔道表面與催化劑發生反應，反應過程遵循Eley-Rideal機理[24]，本質上是兩種氣態物質與固態催化劑的吸附氧化過程，催化劑能加快化學反應速度并控制化學反應方向，將氣態混合物中的有害成分轉化為無害物。在SCR催化器中涉及到NO、NO2、NH3、O2、H2O、N2、CO2等物質反應，具體反應機理如下：

SCR反應的反應速率主要由基于化學反應動力學的Arrhenius公式確定[25]。

式（4）為標準SCR反應，由于排氣中的NO與NO2的比例一般在9∶1以上，此時SCR 催化器中主要發生標準SCR反應。式（5）稱為快速SCR反應，此反應可以在較低溫度下進行，并且在較低溫度下反應速率是標準SCR反應的17倍[26]。提高NOx中NO2的比例可以使SCR在較低溫度下發生快速SCR反應，有利于提高NOx轉化率[27]。當（NO2）/（NOX）比例繼續增大時，反應所示的緩慢SCR反應占主導作用，NOX轉化效率降低，式（6）為慢速SCR反應。式（7）為NH3氧化反應。式（8）中，為反應的速率常數，s-l；為指前因子，s-1；a為反應活化能，kJ/mol；為摩爾氣體常數，8.314 J/(mol·K)；1為熱力學溫度，K。

1.3 SCR性能評價指標

SCR一維性能參數主要以NOX轉化效率與NH3逃逸率為評價指標，定義式如式（9）、（10）所示[28]。

NOX轉化效率定義為

NH3逃逸率定義為

2 模型建立及驗證

2.1 試驗設備

試驗所用的發動機是一臺滿足國五排放標準的高壓共軌直列四缸柴油機，其主要技術參數如表1所示，SCR主要參數如表2所示。臺架設備儀器主要有AVL PUMA OPEN V1.5測控系統，AVL Dynoroad 202/12交流電力測功機，AVL 753C柴油質量流量計，AVL 735S柴油溫控系統，AVL FTIR i60，AVL AMA i60。圖2為發動機臺架布置示意圖，將AVL AMA i60置于SCR前端監測SCR入口的排放物濃度，另將AVL FTIR置于SCR后端監測SCR出口的NOX排放及NH3排放。

表1 發動機基本參數

表2 SCR主要結構參數

2.2 模型建立

利用GT-POWER建立SCR一維模型，該模型分為5個部分，進口邊界模塊、出口邊界模塊、SCR催化器模塊、化學反應模塊、監控NOX轉化效率及NH3逃逸率模塊。根據試驗采集的發動機排放數據對模型進行參數設置，包括排氣溫度、排氣流量與排氣組分，載體參數根據SCR實際參數進行設置，SCR出口邊界設置相應的大氣壓力、大氣溫度與氣體組分，圖3為SCR一維模型圖。

2.3 模型驗證

通過臺架試驗采集的3個工況試驗值與仿真對比驗證SCR一維仿真模型的準確性。表3為A、B、C三個工況的邊界條件，圖4為3個工況下NOX轉化效率與NH3逃逸率的試驗值與仿真值對比，可以看出NOX轉化效率與NH3逃逸率的仿真值與試驗值在曲線規律及數值上都有良好的跟隨性和一致性，各點工況誤差均在5%以內，說明模型具有較高的準確性。

表3 仿真邊界條件

3 單因素影響研究

為研究單因素對SCR性能的影響，在眾多因素中選取排氣溫度、排氣流量、氨氮比3個影響較大的因素進行研究。

3.1 排氣溫度的影響

通過一維仿真模型研究不同溫度對SCR性能的影響，在排氣流量200 kg/h，氨氮比為0.9時選取溫度150～500 ℃進行研究。圖5為不同排氣溫度下NOX轉化效率與NH3逃逸率，排氣溫度在150～250 ℃時，NOX轉化效率較低，且隨著溫度升高，NOX轉化效率呈現急劇上升的趨勢，而NH3逃逸率呈現相反的趨勢。低溫時由于催化劑活性較低，反應速率較慢，導致NH3參與反應較少，造成NH3較高的逃逸。因此，低溫時可通過提高排氣溫度，或者適當降低尿素噴射量，以此提高NOX轉化效率，降低NH3逃逸率。從圖中可看出在250～450 ℃為SCR最佳反應區間，此時SCR反應器內發生快反應，以式（5）為主導，在此區間內，NOX轉化效率均在80%以上，而NH3逃逸率均在5%以內。超過450 ℃時，反應以式（7）為主導，NH3與O2發生氧化反應，參與NOX反應的NH3變少，NOX轉化效率下降[29]。

圖5 不同排氣溫度下NOX轉化效率與NH3逃逸率

3.2 排氣流量的影響

選取排氣溫度250 ℃，氨氮比為0.9，排氣流量選取100～400 kg/h進行研究。圖6為不同排氣流量下NOX轉化效率與NH3逃逸率圖，排氣流量在100～200 kg/h時，NOX轉化效率保持在85%左右，NH3逃逸率保持在5%以內。排氣流量高于200 kg/h時，排氣流量每增加50 kg/h，NOX轉化效率平均下降3%，NH3逃逸率平均增加4%；排氣流量達到400 kg/h時，此時NOX轉化效率下降到74%，NH3逃逸率上升到16%。排氣流量加大時，NOX與NH3在催化劑上停留時間變短，因此，排氣流量加大導致NOX轉化效率降低，NH3逃逸率升高。

圖6 不同排氣流量下NOX轉化效率與NH3逃逸率

3.3 氨氮比的影響

選取排氣溫度250 ℃，排氣流量選取200 kg/h，氨氮比選取0.5～1.2進行研究。圖7為不同氨氮比下NOX轉化效率與NH3逃逸率圖，從圖中看出氨氮比提高，NOX轉化效率提高，同時NH3逃逸率也隨之升高。原因是氨氮比越高，參與NOX反應的NH3量越多，則NOX還原成N2和H2O較為徹底，因此NOX轉化效率提高，而氨氮比增大使NOX降低，參數NH3反應的NOX量不夠，從而造成NH3逃逸率變大。氨氮比在0.9以上能保持較高的轉化效率，氨氮比在1.2時轉化效率最高能達到95%，而氨氮比在0.9以下能保持較低的NH3逃逸率。綜合來看，為了保持較高的NOX轉化效率，較低的NH3逃逸率，氨氮比在0.9時最佳，此時NOX轉化效率為86.8%，NH3逃逸率為1.97%。

圖7 不同氨氮比下NOX轉化效率與NH3逃逸率

4 響應曲面優化

響應曲面法（Response Surface Methodology，RSM）是一種解決多變量因素之間的交互作用的試驗設計方法，通過試驗得到基礎數據，并最終通過多元回歸方程來擬合因子與所關心的響應值之間的函數關系，建立連續變量響應曲面模型，從而對因子之間的交互作用進行評價和分析。而采用的Box-Behnken試驗設計是響應曲面設計的方法之一，它可以評價指標和因素間的非線性關系的一種試驗設計方法。該試驗設計主要是適用于3～7個范圍內的試驗設計，其主要的優點試驗次數少，效率高，且所有因子不會同時處于高水平。

4.1 試驗設計

在單因素試驗的基礎上，采用響應面分析法中的Box-Behnken Design建立數學模型，柴油機轉速工作范圍為800～3 200 r/min，不同轉速下影響規律一致，故選取了中間轉速2 000r/min進行響應曲面分析。以NOX轉化效率與NH3逃逸率為指標進行優化試驗設計，選取排氣溫度（1）、排氣流量（2）、氨氮比（3）設計三因素三水平響應面試驗，利用GT-POWER算出計算結果。試驗因素及水平見表4，計算結果見表5。

表4 試驗因素水平表

表5 計算結果

4.2 方差分析

在響應曲面模型中，二階多項式擬合足夠表達試驗設計中的真實響應參數[30]，結合對表4進行回歸方差分析（Analysis of Variance，ANOVA），對因子及響應參數之間進行二次多項式擬合，最終得到NOX轉化效率（1）回歸方程與NH3逃逸率（2）回歸方程如式（11）、式（12）所示。

1=-150.67+0.711+0.092+177.533+1.55×10-41+

0.1713-0.0623-1.26×10-312-2.12×10-422-

95.6932（11）

2=178.61-0.921-0.062-33.73-1.32×10-412-

0.0813+0.0723+1.37×10-312+1.22×10-422+

34.2932（12）

由表6知，建立的NOX轉化效率模型回歸顯著（<0.000 1），模型的決定系數2為0.954 1，接近于1，說明模型與仿真數據具有較高的可信度，可用于NOX轉化效率的預測。建立的NH3逃逸率模型回歸顯著（<0.000 1），模型的決定系數2為0.987 5，接近于1，說明模型與仿真數據具有較高的可信度，可用于NH3逃逸率的預測。其NOX轉化效率模型與NH3逃逸率模型的計算值與預測值對比如圖8所示。

表6 方差分析結果

4.3 結果與討論

4.3.1 交互作用對NOX轉化效率的影響

圖9為排氣溫度、排氣流量與氨氮比三因素交互作用對NOX轉化效率的3D響應曲面圖。

排氣溫度與排氣流量在氨氮比為零水平時對NOX轉化效率的影響響應曲面如圖9a所示。當氨氮比水平為0.9時，排氣溫度占主導因素，排氣溫度在250～450 ℃時，NOX轉化效率保持較高水平，均在80%以上；300～400 ℃時NOX轉化效率能達到90%以上。相比排氣溫度而言，排氣流量對NOX轉化效率的影響趨勢較小，NOX轉化效率的高效率區向低排氣流量下傾斜?？傮w來看，高排氣溫度低排氣流量是NOX轉化效率的最佳選擇，在此試驗中，排氣溫度在350 ℃，排氣流量在200 kg/h，氨氮比為1.0時NOX轉化效率最佳，達到96.4%。

排氣溫度與氨氮比在排氣流量為零水平時對NOX轉化效率的影響響應曲面如圖9b所示。當排氣流量水平為250 kg/h時，NOX轉化效率高效率區向高排氣溫度高氨氮比方向傾斜。因此，調控合適的排氣溫度，選擇高的氨氮比能使NOX轉化效率保持較高水平，在此試驗中，排氣溫度在350 ℃，氨氮比在1.0以上NOX轉化效率能高達到95%，但為了兼顧NH3逃逸率氨氮比也應合理選擇。

排氣流量與氨氮比在排氣溫度為零水平時對NOX轉化效率的影響響應曲面如圖9c所示。當排氣溫度水平為300 ℃時，NOX轉化效率的高效率區域向低排氣流量高氨氮比方向傾斜。在低排氣流量時，NOX在催化劑上停留時間久，高的氨氮比保證足夠的NH3與NOX反應，從而提高NOX轉化效率。但是較高的氨氮比會導致NH3逃逸率升高，因此，氨氮比應適宜。在此試驗中，排氣流量在300 kg/h以下，氨氮比在0.9以上為NOX轉化效率的高效率區。

4.3.2 交互作用對NH3逃逸率的影響

圖10為排氣溫度、排氣流量與氨氮比三因素交互作用對NH3逃逸率的3D響應曲面圖。

排氣溫度與排氣流量在氨氮比為零水平時對NH3逃逸率的影響響應曲面如圖10a所示。當氨氮比水平為0.9時，排氣溫度仍然占主導因素，高溫時，由于化學反應速率加快，NOX與NH3反應完全，參與轉化的NH3多，因此在250 ℃以上時，NH3逃逸率低，均在15%以下。在200 ℃以下，NH3逃逸率較高，均達到了40%以上。低溫時，排氣流量越大，NH3逃逸率越高，原因是排氣流量大時，NH3附著在催化劑上的量變少，以及附著時間變短，導致參與反應的NH3變少?？傮w來看，低排氣溫度高排氣流量造成NH3逃逸率升高，為優化NH3逃逸率，應提高其排氣溫度，降低排氣流量。在此試驗中，為保證高效的NOX轉化效率，排氣溫度在350 ℃，排氣流量在200 kg/h時NOX轉化效率最佳，NH3逃逸率也能保持較低水平。

排氣溫度與氨氮比在排氣流量為零水平時對NH3逃逸率的影響響應曲面如圖10b所示。當排氣流量水平為250 kg/h時，低排氣溫度高氨氮比使得NH3逃逸率升高。低溫時由于催化劑活性不足，再高的氨氮比也導致化學反應不完全，從而造成較高的NH3逃逸率。因此，調控合適的排氣溫度與氨氮比能降低NH3逃逸率。在此試驗中，排氣溫度在300 ℃以上，氨氮比在0.9以上時能使NOX轉化效率保持較高水平，同時NH3逃逸率也保持較低水平。當排氣溫度為350 ℃，排氣流量為200 kg/h，氨氮比為1.0時，NH3逃逸率僅0.5%。

排氣流量與氨氮比在排氣溫度為零水平時對NH3逃逸率的影響響應曲面如圖10c所示。當排氣溫度水平為300 ℃時，選擇低的排氣流量低氨氮比能使NH3逃逸率保持最低。

5 結 論

1）排氣溫度對SCR性能影響較大，隨排氣溫度升高NOX轉化效率呈現先升高后下降的趨勢，NH3逃逸率呈現相反趨勢。排氣溫度250～450 ℃為NOX最佳轉化效率區間，NOX轉化效率均在80%以上，NH3逃逸率均在5%以內。

2）排氣流量增加使NOX轉化效率下降，NH3逃逸率上升，排氣流量在200 kg/h以上尤為明顯，排氣流量每增加50 kg/h，NOX轉化效率平均下降3%，NH3逃逸率平均增加4%。

3）氨氮比增加使得NOX轉化效率提升，同時NH3逃逸率升高，氨氮比在0.9以上能使NOX轉化效率保持較高水平，氨氮比在0.9以下能保證NH3逃逸率保持較低水平。

4）由響應曲面優化結果得出：高排氣溫度低排氣流量配合適宜的氨氮比能使NOX轉化效率保持高效率區，NH3逃逸率保持低水平區域。當排氣溫度在350 ℃，排氣流量在200 kg/h，氨氮比在1.0時，SCR性能最佳NOX轉化效率達到96.4%，NH3逃逸率僅0.5%。

[1] 孟忠偉，晉兆祥，陳超，等. DOC/SCR組合方式及關鍵參數對NOX轉化效率的影響[J]. 環境工程學報，2017，11(10)：5483-5489.

Meng Zhongwei, Jin Zhaoxiang, Chen Chao, et al. Influence of combinations and key parameters on NOXconversion efficiency in DOC/SCR system[J]. Chinese Journal of Environmental Engineering, 2017, 11(10): 5483-5489. (in Chinese with English abstract)

[2] In L D, Soo Y Y, Junhong P, et al. NOXconversion efficiency of SCR diesel vehicle under cold start condition[J]. Journal of ILASS-Korea, 2018, 23(4): 244-253.

[3] Vignesh R, Ashok B. Critical interpretative review on current outlook and prospects of selective catalytic reduction system for De-NOXstrategy in compression ignition engine[J]. Fuel, 2020, 276(9): 117996.

[4] 帥石金，劉世宇，馬驍，重型柴油車滿足近零排放法規的技術分析[J]. 汽車安全與節能學報，2019，10(1)：16-31.

Shuai Shijin, Liu Shiyu, Ma Xiao, et al. Technology analysis of heavy-duty diesel vehicles to meet nearzero emission regulations[J]. Automotive Safety and Energy, 2019, 10(1): 16-31. (in Chinese with English abstract)

[5] 薛文博，許艷玲，唐曉龍，等. 中國氨排放對PM2.5污染的影響[J]. 中國環境科學，2016，36(12)：3531-3539.

Xue Wenbo, Xu Yanling, Tang Xiaolong, et al. Impacts of ammonia emission on PM2.5pollution in China[J]. China Environmental Science, 2016, 36(12): 3531-3539. (in Chinese with English abstract)

[6] Ho K F, Ho S S H, Huang R J, et al. Chemical composition and bioreactivity of PM2.5during 2013 haze events in China[J]. Atmospheric Environment, 2016, 126(2): 162-170.

[7] Feng L, Liao W. Legislation, plans, and policies for prevention and control of air pollution in China: Achievements, challenges, and improvements[J]. Journal of Cleaner Production, 2016, 112(2): 1549-1558.

[8] Oh J, Lee K. Spray characteristics of a urea solution injector and optimal mixer location to improve droplet uniformity and NOx conversion efficiency for selective catalytic reduction[J]. Fuel, 2014, 119(3): 90-97.

[9] Chu X F, Chen L S, Zhou J W, et al. Experimental study of cu&Fe molecular sieve catalyst for urea-SCR after-treatment system[J]. Advanced Materials Research, 2013, 2694(9): 1843-1846.

[10] Byeongkil S, Hwan C H, Sun C J, et al. Physico-chemical property and catalytic activity of a CeO2-Doped MnOx–TiO2catalyst with SO2resistance for low-temperature NH3-SCR of NOX[J]. Journal of Nanoscience and Nanotechnology, 2016, 16(5): 4370-4376.

[11] Capetillo A, Ibarra F. Multiphase injector modelling for automotive SCR systems: A full factorial design of experiment and optimization[J]. Computers & Mathematics with Applications, 2017, 74(1): 188-200.

[12] Vedagiri P, Martin L J, Varuvel E G. Characterization of urea SCR using Taguchi technique and computational methods[J]. Environmental Science & Pollution Research, 2020(1): 11988–11999.

[13] 王國仰，祁金柱，劉世宇，等. 基于多目標遺傳算法的SCR系統氨覆蓋率優化[J]. 汽車工程，2020，42(3)：279-285.

Wang Guoyang, Qi Jinzhu, Liu Shiyu, et al. Optimization of ammonia coverage ratio in selective catalytic reduction system using a multi-objective genetic algorithm[J]. Automotive Engineering, 2020, 42(3): 279-285. (in Chinese with English abstract)

[14] Wang G, Qi J, Liu S, et al. Zonal control for selective catalytic reduction system using a model-based multi-objective genetic algorithm[J]. International Journal of Engine Research, 2021, 22(3): 911-920.

[15] Wei L, Yan F, Hu J, et al. NOx conversion efficiency optimization based on NSGA-II and state-feedback nonlinear model predictive control of selective catalytic reduction system in diesel engine[J]. Applied Energy, 2017. 206(11): 959-971

[16] 仇滔，宋鑫，雷艷，等. 柴油機NOX選擇性催化還原催化箱出口溫度特性及其計算模型[J]. 農業工程學報，2016，32(6)：89-94.

Qiu Tao, Song Xin, Lei Yan, et al. Outlet temperature characteristics for NOXcatalyst container of diesel engine SCR system and its calculation model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(6): 89-94. (in Chinese with English abstract)

[17] 李軍，李志軍，陶漢國，等. 基于ESC測試循環的SCR催化劑動態轉化效率溫度窗口研究[J]. 內燃機工程，2017，38(4)：41-46.

Li Jun, Li Zhijun, Tao Hanguo, et al. Research on the temperature window of SCR catalyst dynamic conversion efficiency on the ESC test[J]. Chinese Internal Combustion Engine Engineering, 2017, 38(4): 41-46. (in Chinese with English abstract)

[18] Liu B, Yan F, Hu J, et al. Modeling and multi-objective optimization of NOXconversion efficiency and NH3slip for a diesel engine[J]. Sustainability, 2016, 8(5): 478-478.

[19] Pio F, Isabella N, Enrico T. Removal of NOXfrom diesel exhausts: The new “Enhanced NH3-SCR” reaction[J]. Sae International Journal of Fuels & Lubricants, 2010, 3(1): 654-663.

[20] Colombo M, Koltsakis G, Nova I, et al. Modelling the ammonia adsorption–desorption process over an Fe–zeolite catalyst for SCR automotive applications[J]. Catalysis Today, 2012, 188(1): 42-52.

[21] 蘇慶運，馮廷智，鄭貫宇，等. 柴油機SCR催化轉化器的數值模擬[J]. 內燃機學報，2017，35(1)：53-59.

Su Qingyun, Feng Tingzhi, Zheng Guanyu, et al. Numerical simulation of diesel SCR system[J]. Transactions of CSICE 2017, 35(1): 53-59. (in Chinese with English abstract)

[22] Barman J, Aroras S, Khan R, et al. Temperature based model approach to optimize SCR calibration for BSIV norms[C]//The 12th International Conference on Automotive Engineering. 2016.

[23] Hirose S, Yamamoto H, Suenobu H, et al. Development of high porosity cordierite honeycomb substrate for SCR application to realize high NOXconversion efficiency and system compactness[J]. SAE International Journal of Materials and Manufacturing, 2014, 7(3): 682-687.

[24] 辛喆，王順喜，張寅，等. Urea-SCR催化器壓力損失及其對柴油機性能的影響[J]. 農業工程學報，2011，27(8)：169-173.

Xin Zhe, Wang Shunxi, Zhang Yin, et al. Pressure loss of urea-SCR converter and its influence on diesel engine performance[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(8): 169-173. (in Chinese with English abstract)

[25] 譚丕強，杜加振，胡志遠，等. 關鍵運行參數對柴油機SCR系統性能的影響[J]. 化工學報，2014，65(10)：4063-4070.

Tan Piqiang, Du Jiazhen, Hu Zhiyuan, et al. Effects of key operating parameters on SCR performance of diesel engines[J]. CIESC Journal, 2014, 65(10): 4063-4070. (in Chinese with English abstract)

[26] 蘇慶運，劉軒江，佟德輝，等. 柴油機Fe-Cu分子篩SCR催化劑的數值模擬[J]. 內燃機學報，2020，38(2)：161-168.

Su Qingyun, Liu Xuanjiang, Tong Dehui, et al, Numerical simulation of diesel combined Fe-Cu zeolite SCR Catalyst [J]. Transactions of CSICE, 2020, 38(2): 161-168. (in Chinese with English abstract)

[27] 張東劍，馬其華，張璐璐，等. 柴油機Urea-SCR系統內流場與化學反應耦合研究[J]. 內燃機工程，2018，39(6)：59-66.

Zhang Dongjian, Ma Qihua, Zhang Lulu, et al. Coupling research of flow field and chemical reaction in a diesel urea-SCR system[J]. Chinese Internal Combustion Engine Engineering, 2018, 39(6): 59-66. (in Chinese with English abstract)

[28] 胡廣地，馬盧平，劉叢志，等. 柴油機DOC+SCR系統數值仿真與結構優化[J]. 內燃機學報，2019，37(2)：155-163.

Hu Guangdi, Ma Luping, Liu Congzhi, et al. Numerical simulation and structural optimization of DOC+ SCR system for diesel engine[J]. Transactions of CSICE, 2019, 37(2): 155-163. (in Chinese with English abstract)

[29] 辛喆，張寅，王順喜，等. 柴油機Urea-SCR催化器轉化效率影響因素研究[J]. 農業機械學報，2011，42(9)：30-34.

Xin Zhe, Zhang Yin, Wang Shunxi, et al. Influence factors of DeNOXbehavior in urea SCR catalytic converter of diesel[J]. Transactions of The Chinese Society for Agricultural Machinery, 2011, 42(9): 30-34. (in Chinese with English abstract)

[30] 王俊，申立中，楊永忠，等. 基于響應曲面法的非道路用高壓共軌柴油機設計點優化標定[J]. 農業工程學報，2017，33(3)：31-39.

Wang Jun, Shen Lizhong, Yang Yongzhong, et al. Optimizing calibration of design points for non-road high pressure common rail diesel engine base on response surface methodology[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(3): 31-39. (in Chinese with English abstract)

Prediction of SCR performance of diesel engine based on response surface methodology

Nie Xuexuan1, Bi Yuhua1※, Shen Lizhong1, Wang Peng1, Yan Jie1, Peng Yiyuan2

(1.,,650500,; 2..,.,650500,)

Selective Catalytic Reduction (SCR) is applied to a diesel engine for the promising solution of NOXemission. In this study, a diesel engine test bench with an SCR system was developed to clarify the large difference in the NOXemission of diesel engine SCR system under various exhaust conditions. A performance test of the SCR system was also used to determine the single and multi-factor interaction. An SCR model was established using GT-POWER software. A systematic analysis was made to explore the influence of exhaust temperature, exhaust mass flow, and ammonia nitrogen ratio on the SCR performance in the heavy-duty engine. Box-Behnken design and response surface method (RSM) was used to simulate the diesel engine SCR system. An RSM optimization was carried out with the NOXconversion efficiency and NH3slip rate as optimization objectives, where the exhaust temperature, exhaust mass flow, and ammonia nitrogen ratio were variable factors. The results showed that the NOXconversion efficiency increased in the range of 150-250 °C, while a high level was then observed in the range of 250-450 °C, finally to decline after 450 °C. There was an opposite trend for the influence of exhaust temperature on NH3slip rate. Specifically, the NH3slip rate remained at a low level, all within 5% above 250 °C. The NOXconversion efficiency decreased, but the NH3slip rate increased, with the increase of exhaust mass flow, especially when the exhaust mass flow was above 200 kg / h. When the exhaust mass flow increased by 50 kg / h, the NOXconversion efficiency decreased by 3%, and the NH3slip rate increased by 4%. The high ammonia nitrogen ratio contributed to improving the NOXconversion efficiency and NH3slip rate. Particularly, the NOXconversion efficiency maintained a high level, when the ammonia nitrogen ratio was above 0.9. Nevertheless, the NH3slip rate maintained a low level, when the ammonia nitrogen ratio was below 0.9. It inferred that the appropriate ammonia nitrogen ratio was expected to optimize the SCR performance. In the response surface optimization, a high exhaust temperature and low exhaust mass flow with a suitable ammonia nitrogen ratio can contribute to the NOXconversion efficiency in the high level, while the NH3slip rate in the low level. An optimal NOXconversion efficiency of SCR performance reached 96.4%, and the NH3slip rate was only 0.5%, when the exhaust temperature was 350 °C, while the exhaust mass flow rate was 200 kg / h, and the ammonia nitrogen ratio was 1.0. Consequently, an optimal combination of NOXconversion efficiency and NH3slip rate can be achieved under the optimization of exhaust temperature, exhaust flow and ammonia nitrogen ratio. This finding can provide effective guidance for urea control in an SCR system under different operating conditions of a diesel engine.

diesel; response surface method; SCR; NOx conversion efficiency; ammonia slip rate

10.11975/j.issn.1002-6819.2021.10.008

TK427

A

1002-6819(2021)-10-0064-09

聶學選，畢玉華，申立中，等. 基于響應曲面法的柴油機SCR性能預測[J]. 農業工程學報，2021，37(10)：64-72.doi：10.11975/j.issn.1002-6819.2021.10.008 http://www.tcsae.org

Nie Xuexuan, Bi Yuhua, Shen Lizhong, et al. Prediction of SCR performance of diesel engine based on response surface methodology[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(10): 64-72. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2021.10.008 http://www.tcsae.org

2021-01-25

2021-04-23

國家自然科學基金重點項目（52066008）；云南省基礎研究重點項目（2018FA030）；云南省重大科技專項計劃項目（2018ZE001）

聶學選，研究方向為內燃機排放控制。Email：1838813119@qq.com

畢玉華，教授，博士生導師，研究方向為內燃機現代設計和排放控制。Email：815569621@qq.com