[发明专利]一种动力电池的峰值功率预测方法

权利要求书

1.一种动力电池的峰值功率预测方法，其特征是：它由以下步骤实现：

步骤A、电池模型参数在线估计的步骤，具体为：

步骤A1、在对二次电池建模时，由于电动汽车运行工况频率特性，因此电池电化学阻抗谱模型中的中频率的阻抗特性由常用的纯阻性元件R和常相位元件Q并联电路简化为纯阻性元件R来描述，得到简化后的电池电化学阻抗谱等效电路模型；

该简化后的电化学阻抗谱等效电路模型包括开路电压OCV_e、欧姆内阻R_o和韦伯阻抗Z_W；

步骤A2、根据步骤A1获得的简化后的电化学阻抗谱等效电路模型建立分数阶卡尔曼滤波器所需的状态方程与观测方程，具体为：

取流经二次电池的总电流I_L在放电时为正值，数据采样周期为1s；

Δr=drdtr,r>0]]>

其中△^r为微分算子，r为微分阶数，当r为小数时，△^r表示分数阶微分算子，当r为整数时，△^r为整数微分算子；

取分数阶元件Z_W是两端电压为U_W的状态量，有：

Δ0.5UW=1WIL=XWIL]]>

对于电池模型参数，扩散参数X_W、开路电压OCV_e和欧姆内阻R_o随着电池荷电状态(SOC)的变化是缓慢的，因此：

Δ1XW≈0Δ1OCVe≈0Δ1Ro≈0]]>

将上述四个方程改写为矩阵形式，获得分数阶联合卡尔曼滤波器的状态方程：

Δ0.5111UWXWOCVeRo=0IL00000000000000UWXWOCVeRo;]]>

取U_L为系统的观测量，则有：

U_L＝OCV_e-I_LR_o-U_W

I_L表示与流经电池的总电流；

取：

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=diag[0.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]]>

步骤A3、利用步骤A2构建的分数阶卡尔曼滤波器所需的状态方程与观测方程，对状态、参数及协方差矩阵按照分数阶联合卡尔曼滤波算法进行时间更新和测量更新：

具体为：

初始化：

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]]>

其中，Q_k是噪声w_k的协方差，为k时刻状态和模型参数x_k的预测值，为k-1时刻状态和模型参数x_k-1的修正值，为k时刻x的噪声协方差矩阵P_k的预测值，为k-1时刻x的噪声协方差矩阵P_k-1的修正值；

状态、参数及协方差矩阵的测量更新：

Lk=Pk-(Ck)T[CkPk-(Ck)T+Rk]-1]]>

x^k+=x^k-+Lkx[yk-g(x^k-,IL,k)]]]>

Pk+=(I-LkCk)Pk-]]>

其中，R_k是噪声v_k的协方差，L_k是k时刻卡尔曼滤波器增益大小；

步骤A4、采集电池的端电压U_L和流经二次电池的总电流I_L，利用步骤A1获得的简化后的电化学阻抗谱等效电路模型和步骤A3更新后的分数阶卡尔曼滤波器所需的状态方程与观测方程，得到开路电压OCV_e、欧姆内阻R_o、扩散参数X_W的估计值，将获得的开路电压OCV_e、欧姆内阻R_o、扩散参数X_W的估计值作为的电池的估计结果，完成基于分数阶联合卡尔曼滤波的二次电池简化阻抗谱模型参数在线估计；

步骤B、根据步骤A获得的电池模型参数在线估计结果，进行功率预测的步骤：

步骤B1、推导分数阶元件的韦伯阻抗Z_W是否为线性时不变元件：

设分数阶元件的两端电压U_W的初始值为0，当施加一个幅值为I_L的阶跃电流激励时，分数阶元件的两端电压U_W在Ns后的电压响应情况计算如下，N为正数：

设电池模型的开路电压OCV_e、欧姆内阻R_o、扩散参数X_W在功率预测过程中数值大小是不变的；

则U_W在k＝1s时的电压响应为：

U^W,1=X^W,0IL-(-1)10.51UW,0]]>

其中，符号^代表预测值；

U_W在k≥2时的电压响应为：

U^W,k+1=X^W,0IL-Σj=1L(-1)j0.5jU^W,k+1-j+w]]>

计算得：

U^W,1=X^WIL]]>

U^W,2=X^W,0IL-((-1)10.51U^W,1+(-1)20.52UW,0)]]>

计算后得：

U^W,2=1.5X^W,0IL]]>

由此递推得到10s和60s处的电压响应为：

U^W,10=3.524X^W,0IL]]>

U^W,60=8.722X^W,0IL]]>

由此，当电池模型参数X_W不变或者缓慢变化时，输出值U_W与I_L为线性关系，推断分数阶元件Z_W为线性时不变元件；

则，分数阶元件的功率预测方法为：

为了预测分数阶元件在k+ΔTs处的电压响应将该电压响应分为零状态响应和零输入响应

U^W,k+ΔT=U^W,k+ΔTzs+U^W,k+ΔTzi]]>

其中零状态响应为：

U^W,k+ΔTzs=aX^W,kImax]]>

对于k+10s时刻，a＝3.524；对于k+60s时刻，a＝8.722；

零输入响应由k时刻之前的数据决定，取分数阶微分的时间记忆长度L＝60；

k+1时刻的零输入响应为：

U^W,k+1zi=-Σj=1L(-1)j0.5jU^W,k+1-jzi]]>

由此递推可得到k+ΔT时刻分数阶元件处零输入电压响应：

U^W,k+ΔTzi=Σj=1Lbk+1-jU^W,k+1-j]]>

其中，是电池模型分数阶元件Z_W两端端电压在(k+1-L)～k时刻之间的估计值，当预测k+10时刻时，b＝α，α为一组常系数矩阵；

当预测k+60时刻时，b＝β，β为另外一组常系数矩阵；

由此，可得到电池恒流放电电流为I_max时的电池在k+ΔT时刻的端电压预测值：

放电峰值功率预测的方法为：

若电池荷电状态是电池极限工作状态的限制条件，SoC_min是电池放电终止荷电状态的最小值，则得到此时的最大放电电流值：

Imax,SoC=(SoCk-SoCmin)×CapacityΔT]]>

Capacity为电池容量值，单位为安培*小时(Ah)，SoC_k为k时刻的电池荷电状态，SoC_min为电池放电限定的最小荷电状态；

若电池端电压U_L为电池极限工作状态的限制条件，若此时以最大放电电流I_max对电池放电，电池模型参数OCV_e在k+△T时刻的预测值为：

上式中，为通过分数阶联合卡尔曼滤波算法计算得到的电池模型参数OCVe在k时刻的估计值，为OCV_e在k+△T时刻的预测值；I_max是在端电压作为电池极限工作状态限制条件时的最大放电电流值，是待求解值，是电池的最大放电电流；OCV_k是在k时刻的电池开路电压值；为电池在k+△T时刻，假定以恒流放电时对应的电池开路电压值；由于多数情况下，OCV的变化都较小并且缓慢，因此可认为在△T时段内的开路电压变化量与放电电流成线性关系；

进而推算，假设电池在放电电流为I_max时，电池端电压的推算值U_L,k+△T：

其中，端电压的推算值是k时刻开路电压估计值、△T时段内开路电压变化值、欧姆内阻电压差、分数阶元件Z_W的零状态响应电压值、零输入响应电压值之和；

由此推出，当以端电压U_L作为电池的极限工作状态的约束条件时，电池的最大工作电流为：

k+ΔT时刻达到峰值放电功率时，综合考虑上述限制条件，最大放电电流值为：

Imax=min(Imax,SoC,Imax,UL,Imaxlim)]]>

k+ΔT时刻的峰值放电功率为：

电池的回馈电流峰值功率的预测方法与上述放电功率的预测方法同理：

设I_L为正值；

若SoC为电池极限工作状态的限制条件，SoC_max是电池的最大荷电状态值，则得到此时的最小回馈电流值：

Imin,SoC=(SoCk-SoCmax)×CapacityΔT]]>

若U_L为电池极限工作状态的限制条件，则得到此时的最小回馈电流值：

k+ΔT时刻达到峰值回馈功率时，综合考虑上述限制条件，最小回馈电流值为：

Imin=max(Imin,SoC,Imin,UL,Iminlim)]]>

由此，得到电池电流回馈的峰值功率：

完成动力电池的峰值功率预测。

2.根据权利要求1所述的一种动力电池的峰值功率预测方法，其特征在于步骤A中，

由于方法中使用的电池模型基于电池的电化学阻抗谱测试数据，电池模型参数有明确的物理意义，欧姆内阻R_o的物理意义为：

R_o≈R_Ω+R_SEI+R_ct

其中，R_Ω为高频欧姆阻抗，R_SEI为SEI膜阻抗，R_ct为电荷转移阻抗；

另外，韦伯阻抗由以下公式来定义：

ZW=1W(jw)0.5]]>

其中，W为离子扩散系数，为了便于阻抗参数在线估计，取：

XW=1W]]>

得到：

ZW=XW(jw)0.5.]]>