《基于数理统计方法的锂电池电解液电导率优化设计.docx》由会员分享,可在线阅读,更多相关《基于数理统计方法的锂电池电解液电导率优化设计.docx(14页珍藏版)》请在第一文库网上搜索。
1、基于数理统计方法的锂电池电解液电导率优化设计摘要:离子电导率是评估锂离子电池电解液性能的重要特征参数,直接影响着电池的低温、倍率等性能,在电解液的设计中极具指导价值。传统的电解液研发模式主要基于经验和实验的“试错法”,存在变量多、实验成本高、开发周期长等问题。针对以上问题,本工作提出了一种结合空间填充混料设计与高斯过程回归的电导率优化设计方法,以包括环状碳酸酯(EC)及不同种类线性碳酸酯、竣酸酯的电解液溶剂组成作为模型的输入,电导率作为模型的输出,并运用最大似然估计求解超参数;通过后续实验验证了模型的有效性,并可预测满足任意电导率要求的电解液溶剂配方。关键词:电导率;空间填充混料设计;高斯过程
2、回归;锂离子电池;电解液目录前言11 .实验材料和方法21.1. 分析测试仪器21.2. 电解液配制22 .实验设计32. 1.混料设计33. 2.空间填充设计34. 3.空间填充混料设计43 .数据处理与分析61. 1.高斯过程回归模型63. 2.模型结果分析64 .结论13前言电解液被称为“锂电池的血液”,其作用是在正负极间传输锂离子,对电池的能量密度、循环寿命、安全性能、高低温性能有直接影响。电解液的重要指标之一是电导率,测量电导率随温度和导电盐浓度的变化有助于确定存储设备的极限工作条件(例如,最高、最低工作温度,最大充电速率)。而极限工作条件的研究对快充技术、深空探测、极地、航天技术等
3、特殊应用场景具有重要意义。早期研究通常采用经验或半经验方法,将一类介电常数大的溶剂和另一类黏度低的溶剂通过调整混合比,以获得电导率高、电化学性能优良的电解液。但电解液的组成相对复杂,变量较多,常规的改变单个变量的实验方法不仅无法考虑到变量之间的相互作用,且存在实验成本高、研究周期长的问题,基于数理统计方法的试验设计能够科学有效地提升研发效率。试验设计通过对试验进行合理安排,以较小的试验规模、较短的试验周期和较低的试验成本,获得理想的试验结果以及得出科学的结论。常见的试验设计有经典筛选设计、响应曲面设计、田口设计、正交超饱和设计和混料设计等。空间填充设计能够均匀分布设计点,可极大地降低试验的次数
4、。而混料设计的优势在于研究变量之间是相互关联的,指定变量之和等于1,满足了电解液研发对溶剂比例影响的探究。结合了两种设计方法的空间填充混料设计能同时进行多变量分析,提高实验效率,然而这种设计带来了多维的实验数据,传统的数据拟合模型如最小二乘法无法全面归纳实验数据中的科学规律。作为一种能根据训练数据的量调整模型表达能力的非参数方法,高斯过程回归在进行此类数据处理时具有显著优势。基于此,本工作以环状碳酸酯(EC)及不同种类线性碳酸酯、竣酸酯等常见的电解液溶剂体系为研究对象,按照空间填充混料设计得到的电解液溶剂组成配制电解液进行电导率的测试,以高斯过程回归对数据进行分析处理,得到预测电导率的数学模型
5、,并对其进行有效性的评价。1 .实睑材料和方法1.1. 分析测试仪器电解液样品分别在10、25和45下手套箱中恒温静置4h再用雷磁电导率仪(DDBJ-350F)测试离子电导率。工作电极为伯电极,电极常数为1。1.2.电解液配制在充满氧气气氛的手套箱中将溶剂碳酸乙烯酯(EC)、碳酸甲乙酯(EMC)、碳酸二甲酯(DMC)、碳酸二乙酯(DEC)、乙酸乙酯(EA)、丙酸乙酯(EP)、丙酸甲酯(MP)按照实验设计要求的质量比配成混合溶剂,然后将锂盐六氟磷酸锂1iPF6溶解于上述溶液中配制成ImoIkg的1iPFe电解液。手套箱中水、氧含量均控制在0.00001%以下。2 .实验设计1. 1.混料设计混料
6、设计保证了指定变量比例之和等于U本设计中存在3组混料设计。行业设计经验可确定环状碳酸酯、线性碳酸酯与陵酸酯的组成区间,将其质量分数分别设置为0.10.3、0.20.6以及0.10.40电解液中常用的环状碳酸酯包括碳酸乙烯酯(Ee)和碳酸丙烯酯(PC),然而PC对天然石墨的结构有一定的破坏作用,在石墨类电极中用量通常较低。考虑到更广泛的适用性,在本工作中环状碳酸酯只选择EC作为研究对象。因此,第一组混料设计是存在着线性约束的混料设计;第二组混料设计是3种线性碳酸酯EMC、DMC、DEC的组合;第三组混料设计是3种线性碳酸酯EA、MP、EP的组合。各溶剂在质量分数上满足如下关系:环状碳酸酯+线性碳
7、酸酯+竣酸=1DMC+DEC+EMC=1EA+MP+EP=12. 2.空间填充设计空间填充设计的算法最初会在指定设计区域内生成大量随机点,随后使用快速Ward算法将这些点聚类到个数与指定的试验次数相同的聚类中。使用MaXPrO(最大投影)最优性准则来获取最终设计点,对于P个变量以及次试验则有:CMaXPr。1:=id-Xk)式中,X为设计区域内的点的位置;Xjk是指第i次试验变量k的点的位置;Xjk是指第j次试验变量k的点的位置。因此最小化JMaXPr。可使任意两个设计点间距离最大。2.3.空间填充混料设计综合以上两种设计的思想,本次设计在质量分数满足混料设计的限制下以空间填充设计的算法均匀分
8、布10个实验点。图1(a)展示了电解液溶剂种类组成质量配比的三元图,图1(b)则将线性碳酸酯和假酸酯的组成质量配比集成在一张图中。环状碳酸酯(EC)(a)电解液溶剂种类EMC(EA)(b)线性碳酸酯组成和殁酸酯组成图1空间填充混料设计的质量分数三元图3.数据处理与分析1. 1.高斯过程回归模型高斯过程回归在处理高维度、小样本、非线性等问题上具有很强的自适应功能,与其他机器学习的方法比较,具有易实现、超参数自适应获取、输出具有置信区间等优点。如果把数据点确定为高斯过程的采样点,且输入向量通过隐藏函数构成的集合可由一个均值向量(X)和核函数k(x,X)决定,其中X、X为随机变量,则这一过程为高斯过
9、程回归。因此高斯过程回归定义为:f(x)GP(%,k(x,x)为了表达上的简洁,通过数据预处理使均值函数为0。考虑模型,其中为输入向量,为函数值,为考虑噪声后的观测值,假设噪声;给定一个新的样本输入,则相对应的预测值为,则可得观测值和预测值的联合先验分布为:1NoiN,S1oJ/n式中,。为维单位,(x,x)+k(x,)k(,x)k()矩阵元素k(m,)用来度量m和。之间的相关性。通过贝叶斯推断计算可出f的后验分布为:尸(fx,y,x)(f,cov(f)其中:T=k(x.,x)k(x,x)+1n1ycov(f)=k(x:x)-k(x:x)k(x,x)+ck(x,x)式中,预测值的均值,cv()
10、为预测值的方差。对于式(7)中的预测值r,本工作选用经典的高斯核函数的变形来确定高斯过程的协方差函数,函数表达式中无法确定的未知量称之为超参数。通过极大似然估计来求解核函数的超参数。3. 2.模型结果分析表1为空间填充混料设计的10个实验点的溶剂质量分数组成。图2为10个实验点的电导率实测值与预测值数据,图中的数据点落在y=x沿线,意味着高斯过程回归模型是真实函数的良好近似。表1空间填充混料设计溶剂组成-环状碳酸酹(EC)线性成陂心一一而EMCDECDMCEAMPEP0.2690.5500.1810.0430.9530.0040.1890.2020.6090.1290.5290.3420.14
11、00.2750.5850.0120.0880.90.2790.4220.2990.3720.4160.2120.4700.0150.5150.1660.5070.3270.9850.0130.0020.8470.0620.0910.2970.3050.3980.0220.3230.6550.6090.3550.0360.1560.60.2440.1140.8190.0670.5400.4390.0210.1180.5890.2930.5600.0390.4010.2600.3490.3910.2780.5660.1560.1930.0170.7900.0480.8990.0530.240.36
12、00.3920.7850.1880.0270.0820.640.2340.2110.4790.3100.0180.7280.2540.10.8300.169E0se-1614121086446810121416实测电导率/(mscm)图2实验点电导率实测值与预测值图3为电导率的变量影响图,可以直观看出每个模型输入变量在不同温度下对电导率的影响,着色带代表了95%置信区间。只有在同一组混料设计中的变量才能进行影响的比较,即不能将不属于同一混料设计的EMC和EA进行影响的比较。(b)线性碳酸酯组成OOOTOOOTP1OCN1210815131190.10.50.90.10.50.90.10.50.
13、9EAMPEP(C)竣酸酯组成图3电导率的变量影响图在第一组混料设计中,环状碳酸酯(Ee)和竣酸酯在研究范围内含量增大对电导率都有正效应,但不同温度下两者的正效应大小不同。在温度为10。C时,竣酸酯的正效应大于环状碳酸酯(EC);在温度为25C和45C时,环状碳酸酯(EC)的正效应显著大于较酸酯。且随着温度的升高,环状碳酸酯(EC)的正效应是增大的,而较酸酯的正效应是减小的。对于线性碳酸酯,在3个温度下,研究范围内含量增大对电导率都有负效应。在第二组线性碳酸酯的混料设计中,3个温度下只有DMC在研究范围内含量增大对电导率有正效应,而EMC和DEC研究范围内含量增大对电导率都有负效应。在第三组较
14、酸酯的混料设计中,3个温度下EA和MP在研究范围内含量增大都对电导率有正效应;在温度为-10和25时,MP的正效应大于EA;在温度为45C时,两者的正效应大致相同,总的来说MP对电导率的提升效果略优于EA,而EA则优于对电导率有负效应的EP。由以上分析可知电解液溶剂变量众多且影响复杂,难以直接观察到3个温度下电导率均最大的最优电导率设计。可以将10、25C、45C电导率按照指定权重加和后得到总的电导率函数,对总电导率函数归一化处理使算法收敛加快,进而求解使归一化后总电导率函数最大的各变量的值。本工作将3个温度下电导率权重均设为1,最终得到最优电导率设计的溶剂组成为:环状碳酸酯(EC):0.3、线性碳酸酯:0.3、竣酸酯:0.4;EMC:0、DEC:0、DMC:1;EA:0、MP:1、EP:0,即EC:0.3、DMC:0.3、MP:0.4o如表2所示,实验结果证明了最优电导率设计的预测较为准确。表2最佳电导率设计的模型性能温度/预测值/(mscm)实测值/(mscm)模型误差/(mscm)-1010.0610.43-0.372514.8714.420.454517.3116.950.36在实际应用中电导率的优化只是电解液设计的初