劉興龍，張 娜，徐 江，姜百超，劉莉莉

(1.中國(guó)能源建設(shè)集團(tuán)遼寧電力勘測(cè)設(shè)計(jì)院有限公司，遼寧 沈陽(yáng) 110179；2.國(guó)網(wǎng)遼寧省電力有限公司經(jīng)濟(jì)技術(shù)研究院，遼寧 沈陽(yáng) 110015)

深入理解節(jié)能發(fā)電調(diào)度[1-3]中的每一個(gè)環(huán)節(jié)至關(guān)重要，而發(fā)電組合方案的制定是其中的重中之重。在發(fā)電組合方案的制定過(guò)程中，只考慮發(fā)電機(jī)組的供電煤耗是非常片面的,不能真正達(dá)到能耗最低的目標(biāo)，而必須要對(duì)發(fā)電機(jī)組的供電煤耗率進(jìn)行網(wǎng)損修正。網(wǎng)損修正過(guò)程中引起的微小變化都可能會(huì)影響發(fā)電機(jī)組上網(wǎng)發(fā)電與否，關(guān)系到發(fā)電企業(yè)的重大利益。因此，網(wǎng)損修正系數(shù)的計(jì)算關(guān)系重大，必須要準(zhǔn)確可信，使各方滿(mǎn)意。而網(wǎng)損修正系數(shù)的精確計(jì)算非常復(fù)雜，難以應(yīng)用于實(shí)際，因此，找到網(wǎng)損修正系數(shù)快速準(zhǔn)確的計(jì)算方法，是一個(gè)亟需解決的關(guān)鍵問(wèn)題。

1 節(jié)能發(fā)電調(diào)度中發(fā)電組合方案區(qū)域內(nèi)優(yōu)化原則和流程

1.1 機(jī)組發(fā)電組合方案區(qū)域內(nèi)優(yōu)化的原則

《節(jié)能發(fā)電調(diào)度辦法》中規(guī)定，節(jié)能發(fā)電調(diào)度下，按以下方式對(duì)次日發(fā)電組合方案進(jìn)行區(qū)域內(nèi)優(yōu)化。

區(qū)域電力調(diào)度機(jī)構(gòu)在各省(區(qū)、市)機(jī)組次日發(fā)電組合方案的基礎(chǔ)上，依據(jù)本區(qū)域內(nèi)各省(區(qū)、市)排序表、各機(jī)組申報(bào)的可調(diào)發(fā)電能力、跨省輸電聯(lián)絡(luò)線(xiàn)的輸送電能力和網(wǎng)損，進(jìn)一步優(yōu)化調(diào)整本區(qū)域內(nèi)發(fā)電機(jī)組的啟停機(jī)方式。即：進(jìn)一步對(duì)各省(區(qū)、市)邊際機(jī)組(被調(diào)用的最后1臺(tái)機(jī)組)考慮網(wǎng)損因素后的供電煤耗率(簡(jiǎn)稱(chēng)邊際供電煤耗率)進(jìn)行比較，對(duì)邊際供電煤耗率較高的省(區(qū)、市)依次調(diào)整安排停機(jī)，對(duì)邊際供電煤耗率較低的省(區(qū)、市)依次調(diào)整安排啟機(jī)，直至區(qū)域中各省(區(qū)、市)的邊際供電煤耗率趨同，或跨省(區(qū)、市)輸電聯(lián)絡(luò)線(xiàn)達(dá)到輸送容量的極限。

1.2 機(jī)組發(fā)電組合方案區(qū)域內(nèi)優(yōu)化過(guò)程

按照節(jié)能發(fā)電調(diào)度辦法，發(fā)電組合方案區(qū)域內(nèi)優(yōu)化過(guò)程如下。

a. 對(duì)全網(wǎng)的發(fā)電機(jī)組按煤耗率由小到大分省進(jìn)行排序。

b. 根據(jù)年平均方式下各省的總負(fù)荷需求，按照排序表，綜合考慮備用要求，確定各省所需發(fā)電容量，并分省確定次日機(jī)組發(fā)電組合的方案。

c. 以年平均運(yùn)行方式為基礎(chǔ)，利用對(duì)潮流計(jì)算中雅克比矩陣進(jìn)行轉(zhuǎn)置求逆的方法[4-6]，計(jì)算各發(fā)電機(jī)組的網(wǎng)損微增率。

d. 利用(c)中得到的網(wǎng)損微增率對(duì)機(jī)組的供電煤耗率進(jìn)行修正。

e. 在(b)所確定的發(fā)電組合方案的基礎(chǔ)上，根據(jù)《辦法》(四)的原則，考慮網(wǎng)損因素，對(duì)區(qū)域內(nèi)發(fā)電機(jī)組的啟停機(jī)方式進(jìn)行優(yōu)化調(diào)整，即：進(jìn)一步對(duì)各省的邊際機(jī)組考慮了網(wǎng)損因素后的供電煤耗率(利用網(wǎng)損微增率修正以后的供電煤耗率)進(jìn)行比較，對(duì)邊際供電煤耗率較高的省依次調(diào)整安排停機(jī)，對(duì)邊際供電煤耗率較低的省依次調(diào)整安排啟機(jī)，直至各省的邊際供電煤耗率趨同或跨省輸電聯(lián)絡(luò)線(xiàn)達(dá)到輸送容量的極限。

2 對(duì)機(jī)組供電煤耗率進(jìn)行網(wǎng)損修正的必要性

按照傳統(tǒng)的調(diào)度規(guī)則，發(fā)電機(jī)組每年發(fā)電量均由政府下達(dá)指標(biāo)確定，其供電煤耗率僅作為日前機(jī)組組合和負(fù)荷分配等短期調(diào)度的依據(jù)，并不影響發(fā)電企業(yè)年度電量指標(biāo)。實(shí)行節(jié)能調(diào)度后，供電煤耗率成為發(fā)電機(jī)組能否上網(wǎng)發(fā)電的唯一依據(jù)。而利用網(wǎng)損微增率對(duì)發(fā)電機(jī)組供電煤耗率進(jìn)行修正，這其中引起的微小變化都可能會(huì)影響發(fā)電機(jī)組上網(wǎng)發(fā)電與否，關(guān)系到發(fā)電企業(yè)的重大利益。

節(jié)能調(diào)度的目標(biāo)是全網(wǎng)整體供電煤耗率最低。由于各發(fā)電機(jī)組電氣位置不同，其輸出功率在系統(tǒng)中引起的損耗各異，在發(fā)電調(diào)度中如直接利用發(fā)電機(jī)組的耗量進(jìn)行調(diào)度，所得到的結(jié)果不是最優(yōu)的，即全網(wǎng)的供電煤耗率不是最低。只有按照網(wǎng)損修正后的耗量進(jìn)行調(diào)度才能使全網(wǎng)的能量消耗最少，所以，在進(jìn)行發(fā)電調(diào)度前必須要對(duì)各個(gè)發(fā)電機(jī)組的供電煤耗率進(jìn)行公平合理的修正。

3 煤耗率網(wǎng)損修正系數(shù)的準(zhǔn)確計(jì)算

3.1 供電煤耗率網(wǎng)損修正數(shù)學(xué)模型

發(fā)電機(jī)組供電煤耗率的網(wǎng)損修正表達(dá)式如下：

(1)

網(wǎng)損微增系數(shù)是本文中提出的一個(gè)概念，表示機(jī)組一天內(nèi)網(wǎng)損微增率的平均水平。

3.2 供電煤耗率網(wǎng)損修正系數(shù)的準(zhǔn)確計(jì)算模型及計(jì)算方法

根據(jù)網(wǎng)損修正的數(shù)學(xué)模型，從中可以得到網(wǎng)損修正系數(shù)的計(jì)算公式為

ρ=1+(?Ploss/?Pi)AVE

(2)

式中：ρ為發(fā)電機(jī)組i的網(wǎng)損修正系數(shù)；(?Ploss/?Pi)AVE為發(fā)電機(jī)組i的網(wǎng)損微增系數(shù)。

由式(2)可以看出，網(wǎng)損微增系數(shù)是唯一的一個(gè)變量，因此如何計(jì)算網(wǎng)損修正系數(shù)的關(guān)鍵就是如何計(jì)算網(wǎng)損微增系數(shù)。

根據(jù)數(shù)學(xué)復(fù)合偏導(dǎo)公式，可得[7-8]：

(3)

式(3)中的方陣正為潮流計(jì)算中的雅可比矩陣。

而系統(tǒng)的有功損耗為

(4)

可得：

(5)

(6)

網(wǎng)損修正系數(shù)的精確計(jì)算步驟如下：

a. 根據(jù)潮流計(jì)算的結(jié)果得到雅克比矩陣，并對(duì)其轉(zhuǎn)置求逆；

b. 根據(jù)式(5)和式(6)得到相應(yīng)的列向量；

c. 用①得到的矩陣和②得到的列向量相乘，最終得到總有功損耗對(duì)各節(jié)點(diǎn)有功功率的偏導(dǎo)，從而提取出總有功損耗對(duì)機(jī)組有功的偏導(dǎo)，即發(fā)電機(jī)組的網(wǎng)損微增率。

d. 計(jì)算機(jī)組各時(shí)段網(wǎng)損微增率的平均值，即為網(wǎng)損微增系數(shù)。

4 網(wǎng)損修正系數(shù)的快速計(jì)算

4.1 網(wǎng)損修正系數(shù)快速計(jì)算方法

根據(jù)網(wǎng)損修正系數(shù)的準(zhǔn)確計(jì)算方法，結(jié)合經(jīng)驗(yàn)，提出以下3種快速計(jì)算方法，用來(lái)簡(jiǎn)便計(jì)算發(fā)電機(jī)組的網(wǎng)損微增系數(shù)，進(jìn)而得到發(fā)電機(jī)組的網(wǎng)損修正系數(shù)。

方法1：將總負(fù)荷等于平均總負(fù)荷時(shí)段的網(wǎng)損微增率取算術(shù)平均值作為網(wǎng)損微增系數(shù)。

方法2：將日平均方式下計(jì)算得到的網(wǎng)損微增率作為網(wǎng)損微增系數(shù)。

方法3：將峰谷方式下計(jì)算得到的網(wǎng)損微增率加權(quán)平均作為網(wǎng)損微增系數(shù)。

4.2 網(wǎng)損修正系數(shù)快速計(jì)算方法比較及驗(yàn)證

首先在我國(guó)某區(qū)域電網(wǎng)的節(jié)點(diǎn)負(fù)荷曲線(xiàn)中總結(jié)出7種典型曲線(xiàn)，并按照一定的比例分配給IEEE30個(gè)節(jié)點(diǎn)算例[9]的各個(gè)節(jié)點(diǎn)上，得到一個(gè)30個(gè)節(jié)點(diǎn)系統(tǒng)96個(gè)時(shí)段的基礎(chǔ)數(shù)據(jù)。默認(rèn)6臺(tái)機(jī)組全部上網(wǎng)發(fā)電，按照等微增率原則，采用夾板二分法，并綜合考慮機(jī)組爬坡率及最小技術(shù)出力的約束，對(duì)發(fā)電機(jī)組進(jìn)行負(fù)荷分配，從而得到算例系統(tǒng)中6臺(tái)機(jī)組所在節(jié)點(diǎn)的發(fā)電有功功率。選取其中的1號(hào)節(jié)點(diǎn)為平衡節(jié)點(diǎn)，對(duì)系統(tǒng)96個(gè)時(shí)段中的每個(gè)時(shí)段都分別進(jìn)行一次潮流計(jì)算，計(jì)算結(jié)果用來(lái)計(jì)算每臺(tái)機(jī)組每個(gè)時(shí)段的網(wǎng)損微增率。

利用3.2節(jié)中的方法計(jì)算每臺(tái)機(jī)組網(wǎng)損微增系數(shù)及網(wǎng)損修正系數(shù)的準(zhǔn)確值，見(jiàn)表1。

分別按照3種快速計(jì)算方法計(jì)算算例系統(tǒng)中6臺(tái)機(jī)組的網(wǎng)損微增系數(shù)及網(wǎng)損修正系數(shù)，結(jié)果見(jiàn)表2—表4所示。

表1 各機(jī)組網(wǎng)損修正系數(shù)準(zhǔn)確值

表2 方法1網(wǎng)損修正系數(shù)計(jì)算結(jié)果

表3 方法2網(wǎng)損修正系數(shù)計(jì)算結(jié)果

表4 方法3網(wǎng)損修正系數(shù)計(jì)算結(jié)果

比較3種計(jì)算方法，可以看到，按照方法2和方法3，準(zhǔn)確率很低，不能應(yīng)用于實(shí)踐當(dāng)中。而按照方法1，即將總負(fù)荷等于總負(fù)荷平均值的時(shí)段的網(wǎng)損微增率取算術(shù)平均作為網(wǎng)損微增系數(shù)，與網(wǎng)損微增系數(shù)的準(zhǔn)確值最為接近，且誤差在10%以?xún)?nèi)，可以作為網(wǎng)損微增系數(shù)的快速計(jì)算方法。

4.3 網(wǎng)損修正系數(shù)快速計(jì)算方法原理分析

對(duì)于每臺(tái)機(jī)組，分別將時(shí)段作為橫坐標(biāo)，將網(wǎng)損微增率作為縱坐標(biāo)，作出其96個(gè)時(shí)段的網(wǎng)損微增率曲線(xiàn)，并畫(huà)出網(wǎng)損微增系數(shù)準(zhǔn)確值橫線(xiàn)，與前一曲線(xiàn)交于幾個(gè)點(diǎn)，如圖1所示。作出96點(diǎn)的總負(fù)荷曲線(xiàn)，并作出平均負(fù)荷的直線(xiàn)，與其相交，如圖2所示。對(duì)照?qǐng)D1和圖2，發(fā)現(xiàn)每臺(tái)機(jī)組的網(wǎng)損微增率曲線(xiàn)與網(wǎng)損微增系數(shù)準(zhǔn)確值橫線(xiàn)的交點(diǎn)及總負(fù)荷曲線(xiàn)與平均負(fù)荷線(xiàn)的交點(diǎn)非常接近。

圖1 各機(jī)組的網(wǎng)損微增率變化曲線(xiàn)

圖2 系統(tǒng)總負(fù)荷變化曲線(xiàn)

經(jīng)過(guò)分析，按照方法1，即將總負(fù)荷等于總負(fù)荷平均值時(shí)段的網(wǎng)損微增率取算術(shù)平均作為網(wǎng)損微增系數(shù)具有較高的準(zhǔn)確率，是有一定依據(jù)的。在負(fù)荷分配的過(guò)程中，采用的是等微增率的原則，而無(wú)論是按照這種原則，還是按照等比例原則，得到的結(jié)果都將是在負(fù)荷高峰時(shí)段，各臺(tái)機(jī)組的出力都相應(yīng)較大，而在負(fù)荷低谷時(shí)段，各臺(tái)機(jī)組的出力也相應(yīng)較小，這就造成了各臺(tái)機(jī)組的發(fā)電功率曲線(xiàn)，以及總負(fù)荷曲線(xiàn)在一定意義上具有某種相似性。而理論及經(jīng)驗(yàn)均表明，對(duì)于同一臺(tái)機(jī)組，在趨勢(shì)上，網(wǎng)損微增率的大小與其出力的大小成正比關(guān)系，即當(dāng)其出力較大時(shí)，網(wǎng)損微增率也較大。因此，各臺(tái)機(jī)組的網(wǎng)損微增率曲線(xiàn)也具有某種相似性，并且與總負(fù)荷曲線(xiàn)存在某種對(duì)應(yīng)關(guān)系。經(jīng)過(guò)驗(yàn)證，出現(xiàn)圖1及圖2中的對(duì)應(yīng)關(guān)系，即每臺(tái)機(jī)組的網(wǎng)損微增率曲線(xiàn)與網(wǎng)損微增系數(shù)準(zhǔn)確值橫線(xiàn)的交點(diǎn)與總負(fù)荷曲線(xiàn)與平均負(fù)荷線(xiàn)的交點(diǎn)非常接近，即將總負(fù)荷等于平均總負(fù)荷時(shí)段的網(wǎng)損微增率取算術(shù)平均作為網(wǎng)損微增系數(shù)，準(zhǔn)確率非常高，為網(wǎng)損修正系數(shù)的計(jì)算提供了一種簡(jiǎn)單可行的方式。

5 結(jié)論

a. 由發(fā)電機(jī)組煤耗率的網(wǎng)損修正模型，引出了網(wǎng)損微增系數(shù)這一概念，并將網(wǎng)損修正系數(shù)的計(jì)算問(wèn)題轉(zhuǎn)化為網(wǎng)損微增系數(shù)的計(jì)算問(wèn)題。

b. 給出了網(wǎng)損微增系數(shù)準(zhǔn)確值的計(jì)算方法，在明確了實(shí)用過(guò)程中準(zhǔn)確計(jì)算網(wǎng)損微增系數(shù)的不可行性的基礎(chǔ)上，提出了網(wǎng)損微增系數(shù)需要一種簡(jiǎn)便可行的快速計(jì)算方法。

c. 提出了將總負(fù)荷等于總負(fù)荷平均值時(shí)段的網(wǎng)損微增率取算術(shù)平均作為網(wǎng)損微增系數(shù)，將日平均方式下計(jì)算得到的網(wǎng)損微增率作為網(wǎng)損微增系數(shù)，將峰谷方式下計(jì)算得到的網(wǎng)損微增率加權(quán)平均作為網(wǎng)損微增系數(shù)，這3種網(wǎng)損微增系數(shù)的簡(jiǎn)便計(jì)算方法。以IEEE30個(gè)節(jié)點(diǎn)的算例為例，對(duì)3種計(jì)算方法比較驗(yàn)證得到，將總負(fù)荷等于總負(fù)荷平均值時(shí)段的網(wǎng)損微增率取算術(shù)平均作為網(wǎng)損微增系數(shù)誤差在10%以?xún)?nèi)，可以作為網(wǎng)損微增系數(shù)的簡(jiǎn)便計(jì)算方法，進(jìn)而得到了網(wǎng)損修正系數(shù)的快速計(jì)算方法。

參考文獻(xiàn)：

[1] 中國(guó)國(guó)務(wù)院．國(guó)務(wù)院關(guān)于印發(fā)節(jié)能減排綜合性工作方案的通知[Z]．北京，2007.

[2] 國(guó)家發(fā)展改革委員會(huì)．節(jié)能發(fā)電調(diào)度改革試點(diǎn)將于近期啟動(dòng)[Z]．發(fā)改委網(wǎng)站，2007.

[3] 中國(guó)國(guó)務(wù)院．國(guó)務(wù)院辦公廳關(guān)于轉(zhuǎn)發(fā)發(fā)展改革委等部門(mén)節(jié)能發(fā)電調(diào)度辦法(試行)的通知[Z]．北京，2007.

[4] 楊淑英．電力系統(tǒng)概論[M]．北京：中國(guó)電力出版社，2007.

[5] 于爾鏗. 能量管理系統(tǒng)[M]. 北京：科學(xué)技術(shù)出版社，2001.

[6] 李文沅. 電力系統(tǒng)安全經(jīng)濟(jì)運(yùn)行—模型與方法[M]. 重慶：重慶大學(xué)出版社，1989.

[7] 陳 懇，李小銳，徐 敏．網(wǎng)損微增率新解法與轉(zhuǎn)置雅可比矩陣法用于有功優(yōu)化計(jì)算的比較[J]．中國(guó)電機(jī)工程學(xué)報(bào),2000,20(7):34-36.

[8] 張伯明，陳壽孫．高等電力網(wǎng)絡(luò)分析[M]．北京：清華大學(xué)出版社．1996.

[9] J Mutale,G Strbac,S Curcic and N Jenkins. Allocation of losses in distribution systems with embedded generation. IEE Proc-Gener. Trans. Distrib, 2000,147(1):7-14.