CN112989612B - 一种基于麦考密克包络的电热综合能源系统线性建模方法 - Google Patents
一种基于麦考密克包络的电热综合能源系统线性建模方法 Download PDFInfo
- Publication number
- CN112989612B CN112989612B CN202110292002.2A CN202110292002A CN112989612B CN 112989612 B CN112989612 B CN 112989612B CN 202110292002 A CN202110292002 A CN 202110292002A CN 112989612 B CN112989612 B CN 112989612B
- Authority
- CN
- China
- Prior art keywords
- representing
- power
- output
- energy system
- envelope
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005485 electric heating Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 238000009826 distribution Methods 0.000 claims description 19
- 238000005259 measurement Methods 0.000 claims description 15
- 239000000243 solution Substances 0.000 claims description 15
- 238000010438 heat treatment Methods 0.000 claims description 14
- 230000008878 coupling Effects 0.000 claims description 13
- 238000010168 coupling process Methods 0.000 claims description 13
- 238000005859 coupling reaction Methods 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000001172 regenerating effect Effects 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 5
- 238000009833 condensation Methods 0.000 claims description 4
- 230000005494 condensation Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 24
- 238000012360 testing method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 239000000758 substrate Substances 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于麦考密克包络的电热综合能源系统线性建模方法,包括:建立电热综合能源系统;获取电热综合能源系统的拓扑结构和基本元件信息并建立基本元件模型;基于基本元件模型构建系统目标函数及约束条件,根据麦考密克包络的线性化策略建立电热综合能源系统线性化的联合优化调度模型,实现电热综合能源系统的高速求解,完成线性建模。本发明通过建立线性化的联合优化调度模型,热力系统采用变质量流量和变温度运行方式,给整个系统带来灵活性,同时模型通过麦考密克包络的线性化方法,保证了模型的高效求解;本发明可以用于给调度机构和系统规划提供参考,对于电热综合能源系统的经济运行具有重要意义。
Description
技术领域
本发明涉及多能源系统运行以及规划的技术领域,尤其涉及一种基于麦考密克包络的电热综合能源系统线性建模方法。
背景技术
在向绿色可持续社会转型的过程中,可再生能源发电在全球得到迅速发展。然而,在供暖期,大量热电联产机组采用“以热定电”的方式运行,限制了其自身的调峰能力,导致电力系统需要深度调峰时,整个电力系统的调峰能力不足,只能压缩风电上网的空间,造成了大量的弃风。高比例可再生能源的消纳成为亟待解决的问题。
鉴于电能和热能两种能源的紧密联系,从整体能源消耗的角度来说,热力系统与电力系统的联合运行,是解决高比例可再生能源消纳问题的关键;电热综合能源系统增加了可再生能源的渗透率,促进了新能源的消纳,可以实现多种能源的优势互补,提高能源的利用效率。
在以往的研究中,热力系统大多采用定质量流量定温度、变质量流量定温度、定质量流量变温度的运行方式,鲜有采用变质量流量和变温度的运行方式,变质量流量和变温度的运行方式使系统运行更加灵活;同时,也引入了由质量流量速率和温度的乘积带来的非线性项,使得潮流的计算更加复杂。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有存在的问题,提出了本发明。
因此,本发明解决的技术问题是:现有方案大多采用定质量流量定温度、变质量流量定温度、定质量流量变温度的运行方式,灵活性差,变质量流量和变温度在提升系统运行灵活性的同时而带来求解困难的问题。
为解决上述技术问题,本发明提供如下技术方案:基于传统配电网及运行方式为变质量流量和变温度的热力系统建立电热综合能源系统;获取所述电热综合能源系统的拓扑结构和基本元件信息并建立基本元件模型;基于所述基本元件模型构建系统目标函数及约束条件,根据麦考密克包络的线性化策略建立电热综合能源系统线性化的联合优化调度模型,实现所述电热综合能源系统的高速求解,完成线性建模。
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述电热综合能源系统的拓扑结构包括多个节点通过多条边连接构成,所述节点包括电力节点、热力节点、耦合节点,所述边包括电力线路、热力管道。
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述目标函数包括,
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述和/>的定义分别包括,
定义运行计量值为:
其中,h表示燃气机组的序号,表示燃气机组的输出热功率,/>表示第h个燃气机组的m阶计量值因子;
热电联产机组的运行计量值是关于它输出电功率和热功率的双变量函数,表示为:
其中,表示热电联产机组的m阶计量值因子;
电锅炉将电能转化为热能,产生一定数量的热蒸汽,热水等热媒,表示为:
其中,f表示电锅炉的序号,表示电锅炉的输出热功率,ηf表示电锅炉的输出效率,Γ表示电价;
燃气锅炉运行计量值是关于其输出热功率的二次函数,表示为:
其中,g表示燃气锅炉的序号,表示燃气锅炉的输出热功率,/>表示第g个燃气锅炉的m阶计量值因子;
分布式风电的弃风计量值与弃风量成正比,表示为:
其中,n表示分布式风电的序号,Wgn,t和分别表示分布式风电的实际输出电功率和预测发电功率,σ表示弃风的罚计量值。
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述约束条件包括热网管道约束、配电网约束、电源热源输出约束。
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述热网管道约束包括,
Φi=cmi(Ts,i-Tr,i)
Βf·Pl=0
其中,Φl表示热源/热负荷节点i的热功率,c表示管道工质的比热容,mi表示连接供热管网和回热管网的管道中的质量流量速率,Ts,i和Tr,i分别表示供热管网和回热管网的相对温度,Ts,i≥Tr,i恒成立,Tin和Tout分别表示管道的入口温度和出口温度,Tam表示环境温度,λ表示管道单位长度的热传导系数,L表示管道的长度,ml表示管道的质量流量速率,和/>分别表示其入口和出口与节点i相连接的管道集合,/>表示管道k的质量流量速率,Tk,out表示管道k的出口温度,Tj,in表示管道j的入口温度,Pl表示沿管道的压力损失,μ表示流量粘度,D表示管道的内径,/>和/>分别表示热源和热负荷的质量流量速率,Βf表示基本回路矩阵,可以写作Βf=[bhk],关联矩阵Α=[aik]定义如下:
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述配电网约束包括,
其中,Pij和Qij分别表示配电线路ij上传输的有功和无功功率,和/>分别表示节点j上由发电机组发出的有功和无功功率,/>和/>分别表示节点j上由负荷吸收的有功和无功功率,rij和xij分别表示配电线路ij上的电阻和电抗,Uj和U0分别表示节点j和平衡节点上的电压,π(j)表示节点注入功率来自于节点j的节点集合。
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:所述电源热源输出约束包括,燃气机组输出约束:
热电联产机组输出约束:
其中,e表示CHP机组的序号,cm,e表示CHP机组背压运行状态的热电转换效率,cv,e表示CHP机组运行在抽汽凝气运行状态的输出电功率,和/>分别表示CHP机组最大和最小输出电功率,/>表示CHP机组的最大输出热功率,/>和/>分别表示CHP机组的输出电功率和输出热功率;
燃气锅炉的输出约束为:
分布式风电的输出约束为:
电锅炉的输出约束为:
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:采用所述麦考密克包络的线性化策略将非线性约束线性化包括,根据式定义沿管道的热损失为:
由于0<λL/cml<<1.使用等价无穷小变换e-x=1-x,所述热损失近似为:
作为本发明所述的基于麦考密克包络的电热综合能源系统线性建模方法的一种优选方案,其中:采用所述麦考密克包络的线性化策略将非线性约束线性化进一步包括,基于所述麦考密克包络来保证凸性并且使得边界充分紧,则所述热管道模型中的Φi式定义为:
Φi=c(ws,i-wr,i)
ws,i≥mi Ts,i+mi Ts,i -miTs,i
wr,i≥mi Tr,i+mi Tr,i -miTr,i
同理,将式替换为线性约束。
本发明的有益效果:本发明通过基于麦考密克包络的电热综合能源系统线性化建模方法,建立了线性化的联合优化调度模型,热力系统采用变质量流量和变温度运行方式,给整个系统带来了灵活性,同时模型通过麦考密克包络的线性化方法,保证了模型的高效求解;本发明可以用于给调度机构和系统规划提供参考,对于电热综合能源系统的经济运行具有重要意义。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明一个实施例提供的一种基于麦考密克包络的电热综合能源系统线性建模方法的基本流程示意图;
图2为本发明一个实施例提供的一种基于麦考密克包络的电热综合能源系统线性建模方法的基于麦考密克包络策略的原理图;
图3为本发明一个实施例提供的一种基于麦考密克包络的电热综合能源系统线性建模方法的电热综合能源系统结构的示意图;
图4为本发明一个实施例提供的一种基于麦考密克包络的电热综合能源系统线性建模方法的实验结果对比图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
参照图1~2,为本发明的一个实施例,提供了一种基于麦考密克包络的电热综合能源系统线性建模方法,包括:
S1:基于传统配电网及运行方式为变质量流量和变温度的热力系统建立电热综合能源系统;需要说明的是,
电热综合能源系统中的热力系统采用变质量流量和变温度的运行方式,以往的热力系统大多采用定质量流量定温度、变质量流量定温度、定质量流量变温度的运行方式,鲜有采用变质量流量和变温度的运行方式,变质量流量和变温度的运行方式使系统运行更加灵活的同时,也引入了由质量流量速率和温度的乘积带来的非线性项,使得潮流的计算更加复杂,在本实施例中,质量流量速率、温度均为未知的代求量。
S2:获取电热综合能源系统的拓扑结构和基本元件信息并建立基本元件模型;需要说明的是,
电热综合能源系统拓扑结构由多个节点通过多条边连接构成,节点包括电力节点、热力节点、耦合节点三类,边包括电力线路、热力管道的两类;
具体的,电力节点之间或者电力节点和耦合节点之间通过电力线路连接或者不连接,热力节点之间或者热力节点和耦合节点之间通过天然气管道连接或者不连接;节点上有设备和负荷,设备包含非热电联产机组、热电联产机组、电锅炉、燃气锅炉、热泵、分布式风电六类,非热电联产机组位于电力节点/耦合节点上,热电联产机组、电锅炉、热泵仅存在于耦合节点上,燃气锅炉位于热力节点/耦合节点上;负荷包括电力负荷和热力负荷,电力负荷位于电力节点/耦合节点上,热力负荷位于热力节点/耦合节点上。
基本元件模型主要包括机组、风电、锅炉、网络等。
S3:基于基本元件模型构建系统目标函数及约束条件,根据麦考密克包络的线性化策略建立电热综合能源系统线性化的联合优化调度模型,实现电热综合能源系统的高速求解,完成线性建模;需要说明的是,
建立电热综合能源系统联合优化调度模型,其特征在于,模型主要包括以下几部分:
建立以下的目标函数:
其中,和/>的定义分别包括,
非热电联产机组-燃气机组:非热电联产机组在本模型中主要是传统的火电机组,本实施例以燃气机组为例,其运行计量值为:
其中,h表示燃气机组的序号,表示燃气机组的输出热功率,/>表示第h个燃气机组的m阶计量值因子;
热电联产机组:热电联产机组的运行计量值是关于它输出电功率和热功率的双变量函数,表示为:
其中,表示热电联产机组的m阶计量值因子;
电锅炉:电锅炉将电能转化为热能,产生一定数量的热蒸汽,热水等热媒,表示为:
其中,f表示电锅炉的序号,表示电锅炉的输出热功率,ηf表示电锅炉的输出效率,Γ表示电价;
燃气锅炉:燃气锅炉运行计量值是关于其输出热功率的二次函数,表示为:
其中,g表示燃气锅炉的序号,表示燃气锅炉的输出热功率,/>表示第g个燃气锅炉的m阶计量值因子;
分布式风电:分布式风电的弃风计量值与弃风量成正比,表示为:
其中,n表示分布式风电的序号,Wgn,t和分别表示分布式风电的实际输出电功率和预测发电功率,σ表示弃风的罚计量值。
同时建立如下约束条件:
热网管道约束包括,
Φi=cmi(Ts,i-Tr,i)
Βf·Pl=0
其中,Φl表示热源/热负荷节点i的热功率,c表示管道工质的比热容,mi表示连接供热管网和回热管网的管道中的质量流量速率,Ts,i和Tr,i分别表示供热管网和回热管网的相对温度,Ts,i≥Tr,i恒成立,Tin和Tout分别表示管道的入口温度和出口温度,Tam表示环境温度,λ表示管道单位长度的热传导系数,L表示管道的长度,ml表示管道的质量流量速率,和/>分别表示其入口和出口与节点i相连接的管道集合,/>表示管道k的质量流量速率,Tk,out表示管道k的出口温度,Tj,in表示管道j的入口温度,Pl表示沿管道的压力损失,μ表示流量粘度,D表示管道的内径,/>和/>分别表示热源和热负荷的质量流量速率,关联矩阵Α=[aik]定义如下:
Βf表示基本回路矩阵,可以写作Βf=[bhk],当bhk=1时表示管道k中的流速方向与回路h相一致;当bhk=-1时表示管道k中的流速方向与回路h相反;在区域供热网络中,质量流量是由循环水泵来维持的,鉴于循环水泵只消耗很少的一部分电量,可以假设循环水泵的费用可以忽略不计。
配电网约束包括,
其中,Pij和Qij分别表示配电线路ij上传输的有功和无功功率,和/>分别表示节点j上由发电机组发出的有功和无功功率,/>和/>分别表示节点j上由负荷吸收的有功和无功功率,rij和xij分别表示配电线路ij上的电阻和电抗,Uj和U0分别表示节点j和平衡节点上的电压,π(j)表示节点注入功率来自于节点j的节点集合。
电源热源输出约束包括,
燃气机组输出约束:
抽汽凝汽式机组从汽轮机中间级抽出一部分蒸汽,供给热负荷,实现在发电的同时还供热。抽汽凝汽式机组的电热耦合特性,热电联产机组输出约束:
其中,e表示CHP机组的序号,cm,e表示CHP机组背压运行状态的热电转换效率,cv,e表示CHP机组运行在抽汽凝气运行状态的输出电功率,和/>分别表示CHP机组最大和最小输出电功率,/>表示CHP机组的最大输出热功率,/>和/>分别表示CHP机组的输出电功率和输出热功率;
燃气锅炉的输出约束为:
分布式风电的输出约束为:
电锅炉的输出约束为:
进一步的,观察电热综合能源系统,目标函数是二次的,但约束条件中式Φi=cmi(Ts,i-Tr,i)和式中由质量流量和温度的乘积而引入了双线性函数项,式/>引入了指数项。因此,本发明采用麦考密克包络的线性化策略将非线性约束线性化。
具体包括,等价无穷小变换:根据式定义沿管道的热损失为:
在实践中,由于0<λL/cml<<1.使用等价无穷小变换e-x=1-x,热损失近似为:
式Φi=cmi(Ts,i-Tr,i)和中的双线性函数项是由于热力系统采用了变质量流量和变温度的运行方式引起的,可以通过麦考密克包络来保证凸性并且使得边界充分紧,使用麦考密克包络,热管道模型中的Φi式定义为:
Φi=c(ws,i-wr,i)
ws,i≥mi Ts,i+mi Ts,i -miTs,i
wr,i≥mi Tr,i+mi Tr,i -miTr,i
同理,将式替换为线性约束。电热综合能源系统联合优化模型被改造成了二次规划问题,可以利用现有的成熟求解器对其进行求解。
本发明方法通过基于麦考密克包络的电热综合能源系统线性化建模方法,建立了线性化的联合优化调度模型,热力系统采用变质量流量和变温度运行方式,给整个系统带来了灵活性,同时模型通过麦考密克包络的线性化方法,保证了模型的高效求解。
实施例2
参照图3~4为本发明另一个实施例,为对本方法中采用的技术效果加以验证说明,本实施例采用传统技术方案与本发明方法进行对比测试,以科学论证的手段对比试验结果,以验证本方法所具有的真实效果。
结合实施案例及附图对本发明方法作进一步说明:如图3所示,电热综合能源系统主要由热力系统与电力系统通过热电联产机组和电锅炉耦合而成,热电联产机组采用抽汽凝汽式机组,同时生产电能和热能,输出热功率与输出电功率有强耦合特性。电锅炉将电能转化为热能,产生一定数量的热蒸汽,热水等热媒,其输入端连入电力系统作为电负荷,电锅炉的输出端连入热力系统作为热源。热力系统主要由热源、热网和热负荷3部分组成,分别承担着热力生产、传输配送和使用的职能,其中,热网将热源与热负荷连接起来,并将热源产生的热量通过管道工质输送到热用户。在我国目前最常用的管道工质是热水,热网一般分为输热网和配热网,中间以换热站连接,换热站作为连接输热网和配热网的接口,起到交换热能的作用,在输热网中它相当于是热负荷,在配热网中它相当于是热源。
按照本发明内容的完整方法实施的实施例情况如下:处理获得某测试案例中由IEEE改造的电热综合能源测试系统。
首先,初始化该电热综合能源系统的参数。该电热综合能源测试结构示意图如图3所示,其中,浅色线表示热力网络,深色线表示电力网络,WG表示分布式风电,GT表示燃气机组,GB表示燃气锅炉,CHP表示热电联产机组,EB表示电锅炉。包含的电力系统和热力系统的原始数据来源于公开文献,具体可以在Liu,J.Wu,N.Jenkins,and A.Bagdanavicius,"Combined analysis of electricity and heat networks,"Applied Energy,vol.162,pp.1238-1250,2016.的附录A中查到。在电力系统的3和12节点分别设置两个1.5MVar和2MVar的静止无功补偿器来补偿无功功率和维持电压,同时设置两个案例来进行分析比较,一种是电热综合能源系统运行在传统的定质量流量和变温度运行方式下,另外一种是本模型中电热综合能源系统所采用的变质量流量和变温度的运行方式。
实施过程具体如下:
1、建立电热综合能源系统所包含的基本元件模型,其中电力系统是传统的配电网,热力系统采用变质量流量和变温度的运行方式;
2、其次,建立电热综合能源系统联合优化调度模型;
3、最后,采用基于麦考密克包络的线性化方法来保证高速求解;
4、将数值结果与常规的定质量流量和变温度结果进行比较。
根据本实施例提出的方法得到该系统采用变质量流量和变温度运行方式与传统的采用定质量流量变温度运行方式的风电消纳率对比图如图4所示,可以看到,与定质量流量和变温度相比,变质量流量和变温度运行方式为系统带来了灵活性,使得系统可以更好的消纳弃风,带来高额的经济收益,同时,线性化建模使得高效求解得到保证,整个求解速度约0.8秒,传统的基于IPOPT的非线性化求解方法在测试系统比本发明测试系统小的情况下,求解速度为96s,本发明方法相较于传统方法大大提升了求解速度,体现本发明方法的高效求解。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (8)
1.一种基于麦考密克包络的电热综合能源系统线性建模方法,其特性在于,包括:
基于传统配电网及运行方式为变质量流量和变温度的热力系统建立电热综合能源系统;
获取所述电热综合能源系统的拓扑结构和基本元件信息并建立基本元件模型;
基于所述基本元件模型构建系统目标函数及约束条件,根据麦考密克包络的线性化策略建立电热综合能源系统线性化的联合优化调度模型,实现所述电热综合能源系统的高速求解,完成线性建模;
采用所述麦考密克包络的线性化策略将非线性约束线性化包括,
根据式定义沿管道的热损失为:
由于0<λL/cml<<1.使用等价无穷小变换e-x=1-x,所述热损失近似为:
其中,Φi表示热源/热负荷节点i的热功率,c表示管道工质的比热容,ml表示管道的质量流量速率,Tin和Tout分别表示管道的入口温度和出口温度,Tam表示环境温度,λ表示管道单位长度的热传导系数,L表示管道的长度;
采用所述麦考密克包络的线性化策略将非线性约束线性化进一步包括,
基于所述麦考密克包络来保证凸性并且使得边界充分紧,则热管道模型中的Φi式定义为:
Φi=c(ws,i-wr,i)
ws,i≥mi Ts,i+mi Ts, i-miTs, i
wr,i≥mi Tr,i+mi Tr,i -miTr,i
其中,mi表示连接供热管网和回热管网的管道中的质量流量速率,m i和表示连接供热管网和回热管网的管道中的最小和最大质量流量速率,Ts,i和Tr,i分别表示供热管网和回热管网的相对温度,Ts,i≥Tr,i恒成立,Ts,i 和/>分别表示供热管网的最小和最大相对温度,Tr,i 和/>分别表示回热管网的最小和最大相对温度;
同理,将式替换为线性约束,其中,/>和/>分别表示其入口和出口与节点i相连接的管道集合,k表示管道k,d表示管道d。
2.如权利要求1所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述电热综合能源系统的拓扑结构包括多个节点通过多条边连接构成,所述节点包括电力节点、热力节点、耦合节点,所述边包括电力线路、热力管道。
3.如权利要求1或2所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述目标函数包括,
其中,表示为定义运行计量值,/>表示为热点联产机组的运行计量值,/>表示为电锅炉将电能转换成热能的值,/>表示为热锅炉运行计量值,/>表示为分布式风电的弃风计量值。
4.如权利要求3所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述和/>的定义分别包括,
定义运行计量值为:
其中,h表示燃气机组的序号,表示燃气机组的输出电功率,/>表示第h个燃气机组的m阶计量值因子;
热电联产机组的运行计量值是关于它输出电功率和热功率的双变量函数,表示为:
其中,表示热电联产机组的m阶计量值因子,/>和/>分别表示CHP机组的输出电功率和输出热功率,cv,e表示CHP机组运行在抽汽凝气运行状态的输出电功率;
电锅炉将电能转化为热能,产生的热媒,表示为:
其中,f表示电锅炉的序号,表示电锅炉的输出热功率,ηf表示电锅炉的输出效率,Γ表示电价;
燃气锅炉运行计量值是关于其输出热功率的二次函数,表示为:
其中,g表示燃气锅炉的序号,表示燃气锅炉的输出热功率,/>表示第g个燃气锅炉的m阶计量值因子;
分布式风电的弃风计量值与弃风量成正比,表示为:
其中,n表示分布式风电的序号,Wgn,t和分别表示分布式风电的实际输出电功率和预测发电功率,σ表示弃风的罚计量值。
5.如权利要求4所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述约束条件包括热网管道约束、配电网约束、电源热源输出约束。
6.如权利要求5所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述热网管道约束包括,
Φi=cmi(Ts,i-Tr,i)
Βf·Pl=0
其中,Φi表示热源/热负荷节点i的热功率,c表示管道工质的比热容,mi表示连接供热管网和回热管网的管道中的质量流量速率,m i和表示连接供热管网和回热管网的管道中的最小和最大质量流量速率,Ts,i和Tr,i分别表示供热管网和回热管网的相对温度,Ts,i≥Tr,i恒成立,Ts,i 和/>分别表示供热管网的最小和最大相对温度,Tr,i 和/>分别表示回热管网的最小和最大相对温度,Tin和Tout分别表示管道的入口温度和出口温度,Tin 和/>分别表示管道的最小和最大入口温度,Tout 和/>分别表示管道的最小和最大出口温度,Tam表示环境温度,λ表示管道单位长度的热传导系数,L表示管道的长度,ml表示管道的质量流量速率,ml 和/>分别表示管道的最小和最大质量流量速率,/>和/>分别表示其入口和出口与节点i相连接的管道集合,/>表示管道k的质量流量速率,Tk,out表示管道k的出口温度,Td,in表示管道d的入口温度,Pl表示沿管道的压力损失,μ表示流量粘度,D表示管道的内径,和/>分别表示热源和热负荷的质量流量速率,Βf表示基本回路矩阵,写作Βf=[bhk],Βf表示基本回路矩阵,可以写作Βf=[bhk],当bhk=1时表示管道k中的流速方向与回路h相一致;当bhk=-1时表示管道k中的流速方向与回路h相反;
关联矩阵Α=[aik]定义如下:
7.如权利要求6所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述配电网约束包括,
其中,Pij和Qij分别表示配电线路ij上传输的有功和无功功率,和/>分别表示节点j上由发电机组发出的有功和无功功率,/>和/>分别表示节点j上由负荷吸收的有功和无功功率,rij和xij分别表示配电线路ij上的电阻和电抗,Uj和U0分别表示节点j和平衡节点上的电压,π(j)表示节点注入功率来自于节点j的节点集合。
8.如权利要求7所述的基于麦考密克包络的电热综合能源系统线性建模方法,其特征在于:所述电源热源输出约束包括,
燃气机组输出约束:
其中,和/>表示燃气机组的最小和最大电功率;
热电联产机组输出约束:
其中,e表示CHP机组的序号,cm,e表示CHP机组背压运行状态的热电转换效率,cv,e表示CHP机组运行在抽汽凝气运行状态的输出电功率,和/>分别表示CHP机组最大和最小输出电功率,/>表示CHP机组的最大输出热功率,/>和/>分别表示CHP机组的输出电功率和输出热功率;
燃气锅炉的输出约束为:
其中,表示燃气锅炉的最大热功率;
分布式风电的输出约束为:
电锅炉的输出约束为:
其中,表示电锅炉的最大热功率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110292002.2A CN112989612B (zh) | 2021-03-18 | 2021-03-18 | 一种基于麦考密克包络的电热综合能源系统线性建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110292002.2A CN112989612B (zh) | 2021-03-18 | 2021-03-18 | 一种基于麦考密克包络的电热综合能源系统线性建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112989612A CN112989612A (zh) | 2021-06-18 |
CN112989612B true CN112989612B (zh) | 2023-09-19 |
Family
ID=76333428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110292002.2A Active CN112989612B (zh) | 2021-03-18 | 2021-03-18 | 一种基于麦考密克包络的电热综合能源系统线性建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112989612B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113836739B (zh) * | 2021-09-30 | 2023-06-02 | 沈阳工程学院 | 一种计及供热管网多损耗的电热联合系统潮流计算方法 |
CN115906411B (zh) * | 2022-10-24 | 2024-06-04 | 国网江苏省电力有限公司苏州供电分公司 | 考虑全动态的电热综合能源系统最优能流建模方法及系统 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101213679A (zh) * | 2005-06-28 | 2008-07-02 | Bsst有限责任公司 | 用于可变热功率源的热电发电机 |
CN106845701A (zh) * | 2017-01-11 | 2017-06-13 | 东南大学 | 一种基于热网和房屋热惯性的综合能源系统优化方法 |
CN106999060A (zh) * | 2014-08-11 | 2017-08-01 | 伊利诺伊大学评议会 | 用于分析温度特性和热传送特性的表皮器件 |
CN107067116A (zh) * | 2017-04-26 | 2017-08-18 | 燕山大学 | 一种多区域电热综合系统经济环境联合调度求解方法 |
CN109787259A (zh) * | 2019-01-23 | 2019-05-21 | 西安交通大学 | 一种基于新能源随机波动的多类型储能联合规划方法 |
CN109978625A (zh) * | 2019-03-28 | 2019-07-05 | 河海大学 | 一种计及电热气网络的综合能源系统多目标运行优化方法 |
CN110032828A (zh) * | 2019-05-10 | 2019-07-19 | 四川大学 | 一种计及需求响应柔性调节的配电网两阶段分布鲁棒dg优化配置线性方法 |
CN110263966A (zh) * | 2019-05-06 | 2019-09-20 | 天津大学 | 考虑动态传热过程的电-热综合能源系统优化调度方法 |
CN110348602A (zh) * | 2019-06-06 | 2019-10-18 | 国网浙江省电力有限公司经济技术研究院 | 计及天然气管网和热力管网特性的综合能源系统优化方法 |
CN111222257A (zh) * | 2020-02-27 | 2020-06-02 | 东南大学 | 基于凸优化的电热水多能流协同调度方法 |
CN111428351A (zh) * | 2020-03-11 | 2020-07-17 | 国网辽宁省电力有限公司大连供电公司 | 基于前推回代法的电-热综合能源系统潮流计算方法 |
AU2020100429A4 (en) * | 2020-03-20 | 2020-09-10 | Southeast University | A dynamic optimal energy flow computing method for the combined heat and power system |
CN111815068A (zh) * | 2020-07-22 | 2020-10-23 | 谢仕炜 | 二阶段约束边界紧缩算法求解城市综合能源网络优化方法 |
CN111950122A (zh) * | 2020-07-08 | 2020-11-17 | 国网(苏州)城市能源研究院有限责任公司 | 园区综合能源系统运行优化方法 |
-
2021
- 2021-03-18 CN CN202110292002.2A patent/CN112989612B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101213679A (zh) * | 2005-06-28 | 2008-07-02 | Bsst有限责任公司 | 用于可变热功率源的热电发电机 |
CN106999060A (zh) * | 2014-08-11 | 2017-08-01 | 伊利诺伊大学评议会 | 用于分析温度特性和热传送特性的表皮器件 |
CN106845701A (zh) * | 2017-01-11 | 2017-06-13 | 东南大学 | 一种基于热网和房屋热惯性的综合能源系统优化方法 |
CN107067116A (zh) * | 2017-04-26 | 2017-08-18 | 燕山大学 | 一种多区域电热综合系统经济环境联合调度求解方法 |
CN109787259A (zh) * | 2019-01-23 | 2019-05-21 | 西安交通大学 | 一种基于新能源随机波动的多类型储能联合规划方法 |
CN109978625A (zh) * | 2019-03-28 | 2019-07-05 | 河海大学 | 一种计及电热气网络的综合能源系统多目标运行优化方法 |
CN110263966A (zh) * | 2019-05-06 | 2019-09-20 | 天津大学 | 考虑动态传热过程的电-热综合能源系统优化调度方法 |
CN110032828A (zh) * | 2019-05-10 | 2019-07-19 | 四川大学 | 一种计及需求响应柔性调节的配电网两阶段分布鲁棒dg优化配置线性方法 |
CN110348602A (zh) * | 2019-06-06 | 2019-10-18 | 国网浙江省电力有限公司经济技术研究院 | 计及天然气管网和热力管网特性的综合能源系统优化方法 |
CN111222257A (zh) * | 2020-02-27 | 2020-06-02 | 东南大学 | 基于凸优化的电热水多能流协同调度方法 |
CN111428351A (zh) * | 2020-03-11 | 2020-07-17 | 国网辽宁省电力有限公司大连供电公司 | 基于前推回代法的电-热综合能源系统潮流计算方法 |
AU2020100429A4 (en) * | 2020-03-20 | 2020-09-10 | Southeast University | A dynamic optimal energy flow computing method for the combined heat and power system |
CN111950122A (zh) * | 2020-07-08 | 2020-11-17 | 国网(苏州)城市能源研究院有限责任公司 | 园区综合能源系统运行优化方法 |
CN111815068A (zh) * | 2020-07-22 | 2020-10-23 | 谢仕炜 | 二阶段约束边界紧缩算法求解城市综合能源网络优化方法 |
Non-Patent Citations (3)
Title |
---|
Co-optimization planning of integrated electricity and district heating systems based on improved quadratic convex relaxation;Xu Wang等;《Applied Energy》;20210122;第285卷;第1-18页 * |
含0-1变量的电力与能源系统优化调度方法研究;陈鸿琳;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20200115(第1期);第C042-217页 * |
考虑新能源发电不确定性的配电网重构策略研究;罗琳;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20200715(第7期);第C042-536页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112989612A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106447113B (zh) | 一种基于运行优化模型的多区域综合能源系统运行方法 | |
CN112989612B (zh) | 一种基于麦考密克包络的电热综合能源系统线性建模方法 | |
CN109492325B (zh) | 一种基于扩展能源集线器的多能耦合系统潮流分析方法 | |
CN110263966A (zh) | 考虑动态传热过程的电-热综合能源系统优化调度方法 | |
CN110502791B (zh) | 基于能源集线器的综合能源系统稳态建模方法 | |
CN106844895B (zh) | 一种冷热电联供微网能量流的解耦计算方法 | |
CN113792969A (zh) | 考虑气网动态特性和电转气综合能源系统优化调度方法 | |
CN109359839A (zh) | 一种综合能源系统中区域热网的节点热价计算方法 | |
CN111428962B (zh) | 一种基于热网动态特性端口等值模型的电热协同调度方法 | |
CN114077934A (zh) | 一种综合能源微网互联系统及其调度方法 | |
CN110298556B (zh) | 一种基于能量价值的多能协同园区能源利用效率控制方法 | |
CN113255105B (zh) | 双向耦合网络结构的电气热综合能源系统潮流计算方法 | |
Man et al. | State estimation for integrated energy system containing electricity, heat and gas | |
Tang et al. | Optimal dispatch of integrated energy campus microgrids considering the time-delay of pipelines | |
CN113761727B (zh) | 一种含分布式电热泵的热电联合系统优化调度模型构建方法 | |
CN112906220B (zh) | 综合能源微网园区系统状态的估计方法 | |
Zhao et al. | Optimal scheduling method for electrical-thermal integrated energy system considering heat storage characteristics of heating network | |
CN112421635B (zh) | 一种基于多能微网群网内网外多能流潮流计算的方法 | |
KR101515374B1 (ko) | 열병합발전 시스템으로 구축된 개별 열/전기 에너지 및 중앙 열/전기 에너지 공급 연계 시스템, 및 이를 운용하는 방법 | |
CN112348343B (zh) | 一种考虑不确定性的多能流配电网运行成本评估方法 | |
Zhao et al. | A multi-objective power flow optimization model of electric and thermal distribution network using the power flow method | |
Chen et al. | The application of energy network theory in the analysis of district electricity and heating system | |
CN115860325A (zh) | 一种基于效能转换模型的多元电气互联网络分析方法 | |
CN111969588B (zh) | 一种多能流系统最优能流的求解方法、存储介质及设备 | |
Chen et al. | A joint solution framework for the optimal operation and evaluation of integrated energy systems based on exergy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |