汪明先， 吳建德,2， 王曉東,2

(1.昆明理工大學(xué) 信息工程與自動(dòng)化學(xué)院，云南 昆明 650500；2.云南省礦物管道輸送工程技術(shù)研究中心，云南 昆明 650500)

PSO-ELM的漿體管道臨界淤積流速預(yù)測(cè)模型研究*

汪明先1， 吳建德1,2， 王曉東1,2

(1.昆明理工大學(xué) 信息工程與自動(dòng)化學(xué)院，云南 昆明 650500；2.云南省礦物管道輸送工程技術(shù)研究中心，云南 昆明 650500)

針對(duì)漿體管道臨界淤積流速預(yù)測(cè)難度大、精度低等問(wèn)題，提出了粒子群優(yōu)化—極限學(xué)習(xí)機(jī)(PSO-ELM)的臨界淤積流速預(yù)測(cè)模型。利用PSO算法對(duì)ELM模型參數(shù)輸入權(quán)值和隱元偏置進(jìn)行優(yōu)化，應(yīng)用優(yōu)化得到的ELM模型對(duì)預(yù)測(cè)集進(jìn)行預(yù)測(cè)。通過(guò)實(shí)驗(yàn)仿真得到預(yù)測(cè)結(jié)果的最大誤差為5.73 %，預(yù)測(cè)效果優(yōu)于常規(guī)的ELM模型和反向傳播(BP)神經(jīng)網(wǎng)絡(luò)模型。

漿體管道; 臨界淤積流速; 粒子群優(yōu)化算法; 極限學(xué)習(xí)機(jī)

0 引 言

近年來(lái)，漿體管道輸送獲得迅猛發(fā)展，已經(jīng)成為第五大運(yùn)輸手段，是一種高效率、低成本、減少土地占用、無(wú)污染的固體物料輸送技術(shù)[1]。然而，能否確定某一實(shí)際輸送管道合適的臨界淤積流速值直接關(guān)系到其能否安全運(yùn)行，并在相同條件下能否取得最大的經(jīng)濟(jì)效益[2]。目前應(yīng)用較多的漿體管道臨界淤積流速預(yù)測(cè)方法有經(jīng)驗(yàn)半經(jīng)驗(yàn)公式法[3]、反向傳播(back propagation,BP)神經(jīng)網(wǎng)絡(luò)法[4]和支持向量機(jī)(SVM)預(yù)測(cè)模型[5]，在某一程度上來(lái)說(shuō)也獲得了良好的效果，但其均有無(wú)法避免的缺陷，如經(jīng)驗(yàn)公式往往考慮的因素比較單一，只適用于某些具體條件，由此導(dǎo)致計(jì)算出來(lái)的結(jié)果與實(shí)際所測(cè)得的數(shù)據(jù)之間有很大的誤差；神經(jīng)網(wǎng)絡(luò)所需的數(shù)據(jù)量較大，經(jīng)常會(huì)遇到過(guò)學(xué)習(xí)、陷入局部極值等問(wèn)題[6]，在實(shí)際應(yīng)用中為解決該問(wèn)題，文獻(xiàn)[7～9]提出了粒子群優(yōu)化BP(particle swarm optimization,PSO-BP)算法、改進(jìn)的BP算法用于研究多元有害氣體的檢測(cè)及傳感器環(huán)境溫度補(bǔ)償問(wèn)題，但其無(wú)法解決訓(xùn)練速度慢、預(yù)測(cè)效率一般等問(wèn)題； SVM預(yù)測(cè)模型的參數(shù)選擇較為復(fù)雜，從而嚴(yán)重影響了它的預(yù)測(cè)效率，文獻(xiàn)[10]用人工魚(yú)群算法優(yōu)化SVM的懲罰參數(shù)和核函數(shù)參數(shù)，在一定程度上避免了參數(shù)選擇的盲目性，提高了模型的識(shí)別精度，但仍然無(wú)法從根本上解決該問(wèn)題。此時(shí)，極限學(xué)習(xí)機(jī)(extreme learning machine,ELM)預(yù)測(cè)模型[11]應(yīng)運(yùn)而生。

本文在前人實(shí)驗(yàn)研究的基礎(chǔ)上，通過(guò)理論分析研究，確定了管道直徑、粒料平均直徑、物料密度與液體密度比和漿體體積濃度為影響漿體管道臨界淤積流速的主要因素，建立將其作為輸入因子，臨界淤積流速值為輸出因子的ELM預(yù)測(cè)模型，并針對(duì)ELM的不足，用粒子群算法對(duì)ELM中網(wǎng)絡(luò)的輸入層權(quán)值和隱單元閾值進(jìn)行優(yōu)選，進(jìn)而提出了PSO-ELM模型對(duì)大直徑漿體輸送管道的臨界淤積流速進(jìn)行了研究。

1 材料與方法

1.1 臨界淤積流速影響因素分析

從前人的實(shí)驗(yàn)中可以發(fā)現(xiàn)影響漿體管道臨界流速的因素非常多，其中主要包括顆粒的大小、密度、形狀、粒度分布、漿體濃度、管道運(yùn)輸直徑、外部溫度等等[12]。Duckworth 等人采用比水輕和比水重的物質(zhì)分別進(jìn)行了實(shí)驗(yàn)，觀察到的實(shí)驗(yàn)現(xiàn)象為：在固體顆粒比水重的情況下，臨界流速隨著固體密度的加大而增大；在水比固體顆粒重的情況下，臨界流速隨著固體密度的減小而加快。通過(guò)濃度和臨界流速關(guān)系Kazanskij公式可以總結(jié)出這樣的結(jié)論：對(duì)于某一具體管徑來(lái)說(shuō)，通常情況下，漿體輸送管道的臨界淤積流速與漿體濃度呈正相關(guān)關(guān)系時(shí)，其漿體濃度處于低濃度范圍；反之，但當(dāng)濃度超過(guò)一定范圍時(shí)，呈負(fù)相關(guān)關(guān)系，且易形成層流狀態(tài)，阻力也快速增大。杜蘭德公式[13]、Wasp公式[14]中均表明臨界淤積流速與管徑D的1/3方成正比，汪東等人根據(jù)試驗(yàn)數(shù)據(jù)繪制出臨界流速與管徑的次方關(guān)系圖提出，在其他條件相同時(shí)，臨界流速正比于管徑的 1/3 次方[15]。漿體中顆粒在合理比例時(shí)，細(xì)顆粒間容易形成絮網(wǎng)結(jié)構(gòu)，在阻尼作用的影響下，大大降低了粗顆粒的沉降速度，且不易形成滑、跳移狀態(tài)，減少阻力損失。

基于以上分析可以看出，管道直徑D、粒料平均直徑d、物料密度與液體密度的比值S=ρs/ρ和漿體體積濃度Cv為影響淤積臨界流速的主要因素。因此，漿體在水力固體運(yùn)輸管道中淤積臨界流速可以表示為

(1)

1.2 基于PSO-ELM的預(yù)測(cè)建模

本文所提出的PSO-ELM算法，其基本思想是運(yùn)用粒子群優(yōu)化對(duì)單隱層前饋網(wǎng)絡(luò)的輸入層權(quán)值和隱元偏置進(jìn)行優(yōu)選，進(jìn)而將其優(yōu)選所得的輸入權(quán)值和隱元偏置應(yīng)用于漿體管道臨界淤積流速的預(yù)測(cè)。算法的基本步驟如下：

1)隨機(jī)產(chǎn)生種群：設(shè)定粒子的群體規(guī)模m，種群的個(gè)體維度D，該個(gè)體維度取決于ELM模型的輸入?yún)?shù)，學(xué)習(xí)因子的初始值c1ini和c2ini，最大迭代次數(shù)kmax，實(shí)驗(yàn)結(jié)束條件最小誤差。

4)粒子速度更新與位置更新：速度調(diào)整規(guī)則如下

(2)

5)全局求解最佳適應(yīng)度：不斷比較各粒子當(dāng)前的適應(yīng)度f(wàn)(xi)與其歷史最優(yōu)適應(yīng)度f(wàn)(Pbest)，并結(jié)合結(jié)束條件確定其值大小。將結(jié)束的條件設(shè)置為尋優(yōu)達(dá)到最大迭代次數(shù)kmax或評(píng)價(jià)值小于指定的精度,該精度設(shè)為0.001。

6)將最優(yōu)解ai和bi代入ELM模型。

7)用測(cè)試數(shù)據(jù)對(duì)ELM進(jìn)行訓(xùn)練，然后求解。

1.3 數(shù)據(jù)分析方法

采用均方根誤差(RMSE)、決定系數(shù)(R2)、最大相對(duì)誤差、平均相對(duì)誤差等指標(biāo)評(píng)價(jià)各模型計(jì)算結(jié)果與標(biāo)準(zhǔn)值之間的誤差以及一致性，計(jì)算公式如下

(3)

(4)

2 實(shí)驗(yàn)與分析

2.1 實(shí)驗(yàn)數(shù)據(jù)

結(jié)合國(guó)內(nèi)外有關(guān)鐵精礦、煤漿、硫精礦、磷精礦、電廠灰渣、尾礦等漿體管道整理得到的82組數(shù)據(jù)[17]，現(xiàn)將其中的77個(gè)樣本構(gòu)成訓(xùn)練集，用來(lái)進(jìn)行ELM模型的訓(xùn)練，考慮到樣本總體組成比較復(fù)雜和磷精礦管道實(shí)驗(yàn)數(shù)據(jù)過(guò)少等問(wèn)題，故在測(cè)試數(shù)據(jù)的選取時(shí)，不考慮磷精礦，隨機(jī)在其它 5類(lèi)漿體管道中抽取5個(gè)樣本作為測(cè)試數(shù)據(jù)，建立以管道直徑、物料平均粒徑、漿體體積濃度、物料比為輸入因子，輸出變量為臨界淤積流速值 的ELM模型，對(duì)漿體管道臨界淤積流速進(jìn)行預(yù)測(cè)。

2.2 預(yù)測(cè)結(jié)果分析

假定隱含層節(jié)點(diǎn)數(shù) ,在確定了ELM模型中的隱含層節(jié)點(diǎn)數(shù)后，采用PSO對(duì)模型隱含層輸入權(quán)值ai和偏置向量bi進(jìn)行優(yōu)化選擇，以確定最優(yōu)的ELM模型。PSO的各參數(shù)如下：最大迭代次數(shù)kmax設(shè)為100，粒子數(shù)m設(shè)為20。最大慣性權(quán)重wmax取0.9；最小慣性權(quán)重wmin取0.4。學(xué)習(xí)因子c1初始值取1.68，c2初始值取1.3，學(xué)習(xí)因子c1,c2的最終迭代值c1fin和c2fin均取0.5。經(jīng)PSO算法尋優(yōu)后，得到最優(yōu)適應(yīng)度f(wàn)(gbesti)=0.009 335,其對(duì)應(yīng)的ELM參數(shù)ai,bi值即為最優(yōu)參數(shù)，此時(shí)

a=

圖1 適應(yīng)度隨迭代次數(shù)的變化曲線

根據(jù)優(yōu)化得到的PSO-ELM模型對(duì)預(yù)測(cè)集進(jìn)行預(yù)測(cè)，其預(yù)測(cè)結(jié)果及真實(shí)值與預(yù)測(cè)值之間的差值圖如圖2、圖3所示。

從圖2、圖3可知，PSO-ELM模型的預(yù)測(cè)結(jié)果與真實(shí)值擬合效果較好，最大差值為0.085 2，且各差值在區(qū)間[-0.1 0.1]內(nèi)波動(dòng)，這說(shuō)明PSO-ELM模型的預(yù)測(cè)值的誤差是隨機(jī)性的，不存在系統(tǒng)誤差，也從側(cè)面證明，本模型所選擇的輸入、輸出參數(shù)是正確可行的?，F(xiàn)將該模型與ELM模型、BP神經(jīng)網(wǎng)絡(luò)模型進(jìn)行對(duì)比研究，利用上文中所確定的訓(xùn)練集、驗(yàn)證集對(duì)其進(jìn)行訓(xùn)練，并產(chǎn)生預(yù)測(cè)結(jié)果，對(duì)比結(jié)果見(jiàn)圖4及表1。

圖2 真實(shí)值與預(yù)測(cè)值曲線圖

圖3 真實(shí)值與預(yù)測(cè)值之間的差值圖

圖4 三種模型的預(yù)測(cè)結(jié)果對(duì)比

由表1可以看出，采用PSO-ELM方法的5個(gè)樣本中，相對(duì)誤差最大值僅為5.73 %，小于ELM的最大相對(duì)誤差7.69 %。與此同時(shí)，分析二者的誤差可知，采用PSO-ELM方法得到的預(yù)測(cè)結(jié)果的相對(duì)誤差與ELM相比均有相當(dāng)程度上的減少，這說(shuō)明運(yùn)用PSO算法對(duì)優(yōu)化ELM的輸入權(quán)值和隱元偏置進(jìn)行優(yōu)選這一方法是可行的，所優(yōu)選得到的輸入權(quán)值和隱元偏置能夠提高網(wǎng)絡(luò)的泛化能力，進(jìn)而能夠提高該模型的預(yù)測(cè)精度。

表1 臨界淤積流速預(yù)測(cè)結(jié)果

2.3 模型精度比較分析

為顯示本文所建預(yù)測(cè)模型的優(yōu)勢(shì)，定義RMSE、R2、最大相對(duì)誤差、平均相對(duì)誤差作為評(píng)價(jià)指標(biāo)評(píng)價(jià)PSO-ELM模型、ELM模型、BP神經(jīng)網(wǎng)絡(luò)模型的計(jì)算結(jié)果的精度如表2所示。

從表2可以看出，BP神經(jīng)網(wǎng)絡(luò)模型由于樣本數(shù)據(jù)少，臨界淤積流速預(yù)測(cè)值誤差大，最大誤差可達(dá)11.8 %；PSO-ELM預(yù)測(cè)模型較常規(guī)的ELM預(yù)測(cè)模型預(yù)測(cè)精度更高，從模型決定系數(shù)R2分析可知，該模型吻合度更好，但在尋優(yōu)過(guò)程中消耗了更多的時(shí)間?？梢?jiàn)，本文提出的PSO-ELM模型方法在漿體輸送管道臨界淤積流速的預(yù)測(cè)方面較常規(guī)ELM模型有著更高的預(yù)測(cè)精度。

表2 幾種模型預(yù)測(cè)結(jié)果評(píng)價(jià)指標(biāo)比較

3 結(jié)束語(yǔ)

1)根據(jù)ELM原理，確定以管徑、物料平均粒徑、漿體體積濃度、物料比為輸入因子，臨界淤積流速值為輸出因子，結(jié)合PSO，建立了PSO-ELM預(yù)測(cè)模型，避免了參數(shù)選擇的隨機(jī)性對(duì)模型精度的削弱。預(yù)測(cè)結(jié)果表明：PSO-ELM模型精度高，且誤差在6 %以?xún)?nèi)。

2)利用BP神經(jīng)網(wǎng)絡(luò)和常規(guī)的ELM模型對(duì)結(jié)果進(jìn)行預(yù)測(cè)，比較三者之間的預(yù)測(cè)效果，比較結(jié)果表明：PSO-ELM模型的預(yù)測(cè)精度要優(yōu)于ELM模型和BP神經(jīng)網(wǎng)絡(luò)模型，PSO-ELM模型被證明是一種可行的臨界淤積流速預(yù)測(cè)方法。

3)為進(jìn)一步提高粒子群算法的全局搜索能力，減少無(wú)效迭代次數(shù)，以更快的速度求得最優(yōu)適應(yīng)度，擬對(duì)影響粒子群搜索性能的慣性權(quán)重優(yōu)化進(jìn)行深入研究，以期提出改進(jìn)的PSO-ELM臨界淤積流速預(yù)測(cè)模型。

4)利用粒子群算法對(duì)極限學(xué)習(xí)機(jī)的輸入權(quán)值和隱元偏置的優(yōu)選過(guò)程是一個(gè)隨機(jī)過(guò)程，需要花費(fèi)一定的時(shí)間，在后階段將具體考慮如何用較少的時(shí)間來(lái)實(shí)現(xiàn)該優(yōu)選過(guò)程。

[1] 費(fèi)祥俊.漿體的物理特性與管道輸送流速[J].管道技術(shù)與設(shè)備,2000(1):1-4.

[2] 費(fèi)祥俊.漿體輸送管道的不淤流速研究[J].煤炭學(xué)報(bào),1997,22(5):532-536.

[3] Durand R.Basic relationships of the transportation of solids in pipes experiment[C]∥Research Proc Intern Assoc Hydr Res the 5th Congr,Minneapolis,1953:89-93.

[4] 韓 旭,張奇志,佟慶理,等.用人工神經(jīng)網(wǎng)絡(luò)法研究大直徑漿體輸送管道的淤積臨界流速[J].中國(guó)有色金屬學(xué)報(bào),1997,7(4):26-29.

[5] 李 杰,楚 恒,朱維樂(lè),等.基于支持向量機(jī)和遺傳算法的紋理識(shí)別[J].四川大學(xué)學(xué)報(bào),2005,37(4):104-108.

[6] 高海兵,高 亮,周 馳,等.基于粒子群優(yōu)化的神經(jīng)網(wǎng)絡(luò)訓(xùn)練算法研究[J].電子學(xué)報(bào),2004,32(9):1572-1574.

[7] 龔雪飛,徐 景,孫壽通,等.PSO-BP神經(jīng)網(wǎng)絡(luò)在多元有害氣體檢測(cè)中的應(yīng)用[J].傳感器與微系統(tǒng),2015,34(2):154-156.

[8] 趙 斌,王建華.基于PSO-BP的紅外溫度傳感器環(huán)境溫度補(bǔ)償[J].傳感器與微系統(tǒng),2015,34(2):47-49.

[9] 康國(guó)煉,楊遂軍,葉樹(shù)亮.改進(jìn)BP算法在熱流傳感器溫度補(bǔ)償中的應(yīng)用[J].傳感器與微系統(tǒng),2016,35(2):154-156.

[10] 洪 潔,王 璐,汪 超,等.基于人工魚(yú)群算法優(yōu)化SVM的手部動(dòng)作sEMG識(shí)別[J].傳感器與微系統(tǒng),2016,35(2):23-25.

[11] Huang Guangbin,Zhu Qinyu,Siew Chee Kheong.Extreme learning machine:Theory and applications[J].Neurocomputing,2006,70(1/3):489-501.

[12] 汪 東,許振良,孟慶華.漿體管道輸送臨界流速的影響因素及計(jì)算分析[J].管道技術(shù)與設(shè)備,2004(6):1-2.

[13] Durand R.The hydraulic transportation of coal and solid materials in pipes[C]∥Colloq of National Coal Board,London,1952：39-52．

[14] Wasp E J,Kenny J P,Gandhi R L.Solid-liquid flow-Slurry pipeline transition[M].Rockport,MA:Trans Tech Publ,1977.

[15] 汪 東.非均質(zhì)流水平管道臨界流速研究[J].礦冶工程,2008(4):20-22.

[16] 王小川,史 峰,郁 磊,等.Matlab神經(jīng)網(wǎng)絡(luò)43個(gè)案例分析[M].北京:北京航空航天大學(xué)出版社,2013:243-255.

[17] 丁宏達(dá).物料粒徑對(duì)漿體管道輸送臨界流速的影響[J].水力采煤與管道運(yùn)輸,2011(4):1-6.

Study on prediction model of critical deposition velocity in slurry pipeline based on PSO-ELM*

WANG Ming-xian1, WU Jian-de1,2, WANG Xiao-dong1,2

(1.Faculty of Information Engineering and Automation,Kunming University of Science and Technology,Kunming 650500,China; 2.Engineering Research Center for Mineral Pipeline Transportation of Yunnan,Kunming 650500,China)

Aiming at problems of great difficulties and low precision in predicting critical deposition velocity of slurry pipeline,a prediction model for critical deposition velocity based on particle swarm optimization and extreme learning machine(PSO-ELM)is proposed.The method using PSO algorithm to optimize the ELM model parameters such as the input weight and hidden element bias，the optimized ELM model is used to predict prediction set.The maximum error of predicted result is 5.73 % by experimental simulation,which shows that the prediction effect is superior to the conventional ELM and back propagation(BP)neural network model.

slurry pipeline; critical deposition velocity; particle swarm optimization(PSO)algorithm; extreme learning machine(ELM)

10.13873/J.1000—9787(2017)03—0066—04

2016—04—13

國(guó)家自然科學(xué)基金資助項(xiàng)目(51169007)；云南省科技計(jì)劃項(xiàng)目(2013DH034)；云南省中青年學(xué)術(shù)和技術(shù)帶頭人后備人才培養(yǎng)計(jì)劃項(xiàng)目(2011CI017)

TP 391; X 936

A

1000—9787(2017)03—0066—04

汪明先(1990-)，男，碩士研究生，研究方向?yàn)殚L(zhǎng)距離漿體管道輸送控制與智能控制算法研究。