一种基于HSMM和经验模型的燃料电池故障预测方法与流程

本发明属于燃料电池
技术领域：
，更为具体地讲，涉及一种基于hsmm和经验模型的燃料电池故障预测方法。
背景技术：
：燃料电池是一个通过电化学反应把燃料中化学能量直接转化为电能的装置，它具有环境友好、能量转换率高等优点，然而，可靠性和持久性是阻碍燃料电池发展的一个主要原因。故障预测技术可以在燃料电池故障早期及时的诊断它的健康状态、预测其剩余可用寿命，使得我们可以及时的对它进行维修，以此延长燃料电池寿命。近几十年来，故障预测技术在解决燃料电池的寿命相对较短的问题上表现出了很大的潜力。基于隐式半马尔科夫模型(hsmm)的故障预测方法是一种基于数据的预测方法，不需要对燃料电池的内部结构和物理退化过程进行建模，直接对系统的退化数据进行分析即可预测其剩余寿命。然而当燃料电池的电压退化趋势与hsmm全寿命训练数据的电压退化趋势相似度较低时，基于hsmm的预测结果较差。为了解决这个问题，本发明提出了将hsmm与经验模型相结合的混合故障预测方法，若燃料电池的当前运行电压退化趋势与hsmm训练电压退化趋势相似度较高，则利用hsmm对燃料电池的剩余寿命进行预测；若两者的相似度较低，则利用经验模型估计燃料电池剩余运行寿命。技术实现要素：本发明的目的在于克服现有技术的不足，提供一种基于hsmm和经验模型的燃料电池故障预测方法，根据燃料电池的退化数据，对燃料电池的剩余寿命进行预测。为实现上述发明目的，本发明一种基于hsmm和经验模型的燃料电池故障预测方法，其特征在于，包括以下步骤：(1)、设置初始状态概率矩阵π0、初始状态转移概率矩阵a0、初始观测值概率矩阵b0、初始状态持续时间概率矩阵p0，以及隐含状态数n和观测值个数m，构建hsmm模型λ0＝(n,m,a0,b0,π0,p0)；(2)、采集一组燃料电池的全寿命电压数据作为训练数据，利用welch-baum算法更新hsmm模型中a0、b0、π0和p0的值，得到更新后的hsmm模型λ＝(n,m,a,b,π,p)；(3)、采集当前燃料电池的一段电压数据作为测试数据，将测试数据输入至训练好的hsmm模型中，利用前向算法计算当系统处于n种不同隐含状态时的p(o|λk),k＝1,2,…,n；其中，λ1＝(n,m,a,b,π1,p)，λ2＝(n,m,a,b,π2,p)，……，λn＝(n,m,a,b,πn,p)；αt(i)＝p(o1，o2…ot，st＝i|st+1≠i，λk)其中，t＝1，2，…，t，表示t个观测时刻，o＝(o1,o2,o3,…ot)表示从初始时刻t＝1到t时刻的测试数据，st＝i表示在t时刻时系统处于i状态；(4)、从n个不同的p(o|λk)中选出最大值p(o|λm)；p(o|λm)＝max(p(o|λ1),p(o|λ2),…,p(o|λn)),1≤m≤n将p(o|λm)所对应的hsmm模型标记为当前燃料电池的模型λm＝(n,m,a,b,πm,p)，当前燃料电池的健康状态处于隐含状态m；(5)、从训练数据中提取处于隐含状态m且与测试数据相同长度的观测数据on，n表示观测数据组数，将观测数据on输入至λm＝(n,m,a,b,πm,p)，通过前向后向算法分别计算出n个不同的p(oj|λm)，j＝1，2，…，n；将p(o|λm)分别与p(oj|λm)进行一一比较，选择大小最相近的一组观测数据，以该组观测数据离状态转移点的时间作为当前测试数据距离状态转移点的时间ttransfer；(6)、利用hsmm算法，计算燃料电池的剩余寿命rulhsmm：其中，amμ表示从状态m转移到状态μ的概率；(7)、根据测试数据建立经验模型：其中，a,b,c,d为待估计参数，可以通过梯度下降法估算；(8)、利用经验模型，计算燃料电池剩余寿命rulempirialrulempirial＝arg(vfailure(t))-tcurrent其中，arg(vfailure(t))表示经验模型到预设电压阈值的时间，tcurrent为测试数据的当前时间；(9)、计算训练数据和测试数据的电压梯度值gradζ(t)：gradζ(t)＝vζ(t)-vζ(t-1)其中，ζ＝1,2，v1(t)和grad1(t)分别表示t时刻训练数据的电压值和电压梯度值；v2(t)和grad2(t)分别表示t时刻测试数据的电压值和电压梯度值；(10)、计算指标β：其中，tcurrent表示当前时间；tdef为预设值，用来决定用几个单位的电压梯度值来计算相似度；(11)、估算燃料电池的剩余寿命rulfinal：rulfinal＝α·rulhsmm+(1-α)·rulempirial本发明的发明目的是这样实现的：本发明基于hsmm和经验模型的燃料电池故障预测方法，先采集一组燃料电池的全寿命电压数据，作为训练数据，再通过welch-baum算法利用训练数据训练hsmm模型；再采集一组燃料电池的电压退化数据，作为测试数据；将测试数据输入hsmm中，通过前向算法估计当前健康状态，并根据状态持续时间计算其剩余寿命；根据测试数据建立经验模型，并估计参数及预测未来走势，根据未来电压走势计算其剩余寿命；以训练数据和测试数据的电压梯度值的相似度作为标准，将hsmm得到的剩余寿命和经验模型计算出的剩余寿命结合起来，估算出燃料电池的最终剩余寿命。同时，本发明基于hsmm和经验模型的燃料电池故障预测方法还具有以下有益效果：(1)、故障预测所需数据仅为燃料电池的输出电压数据，对燃料电池内部的数据没有要求，因此在采集数据方面成本较低，较方便，也不会在采集过程中对燃料电池的内部结构造成影响。(2)、该方法相较于基于经验模型的故障预测方法，在预测速度上有很大的提高。(3)、本发明相较于目前提出的其它燃料电池故障预测方法，有较高的预测精度。附图说明图1是本发明基于hsmm和经验模型的燃料电池故障预测方法流程图；图2是hsmm训练数据状态划分结果；图3是hsmm健康状态诊断结果；图4是基于hsmm剩余寿命预测结果；图5是经验模型各个时刻预测的电压值；图6是经验模型预测的剩余寿命图7是基于hsmm和经验模型的燃料电池的预测剩余寿命。具体实施方式下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。实施例图1是本发明基于hsmm和经验模型的燃料电池故障预测方法流程图。在本实施例中，如图1所示，本发明一种基于hsmm和经验模型的燃料电池故障预测方法，包括以下步骤：s1、随机设置初始状态概率矩阵π0、初始状态转移概率矩阵a0、初始观测值概率矩阵b0、初始状态持续时间概率矩阵p0，以及隐含状态数n为3和观测值个数m为10，构建hsmm模型λ0＝(n,m,a0,b0,π0,p0)。pai0:[3×1double]a0:[3×3double]b0:[3×10double]p0:[3×300double]s2、采集一组燃料电池的全寿命电压数据作为训练数据，利用welch-baum算法更新hsmm模型中a0、b0、π0和p0的值，得到更新后的hsmm模型λ＝(n,m,a,b,π,p)；welch-baum算法是极大似然估计的一个应用，它采用了一种多次迭代优化的算法，用已知的训练数据和hsmm模型构造函数，该函数中包括了hsmm中所有的六个参数n,m,a0,b0,π0,p0，然后令各偏导为零，得到新的hsmm各参数的估计值，用新旧hsmm模型参数之间的函数关系反复进行迭代运算，直到hsmm的各个参数估计都稳定不再发生明显变化为止，得到更新后的hsmm模型λ＝(n,m,a,b,π,p)。其中更新后的模型中，另外，还有一个3*10的矩阵b和一个3*300的矩阵p，由于矩阵过大在此处就没有展示。根据状态持续时间概率矩阵p，利用公式：计算得到各个状态的平均持续时间为：根据训练结果，可作出如图2所示的状态划分图。s3、沿用上面的训练结果，采集当前燃料电池的一段电压数据作为测试数据，将测试数据输入至3个不同的hsmm模型：λ1＝(n,m,a,b,π1,p)，λ2＝(n,m,a,b,π2,p)，λ3＝(n,m,a,b,π3,p)，hsmm模型通过前向算法分别计算出3个不同的p(o|λk),k＝1,2,3，如表1所示。表1是健康状态诊断时各样本的p(o|λ)值。正常状态退化状态1退化状态2样本1-inf-inf-0.2344样本2-0.1860-5.3669-5.2402样本3-inf-0.1846-5.5593样本4-inf-inf-0.2511样本5-inf-0.1503-4.8085样本6-inf-inf-0.1732样本7-0.0700-2.6605-2.5971样本8-inf-0.1735-5.5482样本9-inf-inf-0.1876样本10-inf-0.1120-4.4118样本11-inf-inf-0.1375样本12-0.1001-3.8007-3.7101样本13-inf-0.0560-2.2059样本14-inf-inf-0.0952表1s4、以第二个样本为例，p(o|λ1)为-0.1860，p(o|λ2)也为-5.3669，p(o|λ3)为-5.2402。因p(o|λ1)的值最大，我们认为此时的燃料电池系统处于第一个隐含状态。在本实施例中，我们取了14组数据来做当前健康状态的估计，估计结果如图3所示，从图3中可以看出，对于实验所选取的14个样本点，hsmm都能精确的估计其健康状态。s5、以第二个样本为例，该样本处于第一个隐含状态，根据上述
发明内容中第五步描述的方法，计算得到的ttransfer为20。s6、利用hsmm算法，计算燃料电池的剩余寿命rulhsmm。以第二个样本为例，该样本处于第一个隐含状态，根据上述式子，可计算得到：rul3＝72hrul2＝44h+1*rul3＝116hrulhsmm＝20h+1*rul2+0*rul3＝136h根据该方法计算所得的rulhsmm如图4所示。图中，横坐标表示的是进行寿命预测时的时刻；纵坐标代表各时刻寿命预测的结果。粗体实线表示的是hsmm对于阳极氯气浓度为的燃料电池寿命预测曲线；粗体虚线为各个时刻所对应的真实寿命值；上侧细虚线为真实寿命值的置信上限；下侧细点线为真实寿命值的置信下限，两者的置信概率均为15％。s7、在本实施例中，根据测试数据设置50组不同的初始值，限制voltage(t)是关于时间t的单调递减函数，并且限制在一定时间内电压不会退化到失效，再用梯度下降法分别估计出参数a,b,c,d的值；然后根据估计出的参数值，分别预测燃料电池电压在未来的走势，取这50组不同预测曲线中的中位数所对应的模型，作为我们所求的经验模型。以燃料电池前10个小时运行的电压数据为例，根据梯度下降法训练得到的经验模型模型为：s8、利用经验模型，计算燃料电池剩余寿命rulempirialrulempirial＝arg(vfailure(t))-tcurrent其中，arg(vfailure(t))表示经验模型到预设电压阈值的时间，tcurrent为测试数据的当前时间；以燃料电池前10个小时运行的电压数据为例，我们预设的电压阈值为0.89v，测试数据的当前时间为10小时，根据上述式子计算可得：rulempirial＝144h-10h＝134h在10、20、……140时刻所选的经验模型电压预测结果如图5所示。经验模型的剩余寿命预测结果如图6所示。s9、根据上述
发明内容中第九步展示的电压梯度计算公式，分别计算hsmm训练数据和当前测试数据的电压梯度值。s10、同样以燃料电池前10个小时运行的电压数据为例，设置tdef＝10，根据上述
发明内容中第十步展示的β计算公式，可得β＝0.26。其中，指标β用来衡量训练数据和测试数据之间的相似度，β值越小，则相似度越高；β值越大，相似度越低。s11、同样以燃料电池前10个小时运行的电压数据为例，根据上述
发明内容中第十一步展示的α计算公式及rulfinal，此时由于β＝0.26，那么α＝1，因此rulfinal＝rulhsmm＝136h。根据该结合方法得到的最终剩余寿命结果如图7所示。尽管上面对本发明说明性的具体实施方式进行了描述，以便于本
技术领域：
的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本
技术领域：
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。当前第1页12