CN114156948A - 工业园区氢能综合利用站的能量管理方法 - Google Patents

工业园区氢能综合利用站的能量管理方法 Download PDF

Info

Publication number
CN114156948A
CN114156948A CN202210027215.7A CN202210027215A CN114156948A CN 114156948 A CN114156948 A CN 114156948A CN 202210027215 A CN202210027215 A CN 202210027215A CN 114156948 A CN114156948 A CN 114156948A
Authority
CN
China
Prior art keywords
formula
power
time
hydrogen
output
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
Application number
CN202210027215.7A
Other languages
English (en)
Other versions
CN114156948B (zh
Inventor
毕锐
王孝淦
袁华凯
朱正轩
吴红斌
王磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202210027215.7A priority Critical patent/CN114156948B/zh
Publication of CN114156948A publication Critical patent/CN114156948A/zh
Application granted granted Critical
Publication of CN114156948B publication Critical patent/CN114156948B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • 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
    • 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/28—Arrangements for balancing of the load in a network by storage of energy
    • 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/28—Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • 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]
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • 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
    • Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B10/00—Integration of renewable energy sources in buildings
    • Y02B10/10—Photovoltaic [PV]
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/30—Hydrogen technology
    • Y02E60/36—Hydrogen production from non-carbon containing sources, e.g. by water electrolysis
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00—Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin
    • 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
    • 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/20—Climate change mitigation technologies for sector-wide applications using renewable energy
    • 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
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/40—Fuel cell technologies in production processes
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Educational Administration (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Fuel Cell (AREA)

Abstract

本发明公开了一种工业园区氢能综合利用站的能量管理方法，其步骤包括，1，确定工业园区综合利用电站的主要设备，搭建工业园区氢能综合利用站的系统结构图；2，建立电解槽、储氢罐、燃料电池、燃气锅炉和风光的数学模型，以及电母线和热母线的功率平衡约束方程；3，以工业园区碳排放量最小为约束，计及风光在不同时段的出力，考虑氢能综合利用站的工作状态和能量管理方案。本发明能够更加合理的对氢能综合利用站的能流进行控制，在保持供能可靠性的条件下，优化能量流动和充分实现能量的梯级利用，使工业园区处于可靠经济的供能运行状态，从而能够提高能量的利用效率。

Description

工业园区氢能综合利用站的能量管理方法
技术领域
本发明涉及工业园区综合能源能量管理技术领域，具体说是一种工业园区氢能综合利用站能量管理策略。
背景技术
随着能源危机的逐渐加剧，在双碳政策大力推广和能源互联网不断发展的背景下，以电能为核心的能源网络结构持续向低碳化转型，分布式能源发电越来越受到重视。但是由于新能源出力的不确定性和随机性，出现严重的可再生能源的弃风弃光现象并且对系统的稳定提出了很大的挑战。氢能作为新一代清洁能源，其对可再生能源的消纳能力提高和实现电力系统向低碳化转型有着不可忽视的作用。工业园区作为低碳化转型的重要对象，对工业园区氢能综合利用站进行能量管理提高能源的利用效率势在必行。
目前，工业园区的能量管理主要集中于对分布式电源、储能的能量控制策略，。尽管有研究考虑氢能在工业园区中能量管理的影响，但是其主要考虑燃料电池产电的影响，很少考虑燃料电池产生的余热。目前，由于光伏和风电的边际成本进一步降低，未来工业园区内的供电将主要是光风可再生能源，为了进一步提高风光的消纳能力，未来氢能将会广泛推广应用于工业园区中。然而仅仅利用氢能综合系统的供电能力提高风光的消纳能力，无法将氢能综合利用系统利用最大化。
发明内容
本发明为促进工业园区向低碳转型，提供一种工业园区氢能综合利用站的能量管理方法，以期能够更加合理的使用氢能和可再生能源，实现工业园区内供能的可靠以及能量的梯级利用，优化工业园区的功率流动和提高能源的利用率，从而使工业园区处于高效经济运行状态。
为实现上述目的，本发明采用以下技术方案：
本发明一种工业园区氢能综合利用站的能量管理方法的特点包括以下步骤：
1)，搭建工业园区氢能综合利用站的系统结构，包括：
制氢过程设备：电解槽；
储氢设备：储氢罐，用于为燃料电池提供氢能，同时作为加氢站为氢燃料汽车提供燃料；
用氢设备：燃料电池，包括：固体氧化物燃料电池和质子交换膜燃料电池；
燃气锅炉：用于提供热能；
由光伏风电和外电网提供电功率；
2)，基于工业园区氢能综合利用站的系统结构，建立氢能综合利用站的精细化数学模型，包括：电解槽、燃料电池、储氢罐、燃气锅炉和光伏风电的数学模型；
3)，考虑风光出力不同时段的氢能综合利用站的工作状态，制定能量管理方案。
本发明所述的一种工业园区氢能综合利用站的能量管理方法的特点在于，所述步骤2包括：
步骤2.1、建立电解槽的数学模型：
步骤2.1.1利用式(1)得到电解槽的电压特性：
Figure BDA0003464565410000021
式(1)中，UEL(t)是t时刻电解槽的电压；urev是电解槽的开路电压；r1、r2分别为电解槽的欧姆特性参数，s1、s2、s3、s4、s5和s6分别为电解槽的电极过电压参数；TEL为电解槽的工作温度；AEL为电解槽的面积；IEL(t)为t时刻电解槽的直流电流；
步骤2.1.2利用式(2)和式(3)分别建立电解槽的效率和功率数学模型：
Figure BDA0003464565410000022
式(2)中，a1～a5为五个法拉第效率系数；IEL为直流电流；
PEL_H2＝PELηEL (3)
式(3)中，PEL为电解槽从电母线电解制氢所需的输入功率；PEL_H2为电解槽输出的氢功率；ηEL为电解槽的效率；
步骤2.1.2利用式(4)建立电解槽的产氢速率模型：
Figure BDA0003464565410000023
式(4)中，nEL为产氢的速率；HHV为氢气的高热值；PEL_H2为电解槽输出的氢功率；
步骤2.2建立质子交换膜燃料电池的数学模型；
步骤2.2.1利用式(5)建立单电池的输出电压Vcell的数学模型：
Vcell＝Vn-Vact-Vohm-Vcon (5)
式(5)中，Vn、Vact、Vohm、Vcon分别为能斯特电压、活化极化过电压、欧姆极化过电压和浓差极化过电压；
步骤2.2.2利用式(6)～式(9)建立能斯特电压、活化极化过电压、欧姆极化过电压损耗和浓差极化过电压的数学模型：
利用式(6)建立燃料电池能斯特电压的数学模型：
Figure BDA0003464565410000031
式(6)中，Vn为能斯特电压；TFC、
Figure BDA0003464565410000032
分别为燃料电池的工作温度、氢气分压、氧气分压和水蒸气分压；
利用式(7)建立活化极化过电压的数学模型：
Figure BDA0003464565410000033
式(7)中，Vact为活化极化电压；ξ1、ξ2、ξ3、ξ4分别为四个经验参数；
Figure BDA0003464565410000034
J分别为电池的电流密度和气液界面的氧溶解浓度；
Figure BDA0003464565410000035
为输入的氧气分压；
利用式(8)建立欧姆极化过电压的数学模型：
Figure BDA0003464565410000036
式(8)中，Vohm为欧姆极化过电压；Rm为质子通过膜电极的等效阻抗；Rc电子通过极板的等效阻抗；ρ为膜阻率；β为膜的厚度；A为膜的有效面积；
利用式(9)建立浓差极化过电压的数学模型：
Figure BDA0003464565410000037
式(9)中，Vcon为浓差极化过电压；λ为温度相关系数；ω为气体在催化层中的增长率；
步骤2.2.3建立质子交换膜燃料电池的输出功率模型：
利用式(10)得到燃料电池的工作效率ηFC：
Figure BDA0003464565410000038
式(10)中，μ为燃料利用率；Vcell为电池的输出电压；
利用式(11)建立燃料电池的输出电功率和热功率数学模型：
Figure BDA0003464565410000041
式(11)中，PFC为燃料电池输出的电功率；
Figure BDA0003464565410000042
为储氢罐输出的功率；QFC为燃料电池输出的热功率；
利用式(12)建立燃料电池的耗氢速率数学模型：
Figure BDA0003464565410000043
式(12)中，nFC为燃料电池的耗氢速率；LHV为氢气的低热值；
Figure BDA0003464565410000044
为储氢罐输出的功率；
步骤2.3建立光伏风电以及燃气锅炉的数学模型；
步骤2.3.1利用式(13)建立光伏系统的数学模型：
Figure BDA0003464565410000045
式(13)中，Cpv为光伏系统的额定功率；S(t)为t时刻的辐照度；Sstd为标准辐照度；k为功率温度系数；Tc(t)为t时刻的光伏发电板表面温度；Tref为参考温度；
步骤2.3.2利用式(14)建立风电的数学模型：
Figure BDA0003464565410000046
式(14)中，PWT(t)为风电在t时刻的出力；PN为风机的额定功率；v(t)为t时刻的风速；vin(t)为t时刻风机的切入风速；vN为额定功率时的风速；vout风机的切出风速；
步骤2.3.3通过式(15)得到燃气锅炉的出力模型：
QGB(t)＝ηGBV(t)LV (15)
式(15)中，ηGB为燃气锅炉的效率；V(t)为t时刻消耗的天然气的体积；LV为天然气的低热值；QGB(t)为t时刻燃气锅炉的输出功率；
步骤2.4利用式(16)建立储氢罐的数学模型：
Figure BDA0003464565410000051
式(16)中，PHT(t)为t时刻储氢罐的压强；EHT(t)为t时刻的储氢量；R为阿伏伽德罗常数；
Figure BDA0003464565410000052
为标准开氏温度；VHT为储氢罐的体积；nEL(t)为t时刻电解槽的产氢量；nFC(t)为t时刻燃料电池的耗氢量；EHT(t-1)为t-1时刻的储氢量；Δt为时间间隔；PN为额定压强；
步骤2.5利用式(17)建立氢能综合利用站的碳排放数学模型：
Figure BDA0003464565410000053
式(17)中，βe、βh分别表示单位电量、天然气量对应的碳排放；Pe,graid(t)表示t时刻向电网购买的电量；V(t)表示t时刻消耗的天然气体积；LV表示天然气的低热值；
步骤2.6在氢能综合利用站中，当电解槽工作时，风光发电和市电向电解槽供电，从而利用式(18)得到电母线t时刻的功率平衡约束：
PPV(t)+PWT(t)+δPgrid(t)-Pload(t)＝PEL(t) (18)
式(18)中，PPV(t)为t时刻的光伏出力；PWT(t)为t时刻的风机出力；Pgrid(t)为t时刻与电网交互功率；Pload(t)为t时刻的负荷功率；PEL(t)为t时刻电解槽的输入功率；δ表示0-1变量；
利用式(19)得到热母线的平衡约束：
QGB(t)＝Qload(t) (19)
式(19)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；
步骤2.7在氢能综合利用站中，当电解槽不工作，由燃料电池、风光发电和市电向电负荷供能；由燃气锅炉和燃料电池向热负荷供能，从而利用式(20)得到电母线功率平衡约束：
PPV(t)+PWT(t)+δPgrid(t)+PFC(t)＝Pload(t) (20)
式(20)中，PFC(t)为t时刻燃料电池输出的电功率；
利用式(23)得到热母线平衡约束：
δQGB(t)+QFC(t)＝Qload(t) (21)
式(23)中，QFC(t)为燃料电池输出的热功率；
步骤2.8，利用式(22)建立以碳排量
Figure BDA0003464565410000061
最小为约束的目标函数：
Figure BDA0003464565410000062
当风光处于出力较大时段时，按如下步骤制定氢能综合利用站工作状态的能量管理方案：
步骤3.1a，利用式(23)判断氢能综合利用站的工作状态，若满足式(23)，则执行步骤3.2a-步骤3.5a，否则，执行步骤4.1；
PPV(t)+PWT(t)≥Pload(t) (23)
步骤3.2a，利用式(24)判断储氢罐是否满足约束，若满足约束，则执行步骤3.3a，否则，电解槽不工作：
SOHT(t)≤SOHTmax (24)
式(24)中，SOHT(t)为t时刻的储氢罐状态；SOHTmax表示储氢罐最大储氢状态量；
步骤3.3a，利用式(25)得到t时刻电解槽的输入功率PEL(t)：
Figure BDA0003464565410000063
式(25)中，CEL为电解槽的容量；
Figure BDA0003464565410000064
为储氢罐的最大储氢容量；EHT(t)为t时刻储氢罐的储氢量；
步骤3.4a，利用式(26)得到满足热负荷的t时刻燃气锅炉出力QGB(t)：
QGB(t)＝Qload(t) (26)
式(26)中，Qload(t)为t时刻热负荷的功率；
步骤3.5a，利用式(27)得到工业园区与外电网购电的交互功率Pgrid(t)：
Pgrid(t)+PPV(t)+PWT(t)＝Pload(t)+PEL(t) (27)。
当风光出力处于较小时段时，按如下步骤制定氢能综合利用站工作状态的能量管理方案：
步骤3.1b，利用式(28)判断储氢罐是否满足约束，若满足约束，则执行步骤3.2b，否则，燃料电池不工作；
SOHT(t)≤SOHTmax (28)
式(29)中，SOHT(t)为t时刻的储氢罐状态；SOHTmax表示储氢罐最大储氢量。
步骤3.2b，利用式(29)得到t时刻燃料电池输出的电功率PFC(t)：
Figure BDA0003464565410000071
式(30)中，CFC为燃料电池的容量；
Figure BDA0003464565410000072
为储氢罐的最大储氢容量；EHT(t)为t时刻储氢罐的储氢量；ηFC为燃料电池的工作效率；
步骤3.3b，利用式(30)得到t时刻燃料电池的输出热功率QFC(t)：
Figure BDA0003464565410000073
式(31)中，PFC(t)为t时刻燃料电池输出的电功率；
步骤3.4b，利用式(31)得到t时刻燃气锅炉输出热功率QGB(t)：
Figure BDA0003464565410000074
式(32)中，Qload(t)为t时刻热负荷的功率；QFC(t)为t时刻燃料电池输出的热功率；
步骤3.5b，利用式(32)得到t时刻工业园区与外电网的交互功率Pgrid(t)：
Pgrid(t)＝max{Pload(t)-PPV(t)-PWT(t)-PFC(t),0} (32)
步骤3.6b，利用式(33)得到工业园区所需要的热能：
δQGB(t)+QFC(t)＝Qload(t) (33)
式(33)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；QFC(t)为t时刻燃料电池输出的热功率。
与现有技术相比，本发明的有益效果为在于：
(1)本发明提出工业园区氢能综合利用站的系统结构，其包括光伏风电单元、燃气锅炉与氢储能系统；在高比例可再生能源渗透的背景下，通过考虑风光出力特性，提出氢能综合利用站的能量管理方案，通过配置氢储能系统有效提高了风光的消纳能力，并在满足工业园区安全、可靠运行的基础上，实现了工业园区内能量的合理流动以及设备的有效运行。
(2)本发明的能量管理方法利用燃料电池的热电耦合特性，充分利用其余热满足热负荷的需求，实现了氢能综合利用站的能量梯级利用；在计及风光出力特性的基础上，实现了氢能综合利用站在风光出力不同时段的灵活运行，并通过约束电解槽和燃料电池的运行状态，进一步提高了氢能综合利用站的能源利用效率。
附图说明
图1为本发明工业园区氢能综合利用站的系统结构图；
图2为本发明工业园区氢能综合利用站的能量管理方法的流程图；
图3为本发明工业园区氢能综合利用站的能量管理方法的示意图。
具体实施方式
本实施例中，如图2所示，一种工业园区氢能综合利用站的能量管理方法包括以下步骤；
1)，搭建工业园区氢能综合利用站的系统结构框图；
如图1所示工业园区存在制氢、储氢和用氢为一体的氢能综合能源一体站，制氢过程设备为电解槽，由风光可再生能源和外电网提供电功率，由燃气锅炉提供热能；储氢设备为储氢罐，可以为燃料电池提供氢能亦可以作为加氢站为氢燃料汽车提供燃料；用氢设备为燃料电池，常用的燃料电池主要为固体氧化物燃料电池(SOFC)和质子交换膜燃料电池(PEMFC)。
2)，基于氢能综合利用站的系统结构框图，建立氢能综合利用站设备的精细化数学模型，其中包括电解槽、燃料电池、储氢罐、燃气锅炉和光伏风电的数学模型；
步骤2.1、建立电解槽的数学模型：
步骤2.1.1利用式(1)得到电解槽的电压特性为：
Figure BDA0003464565410000081
式(1)中，UEL(t)是t时刻电解槽的电压；urev是电解槽的开路电压；r1、r2分别为电解槽的欧姆特性参数，s1、s2、s3、s4、s5和s6分别为电解槽的电极过电压参数；TEL为电解槽的工作温度；AEL为电解槽的面积；IEL(t)为t时刻电解槽的直流电流；
步骤2.1.2利用式(2)、式(3)建立电解槽的效率和功率数学模型：
Figure BDA0003464565410000091
式(2)中，a1～a5为法拉第效率系数；IEL为直流电流；AEL为电解槽面积。
PEL_H2＝PELηEL (3)
式(3)中，PEL为电解槽从电母线电解制氢所需的输入功率；PEL_H2为电解槽输出的氢功率；ηEL为电解槽的效率。
步骤2.1.2利用式(4)建立电解槽的产氢速率模型：
Figure BDA0003464565410000092
式(4)中，nEL为产氢的速率；HHV为氢气的高热值；PEL_H2为电解槽输出氢功率。
步骤2.2建立质子交换膜燃料电池的数学模型；
步骤2.2.1利用式(5)建立单电池的输出电压Vcell的数学模型：
Vcell＝Vn-Vact-Vohm-Vcon (5)
式(5)中，Vn、Vact、Vohm、Vcon分别为能斯特电压、活化极化过电压损耗、欧姆极化过电压损耗和浓差极化过电压损耗。
步骤2.2.2利用式(6)～(9)建立能斯特电压、活化极化过电压损耗、欧姆极化过电压损耗和浓差极化过电压损耗数学模型：
利用式(6)建立燃料电池能斯特电压的数学模型：
Figure BDA0003464565410000093
式(6)中，Vn为能斯特电压；TFC、
Figure BDA0003464565410000094
分别为燃料电池工作温度、氢气分压、氧气分压和水蒸气分压，由于氢燃料电池中的产物为液态水，
Figure BDA0003464565410000095
为1.
利用式(7)建立活化极化过电压数学模型：
Figure BDA0003464565410000096
式(7)中，Vact为活化极化电压；ξ1、ξ2、ξ3、ξ4分别为经验参数；
Figure BDA0003464565410000101
J分别为电池的电流密度和气液界面的氧溶解浓度；
Figure BDA0003464565410000102
为输入的氧气分压；TFC为燃料电池工作温度。
利用公式(8)建立欧姆极化过电压数学模型：
Figure BDA0003464565410000103
式(8)中，Vohm为欧姆极化过电压；Rm为质子通过膜电极的等效阻抗；Rc电子通过极板的等效阻抗；ρ为膜阻率；β为膜的厚度；A为膜的有效面积。由于Rm远远大于Rc，因此忽略Rc。
利用式(9)建立浓差极化过电压数学模型：
Vcon＝λeωJ (9)
式(9)中，Vcon为浓差极化过电压；λ为温度相关系数；ω为气体在催化层中的增长率。
步骤2.2.3建立质子交换膜燃料电池的输出功率模型：
由于氢燃料电池由多个单电池组合串联而成，利用式(10)可以得到燃料电池的利用效率：
Figure BDA0003464565410000104
式(10)中，ηfc为燃料电池工作效率；HHV为氢气的高热值；μ为燃料利用率；Vcell为电池的输出电压。
利用式(11)建立燃料电池的输出电功率和热功率数学模型：
Figure BDA0003464565410000105
式(11)中，PFC为燃料电池输出的电功率；
Figure BDA0003464565410000106
为储氢罐输出的功率；ηFC为燃料电池的工作效率；QFC为燃料电池输出的热功率。
利用式(12)建立燃料电池的耗氢速率数学模型：
Figure BDA0003464565410000107
式(12)中，nFC为燃料电池的耗氢速率；LHV为氢气的低热值；
Figure BDA0003464565410000108
为储氢罐输出的功率。
步骤2.3建立风光可再生能源以及燃气锅炉的数学模型模型；
步骤2.3.1利用式(13)建立光伏的数学模型：
Figure BDA0003464565410000111
式(13)中，Cpv为光伏系统的额定功率；S(t)为t时刻的辐照度；Sstd标准辐照度；k为功率温度系数；Tc(t)为t时刻的光伏发电板表面温度；Tref为参考温度。
步骤2.3.2利用式(14)建立风电的数学模型：
Figure BDA0003464565410000112
式(14)中，PWT(t)为风电在t时刻的出力；PN为风机的额定功率；v(t)为t时刻的风速；vin(t)为t时刻风机的切入风速；vN为额定功率时的风速；vout风机的切出风速。
步骤2.3.3通过式(15)得到燃气锅炉的出力模型：
QGB(t)＝ηGBV(t)LV (15)
式(15)中，ηGB为燃气锅炉的效率；V(t)为t时刻消耗的天然气的体积；LV为天然气的低热值；QGB(t)为t时刻燃气锅炉输出功率。
步骤2.4利用式(16)建立氢储能的数学模型：
Figure BDA0003464565410000113
式(16)中，PHT(t)为t时刻储氢罐压强；EHT(t)为t时刻储氢量；R为阿伏伽德罗常数；
Figure BDA0003464565410000114
为标准开氏温度；VHT为储氢罐的体积；nEL(t)为t时刻电解槽的产氢量；nFC(t)为燃料电池耗氢量；EHT(t-1)为t-1时刻的储氢量；Δt为时间间隔；PN为额定压强。
步骤2.5利用式(17)建立氢能综合利用站的碳排放数学模型：
Figure BDA0003464565410000121
式(17)中，βe、βh分别表示单位电量、天然气量对应的碳排放；Pe,graid(t)表示t时刻向电网购买的电量；V(t)表示t时刻消耗的天然气体积；LV表示天然气的低热值；Δt表示时间间隔。
步骤2.6为避免燃料电池和电解槽同时工作时燃料电池产生的电能制氢降低氢储能系统的工作效率。氢能综合利用站有两种工作方式，一种是当电解槽开始工作时，此时风光发电单元和市电向电解槽供电，利用式(18)可得到电母线此时的功率平衡约束：
PPV(t)+PWT(t)+δPgrid(t)-Pload(t)＝PEL(t) (18)
式(18)中，PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pgrid(t)为t时刻与电网交互功率；Pload(t)为t时刻负荷功率；PEL(t)为t时刻电解槽的输入功率；δ表示0-1变量。
此时，燃料电池不工作，利用式(19)可以得到热母线的平衡约束为：
QGB(t)＝Qload(t) (19)
式(19)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率。
步骤2.7另一种为燃料电池供能模式，此时电解槽不工作，此时电负荷由燃料电池、风光发电单元和市电参与供能；即热负荷则由燃气锅炉和燃料电池供能，此时电母线功率平衡约束为：
PPV(t)+PWT(t)+δPgrid(t)+PFC(t)＝Pload(t) (20)
式(20)中，PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pgrid(t)为t时刻与电网交互功率；Pload(t)为t时刻负荷功率；PFC(t)为t时刻燃料电池输出的电功率；δ表示0-1变量。
利用式(21)热母线平衡约束为：
δQGB(t)+QFC(t)＝Qload(t) (21)
式(21)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；δ表示0-1变量；QFC(t)为燃料电池输出热功率。
步骤2.8，在考虑工业园区氢能综合利用站的能量管理方案时，要以碳排量
Figure BDA0003464565410000122
最小为约束，即：
Figure BDA0003464565410000131
3)，考虑风光不同出力情况下的氢能综合利用站的工作状态，其中包括：主要考虑风光出力较大时段、风光出力较小时段的氢能综合利用站工作状态，制定能量管理方案，其过程如图3所示。
当风光处于出力较大时段时，考虑氢能综合利用站工作状态的能量管理方案为：
步骤3.1，利用式(23)判断氢能综合利用站的工作状态，若满足则执行步骤3.2-步骤3.5，否则执行步骤4.1；
PPV(t)+PWT(t)≥Pload(t) (23)
式(23)中，PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pload(t)为t时刻负荷功率。步骤3.2，利用式(24)判断储氢罐是否满足约束，若满足则继续执行：否则电解槽不工作；
SOHT(t)≤SOHTmax (24)
式(24)中，SOHT(t)为t时刻的储氢罐状态；SOHTmax表示储氢罐最大储氢状态量。
步骤3.3，利用式(25)得到此时电解槽的输入功率：
Figure BDA0003464565410000132
式(25)中，PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pload(t)为t时刻负荷功率；PEL(t)为t时刻电解槽的输入电功率；CEL为电解槽容量；
Figure BDA0003464565410000133
为储氢罐最大储氢容量；EHT(t)为t时刻储氢罐储氢量；ηEL为电解槽工作效率；Δt为时间间隔。
步骤3.4，利用式(26)可以得到满足热负荷的燃气锅炉出力为：
QGB(t)＝Qload(t) (26)
式(26)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率。
步骤3.5，利用式(27)得到工业园区与外电网购电的交互功率为：
Pgrid(t)+PPV(t)+PWT(t)＝Pload(t)+PEL(t) (27)
式(27)中，PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pload(t)为t时刻负荷功率；PEL(t)为t时刻电解槽的输入电功率；Pgrid(t)为t时刻与电网交互功率。
当风光出力处于较小时段时，考虑氢能综合利用站工作状态的能量管理方案为：
步骤4.1，利用式(29)判断储氢罐是否满足约束；否则燃料电池不工作；
SOHT(t)≤SOHTmax (28)
式(28)中，SOHT(t)为t时刻的储氢罐状态；SOHTmax表示储氢罐最大储氢量。
步骤4.2，此时燃料电池工作而电解槽不工作，利用式(30)可以得到燃料电池输出电功率为：
Figure BDA0003464565410000141
式(29)中，PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pload(t)为t时刻负荷功率；PFC(t)为t时刻电解槽的输入电功率；CFC为燃料电池的容量；
Figure BDA0003464565410000142
为储氢罐最大储氢量；EHT(t)为t时刻储氢罐储氢量；ηFC为燃料电池工作效率；Δt为时间间隔。
步骤4.3，利用式(31)可以得到燃料电池的输出热功率为：
Figure BDA0003464565410000143
式(30)中，ηFC为燃料电池的工作效率；QFC(t)为t时刻燃料电池输出的热功率；PFC(t)为t时刻燃料电池输出电功率。
步骤4.3，利用式(31)可以得到燃气锅炉输出热功率：
Figure BDA0003464565410000144
式(31)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；QFC(t)为t时刻燃料电池输出的热功率。
步骤4.4,利用式(32)可以得到工业园区与外电网的交互功率为：
Pgrid(t)＝max{Pload(t)-PPV(t)-PWT(t)-PFC(t),0} (32)
式(32)中，PFC(t)为t时刻燃料电池输出电功率；PPV(t)为t时刻光伏出力；PWT(t)为t时刻风机出力；Pload(t)为t时刻负荷功率；Pgrid(t)为t时刻与电网交互功率。
步骤4.5,工业园区所需要的热能由燃气锅炉和燃料电池提供：
δQGB(t)+QFC(t)＝Qload(t) (33)
式(33)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；QFC(t)为t时刻燃料电池输出的热功率；δ为0-1变量。

Claims (4)

1.一种工业园区氢能综合利用站的能量管理方法，其特征包括以下步骤：
1)，搭建工业园区氢能综合利用站的系统结构，包括：
制氢过程设备：电解槽；
储氢设备：储氢罐，用于为燃料电池提供氢能，同时作为加氢站为氢燃料汽车提供燃料；
用氢设备：燃料电池，包括：固体氧化物燃料电池和质子交换膜燃料电池；
燃气锅炉：用于提供热能；
由光伏风电和外电网提供电功率；
2)，基于工业园区氢能综合利用站的系统结构，建立氢能综合利用站的精细化数学模型，包括：电解槽、燃料电池、储氢罐、燃气锅炉和光伏风电的数学模型；
3)，考虑风光出力不同时段的氢能综合利用站的工作状态，制定能量管理方案。
2.根据权利要求1所述的一种工业园区氢能综合利用站的能量管理方法，其特征在于，所述步骤2包括：
步骤2.1、建立电解槽的数学模型：
步骤2.1.1利用式(1)得到电解槽的电压特性：
Figure FDA0003464565400000011
式(1)中，UEL(t)是t时刻电解槽的电压；urev是电解槽的开路电压；r1、r2分别为电解槽的欧姆特性参数，s1、s2、s3、s4、s5和s6分别为电解槽的电极过电压参数；TEL为电解槽的工作温度；AEL为电解槽的面积；IEL(t)为t时刻电解槽的直流电流；
步骤2.1.2利用式(2)和式(3)分别建立电解槽的效率和功率数学模型：
Figure FDA0003464565400000012
式(2)中，a1～a5为五个法拉第效率系数；IEL为直流电流；
PEL_H2＝PELηEL (3)
式(3)中，PEL为电解槽从电母线电解制氢所需的输入功率；PEL_H2为电解槽输出的氢功率；ηEL为电解槽的效率；
步骤2.1.2利用式(4)建立电解槽的产氢速率模型：
Figure FDA0003464565400000021
式(4)中，nEL为产氢的速率；HHV为氢气的高热值；PEL_H2为电解槽输出的氢功率；
步骤2.2建立质子交换膜燃料电池的数学模型；
步骤2.2.1利用式(5)建立单电池的输出电压Vcell的数学模型：
Vcell＝Vn-Vact-Vohm-Vcon (5)
式(5)中，Vn、Vact、Vohm、Vcon分别为能斯特电压、活化极化过电压、欧姆极化过电压和浓差极化过电压；
步骤2.2.2利用式(6)～式(9)建立能斯特电压、活化极化过电压、欧姆极化过电压损耗和浓差极化过电压的数学模型：
利用式(6)建立燃料电池能斯特电压的数学模型：
Figure FDA0003464565400000022
式(6)中，Vn为能斯特电压；TFC、
Figure FDA0003464565400000023
分别为燃料电池的工作温度、氢气分压、氧气分压和水蒸气分压；
利用式(7)建立活化极化过电压的数学模型：
Figure FDA0003464565400000024
式(7)中，Vact为活化极化电压；ξ1、ξ2、ξ3、ξ4分别为四个经验参数；
Figure FDA0003464565400000025
J分别为电池的电流密度和气液界面的氧溶解浓度；
Figure FDA0003464565400000026
为输入的氧气分压；
利用式(8)建立欧姆极化过电压的数学模型：
Figure FDA0003464565400000027
式(8)中，Vohm为欧姆极化过电压；Rm为质子通过膜电极的等效阻抗；Rc电子通过极板的等效阻抗；ρ为膜阻率；β为膜的厚度；A为膜的有效面积；
利用式(9)建立浓差极化过电压的数学模型：
Figure FDA0003464565400000031
式(9)中，Vcon为浓差极化过电压；λ为温度相关系数；ω为气体在催化层中的增长率；
步骤2.2.3建立质子交换膜燃料电池的输出功率模型：
利用式(10)得到燃料电池的工作效率ηFC：
Figure FDA0003464565400000032
式(10)中，μ为燃料利用率；Vcell为电池的输出电压；
利用式(11)建立燃料电池的输出电功率和热功率数学模型：
Figure FDA0003464565400000033
式(11)中，PFC为燃料电池输出的电功率；
Figure FDA0003464565400000034
为储氢罐输出的功率；QFC为燃料电池输出的热功率；
利用式(12)建立燃料电池的耗氢速率数学模型：
Figure FDA0003464565400000035
式(12)中，nFC为燃料电池的耗氢速率；LHV为氢气的低热值；
Figure FDA0003464565400000036
为储氢罐输出的功率；
步骤2.3建立光伏风电以及燃气锅炉的数学模型；
步骤2.3.1利用式(13)建立光伏系统的数学模型：
Figure FDA0003464565400000037
式(13)中，Cpv为光伏系统的额定功率；S(t)为t时刻的辐照度；Sstd为标准辐照度；k为功率温度系数；Tc(t)为t时刻的光伏发电板表面温度；Tref为参考温度；
步骤2.3.2利用式(14)建立风电的数学模型：
Figure FDA0003464565400000041
式(14)中，PWT(t)为风电在t时刻的出力；PN为风机的额定功率；v(t)为t时刻的风速；vin(t)为t时刻风机的切入风速；vN为额定功率时的风速；vout风机的切出风速；
步骤2.3.3通过式(15)得到燃气锅炉的出力模型：
QGB(t)＝ηGBV(t)LV (15)
式(15)中，ηGB为燃气锅炉的效率；V(t)为t时刻消耗的天然气的体积；LV为天然气的低热值；QGB(t)为t时刻燃气锅炉的输出功率；
步骤2.4利用式(16)建立储氢罐的数学模型：
Figure FDA0003464565400000042
式(16)中，PHT(t)为t时刻储氢罐的压强；EHT(t)为t时刻的储氢量；R为阿伏伽德罗常数；
Figure FDA0003464565400000043
为标准开氏温度；VHT为储氢罐的体积；nEL(t)为t时刻电解槽的产氢量；nFC(t)为t时刻燃料电池的耗氢量；EHT(t-1)为t-1时刻的储氢量；Δt为时间间隔；PN为额定压强；
步骤2.5利用式(17)建立氢能综合利用站的碳排放数学模型：
Figure FDA0003464565400000044
式(17)中，βe、βh分别表示单位电量、天然气量对应的碳排放；Pe,graid(t)表示t时刻向电网购买的电量；V(t)表示t时刻消耗的天然气体积；LV表示天然气的低热值；
步骤2.6在氢能综合利用站中，当电解槽工作时，风光发电和市电向电解槽供电，从而利用式(18)得到电母线t时刻的功率平衡约束：
PPV(t)+PWT(t)+δPgrid(t)-Pload(t)＝PEL(t) (18)
式(18)中，PPV(t)为t时刻的光伏出力；PWT(t)为t时刻的风机出力；Pgrid(t)为t时刻与电网交互功率；Pload(t)为t时刻的负荷功率；PEL(t)为t时刻电解槽的输入功率；δ表示0-1变量；
利用式(19)得到热母线的平衡约束：
QGB(t)＝Qload(t) (19)
式(19)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；
步骤2.7在氢能综合利用站中，当电解槽不工作，由燃料电池、风光发电和市电向电负荷供能；由燃气锅炉和燃料电池向热负荷供能，从而利用式(20)得到电母线功率平衡约束：
PPV(t)+PWT(t)+δPgrid(t)+PFC(t)＝Pload(t) (20)
式(20)中，PFC(t)为t时刻燃料电池输出的电功率；
利用式(23)得到热母线平衡约束：
δQGB(t)+QFC(t)＝Qload(t) (21)
式(23)中，QFC(t)为燃料电池输出的热功率；
步骤2.8，利用式(22)建立以碳排量
Figure FDA0003464565400000051
最小为约束的目标函数：
Figure FDA0003464565400000052
3.根据权利要求2所述的一种工业园区氢能综合利用站的能量管理方法，其特征在于，当风光处于出力较大时段时，按如下步骤制定氢能综合利用站工作状态的能量管理方案：
步骤3.1a，利用式(23)判断氢能综合利用站的工作状态，若满足式(23)，则执行步骤3.2a-步骤3.5a，否则，执行步骤4.1；
PPV(t)+PWT(t)≥Pload(t) (23)
步骤3.2a，利用式(24)判断储氢罐是否满足约束，若满足约束，则执行步骤3.3a，否则，电解槽不工作：
SOHT(t)≤SOHTmax (24)
式(24)中，SOHT(t)为t时刻的储氢罐状态；SOHTmax表示储氢罐最大储氢状态量；
步骤3.3a，利用式(25)得到t时刻电解槽的输入功率PEL(t)：
Figure FDA0003464565400000053
式(25)中，CEL为电解槽的容量；
Figure FDA0003464565400000061
为储氢罐的最大储氢容量；EHT(t)为t时刻储氢罐的储氢量；
步骤3.4a，利用式(26)得到满足热负荷的t时刻燃气锅炉出力QGB(t)：
QGB(t)＝Qload(t) (26)
式(26)中，Qload(t)为t时刻热负荷的功率；
步骤3.5a，利用式(27)得到工业园区与外电网购电的交互功率Pgrid(t)：
Pgrid(t)+PPV(t)+PWT(t)＝Pload(t)+PEL(t) (27)。
4.根据权利要求2中所述的一种工业园区氢能综合利用站的能量管理方法，其特征在于，当风光出力处于较小时段时，按如下步骤制定氢能综合利用站工作状态的能量管理方案：
步骤3.1b，利用式(28)判断储氢罐是否满足约束，若满足约束，则执行步骤3.2b，否则，燃料电池不工作；
SOHT(t)≤SOHTmax (28)
式(29)中，SOHT(t)为t时刻的储氢罐状态；SOHTmax表示储氢罐最大储氢量；
步骤3.2b，利用式(29)得到t时刻燃料电池输出的电功率PFC(t)：
Figure FDA0003464565400000062
式(30)中，CFC为燃料电池的容量；
Figure FDA0003464565400000063
为储氢罐的最大储氢容量；EHT(t)为t时刻储氢罐的储氢量；ηFC为燃料电池的工作效率；
步骤3.3b，利用式(30)得到t时刻燃料电池的输出热功率QFC(t)：
Figure FDA0003464565400000064
式(31)中，PFC(t)为t时刻燃料电池输出的电功率；
步骤3.4b，利用式(31)得到t时刻燃气锅炉输出热功率QGB(t)：
Figure FDA0003464565400000065
式(32)中，Qload(t)为t时刻热负荷的功率；QFC(t)为t时刻燃料电池输出的热功率；
步骤3.5b，利用式(32)得到t时刻工业园区与外电网的交互功率Pgrid(t)：
Pgrid(t)＝max{Pload(t)-PPV(t)-PWT(t)-PFC(t),0} (32)
步骤3.6b，利用式(33)得到工业园区所需要的热能：
δQGB(t)+QFC(t)＝Qload(t) (33)
式(33)中，QGB(t)为t时刻燃气锅炉的出力；Qload(t)为t时刻热负荷的功率；QFC(t)为t时刻燃料电池输出的热功率。
CN202210027215.7A 2022-01-11 2022-01-11 工业园区氢能综合利用站的能量管理方法 Active CN114156948B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210027215.7A CN114156948B (zh) 2022-01-11 2022-01-11 工业园区氢能综合利用站的能量管理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210027215.7A CN114156948B (zh) 2022-01-11 2022-01-11 工业园区氢能综合利用站的能量管理方法

Publications (2)

Publication Number Publication Date
CN114156948A true CN114156948A (zh) 2022-03-08
CN114156948B CN114156948B (zh) 2023-07-28

Family

ID=80449719

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210027215.7A Active CN114156948B (zh) 2022-01-11 2022-01-11 工业园区氢能综合利用站的能量管理方法

Country Status (1)

Country Link
CN (1) CN114156948B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019134532A1 (zh) * 2018-01-02 2019-07-11 佛山科学技术学院 一种含风光可再生能源的热电联供系统运行方法
CN111242806A (zh) * 2020-02-19 2020-06-05 武汉理工大学 一种考虑不确定性的电-热-氢多能量系统的规划方法
CN111738503A (zh) * 2020-06-15 2020-10-02 国网安徽省电力有限公司经济技术研究院 以氢能为核心的综合能源微网日前运行调度方法及系统
CN112464477A (zh) * 2020-11-27 2021-03-09 国网山东省电力公司青岛供电公司 计及需求响应的多能耦合综合能源运行仿真方法
CN112736939A (zh) * 2020-12-11 2021-04-30 国网江苏省电力有限公司南通供电分公司 掺氢天然气综合能源系统制氢储氢装置优化容量配置方法
CN113872240A (zh) * 2021-09-30 2021-12-31 国网上海市电力公司 一种风光氢综合能源系统容量配置方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019134532A1 (zh) * 2018-01-02 2019-07-11 佛山科学技术学院 一种含风光可再生能源的热电联供系统运行方法
CN111242806A (zh) * 2020-02-19 2020-06-05 武汉理工大学 一种考虑不确定性的电-热-氢多能量系统的规划方法
CN111738503A (zh) * 2020-06-15 2020-10-02 国网安徽省电力有限公司经济技术研究院 以氢能为核心的综合能源微网日前运行调度方法及系统
CN112464477A (zh) * 2020-11-27 2021-03-09 国网山东省电力公司青岛供电公司 计及需求响应的多能耦合综合能源运行仿真方法
CN112736939A (zh) * 2020-12-11 2021-04-30 国网江苏省电力有限公司南通供电分公司 掺氢天然气综合能源系统制氢储氢装置优化容量配置方法
CN113872240A (zh) * 2021-09-30 2021-12-31 国网上海市电力公司 一种风光氢综合能源系统容量配置方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
曾鸣 等: "综合能源系统建模及效益评价体系综述与展望", vol. 42, no. 06, pages 1697 - 1708 *

Also Published As

Publication number Publication date
CN114156948B (zh) 2023-07-28

Similar Documents

Publication Publication Date Title
Xing et al. Modeling and operation of the power-to-gas system for renewables integration: a review
Lajnef et al. Modeling, control, and simulation of a solar hydrogen/fuel cell hybrid energy system for grid‐connected applications
Feroldi et al. Energy management of a hybrid system based on wind–solar power sources and bioethanol
CN106654319B (zh) 一种基于变异粒子群和差分进化混合算法的pemfc系统温度建模方法
Wang et al. Energy and exergy analysis of a proton exchange membrane water electrolysis system without additional internal cooling
Li et al. Modeling and simulation of hydrogen energy storage system for power-to-gas and gas-to-power systems
CN112271752A (zh) 一种分布式能源系统
Luta et al. Energy management system for a remote renewable fuel cell system
CN117863980A (zh) 一种燃料电池/锂电池混合动力系统瞬时功率匹配及综合热管理方法
CN203907723U (zh) 一种燃料电池建筑发电供暖系统
CN206686115U (zh) 一种基于燃料电池的多能互补热电联产系统
CN105811443A (zh) 基于甲醇水重整制氢发电系统的削峰填谷供电系统及方法
CN116191485A (zh) 一种基于状态机的综合能源系统控制方法
CN114156948B (zh) 工业园区氢能综合利用站的能量管理方法
Ali et al. Intelligent hybrid energy system and grid integration using microcontrollers
Lin et al. Operational characteristics of an integrated island energy system based on multi-energy complementarity
Awad et al. Hybrid PV/FC system design and simulation
Yang et al. Theoretical and Experimental Verification of Wind Powered Hydrogen Storage Energy System
Pierre et al. Implementation of an Advanced PEM Hydrogen Storage System Based Cogeneration Using Photovoltaic System in a Building
CN112907147B (zh) 一种基于最佳工况点的分布式能源系统运行方法及系统
CN103939969A (zh) 一种燃料电池建筑发电供暖系统
Takubo et al. A Hydrogen Management Method in Residential Distributed Generation Systems with Photovoltaic Cells and Hydrogen‐Storage Type Fuel Cells
CN213461148U (zh) 一种分布式能源系统
Haddad et al. An investigation on coupling fuel cell and photovoltaic systems for power generation
CN112751341B (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