李國(guó)慶，劉 玢，陳厚合

（東北電力大學(xué) 電氣工程學(xué)院，吉林 吉林 132012）

0 引言

由于風(fēng)電的隨機(jī)性、波動(dòng)性和間歇性［1-2］，風(fēng)電并網(wǎng)給電力系統(tǒng)區(qū)域間可用輸電能力ATC（Available Transfer Capability）［3-5］計(jì)算的可信性帶來(lái)巨大挑戰(zhàn)［6-8］。在ATC的計(jì)算中，容量效益裕度CBM（Capacity Benefit Margin）是指負(fù)荷所在處的主輸電系統(tǒng)為供電服務(wù)公司留取的一部分線路傳輸功率，以保證電力系統(tǒng)運(yùn)行的可靠性水平［9］。CBM留取的合理性和準(zhǔn)確性將直接影響所發(fā)布ATC值的可信度，進(jìn)而影響系統(tǒng)運(yùn)行的可靠性和經(jīng)濟(jì)性。隨著大規(guī)模風(fēng)電并網(wǎng)運(yùn)行，為了確保系統(tǒng)ATC評(píng)估的可信性和準(zhǔn)確性，需要對(duì)區(qū)域電網(wǎng)CBM的計(jì)算模型及多區(qū)域電網(wǎng)輸電斷面CBM的求取模型進(jìn)行更為深入的研究［10-11］。

當(dāng)前，電力系統(tǒng)對(duì)于CBM的研究尚不夠完善，關(guān)于新能源并網(wǎng)對(duì)CBM留取的影響、CBM留取對(duì)系統(tǒng)ATC可信度的影響及對(duì)系統(tǒng)運(yùn)行可靠性和經(jīng)濟(jì)性的影響尚未考慮在內(nèi)。文獻(xiàn)［12-13］采用傳統(tǒng)的確定性方法計(jì)算CBM，主要考慮傳統(tǒng)機(jī)組停機(jī)故障等不確定性因素，其發(fā)生概率較小，因此CBM大小一般取系統(tǒng)內(nèi)最大發(fā)電出力的一個(gè)倍數(shù)，或者取最大輸電能力TTC（Total Transfer Capability）的一個(gè)固定百分比。此類(lèi)計(jì)算方法并不能很好地反映系統(tǒng)變化對(duì)CBM的影響，特別是風(fēng)電并網(wǎng)后，該方法將難以應(yīng)對(duì)風(fēng)電波動(dòng)性的影響。文獻(xiàn)［14］按照與各區(qū)域缺電時(shí)間期望LOLE（Loss Of Load Expectation）指標(biāo)大小成反比的方法進(jìn)行輸電斷面CBM分配，此類(lèi)方法雖然步驟簡(jiǎn)單且計(jì)算速度快，但未考慮系統(tǒng)運(yùn)行的經(jīng)濟(jì)性，只是簡(jiǎn)單地按照LOLE指標(biāo)對(duì)CBM進(jìn)行分配。文獻(xiàn)［10］采用概率性方法對(duì)CBM進(jìn)行計(jì)算，并以各輸電斷面CBM之和最小為目標(biāo)對(duì)各輸電斷面CBM進(jìn)行經(jīng)濟(jì)分配，此類(lèi)方法可使系統(tǒng)的ATC值最大化，但并未考慮系統(tǒng)留取CBM的成本問(wèn)題和系統(tǒng)運(yùn)行經(jīng)濟(jì)性。

因此，本文首先針對(duì)大規(guī)模風(fēng)電并網(wǎng)對(duì)系統(tǒng)運(yùn)行可靠性的影響，引入可表征系統(tǒng)運(yùn)行可靠性的LOLE指標(biāo)［15］，通過(guò)計(jì)算系統(tǒng)發(fā)電裕度建立其與LOLE指標(biāo)的對(duì)應(yīng)關(guān)系；然后，綜合考慮風(fēng)電場(chǎng)出力隨機(jī)性、負(fù)荷波動(dòng)、常規(guī)機(jī)組故障等各種不確定性因素對(duì)系統(tǒng)發(fā)電裕度的影響，建立各不確定性因素與系統(tǒng)LOLE指標(biāo)的內(nèi)在聯(lián)系。對(duì)于區(qū)域電網(wǎng)，本文根據(jù)系統(tǒng)對(duì)可靠性指標(biāo)的不同要求，根據(jù)系統(tǒng)輸電斷面CBM大小與LOLE指標(biāo)的對(duì)應(yīng)關(guān)系，建立CBM計(jì)算模型；對(duì)于多區(qū)域電網(wǎng)，以區(qū)域電網(wǎng)CBM計(jì)算模型為基礎(chǔ)，構(gòu)建以經(jīng)濟(jì)性為目標(biāo)、滿足可靠性約束的送電區(qū)域各機(jī)組發(fā)電裕度分配及多區(qū)域電網(wǎng)各輸電斷面CBM留取的數(shù)學(xué)模型，力求在保證風(fēng)電并網(wǎng)后系統(tǒng)運(yùn)行可靠性的基礎(chǔ)上，提升系統(tǒng)運(yùn)行的經(jīng)濟(jì)性；最后，以IEEE 30節(jié)點(diǎn)和IEEE 118節(jié)點(diǎn)系統(tǒng)為例進(jìn)行分析驗(yàn)證。

1 區(qū)域電網(wǎng)CBM模型及計(jì)算

大規(guī)模風(fēng)電并網(wǎng)后，風(fēng)電出力的波動(dòng)性和隨機(jī)性增加了系統(tǒng)運(yùn)行的不確定性。根據(jù)CBM的物理意義，其為從互聯(lián)區(qū)域輸入電能提供了必要條件，當(dāng)系統(tǒng)運(yùn)行可靠性降低時(shí)，能夠通過(guò)輸電斷面預(yù)留的CBM從其互聯(lián)電網(wǎng)輸入電能，從而改善系統(tǒng)運(yùn)行的可靠性。因此，電力系統(tǒng)輸電斷面CBM的大小與系統(tǒng)要求的可靠性水平相關(guān)。本文基于系統(tǒng)對(duì)可靠性指標(biāo)的不同要求建立區(qū)域電網(wǎng)CBM計(jì)算模型以保障系統(tǒng)運(yùn)行的可靠性。

1.1 區(qū)域電網(wǎng)發(fā)電裕度建模

風(fēng)電所特有的電源特性成為影響系統(tǒng)運(yùn)行可靠性的關(guān)鍵因素。本文利用蒙特卡洛模擬法模擬風(fēng)電場(chǎng)的隨機(jī)因素，考慮風(fēng)電場(chǎng)風(fēng)速的隨機(jī)性，并根據(jù)風(fēng)力發(fā)電機(jī)組的輸出功率曲線進(jìn)行風(fēng)電機(jī)組出力預(yù)測(cè)，計(jì)算風(fēng)電場(chǎng)各時(shí)段綜合輸出功率。

根據(jù)系統(tǒng)發(fā)電量與負(fù)荷的關(guān)系，確定系統(tǒng)的發(fā)電裕度M，即系統(tǒng)的總發(fā)電量超出系統(tǒng)負(fù)荷的部分。借助M表征該區(qū)域的發(fā)電量是否足夠支撐其負(fù)荷，當(dāng)系統(tǒng)發(fā)電量不足時(shí)，通過(guò)輸電斷面上所留取的CBM，使得該區(qū)域可以從與其互聯(lián)的送電區(qū)域輸入電能，以滿足系統(tǒng)可靠性要求。

根據(jù)系統(tǒng)中常規(guī)機(jī)組出力模型（C）和風(fēng)電機(jī)組出力模型（W），求出系統(tǒng)的總出力G=W+C。由于G、W、C均為離散的隨機(jī)變量，故可利用卷積公式求取G的概率質(zhì)量函數(shù)：

其中，PG（W+C=z）為G數(shù)組的概率質(zhì)量函數(shù)；k為自變量；PW（W=z-k）為 W=z-k 時(shí)的概率；PC（C=k）為C=k時(shí)的概率。

考慮系統(tǒng)中負(fù)荷的波動(dòng)性，設(shè)負(fù)荷模型為L(zhǎng)，則系統(tǒng)的發(fā)電裕度為M=G-L。系統(tǒng)的總出力與總負(fù)荷是相互獨(dú)立的隨機(jī)變量，所以仍然利用卷積公式求取M的概率質(zhì)量函數(shù)：

其中，設(shè) G-L=z，PM（G-L=z）為 M 數(shù)組的概率質(zhì)量函數(shù)；k為自變量；PG（G=z+k）為 G=z+k時(shí)的概率；PL（L=k）為 L=k 時(shí)的概率。

求取M的概率質(zhì)量函數(shù)時(shí)，以式（2）為目標(biāo)，綜合考慮系統(tǒng)的有功平衡約束：

機(jī)組容量約束：

機(jī)組爬坡率約束：

其中，N為常規(guī)機(jī)組臺(tái)數(shù)；S為風(fēng)電機(jī)組臺(tái)數(shù)；L為負(fù)荷預(yù)測(cè)值表示風(fēng)電機(jī)組總出力值表示常規(guī)機(jī)組總出力值，pi（t）為第i臺(tái)常規(guī)機(jī)組的出力；pimin和pimax分別為機(jī)組i的最小和最大出力限值；RDamp，i和 RUamp，i分別為機(jī)組 i的有功出力下降速率和上升速率，單位為MW／min。

1.2 區(qū)域電網(wǎng)CBM計(jì)算

對(duì)某一確定的區(qū)域電網(wǎng)，當(dāng)風(fēng)電并網(wǎng)運(yùn)行后，根據(jù)上述模型，可求取風(fēng)電場(chǎng)出力模型W、常規(guī)機(jī)組出力模型C及負(fù)荷模型L。由式（1）、（2）以及約束條件式（3）—（5），可求得M的概率質(zhì)量函數(shù)并繪制M的概率分布函數(shù)曲線，如圖1中實(shí)曲線所示。LOLE指標(biāo)是指某個(gè)時(shí)段系統(tǒng)容量小于日最高負(fù)荷天數(shù)的期望值［16］。對(duì)M的負(fù)值區(qū)域的概率進(jìn)行累積分布處理，可求出此時(shí)M與LOLE的對(duì)應(yīng)關(guān)系。當(dāng)系統(tǒng)對(duì)LOLE指標(biāo)要求不同時(shí)，M的概率分布發(fā)生變化，進(jìn)而輸電斷面留取的CBM相應(yīng)變化。經(jīng)過(guò)計(jì)算，得到當(dāng)LOLE為α?xí)r，CBM為0。若系統(tǒng)對(duì)可靠性要求提高，LOLE為β且α>β時(shí)，M的概率分布發(fā)生變化，如圖1中虛曲線所示，進(jìn)而經(jīng)過(guò)計(jì)算得到輸電斷面留取φ MW CBM，當(dāng)系統(tǒng)為受電區(qū)域時(shí)便可以從其互聯(lián)的送電區(qū)域輸入不超過(guò)φ MW的電能，彌補(bǔ)系統(tǒng)發(fā)電量不足，供給負(fù)荷。

圖1 某時(shí)刻系統(tǒng)發(fā)電裕度的概率分布Fig.1 Probability distribution of generation margin at a particular moment

因此，根據(jù)不同系統(tǒng)對(duì)于LOLE指標(biāo)的不同要求，可求出該時(shí)段輸電斷面需要為本區(qū)域留取的CBM的大小。根據(jù)LOLE與CBM之間的對(duì)應(yīng)關(guān)系，可繪制其關(guān)系曲線，以CBM取值為橫坐標(biāo)，LOLE大小為縱坐標(biāo)，如圖2所示，以便在之后的實(shí)用計(jì)算中應(yīng)用該曲線根據(jù)系統(tǒng)對(duì)可靠性指標(biāo)的不同要求求取對(duì)應(yīng)的CBM的大小。若系統(tǒng)要求的可靠性指標(biāo)LOLE為θ h／a，則通過(guò)該曲線中所標(biāo)注的點(diǎn)可以求出，此時(shí)CBM為γ MW。即為滿足系統(tǒng)可靠性要求，該時(shí)段系統(tǒng)輸電斷面需要為其留取的CBM大小為γMW。

圖2 某時(shí)刻CBM-LOLE關(guān)系曲線Fig.2 CBM-LOLE curve at a particular moment

2 多區(qū)域電網(wǎng)CBM求取模型

當(dāng)互聯(lián)電網(wǎng)存在多個(gè)區(qū)域電網(wǎng)時(shí)（以區(qū)域A、B、C 3區(qū)電網(wǎng)為例），如果區(qū)域B、C同時(shí)為送電區(qū)域、區(qū)域A為受電區(qū)域，區(qū)域B、C至區(qū)域A的輸電斷面留取的CBM可分別記為指為滿足受電區(qū)域可靠性指標(biāo)，B或C送電區(qū)域單獨(dú)作用時(shí)輸電斷面為受電區(qū)域留取的CBM，由前文方法可分別計(jì)算得到B、C兩區(qū)域至A區(qū)輸電斷面留取的CBM之和，即該代數(shù)和已超過(guò)區(qū)域A的實(shí)際需求，造成一定程度的經(jīng)濟(jì)浪費(fèi)。故本文綜合考慮B、C送電區(qū)域協(xié)同向受電區(qū)域A輸送電能，提高受電區(qū)域可靠性并減少電能浪費(fèi)。此時(shí)系統(tǒng)輸電斷面中留取CBM的大小等于送電區(qū)域各機(jī)組為受電區(qū)域預(yù)留發(fā)電裕度的總和。因此，以送電區(qū)域各機(jī)組發(fā)電裕度獲取經(jīng)濟(jì)性最佳為目標(biāo)得到各時(shí)段送電區(qū)域各機(jī)組為受電區(qū)域預(yù)留發(fā)電裕度的大小，進(jìn)而求得各送電區(qū)域至受電區(qū)域輸電斷面上CBM的大小，使系統(tǒng)運(yùn)行可靠性得到保障且經(jīng)濟(jì)性提高。

2.1 目標(biāo)函數(shù)

以所有送電區(qū)域各機(jī)組為受電區(qū)域所留取發(fā)電裕度獲取總成本最小為多區(qū)域電網(wǎng)各輸電斷面CBM模型的目標(biāo)函數(shù)：

其中，n為可調(diào)度的機(jī)組總數(shù)；qi為機(jī)組i為受電區(qū)域預(yù)留的發(fā)電量，即通過(guò)輸電斷面上的CBM向受電區(qū)域傳輸功率的大??；Fi（qi）為機(jī)組i所預(yù)留發(fā)電裕度成本，以二次函數(shù)表示為其中，ai、bi、ci為給定常數(shù)。

2.2 約束條件

首先，各機(jī)組出力要在其可行范圍內(nèi)；其次，受電區(qū)域可靠性未達(dá)到系統(tǒng)要求時(shí)，送電區(qū)域應(yīng)在滿足其自身可靠性要求的基礎(chǔ)上，通過(guò)輸電斷面CBM將各機(jī)組留取的發(fā)電裕度輸送至受電區(qū)域，以彌補(bǔ)受電區(qū)域可靠性的不足，故約束條件如下：

其中，Pi為機(jī)組i的最大技術(shù)出力；pi為機(jī)組i現(xiàn)有的有功出力；CA為此時(shí)段受電區(qū)域A的發(fā)電不足量，為滿足系統(tǒng)可靠性要求所需從互聯(lián)區(qū)域輸入功率的大小，即輸電斷面需要為區(qū)域A留取CBM的大?。籦、c分別為此時(shí)段區(qū)域B、C機(jī)組中未滿發(fā)的機(jī)組數(shù)量；CB、CC為此時(shí)段送電區(qū)域B、C中各機(jī)組發(fā)電裕度大小的總和，即在滿足自身系統(tǒng)可靠性約束后可以通過(guò)輸電斷面中留取的CBM向區(qū)域A輸送電能的大小。

約束條件中，式（7）表示送電區(qū)域機(jī)組提供的通過(guò)輸電斷面中的CBM向受電區(qū)域傳輸?shù)墓β实拇笮≡谄淇尚蟹秶鷥?nèi)，即在保證其所在區(qū)域安排的發(fā)電量后可以提供的有功功率的范圍；式（8）表示送電區(qū)域向受電區(qū)域輸送一部分功率后能夠使受電區(qū)域滿足可靠性要求；式（9）、（10）表示送電區(qū)域通過(guò)輸電斷面中CBM為受電區(qū)域輸送功率后，其自身運(yùn)行仍然滿足可靠性水平要求。

3 模型的求解

遺傳算法GA（Genetic Algorithm）模擬物競(jìng)天擇的生物進(jìn)化過(guò)程，包括遺傳、突變、自然選擇以及雜交等。該方法已被引入廣泛的工程問(wèn)題中，進(jìn)而快速發(fā)展成一種“自適應(yīng)啟發(fā)式概率性迭代式全局搜索算法”。采用GA求解多區(qū)域電網(wǎng)送電區(qū)域各機(jī)組發(fā)電裕度經(jīng)濟(jì)分配模型流程如圖3所示。

圖3 GA流程圖Fig.3 Flowchart of GA

4 算例分析

4.1 IEEE 30節(jié)點(diǎn)系統(tǒng)算例分析

首先以IEEE 30節(jié)點(diǎn)3區(qū)系統(tǒng)為例進(jìn)行分析，該系統(tǒng)常規(guī)機(jī)組裝機(jī)容量為435 MW，區(qū)域A、B、C的節(jié)點(diǎn)16、3、28各加入10 MW的風(fēng)電場(chǎng)（考慮到風(fēng)電的波動(dòng)性，本算例選取風(fēng)電在電網(wǎng)中的貢獻(xiàn)率不超過(guò)10%），系統(tǒng)結(jié)構(gòu)如圖4所示。

采用蒙特卡洛抽樣24 h風(fēng)電場(chǎng)出力模型W、常規(guī)機(jī)組出力模型C及負(fù)荷模型L，取LOLE<1d／10 a（或LOLE<2.4 h／a），即10 a中負(fù)荷超過(guò)最大發(fā)電容量的時(shí)間累計(jì)最多不超過(guò)1d。由式（1）和（2）及CBMLOLE曲線分別計(jì)算各時(shí)段輸電斷面CBM值如圖5所示，CBM大于0代表系統(tǒng)滿足可靠性要求，各機(jī)組可通過(guò)輸電斷面上的CBM向外界提供電能；CBM小于0代表系統(tǒng)可靠性未滿足要求，需在輸電斷面為其預(yù)留CBM從互聯(lián)電網(wǎng)輸入電能。

圖4 IEEE 30節(jié)點(diǎn)系統(tǒng)圖Fig.4 IEEE 30-bus system

圖5 IEEE 30節(jié)點(diǎn)系統(tǒng)輸電斷面CBM留取情況Fig.5 Calculated CBM of transmission section of IEEE 30-bus system

圖6 區(qū)域A的發(fā)電量和負(fù)荷情況Fig.6 Power generation and load of region A

圖7 區(qū)域B的發(fā)電量和負(fù)荷情況Fig.7 Power generation and load of region B

圖8 區(qū)域C的發(fā)電量和負(fù)荷情況Fig.8 Power generation and load of region C

區(qū)域A、B、C的10 MW風(fēng)電場(chǎng)出力預(yù)測(cè)、常規(guī)機(jī)組出力及負(fù)荷預(yù)測(cè)數(shù)據(jù)分別如圖6、7、8所示。同樣采用本文所建區(qū)域電網(wǎng)CBM計(jì)算模型，根據(jù)區(qū)域A、B、C的出力及負(fù)荷情況計(jì)算各時(shí)段各區(qū)域輸電斷面上的CBM，計(jì)算結(jié)果如圖9所示。由圖9可見(jiàn)，區(qū)域A發(fā)電量不足，可靠性不滿足要求，需要輸電斷面為其預(yù)留CBM從互聯(lián)電網(wǎng)輸入電能；而區(qū)域B、C可靠性滿足要求，各機(jī)組可通過(guò)輸電斷面上的CBM向外界提供電能。限于篇幅，以第13個(gè)時(shí)段為例進(jìn)行分析：該時(shí)段，由于受電區(qū)域A系統(tǒng)可靠性未滿足要求，需從外界輸入37.2 MW功率，因此輸電斷面共需要為其留取CBM 37.2MW；送電區(qū)域B和C中各機(jī)組可分別通過(guò)輸電斷面向互聯(lián)區(qū)域輸送的功率總和分別為81.9789 MW和13.3401 MW。

圖9 IEEE 30節(jié)點(diǎn)系統(tǒng)各區(qū)域輸電斷面CBM留取情況Fig.9 Calculated CBM of three regional transmission sections of IEEE 30-bus system

采用本文送電區(qū)域各機(jī)組發(fā)電裕度經(jīng)濟(jì)分配模型，可將遺傳算法種群大小設(shè)為Nsize=100，迭代終止次數(shù)為200。根據(jù)各區(qū)域該時(shí)段發(fā)電情況及輸電斷面CBM值，計(jì)算可得到第13個(gè)時(shí)段系統(tǒng)送電區(qū)域各機(jī)組為受電區(qū)域預(yù)留發(fā)電裕度分配情況如表1所示。

各送電區(qū)域發(fā)電機(jī)組為受電區(qū)域預(yù)留發(fā)電裕度的總和即為該送電區(qū)域至受電區(qū)域輸電斷面留取CBM的大小，則各送電區(qū)域至受電區(qū)域輸電斷面CBM優(yōu)化結(jié)果見(jiàn)表 2，算法收斂過(guò)程見(jiàn)圖10。

表1 IEEE 30節(jié)點(diǎn)系統(tǒng)第13個(gè)時(shí)段送電區(qū)域機(jī)組發(fā)電裕度Table 1 Generation margin of units in feeding regions of IEEE 30-bus system in 13th period

表2 IEEE 30節(jié)點(diǎn)系統(tǒng)第13個(gè)時(shí)段輸電斷面留取的優(yōu)化CBMTable 2 Optimized CBM of transmission sections of IEEE 30-bus system in 13th period

圖10 IEEE 30節(jié)點(diǎn)系統(tǒng)第13個(gè)時(shí)段所需費(fèi)用最小化收斂過(guò)程Fig.10 Convergence process of cost minimization for IEEE 30-bus system in 13th period

若采用傳統(tǒng)的確定性CBM留取方法，即簡(jiǎn)單地將輸電斷面的CBM直接留取為一個(gè)確定的值，則本系統(tǒng)輸電斷面共需留取CBM為最大發(fā)電單元出力40 MW，且傳統(tǒng)方法直接按照與區(qū)域B、C LOLE指標(biāo)大小成反比的方法對(duì)CBM進(jìn)行分配。將本文概率性方法與傳統(tǒng)確定性方法進(jìn)行比較，結(jié)果如表3所示。由計(jì)算結(jié)果可見(jiàn)，傳統(tǒng)方法是將CBM直接取為一個(gè)較大值，來(lái)保證系統(tǒng)運(yùn)行可靠性，但無(wú)法保證經(jīng)濟(jì)性；而本文所構(gòu)建模型是在滿足系統(tǒng)要求的可靠性水平的基礎(chǔ)上將CBM取為一個(gè)較小的值，且經(jīng)濟(jì)性較好，更為合理。

表3 CBM計(jì)算和分配結(jié)果對(duì)比Table 3 Comparison of CBM calculation and allocation results

4.2 IEEE 118節(jié)點(diǎn)系統(tǒng)算例分析

IEEE 118節(jié)點(diǎn)系統(tǒng)可劃分為3個(gè)區(qū)域［17］。該系統(tǒng)總裝機(jī)容量為5790 MW。假設(shè)在系統(tǒng)中新增500 MW風(fēng)電，平均分配在A、B、C3個(gè)區(qū)域的節(jié)點(diǎn)39、67和102。采用區(qū)域電網(wǎng)輸電斷面留取CBM模型，可以求得當(dāng)LOLE為2.4 h／a時(shí)，該118節(jié)點(diǎn)系統(tǒng)各時(shí)段輸電斷面共為其留取CBM的大小，計(jì)算結(jié)果見(jiàn)圖11。

圖11 IEEE 118節(jié)點(diǎn)系統(tǒng)各時(shí)段輸電斷面CBM留取情況Fig.11 Calculated CBM of three regional transmission sections of IEEE 118-bus system

同樣，可求得A、B、C 3區(qū)域各時(shí)段輸電斷面留取CBM的情況。以第16個(gè)時(shí)段為例，此時(shí)區(qū)域A可靠性未滿足要求，為受電區(qū)域，需要輸電斷面留取201.1714 MW的CBM。區(qū)域B和C可靠性已滿足要求，為送電區(qū)域，各機(jī)組可分別通過(guò)輸電斷面上的CBM向外界輸送功率112 MW和534.8265 MW。由送電區(qū)域各機(jī)組發(fā)電裕度經(jīng)濟(jì)分配模型計(jì)算得，此時(shí)段該系統(tǒng)區(qū)域B、C各機(jī)組通過(guò)輸電斷面為區(qū)域A預(yù)留發(fā)電裕度經(jīng)濟(jì)分配結(jié)果如表4所示。進(jìn)而可分別求得由區(qū)域B、C至區(qū)域A輸電斷面上CBM的優(yōu)化結(jié)果如表5所示。圖12為經(jīng)濟(jì)優(yōu)化計(jì)算目標(biāo)函數(shù)迭代300次的收斂過(guò)程。

若采用傳統(tǒng)方法對(duì)IEEE 118節(jié)點(diǎn)3區(qū)系統(tǒng)進(jìn)行CBM計(jì)算，可以得到此時(shí)輸電斷面為區(qū)域A共留取CBM為450 MW，明顯多于本文方法所求得的201.1714 MW，故其所需經(jīng)濟(jì)成本必然高于本文方法的6.2088×105$，對(duì)比結(jié)果如表6所示。

表4 系統(tǒng)第16個(gè)時(shí)段送電區(qū)域發(fā)電裕度經(jīng)濟(jì)分配結(jié)果Table 4 Results of generation margin economic allocation for feeding regions of system in 16th period

表5 IEEE 118節(jié)點(diǎn)系統(tǒng)第16個(gè)時(shí)段各輸電斷面CBM的優(yōu)化結(jié)果Table 5 Optimized CBM of transmission sections of IEEE 118-bus system in 16th period

圖12 系統(tǒng)第16個(gè)時(shí)段費(fèi)用最小化收斂過(guò)程Fig.12 Convergence process of cost minimization for system in 16th period

表6 CBM計(jì)算和分配結(jié)果對(duì)比Table 6 Comparison of CBM calculation and allocation results

由分析可見(jiàn)，本文所構(gòu)建模型既保證了各時(shí)段系統(tǒng)運(yùn)行的可靠性水平，又避免了資源浪費(fèi)，實(shí)現(xiàn)了系統(tǒng)運(yùn)行的經(jīng)濟(jì)性。因此，在電力系統(tǒng)區(qū)域間ATC的評(píng)估與決策中，采用基于可靠性指標(biāo)并兼顧系統(tǒng)運(yùn)行經(jīng)濟(jì)性的CBM模型是合理和有效的。

5 結(jié)論

本文提出了基于可靠性指標(biāo)的CBM概率性計(jì)算模型，并對(duì)互聯(lián)系統(tǒng)提出了以送電區(qū)域各機(jī)組預(yù)留發(fā)電裕度獲取經(jīng)濟(jì)性最佳為目標(biāo)的CBM分配模型。本文方法具有以下特點(diǎn)。

a.不再按照傳統(tǒng)的確定性方法將輸電斷面留取的CBM取為一個(gè)固定值，而是采用概率性方法進(jìn)行CBM計(jì)算，根據(jù)系統(tǒng)發(fā)電裕度M，并引入LOLE指標(biāo)，不同系統(tǒng)對(duì)于可靠性水平的要求不同時(shí)，其取值不同，通常系統(tǒng)所要求的LOLE取值為2.4 h／a。充分考慮風(fēng)電的波動(dòng)性，通過(guò)繪制CBM-LOLE曲線并結(jié)合系統(tǒng)的可靠性要求，求取CBM的值。采用該方法計(jì)算后，輸電斷面留取CBM的大小有所降低，但仍然滿足系統(tǒng)運(yùn)行可靠性要求，使得CBM的留取更加科學(xué)、可靠，CBM的合適取值減少了電力系統(tǒng)失負(fù)荷的概率，使發(fā)布的ATC的值更加可靠，為系統(tǒng)創(chuàng)造了效益。

b.對(duì)互聯(lián)系統(tǒng)輸電斷面進(jìn)行CBM計(jì)算，不再采用未考慮系統(tǒng)運(yùn)行經(jīng)濟(jì)性的傳統(tǒng)方法只按照與各送電區(qū)域LOLE指標(biāo)大小成反比的方法進(jìn)行經(jīng)濟(jì)分配，而是以送電區(qū)域各機(jī)組預(yù)留發(fā)電裕度獲取經(jīng)濟(jì)性最優(yōu)為目標(biāo)，滿足可靠性為約束，得到送電區(qū)域各機(jī)組預(yù)留發(fā)電裕度的分配方案，進(jìn)而得到各輸電斷面CBM的值，確定最佳的運(yùn)行方式，使系統(tǒng)運(yùn)行的可靠性與經(jīng)濟(jì)性相結(jié)合。

參考文獻(xiàn)：

［1］WANG C，LU Z，QIAO Y.A consideration of the wind power benefits in day-ahead scheduling of wind-coal intensive power systems［J］.IEEE Transactions on Power Systems，2013，28（1）：236-245.

［2］DING Y，WANG P，GOEL L，et al.Long-term reserve expansion of power systems with high wind power penetration using universal generating function methods［J］.IEEE Transactions on Power Systems，2011，26（2）：766-774.

［3］NERC.Available transfer capability definitions and determination［R］.New Jersey，USA：North American Electricity Reliability Council，1996.

［4］陳厚合，李國(guó)慶，姜濤．計(jì)及靜態(tài)電壓穩(wěn)定約束的交直流系統(tǒng)可用輸電能力［J］．電網(wǎng)技術(shù)，2012，36（2）：75-81．CHEN Houhe，LI Guoqing，JIANG Tao.Available transfer capability of hybrid AC／DC power system considering constraint of static voltage stability［J］.Power System Technology，2012，36（2）：75-81.

［5］SALIM N A，OTHMAN M M，SERWAN M S，et al.Determination of available transfer capability with implication of cascading collapse uncertainty［J］.IET Generation Transmission&Distribution，2014，8（4）：705-715.

［6］LUO G，CHEN J，CAI D，et al.Probabilistic assessment of available transfer capability considering spatial correlation in wind power integrated system［J］.IET Generation Transmission&Distribution，2013，7（12）：1527-1535.

［7］SHAYESTEH E，HOBBS B F，SODER L，et al.ATC-based system reduction for planning power systems with correlated wind and loads［J］.IEEE Transactions on Power Systems，2015，30（1）：429-438.

［8］ZHOU Z，BOTTERUD A.Dynamic scheduling of operating reserve in co-optimized electricity markets with wind power［J］.IEEE Transactions on Power Systems，2014，29（1）：160-171.

［9］NERC.Transmission capability margins and their use in ATC determination［R］.New Jersey，USA：North American Electricity Reliability Council，1999.

［10］RAMEZANI M，HAGHIFAM M，SINGH C，et al.Determination of capacity benefit margin in multiarea power systems using particle swarm optimization［J］.IEEE Transactions on Power Systems，2009，24（2）：631-641.

［11］AHMADI-KHATIR A，CONEJO A J，CHERKAOUUI R.Multiarea unit scheduling and reserve allocation under wind power uncertainty［J］.IEEE Transactions on Power Systems，2014，29（4）：1701-1710.

［12］EJEBE G C，WAIGHT J G，SANOTS-NIETO M，et al.Fast calculation of linear available transfer capability［J］.IEEE Transactions on Power Systems，1999，15（3）：255-260.

［13］LYON J D，HEDMAN K W，ZHANG M.Reserve requirements

（continued on page 14）（continued from page 6）to efficiently manage intra-zonal congestion［J］.IEEE Transactions on Power Systems，2014，29（1）：251-258.

［14］OU Y，SINGH C.Assessment of available transfer capability and margins［J］.IEEE Power Engineering Review，2002，22（5）：69.

［15］BENIDRIS M，MITRA J.Reliability and sensitivity analysis of composite power systems under emission constraints［J］.IEEE Transactions on Power Systems，2014，29（1）：404-412.

［16］MATOS M A，BESSA R J.Setting the operating reserve using probabilistic wind power forecasts［J］.IEEE Transactions on Power Systems，2011，26（2）：594-603.

［17］MEHRJERDI H，LEFEBVRE S，SAAD M，et al.A decentralized control of partitioned power networks for voltage regulation and prevention against disturbance propagation［J］.IEEE Transactions on Power Systems，2013，28（2）：1461-1469.