CN116387569A - 一种基于svm的全钒液流电池泵故障检测方法 - Google Patents
一种基于svm的全钒液流电池泵故障检测方法 Download PDFInfo
- Publication number
- CN116387569A CN116387569A CN202310160420.5A CN202310160420A CN116387569A CN 116387569 A CN116387569 A CN 116387569A CN 202310160420 A CN202310160420 A CN 202310160420A CN 116387569 A CN116387569 A CN 116387569A
- Authority
- CN
- China
- Prior art keywords
- pump
- classification
- data
- redox flow
- hyperplane
- 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.)
- Pending
Links
- 229910052720 vanadium Inorganic materials 0.000 title claims abstract description 47
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 238000012549 training Methods 0.000 claims abstract description 26
- 238000012795 verification Methods 0.000 claims abstract description 24
- 238000012706 support-vector machine Methods 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 230000006870 function Effects 0.000 claims description 67
- 238000000034 method Methods 0.000 claims description 45
- 230000009977 dual effect Effects 0.000 claims description 21
- 230000008859 change Effects 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000004146 energy storage Methods 0.000 abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 3
- 238000003745 diagnosis Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000002572 peristaltic effect Effects 0.000 description 5
- 238000007599 discharging Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012983 electrochemical energy storage Methods 0.000 description 2
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 239000007772 electrode material Substances 0.000 description 1
- 239000003792 electrolyte Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 229910052744 lithium Inorganic materials 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- LEONUFNNVUYDNQ-UHFFFAOYSA-N vanadium atom Chemical compound [V] LEONUFNNVUYDNQ-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/18—Regenerative fuel cells, e.g. redox flow batteries or secondary fuel cells
- H01M8/184—Regeneration by electrochemical means
- H01M8/188—Regeneration by electrochemical means by recharging of redox couples containing fluids; Redox flow type batteries
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04313—Processes for controlling fuel cells or fuel cell systems characterised by the detection or assessment of variables; characterised by the detection or assessment of failure or abnormal function
- H01M8/0438—Pressure; Ambient pressure; Flow
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04313—Processes for controlling fuel cells or fuel cell systems characterised by the detection or assessment of variables; characterised by the detection or assessment of failure or abnormal function
- H01M8/04664—Failure or abnormal function
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04694—Processes for controlling fuel cells or fuel cell systems characterised by variables to be controlled
- H01M8/04746—Pressure; Flow
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04992—Processes for controlling fuel cells or fuel cell systems characterised by the implementation of mathematical or computational algorithms, e.g. feedback control loops, fuzzy logic, neural networks or artificial intelligence
-
- 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/50—Fuel cells
Abstract
本发明涉及一种基于SVM的全钒液流电池泵故障检测方法,适用于电网储能系统,该方法不需要任何流量传感器,可对全钒液流电池泵的正极泵故障、负极泵故障及双侧泵故障进行分类。从电池状态曲线中提取特征参数,对支持向量机的参数进行了优化,最后通过支持向量机训练得到故障预测结果,实现对全钒液流电池泵的正极泵故障、负极泵故障及双侧泵故障进行故障分类。采用支持向量机算法检测全钒液流电池水泵在实际运行中的运行状态,并利用交叉验证算法优化支持向量机中的参数,能够很好地解决全钒液流电池泵故障检测的技术问题。
Description
技术领域
本发明涉及全钒液流电池技术领域,尤其涉及一种基于SVM的全钒液流电池泵故障检测方法。
背景技术
随着应用市场对储能系统可靠性和容量稳定性的需求逐年增加,用于电网应用的储能系统的部署速度正在加快。众所周知,电化学储能技术是发展固定式电力储能系统的关键一环。在众多的电化学储能技术中,全钒液流电池因具有安全性高、输出功率和容量相互独立、过载能力强、深放电能力强和循环寿命长等优点,而具有广泛的商业应用价值。
电池的故障诊断在现代工业系统中发挥着越来越重要的作用。随着大数据处理技术的快速发展,数据驱动方法在储能电池故障诊断领域得到了日益广泛的应用。在这些电池的故障诊断方法中,常用的方法有支持向量机(Support Vector Machine,SVM)。在与故障诊断相关的应用中,故障样本的数量通常比正常样本的数量小得多。在这种样本数不平衡的情况下,SVM只需少量学习样本就可以提供更可靠的结果。此外,SVM的结果通常可以用一小部分样本(即支持向量)表示,这使得SVM可以在执行较小的计算量下保证全局最优解。
SVM主要处理的问题的类型有两类,分别是线性可分问题和线性不可分问题。线性可分问题相对简单,SVM可以直接对其求解,并寻找最优的分类效果;而线性不可分问题的计算相对复杂,需要在分类之前使用一个适当的转换方式将低维的线性不可分问题转化为高维的线性可分问题之后再对其进行求解,因此转换方式的选取在解决线性不可分问题上至关重要。
在现有的研究中,SVM方法主要用于锂电池的故障诊断,而针对全钒液流电池的研究大多集中在电池性能分析、电极材料优化和流动路径优化上,很少有研究集中在全钒液流电池的故障诊断上。
发明内容
本发明的目的是提供一种基于SVM的全钒液流电池泵故障检测方法,将SVM方法应用于全钒液流电池的故障诊断,该方法不需要任何流量传感器,可对全钒液流电池泵的正极泵故障、负极泵故障及双侧泵故障进行分类。采用支持向量机算法检测全钒液流电池水泵在实际运行中的运行状态,并利用交叉验证算法优化支持向量机中的参数,能够很好地解决全钒液流电池泵故障检测的技术问题。
为实现上述目的,本发明所设计的一种基于SVM的全钒液流电池泵故障检测方法,包括以下步骤:
步骤1:对全钒液流电池分别进行正极泵、负极泵和双侧泵的模拟故障实验,泵故障模拟方式具体为,充放电过程中停泵数次,停泵过程即视为发生了泵故障,且充放电总时间大于停泵总时间,之后再将电池充满或充分放电;记录全钒液流电池分别在正极泵故障、负极泵故障或双侧泵故障的情况下的实验数据;
步骤2:从步骤1的实验数据中提取电流与电压信号的变化曲线,在全钒液流电池发生泵故障的情况下,电压信号的变化是最直接和最明显的,为了消除电流变化对电压的影响,引入电阻变化率ΔR以直观地反映电池的状态,并用式(1)表示:
其中Ut为某一时刻采集的电压值;Ut-1为某一时刻的前一秒内采集的电压值;It为某一时刻采集的电流值;
步骤3:将电池的荷电状态SOC与电阻变化率ΔR作为样本特征输入支持向量机中训练,利用对偶函数和拉格朗日优化算法,结合核函数,建立故障分类器模型,并用式(2)、式(3)表示:
式(2)为非线性SVM的分类决策函数,式(3)为径向基核函数,核函数将有限维数据映射到更高维空间,用于解决线性不可分割问题;
其中f(x)为分类决策函数,x为样本数据,xi为样本数据点,yi为数据类别,ai为拉格朗日乘子,满足条件ai≥0;b为偏移项;K(x,xi)为径向基核函数,满足Mercer定理;γ为径向基核函数的特征参数,反映了训练样本的特征;
步骤4:结合交叉验证和网格搜索的方法,针对步骤3的故障分类器模型中的特征参数进行优化;交叉验证方法将原始数据划分为训练集和验证集,训练分类器后,验证集的准确性可以反映分类器的准确性;网格搜索方法在特定范围内构建网格,并结合交叉验证方法得到网格中分类精度最高的惩罚因子C和核函数参数γ;
步骤5:将优化后的惩罚因子C和核函数参数γ输入到步骤3的故障分类器模型中进行训练,得到最终的故障分类结果。
作为优选方案,所述步骤3包括通过故障分类器模型找到相应的超平面,根据超平面对数据进行分类。
作为优选方案,通过故障分类器模型找到相应的超平面,具体包括以下步骤,
步骤301:将从步骤2的实验数据中提取的样本数据作为训练数据集(xi,yi;i=1,...,n)输入故障分类器模型;
步骤302:故障分类器模型按照区间最大化的原理找到相应的超平面,由超平面对这些数据进行分类,超平面的方程由式(3-1)表示:
wTx+b=0 (3-1)
其中w为参数向量;T代表转置;x为样本数据;b为偏移项,决定了超平面和原点之间的距离;
步骤303:确定超平面后,用|wTx+b|表示从点x到超平面的距离,通过观察wTx+b的符号与数据类别y的符号是否一致来判断分类是否正确,即可以用y(wTx+b)的正负性来判定或表示分类的正确性,定义函数间隔为:
步骤304:对于一个样本数据点x0,x0为x垂直投影到超平面上的对应点;w为垂直于超平面的一个向量;γ为样本x到分类间隔的距离;则x可以表示为将x0代入超平面的方程,即可算出样本x到分类间隔的距离γ,由式(3-4)表示;令γ乘上对应的数据类别y,即可得出几何间隔由式(3-5)表示:
其中||w||表示的是范数。
作为优选方案,根据超平面对数据进行分类,具体包括以下步骤,
其中,令函数间隔根据式(3-7)、(3-8)可得式(3-9),式(3-9)即为最大间隔分类器的目标函数;由于求1/||w||的最大值相当于求1/2||w||2的最小值,故将式(3-9)转化为等价问题表达式(3-10);
步骤306:通过拉格朗日对偶性,将目标函数变换到对偶变量的优化问题,即通过求解与原问题等价的对偶问题得到原始问题的最优解,通过给每一个约束条件加上一个拉格朗日乘子α定义拉格朗日函数:
其中,p*表示这个问题的最优值,且和最初的问题等价;在要求约束条件得到满足的情况下最小化1/2||w||2实际上等价于直接最小化θ(w),则目标函数可由式(3-13)表示;将式(3-13)求解最小和最大的顺序交换,得到原始问题的对偶问题,由式(3-14)表示,d*为对偶问题的最优值,且有d*≤p*;
步骤307:通过求解对偶问题来间接地求解原始问题,结合核函数,得到最终的分类决策函数式(2):
其中,核函数相当于将最初的分类函数式(3-15)映射成式(3-16);而其中的α可以通过求解式(3-17)中的对偶问题得到。
作为优选方案,所述步骤4具体包括以下步骤:
步骤401:为了排除数据中噪声的影响,结合软区间最大化原理,SVM允许数据点在一定程度上偏离超平面:
式(4-1)为非线性支持向量机的分类决策函数;其中ξi为松弛变量;C为惩罚因子,用于维持分离误差和计算复杂度之间适当的平衡;
步骤402:将原始数据随机分为两组,一组作为训练集,另一组作为验证集,利用训练集训练分类器,然后利用验证集验证模型,记录最后的分类准确率;
步骤403:将待搜索参数惩罚因子C和核函数参数γ对应的取值范围看作一个空间,并将此参数空间按照规定的步长划分成网格,随后在每个网格点处测试样本,通过遍历网格中所有的点来寻找最优参数。
作为优选方案,步骤3中样本数据为全钒液流电池的个数,数据类别包括实际分类、预测分类,实际分类与预测分类符号重合即为分类正确,二者不重合即为分类不正确;预测分类准确率的计算方法为(正确分类样本数/总样本数)*100%。
作为优选方案,所述步骤1中,以充电过程为例进行泵故障模拟实验,分别在三种电流下对电池进行充电,且每种电流下分别停泵两次,即先充4分钟,停泵30秒;再充4分钟,再停泵30秒;之后将电池充满。
本发明的有益效果:
该检测方法将基于网格搜索的SVM算法应用于全钒液流电池的故障诊断,因网格搜索是一种优化算法的方法,在SVM算法中应用比较普遍。该方法适用于电网储能系统,该方法不需要任何流量传感器,可对全钒液流电池泵的正极泵故障、负极泵故障及双侧泵故障进行分类。从电池状态曲线中提取特征参数,对支持向量机的参数进行了优化,最后通过支持向量机训练得到故障预测结果,实现对全钒液流电池泵的正极泵故障、负极泵故障及双侧泵故障进行故障分类。采用支持向量机算法检测全钒液流电池水泵在实际运行中的运行状态,并利用交叉验证算法优化支持向量机中的参数,能够很好地解决全钒液流电池泵故障检测的技术问题。
附图说明
图1为本发明的流程图;
图2为本发明参数优化结果示意图;
图3为本发明预测分类示意图。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚,下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。可以理解的是,此处所描述的具体实施例仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、仅用于描述目的,而不能理解为指示或暗示相对重要性。其中,术语“第一位置”和“第二位置”为两个不同的位置。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
本发明涉及一种基于SVM的全钒液流电池泵故障检测方法,适用于电网储能系统,该方法不需要任何流量传感器,从电池状态曲线中提取特征参数,对支持向量机的参数进行了优化,最后通过支持向量机训练得到故障预测结果,实现对全钒液流电池泵的正极泵故障、负极泵故障及双侧泵故障进行故障分类。
如图1所示,本发明提出的一种基于SVM的全钒液流电池泵故障检测方法,按如下步骤进行。
步骤1:对全钒液流电池分别进行正极泵、负极泵和双侧泵的模拟故障实验:
(i)先用电流It对电池充分放电;
(ii)再用电流It对电池充电;
(iii)将正极泵关闭,模拟正极泵故障;继续对电池进行充电,直至将电池充满;
(iv)将负极泵关闭,模拟负极泵故障;继续对电池进行充电,直至将电池充满;
(v)将双侧泵同时关闭,模拟双侧泵故障;继续对电池进行充电,直至将电池充满;
(vi)记录全钒液流电池分别在正极泵故障、负极泵故障或双侧泵故障的情况下的实验数据;
本实例中全钒液流电池的实验参数设置如表1所示,以充电过程为例进行泵故障模拟实验,分别在It为500mA、1000mA、1500mA这三种电流下对电池进行充电,且每种电流下分别停泵两次,即先充4分钟,停泵30秒;再充4分钟,再停泵30秒;之后将电池充满。放电过程和充电过程正好相反,不再赘述。在进行全钒液流电池泵故障模拟实验中,充电时,停泵过程电池电压明显升高;放电时,停泵过程电池电压明显降低。
表1全钒液流电池的实验参数设置
实验参数 | 设置值 |
初始电解液体积/mL | 20mL |
充放电电流/mA | 500、1000、1500 |
放电截止电压/V | 1 |
充电截止电压/V | 1.7 |
蠕动泵流速/mL/min | 80 |
步骤2:从步骤1的实验数据中提取电流与电压信号的变化曲线。由于电压信号受全钒液流电池的电化学性质、温度、流速和材料性质的综合影响。因此,电压信号是全钒液流电池健康检测指标中的关键信号之一。在发生故障的情况下,电压信号的变化是最直接和最明显的。因为电流变化也会导致电压变化,电池里面有内阻,电流升高,电阻上的分压就会变大,导致整体电压也会升高。因此,为了消除电流变化对电压的影响,引入电阻变化率ΔR以直观地反映电池的状态,并用式(1)表示:
其中Ut为某一时刻采集的电压值;Ut-1为某一时刻的前一秒内采集的电压值;It为某一时刻采集的电流值。
步骤3:将电池的荷电状态SOC与电阻变化率ΔR作为样本特征输入支持向量机中训练,利用对偶函数和拉格朗日优化算法,结合核函数,建立故障分类器模型,并用式(2)-式(3)表示:
式(2)为非线性SVM的分类决策函数,式(3)为径向基核函数,核函数将有限维数据映射到更高维空间,用于解决线性不可分割问题;
其中f(x)为分类决策函数,x为样本数据,xi为数据点,yi为数据类别,ai为拉格朗日乘子,满足条件ai≥0;b为偏移项;K(x,xi)为径向基核函数,满足Mercer定理;γ为径向基核函数的特征参数,反映了训练样本的特征;
步骤301:将从步骤2的实验数据中提取的样本数据作为训练数据集(xi,yi;i=1,...,n)输入故障分类器模型。
步骤302:故障分类器模型按照区间最大化的原理找到相应的超平面,由超平面对这些数据进行分类,超平面的方程由式(3-1)表示:
wTx+b=0 (3-1)
其中w为参数向量;T代表转置;x为样本数据;b为偏移项,决定了超平面和原点之间的距离;
步骤303:确定超平面后,用|wTx+b|表示从点x到超平面的距离,通过观察wTx+b的符号与数据类别y的符号是否一致来判断分类是否正确,即可以用y(wTx+b)的正负性来判定或表示分类的正确性,定义函数间隔为:
步骤304:对于一个样本数据点x0,x0为x垂直投影到超平面上的对应点;w为垂直于超平面的一个向量;γ为样本x到分类间隔的距离;则x可以表示为将x0代入超平面的方程,即可算出样本x到分类间隔的距离γ,由式(3-4)表示;令γ乘上对应的数据类别y,即可得出几何间隔由式(3-5)表示:
其中||w||表示的是范数;
其中,令函数间隔根据式(3-7)、(3-8)可得式(3-9),式(3-9)为最大间隔分类器的目标函数;由于求1/||w||的最大值相当于求1/2||w||2的最小值,故将式(3-9)转化为等价问题表达式(3-10);
步骤306:通过拉格朗日对偶性,将目标函数变换到对偶变量的优化问题,即通过求解与原问题等价的对偶问题得到原始问题的最优解。通过给每一个约束条件加上一个拉格朗日乘子α定义拉格朗日函数:
其中,p*表示这个问题的最优值,且和最初的问题等价;在要求约束条件得到满足的情况下最小化1/2||w||2实际上等价于直接最小化θ(w),则目标函数可由式(3-13)表示;将式(3-13)求解最小和最大的顺序交换,得到原始问题的对偶问题,由式(3-14)表示,d*为对偶问题的最优值,且有d*≤p*;
步骤307:通过求解对偶问题来间接地求解原始问题,结合核函数,得到最终的分类决策函数式(2):
其中,核函数相当于将最初的分类函数式(3-15)映射成式(3-16);而其中的α可以通过求解式(3-17)中的对偶问题得到;
步骤4:结合交叉验证和网格搜索的方法,针对步骤3的故障分类器模型中的特征参数进行优化;交叉验证方法将原始数据划分为训练集和验证集,训练分类器后,验证集的准确性可以反映分类器的准确性;网格搜索方法在特定范围内构建网格,并结合交叉验证方法得到网格中分类精度最高的惩罚因子C和核函数参数γ;
步骤401:为了排除数据中噪声的影响,结合软区间最大化原理,SVM允许数据点在一定程度上偏离超平面:
式(4-1)为非线性支持向量机的分类决策函数;其中ξi为松弛变量;C为惩罚因子,用于维持分离误差和计算复杂度之间适当的平衡。
步骤402:将原始数据随机分为两组,一组作为训练集,另一组作为验证集,利用训练集训练分类器,然后利用验证集验证模型,记录最后的分类准确率;
步骤403:将待搜索参数惩罚因子C和核函数参数γ对应的取值范围看作一个空间,并将此参数空间按照规定的步长划分成网格,随后在每个网格点处测试样本,通过遍历网格中所有的点来寻找最优参数。
本实施例中,惩罚因子C和核函数参数γ的取值范围设定为[2-8,28];步长间隔大小设为1;最终获得的最佳特征参数如图2所示,惩罚因子C=0.25,核函数参数γ=256;在交叉验证方法下的最佳分类精度为99.3243%;
步骤5:将优化后的惩罚因子C和核函数参数γ输入到步骤3的故障分类器模型中进行训练,得到最终的故障分类结果,如图3所示。
图3为在正极泵故障情况下的分类准确率示意图,图中横坐标为样本数据x,样本数据0-140分别代表样本个数,即全钒液流电池的个数。纵坐标为数据类别y,图中“○”代表样本点的实际分类,“*”代表预测分类,实际分类与预测分类符号重合即为分类正确,二者不重合即为分类不正确;标签“1”代表蠕动泵的正常运行状态,标签“-1”代表蠕动泵的故障状态;从图中可以看出,在蠕动泵正常运行状态下,重合率为100%;在蠕动泵处于故障状态下,有部分不重合。
预测分类准确率的计算方法为:(正确分类样本数/总样本数)*100%。输入实验参数,根据SVM算法得知,在正极泵故障情况下的分类准确率为99.1935%。在双侧泵故障和负极泵故障情况下的分类准确率示意图与图3相似,不再示意。同样,根据SVM算法得知在双侧泵故障和负极泵故障情况下的分类准确率分别为100%和98.3871%。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (7)
1.一种基于SVM的全钒液流电池泵故障检测方法,其特征在于,包括以下步骤:
步骤1:对全钒液流电池分别进行正极泵、负极泵和双侧泵的模拟故障实验,泵故障模拟方式具体为,充放电过程中停泵数次,停泵过程即视为发生了泵故障,且充放电总时间大于停泵总时间,之后再将电池充满或充分放电;记录全钒液流电池分别在正极泵故障、负极泵故障或双侧泵故障的情况下的实验数据;
步骤2:从步骤1的实验数据中提取电流与电压信号的变化曲线,在全钒液流电池发生泵故障的情况下,电压信号的变化是最直接和最明显的,为了消除电流变化对电压的影响,引入电阻变化率ΔR以直观地反映电池的状态,并用式(1)表示:
其中Ut为某一时刻采集的电压值;Ut-1为某一时刻的前一秒内采集的电压值;It为某一时刻采集的电流值;
步骤3:将电池的荷电状态SOC与电阻变化率ΔR作为样本特征输入支持向量机中训练,利用对偶函数和拉格朗日优化算法,结合核函数,建立故障分类器模型,并用式(2)、式(3)表示:
式(2)为非线性SVM的分类决策函数,式(3)为径向基核函数,核函数将有限维数据映射到更高维空间,用于解决线性不可分割问题;
其中f(x)为分类决策函数,x为样本数据,xi为样本数据点,yi为数据类别,ai为拉格朗日乘子,满足条件ai≥0;b为偏移项;K(x,xi)为径向基核函数,满足Mercer定理;γ为径向基核函数的特征参数,反映了训练样本的特征;
步骤4:结合交叉验证和网格搜索的方法,针对步骤3的故障分类器模型中的特征参数进行优化;交叉验证方法将原始数据划分为训练集和验证集,训练分类器后,验证集的准确性可以反映分类器的准确性;网格搜索方法在特定范围内构建网格,并结合交叉验证方法得到网格中分类精度最高的惩罚因子C和核函数参数γ;
步骤5:将优化后的惩罚因子C和核函数参数γ输入到步骤3的故障分类器模型中进行训练,得到最终的故障分类结果。
2.根据权利要求1所述的基于SVM的全钒液流电池泵故障检测方法,其特征在于:所述步骤3包括通过故障分类器模型找到相应的超平面,根据超平面对数据进行分类。
3.根据权利要求2所述的基于SVM的全钒液流电池泵故障检测方法,其特征在于:通过故障分类器模型找到相应的超平面,具体包括以下步骤,
步骤301:将从步骤2的实验数据中提取的样本数据作为训练数据集(xi,yi;i=1,…,n)输入故障分类器模型;
步骤302:故障分类器模型按照区间最大化的原理找到相应的超平面,由超平面对这些数据进行分类,超平面的方程由式(3-1)表示:
wTx+b=0 (3-1)
其中w为参数向量;T代表转置;x为样本数据;b为偏移项,决定了超平面和原点之间的距离;
步骤303:确定超平面后,用|wTx+b|表示从点x到超平面的距离,通过观察wTx+b的符号与数据类别y的符号是否一致来判断分类是否正确,即可以用y(wTx+b)的正负性来判定或表示分类的正确性,定义函数间隔为:
步骤304:对于一个样本数据点x0,x0为x垂直投影到超平面上的对应点;w为垂直于超平面的一个向量;γ为样本x到分类间隔的距离;则x可以表示为将x0代入超平面的方程,即可算出样本x到分类间隔的距离γ,由式(3-4)表示;令γ乘上对应的数据类别y,即可得出几何间隔由式(3-5)表示:
其中||w||表示的是范数。
4.根据权利要求2所述的基于SVM的全钒液流电池泵故障检测方法,其特征在于:根据超平面对数据进行分类,具体包括以下步骤,
s.t.yi(wxi+b)≥1,i=1,2,…,n
其中,令函数间隔根据式(3-7)、(3-8)可得式(3-9),式(3-9)即为最大间隔分类器的目标函数;由于求1/||w||的最大值相当于求1/2||w||2的最小值,故将式(3-9)转化为等价问题表达式(3-10);
步骤306:通过拉格朗日对偶性,将目标函数变换到对偶变量的优化问题,即通过求解与原问题等价的对偶问题得到原始问题的最优解,通过给每一个约束条件加上一个拉格朗日乘子α定义拉格朗日函数:
其中,p*表示这个问题的最优值,且和最初的问题等价;在要求约束条件得到满足的情况下最小化1/2||w||2实际上等价于直接最小化θ(w),则目标函数可由式(3-13)表示;将式(3-13)求解最小和最大的顺序交换,得到原始问题的对偶问题,由式(3-14)表示,d*为对偶问题的最优值,且有d*≤p*;
步骤307:通过求解对偶问题来间接地求解原始问题,结合核函数,得到最终的分类决策函数式(2):
其中,核函数相当于将最初的分类函数式(3-15)映射成式(3-16);而其中的α可以通过求解式(3-17)中的对偶问题得到。
5.根据权利要求1所述的基于SVM的全钒液流电池泵故障检测方法,其特征在于:所述步骤4具体包括以下步骤:
步骤401:为了排除数据中噪声的影响,结合软区间最大化原理,SVM允许数据点在一定程度上偏离超平面:
式(4-1)为非线性支持向量机的分类决策函数;其中ξi为松弛变量;C为惩罚因子,用于维持分离误差和计算复杂度之间适当的平衡;
步骤402:将原始数据随机分为两组,一组作为训练集,另一组作为验证集,利用训练集训练分类器,然后利用验证集验证模型,记录最后的分类准确率;
步骤403:将待搜索参数惩罚因子C和核函数参数γ对应的取值范围看作一个空间,并将此参数空间按照规定的步长划分成网格,随后在每个网格点处测试样本,通过遍历网格中所有的点来寻找最优参数。
6.根据权利要求1所述的基于SVM的全钒液流电池泵故障检测方法,其特征在于:步骤3中样本数据为全钒液流电池的个数,数据类别包括实际分类、预测分类,实际分类与预测分类符号重合即为分类正确,二者不重合即为分类不正确;预测分类准确率的计算方法为(正确分类样本数/总样本数)*100%。
7.根据权利要求1所述的基于SVM的全钒液流电池泵故障检测方法,其特征在于:所述步骤1中,以充电过程为例进行泵故障模拟实验,分别在三种电流下对电池进行充电,且每种电流下分别停泵两次,即先充4分钟,停泵30秒;再充4分钟,再停泵30秒;之后将电池充满。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310160420.5A CN116387569A (zh) | 2023-02-23 | 2023-02-23 | 一种基于svm的全钒液流电池泵故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310160420.5A CN116387569A (zh) | 2023-02-23 | 2023-02-23 | 一种基于svm的全钒液流电池泵故障检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116387569A true CN116387569A (zh) | 2023-07-04 |
Family
ID=86960504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310160420.5A Pending CN116387569A (zh) | 2023-02-23 | 2023-02-23 | 一种基于svm的全钒液流电池泵故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116387569A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117291124A (zh) * | 2023-11-24 | 2023-12-26 | 艾博特瑞能源科技(苏州)有限公司 | 一种全钒液流电池性能预测方法 |
-
2023
- 2023-02-23 CN CN202310160420.5A patent/CN116387569A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117291124A (zh) * | 2023-11-24 | 2023-12-26 | 艾博特瑞能源科技(苏州)有限公司 | 一种全钒液流电池性能预测方法 |
CN117291124B (zh) * | 2023-11-24 | 2024-02-06 | 艾博特瑞能源科技(苏州)有限公司 | 一种全钒液流电池性能预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Battery fault diagnosis for electric vehicles based on voltage abnormality by combining the long short-term memory neural network and the equivalent circuit model | |
CN110752410B (zh) | 一种快速分选和重组退役锂电池的方法 | |
CN111090050A (zh) | 一种基于支持向量机和k均值的锂电池故障诊断方法 | |
CN108896914A (zh) | 一种锂电池健康状况的梯度提升树建模与预测方法 | |
US11830991B2 (en) | Method and apparatus based on neural network for pre-diagnosing defect and fire in battery cell | |
CN107422266B (zh) | 一种大容量电池储能系统的故障诊断方法及装置 | |
CN112630660B (zh) | 一种基于支持向量机的电池故障识别方法 | |
CN110376530A (zh) | 电池内短路检测装置及方法 | |
CN114899457B (zh) | 一种针对质子交换膜燃料电池系统的故障检测方法 | |
CN116387569A (zh) | 一种基于svm的全钒液流电池泵故障检测方法 | |
CN112287980B (zh) | 基于典型特征向量的动力电池筛选方法 | |
CN114565008A (zh) | 基于长短时记忆自编码器的锂电池异常检测方法 | |
CN112305441B (zh) | 一种集成式聚类下的动力电池健康状态评估方法 | |
CN115327403A (zh) | 基于新能源汽车大数据的动力电池安全性评价方法及装置 | |
CN113866641A (zh) | 一种锂离子电池的故障检测方法及装置 | |
CN114460481A (zh) | 基于Bi-LSTM和注意力机制的储能电池热失控预警方法 | |
Liu et al. | Online diagnosis and prediction of power battery voltage comprehensive faults for electric vehicles based on multi-parameter characterization and improved K-means method | |
Quan et al. | An enhanced fault diagnosis method for fuel cell system using a kernel extreme learning machine optimized with improved sparrow search algorithm | |
CN116449218B (zh) | 一种锂电池健康状态的估计方法 | |
CN115825756B (zh) | 分散式储能电站故障多级诊断方法及系统 | |
CN116893366A (zh) | 用于运行系统以识别设备的电蓄能器的异常的方法和装置 | |
CN115656837A (zh) | 一种串联型电池故障预测方法 | |
CN113391214A (zh) | 一种基于电池充电电压排名变化的电池微故障诊断方法 | |
CN114595742A (zh) | 一种燃料电池故障数据采样方法及系统 | |
CN112305442A (zh) | 一种基于kmeans聚类的动力电池SOH快速估计方法 |
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 |