CN109239602A - 一种动力电池的欧姆内阻的估算方法 - Google Patents

Abstract

本发明公开了一种动力电池的欧姆内阻的估算方法，包括：建立反映欧姆内阻变化的离散化计算公式，建立动力电池的PNGV等效电路模型，计算得到动力电池的开路电压U_OC和欧姆内阻的函数关系式，然后根据反映欧姆内阻变化的离散化计算公式建立动力电池的欧姆内阻的状态方程，根据动力电池的开路电压U_OC和欧姆内阻的函数关系式建立所述动力电池的欧姆内阻的观测方程，并根据粒子滤波算法进行迭代得到所述动力电池的欧姆内阻。本发明提出的动力电池的欧姆内阻的估算方法，能够精确地估算出动力电池的欧姆内阻。

Description

一种动力电池的欧姆内阻的估算方法

技术领域

本发明涉及新能源汽车的电池技术领域，尤其涉及一种动力电池的欧姆内阻的估算方法。

背景技术

锂电子电池由于其高能量比和高功率比成为目前电动汽车的常用电池种类。在纯电动汽车上，电池不止作为汽车的动力来源，同时也为电空调、电制动、电转向三电系统供应能源，决定着电动汽车的续航以及驾乘品质。欧姆内阻主要是由电极材料、电解液、隔膜的电阻及各部分零件的接触电阻组成。电池在使用初期初期的内阻较小，但在长期使用后，由于电池内部电解液的枯竭和电池内部化学物质活性的降低，内阻会逐渐变大，直到电池内部的电量无法被正常释放出来，此时电池寿命终止。当电池内阻达到一定程度的时候，会严重影响到电池的效率，然而电池的欧姆内阻无法用常规方法精确测量，所以对于欧姆内阻的估算研究显得极为重要。

以上背景技术内容的公开仅用于辅助理解本发明的构思及技术方案，其并不必然属于本专利申请的现有技术，在没有明确的证据表明上述内容在本专利申请的申请日已经公开的情况下，上述背景技术不应当用于评价本申请的新颖性和创造性。

发明内容

为了解决上述技术问题，本发明提出一种动力电池的欧姆内阻的估算方法，能够精确地估算出动力电池的欧姆内阻。

为了达到上述目的，本发明采用以下技术方案：

本发明公开了一种动力电池的欧姆内阻的估算方法，包括以下步骤：

S1：建立反映欧姆内阻变化的离散化计算公式：R_0(k)＝R_0(k-1)+r_(k-1)，其中R₀是所述动力电池的欧姆内阻，k表示第k次采样周期，r表示误差干扰矩阵；

S2：建立所述动力电池的PNGV等效电路模型，计算得到所述动力电池的开路电压U_OC和欧姆内阻的函数关系式：U_OC(k)＝f(R₀)；

S3：根据步骤S1中的反映欧姆内阻变化的离散化计算公式建立所述动力电池的欧姆内阻的状态方程：x_k＝x_k-1+r_k-1+w_k，根据步骤S2中的所述动力电池的开路电压U_OC和欧姆内阻的函数关系式建立所述动力电池的欧姆内阻的观测方程：z_k＝f(x_k)+v_k，并根据粒子滤波算法进行迭代得到所述动力电池的欧姆内阻，其中，x_k是第k次采样周期的欧姆内阻，z_k是第k次采样周期的系统观测值，w_k为系统噪声，v_k为观测噪声。

优选地，步骤S2具体包括：

S21：建立所述动力电池的PNGV等效电路模型，数学表达式为：

其中，U_OC是所述动力电池的开路电压，C_b是所述动力电池的等效电容，I是所述动力电池的负载电流，U是所述动力电池的端电压，R_P是所述动力电池的极化内阻，C_p是所述动力电池内部的极化电容，U_P是流经极化内阻的电压；

S22：根据所述动力电池的PNGV等效电路模型的数学表达式，得到对应的离散状态空间整理方程为：

U(k)＝a₁U(k-1)+a₂U(k-2)+a₃I(k)+a₄I(k-1)+a₅

其中，

θ＝1/(C_pR_p)，T为采样周期；

S23：测试得到所述动力电池的U(k-2)、U(k-1)、U(k)、I(k-1)、I(k)，并代入到步骤S22中的离散状态空间整理方程，计算得到不同温度及SOC(State of Charge，蓄电池剩余容量)下的R₀、R_P、C_b、C_p；

S24：将步骤S23计算得到的不同温度及SOC下的R₀、R_P、C_b、C_p输入到步骤S21中的所述动力电池的PNGV等效电路模型的数学表达式，拟合得到所述动力电池的开路电压U_OC和欧姆内阻的函数关系式：U_OC(k)＝f(SOC)。

优选地，步骤S22中根据所述动力电池的PNGV等效电路模型的数学表达式，得到对应的离散状态空间整理方程具体包括：

S221：以所述动力电池的PNGV等效电路模型的两电容电压值为状态变量，以所述动力电池的端电压为输出变量，将步骤S21中的数学表达式转换为对应的空间状态方程为：

S222：将步骤S221中的空间状态方程转换为离散状态空间方程为：

S223：将步骤S222中的离散状态空间方程进行整理得到对应的离散状态空间整理方程为：U(k)＝a₁U(k-1)+a₂U(k-2)+a₃I(k)+a₄I(k-1)+a₅。

优选地，步骤S23中采用HPPC测试方法测试得到所述动力电池的U(k-2)、U(k-1)、U(k)、I(k-1)、I(k)。

优选地，步骤S3中根据粒子滤波算法进行迭代得到所述动力电池的欧姆内阻具体包括：

S31：初值化：在先验概率密度p(x₀)中分布采样产生粒子群将每个粒子的权值统一设为i代表是第i个粒子，N_s为粒子总数；

S32：更新：第k次采样周期的粒子x_k的后验概率分布p(x_k|z_k)用带权样本来描述，是粒子集的对应权重，带权样本根据重要密度函数q(x_0:k-1|z_1:k)获得；

S34：预测估计：根据步骤S32得到状态估计公式：方差估计公式：并根据所述动力电池的欧姆内阻的状态方程预测下一次采样周期的欧姆内阻：

S35：检测k是否迭代完成，如果否，则令k＝k+1，并返回步骤S32，如果是，则结束迭代，得到所述动力电池的欧姆内阻。

优选地，步骤S3还包括：

S33：重采样：在步骤S32之后，判断是否小于阈值N_threshold，如果是，则进行重采样，将带权样本映射成权值相同的新样本

优选地，步骤S32中带权样本根据重要密度函数q(x_0:k-1|z_1:k)获得的具体步骤包括：

S321：在第k次采样周期中，采样粒子集由重要密度函数得到：

S322：计算重要性权值：

将后验概率分解为：

将重要密度函数分解为：q(x_0:k|z_1:k)＝q(x_k|x_0:k-1|z_1:k)q(x_0:k-1|z_1:k-1)，

得到重要性权值更新公式为：

S323：对重要性权值进行归一化：

S324：得到x_k的最优估计值为：

与现有技术相比，本发明的有益效果在于：本发明提出的动力电池的欧姆内阻的估算方法采用粒子滤波算法进行迭代估算，其中状态方程根据反映欧姆内阻变化的离散化计算公式来建立，观测方程采用通过PNGV等效电路模型拟合出的动力电池的开路电压和欧姆内阻的函数关系式；将粒子滤波和PNGV等效电路模型进行结合，既利用了粒子滤波算法提高算法的精准度和在非线性非高斯系统强表达能力，也吸收了PNGV模型阶数较低计算相对容易，且物理意义清晰，模型参数辨别简单的优势，达到精确估算欧姆内阻的目的。

附图说明

图1是本发明优选实施例的动力电池的欧姆内阻估算方法的流程示意图；

图2是本发明优选实施例的动力电池的PNGV等效电路模型示意图。

具体实施方式

下面对照附图并结合优选的实施方式对本发明作进一步说明。

如图1所示，本发明优选实施例提出一种结合粒子滤波算法和PNGV等效电路模型的动力电池的欧姆内阻的估算方法，首先根据PNGV模型建立等效电路，然后建立动力电池的欧姆内阻的状态方程和观测方程，再用粒子滤波算法进行迭代，算出更为精准的动力电池的欧姆内阻估算值；该动力电池的欧姆内阻的估算方法包括以下步骤：

S1：建立反映欧姆内阻变化的离散化计算公式：R_0(k)＝R_0(k-1)+r_(k-1)，其中R₀是所述动力电池的欧姆内阻，k表示第k次采样周期，r表示误差干扰矩阵；

其中电池的充放电效率可以通过充放电试验测得。

S2：建立动力电池的PNGV等效电路模型，计算得到动力电池的开路电压U_OC和欧姆内阻的函数关系式：U_OC(k)＝f(R₀)，其中k表示第k次采样周期；

其中，步骤S2具体包括：

S21：如图2所示，建立动力电池的PNGV等效电路模型，数学表达式为：

其中，U_OC是动力电池的开路电压，C_b是动力电池的等效电容，描述动力电池开路电压随负载电流积分变化，R₀是动力电池的欧姆内阻，I是动力电池的负载电流，放电为正，充电为负，U是动力电池的端电压，R_P是动力电池内部由于电解液浓度差造成的极化内阻，C_p是动力电池内部由于电解液浓度差造成的极化电容，U_P是流经极化内阻的电压；

S22：根据动力电池的PNGV等效电路模型的数学表达式，得到对应的离散状态空间整理方程，具体包括：

S221：以动力电池的PNGV等效电路模型的两电容电压值为状态变量，以动力电池的端电压为输出变量，将步骤S21中的数学表达式转换为对应的空间状态方程为：

S222：将步骤S221中的空间状态方程转换为离散状态空间方程为：

S223：将步骤S222中的离散状态空间方程进行整理得到对应的离散状态空间整理方程为：

U(k)＝a₁U(k-1)+a₂U(k-2)+a₃I(k)+a₄I(k-1)+a₅，

其中，

θ＝1/(C_pR_p)，T为采样周期，k表示第k次采样周期；

S23：测试得到动力电池的U(k-2)、U(k-1)、U(k)、I(k-1)、I(k)，并代入到步骤S22中的离散状态空间整理方程，计算得到不同温度及SOC下的R₀、R_P、C_b、C_p；

在本实施例中，采用HPPC测试方法测试得到动力电池的U(k-2)、U(k-1)、U(k)、I(k-1)、I(k)；具体地，在常温下恒定电流100A(1C)放电，在不同的SOC点进行复合脉冲实验(对于100AH的电池容量，100A电流放电6分钟，可放出10AH，即SOC下降0.1)，放电后静置1h，进入下一个HPPC循环，放电前测量电池电压，以得到对应SOC的电池开路电压近似值。

S24：将步骤S23计算得到的不同温度及SOC下的R₀、R_P、C_b、C_p输入到步骤S21中的动力电池的PNGV等效电路模型的数学表达式，拟合得到动力电池的开路电压U_OC和欧姆内阻的曲线(函数关系式)：U_OC(k)＝f(R₀)。

S3：根据步骤S1中的反映欧姆内阻变化的离散化计算公式建立动力电池的欧姆内阻的状态方程：x_k＝x_k-1+r_k-1+w_k，根据步骤S2中的动力电池的开路电压U_OC和欧姆内阻的函数关系式建立动力电池的欧姆内阻的观测方程：z_k＝f(x_k)+v_k，并根据粒子滤波算法进行迭代得到动力电池的欧姆内阻，其中，x_k是第k次采样周期的状态值(欧姆内阻)，z_k是第k次采样周期的系统观测值，w_k为系统噪声，v_k为观测噪声，w_k和v_k是相互独立的随机噪声。

其中，根据粒子滤波算法进行迭代得到动力电池的欧姆内阻具体包括：

S31：初值化：在先验概率密度p(x₀)中分布采样产生粒子群将每个粒子的权值统一设为i代表是第i个粒子，N_s为粒子总数；

S32：更新：第k次采样周期的粒子x_k的后验概率分布p(x_k|z_k)用带权样本来描述，是粒子集的对应权重，带权样本根据重要密度函数q(x_0:k-1|z_1:k)获得；

在本实施例中，带权样本根据重要密度函数q(x_0:k-1|z_1:k)获得具体包括：

S321：在第k次采样周期中，采样粒子集由重要密度函数得到：

S322：计算重要性权值：

将后验概率分解为：

将重要密度函数分解为：q(x_0:k|z_1:k)＝q(x_k|x_0:k-1|z_1:k)q(x_0:k-1|z_1:k-1)，

得到重要性权值更新公式为：

S323：对重要性权值进行归一化：

S324：得到x_k的最优估计值为：

S33：重采样：判断是否小于阈值N_threshold，如果是(也即则进行重采样，将带权样本映射成权值相同的新样本如果否则直接进入步骤S34；通过重采样步骤，可防止粒子退化。

S34：预测估计：根据步骤S32得到状态估计公式：方差估计公式：(该公式中的上标T为向量转置符合)，并根据动力电池的欧姆内阻的状态方程预测下一次采样周期的状态值(欧姆内阻)：

S35：检测k是否迭代完成，如果否，则令k＝k+1，并返回步骤S32，如果是，则结束迭代，得到动力电池的欧姆内阻。

本发明优选实施例提出的动力电池的欧姆内阻的估算方法采用粒子滤波算法进行迭代估算，其中状态方程采用反映欧姆内阻变化的离散化计算公式，观测方程采用通过PNGV等效电路模型拟合出的动力电池的开路电压和欧姆内阻的函数关系式；将粒子滤波和PNGV等效电路模型进行结合，既利用了粒子滤波算法提高算法的精准度和在非线性非高斯系统强表达能力，也吸收了PNGV模型阶数较低计算相对容易，且物理意义清晰，模型参数辨别简单的优势，达到精确估算欧姆内阻的目的。其中如果将PNGV等效电路模型的等效方程直接作为粒子滤波算法的观测方程的话，模型会比较复杂，而且PNGV等效电路模型的参数较多容易引起观测系统不稳定，而在本发明优选实施例中，首先通过PNGV等效电路模型估算欧姆内阻，得到动力电池的开路电压与欧姆内阻的拟合曲线，简化观测维度增强稳定性，再利用粒子滤波算法的强表达能力提高估算精准性。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干等同替代或明显变型，而且性能或用途相同，都应当视为属于本发明的保护范围。

Claims (7)

1.一种动力电池的欧姆内阻的估算方法，其特征在于，包括以下步骤：

S1：建立反映欧姆内阻变化的离散化计算公式：R_0(k)＝R_0(k-1)+r_(k-1)，其中R₀是所述动力电池的欧姆内阻，k表示第k次采样周期，r表示误差干扰矩阵；

S2：建立所述动力电池的PNGV等效电路模型，计算得到所述动力电池的开路电压U_OC和欧姆内阻的函数关系式：U_OC(k)＝f(R₀)；

S3：根据步骤S1中的反映欧姆内阻变化的离散化计算公式建立所述动力电池的欧姆内阻的状态方程：x_k＝x_k-1+r_k-1+w_k，根据步骤S2中的所述动力电池的开路电压U_OC和欧姆内阻的函数关系式建立所述动力电池的欧姆内阻的观测方程：z_k＝f(x_k)+v_k，并根据粒子滤波算法进行迭代得到所述动力电池的欧姆内阻，其中，x_k是第k次采样周期的欧姆内阻，z_k是第k次采样周期的系统观测值，w_k为系统噪声，v_k为观测噪声。

2.根据权利要求1所述的动力电池的欧姆内阻的估算方法，其特征在于，步骤S2具体包括：

S21：建立所述动力电池的PNGV等效电路模型，数学表达式为：

其中，U_OC是所述动力电池的开路电压，C_b是所述动力电池的等效电容，I是所述动力电池的负载电流，U是所述动力电池的端电压，R_P是所述动力电池的极化内阻，C_p是所述动力电池内部的极化电容，U_P是流经极化内阻的电压；

S22：根据所述动力电池的PNGV等效电路模型的数学表达式，得到对应的离散状态空间整理方程为：

U(k)＝a₁U(k-1)+a₂U(k-2)+a₃I(k)+a₄I(k-1)+a₅

其中，

θ＝1/(C_pR_p)，T为采样周期；

S23：测试得到所述动力电池的U(k-2)、U(k-1)、U(k)、I(k-1)、I(k)，并代入到步骤S22中的离散状态空间整理方程，计算得到不同温度及SOC下的R₀、R_P、C_b、C_p；

S24：将步骤S23计算得到的不同温度及SOC下的R₀、R_P、C_b、C_p输入到步骤S21中的所述动力电池的PNGV等效电路模型的数学表达式，拟合得到所述动力电池的开路电压U_OC和欧姆内阻的函数关系式：U_OC(k)＝f(R₀)。

3.根据权利要求2所述的动力电池的欧姆内阻的估算方法，其特征在于，步骤S22中根据所述动力电池的PNGV等效电路模型的数学表达式，得到对应的离散状态空间整理方程具体包括：

S221：以所述动力电池的PNGV等效电路模型的两电容电压值为状态变量，以所述动力电池的端电压为输出变量，将步骤S21中的数学表达式转换为对应的空间状态方程为：

S222：将步骤S221中的空间状态方程转换为离散状态空间方程为：

S223：将步骤S222中的离散状态空间方程进行整理得到对应的离散状态空间整理方程为：U(k)＝a₁U(k-1)+a₂U(k-2)+a₃I(k)+a₄I(k-1)+a₅。

4.根据权利要求2所述的动力电池的欧姆内阻的估算方法，其特征在于，步骤S23中采用HPPC测试方法测试得到所述动力电池的U(k-2)、U(k-1)、U(k)、I(k-1)、I(k)。

5.根据权利要求1所述的动力电池的欧姆内阻的估算方法，其特征在于，步骤S3中根据粒子滤波算法进行迭代得到所述动力电池的欧姆内阻具体包括：

S31：初值化：在先验概率密度p(x₀)中分布采样产生粒子群将每个粒子的权值统一设为i代表是第i个粒子，N_s为粒子总数；

S32：更新：第k次采样周期的粒子x_k的后验概率分布p(x_k|z_k)用带权样本来描述，是粒子集的对应权重，带权样本根据重要密度函数q(x_0:k-1|z_1:k)获得；

S34：预测估计：根据步骤S32得到状态估计公式：方差估计公式：并根据所述动力电池的欧姆内阻的状态方程预测下一次采样周期的欧姆内阻：

S35：检测k是否迭代完成，如果否，则令k＝k+1，并返回步骤S32，如果是，则结束迭代，得到所述动力电池的欧姆内阻。

6.根据权利要求5所述的动力电池的欧姆内阻的估算方法，其特征在于，步骤S3还包括：

S33：重采样：在步骤S32之后，判断是否小于阈值N_threshold，如果是，则进行重采样，将带权样本映射成权值相同的新样本

7.根据权利要求5所述的动力电池的欧姆内阻的估算方法，其特征在于，步骤S32中带权样本根据重要密度函数q(x_0:k-1|z_1:k)获得的具体步骤包括：

S321：在第k次采样周期中，采样粒子集由重要密度函数得到：

S322：计算重要性权值：

将后验概率分解为：

将重要密度函数分解为：q(x_0:k|z_1:k)＝q(x_k|x_0:k-1|z_1:k)q(x_0:k-1|z_1:k-1)，

得到重要性权值更新公式为：

S323：对重要性权值进行归一化：

S324：得到x_k的最优估计值为：