（國網(wǎng)陜西省電力公司電力科學(xué)研究院，陜西 西安 710100）

統(tǒng)一潮流控制器（UPFC）是目前柔性交流輸電系統(tǒng)（flexible AC transmission system，F(xiàn)ACTS）中功能最強大、綜合性能最好的裝置[1-2]。由于UPFC集合了靜止同步補償器（static synchronous compensator，STATCOM）和靜止同步串聯(lián)補償器（static synchronous series compensator，SSSC）的優(yōu)點，在互聯(lián)電網(wǎng)潮流控制、穩(wěn)定電壓、提高電網(wǎng)穩(wěn)定性方面均表現(xiàn)出良好的應(yīng)用效果，因此其控制策略的研究是當(dāng)前電力系統(tǒng)研究的熱點[3]。文獻[4]通過將UPFC的有功和無功潮流解耦設(shè)計了UPFC并聯(lián)側(cè)的控制器，并通過三環(huán)解耦設(shè)計了其串聯(lián)側(cè)的控制器，但是這種控制方式破壞了UPFC的非線性特性。文獻[5]利用暫態(tài)能量函數(shù)法和模糊滑?？刂品ㄔO(shè)計了UPFC的控制器來抑制系統(tǒng)的低頻振蕩，但控制器設(shè)計基于的模型是系統(tǒng)等效模型。文獻[6]通過構(gòu)造UPFC的Lyapunov函數(shù)和坐標變換完成了UPFC控制器設(shè)計，但是在控制器設(shè)計過程中仍進行了部分線性化處理。由于UPFC裝置及電力系統(tǒng)都是非線性系統(tǒng)，因此需要利用非線性方法設(shè)計系統(tǒng)的穩(wěn)定控制器。

Hamilton系統(tǒng)理論上可以根據(jù)系統(tǒng)的內(nèi)部結(jié)構(gòu)特征，在不經(jīng)過任何線性化處理的基礎(chǔ)上完成控制器的設(shè)計，因此在非線性系統(tǒng)控制方面具有廣闊的應(yīng)用前景[7]。文獻[8]利用廣義Hamilton理論設(shè)計了TCSC（thyristor controlled series compensator）的控制器并在仿真平臺上搭建實際電網(wǎng)模型驗證了其控制效果。文獻[9]通過無源控制方法設(shè)計了電網(wǎng)穩(wěn)定控制器，并在設(shè)計中考慮了相關(guān)元件的轉(zhuǎn)移電導(dǎo)。文獻[10]在考慮發(fā)電機廣域勵磁信號存在時滯的前提下，通過偽廣義Hamilton理論設(shè)計了STATCOM與時滯勵磁的協(xié)調(diào)控制器。但是，利用Hamilton系統(tǒng)理論設(shè)計UPFC控制器的研究還很欠缺。

本文在考慮系統(tǒng)相關(guān)元件轉(zhuǎn)移電導(dǎo)的情況下根據(jù)Hamilton系統(tǒng)理論設(shè)計了多機勵磁和UPFC的協(xié)調(diào)控制策略。在定義出UPFC對電網(wǎng)影響的前提下，建立了包含UPFC裝置的多機電力系統(tǒng)數(shù)學(xué)模型。進而將此系統(tǒng)表示成含參數(shù)擾動的廣義Hamilton系統(tǒng)形式。根據(jù)邊界函數(shù)法和L2干擾抑制控制法設(shè)計了相應(yīng)的協(xié)調(diào)控制器。系統(tǒng)動態(tài)響應(yīng)結(jié)果表明了本文所提非線性協(xié)調(diào)控制策略的正確性和良好應(yīng)用效果。

1 廣義Hamilton理論L2干擾抑制控制法

考慮如下的含有擾動項的仿射非線性系統(tǒng)：

式中：x為系統(tǒng)的n維狀態(tài)變量，x∈ Rn；f(x)為n維微分流形；u為m維控制輸入向量，u∈Rm；w為系統(tǒng)的s維擾動項，w∈ Rs；g(x)，G(x)為適當(dāng)階數(shù)的系數(shù)矩陣。

若存在一個標量函數(shù)H(x)，則能將式（1）所示的系統(tǒng)表示為

式（2）中的∑(x)在x的定義域上是一致有界的，并且在系統(tǒng)的平衡點x0處的取值為零，此式中其余變量的含義可以參考文獻[11]。根據(jù)邊界函數(shù)理論[12]，∑(x)的表達式需要滿足下式：

式中：ρ(x)是一個正定函數(shù)。

由于系統(tǒng)中含有參數(shù)擾動，則考慮如下的懲罰函數(shù)：

式中：h為權(quán)重矩陣。

當(dāng)式（2）所示系統(tǒng)滿足以下三點要求時：

1）式（2）所示系統(tǒng)是零狀態(tài)可檢的；

2）系統(tǒng)的Hamilton函數(shù)H(x)為非負定函數(shù)；

3）給定的干擾抑制水平γ（大于零）滿足：

則可令控制變量u滿足下式，從而使系統(tǒng)在平衡點x0處漸近穩(wěn)定：

式中：λ為滿足R(x)-λ/2gˉgˉT≥ 0的非負數(shù)；I為單位對角矩陣。

上述控制策略的推導(dǎo)過程詳見文獻[13]。

2 含UPFC裝置的多機電力系統(tǒng)廣義Hamilton系統(tǒng)實現(xiàn)

本文所研究的n機系統(tǒng)包含m個負荷節(jié)點和一臺UPFC裝置。其中發(fā)電機采用三階模型，1,2,…,n為其節(jié)點；n+1和n+2為安裝有UPFC的母線節(jié)點；n+3,…,n+m+2為負荷節(jié)點。

2.1 UPFC的動態(tài)模型

UPFC裝置可以看成是由一臺STATCOM和一臺SSSC合并而成的，其模型結(jié)構(gòu)圖如圖1所示。

UPFC的非線性動態(tài)方程[14]如下：

式中，id1，iq1分別為UPFC并聯(lián)部分注入的d軸和q軸電流；id2，iq2分別為UPFC串聯(lián)部分注入的d軸和q軸電流；Un+1∠θn+1，Un+2∠θn+2分別為UPFC并聯(lián)部分和串聯(lián)部分接入母線的電壓；Rs1，Ls1分別為并聯(lián)耦合變壓器的等效電阻和漏電感；Rs2，Ls2分別為串聯(lián)耦合變壓器的等效電阻和漏電感；Cdc為直流側(cè)電容器的電容；Rdc為兩個換流器的等效切換損失；Udc為直流側(cè)電容電壓；ωs為系統(tǒng)角頻率；ζ1,ζ2,ζ3,ζ4為UPFC的等效有界輸入擾動；k1，k2分別為并聯(lián)側(cè)換流器調(diào)制比和串聯(lián)側(cè)換流器調(diào)制比；α1，α2分別為并聯(lián)側(cè)換流器相位移和串聯(lián)側(cè)換流器相位移；其中的ωs為有名值，其余變量均為標幺值。

定義UPFC的控制輸入變量為u1，u2，u3和u4，其表達式為

圖1 UPFC的電路結(jié)構(gòu)圖Fig.1 Circuit structure diagram of UPFC

2.2 多機電力系統(tǒng)模型

將由于UPFC裝置的存在，電力系統(tǒng)中發(fā)電機輸出的附加電磁功率和d軸電流定義為ΔPei和ΔIdi。則在n機電力系統(tǒng)中第i臺發(fā)電機的三階動態(tài)方程為

式中：下標i為第i臺發(fā)電機的相應(yīng)變量，i=1,2,…,n；ωi為發(fā)電機角速度；ω0為是發(fā)電機同步角速度；xdi為發(fā)電機的d軸同步電抗；x′di為發(fā)電機的d軸暫態(tài)電抗；T′d0i為發(fā)電機的d軸暫態(tài)開路時間常數(shù)；E′qi為發(fā)電機的q軸暫態(tài)電勢；Gij，Bij分別為發(fā)電機i與發(fā)電機j間的互電導(dǎo)及互電納；δij為E′qi的復(fù)矢量與E′qj的復(fù)矢量之間的相角差；Di，Mi分別發(fā)電機的阻尼系數(shù)和慣性時間常數(shù)；Pmi為固定的機械輸入功率；ufi為發(fā)電機勵磁輸入控制量；Efis為平衡狀態(tài)下勵磁電壓值；wi1，wi2為發(fā)電機的有界輸入擾動。

2.3 系統(tǒng)的廣義Hamilton形式實現(xiàn)

在得到UPFC裝置和第i臺發(fā)電機的數(shù)學(xué)模型后，聯(lián)立式（5）和式（7）就可以得到含有UPFC裝置的多機電力系統(tǒng)的數(shù)學(xué)模型表達式。根據(jù)實際電力系統(tǒng)將x=[δ1,ω1,E′q1,…,δn,ωn,E′qn,id1,iq1,id2,iq2,Udc]T作為其狀態(tài)變量，系統(tǒng)的初始平衡點為x0。為了簡便，后續(xù)表述定義：

構(gòu)造上述系統(tǒng)的Hamilton能量函數(shù)為

此能量函數(shù)H(x)具有明確的物理意義，代表了包含UPFC裝置附加能量在內(nèi)的系統(tǒng)總能量。

計算能量函數(shù)H(x)對各個狀態(tài)變量的偏導(dǎo)數(shù)，則多機系統(tǒng)的動態(tài)方程可表示為

式中：u為多機電力系統(tǒng)的控制輸入量；u10，u20，u30，u40分別為UPFC在系統(tǒng)平衡狀態(tài)下的控制變量初始值；id10，iq10，id20，iq20為UPFC在系統(tǒng)平衡狀態(tài)下的初始值；w為多維擾動矩陣。

3 發(fā)電機勵磁與UPFC裝置的協(xié)調(diào)控制器設(shè)計

將裝設(shè)有UPFC的多機電力系統(tǒng)表示成廣義Hamilton形式后，要想利用式（4）所示的擴展L2干擾抑制控制法設(shè)計系統(tǒng)的協(xié)調(diào)控制器，需要首先驗證一下系統(tǒng)的幾個性質(zhì)：

1）當(dāng)控制變量u=0時，很明顯多機系統(tǒng)是零狀態(tài)可檢的。

2）從H(x)的表達式可以看出，在x0處，H(x)的值為零，并且?xH(x)|x=x0=0。要想保證H(x)為非負定函數(shù)，就要求其Hessian矩陣正定。但是像H(x)這樣一個復(fù)雜函數(shù)的Hessian矩陣的正定條件很難用具體表達式給出，故通常采用計算機數(shù)值計算法進行判定[15]。

3）通過對下式進行計算來確定系統(tǒng)的干擾抑制水平γ：

并且從∑(x)的表達式可以驗證∑(x0)=0且在x的定義域空間中一致有界。根據(jù)式（3）可定義出ρ(x)=[ρ1i(x),ρ2(x)]T，只要找到合適的ρ(x)就可以根據(jù)式（4）來設(shè)計多機勵磁與UPFC裝置的非線性協(xié)調(diào)控制策略。

多機系統(tǒng)中第i臺發(fā)電機勵磁的控制策略為

式中：hi為特定系數(shù)。

UPFC裝置的控制策略為

式中：h2為待定系數(shù)。

由式（6）可將u10，u20，u30，u40轉(zhuǎn)化成UPFC裝置的觸發(fā)控制變量：

4 仿真分析與驗證

為驗證本文所提的多機勵磁與UPFC的協(xié)調(diào)控制策略的有效性，在PSCAD仿真環(huán)境中搭建如圖2所示的包含UPFC裝置的四機兩區(qū)域仿真模型。此仿真模型中所有發(fā)電機的額定容量為900 MV·A，額定電壓是20 kV。系統(tǒng)輸出電壓的額定值為230 kV。以額定值為基準值的發(fā)電機和UPFC裝置的具體仿真參數(shù)分別如表1和表2所示。根據(jù)系統(tǒng)在平衡點處的初始值和第1節(jié)所述控制理論的約束條件通過多次仿真優(yōu)化法確定控制器參數(shù)，取λ1i=λ2=0.1，hi=h2=2，γ=2，邊界函數(shù)取通過PSCAD/EMTDC仿真平臺中測量單元的測量值來模擬實際系統(tǒng)中PMU裝置的測量值，為本仿真提供所需的發(fā)電機和UPFC裝置的實時運行參數(shù)。

圖2 四機兩區(qū)域系統(tǒng)結(jié)構(gòu)圖Fig.2 The system structure of two-area including four generators

表1 發(fā)電機仿真參數(shù)Tab.1 Generator simulation parameters

表2 UPFC裝置仿真參數(shù)Tab.2 UPFC device simulation parameters

本仿真中系統(tǒng)初始運行于穩(wěn)定平衡點，設(shè)定的系統(tǒng)故障情況如下：1）UPFC裝置串聯(lián)母線右側(cè)在2 s時發(fā)生三相瞬時性接地短路故障，0.1 s后故障切除，系統(tǒng)恢復(fù)正常；2）母線7處的負荷在2 s時忽然減少為原來的50%，并在2.5 s時恢復(fù)正常負荷。為了充分驗證本文所提非線性控制策略的有效性，將其與由發(fā)電機的AVR/PSS和UPFC的PI控制組成的分散控制進行對比仿真。

故障1的設(shè)定用來驗證非線性協(xié)調(diào)控制器的有效性；故障2的設(shè)定用來驗證非線性協(xié)調(diào)控制器的魯棒性。圖3～圖8分別展示了系統(tǒng)聯(lián)絡(luò)線傳輸功率、發(fā)電機功角差及母線9電壓在兩種故障情況下的動態(tài)響應(yīng)曲線。

圖3 G1和G3的功角差Δδ13的響應(yīng)曲線（故障1）Fig.3 Response curves of the power angle difference Δδ13betweenG1andG3（fault 1）

圖4 聯(lián)絡(luò)線傳輸功率P7-8的響應(yīng)曲線（故障1）Fig.4 Response curves of tie line transmission powerP7-8（fault 1）

圖5 母線9電壓的響應(yīng)曲線（故障1）Fig.5 Response curves of bus 9 voltage（fault 1）

圖6 G1和G3的功角差Δδ13的響應(yīng)曲線（故障2）Fig.6 Response curves of the power angle difference Δδ13betweenG1andG3（fault 2）

圖7 聯(lián)絡(luò)線傳輸功率P7-8的響應(yīng)曲線（故障2）Fig.7 Response curves of tie line transmission powerP7-8（fault 2）

圖8 母線9電壓的響應(yīng)曲線（故障2）Fig.8 Response curves of bus 9 voltage（fault 2）

從圖3～圖8中可以看出，與傳統(tǒng)分散控制相比，本文設(shè)計的非線性協(xié)調(diào)控制器具有較快的響應(yīng)速度，并且能夠以較小的超調(diào)量在較短的時間內(nèi)使系統(tǒng)迅速恢復(fù)穩(wěn)定，從而驗證了本控制器的有效性和魯棒性。

5 結(jié)論

本文利用廣義Hamilton理論的擴展L2干擾抑制控制法設(shè)計了多臺發(fā)電機勵磁與UPFC裝置的協(xié)調(diào)控制器，解決了具有參數(shù)擾動的含UPFC裝置的多機電力系統(tǒng)的穩(wěn)定性問題，同時系統(tǒng)動態(tài)方程的建立和Hamilton函數(shù)的構(gòu)造為研究多臺UPFC裝置和多臺發(fā)電機的協(xié)調(diào)控制問題提供了基礎(chǔ)。仿真結(jié)果驗證了所提非線性控制器的有效性和魯棒性。本文中的一些控制器參數(shù)是通過多次仿真優(yōu)化得到的。下一步工作的重點是，研究一種適用于確定非線性控制器參數(shù)的簡便優(yōu)化方法。