考慮網(wǎng)絡(luò)損耗的碳排放流理論改進(jìn)與完善

馮 欣 ，楊 軍

（1.武漢大學(xué) 電氣工程學(xué)院，湖北 武漢 430072；2.國(guó)網(wǎng)蘇州供電公司，江蘇 蘇州 215004）

0 引言

隨著世界工業(yè)經(jīng)濟(jì)的發(fā)展、人口的劇增、人類(lèi)生產(chǎn)生活方式的愈加無(wú)節(jié)制，世界氣候面臨越來(lái)越嚴(yán)重的問(wèn)題，二氧化碳排放量越來(lái)越大，地球臭氧層正遭受前所未有的危機(jī)，人類(lèi)的生存環(huán)境和健康安全已經(jīng)受到嚴(yán)重危害。為了應(yīng)對(duì)當(dāng)今氣候變化與能源危機(jī)，減輕溫室效應(yīng)對(duì)全球的影響，低碳發(fā)展是社會(huì)可持續(xù)發(fā)展的必經(jīng)之路［1-6］。低碳之路，電力先行。電力系統(tǒng)作為碳排放的主要來(lái)源之一，必將在低碳電力領(lǐng)域做出重大貢獻(xiàn)。文獻(xiàn)［7-8］從國(guó)家層面進(jìn)行碳排放估算，認(rèn)為發(fā)展低碳電力必須對(duì)電力系統(tǒng)碳排放有深刻的認(rèn)識(shí)，分析其碳排放結(jié)構(gòu)，了解構(gòu)成電力碳排放的不同原因［9-13］，分析其影響電力系統(tǒng)相關(guān)排放的機(jī)理。

在當(dāng)前電力系統(tǒng)低碳減排的發(fā)展趨勢(shì)下［14］，電力系統(tǒng)碳排放流理論的引入為電力系統(tǒng)低碳發(fā)展的研究帶來(lái)了新的思路和方法［15-17］。 文獻(xiàn)［15］提出了將碳排放視為依附于潮流而存在的虛擬網(wǎng)絡(luò)流的思想，將碳排放與電力潮流分析相結(jié)合，建立了電力系統(tǒng)碳排放流的理念和理論框架。文獻(xiàn)［16］將潮流計(jì)算與電力系統(tǒng)碳排放流計(jì)算相結(jié)合，進(jìn)一步提出了碳排放流計(jì)算的基本方法。文獻(xiàn)［17］分析了碳排放流從發(fā)電側(cè)流向用電側(cè)的具體過(guò)程，研究了碳排放流在電網(wǎng)中的分布特性和機(jī)理。碳排放流理論的提出為電力系統(tǒng)碳排放量的計(jì)算與分析提供了新的思路和手段，具有良好發(fā)展前景。

但是，文獻(xiàn)［15-17］中提出的碳排放流理論只能應(yīng)用在忽略網(wǎng)損的情況下，無(wú)法應(yīng)用于實(shí)際有損電力系統(tǒng)，需要進(jìn)一步完善。本文考慮實(shí)際電力系統(tǒng)中的網(wǎng)絡(luò)損耗，通過(guò)計(jì)算負(fù)荷節(jié)點(diǎn)等效負(fù)荷需求和線(xiàn)路等效傳輸功率，將網(wǎng)損分?jǐn)偟礁髫?fù)荷，從而將有損網(wǎng)絡(luò)轉(zhuǎn)換成無(wú)損網(wǎng)絡(luò)，對(duì)碳排放流理論進(jìn)行了補(bǔ)充和完善，使其能夠適用于實(shí)際電力系統(tǒng)。所提方法能夠追蹤電力系統(tǒng)中碳排放流的具體流向，回溯負(fù)荷側(cè)碳流來(lái)源，為碳交易市場(chǎng)下負(fù)荷側(cè)用戶(hù)應(yīng)當(dāng)繳納的碳配額計(jì)算［18-19］提供參考依據(jù)，有利于促進(jìn)用戶(hù)參與節(jié)能減排行動(dòng)，推動(dòng)低碳電力的發(fā)展。

1 未考慮網(wǎng)損的碳排放流理論基礎(chǔ)

碳排放流基本理論［15-17］介紹了碳排放流在電力系統(tǒng)中的分布機(jī)理與特性，提出了用機(jī)組-節(jié)點(diǎn)碳流關(guān)聯(lián)矩陣RU-N、機(jī)組-支路碳流關(guān)聯(lián)矩陣RU-B和機(jī)組-負(fù)荷碳流關(guān)聯(lián)矩陣RU-L3個(gè)關(guān)聯(lián)矩陣來(lái)描述電力系統(tǒng)中的碳排放流分布。其中，機(jī)組-節(jié)點(diǎn)碳流關(guān)聯(lián)矩陣表示系統(tǒng)中所有發(fā)電機(jī)組的碳排放流注入對(duì)流入系統(tǒng)中任一節(jié)點(diǎn)碳流率的貢獻(xiàn)情況，節(jié)點(diǎn)碳流率的物理意義為發(fā)電側(cè)為維持節(jié)點(diǎn)有功功率每單位時(shí)間產(chǎn)生的碳排放量，RU-N可由下式計(jì)算：

其中，EG為發(fā)電機(jī)組碳排放強(qiáng)度向量；PN為節(jié)點(diǎn)有功通量矩陣；PB為支路潮流分布矩陣；PG為機(jī)組注入分布矩陣。

機(jī)組-支路碳流關(guān)聯(lián)矩陣RU-B，k表示系統(tǒng)中第k臺(tái)發(fā)電機(jī)組的碳排放流注入對(duì)系統(tǒng)中所有支路碳流率的貢獻(xiàn)情況，支路碳流率的物理意義為單位時(shí)間內(nèi)系統(tǒng)為維持有功潮流而在發(fā)電廠產(chǎn)生的碳排放量，RU-B，k可由下式計(jì)算：

其中，為N維行向量，其第k個(gè)元素為1，其余元素為0；eGk為第k臺(tái)發(fā)電機(jī)的碳排放強(qiáng)度。

機(jī)組－負(fù)荷碳流關(guān)聯(lián)矩陣RU-L表示系統(tǒng)中的各負(fù)荷節(jié)點(diǎn)碳流率來(lái)自于系統(tǒng)中各機(jī)組的分量，負(fù)荷節(jié)點(diǎn)碳流率的物理意義為發(fā)電側(cè)為供應(yīng)節(jié)點(diǎn)負(fù)荷每單位時(shí)間產(chǎn)生的碳排放量，RU-L可由下式計(jì)算：

其中，ζN為N維行向量，其所有元素都為1；PL為負(fù)荷分布矩陣。

上述關(guān)聯(lián)矩陣的計(jì)算前提是無(wú)損網(wǎng)絡(luò)，但實(shí)際電力系統(tǒng)網(wǎng)絡(luò)是有損的，因此該方法并不適用于實(shí)際電力系統(tǒng)。為使碳排放流理論適用于實(shí)際電力系統(tǒng)，本文對(duì)其進(jìn)行了改進(jìn)，將實(shí)際有損網(wǎng)絡(luò)轉(zhuǎn)換成無(wú)損網(wǎng)絡(luò)［20］，并對(duì)上述3個(gè)關(guān)聯(lián)矩陣進(jìn)行了重新計(jì)算。

2 無(wú)損網(wǎng)絡(luò)轉(zhuǎn)換

當(dāng)系統(tǒng)為有損網(wǎng)絡(luò)時(shí)，則需要考慮到網(wǎng)損，可通過(guò)某種轉(zhuǎn)換將有損網(wǎng)絡(luò)等效為無(wú)損網(wǎng)絡(luò)。隨著碳交易市場(chǎng)的發(fā)展，文獻(xiàn)［18］提出消費(fèi)者購(gòu)買(mǎi)電力產(chǎn)品時(shí)除了應(yīng)當(dāng)支付相應(yīng)的現(xiàn)金，還需要提交與這些能源商品、服務(wù)的碳排放相當(dāng)?shù)奶寂漕~。電力系統(tǒng)碳排放的本質(zhì)來(lái)源是負(fù)荷側(cè)的消費(fèi)需要，需求導(dǎo)致生產(chǎn)，所以電力負(fù)荷才是碳排放的真正源頭［15］。因此，電力系統(tǒng)網(wǎng)絡(luò)損耗的產(chǎn)生也是由消費(fèi)者需求產(chǎn)生的?；诖?，在碳交易市場(chǎng)背景下，本文將網(wǎng)損等效到負(fù)荷側(cè)，把由網(wǎng)損引起的碳排放分?jǐn)偨o負(fù)荷側(cè)消費(fèi)者，為碳交易市場(chǎng)下負(fù)荷側(cè)用戶(hù)繳納的碳配額計(jì)算提供相應(yīng)的參考依據(jù)，促使負(fù)荷側(cè)用戶(hù)參與節(jié)能行動(dòng)，推動(dòng)低碳電力的發(fā)展。

令為等效無(wú)損網(wǎng)絡(luò)下節(jié)點(diǎn)i的有功潮流，為等效無(wú)損網(wǎng)絡(luò)下支路i-j流過(guò)的有功功率，則根據(jù)文獻(xiàn)［20］有：

其中，為有功潮流直接流入節(jié)點(diǎn)i的節(jié)點(diǎn)集合；PGi為接入節(jié)點(diǎn)i的發(fā)電機(jī)組出力。由于那么可以表示為其中文獻(xiàn)［20］中提出，由于線(xiàn)路有功損耗相對(duì)較小，則可假設(shè)其中，Pj-i為實(shí)際有損網(wǎng)絡(luò)中支路j-i從節(jié)點(diǎn)j流向節(jié)點(diǎn)i的有功潮流，Pj為實(shí)際有損網(wǎng)絡(luò)中節(jié)點(diǎn) j流過(guò)的有功潮流，則式（4）可以寫(xiě)成：

亦可寫(xiě)成：

其中，Pg表示等效為無(wú)損網(wǎng)絡(luò)后電網(wǎng)各節(jié)點(diǎn)的有功潮流向量；Au為n×n階矩陣，可表示為式（7）。

由于Au和PG已知，則可計(jì)算出Pg。由Pg可計(jì)算出等效為無(wú)損網(wǎng)絡(luò)后網(wǎng)絡(luò)中支路有功潮流和各節(jié)點(diǎn)的有功負(fù)荷。線(xiàn)路i-l的等效有功潮流為：

其中，為從節(jié)點(diǎn)i直接流入有功潮流的節(jié)點(diǎn)集合。同樣，可求出節(jié)點(diǎn)i的等效有功負(fù)荷為：

其中，包含了節(jié)點(diǎn)i的實(shí)際負(fù)荷以及該負(fù)荷分?jǐn)偟木W(wǎng)損。

因此，聯(lián)立式（7）—（9），可將有損網(wǎng)絡(luò)等效為無(wú)損網(wǎng)絡(luò)，使得網(wǎng)損被分?jǐn)偟礁髫?fù)荷處，再將其應(yīng)用于碳排放流理論的計(jì)算。

3 仿真分析與校驗(yàn)

3.1 網(wǎng)損轉(zhuǎn)移計(jì)算

為了更清楚地展示有損網(wǎng)絡(luò)轉(zhuǎn)換成無(wú)損網(wǎng)絡(luò)的過(guò)程，及網(wǎng)損如何被分?jǐn)偟礁髫?fù)荷處，現(xiàn)通過(guò)2機(jī)4節(jié)點(diǎn)系統(tǒng)和IEEE 30節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)來(lái)進(jìn)行說(shuō)明。

3.1.1 2機(jī)4節(jié)點(diǎn)系統(tǒng)

圖1所示為一個(gè)4節(jié)點(diǎn)交流系統(tǒng)［20］（圖中功率流的單位為MW），節(jié)點(diǎn)1和節(jié)點(diǎn)2處分別裝有發(fā)電機(jī)組，發(fā)電機(jī)有功出力分別為PG1=400 MW和PG2=114 MW；節(jié)點(diǎn)3和節(jié)點(diǎn)4為負(fù)荷節(jié)點(diǎn)，負(fù)荷分別為PL3=300 MW和PL4=200 MW。計(jì)算可知，該網(wǎng)絡(luò)系統(tǒng)總網(wǎng)損為（225-218）+（83-82）+（173-171）+（60-59）+（115-112）=14（MW）。

圖1 4節(jié)點(diǎn)系統(tǒng)穩(wěn)態(tài)潮流分布Fig.1 Steady power-flow of 4-bus system

則根據(jù)式（7）可以求出矩陣Au為：

將 Au代入式（6）中有：

求解上式可得到：

當(dāng)Pg被求解出來(lái)后，則可以根據(jù)式（9）求出網(wǎng)絡(luò)等效情況下負(fù)荷節(jié)點(diǎn)3和負(fù)荷節(jié)點(diǎn)4處的等效負(fù)荷分別為：

網(wǎng)絡(luò)未進(jìn)行等效時(shí)，負(fù)荷節(jié)點(diǎn)3和負(fù)荷節(jié)點(diǎn)4處的實(shí)際負(fù)荷分別為PL3=300 MW和PL4=200 MW；網(wǎng)絡(luò)等效后負(fù)荷節(jié)點(diǎn)3和負(fù)荷節(jié)點(diǎn)4所分?jǐn)偩W(wǎng)損分別為309.76-300=9.76（MW）和 204.24-200=4.24（MW），2個(gè)負(fù)荷節(jié)點(diǎn)所分?jǐn)偩W(wǎng)損之和為9.76+4.24=14（MW），等于系統(tǒng)的總網(wǎng)損。這表明經(jīng)過(guò)文中的潮流追蹤算法計(jì)算后網(wǎng)絡(luò)中的網(wǎng)損被分?jǐn)偟搅烁髫?fù)荷處。

同時(shí)，根據(jù)式（9）可以追蹤到各負(fù)荷來(lái)自于各臺(tái)發(fā)電機(jī)組的功率如表1所示。

由表1可發(fā)現(xiàn)，網(wǎng)絡(luò)等效成無(wú)損網(wǎng)絡(luò)后，發(fā)電機(jī)發(fā)出的有功功率全部都由負(fù)荷消納。負(fù)荷節(jié)點(diǎn)3來(lái)自于發(fā)電機(jī) G1和G2的功率分別為276.32 MW和33.44 MW，其分?jǐn)偟木W(wǎng)損為9.76 MW。負(fù)荷節(jié)點(diǎn)4來(lái)自于發(fā)電機(jī)G1和G2的功率分別為123.68 MW和80.56 MW，其分?jǐn)偟木W(wǎng)損為4.24 MW。所有負(fù)荷分?jǐn)偟目偩W(wǎng)損為14 MW，等于系統(tǒng)的總網(wǎng)損。

表1 4節(jié)點(diǎn)系統(tǒng)電力分布Table1 Power distribution of 4-bus system

3.1.2 IEEE 30節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)

本文進(jìn)一步利用IEEE 30節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)來(lái)說(shuō)明網(wǎng)損是如何轉(zhuǎn)移到負(fù)荷側(cè)的。IEEE 30節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)中各機(jī)組有功出力、各用戶(hù)負(fù)荷及系統(tǒng)穩(wěn)態(tài)有功潮流分布如圖2所示，圖中功率流單位為MW。

可根據(jù)圖2中各線(xiàn)路首末兩端的功率差求得系統(tǒng)的總網(wǎng)損為2.36 MW。同理根據(jù)式（4）—（9），可求得系統(tǒng)各節(jié)點(diǎn)的等效負(fù)荷。系統(tǒng)各節(jié)點(diǎn)的實(shí)際負(fù)荷和等效負(fù)荷如表2所示，表2中同時(shí)列出了各個(gè)節(jié)點(diǎn)分?jǐn)偟木W(wǎng)損，且所有負(fù)荷節(jié)點(diǎn)分?jǐn)偟木W(wǎng)損之和為2.36 MW，等于系統(tǒng)的總網(wǎng)損。這表明網(wǎng)絡(luò)等效后系統(tǒng)中的網(wǎng)損都被分?jǐn)偟搅烁髫?fù)荷處。同時(shí)追蹤網(wǎng)絡(luò)等效后各負(fù)荷來(lái)自于各臺(tái)機(jī)組的功率，可發(fā)現(xiàn)發(fā)電機(jī)發(fā)出的有功功率全部都由負(fù)荷消納。

3.2 等效網(wǎng)絡(luò)下碳流理論的計(jì)算

為了展示考慮網(wǎng)絡(luò)損耗的碳排放流理論在實(shí)際電力系統(tǒng)中的應(yīng)用，本文采用如圖2所示IEEE 30節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)進(jìn)行示例分析。

圖2 IEEE 30節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)穩(wěn)態(tài)潮流分布Fig.2 Steady power-flow of IEEE 30-bus standard test system

在實(shí)際有損網(wǎng)絡(luò)和忽略網(wǎng)損2種情況下，本文應(yīng)用碳排放理論分別計(jì)算碳排放流各相關(guān)矩陣。設(shè)定機(jī)組G1、G2和G6為燃煤機(jī)組，其碳排放強(qiáng)度為0.875 kg/（kW·h）（即生產(chǎn)每 kW·h 電量排放 CO20.875 kg）；機(jī)組G3和機(jī)組G5為燃?xì)鈾C(jī)組，其碳排放強(qiáng)度為 0.525 kg /（kW·h）［16］；機(jī)組 G4為水輪機(jī)，其碳排放強(qiáng)度為0。

表2 各節(jié)點(diǎn)處負(fù)荷及所分?jǐn)偟木W(wǎng)損Table 2 Load and allocated power-loss of each node MW

可分別計(jì)算在無(wú)損網(wǎng)絡(luò)和有損網(wǎng)絡(luò)中各個(gè)節(jié)點(diǎn)處碳流率的注入量來(lái)自于系統(tǒng)中各個(gè)發(fā)電機(jī)的分量。各節(jié)點(diǎn)處的總碳流率如圖3所示，可發(fā)現(xiàn)除了在節(jié)點(diǎn)5、9、11、13處外，其余節(jié)點(diǎn)的總碳流率在考慮網(wǎng)損情況下都要比忽略網(wǎng)損情況下大。這是因?yàn)樵谟袚p情況下，網(wǎng)絡(luò)中線(xiàn)路損耗要消耗一部分碳排放流，但是由于節(jié)點(diǎn)5、11和13分別和發(fā)電機(jī)組G3、G5、G6在同一條母線(xiàn)上，有損和無(wú)損情況下各發(fā)電機(jī)組向相應(yīng)節(jié)點(diǎn)處的注入碳流率一定；節(jié)點(diǎn)9和節(jié)點(diǎn)11之間的線(xiàn)路電阻為零，相當(dāng)于節(jié)點(diǎn)9和節(jié)點(diǎn)11處于同一母線(xiàn)，因此并沒(méi)有線(xiàn)路損耗，所以在有損和無(wú)損2種情況下，節(jié)點(diǎn)5、9、11、13的總注入碳流率相等。另外，由于節(jié)點(diǎn)8處安裝的是水輪機(jī)組，因此其注入的總碳流率為零。

計(jì)算無(wú)損網(wǎng)絡(luò)和有損網(wǎng)絡(luò)2種情況下各臺(tái)發(fā)電機(jī)組對(duì)各條支路的碳流率貢獻(xiàn)可發(fā)現(xiàn)，有損情況下大多數(shù)支路的總碳流率比無(wú)損情況下大，這是因?yàn)橛袚p情況下考慮了線(xiàn)路損耗。只有小部分支路是相等或者較小，這部分支路的支路碳流率如圖4所示。從圖4中可發(fā)現(xiàn)，支路11-9和13-12處有損和無(wú)損情況下支路碳流率相等，這是因?yàn)檫@2條支路之間的電阻為零，支路損耗為零，節(jié)點(diǎn)11和節(jié)點(diǎn)9、節(jié)點(diǎn)13和節(jié)點(diǎn)12相當(dāng)于在同一條母線(xiàn)上，因此其支路碳流率相等。在支路2-5、支路5-7、支路9-6、支路10-17、支路27-25處，有損網(wǎng)絡(luò)下的支路碳流率要比無(wú)損情況下小，是由于其他支路的碳流率增大導(dǎo)致這5條支路的碳流率減小。將該結(jié)果與潮流結(jié)果對(duì)比發(fā)現(xiàn)，與潮流結(jié)果相符。根據(jù)潮流結(jié)果發(fā)現(xiàn)，在有損情況下，上述5條支路流過(guò)的潮流也比無(wú)損情況下小，原因主要是其他支路增加的潮流補(bǔ)償了這部分支路減少的潮流。

圖3 各節(jié)點(diǎn)流入總碳流率Fig.3 Gross carbon emission flow rate of each node

圖4 部分支路流過(guò)的總碳流率Fig.4 Gross carbon emission flow rate of some branches

計(jì)算無(wú)損網(wǎng)絡(luò)和有損網(wǎng)絡(luò)2種情況下各發(fā)電機(jī)組對(duì)各支路的碳流率貢獻(xiàn)可發(fā)現(xiàn)，2種情況下各機(jī)組對(duì)支路8-6和支路8-28注入的碳流率均為零。這是因?yàn)橹?-6和支路8-28連接的是水輪機(jī)G4，支路功率完全由G4供應(yīng)，而G4并沒(méi)有碳排放，所以其碳流率都為零，因此這2條支路中的線(xiàn)損所引起的碳流率也為零。

因此，當(dāng)線(xiàn)路連接的是水輪機(jī)或風(fēng)機(jī)等沒(méi)有碳排放的機(jī)組時(shí)，該線(xiàn)路的功率由相應(yīng)機(jī)組提供，則該支路流過(guò)的碳流率為零，支路線(xiàn)損所引起的碳流率也為零。

圖5 各負(fù)荷處總碳流率Fig.5 Gross carbon emission flow rate of each load

計(jì)算有損和無(wú)損網(wǎng)絡(luò)2種情況下的負(fù)荷節(jié)點(diǎn)碳流率，相應(yīng)各節(jié)點(diǎn)處的負(fù)荷總碳流率如圖5所示。從圖5中可發(fā)現(xiàn)，在有損情況下，各節(jié)點(diǎn)負(fù)荷的總碳流率比無(wú)損情況下大。這是因?yàn)樵谟袚p網(wǎng)絡(luò)中，線(xiàn)路損耗轉(zhuǎn)移到了負(fù)荷側(cè)，使得各負(fù)荷節(jié)點(diǎn)等效的負(fù)荷需求比實(shí)際大，其所需的碳流率也相應(yīng)地增大，線(xiàn)路損耗成功地分?jǐn)偟礁髫?fù)荷處。計(jì)算考慮網(wǎng)損和忽略網(wǎng)損2種情況下的各節(jié)點(diǎn)負(fù)荷處總碳流率的相對(duì)偏差率，可發(fā)現(xiàn)節(jié)點(diǎn)26處相對(duì)偏差率最大，為6.83%。

本文對(duì)IEEE 30節(jié)點(diǎn)系統(tǒng)進(jìn)行仿真，考慮到實(shí)際大電網(wǎng)中復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)和長(zhǎng)線(xiàn)路實(shí)際情況，則考慮網(wǎng)損比忽略網(wǎng)損時(shí)的負(fù)荷節(jié)點(diǎn)碳流率大?？梢?jiàn)，實(shí)際電網(wǎng)中網(wǎng)絡(luò)損耗引起的碳排放不能被忽略，因此本文所提改進(jìn)方法實(shí)用價(jià)值較好，能為碳交易市場(chǎng)下負(fù)荷側(cè)應(yīng)當(dāng)交納的碳配額計(jì)算提供參考依據(jù)。各機(jī)組注入電網(wǎng)的總碳流率如圖6所示。由圖6可見(jiàn)，有損和無(wú)損2種情況下除了平衡機(jī)組G1注入電網(wǎng)的碳流率不同之外，其余機(jī)組發(fā)出的碳流率均相等，符合潮流計(jì)算結(jié)果，說(shuō)明有損網(wǎng)絡(luò)和無(wú)損網(wǎng)絡(luò)都不影響除了平衡節(jié)點(diǎn)外的發(fā)電機(jī)組出力。

圖6 各機(jī)組注入電網(wǎng)的碳流率Fig.6 Rate of carbon emission flow injected into power grid by each unit

4 結(jié)論

本文考慮到實(shí)際電力系統(tǒng)中的網(wǎng)絡(luò)損耗，通過(guò)計(jì)算線(xiàn)路等效傳輸功率和負(fù)荷節(jié)點(diǎn)等效負(fù)荷需求，將網(wǎng)損分?jǐn)偟礁髫?fù)荷，從而將有損網(wǎng)絡(luò)轉(zhuǎn)換成無(wú)損網(wǎng)絡(luò)，對(duì)碳排放流理論進(jìn)行了完善，使其能夠適用于實(shí)際有損電力系統(tǒng)。IEEE 30節(jié)點(diǎn)系統(tǒng)仿真結(jié)果表明，忽略網(wǎng)損時(shí)的碳流相關(guān)計(jì)算結(jié)果相比于實(shí)際有損電網(wǎng)存在一定偏差，而本文提出的考慮網(wǎng)絡(luò)損耗的碳排放流理論能夠清晰描述電力系統(tǒng)中碳排放流的流向和分布，可以回溯負(fù)荷側(cè)的碳流來(lái)源，為在碳交易市場(chǎng)背景下負(fù)荷側(cè)用戶(hù)應(yīng)當(dāng)繳納的碳配額計(jì)算提供參考依據(jù)，有利于促使用戶(hù)參與節(jié)能減排行動(dòng)，推動(dòng)低碳電力的發(fā)展。

［1］康重慶，陳啟鑫，夏清．低碳電力技術(shù)的研究展望［J］．電網(wǎng)技術(shù)，2009，33（2）：1-7．KANG Chongqing，CHEN Qixin，XIA Qing．Prospectsoflow carbon eletricity［J］．Power System Technology，2009，33（2）：1-7．

［2］曾鳴，薛松，魏陽(yáng)，等．低碳智慧島能源體系設(shè)計(jì)及經(jīng)濟(jì)可靠性評(píng)價(jià)模型［J］．電力自動(dòng)化設(shè)備，2012，32（12）：58-63．ZENG Ming，XUE Song，WEI Yang，etal．Energysystem for low-carbon-smart islands and economy-reliability assessment model［J］．Electric Power Automation Equipment，2012，32（12）：58-63．

［3］辛建波，范瑞祥，鄭蜀江，等．智能配用電園區(qū)低碳效益仿真與評(píng)價(jià)系統(tǒng)及其在江西的應(yīng)用［J］．電力系統(tǒng)保護(hù)與控制，2014，42（7）：86-90．XIN Jianbo，F(xiàn)AN Ruixiang，ZHENG Shujiang，et al．Low-carbon benefit simulation and evaluation system for smart distribution grid and its application in Jiangxi Province［J］．Power System Protection and Control，2014，42（7）：86-90．

［4］朱澤磊，周京陽(yáng)，韓忠旭，等．基于多階段不同政策機(jī)制的低碳電力調(diào)度模型［J］．電網(wǎng)技術(shù)，2013，37（2）：287-294．ZHU Zelei，ZHOU Jingyang，HAN Zhongxu，et al．A low-carbon power dispatch model based on policy mechanisms corresponding to different stages［J］．Power System Technology，2013，37（2）：287-294．

［5］文旭，王俊梅，郭琳，等．計(jì)及污染氣體排放風(fēng)險(xiǎn)的多目標(biāo)隨機(jī)動(dòng)態(tài)環(huán)境經(jīng)濟(jì)調(diào)度模型［J］．電力自動(dòng)化設(shè)備，2015，35（5）：131-138．WEN Xu，WANG Junmei，GUO Lin，et al． Multi-objective stochastic and dynamic model of environmental and economic dispatch considering gas pollution emission risk［J］．Electric Power Automation Equipment，2015，35（5）：131-138．

［6］龍軍，鄭斌，郭小璇，等．一種求解環(huán)境經(jīng)濟(jì)發(fā)電調(diào)度的交互式多目標(biāo)優(yōu)化方法［J］．電力自動(dòng)化設(shè)備，2013，33（5）：83-88．LONG Jun，ZHENG Bin，GUO Xiaoxuan，et al．Interactive multiobjective optimization ofenvironmentaland economic power generation scheduling［J］．Electric Power Automation Equipment，2013，33（5）：83-88．

［7］張斌．2020年我國(guó)能源電力消費(fèi)及碳排放強(qiáng)度情景分析［J］．中國(guó)能源，2009（3）：28-31．ZHANG Bin．Scenario analyses of China’s energy and power use and carbon emission in 2020［J］．Energy of China，2009（3）：28-31．

［8］CHUNG W，TOHNO S，SHIM S Y．An estimation of energy and GHG emission intensity caused by energy consumption in Korea：an energy IO approach［J］．Applied Energy，2009，86（10）：1902-1914．

［9］IPCC．Intergovernmental panel for climate change：the fourth assessment report［M］．Cambridge，UK：Cambridge University Press，2007．

［10］GRUBB M，JAMASB T，POLLITT M G.Delivering a low carbon electricity system［M］．Cambridge，UK：Cambridge University Press，2008．

［11］JAMASB T，NUTTALL W J，POLLITT M G.Future electricity technologies and systems［M］．Cambridge，UK：Cambridge University Press，2008．

［12］張坤民，潘家華，崔大鵬．低碳經(jīng)濟(jì)論［M］．北京：中國(guó)環(huán)境科學(xué)出版社，2008．

［13］魏一鳴，劉蘭翠，范英，等．中國(guó)能源報(bào)告（2008）：碳排放研究［M］．北京：科學(xué)出版社，2008．

［14］陳曉科，周天睿，李欣，等．電力系統(tǒng)的碳排放結(jié)構(gòu)分解與低碳目標(biāo)貢獻(xiàn)分析［J］．電力系統(tǒng)自動(dòng)化，2012，36（2）：18-25．CHEN Xiaoke，ZHOU Tianrui，LI Xin，et al． Structure identification of CO2emission for power system and analysis of its low-carbon contribution［J］．Automation of Electric Power System，2012，36（2）：18-25．

［15］周天睿，康重慶，徐乾耀，等．電力系統(tǒng)碳排放流分析理論初探［J］．電力系統(tǒng)自動(dòng)化，2012，36（7）：38-43．ZHOU Tianrui，KANG Chongqing，XU Qianyao，et al．Preliminary theoretical investigation on power system carbon emission flow［J］．Automation of Electric Power Systems，2012，36（7）：38-43．

［16］周天睿，康重慶，徐乾耀，等．電力系統(tǒng)碳排放流的基本計(jì)算方法初探［J］．電力系統(tǒng)自動(dòng)化，2012，36（11）：44-49．ZHOU Tianrui，KANG Chongqing，XU Qianyao，et al．Preliminary investigation on a method for carbon emission flow calculation of power system［J］．Automation of Electric Power Systems，2012，36（11）：44-49．

［17］周天睿，康重慶，徐乾耀，等．碳排放流在電力網(wǎng)絡(luò)中分布的特性與機(jī)理分析［J］．電力系統(tǒng)自動(dòng)化，2012，36（15）：39-44．ZHOU Tianrui，KANG Chongqing，XU Qianyao，et al．Analysis on distribution characteristics and mechanisms of carbon emission flow in electric power network［J］．Automation of Electric Power Systems，2012，36（15）：39-44．

［18］范進(jìn).基于個(gè)人碳交易行為模型的電力消費(fèi)選擇研究［D］.合肥：中國(guó)科學(xué)技術(shù)大學(xué)，2012.FAN Jin.The electricity consumption choice based on the model of personal carbon trading behavior［D］.Hefei：University of Science and Technology of China，2012.

［19］宋旭東，莫娟，向鐵元．電力行業(yè)碳排放權(quán)的初始分配機(jī)制［J］．電力自動(dòng)化設(shè)備，2013，33（1）：44-49．SONG Xudong，MO Juan，XIANG Tieyuan．Initialallocation mechanism of carbon emission permit in electric power industry［J］．Electric Power Automation Equipment，2013，33（1）：44-49．

［20］BIALEK J．Tracing the flow of electricity［J］．Generation，Transmission and Distribution，IEE Proceedings，1996，143 （4）：313-320．