(57)【要約】  【課題】電流積算による残存容量と開放電圧に基づく残存容量との双方の利点を生かして精度高く残存容量を求めつつ、低温充電時の電圧ヒステリシスの影響による精度低下を防止する。  【解決手段】セル温度が基準値より高い場合或いはセル温度が基準値以下でも発電状態でない場合には、通常のウェイトｗを算出し（Ｓ９）、電流積算による残存容量ＳＯＣｃと開放電圧の推定による残存容量ＳＯＣｖとを重み付け合成して残存容量ＳＯＣを算出する。一方、セル温度が基準値以下且つ発電状態である場合には、通常時よりも残存容量ＳＯＣｃの重みを大きくした低温発電用のウェイトｗを算出し（Ｓ１１）、最終的な残存容量ＳＯＣを合成・算出する。これにより、残存容量ＳＯＣｃ，ＳＯＣｖの双方の利点を生かして精度高く残存容量を求めつつ、低温充電時の電圧ヒステリシスの影響による残存容量ＳＯＣの精度低下を防止することができる。【選択図】図７
2006038494.tif 000002
【特許請求の範囲】
【請求項１】
蓄電デバイスの充放電電流の積算値に基づいて第１の残存容量を算出する第１の演算手段と、
上記蓄電デバイスの内部インピーダンスから推定した開放電圧に基づいて第２の残存容量を算出する第２の演算手段と、
上記第１の残存容量と上記第２の残存容量とを上記蓄電デバイスの使用状況に応じて設定したウェイトを用いて重み付け合成した第３の残存容量を、上記蓄電デバイスの最終的な残存容量として算出する第３の演算手段と、
上記蓄電デバイスが基準温度以下で充電状態にあるとき、上記ウェイトの値を上記第１の残存容量の重みを大きくする方向に調整するウェイト調整手段とを備えたことを特徴とする蓄電デバイスの残存容量演算装置。
【請求項２】
上記ウェイト調整手段は、
上記蓄電デバイスへの充電量に応じて上記ウェイトの値を調整することを特徴とする請求項１記載の蓄電デバイスの残存容量演算装置。
【請求項３】
上記ウェイト調整手段は、
上記ウェイトを、上記蓄電デバイスの充放電電流の電流変化率に基づいて設定することを特徴とする請求項１又は２記載の蓄電デバイスの残存容量演算装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、二次電池や電気化学キャパシタ等の蓄電デバイスの残存容量を演算する蓄電デバイスの残存容量演算装置に関する。
【背景技術】
【０００２】
近年、ニッケル水素電池やリチウムイオン電池等の二次電池、電気二重層キャパシタ等の電気化学キャパシタといった蓄電デバイスの小型軽量化・高エネルギー密度化が進み、携帯型の情報通信機器から電気自動車やハイブリッド自動車等の電源として活発に利用されている。
【０００３】
このような蓄電デバイスを有効に活用するには、その残存容量を正確に把握することが重要であり、従来から、蓄電デバイスの充放電電流を積算して残存容量を求める技術や、開放電圧に基づいて残存容量を求める技術が知られている。
【０００４】
例えば、特許文献１には、電気自動車の車両停止時の電池電圧から求めた開放電圧により停止時残存容量を求めると共に、電池の放電電流の積算値に基づいて放電電気容量を検出し、この放電電気容量と停止時残存容量とから満充電容量を算出し、この満充電容量と放電電気容量とから残存容量を求める技術が開示されている。
【０００５】
また、特許文献２には、リチウムイオン電池のような電池容量と電池電圧とに直線的な比例関係があるものにおいて、任意の時間のあいだ放電または充電したときの電流積算量と、放電または充電前の電圧、放電または充電後の電圧より、残存容量を求める技術が開示されている。
【０００６】
更に、特許文献３には、電池の充放電電流を積分して求めた残存容量と、電池の開放端子電圧に基づいて推定した残存容量との差の変化率に基づいて、残存容量の演算方法を補正する技術が開示されている。
【特許文献１】特開平６−２４２１９３号公報
【特許文献２】特開平８−１７９０１８号公報
【特許文献３】特開平１１−２２３６６５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、充放電電流を積算して残存容量を求める技術と開放電圧の推定値に基づいて残存容量を求める技術とは、それぞれに一長一短があり、前者は、突入電流等の負荷変動に強く、安定した残存容量が得られる反面、誤差が累積し易い（特に、高負荷継続時には誤差が大きくなる）という欠点があり、また、後者は、通常の使用時において、正確な値を求めることができる反面、短時間で負荷が大きく変動した場合に演算値が変動しやすいという欠点がある。
【０００８】
従って、特許文献１，２，３のように、単に、両者の技術を組合わせただけでは、電流積算による誤差の累積を排除することは困難である。特に、ハイブリッド車等のように充放電が連続する状態では、残存容量の演算精度が低下したり、残存容量の演算値が急激に変化するといった事態が生じる虞があり、安定した精度を確保することは困難である。
【０００９】
しかも、特許文献１，２，３に開示の技術では、低温時に顕著となる電池の電圧ヒステリシスの影響を考慮しておらず、特に、低温の充電時には、充電に伴って電池の端子電圧が急上昇して開放電圧の推定精度が悪化し、充放電電流の積算値に基づく残存容量と、開放電圧の推定値に基づく残存容量とを一義的に組合わせただけでは、残存容量の精度低下を回避することは困難である。
【００１０】
本発明は上記事情に鑑みてなされたもので、電流積算による残存容量と開放電圧に基づく残存容量との双方の利点を生かして精度高く残存容量を求めつつ、低温充電時の電圧ヒステリシスの影響による精度低下を防止することのできる蓄電デバイスの残存容量演算装置を提供することを目的としている。
【課題を解決するための手段】
【００１１】
上記目的を達成するため、本発明による蓄電デバイスの残存容量演算装置は、蓄電デバイスの充放電電流の積算値に基づいて第１の残存容量を算出する第１の演算手段と、上記蓄電デバイスの内部インピーダンスから推定した開放電圧に基づいて第２の残存容量を算出する第２の演算手段と、上記第１の残存容量と上記第２の残存容量とを上記蓄電デバイスの使用状況に応じて設定したウェイトを用いて重み付け合成した第３の残存容量を、上記蓄電デバイスの最終的な残存容量として算出する第３の演算手段と、上記蓄電デバイスが基準温度以下で充電状態にあるとき、上記ウェイトの値を上記第１の残存容量の重みを大きくする方向に調整するウェイト調整手段とを備えたことを特徴とする。
【００１２】
その際、蓄電デバイスが基準温度以下の充電状態では、ウェイトの値を充電量に応じて調整することが望ましい。また、ウェイトは、蓄電デバイスの充放電電流の電流変化率に基づいて設定することが望ましい
【発明の効果】
【００１３】
本発明の蓄電デバイスの残存容量演算装置は、電流積算による残存容量と開放電圧に基づく残存容量との双方の利点を生かして精度高く残存容量を求めつつ、低温充電時の電圧ヒステリシスの影響による精度低下を防止することができる。
【発明を実施するための最良の形態】
【００１４】
以下、図面を参照して本発明の実施の形態を説明する。図１〜図１２は本発明の実施の一形態に係わり、図１はハイブリッド車への適用例を示すシステム構成図、図２はバッテリ残存容量の推定アルゴリズムを示すブロック図、図３は等価回路モデルを示す回路図、図４は電流の移動平均処理無しの場合の残存容量を示す説明図、図５は電流の移動平均処理有りの場合の残存容量を示す説明図、図６は実車走行時の残存容量演算結果を示す説明図、図７はバッテリ残存容量推定処理のフローチャート、図８は電流容量テーブルの説明図、図９はインピーダンステーブルの説明図、図１０は残存容量テーブルの説明図、図１１は通常用ウェイトテーブルの説明図、図１２は低温発電用ウェイトテーブルの説明図である。
【００１５】
図１は、本発明をエンジンとモータとを併用して走行するハイブリッド車両（ＨＥＶ）に適用した例を示し、同図において、符号１は、ＨＥＶの電源ユニットである。この電源ユニット１には、蓄電デバイスとして例えば複数のセルを封止した電池パックを複数個直列に接続して構成されるバッテリ２と、バッテリ２の残存容量の演算、バッテリ２の冷却や充電の制御、異常検出及び異常検出時の保護動作等のエネルギーマネージメントを行う演算ユニット（演算ＥＣＵ）３とが１つの筐体内にパッケージされている。
【００１６】
尚、本形態においては、蓄電デバイスとしてリチウムイオン二次電池を例に取って説明するが、本発明による残存容量の演算手法は、電気化学キャパシタやその他の二次電池にも適用可能である。
【００１７】
演算ＥＣＵ３は、マイクロコンピュータ等から構成され、電圧センサ４で測定したバッテリ２の端子電圧Ｖ、電流センサ５で測定したバッテリ２の充放電電流Ｉ、温度センサ６で測定したバッテリ２の温度（セル温度）Ｔに基づいて、所定時間ｔ毎に充電状態(State of charge;SOC）すなわち残存容量ＳＯＣ(t)を演算する。この残存容量ＳＯＣ(t)は、電源ユニット１の演算ＥＣＵ３から、例えばＣＡＮ（Controller Area Network）通信等を介してＨＥＶ制御用電子制御ユニット（ＨＥＶ制御用ＥＣＵ）１０に出力され、車両制御用の基本データ、バッテリ残量や警告用の表示用データ等として使用される。
【００１８】
尚、後述するように、演算ＥＣＵ３は、残存容量ＳＯＣの演算に際し、周期的な演算における１演算周期前のデータ（後述する電流積算による残存容量演算の際のベース値）ＳＯＣ(t-1)を用いており、また、バッテリ２がモータの発電によって充電状態にあることを示す発電信号の入力により、低温時に顕著となるバッテリ２の電圧ヒステリシスの影響を低減し、演算精度を確保するようにしている。
【００１９】
ＨＥＶ制御用ＥＣＵ１０は、同様にマイクロコンピュータ等から構成され、運転者からの指令に基づいて、ＨＥＶの運転、その他、必要な制御を行う。すなわち、ＨＥＶ制御用ＥＣＵ１０は、電源ユニット１からの信号や図示しないセンサ・スイッチ類からの信号により、車両の状態を検出し、バッテリ２の直流電力を交流電力に変換してモータ１５を駆動するインバータ２０を初めとして、図示しないエンジンや自動変速機等を、専用の制御ユニットを介して或いは直接的に制御する。
【００２０】
演算ＥＣＵ３における残存容量ＳＯＣの演算は、図２に示す推定アルゴリズムに従って実行される。このＳＯＣ推定アルゴリズムでは、バッテリ２で測定可能なパラメータ、すなわち、端子電圧Ｖ、電流Ｉ、温度Ｔを用い、第１〜３の演算手段としての機能により、電流積算に基づく第１の残存容量としての残存容量ＳＯＣｃと、バッテリ開放電圧の推定値に基づく第２の残存容量としての残存容量ＳＯＣｖとを並行して演算する。そして、残存容量ＳＯＣｃ，ＳＯＣｖを、バッテリ２の使用状況に応じて随時変化させるウェイト（重み係数）ｗを用いて重み付け合成した第３の残存容量としての残存容量ＳＯＣを、バッテリ２の残存容量として算出する。
【００２１】
この場合、バッテリ２には、低温時に大きな電圧ヒステリシスが発生し、特に、低温状態での充電時には、バッテリ２の端子電圧が急上昇して開放電圧の推定値を正確に求めることができず、残存容量ＳＯＣｖの精度低下が懸念される。従って、演算ＥＣＵ３には、更に、ウェイト調整手段としての機能が備えられており、バッテリ２が基準温度以下の低温時で充電状態にあるとき、電流積算に基づく残存容量ＳＯＣｃの重みを大きくする方向にウェイトｗの値を調整することで、低温時に顕著となる電圧ヒステリシスの影響を回避するようにしている。このウェイト調整手段としての機能については、後述する。
【００２２】
一般的に、バッテリの残存容量を算出する技術としては、バッテリ電流の積算値に基づて残存容量を求める技術と、バッテリの開放電圧に基づいて残存容量を求める技術とがあり、それぞれに一長一短がある。前者は、突入電流等の負荷変動に強く、安定した残存容量が得られる反面、電流誤差が累積し易い（特に、高負荷継続時には誤差が大きくなる）という欠点がある。また、後者は、電流が安定している領域では、正確な値を求めることができる反面、短時間で負荷が大きく変動した場合やバッテリの電圧ヒステリシスが大きくなる低温時には、開放電圧を推定する際のインピーダンスを正確に求めることができず、残存容量の算出値が振動し易いという欠点がある。
【００２３】
従って、本ＳＯＣ推定アルゴリズムでは、電流Ｉを積算して求めた残存容量ＳＯＣｃ(t)と、バッテリ開放電圧の推定値から求めた残存容量ＳＯＣｖ(t)とを、バッテリ２の使用状況に応じて随時変化させるウェイト（重み係数）ｗにより重み付けして合成することにより、残存容量ＳＯＣｃ，ＳＯＣｖ双方の欠点を打消して互いの利点を最大限に引き出すようにしている。ウェイトｗは、ｗ＝０〜１の間で変化させ、合成後の最終的な残存容量ＳＯＣ(t)は、以下の（１）式で与えられる。  
ＳＯＣ(t)＝ｗ・ＳＯＣｃ(t)＋(１−ｗ)・ＳＯＣｖ(t)…（１）
【００２４】
ウェイトｗは、現在のバッテリの使用状況を的確に表すことのできるパラメータを用いて決定する必要があり、そのパラメータとしては、単位時間当たりの電流の変化率や残存容量ＳＯＣｃ，ＳＯＣｖの間の偏差等を用いることが可能である。単位時間当たりの電流変化率は、バッテリの負荷変動を直接的に反映しているが、単なる電流変化率では、スパイク的に発生する電流の急激な変化の影響を受けてしまう。
【００２５】
従って、本形態においては、スパイク的な電流の変化の影響を防止するため、所定のサンプリング数の単純平均、移動平均、加重平均等の処理を施した電流変化率を用いるようにしており、特に、電流の遅れを考慮した場合、バッテリの充放電状態の変化に対して、過去の履歴を過剰となることなく適切に反映することのできる移動平均を用いてウェイトｗを決定するようにしている。
【００２６】
この電流Ｉの移動平均値に基づいてウェイトｗを決定することにより、電流Ｉの移動平均値が大きいときには、電流積算に基づく残存容量ＳＯＣｃのウェイトを高くして開放電圧の推定値に基づく残存容量ＳＯＣｖのウェイトを下げ、負荷変動の影響を電流積算によって正確に反映すると共に、開放電圧推定時の振動を防止することができる。逆に、電流Ｉの移動平均値が小さいときには、電流積算に基づく残存容量ＳＯＣｃのウェイトを下げ、開放電圧の推定値に基づく残存容量ＳＯＣｖのウェイトを高くすることにより、電流積算時の誤差の累積による影響を回避し、開放電圧の推定により正確な残存容量を算出することができる。
【００２７】
すなわち、電流Ｉの移動平均は、電流の高周波成分に対するローパスフィルタとなり、この移動平均のフィルタリングにより、走行中の負荷変動で発生する電流のスパイク成分を、遅れ成分を助長することなく除去することができる。これにより、バッテリ状態をより的確に把握することができ、残存容量ＳＯＣｃ，ＳＯＣｖ双方の欠点を打消して互いの利点を最大限に引き出し、残存容量の推定精度を大幅に向上することができる。
【００２８】
以上の電流Ｉの移動平均によるウェイトｗは、バッテリ２の電圧ヒステリシスの影響を無視できる通常の状態では、バッテリの使用状況に応じて残存容量ＳＯＣｃ，ＳＯＣｖの重み付けを適切に設定することができるが、電圧ヒステリシスが顕著となる低温状態では、ウェイトｗの値を調整する必要がある。すなわち、低温時には、バッテリ２の電圧ヒステリシスの影響が大きくなり、特に、発電時（バッテリ充電時）には、発電（充電）開始と共にバッテリ２の端子電圧が急上昇し、過渡的な電圧変動が発生する。従って、この過渡的な電圧変動によってバッテリのインピーダンスを正確に求めることが困難となり、開放電圧の推定値に基づく残存容量ＳＯＣｖの精度が悪化する虞がある。
【００２９】
このため、演算ＥＣＵ３は、ウェイト調整手段としての機能により、低温発電時（低温充電時）には、ウェイトｗの値を調整して電流積算のウェイトを高くし、開放電圧推定のウェイトを下げることにより、残存容量ＳＯＣｖの精度低下を補償する。具体的には、バッテリ温度Ｔと電流変化率ΔＩ／Δｔ（詳細には、低温になる程、バッテリの内部インピーダンスが増加して電流変化率が小さくなるため、温度補正した補正後電流変化率ｋΔＩ／Δｔを用いる）とをパラメータとしてウェイトｗを格納するテーブルを、通常時に使用する通常用ウェイトテーブルと低温発電時に使用する低温発電用ウェイトテーブルとの２つのテーブル（図１１，図１２参照）に分けて作成しておき、バッテリ温度（セル温度）及び発電状態に応じて使い分ける。
【００３０】
低温発電用ウェイトテーブルには、通常用ウェイトテーブルに比較して、テーブルの同じアドレス（参照位置）には、電流積算に基づく残存容量ＳＯＣｃの重みを大きくするようなウェイトｗの値が格納されている。低温発電時には、この低温発電用ウェイトテーブルを用いることにより、電圧ヒステリシスによる残存容量ＳＯＣｖの精度低下を補償し、電圧ヒステリシスの影響を受けにくい残存容量ＳＯＣｃの重みを大きくして残存容量ＳＯＣの精度低下を防止することができる。
【００３１】
更に、本ＳＯＣ推定アルゴリズムの特徴として、電池理論に基づいてバッテリ内部状況を電気化学的に把握し、バッテリ開放電圧に基づく残存容量ＳＯＣｖの演算精度の向上を図っている。以下、本推定アルゴリズムによる残存容量ＳＯＣｃ，ＳＯＣｖの演算について詳述する。
【００３２】
先ず、電流積算による残存容量ＳＯＣｃは、以下の（２）式に示すように、ウェイトｗを用いて合成した残存容量ＳＯＣをベース値として、所定時間毎に電流Ｉを積算して求められる。  
ＳＯＣｃ(t)＝ＳＯＣ(t-1)−∫[(１００ηＩ／Ａｈ)＋ＳD]ｄｔ／３６００…（２）
但し、η  ：電流効率
Ａｈ：電流容量（温度による変数）
ＳD ：自己放電率
【００３３】
（２）式における電流効率η及び自己放電率ＳDは、それぞれ定数と見なすことができるが（例えば、η＝１、ＳD＝０）、電流容量Ａｈは、温度に依存して変化する。従って、この電流積算による残存容量ＳＯＣｃの算出に際しては、温度によるセル容量の変動を関数化して算出した電流容量Ａｈを用いている。
【００３４】
また、（２）式による残存容量ＳＯＣｃ(t)の演算は、具体的には演算ＥＣＵ３における離散時間処理によって実行され、１演算周期前の合成残存容量ＳＯＣ(t-1)を、電流積算のベース値（初期値）として入力している（図２のブロック図における遅延演算子Ｚ-1）。従って、誤差が累積したり、発散することがなく、万一、初期値が真値と大きく異なっていても、所定の時間経過後（例えば、数分後）には、真値に収束させることができる。
【００３５】
一方、開放電圧の推定に基づく残存容量ＳＯＣｖを求めるには、先ず、図３に示す等価回路モデルを用いてバッテリの内部インピーダンスＺを求める。この等価回路は、抵抗分Ｒ１〜Ｒ３、容量分Ｃ１，ＣＰＥ１，ＣＰＥ２（但し、ＣＰＥ１，ＣＰＥ２は二重層容量分）の各パラメータを、直列及び並列に組合わせた等価回路モデルであり、交流インピーダンス法における周知のCole-Coleプロットをカーブフィッティングすることにより、各パラメータを決定する。
【００３６】
これらのパラメータから求められるインピーダンスＺは、バッテリの温度や電気化学的な反応速度、充放電電流の周波数成分によって大きく変化する。従って、インピーダンスＺを決定するパラメータとして、前述の単位時間当たりの電流Ｉの移動平均値を周波数成分の置き換えとして採用し、電流Ｉの移動平均値と温度Ｔとを条件とするインピーダンス測定を行ってデータを蓄積した後、温度Ｔと単位時間当たりの電流Ｉの移動平均値とに基づいてインピーダンスＺのテーブル（後述する図９のインピーダンステーブル）を作成する。そして、このテーブルを利用してインピーダンスＺを求め、このインピーダンスＺと、実測した端子電圧Ｖと電流Ｉとから、以下の（３）式を用いて開放電圧Ｖｏの推定値を求める。  
Ｖ＝Ｖｏ−Ｉ・Ｚ…（３）
【００３７】
開放電圧Ｖｏの推定後は、バッテリ内の電気化学的な関係に基づいて残存容量ＳＯＣｖを演算する。具体的には、平衡状態での電極電位とイオンの活量との関係を記述した周知のネルンストの式を適用し、開放電圧Ｖｏと残存容量ＳＯＣｖとの関係を表すと、以下の（４）式を得ることができる。  
Ｖｏ＝Ｅ＋[(Ｒｇ・Ｔ／Ｎｅ・Ｆ)×ｌｎＳＯＣｖ／(１００−ＳＯＣｖ)]＋Ｙ…（４）
但し、Ｅ  ：標準電極電位（本形態のリチウムイオン電池では、Ｅ＝３．７４５）
Ｒｇ：気体定数（８．３１４Ｊ／ｍｏｌ−Ｋ）
Ｔ  ：温度（絶対温度Ｋ）
Ｎｅ：イオン価数（本形態のリチウムイオン電池では、Ｎｅ＝１）
Ｆ  ：ファラデー定数（９６４８５Ｃ／ｍｏｌ）
【００３８】
尚、（４）式におけるＹは補正項であり、常温における電圧−ＳＯＣ特性をＳＯＣの関数で表現したものである。ＳＯＣｖ＝Ｘとすると、以下の（５）式に示すように、ＳＯＣの三次関数で表すことができる。  
Ｙ＝−１０-6Ｘ3＋９・１０-5Ｘ2＋０．０１３Ｘ−０．７３１１…（５）
【００３９】
以上の（４）式により、残存容量ＳＯＣｖには、開放電圧Ｖｏのみならず温度Ｔとの間にも強い相関性があることがわかる。この場合、開放電圧Ｖｏと温度Ｔとをパラメータとして、直接、（４）式を用いて残存容量ＳＯＣｖを算出することも可能であるが、実際には、使用する電池特有の充放電特性や使用条件等に対する考慮が必要となる。
【００４０】
従って、以上の（４）式の関係から実際の電池の状態を把握する場合には、常温でのＳＯＣ−Ｖｏ特性を基準として、各温度域での充放電試験或いはシミュレーションを行い、実測データを蓄積する。そして、蓄積した実測データから開放電圧Ｖｏと温度Ｔとをパラメータとする残存容量ＳＯＣｖのテーブル（後述する図１０の残存容量テーブル）を作成しておき、このテーブルを利用して残存容量ＳＯＣｖを求める。そして、前述の（１）式に示したように、電流積算に基づく残存容量ＳＯＣｃと開放電圧Ｖｏの推定値に基づく残存容量ＳＯＣｖとがウェイトｗを用いて重み付け合成され、最終的な残存容量ＳＯＣが算出される。
【００４１】
ここで、残存容量の演算における電流の移動平均処理の有無の影響を比較すると、電流の移動平均処理を行うことなく残存容量ＳＯＣｖを算出した場合には、図４に示すように、電流のスパイク成分の影響を受けて局所的な残存容量ＳＯＣｖの急激な変化が発生し、最終的な合成残存容量ＳＯＣの精度を低下させる原因となる。これに対し、電流の移動平均処理を行って残存容量ＳＯＣｖを算出した場合には、図５に示すように、残存容量ＳＯＣｖから電流のスパイク成分の影響が除去され、比較的負荷変動が小さい条件下での残存容量を正確に把握することが可能となる。
【００４２】
実走行時の残存容量の演算結果は、図６に示され、比較的アップダウンの多い走行条件でセル温度が略４５°Ｃの状態において、電流積算による残存容量ＳＯＣｃと合成後の残存容量ＳＯＣの変化が示されている。図６に示される経過時間１５００ｓｅｃ付近までのバッテリの充放電が繰返される状態においては、電流積算による残存容量ＳＯＣｃの演算結果が合成後の残存容量ＳＯＣに良好に反映されている。また、経過時間１５００ｓｅｃ以後、バッテリへの充電量が増加傾向にある状態において、電流積算による残存容量ＳＯＣｃの上昇が鈍化して誤差が拡大する傾向にあるが、開放電圧の推定による残存容量ＳＯＣｖ（図示せず）が合成後の残存容量ＳＯＣに重みを増して反映され、充電量の増加に応じて合成後の残存容量ＳＯＣが上昇し、精度良く残存容量の変化を捉えている。
【００４３】
次に、以上のＳＯＣ推定アルゴリズムに従った残存容量ＳＯＣｃ，ＳＯＣｖの演算及び合成処理について、図７のフローチャートを用いて説明する。
【００４４】
図７のフローチャートは、電源ユニット１の演算ＥＣＵ３におけるバッテリ残存容量推定の基本的な処理を示すものであり、同図においては、説明の都合上、電流積算による残存容量ＳＯＣｃの演算に続いて開放電圧Ｖｏの推定による残存容量ＳＯＣｖの演算を行うようにしているが、実際には、残存容量ＳＯＣｃ，ＳＯＣｖの演算は、並行して実行される。
【００４５】
図７のバッテリ残存容量推定処理は、所定時間毎（例えば、０．１ｓｅｃ毎）に実行され、先ず、ステップＳ１において、バッテリ２の端子電圧Ｖ、電流Ｉ、温度Ｔ、及び、前回の演算処理時に推定・合成した残存容量ＳＯＣ(t-1)のデータ入力の有無を調べる。尚、端子電圧Ｖは複数の電池パックの平均値、電流Ｉは複数の電池パックの電流の総和を取り、それぞれ、例えば０．１ｓｅｃ毎にデータを取得するものとする。また、温度Ｔは、例えば１０ｓｅｃ毎に取得するものとする。
【００４６】
その結果、ステップＳ１において新たなデータ入力がない場合には、そのまま本処理を抜け、新たなデータ入力がある場合、ステップＳ１からステップＳ２へ進んで、バッテリ電流容量を、図８に示す電流容量テーブルを参照して演算する。この電流容量テーブルは、温度Ｔをパラメータとして、所定の基準とする定格容量（例えば、１つの電池パック内の所定セル数を基準単位とした場合の定格電流容量）に対する容量比Ａｈ’を格納したものであり、常温（２５°Ｃ）における容量比Ａｈ’（＝１．００）に対し、低温になる程、電流容量が減少するため、容量比Ａｈ’の値が大きくなる。この電流容量テーブルから参照した容量比Ａｈ’を用い、計測対象毎の温度Ｔにおける電流容量Ａｈを算出する。
【００４７】
次に、ステップＳ３へ進み、電流容量テーブルから求めた電流容量Ａｈ、電流Ｉの入力値、１演算周期前の合成残存容量ＳＯＣ(t-1)を用い、前述の（２）式に従って、電流積算による残存容量ＳＯＣｃ(t)を算出する。更に、ステップＳ４において、電流Ｉを移動平均して単位時間当りの電流変化率ΔＩ／Δｔを取得する。この移動平均は、例えば、電流Ｉのサンプリングを０．１ｓｅｃ毎、電流積算の演算周期を０．５ｓｅｃ毎とした場合、５個のデータを移動平均する。
【００４８】
続くステップＳ５では、バッテリ等価回路のインピーダンスＺを、図９に示すインピーダンステーブルを参照して演算し、得られたインピーダンスＺからバッテリ２の開放電圧Ｖｏを推定する。このインピーダンステーブルは、電流変化率ΔＩ／Δｔ（単位時間当たりの電流Ｉの移動平均値）を温度補正した補正後電流変化率ｋΔＩ／Δｔと温度Ｔとをパラメータとして、等価回路のインピーダンスＺを格納したものであり、概略的には、補正後電流変化率ｋΔＩ／Δｔが同じ場合には、温度Ｔが低くなる程、インピーダンスＺが増加し、同じ温度では、補正後電流変化率ｋΔＩ／Δｔが小さくなる程、インピーダンスＺが増加する傾向を有している。
【００４９】
その後、ステップＳ６へ進み、電圧−ＳＯＣ特性の演算を行い、残存容量ＳＯＣｖを算出する。すなわち、温度Ｔと推定した開放電圧Ｖｏとをパラメータとして、図１０に示す残存容量テーブルを参照し、残存容量ＳＯＣｖを算出する。この残存容量テーブルは、前述したように、ネルンストの式に基づいてバッテリ内の電気化学的な状態を把握して作成したテーブルであり、概略的には、温度Ｔ及び開放電圧Ｖｏが低くなる程、残存容量ＳＯＣｖが小さくなり、温度Ｔ及び開放電圧Ｖｏが高くなる程、残存容量ＳＯＣｖが大きくなる傾向を有している。
【００５０】
尚、図９，１０においては、通常の条件下で使用される範囲のデータを示し、他の範囲のデータは記載を省略してある。
【００５１】
その後、ステップＳ７へ進み、セル温度が基準値（例えば、１０°Ｃ）以下か否かを判定する。この判定は、セル温度が低い状態で電圧ヒステリシスが大きくなり、開放電圧Ｖｏの推定が不安定となって誤差が拡大することを回避するためのものであり、セル温度Ｔが基準値以下の場合、更に、ステップＳ８で、バッテリを充電するための発電状態にあるか否かを調べる。
【００５２】
そして、ステップＳ７においてセル温度が基準値より高い場合、或いはステップＳ８においてセル温度が基準値以下であっても発電状態でない場合には、電圧ヒステリシスの影響を無視できる通常状態と判断してステップＳ９へ進み、図１１に示す通常用ウェイトテーブルを参照してウェイトｗを算出する。通常用ウェイトテーブルは、補正後電流変化率ｋΔＩ／Δｔをパラメータとする一次元テーブルであり、概略的には、補正後電流変化率ｋΔＩ／Δｔが小さくなる程、すなわち、バッテリ負荷変動が小さい程、ウェイトｗの値を小さくして電流積算による残存容量ＳＯＣｃの重みを小さくする特性に設定されている。
【００５３】
そして、通常用ウェイトテーブルからウェイトｗを算出した後、ステップＳ９からステップＳ１２へ進み、前述の（１）式に従って、電流積算による残存容量ＳＯＣｃと開放電圧Ｖｏの推定による残存容量ＳＯＣｖとをウェイトｗを用いて重み付けし、最終的な残存容量ＳＯＣ(t)を合成して算出することにより、１サイクルの本演算処理を終了する。
【００５４】
一方、ステップＳ７，Ｓ８において、セル温度が基準値以下、且つ発電状態である場合には、ステップＳ７，Ｓ８を経てステップＳ１０へ進み、バッテリ２における入出力可能な最大電力で示される入出力可能パワー量等に基づいて発電量を算出する。次いで、ステップＳ１１へ進み、図１２に示す低温発電用ウェイトテーブルを参照してウェイトｗを算出する。
【００５５】
低温発電用ウェイトテーブルは、通常用ウェイトテーブルと同様の一次元テーブルであるが、補正後電流変化率ｋΔＩ／Δｔを、発電量に応じた係数ＰＷ（ＰＷ≧１．０）で更に補正したパラメータ（ＰＷ・ｋΔＩ／Δｔ）を用いたテーブルとして構成され、基準温度以下の低温域において、相対的に通常用ウェイトテーブルよりも大きい値のウェイトｗが格納されている。係数ＰＷは、発電量に応じてウェイトｗの値を調整するためのものであり、発電量が大きい程、係数ＰＷの値を大きくしてテーブルから参照されるウェイトｗを１に近づける方向に調整する。
【００５６】
そして、ステップＳ１１で低温発電用ウェイトテーブルからウェイトｗを算出した後は、ステップＳ１２へ進み、前述の（１）式に従って、電流積算による残存容量ＳＯＣｃと開放電圧Ｖｏの推定による残存容量ＳＯＣｖとをウェイトｗを用いて重み付けし、最終的な残存容量ＳＯＣ(t)を合成して算出することにより、１サイクルの本演算処理を終了する。この低温発電時の最終的な残存容量ＳＯＣは、通常時よりも電流積算による残存容量ＳＯＣｃの重みが大きくされており、電圧ヒステリシスの影響による残存容量ＳＯＣの精度低下が防止される。
【００５７】
尚、本形態においては、通常用と低温発電用との２つのテーブルを選択的に用いる例について説明しているが、補正後電流変化率ｋΔＩ／Δｔに、電圧ヒステリシスの温度条件によって可変される係数Ｋ１を乗算したパラメータＫ１・ｋΔＩ／Δｔに基づいて作成した１つのウェイトテーブルを用いるようにしても良い。係数Ｋ１は、通常状態ではＫ１＝１とされて通常用ウェイトテーブルと等価になり、低温発電時には、Ｋ１＞１．０として発電量に応じてＫ１の値を可変することで、発電用ウェイトテーブルと等価となる。
【００５８】
以上のように、電流積算による残存容量ＳＯＣｃと開放電圧の推定値に基づく残存容量ＳＯＣｖとを用いて残存容量を演算する際に、バッテリの使用状況に応じて設定したウェイトｗを用いて互いの重み付けを最適化して演算精度を均一化することができる。しかも、バッテリの電圧ヒステリシスの影響が顕著となる低温発電時には、ウェイトｗの値を調整して最適化し、電流積算に基づく残存容量ＳＯＣｃの重みを大きくすることにより、電圧ヒステリシスによる残存容量ＳＯＣｖの精度低下に起因する合成残存容量ＳＯＣの精度低下を防止することができ、常時、正確にバッテリ（蓄電デバイス）の残存容量を求めることができる。
【図面の簡単な説明】
【００５９】
【図１】ハイブリッド車への適用例を示すシステム構成図
【図２】バッテリ残存容量の推定アルゴリズムを示すブロック図
【図３】等価回路モデルを示す回路図
【図４】電流の移動平均処理無しの場合の残存容量を示す説明図
【図５】電流の移動平均処理有りの場合の残存容量を示す説明図
【図６】実車走行時の残存容量演算結果を示す説明図
【図７】バッテリ残存容量推定処理のフローチャート
【図８】電流容量テーブルの説明図
【図９】インピーダンステーブルの説明図
【図１０】残存容量テーブルの説明図
【図１１】通常用ウェイトテーブルの説明図
【図１２】低温発電用ウェイトテーブルの説明図
【符号の説明】
【００６０】
１  電源ユニット
２  バッテリ
３  演算ユニット（第１，第２，第３の演算手段、ウェイト調整手段）
ＳＯＣｃ  残存容量（第１の残存容量）
ＳＯＣｖ  残存容量（第２の残存容量）
ＳＯＣ  残存容量（第３の残存容量）
Ｉ  充放電電流
Ｖｏ  開放電圧
Ｚ  インピーダンス
ｗ  ウェイト
代理人      弁理士      伊    藤      進
【図１】
2006038494.tif 000003
【図２】
2006038494.tif 000004
【図３】
2006038494.tif 000005
【図４】
2006038494.tif 000006
【図５】
2006038494.tif 000007
【図６】
2006038494.tif 000008
【図７】
2006038494.tif 000009
【図８】
2006038494.tif 000010
【図９】
2006038494.tif 000011
【図１０】
2006038494.tif 000012
【図１１】
2006038494.tif 000013
【図１２】
2006038494.tif 000014