朱良濤 程琪 馬平

摘要： 為了更好的應用電力系統(tǒng)同步相量測量裝置（phasor measurement unit，PMU）測量數(shù)據(jù)評估系統(tǒng)的電壓穩(wěn)定性，本文提出利用戴維南等值參數(shù)求取電壓穩(wěn)定裕度的改進方法。加入初值優(yōu)選環(huán)節(jié)，全面考慮戴維南等值參數(shù)和節(jié)點電壓電流不斷變化的關系，把每一時刻求出的戴維南等值參數(shù)作為下一時刻的計算初值，逐步迭代求解?；趩蝹€電壓穩(wěn)定裕度指標估算的局限性，根據(jù)系統(tǒng)本身的特性，求出了阻抗模裕度指標和角度裕度指標的加權平均值表征電壓穩(wěn)定裕度。通過對IEEE39節(jié)點系統(tǒng)的仿真分析，證明了其準確性和有效性。該研究能準確的反應系統(tǒng)的電壓穩(wěn)定性，提高表征的精度和準確性，具有一定的實際應用價值。

關鍵詞： 電壓穩(wěn)定裕度； 戴維南等值； 全微分； 初值優(yōu)選； 阻抗模裕度； 角度裕度

中圖分類號： TM712文獻標識碼： A

收稿日期： 20170602； 修回日期： 20171020

基金項目： 國家自然科學基金資助項目（51477078）

作者簡介： 朱良濤（1991），男，碩士研究生，主要研究方向為電力系統(tǒng)分析與控制。

通訊作者： 馬平（1973），女，副教授，碩士生導師，主要研究方向為電力系統(tǒng)分析與控制。Email： 814967184@qq.com近年來，電壓不穩(wěn)定因素造成的國內(nèi)外停電大事故對社會經(jīng)濟造成了巨大影響，電壓穩(wěn)定問題受到高度重視，研究電壓穩(wěn)定的方法也越來越成熟。電壓穩(wěn)定裕度作為衡量電壓穩(wěn)定程度的重要指標，主要求取方法有靈敏度法[13]、連續(xù)潮流法與遺傳算法相結合的方法[4]、相量法等，但這些方法存在誤差過大、計算局限性較大、臨界點潮流不收斂的問題。如Sui H等人[5]在已知系統(tǒng)電壓分布的條件下提出了一種根據(jù)傳輸路徑對電壓進行穩(wěn)定分析的方法；謝武忠等人[6]提出了一種非線性規(guī)劃模型來求解電壓穩(wěn)定臨界點，但受外界因素影響較大。在對電壓進行穩(wěn)定性分析過程中，戴維南等值方法由于直觀簡便，并且可準確簡化網(wǎng)絡的優(yōu)勢被眾多學者青睞。傳統(tǒng)戴維南等值方法是假設在相鄰兩個采樣點之間其等值參數(shù)不變，但在電力系統(tǒng)實際運行時，因網(wǎng)絡拓撲結構、運行模式、無功功率等因素的變化，戴維南等值參數(shù)也會改變，因此關鍵在于求取戴維南等值參數(shù)。李卓藝等人[79]通過建立與戴維南等值有關的數(shù)學模型來求解其等值參數(shù)；劉松明等人[1012]通過改進戴維南等值模型求解戴維南等值參數(shù)；李來福等人[13]對戴維南等值參數(shù)求解過程中的參數(shù)漂移問題進行了分析；湯涌等人[14]提出的基于全微分的戴維南等值參數(shù)計算方法，克服了兩個時步間戴維南等值參數(shù)求解帶來的誤差，提高了準確性和速度；韓巍等人[1720]從不同角度對電壓穩(wěn)定裕度和電壓穩(wěn)定性進行了分析判斷。因此，本文在文獻[14]的基礎上，對用戴維南等值參數(shù)求取電壓穩(wěn)定裕度的方法做了改進。文獻[14]中，在等值點負荷存在擾動的條件下，假設滿足系統(tǒng)等值電勢和阻抗不變求取等值參數(shù)的初值，但實際系統(tǒng)中很難滿足這樣的運行條件，即使?jié)M足，也很難判斷PMU實測的哪兩組數(shù)據(jù)滿足上述要求，因此通常假設一個數(shù)據(jù)窗中的前兩組滿足上述要求，求得參數(shù)的初值，但這樣的初值會嚴重影響計算結果的精度。針對該方法對初值的高度敏感性，本文考慮輸電系統(tǒng)運行時電抗遠大于電阻[15]的實際情況選取起始初值，并根據(jù)系統(tǒng)內(nèi)部等值參數(shù)和外部負荷的不斷變化關系，每一時刻的計算初值均采用前一時刻得到的等值參數(shù)，層層迭代，快速準確的求取戴維南等值參數(shù)。同時，利用求得的等值參數(shù)，采用阻抗模裕度指標和角度裕度指標的加權平均值計算電壓的穩(wěn)定裕度。通過對IEEE39節(jié)點系統(tǒng)的仿真分析，證明了其準確性和有效性。

1改進的全微分跟蹤算法

1.1傳統(tǒng)戴維南等值方法

戴維南等值就是將系統(tǒng)等值成一個電壓源和一個阻抗串聯(lián)的兩節(jié)點系統(tǒng)，戴維南等值系統(tǒng)如圖1所示。

圖1戴維南等值系統(tǒng)圖1中，k和Zk表示戴維南等值電勢和等值阻抗；k和k表示負荷母線的電壓和電流，可通過PMU裝置測得。其中，k=Zkk+k，若假設k和Zk在相鄰的兩個時刻不變，通過采樣點k和k+1的數(shù)據(jù)，可得方程組

k=Zkk+k

k=Zkk+1+k+1（1）

由式（1）解得k和Zk的值，對戴維南等值參數(shù)進行實時跟蹤求解。

1.2改進的全微分跟蹤算法

文獻[14]所述的全微分方法，由圖1根據(jù)潮流方程求得戴維南等值參數(shù)及負荷母線電壓幅值，并取全微分，聯(lián)立可得

U2k-UkEkr+PkRk+QkXk=0

UkEki-PkXk+QkRk=0

U2k+1-Uk+1Ek+1r+Pk+1Rk+1+Qk+1Xk+1=0

Uk+1Ek+1i-Pk+1Xk+1+Qk+1Rk+1=0

Pk+1-Pk=fP（Uk，Uk+1，Rk，Xk，Rk+1，Xk+1，Ekr，Eki，Ek+1r，Ek+1i）

Qk+1-Qk=fQ（Uk，Uk+1，Rk，Xk，Rk+1，Xk+1，Ekr，Eki，Ek+1r，Ek+1i）（2）

方程組（2）包括6個方程和8個未知數(shù)，依據(jù)傳統(tǒng)戴維南等值方法，假設在等值點負荷擾動時，系統(tǒng)等值電勢和阻抗不變，利用式（1）求得其中兩個未知數(shù)的初值，即R0+jX0初值，此初值的精度會高度影響后續(xù)時刻等值參數(shù)的精度。

在參數(shù)求取過程中，本文不設參考相位，直接運用真實值進行推導。初值的選取很重要，為減小初值對后續(xù)時刻等值參數(shù)的影響，本文對該問題做深入探討。采用PMU測量數(shù)據(jù)，對于節(jié)點k，有

k=Zkk+k（3）

其中，=Er+jEi；Z=R+jX；=Ir+jIi和=Ur+jUi。

分離式（3），可得

Er=RIr-XIi+Ur（4）

Ei=RIi+XIr+Ui（5）

整理得

Ir=REr+XEi-UiX-UrRR2+X2（6）

Ii=REi-XEr+UrX-UiRR2+X2（7）

式（6）和式（7）中，Ir、Ii對Er、Ei、R、X和U求全微分，整理得

dIr=IrErdEr+IrEidEi+IrRdR+IrXdX+IrUrdUr+IrUidUi（8）

dIi=IiErdEr+IiEidEi+IiRdR+IiXdX+IiUrdUr+IiUidUi（9）

在每一步電氣量合理波動情況下，式（8）和式（9）可表示為差分方程

Ir（k+1）-Irk=RkR2k+X2k（Er（k+1）-Erk）+XkR2k+X2k×（Ei（k+1）-Eik）+（（Erk-Urk）（R2k+X2k）R2k+X2k2-

2Rk（RkErk+XkEik-UikXk-UrkRk）R2k+X2k2）×（R（k+1）-Rk）+（（Eik-Uik）（R2k+X2k）R2k+X2k2-

2Xk（RkErk+XkEik-UikXk-UrkRk）R2k+X2k2）×（X（k+1）-Xk）-RkR2k+X2k（Ur（k+1）-Urk）-

XkR2k+X2k（Ui（k+1）-Uik）=F1（10）

Ii（k+1）-Iik=-XkR2k+X2k（Er（k+1）-Erk）+RkR2k+X2k×（Ei（k+1）-Eik）+（（Eik-Uik）（R2k+X2k）R2k+X2k2-

2Rk（RkEik-XkErk+UrkXk-UikRk）R2k+X2k2）×（R（k+1）-Rk）+（（Urk-Erk）（R2k+X2k）R2k+X2k2-

2Xk（-XkErk+RkEik+UrkXk-UikRk）R2k+X2k2）×（X（k+1）-Xk）+XkR2k+X2k（Ur（k+1）-Urk）-

RkR2k+X2k（Ui（k+1）-Uik）=F2（11）

由式（4）和式（5）及式（10）和式（11）得方程組為

Er（k+1）=R（k+1）Ir（k+1）-X（k+1）Ii（k+1）+Ur（k+1）

Ei（k+1）=R（k+1）Ii（k+1）-X（k+1）Ir（k+1）+Ui（k+1）

Ir（k+1）-Irk=F1

Ii（k+1）-Iik=F2（12）

方程組（12）包含8個未知數(shù)和4個方程，須給定其中4個未知數(shù)的值或再列寫4個方程才能求解。本文采用第1種方法，即先給定其中4個未知數(shù)的值。為避免傳統(tǒng)的戴維南等值方法存在的參數(shù)漂移問題，在一個數(shù)據(jù)窗內(nèi)，經(jīng)過初值優(yōu)選得到4個未知數(shù)的起始初值Er0、Ei0、R0、X0，利用起始初值Er0、Ei0、R0、X0和測量數(shù)據(jù)Ui1、Ur1、Ii1、Ir1通過式（12）求取一組戴維南等值參數(shù)Er1、Ei1、R1、X1，進而把Er1、Ei1、R1、X1作為下一組測量數(shù)據(jù)的計算初值計算Ei2、Er2、R2、X2。以此類推，該時刻得到的等值參數(shù)作為求解下一時刻的初值，逐步迭代計算，充分考慮兩個時刻間每一個參數(shù)的變化。

1.3戴維南等值參數(shù)的初值優(yōu)選

在實際電力系統(tǒng)中，負荷變化會導致母線電壓的變化，戴維南等值阻抗和等值電勢也會發(fā)生變化。但如果一個數(shù)據(jù)窗內(nèi)時間相隔太近，很容易出現(xiàn)數(shù)據(jù)參數(shù)漂移問題，本文對采樣數(shù)據(jù)進行優(yōu)選，從而得到理想初值。

1）在同一數(shù)據(jù)窗口進行采樣時，電壓和電流的幅值差滿足一定條件，即

（k+1）-k>γ1Uk

（k+1）-k>γ2Ik（13）

式中，γ1和γ2由系統(tǒng)本身決定。

2）由條件1）篩選的兩組數(shù)據(jù)，得到的等值阻抗要滿足電抗遠大于電阻，即X>>R。等值阻抗為

Zk=Rk+jXk=-（k+1）-k（k+1）-k=-ΔkΔk=Zk∠φZk=Zkcos φZk+jZksin φZk（14）

其中，φZk=-（Δφuk-Δφik）為等值阻抗角。若X>>R，有

sin φZkcos φZk>>1（15）

將經(jīng)過優(yōu)選得到的兩個采樣點數(shù)據(jù)代入方程組（1），即可求得優(yōu)化后的起始初值。

2靜態(tài)電壓穩(wěn)定裕度指標

2.1阻抗模裕度指標

在圖1所示系統(tǒng)中，靜態(tài)電壓穩(wěn)定臨界點條件是Zeq=Zl，電壓穩(wěn)定裕度指標用阻抗模裕度表示為

MZ=Zl-ZeqZl（16）

當MZ>0時，系統(tǒng)運行在穩(wěn)定域；當MZ=0時，達到電壓穩(wěn)定臨界點。說明MZ的值越小，表征裕度越低。

2.2角度表示的電壓穩(wěn)定裕度指標

在靜態(tài)電壓穩(wěn)定臨界點處，圖1所示系統(tǒng)兩端電壓相角差滿足[16]

δcr=-05Φ+05α， -α≤Φ≤π/2

05Φ+15α， -π/2≤Φ≤-α（17）

其中，Φ為負荷功率因數(shù)角；α為等值阻抗角；δcr為節(jié)點臨界電壓相角。角度表示的電壓穩(wěn)定裕度為

Mδ=δcr-δδcr（18）

其中，δ=δ1-δ2。Mδ=0即為電壓穩(wěn)定臨界點，電壓穩(wěn)定裕度隨Mδ同向變化。

2.3靜態(tài)電壓平均穩(wěn)定裕度

阻抗模裕度指標和角度表示的電壓穩(wěn)定裕度指標是以阻抗?；蚪嵌确矫娴膯蝹€指標來衡量穩(wěn)定裕度，M=σ1MZ+σ2Mδ（19）

其中，σ1和σ2為權重系數(shù)，其值根據(jù)系統(tǒng)自身特性選取。

3算例分析

本文對IEEE39節(jié)點系統(tǒng)的節(jié)點8進行分析，設該節(jié)點裝有PMU測量裝置。一個數(shù)據(jù)窗口有8個采樣點，采樣點數(shù)據(jù)如表1所示。

運用13的優(yōu)選方法對測量數(shù)據(jù)進行優(yōu)選，其中γ1=0001，γ2=001。對采樣點6和7進行優(yōu)選，優(yōu)化初值和數(shù)據(jù)窗內(nèi)第1和第2組數(shù)據(jù)計算的未優(yōu)選戴維南參數(shù)的起始初值如表2所示。

分別以表2中的兩組初值作為表1中采樣點1對應的戴維南等值參數(shù)，應用式（12）迭代，求出采樣點2～8的戴維南等值參數(shù)，戴維南等值參數(shù)如表3所示。

由表2和表3可以看出，優(yōu)選的初值與電力系統(tǒng)的實際運行情況相符，滿足X>>R，且電阻為正數(shù)；未優(yōu)選初值中電阻為負數(shù)，這與系統(tǒng)的實際情況嚴重不符；采用優(yōu)選的初值計算出的后續(xù)采樣點對應的阻抗參數(shù)，即穩(wěn)定也都符合電力系統(tǒng)的實際情況；采用未優(yōu)選的初值計算出的后續(xù)采樣點的阻抗參數(shù)中，電阻均為負數(shù)，嚴重偏離系統(tǒng)的實際情況，表明初值的優(yōu)選對全微分方法非常重要，初值選取的合理是該方法有效性的重要保障。

根據(jù)系統(tǒng)特性，本文權重系數(shù)選取σ1=072，σ2=028，求解電壓穩(wěn)定裕度指標如表4所示。由表4可以看出，電壓穩(wěn)定裕度趨于穩(wěn)定。因此，本文方法具有有效性。

4結束語

本文利用改進的戴維南等值方法對要研究區(qū)域進行等值分析，并根據(jù)求得的等值參數(shù)計算電壓穩(wěn)定裕度。用戴維南等值方法對電力系統(tǒng)進行等值關鍵在于求取其等值參數(shù)。本文對原有全微分方法做了改進，通過初值優(yōu)選得到起始初值，而后每一時刻的計算初值都采用上一時刻得到的等值參數(shù)，全面考慮到戴維南等值參數(shù)的幅值和相角隨著時間變化的實際情況，減小了計算時間和計算誤差，通過檢驗證明了其有效性?；趩蝹€電壓穩(wěn)定裕度指標估算的局限性，根據(jù)系統(tǒng)本身的特性，本文以阻抗模裕度指標和角度裕度指標的加權平均值來表征電壓穩(wěn)定裕度，所得結果更能準確的反應系統(tǒng)的電壓穩(wěn)定性，提高了表征的精度和準確性。

參考文獻：

[1]李國慶， 姜濤， 徐秋蒙， 等. 基于局部電壓穩(wěn)定指標的裕度靈敏度分析及應用[J]. 電力自動化設備， 2012， 32（4）： 15.

[2]代飛， 崔挺， 徐箭， 等. 基于綜合靈敏度分析的電壓校正控制[J]. 電力自動化設備， 2011， 31（12）： 1520.

[3]張小兵， 吳政球， 李連偉， 等. 基于擬合方法的N1靜態(tài)電壓穩(wěn)定裕度計算[J]. 電力系統(tǒng)保護與控制， 2010， 38（14）： 2327.

[4]李娟， 陳繼軍， 司雙. 連續(xù)潮流與免疫遺傳算法結合的靜態(tài)電壓穩(wěn)定裕度計算[J]. 電力系統(tǒng)保護與控制， 2010， 38（18）： 2427， 32.

[5]Sui H， Zhao J. Voltage Stability Index Based on Synchronized Phasor Measurement[J]. Electric Power Automation Equipment， 2008， 28（3）： 59.

[6]謝武忠， 胡剛， 劉艷， 等. 一種計算電力系統(tǒng)電壓穩(wěn)定臨界點新模型[J]. 中國電力， 2011， 44（8）： 14.

[7]李卓藝， 劉俊勇， 劉友波， 等. 系統(tǒng)持續(xù)擾動下的電力系統(tǒng)戴維南等值參數(shù)在線識別[J]. 中國電機工程學報， 2015， 35（8）： 19001908.

[8]李連偉， 吳政球， 鐘浩. 基于節(jié)點戴維南等值的靜態(tài)電壓穩(wěn)定裕度快速求解[J]. 中國電機工程學報， 2010， 30（4）： 7983.

[9]付正鑫， 盧繼平， 喬梁， 等. 基于戴維南參數(shù)等值的電壓穩(wěn)定性分析[J]. 中國電力， 2014， 47（5）： 4447.

[10]劉明松， 張伯明， 姚良忠， 等. 基于PMU和改進戴維南等值模型的電壓穩(wěn)定在線監(jiān)視[J]. 電力系統(tǒng)自動化， 2009， 33（10）： 610.

[11]李衛(wèi)星， 牟曉明， 李志民. 電力系統(tǒng)戴維南等值參數(shù)的解析與思考[J]. 中國電機工程學報， 2012， 32（S1）： 2832.

[12]牟善科， 丁濤， 顧偉， 等. 基于偏差校正的戴維南等值參數(shù)在線跟蹤改進算法[J]. 電力系統(tǒng)保護與控制， 2011， 39（2）： 2328.

[13]李來福， 于繼來， 柳焯. 戴維南等值跟蹤的參數(shù)漂移問題研究[J]. 中國電機工程學報， 2005， 25（20）： 15.

[14]湯涌， 孫華東， 易俊， 等. 基于全微分的戴維南等值參數(shù)跟蹤算法[J]. 中國電機工程學報， 2009， 29（13）： 4853.

[15]Corsi S， Taranto G N. A RealTime Voltage Instability Identification Algorithm Based on Local Phasor Measurements[J]. IEEE Transactions on Power Systems， 2008， 23（3）： 12711279.

[16]劉道偉， 韓學山， 王勇， 等. 在線電力系統(tǒng)靜態(tài)穩(wěn)定域的研究及應用[J]. 中國電機工程學報， 2009， 29（34）： 4249.

[17]韓巍， 于汀， 薛安成， 等. 采用關鍵PMU量測的含風電系統(tǒng)電壓穩(wěn)定快速評估方法[J]. 電網(wǎng)技術， 2014， 38（s1）： 610.

[18]Mousavi O A， Bozorg M， Cherkaoui R. Preventive Reactive Power Management for Improving Voltage Stability Margin[J]. Electric Power Systems Research， 2013， 96（1）： 3646.

[19]Hu F， Sun K， Rosso A D， et al. MeasurementBased RealTime Voltage Stability Monitoring for Load Areas[J]. IEEE Transactions on Power Systems， 2016， 31（4）： 27872798.

[20]劉明波， 袁康龍， 楊柳青， 等. 兩種計算靜態(tài)電壓穩(wěn)定裕度方法的比較[J]. 華南理工大學學報：自然科學版， 2013， 41（4）： 17.