張立峰, 肖 凱, 華回春
(1. 華北電力大學(xué) 自動(dòng)化系, 河北 保定 071003; 2. 華北電力大學(xué) 數(shù)理系, 河北 保定 071003)
隨著我國(guó)工業(yè)的不斷發(fā)展,對(duì)兩相流監(jiān)測(cè)日趨重要。兩相流系統(tǒng)在電力、化工、石油等工業(yè)部門廣泛存在,兩相流流型辨識(shí)也成為兩相流在線監(jiān)測(cè)的基礎(chǔ),是兩相流測(cè)量的重要研究方向[1]。兩相流各相間存在速度差且邊界尚不明確,其研究難度高于單向流[2]。兩相流流型變化直接決定相含率、流量及混合流速等參數(shù)變化,因此需對(duì)兩相流的流型進(jìn)行實(shí)時(shí)監(jiān)測(cè)[3]。電阻層析成像(electrical resistance tomography,ERT)技術(shù)是一種兩相流參數(shù)實(shí)時(shí)監(jiān)測(cè)技術(shù),其物理基礎(chǔ)是:根據(jù)不同介質(zhì)的電導(dǎo)率特性,通過測(cè)量敏感場(chǎng)邊界信號(hào),反推出內(nèi)部的電導(dǎo)率分布,通過電導(dǎo)率圖像判斷出物質(zhì)內(nèi)部的介質(zhì)分布狀況[4]。該方法具有非侵入、結(jié)構(gòu)簡(jiǎn)單及實(shí)時(shí)性高等特點(diǎn)。流型辨識(shí)則是基于ERT來滿足區(qū)分不同流型要求的一種技術(shù),其原理是依據(jù)ERT傳感器內(nèi)不同流型的電導(dǎo)率分布不同,通過測(cè)量ERT傳感器電極上的電壓得到管道內(nèi)部流型的特征,根據(jù)獨(dú)立邊界電壓測(cè)量值進(jìn)行圖像重建然后進(jìn)行流型辨識(shí),或直接根據(jù)獨(dú)立邊界電壓測(cè)量值進(jìn)行流型辨識(shí)。申超群等將ERT測(cè)量數(shù)據(jù)用模糊聚類與神經(jīng)網(wǎng)絡(luò)相結(jié)合的方法對(duì)4種流型進(jìn)行辨識(shí)[5];陳德運(yùn)等采用主成分分析方法對(duì)ERT系統(tǒng)中的邊界測(cè)量電壓數(shù)據(jù)進(jìn)行特征提取,然后作為支持向量機(jī)的輸入,對(duì)油水兩相流的均勻流、層流、環(huán)狀流、核心流四種流型進(jìn)行流型辨識(shí)[6];朱波則在ERT特征提取上采用了小波分析方法提取特征然后與神經(jīng)網(wǎng)絡(luò)進(jìn)行結(jié)合對(duì)流型進(jìn)行辨識(shí),取得了較好的效果[7];許燕斌等用獨(dú)立成分分析方法對(duì)氣液兩相流的ERT數(shù)據(jù)進(jìn)行了特征提取,并通過實(shí)驗(yàn)驗(yàn)證了其有效性[8]。
AdaBoost(Adaptive Boosting)算法基本原理是以上一次弱分類器的分類結(jié)果為依據(jù)來對(duì)下一個(gè)弱分類器調(diào)整樣本權(quán)重和弱分類器權(quán)重,最后將多個(gè)弱分類器組合為一個(gè)強(qiáng)分類器以此來達(dá)到提升單個(gè)弱分類器分類效果的目的[9,10]。由于AdaBoost具有可以將不同的分類算法作為弱分類器且能有效改善過擬合問題等優(yōu)點(diǎn)被廣泛應(yīng)用于各類研究之中。本文將以一維卷積神經(jīng)網(wǎng)絡(luò)作為弱分類器的AdaBoost算法結(jié)合形成1D-CNN-AdaBoost算法并用實(shí)驗(yàn)測(cè)得的ERT數(shù)據(jù)對(duì)泡狀流、彈狀流、段塞流及混狀流4種流型進(jìn)行流型辨識(shí)。
電阻層析成像系統(tǒng)主要是由傳感器、數(shù)據(jù)采集系統(tǒng)及圖像重建計(jì)算機(jī)構(gòu)成。其原理如圖1所示[11]。
圖1 電阻層析成像系統(tǒng)原理圖Fig.1 Schematic diagram of electrical resistance tomography system
由計(jì)算機(jī)發(fā)出控制信號(hào)經(jīng)采集系統(tǒng)對(duì)ERT傳感器的k個(gè)電極采用相鄰激勵(lì)模式進(jìn)行激勵(lì),并用相鄰測(cè)量模式采集出M個(gè)獨(dú)立邊界電壓測(cè)量值[12]。
(1)
本文采用16電極ERT傳感器,可得120個(gè)獨(dú)立測(cè)量值。將測(cè)量數(shù)據(jù)傳回計(jì)算機(jī),采用圖像重建算法即可獲得管道截面電導(dǎo)率分布的重建圖像[13]。
一維卷積神經(jīng)網(wǎng)絡(luò)(one-dimensional convolutional neural network,1D-CNN)是一種前饋神經(jīng)網(wǎng)絡(luò)。該網(wǎng)絡(luò)的基本工作原理是根據(jù)所選取的卷積核在輸入樣本上以特定步長(zhǎng)滑動(dòng)來選擇樣本的部分特征值,然后將其組合起來,形成更高級(jí)的特征,以此來進(jìn)行預(yù)測(cè)[14]。
本文使用1D-CNN結(jié)構(gòu)圖如圖2所示。1D-CNN中的網(wǎng)絡(luò)結(jié)構(gòu)經(jīng)過多次調(diào)整最終使用7層Conv1D來提取特征值,每2層Conv1D后添加1層池化層(MaxPooling1D)來保留主要特征,減少計(jì)算量。每層卷積層使用雙曲正切函數(shù)tanh作為卷積神經(jīng)網(wǎng)絡(luò)的激活函數(shù)來提高神經(jīng)網(wǎng)絡(luò)對(duì)模型的表達(dá)能力,其中tanh函數(shù)可表示為:
(2)
圖2 一維卷積神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.2 1D-CNN structure diagram
1D-CNN-AdaBoost算法是將1D-CNN作為Ada-Boost算法的弱分類器的一種算法,是一種自適應(yīng)增強(qiáng)算法,該算法由多個(gè)1D-CNN分類器組成,核心思想是依次訓(xùn)練多個(gè)1D-CNN,根據(jù)前1個(gè)1D-CNN分類器的結(jié)果來調(diào)整下1個(gè)1D-CNN分類器的樣本權(quán)值,最終根據(jù)每個(gè)1D-CNN分類器的分類效果決定其在強(qiáng)分類器中的權(quán)重。其算法流程圖如圖3所示。
圖3 1D-CNN-AdaBoost算法流程圖Fig.3 1D-CNN-AdaBoost Algorithm flowchart
1D-CNN-AdaBoost算法計(jì)算過程如下:
1) 輸入訓(xùn)練樣本;
2) 初始化權(quán)值分布Di;
(3)
式中:wi對(duì)應(yīng)第i個(gè)樣本的權(quán)值;N代表樣本數(shù)量。
3) 訓(xùn)練弱分類器,共進(jìn)行T次迭代t=1,…,T;
任意選擇一個(gè)卷積神經(jīng)網(wǎng)絡(luò),作為第t個(gè)弱分類器Ht并計(jì)算該弱分類器在權(quán)值分布Dt上的分類誤差et為:
(4)
式中:wt(i)對(duì)應(yīng)第t次迭代中第i個(gè)被錯(cuò)誤分類樣本的權(quán)值。
計(jì)算當(dāng)前卷積神經(jīng)網(wǎng)絡(luò)在最終強(qiáng)分類器中的權(quán)重αt:
(5)
更新下一次迭代的樣本權(quán)值Dt+1:
(6)
4) 將各個(gè)卷積神經(jīng)網(wǎng)絡(luò)弱分類器Ht及其權(quán)重αt組合為1個(gè)強(qiáng)分類器Hfinal;
(7)
通過符號(hào)函數(shù)sign的作用,得到最終的強(qiáng)分類器:
(8)
5) 將測(cè)試樣本輸入強(qiáng)分類器并輸出分類結(jié)果。
使用華北電力大學(xué)先進(jìn)測(cè)量實(shí)驗(yàn)室的可移動(dòng)氣水兩相流實(shí)驗(yàn)裝置進(jìn)行實(shí)驗(yàn)研究,裝置如圖4所示。
圖4 實(shí)驗(yàn)裝置示意圖Fig.4 Schematic diagram of experimental setup
裝置主要由水路循環(huán)、氣路循環(huán)及電阻層析成像3部分組成。儲(chǔ)氣罐里的氣體經(jīng)過壓力表、氣體渦輪流量計(jì)、壓力表以及浮子流量計(jì)到達(dá)ERT傳感器, 同時(shí)水箱里的水由水泵抽取經(jīng)過電磁流量計(jì)、閥門到達(dá)ERT,與氣體混合形成氣水兩相流;在水流量不變的情況下,通過調(diào)節(jié)氣相路的閥門開度進(jìn)而改變氣量可以得到不同的氣水兩相流流型;再由ERT將不同流型的信息進(jìn)行采集。
本次實(shí)驗(yàn)共采集泡狀流、彈狀流、段塞流、混狀流各1 500個(gè)樣本,共計(jì)6 000個(gè)樣本。每個(gè)樣本由120個(gè)獨(dú)立測(cè)量電壓值組成,構(gòu)成120個(gè)特征。泡狀流、彈狀流、段塞流、混狀流分別如圖5所示。
圖5 實(shí)驗(yàn)流型Fig.5 Experimental flow pattern
華北電力大學(xué)先進(jìn)測(cè)量實(shí)驗(yàn)室的可移動(dòng)氣水兩相流實(shí)驗(yàn)裝置實(shí)物圖如圖6所示。
圖6 實(shí)驗(yàn)裝置實(shí)物圖Fig.6 Physical image of experimental device
實(shí)驗(yàn)獲得的泡狀流、彈狀流、段塞流、混狀流各1 500個(gè)樣本共計(jì)6 000個(gè)樣本進(jìn)行訓(xùn)練集和測(cè)試集劃分,從6 000個(gè)樣本中隨機(jī)抽取1 200個(gè)數(shù)據(jù)作為測(cè)試集,剩余4 800個(gè)數(shù)據(jù)作為訓(xùn)練集。然后將訓(xùn)練樣本以每批次100個(gè)樣本作為一維卷積神經(jīng)網(wǎng)絡(luò)的輸入。損失函數(shù)選擇交叉熵?fù)p失函數(shù),學(xué)習(xí)率設(shè)置為0.2,并將AdaBoost的弱分類器個(gè)數(shù)設(shè)置為5。
圖7 1D-CNN流型辨識(shí)結(jié)果Fig.7 1D-CNN flow pattern recognition results
為了驗(yàn)證AdaBoost算法對(duì)1D-CNN的流型辨識(shí)率的提升效果,本文對(duì)比了1D-CNN和1D-CNN-AdaBoost對(duì)流型辨識(shí)的結(jié)果。其中1D-CNN對(duì)4種流型的辨識(shí)結(jié)果的混淆矩陣如圖7所示,圖7中分別對(duì)應(yīng)泡狀流、彈狀流、段塞流及混狀流,右側(cè)顏色條由0.0到1.0分別對(duì)應(yīng)著分類精度,顏色越深即分類精度越高。由圖7可見單純使用1D-CNN對(duì)段塞流的辨識(shí)效果最佳,其辨識(shí)正確率可以達(dá)到100%,對(duì)于彈狀流和混狀流的辨識(shí)正確率也分別能達(dá)到90%和91%,但是對(duì)于泡狀流的辨識(shí)正確率僅為84%。平均辨識(shí)正確率為91.2%。
采用1D-CNN-AdaBoost進(jìn)行4種流型辨識(shí)結(jié)果的混淆矩陣如圖8所示。由圖8可見,將1D-CNN與AdaBoost結(jié)合后的1D-CNN-AdaBoost算法對(duì)泡狀流、彈狀流、混狀流的辨識(shí)正確率在原有1D-CNN算法上均有了明顯的提高,尤其對(duì)泡狀流的辨識(shí)正確率由原來的84%提高到了95%。平均辨識(shí)率也由原來的91.2%提升到了97%,驗(yàn)證了AdaBoost算法對(duì)1D-CNN的流型辨識(shí)率的提升效果。
圖8 1D-CNN-AdaBoost流型辨識(shí)結(jié)果Fig.8 1D-CNN-AdaBoost flow pattern recognition results
為了探究1D-CNN-AdaBoost與其他流型辨識(shí)算法的流型辨識(shí)效果,本文對(duì)比了該算法與BP、支持向量機(jī)(support vector machines, SVM)以及決策樹3種算法的型辨識(shí)率。結(jié)果如表1所示。由表1可見1D-CNN-AdaBoost算法辨識(shí)精度優(yōu)于BP、SVM及決策樹算法,對(duì)泡狀流、彈狀流、段塞流、混狀流4種流型的辨識(shí)正確率均能達(dá)到95%以上,平均辨識(shí)精度可達(dá)到97%。
表1 各算法流型辨識(shí)率Tab.1 Flow pattern recognition rate of each algorithm (%)
本文研究了AdaBoost算法對(duì)1D-CNN的流型辨識(shí)率的提升效果,并通過實(shí)驗(yàn)驗(yàn)證了其有效性。最終選擇了1D-CNN作為AdaBoost算法的弱分類器,形成了1D-CNN-AdaBoost算法。然后用實(shí)驗(yàn)選取4種流型共6 000個(gè)電阻層析成像邊界電壓測(cè)量值樣本,采用1D-CNN-AdaBoost算法進(jìn)行流型辨識(shí),并與BP、SVM及決策樹流型辨識(shí)算法進(jìn)行比較。實(shí)驗(yàn)結(jié)果表明,采用1D-CNN-AdaBoost算法對(duì)泡狀流、彈狀流、段塞流及混狀流的辨識(shí)正確率均高于其它算法,平均辨識(shí)正確率達(dá)到了97%,為流型辨識(shí)提供了一種新的方法。