WO2011155017A1 - 充電状態推定装置 - Google Patents

充電状態推定装置 Download PDF

Info

Publication number
WO2011155017A1
WO2011155017A1 PCT/JP2010/059626 JP2010059626W WO2011155017A1 WO 2011155017 A1 WO2011155017 A1 WO 2011155017A1 JP 2010059626 W JP2010059626 W JP 2010059626W WO 2011155017 A1 WO2011155017 A1 WO 2011155017A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
charge
battery
state
calculation unit
Prior art date
Application number
PCT/JP2010/059626
Other languages
English (en)
French (fr)
Inventor
吉岡　省二
朗子 田渕
啓太 畠中
英俊 北中
Original Assignee
三菱電機株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to KR1020127027258A priority Critical patent/KR101338639B1/ko
Priority to US13/640,056 priority patent/US8975897B2/en
Priority to EP10852857.1A priority patent/EP2579059B1/en
Priority to BR112012027318A priority patent/BR112012027318A2/pt
Priority to JP2012519152A priority patent/JP5511951B2/ja
Priority to CN201080067268.7A priority patent/CN102918411B/zh
Priority to RU2012155314/28A priority patent/RU2565339C2/ru
Priority to MX2012012939A priority patent/MX2012012939A/es
Priority to PCT/JP2010/059626 priority patent/WO2011155017A1/ja
Priority to AU2010354957A priority patent/AU2010354957B2/en
Publication of WO2011155017A1 publication Critical patent/WO2011155017A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/389—Measuring internal impedance, internal conductance or related variables
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/396—Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery

Definitions

  • the present invention relates to a state-of-charge estimating apparatus that repeatedly performs charging and discharging, for example, estimating a state of charge (SOC) of a battery such as a lithium ion secondary battery.
  • SOC state of charge
  • Patent Document 1 uses a simplified equivalent circuit model of a battery when estimating an open circuit voltage. In order to obtain an open circuit voltage by calculation, the voltage, current, temperature, resistance, capacitance, etc. Although the control information and the parameter information representing the battery characteristics are used, the same battery parameter is used for the constant voltage control and the constant current control.
  • the charge state value can be estimated from the open circuit voltage using an equivalent circuit model from the measured current and voltage, but the internal resistance that greatly affects the accuracy of the estimated charge state value is the current, Since it is composed of a diffusion resistance or an electrode reaction resistance that exhibits nonlinear characteristics with respect to voltage, it is difficult to estimate a state of charge using a conventional equivalent circuit model described by a linear resistance.
  • the present invention has been made in view of the above, and an object of the present invention is to obtain a state-of-charge estimation device capable of accurately estimating a state of charge and a deterioration state of a battery.
  • the present invention is a charging state estimation device that is connected to a power storage device in which a plurality of batteries are connected and estimates a charging state value indicating a remaining capacity of the power storage device. And based on a battery capacity, a previous value of the estimated state of charge, and a current value flowing in and out between the power storage device and a current control device that controls a charge / discharge amount of the power storage device.
  • a first calculation unit that calculates the current value of the calculated charge state estimated value as a first charge state estimated value, and at the time of constant current control that charges the power storage device with a constant current, an equivalent circuit model of the battery, While calculating the charging state estimated value calculated based on the voltage value of the battery as the current value of the second charging state estimated value, at the time of constant voltage control for charging the power storage device with a constant voltage, the battery Equivalent circuit model
  • a second calculation unit that calculates a current value of the estimated state of charge calculated in consideration of a resistance change of the battery based on the voltage value of the battery and the second state of charge as the second state of charge estimated value; And a correction calculation unit that periodically corrects the state of charge estimated value with the second state of charge estimated value.
  • the open circuit voltage and the charging state of the battery are estimated from the equivalent circuit model including the non-linear resistance at the time of constant voltage control, the charging state of the battery In addition, there is an effect that the deterioration state can be accurately estimated.
  • FIG. 1 is a configuration diagram of a power storage device to which a charge state estimation device according to an embodiment of the present invention is applied.
  • FIG. 2 is a configuration diagram of an electric power storage system to which the charge state estimation device according to the embodiment of the present invention is applied.
  • FIG. 3 is a flowchart for explaining the SOC estimation process performed by the state-of-charge estimation apparatus according to the embodiment of the present invention.
  • FIG. 4 is a configuration diagram of the first arithmetic unit according to the present embodiment.
  • FIG. 5 is a diagram showing a distributed constant system equivalent circuit model of the electricity storage device applicable to the second calculation according to the present embodiment.
  • FIG. 6 is a diagram showing a lumped-constant equivalent circuit model of the electricity storage device applicable to the second calculation according to the present embodiment.
  • FIG. 1 is a configuration diagram of a power storage device to which a charge state estimation device according to an embodiment of the present invention is applied.
  • FIG. 2 is a configuration diagram of an electric power storage system to which the
  • FIG. 7A is a diagram illustrating a relationship between the second calculation unit 31 and the data bank 202.
  • FIG. 7B is a diagram for explaining the operation of the second estimation process by the second calculation unit 31.
  • FIG. 8A shows the configuration of the deterioration amount calculation unit 35.
  • FIG. 8-2 is a diagram for explaining operations of the voltage simulation unit 231 and the resistance / capacitance calculation unit 232.
  • FIG. 9 is a diagram for explaining the relationship between the product of the resistance and the capacitance of the capacitor and the battery capacity.
  • FIG. 10 is a diagram for explaining the operation of the resistance value calculation unit according to the present embodiment.
  • FIG. 11 is a diagram for explaining a change in resistance during constant voltage control.
  • FIG. 12 is a diagram showing the SOC calculation accuracy in a deteriorated battery.
  • FIG. 1 is a configuration diagram of a power storage device 1 to which a charge state estimation device according to an embodiment of the present invention is applied.
  • the power storage device 1 includes battery modules 11-1 to 1n-m in which k unit cells are connected in series.
  • the battery module 11-1 is the first battery module in a row, and is obtained by connecting the cells 1-1 to 1-k in series.
  • the battery module 12-1 is the second battery module in the first row
  • the battery module 1n-1 is the nth battery module in the first row.
  • the battery module 11-2 is the first battery module in the second row
  • the battery module 12-2 is the second battery module in the second row
  • the battery module 1n-2 is the second battery module in the second row.
  • the battery module 11-m is the first battery module in the m row
  • the battery module 12-m is the second battery module in the m row
  • the battery module 1n-m is the nth battery module in the m row. Therefore, the total number of cells is n ⁇ m ⁇ k.
  • Each of the cells 1-1 to 1-k is a power storage device that can be repeatedly charged and discharged, such as a lithium ion secondary battery, and the charge state value can be observed from the open circuit voltage value.
  • Nickel metal hydride batteries, lead storage batteries, electric double layer capacitors, lithium ion capacitors, and the like can also be used as power storage devices constituting the power direct storage device.
  • the power storage device 1 may be provided with a breaker, a battery monitoring device, and the like, which are omitted here.
  • the terminal voltage of the entire power storage device is the total voltage Vall
  • the total current of the charge / discharge current is Iall
  • the charging direction is positive.
  • the voltage caused by the resistance component of the conductor or cable used for connection between the terminals is added to Vall.
  • the voltage at zero current is the open circuit voltage.
  • the open circuit voltage and the SOC have a one-to-one relationship under a constant temperature environment, and the SOC generally indicates a monotonically increasing function with respect to the open circuit voltage.
  • FIG. 2 is a configuration diagram of a power storage system to which the state-of-charge estimation device according to the embodiment of the present invention is applied.
  • the power storage device 1 is controlled by the current control device 3 based on a command from the controller 2.
  • the controller 2 includes an A / D converter 203, a data bank 202, and a parameter calculation unit 201.
  • the controller 2 is a charge state estimation device according to the present embodiment.
  • the A / D converter 203 includes analog signals sent from the total voltage sensor 4a, the total current detection sensor 4b, and the temperature sensor 4c in the power storage system, and each cell voltage sensor (not shown) in the power storage device 1.
  • the analog signal 4d to be sent is converted into a digital signal.
  • the data bank 202 stores battery data and stores an expression describing the relationship between the open circuit voltage and the SOC, or a data table indicating this relationship. Since the relationship between the open circuit voltage and the SOC differs depending on the type of the electricity storage device to be used and the electrode material types of the positive electrode and the negative electrode, it is obtained by input or measurement in advance. As a method for obtaining it, open circuit voltage data is obtained when a known constant current is applied for a certain period of time and the SOC is changed at a constant interval in an environment controlled at a constant temperature. It is desirable to measure the open circuit voltage after maintaining a zero current state for several hours. The open circuit voltage can also be obtained from the state in which the power storage device 1 is incorporated in the system, and is executed at the beginning of the control program, or the program is executed using a long unused period or maintenance period. It can also be done manually.
  • the method of obtaining from the fully discharged state is a method of obtaining an open circuit voltage with respect to the SOC by charging at a constant SOC interval with the state of charge at the time of complete discharge discharged to the lower limit voltage of the battery being zero.
  • the method from the fully charged state is a method in which the SOC is changed at a constant SOC interval by setting the SOC in a state of being charged by constant current and constant voltage charging up to the upper limit voltage of the battery as 100% and changing the SOC. Comparing the two results in a more stable state when measured from the state of charge, and the accuracy of the data is higher.
  • the current control device 3 controls the charge / discharge amount of the power storage device 1 according to a command from the parameter calculation unit 201 in the controller 2.
  • the current control device 3 outputs the power supplied from the AC grid 7 via the power conversion device 6 and the power stored in the power storage device 1 to the motor 5.
  • the current control device 3 stores the power generated by the motor / generator 5 in the power storage device 1, and for the power exceeding the capacity of the power storage device 1, the AC grid 7 is connected via the power conversion device 6.
  • the AC grid 7 is, for example, an AC power supply network that supplies power to the railway vehicle.
  • the parameter calculation unit 201 includes a first calculation unit 30, a second calculation unit 31, and a correction calculation unit 32 as illustrated in FIG. 2.
  • the first computing unit 30 is configured to calculate a first charge state estimated value (hereinafter simply referred to as “(first SOC)”) based on the battery capacity (Ah) and the integrated value of the energized electricity during battery charging / discharging. Is calculated.
  • the second calculation unit 31 uses the open circuit voltage calculated from the equivalent circuit model having the resistance component and the capacitor component, and the relationship between the open circuit voltage and the charge state estimated value, to calculate a second charge state estimated value (hereinafter referred to as the second charged state estimated value Simply referred to as “(second SOC)”.
  • the second calculation unit 31 calculates the charge state estimated value from the open circuit voltage during the second estimation process for estimating the charge state estimated value in a current zero state and the constant current control in which the voltage changes.
  • a third estimation process for estimation and a fourth estimation process for estimating the charge state estimation value from the open circuit voltage during constant voltage control controlled at a constant voltage are executed.
  • the correction calculation unit 32 periodically corrects the SOC estimated by the first calculation unit 30 with the SOC estimated by the second calculation unit 31.
  • FIG. 3 is a flowchart for explaining the SOC estimation process performed by the state-of-charge estimation apparatus according to the embodiment of the present invention.
  • the flowchart shown in FIG. 3 is repeated at the same time interval as the interval for acquiring data.
  • This time interval means an acquisition interval of each data of current, voltage, and temperature, and an interval of milliseconds to several minutes is appropriate.
  • the SOC (N-1), the electron / ion resistance R0 (N-1), the electrode reaction resistance R (N-1), and the capacitance component obtained in the previous calculation flow C (N-1) and the like are recorded (step S10).
  • the second calculation unit 31 first switches the second to fourth estimation processing methods depending on whether the constant current control or constant voltage control is being performed. For example, in the open circuit state where the current is zero (Step S12, No to S13, Yes), the second calculation unit 31 executes the second estimation process (Step S14). Further, the second calculation unit 31 executes the third estimation process (step S15) during the constant current control in which the voltage changes (step S12, No to S13, No). Furthermore, the 2nd calculating part 31 performs a 4th estimation process at the time of the constant voltage control controlled by a fixed voltage (step S12, Yes) (step S16). The correction calculation unit 32 corrects the first SOC with the second SOC.
  • FIG. 4 is a configuration diagram of the first arithmetic unit 30 according to the present embodiment.
  • the first calculation unit 30 includes, as main components, an average current calculation unit 211, an energized electricity amount calculation unit 212, an SOC change amount calculation unit 213, and an SOC calculation unit 214.
  • the total current Iall flowing through the power storage device 1 is input to the average current calculation unit 211.
  • the average current calculation unit 211 multiplies the total current Iall by a predetermined gain to obtain an average current value. That is, the average current value is obtained by dividing the total current Iall by the parallel number m of the battery modules 11-1 to 1nm.
  • the energization electricity calculation unit 212 integrates the energization current during charging and discharging with the calculation cycle as ⁇ t, and calculates the electricity after energization for a predetermined time.
  • the SOC change amount calculation unit 213 divides the amount of electricity (coulomb) obtained by the energized electricity amount calculation unit 212 by the battery capacity (Ah) and 3600 (s), and multiplies by 100 to obtain the change amount ⁇ SOC (%). calculate.
  • the SOC calculation unit 214 adds the amount of change ⁇ SOC to the SOC (N ⁇ 1) obtained in the previous calculation flow at the time of charging, and subtracts the amount of change ⁇ SOC from the SOC (N ⁇ 1) at the time of discharging.
  • the current value SOC (N) of the value is obtained.
  • the SOC (N ⁇ 1) is stored in the data bank 202, and the SOC calculation unit 214 estimates the SOC (N) using the SOC (N ⁇ 1) from the data bank 202. This is the simplest method for determining the SOC, but often includes errors in current measurements.
  • the current value information of the battery capacity (Ah) is changed and input to the data bank 202.
  • the calculation period ⁇ t may be generated inside the parameter calculation unit 201 or may be generated outside the parameter calculation unit 201.
  • FIG. 5 is a diagram illustrating a distributed constant system equivalent circuit model of the electricity storage device applicable to the second arithmetic unit 31 according to the present embodiment.
  • the second computing unit 31 estimates the SOC, resistance, and capacitance of the capacitor based on the equivalent circuit model. Strictly speaking, the second calculation unit 31 obtains an open circuit voltage by fitting a measured value of current and a measured value of voltage to a numerical model discretized based on the distributed constant system equivalent circuit shown in FIG. .
  • the calculation formula includes the negative terminal resistance 8a, the positive terminal resistance 8b, the negative electrode layer electronic resistance 9a, the positive electrode layer electronic resistance 9b, and the negative electrode interface capacitor of the distributed constant equivalent circuit model shown in FIG.
  • Capacitance 10a Capacitance capacitance 10b of positive electrode interface, Resistance component 11a of negative electrode interface, Resistance component 11b of positive electrode interface, Potential difference 12a generated at negative electrode interface, Potential difference 12b generated at positive electrode interface, Resistance of electrolyte in negative electrode 13a, the resistance 13b of the electrolyte in the positive electrode, and the resistance 14 of the electrolyte in the separator.
  • a plurality of resistance values, capacitance values, and open circuit voltages are obtained by fitting the actually measured current value and voltage temporal change measurement data and calculated values.
  • the capacitance 10a of the capacitor at the negative electrode interface and the capacitance 10b of the capacitor at the positive electrode interface are due to the electric double layer formed at the electrode active material / electrolyte interface, and are expressed in units of farads.
  • the capacitance 10a of the capacitor at the negative electrode interface and the capacitance 10b of the capacitor at the positive electrode interface are proportional to the surface area of the electrode active material, and also vary depending on the properties of the electrolyte and the electrode potential.
  • the resistance component 11a at the negative electrode interface and the resistance component 11b at the positive electrode interface are resistances when charge carriers change from ions to electrons and from electrons to ions, and are resistances generated at the electrode / electrolyte interface, Includes diffused resistance.
  • the distributed constant system equivalent circuit model shown in FIG. 5 represents one single cell (eg, 1-1) shown in FIG. 1, and the current I flowing through this distributed constant system equivalent circuit model This is a value obtained by dividing the current Iall by the parallel number of the battery modules 11-1 to 1nm.
  • FIG. 6 shows a lumped parameter equivalent circuit model that can save calculation resources by shortening the calculation cycle.
  • FIG. 6 is a diagram illustrating a lumped-constant equivalent circuit model of an electricity storage device applicable to the second arithmetic unit 31 according to the present embodiment.
  • the lumped-constant equivalent circuit model is a simplification of the distributed-constant equivalent circuit model of FIG. 5.
  • the current I flowing through the lumped-constant equivalent circuit model is the total current Iall as shown in FIG. The value is divided by the parallel number of 1 to 1 nm.
  • Each element in this model is configured as follows. That is, the lumped-constant equivalent circuit model shown in FIG.
  • resistor 15 related to electrode reaction
  • capacitance 18 of a capacitor generated at the electrode interface
  • resistor 16 A resistance 16 relating to electrons / ions
  • electromotive force portion 17 A resistance 16 relating to an open circuit voltage
  • the value of the resistor 15 is R
  • the capacitance value of the capacitance 18 is C
  • the value of the resistor 16 is R0 (electron / ion resistance).
  • the current flowing through the resistor 15 is I1
  • the current flowing through the capacitance 18 is I2
  • the sum of these is the current I. Since the voltage applied to the capacitor (capacitance 18) is equal to the voltage across the resistor 15 through which the current I1 flows, and the time change of the charge Q stored in the capacitor (capacitance 18) corresponds to the current I2.
  • a differential equation relating to the charge Q in equation (1) is obtained.
  • Equation (2) is obtained.
  • the charge Q (N) at the time N can be expressed as the following equation (2) by using the charge Q (N-1) and the current I at the previous time (N-1).
  • the current I2 flowing through the capacitor can be expressed by the change with time of the charge Q shown in the equation (3).
  • the open circuit voltage Voc is calculated from the equation (4).
  • FIG. 7 is a diagram for explaining the configuration and operation of the second arithmetic unit 31 according to the present embodiment.
  • FIG. 7A is a diagram illustrating the relationship between the second calculation unit 31 and the data bank 202
  • FIG. 7B is a diagram for explaining the operation of the second estimation process performed by the second calculation unit 31.
  • FIG. 7A is a diagram illustrating the relationship between the second calculation unit 31 and the data bank 202
  • FIG. 7B is a diagram for explaining the operation of the second estimation process performed by the second calculation unit 31.
  • the resistance value calculation unit 36 includes a deterioration amount calculation unit 35, a diffusion species concentration calculation unit 225, a current value calculation unit 226, and a diffusion resistance calculation unit 227. It is comprised.
  • the voltage change determination part 221 determines the presence or absence of the change of the battery voltage V in the calculation period ⁇ t. That is, the voltage change determination unit 221 performs the determination in step S12 in FIG.
  • the energization determining unit 222 determines whether or not the current I is zero when the battery voltage V is changing (No in step S12). That is, the energization determination unit 222 performs the determination in step S13 in FIG.
  • the open circuit voltage calculation unit 224 calculates the open circuit voltage Voc from the battery voltage V and the value R0 (N ⁇ 1) of the resistor 16 stored in the data bank 202 using the equation (4). It is assumed that the value R 0 (N ⁇ 1) of the resistor 16 is stored in the data bank 202.
  • the SOC calculation unit 214 takes in the open circuit voltage Voc from the open circuit voltage calculation unit 224 and estimates the SOC (N) using the relationship between the open circuit voltage Voc stored in the data bank 202 and the SOC.
  • the open circuit voltage calculation unit 224 calculates the battery voltage V, the resistor 15, and the resistor 16 from the equation (4). , And the current I is zero, and the open circuit voltage Voc is obtained.
  • the SOC calculation unit 214 takes in the open circuit voltage Voc from the open circuit voltage calculation unit 224, and estimates the second SOC using the relationship between the open circuit voltage Voc and the SOC stored in the data bank 202.
  • the open circuit voltage calculation unit 224 obtains the open circuit voltage Voc based on the battery voltage V, the resistor 15, the resistor 16, and the current I from the equation (4).
  • the SOC calculation unit 214 estimates the second SOC using the relationship between the open circuit voltage Voc and the SOC. The accuracy of the second estimation process for obtaining the SOC without passing current is higher.
  • FIG. 8 is a configuration diagram of the deterioration amount calculation unit 35 according to the present embodiment.
  • FIG. 8A shows the configuration of the deterioration amount calculation unit 35
  • FIG. 8B is a diagram for explaining the operation of the voltage simulation unit 231 and the resistance / capacitance calculation unit 232.
  • FIG. 9 is a diagram for explaining the relationship between the product of the resistance and the capacitance of the capacitor and the battery capacity.
  • the deterioration amount calculation unit 35 shown in FIG. 8A includes a voltage simulation unit 231 and a resistance / capacitance calculation unit 232 as main components.
  • the deterioration amount calculation unit 35 estimates the deterioration state of the battery from the current and voltage data for a certain period of time based on changes in the resistance 15, the resistance 16, and the capacitance 18.
  • the deterioration amount calculation unit 35 estimates the progress of the battery deterioration state by sequentially calculating the resistance component and the capacitance component of the capacitor in the assumed control mode, for example, a plurality of voltage data By fitting the calculated values, the values of the resistance 15, the resistance 16, and the capacitance 18 are set as the electrode reaction resistance R (N), the electron / ion resistance R0 (N), and the capacitance component C (N), respectively. Identify.
  • the voltage data to be fitted is p data from time t1 to tp, and the individual data time interval is ⁇ t.
  • the voltage simulation unit 231 uses the charge accumulated in the capacitor 18 immediately before time t1 as Q (k ⁇ 1), the electrode reaction resistance R (N), the electron / ion resistance R0 (N), and the capacitance Substituting the previous value for each of the components C (N), substituting the values used when calculating the SOC in FIG. 7-1 for the current I and the open circuit voltage Voc, and using the equation (5), the voltage for the time t1 Vk is calculated.
  • the charge Q (N) in the equation (5) is calculated using the equation (2).
  • the voltage simulation unit 231 calculates from t2 to tp.
  • the resistance / capacitance calculation unit 232 compares the calculated value (voltage Vk) calculated by the voltage simulation unit 231 with the actual measurement data (Vp), and the sum ⁇ V of the difference between the actual measurement value and the calculation value is determined as the determination value. If it is less than the product of ⁇ and the number of data p, it is regarded as a match. On the other hand, when the sum ⁇ V of the difference between the actual measurement value and the calculated value is equal to or greater than the product of the determination value ⁇ and the number of data p, the resistance / capacitance calculation unit 232 determines the electrode reaction resistance R (N), the electron / ion resistance R0 ( N), changing the value of the capacitance component C (N).
  • the resistance / capacitance calculation unit 232 increases the electrode reaction resistance R (N), the electron / ion resistance R0 (N), and the capacitance component C (N) if the discharge voltage is lower than the calculated value, and vice versa. Make it smaller.
  • the resistance / capacity calculation unit 232 calculates the battery capacity using the relationship between the product and the capacity of the deteriorated battery (see FIG. 9). However, since the battery deterioration is not a reaction that proceeds rapidly, it is not necessary to calculate the deterioration in the data acquisition cycle, and it is sufficient to have a single frequency immediately after starting or immediately before stopping.
  • the resistor 15 at the time of constant voltage control rises with time and does not take a constant value.
  • the reason why the value of the resistor 15 is increased is that the resistance value is proportional to the reciprocal of the lithium ion concentration C Li (N) in the electrode active material, as shown in the Nernstein-Einstein equation (6). In constant voltage control, this is due to a decrease in the number of movable lithium ions.
  • D Li is a diffusion constant in the lithium ion active material
  • T is a module temperature
  • A is a proportionality constant.
  • FIG. 10 is a diagram for explaining the operation of the resistance value calculation unit 36 according to the present embodiment
  • FIG. 11 is a diagram for explaining a change in resistance during constant voltage control.
  • the diffusion species concentration calculation unit 225 calculates the equation (7) based on the boundary condition from C Li (N ⁇ 1), which is the concentration of mobile lithium ions in the active material immediately before entering constant voltage control, and the diffusion constant D Li. It solves and calculates lithium ion density
  • the current value calculation unit 226 calculates the current value from the change in the battery voltage V. Based on the current value calculated by the current value calculation unit 226 and the voltage at the time of constant current control, the diffusion resistance calculation unit 227 calculates the resistance 15 of the resistor 15 when the voltage is not constant (CC charging: constant ⁇ current charging). Estimate the value. Further, the diffusion resistance calculation unit 227 estimates the value of the resistance 15 when the voltage is constant (CV charging: constant voltage charging) using the equation (6).
  • the value of the resistance 15 relating to the electrode reaction does not change during CC charging, but increases as shown in FIG. 11 during CV charging. This increase in resistance is temporary, occurring between a few seconds and a few hours under constant voltage control conditions.
  • the second arithmetic unit 31 also executes a process of excluding the increase in resistance at a constant voltage as shown in FIG. 11 from the calculation for estimating capacity deterioration. To do.
  • the short-term resistance increase rate calculated here that is, the value obtained by dividing the resistance value increase amount in FIG. 11 by the constant voltage control time indicates the deterioration of the battery related to mass transfer.
  • the 2nd calculating part 31 concerning this Embodiment estimates the battery life by long-term deterioration using the fact that this value becomes large with the progress of deterioration.
  • FIG. 12 is a diagram showing the SOC calculation accuracy in a deteriorated battery.
  • Comparative Example 1 is a transition of the first SOC calculated by the first calculation unit 30 based on integration of energized electricity. Since the current value used for the calculation generally includes an error, the SOC deviates from the true value.
  • Example 1 shows the transition of the SOC estimated by the state-of-charge estimation device according to the present embodiment.
  • the charging / discharging start point 19 and the charging / discharging end point 20 are correct SOCs obtained after being left in an open circuit for 3 days or more.
  • the current value is detected and integrated smaller during charging than when actually charged, and is detected and integrated greatly during discharging, resulting in a lower value than the true value over the entire charge / discharge.
  • the SOC estimation method according to Comparative Example 2 since the SOC is calculated without considering the capacity reduction due to deterioration, both charging and discharging are calculated to be larger than the true value of the SOC, and therefore, the SOC changes excessively.
  • the deterioration amount calculation unit 35 considers a decrease in capacity due to battery deterioration, and the resistance value calculation unit 36 uses the resistance value calculation unit 36 to estimate the increase in resistance at a constant voltage. Since it is excluded from the calculation, it is possible to improve the SOC calculation accuracy.
  • the state-of-charge estimation device includes the battery capacity, the previous value SOC (N ⁇ 1) of the state of charge estimation value, and the current that controls the charge / discharge amount of the power storage device 1.
  • a first calculation unit 30 that calculates a current value of a charge state estimated value calculated based on a current flowing in and out between the control device 3 and the power storage device 1 as a first charge state estimated value; At the time of current control, the estimated state of charge calculated based on the equivalent circuit model of the battery and the voltage of the battery is calculated as the current value of the second state of charge estimated value, while at the time of constant voltage control, the equivalent circuit model of the battery is calculated.
  • the third SOC at the time of current control and the fourth SOC at the time of constant voltage control can be estimated, and the state of charge and deterioration of the battery can be accurately estimated as compared with the prior art. .
  • the equivalent circuit model of the battery used in the second arithmetic unit 31 is a lumped constant equivalent circuit model having one parallel circuit of a capacitor and a nonlinear resistor, or a plurality of parallel circuits of a capacitor and a nonlinear resistor. Since it is configured with the distributed constant system equivalent circuit model that it has, it can save calculation resources by using the lumped constant system equivalent circuit model, and it can perform strict calculations by using the distributed constant system equivalent circuit model. It is also possible to do this.
  • the second calculation unit 31 includes the charge Q (N ⁇ 1) stored in the capacitor, the capacitance value C (N ⁇ 1) of the capacitor, and the resistance value (R) of the battery. (N-1), R0 (N-1)) and the voltage V of the battery are used to calculate the second charge state estimated value, so that the second charge state estimated value can be obtained with a simple configuration. Can be estimated.
  • the second calculation unit 31 also includes a plurality of continuous voltage measurement values (Vp) measured after the current changes, a charge Q (N ⁇ 1), a capacitance value C (N ⁇ 1), and a resistance.
  • Vp continuous voltage measurement values
  • the sum ⁇ V of the difference between the calculated voltage value (Vk) calculated based on the values R (N ⁇ 1) and R0 (N ⁇ 1) is a predetermined value (product of the judgment value ⁇ and the number of data p) Until the capacitance value C (N) and the resistance values R (N) and R0 (N) are calculated, the sum ⁇ V of the difference between the voltage measurement value (Vp) and the voltage calculation value (Vk) is a predetermined value.
  • the deterioration amount calculation unit 35 that outputs the capacitance C (N) and the resistance values R (N) and R0 (N) when they coincide with each other as an index for estimating the deterioration state of the battery is provided. Compared to technology, it is possible to estimate the state of charge value by the equivalent circuit model with higher accuracy.
  • the first calculation unit 30 includes the deterioration amount calculation unit 35. It is possible to accurately calculate the first charge state estimated value using the battery capacity calculated in (1).
  • the second calculation unit 31 uses the movable lithium ion concentration C Li (N ⁇ 1) and the diffusion constant D Li in the lithium ion active material to determine the resistance value R (N ) Is calculated, the increase in resistance at a constant voltage can be excluded from the calculation for estimating the capacity deterioration, and the SOC calculation accuracy can be improved.
  • the present invention can be applied to a state-of-charge estimation device that estimates SOC in a storage battery such as a secondary battery, and is particularly useful as an invention that can improve the estimation accuracy of SOC.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Secondary Cells (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Abstract

　電池容量と、充電状態推定値の前回値と、電力貯蔵装置１の充放電量を制御する電流制御装置３と電力貯蔵装置１との間で流出入する電流と、に基づいて演算した充電状態推定値の今回値を第１の充電状態推定値として算出する第１の演算部３０と、定電流制御時には、電池の等価回路モデルと電池の電圧とに基づいて演算した充電状態推定値を第２の充電状態推定値の今回値として算出する一方で、定電圧制御時には、電池の等価回路モデルと電池の電圧とに基づいて、電池の抵抗変化を考慮して演算した充電状態推定値の今回値を第２の充電状態推定値として算出する第２の演算部３１と、第１の充電状態推定値を第２の充電状態推定値で定期的に補正する補正演算部３２と、を備える充電状態推定装置。

Description

充電状態推定装置
　本発明は、充電と放電を繰り返し行う、例えばリチウムイオン二次電池等の電池の充電状態値（ＳＯＣ:State Of Charge）の推定を行う充電状態推定装置に関するものである。
　自動車や鉄道車両システムにおいて、電池を電力バッファとして、車両運行時の駆動モータへの電力供給、減速または停止時の電気エネルギー回収を行う蓄電システムでは、車両運行時における電池の充電状態を正確に把握する必要がある。
　従来、充電状態値を推定する方式としては、充放電電流値の積算値から計算する方式と、電池の開回路電圧から推定する方式が知られており、これらを併用する従来技術では、電池が充放電している際には電流積算から充電状態値を推定し、待機時には開回路電圧を推定することで充電状態値を推定するという方法が採用されていた（例えば下記特許文献１）。また、充電状態値は開回路電圧からは推定可能であるが、充放電を行っている閉回路において開回路電圧を予測する技術が必要となる。これについては、閉回路時の通電電流と電池電圧の測定値から、電池の等価回路モデルを用いて内部抵抗を算出する方法がある（例えば下記特許文献１、２、３）。
特開２００８－１９９７２３号公報 特開平８－１４０２７０号公報 特開２００３－７５５１８号公報
　蓄電システムを搭載した鉄道車両などの走行可能距離は、充電状態値から推定されるため、蓄電システムにおける電池充電状態値の正確な把握は、機器の効率的運用を可能にする。直並列した多数の電池で構成される電力貯蔵装置の充電状態値を推定する場合、電力貯蔵装置を構成する電池すべての電圧、電流、充電状態値を測定または推定することが望ましい。そして、電池の充電状態は、一般に開回路電圧と一対一の関係にあるため、開回路電圧から推定することが可能である。
　特許文献１は、開回路電圧の推定に際して電池の単純化した等価回路モデルを用いているが、開回路電圧を演算で求めるために、電圧、電流の他、温度、抵抗、静電容量などの制御情報と電池特性を表すパラメータの情報とを用いているが、定電圧制御と定電流制御とで同じ電池パラメータを利用しているため、すべての制御条件に対応できないという課題があった。
　特許文献２、３についても同様に、測定電流と電圧から等価回路モデルを用いて開回路電圧から充電状態値を推定できるが、推定する充電状態値の精度に大きく影響する内部抵抗は、電流、電圧に対して非線形特性を示す拡散抵抗や電極反応抵抗で構成されるため、線形抵抗で記述された従来の等価回路モデルによる充電状態値の推定は困難である。
　本発明は、上記に鑑みてなされたものであって、電池の充電状態および劣化状態を精度よく推定することが可能な充電状態推定装置を得ることを目的とする。
　上述した課題を解決し、目的を達成するために、本発明は、電池を複数接続した電力貯蔵装置に接続され、前記電力貯蔵装置の残存容量を示す充電状態値を推定する充電状態推定装置であって、電池容量と、前記充電状態推定値の前回値と、前記電力貯蔵装置の充放電量を制御する電流制御装置と前記電力貯蔵装置との間で流出入する電流値と、に基づいて演算した充電状態推定値の今回値を第１の充電状態推定値として算出する第１の演算部と、前記電力貯蔵装置を一定の電流で充電する定電流制御時には、前記電池の等価回路モデルと前記電池の電圧値とに基づいて演算した充電状態推定値を第２の充電状態推定値の今回値として算出する一方で、前記電力貯蔵装置を一定の電圧で充電する定電圧制御時には、前記電池の等価回路モデルと前記電池の電圧値とに基づいて、電池の抵抗変化を考慮して演算した充電状態推定値の今回値を第２の充電状態推定値として算出する第２の演算部と、前記第１の充電状態推定値を第２の充電状態推定値で定期的に補正する補正演算部と、を備えたことを特徴とする。
　この発明によれば、通電時の温度、電流変化、電圧変化から電池の開回路電圧および充電状態を、定電圧制御時に非線形抵抗を含む等価回路モデルから推定するようにしたので、電池の充電状態および劣化状態を精度よく推定することができるという効果を奏する。
図１は、本発明の実施の形態にかかる充電状態推定装置が適用される電力貯蔵装置の構成図である。 図２は、本発明の実施の形態にかかる充電状態推定装置が適用される電力貯蔵システムの構成図である。 図３は、本発明の実施の形態にかかる充電状態推定装置によるＳＯＣ推定処理を説明するためのフローチャートである。 図４は、本実施の形態にかかる第１の演算部の構成図である。 図５は、本実施の形態にかかる第２の演算に適用可能な蓄電デバイスの分布定数系等価回路モデルを示す図である。 図６は、本実施の形態にかかる第２の演算に適用可能な蓄電デバイスの集中定数系等価回路モデルを示す図である。 図７－１は、第２の演算部３１とデータバンク２０２との関係を示す図である。 図７－２は、第２の演算部３１による第２の推定処理の動作を説明するための図である。 図８－１は、劣化量算出部３５の構成を示すものである。 図８－２は、電圧シミュレーション部２３１および抵抗／容量算出部２３２の動作を説明するための図である。 図９は、抵抗とコンデンサの静電容量との積と、電池容量との関係を説明するための図である。 図１０は、本実施の形態にかかる抵抗値算出部の動作を説明するための図である。 図１１は、定電圧制御時における抵抗の変化を説明するための図である。 図１２は、劣化した電池におけるＳＯＣ演算精度を示す図である。
　以下に、本発明にかかる充電状態推定装置の実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
実施の形態．
　図１は、本発明の実施の形態にかかる充電状態推定装置が適用される電力貯蔵装置１の構成図である。
　電力貯蔵装置１は、ｋ個の単電池を直列に接続した電池モジュール１１－１～１ｎ－ｍを有して構成されている。例えば、電池モジュール１１－１は、１列１番目の電池モジュールであり、単電池１－１～１－ｋを直列接続したものである。以下同様に、電池モジュール１２－１は、１列２番目の電池モジュール、電池モジュール１ｎ－１は、１列ｎ番目の電池モジュールである。電池モジュール１１－２は、２列１番目の電池モジュール、電池モジュール１２－２は、２列２番目の電池モジュール、電池モジュール１ｎ－２は、２列ｎ番目の電池モジュールである。電池モジュール１１―ｍは、ｍ列１番目の電池モジュール、電池モジュール１２―ｍは、ｍ列２番目の電池モジュール、電池モジュール１ｎ―ｍは、ｍ列ｎ番目の電池モジュールである。従って、単電池の総数はｎ×ｍ×ｋ個である。
　各単電池１－１～１－ｋは、リチウムイオン二次電池等の繰り返し充放電が可能な蓄電デバイスであり、開回路電圧値から充電状態値が観測可能なものである。ニッケル水素電池、鉛蓄電池、電気二重層キャパシタ、リチウムイオンキャパシタなども本電力直貯蔵装置を構成する蓄電デバイスとして使用可能である。電力貯蔵装置１には他にブレーカ、電池監視装置等が設けられることがあるが、ここでは省略する。
　電力貯蔵装置全体の端子間電圧を総電圧Ｖａｌｌ、充放電電流の総電流をＩａｌｌとし、充電方向を正とする。なお、充放電中は端子間の接続に用いられる導体やケーブルの抵抗成分に起因する電圧がＶａｌｌに加算される。
　通電していない充放電前の充電状態値を求めるため、電流ゼロの時の電圧を開回路電圧とする。リチウムイオン二次電池などの蓄電デバイスでは、一定温度環境下において、開回路電圧とＳＯＣは一対一の関係にあり、ＳＯＣは開回路電圧に関して一般に単調増加関数を示す。
　以下、電力貯蔵装置１を用いた電力貯蔵システムを説明した上で、本実施の形態にかかる充電状態推定装置の構成を説明する。
　図２は、本発明の実施の形態にかかる充電状態推定装置が適用される電力貯蔵システムの構成図である。電力貯蔵装置１は、コントローラ２からの指令に基づいて電流制御装置３によって制御される。このコントローラ２は、Ａ／Ｄコンバータ２０３とデータバンク２０２とパラメータ演算部２０１とを有して構成されている。なお、コントローラ２は、本実施の形態にかかる充電状態推定装置である。
　Ａ／Ｄコンバータ２０３は、電力貯蔵システム内の総電圧センサ４ａ、総電流検出センサ４ｂ、および温度センサ４ｃから送られるアナログ信号と、電力貯蔵装置１内の各セル電圧センサ（図示せず）から送られるアナログ信号４ｄとをデジタル信号に変換する。
　データバンク２０２には、電池データを格納すると共に、開回路電圧とＳＯＣの関係を記述した式、あるいはこの関係を示すデータテーブルなどが格納されている。開回路電圧とＳＯＣの関係は、使用する蓄電デバイスの種類、正極、負極の電極材料種によって異なるため、予めインプットするか測定で求める。求める方法としては、一定の温度に管理された環境において、既知の定電流を一定時間通電させてＳＯＣを一定間隔で変化させたときの開回路電圧データを取得する。開回路電圧は、電流ゼロの状態を数時間程度保持した後に測定することが望ましい。また、開回路電圧は、電力貯蔵装置１をシステムに組み込んだ状態から求めることも可能であり、制御プログラムの先頭で実行し、または長期の未使用期間、メンテナンス期間を利用してプログラムの実行をマニュアルで行うことも可能である。
　開回路電圧とＳＯＣとの関係を求める方法には、完全放電状態から求める方法と、満充電状態から求める方法とがある。完全放電状態から求める方法は、電池の下限電圧まで放電した完全放電時の充電状態値をゼロとして、一定ＳＯＣ間隔で充電してＳＯＣに対する開回路電圧を求める方法である。また、満充電状態からの方法は、電池の上限電圧まで定電流定電圧充電で充電した状態のＳＯＣを１００％として一定ＳＯＣ間隔で放電し、ＳＯＣを変化させる方法である。両者を比較すると充電状態から測定した方が安定した状態を得られ、データの精度は高い。
　電流制御装置３は、コントローラ２内のパラメータ演算部２０１からの指令によって電力貯蔵装置１の充放電量を制御する。電力が不足する場合、電流制御装置３は、交流グリッド７から電力変換装置６を介して供給される電力と、電力貯蔵装置１に貯蔵された電力と、を合わせてモータ５へ出力する。また、電流制御装置３は、モータ／発電機５で発電された電力を電力貯蔵装置１に貯蔵し、電力貯蔵装置１の容量をオーバーする電力に関しては、電力変換装置６を介して交流グリッド７に回生させる。なお、交流グリッド７は、例えば、鉄道車両に電力を供給する交流電力供給網である。
　パラメータ演算部２０１は、図２に示すように、第１の演算部３０、第２の演算部３１、および補正演算部３２を有して構成されている。第１の演算部３０は、電池容量（Ａｈ）と電池充放電時の通電電気量の積算値とに基づいて第１の充電状態推定値（以下単に「（第１のＳＯＣ）」と称する）を演算する。第２の演算部３１は、抵抗成分とコンデンサ成分を有する等価回路モデルから算出した開回路電圧と、開回路電圧と充電状態推定値の関係と、を用いて第２の充電状態推定値（以下単に「（第２のＳＯＣ）」と称する）を演算する。具体的には、第２の演算部３１は、電流ゼロの状態において充電状態推定値を推定する第２の推定処理と、電圧が変化する定電流制御時において開回路電圧から充電状態推定値を推定する第３の推定処理と、一定の電圧で制御する定電圧制御時において開回路電圧から充電状態推定値を推定する第４の推定処理とを実行する。補正演算部３２は、第１の演算部３０で推定されたＳＯＣを第２の演算部３１で推定されたＳＯＣで定期的に補正する。
　次に、図３を用いて、パラメータ演算部２０１で実行されるＳＯＣの推定動作を概説する。図３は、本発明の実施の形態にかかる充電状態推定装置によるＳＯＣ推定処理を説明するためのフローチャートである。図３に示されるフローチャートは、データを取得する間隔と同じ時間間隔で繰り返される。この時間間隔は電流、電圧、温度の各データの取得インターバルを意味し、ミリ秒間隔から数分間隔が適当である。
　データバンク２０２には、例えば、１つ前の計算フローで得たＳＯＣ（Ｎ－１）、電子・イオン抵抗Ｒ０（Ｎ－１）、電極反応抵抗Ｒ（Ｎ－１）、および静電容量成分Ｃ（Ｎ－１）などが記録されている（ステップＳ１０）。
　第２の演算部３１は、まず、定電流制御時であるか定電圧制御時であるかによって第２～第４の推定処理方法を切り替える。例えば、第２の演算部３１は、電流ゼロの開回路状態では（ステップＳ１２、Ｎｏ～Ｓ１３、Ｙｅｓ）、第２の推定処理を実行する（ステップＳ１４）。また、第２の演算部３１は、電圧が変化する定電流制御時には（ステップＳ１２、Ｎｏ～Ｓ１３、Ｎｏ）、第３の推定処理を実行する（ステップＳ１５）。さらに、第２の演算部３１は、一定の電圧で制御する定電圧制御時には（ステップＳ１２、Ｙｅｓ）、第４の推定処理を実行する（ステップＳ１６）。補正演算部３２は、第１のＳＯＣを第２のＳＯＣで補正する。
　データバンク２０２には、第１の演算部３０および第２の演算部３１で演算されたＳＯＣ（Ｎ）、電子・イオン抵抗Ｒ０（Ｎ）、電極反応抵抗Ｒ（Ｎ）、および静電容量成分Ｃ（Ｎ）が記録される（ステップＳ１８）。
　次に、図４～図１１を用いて、第１の演算部３０および第２の演算部３１の構成および動作を詳細に説明する。
　まず、第１の演算部３０の構成および動作を説明する。図４は、本実施の形態にかかる第１の演算部３０の構成図である。第１の演算部３０は、主たる構成として、平均電流算出部２１１、通電電気量算出部２１２、ＳＯＣ変化量算出部２１３、およびＳＯＣ算出部２１４を有して構成されている。
　電力貯蔵装置１を流れる総電流Ｉａｌｌは、平均電流算出部２１１に入力され、平均電流算出部２１１では、総電流Ｉａｌｌに所定のゲインを乗じて電流平均値を求める。すなわち、総電流Ｉａｌｌを電池モジュール１１－１～１ｎ－ｍの並列数ｍで割ることで、電流の平均値を得ている。通電電気量算出部２１２は、演算周期をΔｔとして、充電および放電時の通電電流を積分し、一定時間通電後の電気量を算出する。ＳＯＣ変化量算出部２１３は、通電電気量算出部２１２で得られた電気量（クーロン）を電池容量（Ａｈ）および３６００（ｓ）で除して、１００を乗じて変化量ΔＳＯＣ（％）を算出する。
　ＳＯＣ算出部２１４は、充電時には１つ前の計算フローで得られたＳＯＣ（Ｎ－１）に変化量ΔＳＯＣを加算し、放電時にはＳＯＣ（Ｎ－１）から変化量ΔＳＯＣ減算することで充電状態値の今回値ＳＯＣ（Ｎ）を得る。なお、ＳＯＣ（Ｎ－１）は、データバンク２０２に格納されており、ＳＯＣ算出部２１４は、データバンク２０２からのＳＯＣ（Ｎ－１）を用いてＳＯＣ（Ｎ）を推定する。これはＳＯＣを求める最も簡単な方法であるが、電流測定値の誤差を含むことが多い。また、電池の長期使用による劣化によって、ＳＯＣの１００％容量値の低下も誤差要因となるため、データバンク２０２に電池容量（Ａｈ）の現在値情報を変更入力する。なお、演算周期Δｔは、パラメータ演算部２０１の内部で生成されたものでもよいし、パラメータ演算部２０１の外部で生成されたものでもよい。
　次に、第２の演算部３１に関して説明する。図５は、本実施の形態にかかる第２の演算部３１に適用可能な蓄電デバイスの分布定数系等価回路モデルを示す図である。第２の演算部３１は、等価回路モデルに基づいてＳＯＣ、抵抗、コンデンサの静電容量を推定する。厳密には、第２の演算部３１は、図５に示す分布定数系等価回路に基づいて離散化した数値モデルに、電流の計測値および電圧の計測値をフィッティングさせることで開回路電圧を求める。計算式には、図５に示される分布定数系等価回路モデルの負極端子の抵抗８ａ、正極端子の抵抗８ｂ、負極電極層の電子抵抗９ａ、正極電極層の電子抵抗９ｂ、負極界面のコンデンサの静電容量１０ａ、正極界面のコンデンサの静電容量１０ｂ、負極界面の抵抗成分１１ａ、正極界面の抵抗成分１１ｂ、負極界面に発生する電位差１２ａ、正極界面に発生する電位差１２ｂ、負極内電解質の抵抗１３ａ、正極内電解質の抵抗１３ｂ、およびセパレータ内電解液の抵抗１４が含まれる。実際に計測される電流値および電圧の時間変化実測データと、計算値と、をフィッティングさせることで複数の抵抗値、静電容量値、および開回路電圧を求める。
　負極界面のコンデンサの静電容量１０ａおよび正極界面のコンデンサの静電容量１０ｂは、電極活物質と電解液界面に形成される電気二重層によるものであり、ファラッドの単位で表される。負極界面のコンデンサの静電容量１０ａおよび正極界面のコンデンサの静電容量１０ｂは、電極活物質の表面積に比例し、電解液の性状や電極電位によっても変化する。
　負極界面の抵抗成分１１ａおよび正極界面の抵抗成分１１ｂは、イオンから電子、電子からイオンへ電荷担体が変化する時の抵抗であり、電極と電解液界面に発生する抵抗であり、電荷移動抵抗と拡散抵抗を含む。
　このような分布定数系等価回路モデルによれば、厳密な計算が可能であるが、モデル規模が大きくなるというデメリットがあるため、計算周期を長くとる必要が生じ、電流、電圧の速い変化に対する応答が困難になる。また、多くの計算リソースが必要になるという問題も発生する。なお、図５に示される分布定数系等価回路モデルは、図１に示される１つの単電池（例えば１－１）を示すものであり、この分布定数系等価回路モデルに流れる電流Ｉは、総電流Ｉａｌｌを電池モジュール１１－１～１ｎ－ｍの並列数で除した値である。
　一方、計算周期を短縮して計算リソースが節約できる集中定数系等価回路モデルを図６に示す。図６は、本実施の形態にかかる第２の演算部３１に適用可能な蓄電デバイスの集中定数系等価回路モデルを示す図である。集中定数系等価回路モデルは、図５の分布定数系等価回路モデルを簡素化したものであり、この集中定数系等価回路モデルに流れる電流Ｉは、上述同様に、総電流Ｉａｌｌを電池モジュール１１－１～１ｎ－ｍの並列数で除した値である。このモデルにおける各要素は、以下のように構成されている。すなわち、図６に示される集中定数系等価回路モデルは、電極反応に関する抵抗１５（以下単に「抵抗１５」と称する）、電極界面に発生するコンデンサの静電容量１８（以下単に「静電容量１８」と称する）、電子・イオンに関する抵抗１６（以下単に「抵抗１６」と称する）、開回路電圧に相当する起電力部１７とを有し、抵抗１５と静電容量１８は並列に接続されている。
　抵抗１５の値はＲ、静電容量１８の静電容量値はＣ、抵抗１６の値はＲ０（電子・イオン抵抗）とする。抵抗１５を流れる電流をＩ１とし、静電容量１８を流れる電流をＩ２とすると、これらの和は電流Ｉである。コンデンサ（静電容量１８）にかかる電圧は、電流Ｉ１が流れる抵抗１５の両端電圧に等しく、またコンデンサ（静電容量１８）に蓄えられる電荷Ｑの時間変化が電流Ｉ２に相当することから、（１）式の電荷Ｑに関する微分方程式が得られる。
Figure JPOXMLDOC01-appb-M000001
　計算周期をΔｔとし、（１）式を電荷Ｑに関して離散化すると（２）式が得られる。時刻Ｎにおける電荷Ｑ（Ｎ）は、一つ前の時刻（Ｎ－１）における電荷Ｑ（Ｎ－１）と電流Ｉとを用いることで、（２）式のように表すことができる。
Figure JPOXMLDOC01-appb-M000002
　また、コンデンサ（静電容量１８）を流れる電流Ｉ２は、（３）式に示す電荷Ｑの時間変化で表すことができる。
Figure JPOXMLDOC01-appb-M000003
　（１）式と、（２）式と、電池電圧Ｖとを用いて、（４）式から開回路電圧Ｖｏｃを計算する。
Figure JPOXMLDOC01-appb-M000004
　図７は、本実施の形態にかかる第２の演算部３１の構成および動作を説明するための図である。図７－１は、第２の演算部３１とデータバンク２０２との関係を示す図であり、図７－２は、第２の演算部３１による第２の推定処理の動作を説明するための図である。
　図７－１に示される第２の演算部３１は、電圧変化判定部２２１と、通電判定部２２２と、電荷／電流値算出部２２３と、開回路電圧算出部２２４と、ＳＯＣ算出部２１４と、抵抗値算出部３６とを有して構成され、さらに抵抗値算出部３６は、劣化量算出部３５と、拡散種濃度算出部２２５と、電流値算出部２２６と、拡散抵抗算出部２２７とを有して構成されている。
　電圧変化判定部２２１は、演算周期Δｔにおける電池電圧Ｖの変化の有無を判定する。すなわち、電圧変化判定部２２１は、図３におけるステップＳ１２の判定を行う。
　通電判定部２２２は、電池電圧Ｖが変化している場合（ステップＳ１２、Ｎｏであるとき）、電流Ｉがゼロであるか否かを判定する。すなわち、通電判定部２２２は、図３におけるステップＳ１３の判定を行う。
　電荷／電流値算出部２２３は、計算の初期値としてコンデンサ１８に蓄電されていた電荷Ｑ（Ｎ－１）と、電極反応抵抗Ｒ（Ｎ－１）と、静電容量値Ｃ（Ｎ－１）と、電流Ｉと、演算周期Δｔとから、（２）式を用いて電荷Ｑ（Ｎ）を求めると共に、（３）式を用いて電流Ｉ２を求める。なお、電荷Ｑ（Ｎ－１）、電極反応抵抗Ｒ（Ｎ－１）、および静電容量値Ｃ（Ｎ－１）は、データバンク２０２に格納されているものとする。さらに、電荷／電流値算出部２２３は、Ｉ＝Ｉ１＋Ｉ２の関係から電流Ｉ１を求める。
　開回路電圧算出部２２４は、電池電圧Ｖと、データバンク２０２に格納されている抵抗１６の値Ｒ０（Ｎ－１）とから（４）式を用いて開回路電圧Ｖｏｃを算出する。なお、抵抗１６の値Ｒ０（Ｎ－１）は、データバンク２０２に格納されているものとする。
　ＳＯＣ算出部２１４は、開回路電圧算出部２２４からの開回路電圧Ｖｏｃを取り込み、データバンク２０２に格納されている開回路電圧ＶｏｃとＳＯＣとの関係を用いて、ＳＯＣ（Ｎ）を推定する。
　以下、第２の演算部３１による第２の推定処理および第３の推定処理を説明する。通電状態から電流ゼロのとき、分極緩和現象によって大きく電圧が変化し、その後徐々にある一定電圧に漸近する。開回路電圧算出部２２４は、電池電圧Ｖと、データバンク２０２に格納されている抵抗１６の値Ｒ０（Ｎ－１）とに基づいて、（４）式から電池電圧Ｖ、抵抗１５、抵抗１６、および電流Ｉをゼロとして、開回路電圧Ｖｏｃを求める。ＳＯＣ算出部２１４は、開回路電圧算出部２２４からの開回路電圧Ｖｏｃを取り込み、データバンク２０２に格納されている開回路電圧ＶｏｃとＳＯＣとの関係を用いて、第２のＳＯＣを推定する。
　電流値一定のとき、開回路電圧算出部２２４は、（４）式から電池電圧Ｖ、抵抗１５、抵抗１６、および電流Ｉに基づいて開回路電圧Ｖｏｃを求める。ＳＯＣ算出部２１４は、上述同様に、開回路電圧ＶｏｃとＳＯＣとの関係を用いて第２のＳＯＣを推定する。なお、電流を流さずにＳＯＣを求める第２の推定処理の方が精度は高い。
　次に、図８および図９を用いて、劣化量算出部３５の構成および動作を説明する。図８は、本実施の形態にかかる劣化量算出部３５の構成図である。図８－１は、劣化量算出部３５の構成を示すものであり、図８－２は、電圧シミュレーション部２３１および抵抗／容量算出部２３２の動作を説明するための図である。また、図９は、抵抗とコンデンサの静電容量との積と、電池容量との関係を説明するための図である。
　図８－１に示される劣化量算出部３５は、主たる構成として、電圧シミュレーション部２３１と抵抗／容量算出部２３２を有して構成されている。劣化量算出部３５は、ある一定時間の電流、電圧データから抵抗１５、抵抗１６、静電容量１８の変化によって電池の劣化状態を推定するものである。具体的には、劣化量算出部３５は、想定される制御モードにおける抵抗成分、コンデンサの静電容量成分を逐次算出することにより、電池劣化状態の進行度を推定し、例えば、複数の電圧データと計算値をフィッティングさせることで、抵抗１５、抵抗１６、静電容量１８の値それぞれを、電極反応抵抗Ｒ（Ｎ）、電子・イオン抵抗Ｒ０（Ｎ）、静電容量成分Ｃ（Ｎ）として特定する。
　フィッティングさせる電圧データは、時間ｔ１からｔｐまでのｐ個であり、個々のデータ時間間隔はΔｔである。参照するデータ数が多いほど計算精度は向上するが、計算負荷を考慮すると１０個以内のデータ数が好ましい。以下、実際の計算手順を示す。電圧シミュレーション部２３１は、時間ｔ１の１つ前のコンデンサ１８に蓄えられた電荷をＱ（ｋ－１）とし、電極反応抵抗Ｒ（Ｎ）、電子・イオン抵抗Ｒ０（Ｎ）、および静電容量成分Ｃ（Ｎ）のそれぞれに直前の値を代入し、電流Ｉおよび開回路電圧Ｖｏｃに図７－１のＳＯＣ算出時に使用した値を代入して、（５）式を用いて時間ｔ１に対する電圧Ｖｋを算出する。なお、（５）式内の電荷Ｑ（Ｎ）は（２）式を用いて計算する。
Figure JPOXMLDOC01-appb-M000005
　同様に、電圧シミュレーション部２３１は、ｔ２からｔｐまでを計算する。
　抵抗／容量算出部２３２は、電圧シミュレーション部２３１で算出されたこれらの計算値（電圧Ｖｋ）と実測データ（Ｖｐ）とを比較して、実測値と計算値の差の総和ΣΔＶが、判定値εとデータ数ｐとの積未満の場合には一致と見なす。一方、実測値と計算値の差の総和ΣΔＶが、判定値εとデータ数ｐとの積以上の場合、抵抗／容量算出部２３２は、電極反応抵抗Ｒ（Ｎ）、電子・イオン抵抗Ｒ０（Ｎ）、静電容量成分Ｃ（Ｎ）の値を変更する。例えば、抵抗／容量算出部２３２は、放電電圧が計算値より低ければ電極反応抵抗Ｒ（Ｎ）、電子・イオン抵抗Ｒ０（Ｎ）、静電容量成分Ｃ（Ｎ）を大きく、逆に高ければ小さくする。
　フィッティングパラメータには、抵抗１５と静電容量１８との積を用いると便利である。抵抗／容量算出部２３２は、この積と劣化した電池の容量との関係（図９参照）を用いて、電池容量を算出する。ただし、電池劣化は、急激に進行する反応ではないため、劣化の計算をデータ取得周期で行う必要はなく、起動直後または停止直前に１回の頻度で十分である。
　次に、第４の推定処理に関して説明する。定電圧制御時の抵抗１５は、時間経過とともに上昇して一定の値をとらない。抵抗１５の値が上量する理由は、ネルンストアインシュタインの式である（６）式に示されるように、この抵抗値が電極活物質内のリチウムイオン濃度ＣＬｉ（Ｎ）の逆数に比例するためであり、定電圧制御では可動リチウムイオン数の減少に起因する。
Figure JPOXMLDOC01-appb-M000006
　ここで、ＤＬｉはリチウムイオンの活物質内の拡散定数、Ｔはモジュール温度、Ａは比例定数である。リチウムイオン濃度ＣＬｉ（Ｎ）の時間変化は、拡散方程式である（７）式を含む物質移動式を数値的に解き、抵抗１５を算出する。
Figure JPOXMLDOC01-appb-M000007
　以下、図１０および図１１を用いて抵抗値算出部３６の構成および動作を説明する。図１０は、本実施の形態にかかる抵抗値算出部３６の動作を説明するための図であり、図１１は、定電圧制御時における抵抗の変化を説明するための図である。
　拡散種濃度算出部２２５は、定電圧制御に入る直前の活物質内の可動リチウムイオン濃度であるＣＬｉ（Ｎ－１）と拡散定数ＤＬｉとから、（７）式を境界条件の元で解いて、演算周期Δｔ時間後のリチウムイオン濃度ＣＬｉ（Ｎ）を算出する。
　電流値算出部２２６は、電池電圧Ｖの変化から電流値を逆算する。拡散抵抗算出部２２７は、電流値算出部２２６で算出された電流値と、定電流制御時における電圧と、に基づいて、電圧が一定ではない場合（ＣＣ充電： constant current充電）における抵抗１５の値を推定する。さらに拡散抵抗算出部２２７は、（６）式を用いて、電圧が一定の場合（ＣＶ充電：constant voltage充電）における抵抗１５の値を推定する。
　電極反応に関する抵抗１５の値は、ＣＣ充電時においては変化しないが、ＣＶ充電時においては図１１のように増加する。この抵抗増加は、定電圧制御条件下における数秒から数時間の間に発生する一時的なものである。第２の演算部３１は、電池の劣化による抵抗増加をリアルタイムに評価するときに、図１１に示すような定電圧時の抵抗増加を、容量劣化推定のための計算から除外するという処理も実行する。
　また、ここで算出した短期的な抵抗の増加速度、すなわち、図１１の抵抗値上昇量を定電圧制御時間で除した値は、物質移動に関わる電池の劣化を示している。このように、本実施の形態にかかる第２の演算部３１は、劣化の進行と共にこの値が大きくなることを利用して、長期的劣化による電池寿命を予測する。
　図１２は、劣化した電池におけるＳＯＣ演算精度を示す図である。比較例１は、通電電気量の積算による第１の演算部３０で算出された第１のＳＯＣの推移である。計算に用いる電流値には、一般に誤差が含まれるので、ＳＯＣが真値からずれることとなる。
　比較例２は、この誤差を、等価回路モデルによる従来のＳＯＣ推定処理で補ったものであり、フレッシュな初期の電池においては、良好なＳＯＣ推定値が得られる。しかしながら、長期間使用して劣化が進行したリチウムイオン電池のＳＯＣ推定においては、大きな誤差が生じる。
　実施例１は、本実施の形態にかかる充電状態推定装置によって推定されたＳＯＣの推移を示すものである。
　充放電開始点１９と充放電終了点２０においては、開回路で３日以上放置して求めた正しいＳＯＣである。ここで、比較例１によるＳＯＣ推定方法では、電流値が実際より充電時に小さく検知積算され、放電時に大きく検知積算されて、充放電全体にわたり真値を下まわる結果となる。一方、比較例２によるＳＯＣ推定方法では、劣化による容量低下を考慮しないままＳＯＣを算出するため、充電、放電ともにＳＯＣの真値より大きく算出されるため、過剰に変化することになる。
　実施例１によるＳＯＣ推定方法では、劣化量算出部３５によって、電池劣化による容量低下を考慮していることに加え、抵抗値算出部３６によって、定電圧時の抵抗増加を容量劣化推定のための計算から除外するようにしたので、ＳＯＣ演算精度を向上させることが可能である。
　以上に説明したように、本実施の形態にかかる充電状態推定装置は、電池容量と、充電状態推定値の前回値ＳＯＣ（Ｎ－１）と、電力貯蔵装置１の充放電量を制御する電流制御装置３と電力貯蔵装置１との間で流出入する電流と、に基づいて演算した充電状態推定値の今回値を第１の充電状態推定値として算出する第１の演算部３０と、定電流制御時には、電池の等価回路モデルと電池の電圧とに基づいて演算した充電状態推定値を第２の充電状態推定値の今回値として算出する一方で、定電圧制御時には、電池の等価回路モデルと電池の電圧とに基づいて、電池の抵抗変化を考慮して演算した充電状態推定値の今回値を第２の充電状態推定値として算出する第２の演算部３１と、第１の充電状態推定値を第２の充電状態推定値で定期的に補正する補正演算部３２と、を備えるようにしたので、電池容量と電池充放電時の通電電気量の積算値とに基づく第１のＳＯＣと、電流ゼロの状態における第２のＳＯＣと、定電流制御時における第３のＳＯＣと、定電圧制御時における第４のＳＯＣとを推定することができ、従来技術に比して電池の充電状態および劣化状態を精度よく推定することが可能である。
　また、第２の演算部３１で用いられる電池の等価回路モデルは、コンデンサと非線形抵抗の並列回路を１つ有してなる集中定数系等価回路モデル、または、コンデンサと非線形抵抗の並列回路を複数有してなる分布定数系等価回路モデルで構成するようにしたので、集中定数系等価回路モデルを用いて計算リソースが節約することもできるし、分布定数系等価回路モデルを用いて厳密な計算を行うことも可能である。
　また、本実施の形態にかかる第２の演算部３１は、コンデンサに蓄えられた電荷Ｑ（Ｎ－１）と、コンデンサの静電容量値Ｃ（Ｎ－１）と、電池の抵抗値（Ｒ（Ｎ－１）、Ｒ０（Ｎ－１））と、電池の電圧Ｖと、に基づいて第２の充電状態推定値を算出するようにしたので、簡単な構成で第２の充電状態推定値を推定可能である。
　また、第２の演算部３１は、電流が変化した後に計測された連続する複数の電圧計測値（Ｖｐ）と、電荷Ｑ（Ｎ－１）と静電容量値Ｃ（Ｎ－１）と抵抗値Ｒ（Ｎ－１）、Ｒ０（Ｎ－１）とに基づいて算出された電圧計算値（Ｖｋ）と、の差の総和ΣΔＶが所定の値（判定値εとデータ数ｐとの積）になるまで静電容量値Ｃ（Ｎ）および抵抗値Ｒ（Ｎ）、Ｒ０（Ｎ）を計算し、電圧計測値（Ｖｐ）と電圧計算値（Ｖｋ）と差の総和ΣΔＶが所定の値と一致したときの静電容量Ｃ（Ｎ）と抵抗値Ｒ（Ｎ）、Ｒ０（Ｎ）とを、電池の劣化状態を推定する指標として出力する劣化量算出部３５を有するようにしたので、従来技術に比して、等価回路モデルによる充電状態値を精度よく推定することが可能である。
　また、静電容量Ｃ（Ｎ）と抵抗値Ｒ（Ｎ）、Ｒ０（Ｎ）との積に基づいて電池容量を算出するようにしたので、第１の演算部３０は、劣化量算出部３５で算出された電池容量を用いて第１の充電状態推定値を精度よく算出することが可能である。
　また、第２の演算部３１は、可動リチウムイオン濃度ＣＬｉ（Ｎ－１）とリチウムイオンの活物質内の拡散定数ＤＬｉとに基づいて、定電圧制御時における電池の抵抗値Ｒ（Ｎ）を算出する抵抗値算出部３６を有するようにしたので、定電圧時の抵抗増加を容量劣化推定のための計算から除外することができ、ＳＯＣ演算精度を向上させることが可能である。
　以上のように、本発明は、二次電池等の蓄電池におけるＳＯＣを推定する充電状態推定装置に適用可能であり、特に、ＳＯＣの推定精度を向上させることができる発明として有用である。
　１　電力貯蔵装置
　１－１、１－２、１－ｋ　単電池
　２　コントローラ
　３　電流制御装置
　４ａ　総電圧センサ
　４ｂ　総電流検出センサ
　４ｃ　温度センサ
　４ｄ　電圧センサから送られるアナログ信号
　５　モータ／発電機
　６　電力変換装置
　７　交流グリッド
　８ａ　負極端子の抵抗
　８ｂ　正極端子の抵抗
　９ａ　負極電極層の電子抵抗
　９ｂ　正極電極層の電子抵抗
　１０ａ　負極界面のコンデンサの静電容量
　１０ｂ　正極界面のコンデンサの静電容量
　１１－１、１２－１、１ｎ－１、１１－２、１２－２、１ｎ－２、１１―ｍ、１２―ｍ、１ｎ―ｍ　電池モジュール
　１１ａ　負極界面の抵抗成分
　１１ｂ　正極界面の抵抗成分
　１２ａ　負極界面に発生する電位差
　１２ｂ　正極界面に発生する電位差
　１３ａ　負極内電解質の抵抗
　１３ｂ　正極内電解質の抵抗
　１４　セパレータ内電解液の抵抗
　１５　電極反応に関する抵抗
　１６　電子・イオンに関する抵抗
　１７　開回路電圧に相当する起電力部
　１８　電極界面に発生するコンデンサの静電容量
　１９　充放電開始点
　２０　充放電終了点
　３０　第１の演算部
　３１　第２の演算部
　３２　補正演算部
　３５　劣化量算出部
　３６　抵抗値算出部
　２０１　パラメータ演算部
　２０２　データバンク部
　２０３　Ａ／Ｄコンバータ
　２１１　平均電流算出部
　２１２　通電電気量算出部
　２１３　ＳＯＣ変化量算出部
　２１４　ＳＯＣ算出部
　２２１　電圧変化判定部
　２２２　通電判定部
　２２３　電荷／電流値算出部
　２２４　開回路電圧算出部
　２２５　拡散種濃度算出部
　２２６　電流値算出部
　２２７　拡散抵抗算出部
　２３１　電圧シミュレーション部
　２３２　抵抗／容量算出部
　Ｉａｌｌ　総電流
　Ｔ　モジュール温度
　Ｖａｌｌ　総電圧
　Ｖｏｃ　開回路電圧

Claims (8)

  1. 　電池を複数接続した電力貯蔵装置に接続され、前記電力貯蔵装置の残存容量を示す充電状態値を推定する充電状態推定装置であって、
    　電池容量と、前記充電状態推定値の前回値と、前記電力貯蔵装置の充放電量を制御する電流制御装置と前記電力貯蔵装置との間で流出入する電流と、に基づいて演算した充電状態推定値の今回値を第１の充電状態推定値として算出する第１の演算部と、
    　前記電力貯蔵装置を一定の電流で充電する定電流制御時には、前記電池の等価回路モデルと前記電池の電圧とに基づいて演算した充電状態推定値を第２の充電状態推定値の今回値として算出する一方で、前記電力貯蔵装置を一定の電圧で充電する定電圧制御時には、前記電池の等価回路モデルと前記電池の電圧とに基づいて、電池の抵抗変化を考慮して演算した充電状態推定値の今回値を第２の充電状態推定値として算出する第２の演算部と、
    　前記第１の充電状態推定値を第２の充電状態推定値で定期的に補正する補正演算部と、
    　を備えたことを特徴とする充電状態推定装置。
  2. 　前記電池の等価回路モデルは、コンデンサと非線形抵抗の並列回路を１つ有してなる等価回路モデル、または、コンデンサと非線形抵抗の並列回路を複数有してなる等価回路モデルで構成されていること、を特徴とする請求項１に記載の充電状態推定装置。
  3. 　前記第２の演算部は、前記コンデンサに蓄えられた電荷と、前記コンデンサの静電容量値と、前記電池の抵抗値と、前記電池の電圧とに基づいて前記第２の充電状態推定値を算出すること、を特徴とする請求項１に記載の充電状態推定装置。
  4. 　前記第２の演算部は、電流が変化した後に計測された連続する複数の電圧計測値と、前記電荷と前記静電容量値と前記抵抗値とに基づいて算出された電圧計算値と、の差の総和が所定の値になるまで前記静電容量値および前記抵抗値を計算し、前記電圧計測値と前記電圧計算値と差の総和が前記所定の値と一致したときの前記コンデンサの静電容量と前記電池の抵抗値とを、電池の劣化状態を推定する指標として出力する劣化量算出部を有すること、を特徴とする請求項１に記載の充電状態推定装置。
  5. 　前記劣化量算出部は、前記コンデンサの静電容量と前記電池の抵抗値との積に基づいて電池容量を算出すること、を特徴とする請求項４に記載の充電状態推定装置。
  6. 　前記第２の演算部は、可動リチウムイオン濃度とリチウムイオンの活物質内の拡散定数とに基づいて、前記定電圧制御時における電池の抵抗値を算出する抵抗算出部を有すること、を特徴とする請求項１に記載の充電状態推定装置。
  7. 　前記第２の演算部は、前記劣化量算出部で算出された静電容量および抵抗値を用いて第２の充電状態推定値を算出すること、を特徴とする請求項１に記載の充電状態推定装置。
  8. 　前記第２の演算部は、前記拡散抵抗算出部で算出された電池の抵抗値を用いて前記第２の充電状態推定値を算出すること、を特徴とする請求項１に記載の充電状態推定装置。
     
PCT/JP2010/059626 2010-06-07 2010-06-07 充電状態推定装置 WO2011155017A1 (ja)

Priority Applications (10)

Application Number Priority Date Filing Date Title
KR1020127027258A KR101338639B1 (ko) 2010-06-07 2010-06-07 충전 상태 추정 장치
US13/640,056 US8975897B2 (en) 2010-06-07 2010-06-07 State-of-charge estimating apparatus
EP10852857.1A EP2579059B1 (en) 2010-06-07 2010-06-07 Charge status estimation apparatus
BR112012027318A BR112012027318A2 (pt) 2010-06-07 2010-06-07 aparelho de estimativa do estado de carga
JP2012519152A JP5511951B2 (ja) 2010-06-07 2010-06-07 充電状態推定装置
CN201080067268.7A CN102918411B (zh) 2010-06-07 2010-06-07 充电状态推定装置
RU2012155314/28A RU2565339C2 (ru) 2010-06-07 2010-06-07 Устройство оценки состояния заряда
MX2012012939A MX2012012939A (es) 2010-06-07 2010-06-07 Aparato de estimacion de estado de carga.
PCT/JP2010/059626 WO2011155017A1 (ja) 2010-06-07 2010-06-07 充電状態推定装置
AU2010354957A AU2010354957B2 (en) 2010-06-07 2010-06-07 Charge status estimation apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2010/059626 WO2011155017A1 (ja) 2010-06-07 2010-06-07 充電状態推定装置

Publications (1)

Publication Number Publication Date
WO2011155017A1 true WO2011155017A1 (ja) 2011-12-15

Family

ID=45097649

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2010/059626 WO2011155017A1 (ja) 2010-06-07 2010-06-07 充電状態推定装置

Country Status (10)

Country Link
US (1) US8975897B2 (ja)
EP (1) EP2579059B1 (ja)
JP (1) JP5511951B2 (ja)
KR (1) KR101338639B1 (ja)
CN (1) CN102918411B (ja)
AU (1) AU2010354957B2 (ja)
BR (1) BR112012027318A2 (ja)
MX (1) MX2012012939A (ja)
RU (1) RU2565339C2 (ja)
WO (1) WO2011155017A1 (ja)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012251806A (ja) * 2011-05-31 2012-12-20 Toshiba Corp 算出方法、算出プログラム、算出システム、算出装置
JP2014035283A (ja) * 2012-08-09 2014-02-24 Denso Corp 電源システム
US20140139230A1 (en) * 2012-11-19 2014-05-22 Infineon Technologies Ag Method and apparatus for determining a charge state
CN104115017A (zh) * 2012-02-17 2014-10-22 丰田自动车株式会社 电池系统以及劣化判别方法
EP2762907A4 (en) * 2012-10-26 2016-01-20 Lg Chemical Ltd DEVICE AND METHOD FOR ESTIMATING THE LOADING CONDITION OF A BATTERY
EP2952919A4 (en) * 2012-12-04 2016-12-14 Lg Chemical Ltd METHOD AND DEVICE FOR ESTIMATING PARAMETERS FOR SECONDARY BATTERY
US9525300B2 (en) 2013-01-18 2016-12-20 Toyota Jidosha Kabushiki Kaisha Electricity storage system
JPWO2015008620A1 (ja) * 2013-07-15 2017-03-02 古河電気工業株式会社 二次電池状態検出装置および二次電池状態検出方法
JP2017195698A (ja) * 2016-04-20 2017-10-26 スズキ株式会社 電池管理装置および電池管理方法
JP2018040676A (ja) * 2016-09-07 2018-03-15 日立化成株式会社 車載用蓄電デバイスモデルのパラメータ同定方法及びパラメータ同定装置
US10209319B2 (en) 2013-02-01 2019-02-19 Toyota Jidosha Kabushiki Kaisha State of deterioration or state of charges estimating apparatus for secondary battery
JP2019100897A (ja) * 2017-12-05 2019-06-24 日立化成株式会社 デバイス状態検知装置、電源システムおよび自動車
JP2019146302A (ja) * 2018-02-16 2019-08-29 株式会社半導体エネルギー研究所 二次電池の充電状態推定装置及び異常検出装置、及び二次電池の管理システム
WO2020208762A1 (ja) * 2019-04-11 2020-10-15 三菱電機株式会社 蓄電池の特性推定装置および特性推定方法
US10807494B2 (en) 2016-03-09 2020-10-20 Vehicle Energy Japan Inc. Battery management system, battery system and hybrid vehicle control system
JP7487125B2 (ja) 2021-02-10 2024-05-20 株式会社豊田自動織機 燃料電池システム

Families Citing this family (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012143996A1 (ja) * 2011-04-18 2012-10-26 日立ビークルエナジー株式会社 蓄電装置
US9213070B2 (en) * 2011-05-31 2015-12-15 Kabushiki Kaisha Toshiba Calculation method, calculation system, and calculation apparatus
JP5291845B1 (ja) * 2012-01-26 2013-09-18 カルソニックカンセイ株式会社 電池の状態推定装置
KR101327591B1 (ko) * 2012-07-02 2013-11-12 엘에스산전 주식회사 인버터의 직류링크 커패시터 진단장치
WO2014024226A1 (ja) * 2012-08-09 2014-02-13 トヨタ自動車株式会社 蓄電システム
US9182449B2 (en) * 2012-10-12 2015-11-10 GM Global Technology Operations LLC Method and system for estimating battery capacity in a vehicle
US9841464B2 (en) * 2012-12-26 2017-12-12 Mitsubishi Electric Corporation Life prediction apparatus for electrical storage device and life prediction method for electrical storage device
CA2897702C (en) * 2013-03-18 2017-11-21 Mitsubishi Electric Corporation Apparatus and method for estimating power storage device electric energy
WO2014167602A1 (ja) * 2013-04-11 2014-10-16 トヨタ自動車株式会社 電池システム
US9347998B2 (en) * 2013-04-17 2016-05-24 Allegro Microsystems, Llc System and method for measuring battery voltage
KR102248599B1 (ko) 2014-05-20 2021-05-06 삼성에스디아이 주식회사 배터리의 충전방법 및 이를 위한 배터리 관리 시스템
US9193273B1 (en) * 2014-06-15 2015-11-24 Efficient Drivetrains, Inc. Vehicle with AC-to-DC inverter system for vehicle-to-grid power integration
US9987942B2 (en) * 2014-09-03 2018-06-05 Ford Global Technologies, Llc Method of operating vehicle powertrain based on prediction of how different chemical type batteries connected in parallel will operate to output demanded current
JP6262686B2 (ja) * 2015-04-27 2018-01-17 ファナック株式会社 平滑コンデンサの寿命予測手段を有するモータ制御装置
KR102032229B1 (ko) * 2015-11-18 2019-10-16 현대일렉트릭앤에너지시스템(주) 배터리 수명 상태 추정 시스템 및 그 방법
JP6647111B2 (ja) 2016-03-29 2020-02-14 古河電気工業株式会社 二次電池劣化推定装置および二次電池劣化推定方法
JP6339618B2 (ja) * 2016-03-29 2018-06-06 古河電気工業株式会社 二次電池劣化推定装置および二次電池劣化推定方法
JP5980459B1 (ja) * 2016-03-30 2016-08-31 本田技研工業株式会社 電源装置、該電源装置を有する輸送機器、蓄電部の充電率と開放端電圧の相関情報を推定する推定方法、および該相関情報を推定するためのプログラム
CN105807168B (zh) * 2016-05-11 2018-06-08 大连理工大学 一种用于修正超级电容器soc估算的方法
CN110167783B (zh) 2017-01-09 2022-09-20 沃尔沃卡车集团 一种用于确定电池组的充电状态的方法和装置
JP6690584B2 (ja) * 2017-03-10 2020-04-28 トヨタ自動車株式会社 電池状態推定装置
CN107271906B (zh) * 2017-05-31 2019-10-18 宁德时代新能源科技股份有限公司 电池包健康度估算方法和装置
KR102203245B1 (ko) 2017-11-01 2021-01-13 주식회사 엘지화학 배터리 soc 추정 장치 및 방법
KR20190107888A (ko) * 2018-03-13 2019-09-23 한국전자통신연구원 제로 에너지 타운 피크 전력 관리 방법 및 장치
HRP20230478T1 (hr) * 2018-03-28 2023-07-21 Toyo System Co., Ltd. Uređaj za određivanje stanja degradacije i postupak određivanja stanja degradacije
EP3579007B1 (en) * 2018-06-07 2022-08-17 Rolls-Royce Deutschland Ltd & Co KG Method and apparatus for estimating a state of charge of a battery
CN111239613B (zh) * 2018-11-29 2022-07-05 宏碁股份有限公司 电池电量估计方法与电子装置
CN113196523A (zh) * 2018-12-18 2021-07-30 松下知识产权经营株式会社 电池状态估计装置、电池状态估计方法以及电池系统
CN113013507B (zh) * 2019-12-19 2024-04-09 横河电机株式会社 二次电池管理系统及二次电池管理方法
JP7191873B2 (ja) * 2020-01-17 2022-12-19 株式会社東芝 充放電制御装置、充放電システム、充放電制御方法及び充放電制御プログラム
JP7387660B2 (ja) * 2021-02-10 2023-11-28 株式会社東芝 電池の診断方法、電池の診断装置、電池の診断システム、電池搭載機器及び電池の診断プログラム
JP2022139620A (ja) * 2021-03-12 2022-09-26 株式会社東芝 電池の診断方法、電池の診断装置、電池の診断システム、電池搭載機器及び電池の診断プログラム
US20230358815A1 (en) * 2022-05-09 2023-11-09 Samsung Electronics Co., Ltd. Method and device with estimating battery state

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08140270A (ja) 1994-11-04 1996-05-31 Mitsubishi Electric Corp 二次電池のパラメータ測定方法ならびにそれを用いた二次電池の充放電制御方法および寿命予測方法、ならびに、二次電池の充放電制御装置およびそれを用いた電力貯蔵装置
JP2002238106A (ja) * 2001-02-14 2002-08-23 Denso Corp ハイブリッド車の電池状態制御方法
JP2003075518A (ja) 2001-09-05 2003-03-12 Nissan Motor Co Ltd 二次電池の充電率推定装置
JP2008199723A (ja) 2007-02-09 2008-08-28 Railway Technical Res Inst バッテリーの残容量推定装置
JP2010071703A (ja) * 2008-09-17 2010-04-02 Calsonic Kansei Corp 組電池の残量推定方法及び装置

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1081499B1 (en) * 1998-05-28 2003-07-09 Toyota Jidosha Kabushiki Kaisha Means for estimating charged state of battery and method for estimating degraded state of battery
JP3879278B2 (ja) * 1998-11-10 2007-02-07 日産自動車株式会社 ハイブリッド車の充電量演算方法および充電量演算装置
DE10203810A1 (de) * 2001-06-29 2003-01-16 Bosch Gmbh Robert Verfahren zur Ermittlung des Ladezustands und/oder der Leistungsfähigkeit eines Ladungsspeichers
JP3714246B2 (ja) * 2001-12-18 2005-11-09 日産自動車株式会社 二次電池の充電率推定装置
JP3852371B2 (ja) * 2002-05-24 2006-11-29 日産自動車株式会社 二次電池の充電率推定装置
JP3689084B2 (ja) * 2002-12-11 2005-08-31 三菱電機株式会社 バッテリ充電状態演算装置およびバッテリ充電状態演算方法
US7321220B2 (en) * 2003-11-20 2008-01-22 Lg Chem, Ltd. Method for calculating power capability of battery packs using advanced cell model predictive techniques
JP2006129588A (ja) * 2004-10-28 2006-05-18 Sanyo Electric Co Ltd 二次電池の電力制御方法及び電源装置
CN101116003B (zh) * 2005-01-27 2010-04-07 松下电动车辆能源股份有限公司 二次电池充/放电电量、极化电压及soc的估计方法和装置
RU2283504C1 (ru) * 2005-06-28 2006-09-10 Федеральное государственное унитарное предприятие "Производственно-конструкторское предприятие "ИРИС" Автоматизированная система контроля и диагностики аккумуляторных батарей
JP4830382B2 (ja) * 2005-07-19 2011-12-07 日産自動車株式会社 二次電池の充電率推定装置
JP4984527B2 (ja) * 2005-12-27 2012-07-25 トヨタ自動車株式会社 二次電池の充電状態推定装置および充電状態推定方法
KR100859688B1 (ko) * 2006-10-12 2008-09-23 삼성에스디아이 주식회사 배터리 관리 시스템 및 그의 구동 방법
JP4706648B2 (ja) * 2007-03-06 2011-06-22 トヨタ自動車株式会社 電動車両、充電状態推定方法および充電状態推定方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読取可能な記録媒体
JP4703593B2 (ja) * 2007-03-23 2011-06-15 株式会社豊田中央研究所 二次電池の状態推定装置
JP5036416B2 (ja) * 2007-06-15 2012-09-26 トヨタ自動車株式会社 電源システムおよびそれを備えた車両、ならびに充放電制御方法
JP4544273B2 (ja) * 2007-06-20 2010-09-15 トヨタ自動車株式会社 車両用電源装置および車両用電源装置における蓄電装置の充電状態推定方法
US7830119B2 (en) * 2007-08-29 2010-11-09 Gm Global Technology Operations, Inc. Adaptive battery estimator and method
US7957921B2 (en) * 2008-02-19 2011-06-07 GM Global Technology Operations LLC Model-based estimation of battery hysteresis
JP5358990B2 (ja) * 2008-03-21 2013-12-04 株式会社明電舎 電力変換装置
JP2010019595A (ja) * 2008-07-08 2010-01-28 Fuji Heavy Ind Ltd 蓄電デバイスの残存容量演算装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08140270A (ja) 1994-11-04 1996-05-31 Mitsubishi Electric Corp 二次電池のパラメータ測定方法ならびにそれを用いた二次電池の充放電制御方法および寿命予測方法、ならびに、二次電池の充放電制御装置およびそれを用いた電力貯蔵装置
JP2002238106A (ja) * 2001-02-14 2002-08-23 Denso Corp ハイブリッド車の電池状態制御方法
JP2003075518A (ja) 2001-09-05 2003-03-12 Nissan Motor Co Ltd 二次電池の充電率推定装置
JP2008199723A (ja) 2007-02-09 2008-08-28 Railway Technical Res Inst バッテリーの残容量推定装置
JP2010071703A (ja) * 2008-09-17 2010-04-02 Calsonic Kansei Corp 組電池の残量推定方法及び装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2579059A4 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012251806A (ja) * 2011-05-31 2012-12-20 Toshiba Corp 算出方法、算出プログラム、算出システム、算出装置
CN104115017A (zh) * 2012-02-17 2014-10-22 丰田自动车株式会社 电池系统以及劣化判别方法
JP2014035283A (ja) * 2012-08-09 2014-02-24 Denso Corp 電源システム
EP2762907A4 (en) * 2012-10-26 2016-01-20 Lg Chemical Ltd DEVICE AND METHOD FOR ESTIMATING THE LOADING CONDITION OF A BATTERY
US20140139230A1 (en) * 2012-11-19 2014-05-22 Infineon Technologies Ag Method and apparatus for determining a charge state
EP2952919A4 (en) * 2012-12-04 2016-12-14 Lg Chemical Ltd METHOD AND DEVICE FOR ESTIMATING PARAMETERS FOR SECONDARY BATTERY
US9720047B2 (en) 2012-12-04 2017-08-01 Lg Chem, Ltd. Apparatus and method for estimating parameter of secondary battery
US9525300B2 (en) 2013-01-18 2016-12-20 Toyota Jidosha Kabushiki Kaisha Electricity storage system
US10209319B2 (en) 2013-02-01 2019-02-19 Toyota Jidosha Kabushiki Kaisha State of deterioration or state of charges estimating apparatus for secondary battery
JPWO2015008620A1 (ja) * 2013-07-15 2017-03-02 古河電気工業株式会社 二次電池状態検出装置および二次電池状態検出方法
US10807494B2 (en) 2016-03-09 2020-10-20 Vehicle Energy Japan Inc. Battery management system, battery system and hybrid vehicle control system
JP2017195698A (ja) * 2016-04-20 2017-10-26 スズキ株式会社 電池管理装置および電池管理方法
JP2018040676A (ja) * 2016-09-07 2018-03-15 日立化成株式会社 車載用蓄電デバイスモデルのパラメータ同定方法及びパラメータ同定装置
JP2019100897A (ja) * 2017-12-05 2019-06-24 日立化成株式会社 デバイス状態検知装置、電源システムおよび自動車
JP7038530B2 (ja) 2017-12-05 2022-03-18 昭和電工マテリアルズ株式会社 デバイス状態検知装置、電源システムおよび自動車
JP2019146302A (ja) * 2018-02-16 2019-08-29 株式会社半導体エネルギー研究所 二次電池の充電状態推定装置及び異常検出装置、及び二次電池の管理システム
WO2020208762A1 (ja) * 2019-04-11 2020-10-15 三菱電機株式会社 蓄電池の特性推定装置および特性推定方法
JPWO2020208762A1 (ja) * 2019-04-11 2021-10-28 三菱電機株式会社 蓄電池の特性推定装置および特性推定方法
CN113661399A (zh) * 2019-04-11 2021-11-16 三菱电机株式会社 蓄电池的特性推测装置以及特性推测方法
JP7069412B2 (ja) 2019-04-11 2022-05-17 三菱電機株式会社 蓄電池の特性推定装置および特性推定方法
CN113661399B (zh) * 2019-04-11 2024-04-19 三菱电机株式会社 蓄电池的特性推测装置以及特性推测方法
JP7487125B2 (ja) 2021-02-10 2024-05-20 株式会社豊田自動織機 燃料電池システム

Also Published As

Publication number Publication date
CN102918411B (zh) 2015-02-25
RU2565339C2 (ru) 2015-10-20
MX2012012939A (es) 2013-01-17
EP2579059A1 (en) 2013-04-10
EP2579059A4 (en) 2013-05-29
AU2010354957A1 (en) 2012-12-13
CN102918411A (zh) 2013-02-06
AU2010354957B2 (en) 2014-04-17
JPWO2011155017A1 (ja) 2013-08-01
BR112012027318A2 (pt) 2016-08-02
EP2579059B1 (en) 2014-04-02
US20130027047A1 (en) 2013-01-31
RU2012155314A (ru) 2014-07-20
JP5511951B2 (ja) 2014-06-04
US8975897B2 (en) 2015-03-10
KR20120139818A (ko) 2012-12-27
KR101338639B1 (ko) 2013-12-06

Similar Documents

Publication Publication Date Title
JP5511951B2 (ja) 充電状態推定装置
Lai et al. Co-estimation of state of charge and state of power for lithium-ion batteries based on fractional variable-order model
CN110914696B (zh) 用于在电池的操作期间估计电池开路池格电压、充电状态以及健康状态的方法和系统
Wang et al. Modeling and state-of-charge prediction of lithium-ion battery and ultracapacitor hybrids with a co-estimator
Xiong et al. A robust state-of-charge estimator for multiple types of lithium-ion batteries using adaptive extended Kalman filter
Wei et al. Adaptive estimation of state of charge and capacity with online identified battery model for vanadium redox flow battery
Wang et al. Online dynamic equalization adjustment of high-power lithium-ion battery packs based on the state of balance estimation
CN107076801B (zh) 电池管理系统中的电化学阻抗谱
Li et al. A comparative study of state of charge estimation algorithms for LiFePO4 batteries used in electric vehicles
Wang et al. Online state of charge estimation for the aerial lithium-ion battery packs based on the improved extended Kalman filter method
US20190288344A1 (en) Battery Control Device
KR102572652B1 (ko) 배터리의 충전상태를 추정하는 방법
Wang et al. Battery available power prediction of hybrid electric vehicle based on improved Dynamic Matrix Control algorithms
Jantharamin et al. A new dynamic model for lead-acid batteries
Qiu et al. Battery hysteresis modeling for state of charge estimation based on Extended Kalman Filter
CN110632520A (zh) 一种动力电池soc的估算装置及其估算方法
KR101268082B1 (ko) 분극전압과 개로전압을 이용한 배터리 잔존용량 추정방법
CN114503392A (zh) 涉及多个电池的判定装置、蓄电系统、判定方法和判定程序
EP3920308A1 (en) Battery control device
US20240044996A1 (en) Systems, methods, and devices for power rating estimation in energy storage systems
Lystianingrum et al. State of health and life estimation methods for supercapacitors
JP7174327B2 (ja) 二次電池の状態判定方法
JP2024501600A (ja) バッテリの充電状態および健全状態を推定する方法およびシステム
Kustiman et al. Battery state of charge estimation based on coulomb counting combined with recursive least square and pi controller
dos Santos et al. Electric Vehicles Batteries Modeling Analysis Based on a Multiple Layered Perceptron Identification Approach

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 201080067268.7

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 10852857

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2012519152

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 13640056

Country of ref document: US

ENP Entry into the national phase

Ref document number: 20127027258

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: MX/A/2012/012939

Country of ref document: MX

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 10269/CHENP/2012

Country of ref document: IN

ENP Entry into the national phase

Ref document number: 2010354957

Country of ref document: AU

Date of ref document: 20100607

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2010852857

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2012155314

Country of ref document: RU

Kind code of ref document: A

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112012027318

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 112012027318

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20121024