吳擁勛，江岳文，溫步瀛

（福州大學(xué) 電氣工程與自動化學(xué)院，福建 福州 350108）

為了應(yīng)對全球氣候變化，緩解資源環(huán)境約束，被稱為綠色清潔能源的風(fēng)能越來越受到世界各國的重視。風(fēng)電良好的節(jié)能減排效益決定了對含有風(fēng)電場的系統(tǒng)進行隨機生產(chǎn)模擬、計算其可靠性指標(biāo)和動態(tài)費用的必要性。

隨機生產(chǎn)模擬是一種考慮機組隨機故障和系統(tǒng)負(fù)荷隨機性，從而計算出最優(yōu)運行方式下各發(fā)電機組在一段時間內(nèi)的發(fā)電量以及系統(tǒng)可靠性指標(biāo)的算法[1]。目前，國內(nèi)外提出的隨機生產(chǎn)模擬主要有：等效電量函數(shù)法[2-4]、半不變量法（累積量法）[5]、分段直線逼近法[6]、分塊法[7]。這些算法的核心都是將時序負(fù)荷曲線轉(zhuǎn)化為持續(xù)負(fù)荷曲線，通過增大持續(xù)負(fù)荷來體現(xiàn)機組的隨機停運。文獻[8-9]提出了一種基于時序負(fù)荷曲線的電力系統(tǒng)概論性生產(chǎn)模擬方法。以上文獻的方法適用于不含風(fēng)電機組的系統(tǒng)。由于風(fēng)電的反調(diào)峰特性，給電力系統(tǒng)調(diào)頻、調(diào)峰帶來了難度[10-12]，因此，在對含有風(fēng)電機組的系統(tǒng)使用隨機生產(chǎn)模擬進行模擬時，就需要考慮風(fēng)電場的影響[13-14]。

文獻[2-7]將時序負(fù)荷曲線轉(zhuǎn)化為持續(xù)負(fù)荷曲線的同時丟失了與時序相關(guān)的信息和約束，在全面分析生產(chǎn)成本方面已難以滿足電力系統(tǒng)的實際需求，因此文獻[15]提出等效電量頻率法。該方法在一定程度計及了負(fù)荷的時間信息，通過修正后的等效負(fù)荷頻率曲線來計算機組期望開機次數(shù)，既可以得到傳統(tǒng)算法計算出的機組生產(chǎn)電量和可靠性指標(biāo)，還可以計算與機組啟停、維護相關(guān)的動態(tài)費用。文獻[15]建立考慮計及風(fēng)電機組的隨機生產(chǎn)模擬改進算法，但在計算機組啟停次數(shù)時可以繼續(xù)深入研究。本文針對文獻[15]中的等效電量頻率法機組啟停次數(shù)計算方法進行改進，將除第i臺機組外的其他機組故障對第i臺機組的影響納入到機組啟停次數(shù)計算中，使等效電量頻率法中計算的機組啟停次數(shù)考慮得更周全，并在算例中以改進后的等效電量頻率法與原來的等效電量頻率法計算結(jié)果進行比較。同時，由于本文考慮了應(yīng)變分量f（i）3（x），優(yōu)化了機組組合和啟停次數(shù)期望值，系統(tǒng)的失負(fù)荷概率和缺電成本都將得到改善。

1 等效電量頻率法的基本原理

1.1 等效電量函數(shù)法簡述

等效持續(xù)負(fù)荷曲線（ELDC）是隨機生產(chǎn)模擬技術(shù)的重要概念。如圖1所示，t=F（x）為一條持續(xù)負(fù)荷曲線，其橫坐標(biāo)表示系統(tǒng)的負(fù)荷，縱坐標(biāo)表示負(fù)荷的持續(xù)時間。圖中，T為研究周期；PLMAX為系統(tǒng)最大負(fù)荷；曲線上任一點（xi，ti）表示負(fù)荷大于或等于xi的持續(xù)時間為ti；假設(shè)系統(tǒng)總裝機容量為CS，Ci為第i臺機組的額定容量。從圖1可知，在安排完所有機組后，陰影部分即為電力不足期望值（EEENS）。

圖1 持續(xù)負(fù)荷曲線Fig. 1 Durative load curve

由文獻[1]可知，電量函數(shù)：

安排第i臺機組后，考慮機組隨機停運的電量函數(shù)：

第i臺發(fā)電機組所生產(chǎn)的電量：

假設(shè)系統(tǒng)中一共有N臺發(fā)電機組，當(dāng)N臺發(fā)電機全部安排完畢，電量函數(shù)為E（N）（J）。這時系統(tǒng)電量不足期望值：

失負(fù)荷概率：

1.2 負(fù)荷頻率曲線的含義

負(fù)荷轉(zhuǎn)移頻率代表的是在研究周期T內(nèi)，負(fù)荷水平x向上方向穿過時序負(fù)荷曲線的平均次數(shù)σ與研究周期T的比值，則負(fù)荷轉(zhuǎn)移頻率為fL（x）=σ/T。

負(fù)荷頻率曲線形象反映了研究周期T內(nèi)負(fù)荷水平到達(dá)x的次數(shù)。由于安排發(fā)電機組開機投入運行和停機都是根據(jù)負(fù)荷的變化來進行的，每當(dāng)負(fù)荷達(dá)到新的負(fù)荷水平時，就會安排相應(yīng)機組承擔(dān)增加的負(fù)荷。反之，負(fù)荷下降時，相應(yīng)的機組就會停機備用。因此，負(fù)荷為x處的機組在研究周期T內(nèi)的啟停機次數(shù)與fL（x）是對應(yīng)的。

1.3 負(fù)荷頻率曲線的建立

本文采用文獻[16]中形成負(fù)荷頻率表的方法，求取負(fù)荷轉(zhuǎn)移頻率fL（x）。

設(shè)負(fù)荷xi出現(xiàn)是時間為tij（j為負(fù)荷出現(xiàn)變化的時間），因此，點（j，xi）的增量頻率為fij=1/T。負(fù)荷xi的總增量頻率為：

負(fù)荷xi的累積頻率為：

由累積頻率可得負(fù)荷轉(zhuǎn)移頻率：

1.4 修正等效負(fù)荷頻率曲線

設(shè)1.3節(jié)求得的負(fù)荷頻率曲線為f（0）L（x），安排完第i－1臺發(fā)電機組后等效負(fù)荷頻率曲線為f（i-1）L（x），則安排完第i臺發(fā)電機組的等效負(fù)荷頻率曲線為：

循環(huán)分量f（i）2（x）反映了第i臺機組按照計劃開機、運行的過程中，隨時可能因為自身故障而引起自身計劃外的啟停[16]。因此，循環(huán)分量：

應(yīng)變分量f（i）3（x）分析的內(nèi)容與循環(huán)分量f（i）2（x）一樣，反映的都是機組在運行過程中，隨機可能發(fā)生停運而引起計劃外的啟停。然而，循環(huán)分量是關(guān)于機組自身故障引起的啟停次數(shù)變化的分析，與之不同的是，應(yīng)變分量反映的是開機順序中由于上一臺機組乃至上n臺機組的故障停運而引起的下一臺機組的啟停次數(shù)增加。例如，第i－1臺機組在計劃的運行過程中發(fā)生故障，在不額外安排備用機組的情況下，其出力由第i臺機組承擔(dān)，因此第i臺機組的開機次數(shù)加1（當(dāng)?shù)趇臺機組的出力不足以替代第i－1臺機組的作用時，第i+1臺機組，…，第i+N臺機組相繼開機替代，這時相應(yīng)的機組開機次數(shù)都要加1）；當(dāng)故障修復(fù)后，第i－1臺機組重新啟動運行，此時，負(fù)荷不變或者負(fù)荷變小的情況下，第i臺機組停機次數(shù)加1，反之，第i臺機組繼續(xù)運行。由于考慮了除本臺機組外的其他機組故障影響，因此，應(yīng)變分量：

其中：

τi是平均循環(huán)周期，亦稱作機組平均故障間隔時間。設(shè)Mi－l+1作為第i臺的狀態(tài)決策值。設(shè)第i－l+1臺機組發(fā)生故障，Ci－l+1為第i－l+1臺機組額定容量，當(dāng)Ci－l+1≤Ci－l+2，Mi－l+1=0，表明第i－l+1臺機組因故障而停止運行所需要代替的那部分發(fā)電量由第i－l+2臺機組替代發(fā)電，而不需由第i臺機組開機來代替發(fā)電；當(dāng)Ci－l+2

將式（10）—（12）進行離散化，可得

根據(jù)等效電量函數(shù)的定義，可得

第i臺機組的期望開機次數(shù)為：

本文將風(fēng)電場時序出力曲線從時序負(fù)荷曲線中分離得到凈負(fù)荷曲線，通過凈負(fù)荷曲線計算負(fù)荷頻率曲線，而負(fù)荷頻率曲線是關(guān)于負(fù)荷上升或者下降的次數(shù)，并且發(fā)現(xiàn)機組啟停次數(shù)與負(fù)荷頻率有一定關(guān)聯(lián)性，因此可以使用負(fù)荷頻率曲線計算機組啟停次數(shù)。 由于承擔(dān)負(fù)荷時盡量先使用風(fēng)電，因此風(fēng)電機組始終處于運行狀態(tài)，并且風(fēng)電出力的不確定性，因此在文中沒有將風(fēng)電機組納入等小負(fù)荷頻率曲線的修正，因此也就沒有納入機組啟停機次數(shù)計算中。

設(shè)系統(tǒng)中共有N臺機組，在文獻[15]中采用新指標(biāo)機組單位容量啟停頻率IFGSUC來綜合評價負(fù)荷波動和風(fēng)電場對系統(tǒng)常規(guī)機組啟停造成的影響。

2 動態(tài)費用分析

2.1 常規(guī)機組動態(tài)費用計算

常規(guī)機組（不含風(fēng)電機組）的動態(tài)費用Ct可按下式計算（對于含核電機組的系統(tǒng)，為保證其安全及經(jīng)濟性，在較短的周期內(nèi)核電機組一般不啟停，將其啟停費用置0處理[16]）：

式中，Cfuel為燃料費用；Co&m為運行和維護費用；Cuec為缺電成本；Cenvi為環(huán)境成本；Ctsu為汽輪機啟動費用；Cbsu為鍋爐啟動費用；Csd停機費用。

系統(tǒng)中N臺汽輪機組的啟動費用為：

式中，Ctsu，i為第i臺機組的單位發(fā)電量每次的汽輪機啟動費用。

計算Cbsu前要先計算機組的平均離線時間：

鍋爐啟動總費用為：

式中，Cbsu，i為第i臺機組的單位發(fā)電量每次的鍋爐啟動費用。

另外，可以根據(jù)實際負(fù)荷變化情況來選擇鍋爐應(yīng)該處于熱備用還是冷備用[17]。其他費用計算同文獻[15]。

2.2 風(fēng)電機組動態(tài)費用計算

計及含風(fēng)電機組的系統(tǒng)的生產(chǎn)成本Cz:

式中，Cad為風(fēng)電的可避免費用，主要是由于風(fēng)電接入減少了火電機組所造成的燃料費用、運行費用和環(huán)境成本。

在增加應(yīng)變分量f（i）3（x）后，各機組的啟停次數(shù)都會有相應(yīng)的改變。由于考慮了當(dāng)?shù)趇－l+1臺機組故障后的出力將由發(fā)電計劃中開機順序在第i－l+1之后的機組承擔(dān)，因此當(dāng)有機組故障時，負(fù)責(zé)替代出力的機組的開機和停機次數(shù)必然增加，這將在一定程度上降低系統(tǒng)的缺電期望值，而每臺機組的啟停費用也將會有所變化。

3 算例分析

本文所采用算例與文獻[15]相同，系統(tǒng)為EPRI 36機組，系統(tǒng)裝機容量為8 800 MW。負(fù)荷曲線采用IEEE RTS負(fù)荷數(shù)據(jù)[19]，年最大負(fù)荷為8 502 MW，模擬時間為1月份前30 d，共720 h。機組可靠性數(shù)據(jù)和運行維護費用見文獻[18]，機組啟停參數(shù)和費用見文獻[19]，環(huán)境成本數(shù)據(jù)見文獻[20]。風(fēng)電裝機共有1 200 MW，占據(jù)系統(tǒng)總裝機容量的12%。算例主要是對加入風(fēng)電機組后的系統(tǒng)進行隨機生產(chǎn)模擬，進行成本計算后，在與文獻[15]中的計算結(jié)果進行比較。模擬結(jié)果見表1。

表1 隨機生產(chǎn)模擬結(jié)果Tab. 1 Probabilistic production simulation result

從費用上看，缺電成本下降了7.6%，啟停費用增加了0.53%，燃料費用下降了0.08%，系統(tǒng)的機組單位容量啟停頻率指標(biāo)IFGSUC略微增大，本文模擬的總費用下降0.04%。從失負(fù)荷概率（LOLP）上看，改進等效電量頻率法后，本文的失負(fù)荷概率（LOLP）下降了0.35%。雖然機組的啟停費用增加，但是系統(tǒng)的電量不足期望值（EEENS）和失負(fù)荷概率（LOLP）有所下降，并且系統(tǒng)下降的缺電成本大于系統(tǒng)增加的啟停費用，因此改進后的算法優(yōu)化了系統(tǒng)機組組合，提高了發(fā)電系統(tǒng)的可靠性和經(jīng)濟性。

4 結(jié)論

本文充分挖掘了等效電量頻率法體現(xiàn)負(fù)荷時序特性的潛力，最后通過EPRI 36機組系統(tǒng)的隨機生產(chǎn)模擬結(jié)果可以得出以下的結(jié)論：

1）本文改進了等效電量頻率法，提出了應(yīng)變分量的概念，考慮發(fā)電計劃中開機順序在第i臺機組前的機組因故障造成后續(xù)機組啟停次數(shù)改變的影響。

2）缺電期望值、缺電頻率相較于原始的等效電量頻率法有所下降，說明改進后的等效電量頻率法相對提高了系統(tǒng)的可靠性。

3）改進后的等效電量頻率法雖然增加了機組的啟停的次數(shù)和啟停費用，但是同時也降低了系統(tǒng)的缺電期望值，最后總費用還有一定減少，說明改進后的等效電量頻率法相對提高了系統(tǒng)的經(jīng)濟性。

[1] 王錫凡. 電力系統(tǒng)隨機生產(chǎn)模擬的等效電量函數(shù)法[J].西安交通大學(xué)學(xué)報，1984，18（6）：13-26.WANG Xifan. EEF approach to power system probabilistic modeling[J]. Journal of Xi’an Jiaotong University，1984，18（6）: 13-26（in Chinese）.

[2] WANG Xifan.Equivalent energy function approach to power system probabilistic modeling[J]. IEEE Trans on Power Systems，1988，3（3）: 823-829.

[3] 周景宏，胡兆光，田建偉，等. 含能效電廠的電力系統(tǒng)生產(chǎn)模擬[J]. 電力系統(tǒng)自動化，2010，34（18）: 27-31.ZHOU Jinghong，HU Zhaoguang，TIAN Jianwei，et al.Power system production simulation including efficiency power plant[J]. Automation of Electric Power Systems，2010，34（18）: 27-31（in Chinese）.

[4] 李林川，王錫凡，王秀麗. 基于等效電量函數(shù)法的互聯(lián)電力系統(tǒng)隨機生產(chǎn)模擬[J]. 中國電機工程學(xué)報，1996，16（3）: 180-184.LI Linchuan，WANG Xifan，WANG Xiuli. Probabilistic modeling for interconnected power systems based on the equivalent energy function approach[J]. Proceedings of the CSEE，1996，16（3）: 180-184（in Chinese）.

[5] 陳剛，相年德，陳雪青. 一種基于負(fù)荷分解的隨機生產(chǎn)模擬新方法[J]. 中國電機工程學(xué)報，1992，12（3）: 47-52.CHEN Gang，XIANG Niande，CHEN Xueqing. A new method for probabilistic production simulation based on load decomposition[J]. Proceedings of the CSEE，1992，12（3）: 47-52（in Chinese）.

[6] BOOTH R R. Power system simulation model based on probability analysis[J]. IEEE Transactions on Power Systems，1972，91（1）：62-69.

[7] SCHENK K F，MISRA R B，VASSOS S，et al. A new method for the evaluating of expected energy generation and loss of probability[J]. IEEE Transactions on Power Systems，1984，103（2）：294-303.

[8] 夏清，王少軍，相年德. 時序負(fù)荷曲線下電力系統(tǒng)概率性生產(chǎn)模擬[J]. 中國電機工程學(xué)報，1994，14（3）：21-28.XIA Qing，WANG Shaojun，XIANG Niande. The probabilistic power system production simulation based on chronological load curve[J]. Proceedings of the CSEE，1994，14（3）：21-28（in Chinese）.

[9] 夏清，王少軍，相年德，等. 時序負(fù)荷曲線下電力系統(tǒng)概率性水電生產(chǎn)模擬[J]. 中國電機工程學(xué)報，1998，18（6）：429-433.XIA Qing，WANG Shaojun，XIANG Niande，et al. Hydro station’s probabilistic production simulation[J]. Proceedings of the CSEE，1998，18（6）: 429-433（in Chinese）.

[10] 洪翠，溫步瀛，陳群，等. 風(fēng)電場出力特性及其不確定性的對策分析[J]. 電網(wǎng)與清潔能源，2012，28（11）: 65-68.HONG Cui，WEN Buying，CHEN Qun，et al. Analysis of wind farm output characteristics and solutions to its uncertainty[J]. Power System and Clean Energy，2012，28（11）:65-68（in Chinese）.

[11] 姚金雄，張世強. 基于調(diào)峰能力分析的電網(wǎng)風(fēng)電接納能力研究[J]. 電網(wǎng)與清潔能源，2010，26（7）: 25-28.YAO Jinxiong，ZHANG Shiqiang. Analysis on capacity of wind power integration into grid based on peak load regulation[J]. Power System and Clean Energy，2010，26（7）: 25-28（in Chinese）.

[12] 魏磊，張琳，姜寧，等.包含風(fēng)電的電力系統(tǒng)調(diào)峰能力計算方法探討[J]. 電網(wǎng)與清潔能源，2010，26（8）: 59-63.WEI Lei，ZHANG Lin，JIANG Ning，et al. Discussion on computing methods for peak load regulation ability of power system connected with large-scale wind farms[J].Power System and Clean Energy，2010，26（8）: 59-63（in Chinese）.

[13] 范李平，楊力森，武粉桃. 風(fēng)電場并網(wǎng)對電力系統(tǒng)穩(wěn)定性影響[J]. 電網(wǎng)與清潔能源，2009，25（6）: 58-61.FAN Liping，YANG Lisen，WU Fentao. Influence of wind farm interconnected to power grid on stability of power system[J]. Power System and Clean Energy，2009，25（6）:58-61（in Chinese）.

[14] 范高鋒，趙海翔，戴慧珠. 大規(guī)模風(fēng)電對電力系統(tǒng)的影響和應(yīng)對策略[J]. 電網(wǎng)與清潔能源，2008，24（7）: 44-48.FAN Gaofeng，ZHAO Haixiang，DAI Huizhu. The impact and countermeasure of large scale wind power on power system[J]. Power System and Clean Energy，2008，24（7）:44-48（in Chinese）.

[15] 張節(jié)譚，程浩忠，胡澤春，等. 含風(fēng)電場的電力系統(tǒng)隨機生產(chǎn)模擬[J]. 中國電機工程學(xué)報，2009，29（28）：34-39.ZHANG Jietan，CHENG Haozhong，HU Zechun，et al.Power system probabilistic production simulation including wind farms[J]. Proceedings of the CSEE，2009，29（28）:34-39（in Chinese）.

[16] 曲翀，王秀麗，謝紹宇，等. 含風(fēng)電電力系統(tǒng)隨機生產(chǎn)模擬的改進算法[J]. 西安交通大學(xué)學(xué)報，2012，46（6）：116-117.QU Chong，WANG Xiuli，XIE Shaoyu，et al. A improved algorithm for probabilistic production simulation of power system with wind power[J]. Journal of Xi’an Jiaotong University，2012，46（6）: 116-117（in Chinese）.

[17] MALIK A S. An application of frequency and duration approach in generation planning[J]. IEEE Trans on Power Systems，1997，12（3）: 1076-1084.

[18] The Reliability Test System Task Force of the Application of Probability Methods Subcommittee. IEEE reliability test system[J]. IEEE Transaction on Power Apparatus and Systems，1979，98（6）: 2047-2054.

[19] MALIK A S，CORY B J. An application of FD approach in generation planning[J]. IEEE Transactions on Power Systems，1989，4（2）: 419-425.

[20] 陸華，周浩. 發(fā)電廠的環(huán)境成本分析[J]. 環(huán)境保護，2004（4）: 51-54.LU Hua，ZHOU Hao. Environmental cost analysis of power plants[J]. Environmental Protection，2004（4）: 51-54（in Chinese）.