CN114734881A - 一种用于氢电混合动力系统的容量优化配置方法 - Google Patents

一种用于氢电混合动力系统的容量优化配置方法 Download PDF

Info

Publication number
CN114734881A
CN114734881A CN202210492230.9A CN202210492230A CN114734881A CN 114734881 A CN114734881 A CN 114734881A CN 202210492230 A CN202210492230 A CN 202210492230A CN 114734881 A CN114734881 A CN 114734881A
Authority
CN
China
Prior art keywords
power
fuel cell
particle
hydrogen
power battery
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
CN202210492230.9A
Other languages
English (en)
Other versions
CN114734881B (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.)
Suzhou Huachang Energy Technology Co ltd
Original Assignee
Suzhou Huachang Energy Technology Co ltd
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 Suzhou Huachang Energy Technology Co ltd filed Critical Suzhou Huachang Energy Technology Co ltd
Priority to CN202210492230.9A priority Critical patent/CN114734881B/zh
Publication of CN114734881A publication Critical patent/CN114734881A/zh
Application granted granted Critical
Publication of CN114734881B publication Critical patent/CN114734881B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L58/00Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
    • B60L58/40Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for controlling a combination of batteries and fuel cells
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L50/00Electric propulsion with power supplied within the vehicle
    • B60L50/50Electric propulsion with power supplied within the vehicle using propulsion power supplied by batteries or fuel cells
    • B60L50/75Electric propulsion with power supplied within the vehicle using propulsion power supplied by batteries or fuel cells using propulsion power supplied by both fuel cells and batteries
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L58/00Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
    • B60L58/10Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries
    • B60L58/12Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries responding to state of charge [SoC]
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60LPROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L58/00Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
    • B60L58/30Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling fuel cells
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/30Hydrogen technology
    • Y02E60/50Fuel cells

Abstract

本发明公开了一种用于氢电混合动力系统的容量优化配置方法,包括根据氢电混合动力系统的约束条件,确定燃料电池额定功率P0的可行域和动力电池额定容量Q0的可行域;建立目标函数
Figure DDA0003631537680000011
确定P0和Q0在各自可行域内使目标函数最小的P0和Q0的值,并根据其分配氢电混合动力系统的燃料电池额定功率和动力电池额定容量,其方法包括设置粒子群组数为和迭代次数为M,计算每组粒子群的目标函数值并比较得到其中的最小值,从而更新局部最小值和全局最小值,并得到当前全局最小值所对应的粒子位置
Figure DDA0003631537680000012
判断当前迭代次数j是否达到M,若是,则令
Figure DDA0003631537680000013
Figure DDA0003631537680000014
本发明能够大幅降低系统运行过程中的损耗,提高系统的经济性。

Description

一种用于氢电混合动力系统的容量优化配置方法
技术领域
本发明涉及氢电混合动力技术领域,特别涉及一种用于氢电混合动力系统的容量优化配置方法。
背景技术
化石能源的日渐匮乏以及尾气排放问题使得新型清洁能源逐步进入人们的视野,其中,质子交换膜燃料电池(PEMFC)具有零排放、功率密度高、启动速度快、工作温度低以及能量转化效率高等诸多优势,在交通领域各类动力系统中得到广泛关注和应用。
目前,PEMFC应用于交通领域时,由于运行工况的随机性,其负载功率中含有大量高频部分,而PEMFC动态响应较慢,且频繁的变载将导致PEMFC寿命的快速衰减,因此通常配合动力电池而共同组成氢电混合动力系统,从而PEMFC和动力电池联合供电以满足动力系统的负载需求。在氢电混合动力系统中,燃料电池和动力电池的容量配比对系统动力性能以及运行经济性有着显著的影响,不合理的容量配置不仅会增加系统的运行氢耗,还会加速燃料电池以及动力电池的寿命衰减,使得系统运行成本大幅增加,同时,氢电混合动力系统的功率分配策略以容量配置为基础,使用不合理的容量配置将导致功率分配策略优化困难。
因此,寻求一种氢电混合动力系统的容量优化配置方法对提高系统经济性具有重要意义。
发明内容
本发明的目的是提供一种降低系统运行损耗的氢电混合动力系统的容量优化配置方法。
为达到上述目的,本发明采用的技术方案如下:
一种用于氢电混合动力系统的容量优化配置方法,其包括:
根据所述氢电混合动力系统的约束条件,确定燃料电池额定功率P0的可行域[P0_min,P0_max]和动力电池额定容量Q0的可行域[Q0_min,Q0_max];
基于系统氢耗
Figure BDA0003631537660000012
燃料电池寿命损耗Lfc_loss及动力电池损耗Lbat_loss建立目标函数
Figure BDA0003631537660000011
c1为氢气单价,c2为燃料电池单价,c3为动力电池单价;
确定P0和Q0在各自可行域内使所述目标函数最小的P0和Q0的值,并根据其分配所述氢电混合动力系统的燃料电池额定功率和动力电池额定容量,其方法包括:
设置粒子群组数为N和迭代次数为M,并初始化局部最小值
Figure BDA0003631537660000021
为正无穷、全局最小值
Figure BDA0003631537660000022
为正无穷以及当前迭代次数j为0,同时初始化粒子初始位置
Figure BDA0003631537660000023
和粒子速度
Figure BDA0003631537660000024
计算每组粒子群的目标函数值并比较得到其中的最小值
Figure BDA0003631537660000025
从而根据所述
Figure BDA0003631537660000026
更新
Figure BDA0003631537660000027
Figure BDA0003631537660000028
并得到当前全局最小值所对应的粒子位置
Figure BDA0003631537660000029
判断当前迭代次数j是否达到M,若是,则令
Figure BDA00036315376600000210
Figure BDA00036315376600000211
否则,令j=j+1,更新粒子位置和速度,并重新计算每组粒子群的目标函数值,直至当前迭代次数j达到M。
进一步地,所述氢电混合动力系统的约束条件包括系统平均功率需求约束条件、系统装配空间限制条件、动力电池SOC波动限制条件及系统瞬时需求功率约束条件;
所述系统平均功率需求约束条件为:
Figure BDA00036315376600000212
式中,P0为燃料电池额定功率,
Figure BDA00036315376600000213
为系统平均功率;
所述系统装配空间限制条件为:
Figure BDA00036315376600000214
Figure BDA00036315376600000215
式中,P0为燃料电池额定功率,ρV_fc为燃料电池比功率,Vconf_fc为系统预留燃料电池装配空间,Q0为动力电池额定容量,ρV_bat为动力电池比能量,Vconf_bat为系统预留动力电池装配空间;
所述动力电池SOC波动限制条件为:
Figure BDA00036315376600000216
ΔEmax=Emax-Emin
Figure BDA00036315376600000217
式中,Q0为动力电池额定容量,ηdischar为动力电池平均放电效率,SOCH为动力电池允许的SOC最高值,SOCL为动力电池允许的SOC最低值,ΔEmax为波动能量的最大变化值,Emax为测试工况下的波动能量E的最大值,Emin为测试工况下的波动能量E的最小值,E为测试工况下的波动能量,k为历史数据总时长,Preq为系统需求功率,
Figure BDA00036315376600000318
为系统平均功率;
所述系统瞬时需求功率约束条件为:
Figure BDA0003631537660000031
式中,P0为燃料电池额定功率,Q0为动力电池额定容量,ηdischar为动力电池平均放电效率,Ic_max为动力电池最大放电倍率,Preq_max为需求功率最大值。
进一步地,所述容量优化配置方法还包括:根据下式初始化粒子初始位置
Figure BDA0003631537660000032
Figure BDA0003631537660000033
Figure BDA0003631537660000034
式中,P0_min和P0_max来自于燃料电池额定功率的可行域[P0_min,P0_max],Q0_min和Q0_max来自于动力电池额定容量的可行域[Q0_min,Q0_max],rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群;
根据下式初始化粒子速度
Figure BDA0003631537660000035
Figure BDA0003631537660000036
Figure BDA0003631537660000037
式中,vmin为粒子最小速度,vmax为粒子最大速度,rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群。
进一步地,所述容量优化配置方法还包括:根据
Figure BDA0003631537660000038
的值配置所述燃料电池额定功率P0,且根据
Figure BDA0003631537660000039
的值配置所述动力电池额定容量Q0,并基于k个小时的历史需求功率数据运行氢电混合动力系统模型以得到k个小时的系统氢耗
Figure BDA00036315376600000310
燃料电池寿命损耗Lfc_loss、动力电池寿命损耗Lbat_loss,从而计算每组粒子群的目标函数值。
进一步地,所述更新
Figure BDA00036315376600000311
Figure BDA00036315376600000312
的方法包括:
Figure BDA00036315376600000313
则令
Figure BDA00036315376600000314
并得到当前局部最小值所对应的粒子位置
Figure BDA00036315376600000315
Figure BDA00036315376600000316
则令
Figure BDA00036315376600000317
并得到当前全局最小值所对应的粒子位置混合动力系统的燃料电池额定功率和动力电池额定容量,其方法包括:
设置粒子群组数为N和迭代次数为M,并初始化局部最小值
Figure BDA0003631537660000041
为正无穷、全局最小值
Figure BDA0003631537660000042
为正无穷以及当前迭代次数j为0,同时初始化粒子初始位置
Figure BDA0003631537660000043
和粒子速度
Figure BDA0003631537660000044
计算每组粒子群的目标函数值并比较得到其中的最小值
Figure BDA0003631537660000045
从而根据所述
Figure BDA0003631537660000046
更新
Figure BDA0003631537660000047
Figure BDA0003631537660000048
并得到当前全局最小值所对应的粒子位置
Figure BDA0003631537660000049
判断当前迭代次数j是否达到M,若是,则令
Figure BDA00036315376600000410
Figure BDA00036315376600000411
否则,令j=j+1,更新粒子位置和速度,并重新计算每组粒子群的目标函数值,直至当前迭代次数j达到M。
进一步地,所述氢电混合动力系统的约束条件包括系统平均功率需求约束条件、系统装配空间限制条件、动力电池SOC波动限制条件及系统瞬时需求功率约束条件;
所述系统平均功率需求约束条件为:
Figure BDA00036315376600000412
式中,P0为燃料电池额定功率,
Figure BDA00036315376600000413
为系统平均功率;
所述系统装配空间限制条件为:
Figure BDA00036315376600000414
Figure BDA00036315376600000415
式中,P0为燃料电池额定功率,ρV_fc为燃料电池比功率,Vconf_fc为系统预留燃料电池装配空间,Q0为动力电池额定容量,ρV_bat为动力电池比能量,Vconf_bat为系统预留动力电池装配空间;
所述动力电池SOC波动限制条件为:
Figure BDA00036315376600000416
ΔEmax=Emax-Emin
Figure BDA00036315376600000417
式中,Q0为动力电池额定容量,ηdischar为动力电池平均放电效率,SOCH为动力电池允许的SOC最高值,SOCL为动力电池允许的SOC最低值,ΔEmax为波动能量的最大变化值,Emax为测试工况下的波动能量E的最大值,Emin为测试工况下的波动能量E的最小值,E为测试工况下的波动能量,k为历史数据总时长,Preq为系统需求功率,
Figure BDA0003631537660000051
为系统平均功率;
所述系统瞬时需求功率约束条件为:
Figure BDA0003631537660000052
式中,P0为燃料电池额定功率,Q0为动力电池额定容量,ηdischar为动力电池平均放电效率,Ic_max为动力电池最大放电倍率,Preq_max为需求功率最大值。
进一步地,所述容量优化配置方法还包括:根据下式初始化粒子初始位置
Figure BDA0003631537660000053
Figure BDA0003631537660000054
Figure BDA0003631537660000055
式中,P0_min和P0_max来自于燃料电池额定功率的可行域[P0_min,P0_max],Q0_min和Q0_max来自于动力电池额定容量的可行域[Q0_min,Q0_max],rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群;
根据下式初始化粒子速度
Figure BDA0003631537660000056
Figure BDA0003631537660000057
Figure BDA0003631537660000058
式中,vmin为粒子最小速度,vmax为粒子最大速度,rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群。
进一步地,所述容量优化配置方法还包括:根据
Figure BDA0003631537660000059
的值配置所述燃料电池额定功率P0,且根据
Figure BDA00036315376600000510
的值配置所述动力电池额定容量Q0,并基于k个小时的历史需求功率数据运行氢电混合动力系统模型以得到k个小时的系统氢耗
Figure BDA00036315376600000511
燃料电池寿命损耗Lfc_loss、动力电池寿命损耗Lbat_loss,从而计算每组粒子群的目标函数值。
进一步地,所述更新
Figure BDA00036315376600000512
Figure BDA00036315376600000513
的方法包括:
Figure BDA00036315376600000514
则令
Figure BDA00036315376600000515
并得到当前局部最小值所对应的粒子位置
Figure BDA00036315376600000516
Figure BDA00036315376600000517
则令
Figure BDA00036315376600000518
并得到当前全局最小值所对应的粒子位置
Figure BDA0003631537660000061
进一步地,所述容量优化配置方法还包括:若当前迭代次数没有达到M,则令j=j+1,并根据下式更新粒子位置和速度:
Figure BDA0003631537660000062
Figure BDA0003631537660000063
式中,
Figure BDA0003631537660000064
为第j次迭代下第i组粒子群的粒子位置,
Figure BDA0003631537660000065
为第j-1次迭代下第i组粒子群的粒子位置,
Figure BDA0003631537660000066
为第j-1次迭代下第i组粒子群的粒子速度,
Figure BDA0003631537660000067
为第j次迭代下第i组粒子群的粒子速度,w为权重系数,c1、c2为加速度常数,
Figure BDA0003631537660000068
为第j-1次迭代下当前局部最小值所对应的粒子位置,r1、r2为随机函数,xglobal_min为当前全局最小值所对应的粒子位置。
进一步地,所述权重系数w和所述加速度常数c1、c2的取值被配置为令所述目标函数在所述迭代次数内收敛。
进一步地,所述系统氢耗
Figure BDA0003631537660000069
根据下式确定:
Figure BDA00036315376600000610
式中,k为系统运行时长,Pfc为燃料电池输出功率,
Figure BDA00036315376600000611
为燃料电池节电压平均值。
进一步地,所述燃料电池寿命损耗Lfc_loss根据下式确定:
Figure BDA00036315376600000612
式中,l为工况修正系数,k为系统运行时长,n1为k小时内燃料电池启停次数,V0为燃料电池寿命到达极限时相对额定功率下的平均节电压,a为燃料电池V0的衰减比,V1为燃料电池启停时的节电压衰退率,V2为燃料电池变载时的节电压衰退率,U1为燃料电池高负载运行下的节电压衰退率,U2为燃料电池低负载运行下的节电压衰退率,t1为k小时内燃料电池高负载运行时间,t2为k小时内燃料电池低负载运行时间,ΔPfc为燃料电池功率变化率。
进一步地,所述动力电池损耗Lbat_loss根据下式确定:
Figure BDA00036315376600000613
式中,α、β为前指因子,SOC为动力电池的荷电状态,R为气体常数,T为运行温度,Ic为充放电倍率,Ah为动力电池安时通量,Q0为动力电池寿命到达极限时的相对额定容量,b为动力电池Q0的衰减比。
本发明具有的优点:在系统的功率约束和能量约束条件下,通过二维粒子群优化方法得到燃料电池额定功率和动力电池额定容量的最优配置方案,能够大幅降低系统运行过程中的损耗,提高系统的经济性,并且,基于系统氢耗、燃料电池损耗、动力电池损耗以及各自市场价格构建优化目标函数,将系统主要损耗统一为经济损耗,从而实现氢电混合动力系统的多目标优化。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的氢电混合动力系统的示意框图;
图2是本发明实施例提供的用于氢电混合动力系统的容量优化配置方法的逻辑示意图;
图3是本发明实施例提供的氢电混合动力系统的历史需求功率Preq曲线图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,更清楚地了解本发明的目的、技术方案及其优点,以下结合具体实施例并参照附图对本发明实施例中的技术方案进行清楚、完整地描述。需要说明的是,附图中未绘示或描述的实现方式,为所属技术领域中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应的值。显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。除此,本发明的说明书和权利要求书中的术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、装置、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
在本发明的一个实施例中,提供了一种用于氢电混合动力系统的容量优化配置方法。
如图1所示,氢电混合动力系统包括燃料电池、动力电池、DC/DC变换器、DC/AC变换器、电动机等,其中,燃料电池和动力电池协调为系统供电;DC/DC变换器为单向DC/DC转换器,其用于燃料电池至直流母线的电压变换;DC/AC变换器用于连接直流母线与电动机;电动机为三相交流电机,用于为系统提供动力。
本容量优化配置方法包括以下步骤:
S1、根据氢电混合动力系统的约束条件,确定燃料电池额定功率P0的可行域[P0_min,P0_max]和动力电池额定容量Q0的可行域[Q0_min,Q0_max]。
具体地,氢电混合动力系统的约束条件包括以下四个方面:
(1)系统平均功率需求约束条件,即系统所配燃料电池额定功率应大于系统平均功率:
Figure BDA0003631537660000081
式中,P0为燃料电池额定功率;
Figure BDA0003631537660000082
为系统平均功率。
(2)系统装配空间限制条件:
Figure BDA0003631537660000083
Figure BDA0003631537660000084
式中,P0为燃料电池额定功率;ρV_fc为燃料电池比功率,单位为kW/L;Vconf_fc为系统预留燃料电池装配空间,单位为L;Q0为动力电池额定容量;ρV_bat为动力电池比能量,单位为kW/L;Vconf_bat为系统预留动力电池装配空间,单位为L。
(3)动力电池SOC波动限制条件:
Figure BDA0003631537660000085
ΔEmax=Emax-Emin
Figure BDA0003631537660000086
式中,Q0为动力电池额定容量;ηdischar为动力电池平均放电效率;SOCH为动力电池允许的SOC最高值;SOCL为动力电池允许的SOC最低值;ΔEmax为波动能量的最大变化值,Emax为测试工况下的波动能量E的最大值,Emin为测试工况下的波动能量E的最小值,E为测试工况下的波动能量,k为历史数据总时长,单位为h;Preq为系统需求功率;
Figure BDA0003631537660000091
为系统平均功率。
(4)系统瞬时需求功率约束条件:
Figure BDA0003631537660000092
式中,P0为燃料电池额定功率;Q0为动力电池额定容量;ηdischar为动力电池平均放电效率;Ic_max为动力电池最大放电倍率,单位为C;Preq_max为需求功率最大值,单位为kW。
根据上述四个约束条件可得出燃料电池额定功率的可行域边界P0_min和P0_max以及动力电池额定容量的可行域边界Q0_min和Q0_max的取值如下:
Figure BDA0003631537660000093
在本实施例中,结合上式并根据各项参数计算得到:
Figure BDA0003631537660000094
S2、基于系统氢耗
Figure BDA0003631537660000095
燃料电池寿命损耗Lfc_loss及动力电池损耗Lbat_loss建立优化目标函数:
Figure BDA0003631537660000096
式中,P0为燃料电池额定功率;Q0为动力电池额定容量;c1为氢气单价,单位为元/kg;c2为燃料电池单价,单位为元/kW;c3为动力电池单价,单位为元/kWh。
在本实施例中,c1取60元/kg,c2取11000元/kW,c3取7500元/kWh。
S3、通过粒子群寻优算法确定P0在[P0_min,P0_max]范围内和Q0在[Q0_min,Q0_max]范围内使目标函数最小的P0和Q0的值,如图2所示,S3主要包括以下几步:
S31、初始化粒子种群大小N=50;迭代次数M=20,初始化迭代次数j=0,初始化局部最小值
Figure BDA0003631537660000097
和全局最小值
Figure BDA0003631537660000098
均为正无穷。在本发明的一个实施例中,N≥20,M≥10,具体数值根据实际情况选择,不以此限定本发明的保护范围。
并且,随机分布粒子初始位置
Figure BDA0003631537660000101
至可行域范围内:
Figure BDA0003631537660000102
Figure BDA0003631537660000103
式中,P0_min和P0_max来自于燃料电池额定功率的可行域[P0_min,P0_max],Q0_min和Q0_max来自于动力电池额定容量的可行域[Q0_min,Q0_max],rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群。
以及,初始化粒子初始位置
Figure BDA0003631537660000104
至设定速度范围内:
Figure BDA0003631537660000105
Figure BDA0003631537660000106
式中,vmin为粒子最小速度,vmax为粒子最大速度,rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群。在本实施例中,vmin取值为-2,vmax取值为+2。
S32、按照
Figure BDA0003631537660000107
的值配置燃料电池额定功率P0,且按照
Figure BDA0003631537660000108
的值配置动力电池额定容量Q0,并选择k个小时(在本实施例中,k取12)的历史需求功率数据(如图3所示),通过运行氢电混合动力系统模型得到k个小时的系统氢耗
Figure BDA0003631537660000109
燃料电池寿命损耗Lfc_loss、动力电池寿命损耗Lbat_loss,从而计算每组粒子群的目标函数值
Figure BDA00036315376600001010
然后从这些目标函数值中选择最小的值作为下一步更新的基础参数,具体过程见S33。
系统氢耗
Figure BDA00036315376600001011
根据下式确定:
Figure BDA00036315376600001012
式中,k为系统运行时长,单位为h;Pfc为燃料电池输出功率,Pfc由系统功率分配策略给出,单位为kW;
Figure BDA00036315376600001013
为燃料电池节电压平均值,
Figure BDA00036315376600001014
由燃料电池模型给出,单位为V。
燃料电池寿命损耗Lfc_loss根据下式确定:
Figure BDA00036315376600001015
式中,l为工况修正系数;k为系统运行时长;n1为k小时内燃料电池启停次数;V0为燃料电池寿命到达极限时相对额定功率下的平均节电压;a为燃料电池V0的衰减比;V1为燃料电池启停时的节电压衰退率;V2为燃料电池变载时的节电压衰退率;U1为燃料电池高负载运行下的节电压衰退率;U2为燃料电池低负载运行下的节电压衰退率;t1为k小时内燃料电池高负载运行时间;t2为k小时内燃料电池低负载运行时间;ΔPfc为燃料电池功率变化率。
动力电池损耗Lbat_loss根据下式确定:
Figure BDA0003631537660000111
式中,α、β为前指因子,α、β与动力电池类型相关;,SOC为动力电池的荷电状态,R为气体常数,T为运行温度,Ic为充放电倍率,Ah为动力电池安时通量,Q0为动力电池寿命到达极限时的相对额定容量,b为动力电池Q0的衰减比。
S33、更新
Figure BDA0003631537660000112
Figure BDA0003631537660000113
则令
Figure BDA0003631537660000114
并得到当前局部最小值所对应的粒子位置
Figure BDA0003631537660000115
由于初始化迭代次数j=0,初始化局部最小值
Figure BDA0003631537660000116
为正无穷,此时,若S32计算得到的
Figure BDA0003631537660000117
小于局部最小值
Figure BDA0003631537660000118
的初始值即正无穷,则将S32计算得到的
Figure BDA0003631537660000119
赋给局部最小值,从而当前
Figure BDA00036315376600001110
的值为S32计算得到的
Figure BDA00036315376600001111
并保存当前局部最小值
Figure BDA00036315376600001112
的值所对应的粒子位置,即
Figure BDA00036315376600001113
从而为S34中更新全局最小值做准备。
S34、更新
Figure BDA00036315376600001114
Figure BDA00036315376600001115
则令
Figure BDA00036315376600001116
并得到当前全局最小值所对应的粒子位置
Figure BDA00036315376600001117
由于初始化全局最小值
Figure BDA00036315376600001118
为正无穷,此时迭代次数j=0,若S33更新后的
Figure BDA00036315376600001119
小于全局最小值
Figure BDA00036315376600001120
的初始值即正无穷,则将S33更新后的
Figure BDA00036315376600001121
赋给全局最小值,从而当前
Figure BDA00036315376600001122
的值为S33更新后的
Figure BDA00036315376600001123
也就是S32计算得到的
Figure BDA00036315376600001124
并保存当前全局最小值
Figure BDA00036315376600001125
所对应的粒子位置,即
Figure BDA00036315376600001126
S35、判断当前迭代次数j是否达到M即20:
若j达到M,则将S34中的
Figure BDA00036315376600001127
作为本次目标函数的最优解,令
Figure BDA00036315376600001128
Figure BDA00036315376600001129
结束迭代,根据P0和Q0来对系统进行优化管理。
若j没有达到M,则令j=j+1,按照下式更新粒子位置和速度,并返回S32以重新计算每组粒子群的目标函数值,直至当前迭代次数j达到M:
Figure BDA00036315376600001130
Figure BDA00036315376600001131
式中,
Figure BDA00036315376600001132
为第j次迭代下第i组粒子群的粒子位置;
Figure BDA00036315376600001133
为第j-1次迭代下第i组粒子群的粒子位置;
Figure BDA0003631537660000121
为第j-1次迭代下第i组粒子群的粒子速度;
Figure BDA0003631537660000122
为第j次迭代下第i组粒子群的粒子速度;w为权重系数,w用于调整搜索范围;c1、c2为加速度常数即粒子的最大加速度;
Figure BDA0003631537660000123
为第j-1次迭代下当前局部最小值所对应的粒子位置;r1、r2为随机函数,r1=r2=rand(0,1);xglobal_min为当前全局最小值所对应的粒子位置。需要注意的是,权重系数w和加速度常数c1、c2的取值需要令目标函数在迭代次数内实现收敛。
S4、根据S3确定的本次目标函数的最优解来分配燃料电池额定功率和动力电池额定容量,完成系统容量的分配。在本实施例中,寻优结果为P0=38kW,Q0=47kWh。
本发明通过建立关于系统氢耗、燃料电池寿命衰减以及动力电池衰减的目标函数,并以该目标函数为优化目标,结合粒子群优化方法寻求在系统约束下的最优容量配置方案,根据本容量优化配置方法可大幅降低系统的运行成本。
以上所述仅为本发明的优选实施例,并非因此限制其专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (10)

1.一种用于氢电混合动力系统的容量优化配置方法,其特征在于,包括:
根据所述氢电混合动力系统的约束条件,确定燃料电池额定功率P0的可行域[P0_min,P0_max]和动力电池额定容量Q0的可行域[Q0_min,Q0_max];
基于系统氢耗
Figure FDA0003631537650000011
燃料电池寿命损耗Lfc_loss及动力电池损耗Lbat_loss建立目标函数
Figure FDA0003631537650000012
c1为氢气单价,c2为燃料电池单价,c3为动力电池单价;
确定P0和Q0在各自可行域内使所述目标函数最小的P0和Q0的值,并根据其分配所述氢电混合动力系统的燃料电池额定功率和动力电池额定容量,其方法包括:
设置粒子群组数为N和迭代次数为M,并初始化局部最小值
Figure FDA0003631537650000013
为正无穷、全局最小值
Figure FDA0003631537650000014
为正无穷以及当前迭代次数j为0,同时初始化粒子初始位置
Figure FDA0003631537650000015
和粒子速度
Figure FDA0003631537650000016
计算每组粒子群的目标函数值并比较得到其中的最小值
Figure FDA0003631537650000017
从而根据所述
Figure FDA0003631537650000018
更新
Figure FDA0003631537650000019
Figure FDA00036315376500000110
并得到当前全局最小值所对应的粒子位置
Figure FDA00036315376500000111
判断当前迭代次数j是否达到M,若是,则令
Figure FDA00036315376500000112
Figure FDA00036315376500000113
否则,令j=j+1,更新粒子位置和速度,并重新计算每组粒子群的目标函数值,直至当前迭代次数j达到M。
2.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,所述氢电混合动力系统的约束条件包括系统平均功率需求约束条件、系统装配空间限制条件、动力电池SOC波动限制条件及系统瞬时需求功率约束条件;
所述系统平均功率需求约束条件为:
Figure FDA00036315376500000114
式中,P0为燃料电池额定功率,
Figure FDA00036315376500000115
为系统平均功率;
所述系统装配空间限制条件为:
Figure FDA00036315376500000116
Figure FDA00036315376500000117
式中,P0为燃料电池额定功率,ρV_fc为燃料电池比功率,Vconf_fc为系统预留燃料电池装配空间,Q0为动力电池额定容量,ρV_bat为动力电池比能量,Vconf_bat为系统预留动力电池装配空间;
所述动力电池SOC波动限制条件为:
Figure FDA0003631537650000021
ΔEmax=Emax-Emin
Figure FDA0003631537650000022
式中,Q0为动力电池额定容量,ηdischar为动力电池平均放电效率,SOCH为动力电池允许的SOC最高值,SOCL为动力电池允许的SOC最低值,ΔEmax为波动能量的最大变化值,Emax为测试工况下的波动能量E的最大值,Emin为测试工况下的波动能量E的最小值,E为测试工况下的波动能量,k为历史数据总时长,Preq为系统需求功率,
Figure FDA0003631537650000023
为系统平均功率;
所述系统瞬时需求功率约束条件为:
Figure FDA0003631537650000024
式中,P0为燃料电池额定功率,Q0为动力电池额定容量,ηdischar为动力电池平均放电效率,Ic_max为动力电池最大放电倍率,Preq_max为需求功率最大值。
3.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,还包括:根据下式初始化粒子初始位置
Figure FDA0003631537650000025
Figure FDA0003631537650000026
Figure FDA0003631537650000027
式中,P0_min和P0_max来自于燃料电池额定功率的可行域[P0_min,P0_max],Q0_min和Q0_max来自于动力电池额定容量的可行域[Q0_min,Q0_max],rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群;
根据下式初始化粒子速度
Figure FDA0003631537650000028
Figure FDA0003631537650000029
Figure FDA00036315376500000210
式中,vmin为粒子最小速度,vmax为粒子最大速度,rand(0,1)为随机函数,i为N组粒子群中的第i组粒子群。
4.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,还包括:
根据
Figure FDA0003631537650000031
的值配置所述燃料电池额定功率P0,且根据
Figure FDA0003631537650000032
的值配置所述动力电池额定容量Q0,并基于k个小时的历史需求功率数据运行氢电混合动力系统模型以得到k个小时的系统氢耗
Figure FDA0003631537650000033
燃料电池寿命损耗Lfc_loss、动力电池寿命损耗Lbat_loss,从而计算每组粒子群的目标函数值。
5.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,所述更新
Figure FDA0003631537650000034
Figure FDA0003631537650000035
的方法包括:
Figure FDA0003631537650000036
则令
Figure FDA0003631537650000037
并得到当前局部最小值所对应的粒子位置
Figure FDA0003631537650000038
Figure FDA0003631537650000039
则令
Figure FDA00036315376500000310
并得到当前全局最小值所对应的粒子位置
Figure FDA00036315376500000311
6.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,还包括:若当前迭代次数没有达到M,则令j=j+1,并根据下式更新粒子位置和速度:
Figure FDA00036315376500000312
Figure FDA00036315376500000313
式中,
Figure FDA00036315376500000314
为第j次迭代下第i组粒子群的粒子位置,
Figure FDA00036315376500000315
为第j-1次迭代下第i组粒子群的粒子位置,
Figure FDA00036315376500000316
为第j-1次迭代下第i组粒子群的粒子速度,
Figure FDA00036315376500000317
为第j次迭代下第i组粒子群的粒子速度,w为权重系数,c1、c2为加速度常数,
Figure FDA00036315376500000318
为第j-1次迭代下当前局部最小值所对应的粒子位置,r1、r2为随机函数,xglobal_min为当前全局最小值所对应的粒子位置。
7.根据权利要求6所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,所述权重系数w和所述加速度常数c1、c2的取值被配置为令所述目标函数在所述迭代次数内收敛。
8.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,所述系统氢耗
Figure FDA0003631537650000041
根据下式确定:
Figure FDA0003631537650000042
式中,k为系统运行时长,Pfc为燃料电池输出功率,
Figure FDA0003631537650000043
为燃料电池节电压平均值。
9.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,所述燃料电池寿命损耗Lfc_loss根据下式确定:
Figure FDA0003631537650000044
式中,l为工况修正系数,k为系统运行时长,n1为k小时内燃料电池启停次数,V0为燃料电池寿命到达极限时相对额定功率下的平均节电压,a为燃料电池V0的衰减比,V1为燃料电池启停时的节电压衰退率,V2为燃料电池变载时的节电压衰退率,U1为燃料电池高负载运行下的节电压衰退率,U2为燃料电池低负载运行下的节电压衰退率,t1为k小时内燃料电池高负载运行时间,t2为k小时内燃料电池低负载运行时间,ΔPfc为燃料电池功率变化率。
10.根据权利要求1所述的用于氢电混合动力系统的容量优化配置方法,其特征在于,所述动力电池损耗Lbat_loss根据下式确定:
Figure FDA0003631537650000045
式中,α、β为前指因子,SOC为动力电池的荷电状态,R为气体常数,T为运行温度,Ic为充放电倍率,Ah为动力电池安时通量,Q0为动力电池寿命到达极限时的相对额定容量,b为动力电池Q0的衰减比。
CN202210492230.9A 2022-05-07 2022-05-07 一种用于氢电混合动力系统的容量优化配置方法 Active CN114734881B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210492230.9A CN114734881B (zh) 2022-05-07 2022-05-07 一种用于氢电混合动力系统的容量优化配置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210492230.9A CN114734881B (zh) 2022-05-07 2022-05-07 一种用于氢电混合动力系统的容量优化配置方法

Publications (2)

Publication Number Publication Date
CN114734881A true CN114734881A (zh) 2022-07-12
CN114734881B CN114734881B (zh) 2023-12-12

Family

ID=82284689

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210492230.9A Active CN114734881B (zh) 2022-05-07 2022-05-07 一种用于氢电混合动力系统的容量优化配置方法

Country Status (1)

Country Link
CN (1) CN114734881B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003163013A (ja) * 2001-11-27 2003-06-06 Nippon Telegr & Teleph Corp <Ntt> 燃料電池システム制御法
CN106887841A (zh) * 2017-03-23 2017-06-23 东北大学 一种含电动汽车微电网容量配置的多种群遗传粒子群优化方法
CN110190628A (zh) * 2019-06-12 2019-08-30 西南交通大学 一种基于综合能源的复线全并联at牵引供电系统及方法
CN110834625A (zh) * 2019-11-11 2020-02-25 常熟理工学院 一种自适应异步粒子群的双电耦合燃料电池汽车能效优化方法
CN111162295A (zh) * 2019-12-30 2020-05-15 电子科技大学 一种考虑退化的燃料电池混合系统的能量管理方法
US20210224351A1 (en) * 2020-01-21 2021-07-22 Kinaxis Inc. Method and system for optimizing an objective having discrete constraints
CN113644674A (zh) * 2021-09-06 2021-11-12 华北电力大学 基于量子粒子群的氢混合储能容量配置系统、方法及应用

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003163013A (ja) * 2001-11-27 2003-06-06 Nippon Telegr & Teleph Corp <Ntt> 燃料電池システム制御法
CN106887841A (zh) * 2017-03-23 2017-06-23 东北大学 一种含电动汽车微电网容量配置的多种群遗传粒子群优化方法
CN110190628A (zh) * 2019-06-12 2019-08-30 西南交通大学 一种基于综合能源的复线全并联at牵引供电系统及方法
CN110834625A (zh) * 2019-11-11 2020-02-25 常熟理工学院 一种自适应异步粒子群的双电耦合燃料电池汽车能效优化方法
CN111162295A (zh) * 2019-12-30 2020-05-15 电子科技大学 一种考虑退化的燃料电池混合系统的能量管理方法
US20210224351A1 (en) * 2020-01-21 2021-07-22 Kinaxis Inc. Method and system for optimizing an objective having discrete constraints
CN113644674A (zh) * 2021-09-06 2021-11-12 华北电力大学 基于量子粒子群的氢混合储能容量配置系统、方法及应用

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姚海涛: "含氢储的混合储能在风光互补发电系统中容量优化研究", 电气开关 *

Also Published As

Publication number Publication date
CN114734881B (zh) 2023-12-12

Similar Documents

Publication Publication Date Title
CN109606137B (zh) 融合成本寿命因素的多源电驱动系统经济性优化方法
Etxeberria et al. Hybrid energy storage systems for renewable energy sources integration in microgrids: A review
CN108987770A (zh) 一种多堆燃料电池发电系统的协调优化控制方法
CN104716644A (zh) 一种可再生能源冷热电微网系统及控制方法
WO2014153946A1 (zh) 独立微网系统的优化方法
CN110758122B (zh) 一种电-电混合电源系统的燃料电池效率优化方法
CN112060982B (zh) 一种动态平衡的燃料电池无人机能量管理方法
CN109713666B (zh) 一种基于K-means聚类的电力市场下分布式储能经济性调控方法
CN106099965A (zh) 交流微电网并网状态下复杂混合储能系统的协调控制方法
CN113555590B (zh) 多燃料电池模块发电系统及其控制方法
CN106911136B (zh) 基于温度和功率控制平抑分布式能源功率波动的方法及系统
CN110752629B (zh) 一种交直流混合家庭微电网能量优化管理方法
CN114889498B (zh) 一种氢电混合动力系统的功率优化分配方法
CN115907240B (zh) 计及互补互济运行特性的电网多类型调峰资源规划方法
CN111525624A (zh) 基于蓄电池储能系统的家庭分布式能源调度方法
CN115549159A (zh) 一种考虑调峰成本的大规模电动汽车分群调度策略
CN113270884B (zh) 一种基于v2g微网系统的电力能源管理系统和方法
CN108564230A (zh) 一种户用分布式能源管理方法及系统
CN111231713A (zh) 一种电动汽车充放电系统及控制方法
CN204407890U (zh) 一种可再生能源冷热电微网系统
CN110470993B (zh) 一种启停电池soc算法
CN114734881A (zh) 一种用于氢电混合动力系统的容量优化配置方法
CN115107538B (zh) 汽车的能量管理方法及装置
CN111509716A (zh) 电网柔性负荷控制方法、装置、计算机设备和存储介质
CN116361920A (zh) 基于nsgaⅱ算法的氢燃料电池混动汽车能量处理方法

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