质子交换膜燃料电池退化预测方法

汪建锋1 王荣杰1,2 林安辉1 王亦春1 张 博1

(1. 集美大学轮机工程学院 厦门 361021 2. 电工材料电气绝缘国家重点实验室(西安交通大学) 西安 710049)

摘要 耐久度是制约质子交换膜燃料电池大规模应用的主要障碍之一,性能退化预测技术可以有效提高质子交换膜燃料电池的耐久度。该文提出一种结合小波阈值去噪方法的正则化堆叠长短期记忆网络的性能退化预测方法。通过小波阈值去噪法,获得消除噪声和尖峰后的平滑数据。针对退化数据不确定性和高度非线性导致的特征难以提取问题,引入了正则化堆叠长短期记忆网络模型,该模型通过引入参数优化算法有效地避免了过拟合风险,提高了预测精度和可靠性。为验证该方法的有效性,采用两种不同工况下的质子交换膜燃料电池老化数据进行验证。结果表明,所提方法在稳态工况下的最大误差为0.016 3 V,误差区间在0.5%以内;动态工况下的最大误差为0.006 4 V,误差区间在0.2%以内。

关键词:质子交换膜燃料电池 性能退化预测 小波阈值去噪 长短期记忆网络

0 引言

质子交换膜燃料电池(Proton Exchange Membrane Fuel Cell, PEMFC)由于其零污染、能量转换效率高、低噪声和快速启动的优点[1-2],被认为是最有前途的新能源之一。目前PEMFC已经被应用于混合动力船舶、电动汽车以及微电网等方向[3-7]。然而,质子交换膜燃料电池退化带来的耐久性和高维护成本严重阻碍了PEMFC在各领域的大规模应用[8]。因此,M. Jouin等在2013年首次提出了燃料电池寿命预测和健康管理技术(Prognostics and Health Management, PHM)的概念[9],其中退化预测为PEMFC预知和健康管理的关键。

PEMFC的退化问题涉及多物理和多尺度的问题,是一个非常复杂的非线性过程,整个电堆当中涉及材料和应用中的各种因素都有可能影响PEMFC的降解。因此,准确预测PEMFC性能退化是一个很大的挑战。PEMFC系统的退化受启停条件、启动方式、可变负载、高载、过载、低负载和空载等多种状态的混合影响。每种操作状态都会因为其独特的动态特性对PEMFC的耐久度产生不利影响,且复杂的操作环境会加速PEMFC的老化。PEMFC在典型运行条件下的降解机理如图1所示。

width=228,height=144.75

图1 PEMFC在典型运行条件下的降解机理

Fig.1 Degradation mechanism of PEMFC under typical operating conditions

PEMFC退化指标主要有:电压、电流、内阻以及电化学活性面积。质子交换膜燃料电池的退化将会导致电化学活性面积下降以及电压、电流、内阻等指标的变化[10-11]。相比之下,电压老化数据能很好地反映燃料电池内部的退化趋势,并且易于获取,因此本文采用电压指标作为性能退化预测的健康指标。现有的主流质子交换膜燃料电池退化预测方法主要分为三大类:基于物理模型的方法[12-13]、基于数据驱动的方法[14-15]和基于混合模型的方法[16-18]

基于物理模型的方法主要采用特定的物理模型或半经验退化模型对燃料电池进行退化预测,目前这类方法以滤波算法和半经验模型为主,典型的代表方法有粒子滤波(Particle Filter, PF)算法[12]、无迹卡尔曼滤波(Unscented Kalman Filter, UKF)算法[19]、扩展卡尔曼滤波(Extended Karman Filter, EKF)算法[20]和半经验模型[21]等,这些方法虽然能对PEMFC的退化过程进行预测,但由于未完全了解PEMFC内部降解机理,使得其非线性处理能力有限且计算效率较低,因此针对PEMFC老化过程的高度非线性特性难以进行有效的预测。

基于数据驱动的预测方法通过学习大量的历史数据来建立老化模型,不需要对PEMFC的材料特性、加载条件、失效机制和老化机理等进行深入的了解[22],因此模型搭建难度较低。基于数据驱动的方法主要是基于人工智能算法[23],文献[24]提出了一种基于遗传优化算法的堆叠状态回升网络(Stacked Echo State Network Based on Genetic Algorithm Optimization, S-ESN-GAO)的预测方法,该方法通过遗传学算法对模型的超参数进行优化,所提出的方法有更强的泛化能力和更高的预测精度。文献[25]提出了一种堆叠的长短期记忆网络(Stack Long Short-Term Memory, S-LSTM)的预测方法,通过堆叠和dropout操作获得了不错的预测效果。文献[23]提出了一种基于网格长短期记忆递归神经网络(Grid Long Short-Term Memory, G-LSTM)的预测方法,在传统递归神经网络的基础上加入了长短期记忆网络单元避免梯度爆炸和梯度消失,然后采用不同工况条件下的燃料电池数据进行验证,得到的结果有较高的准确度和泛化能力。上述方法有效地实现了PEMFC的性能退化预测,但预测精度还有待提高。

基于混合模型的预测方法采用模型驱动和数据驱动的混合方法[14]。文献[21]提出了一种基于自回归滑动平均模型(Autoregressive and Moving Average, ARMA)和延迟神经网络(Time Delay Neural Network, TDNN)的燃料电池预测方法,在该方法中ARMA用于滤波,TDNN用于预测。文献[26]提出了一种结合深度信念网络(Deep Belief Network, DBN)和极限学习机(Extreme Learning Machine, ELM)的质子交换膜燃料电池退化预测方法,在一定程度上提高了预测精度和可靠性。文献[27]提出了一种基于最小绝对值收敛套索算法和回声状态网络(Least Absolute Shrinkage and Selection Operator-Echo State Network, LASSO-ESN)的退化预测方法,采用LASSO评估输入参数对质子交换膜燃料电池性能的贡献,通过对输入参数进行迭代更新的策略实现退化预测,并有较高的准确度。混合模型的方法在特定的条件下性能优越,但其模型的搭建更加复杂,同时训练时间较长[28]。目前,基于混合模型的预测方法在预测精度上比基于物理模型的预测方法更高,但是其模型搭建的复杂度和计算复杂度也随着预测精度的提高而增加,因此并不方便部署和实现在线预测。

对比基于物理模型的方法和基于混合模型的方法,基于数据驱动的方法可以灵活地表示老化电压数据中观察到的退化特征,且无需对燃料电池的先验知识有深入的了解[29];此外,基于数据驱动的方法易于部署,计算复杂度较低,更适合于实际的在线应用。LSTM作为一种循环神经网络,在长时间序列预测问题上取得了很好的成果,其能有效避免梯度爆炸和梯度消失的优点更适合处理数千小时持续运行下PEMFC的性能退化预测问题。本文采用了一种小波阈值去噪(Wavelet Threshold Denoising, WTD)的方法对PEMFC的老化数据进行平滑去噪处理,采用正则化对点长短期记忆网络(Regularization Stack Long Short-Term Memory, RS-LSTM)模型进行预测。采用两种不同工况下的PEMFC老化数据进行验证,将老化数据分为五种不同的样本大小对WTD-RS-LSTM进行验证。此外还设置了不对老化数据进行WTD平滑去噪处理的RS-LSTM预测实验。以上实验结果表明,结合小波阈值去噪方法的正则化堆叠长短期记忆网络(Regularization Stack Long Short-Term Memory combined with Wavelet Threshold Denoising, WTD-RS-LSTM)方法可以有效实现对PEMFC性能退化的高准确度预测。

1 问题描述

1.1 测试平台说明

本文所使用的1 kW燃料电池测试平台如图2所示,测试台由五片单电池组成,活化面积为100 cm2,电流密度为width=21.45,height=12width=48,height=14.55。测试台来自法国FCLAB,稳态工况下(FC1)和动态工况下(FC2)的老化数据来自文献[30]。其原理如图3所示。

width=224.25,height=108

图2 PEMFC系统测试平台[30]

Fig.2 PEMFC test bench[30]

width=215.25,height=83.25

图3 测试台原理

Fig.3 The principal scheme of the test bench

测试台主要由加湿系统、供气子系统、氢气子系统、电子负载、冷却系统和PEMFC电堆组成。为了避免燃料电池堆的“饥饿”现象,测试台采用压力和流量控制器来调节气体的供应。

另外,在这个系统中增加了一个额外的冷却水子系统用来将堆栈的内部温度保持在可接受的范围内。此外,为了有效地防止氢气泄漏和控制电堆阴阳极之间的压力差,系统中还增加了氢气循环泵和背压阀组件。测试台适用于1 kW燃料电池电堆[30],在额定功率下的物理参数控制范围见表1。

表1 测试台参数控制范围

Tab.1 Test bench parameter control range

参数控制范围 冷却温度/℃20~80 冷却流量/(L/min)0~10 气体温度/℃20~80 气体相对湿度(%)空气流量/(L/min)氢气流量/(L/min)0~1000~1000~30 气体压力/Pa燃料电池电流/A0~2×1050~300

1.2 退化特性分析

通过实际的PEMFC老化数据更有助于分析电压的老化特性,从而在下一步建立更准确的退化预测模型。FC1电压退化-恢复现象如图4所示,从实际的电压老化曲线可以观察到PEMFC电压老化趋势具有以下的特征。

width=218.25,height=129.75

图4 FC1电压退化-恢复现象

Fig.4 FC1 voltage degradation-recovery phenomenon

1)高度非线性。在整个质子交换膜燃料电池退化当中,由于其恢复现象,质子交换膜燃料电池电压衰退趋势总是呈现非线性的特征。衰退趋势在一个小周期中呈现单调的退化趋势,但在进入下一个小周期时发生电压恢复现象,使得整个衰退过程呈现高度的非线性趋势。

2)周期性。电压退化和恢复现象总是在每个小周期重复发生。

3)随机性。燃料电池运行时进/出口氢气温度进/排气温度冷却水入/出口温度参数的随机变化将导致电压呈现随机性的变化,这在短周期内更为明显。此外,每次电压恢复的幅度不同,这也证明了退化数据具有随机性。

1.3 预测当中的问题

根据以上对PEMFC退化数据的分析,现有的PEMFC性能退化预测方法存在两大问题:

1)特征提取能力低。PEMFC电压老化数据具有高度的非线性、周期性以及压降和电压恢复的随机性特点,老化数据由于长期的老化实验使得其数据量非常庞大。以上特性使得老化电压数据的特征提取变得更加困难,传统算法难以完全提取其有效特征。

2)缺乏泛化能力。现有的多数算法只能单独预测FC1数据或者FC2数据,即算法只能较为准确地预测到稳态或者暂态工况下的退化趋势,而不能做到二者皆可。且它们的预测性能易受训练数据量的影响,当其训练数据量发生改变时(通常是减小),其预测精度将会发生下降。

本文旨在设计一个具有高精度、高可靠性的PEMFC性能退化预测方法。WTD对原始数据进行预处理可以较大程度地降低后续网络特征提取难度;LSTM网络具有避免梯度爆炸和梯度消失的优点,堆叠的方式可以进一步加强网络的特征提取能力,对参数进行优化和动态调整学习率的方式可以有效防止预测模型发生过拟合。因此本文提出WTD-RS-LSTM的预测方法以实现PEMFC性能退化的高准确度和高可靠性预测。

2 基于WTD-RS-LSTM的融合预测框架

2.1 小波阈值去噪

WTD具有突出时频定位、多尺度和多分辨率的特性[31],该方法能在去除原始数据中存在的噪声和尖峰的同时保持信号的有效特征,因此采用该方法对PEMFC的老化数据进行平滑处理[32],平滑过程如图5所示。

width=228,height=105.75

图5 小波阈值去噪原理

Fig.5 Wavelet threshold denoising schematic

WTD包括三个阶段:小波分解、阈值处理和信号重构,通过对原始信号进行小波分解得到不同的小波系数,然后根据信号对应的属性选择适当的阈值和阈值处理函数,最后通过数据重构获得去噪后的信号。小波分解将母小波函数变换成与原始信号紧密匹配的基函数,在此将本文需要平滑的时间序列定义为width=21.45,height=14.55,其表达式为

width=125.15,height=17.15 (1)

width=164.55,height=33.45 (2)

式中,tab分别为时间序列中的时刻、缩放因子和沿时间轴的位移因子;width=39.45,height=29.15为小波函数。

母小波函数width=30,height=16.3和变换后的小波基函数width=36,height=16.3表示为

width=205.7,height=33.45(3)

原始信号width=21.45,height=14.55变换后得到的width=36,height=17.15可表示为

width=86.55,height=29.15 (4)

式中,width=9.45,height=10.3为信号分解层数。信号分解后的近似部分保存在width=13.7,height=16.3当中,噪声部分保存在width=13.7,height=16.3中,阈值处理过程分为软阈值处理和硬阈值处理,在本文中采用软阈值处理,将小波系数的绝对值低于阈值的置零,高于阈值的作系数收缩处理,表达式为

width=133.7,height=36.85(5)

式中,width=11.15,height=10.3为小波系数;width=10.3,height=11.15为阈值。令width=24,height=14.55为将width=36,height=16.3中各元素绝对值按从小到大排列后取二次方的信号序列,表达式见式(6),则width=10.3,height=11.15可表示为式(7)。

width=171.45,height=25.7 (6)

width=118.3,height=17.15 (7)

WTD的原理为削弱信号中的无用部分,增强信号中的有用部分,通过将信号进行小波分解后利用小波系数的大小区分是否为噪声信号,从而达到正常信号和噪声信号的区分。

2.2 长短期记忆网络

LSTM是一种特殊的循环神经网络(Recurrent Neural Netuork, RNN),解决了RNN在长序列训练过程中的梯度爆炸和梯度消失的问题,能在长时间序列训练中取得更好的效果。LSTM的网络结构如图6所示,该单元由遗忘门、输入门、输出门和细胞状态更新组成,由以下复合函数实现。

width=126,height=17.15 (8)

width=97.7,height=15.45 (9)

width=118.3,height=18 (10)

width=66.85,height=17.15(11)

width=111.45,height=15.45 (12)

width=58.3,height=15.45 (13)

width=221.25,height=83.25

图6 LSTM存储单元结构

Fig.6 LSTM memory cell structure

式中,σ为Sigmoid函数;width=52.3,height=12分别为遗忘门、细胞状态更新、输入门和输出门;width=42.85,height=13.7分别为输入权重、循环权重和偏置函数;C为细胞状态;width=13.7,height=15.45为输出值。tanh函数的作用是缓解梯度消失并加速收敛,其函数表示为

width=85.7,height=25.7 (14)

2.3 正则化堆叠长短期记忆网络

高度非线性和不确定性是PEMFC电堆电压退化数据特征难以提取的重要因素,为克服这些问题,提高PEMFC性能退化预测的准确率,提出一种RS-LSTM的预测模型,模型结构如图7所示。该模型为了解决PEMFC电堆电压退化数据当中存在的特征难以提取的问题,采用堆叠的方法提高模型的学习能力。为防止过拟合,在模型的参数更新环节添加参数优化算法对模型当中的参数进行优化。在该网络结构中,width=6.85,height=11.15时步的权重梯度和偏置向量的梯度表示为

width=108.85,height=30 (15)

width=93.45,height=30 (16)

式中,width=10.3,height=16.3t时间步长的误差项。权重矩阵和偏置向量的最终梯度可以表示为

width=64.3,height=30 (17)

width=49.7,height=30 (18)

式中,j为单元数量。

width=468,height=204

图7 RS-LSTM结构

Fig.7 RS-LSTM structure diagram

网络在训练过程当中将产生大量的参数,堆叠的LSTM使得参数数量成倍地增长,这导致模型具有强学习能力的同时,也容易造成过拟合。因此对模型中的参数进行优化,具体的参数优化方法目标函数可以表示为

width=228,height=29.15 (19)

式中,width=21.45,height=17.15为损失函数;N为样本数量;width=21.45,height=17.15为待学习的神经网络;width=7.7,height=11.15为范数函数;width=10.3,height=11.15为函数width=21.45,height=17.15width=21.45,height=14.55的参数;width=10.3,height=11.15为正则化系数。width=21.45,height=12.85范数在零点不可导,因此近似为

width=74.55,height=29.15 (20)

式中,D为参数数量;width=9.45,height=10.3为一个非常小的常数。参数优化的函数可以表示为

width=226.3,height=45.45 (21)

2.4 基于融合方法的预测实现

为了提高模型的预测准确度,提出了一种基于WTD-RS-LSTM的PEMFC性能退化预测方法,预测框架如图8所示。在提出的框架当中,电堆电压被作为老化指标。WTD方法可以重构出更适合进行预测的数据,重构的输出数据经过堆叠的LSTM进行特征提取和训练,网络训练得到的参数在经过优化算法后,可以防止过拟合。选取方均根误差(Root Mean Square Error, RMSE)和平均绝对百分比误差(Mean Absolute Percentage Error, MAPE)两个不同的性能指标来评估模型的预测性能,计算公式分别表示为

width=129.45,height=32.55 (22)

width=143.15,height=35.15 (23)

式中,width=53.15,height=18分别为实际电压值和预测电压值;T为整个预测时间当中的采样点总数。RMSE和MAPE的值越小,说明模型性能越好。

width=474.75,height=258.75

图8 WTD-RS-LSTM的PEMFC性能退化预测框架

Fig.8 WTD-RS-LSTM PEMFC performance degradation prediction framework

具体的实验步骤总结如下:

(1)对原始老化数据进行预处理,随后FC1和FC2数据被划分为不同长度的训练集和验证集,两种工况下不同长度的验证集用来验证模型的泛化能力和通用性。

(2)将训练集输入模型当中,训练模型参数初始化。

(3)模型通过前向传播和反向传播计算神经单元的参数,参数通过优化算法被优化,参数优化算法的目标函数为式(19)。

(4)Warmup策略动态调整网络的学习率,Adam优化器改进模型参数。

(5)判断模型的性能是否达到预设指标,如果达到则输出训练好的模型参数;否则重复步骤(3)和步骤(4)。

(6)调整模型的超参数,直到达到最优性能后,输出训练好的预测模型。

3 结果与讨论

3.1 数据预处理

本文采用的数据集来自法国FCLAB实验室在IEEE PHM公开挑战赛上提供的数据。FC1老化数据由1 kW燃料电池电堆在70 A的恒定电流下运行1 154 h获取;FC2老化数据由1 kW燃料电池电堆在频率为5 Hz、振荡为7~70 A动态电流下运行 1 120 h获取。系统采集的数据变量见表2,FC1老化数据集包含143 862组老化数据、FC2老化数据集包含127 371组老化数据。老化指标的老化曲线如图9所示,FC1和FC2电堆的长期老化测试中的电压老化数据中都有非常明显的老化趋势,随时间推移逐渐下降。此外,电压老化数据中包含大量的噪声和尖峰,大量的数据点也会造成计算时间的延长。因此,在本研究中,对原始电压数据进行数据预处理,预处理的方法采用WTD方法。将电压老化数据表示为式(1),原始信号与母小波函数匹配生成的小波基函数表示为式(3)。通过多次的实验对比,确定将WTD算法的分解层数设置为6层时能得到最后的平滑效果,此时小波基函数可以表示为

表2 数据采集系统监测变量

Tab.2 The data acquisition system monitors the variables

参数物理意义 单电池 Utot/VI/AJ/(A/cm2)电堆电压电流电流密度 进/出口氢气温度 Tin/out_AIR/℃进/排气温度 Tin/out_WAT/℃冷却水入/出口温度 进/排气压力 Pin/out_AIR/mbar氢气入/出口压力 /(L/min)氢气出/入口流量 Din/out_AIR/(L/min)进/排气流量 D_WAT/(L/min)冷却水流量 (%)氢气湿度

注:1bar=10-5Pa

width=201.75,height=128.25

width=192,height=125.25

图9 FC1和FC2平滑后曲线

Fig.9 FC1 and FC2 smoothed post-curves

width=86.55,height=29.15 (24)

WTD方法中的width=24,height=14.55序列可以表示为

width=184.3,height=35.15(25)

式中,对于FC1和FC2,N分别为143 862和127 371。FC1和FC2平滑后曲线如图9所示。由图9可知,重构后的FC1和FC2电压退化数据,不仅保留了原始数据当中的主要特征,而且有效去除了噪声和尖峰,电压数据被明显地平滑处理了。该项操作降低了不确定性对预测结果的干扰。

3.2 退化预测

本文的仿真实验平台为AMD Ryzen 7 5800 8-Core Processor@3.40 GHz处理器和16 GB内存,运行环境为Python3.9+tensorflow2.8.0。

在本节中,FC1和FC2数据被用于训练WTD-RS-LSTM并验证预测结果。设置了五种不同长度的训练数据集,其长度分别550 h、700 h、800 h、900 h和1 000 h,其余的数据集作为验证集对算法进行验证。WTD-RS-LSTM对PEMFC在两种不同工况下的五种不同训练状态的预测结果如图10和图11所示,评价指标RMSE和MAPE计算结果见表3。预测误差如图12和图13所示。根据图12和图13可以看出,所提方法在稳态工况下的最大误差为0.016 3 V,误差区间在0.5%以内;动态工况下的最大误差为0.0064 V,误差区间在0.2%以内。分析表3中不同情况下的预测性能指标,该算法在稳态情况下和动态情况下均取得较好的预测准确度,同时动态工况下的预测性能要优于稳态情况下的预测性能。对不同的训练长度,动态训练长度大于550 h后,数据集的退化-恢复现象变得更加剧烈。

width=224.25,height=128.25

图10 稳态工况下在不同训练长度下的预测结果

Fig.10 Steady-state predictions with different training lengths

width=222,height=126.75

图11 动态工况下在不同训练长度下的预测结果

Fig.11 Dynamic-state predictions with different training lengths

表3 WTD-RS-LSTM预测性能指标

Tab.3 WTD-RS-LSTM predictive performance indicators

训练长度/h稳态动态 RMSE/VMAPE(%)RMSE/VMAPE(%) 550 h0.002 060.000 490.000 750.000 20 700 h0.002 030.000 480.000 980.000 22 800 h0.001 370.000 390.001 670.000 32 900 h0.000 790.000 230.001 520.000 32 1 000 h0.000 910.000 220.001 250.000 29

width=225.75,height=128.25

图12 稳态情况下的误差曲线

Fig.12 Error curves under steady-state conditions

width=231.75,height=128.25

图13 暂态情况下的误差曲线

Fig.13 Error curves under dynamic-state conditions

3.3 对比试验

为了验证本文所提方法的预测性能,将预测结果与RVM、LSTM、EKF-LSTM、CNN、RNN和RS-LSTM方法的预测结果进行对比,并比较不同训练数据集长度下各个模型预测性能指标。稳态工况不同训练长度下预测的各项性能指标结果见表4;动态工况不同训练长度下预测的各项结果见表5,对比如图14所示。通过对550 h、700 h和800 h三种不同训练数据集长度下的预测结果进行比较,可以看出,随着训练长度的增加,模型的预测性能也得到提高。RVM、LSTM、CNN、RNN和EKF-LSTM在稳态情况下可获得较好的预测性能,但以上这五种方法对动态情况下的PEMFC进行预测时则达不到与稳态工况预测时相同的预测性能;RS-LSTM方法在稳态情况下获得较高准确度的同时,在动态的情况下会得到更好的预测性能;此外,添加了WTD方法的WTD-RS-LSTM方法在稳态和动态的情况下均能获得最好的预测性能,且在更难以预测的动态情况下能获得比稳态情况下更好的预测性能,说明了所提方法具有更好的特征学习能力。预测结果说明,无论是在稳态还是动态工况下,所提出的方法都有更好的预测性能和鲁棒性。

同时,根据图14可以看出,RS-LSTM在对PEMFC进行预测时对比其他网络均具有更好的预测准确度,这点在更加复杂的动态工况下的PEMFC上得到了更明显的体现。在通过添加WTD的方法对原始数据进行预处理后,稳态和动态情况下的预测精度均得到了较大的提升,验证了所提方法的有效性。

表4 稳态工况下不同训练长度下的性能分析

Tab.4 Performance analysis under steady-state working conditions under different training lengths

训练长度/hRVM[33]LSTM[33]EKF-LSTM[33]CNNRNNRS-LSTMWTD-RS-LSTM RMSE/V5500.017 00.016 80.009 20.027 60.018 30.012 850.002 06 7000.018 60.012 10.011 00.001 90.007 10.008 910.002 03 8000.012 80.010 70.010 40.008 80.006 90.006 100.001 37 MAPE(%)5500.004 40.004 10.002 20.006 40.005 20.003 510.000 49 7000.005 20.003 00.002 20.005 50.001 80.002 370.000 48 8000.003 40.002 20.002 20.002 30.001 70.001 430.000 39

表5 动态工况下不同训练长度下的性能分析

Tab.5 Performance analysis under different training lengths under variable working conditions

训练长度/hRVM[33]LSTM[33]EKF-LSTM[33]CNNRNNRS-LSTMWTD-RS-LSTM RMSE/V5500.032 10.025 10.021 60.013 50.006 10.003 020.000 75 7000.033 10.033 50.019 50.014 00.022 00.003 030.000 98 8000.032 10.023 90.026 20.010 90.009 70.004 020.001 67 MAPE(%)5500.008 60.006 60.005 50.003 10.001 50.000 720.000 20 7000.008 00.008 70.004 60.003 40.005 40.000 720.000 22 8000.008 10.007 10.006 10.002 90.002 50.000 910.000 32

width=437.25,height=261.75

图14 稳态和动态情况下的RMSE和MAPE

Fig.14 RMSE and MAPE in steady and dynamic state

4 结论

本文提出了一种新的基于数据驱动的PEMFC性能退化预测方法,该方法通过WTD对原始数据进行平滑去噪,然后使用RS-LSTM预测PEMFC性能退化曲线。该方法在PEMFC性能退化预测当中能准确提取老化电压数据中的高度非线性特征并防止过拟合。利用所提出的方法在不同长度训练集的情况下,对稳态和动态情况下的PEMFC进行性能退化预测,验证了所提方法对比传统的预测方法具有更高的准确度和可靠性。与传统方法相比,所建立的退化预测模型可以在各种训练样本大小取得更高的预测精度和更稳定的性能,且易于部署,适合在线运用。可以确定,所提出的预测模型对于PEMFC性能退化预测具有很高的潜力。

参考文献

[1] 陈锦洲, 林飞, 何洪文, 等. 质子交换膜燃料电池/电解槽系统建模及负荷追踪策略[J]. 电工技术学报, 2020, 35(增刊2): 636-643. Chen Jinzhou, Lin Fei, He Hongwen, et al. Proton exchange membrane fuel cell/electrolyzer hybrid power system modeling and load tracking strategy[J]. Transactions of China Electrotechnical Society, 2020, 35(S2): 636-643.

[2] 李梓丘, 乔颖, 鲁宗相. 海上风电-氢能系统运行模式分析及配置优化[J]. 电力系统自动化, 2022, 46(8): 104-112. Li Ziqiu, Qiao Ying, Lu Zongxiang, et al. Operation mode analysis and configuration optimization of offshore wind-hydrogen system[J]. Automation of Electric Power Systems, 2022, 46(8): 104-112.

[3] Song Yuxi, Zhang Caizhi, Ling Chunyu, et al. Review on current research of materials fabrication and application for bipolar plate in proton exchange membrane fuel cell[J]. International Journal of Hydrogen Energy, 2020, 45(54): 29832-29847.

[4] 张雪霞, 黄平, 蒋宇, 等. 动态机车工况下质子交换膜燃料电池电堆衰退性能分析[J]. 电工技术学报, 2022, 37(18): 4798-4806. Zhang Xuexia, Huang Ping, Jiang Yu, et al. Degradation performance analysis of proton exchange membrane fuel cell stack under dynamic locomotive conditions[J]. Transactions of China Electrotechnical Society, 2022, 37(18): 4798-4806.

[5] 高锋阳, 高翾宇, 张浩然, 等. 全局与瞬时特性兼优的燃料电池有轨电车能量管理策略[J]. 电工技术学报, 2023, 38(21): 5923-5938. Gao Fengyang, Gao Xuanyu, Zhang Haoran, et al. Management strategy for fuel cell trams with both global and transient characteristics[J]. Transactions of China Electrotechnical Society. 2023, 38(21): 5923-5938.

[6] 马小勇, 王议锋, 王萍, 等. 燃料电池用交错并联型Boost变换器参数综合设计方法[J]. 电工技术学报, 2022, 37(2): 397-408. Ma Xiaoyong, Wang Yifeng, Wang Ping, et al. Comprehensive parameter design method of interleavedBoost converter for fuel cell applications[J]. Transactionsof China Electrotechnical Society, 2022, 37(2): 397-408.

[7] 唐钧涛, 戚志东, 裴进, 等. 基于电荷泵的燃料电池有源网络升压变换器[J]. 电工技术学报, 2022, 37(4): 905-917. Tang Juntao, Qi Zhidong, Pei Jin, et al. An active network DC-DC Boost converter with a charge pump employed in fuel cells[J]. Transactions of China Electrotechnical Society, 2022, 37(4): 905-917.

[8] Taghiabadi M M, Zhiani M, Silva V, et al. Effect of MEA activation method on the long-term performance of PEM fuel cell[J]. Applied Energy, 2019, 242: 602-611.

[9] Jouin M, Gouriveau R, Pera M, et al. Prognostics and health management of PEMFC–state of the art and remaining challenges[J]. International Journal of Hydrogen Energy, 2013, 38(35): 15307-15317.

[10] Liu Hao, Chen Jian, Hissel D, et al. Prognostics methods and degradation indexes of proton exchange membrane fuel cells: a review[J]. Renewable and Sustainable Energy Reviews, 2020, 123: 109721.

[11] Pei Pucheng, Chen Dongfang, Wu Zhiyao, et al. Nonlinear methods for evaluating and online predicting the lifetime of fuel cells[J]. Applied Energy, 2019, 254: 113730.

[12] Ou Mingyang, Zhang Ruofan, Shao Zhifang, et al. A novel approach based on semi-empirical model for degradation prediction of fuel cells[J]. Journal Of Power Sources, 2020, 488: 229435-229445.

[13] Zhang Xinfeng, Yang Daijun, Luo Minghui, et al. Load profile based empirical model for the lifetime prediction of an automotive PEM fuel cell[J]. International Journal of Hydrogen Energy, 2017, 42(16): 11868-11878.

[14] Ibrahim M, Steiner Y, Jemei S, et al. Wavelet-based approach for online fuel cell remaining useful lifetime prediction[J]. IEEE Transactions Industrial Electronics, 2016, 63(8): 5057-5068.

[15] Liu Hui, Liu Zhenyu, Jia Weiqiang, et al. Remaining useful life prediction using a novel feature-attention-based end-to-end approach[J]. IEEE Transactions on Industrial Informatics, 2016, 17(2): 1197-1207.

[16] Cheng Yujie, Zerhouni N, et al. A hybrid remaining useful life prognostic method for proton exchange membrane fuel cell[J]. International Journal of Hydrogen Energy, 2018, 43(27): 12314-12327.

[17] Ma Jian, Liu Xue, Zou Xinyu, et al. Degradation prognosis for proton exchange membrane fuel cell based on hybrid transfer learning and intercell differences[J]. ISA Transactions, 2021, 113: 149-165.

[18] Zhou Daming, Gao Fei, Elena B, et al. Degradation prediction of PEM fuel cell using a moving window based hybrid prognostic approach[J]. Energy, 2017, 138: 1175-1186.

[19] Chen K, Laghrouche S. Fuel cell health prognosis using unscented Kalman filter: Postal fuel cell electric vehicles case study[J]. International Journal of Hydrogen Energy, 2019, 44(3): 1930-1939.

[20] Bressel M, Hilairet M, Hissel D, et al. Extended Kalman filter for prognostic of proton exchange membrane fuel cell[J]. Applied Energy, 2016, 164: 220-227.

[21] Zhou Daming, Al-Durra A, Zhang Ke, et al. Online remaining useful lifetime prediction of proton exchange membrane fuel cells using a novel robust methodology[J]. Journal of Power Sources, 2018, 399: 314-328.

[22] Yin Shen, Xie Xiaochen, Lam J, et al. An improved incre-mental learning approach for KPI prognosis of dynamic fuel cell system[J]. IEEE Transactions on Cybernetics, 2016, 46(12): 3135-3144.

[23] Ma Rui, Yang Tao, Elena B, et al. Data-driven proton exchange membrane fuel cell degradation predication through deep learning method[J]. Applied Energy, 2018, 231: 102-115.

[24] Deng Zhihua, Chen Qihong, Zhang Liyan, et al. Degradation prediction of PEMFCs using stacked echo state network based on genetic algorithm optimization[J]. IEEE Transactions on Transportation Electrification, 2022, 8(1): 1454-1466.

[25] Fu-Kwun W, Cheng Xiaobin, Kai-Chun H, et al. Stacked long short-term memory model for proton exchange membrane fuel cell systems degradation[J]. Power Sourses, 2020, 488(22): 7591-7598.

[26] Xie Yucen, Zou Jianxiao, Li Zhongliang, et al. A novel deep belief network and extreme learning machine based performance degradation prediction method for proton exchange membrane fuel cell[J]. IEEE Access, 2020, 17: 6661-6675.

[27] He Kai, Mao L, Yu Jianbo, et al. Long-term performance prediction of PEMFC based on LASSO-ESN[J]. IEEE Transactions on Instrumentation and Measurement, 2021, 70: 1-11.

[28] Liu Jiawei, Li Qi. Remaining useful life prediction of PEMFC based on long short-term memory recurrent neural networks[J]. International Journal of Hydrogen Energy, 2018, 44(11): 5470-5480.

[29] Jouin M, Gouriveau R, Hissel D, et al. Prognostics of PEM fuel cell in a particle filtering framework[J]. International Journal of Hydrogen Energy, 2014, 39: 481-494.

[30] FCLAB Research. IEEE PHM 2014 Data challenge [Z/OL]. 2014.

[31] Wang Fuyu, Cen Jian, Yu Zongwei, et al. Research on a hybrid model for cooling load prediction based on wavelet threshold denoising and deep learning: a study in China[J]. Energy Reports, 2022, 8: 10950-10962.

[32] Peng Shanbi, Chen Ruolei, Yu Bin, et al. Daily natural gas load forecasting based on the combination of long short term memory, local mean decomposition, and wavelet threshold denoising algorithm[J]. Journal of Natural Gas Science and Engineering, 2021, 95: 104175.

[33] Ma Rui, Xie Renyou, Xu Liangcai, et al. A hybrid prognostic method for PEMFC with aging parameter prediction[J]. IEEE Transactions on Transportation Electrification, 2021, 7(4): 2318-2331.

Degradation Prediction Method of Proton Exchange Membrane Fuel Cell

Wang Jianfeng1 Wang Rongjie1,2 Lin Anhui1 Wang Yichun1 Zhang Bo1

(1. School of Marine Engineering Jimei University Xiamen 361021 China 2. State Key Laboratory of Electrical Insulation and Power Equipment Xi’an Jiaotong University Xi’an 710049 China)

Abstract Durability is one of the main obstacles to the large-scale application of proton exchange membrane fuel cell (PEMFC). Performance degradation prediction technology can effectively improve the durability of PEMFC. Through the study of PEMFC aging data, it is found that the actual PEMFC aging data is highly nonlinear, periodic and random, which makes it difficult for the prediction algorithm to extract the features effectively. In addition, in the problem of degradation prediction, the prediction algorithm needs to predict the degradation of PEMFC under different working conditions, which requires the prediction algorithm to have stronger generalization ability. To solve the above problems, a performance degradation prediction method of regularization stack long short-term memory combined with wavelet threshold denoising method (WTD-RS-LSTM) method is proposed. Firstly, the WTD method is used to process the original data, and the smooth data after eliminating noise and spikes is obtained by wavelet decomposition, threshold processing and data reconstruction. Then the RS-LSTM model is introduced to solve the problem of feature extraction caused by uncertainty and high nonlinearity of degraded data. The generalization ability of the model is improved by introducing parameter optimization algorithm. The model is stacked to enhance its learning ability. For increase the reliability of the model, Warmup strategy was used to dynamically adjust the learning rate of the network. Through the above operations, the overfitting phenomenon which may occur in the training of the model is effectively avoided, and the prediction accuracy and reliability of the prediction algorithm are improved. For verify the effectiveness of the proposed method, PEMFC aging data under two different working conditions are used for verification. The datasets under different working conditions are divided into five different lengths of training sets and test sets to train and test the proposed algorithm. The verification results show that under steady-state conditions, the maximum error of the proposed method is 0.016 3 V, and the error interval is within 0.5%. The prediction performance increases with the training length, and the best prediction performance is obtained at the training length of 1 000 h, when the RMSE and MAPE are 0.000 91 and 0.000 22%, respectively. Under dynamic conditions, the maximum error is 0.006 4 V and the error interval is within 0.2%. The best performance was achieved when the training length was 550 h, when the RMSE and MAPE are 0.000 75 and 0.000 20%, respectively. According to the above experimental results and the comparison with the existing traditional algorithms, the following conclusions are drawn: (1) the proposed method can make more accurate PEMFC degradation prediction under different working conditions and different training lengths, and has stronger generalization ability; (2) Comparing the prediction accuracy of the two conditions under different training lengths, it is found that the prediction of PEMFC degradation under dynamic conditions by the proposed method is better than that under steady-state conditions. Therefore, the proposed method has stronger prediction ability under dynamic conditions. (3) The proposed method has a simple structure, easy to deploy and is suitable for online application; (4) The aging of PEMFC under dynamic conditions will produce more randomness, which will have a great impact on the stability of the prediction algorithm.

Keywords: Proton exchange membrane fuel cell (PEMFC), degradation prediction, wavelet threshold denoising, long short-term memory (LSTM)

DOI: 10.19595/j.cnki.1000-6753.tces.230326

中图分类号:TM911; TK91

国家自然科学基金(51879118)、福建省自然科学基金(2020J01688)、电力设备电气绝缘国家重点实验室基金(EIPE23202)和福建省中青年教师教育科研项目(JAT220173)资助。

收稿日期 2023-03-17

改稿日期 2023-03-31

作者简介

汪建锋 男,1999年生,硕士研究生,研究方向为质子交换膜燃料电池寿命预测。E-mail: 1307135971@qq.com

王荣杰 男,1981年生,博士,教授,研究方向为智能信息处理和电力电路故障诊断。E-mail: roger811207@163.com(通信作者)

(编辑 郭丽军)