WO2019017411A1 - 推定装置、蓄電装置、推定方法、及びコンピュータプログラム - Google Patents

推定装置、蓄電装置、推定方法、及びコンピュータプログラム Download PDF

Info

Publication number
WO2019017411A1
WO2019017411A1 PCT/JP2018/026994 JP2018026994W WO2019017411A1 WO 2019017411 A1 WO2019017411 A1 WO 2019017411A1 JP 2018026994 W JP2018026994 W JP 2018026994W WO 2019017411 A1 WO2019017411 A1 WO 2019017411A1
Authority
WO
WIPO (PCT)
Prior art keywords
characteristic
charge
storage
discharge
storage amount
Prior art date
Application number
PCT/JP2018/026994
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
Priority claimed from JP2018133649A external-priority patent/JP6477959B2/ja
Application filed by 株式会社Ｇｓユアサ filed Critical 株式会社Ｇｓユアサ
Priority to EP18834355.2A priority Critical patent/EP3657189A4/en
Priority to US16/631,614 priority patent/US11237214B2/en
Publication of WO2019017411A1 publication Critical patent/WO2019017411A1/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/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]
    • 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
    • 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
    • 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
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
    • 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 an estimation device, a power storage device including the estimation device, an estimation method of an internal state, and a computer program.
  • LiCoO 2 LiCoO 2
  • LiMeO 2 LiMeO 2
  • an Li 2 MnO 3 -based active material which is a lithium excess type is being studied.
  • This material has the property of hysteresis which causes a difference in voltage and electrochemical characteristics between SOC-OCV (Open Circuit Voltage) at the time of charge and discharge for the same SOC (State Of Charge). .
  • SOC-OCV Open Circuit Voltage
  • SOC-OCV curve can not be uniquely determined, it is also difficult to predict the dischargeable energy at a certain point.
  • the material of the lithium excess type has a property of voltage drop (hereinafter referred to as VF) in which the SOC-OCP (Open Circuit Potential) curve of the positive electrode changes substantially over the entire area by repetition of charge and discharge. Since the value of the average discharge potential decreases, it is necessary to estimate not only the dischargeable capacity but also the dischargeable power amount as the current state of health (SOH). Since the SOC-OCV curve shape of the battery cell based on the unipolar SOC-OCP curve (hereinafter, also simply referred to as "cell”) changes significantly due to deterioration even if the latest charge / discharge history is the same, OCV The law can not be adopted.
  • SOC-OCP Open Circuit Potential
  • the conditions under which the latest charge and discharge history is the same include, for example, charge after passing through a complete discharge state.
  • the SOC-OCV curve shape of the cell is significantly changed because the unipolar SOC-OCP curve changes according to the deterioration.
  • the SOC-OCV of the storage element also changes depending on the state of the counter electrode, in order to control the storage element, it is necessary to estimate the internal state of the storage element including the state of the counter electrode. That is, the SOC-OCP characteristics of the positive electrode at current full charge and complete discharge, state of charge (proportion of SOC-OCP of positive electrode), degree of contraction of positive electrode, degree of contraction of negative electrode, capacity balance between positive electrode and negative electrode It is necessary to estimate parameters such as the amount of deviation and polarization. And it is necessary to estimate SOC correctly.
  • the current SOH and SOC estimation techniques for non-aqueous electrolyte secondary batteries make it difficult to use an active material having VF and hysteresis properties.
  • the measurement unit of the deterioration determination device disclosed in Patent Document 1 measures the capacity versus open-circuit voltage characteristic of the lithium ion secondary battery.
  • the determination unit can set a parameter for specifying the open-circuit voltage characteristic, and uses the parameter to specify the open-circuit voltage characteristic that substantially matches the open-circuit voltage characteristic measured by the measurement unit. judge.
  • the deterioration determination device of Patent Document 1 is for determining the deterioration state due to wear and lithium deposition, and can not estimate the internal state of the storage element including the active material that generates VF.
  • the present invention provides an estimation device capable of estimating the internal state of a storage element having a single pole whose storage capacity-potential characteristic changes due to repeated charge / discharge, a storage device provided with the estimation device, an estimation method of the internal state, and a computer program Intended to be provided.
  • the storage amount means the SOC, the amount of power that can be released, and the like.
  • An estimation device is a storage element having a single electrode including an active material in which a first characteristic that is a storage amount-potential charge characteristic and a second characteristic that is a storage amount-potential discharge characteristic change with repetition of charge and discharge. Estimate the internal state.
  • the estimation device changes according to repetition of charge and discharge, an acquisition unit that acquires a part of information of a third characteristic that is a storage amount-voltage charge characteristic of the storage element and a fourth characteristic that is a storage amount-voltage discharge characteristic.
  • a plurality of storage amount characteristics which are at least one of the first characteristic, the second characteristic, the third characteristic, the fourth characteristic, and V-dQ / dV, are stored corresponding to the change of the characteristic value, or the characteristic value
  • a first estimation unit configured to estimate an internal state of the storage element based on the information and the storage amount characteristic or the function.
  • the internal state of the storage element having a single pole in which the storage amount-potential characteristic changes due to repetition of charge and discharge can be favorably estimated.
  • FIG. 5 is a conceptual diagram showing a relationship between a potential range of a positive electrode corresponding to a predetermined voltage range and a range of charge amount corresponding to each deterioration state in each potential range.
  • FIG. 2 is a perspective view showing an example of a power storage device. It is a perspective view which shows the other example of an electrical storage apparatus. It is an exploded perspective view of a battery module. It is a block diagram of a battery module. It is a flowchart which shows the procedure of the estimation process of the parameter and storage amount characteristic by CPU.
  • FIG. 6 is a graph showing the relationship between the degree of deterioration k ′ and the amount of charge Q in the potential range P1 to P2V. It is a graph which shows discharge DOD-V of a positive electrode, discharge DOD-V of a battery, and DOD-V of a negative electrode which were calculated using the parameter calculated
  • the figure shows the result of finding the difference between the measured SOC-V and the second SOC-V calculated by calculating the parameters according to the above-mentioned method when charging and discharging 2.0 to 4.6 V per battery is repeated 500 cycles.
  • It is a flowchart which shows the procedure of the SOC estimation process by CPU. It is a flowchart which shows the procedure of the SOC estimation process by CPU.
  • An estimation device includes a storage element having a single electrode including an active material in which a first characteristic that is a storage amount-potential charge characteristic and a second characteristic that is a storage amount-potential discharge characteristic change with repetition of charge and discharge. Estimate the internal state.
  • the estimation device includes an acquisition unit for acquiring information of a part of the fourth characteristic, which is the third characteristic and / or the storage characteristic and voltage discharge characteristic of the storage element, and the charge and discharge are repeated.
  • a plurality of storage amount characteristics which are at least one of the first characteristic, the second characteristic, the third characteristic, the fourth characteristic, and V-dQ / dV, are stored corresponding to the change of the changing feature value, or
  • a first estimation unit configured to estimate an internal state of the storage element based on the information and the storage amount characteristic or the function.
  • V ⁇ dQ / dV is a relationship between the potential of the single pole or the voltage V of the storage element, and a differential value dQ / dV obtained by differentiating charge amount of charge or discharge capacity Q by V.
  • the plurality of stored charge amount characteristics reflect the deterioration based on the unipolar VF according to the feature value.
  • the information corresponds to a feature value.
  • the current internal state of the storage element including the state of the counter electrode, can be favorably estimated by a simple method.
  • the acquisition unit acquires information of a third characteristic and / or a fourth characteristic in a predetermined voltage range of the storage element, obtains a feature value from the information, and based on the feature value,
  • a second estimation unit configured to estimate a third characteristic or a fourth characteristic at the time of acquiring the information with reference to the storage amount characteristic or the function, the first estimation unit estimating the information and the information
  • the parameter indicating the internal state of the storage element may be estimated based on the third characteristic or the fourth characteristic.
  • the storage amount characteristic calculated based on the feature value and the stored storage amount characteristic, or based on the storage amount characteristic calculated based on the function, and the information measured in a predetermined voltage range, for example Parameters are obtained by regression calculation.
  • the current internal state of the storage element can be estimated with high accuracy, including the state of the opposite electrode, by a simple method. Even in a complicated use environment, the internal state of the storage element can be monitored with high accuracy.
  • the characteristic value may be at least one of a charge quantity, a discharge capacity, and an average discharge potential in the predetermined voltage range.
  • the voltage range of the storage element corresponds to a potential range in which there is a linear relationship between the amount of charge electricity or discharge capacity and the average discharge potential, and the potential difference with the counter electrode (voltage of the storage element) does not change before and after deterioration. It is set as a predetermined voltage range.
  • the charge quantity of electricity or the discharge capacity in a predetermined voltage range is used as the feature value.
  • a plurality of storage amount characteristics are stored in association with the feature value corresponding to the degree of deterioration, or stored as a function of the feature value.
  • the current storage amount characteristic can be accurately estimated using this. Even when the average discharge potential is used, the storage amount characteristic can be estimated with high accuracy.
  • the current third characteristic and the fourth characteristic can be accurately estimated based on the current storage amount-potential characteristic, the storage amount-voltage characteristic, or V-dQ / dV.
  • the parameters include a first characteristic, a second characteristic, a ratio of the first characteristic, a degree of contraction of the single electrode, a degree of contraction of the opposite electrode, and a displacement of a capacity balance between the single electrode and the opposite electrode It may be at least one selected from the amount and the polarization.
  • the second estimation unit estimates V-dQ / dV by interpolation calculation with reference to the V-dQ / dV, and the parameter is a processing coefficient required for the interpolation calculation. May be included.
  • the third characteristic or the fourth characteristic is obtained, so that the calculation accuracy and the calculation speed are good. Then, by using the processing coefficient necessary for the interpolation calculation as a parameter, the third characteristic or the fourth characteristic can be corrected with high accuracy by the parameter obtained by the regression calculation.
  • the first estimation unit is configured to minimize the difference between the information and the third characteristic or the fourth characteristic estimated by the second estimation unit by optimization calculation. Parameters may be determined.
  • the current parameter can be determined with high accuracy by optimization calculation, and the current internal state of the storage element can be detected with high accuracy.
  • the above-described estimation apparatus may further include a third estimation unit that estimates a third characteristic or a fourth characteristic using the parameter estimated by the first estimation unit.
  • the active material causes two or more electrochemical reactions according to the transition of charge and discharge, and the hysteresis between the third characteristic and the fourth characteristic shown when one electrochemical reaction occurs is
  • the storage amount is smaller than the hysteresis when another electrochemical reaction occurs, and the acquired voltage based on the third characteristic or the fourth characteristic estimated by the third estimation unit and the charge / discharge history of the storage element
  • a fourth estimation unit may be provided to estimate a third characteristic or a fourth characteristic for reference when estimating.
  • the effect of hysteresis is accurate based on the third characteristic or the fourth characteristic, parameter according to the current internal state of the storage element, and the charge / discharge history of the storage element.
  • the third or fourth characteristic for reference can be estimated.
  • the above-described estimation apparatus may further include a fifth estimation unit configured to estimate a storage amount based on the history of charge and discharge, the third characteristic or the fourth characteristic for reference, and the acquired voltage.
  • the storage amount of the storage element having the active material having the property of VF and the storage amount-voltage characteristic exhibiting hysteresis since the voltage is used, the storage amount is not limited to the SOC, and the amount of current energy stored in the storage element, such as the amount of power, can be estimated. Based on the charge / discharge characteristics, dischargeable energy up to SOC 0% and charge energy required up to SOC 100% can be predicted. The present remaining power and storable power can be estimated. Therefore, it is possible to accurately perform balancing in the case of using a plurality of storage elements, control of regeneration reception, estimation of a traveling distance in the case of mounting the storage elements, and the like.
  • a power storage device includes a power storage element and any one of the estimation devices described above.
  • the storage amount of the storage element can be accurately estimated even in a complicated use environment.
  • An estimation method includes a storage element having a single electrode including an active material in which a first characteristic, which is a storage amount-potential charge characteristic, and a second characteristic, which is a storage amount-potential discharge characteristic, change with repetition of charge and discharge. Estimate the internal state.
  • the estimation method includes a first characteristic, a second characteristic, a third characteristic which is a storage amount-voltage charge characteristic of the storage element, a storage amount-voltage discharge characteristic, corresponding to a change in the feature value which changes due to repetition of charge and discharge.
  • a plurality of fourth characteristics, and at least one of V-dQ / dV, is stored in the storage unit as a plurality or stored as a function of the feature value;
  • the partial information of three characteristics and / or the fourth characteristic is acquired, and the internal state of the storage element is estimated based on the information and the storage amount characteristic or the function.
  • the current internal state of the storage element can be estimated with high accuracy by detecting the states of the single pole and the opposite pole even in a complicated use environment by a simple method.
  • a computer program includes a storage element having a single electrode including an active material in which a first characteristic, which is a storage amount-potential charge characteristic, and a second characteristic, which is a storage amount-potential discharge characteristic, change with repetition of charge and discharge.
  • the computer for estimating the internal state acquires a third characteristic which is a storage amount-voltage charge characteristic of the storage element and / or a part of information of a fourth characteristic which is a storage amount-voltage discharge characteristic, a first characteristic, Referring to a table in which a plurality of storage amount characteristics, which are at least one of a second characteristic, a third characteristic, a fourth characteristic and V-dQ / dV, are stored corresponding to changes in feature values, or A process of estimating the internal state of the storage element based on the information and the storage amount characteristic or the function is performed with reference to the function of the feature value.
  • the single electrode of the electrode body of the storage element includes an active material having the property of VF and having the storage amount-potential characteristic hysteresis.
  • the active material has the property of VF
  • the shapes of the unipolar SOC-OCP curve and the SOC-OCV curve of the storage element change due to repetition of charge and discharge.
  • An electric storage element containing this active material flows a minute current, and is charged when fully charged to fully charged (full charge) and when fully charged or fully discharged (fully discharged). It has a hysteresis whose maximum potential difference between V curves is 100 mV or more.
  • FIG. 1 is a graph showing an example of the SOC-OCP of the positive electrode.
  • the horizontal axis is SOC (%), and the vertical axis is the potential E as OCP (VvsLi / Li + : Li / Li + ).
  • the charge / discharge curve before deterioration is shown by a broken line, and the charge / discharge curve after deterioration is shown by a solid line. As shown in FIG. 1, the deterioration causes VF, and the charge / discharge curve shifts downward.
  • the hysteresis is not provided, and the unipolar SOC-OCP curve does not change by repetition of charge and discharge.
  • the shape of the SOC-OCV curve of the storage element changes due to the repetition of charge and discharge due to the deterioration of the single pole (the reduction of the curve) or the increase of the amount of displacement of the capacity balance between both electrodes.
  • the internal states of the current storage element, single pole, and counter pole are estimated.
  • the charge amount-potential charge characteristic (first characteristic) of the single pole or its V-dQ / dV, the charge amount-potential discharge characteristic (second characteristic) of the single pole or its V-dQ / DV dV above, referred to as k
  • ratio of charge amount-potential charge characteristic (first characteristic) of unipolar Rcha described later
  • degree of contraction of positive electrode m
  • degree of contraction of negative electrode n
  • positive electrode At least one selected from the displacement amount ( ⁇ Qim) of the capacity balance between the negative electrodes, the polarization ( ⁇ V), and the processing coefficient (q) of interpolation calculation described later can be mentioned.
  • the k can also be expressed as the degree of deterioration k ′ (a state without deterioration is 1.0).
  • m and n respectively indicate the degree of contraction in the x-axis direction of the storage amount-V curve of the positive electrode and the storage amount-V curve of the negative electrode.
  • ⁇ V indicates the degree of contraction in the y-axis direction of the storage amount-V curve.
  • the information processing unit 60 stores a plurality of storage capacity characteristics according to the deterioration, or stores them as a function of the feature value.
  • the storage amount characteristic a single charge SOC-V (first characteristic) or a discharge SOC-V (second characteristic), a charge SOC-V (third characteristic) of the storage element or a discharge SOC-V (fourth characteristic) And unipolar charge V-dQ / dV or discharge V-dQ / dV, charge V-dQ / dV of storage element or discharge V-dQ / dV, and the like.
  • the characteristic value includes the amount of charge, the discharge capacity, the average discharge potential, the absolute value of the dQ / dV or dV / dQ peak, or the peak intensity ratio thereof.
  • the peak intensity ratio of dQ / dV the ratio of the peak integrated intensity of dQ / dV in the predetermined potential range a and the peak integrated intensity in the potential range b can be mentioned. By taking this ratio, it is possible to estimate the current charge amount characteristic regardless of the decrease in the degree m of contraction of the positive electrode. Among them, the amount of charge and the discharge capacity are preferable.
  • the charge / discharge curve shape changes continuously and uniquely according to the change (deterioration) of the feature value.
  • LiMeO 2 -Li respect 2 MnO 3 system of the active material along with the repetition of charge and discharge, the crystal structures have been reported to vary (Journal of Power Sources, vol.229 ( 2013), pp239-248). It is considered that the shape of the charge / discharge curve changes with the change of the crystal structure. It is suggested from the results of the paper that changes in crystal structure occur continuously in the short-term one temperature level charge / discharge cycle. And, from the report that the crystal structure has changed from layered to spinel-like crystals, it is presumed that the way of change is one way.
  • the crystal structure changes continuously and uniquely at one temperature level in the short term. From this report, the present inventors considered that the charge / discharge curve shape changes continuously and uniquely according to the change of the crystal structure in a long-term and in any use history. From the experimental results described later, it was confirmed that the charge / discharge curve shape changes continuously and uniquely in the long run even if the usage history is different.
  • the charge and discharge curve shape of a single pole does not change with repetition of charge and discharge.
  • the charge / discharge curve shape of the cell changes individually, that is, unambiguously, due to the repetition of charge / discharge.
  • the storage amount characteristic of the single pole changes continuously and unambiguously with the change of the feature value, so that the transition of the change of the storage amount characteristic with respect to the change of the feature value is partially stored.
  • the current first SOC-V third characteristic or fourth characteristic
  • the parameter can be accurately estimated. That is, first, a plurality of the storage capacity characteristics described above are stored in the table in response to the change according to the deterioration of the feature value. Where a unipolar SOC-V or V-dQ / dV is stored, a counter-polar SOC-V or V-dQ / dV is also stored.
  • SOC-V is a characteristic obtained when charged or discharged by a minute current.
  • the unipolar SOC-V may be SOC-OCP.
  • the SOC-V of the storage element may be SOC-OCV.
  • the counter electrode is a negative electrode and made of, for example, a material such as graphite, hard carbon, or Ti, and there is no potential drop with time, SOC-V or V-dQ / dV is stored as it is when it is new.
  • potential drops occur with time, a plurality of potential drops are stored according to the deterioration.
  • parameters relating to the negative electrode are preferably provided in addition to the ⁇ Qim and the n. Instead of storing a plurality of storage capacity characteristics, a function of the feature value may be stored.
  • the CPU 62 acquires the actually measured charge SOC-V or the actually measured discharge SOC-V (a part of information of the third characteristic or the fourth characteristic) by the constant current of the current predetermined voltage range.
  • the voltage may be OCV.
  • the voltage range may be set from the potential range by determining the unipolar potential range in consideration of the potential of the counter electrode assumed. Feature values are extracted from the measured SOC-V.
  • the characteristic value is the charge quantity or the discharge capacity
  • the numerical value in the predetermined voltage range is extracted. If the characteristic value is the average discharge potential, it is determined by dividing the integrated power (Wh) at the time of discharge by the discharge capacity (Ah).
  • unipolar potential range corresponding to the predetermined voltage range there is a linear relationship between the amount of charge or discharge capacity and the mean discharge potential of the unipolar, and before and after deterioration, the potential difference with the counter electrode (storage element voltage) It is preferable to select the range which does not change. It is preferable to select a wide potential range on the high potential side, in which the charge amount characteristic does not significantly change due to the charge and discharge history.
  • the battery voltage is converted to a unipolar potential in consideration of the currently estimated potential of the counter electrode.
  • FIG. 2 is a conceptual diagram showing the relationship between the potential range of the positive electrode corresponding to the predetermined voltage range and the range of the charge amount corresponding to each deterioration state in each potential range.
  • the potential range becomes narrower in the order of a, b and c.
  • the potential range narrows, the range of the charge amount of electricity narrows. That is, the error increases with the reduction of the potential range to be used.
  • the potential range is wide, it takes time and effort to acquire the amount of charge electricity. Therefore, it is preferable to set an appropriate potential range in consideration of the balance between the estimation accuracy and the ease of measurement.
  • the characteristic value is the charge quantity or the discharge capacity, k can also be determined with high accuracy.
  • the CPU 62 performs optimization calculation based on the measured charge SOC-V (or the measured discharge SOC-V) and the first charge SOC-V (or the first discharge SOC-V) to obtain a parameter.
  • the parameters estimate the internal state. Although Rcha, q, k or k 'is not essential as an initial value of the parameter, it is preferable to input as an initial value from the viewpoint of calculation accuracy and calculation speed.
  • the CPU 62 calculates the second charge SOC-V or the second discharge SOC-V (third characteristic or fourth characteristic) based on the parameter.
  • FIG. 3 shows an example of a power storage device.
  • Power storage device 50 includes a plurality of power storage elements 200, monitoring device 100, and a storage case 300 for storing them.
  • Power storage device 50 may be used as a power source of an electric vehicle (EV) or a plug-in hybrid electric vehicle (PHEV).
  • the storage element 200 is not limited to a square cell, and may be a cylindrical cell or a pouch cell.
  • the monitoring device 100 may be a circuit board disposed to face the plurality of storage elements 200. Monitoring device 100 monitors the state of storage element 200.
  • the monitoring device 100 may be an estimation device. Alternatively, a computer or server wired or wirelessly connected to the monitoring apparatus 100 may execute a method of estimating the internal state based on the information output by the monitoring apparatus 100.
  • FIG. 4 shows another example of the power storage device.
  • the power storage device (hereinafter referred to as a battery module) 1 may be a 12 volt power source or a 48 volt power source suitably mounted on an engine vehicle.
  • FIG. 4 is a perspective view of the battery module 1 for 12 V power supply
  • FIG. 5 is an exploded perspective view of the battery module 1
  • FIG. 6 is a block diagram of the battery module 1.
  • the battery module 1 has a rectangular parallelepiped case 2.
  • a plurality of bus bars 4 a plurality of bus bars 4
  • a BMU (Battery Management Unit) 6 Battery Management Unit
  • the battery 3 includes a rectangular parallelepiped case 31 and a pair of terminals 32 and 32 provided on one side of the case 31 and having different polarities.
  • the case 31 accommodates an electrode body 33 in which a positive electrode plate, a separator, and a negative electrode plate are stacked.
  • At least one of the positive electrode active material of the positive electrode plate of the electrode body 33 and the negative electrode active material of the negative electrode plate has the properties of VF and hysteresis.
  • the positive electrode active material LiMeO 2 -Li 2 MnO 3 solid solution described above, Li 2 O-LiMeO 2 solid solution, Li 3 NbO 4 -LiMeO 2 solid solution, Li 4 WO 5 -LiMeO 2 solid solution, Li 4 TeO 5 -LiMeO 2 solid solution, Li 3 SbO 4 -LiFeO 2 solid solution, Li 2 RuO 3 -LiMeO 2 solid solution, and a Li-excess active material such as Li 2 RuO 3 -Li 2 MeO 3 solid solution.
  • the negative electrode active material examples include hard carbon, metals such as hard carbon, Si, Sn, Cd, Zn, Al, Bi, Pb, Ge, Ag, or alloys thereof, chalcogenides containing these, and the like.
  • chalcogenide is SiO.
  • the technology of the present invention is applicable as long as at least one of the positive electrode active material and the negative electrode active material is included.
  • Case 2 is made of synthetic resin.
  • the case 2 includes a case body 21, a lid 22 closing the opening of the case body 21, a BMU accommodating portion 23 provided on the outer surface of the lid 22, a cover 24 covering the BMU accommodating portion 23, and an inner lid 25 and a partition plate 26.
  • the inner lid 25 and the partition plate 26 may not be provided.
  • the battery 3 is inserted between the partition plates 26 of the case main body 21.
  • a plurality of metal bus bars 4 are mounted on the inner lid 25.
  • the inner cover 25 is disposed on the terminal surface on which the terminals 32 of the battery 3 are provided, and the adjacent terminals 32 of the adjacent batteries 3 are connected by the bus bar 4 and the batteries 3 are connected in series.
  • the BMU accommodating portion 23 has a box shape, and has a rectangular projecting portion 23 a at the center of one long side. A pair of external terminals 5 and 5 made of a metal such as a lead alloy and having different polarities is provided on both sides of the protrusion 23 a in the lid 22.
  • the BMU 6 is formed by mounting an information processing unit 60, a voltage measurement unit 8, and a current measurement unit 9 on a substrate. The BMU 6 is accommodated in the BMU accommodating portion 23 and the BMU accommodating portion 23 is covered with the cover 24, whereby the battery 3 and the BMU 6 are connected.
  • the information processing unit 60 includes a CPU 62 and a memory 63.
  • the memory 63 stores an internal state estimation program (hereinafter referred to as a program) 63a according to the present embodiment and a table 63b.
  • the program 63a is provided in the state of being stored in a computer readable recording medium 70 such as a CD-ROM, a DVD-ROM, a USB memory, etc., and is stored in the memory 63 by being installed in the BMU 6.
  • the program 63a may be obtained from an external computer (not shown) connected to a communication network and stored in the memory 63.
  • the table 63b stores, as an example, a plurality of positive electrode SOC-Vs or their V-dQ / dV according to deterioration, a negative electrode SOC-V or their V-dQ / dVs, or a function of a feature value .
  • the memory 63 and the CPU 62 functioning as the estimation unit are not limited to those installed in the BMU 6. These may be mounted on an external device, and when the actually measured charge SOC-V or the actually measured discharge SOC-V in a predetermined voltage range is acquired, the internal state may be estimated and the result may be passed to the BMU 6.
  • the table 63b may store the SOC-V of the battery 3 or its V-dQ / dV.
  • the storage amount characteristic stored in the table 63b will be described by taking a specific example. For each battery 3, each No. 1 battery was tested under the conditions of voltage range, cycle number, and test temperature shown in Table 1 below. Cycle test was conducted.
  • the conditions of the confirmation test for determining the SOC-V of the positive electrode are as follows. ⁇ Anode: Li metal ⁇ Test rate: charge 0.1CA, discharge 0.1CA ⁇ Test temperature: 25 ° C Thereby, each No.
  • the SOC-V or V-dQ / dV of the positive electrode is determined as the storage capacity characteristic for the test of (1).
  • the storage amount characteristic is stored in the table 63 b in association with the charge amount or discharge capacity of the predetermined voltage range, or the average discharge potential.
  • the charge quantity characteristic of the single pole in a deteriorated state is acquired, and this is arranged in the order of the feature value, and the feature value and the charge quantity characteristic are associated. It is confirmed that the storage capacity characteristic changes continuously and uniquely in the long run even if the usage history is different.
  • the SOC-V of the active material of the opposite electrode is stored in the table 63 b.
  • the CPU 62 executes an internal state estimation process described later according to the program read from the memory 63.
  • the voltage measurement unit 8 is connected to both ends of the battery 3 via a voltage detection line, and measures the voltage of each battery 3 at predetermined time intervals.
  • the current measuring unit 9 measures the current flowing through the battery 3 via the current sensor 7 at predetermined time intervals.
  • the external terminals 5 and 5 of the battery module 1 are connected to a load 11 such as a starter motor for starting the engine and electrical components.
  • the ECU (Electronic Control Unit) 10 is connected to the BMU 6 and the load 11.
  • FIG. 7 is a flowchart showing a procedure of estimation processing of parameters and storage capacity characteristics by the CPU 62.
  • the CPU 62 repeats the processing from S1 at predetermined intervals.
  • the CPU 62 acquires the measured storage amount characteristic by the constant current of the predetermined voltage range (S1).
  • the CPU 62 obtains, for example, SOC-V of charge or discharge.
  • the CPU 62 extracts a feature value from the measured storage capacity characteristic (S2).
  • the CPU 62 extracts the amount of charge in a predetermined voltage range, for example, from the SOC-V.
  • the CPU 62 calculates the first storage capacity characteristic of the battery 3 based on the acquired feature value. (S3).
  • the CPU 62 performs optimization calculation to calculate a parameter based on the actual storage amount characteristic and the first storage amount characteristic (S4).
  • the design variables are the parameters described above, and the objective function is the residual sum of squares f (x).
  • the CPU 62 calculates the storage amount characteristic again based on the parameter, and obtains the second storage amount characteristic (S5).
  • the CPU 62 stores the calculated parameter and the second storage capacity characteristic in the table 63b (S6), and ends the process.
  • the CPU 62 acquires the SOC-V of the battery 3 in the voltage range C1V to C2V (S1).
  • the CPU 62 extracts the charge quantity Q (x) when the voltage range C1V to C2V, that is, the potential range of the positive electrode is P1V to P2V, as the feature value from the SOC-V (S2).
  • This charge quantity of electricity is defined as QinP1-P2V.
  • the V-dQ / dV to be obtained can be obtained by the following equation.
  • dQ / dV dQ / dV (i) ⁇ q + dQ / dV (i ⁇ 1) ⁇ (1 ⁇ q)
  • the dQ / dV of i-1 here no.
  • FIG. 8 is a graph showing the relationship between the degree of deterioration k 'and the amount of charge Qin P1-P2V in the potential range P1-P2V.
  • the horizontal axis is the amount of charge (%), and the vertical axis is the degree of deterioration k '.
  • the solid line graph is at the time of charging, and the broken line graph is at the discharging time. As shown in FIG. 8, a linear correlation is seen between Q and k '.
  • the unipolar potential range has a linear relationship between the amount of charge or discharge capacity and the mean discharge potential of the unipolar electrode, and a range in which the potential difference with the counter electrode does not change before and after deterioration is selected. It may change depending on the state of the opposite electrode. Since the SOC-V of the battery 3 changes according to the state of the counter electrode, it is necessary to estimate the internal state including the state of the counter electrode in order to control the battery 3.
  • the SOC-V of the positive electrode is determined from Q (x) and Rcha, the first SOC-V of the battery 3 is determined, and the optimization calculation is performed based on the measured SOC-V.
  • the CPU 62 converts V-dQ / dV of the positive electrode determined by Q (x) into SOC-V of the positive electrode. Then, the CPU 62 calculates the SOC-V (first storage amount characteristic) of the battery 3 based on the SOC-V of the positive electrode and the SOC-V of the negative electrode (S3). At this time, as parameters, a processing coefficient q 1 of interpolation calculation, SOC-V (k 1 ) of positive electrode, ratio Rcha 1 of charging SOC-V, degree of deterioration k ′ 1 , degree m 1 of contraction of positive electrode, and An amount of deviation ⁇ Qim 1 of the capacity balance is obtained.
  • the CPU 62 calculates a parameter indicating the current internal state of the battery 3.
  • the CPU 62 obtains the residual square sum f (x) of the objective function by the following equation.
  • f (x) ⁇ C1V C2V (V (calculated value)-V (measured value)) 2
  • the voltage range may be set to a range different from C1V to C2V. For each of the first storage amount characteristic estimation and the parameter estimation, it is preferable to select a voltage range in which an optimum result can be obtained.
  • Parameters k 2 , Rcha 2 , k ′ 2 , m 2 , and ⁇ Qim 2 at which the error between the measured SOC-V acquired in S1 and the first SOC-V calculated in S3 is minimized are determined.
  • the CPU 62 calculates the charge SOC-V or the discharge SOC-V (second storage capacity characteristic) of the battery 3 using the obtained parameters (S5).
  • FIG. 9 is a graph showing the discharge DOD-V of the positive electrode, the discharge DOD-V of the battery 3 and the DOD-V of the negative electrode calculated using the parameters determined in S4.
  • the horizontal axis is capacity (%), and the vertical axis is potential or voltage (V).
  • e is the positive electrode DOD-V
  • f is the negative electrode DOD-V
  • g is the battery 3 DOD-V. represented by the graph of the e and f, displacement amount DerutaQim of capacity balance a capacitance difference between the positive and negative electrodes, a DerutaQim 2 obtained.
  • the internal resistance (R) can be calculated from the polarization ( ⁇ V) and the current at the time of acquiring the measured storage amount characteristic.
  • FIG. 11 is a graph showing the result of finding the difference between the measured value and the calculated value of V-dQ / dV of the positive electrode when charge / discharge of 2.0 to 4.6 V per battery 3 is repeated 50 cycles.
  • the horizontal axis is the potential (VvsLi / Li + : Li / Li + ), and the vertical axis is the difference (%).
  • h represents the difference between the charge V-dQ / dV and i the discharge V-dQ / dV.
  • the measured value of V-dQ / dV of the positive electrode is obtained by disassembling the battery 3 which has been repeatedly charged and discharged for 50 cycles, and determining V-dQ / dV of only the positive electrode.
  • Difference ⁇ dQ / dV (calculated value) ⁇ dQ / dV (measured value) ⁇ / (maximum value of absolute value of dQ / dV (measured value)) ⁇ 100 From FIG. 11, it can be seen that the difference is small.
  • FIG. 13 is a graph showing the result of finding the difference between the measured value and the calculated value of V-dQ / dV of the positive electrode when charge / discharge of 2.0 to 4.6 V per battery 3 is repeated 500 cycles.
  • the horizontal axis is the potential (VVvsLi / Li + : Li / Li + ), and the vertical axis is the difference (%).
  • j is the difference between the charge V-dQ / dV and the discharge V-dQ / dV.
  • the measured value of V-dQ / dV of the positive electrode is obtained by disassembling the battery 3 which has been repeatedly charged and discharged for 50 cycles, and determining V-dQ / dV of only the positive electrode.
  • Difference ⁇ dQ / dV (calculated value) ⁇ dQ / dV (measured value) ⁇ / (maximum value of absolute value of dQ / dV (measured value)) ⁇ 100 It can be seen from FIG. 13 that the difference is small.
  • 14 and 15 are flowcharts showing the procedure of the SOC estimation process by the CPU 62.
  • the CPU 62 repeats the processing from S11 at predetermined intervals.
  • the second storage capacity characteristic is SOC-OCV.
  • the voltage with which the oxidation amount and the reduction amount of the reaction generating the hysteresis are small is previously determined by experiments and is set as the threshold value V1.
  • the voltage acquired after the voltage becomes nobler than V1 is set to the upper reference voltage (Vup).
  • Vup is updated when the acquired voltage is greater than the previously acquired voltage.
  • the acquired voltage is set to the lower reference voltage (Vlow) after the voltage goes below V1.
  • Vlow is updated when the acquired voltage is smaller than the previously acquired voltage.
  • the CPU 62 acquires the voltage and current between the terminals of the battery 3 (S11).
  • the threshold V1 and the upper reference voltage Vup are indicated by OCV
  • the correction value to the OCV can be obtained by estimating the voltage when the current is zero using a regression line from a plurality of voltage and current data. If the amount of current flowing through the battery 3 is as small as the dark current (small current), the acquired voltage is regarded as OCV.
  • the CPU 62 determines whether the absolute value of the current is equal to or more than the pause threshold (S12).
  • the rest threshold is set to determine whether the state of the battery 3 is in a charged state, a discharged state, or a rest state.
  • the process proceeds to S22.
  • the CPU 62 determines whether the current is larger than 0 (S13). If the current is greater than 0, it can be determined that the state of the battery 3 is in the charged state. When the CPU 62 determines that the current is not larger than 0 (S13: NO), the process proceeds to S18.
  • the CPU 62 determines whether the acquired voltage is larger than Vup stored in the memory 63 last time (S15). When the CPU 62 determines that the voltage is not larger than the previous Vup (S15: NO), the process proceeds to S17.
  • the memory 63 updates the voltage to Vup (S16).
  • the CPU 62 estimates SO by current integration (S17), and ends the process.
  • the CPU 62 determines whether the voltage is less than V1 (S18). When the CPU 62 determines that the voltage is not less than V1 (S18: NO), the process proceeds to S21. When the CPU 62 determines that the voltage is less than V1 (S18: YES), it determines whether the acquired voltage is smaller than the lower reference voltage Vlow stored in the memory 63 last time (S19). When the CPU 62 determines that the voltage is not smaller than the previous Vlow (S19: NO), the process proceeds to S21. When the CPU 62 determines that the voltage is smaller than the previous Vup (S19: YES), the memory 63 updates the voltage to Vlow (S20) The CPU 62 estimates the SOC by current integration (S21), and ends the process.
  • the CPU 62 determines whether the set time has elapsed (S22).
  • the set time is an experimentally determined time sufficient to regard the acquired voltage as an OCV.
  • the CPU 62 determines whether or not the time has been exceeded, based on the number of acquisition times and the acquisition interval of the current after determining that the apparatus is in the pause state. This makes it possible to estimate the SOC with higher accuracy in the resting state.
  • the CPU 62 estimates the SOC by current integration (S23), and ends the process.
  • the acquired voltage can be regarded as an OCV.
  • the CPU 62 acquires the most recent second storage capacity characteristic from the table 63b (S24).
  • the second storage amount characteristic is the charge SOC-OCV and the discharge SOC-OCV of the battery 3.
  • the CPU 62 calculates a voltage reference storage amount characteristic based on the acquired second storage amount characteristic (S25).
  • the CPU 62 calculates a voltage reference charge SOC-OCV or a voltage reference discharge SOC-OCV based on the charge SOC-OCV, the discharge SOC-OCV, and the Vup.
  • the CPU 62 uses, for example, the charge SOC-OCV, the discharge SOC-OCV, and the Vup in consideration of the oxidation amount and the reduction amount of the reaction that causes the hysteresis, the voltage reference charge SOC-OCV or the voltage reference discharge SOC- Calculate OCV. Vup can be corrected based on the current Rcha.
  • the CPU 62 reads the SOC corresponding to the voltage acquired in S1 in the voltage reference charge SOC-OCV or the voltage reference SOC-OCV to estimate the SOC (S26), and ends the processing.
  • a plurality of storage amount-potential characteristics, storage amount-voltage characteristics, or V-dQ / dV that reflect VF due to deterioration of a single pole are stored, and the storage amount characteristics and characteristics
  • the first storage amount-voltage characteristic is calculated based on the value.
  • the parameter is determined by optimization calculation based on the first storage amount-voltage characteristic and the measured storage amount-voltage characteristic acquired in the predetermined voltage range.
  • the parameter can detect the unipolar state and the counterpolar state corresponding to VF in a simple manner. Therefore, the current internal state of the battery 3 can be estimated with high accuracy.
  • the number of data of the storage amount characteristic stored in the table 63b may be small.
  • the internal state can be monitored with high accuracy even in a complicated use environment.
  • the voltage reference storage voltage-voltage characteristics are accurately based on the second storage amount-potential characteristic according to the current deterioration state of the single pole and the charge / discharge history of the storage element. It can be estimated.
  • the amount of stored electricity can be relatively easily estimated by using the knowledge of the behavior of the hysteresis in combination with the storage element including the active material having a high capacity and the property of VF. Since the voltage is used, the storage amount is not limited to the SOC, and the amount of current energy stored in the storage element, such as the amount of power, can be estimated. Based on the charge / discharge characteristics, dischargeable energy up to SOC 0% and charge energy required up to SOC 100% can be predicted.
  • the storage capacity-potential characteristic or V-dQ / may also be applied to the case where the negative electrode includes the active material having VF and hysteresis.
  • dV can be estimated.
  • the estimation of the storage amount by voltage reference is not limited to the case of stopping, but may be performed in real time during charging or discharging.
  • the current OCV is calculated from the acquired voltage and current.
  • the calculation of the OCV can be obtained by, for example, estimating the voltage when the current is zero using regression lines from data of a plurality of voltages and currents. Further, when the current is small as dark current, the acquired voltage V can be read as OCV.
  • the estimation device and the storage amount estimation device according to the present invention are not limited to the case of being applied to a lithium ion secondary battery for vehicle use, and may be applied to other power storage devices such as railway regenerative power storage devices and solar power generation systems. Applicable In a power storage device in which a minute current flows, the voltage between the positive electrode terminal and the negative electrode terminal of the power storage element can be regarded as OCV.
  • the storage element is not limited to a lithium ion secondary battery, and may be another secondary battery or an electrochemical cell having VF and hysteresis characteristics.
  • a CMU Cell Monitoring Unit
  • the storage amount estimation device may be part of a storage module in which the monitoring device 100 and the like are incorporated.
  • the storage amount estimation device may be configured separately from the storage element and the storage module, and may be connected to the storage module including the storage element for which the storage amount is to be estimated when the storage amount is estimated.
  • the heat storage amount estimation device may remotely monitor the storage element or the storage module.
  • the present invention is applicable to estimation of the internal state of a storage element such as a lithium ion secondary battery.

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)
  • Power Engineering (AREA)
  • Secondary Cells (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Abstract

推定装置６は、蓄電素子３の蓄電量－電圧充電特性である第３特性及び／又は蓄電量－電圧放電特性である第４特性の一部の情報を取得する取得部６２と、充放電の繰り返しにより変化する特徴値の変化に対応して、第１特性、第２特性、蓄電量－電圧充電特性である第３特性、蓄電量－電圧放電特性である第４特性、及び、Ｖ－ｄＱ／ｄＶの少なくともいずれかである蓄電量特性を複数記憶し、又は前記特徴値の関数として記憶する記憶部６３と、前記情報及び前記蓄電量特性に基づいて、蓄電素子３の内部状態を推定する第１推定部６２とを備える。

Description

推定装置、蓄電装置、推定方法、及びコンピュータプログラム
　本発明は、推定装置、該推定装置を含む蓄電装置、内部状態の推定方法、及びコンピュータプログラムに関する。
　電気自動車、ハイブリッド車等に用いられる車両用の二次電池や、電力貯蔵装置、太陽光発電システム等に用いられる産業用の二次電池においては、高容量化が求められている。これまで様々な検討と改良が行われてきて、電極構造等の改良のみで更なる高容量化を実現することは困難な傾向にある。その為、現行の材料より高容量である正極材料の開発が進められている。
　従来、リチウムイオン二次電池等の非水電解質二次電池用の正極活物質として、α－ＮａＦｅＯ２型結晶構造を有するリチウム遷移金属複合酸化物が検討され、ＬｉＣｏＯ２を用いた非水電解質二次電池が広く実用化されていた。ＬｉＣｏＯ２の放電容量は１２０～１３０ｍＡｈ／ｇ程度であった。
　リチウム遷移金属複合酸化物をＬｉＭｅＯ２（Ｍｅは遷移金属）で表したとき、ＭｅとしてＭｎを用いることが望まれてきた。ＭｅとしてＭｎを含有させた場合、Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５を超える場合には、充電をするとスピネル型へと構造変化が起こり、結晶構造が維持できない為、充放電サイクル性能が著しく劣る。
　Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５以下であり、Ｍｅに対するＬｉのモル比Ｌｉ／Ｍｅが略１であるＬｉＭｅＯ２型活物質が種々提案され、実用化されている。リチウム遷移金属複合酸化物であるＬｉＮｉ１／２Ｍｎ１／２Ｏ２及びＬｉＮｉ１／３Ｃｏ１／３Ｍｎ１／３Ｏ２等を含有する正極活物質は１５０～１８０ｍＡｈ／ｇの放電容量を有する。
　ＬｉＭｅＯ２型活物質に対し、Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５を超え、遷移金属（Ｍｅ）の比率に対するＬｉの組成比率Ｌｉ／Ｍｅが１より大きいリチウム遷移金属複合酸化物を含む、いわゆるリチウム過剰型活物質も知られている。
　上述の高容量の正極材料として、リチウム過剰型であるＬｉ2 ＭｎＯ3 系の活物質が検討されている。この材料は、同一のＳＯＣ（State Of Charge）に対して、充電時及び放電時の各ＳＯＣ－ＯＣＶ（Open Circuit Voltage）間に、電圧や電気化学的特性の差が生じる、ヒステリシスという性質を有する。
　ヒステリシスを有する場合、ＳＯＣに対して電圧が一義的に決まらない為、ＳＯＣ－ＯＣＶに基づいてＳＯＣを推定するＯＣＶ法によるＳＯＣの推定は困難である。ＳＯＣ－ＯＣＶ曲線が一義的に決まらない為、ある時点での放電可能エネルギーを予測することも困難である。
　リチウム過剰型の材料は、充放電の繰り返しにより正極のＳＯＣ－ＯＣＰ（Open Circuit Potential）曲線が略全域に亘って変化する、電位降下（Voltage Fade、以下ＶＦという）という性質を有する。平均放電電位の値が減少するため、現時点のＳＯＨ（State of Health）として放電可能容量だけでなく、放電可能電力量を推定する必要がある。直近の充放電の履歴が同一であっても、劣化により単極のＳＯＣ－ＯＣＰ曲線に基づく電池セル（以下、単に「セル」ともいう。）のＳＯＣ－ＯＣＶ曲線形状が大幅に変わるため、ＯＣＶ法は採用できない。直近の充放電履歴が同一の条件とは、例えば、完全放電状態を経由した後の充電が挙げられる。完全放電状態を経由した後の充電において、劣化に応じて単極のＳＯＣ－ＯＣＰ曲線が変わる為、セルのＳＯＣ－ＯＣＶ曲線形状が大幅に変わってしまう。
　二次電池の充放電電流を積算して決定する電流積算法によりＳＯＣを推定する場合、電流積算が長期継続されると、電流センサの計測誤差が蓄積する。また、電池容量は経時的に小さくなる。その為、電流積算法によって推定されるＳＯＣは、その推定誤差が経時的に大きくなる。従来、電流積算を長期継続した場合にＯＣＶ法によりＳＯＣを推定して、誤差の蓄積をリセットするＯＣＶリセットが行われている。
　ＶＦ及びヒステリシスを有する電極材料を用いた蓄電素子においても、電流積算を継続すると誤差が蓄積する。しかし、ＳＯＣに対して電圧が一義的に決まらない為、ＯＣＶリセットを行うことは困難である。
　また、対極の状態によっても蓄電素子のＳＯＣ－ＯＣＶが変化するので、蓄電素子を制御するためには対極の状態も含めた蓄電素子の内部状態を推定する必要がある。
　即ち、現在の満充電時及び完全放電時の正極のＳＯＣ－ＯＣＰ特性、充電状態（正極のＳＯＣ－ＯＣＰの割合）、正極の収縮の度合、負極の収縮の度合、正極及び負極間の容量バランスのずれ量、及び分極等のパラメータを推定する必要がある。そして、正確にＳＯＣを推定する必要がある。
　現行の非水電解質二次電池のＳＯＨ及びＳＯＣの推定技術では、ＶＦ及びヒステリシスの性質を有する活物質を使いこなすことは困難である。
　特許文献１に開示の劣化判定装置の測定部は、リチウムイオン二次電池の容量対開放電圧特性を測定する。判定部は、開放電圧特性を特定するためのパラメータを設定可能であり、測定部によって測定された開放電圧特性と略一致する開放電圧特性を特定するパラメータを用いて摩耗およびリチウム析出による劣化状態を判定する。
特開２０１１－２２０９１７号公報
　特許文献１の劣化判定装置は、摩耗およびリチウム析出による劣化状態を判定するものであり、ＶＦを生じる活物質を含む蓄電素子の内部状態を推定することはできない。
　本発明は、蓄電量－電位特性が充放電の繰り返しにより変化する単極を持つ蓄電素子の内部状態を推定できる推定装置、該推定装置を備える蓄電装置、内部状態の推定方法、及びコンピュータプログラムを提供することを目的とする。
　ここで、蓄電量とは、ＳＯＣ、電力放出可能量等を意味する。
　本発明に係る推定装置は、蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定する。推定装置は、前記蓄電素子の蓄電量－電圧充電特性である第３特性及び蓄電量－電圧放電特性である第４特性の一部の情報を取得する取得部と、充放電の繰り返しにより変化する特徴値の変化に対応して、第１特性、第２特性、第３特性、第４特性、及び、Ｖ－ｄＱ／ｄＶの少なくともいずれかである蓄電量特性を複数記憶し、又は前記特徴値の関数として記憶する記憶部と、前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する第１推定部とを備える。
　上記構成により、蓄電量－電位特性が充放電の繰り返しにより変化する単極を持つ蓄電素子の内部状態を良好に推定できる。
正極のＳＯＣ－ＯＣＰの一例を示すグラフである。 所定電圧範囲に対応する正極の電位範囲と、各電位範囲における、各劣化状態に対応する充電電気量の範囲との関係を示す概念図である。 蓄電装置の一例を示す斜視図である。 蓄電装置の他の例を示す斜視図である。 電池モジュールの分解斜視図である。 電池モジュールのブロック図である。 ＣＰＵによるパラメータ及び蓄電量特性の推定処理の手順を示すフローチャートである。 電位範囲Ｐ１－Ｐ２Ｖにおける劣化程度ｋ′と充電電気量Ｑとの関係を示すグラフである。 Ｓ４で求めたパラメータを用いて算出した正極の放電ＤＯＤ－Ｖ、電池の放電ＤＯＤ－Ｖ、及び負極のＤＯＤ－Ｖを示すグラフである。 電池につき２．０－４．６Ｖの充放電を５０サイクル繰り返した場合に、実測したＳＯＣ－Ｖと、上述の方法によりパラメータを求めて算出した第２ＳＯＣ－Ｖとの差異を求めた結果を示すグラフである。 電池につき２．０－４．６Ｖの充放電を５０サイクル繰り返した場合の、正極のＶ－ｄＱ／ｄＶの実測値と算出値との差異を求めた結果を示すグラフである。 電池につき２．０－４．６Ｖの充放電を５００サイクル繰り返した場合に、実測したＳＯＣ－Ｖと、上述の方法によりパラメータを求めて算出した第２ＳＯＣ－Ｖとの差異を求めた結果を示すグラフである。 電池につき２．０－４．６Ｖの充放電を５００サイクル繰り返した場合の、正極のＶ－ｄＱ／ｄＶの実測値と算出値との差異を求めた結果を示すグラフである。 ＣＰＵによるＳＯＣ推定処理の手順を示すフローチャートである。 ＣＰＵによるＳＯＣ推定処理の手順を示すフローチャートである。
（実施形態の概要）
　実施形態に係る推定装置は、蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定する。推定装置は、前記蓄電素子の蓄電量－電圧充電特性である第３特性及び／又は蓄電量－電圧放電特性である第４特性の一部の情報を取得する取得部と、充放電の繰り返しにより変化する特徴値の変化に対応して、第１特性、第２特性、第３特性、第４特性、及び、Ｖ－ｄＱ／ｄＶの少なくともいずれかである蓄電量特性を複数記憶し、又は前記特徴値の関数として記憶する記憶部と、前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する第１推定部とを備える。
　ここで、Ｖ－ｄＱ／ｄＶは、前記単極の電位又は前記蓄電素子の電圧Ｖと、充電電気量若しくは放電容量Ｑを前記Ｖで微分した微分値ｄＱ／ｄＶとの関係である。
　上記構成によれば、複数記憶された蓄電量特性は、単極のＶＦに基づく劣化を特徴値に応じて反映する。前記情報は特徴値に対応する。前記情報と記憶された蓄電量特性とに基づいて、簡便な方法で、対極の状態も含め、蓄電素子の現在の内部状態が良好に推定できる。
　上述の推定装置において、前記取得部は、前記蓄電素子の所定電圧範囲における、第３特性及び／又は第４特性の情報を取得し、前記情報から特徴値を求め、該特徴値に基づいて、前記蓄電量特性又は前記関数を参照して、前記情報を取得した時点の第３特性又は第４特性を推定する第２推定部を備え、前記第１推定部は、前記情報と、推定した前記第３特性又は前記第４特性とに基づき、前記蓄電素子の内部状態を示すパラメータを推定してもよい。
　上記構成によれば、特徴値及び記憶された蓄電量特性に基づいて算出した蓄電量特性、又は前記関数に基づいて算出した蓄電量特性、並びに所定電圧範囲で実測した前記情報に基づいて、例えば回帰計算によりパラメータが求められる。簡便な方法で、対極の状態も含め、蓄電素子の現在の内部状態が高精度に推定できる。
　複雑な使用環境下においても、高精度に蓄電素子の内部状態を監視することができる。
　上述の推定装置において、前記特徴値は、前記所定電圧範囲における充電電気量、放電容量、及び、平均放電電位の少なくともいずれかであってもよい。
　充電電気量又は放電容量と平均放電電位との間に直線関係があり、劣化の前後で、対極との電位差（蓄電素子の電圧）が変わらない電位範囲に対応する蓄電素子の電圧範囲が、前記所定の電圧範囲とされる。所定の電圧範囲における充電電気量又は放電容量が特徴値として用いられる。劣化の程度に応じた該特徴値に対応付けて蓄電量特性が複数記憶され、又は特徴値の関数として記憶される。これを用いて精度良く、現在の蓄電量特性が推定できる。平均放電電位を用いた場合も、同様に蓄電量特性が精度良く推定できる。現在の蓄電量－電位特性、蓄電量－電圧特性、又はＶ－ｄＱ／ｄＶに基づき、精度良く現在の第３特性、第４特性を推定できる。
　上述の推定装置において、前記パラメータは、第１特性、第２特性、第１特性の割合、前記単極の収縮の度合、対極の収縮の度合、前記単極及び前記対極間の容量バランスのずれ量、及び分極から選択される少なくとも１種であってもよい。
　但し、前記第１特性の割合は下記式で表される。
　（第１特性の割合）＝（前記情報を取得したときに、前記記憶部に記憶されている前記蓄電素子の到達最大電圧に対応する前記単極の電位における、第１特性及び第２特性間の蓄電量の差）／（前記蓄電量の差の最大値）
　上述のパラメータにより、劣化によりＶＦが生じる単極の状態、及び対極の状態を検知することができ、蓄電素子の現在の内部状態が良好に検知できる。
　上述の推定装置において、前記第２推定部は、前記Ｖ－ｄＱ／ｄＶを参照して、内挿計算によりＶ－ｄＱ／ｄＶを推定し、前記パラメータは、前記内挿計算に必要な処理係数を含んでもよい。
　上記構成によれば、内挿計算により現在のＶ－ｄＱ／ｄＶを推定した後、第３特性又は第４特性を求めるので、算出精度及び算出速度が良好である。そして、内挿計算に必要な処理係数をパラメータとすることで、回帰計算により求められたパラメータにより、第３特性又は第４特性を精度良く補正することができる。
　上述の推定装置において、前記第１推定部は、最適化計算により、前記情報と、前記第２推定部により推定した前記第３特性又は前記第４特性との差が最小になるように、前記パラメータを求めてもよい。
　上記構成によれば、最適化計算により現在のパラメータを精度良く求めることができ、精度良く現在の蓄電素子の内部状態が検知できる。
　上述の推定装置において、前記第１推定部により推定した前記パラメータを用いて、第３特性又は第４特性を推定する第３推定部を備えてもよい。
　上記構成によれば、現在の蓄電素子の内部状態を精度良く表すパラメータを用い第３特性又は第４特性を補正して推定することができる。
　上述の推定装置において、前記活物質は、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、第３特性及び第４特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さく、前記第３推定部により推定した第３特性又は第４特性、及び前記蓄電素子の充放電の履歴に基づいて、取得した電圧により蓄電量を推定するときの参照のための第３特性又は第４特性を推定する第４推定部を備えてもよい。
　上記構成によれば、現在の蓄電素子の内部状態に応じた第３特性又は第４特性、パラメータ、及び蓄電素子の充放電の履歴に基づき、劣化によるＶＦの影響に加え、ヒステリシスの影響を精度良く反映して、参照のための第３特性又は第４特性が推定できる。
　上述の推定装置において、前記充放電の履歴、参照のための前記第３特性又は前記第４特性、並びに取得した電圧に基づいて、蓄電量を推定する第５推定部を備えてもよい。
　上記構成によれば、ＶＦの性質を有し、蓄電量－電圧特性がヒステリシスを示す活物質を有する蓄電素子の蓄電量が良好に容易に推定できる。
　電圧を用いるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量が推定できる。充放電特性に基づいて、ＳＯＣ０％までの放電可能なエネルギー、及びＳＯＣ１００％までに必要な充電エネルギーを予測することができる。現在の残存電力量と貯蔵可能電力量とが推定できる。
　従って、複数の蓄電素子を用いる場合のバランシング、回生受け入れの制御、蓄電素子を車載した場合の走行距離の推定等を精度良く行うことができる。
　実施形態に係る蓄電装置は、蓄電素子と、上述のいずれかの推定装置とを備える。
　上記構成においては、蓄電素子の蓄電量が、複雑な使用環境下においても、精度良く推定できる。
　実施形態に係る推定方法は、蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定する。推定方法は、充放電の繰り返しにより変化する特徴値の変化に対応して、第１特性、第２特性、前記蓄電素子の蓄電量－電圧充電特性である第３特性、蓄電量－電圧放電特性である第４特性、及び、Ｖ－ｄＱ／ｄＶ、の少なくともいずれかである蓄電量特性を、記憶部に、複数記憶し、又は前記特徴値の関数として記憶してあり、前記蓄電素子の第３特性及び／又は第４特性の一部の情報を取得し、前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する。
　上記構成によれば、簡便な方法で、複雑な使用環境下における場合でも、単極及び対極の状態を検知して、高精度に蓄電素子の現在の内部状態が推定できる。
　実施形態に係るコンピュータプログラムは、蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定するコンピュータに、前記蓄電素子の蓄電量－電圧充電特性である第３特性及び／又は蓄電量－電圧放電特性である第４特性の一部の情報を取得し、第１特性、第２特性、第３特性、第４特性、及び、Ｖ－ｄＱ／ｄＶ、の少なくともいずれかである蓄電量特性が、特徴値の変化に対応して、複数記憶されたテーブルを参照し、又は前記特徴値の関数を参照し、前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する、処理を実行させる。
　以下、本発明を、実施形態を示す図面に基づいて具体的に説明する。
　上述したように、本実施形態に係る蓄電素子の電極体の単極は、ＶＦの性質を有し、蓄電量－電位特性がヒステリシスを有する活物質を含む。
　活物質がＶＦの性質を有する場合、充放電の繰り返しにより、単極のＳＯＣ－ＯＣＰ曲線、及び、蓄電素子のＳＯＣ－ＯＣＶ曲線の形状が変化する。この活物質を含む蓄電素子は、微小電流を流して、完全放電状態から満充電状態に充電したとき（フル充電）、及び満充電状態から完全放電状態に放電したとき（フル放電）のＳＯＣ－Ｖ曲線間の最大の電位差が１００ｍＶ以上であるヒステリシスを有する。
　図１は、正極のＳＯＣ－ＯＣＰの一例を示すグラフである。横軸はＳＯＣ（％）、縦軸はＯＣＰとしての電位Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+）である。劣化前の充放電曲線を破線で、劣化後の充放電曲線を実線で示す。
　図１に示すように、劣化によりＶＦが生じ、充放電曲線は下側にシフトする。
　ＶＦの性質を有さない活物質の場合、前記ヒステリシスは有さず、単極のＳＯＣ－ＯＣＰ曲線は充放電の繰り返しによっては変化しない。単極の劣化（曲線の縮小）又は両極間の容量バランスのずれ量の拡大によって、蓄電素子のＳＯＣ－ＯＣＶ曲線の形状は充放電の繰り返しにより変化する。
　本実施の形態においては、現在の蓄電素子、単極、及び対極の内部状態を推定する。内部状態を示すパラメータとしては、単極の蓄電量－電位充電特性（第１特性）又はそのＶ－ｄＱ／ｄＶ、単極の蓄電量－電位放電特性（第２特性）又はそのＶ－ｄＱ／ｄＶ（以上、ｋという）、単極の蓄電量－電位充電特性（第１特性）の割合（後述するＲcha）、正極の収縮の度合（ｍ）、負極の収縮の度合（ｎ）、正極及び負極間の容量バランスのずれ量（ΔＱim）、分極（ΔＶ）、及び後述する内挿計算の処理係数（ｑ）から選択される少なくとも１種以上が挙げられる。前記ｋは劣化の程度ｋ′（劣化がない状態を１．０とする）で表すこともできる。ｍ及びｎは夫々、正極の蓄電量－Ｖ曲線、負極の蓄電量－Ｖ曲線のｘ軸方向の収縮の度合を示す。ΔＶは蓄電量－Ｖ曲線のｙ軸方向の収縮の度合を示す。
　以下、蓄電量がＳＯＣである場合につき説明する。
（１）まず、後述する情報処理部６０は、劣化に応じて複数の蓄電量特性を記憶し、又は特徴値の関数として記憶する。
　蓄電量特性としては、単極の充電ＳＯＣ－Ｖ（第１特性）又は放電ＳＯＣ－Ｖ（第２特性）、蓄電素子の充電ＳＯＣ－Ｖ（第３特性）又は放電ＳＯＣ－Ｖ（第４特性）、単極の充電Ｖ－ｄＱ／ｄＶ又は放電Ｖ－ｄＱ／ｄＶ、蓄電素子の充電Ｖ－ｄＱ／ｄＶ又は放電Ｖ－ｄＱ／ｄＶ、等が挙げられる。
　充放電の繰り返しにより変化する特徴値と、上述の蓄電量特性との間に、相関関係がある。
　特徴値としては、充電電気量、放電容量、平均放電電位、ｄＱ／ｄＶ又はｄＶ／ｄＱピークの絶対値またはそのピーク強度比等が挙げられる。ｄＱ／ｄＶのピーク強度比の例としては、所定の電位範囲ａにおけるｄＱ／ｄＶのピーク積分強度と電位範囲ｂにおけるピーク積分強度との比が挙げられる。この比を取ることで、正極の収縮の度合ｍの減少に関わらず、現時点の蓄電量特性を推定することが可能である。
　中でも、充電電気量、放電容量が好ましい。
　ＶＦの性質を有する活物質の場合、特徴値の変化（劣化）に対応して、充放電曲線形状が連続的及び一義的に変化すると考えられる。ＬｉＭｅＯ２-Ｌｉ２ＭｎＯ３系の活物質に関して、充放電の繰り返しに伴い、結晶構造が変化することが報告されている（Journal of Power Sources, vol.229(2013), pp239-248）。結晶構造の変化に伴い、充放電曲線の形状が変化すると考えられる。短期的で一つの温度水準の充放電サイクルにおいて、結晶構造の変化は連続的に生じることが論文の結果から示唆される。そして、結晶構造が層状からスピネル類縁結晶に変化したという報告から、変化の仕方は一通りであると推察される。即ち、短期的で一つの温度水準において、結晶構造が連続的かつ一義的に変化している。この報告から、長期的かつ、いかなる使用履歴においても、結晶構造の変化に従って充放電曲線形状が連続的かつ一義的に変化すると、本願発明者等は考えた。後述する実験結果から、長期的に、使用履歴が異なっても、充放電曲線形状が連続的かつ一義的に変化することが確認された。
　ＶＦの性質を有さない活物質の場合、単極の充放電曲線形状は充放電の繰り返しによっては変化しない。上述のように単極の劣化又は容量バランスのずれ量の拡大によって、個別に、即ち非一義的に、セルの充放電曲線形状は充放電の繰り返しにより変化する。
　本実施の形態の場合、特徴値の変化に対して、単極の蓄電量特性が連続的かつ一義的に変化するので、特徴値の変化に対する蓄電量特性の変化の推移を一部記憶しておくことで、現在の第１ＳＯＣ－Ｖ（第３特性又は第４特性）が精度良く推定でき、これに基づき、パラメータが精度良く推定できる。
　即ち、まず、特徴値の劣化に従う変化に対応して、上述のいずれかの蓄電量特性が複数、テーブルに記憶される。単極のＳＯＣ－Ｖ又はＶ－ｄＱ／ｄＶが記憶される場合、対極のＳＯＣ－Ｖ又はＶ－ｄＱ／ｄＶも記憶される。ＳＯＣ－Ｖは、微小電流により充電又は放電したときに得られる特性である。単極のＳＯＣ－ＶはＳＯＣ－ＯＣＰであってもよい。蓄電素子のＳＯＣ－ＶはＳＯＣ－ＯＣＶであってもよい。
　対極が負極であり、例えばグラファイト、ハードカーボン、Ｔｉ等の材料からなり、経時的な電位降下が生じない場合、ＳＯＣ－Ｖ又はＶ－ｄＱ／ｄＶは新品時のものが記憶される。経時的に電位降下が生じる場合、劣化に応じて複数記憶される。また、負極に係るパラメータも、前記ΔＱim、前記ｎ以外に設けるのが好ましい。
　複数の蓄電量特性を記憶する代わりに、特徴値の関数が記憶されてもよい。
（２）ＣＰＵ６２は、現在の所定電圧範囲の定電流による実測充電ＳＯＣ－Ｖ又は実測放電ＳＯＣ－Ｖ（第３特性又は第４特性の一部の情報）を取得する。電圧はＯＣＶであってもよい。電圧範囲は、想定される対極の電位を考慮して単極の電位範囲を求め、該電位範囲から設定してもよい。
　実測ＳＯＣ－Ｖから特徴値が抽出される。特徴値が充電電気量又は放電容量である場合、前記所定電圧範囲における数値が抽出されることになる。特徴値が、平均放電電位である場合、放電時の積算電力（Ｗｈ）を放電容量（Ａｈ）で除して求める。
　前記所定電圧範囲に対応する単極の電位範囲として、充電電気量又は放電容量と単極の平均放電電位との間に直線関係があり、劣化の前後で、対極との電位差（蓄電素子電圧）が変わらない範囲を選択するのが好ましい。充放電履歴によって蓄電量特性が大きく変化しない、高電位側における広い電位範囲を選択するのが好ましい。現時点において推定される対極の電位を考慮して、電池電圧が単極の電位に変換される。
　図２は、前記所定電圧範囲に対応する正極の電位範囲と、各電位範囲における、各劣化状態に対応する充電電気量の範囲との関係を示す概念図である。電位範囲はａ、ｂ、ｃの順に狭くなる。電位範囲が狭くなった場合、充電電気量の範囲が狭くなる。即ち、使用する電位範囲の縮小に伴い、誤差が増大する。一方、電位範囲が広い場合、充電電気量の取得に時間及び労力が要される。従って、推定精度及び測定の容易さのバランスを考慮して、適切な電位範囲を設定するのが好ましい。
（３）ＣＰＵ６２は、現在の特徴値に基づき、記憶された蓄電量特性又は関数を参照して、蓄電素子の第１充電ＳＯＣ－Ｖ又は第１放電ＳＯＣ－Ｖ（第３特性又は第４特性）を算出する。
　まず、記憶された蓄電量特性又は関数に基づき前記ｋ（単極のＳＯＣ－Ｖ又はＶ－ｄＱ／ｄＶ）を求めた後、蓄電素子の第１ＳＯＣ－Ｖを求めるのが好ましい。この場合、計算精度及び計算速度が良好である。特徴値が充電電気量又は放電容量である場合、精度良くｋを求めることもできる。
（４）ＣＰＵ６２は、実測充電ＳＯＣ－Ｖ（又は実測放電ＳＯＣ－Ｖ）と、第１充電ＳＯＣ－Ｖ（又は第１放電ＳＯＣ－Ｖ）とに基づき、最適化計算を行ってパラメータを求める。パラメータにより内部状態が推定される。パラメータの初期値として、Ｒcha、ｑ、ｋ又はｋ′は必須ではないが、計算精度及び計算速度の観点から初期値として入力するのが好ましい。
（５）ＣＰＵ６２は、パラメータに基づいて第２充電ＳＯＣ－Ｖ又は第２放電ＳＯＣ－Ｖ（第３特性又は第４特性）を算出する。
（実施形態１）
　以下、実施形態１として、車両に搭載される蓄電装置が例示される。
　図３は、蓄電装置の一例を示す。蓄電装置５０は、複数の蓄電素子２００と、監視装置１００と、それらを収容する収容ケース３００とを備えている。蓄電装置５０は、電気自動車（ＥＶ）や、プラグインハイブリッド電気自動車（ＰＨＥＶ）の動力源として使用されてもよい。
　蓄電素子２００は、角形セルに限定されず、円筒形セルやパウチセルであってもよい。
　監視装置１００は、複数の蓄電素子２００と対向して配置される回路基板であってもよい。監視装置１００は、蓄電素子２００の状態を監視する。監視装置１００が、推定装置であってもよい。代替的に、監視装置１００と有線接続または無線接続されるコンピュータやサーバが、監視装置１００が出力する情報に基づいて内部状態を推定する方法を実行してもよい。
　図４は、蓄電装置の他の例を示す。蓄電装置（以下、電池モジュールという）１は、エンジン車両に好適に搭載される、１２ボルト電源や、４８ボルト電源であってもよい。図４は１２Ｖ電源用の電池モジュール１の斜視図、図５は電池モジュール１の分解斜視図、図６は電池モジュール１のブロック図である。
　電池モジュール１は直方体状のケース２を有する。ケース２に複数のリチウムイオン二次電池（以下、電池という）３、複数のバスバー４、ＢＭＵ（Battery Management Unit）６、電流センサ７が収容される。
　電池３は、直方体状のケース３１と、ケース３１の一側面に設けられた、極性が異なる一対の端子３２，３２とを備える。ケース３１には、正極板、セパレータ、及び負極板を積層した電極体３３が収容されている。
　電極体３３の正極板が有する正極活物質及び負極板が有する負極活物質の少なくとも一方は、ＶＦ及びヒステリシスの性質を有する。
　正極活物質としては、上述のＬｉＭｅＯ２-Ｌｉ２ＭｎＯ３固溶体、Ｌｉ2Ｏ－ＬｉＭｅＯ2固溶体、Ｌｉ3ＮｂＯ4 －ＬｉＭｅＯ2固溶体、Ｌｉ4 ＷＯ5 －ＬｉＭｅＯ2固溶体、Ｌｉ4 ＴｅＯ5 －ＬｉＭｅＯ2固溶体、Ｌｉ3ＳｂＯ4 －ＬｉＦｅＯ2固溶体、Ｌｉ2ＲｕＯ3 －ＬｉＭｅＯ2固溶体、Ｌｉ2ＲｕＯ3 －Ｌｉ2 ＭｅＯ3 固溶体等のＬｉ過剰型活物質が挙げられる。負極活物質としては、ハードカーボン、Ｓｉ、Ｓｎ、Ｃｄ、Ｚｎ、Ａｌ、Ｂｉ、Ｐｂ、Ｇｅ、Ａｇ等の金属若しくは合金、又はこれらを含むカルコゲン化物等が挙げられる。カルコゲン化物の一例として、ＳｉＯが挙げられる。本発明の技術は、これらの正極活物質及び負極活物質の少なくとも一方が含まれていれば適用可能である。
　ケース２は合成樹脂製である。ケース２は、ケース本体２１と、ケース本体２１の開口部を閉塞する蓋部２２と、蓋部２２の外面に設けられたＢＭＵ収容部２３と、ＢＭＵ収容部２３を覆うカバー２４と、中蓋２５と、仕切り板２６とを備える。中蓋２５や仕切り板２６は、設けられなくてもよい。
　ケース本体２１の各仕切り板２６の間に、電池３が挿入されている。
　中蓋２５には、複数の金属製のバスバー４が載置されている。電池３の端子３２が設けられている端子面に中蓋２５が配置されて、隣り合う電池３の隣り合う端子３２がバスバー４により接続され、電池３が直列に接続されている。
　ＢＭＵ収容部２３は箱状をなし、一長側面の中央部に、外側に角型に突出した突出部２３ａを有する。蓋部２２における突出部２３ａの両側には、鉛合金等の金属製で、極性が異なる一対の外部端子５，５が設けられている。ＢＭＵ６は、基板に情報処理部６０、電圧計測部８、及び電流計測部９を実装してなる。ＢＭＵ収容部２３にＢＭＵ６を収容し、カバー２４によりＢＭＵ収容部２３を覆うことにより、電池３とＢＭＵ６とが接続される。
　図６に示すように、情報処理部６０は、ＣＰＵ６２と、メモリ６３とを備える。
　メモリ６３には、本実施形態に係る内部状態の推定プログラム（以下、プログラムという）６３ａと、テーブル６３ｂとが記憶されている。プログラム６３ａは、例えば、ＣＤ－ＲＯＭやＤＶＤ－ＲＯＭ、ＵＳＢメモリ等のコンピュータ読み取り可能な記録媒体７０に格納された状態で提供され、ＢＭＵ６にインストールすることによりメモリ６３に格納される。代替的に、通信網に接続されている図示しない外部コンピュータからプログラム６３ａを取得し、メモリ６３に記憶させてもよい。
　テーブル６３ｂには、一例として、劣化に応じた複数の正極のＳＯＣ－Ｖ又はそのＶ－ｄＱ／ｄＶ、負極のＳＯＣ－Ｖ又はそのＶ－ｄＱ／ｄＶ、又は特徴値の関数が格納されている。メモリ６３、及び推定部（特に第１推定部）として機能するＣＰＵ６２は、ＢＭＵ６に搭載されている場合に限定されない。これらを外部装置に搭載し、所定電圧範囲の実測充電ＳＯＣ－Ｖ又は実測放電ＳＯＣ－Ｖを取得したときに、内部状態を推定し、ＢＭＵ６へ結果を渡すことにしてもよい。
　テーブル６３ｂには、電池３のＳＯＣ－Ｖ又はそのＶ－ｄＱ／ｄＶを格納することにしてもよい。
　テーブル６３ｂに記憶する蓄電量特性について、具体例を挙げて説明する。
　電池３につき、下記の表１に示す電圧範囲、サイクル数、及び試験温度の条件で各Ｎｏ．のサイクル試験を行った。
Figure JPOXMLDOC01-appb-T000001
　充放電の条件は、以下の通りである。
・負極：グラファイト
・試験レート：充電０．５ＣＡ、放電１．０ＣＡ
　正極のＳＯＣ－Ｖを求めるための確認試験の条件は、以下の通りである。
・負極：Ｌｉ金属 
・試験レート：充電０．１ＣＡ、放電０．１ＣＡ
・試験温度：２５℃
　これにより、各Ｎｏ．の試験につき、蓄電量特性として、正極のＳＯＣ－Ｖ、又はＶ－ｄＱ／ｄＶが求められる。この蓄電量特性が所定電圧範囲の充電電気量若しくは放電容量、又は平均放電電位と対応付けて、テーブル６３ｂに格納される。各試験により、劣化した状態の単極の蓄電量特性を取得し、これを特徴値の順に並べて、特徴値と蓄電量特性とが対応付けられる。長期的に、使用履歴が異なっても、蓄電量特性が連続的かつ一義的に変化することが確認される。
　対極の活物質のＳＯＣ－Ｖがテーブル６３ｂに格納される。
　ＣＰＵ６２はメモリ６３から読み出したプログラムに従って、後述する内部状態推定処理を実行する。
　電圧計測部８は、電圧検知線を介して電池３の両端に夫々接続されており、各電池３の電圧を所定時間間隔で測定する。
　電流計測部９は、電流センサ７を介して電池３に流れる電流を所定時間間隔で計測する。
　電池モジュール１の外部端子５，５は、エンジン始動用のスターターモータ及び電装品等の負荷１１に接続されている。
　ＥＣＵ（Electronic Control Unit）１０は、ＢＭＵ６及び負荷１１に接続されている。
　以下、本実施形態に係る内部状態推定方法について説明する。
　図７は、ＣＰＵ６２によるパラメータ及び蓄電量特性の推定処理の手順を示すフローチャートである。
　ＣＰＵ６２は、所定の間隔でＳ１からの処理を繰り返す。
　ＣＰＵ６２は、前記所定電圧範囲の定電流による実測蓄電量特性を取得する（Ｓ１）。ＣＰＵ６２は、例えば充電又は放電のＳＯＣ－Ｖを取得する。
　ＣＰＵ６２は、実測蓄電量特性から特徴値を抽出する（Ｓ２）。ＣＰＵ６２は、例えばＳＯＣ－Ｖから所定電圧範囲の充電電気量を抽出する。
　ＣＰＵ６２は、取得した特徴値に基づき、電池３の第１蓄電量特性を算出する。（Ｓ３）。
　ＣＰＵ６２は、実測蓄電量特性と第１蓄電量特性とに基づき、最適化計算を行ってパラメータを算出する（Ｓ４）。設計変数は上述のパラメータであり、目的関数は残差二乗和ｆ（ｘ）である。
　ＣＰＵ６２は、パラメータに基づいて再度蓄電量特性を算出し、第２蓄電量特性を得る（Ｓ５）。
　ＣＰＵ６２は、算出したパラメータ及び第２蓄電量特性をテーブル６３ｂに記憶し（Ｓ６）、処理を終了する。
　以下、具体的に説明する。まず、上述のＳ１～Ｓ３について説明する。
（１）ＣＰＵ６２は、電圧範囲Ｃ１Ｖ～Ｃ２Ｖにおける電池３のＳＯＣ－Ｖを取得する（Ｓ１）。ＣＰＵ６２は、該ＳＯＣ－Ｖから特徴値として、電圧範囲Ｃ１Ｖ～Ｃ２Ｖ、即ち正極の電位の範囲がＰ１Ｖ～Ｐ２Ｖであるときの充電電気量Ｑ(x)を抽出する（Ｓ２）。この充電電気量がＱinP1-P2Vと定義される。
（２）テーブル６３ｂには、表１のＮｏ．１～１５に対し、Ｑin P1-P2Vと対応付けて、正極のＶ－ｄＱ／ｄＶが記憶されているとする。また、新品時の負極のＶ－ｄＱ／ｄＶも記憶されているとする。
　例えば取得した特徴値Ｑ(x)in P1-P2Vが、表１のＮｏ．１２及びＮｏ．１５夫々のＱinP1-P2Vの間にある場合、ＣＰＵ６２は、Ｎｏ．１２及びＮｏ．１５夫々のＶ－ｄＱ／ｄＶを取得する。ＣＰＵ６２は、この２つのＶ－ｄＱ／ｄＶを用い、内挿計算を行って、特徴値に対応するＶ－ｄＱ／ｄＶを算出する。内挿計算として天秤法が挙げられる。
　求めるＶ－ｄＱ／ｄＶは下記の式により求められる。
　ｄＱ／ｄＶ＝ｄＱ／ｄＶ(i)×ｑ＋ｄＱ／ｄＶ(i-1)×（１－ｑ）
但し、ｄＱ／ｄＶ(i)：Ｎｏ．ｉのｄＱ／ｄＶ、ここではＮｏ．１５のｄＱ／ｄＶ
　　　ｄＱ／ｄＶ(i-1)：Ｎｏ．i-1のｄＱ／ｄＶ、ここではＮｏ．１２のｄＱ／ｄＶ
　　　ｑ：内挿計算の処理係数、０≦ｑ≦１
　　　　　ｑ＝（Ｑ(i-1)－Ｑ(x)）／（Ｑ(i-1)－Ｑ(i)）
　　　　　Ｑ(i)：Ｎｏ．ｉのＱin P1-P2V、ここではＮｏ．１５のＱin P1-P2V
　　　　　Ｑ(i-1) ：Ｎｏ．i-1のＱin P1-P2V、ここではＮｏ．１２のＱin P1-P2V
　上述した通り、正極の活物質はヒステリシスの性質を有するので、Ｑin P1-P2Vは最大で数十％の変動がある。従って、Ｑ(x)は下記のように補正するのが好ましい。
　図８は、電位範囲Ｐ１－Ｐ２Ｖにおける前記劣化程度ｋ′と充電電気量Ｑin P1-P2Vとの関係を示すグラフである。横軸は充電電気量（％）、縦軸は劣化程度ｋ′である。図８中、実線のグラフは充電時、破線のグラフは放電時のグラフである。
　図８に示すように、Ｑとｋ′との間には直線的な相関が見られる。一方、充電電気量はヒステリシスにより、充電状態によって変動する。図８の場合、最大で２０％変動している。
　Ｑ(x)が充電電気量である場合、補正値Ｑ(x),chaは以下の式により算出される。
　Ｑ(x),cha＝ｎ×（１００＋ΔＱox,max×Ｒcha）/１００×Ｑ(Ｒcha＝０)
　但し、
　ｎ：正極有効度（正極の収縮の度合）、０≦ｎ≦１
　ΔＱox,max：変動量、ΔＱoxの最大値（％）
　ΔＱox＝ΔＳＯＣmax－ΔＳＯＣ
　ΔＳＯＣ：Ｖｕｐ′におけるフル放電ＳＯＣ－Ｖ及びフル充電ＳＯＣ－Ｖ間のＳＯＣの差
　Ｖｕｐ′：電池３のＳＯＣ－Ｖを取得した時点でメモリ６３に記憶されている、後述する電池３の到達最大電圧（上側基準電圧）Ｖｕｐに対応する正極の電位
　ΔＳＯＣmax：ΔＳＯＣの最大値
　Ｒcha：正極の充電ＳＯＣ－Ｖの割合、Ｒcha＝（ΔＱox,max－ΔＱox）/ΔＱox,max＝ΔＳＯＣ/ΔＳＯＣmax、０≦Ｒcha≦１
　Ｑ(Ｒcha＝０)：放電容量
　Ｑ(x)が放電容量である場合、補正値Ｑ(x),disは以下の式により算出される。
　Ｑ(x),dis＝ｎ×Ｑ(x)
　劣化程度ｋ′の推定はＲchaの精度に依存する。前記単極の電位範囲は、充電電気量又は放電容量と単極の平均放電電位との間に直線関係があり、劣化の前後で対極との電位差が変わらない範囲を選択するが、該範囲は対極の状態により変化する可能性がある。対極の状態により電池３のＳＯＣ－Ｖが変化するので、電池３の制御のためには対極の状態も含めた内部状態を推定する必要がある。
　Ｑ(x)及びＲchaから正極のＳＯＣ－Ｖを求め、電池３の第１ＳＯＣ－Ｖを求め、実測ＳＯＣ－Ｖに基づき最適化計算を行うことで、良好に内部状態が推定できる。
（３）ＣＰＵ６２は、Ｑ(x)により求めた正極のＶ－ｄＱ／ｄＶを正極のＳＯＣ－Ｖに変換する。そして、ＣＰＵ６２は、正極のＳＯＣ－Ｖ及び負極のＳＯＣ－Ｖに基づいて電池３のＳＯＣ－Ｖ（第１蓄電量特性）を算出する（Ｓ３）。このときに、パラメータとして、内挿計算の処理係数ｑ1 、正極のＳＯＣ－Ｖ（ｋ1 ）、充電ＳＯＣ－Ｖの割合Ｒcha1 、劣化程度ｋ′1 、正極の収縮の度合ｍ1 、及び容量バランスのずれ量ΔＱim1 が得られている。
（４）上述のＳ４で、ＣＰＵ６２は電池３の現在の内部状態を示すパラメータを算出する。
　ＣＰＵ６２は、下記の式により、目的関数の残差二乗和ｆ（ｘ）を求める。
　ｆ（ｘ）＝ΣC1V  C2V（Ｖ(算出値)－Ｖ(実測値)）2 
　ここで、電圧範囲はＣ１Ｖ～Ｃ２Ｖと異なる範囲に設定してもよい。第１蓄電量特性の推定、パラメータの推定夫々につき、最適な結果が得られる電圧範囲を選択するのが好ましい。Ｓ１で取得した実測ＳＯＣ－Ｖと、Ｓ３で算出した第１ＳＯＣ－Ｖとの誤差が最小になるときのパラメータｋ2 、Ｒcha2 、ｋ′2 、ｍ2 、及びΔＱim2が求められる。
　ＣＰＵ６２は、求められたパラメータを用いて、電池３の充電ＳＯＣ－Ｖ又は放電ＳＯＣ－Ｖ（第２蓄電量特性）を算出する（Ｓ５）。
　図９は、Ｓ４で求めたパラメータを用いて算出した正極の放電ＤＯＤ－Ｖ、電池３の放電ＤＯＤ－Ｖ、及び負極のＤＯＤ－Ｖを示すグラフである。横軸は容量（％）、縦軸は電位又は電圧（Ｖ）である。図中、ｅは正極のＤＯＤ－Ｖ、ｆは負極のＤＯＤ－Ｖ、ｇは電池３のＤＯＤ－Ｖである。
　ｅとｆとのグラフにより示される、正極と負極との容量差である容量バランスのずれ量ΔＱimが、求めたΔＱim2である。
　同様にして、他のパラメータである負極の収縮の度合（ｎ）、及び分極（ΔＶ）、を設計変数として、現在の電池３の内部状態を示すパラメータが求められる。分極（ΔＶ）と実測蓄電量特性を取得した時点の電流とから内部抵抗（Ｒ）を算出することができる。
　図１０は、電池３につき２．０－４．６Ｖの充放電を５０サイクル繰り返した場合に、実測したＳＯＣ－Ｖと、上述の方法によりパラメータを求めて算出した第２ＳＯＣ－Ｖとの差異を求めた結果を示すグラフである。横軸は容量（％）、縦軸は差異（Ｖ）である。
　差異は下記の式により求められる。
　差異＝｛電池３の電圧（算出値）｝－｛電池３の電圧（実測値）｝
　図１０より、差異は小さいことが分かる。
　図１１は、電池３につき２．０－４．６Ｖの充放電を５０サイクル繰り返した場合の、正極のＶ－ｄＱ／ｄＶの実測値と算出値との差異を求めた結果を示すグラフである。横軸は電位（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+）、縦軸は差異（％）である。図中、ｈは充電Ｖ－ｄＱ／ｄＶ、ｉは放電Ｖ－ｄＱ／ｄＶの差異を示す。
　正極のＶ－ｄＱ／ｄＶの実測値は、充放電を５０サイクル繰り返した電池３を解体し、正極のみのＶ－ｄＱ／ｄＶを求めたものである。
　差異は下記の式により求められる。
　差異＝｛ｄＱ／ｄＶ（算出値）－ｄＱ／ｄＶ（実測値）｝/（ｄＱ／ｄＶの絶対値の最大値（実測値））×１００
　図１１より、差異は小さいことが分かる。
　図１２は、電池３につき２．０－４．６Ｖの充放電を５００サイクル繰り返した場合に、実測したＳＯＣ－Ｖと、上述の方法によりパラメータを求めて算出した第２ＳＯＣ－Ｖとの差異を求めた結果を示すグラフである。横軸は容量（％）、縦軸は差異（Ｖ）である。
　差異は下記の式により求められる。
　差異＝｛電池３の電圧（算出値）｝－｛電池３の電圧（実測値）｝
　図１２より、差異は小さいことが分かる。
　図１３は、電池３につき２．０－４．６Ｖの充放電を５００サイクル繰り返した場合の、正極のＶ－ｄＱ／ｄＶの実測値と算出値との差異を求めた結果を示すグラフである。横軸は電位（ＶＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+）、縦軸は差異（％）である。図中、ｊは充電Ｖ－ｄＱ／ｄＶ、ｋは放電Ｖ－ｄＱ／ｄＶの差異を示す。
　正極のＶ－ｄＱ／ｄＶの実測値は、充放電を５０サイクル繰り返した電池３を解体し、正極のみのＶ－ｄＱ／ｄＶを求めたものである。
　差異は下記の式により求められる。
　差異＝｛ｄＱ／ｄＶ（算出値）－ｄＱ／ｄＶ（実測値）｝/（ｄＱ／ｄＶの絶対値の最大値（実測値））×１００
　図１３より、差異は小さいことが分かる。
　以下、直近に記憶した第２ＳＯＣ－Ｖを用いてＳＯＣを推定する場合について説明する。
　図１４及び図１５は、ＣＰＵ６２によるＳＯＣ推定処理の手順を示すフローチャートである。ＣＰＵ６２は、所定の間隔でＳ１１からの処理を繰り返す。以下、第２蓄電量特性がＳＯＣ－ＯＣＶである場合につき説明する。
　予め実験により、ヒステリシスを生じる反応の酸化量及び還元量が小さい電圧が求められ、閾値Ｖ１とされる。電圧がＶ１よりも貴になった後に取得した電圧が上側基準電圧（Ｖｕｐ）に設定される。Ｖｕｐは、取得した電圧が前回取得した電圧より大きい場合に更新される。電圧がＶ１よりも卑になった後に取得した電圧が下側基準電圧（Ｖｌｏｗ）に設定される。Ｖｌｏｗは、取得した電圧が前回取得した電圧より小さい場合に更新される。
　ＣＰＵ６２は、電池３の端子間の電圧及び電流を取得する（Ｓ１１）。閾値Ｖ１及び上側基準電圧ＶｕｐをＯＣＶで示している場合、電池３の電流量が大きいとき、取得した電圧をＯＣＶに補正する必要がある。ＯＣＶへの補正値は、複数の電圧及び電流のデータから回帰直線を用いて、電流がゼロである場合の電圧を推定すること等により得られる。電池３を流れる電流量が暗電流程度に小さい（微小電流である）場合、取得した電圧がＯＣＶとみなされる。
　ＣＰＵ６２は、電流の絶対値が休止閾値以上であるか否かを判定する（Ｓ１２）。休止閾値は、電池３の状態が充電状態又は放電状態と、休止状態とのいずれであるかを判定する為に設定される。ＣＰＵ６２は電流の絶対値が休止閾値以上でないと判定した場合（Ｓ１２：ＮＯ）、処理をＳ２２へ進める。
　ＣＰＵ６２は、電流の絶対値が休止閾値以上であると判定した場合（Ｓ１２：ＹＥＳ）、電流が０より大きいか否かを判定する（Ｓ１３）。電流が０より大きい場合、電池３の状態は充電状態であると判定できる。ＣＰＵ６２は電流が０より大きくないと判定した場合（Ｓ１３：ＮＯ）、処理をＳ１８へ進める。
　ＣＰＵ６２は電流が０より大きいと判定した場合（Ｓ１３：ＹＥＳ）、電圧がＶ１以上であるか否かを判定する（Ｓ１４）。ＣＰＵ６２は電圧がＶ１以上でないと判定した場合（Ｓ１４：ＮＯ）、処理をＳ１７へ進める。
　ＣＰＵ６２は電圧がＶ１以上であると判定した場合（Ｓ１４：ＹＥＳ）、取得した電圧が前回メモリ６３に記憶されたＶｕｐより大きいか否かを判定する（Ｓ１５）。ＣＰＵ６２は電圧が前回のＶｕｐより大きくないと判定した場合（Ｓ１５：ＮＯ）、処理をＳ１７へ進める。
　ＣＰＵ６２は電圧が前回のＶｕｐより大きいと判定した場合（Ｓ１６：ＹＥＳ）、メモリ６３において、電圧をＶｕｐに更新する（Ｓ１６）。
　ＣＰＵ６２は、電流積算によりＳＯを推定し（Ｓ１７）、処理を終了する。
　ＣＰＵ６２は電流が０より小さく、電池３の状態が放電状態であると判定した場合（Ｓ１３：ＮＯ）、電圧がＶ１未満であるか否かを判定する（Ｓ１８）。ＣＰＵ６２は電圧がＶ１未満でないと判定した場合（Ｓ１８：ＮＯ）、処理をＳ２１へ進める。
　ＣＰＵ６２は電圧がＶ１未満であると判定した場合（Ｓ１８：ＹＥＳ）、取得した電圧が前回メモリ６３に記憶された下側基準電圧Ｖｌｏｗより小さいか否かを判定する（Ｓ１９）。
　ＣＰＵ６２は電圧が前回のＶｌｏｗより小さくないと判定した場合（Ｓ１９：ＮＯ）、処理をＳ２１へ進める。
　ＣＰＵ６２は電圧が前回のＶｕｐより小さいと判定した場合（Ｓ１９：ＹＥＳ）、メモリ６３において、電圧をＶｌｏｗに更新する（Ｓ２０）
　ＣＰＵ６２は、電流積算によりＳＯＣを推定し（Ｓ２１）、処理を終了する。
　ＣＰＵ６２は電流の絶対値が休止閾値未満であり、電池３の状態が休止状態であると判定した場合（Ｓ１２：ＮＯ）、設定時間が経過したか否かを判定する（Ｓ２２）。設定時間は、取得した電圧をＯＣＶとみなす為に十分な、実験により求めた時間である。ＣＰＵ６２は、休止状態であると判定してからの電流の取得回数及び取得間隔に基づき、前記時間を超えたか否かを判定する。これにより、休止状態において、より高精度にＳＯＣを推定することができる。
　ＣＰＵ６２は設定時間が経過していないと判定した場合（Ｓ２２：ＮＯ）、電流積算によりＳＯＣを推定し（Ｓ２３）、処理を終了する。
　ＣＰＵ６２は設定時間が経過したと判定した場合（Ｓ２２：ＹＥＳ）、取得した電圧はＯＣＶとみなすことができる。
　ＣＰＵ６２は、直近の第２蓄電量特性をテーブル６３ｂから取得する（Ｓ２４）。ここで、第２蓄電量特性は電池３の充電ＳＯＣ－ＯＣＶ及び放電ＳＯＣ－ＯＣＶである。なお、最後に実測蓄電量特性を取得した日から期間が空いた場合、取得後から現時点までの履歴を考慮して、推定した第２蓄電量特性を補正するのが好ましい。
　ＣＰＵ６２は、取得した第２蓄電量特性に基づいて電圧参照用蓄電量特性を算出する（Ｓ２５）。ＣＰＵ６２は、充電ＳＯＣ－ＯＣＶ、放電ＳＯＣ－ＯＣＶ、及びＶｕｐに基づいて電圧参照用充電ＳＯＣ－ＯＣＶ又は電圧参照用放電ＳＯＣ－ＯＣＶを算出する。ＣＰＵ６２は、例えば、ヒステリシスを生じる反応の酸化量，還元量を考慮して、充電ＳＯＣ－ＯＣＶ、放電ＳＯＣ－ＯＣＶ、及びＶｕｐを用いて、電圧参照用充電ＳＯＣ－ＯＣＶ又は電圧参照用放電ＳＯＣ－ＯＣＶを算出する。Ｖｕｐは現在のＲchaに基づいて補正することができる。
　ＣＰＵ６２は、電圧参照用充電ＳＯＣ－ＯＣＶ又は電圧参照用ＳＯＣ－ＯＣＶにおいて、Ｓ１で取得した電圧に対応するＳＯＣを読み取ってＳＯＣを推定し（Ｓ２６）、処理を終了する。
　なお、ＣＰＵ６２が電圧計測部８から取得する電圧は、電流により多少変動するので、実験により補正係数を求めて電圧を補正することもできる。
　以上のように、本実施形態においては、単極の劣化によるＶＦを反映する複数の蓄電量－電位特性、蓄電量－電圧特性、又はＶ－ｄＱ／ｄＶが記憶され、この蓄電量特性及び特徴値に基づいて第１蓄電量－電圧特性が算出される。第１蓄電量－電圧特性、及び所定電圧範囲で取得した実測蓄電量－電圧特性に基づき、最適化計算により、パラメータが求められる。簡便な方法で、パラメータにより、ＶＦに対応する単極の状態、及び対極の状態が検知できる。従って、電池３の現在の内部状態が高精度に推定できる。
　テーブル６３ｂに記憶する蓄電量特性のデータ数は少なくて済む。
　複雑な使用環境下においても、高精度に内部状態を監視することができる。
　活物質がヒステリシスを有する場合に、現在の単極の劣化状態に応じた第２蓄電量－電位特性、及び蓄電素子の充放電の履歴に基づいて、精度良く電圧参照用蓄電量－電圧特性が推定できる。高容量であり、ＶＦの性質を有する活物質を含む蓄電素子につき、ヒステリシスの挙動の知見を併用することで、蓄電量が比較的容易に推定できる。
　電圧を用いるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量が推定できる。充放電特性に基づいて、ＳＯＣ０％までの放電可能なエネルギー、及びＳＯＣ１００％までに必要な充電エネルギーを予測することができる。現時点での残存電力量と貯蔵可能電力量とが推定できる。
　従って、複数の蓄電素子を用いる場合のバランシング、回生受け入れの制御、蓄電素子を車載した場合の走行距離の推定等を精度良く行うことができる。
　本発明は上述した実施の形態の内容に限定されるものではなく、請求項に示した範囲で種々の変更が可能である。即ち、請求項に示した範囲で適宜変更した技術的手段を組み合わせて得られる実施形態も本発明の技術的範囲に含まれる。
　実施形態１においては、正極がＶＦ及びヒステリシスを有する活物質を含む場合につき説明したが、負極がＶＦ及びヒステリシスを有する活物質を含む場合も、同様にして蓄電量－電位特性又はＶ－ｄＱ／ｄＶを推定することができる。
　本発明に係る電圧参照による蓄電量の推定は休止時に行う場合に限定されず、充電時又は放電時にリアルタイムに行ってもよい。この場合、取得した電圧及び電流から現在のＯＣＶが算出される。ＯＣＶの算出は、複数の電圧及び電流のデータから回帰直線を用いて、電流がゼロである場合の電圧を推定すること等により得られる。また、電流が暗電流のように小さい場合は、取得した電圧ＶをＯＣＶに読み替えることもできる。
　本発明に係る推定装置及び蓄電量推定装置は、車載用のリチウムイオン二次電池に適用される場合に限定されず、鉄道用回生電力貯蔵装置、太陽光発電システム等の他の蓄電装置にも適用できる。微小電流が流れる蓄電装置においては、蓄電素子の正極端子・負極端子間の電圧をＯＣＶとみなすことができる。
　蓄電素子は、リチウムイオン二次電池に限定されるものではなく、ＶＦ及びヒステリシス特性を有する他の二次電池や電気化学セルであってもよい。
　監視装置１００又はＢＭＵ６が蓄電量推定装置である場合を例示した。代替的に、ＣＭＵ（Cell Monitoring Unit）が蓄電量推定装置でよい。蓄電量推定装置は、監視装置１００等が組み込まれた蓄電モジュールの一部であってもよい。蓄電量推定装置は、蓄電素子や蓄電モジュールとは別個に構成されて、蓄熱量推定対象の蓄電素子を含む蓄電モジュールに、蓄熱量の推定時に接続されてもよい。蓄熱量推定装置は、蓄電素子や蓄電モジュールを遠隔監視してもよい。
　本発明は、リチウムイオン二次電池等の蓄電素子の内部状態の推定に適用できる。
　１、５０　電池モジュール（蓄電装置）
　２　ケース
　２１　ケース本体
　２２　蓋部
　２３　ＢＭＵ収容部
　２４　カバー
　２５　中蓋
　２６　仕切り板
　３、２００　電池（蓄電素子）
　３１　ケース
　３２　端子
　３３　電極体
　４　バスバー
　５　外部端子
　６　ＢＭＵ（推定装置）
　６０　情報処理部 
　６２　ＣＰＵ（取得部、第１推定部、第２推定部、第３推定部、第４推定部、第５推定部）
　６３　メモリ（記憶部）
　６３ａ　プログラム
　６３ｂ　テーブル
　７　電流センサ
　８　電圧計測部
　９　電流計測部
　１０　ＥＣＵ
　１００　監視装置（推定装置）
　３００　収容ケース

Claims (12)

  1. 　蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定する推定装置であって、
    　前記蓄電素子の蓄電量－電圧充電特性である第３特性及び／又は蓄電量－電圧放電特性である第４特性の一部の情報を取得する取得部と、
    　充放電の繰り返しにより変化する特徴値の変化に対応して、第１特性、第２特性、第３特性、第４特性、及び、Ｖ－ｄＱ／ｄＶの少なくともいずれかである蓄電量特性を複数記憶し、又は前記特徴値の関数として記憶する記憶部と、
    　前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する第１推定部とを備える、推定装置。
  2. 　前記取得部は、前記蓄電素子の所定電圧範囲における、第３特性及び／又は第４特性の情報を取得し、
    　前記情報から特徴値を求め、該特徴値に基づいて、前記蓄電量特性又は前記関数を参照して、前記情報を取得した時点の第３特性又は第４特性を推定する第２推定部を備え、
    　前記第１推定部は、
    　前記情報と、推定した前記第３特性又は前記第４特性とに基づき、前記蓄電素子の内部状態を示すパラメータを推定する、請求項１に記載の推定装置。
  3. 　前記特徴値は、前記所定電圧範囲における充電電気量、放電容量、及び、平均放電電位の少なくともいずれかである、請求項２に記載の推定装置。
  4. 　前記パラメータは、第１特性、第２特性、第１特性の割合、前記単極の収縮の度合、対極の収縮の度合、前記単極及び前記対極間の容量バランスのずれ量、及び分極から選択される少なくとも１種である、請求項２又は３に記載の推定装置。
    　但し、前記第１特性の割合は下記式で表される。
    　（第１特性の割合）＝（前記情報を取得したときに、前記記憶部に記憶されている前記蓄電素子の到達最大電圧に対応する前記単極の電位における、第１特性及び第２特性間の蓄電量の差）/（前記蓄電量の差の最大値）
  5. 　前記第２推定部は、前記Ｖ－ｄＱ／ｄＶを参照して、内挿計算によりＶ－ｄＱ／ｄＶを推定し、
    　前記パラメータは、前記内挿計算に必要な処理係数を含む、請求項２から４までのいずれか１項に記載の推定装置。
  6. 　前記第１推定部は、最適化計算により、前記情報と、前記第２推定部により推定した前記第３特性又は前記第４特性との差が最小になるように、前記パラメータを求める、請求項２から５までのいずれか１項に記載の推定装置。
  7. 　前記第１推定部により推定した前記パラメータを用いて、第３特性又は第４特性を推定する第３推定部を備える、請求項２から６までのいずれか１項に記載の推定装置。
  8. 　前記活物質は、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、第３特性及び第４特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さく、
    　前記第３推定部により推定した第３特性又は第４特性、及び前記蓄電素子の充放電の履歴に基づいて、取得した電圧により蓄電量を推定するときの参照のための第３特性又は第４特性を推定する第４推定部　を備える、請求項７に記載の推定装置。
  9. 　前記充放電の履歴、参照のための前記第３特性又は前記第４特性、並びに取得した電圧に基づいて、蓄電量を推定する第５推定部
    　を備える、請求項８に記載の推定装置。
  10. 　蓄電素子と、
    　請求項１から９までのいずれか１項に記載の推定装置と
    　を備える、蓄電装置。
  11. 　蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定する推定方法であって、
    　充放電の繰り返しにより変化する特徴値の変化に対応して、第１特性、第２特性、前記蓄電素子の蓄電量－電圧充電特性である第３特性、蓄電量－電圧放電特性である第４特性、及び、Ｖ－ｄＱ／ｄＶ、の少なくともいずれかである蓄電量特性を、記憶部に、複数記憶し、又は前記特徴値の関数として記憶してあり、
    　前記蓄電素子の第３特性及び／又は第４特性の一部の情報を取得し、
    　前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する、推定方法。
  12. 　蓄電量－電位充電特性である第１特性及び蓄電量－電位放電特性である第２特性が充放電の繰り返しにより変化する活物質を含む単極を有する蓄電素子の内部状態を推定するコンピュータに、
    　前記蓄電素子の蓄電量－電圧充電特性である第３特性及び／又は蓄電量－電圧放電特性である第４特性の一部の情報を取得し、
    　第１特性、第２特性、第３特性、第４特性、及び、Ｖ－ｄＱ／ｄＶ、の少なくともいずれかである蓄電量特性が、特徴値の変化に対応して、複数記憶されたテーブルを参照し、又は前記特徴値の関数を参照し、
    　前記情報、及び、前記蓄電量特性又は前記関数に基づいて、前記蓄電素子の内部状態を推定する、
    　処理を実行させる、コンピュータプログラム。
PCT/JP2018/026994 2017-07-19 2018-07-19 推定装置、蓄電装置、推定方法、及びコンピュータプログラム WO2019017411A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP18834355.2A EP3657189A4 (en) 2017-07-19 2018-07-19 ESTIMATION DEVICE, ENERGY STORAGE DEVICE, ESTIMATION METHOD, AND COMPUTER PROGRAM
US16/631,614 US11237214B2 (en) 2017-07-19 2018-07-19 Estimation device, energy storage apparatus, estimation method, and computer program

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2017140257 2017-07-19
JP2017-140257 2017-07-19
JP2018-133649 2018-07-13
JP2018133649A JP6477959B2 (ja) 2017-07-19 2018-07-13 推定装置、蓄電装置、推定方法、及びコンピュータプログラム

Publications (1)

Publication Number Publication Date
WO2019017411A1 true WO2019017411A1 (ja) 2019-01-24

Family

ID=65015172

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/026994 WO2019017411A1 (ja) 2017-07-19 2018-07-19 推定装置、蓄電装置、推定方法、及びコンピュータプログラム

Country Status (2)

Country Link
US (1) US11237214B2 (ja)
WO (1) WO2019017411A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021191993A1 (ja) * 2020-03-24 2021-09-30 Tdk株式会社 二次電池の制御装置、二次電池の制御システム、二次電池パック及び二次電池の制御方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11585862B2 (en) * 2020-02-28 2023-02-21 Denso Corporation Battery deterioration prediction system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002286818A (ja) * 2001-03-26 2002-10-03 Toyota Motor Corp バッテリ容量判定装置
JP2011220917A (ja) 2010-04-13 2011-11-04 Toyota Motor Corp リチウムイオン二次電池の劣化判定装置および劣化判定方法
WO2014045706A1 (ja) * 2012-09-18 2014-03-27 株式会社豊田自動織機 充電状態推定装置及び充電状態推定方法
WO2014156265A1 (ja) * 2013-03-29 2014-10-02 日立オートモティブシステムズ株式会社 電池制御装置
WO2017010475A1 (ja) * 2015-07-13 2017-01-19 三菱電機株式会社 リチウムイオン電池の充電状態推定方法およびリチウムイオン電池の充電状態推定装置
WO2017098686A1 (ja) * 2015-12-10 2017-06-15 ソニー株式会社 電池パック、蓄電装置及び劣化検出方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8531158B2 (en) 2010-11-01 2013-09-10 GM Global Technology Operations LLC Method and apparatus for assessing battery state of health
US9742042B2 (en) 2013-11-23 2017-08-22 Hrl Laboratories, Llc Voltage protection and health monitoring of batteries with reference electrodes
US20170033572A1 (en) * 2015-07-27 2017-02-02 Robert Bosch Gmbh Capacity estimation in a secondary battery
JP6406468B1 (ja) 2017-03-29 2018-10-17 株式会社Ｇｓユアサ 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム
CN110476073B (zh) 2017-03-29 2021-10-29 株式会社杰士汤浅国际 蓄电量估计装置、蓄电模块、蓄电量估计方法、和记录介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002286818A (ja) * 2001-03-26 2002-10-03 Toyota Motor Corp バッテリ容量判定装置
JP2011220917A (ja) 2010-04-13 2011-11-04 Toyota Motor Corp リチウムイオン二次電池の劣化判定装置および劣化判定方法
WO2014045706A1 (ja) * 2012-09-18 2014-03-27 株式会社豊田自動織機 充電状態推定装置及び充電状態推定方法
WO2014156265A1 (ja) * 2013-03-29 2014-10-02 日立オートモティブシステムズ株式会社 電池制御装置
WO2017010475A1 (ja) * 2015-07-13 2017-01-19 三菱電機株式会社 リチウムイオン電池の充電状態推定方法およびリチウムイオン電池の充電状態推定装置
WO2017098686A1 (ja) * 2015-12-10 2017-06-15 ソニー株式会社 電池パック、蓄電装置及び劣化検出方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JOURNAL OF POWER SOURCES, vol. 229, 2013, pages 239 - 248

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021191993A1 (ja) * 2020-03-24 2021-09-30 Tdk株式会社 二次電池の制御装置、二次電池の制御システム、二次電池パック及び二次電池の制御方法

Also Published As

Publication number Publication date
US20200182938A1 (en) 2020-06-11
US11237214B2 (en) 2022-02-01

Similar Documents

Publication Publication Date Title
JP6477957B2 (ja) 推定装置、蓄電装置、推定方法、及びコンピュータプログラム
JP6406468B1 (ja) 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム
JP6409208B1 (ja) 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム
JP6443577B2 (ja) 管理装置、蓄電モジュール、管理方法、及びコンピュータプログラム
WO2015174279A1 (ja) 二次電池監視装置、および、二次電池の電池容量の予測方法
JP7115345B2 (ja) 推定装置、蓄電装置、推定方法、及びコンピュータプログラム
WO2019240225A1 (ja) 推定装置、蓄電装置、推定方法、及びコンピュータプログラム
WO2018181620A1 (ja) 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム
WO2019017411A1 (ja) 推定装置、蓄電装置、推定方法、及びコンピュータプログラム
WO2020241448A1 (ja) 推定装置、推定方法及びコンピュータプログラム
JP2022101196A (ja) 劣化推定装置、劣化推定方法及びコンピュータプログラム
WO2019017183A1 (ja) 推定装置、蓄電装置、推定方法、及びコンピュータプログラム
JP2021058017A (ja) 推定装置、推定方法及びコンピュータプログラム
JP6406469B1 (ja) 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム
WO2018181624A1 (ja) 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム
JP6406470B1 (ja) 管理装置、蓄電モジュール、管理方法、及びコンピュータプログラム
WO2018221514A1 (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: 18834355

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2018834355

Country of ref document: EP

Effective date: 20200219