于繼來,郭鈺鋒,盛玉和,張偉,王松巖
(1. 哈爾濱工業(yè)大學(xué) 電氣工程及自動化學(xué)院,黑龍江 哈爾濱 150001;2. 國網(wǎng)黑龍江省電力有限公司電力科學(xué)研究院,黑龍江 哈爾濱 150030)
電力市場是一種電能交易合同必須緊密依托電網(wǎng)輸配電功能才能真正付諸實施的特殊市場。此特性決定了由發(fā)用電商基于某種交易模式及其規(guī)則形成的電能交易合同,在局限于已有電力系統(tǒng)路網(wǎng)結(jié)構(gòu)并受到電網(wǎng)實時運行方式變化的干擾下,必然會出現(xiàn)不以交易主體主觀意志或商業(yè)意愿為轉(zhuǎn)移、與真正電能分配格局產(chǎn)生較大差異甚至迥異的現(xiàn)象。此現(xiàn)象背后的內(nèi)因如若揭示不明,差異程度如若不能客觀定量評估,則勢必造成對部分能量流-經(jīng)濟流分析結(jié)果的錯亂,由此形成不同主體間的長期不合理交叉補貼問題,影響市場公正性[1-3]。
在吸收和借鑒前人思想的基礎(chǔ)上,本文作者認為上述問題需要先后從2個關(guān)鍵點上取得突破:一是如何基于電力系統(tǒng)物理平臺提供的實時狀態(tài)信息,厘清發(fā)電商的上網(wǎng)電能是如何實時地配送給不同購電商的?或購電商的下網(wǎng)電能是如何實時地從多個發(fā)電商處吸納的?二是如何基于市場平臺提供的交易信息,將本已由合同或政策固化的市場均衡價格信號,折算為能夠跟隨電網(wǎng)運行方式改變的實時價格信號?
第1突破點實際上就是實時能流的正確解析問題,第2突破點就是實時價格的正確計算問題。
關(guān)于第1突破點,近20多年來,眾多學(xué)者從不同角度,圍繞電力網(wǎng)絡(luò)源流電氣分配關(guān)系問題,進行了大量研究。就大類而言,一類研究基于交流電路相量體系。其中,無論是側(cè)重物理解析[4-6],還是側(cè)重數(shù)學(xué)解析[7-9],一直受到有功、無功分量非線性耦合的困擾,尚未取得共識性[10-11]。另一類研究基于交流電路瞬時值體系[12],其在完善交流節(jié)點瞬時電流和瞬時功率比例共享定理證明過程的基礎(chǔ)上[13-14],提出了能夠正確解析交直流電力網(wǎng)絡(luò)源流電能分配關(guān)系的新算法并開發(fā)出相應(yīng)程序[15],此對促進第1點的突破,提供了有力的技術(shù)支持。
關(guān)于第2突破點,具有實用性的研究最早可追溯到20世紀80年代末文獻[16-17]提出了節(jié)點電價概念:在滿足電力系統(tǒng)運行約束和市場主體資源約束的條件下,某節(jié)點增加單位負荷需求時的邊際成本。節(jié)點邊際電價已在國內(nèi)外多個電力市場實踐中得到很好應(yīng)用,其對優(yōu)化電網(wǎng)運行方式、整合市場資源配給效益,發(fā)揮了積極作用。文獻[18]大力倡議將節(jié)點電價應(yīng)用于中國現(xiàn)貨市場并分析了其與市場均衡價格的關(guān)系。文獻[19]針對現(xiàn)貨電力市場,結(jié)合節(jié)點邊際電價和分區(qū)輸電電價概念,探討了能夠提供位置信號的電價機制問題。然而,節(jié)點邊際電價一般基于最優(yōu)潮流(optimal power flow, OPF)計算得出[17,20],技術(shù)實現(xiàn)上雖無太大困難,但在理論和應(yīng)用上,尚存在下述需特別注意的問題。
(1)節(jié)點邊際價格α只反映特定存量負荷需求Db(基點工況)再增加1個單位時,成本由原先的Cb增加的程度ΔC(=α),不能完全用作對變化后總需求Db+1的下網(wǎng)費進行實時核算的價格。若設(shè)實時需求為D(t),核算價格為β(t),則下網(wǎng)費為二者之積。對Db+1,下網(wǎng)價應(yīng)為(Cb+ΔC) /(Db+1),它只在α=Cb/Db條件下,才等于α。然而,對實時變化的電網(wǎng)運行狀態(tài),理論上此條件幾乎不成立,且條件偏離程度,與Db、Cb和α均有關(guān)。
(2)當α由OPF生成時,其為優(yōu)化系統(tǒng)、特別是對調(diào)控阻塞,提供非常重要的信號。但絕大多數(shù)OPF算法只取性能目標改善最大的方向?qū)?yōu),此方向與電網(wǎng)實時運行狀態(tài)實際變化方向很難完全吻合,因此,α也就不能真正代表邊際成本。
(3)應(yīng)用OPF計算時,從為系統(tǒng)性能調(diào)控提供最佳決策的角度,一般給出算法迭代若干次收斂后的最優(yōu)解。如此給出的邊際成本信號,更多地反映最優(yōu)需求點Dopt處的邊際成本信息,而不是基于Db、根據(jù)實際變化方向增加1個單位需求時的邊際成本。實際上,Dopt與Db可能離得較遠。由問題(1)知,在Db處,β(t)與α就存在偏差,如若換成Dopt,偏差可能更大。
因此,基于反映電網(wǎng)真實運行狀況的實時潮流方式(由在線潮流計算或在線量測系統(tǒng)采集獲得),而非OPF方式計算的實時價格β(t),才是真正符合第2突破點內(nèi)涵要求的信號。
綜上,本文在電能追蹤程序[15]正確提供能流信息基礎(chǔ)上,重點針對第2突破點,從實時能量流-經(jīng)濟流耦合角度,探討實時電價的正確計算問題。
基于瞬時值體系的電能追蹤方法[15],沿“解析瞬時電流成分解析瞬時功率成分積分計算周期電能成分”的路線,并考慮兼容現(xiàn)有關(guān)口式電能計量模式,先求出所有交/直流輸變電支路送端源與受端流間的電能分配關(guān)系,并記由送端源i與受端流j構(gòu)成的支路源-流對的直達能量成分信息向量為
圖1側(cè)重示意了支路源-流對利用交/直流支路進行能量配送的情況:在支路l送端(節(jié)點Bs內(nèi)側(cè)),產(chǎn)自源i的電能配裝在“i號貨卡”上(“發(fā)貨方”為源i),其中的Aij“電能包裹”標簽上的“收貨方”為流j,因中途損耗一部分,當其運抵受端(節(jié)點Br內(nèi)側(cè))時被流j簽收的“電能包裹”為aij。對流j,它除了可能簽收來自源i配送的包裹aij外,還可能收到來自其他送端源配送的包裹 a1j、···、ai-1,j、ai+1,j、···、amj,這些系列包裹a·j的累加量就是流j在受端(Br內(nèi)側(cè))總簽收的來自送端源集合的電能。
圖1 支路源-流對能量配送示意Fig. 1 Schematic diagram of energy distribution for branch source-flow pair
在掌握所有支路源-流對直達能量成分信息向量的基礎(chǔ)上,電能追蹤方法[15]可以基于任意網(wǎng)絡(luò)源-流對
圖2 網(wǎng)絡(luò)源-流對能量配送示意Fig. 2 Schematic diagram of energy distribution for the network source-flow pair
針對集合Ω匯總?cè)绫?所示的各子鏈直達能量信息,即可獲得
表1 子鏈直達能量成分信息構(gòu)成Table 1 Through-energy information of a sub-chain
電源和負荷設(shè)備單元一旦并網(wǎng)運行,就定會產(chǎn)生不同程度的輸配電業(yè)務(wù)需求,并接收電網(wǎng)平臺提供的中間服務(wù)。理論上,電網(wǎng)實時狀態(tài)的變化,意味著為市場主體提供服務(wù)的具體內(nèi)容隨之變化,也說明依托電網(wǎng)的能量流和經(jīng)濟流緊密耦合。因此,在分析電網(wǎng)源流間經(jīng)濟流分布信息時,需綜合考慮電網(wǎng)平臺提供的能流信息和市場平臺提供的交易信息間的耦合關(guān)系。一般的,能量流和經(jīng)濟流耦合的形式及其主要誘導(dǎo)因素如下。
(1)單個網(wǎng)絡(luò)源-流對間的能量流和經(jīng)濟流耦合。
由第1節(jié)知,單個網(wǎng)絡(luò)源-流對
由此可知,經(jīng)濟流主要受市場主體間的實時商業(yè)供求關(guān)系的影響,能量流主要受電網(wǎng)源-流對間的實時電氣配送關(guān)系的影響,且二者均隨電網(wǎng)實時狀態(tài)而變。單個網(wǎng)絡(luò)源-流對
(2)多個網(wǎng)絡(luò)源-流對間的能量流和經(jīng)濟流耦合。
對同一電網(wǎng)內(nèi)的多個網(wǎng)絡(luò)源-流對
綜合上述2類耦合形式可知,其原發(fā)因素均因電網(wǎng)實時狀態(tài)改變而引起,且針對第1類耦合形式分析的結(jié)果,是進一步分析多個網(wǎng)絡(luò)源-流對間的能量流和經(jīng)濟流耦合關(guān)系的基礎(chǔ)(因為第2類耦合關(guān)系可分解成多個第1類耦合關(guān)系進行分析)。
下文首先分析與單個網(wǎng)絡(luò)源-流對某一路徑子鏈直達能量成分關(guān)聯(lián)的經(jīng)濟流分布信息,然后再結(jié)合典型交易模式闡述與其關(guān)聯(lián)的多個網(wǎng)絡(luò)源-流對的經(jīng)濟流分布量化分析問題。
對實際電網(wǎng),基于第1節(jié)電能追蹤應(yīng)用模塊(只需調(diào)用現(xiàn)有EMS/DMS的在線潮流模塊計算結(jié)果)提供的能流信息和市場平臺提供的交易信息,完全可以應(yīng)用上述關(guān)系式在線計算用于核算關(guān)聯(lián)子鏈各種費用成分的實時電價。
繼續(xù)討論子鏈能耗特性對實時電價的影響以及2種計費方式的主要特征如下。
(1)鏈直達能量傳輸效率η們越高,意味著傳輸能耗越低。對mode I,在鏈首能量成分一定時,由式(8)提供的等效過網(wǎng)價格和由式(11)提供的等效政府性基金價格就會越高,2種價格隨η單調(diào)抬升的效應(yīng)傳導(dǎo)到用戶端,使得用戶等效購電總價也隨之升高(符合成本平衡回收原理);對mode II,在鏈末能量成分一定時,由式(8)提供的等效過網(wǎng)價格和由式(3)提供的等效上網(wǎng)電價就會越低,2種價格隨η單調(diào)走低的效應(yīng)傳導(dǎo)到用戶端,使得用戶等效購電總價cjnet也隨之降低。由此可見,和cjnet具有完全相反的隨η單調(diào)變化的特性。然而,必須透過此種表面“反調(diào)”特性注意到:雖然2種計費方式采用的能量值(鏈首電源端計量值和鏈末用戶端計量值)和總價格有所不同,但二者核算后的用戶總支出及其他主體(電源、電網(wǎng)、政府)收入均對應(yīng)相等。這說明,本文方法在考慮電網(wǎng)平臺能流信息和市場平臺交易信息耦合特性的基礎(chǔ)上推導(dǎo)出的實時電價模型及利用其核算后的各種費用成分,具有鏈首鏈末等值性(形象化的理解就是:分別由鏈首和鏈末兩撥“會計”核算的結(jié)果完全相同)。實際上,此結(jié)論還可從理論上推廣為:鏈首和鏈末、鏈各路段段首和段末(乃至鏈上的任一點)均可作為費用核算點,身處其中的“會計”只需利用所在點計量的直達能量成分和折算至所在點的實時價格進行核算,必能給出相同的財務(wù)結(jié)果(即不同核算點的“會計”不會因為“辦公地點”的差異而得出不同的財務(wù)結(jié)果)。此特性很重要,屬于電力與交易耦合系統(tǒng)的基本特性,本文將其歸納為如下的核算一致性原理。
核算一致性原理:對源-流對
原理的正確性不難由式(12)等價格關(guān)系式形成過程證明。限于篇幅,不再贅述。
上述原理可拓展至源-流對
(3)當線損費用獨立核算(不含在輸配電價中)且其費用需要由電源和負荷分攤時,表1基于物理系統(tǒng)給出的子鏈及其各路段能耗向源inet和流jnet分攤的結(jié)果,才會起直接作用。然而,在我國目前執(zhí)行中發(fā) 〔2015〕 9號文規(guī)定用戶購電價格由上網(wǎng)電價、輸配電價(含線損)和政府性基金三部分組成的費用核算機制下,由物理系統(tǒng)得出的能耗分攤信息,將不再用于除(1) ~(3)以外的其他相關(guān)經(jīng)濟流分析(文件已規(guī)定將線損費用統(tǒng)一核算到用戶端)??紤]到中國目前電力市場改革政策,本文暫不討論線損費用獨立核算問題。
就源-流對
式中:上標中的I和II分別對應(yīng)mode I和II,對應(yīng)的電能分別為集合Ω在源端配送的和在流端接收的總電能為
對mode I和II,式(17)中的4項資金流均對應(yīng)相等,仍然遵守核算一致性原理。
在實際應(yīng)用時,如果需要同時分析多個網(wǎng)絡(luò)源-流對
第2節(jié)分析了任意多個網(wǎng)絡(luò)源-流對
當發(fā)電商inet和購電商jnet間只簽有雙邊直接交易電能合同時,可以直接采用2.2節(jié)方法進行經(jīng)濟流及其財務(wù)核算。
更一般的,若發(fā)電商inet和NL個購電商同時簽有雙邊直接交易合同,且在某運行方式時,inet和NL個購電商所在電網(wǎng)節(jié)點間存在配送直達能量成分的路徑鏈集合群 Φ={Ω1, Ω2, ··· , ΩNL},則inet和第k個購電商間的經(jīng)濟流及其對應(yīng)的實時價格信號,只需依據(jù)子鏈集合Ωk提供的信息,采用2.2節(jié)方法計算獲得。在此基礎(chǔ)上,匯總集合群Φ的結(jié)果,可得發(fā)電商inet網(wǎng)售電總收入額為
與其對應(yīng)的平均網(wǎng)售電價格為
相對照地,若購電商jnet和NG個發(fā)電商同時簽有雙邊直接交易合同,且在某運行方式時,jnet和NG個發(fā)電商所在電網(wǎng)節(jié)點間存在配送直達能量成分的路徑鏈集合群 Θ={Ω1, Ω2, ··· , ΩNG},則jnet和第k個發(fā)電商間的經(jīng)濟流及其對應(yīng)的實時價格信號,同樣只需依據(jù)子鏈集合Ωk提供的信息,采用2.2節(jié)方法計算獲得。在此基礎(chǔ)上,匯總集合群Θ的結(jié)果,可得購電商jnet網(wǎng)購電總支出額為
與其對應(yīng)的平均網(wǎng)購電價格為
設(shè)發(fā)電商群{inet}(含NG個發(fā)電商)與某購電商群{jnet}(含NL個購電商)簽有雙多邊交易合同,且在某運行方式時,群{inet}和群{jnet}所在電網(wǎng)節(jié)點間存在配送直達能量成分的路徑鏈集合群, ··· ,},則群{inet}中第 kG個發(fā)電商和群{jnet}中第kL個購電商間的經(jīng)濟流及其對應(yīng)的實時價格信號,可依據(jù)子鏈集合ΩkGkL提供的信息,采用2.2節(jié)方法計算獲得。
依托新英格蘭10機39節(jié)點系統(tǒng)[21]模擬給出計算條件:系統(tǒng)頻率設(shè)為50 Hz,基準功率為100 MW;33號節(jié)點機組G33作為待分析網(wǎng)絡(luò)源(發(fā)電商),16號節(jié)點負荷L16作為待分析網(wǎng)絡(luò)流(購電商),二者構(gòu)成一對商業(yè)伙伴{G33,L16};設(shè)系統(tǒng)某運行方式持續(xù)1 h,該方式下與{G33, L16}關(guān)聯(lián)的所有電能配送流向及其路徑信息如圖3和表2所示。
圖3 商業(yè)伙伴{G33, L16}關(guān)聯(lián)的電能配送流向Fig. 3 Energy flow directions associated with {G33, L16}
由圖3和表2可見,在源-流對
另外,基于表2模擬各路段過網(wǎng)電價,列于表3。同時模擬政府性基金電價為14元/MWh;G33上網(wǎng)電價為350元/(MW·h);G35和G36的上網(wǎng)電價分別為360元/(MW·h)和 370元/(MW·h)。
表2 與{G33, L16}關(guān)聯(lián)的電能配送路徑鏈Table 2 Energy distribution chains for {G33, L16}
表3中,33-19、35-22、36-23為電廠主變,不納入過網(wǎng)費計算范疇,且相關(guān)發(fā)電商的上網(wǎng)費以網(wǎng)供電能(電廠高壓母線側(cè))為準核算,即對相應(yīng)路徑鏈,應(yīng)取次路段首端(即首路段末端)的能量。
表3 各路段過網(wǎng)電價模擬值Table 3 Simulated transmission price of branch sections
4.2.1 源-流對
由表2知,
表4 1號鏈能流信息Table 4 Energy flow information of the 1st chain
由表5可見,當核算上網(wǎng)電費時,若采用負荷側(cè)的電能計量值 188.432133 MW·h,則對應(yīng)的結(jié)算價格須采用考慮電網(wǎng)實時運行狀態(tài)因素折算后的價格352.782元/(MW·h)。在此條件下,考慮價格聯(lián)動機制和成本平衡回收原理,折算出的下網(wǎng)電價分別為 463.890元/(MW·h)(mode I)和467.577元/(MW·h)(mode II);負荷側(cè)核算的政府性基金價格為 14.000元/(MW·h)(mode II),折算到 mode I的價格為 13.890元/(MW·h)。類似地,可得子鏈及其各路段在mode I和II下的過網(wǎng)價。表5最下面兩行表明,無論mode I還是II,各項費用均對應(yīng)相等,滿足核算一致性原理。
表5 1號鏈經(jīng)濟流信息Table 5 Economic flow information of the 1st chain
相比較的,若直接按中發(fā) 〔2015〕 9號文規(guī)定,則下網(wǎng)電價 = 350 + 14 + 100 = 464元/(MW·h),由此核算的下網(wǎng)費 = 464 × 188.432133 = 87432.509712元,相對于表 5中的 88106.466 元,少了 673.956288元(約7.649‰),說明其不完全滿足核算一致性原理。
4.2.2 其他關(guān)聯(lián)源-流對間經(jīng)濟流計算結(jié)果
類似于4.2.1節(jié),可計算出G33與OD和L16與OG間的經(jīng)濟流分布。限于篇幅,下面只重點以核算L16的下網(wǎng)費為例進行分析。
L16除了通過網(wǎng)絡(luò)源-流對
表6 8 ~10號鏈的能流信息Table 6 Energy flow information of the 8th to 10th chain
由表6能流信息和4.1節(jié)提供的價格信息,可以計算出負荷L16在8 ~10號鏈末端的下網(wǎng)電價(以 mode II為例)分別為:567.587元/(MW·h)、557.915元/(MW·h)和 507.103元/(MW·h);再乘以下網(wǎng)電能,可得對應(yīng)的下網(wǎng)費分別為:69412.389、568.306、7071.036元。與通過1號鏈購電的下網(wǎng)費合計為88106.466元,可得L16從發(fā)電商G33、G35和G36購電的總費用為165158.197元。除以L16從3個發(fā)電商所購總電能(負荷側(cè)計量)325.688573 MW·h,可得 L16的平均下網(wǎng)電價為507.1047元/(MW·h)。
4.2.3 結(jié)果討論
由基本理論和算例結(jié)果可知,只要能量管理系統(tǒng)提供相應(yīng)潮流狀態(tài)(實時潮流狀態(tài)或目標調(diào)控后的最優(yōu)潮流狀態(tài)),市場平臺提供有關(guān)對象的合同簽約價格,則均可精細化分析出電網(wǎng)不同網(wǎng)絡(luò)源-流對間或不同市場交易對象間的所有相關(guān)能量流和經(jīng)濟流分布結(jié)果。這些結(jié)果,可應(yīng)用于指導(dǎo)發(fā)電商/購電商優(yōu)化其買賣行為或指導(dǎo)電網(wǎng)公司優(yōu)化其運行方式。如可將4.2.2節(jié)計算結(jié)果應(yīng)用于指導(dǎo)如何進一步優(yōu)化或降低L16的平均單位購電成本。因L16從G33購電價格最低,從G35購電價格最高,若L16適當降低從G35購電的量、同時提高從G33購電量,則L16的平均實時下網(wǎng)電價必然會有所降低,從而節(jié)省購電費用。此外能量流和經(jīng)濟流耦合分布信息計算結(jié)果,可進一步應(yīng)用于分析成本分攤和交叉補貼問題[22-25]及更為復(fù)雜的技術(shù)-經(jīng)濟協(xié)同下的多目標優(yōu)化問題。
本文通過對電力網(wǎng)絡(luò)源流間的能量流與經(jīng)濟流的耦合分析,得出結(jié)論:(1)不同網(wǎng)絡(luò)源流間的能量流分布決定于電網(wǎng)平臺提供的系統(tǒng)運行信息,不同商業(yè)伙伴間的經(jīng)濟流分布決定于市場平臺提供的電能交易信息,二者因依托于共同的電網(wǎng)而產(chǎn)生耦合關(guān)系;(2)考慮能量流-經(jīng)濟流耦合特性,基于上網(wǎng)能量計費方式和下網(wǎng)能量計費方式分別核算各項費用的電價成分,因電網(wǎng)運行狀態(tài)的時變性也具有了時變特征,且兩種計費方式滿足核算一致性原理;(3)本文方法可進一步應(yīng)用于交叉補貼機理揭示、非政策性補貼費用的清退與返還規(guī)則設(shè)計,以及多個技術(shù)-經(jīng)濟目標的協(xié)同優(yōu)化等問題。