锂离子电池数字孪生工作流中多阶等效电路模型与经验老化耦合的基准评测

《Journal of Energy Storage》：Benchmarking of multi-order equivalent-circuit models with empirical ageing for Li-ion battery digital-twin workflows

　　

准确且计算高效的电池模型对于电池管理系统（BMS）与电池数字孪生工作流至关重要。然而模型选择需谨慎权衡动态电压保真度、计算开销、能量计量与老化表征之间的关系。研究人员提出一种可复现的基于MATLAB的基准评测工作流，用于在共同参数化与老化假设下比较电池包尺度的等效电路模型（ECM）。该工作采用SAFT型RC模型以及一阶、二阶、三阶Thevenin模型，基于开放获取的荷电状态（SoC）–健康状态（SoH）依赖参数曲面、经验日历–循环老化定律，并通过DATTES工具箱实现自动化后处理。基准对象为64s29p磷酸铁锂（LiFePO?）电池包，在固定25 °C下施加重复的24 h合成电流工况，等效运行三年。结果表明高阶ECM对三阶Thevenin内部参考的电压一致性更优，但超越二阶后的边际收益递减；在所考察的模型族、工况与MATLAB实现中，二阶Thevenin模型在内部参考电压保真度与单步计算成本间取得最有利折衷。老化输出在各ECM间呈现相同的SoH与内阻增长轨迹，因为所采用老化律由共同日应力指标驱动，而非ECM特定的极化态或瞬态过电位历史；因此这些老化结果应解释为所实现模块化老化耦合策略的属性，而非ECM阶数不影响真实电池退化的普遍证明。该研究为受控ECM比较提供透明可复现工作流，而实验验证、电热耦合及估计器在环测试在预测性部署前仍属必要。

研究背景方面，全球电动交通、分布式储能与电气化趋势提升了对锂离子（Li-ion）电池准确、可解释且计算高效模型的需求。嵌入式BMS与电池数字孪生工作流要求模型在长周期下数值稳健且计算可行。等效电路模型（ECM）因物理可解释性与计算效率的平衡而被广泛采用，其中SAFT型RC网络与基于电压的Thevenin模型分别代表表面–本体电荷重分布与多时间常数极化动力学。已有研究指出一阶与二阶Thevenin常在精度–复杂度间具优势，但将ECM族、经验老化律、开放参数曲面与标准化后处理集成进透明可复现的长周期包级基准仍有限。本研究即在此缺口下开展，旨在以受控方式量化ECM族与阶数在相同参数化、老化律与工况下对电压响应、能量吞吐、往返效率、退化指标及计算开销的影响。论文发表于《Journal of Energy Storage》。

关键技术方法上，研究人员搭建统一MATLAB仿真框架，对64s29p LiFePO?电池包在固定25 °C下重复合成24 h电流工况（含夜间充电、晨间放电、午间脉冲与待机噪声），以60 s步长模拟约150万步等效三年。模型涵盖SAFT型2C–2R约化RC模型与一阶至三阶Thevenin模型；电气参数取自开放SoC–SoH参数曲面（Seger等第二寿命数据集，SoH域0.5–0.8，插值坐标钳制于0.8），RC–SAFT内部C_c、C_b、R_e由一阶Thevenin时间常数τ₁=R₁C₁解析匹配，容量拆分因子α=0.20做鲁棒性扫描；老化采用Batteries2020经验日历–循环律，日应力（存储时长、DoD、Mid-SoC、温度）驱动容量衰减ΔC与内阻增量ΔR_int，经DATTES工具箱自动后处理能效、电压RMSE/NRMSE/TVE等指标。

研究结果部分，2. 建模框架与方法学中，2.1节给出从参数初始化到后处理的统一流程；2.2节定义RC–SAFT状态向量x=[SoC;V_c;V_b]与Thevenin各阶极化电压离散更新式V_p,i,k+1=α_iV_p,i,k+R_i(1?α_i)I_k，α_i=exp(?T_s/(R_iC_i))；2.3节说明双线性插值边界处理使所有ECM共享V_OC与R_int曲面，三阶Thevenin第三支路由二阶参数解析派生仅作内部高阶参考；2.4节耦合日历老化率k_cal(T,SoC)=A_calexp(?E_a/RT)/(a+b·SoC)与循环老化ΔC_cyc=N_cycles·k_cyc(DoD,T,Mid-SoC)，日容量损失ΔC_day=ΔC_cal+ΔC_cyc，内阻增量ΔR_int,day=k_rint·ΔC_day；2.5–2.6节强调60 s定步长指数更新数值稳定、全链路可复现。

1. 3.
   结果与性能分析中，3.1节给出包技术参数（单体60 Ah、标称204.8 V、11.7 kWh可用能）与合成24 h工况；3.2节电压动力学显示第1095天末端工况下RC–SAFT在高功率暂态与弛豫段系统性低估，Th1仍有延迟整定，Th2最贴近Th3内部参考，Th3因第三支路非独立辨识故增量改善有限；3.3节能量与效率表明总处理能量序为RC–SAFT>Th1>Th2>Th3（Th3约24.81 MWh、RC–SAFT约25.09 MWh），平均往返效率RC–SAFT 92.51%、Th1 91.53%、Th2 89.90%、Th3 89.09%，低阶因未充分解析长尾弛豫而表观更高效，α扫描证实RC–SAFT能耗对拆分因子不敏感（效率波动0.0005 pp）；3.4节老化轨迹所有ECM完全一致，终态SoH=83.47%、ΔR_int=+19.83%，敏感性分析显示循环损失尺度系数主导SoH（±20%扰动致SoH 79.84–87.10%），k_rint仅缩放电阻增长，日历系数与SoC系数影响微弱；3.5节内阻增长与效率退化在等温下呈近线性（ΔR_int每增1%效率降约0.09 pp）；3.6节计算性能以每步微基准为主指标，RC–SAFT相对耗时1.00、Th1 0.48、Th2 0.50、Th3 0.54，RC–SAFT因需重算等效电容与耦合电阻反重于其概念阶数，Th族随阶数单调增但绝对差小；3.7节以Th3为内部参考的RMSE/NRMSE/TVE确认Th2相较RC–SAFT与Th1显著拉近与高阶参考差距。

讨论与结论部分，研究人员明确指出本研究非论证某ECM拓扑普适最优，而是把已建立模型结构、开放SoC–SoH参数曲面、经验日历–循环老化与DATTES后处理集成进单一MATLAB基准工作流，使相同工况下模型族行为可比。核心结论为：在所用参数曲面、老化律、合成24 h工况、64s29p LiFePO?包与MATLAB实现内，二阶Thevenin模型在内部参考电压保真度与单步计算成本间提供最有利折衷，超越二阶的边际电压收益不抵复杂度；各ECM老化轨迹相同是模块化老化耦合（由共同日应力而非ECM极化态驱动）的方法论属性，不得外推为“ECM阶数不影响真实退化”；该基准工作流本身透明可复现，但绝对电压精度、效率值与SoH预测均非实包实验验证，电热耦合、多化学体系、估计器在环与实测动态IV校验属部署前必要扩展。