上海交通大学学报, 2023, 57(4): 422-431 doi: 10.16183/j.cnki.jsjtu.2021.513

新型电力系统与综合能源

考虑经济性与碳排放的电-气综合能源系统多目标规划

朱海南,1, 王娟娟1, 陈兵兵1, 张厚望2, 陈健2, 吴秋伟3

1.国网山东省电力公司潍坊供电公司, 山东 潍坊 261000

2.山东大学 电气工程学院, 济南 250061

3.丹麦科技大学 电气工程学院,丹麦 灵比 2800

Multi-Objective Planning of Power-Gas Integrated Energy System Considering Economy and Carbon Emission

ZHU Hainan,1, WANG Juanjuan1, CHEN Bingbing1, ZHANG Houwang2, CHEN Jian2, WU Qiuwei3

1. State Grid Weifang Power Supply Company, Weifang 261000, Shandong, China

2. School of Electrical Engineering, Shandong University, Jinan 250061, China

3. Department of Electrical Engineering, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark

责任编辑: 王一凡

收稿日期: 2021-12-16   修回日期: 2022-01-6   接受日期: 2022-02-7  

基金资助: 国网山东省电力公司科技项目(520604200003)
国家重点研发计划(2018YFA0702200)

Received: 2021-12-16   Revised: 2022-01-6   Accepted: 2022-02-7  

作者简介 About authors

朱海南(1987-),高级工程师,博士,从事电网运行与控制研究;E-mail:hainanzhu@mail.sdu.edu.cn.

摘要

为加速电-气系统快速、经济的低碳转型,构建了一种综合考虑经济性成本与碳排放量的电-气综合能源系统多目标随机优化规划模型.首先建立电-气网络与相关设备的数学模型,并运用场景法表征电、气负荷与光伏出力的不确定性.其次建立综合考虑系统经济性成本和碳排放量两个指标的混合整数二次约束规划(MIQCP)模型,对电网馈线、气网管道、变电站、配气站、燃气机组、电转气装置、光伏及储能装置进行统筹规划.最后,构建算例验证模型的可行性及有效性.结果表明:在不同的目标函数权重选择下,模型可以充分考虑电-气网络线路与多种综合能源设备间的耦合关系,获得整体最优的规划方案.

关键词: 综合能源; 多目标; 扩展规划; 低碳转型

Abstract

In order to accelerate the rapid and economic low-carbon transformation of the power-gas system, a multi-objective stochastic optimization programming model for the whole equipment of the power-gas system was established, which comprehensively considered the economic cost and carbon emissions. First, the mathematical model of the electric-gas network and related equipment was established, and the uncertainty characteristics of the electric and gas loads and photovoltaic output were analyzed by using the scenario method. Next, a mixed-integer quadratically constrained programming (MIQCP) model considering the economic cost and carbon emissions of the system was established. An overall planning was made for power feeders, gas network pipelines, substations, gas distribution stations, gas units, power-to-gas devices, photovoltaic, and energy storage devices. Finally, a numerical example was built to verify the feasibility and effectiveness of the model. The results show that the model can fully consider the coupling relationship between power-gas network lines and a variety of comprehensive energy equipment under different weight choices of objective function, and obtain the overall optimal planning scheme.

Keywords: integrated energy; multi-objective; expansion planning; low-carbon transformation

PDF (1667KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

朱海南, 王娟娟, 陈兵兵, 张厚望, 陈健, 吴秋伟. 考虑经济性与碳排放的电-气综合能源系统多目标规划[J]. 上海交通大学学报, 2023, 57(4): 422-431 doi:10.16183/j.cnki.jsjtu.2021.513

ZHU Hainan, WANG Juanjuan, CHEN Bingbing, ZHANG Houwang, CHEN Jian, WU Qiuwei. Multi-Objective Planning of Power-Gas Integrated Energy System Considering Economy and Carbon Emission[J]. Journal of Shanghai Jiaotong University, 2023, 57(4): 422-431 doi:10.16183/j.cnki.jsjtu.2021.513

近年来,随着环境问题的日益凸显,传统能源带来的污染问题逐渐显露,实现能源生产和使用的低碳化已成为人类社会可持续发展的必然趋势[1].电-气网络是能源供应的重要环节,实现电-气系统快速、经济的低碳转型是推进我国经济社会发展动力转换的重要引擎,是保障国家能源安全的重要举措.综合考虑经济性与碳排放量,对电-气网络线路及相关设备进行统筹规划,保证规划网络的经济性与灵活性,提高光伏等清洁能源的消纳能力,具有重大研究意义.

电-气系统联合规划以电网规划为基础,在电-气网络耦合日益密切的当下逐渐成为国内外研究的热点问题.文献[2]以电-气网络综合投资与运行成本最小为目标,对电-气网络线路进行联合扩展规划.文献[3]建立了计及传输线重构的分布鲁棒电-气综合能源系统扩展规划模型,以系统总投资成本和包括弃负荷成本在内的总运行成本最低为目标,对电力线路与天然气管道进行规划.文献[4-5]建立了考虑风力发电(简称风电)出力随机性的规划模型,模型聚焦于风电场、电转气装置等综合能源设备的选址定容问题,未考虑电-气网络线路的规划问题.文献[6]建立了考虑风电不确定性的综合能源系统扩展规划模型,在风电装机容量确定的前提下对电网线路与耦合设备进行规划.文献[7]建立了考虑需求侧响应的电-气系统规划模型,以建设、运营总成本和需求侧响应补偿最小为目标,对电-气网络线路和燃气机组、电转气装置进行规划.文献[8]所提的扩展规划模型考虑全球电网的廉价燃料获取和输电瓶颈问题,以低碳经济为导向,实现对电-气网络线路与燃气机组、压气泵的规划.总结上述文献可知,目前电-气网络规划模型部分聚焦于电-气网络线路、拓扑的扩展规划问题[2-3],部分聚焦于综合能源设备的选址定容问题[4-5],对同时进行线路扩展规划与设备选址定容的研究也局限于少数几种设备类型[6-8].电-气综合能源系统内部各设备相互耦合,如果规划方案仅考虑其中线路或部分类型设备的配置,势必会引起资源的不合理配置.

总结目前电-气综合能源系统规划中常见的待规划设备类型,实现了对电网馈线、气网管道以及变压器、配气站、燃气轮机、电转气装置、光伏及储能装置等多种类型设备的统筹规划.另外,规划模型综合考虑了系统的经济性成本与碳排放量,决策者可以根据需求自由选择两个目标的权重系数,获得符合需求的规划方案.

1 数学模型

1.1 综合能源系统设备数学模型

1.1.1 燃气机组

燃气机组可以在电网电能不足时通过燃烧天然气产生电能,实现天然气网向电网的能量传递.相较于传统燃煤机组,燃气机组产生单位电能时碳排放量更低,具有良好的环境效益.燃气机组数学模型如下:

Pi,d,hgt=gTgPg,i,d,hgt
Qi,d,hgt=gTgηgPg,i,d,hgt
Fm,d,hgt=gTgPg,i,d,hgtζgHgas
0Pg,i,d,hgtgTgZg,igtPgmax
gTgZg,igt1

式中: iNg, Ng为配电网中备选燃气机组节点集合; Tg为备选燃气机组类型集合; Pi,d,hgtQi,d,hgt分别表示配电网节点i处燃气机组在典型日d中时刻h输出的有功和无功功率; Fm,d,hgt表示气网节点m处燃气机组在典型日d中时刻h消耗的天然气体积; Pg,i,d,hgt表示配电网节点i处g类型燃气机组在典型日d中时刻h输出的有功功率; Zg,igt为0-1变量,表示配电网节点i处是否新建g类型燃气机组; Pgmaxζgηg分别为g类型燃气机组的额定容量、气转电效率、无功系数;Hgas为天然气热值.

式(1)~(4)表示燃气机组能量转换过程,式(5)为燃气机组建设逻辑约束.式(1)表示节点i处燃气机组的有功功率.式(2)、式(3)分别表示节点i处燃气机组有功与无功、有功与气量的转换关系. 式(4)保证各类型机组发出的有功功率小于建设的额定功率.式(5)保证节点i最多建设一种类型的备选燃气机组.

1.1.2 电转气装置

电转气装置可以将电网中的过剩电能转换为天然气,实现电网向天然气网的能量传递.在电转气装置运行过程中,产生单位体积的天然气需要单位CO2气体作为原料,因此该装置具有一定的固碳作用.电转气装置数学模型如下:

Pi,d,hp2g=pTpgPp,i,d,hp2g
Fm,d,hp2g=pTpgζpPp,i,d,hp2gHgas
0Pp,i,d,hp2gpTpgZp,ip2gPpmax
pTpgZp,ip2g1

式中: iNpg, Npg为配电网中备选电转气装置节点集合; Tpg为备选电转气装置类型集合; Pi,d,hp2g表示配电网节点i处电转气装置在典型日d中时刻h消耗的有功功率; Fm,d,hp2g表示气网节点m处电转气装置在典型日d中时刻h产生的天然气体积; Pp,i,d,hp2g表示配电网节点i处p类型电转气装置在典型日d中时刻h消耗的有功功率; Zp,ip2g为0-1变量,表示配电网节点i处是否新建p类型电转气装置; Ppmaxζp分别为p类型电转气装置的额定容量、电转气效率.

式(6)表示节点i处电转气装置的有功功率.式(7)表示节点i处电转气装置有功功率与天然气流量的转换关系.式(8)保证各类型电转气装置发出的有功功率小于各自的额定功率.式(9)保证节点i最多建设一种类型的备选电转气装置.

1.1.3 储能装置

在电-气综合能源系统中,储能装置通过在不同时刻的充放电调节网络供需关系,提高系统光伏利用率.由于储能装置充放电效率小于1,所以最优解将自动满足充放电互补松弛[9],基于此建立的不含充放电标志位的储能装置数学模型如下:

eTeZe,iessEeminEi,d,hesseTeZe,iessEemax
0Pi,d,hess,chaeTeZe,iessPemax
0Pi,d,hess,diseTeZe,iessPemax
Ei,d,hess=Ei,d,h-1ess+μPi,d,hess,cha-Pi,d,hess,dis/μ
Ei,d,24ess=Ei,d,0ess
eTeZe,iess1

式中: iNes, Nes为配电网中备选储能装置节点集合; Te为备选储能装置类型集合; Ze,iess为0-1变量,表示配电网节点i处是否新建e类型储能装置; Pemax为e类型储能装置的功率上限; EeminEemax分别为e类型储能装置的最小容量、最大容量; Ei,d,hessPi,d,hess,chaPi,d,hess,dis分别表示配电网节点i处储能装置在典型日d中时刻h时的电量、充电功率、放电功率;μ为储能装置的充放电效率.

式(10)为储能装置容量约束.式(11)、式(12)分别表示充电功率、放电功率约束.式(13)为储能装置电量变化与充、放电功率间的关系约束.式(14)保证典型日内储能装置电量的初、末状态相同.式(15)为储能装置建设逻辑约束.

1.1.4 分布式光伏

分布式光伏建设与运行约束如下:

vTvZv,ipv1
0Pi,d,hpvvTvZv,ipvPvmaxP^d,hpv

式中: iNv, Nv为配电网中备选光伏电站节点集合; Tv为备选光伏电站类型集合; Zv,ipv为0-1变量,表示配电网节点i处是否新建v类型光伏电站; Pi,d,hpv表示配电网节点i处光伏电站在典型日d中时刻h的输出功率; Pvmax表示v类型光伏电站的装机容量; P^d,hpv表示单位容量光伏电站在典型日d中时刻h的输出功率.

式(16)为光伏电站建设逻辑约束,式(17)为光伏电站各时刻输出的有功功率约束.

1.2 配电网模型

1.2.1 变电站

变电站建设与运行约束如下:

(Pi,d,hsub)2+(Qi,d,hsub)2sTsZs,isub(Ssmax)2
Zs,isub=Zn,s,isub+Zc,s,isub
sTsZs,isub1

式中: iNs, Ns为配电网中变电站节点集合; Ts为备选变电站类型集合; Ssmax表示s类型变电站的额定容量; Pi,d,hsubQi,d,hsub分别表示节点i处变电站在典型日d中时刻h输出的有功和无功功率; Zn,s,isubZc,s,isub为0-1变量,分别表示配电网备选变电站节点i处是否存在新建的、扩建的s类型变电站; Zs,isub为0-1变量,表示规划后节点i处是否存在s类型变电站.

式(18)为变电站有功功率、无功功率与其额定容量的关系约束.式(19)、式(20)保证在配电网节点i处仅存在一种类型的变电站.

1.2.2 配电网潮流

配电网潮流模型采用适合于辐射状网络的线性交流DistFlow模型[10],模型如下:

Vi,d,h-Vj,d,h=cTc(rcPc,ij,d,hcod+xcQc,ij,d,hcod)LijVref+bij,d,h
VminVi,d,hVmax
jiBePji,d,hbra+Pi,d,hgt+Pi,d,hpv+Pi,d,hess,dis=ijBePij,d,hbra+Pi,d,hp2g+Pi,d,hess,cha+PiloadS^d,h
jiNeQji,d,hbra-ijNeQij,d,hbra+Qi,d,hgt+QiloadS^d,h=0
Pij,d,hbra=cTcPc,ij,d,hcod
Qij,d,hbra=cTcQc,ij,d,hcod
Pc,ij,d,hcodZc,ijcodPcod,cmax
Qc,ij,d,hcodZc,ijcodQcod,cmax
bij,d,h(Vmax-Vmin)(1-cTcZc,ijcod)
Zc,ijcod=Zn,c,ijcos+Zc,c,ijcod
cTcZc,ijcod1

式中: ijBeiNecTc, Be为配电网中备选线路集合; Ne为配电网节点集合; Tc为备选传输线类型集合; Vi,d,h表示节点i在典型日d中时刻h的电压幅值; Vref为参考节点电压; Pc,ij,d,hcodQc,ij,d,hcod分别表示配电网线路ij处在典型日d中h时刻c类型馈线流过的有功、无功功率; Pij,d,hbraQij,d,hbra分别表示在典型日d中h时刻流过线路ij的有功、无功功率; PiloadQiload分别表示节点i的有功、无功功率基值; S^d,h表示在典型日d中h时刻的归一化电负荷; bij,d,h为辅助松弛变量; rcxc分别表示e类型导线的单位长度电阻值、电抗值;Lij 表示线路ij的长度; VmaxVmin分别表示节点电压最大值与最小值; Zc,ijcod为0?1变量,表示规划后配电网线路ij处是否存在c类型馈线; Zn,c,ijcodZc,c,ijcod为0-1变量,分别表示配电网线路ij处是否存在新建的、扩建的c类型馈线; Pcod,cmaxQcod,cmax表示流过c类型馈线的最大有功、无功功率.

式(21)~(29)为配电网运行约束,式(30)~(31)为馈线建设逻辑约束.式(21)为简化的DistFlow潮流方程.式(22)为节点电压上下限约束.式(23)、式(24)为节点有功、无功功率平衡约束.式(25)、式(26)表示流过线路的有功、无功功率与各类型馈线中有功、无功功率的关系.式(27)、式(28)为各类型馈线有功、无功功率上下限约束.式(29)约束了变量bij,d,h的取值,当线路ij处不建设馈线时, cTcZc,ijcod=0,bij,d,h为在[Vmin-Vmax,Vmax-Vmin] 范围内随意变动的辅助变量,此时代表线路两端电压与流过线路功率间关系的约束式(21)被松弛;反之, 当线路ij存在馈线时, bij,d,h=0,式(21)正常存在.式(30)、式(31)保证在配电网线路ij处仅存在一种类型的馈线.

1.2.3 放射状拓扑

为保证规划后配电网为放射状拓扑,需要引入放射状拓扑约束.采用潮流剥离配电网规划方法,通过引入非负连续辅助变量yij+yij-, 保证规划后配电网为放射状拓扑[11].建立模型如下:

yij+0, yij-0
yij++yij-1
cTcZc,ijcod=yij++yij-
kiBeyki++ijBeyij-=1-sTsZs,isub1, iNs
kiBeyki++ijBeyij-=1, iNe-Ns

式中: ijBe.yij+yij-为虚拟潮流方向辅助变量,通过式(32)~(36)的共同作用, yij+yij-的值被约束为0或1. yij+=1代表存在由节点i流向节点j的虚拟潮流, yij+=0代表不存在; yij-=1代表存在由节点j流向节点i的虚拟潮流, yij-=0代表不存在.

式(32)、式(33)约束连续辅助变量yij+yij-的取值范围.式(34)保证当线路ij处存在馈线时yij+yij-之和为1,反之为0.式(35)保证潮流不从配电网其他节点倒流向变电站.式(36)保证配电网除变电站外的任意节点都有且仅有一个父节点.

1.3 天然气网模型

1.3.1 配气站

配气站建设与运行约束如下:

0Fm,d,hgataTaFamaxZa,mgat
Za,mgat=Zn,a,mgat+Zc,a,mgat
aTaZa,mgat1

式中: mNa, Na为天然气网中备选配气站节点集合; Ta为备选配气站类型集合; Famax表示a类型配气站的额定容量; Fm,d,hgat表示天然气网节点m处配气站在典型日d中时刻h输出的气量; Zn,a,mgatZc,a,mgat为0-1变量,分别表示节点m处是否存在新建的、扩建的a类型配气站; Za,mgat为0-1变量,表示规划后节点m处是否存在a类型配气站.

式(37)为配气站流量上下限约束,式(38)、式(39)为配气站建设逻辑约束.

1.3.2 天然气网

忽略天然气管道Weymouth方程、压缩机方程,建立如下气网模型[12]:

mnBtFmn,d,hbra+Fm,d,hp2g=nmBtFmn,d,hbra+Fm,d,hgt+FmloadF^d,h
Ff,mn,d,hpipZf,mnpipFfmax
Fmn,d,hbra=fTfFf,mn,d,hpip
Zf,mnpip=Zn,f,mnpip+Zc,f,mpip
fTfZf,mnpip1

式中: mNtmnBtfTf, Nt为气网中备选管道集合,n为天然气网节点序号; Bt为气网节点集合; Tf为备选管道类型集合; Ff,mn,d,hpip表示气网线路mn处在典型日d中h时刻f类型管道流过的流量; Fmn,d,hbra表示在典型日d中h时刻流过线路mn的气量; Fmload表示节点m的气负荷基值; F^d,h表示在典型日d中h时刻的归一化气负荷; Zf,mnpip为0-1变量,表示规划后气网线路mn处是否存在f类型管道; Zn,f,mnpipZc,f,mpip为0-1变量,分别表示气网线路mn处是否存在新建的、扩建的f类型管道;Ffmax表示f类型管道最大传输容量.

式(40)~(42)为天然气网运行约束,式(43)~(44)为馈线建设逻辑约束.

1.3.3 放射状拓扑

与电网放射状拓扑约束相同,将ymn+ymn-设置为连续变量,获得气网放射状网络拓扑约束如下:

ymn+0, ymn-0, mnBt
fTfZf,mnpip=ymn++ymn-, mnBt
ymn++ymn-1, mnBt
lmBtylm++mnBtymn-=1-aTaZa,mgat,mNa
lmBtylm++mnNeymn-=1, mNt-Na

2 规划模型

2.1 目标函数

目标函数由系统经济性成本f1与碳排放量f2两部分组成.经济性成本f1包括设备投资成本Cinv与系统运行成本Cope两部分,其中投资成本折算至一年,运行成本为一年的购电成本与购气成本;碳排放量f2包括传统煤炭发电产生的碳排放、燃气机组发电产生的碳排放以及电转气装置运行过程中减少的碳排放.

f1f2的具体表达形式如下:

f1=Cinv+Cope

Cinv=mnBtfTfκ(Cf,npipLmnZn,f,mnpip+Cf,cpipLmnZc,f,mnpip)+

ijBecTcκ(Cc,ncodLijZn,c,ijcod+Cc,ccodLijZc,c,ijcod)+

mNtaTaκ(Ca,ngatZn,a,mgat+Ca,cgatZc,a,mgat)+

iNssTsκ(Cs,nsubZn,s,isub+Cs,csubZc,s,isub)+

iNggTgκCggtZg,igt+iNpgpTpgκCpp2gZp,ip2g+

iNvvTvκCvpvZv,ipv+iNeseTeκCeessZe,iess
Cope=CeiNsd(ndhPi,d,hsub)+CamNad(ndhFm,d,hgat)
f2=WcηciNsd(ndhPi,d,hsub)+WgtiNgd(ndhPi,d,hgt)-Wp2giNpgd(ndhPi,d,hp2g)
κ=r1-(1+r)-T

式中: Cf,npipCf,cpip分别表示单位长度f类型天然气管道新建、扩建成本; Cc,ncodCc,ccod分别表示单位长度c类型馈线新建、扩建成本; Ca,ngatCa,cgat分别表示a类型配气站新建、扩建成本; Cs,nsubCs,csub分别表示s类型变电站新建、扩建成本; CggtCpp2gCvpvCeess分别表示g类型燃气机组、p类型电转气装置、v类型光伏电站、e类型储能装置建设成本; CeCa分别表示单位电价、单位气价; nd表示典型日d在一年中出现的时间; ηc为变电站从外界电网获得的电能中传统煤炭发电所占的比例,本文取1;WcWgt分别表示传统煤炭发电、燃气机组发电生产单位电能时产生的CO2气体质量;Wp2g表示电转气装置产生单位天然气时消耗的CO2气体质量; κ表示设备的年投资成本折算系数; r表示年利率;T表示规划年限.

经济性成本f1与碳排放量f2由于单位、数量级上的差异无法直接相加,为了协调两个目标函数,本文首先进行模糊化处理将目标函数归一化,分别表示出各个目标函数对各自解的满意度,再通过加权满意度指标法求得最优的折中解,最终的目标函数如下:

minω1f1-f1,minf1,max-f1,min+ω2f2-f2,minf2,max-f2,min

式中: ω1ω2为权重系数,满足ω1+ω2=1,决策者可根据需求确定二者权重;f1,maxf1,minf2,maxf2,min分别表示f1f2的最大值、最小值,其中f1,minf2,max为规划目标仅考虑系统经济性成本f1时规划期内系统的经济性成本、碳排放量,f1,maxf2,min为规划目标仅考虑系统碳排放量f2时规划期内系统的经济性成本、碳排放量.

2.2 约束条件

约束条件包括设备约束式(1)~(17),配电网约束式(18)~(36)及气网约束式(37)~(49)等3个部分,其中式(18)为二次约束,其余均为线性约束.

本文所建立规划模型为混合整数二次约束规划(MIQCP)模型,可调用商用求解器Gurobi求解.

3 算例分析

3.1 算例介绍

为验证模型的有效性,构建电-气网络联合扩展规划算例,待扩建电-气综合能源系统[13]图1所示.

图1

图1   电-气综合能源系统[13]

Fig.1   Integrated power and natural gas system[13]


备选配电网馈线、天然气管道参数分别如表1表2所示,网络中原有线路均为1型线路.配电网节点22、23为变电站节点,为原有1型变电站;配气网节点26为配气站节点,为原有1型配气站;备选变电站与配气站参数分别如表3表4所示.配电网节点6、18分别与配气网节点30、3耦合,为备选燃气机组与电转气装置节点;备选燃气机组、电转气装置参数分别如表5表6所示.配电网节点2、4、11为备选光伏与储能装置节点;备选光伏电站、储能装置参数分别如表7表8所示.各备选设备以容量为依据进行分类,型号越大,对应容量越大,即2型设备容量大于1型,3型大于2型.单位电价Ce、单位气价Ca分别取0.5 元/(kW·h)、3.95 元/m3.碳排放量系数[14]WcWgtWp2g分别取0.623 kg/(kW·h),0.184 kg/(kW·h),1.96 kg/m3.电压上下限VmaxVmin分别取1.05(p.u.)、0.95(p.u.).规划年限T=10 a,年利率r=5%.

表1   备选配电网馈线参数

Tab.1  Parameters of alternate feeder

类型Cc,ccod/
(元·km-1)
Cc,ncod/
(元·km-1)
Pcod,cmax/MWQcod,cmax/Mvarrc/(Ω·km-1)xc/(Ω·km-1)
111.526.443.390.822 30.567
211.52×10422.26×10413.67.160.4110.283 5

新窗口打开| 下载CSV


表2   备选天然气网管道参数

Tab.2  Parameters of alternate piping

类型Cf,cpip/
(元·km-1)
Cf,npip/
(元·km-1)
Ffmax/(m3·h-1)
114.88420
235.65×10437.14×1042 600

新窗口打开| 下载CSV


表3   备选变电站参数

Tab.3  Parameters of alternate transformer substation

类型Cs,csub/元Cs,nsub/元Ssmax/( MV·A)
1198.0×10410
2185.7×104297.1×10415
3267.3×104445.7×10420

新窗口打开| 下载CSV


表4   备选配气站参数

Tab.4  Parameters of alternate gas distributing station

类型Ca,cgat/元Ca,ngat/元Famax/(m3·h-1)
1399.0×1043 200
2331.4×104533.7×1044 200
3468.0×104668.5×1045 200

新窗口打开| 下载CSV


表5   备选燃气机组参数

Tab.5  Parameters of alternate gas generator

类型ζgηgCggt/(元· MW-1)Pgmax/MW
10.350.3440.24×1041.08
20.400.3483.46×1042.62

新窗口打开| 下载CSV


表6   备选电转气装置参数

Tab.6  Parameters of alternate power to gas device

类型ζpCpp2g/元Ppmax/MW
10.60666.4×1041
20.621 928.0×1043

新窗口打开| 下载CSV


表7   备选光伏电站参数

Tab.7  Parameters of alternate PV

类型Ppmax/MWCvpv/元Np
1102 704×1042、4、11
2154 056×1042、4、11
3205 406×1042、4、11

新窗口打开| 下载CSV


表8   备选储能装置参数

Tab.8  Parameters of alternate energy storage device

类型Eemax/
(MW·h)
Pemax/
MW
Ceess/
Nesμ
151.75883.2×1042、4、110.95
272.451 236.0×1042、4、110.95
3103.51 766.0×1042、4、110.95

新窗口打开| 下载CSV


利用场景法表征电负荷、气负荷与光伏出力的不确定性特征[15-16],夏季典型日、冬季典型日及春秋过渡季典型日的数据如图2~4所示.其中夏季、冬季、春秋过渡季典型日的权重分别取0.25、0.25、0.50,即夏季、冬季、春秋过渡季典型日在一年中出现的时间nd分别取91、91、183 d.

图2

图2   归一化电负荷变化

Fig.2   Curves of normalized electrical load


图3

图3   归一化气负荷变化

Fig.3   Curves of normalized gas load


图4

图4   单位容量光伏出力变化

Fig.4   Variation curves of photovoltaic output per unit capacity


为验证模型有效性,设置如下对照算例:

算例1,兼顾系统经济成本与碳排放量,权重系数取ω1=0.5,ω2=0.5.

算例2,只考虑系统经济性成本,权重系数取ω1=1,ω2=0.

算例3,只考虑系统碳排放量,权重系数取ω1=0,ω2=1.

3.2 结果分析

各算例中综合能源设备选址定型结果如表9所示,表中e、c、n分别代表原有设备、扩建设备与新建设备,如e1代表原有1型设备.

表9   各方案设备选型比较

Tab.9  Equipment selection of each case

设备名称位置设备类型
算例1算例2算例3
光伏电站2、4、11n3、n1、n3n1、n1、n2n3、n3、n3
储能装置2、4、11n3、n3、n30、0、0n3、n3、n3
燃气机组6、18n2、n2n1、n1n2、n2
电转气6、180、00、0n2、n2
变电站22、23c2、e1c2、e1c3、c3
配气站26c2c2c3

新窗口打开| 下载CSV


对比算例1与算例2规划结果,可以发现算例1中新建了两处3型光伏电站与1处1型光伏电站,大于算例2中光伏电站的建设容量.这是由于算例1兼顾系统碳排放量,所以倾向于建设更多的光伏电站增加系统清洁能源渗透率.与此相应,算例1中建设了3处3型储能装置,在光伏发电过剩时吸收多余电能,在光伏出力较少时发出电能,以此来平抑光伏出力的波动性,提高光伏的电能利用率,而光伏电站容量较小的算例2中没有建设储能装置.

另外,算例1在节点6、18两处分别规划了2型燃气机组,而算例2仅规划了1型燃气机组.这是由于在电-气网络运行过程中,燃气机组有两方面的作用:一是在电网电能不足、气网天然气过盛时将天然气转换为电能,提高电-气网络内设备利用率,降低系统整体的经济性成本;二是替代传统化石能源发电,降低系统的碳排放量.算例2仅考虑到燃气机组在降低系统经济性成本上的作用,未考虑到其在降低碳排放量方面产生的效益,因此规划结果中燃气机组容量较算例1较小.

各算例中配电网、天然气网线路扩展规划结果如图5所示.不同算例中线路扩展规划结果呈现出较大的差异性,这是不同权重系数下综合能源设备建设容量的差异导致的.在配电网与天然气网中,各综合能源设备扮演着分布式源、荷的角色,其容量的不同导致网络中源、荷分布的差异.这种差异会影响网络最优潮流分布,进而使线路扩展规划呈现不同结果.

图5

图5   各算例线路扩展规划结果

Fig.5   Expansion planning results of cases


各算例规划方案的经济性成本及碳排放量如表10所示.算例1所获得的规划方案年投资、运行产生的总经济成本为1.678×108元,较之仅考虑经济成本的算例2增加了1.520 ×107元;而算例1中系统年碳排放量为4.378×107 kg,较之算例2减少了2.403×107 kg.简而言之,算例1中规划方案以每年增加1.521×109元成本为代价,减少了CO2气体的排放2.403×107 kg.

表10   各方案成本比较

Tab.10  Cost of each case

算例投资成本/元运行成本/元总经济成本/元碳排放量/kg
算例13.422×1071.336×1081.678×1084.378×107
算例21.854×1071.341×1081.526×1086.781×107
算例34.724×1071.395×1081.867×1083.517×107

新窗口打开| 下载CSV


对比算例1与算例3规划结果,可以发现仅考虑碳排放量的算例3的规划方案中所有设备均进行了最大容量的建设.算例3所获得的规划方案相较于算例1每年增加了1.890×107元的成本,减少了8.610×106 kg CO2气体的排放.

通过以上分析可以看出,电-气系统线路与各种综合能源设备在规划运行过程中相互耦合、相互影响,在规划过程中对线路及多种设备进行统一、合理的规划才能获得最优的规划结果.本文所提出的模型可以在不同的目标函数权重选择下,综合考虑系统经济性与碳排量,实现对电-气网络线路以及变电站、配气站、光伏电站、储能装置、燃气机组、电转气装置的统筹配置.

4 结语

合理的规划方案是提高电-气综合能源系统经济性、降低系统整体碳排放量的基础.提出一种综合考虑经济性成本与碳排放量的电-气综合能源系统多目标随机优化规划模型,实现对电网馈线、气网管道以及变压器、配气站、燃气轮机、电转气装置、光伏及储能装置等多种类型综合能源设备的统筹规划.算例结果表明:① 通过经济性成本与碳排量二者权重系数的选择可以获得不同规划要求下的电-气综合能源系统规划方案;② 规划模型可以统筹考虑电-气网络线路与多类型综合能源设备的耦合关系,获得待规划系统的最优规划方案.但是,本文对气网的建模不够精确,且没有考虑热网模型,因此考虑更精细的气网模型并将热网引入规划模型将是下一步的工作方向.

参考文献

奚剑明, 吴瀚然.

供给侧改革中资源与能源行业去产能政策效果研究

[J]. 工业技术经济, 2021, 40(1): 113-119.

DOI:10.3969/j.issn.1004-910X.2021.01.013      [本文引用: 1]

中国经济历经复苏和快速增长阶段后,部分产业出现产能过剩现象,产业结构的不均衡变化导致了GDP增速放缓。2016年初中国政府提出供给侧结构性改革,旨在通过 “三去一降一补”政策调整产业结构,刺激经济发展。本文构建了一个包含资源与能源产业的一般均衡模型,量化分析中国供给侧结构性改革的政策效果。研究表明,在去产能政策削减上游资源产业的产出后,政策冲击会传导至下游能源产业,导致化石能源相对于绿色能源产量下降14.04%,能源结构实现优化。随后的政策对比实验引入了新能源补贴、碳排放税征收及资源税收缴等财政手段,与去产能政策进行效果对比。实验结果显示,3种政策在不同程度上实现了能源结构调整,其中征收碳税的政策使化石能源占比降幅最大,但同时会使资本利率上升,加重企业融资成本;新能源补贴与资源税的政策也可以实现调控目标,但与去产能政策相比效果更弱。

XI Jianming, WU Hanran.

Research on the Effect of supply side reform and the de-capacity policy in resource and energy industry

[J]. Journal of Industrial Technological Economics, 2021, 40(1): 113-119.

[本文引用: 1]

UNSIHUAY C, MARANGON-LIMA J W, SOUZA A C Z D.

Integrated power generation and natural gas expansion planning

[C]// Hans-Björn Püttgen. 2007 IEEE Lausanne Power Tech. New York, USA: IEEE, 2007: 1-6.

[本文引用: 2]

魏震波, 郭毅, 魏平桉, .

考虑传输线重构的电气综合能源系统分布鲁棒扩展规划模型

[J]. 电力自动化设备, 2021, 41(2): 16-23.

[本文引用: 2]

WEI Zhenbo, GUO Yi, WEI Ping’an, et al.

Distribution robust expansion planning model for integrated natural gas and electric power systems considering transmission switching

[J]. Electric Power Automation Equipment, 2021, 41(2): 16-23.

[本文引用: 2]

薛友, 李杨, 高滢, .

计及风电出力随机特性的电-气综合能源系统随机优化

[J]. 电力建设, 2018, 39(12): 2-12.

DOI:10.3969/j.issn.1000-7229.2018.12.001      [本文引用: 2]

不断成熟的电转气 (power to gas, P2G) 技术和燃气机组为电力系统和天然气系统的双向耦合和闭环运行, 进而形成电-气综合能源系统 (integrated electricity-gas energy system, IEGES) 奠定了基础。IEGES的发展有助于提升电力系统和天然气系统运行的灵活性, 并为消纳风电等间歇性可再生能源发电提供新途径。在此背景下, 研究了含P2G装置和燃气机组的IEGES的随机优化策略。首先, 构建了计及风电出力随机特性的IEGES随机优化模型, 以运行成本最小化为优化目标, 并考虑了电力系统和计及管存影响的天然气系统的相关运行约束。其中, 利用快速搜索密度聚类算法对历史风速数据进行场景聚类, 以处理风电出力的随机特性。然后, 采用AMPL/IPOPT商业求解器求解所建立的随机优化模型。最后, 以修改的IEEE 39节点电力系统和比利时20节点天然气系统所构成的IEGES为例, 对所构建的优化模型和采用的求解方法进行说明。

XUE You, LI Yang, GAO Ying, et al.

Stochastic optimization in integrated electricity-gas energy systems considering stochastic characteristics of wind power outputs

[J]. Electric Power Construction, 2018, 39(12): 2-12.

DOI:10.3969/j.issn.1000-7229.2018.12.001      [本文引用: 2]

 The power to gas (P2G) technology and natural gas generation units (NGGUs) facilitate the bidirectional coupling and hence closed-loop operation between a power system and a natural gas system concerned, and the development of an integrated electricity-gas energy system (IEGES). IEGES can help promote the operation flexibility of the power system and natural gas system concerned and offer a potential alternative to accommodate intermittent generation from wind power and other renewable energy sources. Given this background, stochastic optimization strategies in IEGES with P2G facilities and NGGUs are addressed. First, a stochastic optimization model of IEGES with stochastic characteristics of wind power output is presented, with an objective of minimizing the overall operation costs, and the operation constraints in both the power system and natural gas system. Among them, by employing a fast searching density clustering algorithm, historical wind speed data are clustered so as to simulate the stochastic characteristics of wind power outputs. Then, the commercial solver AMPL/IPOPT is used to solve the presented stochastic optimization model. Finally, a sample IEGES with a modified version of the IEEE 39-node power system and the Belgium 20-node natural gas system included is employed to demonstrate the proposed method.

王芃, 刘伟佳, 林振智, .

基于场景分析的风电场与电转气厂站协同选址规划

[J]. 电力系统自动化, 2017, 41(6): 20-29.

[本文引用: 2]

WANG Peng, LIU Weijia, LIN Zhenzhi, et al.

Scenario analysis based collaborative site selection planning of wind farms and power-to-gas plants

[J]. Automation of Electric Power Systems, 2017, 41(6): 20-29.

[本文引用: 2]

李哲, 王成福, 梁军, .

计及风电不确定性的电-气-热综合能源系统扩展规划方法

[J]. 电网技术, 2018, 42(11): 3477-3487.

[本文引用: 2]

LI Zhe, WANG Chengfu, LIANG Jun, et al.

Expansion planning method of integrated energy system considering uncertainty of wind power

[J]. Power System Technology, 2018, 42(11): 3477-3487.

[本文引用: 2]

GAO Y, WEN F S, WANG K, et al.

Optimal collaborative planning with demand side management in integrated gas-electricity energy systems

[C]// 2018 IEEE Power & Energy Society General Meeting. New York, USA: IEEE, 2018: 1-5.

[本文引用: 2]

JING Q, DONG Z Y, ZHAO J H, et al.

Low carbon oriented expansion planning of integrated gas and power systems

[J]. IEEE Transactions on Power Systems, 2015, 30(2): 1035-1046.

DOI:10.1109/TPWRS.2014.2369011      URL     [本文引用: 2]

李捷, 余涛, 潘振宁.

基于强化学习的增量配电网实时随机调度方法

[J]. 电网技术, 2020, 44(9): 3321-3332.

[本文引用: 1]

LI Jie, YU Tao, PAN Zhenning, et al.

Real-time stochastic dispatch method for incremental distribution network based on reinforcement learning

[J]. Power System Technology, 2020, 44(9): 3321-3332.

[本文引用: 1]

BARAN M, WU F.

Optimal sizing of capacitors placed on a radial distribution system

[J]. IEEE Transactions on Power Delivery, 1989, 4(1): 735-743.

DOI:10.1109/61.19266      URL     [本文引用: 1]

ZHOU X Z, GUO C X, WANG Y F, et al.

Optimal expansion co-planning of reconfigurable electricity and natural gas distribution systems incorporating energy hubs

[J]. Energies, 2017, 10(1): 124.

DOI:10.3390/en10010124      URL     [本文引用: 1]

ZHANG X P, SHAHIDEHPOUR M, ALABDULWAHAB A, et al.

Optimal expansion planning of energy hub with multiple energy infrastructures

[J]. IEEE Transactions on Smart Grid, 2017, 6(5): 2302-2311.

DOI:10.1109/TSG.2015.2390640      URL     [本文引用: 1]

SALDARRIAGA C A, HINCAPIE R A, SALAZAR H, et al.

A holistic approach for planning natural gas and electricity distribution networks

[J]. IEEE Transactions on Power Systems, 2013, 28(4): 4052-4063.

DOI:10.1109/TPWRS.2013.2268859      URL     [本文引用: 3]

栗然, 申雪, 钟超, .

考虑环境效益的分布式电源多目标规划

[J]. 电网技术, 2014, 38(6): 1471-1478.

[本文引用: 1]

LI Ran, SHEN Xue, ZHONG Chao, et al.

Multi-objective planning of distributed generation considering environmental benefit

[J]. Power System Technology, 2014, 38(6): 1471-1478.

[本文引用: 1]

胡源, 薛松, 杨素, .

综合能源背景下的配电网多场景规划

[J]. 中国电力, 2021, 54(4): 175-184.

[本文引用: 1]

HU Yuan, XUE Song, YANG Su, et al.

Multi-scenario planning of distribution network under integrated energy background

[J]. Electric Power, 2021, 54(4): 175-184.

[本文引用: 1]

SAROJ, KAVITA.

Review: Study on simple k mean and modified K mean clustering technique

[J]. International Journal of Computer Science Engineering and Technology, 2016, 6(7): 279-281.

[本文引用: 1]

/