一种氢电混合动力系统的功率分配方法

1.本发明属于燃料电池汽车及船舶技术领域，具体涉及一种氢电混合动力系统的功率分配方法。


背景技术：

2.化石能源的日渐匮乏以及尾气排放问题使得新型清洁能源逐步进入人们的视野。质子交换膜燃料电池(pemfc)具有零排放、功率密度高、启动速度快、工作温度低以及能量转化效率高等诸多优势，在交通领域得到广泛关注和应用。
3.质子交换膜燃料电池是一种以氢气为燃料、氧气为氧化剂的电化学发电装置。氢气在电堆阳极发生氧化反应失去电子转化为质子，通过质子交换膜到达阴极并与阴极产生的氧离子生成产物水。质子交换膜燃料电池具有诸多优势，但由于质子交换燃料电池动态响应慢，而实际交通工况的需求功率变化快，单一的燃料电池难以匹配动力系统所需功率，因此常采用燃料电池 辅助能源的方式为动力系统提供动力。这种氢电混合动力系统既克服了燃料电池动态响应慢，纹波耐受能力弱的缺点，又克服了纯电动汽车续航里程短的问题。但多级能量源的设计模式也引入了新的问题，即系统的功率分配问题，不合理的功率分配策略不仅会导致系统氢耗的增加，还会大幅加快燃料电池的性能衰减和辅助电源的寿命衰减，因此寻求一种优化的功率分配策略对提高系统整体性能、降低系统成本具有重要意义。


技术实现要素：

4.针对上述现有技术中存在的问题，本发明提出了一种氢电混合动力系统的功率分配方法，根据电动机需求功率分配燃料电池承担低频功率部分，同时采用soc闭环控制策略给燃料电池输出功率叠加一个微调量，使辅助电源控制在正常范围内，进而保证动力系统的可靠运行。
5.本发明具体技术方案如下：
6.一种氢电混合动力系统的功率分配方法，其特征在于，包括以下步骤：
7.步骤1：设置一个容量为2n的实时双向队列空间，电动机实时需求功率p
req
依次填充实时双向队列空间，并采用滑动窗口的形式不断更新被完全填充的实时双向队列空间，获得更新后实时双向队列空间中填充的数据；其中，n≥4；
8.步骤2：对更新后实时双向队列空间中填充的数据进行n层小波分解和小波重构，得到电动机实时需求功率p
req
的低频分量pn，作为燃料电池的基量；其中，小波分解的层数n的获取步骤为：
9.步骤2.1：选取一个被电动机历史需求功率p
req
′
完全填充的容量为2m的历史双向队列空间，初始化小波分解的层数m＝1；其中，m≥10；
10.步骤2.2：对历史双向队列空间填充的数据进行m层小波分解和小波重构，得到重构信号ym；
11.步骤2.3：对重构信号ym进行快速傅里叶变换，得到[0,fs/2]频率区间内各频点的幅值a
k,m
和频率f
k,m
，计算各频点的幅值a
k,m
与频率f
k,m
的乘积a
k,m
×fk,m
；其中，fs为电动机历史需求功率p
req
′
的数据采样频率；
[0012]
步骤2.4：若在[0,fs/2]频率区间内均满足a
k,m
×fk,m
≤ac×
fc，则选择当前层数m的值为小波分解的层数n；否则，令m＝m 1，转回步骤2.2；其中，fc和ac分别为截止频率和截止频率所对应的幅值，获取步骤为：
[0013]
步骤2.4.1：获取燃料电池的极限变载速率a；
[0014]
步骤2.4.2：对历史双向队列空间填充的数据进行傅里叶级数展开并按谐波频率升序排列，得到：
[0015]
x(t)＝b0 b
1 sin(2πf1t ψ1) b
2 sin(2πf2t ψ2) ... b
k sin(2πfkt ψk) ...
[0016]
其中，x(t)为历史双向队列空间填充的数据；b0为常数项；b1为一次谐波系数，f1为一次谐波频率，ψ1为一次谐波的相位偏移量；b2为二次谐波系数，f2为二次谐波频率，ψ2为二次谐波的相位偏移量；bk为k次谐波系数，fk为k次谐波频率，ψk为k次谐波的相位偏移量；
[0017]
步骤2.4.3：初始化k＝1；
[0018]
步骤2.4.4：根据sin函数特性求取k次谐波分量的斜率最大值r＝2πb
kfk
cos0＝2πb
kfk
；
[0019]
步骤2.4.5：若r《a，则令k＝k 1，转回步骤2.4.4；否则，令ac＝bk，fc＝fk；
[0020]
步骤3：实时采集辅助电源的soc(电量)值soc
battery
，设定soc期望值soc
ref
，根据soc
battery
和soc
ref
计算pi(比例积分)控制器的输出量u：
[0021][0022]
其中，k
p
为比例增益；ti为积分常数；τ为当前时间；
[0023]
步骤4：对输出量u进行截止频率为fc的低通滤波，低通滤波后输出p
lowpass
，作为燃料电池的微调量，使辅助电源的soc值控制在正常范围内；
[0024]
步骤5：对燃料电池的基量pn和微调量p
lowpass
进行求和，并将求和结果限幅在[p
min
,p
max
]区间范围内，得到燃料电池的输出功率p
fc
；其中，p
min
为燃料电池的最低输出功率；p
max
为燃料电池的最高输出功率；
[0025]
步骤6：计算辅助电源的输出功率p
bat
＝p
req-p
fc
。
[0026]
进一步地，步骤2中小波分解采用haar(哈尔小波)小波为母小波，haar小波的函数ψ(x)为：
[0027][0028]
进一步地，步骤2.4.1中获取燃料电池的极限变载速率a的方法为：固定燃料电池的变载时间为1s，以0.1kw为间隔逐步增大变载步长，当燃料电池的实际输出功率无法实现变载步长的要求时，记录此时的变载步长a
′
kw，则燃料电池的极限变载速率a＝a
′
kw/s。
[0029]
进一步地，步骤3中k
p
的取值方法为：
[0030]
若soc
battery
≤soc0，则k
p
＝k
pmax
；
[0031]
若soc0＜soc
battery
＜soc1，则
[0032]
若soc1≤soc
battery
≤soc2，则k
p
＝0；
[0033]
若soc2＜soc
battery
＜soc3，则
[0034]
若soc
battery
≥soc3，则k
p
＝k
pmax
；
[0035]
其中，soc0为辅助电源出现过放现象时的soc临界值；soc3为辅助电源出现过充现象时的soc临界值；soc1为辅助电源正常工作时的soc下限；soc2为辅助电源正常工作时的soc上限；k
pmax
为比例增益极限值。
[0036]
进一步地，k
pmax
的获取方式为：以soc期望值soc
ref
为目标，在soc
battery
≤soc0和soc
battery
≥soc3两种情况下调节pi控制器参数，当两种情况下超调量和响应时间均满足氢电混合动力系统要求时，此时pi控制器调控的比例增益的值为k
pmax
，积分常数的值为ti。
[0037]
进一步地，步骤5中燃料电池的输出功率p
fc
的公式为：
[0038][0039]
本发明的有益效果为：
[0040]
1、本发明提出了一种氢电混合动力系统的功率分配方法，根据电动机需求功率分配燃料电池承担低频功率部分，辅助电源承担高频功率部分，采用soc闭环控制策略给燃料电池输出功率叠加一个微调量，使辅助电源控制在正常范围内，进而保证动力系统的可靠运行；
[0041]
2、本发明采用小波分解和重构算法提取电动机需求功率的低频成分，并运用傅里叶级数分解的方法将小波分解层数和低通滤波的截止频率与燃料电池的极限变载速率相关联，实现燃料电池的功率匹配；
[0042]
3、优选地，本发明通过设计比例增益跟随辅助电源soc值变化的pi控制器，在燃料电池输出基量上叠加一个相对微弱的调控量，实现了辅助电源的电量平衡；相比传统pi控制器，本发明中设计的pi控制器在soc值偏离正常区间时可更快地进行调控，而在正常区间内时不调控，减小了控制难度。
附图说明
[0043]
图1为本发明实施例1中氢电混合动力系统的拓扑图；
[0044]
图2为本发明实施例1提出的氢电混合动力系统的功率分配方法的流程图；
[0045]
图3为本发明实施例1中电动机实时需求功率p
req
的曲线图；
[0046]
图4为本发明实施例1中燃料电池的基量pn的曲线图；
[0047]
图5为本发明实施例1中燃料电池的输出功率p
fc
的曲线图；
[0048]
图6为本发明实施例1中燃料电池的输出功率p
fc
变化率的曲线图；
[0049]
图7为本发明实施例1中辅助电源的soc曲线图。
具体实施方式
[0050]
为使本发明的目的、技术方案和优点更加清晰，结合以下具体实施例，并参照附图，对本发明做进一步的说明。
[0051]
下述非限制性实施例可以使本领域的普通技术人员更全面的理解本发明，但不以任何方式限制本发明。
[0052]
实施例1
[0053]
本实施例基于如图1所示的氢电混合动力系统，提出了一种氢电混合动力系统的功率分配方法。
[0054]
所述氢电混合动力系统包括能量管理系统、燃料电池、辅助电源、dc/dc模块、dc/ac模块和电动机；其中，能量管理系统根据功率分配方法控制燃料电池和辅助电源进行功率输出；燃料电池的最大功率为50kw，最大变载能力为3kw/s；辅助电源的容量为138ah，电压范围为570～760v；dc/dc模块为单向dc/dc转换器，用于燃料电池至直流母线的电压变换；dc/ac模块为双向dc/ac转换器，用于连接直流母线与电动机；电动机包括三相交流电机和电机控制器。
[0055]
所述功率分配方法如图2所示，具体包括以下步骤：
[0056]
步骤1：设置一个容量为29＝512的实时双向队列空间，电动机实时需求功率p
req
依次填充实时双向队列空间，电动机实时需求功率p
req
的曲线如图3所示，并采用滑动窗口的形式不断更新被完全填充的实时双向队列空间，即实时双向队列空间被完全填充时，删除实时双向队列空间的第1个位置填充的数据，之后各位置填充的数据往前滑动一个位置，最新的电动机实时需求功率p
req
填充至最后一个位置，获得更新后实时双向队列空间，填充的数据为[a1,a2,...,a
512
]；其中，a1、a2、
…
、a
512
分别为双向队列空间中第1个、第2个、
…
、第512个位置处填充的数据。
[0057]
步骤2：对进行n层小波分解，得到第n层分解的近似信号xn，小波重构后得到电动机实时需求功率p
req
的低频分量pn，作为燃料电池的基量，如图4所示；其中，小波分解采用haar小波为母小波，haar小波的函数ψ(x)为：
[0058][0059]
小波分解的层数n的获取步骤具体为：
[0060]
步骤2.1：选取一个被电动机历史需求功率p
req
′
完全填充的容量为2
10
＝1024的历史双向队列空间，初始化小波分解的层数m＝1；
[0061]
步骤2.2：利用haar滤波器对历史双向队列空间填充的数据进行m层小波分解，并对第m层分解的近似信号进行小波重构，得到重构信号ym；
[0062]
步骤2.3：对重构信号ym进行快速傅里叶变换，得到[0,fs/2]频率区间内各频点的幅值a
k,m
和频率f
k,m
，计算各频点的幅值a
k,m
与频率f
k,m
的乘积a
k,m
×fk,m
；其中，fs＝1hz为电动机历史需求功率p
req
′
的数据采样频率；
[0063]
步骤2.4：若在[0,fs/2]频率区间内均满足a
k,m
×fk,m
≤ac×
fc，则选择当前层数m的值为小波分解的层数n，在本实施例中n＝6；否则，令m＝m 1，转回步骤2.2；其中，fc和ac分别
为截止频率和截止频率所对应的幅值，获取步骤为：
[0064]
步骤2.4.1：获取燃料电池的极限变载速率a＝3kw/s；
[0065]
步骤2.4.2：对历史双向队列空间填充的数据进行傅里叶级数展开并按谐波频率升序排列，得到：
[0066]
x(t)＝b0 b
1 sin(2πf1t ψ1) b
2 sin(2πf2t ψ2) ... b
k sin(2πfkt ψk) ...
[0067]
其中，x(t)为历史双向队列空间填充的数据；b0为常数项；b1为一次谐波系数，f1为一次谐波频率，ψ1为一次谐波的相位偏移量；b2为二次谐波系数，f2为二次谐波频率，ψ2为二次谐波的相位偏移量；bk为k次谐波系数，fk为k次谐波频率，ψk为k次谐波的相位偏移量；
[0068]
步骤2.4.3：初始化k＝1；
[0069]
步骤2.4.4：根据sin函数特性求取k次谐波分量的斜率最大值r＝2πb
kfk
cos0＝2πb
kfk
；
[0070]
步骤2.4.5：若r《3，则令k＝k 1，转回步骤2.4.4；否则，令ac＝bk，fc＝fk；在本实施例中ac＝28kw，fc＝0.016hz；
[0071]
步骤3：实时采集辅助电源的soc值soc
battery
，设定soc期望值soc
ref
＝75％，根据soc
battery
和soc
ref
计算pi控制器的输出量u：
[0072][0073]
其中，k
p
为比例增益；ti为积分常数；τ为氢电混合动力系统的当前时间；
[0074]kp
的取值方法为：
[0075]
若soc
battery
≤soc0，则k
p
＝k
pmax
；
[0076]
若soc0＜soc
battery
＜soc1，则
[0077]
若soc1≤soc
battery
≤soc2，则k
p
＝0；
[0078]
若soc2＜soc
battery
＜soc3，则
[0079]
若soc
battery
≥soc3，则k
p
＝k
pmax
；
[0080]
其中，soc0＝30％为辅助电源出现过放现象时的soc临界值；soc3＝90％为辅助电源出现过充现象时的soc临界值；soc1＝50％为辅助电源正常工作时的soc下限；soc2＝70％为辅助电源正常工作时的soc上限；k
pmax
为比例增益极限值，获取方式为：以soc期望值soc
ref
＝75％为目标，在soc
battery
＝30％和soc
battery
＝90％两种情况下调节pi控制器参数，当两种情况下超调量低于20％，响应时间低于2h时，此时pi控制器调控的比例增益的值为k
pmax
＝0.6，积分常数的值为ti＝0.003。
[0081]
步骤4：利用低通滤波器对输出量u进行截止频率为fc＝0.016hz的低通滤波，低通滤波后输出p
lowpass
，作为燃料电池的微调量，使辅助电源的soc值控制在正常范围内；
[0082]
步骤5：对燃料电池的基量pn和微调量p
lowpass
进行求和，并利用限幅器将求和结果限幅在[p
min
,p
max
]区间范围内，得到燃料电池的输出功率p
fc
，具体为：
[0083][0084]
其中，p
min
＝0为燃料电池的最低输出功率；p
max
＝40kw为燃料电池的最高输出功率；
[0085]
燃料电池的输出功率p
fc
的曲线如图5所示，任意时刻输出功率p
fc
的变化率如图6所示，燃料电池输出功率p
fc
的变化率均小于燃料电池极限变载能力。
[0086]
步骤6：计算辅助电源的输出功率p
bat
＝p
req-p
fc
，整个循环工况下的辅助电源soc曲线如图7所示。
[0087]
在本实施例中，电动机实时需求功率p
req
、辅助电源的soc值soc
battery
、燃料电池的基量pn、输出量u、微调量p
lowpass
、燃料电池输出功率p
fc
、辅助电源的输出功率p
bat
均为与时间t相关的参数。
[0088]
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。