基于Wirtinger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法

【專利摘要】本發(fā)明公開了一種基于Wirtinger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法，用于分析電力系統(tǒng)所能承受的最大時(shí)滯穩(wěn)定裕度。該方法的具體步驟如下：首先，建立考慮時(shí)滯影響的電力系統(tǒng)模型。然后，針對(duì)所建模型構(gòu)建Lyapunov泛函，在泛函的求導(dǎo)過程中通過采用Wirtinger不等式進(jìn)行放縮，以減少判據(jù)的保守性。最后將所得判據(jù)用一組線性矩陣不等式(LMI)表示。

【專利說明】

基于W i rt i nger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法

技術(shù)領(lǐng)域

[0001 ]本發(fā)明涉及一種基于Wirtinger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法，適用于 解決互聯(lián)電力系統(tǒng)廣域控制策略中的延時(shí)問題，屬于電力系統(tǒng)技術(shù)領(lǐng)域。

【背景技術(shù)】

[0002] 隨著電網(wǎng)規(guī)模的不斷擴(kuò)大，單純的依靠本地控制信息已經(jīng)無法滿足電網(wǎng)性能的要 求。因此電網(wǎng)穩(wěn)定控制策略逐步從本地控制趨向于全局控制，特別是近年來廣域測(cè)量系統(tǒng) (WAMS)的快速發(fā)展，使電力系統(tǒng)全局控制策略的實(shí)現(xiàn)成為可能。

[0003] 在此基礎(chǔ)上，現(xiàn)代電力系統(tǒng)控制器將從全局角度出發(fā)，利用本地信號(hào)和遠(yuǎn)端信號(hào) 作為反饋信號(hào)設(shè)計(jì)廣域控制器。由于遠(yuǎn)端信號(hào)的引入，信號(hào)的延時(shí)將變得不可避免，已有研 究表明，即使很小的時(shí)滯都可能對(duì)電力系統(tǒng)穩(wěn)定性產(chǎn)生影響。因此考慮電力系統(tǒng)所能承受 的最大時(shí)滯，對(duì)于電網(wǎng)的安全穩(wěn)定運(yùn)行具有十分重要的意義。

[0004] 目前關(guān)于時(shí)滯系統(tǒng)的研究，具有多種求解方法，常采用的求解方法是構(gòu)造 Lyapunov泛函，基于Lyapunov穩(wěn)定性理論，得到系統(tǒng)穩(wěn)定判據(jù)，最后借助線性矩陣不等式 (LMI)來求解時(shí)滯穩(wěn)定裕度。但是由上述方法得到的判據(jù)，實(shí)際是對(duì)求導(dǎo)后的Lyapunov泛函 放縮所得到的充分條件，因此造成所得判據(jù)具有一定的保守性。所以降低保守性問題成為 該方法研究的重點(diǎn)，也是難點(diǎn)之一。

【發(fā)明內(nèi)容】

[0005] 發(fā)明目的：針對(duì)現(xiàn)有技術(shù)中存在的問題，本發(fā)明提出一種基于Wirtinger不等式的 時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法，首先構(gòu)造全新的Lyapunov泛函，將時(shí)滯下限不為零考慮進(jìn) 判據(jù)中，然后利用Wirtinger不等式放縮技巧，降低結(jié)果的保守性。

[0006] 技術(shù)方案:一種基于Wirtinger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法，包括如下 步驟： j x(t) = Ax(t) + A,)X{t - ci(t))

[0007] (1)建立包含廣域控制回路的時(shí)滯電力系統(tǒng)模型],、，、' \式中：

[x(j) = (p(j) /eH/{0,0] hi彡d(tXh2,辦)，y彡1。其中0彡hi

[0008] (2)給定穩(wěn)定判定條件：

[0009] 若存在正定矩陣P G R4nX4n;正定矩陣Qi G RnXn，i = 1，2，3;正定矩陣Zj G RnXn，j = 1， 2;矩陣XkGR2nX2n，k=l，2使下列矩陣不等式成立，則時(shí)滯電力系統(tǒng)是漸近穩(wěn)定的：

[0012] 其中:(1)。= //e((7/ /)(72) + Q + 4 + 4 + 勻允2

[0013] r：!-^ Gf Gl G[]r

[0014] G3 =ei_e2，G4=ei+e2_2es，Gs= e2_e4，G6= e2+e4_2e7

[0015] r2=[Gf Ggr 0i

[0016] Gy = e3_e2，G8 = e3+e2_2e6，G9 = e2_e4，Gio = e2+e4_2e7

[0017] Qj -diag^O.-Q^J> Q： = ^^g(〇2^°2,；'Q： °.i/,);

[0018] Q,,=diag{0,,-{\-d)0,,Q5tt) _ Z^h2 xdiagiZ^Q.J + Uh-l\)^diag{Z^0,:)

[0019] 之=c//tzg(Z丨，3Z|)，之2 =r"ag(Z2,3Z2)

[0020] G^\e[ d(!)el {d{l)-h{)el {h.-dU^J 一 T

[0021 ] (?2 = A^C-, ~ (1 - cl)~ ~ cl)c!j -+ (1- d)C^

[0022] ei=[0i-i I O7-1], 1 = 1, , 7

[0023] 對(duì)于矩陣A，He(A)=A+AT。其中I代表單位矩陣。

[0024] (3)利用Matlab中的線性矩陣(LMI)工具箱判斷給定時(shí)滯d(t)是否滿足步驟(2)給 出的判定條件，若滿足，則可判定在延時(shí)d(t)條件下的時(shí)滯電力系統(tǒng)是漸近穩(wěn)定的。

[0025] 時(shí)滯電力系統(tǒng)模型乂(〇 =也的+ .^^_

【附圖說明】

[0026] 圖1為本發(fā)明所提的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法流程圖；

[0027]圖2為包含WADC的時(shí)滯電力系統(tǒng)；

[0028]圖3為四機(jī)兩區(qū)域電力系統(tǒng)；

[0029]圖4為d(t)=0ms情況下四機(jī)兩區(qū)域電力系統(tǒng)響應(yīng)圖；

[0030]圖5為d(t) = 110ms情況下四機(jī)兩區(qū)域電力系統(tǒng)響應(yīng)圖。

【具體實(shí)施方式】

[0031]下面結(jié)合具體實(shí)施例，進(jìn)一步闡明本發(fā)明，應(yīng)理解這些實(shí)施例僅用于說明本發(fā)明 而不用于限制本發(fā)明的范圍，在閱讀了本發(fā)明之后，本領(lǐng)域技術(shù)人員對(duì)本發(fā)明的各種等價(jià) 形式的修改均落于本申請(qǐng)所附權(quán)利要求所限定的范圍。

[0032] (1)電力系統(tǒng)時(shí)滯模型的建立

[0033]通常情況下，電力系統(tǒng)可由一組微分代數(shù)方程描述，在系統(tǒng)運(yùn)行點(diǎn)附近對(duì)其線性 化，最終系統(tǒng)可表示為：

[i, (/) = A,x, (/) + _ KU ⑴ | j(，) = ChV|(0

[0035] 式中：xi為系統(tǒng)狀態(tài)變量;Ai為系統(tǒng)狀態(tài)矩陣;Bi為系統(tǒng)輸入矩陣;Ci為系統(tǒng)輸出矩 陣;u為系統(tǒng)控制輸入，y為系統(tǒng)控制輸出。

[0036] 系統(tǒng)的控制策略采用基于WAMS的廣域阻尼控制，其中廣域阻尼控制器(WADC)的控 制輸入信號(hào)含有系統(tǒng)遠(yuǎn)端信號(hào)，其狀態(tài)方程可表示為： f.v {t) = Ax U) + Bu (t)

[0037] ^ f ( (2) UV⑴=(、弋⑴+ Z). W?

[0038] 式中：x。為控制器的狀態(tài)變量;A。為控制系統(tǒng)的狀態(tài)矩陣;B。為控制系統(tǒng)的輸入矩 陣;C。為控制系統(tǒng)的輸出矩陣;u。為控制系統(tǒng)輸入，y。為系控制系統(tǒng)輸出。D。為標(biāo)量，反映了 輸出y。與輸入u。之間的直接關(guān)聯(lián)。

[0039] 圖2給出了包含廣域控制回路的時(shí)滯電力系統(tǒng)結(jié)構(gòu)關(guān)系圖，其中的變延時(shí)d(t)是 由于遠(yuǎn)端信號(hào)傳輸所引起的傳輸延時(shí)。根據(jù)圖2所示的關(guān)系圖，可以得到：

[0040] (3) |w(/) = V(.(〇

[0041 ]其中：hi彡d(tXh2,

[0042]進(jìn)一步，時(shí)滯電力系統(tǒng)模型可表示為：

[0043] i(/) = Ax{() + Adx{t-d{t)) (4). T 「4 B.C 1 「5,1) C, 0]

[0044] 式中，4= 1 e 1 L 丨 c」 L〇 A J L sc\

[0045] (2)電力系統(tǒng)延時(shí)依賴穩(wěn)定性判據(jù)方法

[0046]求解時(shí)滯穩(wěn)定判據(jù)過程中，常用的放縮技巧是利用Jensen不等式。這種方法雖然 可行但增加了結(jié)果的保守性。本發(fā)明所提方法采用一種全新的不等式-Wirtinger不等式進(jìn) 行放縮，可大大降低所得結(jié)果的保守性。首先，給出本發(fā)明所提方法用到的兩個(gè)重要的引 理。

[0047]引理1:對(duì)于給定的正定矩陣M>0,以下不等式對(duì)于在區(qū)間[a，b]上連續(xù)可微函數(shù)x 都成立：

[0051 ]引理2:對(duì)于給定的正定矩陣R>0，矩陣I，W2和標(biāo)量a G (〇，1)，定義對(duì)于所有的I， 函數(shù) 0 (a，R):

[0053] 如果存在矩陣乂，使A X >0,那么以下不等式成立： -* 'R XTW￡~

[0054] min 0(a., R)> " (仙 Lmi」L*及」L:i^i

[0055] 構(gòu)造如下Lyapunov泛函： V ~ if Pq i- [ x1 (s)Q,x(s)cls + I x! (s)()2x(s)cis + [ x1 (s)Q,x(s)ds r n h-h h--h2

[0056] 〇 , _h t (5) + J [ i：1 (a0Z + | ^ j" x1 (s)Z lx(.s)dsci0

[0057] 式中：PGR4nX4n;QiGRnXn，i = l，2,3;ZjGRnXn，j = l，2，

[0058] = x1 (t) | ^xr(s)d$ I . xJ (syds

[0059] SeiGRM'eftOi-丨I 〇7-」，1 = 1，…，7為分塊坐標(biāo)矩陣，可得如下穩(wěn)定性判據(jù)：

[0060] 判據(jù):若存在正定矩陣P G R4nX4n;正定矩陣Qi G RnXn，i = 1，2，3 ;正定矩陣Zj G RnXn， j = l，2;矩陣XkGR2nX2n，k=l，2使下列矩陣不等式成立，則時(shí)滯電力系統(tǒng)是漸近穩(wěn)定的：

(6)

[0063] 其中：c丨)(1 = /"，((，丫/%) + Q + 4 + 4 + (7:丨f

[0064] ri: = [Gf GTA Gl Gl^

[0065] 〇3 =ei_e2，G4=ei+e2_2es，Gs= e2_e4，G6= e2+e4_2e7

[0066] r2: = [(3f Gj G^J

[0067] 〇7 = e3_e2，G8 = e3+e2_2e6，G9 = e2_e4，Gio = e2+e4_2e7

[0068] Q = , 4 =t/Vag(a,02;,,

[0069] (), -^K2',^°5") 7 Z ^h? y.diag(Zi,0,ii) + (h2-h^xdiagiZ,,%^

[0070] 氧-―g(Z丨，耳)，之二-g(Z2,3Z2)

[0071] 0；=[< (/(0< (d ⑴-(/z?-

[0072] Gj - - (1.-一(1 一 -4+:(1 -

[0073] ei=[0i-i I 07-i] ,i = l ,7 [0074]那么系統(tǒng)(4)漸近穩(wěn)定。

[0075] 對(duì)于矩陣A，He(A)=A+AT。其中I代表單位矩陣。

[0076]證明：

[0077] 對(duì)于判據(jù)中的Lyapunov泛函進(jìn)行求導(dǎo)可得：

[0081 ] 然后對(duì)求導(dǎo)后的Lyapunov泛函(7)的最后兩項(xiàng)利用引用1,2處理，以

[0083]針對(duì)式(9)分別使用引理1可得： 處理方法相同，具體操作如下： (9)

[0085] 其中：

[0086] 111 = 63^1 ^12=64^1 ^21 = G5Cl^22 = G6ll [0087] 針對(duì)上式，使用引理2可得：

[0090]另一項(xiàng)-「了 辦采用同樣的處理方法：

[0092] 將處理后的兩項(xiàng)代入到求導(dǎo)后的Lyapunov泛函中去可得：

[0096]系統(tǒng)漸近穩(wěn)定，需要使所求Lyapunov泛函導(dǎo)數(shù)小于0,判據(jù)得證。

[0097]注：判據(jù)中的不等式(6)依賴于d (t)和,無法直接使用LMI工具箱求解。但 是不等式(6)是關(guān)于d(t)和的凸函數(shù)，所以只要使上述不等式在cKUzln， d(t) = -jii ： d(t) - hx . d(t) = j.i \ d{l) -h2 , d{t) = -pi ； d{l)二 h2、?'(〇 = "上都成立即 可。

[0098]下面介紹本發(fā)明的一個(gè)實(shí)施例：

[0099]四機(jī)兩區(qū)域電力系統(tǒng)如圖3所示，1號(hào)發(fā)電機(jī)上安裝有廣域阻尼控制器，選擇《13作 為控制器反饋信號(hào)。廣域阻尼控制器常規(guī)超前滯后WADC，如下式所示：

(10)

[0101] 其中：Tw=10s，Ti = 0.324s，T2 = 0.212s

[0102] 在111 = 0^ = 0情況下，利用本發(fā)明所提方法得到四機(jī)兩區(qū)域電力系統(tǒng)在1=10與 Ka = 22時(shí)系統(tǒng)的穩(wěn)定裕度分別為346.4ms與101.0ms。

[0103]設(shè)置系統(tǒng)在母線3處發(fā)生三相短路故障，持續(xù)200ms。通過圖4中無時(shí)滯條件下的系 統(tǒng)響應(yīng)，可以看出廣域控制器WADC優(yōu)化了系統(tǒng)的性能，消除了內(nèi)部振蕩。圖5中給出了系統(tǒng) 中存在110ms延時(shí)情況下的系統(tǒng)響應(yīng)，可以看出延時(shí)110ms情況下，K a= 10系統(tǒng)是穩(wěn)定的，而 Ka = 22系統(tǒng)是不穩(wěn)定的，呈現(xiàn)發(fā)散趨勢(shì)。這也符合本發(fā)明所提方法求得Ka = 22條件下的時(shí)滯 穩(wěn)定裕度101 .〇ms。

【主權(quán)項(xiàng)】

1. 一種基于Wirtinger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法，其特征在于，包括如下 步驟： (1) 建立包含廣域控制回路的時(shí)滯電力系統(tǒng)模型?中：In < d (t)^h2, d{t) < μ -μ^Ι ；(2) 給定穩(wěn)定判定條件： 若存在正定矩陣P e R4nX4n;正定矩陣& e RnXn，i = 1，2，3;正定矩陣^ e RnXn，j = 1，2;矩 陣XkeR2nX2n，k=l，2使下列矩陣不等式成立，則時(shí)滯電力系統(tǒng)是漸近穩(wěn)定的：ei=[Oi-1 I 〇7-1]，1 = 1，···，7 對(duì)于矩陣A，He(A)=A+AT。其中I代表單位矩陣。 (3) 利用Matlab中的線性矩陣(LMI)工具箱判斷給定時(shí)滯d(t)是否滿足步驟(2)給出的 判定條件，若滿足，則可判定在延時(shí)d(t)條件下的時(shí)滯電力系統(tǒng)是漸近穩(wěn)定的。2. 如權(quán)利要求1所述的基于Wirtinger不等式的時(shí)滯電力系統(tǒng)穩(wěn)定性判定方法，其特征 在于，時(shí)滯電力系統(tǒng)模型式中：χ=[χι X。]1"， XI為系統(tǒng)狀態(tài)變量;X。為控制器的狀態(tài)變量;Αι為系統(tǒng) 狀態(tài)矩陣;Βι為系統(tǒng)輸入矩陣;A。為控制系統(tǒng)的狀態(tài)矩陣;C。為控制系統(tǒng)的輸出矩陣;Ci為系 統(tǒng)輸出矩陣;B。為控制系統(tǒng)的輸入矩陣。

【文檔編號(hào)】H02J3/00GK105958476SQ201610299335

【公開日】2016年9月21日

【申請(qǐng)日】2016年5月6日

【發(fā)明人】孫永輝, 李寧, 衛(wèi)志農(nóng), 孫國(guó)強(qiáng), 張世達(dá), 郭敏, 秦晨

【申請(qǐng)人】河海大學(xué)

總結(jié)

以上是生活随笔為你收集整理的时滞电力系统matlab,基于Wirtinger不等式的时滞电力系统稳定性判定方法的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。