A method of quickly obtaining all cell capacities in battery pack
Technical field
The present invention relates to a kind of methods for quickly obtaining all cell capacities in battery pack.
Background technique
The present invention is related with the capacity of single batteries all in battery pack, is especially connected in series institute's structure by multiple single batteries
At battery pack in all single batteries capacity.Lithium ion battery is because it is high with energy, cell voltage is high, operating temperature
The advantages that range is wide, storage life is long obtains in fields such as electric car field, electrochemical energy storage field, 3C electronic products
To being widely applied.The available capacity of battery is related with the battery last working time, battery resistance value and battery moment charge and discharge electric energy
Power is closely related, in addition in groups of battery pack of connecting, the close phase of consistency of the capacity distribution and battery pack of single battery
It closes.
Capacity of lithium ion battery is one of the important performance indexes for measuring battery performance, it indicates (to put under certain condition
Electric rate, temperature, final voltage etc.) electricity (discharge test can be done with JS-150D) released of battery, i.e. the capacity of battery, usually
It (referred to as, is indicated as unit of ampere-hour with AH, 1Ah=3600C).Lithium ion battery in actual use,
Substantially there is no batteries to be completely full of the case where putting, it is difficult to the total capacity of battery is known, to influence battery and nuclear power state
(SOC) estimation, the accurate real-time calculating to battery capacity, for battery safety management, there are significances.
Chinese invention patent (the patent No.: CN109342955A, a kind of patent name: reckoning side of capacity of lithium ion battery
Method and system), which obtains the lithium battery right by the way that lithium ion battery is carried out discharge test under n experimental temperature
Discharge capacity under experimental temperature establishes model in conjunction with Arrhenius formula, and the battery extrapolated under arbitrary temp point holds
Magnitude.This method needs to carry out n discharge test in advance, complicated for operation, and this method can only extrapolate the battery in not equality of temperature
Under degree for capability value, accurate calculation cannot be carried out to the capability value after cell degradation, practicability is poor.
Chinese invention patent (the patent No.: 108152747 A of CN, patent name: the detection method and dress of accumulator capacity
Set), the patent is by utilizing the battery discharge parameter in storage battery power supply pitch process in acquisition wind power generating set
Measured value;The functional relation using discharge parameter building is tested according to the measured value and by battery, detects institute
State battery current capacities.This method needs constructing measured value and battery capacity using a large amount of battery test of preceding progress
Relational database the measured value of needs is calculated, by looking into actual use by grabbing suitable discharge parameter
Corresponding capacity in ready database in advance is taken, battery current capacities are evaluated.This method needs opening relationships in advance
Database, practice trouble are expensive;And due to the nonlinear characteristic of battery, in implementation process, accumulator capacity
Whether also set up with the corresponding relationship of corresponding measured value, accuracy as a result lacks data and supports.
Chinese invention patent (the patent No.: CN108732508, a kind of patent name: real-time estimation of capacity of lithium ion battery
Method), this method is tested by cell durability, obtains charge data --- the real time capacity VS voltage in normal capacity test
Then curve carries out data to curve and fits within carry out difference.Peak position in difference result curve and size and battery
Residual capacity have correlation, comparison initial situation and certain moment in the case of difference curves, it can for battery capacity into
Row estimation.This method needs to carry out durability test to battery, and the actual operation of this method is poor;The peak of difference result curve
Value size is related with data screening processing method, and error is larger.
Chinese patent (the patent No.: CN109031153, a kind of patent name: health status On-line Estimation of lithium ion battery
Method), this method is tested by short-term cycle life, as the initial parameter of model training, using capacity increment analytic approach,
Multiple characteristic parameters are extracted from incremental capacity plot, form set of characteristic parameters, the health status of battery is portrayed, and
It is exported the value of characteristic parameter as the model of multi output Gaussian process regression model, completes commenting online to battery capacity
Estimate.This method needs to carry out cycle life test, and there are difficulty for practice.
Summary of the invention
The object of the present invention is to provide a kind of charging curves based in normal battery operation, and characteristic value ginseng is extracted after processing
Number, brings model into and is calculated, and completes the capability value for calculating all single batteries in battery pack in real time.
Object of the invention can be realized by the following technical scheme:
A method of all cell capacities in battery pack quickly being obtained, it includes following steps:
S1, it is directed to target lithium ion battery, obtains battery datum curve data, battery nominal capacity Capinitial;
S2, lithium ion battery datum curve data are handled, the capacity or charging for recording characteristic value position cells hold
Amount;
S3, the coefficient of relationship for calculating capacity and total capacity between battery characteristics value;
S4, battery pack charging curve is handled, records the charging capacity of characteristic value position cells;
S5, capacity between the battery characteristics value that S4 step obtains, and the coefficient of relationship according to S3 step are calculated, by
One calculates the capacity of all single batteries in battery pack.
Wherein, battery datum curve data described in the step S1 can be from producer and obtain battery SOC-OCV song
Line number evidence is also possible to history charging curve data.
Wherein, lithium ion battery datum curve data are handled described in the step S2, records characteristic value position
The capacity or charging capacity of battery: it according to battery material type, takes SOC to be greater than 20% data ferric phosphate lithium cell, asks
Incremental capacity plot is taken out, calculation formula isCharacteristic value is extracted, characteristic value position is that capacity increases
The maximum position in curve is measured, and records characteristic value position corresponding capacity Q or charging capacity Q ';For ternary material electricity
Pond takes data of the SOC greater than 20% to seek d2Q/dV2Curve, calculation formula areExtract characteristic value, feature
Value position isPosition, and record characteristic value position corresponding capacity Q or charging capacity Q ';Q=SOC*
Capinitial;
Wherein: Q is the capacity of battery, and Q ' is charging capacity, and dQ is the differential of capacity, d2Q is the second-order differential of capacity, Δ
QkThe difference of capacity between neighbouring sample point, V are the voltage of battery, and dV is the differential of voltage, dV2For the second-order differential of voltage, Δ
VkThe difference of voltage between neighbouring sample point, for each sampled point k, Δ Qk=Qk-Qk-1, Δ Vk=Vk-Vk-1, Δ Qk-1=
Qk-1-Qk-2, Δ Vk-1=Vk-1-Vk-2。
Wherein, the characteristic value of the battery datum curve has 2, i.e. the First Eigenvalue and Second Eigenvalue, records the 1st
The capacity Q of characteristic value position1Or charging capacity Q '1, the capacity Q of the 2nd characteristic value position2Or charging capacity Q '2。
Wherein, the coefficient of relationship g=Δ Q/ of capacity and total capacity between battery characteristics value is calculated in the step S3
Capinitial;Wherein Δ Q is the 1st characteristic value to the capacity difference between the 2nd characteristic value, Δ Q=| Q2-Q1|;If using going through
History charging curve data, then Δ Q=| Q '2-Q‘1|。
Wherein, the step S4 handles battery pack charging curve, record characteristic value position cells charging capacity tool
Body are as follows: extract the data for meeting Δ V >=X condition in battery pack charging process.According to battery material type, for LiFePO4
Battery takes data of the SOC greater than 20% using capacity increment analytic approach, by extracting characteristic value to incremental capacity plot, and records
The corresponding charging capacity in characteristic value position;For ternary material battery, data of the SOC greater than 20% is taken to seek d2Q/dV2Curve,
Characteristic value is extracted, and records the corresponding charging capacity in characteristic value position;The corresponding charging capacity of characteristic value includes the 1st characteristic value
The charging capacity Q ' of position1jWith the charging capacity Q ' of the 2nd characteristic value position2j, in which: Q1jFor the 1st characteristic value position jth
Number cell capacity, Q2jFor the 2nd characteristic value position jth cell capacity, Q '1jFor the 1st characteristic value position jth
Number monomer charging capacity, Q '2jFor the 2nd characteristic value position jth monomer charging capacity.
Wherein, the S5 step calculates capacity between the battery characteristics value that S4 step obtains, and according to S3 step
Coefficient of relationship, calculate the capacity of all single batteries in battery pack one by one specifically: calculate jth single battery the 1st
Characteristic value is to the capacity difference between the 2nd characteristic value, Δ Qj=| Q2j-Q1j|=| Q '2j-Q‘1j|, further according to formula Capj=Δ
Qj/ g obtains jth cell capacity, then calculates other monomers battery capacity in battery pack one by one;Wherein, CapjIt is
J cell capacity, g are the coefficient of relationship of capacity and total capacity between obtained battery characteristics value in claim 5.
Wherein, the value range of the X is 1mV≤X≤10mV.
Wherein, the battery pack charging curve can be the last charging curve, is also possible to nearest 10 times and fills
Charging curve is extracted in electric curve to be combined with that time recently, is become one and is included 2 characteristic value point curves.
Wherein, charging curve is extracted in nearest 10 charging curves to be combined with that time recently, is chosen first
In nearest 10 charging curves as worst curve charging current and environment temperature differs curve less than 5 degree, in residue
The curve comprising two characteristic value position curves can be combined into that nearest secondary curve by choosing again in curve, if there is multiple songs
Line is met the requirements, then preferential selection environment temperature is immediate;If environment temperature is the same, it is bent preferentially to choose newest charging
The upper curve closest to that last secondary charging of line, i.e. time.
The beneficial effects of the present invention are: 1. the method for the present invention are applicable in battery pack normal charging process, do not influence electricity
Pond work input and output;2. needing to get battery SOC-OCV curve data or history charging curve, a battery in advance
Nominal capacity does not need additionally to test battery parameter;3. the total capacity of all single batteries can all obtain in real time in battery pack;4.
It does not need battery and carries out the full total capacity full of all single batteries of putting that you can get it.
Detailed description of the invention
Fig. 1 is flow chart of the invention;
Fig. 2 is the incremental capacity plot of ferric phosphate lithium ion battery datum curve;
Fig. 3 is the incremental capacity plot of ferric phosphate lithium ion battery charge data;
Fig. 4 is the incremental capacity plot and d of ternary lithium ion battery datum curve2Q/dV2Curve;
Fig. 5 is the incremental capacity plot and d of ternary lithium ion battery charge data2Q/dV2Curve;
Fig. 6 is all cell capacity histograms in battery pack;
Fig. 7 is the Error Graph of the present invention calculated capacity and true capacity.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
As shown in Figure 1, a kind of method for quickly obtaining all cell capacities in battery pack, the middle battery of battery pack can
To be ferric phosphate lithium cell or ternary material battery;Battery pack can be by multiple battery cores it is first in parallel after concatenated battery pack system again
System.For target lithium ion battery, battery datum curve data, battery nominal capacity Cap are obtainedinitial.The lithium ion battery
Datum curve data can be from producer and obtain battery SOC-OCV curve data, be also possible to history charging curve data.From factory
Family obtains battery SOC-OCV curve data, and SOC value density includes 0.1% and 5% between 0.1%~5%.To curve number
According to being handled, according to battery material type, data of the SOC greater than 20% are taken to utilize capacity increment point ferric phosphate lithium cell
Analysis method by extracting characteristic value to incremental capacity plot, and records characteristic value position corresponding capacity Q or charging capacity Q ';It is right
In ternary material battery, data of the SOC greater than 20% is taken to seek d2Q/dV2Curve extracts characteristic value, and records characteristic value position
Corresponding capacity Q or charging capacity Q '.Q=SOC*Capinitial。
The datum curve data of ferric phosphate lithium cell are handled, the capacity increment using SOC as abscissa is sought out
Curve, calculation formula areCharacteristic value is extracted, characteristic value position is the pole in incremental capacity plot
Big value position;The datum curve data of ternary battery are handled, the d using voltage as abscissa is sought2Q/dV2Curve, meter
Calculating formula isCharacteristic value is extracted, characteristic value position is Position.For LiFePO4 or
Person's ternary battery is greater than 20% in SOC, and the section less than 100%, the characteristic value of curve is all 2.Record the 1st feature
It is worth the capacity Q of position1Or charging capacity Q '1, the capacity Q of the 2nd characteristic value position2Or charging capacity Q '2.Wherein: Q is battery
Capacity, Q ' be charging capacity, dQ be capacity differential, d2Q is the second-order differential of capacity, Δ QkThe capacity between neighbouring sample point
Difference, V be battery voltage, dV be voltage differential, dV2For the second-order differential of voltage, Δ VkThe voltage between neighbouring sample point
Difference, for each sampled point k, Δ Qk=Qk-Qk-1, Δ Vk=Vk-Vk-1, Δ Qk-1=Qk-1-Qk-2, Δ Vk-1=Vk-1-
Vk-2;
Calculate the coefficient of relationship of capacity and total capacity between battery characteristics value, g=Δ Q/Capinitial;Wherein Δ Q is the 1st
A characteristic value is to the capacity difference between the 2nd characteristic value, Δ Q=| Q2-Q1|;If using history charging curve data, Δ Q
=| Q '2-Q‘1|。
Extract the data for meeting Δ V >=X condition in battery pack charging process.Data include the voltage of all single batteries
Value, the charging time of battery pack and the charged capacity value of battery pack etc..The value range of X is 1mV≤X≤10mV.According to battery
Material type takes data of the SOC greater than 20% using capacity increment analytic approach ferric phosphate lithium cell, by capacity increment
Curve extracts characteristic value, and records the corresponding charging capacity in characteristic value position;For ternary material battery, SOC is taken to be greater than 20%
Data seek d2Q/dV2Curve extracts characteristic value, and records the corresponding charging capacity in characteristic value position.Record all batteries
The charging capacity Q ' of 1st characteristic value position of curve1j, record the charging of the 2nd characteristic value position of all battery curves
Capacity Q '2j, in which: Q1jFor the 1st characteristic value position jth cell capacity, Q2jFor the 2nd characteristic value position jth list
Body battery capacity, Q '1jFor the 1st characteristic value position jth monomer charging capacity, Q '2jFor the 2nd characteristic value position jth list
Bulk charging capacity.
The 1st characteristic value of jth single battery is calculated to the capacity difference between the 2nd characteristic value, Δ Qj=| Q2j-Q1j|;
ΔQjThe 1st characteristic value is also equal to the 2nd characteristic value charging capacity.If certain primary charging curve contains the 1st simultaneously
Characteristic value and the 2nd characteristic value, then Δ Qj=| Q '2j-Q‘1j|.If single charge curve does not include two characteristic value positions, can
It is combined with extracting charging curve from nearest 10 charging curves with that time recently, becomes one comprising 2
Characteristic value point curve.Two curves being combined should meet, and charging current is the same, and charging ambient temperature is differed less than 5 degree
Requirement, otherwise without combination.
Calculate the corresponding total capacity of jth single battery, calculation formula are as follows:
Capj=Δ Qj/g;
Thus the capacity of each single battery in battery pack is calculated one by one.
Example 1:
Target lithium ion battery is the ferric phosphate lithium cell of CATL, nominal capacity Capinitial=180Ah, battery benchmark
Curve data uses SOC-OCV relation curve data, and data break is Δ SOC=2%.Each corresponding appearance of OCV point
Magnitude passes through formula: Q=SOC*180 is calculated;
Take 5 points three times smothing filtering method (to carry out each 2 data before and after smothing filtering position by choosing, in total 5
A data, are fitted using 3 rank multinomials, the numerical value after seeking smothing filtering) seek SOC-OCV relation data capacity increase
Curve is measured, obtained incremental capacity plot is as shown in dQ/dV dotted line in Fig. 2, and Fig. 2 is using SOC as abscissa, and left ordinate is
Battery open circuit voltage OCV, right ordinate are dQ/dV, and the 1st characteristic value position is the A point position in Fig. 2, and corresponding SOC is
55.6%, capacity Q1=100.08Ah, the 2nd characteristic value position is the B point position in Fig. 2 in figure, and corresponding SOC is 84.9%,
Capacity Q2=152.82Ah;
Calculate the coefficient of relationship of capacity and total capacity between battery characteristics value:
G=Δ Q/Capinitial=| Q2-Q1|/Capinitial=(152.82-100.08)/180=0.293
The battery pack of target lithium ion battery assembling, the commercialization energy-storage system being composed in series by 240 single batteries, figure
3 charge data condition is battery pack with 0.35C, and 63A electric current carries out constant-current charge to any one monomer battery voltage and reaches
3.6V, that is, stop charging, and the battery pack SOC for early period of charging is 20%.
The extraction conditions of charge data are Δ V >=2mV, extract the voltage value and charge capacity value of each single battery,.
Voltage value and charge capacity value to battery using 5 points three times smooth filtering method seek incremental capacity plot, the in battery pack
For the incremental capacity plot of No. 1 single battery as shown in figure 3, Fig. 3 is using SOC as abscissa, left ordinate is cell voltage, and the right side is vertical
Coordinate is dQ/dV, and the 1st characteristic value position is the C point position in Fig. 3 in figure, and corresponding charged capacity value is Q '11=67.66Ah,
The 2nd characteristic value position is the D point position in Fig. 3 in figure, and corresponding charged capacity value is Q '21=119.2Ah;
Calculate the total capacity of No. 1 single battery of battery pack:
Cap1=Δ Q1/ g=| Q '21-Q‘11|/g=(119.2-67.06)/0.293=177.9Ah
Same calculation method calculates the total capacity of No. 2 single batteries in battery pack to total appearance of No. 240 single batteries
Amount, and all single battery total capacities and it is depicted as histogram, Fig. 6 is that all cell capacities are straight in battery pack
Square distribution map.
And the capability value calculated using the inventive method is compared with battery true capacity value, volume error rate
Such as Fig. 7, show that maximum capacity error is 3.86% in figure.
Example 2:
Target lithium ion battery is 21700 ternary batteries of power mind, nominal capacity Capinitial=4.5Ah, battery base
Directrix curve data use SOC-OCV relation curve data, and data break is Δ SOC=3%.Each OCV point is corresponding
Capability value passes through formula: Q=SOC*4.5 is calculated;
Take 5 points three times smooth filtering method seek the incremental capacity plot and d of SOC-OCV relation data2Q/dV2It is bent
Line, obtained incremental capacity plot is as shown in the dQ/dV dotted line in Fig. 4, d2Q/dV2As indicated with broken lines in fig. 4, Fig. 4 is curve
Using open-circuit voltage OCV as abscissa, left ordinate is battery SOC, and right 1st ordinate is dQ/dV, and right 2nd ordinate is
d2Q/dV2, the 1st characteristic value position is the E point position in Fig. 4 in figure, and corresponding SOC value is 39.8%, capacity Q1=1.791Ah,
The 2nd characteristic value position is the F point position in Fig. 4 in figure, and corresponding SOC value is 76.4%, capacity Q2=3.438Ah;
Calculate the coefficient of relationship of capacity and total capacity between battery characteristics value:
G=Δ Q/Capinitial=| Q2-Q1|/Capinitial=(3.438-1.791)/4.5=0.336
With 1C, 4.5A electric current carries out constant-current charge to voltage and reaches 4.2V target lithium ion battery, that is, stops charging, charging
The battery pack SOC of early period is 15%.
The extraction conditions of charge data are Δ V >=3mV, extract the voltage value and charge capacity value of battery,.To voltage value and
Charge capacity value using 5 points three times smooth filtering method seek incremental capacity plot and d2Q/dV2Curve, the charging number of the battery
According to incremental capacity plot and d2Q/dV2For curve as shown in figure 5, Fig. 5 is using voltage as abscissa, left ordinate is battery SOC,
Right 1st ordinate is dQ/dV, and right 2nd ordinate is d2Q/dV2, the 1st characteristic value position is the H point in Fig. 5 in figure
It sets, corresponding charged capacity value is Q '11=0.945Ah, the 2nd characteristic value position is the M point position in Fig. 5 in figure, corresponding to fill
Capacitance value is Q '21=2.601Ah;
Calculate the total capacity of No. 1 single battery of battery pack:
Cap1=Δ Q1/ g=| Q '21-Q‘11|/g=(2.443-0.945)/0.336=4.46Ah
Same calculation method calculates the total capacity of other single batteries in battery pack.
Example 3:
Target lithium ion battery is the Soft Roll ternary battery of state's pavilion high-tech, and 15Ah two is simultaneously rear to connect, and nominal capacity is
Capinitial=30Ah, battery datum curve data use history curve data, charging current 6A, and data are extracted
Condition is Δ V >=5mV;
Take 5 points three times smooth filtering method acquire the incremental capacity plot and d of history charging curve data2Q/dV2It is bent
Line, and using charging voltage as abscissa, left ordinate is battery SOC, and right 1st ordinate is dQ/dV, right 2nd ordinate
It is d2Q/dV2, acquire the corresponding charging capacity Q ' in the 1st characteristic value position1=16.53Ah, the 2nd characteristic value position are corresponding
Charging capacity Q '2=26.13Ah;
Calculate the coefficient of relationship of capacity and total capacity between battery characteristics value:
G=Δ Q/Capinitial=| Q '2-Q‘1|/Capinitial=(26.13-16.53)/30=0.32
The battery pack of target lithium ion battery assembling forms battery pack, target lithium ion by 2 and 16 strings by 32 battery cores
With 0.5C, 15A electric current carries out constant-current charge to any monomer battery voltage and reaches 4.2V battery, that is, stops charging, charges early period
Battery pack SOC be 13%.
The extraction conditions of charge data are Δ V >=1mV, extract the voltage value and charge capacity value of each single battery,
After charging complete, after weeding out data of the SOC less than 15%,.Voltage value and charge capacity value to each single battery are adopted
With 5 points three times smooth filtering method seek incremental capacity plot and d2Q/dV2Curve records the 1st spy of each single battery
The charged capacity value of the charged capacity value of value indicative position and the 2nd characteristic value position, wherein the 1st feature of No. 1 single battery
It is Q ' that value position, which corresponds to charged capacity value,11=0.945Ah, the 2nd corresponding charged capacity value in characteristic value position are Q '21=
2.601Ah;
Calculate the total capacity of No. 1 single battery in battery pack:
Cap1=Δ Q1/ g=| Q '21-Q‘11|/g=(21.56-11.576)/0.32=31.2Ah
Same calculation method calculates in battery pack No. 2 single batteries to the total capacity of No. 16 single batteries.
Embodiment 4:
Case study on implementation 4 is similar with case study on implementation 3, unlike, when handling charging curve, single charge curve does not include
Two characteristic value positions, therefore extract charging curve from nearest 10 charging curves and be combined with that time recently, make it
It include 2 characteristic value point curves as one.It is chosen in nearest 10 charging curves first with worst curve charging current one
Sample and environment temperature differ the curve less than 5 degree, and packet can be combined into that nearest secondary curve by choosing again in residual curve
Curve containing two characteristic value position curves, if there is multiple curves are met the requirements, then preferential selection environment temperature is immediate;
If environment temperature is the same, newest charging curve (i.e. the time upper curve closest to that last secondary charging) is preferentially chosen.
Data extraction conditions can be the data of Δ V >=10mV condition.Aforementioned formula is recycled to calculate monomer electricity in battery pack one by one
Tankage.
The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention.The technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the above embodiments and description only describe this
The principle of invention, without departing from the spirit and scope of the present invention, various changes and improvements may be made to the invention, these changes
Change and improvement all fall within the protetion scope of the claimed invention.The claimed scope of the invention by appended claims and its
Equivalent thereof.