CN114676897A - 一种含chp-p2g-氢能的园区综合能源系统优化调度方法 - Google Patents
一种含chp-p2g-氢能的园区综合能源系统优化调度方法 Download PDFInfo
- Publication number
- CN114676897A CN114676897A CN202210260006.7A CN202210260006A CN114676897A CN 114676897 A CN114676897 A CN 114676897A CN 202210260006 A CN202210260006 A CN 202210260006A CN 114676897 A CN114676897 A CN 114676897A
- Authority
- CN
- China
- Prior art keywords
- power
- chp
- unit
- energy
- cost
- 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.)
- Withdrawn
Links
- 239000001257 hydrogen Substances 0.000 title claims abstract description 133
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 20
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 151
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 107
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 84
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 84
- 239000003345 natural gas Substances 0.000 claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 239000007789 gas Substances 0.000 claims description 81
- 238000003860 storage Methods 0.000 claims description 56
- 230000009194 climbing Effects 0.000 claims description 26
- 238000006243 chemical reaction Methods 0.000 claims description 19
- 230000005611 electricity Effects 0.000 claims description 19
- 238000012423 maintenance Methods 0.000 claims description 16
- 238000004519 manufacturing process Methods 0.000 claims description 14
- 230000007613 environmental effect Effects 0.000 claims description 10
- 230000008901 benefit Effects 0.000 claims description 8
- 230000007423 decrease Effects 0.000 claims description 8
- 238000004146 energy storage 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
- 150000002431 hydrogen Chemical class 0.000 claims description 7
- 238000010248 power generation Methods 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 7
- RZGZMLICFFEUIQ-UHFFFAOYSA-N 5-[(1-phenylcyclohexyl)amino]pentanoic acid Chemical compound C=1C=CC=CC=1C1(NCCCCC(=O)O)CCCCC1 RZGZMLICFFEUIQ-UHFFFAOYSA-N 0.000 claims description 4
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims 1
- 230000002457 bidirectional effect Effects 0.000 abstract description 5
- 238000010206 sensitivity analysis Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 9
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 5
- 229910002092 carbon dioxide Inorganic materials 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004148 unit process Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007728 cost analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005868 electrolysis reaction Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/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
- 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
- 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)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种含CHP‑P2G‑氢能的园区综合能源系统优化调度方法,包括以下步骤:步骤1、构建园区综合能源系统;步骤2、构建考虑CHP‑P2G‑氢能的园区综合能源系统优化调度模型;步骤3、构建优化调度模型的约束条件;步骤4、以总成本最小为目标,在MATLAB环境下,以总成本最小为目标,通过Yalmip工具箱调用Cplex求解器进行求解;步骤5、对运行成本和碳排放成本进行灵敏度分析,设置不同的偏好度。本发明可以实现电力网络与天然气网络能量的双向流动,通过电解槽将不确定性的风能转化为氢气,提升风电消纳水平,同时实现能量高效利用,同时权重系数的设置对调度部门具有一定的参考价值。
Description
技术领域
本发明属于园区综合能源系统优化调度技术领域,具体涉及一种含CHP-P2G-氢能的园区综合能源系统优化调度方法以及对系统运行成本和碳排放成本偏好度的设置范围。
背景技术
随着化石能源危机逐渐加剧和日益严重的环境污染问题,完成能源转型和构建现代化能源体系已迫在眉睫。由于风电等可再生能源发电具有较为明显的间歇性和不确定性,造成弃风现象严重。氢作为一种清洁能源,可直接高效利用,例如为氢燃料汽车或某些化学过程提供原料。近年来,随着综合能源系统这一概念的提出,其与可再生能源的相互耦合将极大提高能源利用效率,减少污染排放。因此,多能耦合、多荷交叉的综合能源系统是未来能源供应体系的重要发展趋势。电力网络与天然气网络是当前能源领域最主要的大规模传输载体,因此两者之间的耦合关系一直受到广泛关注。通常借助燃气轮机等设备将天然气转化为电能,而电转气(Power to Gas,P2G)技术可以将电能转化为氢气或天然气,进一步加强了电力网络与天然气网络之间的耦合,实现了能量之间的双向流动,为消纳风电等可再生能源提供了一个新的方案。同时,为提高可再生能源的利用率,通常需要加入储能装置,与其他储能方式相比,储氢在清洁绿色等方面具有明显优势。
现有研究大多为电转天然气,但是由于电转氢气的效率更高,所以其为解决可再生能源消纳问题的重要途径,可以从源头上实现零碳或低碳。为实现能量的高效利用,需合理充分的调控P2G两阶段运行过程。我国三北地区风力资源丰富,但由于其具有随机性以及反调峰特性,如何合理调度使用已成为当务之急。
发明内容
为解决上述技术问题,本发明提出一种含CHP-P2G-氢能的园区综合能源系统优化调度方法,该方法充分利用电制氢过程,实现P2G效率最大化,构建的系统通过热电联产(CHP)机组和P2G,实现电力网络与天然气网络能量的双向流动,优先利用电制氢将具有不确定性和反调峰特性的风能转化为氢能供给氢负荷,提升风电消纳水平,减少碳排放,实现能量高效利用,同时对运行成本和碳排放成本进行灵敏度分析,设置不同偏好度,可为调度部门提供参考,具有一定的实践意义。
本发明采取的技术方案为:
一种含CHP-P2G-氢能的园区综合能源系统优化调度方法及成本偏好度设置,包括以下步骤:
步骤1、构建园区综合能源系统:
该系统的能源侧包括:风能、上级电网、上级气网;
能量转换设备包含热电联产CHP机组、电解槽、甲烷反应器和燃气锅炉;
负荷侧包括:电负荷、氢负荷、热负荷;
CHP机组主要利用甲烷反应器产生的天然气,不足部分从上级气网购得,产生电能和热能供给负荷,电解槽将风电产生的波动性较大的电能转化为氢能,一部分经储氢罐供给氢负荷,其余部分传输至甲烷反应器,甲烷反应器将氢气经与二氧化碳反应生成天然气供气供给CHP机组,形成一个循环。
步骤2、构建考虑CHP-P2G-氢能的园区综合能源系统优化调度模型:
以经济性最优和碳排放最低为目标,引入碳税ctax,将碳排放量转化为碳排放成本,以园区运行成本F1和碳排放成本F2之和最小为目标函数,如式(1)所示:
式(1)中,minF为总目标函数最小,eCHP为CHP机组的碳排放量,eGB为燃气锅炉的碳排放量。
(i)运行成本F1:
运行成本F1包括购能成本cbuy和运维成本co,如式(2)所示。其中,购能成本cbuy包含购电成本和购气成本,运维成本co为各设备在t时段运行维护成本之和。
式(2)中,ce和cg分别为单位购电、购气价格;Pbuye和Pbuyg分别为t时段从上级网络中购买的电量与天然气量;Ki为设备的单位可变运行维护成本;Pi为设备在t时段的出力。
(ii)碳排放成本F2:
碳排放成本F2包括CHP机组碳排放量eCHP和燃气锅炉碳排放量eGB两部分,如式(3)所示。式(3)中,eCHP为CHP机组碳排放量cCHP和碳排放配额cCHPb之差。
PCHP=PCHPe+γPCHPh
式(3)、式(4)中,μGB为燃气锅炉单位CO2排放强度;PGBh为燃气锅炉的产热功率;PCHP为CHP机组在t时刻折算到纯凝工况下的电功率;PCHPe为CHP机组的产电功率;
PCHPh为CHP机组的产热功率;cCHPb为CHP机组的碳排放配额;γ为固定进气量下,CHP机组增加单位热功率时电功率的减少值;α1,α2,α3分别为CHP机组碳排放量系数;β为单位电量碳交易配额。
步骤3、构建优化调度模型的约束条件:
(a)功率平衡约束:
式(5)中,PWT为风电出力;Peup和Pgup分别为系统从上级网络的购电量和购气量;Peload、Phload和Pgload分别为电负荷、热负荷和氢负荷。PCHPe为CHP机组的产电功率;PEC为电解槽的耗电功率;PMR为甲烷反应器的输出功率;PCHPg为CHP机组的耗气功率;PGBg为燃气锅炉的耗气功率;PCHPh为CHP机组的产热功率;PGBh为燃气锅炉的产热功率。
(b)电转气部分约束:
优化调度模型的P2G部分涉及约束如式(6)所示,包含电解槽、甲烷反应器和储氢罐耦合部分。
式(6)中,PEC为电解槽的输出功率;PMRi为甲烷反应器的输入功率;PHSTi为储氢罐的输入功率;PHSTo为储氢罐的输出功率。
(c)储能设备约束:
储能设备为储氢罐,其约束如式(7)所示:
式(7)中,和为储氢容量上、下限;为储氢罐输入功率上限;为储氢罐输入功率下限;为储氢罐输出功率上限;为储氢罐输出功率下限;EHST,t为储氢罐在t时刻的容量;PHSTi为储氢罐在t时刻的输入功率;PHSTo为储氢罐在t时刻的输出功率;ET为储氢罐在末端时刻的容量;E1为储氢罐在初始时刻的容量。
(d)机组设备约束:
该部分包含设备出力上、下限约束,如式(8)所示,以及可控机组爬坡约束,如式(9)所示。
式(8-9)中,为CHP机组电出力上限,为CHP机组电出力下限,为CHP机组热出力上限,为CHP机组热出力下限,为燃气锅炉输入功率上限,为燃气锅炉输入功率下限,为电解槽输入功率上限,为电解槽输入功率下限,为甲烷反应器输入功率上限,为甲烷反应器输入功率下限,为CHP机组的爬坡功率上限,为CHP机组的爬坡功率下限,为燃气锅炉的爬坡功率上限,为燃气锅炉的爬坡功率下限,为电解槽的爬坡功率上限,为电解槽的爬坡功率下限,为甲烷反应器的爬坡功率上限,为甲烷反应器的爬坡功率下限,PCHPe为CHP机组的输出电功率,PCHPh为CHP机组的输出热功率,PGBg为燃气锅炉的耗气功率,PECi为电解槽的输入功率,PMRi为甲烷反应器的输入功率,为CHP机组在t时刻折算到纯凝工况下的电功率,为CHP机组在(t-1)时刻折算到纯凝工况下的电功率,PGBh,t为燃气锅炉在t时刻的产热功率,PGBh,t-1为燃气锅炉在(t-1)时刻的产热功率,PECi,t+1为电解槽在(t+1)时刻的输入功率,PECi,t为电解槽在t时刻的输入功率,PMRi,t+1为甲烷反应器在(t+1)时刻的输入功率,PMRi,t为甲烷反应器在t时刻的输入功率。
步骤4、以总成本最小为目标,通过Cplex求解器进行求解。
本发明一种含CHP-P2G-氢能的园区综合能源系统优化调度方法,技术效果如下:
1)本发明能够有效降低碳排放量,实现电能、氢能、天然气之间的能量流动,充分考虑P2G两阶段转化效率的不同,将风能转化成系统所需氢气,提高系统经济性,实现了P2G中能量的高效利用,具有较优的经济效益和环境效益。
2)氢气不易运输与储存,而本发明创新性的将P2G设备拆分成两阶段分步利用,优先将第一阶段产生的氢气通过储氢罐直接供给氢负荷,不仅减少了经长管拖车运输的危险性,同时也能更加快速灵活的利用氢能,如若有多余,则进行后续的甲烷化反应。
3)通过引入权重系数对系统进行灵敏度分析,形成不同情境下的优化配置方案对比,不同的优化偏好设置会显著影响系统配置结果。本发明所研究的综合能源系统合理调度系数范围为(r1=0.7,r2=0.3)~(r1=0.4,r2=0.6),未来在优化过程中通过考虑对不同方面的关注程度,可选取合适的权重系数,从而实现对经济指标和环境指标的灵活调度控制。
附图说明
图1为电转气运行过程图。
图2为本发明所建立园区综合能源系统的能量流向图。
图3为典型日负荷及风电预测数据曲线图。
图4(1)为场景一机组出力图;
图4(2)为场景二机组出力图;
图4(3)为场景三机组出力图。
图5为电解槽、甲烷反应器出力图。
图6为储氢罐容量变化示意图。
图7为不同权重系数下经济指标与环境指标对比图。
具体实施方式
原理分析:
电转气包含电转氢气和电转天然气两部分,其中电转氢气效率可达75%~85%,电转天然气效率为45%~60%。其两步化学方程式如下所示。
CO2+4H2→CH4+2H2O
电解槽通过利用风能电解水产生氢气和氧气,其中,一部分氢气注入储氢罐,供给氢负荷,另一部分未被利用的氢气则通入甲烷反应器,产生的CH4可供CHP机组和燃气锅炉使用。本发明系统所涉及电转气运行过程图如图1所示。
为充分消纳风电,充分利用P2G两阶段运行过程,本发明构建了包含风电、天然气、氢能的园区综合能源系统,建立了考虑电转气的园区综合能源系统优化调度模型,以经济性最优和碳排放最低为目标函数,提高了园区能量供给的自主性,并通过算例分析验证了该模型所具备的良好性能,最后引入权重系数,对算例进行灵敏度分析与比较,且每种情况下对应的经济和环境因素也不同,从而选取最优区间范围,为决策者提供参考。包括以下步骤:
步骤1、构建园区综合能源系统:
本发明所构建的园区综合能源系统,如图2所示。通过热电联产(CHP)机组和P2G可以实现电力网络与天然气网络能量的双向流动,通过电解槽将不确定性的风能转化为氢气,提升风电消纳水平,同时实现能量高效利用。
该系统的能源侧包括:风能、上级电网、上级气网;
能量转换设备包含热电联产CHP机组、电解槽、甲烷反应器和燃气锅炉;
负荷侧包括:电负荷、氢负荷、热负荷;
CHP机组主要利用甲烷反应器产生的天然气,不足部分从上级气网购得,产生电能和热能供给负荷,电解槽将风电产生的稳定性较差的电能转化为氢气,一部分经储氢罐供给氢负荷,其余部分传输至甲烷反应器,甲烷反应器将氢气经与二氧化碳反应生成天然气供气供给CHP机组,形成一个循环。
步骤2、构建考虑P2G的园区综合能源系统优化调度模型:
为实现经济性最优和碳排放最低,本发明通过引入碳税ctax,将碳排放量转化为碳排放成本,从而使多目标问题转变为单目标优化问题进行求解。
以园区运行成本F1和碳排放成本F2之和最小为目标函数,如式(1)所示:
minF=F1+F2
F1=cbuy+co
式(1)中,minF为总目标函数最小;
eCHP为CHP机组的碳排放量;
eGB为燃气锅炉的碳排放量。
(i)运行成本F1:
运行成本F1包括购能成本cbuy和运维成本co,如式(2)所示。其中,购能成本cbuy包含购电成本和购气成本,运维成本co为各设备在t时段运行维护成本之和。
式(2)中,ce和cg分别为单位购电、购气价格;
Pbuye和Pbuyg分别为t时段从上级网络中购买的电量与天然气量;
Ki为设备的单位可变运行维护成本;
Pi为设备在t时段的出力。
(ii)碳排放成本F2:
碳排放成本F2包括CHP机组碳排放成本eCHP和燃气锅炉碳排放成本eGB两部分,如式(3)所示。式(3)中,eCHP为CHP机组碳排放量cCHP和碳排放配额cCHPb之差。
PCHP=PCHPe+γPCHPh
式(3-4)中,μGB为燃气锅炉单位CO2排放强度;
PGBh为燃气锅炉产热功率;
PCHP为CHP机组在t时刻折算到纯凝工况下的电功率;
PCHPe为CHP机组的输出电功率;
PCHPh为CHP机组的输出热功率;
cCHPb为CHP机组的碳排放配额;
γ为固定进气量下,CHP机组增加单位热功率时电功率的减少值;
α1,α2,α3分别为CHP机组碳排放量系数;
β为单位电量碳交易配额。
步骤3、构建优化调度模型的约束条件:
(a)功率平衡约束:
式(5)中,PWT为风电出力;
Peup和Pgup分别为系统从上级网络的购电量和购气量;
Peload、Phload和Pgload分别为电负荷、热负荷和氢负荷。
PCHPe为CHP机组的输出电功率;
PEC为电解槽的耗电功率;
PMR为甲烷反应器的输出功率;
PCHPg为CHP机组的耗气功率;
PGBg为燃气锅炉的耗气功率;
PCHPh为CHP机组的输出热功率;
PGBh为燃气锅炉的输出功率。
(b)电转气部分约束:
优化调度模型的P2G部分涉及约束如式(6)所示,包含电解槽、甲烷反应器和储氢罐耦合部分。
PEC=PMRi+PHSTi
PHSTo=Pgload (6);
式(6)中,PEC为t时刻电解槽的输出功率;
PMRi为甲烷反应器的输入功率;
PHSTi为储氢罐的输入功率;
PHSTo为储氢罐的输出功率。
(c)储能设备约束:
储能设备为储氢罐,其约束如式(7)所示:
EHST,t为储氢罐在t时刻的容量;
PHSTi为储氢罐在t时刻的输入功率;
PHSTo为储氢罐在t时刻的输出功率;
ET为储氢罐在末端时刻的容量;
E1为储氢罐在初始时刻的容量。
(d)机组设备约束:
该部分包含设备出力上、下限约束,如式(8)所示,以及可控机组爬坡约束,如式(9)所示。
PCHPe为CHP机组的输出电功率;
PCHPh为CHP机组的输出热功率;
PGBg为燃气锅炉的耗气功率;
PECi为电解槽的输入功率;
PMRi为甲烷反应器的输入功率;
PGBh,t为燃气锅炉在t时刻的产热功率;
PGBh,t-1为燃气锅炉在(t-1)时刻的产热功率
PECi,t+1为电解槽在(t+1)时刻的输入功率;
PECi,t为电解槽在t时刻的输入功率;
PMRi,t+1为甲烷反应器在(t+1)时刻的输入功率;
PMRi,t为甲烷反应器在t时刻的输入功率。
步骤4、在MATLAB环境中进行求解:
本发明在MATLAB环境下,以总成本最小为目标,通过Yalmip工具箱调用Cplex求解器进行求解。
步骤5、在不同场景下进行仿真验证:
本发明以我国某省某工业园区为研究对象,对模型进行验证。甘肃省风电装机容量占总装机容量的25.08%以上,市场环境成熟,配套政策优越,适合风电项目的发展。为实现能源供应的经济性和环保性,构建了天然气-风能-氢能园区综合能源系统。
为验证包含P2G的园区综合能源系统与电、热、气、氢各系统独立运行时相比所具有的优势,选取以下三种场景进行分析,如表1所示:
表1场景设置
表1中,“√”代表考虑,“×”代表不考虑。
单位调度时间为1h,电价采用0.5元/(kW·h),天然气单位价格取3元/m3。各类设备功能计划以满足园区负荷为目的。图3展示了该园区典型日负荷预测数据和风电出力基础数据从图3可以看出,从14时和17时为用电高峰期,9时和17时热负荷较高,夜间氢负荷较高,风电出力较大。
各种设备运行参数如表2所示
表2各设备运行参数
表3各场景成本
(1)表3为各场景下成本对比分析。由表3可知,场景1~场景3的总成本依次降低,场景3较场景1降低了24.8%,说明在考虑电转氢气以及加入甲烷反应器的过程中,可以最大发挥系统的经济性。场景1购能成本最高,主要原因为在不包含P2G时,氢负荷供给全部依靠购买所得,而氢气购买价格较高,所以导致购能成本较高,场景2购能成本次之,因为其所需一部分氢气可由风电电解水得到,剩余部分则需从外界购买,场景3够能成本最低,主要原因是其购氢成本较低。场景1~场景3碳排放成本依次降低,说明逐步考虑P2G后,系统的减排效果越来越好。
(2)图4(1)~图4(3)为各场景机组出力图。对于电负荷,在1:00~5:00时,由于负荷需求较低,而此时风电出力较大,场景3中,为有效降低弃风,电解槽在夜间出力升高,以达到消纳弃风和供给此时较高氢负荷的作用,7:00之后,热负荷增大,为满足需求,CHP机组出力增加,同时购电量增大,在17:00时,热负荷较高,CHP机组和燃气锅炉出力较大,在20:00之后用电需求降低,热负荷略有起伏,此时CHP机组出力减少,燃气锅炉出力增加,满足负荷要求。对于氢负荷,场景1和2完全依赖于购氢,价格较为昂贵。对于场景3,在6:00之后电解槽出力随着氢气需求趋于平缓而逐渐下降,在13:00时,由于负荷需求降低,此时电解槽的出力也较低。
(3)本发明构建的含风电-氢能-P2G的园区IES的电解槽、甲烷反应器出力如图5所示,储氢罐容量变化如图6所示。
图5可知,在1:00~7:00时,由于风力资源丰富,且此时氢负荷较大,使得电解槽出力较高,提供给甲烷反应器的氢气减少,在7:00~13:00时,风电出力下降,电解槽输出功率也随之下降,由于其产生的氢气一部分要供给氢负荷,所以甲烷反应器的出力降低,在13:00后风电开始逐渐升高,前期输出功率下降的甲烷反应器电也开始随着电解槽输出功率的上升而逐渐上升,夜间22:00~24:00时,风电资源丰富,此时氢负荷需求处于较低水平,所以甲烷反应器输出功率增加,产生的天然气供CHP机组和燃气锅炉使用。
由图6可知,储氢罐容量在1:00~5:00时逐步下降,是由于此时氢负荷需求较高,而电解槽运行产生氢气需要一段时间,所以消耗了储氢罐原有氢气,在5:00~10:00时,由于氢能需求缓慢下降,此时风电和电解槽出力还在较高水平,所以储氢罐容量逐渐上升,在10:00~18:00时,风电出力较小,电解槽出力较低,供给氢负荷导致储氢罐容量降低,在18:00~24:00时,随着风电增加,电解槽产生的氢气增多,储氢罐容量随之上升。
(4)灵敏度分析:
对系统运行成本F1和碳排放成本F2分别赋予不同权重r1、r2,分析其对系统的影响。其中,权重系数用于表示被给予的关注度,且r1+r2=1。当r1>r2时,意味着对系统的运行成本,即经济性更加关注;当r2>r1时,说明对系统的碳排放成本,即环保性给予了更多的关注。图7为不同权重系数下经济指标、环境指标及总成本的变化曲线。
不同的系数设置会影响最终的优化结果。当r2越大,即决策者更关注环保性时,碳排放成本显著降低,亦则反之。从图7可以得出,当权重系数从极端关注经济效益(r1=1,r2=0)到逐渐关注经济效益(r1=0.7,r2=0.3)时,随着对经济性关注度的减少,运行成本逐渐增加,碳排放成本降低。然而,当过度关注碳排放时,将导致运行成本大幅增加,但此时碳排放的减少非常有限。因此,在坐标轴右侧的系数范围并没有实际参考价值,性价比非常低。综上,合理系数选取范围在(r1=0.7,r2=0.3)~(r1=0.4,r2=0.6)之间。
本发明建立了天然气-风能-氢能优化模型。综合能源系统IES主要由风电、CHP机组、燃气锅炉和P2G组成,所涉及的能量流包括天然气流、电力流、热力流和氢气流。该模型具体包括设备数学建模、目标函数和约束条件。结合算例仿真,得出电解槽、甲烷反应器、储氢罐、CHP机组的组合实现了电、氢气、天然气之间的能量的双向流动,借助P2G两阶段运行,考虑其不同阶段的转换效率,实现了能量更加高效的利用。
Claims (4)
1.一种含CHP-P2G-氢能的园区综合能源系统优化调度方法,其特征在于包括以下步骤:
步骤1、构建园区综合能源系统:
该系统的能源侧包括:风能、上级电网、上级气网;
能量转换设备包含热电联产CHP机组、电解槽、甲烷反应器和燃气锅炉;
负荷侧包括:电负荷、氢负荷、热负荷;
步骤2、构建考虑CHP-P2G-氢能的园区综合能源系统优化调度模型:
以经济性最优和碳排放最低为目标,引入碳税ctax,将碳排放量转化为碳排放成本,以园区运行成本F1和碳排放成本F2之和最小为目标函数,如式(1)所示:
式(1)中,min F为总目标函数最小,eCHP为CHP机组的碳排放量,eGB为燃气锅炉的碳排放量;
(i)运行成本F1:
运行成本F1包括购能成本cbuy和运维成本co,如式(2)所示;其中,购能成本cbuy包含购电成本和购气成本,运维成本co为各设备在t时段运行维护成本之和;
式(2)中,ce和cg分别为单位购电、购气价格;Pbuye和Pbuyg分别为t时段从上级网络中购买的电量与天然气量;Ki为设备的单位可变运行维护成本;Pi为设备在t时段的出力;
(ii)碳排放成本F2:
碳排放成本F2包括CHP机组碳排放量eCHP和燃气锅炉碳排放量eGB两部分,如式(3)所示;式(3)中,eCHP为CHP机组碳排放量cCHP和碳排放配额cCHPb之差;
PCHP=PCHPe+γPCHPh
式(3)、式(4)中,μGB为燃气锅炉单位CO2排放强度;PGBh为燃气锅炉的产热功率;PCHP为CHP机组在t时刻折算到纯凝工况下的电功率;PCHPe为CHP机组的产电功率;
PCHPh为CHP机组的产热功率;cCHPb为CHP机组的碳排放配额;γ为固定进气量下,CHP机组增加单位热功率时电功率的减少值;α1,α2,α3分别为CHP机组碳排放量系数;β为单位电量碳交易配额;
步骤3、构建优化调度模型的约束条件:
(a)功率平衡约束:
式(5)中,PWT为风电出力;Peup和Pgup分别为系统从上级网络的购电量和购气量;Peload、Phload和Pgload分别为电负荷、热负荷和氢负荷;PCHPe为CHP机组的产电功率;PEC为电解槽的耗电功率;PMR为甲烷反应器的输出功率;PCHPg为CHP机组的耗气功率;PGBg为燃气锅炉的耗气功率;PCHPh为CHP机组的产热功率;PGBh为燃气锅炉的产热功率;
(b)电转气部分约束:
优化调度模型的P2G部分涉及约束如式(6)所示,包含电解槽、甲烷反应器和储氢罐耦合部分;
式(6)中,PEC为电解槽的输出功率;PMRi为甲烷反应器的输入功率;PHSTi为储氢罐的输入功率;PHSTo为储氢罐的输出功率;
(c)储能设备约束:
储能设备为储氢罐,其约束如式(7)所示:
式(7)中,和为储氢容量上、下限;为储氢罐输入功率上限;为储氢罐输入功率下限;为储氢罐输出功率上限;为储氢罐输出功率下限;EHST,t为储氢罐在t时刻的容量;PHSTi为储氢罐在t时刻的输入功率;PHSTo为储氢罐在t时刻的输出功率;ET为储氢罐在末端时刻的容量;E1为储氢罐在初始时刻的容量;
(d)机组设备约束:
该部分包含设备出力上、下限约束,如式(8)所示,以及可控机组爬坡约束,如式(9)所示;
式(8-9)中,为CHP机组电出力上限,为CHP机组电出力下限,为CHP机组热出力上限,为CHP机组热出力下限,为燃气锅炉输入功率上限,为燃气锅炉输入功率下限,为电解槽输入功率上限,为电解槽输入功率下限,为甲烷反应器输入功率上限,为甲烷反应器输入功率下限,为CHP机组的爬坡功率上限,为CHP机组的爬坡功率下限,为燃气锅炉的爬坡功率上限,为燃气锅炉的爬坡功率下限,为电解槽的爬坡功率上限,为电解槽的爬坡功率下限,为甲烷反应器的爬坡功率上限,为甲烷反应器的爬坡功率下限,PCHPe为CHP机组的输出电功率,PCHPh为CHP机组的输出热功率,PGBg为燃气锅炉的耗气功率,PECi为电解槽的输入功率,PMRi为甲烷反应器的输入功率,为CHP机组在t时刻折算到纯凝工况下的电功率,为CHP机组在(t-1)时刻折算到纯凝工况下的电功率,PGBh,t为燃气锅炉在t时刻的产热功率,PGBh,t-1为燃气锅炉在(t-1)时刻的产热功率,PECi,t+1为电解槽在(t+1)时刻的输入功率,PECi,t为电解槽在t时刻的输入功率,PMRi,t+1为甲烷反应器在(t+1)时刻的输入功率,PMRi,t为甲烷反应器在t时刻的输入功率;
步骤4、以总成本最小为目标,通过Cplex求解器进行求解。
2.根据权利要求1所述一种含CHP-P2G-氢能的园区综合能源系统优化调度方法,其特征在于:对系统运行成本F1和碳排放成本F2分别赋予不同权重r1、r2,分析其对系统的影响:
其中,权重系数用于表示被给予的关注度,且r1+r2=1;
当r1>r2时,意味着对系统的运行成本,即经济性更加关注;
当r2>r1时,说明对系统的碳排放成本,即环保性给予了更多的关注。
3.根据权利要求2所述一种含CHP-P2G-氢能的园区综合能源系统优化调度方法,其特征在于:不同的系数设置会影响最终的优化结果;当r2越大,即决策者更关注环保性时,碳排放成本显著降低,亦则反之;当权重系数从极端关注经济效益(r1=1,r2=0)到逐渐关注经济效益(r1=0.7,r2=0.3)时,随着对经济性关注度的减少,运行成本逐渐增加,碳排放成本降低;然而,当过度关注碳排放时,将导致运行成本大幅增加,但此时碳排放的减少非常有限;合理系数选取范围在(r1=0.7,r2=0.3)~(r1=0.4,r2=0.6)之间。
4.考虑CHP-P2G-氢能的园区综合能源系统优化调度模型,其特征在于:
该模型以经济性最优和碳排放最低为目标,引入碳税ctax,将碳排放量转化为碳排放成本,以园区运行成本F1和碳排放成本F2之和最小为目标函数,如式(1)所示:
式(1)中,min F为总目标函数最小,eCHP为CHP机组的碳排放量,eGB为燃气锅炉的碳排放量;
(i)运行成本F1:
运行成本F1包括购能成本cbuy和运维成本co,如式(2)所示;其中,购能成本cbuy包含购电成本和购气成本,运维成本co为各设备在t时段运行维护成本之和;
式(2)中,ce和cg分别为单位购电、购气价格;Pbuye和Pbuyg分别为t时段从上级网络中购买的电量与天然气量;Ki为设备的单位可变运行维护成本;Pi为设备在t时段的出力;
(ii)碳排放成本F2:
碳排放成本F2包括CHP机组碳排放量eCHP和燃气锅炉碳排放量eGB两部分,如式(3)所示;式(3)中,eCHP为CHP机组碳排放量cCHP和碳排放配额cCHPb之差;
PCHP=PCHPe+γPCHPh
式(3)、式(4)中,μGB为燃气锅炉单位CO2排放强度;PGBh为燃气锅炉的产热功率;PCHP为CHP机组在t时刻折算到纯凝工况下的电功率;PCHPe为CHP机组的产电功率;PCHPh为CHP机组的产热功率;cCHPb为CHP机组的碳排放配额;γ为固定进气量下,CHP机组增加单位热功率时电功率的减少值;α1,α2,α3分别为CHP机组碳排放量系数;β为单位电量碳交易配额。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210260006.7A CN114676897A (zh) | 2022-03-16 | 2022-03-16 | 一种含chp-p2g-氢能的园区综合能源系统优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210260006.7A CN114676897A (zh) | 2022-03-16 | 2022-03-16 | 一种含chp-p2g-氢能的园区综合能源系统优化调度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114676897A true CN114676897A (zh) | 2022-06-28 |
Family
ID=82073666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210260006.7A Withdrawn CN114676897A (zh) | 2022-03-16 | 2022-03-16 | 一种含chp-p2g-氢能的园区综合能源系统优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114676897A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115619006A (zh) * | 2022-09-23 | 2023-01-17 | 河海大学 | 一种计及辅助服务的电-气-氢混联综合能源系统优化调度方法 |
-
2022
- 2022-03-16 CN CN202210260006.7A patent/CN114676897A/zh not_active Withdrawn
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115619006A (zh) * | 2022-09-23 | 2023-01-17 | 河海大学 | 一种计及辅助服务的电-气-氢混联综合能源系统优化调度方法 |
CN115619006B (zh) * | 2022-09-23 | 2023-12-05 | 河海大学 | 一种计及辅助服务的电-气-氢混联综合能源系统优化调度方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108154309B (zh) | 计及冷热电多负荷动态响应的能源互联网经济调度方法 | |
CN111639824B (zh) | 一种含电转气的区域综合能源系统热电优化调度方法 | |
CN108832656B (zh) | 基于电转气和可再生能源利用微能源网多目标规划方法 | |
CN113315242B (zh) | 基于氢能经济促进弃风消纳的虚拟弃风-制氢联合体 | |
CN108665188B (zh) | 一种基于优化模型的园区多能源主体综合配比方法 | |
CN113850474B (zh) | 一种热电氢多能流综合能源系统及其优化调度方法 | |
CN112990523A (zh) | 基于多目标模型预测控制的区域综合能源系统分层优化运行方法 | |
CN113078687B (zh) | 孤岛多能互补电-气耦合系统能量优化调度方法 | |
CN113159407A (zh) | 基于区域综合能源系统的多能储存模块容量优化配置方法 | |
CN110518865A (zh) | 一种风电-光电耦合的工业规模稳定供氢系统 | |
CN116093949A (zh) | 考虑两阶段p2g混合储能与碳势控制需求响应优化方法 | |
CN110957722A (zh) | 一种含电转气设备的微型能源网日前优化调度方法 | |
CN114676897A (zh) | 一种含chp-p2g-氢能的园区综合能源系统优化调度方法 | |
Li et al. | Planning model of integrated energy system considering P2G and energy storage | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法 | |
CN116596123A (zh) | 一种园区综合能源系统低碳优化调度方法 | |
CN116205458A (zh) | 考虑不确定性的水电制氢的综合能源系统优化调度模型建立方法 | |
CN115495906A (zh) | 高比例可再生能源接入的配网多能互补协同优化方法 | |
CN114330835A (zh) | 一种综合能源微网中电/热混合储能系统的优化配置方法 | |
Li et al. | Low-carbon economic dispatch of integrated energy system considering p2g and carbon trading mechanism | |
Liu et al. | Day-ahead economic dispatch of integrated energy system including power to gas | |
CN114362152B (zh) | 综合能源系统多时间尺度调度方法 | |
Tao et al. | Research on Power Storage Optimization Operation Strategy for Wind-Photovoltaic Power Coupled Hydrogen Production | |
CN109375588B (zh) | 一种考虑广义需求响应的工厂多能互补优化调度方法 | |
Qiu et al. | IES Low-Carbon Economic Dispatch Considering Natural Gas Pressure Energy Power Generation under P2G-CCS |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220628 |