趙 淵 吳小平 謝開貴

（1. 重慶大學輸配電裝備及系統(tǒng)安全與新技術國家重點實驗室 重慶 400044 2. 廣東電網公司佛山供電局 佛山 528000）

1 引言

近年來，國內外電網多次發(fā)生因遭受嚴重有功缺額致使頻率快速下跌，甚至導致頻率崩潰的嚴重電力事故[1]，系統(tǒng)保持頻率穩(wěn)定的能力正面臨嚴峻挑戰(zhàn)。擾動事故下頻率的動態(tài)行為及其穩(wěn)定一直是穩(wěn)定研究的薄弱環(huán)節(jié)，直至2004年IEEE/CIGRE穩(wěn)定定義聯(lián)合工作組才對頻率穩(wěn)定給出了定義：頻率穩(wěn)定是指電力系統(tǒng)發(fā)生突然的有功缺額后，發(fā)電和負荷需求出現(xiàn)不平衡，系統(tǒng)頻率保持或恢復到允許范圍內而不發(fā)生頻率崩潰的能力[2,3]。

目前頻率穩(wěn)定研究存在以下不足：①對頻率穩(wěn)定中的機-網協(xié)調問題缺乏足夠重視。當?shù)皖l運行時間超過機組低頻保護整定值時，機組被切除，功率缺額的加大將使頻率下降更嚴重，有可能形成連鎖切機，導致頻率崩潰，1999年臺灣大停電便是機-網嚴重不協(xié)調的典型例子。②在頻率動態(tài)計算和穩(wěn)定分析中缺乏旋轉備用影響的量化分析。若能在低頻減載配置和機組異常頻率保護中考慮旋轉備用的作用，則能減少過控現(xiàn)象。③缺乏從頻率動態(tài)角度來評價頻率的穩(wěn)定性。對電力系統(tǒng)頻率穩(wěn)定研究最有價值的頻率數(shù)據并非系統(tǒng)穩(wěn)態(tài)頻率，而是系統(tǒng)遭受擾動時的動態(tài)頻率，文獻[4]提出的頻率穩(wěn)定計算方法只能得到系統(tǒng)穩(wěn)態(tài)頻率，而不能求得頻率動態(tài)過程中的最低頻率或某低頻值下的運行時間等運行人員關注的動態(tài)信息。④缺乏有效方法和概率指標來計及不確定因素對頻率穩(wěn)定的影響。正確考慮系統(tǒng)運行工況的不確定性，可減少確定性方法造成的欠控或過控損失[5]，概率指標可從概率統(tǒng)計角度對系統(tǒng)保持頻率穩(wěn)定的能力進行定量描述。因此將可靠性評估方法引入頻率穩(wěn)定研究中以全面客觀地反映頻率穩(wěn)定的安全裕度水平已成為迫切需要。

基于上述不足，本文推導了計及旋轉備用影響的頻率動態(tài)特性，描述了低頻減載作用下頻率穩(wěn)定的分析過程，并判斷機組低頻保護是否動作；為實現(xiàn)頻率穩(wěn)定的風險評估，對發(fā)電方式、網絡拓撲和負荷水平分別建立了隨機模型；構建了評價頻率穩(wěn)定可靠性的指標體系，最后基于非序貫蒙特卡洛仿真對IEEE RTS79測試系統(tǒng)進行了驗算。

2 基于頻率動態(tài)特性的頻率穩(wěn)定分析模型

頻率動態(tài)特性是指頻率在擾動事故下隨時間的動態(tài)變化過程。在動態(tài)過程中，頻率或會重新恢復至期望值，或發(fā)生失穩(wěn)，因此可根據頻率動態(tài)過程來判斷頻率的穩(wěn)定性。由于頻率穩(wěn)定重點關注全網頻率的平均變化過程，而單機模型能較好反映系統(tǒng)平均頻率的動態(tài)過程，文獻[6-8]指出單機模型的計算結果與實際系統(tǒng)記錄較為吻合，作為頻率動態(tài)特性分析和低頻減載配置依據是適宜的。

2.1 考慮旋轉備用的電力系統(tǒng)頻率動態(tài)特性

標幺值下的系統(tǒng)運動方程為

由于 dω*/dt=dΔω*/dt=dΔf*/dt，并考慮到發(fā)電機組慣性較大，一般電角速度的變化和轉矩不是太大，則轉矩標幺值可近似認為等于功率標幺值[9]，即Tm*-Te*=Pm*-Pe*，單機帶綜合負荷時，發(fā)電機輸出電磁功率 Pe*即所帶負荷有功 PL*，故式（1）可寫成

式中，TJ為系統(tǒng)慣性時間常數(shù)；Pm*、PL*分別為發(fā)電機機械輸入功率和負荷有功功率標幺值；Δf*為以系統(tǒng)額定頻率fN為基準的頻率偏差標幺值。

負荷功率PL*可表示為

式中，kL*為負荷的功率頻率調節(jié)效應系數(shù)。

原動機調速器動作有一定延時，調節(jié)閥門開度以增加機械功率也需要時間，故ΔPm*滯后于Δf*的變化，兩者之間為時間函數(shù)關系。穩(wěn)定后，發(fā)電機機械功率變化對頻率變化之比為機組的標幺單位調節(jié)功率kG*。為簡化計算，ΔPm*滯后于Δf*的關系可用一階慣性環(huán)節(jié)表示

式中，TR為慣性環(huán)節(jié)的時間常數(shù)，則式（2）可寫為

令ΔP0*=Pm0*－PL0*為故障瞬間初始功率缺額標幺值。對式（4）和式（5）進行頻域分析，可得Δf*的頻域表達式為

從而Δf*的時域表達式為

考慮備用時頻率動態(tài)特性有名值表達式為

式中，f0為發(fā)生功率缺額瞬間時的初始頻率，本文假設發(fā)生功率缺額瞬間，初始頻率等于額定頻率fN；ks*為系統(tǒng)功率頻率調節(jié)效應系數(shù)。

由于機組出力變化滯后系統(tǒng)頻率變化，即故障瞬間機組出力不變，則故障瞬間有功缺額可歸納為兩種情況[10]：①僅有發(fā)電機故障引起的功率缺額：此時運行機組丟失的出力即是故障瞬間的功率缺額；②輸電線路故障引起的功率缺額：某些重要輸電線路故障，可能導致發(fā)電廠出現(xiàn)“窩電”現(xiàn)象，例如只有單回送出線路的發(fā)電廠，線路故障導致機組出力完全受阻，而有多回送出線路的發(fā)電廠，部分線路故障導致穩(wěn)定極限降低，發(fā)電廠出力受限。此時的功率缺額等于受阻或受限機組為保證暫態(tài)和靜態(tài)穩(wěn)定而減少的出力。

2.2 低頻減載作用下的頻率穩(wěn)定分析模型

低頻減載要保證各輪動作的選擇性，即前一輪動作后如果不能阻止頻率下降，后一輪才動作[11]，動作原理如圖 1，圖中 fqi(i=1,2,…,n)為第 i輪的整定啟動頻率，第 i輪在 tqi時刻啟動，經過時延 tdi后在ti時刻動作，ti時刻的頻率用f (ti)表示。

頻率變化過程可逐段表示，設ti-1時第i-1輪動作切除部分負荷，則ti-1至ti時段內頻率變化公式為

圖1 低頻減載裝置的工作原理Fig.1 Principle of under frequency load shedding

如第i-1輪動作后，功率缺額完全得到補償，即Pmi-1*-PLi-1*＞0，則系統(tǒng)頻率將恢復。

計算各輪動作時間是分析頻率動態(tài)過程的關鍵，文獻[8,12]采用了近似方法，即認為頻率的變化率在研究的各時段內不變，實際上頻率的變化率隨時間在不斷變化，如圖2所示，第i-1輪動作后，如仍不能制止頻率下降，則應計算第 i輪啟動時間tqi，近似計算的第 i輪啟動時間為qit′，比準確時間tqi提前，當頻率下降速度較慢時，將產生較大誤差。由于很難從式（11）得到當頻率變化Δfd*時所需時間Δt的精確表達式，考慮到功率缺額較大時，低頻減載各輪動作很快，此時機組出力變化不大，所以可用無備用時的頻率動態(tài)特性來計算Δt。

圖2 低頻減載各輪動作時間計算示意圖Fig.2 Action time for under frequency load shedding

令式（5）中ΔPm*(t)=0，可得無備用下頻率偏差標幺值表達式

則第i輪動作時間ti的表達式為

式中，Δtqi為第i-1輪動作時間ti-1至第i輪啟動時間tqi的時間間隔；fqi為第i輪的啟動頻率；tdi為第i輪的動作時延；PLi-1*為第 i-1輪動作切負荷后剩余的負荷功率，計算如式（16）所示。

式（16）右邊第二項是因負荷自身功率頻率調節(jié)效應而減少的負荷有功；Δf*(ti-1)為第 i-1輪在 ti-1時刻切負荷時的頻率偏差，Δf*(ti-1)=[ f(ti-1)-f0]/fN；Cm*為第m輪切除的負荷量，為簡化計算模型，式（16）中未考慮已切負荷隨頻率的變化。需指出的是，若低頻減載各輪動作后的機組機械功率 Pmi*和負荷功率PLi*不易求得，可用頻率動態(tài)過程另一計算公式[13]：

式中，Ci*為第i輪切除的負荷量。

在頻率恢復或切負荷過程中，若出現(xiàn)頻率懸浮現(xiàn)象，則低頻減載特殊輪將再切除部分負荷，此時的頻率動態(tài)計算與基本輪類似；若機組低頻運行時間超出低頻保護整定值，則切除該機組，如切機時低頻減載已動作至第i輪（i＜n），則需判斷切機后第i+1～n輪是否啟動，計算過程與基本輪類似。

2.3 各抽樣系統(tǒng)狀態(tài)下計算參數(shù)的等值

用指數(shù)分布模擬發(fā)電機和線路運行狀態(tài)的不確定性，且只考慮元件正常和停運兩個狀態(tài)，使用非序貫蒙特卡洛抽樣得到元件和系統(tǒng)狀態(tài)[14]。

取故障后瞬間運行發(fā)電機總有功出力Si為功率基準，其下標 i表示第 i個系統(tǒng)狀態(tài)。系統(tǒng)等效慣性時間常數(shù)TJ、負荷等效頻率調節(jié)效應系數(shù)kL*和發(fā)電機等效單位調節(jié)功率kG*按以下各式分別計算。

式中，K為系統(tǒng)狀態(tài) i下處于運行狀態(tài)的發(fā)電機總數(shù)；MVAk為第 k臺發(fā)電機的額定容量；TJ0k為第 k臺發(fā)電機以自身額定容量為基準的慣性常數(shù)；kL0*為以系統(tǒng)額定負荷 PLN為基準的負荷頻率調節(jié)系數(shù)；R為運行中具有旋轉備用的機組總數(shù)；kG0r*為第r臺發(fā)電機以自身額定容量為基準的單位調節(jié)功率。

低頻減載的整定是事先人為假定一個比較嚴重的功率缺額，如功率缺額為運行發(fā)電機總有功出力的40%。在分析系統(tǒng)狀態(tài)i時，由于出現(xiàn)了機組故障，因此系統(tǒng)狀態(tài) i下的運行發(fā)電機總有功出力將小于基態(tài)下的有功出力，在頻率動態(tài)分析時，低頻減載的整定值需在新的功率基準下重新計算標幺值，這時需將各輪整定值以系統(tǒng)狀態(tài)i下的Si為功率基準進行折算，如式（21）。

3 頻率穩(wěn)定可靠性評估指標體系

本文從頻率失穩(wěn)指標、期望切負荷量指標和頻率動態(tài)指標三大類構建了如下風險指標。

3.1 頻率失穩(wěn)指標

（1）頻率失穩(wěn)概率。頻率穩(wěn)定性還沒有嚴格判據，目前可使用兩種判別方法：①若頻率低于某一低頻值的時間超過規(guī)定時間，如低于47.5Hz的時間超過5s，即可判別頻率失穩(wěn)。②若頻率下降至某一低頻值以下，則判別頻率失穩(wěn)，如低于 47Hz時系統(tǒng)頻率基本上已無法挽回，終將導致頻率失穩(wěn)。本文采用第二種方法。

式中，F(xiàn)為系統(tǒng)頻率失穩(wěn)的故障狀態(tài)集合；Pi為系統(tǒng)狀態(tài)i發(fā)生的概率。

（2）累計頻率失穩(wěn)次數(shù)

式中，fi為系統(tǒng)進入失穩(wěn)狀態(tài) i的年平均次數(shù)；vk為第k個元件從狀態(tài)i離開的轉移率，如其處于工作狀態(tài)，則vk取失效率，否則取修復率；N為元件總數(shù)。

（3）頻率失穩(wěn)時間的概率分布。頻率失穩(wěn)時間是指故障瞬間至頻率下降至最低允許頻率所經歷的時間。此處定義的概率分布是離散概率分布。對失穩(wěn)時間概率分布計算說明如下：首先將失穩(wěn)時間按一定時間間隔平均分為多個區(qū)間，然后對每一個導致頻率失穩(wěn)的嚴重故障狀態(tài)分別計算其失穩(wěn)時間，并將此失穩(wěn)時間及該狀態(tài)發(fā)生的概率歸類至某一個區(qū)間。下文中其他指標的概率分布計算思路與此類似。

3.2 各輪期望切負荷量指標

期望切負荷量是低頻減載裝置為防止頻率快速下降而動作切除的負荷，計算公式為

式中，n為低頻減載的總輪數(shù)；Lcm為第m輪的期望負荷削減量，MW；S為造成系統(tǒng)低頻切負荷的故障集合；Pi為系統(tǒng)狀態(tài)i出現(xiàn)的概率；為系統(tǒng)狀態(tài)i下第m輪是否動作的判斷標志，若第m輪動作，則為1，反之為0；為狀態(tài)i下第m輪負荷切除量的標幺值。

3.3 頻率動態(tài)指標

3.3.1 頻率偏移安全裕度概率分布

頻率偏移安全裕度與人為設定的頻率限值 fsm或時間限值tsm有關，可用圖3來定義和解釋。

圖3 從頻率限值和時間限值角度計算頻率偏移安全裕度Fig.3 Frequency deviation safety margin from the perspective of frequency limit and time limit

（1）從頻率限值的角度。根據曲線1和2定義的頻率偏移安全裕度

式中，ε 為從頻率限值角度計算的頻率偏移安全裕度標幺值，表征最低頻率fmin與頻率限值fsm的相差程度。當功率缺額不大時，頻率能較快恢復，此時最低頻率高于頻率限值，裕度ε＞0表示安全；當功率缺額很大，頻率下降至頻率限值以下，并且可能由于切負荷不足，使頻率無法恢復，這種惡劣情況發(fā)生的概率很小，此時裕度ε＜0表示不安全。

（2）從時間限值的角度。根據曲線3定義的頻率偏移安全裕度，可表征低頻運行時間tda與時間限值tsm的相差程度。

頻率偏移安全裕度為評價頻率穩(wěn)定提供了新的視角。需指出的是，裕度指標與選定的頻率限值和時間限值有很大關系，不同限值下的裕度指標相差較大，同一個系統(tǒng)可人為選定多對頻率限值和時間限值，以考察系統(tǒng)在不同限值上的安全裕度。

3.3.2 系統(tǒng)最低頻率和穩(wěn)態(tài)頻率概率分布

頻率動態(tài)變化過程中的最低頻率和穩(wěn)態(tài)頻率是頻率穩(wěn)定分析應重點關注的指標。各指標的概率分布反映了頻率穩(wěn)定可靠性水平的內在分布規(guī)律。

4 算例分析

對IEEE-RTS79[15]可靠性測試系統(tǒng)進行了評估，因文獻[15]未給出機組慣性時間常數(shù)、機組單位調節(jié)功率、負荷頻率調節(jié)效應系數(shù)等相關參數(shù)，本文參照文獻[8,11,15,16]作如下假設：

（1）系統(tǒng)中額定容量小于100MW的發(fā)電機單位調節(jié)功率標幺值取17，慣性時間常數(shù)取5.8s；大于等于100MW且小于200MW的發(fā)電機分別取20s和8.1s；大于200MW的發(fā)電機分別取25s和9.3s。

（2）負荷頻率調節(jié)效應取2.0；一階慣性環(huán)節(jié)時間常數(shù)取8s；頻率偏移安全裕度對應的頻率限值和時間限值分別取48.5Hz和10s；系統(tǒng)最低允許頻率取 47.0Hz；為簡化計算，僅對 RTS79系統(tǒng)中 1臺 350MW 發(fā)電機配置機組低頻跳閘保護，整定參數(shù)為[48.25Hz，20s]；系統(tǒng)低頻減載的整定見表1。

為分析旋轉備用和低頻切機對頻率穩(wěn)定的影響，規(guī)定：無備用時為情形 1，無備用且配置低頻切機為情形 2，有 10%旋轉備用為情形 3，有 10%旋轉備用且配置低頻切機為情形4。表 2為系統(tǒng)年峰荷2 850MW下，頻率失穩(wěn)風險指標的計算結果。

表3為低頻減載各輪期望切負荷量的計算結果。

表1 低頻減載整定表Tab.1 Setting of under frequency load shedding

表2 年峰荷下頻率失穩(wěn)指標計算結果Tab.2 Results of frequency instability for peak load

表3 年峰荷下各輪期望切負荷量指標計算結果Tab.3 Results of expected load shedding for peak load（單位：MW）

由表可見：

（1）旋轉備用能明顯改善頻率的穩(wěn)定性。無備用時，頻率失穩(wěn)概率為0.003 6，累計頻率失穩(wěn)次數(shù)為 2.99次/a，期望切負荷總量為 151.75MW；當投入10%旋轉備用時，失穩(wěn)概率下降為0.002 5，是無備用時的69.44%，累計失穩(wěn)次數(shù)為2.04次/a，是無備用時的68.23%；期望切負荷總量為113.86MW，是無備用時的75.03%。

（2）低頻切機加重了頻率失穩(wěn)概率和累計失穩(wěn)次數(shù)的嚴重程度，對比情形 4和情形3，嚴重程度分別上升了 36.0%和 38.5%，但期望切負荷總量僅增加 6.40%，并且前四輪和特殊輪切負荷量相等。這是因為：

（1）當功率缺額比較嚴重時，系統(tǒng)頻率雖然可以在前幾輪動作后開始恢復，但恢復速度有時較慢，可能在恢復至期望值之前，低頻運行時間已超出機組低頻跳閘整定值，機組被切除，頻率又開始下降，使最末幾輪啟動，換言之，低頻切機只影響最末幾輪的動作，而不影響前幾輪，且最末幾輪動作概率較小，所以使期望切負荷總量變化不大。

（2）由于整定值的原因，低頻切機必然動作于特殊輪之后，所以低頻切機不影響特殊輪的動作。

圖4為有無備用時頻率失穩(wěn)。由圖可見：

（1）旋轉備用的投入能延緩系統(tǒng)頻率失穩(wěn)。無備用時導致頻率失穩(wěn)的最小功率缺額為926.37MW，對應的失穩(wěn)時間為10.51s，有10%備用時頻率失穩(wěn)的最小功率缺額為 961.70MW，對應的失穩(wěn)時間為8.96s，此時無備用對應的失穩(wěn)時間為 4.82s，比有備用時快了4.14s。

（2）隨著功率缺額加大，兩者的失穩(wěn)時間逐漸接近。這是因為當功率缺額很嚴重，且旋轉備用還沒來得及發(fā)揮作用時，頻率的快速下降已經使頻率失穩(wěn)。

圖4 有無備用時頻率失穩(wěn)時間Fig.4 Frequency instability time under different spinning reserve situation

圖5為頻率偏移安全裕度的概率分布，由圖可見：

（1）對大部分故障狀態(tài)，頻率動態(tài)過程中的最低頻率高于選定的頻率限值，或低頻運行時間小于選定的時間限值，對應的頻率安全裕度為正，只有很少部分的嚴重故障使得安全裕度為負。

（2）備用的投入能改善頻率偏移安全裕度的概率分布。安全裕度在0.45以上時，有備用對應的概率分布較無備用時稍大，分別為0.558 7和0.488 5，安全裕度為負時，有無備用對應的概率分別為0.101 1和0.159 6。

圖6為有無備用時最低頻率概率分布。由圖可見，備用的投入使最低頻率分布較無備用時在高頻處概率較大，在低頻處概率較小，這對于改善系統(tǒng)頻率穩(wěn)定性有明顯效果。

圖5 頻率限值和時間限值下頻率偏移安全裕度概率分布Fig.5 Probability distribution of frequency deviation safety margin under frequency limit and time limit

圖7為有無備用時穩(wěn)態(tài)頻率概率分布。由圖可見，有備用時的穩(wěn)態(tài)頻率雖然在低頻處（＜49.4Hz）概率較低，但其超調和過控現(xiàn)象比較嚴重，穩(wěn)態(tài)頻率高于50.2Hz的概率為0.092 1，需引起特別關注。

圖7 有無備用時系統(tǒng)穩(wěn)態(tài)頻率概率分布Fig.7 Probability distribution of system steady state frequency under different spinning reserve situation

圖8給出了備用容量在5%～15%之間變化時系統(tǒng)指標相對無備用時的下降程度變化曲線。由圖可見，頻率穩(wěn)定性隨備用容量的增加得到不斷改善。當備用容量為15%時，頻率失穩(wěn)概率、累計失穩(wěn)次數(shù)和期望切負荷總量的嚴重程度分別下降47.22%、50.03%和33.20%。

圖8 系統(tǒng)指標隨旋轉備用變化曲線Fig.8 System indices with respect to spinning reserve variation

機組低頻保護的切機容量和切機時間可根據實際情況進行調整，以切機容量為例，假設切機容量可連續(xù)調整，圖9為低頻切機容量在350～550MW之間變化時指標相對無切機時的上升程度變化曲線。由圖可見，失穩(wěn)概率和累計失穩(wěn)次數(shù)隨切機容量增加而快速上升，而期望切負荷總量變化相對緩慢。

圖9 系統(tǒng)指標隨低頻切機容量變化曲線Fig.9 System indices with respect to generator tripping

上述算例使用系統(tǒng)年峰值負荷，未考慮負荷變化，得到的是年度化指標。年度指標則是考慮整年負荷變化計算得出的指標，由于年度指標能計及負荷隨機變化，在進行可靠性評估的成本—效益分析時更有說服力。下面在頻率穩(wěn)定可靠性評估中考慮年度負荷的變化，并計算了頻率穩(wěn)定年度指標。

利用均值聚類技術產生年負荷的多層水平模型，每一層負荷可按計算年度化指標時的恒定負荷水平求得該層負荷水平下的頻率穩(wěn)定指標，總年度指標可通過每層負荷下的指標加權該層負荷水平發(fā)生的概率求和獲得。如需計及負荷變化的不確定性，可使用正態(tài)分布來模擬。

針對RTS79系統(tǒng)選擇3種不同的負荷層數(shù)進行計算和比較，并采用七分段正態(tài)分布來考慮負荷的不確定性。有備用時頻率穩(wěn)定年度指標見表 4，其中N、Y分別表示不計及和計及負荷的不確定性。由表可見：

（1）各負荷模型所對應的指標接近，但都比年峰荷下的年度化指標小得多。

（2）計及負荷不確定性時，指標嚴重程度略有上升。

表4 不同層數(shù)負荷模型下系統(tǒng)指標的比較Tab.4 Comparison of system indices under different multiple load models

表5列出了計及負荷不確定性時，20層負荷模型中各級負荷對年度指標的貢獻程度，即各負荷水平下的指標占年度總指標的百分比。由表可見：

表5 20層負荷模型下各級負荷對指標的貢獻Tab.5 The contribution to the indices of every level load level for 20 grade load model

（1）占年峰荷百分數(shù)較高的前面幾個負荷級對總指標的貢獻較大。前6個負荷級（≥83.44%）對頻率失穩(wěn)概率指標和累計頻率失穩(wěn)次數(shù)指標的貢獻分別為90.22%和87.87%，前7個負荷級（≥80.66%）對期望切負荷總量指標的貢獻為86.93%。

（2）20個負荷級中，第4個負荷級（88.88%）對失穩(wěn)概率指標、累計失穩(wěn)次數(shù)指標和期望切負荷總量指標的貢獻都最大，分別為24.52%、24.63%和20.23%，這是因為該級負荷水平較高，且其發(fā)生的概率相對其他高級負荷也較大。

（3）當負荷水平較低（＜73.00%）時，對年度指標的貢獻很小，甚至可以忽略，因此對低級負荷可不必計算其指標，以節(jié)省計算量。

5 結論

本文提出的頻率失穩(wěn)概率評估考慮了旋轉備用、低頻減載和低頻切機對頻率動態(tài)特性的影響，同時考慮了發(fā)輸電元件和負荷水平的不確定性，基于所構建的頻率穩(wěn)定可靠性評估指標體系，對IEEE RTS79系統(tǒng)進行了驗算，并分別對旋轉備用容量和低頻切機容量變化對頻率穩(wěn)定性的影響、考慮負荷變化時的頻率穩(wěn)定年度指標進行了分析計算，驗證了該方法的有效性。本文所提方法能更加準確合理地表征系統(tǒng)頻率穩(wěn)定的可靠性水平，可為低頻減載的配置和機-網的協(xié)調問題提供概率角度的參考。

[1]Horne J, Flynn D, Littler T. Frequency stability issues for islanded power systems[C]. Proceedings of 2004 IEEE PES Power Systems Conference and Expositon,New York, USA, 2004, 299-306.

[2]IEEE/CIGRE joint task force on stability terms and definitions. Definition and classification of power system stability[J]. IEEE Transactions on Power System, 2004, 19(2): 1387-1401.

[3]孫華東, 湯涌, 馬世英. 電力系統(tǒng)穩(wěn)定的定義與分類述評[J]. 電網技術, 2006, 30(17): 31-35.Sun Huadong, Tang Yong, Ma Shiying. A commentary on definition and classification of power system stability[J]. Power System Technology, 2006,30(17): 31-35.

[4]蔡澤祥, 徐志勇, 申洪. 電力系統(tǒng)頻率穩(wěn)定分析的直接法[J]. 電力系統(tǒng)及其自動化學報, 1999, 11(5-6):13-17.Cai Zexiang, Xu Zhiyong, Shen Hong. A direct method for frequency stability analysis of power system[J]. Proceedings of the CSU-EPSA, 1999,11(5-6): 13-17.

[5]薛禹勝, 任先成, Wu Q H, 等. 關于低頻低壓切負荷決策優(yōu)化協(xié)調的評述[J]. 電力系統(tǒng)自動化, 2009,33(9): 100-107.Xue Yusheng, Ren Xiancheng, Wu Q H, et al. A review on optimization and coordination of under frequency/voltage load shedding[J]. Automation of Electric Power Systems, 2009, 33(9): 100-107.

[6]Pasand M S, Seyedi H. New centralized adaptive under frequency load shedding algorithms[C]. Large Engineering System Conference on Power Engineering, Montreal, Canada, 2007, 44-48.

[7]Chin V, Dong Z Y, Saha T K, et al. Adaptive and optimal under frequency load shedding[C].Australasian Universities Power Engineering Conference, Sydney, Australia, 2008, 1-6.

[8]DL 428—91. 電力系統(tǒng)自動低頻減負荷技術規(guī)定[S],1991.

[9]李光琦. 電力系統(tǒng)暫態(tài)分析[M]. 北京: 中國電力出版社, 2007.

[10]Li W Y. Risk assessment of power systems-models,methods, and applications[M]. USA and Canada:IEEE Press and Wiley&Sons, 2005.

[11]楊冠城. 電力系統(tǒng)自動裝置原理[M]. 北京: 水利電力出版社, 1995.

[12]Shrestha G B, Lee K C. Load shedding schedules considering probabilistic outages[C]. The 7th International Power Engeineering Conference,Singapore, 2005, 2: 950-955.

[13]Imai S, Yasuda T. UFLS program to ensure stable island operation[C]. Power Systems Conference and Exposition, New York, USA, 2004, 1-6.

[14]趙淵, 周家啟, 劉志宏. 大電網可靠性的序貫和非序貫蒙特卡洛仿真的收斂性分析及比較[J]. 電工技術學報, 2009, 24(11): 127-133.Zhao Yuan, Zhou Jiaqi, Liu Zhihong. Convergence analysis and comparison of sequential and nonsequential monte-carlo simulation for bulk power system reliability assessment[J]. Transactions of China Electrotechnical Society, 2009, 24(11):127-133.

[15]IEEE Task Force. IEEE reliability test system[J].IEEE Transaction on Power Apparatus and Systems,1979, 98(6): 2047-2054.

[16]陳珩. 電力系統(tǒng)穩(wěn)態(tài)分析[M]. 北京: 中國電力出版社, 2007.