汪醒鵬

（湖南省水利水電科學(xué)研究所 長(zhǎng)沙市 410007）

范 哲

（益陽(yáng)市柘溪水電站 益陽(yáng)市 413500）

水輪發(fā)電機(jī)組承擔(dān)著將水能轉(zhuǎn)化為旋轉(zhuǎn)機(jī)械能，再將旋轉(zhuǎn)機(jī)械能轉(zhuǎn)化為電能的任務(wù)。與風(fēng)能、水能、熱能相比，電能對(duì)存儲(chǔ)運(yùn)送介質(zhì)有特殊要求，電能不能大量存儲(chǔ)，因此電能的生產(chǎn)和消費(fèi)必須同時(shí)進(jìn)行。GB/T 15945-2008[1]中規(guī)定：電力系統(tǒng)正常頻率范圍為50±0.2 Hz，頻率偏差過(guò)大不僅直接影響發(fā)電機(jī)組自身穩(wěn)定運(yùn)行，造成電網(wǎng)解列，而且導(dǎo)致以電動(dòng)機(jī)為動(dòng)力設(shè)備的生產(chǎn)制造廠家產(chǎn)生大量廢品。因此，水輪機(jī)調(diào)節(jié)系統(tǒng)必須根據(jù)用戶負(fù)荷的變化不斷調(diào)節(jié)水輪發(fā)電機(jī)的輸出功率，水輪機(jī)調(diào)節(jié)系統(tǒng)的優(yōu)劣性決定電能品質(zhì)。水輪發(fā)電機(jī)PID 調(diào)速器參數(shù)優(yōu)化的方法有：下降梯度法[2]、神經(jīng)網(wǎng)絡(luò)[3]、遺傳算法[4]、模擬退火算法[5]等，這些方法對(duì)提高PID 調(diào)速器的調(diào)節(jié)品質(zhì)均有重要的理論價(jià)值，但是下降梯度法要求提取精確的函數(shù)，處理能力有限；神經(jīng)網(wǎng)絡(luò)易于陷入局部最優(yōu)；GA 算法進(jìn)化速度慢，易于早熟。PSO 算法[6]具有并行處理、尋優(yōu)效率高的特點(diǎn)，與結(jié)構(gòu)簡(jiǎn)單、魯棒性強(qiáng)的PID 調(diào)速器相適應(yīng)。

基本粒子群算法在處理水輪機(jī)PID 調(diào)速器時(shí)會(huì)出現(xiàn)局部最優(yōu)或收斂緩慢等問(wèn)題，因此本文提出了一種改進(jìn)的自適應(yīng)PSO 算法，該算法根據(jù)函數(shù)尋優(yōu)策略改進(jìn)慣性權(quán)重對(duì)搜索能力的影響，提高求解精度并加快收斂速率。在此基礎(chǔ)上引入微分進(jìn)化算子[7]選取部分隨機(jī)粒子進(jìn)行變異處理，從而提高了種群的多樣性，擴(kuò)大搜索解空間的區(qū)域，避免早熟和局部最優(yōu)，并利用matlab 對(duì)該模型進(jìn)行仿真優(yōu)化，仿真結(jié)果驗(yàn)證了該方法的正確性和有效性，為水輪機(jī)調(diào)節(jié)系統(tǒng)提供一種可行的控制方案。

1 水輪機(jī)調(diào)節(jié)控制系統(tǒng)matlab模型

水輪機(jī)控制系統(tǒng)主要由PID 控制模塊和水輪機(jī)被控模塊組成，控制模塊[8]系統(tǒng)結(jié)構(gòu)在matlab 建立如圖1 所示模型。其中，控制模塊中的Kp、Ki、Kd 為simulink 系統(tǒng)中的待優(yōu)化參數(shù)，分別表示比例常數(shù)、積分常數(shù)和微分常數(shù)。被控模塊Subsystem 中封裝了水輪機(jī)調(diào)速系統(tǒng)被控對(duì)象，包括機(jī)械液壓子模塊[9]（MHS）、水輪機(jī)及引水系統(tǒng)子模塊[10]（WDS）和發(fā)電機(jī)及負(fù)荷子模塊[11]（GS）。Subsystem 模塊、MHS 模塊、WDS 模塊、GS 模塊的matlab 建模分別如圖2、圖3、圖4 和圖5 所示。

圖1 水輪機(jī)控制系統(tǒng)數(shù)學(xué)模型

圖2 水輪機(jī)調(diào)節(jié)系統(tǒng)被控模塊

圖3 機(jī)械液壓模塊MHS

圖4 水輪機(jī)及引水系統(tǒng)模塊WDS

圖5 發(fā)電機(jī)及負(fù)荷模塊GS

2 改進(jìn)的PSO優(yōu)化算法

2.1 基本PSO 優(yōu)化算法

與魚群算法、Ant 算法相同，PSO 算法是從生物群體智能獲得啟示，利用并行計(jì)算的智能優(yōu)化方法。Kennedy 和Eberhart 在研究鳥(niǎo)類捕食過(guò)程時(shí)，模擬其搜索解區(qū)域方式，將粒子抽象為潛在解，使每個(gè)粒子綁定一個(gè)適應(yīng)度函數(shù)來(lái)決定其適應(yīng)度值。粒子的速度決定粒子移動(dòng)的方向和距離，速度隨自身及其他粒子的移動(dòng)經(jīng)驗(yàn)進(jìn)行動(dòng)態(tài)調(diào)整，從而實(shí)現(xiàn)個(gè)體在可解空間中的尋優(yōu)。粒子在求解空間中運(yùn)動(dòng)，通過(guò)追蹤個(gè)體極值Pbest 和群體極值Gbest 更新個(gè)體位置。

基本粒子群算法思路： 存在一個(gè)n 維解區(qū)域搜索空間，隨機(jī)分布m 個(gè)粒子(X1，X2,……,Xm),每個(gè)粒子均代表該空間中可能解，數(shù)學(xué)描述為n 維向量Xi=（Xi1，Xi2，……,Xin），根據(jù)計(jì)算每個(gè)Particle 的個(gè)體極值Pid=(Pi1,Pi2,……,Pin)與全體極值Pg=(Pg1,Pg2,……,Pgn), 利用速度調(diào)整公式實(shí)時(shí)跟新自身位置最終找出最優(yōu)解。速度調(diào)整公式為：

式中c1和c2取[0,2]，r1和r2為程序自動(dòng)隨機(jī)數(shù)。

自身位置更新公式為：

基本粒子群算法參數(shù)依賴性小，計(jì)算效率高，但是在求解水輪機(jī)PID 調(diào)速器這類復(fù)雜模型過(guò)程中，易出現(xiàn)早熟和局部收斂，導(dǎo)致所得結(jié)果精度較低或者收斂速度過(guò)慢。

2.2 改進(jìn)的PSO 優(yōu)化策略

（1）自適應(yīng)權(quán)重因子調(diào)整策略。

從基本粒子群算法速度調(diào)整公式可以看出，權(quán)重因子w 直接影響速度梯度，從而控制當(dāng)前粒子在n 維空間中搜索解的能力。若w 數(shù)值偏小，粒子的記憶能力差，粒子移動(dòng)緩慢，降低求解效率，若w 數(shù)值偏大，步長(zhǎng)變大，利于搜索全體空間，但精度降低，不利于求得最優(yōu)解，因此本文將采用一種自適應(yīng)調(diào)整權(quán)重因子的方法，其表達(dá)式為：

式中 k，k+1 表示循環(huán)迭代次數(shù)，大量試驗(yàn)經(jīng)驗(yàn)表明：Wmin取0.85，Wmax取2.2 時(shí)，收斂速度快，所求精度高。

通過(guò)限定權(quán)重因子的最小值并引入迭代次數(shù)增強(qiáng)其自適應(yīng)能力后，在算法尋優(yōu)前期，權(quán)重因子偏大，步長(zhǎng)大，擴(kuò)展搜索全域的能力，加快收斂速度。在尋優(yōu)后期，權(quán)重因子偏小，利于在最優(yōu)解附近微調(diào)整，從而利于調(diào)整精度。

（2）微分進(jìn)化調(diào)整策略。

執(zhí)行新的權(quán)重因子更新公式的同時(shí)，引入微分進(jìn)化操作，在每一代的粒子中隨機(jī)選取部分個(gè)體，進(jìn)行變異進(jìn)化，公式為：

式中 k，k+1——粒子進(jìn)化代數(shù)；

F——向量差分因子，其大小控制粒子的放大系數(shù)，本文取0.9。

2.3 改進(jìn)PSO 算法的適應(yīng)度函數(shù)

傳統(tǒng)的PSO 算法優(yōu)化水輪機(jī)PID 調(diào)速器參數(shù)中，引用時(shí)間絕對(duì)誤差乘積指標(biāo)ITAE 指標(biāo)作為粒子的適應(yīng)度函數(shù)，其表達(dá)式為：

式中 t——時(shí)間；

ts——積分時(shí)間上限；

e（t）——誤差。

為保證水輪機(jī)調(diào)速控制器具有滿意的動(dòng)態(tài)響應(yīng)的特性，本文在ITAE 的基礎(chǔ)上，綜合考慮最大轉(zhuǎn)速偏差相對(duì)值Xmax及控制量平方的積分，使用帶加權(quán)系數(shù)的ITAE 作為粒子的適應(yīng)度函數(shù)，其表達(dá)式為：

3 matlab中算法執(zhí)行過(guò)程

在matlab 的simulink 模塊建立本文圖1～圖5所示數(shù)學(xué)模型，在M 文件中編寫改進(jìn)PSO 算法程序，二者交互過(guò)程如圖6 所示。

改進(jìn)PSO 算法優(yōu)化水輪機(jī)PID 調(diào)速器算法步驟如下：

①以隨機(jī)方式賦值給每個(gè)Particle 原始位置和原始速度。②根據(jù)適應(yīng)度函數(shù)計(jì)算出粒子適應(yīng)度值。③將每個(gè)粒子的適應(yīng)度值與其歷史位置對(duì)比，優(yōu)異者作為Pbest。④將粒子的適應(yīng)度值與整個(gè)粒子群的最優(yōu)歷史位子作對(duì)比，優(yōu)異者作為整個(gè)粒子群的Gbest。⑤根據(jù)式（1）和式（3）更新每個(gè)Particle 的移動(dòng)方向和速度。根據(jù)式（2）更新粒子當(dāng)前位置。⑥根據(jù)式（4）隨機(jī)選取部分粒子進(jìn)行變異操作。⑦根據(jù)式（6）輸出性能指標(biāo)，若滿足條件跳出循環(huán)，否則轉(zhuǎn)第②步。

圖6 算法流程圖

4 工程實(shí)例分析

4.1 計(jì)算數(shù)據(jù)

已知湖南農(nóng)村某小水電站1#混流式水輪機(jī)裝機(jī)容量700 kW，額定水頭97.9 m，最小水頭96 m，最大水頭99 m，額定流量0.71 m3/s，額定轉(zhuǎn)速1 000 r/min，采用單機(jī)孤網(wǎng)帶負(fù)荷模式運(yùn)行。改進(jìn)PSO 算法的試驗(yàn)參數(shù)設(shè)置范圍如本文2.1，2.2，2.3 所述。

4.2 空載工況

將模型置于空載條件下，對(duì)試驗(yàn)仿真機(jī)加5%的頻率擾動(dòng)，研究機(jī)組轉(zhuǎn)速過(guò)渡過(guò)程，其轉(zhuǎn)速偏差結(jié)果如圖7 所示: 藍(lán)紅線分別代表機(jī)組經(jīng)過(guò)改進(jìn)PSO算法優(yōu)化前后的轉(zhuǎn)速響應(yīng)曲線，經(jīng)過(guò)改進(jìn)PSO 算法優(yōu)化過(guò)的水輪機(jī)PID 調(diào)速器響應(yīng)曲線超調(diào)量大幅度減小，調(diào)整時(shí)間也一半，整個(gè)過(guò)程有明顯改進(jìn)。機(jī)組輸出性能指標(biāo)結(jié)果如圖8 所示：ITAE 進(jìn)入穩(wěn)定的時(shí)間加快，且穩(wěn)定值下降幅度較大。

為驗(yàn)證改進(jìn)權(quán)重慣性因子對(duì)PSO 算法的影響，將改進(jìn)PSO 算法與基本PSO 算法的尋優(yōu)收斂曲線進(jìn)行對(duì)比，結(jié)果如圖9 所示：在迭代次數(shù)小于20 時(shí)，紅色收斂曲線的下降幅度大于藍(lán)色線，說(shuō)明在算法早期，改進(jìn)PSO 算法收斂速度較快，當(dāng)?shù)螖?shù)大于20 時(shí)，紅色收斂曲線線不斷作微調(diào)整，說(shuō)明改進(jìn)PSO 算法精度作微調(diào)整。試驗(yàn)結(jié)果與權(quán)重慣性因子特性相吻合。

圖7 轉(zhuǎn)速偏差結(jié)果

圖8 機(jī)組輸出性能指標(biāo)結(jié)果

圖9 尋優(yōu)收斂曲線對(duì)比圖

圖10 10%試驗(yàn)負(fù)載下轉(zhuǎn)速偏差結(jié)果

圖11 輸出特性指標(biāo)

圖12 尋優(yōu)收斂曲線對(duì)比圖

4.3 負(fù)荷工況

為研究水輪機(jī)調(diào)節(jié)的動(dòng)態(tài)特性，將機(jī)組置于10%試驗(yàn)負(fù)荷下，其機(jī)組轉(zhuǎn)速偏差結(jié)果和輸出特性指標(biāo)ITAE 分別如圖10，圖11 所示，結(jié)果表明：經(jīng)過(guò)PSO 算法優(yōu)化過(guò)的PID 水輪機(jī)調(diào)速器在孤網(wǎng)10%負(fù)荷擾動(dòng)下，機(jī)組轉(zhuǎn)速偏差的過(guò)渡過(guò)程比未經(jīng)優(yōu)化的機(jī)組過(guò)渡過(guò)程有明顯改進(jìn)；超調(diào)量減少、調(diào)節(jié)時(shí)間縮短、機(jī)組轉(zhuǎn)速偏差的ITAE 指標(biāo)下降。

同理，將改進(jìn)PSO 算法與基本PSO 算法的尋優(yōu)收斂曲線做比較，結(jié)果如圖12 所示:在整個(gè)過(guò)程中，兩種算法的調(diào)準(zhǔn)精度相差較小，但改進(jìn)PSO 算法的收斂速度較快。優(yōu)化性能較好。

5 結(jié) 論

為了提高傳統(tǒng)PSO 算法對(duì)PID 型水輪機(jī)調(diào)速器Ki,Kp,Kd 三個(gè)參數(shù)的優(yōu)化質(zhì)量，本文選擇自適應(yīng)調(diào)整權(quán)重慣性因子，該方法綜合考慮了實(shí)際粒子群算法迭代過(guò)程中的收斂步長(zhǎng)，平衡收斂速度和精確度的影響，避免陷入局部最優(yōu)的同時(shí)提高了計(jì)算效率，同時(shí)引入微分進(jìn)化因子，增加了全局搜索能力，確保了迭代過(guò)程中出現(xiàn)結(jié)果接近而早熟，最后通過(guò)5%頻率擾動(dòng)機(jī)組仿真試驗(yàn)和10%負(fù)荷擾動(dòng)試驗(yàn)證明了改進(jìn)PSO 算法優(yōu)于基本PSO 算法，經(jīng)過(guò)改進(jìn)PSO 算法優(yōu)化過(guò)的PID 調(diào)速器也具有較好的性能。

[1] 程遠(yuǎn)楚，張江濱.水輪機(jī)自動(dòng)調(diào)節(jié)[M].北京：中國(guó)水利水電出版社，2010.

[2] 湯儀平，金福江. 最速下降法和共軛梯度的混合算法及全局收斂[J]. 華僑大學(xué)學(xué)報(bào)（自然科學(xué)版），2007，（2）:124-126.

[3] 侯樹(shù)文，王琳琳，盧家濤，等. BP 網(wǎng)絡(luò)模糊PID 控制在水輪機(jī)調(diào)速器中的應(yīng)用[J]. 人民長(zhǎng)江，2008，39（2）: 79-80.

[4] 孫美鳳，張永華，張俊紅. 水輪機(jī)調(diào)速器參數(shù)自整定遺傳算法的研究[J]. 水力發(fā)電，2007，（5）.

[5] 王勇，劉厚林，袁壽其，等.離心泵非設(shè)計(jì)工況空化振動(dòng)噪聲的試驗(yàn)測(cè)試[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2012，28（2）：35-38.

[6] 魚瑞文，周佩華，錢苗旺，等. 基于模糊復(fù)合最優(yōu)粒子群算法的礦用電機(jī)控制系統(tǒng)參數(shù)優(yōu)化[J]. 煤礦機(jī)械，2012，33（3）.

[7] 劉自發(fā)，張建華. 一種求解電力經(jīng)濟(jì)負(fù)荷分配問(wèn)題的改進(jìn)微分進(jìn)化算法[J]. 中國(guó)電機(jī)工程學(xué)報(bào)，2008，（10）:100-105.

[8] 魏守平.水輪機(jī)調(diào)節(jié)系統(tǒng)仿真[M].武漢：華中科技大學(xué)出版社，2011.

[9] 魏守平.水輪機(jī)調(diào)節(jié)系統(tǒng)的MATLAB 仿真模型[J]. 水電自動(dòng)化與大壩監(jiān)測(cè)，2009，（4）.

[10] 曹程杰，莫岳平，劉丹丹. 基于改進(jìn)遺傳算法的水輪機(jī)PID 調(diào)速器最優(yōu)參數(shù)整定[J].微計(jì)算機(jī)信息，2009，（34）.

[11] 羅旋，魏守平.基于MATLAB 的水輪機(jī)調(diào)節(jié)系統(tǒng)仿真[J].西北水電，2005，（3）.