一种动力电池的峰值功率预测方法
【专利摘要】一种动力电池的峰值功率预测方法，涉及动力电池的峰值功率预测技术。本发明是为了提高动力电池的峰值功率预测的准确性。本发明的方法包括两部分：A基于简化的电化学阻抗谱等效电路模型和分数阶联合卡尔曼滤波的参数在线估计、B基于零状态响应和零输入响应分解的电池峰值功率预测方法。本发明选用了包含分数阶元件的简化阻抗谱模型作为电池峰值功率预测的参考模型，因此本发明提出的基于此模型的动力电池峰值功率预测方法不仅可以准确的预测电池的短时峰值输出功率，也可以准确预测电池在较长时间段内的峰值功率输出能力。本发明适用于电动汽车在线峰值功率预测。
【专利说明】一种动力电池的峰值功率预测方法
【技术领域】
[0001] 本发明涉及动力电池的峰值功率预测技术。
【背景技术】
[0002] 国内专利中公开的有关电池峰值功率预测方法的内容较少，国内外的论文期刊中 有一些动力电池峰值功率预测的方法。
[0003] 这些方法大都采用电池的一阶RC等效电路模型，以S0C、端电压、电流作为峰值功 率预测的限制条件，早期的功率预测方法采用离线电池模型参数数据完成对电池输出/回 馈峰值功率的预测，但是电池随着老化或者温度、工况等环境的变化，模型参数也会随之改 变，因此，基于离线模型参数数据的功率估计方法可靠性较差；
[0004] 最新的发表的论文中，电池模型参数数据由参数估计方法在线得到，可以有效提 高估算的准确性，但是由于一阶RC模型自身的误差较大，功率预测只能局限在较短的时间 内是有效的。
[0005] 为了减少状态方程在用于功率预测时的离散化误差，论文中采用较小的离散化时 间间隔（如取时间间隔为Is或者0. Is)，通过递推计算的方法，得到电池在当前时刻之后 l〇s处的峰值功率预测值或者60s处的预测值，这种方法的缺点是存在大量的递推计算，实 时性差。另外这种方法也无法改善由于模型自己误差造成的功率估计不准确的问题。
【发明内容】
[0006] 本发明是为了提高动力电池的峰值功率预测的可靠性，减少计算量，从而提供一 种动力电池的峰值功率预测方法。
[0007] -种动力电池的峰值功率预测方法，它由以下步骤实现：
[0008] 步骤A、电池模型参数在线估计的步骤，具体为：
[0009] 步骤A1、在对二次电池建模时，由于电动汽车运行工况频率特性，因此电池电化学 阻抗谱模型中的中频率的阻抗特性可以由常用的纯阻性元件R和常相位元件Q并联电路简 化为纯阻性元件R来描述，得到简化后的电池电化学阻抗谱等效电路模型；
[0010] 该简化后的电化学阻抗谱等效电路模型包括开路电压〇CVe、欧姆内阻R。和韦伯阻 抗z w;
[0011] 步骤A2、根据步骤A1获得的简化后的电化学阻抗谱等效电路模型建立分数阶卡 尔曼滤波器所需的状态方程与观测方程，具体为：
[0012] 取流经二次电池的总电流L在放电时为正值，数据采样周期为Is ;
[0013]
【权利要求】
1. 一种动力电池的峰值功率预测方法，其特征是：它由以下步骤实现： 步骤A、电池模型参数在线估计的步骤，具体为： 步骤A1、在对二次电池建模时，由于电动汽车运行工况频率特性，因此电池电化学阻抗 谱模型中的中频率的阻抗特性由常用的纯阻性元件R和常相位元件Q并联电路简化为纯阻 性元件R来描述，得到简化后的电池电化学阻抗谱等效电路模型； 该简化后的电化学阻抗谱等效电路模型包括开路电压OCV^欧姆内阻R。和韦伯阻抗Zw ; 步骤A2、根据步骤Al获得的简化后的电化学阻抗谱等效电路模型建立分数阶卡尔曼 滤波器所需的状态方程与观测方程，具体为： 取流经二次电池的总电流k在放电时为正值，数据采样周期为Is ;
其中A 为微分算子，r为微分阶数，当r为小数时，A 表示分数阶微分算子，当！为 整数时，A 为整数微分算子； 取分数阶元件Zw是两端电压为Uw的状态量，有：
对于电池模型参数，扩散参数Xw、开路电压OCVe和欧姆内阻R。随着电池荷电状态（SOC) 的变化是缓慢的，因此：
将上述四个方程改写为矩阵形式，获得分数阶联合卡尔曼滤波器的状态方程：
取队为系统的观测量，则有： Ul = OCV6-IlR0-Uw 込表示与流经电池的总电流； 取：
获得分数阶联合卡尔曼滤波器的观测方程：
该方程离散化后，有：
其中，《，V分别表示系统的状态噪声和观测噪声； 根据分数阶微分的级数定义（又称为Griinwald-Letnikov分数阶微分定义）：
根据分数阶微分的级数定义式，其中(-的计算量将随着时间的增加而 戶 不断增大，这种情况不适合工程应用，为此，将上式改写为下面的形式：
步骤A3、利用步骤A2构建的分数阶卡尔曼滤波器所需的状态方程与观测方程，对状 态、参数及协方差矩阵按照分数阶联合卡尔曼滤波算法进行时间更新和测量更新： 具体为： 初始化：
其中，E[x]表示X的数学期望，在方法计算时为经验预设值，&表示X在初始时刻（k =〇)的估计值，A+表示X在初始时刻（k = 0)的噪声协方差的估计值； 状态、参数及协方差矩阵的时间更新：
其中，Qk是噪声wk的协方差，右为k时刻状态和模型参数X k的预测值，为k-1时刻 状态和模型参数Xlrt的修正值，I为k时刻X的噪声协方差矩阵Pk的预测值,i二为k-1时 刻X的噪声协方差矩阵Plrt的修正值； 状态、参数及协方差矩阵的测量更新：
其中，Rk是噪声vk的协方差，Lk是k时刻卡尔曼滤波器增益大小； 步骤A4、采集电池的端电压队和流经二次电池的总电流L利用步骤Al获得的简化后 的电化学阻抗谱等效电路模型和步骤A3更新后的分数阶卡尔曼滤波器所需的状态方程与 观测方程，得到开路电压OCV^欧姆内阻R。、扩散参数X w的估计值，将获得的开路电压OCV^ 欧姆内阻R。、扩散参数Xw的估计值作为的电池的估计结果，完成基于分数阶联合卡尔曼滤 波的二次电池简化阻抗谱模型参数在线估计； 步骤B、根据步骤A获得的电池模型参数在线估计结果，进行功率预测的步骤： 步骤B1、推导分数阶元件的韦伯阻抗Zw是否为线性时不变元件： 设分数阶元件的两端电压Uw的初始值为0,当施加一个幅值为L的阶跃电流激励时， 分数阶元件的两端电压Uw在Ns后的电压响应情况计算如下，N为正数： 设电池模型的开路电压〇CVe、欧姆内阻R。、扩散参数Xw在功率预测过程中数值大小是 不变的； 则仏在1^ = Is时的电压响应为：
其中，符号'代表预测值； 仏在1^彡2时的电压响应为：
计算得：
计算后得：
由此递推得到IOs和60s处的电压响应为：
由此，当电池模型参数Xw不变或者缓慢变化时，输出值Uw与k为线性关系，推断分数 阶元件Zw为线性时不变元件； 贝1J，分数阶元件的功率预测方法为： 为了预测分数阶元件在k+ATs处的电压响应f/fn)[+Ar，将该电压响应分为零状态响应 ^U\k+AT和零输入响应 .k+AT:
其中零状态响应为：
对于k+10s时刻，a= 3. 524 ;对于k+60s时刻，a=8. 722 ; 零输入响应由k时刻之前的数据决定，取分数阶微分的时间记忆长度L= 60 ;k+1时刻的零输入响应为：
由此递推可得到k+AT时刻分数阶元件处零输入电压响应：
其中，是电池模型分数阶元件2"两端端电压在（k+1-L)?k时刻之间的估计 值，当预测k+10时刻时，b=a，a为一组常系数矩阵； 当预测k+60时刻时，b= @@为另外一组常系数矩阵； 由此，可得到电池恒流放电电流为Imax时的电池在k+AT时刻的端电压预测值：
放电峰值功率预测的方法为： 若电池荷电状态是电池极限工作状态的限制条件，SoCniin是电池放电终止荷电状态的 最小值，则得到此时的最大放电电流值：
Capacity为电池容量值，单位为安培*小时（Ah)，SoCk为k时刻的电池荷电状态，SoCmin 为电池放电限定的最小荷电状态； 若电池端电压队为电池极限工作状态的限制条件，若此时以最大放电电流Iniax对电池 放电，电池模型参数OCVe在k+AT时刻的预测值为：
上式中为通过分数阶联合卡尔曼滤波算法计算得到的电池模型参数OCVe在k 时刻的估计值，为OCVe在k+ A T时刻的预测值；Imax是在端电压作为电池极限工 作状态限制条件时的最大放电电流值，是待求解值，/=x是电池的最大放电电流；OCVk是在 k时刻的电池开路电压值；为电池在k+aT时刻，假定以恒流放电时对应的 电池开路电压值；由于多数情况下，OCV的变化都较小并且缓慢，因此可认为在AT时段内 的开路电压变化量与放电电流成线性关系； 进而推算，假设电池在放电电流为Imax时，电池端电压的推算值Uuk+AT :
其中，端电压的推算值是k时刻开路电压估计值、A T时段内开路电压变化值、欧姆内 阻电压差、分数阶元件Zw的零状态响应电压值、零输入响应电压值之和； 由此推出，当以端电压队作为电池的极限工作状态的约束条件时，电池的最大工作电 流为：
k+AT时刻达到峰值放电功率时，综合考虑上述限制条件，最大放电电流值为：
k+AT时刻的峰值放电功率为：
电池的回馈电流峰值功率的预测方法与上述放电功率的预测方法同理： 设込为正值； 若SoC为电池极限工作状态的限制条件，SoCmax是电池的最大荷电状态值，则得到此时 的最小回馈电流值：
若队为电池极限工作状态的限制条件，则得到此时的最小回馈电流值：
k+AT时刻达到峰值回馈功率时，综合考虑上述限制条件，最小回馈电流值为：
由此，得到电池电流回馈的峰值功率：
完成动力电池的峰值功率预测。
2.根据权利要求1所述的一种动力电池的峰值功率预测方法，其特征在于步骤A中， 由于方法中使用的电池模型基于电池的电化学阻抗谱测试数据，电池模型参数有明确 的物理意义，欧姆内阻R。的物理意义为： R0?Rn+RsEi+Rct 其中，Rn为高频欧姆阻抗，Rsei为SEI膜阻抗，R。,为电荷转移阻抗； 另外，韦伯阻抗由以下公式来定义：
其中，W为离子扩散系数，为了便于阻抗参数在线估计，取：
得到： 〇
【文档编号】G01R21/00GK104267354SQ201410592570
【公开日】2015年1月7日 申请日期:2014年10月29日 优先权日:2014年10月29日
【发明者】朱春波, 李晓宇, 魏国, 裴磊 申请人:哈尔滨工业大学