CN117436672B - 考虑等效循环寿命和温控负荷的综合能源运行方法及系统 - Google Patents
考虑等效循环寿命和温控负荷的综合能源运行方法及系统 Download PDFInfo
- Publication number
- CN117436672B CN117436672B CN202311754468.5A CN202311754468A CN117436672B CN 117436672 B CN117436672 B CN 117436672B CN 202311754468 A CN202311754468 A CN 202311754468A CN 117436672 B CN117436672 B CN 117436672B
- Authority
- CN
- China
- Prior art keywords
- power
- cost
- chp
- constraint
- load
- 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
- 238000000034 method Methods 0.000 title claims abstract description 131
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 55
- 230000004044 response Effects 0.000 claims abstract description 53
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 50
- 238000005457 optimization Methods 0.000 claims abstract description 26
- 238000013486 operation strategy Methods 0.000 claims abstract description 15
- AMXOYNBUYSYVKV-UHFFFAOYSA-M lithium bromide Chemical compound [Li+].[Br-] AMXOYNBUYSYVKV-UHFFFAOYSA-M 0.000 claims description 86
- 239000007789 gas Substances 0.000 claims description 78
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 78
- 238000010248 power generation Methods 0.000 claims description 73
- 238000004146 energy storage Methods 0.000 claims description 54
- 239000002918 waste heat Substances 0.000 claims description 52
- 239000003345 natural gas Substances 0.000 claims description 39
- 238000005057 refrigeration Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 17
- 238000007599 discharging Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 16
- 230000005611 electricity Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 12
- 230000020169 heat generation Effects 0.000 claims description 12
- 238000010276 construction Methods 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 238000012423 maintenance Methods 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000001816 cooling Methods 0.000 claims description 6
- 230000002222 downregulating effect Effects 0.000 claims description 6
- 238000011084 recovery Methods 0.000 claims description 6
- 238000003303 reheating Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 3
- 102100031077 Calcineurin B homologous protein 3 Human genes 0.000 claims 4
- 101000777270 Homo sapiens Calcineurin B homologous protein 3 Proteins 0.000 claims 4
- 101001004953 Homo sapiens Lysosomal acid lipase/cholesteryl ester hydrolase Proteins 0.000 claims 4
- 102100026001 Lysosomal acid lipase/cholesteryl ester hydrolase Human genes 0.000 claims 4
- 230000003321 amplification Effects 0.000 claims 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims 2
- 230000007246 mechanism Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- GELKBWJHTRAYNV-UHFFFAOYSA-K lithium iron phosphate Chemical compound [Li+].[Fe+2].[O-]P([O-])([O-])=O GELKBWJHTRAYNV-UHFFFAOYSA-K 0.000 description 1
- 230000035515 penetration 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/06313—Resource planning in a project environment
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Finance (AREA)
- Tourism & Hospitality (AREA)
- Game Theory and Decision Science (AREA)
- Accounting & Taxation (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Biodiversity & Conservation Biology (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
- Feedback Control In General (AREA)
Abstract
一种考虑等效循环寿命和温控负荷的综合能源运行方法,先基于等效循环寿命、温控负荷、综合需求响应构建综合能源系统优化运行模型,再基于构建的模型进行计算,输出优化后的综合能源系统运行策略,该综合能源系统优化运行模型的目标函数为系统运行成本最小,约束条件包括功率平衡约束、设备出力约束、基线负荷约束。本发明能够合理调用多种可调性资源,缓解电能供需平衡压力;考虑等效循环寿命,能够更精准的计算电储能系统EES、电动汽车EV的运行成本;考虑了阶梯碳交易机制,能够精准计算系统低碳运行成本;同时引入考虑需求响应的现货市场出清收益,提高系统运行经济性并激励系统及用户提供需求响应服务。
Description
技术领域
本发明属于综合能源系统调度自动化技术领域,具体涉及一种考虑等效循环寿命和温控负荷的综合能源运行方法及系统。
背景技术
随着电力行业向低碳、可持续发展迈进,新能源在电网中的渗透比例不断上升,电源侧的电量平衡能力日益不足。在极端气候等不利因素的影响下,容易出现源荷两侧电量不平衡的情况,造成的电量缺口严重影响系统的安全稳定运行。
综合能源系统(IES)作为大工业用户的用能管理主体,负责安排系统内供能设备的运行计划,并为系统内用户提供各类能量流。IES通过多能流的耦合互补、协调管理和系统内综合需求响应(IDR)资源的合理调用,能够降低系统整体的运行成本,并具备向电网提供DR服务的能力。电储能系统(EES)、电动汽车(EV)和温控负荷是IES内的优质灵活性资源。
一方面,现有的IES运行策略研究主要对EES和EV的运行成本进行基于交互功率或放电深度的定值计算,但得到的结果缺乏精确性,这将影响IES在现货市场下的运行策略;另一方面,现有研究主要集中于利用DR和灵活性资源降低IES的运行成本,但未考虑IES可以调用系统内电储能设备、温控负荷额外的功率和容量,参与现货市场,为电网提供DR服务;还有,员工私有的EV在白天大量被闲置,具有极大的调度潜力。因此,迫切需要建立一种能够合理调用EES、EV和温控负荷等优质灵活性资源,考虑等效循环寿命、综合能源需求响应的综合能源系统优化运行模型,以优化IES运行策略,为电网提供优质需求响应服务。
发明内容
本发明的目的是针对现有技术存在的上述问题,提供一种能够合理调用EES、EV、温控负荷等优质灵活性资源,并且考虑等效循环寿命、综合能源需求响应,以优化IES运行策略的考虑等效循环寿命和温控负荷的综合能源运行方法及系统。
为实现以上目的,本发明的技术方案如下:
第一方面,本发明提出一种考虑等效循环寿命和温控负荷的综合能源运行方法,所述综合能源运行方法包括以下步骤:
S1、基于等效循环寿命、温控负荷、综合需求响应构建综合能源系统优化运行模型;所述综合能源系统优化运行模型的目标函数为系统运行成本最小,约束条件包括功率平衡约束、设备出力约束、基线负荷约束;
S2、基于构建的综合能源系统优化运行模型进行计算,输出优化后的综合能源系统运行策略。
所述目标函数为:
;
;
;
;
;
;
;
上式中,为系统运行成本;/>、/>、/>、/>、/>分别为能源购置成本、设备运行成本、多类型需求响应成本、碳排放成本、现货市场出清收益;/>、/>分别为综合能源系统IES向电网购电的价格和功率;/>、/>分别为综合能源系统IES从气网购气的价格、气网输入综合能源系统IES的天然气功率;/>、/>分别为热电联产机组CHP的启、停成本;/>、/>、/>、/>、/>、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB、溴化锂制冷系统LBRS、卡琳娜循环发电系统KLN、热储能系统TES、冷储能系统CES的运行成本;/>、/>分别为需求响应前、需求响应后的用户购电成本;/>、/>、/>、/>分别为冷、热、电、气四种能源对应的功率;/>为/>时刻的综合需求响应DR资源出清价格;/>为/>时刻综合能源系统IES第m段投标的中标量;/>为规划期限;/>为投标段数;/>为温度调整成本;/>为电动汽车EV参与综合需求响应DR调用的补偿成本;/>为电储能系统EES的循环寿命损耗成本;/>为碳交易基价;/>为碳排放交易额;/>为碳交易价格增长率;/>为阶梯式碳交易区间长度;/>为电动汽车EV参与综合能源系统IES内综合需求响应DR调用的经济补偿系数;/>为电动汽车EV的循环寿命损耗成本。
所述电储能系统EES的循环寿命损耗成本、电动汽车EV的循环寿命损耗成本基于以下公式计算得到:
;
;
;
;
上式中,、/>、/>分别为电储能系统EES的单位功率投资成本、运维成本、单位容量投资成本;/>、/>分别为电储能系统EES的额定功率、额定容量;/>、/>分别为电动汽车EV电池单位容量成本、电动汽车EV额定容量;/>、/>分别为电储能系统EES、电动汽车EV的电池寿命损耗成本;/>、/>分别为电储能系统EES、电动汽车EV在100%放电深度下的等效循环次数;/>、/>分别为电储能系统EES、电动汽车EV的线性化等效循环次数;
所述温度调整成本基于温控负荷等效热参数模型得到,所述温控负荷等效热参数模型为:
;
;
;
;
上式中,、/>分别为室内温度下限、上限;/>、/>、/>分别为/>时刻的室内设定温度、室内下调温度、室内上调温度;/>为/>时刻的室内温度;/>为整日调整温度累计值;/>为温度调整惩罚系数;/>为惩罚增幅系数;/>为单位区间。
所述功率平衡约束包括电功率平衡约束、热功率平衡约束、冷功率平衡约束、天然气平衡约束;其中,所述电功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、光伏系统PV、卡琳娜循环发电系统KLN生产的电功率;/>、/>、/>分别为电锅炉CACS、电锅炉EB、工业制冷机IR消耗的电功率;/>、/>分别为电储能系统EES的充、放电功率;/>、/>为电动汽车EV的充、放电功率;/>为/>时刻的生活类负荷功率;/>为/>时刻的工业类负荷功率;
所述热功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB生产的热功率;/>、/>分别为热储能系统TES的充、放热功率;/>为/>时刻工业类负荷产生的冷功率;
所述冷功率平衡约束为:
;
上式中,、/>分别为工业制冷机IR、溴化锂制冷系统LBRS生产的冷功率;、/>分别为电锅炉CACS充、放冷功率;/>为/>时刻工业类负荷产生的冷功率;
所述天然气平衡约束为:
;
上式中,、/>分别为热电联产机组CHP、燃气锅炉GB耗费的天然气功率;为/>时刻工业类负荷产生的天然气功率;
所述设备出力约束包括热电联产机组约束、燃汽锅炉约束、余热锅炉约束、溴化锂制冷机约束、卡琳娜循环发电系统约束,所述热电联产机组约束为:
;
上式中,为热电联产机组CHP最大产热功率时对应的发电功率;/>、、/>分别为热电联产机组CHP的最小发电功率、最大发电功率和发电效率;/>、/>分别为热电联产机组CHP的热电耦合出力区间的下限斜率、上限斜率;/>为热电联产机组CHP的启停变量;/>为天然气热值;/>为热电联产机组CHP热电联产过程中产生的余热功率;/>为热电联产机组CHP最大产热功率;
所述燃汽锅炉约束为:
;
上式中,为燃气锅炉GB的产热效率;/>为燃气锅炉GB的最大产热功率;
所述余热锅炉约束为:
;
上式中,为余热锅炉利用高温热气残余热量再加热后输出的热功率;/>为余热锅炉热能回收效率;/>为/>时刻燃气锅炉GB产生的余热功率;
所述溴化锂制冷机约束为:
;
;
上式中,为溴化锂制冷系统LBRS从系统余热功率吸收的热功率;/>为余热利用系数;/>为溴化锂制冷系统LBRS的热冷转换效率;/>为/>时刻光伏系统产热功率;/>为/>时刻溴化锂制冷系统产热功率;
所述卡琳娜循环发电系统约束为:
;
上式中,为卡琳娜循环发电系统KLN从系统余热功率吸收的热功率;/>为卡琳娜循环发电系统KLN的热电转换效率;
所述基线负荷约束为:
;
;
上式中,为综合能源系统IES未提供综合需求响应DR服务的综合能源系统IES基线负荷,对应综合能源系统IES的日常生产计划。
第二方面,本发明提出考虑等效循环寿命和温控负荷的综合能源运行系统,所述运行系统包括模型构建模块、优化计算模块;
所述模型构建模块用于基于等效循环寿命、温控负荷、综合需求响应构建综合能源系统优化运行模型,所述综合能源系统优化运行模型的目标函数为系统运行成本最小,约束条件包括功率平衡约束、设备出力约束、基线负荷约束;
所述优化计算模块用于基于综合能源系统优化运行模型进行计算,输出优化后的综合能源系统运行策略。
所述模型构建模块用于如下目标函数:
所述目标函数为:
;
;
;
;
;
;
上式中,为系统运行成本;/>、/>、/>、/>、/>分别为能源购置成本、设备运行成本、需求响应成本、碳排放成本、现货市场出清收益;/>、/>分别为综合能源系统IES向电网购电的价格和功率;/>、/>分别为综合能源系统IES从气网购气的价格、气网输入综合能源系统IES的天然气功率;/>、/>分别为热电联产机组CHP的启、停成本;/>、/>、/>、/>、/>、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB、溴化锂制冷系统LBRS、卡琳娜循环发电系统KLN、热储能系统TES、冷储能系统CES的运行成本;/>、/>分别为需求响应前、需求响应后的用户购电成本;/>、/>、/>、/>分别为冷、热、电、气四种能源对应的功率;/>为/>时刻的综合需求响应DR资源出清价格;/>为/>时刻综合能源系统IES第m段投标的中标量;/>为规划期限;/>为投标段数;/>为温度调整成本;/>为电动汽车EV参与综合需求响应DR调用的补偿成本;/>为电储能系统EES的循环寿命损耗成本;/>为碳交易基价;/>为碳排放交易额;/>为碳交易价格增长率;/>为阶梯式碳交易区间长度;/>为电动汽车EV参与综合能源系统IES内综合需求响应DR调用的经济补偿系数;/>为电动汽车EV的循环寿命损耗成本。
所述电储能系统EES的循环寿命损耗成本、电动汽车EV的循环寿命损耗成本基于以下公式计算得到:
;
;
;
;
上式中,、/>、/>分别为电储能系统EES的单位功率投资成本、运维成本、单位容量投资成本;/>、/>分别为电储能系统EES的额定功率、额定容量;/>、/>分别为电动汽车EV电池单位容量成本、电动汽车EV额定容量;/>、/>分别为电储能系统EES、电动汽车EV的电池寿命损耗成本;/>、/>分别为电储能系统EES、电动汽车EV在100%放电深度下的等效循环次数;/>、/>分别为电储能系统EES、电动汽车EV的线性化等效循环次数;
所述温度调整成本基于温控负荷等效热参数模型得到,所述温控负荷等效热参数模型为:
;
;
;
;
上式中,、/>分别为室内温度下限、上限;/>、/>、/>分别为/>时刻的室内设定温度、室内下调温度、室内上调温度;/>为/>时刻的室内温度;/>为整日调整温度累计值;/>为温度调整惩罚系数;/>为惩罚增幅系数;/>为单位区间。
所述功率平衡约束包括电功率平衡约束、热功率平衡约束、冷功率平衡约束、天然气平衡约束;其中,所述电功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、光伏系统PV、卡琳娜循环发电系统KLN生产的电功率;/>、/>、/>分别为电锅炉CACS、电锅炉EB、工业制冷机IR消耗的电功率;/>、/>分别为电储能系统EES的充、放电功率;/>、/>为电动汽车EV的充、放电功率;/>为/>时刻的生活类负荷功率;/>为/>时刻的工业类负荷功率;
所述热功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB生产的热功率;/>、/>分别为热储能系统TES的充、放热功率;/>为/>时刻工业类负荷产生的冷功率;
所述冷功率平衡约束为:
;
上式中,、/>分别为工业制冷机IR、溴化锂制冷系统LBRS生产的冷功率;、/>分别为电锅炉CACS充、放冷功率;/>为/>时刻工业类负荷产生的冷功率;
所述天然气平衡约束为:
;
上式中,、/>分别为热电联产机组CHP、燃气锅炉GB耗费的天然气功率;为/>时刻工业类负荷产生的天然气功率;
所述设备出力约束包括热电联产机组约束、燃汽锅炉约束、余热锅炉约束、溴化锂制冷机约束、卡琳娜循环发电系统约束,所述热电联产机组约束为:
;
上式中,为热电联产机组CHP最大产热功率时对应的发电功率;/>、、/>分别为热电联产机组CHP的最小发电功率、最大发电功率和发电效率;/>、/>分别为热电联产机组CHP的热电耦合出力区间的下限斜率、上限斜率;/>为热电联产机组CHP的启停变量;/>为天然气热值;/>为热电联产机组CHP热电联产过程中产生的余热功率;/>为热电联产机组CHP最大产热功率;
所述燃汽锅炉约束为:
;
上式中,为燃气锅炉GB的产热效率;/>为燃气锅炉GB的最大产热功率;
所述余热锅炉约束为:
;
上式中,为余热锅炉利用高温热气残余热量再加热后输出的热功率;/>为余热锅炉热能回收效率;/>为/>时刻燃气锅炉GB产生的余热功率;
所述溴化锂制冷机约束为:
;
;
上式中,为溴化锂制冷系统LBRS从系统余热功率吸收的热功率;/>为余热利用系数;/>为溴化锂制冷系统LBRS的热冷转换效率;/>为/>时刻光伏系统产热功率;/>为/>时刻溴化锂制冷系统产热功率;
所述卡琳娜循环发电系统约束为:
;
上式中,为卡琳娜循环发电系统KLN从系统余热功率吸收的热功率;/>为卡琳娜循环发电系统KLN的热电转换效率;
所述基线负荷约束为:
;
;
上式中,为综合能源系统IES未提供综合需求响应DR服务的综合能源系统IES基线负荷,对应综合能源系统IES的日常生产计划。
第三方面,本发明考虑等效循环寿命和温控负荷的综合能源运行设备,所述运行设备包括存储器和处理器;
所述存储器,用于存储计算机程序代码,并将所述计算机程序代码传输给所述处理器;
所述处理器,用于根据所述计算机程序代码中的指令执行前述的方法。
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现前述的方法。
与现有技术相比,本发明的有益效果为:
本发明一种考虑等效循环寿命和温控负荷的综合能源运行方法,构建了现货市场下含电储能系统EES、电动汽车EV、温控负荷的综合能源系统优化运行模型,能够合理调用综合能源系统IES内多种可调性资源,为IES提供最优DR运行策略,缓解电能供需平衡压力,提高综合能源系统IES的整体经济性;该综合能源系统优化运行模型的目标函数为系统运行成本最小,在系统运行成本的计算中,考虑等效循环寿命,能够更精准的计算电储能系统EES、电动汽车EV的运行成本;考虑了阶梯碳交易机制,能够精准计算系统低碳运行成本;同时,引入考虑需求响应的现货市场出清收益,能够提高系统运行经济性并激励系统及用户提供需求响应服务。
附图说明
图1为本发明所述运行方法的流程图。
图2为本发明所述运行系统的结构框图。
图3为本发明所述运行设备的结构框图。
具体实施方式
下面结合具体实施方式以及附图对本发明作进一步详细的说明。
本发明所面对的综合能源系统IES内部存在电、热、冷、气4种能源形式和多类能源设备,能源设备的主要组成为:热电联产机组CHP、燃气锅炉GB、光伏系统PV、电锅炉EB、CACS、卡琳娜循环发电系统KLN、溴化锂制冷系统LBRS、工业制冷机IR、电储能系统EES、冷储能系统CES、热储能系统TES。综合能源系统IES通过设备间的能量转换和系统内的能流互补,实现低碳经济运行,并通过调整运行策略来削减购电功率,从而为电网提供DR服务。
实施例1:
参见图1,本发明一种考虑等效循环寿命和温控负荷的综合能源运行方法,依次按照以下步骤进行:
S1、构建综合能源系统优化运行模型;所述综合能源系统优化运行模型的目标函数为系统运行成本最小,所述目标函数具体为:
;
;
;
;
;
;
;
上式中,为系统运行成本;/>、/>、/>、/>、/>分别为能源购置成本、设备运行成本、多类型需求响应成本、碳排放成本、现货市场出清收益;/>、/>分别为综合能源系统IES向电网购电的价格和功率;/>、/>分别为综合能源系统IES从气网购气的价格、气网输入综合能源系统IES的天然气功率;/>、/>分别为热电联产机组CHP的启、停成本;/>、/>、/>、/>、/>、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB、溴化锂制冷系统LBRS、卡琳娜循环发电系统KLN、热储能系统TES、冷储能系统CES的运行成本;/>、/>分别为需求响应前、需求响应后的用户购电成本;/>、/>、/>、/>分别为冷、热、电、气四种能源对应的功率;/>为/>时刻的综合需求响应DR资源出清价格;/>为/>时刻综合能源系统IES第m段投标的中标量;/>为规划期限;/>为投标段数;/>为温度调整成本;/>为电动汽车EV参与综合需求响应DR调用的补偿成本;/>为电储能系统EES的循环寿命损耗成本;/>为碳交易基价;/>为碳排放交易额;/>为碳交易价格增长率;/>为阶梯式碳交易区间长度;/>为电动汽车EV参与综合能源系统IES内综合需求响应DR调用的经济补偿系数;/>为电动汽车EV的循环寿命损耗成本;
所述模型的约束条件包括功率平衡约束、设备出力约束、基线负荷约束;所述功率平衡约束包括电功率平衡约束、热功率平衡约束、冷功率平衡约束、天然气平衡约束;所述设备出力约束包括热电联产机组约束、燃汽锅炉约束、余热锅炉约束、溴化锂制冷机约束、卡琳娜循环发电系统约束,其中,
所述电功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、光伏系统PV、卡琳娜循环发电系统KLN生产的电功率;/>、/>、/>分别为电锅炉CACS、电锅炉EB、工业制冷机IR消耗的电功率;/>、/>分别为电储能系统EES的充、放电功率;/>、/>为电动汽车EV的充、放电功率;/>为/>时刻的生活类负荷功率;/>为/>时刻的工业类负荷功率;
所述热功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB生产的热功率;/>、/>分别为热储能系统TES的充、放热功率;/>为/>时刻工业类负荷产生的冷功率;
所述冷功率平衡约束为:
;
上式中,、/>分别为工业制冷机IR、溴化锂制冷系统LBRS生产的冷功率;、/>分别为电锅炉CACS充、放冷功率;/>为/>时刻工业类负荷产生的冷功率;
所述天然气平衡约束为:
;
上式中,、/>分别为热电联产机组CHP、燃气锅炉GB耗费的天然气功率;为/>时刻工业类负荷产生的天然气功率;
所述热电联产机组约束为:
;
上式中,为热电联产机组CHP最大产热功率时对应的发电功率;/>、、/>分别为热电联产机组CHP的最小发电功率、最大发电功率和发电效率;/>、/>分别为热电联产机组CHP的热电耦合出力区间的下限斜率、上限斜率;/>为热电联产机组CHP的启停变量;/>为天然气热值;/>为热电联产机组CHP热电联产过程中产生的余热功率;/>为热电联产机组CHP最大产热功率;
所述燃汽锅炉约束为:
;
上式中,为燃气锅炉GB的产热效率;/>为燃气锅炉GB的最大产热功率;
所述余热锅炉约束为:
;
上式中,为余热锅炉利用高温热气残余热量再加热后输出的热功率;/>为余热锅炉热能回收效率;/>为/>时刻燃气锅炉GB产生的余热功率;
所述溴化锂制冷机约束为:
;
;/>
上式中,为溴化锂制冷系统LBRS从系统余热功率吸收的热功率;/>为余热利用系数;/>为溴化锂制冷系统LBRS的热冷转换效率;/>为/>时刻光伏系统产热功率;/>为/>时刻溴化锂制冷系统产热功率;
所述卡琳娜循环发电系统约束为:
;
上式中,为卡琳娜循环发电系统KLN从系统余热功率吸收的热功率;/>为卡琳娜循环发电系统KLN的热电转换效率;
根据用户基线负荷认定IES提供DR服务的实施效果,为避免综合能源系统IES参与综合需求响应DR服务导致其余时段的负荷高升,造成非高峰时段电网电量不平衡,设置基线负荷约束;所述基线负荷约束为:
;
;
上式中,为综合能源系统IES未提供综合需求响应DR服务的综合能源系统IES基线负荷,对应综合能源系统IES的日常生产计划;
上述综合能源系统优化运行模型中,由于现阶段国内的配额方式主要为无偿分配,碳排放交易额的计算公式如下:
;
上式中,为实际碳排放量;/>、/>分别为热电联产机组CHP、燃气锅炉GB的初始碳排放分配量;/>、/>分别为单位电功率、单位热功率的碳排放配量;/>为单位天然气燃烧产生的实际碳排放量;
上述综合能源系统优化运行模型中,电储能系统EES的循环寿命损耗成本、电动汽车EV的循环寿命损耗成本/>的计算公式及推导过程如下:/>
步骤A1、引入布尔变量构建线性化计算模型,将等效循环次数线性化,得到电储能系统EES的线性化等效循环次数/>、电动汽车EV的线性化等效循环次数/>;
电储能系统EES的等效循环次数为:
;
上式中,为电储能系统EES在100%放电深度下的等效循环次数;/>为电储能系统EES在标准放电深度/>下的等效循环寿命;/>为放电深度;/>为雨流计数法拟合指数;
电动汽车EV的等效循环次数为:
;
上式中,为电动汽车EV在100%放电深度下的等效循环次数;/>为电动汽车EV在标准放电深度/>下的等效循环寿命;/>、/>、/>均为雨流计数法拟合指数;、/>分别为放电循环的起始荷电状态和标准荷电状态;
因等效循环次数为指数函数,以为决策变量,选择分段点/>,所述布尔变量/>为:
;
上式中,、/>均为分段线性化函数参数,线性化的拟合程度取决于分段点的选取,分段点越多,取点位置越合适,则拟合程度越高;
电储能系统EES的线性化计算模型为:
;
上式中,、/>分别为辅助端点变量、辅助分段变量;/>为分段数量;
电动汽车EV的线性化计算模型构建流程与电储能系统EES一致;
步骤A2、调节DOD,通过计算充放电状态切换次数和状态切换前的累计电量来精确测算电化学电池的运行周期下的充放电深度,根据以下公式计算电储能系统EES的循环寿命损耗成本、电动汽车EV的循环寿命损耗成本/>:
;
;
;
;
上式中,、/>、/>分别为电储能系统EES的单位功率投资成本、运维成本、单位容量投资成本;/>、/>分别为电储能系统EES的额定功率、额定容量;/>、/>分别为电动汽车EV电池单位容量成本、电动汽车EV额定容量;/>、/>分别为电储能系统EES、电动汽车EV的电池寿命损耗成本;/>、/>分别为电储能系统EES、电动汽车EV的线性化等效循环次数;
上述综合能源系统优化运行模型中,温度调整成本基于温控负荷等效热参数模型得到,所述温控负荷等效热参数模型的表达式如下:
;
;
;
;
上式中,、/>分别为室内温度下限、上限;/>、/>、/>分别为/>时刻的室内设定温度、室内下调温度、室内上调温度;/>为/>时刻的室内温度;/>为整日调整温度累计值;/>为温度调整惩罚系数;/>为惩罚增幅系数;/>为单位区间;
S2、基于构建的综合能源系统优化运行模型进行计算,输出优化后的综合能源系统运行策略。
下面为验证本发明所提考虑等效循环寿命和温控负荷的综合能源运行方法的有效性,在综合能源系统上分别利用以下4种方法进行优化,方法1为本发明所提方法,方法2为考虑交换功率成本和温控负荷的综合能源系统运行方法,方法3为将系统中的室外气温上升2℃后再运行本发明所提方法,方法4为考虑等效循环寿命和温控负荷,不考虑低碳的综合能源系统运行方法;
参数设置:天然气热值为3.6万;购买天然气价格为3.5元/>;电储能系统采用磷酸铁锂电池;单位功率投资成本为2200元/kW;单位功率维护成本200元/kW;单位容量投资成本为1500元/(kW•h);冷、热、电能源的单位需求响应补偿设为2元/kW;
(1)运行由4种方法输出的优化后综合能源系统运行策略,计算各项成本并进行对比分析,各项成本的计算结果如表1所示,DR中标量和出清价格为800kW•h、1.73元/(kW•h):
表1、4种方法的各项成本
由表1可知,对比方法1和方法2,方法2的总运行成本较方法1更低,这是因为方法2中对电储能系统EES、电动汽车EV调用成本的计算不够准确,会影响IES中其他能源设备的运行结果,方法1在成本计算中充分考虑电化学电池的循环寿命更为合理;
对比方法1和方法3,更高的室外温度将增加综合能源系统IES的各项运行成本,方法1考虑了温控负荷的影响,能够反映室内外温度、辐照强度、等效热源等对综合能源系统IES运行策略的影响,从而提高了模型结果的准确性;
对比方法1和方法4,方法4未考虑阶梯碳交易成本时,碳排放成本上升至1103.6元,可见方法1在考虑阶梯碳交易后,通过减少高排碳设备的使用成功降低了碳排放,这不仅有助于环境保护和解决电量不平衡问题,同时降低了整体运行成本,促进实现“双碳”目标。
(2)在不同DR中标量和出清价格下运行方法1,计算调用DR资源成本以及系统总成本,计算结果如表2所示:
表2、不同DR中标量和出清价格所对应的调用DR资源成本以及系统总成本
/>
由表2可知,随着综合能源运行系统IES调用内部DR资源成本上升,但系统总运行成本减小,证明本发明所提方法能够激励综合能源运行系统IES为电网提供DR服务。
实施例2:
参见图2,一种考虑等效循环寿命和温控负荷的综合能源运行系统,包括模型构建模块、优化计算模块;所述模型构建模块用于基于等效循环寿命、温控负荷、综合需求响应构建如下的综合能源系统优化运行模型,所述综合能源系统优化运行模型包括目标函数及约束条件:
目标函数为:
;
;
;
;
;
;
上式中,为系统运行成本;/>、/>、/>、/>、/>分别为能源购置成本、设备运行成本、需求响应成本、碳排放成本、现货市场出清收益;/>、/>分别为综合能源系统IES向电网购电的价格和功率;/>、/>分别为综合能源系统IES从气网购气的价格、气网输入综合能源系统IES的天然气功率;/>、/>分别为热电联产机组CHP的启、停成本;/>、/>、/>、/>、/>、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB、溴化锂制冷系统LBRS、卡琳娜循环发电系统KLN、热储能系统TES、冷储能系统CES的运行成本;/>、/>分别为需求响应前、需求响应后的用户购电成本;/>、/>、/>、/>分别为冷、热、电、气四种能源对应的功率;/>为/>时刻的综合需求响应DR资源出清价格;/>为/>时刻综合能源系统IES第m段投标的中标量;/>为规划期限;/>为投标段数;/>为温度调整成本;/>为电动汽车EV参与综合需求响应DR调用的补偿成本;/>为电储能系统EES的循环寿命损耗成本;/>为碳交易基价;/>为碳排放交易额;/>为碳交易价格增长率;/>为阶梯式碳交易区间长度;/>为电动汽车EV参与综合能源系统IES内综合需求响应DR调用的经济补偿系数;/>为电动汽车EV的循环寿命损耗成本;
基于以下公式计算电储能系统EES的循环寿命损耗成本、电动汽车EV的循环寿命损耗成本/>:
;
;
;
;
上式中,、/>、/>分别为电储能系统EES的单位功率投资成本、运维成本、单位容量投资成本;/>、/>分别为电储能系统EES的额定功率、额定容量;/>、/>分别为电动汽车EV电池单位容量成本、电动汽车EV额定容量;/>、/>分别为电储能系统EES、电动汽车EV的电池寿命损耗成本;/>、/>分别为电储能系统EES、电动汽车EV在100%放电深度下的等效循环次数;/>、/>分别为电储能系统EES、电动汽车EV的线性化等效循环次数;
基于温控负荷等效热参数模型计算温度调整成本,所述温控负荷等效热参数模型为:
;
;
;/>
;
上式中,、/>分别为室内温度下限、上限;/>、/>、/>分别为/>时刻的室内设定温度、室内下调温度、室内上调温度;/>为/>时刻的室内温度;/>为整日调整温度累计值;/>为温度调整惩罚系数;/>为惩罚增幅系数;/>为单位区间;
约束条件包括低碳经济运行约束条件、功率平衡约束、设备出力约束、基线负荷约束;所述功率平衡约束包括电功率平衡约束、热功率平衡约束、冷功率平衡约束、天然气平衡约束;所述设备出力约束包括热电联产机组约束、燃汽锅炉约束、余热锅炉约束、溴化锂制冷机约束、卡琳娜循环发电系统约束;
所述电功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、光伏系统PV、卡琳娜循环发电系统KLN生产的电功率;/>、/>、/>分别为电锅炉CACS、电锅炉EB、工业制冷机IR消耗的电功率;/>、/>分别为电储能系统EES的充、放电功率;/>、/>为电动汽车EV的充、放电功率;/>为/>时刻的生活类负荷功率;/>为/>时刻的工业类负荷功率;
所述热功率平衡约束为:
;
上式中,、/>、/>分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB生产的热功率;/>、/>分别为热储能系统TES的充、放热功率;/>为/>时刻工业类负荷产生的冷功率;
所述冷功率平衡约束为:
;
上式中,、/>分别为工业制冷机IR、溴化锂制冷系统LBRS生产的冷功率;、/>分别为电锅炉CACS充、放冷功率;/>为/>时刻工业类负荷产生的冷功率;
所述天然气平衡约束为:
;
上式中,、/>分别为热电联产机组CHP、燃气锅炉GB耗费的天然气功率;为/>时刻工业类负荷产生的天然气功率;
所述热电联产机组约束为:
;
上式中,为热电联产机组CHP最大产热功率时对应的发电功率;/>、、/>分别为热电联产机组CHP的最小发电功率、最大发电功率和发电效率;/>、/>分别为热电联产机组CHP的热电耦合出力区间的下限斜率、上限斜率;/>为热电联产机组CHP的启停变量;/>为天然气热值;/>为热电联产机组CHP热电联产过程中产生的余热功率;/>为热电联产机组CHP最大产热功率;
所述燃汽锅炉约束为:
;
上式中,为燃气锅炉GB的产热效率;/>为燃气锅炉GB的最大产热功率;
所述余热锅炉约束为:
;
上式中,为余热锅炉利用高温热气残余热量再加热后输出的热功率;/>为余热锅炉热能回收效率;/>为/>时刻燃气锅炉GB产生的余热功率;
所述溴化锂制冷机约束为:
;
;/>
上式中,为溴化锂制冷系统LBRS从系统余热功率吸收的热功率;/>为余热利用系数;/>为溴化锂制冷系统LBRS的热冷转换效率;/>为/>时刻光伏系统产热功率;/>为/>时刻溴化锂制冷系统产热功率;
所述卡琳娜循环发电系统约束为:
;
上式中,为卡琳娜循环发电系统KLN从系统余热功率吸收的热功率;/>为卡琳娜循环发电系统KLN的热电转换效率;
所述基线负荷约束为:
;
;
上式中,为综合能源系统IES未提供综合需求响应DR服务的综合能源系统IES基线负荷,对应综合能源系统IES的日常生产计划;
所述优化计算模块用于基于综合能源系统优化运行模型进行计算,输出优化后的综合能源系统运行策略。
实施例3:
参考图3,一种考虑等效循环寿命和温控负荷的综合能源运行设备,包括存储器和处理器;所述存储器,用于存储计算机程序代码,并将所述计算机程序代码传输给所述处理器;所述处理器,用于根据所述计算机程序代码中的指令执行实施例1所述的方法。
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现实施例1所述的方法。
Claims (6)
1.考虑等效循环寿命和温控负荷的综合能源运行方法,其特征在于:
所述综合能源运行方法包括以下步骤:
S1、构建综合能源系统优化运行模型;所述综合能源系统优化运行模型的目标函数为系统运行成本最小,所述系统运行成本为系统总成本与现货市场出清收益的差值,所述系统总成本包括能源购置成本、设备运行成本、需求响应成本、碳排放成本,其中,所述设备运行成本基于电储能系统EES的循环寿命损耗成本得到,所述需求响应成本基于温度调整成本、电动汽车EV的循环寿命损耗成本得到;所述综合能源系统优化运行模型的约束条件包括功率平衡约束、设备出力约束、基线负荷约束;
所述目标函数为:
min W=min(Wene+Wope+WDR+Wcarbon-Wearn);
上式中,W为系统运行成本;Wene、Wope、WDR、Wcarbon、Wearn分别为能源购置成本、设备运行成本、需求响应成本、碳排放成本、现货市场出清收益;cele,t、Ebuy,t分别为综合能源系统IES向电网购电的价格和功率;cgas、Gbuy,t分别为综合能源系统IES从气网购气的价格、气网输入综合能源系统IES的天然气功率;分别为热电联产机组CHP的启、停成本;分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB、溴化锂制冷系统LBRS、卡琳娜循环发电系统KLN、热储能系统TES、冷储能系统CES的运行成本;/>Wuser分别为需求响应前、需求响应后的用户购电成本;分别为冷、热、电、气四种能源对应的功率;/>为t时刻的综合需求响应DR资源出清价格;/>为t时刻综合能源系统IES第m段投标的中标量;T为规划期限;M为投标段数;Wadj为温度调整成本;/>为电动汽车EV参与综合需求响应DR调用的补偿成本;/>为电储能系统EES的循环寿命损耗成本;λc为碳交易基价;Ctr为碳排放交易额;αc为碳交易价格增长率;d为阶梯式碳交易区间长度;ωEV为电动汽车EV参与综合能源系统IES内综合需求响应DR调用的经济补偿系数;/>为电动汽车EV的循环寿命损耗成本;
所述电储能系统EES的循环寿命损耗成本电动汽车EV的循环寿命损耗成本/>基于以下公式计算得到:
上式中,mp、my、ms分别为电储能系统EES的单位功率投资成本、运维成本、单位容量投资成本;分别为电储能系统EES的额定功率、额定容量;me、/>分别为电动汽车EV电池单位容量成本、电动汽车EV额定容量;/>分别为电储能系统EES、电动汽车EV的电池寿命损耗成本;/>分别为电储能系统EES、电动汽车EV在100%放电深度下的等效循环次数;/>分别为电储能系统EES、电动汽车EV的线性化等效循环次数;
所述温度调整成本Wadj基于温控负荷等效热参数模型得到,所述温控负荷等效热参数模型为:
上式中,分别为室内温度下限、上限;/>Tdown,t、Tup,t分别为t时刻的室内设定温度、室内下调温度、室内上调温度;/>为t时刻的室内温度;Tadj为整日调整温度累计值;λ为温度调整惩罚系数;ω为惩罚增幅系数;l为单位区间;
S2、基于构建的综合能源系统优化运行模型进行计算,输出优化后的综合能源系统运行策略。
2.根据权利要求1所述的考虑等效循环寿命和温控负荷的综合能源运行方法,其特征在于:
所述功率平衡约束包括电功率平衡约束、热功率平衡约束、冷功率平衡约束、天然气平衡约束;其中,所述电功率平衡约束为:
上式中,ECHP,t、EPV,t、EKLN,t分别为热电联产机组CHP、光伏系统PV、卡琳娜循环发电系统KLN生产的电功率;ECACS,t、EEB,t、EIR,t分别为电锅炉CACS、电锅炉EB、工业制冷机IR消耗的电功率;EEESC,t、EEESD,t分别为电储能系统EES的充、放电功率;EEVC,t、EEVD,t为电动汽车EV的充、放电功率;Euser,t为t时刻的生活类负荷功率;Eload,t为t时刻的工业类负荷功率;
所述热功率平衡约束为:
HCHP,t+HGB,t+HTESD,t+HEB,t=Hload,t+HTESC,t;
上式中,HCHP,t、HGB,t、HEB,t分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB生产的热功率;HTESD,t、HTESC,t分别为热储能系统TES的充、放热功率;Hload,t为t时刻工业类负荷产生的冷功率;
所述冷功率平衡约束为:
LIR,t+LLBRS,t+LCESD,t=Lload,t+LCESC,t;
上式中,LIR,t、LLBRS,t分别为工业制冷机IR、溴化锂制冷系统LBRS生产的冷功率;LCESC,t、LCESD,t分别为电锅炉CACS充、放冷功率;Lload,t为t时刻工业类负荷产生的冷功率;
所述天然气平衡约束为:
Gbuy,t=GCHP,t+GGB,t+Gload,t;
上式中,GCHP,t、GGB,t分别为热电联产机组CHP、燃气锅炉GB耗费的天然气功率;Gload,t为t时刻工业类负荷产生的天然气功率;
所述设备出力约束包括热电联产机组约束、燃汽锅炉约束、余热锅炉约束、溴化锂制冷机约束、卡琳娜循环发电系统约束,所述热电联产机组约束为:
上式中,ECHP,0为热电联产机组CHP最大产热功率时对应的发电功率;ECHP,min、ECHP,max、ηCHP分别为热电联产机组CHP的最小发电功率、最大发电功率和发电效率;kdown、kup分别为热电联产机组CHP的热电耦合出力区间的下限斜率、上限斜率;uCHP,t为热电联产机组CHP的启停变量;Qgas为天然气热值;为热电联产机组CHP热电联产过程中产生的余热功率;HCHP,max为热电联产机组CHP最大产热功率;
所述燃汽锅炉约束为:
上式中,ηGB为燃气锅炉GB的产热效率;HGB,max为燃气锅炉GB的最大产热功率;
所述余热锅炉约束为:
上式中,为余热锅炉利用高温热气残余热量再加热后输出的热功率;ηYR,h为余热锅炉热能回收效率;/>为t时刻燃气锅炉GB产生的余热功率;
所述溴化锂制冷机约束为:
LLBRS,t=ηLBRSHLBRS,t;
上式中,HLBRS,t为溴化锂制冷系统LBRS从系统余热功率吸收的热功率;β为余热利用系数;ηLBRS为溴化锂制冷系统LBRS的热冷转换效率;HPV,t为t时刻光伏系统产热功率;HLBRS,t为t时刻溴化锂制冷系统产热功率;
所述卡琳娜循环发电系统约束为:
上式中,HKLN,t为卡琳娜循环发电系统KLN从系统余热功率吸收的热功率;ηKLN为卡琳娜循环发电系统KLN的热电转换效率;
所述基线负荷约束为:
上式中,为综合能源系统IES未提供综合需求响应DR服务的综合能源系统IES基线负荷,对应综合能源系统IES的日常生产计划。
3.考虑等效循环寿命和温控负荷的综合能源运行系统,其特征在于:
所述运行系统包括模型构建模块、优化计算模块;
所述模型构建模块用于构建综合能源系统优化运行模型;所述综合能源系统优化运行模型的目标函数为系统运行成本最小,所述系统运行成本为系统总成本与现货市场出清收益的差值,所述系统总成本包括能源购置成本、设备运行成本、需求响应成本、碳排放成本,其中,所述设备运行成本基于电储能系统EES的循环寿命损耗成本得到,所述需求响应成本基于温度调整成本、电动汽车EV的循环寿命损耗成本得到;所述综合能源系统优化运行模型的约束条件包括功率平衡约束、设备出力约束、基线负荷约束;
所述模型构建模块用于构建如下目标函数:
所述目标函数为:
minW=min(Wene+Wope+WDR+Wcarbon-Wearn);
上式中,W为系统运行成本;Wene、Wope、WDR、Wcarbon、Wearn分别为能源购置成本、设备运行成本、需求响应成本、碳排放成本、现货市场出清收益;cele,t、Ebuy,t分别为综合能源系统IES向电网购电的价格和功率;cgas、Gbuy,t分别为综合能源系统IES从气网购气的价格、气网输入综合能源系统IES的天然气功率;分别为热电联产机组CHP的启、停成本;分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB、溴化锂制冷系统LBRS、卡琳娜循环发电系统KLN、热储能系统TES、冷储能系统CES的运行成本;/>Wuser分别为需求响应前、需求响应后的用户购电成本;分别为冷、热、电、气四种能源对应的功率;/>为t时刻的综合需求响应DR资源出清价格;/>为t时刻综合能源系统IES第m段投标的中标量;T为规划期限;M为投标段数;Wadj为温度调整成本;/>为电动汽车EV参与综合需求响应DR调用的补偿成本;/>为电储能系统EES的循环寿命损耗成本;λc为碳交易基价;Ctr为碳排放交易额;αc为碳交易价格增长率;d为阶梯式碳交易区间长度;ωEV为电动汽车EV参与综合能源系统IES内综合需求响应DR调用的经济补偿系数;/>为电动汽车EV的循环寿命损耗成本;
所述电储能系统EES的循环寿命损耗成本电动汽车EV的循环寿命损耗成本/>基于以下公式计算得到:
上式中,mp、my、ms分别为电储能系统EES的单位功率投资成本、运维成本、单位容量投资成本;分别为电储能系统EES的额定功率、额定容量;me、/>分别为电动汽车EV电池单位容量成本、电动汽车EV额定容量;/>分别为电储能系统EES、电动汽车EV的电池寿命损耗成本;/>分别为电储能系统EES、电动汽车EV在100%放电深度下的等效循环次数;/>分别为电储能系统EES、电动汽车EV的线性化等效循环次数;
所述温度调整成本Wadj基于温控负荷等效热参数模型得到,所述温控负荷等效热参数模型为:
上式中,分别为室内温度下限、上限;/>Tdown,t、Tup,t分别为t时刻的室内设定温度、室内下调温度、室内上调温度;/>为t时刻的室内温度;Tadj为整日调整温度累计值;λ为温度调整惩罚系数;ω为惩罚增幅系数;l为单位区间;
所述优化计算模块用于基于综合能源系统优化运行模型进行计算,输出优化后的综合能源系统运行策略。
4.根据权利要求3所述的考虑等效循环寿命和温控负荷的综合能源运行系统,其特征在于:
所述功率平衡约束包括电功率平衡约束、热功率平衡约束、冷功率平衡约束、天然气平衡约束;其中,所述电功率平衡约束为:
上式中,ECHP,t、EPV,t、EKLN,t分别为热电联产机组CHP、光伏系统PV、卡琳娜循环发电系统KLN生产的电功率;ECACS,t、EEB,t、EIR,t分别为电锅炉CACS、电锅炉EB、工业制冷机IR消耗的电功率;EEESC,t、EEESD,t分别为电储能系统EES的充、放电功率;EEVC,t、EEVD,t为电动汽车EV的充、放电功率;Euser,t为t时刻的生活类负荷功率;Eload,t为t时刻的工业类负荷功率;
所述热功率平衡约束为:
HCHP,t+HGB,t+HTESD,t+HEB,t=Hload,t+HTESC,t;
上式中,HCHP,t、HGB,t、HEB,t分别为热电联产机组CHP、燃气锅炉GB、电锅炉EB生产的热功率;HTESD,t、HTESC,t分别为热储能系统TES的充、放热功率;Hload,t为t时刻工业类负荷产生的冷功率;
所述冷功率平衡约束为:
LIR,t+LLBRS,t+LCESD,t=Hload,t+LCESC,t;
上式中,LIR,t、LLBRS,t分别为工业制冷机IR、溴化锂制冷系统LBRS生产的冷功率;LCESC,t、LCESD,t分别为电锅炉CACS充、放冷功率;Lload,t为t时刻工业类负荷产生的冷功率;
所述天然气平衡约束为:
Gbuy,t=GCHP,t+GGB,t+Gload,t;
上式中,GCHP,t、GGB,t分别为热电联产机组CHP、燃气锅炉GB耗费的天然气功率;Gload,t为t时刻工业类负荷产生的天然气功率;
所述设备出力约束包括热电联产机组约束、燃汽锅炉约束、余热锅炉约束、溴化锂制冷机约束、卡琳娜循环发电系统约束,所述热电联产机组约束为:
上式中,ECHP,0为热电联产机组CHP最大产热功率时对应的发电功率;ECHP,min、ECHP,max、ηCHP分别为热电联产机组CHP的最小发电功率、最大发电功率和发电效率;kdown、kup分别为热电联产机组CHP的热电耦合出力区间的下限斜率、上限斜率;uCHP,t为热电联产机组CHP的启停变量;Qgas为天然气热值;为热电联产机组CHP热电联产过程中产生的余热功率;HCHP,max为热电联产机组CHP最大产热功率;
所述燃汽锅炉约束为:
上式中,ηGB为燃气锅炉GB的产热效率;HGB,max为燃气锅炉GB的最大产热功率;
所述余热锅炉约束为:
上式中,为余热锅炉利用高温热气残余热量再加热后输出的热功率;ηYR,h为余热锅炉热能回收效率;/>为t时刻燃气锅炉GB产生的余热功率;
所述溴化锂制冷机约束为:
LLBRS,t=ηLBRSHLBRS,t;
上式中,HLBRS,t为溴化锂制冷系统LBRS从系统余热功率吸收的热功率;β为余热利用系数;ηLBRS为溴化锂制冷系统LBRS的热冷转换效率;HPV,t为t时刻光伏系统产热功率;HLBRS,t为t时刻溴化锂制冷系统产热功率;
所述卡琳娜循环发电系统约束为:
上式中,HKLN,t为卡琳娜循环发电系统KLN从系统余热功率吸收的热功率;ηKLN为卡琳娜循环发电系统KLN的热电转换效率;
所述基线负荷约束为:
上式中,为综合能源系统IES未提供综合需求响应DR服务的综合能源系统IES基线负荷,对应综合能源系统IES的日常生产计划。
5.考虑等效循环寿命和温控负荷的综合能源运行设备,其特征在于:
所述运行设备包括存储器和处理器;
所述存储器,用于存储计算机程序代码,并将所述计算机程序代码传输给所述处理器;
所述处理器,用于根据所述计算机程序代码中的指令执行如权利要求1或2所述的方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1或2所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311754468.5A CN117436672B (zh) | 2023-12-20 | 2023-12-20 | 考虑等效循环寿命和温控负荷的综合能源运行方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311754468.5A CN117436672B (zh) | 2023-12-20 | 2023-12-20 | 考虑等效循环寿命和温控负荷的综合能源运行方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117436672A CN117436672A (zh) | 2024-01-23 |
CN117436672B true CN117436672B (zh) | 2024-03-12 |
Family
ID=89558596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311754468.5A Active CN117436672B (zh) | 2023-12-20 | 2023-12-20 | 考虑等效循环寿命和温控负荷的综合能源运行方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117436672B (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013039554A1 (en) * | 2011-09-16 | 2013-03-21 | Narayam Amit | A system and a method for optimization and management of demand response and distribute energy resources |
CN112365021A (zh) * | 2020-09-21 | 2021-02-12 | 华北电力大学 | 一种基于混合储能的区域综合能源系统规划优化方法 |
CN113517690A (zh) * | 2021-07-02 | 2021-10-19 | 东北电力大学 | 含电动汽车充电站的社区综合能源系统双层调度方法 |
CN113627720A (zh) * | 2021-06-30 | 2021-11-09 | 四川大学 | 基于电能替代的综合园区风光储系统配置方法 |
CN114022046A (zh) * | 2021-11-30 | 2022-02-08 | 国网四川省电力公司电力科学研究院 | 一种考虑交通均衡的综合能源系统优化调度方法 |
CN114066264A (zh) * | 2021-11-18 | 2022-02-18 | 国网山东综合能源服务有限公司 | 计及综合需求响应的综合能源系统运行优化方法及系统 |
CN114662330A (zh) * | 2022-03-31 | 2022-06-24 | 华北电力大学 | 计及碳交易机制和需求响应的综合能源系统模型构建方法 |
CN115438851A (zh) * | 2022-08-31 | 2022-12-06 | 南京理工大学 | 一种考虑极端情况的源网荷储协调运行方法 |
CN115577909A (zh) * | 2022-09-13 | 2023-01-06 | 四川大学 | 考虑价格型需求响应和v2g的园区综合能源系统调度方法 |
CN115879651A (zh) * | 2023-02-21 | 2023-03-31 | 国网天津市电力公司城西供电分公司 | 考虑电动汽车参与的综合能源系统低碳优化方法及装置 |
WO2023082697A1 (zh) * | 2021-11-15 | 2023-05-19 | 中国电力科学研究院有限公司 | 电-热综合能源系统协调优化方法、系统、设备、介质及程序 |
CN116167573A (zh) * | 2022-12-28 | 2023-05-26 | 国网河北省电力有限公司雄安新区供电公司 | 一种高铁站综合能源需求响应优化调度方法及相关装置 |
CN116316567A (zh) * | 2023-02-24 | 2023-06-23 | 国网甘肃省电力公司经济技术研究院 | 一种阶梯碳交易机制下综合能源需求响应优化方法 |
CN116542370A (zh) * | 2023-04-23 | 2023-08-04 | 华北电力大学 | 一种考虑碳捕集与碳交易的园区低碳经济运行方法 |
CN116822683A (zh) * | 2023-02-18 | 2023-09-29 | 南昌大学 | 一种基于碳交易计及不确定性的综合能源系统优化运行方法 |
-
2023
- 2023-12-20 CN CN202311754468.5A patent/CN117436672B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013039554A1 (en) * | 2011-09-16 | 2013-03-21 | Narayam Amit | A system and a method for optimization and management of demand response and distribute energy resources |
CN112365021A (zh) * | 2020-09-21 | 2021-02-12 | 华北电力大学 | 一种基于混合储能的区域综合能源系统规划优化方法 |
CN113627720A (zh) * | 2021-06-30 | 2021-11-09 | 四川大学 | 基于电能替代的综合园区风光储系统配置方法 |
CN113517690A (zh) * | 2021-07-02 | 2021-10-19 | 东北电力大学 | 含电动汽车充电站的社区综合能源系统双层调度方法 |
WO2023082697A1 (zh) * | 2021-11-15 | 2023-05-19 | 中国电力科学研究院有限公司 | 电-热综合能源系统协调优化方法、系统、设备、介质及程序 |
CN114066264A (zh) * | 2021-11-18 | 2022-02-18 | 国网山东综合能源服务有限公司 | 计及综合需求响应的综合能源系统运行优化方法及系统 |
CN114022046A (zh) * | 2021-11-30 | 2022-02-08 | 国网四川省电力公司电力科学研究院 | 一种考虑交通均衡的综合能源系统优化调度方法 |
CN114662330A (zh) * | 2022-03-31 | 2022-06-24 | 华北电力大学 | 计及碳交易机制和需求响应的综合能源系统模型构建方法 |
CN115438851A (zh) * | 2022-08-31 | 2022-12-06 | 南京理工大学 | 一种考虑极端情况的源网荷储协调运行方法 |
CN115577909A (zh) * | 2022-09-13 | 2023-01-06 | 四川大学 | 考虑价格型需求响应和v2g的园区综合能源系统调度方法 |
CN116167573A (zh) * | 2022-12-28 | 2023-05-26 | 国网河北省电力有限公司雄安新区供电公司 | 一种高铁站综合能源需求响应优化调度方法及相关装置 |
CN116822683A (zh) * | 2023-02-18 | 2023-09-29 | 南昌大学 | 一种基于碳交易计及不确定性的综合能源系统优化运行方法 |
CN115879651A (zh) * | 2023-02-21 | 2023-03-31 | 国网天津市电力公司城西供电分公司 | 考虑电动汽车参与的综合能源系统低碳优化方法及装置 |
CN116316567A (zh) * | 2023-02-24 | 2023-06-23 | 国网甘肃省电力公司经济技术研究院 | 一种阶梯碳交易机制下综合能源需求响应优化方法 |
CN116542370A (zh) * | 2023-04-23 | 2023-08-04 | 华北电力大学 | 一种考虑碳捕集与碳交易的园区低碳经济运行方法 |
Non-Patent Citations (3)
Title |
---|
Optimal Scheduling Method of Park Integrated Energy System Considering Carbon Transaction Cost;K. Ding, H. Ma, Y. Wu, Y. Liu, S. Liu and H. Dong;2022 IEEE 6th Conference on Energy Internet and Energy System Integration (EI2);20230510;504-509 * |
基于需求响应的区域微电网优化调度;徐斌;徐斌;刘红新;马骏;丁倩;;电力科学与技术学报;20180328(01);134-142 * |
工业园区综合能源系统低碳经济优化运行模型;李嘉祺 等;高电压技术;20220831;3190-3200 * |
Also Published As
Publication number | Publication date |
---|---|
CN117436672A (zh) | 2024-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | Optimal operation of regional integrated energy system considering demand response | |
CN110288152B (zh) | 考虑电/热柔性负荷的区域综合能源系统储能配置方法 | |
CN109659927B (zh) | 一种考虑储能参与度的综合能源微网储能容量配置方法 | |
CN109523052B (zh) | 一种考虑需求响应和碳交易的虚拟电厂优化调度方法 | |
CN105337303B (zh) | 一种含热泵的热电联产型微电网容量优化配置方法 | |
CN113095791B (zh) | 一种综合能源系统运行方法及系统 | |
CN108229025A (zh) | 一种冷热电联供型多微网主动配电系统经济优化调度方法 | |
CN110661254A (zh) | 一种区域综合能源系统冷热电互补效益量化方法 | |
CN108832656A (zh) | 基于电转气和可再生能源利用微能源网多目标规划方法 | |
CN105225022A (zh) | 一种热电联产型微电网的经济性优化运行方法 | |
CN112258021B (zh) | 用于家庭燃料电池热电联供建筑的能效评估方法及系统 | |
CN110826815A (zh) | 一种考虑综合需求响应的区域综合能源系统运行优化方法 | |
CN111416388B (zh) | 一种高比例新能源系统储能需求和储能时长确定方法 | |
CN113762708A (zh) | 一种考虑多目标协同的园区级综合能源系统规划方法 | |
CN111738503A (zh) | 以氢能为核心的综合能源微网日前运行调度方法及系统 | |
CN111382902A (zh) | 基于运行效益增量的区域综合能源系统储能优化配置方法 | |
Han et al. | Economic evaluation of micro-grid system in commercial parks based on echelon utilization batteries | |
CN113344249A (zh) | 基于区块链的冷热电联供多微网优化调度方法和系统 | |
CN115864376A (zh) | 一种多时间尺度电力保供分析方法及系统 | |
CN116341847A (zh) | 一种变工况下综合能源系统碳排放的调度方法及系统 | |
CN115358519A (zh) | 一种虚拟电厂优化调度方法及装置 | |
CN115099007A (zh) | 基于综合成本-能耗曲线的综合能源系统优化运行方法 | |
CN112883630B (zh) | 用于风电消纳的多微网系统日前优化经济调度方法 | |
CN112421642B (zh) | 一种综合能源系统ies可靠性评估方法和系统 | |
CN113885367A (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 |