Summary of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide a kind of diagnosis quickly, accurate
The diagnostic method of the battery internal short-circuit of identification.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of diagnostic method of battery internal short-circuit is wrapped to obtain the internal short-circuit resistance value of battery pack internal short-circuit battery cell
Include following steps:
1) charge state difference of each battery cell of internal battery pack is obtained;
2) the difference electricity of each battery cell is calculated according to the charge state difference of battery cell in battery pack;
3) the change rate L of difference electricity is obtained using linear regression to get battery cell is arrived in time of measuring section
Average drain currents Ideplete;
4) according to average drain currents IdepleteAnd the average terminal voltage value of battery pack obtains battery pack in time of measuring section
The near short circuit resistance value of interior each battery cell;
5) by the near short circuit resistance value of each battery cell respectively compared with the short-circuit resistance value threshold value of setting, if near short circuit
Resistance value is greater than short-circuit resistance value threshold value, then determines that the battery cell is normal monomer, if near short circuit resistance value is less than or equal to short circuit
Resistance value threshold value then determines the battery cell for short-circuit monomer.
The step 1) specifically includes the following steps:
11) overall permanence of average cell model battery eliminator group is used, it is average with battery pack current value and battery pack
Voltage value estimates according to extended Kalman filter EKF algorithm the average lotus of battery pack as input value in a manner of high frequency
Electricity condition;
12) using the difference in difference battery model battery eliminator group between single battery characteristic and battery pack overall permanence
It is different, with the difference and average cell model of the average terminal voltage of battery pack current value I, the end voltage of battery cell and battery pack
The average state of charge of estimation is estimated in a manner of low frequency as input value, and according to extended Kalman filter EKF algorithm
The charge state difference of each battery cell.
In the step 2), the calculating formula of difference electricity are as follows:
ΔCk=C Δ SOCk
Wherein, C is the current capacities of battery cell in battery pack, Δ SOCkIt is that k moment battery cell SOC and battery pack are flat
Difference state-of-charge between equal SOC, Δ CkDifference electricity between k moment battery cell electricity and battery pack average electricity
Amount.
In the step 3), the change rate L calculating formula of difference electricity are as follows:
Wherein, Δ C(n)For value of the linear regression straight line in time of measuring segment at rear cut off, Δ C(1)It is linear
The value of regression straight line starting point in time of measuring segment, t(n)For the terminating point of time of measuring segment, t(1)For measurement
The starting point of time interval section.
In the step 4), the calculating formula of the near short circuit resistance value of each battery cell in battery pack are as follows:
Wherein, RISCFor the near short circuit resistance value of battery cell each in battery pack, UmeanIt is electric in time of measuring segment
The average terminal voltage value of pond group,For i-th of the battery pack terminal voltage value acquired in time of measuring segment, n is when measuring
Between the battery pack terminal voltage value sum that acquires in segment.
In the step 5), the short-circuit resistance value threshold value set is 200 Ω.
In the step 11), average cell model is Order RC model.
The expression formula of the average cell model are as follows:
Umean=Uoc(SOCmean)-IR0-UD-UT
Wherein, UmeanFor the end voltage of average cell model, UocFor the voltage source of average cell group model, SOCmeanFor electricity
The average state-of-charge of pond group, I are battery pack current value, R0For internal resistance, UDFor the distribution voltage of activation polarization internal resistance, UTIt is dense
The distribution voltage of poor polarization resistance.
In the step 12), difference battery model is Rint model.
The expression formula of the difference battery model are as follows:
Wherein,For the terminal voltage value of difference battery,For in SOCmeanNear, i-th single battery and average
The open-circuit voltage U of batteryocDifference, Δ SOCiFor the average state-of-charge SOC of i-th of single battery and battery packmeanDifference,
I is battery pack current value, Δ RiIt is averaged the difference internal resistance of internal resistance for i-th internal resistance of single cell and battery pack.
Compared with prior art, the invention has the following advantages that
One, diagnosis is quick: the present invention can quantitatively be diagnosed under limited hardware condition with the calculation amount of very little online
There are the single batteries of internal short-circuit for internal battery pack, and obtain the near short circuit resistance value of internal short-circuit battery cell in battery pack, energy
It is enough quantitatively to be prejudged in advance when internal short-circuit sign is small, the probability that thermal runaway occurs for battery is reduced, the usability of battery pack is improved
Energy.
Two, accurate identification: the present invention uses average cell model and difference battery model, in conjunction with Extended Kalman filter
EKF algorithm realizes the identification of each battery cell charge state difference in battery pack, and then realizes the internal short-circuit failure of battery pack
Diagnosis, improves dynamic property, safety and the durability of battery pack.
Embodiment
The difference that the present invention picks out the state-of-charge SOC of internal battery pack battery cell by dividing models coupling EKF,
The difference state-of-charge Δ SOC obtained using this method is diagnosed to be the near short circuit resistance value of internal short-circuit battery cell, can with compared with
Low calculation amount online in quantitative diagnosis battery pack internal short-circuit battery cell near short circuit resistance value.Internal short-circuit battery list occurs
Know from experience internal self-contained circuit, constantly consumes the electricity of battery cell, in turn resulting in battery cell charge state difference Δ SOC has
The trend constantly increased.
The present invention diagnose battery cell internal short-circuit process as shown in Figure 1, specifically includes the following steps:
The operational data of S1, the battery pack measured by sensor utilize frequency dividing in conjunction with the battery basic parameter of battery pack
Model, including average cell model and difference battery model pick out internal battery pack electricity in conjunction with extended Kalman filter EKF
The difference of the state-of-charge SOC of pond monomer;
Step S1 specifically includes the following steps:
1) battery pack operational data is acquired using current sensor, voltage sensor and temperature sensor;
2) overall permanence using Order RC model as shown in Figure 2 as average cell model battery eliminator group, here
A big battery will be regarded as by single battery battery pack in series, with battery pack current value I that current sensor measures with
And the average voltage level U' that voltage sensor measuresmeanBattery is estimated in a manner of high frequency as input value, and in conjunction with EKF algorithm
The average state-of-charge SOC of groupmean, and by SOCmeanEstimated result be output to controller, controller controls battery pack charge and discharge
Electricity.
In Fig. 2, R0Represent the internal resistance of battery pack;The R being connected in parallelDCDAnd RTCTRespectively represent the activation polarization of battery pack
Internal resistance and concentration polarization internal resistance;UDAnd UTRespectively R in modelDCDAnd RTCTDistribution voltage;UocFor average cell group model
Voltage source represents the open-circuit voltage of battery pack, in the state of the equilibrium, UocWith the state-of-charge SOC of battery packmeanThere is one-to-one correspondence
Relationship;UmeanFor the end voltage of average cell model.I is the size of current of battery pack, because recognized battery pack is by more
A single battery is connected in series, thus by the electric current of battery pack and all monomers be it is equal, measured by current sensor.
The Order RC model of step 2), each parameter relationship formula are as follows:
Umean=Uoc(SOCmean)-IR0-UD-UT (1)
Wherein, SOCmeanIndicate the average state-of-charge of battery pack, i.e., all single batteries is averaged in expression battery pack
SOC;UocRepresent the open-circuit voltage of battery pack, in the state of the equilibrium, UocValue and SOCmeanValue have one-to-one relationship;
I is the size of current of battery pack, and the battery pack of the program is connected in series by multiple single batteries, and battery pack and all lists are passed through
The electric current of body is equal;R0Indicate the ohmic internal resistance of battery pack;UDAnd UTIndicate the distribution voltage of the polarization resistance of battery pack.
UmeanThe end voltage of average cell model, theoretically with monomers all in battery pack end voltage(being measured by voltage sensor)
Average value U'meanIt is equal (error to be normally present, because in order to reduce calculation amount in existing engineering practice, all
Using better simply equivalent-circuit model, generally can not simulated battery completely genuine property, and only approximate equivalent mould
Type.), U'meanAs shown in formula (2):
In formula,The end voltage of i-th of single battery is represented, N indicates the number of monomer in the battery pack.
EKF algorithm, state equation and output equation are respectively (3) and (4) formula:
xk+1=f (xk,uk)+wk (3)
yk=g (xk,uk)+vk (4)
Wherein, f (xk,uk) it is function of state, g (xk,uk) it is measurement functions, xkFor state vector, ukFor input value, step
U in S2kFor the battery pack current value that current sensor measures, ykFor output valve, i.e. model estimate value, y in step S2kFor model
Output voltage, wk, vkIt is mean value is 0, variance random process noise, collective effect determines kalman gain in EKF algorithm
It is worth, here wk, vkVariance take Var (wk)=1e-8, Var (vk)=0.01^2.The size of the two variances indicates external and perturbs
Or parameter error (including model and system parameter error) is worth smaller, table to the influence degree of model and system output result
Show that respective value is with a high credibility, the two parameter values are bigger, indicate that corresponding output result distortion level is big.So general Var
(wk) bigger, Var (vk) smaller, then it represents that it is big that input value is disturbed degree, and it is small that system output reduced value is disturbed degree, then institute
Reduced value is more reliable, and system more believes that actually measured reduced parameter value, corresponding kalman gain value are also bigger.In turn,
Var(wk) smaller, Var (vk) bigger, then input value is more reliable, and system more believes the parameter estimated by model measured, corresponding
Kalman gain value is also smaller.General estimation early period takes Var (wk) larger, Var (vk) smaller, so that estimated value restrains as early as possible.
Parameter matrix is obtained using first order Taylor formula linearization equations (3) and (4)
State vector x in EKF algorithmkIt may be expressed as:
xk=[SOCmean,k,UD,k,UT,k]T (7)
Wherein, SOCmean,kIndicate the average state-of-charge of k node time instance battery pack;UD,kAnd UT,kIndicate k node time instance electricity
The distribution voltage of the polarization resistance of pond group.
g(xk,uk)=Uoc(SOCmean,k)-IkR0-UD,k-UT,k (9)
In formula, f (xk,uk) be average cell function of state;g(xk,uk) be average cell measurement functions, i.e., it is average
The terminal voltage value of the average cell of battery model estimation;Δ t is the sampling time, takes 1s, i.e. average cell model estimates battery pack
The frequency of average SOC is 1HZ.τD, τTRespectively represent the activation polarization internal resistance R of battery packDCDWith concentration polarization internal resistance RTCT;
SOCmean,kFor the average state-of-charge of k node time instance battery pack, i.e., the average SOC of all single batteries in battery pack;UD,kWith
UT,kFor the distribution voltage of the polarization resistance of k node time instance battery pack;IkThe current value measured for k node time instance current sensor;
Uoc(SOCmean,k) indicate k node time instance battery pack open-circuit voltage;η is coulombic efficiency, and when electric discharge takes 1, less than 1 when charging, this
In take 0.99;R0For the internal resistance of average cell.
F (x in formula (8)k,uk) and (9) in g (xk,uk) to state vector xkPartial differential is carried out, and is taken respectivelyWithObtain the parameter matrix of average cell model are as follows:
3) using Rint model as shown in Figure 3 as single battery characteristic in the equivalent battery pack of difference battery model with
Difference between battery pack overall permanence, the battery pack current value I measured with current sensor, the monomer that voltage sensor measures
The difference value of voltage and average voltage levelAs shown in formula (13) and average cell model estimation average state-of-charge
SOCmeanValue is used as input value, and single battery state-of-charge SOC is estimated in a manner of low frequency in conjunction with EKF algorithmiIt is flat with battery pack
Equal state-of-charge SOCmeanBetween discrepancy delta SOCi, and the Δ SOC that difference battery model is estimatediValue is output to controller.
The Rint model of step S3, each parameter relationship formula are as follows:
ΔSOCiIndicate the average state-of-charge SOC of i-th of single battery and battery packmeanDifference;It represents
SOCmeanNear, the open-circuit voltage U of i-th of single battery and average cellocDifference, in the state of the equilibrium,Value and Δ
SOCiValue have one-to-one relationship;ΔRiIndicate that i-th of internal resistance of single cell and battery pack are averaged in the difference of internal resistance
Resistance;For the end voltage of difference battery, theoretically with the end voltage of i-th of single batteryWith the average end electricity of battery pack
Press UmeanDifferenceEqual (it is normally present error, because not accounting for the difference of polarization resistance in differential pattern, and
Consider that the difference of polarization resistance can greatly increase the complexity of calculating, applied model is approximate difference battery model.),
As shown in formula (13):
The EKF algorithm of step 3), state equation and output equation such as formula (3) and formula (4).
State vector may be expressed as: in the EKF algorithm of step 3)
Indicate that i-th of monomer is averaged state-of-charge SOC in k node time instance battery SOC and battery packmeanDifference
It is different.
In formula,For the function of state of difference battery;For the measurement functions of difference battery, i.e. difference battery
The end voltage of the difference battery of model estimation;Indicate the average state-of-charge SOC in equilibrium statemean
Near, i-th of single battery open-circuit voltage and battery pack average open-circuit voltage Uoc(SOCmean) difference;ΔRiIt indicates i-th
The internal resistance of single battery and battery pack are averaged the difference internal resistance of internal resistance;IkThe current value measured for k node time instance current sensor.
F (x in formula (15)k) and (16) in g (xk) to state vector xkPartial differential is carried out, and is taken respectivelyWith
It should be pointed out that in SOCmeanNeighbouring Δ SOC- Δ UocBetween have one-to-one relationship, formula (18) indicate
In SOCmeanNeighbouring Δ UocTo the derivative value of Δ SOC, can also indicate are as follows:
The frequency of difference battery model estimation Δ SOC takes 0.01HZ in step 3).
Difference state-of-charge Δ between S2, the SOC that is averaged using the obtained battery cell SOC of step S1 and battery pack
The difference of battery cell electricity Yu battery pack average electricity is calculated in conjunction with the capacity parameter of battery cell in the estimated value of SOC
Electricity Δ C;
Obtain battery cell difference electricity Δ C's by the estimated value of battery cell difference state-of-charge Δ SOC in step S2
Calculation formula are as follows:
ΔCk=C Δ SOCk
In formula, C indicates the current capacities of battery cell, i.e. total electricity of present battery monomer when fully charged, unit is
Ah, in the case where not considering the variation of the battery capacity as caused by the factors such as durability, temperature, it is believed that C and the battery
The initial capacity of monomer is consistent;ΔSOCkIndicate the difference lotus that k node time instance battery cell SOC and battery pack are averaged between SOC
Electricity condition;ΔCkIndicate the difference electricity between k node time instance battery cell electricity and battery pack average electricity.
S3, by the estimated value of battery cell difference electricity Δ C, obtain the change rate of Δ C using linear regression to get to poor
Average drain currents of the different battery cell within this time;
In step s3, it obtains utilizing least-squares algorithm in the calculation method of the change rate of Δ C using linear regression;
In step S3 the change rate of the Δ C of difference battery cell and within this time the battery cell average drain currents
It is unified numerical value, relationship are as follows:
The change rate of L expression Δ C;ΔC(n)After indicating linear regression straight line in measured time interval section at cut off
Value;ΔC(1)Indicate value of the linear regression straight line in measured time interval section at its real point;t(n)Indicate the measured time
The terminating point of segment;t(1)Indicate the starting point of measured time interval section;IdepleteIndicate difference battery cell measured
Leakage current in time interval section.
S4, using the average drain currents of difference battery cell, be diagnosed to be internal short-circuit in conjunction with the average voltage in this period
The near short circuit resistance value of battery cell;
Estimate that the internal short-circuit resistance value of internal short-circuit battery cell utilizes Ohm's law in step S4:
RISCFor the internal short-circuit resistance value of estimated internal short-circuit battery cell;IdepleteIt is difference battery cell measured
Leakage current in time interval section;UmeanFor the average terminal voltage value of battery pack in measured time interval section, may be expressed as:
For i-th of battery pack terminal voltage value that voltage sensor in time of measuring segment acquires, n is time of measuring
The battery pack terminal voltage value sum that voltage sensor acquires in segment.
Whether S5, the battery cell near short circuit resistance value estimated by step S4 are greater than default short-circuit resistance value, are to go to
Otherwise step S6 goes to step S7;
In this embodiment, the preset value in step S5 may be set to 200 Ω.
S6, battery management system determine that the battery cell is normal monomer;
S7, battery management system determine the battery cell for internal short-circuit monomer.
According to the above process, the present invention picks out the charged of internal battery pack battery cell by dividing models coupling EKF
The difference of state SOC, the difference state-of-charge Δ SOC obtained using this method are diagnosed to be the near short circuit of internal short-circuit battery cell
Resistance value, can with the near short circuit resistance value of internal short-circuit battery cell in lower calculation amount online quantitative diagnosis battery pack,
Internal short-circuit sign is diagnosable when small to be arrived, and then provides reliable reference frame for battery failures diagnosis, reduces battery and heat mistake occurs
The probability of control, and then improve the reliability that battery pack uses.