JP5036662B2 - Secondary battery monitoring device and secondary battery system - Google Patents

Secondary battery monitoring device and secondary battery system Download PDF

Info

Publication number
JP5036662B2
JP5036662B2 JP2008225649A JP2008225649A JP5036662B2 JP 5036662 B2 JP5036662 B2 JP 5036662B2 JP 2008225649 A JP2008225649 A JP 2008225649A JP 2008225649 A JP2008225649 A JP 2008225649A JP 5036662 B2 JP5036662 B2 JP 5036662B2
Authority
JP
Japan
Prior art keywords
battery
value
resistance
estimated
unit
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.)
Expired - Fee Related
Application number
JP2008225649A
Other languages
Japanese (ja)
Other versions
JP2010060406A (en
Inventor
也寸志 天野
修二 戸村
孝治 梅野
哲矢 淵本
毅 竹本
伸烈 芳賀
勇二 西
Original Assignee
トヨタ自動車株式会社
株式会社豊田中央研究所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by トヨタ自動車株式会社, 株式会社豊田中央研究所 filed Critical トヨタ自動車株式会社
Priority to JP2008225649A priority Critical patent/JP5036662B2/en
Publication of JP2010060406A publication Critical patent/JP2010060406A/en
Application granted granted Critical
Publication of JP5036662B2 publication Critical patent/JP5036662B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/70Energy storage systems for electromobility, e.g. batteries

Description

本発明は、二次電池の監視装置に係り、特に大電流域で発生する電池抵抗の上昇度合いを推定する二次電池監視装置に関する。   The present invention relates to a secondary battery monitoring apparatus, and more particularly to a secondary battery monitoring apparatus that estimates the degree of increase in battery resistance that occurs in a large current region.

充放電可能な二次電池がハイブリッド電気車両(HEV)に搭載される場合、例えば、二次電池は電動機を車両駆動力源として駆動する一方で、この電動機の回生発電時の発電電力やエンジンの回転に伴って発生する発電機の発電電力等によって、二次電池が充電される。二次電池の状態は使用状況によって種々に変化するため、従来より二次電池の状態の監視、さらには監視結果に基づいて電池性能の劣化を防止することが行われている。   When a chargeable / dischargeable secondary battery is mounted on a hybrid electric vehicle (HEV), for example, the secondary battery drives an electric motor as a vehicle driving force source, while the electric power generated during regenerative power generation of the electric motor or the engine The secondary battery is charged by the power generated by the generator generated by the rotation. Since the state of the secondary battery changes variously depending on the use state, conventionally, the state of the secondary battery is monitored, and further, deterioration of the battery performance is prevented based on the monitoring result.

例えば特開平11−187577号公報(特許文献1)には、電池の使用環境および電池の状態に応じた適切な電力で充放電を行なうために、電池温度が所定温度以下または所定温度以上のときには、充電電力上限値および放電電力上限値を常温時に比して小さくする充放電制御装置が開示されている。この装置によれば、過充放電による二次電池の劣化を防止できると述べられている。   For example, in Japanese Patent Laid-Open No. 11-187777 (Patent Document 1), in order to perform charging / discharging with an appropriate electric power according to the use environment of the battery and the state of the battery, when the battery temperature is equal to or lower than a predetermined temperature. A charging / discharging control device is disclosed in which the charging power upper limit value and the discharging power upper limit value are made smaller than those at normal temperature. According to this device, it is stated that deterioration of the secondary battery due to overcharge / discharge can be prevented.

また、特開2004−31170号公報(特許文献2)には、分極の影響度合いと、電流と電圧との関係の特性とを考慮して、二次電池の内部抵抗を高精度に算出することができる二次電池の内部抵抗検出装置が開示されている。この装置では内部抵抗検出に、分極状態(分極の影響度合い)を表す指数として分極指数が用いられる。当該分極指数は、電極近傍の溶液濃度を電気量で表現したものである。当該分極指数によって、充放電による電極近傍の溶液濃度変化および拡散による解消分を考慮して内部抵抗を検出する。   Japanese Patent Application Laid-Open No. 2004-31170 (Patent Document 2) calculates the internal resistance of a secondary battery with high accuracy in consideration of the degree of influence of polarization and the characteristics of the relationship between current and voltage. An internal resistance detection device for a secondary battery that can be used is disclosed. In this apparatus, a polarization index is used as an index representing the polarization state (degree of influence of polarization) for detecting internal resistance. The polarization index represents the solution concentration in the vicinity of the electrode in terms of electricity. Based on the polarization index, the internal resistance is detected in consideration of the solution concentration change near the electrode due to charge and discharge and the elimination due to diffusion.

さらに、特許文献2には、内部抵抗を高精度に検出し、内部抵抗の増加によって二次電池の劣化の検出を行なうことにより、エンジンの始動性の低下等を未然かつ確実に防止可能であると述べられている。   Furthermore, in Patent Document 2, it is possible to prevent deterioration of engine startability and the like in advance by reliably detecting internal resistance and detecting deterioration of the secondary battery by increasing the internal resistance. It is stated.

また、特開2000−123886号公報(特許文献3)では、履歴現象による影響とは係わりなく、車両用二次電池の満充電状態を判定することが可能な満充電判定装置が開示されている。特に、バッテリ内に生じる分極による電圧変化と、分極の大きさを示す分極指標との関係を示す一次直線の勾配を算出するとともに、算出した勾配と、予め定められた満充電状態の電圧変化と分極指標との関係を示す一次直線の勾配との比較判定によって、満充電状態を判定することが開示されている。また、分極による電圧変化と電極表面での電解液の濃度変化とが関連することが記載されている。   Japanese Patent Laid-Open No. 2000-123886 (Patent Document 3) discloses a full charge determination device capable of determining the full charge state of a vehicular secondary battery irrespective of the influence of the history phenomenon. . In particular, the gradient of the linear line indicating the relationship between the voltage change caused by the polarization generated in the battery and the polarization index indicating the magnitude of the polarization is calculated, and the calculated gradient and the voltage change in a fully charged state are determined in advance. It is disclosed that the fully charged state is determined by comparison with a gradient of a linear straight line indicating a relationship with the polarization index. Further, it is described that a voltage change due to polarization is related to a change in the concentration of the electrolytic solution on the electrode surface.

さらに、特開2006−42497号公報(特許文献4)には、二次電池内部でのイオン濃度分布推定に基づくローカルSOCに応じて、一部の動力源が二次電池からの電力によって駆動力を発生する動力出力装置において、複数の動力源の間での駆動力配分制御を適正化することが開示されている。具体的には、二次電池の電極でのイオン濃度分布推定に基づくローカルSOCと、二次電池全体マクロで見た全体SOCとの比較に従って、駆動力配分制御が修正される。   Furthermore, in Japanese Patent Application Laid-Open No. 2006-42497 (Patent Document 4), according to the local SOC based on the estimation of the ion concentration distribution inside the secondary battery, some power sources are driven by the power from the secondary battery. In the power output apparatus that generates the power, it is disclosed to optimize the driving force distribution control among a plurality of power sources. Specifically, the driving force distribution control is corrected according to the comparison between the local SOC based on the estimation of the ion concentration distribution at the electrode of the secondary battery and the overall SOC as seen in the secondary battery overall macro.

特開平11−187577号公報JP-A-11-187777 特開2004−31170号公報JP 2004-31170 A 特開2000−123886号公報JP 2000-123886 A 特開2006−42497号公報JP 2006-42497 A

一般に、二次電池の充放電を大電流で繰り返すと、そのような大電流域において電池抵抗(内部抵抗)が上昇する場合がある。電池抵抗の上昇は電池性能の低下を招くので、当該抵抗上昇現象の発生度合い、換言すれば電池抵抗の上昇度合いを把握することが望まれる。しかしながら、上記特許文献1〜4にはかかる電池抵抗上昇度合いの検出について述べられていない。   In general, when charging / discharging of a secondary battery is repeated with a large current, the battery resistance (internal resistance) may increase in such a large current region. Since an increase in battery resistance causes a decrease in battery performance, it is desirable to grasp the degree of occurrence of the resistance increase phenomenon, in other words, the degree of increase in battery resistance. However, the above Patent Documents 1 to 4 do not describe the detection of the battery resistance increase degree.

ここで、上記の大電流域における電池抵抗の上昇は、非可逆的な経年劣化による抵抗上昇とは異なり、可逆的であることが確認された。すなわち、当該電池抵抗上昇現象は、放電過多が原因の場合は充電によって解消可能であり、逆に充電過多が原因の場合は放電によって解消可能であることが、実験から確認された。このため、大電流域での電池抵抗上昇度合いを取得するためには、経年劣化による抵抗上昇とは異なる新たな手法が必要である。   Here, it was confirmed that the battery resistance increase in the above-described large current region is reversible, unlike the resistance increase due to irreversible aging. That is, it has been confirmed from experiments that the battery resistance increase phenomenon can be eliminated by charging when it is caused by excessive discharge, and can be eliminated by discharging when it is caused by excessive charging. For this reason, in order to acquire the degree of battery resistance increase in a large current region, a new method different from the resistance increase due to aging is necessary.

本発明は、大電流域における電池抵抗上昇の度合いを推定可能な二次電池監視装置及び二次電池システムを提供することを目的とする。 An object of the present invention is to provide a secondary battery monitoring device and a secondary battery system capable of estimating the degree of battery resistance increase in a large current region.

本発明に係る二次電池の監視装置は、二次電池の電池電圧と電池電流を電池状態量として測定する測定部と、バトラー・ボルマーの式において、負極中および正極中における電気化学反応を所定のものと仮定した第1モデル式、活物質拡散モデル式、及び前記電池状態量の一方に基づき、推定電池状態量を取得する取得部と、前記推定電池状態量と前記電池状態量との間の電流推定誤差、及び電極間電解液濃度差に基づいて前記二次電池の電池抵抗の上昇度合いを推定する電池抵抗上昇度合い推定部と、を備える。
また、本発明に係る二次電池の監視装置において、前記取得部は、前記第1モデル式において、前記負極と前記正極の活物質モデルを別々とせずにまとめて単独の素子として簡略化した第2モデル式に基づいて、推定電池状態量を取得する取得部である。 Further, in the secondary battery monitoring device according to the present invention, the acquisition unit simplifies the active material model of the negative electrode and the positive electrode as a single element in the first model formula without separating them. This is an acquisition unit that acquires the estimated battery state based on the two model formulas. Monitoring device for a secondary battery according to the present invention includes a measuring unit for measuring a battery voltage and battery current of the secondary battery as batteries state quantity, in the formula of the Butler-Volmer, the electrochemical reaction at the negative electrode and the positive electrode Based on one of the first model formula, the active material diffusion model formula, and the battery state quantity assumed to be predetermined, an acquisition unit that obtains an estimated battery state quantity, and the estimated battery state quantity and the battery state quantity current estimation error between, and a battery resistance increasing degree estimation unit for estimating the degree of increase battery resistance of the secondary battery based on the inter-electrode electrolytic solution concentration differences, Ru comprising a. Monitoring device for a secondary battery according to the present invention includes a measuring unit for measuring a battery voltage and battery current of the secondary battery as batteries state quantity, in the formula of the Butler-Volmer, the electrical reaction at the negative electrode and the positive electrode Based on one of the first model formula, the active material diffusion model formula, and the battery state quantity assumed to be predetermined, an acquisition unit that obtains an estimated battery state quantity, and the estimated battery state quantity and the battery state quantity current estimation error between, and a battery resistance increasing degree estimation unit for estimating the degree of increase battery resistance of the secondary battery based on the inter-electrode electrolytic solution concentration differences, Ru comprising a.
Further, in the monitoring apparatus for a secondary battery according to the present invention, the acquisition unit may be simplified as a single element by combining the active material models of the negative electrode and the positive electrode in the first model formula without separating them. It is an acquisition part which acquires an estimated battery state quantity based on 2 model formulas. Further, in the monitoring apparatus for a secondary battery according to the present invention, the acquisition unit may be simplified as a single element by combining the active material models of the negative electrode and the positive electrode in the first model formula without separating them. It is an acquisition part which acquires an estimated battery state quantity based on 2 model formulas.

ここで、前記電池抵抗上昇度合い推定部は、前記一方の電池状態量についての推定値と測定値とともに、前記二次電池について推定された電極間の電解液塩濃度差と、前記電池電流の測定値の移動平均値との少なくとも一方に基づいて、前記電池抵抗上昇度合いを推定することが好ましい。 Here, the battery resistance increase degree estimation unit measures the difference between the electrolyte salt concentration estimated for the secondary battery and the battery current together with the estimated value and the measured value for the one battery state quantity. The battery resistance increase degree is preferably estimated based on at least one of the moving average values.

また、前記二次電池の電池温度を測定する温度測定部をさらに備え、前記電池状態量推定部は、前記他方の電池状態量の測定値とともに前記電池温度の測定値に基づき、前記一方の電池状態量を前記電池モデルに従って推定することが好ましい。   The battery further includes a temperature measurement unit that measures a battery temperature of the secondary battery, wherein the battery state quantity estimation unit is configured based on the measured value of the battery temperature together with the measured value of the other battery state quantity. It is preferable to estimate the state quantity according to the battery model.

また、前記電池抵抗上昇度合い推定部は、前記電池抵抗上昇度合いとして、前記一方の電池状態量についての測定値と推定値との比もしくは差、前記比もしくは差の増加率であって前記二次電池について推定された電極間の電解液塩濃度差に対する増加率、または、前記比もしくは差と前記推定された電解液塩濃度差との一次近似式において前記推定電解液塩濃度差を0とした場合の値を算出することが好ましい。 In addition, the battery resistance increase degree estimation unit may calculate, as the battery resistance increase degree, a ratio or difference between a measured value and an estimated value for the one battery state quantity , an increase rate of the ratio or difference, The rate of increase with respect to the electrolyte salt concentration difference between the electrodes estimated for the battery, or the estimated electrolyte salt concentration difference in the linear approximation of the ratio or difference and the estimated electrolyte salt concentration difference is set to 0 It is preferable to calculate the case value.

また、前記電池抵抗上昇度合いが予め設定された通知設定値よりも大きい場合に、前記電池抵抗の上昇を通知する通知部をさらに備えることが好ましい。   In addition, it is preferable that the information processing apparatus further includes a notification unit that notifies the battery resistance increase when the battery resistance increase degree is greater than a preset notification setting value.

また、前記電池抵抗上昇度合いが予め設定された使用制限設定値よりも大きい場合に、充電電力量、放電電力量、充電電流量、放電電流量、電池電圧上限値、電池電圧下限値の少なくともいずれか1を制限する使用制限部をさらに備えることが好ましい。 Further, when the battery resistance increasing degree is greater than a preset use limitation set value, the amount of charge power, amount of discharge power, the amount of charge current, discharge current, battery voltage upper limit, at least for batteries voltage lower limit It is preferable to further include a use restriction unit that restricts any one .

また、前記電池抵抗上昇度合いが予め設定された充電率制御中心調整設定値よりも大きい場合に、充電率制御の制御中心値を調整する充電率制御部をさらに備えることが好ましい。   Moreover, it is preferable to further include a charge rate control unit that adjusts the control center value of the charge rate control when the battery resistance increase degree is larger than a preset charge rate control center adjustment set value.

また、前記二次電池の診断を行う診断部をさらに備え、前記診断部は、前記電池抵抗上昇度合いが予め設定された電圧下限値診断禁止設定値よりも大きい場合に、前記電池電圧の下限値についての故障診断を行わないことが好ましい。   In addition, a diagnostic unit that diagnoses the secondary battery is further provided, and the diagnostic unit has a lower limit value of the battery voltage when the battery resistance increase degree is greater than a preset voltage lower limit value diagnosis prohibition set value. It is preferable not to perform a fault diagnosis for.

また、前記電池抵抗上昇度合いが予め設定された充放電要求設定値よりも大きい場合に、前記二次電池の強制的な充電または放電を行う強制的充放電制御部をさらに備えることが好ましい。   In addition, it is preferable to further include a forcible charge / discharge control unit that forcibly charges or discharges the secondary battery when the battery resistance increase degree is larger than a preset charge / discharge request set value.

また、前記二次電池および前記監視装置がエンジン搭載車両に適用される場合、前記強制的充放電制御部は、エンジンに対して要求するエンジン要求パワーを増加させて、エンジン出力によって前記強制的な充電または放電を行うことが好ましい。   In addition, when the secondary battery and the monitoring device are applied to an engine-equipped vehicle, the forced charge / discharge control unit increases the required engine power required for the engine, and the forced output according to the engine output. It is preferable to charge or discharge.

また、前記二次電池が複数個接続されて電池パックが構成されており、前記強制的充放電制御部は、前記電池パックの温度が予め定められた定常状態にない場合には、前記強制的な充電または放電を行わないことが好ましい。   In addition, a plurality of the secondary batteries are connected to form a battery pack, and the forced charge / discharge control unit performs the forced charge / discharge control when the temperature of the battery pack is not in a predetermined steady state. It is preferable not to perform proper charging or discharging.

また、前記電池モデルにおける抵抗成分のパラメータについて前記二次電池の経時劣化に伴う変化を推定し、当該経時劣化推定結果を前記電池モデルに反映させるパラメータ変化推定部をさらに備えることが好ましい。   Further, it is preferable to further include a parameter change estimation unit that estimates a change of the resistance component in the battery model due to deterioration of the secondary battery with time and reflects the result of deterioration with time in the battery model.

また、前記電池抵抗上昇度合い推定部は、前記パラメータ変化推定部による前記経時劣化推定結果に基づいて前記電池抵抗上昇度合いの推定を実行するか否かを判定することが好ましい。   Moreover, it is preferable that the said battery resistance increase degree estimation part determines whether the estimation of the said battery resistance increase degree is performed based on the said temporal deterioration estimation result by the said parameter change estimation part.

また、前記電池状態量推定部は、前記電池電圧の測定値に基づき前記電池モデルに従って前記電池電流を推定し、前記電池抵抗上昇度合い推定部は、前記電池電流の推定値と測定値とに基づいて前記電池抵抗上昇度合いを推定することが好ましい。   The battery state quantity estimating unit estimates the battery current according to the battery model based on the measured value of the battery voltage, and the battery resistance increase degree estimating unit is based on the estimated value and measured value of the battery current. It is preferable to estimate the degree of battery resistance increase.

また、前記電池状態量推定部は、前記電池電流の測定値に基づき前記電池モデルに従って前記電池電圧を推定し、前記電池抵抗上昇度合い推定部は、前記電池電圧の推定値と測定値とに基づいて前記電池抵抗上昇度合いを推定することが好ましい。
また、本発明に係る二次電池システムは、二次電池の電池電圧と電池電流を電池状態量として測定する測定部と、バトラー・ボルマーの式において、負極中および正極中における電気化学反応を所定のものと仮定した第1モデル式、活物質拡散モデル式、及び前記電池状態量の一方に基づき、推定電池状態量を取得する取得部と、前記推定電池状態量が前記電池状態量よりも所定値以上大きい場合に所定の処理を行う処理部と、を備える。 Further, in the secondary battery system according to the present invention, a measuring unit that measures the battery voltage and the battery current of the secondary battery as the battery state amount, and the Butler-Bolmer formula, determine the electrochemical reaction in the negative electrode and the positive electrode. An acquisition unit that acquires an estimated battery state amount based on one of the first model formula, the active material diffusion model formula, and the battery state amount, and the estimated battery state amount is determined more than the battery state amount. A processing unit that performs a predetermined process when the value is larger than the value is provided.
また、前記処理部は、電極間電解液濃度差が所定値よりも大きい場合に前記所定の処理を行うことが好ましい。 Further, it is preferable that the processing unit performs the predetermined processing when the difference in the concentration of the electrolytic solution between the electrodes is larger than the predetermined value. The battery state quantity estimating unit estimates the battery voltage according to the battery model based on the measured value of the battery current, and the battery resistance increase degree estimating unit is based on the estimated value and measured value of the battery voltage. It is preferable to estimate the degree of battery resistance increase. The battery state quantity estimating unit estimates the battery voltage according to the battery model based on the measured value of the battery current, and the battery resistance increase degree estimating unit is based on the estimated value and measured value of the battery voltage. It is preferred. to estimate the degree of battery resistance increase.
In addition, the secondary battery system according to the present invention includes a measurement unit that measures the battery voltage and battery current of the secondary battery as battery state quantities, and the electrochemical reaction in the negative electrode and the positive electrode in the Butler-Volmer formula. An acquisition unit that acquires an estimated battery state quantity based on one of the first model formula, the active material diffusion model formula, and the battery state quantity, and the estimated battery state quantity is more predetermined than the battery state quantity And a processing unit that performs a predetermined process when the value is larger than the value. In addition, the secondary battery system according to the present invention includes a measurement unit that measures the battery voltage and battery current of the secondary battery as battery state quantities, and the electrical reaction in the negative electrode and the positive electrode in the Butler-Volmer formula. An acquisition unit that acquires an estimated battery state quantity based on one of the first model formula, the active material diffusion model formula, and the battery state quantity, and the estimated battery state quantity is more predetermined than the battery state quantity And a processing unit that performs a predetermined process when the value is larger than the value.
Moreover, it is preferable that the said process part performs the said predetermined process, when the electrolyte solution concentration difference between electrodes is larger than a predetermined value. Moreover, it is preferred that the said process part performs the said predetermined process, when the electrolyte solution concentration difference between electrodes is larger than a predetermined value.

上記構成によれば、大電流域での電池抵抗上昇の度合いを推定することができる。 According to the above configuration, the degree of battery resistance increase in a large current region can be estimated.

まず、図1に充放電可能な二次電池10の概略構成図を示す。以下の説明では二次電池10としてリチウムイオン電池を例示するが、二次電池10は当該例示に限定されるものではない。 First, the schematic block diagram of the secondary battery 10 which can be charged / discharged in FIG. 1 is shown. In the following description, a lithium ion battery is exemplified as the secondary battery 10, but the secondary battery 10 is not limited to the illustration.

図1に示すように、二次電池10は、負極12と、セパレータ14と、正極15とを含む。電極12,15はリチウムイオン電池に用いられる各種材料で構成可能であり、また、セパレータ14は例えば負極12と正極15との間に設けられた樹脂に電解液を浸透させることによって構成可能である。セパレータ14は以下の説明から理解されるように「イオン伝導体」と呼ぶことが可能である。なお、図1では座標軸xを電極12,15の厚み方向、換言すれば電極12,15の積層方向に規定している。   As shown in FIG. 1, the secondary battery 10 includes a negative electrode 12, a separator 14, and a positive electrode 15. The electrodes 12 and 15 can be made of various materials used for lithium ion batteries, and the separator 14 can be made by, for example, infiltrating an electrolyte into a resin provided between the negative electrode 12 and the positive electrode 15. . The separator 14 can be referred to as an “ionic conductor” as will be understood from the following description. In FIG. 1, the coordinate axis x is defined in the thickness direction of the electrodes 12 and 15, in other words, in the stacking direction of the electrodes 12 and 15.

負極12および正極15の各々は、球状の活物質18の集合体で構成される。放電時には、負極12の活物質18の界面上では、リチウムイオンLi+および電子e-を放出する化学反応が行なわれる。一方、正極15の活物質18の界面上では、リチウムイオンLi+および電子e-を吸収する化学反応が行なわれる。 Each of the negative electrode 12 and the positive electrode 15 is composed of an aggregate of spherical active materials 18. During discharge, a chemical reaction that releases lithium ions Li + and electrons e is performed on the interface of the active material 18 of the negative electrode 12. On the other hand, a chemical reaction that absorbs lithium ions Li + and electrons e is performed on the interface of the active material 18 of the positive electrode 15.

負極12には電子e-を吸収する集電体13が設けられ、正極15には電子e-を放出する集電体16が設けられる。なお、集電体は電流コレクタとも呼ばれる。負極の集電体13は例えば銅で構成され、正極の集電体16は例えばアルミニウムで構成される。集電体13には負極端子11nが設けられ、集電体16には正極端子11pが設けられる。セパレータ14を介したリチウムイオンLi+の授受によって、二次電池10では充放電が行われ、充電電流Ib(<0)または放電電流Ib(>0)が生じる。 The negative electrode 12 is provided with a current collector 13 that absorbs electrons e , and the positive electrode 15 is provided with a current collector 16 that emits electrons e . The current collector is also called a current collector. The negative electrode current collector 13 is made of, for example, copper, and the positive electrode current collector 16 is made of, for example, aluminum. The current collector 13 is provided with a negative electrode terminal 11n, and the current collector 16 is provided with a positive electrode terminal 11p. The secondary battery 10 is charged and discharged by the exchange of lithium ions Li + through the separator 14, and a charging current Ib (<0) or a discharging current Ib (> 0) is generated.

上記のように、放電時には、負極12から放出されたリチウムイオンLi+は、拡散および泳動によって正極15へ移動し、正極15に吸収される。このとき、電解液内におけるリチウムイオンLi+の拡散に遅れが生じると、負極12内の電解液中ではリチウムイオンLi+濃度(すなわち電解液塩濃度)が増加する一方、正極15内の電解液中ではリチウムイオンLi+濃度が減少する。かかる様子を図2に示す。なお、図2中に示した平均塩濃度とは、例えば二次電池10の長時間の放置によって電解液塩濃度が二次電池10内全体において均一になった場合の当該電解液塩濃度に対応する。 As described above, at the time of discharging, lithium ions Li + released from the negative electrode 12 move to the positive electrode 15 by diffusion and migration and are absorbed by the positive electrode 15. At this time, if a delay occurs in the diffusion of lithium ions Li + in the electrolytic solution, the lithium ion Li + concentration (that is, the electrolyte salt concentration) increases in the electrolytic solution in the negative electrode 12, while the electrolytic solution in the positive electrode 15. Among them, the lithium ion Li + concentration decreases. Such a situation is shown in FIG. The average salt concentration shown in FIG. 2 corresponds to the electrolyte salt concentration when the electrolyte salt concentration becomes uniform throughout the secondary battery 10 due to, for example, leaving the secondary battery 10 for a long time. To do.

ここで、図3に、電解液塩濃度と反応抵抗との関係を示す。なお、反応抵抗は、活物質18の界面での反応電流発生時に等価的に電気抵抗として作用する抵抗であり、換言すれば電極表面におけるリチウムイオンLi+の出入りに関する抵抗成分である。反応抵抗は電荷移動抵抗とも呼ばれる。 Here, FIG. 3 shows the relationship between the electrolyte salt concentration and the reaction resistance. The reaction resistance is a resistance that acts equivalently as an electric resistance when a reaction current is generated at the interface of the active material 18, in other words, a resistance component related to the entry and exit of lithium ions Li + on the electrode surface. Reaction resistance is also called charge transfer resistance.

図3の特性図によれば、反応抵抗は電解液塩濃度の関数であることが分かる。特に、あるしきい値塩濃度cthを境界にして、電解液塩濃度が高い領域では電解液塩濃度の変化に対して反応抵抗の変化は小さく、他方、電解液塩濃度が低い領域では電解液塩濃度が低いほど反応抵抗値が大きくなる。すなわち、電解液塩濃度が低い領域では、電解液塩濃度が高い領域と比較して、電解液塩濃度に対する反応抵抗値の変化率が大きい。 According to the characteristic diagram of FIG. 3, it can be seen that the reaction resistance is a function of the electrolyte salt concentration. In particular, with a certain threshold salt concentration c th as a boundary, in the region where the electrolyte salt concentration is high, the change in the reaction resistance is small with respect to the change in the electrolyte salt concentration, while in the region where the electrolyte salt concentration is low The reaction resistance value increases as the liquid salt concentration decreases. That is, in the region where the electrolyte salt concentration is low, the rate of change of the reaction resistance value with respect to the electrolyte salt concentration is larger than in the region where the electrolyte salt concentration is high.

このため、図2および図3に鑑みれば、放電時に正極15内での電解液塩濃度が減少した場合であっても、正極15内の電解液塩濃度が高いときには、反応抵抗の低下はほとんど生じないことが分かる。これに対し、正極15内の電解液塩濃度が低いときには、放電時における正極15内での電解液塩濃度の低下は、反応抵抗の増加を招くことが分かる。   Therefore, in view of FIG. 2 and FIG. 3, even when the electrolyte salt concentration in the positive electrode 15 decreases during discharge, when the electrolyte salt concentration in the positive electrode 15 is high, there is almost no decrease in reaction resistance. It turns out that it does not occur. On the other hand, when the electrolyte salt concentration in the positive electrode 15 is low, it can be seen that a decrease in the electrolyte salt concentration in the positive electrode 15 during discharge causes an increase in reaction resistance.

このような反応抵抗増加の要因として、例えば平均電解液塩濃度が減少することによって正極内電解液塩濃度がしきい値塩濃度cthを下回る場合(図4(a)参照)や、例えば放電が繰り返されて累積的に正極内電解液塩濃度が低下することによって正極内電解液塩濃度がしきい値塩濃度cthを下回る場合(図4(b)参照)が考えられる。 As a cause of such an increase in reaction resistance, for example, when the average electrolyte salt concentration decreases to cause the electrolyte salt concentration in the positive electrode to fall below the threshold salt concentration c th (see FIG. 4A), for example, discharge Is repeated and the electrolyte salt concentration in the positive electrode is cumulatively decreased, so that the electrolyte salt concentration in the positive electrode is lower than the threshold salt concentration c th (see FIG. 4B).

なお、上記では反応抵抗の上昇が、放電時に正極15内の電解液塩濃度が低下することによって発現する場合を例示したが、充電時にも負極12内の電解液塩濃度が低下することによって同様に発現する。   In addition, although the case where the increase in reaction resistance is expressed by the decrease in the electrolyte salt concentration in the positive electrode 15 at the time of discharging has been exemplified above, the same phenomenon is caused by the decrease in the electrolyte salt concentration in the negative electrode 12 at the time of charging. Expressed in

反応抵抗と、電極12,15での電子e-の移動に対する純電気的な抵抗(純抵抗)とを併せたものが、二次電池10をマクロに見た場合の電池抵抗(内部抵抗)における直流抵抗成分に相当する。このため、本実施の形態では、上述の反応抵抗が上昇する現象を電池抵抗の上昇現象として検出して監視するための装置を説明する。 The combination of the reaction resistance and the pure electrical resistance (pure resistance) against the movement of electrons e at the electrodes 12 and 15 is the battery resistance (internal resistance) when the secondary battery 10 is viewed macroscopically. Corresponds to DC resistance component. For this reason, in the present embodiment, an apparatus for detecting and monitoring the phenomenon in which the reaction resistance increases as a battery resistance increase phenomenon will be described.

図5に、実施の形態に係る二次電池の監視装置5を説明する概略構成図を示す。なお、図5には、説明のために、二次電池10および接続対象20も併記している。 In FIG. 5, the schematic block diagram explaining the monitoring apparatus 5 of the secondary battery which concerns on embodiment is shown. In FIG. 5, the secondary battery 10 and the connection target 20 are also shown for explanation.

なお、図5では1個の二次電池10を例示しているが、複数個の二次電池10が直列接続された電池パックが利用される場合もある。 Although FIG. 5 illustrates one secondary battery 10, a battery pack in which a plurality of secondary batteries 10 are connected in series may be used.

また、二次電池10が接続される接続対象20は、例えば二次電池10からの出力電力によって駆動される負荷であり、また、発電・給電要素を含むことも可能である。一例として二次電池10がハイブリッド電気車両(Hybrid Electrical Vehicle:HEV)に適用される場合、接続対象20には例えばモータ・ジェネレータ(M/G)と、インバータと、エンジンとが含まれる。M/Gは、二次電池10から電力が供給されるときはモータとして機能し(すなわち負荷となる)、例えば車両の制動時(減速時)やエンジンによって駆動されるときには発電機として機能する回転電機である。M/Gと二次電池10との間での電力の授受はインバータを介して行われる。   Further, the connection target 20 to which the secondary battery 10 is connected is, for example, a load driven by output power from the secondary battery 10, and may include a power generation / feeding element. As an example, when the secondary battery 10 is applied to a hybrid electric vehicle (HEV), the connection target 20 includes, for example, a motor / generator (M / G), an inverter, and an engine. M / G functions as a motor (i.e., becomes a load) when power is supplied from the secondary battery 10, for example, rotation that functions as a generator when the vehicle is braked (during deceleration) or driven by an engine. Electric. Power is exchanged between the M / G and the secondary battery 10 via an inverter.

図5に例示の監視装置5は、電流測定部30と、電圧測定部32と、温度測定部34と、監視制御部50とを含んでいる。 The monitoring device 5 illustrated in FIG. 5 includes a current measurement unit 30, a voltage measurement unit 32, a temperature measurement unit 34, and a monitoring control unit 50.

各測定部30,32,34は、例えば一般的なセンサ、測定器等で構成することが可能である。 Each of the measurement units 30, 32, and 34 can be configured by, for example, a general sensor, a measuring instrument, or the like.

電流測定部30は、二次電池10と接続対象20との間に接続されており、電池電流Ibを測定する。電池電流Ibの向き、換言すればその電流値Ibの符号(正負)から、当該電流Ibが二次電池10から流出する電流(放電電流)であるか、二次電池10へ流入する電流(充電電流)であるかを判別することが可能である。電池電流Ibの符号と放電電流/充電電流との対応関係は、電流測定器30を構成するセンサ等の端子極性の設定によって異なるが、ここでは放電時の電池電流Ibは正(>0)であり、充電時の電池電流Ibは負(<0)であるものとする。   The current measuring unit 30 is connected between the secondary battery 10 and the connection target 20 and measures the battery current Ib. From the direction of the battery current Ib, in other words, the sign (positive / negative) of the current value Ib, the current Ib is a current flowing out of the secondary battery 10 (discharge current) or a current flowing into the secondary battery 10 (charging) Current). The correspondence between the sign of the battery current Ib and the discharge current / charge current differs depending on the terminal polarity setting of the sensor or the like constituting the current measuring device 30, but here the battery current Ib during discharge is positive (> 0). It is assumed that the battery current Ib during charging is negative (<0).

電圧測定部32は、二次電池10の正極端子と負極端子との間に接続されており、両端子間の電圧を電池電圧Vbとして測定する。   The voltage measuring unit 32 is connected between the positive terminal and the negative terminal of the secondary battery 10 and measures the voltage between both terminals as the battery voltage Vb.

なお、電圧測定部32と上記の電流測定部30とを総称して測定部と呼ぶとき、当該測定部は、それぞれが電池状態量に当たる電池電圧Vbと電池電流Ibとを測定する要素である。   In addition, when the voltage measuring unit 32 and the current measuring unit 30 are collectively referred to as a measuring unit, the measuring unit is an element that measures the battery voltage Vb and the battery current Ib corresponding to the battery state quantities.

温度測定部34は、例えば温度検知部が二次電池10の外面に取り付けられることにより、電池温度Tbを測定する。   The temperature measurement unit 34 measures the battery temperature Tb, for example, by attaching a temperature detection unit to the outer surface of the secondary battery 10.

測定部30,32,34は、監視制御部50に接続されており、当該監視制御部50へ測定値Ib,Vb,Tbを送出する。監視制御部50は、受信した測定値Ib,Vb,Tbに基づいて二次電池10の電池抵抗が上昇する現象を監視し、必要に応じて二次電池10や接続対象20を制御する。   The measurement units 30, 32, and 34 are connected to the monitoring control unit 50 and send measured values Ib, Vb, and Tb to the monitoring control unit 50. The monitoring control unit 50 monitors the phenomenon in which the battery resistance of the secondary battery 10 increases based on the received measurement values Ib, Vb, and Tb, and controls the secondary battery 10 and the connection target 20 as necessary.

監視制御部50は、例えばマイクロコンピュータ、記憶装置等を含んで構成することが可能であり、例えばHEVに搭載される各種ECU(Electric Control Unit)の1つまたは複数を利用して構成することが可能である。この場合、監視制御部50の後述する各部(図16参照)をソフトウェアによって実現可能である。具体的には、監視制御部50の後述する各部の機能を、上記マイクロプロセッサが所定のプログラムを実行することにより実現される。当該プログラムは記憶装置に格納される。記憶装置は、例えばROM(Read Only Memory)、RAM(Random Access Memory)、書き換え可能な不揮発性メモリ(EPROM(Erasable Programmable ROM)等)、ハードディスク装置、DVD(Digital Versatile Disc)等の1つまたは複数で構成可能である。記憶装置は、各種のデータ、情報等も格納し、また、プログラムを実行するための作業領域を提供する。なお、監視制御部50の後述する各部の一部または全部をハードウェアによって実現することも可能である。監視制御部50については後に詳述する。   The monitoring control unit 50 can be configured to include, for example, a microcomputer, a storage device, and the like. For example, the monitoring control unit 50 can be configured using one or more of various ECUs (Electric Control Units) mounted on the HEV. Is possible. In this case, each part (refer FIG. 16) mentioned later of the monitoring control part 50 is realizable by software. Specifically, the function of each unit described later of the monitoring control unit 50 is realized by the microprocessor executing a predetermined program. The program is stored in a storage device. The storage device is one or more of, for example, ROM (Read Only Memory), RAM (Random Access Memory), rewritable nonvolatile memory (EPROM (Erasable Programmable ROM), etc.), hard disk device, DVD (Digital Versatile Disc), etc. Can be configured. The storage device stores various data, information, and the like, and provides a work area for executing the program. It should be noted that some or all of the units described below of the monitoring control unit 50 can be realized by hardware. The monitoring controller 50 will be described in detail later.

ここで、図1に示した概略電池モデルに適用される基礎的な電池モデル式について説明する。当該基礎的な電池モデル式は上記非特許文献1に開示されるものである。当該基礎的な電池モデル式は、以下の(1)〜(11)式からなる基礎方程式により示される。なお、電池モデル式で用いられる変数および定数の一覧表を図6に示す。   Here, a basic battery model formula applied to the schematic battery model shown in FIG. 1 will be described. The basic battery model formula is disclosed in Non-Patent Document 1 above. The basic battery model equation is represented by a basic equation consisting of the following equations (1) to (11). A list of variables and constants used in the battery model equation is shown in FIG.

なお、以下に説明するモデル式中の変数および定数について、添字eは電解液中の値であることを示し、sは活物質中の値であることを示す。また、添字jは正極および負極を区別するものであり、j=1は正極における値を示し、j=2は負極における値を示すものとする。なお、正極および負極での変数または定数を包括的に表記する場合には、添字jを省略して表記することもある。また、時間の関数であることを示す(t)の表記、電池温度の依存性を示す(T)の表記、または、局所充電率(SOC:State Of Charge)θの依存性を示す(θ)等について、表記を省略することもある。また、変数または定数に付された記号#は、平均値を表すものとする。   For the variables and constants in the model formulas described below, the subscript e indicates a value in the electrolytic solution, and s indicates a value in the active material. The subscript j distinguishes between the positive electrode and the negative electrode, j = 1 indicates the value at the positive electrode, and j = 2 indicates the value at the negative electrode. In addition, when the variable or constant in a positive electrode and a negative electrode is described comprehensively, it may describe by omitting the subscript j. In addition, the notation (t) indicating that it is a function of time, the notation (T) indicating the dependency of the battery temperature, or the dependency of the state of charge (SOC) θ (θ) In some cases, the notation is omitted. The symbol # attached to a variable or constant represents an average value.

上記(1),(2)式は、電極(活物質)における電気化学反応を示す式であり、バトラー・ボルマーの式と呼ばれる。   The above formulas (1) and (2) are formulas indicating an electrochemical reaction in the electrode (active material), and are called Butler-Volmer formulas.

また、電解液中のリチウムイオン濃度保存則に関する式として、(3)式が成立する。一方、活物質内のリチウム濃度保存則に関する式として、(4)式の拡散方程式および(5)および(6)式に示す境界条件式が適用される。(5)式は、活物質の中心部における境界条件を示し、(6)式は、活物質の電解液との界面(以下、単に「界面」とも呼ぶ)における境界条件を示す。   In addition, the equation (3) is established as an equation relating to the conservation law of lithium ion concentration in the electrolytic solution. On the other hand, as an equation relating to the law of conservation of lithium concentration in the active material, the diffusion equation (4) and the boundary condition equation shown in equations (5) and (6) are applied. Equation (5) represents the boundary condition at the center of the active material, and equation (6) represents the boundary condition at the interface between the active material and the electrolyte (hereinafter also simply referred to as “interface”).

ここで、活物質界面における局所的なリチウム濃度分布である局所SOCθjは(7)式で定義される。(7)式中のCsejは、(8)式に示されるように、正極および負極の活物質界面でのリチウム濃度を示している。また、csj,maxは、活物質内での限界リチウム濃度を示している。 Here, the local SOC θ j that is a local lithium concentration distribution at the active material interface is defined by equation (7). C sej in the equation (7) indicates the lithium concentration at the active material interface between the positive electrode and the negative electrode, as shown in the equation (8). Further, c sj, max indicates the limit lithium concentration in the active material.

さらに、電解液中の電荷保存則に関する式として(9)式が成立し、活物質中の電荷保存則に関する式として(10)式が成立する。また、活物質界面での電気化学反応式として、電流密度I(t)と、反応電流密度jj Liとの関係を示す(11)式が成立する。 Further, the equation (9) is established as an equation relating to the charge conservation law in the electrolytic solution, and the equation (10) is established as an equation relating to the charge conservation law in the active material. Further, as an electrochemical reaction formula at the active material interface, formula (11) showing a relationship between the current density I (t) and the reaction current density j j Li is established.

上記(1)〜(11)式の基礎方程式は、非特許文献1に開示されたものであるが、実機に搭載されて二次電池の状態をオンラインで推定する場合、これらの電池モデル式をそのまま適用することは、監視制御部50の演算負荷および演算時間の制約上、極めて困難である。このため、監視制御部50では上記電池モデル式を以下に説明するように簡易化して利用する。   The basic equations of the above formulas (1) to (11) are those disclosed in Non-Patent Document 1. However, when the state of the secondary battery is estimated on-line when mounted on an actual machine, these battery model formulas are expressed as follows. It is extremely difficult to apply it as it is because of restrictions on the computation load and computation time of the monitoring controller 50. For this reason, the supervisory control unit 50 uses the battery model formula in a simplified manner as described below.

図7に電池モデルを説明する概念図を示す。本実施の形態では、負極12中および正極15中における電気化学反応を一様なものと仮定する。すなわち、各電極12,15においてx方向(図1参照)には反応が均一に生じるものと仮定する。また、負極12内および正極15内で各活物質の反応が均一と仮定するので、負極12および正極15のそれぞれにおいて1個の活物質モデルを取扱うことにする。この結果、図1に示した二次電池の概略構造は、図7に示すようにモデリングされる。   FIG. 7 is a conceptual diagram for explaining the battery model. In the present embodiment, it is assumed that the electrochemical reaction in the negative electrode 12 and the positive electrode 15 is uniform. That is, it is assumed that the reaction occurs uniformly in the x direction (see FIG. 1) in each of the electrodes 12 and 15. In addition, since it is assumed that the reaction of each active material is uniform in the negative electrode 12 and the positive electrode 15, one active material model is handled in each of the negative electrode 12 and the positive electrode 15. As a result, the schematic structure of the secondary battery shown in FIG. 1 is modeled as shown in FIG.

放電時には、負極活物質モデル18nの表面での電極反応により、活物質モデル18n内のリチウム原子Liが、電子e-の放出によりリチウムイオンLi+となってセパレータ14中の電解液に放出される。一方、正極活物質モデル18pの表面の電極反応では、電解液中のリチウムイオンLi+が取込まれて電子e-を吸収する。これにより、正極活物質モデル18pの内部にリチウム原子Liが取り込まれる。負極活物質18nからのリチウムイオンLi+の放出および正極活物質18pでのリチウムイオンLi+の取り込みによって、正極側集電体16から負極側集電体13に向けて電流が流れる。 During discharge, due to an electrode reaction on the surface of the negative electrode active material model 18n, lithium atoms Li in the active material model 18n are released into the electrolyte in the separator 14 as lithium ions Li + by the release of electrons e −. . On the other hand, in the electrode reaction on the surface of the positive electrode active material model 18p, lithium ions Li + in the electrolytic solution are taken in and absorb the electrons e . As a result, lithium atoms Li are taken into the positive electrode active material model 18p. A current flows from the positive electrode side current collector 16 toward the negative electrode side current collector 13 by the release of the lithium ion Li + from the negative electrode active material 18n and the uptake of the lithium ion Li + by the positive electrode active material 18p.

反対に、二次電池の充電時には、負極活物質モデル18nの表面での電極反応により電解液中のリチウムイオンLi+が取り込まれ、正極活物質モデル18pの表面での電極反応では電解液ヘリチウムイオンLi+が放出される。 On the other hand, when the secondary battery is charged, lithium ions Li + in the electrolytic solution are taken in by the electrode reaction on the surface of the negative electrode active material model 18n, and in the electrode reaction on the surface of the positive electrode active material model 18p, Li + is released.

電池モデル式では、充放電時における活物質モデル18p(j=1)および活物質モデル18n(j=2)の表面での電極反応、活物質モデル18p,18nの内部でのリチウムの拡散(径方向)および電解液中のリチウムイオンの拡散(濃度分布)や、各部位での電位分布や温度分布がモデリングされる。   In the battery model formula, the electrode reaction on the surface of the active material model 18p (j = 1) and the active material model 18n (j = 2) during charging / discharging, the diffusion of lithium (diameter inside the active material models 18p, 18n) Direction) and the diffusion (concentration distribution) of lithium ions in the electrolyte, and the potential distribution and temperature distribution at each site.

各活物質モデル18p,18n内でのリチウム濃度csは、図8に示すように、周方向での位置依存性はないものと仮定して、半径方向の座標r(r:活物質18の中心からの距離、rs:活物質18の半径)上での関数として表すことができる。この活物質モデルは、界面での電気化学反応に伴う、活物質18内部でのリチウム拡散現象を推定するためのものであり、径方向にN分割(N:2以上の自然数)された各領域(k=1〜N)について、リチウム濃度cs,k(t)が、後述する拡散方程式に従って推定される。なお、図8では活物質18を等間隔でN分割する場合を例示しているが、不等間隔で分割することも可能である。 As shown in FIG. 8, the lithium concentration c s in each of the active material models 18p and 18n is assumed to have no position dependency in the circumferential direction. Distance from the center, r s : radius of the active material 18). This active material model is for estimating the lithium diffusion phenomenon inside the active material 18 due to the electrochemical reaction at the interface, and is divided into N regions (N: a natural number of 2 or more) in the radial direction. For (k = 1 to N), the lithium concentration c s, k (t) is estimated according to the diffusion equation described below. Although FIG. 8 illustrates the case where the active material 18 is divided into N at equal intervals, it is also possible to divide the active material 18 at unequal intervals.

これらの仮定の結果、非特許文献1における基礎方程式の(1)〜(6)式および(8)式は、下記(1’)〜(6’)式および(8’)式に簡易化される。   As a result of these assumptions, the equations (1) to (6) and (8) of the basic equations in Non-Patent Document 1 are simplified to the following equations (1 ′) to (6 ′) and (8 ′). The

特に(3’)式では、電解液塩濃度を時間に対して不変と仮定することによって、cej(t)が一定値であると仮定する。また、各電極12,15について1つずつ定義される活物質モデル18n,18pに対しては、拡散方程式(4)〜(6)式が極座標方向の分布のみを考慮して、拡散方程式(4’)〜(6’)式に変形される。なお、(8’)式での活物質界面でのリチウム濃度csejは、図8に示したN分割領域のうちの最外周の領域におけるリチウム濃度csj(t)に対応する。 In particular, in equation (3 ′), it is assumed that c ej (t) is a constant value by assuming that the electrolyte salt concentration remains unchanged with respect to time. For the active material models 18n and 18p defined for each of the electrodes 12 and 15, the diffusion equations (4) to (6) are considered in consideration of only the distribution in the polar coordinate direction and the diffusion equation (4 ') To (6'). Note that the lithium concentration c sej at the active material interface in the equation (8 ′) corresponds to the lithium concentration c sj (t) in the outermost region of the N-divided region shown in FIG.

さらに、電解液中の電荷保存則に関する(9)式は、電極中の電気化学反応を一様とし、電解液塩濃度を時間に対して不変と仮定する(3’)式により、下記(12)式に簡易化される。すなわち、電解液の変位φejは、xの二次関数により近似される。そして、過電圧ηj#の算出に用いる電解液中の平均電位φej#は、(12)式を電極厚さLjで積分した下記(13)式によって求められる。 Furthermore, the equation (9) relating to the law of conservation of electric charge in the electrolytic solution is expressed by the following equation (12) based on the equation (3 ′) assuming that the electrochemical reaction in the electrode is uniform and the electrolyte salt concentration is unchanged with respect to time. ). That is, the displacement φ ej of the electrolytic solution is approximated by a quadratic function of x. The overvoltage eta j # average potential phi ej # in the electrolytic solution used for the calculation of is determined by the following expression (13) obtained by integrating the electrode thickness L j (12) below.

負極12については、(12)式より下記(14)式が成立するため、電解液平均電位φe2#およびセパレータ14との境界における電解液電位との電位差は、下記(15)式で示されることとなる。同様に、正極15についても、電解液平均電位φe1#およびセパレータ14との境界における電解液電位との電位差は、下記(16)式で示される。 For the negative electrode 12, since the following formula (14) is established from the formula (12), the potential difference between the electrolyte average potential φ e2 # and the electrolyte potential at the boundary with the separator 14 is represented by the following formula (15). It will be. Similarly, for the positive electrode 15, the potential difference between the electrolyte average potential φ e1 # and the electrolyte potential at the boundary with the separator 14 is expressed by the following equation (16).

同様に、活物質中の電荷保存則に関する(10)式についても、下記(17)式に簡易化される。すなわち、活物質φsjについても、xの二次関数により近似される。そして、過電圧ηj#の算出に用いる活物質中の平均電位φsj#は、(17)式を電極厚さLjで積分した下記(18)式によって求められる。このため、電解液中の電位と同様にして、正極15における、活物質平均電位φs1#および集電体16との境界における活物質電位との電位差は、下記(19)式で示される。負極12についても同様に、(20)式が成立する。 Similarly, the formula (10) related to the law of conservation of charge in the active material is simplified to the following formula (17). That is, the active material φ sj is also approximated by a quadratic function of x. The average potential φ sj # in the active material used for calculating the overvoltage η j # is obtained by the following equation (18) obtained by integrating the equation (17) with the electrode thickness L j . For this reason, in the same manner as the potential in the electrolytic solution, the potential difference between the active material average potential φ s1 # and the active material potential at the boundary with the current collector 16 in the positive electrode 15 is expressed by the following equation (19). Similarly, the formula (20) is established for the negative electrode 12.

図9に、二次電池の端子電圧V(t)と、上述のように求めた各平均電位との関係が示す。なお、図9中において、セパレータ14では反応電流密度jj Li=0であるため、セパレータ14での電圧降下は、電流密度I(t)に比例し、Ls/Ks eff・I(t)となる。 FIG. 9 shows the relationship between the terminal voltage V (t) of the secondary battery and each average potential obtained as described above. In FIG. 9, since the reaction current density j j Li = 0 in the separator 14, the voltage drop at the separator 14 is proportional to the current density I (t), and L s / K s eff · I (t )

さらに、上述のように、各電極中における電気化学反応を一様と仮定したことにより、極板単位面積当たりの電流密度I(t)と反応電流密度(リチウム生成量)jj Liとの間には、下記(21)式が成立することとなる。 Furthermore, as described above, by assuming that the electrochemical reaction in each electrode is uniform, the current density I (t) per electrode plate unit area and the reaction current density (lithium generation amount) j j Li The following equation (21) is established.

したがって、図9に示された電位関係および上記(21)式より、電池電圧V(t)について下記(22)式が成立する。なお、(22)式は、図9に示された(23)式の電位関係を前提とするものである。   Therefore, from the potential relation shown in FIG. 9 and the above equation (21), the following equation (22) is established for the battery voltage V (t). The expression (22) is based on the potential relationship of the expression (23) shown in FIG.

次に、平均過電圧η#(t)を算出する。jj Liを一定とし、かつ、バトラー・ボルマーの関係式において、αajcj=0.5(すなわち充放電効率は同一)とすると、下記(24)式が成立する。(24)式を逆変換することにより、平均過電圧η#(t)は、下記(25)式により求められる。 Next, an average overvoltage η # (t) is calculated. When j j Li is constant and α aj = α cj = 0.5 (that is, the charge / discharge efficiency is the same) in the Butler-Volmer relational expression, the following expression (24) is established. By inversely transforming the equation (24), the average overvoltage η # (t) is obtained by the following equation (25).

そして、図9に従って平均電位φs1およびφs2を求めて(22)式に代入し、さらに、(25)式より求めた平均過電圧η1#(t)およびη2#(t)を(22)式に代入する。この結果、簡易化された電気化学反応式である(1’)、(21)式および電圧関係式(2’)に基づいて、電気化学反応モデル式に従った下記の電圧−電流関係モデル式(M1a)が導出される。 Then, the average potentials φ s1 and φ s2 are obtained according to FIG. 9 and substituted into the equation (22), and the average overvoltages η 1 # (t) and η 2 # (t) obtained from the equation (25) are expressed as (22) Substitute in the expression. As a result, based on the simplified electrochemical reaction formulas (1 ′), (21) and the voltage relation formula (2 ′), the following voltage-current relation model formula according to the electrochemical reaction model formula is shown. (M1a) is derived.

さらに、リチウム濃度保存則(拡散方程式)(4’)式および境界条件式(5’),(6’)式によって、正極15および負極12のそれぞれでの活物質モデル18p,18nについての活物質拡散モデル式(M2a)が求められる。   Furthermore, the active material for the active material models 18p and 18n in the positive electrode 15 and the negative electrode 12 respectively according to the lithium concentration conservation law (diffusion equation) (4 ′) and boundary condition equations (5 ′) and (6 ′). A diffusion model formula (M2a) is obtained.

なお、モデル式(M1a)の右辺第1項は活物質表面での反応物質(リチウム)濃度により決定される開放電圧(OCV)を示し、右辺第2項は過電圧(η1#−η2#)を示し、右辺第3項は電池電流による電圧降下を示している。すなわち、二次電池の直流純抵抗が(M1a)式中のRd(T)で表される。 The first term on the right side of the model formula (M1a) indicates the open circuit voltage (OCV) determined by the concentration of the reactant (lithium) on the active material surface, and the second term on the right side indicates the overvoltage (η 1 # −η 2 # ) And the third term on the right side shows a voltage drop due to the battery current. That is, the DC pure resistance of the secondary battery is represented by Rd (T) in the formula (M1a).

また(M2a)式において、反応物質であるリチウムの拡散速度を規定するパラメータとして用いられる拡散係数Ds1およびDs2は温度依存性を有する。このため、拡散係数Ds1およびDs2は、温度測定部34によって検出された電池温度に対して可変に設定される可変パラメータとして、例えば図10に示すようなマップに従って、設定される。図10に示されるように、拡散係数Ds1およびDs2は、電池温度の低温時に相対的に低下する一方で、温度上昇に応じて相対的に上昇する特性を有する。したがって、このような温度依存性を反映して、温度測定部34の測定値Tbに従って検知される電池温度Tに応じて、(M2a)式中の拡散係数Ds1およびDs2を電池パラメータ値として設定するマップが予め作成される。 In the equation (M2a), the diffusion coefficients D s1 and D s2 used as parameters for defining the diffusion rate of lithium as a reactant have temperature dependence. Therefore, the diffusion coefficients D s1 and D s2 are set as variable parameters that are variably set with respect to the battery temperature detected by the temperature measuring unit 34, for example, according to a map as shown in FIG. As shown in FIG. 10, the diffusion coefficients D s1 and D s2 have a characteristic of relatively decreasing when the battery temperature is low while relatively increasing with increasing temperature. Therefore, reflecting such temperature dependency, the diffusion coefficients D s1 and D s2 in the equation (M2a) are used as battery parameter values according to the battery temperature T detected according to the measurement value Tb of the temperature measurement unit 34. A map to be set is created in advance.

なお、拡散係数Ds1およびDs2は、必要に応じて、温度依存性だけでなく局所SOCθ依存性を有してもよく、その場合、検知される電池温度Tと推定される局所SOCθに応じて、これらの拡散係数をパラメータ値として設定する二次元マップが予め作成される。 Note that the diffusion coefficients D s1 and D s2 may have not only temperature dependency but also local SOC θ dependency as necessary, in which case, depending on the detected battery temperature T and the estimated local SOC θ. Thus, a two-dimensional map for setting these diffusion coefficients as parameter values is created in advance.

同様に、(M1a)式中の開放電圧U1およびU2についても、図11(a),(b)に示すように、局所SOCθの低下に伴って上昇または低下するような依存性を有する。したがって、このような局所SOC依存性を反映して、推定された局所的SOCθに応じて、(M1a)式の開放電圧U1およびU2を設定するマップが予め設定される。なお、図11(a)には、正極(j=1)における開放電圧U1の局所SOCθ依存性が示され、図11(b)には、負極(j=2)における開放電圧U2の局所SOCθ依存性が示される。 Similarly, as shown in FIGS. 11A and 11B, the open circuit voltages U 1 and U 2 in the equation (M1a) have a dependency that increases or decreases as the local SOC θ decreases. . Therefore, a map for setting open-circuit voltages U 1 and U 2 in the equation (M1a) is set in advance according to the estimated local SOC θ, reflecting such local SOC dependency. 11A shows the local SOC θ dependence of the open circuit voltage U 1 at the positive electrode (j = 1), and FIG. 11B shows the open circuit voltage U 2 at the negative electrode (j = 2). Local SOCθ dependence is shown.

さらに、図示は省略するが、(M1a)式中の交換電流密度i01およびi02についても、局所SOCθおよび電池温度Tに対する依存性を有するので、θおよびTを引数とする二次元マップを予め準備して、現時点での局所SOCθおよび電池温度Tに応じて交換電流密度i01およびi02が電池パラメータ値として決定される。 Furthermore, although illustration is omitted, since the exchange current densities i 01 and i 02 in the equation (M1a) also have a dependency on the local SOC θ and the battery temperature T, a two-dimensional map with θ and T as arguments is previously stored. In preparation, the exchange current densities i 01 and i 02 are determined as battery parameter values in accordance with the current local SOC θ and battery temperature T.

同様に、直流純抵抗Rdも温度依存性を有するので、温度依存性を反映して、温度測定部34の検出値Tbに従って検知される電池温度Tに応じて、直流純抵抗Rdを電池パラメータ値として設定するマップが予め作成される。   Similarly, since the DC pure resistance Rd also has temperature dependence, the DC pure resistance Rd is changed to the battery parameter value according to the battery temperature T detected according to the detection value Tb of the temperature measurement unit 34, reflecting the temperature dependence. A map to be set as is created in advance.

なお、上述したマップ群については、二次電池10についての、周知の交流インピーダンス測定等の実験結果に基づいて作成することが可能である。   In addition, about the map group mentioned above, it is possible to create based on experimental results, such as a known alternating current impedance measurement, about the secondary battery 10.

ここで、上記では負極12および正極15のそれぞれに別個の活物質モデルを使用する方法を例示したが(図7参照)、これに代えて、負極12および正極15に共通の活物質モデルを適用して、さらなる演算負荷の軽減を図るモデル化も可能である。この場合には、負極12および正極15の活物質モデル18n,18pをまとめて単独の素子として扱うため、下記(26)式に示すような式の置き換えが必要となる。(26)式では、正極および負極の区別を示す添字jが削除される。   Here, the method of using a separate active material model for each of the negative electrode 12 and the positive electrode 15 is illustrated above (see FIG. 7), but instead, a common active material model is applied to the negative electrode 12 and the positive electrode 15. Thus, it is possible to model to further reduce the calculation load. In this case, since the active material models 18n and 18p of the negative electrode 12 and the positive electrode 15 are collectively handled as a single element, it is necessary to replace the equation as shown in the following equation (26). In the equation (26), the suffix j indicating the distinction between the positive electrode and the negative electrode is deleted.

この結果、モデル式(M1a)および(M2a)をさらに簡易化したモデルとして、下記(M1b)式および(M2b)式が得られる。また、このように正極および負極を共通化した活物質モデルを適用した電池モデル式において、電流密度I(t)と反応電流密度jj Liとの関係式は(21)式に代えて(21’)式が適用される。 As a result, the following equations (M1b) and (M2b) are obtained as models obtained by further simplifying the model equations (M1a) and (M2a). Further, in the battery model formula to which the active material model in which the positive electrode and the negative electrode are shared is applied in this way, the relational expression between the current density I (t) and the reaction current density j j Li is replaced with the formula (21) (21 ') The formula is applied.

ところで、上述の(M1a)式では、過電圧を示す右辺第2項にarcsinh項が存在するため、非線形式を解く必要が生じる。このため、(M1a)式の演算には繰返し計算が必要となり、演算負荷が増大する他、演算の安定性を損なう可能性もある。したがって、他の変形例として、(M1a)式中のarcsinh項を一次近似(線形近似)した電圧−電流関係モデル式(M1c)式を導出することが可能である。   By the way, in the above equation (M1a), since the arcsinh term exists in the second term on the right side indicating overvoltage, it is necessary to solve the nonlinear equation. For this reason, the calculation of equation (M1a) requires repeated calculation, which increases the calculation load and may impair the stability of the calculation. Therefore, as another modified example, it is possible to derive a voltage-current relationship model equation (M1c) that is a first-order approximation (linear approximation) of the arcsinh term in the equation (M1a).

なお、(M1c)式では、線形近似の結果、右辺第2項も電流密度I(t)と反応抵抗Rrとの積で示されることとなる。反応抵抗Rrは、上記(27)式に示されるように、交換電流密度i01,i02についての局所SOCθおよび電池温度Tの依存性を有する。したがって、電圧−電流関係モデル式(M1c)式を用いる場合には、反応抵抗Rr(θi,T)についても、局所SOCθおよび電池温度Tに対する依存性を反映したマップを実験結果等に基づき予め作成して、電池パラメータの1つとしてマップ参照により適宜設定することが必要となる。 In the equation (M1c), as a result of linear approximation, the second term on the right side is also represented by the product of the current density I (t) and the reaction resistance Rr. The reaction resistance Rr is dependent on the local SOC θ and the battery temperature T with respect to the exchange current densities i 01 and i 02 as shown in the above equation (27). Therefore, when the voltage-current relationship model equation (M1c) is used, a map reflecting the dependence on the local SOC θ and the battery temperature T for the reaction resistance Rr (θ i , T) is preliminarily determined based on experimental results. It is necessary to create and appropriately set one of the battery parameters by referring to the map.

この結果、電圧−電流関係モデル式(M1c)式は、上記(28)式に示すように電池電圧V(t)、電池電流密度I(t)、電池パラメータとしてのRr(θ,T)およびRd(T)、ならびに開放電圧U1,U2の間に成立する線形モデル式として表現されることとなる。 As a result, the voltage-current relationship model equation (M1c) is expressed by the battery voltage V (t), the battery current density I (t), Rr (θ, T) as the battery parameters, as shown in the above equation (28). It will be expressed as a linear model equation established between Rd (T) and the open circuit voltages U 1 and U 2 .

同様に、上述の(M1b)式についても、右辺第2項のarcsinh項を線形近似することによって、下記(M1d)式が得られる。(M1d)式についても、(M1c)式と同様に線形モデル式として表現される。   Similarly, the following equation (M1d) can be obtained by linearly approximating the arcsinh term of the second term on the right side of the above equation (M1b). The expression (M1d) is also expressed as a linear model expression in the same manner as the expression (M1c).

上記電池モデルに基づき、監視制御部50(図5参照)は、上記電圧−電流関係モデル式(M1a)〜(M1d)のいずれかと、これに対応する活物質拡散モデル(M2a)または(M2b)との組合わせによって実現される、電池状態量推定部52を含んでいる。ここでは、電池状態量推定部52が、電圧測定部32によって測定された電池電圧Vbに基づき、上記電池モデルに従って電池電流密度I(t)を推定する場合を例示するが、上記電圧−電流関係モデル式(M1a)〜(M1d)に従えば、電流測定部30によって測定された電池電流Ibに基づいて(測定電池電流値Ibを電極表面積で除算することによって電池電流密度を算出可能である)、電池電圧V(t)を推定するように電池状態量推定部52を構成することも可能である。   Based on the battery model, the supervisory control unit 50 (see FIG. 5) can select one of the voltage-current relationship model equations (M1a) to (M1d) and the corresponding active material diffusion model (M2a) or (M2b). The battery state quantity estimation part 52 implement | achieved by the combination with these is included. Here, a case where the battery state quantity estimation unit 52 estimates the battery current density I (t) according to the battery model based on the battery voltage Vb measured by the voltage measurement unit 32 is illustrated. According to the model formulas (M1a) to (M1d), based on the battery current Ib measured by the current measuring unit 30 (the battery current density can be calculated by dividing the measured battery current value Ib by the electrode surface area). It is also possible to configure the battery state quantity estimation unit 52 so as to estimate the battery voltage V (t).

図12は電池状態量推定部52を説明する概略構成図である。図12の例示では、電池状態量推定部52は、拡散推定部100と、開放電圧推定部110と、電流推定部120と、電池パラメータ値設定部130と、境界条件設定部140とを含む。   FIG. 12 is a schematic configuration diagram illustrating the battery state quantity estimation unit 52. In the example of FIG. 12, the battery state quantity estimation unit 52 includes a diffusion estimation unit 100, an open-circuit voltage estimation unit 110, a current estimation unit 120, a battery parameter value setting unit 130, and a boundary condition setting unit 140.

拡散推定部100は、活物質拡散モデル式(M2a)または(M2b)式により、境界条件設定部140によって(5’)および(6’)式に従って設定された境界条件に基づいて、活物質内部でのリチウム濃度分布を、例えば差分形式により逐次演算する。また、拡散推定部100は、推定されたリチウム濃度分布に基づき、最外周の領域におけるリチウム濃度を活物質界面でのリチウム濃度csejとして、(7)式に従って局所SOCθを設定する。 Based on the boundary conditions set by the boundary condition setting unit 140 according to the equations (5 ′) and (6 ′) using the active material diffusion model equation (M2a) or (M2b), The lithium concentration distribution at is sequentially calculated, for example, in a differential format. In addition, based on the estimated lithium concentration distribution, diffusion estimation unit 100 sets the local SOC θ according to the equation (7), with the lithium concentration in the outermost region as the lithium concentration c sej at the active material interface.

開放電圧推定部110は、図11に示された特性を反映したマップに従い、正極および負極それぞれの開放電圧、または正極および負極を合成した開放電圧を求める。なお、図12中では、これらを包括的に開放電圧U(θ)と表記している。   The open-circuit voltage estimation unit 110 obtains the open-circuit voltage of each of the positive electrode and the negative electrode or the open-circuit voltage obtained by combining the positive electrode and the negative electrode according to a map reflecting the characteristics shown in FIG. In FIG. 12, these are comprehensively expressed as an open circuit voltage U (θ).

電池パラメータ値設定部130は、温度測定部34の検出値Tbに従って検知される電池温度T、および、拡散推定部100による推定に基づく現在の局所SOCθに応じて、使用する電池モデル式中の電池パラメータを設定する。上述のように、拡散推定部100で用いられる、モデル式(M2a),(M2b)中の拡散係数Ds1,Ds2,Dsが電池温度Tに応じて設定される他、モデル式(M1a)〜(M1d)中の直流純抵抗Rd、モデル式(M1a),(M1b)中の交換電流密度i01,i02、モデル式(M1c),(M1d)中の反応抵抗Rr等が、電池パラメータ値設定部130により設定される。 The battery parameter value setting unit 130 uses the battery temperature T detected according to the detection value Tb of the temperature measurement unit 34 and the current local SOC θ based on the estimation by the diffusion estimation unit 100. Set the parameters. As described above, the diffusion coefficients D s1 , D s2 , and D s in the model equations (M2a) and (M2b) used in the diffusion estimation unit 100 are set according to the battery temperature T, and the model equation (M1a ) To (M1d) DC pure resistance Rd, exchange current densities i 01 and i 02 in model equations (M1a) and (M1b), reaction resistances Rr in model equations (M1c) and (M1d), etc. It is set by the parameter value setting unit 130.

電流推定部120は、上述の(M1a)〜(M1d)式に基づいて、電池電流密度I(t)を算出するための下記(M3a)〜(M3d)式のいずれかに従って、開放電圧推定部110によって推定された開放電圧U(θ)と、電圧測定部34の検出値Vbに従って検知される現在の電池電圧V(t)とを代入して、電池電流密度I(t)を算出する。   The current estimation unit 120 is based on any one of the following formulas (M3a) to (M3d) for calculating the battery current density I (t) based on the above formulas (M1a) to (M1d). The battery current density I (t) is calculated by substituting the open-circuit voltage U (θ) estimated by 110 and the current battery voltage V (t) detected according to the detection value Vb of the voltage measurement unit 34.

例えば、(M3a)式は、(M1a)式を電池電流密度I(t)について解いたものに相当する。非線形方程式である(M3a)式をニュートン法で解く場合には、I(t)の初期値を仮定し、さらに、電池電圧V(t)、開放電圧U(θ)および電池パラメータ値を代入してI(t)を算出して、算出したI(t)と仮定したI(t)とが一致するまで反復計算を行なうことにより解くことができる。   For example, the equation (M3a) corresponds to a solution of the equation (M1a) with respect to the battery current density I (t). When solving the equation (M3a) which is a nonlinear equation by the Newton method, an initial value of I (t) is assumed, and further, the battery voltage V (t), the open circuit voltage U (θ) and the battery parameter value are substituted. I (t) is calculated, and the calculation can be solved by performing iterative calculation until the calculated I (t) matches the assumed I (t).

同様に、(M1b)式を用いた電池モデル式では、下記(M3b)式を(M3a)式と同様の手法で解くことにより、電池電流密度I(t)を算出できる。   Similarly, in the battery model equation using the equation (M1b), the battery current density I (t) can be calculated by solving the following equation (M3b) by the same method as the equation (M3a).

また、線形近似された(M1c),(M1d)式を電池モデル式では、下記(M3c),(M3d)式に従って、非線形方程式のような反復計算を行なうことなく、電池電圧V(t)、開放電圧U(θ)および電池パラメータ値を代入して一意に電池電流密度I(t)を算出できる。これにより、演算負荷の軽減および、演算の安定性向上が図られる。   Further, in the battery model equation, the linearly approximated equations (M1c) and (M1d) are expressed in accordance with the following equations (M3c) and (M3d) without performing iterative calculation like a nonlinear equation, and the battery voltage V (t), The battery current density I (t) can be uniquely calculated by substituting the open circuit voltage U (θ) and the battery parameter value. Thereby, calculation load can be reduced and calculation stability can be improved.

次に、境界条件設定部140は、上記(21)式または(21’)式に従って、演算された電流密度I(t)を反応電流密度(リチウム生成量)jj Liに変換して、(6’)式に従って活物質拡散モデル式(M2a),(M2b)式の境界条件を更新する。 Next, the boundary condition setting unit 140 converts the calculated current density I (t) into the reaction current density (lithium generation amount) j j Li according to the above formula (21) or (21 ′), The boundary conditions of the active material diffusion model equations (M2a) and (M2b) are updated according to the equation 6 ′).

このように、簡易化された電気化学反応モデル式の採用によって演算負荷が軽減されるとともに、二次電池10の内部状態を電気化学反応に基づいて高精度に推定可能な電池モデル式を実現できる。このため、かかる電池モデル式によれば、実機搭載に適したモデルを提供可能である。   As described above, the calculation load is reduced by adopting the simplified electrochemical reaction model formula, and the battery model formula capable of estimating the internal state of the secondary battery 10 with high accuracy based on the electrochemical reaction can be realized. . For this reason, according to such a battery model formula, a model suitable for mounting on an actual machine can be provided.

図13に電池状態量推定部52による処理を説明するフローチャートを例示する。図13に示される一連の処理は、監視制御部50において所定演算周期ごとに呼び出されて実行される。   FIG. 13 illustrates a flowchart for explaining processing by the battery state quantity estimation unit 52. A series of processing shown in FIG. 13 is called and executed by the monitoring control unit 50 every predetermined calculation cycle.

ステップS100では、電池状態量推定部52は、電圧測定部32による検出値Vbに基づき電池電圧V(t)を取得し、ステップS110において、温度測定部34による検出値Tbに基づき現在の電池温度Tを取得する。   In step S100, the battery state quantity estimation unit 52 acquires the battery voltage V (t) based on the detection value Vb by the voltage measurement unit 32, and in step S110, the current battery temperature based on the detection value Tb by the temperature measurement unit 34. Get T.

ステップS120では、電池状態量推定部52は、拡散モデル式(M2a),(M2b)式による前回の演算時におけるリチウム濃度分布に基づき、活物質表面のリチウム濃度に基づいて局所的SOCθを算出する。ステップS120の処理は拡散推定部100の機能に相当する。   In step S120, the battery state quantity estimation unit 52 calculates the local SOC θ based on the lithium concentration on the surface of the active material, based on the lithium concentration distribution at the previous calculation by the diffusion model equations (M2a) and (M2b). . The process of step S120 corresponds to the function of the diffusion estimation unit 100.

ステップS130では、電池状態量推定部52は、図11に示した特性に従うマップの参照により、負極12および正極15ごとに、または正負極共通に、局所的SOCθから開放電圧U(θ)値を算出する。ステップS130の処理は開放電圧推定部110の機能に相当する。   In step S130, the battery state quantity estimation unit 52 calculates the open circuit voltage U (θ) value from the local SOC θ for each of the negative electrode 12 and the positive electrode 15 or in common to the positive and negative electrodes by referring to the map according to the characteristic shown in FIG. calculate. The process of step S130 corresponds to the function of the open-circuit voltage estimation unit 110.

ステップS140では、電池状態量推定部52は、電流推定部120の機能により、電池電圧V(t)と、開放電圧U(θ)と、電池温度Tに応じて設定された電池パラメータ値とに応じて、電圧−電流関係モデル式(M1a)〜(M1d)のいずれかに従って、電池電流密度I(t)の推定値を算出する。   In step S140, the battery state quantity estimation unit 52 uses the function of the current estimation unit 120 to change the battery voltage V (t), the open circuit voltage U (θ), and the battery parameter value set according to the battery temperature T. Accordingly, an estimated value of battery current density I (t) is calculated according to any one of voltage-current relationship model formulas (M1a) to (M1d).

さらに、電池状態量推定部52は、ステップS150において、推定された電池電流密度I(t)から、(21)式または(21’)式に基づいて、反応電流密度(リチウム生成量)jj Liを算出するとともに、算出した反応電流密度を用いて拡散モデル方程式(M2a),(M2b)の活物質界面における境界条件(活物質界面)を設定する。ステップS150の処理は境界条件設定部140の機能に相当する。 Further, in step S150, the battery state quantity estimation unit 52 determines the reaction current density (lithium generation amount) j j from the estimated battery current density I (t) based on the formula (21) or (21 ′). While calculating Li , the boundary condition (active material interface) in the active material interface of diffusion model equations (M2a) and (M2b) is set using the calculated reaction current density. The processing in step S150 corresponds to the function of the boundary condition setting unit 140.

そして、ステップS160では、電池状態量推定部52は、拡散方程式モデル(M2a),(M2b)に従って、負極12および正極15ごとに、または正負極共通に設定された活物質モデル内のリチウム濃度分布を計算し、各領域のリチウム濃度推定値を更新する。ステップS160の処理は拡散推定部100の機能に相当する。なお、上述のように、この際演算された最外周の分割領域におけるリチウム濃度は、次回の演算実行時に、活物質表面での局所的SOCθの算出に用いられる(ステップS120参照)。   In step S160, the battery state quantity estimating unit 52 determines the lithium concentration distribution in the active material model set for each of the negative electrode 12 and the positive electrode 15 or in common with the positive and negative electrodes according to the diffusion equation models (M2a) and (M2b). And update the estimated lithium concentration in each region. The process of step S160 corresponds to the function of the diffusion estimation unit 100. As described above, the lithium concentration in the outermost divided region calculated at this time is used to calculate the local SOC θ on the active material surface when the next calculation is performed (see step S120).

なお、上記の電池モデル式は、全電池電流が活物質18を流れて電気化学反応に関与するとの前提で導出されている。しかしながら、実際には、特に低温時等において、電解液と活物質との界面に電気二重層キャパシタが生じ、その影響によって全電池電流が、電気化学反応に関与する電気化学反応電流と、キャパシタを流れるキャパシタ電流とに分流される場合がある。そのような場合には例えば、電気二重層キャパシタを通過するキャパシタ電流成分を電気化学反応電流成分と分離するように、電池モデル式を構成するのが、より好ましい。   The above battery model equation is derived on the assumption that the entire battery current flows through the active material 18 and participates in the electrochemical reaction. However, in reality, an electric double layer capacitor is generated at the interface between the electrolytic solution and the active material, particularly at low temperatures, and the total battery current is affected by the electrochemical reaction current involved in the electrochemical reaction and the capacitor. In some cases, the current is shunted to the flowing capacitor current. In such a case, for example, it is more preferable to configure the battery model formula so that the capacitor current component passing through the electric double layer capacitor is separated from the electrochemical reaction current component.

ここで、上記のように、非特許文献1に開示される基礎的な電池モデルでは、電極12,15の表面におけるリチウムイオンLi+の反応、電極12,15の活物質18内におけるリチウムイオンLi+の拡散、および電解液内でのリチウムイオンLi+の拡散がモデル化されている。これに対し、電池状態量推定部52に適用される簡易化電池モデルは、上記基礎的な電池モデルにおいて、電極厚さ方向の反応は一様であるとする仮定、および、電極12,15内でのリチウムイオンLi+の濃度は一定であるとする仮定の下で構成されている。 Here, as described above, in the basic battery model disclosed in Non-Patent Document 1, the reaction of lithium ions Li + on the surfaces of the electrodes 12 and 15 and the lithium ions Li in the active material 18 of the electrodes 12 and 15 are performed. + Diffusion and lithium ion Li + diffusion in the electrolyte are modeled. On the other hand, in the simplified battery model applied to the battery state quantity estimation unit 52, in the basic battery model, it is assumed that the reaction in the electrode thickness direction is uniform. It is configured under the assumption that the concentration of lithium ion Li + in is constant.

電解液中のリチウムイオンLi+の濃度、すなわち電解液塩濃度が十分に高い場合には、簡易化電池モデルの上記仮定を満足することは可能であり、その結果、電池電流密度I(t)を精度良く推定可能である。これは、上記のように、電解液塩濃度が十分に高い場合には、充放電によって電極内電解液塩濃度が変化したとしても、当該塩濃度変化が反応抵抗に及ぼす影響が小さいからである(図2および図3参照)。 If the lithium ion Li + concentration in the electrolyte, ie, the electrolyte salt concentration, is sufficiently high, it is possible to satisfy the above assumption of the simplified battery model and, as a result, the battery current density I (t) Can be estimated with high accuracy. This is because, as described above, when the electrolyte salt concentration is sufficiently high, even if the electrolyte salt concentration in the electrode changes due to charge / discharge, the influence of the salt concentration change on the reaction resistance is small. (See FIGS. 2 and 3).

これに対し、簡易化電池モデルにおける上記仮定は電極内の電解液塩濃度が低い場合に生じる反応抵抗の上昇(図3参照)が考慮されていないので、簡易化電池モデルによって推定された電池電流密度I(t)と、実際に測定された電池電流Ibに基づく電流密度との間には誤差が生じる。   On the other hand, the assumption in the simplified battery model does not take into account the increase in reaction resistance (see FIG. 3) that occurs when the electrolyte salt concentration in the electrode is low, so the battery current estimated by the simplified battery model is not considered. There is an error between the density I (t) and the current density based on the actually measured battery current Ib.

かかる点に鑑みれば、当該誤差に基づいて電池抵抗の上昇の度合い(指標)を推定することが可能である。例えば、電解液塩濃度(ここではリチウムイオン濃度)の拡散方程式を簡易化することにより、電極内の電解液における塩濃度変化は次の(51),(52)式によって推定することができる。   In view of this point, it is possible to estimate the degree of increase (index) of battery resistance based on the error. For example, by simplifying the diffusion equation of the electrolyte salt concentration (here, lithium ion concentration), the salt concentration change in the electrolyte solution in the electrode can be estimated by the following equations (51) and (52).

(51),(52)式において、Δceは負極内電解液塩濃度と正極内電解液塩濃度との差(図2参照)であり、Deffは電解液の有効拡散係数であり、εeは電解液の体積分率であり、t+ 0はリチウムイオンLi+の輸率であり、Fはファラデー定数であり、Δtは当該推定処理の時間間隔(時間刻み)であり、Δxは拡散距離(図2参照)であり、Tは電池温度、Iは電池電流密度である。なお、塩濃度差Δceは、例えば放電時については、図2に示すように、負極での塩濃度増加量(平均塩濃度から放電時塩濃度への増加量)と、正極での塩濃度減少量(平均塩濃度から放電時塩濃度への減少量)との合計として把握される。かかる点に鑑みれば、塩濃度差Δceは平均塩濃度からの変化量であると捉えることも可能である。 (51), in equation (52), .DELTA.c e is the difference between the negative electrode in the electrolyte salt concentration and the positive electrode in the electrolyte salt concentration (see FIG. 2), D eff is the effective diffusion coefficient of the electrolyte solution, epsilon e is the volume fraction of the electrolyte, t + 0 is the transport number of lithium ion Li + , F is the Faraday constant, Δt is the time interval (time step) of the estimation process, and Δx is the diffusion The distance (see FIG. 2), T is the battery temperature, and I is the battery current density. Incidentally, the salt concentration difference .DELTA.c e, for example the time of discharge, as shown in FIG. 2, with the salt concentration increase at the negative electrode (the amount of increase in the average salt concentration to the discharge time of the salt concentration), the salt concentration of the positive electrode It is grasped as the sum of the amount of decrease (the amount of decrease from the average salt concentration to the salt concentration during discharge). In view of this point, the salt concentration difference .DELTA.c e can be also be understood that the amount of change from the average salt concentration.

式(51)、(52)によって推定された電極間電解液塩濃度差Δceと、電流推定誤差(Im−Ir)(Imは推定電流密度、Irは測定電流Ibに基づく電流密度)との相関を図14に示す。図14によれば、電極間電解液塩濃度差Δceが大きい領域において、電流推定誤差が大きくなる現象が認められる。 The interelectrode electrolyte salt concentration difference Δc e estimated by the equations (51) and (52) and the current estimation error (Im−Ir) (Im is the estimated current density, Ir is the current density based on the measured current Ib) The correlation is shown in FIG. According to FIG 14, in the region between the electrodes electrolyte salt concentration difference .DELTA.c e is large, a phenomenon is observed that the current estimation error increases.

したがって、電極間電解液塩濃度差Δceが大きい場合における電流推定誤差(Im−Ir)の値を、電池抵抗上昇度合いとして利用することができる。なお、電極間電解液塩濃度差Δceが大きいという条件は、例えば、Δceの値が予め設定された所定値以上であるという条件、または、予め設定された所定範囲内に存在するという条件として設定することが可能である。 Therefore, the value of the current estimation error (Im-Ir) when the inter-electrode electrolyte salt concentration difference .DELTA.c e is large, can be utilized as a battery resistance increasing degree. The condition that the condition that the inter-electrode electrolyte salt concentration difference .DELTA.c e is large, for example, the condition that the value of .DELTA.c e is preset predetermined value or more, or is present within a preset predetermined range Can be set as

また、電流推定誤差(Im−Ir)すなわちImとIrとの差分の代わりに、推定電流密度Imと測定電流に基づく電流密度Irとの比の値によっても、電池抵抗上昇度合いを表すことは可能である。かかる点に基づき、抵抗上昇率Dh(=Im/Ir)と電極間電解液塩濃度差Δceとの相関を図15の模式図に示す。図15においても、電極間電解液塩濃度差Δceが大きい領域において、電池抵抗が高くなる現象が認められる。 Further, instead of the current estimation error (Im-Ir), that is, the difference between Im and Ir, the degree of battery resistance increase can be expressed by the value of the ratio between the estimated current density Im and the current density Ir based on the measured current. It is. Based on this point, shown in the schematic diagram of FIG. 15 and the resistance increase rate Dh (= Im / Ir) correlations between the inter-electrode electrolyte salt concentration difference .DELTA.c e. Also in FIG. 15, in the region between the electrodes electrolyte salt concentration difference .DELTA.c e is large, a phenomenon is observed that cell resistance increases.

ここで、電極間電解液塩濃度差Δceが大きい領域において電流推定誤差(Im−Ir)が発生するのは、電極内電解液塩濃度の低下により発生する電池抵抗上昇分が実際の二次電池10と電池モデルとで異なるからであると解釈することができる。その一方で、電池抵抗上昇に起因する電圧変化量ΔVは、実際の二次電池10と電池モデルとで等しい。このため、実際に発現する電池抵抗の増加分をRrとし、電池モデルにおける電池抵抗の増加分をRmとすると、次の(53)式が成り立ち、当該(53)式から(54)式が得られる。 Here, the current estimation error (Im-Ir) is generated in the region between the electrodes electrolyte salt concentration difference .DELTA.c e is large, two battery resistance rise caused by lowering of the electrode in the electrolytic solution salt concentration actual primary It can be interpreted that the battery 10 and the battery model are different. On the other hand, the voltage change amount ΔV due to the battery resistance increase is equal between the actual secondary battery 10 and the battery model. For this reason, assuming that the actual increase in battery resistance is Rr and the increase in battery resistance in the battery model is Rm, the following equation (53) is established, and equation (54) is obtained from the equation (53). It is done.

したがって、推定電流密度Imと測定電流に基づく電流密度Irとの比は、電池抵抗上昇率Dhを表すことになる。 Therefore, the ratio between the current density Ir based on the estimated current density Im and the measured current will represent a battery resistance increase rate D h.

また、電極間電解液塩濃度差Δceに対する電流推定誤差(Im−Ir)の増加率、換言すれば電極間電解液塩濃度差Δeと電流推定誤差(Im−Ir)とについての一次近似式における直線の傾きによっても、電池抵抗上昇度合いを表すことが可能である。同様に、電極間電解液塩濃度差Δceに対する電流密度比(Im/Ir)の増加率によっても、電池抵抗上昇度合いを表すことが可能である。 The rate of increase of the current estimation error with respect to the inter-electrode electrolyte salt concentration difference Δc e (Im-Ir), i.e. the inter-electrode electrolyte salt if density difference delta e and the current estimation error and (Im-Ir) first-order approximation of the The degree of increase in battery resistance can also be expressed by the slope of the straight line in the equation. Similarly, the rate of increase in current density ratio between the electrodes electrolyte salt concentration difference Δc e (Im / Ir), it is possible to represent the battery resistance increases degree.

また、電極間電解液塩濃度差Δceと電流推定誤差(Im−Ir)とについての一次近似式、または、電極間電解液塩濃度差Δceと電流密度比(Im/Ir)とについての一次近似式(図15中の一点鎖線を参照)において、Δce=0とした場合の値(いわゆる切片値)によっても、電池抵抗上昇度合いを表すことが可能である。 Further, a linear approximation formula for the interelectrode electrolyte salt concentration difference Δc e and the current estimation error (Im−Ir), or the interelectrode electrolyte salt concentration difference Δc e and the current density ratio (Im / Ir) In the primary approximate expression (see the alternate long and short dash line in FIG. 15), the battery resistance increase degree can also be expressed by a value (so-called intercept value) when Δc e = 0.

なお、上記ではImを推定電流密度とし、Irを測定電流に基づく電流密度としたが、Imを推定された電流密度に電極表面積を乗算して得られる電流値とし、Irを測定電流値Ibとしても、上記説明はあてはまる。   In the above, Im is an estimated current density and Ir is a current density based on the measured current. However, Im is a current value obtained by multiplying the estimated current density by the electrode surface area, and Ir is a measured current value Ib. However, the above description applies.

また、上記では電池抵抗の上昇現象が発生しうる領域、換言すれば当該現象の発現を監視するのに好適な領域を、電極間電解液塩濃度差Δceが大きいという上記条件で設定したが、当該条件に代えて、または当該条件に加えて、充放電電流値が大きい領域という条件を用いてもよい。充放電電流値が大きい領域という条件は、例えば、測定電流値が予め設定された所定値以上であるという条件、または、予め設定された所定範囲内に存在するという条件として設定することが可能である。また、測定電流値に代えて、当該測定値に基づいて算出される電流密度、または、測定電流値もしくは電流密度の移動平均値を利用して同様の条件を設定することも可能である。 A region rising phenomenon of the battery resistance may occur in the above, a suitable region to monitor the expression of the phenomenon in other words, has been set in the condition that a large inter-electrode electrolyte salt concentration difference .DELTA.c e Instead of the condition or in addition to the condition, a condition of a region where the charge / discharge current value is large may be used. The condition that the charge / discharge current value is large can be set, for example, as a condition that the measured current value is greater than or equal to a preset predetermined value, or a condition that the measured current value is within a preset predetermined range. is there. Further, instead of the measured current value, it is possible to set the same condition by using the current density calculated based on the measured value or the measured current value or the moving average value of the current density.

上記考察に基づいた監視制御部50の概略構成を図16に例示する。監視制御部50は、上記の電池状態量推定部52に加え、電池抵抗上昇度合い推定部54を含んでいる。電池抵抗上昇度合い推定部54は、ここでは、電池状態量推定部52の電流推定部120(図12参照)が推定した電流密度I(t)(上記の推定電流密度Imに相当する)と、電流測定部30が測定した電池電流Ibとを取得し、これらの値I(t),Ibに基づいて電池抵抗上昇度合いDgを推定する。例えば、電池抵抗上昇度合い推定部54は、電池電流測定値Ibから電流密度Irを算出し、電流密度Im,Irを用いて上記の各種形式で表現可能な電池抵抗上昇度合いDgを算出する。 FIG. 16 illustrates a schematic configuration of the monitoring control unit 50 based on the above consideration. The monitoring control unit 50 includes a battery resistance increase degree estimation unit 54 in addition to the battery state quantity estimation unit 52 described above. Here, the battery resistance increase degree estimation unit 54 includes a current density I (t) (corresponding to the estimated current density Im described above) estimated by the current estimation unit 120 (see FIG. 12) of the battery state quantity estimation unit 52, and The battery current Ib measured by the current measuring unit 30 is acquired, and the battery resistance increase degree Dg is estimated based on these values I (t) and Ib. For example, the battery resistance increase degree estimation unit 54 calculates the current density Ir from the battery current measurement value Ib, and calculates the battery resistance increase degree Dg that can be expressed in the various forms described above using the current densities Im and Ir.

また、電池抵抗上昇度合い推定部54は、さらに、電極間電解液塩濃度差Δceの推定((51),(52)式参照)と、測定電池電流Ibまたは電流密度Irの移動平均値の算出との少なくとも1つを実行し、当該実行結果にも基づいて電池抵抗上昇度合いDgを推定することが好ましい。より具体的には、電極間電解液塩濃度差Δceと電流Ib,Irの移動平均値との一方または両方によって電池抵抗の上昇現象が発生する領域(上記の塩濃度差Δceが大きい領域または/および電流Ib,Irの大きい領域)を設定し、電池抵抗上昇度合いDgの算出に用いる値Im,Ir等が当該設定領域内のデータであるか否かを判定するのが好ましい。 Also, the battery resistance increasing degree estimation unit 54, further, estimates the inter-electrode electrolyte salt concentration difference Δc e ((51), ( 52) reference expression) and, moving average value of the measured battery current Ib or current density Ir performing at least one of the calculation, it is preferable to estimate the battery resistance increases the degree D g also based on the execution result. More specifically, the inter-electrode electrolyte salt concentration difference .DELTA.c e and the current Ib, Ir, or both rising phenomenon of the battery resistance is generated by the region (region above the salt concentration difference .DELTA.c e is large and the moving average of the or / and the current Ib, and set a large area) of Ir, a value used to calculate the battery resistance increasing degree D g Im, Ir or the like is preferable to determine whether the data of the set area.

このように監視領域を設定することにより、電池抵抗上昇度合いDgの信頼性が向上する。すなわち、図14に示されるように電池抵抗上昇度合いDg(図14では電流推定誤差(Im−Ir)が図示される)には、任意の塩濃度差Δceについて、ある程度のばらつきが生じうる。このため、電池抵抗上昇現象の監視に適した領域を予め選定することによって、電池抵抗上昇度合いDgの信頼性を向上させることができる。また、かかる監視領域を設けずに、換言すれば電極間電解液塩濃度差Δceとは無関係に、電池抵抗上昇度合いDgを継続的に算出する構成に比べて、演算負荷を軽減することができる。 By setting the monitoring area, thereby improving the reliability of the battery resistance increases the degree D g. That is, as shown in FIG. 14, the battery resistance increase degree D g (the current estimation error (Im−Ir) is shown in FIG. 14) may vary to some extent with respect to an arbitrary salt concentration difference Δc e. . Therefore, by pre-selecting a region suitable for monitoring the battery resistance increases phenomenon, it is possible to improve the reliability of the battery resistance increases the degree D g. Further, without providing such a monitoring region, in other words, regardless of the interelectrode electrolyte salt concentration difference Δc e , the calculation load can be reduced as compared with the configuration in which the battery resistance increase degree D g is continuously calculated. Can do.

このとき、監視領域設定に電流Ib,Irの移動平均値を用いる場合には、ノイズによる電流Ib,Irの変動が抑制されるので、領域判定の精度が高くなる。また、上記のように測定電流値Ibまたはこれに基づく電流密度Irによっても、領域判定をすることは可能であり、この場合には移動平均値を用いる構成と比較して、演算負荷が小さくて済む。   At this time, when the moving average values of the currents Ib and Ir are used for monitoring region setting, fluctuations in the currents Ib and Ir due to noise are suppressed, so that the accuracy of region determination is increased. Further, as described above, the region can also be determined by the measured current value Ib or the current density Ir based on the measured current value Ib. In this case, the calculation load is small compared to the configuration using the moving average value. That's it.

図17に、電池抵抗上昇度合いの推定に関して、監視制御部50による処理を説明するフローチャートを例示する。但し、電池抵抗上昇度合いの推定処理はかかる例示に限定されるものではない。図17に示される一連の処理は、監視制御部50において所定演算周期ごとに呼び出されて実行される。   FIG. 17 illustrates a flowchart for explaining processing by the monitoring control unit 50 regarding the estimation of the battery resistance increase degree. However, the battery resistance increase degree estimation process is not limited to this example. A series of processing shown in FIG. 17 is called and executed by the monitoring control unit 50 every predetermined calculation cycle.

ステップS200では、監視制御部50は、電流測定部30が測定した電池電流Ibと、電圧測定部32が測定した電池電圧Vbと、温度測定部34が測定した電池温度Tbとを取得する。なお、測定部30,32,34による測定は同じサンプリングタイミングで行われる。   In step S200, the monitoring control unit 50 acquires the battery current Ib measured by the current measurement unit 30, the battery voltage Vb measured by the voltage measurement unit 32, and the battery temperature Tb measured by the temperature measurement unit 34. Note that the measurement by the measurement units 30, 32, and 34 is performed at the same sampling timing.

ステップS210では、監視制御部50は、測定値Vb,Tbに基づき上記の電池モデルに従って電池電流密度Imを推定する。ステップS210での処理は電池状態量推定部52の上記機能に相当する(図12および図13参照)。   In step S210, the monitoring controller 50 estimates the battery current density Im according to the battery model based on the measured values Vb and Tb. The process in step S210 corresponds to the above function of the battery state quantity estimation unit 52 (see FIGS. 12 and 13).

ステップS220において、監視制御部50は、上記(51),(52)式に基づいて電極間電解液塩濃度差Δceを算出する。(51)式中の電池電流密度I(t)として、ここでは測定電流値Ibに基づいて算出される電流密度Irを用いるものとする。 In step S220, the monitoring control unit 50, the (51), calculates the inter-electrode electrolyte salt concentration difference .DELTA.c e based on the equation (52). As the battery current density I (t) in the equation (51), here, the current density Ir calculated based on the measured current value Ib is used.

そして、監視制御部50は、ステップS230において、塩濃度差Δceの推定値が予め設定された所定範囲内にあるという条件(ce1<Δce<Δce2)と、測定値に基づく電池電流密度Irが予め設定された所定値I0以上であるという条件(Ir>I0)との両方を満足するか否かを判定する。監視制御部50は、両条件を満足すると判定した場合には、ステップS235において、上記の電池抵抗上昇率Dhの推定に用いるデータを不図示のRAM等に保存する。 In step S230, the monitoring control unit 50 determines that the estimated value of the salt concentration difference Δc e is within a predetermined range (c e1 <Δc e <Δc e2 ) and the battery current based on the measured value. It is determined whether or not both of the conditions (Ir> I0) that the density Ir is equal to or higher than a predetermined value I0 set in advance are satisfied. Monitoring control unit 50, when it is determined that satisfies the both conditions in step S235, stores the data used to estimate the above-mentioned battery resistance increase rate D h in the RAM (not shown) or the like.

ここでは、電池抵抗上昇率Dhを最小自乗法によって推定する場合を例示する。この場合、(54)式から下記(55)式が得られることに鑑み、(56)式で与えられる中間変数Sx,Sy,Sxx,Sxyを算出して保存する。 Here, a case where the battery resistance increase rate D h is estimated by the method of least squares is illustrated. In this case, considering that the following equation (55) is obtained from equation (54), intermediate variables S x , S y , S xx , and S xy given by equation (56) are calculated and stored.

また、ステップS235では、監視制御部50は、中間変数Sx,Sy,Sxx,Sxyの組(または群)の取得回数に対応するカウンタ値Mの値を1カウント増加する(インクリメントする)。ステップS235の実行後、処理は上記ステップS200へ移行する。 In step S235, the monitoring controller 50 increments (increments) the value of the counter value M corresponding to the number of acquisitions (or groups) of the intermediate variables S x , S y , S xx , S xy by one. ). After execution of step S235, the process proceeds to step S200.

他方、監視制御部50は、上記ステップS230において両条件を同時には満足しないと判定した場合には、ステップS240においてカウンタ値Mが、予め設定された設定値M0よりも大きいという条件(M>M0)を満足するか否かを判定する。なお、設定値M0は、最小自乗法による推定に利用するデータ数(データ組数)を規定するものである。設定値Mが大きいほど推定精度は向上する一方で、設定値Mが小さいほど演算負荷を小さくすることができるので、例えば監視制御部50を構成するマイクロプロセッサの処理能力に応じて適宜、設定すればよい。   On the other hand, if the monitoring control unit 50 determines in step S230 that both conditions are not satisfied at the same time, the condition that the counter value M is larger than the preset setting value M0 in step S240 (M> M0). ) Is satisfied. The set value M0 defines the number of data (number of data sets) used for estimation by the least square method. As the setting value M is larger, the estimation accuracy is improved. On the other hand, the smaller the setting value M is, the smaller the calculation load is. Therefore, for example, it is set as appropriate according to the processing capability of the microprocessor constituting the monitoring control unit 50. That's fine.

ステップS240での判定の結果、M>M0という条件を満たさない場合、監視制御部50の処理は上記ステップS200へ移行する。これに対し、M>M0という条件を満たす場合、監視制御部50は、ステップS250において、上記ステップS235において逐次更新された中間変数Sx,Sy,Sxx,Sxyを下記(57)式に代入することによって、電池抵抗上昇率Dhを推定する。 As a result of the determination in step S240, if the condition of M> M0 is not satisfied, the process of the monitoring control unit 50 proceeds to step S200. On the other hand, when the condition of M> M0 is satisfied, the monitoring control unit 50 sets the intermediate variables S x , S y , S xx , S xy sequentially updated in step S235 in step S250 to the following formula (57). By substituting into, the battery resistance increase rate D h is estimated.

なお、(57)式においてMを2倍しているのは、(55)式の近似直線における原点との傾きを求めるためである。   The reason why M is doubled in the equation (57) is to obtain the inclination of the approximate straight line of the equation (55) with respect to the origin.

また、監視制御部50は、電池抵抗上昇率Dhの推定ばらつきによる影響を抑制するために、ステップS260において電池抵抗上昇率Dhに緩変化処理を施す。緩変化処理後の電池抵抗上昇率には符号Dhfを用いることにする。緩変化処理は次の(58)式に従って行うことが可能である。なお、(58)式中のkは緩変化処理係数であり、0<k<1の範囲から予め設定される。 The monitoring control unit 50, in order to suppress the effect of the estimated variation in battery resistance increase rate D h, subjected to gradual change process to the battery resistance increase rate D h in step S260. The symbol D hf is used as the battery resistance increase rate after the slow change process. The slow change process can be performed according to the following equation (58). Note that k in the equation (58) is a slow change processing coefficient, and is preset from the range of 0 <k <1.

その後、監視制御部50は、ステップS270において中間変数Sx,Sy,Sxx,Sxy、カウンタ値M等をクリアして、次回の推定処理に備える。 Thereafter, the supervisory control unit 50 clears the intermediate variables S x , S y , S xx , S xy , the counter value M, and the like in step S270 to prepare for the next estimation process.

上記ステップS220〜S270の処理は電池抵抗上昇度合い推定部54の機能に相当する。   The processes in steps S220 to S270 correspond to the function of the battery resistance increase degree estimation unit 54.

上記手法を用いて、実際の電動車両における電池電流、電池電圧および電池温度の時系列データから、電池抵抗上昇率Dhfを推定した結果を図18に示す。図18によれば、電池抵抗上昇率Dhfが時間経過とともに増加している様子が確認できる(図18中の矢印参照)。 FIG. 18 shows a result of estimating the battery resistance increase rate D hf from the time series data of the battery current, the battery voltage, and the battery temperature in the actual electric vehicle using the above method. According to FIG. 18, it can be confirmed that the battery resistance increase rate D hf increases with time (see the arrow in FIG. 18).

図19および図20に、電池電流の時間変化(各図の(a)参照)、電極間電解液塩濃度差Δceの時間変化(各図の(b)参照)、および、電極間電解液塩濃度差Δceと電流推定誤差(Im−Ir)との相関図(各図の(c)参照)を示す。図19と図20とは各種データの測定時期が異なり、図20の方が後の時間におけるものである。図19に示されるように測定初期では、塩濃度差Δceの大きい領域であっても電流推定誤差(Im−Ir)はほとんど生じていない。これに対し、図20に示されるように測定時間の後期では、塩濃度差Δceの大きい領域において電流推定誤差(Im−Ir)が大きくなっている。 19 and 20, (see (a) of each figure) the time change of the battery current, (see (b) of each figure) time variation of the inter-electrode electrolyte salt concentration difference .DELTA.c e, and, the inter-electrode electrolyte solution correlation diagram between the salt concentration difference .DELTA.c e and the current estimation error (Im-Ir) (in each figure (c) refer) shows a. 19 and 20 differ in the measurement timing of various data, and FIG. 20 is at a later time. The initial measurement as shown in Figure 19, a large area even when the current estimation error of the salt concentration difference Δc e (Im-Ir) hardly occurs. In contrast, in the late measurement time, as shown in FIG. 20, the current estimation error (Im-Ir) is increased in the region of high salt concentration difference .DELTA.c e.

電池抵抗上昇度合いDgを利用して種々の制御や処理を行うことが可能であり、かかる制御等の例を上記図16を参照して説明する。 Using the battery resistance increases the degree D g it is possible to perform various control and processing, an example of such control and the like will be described with reference to FIG 16.

図16の例では、電池抵抗上昇度合いDgが、予め設定された値(通知設定値)よりも大きい場合、通知部56によって電池抵抗の上昇を通知する。通知部56による通知は例えば視覚的、聴覚的に行うことが可能であり、例えば、アラーム音や音声等を発生する機器、ランプ、表示デバイス等によって通知部56を構成可能である。通知部56は、監視装置5用に設けてもよいし、例えば車両に搭載された機器と兼用にしてもよい。かかる構成によれば、二次電池10のメンテナンスの必要性をユーザや保守担当者等に知らせることができる。 In the example of FIG. 16, when the battery resistance increase degree Dg is larger than a preset value (notification setting value), the notification unit 56 notifies the battery resistance increase. The notification by the notification unit 56 can be performed, for example, visually or audibly. For example, the notification unit 56 can be configured by a device, a lamp, a display device, or the like that generates an alarm sound or a sound. The notification unit 56 may be provided for the monitoring device 5, or may be shared with a device mounted on the vehicle, for example. According to such a configuration, it is possible to notify the user, maintenance staff, etc. of the necessity of maintenance of the secondary battery 10.

また、図16の例では、監視制御部50は、電池抵抗上昇度合いDgが、予め設定された値(使用制限設定値)よりも大きい場合に、二次電池10の使用を制限する使用制限部58を含んでいる。使用制限部58による電池使用制限は、例えば充電電力量、放電電力量、充電電流量、放電電流量、電池電圧上限値、電池電圧下限値等を制限することによって行われる。充電電力量等の制限は公知の各種手法を利用することによって行うことが可能であり、使用制限部58は必要に応じて二次電池10や接続対称20(図5参照)を制御する。かかる構成によれば、電池抵抗上昇現象の進行を抑制することができる。 In the example of FIG. 16, the monitoring control unit 50 restricts use of the secondary battery 10 when the battery resistance increase degree Dg is larger than a preset value (usage restriction setting value). Part 58 is included. The battery use restriction by the use restriction unit 58 is performed by restricting, for example, a charge power amount, a discharge power amount, a charge current amount, a discharge current amount, a battery voltage upper limit value, a battery voltage lower limit value, and the like. The amount of charging power can be limited by using various known methods, and the use limiting unit 58 controls the secondary battery 10 and the connection symmetry 20 (see FIG. 5) as necessary. According to such a configuration, the progress of the battery resistance increase phenomenon can be suppressed.

また、図16の例では、監視制御部50はSOC制御部60を含んでいる。例えばHEVに搭載される二次電池については、モータを駆動可能な電池容量を常時確保するとともに、回生電力を受け入れられるだけの空き容量を常時残しておくという制御が採用される場合がある。この場合、二次電池のSOCは電池容量の例えば40%〜80%の範囲内、すなわち60%を中心として±20%の範囲になるように制御される。監視制御部50のSOC制御部60は例えば、このようなSOC制御を行うとともに、電池抵抗上昇度合いDgが、予め設定された値(SOC制御中心調整設定値)よりも大きい場合にはSOC制御中心値を調整(変更)する。 In the example of FIG. 16, the monitoring control unit 50 includes an SOC control unit 60. For example, with respect to a secondary battery mounted on an HEV, there may be a case where control is performed such that a battery capacity capable of driving a motor is always secured and a free capacity that can accept regenerative power is always left. In this case, the SOC of the secondary battery is controlled to be within a range of, for example, 40% to 80% of the battery capacity, that is, within a range of ± 20% centering on 60%. SOC control unit 60 of the monitoring control unit 50 is, for example, performs such SOC control, battery resistance increasing degree D g is, SOC control is larger than a predetermined value (SOC control center adjustment setting) Adjust (change) the center value.

このとき、SOC制御中心値を、予め設定された通常設定値よりも下げることによって、上記通常設定値での制御時と比較して、より多くの電力を充電させることが可能になる。これにより、速やかに電池抵抗を正常状態に回復させることができる。また、SOC制御中心値を通常設定値よりも上げることによって、上記通常設定値での制御時と比較して、より多くの電力を放電させることができる。これにより、速やかに電池抵抗を正常状態に回復させることが出来ると共に、放電時に電池電圧が下限電圧に到達するのを抑制することができる。   At this time, by lowering the SOC control center value than the normal set value set in advance, it becomes possible to charge more electric power as compared with the control at the normal set value. Thereby, the battery resistance can be promptly restored to a normal state. Further, by raising the SOC control center value above the normal set value, it is possible to discharge more electric power compared to the control at the normal set value. Thereby, the battery resistance can be quickly restored to a normal state, and the battery voltage can be prevented from reaching the lower limit voltage during discharging.

また、図16の例では、監視制御部50は、二次電池10について各種の診断を行う診断部62を含んでいる。なお、例えばHEVには補機を含む各種部品が正常に動作するか否かを診断・検出するダイアグ機能を有するECUが設けられており、二次電池10がHEVに搭載される場合には当該ECUを利用して診断部62を構成することも可能である。監視制御部50の診断部62は、電池抵抗上昇度合いDgが、予め設定された値(電圧下限値診断禁止設定値)よりも大きい場合に、二次電池10の電池電圧の下限値についての故障診断判定(ダイアグ判定)は行わないように、構成されている。かかる構成によれば、放電時に電池電圧が下限値まで低下した場合であっても、経年劣化による電池抵抗上昇が発生したとする誤った判定が行われるのを防止することができる。 In the example of FIG. 16, the monitoring control unit 50 includes a diagnosis unit 62 that performs various diagnoses on the secondary battery 10. For example, the HEV is provided with an ECU having a diagnosis function for diagnosing and detecting whether or not various parts including the auxiliary machine operate normally. When the secondary battery 10 is mounted on the HEV, It is also possible to configure the diagnosis unit 62 using the ECU. When the battery resistance increase degree Dg is larger than a preset value (voltage lower limit value diagnosis prohibition setting value), the diagnosis unit 62 of the monitoring control unit 50 determines the lower limit value of the battery voltage of the secondary battery 10. Failure diagnosis determination (diagnosis determination) is not performed. According to such a configuration, it is possible to prevent an erroneous determination that an increase in battery resistance due to aging has occurred even when the battery voltage drops to the lower limit value during discharge.

また、図16の例では、監視制御部50は、二次電池10の充電または放電を強制的に行う強制的充放電制御部64を含んでいる。当該強制的充放電制御部64は、電池抵抗上昇度合いDgが、予め設定された値(充放電要求設定値)よりも大きい場合に、上記の強制的な充電または放電を行う。これにより、電解液塩濃度を調整して、電池抵抗上昇現象を低減・解消することができる。このとき、上記の使用制限部58による電池使用制限に比べて、速やかな電解液塩濃度の調整が可能である。 In the example of FIG. 16, the monitoring control unit 50 includes a forced charge / discharge control unit 64 that forcibly charges or discharges the secondary battery 10. The forced charge / discharge control unit 64 performs the above-described forced charge or discharge when the battery resistance increase degree Dg is larger than a preset value (charge / discharge request set value). As a result, the electrolyte salt concentration can be adjusted to reduce or eliminate the battery resistance increase phenomenon. At this time, as compared with the battery use restriction by the use restriction unit 58, the electrolyte salt concentration can be quickly adjusted.

二次電池10がHEVに搭載される例の場合、強制的な放電は例えば、強制的充放電制御部64の制御の下、二次電池10を用いたモータ駆動走行に切り換えることによって行うことが可能である。また、強制的な充電は例えば、強制的充放電制御部64がエンジンの駆動によってM/Gを発電機として機能させることにより行うことが可能である。また、強制的な充放電は、強制的充放電制御部64が、エンジンに対して要求するパワー(エンジン要求パワー)を増加または減少させることによって行うことが可能である。   In the example in which the secondary battery 10 is mounted on the HEV, forcible discharge can be performed by switching to motor-driven running using the secondary battery 10 under the control of the forced charge / discharge control unit 64, for example. Is possible. Further, forcible charging can be performed, for example, by causing the M / G to function as a generator by the forced charge / discharge control unit 64 by driving the engine. The forced charge / discharge can be performed by increasing or decreasing the power (engine required power) required by the forced charge / discharge control unit 64 for the engine.

例えばHEVにおいては、一般に、複数個の二次電池10が直列接続された電池パックが用いられる。かかる電池パックに対して強制的な充放電を行うと、強制的な充放電を要しない二次電池10では、換言すれば電解液塩濃度が所定の正常範囲にある二次電池10では、電解液塩濃度が正常範囲から逸脱してしまうことがある。かかる逸脱は、電池温度が定常状態にない場合に起こりやすい。このため、電池パックに対する強制的な充放電は、各二次電池10の温度に基づいて行うことが可能である。しかし、電池パック中の全ての二次電池10に温度センサ等を設けるのは、装置の複雑化、コスト増加を招いてしまう。   For example, in HEV, a battery pack in which a plurality of secondary batteries 10 are connected in series is generally used. When the battery pack is forcibly charged and discharged, in the secondary battery 10 that does not require forcible charge and discharge, in other words, in the secondary battery 10 in which the electrolyte salt concentration is in a predetermined normal range, The liquid salt concentration may deviate from the normal range. Such deviation is likely to occur when the battery temperature is not in a steady state. For this reason, forcible charging / discharging of the battery pack can be performed based on the temperature of each secondary battery 10. However, providing temperature sensors or the like for all the secondary batteries 10 in the battery pack leads to complication of the apparatus and an increase in cost.

これらの点に鑑みれば、電池パックに温度測定器を設け、当該温度測定器による測定温度に基づいて強制的な充放電を制御するのが好ましい。より具体的には、強制的充放電制御部64は、電池パックの温度が予め設定された温度範囲内にある場合(すなわち予め定められた定常状態にある場合)には、各二次電池10の温度は定常状態にあるものと判定して、電池パックに対する強制的な充放電を許容する。逆に、強制的充放電制御部64は、電池パックの温度が定常状態にない場合には、強制的な充放電を行わない。かかる構成によれば、電解液塩濃度が正常範囲から逸脱するのを防止することができ、しかもそのような構成を簡便、低コストで実現することができる。   In view of these points, it is preferable to provide a temperature measuring device in the battery pack and control forced charge / discharge based on the temperature measured by the temperature measuring device. More specifically, the forcible charge / discharge control unit 64 determines that each secondary battery 10 is in a state where the temperature of the battery pack is within a preset temperature range (that is, in a predetermined steady state). Is determined to be in a steady state, and forcible charging / discharging of the battery pack is permitted. Conversely, the forced charge / discharge control unit 64 does not perform forced charge / discharge when the temperature of the battery pack is not in a steady state. According to such a configuration, the electrolyte salt concentration can be prevented from deviating from the normal range, and such a configuration can be realized simply and at low cost.

なお、電池パックの温度を測定する温度測定器を、温度測定部34(図5参照)として利用することも可能である。   Note that a temperature measuring device that measures the temperature of the battery pack can be used as the temperature measuring unit 34 (see FIG. 5).

なお、上記説明では通知部56が電池抵抗上昇度合いDgを受信して通知設定値との比較判定処理を行う構成を例示したが、電池抵抗上昇度合い推定部54が当該比較判定処理を行い、その判定結果に基づいて通知部56を制御するように構成することも可能である。かかる点は、使用制限部58、SOC制御部60、診断部62、および強制的充放電制御部64についても同様である。また、通知設定値、使用制限設定値、SOC制御中心調整設定値、電圧下限値診断禁止設定値、および充放電要求設定値は、全てが異なる値であってもよいし、一部について同じ値であってもよい。 Incidentally, in the above description, but the notification portion 56 illustrates the configuration for comparing processing for determining whether the notification settings to receive the battery resistance increases the degree D g, the battery resistance increases degree estimation unit 54 performs the comparison judgment processing, The notification unit 56 may be configured to be controlled based on the determination result. The same applies to the use restriction unit 58, the SOC control unit 60, the diagnosis unit 62, and the forced charge / discharge control unit 64. Further, the notification set value, the use limit set value, the SOC control center adjustment set value, the voltage lower limit value diagnosis prohibition set value, and the charge / discharge request set value may all be different values, or the same value for some of them. It may be.

図21に、監視制御部50が、上記電池モデルで使用するパラメータを二次電池10の経時劣化に応じて補償するパラメータ変化推定部66を含む場合の構成を例示する。なお、図21では電池抵抗上昇度合い推定部54等の図示は省略している。パラメータ変化推定部66は、上記電池モデルにおける抵抗成分のパラメータについて二次電池10の経時劣化に伴う変化を推定し、当該経時劣化推定結果を電池モデルに反映させる。パラメータ変化推定部66の構成例は後述する。かかる構成によれば、経年劣化に伴う電池抵抗(反応抵抗(直流抵抗)および拡散抵抗)の変化の影響を排除して、上記の大電流域での電池抵抗上昇現象を精度良く検出することができる。   FIG. 21 illustrates a configuration in a case where the monitoring control unit 50 includes a parameter change estimation unit 66 that compensates for parameters used in the battery model according to deterioration with time of the secondary battery 10. In FIG. 21, the battery resistance increase degree estimation unit 54 and the like are not shown. The parameter change estimation unit 66 estimates a change associated with the deterioration of the secondary battery 10 with respect to the parameter of the resistance component in the battery model, and reflects the result of the deterioration with time in the battery model. A configuration example of the parameter change estimation unit 66 will be described later. According to such a configuration, it is possible to accurately detect the battery resistance increase phenomenon in the above-described large current range by eliminating the influence of changes in battery resistance (reaction resistance (DC resistance) and diffusion resistance) due to aging. it can.

パラメータ変化推定部66は、例えば経時劣化の推定結果を記憶装置に蓄積することによって、経時劣化の傾向を学習することが可能である。これに対し、パラメータ変化推定部66による経時劣化の学習が不十分である場合、電池抵抗上昇度合いDgの推定を行わないようにしてもよい。 The parameter change estimation unit 66 can learn the tendency of deterioration over time by, for example, accumulating the estimation result of deterioration over time in a storage device. In contrast, if the learning of the time degradation due to parameter change estimation unit 66 is insufficient, it may not be performed to estimate the battery resistance increasing degree D g.

例えば、電池抵抗上昇度合いDgの推定に用いる今回の電池温度測定値に対して十分な推定学習実績が無い場合、経時劣化による影響を精度良く推定し排除することは難しいことがある。かかる点は、電池電流または電池電圧の測定値について十分な推定学習実績が無い場合も同様である。このため、例えば経時劣化推定結果の蓄積データ数が予め設定された所定数に満たない場合には電池抵抗上昇度合いDgの推定を行わないことによって、上記電池抵抗上昇現象の発生を誤って検出するのを防止することができる。 For example, when there is not sufficient estimated learning results for the current battery temperature measurement value used for estimating the battery resistance increase degree Dg , it may be difficult to accurately estimate and eliminate the influence of deterioration over time. The same applies to the case where there is no sufficient estimated learning record for the measured value of the battery current or the battery voltage. Thus, for example, by not performing the estimation of battery resistance increasing degree D g in the case of storing the number of data time degradation estimation result is less than a predetermined number set in advance, falsely detect the occurrence of the battery resistance increases phenomenon Can be prevented.

また、例えば、電池抵抗上昇度合いDgの推定に用いる今回の電池温度測定値に対して経時劣化推定結果の過去のデータが大きなばらつきを有する場合も、経時劣化による影響を精度良く推定し排除することは難しいことがある。かかる点は、電池電流または電池電圧の測定値に対して経時劣化推定結果データのばらつきが大きい場合も同様である。このため、例えば経時劣化推定結果の蓄積データの全数が予め設定された所定範囲内に存在しない場合、または、例えば当該蓄積データが上記所定範囲内に収束していない場合には電池抵抗上昇度合いDgの推定を行わないことによって、上記誤検出を防止することができる。なお、蓄積データが収束しているか否かは、当該蓄積データのうちで上記所定範囲内に存在するデータの個数が、予め設定された割合以上であるか否かによって判定することが可能である。 Further, for example, even if the historical data of the time-dependent deterioration estimation results for the current battery temperature measurements used to estimate the battery resistance increasing degree D g has a large variation, accurately estimated to eliminate the influence of degradation over time It can be difficult. The same applies to the case where the variation of the time degradation estimation result data is large with respect to the measured value of the battery current or the battery voltage. For this reason, for example, when the total number of accumulated data of the estimation results of deterioration with time does not exist within a predetermined range set in advance, or when the accumulated data does not converge within the predetermined range, for example, the battery resistance increase degree D By not estimating g, the erroneous detection can be prevented. Note that whether or not the accumulated data has converged can be determined by whether or not the number of pieces of data existing within the predetermined range of the accumulated data is equal to or greater than a preset ratio. .

すなわち、電池抵抗上昇度合い推定部54は、パラメータ変化推定部66による経時劣化推定結果に基づいて、電池抵抗上昇度合いDgの推定を実行するか否かを判定するのが好ましい。 That is, the battery resistance increasing degree estimation unit 54, based on the time degradation estimation result by the parameter change estimating unit 66, preferably determine whether to execute the estimation of battery resistance increasing degree D g.

以下に、電池モデル中のパラメータの経時劣化について説明する。まず、負極12および正極15での電子e-の移動に対する純電気的な抵抗(純抵抗)Rdと、活物質界面での反応電流発生時に等価的に電気抵抗として作用する反応抵抗(電荷移動抵抗)Rrとを合わせたものが、二次電池10をマクロに見た場合の直流抵抗Raに相当する。このとき、直流抵抗Raと、電池電圧Vと、電池温度Tと、電池電流密度Iと、開放電圧(OCV)Uとの間には(71)式が成り立つ。なお、(71)式は例えば(M1d)式に対応する。 Hereinafter, deterioration with time of parameters in the battery model will be described. First, a pure electrical resistance (pure resistance) Rd against the movement of electrons e at the negative electrode 12 and the positive electrode 15 and a reaction resistance (charge transfer resistance) that acts as an electrical resistance equivalently when a reaction current is generated at the active material interface. ) a combination of the Rr corresponds to the DC resistance R a when viewed rechargeable battery 10 to the macro. At this time, the DC resistance R a, the battery voltage V, and the battery temperature T, a battery current density I, between the open circuit voltage (OCV) U (71) expression holds. The expression (71) corresponds to, for example, the expression (M1d).

開放電圧Uは、上記図11に示すように、活物質表面における局所的SOCθに依存して変化する特性を有する。したがって、二次電池10の初期状態において局所的SOCθと開放電圧Uとの関係を測定することにより、図11に示す特性を局所SOCθの変化に対する開放電圧U(θ)の変化特性を予め記憶する特性マップを作成することができる。   As shown in FIG. 11, the open circuit voltage U has a characteristic that changes depending on the local SOC θ on the surface of the active material. Therefore, by measuring the relationship between local SOC θ and open circuit voltage U in the initial state of secondary battery 10, the characteristics shown in FIG. 11 are stored in advance as the characteristics of change in open circuit voltage U (θ) with respect to changes in local SOC θ. A characteristic map can be created.

また、直流抵抗Raは、図22に示すように、局所SOCθおよび電池温度Tの変化に応じて変化する特性を有する。すなわち、直流抵抗Raは、局所SOCθおよび電池温度Tの関数として与えられる。したがって、二次電池10の初期状態における実測実験結果に基づき、図22に例示した特性に従って、局所的SOCθと電池温度Tとの組合わせに対応して直流抵抗Raの値を決定する特性マップ(直流抵抗マップ)を作成することができる。 Also, DC resistance R a, as shown in FIG. 22, has a characteristic that varies in response to changes in the local SOCθ and battery temperature T. That is, DC resistance Ra is given as a function of local SOC θ and battery temperature T. Thus, based on the measured experimental results in the initial state of the secondary battery 10, according to the characteristics illustrated in FIG. 22, the characteristic map for determining the value of DC resistance R a corresponding to the combination of local SOCθ and battery temperature T (DC resistance map) can be created.

直流抵抗Raは、上記のように電池温度Tおよび局所SOCθに応じて変化するとともに、電池使用に伴う経年的な劣化によっても変化する。このため、直流抵抗マップに記憶された初期状態(代表的には新品時)における直流抵抗Raと、実際の直流抵抗Raとの間には、経年劣化に伴って、差異が生じる。かかる直流抵抗Raの差異は、上記の電池抵抗上昇度合いに推定誤差を発生させることになる。 DC resistance R a, together with changes according to the battery temperature T and the local SOCθ as described above, also vary according to secular deterioration accompanying the battery used. Therefore, the DC resistance R a in the initial state stored in the DC resistance map (when typically new), between the actual DC resistance R a, with the aging, differences may occur. Such a difference in the DC resistance Ra causes an estimation error in the battery resistance increase degree.

なお、初期状態とは新品時に限定されるものではなく、例えば、新品時と予想される最大劣化時との間の中間的な状態に対応して、初期状態を定義してもよい。このようにすると、以下の説明で明らかになるように、初期状態パラメータ値に対する変化率推定値の範囲を狭くできるので、推定精度の向上を図ることができる。   Note that the initial state is not limited to a new product. For example, the initial state may be defined corresponding to an intermediate state between a new product and an expected maximum deterioration. This makes it possible to narrow the range of the rate of change estimated value with respect to the initial state parameter value, as will be apparent from the following description, so that the estimation accuracy can be improved.

ここでは、上記パラメータ変化推定部66(図21参照)が、直流抵抗Raについて、初期状態パラメータ値に対する変化率grを推定することによって、直流抵抗の経年変化を推定する場合を説明する。推定した直流抵抗変化率grを電池モデルに反映することによって、経年変化により直流抵抗が増加した場合においても、精度よく電池抵抗上昇度合いが推定されることとなる。 Here, the parameter change estimating unit 66 (see FIG. 21) is, for DC resistance R a, by estimating the rate of change gr with respect to the initial state parameter value, a case of estimating the aging of DC resistance. By reflecting the estimated DC resistance change rate gr in the battery model, even when the DC resistance increases due to secular change, the battery resistance increase degree can be accurately estimated.

直流抵抗変化率grを次の(72)式で定義する。なお、Ranは初期状態での直流抵抗Raのことである。 The DC resistance change rate gr is defined by the following equation (72). Incidentally, R an, is that the DC resistance R a in the initial state.

かかる直流抵抗変化率grを推定することにより、二次電池の使用に伴う経年的な直流抵抗Raの変化を、電池温度Tbや局所SOCθの変化による直流抵抗の変化と切離して推定することができる。 By estimating such DC resistance change rate gr, a change in the secular DC resistance R a with the use of the secondary battery, it is estimated separately from the change in DC resistance due to the change in battery temperature Tb and local SOCθ it can.

図21を参照して、パラメータ変化推定部66による直流抵抗変化率grの推定学習を説明する。   With reference to FIG. 21, the estimation learning of the DC resistance change rate gr by the parameter change estimation unit 66 will be described.

図21において、パラメータ特性マップ220は、電池モデル式中のパラメータである拡散係数Dsおよび直流抵抗Raについて、初期状態における実測結果に基づいた特性マップを格納している。すなわち、電池温度Tbおよび局所SOCθ等の時々刻々変化する電池状態に対応して、現時点での電池状態に対応する拡散係数Dsおよび直流抵抗Ra(Ran)を読出可能なように構成されている。パラメータ特性マップ220は例えば監視制御部50を構成する上記記憶装置に記憶される。 In Figure 21, the parameter characteristic map 220, the diffusion coefficient D s and DC resistance R a is a parameter in the battery model equation, and stores the characteristic map based on the actual measurement result in the initial state. That is, it is configured such that the diffusion coefficient D s and the DC resistance R a (R an ) corresponding to the current battery state can be read in response to the battery state that changes every moment such as the battery temperature Tb and the local SOC θ. ing. The parameter characteristic map 220 is stored in, for example, the storage device that constitutes the monitoring control unit 50.

パラメータ変化推定部66は、例えば最小自乗法を用いて(72)式に示された直流抵抗変化率grの推定値gr’を算出する。ここでは最小自乗法として後述の(74)〜(78)式による忘却要素付きの逐次最小自乗法を例示するが、他の逐次最小自乗法や、一括最小自乗法等の他方式の最小自乗法を適用してもよい。以下に、忘却係数付きの逐次最小自乗法について説明する。   The parameter change estimation unit 66 calculates the estimated value gr ′ of the DC resistance change rate gr shown in the equation (72) using, for example, the least square method. Here, as the least square method, a sequential least square method with a forgetting element according to the following equations (74) to (78) is illustrated, but other sequential least square methods, other methods such as the least square method, and other methods such as the least square method are used. May be applied. Below, the sequential least square method with a forgetting factor is demonstrated.

逐次最小自乗法によれば、下記の(73)式で示す線形回帰モデルで表されるシステムにおいて、(73)式中のパラメータΘは、(74)〜(76)式で示される時間更新式を、(77),(78)式の初期条件により逐次演算することによって推定される。各式においてパラメータΘの推定値は、Θ’で示されている。   According to the sequential least square method, in the system represented by the linear regression model represented by the following equation (73), the parameter Θ in the equation (73) is a time update equation represented by the equations (74) to (76). Is estimated by sequentially calculating according to the initial conditions of equations (77) and (78). In each equation, the estimated value of the parameter Θ is indicated by Θ ′.

(74),(76)式においてλは忘却係数であり、通常はλ<1.0である。また、Pは共分散行列であり、(78)式の初期値P(0)は単位行列Iの対角要素に定数γを乗じた行列とし、γには通常は102〜103程度の大きな値を用いる。パラメータΘ’の初期値Θ’0は通常はゼロベクトルとされる。 In equations (74) and (76), λ is a forgetting factor, and usually λ <1.0. Further, P is a covariance matrix, and the initial value P (0) of the equation (78) is a matrix obtained by multiplying the diagonal element of the unit matrix I by a constant γ, and γ is usually about 10 2 to 10 3 . Use a large value. The initial value Θ ′ 0 of the parameter Θ ′ is normally a zero vector.

このような忘却要素付きの逐次最小自乗法によれば、直流抵抗Raの変化率grを以下のようにして推定可能である。 According to the recursive least square method with forgetting factor, which is the change rate gr of DC resistance R a can be estimated as follows.

新品状態から経年変化(劣化)した二次電池の直流抵抗Raは、(72)式の定義により、Ra=gr・Ranと表せるので、これを(71)式に代入し、さらに(73)式の形に書き直すと、電池モデル式に基づく線形回帰モデル式として(79)式が得られる。 The DC resistance R a of the secondary battery that has been aged (deteriorated) from the new state can be expressed as R a = gr · R an according to the definition of the equation (72). When rewritten in the form of the equation (73), the equation (79) is obtained as a linear regression model equation based on the battery model equation.

(79)式において、左辺の開放電圧U(θ)は、開放電圧推定部110(図12参照)により今回推定された値を用い、Vには測定された電池電圧Vbを用いることにより、Yを計算することができる。また、右辺については、電池状態量推定部52が今回の演算周期での電池温度Tbおよび局所的SOCθを引数としてパラメータ特性マップ220を参照することにより、直流抵抗の初期状態時パラメータ値Ranが求められる。また、電池電流密度Iとしては、電流測定部30による現在の測定値Ibから算出した単位極板面積当たりの電流値を代入することにより、Zを計算することができる。 In equation (79), the open-circuit voltage U (θ) on the left side is the value estimated this time by the open-circuit voltage estimation unit 110 (see FIG. 12), and the measured battery voltage Vb is used for V. Can be calculated. For the right side, the battery state quantity estimating unit 52 refers to the parameter characteristic map 220 with the battery temperature Tb and the local SOC θ in the current calculation cycle as arguments, so that the parameter value R an in the initial state of the DC resistance is Desired. As the battery current density I, Z can be calculated by substituting the current value per unit electrode plate area calculated from the current measurement value Ib by the current measurement unit 30.

このように演算したYおよびZを用いて(74)〜(78)式の忘却要素付き逐次最小自乗法により、推定パラメータΘとして、直流抵抗Raの変化率grを逐次推定することが可能となる。 Such using computed Y and Z in (74) - (78) below the forgetting factor with recursive least square method, as the estimated parameter theta, can be sequentially estimated change rate gr of DC resistance R a Become.

パラメータ変化推定部66は、推定された直流抵抗変化率gr’を電池モデル式に反映することによって、直流抵抗Raの経年変化を反映する。具体的には、電池モデルにおいて、Raに代えて、gr’×Ranを適用する。これにより、電池モデルによる推定精度を向上させることができる。 Parameter change estimation unit 66, by reflecting in the battery model equation the estimated DC resistance change rate gr ', to reflect the aging of the DC resistance R a. Specifically, in the battery model, gr ′ × R an is applied instead of R a . Thereby, the estimation precision by a battery model can be improved.

なお、推定した直流抵抗変化率gr’に代えて、当該推定値gr’を所定の時定数で平滑処理を行なったものを電池モデル式に反映させるようにしてもよい。これによれば、外乱等による直流抵抗変化率の変動を除去することができる。   Instead of the estimated DC resistance change rate gr ', the estimated value gr' that has been smoothed with a predetermined time constant may be reflected in the battery model equation. According to this, it is possible to eliminate fluctuations in the DC resistance change rate due to disturbance or the like.

上記のように直流抵抗Raは電池温度および活物質表面の局所SOCに依存し、特に電池温度の変化に対して大きく変化するため、仮に直流抵抗Raを直接推定しようとすると、温度変化に対して十分な追従速度が必要となる。推定におけるパラメータ追従速度を大きくすると、ノイズや外乱等によるパラメータ変動が大きくなる傾向があり、パラメータを安定的に推定できないという問題点がある。 The DC resistance R a as described above depending on local SOC of the battery temperature and the surface of the active material, in particular for changes greatly with respect to changes in battery temperature, if you try to directly estimate the DC resistance R a, the temperature change On the other hand, sufficient tracking speed is required. When the parameter tracking speed in estimation is increased, there is a tendency that parameter fluctuations due to noise, disturbance, etc. tend to increase, and parameters cannot be estimated stably.

しかしながら、上記手法によれば、直流抵抗Raの初期状態におけるパラメータ値Ranに対する変化率のみを逐次推定すればよい。経時変化により生じる初期状態に対するパラメータ値の変化速度は、一般に電池状態変化によるパラメータ値変化に比較して非常に遅い(時定数が大きい)ため、推定における追従速度を大きくする必要がない。したがって、上記問題点を回避することができる。 However, according to the above method, only the rate of change of the DC resistance R a with respect to the parameter value R an in the initial state needs to be estimated sequentially. The change speed of the parameter value with respect to the initial state caused by the change with time is generally very slow (the time constant is large) compared to the change of the parameter value due to the battery state change. Therefore, the above problem can be avoided.

また、忘却係数付きの逐次最小自乗法によれば、過去に測定された電池データよりも現在の電池データの方に、より重み付けを大きくしてパラメータ推定を実行する。このため、経時変化による直流抵抗変化率が仮に電池温度や充電率等に依存して変化したとしても、この変化に追従して各電池状態に応じた直流抵抗変化率を正確に推定することが可能となる。   In addition, according to the sequential least square method with a forgetting factor, parameter estimation is performed with greater weighting on current battery data than on battery data measured in the past. For this reason, even if the DC resistance change rate due to changes over time changes depending on the battery temperature, the charging rate, etc., it is possible to accurately estimate the DC resistance change rate according to each battery state following this change. It becomes possible.

なお、上記では(73)式に基づいて反応抵抗と純抵抗とを一体に取扱った直流抵抗Raを用いたが、反応抵抗と純抵抗とを分離してそれぞれ独立にパラメータ推定(変化率推定)を行うことも可能である。また、電気二重層キャパシタを通過するキャパシタ電流成分を考慮してパラメータ推定(変化率推定)を行うことも可能である。 In the above (73) the reaction resistance and the pure resistance and the was used DC resistance R a which dealt integrally parameter estimation (change ratio estimating independently by separating the reaction resistance and the pure resistance based on the formula ) Is also possible. It is also possible to perform parameter estimation (change rate estimation) in consideration of the capacitor current component passing through the electric double layer capacitor.

ところで、二次電池が劣化すると、反応抵抗や純抵抗といった直流抵抗が増加するだけでなく、活物質内の反応関与物質の拡散速度が低下して(すなわち拡散係数が低下して)、いわゆる拡散抵抗も増加する。拡散抵抗の増加は、特に大電流での充放電を継続する場合に電池性能や電流−電圧特性に大きな影響を及ぼす。このため、大電流充放電による電池抵抗上昇度合いを精度よく推定するためには、拡散抵抗の変化、すなわち活物質における拡散係数の変化を推定することが、より好ましい。   By the way, when the secondary battery deteriorates, not only the direct current resistance such as reaction resistance and pure resistance increases, but also the diffusion rate of the reaction participating substance in the active material decreases (that is, the diffusion coefficient decreases), so-called diffusion. Resistance also increases. The increase in diffusion resistance has a great effect on battery performance and current-voltage characteristics, particularly when charging / discharging with a large current is continued. For this reason, in order to accurately estimate the degree of increase in battery resistance due to large current charge / discharge, it is more preferable to estimate the change in diffusion resistance, that is, the change in the diffusion coefficient in the active material.

但し、拡散係数Dsは、上記の直流抵抗Raとは異なり、(71),(79)式のような線形モデル式を構成することができない。なぜなら、拡散係数Dsは、(M2a)式等の拡散方程式中のパラメータ値だからである。このため、最小自乗法を適用することができない。 However, the diffusion coefficient D s, unlike the above DC resistance R a, (71), it is impossible to constitute a linear model expression such as (79) below. This is because the diffusion coefficient D s is a parameter value in a diffusion equation such as (M2a). For this reason, the method of least squares cannot be applied.

図23に、パラメータ変化推定部66による拡散係数の推定学習構成を説明するブロック図を例示する。なお、図21では電池抵抗上昇度合い推定部54等の図示は省略している。図23の例では、パラメータ変化推定部66は、推定用データ保存部260、評価関数演算範囲判定部270と、変化率演算部300とを含んでいる。また、変化率演算部300は、推定処理制御部310と、電池モデル演算部320と、評価関数演算部330とを含んでいる。   FIG. 23 illustrates a block diagram for explaining a diffusion coefficient estimation learning configuration by the parameter change estimation unit 66. In FIG. 21, the battery resistance increase degree estimation unit 54 and the like are not shown. In the example of FIG. 23, the parameter change estimation unit 66 includes an estimation data storage unit 260, an evaluation function calculation range determination unit 270, and a change rate calculation unit 300. The change rate calculation unit 300 includes an estimation process control unit 310, a battery model calculation unit 320, and an evaluation function calculation unit 330.

拡散係数Dsについても、下記(91)式に従って、初期状態パラメータ値に対する変化率としてgdを定義する。なお、Dsnは、拡散係数Dsの初期状態パラメータ値である For the diffusion coefficient D s , gd is defined as the rate of change with respect to the initial state parameter value according to the following equation (91). D sn is an initial state parameter value of the diffusion coefficient D s.

推定用データ保存部260は、所定のデータ保存開始条件の成立に応答してデータ保存期間を開始し、データ保存期間中は測定部30,32,34によって順次測定される電池データIb,Vb,Tbを時系列的なデータ配列として保存する。また、推定用データ保存部260は、所定のデータ保存終了条件の成立に応答してデータ保存期間を終了させる。データ保存期間では、各タイミングにおいて、当該タイミングにおける電池電流Ib、電池電圧Vb、電池温度Tbが互いに対応付けられて記憶される。   The estimation data storage unit 260 starts a data storage period in response to the establishment of a predetermined data storage start condition, and battery data Ib, Vb, which are sequentially measured by the measurement units 30, 32, and 34 during the data storage period. Tb is stored as a time-series data array. Further, the estimation data storage unit 260 ends the data storage period in response to establishment of a predetermined data storage end condition. In the data storage period, at each timing, the battery current Ib, the battery voltage Vb, and the battery temperature Tb at the timing are stored in association with each other.

評価関数演算範囲判定部270は、推定用データ保存部260によって電池データが記憶されるデータ保存期間に、拡散抵抗の影響が電池電圧に大きく表れる時間的範囲が含まれているかどうかを判定する。そして、そのような範囲が含まれている場合には、評価関数演算範囲判定部270は、当該範囲の電池データを用いて評価関数の演算を行ない、拡散係数変化率の推定を行なうように評価関数演算範囲として決定する。   The evaluation function calculation range determination unit 270 determines whether or not the data storage period in which the battery data is stored by the estimation data storage unit 260 includes a time range in which the influence of the diffusion resistance appears largely in the battery voltage. If such a range is included, the evaluation function calculation range determination unit 270 performs an evaluation function calculation using the battery data in the range, and evaluates the diffusion coefficient change rate. Determine as the function operation range.

図24に概略的に示されるように、推定用データ保存部260は、時刻t0から始まった電池緩和状態が終了した後、電流が流れ始める時刻t1から一定時間だけ前からデータ保存期間を開始させる。すなわち、電池データ保存期間は、二次電池10が緩和した状態が所定時間以上継続した後に、緩和状態が終了したと判断される場合に開始する。これにより、測定された電池電流Ibの絶対値が所定値以下であり、かつ、電池モデルによって求められた活物質モデル内の各分割区間における最大リチウム濃度差の絶対値が所定値以下である状態が、所定時間以上継続した場合に、電池が緩和した状態であると判定することができる。また、緩和状態の終了は、上記の電流絶対値条件または最大リチウム濃度差条件が不成立になった場合に判定することができる。   As schematically shown in FIG. 24, the estimation data storage unit 260 starts the data storage period from a time t1 when current starts to flow after the battery relaxation state started from time t0 ends. . That is, the battery data storage period starts when it is determined that the relaxation state has ended after the secondary battery 10 has been relaxed for a predetermined time or longer. As a result, the absolute value of the measured battery current Ib is not more than a predetermined value, and the absolute value of the maximum lithium concentration difference in each divided section in the active material model determined by the battery model is not more than the predetermined value. However, when the battery continues for a predetermined time or more, it can be determined that the battery is in a relaxed state. The end of the relaxed state can be determined when the current absolute value condition or the maximum lithium concentration difference condition is not satisfied.

そして、評価関数演算範囲判定部270は、データ保存期間中には上述した評価関数演算範囲が存在するかどうかの判定を逐次行ない、例えば時刻t2〜t3の範囲を評価関数演算範囲T1cとして決定する。このとき、評価関数演算範囲の開始条件は、例えば、電池モデル式によって求められた活物質界面のリチウム濃度と、活物質内の平均リチウム濃度との差が所定値以上であり、かつ、電池電流が充電状態(Ib<0)となったことにより成立する。他方、評価関数演算範囲の終了条件は、例えば、上記開始条件の後、二次電池の充放電が大きく変化することなく二次電池内部でのリチウム拡散が進行した状態、すなわち拡散係数の推定に適した状態が所定時間以上継続した場合に成立する。   Then, the evaluation function calculation range determination unit 270 sequentially determines whether or not the above-described evaluation function calculation range exists during the data storage period, for example, determines the range from time t2 to t3 as the evaluation function calculation range T1c. . At this time, the start condition of the evaluation function calculation range is, for example, that the difference between the lithium concentration at the active material interface obtained by the battery model equation and the average lithium concentration in the active material is a predetermined value or more, and the battery current Is established when the charging state is reached (Ib <0). On the other hand, the end condition of the evaluation function calculation range is, for example, an estimation of the diffusion coefficient after the start condition, in which the lithium diffusion inside the secondary battery has progressed without significant change in the charge / discharge of the secondary battery. It is established when a suitable state continues for a predetermined time or more.

そして、データ保存期間の開始から一定時間が経過した時刻t4では、データ保存期間が終了されて、変化率演算部300による拡散係数変化率の推定処理が開始される。   Then, at a time t4 when a certain time has elapsed from the start of the data storage period, the data storage period ends, and the diffusion coefficient change rate estimation processing by the change rate calculation unit 300 is started.

変化率演算部300において、電池モデル演算部320は、推定処理制御部310によって設定された拡散係数変化率の候補値を用いて拡散係数Dsを設定し、かつ、データ保存期間の電池温度Tbおよび電池電圧Vbのデータ配列を用いて電池モデル式に従った演算を行うことにより、当該データ保存期間に対応した単位極板面積当たりの電池電流の推定値を算出する。なお、かかる電池電流の推定処理を、電池状態量推定部52の機能を利用して行ってもよい。 In the change rate calculation unit 300, the battery model calculation unit 320 sets the diffusion coefficient D s using the diffusion coefficient change rate candidate value set by the estimation processing control unit 310, and the battery temperature Tb during the data storage period. Then, an estimated value of the battery current per unit electrode plate area corresponding to the data storage period is calculated by performing an operation according to the battery model equation using the data array of the battery voltage Vb. The battery current estimation process may be performed using the function of the battery state quantity estimation unit 52.

変化率演算部300の評価関数演算部330は、評価関数演算範囲判定部270により決定された評価関数演算範囲に対応した期間内における、各時点での電池モデル演算部320による単位極板面積当たりの電池電流の推定値と、データ配列中の電池電流(測定値)Ibを両面極板面積で割ることにより算出した単位極板直積当たりの電池電流(測定値)との誤差に基づいて評価関数Jを算出する。   The evaluation function calculation unit 330 of the rate of change calculation unit 300 per unit plate area by the battery model calculation unit 320 at each time point within a period corresponding to the evaluation function calculation range determined by the evaluation function calculation range determination unit 270. An evaluation function based on an error between an estimated value of the battery current and a battery current (measured value) per unit electrode plate calculated by dividing the battery current (measured value) Ib in the data array by the double-sided electrode plate area J is calculated.

評価関数Jは、次の(92)式で与えられる。なお、(92)式において、Ite(n)は電池モデルに従って推定された単位極板面積当たりの電池電流であり、Ib(n)は電池電流の測定値Ibを単位極板面積に換算したものである。また、NsおよびNeは、推定用データ保存部260に記憶される時系列の電池データDT(n)の配列番号(時間軸)を示す変数nについて、評価関数演算範囲の開始タイミングおよび終了タイミングにそれぞれ対応する。 The evaluation function J is given by the following equation (92). In Equation (92), I te (n) is a battery current per unit electrode plate area estimated according to the battery model, and Ib (n) is a measured value Ib of the battery current converted into a unit electrode plate area. Is. Ns and Ne are the start timing and end timing of the evaluation function calculation range for the variable n indicating the array number (time axis) of the time-series battery data DT (n) stored in the estimation data storage unit 260. Each corresponds.

変化率演算部300の推定処理制御部310は、拡散係数変化率の候補値(例えば下記のGSM法により設定される)を複数個切換えて設定するとともに、各候補値に対する評価関数Jの比較に基づいて、拡散係数変化率の推定値gd’を算出する。推定処理制御部310による拡散係数変化率推定演算処理は、保存した電池データを用いて、評価関数演算範囲T1cにおいて評価関数Jが最小となる変化率を、所定の繰り返し回数探索することにより設定する。評価関数Jが最小となるような拡散係数変化率の探索手法として、例えば公知のGSM法(黄金分割法)を用いることが可能である。   The estimation processing control unit 310 of the change rate calculation unit 300 sets a plurality of candidate values (for example, set by the following GSM method) of the diffusion coefficient change rate, and compares the evaluation function J for each candidate value. Based on this, an estimated value gd ′ of the diffusion coefficient change rate is calculated. The diffusion coefficient change rate estimation calculation processing by the estimation processing control unit 310 is set by searching for a change rate at which the evaluation function J is minimum in the evaluation function calculation range T1c by using the stored battery data, by searching a predetermined number of times. . As a technique for searching for a diffusion coefficient change rate that minimizes the evaluation function J, for example, a known GSM method (golden section method) can be used.

推定処理制御部310は、所定の探索繰り返し回数が終了したときの探索範囲の中心値を拡散係数変化率推定値gd’として採用する。そして、推定処理制御部310は、算出された拡散係数変率の推定値gd’を電池モデル式に反映する。具体的には、電池モデルにおいて、Dsに代えて、gd’×Dsnを適用する。これにより、電池モデルによる推定精度を向上させることができる。 The estimation processing control unit 310 employs the center value of the search range when the predetermined number of search repetitions is completed as the diffusion coefficient change rate estimated value gd ′. Then, the estimation process control unit 310 reflects the calculated diffusion coefficient variation estimated value gd ′ in the battery model equation. Specifically, in the battery model, gd ′ × D sn is applied instead of D s . Thereby, the estimation precision by a battery model can be improved.

なお、パラメータ変化推定部66は、上記例示の構成に限定されるものではなく、各種の手法に応じて構成してもよい。   The parameter change estimation unit 66 is not limited to the configuration illustrated above, and may be configured according to various methods.

ここで、例えば電池モデルにおいて電池温度を予め一定値(定数)に設定することにより、電池温度の測定値を用いずに電池電流を推定するように、電池状態量推定部52を構成してもよい。このとき、電池温度についての当該一定値は、例えば使用状態すなわち充放電を行っている状態の二次電池10について予め測定した電池温度の代表的な値(平均値等)に設定することが可能である。これに対し、上記のように電池温度の測定値も利用することによって、高い推定精度が得られる。   Here, for example, by setting the battery temperature to a constant value (constant) in advance in the battery model, the battery state quantity estimating unit 52 may be configured to estimate the battery current without using the measured value of the battery temperature. Good. At this time, the constant value for the battery temperature can be set to a representative value (average value or the like) of the battery temperature measured in advance for the secondary battery 10 in a use state, that is, in a state of being charged / discharged. It is. On the other hand, high estimation accuracy is obtained by using the measured value of the battery temperature as described above.

また、上記のように、電圧−電流関係モデル式(M1a)〜(M1d)に従い、電池電流Ibの測定値に基づいて電池電圧V(t)を推定するように電池状態量推定部52を構成することも可能である。この場合も、電池電圧の推定値と測定値とに基づいて、例えば電池電圧の推定値と測定値との差や比等を算出することによって、大電流域における電池抵抗の上昇度合いを推定することができ、このようにして推定された上昇度合いは、電池電流に基づく上記の電池抵抗上昇度合いに代えてまたは加えて、監視制御部50で利用することができる。   Further, as described above, the battery state quantity estimating unit 52 is configured to estimate the battery voltage V (t) based on the measured value of the battery current Ib according to the voltage-current relationship model equations (M1a) to (M1d). It is also possible to do. Also in this case, the degree of increase in battery resistance in a large current region is estimated by calculating, for example, the difference or ratio between the estimated value and measured value of the battery voltage based on the estimated value and measured value of the battery voltage. The degree of increase estimated in this way can be used by the monitoring control unit 50 instead of or in addition to the above-described degree of increase in battery resistance based on the battery current.

二次電池を説明する概略構成図である。 It is a schematic block diagram explaining a secondary battery. 二次電池における電解液塩濃度分布を説明する概念図である。 It is a conceptual diagram explaining the electrolyte solution salt density | concentration distribution in a secondary battery. 電解液塩濃度と反応抵抗との関係を示す特性図である。 It is a characteristic view which shows the relationship between electrolyte solution salt concentration and reaction resistance. 電極内電解液塩濃度の低下を例示する概念図である。 It is a conceptual diagram which illustrates the fall of the electrolyte solution salt density | concentration in an electrode. 実施の形態に係る二次電池監視装置を説明する概略構成図である。 It is a schematic block diagram explaining the secondary battery monitoring apparatus which concerns on embodiment. 電池モデル式で用いられる変数等の一覧を示す図である。 It is a figure which shows the list of variables etc. which are used with a battery model type | formula. 電池モデルを説明する概念図である。 It is a conceptual diagram explaining a battery model. 極座標で示された活物質モデルを示す概念図である。 It is a conceptual diagram which shows the active material model shown by the polar coordinate. 二次電池の端子電圧と各種平均電位との関係を示す概念図である。 It is a conceptual diagram which shows the relationship between the terminal voltage of a secondary battery, and various average electric potential. 拡散係数の温度依存性を説明する概念図であるIt is a conceptual diagram explaining the temperature dependence of a diffusion coefficient. 開放電圧と局所SOCとの関係を示す概念図である。 It is a conceptual diagram which shows the relationship between an open circuit voltage and local SOC. 実施の形態に係る電池状態量推定部を説明する概略構成図である。 It is a schematic block diagram explaining the battery state quantity estimation part which concerns on embodiment. 実施の形態に係る電池状態量推定部での処理を例示するフローチャートである。 It is a flowchart which illustrates the process in the battery state quantity estimation part which concerns on embodiment. 実施の形態について電極間電解液塩濃度と電流推定誤差との間の相関図である。 It is a correlation diagram between electrolyte solution salt concentration and electric current estimation error about embodiment. 実施の形態について電極間電解液塩濃度と抵抗上昇率との間の相関図である。 It is a correlation figure between electrolyte solution salt concentration and resistance increase rate about embodiment. 実施の形態に係る監視制御部を説明する概略構成図である。 It is a schematic block diagram explaining the monitoring control part which concerns on embodiment. 実施の形態に係る監視制御部での処理を例示するフローチャートである。 It is a flowchart which illustrates the process in the monitoring control part which concerns on embodiment. 実施の形態について電池抵抗上昇度合いの推定結果を示す図である。 It is a figure which shows the estimation result of a battery resistance raise degree about embodiment. 実施の形態について二次電池の状態を測定した結果を示す図である。 It is a figure which shows the result of having measured the state of the secondary battery about embodiment. 実施の形態について二次電池の状態を測定した結果を示す図である。 It is a figure which shows the result of having measured the state of the secondary battery about embodiment. 実施の形態に係るパラメータ変化推定部による推定学習を説明するブロック図である。 It is a block diagram explaining the estimation learning by the parameter change estimation part which concerns on embodiment. 局所SOCおよび電池温度の変化に対する直流抵抗の変化特性を示す概念図である。 It is a conceptual diagram which shows the change characteristic of direct-current resistance with respect to the change of local SOC and battery temperature. 実施の形態に係るパラメータ変化推定部による他の推定学習を説明するブロック図である。 It is a block diagram explaining the other estimation learning by the parameter change estimation part which concerns on embodiment. 実施の形態についてデータ保存期間および評価関数演算範囲を説明する概念図である。 It is a conceptual diagram explaining the data preservation | save period and evaluation function calculation range about embodiment.

符号の説明Explanation of symbols

5 二次電池の監視装置、10 二次電池、12 負極、15 正極、30 電流測定部、32 電圧測定部、34 温度測定部、50 監視制御部、52 電池状態量推定部、54 電池抵抗上昇度合い推定部、56 通知部、58 使用制限部、60 充電率制御部、62 診断部、64 強制的充放電制御部、66 パラメータ変化推定部、V,Vb 電池電圧、I,Ib,Im,Ir 電池電流密度(電池電流)、T,Tb 電池温度、Dg 電池抵抗上昇度合い、Δce 電解液塩濃度差。 5 Secondary battery monitoring device, 10 Secondary battery, 12 Negative electrode, 15 Positive electrode, 30 Current measurement unit, 32 Voltage measurement unit, 34 Temperature measurement unit, 50 Monitoring control unit, 52 Battery state quantity estimation unit, 54 Battery resistance increase Degree estimation unit, 56 notification unit, 58 use restriction unit, 60 charge rate control unit, 62 diagnosis unit, 64 forced charge / discharge control unit, 66 parameter change estimation unit, V, Vb battery voltage, I, Ib, Im, Ir battery current density (battery current), T, Tb battery temperature, D g battery resistance increasing degree, .DELTA.c e electrolyte salt concentration difference.

Claims (18)

  1. 二次電池の電池電圧と電池電流を電池状態量として測定する測定部と、
    バトラー・ボルマーの式において、負極中および正極中における電気化学反応を所定のものと仮定した第1モデル式、活物質拡散モデル式、及び前記電池状態量の一方に基づき、推定電池状態量を取得する取得部と、
    前記推定電池状態量と前記電池状態量との間の電流推定誤差、及び電極間電解液濃度差に基づいて前記二次電池の電池抵抗の上昇度合いを推定する電池抵抗上昇度合い推定部と、
    を備える二次電池の監視装置。
    A measurement unit that measures a battery voltage and battery current of the secondary battery as batteries state quantity,
    In the Butler-Volmer equation, the estimated battery state quantity is obtained based on one of the first model expression, the active material diffusion model expression, and the battery state quantity assuming that the electrochemical reaction in the negative electrode and the positive electrode is predetermined. An acquisition unit to In the Butler-Volmer equation, the estimated battery state quantity is obtained based on one of the first model expression, the active material diffusion model expression, and the battery state quantity assuming that the electrochemical reaction in the negative electrode and the positive electrode is predetermined . An acquisition unit to
    A battery resistance increase degree estimation unit that estimates a degree of increase in battery resistance of the secondary battery based on a current estimation error between the estimated battery state quantity and the battery state quantity, and an interelectrode electrolyte concentration difference ; A battery resistance increase degree estimation unit that estimates a degree of increase in battery resistance of the secondary battery based on a current estimation error between the estimated battery state quantity and the battery state quantity, and an interelectrode electrolyte concentration difference ;
    Monitoring device for a secondary battery Ru comprising a. Monitoring device for a secondary battery Ru comprising a.
  2. 前記取得部は、
    前記第1モデル式において、前記負極と前記正極の活物質モデルを別々とせずにまとめて単独の素子として簡略化した第2モデル式に基づいて、推定電池状態量を取得する取得部である、請求項1記載の二次電池の監視装置。 In the first model formula, it is an acquisition unit that acquires an estimated battery state amount based on the second model formula that simplifies the active material models of the negative electrode and the positive electrode as a single element without separating them. The secondary battery monitoring device according to claim 1 . The acquisition unit The acquisition unit
    In the first model formula, the acquisition unit obtains an estimated battery state quantity based on a second model formula obtained by simplifying the active material models of the negative electrode and the positive electrode together as a single element without being separated. The monitoring apparatus for a secondary battery according to claim 1 . In the first model formula, the acquisition unit obtains an estimated battery state quantity based on a second model formula obtained by simplifying the active material models of the negative electrode and the positive electrode together as a single element without being separated. The monitoring apparatus for a secondary battery according to claim 1 .
  3. 前記電池抵抗上昇度合い推定部は、
    前記一方の電池状態量についての推定値と測定値とともに、前記二次電池について推定された電極間の電解液塩濃度差と、前記電池電流の測定値の移動平均値との少なくとも一方に基づいて、前記電池抵抗上昇度合いを推定する請求項1または2記載の二次電池の監視装置。 Based on at least one of the estimated difference in electrolyte salt concentration between the electrodes for the secondary battery and the moving average of the measured values ​​of the battery current, as well as the estimated and measured values ​​for one of the battery states. The secondary battery monitoring device according to claim 1 or 2 , wherein the degree of increase in battery resistance is estimated. The battery resistance increase degree estimation unit The battery resistance increase degree estimation unit
    Based on the estimated value and the measured value for the one battery state quantity, and at least one of the electrolyte salt concentration difference between the electrodes estimated for the secondary battery and the moving average value of the measured value of the battery current The secondary battery monitoring device according to claim 1 , wherein the battery resistance increase degree is estimated. Based on the estimated value and the measured value for the one battery state quantity, and at least one of the electrolyte salt concentration difference between the electrodes estimated for the secondary battery and the moving average value of the measured value of the battery current The secondary battery monitoring device according to claim 1 , which the battery resistance increase degree is estimated.
  4. 前記二次電池の電池温度を測定する温度測定部をさらに備え、
    前記電池状態量推定部は、前記他方の電池状態量の測定値とともに前記電池温度の測定値に基づき、前記一方の電池状態量を前記電池モデルに従って推定する請求項1から3のいずれか1項に記載の二次電池の監視装置。
    A temperature measuring unit for measuring a battery temperature of the secondary battery;
    The said battery state quantity estimation part estimates one said battery state quantity according to the said battery model based on the measured value of the said battery temperature with the measured value of said other battery state quantity, The any one of Claim 1 to 3 The monitoring apparatus of the secondary battery as described in .
  5. 記電池抵抗上昇度合い推定部は、前記電池抵抗上昇度合いとして、
    前記一方の電池状態量についての測定値と推定値との比もしくは差、前記比もしくは差の増加率であって前記二次電池について推定された電極間の電解液塩濃度差に対する増加率、または、
    前記比もしくは差と前記推定された電解液塩濃度差との一次近似式において前記推定電解液塩濃度差を0とした場合の値を算出する請求項1ないし4のいずれか1項に記載の二次電池の監視装置。
    Before SL battery resistance increasing degree estimating unit, as the battery resistance increasing degree,
    The ratio or difference between the measured value and the estimated value for the one battery state quantity, the rate of increase of the ratio or difference and the rate of increase relative to the difference in electrolyte salt concentration between the electrodes estimated for the secondary battery, or , The ratio or difference between the measured value and the estimated value for the one battery state quantity, the rate of increase of the ratio or difference and the rate of increase relative to the difference in electrolyte salt concentration between the electrodes estimated for the secondary battery, or,
    5. The value according to claim 1, wherein a value when the estimated electrolyte salt concentration difference is set to 0 in a first-order approximation formula of the ratio or difference and the estimated electrolyte salt concentration difference is calculated. Secondary battery monitoring device. 5. The value according to claim 1, which a value when the estimated electrolyte salt concentration difference is set to 0 in a first-order approximation formula of the ratio or difference and the estimated electrolyte salt concentration difference is calculated. Secondary battery monitoring device.
  6. 記電池抵抗上昇度合いが予め設定された通知設定値よりも大きい場合に、前記電池抵抗の上昇を通知する通知部をさらに備える請求項1ないし5のいずれか1項に記載の二次電池の監視装置。 If prior SL cell resistance increasing degree is greater than a preset notification settings, the secondary battery according to any one of claims 1 to 5, further comprising a notification unit for notifying an increase in the battery resistance Monitoring device.
  7. 記電池抵抗上昇度合いが予め設定された使用制限設定値よりも大きい場合に、充電電力量、放電電力量、充電電流量、放電電流量、電池電圧上限値、電池電圧下限値の少なくともいずれか1を制限する使用制限部をさらに備える請求項1ないし6のいずれか1項に記載の二次電池の監視装置。 If prior SL cell resistance increasing degree is greater than a preset use limitation set value, the amount of charge power, amount of discharge power, the amount of charge current, discharge current, battery voltage upper limit value, at least one of the battery voltage lower limit The monitoring apparatus for a secondary battery according to claim 1 , further comprising a use restriction unit that restricts 1 .
  8. 前記電池抵抗上昇度合いが予め設定された充電率制御中心調整設定値よりも大きい場合に、充電率制御の制御中心値を調整する充電率制御部をさらに備える請求項1ないし7のいずれか1項に記載の二次電池の監視装置。 The charge rate control part which adjusts the control center value of charge rate control when the said battery resistance raise degree is larger than the preset charge rate control center adjustment set value, The any one of Claim 1 thru | or 7 further provided. The monitoring apparatus of the secondary battery as described in .
  9. 前記二次電池の診断を行う診断部をさらに備え、
    前記診断部は、前記電池抵抗上昇度合いが予め設定された電圧下限値診断禁止設定値よりも大きい場合に、前記電池電圧の下限値についての故障診断を行わない請求項1ないし8のいずれか1項に記載の二次電池の監視装置。 Any one of claims 1 to 8 in which the diagnosis unit does not perform a failure diagnosis for the lower limit value of the battery voltage when the degree of increase in battery resistance is larger than a preset voltage lower limit value diagnosis prohibition setting value. The secondary battery monitoring device described in the section . A diagnostic unit for diagnosing the secondary battery; A diagnostic unit for diagnosing the secondary battery;
    9. The diagnosis unit according to claim 1, wherein the diagnosis unit does not perform a failure diagnosis for the lower limit value of the battery voltage when the battery resistance increase degree is larger than a preset voltage lower limit value diagnosis prohibition set value. The monitoring apparatus of the secondary battery as described in a term . 9. The diagnosis unit according to claim 1, wherein the diagnosis unit does not perform a failure diagnosis for the lower limit value of the battery voltage when the battery resistance increase degree is larger than a preset voltage lower limit value diagnosis prohibition set value. monitoring apparatus of the secondary battery as described in a term .
  10. 前記電池抵抗上昇度合いが予め設定された充放電要求設定値よりも大きい場合に、前記二次電池の強制的な充電または放電を行う強制的充放電制御部をさらに備える請求項1ないし9のいずれか1項に記載の二次電池の監視装置。 The forced charging / discharging control part which performs the forced charge or discharge of the said secondary battery when the said battery resistance raise degree is larger than the preset charging / discharging request | requirement setting value. or the secondary battery monitoring device according to item 1.
  11. 前記二次電池および前記監視装置がエンジン搭載車両に適用される場合、
    前記強制的充放電制御部は、エンジンに対して要求するエンジン要求パワーを増加させて、エンジン出力によって前記強制的な充電または放電を行う請求項10に記載の二次電池の監視装置。 The secondary battery monitoring device according to claim 10, wherein the forced charge / discharge control unit increases the engine required power required for the engine and performs the forced charge or discharge by the engine output . When the secondary battery and the monitoring device are applied to an engine-equipped vehicle, When the secondary battery and the monitoring device are applied to an engine-equipped vehicle,
    The secondary battery monitoring apparatus according to claim 10, wherein the forced charge / discharge control unit increases the required engine power required for the engine and performs the forced charge or discharge according to an engine output . The secondary battery monitoring apparatus according to claim 10, wherein the forced charge / discharge control unit increases the required engine power required for the engine and performs the forced charge or discharge according to an engine output .
  12. 前記二次電池が複数個接続されて電池パックが構成されており、
    前記強制的充放電制御部は、前記電池パックの温度が予め定められた定常状態にない場合には、前記強制的な充電または放電を行わない請求項10または11に記載の二次電池の監視装置。 The monitoring of the secondary battery according to claim 10 or 11, wherein the forced charge / discharge control unit does not perform the forced charge or discharge when the temperature of the battery pack is not in a predetermined steady state. apparatus. A plurality of secondary batteries are connected to form a battery pack, A plurality of secondary batteries are connected to form a battery pack,
    The secondary battery monitoring according to claim 10 or 11, wherein the forced charge / discharge control unit does not perform the forced charge or discharge when the temperature of the battery pack is not in a predetermined steady state. apparatus. The secondary battery monitoring according to claim 10 or 11, wherein the forced charge / discharge control unit does not perform the forced charge or discharge when the temperature of the battery pack is not in a predetermined steady state .
  13. 前記電池モデルにおける抵抗成分のパラメータについて前記二次電池の経時劣化に伴う変化を推定し、当該経時劣化推定結果を前記電池モデルに反映させるパラメータ変化推定部をさらに備える請求項1ないし12のいずれか1項に記載の二次電池の監視装置。 The parameter change estimation unit according to any one of claims 1 to 12, further comprising: a parameter change estimation unit that estimates a change of the resistance component in the battery model due to deterioration with time of the secondary battery and reflects the result of deterioration with time in the battery model. 2. The secondary battery monitoring device according to item 1 .
  14. 前記電池抵抗上昇度合い推定部は、前記パラメータ変化推定部による前記経時劣化推定結果に基づいて前記電池抵抗上昇度合いの推定を実行するか否かを判定する請求項13に記載の二次電池の監視装置。 The secondary battery monitoring according to claim 13, wherein the battery resistance increase degree estimation unit determines whether or not to perform the estimation of the battery resistance increase degree based on the time degradation estimation result by the parameter change estimation unit. apparatus.
  15. 記電池状態量推定部は、前記電池電圧の測定値に基づき前記電池モデルに従って前記電池電流を推定し、
    前記電池抵抗上昇度合い推定部は、前記電池電流の推定値と測定値とに基づいて前記電池抵抗上昇度合いを推定する請求項1ないし14のいずれか1項に記載の二次電池の監視装置。 The secondary battery monitoring device according to any one of claims 1 to 14, wherein the battery resistance increase degree estimation unit estimates the battery resistance increase degree based on the estimated value and the measured value of the battery current . Before SL battery state estimating unit, the battery current is estimated according to the battery model based on the measured value of the battery voltage, Before SL battery state estimating unit, the battery current is estimated according to the battery model based on the measured value of the battery voltage,
    The secondary battery monitoring apparatus according to claim 1, wherein the battery resistance increase degree estimation unit estimates the battery resistance increase degree based on an estimated value and a measured value of the battery current . The secondary battery monitoring apparatus according to claim 1, wherein the battery resistance increase degree estimation unit estimates the battery resistance increase degree based on an estimated value and a measured value of the battery current .
  16. 前記電池状態量推定部は、前記電池電流の測定値に基づき前記電池モデルに従って前記電池電圧を推定し、The battery state quantity estimation unit estimates the battery voltage according to the battery model based on the measured value of the battery current,
    前記電池抵抗上昇度合い推定部は、前記電池電圧の推定値と測定値とに基づいて前記電池抵抗上昇度合いを推定する請求項1ないし14のいずれか1項に記載の二次電池の監視装置。 The secondary battery monitoring apparatus according to claim 1, wherein the battery resistance increase degree estimation unit estimates the battery resistance increase degree based on an estimated value and a measured value of the battery voltage.
  17. 二次電池の電池電圧と電池電流を電池状態量として測定する測定部と、 A measurement unit for measuring the battery voltage and battery current of the secondary battery as a battery state quantity;
    バトラー・ボルマーの式において、負極中および正極中における電気化学反応を所定のものと仮定した第1モデル式、活物質拡散モデル式、及び前記電池状態量の一方に基づき、推定電池状態量を取得する取得部と、In the Butler-Volmer equation, the estimated battery state quantity is obtained based on one of the first model expression, the active material diffusion model expression, and the battery state quantity assuming that the electrochemical reaction in the negative electrode and the positive electrode is predetermined. An acquisition unit to
    前記推定電池状態量が前記電池状態量よりも所定値以上大きい場合に所定の処理を行う処理部と、A processing unit that performs a predetermined process when the estimated battery state quantity is larger than the battery state quantity by a predetermined value;
    を備える二次電池システム。A secondary battery system comprising:
  18. 前記処理部は、電極間電解液濃度差が所定値よりも大きい場合に前記所定の処理を行う請求項17記載の二次電池システム。The secondary battery system according to claim 17, wherein the processing unit performs the predetermined processing when a difference in concentration of electrolyte solution between electrodes is larger than a predetermined value.
JP2008225649A 2008-09-03 2008-09-03 Secondary battery monitoring device and secondary battery system Expired - Fee Related JP5036662B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008225649A JP5036662B2 (en) 2008-09-03 2008-09-03 Secondary battery monitoring device and secondary battery system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008225649A JP5036662B2 (en) 2008-09-03 2008-09-03 Secondary battery monitoring device and secondary battery system

Publications (2)

Publication Number Publication Date
JP2010060406A JP2010060406A (en) 2010-03-18
JP5036662B2 true JP5036662B2 (en) 2012-09-26

Family

ID=42187357

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008225649A Expired - Fee Related JP5036662B2 (en) 2008-09-03 2008-09-03 Secondary battery monitoring device and secondary battery system

Country Status (1)

Country Link
JP (1) JP5036662B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016160667A (en) * 2015-03-03 2016-09-05 日立建機株式会社 Hybrid type construction machine

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012046265A1 (en) * 2010-10-05 2012-04-12 トヨタ自動車株式会社 Estimation device and estimation method
JP5704014B2 (en) * 2011-08-02 2015-04-22 トヨタ自動車株式会社 Secondary battery monitoring device
EP2746796B1 (en) 2011-08-03 2016-07-20 Toyota Jidosha Kabushiki Kaisha Device for estimating state of deterioration of secondary battery and method for estimating state of deterioration
JP5770563B2 (en) * 2011-08-23 2015-08-26 トヨタ自動車株式会社 Vehicle power supply system
JP5687584B2 (en) * 2011-08-23 2015-03-18 トヨタ自動車株式会社 Lithium-ion battery condition measurement device
JP5989320B2 (en) * 2011-10-07 2016-09-07 トヨタ自動車株式会社 Vehicle control apparatus and remaining capacity estimation method
JP2013108937A (en) * 2011-11-24 2013-06-06 Toyota Motor Corp Battery state estimation device
JP5586114B2 (en) * 2011-11-25 2014-09-10 トヨタ自動車株式会社 Battery state estimation device
DE112012005901T5 (en) 2012-02-17 2014-11-13 Toyota Jidosha Kabushiki Kaisha Battery system and deterioration determination method
JP5704108B2 (en) * 2012-03-30 2015-04-22 トヨタ自動車株式会社 Battery system and estimation method
JP5768754B2 (en) * 2012-03-30 2015-08-26 トヨタ自動車株式会社 Estimation apparatus and estimation method
JP5783122B2 (en) 2012-04-11 2015-09-24 トヨタ自動車株式会社 Battery state estimation device
JP5831631B2 (en) * 2012-05-15 2015-12-09 トヨタ自動車株式会社 Battery system and method for determining polarization of secondary battery
WO2014109271A1 (en) * 2013-01-14 2014-07-17 Semiconductor Energy Laboratory Co., Ltd. Electrochemical device
JP2014158414A (en) 2013-01-21 2014-08-28 Semiconductor Energy Lab Co Ltd Vehicle having power storage body
JP2014153269A (en) * 2013-02-12 2014-08-25 Toyota Motor Corp Secondary battery inspection method
TWI619329B (en) * 2013-03-01 2018-03-21 半導體能源研究所股份有限公司 Power storage system
JP5862631B2 (en) 2013-10-08 2016-02-16 トヨタ自動車株式会社 Power storage system
JP6151163B2 (en) * 2013-12-06 2017-06-21 株式会社東芝 Battery state calculation device and battery state calculation method
JP6102734B2 (en) 2013-12-27 2017-03-29 トヨタ自動車株式会社 Battery system
JP6154352B2 (en) 2014-04-21 2017-06-28 トヨタ自動車株式会社 Battery system
KR20160011028A (en) 2014-07-21 2016-01-29 삼성전자주식회사 Method and device to detect abnormal state of battery
JP6213511B2 (en) 2015-03-25 2017-10-18 トヨタ自動車株式会社 Electric vehicle and control method thereof
JP6428743B2 (en) 2016-10-26 2018-11-28 トヨタ自動車株式会社 Automobile
JP2018185260A (en) * 2017-04-27 2018-11-22 トヨタ自動車株式会社 Charge rate estimation method and onboard battery system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11346444A (en) * 1998-06-02 1999-12-14 Toyota Motor Corp Estimating method of battery charged condition
JP2006226788A (en) * 2005-02-16 2006-08-31 Fuji Heavy Ind Ltd Battery management system
JP4802945B2 (en) * 2006-08-31 2011-10-26 トヨタ自動車株式会社 Secondary battery control system and hybrid vehicle equipped with the same
JP2008058278A (en) * 2006-09-04 2008-03-13 Toyota Motor Corp Internal state estimating device for secondary cell, internal state estimation method for secondary cell, program, and recording medium

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016160667A (en) * 2015-03-03 2016-09-05 日立建機株式会社 Hybrid type construction machine
EP3064652A1 (en) 2015-03-03 2016-09-07 Hitachi Construction Machinery Co., Ltd. Hybrid construction machine
KR20160107092A (en) 2015-03-03 2016-09-13 히다찌 겐끼 가부시키가이샤 Hybrid construction machine
US9744962B2 (en) 2015-03-03 2017-08-29 Hitachi Construction Machinery Co., Ltd. Hybrid construction machine

Also Published As

Publication number Publication date
JP2010060406A (en) 2010-03-18

Similar Documents

Publication Publication Date Title
US20200088806A1 (en) Systems and methods for determining battery state of charge
US10300798B2 (en) Determining battery DC impedance
EP2700966B1 (en) Apparatus and method for estimating battery state
DE102011117249B4 (en) A method and apparatus for assessing battery health
US20160011274A1 (en) Battery life estimation method and battery life estimation apparatus
EP1429152B1 (en) Method for predicting the internal resistance of a battery and a monitoring device for batteries
US7012434B2 (en) Method for determining the amount of charge which can be drawn from a storage battery and monitoring device
US7362076B2 (en) Remaining capacity calculating device and method for electric power storage
CN107076801B (en) Electrochemical impedance spectroscopy in battery management systems
US9537325B2 (en) Battery state estimation system, battery control system, battery system, and battery state estimation method
US9037426B2 (en) Systems and methods for determining cell capacity values in a multi-cell battery
JP5656415B2 (en) Secondary battery state determination device and control device
US7630842B2 (en) Secondary battery charge/discharge electricity amount estimation method and device, secondary battery polarization voltage estimation method and device and secondary battery remaining capacity estimation method and device
DE102011117098B4 (en) Method and device for assessing the battery termination status
US9071072B2 (en) Available charging/discharging current calculation method and power supply device
DE10158029B4 (en) Method for calculating the dynamic state of charge in a battery
JP4649101B2 (en) Secondary battery status detection device and status detection method
US8887872B2 (en) Method of determining state of charge of energy storage system
US9187007B2 (en) Online battery capacity estimation
KR100970841B1 (en) Apparatus and Method for estimating battery&#39;s state of health based on battery voltage variation pattern
EP1960800B1 (en) Determination of ir-free voltage in hybrid vehicle applications
US7136762B2 (en) System for calculating remaining capacity of energy storage device
KR101399388B1 (en) Apparatus and Method for estimating the life span of battery
US6362598B2 (en) Method for determining the state of charge and loading capacity of an electrical storage battery
US8521497B2 (en) Systems, methods and computer-readable media for modeling cell performance fade of rechargeable electrochemical devices

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100719

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20111125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120221

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120419

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120605

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120703

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150713

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150713

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees