RBF-BP網(wǎng)絡(luò)在火電廠氮氧化物濃度檢測(cè)的應(yīng)用*

嚴(yán) 玥, 嚴(yán) 實(shí), 江 赟

(1.重慶工商大學(xué) 檢測(cè)控制集成系統(tǒng)重慶市市級(jí)工程實(shí)驗(yàn)室,重慶 400067;2.重慶川儀分析儀器有限公司，重慶 400060)

?

應(yīng)用技術(shù)

RBF-BP網(wǎng)絡(luò)在火電廠氮氧化物濃度檢測(cè)的應(yīng)用*

嚴(yán) 玥1, 嚴(yán) 實(shí)2, 江 赟1

(1.重慶工商大學(xué) 檢測(cè)控制集成系統(tǒng)重慶市市級(jí)工程實(shí)驗(yàn)室,重慶 400067;2.重慶川儀分析儀器有限公司，重慶 400060)

在分析影響火電廠氮氧化物濃度檢測(cè)精度的多種原因基礎(chǔ)上，提出建立RBF-BP神經(jīng)網(wǎng)絡(luò)模型進(jìn)行改進(jìn)的方法，并詳細(xì)說(shuō)明了樣本數(shù)據(jù)、神經(jīng)網(wǎng)絡(luò)構(gòu)成、訓(xùn)練函數(shù)、回歸因子、動(dòng)量因子等關(guān)鍵技術(shù)。以2015年實(shí)際火電廠檢測(cè)數(shù)據(jù)為依據(jù)，對(duì)建立的RBF-BP神經(jīng)網(wǎng)絡(luò)進(jìn)行了仿真實(shí)驗(yàn)，并與單一RBF仿真實(shí)驗(yàn)進(jìn)行對(duì)比分析，實(shí)驗(yàn)表明平均相對(duì)誤差為0.3 %～1 %，證明了該方法的有效性。

干擾; 濃度檢測(cè); 徑向基函數(shù); 反向傳播; 氮氧化物

0 引 言

目前火電廠尾氣排放在線監(jiān)測(cè)儀通常采用紅外線吸收法以及傅里葉吸收法原理。由于存在大氣溫度、樣氣溫度、燃煤質(zhì)量、電器干擾以及樣氣交叉等干擾因素的影響，表現(xiàn)出干擾因素多樣化、形成因素不穩(wěn)定的特點(diǎn)。采用硬件補(bǔ)償存在硬件漂移和調(diào)試?yán)щy等問(wèn)題，且也無(wú)法做到全額補(bǔ)償；而采用信息融合技術(shù)的軟件補(bǔ)償方式不需要調(diào)試，避免對(duì)各種因素的定性定量評(píng)估以及數(shù)學(xué)模型，補(bǔ)償效果好，因而逐漸興起，并在相關(guān)領(lǐng)域展現(xiàn)了較好的研究和應(yīng)用價(jià)值[1～3]。

從實(shí)際工程角度上，目前國(guó)內(nèi)大型火力發(fā)電廠通常使用單一傳感器，測(cè)量經(jīng)轉(zhuǎn)換后的NO濃度作為氮氧化物濃度測(cè)量值。一般監(jiān)測(cè)精度只能控制在1 %以?xún)?nèi)。通常分析儀器整重約200 kg，體積約0.96 m3,安裝周期需要大約一周。受到施工安置難度、成本以及實(shí)際工作環(huán)境相對(duì)惡劣等(排放煙囪高度為180～210 m,樣氣采集高度90 m)原因約束，很難再通過(guò)增加儀器硬件或者使用高精度儀器的方式提高實(shí)際測(cè)量精度。

本文不再以硬件方式改進(jìn)儀器最小量程，而是在現(xiàn)有氣體分析儀器元件敏感范圍內(nèi)，根據(jù)信息融合技術(shù)，采用反向傳播(back propagation,BP)和徑向基函數(shù)(radial-basis-function,RBF)神經(jīng)網(wǎng)絡(luò)結(jié)合的方法，以實(shí)際火電廠工業(yè)排放數(shù)據(jù)為訓(xùn)練樣本，設(shè)計(jì)出合理的神經(jīng)網(wǎng)絡(luò)，完成網(wǎng)絡(luò)訓(xùn)練，并模擬儀器實(shí)際工作環(huán)境設(shè)置測(cè)試數(shù)據(jù)進(jìn)行計(jì)算機(jī)仿真，將單一RBF網(wǎng)絡(luò)和RBF-BP網(wǎng)絡(luò)以相對(duì)誤差、殘差、迭代耗時(shí)等對(duì)比以分析其有效性。

1 干擾因素分析和解決方案提出

由于火電實(shí)際工藝流程中無(wú)法做到相對(duì)穩(wěn)定，對(duì)氮氧化物濃度檢測(cè)的干擾主要表現(xiàn)在其濃度輸出值(p)不僅只決定于一個(gè)目標(biāo)參量(吸收的紅外線能量e)，而是與非目標(biāo)參量如樣氣流量fr，粉塵含量d，溫度t，含水量w以及電器噪音f，各種有機(jī)物o(如：焦油、苯萘)等有關(guān)的多元函數(shù)，即

P=f(fr,d,t,w,n,o,…)

(1)

在非線性映射的研究過(guò)程中，充分利用傳感器在不同時(shí)間或者空間上的冗余互補(bǔ)信息，依據(jù)某種準(zhǔn)則來(lái)進(jìn)行分析、支配和使用，以獲得比當(dāng)前單一數(shù)據(jù)更為充分的信息是基于傳感數(shù)據(jù)融合的核心[4]。在神經(jīng)網(wǎng)絡(luò)、專(zhuān)家系統(tǒng)、遺傳算法和模糊系統(tǒng)等數(shù)據(jù)融合的研究領(lǐng)域中，神經(jīng)網(wǎng)絡(luò)的一致性逼近原則認(rèn)為：如果能夠確保足夠多的淺層網(wǎng)絡(luò)的隱層神經(jīng)元個(gè)數(shù)或者足夠深的深度網(wǎng)絡(luò)，就一定能找到一個(gè)淺層網(wǎng)絡(luò)和一個(gè)深層網(wǎng)絡(luò)已任意精度逼近任意一個(gè)非線性映射[5]。

由于地理環(huán)境及煤炭質(zhì)量等多種干擾因素并不具有橫向可比性，因此可以引入“時(shí)間參數(shù)”，使用同一火電廠的實(shí)時(shí)檢測(cè)數(shù)據(jù)為訓(xùn)練數(shù)據(jù)，尋找合適的網(wǎng)絡(luò)結(jié)構(gòu)并解決訓(xùn)練函數(shù)，動(dòng)量因子、隱藏層設(shè)計(jì)等關(guān)鍵問(wèn)題，以得到合理的神經(jīng)網(wǎng)絡(luò)，并用模擬現(xiàn)實(shí)環(huán)境的測(cè)試樣本進(jìn)行分析和評(píng)估。

2 關(guān)鍵設(shè)計(jì)說(shuō)明

2.1 氣體分析儀器說(shuō)明

以川儀P200—UV為分析檢測(cè)單元，分析單元結(jié)構(gòu)如圖1所示。儀器(FS,滿(mǎn)量程)情況下誤差相關(guān)指標(biāo)為：線性誤差不超過(guò)±2 %，零點(diǎn)、量程漂移不超過(guò)5 %，輸出波動(dòng)不大于1 %。

圖1 P200—UV儀器結(jié)構(gòu)圖Fig 1 Structure of P200—UV

2.2 訓(xùn)練過(guò)程

BP神經(jīng)網(wǎng)絡(luò)的自學(xué)習(xí)能力強(qiáng)，實(shí)現(xiàn)相對(duì)容易，已經(jīng)在識(shí)別與去噪、回歸、分類(lèi)等領(lǐng)域得到了相當(dāng)廣泛的運(yùn)用[6～8]。但其學(xué)習(xí)算法容易陷入“局部最優(yōu)”，且網(wǎng)絡(luò)初始值的選取、隱藏層層數(shù)的選取很大程序上依賴(lài)于經(jīng)驗(yàn)。而RBF具有較高的收斂速度并在儀器儀表的很多領(lǐng)域中都有很廣泛的應(yīng)用價(jià)值[9,10]。因此采用BP與RBF相結(jié)合的方法，以提高其泛化性能、克服非標(biāo)準(zhǔn)訓(xùn)練樣本集的錯(cuò)誤引導(dǎo)。

神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過(guò)程可簡(jiǎn)述為：訓(xùn)練樣本進(jìn)行RBF神經(jīng)網(wǎng)絡(luò)訓(xùn)練，得到一個(gè)相對(duì)穩(wěn)定合理的網(wǎng)絡(luò)結(jié)構(gòu)，并用測(cè)試數(shù)據(jù)進(jìn)行預(yù)估。再將同樣的訓(xùn)練樣本經(jīng)由已經(jīng)訓(xùn)練好的RBF網(wǎng)絡(luò)的結(jié)果作為BP網(wǎng)絡(luò)的輸入并對(duì)其進(jìn)行訓(xùn)練，最后以同樣的測(cè)試數(shù)據(jù)進(jìn)行仿真并進(jìn)行分析。

2.3 相關(guān)數(shù)據(jù)說(shuō)明

1)訓(xùn)練樣本

神經(jīng)網(wǎng)絡(luò)樣本輸入向量包括樣氣溫度、流量、含水量、粉塵顆粒以及儀器輸出電壓組成。以火電廠現(xiàn)場(chǎng)得到的80組真實(shí)數(shù)據(jù)作為訓(xùn)練樣本。在現(xiàn)場(chǎng)環(huán)境中經(jīng)由預(yù)處理裝置[11]能夠保證樣氣在進(jìn)入分析儀器前能夠達(dá)到的相關(guān)性能指標(biāo)為：溫度：3～5 ℃；粉塵顆粒：0～0.5 μm；流量40～60 L/H；水分含量<0.8(現(xiàn)場(chǎng)工藝環(huán)境無(wú)法做到零水分)。但由于現(xiàn)場(chǎng)未安裝相應(yīng)傳感器，并不能確定其精確值，故采用二維插值方法生成相關(guān)區(qū)域中的值，其Mathlab相關(guān)代碼如下

N=size(x,2);

[xx0,yy0]=meshgrid(1∶N,1∶6);

[xx1,yy1]=meshgrid(linspace(1,N,100),1∶6);

XX=interp2(xx0,yy0,x,xx1,yy1,’cubic’)

式中 x為輸入向量與對(duì)應(yīng)的目標(biāo)輸出合并而成的6×80的矩陣，XX為經(jīng)過(guò)插值得到的訓(xùn)練樣本。

2)測(cè)試樣本

在同等環(huán)境下20組標(biāo)準(zhǔn)濃度氣體，其中溫度、流量、含水量、粉塵顆粒仍然以二維插值方法在同等范圍內(nèi)生成作為測(cè)試用例。

3)神經(jīng)網(wǎng)絡(luò)輸出

輸出結(jié)果為濃度標(biāo)量，即網(wǎng)絡(luò)輸出層神經(jīng)元節(jié)點(diǎn)個(gè)數(shù)為1。

2.4 神經(jīng)網(wǎng)絡(luò)各層設(shè)計(jì)

1)樣本數(shù)據(jù)：由前所述，輸入向量定義為5×80的矩陣，目標(biāo)輸出值為1×80行向量。

2)RBF網(wǎng)絡(luò)設(shè)計(jì)

第一層為輸入層，兩個(gè)隱含層，以及一個(gè)輸出層，其中回歸因子的選擇遵循以下算法：

對(duì)于i≠j,m和wj正交，y(t)的能量為

(2)

移去其均值后，y為期望輸出的矢量，則y(t)的方差為

(3)

(4)

對(duì)于可選擇的若干個(gè)回歸因子，每個(gè)回歸因子對(duì)對(duì)應(yīng)一個(gè)誤差壓縮比，從誤差壓縮比中選擇最大的一個(gè)，這個(gè)誤差壓縮比對(duì)應(yīng)的回歸因子就是最終選擇的回歸因子。

不進(jìn)行閥值自適應(yīng)調(diào)整，使用均方誤差(mean square error，MSE)為誤差容限單位，當(dāng)訓(xùn)練誤差仍大于目標(biāo)誤差時(shí)，RBF將以逐漸增加神經(jīng)元使訓(xùn)練誤差減小，若達(dá)到最大神經(jīng)元個(gè)數(shù)極限仍為達(dá)到目標(biāo)誤差時(shí)，網(wǎng)絡(luò)不收斂；

3)BP網(wǎng)絡(luò)設(shè)計(jì)

增加隱含層可以明顯減少網(wǎng)絡(luò)規(guī)模，因此設(shè)定BP網(wǎng)絡(luò)隱含層為10，隱含層使用Sigmoid函數(shù)，在傳統(tǒng)LMS(Least Mean Square Algorithm，最小均方誤差算法)[12]的權(quán)值更新階段引入動(dòng)量因子α(0<α<1)，即

(5)

節(jié)點(diǎn)間傳遞函數(shù)為硬限幅傳遞函數(shù)hardlim；訓(xùn)練函數(shù)為動(dòng)量反傳和動(dòng)態(tài)自適應(yīng)學(xué)習(xí)率的梯度下降traingdx。使用MSE作為誤差容限單位。BP網(wǎng)絡(luò)將以增加迭代次數(shù)使訓(xùn)練誤差減小，若達(dá)到指定迭代次數(shù)極限還未達(dá)到指定誤差，則網(wǎng)絡(luò)不收斂，實(shí)驗(yàn)失敗。

3 實(shí) 驗(yàn)

用前述訓(xùn)練樣本輸入單一RBF網(wǎng)絡(luò)，設(shè)置誤差容限為1×10-8，最大神經(jīng)元個(gè)數(shù)極限為100，其訓(xùn)練過(guò)程及結(jié)果如圖2所示。此時(shí)RBF網(wǎng)絡(luò)誤差率呈不斷下降趨勢(shì)，迭代至72次后低于指定誤差訓(xùn)練結(jié)束，表現(xiàn)出了較好的收斂速度，此時(shí)RBF平均相對(duì)誤差為4.14 %，但最大相對(duì)誤差為19.8 %。由此可見(jiàn)單一RBF對(duì)訓(xùn)練樣本數(shù)據(jù)表現(xiàn)出了極不穩(wěn)定的性能。

圖2 誤差性能曲線Fig 2 Error performance curve

將上述 RBF網(wǎng)絡(luò)得到輸出作為BP網(wǎng)絡(luò)輸入，采用feedforwardnet函數(shù)創(chuàng)建BP網(wǎng)絡(luò)，設(shè)置誤差容限同樣為1×10-8，隱含層設(shè)置為10，迭代次數(shù)極限為1 000，進(jìn)行訓(xùn)練，RBF-BP網(wǎng)絡(luò)的迭代次數(shù)為13次，遠(yuǎn)低于指定迭代次數(shù)，達(dá)到訓(xùn)練目的完成訓(xùn)練。為了更加直觀反應(yīng)訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)測(cè)試數(shù)據(jù)的準(zhǔn)確性，圖4給出給出了單一RBF，RBF-BP組合神經(jīng)網(wǎng)絡(luò)的20組組測(cè)試數(shù)據(jù)的測(cè)試結(jié)果以及真實(shí)值。由圖3可知，單一RBF輸出對(duì)于部分測(cè)試樣本沒(méi)有明顯誤差，但在部分樣本上表現(xiàn)出了明顯差異，且該差異與氣體濃度值大小無(wú)明顯關(guān)系，很難加以?xún)?yōu)化，預(yù)測(cè)精度不如RBF-BP。

圖3 RBF和RBF-BP測(cè)試結(jié)果比較Fig 3 Comparison of the results of the RBF and RBF-BP

同時(shí)RBF-BP組合神經(jīng)預(yù)測(cè)結(jié)果的輸出向量與目標(biāo)向量的誤差區(qū)間相對(duì)較窄，表現(xiàn)出了較好的穩(wěn)定性，該特性可以通過(guò)RBF,RBF-BP的殘差對(duì)比得到，如圖4所示。此時(shí)，RBF-BP的平均相對(duì)誤差為0.265 7 %，最大相對(duì)誤差僅為1.00 %。優(yōu)于現(xiàn)有儀器精度要求。

圖4 RBF和RBF-BP測(cè)試殘差比較Fig 4 Comparison of the residual of the RBF and RBF-BP

4 結(jié) 論

本文從影響氮氧化物濃度測(cè)試精度的多種不確定因素出發(fā)，設(shè)計(jì)RBF-BP組合神經(jīng)網(wǎng)絡(luò)，該網(wǎng)絡(luò)結(jié)合單一RBF和BP神經(jīng)網(wǎng)絡(luò)的優(yōu)點(diǎn)，采用實(shí)際工程環(huán)境數(shù)據(jù)對(duì)其進(jìn)行網(wǎng)絡(luò)訓(xùn)練和仿真，以較快的速度和較高的精度提升氮氧化物濃度測(cè)試精度，測(cè)試數(shù)據(jù)顯示了其有效性。對(duì)該網(wǎng)絡(luò)的硬件化可用于改進(jìn)現(xiàn)有氮氧化物濃度分析儀的檢測(cè)精度以及實(shí)際工作環(huán)境中的分析儀器自動(dòng)校準(zhǔn)。

[1] 李玉軍,湯曉君,劉君華.粒子群優(yōu)化算法在改善傳感器穩(wěn)定性中的應(yīng)用[J].儀器儀表學(xué)報(bào)，2010，31(8):1756-1761.

[2] 嚴(yán)文娟,張 晶,胡廣芹,等.BP神經(jīng)網(wǎng)絡(luò)用于肝炎患者舌診近紅外光譜的研究[J].光譜學(xué)與光譜分析，2010,30(10):2628-2630.

[3] 陳曉峰,龍長(zhǎng)江,牛智有,等.基于潛在語(yǔ)義分析與NIR的中藥材分類(lèi)研究[J].光學(xué)學(xué)報(bào)，2014，34(9):288-293.[4] 周觀民,李榮會(huì).基于神經(jīng)網(wǎng)絡(luò)的傳感器網(wǎng)絡(luò)數(shù)據(jù)融合技術(shù)研究[J].計(jì)算機(jī)仿真,2011,28(10):118-120.

[5] 張 蕾,章 毅.大數(shù)據(jù)分析的無(wú)線神經(jīng)網(wǎng)絡(luò)方法[J].計(jì)算機(jī)研究與發(fā)展,2016,53(1):68-79.

[6]LarochelleH，LaulyS.Aneuralautoregressivetopicmodel[C]∥AdvancesinNeuralInformationProcessingSystems.Nevada，UnitedStates:MITPress，2012:2717-2725.

[7]TangY,SalakhutdinovR,HintonG.RobustBoltzmannMachinesforRecognitionandDenoising[C]∥ComputerVisionandPatternRecognition(CVPR),2012IEEEConferenceon.USA:IEEE,2012:2264-2271.

[8] 洪 磊,龔雪飛,孫壽通,等.Adaboost集成BP神經(jīng)網(wǎng)絡(luò)在傳感器陣列檢測(cè)系統(tǒng)中的應(yīng)用,傳感器與微系統(tǒng),2015,34(4):148-154.

[9]KevinHoAuthorVitae.TrainingRBFnetworktotoleratesinglenodefault[J].Neurocomputing,2011,74(6):1046-1052.

[10] 李德霞,黃廷磊,林 科,等.基于RBF神經(jīng)網(wǎng)絡(luò)算法在重金屬濃度測(cè)量中的應(yīng)用[J].儀器技術(shù)與傳感器,2015(10):85-87.

[11] 嚴(yán) 玥,嚴(yán) 實(shí),楊永斌,等.Adaboost集成BP神經(jīng)網(wǎng)絡(luò)在火電廠SO2濃度檢測(cè)中的應(yīng)用[J].傳感器與微系統(tǒng),2016,35(9):148-151.

[12] 陳 明.Mathlab神經(jīng)網(wǎng)絡(luò)原理與實(shí)例精解[M].北京:清華大學(xué)出版社,2013:156-172.

Application of RBF-BP neural network in the detection of nitrogen oxide concentration in thermal power plant*

YAN Yue1, YAN Shi2, JIANG Yun1

(1.Chongqing Engineering Laboratory for Detection,Control and Integrated System,Chongqing Technology and Business University,Chongqing 400067,China;2.Chongqing Chuanyi Analyzer Co Ltd,Chongqing 400060,China)

Aiming at detection accuracy of nitrogen oxide concentration may interfered by all kinds of factors,the method of establishing the RBF-BP network is proposed，and then the sample data,Neural network structure, training function,regression factor,momentum factor are described.According from A Large thermal power plant’s real numbers in 2015,the computer simulation and analysis show that this method can effectively improve the accuracy.The overall-average deviation was 0.3 %～1 %.

interferes; concentration detection; RBF;BP; nitrogen oxides

10.13873/J.1000—9787(2016)12—0149—03

TP 212.2； TN 911.72

B

1000—9787(2016)12—0149—03

嚴(yán) 玥(1976-)，女，滿(mǎn)，重慶人，碩士，講師，主要從事信號(hào)處理及運(yùn)用、傳感器技術(shù)及軟件方向研究工作。