周 鑫， 吳 佳

（浙江浙能臺(tái)州第二發(fā)電有限責(zé)任公司，浙江 臺(tái)州 317109）

火電廠SCR煙氣脫硝系統(tǒng)建模與噴氨量最優(yōu)控制

周 鑫， 吳 佳

（浙江浙能臺(tái)州第二發(fā)電有限責(zé)任公司，浙江 臺(tái)州 317109）

SCR（選擇性催化還原）是發(fā)電廠目前普遍采用的煙氣脫硝方法，脫硝系統(tǒng)的噴氨量不僅影響煙氣脫硝的效率，過(guò)量噴氨也會(huì)造成氨逃逸率升高，導(dǎo)致環(huán)境的二次污染。SCR系統(tǒng)反應(yīng)機(jī)理復(fù)雜，具有非線性、大慣性等特點(diǎn)，傳統(tǒng)PID控制方法無(wú)法實(shí)現(xiàn)噴氨量的精確控制。將KPLS（核偏最小二乘方法）與GA（遺傳算法）結(jié)合，提出了GA-KPLS建模方法，并建立了SCR系統(tǒng)模型。仿真結(jié)果表明，模型具有較好的學(xué)習(xí)及泛化能力。將SCR模型應(yīng)用于模型預(yù)測(cè)控制方法中，實(shí)時(shí)計(jì)算最優(yōu)噴氨量，實(shí)現(xiàn)了對(duì)噴氨量的精確控制。實(shí)驗(yàn)結(jié)果表明，此方法與傳統(tǒng)PID控制方法相比，顯著提高了脫硝率，同時(shí)降低了氨逃逸率。

選擇性催化還原；煙氣脫硝；核偏最小二乘；遺傳算法；最優(yōu)控制

0 引言

發(fā)電廠鍋爐燃燒產(chǎn)生的NOX排放到大氣中會(huì)對(duì)環(huán)境造成嚴(yán)重污染[1]。降低NOX排放的措施分為2種：一是燃燒中減排，即采用分級(jí)燃燒或改進(jìn)燃燒器等方法使燃煤產(chǎn)生的NOX量減少；二是燃燒后減排，即通過(guò)加裝脫硝裝置對(duì)尾部煙氣進(jìn)行脫硝處理。SCR（選擇性催化還原方法）煙氣脫硝系統(tǒng)反應(yīng)機(jī)理復(fù)雜，脫硝效率受到噴氨量、反應(yīng)溫度、煙氣速度、催化劑活性等因素影響，其中噴氨量對(duì)脫硝效果影響最大，是關(guān)鍵的可調(diào)因素。噴氨量過(guò)少會(huì)導(dǎo)致煙氣NOX含量超標(biāo)；噴氨量過(guò)多，不僅會(huì)提高氨逃逸率造成二次污染，同時(shí)也會(huì)增加脫硝成本。

隨著預(yù)測(cè)控制、核偏最小二乘、神經(jīng)網(wǎng)絡(luò)以及支持向量機(jī)等技術(shù)的發(fā)展和在工業(yè)上的成功應(yīng)用，可以充分利用發(fā)電廠歷史運(yùn)行數(shù)據(jù)建立SCR煙氣脫硝系統(tǒng)模型，同時(shí)采用先進(jìn)控制算法，實(shí)現(xiàn)噴氨量的最優(yōu)控制，在降低NOX排放的同時(shí)減少脫硝成本。采用KPLS（核偏最小二乘）方法建立SCR煙氣脫硝系統(tǒng)模型，利用遺傳算法優(yōu)化核函數(shù)參數(shù)，從而提高模型的學(xué)習(xí)和泛化能力，將SCR模型與預(yù)測(cè)控制方法相結(jié)合，設(shè)計(jì)出基于GA-KPLS的噴氨量最優(yōu)控制系統(tǒng)。利用發(fā)電廠實(shí)際運(yùn)行數(shù)據(jù)進(jìn)行仿真驗(yàn)證，實(shí)驗(yàn)結(jié)果表明GAKPLS噴氨量最優(yōu)控制系統(tǒng)在顯著降低NOX排放量的同時(shí)減少了氨逃逸率。

1 SCR噴氨煙氣脫硝系統(tǒng)裝置

某1 000 MW超超臨界機(jī)組SCR煙氣脫硝系統(tǒng)結(jié)構(gòu)如圖1所示。

圖1 SCR煙氣脫硝系統(tǒng)示意

該SCR噴氨煙氣脫硝系統(tǒng)采用TiO2作為催化劑。稀釋空氣與來(lái)自氨氣站的氨氣混合，經(jīng)噴氨控制閥調(diào)整噴氨流量后通過(guò)噴嘴噴出，并與煙氣充分混合?；旌蠚怏w在TiO2的催化作用下發(fā)生選擇性催化還原反應(yīng)，煙氣中的NOX與NH3發(fā)生反應(yīng)，生成無(wú)害的氮?dú)夂退?，從而達(dá)到煙氣脫硝的目的[2]。SCR脫硝系統(tǒng)主要反應(yīng)過(guò)程如下：

實(shí)際運(yùn)行中同時(shí)會(huì)發(fā)生一些副反應(yīng)并產(chǎn)生一些有害副產(chǎn)物：

過(guò)量噴氨不僅會(huì)增加SCR系統(tǒng)中氨氣濃度，提高副反應(yīng)速度，導(dǎo)致重新生成NOX，影響脫硝效率，同時(shí)，過(guò)量的NH3與SO3反應(yīng)生成硫酸氫氨等有害副產(chǎn)物，這些副產(chǎn)物會(huì)導(dǎo)致催化劑失活或者堵塞催化劑孔板，從而降低脫硝效率，還會(huì)腐蝕管道影響機(jī)組的安全運(yùn)行，同時(shí)，氨逃逸也會(huì)造成二次污染[2]。

機(jī)組運(yùn)行工況的不斷變化會(huì)導(dǎo)致煙氣中NOX濃度產(chǎn)生劇烈波動(dòng)，傳統(tǒng)PID（比例-積分-微分）控制效果并不理想。因此，必須設(shè)計(jì)更優(yōu)的自動(dòng)控制系統(tǒng)控制噴氨量，以最合理的噴氨量保證脫硝效率，同時(shí)減少不利副產(chǎn)物的生成，避免二次污染。

2 SCR煙氣脫硝系統(tǒng)模型

2.1 GA-KPLS模型

KPLS是將核函數(shù)技術(shù)與偏最小二乘算法結(jié)合，利用非線性函數(shù)提高輸入變量的維度，把輸入空間映射到高維特征空間并應(yīng)用偏最小二乘求解，從而構(gòu)建出一種解決非線性問(wèn)題的算法。核函數(shù)是核偏最小二乘算法中的關(guān)鍵因素，但如何選擇、確定核函數(shù)目前尚無(wú)理論參考。通常認(rèn)為徑向基函數(shù)一般情況下優(yōu)于其他常用的核函數(shù)，因此選用徑向基函數(shù)作為核函數(shù)。

徑向基核函數(shù)中包含一個(gè)參數(shù)σ，它的確定尚無(wú)統(tǒng)一方法，采用GA（遺傳算法）對(duì)其進(jìn)行尋優(yōu)。遺傳算法按照適者生存和優(yōu)勝劣汰的原理，經(jīng)過(guò)迭代計(jì)算求取最優(yōu)近似解，具有良好的全局搜索能力且計(jì)算速度快。

設(shè)變量X∈Rn×p，Y∈Rn，p為自變量個(gè)數(shù)，n為觀測(cè)樣本個(gè)數(shù)，對(duì)數(shù)據(jù)進(jìn)行歸一化處理。KPLS算法描述如下：

（1）計(jì)算核矩陣K。核矩陣元素計(jì)算公式為：

（2）隨機(jī)初始化Y的得分向量u。

（3）計(jì)算X的得分向量t并歸一化。

（4）計(jì)算權(quán)值向量c。

（5）計(jì)算Y的得分向量u并進(jìn)行歸一化。

（6）重復(fù)步驟（2）—（5）直至t收斂。

（7）縮小K與Y，重復(fù)步驟（2）—（5）直至提取p個(gè)t與u。

訓(xùn)練樣本擬合公式為：

式中：T與U為得分向量t與u構(gòu)成的矩陣。

預(yù)測(cè)樣本擬合公式為：

式中：Kt為預(yù)測(cè)核矩陣；Yt為模型預(yù)測(cè)結(jié)果。

2.2 仿真試驗(yàn)與結(jié)果分析

根據(jù)工程實(shí)際和理論分析，選取SCR入口NOX濃度、噴氨量、煙氣溫度和機(jī)組負(fù)荷4個(gè)影響脫硝效率的因素作為SCR煙氣脫硝系統(tǒng)模型的輸入量，SCR出口NOX濃度作為輸出量。從運(yùn)行歷史數(shù)據(jù)中選取150組作為樣本，其中100組作為訓(xùn)練樣本，50組作為預(yù)測(cè)樣本。采用遺傳算法對(duì)參數(shù)σ進(jìn)行尋優(yōu)，尋優(yōu)結(jié)果為σ=0.01。模型訓(xùn)練及預(yù)測(cè)效果如圖2與3所示。

圖2 KPLS模型訓(xùn)練效果

圖3 KPLS模型預(yù)測(cè)效果

模型精度評(píng)價(jià)指標(biāo)采用RMSE（均方根誤差）和MAPE（平均絕對(duì)百分比誤差），計(jì)算公式如下：

式中：yi為測(cè)量值為預(yù)測(cè)值；n為預(yù)測(cè)樣本的個(gè)數(shù)。

SCR系統(tǒng)KPLS模型的訓(xùn)練及預(yù)測(cè)精度分別為RMSET=3.958 4 mg/m3，RMSEP=4.499 5 mg/m3，MAPET=5.41%，MAPEP=5.71%，其中下標(biāo)T和下標(biāo)P分別代表訓(xùn)練和預(yù)測(cè)評(píng)價(jià)指標(biāo)。從圖2與3中可以看出，擬合值和預(yù)測(cè)值分布在理想直線附近，說(shuō)明擬合及預(yù)測(cè)效果較好，模型具有較強(qiáng)的學(xué)習(xí)能力和泛化能力。

3 噴氨量最優(yōu)控制的仿真

3.1 模型預(yù)測(cè)控制

模型預(yù)測(cè)控制是一種智能控制算法，根據(jù)模型預(yù)測(cè)的輸出相應(yīng)調(diào)整輸入量，通過(guò)極小化目標(biāo)函數(shù)求解輸入量的變化量，從而使實(shí)際輸出與設(shè)定值保持一致。模型預(yù)測(cè)控制算法被轉(zhuǎn)化為求解非線性優(yōu)化問(wèn)題：

式中：yp為模型輸出；yr為設(shè)定值；Δu為輸入量的變化量；N與M分別為預(yù)測(cè)長(zhǎng)度和控制長(zhǎng)度；q與r分別為輸入輸出的權(quán)重系數(shù)。

針對(duì)上面的非線性優(yōu)化問(wèn)題，采用PSO（粒子群優(yōu)化）算法求取最優(yōu)解。PSO算法具有簡(jiǎn)單、容易實(shí)現(xiàn)、調(diào)整參數(shù)少的優(yōu)點(diǎn)?；煦缡欠蔷€性系統(tǒng)中普遍存在的現(xiàn)象?；煦邕\(yùn)動(dòng)具有內(nèi)在隨機(jī)性、非規(guī)則有序性和遍歷性等特點(diǎn)。采用Logistic映射xn+1=4xn（1-xn）產(chǎn)生混沌變量，其中n為迭代次數(shù)；初始值x0∈（0，1）且x0≠0.25，0.5，0.75。采用這種方法生成PSO算法的初始種群，可以使種群內(nèi)粒子分布較為均勻，具有多樣性、遍歷性的特點(diǎn)，有利于全局搜索。

基于以上建模方法以及預(yù)測(cè)控制算法，提出一種SCR噴氨量最優(yōu)控制方法。首先應(yīng)用GAKPLS方法建立SCR系統(tǒng)模型，再利用模型以及粒子群算法對(duì)噴氨量進(jìn)行優(yōu)化，通過(guò)改變噴氨量來(lái)提高脫硝效率，達(dá)到最優(yōu)控制目的。控制算法流程如圖4所示。

3.2 實(shí)驗(yàn)結(jié)果和分析

以某1 000 MW超超臨界鍋爐SCR煙氣脫硝系統(tǒng)作為研究對(duì)象，選取150組運(yùn)行歷史數(shù)據(jù)作為訓(xùn)練及預(yù)測(cè)樣本，其中訓(xùn)練樣本個(gè)數(shù)為100，預(yù)測(cè)樣本個(gè)數(shù)為50。對(duì)訓(xùn)練及預(yù)測(cè)樣本進(jìn)行預(yù)處理，利用遺傳算法對(duì)核函數(shù)參數(shù)進(jìn)行尋優(yōu)，建立SCR煙氣脫硝系統(tǒng)的KPLS模型。噴氨量最優(yōu)控制采用粒子群算法求解非線性優(yōu)化問(wèn)題。其中，設(shè)定慣性因子 ω=0.729，學(xué)習(xí)因子 c1=c2=1.494 45，種群大小為15，迭代次數(shù)為50，噴氨量改變量的范圍設(shè)為0～100 kg/h，根據(jù)國(guó)家相關(guān)規(guī)定SCR出口NOX濃度設(shè)定值設(shè)為50 mg/m3。仿真實(shí)驗(yàn)結(jié)果如圖5與6所示。

圖4 噴氨量最優(yōu)控制算法流程

圖5 SCR出口NOX濃度

SCR脫硝系統(tǒng)反應(yīng)機(jī)理復(fù)雜，具有大慣性、非線性等特征，傳統(tǒng)PID方法無(wú)法取得良好的控制效果。從圖5—6可知，采用傳統(tǒng)PID控制時(shí)平均脫硝率只有61.45%，SCR出口NOX平均濃度為71.630 6 mg/m3。相比之下，模型預(yù)測(cè)控制不需要深入了解被控對(duì)象的內(nèi)部機(jī)理，有效提高了系統(tǒng)的魯棒性。同時(shí)，模型預(yù)測(cè)控制對(duì)噴氨量進(jìn)行滾動(dòng)尋優(yōu)，具有較好的動(dòng)態(tài)控制效果。對(duì)噴氨量進(jìn)行最優(yōu)控制后，平均脫硝率為73.13%，SCR出口NOX平均濃度為49.578 5 mg/m3，脫硝效果有了顯著提升。

由于模型預(yù)測(cè)控制的性能指標(biāo)中引入了噴氨量的變化量，使得SCR出口NOX濃度在跟蹤設(shè)定值的同時(shí)盡量減少噴氨量的變化，實(shí)現(xiàn)噴氨量的精確控制，避免閥門頻繁調(diào)整帶來(lái)的損耗及安全隱患。仿真實(shí)驗(yàn)結(jié)果如圖7與8所示。

圖6 SCR脫硝率

圖7 噴氨量

圖8 氨逃逸率

傳統(tǒng)PID控制方式僅利用SCR出口NOX濃度進(jìn)行反饋控制，很難解決工況波動(dòng)情況下SCR系統(tǒng)的非線性和大慣性問(wèn)題，結(jié)果造成過(guò)量噴氨，導(dǎo)致氨氣逃逸，對(duì)環(huán)境造成了二次污染。采用傳統(tǒng)PID控制時(shí)，平均噴氨量為47.266 4 kg/h，平均氨逃逸率為10.79%。采用最優(yōu)控制后，噴氨量得到了精確控制，平均噴氨量為55.598 0 kg/h，平均氨逃逸率為9.89%。平均氨逃逸率有所下降，說(shuō)明過(guò)量噴氨的情況得到了改善。同時(shí)，平均噴氨量有所上升，因?yàn)槊撓趼视辛孙@著提高，脫硝反應(yīng)所需的氨氣量也隨之上升，導(dǎo)致了噴氨量的增加。

4 結(jié)語(yǔ)

利用GA-KPLS方法建立了SCR煙氣脫硝系統(tǒng)模型，模型具有較好的學(xué)習(xí)及泛化能力，充分反映了SCR系統(tǒng)的非線性過(guò)程。將模型與預(yù)測(cè)控制方法相結(jié)合，利用PSO算法對(duì)噴氨量進(jìn)行最優(yōu)控制。仿真實(shí)驗(yàn)結(jié)果表明，模型預(yù)測(cè)控制在顯著提高脫硝率的同時(shí)減小了氨逃逸率，避免了過(guò)量噴氨，實(shí)現(xiàn)了噴氨量的最優(yōu)控制。模型預(yù)測(cè)噴氨量最優(yōu)控制具有一定的工程應(yīng)用價(jià)值，為噴氨量的精確控制提供了一種有效方法。

[1]王春林，周昊，李國(guó)能，等.基于遺傳算法和支持向量機(jī)的低NOX燃燒優(yōu)化[J].中國(guó)電機(jī)工程學(xué)報(bào)，2007（11）∶40-44.

[2]周洪煜.超臨界鍋爐煙氣脫硝噴氨量混結(jié)構(gòu)-徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)最優(yōu)控制[J].中國(guó)電機(jī)工程學(xué)報(bào)，2011（5）∶108-113.

[3]閻綱.基于最小二乘支持向量機(jī)的鋁電解過(guò)程建模與控制研究[D].長(zhǎng)沙：中南大學(xué)，2012.

[4]杜振，錢徐悅，何勝，等.燃煤電廠煙氣SCR脫硝成本分析與優(yōu)化[J].中國(guó)電力，2013，46（10）∶124-128.

[5]梁川，沈越.1000MW機(jī)組SCR煙氣脫硝系統(tǒng)優(yōu)化運(yùn)行[J].中國(guó)電力，2012，45（1）∶41-44.

[6]董建勛，王松嶺，李永華，等.選擇性催化還原煙氣脫硝過(guò)程數(shù)學(xué)模擬研究[J].熱能動(dòng)力工程，2007，22（5）∶569-573.

[7]周響球,楊晨.選擇性催化還原脫硝反應(yīng)器數(shù)學(xué)模型及仿真[J].重慶大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，30（6）∶39-43.

[8]曹志勇，譚城軍，李建中，等.燃煤鍋爐SCR煙氣脫硝系統(tǒng)噴氨優(yōu)化調(diào)整試驗(yàn)[J].中國(guó)電力，2011，44（11）∶55-58.

[9]NOVA I，LIETTI L，TRONCONI E，et al.Transient response method applied to the kinetic analysis of the De-NOX-SCR reaction[J].Chemical Engineering Science，2001，56（4）∶1229-1237.

[10]趙乾.SCR煙氣脫硝系統(tǒng)模擬優(yōu)化及噴氨量最優(yōu)控制[D].重慶：重慶大學(xué)，2012.

[11]LIETTI L，NOVA I，CAMURRI S，et al.Dynamics of the SCR-DeNOXreaction by the transient-response method[J]. Aiche Journal，1997，43（10）∶2559-2570.

[12]周荔丹.廣義預(yù)測(cè)控制算法及應(yīng)用研究[D].長(zhǎng)沙：湖南大學(xué)，2002.

[13]黃山，蔣鷺，王天才，等.神經(jīng)網(wǎng)絡(luò)與遺傳算法結(jié)合的球團(tuán)豎爐燃燒優(yōu)化[J].東南大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，42（1）∶88-93.

[14]董建勛.燃煤電廠SCR煙氣脫硝試驗(yàn)研究及數(shù)學(xué)模型建立[D].保定：華北電力大學(xué)，2007.

[15]穆朝絮，張瑞民，孫長(zhǎng)銀.基于粒子群優(yōu)化的非線性系統(tǒng)最小二乘支持向量機(jī)預(yù)測(cè)控制方法[J].控制理論與應(yīng)用，2010，27（2）∶164-168.

（本文編輯：陸 瑩）

Modeling of SCR Flue Gas Denitration System and Optimal Control of Spraying Ammonia Flow in Thermal Power Plant

ZHOU Xin，WU Jia
（Zhejiang Zheneng Taizhou Second Electric Power Generation Co.，Ltd.，Taizhou Zhejiang 317109，China）

Selective catalytic reduction（SCR）method is usually used for flue gas denitration in thermal power plant.Spraying ammonia flow of denitration system can affect the efficiency of flue gas denitration And excess ammonia spraying results in higher rates of ammonia escape which cause secondary pollution of the environment.Reaction mechanism of SCR system is very complex and it has the characteristics of nonlinearity and large inertia.Therefore，it is difficult for the traditional PID control methods to achieve precise control of the amount of ammonia injection.Combining the kernel partial least squares（KPLS）and genetic algorithm（GA），GA-KPLS modeling method is proposed and the SCR system model is established.The simulation results show that the learning and generalization abilities of the model are both better.In order to precisely control the amount of ammonia spraying，model predictive control method is used to calculate the real-time optimal amount of ammonia spraying.Experimental results show that compared with traditional PID control this method significantly improves the denitration rate as well as reduces the ammonia escape rate.

selective catalytic reduction；flue gas denitration；kernel partial least squares；genetic algorithm；optimal control

X701.3

B

1007-1881（2015）11-0015-05

2015-09-17

周 鑫（1990），男，助理工程師，從事發(fā)電廠機(jī)組運(yùn)行工作。