電力系統(tǒng)穩(wěn)定器的混合差分進(jìn)化算法設(shè)計(jì)研究

摘要：本文采用混合差分進(jìn)化算法設(shè)計(jì)微電網(wǎng)中的穩(wěn)定器。首先，對(duì)單發(fā)電機(jī)對(duì)無線匯流排系統(tǒng)的穩(wěn)定器進(jìn)行研究，變化發(fā)電機(jī)有功、無功功率及輸電線阻抗，采用混合差分進(jìn)化算法，指定不同目標(biāo)函數(shù)極點(diǎn)使穩(wěn)定器工作于復(fù)平面左半部分，以求優(yōu)良好的動(dòng)態(tài)穩(wěn)定性能。然后，再延伸到多機(jī)和復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)的電網(wǎng)中。

　　隨著用電量不斷快速增長(zhǎng)，電力系統(tǒng)穩(wěn)定度下降，為改善穩(wěn)定性加裝穩(wěn)定器^[1]。設(shè)計(jì)穩(wěn)定器的關(guān)鍵是如何調(diào)整參數(shù)以獲得最佳動(dòng)態(tài)特性。過去用過的方法有模糊控制法、分散模態(tài)控制法、多目標(biāo)基因法等，而混合差分進(jìn)化法才是解決線性與非線性最小化問題的最好方法^[2-4]。

　　本文使用混合差分進(jìn)化算法，設(shè)計(jì)單機(jī)與無限匯流排系統(tǒng)中的穩(wěn)定器^[5-7]，然后再將穩(wěn)定器擴(kuò)展使用于多機(jī)和復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)的系統(tǒng)中去。

1 系統(tǒng)模型

　　圖1為單發(fā)電機(jī)與無限匯流排系統(tǒng)，包括發(fā)電機(jī)及動(dòng)態(tài)勵(lì)磁系統(tǒng)，發(fā)電機(jī)使用交、直雙軸模型，其參數(shù)如表1所示^[8]，為系統(tǒng)狀態(tài)參數(shù)，ΔE'd、ΔE'q為d軸q軸的暫態(tài)電壓，Δω與Δδ為角速度和角度，ΔE_FD為勵(lì)磁電壓，ΔV_s為系統(tǒng)內(nèi)部穩(wěn)定器輸出。利用混合為變數(shù)向量的目標(biāo)函數(shù)。此處變數(shù)向量為發(fā)電機(jī)電力系統(tǒng)參數(shù)。

　　圖2為靜態(tài)勵(lì)磁系統(tǒng)的方框圖，參數(shù)如表2所示^[8]。勵(lì)磁系統(tǒng)及發(fā)電機(jī)方程式在某工作點(diǎn)線性化，表示如下：

(1)

　　其中，為系統(tǒng)狀態(tài)參數(shù)，ΔE'd，ΔE'q為d軸q軸的暫態(tài)電壓，Δω與Δδ為角速度和角度，ΔE_FD為勵(lì)磁電壓，ΔV_s為系統(tǒng)內(nèi)部穩(wěn)定器輸出。

　　用混合差分進(jìn)化算法設(shè)計(jì)穩(wěn)定器時(shí)，應(yīng)在一些限制條件下求目標(biāo)函數(shù)最優(yōu)化，即：

　　目標(biāo)函數(shù)： (2)

　　限制條件： (3)

(4)

　　式(2)中，為變數(shù)向量的目標(biāo)函數(shù)。此處變數(shù)向量為發(fā)電機(jī)電力系統(tǒng)參數(shù)。

　　穩(wěn)定器以發(fā)電機(jī)速度偏差為輸入信號(hào)見圖3，其傳遞函數(shù)如下：

(5)

　　若時(shí)間常數(shù)T₅為固定值5 sec，其余參數(shù)K_S、T₁、T₂、T₃和T₄可調(diào)，結(jié)合發(fā)電機(jī)與勵(lì)磁系統(tǒng)可得線性方程為：

(6)

　　其中為系統(tǒng)狀態(tài)變數(shù)，ΔV₁、ΔV₂為電力系統(tǒng)穩(wěn)定器內(nèi)部電壓，ΔV_pss為輸出信號(hào)。

2 穩(wěn)定器設(shè)計(jì)

2.1 穩(wěn)定器設(shè)計(jì)原理

　　為增加系統(tǒng)阻尼，在發(fā)電機(jī)輸出有效功率為P，無效功率為Q，輸電線路電抗X_e變化的條件下，將系統(tǒng)機(jī)電模式于指定的復(fù)平面內(nèi)，定義三種目標(biāo)函數(shù)如下：

(7)

(8)

(9)

　　式中N_i為P范圍內(nèi)總點(diǎn)數(shù)，N_j為Q范圍內(nèi)總點(diǎn)數(shù)，N_k為X_e范圍內(nèi)總點(diǎn)數(shù)，相位超前或落后補(bǔ)償型電力系統(tǒng)穩(wěn)定器的參數(shù)限制為^[9]：

(10)

　　若(7)~(9)中目標(biāo)函數(shù)為零，結(jié)果將使機(jī)電模式實(shí)部皆位于垂線左邊，阻尼比區(qū)域內(nèi)。設(shè)定發(fā)電機(jī)輸出有效功率為P，無效功率Q及輸電線路阻抗X_e取值變化范圍如下^[8]：

2.2 混合差分進(jìn)化法穩(wěn)定器設(shè)計(jì)流程

　　混合差分進(jìn)化法穩(wěn)定器設(shè)計(jì)流程見圖4。用隨機(jī)方式取得差異向量產(chǎn)生一次進(jìn)化中的突變，擴(kuò)大搜尋范圍，通過交配產(chǎn)生進(jìn)化中的下一代，父子兩代經(jīng)選擇環(huán)節(jié)擇其性能優(yōu)者繼續(xù)向前進(jìn)化，直到目標(biāo)函數(shù)獲得最優(yōu)。這就是差分進(jìn)化法^[10-12] 。

　　單純差分進(jìn)化法有缺點(diǎn)。遇到多代子不如其父，目標(biāo)函數(shù)收斂將很慢;又有時(shí)會(huì)錯(cuò)誤地收斂到局部最小解;所以，圖4中增加了遷移及加速環(huán)節(jié)。為此更名為混合差分進(jìn)化法，其詳見后 ^[13-15] 。

　　(1)初始化：以隨機(jī)方式產(chǎn)生N_p個(gè)族群，在每一個(gè)解中有D個(gè)變數(shù)，此變數(shù)均勻分布在整個(gè)求解空間內(nèi)，表示如下：

(11)

　　式(11)中為均勻分布隨機(jī)數(shù)，、為變數(shù)X_j 的最大最小值。

　　(2)突變：利用隨機(jī)方式得到差異向量以產(chǎn)生一個(gè)擾動(dòng)向量，擴(kuò)大搜索范圍，在突變過程中，第G+1代突變向量為：

(12)

　　其中，F為常數(shù)，、為隨機(jī)產(chǎn)生的兩個(gè)向量。

　　(3)交配：交配過程是第G代向量與第G+1代突變向量交換、混合成為G+1實(shí)驗(yàn)向量，表達(dá)式如下:

(13)

　　(4)選擇：如果下一代更好，則替換低G代，相反則保留第G代。首先做競(jìng)爭(zhēng)，然后選擇最佳個(gè)體：

(14)

(15)

　　將以上步驟反復(fù)進(jìn)行，直到有合適結(jié)果為止。針對(duì)差分進(jìn)化法缺點(diǎn)，采用混合差分進(jìn)化法進(jìn)行改善，并將演算步驟代入遷移及加速過程，提高收斂速度及克服局部最小解^[13-15]。