一种快速得到电池包中所有单体电池容量的方法
技术领域
本发明涉及一种快速得到电池包中所有单体电池容量的方法。
背景技术
本发明与电池组中所有单体电池的容量有关,特别是由多个单体电池相串联所构成的电池组中的所有单体电池的容量。锂离子电池因其具有能量高、电池电压高、工作温度范围宽、贮存寿命长等优点,已经在电动汽车领域、电化学储能领域、3C电子产品等领域得到广泛的应用。电池的有效容量与电池持续工作时间有关,电池阻值与电池瞬间充放电能力密切相关,另外在串联成组的电池包中,单体电池的容量分布与电池包的一致性紧密相关。
锂离子电池容量是衡量电池性能的重要性能指标之一,它表示在一定条件下(放电率、温度、终止电压等)电池放出的电量(可用JS-150D做放电测试),即电池的容量,通常以安培·小时为单位(简称,以A·H表示,1A·h=3600C)。实际使用过程中的锂离子电池,基本不会存在电池满充满放的情况,难以知道电池的总容量,从而影响电池和核电状态(SOC)的估计,对电池容量的精准实时计算,对于电池安全管理存在重要意义。
中国发明专利(专利号:CN109342955A,专利名称:一种锂离子电池容量的推算方法及***),该专利通过将锂离子电池在n个实验温度下进行放电测试,得出该锂电池在对于实验温度下的放电容量,结合阿伦尼乌斯公式建立模型,推算出任意温度点下的电池容量值。该方法需要提前进行n个放电测试,操作复杂,且该方法只能推算出该电池在不同温度下对于的容量值,不能对电池老化后的容量值进行精准计算,实用性较差。
中国发明专利(专利号:CN 108152747 A,专利名称:蓄电池容量的检测方法与装置),该专利通过在获取风力发电机组在利用蓄电池供电变桨过程中的蓄电池放电参数的测量值;根据所述测量值以及通过蓄电池试验利用所述放电参数构建的函数关系,检测所述蓄电池当前容量。该方法需要在使用前进行大量的蓄电池试验构造出测量值与电池容量的关系数据库,在实际使用过程中通过抓取合适的放电参数,计算出需要的测量值,通过查取提前准备好的数据库中对应的容量,评估出蓄电池当前容量。该方法需要提前建立关系数据库,实际运用麻烦,成本昂贵;而且由于电池的非线性特性,在实施过程中,蓄电池容量与对应测量值的对应关系是否还成立,结果的准确度缺乏数据支持。
中国发明专利(专利号:CN108732508,专利名称:一种锂离子电池容量的实时估计方法),该方法通过电池耐久性试验,获取标准容量测试中的充电数据——实时容量VS电压曲线,然后对曲线进行数据拟合在进行差分。差分结果曲线中的峰值位置与大小与电池的剩余容量具有相关性,对比初始情况和某时刻情况下的差分曲线,即可以对于电池容量进行估计。该方法需要对电池进行耐久性试验,这种方法的实际操作性差;差分结果曲线的峰值大小与数据筛选处理方法有关,误差较大。
中国专利(专利号:CN109031153,专利名称:一种锂离子电池的健康状态在线估计方法),该方法通过短期的循环寿命测试,作为模型训练的初始参数,利用容量增量分析法,从容量增量曲线中提取多个特征参数,形成特征参数集合,对电池的健康状态进行刻画,并将特征参数的值作为多输出高斯过程回归模型方法的模型输出,完成对电池容量的在线评估。这种方法需要进行循环寿命测试,实际运用存在难度。
发明内容
本发明的目的是提供一种基于电池正常工作中的充电曲线,处理后提取特征值参数,带入模型进行计算,完成实时计算电池包中所有单体电池的容量值。
本发明的目的可通过下列技术方案来实现:
一种快速得到电池包中所有单体电池容量的方法,其包含如下步骤:
S1、针对目标锂离子电池,获取电池基准曲线数据、电池标称容量Capinitial;
S2、对锂离子电池基准曲线数据进行处理,记录特征值位置电池的容量或充电容量;
S3、计算电池特征值之间容量与总容量的关系系数;
S4、对电池包充电曲线进行处理,记录特征值位置电池的充电容量;
S5、计算S4步骤得到的电池特征值之间容量,并根据S3步骤中所述的关系系数,逐一计算出电池包中所有单体电池的容量。
其中,所述步骤S1中所述的电池基准曲线数据可以是从厂家获取电池SOC-OCV曲线数据,也可以是历史充电曲线数据。
其中,所述步骤S2中所述对锂离子电池基准曲线数据进行处理,记录特征值位置电池的容量或充电容量:根据电池材料类型,对于磷酸铁锂电池取SOC大于20%的数据,求取出容量增量曲线,计算公式为提取出特征值,特征值位置是容量增量曲线中的极大值位置,并记录特征值位置对应的容量Q或充电容量Q′;对于三元材料电池,取SOC大于20%的数据求取d2Q/dV2曲线,计算公式为提取特征值,特征值位置是的位置,并记录特征值位置对应的容量Q或充电容量Q′;Q=SOC*Capinitial;
其中:Q为电池的容量,Q′为充电容量,dQ为容量的微分,d2Q为容量的二阶微分,ΔQk为相邻采样点间容量的差值,V为电池的电压,dV为电压的微分,dV2为电压的二阶微分,ΔVk为相邻采样点间电压的差值,对于每个采样点k,ΔQk=Qk-Qk-1,ΔVk=Vk-Vk-1,ΔQk-1=Qk-1-Qk-2,ΔVk-1=Vk-1-Vk-2。
其中,所述电池基准曲线的特征值有2个,即第一特征值和第二特征值,记录第1个特征值位置的容量Q1或充电容量Q‘1,第2个特征值位置的容量Q2或充电容量Q‘2。
其中,所述步骤S3中计算电池特征值之间容量与总容量的关系系数g=ΔQ/Capinitial;其中ΔQ为第1个特征值到第2个特征值之间的容量差,ΔQ=|Q2-Q1|;如果采用历史充电曲线数据,则ΔQ=|Q‘2-Q‘1|。
其中,所述步骤S4对电池包充电曲线进行处理,记录特征值位置电池充电容量具体为:提取在电池包充电过程中满足ΔV≥X条件的数据。根据电池材料类型,对于磷酸铁锂电池取SOC大于20%的数据利用容量增量分析法,通过对容量增量曲线提取特征值,并记录特征值位置对应的充电容量;对于三元材料电池,取SOC大于20%的数据求取d2Q/dV2曲线,提取特征值,并记录特征值位置对应的充电容量;特征值对应的充电容量包括第1个特征值位置的充电容量Q‘1j和第2个特征值位置的充电容量Q‘2j,其中:Q1j为第1个特征值位置第j号单体电池容量,Q2j为第2个特征值位置第j号单体电池容量,Q‘1j为第1个特征值位置第j号单体充电容量,Q‘2j为第2个特征值位置第j号单体充电容量。
其中,所述S5步骤计算S4步骤得到的电池特征值之间容量,并根据S3步骤中所述的关系系数,逐一计算出电池包中所有单体电池的容量具体为:计算第j号单体电池第1个特征值到第2个特征值之间的容量差,ΔQj=|Q2j-Q1j|=|Q‘2j-Q‘1j|,再根据公式Capj=ΔQj/g得出第j号单体电池容量,然后逐一计算出电池包中其他单体电池容量;其中,Capj为第j号单体电池容量,g为权利要求5中所得到的电池特征值之间容量与总容量的关系系数。
其中,所述X的取值范围为1mV≤X≤10mV。
其中,所述电池包充电曲线可以是最近一次的充电曲线,也可以是最近10次的充电曲线中抽取充电曲线跟最近那次进行组合,使其成为一根包含2个个特征值点曲线。
其中,所述最近10次的充电曲线中抽取充电曲线跟最近那次进行组合,首先选取最近10次充电曲线中跟最次曲线充电电流一样以及环境温度相差小于5度的曲线,在剩余曲线中再选取跟最近那次曲线可以组合成包含两个特征值位置曲线的曲线,如果有多个曲线满足要求,则优先选取环境温度最接近的;如果环境温度一样,则优先选取最新的充电曲线,即时间上最接近最后那次充电的曲线。
本发明的有益效果在于:1.本发明方法在电池组正常充电过程中适用,不影响电池工作输入输出;2.只需要提前获取到电池SOC-OCV曲线数据或者一条历史充电曲线、电池标称容量,不需要额外测试电池参数;3.电池组中所有单体电池的总容量都可实时得出;4.不需要电池进行满充满放即可得出所有单体电池的总容量。
附图说明
图1是本发明的流程图;
图2是磷酸铁锂锂离子电池基准曲线的容量增量曲线;
图3是磷酸铁锂锂离子电池充电数据的容量增量曲线;
图4是三元锂离子电池基准曲线的容量增量曲线和d2Q/dV2曲线;
图5是三元锂离子电池充电数据的容量增量曲线和d2Q/dV2曲线;
图6是电池包中所有单体电池容量直方分布图;
图7是本发明计算出的容量与真实容量的误差图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
如图1所示,一种快速得到电池包中所有单体电池容量的方法,电池包的中电池可以是磷酸铁锂电池或三元材料电池;电池包可以是由多个电芯先并联后再串联的电池包***。针对目标锂离子电池,获取电池基准曲线数据、电池标称容量Capinitial。该锂离子电池基准曲线数据可以是从厂家获取电池SOC-OCV曲线数据,也可以是历史充电曲线数据。从厂家获取电池SOC-OCV曲线数据,SOC取值密度在0.1%~5%之间,包含0.1%和5%。对曲线数据进行处理,根据电池材料类型,对于磷酸铁锂电池取SOC大于20%的数据利用容量增量分析法,通过对容量增量曲线提取特征值,并记录特征值位置对应的容量Q或充电容量Q′;对于三元材料电池,取SOC大于20%的数据求取d2Q/dV2曲线,提取特征值,并记录特征值位置对应的容量Q或充电容量Q′。Q=SOC*Capinitial。
对于磷酸铁锂电池的基准曲线数据进行处理,求取出以SOC为横坐标的容量增量曲线,计算公式为提取出特征值,特征值位置是容量增量曲线中的极大值位置;对于三元电池的基准曲线数据进行处理,求取以电压为横坐标的d2Q/dV2曲线,计算公式为提取出特征值,特征值位置是 的位置。对于磷酸铁锂或者三元电池,在SOC大于20%,小于100%的区间,其曲线的特征值都为2个。记录第1个特征值位置的容量Q1或充电容量Q‘1,第2个特征值位置的容量Q2或充电容量Q‘2。其中:Q为电池的容量,Q′为充电容量,dQ为容量的微分,d2Q为容量的二阶微分,ΔQk为相邻采样点间容量的差值,V为电池的电压,dV为电压的微分,dV2为电压的二阶微分,ΔVk为相邻采样点间电压的差值,对于每个采样点k,ΔQk=Qk-Qk-1,ΔVk=Vk-Vk-1,ΔQk-1=Qk-1-Qk-2,ΔVk-1=Vk-1-Vk-2;
计算电池特征值之间容量与总容量的关系系数,g=ΔQ/Capinitial;其中ΔQ为第1个特征值到第2个特征值之间的容量差,ΔQ=|Q2-Q1|;如果采用历史充电曲线数据,则ΔQ=|Q‘2-Q‘1|。
提取在电池包充电过程中满足ΔV≥X条件的数据。数据包括所有单体电池的电压值,电池包的充电时间和电池包的充电容量值等。X的取值范围为1mV≤X≤10mV。根据电池材料类型,对于磷酸铁锂电池取SOC大于20%的数据利用容量增量分析法,通过对容量增量曲线提取特征值,并记录特征值位置对应的充电容量;对于三元材料电池,取SOC大于20%的数据求取d2Q/dV2曲线,提取特征值,并记录特征值位置对应的充电容量。记录下所有电池曲线的第1个特征值位置的充电容量Q‘1j,记录下所有电池曲线的第2个特征值位置的充电容量Q‘2j,其中:Q1j为第1个特征值位置第j号单体电池容量,Q2j为第2个特征值位置第j号单体电池容量,Q‘1j为第1个特征值位置第j号单体充电容量,Q‘2j为第2个特征值位置第j号单体充电容量。
计算第j号单体电池第1个特征值到第2个特征值之间的容量差,ΔQj=|Q2j-Q1j|;ΔQj也等于第1个特征值到第2个特征值充电容量。如果某一次的充电曲线同时包含了第1特征值和第2特征值,则ΔQj=|Q‘2j-Q‘1j|。如果单次充电曲线不包含两个特征值位置,则可以从最近10次的充电曲线中抽取充电曲线跟最近那次进行组合,使其成为一根包含2个个特征值点曲线。进行组合的两根曲线应该满足,充电电流一样,充电环境温度相差小于5度的要求,否则不进行组合。
计算第j号单体电池对应的总容量,计算公式为:
Capj=ΔQj/g;
由此逐一计算出电池包里每个单体电池的容量。
实例1:
目标锂离子电池是CATL的磷酸铁锂电池,标称容量为Capinitial=180Ah,电池基准曲线数据使用的是SOC-OCV关系曲线数据,数据间隔是ΔSOC=2%。每一个OCV点对应的容量值通过公式:Q=SOC*180计算得出;
采取五点三次平滑滤波法(通过选取要进行平滑滤波位置前后各2个数据,总共5个数据,采用3阶多项式进行拟合,求取平滑滤波后的数值)求取SOC-OCV关系数据的容量增量曲线,得到的容量增量曲线如图2中dQ/dV点画线所示,图2是以SOC为横坐标,左纵坐标是电池开路电压OCV,右纵坐标是dQ/dV,第1个特征值位置为图2中的A点位置,对应的SOC是55.6%,容量Q1=100.08Ah,图中第2个特征值位置为图2中的B点位置,对应SOC是84.9%,容量Q2=152.82Ah;
计算电池特征值之间容量与总容量的关系系数:
g=ΔQ/Capinitial=|Q2-Q1|/Capinitial=(152.82-100.08)/180=0.293
目标锂离子电池组装的电池包,由240个单体电池串联组成的商业化储能***,图3的充电数据条件是电池包在以0.35C,63A电流进行恒流充电至任一个单体电池电压达到3.6V,即停止充电,充电前期的电池包SOC为20%。
充电数据的提取条件是ΔV≥2mV,提取每一个单体电池的电压值和充电电量值,。对电池的电压值和充电电量值采用五点三次平滑滤波方法求取容量增量曲线,电池包中第1号单体电池的容量增量曲线如图3所示,图3是以SOC为横坐标,左纵坐标是电池电压,右纵坐标是dQ/dV,图中第1个特征值位置为图3中的C点位置,对应充电容量值为Q‘11=67.66Ah,图中第2个特征值位置为图3中的D点位置,对应的充电容量值为Q‘21=119.2Ah;
计算电池包1号单体电池的总容量:
Cap1=ΔQ1/g=|Q‘21-Q‘11|/g=(119.2-67.06)/0.293=177.9Ah
同样的计算方法计算出电池包中2号单体电池的总容量到240号单体电池的总容量,并把所有的单体电池总容量并绘制成直方分布图,图6是电池包中所有单体电池容量直方分布图。
并把采用该发明方法计算出来的容量值与电池真实容量值进行比较,容量误差率如图7,图中显示最大容量误差为3.86%。
实例2:
目标锂离子电池是力神的21700三元电池,标称容量为Capinitial=4.5Ah,电池基准曲线数据使用的是SOC-OCV关系曲线数据,数据间隔是ΔSOC=3%。每一个OCV点对应的容量值通过公式:Q=SOC*4.5计算得出;
采取五点三次平滑滤波方法求取SOC-OCV关系数据的容量增量曲线和d2Q/dV2曲线,得到的容量增量曲线如图4中的dQ/dV点画线所示,d2Q/dV2曲线如图4的虚线所示,图4是以开路电压OCV为横坐标,左纵坐标是电池SOC,右第1个纵坐标是dQ/dV,右第2个纵坐标是d2Q/dV2,图中第1个特征值位置为图4中的E点位置,对应SOC值为39.8%,容量Q1=1.791Ah,图中第2个特征值位置为图4中的F点位置,对应SOC值为76.4%,容量Q2=3.438Ah;
计算电池特征值之间容量与总容量的关系系数:
g=ΔQ/Capinitial=|Q2-Q1|/Capinitial=(3.438-1.791)/4.5=0.336
目标锂离子电池以1C,4.5A电流进行恒流充电至电压达到4.2V,即停止充电,充电前期的电池包SOC为15%。
充电数据的提取条件是ΔV≥3mV,提取电池的电压值和充电电量值,。对电压值和充电电量值采用五点三次平滑滤波方法求取容量增量曲线和d2Q/dV2曲线,该电池的充电数据的容量增量曲线和d2Q/dV2曲线如图5所示,图5是以电压为横坐标,左纵坐标是电池SOC,右第1个纵坐标是dQ/dV,右第2个纵坐标是d2Q/dV2,图中第1个特征值位置为图5中的H点位置,对应充电容量值为Q‘11=0.945Ah,图中第2个特征值位置为图5中的M点位置,对应的充电容量值为Q‘21=2.601Ah;
计算电池包1号单体电池的总容量:
Cap1=ΔQ1/g=|Q‘21-Q‘11|/g=(2.443-0.945)/0.336=4.46Ah
同样的计算方法计算出电池包中其它单体电池的总容量。
实例3:
目标锂离子电池是国轩高科的软包三元电池,15Ah两并后串联,标称容量为Capinitial=30Ah,电池基准曲线数据使用的是历史曲线数据,充电电流大小为6A,数据提取条件是ΔV≥5mV;
采取五点三次平滑滤波方法求得历史充电曲线数据的容量增量曲线和d2Q/dV2曲线,并以充电电压为横坐标,左纵坐标是电池SOC,右第1个纵坐标是dQ/dV,右第2个纵坐标是d2Q/dV2,求得第1个特征值位置对应的充电容量Q‘1=16.53Ah,第2个特征值位置对应的充电容量Q‘2=26.13Ah;
计算电池特征值之间容量与总容量的关系系数:
g=ΔQ/Capinitial=|Q‘2-Q‘1|/Capinitial=(26.13-16.53)/30=0.32
目标锂离子电池组装的电池包,由32个电芯通过2并16串组成电池包,目标锂离子电池以0.5C,15A电流进行恒流充电至任一单体电池电压达到4.2V,即停止充电,充电前期的电池包SOC为13%。
充电数据的提取条件是ΔV≥1mV,提取每一个单体电池的电压值和充电电量值,充电完成后,剔除掉SOC小于15%的数据后,。对每一个单体电池的电压值和充电电量值采用五点三次平滑滤波方法求取容量增量曲线和d2Q/dV2曲线,记录每一个单体电池第1个特征值位置的充电容量值和第2个特征值位置的充电容量值,其中第1号单体电池第1个特征值位置对应充电容量值为Q‘11=0.945Ah,第2个特征值位置对应的充电容量值为Q‘21=2.601Ah;
计算电池包中1号单体电池的总容量:
Cap1=ΔQ1/g=|Q‘21-Q‘11|/g=(21.56-11.576)/0.32=31.2Ah
同样的计算方法计算出电池包中2号单体电池到16号单体电池的总容量。
实施例4:
实施案例4与实施案例3类似,不同的是,在处理充电曲线时,单次充电曲线不包含两个特征值位置,因此从最近10次的充电曲线中抽取充电曲线跟最近那次进行组合,使其成为一根包含2个个特征值点曲线。首先选取最近10次充电曲线中跟最次曲线充电电流一样以及环境温度相差小于5度的曲线,在剩余曲线中再选取跟最近那次曲线可以组合成包含两个特征值位置曲线的曲线,如果有多个曲线满足要求,则优先选取环境温度最接近的;如果环境温度一样,则优先选取最新的充电曲线(即时间上最接近最后那次充电的曲线)。数据提取条件可以是ΔV≥10mV条件的数据。再利用前述公式逐一计算出电池包中单体电池容量。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。