CN107832979A - 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 - Google Patents
一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 Download PDFInfo
- Publication number
- CN107832979A CN107832979A CN201711280808.XA CN201711280808A CN107832979A CN 107832979 A CN107832979 A CN 107832979A CN 201711280808 A CN201711280808 A CN 201711280808A CN 107832979 A CN107832979 A CN 107832979A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- munder
- mtd
- 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 23
- 230000010354 integration Effects 0.000 title claims abstract description 20
- 238000005457 optimization Methods 0.000 title claims abstract description 20
- 238000004146 energy storage Methods 0.000 claims abstract description 34
- 239000002918 waste heat Substances 0.000 claims abstract description 16
- 239000000446 fuel Substances 0.000 claims abstract description 7
- 230000005611 electricity Effects 0.000 claims description 34
- 239000007789 gas Substances 0.000 claims description 28
- 238000004378 air conditioning Methods 0.000 claims description 17
- 238000001816 cooling Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 13
- AMXOYNBUYSYVKV-UHFFFAOYSA-M lithium bromide Chemical compound [Li+].[Br-] AMXOYNBUYSYVKV-UHFFFAOYSA-M 0.000 claims description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 238000010438 heat treatment Methods 0.000 claims description 11
- 238000010521 absorption reaction Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000010795 Steam Flooding Methods 0.000 claims description 6
- 239000003345 natural gas Substances 0.000 claims description 6
- 238000011084 recovery Methods 0.000 claims description 6
- 238000002485 combustion reaction Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 claims description 5
- 239000000567 combustion gas Substances 0.000 claims description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 3
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 239000003546 flue gas Substances 0.000 claims description 3
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 241000196324 Embryophyta Species 0.000 claims 1
- 238000005265 energy consumption Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005057 refrigeration Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000005619 thermoelectricity Effects 0.000 description 1
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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Coloring Foods And Improving Nutritive Qualities (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法。该方法首先搭建考虑能量梯级利用的工厂综合能源系统供能架构;以运行维护成本、购电成本、购热成本以及燃料成本构成的日运行费用最小为优化目标,基于冷热电功率平衡约束、设备物理约束和储能设备约束,建立考虑能量梯度利用的优化模型;对优化模型进行求解,得到工厂最优供能策略,实现工厂的自趋优调度。本发明考虑了冷热电多能耦合,实现了多种能源协同互补,引导用户制定合理的用能方案,提高了用户侧的用能效率,减少用户的用能成本;充分利用工厂中的循环余热,优化能量的梯级利用,可显著降低系统的日运行费用,较传统优化策略更贴近工程实际。
Description
技术领域
本发明涉及一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,属于综合能源及电力需求响应领域。
背景技术
能源是社会和经济发展的动力和基础。以智能电网为核心的能源系统革命在全世界范围内推动了电力系统的信息化和智能化。综合能源系统(integrated energy system,IES)是能源互联网的重要物理载体,是实现多能源互补、能量梯级利用等技术的关键。工业园区负荷需求量大、负荷类型多样,各种能量转换过程大多涉及热的利用。只有综合考虑能量的数量与品质两方面品质,才能对能量进行充分利用。因此有必要对工厂进行用能优化管理,实现能量梯度利用,提升工厂的经济效益和能源利用率。
发明内容
为了克服上述现有技术的不足,本发明提供了一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,该方法首先搭建考虑能源梯度利用的工厂综合能源系统的供能架构;以用户的日运行费用最小为目标,在冷热电平衡约束和设备多种设备约束的条件下,构建考虑工厂综合能源系统经济优化调度模型,通过调度各个设备的运行状态与出力,降低工厂的日运行费用;该方法可以应用在不同类型的工业园区综合能源系统中。
本发明的目的是通过以下技术方案来实现的:一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,该方法包括以下步骤:
步骤1:将热细分为高、中、低三种品位的热能,考虑其对口的温度利用区间及相应的利用技术,搭建考虑能量梯级利用的工厂综合能源系统供能架构,该架构具体为:
工厂通过集中式电力母线与公共电力交换电力,采用“自发自用、余量上网”的运行机制,优先满足本地各类负荷需求,同时允许将富裕电量输送到配电系统;同时,工厂内部无燃气生产,与燃气公司之间只存在单向购买行为;工业用户可向园区内的大型热电联产(combined heat and power,CHP)购买工业用汽,以满足工厂的蒸汽负荷需求;微型燃气轮机通过天然气燃烧释放出的高品位热进行发电,排出的高温烟气通过余热锅炉加以利用,余热锅炉通常被布置成两级,第一级产生中品位热蒸汽供吸收式溴化锂制冷机进行制冷或制热,第二级产生低品位热供应热水;燃气锅炉通过天然气燃烧产生的高品位热蒸汽以及从园区CHP中购买的高品位蒸汽用于满足工厂内的蒸汽驱动设备需求;蒸汽驱动设备使用高品位热蒸汽后,对蒸汽进行回收利用,产生中品位热和低品位热蒸汽,同样用于吸收式溴化锂制冷机制冷/制热以及供应热水;在供冷/热系统中,户用空调可提供空间冷热负荷,作为空间冷热负荷的调峰设备使用;冰蓄冷装置在制冷量满足当前所需且电价较低时储存冷水,待需要时释放冷量;
步骤2:以运行维护成本、购电成本、购热成本以及燃料成本构成的日运行费用最小为优化目标,基于冷热电功率平衡约束、设备物理约束和储能设备约束,建立考虑能量梯度利用的优化模型;其中,热功率平衡约束如下:
a)高品位热平衡约束:
b)中品位热平衡约束:
c)低品位热平衡约束:
式中:和分别为系统的高、中、低品位热功率;表示购热功率;表示第i个燃气锅炉产生的高品位热功率,和表示第i个蒸汽驱动设备的蒸汽需求热功率、中品位热回收功率和低品位热回收热功率;和分别表示第i个余热锅炉的中品位热回收热功率和低品位热回收热功率;表示第i个吸收式制冷机的耗热功率;表示热水负荷;
步骤3:对优化模型进行求解,得到工厂最优供能策略。
进一步地,所述步骤2中的优化目标为全天的运行费用:
Min CATC=COM+CES+CF+CH
1)运行维护成本:
式中:ξOM.i为设备i单位输出功率的运行维护费用;表示第i个设备在时段t的输出功率;T为单位时段长度;
2)购电成本:
式中:和分别为时段t的购电价格和购电功率;和分别为时段t的售电价格和售电功率;
3)燃料成本:
式中:和分别为时段t第i个燃气轮机和第i个燃气锅炉的燃气消耗速率;为气价;
4)购热成本:
式中,和分别为购热价格和购热功率;
约束条件:
1)电功率平衡约束
a)交流母线总负荷约束:
式中:为时段t的交流负荷;为第i个燃气轮机的产电功率;为交流电负荷;为交直流转化器的电功率;为第i个户用空调的耗电功率;为冰蓄冷空调系统的耗电功率;
b)交直流转换器效率约束:
式中:为时段t的直流母线总负荷;ηA/D为交流到直流的转换效率;ηD/A为直流到交流的转换效率;
c)直流母线总负荷约束:
式中:为光伏机组输出功率;和分别为电池储能的充电功率与放电功率;为时段t的直流负荷;
d)联络线约束与购售电状态约束:
式中:和分别为系统向电网购电和售电的功率上限;和分别为时段t处于购电和售电的0-1状态变量,取1表示购电,取1表示售电;
2)冷功率平衡约束
式中:和分别为吸收式制冷机、户用空调和冰蓄冷空调的供冷功率,为系统冷负荷。
3)设备运行约束
式中:和分别表示设备i在时段t的输入输出功率;和分别表示设备i在时段t的输出功率上下限;和分别表示设备i在时段t的输入功率上下限;
4)储能设备约束
储能设备需要满足储能状态约束与充放能功率约束;为了保证调度的连续性,调度周期前后,储能设备的储能量应保持一致;
SL.i=ST.i
式中:表示储能设备在时段t的储存量;和分别表示储能设备的最大和最小的储存容量;SL.i和ST.i分别表示储能设备的始状态和终状态;和分别表示储能设备的最大充电和放电功率;和分别表示储能设备在时段t处于充能和放能的0-1状态变量,取1表示充能,取1表示放能,保证设备不能同时充放能。
进一步地,所述步骤3中采用混合整数线性规划法对步骤2建立的模型进行优化求解,输出优化结果为各个设备的运行状态与出力。
与现有技术相比,本发明的有益效果是:一方面,考虑了冷热电多能耦合,实现了多种能源协同互补,引导用户制定合理的用能方案,提高了用户侧的用能效率,减少用户的用能成本;另一方面,充分利用了工厂中的循环余热,优化能量的梯级利用,可显著降低系统的日运行费用,较传统优化策略更贴近工程实际。本发明可适用于不同类型的工业园区综合能源系统中。
附图说明
图1为考虑能量梯级利用的工厂综合能源系统供能架构示意图;
图2为热能温度利用区间示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
本发明提供的一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,该方法包括以下步骤:
步骤1:将热细分为高、中、低三种品位的热能,考虑其对口的温度利用区间及相应的利用技术,搭建考虑能量梯级利用的工厂综合能源系统供能架构,具体为:
在工厂中,热能品位越高,其利用方式越多。经过若干热利用后,余热品位降低,可利用数目下降,需要高品位热驱动的循环子系统无法利用中、低品位余热。传统的工厂综合能源系统架构未充分考虑热能品位区别与能量的梯度利用。考虑能量梯级利用的工厂综合能源系统供能架构如图1所示,包含冷、热、电和气4种能源形式,系统中负荷种类多样、功能设备丰富,其主要设备有微型燃气轮机、光伏电池、余热锅炉、吸收式制冷机、户用空调、燃气锅炉、电池储能、蓄冷装置、蒸汽负荷。该系统通过集中式电力母线和公共电网交换电力,采用“自发自用、余量上网”的运行机制,优先满足本地各类负荷需求,同时允许将富裕电量输送到配电系统。同时,综合能源系统内部无燃气生产,与燃气公司之间只存在单向购买行为。园区内的工业用户可向CHP购买蒸汽,以满足蒸汽负荷需求。微型燃气轮机通过天然气燃烧释放出的高品位热进行发电,排出的高温烟气通过余热锅炉加以利用,余热锅炉通常被布置成两级,第一级产生中品位热蒸汽供吸收式溴化锂制冷机进行制冷或制热,第二级产生低品位热供应热水。燃气锅炉通过天然气燃烧产生的高品位热蒸汽以及从园区CHP中购买的高品位蒸汽用于满足工厂内的蒸汽驱动设备需求。蒸汽驱动设备使用高品位热蒸汽后,对蒸汽进行回收利用,产生中品位热和低品位热蒸汽,同样用于吸收式溴化锂制冷机制冷/制热以及供应热水。在供冷/热系统中,户用空调可提供空间冷热负荷,作为空间冷热负荷的调峰设备使用。冰蓄冷装置在制冷量满足当前所需且电价较低时储存冷水,待需要时释放冷量。
步骤2:以运行维护成本、购电成本、购热成本以及燃料成本构成的日运行费用最小为优化目标,基于冷热电功率平衡约束、设备物理约束和储能设备约束,建立考虑能量梯度利用的优化模型;优化目标为全天的运行费用:
Min CATC=COM+CES+CF+CH
1)运行维护成本:
式中:ξOM.i为设备i单位输出功率的运行维护费用;表示第i个设备在时段t的输出功率;T为单位时段长度;
2)购电成本:
式中:和分别为时段t的购电价格和购电功率;和分别为时段t的售电价格和售电功率;
3)燃料成本:
式中:和分别为时段t第i个燃气轮机和第i个燃气锅炉的燃气消耗速率;为气价;
4)购热成本:
式中,和分别为购热价格和购热功率;
约束条件:
1)电功率平衡约束
a)交流母线总负荷约束:
式中:为时段t的交流负荷;为第i个燃气轮机的产电功率;为交流电负荷;为交直流转化器的电功率;为第i个户用空调的耗电功率;为冰蓄冷空调系统的耗电功率;
b)交直流转换器效率约束:
式中:为时段t的直流母线总负荷;ηA/D为交流到直流的转换效率;ηD/A为直流到交流的转换效率;
c)直流母线总负荷约束:
式中:为光伏机组输出功率;和分别为电池储能的充电功率与放电功率;为时段t的直流负荷;
d)联络线约束与购售电状态约束:
式中:和分别为系统向电网购电和售电的功率上限;和分别为时段t处于购电和售电的0-1状态变量,取1表示购电,取1表示售电,系统不能同时购售电;
2)热功率平衡约束
a)高品位热平衡约束:
b)中品位热平衡约束:
c)低品位热平衡约束:
式中:和分别为系统的高、中、低品位热功率;表示购热功率;表示第i个燃气锅炉产生的高品位热功率,和表示第i个蒸汽驱动设备的蒸汽需求热功率、中品位热回收功率和低品位热回收热功率;和分别表示第i个余热锅炉的中品位热回收热功率和低品位热回收热功率;表示第i个吸收式制冷机的耗热功率;表示热水负荷;
3)冷功率平衡约束
式中:和分别为吸收式制冷机、户用空调和冰蓄冷空调的供冷功率,为系统冷负荷。
4)设备运行约束
式中:和分别表示设备i在时段t的输入输出功率;和分别表示设备i在时段t的输出功率上下限;和分别表示设备i在时段t的输入功率上下限;
5)储能设备约束
储能设备需要满足储能状态约束与充放能功率约束;为了保证调度的连续性,调度周期前后,储能设备的储能量应保持一致;
SL.i=ST.i
式中:表示储能设备在时段t的储存量;和分别表示储能设备的最大和最小的储存容量;SL.i和ST.i分别表示储能设备的始状态和终状态;和分别表示储能设备的最大充电和放电功率;和分别表示储能设备在时段t处于充能和放能的0-1状态变量,取1表示充能,取1表示放能,保证设备不能同时充放能。
步骤3:采用混合整数线性规划法对步骤2建立的模型进行优化求解,输出优化结果为各个设备的运行状态与出力,通过调节系统中各设备的运行方式与工作状态,降低工业用户的运行费用。
考虑能量梯级利用的工厂综合能源系统经济优化调度方法的优势在于:考虑了工厂内部存在的冷热电耦合,实现多能互补与能量梯级利用,提升工厂的能源利用率,减少了运行费用,引导用户选择最优的用能模式。同时,较传统优化策略更贴近工程实际。
上述对实施例的描述是为便于本技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对上述实施例做出各种修改,并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此,本发明不限于上述实施例,本领域技术人员根据本发明的揭示,对于本发明做出的改进和修改都应该在本发明的保护范围之内。
Claims (3)
1.一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,其特征在于,该方法包括以下步骤:
步骤1:将热细分为高、中、低三种品位的热能,考虑其对口的温度利用区间及相应的利用技术,搭建考虑能量梯级利用的工厂综合能源系统供能架构,该架构具体为:
工厂通过集中式电力母线与公共电力交换电力,优先满足本地各类负荷需求,同时允许将富裕电量输送到配电系统;同时,工厂内部无燃气生产,与燃气公司之间只存在单向购买行为;工业用户可向园区内的大型热电联产CHP购买工业用汽,以满足工厂的蒸汽负荷需求;微型燃气轮机通过天然气燃烧释放出的高品位热进行发电,排出的高温烟气通过余热锅炉加以利用,余热锅炉通常被布置成两级,第一级产生中品位热蒸汽供吸收式溴化锂制冷机进行制冷或制热,第二级产生低品位热供应热水;燃气锅炉通过天然气燃烧产生的高品位热蒸汽以及从园区CHP中购买的高品位蒸汽用于满足工厂内的蒸汽驱动设备需求;蒸汽驱动设备使用高品位热蒸汽后,对蒸汽进行回收利用,产生中品位热和低品位热蒸汽,同样用于吸收式溴化锂制冷机制冷/制热以及供应热水;在供冷/热系统中,户用空调可提供空间冷热负荷,作为空间冷热负荷的调峰设备使用;冰蓄冷装置在制冷量满足当前所需且电价较低时储存冷水,待需要时释放冷量;
步骤2:以运行维护成本、购电成本、购热成本以及燃料成本构成的日运行费用最小为优化目标,基于冷热电功率平衡约束、设备物理约束和储能设备约束,建立考虑能量梯度利用的优化模型;其中,热功率平衡约束如下:
a)高品位热平衡约束:
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>h</mi>
<mi>i</mi>
<mi>g</mi>
<mi>h</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>G</mi>
<mi>B</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>h</mi>
<mi>i</mi>
<mi>g</mi>
<mi>h</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>e</mi>
<mi>a</mi>
<mi>m</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
b)中品位热平衡约束:
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
<mi>d</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>e</mi>
<mi>a</mi>
<mi>m</mi>
<mo>.</mo>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
<mi>d</mi>
<mi>l</mi>
<mi>e</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
<mo>.</mo>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
<mi>d</mi>
<mi>l</mi>
<mi>e</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>G</mi>
<mi>B</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>e</mi>
<mi>a</mi>
<mi>m</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
<mi>d</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
c)低品位热平衡约束:
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>e</mi>
<mi>a</mi>
<mi>m</mi>
<mo>.</mo>
<mi>l</mi>
<mi>o</mi>
<mi>w</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
<mo>.</mo>
<mi>l</mi>
<mi>o</mi>
<mi>w</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
<mi>d</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<mo>&lsqb;</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>G</mi>
<mi>B</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>e</mi>
<mi>a</mi>
<mi>m</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&rsqb;</mo>
<mo>+</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>w</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
式中: 和分别为系统的高、中、低品位热功率;表示购热功率;表示第i个燃气锅炉产生的高品位热功率, 和表示第i个蒸汽驱动设备的蒸汽需求热功率、中品位热回收功率和低品位热回收热功率;和分别表示第i个余热锅炉的中品位热回收热功率和低品位热回收热功率;表示第i个吸收式制冷机的耗热功率;表示热水负荷;
步骤3:对优化模型进行求解,得到工厂最优供能策略。
2.根据权利要求1所述的一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,其特征在于,所述步骤2中的优化目标为全天的运行费用:
Min CATC=COM+CES+CF+CH
1)运行维护成本:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>O</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<msub>
<mi>&xi;</mi>
<mrow>
<mi>O</mi>
<mi>M</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
</msub>
<msubsup>
<mi>P</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mi>T</mi>
</mrow>
式中:ξOM.i为设备i单位输出功率的运行维护费用;表示第i个设备在时段t的输出功率;T为单位时段长度;
2)购电成本:
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&xi;</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&xi;</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</mrow>
式中:和分别为时段t的购电价格和购电功率;和分别为时段t的售电价格和售电功率;
3)燃料成本:
<mrow>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<mrow>
<mo>(</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>F</mi>
<mrow>
<mi>G</mi>
<mi>T</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>&xi;</mi>
<mrow>
<mi>g</mi>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>G</mi>
<mi>B</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>&xi;</mi>
<mrow>
<mi>g</mi>
<mi>a</mi>
<mi>s</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</mrow>
式中:和分别为时段t第i个燃气轮机和第i个燃气锅炉的燃气消耗速率;为气价;
4)购热成本:
<mrow>
<msub>
<mi>C</mi>
<mi>H</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<msubsup>
<mi>&xi;</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>h</mi>
</msubsup>
<msubsup>
<mi>H</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mi>T</mi>
</mrow>
式中,和分别为购热价格和购热功率;
约束条件:
1)电功率平衡约束
a)交流母线总负荷约束:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>T</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mo>-</mo>
<mi>l</mi>
<mi>o</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mo>-</mo>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>d</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
<mi>e</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
式中:为时段t的交流负荷;为第i个燃气轮机的产电功率;为交流电负荷;为交直流转化器的电功率;为第i个户用空调的耗电功率;为冰蓄冷空调系统的耗电功率;
b)交直流转换器效率约束:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
<mo>-</mo>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>A</mi>
<mo>/</mo>
<mi>D</mi>
</mrow>
</msub>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>D</mi>
<mo>/</mo>
<mi>A</mi>
</mrow>
</msub>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中:为时段t的直流母线总负荷;ηA/D为交流到直流的转换效率;ηD/A为直流到交流的转换效率;
c)直流母线总负荷约束:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
<mo>-</mo>
<mi>l</mi>
<mi>o</mi>
<mi>a</mi>
<mi>c</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mo>.</mo>
<mi>s</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>E</mi>
<mi>S</mi>
<mo>.</mo>
<mi>r</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
式中:为光伏机组输出功率;和分别为电池储能的充电功率与放电功率;为时段t的直流负荷;
d)联络线约束与购售电状态约束:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>y</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>s</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:和分别为系统向电网购电和售电的功率上限;和分别为时段t处于购电和售电的0-1状态变量,取1表示购电,取1表示售电;
2)冷功率平衡约束
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>d</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
<mi>e</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
式中:和分别为吸收式制冷机、户用空调和冰蓄冷空调的供冷功率,为系统冷负荷。
3)设备运行约束
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:和分别表示设备i在时段t的输入输出功率;和分别表示设备i在时段t的输出功率上下限;和分别表示设备i在时段t的输入功率上下限;
4)储能设备约束
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>L</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>r</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>r</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>w</mi>
<mrow>
<mi>r</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>w</mi>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>w</mi>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>r</mi>
<mo>.</mo>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:表示储能设备在时段t的储存量;和分别表示储能设备的最大和最小的储存容量;SL.i和ST.i分别表示储能设备的始状态和终状态;和分别表示储能设备的最大充电和放电功率;和分别表示储能设备在时段t处于充能和放能的0-1状态变量,取1表示充能,取1表示放能,保证设备不能同时充放能。
3.根据权利要求1所述的一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法,其特征在于:所述步骤3中采用混合整数线性规划法对步骤2建立的模型进行优化求解,输出优化结果为各个设备的运行状态与出力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711280808.XA CN107832979B (zh) | 2017-12-06 | 2017-12-06 | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711280808.XA CN107832979B (zh) | 2017-12-06 | 2017-12-06 | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107832979A true CN107832979A (zh) | 2018-03-23 |
CN107832979B CN107832979B (zh) | 2021-05-11 |
Family
ID=61641643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711280808.XA Active CN107832979B (zh) | 2017-12-06 | 2017-12-06 | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107832979B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108613247A (zh) * | 2018-04-02 | 2018-10-02 | 上海航天智慧能源技术有限公司 | 一种汽水两用燃气锅炉组的热负荷分配方法 |
CN108830743A (zh) * | 2018-05-25 | 2018-11-16 | 天津大学 | 考虑多种蓄冷装置的园区综合能源系统优化调度方法 |
CN109523137A (zh) * | 2018-10-29 | 2019-03-26 | 天津大学 | 考虑楼宇热负荷需求响应的园区综合能源优化调度方法 |
CN109685683A (zh) * | 2018-12-11 | 2019-04-26 | 新奥数能科技有限公司 | 用于泛能站的能源调度模型优化方法、装置、介质及设备 |
CN109740827A (zh) * | 2019-02-14 | 2019-05-10 | 华北电力大学 | 一种基于双层优化的区域综合能源系统规划优化方法 |
CN109764338A (zh) * | 2019-01-16 | 2019-05-17 | 广州环投设计研究院有限公司 | 一种垃圾焚烧发电厂余热梯级利用方法 |
CN110138004A (zh) * | 2019-05-08 | 2019-08-16 | 上海电力学院 | 一种多能互补系统优化运行方法 |
CN110458353A (zh) * | 2019-08-08 | 2019-11-15 | 上海交通大学 | 电热联合微网能量梯级优化方法及系统 |
CN110598313A (zh) * | 2019-09-10 | 2019-12-20 | 国网河北省电力有限公司 | 计及储能全寿命周期运维的综合能源系统优化配置方法 |
CN110619405A (zh) * | 2018-06-19 | 2019-12-27 | 新智数字科技有限公司 | 泛能站混合整数分段非线性优化方法及装置 |
CN110783917A (zh) * | 2019-11-05 | 2020-02-11 | 国网江苏省电力有限公司镇江供电分公司 | 含有新能源消纳的多能源能量枢纽配置方法 |
CN110889600A (zh) * | 2019-11-13 | 2020-03-17 | 国家电网有限公司 | 一种考虑柔性热负荷的区域综合能源系统优化调度方法 |
CN111121340A (zh) * | 2020-01-15 | 2020-05-08 | 安徽普泛能源技术有限公司 | 一种循环利用多品质余热的方法及系统 |
CN111241459A (zh) * | 2020-01-07 | 2020-06-05 | 国网浙江省电力有限公司杭州供电公司 | 一种考虑自趋优的工业园区分布式综合需求响应方法 |
CN112202203A (zh) * | 2020-09-30 | 2021-01-08 | 国网山东省电力公司经济技术研究院 | 一种电热微网能量耦合的梯级利用优化方法 |
CN112861335A (zh) * | 2021-02-01 | 2021-05-28 | 昆明理工大学 | 一种含p2g及储能的综合能源系统低碳经济调度方法 |
CN112883560A (zh) * | 2021-01-28 | 2021-06-01 | 华南理工大学 | 基于用户侧负荷削减响应的多能耦合供能网络的优化方法 |
CN113378409A (zh) * | 2021-07-06 | 2021-09-10 | 国网江苏省电力有限公司营销服务中心 | 综合能源系统多能互补优化调度方法及系统 |
CN113887894A (zh) * | 2021-09-17 | 2022-01-04 | 南京师范大学 | 一种综合能源系统中的高温蒸汽热泵提升能效及消纳风电的方法 |
CN114463130A (zh) * | 2022-02-08 | 2022-05-10 | 河北农业大学 | 基于梯型碳交易机制和需求响应的能源系统调度方法 |
CN116484999A (zh) * | 2023-01-16 | 2023-07-25 | 华北电力大学 | 考虑梯级利用的城镇综合能源系统调度优化方法及其应用 |
CN117151701A (zh) * | 2023-10-31 | 2023-12-01 | 山东欣历能源有限公司 | 一种热电联产的工业余热回收利用系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007228676A (ja) * | 2006-02-22 | 2007-09-06 | Hitachi Ltd | 発電所運用計画作成装置およびその方法 |
CN104734168A (zh) * | 2015-03-13 | 2015-06-24 | 山东大学 | 一种基于电热联合调度的微电网运行优化系统及方法 |
CN105869075A (zh) * | 2016-04-19 | 2016-08-17 | 东南大学 | 一种冷热电联供型微型能源网经济优化调度方法 |
CN107025519A (zh) * | 2017-03-30 | 2017-08-08 | 厦门大学 | 区域分布式冷热电多联产系统混合整数非线性模型优化方法 |
CN107194516A (zh) * | 2017-06-07 | 2017-09-22 | 华北电力大学 | 含多主体的多能互补微电网分布式优化调度方法 |
US20170308968A1 (en) * | 2014-10-21 | 2017-10-26 | Qin Zhou | System, method and apparatus for capacity determination for micro grid and tangible computer readable medium |
-
2017
- 2017-12-06 CN CN201711280808.XA patent/CN107832979B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007228676A (ja) * | 2006-02-22 | 2007-09-06 | Hitachi Ltd | 発電所運用計画作成装置およびその方法 |
US20170308968A1 (en) * | 2014-10-21 | 2017-10-26 | Qin Zhou | System, method and apparatus for capacity determination for micro grid and tangible computer readable medium |
CN104734168A (zh) * | 2015-03-13 | 2015-06-24 | 山东大学 | 一种基于电热联合调度的微电网运行优化系统及方法 |
CN105869075A (zh) * | 2016-04-19 | 2016-08-17 | 东南大学 | 一种冷热电联供型微型能源网经济优化调度方法 |
CN107025519A (zh) * | 2017-03-30 | 2017-08-08 | 厦门大学 | 区域分布式冷热电多联产系统混合整数非线性模型优化方法 |
CN107194516A (zh) * | 2017-06-07 | 2017-09-22 | 华北电力大学 | 含多主体的多能互补微电网分布式优化调度方法 |
Non-Patent Citations (1)
Title |
---|
何仲潇: ""考虑多能协同的工厂综合需求侧响应模型"", 《电力自动化设备》 * |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108613247A (zh) * | 2018-04-02 | 2018-10-02 | 上海航天智慧能源技术有限公司 | 一种汽水两用燃气锅炉组的热负荷分配方法 |
CN108830743B (zh) * | 2018-05-25 | 2021-10-15 | 天津大学 | 考虑多种蓄冷装置的园区综合能源系统优化调度方法 |
CN108830743A (zh) * | 2018-05-25 | 2018-11-16 | 天津大学 | 考虑多种蓄冷装置的园区综合能源系统优化调度方法 |
CN110619405B (zh) * | 2018-06-19 | 2022-12-02 | 新智数字科技有限公司 | 泛能站混合整数分段非线性优化方法及装置 |
CN110619405A (zh) * | 2018-06-19 | 2019-12-27 | 新智数字科技有限公司 | 泛能站混合整数分段非线性优化方法及装置 |
CN109523137A (zh) * | 2018-10-29 | 2019-03-26 | 天津大学 | 考虑楼宇热负荷需求响应的园区综合能源优化调度方法 |
CN109523137B (zh) * | 2018-10-29 | 2022-11-04 | 天津大学 | 考虑楼宇热负荷需求响应的园区综合能源优化调度方法 |
CN109685683A (zh) * | 2018-12-11 | 2019-04-26 | 新奥数能科技有限公司 | 用于泛能站的能源调度模型优化方法、装置、介质及设备 |
CN109764338A (zh) * | 2019-01-16 | 2019-05-17 | 广州环投设计研究院有限公司 | 一种垃圾焚烧发电厂余热梯级利用方法 |
CN109740827A (zh) * | 2019-02-14 | 2019-05-10 | 华北电力大学 | 一种基于双层优化的区域综合能源系统规划优化方法 |
CN110138004A (zh) * | 2019-05-08 | 2019-08-16 | 上海电力学院 | 一种多能互补系统优化运行方法 |
CN110458353A (zh) * | 2019-08-08 | 2019-11-15 | 上海交通大学 | 电热联合微网能量梯级优化方法及系统 |
CN110458353B (zh) * | 2019-08-08 | 2023-04-18 | 上海交通大学 | 电热联合微网能量梯级优化方法及系统 |
CN110598313B (zh) * | 2019-09-10 | 2023-05-16 | 国网河北省电力有限公司 | 计及储能全寿命周期运维的综合能源系统优化配置方法 |
CN110598313A (zh) * | 2019-09-10 | 2019-12-20 | 国网河北省电力有限公司 | 计及储能全寿命周期运维的综合能源系统优化配置方法 |
CN110783917A (zh) * | 2019-11-05 | 2020-02-11 | 国网江苏省电力有限公司镇江供电分公司 | 含有新能源消纳的多能源能量枢纽配置方法 |
CN110889600A (zh) * | 2019-11-13 | 2020-03-17 | 国家电网有限公司 | 一种考虑柔性热负荷的区域综合能源系统优化调度方法 |
CN111241459A (zh) * | 2020-01-07 | 2020-06-05 | 国网浙江省电力有限公司杭州供电公司 | 一种考虑自趋优的工业园区分布式综合需求响应方法 |
CN111241459B (zh) * | 2020-01-07 | 2023-07-21 | 国网浙江省电力有限公司杭州供电公司 | 一种考虑自趋优的工业园区分布式综合需求响应方法 |
CN111121340A (zh) * | 2020-01-15 | 2020-05-08 | 安徽普泛能源技术有限公司 | 一种循环利用多品质余热的方法及系统 |
CN112202203A (zh) * | 2020-09-30 | 2021-01-08 | 国网山东省电力公司经济技术研究院 | 一种电热微网能量耦合的梯级利用优化方法 |
CN112883560B (zh) * | 2021-01-28 | 2022-12-16 | 华南理工大学 | 基于用户侧负荷削减响应的多能耦合供能网络的优化方法 |
CN112883560A (zh) * | 2021-01-28 | 2021-06-01 | 华南理工大学 | 基于用户侧负荷削减响应的多能耦合供能网络的优化方法 |
CN112861335B (zh) * | 2021-02-01 | 2022-11-15 | 昆明理工大学 | 一种含p2g及储能的综合能源系统低碳经济调度方法 |
CN112861335A (zh) * | 2021-02-01 | 2021-05-28 | 昆明理工大学 | 一种含p2g及储能的综合能源系统低碳经济调度方法 |
CN113378409A (zh) * | 2021-07-06 | 2021-09-10 | 国网江苏省电力有限公司营销服务中心 | 综合能源系统多能互补优化调度方法及系统 |
CN113887894A (zh) * | 2021-09-17 | 2022-01-04 | 南京师范大学 | 一种综合能源系统中的高温蒸汽热泵提升能效及消纳风电的方法 |
CN114463130A (zh) * | 2022-02-08 | 2022-05-10 | 河北农业大学 | 基于梯型碳交易机制和需求响应的能源系统调度方法 |
CN116484999A (zh) * | 2023-01-16 | 2023-07-25 | 华北电力大学 | 考虑梯级利用的城镇综合能源系统调度优化方法及其应用 |
CN117151701A (zh) * | 2023-10-31 | 2023-12-01 | 山东欣历能源有限公司 | 一种热电联产的工业余热回收利用系统 |
CN117151701B (zh) * | 2023-10-31 | 2024-02-09 | 山东欣历能源有限公司 | 一种热电联产的工业余热回收利用系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107832979B (zh) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107832979A (zh) | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 | |
CN109919478B (zh) | 一种考虑综合供能可靠性的综合能源微网规划方法 | |
CN107609684A (zh) | 一种基于微电网的综合能源系统经济优化调度方法 | |
CN108717594B (zh) | 一种冷热电联供型多微网系统经济优化调度方法 | |
CN111445090B (zh) | 一种离网型综合能源系统双层规划方法 | |
CN104734168B (zh) | 一种基于电热联合调度的微电网运行优化系统及方法 | |
CN106786793B (zh) | 一种基于鲁棒优化的冷热电联供型微网运行方法 | |
CN106058942B (zh) | 考虑风电不确定性的含电转气和cchp的能量枢纽优化模型 | |
CN103728881B (zh) | 一种多楼宇冷热电联供系统的优化运行方法 | |
CN108229025A (zh) | 一种冷热电联供型多微网主动配电系统经济优化调度方法 | |
CN109146182A (zh) | 计及多种储能的分布式冷热电联供系统的经济调度方法 | |
CN109190859A (zh) | 冷热电联供型多微网系统及其经济优化调度方法 | |
CN109523065A (zh) | 一种基于改进量子粒子群算法的微能源网优化调度方法 | |
CN104616208A (zh) | 一种基于模型预测控制的冷热电联供型微电网运行方法 | |
CN113159407B (zh) | 基于区域综合能源系统的多能储存模块容量优化配置方法 | |
CN106950840A (zh) | 面向电网削峰的综合能源系统分层分布式协调控制方法 | |
CN108008629A (zh) | 一种多种能源互补利用系统的优化运行方法 | |
CN113722895A (zh) | 一种基于多站融合的综合能源系统优化配置方法 | |
CN114066204A (zh) | 一种综合能源系统的一体化优化规划与运行方法及装置 | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法 | |
CN113806952B (zh) | 一种考虑源-荷-储的冷热电综合能源系统及其优化运行方法 | |
CN115081700A (zh) | 基于综合储能技术的数据中心多能协同优化方法及系统 | |
CN109617052B (zh) | 一种大规模电储热单元智能化分层控制方法 | |
CN107565605A (zh) | 一种基于微电网的工厂设备自动趋向优化的方法 | |
CN111008760A (zh) | 一种冷热电联供微电网系统的双层优化方法 |
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 | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 310058 Yuhang Tang Road, Xihu District, Hangzhou, Zhejiang 866 Patentee after: ZHEJIANG University Patentee after: Guangzhou Power Supply Bureau of Guangdong Power Grid Co.,Ltd. Address before: 310058 Yuhang Tang Road, Xihu District, Hangzhou, Zhejiang 866 Patentee before: ZHEJIANG University Patentee before: GUANGZHOU POWER SUPPLY Co.,Ltd. |