CN113343471B - 随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统 - Google Patents
随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统 Download PDFInfo
- Publication number
- CN113343471B CN113343471B CN202110680527.3A CN202110680527A CN113343471B CN 113343471 B CN113343471 B CN 113343471B CN 202110680527 A CN202110680527 A CN 202110680527A CN 113343471 B CN113343471 B CN 113343471B
- Authority
- CN
- China
- Prior art keywords
- fuel cell
- data
- impedance spectrum
- electrochemical
- equation
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Fuel Cell (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及一种随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统,包括数据记录步骤:在燃料电池衰减实验中,以T为时间周期,记录待测数据的数值变化;方程建立步骤:建立待测数据与衰减测试时间的方程;参数计算验证步骤:基于建立的方程计算验证方程的参数;待测数据预估步骤:根据计算得到参数的估计值,预估待测数据随测试时间的变化;燃料电池电化学交流阻抗谱预估步骤:基于预估的待测数据,预估燃料电池电化学交流阻抗谱。本方法通过燃料电池衰减过程中累积的高频、低频测试结果,基于贝叶斯原理及等效电路模型预估交流阻抗谱变化。与现有技术相比,本发明精度高、计算结果可靠性高、同时可以降低燃料电池寿命测试所需时间。
Description
技术领域
本发明涉及燃料电池领域,具体地,涉及随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统。
背景技术
燃料电池在运行过程中类似一个黑盒,电化学交流阻抗谱作为一种有效的电化学测试手段,可以为燃料电池运行状态提供多维的数据,用于评估燃料电池当前运行状态以及剩余寿命预测。
电化学交流阻抗谱在测试过程中需要消耗额外的能源用于提供交流信号的输入,并且长时间的测试会对燃料电池造成一定的不可逆损伤。有效的电化学交流阻抗谱模型对于燃料电池的电化学交流阻抗谱测试至关重要,可以有效地减少测试时间,减低对燃料电池的损伤,提高燃料电池寿命。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统。
根据本发明提供的一种随燃料电池衰减变化的电化学交流阻抗谱预测方法,包括如下步骤:
数据记录步骤:在燃料电池衰减实验中,以T为时间周期,记录待测数据的数值变化;
方程建立步骤:建立待测数据与衰减测试时间的方程;
参数计算验证步骤:基于建立的方程计算验证方程的参数;
待测数据预估步骤:根据计算得到参数的估计值,预估待测数据随测试时间的变化;
燃料电池电化学交流阻抗谱预估步骤:基于预估的待测数据,预估燃料电池电化学交流阻抗谱。
优选地,还包括偏差预估步骤:基于高斯分布的累积分布函数,计算预估值在设定的置信程度下的取值区间,得到概率函数。
优选地,所述待测数据包括高频内阻rohm、高频阻抗弧半径Rhf以及低频阻抗弧半径Rlf。
优选地,待测数据与衰减测试时间的方程为:r=a2x2+a1x+a0+ε,其中ε为噪声均方差,a2,a1,a0为多项式系数,仪器噪声满足高斯分布N(0,ε2),取a=(a2,a1,a0)T,x=(x2,x,1)T,则r=aT*x+ε。
优选地,所述参数计算验证步骤包括:
后验分布、先验分布、似然函数间满足:
均值计算公式为:
方差计算公式为:
优选地,时间周期T为1-100h。
优选地,燃料电池电化学交流阻抗谱预估步骤包括:基于得到预估的待测数据,根据燃料电池的Randel电路,模拟相应的奈奎斯特图。
根据本发明提供的一种随燃料电池衰减变化的电化学交流阻抗谱预测系统,包括如下模块:
数据记录模块:在燃料电池衰减实验中,以T为时间周期,记录待测数据的数值变化;
方程建立模块:建立待测数据与衰减测试时间的方程;
参数计算验证模块:基于建立的方程计算验证方程的参数;
待测数据预估模块:根据计算得到参数的估计值,预估待测数据随测试时间的变化;
燃料电池电化学交流阻抗谱预估模块:基于预估的待测数据,预估燃料电池电化学交流阻抗谱。
与现有技术相比,本发明具有如下的有益效果:
1、本发明精度高、计算结果可靠性高、同时可以降低燃料电池寿命测试所需时间。
2、本发明通过燃料电池衰减过程中累积的高频、低频测试结果,基于贝叶斯原理及等效电路模型预估交流阻抗谱变化。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为随燃料电池衰减变化的电化学交流阻抗谱预测方法步骤流程图。
图2为随燃料电池衰减变化的电化学交流阻抗谱预测方法对应设备原理图。
图4为高频阻抗弧半径Rhf随测试时间t变化示意图。
图5为低频阻抗弧半径Rlf随测试时间t变化示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示,本发明提供了一种基于贝叶斯原理及等效电路模型预估的燃料电池交流阻抗谱衰减预估的方法,包含以下步骤:
(1)燃料电池衰减测试过程中以周期T进行电化学交流阻抗谱测试,记录高频内阻rohm、高频阻抗弧半径Rhf以及低频阻抗弧半径Rlf,一般地,测量周期T取1~100h;
(2)高频内阻rohm、高频阻抗弧半径Rhf以及低频阻抗弧半径Rlf随电池运行时间衰减认为符合多项式关系,特征的,取二次多项式:r=a2x2+a1x+a0,其中r可以为rohm、Rhf、Rlf中的任意一项,x为衰减测试的时间,a为系数。考虑到仪器测量噪声,修正上述二次多项式为r=a2x2+a1x+a0+ε,其中ε为噪声均方差,仪器噪声满足高斯分布N(0,ε2),取a=(a2,a1,a0)T,x=(x2,x,1)T,r=aT*x+ε;
根据贝叶斯原理,后验分布、先验分布、似然函数间满足如下关系:
a的后验分布高斯函数为两高斯分布P(r|a)和P(a)乘积,因此均值计算公式为:
方差计算公式为:
如图2所示,本发明的测试方法的检测设备包括上位机1、函数发生装置2、可编程电子负载3、电流电压采集设备4、被测燃料电池5、电流传感器6和负载7。所述被测燃料电池5与负载7的正负极相连,其中电流传感器6通过被测燃料电池5与负载7的正极电缆。同时被测燃料电池5与可编程电子负载3的正负极相连。电流电压采集设备4电压采集部分正负极与被测燃料电池5正负极相连。信号通过上位机1发送指令到函数发生装置2,驱动可编程电子负载3,对燃料电池进行拉载。电流电压采集设备4采集信号发送给上位机1进行计算,得到高频内阻rohm、高频阻抗弧半径Rhf以及低频阻抗弧半径Rlf。信号通过上位机1发送指令到函数发生装置2,驱动可编程电子负载3,对燃料电池进行拉载。其特征在于拉载方式可以为正弦波、三角波、方波或其叠加波形。
本发明还提供一种随燃料电池衰减变化的电化学交流阻抗谱预测系统,包括数据记录模块:在燃料电池衰减实验中,以T为时间周期,记录待测数据的数值变化;方程建立模块:建立待测数据与衰减测试时间的方程;参数计算验证模块:基于建立的方程计算验证方程的参数;待测数据预估模块:根据计算得到参数的估计值,预估待测数据随测试时间的变化;燃料电池电化学交流阻抗谱预估模块:基于预估的待测数据,预估燃料电池电化学交流阻抗谱。
本发明提供的具体实施例如下:
实施例1
在燃料电池衰减实验中,以50h为时间周期T,记录高频内阻数据rohm随衰减时间t的变化。
取二次多项式:r=a2x2+a1x+a0,其中r为燃料电池高频内阻,x为衰减测试的时间,a为系数。
考虑到仪器测量噪声,修正上述二次多项式为r=a2x2+a1x+a0+ε,其中ε为噪声均方差。
取0-300h高频内阻实测数据,计算得到r=a2x2+a1x+a0+ε的实际函数。
如图3所示,可见随着测试时间增加,用于参数估计修正的数据增加,估计曲线逐渐逼近实测曲线。
实施例2
在燃料电池衰减实验中,以50h为时间周期T,记录高频阻抗弧半径Rhf随衰减时间t的变化。
取二次多项式:r=a2x2+a1x+a0,其中r为高频阻抗弧半径Rhf,x为衰减测试的时间,a为系数。
考虑到仪器测量噪声,修正上述二次多项式为r=a2x2+a1x+a0+ε,其中ε为噪声均方差。
取0-300h高频阻抗弧半径Rhf实测数据,计算得到r=a2x2+a1x+a0+ε的实际函数。
在衰减测试过程中,如上所述,根据0~600h,0~900h等实测数据预估随后的衰减实验中高频阻抗弧半径Rhf随时间变化曲线。
如图4所示,可见随着测试时间增加,用于参数估计修正的数据增加,估计曲线逐渐逼近实测曲线。
实施例3
在燃料电池衰减实验中,以50h为时间周期T,记录低频阻抗弧半径Rlf随衰减时间t的变化。
取二次多项式:r=a2x2+a1x+a0,其中r为低频阻抗弧半径Rlf,x为衰减测试的时间,a为系数。
考虑到仪器测量噪声,修正上述二次多项式为r=a2x2+a1x+a0+ε,其中ε为噪声均方差。
取0-300h低频阻抗弧半径Rlf实测数据,计算得到r=a2x2+a1x+a0+ε的实际函数。
在衰减测试过程中,如上所述,根据0~600h,0~900h等实测数据预估随后的衰减实验中低频阻抗弧半径Rlf随时间变化曲线。
如图5所示,可见随着测试时间增加,用于参数估计修正的数据增加,估计曲线逐渐逼近实测曲线。
实施例4
如实施例1所述。
结果如图6所示。
实施例5
本方法通过燃料电池衰减过程中累积的高频、低频测试结果,基于贝叶斯原理及等效电路模型预估交流阻抗谱变化。与现有技术相比,本发明精度高、计算结果可靠性高、同时可以降低燃料电池寿命测试所需时间。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (8)
1.一种随燃料电池衰减变化的电化学交流阻抗谱预测方法,其特征在于,包括如下步骤:
数据记录步骤:在燃料电池衰减实验中,以T为时间周期,记录待测数据的数值变化;
方程建立步骤:建立待测数据与衰减测试时间的方程;
参数计算验证步骤:基于建立的方程计算验证方程的参数;
待测数据预估步骤:根据计算得到参数的估计值,预估待测数据随测试时间的变化;
燃料电池电化学交流阻抗谱预估步骤:基于预估的待测数据,预估燃料电池电化学交流阻抗谱;
待测数据与衰减测试时间的方程为:r=a2x2+a1x+a0+ε,其中ε为噪声均方差,a2,a1,a0为多项式系数,仪器噪声满足高斯分布N(0,ε2),取a=(a2,a1,a0)T,x=(x2,x,1)T,则r=aT*x+ε。
2.根据权利要求1所述的随燃料电池衰减变化的电化学交流阻抗谱预测方法,其特征在于,还包括偏差预估步骤:基于高斯分布的累积分布函数,计算预估值在设定的置信程度下的取值区间,得到概率函数。
3.根据权利要求1所述的随燃料电池衰减变化的电化学交流阻抗谱预测方法,其特征在于,所述待测数据包括高频内阻rohm、高频阻抗弧半径Rhf以及低频阻抗弧半径Rlf。
6.根据权利要求1所述的随燃料电池衰减变化的电化学交流阻抗谱预测方法,其特征在于,时间周期T为1-100h。
7.根据权利要求1所述的随燃料电池衰减变化的电化学交流阻抗谱预测方法,其特征在于,燃料电池电化学交流阻抗谱预估步骤包括:基于得到预估的待测数据,根据燃料电池的Randel电路,模拟相应的奈奎斯特图。
8.一种随燃料电池衰减变化的电化学交流阻抗谱预测系统,其特征在于,采用权利要求1至7中任一项所述的随燃料电池衰减变化的电化学交流阻抗谱预测方法,包括如下模块:
数据记录模块:在燃料电池衰减实验中,以T为时间周期,记录待测数据的数值变化;
方程建立模块:建立待测数据与衰减测试时间的方程;
参数计算验证模块:基于建立的方程计算验证方程的参数;
待测数据预估模块:根据计算得到参数的估计值,预估待测数据随测试时间的变化;
燃料电池电化学交流阻抗谱预估模块:基于预估的待测数据,预估燃料电池电化学交流阻抗谱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110680527.3A CN113343471B (zh) | 2021-06-18 | 2021-06-18 | 随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110680527.3A CN113343471B (zh) | 2021-06-18 | 2021-06-18 | 随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113343471A CN113343471A (zh) | 2021-09-03 |
CN113343471B true CN113343471B (zh) | 2022-11-15 |
Family
ID=77477375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110680527.3A Active CN113343471B (zh) | 2021-06-18 | 2021-06-18 | 随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113343471B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033752A (zh) * | 2011-09-30 | 2013-04-10 | 吴昌旭 | 电动车电池寿命预测方法以及延长方法 |
CN107037363A (zh) * | 2016-10-28 | 2017-08-11 | 四川普力科技有限公司 | 一种基于状态滤波的电池交流阻抗谱测量方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8125193B2 (en) * | 2004-03-26 | 2012-02-28 | Eaton Power Quality Company | Method of testing an electrochemical device |
CN103076572A (zh) * | 2012-12-27 | 2013-05-01 | 天津力神电池股份有限公司 | 一种根据eis谱快速判断负极涂布工艺对电池循环性能影响的方法 |
US10638947B2 (en) * | 2013-12-16 | 2020-05-05 | Medtronic Minimed, Inc. | Use of electrochemical impedance spectroscopy (EIS) in intelligent diagnostics |
CN112289385B (zh) * | 2020-09-17 | 2022-08-09 | 西南交通大学 | 大功率质子交换膜燃料电池电堆电化学阻抗谱预测方法 |
CN112946489B (zh) * | 2021-01-20 | 2023-01-17 | 北京交通大学 | 一种基于低频eis的快速容量评估方法 |
-
2021
- 2021-06-18 CN CN202110680527.3A patent/CN113343471B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033752A (zh) * | 2011-09-30 | 2013-04-10 | 吴昌旭 | 电动车电池寿命预测方法以及延长方法 |
CN107037363A (zh) * | 2016-10-28 | 2017-08-11 | 四川普力科技有限公司 | 一种基于状态滤波的电池交流阻抗谱测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113343471A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhu et al. | A state of charge estimation method for lithium-ion batteries based on fractional order adaptive extended kalman filter | |
CN108254696B (zh) | 电池的健康状态评估方法及系统 | |
Liu et al. | State-of-health estimation of lithium-ion batteries based on electrochemical impedance spectroscopy: a review | |
US7583059B2 (en) | Apparatus and method for estimating state of charge of battery using neural network | |
US11422194B2 (en) | Battery diagnosis apparatus and battery diagnosis method based on current pulse method | |
CN111707951A (zh) | 一种电池组一致性评估方法及系统 | |
JP7157909B2 (ja) | 電池容量の推定方法、および電池容量の推定装置 | |
CN111965559B (zh) | 一种锂离子电池soh在线估计方法 | |
US11567140B2 (en) | Method and device for the diagnosis of battery cells | |
US20230122362A1 (en) | Apparatus and methods for testing electrochemical systems | |
Lyu et al. | SOH estimation of lithium-ion batteries based on fast time domain impedance spectroscopy | |
JP2011519118A (ja) | 電気化学システムの測定できない特性を推定する方法 | |
CN106950507A (zh) | 一种智能时钟电池用高可靠性寿命评估方法 | |
CN115616428A (zh) | 一种“充—检”一体电动汽车电池状态检测与评价方法 | |
CN114062993A (zh) | 一种基于时间卷积网络的cvt误差状态预测方法 | |
CN112989569A (zh) | 一种带温度约束条件的锂电池传感器故障滤波诊断方法 | |
Ranieri et al. | Electronic module for the thermal monitoring of a Li-ion battery cell through the electrochemical impedance estimation | |
CN113343471B (zh) | 随燃料电池衰减变化的电化学交流阻抗谱预测方法及系统 | |
Pérez et al. | Modelling of Li-ion batteries dynamics using impedance spectroscopy and pulse fitting: EVs application | |
CN107918704A (zh) | 电荷放大器贮存寿命预测方法、装置、存储介质和计算机设备 | |
CN116774051A (zh) | 考虑时频域多维数据特征的电池容量快速估计方法 | |
KR102054050B1 (ko) | 배터리의 충전 상태를 추정하는 방법 및 배터리 관리 장치 | |
CN116449219A (zh) | 一种改进vffrls的锂电池模型参数辨识方法 | |
CN113777497B (zh) | 一种退化电池在线soc、soh联合估计方法、装置、存储介质和电子设备 | |
CN113933725B (zh) | 一种基于数据驱动确定动力电池荷电状态的方法 |
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 |