基于工况测试和简化阻抗谱等效电路模型的电池分选方法技术

基于工况测试和简化阻抗谱等效电路模型的电池分选方法，涉及一种电池分选方法。它是为了提高电池分选的速度和可靠性。本发明专利技术的分选步骤：A)对多个待分选的电池同时进行充电或放电操作；B)实时检测多个电池在充电或放电操作结束前特定时间段具有的起始时刻电压值和结束时刻电压值，从而分别获得多个电池在充电或放电操作结束前特定时间段起始时刻电压值和结束时刻电压值之间差值；C)以电池在充电或放电操作结束前特定时间段具有的起始时刻电压值和结束时刻电压值之间差值为标准，对多个电池进行分选配组，形成多个电池组。本发明专利技术适用于电池分选。

【技术实现步骤摘要】
基于工况测试和简化阻抗谱等效电路模型的电池分选方法
本专利技术涉及一种电池分选方法。
技术介绍
目前，对于电池分选技术主要有如下公开文献：公开号为CN102343334的、名称为《一种动力电池的动态分选方法及其系统》；公开号为CN102886352的、名称为《一种磷酸铁锂子电池配组分选方法》；公开号为CN103560277的、名称为《一种电动汽车退役电池重组分选方法》；公开号为CN103769374的、名称为《锂离子电池的分选方法》；公开号为CN103817091的、名称为《电池分选方法及系统》；但上述这些电池分选方法普遍存在电池的分选速度慢、可靠性低、分选参数少的问题。
技术实现思路
本专利技术是为了提高电池分选的速度和可靠性，从而提供一种基于工况测试和简化阻抗谱等效电路模型的电池分选方法。基于工况测试和简化阻抗谱等效电路模型的电池分选方法，该方法的测试工况为：分选测试的环境为常温、标准大气压；步骤A1、静置1小时，以保证待测试的电池处于近似的电化学平衡态；步骤A2、将电池以1/3C恒流充电，至上限截止电压之后，充电模式转为恒压模式；步骤A3、当电池的充电电流减小到1/20C时，认为电池已充满，停止充电；步骤A4、静置1小时，以保证待测试的电池恢复至近似的电化学平衡态；步骤A5、按照动态应力测试工况进行电池充放电测试；其中，工况的最大放电倍率为2C，动态工况测试的结束条件为电池的端电压≤电池的放电截止电压；步骤A6、静置10min，以消除或减弱电池动态工况测试过程中的部分极化过程；步骤A7、将电池以1/3C恒流放电，至电池端电压≤电池的放电截止电压为止；在该测试工况下，电池分选采用下述方法实现：步骤一、电池特性参数辨识步骤，具体为：步骤B1、利用简化的电化学阻抗谱等效电路模型，该等效电路模型包括开路电压OCVe、欧姆内阻Ro和韦伯阻抗ZW；步骤B2、根据步骤B1获得的等效电路模型建立分数阶卡尔曼滤波器所需的状态方程与观测方程，具体为：取流经电池的总电流IL在放电时为正值，数据采样周期为1s；其中△r为微分算子，r为微分阶数，当r为小数时，△r表示分数阶微分算子，当r为整数时，△r为整数微分算子；取分数阶元件ZW是两端电压为UW的状态量，有：对于扩散参数XW、开路电压OCVe和欧姆内阻Ro随着电池荷电状态(SOC)的变化是缓慢的，因此：将上述四个方程改写为矩阵形式，获得分数阶联合卡尔曼滤波器的状态方程：取UL为系统的观测量，则有：UL＝OCVe-ILRo-UWIL表示与流经电池的总电流；取：获得分数阶联合卡尔曼滤波器的观测方程：该方程离散化后，有：其中，w，v分别表示系统的状态噪声和观测噪声；根据分数阶微分的级数定义(又称为Grünwald-Letnikov分数阶微分定义)：其中，另取：由上式得到分数阶微分方程的离散化递推表达形式：定义：为了减少计算量，做以下修改：步骤B3、利用步骤B2获得的分数阶卡尔曼滤波器所需的状态方程与观测方程，对状态、参数及协方差矩阵进行时间更新和测量更新：具体为：初始化：其中，E[x]表示x的数学期望，在方法计算时为经验预设值，表示x在初始时刻(k＝0)的估计值，表示x在初始时刻(k＝0)的噪声协方差的估计值。状态、参数及协方差矩阵的时间更新：其中，Qk是噪声wk的协方差，为k时刻状态和模型参数xk的预测值，为k-1时刻状态和模型参数xk-1的修正值，为k时刻x的噪声协方差矩阵Pk的预测值，为k-1时刻x的噪声协方差矩阵Pk-1的修正值。状态、参数及协方差矩阵的测量更新：其中，Rk是噪声vk的协方差，Lk是k时刻卡尔曼滤波器增益大小。电池的容量值通过电流与时间的积分求和得到：其中，开始计算容量值的时刻T1为动态应力测试工况的起点，结束计算容量值的时刻T2为恒流放电的终点；由此，可计算得到T1～T2时段内随SOC缓慢变化的电池模型参数OCVe，Ro，XW以及电池容量Ccap。至此，电池的四个分选参数已通过试验和上述方法计算得到。这四个分选参数分别是：Ccap、OCVe，Ro，XW。步骤二、对电池进行分选，具体为：步骤B1、对辨识得到的电池模型参数数据进行降采样，由于在电池动态工况测试过程中，数据记录间隔为1s，测试时间约为4小时，采用等间隔降采样的方式将辨识得到的电池模型参数值点数减小到N点(200个数据点左右)，并适当剔除数据序列中首尾处得不太稳定的参数，所述电池模型参数包括开路电压OCVe、欧姆内阻Ro、扩散参数XW；步骤B2、电池参数一致性的初级筛选：采用粗大误差剔除理论中的莱以特准则对各个电池分选参数中偏离均值较大的电池在初级筛选中剔除，以保证剩余电池具有较高的一致性，具体步骤如下：步骤B21、计算各个电池分选参数的均值；步骤B22、计算各个参数与均值之间的平均距离，得到各个电池模型参数的一致性评价参数。例如，计算电池编号为1、分选参数Ro对应的一致性评价参数计算过程如下：计算电池分选参数Ro的各个采样点的均值其中k代表第k个采样数据点，为1～N，n为电池数量。计算电池编号为1的与均值的平均距离：即为电池编号1、欧姆内阻Ro的一致性评价参数。对于每个待分选电池，一致性评价参数有4个：Ccap、docve、dRo、dXw。步骤B23、计算电池容量Ccap的均值和均方差σC，计算docve、dRo、dXw的均方差步骤B24、按照莱以特准则，剔除一致性评价参数超出3σ要求的电池；步骤B25、将剩余电池作为分选对象，重复步骤B21至步骤B24，直至所有电池达到莱以特准则数据一致性的要求；步骤B26、电池参数一致性初级筛选结束。步骤B3、计算各个电池的综合离散度，具体方法如下：步骤B31、将初级筛选后各个电池的一致性评价参数归一化处理，设4个参数对电池一致性的影响程度是相同的，按照等权重计算电池的综合离散度；步骤B32、按照综合离散度将电池排序，得到电池的一致性排序；步骤B33、将离散度较大的电池从待选电池中剔除；完成电池分选。本专利技术带来的有益效果是：1、该方法的特点是对电池进行一个较短时间的工况试验，记录电流、电压数据，即可根据这些数据计算得到电池近似的恒流放电容量、开路电压、欧姆内阻、扩散极化参数。恒流放电容量反映电池在标准充放电机制下，电池电极可嵌入/脱出的离子量，开路电压反映电池的热力学特性(正、负极电势特性)和正负极离子初始荷电状态(初始SOC)情况，欧姆内阻和扩散极化参数反映电池的动力学特性。四个参数从不同角度反映了电池的特性，它们可以较全面的评价电池之间的一致性。该方法的优点是试验时间短，参数意义明确；2、开路电压、欧姆内阻、扩散极化参数的获取采用分数阶参数在线估计方法，快速有效；3、具体的分选过程包括参数初筛和综合离散度评价两步，初筛采用粗大误差处理理论中的莱以特准则，将待分选电池参数中的偏离其它电池参数较大的电池找出并剔除。电池初筛的目的是避免一致性过差的电池影响电池的分选结果。而经过初筛后的电池可以认为特性基本都是相似的，再通过计算综合离散度并将其排序，可以进一步找出一致性更好的电池编号，该分选方法快速、可靠。附图说明图1是简化的电池阻抗谱等效电路模型示意图；图2是一节电池的分选测试工况试验电压数据图；图3是电池参数及数据处理后电池容量测试结果图；图4是电本文档来自技高网...

【技术保护点】
基于工况测试和简化阻抗谱等效电路模型的电池分选方法，其特征是：该方法的测试工况为：分选测试的环境为常温、标准大气压；步骤A1、静置1小时，以保证待测试的电池处于近似的电化学平衡态；步骤A2、将电池以1/3C恒流充电，至上限截止电压之后，充电模式转为恒压模式；步骤A3、当电池的充电电流减小到1/20C时，认为电池已充满，停止充电；步骤A4、静置1小时，以保证待测试的电池恢复至近似的电化学平衡态；步骤A5、按照动态应力测试工况进行电池充放电测试；其中，工况的最大放电倍率为2C，动态工况测试的结束条件为电池的端电压≤电池的放电截止电压；步骤A6、静置10min，以消除或减弱电池动态工况测试过程中的部分极化过程；步骤A7、将电池以1/3C恒流放电，至电池端电压≤电池的放电截止电压为止；在该测试工况下，电池分选采用下述方法实现：步骤一、电池特性参数辨识步骤，具体为：步骤B1、利用简化的电化学阻抗谱等效电路模型，该等效电路模型包括开路电压OCVe、欧姆内阻Ro和韦伯阻抗ZW；步骤B2、根据步骤B1获得的等效电路模型建立分数阶卡尔曼滤波器所需的状态方程与观测方程，具体为：取流经电池的总电流IL在放电时为正值，数据采样周期为1s；Δr=drdtr,r>0]]>其中Δr为微分算子，r为微分阶数，当r为小数时，Δr表示分数阶微分算子，当r为整数时，Δr为整数微分算子；取分数阶元件ZW是两端电压为UW的状态量，有：Δ0.5UW=1WIL=XWIL]]>对于扩散参数XW、开路电压OCVe和欧姆内阻Ro随着电池荷电状态(SOC)的变化是缓慢的，因此：Δ1XW≈0Δ1OCVe≈0Δ1Ro≈0]]>将上述四个方程改写为矩阵形式，获得分数阶联合卡尔曼滤波器的状态方程：Δ0.5111UWXWOCVeRo=0IL00000000000000UWXWOCVeRo]]>取UL为系统的观测量，则有：UL＝OCVe‑ILRo‑UWIL表示与流经电池的总电流；取：x=UWXWOCVeRo,N=0.5111,y=UL]]>获得分数阶联合卡尔曼滤波器的观测方程：ΔNx=0IL00000000000000xy=-101-ILx]]>该方程离散化后，有：ΔNxk=0IL,k-100000000000000xk-1+wyk=-101-IL,kxk+v]]>其中，w，v分别表示系统的状态噪声和观测噪声；根据分数阶微分的级数定义(又称为Grünwald‑Letnikov分数阶微分定义)：ΔNxk+1=Σj=0k(-1)jNjxk-j]]>其中，Nj=diag0.5j1j1j1j,]]>rj=1forj=0r(r-1)...(r-j+1)/j!forj>0,]]>另取：γj=Nj,]]>由上式得到分数阶微分方程的离散化递推表达形式：定义：Ak-1=∂f(xk-1,IL,k-1)∂xk-1|xk-1=x^k-1+=0IL,k-100000000000000,]]>Ck=∂g(xk,IL,k)∂xk|xk=x^k-=-101-IL,k]]>为了减少计算量，做以下修改：Σj=1k(-1)jγjxk+1-j=Σj=1L(-1)jγjxk+1-j,k≤64,L=kk>64,L=64]]>步骤B3、利用步骤B2获得的分数阶卡尔曼滤波器所需的状态方程与观测方程，对状态、参数及协方差矩阵进行时间更新和测量更新：具体为：初始化：x^0=E[x],P0+=E[(x-x^0)(x-x^0)T]]]>其中，E[x]表示x的数学期望，在方法计算时为经验预设值，表示x在初始时刻(k＝0)的估计值，表示x在初始时刻(k＝0)的噪声协方差的估计值；状态、参数及协方差矩阵的时间更新：x^k-=f(x^k-1+,IL,k-1)]]>Pk-=(Ak-1+γ1)Pk-1+(Ak-1+γ1)T+Q+Σj=2LγjPk-j+γjT]]>其中，Qk是噪声wk的协方差，为k时刻状态和模型参数xk的预测值，为k‑1时刻状态和模型参数xk‑1的修正值，为k时刻x的噪声协方差矩阵Pk的预测值，为k‑1时刻x的噪声协方差矩阵Pk‑1的修正值；状态、参数及协方差矩阵的测量更新：Lk=Pk-(Ck)T[CkPk-(Ck)T+Rk]-1]]>x^k+=x^k-+Lkx[yx-g(x^k-,IL,k)]]]>Pk+=(I-LkCk)Pk-]]>其中，Rk是噪声vk的协方差，Lk是k时刻卡尔曼滤波器增益大小；电池的容量值通...

【技术特征摘要】
1.基于工况测试和简化阻抗谱等效电路模型的电池分选方法，其特征是：该方法的测试工况为：分选测试的环境为常温、标准大气压；步骤A1、静置1小时，以保证待测试的电池处于近似的电化学平衡态；步骤A2、将电池以1/3C恒流充电，至上限截止电压之后，充电模式转为恒压模式；步骤A3、当电池的充电电流减小到1/20C时，认为电池已充满，停止充电；步骤A4、静置1小时，以保证待测试的电池恢复至近似的电化学平衡态；步骤A5、按照动态应力测试工况进行电池充放电测试；其中，工况的最大放电倍率为2C，动态工况测试的结束条件为电池的端电压≤电池的放电截止电压；步骤A6、静置10min，以消除或减弱电池动态工况测试过程中的部分极化过程；步骤A7、将电池以1/3C恒流放电，至电池端电压≤电池的放电截止电压为止；在该测试工况下，电池分选采用下述方法实现：步骤一、电池特性参数辨识步骤，具体为：步骤B1、利用简化的电化学阻抗谱等效电路模型，该等效电路模型包括开路电压OCVe、欧姆内阻Ro和韦伯阻抗ZW；步骤B2、根据步骤B1获得的等效电路模型建立分数阶卡尔曼滤波器所需的状态方程与观测方程，具体为：取流经电池的总电流IL在放电时为正值，数据采样周期为1s；其中△r为微分算子，r为微分阶数，当r为小数时，△r表示分数阶微分算子，当r为整数时，△r为整数微分算子；取分数阶元件ZW是两端电压为UW的状态量，有：对于扩散参数XW、开路电压OCVe和欧姆内阻Ro随着电池荷电状态(SOC)的变化是缓慢的，因此：将上述四个方程改写为矩阵形式，获得分数阶联合卡尔曼滤波器的状态方程：取UL为系统的观测量，则有：UL＝OCVe-ILRo-UWIL表示流经电池的总电流；取：获得分数阶联合卡尔曼滤波器的观测方程：该方程离散化后，有：其中，w，v分别表示系统的状态噪声和观测噪声；根据分数阶微分的级数定义(又称为Grünwald-Letnikov分数阶微分定义)：其中，另取：由上式得到分数阶微分方程的离散化递推表达形式：定义：为了减少计算量，做以下修改：步骤B3、利用步骤B2获得的分数阶卡尔曼滤波器所需的状态方程与观测方程，对状态、参数及协方差矩阵进行时间更新和测量更新：具体为：初始化：其中，E[x]表示x的数学期望，在方法计算时为经验预设值，表示x在初始时刻(k＝0)的估计值，表示x在初始时刻(k＝0)的噪声协方差的估计值；状态、参数及协方差矩阵的时间更新：其中，Qk是噪声wk的协方差，为k时刻状态和模型参数xk的预测值，为k-1时刻状态和模型参...

【专利技术属性】
技术研发人员：朱春波，李晓宇，逯仁贵，王婷婷，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23