少環(huán)配電網(wǎng)節(jié)點(diǎn)邊際容量成本研究

錢佳琦，江岳文

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

?

少環(huán)配電網(wǎng)節(jié)點(diǎn)邊際容量成本研究

錢佳琦，江岳文

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

0　引　言

在現(xiàn)行的電價(jià)體制中，輸電和配電企業(yè)的供電容量成本(主要是輸電線路的容量成本)不能明確的在電價(jià)中體現(xiàn)出來(lái)，而只能轉(zhuǎn)移到電量電價(jià)或其他電價(jià)中去。隨著電力事業(yè)的發(fā)展和改革，輸配電可能進(jìn)行分離，負(fù)荷也在不斷的增大，如不明確這些設(shè)備的容量成本，將嚴(yán)重制約未來(lái)電網(wǎng)的發(fā)展。因此，在實(shí)際工程中，對(duì)供電設(shè)備容量成本確定具有非常重要的意義。文獻(xiàn)[1]提出配電網(wǎng)邊際容量成本(MDCC)，用來(lái)表示降低單位負(fù)荷而引起設(shè)備投資推遲而產(chǎn)生的收益。文章通過(guò)實(shí)例探究擴(kuò)容投資對(duì)MDCC的影響，并提出可以運(yùn)用MDCC指導(dǎo)設(shè)備擴(kuò)容投資規(guī)劃。文獻(xiàn)[2]提出了節(jié)點(diǎn)長(zhǎng)期增量成本(LRIC)的概念，用來(lái)表示負(fù)荷增加引起支路擴(kuò)容投資提前而產(chǎn)生的成本。文獻(xiàn)[3-4]在LRIC的基礎(chǔ)上，提出節(jié)點(diǎn)邊際容量成本(LMCC)，應(yīng)用于分布式電源的配電網(wǎng)規(guī)劃中。但文獻(xiàn)[2-4]均基于有功潮流進(jìn)行計(jì)算研究，忽略了網(wǎng)損和無(wú)功功率的影響。文獻(xiàn)[5]基于交流潮流計(jì)算建立了考慮無(wú)功功率和網(wǎng)損影響的LMCC模型，并提出了支路邊際容量成本(BMCC)的概念，用于評(píng)價(jià)支路的利用率和充裕度，并基于LMCC重點(diǎn)介紹了電網(wǎng)公司產(chǎn)權(quán)的分布式電源規(guī)劃。但對(duì)于LMCC的計(jì)算研究都是基于開環(huán)的輻射型配電網(wǎng)進(jìn)行的。隨著社會(huì)的進(jìn)步，用戶對(duì)用電質(zhì)量要求的不斷提高，尤其為了保證某些專用用戶的電能質(zhì)量，某些配電網(wǎng)會(huì)在實(shí)際運(yùn)行時(shí)閉合某些聯(lián)絡(luò)線，形成少環(huán)結(jié)構(gòu)的配電網(wǎng)絡(luò)。

基于以上文獻(xiàn)，本文根據(jù)少環(huán)配電網(wǎng)的結(jié)構(gòu)特點(diǎn)，提出分流系數(shù)矩陣的概念和計(jì)算方法，結(jié)合推導(dǎo)得到節(jié)點(diǎn)功率對(duì)支路電流的靈敏度，建立靈敏度系數(shù)矩陣，并以此求取線路擴(kuò)容時(shí)間和單位負(fù)荷增量成本矩陣，從而計(jì)算得到配電網(wǎng)LMCC。運(yùn)用IEEE33節(jié)點(diǎn)算例，計(jì)算結(jié)果對(duì)比改進(jìn)后的基于雅克比矩陣的節(jié)點(diǎn)-支路功率靈敏度系數(shù)矩陣計(jì)算LMCC的方法，驗(yàn)證本文方法的有效性。

1　少環(huán)配電網(wǎng)靈敏度系數(shù)矩陣

1.1分流系數(shù)矩陣

如圖1的少環(huán)配電網(wǎng)，由于存在環(huán)網(wǎng)，則無(wú)法通過(guò)形成離散的0-1變量的節(jié)點(diǎn)-支路網(wǎng)絡(luò)關(guān)聯(lián)矩陣來(lái)量化表示節(jié)點(diǎn)負(fù)荷增加對(duì)支路電流的影響。

圖1　少環(huán)配電網(wǎng)示例

因此本文提出建立分流系數(shù)矩陣，表示少環(huán)配電網(wǎng)節(jié)點(diǎn)負(fù)荷變化時(shí)各支路潮流變化的分流情況。

對(duì)于任意網(wǎng)絡(luò)有

(1)

式中：U表示節(jié)點(diǎn)電壓列向量；Z表示節(jié)點(diǎn)阻抗矩陣；I表示節(jié)點(diǎn)注入電流列向量。

展開得

(2)

式中：Ui(i=1,…,N)表示節(jié)點(diǎn)i的電壓；Zij(i,j=1,…,N)表示節(jié)點(diǎn)阻抗矩陣中的元素；Ii(i=1,…,N)表示節(jié)點(diǎn)i的注入電流。

設(shè)支路k的起始節(jié)點(diǎn)為s，末節(jié)點(diǎn)為e，則

(3)

因此可得分流系數(shù)：

(4)

列成矩陣形式即可得分流系數(shù)矩陣：

(5)

式中：b表示支路數(shù)；N表示節(jié)點(diǎn)數(shù)。

因此，通過(guò)分流系數(shù)矩陣，即可將少環(huán)配電網(wǎng)等效成輻射型的配電網(wǎng)。當(dāng)對(duì)圖1中的節(jié)點(diǎn)3進(jìn)行網(wǎng)絡(luò)等效時(shí)，可得等效的輻射型配電網(wǎng)如圖2所示。通過(guò)網(wǎng)絡(luò)等效，原來(lái)的節(jié)點(diǎn)3等效成了節(jié)點(diǎn)3(1)和節(jié)點(diǎn)3(2)。節(jié)點(diǎn)3(1)和節(jié)點(diǎn)3(2)的負(fù)荷大小分別為節(jié)點(diǎn)3的β73倍和β33倍。當(dāng)節(jié)點(diǎn)3負(fù)荷變化ΔS時(shí)，引起支路3和支路7潮流變化分別為β33ΔS和β73ΔS。且根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)圖可知：

(6)

圖2　等效輻射型配電網(wǎng)

依此類推，可根據(jù)節(jié)點(diǎn)的不同，可將少環(huán)配電網(wǎng)等效成相應(yīng)的輻射型配電網(wǎng)，便于后續(xù)進(jìn)行節(jié)點(diǎn)功率-支路電流靈敏度的推導(dǎo)。

而在分析純輻射型的配電網(wǎng)時(shí)，計(jì)算得到的分流系數(shù)矩陣即為只包含0-1變量的上三角矩陣。矩陣的行代表對(duì)應(yīng)支路的下游節(jié)點(diǎn)，其列代表對(duì)應(yīng)節(jié)點(diǎn)的上游支路。如圖1，斷開支路7，得到輻射型配電網(wǎng)，其分流系數(shù)矩陣：

(7)

1.2節(jié)點(diǎn)功率-支路電流靈敏度

已知某配電網(wǎng)某段線路如圖3所示。當(dāng)某段線路下游的節(jié)點(diǎn)i負(fù)荷功率變化ΔSi。忽略負(fù)荷增加時(shí)，相應(yīng)支路網(wǎng)損增加而引起的支路電流增量不等，即假設(shè)給節(jié)點(diǎn)i輸送功率的支路電流都將變化ΔIk?？傻?/p>

(8)

式中：B表示節(jié)點(diǎn)j到節(jié)點(diǎn)i之間所有支路的集合；Uj表示支路k的起始節(jié)點(diǎn)j的線電壓；Ik、Is分別表示支路k和支路s的電流；Si表示支路k下游某節(jié)點(diǎn)i的負(fù)荷功率；SΣ表示支路k下游除i節(jié)點(diǎn)負(fù)荷外，其他節(jié)點(diǎn)負(fù)荷功率與節(jié)點(diǎn)i下游支路網(wǎng)損之和。

圖3　某配電網(wǎng)部分線路示例

化簡(jiǎn)式(8)，則節(jié)點(diǎn)i注入功率對(duì)于支路k電流的靈敏度可以表示為

(9)

當(dāng)ΔSi→0時(shí)，

(10)

根據(jù)式(10)，即可建立節(jié)點(diǎn)功率-支路電流靈敏度矩陣。由于矩陣元素與當(dāng)前網(wǎng)絡(luò)負(fù)荷分布有關(guān)，可表示為

(11)

式中：S0表示在現(xiàn)網(wǎng)絡(luò)結(jié)構(gòu)下負(fù)荷的分布情況。

2　基于雅克比矩陣靈敏度系數(shù)矩陣

本文對(duì)文獻(xiàn)[6]用于電力系統(tǒng)有功安全校正控制而提出的基于雅克比矩陣的節(jié)點(diǎn)-支路潮流靈敏度的計(jì)算方法進(jìn)行一定的改進(jìn)，用于計(jì)算LMCC(后文簡(jiǎn)稱該方法為“Jacobi法”)，并將其作為對(duì)照方法與本文提出方法進(jìn)行對(duì)比。

在極坐標(biāo)下，運(yùn)用雅克比矩陣進(jìn)行潮流求解的方程為

(12)

(13)

式中:[ΔPk]=(ΔPij)b×1、[ΔQk]=(ΔQij)b×1分別代表有功、無(wú)功支路潮流增量矩陣； [ΔP]=(ΔP)(N-1)×1、[ΔQ]=(ΔQ)(N-1)×1分別代表除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的有功、無(wú)功功率增量矩陣。

根據(jù)式(13)可得

(14)

因此，節(jié)點(diǎn)i功率對(duì)支路k功率的靈敏度為

(15)

式中：φi為節(jié)點(diǎn)i的功率因數(shù)角。

寫成矩陣形式可得

(16)

3　配電網(wǎng)LMCC

3.1支路擴(kuò)容時(shí)間

假設(shè)負(fù)荷每年以恒功率因數(shù)等比增長(zhǎng)，因此對(duì)于支路k有

(17)

(18)

化簡(jiǎn)公式(18)可得擴(kuò)容時(shí)間τk：

(19)

式中：Ik，CAP表示支路k的最大載流量；d表示負(fù)荷年增長(zhǎng)率。

3.2LMCC

考慮擴(kuò)容方式為增加同一型號(hào)和容量的供電設(shè)備并聯(lián)運(yùn)行，則貼現(xiàn)值為

(20)

式中：Ck表示支路k的擴(kuò)容成本；Ck,PV表示支路k投資現(xiàn)值；r表示折現(xiàn)率。

根據(jù)公式(19)和(20)可知，節(jié)點(diǎn)i的負(fù)荷以恒定的功率因數(shù)增加ΔSi時(shí)，擴(kuò)容時(shí)間τk將減小，設(shè)備的折現(xiàn)值則增大。因此，節(jié)點(diǎn)i對(duì)于支路k的單位負(fù)荷增長(zhǎng)成本為

(21)

式中：Ck,PV*表示負(fù)荷變化引起擴(kuò)容時(shí)間改變后，支路k新的投資現(xiàn)值。

當(dāng)ΔSi→0時(shí)，

(22)

結(jié)合式(19)、(20)和(22)，可得節(jié)點(diǎn)i對(duì)于支路k的單位負(fù)荷增長(zhǎng)成本年值為

(23)

式中：R表示等年值系數(shù)。

將Δcki列成矩陣形式：

(24)

則節(jié)點(diǎn)i的LMCC為節(jié)點(diǎn)i對(duì)所有支持該節(jié)點(diǎn)功率輸送的支路的單位增量成本之和，即

(25)

4　算例分析

4.1算例介紹

以IEEE33配電網(wǎng)為例，連接兩條聯(lián)絡(luò)線，如圖3所示。根據(jù)原始線路數(shù)據(jù)設(shè)定相應(yīng)的線路投資成本，節(jié)點(diǎn)0為發(fā)電機(jī)節(jié)點(diǎn)。

圖4　IEEE33節(jié)點(diǎn)配電系統(tǒng)

設(shè)負(fù)荷增長(zhǎng)率為3%，折現(xiàn)率為8%，線路投資回報(bào)年限為40 a。

4.2擴(kuò)容時(shí)間驗(yàn)證

分別對(duì)IEEE33配電網(wǎng)輻射型、連接31-8聯(lián)絡(luò)線形成單環(huán)以及連接31-8和25-14聯(lián)絡(luò)線形成雙環(huán)這3種情況下的配電網(wǎng)進(jìn)行支路擴(kuò)容時(shí)間的計(jì)算，對(duì)比本文方法與“Jacobi法”，結(jié)果如圖5所示。

圖5　擴(kuò)容時(shí)間對(duì)比

通過(guò)對(duì)比得出：相對(duì)于“Jacobi法”，用本文方法計(jì)算得到各支路擴(kuò)容時(shí)間的最大誤差分別為3.23%、3.55%和2.26%，都在3%左右。

4.3LMCC計(jì)算結(jié)果驗(yàn)證

分別對(duì)IEEE33配電網(wǎng)輻射型、連接31-8聯(lián)絡(luò)線形成單環(huán)以及連接31-8和25-14聯(lián)絡(luò)線形成雙環(huán)這3種情況下的配電網(wǎng)進(jìn)行LMCC的計(jì)算，對(duì)比本文方法與“Jacobi法”，結(jié)果如圖6所示。

圖6　LMCC對(duì)比

通過(guò)對(duì)比得出：相對(duì)于“Jacobi法”，用本文方法計(jì)算得到各節(jié)點(diǎn)LMCC的最大誤差分別為4.41%、6.31%和6.81%，這些誤差都在允許的范圍之內(nèi)。

5　結(jié)　論

隨著市場(chǎng)化進(jìn)程的不斷推進(jìn)和電力市場(chǎng)的逐步發(fā)展、推廣，基于容量成本的電價(jià)體制也將不斷完善。本文提出的少環(huán)配電網(wǎng)節(jié)點(diǎn)邊際容量成本的研究用于反映配電網(wǎng)供電設(shè)備的容量成本，主要結(jié)論如下：

① 考慮了含有環(huán)網(wǎng)情況下的配電網(wǎng)LMCC的計(jì)算，更能符合現(xiàn)實(shí)生產(chǎn)運(yùn)行的需要，也是適應(yīng)未來(lái)配電網(wǎng)發(fā)展趨勢(shì)的容量成本研究方法。

② 本文提出的分流系數(shù)矩陣，能夠真實(shí)地反映節(jié)點(diǎn)負(fù)荷變化時(shí)各支路潮流變化的分流情況，基于此得到的配電網(wǎng)LMCC也更準(zhǔn)確。由于配電系統(tǒng)的不斷發(fā)展，運(yùn)用能夠體現(xiàn)真實(shí)潮流狀態(tài)的計(jì)算方法是必不可少的。

③ 由于牛拉法在電力系統(tǒng)潮流計(jì)算有著廣泛的適用性，使得本文改進(jìn)后的基于雅克比矩陣的靈敏度計(jì)算方法廣泛適用于電力網(wǎng)絡(luò)的LMCC計(jì)算分析，但其在計(jì)算過(guò)程中進(jìn)行了多次的矩陣計(jì)算，容易造成了誤差的累積。因此，對(duì)于網(wǎng)絡(luò)結(jié)構(gòu)相對(duì)明確的配電系統(tǒng)，更適用于運(yùn)用基于分流系數(shù)矩陣的配電網(wǎng)LMCC計(jì)算方法。

配電網(wǎng)LMCC反映了配電網(wǎng)各節(jié)點(diǎn)對(duì)供電容量的占用情況和供電設(shè)備的利用程度，對(duì)容量電價(jià)制定、分布式電源規(guī)劃以及配電網(wǎng)需求側(cè)管理等多方面都有著十分重要的作用。

[1]GRAYSON H, WOO C K, BRIAN H, et al. Variations in area and time specific marginal capacity costs of electricity distribution[J]. IEEE Trans on Power Systems, 1998, 13 (2): 560-565.

[2]Li F, Tolley D L. Long-Run Incremental cost pricing based on unused capacity[J]. IEEE Trans on Power Systems, 2007, 22(4): 1683-1689.

[3]歐陽(yáng)武，程浩忠，張秀彬等.配電網(wǎng)節(jié)點(diǎn)邊際容量成本及應(yīng)用[J].電力系統(tǒng)自動(dòng)化，2009，33(1)：18-21.

[4]歐陽(yáng)武，程浩忠，張秀彬.基于配電網(wǎng)節(jié)點(diǎn)邊際容量成本的分布式電源規(guī)劃[J].電力系統(tǒng)自動(dòng)化，2009，33(18)：28-31.

[5]任龍霞，雷金勇，甘德強(qiáng).計(jì)及配電網(wǎng)節(jié)點(diǎn)/支路邊際容量成本的分布式電源規(guī)劃[J].電力系統(tǒng)自動(dòng)化，2010，34(21)：34-39.

[6]劉文穎，徐鵬，梁才等.一種計(jì)及交叉權(quán)重的有功校正控制算法[J].電工技術(shù)學(xué)報(bào)，2014，29(8)：281-288.

[7]王承民，蔣傳文.基于支路電流狀態(tài)變量的靈敏度分析方法研究[J].電工技術(shù)學(xué)報(bào)，2006，21(1)：42-46.

[8]荊朝霞，段獻(xiàn)忠，文福拴等.輸電系統(tǒng)固定成本分?jǐn)倖?wèn)題[J].電力系統(tǒng)自動(dòng)化，2003，27(15)：84-89.

[9]王艷松，何新霞.基于網(wǎng)絡(luò)模型的節(jié)點(diǎn)阻抗矩陣的新算法[J].中國(guó)石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2000，24(3)：98-99.

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

[11]車仁飛，李仁俊.一種少環(huán)配電網(wǎng)三相潮流計(jì)算新方法[J].中國(guó)電機(jī)工程學(xué)報(bào)，2003，23(1)：74-79.

[12]顏麗.基于電流分布的電網(wǎng)功率分布因子的計(jì)算方法[D].北京：華北電力大學(xué)，2010.

[13]Woo C K, DEBRA L Z, REN O, et al. Marginal capacity costs of electricity distribution and demand for distributed generation[J]. The Energy Journal, 1995, 16(2): 111-120.

(責(zé)任編輯：林海文)

Research on the Locational Marginal Capacity Cost of Weakly Meshed Distribution Network

QIAN Jiaqi, JIANG Yuewen

(College of Electrical Engineering and Automation, Fuzhou University, Fuzhou 350108,China)

配電網(wǎng)負(fù)荷增加將導(dǎo)致供電設(shè)備擴(kuò)容提前，從而提高擴(kuò)容投資現(xiàn)值。針對(duì)少環(huán)配電網(wǎng)的網(wǎng)絡(luò)結(jié)構(gòu)，提出了能夠反映支路潮流分流情況的分流系數(shù)矩陣的概念及其計(jì)算方法，并結(jié)合推導(dǎo)而得的節(jié)點(diǎn)注入功率-支路電流靈敏度，建立靈敏度系數(shù)矩陣，用于計(jì)算得到可以表示配電網(wǎng)各節(jié)點(diǎn)占用供電容量成本大小的節(jié)點(diǎn)邊際容量成本(locational marginal capacity cost, LMCC)。進(jìn)一步改進(jìn)了基于雅克比矩陣的節(jié)點(diǎn)-支路潮流靈敏度系數(shù)矩陣的計(jì)算方法，用于計(jì)算LMCC，并作為對(duì)照方法在IEEE33配電系統(tǒng)中進(jìn)行測(cè)試，對(duì)比驗(yàn)證了本文方法的有效性。

少環(huán)配電網(wǎng)；分流系數(shù)矩陣；靈敏度系數(shù)矩陣；節(jié)點(diǎn)邊際容量成本

The increasing of the load in the distribution network will advance the expansion time of the power supply equipment, and improve the value of expansion investment. According to the structure of weakly meshed distribution network, the conception and calculation method of current distribution coefficient matrix, which can reflect the distribution of the power flow, are proposed. The sensitivity coefficient matrix is built by combining the sensitivity of node power injection-branch current and current distribution coefficient matrix. And the locational marginal capacity cost (LMCC), which represents the capacity cost of each node, can be calculated. The sensitivity coefficient matrix depending on the Jacobi matrix is improved to calculate the LMCC, and the IEEE33 distribution system is used as test example to verify the feasibility of the method proposed in this paper.

weakly meshed distribution network; current distribution coefficient matrix; sensitivity coefficient matrix; locational marginal capacity cost

1007-2322(2016)03-0043-05

A

TM744

福建省自然科學(xué)基金項(xiàng)目(2013J01176)

2015-03-17

錢佳琦(1991—)，男，碩士研究生，研究方向?yàn)榉植际诫娫窗l(fā)電與電力市場(chǎng)，E-mail：n130127020@fzu.edu.cn；

江岳文(1977—)，女，副教授，研究方向?yàn)殡娏ο到y(tǒng)優(yōu)化運(yùn)行與電力市場(chǎng)。