李红伟, 荆浩婕, 吴 磊, 李婷玉
(西南石油大学 电气信息学院,四川 成都 610500)
近年来,随着储能技术、电转热技术的逐渐成熟,出现了能源互联网、多能源系统(multi-energy system,MES)的理念。MES建模作为MES运行的基础显得尤为重要,其中能够转换、调节和存储多个能源载体的能量枢纽(energy hub,EH)成为MES建模的关键单元[1-6]。含热电联供(combined heat and power,CHP)的EH可以将配电网(power distribution network,PDN)和区域供热网(district heating network,DHN)耦合,从而解决热能易储存、不易运输以及电能易运输、不易储存的问题[3-6]。
虽然已有学者探讨了基于EH模型的MES的优化规划[1-2]、优化运行[3-4]以及碳排放[5-6]等问题,但针对MES标准化建模的研究很少。常雨芳等[7]根据建立的含风力发电的冷热电联供系统,利用能量的耦合关系将EH设定成以矩阵建模为基础的运行方式,但这种建模方式只适用于结构简单的EH。Wang等[8]提出了一种自动计算EH标准化矩阵建模的方法,利用图论将EH设备特性转化为矩阵形式,进而建立不考虑调度因子的能量流线性方程组进行调度。Wang等[9]提出了用于制定EH中能量转换的自动线性化建模方法,并基于EH的耦合矩阵及其秩分析了EH的灵活性。以上研究将EH设备能量效率视为常数,或简化MES的输入和输出为线性关系,而实际工程中设备的运行效率会随运行条件发生变化。
有一些学者考虑了设备能量效率变化对MES的影响。陈晚晴等[10]将EH设备效率和负载率数据拟合为效率修正模型,反映设备负载率对效率的影响,基于修正模型建立以经济性为目标的MES优化调度模型。李虹等[11]基于MES中EH常用设备效率与负载率之间的函数关系构建EH矩阵模型,从而对MES优化配置进行研究。Huang等[12]提出了基于可变能量效率的EH标准化矩阵建模方法以研究MES的运行优化。Zheng等[13]建立了设备恒定能量效率和设备能量效率变化的2种EH模型,通过比较两者的运行成本及相对误差,证明了设备能量效率变化特性对EH优化运行的影响不能忽略。Chen等[14]为了研究EH设备的能效优化有效性,提出了一种由考虑设备能量效率变化特性的EH、天然气网络和PDN构成的MES调度方法。但以上研究均未考虑可变EH能效模型与电热网相耦合的网络约束。
为研究可变能量效率特性对电热网运行的影响,本文将PDN和DHN与考虑设备能量效率随输入或输出功率变化的EH相结合,建立了电热网优化运行线性模型。首先,对EH设备可变能量效率进行了分段线性化,建立了含CHP、储热设备及电热泵的EH标准化矩阵模型;其次,将所建立的矩阵模型应用到电热网中,使得EH在电热网中发挥节能消纳、提供电能和热能的作用,构建了含EH的电热网优化运行模型;最后,通过算例验证了本文模型的有效性。
本文所研究的EH结构如图1所示,新能源输入为风电;能量转换设备为背压运行模式的CHP和电热泵;在CHP侧配置储热设备实现热电解耦。根据图论将能量设备定义为节点,EH的输入输出被视为特殊的节点,每个节点的输入输出定义为端口,能量流用分支v表示。定义风电输入为节点1;电能输入为节点2;天然气输入为节点3;CHP为节点4;储热设备为节点5;电热泵为节点6;电输出为节点7;热输出为节点8。
图1 EH结构示意图
通过图论将设备特性转换为矩阵形式,此EH的输入向量Vin和输出向量Vout[8]分别为
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
分段线性化后的ΔEg表示为
(11)
储能荷能状态约束[8]为
Eg(t)=Eg(t-1)+ΔEg(t);
(12)
(13)
(14)
进一步引入分线器和集线器,促进分段线性化以后能量流的分离和合并,保证线性化前后结果一致[12]。分段线性化过程相当于将非线性能量设备转换为分线器、线性能量设备、集线器3层结构。
若单输入多输出设备每个输出端口的能量输出与能量输入成比例(如本文所采用的CHP),则以类似于单输入单输出设备的方式分别分段线性化从输入端口到每个输出端口之间的能量转换过程。
图2 分段线性化后的EH
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
节点5扩充端口分支关联矩阵为
(23)
若节点g表示线性设备,则
(24)
Ag(k,b2)=
(25)
(26)
如图2中的CHP将1个输入3个输出各分成了两段,分别以η1,1、η1,2这2个恒定效率表示电效率,η2,1、η2,2、η3,1、η3,2这4个恒定效率表示热效率,η1,s×s=diag(η1,1,η1,2,…,η1,s),可建立节点4的特性矩阵。同理建立节点5特性矩阵为
(27)
线性设备的特性矩阵为
(28)
其中,Hg行数是能量转换过程的个数,列数是端口的个数,矩阵中的元素值由式(28)确定。由此可以得到图2中EH节点6的特性矩阵。
(29)
(30)
(31)
(32)
EH的综合能量流方程可作为能量流约束:
(33)
同时还应满足设备的容量约束:
(34)
电热网通常通过含CHP的EH紧密耦合,因此本文采用水力-热力模型描述DHN潮流分布,采用DistFlow描述PDN潮流模型[16]。
DHN一般由热源、热负荷、供回水网组成。对DHN进行水力分析,水力分析用于描述各管道水质量流量及各节点的注入水流的水质量流量[16-18]:
(35)
(36)
(37)
热力分布描述管道和节点的温度与水质量流量的关系,其热力分布可建模[17-18]为
(38)
(39)
(40)
(41)
(42)
(43)
DHN中热负荷与热源的热功率满足:
(44)
(45)
(46)
PDN一般为辐射式结构,可采用如下线性DistFlow潮流模型对潮流分布进行描述[19]。
(47)
(48)
Uj,t=Ui,t-(rijPij,t+xijQij,t)/U0;
(49)
(50)
本文以天然气价格与分时电价为基准,以全天系统购能总成本最少为目标函数,求解含EH的电热网优化运行问题。目标函数为
(51)
本文将含EH的电热网优化运行问题转换为混合整数线性规划问题,使用MATLAB中YALMIP建模、CPLEX求解器求解[21]。算例测试系统如图3所示,为了维持PDN中各母线的电压质量,在PDN中配置无功补偿装置[17,22]。PDN各母线功率比由MATPOWER33节点PDN标准潮流求得。DHN各节点的热负荷比和质量流量、电负荷、热负荷、风电预测数据等参考文献[17]。本算例根据文献[8]建立节点表和分支表用于EH的标准化输入,设定天然气价格为1.54 元/m3,分时电价如表1所示,EH设备参数如表2所示,其中储热设备的最大储热功率和放热功率均为1 600 kW。
表1 分时电价
表2 EH设备参数
图3 测试系统图
如图4所示,当分段数从2增加到50时,预测的成本迅速增加,同时求解时间略微增加;当分段数从50到100时,运行成本逐渐趋于稳定,但求解时间大幅增加。由于分段数越大,优化结果就越接近实际运行成本,因此,本文以分段数为100的运行成本为参考值计算各个分段数的相对误差。考虑到预测成本的准确性与程序求解的效率,最终选取分段数为50的EH对考虑设备可变能量效率后的电热网情况进行分析,与考虑设备恒定能量效率情况(分段数为0)相比,运行成本相对误差从12.71%(分段数为0)降到了0.03%(分段数为50)。
图4 不同分段数的运行成本和运行时间图
EH的分段数为50时PDN节点电压如图5所示,DHN管道温度如图6所示。电价峰时(11:00)和谷时(05:00)都能维持在额定电压幅值1.00的±5%以内。通过调节回水温度能满足热负荷峰值(05:00)和非峰值(15:00)的热功率需求。在供水温度相同的条件下,回水温度越小,供热量越大,表明分段数为50的EH能保证电热网的正常运行。
图5 PDN节点电压
图6 DHN的管道温度
为了突出本文可变能量效率模型的优势,将恒定效率方案1(分段数为0)与变效率方案2(分段数为50)进行对比,恒定效率视为额定效率。两种方案下的热功率平衡如图7所示。在01:00—06:00、23:00—24:00时电价较低,方案1中CHP和电热泵一起供应热负荷,不考虑实际工程中CHP实际运行效率低于额定效率,达不到预期效果的情况;而方案2从电网购电供给电热泵产生热量作为主要供热设备,CHP不出力。在08:00—22:00时,电价较高,方案2电热泵出力减少,进而降低了购电成本;CHP以更高的热功率供给热负荷,并将一部分热储存到储热设备中。储热设备在13:00—16:00时以相较方案1更低的功率放热进行热负荷供给,减少了热功率损耗,预测情况更接近实际情况。
图7 热功率平衡图
两种方案下的电功率平衡如图8所示,风电利用如图9所示。在01:00—06:00、23:00—24:00时,电价较低,风电发电功率较大,方案1中CHP设备以额定效率运行供应一部分电负荷,浪费了部分风电,并且实际工程中CHP实际运行效率会低于额定效率,不能达到预期效果;而方案2中CHP不出力,增加了风电上网空间,电负荷完全由电网和风电供应。在08:00—22:00电价较高时风电减小,CHP出力增加,和电网、风电一起供电,预测情况更接近实际情况;方案2考虑CHP和储热设备的可变能量效率相比于方案1减小了风电出力偏差,弃风功率从3 314 kW降低到1 667 kW,风电消纳率增加了4.29%。
图8 电功率平衡图
图9 两种方案风电利用图
本文建立的考虑设备可变能量效率的EH能灵活地连接电热网,保证电热网的正常运行。与EH设备恒定能量效率的电热网模型相比,本文建立的考虑EH设备可变能量效率的电热网模型能更准确地描述电热网运行状态,有效提高电热网运行成本预测准确性,减小风电出力偏差。