CN114781922A - 一种电热综合能源系统优化调度方法 - Google Patents
一种电热综合能源系统优化调度方法 Download PDFInfo
- Publication number
- CN114781922A CN114781922A CN202210553323.8A CN202210553323A CN114781922A CN 114781922 A CN114781922 A CN 114781922A CN 202210553323 A CN202210553323 A CN 202210553323A CN 114781922 A CN114781922 A CN 114781922A
- Authority
- CN
- China
- Prior art keywords
- heat
- model
- load
- time period
- power
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005485 electric heating Methods 0.000 title claims abstract description 30
- 230000004044 response Effects 0.000 claims abstract description 27
- 239000000446 fuel Substances 0.000 claims abstract description 20
- 230000005611 electricity Effects 0.000 claims abstract description 17
- 238000012423 maintenance Methods 0.000 claims abstract description 15
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 238000003860 storage Methods 0.000 claims description 44
- 238000005338 heat storage Methods 0.000 claims description 42
- 239000002028 Biomass Substances 0.000 claims description 31
- 238000005457 optimization Methods 0.000 claims description 26
- 238000004146 energy storage Methods 0.000 claims description 14
- 238000009825 accumulation Methods 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000010977 unit operation Methods 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 230000006870 function Effects 0.000 description 14
- 238000010248 power generation Methods 0.000 description 4
- 125000000217 alkyl group Chemical group 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种电热综合能源系统优化调度方法,涉及综合能源系统优化调度技术领域,包括以下步骤:采集电热综合能源系统的运行数据和系统设备参数;建立太阳能集热系统出力模型;建立电负荷需求侧响应模型;建立系统设备运行模型;设定电负荷、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束;建立综合考虑购电费用、燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数;对非线性约束进行线性化变换;建立日前区间电热综合能源系统优化调度模型;基于Pyomo搭建区间最优、劣子模型,调用求解器GLPK进行求解。本发明利用Pyomo‑GLPK的组合对模型进行建模和求解,可以有效的降低系统的运行成本且拥有很好的求解速度。
Description
技术领域
本发明涉及综合能源系统优化调度技术领域,具体为一种电热综合能源系统优化调度方法。
背景技术
在区域综合能源系统内,风电出力及太阳能集热受天气因素的影响,其出力预测具有较强的不确定性;负荷预测也存在预测误差的不确定性。在实际系统决策中,获得精确的概率密度函数往往较困难,然而获得不确定变量的取值范围则相对容易,且所需要的信息也大大减少。
目前对太阳能集热系统的研究大部分为一个独立系统供给热负荷,其并未与其它设备相耦合参与到调度中。GLPK作为一个开源的求解器,被用于能源调度求解的例子较少。
发明内容
本发明针对上述现有技术存在的不足,提供一种电热综合能源系统优化调度方法。
为了解决上述技术问题,本发明提供了如下的技术方案:
一种电热综合能源系统优化调度方法,包括以下步骤:
采集电热综合能源系统的运行数据和系统设备参数;
建立太阳能集热系统出力模型;
建立电负荷需求侧响应模型;
建立系统设备运行模型;
设定电负荷、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束;
建立综合考虑购电费用、燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数;
对非线性约束进行线性化变换;
建立日前区间电热综合能源系统优化调度模型;
基于Pyomo搭建区间最优、劣子模型,调用求解器GLPK进行求解。
作为本发明的进一步技术方案为,所述电热综合能源系统包括电力系统和热力系统;所述系统的运行数据包括电负荷、热负荷功率预测,新能源出力的负荷预测,分时电价;所述系统设备参数包括出力设备的额定功率和转换效率,储能设备的放能、储能功率、效率和损耗率;所述电力系统包括风电场、蓄电池、蓄热式电锅炉和地源热泵,蓄热式电锅炉和地源热泵同时属于热力系统;所述热力系统包括蓄热式电锅炉、地源热泵、生物质锅炉和太阳能集热系统。
作为本发明的进一步技术方案为,所述建立太阳能集热系统出力模型:
V=75L/m2*As (2)
Qs,t=JtAsηd (3)
式中:As为集热器总面积;Qs为一天所需的热负荷;V为储热水箱的容积;JT为当地日平均太阳辐照(MJ/(m2·d));ηcd为采暖集热器平均集热效率;ηL为管路及储热装置热损失率;f为太阳能保证率;Qs,t为t时刻太阳能集热系统产热量;Jt为当地t时间段内的平均太阳辐照量;ηd为使用期集热器平均效率;Ss,t和Ss,t+1为t时间段与下一时间段内的蓄热水箱的蓄热存量;ηsr,in和ηsr,out分别为储、放热效率。
作为本发明的进一步技术方案为,所述建立电负荷需求侧响应模型包括可消减负荷模型和可平移负荷模型;其中,
可削减负荷模型为:
可平移负荷模型为:
Pt shift=mtPshift,t (6)
式中:Pshift,t为时段t的可平移负荷功率;Pshift,t为可平移负荷的额定功率。
作为本发明的进一步技术方案为,所述建立设备运行模型包括蓄电池模型、蓄热式锅炉模型、地源热泵模型和生物质锅炉模型;其中,
蓄电池模型为:
式中:PBattery,t为蓄电池t时间段的出力;和分别为蓄电池t 时间段的输入、输出功率;SBattery,t和SBattery,t+1分别为t时间段与下一时间段内的蓄电量;ηB,in和ηB,out分别为储、放电效率;
蓄热式电锅炉模型为:
QREB,t=ηREBPREB,t (9)
式中:QREB,t为电锅炉在t时间段的放热量;ηREB为蓄热式电锅炉的转换效率;PREB,t为其在t时间段的电功率;ηloss1为蓄热部分的热损失系数;和分别为蓄热式电锅炉的蓄热体在t时间段的蓄热量和放热量;SREB,t 和SREB,t+1分别为t时间段和t+1时间段内的蓄热存量;ηr,in和ηr,out分别为储热、放热效率;
地源热泵模型为:
QHp,t=coph×PHp,t (11)
式中:QHp,t为地源热泵在t时间段的产热功率;coph为地源热泵的制热能效比;PHp,t为其在t时间段的用电功率;
生物质锅炉模型为:
QBb,t=ηBbβBbWBb,t (12)
式中:QBb,t为生物质锅炉在t时间段的产热功率;ηBb为生物质锅炉的热效率;βBb为生物质固化燃料发热值;WBb,t为生物质锅炉在t时间段内的使用燃料重量。
作为本发明的进一步技术方案为,所述设定电负荷出力平衡约束、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束:
电负荷出力平衡约束为:
热负荷出力平衡约束为:
式中:为太阳能集热系统在t时间段的放热量;QREB,t为电锅炉在t 时间段的放热量;和分别为蓄热式电锅炉的蓄热体在t时间段的蓄热量和放热量;QHp,t为地源热泵在t时间段的放热量;Qper,t为用户在t时间段需求的热负荷;Ain,t,Aout,t分别为表示蓄热式电锅炉蓄热部分蓄热、放热状态的0-1变量;
蓄电池出力约束为:
蓄电池储能约束为:
SBattery,t(0)=SBattery,t(T) (17)
蓄热式电锅炉出力约束:
蓄热式电锅炉储能约束:
SREB,t(0)=SREB,t(T) (21)
太阳能集热系统出力约束:
太阳能集热系统储热约束:
Ss,min≤Ss,t≤Ss,max (24)
地源热泵出力约束:
生物质锅炉出力约束:
需求侧响应约束:
作为本发明的进一步技术方案为,所述对非线性约束进行线性化变换,具体为:对(0,1)变量与正数变量相乘的非线性约束进行线性化变换:
X(t)为(0,1)变量,V(t)为正数变量,先假定X(t)一个较大的上限值设为Xmax,增加两个临时变量Y(t),Z(t);线性化的步骤如下:
加入等式约束:
Y(t)=X(t)-Z(t) (30)
加入不等式约束:
这样Y(t)完全等价于V(t)X(t)。
作为本发明的进一步技术方案为,所述建立综合考虑购电费用、物质燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数:
式中:F为系统的总运行成本;Cele为购电费用;Cb为购买生物质燃料的费用;Cwind为弃风惩罚费用;Cyw为能源设备的运维成本;Cbc为柔性负荷补偿成本;Ccut为削减负荷补偿费用;Cshift为平移负荷补偿费用;λele为分时电价,PBuy,t为t时间段系统从电网购买的电量;λBb为购买生物质燃料单位质量价格;λwind为风电机组的弃风惩罚系数,为风电机组在t时间段的预测产电功率;Pwind,t为风电机组在t时间段的实际用电功率;λj,yw为能源设备j的单位运维成本;Pj,t为t时间段设备j的出力;T为调度总时间段;为协议中单位功率电负荷补偿价格;为单位功率负荷平移的补偿价格。
作为本发明的进一步技术方案为,所述建立日前区间电热综合能源系统优化调度模型,具体为:
区间线性规划模型一般形式:
将构建的模型一般形式分解为最优子模型和最劣子模型,通过两阶段分解法分别求得最优值,从而得到最优值取值区间;其中,两阶段分解法包括下限优化求解和上限优化求解,
所述下限优化求解具体为:
目标函数:
约束条件:
上限优化求解具体为:
目标函数:
约束条件:
作为本发明的进一步技术方案为,所述基于Pyomo构建区间最优、劣子模型,调用求解器GLPK进行求解;具体为:
根据Pyomo框架搭建日前区间电热综合能源系统优化调度模型,将系统的运行数据、设备参数代入上述模型使用GLPK进行求解。
本发明的有益效果是:
1、本发明提供一种综合考虑新能源出力不确定性、需求侧响应和储能的日前区间优化调度模型,将太阳能集热系统与系统相耦合进行调度,并利用Pyomo-GLPK的组合对模型进行建模和求解,可以有效的降低系统的运行成本且拥有很好的求解速度。
2、本发明综合考虑需求侧响应和新能源地不确定性,用区间数学进行表示;
3、本发明将太阳能集热系统与其它设备进行耦合对其进行日前调度;
4、本发明模型的搭建和求解使用Pyomo+GLPK的组合,Pyomo一种独立于用于解决问题的底层求解器的建模语言,GLPK是一种开源的线性规划求解器,使得其在工程中应用具有一定的优势。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。
图1是本发明电-热综合能源系统优化调度方法的流程图;
图2是本发明电-热综合能源系统优化调度方法的流程图;
图3是本发明实施例中电-热综合能源系统的结构图;
图4是本发明实施例中电价曲线图;
图5是本发明实施例中电负荷、热负荷、风机发电量和太阳能集热系统集热量预测曲线图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
本发明实施例中,区间线性规划是一种很有效的处理不确定性问题的方法,它是将区间数的理论和方法应用于线性规划中,并且在目标函数或约束条件中含有区间数的一类线性规划。
Pyomo为基于python的开源软件包,用于制定和分析优化模型,GLPK 为一个开源的线性规划求解器,二者相组合可以很好的在工程中运行
参见图1,本发明提供一种电热综合能源系统优化调度方法,包括以下步骤:
步骤1:采集电热综合能源系统的运行数据和系统设备参数;
步骤2:建立太阳能集热系统出力模型;
步骤3:建立电负荷需求侧响应模型;
步骤4:建立系统设备运行模型;
步骤5:设定电负荷、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束;
步骤6:建立综合考虑购电费用、燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数;
步骤7:对非线性约束进行线性化变换;
步骤8:建立日前区间电热综合能源系统优化调度模型;
步骤9:基于Pyomo搭建区间最优、劣子模型,调用求解器GLPK进行求解。
下面通过具体案例,并结合附图,对本发明详细说明。
选取某地的设计算例数据具体分析,算例结构如图2所示,设备包含风电机组、太阳能集热系统、蓄电池、蓄热式电锅炉、地源热泵。以24小时为一个调度周期,单位调度时间为1h,对最优解与最差解模型分别利用 Pyomo-GLPK进行求解。
表1出力设备参数
表2储能设备参数
步骤1,采集综合能源系统的运行数据和系统设备参数;
采集运行数据包括电负荷、热负荷功率预测,新能源出力的负荷预测,分时电价;以及系统设备参数包括出力设备的额定功率和转换效率,储能设备的放能、储能功率、效率和损耗率。
步骤2,建立太阳能集热系统出力模型;
V=75L/m2*As (2)
Qs,t=JtAsηd (3)
式中:As为集热器总面积;Qs为一天所需的热负荷;V为储热水箱的容积;JT为当地日平均太阳辐照(MJ/(m2·d));ηcd为采暖集热器平均集热效率;ηL为管路及储热装置热损失率;f为太阳能保证率;Qs,t为t时刻太阳能集热系统产热量;Jt为当地t时间段内的平均太阳辐照量;ηd为使用期集热器平均效率;Ss,t和Ss,t+1为t时间段与下一时间段内的蓄热水箱的蓄热存量;ηsr,in和ηsr,out分别为储、放热效率。
步骤3,建立电负荷需求侧响应模型:
a.可削减负荷模型:
b.可平移负荷模型:
Pt shift=mtPshift,t (6)
式中:Pshift,t为时段t的可平移负荷功率;Pshift,t为可平移负荷的额定功率。
步骤4,建立设备运行模型;
a.蓄电池模型:
式中:PBattery,t为蓄电池t时间段的出力;和分别为蓄电池t时间段的输入、输出功率;SBattery,t和SBattery,t+1分别为t时间段与下一时间段内的蓄电量;ηB,in和ηB,out分别为储、放电效率。
b.蓄热式电锅炉模型:
QREB,t=ηREBPREB,t (9)
式中:QREB,t为电锅炉在t时间段的放热量;ηREB为蓄热式电锅炉的转换效率;PREB,t为其在t时间段的电功率;ηloss1为蓄热部分的热损失系数;和分别为蓄热式电锅炉的蓄热体在t时间段的蓄热量和放热量;SREB,t和SREB,t+1分别为t时间段和t+1时间段内的蓄热存量;ηr,in和ηr,out分别为储热、放热效率。
c.地源热泵模型:
QHp,t=coph×PHp,t (11)
式中:QHp,t为地源热泵在t时间段的产热功率;coph为地源热泵的制热能效比;PHp,t为其在t时间段的用电功率。
d.生物质锅炉模型:
QBb,t=ηBbβBbWBb,t (12)
式中:QBb,t为生物质锅炉在t时间段的产热功率;ηBb为生物质锅炉的热效率,取0.8;βBb为生物质固化燃料发热值,取5.4(kw/kg);WBb,t为生物质锅炉在t时间段内的使用燃料重量,燃料价格0.7元/kg。
步骤5,设定电负荷、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束;
a.电负荷出力平衡约束:
b.热负荷出力平衡约束:
式中:为太阳能集热系统在t时间段的放热量;QREB,t为电锅炉在t 时间段的放热量;和分别为蓄热式电锅炉的蓄热体在t时间段的蓄热量和放热量;QHp,t为地源热泵在t时间段的放热量;Qper,t为用户在t时间段需求的热负荷;Ain,t,Aout,t分别为表示蓄热式电锅炉蓄热部分蓄热、放热状态的0-1变量。
c.蓄电池出力约束:
d.蓄电池储能约束:
SBattery,t(0)=SBattery,t(T) (17)
e.蓄热式电锅炉出力约束:
f.蓄热式电锅炉储能约束:
SREB,t(0)=SREB,t(T) (21)
g.太阳能集热系统出力约束:
h.太阳能集热系统储热约束:
Ss,min≤Ss,t≤Ss,max (24)
i.地源热泵出力约束:
j.生物质锅炉出力约束:
k.需求侧响应约束:
优选地,步骤6中,对(0,1)变量与正数变量相乘的非线性约束进行线性化变换;
X(t)为(0,1)变量,V(t)为正数变量,先假定X(t)一个较大的上限值设为Xmax。增加两个临时变量Y(t),Z(t)。线性化的步骤如下:
加入等式约束:
Y(t)=X(t)-Z(t) (30)
加入不等式约束:
这样Y(t)完全等价于V(t)X(t)。
步骤7,建立综合考虑购电费用、物质燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数:
式中:F为系统的总运行成本;Cele为购电费用;Cb为购买生物质燃料的费用;Cwind为弃风惩罚费用;Cyw为能源设备的运维成本;Cbc为柔性负荷补偿成本;Ccut为削减负荷补偿费用;Cshift为平移负荷补偿费用;λele为分时电价,PBuy,t为t时间段系统从电网购买的电量;λBb为购买生物质燃料单位质量价格;λwind为风电机组的弃风惩罚系数,为风电机组在t时间段的预测产电功率;Pwind,t为风电机组在t时间段的实际用电功率;λj,yw为能源设备j的单位运维成本;Pj,t为t时间段设备j的出力;T为调度总时间段;为协议中单位功率电负荷补偿价格,取0.2元/kW·h;为单位功率负荷平移的补偿价格,取0.05元/kW·h。参见图3和图4。
优选地,步骤8中,建立日前区间电-热综合能源系统优化调度模型;
区间线性规划模型一般形式:
将构建的模型一般形式分解为最优子模型和最劣子模型,通过两阶段分解法分别求得最优值,从而得到最优值取值区间。
两阶段分解法:
1)下限优化求解
目标函数:
约束条件:
2)上限优化求解
目标函数:
约束条件:
优选地,步骤9中,基于Pyomo搭建区间最优、劣子模型,调用求解器GLPK进行求解;
根据Pyomo框架搭建日前区间电-热综合能源系统优化调度模型,将系统的运行数据、设备参数代入上述模型使用GLPK进行求解,将得到的最优调度结果进行分析。
本实施例选取如下三种场景进行对比分析。
场景1:太阳能集热系统与其它设备独立进行供能,考虑风电、太阳能集热的不确定性为5%,运用Pyomo+GLPK进行求解。
场景2:太阳能集热系统与其它设备独立进行供能,考虑风电、太阳能集热的不确定性为5%,运用Yalmip+CPLEX进行求解。
场景3:太阳能集热系统可由地源热泵辅助储热,考虑风电、太阳能集热的不确定性为5%,运用Pyomo+GLPK进行求解。
表3为调度结果分析,由表3对比结果可以看出场景3可以有效的降低成本和提高对风电的消纳能力。而场景1中运维费用较低,原因在于场景3中增加了对地源热泵和太阳能集热系统的使用且风电利用率高,因而设备的运维费用较高。场景3中太阳能集热系统和其它供能设备耦合使用,虽然较少的增加了设备的使用频率,但增加了风电的消纳能力且降低了系统的运行成本。并且场景1和场景2对比可知,本文采用的求解方式 Pyomo-GLPK,其求解精度和求解速度和常用的求解器CPLEX相差无几。
本发明实施案例综合考虑新能源的不确定性和需求侧响应进行电-热综合能源系统的建模,将太阳能集热系统与其它供能设备相耦合进行日前调度,通过Pyomo-GLPK进行建模并求解计算,对比分析三种调度结果,得到下面结论:
将太阳能集热系统和其它供能设备相互耦合进行日前调度可以在提高新能源的消纳能力的同时有效降低系统的运行成本,且GLPK的求解精度和求解速度不差于CPLEX的同时GLPK因其为开源求解器,使用相对自由,有利于其在工程上的运用。
综上所述,将太阳能集热系统与其它供能设备相耦合进行日前调度可有效降低系统地运行成本和增加新能源的消纳能力,Pyomo-GLPK的组合使用有较好的求解速度和求解结果,更加适用于工程使用。
表3 三种场景下的运行情况
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种电热综合能源系统优化调度方法,其特征在于,包括以下步骤:
采集电热综合能源系统的运行数据和系统设备参数;
建立太阳能集热系统出力模型;
建立电负荷需求侧响应模型;
建立系统设备运行模型;
设定电负荷、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束;
建立综合考虑购电费用、燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数;
对非线性约束进行线性化变换;
建立日前区间电热综合能源系统优化调度模型;
基于Pyomo搭建区间最优、劣子模型,调用求解器GLPK进行求解。
2.根据权利要求1所述的一种电热综合能源系统优化调度方法,其特征在于:
所述电热综合能源系统包括电力系统和热力系统;
所述系统的运行数据包括电负荷、热负荷功率预测,新能源出力的负荷预测,分时电价;
所述系统设备参数包括出力设备的额定功率和转换效率,储能设备的放能、储能功率、效率和损耗率;
所述电力系统包括风电场、蓄电池、蓄热式电锅炉和地源热泵,蓄热式电锅炉和地源热泵同时属于热力系统;
所述热力系统包括蓄热式电锅炉、地源热泵、生物质锅炉和太阳能集热系统。
5.根据权利要求1所述的一种电热综合能源系统优化调度方法,其特征在于:
所述建立设备运行模型包括蓄电池模型、蓄热式锅炉模型、地源热泵模型和生物质锅炉模型;其中,
蓄电池模型为:
式中:PBattery,t为蓄电池t时间段的出力;和分别为蓄电池t时间段的输入、输出功率;SBattery,t和SBattery,t+1分别为t时间段与下一时间段内的蓄电量;ηB,in和ηB,out分别为储、放电效率;
蓄热式电锅炉模型为:
QREB,t=ηREBPREB,t (9)
式中:QREB,t为电锅炉在t时间段的放热量;ηREB为蓄热式电锅炉的转换效率;PREB,t为其在t时间段的电功率;ηloss1为蓄热部分的热损失系数;和分别为蓄热式电锅炉的蓄热体在t时间段的蓄热量和放热量;SREB,t和SREB,t+1分别为t时间段和t+1时间段内的蓄热存量;ηr,in和ηr,out分别为储热、放热效率;
地源热泵模型为:
QHp,t=coph×PHp,t (11)
式中:QHp,t为地源热泵在t时间段的产热功率;coph为地源热泵的制热能效比;PHp,t为其在t时间段的用电功率;
生物质锅炉模型为:
QBb,t=ηBbβBbWBb,t (12)
式中:QBb,t为生物质锅炉在t时间段的产热功率;ηBb为生物质锅炉的热效率;βBb为生物质固化燃料发热值;WBb,t为生物质锅炉在t时间段内的使用燃料重量。
6.根据权利要求1所述的一种电热综合能源系统优化调度方法,其特征在于:
所述设定电负荷出力平衡约束、热负荷出力平衡约束、设备出力约束条件和需求侧响应约束:
电负荷出力平衡约束为:
热负荷出力平衡约束为:
式中:为太阳能集热系统在t时间段的放热量;QREB,t为电锅炉在t时间段的放热量;和分别为蓄热式电锅炉的蓄热体在t时间段的蓄热量和放热量;QHp,t为地源热泵在t时间段的放热量;Qper,t为用户在t时间段需求的热负荷;Ain,t,Aout,t分别为表示蓄热式电锅炉蓄热部分蓄热、放热状态的0-1变量;
蓄电池出力约束为:
蓄电池储能约束为:
SBattery,t(0)=SBattery,t(T) (17)
蓄热式电锅炉出力约束:
蓄热式电锅炉储能约束:
SREB,t(0)=SREB,t(T) (21)
太阳能集热系统出力约束:
太阳能集热系统储热约束:
Ss,min≤Ss,t≤Ss,max (24)
地源热泵出力约束:
生物质锅炉出力约束:
需求侧响应约束:
8.根据权利要求1所述的一种电热综合能源系统优化调度方法,其特征在于:
所述建立综合考虑购电费用、物质燃料购买费用、弃风惩罚费用、设备运维费用和需求侧响应补偿费用的目标函数:
9.根据权利要求1所述的一种电热综合能源系统优化调度方法,其特征在于:所述建立日前区间电热综合能源系统优化调度模型,具体为:
区间线性规划模型一般形式:
将构建的模型一般形式分解为最优子模型和最劣子模型,通过两阶段分解法分别求得最优值,从而得到最优值取值区间;其中,两阶段分解法包括下限优化求解和上限优化求解,
所述下限优化求解具体为:
目标函数:
约束条件:
上限优化求解具体为:
目标函数:
约束条件:
10.根据权利要求1所述的一种电热综合能源系统优化调度方法,其特征在于:所述基于Pyomo构建区间最优、劣子模型,调用求解器GLPK进行求解;具体为:
根据Pyomo框架搭建日前区间电热综合能源系统优化调度模型,将系统的运行数据、设备参数代入上述模型使用GLPK进行求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210553323.8A CN114781922B (zh) | 2022-05-19 | 2022-05-19 | 一种电热综合能源系统优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210553323.8A CN114781922B (zh) | 2022-05-19 | 2022-05-19 | 一种电热综合能源系统优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114781922A true CN114781922A (zh) | 2022-07-22 |
CN114781922B CN114781922B (zh) | 2023-05-09 |
Family
ID=82408516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210553323.8A Active CN114781922B (zh) | 2022-05-19 | 2022-05-19 | 一种电热综合能源系统优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114781922B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013060082A1 (zh) * | 2011-10-23 | 2013-05-02 | 西安交通大学 | 燃气联合循环与太阳能发电联合制热系统及其调度方法 |
JP2014155286A (ja) * | 2013-02-06 | 2014-08-25 | Taisei Corp | 設備運転計画決定方法 |
CN104734168A (zh) * | 2015-03-13 | 2015-06-24 | 山东大学 | 一种基于电热联合调度的微电网运行优化系统及方法 |
CN108932560A (zh) * | 2018-06-13 | 2018-12-04 | 天津大学 | 基于模型预测控制的园区综合能源系统优化调度方法 |
CN109919480A (zh) * | 2019-03-01 | 2019-06-21 | 广东电网有限责任公司 | 一种三层目标的能源互联网规划方法及设备 |
CN110009122A (zh) * | 2018-12-27 | 2019-07-12 | 国网北京市电力公司 | 户用多能互补综合能源利用系统能量优化调度方法及系统 |
CN110138004A (zh) * | 2019-05-08 | 2019-08-16 | 上海电力学院 | 一种多能互补系统优化运行方法 |
CN112069734A (zh) * | 2020-09-04 | 2020-12-11 | 华北电力大学 | 一种电热泵-热电联合系统的优化调控方法及系统 |
CN112685879A (zh) * | 2020-12-16 | 2021-04-20 | 南京国电南自电网自动化有限公司 | 一种区域电热互联能源系统的多目标优化方法 |
WO2021249665A1 (en) * | 2020-06-10 | 2021-12-16 | Eaton Intelligent Power Limited | Method and system for resource management |
CN114123273A (zh) * | 2021-11-12 | 2022-03-01 | 青海综合能源服务有限公司 | 一种风电-光伏-储能联合系统的控制方法及系统 |
US20220074620A1 (en) * | 2020-09-04 | 2022-03-10 | North China Electric Power University | Optimized regulating and controlling method and system for integrated electricity and heat system with heat pumps |
CN114508784A (zh) * | 2022-02-16 | 2022-05-17 | 重庆邮电大学 | 一种多源互补供热系统及其优化控制方法 |
-
2022
- 2022-05-19 CN CN202210553323.8A patent/CN114781922B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013060082A1 (zh) * | 2011-10-23 | 2013-05-02 | 西安交通大学 | 燃气联合循环与太阳能发电联合制热系统及其调度方法 |
JP2014155286A (ja) * | 2013-02-06 | 2014-08-25 | Taisei Corp | 設備運転計画決定方法 |
CN104734168A (zh) * | 2015-03-13 | 2015-06-24 | 山东大学 | 一种基于电热联合调度的微电网运行优化系统及方法 |
CN108932560A (zh) * | 2018-06-13 | 2018-12-04 | 天津大学 | 基于模型预测控制的园区综合能源系统优化调度方法 |
CN110009122A (zh) * | 2018-12-27 | 2019-07-12 | 国网北京市电力公司 | 户用多能互补综合能源利用系统能量优化调度方法及系统 |
CN109919480A (zh) * | 2019-03-01 | 2019-06-21 | 广东电网有限责任公司 | 一种三层目标的能源互联网规划方法及设备 |
CN110138004A (zh) * | 2019-05-08 | 2019-08-16 | 上海电力学院 | 一种多能互补系统优化运行方法 |
WO2021249665A1 (en) * | 2020-06-10 | 2021-12-16 | Eaton Intelligent Power Limited | Method and system for resource management |
CN112069734A (zh) * | 2020-09-04 | 2020-12-11 | 华北电力大学 | 一种电热泵-热电联合系统的优化调控方法及系统 |
US20220074620A1 (en) * | 2020-09-04 | 2022-03-10 | North China Electric Power University | Optimized regulating and controlling method and system for integrated electricity and heat system with heat pumps |
WO2022048127A1 (zh) * | 2020-09-04 | 2022-03-10 | 华北电力大学 | 一种电热泵-热电联合系统的优化调控方法及系统 |
CN112685879A (zh) * | 2020-12-16 | 2021-04-20 | 南京国电南自电网自动化有限公司 | 一种区域电热互联能源系统的多目标优化方法 |
CN114123273A (zh) * | 2021-11-12 | 2022-03-01 | 青海综合能源服务有限公司 | 一种风电-光伏-储能联合系统的控制方法及系统 |
CN114508784A (zh) * | 2022-02-16 | 2022-05-17 | 重庆邮电大学 | 一种多源互补供热系统及其优化控制方法 |
Non-Patent Citations (6)
Title |
---|
FENG ZHOU等: "Enhanced-interval linear programming", 《EUROPEAN JOURNAL OF OPERATIONAL RESEARCH》 * |
PENG WEI等: "The Integration of Wind-Solar-Hydropower Generation in Enabling Economic Robust Dispatch", 《MATHEMATICAL PROBLEMS IN ENGINEERING》 * |
XIAOHONGRAN等: "A framework for uncertainty quantification and economic dispatch model with wind–solar energy", 《INTERNATIONAL JOURNAL OF ELECTRICAL POWER & ENERGY SYSTEMS》 * |
ZHAOYANZHANG等: "Energy management of ultra-short-term optimal scheduling of integrated energy system considering the characteristics of heating network", 《ENERGY》 * |
张沛等: "基于区间线性规划和可逼近理想解排序法的园区型综合能源系统投资决策", 《现代电力》 * |
闻若彤: "面向工业园区的多能源系统优化运行策略", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114781922B (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106372742A (zh) | 考虑电转气多源储能型微网日前最优经济调度方法 | |
CN112086998B (zh) | 孤岛模式下直流微电网的能量优化调度方法 | |
CN115099007B (zh) | 基于综合成本-能耗曲线的综合能源系统优化运行方法 | |
CN113659627A (zh) | 一种含光伏发电和液态空气储能的微电网优化调度方法 | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法 | |
CN117436773B (zh) | 一种含可中断负荷的独立微网资源规划方法和系统 | |
CN117371669B (zh) | 一种考虑碳交易风险成本的园区综合能源系统运行方法 | |
Mohamed et al. | Modelling and environmental/economic power dispatch of microgrid using multiobjective genetic algorithm optimization | |
CN116826867B (zh) | 提升源荷多元调峰主动性的优化调度及费用补偿分摊方法 | |
CN112883630A (zh) | 用于风电消纳的多微网系统日前优化经济调度方法 | |
CN111126675A (zh) | 多能互补微网系统优化方法 | |
CN116805192A (zh) | 考虑最优弃能率的综合能源系统双层规划优化方法及应用 | |
CN114781922A (zh) | 一种电热综合能源系统优化调度方法 | |
Garmabdari et al. | Multi energy system modelling and operation optimisation for university research facility | |
CN116187648A (zh) | 基于热电解耦的虚拟电厂热电联合优化调度方法 | |
CN115659585A (zh) | 计及需求响应的微能网低碳协同调度方法、装置、存储器及设备 | |
US20230119377A1 (en) | Method for the Control of an Energy System, and Associated Device | |
CN115085227A (zh) | 一种微电网源储容量配置方法及装置 | |
Rahmanzadeh et al. | Optimal energy management of microgrid based on fcchp in the presence of electric and thermal loads considering energy storage systems | |
Ma et al. | Multi-Point Layout Planning of Multi-Energy Power Supplies Based on Time-series Production Simulation | |
Cheng et al. | Optimal Economic Dispatch Strategy of a Hybrid Energy Storage Microgrid Based on Genetic Algorithm | |
Zhuo | Control of wind power smoothing with battery energy storage system and thermostatically controlled loads | |
Zhang et al. | Energy Storage Sizing for IES Considering Carbon Trading and Demand Response | |
Li et al. | An Optimal Operation Method of Hybrid AC/DC Distribution Network Under the Goal of Peaking Carbon Dioxide Emissions and Carbon Neutrality | |
Gao et al. | Optimized control strategy of green buildings with PV generation and energy storage integration |
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 |