JP6406468B1 - 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム - Google Patents
蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム Download PDFInfo
- Publication number
- JP6406468B1 JP6406468B1 JP2018062290A JP2018062290A JP6406468B1 JP 6406468 B1 JP6406468 B1 JP 6406468B1 JP 2018062290 A JP2018062290 A JP 2018062290A JP 2018062290 A JP2018062290 A JP 2018062290A JP 6406468 B1 JP6406468 B1 JP 6406468B1
- Authority
- JP
- Japan
- Prior art keywords
- voltage value
- storage amount
- amount
- charge
- soc
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
- G01R31/3842—Arrangements for monitoring battery or accumulator variables, e.g. SoC combining voltage and current measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/42—Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
- H01M10/48—Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Abstract
【課題】蓄電量−電圧値特性がヒステリシスを示す活物質を含む蓄電素子の蓄電量を推定することができる蓄電量推定装置、該蓄電量推定装置を備える蓄電モジュール、蓄電量推定方法、及びコンピュータプログラムを提供する。【解決手段】蓄電量推定装置６は、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、蓄電量−電圧値充電特性及び蓄電量−電圧値放電特性間のヒステリシスが、他の電気化学反応が生じる場合のヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子３の蓄電量を推定する。蓄電量推定装置６は、一の電気化学反応が他の電気化学反応より多く生じる場合に、蓄電量−電圧値放電特性を、蓄電量−電圧値充電特性に向けて移動させて得られる電圧参照蓄電量−電圧値特性に基づき、蓄電量を推定する推定部６２を備える。【選択図】図１５
Description
本発明は、蓄電素子のＳＯＣ（State Of Charge）等の蓄電量を推定する蓄電量推定装
置、該蓄電量推定装置を含む蓄電モジュール、蓄電量推定方法、及びコンピュータプログラムに関する。
置、該蓄電量推定装置を含む蓄電モジュール、蓄電量推定方法、及びコンピュータプログラムに関する。
電気自動車、ハイブリッド車等に用いられる車両用の二次電池や、電力貯蔵装置、太陽光発電システム等に用いられる産業用の二次電池においては、高容量化が求められている。これまで様々な検討と改良が行われてきて、電極構造等の改良のみで更なる高容量化を実現することは困難な傾向にある。その為、現行の材料より高容量である正極材料の開発が進められている。
従来、リチウムイオン二次電池等の非水電解質二次電池用の正極活物質として、α−ＮａＦｅＯ２型結晶構造を有するリチウム遷移金属複合酸化物が検討され、ＬｉＣｏＯ２を用いた非水電解質二次電池が広く実用化されていた。ＬｉＣｏＯ２の放電容量は１２０〜１６０ｍＡｈ／ｇ程度であった。
リチウム遷移金属複合酸化物をＬｉＭｅＯ２（Ｍｅは遷移金属）で表したとき、ＭｅとしてＭｎを用いることが望まれてきた。ＭｅとしてＭｎを含有させた場合、Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５を超える場合には、充電をするとスピネル型へと構造変化が起こり、結晶構造が維持できない為、充放電サイクル性能が著しく劣る。
Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５以下であり、Ｍｅに対するＬｉのモル比Ｌｉ／Ｍｅが略１であるＬｉＭｅＯ２型活物質が種々提案され、実用化されている。リチウム遷移金属複合酸化物であるＬｉＮｉ１／２Ｍｎ１／２Ｏ２及びＬｉＮｉ１／３Ｃｏ１／３Ｍｎ１／３Ｏ２等を含有する正極活物質は１５０〜１８０ｍＡｈ／ｇの放電容量を有する。
リチウム遷移金属複合酸化物をＬｉＭｅＯ２（Ｍｅは遷移金属）で表したとき、ＭｅとしてＭｎを用いることが望まれてきた。ＭｅとしてＭｎを含有させた場合、Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５を超える場合には、充電をするとスピネル型へと構造変化が起こり、結晶構造が維持できない為、充放電サイクル性能が著しく劣る。
Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５以下であり、Ｍｅに対するＬｉのモル比Ｌｉ／Ｍｅが略１であるＬｉＭｅＯ２型活物質が種々提案され、実用化されている。リチウム遷移金属複合酸化物であるＬｉＮｉ１／２Ｍｎ１／２Ｏ２及びＬｉＮｉ１／３Ｃｏ１／３Ｍｎ１／３Ｏ２等を含有する正極活物質は１５０〜１８０ｍＡｈ／ｇの放電容量を有する。
ＬｉＭｅＯ２型活物質に対し、Ｍｅ中のＭｎのモル比Ｍｎ／Ｍｅが０．５を超え、遷移金属（Ｍｅ）の比率に対するＬｉの組成比率Ｌｉ／Ｍｅが１より大きいリチウム遷移金属複合酸化物を含む、いわゆるリチウム過剰型活物質も知られている。
上述の高容量の正極材料として、リチウム過剰型であるＬｉ2 ＭｎＯ3 系の活物質が検討されている。この材料は、充電履歴及び放電履歴に依存して、同一のＳＯＣ（State Of
Charge）に対する電圧値や電気化学的特性が変化する、ヒステリシスという性質を有す
る。
Charge）に対する電圧値や電気化学的特性が変化する、ヒステリシスという性質を有す
る。
二次電池におけるＳＯＣを推定する方法として、二次電池のＯＣＶ（Open Circuit Voltage）とＳＯＣとが一対一対応する相関関係（ＳＯＣ−ＯＣＶ曲線）に基づいてＳＯＣを決定するＯＣＶ法（電圧参照）と、二次電池の充放電電流値を積算してＳＯＣを決定する電流積算法とがある。
ヒステリシスを有する電極材料を用いた場合、ＳＯＣに対して電圧値が一義的に決まらない為、ＯＣＶ法によるＳＯＣの推定が困難である。ＳＯＣ−ＯＣＶ曲線が一義的に決まらない為、ある時点での放電可能エネルギーを予測することも困難である。
電流積算法によってＳＯＣを算出する場合、下記の式（１）を用いる。
ＳＯＣi ＝ＳＯＣi-1 ＋Ｉi ×Δｔi / Ｑ×１００…（１）
ＳＯＣi ：今回のＳＯＣ
ＳＯＣi-1 ：前回のＳＯＣ
Ｉ：電流値
Δｔ：時間間隔
Ｑ：電池容量（available capacity）
ＳＯＣi ＝ＳＯＣi-1 ＋Ｉi ×Δｔi / Ｑ×１００…（１）
ＳＯＣi ：今回のＳＯＣ
ＳＯＣi-1 ：前回のＳＯＣ
Ｉ：電流値
Δｔ：時間間隔
Ｑ：電池容量（available capacity）
電流積算が長期継続されると、電流センサの計測誤差が蓄積する。また、電池容量は経時的に小さくなる為、電流積算法によって推定されるＳＯＣは、その推定誤差が経時的に大きくなる。従来、電流積算を長期継続した場合にＯＣＶ法によりＳＯＣを推定して、誤差の蓄積をリセットするＯＣＶリセットが行われている。
ヒステリシスを有する電極材料を用いた蓄電素子においても、電流積算を継続すると誤差が蓄積する。しかし、ＳＯＣに対して電圧値が一義的に決まらない為、ＯＣＶ法によるＳＯＣの推定を行うこと（ＯＣＶリセットを行うこと）は困難である。
従って、現行のＳＯＣ推定技術ではこのような蓄電素子におけるＳＯＣを精度良く推定することが困難である。
従って、現行のＳＯＣ推定技術ではこのような蓄電素子におけるＳＯＣを精度良く推定することが困難である。
特許文献１に開示の二次電池制御装置においては、充電から放電に切替えた際のＳＯＣである切替時ＳＯＣごとに、放電過程におけるＳＯＣとＯＣＶとの関係を放電時ＯＣＶ情報として記憶する。実際に充電から放電に切替えた際の切替時ＳＯＣと、放電時ＯＣＶ情報とに基づいて、二次電池の放電過程におけるＳＯＣを算出するように、二次電池制御装置は構成されている。
特許文献１の二次電池制御装置においては、充電によって到達した電圧値から、放電時のＳＯＣ−ＯＣＶ曲線を選択し、該ＳＯＣ−ＯＣＶ曲線及び現時点の電圧値に基づいてＳＯＣを推定する。この二次電池制御装置では、充電過程の電圧値に基づいてＳＯＣを推定することができない。複雑なパターンで、充放電が繰り返された場合に、二次電池を高精度に監視することができない。
本発明は、蓄電量−電圧値特性がヒステリシスを示す活物質を含む蓄電素子の蓄電量を推定できる蓄電量推定装置、該蓄電量推定装置を備える蓄電モジュール、蓄電量推定方法、及びコンピュータプログラムを提供することを目的とする。
ここで、蓄電量とは、ＳＯＣ、電力放出可能量等を意味する。
ここで、蓄電量とは、ＳＯＣ、電力放出可能量等を意味する。
本発明に係る蓄電量推定装置は、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定する蓄電量推定装置であって、反応の比率が、他の領域の反応の比率に対して、一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性に基づいて得られる電圧参照蓄電量−電圧値特性を用いて、蓄電量を推定する推定部を備える。
ここで、「一の電気化学反応が生じる場合」とは、「同時に一群として電気化学反応が生じる場合」を含む。「他の電気化学反応が生じる場合」とは、「同時に一群として電気化学反応が生じる場合」を含む。
完全放電状態から満充電状態に充電する、又は満充電状態から完全放電状態に放電する(以下、フル充電又はフル放電という)場合における、蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得してもよい。同一電圧値における両特性の蓄電量の差（Δ蓄電量）は、他の電気化学反応に対応する。一の電気化学反応と他の電気化学反応とは、実質的に独立して生じる。すなわち、他の電気化学反応の反応量はΔ蓄電量に対応し、一の電気化学反応にほとんど影響を与えない。
一の電気化学反応が生じる蓄電量−電圧値放電特性を基準にし、蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させることにより、電圧参照蓄電量−電圧値特性を取得してもよい。即ち、各電圧値を起点とする電圧参照蓄電量−電圧値特性は、一の電気化学反応に基づく、同一の蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させ、各起点の電圧値を上端として切り取った形状を有する。
ここで、「一の電気化学反応が生じる場合」とは、「同時に一群として電気化学反応が生じる場合」を含む。「他の電気化学反応が生じる場合」とは、「同時に一群として電気化学反応が生じる場合」を含む。
完全放電状態から満充電状態に充電する、又は満充電状態から完全放電状態に放電する(以下、フル充電又はフル放電という)場合における、蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得してもよい。同一電圧値における両特性の蓄電量の差（Δ蓄電量）は、他の電気化学反応に対応する。一の電気化学反応と他の電気化学反応とは、実質的に独立して生じる。すなわち、他の電気化学反応の反応量はΔ蓄電量に対応し、一の電気化学反応にほとんど影響を与えない。
一の電気化学反応が生じる蓄電量−電圧値放電特性を基準にし、蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させることにより、電圧参照蓄電量−電圧値特性を取得してもよい。即ち、各電圧値を起点とする電圧参照蓄電量−電圧値特性は、一の電気化学反応に基づく、同一の蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させ、各起点の電圧値を上端として切り取った形状を有する。
上記構成により、高容量で、蓄電量−電圧値特性がヒステリシスを示す活物質を有する蓄電素子の蓄電量が良好に容易に推定できる。
実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、作業量は少ない。
隣り合う電圧参照蓄電量−電圧値特性間に位置すべき電圧参照蓄電量−電圧値特性を内挿計算により求めて補填する場合と比較して、各到達電圧値に対応した電圧参照蓄電量−電圧値特性を、実測した蓄電量−電圧値放電特性から直接的に生成するので、蓄電量推定の精度が向上する。
蓄電素子の劣化に応じて、電圧参照蓄電量−電圧値特性を取得する場合、実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、蓄電素子の使用期間における作業量は少ない。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量を推定できる。充放電曲線に基づいて、ＳＯＣ０％までの放電可能なエネルギー、及びＳＯＣ１００％までに必要な充電エネルギーを予測することができる。
実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、作業量は少ない。
隣り合う電圧参照蓄電量−電圧値特性間に位置すべき電圧参照蓄電量−電圧値特性を内挿計算により求めて補填する場合と比較して、各到達電圧値に対応した電圧参照蓄電量−電圧値特性を、実測した蓄電量−電圧値放電特性から直接的に生成するので、蓄電量推定の精度が向上する。
蓄電素子の劣化に応じて、電圧参照蓄電量−電圧値特性を取得する場合、実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、蓄電素子の使用期間における作業量は少ない。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量を推定できる。充放電曲線に基づいて、ＳＯＣ０％までの放電可能なエネルギー、及びＳＯＣ１００％までに必要な充電エネルギーを予測することができる。
以下、本発明をその実施の形態を示す図面に基づいて具体的に説明する。
（本実施形態の概要）
本実施形態に係る蓄電素子の電極体は、蓄電量−電圧値特性がヒステリシスを有する活物質を含む。
以下、蓄電素子の活物質がＮｉを含むＬｉ過剰型のＬｉＭｅＯ２-Ｌｉ２ＭｎＯ３固溶
体であり、蓄電量がＳＯＣである場合を例として説明する。
図１は、このＬｉ過剰型活物質につき、対極Ｌｉのリチウムセルを用い、電気量と充放電電圧値との関係を求めた結果を示すグラフである。横軸は電気量（ｍＡｈ／ｇ）、縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。ここで、電気量はＳＯＣに対応する。
図１に示すように、ＳＯＣの増加（充電）と、減少（放電）とで電圧値が異なる。即ち同一ＳＯＣに対する電圧値が異なり、ヒステリシスを有する。この活物質の場合、同一ＳＯＣに対する電位差は、高ＳＯＣ領域では低ＳＯＣ領域より小さく、ヒステリシスが小さい。
（本実施形態の概要）
本実施形態に係る蓄電素子の電極体は、蓄電量−電圧値特性がヒステリシスを有する活物質を含む。
以下、蓄電素子の活物質がＮｉを含むＬｉ過剰型のＬｉＭｅＯ２-Ｌｉ２ＭｎＯ３固溶
体であり、蓄電量がＳＯＣである場合を例として説明する。
図１は、このＬｉ過剰型活物質につき、対極Ｌｉのリチウムセルを用い、電気量と充放電電圧値との関係を求めた結果を示すグラフである。横軸は電気量（ｍＡｈ／ｇ）、縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。ここで、電気量はＳＯＣに対応する。
図１に示すように、ＳＯＣの増加（充電）と、減少（放電）とで電圧値が異なる。即ち同一ＳＯＣに対する電圧値が異なり、ヒステリシスを有する。この活物質の場合、同一ＳＯＣに対する電位差は、高ＳＯＣ領域では低ＳＯＣ領域より小さく、ヒステリシスが小さい。
本実施形態においては、ヒステリシスが小さく、後述する電圧参照ＳＯＣ−ＯＣＶ曲線を用いて電圧値からＳＯＣを推定することができる領域を判断し、該領域においてＳＯＣを推定する。
図２は、電気量に対する、Ｘ線吸収分光測定（ＸＡＦＳ測定）によって算出したＬｉ過剰型活物質のＮｉのＫ吸収端エネルギーの推移を示すグラフである。横軸は電気量（ｍＡｈ／ｇ）、縦軸はＮｉのＫ吸収端エネルギーＥ0 （ｅＶ）である。
図３は、充放電時におけるＮｉのＫ吸収端エネルギーの推移を示すグラフである。横軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ）、縦軸はＮｉのＫ吸収端エネルギーＥ0 （ｅＶ）である。
図２は、電気量に対する、Ｘ線吸収分光測定（ＸＡＦＳ測定）によって算出したＬｉ過剰型活物質のＮｉのＫ吸収端エネルギーの推移を示すグラフである。横軸は電気量（ｍＡｈ／ｇ）、縦軸はＮｉのＫ吸収端エネルギーＥ0 （ｅＶ）である。
図３は、充放電時におけるＮｉのＫ吸収端エネルギーの推移を示すグラフである。横軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ）、縦軸はＮｉのＫ吸収端エネルギーＥ0 （ｅＶ）である。
図２に示すように、高ＳＯＣ領域において、充電反応のＮｉのＫ吸収端エネルギー推移が放電反応のエネルギー推移と一致しない。低ＳＯＣ領域において、放電反応のエネルギー推移が充電反応のエネルギー推移と一致しない。即ちヒステリシスを有する、Ｎｉ以外のレドックス反応が主として生じていることが分かる（これをＡの反応とする）。Ａの反応は、高ＳＯＣ領域においては酸化反応であり、低ＳＯＣ領域においては還元反応である。
中間のＳＯＣ領域では、充電反応及び放電反応のＮｉのＫ吸収端エネルギーがＳＯＣに対して略直線的に変化している。
中間のＳＯＣ領域では、充電反応及び放電反応のＮｉのＫ吸収端エネルギーがＳＯＣに対して略直線的に変化している。
図３に示すように、充放電電圧値が略３．７〜４．５Ｖである高ＳＯＣ領域において、ＮｉのＫ吸収端エネルギーは、充電と放電とで略一致している。ＮｉのＫ吸収端エネルギーが一致している場合、即ちＮｉの価数等が等しく、この電圧範囲において、Ｎｉの価数変化と電圧値とが概ね１：１に対応しており、Ｎｉは可逆的に反応していると考えられる。即ち、該ＳＯＣ領域において、ＳＯＣ−ＯＣＰ特性が示すヒステリシスが小さいレドックス反応が主として生じている（これをＢの反応とする）。なお、ＯＣＰは開放電位を意味する。
該ＳＯＣ領域においてＢの反応量はＡの反応量より多く、結果として、低ＳＯＣ領域よりヒステリシスが小さい。
該ＳＯＣ領域においてＢの反応量はＡの反応量より多く、結果として、低ＳＯＣ領域よりヒステリシスが小さい。
本実施形態においては、主としてＢの反応が生じる領域の低い方の電圧値（下限電圧値）を実験により求める。下限電圧値において、ヒステリシスの有無が実質的に切り替わる。Ｂの反応の酸化量及び還元量は小さいと考えられる。そして、電圧値の昇降に基づいて、充電状態又は放電状態が下限電圧値以上の電圧領域に相当する領域にあると判定した場合、到達電圧値に基づいて電圧参照によりＳＯＣを推定する。
なお、ここでは、一例として、Ｎｉの酸化還元反応だけに注目して説明しているが、Ｂの反応はＮｉの酸化還元反応に限定されるものではない。Ｂの反応は、充放電の推移に応じて活物質により生じる一又は一群の反応のうち、蓄電量−電圧値特性のヒステリシスが小さい反応をいう。
なお、ここでは、一例として、Ｎｉの酸化還元反応だけに注目して説明しているが、Ｂの反応はＮｉの酸化還元反応に限定されるものではない。Ｂの反応は、充放電の推移に応じて活物質により生じる一又は一群の反応のうち、蓄電量−電圧値特性のヒステリシスが小さい反応をいう。
次に、蓄電素子の負極がヒステリシスが大きい活物質を含む場合につき説明する。一例として、負極が活物質としてＳｉＯ及びグラファイトを含む場合を挙げる。ＳｉＯの電気化学反応が生じるときのヒステリシスは、グラファイトの電気化学反応が生じるときのヒステリシスより大きい。
図４は、この蓄電素子につき、対極Ｌｉのリチウムセルを用い、電気量と充放電電圧値との関係を求めた結果を示すグラフである。横軸は電気量（ｍＡｈ／ｇ）、縦軸は充放電電圧値（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。ここで、電気量はＳＯＣに対応する。
図４に示すように、充電曲線と、放電曲線とで電圧値が異なる。即ち同一ＳＯＣに対する電圧値が異なり、ヒステリシスを有する。この活物質の場合、同一ＳＯＣに対する電位差は、高ＳＯＣ領域では低ＳＯＣ領域より小さく、ヒステリシスが小さい。
図４は、この蓄電素子につき、対極Ｌｉのリチウムセルを用い、電気量と充放電電圧値との関係を求めた結果を示すグラフである。横軸は電気量（ｍＡｈ／ｇ）、縦軸は充放電電圧値（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。ここで、電気量はＳＯＣに対応する。
図４に示すように、充電曲線と、放電曲線とで電圧値が異なる。即ち同一ＳＯＣに対する電圧値が異なり、ヒステリシスを有する。この活物質の場合、同一ＳＯＣに対する電位差は、高ＳＯＣ領域では低ＳＯＣ領域より小さく、ヒステリシスが小さい。
図５は、ＳＯＣ（又は電圧値）が高くなるのに従い、ヒステリシスが大きい領域と小さい領域とが交互に表れる場合の充放電曲線である。横軸はＳＯＣ（％）、縦軸は電圧値（Ｖ）である。
正極が組成の異なるＬｉ過剰型の活物質を複数含む場合、負極がヒステリシスが大きい活物質を複数含む場合、並びに正極及び負極夫々がヒステリシスが大きい活物質を含む場合、ヒステリシスが大きい領域と小さい領域とが交互に表れる、又は重なって表れることがある。
正極が組成の異なるＬｉ過剰型の活物質を複数含む場合、負極がヒステリシスが大きい活物質を複数含む場合、並びに正極及び負極夫々がヒステリシスが大きい活物質を含む場合、ヒステリシスが大きい領域と小さい領域とが交互に表れる、又は重なって表れることがある。
図５の電圧値がａ〜ｂである領域（２）においては、電圧値がａ以下である領域（１）よりヒステリシスが小さい。領域（２）においては、ヒステリシスが大きいＣの反応とヒステリシスが小さいＤの反応とが生じている。領域（２）においてＤの反応量が多いので、結果として領域（１）よりヒステリシスが小さい。
図５の電圧値がｃ以上である領域（４）においては、電圧値がｂ〜ｃである領域（３）よりヒステリシスが小さい。領域（４）においては、ヒステリシスが大きいＥの反応とヒステリシスが小さいＦの反応とが生じている。領域（４）においてＥの反応量が多いので、結果として領域（３）よりヒステリシスが小さい。
図５の電圧値がｃ以上である領域（４）においては、電圧値がｂ〜ｃである領域（３）よりヒステリシスが小さい。領域（４）においては、ヒステリシスが大きいＥの反応とヒステリシスが小さいＦの反応とが生じている。領域（４）においてＥの反応量が多いので、結果として領域（３）よりヒステリシスが小さい。
領域（２）の下限電圧値ａ、及び領域（４）の下限電圧値ｃを実験により求める。電圧値の昇降に基づいて、充電状態又は放電状態が、下限電圧値ａ以上の電圧領域に相当する領域（２）にあると判定した場合、及び下限電圧値ｃ以上の電圧領域に相当する領域（４）にあると判定した場合、夫々、後述する到達電圧値に基づき、電圧参照によりＳＯＣを推定する。
なお、電圧値の領域は上述したように２又は４の領域に分かれる場合には限定されない。正極又は負極の活物質に応じて、複数の電気化学反応が生じ、ヒステリシスが大きい領域と小さい領域とが交互に表れる場合に、ＳＯＣが相対的に高く、ヒステリシスが小さい領域において、電圧参照によりＳＯＣを推定する。
以下、前記Ａの反応及びＢの反応を有する活物質につき、充放電実験により、前記反応ＢがＳＯＣ４０％以上で多く生じ、これに対応する下限電圧値がＥ0 Ｖであることが求められたとして説明する。
下限電圧値としてＯＣＶを測定できる場合、該下限電圧値は一定であってもよい。下限電圧値としてＣＣＶを測定する場合、蓄電素子の使用に伴う劣化の程度に対応して、下限電圧値を下げる等、更新してもよい。蓄電素子の劣化の原因として、内部抵抗の上昇、容量バランスのずれの増大等が挙げられる。容量バランスのずれとは、例えば、正極における充放電反応以外の副反応の量と、負極における充放電反応以外の副反応の量とに差が生じることによって、正極及び負極のうちの一方が完全には充電されないようになり、正極と負極の、可逆的に電荷イオンが電極から出入りできる容量が相違するようになることをいう。一般的なリチウムイオン電池では、正極における副反応量が、負極における副反応量よりも小さいために、「容量バランスのずれ」が増大すると、負極を完全に充電することができなくなり、蓄電素子から可逆的に取り出せる電気量が減少する。
充電して電圧値が下限電圧値を超えた後、即ち下限電圧値よりも貴になった後、到達した最大の電圧値を到達電圧値とする。
蓄電量推定装置のメモリのテーブルには、下限電圧値から複数の到達電圧値までの複数の電圧参照ＳＯＣ−ＯＣＶ曲線（以下、参照ＳＯＣ−ＯＣＶ曲線という）が格納されている。例えば下限電圧値Ｅ0 Ｖから到達電圧値Ｅ1 ＶまでのＳＯＣ−ＯＣＶ曲線ａ、下限電圧値Ｅ0 Ｖから到達電圧値Ｅ2 ＶまでのＳＯＣ−ＯＣＶ曲線ｂ、下限電圧値Ｅ0 Ｖから到達電圧値Ｅ3 ＶまでのＳＯＣ−ＯＣＶ曲線ｃが格納されている。ここで、Ｅ1 ＞Ｅ2 ＞Ｅ3 である。なお、ＳＯＣ−ＯＣＶ曲線ａ，ｂ，ｃは、後述する比較試験においても参照されるが、図示はしていない。テーブルには、離散的ではなく、全ての到達電圧値に対応したＳＯＣ−ＯＣＶ曲線が連続的に格納されている。連続的には格納せず、ＳＯＣ−ＯＣＶ曲線を離散的に用意し、隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填してもよい。
蓄電量推定装置のメモリのテーブルには、下限電圧値から複数の到達電圧値までの複数の電圧参照ＳＯＣ−ＯＣＶ曲線（以下、参照ＳＯＣ−ＯＣＶ曲線という）が格納されている。例えば下限電圧値Ｅ0 Ｖから到達電圧値Ｅ1 ＶまでのＳＯＣ−ＯＣＶ曲線ａ、下限電圧値Ｅ0 Ｖから到達電圧値Ｅ2 ＶまでのＳＯＣ−ＯＣＶ曲線ｂ、下限電圧値Ｅ0 Ｖから到達電圧値Ｅ3 ＶまでのＳＯＣ−ＯＣＶ曲線ｃが格納されている。ここで、Ｅ1 ＞Ｅ2 ＞Ｅ3 である。なお、ＳＯＣ−ＯＣＶ曲線ａ，ｂ，ｃは、後述する比較試験においても参照されるが、図示はしていない。テーブルには、離散的ではなく、全ての到達電圧値に対応したＳＯＣ−ＯＣＶ曲線が連続的に格納されている。連続的には格納せず、ＳＯＣ−ＯＣＶ曲線を離散的に用意し、隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填してもよい。
参照ＳＯＣ−ＯＣＶ曲線は、以下のようにして求められる。
ＳＯＣが４０％から１００％までの各点のＳＯＣ（％）につき、各ＳＯＣ（％）→４０％→１００％と変化させた場合の放電ＯＣＶ曲線及び充電ＯＣＶ曲線を求める。放電ＯＣＶ曲線は、たとえば、放電方向に微小な電流を通電させて、その際の電圧値を測定することで取得できる。又は、充電状態から各ＳＯＣまで放電し、休止させて電圧値が安定した電圧値を測定することもできる。同様に、充電ＯＣＶは、充電方向で上記の測定を実施すれば取得できる。ＳＯＣ４０％以上においても、前記活物質が僅かにヒステリシスを有している為、放電ＯＣＶ曲線及び充電ＯＣＶ曲線を平均化したＯＣＶ曲線を使用するのが好ましい。放電ＯＣＶ曲線及び充電ＯＣＶ曲線、又はそれらを補正したものを用いてもよい。
なお、始めに放電ＯＣＰ曲線及び充電ＯＣＰ曲線を求めた後、蓄電素子の参照ＳＯＣ−ＯＣＶ曲線に補正してもよい。
ＳＯＣが４０％から１００％までの各点のＳＯＣ（％）につき、各ＳＯＣ（％）→４０％→１００％と変化させた場合の放電ＯＣＶ曲線及び充電ＯＣＶ曲線を求める。放電ＯＣＶ曲線は、たとえば、放電方向に微小な電流を通電させて、その際の電圧値を測定することで取得できる。又は、充電状態から各ＳＯＣまで放電し、休止させて電圧値が安定した電圧値を測定することもできる。同様に、充電ＯＣＶは、充電方向で上記の測定を実施すれば取得できる。ＳＯＣ４０％以上においても、前記活物質が僅かにヒステリシスを有している為、放電ＯＣＶ曲線及び充電ＯＣＶ曲線を平均化したＯＣＶ曲線を使用するのが好ましい。放電ＯＣＶ曲線及び充電ＯＣＶ曲線、又はそれらを補正したものを用いてもよい。
なお、始めに放電ＯＣＰ曲線及び充電ＯＣＰ曲線を求めた後、蓄電素子の参照ＳＯＣ−ＯＣＶ曲線に補正してもよい。
充放電を繰り返した場合に、上記参照ＳＯＣ−ＯＣＶ曲線によるＳＯＣ推定と、従来の電流積算によるＳＯＣ推定との差異を比較した結果について説明する。
図６は、充放電時における、時間に対する電圧値の推移を示すグラフである。横軸は時間（秒）、縦軸は充放電電圧値（ＶｖｓＬｉ／Ｌｉ+ ）である。図６においては微小電流により充放電を実施している為、通電中の電圧値は、ＯＣＶとほぼ同じ値を示すことを確認している。
図６に示すように、１回目の充電が行われて電圧値が下限電圧値Ｅ0 Ｖを超え、Ｅ3 Ｖに到達した後、１回目の放電が行われた。電圧値がＥ0 Ｖに達した後、２回目の充電が行われ、電圧値はＥ1 Ｖに達し、その後、２回目の放電が行われた。
図６は、充放電時における、時間に対する電圧値の推移を示すグラフである。横軸は時間（秒）、縦軸は充放電電圧値（ＶｖｓＬｉ／Ｌｉ+ ）である。図６においては微小電流により充放電を実施している為、通電中の電圧値は、ＯＣＶとほぼ同じ値を示すことを確認している。
図６に示すように、１回目の充電が行われて電圧値が下限電圧値Ｅ0 Ｖを超え、Ｅ3 Ｖに到達した後、１回目の放電が行われた。電圧値がＥ0 Ｖに達した後、２回目の充電が行われ、電圧値はＥ1 Ｖに達し、その後、２回目の放電が行われた。
テーブルには、１回目の到達電圧値としてＥ3 Ｖが記憶される。２回目の充電でＥ3 Ｖを超えた時点で到達電圧値が更新される。１回目の放電及び２回目の充電においてＥ3 Ｖに達するまでの間は、前記ＳＯＣ−ＯＣＶ曲線ｃが使用される。２回目の充電のＥ3 ＶからＥ1 Ｖまでの間はテーブルに記憶された別のＳＯＣ−ＯＣＶ曲線が使用される。２回目の放電のＥ1 Ｖから下限電圧値のＥ0 Ｖまでの間は前記ＳＯＣ−ＯＣＶ曲線ａが使用される。
図７は、初期品の蓄電素子につき、図６に示す充放電を行った場合の、本実施形態の電圧参照によりＳＯＣを推定したときと、従来の電流積算によりＳＯＣを推定したときとの差異を示すグラフである。横軸は時間（秒）、左側の縦軸はＳＯＣ（％）、右側の縦軸は前記差異（％）である。対照としての電流積算によるＳＯＣの推定は、事前に放電容量を確認し、かつ精度の高い電流計を使用している為、前記式（１）中のＱの放電容量及びＩの電流値が正確である。真値に近似していると考えられる。
図中、ｄは電流積算により求めたＳＯＣの推移、ｅは、前記ＳＯＣ−ＯＣＶ曲線ｃ及びａを用いて電圧参照により求めたＳＯＣの推移、ｆは差異である。差異は、（電圧参照により算出したＳＯＣ）−（電流積算により算出したＳＯＣ）により求めた。
図７より差異は略±４％未満であり、小さいことが分かる。
以上より、参照ＳＯＣ−ＯＣＶ曲線によるＳＯＣ推定の精度が良好であることが確認された。
図中、ｄは電流積算により求めたＳＯＣの推移、ｅは、前記ＳＯＣ−ＯＣＶ曲線ｃ及びａを用いて電圧参照により求めたＳＯＣの推移、ｆは差異である。差異は、（電圧参照により算出したＳＯＣ）−（電流積算により算出したＳＯＣ）により求めた。
図７より差異は略±４％未満であり、小さいことが分かる。
以上より、参照ＳＯＣ−ＯＣＶ曲線によるＳＯＣ推定の精度が良好であることが確認された。
しかし、全ての到達電圧値に対応した参照ＳＯＣ−ＯＣＶ曲線を連続的にテーブルに格納する場合、上述の放電ＯＣＶ曲線及び充電ＯＣＶ曲線を実測するので、作業量が膨大である。
参照ＳＯＣ−ＯＣＶ曲線を離散的に用意し、隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填した場合、ＳＯＣ推定の精度が低くなる。
しかも、蓄電素子の劣化に応じて所定時間間隔で電圧参照ＳＯＣ−ＯＣＶ曲線を用意することになるので、蓄電素子の使用期間におけるデータ量は非常に膨大になる。
参照ＳＯＣ−ＯＣＶ曲線を離散的に用意し、隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填した場合、ＳＯＣ推定の精度が低くなる。
しかも、蓄電素子の劣化に応じて所定時間間隔で電圧参照ＳＯＣ−ＯＣＶ曲線を用意することになるので、蓄電素子の使用期間におけるデータ量は非常に膨大になる。
図２及び図３に示すように、高いＳＯＣの領域においては、ヒステリシスの小さいＢの酸化還元反応と、Ａの酸化反応のみが生じ、Ａの還元反応はほとんど生じない。また、ある電位におけるＡの酸化反応が一度生じると、その酸化反応と対となるＡの還元反応が生じるまでは、その電圧においてＡの酸化反応は生じない。従って、一旦、Ａの酸化反応が生じた後において、下限電圧値から到達上限電圧までは、ヒステリシスのないＢの酸化還元反応のみが生じ、その充放電曲線形状は、同じ電圧におけるフル放電時の放電曲線形状と一致する。
フル放電したときの放電ＳＯＣ−ＯＣＶ曲線、及びフル充電したときの充電ＳＯＣ−ＯＣＶ曲線を用意する。放電ＳＯＣ−ＯＣＶ曲線上の或る電圧値におけるＳＯＣと、充電ＳＯＣ−ＯＣＶ曲線上の同一電圧値におけるＳＯＣとの差分をΔＳＯＣとする。
Ａの反応とＢの反応とは実質的に独立して生じる。Ａの反応量はΔＳＯＣに対応し、Ｂの反応に影響を与えない。そこで、Ｂの反応の放電ＳＯＣ−ＯＣＶ曲線を基準にしてΔＳＯＣで補正することにより、容易に参照ＳＯＣ−ＯＣＶ曲線を取得できると考えた。各到達電圧値を起点とする参照ＳＯＣ−ＯＣＶ曲線は、同一のＢの反応に基づく、同一の放電ＳＯＣ−ＯＣＶ曲線を充電ＳＯＣ−ＯＣＶ曲線に向けて移動させて補正し、各起点の到達電圧値を上端として切り取った形状を有する。
フル放電したときの放電ＳＯＣ−ＯＣＶ曲線、及びフル充電したときの充電ＳＯＣ−ＯＣＶ曲線を用意する。放電ＳＯＣ−ＯＣＶ曲線上の或る電圧値におけるＳＯＣと、充電ＳＯＣ−ＯＣＶ曲線上の同一電圧値におけるＳＯＣとの差分をΔＳＯＣとする。
Ａの反応とＢの反応とは実質的に独立して生じる。Ａの反応量はΔＳＯＣに対応し、Ｂの反応に影響を与えない。そこで、Ｂの反応の放電ＳＯＣ−ＯＣＶ曲線を基準にしてΔＳＯＣで補正することにより、容易に参照ＳＯＣ−ＯＣＶ曲線を取得できると考えた。各到達電圧値を起点とする参照ＳＯＣ−ＯＣＶ曲線は、同一のＢの反応に基づく、同一の放電ＳＯＣ−ＯＣＶ曲線を充電ＳＯＣ−ＯＣＶ曲線に向けて移動させて補正し、各起点の到達電圧値を上端として切り取った形状を有する。
図８は、この原理を説明するための説明図である。横軸はＳＯＣ［％］、縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ）である。
図８中、（１）はフル放電ＳＯＣ−ＯＣＶ曲線、（２）はフル充電ＳＯＣ−ＯＣＶ曲線である。放電ＳＯＣ−ＯＣＶ曲線（１）上の或る到達電圧値におけるＳＯＣと、充電ＳＯＣ−ＯＣＶ曲線（２）上の同一到達電圧値におけるＳＯＣとの差分がΔＳＯＣである。このΔＳＯＣ分、放電ＳＯＣ−ＯＣＶ曲線（１）を横軸に平行に、充電ＳＯＣ−ＯＣＶ曲線（２）側にスライドさせることで、ＳＯＣ−ＯＣＶ曲線（３）が得られる。下限電圧値から到達電圧値までの範囲におけるＳＯＣ−ＯＣＶ曲線（３）を参照ＳＯＣ−ＯＣＶ曲線とする。ＳＯＣを推定する場合、該参照ＳＯＣ−ＯＣＶ曲線を用い、取得した電圧値に対応するＳＯＣを読み取る。
図８中、（１）はフル放電ＳＯＣ−ＯＣＶ曲線、（２）はフル充電ＳＯＣ−ＯＣＶ曲線である。放電ＳＯＣ−ＯＣＶ曲線（１）上の或る到達電圧値におけるＳＯＣと、充電ＳＯＣ−ＯＣＶ曲線（２）上の同一到達電圧値におけるＳＯＣとの差分がΔＳＯＣである。このΔＳＯＣ分、放電ＳＯＣ−ＯＣＶ曲線（１）を横軸に平行に、充電ＳＯＣ−ＯＣＶ曲線（２）側にスライドさせることで、ＳＯＣ−ＯＣＶ曲線（３）が得られる。下限電圧値から到達電圧値までの範囲におけるＳＯＣ−ＯＣＶ曲線（３）を参照ＳＯＣ−ＯＣＶ曲線とする。ＳＯＣを推定する場合、該参照ＳＯＣ−ＯＣＶ曲線を用い、取得した電圧値に対応するＳＯＣを読み取る。
図９は、ＳＯＣが４０〜６０％である場合に、実測により得られたＳＯＣ−ＯＣＶ曲線である。横軸はＳＯＣ（％）、縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。
図１０は、ＳＯＣが４０〜６０％である場合に、上述の方法により得られた参照ＳＯＣ−ＯＣＶ曲線である。横軸はＳＯＣ（％）、縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。
図１１は、図９のＳＯＣ−ＯＣＶ曲線ａと図１０の参照ＳＯＣ−ＯＣＶ曲線ｂとの差異を示すグラフである。横軸はＳＯＣ（％）、左側の縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）、右側の縦軸は差異（Ｖ）である。差異は下記の式により求められる。
差異＝｜（参照ＳＯＣ−ＯＣＶ曲線の電圧値）−（実測ＳＯＣ−ＯＣＶ曲線の電圧値）｜
上述の方法により求めた参照ＳＯＣ−ＯＣＶ曲線は、実測により求めた参照ＳＯＣ−ＯＣＶ曲線と一致することが確認された。即ち、下限電圧値から到達上限電圧までの参照ＳＯＣ−ＯＣＶ曲線の形状は、同じ電圧におけるフル放電時の放電曲線の形状と一致する。
上記方法により求めた参照ＳＯＣ−ＯＣＶ曲線を用いることにより、図７に示す結果で証明されたように、ＳＯＣを高精度に推定できる。
図１０は、ＳＯＣが４０〜６０％である場合に、上述の方法により得られた参照ＳＯＣ−ＯＣＶ曲線である。横軸はＳＯＣ（％）、縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）である。
図１１は、図９のＳＯＣ−ＯＣＶ曲線ａと図１０の参照ＳＯＣ−ＯＣＶ曲線ｂとの差異を示すグラフである。横軸はＳＯＣ（％）、左側の縦軸は充放電電圧値Ｅ（ＶｖｓＬｉ／Ｌｉ+ ：Ｌｉ／Ｌｉ+平衡電位を基準にしたときの電位差）、右側の縦軸は差異（Ｖ）である。差異は下記の式により求められる。
差異＝｜（参照ＳＯＣ−ＯＣＶ曲線の電圧値）−（実測ＳＯＣ−ＯＣＶ曲線の電圧値）｜
上述の方法により求めた参照ＳＯＣ−ＯＣＶ曲線は、実測により求めた参照ＳＯＣ−ＯＣＶ曲線と一致することが確認された。即ち、下限電圧値から到達上限電圧までの参照ＳＯＣ−ＯＣＶ曲線の形状は、同じ電圧におけるフル放電時の放電曲線の形状と一致する。
上記方法により求めた参照ＳＯＣ−ＯＣＶ曲線を用いることにより、図７に示す結果で証明されたように、ＳＯＣを高精度に推定できる。
以上の方法により、下限電圧値から各到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を容易に取得することができる。実測するのはフル放電ＳＯＣ−ＯＣＶ曲線及びフル充電ＳＯＣ−ＯＣＶ曲線のみであり、事前に検討すべき作業量が著しく軽減される。
隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填する場合と比較して、各到達電圧値に対応した参照ＳＯＣ−ＯＣＶ曲線を、実測したフル充放電ＳＯＣ−ＯＣＶ曲線から直接的に生成するので、ＳＯＣ推定の精度が向上する。
ＳＯＣ−ＯＣＶ曲線が使用に伴い、変化する場合、蓄電素子の劣化に応じて、参照ＳＯＣ−ＯＣＶ曲線を生成する。この場合、実測するのはフル充放電ＳＯＣ−ＯＣＶ曲線のみであり、蓄電素子の使用期間における作業量が著しく軽減される。ＳＯＣの推定は、現在の電気化学挙動及び使用履歴から推定するのが好ましい。
隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填する場合と比較して、各到達電圧値に対応した参照ＳＯＣ−ＯＣＶ曲線を、実測したフル充放電ＳＯＣ−ＯＣＶ曲線から直接的に生成するので、ＳＯＣ推定の精度が向上する。
ＳＯＣ−ＯＣＶ曲線が使用に伴い、変化する場合、蓄電素子の劣化に応じて、参照ＳＯＣ−ＯＣＶ曲線を生成する。この場合、実測するのはフル充放電ＳＯＣ−ＯＣＶ曲線のみであり、蓄電素子の使用期間における作業量が著しく軽減される。ＳＯＣの推定は、現在の電気化学挙動及び使用履歴から推定するのが好ましい。
（実施形態１）
以下、実施形態１として、車両に搭載される蓄電モジュールを例に挙げて説明する。
図１２は、蓄電モジュールの一例を示す。蓄電モジュール５０は、複数の蓄電素子２００と、監視装置１００と、それらを収容する収容ケース３００とを備えている。蓄電モジュール５０は、電気自動車（ＥＶ）や、プラグインハイブリッド電気自動車（ＰＨＥＶ）の動力源として使用されてもよい。
蓄電素子２００は、角形セルに限定されず、円筒形セルやパウチセルであってもよい。
監視装置１００は、複数の蓄電素子２００と対向して配置される回路基板であってもよい。監視装置１００は、蓄電素子２００の状態を監視する。監視装置１００が、蓄電量推定装置であってもよい。代替的に、監視装置１００と有線接続または無線接続されるコンピュータやサーバが、監視装置１００が出力する情報に基づいて蓄電量推定方法を実行してもよい。
以下、実施形態１として、車両に搭載される蓄電モジュールを例に挙げて説明する。
図１２は、蓄電モジュールの一例を示す。蓄電モジュール５０は、複数の蓄電素子２００と、監視装置１００と、それらを収容する収容ケース３００とを備えている。蓄電モジュール５０は、電気自動車（ＥＶ）や、プラグインハイブリッド電気自動車（ＰＨＥＶ）の動力源として使用されてもよい。
蓄電素子２００は、角形セルに限定されず、円筒形セルやパウチセルであってもよい。
監視装置１００は、複数の蓄電素子２００と対向して配置される回路基板であってもよい。監視装置１００は、蓄電素子２００の状態を監視する。監視装置１００が、蓄電量推定装置であってもよい。代替的に、監視装置１００と有線接続または無線接続されるコンピュータやサーバが、監視装置１００が出力する情報に基づいて蓄電量推定方法を実行してもよい。
図１３は、蓄電モジュールの他の例を示す。蓄電モジュール（以下、電池モジュールという）１は、エンジン車両に好適に搭載される、１２ボルト電源や、４８ボルト電源であってもよい。図１３は１２Ｖ電源用の電池モジュール１の斜視図、図１４は電池モジュール１の分解斜視図、図１５は電池モジュール１のブロック図である。
電池モジュール１は直方体状のケース２を有する。ケース２に複数のリチウムイオン二次電池（以下、電池という）３、複数のバスバー４、ＢＭＵ（Battery Management Unit
）６、電流センサ７が収容される。
電池モジュール１は直方体状のケース２を有する。ケース２に複数のリチウムイオン二次電池（以下、電池という）３、複数のバスバー４、ＢＭＵ（Battery Management Unit
）６、電流センサ７が収容される。
電池３は、直方体状のケース３１と、ケース３１の一側面に設けられた、極性が異なる一対の端子３２，３２とを備える。ケース３１には、正極板、セパレータ、及び負極板を積層した電極体３３が収容されている。
電極体３３の正極板が有する正極活物質及び負極板が有する負極活物質の少なくとも一方は、充放電の推移に応じて２以上の電気化学反応を生じる。一の電気化学反応が生じるときに示す蓄電量−電圧値特性のヒステリシスが、他の電気化学反応が生じるときの前記ヒステリシスより小さい。
正極活物質としては、上述のＬｉＭｅＯ２-Ｌｉ２ＭｎＯ３固溶体、Ｌｉ2Ｏ−ＬｉＭｅＯ2固溶体、Ｌｉ3ＮｂＯ4 −ＬｉＭｅＯ2固溶体、Ｌｉ4 ＷＯ5 −ＬｉＭｅＯ2固溶体、Ｌｉ4 ＴｅＯ5 −ＬｉＭｅＯ2固溶体、Ｌｉ3ＳｂＯ4 −ＬｉＦｅＯ2固溶体、Ｌｉ2ＲｕＯ3 −ＬｉＭｅＯ2固溶体、Ｌｉ2ＲｕＯ3 −Ｌｉ2 ＭｅＯ3 固溶体等のＬｉ過剰型活物質が挙げられる。負極活物質としては、ハードカーボン、Ｓｉ、Ｓｎ、Ｃｄ、Ｚｎ、Ａｌ、Ｂｉ、Ｐｂ、Ｇｅ、Ａｇ等の金属若しくは合金、又はこれらを含むカルコゲン化物等が挙げられる。カルコゲン化物の一例として、ＳｉＯなどが挙げられる。
本発明の技術は、これらの正極活物質及び負極活物質の少なくとも一方が含まれていれば適用可能である。
正極活物質としては、上述のＬｉＭｅＯ２-Ｌｉ２ＭｎＯ３固溶体、Ｌｉ2Ｏ−ＬｉＭｅＯ2固溶体、Ｌｉ3ＮｂＯ4 −ＬｉＭｅＯ2固溶体、Ｌｉ4 ＷＯ5 −ＬｉＭｅＯ2固溶体、Ｌｉ4 ＴｅＯ5 −ＬｉＭｅＯ2固溶体、Ｌｉ3ＳｂＯ4 −ＬｉＦｅＯ2固溶体、Ｌｉ2ＲｕＯ3 −ＬｉＭｅＯ2固溶体、Ｌｉ2ＲｕＯ3 −Ｌｉ2 ＭｅＯ3 固溶体等のＬｉ過剰型活物質が挙げられる。負極活物質としては、ハードカーボン、Ｓｉ、Ｓｎ、Ｃｄ、Ｚｎ、Ａｌ、Ｂｉ、Ｐｂ、Ｇｅ、Ａｇ等の金属若しくは合金、又はこれらを含むカルコゲン化物等が挙げられる。カルコゲン化物の一例として、ＳｉＯなどが挙げられる。
本発明の技術は、これらの正極活物質及び負極活物質の少なくとも一方が含まれていれば適用可能である。
ケース２は合成樹脂製である。ケース２は、ケース本体２１と、ケース本体２１の開口部を閉塞する蓋部２２と、蓋部２２の外面に設けられたＢＭＵ収容部２３と、ＢＭＵ収容部２３を覆うカバー２４と、中蓋２５と、仕切り板２６とを備える。中蓋２５や仕切り板２６は、設けられなくてもよい。
ケース本体２１の各仕切り板２６の間に、電池３が挿入されている。
ケース本体２１の各仕切り板２６の間に、電池３が挿入されている。
中蓋２５には、複数の金属製のバスバー４が載置されている。電池３の端子３２が設けられている端子面に中蓋２５が配置されて、隣り合う電池３の隣り合う端子３２がバスバー４により接続され、電池３が直列に接続されている。
ＢＭＵ収容部２３は箱状をなし、一長側面の中央部に、外側に角型に突出した突出部２３ａを有する。蓋部２２における突出部２３ａの両側には、鉛合金等の金属製で、極性が異なる一対の外部端子５，５が設けられている。ＢＭＵ６は、基板に情報処理部６０、電圧計測部８、及び電流計測部９を実装してなる。ＢＭＵ収容部２３にＢＭＵ６を収容し、カバー２４によりＢＭＵ収容部２３を覆うことにより、電池３とＢＭＵ６とが接続される。
図１５に示すように、情報処理部６０は、ＣＰＵ６２と、メモリ６３とを備える。
メモリ６３には、本実施形態に係るＳＯＣ推定プログラム６３ａと、フルの充電ＳＯＣ−ＯＣＶ曲線及び放電ＳＯＣ−ＯＣＶ曲線が格納されたテーブル６３ｂとが記憶されている。ＳＯＣ−ＯＣＶ曲線はテーブル６３ｂに格納する場合に限定されず、メモリ６３に関数式として格納してもよい。
フルの充電ＳＯＣ−ＯＣＶ曲線及び放電ＳＯＣ−ＯＣＶ曲線は、所定の期間間隔で実測により取得する。電池３の劣化に従って、充放電ＳＯＣ−ＯＣＶ曲線が更新される。なお、電池３の使用の都度、充放電ＳＯＣ−ＯＣＶ曲線を実測してもよい。
メモリ６３には、本実施形態に係るＳＯＣ推定プログラム６３ａと、フルの充電ＳＯＣ−ＯＣＶ曲線及び放電ＳＯＣ−ＯＣＶ曲線が格納されたテーブル６３ｂとが記憶されている。ＳＯＣ−ＯＣＶ曲線はテーブル６３ｂに格納する場合に限定されず、メモリ６３に関数式として格納してもよい。
フルの充電ＳＯＣ−ＯＣＶ曲線及び放電ＳＯＣ−ＯＣＶ曲線は、所定の期間間隔で実測により取得する。電池３の劣化に従って、充放電ＳＯＣ−ＯＣＶ曲線が更新される。なお、電池３の使用の都度、充放電ＳＯＣ−ＯＣＶ曲線を実測してもよい。
ＳＯＣ推定プログラム６３ａは、例えば、ＣＤ−ＲＯＭやＤＶＤ−ＲＯＭ、ＵＳＢメモリ等のコンピュータ読み取り可能な記録媒体７０に格納された状態で提供され、ＢＭＵ６にインストールすることによりメモリ６３に格納される。また、通信網に接続されている図示しない外部コンピュータからＳＯＣ推定プログラム６３ａを取得し、メモリ６３に記憶させることにしてもよい。
ＣＰＵ６２はメモリ６３から読み出したＳＯＣ推定プログラム６３ａに従って、後述するＳＯＣ推定処理を実行する。
ＣＰＵ６２はメモリ６３から読み出したＳＯＣ推定プログラム６３ａに従って、後述するＳＯＣ推定処理を実行する。
電圧計測部８は、電圧検知線を介して電池３の両端に夫々接続されており、各電池３の電圧値を所定時間間隔で測定する。
電流計測部９は、電流センサ７を介して電池３に流れる電流値を所定時間間隔で計測する。
電池モジュール１の外部端子５，５は、エンジン始動用のスターターモータ及び電装品等の負荷１１に接続されている。
ＥＣＵ（Electronic Control Unit）１０は、ＢＭＵ６及び負荷１１に接続されている
。
電流計測部９は、電流センサ７を介して電池３に流れる電流値を所定時間間隔で計測する。
電池モジュール１の外部端子５，５は、エンジン始動用のスターターモータ及び電装品等の負荷１１に接続されている。
ＥＣＵ（Electronic Control Unit）１０は、ＢＭＵ６及び負荷１１に接続されている
。
以下、本実施形態に係るＳＯＣ推定方法について説明する。
図１６及び図１７は、ＣＰＵ６２によるＳＯＣ推定処理の手順を示すフローチャートである。ＣＰＵ６２は、所定の間隔でＳ１からの処理を繰り返す。
ＣＰＵ６２は、電池３の端子間の電圧値及び電流値を取得する（Ｓ１）。後述する下限電圧値及び到達電圧値はＯＣＶであるので、電池３の電流量が大きい場合、取得した電圧値をＯＣＶに補正する必要がある。ＯＣＶへの補正値は、複数の電圧値及び電流値のデータから回帰直線を用いて、電流値がゼロである場合の電圧値を推定すること等により得られる。電池３を流れる電流量が暗電流程度に小さい（微小電流である）場合、取得した電圧値をＯＣＶとみなすことができる。
図１６及び図１７は、ＣＰＵ６２によるＳＯＣ推定処理の手順を示すフローチャートである。ＣＰＵ６２は、所定の間隔でＳ１からの処理を繰り返す。
ＣＰＵ６２は、電池３の端子間の電圧値及び電流値を取得する（Ｓ１）。後述する下限電圧値及び到達電圧値はＯＣＶであるので、電池３の電流量が大きい場合、取得した電圧値をＯＣＶに補正する必要がある。ＯＣＶへの補正値は、複数の電圧値及び電流値のデータから回帰直線を用いて、電流値がゼロである場合の電圧値を推定すること等により得られる。電池３を流れる電流量が暗電流程度に小さい（微小電流である）場合、取得した電圧値をＯＣＶとみなすことができる。
ＣＰＵ６２は、電流値の絶対値が閾値以上であるか否かを判定する（Ｓ２）。閾値は、電池３の状態が充電状態又は放電状態と、休止状態とのいずれであるかを判定する為に設定される。ＣＰＵ６２は電流値の絶対値が閾値以上でないと判定した場合（Ｓ２：ＮＯ）、処理をＳ１３へ進める。
ＣＰＵ６２は、電流値の絶対値が閾値以上であると判定した場合（Ｓ２：ＹＥＳ）、電流値が０より大きいか否かを判定する（Ｓ３）。電流値が０より大きい場合、電池３の状態は充電状態であると判定できる。ＣＰＵ６２は電流値が０より大きくないと判定した場合（Ｓ３：ＮＯ）、処理をＳ９へ進める。
ＣＰＵ６２は電流値が０より大きいと判定した場合（Ｓ３：ＹＥＳ）、電圧値が下限電圧値以上であるか否かを判定する（Ｓ４）。ＣＰＵ６２は電圧値が下限電圧値以上でないと判定した場合（Ｓ４：ＮＯ）、処理をＳ８へ進める。
ＣＰＵ６２は電圧値が下限電圧値以上であると判定した場合（Ｓ４：ＹＥＳ）、電圧参照のフラグをオンにする（Ｓ５）。
ＣＰＵ６２は、取得した電圧値が前回の到達電圧値より大きいか否かを判定する（Ｓ６）。ＣＰＵ６２は電圧値が前回の到達電圧値より大きくないと判定した場合（Ｓ６：ＮＯ）、処理をＳ８へ進める。
ＣＰＵ６２は、取得した電圧値が前回の到達電圧値より大きいか否かを判定する（Ｓ６）。ＣＰＵ６２は電圧値が前回の到達電圧値より大きくないと判定した場合（Ｓ６：ＮＯ）、処理をＳ８へ進める。
ＣＰＵ６２は電圧値が前回の到達電圧値より大きいと判定した場合（Ｓ６：ＹＥＳ）、メモリ６３において、電圧値を到達電圧値に更新する（Ｓ７）。
ＣＰＵ６２は、電流積算によりＳＯＣを推定し（Ｓ８）、処理を終了する。
ＣＰＵ６２は、電流積算によりＳＯＣを推定し（Ｓ８）、処理を終了する。
ＣＰＵ６２は電流値が０より小さく、電池３の状態が放電状態であると判定した場合、Ｓ９において、電圧値が下限電圧値未満であるか否かを判定する（Ｓ９）。ＣＰＵ６２は電圧値が下限電圧値未満でないと判定した場合（Ｓ９：ＮＯ）、処理をＳ１２へ進める。
ＣＰＵ６２は電圧値が下限電圧値未満であると判定した場合（Ｓ９：ＹＥＳ）、電圧参照のフラグをオフする（Ｓ１０）。
ＣＰＵ６２は到達電圧値をリセットする（Ｓ１１）。
ＣＰＵ６２は、電流積算によりＳＯＣを推定し（Ｓ１２）、処理を終了する。
ＣＰＵ６２は電圧値が下限電圧値未満であると判定した場合（Ｓ９：ＹＥＳ）、電圧参照のフラグをオフする（Ｓ１０）。
ＣＰＵ６２は到達電圧値をリセットする（Ｓ１１）。
ＣＰＵ６２は、電流積算によりＳＯＣを推定し（Ｓ１２）、処理を終了する。
ＣＰＵ６２は電流値の絶対値が閾値未満であり、電池３の状態が休止状態であると判定した場合、電圧参照のフラグがオンであるか否かを判定する（Ｓ１３）。ＣＰＵ６２は電圧参照フラグがオンでないと判定した場合（Ｓ１３：ＮＯ）、処理をＳ１６へ進める。
ＣＰＵ６２は電圧参照のフラグがオンであると判定した場合（Ｓ１３：ＹＥＳ）、前回のＳ２で休止状態であると判定してから設定時間が経過したか否かを判定する（Ｓ１４）。設定時間は、取得した電圧値をＯＣＶとみなす為に十分な時間を予め実験により求める。休止状態であると判定してからの電流値の取得回数及び取得間隔に基づき、前記時間を超えたか否かを判定する。これにより、休止状態において、より高精度にＳＯＣを推定することができる。
ＣＰＵ６２は設定時間が経過していないと判定した場合（Ｓ１４：ＮＯ）、処理をＳ１８へ進める。
ＣＰＵ６２はＳ１８において、電流積算によりＳＯＣを推定し、処理を終了する。
ＣＰＵ６２は設定時間が経過していないと判定した場合（Ｓ１４：ＮＯ）、処理をＳ１８へ進める。
ＣＰＵ６２はＳ１８において、電流積算によりＳＯＣを推定し、処理を終了する。
ＣＰＵ６２は設定時間が経過したと判定した場合（Ｓ１４：ＹＥＳ）、取得した電圧値はＯＣＶとみなすことができる。
ＣＰＵ６２は、テーブル６３ｂから充放電ＳＯＣ−ＯＣＶ曲線を読み出す（Ｓ１５）。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線を生成する（Ｓ１６）。ＣＰＵ６２は、放電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第１点が、充電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第２点に重なるように放電ＳＯＣ−ＯＣＶ曲線を横軸に平行に移動させる。即ちΔＳＯＣ分、平行移動させる。そして、下限電圧値から到達電圧値まで曲線を切り取り、参照ＳＯＣ−ＯＣＶ曲線を生成する。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線において、Ｓ１で取得した電圧値に対応するＳＯＣを読み取ってＳＯＣを推定し（Ｓ１７）、処理を終了する。
充放電が繰り返された場合、電圧値が昇降し、即ち充電から放電へ切り換わった変曲点のうち高い変曲点が到達電圧値に設定される。
ＣＰＵ６２は、テーブル６３ｂから充放電ＳＯＣ−ＯＣＶ曲線を読み出す（Ｓ１５）。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線を生成する（Ｓ１６）。ＣＰＵ６２は、放電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第１点が、充電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第２点に重なるように放電ＳＯＣ−ＯＣＶ曲線を横軸に平行に移動させる。即ちΔＳＯＣ分、平行移動させる。そして、下限電圧値から到達電圧値まで曲線を切り取り、参照ＳＯＣ−ＯＣＶ曲線を生成する。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線において、Ｓ１で取得した電圧値に対応するＳＯＣを読み取ってＳＯＣを推定し（Ｓ１７）、処理を終了する。
充放電が繰り返された場合、電圧値が昇降し、即ち充電から放電へ切り換わった変曲点のうち高い変曲点が到達電圧値に設定される。
なお、ＣＰＵ６２が電圧計測部８から取得する電圧値は、電流値により多少変動するので、実験により補正係数を求めて電圧値を補正することもできる。
以上のように、本実施形態においては、ヒステリシスが小さい（ヒステリシスが実質的にない）、下限電圧値から到達電圧値までの範囲において、参照ＳＯＣ−ＯＣＶ曲線、及び現時点の電圧値に基づいてＳＯＣを推定するので、ＳＯＣの推定の精度が良好である。従って、高精度にＯＣＶリセットを行うことができる。
充電及び放電のいずれにおいても、ＳＯＣを推定することができる。電圧値の昇降の変曲点を到達電圧値に設定して参照ＳＯＣ−ＯＣＶ曲線を生成することで、複雑なパターンで充放電が繰り返された場合に、電圧値の履歴のみによってＳＯＣを推定できる。また、取得した電圧値が前回の到達電圧値を超えた場合のみ、到達電圧値を更新することで、充電時の最終の電圧値に基づきＳＯＣ−ＯＣＶ曲線を選択する特許文献１より高精度にＳＯＣを推定できる。
充電及び放電のいずれにおいても、ＳＯＣを推定することができる。電圧値の昇降の変曲点を到達電圧値に設定して参照ＳＯＣ−ＯＣＶ曲線を生成することで、複雑なパターンで充放電が繰り返された場合に、電圧値の履歴のみによってＳＯＣを推定できる。また、取得した電圧値が前回の到達電圧値を超えた場合のみ、到達電圧値を更新することで、充電時の最終の電圧値に基づきＳＯＣ−ＯＣＶ曲線を選択する特許文献１より高精度にＳＯＣを推定できる。
本実施の形態においては、メモリ６３に記憶された到達電圧値に基づき、放電ＳＯＣ−ＯＣＶ曲線をΔＳＯＣ分、平行移動させて参照ＳＯＣ−ＯＣＶ曲線を生成する。全ての到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を一つの放電ＳＯＣ−ＯＣＶ曲線から生成する。到達電圧値毎に実測により参照ＳＯＣ−ＯＣＶ曲線をメモリ６３に格納しておく必要がない。即ち下限電圧値から各到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を容易に取得することができる。実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、作業量が著しく軽減される。
実測した隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填する場合と比較して、ＳＯＣ推定の精度が向上する。
電池３の劣化に応じて、参照ＳＯＣ−ＯＣＶ曲線を取得する場合、実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、電池３の使用期間における作業量が著しく軽減される。
ＳＯＣ−ＯＣＶ曲線から推定できる為、蓄電量としてＳＯＣに限定されるものではなく、電力量等、電池３に蓄えられた現在のエネルギーの量を推定することができる。
実測した隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填する場合と比較して、ＳＯＣ推定の精度が向上する。
電池３の劣化に応じて、参照ＳＯＣ−ＯＣＶ曲線を取得する場合、実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、電池３の使用期間における作業量が著しく軽減される。
ＳＯＣ−ＯＣＶ曲線から推定できる為、蓄電量としてＳＯＣに限定されるものではなく、電力量等、電池３に蓄えられた現在のエネルギーの量を推定することができる。
（実施形態２）
実施形態２においては、リアルタイムにＳＯＣを推定する場合について説明する。ＣＰＵ６２によるＳＯＣ推定処理が異なること以外は、実施形態１と同様の構成を有する。
以下、本実施形態のＣＰＵ６２によるＳＯＣ推定処理について説明する。
実施形態２においては、リアルタイムにＳＯＣを推定する場合について説明する。ＣＰＵ６２によるＳＯＣ推定処理が異なること以外は、実施形態１と同様の構成を有する。
以下、本実施形態のＣＰＵ６２によるＳＯＣ推定処理について説明する。
図１８及び図１９は、ＣＰＵ６２によるＳＯＣ推定処理の手順を示すフローチャートである。ＣＰＵ６２は、所定の間隔でＳ２１からの処理を繰り返す。
ＣＰＵ６２は、電池３の端子間の電圧値及び電流値を取得する（Ｓ２１）。
ＣＰＵ６２は、電流値の絶対値が閾値以上であるか否かを判定する（Ｓ２２）。閾値は、電池３の状態が充電状態又は放電状態と、休止状態とのいずれであるかを判定する為に設定される。ＣＰＵ６２は電流値の絶対値が閾値以上でないと判定した場合（Ｓ２２：ＮＯ）、処理をＳ３３へ進める。
ＣＰＵ６２は、電池３の端子間の電圧値及び電流値を取得する（Ｓ２１）。
ＣＰＵ６２は、電流値の絶対値が閾値以上であるか否かを判定する（Ｓ２２）。閾値は、電池３の状態が充電状態又は放電状態と、休止状態とのいずれであるかを判定する為に設定される。ＣＰＵ６２は電流値の絶対値が閾値以上でないと判定した場合（Ｓ２２：ＮＯ）、処理をＳ３３へ進める。
ＣＰＵ６２は、電流値の絶対値が閾値以上であると判定した場合（Ｓ２２：ＹＥＳ）、電流値が０より大きいか否かを判定する（Ｓ２３）。電流値が０より大きい場合、電池３の状態は充電状態である。ＣＰＵ６２は電流値が０より大きくないと判定した場合（Ｓ２３：ＮＯ）、処理をＳ２９へ進める。
ＣＰＵ６２は電流値が０より大きいと判定した場合（Ｓ２３：ＹＥＳ）、電圧値が前回の到達電圧値より大きいか否かを判定する（Ｓ２４）。前記電圧値は、例えば複数の電圧値及び電流値のデータから回帰直線を用い、電流値がゼロである場合の値を推定することで補正してもよい。ＣＰＵ６２は電圧値が前回の到達電圧値より大きくないと判定した場合（Ｓ２４：ＮＯ）、処理をＳ２６へ進める。
ＣＰＵ６２は電流値が０より大きいと判定した場合（Ｓ２３：ＹＥＳ）、電圧値が前回の到達電圧値より大きいか否かを判定する（Ｓ２４）。前記電圧値は、例えば複数の電圧値及び電流値のデータから回帰直線を用い、電流値がゼロである場合の値を推定することで補正してもよい。ＣＰＵ６２は電圧値が前回の到達電圧値より大きくないと判定した場合（Ｓ２４：ＮＯ）、処理をＳ２６へ進める。
ＣＰＵ６２は電圧値が前回の到達電圧値より大きいと判定した場合（Ｓ２４：ＹＥＳ）、電圧値を到達電圧値に更新する（Ｓ２５）。
ＣＰＵ６２は、電圧値が下限電圧値以上であるか否かを判定する（Ｓ２６）。ＣＰＵ６２は電圧値が下限電圧値以上でないと判定した場合（Ｓ２６：ＮＯ）、電流積算によりＳＯＣを推定し（Ｓ３０）、処理を終了する。
ＣＰＵ６２は、電圧値が下限電圧値以上であるか否かを判定する（Ｓ２６）。ＣＰＵ６２は電圧値が下限電圧値以上でないと判定した場合（Ｓ２６：ＮＯ）、電流積算によりＳＯＣを推定し（Ｓ３０）、処理を終了する。
ＣＰＵ６２は電圧値が下限電圧値以上であると判定した場合（Ｓ２６：ＹＥＳ）、ＣＰＵ６２は、テーブル６３ｂから充放電ＳＯＣ−ＯＣＶ曲線を読み出す（Ｓ２７）。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線を生成する（Ｓ２８）。ＣＰＵ６２は、放電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第１点が、充電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第２点に重なるように放電ＳＯＣ−ＯＣＶ曲線を横軸に平行に移動させる。即ちΔＳＯＣ分、平行移動させる。そして、下限電圧値から到達電圧値まで曲線を切り取り、参照ＳＯＣ−ＯＣＶ曲線を生成する。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線において、現時点のＯＣＶからＳＯＣを読み取ってＳＯＣを推定し（Ｓ２９）、処理を終了する。ＣＰＵ６２は、Ｓ２１で取得した電圧値及び電流値から現時点のＯＣＶを算出する。ＯＣＶの算出は、複数の電圧値及び電流値のデータから回帰直線を用いて、電流値がゼロである場合の電圧値を推定すること等により得られる。また、電流値が暗電流の電流値のように小さい場合は、取得した電圧値をＯＣＶに読み替えることもできる。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線を生成する（Ｓ２８）。ＣＰＵ６２は、放電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第１点が、充電ＳＯＣ−ＯＣＶ曲線上の到達電圧値に対応する第２点に重なるように放電ＳＯＣ−ＯＣＶ曲線を横軸に平行に移動させる。即ちΔＳＯＣ分、平行移動させる。そして、下限電圧値から到達電圧値まで曲線を切り取り、参照ＳＯＣ−ＯＣＶ曲線を生成する。
ＣＰＵ６２は、参照ＳＯＣ−ＯＣＶ曲線において、現時点のＯＣＶからＳＯＣを読み取ってＳＯＣを推定し（Ｓ２９）、処理を終了する。ＣＰＵ６２は、Ｓ２１で取得した電圧値及び電流値から現時点のＯＣＶを算出する。ＯＣＶの算出は、複数の電圧値及び電流値のデータから回帰直線を用いて、電流値がゼロである場合の電圧値を推定すること等により得られる。また、電流値が暗電流の電流値のように小さい場合は、取得した電圧値をＯＣＶに読み替えることもできる。
ＣＰＵ６２は電流値が０より小さく、電池３の状態が放電状態であると判定した場合、Ｓ２９において、電圧値が下限電圧値以上であるか否かを判定する（Ｓ３１）。
ＣＰＵ６２は電圧値が下限電圧値以上であると判定した場合（Ｓ３１：ＹＥＳ）、上記と同様にして電圧参照によりＳＯＣを推定する（Ｓ３２）。
ＣＰＵ６２は電圧値が下限電圧値以上でないと判定した場合（Ｓ３１：ＮＯ）、到達電圧値をリセットする（Ｓ３３）。
ＣＰＵ６２は電流積算によりＳＯＣを推定し（Ｓ３４）、処理を終了する。
ＣＰＵ６２は電圧値が下限電圧値以上であると判定した場合（Ｓ３１：ＹＥＳ）、上記と同様にして電圧参照によりＳＯＣを推定する（Ｓ３２）。
ＣＰＵ６２は電圧値が下限電圧値以上でないと判定した場合（Ｓ３１：ＮＯ）、到達電圧値をリセットする（Ｓ３３）。
ＣＰＵ６２は電流積算によりＳＯＣを推定し（Ｓ３４）、処理を終了する。
ＣＰＵ６２は電流値の絶対値が閾値未満であり、電池３の状態が休止状態であると判定した場合、電圧値が下限電圧値以上であるか否かを判定する（Ｓ３５）。ＣＰＵ６２は電圧値が下限電圧値以上でないと判定した場合（Ｓ３６：ＮＯ）、処理をＳ３６へ進める。
ＣＰＵ６２は電圧値が下限電圧値以上であると判定した場合（Ｓ３６：ＹＥＳ）、前回、Ｓ２２で休止状態であると判定してから設定時間が経過したか否かを判定する（Ｓ３６）。設定時間は、取得した電圧値をＯＣＶとみなす為に十分な時間を予め実験により求める。
ＣＰＵ６２は設定時間が経過していないと判定した場合（Ｓ３６：ＮＯ）、処理をＳ３８へ進める。
ＣＰＵ６２は電流積算によりＳＯＣを推定し（Ｓ３８）、処理を終了する。
ＣＰＵ６２は電圧値が下限電圧値以上であると判定した場合（Ｓ３６：ＹＥＳ）、前回、Ｓ２２で休止状態であると判定してから設定時間が経過したか否かを判定する（Ｓ３６）。設定時間は、取得した電圧値をＯＣＶとみなす為に十分な時間を予め実験により求める。
ＣＰＵ６２は設定時間が経過していないと判定した場合（Ｓ３６：ＮＯ）、処理をＳ３８へ進める。
ＣＰＵ６２は電流積算によりＳＯＣを推定し（Ｓ３８）、処理を終了する。
ＣＰＵ６２は設定時間が経過したと判定した場合（Ｓ３６：ＹＥＳ）、取得した電圧値はＯＣＶとみなすことができ、上記と同様にして電圧参照によりＳＯＣを推定し（Ｓ３７）、処理を終了する。
本実施形態においては、充放電時にリアルタイムにＳＯＣを推定することができる。
ヒステリシスが実質的にない、下限電圧値から到達電圧値までの範囲において、参照ＳＯＣ−ＯＣＶ曲線、及び現時点の電圧値に基づいてＳＯＣを推定する。従って、ＳＯＣの推定の精度が良好である。
メモリ６３に記憶された到達電圧値に基づき、テーブル６３ｂに記憶された放電ＳＯＣ−ＯＣＶ曲線をΔＳＯＣ分、平行移動させて参照ＳＯＣ−ＯＣＶ曲線を生成する。全ての到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を一つの放電ＳＯＣ−ＯＣＶ曲線から生成する。到達電圧値毎に実測により参照ＳＯＣ−ＯＣＶ曲線をメモリ６３に格納しておく必要がない。即ち下限電圧値から各到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を容易に取得することができる。実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、作業量が著しく軽減される。
実測した隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填する場合と比較して、ＳＯＣ推定の精度が向上する。
電池３の劣化に応じて、参照ＳＯＣ−ＯＣＶ曲線を取得する場合、実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、電池３の使用期間における作業量が著しく軽減される。
ＳＯＣ−ＯＣＶ曲線から推定できる為、蓄電量としてＳＯＣに限定されるものではなく、電力量等、電池３に蓄えられた現在のエネルギーの量を推定することができる。
充電及び放電のいずれにおいても、ＳＯＣを推定できる。複雑なパターンで、充放電が繰り返された場合においても、電圧値の履歴のみにより、ＳＯＣを推定できる。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、電池３に蓄えられた現在のエネルギーの量を推定することができる。
ヒステリシスが実質的にない、下限電圧値から到達電圧値までの範囲において、参照ＳＯＣ−ＯＣＶ曲線、及び現時点の電圧値に基づいてＳＯＣを推定する。従って、ＳＯＣの推定の精度が良好である。
メモリ６３に記憶された到達電圧値に基づき、テーブル６３ｂに記憶された放電ＳＯＣ−ＯＣＶ曲線をΔＳＯＣ分、平行移動させて参照ＳＯＣ−ＯＣＶ曲線を生成する。全ての到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を一つの放電ＳＯＣ−ＯＣＶ曲線から生成する。到達電圧値毎に実測により参照ＳＯＣ−ＯＣＶ曲線をメモリ６３に格納しておく必要がない。即ち下限電圧値から各到達電圧値に対応する参照ＳＯＣ−ＯＣＶ曲線を容易に取得することができる。実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、作業量が著しく軽減される。
実測した隣り合うＳＯＣ−ＯＣＶ曲線間に位置すべき曲線を内挿計算により求めて補填する場合と比較して、ＳＯＣ推定の精度が向上する。
電池３の劣化に応じて、参照ＳＯＣ−ＯＣＶ曲線を取得する場合、実測するのはフルの充放電ＳＯＣ−ＯＣＶ曲線のみであり、電池３の使用期間における作業量が著しく軽減される。
ＳＯＣ−ＯＣＶ曲線から推定できる為、蓄電量としてＳＯＣに限定されるものではなく、電力量等、電池３に蓄えられた現在のエネルギーの量を推定することができる。
充電及び放電のいずれにおいても、ＳＯＣを推定できる。複雑なパターンで、充放電が繰り返された場合においても、電圧値の履歴のみにより、ＳＯＣを推定できる。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、電池３に蓄えられた現在のエネルギーの量を推定することができる。
以上のように、蓄電量推定装置は、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定する蓄電量推定装置であって、反応の比率が、他の領域の反応の比率に対して、一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性に基づいて得られる電圧参照蓄電量−電圧値特性を用いて、蓄電量を推定する推定部を備える。
上記構成においては、充電と放電とで、蓄電量に対する電圧値の変化が略一致している、反応の比率が、他の領域の反応の比率に対して、一の電気化学反応側にある領域で、電圧参照蓄電量−電圧値特性に基づいて蓄電量を推定する。
前記一の電気化学反応側にある領域において、放電反応としては、一の電気化学反応のみが生じている。この領域における一の電気化学反応の放電電気量は、一の電気化学反応の充電電気量とみなすことができる。即ち、充電状態、放電状態で蓄電量を推定するときに、同一の電圧参照蓄電量−電圧値特性を使用できる。
蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得する。一の電気化学反応が生じる前記蓄電量−電圧値放電特性を用いることにより、容易に電圧参照蓄電量−電圧値特性を取得できる。
高容量で、蓄電量−電圧値特性がヒステリシスを示す活物質を有する蓄電素子の蓄電量を良好に容易に推定できる。
前記一の電気化学反応側にある領域において、放電反応としては、一の電気化学反応のみが生じている。この領域における一の電気化学反応の放電電気量は、一の電気化学反応の充電電気量とみなすことができる。即ち、充電状態、放電状態で蓄電量を推定するときに、同一の電圧参照蓄電量−電圧値特性を使用できる。
蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得する。一の電気化学反応が生じる前記蓄電量−電圧値放電特性を用いることにより、容易に電圧参照蓄電量−電圧値特性を取得できる。
高容量で、蓄電量−電圧値特性がヒステリシスを示す活物質を有する蓄電素子の蓄電量を良好に容易に推定できる。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量を推定できる。充放電曲線に基づいて、ＳＯＣ０％までの放電可能なエネルギー、及びＳＯＣ１００％までに必要な充電エネルギーを予測することができる。
上述の蓄電量推定装置において、前記一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性を、同一電圧値における前記蓄電量−電圧値充電特性及び前記蓄電量−電圧値放電特性間の蓄電量の差分、前記蓄電量−電圧値充電特性に向けて移動させて得られる電圧参照蓄電量−電圧値特性に基づいて、蓄電量を推定してもよい。
蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得する。同一電圧値における両特性の蓄電量の差（Δ蓄電量）は、他の電気化学反応に対応する。一の電気化学反応と他の電気化学反応とは、実質的に独立して生じる。他の電気化学反応の反応量はΔ蓄電量に対応し、一の電気化学反応に影響を与えない。
上記構成においては、同一電圧値における前記蓄電量−電圧値充電特性及び前記蓄電量−電圧値放電特性間の蓄電量の差分、蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させて電圧参照蓄電量−電圧値特性を得る。即ち、各電圧値を起点とする電圧参照蓄電量−電圧値特性は、同一の一の電気化学反応に基づく、同一の蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させ、各起点の電圧値を上端として切り取った形状を有する。
上記構成によれば、蓄電量の推定の精度がより良好である。
上記構成においては、同一電圧値における前記蓄電量−電圧値充電特性及び前記蓄電量−電圧値放電特性間の蓄電量の差分、蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させて電圧参照蓄電量−電圧値特性を得る。即ち、各電圧値を起点とする電圧参照蓄電量−電圧値特性は、同一の一の電気化学反応に基づく、同一の蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させ、各起点の電圧値を上端として切り取った形状を有する。
上記構成によれば、蓄電量の推定の精度がより良好である。
上述の蓄電量推定装置において、前記蓄電素子の電圧値を取得する電圧取得部と、該電圧取得部が取得した電圧値が、ヒステリシスの有無が切り替わる下限の電圧値よりも高くなった後の電圧値に基づいて到達電圧値を設定する設定部とを備え、前記同一電圧値は、前記到達電圧値であることが好ましい。
上記構成においては、電圧値の昇降に基づき到達電圧値を設定して電圧参照蓄電量−電圧値特性を生成する。複雑なパターンで、充放電が繰り返された場合においても、電圧値の履歴のみにより電圧参照蓄電量−電圧値特性を生成して、蓄電量を推定できる。
上述の蓄電量推定装置において、前記推定部は、前記電圧取得部により取得した電圧値が前記下限の電圧値から前記到達電圧値までの間にあるときに、前記移動蓄電量−電圧値特性を参照して前記蓄電量を推定することが好ましい。
一の電気化学反応が主として生じ、ヒステリシスが小さい、下限電圧値から到達電圧値までの範囲において、他の電気化学反応の反応量は少ない。充放電反応において、一の電気化学反応が主に生じている。上記構成においては、この場合に、電圧参照蓄電量−電圧値特性、及び取得した電圧値に基づいて蓄電量を推定する。従って、蓄電量の推定の精度が良好である。
上述の蓄電量推定装置において、前記設定部は、前記到達電圧値を記憶部に記憶し、前記電圧取得部が取得した電圧値が前記記憶部に前回記憶された到達電圧値より大きい場合に、取得した電圧値を到達電圧値に更新することが好ましい。
更新された到達電圧値に基づいて電圧参照蓄電量−電圧値特性を生成することで、精度良く蓄電量を推定することができる。
上述の蓄電量推定装置において、前記蓄電量はＳＯＣであることが好ましい。
高容量材料についてＳＯＣを推定することで、既存の制御システムへの適用性が向上する。ＳＯＣに基づいて、放電可能エネルギーのような蓄電量を容易に算出することができる。蓄電量推定装置は、ＯＣＶとＳＯＣとが一対一対応しない、ヒステリシスを有する電極材料を用いた蓄電素子の充電状態を、特別なセンサや追加の部品を要することなく、高精度に推定することができる。
蓄電モジュールは、複数の蓄電素子と、上述のいずれかの蓄電量推定装置とを備える。
車両用の蓄電モジュールや産業用の蓄電モジュールは、典型的には複数の蓄電素子が直列に接続される。複数の蓄電素子が直列及び並列に接続される場合もある。蓄電モジュールの性能を発揮するためには、各蓄電素子の蓄電量を精度良く推定し、複数の蓄電素子間で蓄電量のばらつきが生じた時にはバランシング処理を行う必要がある。たとえ各蓄電素子が高容量であっても、複数の蓄電素子間の蓄電量のばらつきを検出できなければ、蓄電モジュールの性能を使い切ることができない。上述の蓄電量推定装置によって各蓄電素子の蓄電量を高精度に推定できるので、蓄電モジュールの性能を最大限発揮することが出来る。蓄電モジュールは、特に高容量に対する要求が高い、ＥＶやＰＨＥＶの動力源として好適に使用される。
本発明に係る蓄電量推定方法は、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定する蓄電量推定方法であって、反応の比率が、他の領域の反応の比率に対して、一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性に基づいて得られる電圧参照蓄電量−電圧値特性を取得し、該電圧参照蓄電量−電圧値特性に基づいて、蓄電量を推定する。
上記方法においては、充電と放電とで、蓄電量に対する電圧値の変化が略一致している、一の電気化学反応が多く生じる場合に、電圧参照蓄電量−電圧値特性に基づいて蓄電量を推定する。
一の電気化学反応が多く生じる領域において、放電反応としては、一の電気化学反応のみが生じている。この領域における一の電気化学反応の放電電気量は、一の電気化学反応の充電電気量とみなすことができる。即ち、充電状態、放電状態で蓄電量を推定するときに、同一の電圧参照蓄電量−電圧値特性を使用できる。
蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得する。一の電気化学反応が生じる前記蓄電量−電圧値放電特性を用いることにより、容易に電圧参照蓄電量−電圧値特性を取得できる。
実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、作業量は少ない。
隣り合う電圧参照蓄電量−電圧値特性間に位置すべき電圧参照蓄電量−電圧値特性を内挿計算により求めて補填する場合と比較して、各到達電圧値に対応して電圧参照蓄電量−電圧値特性参照を生成するので、蓄電量推定の精度が向上する。
蓄電素子の劣化に応じて、電圧参照蓄電量−電圧値特性を取得する場合、実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、蓄電素子の使用期間における作業量は少ない。
高容量で、蓄電量−電圧値特性がヒステリシスを示す活物質を有する蓄電素子の蓄電量を良好に容易に推定できる。
一の電気化学反応が多く生じる領域において、放電反応としては、一の電気化学反応のみが生じている。この領域における一の電気化学反応の放電電気量は、一の電気化学反応の充電電気量とみなすことができる。即ち、充電状態、放電状態で蓄電量を推定するときに、同一の電圧参照蓄電量−電圧値特性を使用できる。
蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得する。一の電気化学反応が生じる前記蓄電量−電圧値放電特性を用いることにより、容易に電圧参照蓄電量−電圧値特性を取得できる。
実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、作業量は少ない。
隣り合う電圧参照蓄電量−電圧値特性間に位置すべき電圧参照蓄電量−電圧値特性を内挿計算により求めて補填する場合と比較して、各到達電圧値に対応して電圧参照蓄電量−電圧値特性参照を生成するので、蓄電量推定の精度が向上する。
蓄電素子の劣化に応じて、電圧参照蓄電量−電圧値特性を取得する場合、実測するのは蓄電量−電圧値放電特性及び蓄電量−電圧値充電特性のみであり、蓄電素子の使用期間における作業量は少ない。
高容量で、蓄電量−電圧値特性がヒステリシスを示す活物質を有する蓄電素子の蓄電量を良好に容易に推定できる。
充電及び放電のいずれにおいても、蓄電量を推定することができる。電圧値の昇降に係る変曲点を到達電圧値に設定して蓄電量−電圧値特性を選択する。複雑なパターンで、充放電が繰り返された場合においても、電圧値の履歴のみにより、蓄電量を推定できる。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量を推定できる。
電圧値を用いることができるので、蓄電量としてＳＯＣに限定されず、電力量等、蓄電素子に蓄えられた現在のエネルギーの量を推定できる。
上述の蓄電量推定方法は、取得した電圧値が、ヒステリシスの有無が切り替わる下限の電圧値よりも高くなった後の電圧値に基づいて到達電圧値を設定し、該到達電圧値における、前記蓄電量−電圧値放電特性及び前記蓄電量−電圧値充電特性間の蓄電量の差分、前記蓄電量−電圧値放電特性を、前記蓄電量−電圧値充電特性側に向けて移動させて前記電圧参照蓄電量−電圧値特性を取得し、取得した電圧値が前記下限の電圧値から前記到達電圧値までの間にあるときに、前記電圧参照蓄電量−電圧値特性を参照して前記蓄電量を推定することが好ましい。
蓄電量−電圧値充電特性、及び蓄電量−電圧値放電特性を取得する。同一電圧値における両特性の蓄電量の差（Δ蓄電量）は、他の電気化学反応に対応する。他の電気化学反応の反応量はΔ蓄電量に対応し、一の電気化学反応に影響を与えない。
そこで、一の電気化学反応が生じる蓄電量−電圧値放電特性を基準にしてΔ蓄電量で補正することにより、容易に電圧参照蓄電量−電圧値特性を取得できる。即ち、各電圧値を起点とする電圧参照蓄電量−電圧値特性は、同一の一の電気化学反応に基づく、同一の蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させ、各起点の電圧値を上端として切り取った形状を有する。
上記方法においては、電圧値の昇降に基づき到達電圧値を設定して電圧参照蓄電量−電圧値特性を生成する。複雑なパターンで、充放電が繰り返された場合においても、電圧値の履歴のみにより電圧参照蓄電量−電圧値特性を生成して、蓄電量を推定できる。
そして、一の電気化学反応が生じ、ヒステリシスが小さい、下限電圧値から到達電圧値までの範囲において、電圧参照蓄電量−電圧値特性、及び現時点の電圧値に基づいて蓄電量を推定する。従って、蓄電量の推定の精度が良好である。
そこで、一の電気化学反応が生じる蓄電量−電圧値放電特性を基準にしてΔ蓄電量で補正することにより、容易に電圧参照蓄電量−電圧値特性を取得できる。即ち、各電圧値を起点とする電圧参照蓄電量−電圧値特性は、同一の一の電気化学反応に基づく、同一の蓄電量−電圧値放電特性を蓄電量−電圧値充電特性に向けて移動させ、各起点の電圧値を上端として切り取った形状を有する。
上記方法においては、電圧値の昇降に基づき到達電圧値を設定して電圧参照蓄電量−電圧値特性を生成する。複雑なパターンで、充放電が繰り返された場合においても、電圧値の履歴のみにより電圧参照蓄電量−電圧値特性を生成して、蓄電量を推定できる。
そして、一の電気化学反応が生じ、ヒステリシスが小さい、下限電圧値から到達電圧値までの範囲において、電圧参照蓄電量−電圧値特性、及び現時点の電圧値に基づいて蓄電量を推定する。従って、蓄電量の推定の精度が良好である。
コンピュータプログラムは、充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定するコンピュータに、前記蓄電量−電圧値放電特性、及び前記蓄電量−電圧値充電特性を取得し、取得した電圧値が、ヒステリシスの有無が切り替わる下限の電圧値よりも高くなったか否かを判定し、前記電圧値が前記下限の電圧値よりも高くなったと判定した場合に、到達電圧値を設定し、設定した前記到達電圧値に対応する、前記蓄電量−電圧値放電特性上の第１点を取得し、前記到達電圧値に対応する、前記蓄電量−電圧値充電特性上の第２点を取得し、前記第１点が前記第２点に重なるように、前記蓄電量−電圧値放電特性を前記蓄電量−電圧値充電特性に向けて移動させて電圧参照蓄電量−電圧値特性を取得し、取得した電圧値が前記下限の電圧値から前記到達電圧値までの間にあるか否かを判定し、前記電圧値が前記下限の電圧値から前記到達電圧値までの間にあると判定した場合に、前記電圧参照蓄電量−電圧値特性を参照して、前記電圧値における前記蓄電量を推定する処理を実行させる。
本発明は上述した実施の形態の内容に限定されるものではなく、請求項に示した範囲で種々の変更が可能である。即ち、請求項に示した範囲で適宜変更した技術的手段を組み合わせて得られる実施形態も本発明の技術的範囲に含まれる。
本発明に係る蓄電量推定装置は、車載用のリチウムイオン二次電池に適用される場合に限定されず、鉄道用回生電力貯蔵装置、太陽光発電システム等の他の蓄電モジュールにも適用できる。微小電流が流れる蓄電モジュールにおいては、蓄電素子の正極端子・負極端子間の電圧値、又は、蓄電モジュールの正極端子・負極端子間の電圧値をＯＣＶとみなすことができる。
蓄電素子は、リチウムイオン二次電池に限定されるものではなく、ヒステリシス特性を有する他の二次電池や電気化学セルであってもよい。
本発明に係る蓄電量推定装置は、車載用のリチウムイオン二次電池に適用される場合に限定されず、鉄道用回生電力貯蔵装置、太陽光発電システム等の他の蓄電モジュールにも適用できる。微小電流が流れる蓄電モジュールにおいては、蓄電素子の正極端子・負極端子間の電圧値、又は、蓄電モジュールの正極端子・負極端子間の電圧値をＯＣＶとみなすことができる。
蓄電素子は、リチウムイオン二次電池に限定されるものではなく、ヒステリシス特性を有する他の二次電池や電気化学セルであってもよい。
監視装置１００又はＢＭＵ６が蓄電量推定装置である場合に限定されない。ＣＭＵ（Cell Monitoring Unit）が蓄電量推定装置であってもよい。蓄電量推定装置は、監視装置１００等が組み込まれた蓄電モジュールの一部であってもよい。蓄電量推定装置は、蓄電素子や蓄電モジュールとは別個に構成されて、蓄熱量推定対象の蓄電素子を含む蓄電モジュールに、蓄熱量の推定時に接続されてもよい。蓄熱量推定装置は、蓄電素子や蓄電モジュールを遠隔監視してもよい。
１、５０ 電池モジュール（蓄電モジュール）
２ ケース
２１ ケース本体
２２ 蓋部
２３ ＢＭＵ収容部
２４ カバー
２５ 中蓋
２６ 仕切り板
３、２００ 電池（蓄電素子）
３１ ケース
３２ 端子
３３ 電極体
４ バスバー
５ 外部端子
６ ＢＭＵ（蓄電量推定装置）
６０ 情報処理部
６２ ＣＰＵ（推定部、電圧取得部、設定部、選択部）
６３ メモリ（記憶部）
６３ａ ＳＯＣ推定プログラム
６３ｂ テーブル
７ 電流センサ
８ 電圧計測部
９ 電流計測部
１０ ＥＣＵ
７０ 記録媒体
１００ 監視装置（蓄電量推定装置）
３００ 収容ケース
２ ケース
２１ ケース本体
２２ 蓋部
２３ ＢＭＵ収容部
２４ カバー
２５ 中蓋
２６ 仕切り板
３、２００ 電池（蓄電素子）
３１ ケース
３２ 端子
３３ 電極体
４ バスバー
５ 外部端子
６ ＢＭＵ（蓄電量推定装置）
６０ 情報処理部
６２ ＣＰＵ（推定部、電圧取得部、設定部、選択部）
６３ メモリ（記憶部）
６３ａ ＳＯＣ推定プログラム
６３ｂ テーブル
７ 電流センサ
８ 電圧計測部
９ 電流計測部
１０ ＥＣＵ
７０ 記録媒体
１００ 監視装置（蓄電量推定装置）
３００ 収容ケース
Claims (9)
- 充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定する蓄電量推定装置であって、
反応の比率が、他の領域の反応の比率に対して、一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性に基づいて得られる電圧参照蓄電量−電圧値特性を用いて、蓄電量を推定する推定部を備える、蓄電量推定装置。 - 前記一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性を、同一電圧値における前記蓄電量−電圧値充電特性及び前記蓄電量−電圧値放電特性間の蓄電量の差分、前記蓄電量−電圧値充電特性に向けて移動させて得られる電圧参照蓄電量−電圧値特性に基づいて、蓄電量を推定する、請求項１に記載の蓄電量推定装置。
- 前記蓄電素子の電圧値を取得する電圧取得部と、
該電圧取得部が取得した電圧値が、ヒステリシスの有無が切り替わる下限の電圧値よりも高くなった後の電圧値に基づいて到達電圧値を設定する設定部と
を備え、
前記同一電圧値は、前記到達電圧値である、請求項２に記載の蓄電量推定装置。 - 前記推定部は、前記電圧取得部により取得した電圧値が前記下限の電圧値から前記到達電圧値までの間にあるときに、前記電圧参照蓄電量−電圧値特性を参照して前記蓄電量を推定する、請求項３に記載の蓄電量推定装置。
- 前記設定部は、
前記到達電圧値を記憶部に記憶し、
前記電圧取得部が取得した電圧値が前記記憶部に前回記憶された到達電圧値より大きい場合に、取得した電圧値を到達電圧値に更新する、請求項３又は４に記載の蓄電量推定装置。 - 蓄電素子と、
請求項１から５までのいずれか１項に記載の蓄電量推定装置と
を備える、蓄電素子モジュール。 - 充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性間のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定する蓄電量推定方法であって、
反応の比率が、他の領域の反応の比率に対して、一の電気化学反応側にある領域で、前記蓄電量−電圧値放電特性に基づいて得られる電圧参照蓄電量−電圧値特性を取得し、
該電圧参照蓄電量−電圧値特性に基づいて、蓄電量を推定する、蓄電量推定方法。 - 取得した電圧値が、ヒステリシスの有無が切り替わる下限の電圧値よりも高くなった後の電圧値に基づいて到達電圧値を設定し、
該到達電圧値における、前記蓄電量−電圧値放電特性及び前記蓄電量−電圧値充電特性間の蓄電量の差分、前記蓄電量−電圧値放電特性を、前記蓄電量−電圧値充電特性側に向けて移動させて前記電圧参照蓄電量−電圧値特性を取得し、
取得した電圧値が前記下限の電圧値から前記到達電圧値までの間にあるときに、前記電圧参照蓄電量−電圧値特性を参照して前記蓄電量を推定する、請求項７に記載の蓄電量推定方法。 - 充放電の推移に応じて２以上の電気化学反応を生じ、一の電気化学反応が生じる場合に示される、充電の推移に応じた蓄電量−電圧値充電特性、及び放電の推移に応じた蓄電量−電圧値放電特性のヒステリシスが、他の電気化学反応が生じる場合の前記ヒステリシスより小さい、活物質を正極及び負極の少なくとも一方に含む蓄電素子の蓄電量を推定するコンピュータに、
前記蓄電量−電圧値放電特性、及び前記蓄電量−電圧値充電特性を取得し、
取得した電圧値が、ヒステリシスの有無が切り替わる下限の電圧値よりも高くなったか否かを判定し、
前記電圧値が前記下限の電圧値よりも高くなったと判定した場合に、到達電圧値を設定し、
設定した前記到達電圧値に対応する、前記蓄電量−電圧値放電特性上の第１点を取得し、
前記到達電圧値に対応する、前記蓄電量−電圧値充電特性上の第２点を取得し、
前記第１点が前記第２点に重なるように、前記蓄電量−電圧値放電特性を前記蓄電量−電圧値充電特性に向けて移動させて電圧参照蓄電量−電圧値特性を取得し、
取得した電圧値が前記下限の電圧値から前記到達電圧値までの間にあるか否かを判定し、
前記電圧値が前記下限の電圧値から前記到達電圧値までの間にあると判定した場合に、前記電圧参照蓄電量−電圧値特性を参照して、前記電圧値における前記蓄電量を推定する
処理を実行させる、コンピュータプログラム。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/JP2018/013053 WO2018181620A1 (ja) | 2017-03-29 | 2018-03-28 | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム |
| US16/498,280 US11181584B2 (en) | 2017-03-29 | 2018-03-28 | Storage amount estimation device, energy storage module, storage amount estimation method, and computer program |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017065005 | 2017-03-29 | ||
| JP2017065005 | 2017-03-29 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP6406468B1 true JP6406468B1 (ja) | 2018-10-17 |
| JP2018169398A JP2018169398A (ja) | 2018-11-01 |
Family
ID=63855243
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018062290A Active JP6406468B1 (ja) | 2017-03-29 | 2018-03-28 | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US11181584B2 (ja) |
| EP (1) | EP3605122A4 (ja) |
| JP (1) | JP6406468B1 (ja) |
| CN (1) | CN110462413B (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11237214B2 (en) | 2017-07-19 | 2022-02-01 | Gs Yuasa International Ltd. | Estimation device, energy storage apparatus, estimation method, and computer program |
| JP6477959B2 (ja) * | 2017-07-19 | 2019-03-06 | 株式会社Ｇｓユアサ | 推定装置、蓄電装置、推定方法、及びコンピュータプログラム |
| JP6477957B2 (ja) * | 2017-07-19 | 2019-03-06 | 株式会社Ｇｓユアサ | 推定装置、蓄電装置、推定方法、及びコンピュータプログラム |
| KR102458526B1 (ko) * | 2018-02-07 | 2022-10-25 | 주식회사 엘지에너지솔루션 | 배터리의 동작 상태에 따라 soc를 추정하는 장치 및 방법 |
| JP7040284B2 (ja) * | 2018-05-23 | 2022-03-23 | トヨタ自動車株式会社 | 二次電池の劣化状態推定方法、劣化状態推定装置、制御方法、及び制御システム |
| CN111856178B (zh) * | 2020-03-31 | 2021-10-08 | 同济大学 | 一种基于锂离子电容器电化学特征的soc分区估计方法 |
| KR102493221B1 (ko) * | 2020-12-04 | 2023-01-30 | 주식회사 에이치에너지 | 적응적 SoC 추정을 통한 ESS 최적 운전 시스템 및 방법 |
Family Cites Families (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4292721B2 (ja) * | 2001-02-14 | 2009-07-08 | 株式会社日本自動車部品総合研究所 | ハイブリッド車の電池状態制御方法 |
| US7570024B2 (en) | 2004-04-06 | 2009-08-04 | Cobasys, Llc | Battery state of charge voltage hysteresis estimator |
| US7795874B1 (en) * | 2007-10-10 | 2010-09-14 | Commissariat A L'energie Atomique | Method of estimation of the state of charge of a lead-acid battery |
| US7957921B2 (en) * | 2008-02-19 | 2011-06-07 | GM Global Technology Operations LLC | Model-based estimation of battery hysteresis |
| JP5866987B2 (ja) * | 2011-11-10 | 2016-02-24 | 日産自動車株式会社 | 二次電池の制御装置およびｓｏｃ検出方法 |
| JP2013158087A (ja) * | 2012-01-27 | 2013-08-15 | Toyota Motor Corp | 蓄電システム及び充電状態推定方法 |
| JP6065561B2 (ja) * | 2012-03-08 | 2017-01-25 | 日産自動車株式会社 | 二次電池の制御装置およびｓｏｃ検出方法 |
| JP2014059206A (ja) * | 2012-09-18 | 2014-04-03 | Toyota Industries Corp | 充電状態推定装置及び充電状態推定方法 |
| JP6221237B2 (ja) * | 2013-01-21 | 2017-11-01 | 株式会社豊田自動織機 | 充電率推定装置および充電率推定方法 |
| JP5994680B2 (ja) * | 2013-02-27 | 2016-09-21 | 株式会社豊田自動織機 | 電池残容量推定方法及び装置 |
| US9500713B1 (en) * | 2013-03-05 | 2016-11-22 | Qnovo Inc. | Circuitry and techniques for determining an SOC of a battery/cell having a silicon-based anode and a fuel gauge therefor |
| US20160254542A1 (en) * | 2013-10-23 | 2016-09-01 | Hitachi, Ltd. | Cathode Active Material for Lithium Ion Secondary Batteries, and Lithium Ion Secondary Battery |
| DE112014006399B4 (de) * | 2014-02-25 | 2019-03-21 | Mitsubishi Electric Corporation | Ladezustandschätzvorrichtung für eine Sekundärbatterie |
| US20150276885A1 (en) * | 2014-03-25 | 2015-10-01 | Stmicroelectronics International N.V. | Method of Computing State of Charge and Battery State of Charge Monitor |
| JP6129411B2 (ja) * | 2014-04-23 | 2017-05-17 | 三菱電機株式会社 | 電池残量推定装置および電池残量推定方法 |
| JP2015230193A (ja) * | 2014-06-04 | 2015-12-21 | ソニー株式会社 | 劣化状態推定装置、充電状態推定装置、ｏｃｖカーブ算出生成装置および蓄電装置 |
| CN104297689B (zh) * | 2014-09-04 | 2017-05-24 | 中兴通讯股份有限公司 | 一种显示电池的电量的方法、装置及电子设备 |
| JP6525648B2 (ja) | 2015-03-10 | 2019-06-05 | 株式会社Ｎｔｔファシリティーズ | 電池容量推定システム、電池容量推定方法および電池容量推定プログラム |
| CN106569143B (zh) * | 2016-11-11 | 2019-11-08 | 常州普莱德新能源电池科技有限公司 | 一种在线计算电芯容量与soh的方法、系统及电动车辆 |
-
2018
- 2018-03-28 JP JP2018062290A patent/JP6406468B1/ja active Active
- 2018-03-28 US US16/498,280 patent/US11181584B2/en active Active
- 2018-03-28 EP EP18774740.7A patent/EP3605122A4/en active Pending
- 2018-03-28 CN CN201880022409.XA patent/CN110462413B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| US20200041574A1 (en) | 2020-02-06 |
| EP3605122A4 (en) | 2020-04-01 |
| CN110462413B (zh) | 2021-12-07 |
| US11181584B2 (en) | 2021-11-23 |
| EP3605122A1 (en) | 2020-02-05 |
| JP2018169398A (ja) | 2018-11-01 |
| CN110462413A (zh) | 2019-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6406468B1 (ja) | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム | |
| JP6409208B1 (ja) | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム | |
| JP6477957B2 (ja) | 推定装置、蓄電装置、推定方法、及びコンピュータプログラム | |
| WO2018181620A1 (ja) | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム | |
| JP6443577B2 (ja) | 管理装置、蓄電モジュール、管理方法、及びコンピュータプログラム | |
| JP7115345B2 (ja) | 推定装置、蓄電装置、推定方法、及びコンピュータプログラム | |
| US11237214B2 (en) | Estimation device, energy storage apparatus, estimation method, and computer program | |
| JP2022101196A (ja) | 劣化推定装置、劣化推定方法及びコンピュータプログラム | |
| JP6406469B1 (ja) | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム | |
| WO2019017183A1 (ja) | 推定装置、蓄電装置、推定方法、及びコンピュータプログラム | |
| WO2018181624A1 (ja) | 蓄電量推定装置、蓄電モジュール、蓄電量推定方法、及びコンピュータプログラム | |
| JP2021058017A (ja) | 推定装置、推定方法及びコンピュータプログラム | |
| JP7375473B2 (ja) | 蓄電量推定装置、蓄電量推定方法及びコンピュータプログラム | |
| JP6406470B1 (ja) | 管理装置、蓄電モジュール、管理方法、及びコンピュータプログラム | |
| WO2018221514A1 (ja) | 管理装置、蓄電モジュール、管理方法、及びコンピュータプログラム | |
| WO2022220214A1 (ja) | 電圧推定方法、電圧推定装置及び電圧推定プログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180725 |
|
| TRDD | Decision of grant or rejection written | ||
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20180725 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20180821 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180903 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6406468 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |