WO2014050073A1 - 電池状態推定装置及び蓄電池システム - Google Patents

電池状態推定装置及び蓄電池システム Download PDF

Info

Publication number
WO2014050073A1
WO2014050073A1 PCT/JP2013/005628 JP2013005628W WO2014050073A1 WO 2014050073 A1 WO2014050073 A1 WO 2014050073A1 JP 2013005628 W JP2013005628 W JP 2013005628W WO 2014050073 A1 WO2014050073 A1 WO 2014050073A1
Authority
WO
WIPO (PCT)
Prior art keywords
timing
voltage
battery
unit
parameter
Prior art date
Legal status (The legal status 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 status listed.)
Ceased
Application number
PCT/JP2013/005628
Other languages
English (en)
French (fr)
Inventor
尚人 右田
石井　洋平
晶彦 山田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sanyo Electric Co Ltd
Original Assignee
Sanyo Electric Co Ltd
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 Sanyo Electric Co Ltd filed Critical Sanyo Electric Co Ltd
Priority to JP2014538171A priority Critical patent/JP6158195B2/ja
Priority to US14/428,665 priority patent/US9846200B2/en
Publication of WO2014050073A1 publication Critical patent/WO2014050073A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

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/3644—Constructional arrangements
    • G01R31/3648—Constructional arrangements comprising digital calculation means, e.g. for performing an algorithm
    • 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/382—Arrangements for monitoring battery or accumulator variables, e.g. 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/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
    • G01R31/3835—Arrangements for monitoring battery or accumulator variables, e.g. SoC involving only voltage measurements
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
    • H01M10/00—Secondary cells; Manufacture thereof
    • H01M10/42—Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
    • H01M10/48—Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10—Energy storage using batteries

Definitions

  • the present invention relates to a battery state estimation device and a storage battery system.
  • the open-circuit voltage of the storage battery can be detected by measuring the terminal voltage of the storage battery when charging / discharging is stopped.
  • the impedance circuit in the internal equivalent circuit of the storage battery also includes a capacity component, the terminal voltage fluctuates transiently for a while (for example, several minutes to several hours) after charging or discharging stops.
  • a conventional method minimizes the coefficient of the approximate expression including a plurality of exponential terms using the voltage value within a predetermined time after charging and discharging are stopped. It calculates using a square method (refer the following patent document 1).
  • the stable open-circuit voltage (open-circuit voltage convergence value) may be estimated by the above-described conventional method, fitting using multiple exponential functions is necessary, so that the parameter setting becomes complicated and the computation load is reduced. heavy.
  • an object of the present invention is to provide a battery state estimation device and a storage battery system capable of accurately estimating a stable open-circuit voltage (open-circuit voltage convergence value) by a simple calculation.
  • the battery state estimation device includes a measurement value acquisition unit that sequentially acquires a measurement value of a battery voltage, which is a terminal voltage of a storage battery, and a reference timing in a non-energization period after the storage battery is charged or discharged.
  • a timing setting unit that sets a second timing that varies only by using the measurement value, a first time difference between the reference timing and the first timing, and a second time difference between the first timing and the second timing, And an open-circuit voltage estimating unit that estimates an open-circuit voltage of the storage battery at a target timing after the second timing.
  • the present invention it is possible to provide a battery state estimation device and a storage battery system capable of accurately estimating a stable open-circuit voltage (open-circuit voltage convergence value) by a simple calculation.
  • 1 is a schematic overall configuration diagram of a storage battery system according to an embodiment of the present invention. It is an internal equivalent circuit schematic of the battery module which concerns on embodiment of this invention. It is a graph which shows the time transition of the terminal voltage and internal voltage of a battery module concerning embodiment of this invention. It is a figure which shows the time fluctuation of the battery voltage during a non-energization period.
  • 1 is an internal block diagram of an open-circuit voltage / charge rate estimation apparatus according to a first embodiment of the present invention. It is an operation
  • FIG. 1 is a schematic overall configuration diagram of a storage battery system 1 according to an embodiment of the present invention.
  • the storage battery system 1 includes a battery module 11, a current sensor 12, a voltage sensor 13, a temperature sensor 14, and a battery control unit 15. You may think that the electric power block 20 is also contained in the component of the storage battery system 1.
  • FIG. 1 is a schematic overall configuration diagram of a storage battery system 1 according to an embodiment of the present invention.
  • the storage battery system 1 includes a battery module 11, a current sensor 12, a voltage sensor 13, a temperature sensor 14, and a battery control unit 15. You may think that the electric power block 20 is also contained in the component of the storage battery system 1.
  • FIG. 1 is a schematic overall configuration diagram of a storage battery system 1 according to an embodiment of the present invention.
  • the storage battery system 1 includes a battery module 11, a current sensor 12, a voltage sensor 13, a temperature sensor 14, and a battery control unit 15. You may think that the electric power block 20 is also contained in the component of the storage battery system 1.
  • the battery module 11 includes one or more storage batteries (in other words, secondary batteries).
  • the storage battery forming the battery module 11 is any type of storage battery, for example, a lithium ion battery or a nickel metal hydride battery.
  • the battery module 11 is formed by a plurality of storage batteries connected in series, but the number of storage batteries forming the battery module 11 may be one. Some or all of the storage batteries included in the battery module 11 may be connected in parallel to each other.
  • the battery module 11 is considered as one storage battery.
  • discharging and charging means discharging and charging of the battery module 11 unless otherwise specified.
  • FIG. 2 is an equivalent circuit in the battery module 11.
  • the equivalent circuit in the battery module 11 can be considered as a series connection circuit of a voltage source VS that outputs the voltage V S and has an internal resistance of zero, and an impedance circuit Z that includes a resistance component and a capacitance component.
  • the power block 20 is connected to the battery module 11.
  • the power block 20 includes a load and a power source.
  • the battery module 11 can supply discharge power to a load in the power block 20 and can be supplied with charging power from a power source in the power block 20.
  • a power conversion circuit (not shown) may be interposed between the battery module 11 and the load and power source in the power block 20.
  • the current sensor 12 is interposed between the battery module 11 and the power block 20 and measures a battery current I which is a current flowing through the battery module 11.
  • the voltage sensor 13 measures a battery voltage V that is a terminal voltage of the battery module 11 (that is, a potential difference between the positive electrode and the negative electrode of the battery module 11).
  • the temperature sensor 14 measures a battery temperature TMP that is the temperature of the battery module 11.
  • the battery temperature TMP is, for example, the surface temperature of the pack that encloses the storage battery in the battery module 11 or the temperature at a specific part in the battery module 11.
  • the symbols I, V, and TMP are also used as symbols indicating the value of the battery current, the value of the battery voltage, and the value of the battery temperature, respectively.
  • the battery current value I, the battery voltage value V, and the battery temperature TMP measured by the sensors 12, 13, and 14 may be referred to as the measured current value I, the measured voltage value V, and the measured temperature TMP, and the description is simplified.
  • the measured values I, V, and T may be referred to for the sake of simplicity.
  • the battery control unit 15 controls charging and discharging of the battery module 11 using battery state data including the measurement values I, V, and TMP.
  • the battery control unit 15 includes a battery state estimation device 30 that estimates an open voltage of the battery module 11 using the battery state data.
  • the battery module 11 Prior to the description of the battery state estimation device 30, a situation as shown in FIG.
  • the battery module 11 is not charged and discharged in the period P0 from the timing t A to the timing t B , and then the battery module 11 is discharged in the period P1 from the timing t B to the timing t C , and thereafter Assume that charging and discharging of the battery module 11 are stopped in the period P2 from the timing t C to the timing t D, and discharging of the battery module 11 is performed again in the period P3 after the timing t D.
  • Each of the periods P1 and P3 is an energization period in which the battery module 11 is discharged or charged.
  • Each of the periods P0 and P2 is a non-energization period in which the battery module 11 is not discharged and charged.
  • charging or discharging of the battery module 11 is also expressed as energization of the battery module 11.
  • the energization period P1 is a discharge period during which discharge is performed, but the energization period P1 may be a charge period during which charging is performed, and in this case, the timing t C is a charging stop timing. .
  • a solid line waveform 301 represents a time transition of the battery voltage V
  • a broken line waveform 302 represents a time transition of the voltage V S (see FIG. 2).
  • the battery voltage V corresponds to the voltage V S minus the voltage drop due to the resistance component of the impedance circuit Z.
  • the battery voltage V gradually increases over a corresponding time (for example, several minutes to several hours) as shown by the waveform 301 due to the presence of the impedance circuit Z, and during the period P2. Converges to a constant stable open-circuit voltage V S —P 2, which is the voltage V S.
  • FIG. 4 is an enlarged view of the waveform 301 during the non-energization period P2.
  • Timing t 0 ⁇ t n including timing t 0 ⁇ t 4 shows a plurality of times between the timing t C and t D (n is an integer).
  • the timing t i + 1 is later than the timing t i, representing the battery voltage V at the timing t i at symbol V i.
  • N is a unit quantity (unit voltage quantity) having a positive predetermined value, for example, several tens of millivolts.
  • Timing t 0 functions as a reference timing.
  • a time length T and a voltage variation parameter (hereinafter also simply referred to as a parameter) R are defined.
  • t i + 1 ⁇ t i represents a time length between timings t i and t i + 1 .
  • T t 1 ⁇ t 0 (A1)
  • R (t 2 ⁇ t 1 ) / (t 1 ⁇ t 0 ) (A2)
  • the voltage variation parameter R is represented by a ratio of the second time length “t 2 -t 1 ” to the first time length “t 1 -t 0 ”.
  • the first time length “t 1 ⁇ t 0 ” is the time length required for the battery voltage V to change by the unit amount N as viewed from the reference timing t 0
  • the second time length “t 2 ⁇ t”. 1 " the time when the battery voltage V as viewed from the reference timing t 0 varies by a unit amount N (starting ie timing t 1) as a starting point, the time required for the battery voltage V is varied by further unit quantity N Length.
  • the energization period P1 before the non-energization period P2 is a discharge period
  • the fluctuation of the battery voltage V during the non-energization period P2 is “increase”. If the energization period P1 is a charging period, the fluctuation of the battery voltage V during the non-energization period P2 is “decrease”.
  • the battery voltage V fluctuates toward the stable open circuit voltage V S — P2 (see FIG. 3), but the fluctuation speed is related to the fact that the impedance circuit Z includes a capacitance component. It gradually decreases with time.
  • the voltage variation parameter R has a value corresponding to the rate of change (decrease rate) of the variation rate of the battery voltage V during the non-energization period P2.
  • t represents an arbitrary timing during the non-energization period P2.
  • t ⁇ t 0 represents an elapsed time from the reference timing t 0 to the timing t.
  • ⁇ V represents the fluctuation amount of the battery voltage V between the timings t 0 and t.
  • V EST (t) represents the battery voltage V at the timing t, which can be estimated by the battery state estimation device 30.
  • the battery state estimation device 30 models the behavior of the battery voltage V after the timing t2 by the equation (A4), and the battery voltage V at the timing t at the time before the timing t based on the estimation equation (A4) " V EST (t) "can be estimated (predicted).
  • V EST (t) a future timing t viewed from the current timing (for example, t 2 ) at which this estimation is performed is referred to as a target timing
  • V EST (t) is also referred to as an estimated open voltage (predicted open voltage).
  • FIG. 5 is an internal block diagram of the open-circuit voltage / charge rate estimation device 50 included in the battery state estimation device 30.
  • the estimation apparatus 50 according to the first embodiment includes each part referred to by reference numerals 51 to 56.
  • the measurement value acquisition unit 51 periodically acquires the measurement current value I, the measurement voltage value V, and the measurement temperature TMP by sampling the output signals of the sensors 12, 13, and 14 at a predetermined sampling period. However, acquisition of the measured temperature TMP can be omitted.
  • the timer 52 measures and acquires the current time, or measures an elapsed time from an arbitrary time. Each part in the battery state estimation device 30 including the estimation device 50 can freely use the acquisition result of the acquisition unit 51 and the acquisition / measurement result of the timer 52.
  • the timing setting unit 53 sets the reference timing, the first timing, and the second timing using the measured voltage value V during the non-energization period after the energization of the battery module 11 is stopped.
  • the non-energization period includes a voltage fluctuation period in which the battery voltage V fluctuates transiently, and the reference timing, the first timing, and the second timing are set during the voltage fluctuation period. Also in the first embodiment, assuming a period P2 as the de-energized period should be carried out estimation of open circuit voltage, it is assumed that the timing t 0 is set to the reference timing (the same applies to other embodiments described later).
  • the timing setting unit 53 sets the reference timing t 0 based on the elapsed time after the energization of the battery module 11 is stopped. More specifically, for example, the time point when a predetermined time has elapsed from the boundary timing t C between the periods P1 and P2 (the timing when the discharge is stopped) is set as the reference timing t 0 .
  • the timing setting unit 53 may set the reference timing t 0 based on the voltage fluctuation amount of the battery voltage V after the energization of the battery module 11 is stopped. More specifically, for example, the measured voltage value V is a predetermined voltage amount when viewed from the voltage value V [t C ], which is the measured voltage value V at the boundary timing t C (timing when discharge is stopped) between the periods P1 and P2. The timing increased by only the reference timing t 0 is set. Alternatively, after the timing t C , the fluctuation amount of the battery voltage V per unit time may be monitored as the monitoring target amount, and the timing at which the monitoring target amount decreases and becomes a predetermined value or less may be set as the reference timing t 0. .
  • the timing setting unit 53 further includes a first timing at which the battery voltage V varies by a predetermined first voltage amount as viewed from the reference timing t 0 and a predetermined first voltage at which the battery voltage V is greater than the first voltage amount as viewed from the reference timing.
  • a second timing that fluctuates by two voltage amounts is set using measurement values V that are sequentially acquired.
  • the first voltage amount is equal to the unit amount N described above, and the second voltage amount is twice the unit amount N (see FIG. 4). Then, the first and second timings are the timings t 1 and t 2 in FIG. 4 , respectively.
  • the parameter deriving unit 54 is based on the time difference (t 1 -t 0 ) between the reference and the first timing set by the timing setting unit 53 and the time difference (t 2 -t 1 ) between the first and second timings.
  • the voltage fluctuation parameter R is derived according to the above equation (A2).
  • Open voltage calculating unit 55 at a timing earlier than and target timing t a second timing t 2 after using the voltage fluctuation parameter R derived by deriving unit 54, the target timing in accordance with the above formula (A4) Estimate the open circuit voltage V EST (t) at t.
  • the timing setting unit 53 holds the measured voltage value V at the reference timing t 0 as the reference voltage value V 0 , and the calculation unit 55 substitutes the held reference voltage value V 0 into the equation (A4). .
  • the value of “T” in the formula (A4) is supplied from the timing setting unit 53 to the calculation unit 55. Further, the calculation unit 55 can recognize the elapsed time (t ⁇ t 0 ) from the reference timing t 0 to the prediction timing t, which is substituted into the equation (A4), using the timer 52.
  • Target timing t is any timing of the timing t 2 after.
  • the target timing t may be a timing at which a predetermined time (for example, one hour or several hours) has elapsed from the timing t C , t 0, or t 2 .
  • a predetermined time for example, one hour or several hours
  • the battery voltage V generally reaches the stable open-circuit voltage V S — P ⁇ b> 2 after about 1 hour has passed since the energization stop, so the timing one hour after the timing t C is set to the target timing t. If set, a practical estimated open circuit voltage V EST (t) is obtained, and V EST (t) functions as an estimated value of the stable open circuit voltage V S — P2.
  • the open circuit voltage V EST (t) is estimated using the expression (A4). Is equivalent to estimating the open circuit voltage of the battery module 11 assuming that the above-mentioned assumption ASM is correct (assuming the proposition ASM is true).
  • the charge rate estimation unit 56 charges the battery module 11 from the estimated open voltage V EST (t) using a predetermined arithmetic expression or table data indicating the relationship between the open voltage of the battery module 11 and the charge rate of the battery module 11. Estimate the rate.
  • the charging rate of the battery module 11 is a ratio of the remaining capacity of the battery module 11 to the full charge capacity of the battery module 11, and is generally called SOC (State of Charge).
  • FIG. 6 shows an operation flowchart relating to the open-circuit voltage estimation process.
  • step S11 periodic acquisition of measurement values I, V, and TMP by the measurement acquisition unit 51 (for example, sequential acquisition with a period of 1 millisecond) is started. However, the acquisition of the measured value TMP can be omitted.
  • the timing setting unit 53 monitors the measurement value I to determine whether or not the energization (charging or discharging) of the battery module 11 that has been performed so far has been stopped, and it is confirmed that the energization has been stopped.
  • the reference timing t 0 is set according to the method described above, and the measured voltage value at the reference timing t 0 is held as the reference voltage value V 0 (steps S12 and S13). Further, thereafter, the timing setting unit 53 compares the sequentially obtained measurement value V with the reference voltage value V 0 , whereby the measurement value V changes from the reference voltage value V 0 by a first voltage amount (for example, N). And the timing at which the measured value V fluctuates from the reference voltage value V 0 by a second voltage amount (for example, 2 ⁇ N) is extracted and set as the first timing (t 1 ) and the second timing (t 2 ), respectively ( Step S14).
  • a voltage fluctuation parameter R is derived from the set reference, first and second timings (step S15), and the open circuit voltage V EST (t) is estimated using the parameter R and the charging rate is estimated (step S15). S16).
  • the parameter R necessary for open-circuit voltage estimation can be obtained only by observing the fluctuation amount of the battery voltage V. Therefore, the stable open-circuit voltage (open-circuit voltage convergence value) V S — P2 can be accurately calculated by simple calculation. Can be estimated.
  • a solid line waveform V REAL1 and a broken line waveform V EST1 represent experimental waveforms of the measured voltage V and the estimated open circuit voltage V EST (t) in P2 during a non-energization period under a certain condition. It can be seen that both waveforms are in good agreement and a good estimate is made.
  • the waveforms V REAL1 and V EST1 in FIG. 7 are obtained under certain temperature conditions (for example, when the battery temperature TMP is 25 ° C.). Depending on the battery temperature TMP, as shown in FIG. The error between the voltage V EST (t) and the actual battery voltage V can be large.
  • the solid line waveform V REAL2 and the broken line waveform V EST2 are the measured voltage V and the estimated open-circuit voltage V EST (t) during the non-energization period under low temperature conditions (for example, the condition where the battery temperature TMP is 0 ° C.). ) Represents the experimental waveform.
  • the waveform V EST2 is a waveform of the estimated open circuit voltage V EST (t) when the parameter correction described later is not performed.
  • V EST the waveform of the estimated open circuit voltage V EST (t) when the parameter correction described later is not performed.
  • the second embodiment a configuration for reducing an open-circuit voltage estimation error will be described.
  • the behavior of the battery voltage V after stopping charging or discharging is affected by factors other than the battery temperature TMP, it is better to consider various factors that affect the behavior.
  • FIG. 9 is an internal block diagram of an open-circuit voltage / charge rate estimating apparatus 50A according to the second embodiment.
  • An estimation device 50A can be provided in the battery state estimation device 30 of FIG. By adding the parameter correction unit 57 to the estimation device 50 of FIG. 5, the estimation device 50A of FIG. 9 is formed. Except for the parameter correction unit 57, the configuration and operation of the estimation device 50A are the same as those of the estimation device 50 described in the first embodiment.
  • the parameter correction unit 57 corrects the voltage fluctuation parameter R derived by the parameter deriving unit 54 based on the correction index.
  • the corrected voltage fluctuation parameter R is represented by the symbol “ RC ”.
  • the open-circuit voltage calculation unit 55 uses the parameter R C as “R” in the equation (A4) and then opens the open-circuit voltage at the target timing t according to the equation (A4). Estimate V EST (t).
  • the correction index includes one or more of the following first to fourth indices.
  • the correction unit 57 may correct the voltage variation parameter R using two, three, or four indexes among the first to fourth indexes.
  • the first index is the battery temperature TMP.
  • the battery temperature TMP as the first index is a measured temperature TMP at a timing based on the time point when the energization of the battery module 11 is stopped. That is, for example, battery temperature TMP of the first index, the timing t C, the timing t C than the predetermined time before the timing, or a measured temperature TMP in the timing after a predetermined time than the timing t C.
  • the second index is the current rate (current amount per unit time) of energization of the battery module 11 during the energization period P1.
  • the third index is a change amount ⁇ SOC [P1] of the charging rate of the battery module 11 during the energization period P1.
  • the estimation device 50A of FIG. 9 or the battery state estimation device 30 of FIG. 1 integrates the measured current value I during the energization period P1, and the charge as the third index is calculated from the integrated value and the full charge capacity of the battery module 11.
  • the rate change amount ⁇ SOC [P1] can be obtained.
  • the fourth index is the charging rate SOC [P2] of the battery module 11 in the non-energization period P2.
  • the charging rate SOC [P2] of the battery module 11 in the non-energizing period P2 is equal to the charging rate of the battery module 11 at the end point of the energizing period P1 (that is, timing t C ).
  • the 1 can determine the charging rate SOC [P0] of the battery module 11 in the non-energization period P0 using the method described in the first embodiment ( That is, it can be obtained from the open circuit voltage estimated during the non-energization period P0), or can be obtained from the open circuit voltage (measured voltage V) measured during the non-energization period P0, and the charge rate SOC [P0] and the above
  • the charging rate SOC [P2] can be obtained using the change amount ⁇ SOC [P1].
  • the behavior of the battery voltage V during the non-energization period is actually investigated while changing the first to fourth indexes in various ways, and the parameter R is converted to the parameter RC according to the correction index using the survey results. It is preferable to create an arithmetic expression or table data to be performed, and hold the obtained arithmetic expression or table data in the parameter correction unit 57.
  • FIG. 10 is an operation flowchart relating to an open-circuit voltage estimation process according to the second embodiment.
  • the operation flowchart of FIG. 10 is the same as that of FIG. 6 except that a parameter correction step S20 is inserted between the parameter R deriving step S15 and the open circuit voltage estimating step S16.
  • the derived parameter R is corrected to the parameter RC (step S20), and the open circuit voltage V EST (t) is estimated using the parameter RC and the charging rate is estimated (step S16).
  • the voltage fluctuation parameter is corrected in consideration of the correction index that affects the behavior of the battery voltage V after the stop of charging or discharging, so that the estimation accuracy of the open-circuit voltage can be improved.
  • FIG. 11 is an internal block diagram of an open-circuit voltage / charge rate estimation apparatus 50B according to the third embodiment.
  • An estimation device 50B can be provided in the battery state estimation device 30 of FIG.
  • the estimation device 50B of FIG. 11 is formed by adding the measurement value holding unit 58 to the estimation device 50 of FIG. Except for the presence of the measurement value holding unit 58 and the contents described later corresponding to the presence, the configuration and operation of the estimation device 50B are the same as those of the estimation device 50 described in the first embodiment.
  • Measured value holding unit 58 after the end of the conduction period P1, the measured value of the battery voltage V is sequentially acquired by the measuring value acquisition unit 51, it holds only a predetermined time T Q min. That is, the measurement value holding unit 58 holds the measurement voltage value V acquired during the holding target period 320 of FIG. 12 in association with each measurement time (that is, the measured voltage value V of the holding target period 320 Hold time series data). Holding period 320, the timing t C, a period from the timing t C until time t CQ has elapsed the predetermined time T Q. Measured value holding unit 58 may together hold a measured voltage value V of the time period exceeding a predetermined time T Q.
  • the timing setting unit 53 does not set the first and second timings t 1 and t 2 with reference to the reference timing t 0 after setting the reference timing t 0.
  • the end timing t CQ of the period 320 is handled as the second timing t 2 (see FIG. 12), and then, based on each measured voltage value V held in the holding unit 58, the reference timing t 2 with the second timing t 2 as a reference. 0 and the first timing t 1 are set.
  • the timing setting unit 53 a second timing voltage value measured voltage value V at t 2 as a timing t CQ with defined as the second timing t 2 Treat as V 2 (see FIG. 4).
  • the setting unit 53 measures the measured voltage value V (that is, V 1 ) that is lower than the voltage value V 2 by the unit amount N and the measured voltage value V that is lower than the voltage value V 2 by twice the unit amount N (that is, V 0 ).
  • the operation after the timings t 0 , t 1 and t 2 are set is as described in the first embodiment.
  • the parameter correction unit 57 of FIG. 9 is added to the estimation device 50B of FIG. 11, and the correction of the parameter R described in the second embodiment and the open circuit voltage V EST (t) based on the corrected parameter R (ie, R C ). An estimation may be performed.
  • the estimated open circuit voltage at that time is set after the lapse of the predetermined time T Q is a timing t 2 is ensured.
  • ⁇ 4th Example A fourth embodiment will be described.
  • the parameter update unit 60 (see FIG. 13) according to the fourth embodiment can be mounted on the estimation device 50, 50A or 50B (FIG. 5, FIG. 9, or FIG. 11).
  • the update unit 60 The operation of the estimation device 50, 50A or 50B on which is mounted will be described.
  • the timing setting unit 53 sets the reference timing t 0 ′ during the non-energization period P2 in addition to the settings of the timings t 0 , t 1 and t 2 , and from the reference timing t 0 ′.
  • the battery voltage V changes by the second voltage amount (for example, 2N) as viewed from the first timing t 1 ′ when the battery voltage V changes by the first voltage amount (for example, N) and the reference timing t 0 ′.
  • Two timings t 2 ′ are set using the measurement values V acquired sequentially.
  • the reference timing t 0 ′ is an arbitrary timing after the reference timing t 0 and may coincide with the timing t 1 or t 2 or may be a predetermined timing after the timing t 2 .
  • the timings t 0 ′, t 1 ′, and t 2 ′ function as update timings.
  • the parameter derivation unit 54 uses the following formula (B1 ) To obtain the update parameter R ′ (see FIG. 14).
  • “t i + 1 ′ ⁇ t i ′” represents a time length between timings t i ′ and t i + 1 ′.
  • R ′ (t 2 ′ ⁇ t 1 ′) / (t 1 ′ ⁇ t 0 ′) (B1)
  • the update parameter R ′ may be corrected to the update parameter R C ′ based on the correction index as in the method of correcting the parameter R to the parameter RC (FIG. 14). reference).
  • the parameter update unit 60 obtains the parameter R UPDATE according to the following formula (B2) or formula (B3).
  • the weighting coefficient w in the formulas (B2) and (B3) may be 1, or may have an arbitrary predetermined value within a range satisfying “0 ⁇ w ⁇ 1”.
  • R UPDATE (1-w) ⁇ R + w ⁇ R ′ (B2)
  • R UPDATE (1-w) ⁇ R C + w ⁇ R C '(B3)
  • Open voltage calculating unit 55 before derivation of parameters R UPDATE, but to estimate the open circuit voltage V EST (t) based on the parameters R, or R C, after deriving the parameters R UPDATE, based on the parameter R UPDATE open Estimate the voltage V EST (t). That is, after deriving the parameter R UPDATE , the calculation unit 55 uses the parameter R UPDATE as “R” in the equation (A4), and then calculates the open circuit voltage V EST (t) at the target timing t according to the equation (A4). presume. As described above, the derivation of the parameter R UPDATE updates the voltage fluctuation parameter used for the open circuit voltage estimation (the parameter R UPDATE corresponds to the updated parameter).
  • the battery voltage V gradually converges to the stable open voltage V S — P2 while being influenced by various factors inside the battery module 11, but rather than the behavior of the battery voltage V in the initial period of the voltage fluctuation period, The behavior of the battery voltage V in the middle or end of the voltage fluctuation period better represents the state of convergence to the stable open-circuit voltage V S —P2. For this reason, the estimation accuracy of the open circuit voltage can be improved by updating the voltage fluctuation parameter as in the fourth embodiment.
  • the operation when the voltage fluctuation parameter is updated only once has been described. According to the same method as described above, a plurality of sets of timings t 0 ′, t 1 ′ and t 2 ′ are set, and the update operation is performed more than once. You can go.
  • FIG. 15 is a specific example of timings t 0 , t 1 and t 2 and timings t 0 ′, t 1 ′ and t 2 ′ assumed in the fifth embodiment.
  • the battery voltages V at the timings t 0 ′, t 1 ′, and t 2 ′ are V 0 ′, (V 0 ′ + N), and (V 0 ′ + 2N), respectively.
  • the parameter deriving unit 54 obtains the voltage fluctuation parameter R and R ′ based on the above (A2) and (B1), and then obtains the voltage fluctuation parameter R J according to the equation (B4).
  • the weighting factor w in the formula (B4) is an arbitrary predetermined value within a range satisfying “0 ⁇ w ⁇ 1”, but may be typically 0.5.
  • the parameter R J may be corrected to the parameter R JC based on the correction index in the same manner as the method of correcting the parameter R to the parameter R C.
  • the calculation unit 55 estimates the open circuit voltage V EST (t) at the target timing t according to the above equation (A4) after using the parameter R J or R JC as “R” in the equation (A4).
  • R J (1 ⁇ w) ⁇ R + w ⁇ R ′ (B4)
  • the fifth embodiment described above can be said to be a form in which the voltage variation parameter is set at a timing of four or more points. As described below, the voltage variation parameter may be set at a timing of four or more points.
  • the timing setting unit 53 sets the timings t 0 , t 1, and t 2 as described above, and further, the third timing at which the battery voltage V fluctuates by a predetermined third voltage amount as viewed from the reference timing t 0 , It sets using the measured value V acquired sequentially.
  • the third voltage amount is larger than the first and second voltage amounts described above.
  • the first and second voltage amounts are 1 and 2 times the unit amount N, respectively, and the third voltage amount is 3 times the unit amount. Then, the third timing coincides with the timing t 3 in FIG.
  • the parameter deriving unit 54 uses t 0 to t 3 as measured values, and calculates the value of the parameter R that accurately satisfies the following expressions (C1) to (C3) using the least square method or the like. To derive.
  • the operation after the derivation of the parameter R is as described in the above embodiments.
  • the parameter R is determined using the timings t 0 to t 3 of four points, but the parameter R may be determined using timings of five points or more (for example, t 0 to t 4 ). By increasing the number of timings for deriving the parameter R, the parameter R is made more appropriate, and the open-circuit voltage estimation accuracy is improved.
  • the variation amount of the battery voltage V from the reference timing t 0 to the i-th timing t i is i times the unit amount N (in other words, the battery voltage between the timings t i and t i + 1 ).
  • the fluctuation amount of the battery voltage V between the timings t i and t i + 1 is a logarithmic function or an exponential function depending on the integer i. It may be expressed. For example, when the fluctuation amount of the battery voltage V between the timings t 0 and t 1 is the unit amount N, the fluctuation amount of the battery voltage V between the timings t 1 and t 2 is “logN” and “N m ”. (M is a predetermined value).
  • the parameter deriving unit 54 and the open-circuit voltage calculating unit 55 regard the assumption ASM as correct (the proposition ASM is true) and derive the parameter R and the open-circuit voltage V EST (t). (The estimation formula of the open circuit voltage V EST (t) can be changed from the above formula (A4)).
  • the estimated open circuit voltage V EST (t) When the voltage fluctuation parameter does not depend on time and does not change, the estimated open circuit voltage V EST (t) according to the above formula (A4) continues to fluctuate with time. Even in that case, if the target timing t is appropriately set (for example, if the timing one hour after the timing t C is set as the target timing t as described in the first embodiment), the estimated open circuit voltage V EST (t) represents the stable open circuit voltage V S — P2 with high accuracy. However, in order to converge the estimated open circuit voltage V EST (t) to the stable open circuit voltage V S — P2, the voltage variation parameter (R, R C , R UPDATE , R J or R JC ) is made a function of time (exponential function, etc.). May be. That is, for example, after derivation of the voltage fluctuation parameter, the value of the voltage fluctuation parameter may be increased with time as shown in FIG.
  • an estimation formula is formed using a logarithmic function, but the form of the open-circuit voltage estimation formula is arbitrary, and a function other than the logarithmic function (for example, an exponential function) ) May be used to form an estimation formula for the open circuit voltage.
  • the target device that is the battery control unit 15, the battery state estimation device 30, or the estimation device 50, 50A, or 50B can be configured by hardware or a combination of hardware and software.
  • Arbitrary specific functions that are all or part of the functions realized in the target device are described as a program, the program is stored in a flash memory that can be mounted on the target device, and the program is executed by the program execution device.
  • the specific function may be realized by executing on a microcomputer (for example, a microcomputer that can be mounted on the target device).
  • the program can be stored and fixed in an arbitrary recording medium (not shown).
  • a recording medium (not shown) for storing and fixing the program may be mounted or connected to a device (such as a server device) different from the target device.
  • an open-circuit voltage estimating unit is formed by the parameter deriving unit 54 and the open-circuit voltage calculating unit 55, or by the parameter deriving unit 54, the parameter correcting unit 57, and the open-circuit voltage calculating unit 55 (FIG. 5, FIG. 9 and FIG. 11), the open-circuit voltage estimation unit may be considered to further include a parameter update unit 60 (see FIG. 13).
  • a device that excludes the charging rate estimation unit 56 from the estimation device 50, 50A, or 50B functions as an open-circuit voltage estimation device. It is also possible to delete the charging rate estimation unit 56 from the estimation device 50, 50A or 50B.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • General Chemical & Material Sciences (AREA)
  • Secondary Cells (AREA)

Description

電池状態推定装置及び蓄電池システム
　本発明は、電池状態推定装置及び蓄電池システムに関する。
　充放電停止時に蓄電池の端子電圧を測定することで蓄電池の開放電圧を検出することができる。但し、蓄電池の内部等価回路におけるインピーダンス回路は容量成分も含むため、充電又は放電が停止した後、暫くの間（例えば数分～数時間）は、端子電圧が過渡的に変動する。端子電圧の収束を待たずに安定開放電圧を推測するべく、或る従来方法では、充放電停止後、所定時間内の電圧値を用いて、複数の指数関数項を含む近似式の係数を最小二乗法を用いて算出する（下記特許文献１参照）。
国際公開２００６／１２９８０２号
　上記従来方法によって安定開放電圧（開放電圧の収束値）を推定することができるかもしれないが、複数の指数関数を用いたフィティングが必要になるため、パラメータの設定が煩雑となり、演算負荷が重い。
　そこで本発明は、簡素な演算により安定開放電圧（開放電圧の収束値）を正確に推定可能な電池状態推定装置及び蓄電池システムを提供することを目的とする。
　本発明に係る電池状態推定装置は、蓄電池の端子電圧である電池電圧の測定値を順次取得する測定値取得部と、前記蓄電池の充電又は放電の停止後の非通電期間中において、基準タイミングを設定するとともに、前記基準タイミングから見て前記電池電圧が所定の第１電圧量だけ変動する第１タイミング及び前記基準タイミングから見て前記電池電圧が前記第１電圧量より大きな所定の第２電圧量だけ変動する第２タイミングを前記測定値を用いて設定するタイミング設定部と、前記基準タイミング及び前記第１タイミング間の第１時間差と、前記第１タイミング及び前記第２タイミング間の第２時間差とに基づき、前記第２タイミングよりも後の対象タイミングにおける前記蓄電池の開放電圧を推定する開放電圧推定部と、を備える。
　本発明によれば、簡素な演算により安定開放電圧（開放電圧の収束値）を正確に推定可能な電池状態推定装置及び蓄電池システムを提供することが可能である。
本発明の実施形態に係る蓄電池システムの概略全体構成図である 本発明の実施形態に係る電池モジュールの内部等価回路図である。 本発明の実施形態に係り、電池モジュールの端子電圧及び内部電圧の時間推移を示すグラフである。 非通電期間中の電池電圧の時間変動を示す図である。 本発明の第１実施例に係る開放電圧／充電率推定装置の内部ブロック図である。 本発明の第１実施例の推定動作に関わる動作フローチャートである。 非通電期間中における測定電圧及び推定開放電圧の実験波形例を表す図である。 非通電期間中における測定電圧及び推定開放電圧の実験波形例を表す図である。 本発明の第２実施例に係る開放電圧／充電率推定装置の内部ブロック図である。 本発明の第２実施例の推定動作に関わる動作フローチャートである。 本発明の第３実施例に係る開放電圧／充電率推定装置の内部ブロック図である。 本発明の第３実施例に係る保持対象期間を示す図である。 本発明の第４実施例に係るパラメータ更新部を示す図である。 本発明の第４実施例の動作に関わるデータの流れを示す図である。 本発明の第５実施例で想定される複数のタイミングの関係を示す図である。 本発明の第７実施例に係る電圧変動パラメータの時間依存性を示す図である。
　以下、本発明の実施形態の例を、図面を参照して具体的に説明する。参照される各図において、同一の部分には同一の符号を付し、同一の部分に関する重複する説明を原則として省略する。尚、本明細書では、記述の簡略化上、情報、信号、物理量、状態量又は部材等を参照する記号又は符号を記すことによって該記号又は符号に対応する情報、信号、物理量、状態量又は部材等の名称を省略又は略記することがある。
　図１は、本発明の実施形態に係る蓄電池システム１の概略全体構成図である。蓄電池システム１は、電池モジュール１１、電流センサ１２、電圧センサ１３、温度センサ１４及び電池制御部１５を備える。電力ブロック２０も蓄電池システム１の構成要素に含まれている、と考えても良い。
　電池モジュール１１は、１以上の蓄電池（換言すれば二次電池）から成る。電池モジュール１１を形成する蓄電池は、任意の種類の蓄電池であり、例えば、リチウムイオン電池、ニッケル水素電池である。図１では、直列接続された複数の蓄電池にて電池モジュール１１が形成されているが、電池モジュール１１を形成する蓄電池の個数は１でも良い。電池モジュール１１に含まれる蓄電池の一部又は全部は、互いに並列接続されていても良い。以下では、電池モジュール１１を１つの蓄電池として捉えて考える。本実施形態において、放電及び充電とは、特に記述なき限り電池モジュール１１の放電及び充電を意味する。図２は、電池モジュール１１内の等価回路である。電池モジュール１１内の等価回路は、電圧ＶSを出力し且つ内部抵抗がゼロの電圧源ＶＳと、抵抗成分と容量成分を含むインピーダンス回路Ｚとの直列接続回路である、と考えることができる。
　電池モジュール１１には電力ブロック２０が接続されている。電力ブロック２０は、負荷及び電力源から成る。電池モジュール１１は、電力ブロック２０内の負荷に対して放電電力を供給することができると共に電力ブロック２０内の電力源から充電電力の供給を受けることができる。電池モジュール１１と、電力ブロック２０内の負荷及び電力源との間に、電力変換回路（不図示）が介在していても良い。
　電流センサ１２は、電池モジュール１１と電力ブロック２０との間に介在し、電池モジュール１１に流れる電流である電池電流Ｉを測定する。電圧センサ１３は、電池モジュール１１の端子電圧（即ち、電池モジュール１１の正極及び負極間の電位差）である電池電圧Ｖを測定する。温度センサ１４は、電池モジュール１１の温度である電池温度ＴＭＰを測定する。電池温度ＴＭＰは、例えば、電池モジュール１１内の蓄電池を包むパックの表面温度、又は、電池モジュール１１内の特定部位における温度である。尚、記号Ｉ、Ｖ、ＴＭＰは、夫々、電池電流の値、電池電圧の値、電池温度の値を示す記号としても用いられる。また、センサ１２、１３及び１４によって測定された電池電流値Ｉ、電池電圧値Ｖ及び電池温度ＴＭＰを、測定電流値Ｉ、測定電圧値Ｖ及び測定温度ＴＭＰと呼ぶこともあるし、記述の簡略化のために測定値Ｉ、Ｖ及びＴと呼ぶこともある。
　電池制御部１５は、測定値Ｉ、Ｖ及びＴＭＰを含む電池状態データを用いて電池モジュール１１の充電及び放電を制御する。また、電池制御部１５は、電池状態データを用いて電池モジュール１１の開放電圧等を推定する電池状態推定装置３０を有している。
　電池状態推定装置３０の説明に先立ち、電池状態推定装置３０の動作説明の具体化のために、図３に示すような状況を想定する。タイミングｔAからタイミングｔBまでの期間Ｐ０において電池モジュール１１の充電及び放電が行われておらず、その後、タイミングｔBからタイミングｔCまでの期間Ｐ１において電池モジュール１１の放電が行われ、その後、タイミングｔCからタイミングｔDまでの期間Ｐ２において電池モジュール１１の充電及び放電が停止され、タイミングｔDより後の期間Ｐ３において再び電池モジュール１１の放電が行われたことを想定する。
　期間Ｐ１及びＰ３の夫々は、電池モジュール１１の放電又は充電が行われている期間である通電期間である。期間Ｐ０及びＰ２の夫々は、電池モジュール１１の放電及び充電が行われていない期間である非通電期間である。以下では、電池モジュール１１の充電又は放電を電池モジュール１１の通電とも表現する。本実施形態では、通電期間Ｐ１が放電の成される放電期間であるが、通電期間Ｐ１は充電の成される充電期間であっても良く、この場合、タイミングｔCは充電の停止タイミングである。
　図３において、実線波形３０１は電池電圧Ｖの時間推移を表しており、破線波形３０２は電圧ＶS（図２参照）の時間推移を表している。放電期間Ｐ１において、電池電圧Ｖは、電圧ＶSから、インピーダンス回路Ｚの抵抗成分による電圧降下を差し引いたものに相当する。放電期間Ｐ１から非通電期間Ｐ２に切り替わると、電池電圧Ｖは、インピーダンス回路Ｚの存在により波形３０１に示す如く相応の時間（例えば数分～数時間）をかけて徐々に上昇し、期間Ｐ２中の電圧ＶSである一定の安定開放電圧ＶS＿Ｐ２に収束する。
　図４は、非通電期間Ｐ２中の波形３０１の拡大図である。タイミングｔ0～ｔ4を含むタイミングｔ0～ｔnは、タイミングｔC及びｔD間の複数の時刻を示す（ｎは整数）。任意の整数ｉについて、タイミングｔi+1はタイミングｔiより遅く、タイミングｔiにおける電池電圧Ｖを記号Ｖiにて表す。また、任意の整数ｉについて等式"Ｖi+1－Ｖi＝Ｎ"が成立する。Ｎは、正の所定値を持つ単位量（単位電圧量）であり、例えば数１０ミリボルトである。タイミングｔ0は、基準タイミングとして機能する。
　式（Ａ１）及び（Ａ２）に示すように、時間長さＴ及び電圧変動パラメータ（以下、単にパラメータとも言う）Ｒを定義する。ここで、"ｔi+1－ｔi"は、タイミングｔi及びｔi+1間の時間長さを表す。
　Ｔ＝ｔ1－ｔ0　　　　　　　　　　　　　…（Ａ１）
　Ｒ＝（ｔ2－ｔ1）／（ｔ1－ｔ0）　　　　…（Ａ２）
　式（Ａ２）に示すように、電圧変動パラメータＲは、第１時間長さ"ｔ1－ｔ0"に対する第２時間長さ"ｔ2－ｔ1"の比にて表される。第１時間長さ"ｔ1－ｔ0"は、基準タイミングｔ0から見て電池電圧Ｖが単位量Ｎだけ変動するのに要する時間長さであり、第２時間長さ"ｔ2－ｔ1"は、基準タイミングｔ0から見て電池電圧Ｖが単位量Ｎだけ変動した時点を起点として（即ちタイミングｔ1を起点として）、電池電圧Ｖが更に単位量Ｎだけ変動するのに要する時間長さである。尚、本実施形態では、非通電期間Ｐ２の前の通電期間Ｐ１が放電期間であることを想定しているため、非通電期間Ｐ２中の電池電圧Ｖの変動は"増加"である。仮に通電期間Ｐ１が充電期間であったならば、非通電期間Ｐ２中の電池電圧Ｖの変動は"減少"である。
　式（Ａ１）及び（Ａ２）より等式"ｔ2－ｔ1＝Ｔ×Ｒ"が成立する（図４も参照）。非通電期間Ｐ２において、電池電圧Ｖは安定開放電圧ＶS＿Ｐ２（図３参照）に向けて変動してゆくが、インピーダンス回路Ｚが容量成分を含んでいることに関連して、その変動速度は時間経過と共に徐々に減少してゆく。タイミングｔiの定義及び式（Ａ２）から明らかなように、電圧変動パラメータＲは、非通電期間Ｐ２中における電池電圧Ｖの変動速度の変化率（減少率）に応じた値を持つ。
　電池状態推定装置３０では、"非通電期間Ｐ２において、電池電圧Ｖが単位量Ｎだけ変動するのに要する時間長さが、電池電圧Ｖが単位量Ｎだけ変動（ここでは増加）するごとにＲ倍増大する"という仮定（換言すれば命題）ＡＳＭをおく。仮定ＡＳＭは、"非通電期間Ｐ２において電池電圧Ｖの変動速度の変化率が不変である"という仮定とも言える。
仮定（命題）ＡＳＭが真であるならば、等式"ｔ3－ｔ2＝Ｔ×Ｒ2"、"ｔ4－ｔ3＝Ｔ×Ｒ3"、・・・及び"ｔn－ｔn―1＝Ｔ×Ｒn-1"が成立すると共に下記式（Ａ３）が成立する。式（Ａ３）から式（Ａ４）が導かれる。
Figure JPOXMLDOC01-appb-M000001
　ここで、"ｔ"は非通電期間Ｐ２中の任意のタイミングを表す。
　 "ｔ－ｔ0"は、基準タイミングｔ0からタイミングｔまでの経過時間を表す。
　 "ΔＶ"は、タイミングｔ0及びｔ間における電池電圧Ｖの変動量を表す。
　"ＶEST（ｔ）"は、電池状態推定装置３０にて推定可能な、タイミングｔにおける電池電圧Ｖを表す。式（Ａ３）から式（Ａ４）を導く際、 関係式"ΔＶ＝ＶEST（ｔ）－Ｖ0" が利用される。
　電池状態推定装置３０は、タイミングｔ２以降の電池電圧Ｖの挙動を式（Ａ４）にてモデル化し、タイミングｔより前の時点において、タイミングｔにおける電池電圧Ｖを、推定式（Ａ４）に基づき"ＶEST（ｔ）"として推定（予測）することができる。尚、以下では、この推定を行う現タイミング（例えばｔ2）から見た未来のタイミングｔを対象タイミングと呼び、"ＶEST（ｔ）"を推定開放電圧（予測開放電圧）とも呼ぶ。
　以下、上述の構成を基本とする、蓄電池システム１（特に電池状態推定装置３０）の具体的な実施例を幾つか説明する。矛盾無き限り、以下に示す複数の実施例の内、任意の２以上の実施例を組み合わせることも可能である。
＜＜第１実施例＞＞
　第１実施例を説明する。図５は、電池状態推定装置３０に内包された開放電圧／充電率推定装置５０の内部ブロック図である。第１実施例に係る推定装置５０は、符号５１～５６によって参照される各部位を備える。
　測定値取得部５１は、センサ１２、１３及び１４の出力信号を所定のサンプリング周期でサンプリングすることにより周期的に測定電流値Ｉ、測定電圧値Ｖ及び測定温度ＴＭＰを取得する。但し、測定温度ＴＭＰの取得は割愛されうる。タイマ５２は、現在時刻を計測及び取得する、又は、任意の時刻からの経過時間を計測する。推定装置５０を含む電池状態推定装置３０内の各部位は、取得部５１の取得結果及びタイマ５２の取得／計測結果を自由に利用可能である。
　タイミング設定部５３は、電池モジュール１１の通電の停止後の非通電期間中において、測定電圧値Ｖを用いて、基準タイミング、第１及び第２タイミングを設定する。非通電期間は、電池電圧Ｖが過渡的に変動する電圧変動期間を含み、電圧変動期間中に基準タイミング、第１及び第２タイミングが設定される。第１実施例でも、開放電圧の推定が行われるべき非通電期間として期間Ｐ２を想定し、タイミングｔ0が基準タイミングに設定されるものとする（後述の他の実施例においても同様）。
　タイミング設定部５３は、例えば、電池モジュール１１の通電が停止してからの経過時間に基づき、基準タイミングｔ0を設定する。より具体的には例えば、期間Ｐ１及びＰ２間の境界タイミングｔC（放電が停止したタイミング）から所定時間が経過した時点を基準タイミングｔ0に設定する。
　或いは例えば、タイミング設定部５３は、電池モジュール１１の通電の停止後における電池電圧Ｖの電圧変動量に基づき、基準タイミングｔ0を設定してもよい。より具体的には例えば、期間Ｐ１及びＰ２間の境界タイミングｔC（放電が停止したタイミング）の測定電圧値Ｖである電圧値Ｖ［ｔC］から見て、測定電圧値Ｖが所定電圧量だけ増加したタイミングを基準タイミングｔ0に設定する。または、タイミングｔC以降、単位時間当たりの電池電圧Ｖの変動量を監視対象量として監視し、監視対象量が減少してきて所定値以下となったタイミングを基準タイミングｔ0に設定しても良い。
　タイミング設定部５３は、更に、基準タイミングｔ0から見て電池電圧Ｖが所定の第１電圧量だけ変動する第１タイミング及び基準タイミングから見て電池電圧Ｖが第１電圧量より大きな所定の第２電圧量だけ変動する第２タイミングを、順次取得される測定値Ｖを用いて、設定する。第１実施例において、第１電圧量は上述の単位量Ｎと等しく、且つ、第２電圧量は単位量Ｎの２倍であるとする（図４参照）。そうすると、第１及び第２タイミングは、夫々、図４のタイミングｔ1及びｔ2である。
　パラメータ導出部５４は、タイミング設定部５３にて設定された基準及び第１タイミング間の時間差（ｔ1－ｔ0）並びに第１及び第２タイミング間の時間差（ｔ2－ｔ1）に基づき、上記式（Ａ２）に従って電圧変動パラメータＲを導出する。
　開放電圧演算部５５は、第２タイミングｔ2以後であって且つ対象タイミングｔよりも前のタイミングにおいて、導出部５４にて導出された電圧変動パラメータＲを用い、上記式（Ａ４）に従って対象タイミングｔにおける開放電圧ＶEST（ｔ）を推定する。タイミング設定部５３は、基準タイミングｔ0での測定電圧値Ｖを基準電圧値Ｖ0として保持しており、演算部５５は、その保持された基準電圧値Ｖ0を式（Ａ４）に代入する。式（Ａ４）中の"Ｔ"の値は、タイミング設定部５３から演算部５５に供給される。また、演算部５５は、式（Ａ４）に代入される、基準タイミングｔ0から予測タイミングｔまでの経過時間（ｔ－ｔ0）を、タイマ５２を利用して認識できる。
　対象タイミングｔは、タイミングｔ2以後の任意のタイミングである。タイミングｔC、ｔ0又はｔ2から所定時間（例えば１時間又は数時間）が経過したタイミングが対象タイミングｔであっても良い。電池モジュール１１の特性にもよるが、通電停止後、１時間程度経過すれば、電池電圧Ｖは概ね安定開放電圧ＶS＿Ｐ２に達するため、タイミングｔCから１時間後のタイミングを対象タイミングｔに設定すれば、実用的な推定開放電圧ＶEST（ｔ）が得られ、ＶEST（ｔ）は安定開放電圧ＶS＿Ｐ２の推定値として機能する。
　上述の説明から理解されるように（式（Ａ４）の導出元の式（Ａ３）からも理解されるように）、式（Ａ４）を用いて開放電圧ＶEST（ｔ）を推定するということは、上記の仮定ＡＳＭが正しいとみなして（命題ＡＳＭが真であるとみなして）、電池モジュール１１の開放電圧を推定していることに等しい。
　充電率推定部５６は、電池モジュール１１の開放電圧と電池モジュール１１の充電率との関係を示す所定の演算式又はテーブルデータを用いて、推定開放電圧ＶEST（ｔ）から電池モジュール１１の充電率を推定する。電池モジュール１１の充電率は、電池モジュール１１の満充電容量に対する電池モジュール１１の残容量の比であり、一般にＳＯＣ（State of Charge）と呼ばれる。
　図６に、開放電圧の推定処理に関わる動作フローチャートを示す。まず、ステップＳ１１において、測定取得部５１による測定値Ｉ、Ｖ及びＴＭＰの周期的な取得（例えば、１ミリ秒を周期とする順次取得）が開始される。但し、測定値ＴＭＰの取得は割愛されうる。その後、タイミング設定部５３は、測定値Ｉを監視することで、今まで行われていた電池モジュール１１の通電（充電又は放電）が停止したか否かを判断し、通電の停止が確認されると、上述した方法に従って基準タイミングｔ0を設定すると共に、基準タイミングｔ0における測定電圧値を基準電圧値Ｖ0として保持する（ステップＳ１２及びＳ１３）。
更にその後、タイミング設定部５３は、順次取得される測定値Ｖを基準電圧値Ｖ0と比較することで、測定値Ｖが基準電圧値Ｖ0から第１電圧量（例えばＮ）だけ変動したタイミング及び測定値Ｖが基準電圧値Ｖ0から第２電圧量（例えば２×Ｎ）だけ変動したタイミングを、夫々、第１タイミング（ｔ1）及び第２タイミング（ｔ2）として抽出、設定する（ステップＳ１４）。設定された基準、第１、第２タイミングから電圧変動パラメータＲが導出され（ステップＳ１５）、当該パラメータＲを用いて開放電圧ＶEST（ｔ）が推定されると共に充電率が推定される（ステップＳ１６）。
　本方法によれば、電池電圧Ｖの変動量を観測するだけで、開放電圧推定に必要なパラメータＲが求まるため、簡素な演算により、安定開放電圧（開放電圧の収束値）ＶS＿Ｐ２を正確に推定することが可能となる。
　図７において、実線波形ＶREAL1及び破線波形ＶEST1は、或る条件下での非通電期間中Ｐ２における測定電圧Ｖ及び推定開放電圧ＶEST（ｔ）の実験波形を表している。両波形は良く一致し、良好な推定が成されていることが分かる。
　＜＜第２実施例＞＞
　但し、図７の波形ＶREAL1及びＶEST1は或る温度条件下（例えば電池温度ＴＭＰが２５℃の条件下）で得られるものであり、電池温度ＴＭＰによっては、図８に示す如く、推定開放電圧ＶEST（ｔ）と実際の電池電圧Ｖとの誤差が大きくなりうる。図８において、実線波形ＶREAL2及び破線波形ＶEST2は、低温度条件下（例えば電池温度ＴＭＰが０℃の条件下）での非通電期間中Ｐ２における測定電圧Ｖ及び推定開放電圧ＶEST（ｔ）の実験波形を表している。但し、波形ＶEST2は、後述のパラメータ補正を行わなかった場合における推定開放電圧ＶEST（ｔ）の波形である。第２実施例では、開放電圧の推定誤差を低減する構成を説明する。また、充電又は放電停止後の電池電圧Ｖの挙動は、電池温度ＴＭＰ以外の要因にも影響されるため、該挙動に影響を与える様々な要因に配慮した方が良い。
　図９は、第２実施例に係る開放電圧／充電率推定装置５０Ａの内部ブロック図である。
図１の電池状態推定装置３０に推定装置５０Ａを設けておくことができる。図５の推定装置５０にパラメータ補正部５７を追加することで図９の推定装置５０Ａが形成される。パラメータ補正部５７を除き、推定装置５０Ａの構成及び動作は、第１実施例で述べた推定装置５０のそれらと同じである。
　パラメータ補正部５７は、パラメータ導出部５４によって導出された電圧変動パラメータＲを、補正用指標に基づき補正する。補正後の電圧変動パラメータＲを記号"ＲC"にて表す。パラメータＲがパラメータＲCへと補正された場合、開放電圧演算部５５は、パラメータＲCを式（Ａ４）の"Ｒ"として用いた上で、上記式（Ａ４）に従って対象タイミングｔにおける開放電圧ＶEST（ｔ）を推定する。
　補正用指標は、以下の第１～第４指標の内、何れか１以上の指標を含む。補正部５７は、第１～第４指標の内の、２つ、３つ又は４つの指標を用いて電圧変動パラメータＲの補正を行っても良い。
　第１指標は、電池温度ＴＭＰである。第１指標としての電池温度ＴＭＰは、電池モジュール１１の通電の停止時点を基準とするタイミングでの測定温度ＴＭＰである。即ち例えば、第１指標としての電池温度ＴＭＰは、タイミングｔC、タイミングｔCより所定時間前のタイミング、又は、タイミングｔCより所定時間後のタイミングにおける測定温度ＴＭＰである。
　第２指標は、通電期間Ｐ１中における電池モジュール１１の通電の電流レート（単位時間当たりの電流量）である。
　第３指標は、通電期間Ｐ１中における電池モジュール１１の充電率の変化量ΔＳＯＣ［Ｐ１］である。図９の推定装置５０Ａ又は図１の電池状態推定装置３０は、通電期間Ｐ１中における測定電流値Ｉを積算し、その積算値と電池モジュール１１の満充電容量とから、第３指標としての充電率の変化量ΔＳＯＣ［Ｐ１］を求めることができる。
　第４指標は、非通電期間Ｐ２における電池モジュール１１の充電率ＳＯＣ［Ｐ２］である。非通電期間Ｐ２における電池モジュール１１の充電率ＳＯＣ［Ｐ２］は、通電期間Ｐ１の終了時点（即ちタイミングｔC）における電池モジュール１１の充電率と等しい。図９の推定装置５０Ａ又は図１の電池状態推定装置３０は、非通電期間Ｐ０における電池モジュール１１の充電率ＳＯＣ［Ｐ０］を、第１実施例で述べた方法を用いて求めることができる（即ち、非通電期間Ｐ０中に推定した開放電圧から求めることができる）、或いは、非通電期間Ｐ０中に実測した開放電圧（測定電圧Ｖ）から求めることができ、充電率ＳＯＣ［Ｐ０］と上記変化量ΔＳＯＣ［Ｐ１］を用いて充電率ＳＯＣ［Ｐ２］を求めることができる。
　非通電期間中における電池電圧Ｖの挙動を第１～第４指標を様々に変化させながら実際に調査しておき、その調査結果を用い、補正用指標に応じてパラメータＲをパラメータＲCに変換する演算式又はテーブルデータを作成し、得られた演算式又はテーブルデータをパラメータ補正部５７に保持させておくと良い。
　図１０は、第２実施例に係る、開放電圧の推定処理に関わる動作フローチャートである。パラメータＲの導出ステップＳ１５と開放電圧の推定ステップＳ１６との間に、パラメータの補正ステップＳ２０が挿入されている点を除き、図１０の動作フローチャートは図６のそれと同様である。導出されたパラメータＲがパラメータＲCへと補正され（ステップＳ２０）、パラメータＲCを用いて開放電圧ＶEST（ｔ）が推定されると共に充電率が推定される（ステップＳ１６）。
　第２実施例では、充電又は放電停止後の電池電圧Ｖの挙動に影響を与える補正用指標を考慮して電圧変動パラメータを補正するので、開放電圧の推定精度の向上が図られる。
　＜＜第３実施例＞＞
　第３実施例を説明する。図１１は、第３実施例に係る開放電圧／充電率推定装置５０Ｂの内部ブロック図である。図１の電池状態推定装置３０に推定装置５０Ｂを設けておくことができる。図５の推定装置５０に測定値保持部５８を追加することで図１１の推定装置５０Ｂが形成される。測定値保持部５８の存在及びその存在に対応する後述の内容を除き、推定装置５０Ｂの構成及び動作は、第１実施例で述べた推定装置５０のそれらと同じである。
　測定値保持部５８は、通電期間Ｐ１の終了後、測定値取得部５１にて順次取得される電池電圧Ｖの測定値を、所定時間ＴQ分だけ保持する。即ち、測定値保持部５８は、図１２の保持対象期間３２０中に取得された測定電圧値Ｖを、各測定時刻と対応付けて保持する（即ち、保持対象期間３２０中の測定電圧値Ｖの時系列データを保持する）。保持対象期間３２０は、タイミングｔCと、タイミングｔCから所定時間ＴQだけ経過したタイミングｔCQまでの期間である。測定値保持部５８は、所定時間ＴQを超える時間分の測定電圧値Ｖを保持していても構わない。
　第３実施例に係るタイミング設定部５３は、基準タイミングｔ0を設定してから基準タイミングｔ0を基準として第１及び第２タイミングｔ1及びｔ2を設定するのではなく、まず、保持対象期間３２０の終了タイミングｔCQを第２タイミングｔ2として取り扱い（図１２参照）、その後、保持部５８に保持されている各測定電圧値Ｖに基づき、第２タイミングｔ2を基準として基準タイミングｔ0及び第１タイミングｔ1を設定する。
　より具体的には例えば、第３実施例に係るタイミング設定部５３は、タイミングｔCQを第２タイミングｔ2として定めると共にタイミングｔCQとしての第２タイミングｔ2での測定電圧値Ｖを電圧値Ｖ2（図４参照）として取り扱う。その後例えば、設定部５３は、電圧値Ｖ2より単位量Ｎだけ低い測定電圧値Ｖ（即ちＶ1）及び電圧値Ｖ2より単位量Ｎの２倍だけ低い測定電圧値Ｖ（即ちＶ0）を保持部５８の保持内容から抽出すると共に、抽出した前者の測定電圧値Ｖ（即ちＶ1）に対応するタイミング及び抽出した後者の測定電圧値Ｖ（即ちＶ0）に対応するタイミングを、夫々、第１タイミングｔ1及び基準タイミングｔ0と定める。
　タイミングｔ0、ｔ1、ｔ2の設定後の動作は、第１実施例で述べた通りである。図１１の推定装置５０Ｂに図９のパラメータ補正部５７を追加し、第２実施例で述べたパラメータＲの補正及び補正後のパラメータＲ（即ちＲC）に基づく開放電圧ＶEST（ｔ）の推定を行うようにしても良い。
　上述の第１実施例では、基準電圧値Ｖ0の設定後、電圧値Ｖ2（＝Ｖ0＋２Ｎ）に相当する測定電圧値Ｖが実際に観測されるまで、タイミングｔ2が確定せず、開放電圧の推定を行うことができない。これに対し、第３実施例の方法によれば、通電停止後、所定時間ＴQを経過した時点がタイミングｔ2として設定されて当該時点において開放電圧の推定が確保される。システムによっては、通電が停止してから或る一定時間内に開放電圧の推定を確保したい（一定時間内に必ず開放電圧を知りたい）という要望もある。第３実施例は、このような要望に応える。
　＜＜第４実施例＞＞
　第４実施例を説明する。第４実施例では、第１～第３実施例と組み合わせ可能なパラメータ更新方法を説明する。推定装置５０、５０Ａ又は５０Ｂ（図５、図９又は図１１）に第４実施例に係るパラメータ更新部６０（図１３参照）を搭載可能であり、以下、第４実施例では、更新部６０が搭載された推定装置５０、５０Ａ又は５０Ｂの動作を説明する。
　第４実施例に係るタイミング設定部５３は、タイミングｔ0、ｔ1及びｔ2の設定に加えて、非通電期間Ｐ２中において、基準タイミングｔ0'を設定するとともに、基準タイミングｔ0'から見て電池電圧Ｖが上記第１電圧量（例えばＮ）だけ変動する第１タイミングｔ1'及び基準タイミングｔ0'から見て電池電圧Ｖが上記第２電圧量（例えば２Ｎ）だけ変動する第２タイミングｔ2'を、順次取得される測定値Ｖを用いて、設定する。基準タイミングｔ0'は、基準タイミングｔ0より後の任意のタイミングであり、タイミングｔ1又はｔ2と一致していても良いし、タイミングｔ2より後の所定タイミングであっても良い。
　タイミングｔ0'、ｔ1'及びｔ2'は更新用タイミングとして機能し、更新用タイミングｔ0'、ｔ1'及びｔ2'が設定されると、パラメータ導出部５４は、下記式（Ｂ１）に従って更新用パラメータＲ'を求める（図１４参照）。"ｔi+1'－ｔi'"は、タイミングｔi'及びｔi+1'間の時間長さを表す。
　Ｒ'＝（ｔ2'－ｔ1'）／（ｔ1'－ｔ0'）　　　　…（Ｂ１）
　パラメータ補正部５７が存在する場合、パラメータＲをパラメータＲCへ補正する方法と同様に、補正用指標に基づき更新用パラメータＲ'を更新用パラメータＲC'へと補正しても良い（図１４参照）。
　更新用パラメータＲ'又はＲC'が求められると、パラメータ更新部６０は、下記式（Ｂ２）又は式（Ｂ３）に従ってパラメータＲUPDATEを求める。式（Ｂ２）及び（Ｂ３）における重み係数ｗは、１であっても良いし、"０＜ｗ＜１"を満たす範囲内の任意の所定値を有しても良い。
　ＲUPDATE＝（１－ｗ）・Ｒ＋ｗ・Ｒ'                …（Ｂ２）
　ＲUPDATE＝（１－ｗ）・ＲC＋ｗ・ＲC'              …（Ｂ３）
　開放電圧演算部５５は、パラメータＲUPDATEの導出前は、パラメータＲ又はＲCに基づいて開放電圧ＶEST（ｔ）を推定するが、パラメータＲUPDATEの導出後は、パラメータＲUPDATEに基づいて開放電圧ＶEST（ｔ）を推定する。即ち、パラメータＲUPDATEの導出後、演算部５５は、パラメータＲUPDATEを式（Ａ４）の"Ｒ"として用いた上で、上記式（Ａ４）に従って対象タイミングｔにおける開放電圧ＶEST（ｔ）を推定する。このように、パラメータＲUPDATEの導出によって、開放電圧推定に用いる電圧変動パラメータが更新される（パラメータＲUPDATEは更新後のパラメータに相当する）。
　通電停止後、電池電圧Ｖは電池モジュール１１内部の様々な要因の影響を受けながら徐々に安定開放電圧ＶS＿Ｐ２に収束してゆくが、電圧変動期間の初期における電池電圧Ｖの挙動よりも、電圧変動期間の中期又は末期における電池電圧Ｖの挙動の方が、安定開放電圧ＶS＿Ｐ２への収束の様子を良く表している。このため、第４実施例の如く電圧変動パラメータを更新することにより、開放電圧の推定精度を向上させることができる。電圧変動パラメータを１回だけ更新するときの動作を説明したが、上述と同様の方法に従って、タイミングｔ0'、ｔ1'及びｔ2'の組を複数組設定し、更新動作を複数回以上行っても良い。
　＜＜第５実施例＞＞
　第１～第４実施例と組み合せて実施される第５実施例を説明する。上述の如く設定されたタイミングｔ0'、ｔ1'及びｔ2'は、電圧変動パラメータの更新に利用されるのではなく、以下のようにパラメータの設定に利用されても良い。図１５は、第５実施例で想定されるタイミングｔ0、ｔ1及びｔ2とタイミングｔ0'、ｔ1'及びｔ2'の具体例である。
図１５の例では、タイミングｔ0'、ｔ1'、ｔ2'における電池電圧Ｖは、夫々、Ｖ0'、（ Ｖ0'＋Ｎ）、 （Ｖ0'＋２Ｎ）である。
　パラメータ導出部５４は、上記（Ａ２）及び（Ｂ１）に基づき電圧変動パラメータＲ及びＲ'を求めた後、式（Ｂ４）に従って電圧変動パラメータＲJを求める。式（Ｂ４）における重み係数ｗは、 "０＜ｗ＜１"を満たす範囲内の任意の所定値であるが、典型的には０．５でよい。パラメータ補正部５７が存在する場合、パラメータＲをパラメータＲCへ補正する方法と同様に、補正用指標に基づきパラメータＲJをパラメータＲJCへと補正しても良い。演算部５５は、パラメータＲJ又はＲJCを式（Ａ４）の"Ｒ"として用いた上で、上記式（Ａ４）に従って対象タイミングｔにおける開放電圧ＶEST（ｔ）を推定する。
　ＲJ＝（１－ｗ）・Ｒ＋ｗ・Ｒ'                …（Ｂ４）
　＜＜第６実施例＞＞
　第１～第４実施例と組み合せて実施される第６実施例を説明する。上述の第５実施例は、４点以上のタイミングで電圧変動パラメータを設定する一形態と言える。以下のように、４点以上のタイミングで電圧変動パラメータを設定しても良い。タイミング設定部５３は、上述の如くタイミングｔ0、ｔ1及びｔ2を設定すると共に、更に、基準タイミングｔ0から見て電池電圧Ｖが所定の第３電圧量だけ変動する第３タイミングを、順次取得される測定値Ｖを用いて、設定する。第３電圧量は、上述の第１及び第２電圧量よりも多い。ここでは、第１、第２電圧量は夫々単位量Ｎの１倍、２倍であって、第３電圧量は単位量の３倍であるとする。そうすると、第３タイミングは、図４のタイミングｔ3と一致する。
　第６実施例に係るパラメータ導出部５４は、ｔ0～ｔ3を実測値として用いて、下記式（Ｃ１）～（Ｃ３）を精度良く満たすパラメータＲの値を、最小二乗法などを用いて導出する。パラメータＲの導出後の動作は、上述の各実施例で述べた通りである。
　ｔ1－ｔ0＝Ｔ　　　　　　　　　　　　　　　…（Ｃ１）
　ｔ2－ｔ1＝Ｔ×Ｒ　　　　　　　　　　　　　…（Ｃ２）
　ｔ3－ｔ2＝Ｔ×Ｒ2　　　　　　　　　　　　　…（Ｃ３）
　ここでは４点のタイミングｔ0～ｔ3を用いてパラメータＲを決定しているが、５点以上のタイミング（例えばｔ0～ｔ4）を用いてパラメータＲを決定しても良い。パラメータＲの導出用のタイミング数を増大させることで、パラメータＲがより適正化され、開放電圧の推定精度の向上が図られる。
　＜＜第７実施例＞＞
　第７実施例では、第１～第６実施例の何れかに適用可能な幾つかの技術を説明する。
　上述の説明では、基準タイミングｔ0から第ｉタイミングｔiまでの電池電圧Ｖの変動量が単位量Ｎのｉ倍であること（換言すれば、タイミングｔi及びｔi+1間における電池電圧Ｖの変動量が常に単位量Ｎであること）を主として想定しているが、タイミングｔi及びｔi+1間における電池電圧Ｖの変動量は、整数ｉに依存する対数関数又は指数関数で表現されていても良い。例えば、タイミングｔ0及びｔ1間における電池電圧Ｖの変動量が単位量Ｎであるとき、タイミングｔ1及びｔ2間における電池電圧Ｖの変動量は"ｌｏｇＮ"や"Ｎm"であっても良い（ｍは所定値）。この場合においても、パラメータ導出部５４及び開放電圧演算部５５は、上記仮定ＡＳＭが正しい（命題ＡＳＭが真である）とみなして、パラメータＲの導出及び開放電圧ＶEST（ｔ）の導出を行えば良い（開放電圧ＶEST（ｔ）の推定式は上記式（Ａ４）から変更されうる）。
　電圧変動パラメータが時間に依存せず変化しない場合、上記式（Ａ４）に従う推定開放電圧ＶEST（ｔ）は時間の経過と共に変動し続ける。その場合であっても、対象タイミングｔを適切に設定すれば（例えば、第１実施例で述べたようにタイミングｔCから１時間後のタイミングを対象タイミングｔに設定すれば）、推定開放電圧ＶEST（ｔ）は安定開放電圧ＶS＿Ｐ２を精度良く表す。しかしながら、推定開放電圧ＶEST（ｔ）を安定開放電圧ＶS＿Ｐ２に収束させるべく、電圧変動パラメータ（Ｒ、ＲC、ＲUPDATE、ＲJ又はＲJC）を時間の関数（指数関数等）にしても良い。即ち例えば、電圧変動パラメータの導出後、図１６に示す如く、電圧変動パラメータの値を時間経過に伴って増大させても良い。
　開放電圧の推定式として例示した上記式（Ａ４）では、対数関数を用いて推定式を形成しているが、開放電圧の推定式の形態は任意であり、対数関数以外の関数（例えば指数関数）を用いて開放電圧の推定式を形成しても良い。
　＜＜変形等＞＞
　本発明の実施形態は、特許請求の範囲に示された技術的思想の範囲内において、適宜、種々の変更が可能である。以上の実施形態は、あくまでも、本発明の実施形態の例であって、本発明ないし各構成要件の用語の意義は、以上の実施形態に記載されたものに制限されるものではない。上述の説明文中に示した具体的な数値は、単なる例示であって、当然の如く、それらを様々な数値に変更することができる。上述の実施形態に適用可能な注釈事項として、以下に、注釈１～注釈３を記す。各注釈に記載した内容は、矛盾なき限り、任意に組み合わせることが可能である。
［注釈１］
　上述の実施形態では、放電停止後の非通電期間Ｐ２中における開放電圧の推定動作を説明しているが、充電停止後の非通電期間中においても同様の開放電圧推定が可能である。
［注釈２］
　電池制御部１５、電池状態推定装置３０又は推定装置５０、５０Ａ若しくは５０Ｂである対象装置を、ハードウェア、或いは、ハードウェアとソフトウェアの組み合わせによって構成することができる。対象装置にて実現される機能の全部又は一部である任意の特定の機能をプログラムとして記述して、該プログラムを対象装置に搭載可能なフラッシュメモリに保存しておき、該プログラムをプログラム実行装置（例えば、対象装置に搭載可能なマイクロコンピュータ）上で実行することによって、その特定の機能を実現するようにしてもよい。上記プログラムは任意の記録媒体（不図示）に記憶及び固定されうる。上記プログラムを記憶及び固定する記録媒体（不図示）は対象装置と異なる機器（サーバ機器等）に搭載又は接続されても良い。
［注釈３］
　パラメータ導出部５４及び開放電圧演算部５５によって、又は、パラメータ導出部５４、パラメータ補正部５７及び開放電圧演算部５５によって、開放電圧推定部が形成されていると考えても良く（図５、図９、図１１参照）、開放電圧推定部は、パラメータ更新部６０を更に含んでいると考えて良い（図１３参照）。推定装置５０、５０Ａ又は５０Ｂから充電率推定部５６を除外した装置は、開放電圧推定装置として機能する。推定装置５０、５０Ａ又は５０Ｂから充電率推定部５６を削除することも可能である。
　　１　蓄電池システム
　１１　電池モジュール
　３０　電池状態推定装置
　５０、５０Ａ、５０Ｂ　推定装置
　５１　測定値取得部
　５３　タイミング設定部
　５４　パラメータ導出部
　５５　開放電圧演算部
　５６　充電率推定部
　５７　パラメータ補正部
　５８　測定値保持部
　６０　パラメータ更新部

Claims (10)

  1. 　蓄電池の端子電圧である電池電圧の測定値を順次取得する測定値取得部と、
    　前記蓄電池の充電又は放電の停止後の非通電期間中において、基準タイミングを設定するとともに、前記基準タイミングから見て前記電池電圧が所定の第１電圧量だけ変動する第１タイミング及び前記基準タイミングから見て前記電池電圧が前記第１電圧量より大きな所定の第２電圧量だけ変動する第２タイミングを前記測定値を用いて設定するタイミング設定部と、
    　前記基準タイミング及び前記第１タイミング間の第１時間差と、前記第１タイミング及び前記第２タイミング間の第２時間差とに基づき、前記第２タイミングよりも後の対象タイミングにおける前記蓄電池の開放電圧を推定する開放電圧推定部と、を備えた電池状態推定装置。
  2. 　前記開放電圧推定部は、前記第１及び第２時間差に基づき、前記非通電期間中における前記電池電圧の変動速度の変化率に応じたパラメータを導出するパラメータ導出部を有し、前記パラメータを用いて前記対象タイミングにおける前記開放電圧を推定する請求項１に記載の電池状態推定装置。
  3. 　前記パラメータは、第１時間長さに対する第２時間長さの比に相当し、
    　前記第１時間長さは、前記基準タイミングから見て前記電池電圧が所定の単位量だけ変動するのに要する時間長さであり、
    　前記第２時間長さは、前記基準タイミングから見て前記電池電圧が前記単位量だけ変動した時点を起点として、前記電池電圧が更に前記単位量だけ変動するのに要する時間長さであり、
    　前記第１電圧量は前記単位量と一致し、前記第２電圧量は前記単位量に応じた電圧量である請求項２に記載の電池状態推定装置。
  4. 　Ｒは、前記比の値を持ち、
    　前記開放電圧推定部は、前記非通電期間において、前記電池電圧が前記単位量だけ変動するのに要する時間長さが、前記電池電圧が前記単位量だけ変動するごとにＲ倍増大するとみなして、前記対象タイミングにおける前記開放電圧を推定する請求項３に記載の電池状態推定装置。
  5. 　前記開放電圧推定部は、前記パラメータ導出部によって導出される前記パラメータを補正用指標に基づき補正するパラメータ補正部を更に有し、補正後のパラメータを用いて前記対象タイミングにおける前記開放電圧を推定し、
    　前記補正用指標は、
    　前記蓄電池の温度、
    　前記蓄電池の充電又は放電が行われていた通電期間中の前記蓄電池の充電又は放電の電流レート、
    　前記通電期間中における前記蓄電池の充電率の変化量、及び、
    　前記非通電期間における前記蓄電池の充電率の内、少なくとも１つを含む請求項２～請求項４の何れかに記載の電池状態推定装置。
  6. 　前記タイミング設定部は、前記蓄電池の充電又は放電が停止してからの経過時間、又は、前記蓄電池の充電又は放電の停止後における前記電池電圧の変動量に基づき、前記基準タイミングを設定する請求項１～請求項５の何れかに記載の電池状態推定装置。
  7. 　前記蓄電池の充電又は放電の停止後、所定時間分以上の前記測定値を保持する測定値保持部を更に備え、
    　前記タイミング設定部は、前記蓄電池の充電又は放電が停止してから前記所定時間だけ経過したタイミングを前記第２タイミングとして設定した後、前記第２タイミングを基準として前記測定値保持部の保持内容を用い前記基準タイミング及び第１タイミングを設定する請求項１～請求項６の何れかに記載の電池状態推定装置。
  8. 　前記タイミング設定部は、前記非通電期間中において、前記基準タイミングより後の更新用基準タイミングを設定するとともに、前記更新用基準タイミングから見て前記電池電圧が前記第１電圧量だけ変動する更新用第１タイミング及び前記更新用基準タイミングから見て前記電池電圧が前記第２電圧量だけ変動する更新用第２タイミングを前記測定値を用いて設定し、
    　前記開放電圧推定部は、前記更新用基準タイミング及び前記更新用第１タイミング間の時間差と、前記更新用第１タイミング及び前記更新用第２タイミング間の時間差とに基づき、前記第１及び第２時間差を用いて導出した前記パラメータを更新するパラメータ更新部を更に有し、更新後のパラメータを用いて前記対象タイミングにおける前記開放電圧を推定する請求項２～請求項４の何れかに記載の電池状態推定装置。
  9. 　前記開放電圧推定部による推定開放電圧に基づき前記蓄電池の充電率を推定する充電率推定部を更に備えた請求項１～請求項８の何れかに記載の電池状態推定装置。
  10. 　蓄電池と、
    　請求項１～請求項９の何れかに記載の電池状態推定装置と、を備えた蓄電池システム。
PCT/JP2013/005628 2012-09-26 2013-09-24 電池状態推定装置及び蓄電池システム Ceased WO2014050073A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2014538171A JP6158195B2 (ja) 2012-09-26 2013-09-24 電池状態推定装置及び蓄電池システム
US14/428,665 US9846200B2 (en) 2012-09-26 2013-09-24 Battery state estimation device and storage battery system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012211933 2012-09-26
JP2012-211933 2012-09-26

Publications (1)

Publication Number Publication Date
WO2014050073A1 true WO2014050073A1 (ja) 2014-04-03

Family

ID=50387511

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/005628 Ceased WO2014050073A1 (ja) 2012-09-26 2013-09-24 電池状態推定装置及び蓄電池システム

Country Status (3)

Country Link
US (1) US9846200B2 (ja)
JP (1) JP6158195B2 (ja)
WO (1) WO2014050073A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017511566A (ja) * 2014-05-12 2017-04-20 ローベルト ボッシュ ゲゼルシャフト ミット ベシュレンクテル ハフツング バッテリの温度を定める方法
JPWO2017221899A1 (ja) * 2016-06-22 2018-11-29 株式会社豊田自動織機 蓄電装置及び蓄電装置の制御方法
WO2025084110A1 (ja) * 2023-10-16 2025-04-24 パナソニックＩｐマネジメント株式会社 電池状態推定システム、電池状態推定方法、電池状態推定プログラム、及び記録媒体
WO2025197446A1 (ja) 2024-03-22 2025-09-25 パナソニックＩｐマネジメント株式会社 電池状態推定システム、電池状態推定方法、及び電池状態推定プログラム

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015166710A (ja) * 2014-03-04 2015-09-24 ソニー株式会社 蓄電部材状態推定装置、電池パック、電動車両、蓄電装置および蓄電部材状態推定方法
US11144106B2 (en) 2015-04-13 2021-10-12 Semiconductor Components Industries, Llc Battery management system for gauging with low power
US10224579B2 (en) 2015-12-31 2019-03-05 Robert Bosch Gmbh Evaluating capacity fade in dual insertion batteries using potential and temperature measurements
US10263447B2 (en) 2016-01-29 2019-04-16 Robert Bosch Gmbh Secondary battery management system
US10686321B2 (en) 2016-01-29 2020-06-16 Robert Bosch Gmbh Secondary battery management
US10243385B2 (en) 2016-01-29 2019-03-26 Robert Bosch Gmbh Secondary battery management system
US9960625B2 (en) 2016-03-31 2018-05-01 Robert Bosch Gmbh Battery management system with multiple observers
US10447046B2 (en) 2016-09-22 2019-10-15 Robert Bosch Gmbh Secondary battery management system with remote parameter estimation
KR102239365B1 (ko) 2017-10-20 2021-04-09 주식회사 엘지화학 배터리 충전 상태 추정 장치
US20200090425A1 (en) * 2018-09-18 2020-03-19 Cambridge Mobile Telematics Inc. Using vehicle electrical system monitored values
KR102856950B1 (ko) * 2019-12-11 2025-09-05 주식회사 엘지에너지솔루션 배터리 퇴화도 진단 장치 및 방법

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004109007A (ja) * 2002-09-19 2004-04-08 Yazaki Corp バッテリの開回路電圧推定方法及び装置
JP2008145199A (ja) * 2006-12-07 2008-06-26 Auto Network Gijutsu Kenkyusho:Kk 充電深度判定方法及び充電深度判定装置
JP2008256436A (ja) * 2007-04-03 2008-10-23 Yazaki Corp 近似式算出装置及びその方法、並びに、バッテリ状態監視装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3869676B2 (ja) * 2000-12-08 2007-01-17 矢崎総業株式会社 車両用バッテリの開回路電圧推定方法及び装置
JP2003307556A (ja) * 2002-02-15 2003-10-31 Yazaki Corp バッテリの開回路電圧推定方法及び装置
JP5170876B2 (ja) 2005-06-03 2013-03-27 古河電気工業株式会社 充電率／残存容量推定方法、電池の状態検知センサ及び電池電源システム
JP5882336B2 (ja) * 2011-08-30 2016-03-09 三洋電機株式会社 バッテリシステム、電動車両、移動体、電力貯蔵装置および電源装置
JP2013083612A (ja) * 2011-10-12 2013-05-09 Mitsumi Electric Co Ltd 電池状態計測方法及び電池状態計測装置
JP5867195B2 (ja) * 2012-03-13 2016-02-24 日産自動車株式会社 二次電池の充電装置および二次電池の充電方法
WO2013141100A1 (ja) * 2012-03-21 2013-09-26 三洋電機株式会社 電池状態推定装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004109007A (ja) * 2002-09-19 2004-04-08 Yazaki Corp バッテリの開回路電圧推定方法及び装置
JP2008145199A (ja) * 2006-12-07 2008-06-26 Auto Network Gijutsu Kenkyusho:Kk 充電深度判定方法及び充電深度判定装置
JP2008256436A (ja) * 2007-04-03 2008-10-23 Yazaki Corp 近似式算出装置及びその方法、並びに、バッテリ状態監視装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017511566A (ja) * 2014-05-12 2017-04-20 ローベルト ボッシュ ゲゼルシャフト ミット ベシュレンクテル ハフツング バッテリの温度を定める方法
JPWO2017221899A1 (ja) * 2016-06-22 2018-11-29 株式会社豊田自動織機 蓄電装置及び蓄電装置の制御方法
WO2025084110A1 (ja) * 2023-10-16 2025-04-24 パナソニックＩｐマネジメント株式会社 電池状態推定システム、電池状態推定方法、電池状態推定プログラム、及び記録媒体
WO2025197446A1 (ja) 2024-03-22 2025-09-25 パナソニックＩｐマネジメント株式会社 電池状態推定システム、電池状態推定方法、及び電池状態推定プログラム

Also Published As

Publication number Publication date
JP6158195B2 (ja) 2017-07-05
US9846200B2 (en) 2017-12-19
US20150234013A1 (en) 2015-08-20
JPWO2014050073A1 (ja) 2016-08-22

Similar Documents

Publication Publication Date Title
JP6158195B2 (ja) 電池状態推定装置及び蓄電池システム
JP5897701B2 (ja) 電池状態推定装置
TWI384246B (zh) 藉電池電壓變化模式估測開路電壓以估測電池電阻特徵之裝置及方法
CN107991623B (zh) 一种考虑温度和老化程度的电池安时积分soc估计方法
EP3002597B1 (en) Battery control device
He et al. Online model-based estimation of state-of-charge and open-circuit voltage of lithium-ion batteries in electric vehicles
TWI708068B (zh) 用於測定鋰硫電池之健康狀態及充電狀態的方法及裝置
JP6441913B2 (ja) バッテリに蓄積された電荷のモニタリング
KR102347014B1 (ko) 축전지 잔량 추정 장치, 축전지 잔량 추정 방법, 및 프로그램
JP5269994B2 (ja) バッテリーの電圧挙動を用いたバッテリーのｓｏｈ推定装置及び方法
JP4473823B2 (ja) 電気エネルギー蓄積器のための複数の部分モデルを用いた状態量およびパラメータの推定装置
JP2014025739A (ja) 電池状態推定装置
AU2022266288B2 (en) Battery management device, and electric power system
JP2014025738A (ja) 残容量推定装置
CN106104907B (zh) 二次电池内部温度推定装置以及二次电池内部温度推定方法
JP6044512B2 (ja) 電池特性学習装置
JP2006275797A (ja) 蓄電デバイスの残存容量演算装置
CN105353313A (zh) 电池荷电状态的估算方法和装置
JP2014068468A (ja) 充電制御装置
JP2011169831A (ja) 電池状態検知装置及び電池状態検知方法
JP2013250071A (ja) 満充電容量検出装置及び蓄電池システム
CN110596604A (zh) 一种基于安时积分法的锂电池soc估计方法
JP2013246088A (ja) 電池の内部抵抗推定方法及びその装置
CN109669138B (zh) 一种精准测定动力铅蓄电池组剩余容量的方法
JP2019060807A (ja) 二次電池の状態診断装置および状態診断方法

Legal Events

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

Ref document number: 13840374

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014538171

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14428665

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13840374

Country of ref document: EP

Kind code of ref document: A1