CN105912799B - 一种液态或半液态金属电池的建模方法 - Google Patents

一种液态或半液态金属电池的建模方法 Download PDF

Info

Publication number
CN105912799B
CN105912799B CN201610268697.XA CN201610268697A CN105912799B CN 105912799 B CN105912799 B CN 105912799B CN 201610268697 A CN201610268697 A CN 201610268697A CN 105912799 B CN105912799 B CN 105912799B
Authority
CN
China
Prior art keywords
battery
soc
equivalent
circuit
capacity
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
Application number
CN201610268697.XA
Other languages
English (en)
Other versions
CN105912799A (zh
Inventor
王康丽
王大磊
蒋凯
程时杰
余丽红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201610268697.XA priority Critical patent/CN105912799B/zh
Publication of CN105912799A publication Critical patent/CN105912799A/zh
Application granted granted Critical
Publication of CN105912799B publication Critical patent/CN105912799B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Abstract

本发明公开了一种液态或半液态金属电池的建模方法,包括以下步骤,根据电池阻抗谱进行拟合,构建电池的阻抗谱拟合电路,根据阻抗谱拟合电路构建电池的等效电路模型;采用对称脉冲对电池进行混合脉冲功率性能测试,获取用于辨识电池的开路电压、电阻、电容参数的测试数据;根据测试数据拟合电池的电动势、欧姆内阻、极化电阻、极化电容、扩散等效电阻和扩散等效电容与SOC的函数关系,辨识出等效电路模型的参数;采用安时法计算电池的SOC;对SOC进行修正,获得修正SOC;根据修正SOC,对等效电路模型的参数进行修正;完成建模;该建模方法能准确模拟液态或半液态金属电池的外特性,为液态或半液态金属电池这类新型电网级储能电池的充放电管理和后续应用奠定基础。

Description

一种液态或半液态金属电池的建模方法
技术领域
[0001] 本发明属于电池储能建模技术领域,更具体地,涉及一种液态或半液态金属电池 的建模方法。
背景技术
[0002] 随着全球能源和环境危机越来越严重,可再生能源成为关注重点,但可再生能源 发电具有间歇性和波动性,将其直接接入电网会降低电能质量,甚至威胁电力系统的安全 性和可靠性。借助大规模储能技术可以有效解决可再生能源间歇性和波动性带来的诸多问 题,实现发电侧与电网侧的隔离,大幅提高电网接纳可再生能源的能力。有效消纳可再生能 源不仅能够削峰填谷,改善供需矛盾,也能削弱发、输、配、用等环节在时间上的统一性,提 高能源利用效率。而液态或半液态金属电池容易放大和生产,储能成本可满足大规模储能 要求,在电网储能规模应用方面具有明显优势。掌握这一类新型电池体系的储能参数,特别 是其输出特性,是非常紧迫和必要的工作。通常通过电池建模来了解不同工况下的电池响 应情况并制定相应的应对措施以有效保证电网供电的安全性。
[0003] 现有的电池模型包括基于阻抗谱分析的内阻法对铅酸电池的建模模型,针对质子 交换膜燃料电池提出的基于电极过程和热力学势的电化学模型,分别对铅酸电池和锂离子 电池建立的等效电路模型,还有将遗传算法与等效电路模型相结合以确定电池最优模型参 数;上述模型主要针对铅酸电池、锂离子电池和燃料电池;液态或半液态金属电池是面向电 网的新型储能电池,其特殊的内部状态以及低电压、大电流的特点使得现有的电池模型对 其并不适用,套用现有的建模方法会造成较大误差。
发明内容
[0004] 针对现有技术的以上缺陷或改进需求,本发明提供了一种液态或半液态金属电池 的建模方法,由此准确模拟液态或半液态金属电池这类新型电池的储能参数和输出特性。
[0005] 为实现上述目的,按照本发明的一个方面,提供了一种液态或半液态金属电池的 建模方法,具体如下:
[0006] (1)通过电化学测试获取液态或半液态金属电池的阻抗谱,并根据获得的阻抗谱 进行拟合,构建所述电池的阻抗谱拟合电路;根据所述阻抗谱拟合电路构建所述电池的等 效电路模型;
[0007]本发明中,根据阻抗谱拟合电路构建得的等效电路模型为二阶Thevenin等效电路 丰旲型;
[0008] 根据该等效电路模型确定建模所需的参数为电池电动势U。。、电池欧姆内阻R〇、并 联的极化电阻Rp和极化电容CP、并联的扩散等效电阻匕和扩散等效电容Cw;
[0009] 由于阻抗谱能够真实反映电池的实际物理构造和化学特性,结合阻抗谱选择的模 型与液态或半液态金属电池实际的物理构造和化学特性的匹配度更高;
[0010] (2)采用对称脉冲对液态或半液态金属电池进行混合脉冲功率性能(HPPC)测试, 获取电池的电流、电压、充电容量、充电时间、放电容量和放电时间的测试数据;这些测试数 据将用于辨识二阶Thevenin等效电路模型的参数;
[0011] (3)根据上述测试数据拟合出液态或半液态金属电池的电动势U。。、欧姆内阻Ro、极 化电阻Rp、极化电容Cp、扩散等效电阻Rw和扩散等效电容。与荷电状态(SOC)的函数关系,辨 识出等效电路模型的参数;
[0012] ⑷采用安时法计算获得液态或半液态金属电池的S0C;
[0013] (5)根据Rakhmatov模型获得修正式;采用修正式对采用安时法获得的SOC进行修 正,获得修正SOC以减小因电池的倍率容量效应和恢复效应产生的容量误差;
[0014] ⑹根据修正S0C,以及在步骤(3)获得的电动势U。。、欧姆内阻Ro、极化电阻心、极化 电容CP、扩散等效电阻1和扩散等效电容CV^SOC的函数关系,对等效电路模型的参数进行 修正。
[0015] 优选地,上述液态或半液态金属电池的建模方法,步骤⑵具体如下:
[0016] (2.1)用额定电流对电池进行N次完整的充放电;N>3;
[0017] (2.2)将电池充满并将电池静置tr时长;
[0018] (2.3)将电池放电0. IT时长并静置tr时长,再对电池施加对称脉冲信号;其中,T为 步骤⑴中每次放电的平均时间;
[0019] 在本步骤中,采用对称脉冲对电池进行HPPC测试,充、放电脉冲幅值相同;
[0020] (2.4)重复步骤(2.3)直至将电池满放;获得用于辨识电池的开路电压、电阻和电 容的测试数据,包括电池电流、电压、充电容量、充电时间、放电容量和放电时间的测试数 据。
[0021] 优选地,上述液态或半液态金属电池的建模方法,其步骤(2.3)中,对称脉冲信号 采用幅值为Ip的短时阶跃信号;11>为0.2C〜0.5C。
[0022] 优选地,上述液态或半液态金属电池的建模方法,其步骤(2.3)中对电池施加对称 脉冲的具体工步为:Ip放电l〇s,静置40s;IP充电10s,静置40s。
[0023] 优选地,上述液态或半液态金属电池的建模方法,步骤⑶具体如下:
[0024] (3.1)采用仿真工具搭建等效电路模型;并在软件工具中对等效电路模型的参数 赋值;
[0025] (3.2)从所述测试数据中获取每个SOC点静置tr时长后的开路电压值,并对电动势 U。。进行赋值;
[0026] (3.3)截取测试数据中脉冲工步获取的时间、电流和电压数据作为仿真工具的输 入;并将电池欧姆内阻Ro、极化电阻Rp、极化电容Cp、扩散等效电阻Rw和扩散等效电容Cw参数 设为仿真工具里的变量;
[0027] (3.4)通过仿真工具进行曲线拟合以辨识参数,在每个施加过脉冲的SOC点获取一 组辨识出的变量值;
[0028] (3.5)采用算式拟合工具对SOC值和与其对应的辨识出的变量值进行算式拟合,获 得等效电路模型的参数与SOC的函数关系;包括电池电动势U。。、欧姆内阻Ro、极化电阻Rp、极 化电容Cp、扩散等效电阻Rw和扩散等效电容(^参数与SOC的函数关系。
[0029] 优选地,上述液态或半液态金属电池的建模方法,其步骤(3.5)中,对开路电压、电 阻、电容参数与SOC的关系采取分段函数拟合;
[0030] 将电池欧姆内阻Ro与SOC的关系分段拟合为:
[0031]
Figure CN105912799BD00061
[0032] 将极化电阻心与30(:的关系分段拟合为:
[0033]
Figure CN105912799BD00062
[0034] 将极化电容心与50(:的关系分段拟合为:
[0035]
Figure CN105912799BD00063
[0036] 将扩散等效电阻IV^SOC的关系分段拟合为:
[0037]
Figure CN105912799BD00064
[0038] 将扩散等效电容CV^SOC的关系分段拟合为:
[0039]
Figure CN105912799BD00065
[0040] 将电池电动势U。。与SOC的关系分段拟合为:
[0041]
Figure CN105912799BD00066
[0042] 其中,f1>2...n (SOC) Wm(SOC) Am(SOC) Jm(SOC)、rHn(SOC)、zHn(SOC) 都是自变量为SOC的函数;
[0043] 其函数式均为其
Figure CN105912799BD00067
中,(30、(31、〇2"_(^、(11、¥1、(12、¥2均为实数拟合系数,在每一个分段内,0彡7彡6,分段数量1彡11 ^Ξ4 ;
[0044]采用分段函数可以提高对电池模型各参数的辨识精度;其中,提高欧姆内阻R0的 辨识精度对液态或半液态金属电池储能特性建模至关重要,液态或半液态金属电池的特性 与锂电、锂硫等固态电池的区别主要体现在欧姆内阻Ro的变化上;提高欧姆内阻Ro的精度可 以提高最终获得的电池模型与实际电池的匹配程度;
[0045] 在本步骤中,仿真工具可采用Matlab/Simulink软件;算式拟合工具可采用Istopt 软件;采用Matlab/Simul ink以及Istopt软件进行分工辨识,相比于现有技术中的采用单一 软件辨识,有利于提高本模型的辨识精度。
[0046] 优选地,上述液态或半液态金属电池的建模方法,步骤⑷具体如下:
[0047] 将充放电电流i (t)作为输入,采用安时法获得
[0048]
Figure CN105912799BD00071
[0049] 其中,to为起始积分时刻,t为当前时刻;i⑴为正时为充电,为负时为放电;Qo为额 定放电容量,SOC⑴为t时刻电池S0C。
[0050] 优选地,上述液态或半液态金属电池的建模方法,步骤⑸具体如下:
[0051] (5.1)通过对Rakhmatov模型进行演绎,获得修正式
[0052]
Figure CN105912799BD00072
[0053] 其中,X为拟合项数,Uk为k时刻之前电池的不可用容量,Ik为k时刻电池电流,ts为 米样时间;81、131、32』2"1上为;1^(1^)的拟合系数;
[0054] 其中,f(ti)是
Figure CN105912799BD00073
.的拟合式:
[0055]
Figure CN105912799BD00074
[0056] β为电池参数,是电极表面活动载流子被补偿的速率,用于衡量电池放电特性;ti 为静置时间;ti为静置前的充、放电时间,m为求和项数;
[0057] (5.2)采用上述修正式对电池SOC进行修正,获得
[0058]
Figure CN105912799BD00075
[0059] 其中,SOC/⑴是指t时刻修正的S0C。
[0060] 对SOC进行修正,可减少因倍率容量效应和恢复效应造成的容量误差。
[0061] Rakhmatov模型为
Figure CN105912799BD00076
步骤(5 · 1)将该模型转 化为更易实现的迭代式,将恒流与变电流的两种情况统一,且去掉了求和符号,大大减少了 计算量,可在很大程度上节省计算处理的时间。
[0062] 优选地,上述液态或半液态金属电池的建模方法,还包括步骤(7),根据电池端电 压对所述等效电路模型与等效电路模型的参数进行验证,具体如下:
[0063] (7.1)根据步骤(6)获得的等效电路模型的参数,获取欧姆内阻Ro两端的欧姆压降 Uo、极化压降Up、扩散等效压降Uw,以及电池端电压Ubat = Ucic-Uo-Up-Uw;
[0064] 其中,UbaA电池端电压,Up为电池极化电阻Rp与极化电容Cp构成的RC环的两端极 化压降,UW为电池扩散电阻1与扩散电容Cw构成的RC环的两端的扩散等效压降;
[0065] (7.2)将电池端电压1]^仿真值与在各工况下测试获得的电池端电压UbaJ|〗试数据 进行比对,通过仿真值与实测值的匹配程度,验证建模获得的所述等效电路模型与等效电 路模型的参数与所述电池的匹配程度;仿真值与实测值越接近,表明建模获得的等效电路 模型与等效电路模型的参数与电池越匹配。
[0066] 总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有 益效果:
[0067] (1)本发明提供的液态或半液态金属电池的建模方法,结合阻抗谱拟合进行等效 电路模型选择,与液态或半液态金属电池的实际物理结构和化学特性相结合;与现有技术 从四种经典模型当中任意选择相比,本发明采用电化学理论测试结果辅助选取合适的等效 电路模型,使模型更加符合电池的实际情况,减少了不必要的模型误差;
[0068] (2)本发明提供的液态或半液态金属电池的建模方法,采用分段函数拟合电池欧 姆内阻Ro、极化电阻心、极化电容Cp、扩散等效电阻Rw和扩散等效电容(^参数,提高了参数辨 识精度;
[0069] 而液态或半液态金属电池与锂电、锂硫等固态电池的区别主要体现在电池欧姆内 阻Ro的变化上;提高电池欧姆内阻Ro的辨识精度,可进一步提高通过本建模方法获得的模型 的精度;
[0070] (3)本发明提供的液态或半液态金属电池的建模方法,在采用安时法获得的SOC基 础上,加入修正式对SOC进行进一步校正,考虑到了电池倍率变化或工步变化对充放电曲线 的影响,减小了实际应用中常见的容量误差;
[0071] 现有的电池建模方法忽略电池的倍率容量效应和恢复效应,对于不同倍率充、放 电以及含有静置工况的充、放电情况,模型曲线误差较大;相较而言,本发明提供的建模方 法,对SOC进行修正后,进一步提高了等效电路模型参数的与电池实际参数的匹配程度;
[0072] (4)本发明提供的液态或半液态金属电池的建模方法,在进行修正时,对现有的 Rakhmatov模型进行了改进,获得更易实现的迭代公式,并将恒流和变电流两种情况统一处 理,具有提高处理的速度和效率的作用;解决了将现有Rakhmatov模型公式直接应用于程序 中,由于计算量大而容易造成程序卡顿的问题;
[0073] (5)本发明提供的液态或半液态金属电池的建模方法,针对液态或半液态金属电 池进行建模,有利于对此类新型电池的管理、测试、生命周期管理、应用条件控制等工作,对 于实现电池在电网储能领域的稳定、可靠工作具有重要意义。
附图说明
[0074] 图1是实施例中液态金属电池阻抗谱拟合获得的曲线图;
[0075] 图2是实施例中液态金属电池的Thevenin等效电路;
[0076] 图3是对实施例中液态金属电池进行HPPC测试获得的电压波形;
[0077] 图4是实施例中进行Simulink拟合获得的效果图(SOC = O.7);
[0078]图5是实施例中Simul ink拟合计算端电压与真实端电压的误差分布图;
[0079] 图6是实施例中获得的SOC-OCV拟合曲线;
[0080] 图7是实施例中液态金属电池建模流程图;
[0081] 图8是实施例中恒流工况下获得的电池端电压仿真值与实测值的对比图;
[0082] 图9是实施例中HPPC工况下获得的电池端电压仿真值与实测值的对比图。
具体实施方式
[0083] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对 本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并 不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要 彼此之间未构成冲突就可以相互组合。
[0084] 本发明针对液态或半液态金属电池提出一种建模方法,旨在准确模拟液态或半液 态金属电池的内部物理结构,真实反映电池外特性;本发明实施例提供的建模方法,具体如 下:
[0085] (1)为液态或半液态金属电池构建合适的等效电路模型;并根据模型确定建模需 要的参数:电池电动势U。。、电池欧姆内阻Ro、并联的极化电阻Rp和极化电容Cp、并联的扩散等 效电阻Rw和扩散等效电容Cw;
[0086] 电池等效电路模型的构建需要根据电池内部构造将电池等效为若干基本电气元 件的组合;常见的电池等效电路模型主要有Rint模型、Thevenin模型、PNGV模型以及GNL模 型;
[0087] 本发明的建模方法基于电池的实际物理构造,并结合了电池的电化学特性;针对 液态或半液态金属电池的结构,开展电化学阻抗谱研究并进行阻抗谱拟合;实施例中对液 态金属电池进行阻抗谱拟合获得的曲线图如图1所示,其中方形为实验数据,圆形为Zview 软件拟合数据;
[0088] 阻抗谱拟合电路包含欧姆内阻Ro、极化电阻心、极化电容Cp以及韦伯阻抗Ws;韦伯 阻抗Ws是由扩散过程所致,可将其等效为多个RC环串联;实施例中,采用一个RC环近似替代 Ws,即扩散等效电阻Rw和扩散等效电容(^,所得模型为如图2所示的二阶Thevenin模型。
[0089] 端口电压表示电池端电压,在电池内部,电池电动势U。。与欧姆内阻Ro和两个RC环 串联,其中一个RC环是由代表电池极化过程的极化电阻心和极化电容心并联组成,另一个RC 环是由代表电池扩散过程的扩散等效电阻Rw和扩散等效电容Cw并联组成。
[0090] (2)采用对称脉冲对液态或半液态金属电池进行HPPC测试,获取用于辨识 Thevenin等效电路模型开路电压、电阻、电容参数的电池实测数据,包括电池电流、电压、充 电容量、充电时间、放电容量和放电时间的测试数据;
[0091] 实施例中,测试对象是容量为20Ah的液态金属电池,其负极材料为Li,正极材料为 Sb基合金;测试设备选用Arbin5V/300A测试仪和蓝电5V/20A测试仪,内置多个独立可编程 通道,通过自由设置恒流、恒压、静置等工步获取响应曲线,完成数据采集;
[0092] 实施例中,根据美国Freedom CAR功率辅助型电池测试手册,对液态或半液态金属 电池进行HPPC测试,通过在放电过程的某些SOC点上使用反馈脉冲来测试其动态功率能力; 其步骤⑵具体如下:
[0093] (2.1)用0.2C的额定电流对电池进行3次完整的充放电;每次放电的平均时间为T;
[0094] (2.2)将电池充满,将电池静置tr时间;实施例中,tr= Ih;
[0095] tr为IOmin〜lh,由于液态或半液态金属电池的内部传质速度比固态电池快,超过 该时长可能使测试结果受自放电影响;
[0096] (2.3)放电0.1T的时长,静置Ih,施加对称脉冲信号;
[0097] 实例中,对称脉冲幅值Ip = O. 3C;施加对称脉冲信号的具体工步为:0.3C放电IOs; 静置40s; 0.3C充电IOs;静置40s;
[0098] (2.4)重复步骤(2.3)直至将电池满放。
[0099] 采用对称脉冲对电池进行测试,充、放电脉冲幅值相同;而标准的HPPC混合功率脉 冲测试工步是采用IC的放电脉冲和0.75C的充电脉冲,是非对称的;在液态金属电池的实际 测试过程中,采用非对称脉冲不容易获取电池极化参数的变化规律,因此选用对称脉冲;
[0100] 另一方面,根据液态或半液态金属电池的大电流放电特性和实际容量选取脉冲幅 值为0.2C-0.5C之间,而非标准的1C,能够保证电池在测试期间的安全稳定运行,避免液态 或半液态金属电池在大倍率(>1〇Α)下由于绝对电流过大可能对电池性能造成的不良影响。
[0101] 图3所示,为对实施例中液态金属电池进行HPPC测试获得的电压波形,纵坐标 Vol tage为电压,横坐标Time为测试时间点;从该图可清楚分辨出恒流放电、静置、脉冲等工 步引起的电压变化。
[0102] (3)根据电池实测数据拟合出Ro、Rp、Cp、Rw、Cw和U。。与SOC (荷电状态)的函数式,即 为辨识出的等效电路模型参数;辨识工具选用Ma tlab/Simulink中的Parameter Estimation,这种方式相比于Excel有更高的灵活性,模型变动只需修改电路图即可;算式 拟合工具选用Istopt;该步骤具体如下:
[0103] (3· 1)在Matlab/Simulink中搭建Thevenin等效电路模型,用powergui将米样时间 设为1.014s;由于测试仪采样设置为ls,但存在0.014s的测量误差,因此在采用时间的设置 里补偿该误差;
[0104] (3.2)用Inport模块表示开路电压U〇c,0utport模块表示电池端电压,电阻、电容等 参数均用字母赋值;譬如,在Ro的赋值区域输入R〇,以便将其设置为变量;
[0105] (3.3)从测试数据中取每个SOC点静置Ih后的开路电压值(OCV)对电动势U。。进行赋 值,BP Ucic = OCV;
[0106] (3.4)截取测试数据中脉冲工步获取的时间、电流和电压数据输入Matlab的工作 空间;导入Parameter Estimation工具箱,其中时间为从0开始的等差数列;将电池欧姆内 阻Ro、极化电阻心、极化电容Cp、扩散等效电阻1^和扩散等效电容(^参数设为变量;
[0107] (3.5)通过工具箱的曲线拟合辨识参数,在每个施加过脉冲的SOC点获取一组辨识 出的变量值;
[0108] (3.6)采用Istopt软件对SOC值和与其对应的辨识出的变量值进行算式拟合,获得 电池欧姆内阻Ro、极化电阻RP、极化电容CP、扩散等效电阻Rw和扩散等效电容Cw参数与SOC的 函数关系式;
[0109] 并采用Istopt软件根据SOC值与步骤(3.3)中对应的电池电动势U。。的值获得U。。与 SOC的函数关系式。
[0110] 在步骤(3.6)中,对最终拟合的开路电压、电阻、电容参数与SOC的表达式采取分段 函数表示;对电池欧姆内阻Ro与SOC的关系分段拟合为:
[0111]
Figure CN105912799BD00111
[0112] 极化电阻心与50(:的关系分段拟合为:
[0113]
Figure CN105912799BD00112
[0114] 极化电容心与50(:的关系分段拟合为:
[0115]
Figure CN105912799BD00113
[0116] 扩散等效电阻IV^SOC的关系分段拟合为:
[0117]
Figure CN105912799BD00114
[0118] 扩散等效电容(^与SOC的关系分段拟合为:
[0119]
Figure CN105912799BD00115
[0120] 电池电动势U。。与SOC的关系分段拟合为:
[0121]
Figure CN105912799BD00116
[0122]其中,f1>2...n (SOC) Wm(SOC) Am(SOC) Jm(SOC)、rHn(SOC)、zHn(SOC) 都是自变量为SOC的函数;
[0123] 其基本形式均为:
Figure CN105912799BD00121
[0124] 其中,(3〇、(31、〇2">、(11、¥1、(12、¥2均为实数拟合系数,在每一个分段内,〇<y<6,分 段数量Kn<4。
[0125] 在实施例中,对于Ro的拟合分为4段进行,S卩η = 4,SOC分为0〜0.3、0.3〜0.88、 0.88〜0.96和0.96〜1共4段;由此提高欧姆内阻Ro的拟合精度;其余参数采用类似的分段 函数表达。
[0126] 采用分段函数可以提高对电池模型各参数的辨识精度;其中,提高欧姆内阻R0的 辨识精度对液态或半液态金属电池储能特性建模来说至关重要;因为液态或半液态金属电 池的特性与锂电、锂硫等固态电池并不相同,在锂离子电池中,随着放电倍率的变化,Ro变 化不大,可认为相同放电状态下Ro为常数,但在液态或半液态金属电池中,放电倍率会影响 内部温度和熔盐体积,若以同样的Ro计算,则会造成较大误差;
[0127] 在本步骤中,相比于传统参数的单一软件辨识,采用Matlab/Simulink以及Istopt 软件进行分工辨识,有利于提高获得的模型的辨识精度。
[0128] 图4所示,为SOC = 0.7时的参数辨识情况;纵坐标Vo I tage为电压,横坐标Time为测 试时间点,Measured是指实测值,Simulated是指仿真值;该图中,虚线(Measured)为从图3 中放大截取的脉冲电压数据,实线(Simulated)为Simulink参数辨识拟合曲线,二者重合度 很高。
[0129] SOC = O.7时的Simulink拟合计算端电压与真实端电压的相对误差,如图5所示;图 5中,纵坐标Error为电压误差,横坐标Time为测试时间点;电压误差不超过±0.2%。
[0130] 图6为拟合的SOC-OCV曲线;纵坐标OCV代表开路电压,同时也是电动势U。。的值,横 坐标为荷电状态S0C,该图中,Measured是指实测值,Fit是指拟合值;这条曲线对应着辨识 出的U。。和SOC的函数关系式。
[0131] ⑷采用安时法计算获得液态或半液态金属电池的S0C;
[0132] 实施例中,步骤⑷具体为将充、放电电流i (t)作为输入,用安时法获取S0C,安时 法是指用电流对时间的积分计算充、放电容量;用安时法获取的SOC关系式如下:
[0133]
Figure CN105912799BD00122
(7)
[0134] 其中,to为起始积分时刻,t为本时刻,Qo为额定放电容量;
[0135] (5)根据Rakhmatov模型获得修正式;采用修正式对采用安时法获得的SOC进行修 正,获得修正SOC以减小因电池的倍率容量效应和恢复效应产生的容量误差;
[0136] 在不同的放电倍率下电池所能放出的容量并不相同,电流越大,放出的容量越少, 称之为电池的倍率容量效应;电池放电一段时间后静置,容量会有少量恢复,称之电池的恢 复效应,本发明所建立的模型较好地补偿了电池的倍率容量效应和恢复效应。
[0137] 步骤⑸具体如下:
[0138] (5.1)根据Rakhmatov模型获取修正式:Rakhmatov模型给出了负载电流i (τ)和电 池寿命L的关系式:
[0139]
Figure CN105912799BD00131
(8)
[0140] 其中α为电池可以输出的最大容量,β为电池参数,表示在电极表面活动载流子被 补偿的速率,它可以用来衡量电池放电特性,不同种类电池的β不同,m为求和项数;
[0141] 当负载电流i (τ)等于恒定电流I,且P2L彡1时:
Figure CN105912799BD00132
(9)
[0142] y
[0143] 式(9)中,1⑴为安时法计算的容量,u⑴不可利用容量;当电池放电时,电池总容 量的一部分供给负荷,另一部分转化为不可利用容量,电流越大,不可利用容量越大,放出 的电量越少;
[0144] 而当电池放电^时间后静置时,u⑴会发生衰减,此时u (t)如以下式(10)所示,不 可用容量转化为可用容量;
[0145]
Figure CN105912799BD00133
(10)
[0146] 其中,ti为静置时间,ti为静置前的充、放电时间;
[0147] 实施例中,恒流公式(9)中u (t)与静置时的式(10)从本质上是相同的,即不管是恒 流还是静置,u⑴都可以采用式(10)表达。
[0148] 当负载电流i (τ)为变电流时,将变电流分解为所有采样时间的叠加;由于ti和ti是 随时间变化的,在每个采样时间内,不可利用容量u⑴理论上都要重新计算一次;这种非迭 代的指数级数运算不仅浪费程序运行时间,也容易造成程序卡顿;为了将非迭代的运算转 化为迭代运算,令放电时间t等于采样时间ts,并将f(ti)用Istopt软件拟合为少数指数项 和的形式:
[0149]
Figure CN105912799BD00134
(11)
[0150] 其中,81、131、32』2^1、1^为;^1^)的拟合系数,1为拟合项数,1越大则公式越精确, 在误差允许的范围内可灵活调整项数X的大小。
[0151] 根据式(11),每进入一个新的采样时间,认为电池在本采样时间静置了零秒,因此 增加21 · f (0) =21 (ai+a2+…+ax)的不可利用容量,除此之外,之前累计的不可利用容量u (t)在本采样时间内静置了。秒,会发生衰减,即第一项是上一时刻的
Figure CN105912799BD00135
ί咅,第二项是上一 时刻的
Figure CN105912799BD00136
倍,以此类推;则迭代公式表示如下:
[0152]
Figure CN105912799BD00137
[0153] 其中,Uk为k时刻之前电池内部由于倍率影响产生的不可利用的容量,Ik为k时刻电 池电流,ts为米样时间;
[0154] (5.2)采用式(11)和式(12)对步骤(4.1)中的基本电池模型的SOC容量进行修正, 其函数式如下:
[0155]
Figure CN105912799BD00141
( 13)
[0156] 对SOC进行修正,可减少因倍率容量效应和恢复效应造成的容量误差。
[0157] 步骤(5.1)将该模型转化为更易实现的迭代式,将恒流与变电流的两种情况统一, 且去掉了求和符号,大大减少了计算量,可在很大程度上节省计算处理的时间。
[0158] ⑹根据修正S0C,以及在步骤(3)获得的电动势U。。、欧姆内阻Ro、极化电阻心、极化 电容CP、扩散等效电阻1和扩散等效电容CV^SOC的函数关系,对等效电路模型的参数进行 修正。
[0159] 实施例中,将修正后的SOC (t)的值代入各参数与SOC的函数式,可以获取对应的修 正后的办為、〇)、1^、(^和1]。。的值。
[0160] 实施例中,根据电池端电压对建立的等效电路模型与等效电路模型的参数进行验 证,具体如下:
[0161] (7.1)根据步骤(6)获得的等效电路模型的参数,获取欧姆内阻Ro两端的欧姆压降 Uo、极化压降UP、扩散等效压降Uw,则电池端电压如下式所示:
[0162] Ubat = Uoc-Uo-Up-Uw (14)
[0163] 其中,1^为电池端电压,Up为Rp与Cp构成的RC环的两端极化压降,Uw为Rw、Cw构成的 RC环的两端的扩散等效压降;
[0164] (7.2)将所述电池端电压1^*仿真值与在各工况下测试获得的电池端电压测试数 据进行比对,通过仿真值与实测值的匹配程度,验证建模获得的所述等效电路模型与等效 电路模型的参数与所述电池的匹配程度;仿真值与实测值越接近,表明建模获得的等效电 路模型与等效电路模型的参数与电池越匹配。
[0165] 综上而言,实施例中,电池建模方法的流程如图7所示,首先构建电池等效电路模 型,再用测试仪获取电池实测数据,用以辨识模型各参数,并以电流为输入,用安时法计算 容量,然后在安时法的基础上对容量进行修正,使得模型各参数的值也得到进一步修正;最 终根据获得的电池端电压,在恒流工况和HPPC工况下验证电路模型及模型参数的准确性。
[0166] 图8所示,是在恒流工况下,对应不同充放电幅值下,本实施例建立的液态金属电 池模型获得的电池端电压仿真值与实测值的对比图,纵坐标Voltage为电压,横坐标Time为 测试时间点;该图中,用三角形、十字形和圆形分别代表〇. 1C、0.2C和0.45C的充放值的电曲 线;仿真曲线均用实线表示,由该图可看出,模型仿真值与真实值紧密重合,经实测,实施例 建立的模型的误差不超过5%,表明本发明建立的模型较好地模拟了真实电池的倍率性能。
[0167] 图9所示,是在HPPC工况下,本实施例建立的液态金属电池模型获得的电池端电压 仿真值与实测值的对比示意图,纵坐标Voltage为电压,横坐标Time为测试时间点;该图中, Measured是指实测值,Simulated是指仿真值;虚线(Measured)为测试仪测得的真实电压曲 线,实线(Simulated)为模型仿真曲线,二者重合较好;将其中SOC = O. 5的脉冲部分放大,实 线和虚线重合度依然很高,误差较小;进一步验证了本发明建立的模型的准确性和实用性。
[0168] 本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以 限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含 在本发明的保护范围之内。

Claims (8)

1. 一种液态或半液态金属电池的建模方法,其特征在于,包括如下步骤: (1) 通过电化学测试获取所述电池的阻抗谱,并根据所述阻抗谱进行拟合,构建所述电 池的阻抗谱拟合电路;根据所述阻抗谱拟合电路构建所述电池的等效电路模型; (2) 采用对称脉冲对所述电池进行混合脉冲功率性能测试,获取所述电池电流、电压、 充电容量、充电时间、放电容量和放电时间的测试数据; (3) 根据所述测试数据拟合所述电池的电动势U。。、欧姆内阻R〇、极化电阻心、极化电容 CP、扩散等效电阻Rw和扩散等效电容CV^SOC的函数关系,辨识出所述等效电路模型的参数; ⑷采用安时法获得所述电池的S0C; ⑸对所述SOC进行修正,获得修正SOC,以补偿所述电池的容量误差; ⑹根据所述修正SOC以及所述函数关系对所述等效电路模型的参数进行修正; 所述步骤⑸包括如下子步骤: (5.1) 通过对Rakhmatov模型进行演绎,获得修正式
Figure CN105912799BC00021
其中,X为拟合项数,Uk为k时刻之前电池内部由于倍率影响产生的不可利用的容量,Ik 为让时刻电池电流,1^为米样时间;31、131、32、&amp;2〜3\、&amp;\为;1^(1^)的拟合系数,111,1{-1表不第一项 k-Ι时刻之前电池内部由于倍率影响产生的不可利用的容量,Ulk-l·表示第二项k_l时刻之前 电池内部由于倍率影响产生的不可利用的容量,Ux,k-ι表示第X项k-ι时刻之前电池内部由于 倍率影响产生的不可利用的容量; 其中,
Figure CN105912799BC00022
的拟合式,
Figure CN105912799BC00023
β为电池参数,ti为静置时间;ti为静置前的充、放电时间,m为求和项数; (5.2) 采用所述修正式对所述电池SOC进行修正,获得
Figure CN105912799BC00024
其中,to为起始积分时刻,t为当前时刻,Qo为电池额定放电容量,SOC⑴为t时刻电池 SOGSOC/⑴是指t时刻修正的SOC,SOC (to)为to时刻电池S0C,i⑴为充放电电流。
2. 如权利要求1所述的建模方法,其特征在于,所述步骤⑵包括如下子步骤: (2.1) 采用额定电流对所述电池进行N次完整的充放电;N多3; (2.2) 将所述电池充满并将电池静置tr时长; (2.3) 将所述电池放电0. IT时长并静置tr时长,再对电池施加对称脉冲信号;其中,T为 步骤(2.1)中每次放电的平均时间; (2.4) 重复步骤(2.3)直至将所述电池满放;获得电池电流、电压、充电容量、充电时间、 放电容量和放电时间的测试数据。
3. 如权利要求2所述的建模方法,其特征在于,所述步骤(2.3)中,对称脉冲信号采用幅 值为Ip的短时阶跃信号;IpSO . 2C〜0.5C。
4. 如权利要求2或3所述的建模方法,其特征在于,所述步骤(2.3)中对电池施加对称脉 冲的具体工步为:Ip放电IOs,静置40s; Ip充电IOs,静置40s。
5. 如权利要求1所述的建模方法,其特征在于,所述步骤⑶包括如下子步骤: (3.1) 采用仿真工具搭建等效电路模型;并在仿真工具中对等效电路模型的参数赋值; (3.2) 从所述测试数据中获取每个SOC点静置tr时长后的开路电压值,并对电动势U。。进 行赋值; (3.3) 以测试数据中脉冲工步获取的时间、电流和电压数据作为仿真工具的输入参数; 并将电池欧姆内阻Ro、极化电阻心、极化电容Cp、扩散等效电阻1和扩散等效电容Cw#数设为 仿真工具里的变量; (3.4) 通过仿真工具进行曲线拟合以辨识参数,在每个施加过脉冲的SOC点获取一组辨 识出的变量值; (3.5) 采用算式拟合工具对SOC值和与其对应的辨识出的变量值进行算式拟合,获得电 池电动势U。。、欧姆内阻Ro、极化电阻心、极化电容Cp、扩散等效电阻Rw和扩散等效电容Cw参数 与SOC的函数关系。
6. 如权利要求5所述的建模方法,其特征在于,所述步骤(3.5)中,采用分段函数拟合欧 姆内阻Ro、极化电阻心、极化电容Cp、扩散等效电阻Rw和扩散等效电容。参数与SOC的函数关 系。
7. 如权利要求1所述的建模方法,其特征在于,所述步骤⑷具体为:采用安时法,将充 放电电流i (t)作为输入,获得
Figure CN105912799BC00031
8. 如权利要求1所述的建模方法,其特征在于,其步骤⑹后,还包括步骤(7):根据电池 端电压对所述等效电路模型与等效电路模型的参数进行验证,具体如下: (7.1) 根据步骤(6)进行修正后获得的等效电路模型的参数,获取欧姆内阻Ro两端的欧 姆压降Uo、极化压降Up、扩散等效压降Uw;并获取电池端电压Ubat = Ucic-Uo-Up-Uw; (7.2) 将所述电池端电压1^*仿真值与在各工况下测试获得的电池端电压1^*实测值进 行比对,通过所述仿真值与实测值的匹配程度,验证建模获得的所述等效电路模型与等效 电路模型的参数与所述电池的匹配程度。
CN201610268697.XA 2016-04-27 2016-04-27 一种液态或半液态金属电池的建模方法 Active CN105912799B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610268697.XA CN105912799B (zh) 2016-04-27 2016-04-27 一种液态或半液态金属电池的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610268697.XA CN105912799B (zh) 2016-04-27 2016-04-27 一种液态或半液态金属电池的建模方法

Publications (2)

Publication Number Publication Date
CN105912799A CN105912799A (zh) 2016-08-31
CN105912799B true CN105912799B (zh) 2018-12-14

Family

ID=56752835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610268697.XA Active CN105912799B (zh) 2016-04-27 2016-04-27 一种液态或半液态金属电池的建模方法

Country Status (1)

Country Link
CN (1) CN105912799B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446480B (zh) * 2016-12-06 2019-03-29 上海动力储能电池系统工程技术有限公司 锂离子电池电极材料开路电势曲线的拟合方法
CN106872905A (zh) * 2017-02-23 2017-06-20 哈尔滨工业大学 一种单体锂离子全电池参数获取方法
CN108663619B (zh) * 2017-03-27 2020-11-17 宁德时代新能源科技股份有限公司 确定电池工作电压曲线的方法、装置和设备
CN107248597B (zh) * 2017-06-23 2019-06-28 华中科技大学 一种液态金属电池的建模方法
CN109273781B (zh) * 2017-07-18 2021-03-09 宁德时代新能源科技股份有限公司 电芯监测方法和电芯监测装置
CN108090244B (zh) * 2017-11-07 2021-09-07 国网江苏省电力公司盐城供电公司 一种并联型锂离子电池系统建模方法
CN111123107A (zh) * 2018-10-30 2020-05-08 北京天诚同创电气有限公司 电池仿真建模方法、装置及电池等效模型
CN110082683A (zh) * 2019-05-09 2019-08-02 合肥工业大学 抑制安时积分soc估计误差的闭环补偿方法
CN110058162B (zh) * 2019-05-24 2020-12-01 中国科学院电工研究所 一种基于线性时不变电池模型结构的参数辨识方法
CN113138340A (zh) * 2020-01-17 2021-07-20 华为技术有限公司 电池等效电路模型的建立方法、健康状态估算方法及装置
CN111474431A (zh) * 2020-04-21 2020-07-31 三峡大学 一种电解电容器分数阶等效电路模型及其参数辨识方法
CN111487535B (zh) * 2020-04-29 2021-05-18 华中科技大学 一种液态金属电池双等效电路模型的参数获取及切换方法
CN111856285B (zh) * 2020-07-06 2021-06-08 大连理工大学 一种电动汽车退役电池组等效模型建模方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102088118A (zh) * 2010-12-28 2011-06-08 深圳市航盛电子股份有限公司 一种电池管理系统、电动车及荷电状态的估算方法
CN102569922A (zh) * 2012-03-05 2012-07-11 同济大学 一种基于单体电池一致性的蓄电池组soc估计改进方法
CN202794475U (zh) * 2012-06-28 2013-03-13 浙江大学城市学院 混合动力车用动力锂电池soc估计器
CN103472397A (zh) * 2013-08-09 2013-12-25 江苏欧力特能源科技有限公司 铅酸电池模型经验公式法中关键参数抗差辨识方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102088118A (zh) * 2010-12-28 2011-06-08 深圳市航盛电子股份有限公司 一种电池管理系统、电动车及荷电状态的估算方法
CN102569922A (zh) * 2012-03-05 2012-07-11 同济大学 一种基于单体电池一致性的蓄电池组soc估计改进方法
CN202794475U (zh) * 2012-06-28 2013-03-13 浙江大学城市学院 混合动力车用动力锂电池soc估计器
CN103472397A (zh) * 2013-08-09 2013-12-25 江苏欧力特能源科技有限公司 铅酸电池模型经验公式法中关键参数抗差辨识方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Research on Grid-connected Technology of Energy Storage System with Liquid Metal Batteries;王大磊等;《Computer Science Research》;20160131;第50卷;第568页 *
纯电动汽车电池管理数据分析系统的研究;纪淼淼;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20150315(第03期);第15-16页 *
采用等效电路的参数自适应电池模型及电池荷电状态估计方法;宁博等;《西安交通大学学报》;20151031;第68-69页 *
锂离子电池特性建模与SOC估算研究;卢杰祥;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130115(第1期);第38-54页 *

Also Published As

Publication number Publication date
CN105912799A (zh) 2016-08-31

Similar Documents

Publication Publication Date Title
CN105912799B (zh) 一种液态或半液态金属电池的建模方法
CN104392080B (zh) 一种锂电池分数阶变阶等效电路模型及其辨识方法
CN102508165B (zh) 一种评价磷酸铁锂电池自放电一致性的方法
CN103926538B (zh) 基于aic准则的变阶数rc等效电路模型及实现方法
CN103632018B (zh) 一种基于Simscape平台的电池建模方法
CN109586373A (zh) 一种电池充电方法和装置
CN105759213A (zh) 一种测量蓄电池剩余容量soc的方法
CN104965179B (zh) 一种锂离子蓄电池的温度组合电路模型及其参数识别方法
CN107991623A (zh) 一种考虑温度和老化程度的电池安时积分soc估计方法
CN104051810B (zh) 一种锂离子储能电池系统soc估算快速修正方法
CN106250576A (zh) 一种基于动态阻抗的锂电池模型的建模方法
CN105954679A (zh) 一种锂电池荷电状态的在线估计方法
CN105445663B (zh) 电池老化程度的检测方法和装置
CN104849672A (zh) 基于等效电路模型的锂电池动态阻抗参数识别方法
CN105974320B (zh) 一种液态或半液态金属电池荷电状态估计方法
CN105093129B (zh) 一种储能电池剩余容量检测方法
CN107121643A (zh) 锂离子电池健康状态联合估计方法
CN103698716A (zh) 一种基于衰减系数的串联电池组可放出电量衰减评价方法
CN109614754A (zh) 一种锂离子电池三维简化模型的仿真方法
CN106124996A (zh) 一种锂离子电池单体的一致性判定方法及装置
CN106093782A (zh) 动态建模的最小二乘支持向量机soc估计方法
Qiuting et al. State of health estimation for lithium-ion battery based on d-ukf
CN105738821B (zh) 准确计算不同温度下电池库仑效率的方法
CN108829911A (zh) 一种开路电压与soc函数关系优化方法
Kim et al. Hysteresis modeling for model-based condition monitoring of lithium-ion batteries

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant