，，，，，

(廣西大學 電氣工程學院,廣西壯族自治區(qū) 南寧 530004)

1 引言

暫態(tài)穩(wěn)定性分析計算是電力系統(tǒng)三大經(jīng)典計算質疑，它是電力系統(tǒng)規(guī)劃、設計、調度運行與控制必不可少的一項重要的計算任務。在正常的穩(wěn)態(tài)運行情況下，各發(fā)電機輸出的電磁轉矩與原動機輸入的機械轉矩之間相互平衡，因而所有發(fā)電機的轉子速度能一直維持恒定。但是，電力系統(tǒng)中不可避免地受到一些大擾動的沖擊，例如，各類短路故障、大容量機組或大負荷的投入和切除等。電力系統(tǒng)的暫態(tài)穩(wěn)定性分析，其主要目的是在大擾動的沖擊時，各發(fā)電機維持同步運行的能力。

電力系統(tǒng)暫態(tài)穩(wěn)定分析方法可以分為：時域仿真法、直接法、人工智能法、動態(tài)安全域方法。把以上方法當中的一種或者幾種方法相互結合形成的混合法等。近年來，有關時域仿真算法許多的改進算法，主要目的是加快仿真速度或者提高計算精度。在數(shù)學上，暫態(tài)穩(wěn)定時域仿真[1-2]可以表述為微分代數(shù)方程的初值問題。有聯(lián)立求解和交替求解這兩類方法。交替求解法是指微分方程組和代數(shù)方程組兩者的求解彼此是獨立但需要交替進行的。但會造成“交接誤差”。由于顯式數(shù)值積分方法的數(shù)值穩(wěn)定性不太好，所以顯式分離求解法的整體數(shù)值穩(wěn)定性也不是很好。因此顯式分離求解法在商業(yè)化的電力系統(tǒng)暫態(tài)穩(wěn)定性計算程序中應用很少。聯(lián)立求解法的基本過程為，先用隱式梯形積分公式把微分方程組進行差分化，它和代數(shù)方程組一起形成聯(lián)立非線性方程組，然后求解此非線性方程組。所以，仿真算法的研究可以分為兩種，一種是常微分數(shù)值積分方法，另一種是對非線性方程組的求解方法。常微分數(shù)值積分方法有改進歐拉法、隱式積分法、龍格-庫塔法等，在此基礎上也有改進方法。文獻[3-5]研究了在電力系統(tǒng)中，中長期的暫態(tài)穩(wěn)定性分析以及變步長變階方法。文獻[6]研究了不誠實牛頓法和變步長技術相結合的方法，提高了系統(tǒng)的仿真速度，驗證算法的正確性和實用性。非線性方程組的求解方法一般用牛頓法迭代求解。文獻[7]基于奇異攝動法，提出求解發(fā)電機狀態(tài)變量和網(wǎng)絡的節(jié)點電壓的方法，計算量比牛頓迭代法要小，從而提高仿真速度。文獻[8]通過引入新的狀態(tài)變量，網(wǎng)絡代數(shù)方程的形式發(fā)生了改變，計算過程中沒有迭代和三角分解，大幅度的提高了計算速度。

本文應用攝動法進行電力系統(tǒng)暫態(tài)計算仿真研究。結合IEEE9及節(jié)點系統(tǒng)，從發(fā)電機的轉速、機端電壓、和仿真時間三個方面。驗證該方法的有效性和實用性。

2 暫態(tài)穩(wěn)定仿真的奇異攝動法

進行電力系統(tǒng)動態(tài)研究時，發(fā)電機和電力系統(tǒng)輸電網(wǎng)絡通?？梢员磉_為：

(1)

0=g(x,U)

(2)

式中：x是同步發(fā)電機的狀態(tài)變量;U是網(wǎng)絡的電壓向量。

采用隱式積分的方法求解微分方程時，一般先把微分方程換成一系列代數(shù)方程，再用求解代數(shù)方程的辦法求出時段終值。隱式梯形積分法的差分方程可以表達為：

(3)

因此，描述電力系統(tǒng)的動態(tài)方程(1)、(2)可以改寫為：

xk+1=xk+εf(xk,Uk)+εf(xk+1,Uk+1)

(4)

0=g(x,U)

(5)

式中：ε=h/2，h為步長。

為簡單起見，式(3)、(4)可寫為：

x=B+εf(x,U)

(6)

0=g(x,U)

(7)

其中，設x=xk+1,U=Uk+1，B=xk+εf(xk,Uk)很顯然在式(4)中ε是一個很小的常數(shù)，可以應用奇異攝動原理進行電力系統(tǒng)仿真計算，把x和U分別展開為ε的冪級數(shù)，令

(8)

(9)

ε0:x0=B

(10)

0=g(x0,U0)

(11)

ε1:x1=f(x0,U0)

(12)

(13)

(14)

(15)

式中：xx1，xU1和UU1都是向量，每一個元素均可由向量x1,U1解出。

因此，x0,U0可由(10)、(11)中求出，x1,U1可由(12)、(13)求出，x2,U2可由(14)、(15)求出。最后x和U可以表達為：

x≈x0+εx1+ε2x2

(16)

U≈U0+εU1+ε2U2

(17)

在計算中如果保留的次冪越高,精度也就會越高，但計算就會隨之增大。實際計算表明，在一般的電力系統(tǒng)仿真計算中，保留前三項就能夠滿足精度要求。如果是在精度要求高時，則需要以保留更多的項來滿足精度要求。

針對x0=xk+εf(xk,Uk)，x1=f(xk,Uk)，則可以得出：x0=xk+εx1

x=x0+εx1+ε2x2+o(ε3)=xk+2εx1+o(ε3)

(18)

U=U0+εU1+ε2U2+o(ε3)=Uk+2εU1+o(ε3)

(19)

綜上可以得出：

ε0:x0=B,0=g(x0,U0)

(20)

(21)

從而可以得出：

x=xk+2εx1,U=Uk+2εU1

(22)

因此，若想求解x、U，只需求解出x0、U0和x1、U1，進一步減少了計算量。

3 電力系統(tǒng)暫態(tài)計算的數(shù)學模型

在電力系統(tǒng)暫態(tài)穩(wěn)定計算中，當選用的元件的數(shù)學模型不同，所得出的計算結果會受直接影響，穩(wěn)定性分析計算的復雜程度也會有很大的影響。因此，選擇合適的數(shù)學模型來表示元件本身，不僅計算過程簡單還可以得出滿足精度要求結果。數(shù)學模型是電力系統(tǒng)暫態(tài)計算中的一個重要的問題。在工程實際問題中，根據(jù)使用場合的不同對發(fā)電機的模型做不同程度的簡化處理。本文選用發(fā)電機三階模型，負荷為恒定阻抗模型。

3.1 發(fā)電機三階實用模型

在實際電力系統(tǒng)，暫態(tài)計算仿真分析中，當要計及勵磁系統(tǒng)的作用時，能選擇的最簡單的模型就是三階模型。它較適用于凸極機。這種實用模型基于如下假定：

(1)不考慮定子d繞組、q繞組的暫態(tài)；

(2)在定子電壓方程中，設在速度變化不大的過渡過程中，其引起的誤差很小；

(3)忽略D繞組、Q繞組，通過補入阻尼項到轉子運動方程中來近似處理。

(23)

3.2 電力網(wǎng)絡數(shù)學模型

電力系統(tǒng)暫態(tài)計算和潮流、短路計算中一樣，用相量把電力網(wǎng)絡的節(jié)點電壓方程可表示成：

YV=I

(24)

式中：分別為電力網(wǎng)絡的節(jié)點電流和節(jié)點電壓組成的列向量；Y表示節(jié)點導納矩陣。上式所描述的網(wǎng)絡方程在形式上為線性方程組，其中的導納矩陣Y僅由電力網(wǎng)絡的結構和參數(shù)所決定。

圖1 第i臺發(fā)電機的向量圖

在電力系統(tǒng)中，為了求解電力網(wǎng)絡方程，需要首先列寫出每一個動態(tài)元件的本身的代數(shù)方程，并對其進行處理，從而才能和網(wǎng)絡方程聯(lián)立求解。

3.3 發(fā)電機與網(wǎng)絡的連接

對于三階同步發(fā)電機模型，在d-q坐標系下，定子電壓方程都可統(tǒng)一表示為：

(25)

4 算例分析

本文的采用電科院3機9節(jié)點測試系統(tǒng)。仿真時間為5 s，設定故障為7號節(jié)點，0.1 s時發(fā)生三相短路故障，步長為0.001 s時，0.15 s時切出故障。

4.1 仿真精度對比分析

如圖2所示，表示利用牛頓法求解時，1、2號發(fā)電機電壓仿真曲線。圖3表示利用攝動法求解時，1、2號發(fā)電機電壓仿真曲線。

由圖2、3可以看出，用兩種方法進行仿真時，攝動法和牛頓法的計算結果曲線基本相同。所以，從這個意義上來說，驗證了它們的精度是相同的。

圖2 牛頓法求解時，1、2號發(fā)電機電壓仿真曲線

圖3 攝動法求解時，1、2號發(fā)電機電壓仿真曲線

表1 不同算法的仿真時間對比

4.2 仿真時間對比分析

由表1可以看出，當系統(tǒng)總仿真時間都為5s，牛頓法計算時間是0.462s，攝動法為0.075s，計算速度提高6倍。證明了攝動法的優(yōu)越性。

5 結論

本文利用攝動法進行電力系統(tǒng)暫態(tài)計算，其計算速度比牛頓法要快，通過仿真，利用算例驗證該方法在保持計算精度的情況下，而提高計算速度，取得較好的計算效果。