JP7685455B2 - 状態診断装置および状態診断システム - Google Patents

Description

本発明は、状態診断装置および状態診断システムに関する。

本技術分野の背景技術として、下記特許文献１の要約には、「劣化判定装置１００の電源部１０７は、交流電源１０９から交流電力を受け付け直流電力に変換し自身に動作電力を供給する。接続部１０８は、組電池６００の測定対象電池と接続し、設定部１０２は、ユーザから、測定対象電池の放電条件情報の入力を受け付ける。放電部１０３は、放電条件情報に基づいて接続された測定対象電池の放電を制御し、放電後の測定対象電池の端子電圧を計測する。演算部１０５は、電池データ記憶部１０４内の放電特性情報と、放電条件情報と、測定対象電池の端子電圧とに基づいて残存容量を算出する。表示部１０６は、残存容量を表示する。」と記載されている。

特開２００６－３００５６１号公報

ところで、上述した技術において、一層適切に電池の状態を把握したいという要望がある。
この発明は上述した事情に鑑みてなされたものであり、バッテリの状態を適切に把握できる状態診断装置および状態診断システムを提供することを目的とする。

上記課題を解決するため本発明の状態診断装置は、計測対象であるバッテリの状態量と、前記状態量から前記バッテリの特性カーブ情報を導出する波形近似モデル情報と、に基づいて前記特性カーブ情報を取得する波形近似モデル処理部と、前記特性カーブ情報における所定条件を満たす特徴点の位置に基づいて、前記バッテリの状態を推定するバッテリ状態推定部と、複数種類の前記バッテリに対応して各種類に応じた前記特性カーブ情報の波形形状を各々決定するための複数の波形パラメータを記憶するデータベースから、計測対象の前記バッテリに対応する前記波形パラメータを選択し、前記波形近似モデル情報に適用する波形パラメータ探索部と、を備え、前記特性カーブ情報は、前記バッテリの充電率に対する端子開放電圧の特性であり、前記特徴点は、前記特性カーブ情報に生じる変曲点であり、前記バッテリ状態推定部は、前記波形パラメータの候補値である候補波形パラメータに基づいて算出した候補特性カーブ情報が、広義単調増加条件を満たすことと、前記特性カーブ情報が所定の変曲点条件を満たすことと、を条件として、前記候補波形パラメータを前記波形パラメータに設定することを特徴とする。

本発明によれば、バッテリの状態を適切に把握できる。

第１実施形態による状態診断システムのブロック図である。 ＳＯＣの変化に対するＯＣＶカーブの変化の一例を示す図である。 バッテリの時間経過に対する充電電流、ＳＯＣ、端子電圧の模式的なグラフである。 波形近似モデル情報により導出される可能性のある各種のＯＣＶカーブを示したものである。 変曲点算出部の動作説明図である。 コンピュータのブロック図である。 データベース作成ルーチンのフローチャートである。 状態推定ルーチンのフローチャート（１／２）である。 状態推定ルーチンのフローチャート（２／２）である。 第２実施形態による状態診断システムのブロック図である。 第３実施形態による状態診断システムのブロック図である。 第４実施形態による状態診断システムのブロック図である。 第５実施形態による状態診断システムのブロック図である。

［実施形態の概要］
後述する実施形態は、主に充放電を行うバッテリの状態診断システムおよび劣化診断システムに関するものであり、特に多種多様なバッテリの劣化診断を行う場合に好適である。これら実施形態の応用例としては、運用中駐車時にバッテリ交換を行うことにより航続距離を伸ばすことが可能な電動車両に付帯する車載ＢＭＳ（Battery Management System：バッテリマネジメントシステム）や、定置充電器などが挙げられる。ＢＭＳにおいて搭載バッテリのＳＯＨ（State of Health：健全性指標）の診断手段を実現するために、後述する実施形態は効果的である。

地球環境保全の有力手段として、陸海空の乗り物等の電動化、例えば電動車両が期待されている。電動車両はエネルギーをバッテリに蓄えてその電力を用いて走行するため、走行時には環境負荷物質を排出しないという特徴がある。しかし、バッテリの単位体積、単位重量当たりのエネルギー密度はガソリンや水素など化学エネルギーと比較すると低く、航続距離が内燃機関車両と比較して短いという課題がある。そのため頻繁に充電を行わざるを得ず、また充電時間も燃料の補充時間と比較すると長く、電動車両利用者の利便性を損ねている。

この課題を解決するための手段として、バッテリ交換サービスが検討されている。バッテリ交換サービスは、バッテリのＳＯＣ（State of Charge：充電率）が少なくなった車両のバッテリを、充電されたバッテリに交換する。これにより、通常なら長時間かかっていた充電時間を、交換という比較的短い時間の作業で代替することができる。交換されたバッテリは、劣化が進んでいない場合には再度満充電され、別の電動車両への交換待ち状態となる。劣化が進んでいた場合には、定置型バッテリや災害時用予備バッテリとして再利用される。転用不可能な程度まで劣化が進んだ場合にはリサイクルされる。

一般にバッテリは劣化進行による性能変化が大きいため、バッテリを有効に運用するためにはＳＯＨ（State of Health：劣化率）の正確な把握を行うことが好ましい。上述した特許文献１の内容を応用すると、バッテリのＳＯＨを推定できると考えられる。すなわち、ＳＯＣとＯＣＶ（Open Circuit Voltage：端子開放電圧）とＳＯＨとの対応関係を予め記憶手段に記憶し、そのデータを用いてＳＯＨを診断できると考えられる。

本実施形態は、特に多種多様なバッテリが搭載されることが予想される電動車両において、一層精度の高いバッテリ状態診断システムを提供するものである。
従来の電動車両にあっては、バッテリの頻繁な交換は考慮されておらず、特定のバッテリの状態診断を実現することを想定していた。バッテリの状態診断装置は車載ＢＭＳ（Battery Management System：バッテリマネジメントシステム）に搭載され、車両とバッテリの種類とは１対１に対応する。このような前提ではバッテリの材料特性や反応式など物理的な知見から導出された数式モデル（以降「物理モデル」と記載する）を詳細に構築して診断精度を向上させることができる。

しかし、バッテリ交換を前提とした場合、想定したものとは別の種類のバッテリが車両と組み合わされる場合がありうる。例えば、正極材料としてリン酸鉄を想定している車載ＢＭＳに対して正極材料がニッケル－コバルト－マンガン三元系を用いるバッテリが組み合わされることも予想される。一般に物理モデルは、精度は非常に高いものの構築に時間を要し、次々に開発される新型バッテリへの対応が困難であるという課題がある。

一方、バッテリ特性の測定結果をそのまま近似式に当てはめて状態診断のよりどころとする手法が代替案として考えられる。これらは物理式によらないためモデルの構築が容易で、新型バッテリにも迅速に対応できるという利点があるが、精度との両立は難しい。
そこで、後述する実施形態の一つの目的は、新型バッテリへの対応が迅速かつ容易で、かつ精度も向上させることができる状態推定システムを提供することである。

電動車両などのバッテリアプリケーションにおいては、バッテリのＳＯＣや寿命などの状態推定が行われる。従来は電気化学的な物理モデルを拠り所としてこれら状態推定が実現されていた。しかし、物理モデルの導出には時間を要し、新型バッテリのサービスインに時間を要する課題があった。後述する実施形態では物理モデルに拠らず、バッテリの特性カーブを直接近似する波形近似モデルを使用することで新型バッテリ導入期間の短縮を図れるという効果がある。一方、波形近似モデルは物理モデルと比較して精度的に不利であると考えられている。それについては十分な精度を確保した波形近似モデルの波形パラメータ（波形を決定するための調整因子）をデータベースに保存しておき、必要に応じて波形パラメータを参照することで、十分な精度を得ることができる。

また、波形近似モデルにより得られた近似波形をどのように判断するかといった判定基準情報も同時にデータベースから与えることもできる。これにより、充電の代わりにバッテリ交換を行うサービスや、劣化バッテリの随時交換といった、バッテリの変更が頻繁に発生する運用方法であっても柔軟に対応できるという効果がある。さらに、後述する実施形態では、前記状態推定装置は前記波形パラメータ候補を前記データベースから取得し、広義単調増加条件（ｘ1<ｘ2ならばｆ（ｘ1）≦ｆ（ｘ2））を用いて好適な波形パラメータを選択する波形パラメータ探索手段（探索部）を有する。

この手法によれば、バッテリの識別ができずデータベースを直接使えない場合に対処できる。バッテリの識別ができない場合には、データベースに格納された前記波形パラメータ候補を順次確認していくことで、最終的に好適な波形パラメータを決定する。波形パラメータの適否を判定するのが前記広義単調増加条件である。バッテリの状態推定にはSOCに対するOCVの波形（OCVカーブ）が用いられるが、OCV値をｘ，その時のSOC値をf(ｘ)とすると、広義単調増加条件（ｘ1<ｘ2ならばｆ（ｘ1）≦ｆ（ｘ2））を満たすという特徴がある。この特徴を用いることで好適な波形パラメータを判定することができ、結果としてバッテリが識別できない場合においても好適な波形パラメータを提供できるという効果がある。

［第１実施形態］
〈第１実施形態の構成〉
図１は、第１実施形態による状態診断システム１０１のブロック図である。
状態診断システム１０１は、劣化診断の対象となるバッテリ１０の状態を取得するものであり、状態診断装置１１と、バッテリ情報サーバ１４（サーバ装置）と、バッテリ運用履歴ストレージ１６と、充電器１７と、を備えている。状態診断装置１１、バッテリ情報サーバ１４およびバッテリ運用履歴ストレージ１６は、通信網１３を介して相互に各種データを送受信する。

通信網１３は、コンピュータ通信網や電話網等によって実現される。バッテリ情報サーバ１４は、各種バッテリ１０の情報を記憶するデータベース１５を備えており、該データベース１５に記憶されている情報を提供する。バッテリ運用履歴ストレージ１６は、状態診断装置１１の診断結果（詳細は後述する）を記録する。バッテリ運用履歴ストレージ１６に記録された情報は、バッテリ１０の交換通知サービス、バッテリ価格決定等、各種のバッテリ関連サービスに利用される。

状態診断装置１１は、通信部１２と、データ加工部２２と、波形近似モデル処理部２４と、バッテリ状態推定部３０と、識別情報取得部４２と、波形パラメータ探索部４４と、充電制御部７８と、を備えている。ここで、波形パラメータ探索部４４は、バッテリ状態診断装置２７を備えている。バッテリ状態推定部３０は、変曲点算出部３４と、変曲点位置比較部３６と、を備えている。

状態診断装置１１は、図示は省略するが、例えば、ＢＭＳ（Battery Management System：バッテリマネジメントシステム）充電器コントローラ、車両統合コントローラ等の装置に搭載することが好ましい。状態診断装置１１を実現するにあたっては、マイクロコントローラに書き込まれたソフトウエアやＰＬＤやＦＰＧＡなどのプログラマブルデバイス、システムオンチップなどによって実現できる。通信部１２は、通信網１３を介してバッテリ情報サーバ１４またはバッテリ運用履歴ストレージ１６と通信する。

バッテリ１０には、各種センサ（図示略）が装着されている。各センサは、バッテリ１０の端子電圧、入出力電流および温度を計測し、その結果をセンサ情報Ｄ２１（状態量）として出力する。なお、状態診断装置１１内の処理に関しては、説明の簡略化のため、温度に対応する処理に関しては、説明を省略する場合がある。状態診断装置１１内のデータ加工部２２は、必要に応じて、センサ情報Ｄ２１を劣化状態診断に用いて好適な指標や物理量であるバッテリ状態量Ｄ２３に変換する。バッテリ状態量Ｄ２３は、例えば、バッテリ１０のＳＯＣ、ＯＣＶ等の離散的な値を含む。

波形近似モデル処理部２４は、波形近似モデル情報Ｄ２８を記憶し、これによって離散的なバッテリ状態量Ｄ２３を加工し、特性カーブ情報Ｄ２５を算出する。波形近似モデル情報Ｄ２８は、ＳＯＣに対するＯＣＶの特性を規定する情報を含み、特性カーブ情報Ｄ２５は、ＳＯＣに対するＯＣＶの特性曲線の情報を含む。

波形近似モデル情報Ｄ２８はスプライン曲線や、Volterra級数によって表現することが考えられる。しかし、波形近似モデル情報Ｄ２８は、このような特定の物理式によらず、波形を近似数式で表現するモデルによって実現することが好ましい。物理モデルを用いないことにより、物理モデルを用いる場合と比較して、多種多様なバッテリの特性カーブ情報Ｄ２５を、より少ない情報量と同定手順で得ることができる。波形近似モデル情報Ｄ２８の具体例については後述する。

バッテリ状態推定部３０は、特性カーブ情報Ｄ２５に基づいて、バッテリ１０の劣化状態を推定する。バッテリ状態推定部３０に含まれる変曲点算出部３４は、通信部１２を介して、バッテリ情報サーバ１４から変曲点探索情報Ｄ３０を受信する。そして、変曲点算出部３４は、変曲点探索情報Ｄ３０を用いて、特性カーブ情報Ｄ２５における一または複数の変曲点の位置を算出し、その結果を変曲点情報Ｄ３４として出力する。

ここで、変曲点探索情報Ｄ３０とは、バッテリ１０の種別に応じて、変曲点の存在範囲、すなわちＳＯＣの範囲等を示す情報である。この変曲点探索情報Ｄ３０を用いることにより、バッテリ１０の劣化を効果的に判断できる変曲点の位置、すなわち変曲点が現れるＳＯＣ等を容易に探索することができる。

例えば、正極にリチウムを含有し、負極に黒鉛を含有するリン酸鉄系リチウムイオン電池においては、下記の事象Ａ１，Ａ２，Ａ３が現れる傾向が強い。
・Ａ１：ＳＯＣが５０％未満の領域に一つの変曲点（低ＳＯＣ側変曲点という）が現れる。
・Ａ２：ＳＯＣが５０％を超える領域に他の一つの変曲点（高ＳＯＣ側変曲点という）が現れる。
・Ａ３：変曲点の数が２個か３個になる。
また、上述した事象Ａ１，Ａ２におけるＳＯＣの領域は、正極、負極、電解液等の詳細構成を限定してゆくと、より狭い範囲にすることができる。そこで、これら事象Ａ１，Ａ２，Ａ３の条件を、変曲点探索情報Ｄ３０にするとよい。

変曲点位置比較部３６は、通信部１２を介して、バッテリ情報サーバ１４から、バッテリ状態の判定情報Ｄ３１を受信する。ここで、判定情報Ｄ３１について説明する。まず、以下の説明において、バッテリ１０の劣化度を表す指標としてはＳＯＨ（State of Health；健全性指標）を用いる。ＳＯＨは、初期の満充電容量（Ａｈ）を１００%とした際の、劣化時の満充電容量（Ａｈ）の割合として定義される。劣化が進むほどＳＯＨは減少する。また、変曲点探索情報Ｄ３０に規定された一または複数の変曲点を、「対象変曲点」と呼ぶ。

対象変曲点の位置は、バッテリ１０の劣化に伴い、初期位置から移動する。判定情報Ｄ３１は、バッテリの劣化度すなわちＳＯＨと、各対象変曲点の初期位置からの移動量と、を対応付けた情報である。すなわち、対象変曲点の移動量が判明すると、判定情報Ｄ３１に基づいて、バッテリ１０のＳＯＨを推定することができる。推定されたＳＯＨは、バッテリ状態情報Ｄ３２として、通信部１２を介して、バッテリ運用履歴ストレージ１６に記憶される。

ところで、バッテリ１０の筐体表面には、例えば、ＲＦＩＤタグやバーコード等（図示せず）が設けられており、タグ識別情報Ｄ４１（識別情報）を提示する。状態診断装置１１における識別情報取得部４２は、ＲＦＩＤタグ等（図示せず）からタグ識別情報Ｄ４１を読み取る。識別情報取得部４２は、タグ識別情報Ｄ４１を、バッテリ１０の製造ナンバや管理番号を特定するバッテリ識別情報Ｄ４３（識別情報）に変換する。バッテリ情報サーバ１４におけるデータベース１５は、このバッテリ識別情報Ｄ４３に基づいて、バッテリ１０の使用材料、仕様、特性データ等を検索し、バッテリ１０の寿命診断に適した変曲点探索情報Ｄ３０を状態診断装置１１に提供する。

しかし、ＲＦＩＤタグ等（図示せず）が損傷した場合、あるいは元々無い場合には、タグ識別情報Ｄ４１およびバッテリ識別情報Ｄ４３は得られない。このような場合の波形パラメータ探索部４４の動作については後述する。また、詳細は後述するが、バッテリ状態推定部３０は、変曲点テーブルＤ７１、ＩＤ・判定情報テーブルＤ７２およびバッテリ仕様情報テーブルＤ７４を作成する。

図２は、ＳＯＣの変化に対するＯＣＶカーブの変化の一例を示す図である。
図２の特性Ｆ２Ａ，Ｆ２Ｂにおいて、何れも横軸がＳＯＣであり、縦軸がＯＣＶである。特性Ｆ２Ａは、ＳＯＣが１００%、すなわち新品同様の状態におけるＯＣＶカーブの例である。特性Ｆ２Ｂは、劣化が進みＳＯＨが８０%になった状態におけるＯＣＶカーブの例である。特性Ｆ２Ａにおいて、ＳＯＣが５０％未満の領域には低ＳＯＣ側変曲点５１ａ（特徴点）が現れており、ＳＯＣが５０％を超える領域には高ＳＯＣ側変曲点５２ａ（特徴点）が現れている。また、特性Ｆ２Ｂには、低ＳＯＣ側変曲点５１ｂ（特徴点）と、高ＳＯＣ側変曲点５２ｂ（特徴点）と、が現れている。ここで、本実施形態において検索する変曲点とは、ＳＯＣの増加に対してＯＣＶが増加し、かつ増加率の微分値が正から負に変化している箇所のことである。

特性Ｆ２Ａ，Ｆ２Ｂの低ＳＯＣ側変曲点５１ａ，５１ｂを比較すると、劣化状態の低ＳＯＣ側変曲点５１ｂは、非劣化状態の低ＳＯＣ側変曲点５１ａより右側（ＳＯＣが大きくなる方向）に移動していることがわかる。これは高ＳＯＣ側変曲点５２ａ，５２ｂについても同様であり劣化状態の高ＳＯＣ側変曲点５２ｂは、非劣化状態の高ＳＯＣ側変曲点５２ａより右側に移動している。このように、変曲点の移動量とＳＯＨの変化とが対応するため、上述したように、変曲点の移動量からＳＯＨを診断することができる。

図３は、バッテリ１０の時間経過に対する充電電流、ＳＯＣ、端子電圧の模式的なグラフである。
図３における充電電流および端子電圧は、センサ情報Ｄ２１（図１参照）の一部である。また、データ加工部２２は、充電電流を積算することによってＳＯＣを取得する。
図２の例では説明の簡単化のため、充電と休止を繰り返す間欠充電を採用した例を示している。図１に示した充電制御部７８は、充電器１７に対して、ＳＯＣが所定値よりも低い場合にはＣＣＣ（Constant Current Charging：定電流充電）、ＳＯＣが所定値以上である場合にはＣＶＣ（Constant Voltage Charging：定電圧充電）を実行させる。充電制御部７８（図１参照）は、ＣＶＣにおいては、充電電流を供給する際のバッテリ１０の端子電圧が所定の電圧Ｖｃになるように充電器１７を制御する。

図３において充電電流が正値である期間や、充電電流を停止した直後の期間には、バッテリ１０に分極現象が生じるため、バッテリ１０の端子電圧はＯＣＶよりも高くなる。一方、充電電流を停止した後に暫く時間が経過すると、バッテリ１０は電気化学反応平衡状態になる。そのため、図示の例では、充電電流が０から正値に立ち上がる直前の端子電圧をＯＣＶであると見做して測定する。すなわち、図３に示すＳＯＣデータ５３のデータ列と、ＯＣＶデータ５４のデータ列と、がバッテリ状態量Ｄ２３（図１参照）に含まれる。

上述したように、波形近似モデル処理部２４（図３参照）は、このバッテリ状態量Ｄ２３と、波形近似モデル情報Ｄ２８と、を用いて、ＳＯＣに対するＯＣＶの特性曲線等の特性カーブ情報Ｄ２５を作成する。

なお、図３では間欠充電を適用した例を示したが、特性カーブ情報Ｄ２５は別の手法によっても取得可能である。例えば電動車両においては走行時の電流を積算してＳＯＣを算出し、車両停止時の端子電圧を測定することで同様のＯＣＶカーブが得られる。ＯＣＶカーブの作成には物理モデルの同定による手法を適用してもよい。但し、本実施形態では、上述した、波形近似モデルによる波形近似モデル情報Ｄ２８（図１参照）を用いて特性カーブ情報Ｄ２５を取得する方法を採用する。

上述したように、波形近似モデル情報Ｄ２８は、波形を近似数式で表現するモデルによって実現する。この手法は対象（バッテリ１０）の特性を物理特性によらず、データから直接導出することができる。このため、新しい特性のバッテリが次々に開発される昨今において物理モデル導出の時間がかからず、結果として新しいバッテリ１０を早期にデータベース１５に反映させることができるという利点がある。波形近似モデル情報Ｄ２８には様々な実現方法があるが、ここではＧＰＲ（Gaussian Process Regression：ガウス過程回帰）を用いる手法を説明する。ＧＰＲは、限られたデータから最も尤もらしい連続波形を生成する数学的手法である。

この手法については広く知られた手法であるため省略し結果のみ示す。この手法の詳細については、例えば「Ｃ.Ｍ.ビショップ著、パターン認識と機械学習、丸善出版(２０１２)」に記載されている。波形形状を決めるのに用いた訓練データをｘ，ｙとする。ここでｘはＳＯＣの事前測定結果、ｙはＯＣＶの事前測定結果に相当する。また、訓練データの個数をＮとする。このときｘ*の各要素に対して予測される出力値ｙ*の平均値は、［数１］によって計算できる。

［数１］において、λは解を安定化させるためのパラメータ、Ｉ_NはＮ×Ｎの単位行列である。ｋはカーネル関数であり例えば下記［数２］のように定義される。

また、［数１］におけるＫは［数２］で定義されるカーネル関数を各要素に配した行列である。

ＧＰＲを適用する利点は、少ない測定点から波形形状を再現できることにある。ＯＣＶカーブ導出の際に取得できるＯＣＶの測定点は、図３に示したように、充電電流が０であるタイミングで取得するため、取得できる条件が限られる。

従って、少ない測定点から可能な限り正確な特性カーブ情報Ｄ２５すなわちＯＣＶ波形を得るためにＧＰＲは好適である。また、ＧＰＲは、その計算過程でホワイトノイズを抑制するため、測定時のノイズやバラツキがある場合であっても安定してＯＣＶカーブが取得できるという利点がある。ここで、［数１］，［数２］において、λ，η_i，θ₀，θ₂，θ₃は本手法の使用者が自由に決定できるパラメータである。また、訓練データの選択にも自由度がある。そのため、波形近似モデル処理部２４は、採用する波形近似モデル情報Ｄ２８、によって異なるＯＣＶカーブを導出する可能性がある。異なるＯＣＶカーブが導出された例を図４に示す。

図４は波形近似モデル情報Ｄ２８により導出される可能性のある各種のＯＣＶカーブを示したものである。
図示の特性Ｆ４Ａ，Ｆ４Ｂ，Ｆ４Ｃにおいて、５５は実測点、ＯＣＶカーブ５６は真のＯＣＶカーブ、ＯＣＶカーブ５７ａ，５７ｂ，５７ｃは導出されたＯＣＶカーブである。
特性Ｆ４Ａは劣化診断に好適なＯＣＶカーブが得られた例である。また、特性Ｆ４Ｂは取得データを忠実に再現しようとしてＯＣＶカーブ５７ｂが乱高下し、特性Ｆ４Ａに示した正しい低ＳＯＣ側変曲点５１と高ＳＯＣ側変曲点５２とが埋もれてしまった例である。また、特性Ｆ４Ｃは、波形の滑らかさを重視しすぎて、特性Ｆ４Ａに示した正しい低ＳＯＣ側変曲点５１と高ＳＯＣ側変曲点５２とが消えてしまった例である。

上述した［数１］、［数２］に含まれるλ，η_i，θ₀，θ₂，θ₃等のパラメータを総称して波形パラメータＤ２６と呼ぶ。図４に示したように、波形パラメータＤ２６の設定状態によっては正しい特性カーブ情報Ｄ２５すなわちＯＣＶ波形が得られないため、本実施形態ではデータベース１５から供給された波形パラメータＤ２６を適用する。波形パラメータ探索部４４は、取得した波形パラメータＤ２６に基づいて、波形近似モデル情報Ｄ２８を決定し、これを波形近似モデル処理部２４に適用する。これにより、波形近似モデル処理部２４においては、正確なＯＣＶカーブ等の特性カーブ情報Ｄ２５を得ることが可能になる。

図５は変曲点算出部３４の動作説明図である。
上述したように、本実施形態では、［数１］および［数２］に基づいて、ＯＣＶカーブ等の特性カーブ情報Ｄ２５を取得する。従って、取得したＯＣＶカーブを二回微分した結果が０になる点から変曲点を得ることも可能である。しかし、図５に示すように、ＯＣＶカーブ５７ｄは１入力１出力のカーブであるため、より単純なグリッドサーチによって変曲点を得ることができる。その一例として、ＯＣＶカーブ５７ｄから高ＳＯＣ側変曲点５２を得る手法について説明する。

高ＳＯＣ側変曲点５２の範囲は既知であるため、その範囲内のＳＯＣ(=ｘ)を幅Δｘで微小分割する。そして分割されたＳＯＣに対してＯＣＶをそれぞれ計算することで離散ＯＣＶ値を得る。離散ＯＣＶデータのｉ番目のサンプル値をｙ_iとすると、二回微分は［数３］のように表せ、これが正から負に変化する点を以て変曲点と判断するとよい。

なお、ＳＯＣ変化に対してＯＣＶがほとんど変化しない領域を変曲点と誤認しないようにするために、ＯＣＶカーブの一回微分が一定値以上であるという条件を付けることも好ましい。また、低ＳＯＣ側変曲点（図示略）についても、高ＳＯＣ側変曲点５２と全く同様に取得できる。
図１に示した変曲点探索情報Ｄ３０は、上述のような変曲点の存在範囲、分割幅Δｘなどの各種条件を変曲点算出部３４に提供する。
以上、図２から図５まで本実施形態による状態診断装置１１についてリン酸鉄系リチウムイオン電池の特性を例として説明したが、本手法はその他のバッテリ、例えば三元系リチウムイオンバッテリに対しても全く同様に適用可能である。

図６は、コンピュータ９８０のブロック図である。
図１に示した状態診断装置１１、バッテリ情報サーバ１４およびバッテリ運用履歴ストレージ１６は、何れも図６に示すコンピュータ９８０を、１台または複数台備えている。
図６において、コンピュータ９８０は、ＣＰＵ９８１と、記憶部９８２と、通信Ｉ／Ｆ（インタフェース）９８３と、入出力Ｉ／Ｆ９８４と、メディアＩ／Ｆ９８５と、を備える。ここで、記憶部９８２は、ＲＡＭ９８２ａと、ＲＯＭ９８２ｂと、ＨＤＤ９８２ｃと、を備える。通信Ｉ／Ｆ９８３は、通信回路９８６に接続される。入出力Ｉ／Ｆ９８４は、入出力装置９８７に接続される。メディアＩ／Ｆ９８５は、記録媒体９８８からデータを読み書きする。ＲＯＭ９８２ｂには、ＣＰＵによって実行される制御プログラム、各種データ等が格納されている。ＣＰＵ９８１は、ＲＡＭ９８２ａに読み込んだアプリケーションプログラムを実行することにより、各種機能を実現する。先に図１において示した、状態診断装置１１、バッテリ情報サーバ１４の内部は、アプリケーションプログラム等によって実現される機能をブロックとして示したものである。

図７はデータベース作成ルーチンのフローチャートである。
図７において処理がステップＳ１０１に進むと、充電制御部７８（図１参照）は、充電器１７を制御することにより、ある特定のバッテリ１０に対して劣化試験を実行する。これにより、バッテリ１０について、データ加工部２２（図１参照）はセンサ情報Ｄ２１を取得し、波形近似が可能となる程度のサンプル数を有するバッテリ状態量Ｄ２３（すなわちＳＯＣ，ＯＣＶ等）を取得する。なおバッテリ１０の温度もバッテリ特性を定める重要な要素であり、センサ情報Ｄ２１は温度も含まれるが、ここでは説明の簡単化のため、温度に関する説明は省略する。

次のステップＳ１０２～１０４では、波形パラメータ探索部４４は、ステップＳ１０１で求めたデータに対して、好適なＯＣＶカーブを作成できる波形パラメータＤ２６を決定する。以下の説明において、波形パラメータＤ２６が決定されるまでの、暫定的な候補となる波形パラメータを、「候補波形パラメータＤ２６Ｃ」と呼ぶ。同様に、候補波形パラメータＤ２６Ｃに基づいて決定された波形近似モデル情報を候補波形近似モデル情報Ｄ２８Ｃと呼ぶ。同様に、候補波形近似モデル情報Ｄ２８Ｃによって得られた特性カーブ情報Ｄ２５を候補特性カーブ情報Ｄ２５Ｃと呼ぶ。なお、図１においては、図示されている符号（例えば「Ｄ２５」）に対して、アルファベット１文字を付した符号（例えば「Ｄ２５Ｃ」）については図示を省略する。

図７のステップＳ１０２では、波形パラメータ探索部４４は、候補波形パラメータＤ２６Ｃの初期値を設定する。初期値が候補波形パラメータＤ２６Ｃになる。次に、処理がステップＳ１０３に進むと、波形パラメータ探索部４４は、候補波形パラメータＤ２６Ｃに基づいて、波形近似モデル処理部２４に対して、候補波形近似モデル情報Ｄ２８Ｃを適用する。これにより、波形近似モデル処理部２４は、バッテリ状態量Ｄ２３と、候補波形近似モデル情報Ｄ２８Ｃと、に基づいて、暫定的なＯＣＶカーブを含む候補特性カーブ情報Ｄ２５Ｃを算出する。

次に、処理がステップＳ１０４に進むと、波形パラメータ探索部４４は、候補特性カーブ情報Ｄ２５ＣにおけるＯＣＶカーブが広義単調増加条件を充足しているか否かを判定する。ここで広義単調増加条件とは、ＳＯＣを変数ｘとし、ＯＣＶをｆ（ｘ）としたとき、ｘ1<ｘ2ならばｆ（ｘ1）≦ｆ（ｘ2）が成立すること、と定義する。

ステップＳ１０４において「Ｎｏ」と判定されると、処理はステップＳ１０５に進む。バッテリ１０は、「ＳＯＣが増加したにもかかわらずＯＣＶが減少することは無い」という物理的性質を有するため、ステップＳ１０４では、この性質を判定条件の一つとして利用している。仮に、この判定条件が満足されない場合は、例えば図４の特性Ｆ４Ｂに示したように、候補特性カーブ情報Ｄ２５Ｃ（ＯＣＶカーブ）が取得データに過敏に反応し、カーブが乱高下している状態であると考えられる。

そのため処理Ｓ１０５において波形パラメータ探索部４４は、候補波形パラメータＤ２６Ｃの感度を低下させる。ここで「感度を低下させる」とは、ＯＣＶカーブがより滑らかになるように候補波形パラメータＤ２６Ｃを調整することをいう。例えば、［数２］のカーネル関数を用いる場合には、η_iを小さくするとよい。η_iは、ある点の測定データが離れた点に与える影響の度合いを表しているため、これによりＯＣＶカーブを滑らかにすることができる。その後、再度ステップＳ１０３、Ｓ１０４の処理が実行され、ステップＳ１０４において「Ｙｅｓ」と判定されるまで、波形パラメータ探索部４４は波形パラメータの感度を下げ続ける。

このようにして広義単調増加条件が充足されると、ステップＳ１０４において「Ｙｅｓ」と判定され、次に処理はステップＳ１０６に進む。ここでは、バッテリ状態推定部３０の変曲点算出部３４は、候補特性カーブ情報Ｄ２５Ｃの候補変曲点情報Ｄ３４Ｃが、「変曲点条件」を充足しているか否かを判定する。ここで、候補変曲点情報Ｄ３４Ｃとは、候補特性カーブ情報Ｄ２５Ｃにおける一または複数の変曲点の位置（ＳＯＣの値）を示す情報である。また、「変曲点条件」とは、「候補変曲点情報Ｄ３４Ｃに含まれる変曲点の数が所定数以上であるか否か」である。

バッテリ１０の材料が決まれば変曲点の数と大まかな位置は所定範囲内に収まると考えられるので、この特徴を波形近似モデルの妥当性判定に用いることができる。ＯＣＶカーブの変曲点の数が多すぎる場合はＳＯＣ範囲の限定により波形近似モデルとして利用可能であるため、問題となるのは例えば図４の特性Ｆ４Ｃのように、変曲点の数が少なすぎる場合である。

候補変曲点情報Ｄ３４Ｃに含まれる変曲点の数が所定数未満であれば、ステップＳ１０６において「Ｎｏ」と判定され、処理はステップＳ１０７に進む。この場合には、変曲点算出部３４は、候補波形パラメータＤ２６Ｃの感度を上昇させる。ここで「感度を上昇させる」とは、ＯＣＶカーブがより取得データを忠実にトレースするように波形パラメータを調整することをいう。［数２］のカーネル関数を用いる場合には、例えばθ₀とθ₃を大きくするとよい。その後、候補変曲点情報Ｄ３４Ｃが変曲点条件を満たすまで、ステップＳ１０３、Ｓ１０４、Ｓ１０６、Ｓ１０７のループが繰り返される。なお、ステップＳ１０３～Ｓ１０７の繰り返し処理は、制約条件付き非線形計画法などを用いることにより自動化することができる。

ステップＳ１０４，Ｓ１０６において共に「Ｙｅｓ」と判定されると、処理はステップＳ１０８に進む。この時点で、バッテリ状態推定部３０は、現時点の候補波形パラメータＤ２６Ｃを、バッテリ１０に対する波形パラメータＤ２６として同定し、現時点の候補変曲点情報Ｄ３４Ｃに基づいて変曲点探索情報Ｄ３０を作成する。そして、バッテリ状態推定部３０は、バッテリ１０のタグ識別情報Ｄ４１と、波形パラメータＤ２６と、変曲点探索情報Ｄ３０と、を対応付けた変曲点テーブルＤ７１を作成する。

次に処理がステップＳ１０９に進むと、バッテリ状態推定部３０は、バッテリ１０のタグ識別情報Ｄ４１と、バッテリ識別情報Ｄ４３と、判定情報Ｄ３１と、を対応付けたＩＤ・判定情報テーブルＤ７２を作成する。

次に処理がステップＳ１１０に進むと、バッテリ状態推定部３０は、バッテリ仕様情報テーブルＤ７４を作成する。ここで、バッテリ仕様情報テーブルＤ７４は、正極材料、負極材料、バッテリ容量、セル接続数など、バッテリ特性に大きく影響するバッテリ１０の仕様（以下、主要仕様と呼ぶ）と、バッテリ１０のタグ識別情報Ｄ４１と、を対応付けたテーブルである。なお、セル接続数とは、多セルモジュールにおける電池セルの直列接続数、並列接続数等である。

次に、処理がステップＳ１１１に進むと、バッテリ状態推定部３０は、変曲点テーブルＤ７１、ＩＤ・判定情報テーブルＤ７２およびバッテリ仕様情報テーブルＤ７４を、それぞれ通信部１２を介してバッテリ情報サーバ１４に送信する。これにより、バッテリ情報サーバ１４は、データベース１５に対して、変曲点テーブルＤ７１、ＩＤ・判定情報テーブルＤ７２およびバッテリ仕様情報テーブルＤ７４を記録する。

ところで、本実施形態では、バッテリ識別情報Ｄ４３が得られる場合は、これに基づいてバッテリ情報サーバ１４から波形パラメータＤ２６を取得できるため、バッテリ仕様情報テーブルＤ７４は必ずしも必要ではない。しかし、バッテリ識別情報Ｄ４３が設定されていない場合や、バッテリ識別情報Ｄ４３が得られない場合も考えられるため、このような場合にバッテリ仕様情報テーブルＤ７４を用いるとよい。なお、バッテリ識別情報Ｄ４３が得られない場合の状態診断装置１１の挙動については後述する。

図８，図９は、状態推定ルーチンのフローチャートである。
図８において処理がステップＳ１２１に進むと、識別情報取得部４２はバッテリ識別情報Ｄ４３を取得できたか否かを判定する。ここで「Ｙｅｓ」と判定されると、処理はステップＳ１２２に進む。

ステップＳ１２２では、識別情報取得部４２は、バッテリ識別情報Ｄ４３に変更が生じたか否かを判定する。例えば、ユーザがバッテリ１０を交換した場合には、今回取得したバッテリ識別情報Ｄ４３は前回のものと異なる。この場合に、ステップＳ１２２では「Ｙｅｓ」と判定される。ステップＳ１２２において「Ｙｅｓ」と判定されると、処理はステップＳ１２３に進む。

ステップＳ１２３においては、識別情報取得部４２は、通信部１２を介して、バッテリ情報サーバ１４に対して、バッテリ識別情報Ｄ４３を送信する。これにより、バッテリ情報サーバ１４は、供給されたバッテリ識別情報Ｄ４３に対応する波形パラメータＤ２６と、変曲点探索情報Ｄ３０と、判定情報Ｄ３１と、を検索し、これらの情報を状態診断装置１１に送信する。次に、処理がステップＳ１２４に進むと、状態診断装置１１は、供給された波形パラメータＤ２６と、変曲点探索情報Ｄ３０と、判定情報Ｄ３１と、を取得する。

次に、処理がステップＳ１２５に進むと、バッテリ状態推定部３０および波形パラメータ探索部４４は、ステップＳ１２４で取得した情報に基づいて、状態診断装置１１の各種データ設定を行う。すなわち、波形パラメータ探索部４４は、波形パラメータＤ２６に基づいて波形近似モデル情報Ｄ２８の状態を設定する。また、バッテリ状態推定部３０は、変曲点探索情報Ｄ３０および判定情報Ｄ３１に基づいて、変曲点算出部３４および変曲点位置比較部３６の状態を設定する。

ステップＳ１２５の処理が終了すると、次に処理はステップＳ１２６に進む。また、上述したステップＳ１２２において「Ｎｏ」と判定された場合も、処理はステップＳ１２６に進む。ステップＳ１２６においては、データ加工部２２は、センサ情報Ｄ２１を取得する。次に、処理がステップＳ１２７に進むと、データ加工部２２は、センサ情報Ｄ２１に基づいて離散的なバッテリ状態量Ｄ２３を算出する。

次に、処理がステップＳ１２８に進むと、波形近似モデル処理部２４は、波形近似モデル情報Ｄ２８と、バッテリ状態量Ｄ２３と、に基づいて、特性カーブ情報Ｄ２５を算出する。次に、処理がステップＳ１２９に進むと、変曲点算出部３４は、特性カーブ情報Ｄ２５に基づいて、変曲点情報Ｄ３４を算出する。

次に、処理がステップＳ１３０に進むと、変曲点位置比較部３６は、変曲点情報Ｄ３４に基づいて、ＳＯＨを含むバッテリ状態情報Ｄ３２を算出する。次に、処理がステップＳ１３１に進むと、通信部１２は、ＳＯＨを含むバッテリ状態情報Ｄ３２を、バッテリ運用履歴ストレージ１６に送信し、本ルーチンの処理が終了する。

一方、上述したステップＳ１２１において「Ｎｏ」と判定された場合、すなわち識別情報取得部４２がバッテリ識別情報Ｄ４３を取得できなかった場合、処理は図９のステップＳ１３６に進む。図９のステップＳ１３６においては、上述したステップＳ１２６（図８参照）と同様に、データ加工部２２は、センサ情報Ｄ２１を取得する。次に、処理がステップＳ１３７に進むと、上述したステップＳ１２７（図８参照）と同様に、データ加工部２２は、センサ情報Ｄ２１に基づいて離散的なバッテリ状態量Ｄ２３を算出する。

次に、処理がステップＳ１４１に進むと、バッテリ状態推定部３０は、バッテリ１０に関する既知情報Ｄ７４Ｐをバッテリ情報サーバ１４に送信する。ここで、既知情報Ｄ７４Ｐとは、上述したバッテリ仕様情報テーブルＤ７４に含まれ得る情報のうち、既に状態診断装置１１が把握している内容を指す。

バッテリ情報サーバ１４は、通信網１３を介して既知情報Ｄ７４Ｐを受信すると、既知情報Ｄ７４Ｐに基づいて、バッテリ１０に対して推定される好適な波形近似モデル情報Ｄ２８を作成するためにデータベース１５から候補となる１件の波形パラメータを探索する。例えば、バッテリ１０の容量が既知であり、これが既知情報Ｄ７４Ｐに含まれていた場合には、バッテリ情報サーバ１４は、バッテリ１０の容量をキーとして検索することにより、候補波形パラメータＤ２６Ｃを検索する。

これにより、容量が不明の場合と比較して、波形パラメータＤ２６の候補数を削減でき、探索時間を低減することができる。また、正極材料、負極材料まで既知である場合には、候補をさらに大幅に絞り込める。バッテリ情報サーバ１４によって探索された、波形パラメータＤ２６の候補を、候補波形パラメータＤ２６Ｃと呼ぶ。そして、バッテリ情報サーバ１４は、探索した１個の候補波形パラメータＤ２６Ｃを、状態診断装置１１に送信する。

バッテリ情報サーバ１４が１個の候補波形パラメータＤ２６Ｃを送信すると、状態診断装置１１において処理はステップＳ１４２に進み、波形パラメータ探索部４４はこの候補波形パラメータＤ２６Ｃを受信する。その後、状態診断装置１１においては、ステップＳ１５４においてループを抜けるか、全ての候補波形パラメータＤ２６Ｃの受信が完了するまでステップＳ１５０～ステップＳ１５８のループを繰り返す。

まず、ステップＳ１５１において、波形パラメータ探索部４４は、候補波形パラメータＤ２６Ｃに基づいて、波形近似モデル処理部２４に対して候補波形近似モデル情報Ｄ２８Ｃを作成させる。これにより、波形近似モデル処理部２４は、バッテリ状態量Ｄ２３と、候補波形近似モデル情報Ｄ２８Ｃと、に基づいて、候補特性カーブ情報Ｄ２５Ｃを算出する。

次に、処理がステップＳ１５２に進むと、上述したステップＳ１０４の処理と同様に、波形パラメータ探索部４４は候補特性カーブ情報Ｄ２５ＣにおけるＯＣＶカーブが広義単調増加条件を充足しているか否かを判定する。

ここで「Ｙｅｓ」と判定されると、処理はステップＳ１５４に進む。ここでは、上述したステップＳ１０６の処理と同様に、バッテリ状態推定部３０の変曲点算出部３４は、候補特性カーブ情報Ｄ２５Ｃの変曲点を示す候補変曲点情報Ｄ３４Ｃが、変曲点条件を充足しているか否かを判定する。

ステップＳ１５２，Ｓ１５４の何れかにおいて「Ｎｏ」と判定されると、処理はステップＳ１５６に進む。ここでは、波形パラメータ探索部４４は、バッテリ情報サーバ１４から、次の新たな候補波形パラメータＤ２６Ｃを受信する。そして、新たな候補波形パラメータＤ２６Ｃを受信すると、波形パラメータ探索部４４において、ステップＳ１５１以降の処理が繰り返される。

一方、ステップＳ１５６において、新たな候補波形パラメータＤ２６Ｃが受信できなかった場合、すなわち、適切な候補波形パラメータＤ２６Ｃが得られないまま、データベース１５に記憶された全ての候補波形パラメータＤ２６Ｃの探索が終了した場合、処理はステップＳ１６０に進む。

ステップＳ１６０においては、バッテリ状態推定部３０は、別診断手法によってバッテリ１０のＳＯＨを推定する。例えば、ＳＯＨは使用開始時からの経過時間の平方根にほぼ比例するが知られているため、これによってＳＯＨを推定することができる。また、過去の充電回数をカウントすることによってもＳＯＨを推定することができる。ステップＳ１６０においては、このような簡易的な手法によってＳＯＨを推定するとよい。

ステップＳ１６０においてバッテリ状態推定部３０がＳＯＨを推定すると、処理はステップＳ１３１（図８参照）に進む。ここでは、上述したように、通信部１２は、ＳＯＨすなわちバッテリ状態情報Ｄ３２を、バッテリ運用履歴ストレージ１６に送信し、本ルーチンの処理が終了する。

また、候補特性カーブ情報Ｄ２５Ｃが広義単調増加条件と、変曲点条件と、を共に満たす場合、ステップＳ１５２，Ｓ１５４において共に「Ｙｅｓ」と判定され、処理はステップＳ１６４に進む。ステップＳ１６４においては、上述したステップＳ１２４と同様に、状態診断装置１１は、供給された波形パラメータＤ２６と、変曲点探索情報Ｄ３０と、判定情報Ｄ３１と、を取得（受信）する。

次に、処理がステップＳ１６５に進むと、上述したステップＳ１２５の処理と同様に、バッテリ状態推定部３０および波形パラメータ探索部４４は、ステップＳ１６４で取得した情報に基づいて、状態診断装置１１の各種データ設定を行う。すなわち、波形パラメータ探索部４４は、波形パラメータＤ２６に基づいて波形近似モデル情報Ｄ２８の状態を設定する。また、バッテリ状態推定部３０は、変曲点探索情報Ｄ３０および判定情報Ｄ３１に基づいて、変曲点算出部３４および変曲点位置比較部３６の状態を設定する。そして、上述したステップＳ１２９，Ｓ１３０，Ｓ１３１が実行され、本ルーチンの処理が終了する。

［第２実施形態］
図１０は、第２実施形態による状態診断システム１０２のブロック図である。なお、以下の説明において、上述した第１実施形態の各部に対応する部分には同一の符号を付し、その説明を省略する場合がある。
本実施形態の状態診断システム１０２においては、波形パラメータ探索部４４がバッテリ情報サーバ１４の内部に設けられている。但し、波形パラメータ探索部４４の構成要素のうち、バッテリ状態診断装置２７は状態診断装置１１に残されている。また、状態診断装置１１には、データ保持部４５が設けられている。

データ保持部４５は、データ加工部２２から出力されたバッテリ状態量Ｄ２３を保持する。さらに、データ保持部４５は、通信部１２および通信網１３を介して、バッテリ状態量Ｄ２３を、バッテリ情報サーバ１４内の波形パラメータ探索部４４に供給する。上述した以外の本実施形態の構成、動作は第１実施形態のもの（図１参照）と同様である。

状態診断装置１１は、例えば車両（図示せず）に設けられた車載用のマイクロコンピュータによって実現することが考えられる。しかし、車載用のマイクロコンピュータでは、波形パラメータ探索部４４を実現するために充分な処理能力が得られない場合も考えられる。このような場合、本実施形態のように、車両の外部に設けられたバッテリ情報サーバ１４に波形パラメータ探索部４４を設けることが好ましい。これにより、一般的に計算能力が高いバッテリ情報サーバ１４を用いて、より短い時間で最適な波形パラメータＤ２６を探索することができる。

但し、本実施形態では、状態診断装置１１とバッテリ情報サーバ１４間の通信が、バッテリ状態診断を行う場合に常に確立していることが必要となる。このため、例えば状態診断装置１１が充電器設備（図示せず）に設けられているような場合に本実施形態を採用することが好ましい。充電器は定置設備であることが多く、この場合は有線で通信路を確立できるため、無線を用いる場合と比較して通信の信頼性を高めることができる。

［第３実施形態］
図１１は、第３実施形態による状態診断システム１０３のブロック図である。なお、以下の説明において、上述した他の実施形態の各部に対応する部分には同一の符号を付し、その説明を省略する場合がある。
図１に示した構成では、車載、あるいは充電器１７の近傍等、バッテリ情報を直接取り扱うことができる箇所に状態診断装置１１を配置することが好ましかった。これに対して、図１１は、状態診断装置１１の機能をネットワーク経由でアプリケーションプログラムによって利用できるようにしたものである。

図１１において、状態診断システム１０３は、データロガー４６と、バッテリ運用履歴ストレージ１６と、バッテリ情報サーバ１４と、データロガー４６と、を備えている。また、アプリケーションサーバ装置４８は、状態診断装置１１を備えている。本実施形態における状態診断装置１１は、データ加工部２２と、波形パラメータ探索部４４と、バッテリ状態推定部３０と、通信部１２ｂと、を備えている。

また、データロガー４６は、通信部１２ａと、データ取得部４７と、を備えている。データ取得部４７は、バッテリ１０のセンサ情報を逐次またはまとめて、通信部１２ａを介してバッテリ情報サーバ１４に送信する。

アプリケーションサーバ装置４８は、アプリケーションサーバ装置側の通信部１２ｂを用いてセンサ情報Ｄ２１を受け取る。また、アプリケーションサーバ装置４８はアプリケーションサーバ装置側の通信部１２ｂを用いて、バッテリ情報サーバ１４からデータベース１５の情報を取得する。アプリケーションサーバ装置４８は、受け取ったセンサ情報Ｄ２１に基づいて、第１実施形態と同様の処理（図８、図９参照）によって、バッテリ１０のＳＯＨ等、バッテリ状態情報Ｄ３２を算出する。

算出されたバッテリ状態情報Ｄ３２は、通信部１２ｂを介して、バッテリ運用履歴ストレージ１６に送出され、バッテリ１０の交換通知サービスやバッテリ価格決定等のバッテリ関連サービスに利用される。本実施形態によれば、通信網１３上にバッテリ診断サービスを提供するアプリケーションサーバ装置４８が数多く存在する場合に好適である。

一般的に診断対象とするバッテリは非常に多いため、アプリケーションサーバ装置４８が少ない場合にはアプリケーションサーバ装置４８の負荷が高くなる。しかし、今後バッテリ診断サービスが普及し、多くのアプリケーションサーバ装置４８が提供されるとみられている。この場合にはアプリケーションサーバ装置４８の一台あたり診断対象とするバッテリの数が少なくなるため、高性能なアプリケーションサーバ装置４８の負荷が軽くなり、十分な性能が得られる。この場合、状態診断装置１１が車載機器として搭載される場合と比較してバージョンアップなどの対応が容易になるという利点がある。

［第４実施形態］
図１２は、第４実施形態による状態診断システム１０４のブロック図である。なお、以下の説明において、上述した他の実施形態の各部に対応する部分には同一の符号を付し、その説明を省略する場合がある。
図１２において、電動車両６０（車両）は、バッテリ１０と、ＢＭＳ（Battery Management System：バッテリマネジメントシステム）６１と、を備えている。そして、ＢＭＳ６１は、状態診断装置１１を備えている。

本実装形態では、電動車両６０において通信網１３を介した通信ができない場合においても、状態診断装置１１は、常にバッテリ１０の運用状態を取得することができる。バッテリ交換時に一回波形パラメータ等が更新できていれば、電波状態の良くない地域を走行している場合においても常に状態の推定を行うことができるという効果がある。また、バッテリマネジメントシステム６１はその目的のため電流、電圧センサを既に搭載しているため、状態診断装置１１を安価に実現できるという利点がある。

［第５実施形態］
図１３は、第５実施形態による状態診断システム１０５のブロック図である。なお、以下の説明において、上述した他の実施形態の各部に対応する部分には同一の符号を付し、その説明を省略する場合がある。
図１３において、電動車両６０はバッテリ１０を備えており、充電器６２は状態診断装置１１を備えている。一般に、充電器６２は定置運用されるため、重量的に余裕があり、状態診断装置１１に対して高性能なマイクロコンピュータを用いることが可能である。また、有線で通信網１３に接続することもでき、通信品質も安定している。その結果、高速にバッテリ１０の状態推定ができるという効果がある。また、推定したバッテリの劣化状態を用いて、充電方法を変更する用途にも適用可能となる効果がある。

［実施形態の効果］
以上のように上述の実施形態によれば、状態診断装置（１１）は、計測対象であるバッテリ（１０）の状態量（Ｄ２１）と、状態量（Ｄ２１）からバッテリ（１０）の特性カーブ情報（Ｄ２５）を導出する波形近似モデル情報（Ｄ２８）と、に基づいて特性カーブ情報（Ｄ２５）を取得する波形近似モデル処理部（２４）と、特性カーブ情報（Ｄ２５）における所定条件を満たす特徴点（５１ａ，５１ｂ，５２ａ，５２ｂ）の位置に基づいて、バッテリ（１０）の状態を推定するバッテリ状態推定部（３０）と、複数種類のバッテリ（１０）に対応して各種類に応じた特性カーブ情報（Ｄ２５）の波形形状を各々決定するための複数の波形パラメータ（Ｄ２６）を記憶するデータベース（１５）から、計測対象のバッテリ（１０）に対応する波形パラメータ（Ｄ２６）を選択し、波形近似モデル情報（Ｄ２８）に適用する波形パラメータ探索部（４４）と、を備える。これにより、計測対象のバッテリ（１０）に対応する波形パラメータ（Ｄ２６）を選択し、波形近似モデル情報（Ｄ２８）に適用できるため、バッテリ（１０）の状態を適切に判定できる。

また、特性カーブ情報（Ｄ２５）は、バッテリ（１０）の充電率（ＳＯＣ）に対する端子開放電圧（ＯＣＶ）の特性であり、バッテリ状態推定部（３０）は、波形パラメータ（Ｄ２６）の候補値である候補波形パラメータ（Ｄ２６Ｃ）に基づいて算出した候補特性カーブ情報（Ｄ２５Ｃ）が、広義単調増加条件を満たすことを条件として、候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）に設定すると一層好ましい。これにより、適切な候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）として設定できる。

また、特徴点（５１ａ，５１ｂ，５２ａ，５２ｂ）は、特性カーブ情報（Ｄ２５）に生じる変曲点であり、バッテリ状態推定部（３０）は、特性カーブ情報（Ｄ２５）が所定の変曲点条件を満たすことをさらなる条件として、候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）に設定すると一層好ましい。これにより、一層適切な候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）として設定できる。

また、変曲点条件は、一つの特徴点（５１ａ，５２ａ）は候補特性カーブ情報（Ｄ２５Ｃ）において充電率（ＳＯＣ）が５０％未満の領域に現れ、他の一つの特徴点（５１ｂ，５２ｂ）は候補特性カーブ情報（Ｄ２５Ｃ）において充電率（ＳＯＣ）が５０％を超える領域に現れ、かつ、全ての特徴点（５１ａ，５１ｂ，５２ａ，５２ｂ）の数は「２」または「３」である、という条件であると一層好ましい。これにより、特にバッテリ（１０）がリチウムイオンバッテリである場合に、一層適切な候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）として設定できる。

また、上述の実施形態は、他の見地によれば、バッテリ（１０）の状態を判定する状態診断システム（１０１）であって、バッテリ（１０）の状態量（Ｄ２１）に基づいて、バッテリ（１０）の状態を推定する状態診断装置（１１）と、通信網（１３）を介して、状態診断装置（１１）との間でデータを送受信するサーバ装置（１４）と、を備え、サーバ装置（１４）は、バッテリ（１０）を特定する識別情報（Ｄ４１，Ｄ４３）または状態量（Ｄ２１）と、バッテリ（１０）の特性カーブ情報（Ｄ２５）の波形形状を決定する波形パラメータ（Ｄ２６）と、を対応付けて記憶するデータベース（１５）を備え、波形近似モデル情報（Ｄ２８）と、状態量（Ｄ２１）と、に基づいてバッテリ（１０）の特性カーブ情報（Ｄ２５）を取得する波形近似モデル処理部（２４）と、特性カーブ情報（Ｄ２５）における所定条件を満たす特徴点（５１ａ，５１ｂ，５２ａ，５２ｂ）の位置に基づいてバッテリ（１０）の状態を推定するバッテリ状態推定部（３０）と、特性カーブ情報（Ｄ２５）の波形形状を決定する波形パラメータ（Ｄ２６）をデータベース（１５）から読み出す波形パラメータ探索部（４４）と、を備える。これにより、この構成においても、計測対象のバッテリ（１０）に対応する波形パラメータ（Ｄ２６）を選択し、波形近似モデル情報（Ｄ２８）に適用できるため、バッテリ（１０）の状態を適切に判定できる。

また、特性カーブ情報（Ｄ２５）は、バッテリ（１０）の充電率（ＳＯＣ）に対する端子開放電圧（ＯＣＶ）の特性であり、バッテリ状態推定部（３０）は、波形パラメータ（Ｄ２６）の候補値である候補波形パラメータ（Ｄ２６Ｃ）に基づいて算出した候補特性カーブ情報（Ｄ２５Ｃ）が、広義単調増加条件を満たすことを条件として、当該候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）に設定すると一層好ましい。これにより、適切な候補波形パラメータ（Ｄ２６Ｃ）を波形パラメータ（Ｄ２６）として設定できる。

また、波形パラメータ探索部（４４）は、サーバ装置（１４）に設けられ、通信網（１３）を介して状態量（Ｄ２１）を受信し、広義単調増加条件に基づいて波形パラメータ（Ｄ２６）を選択する機能を有すると一層好ましい。波形パラメータの探索には、通常は多数回の試行が必要となる。上述した第２実施形態（図９参照）によれば、サーバ装置（１４）が波形パラメータ探索部（４４）を備えるため、データベース（１５）に対して直接的に波形パラメータのやりとりが行えるようになる。これにより、通信によるオーバヘッドを抑制でき、処理時間を一層低減させることができる。

また、波形近似モデル処理部（２４）およびバッテリ状態推定部（３０）は、車両（６０）に搭載されたバッテリマネジメントシステム（６１）によって実現しても好ましい。特に、車両（６０）が電動車両である場合、バッテリマネジメントシステム（６１）が通信網（１３）に対する通信が確保できていない場合においても、バッテリ（１０）の運用状態を取得することができる。例えば、バッテリ（１０）の交換時に一回波形パラメータ等が更新できていれば、電波状態の良くない地域を走行している場合においても常に状態の推定を行うことができるという効果がある。

また、波形近似モデル処理部（２４）およびバッテリ状態推定部（３０）は、バッテリ（１０）を搭載した車両（６０）に接続される充電器（６２）によって実現しても好ましい。一般的に充電器（６２）は定置運用されるため、重量的に余裕があり、高性能なマイクロコンピュータを用いることが可能である。また、有線で通信網（１３）に接続することもでき、通信品質も安定している。その結果、高速にバッテリ（１０）の状態推定ができるという効果がある。また、推定したバッテリの劣化状態を用いて充電方法を変更するという用途にも適用可能となる効果がある。

また、波形近似モデル処理部（２４）およびバッテリ状態推定部（３０）は、サーバ装置（１４）によって実現しても好ましい。これにより、サーバ装置（１４）が有する有する計算機資源を有効利用することが可能となり、計算時間を一層低減することができる。

また、第３実施形態の状態診断システム１０３のように、通信網（１３）に接続されたアプリケーションサーバ装置（４８）をさらに備え、波形近似モデル処理部（２４）およびバッテリ状態推定部（３０）は、アプリケーションサーバ装置（４８）によって実現しても好ましい。特に、アプリケーションサーバ装置（４８）を複数設けることにより、バッテリ（１０）側は、アイドル状態となっているアプリケーションサーバ装置（４８）を随時選択し、バッテリ（１０）の状態推定を行わせることができる。このことにより、状態推定にかかる待ち時間を一層削減でき、平均処理時間を短縮できる。

［変形例］
本発明は上述した実施形態に限定されるものではなく、種々の変形が可能である。上述した実施形態は本発明を理解しやすく説明するために例示したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施形態の構成の一部を他の実施形態の構成に置き換えることが可能であり、ある実施形態の構成に他の実施形態の構成を加えることも可能である。また、各実施形態の構成の一部について削除し、もしくは他の構成の追加・置換をすることが可能である。また、図中に示した制御線や情報線は説明上必要と考えられるものを示しており、製品上で必要な全ての制御線や情報線を示しているとは限らない。実際には、ほとんど全ての構成が相互に接続されていると考えてもよい。上記実施形態に対して可能な変形は、例えば以下のようなものである。

（１）上記実施形態における状態診断装置１１、バッテリ情報サーバ１４等のハードウエアは一般的なコンピュータによって実現できるため、図７～図９に示したフローチャート、その他上述した各種処理を実行するプログラム等を記憶媒体に格納し、または伝送路を介して頒布してもよい。

（２）図７～図９に示した処理、その他上述した各処理は、上記実施形態ではプログラムを用いたソフトウエア的な処理として説明したが、その一部または全部をＡＳＩＣ(Application Specific Integrated Circuit；特定用途向けＩＣ)、あるいはＦＰＧＡ(Field Programmable Gate Array)等を用いたハードウエア的な処理に置き換えてもよい。

（３）上記各実施形態においては、バッテリ情報サーバ１４等は、タグ識別情報Ｄ４１またはバッテリ識別情報Ｄ４３によってバッテリ１０を特定したが、センサ情報Ｄ２１によってバッテリ１０を特定してもよい。

（４）また、上記各実施形態において、特徴点とは「変曲点」であったが、変曲点以外の点を特徴点としてもよい。

（５）状態診断装置１１の機能は、上述した各実施形態の例以外にも、様々な機器によって実現することができる。例えば、波形近似モデル処理部２４およびバッテリ状態推定部３０等を含む状態診断装置１１の機能全体をバッテリ情報サーバ１４によって実現してもよい。

１０ バッテリ
１１ 状態診断装置
１３ 通信網
１４ バッテリ情報サーバ（サーバ装置）
１５ データベース
２４ 波形近似モデル処理部
３０ バッテリ状態推定部
４４ 波形パラメータ探索部
４８ アプリケーションサーバ装置
５１ａ，５２ａ 低ＳＯＣ側変曲点（特徴点）
５１ｂ，５２ｂ 低ＳＯＣ側変曲点（特徴点）
６０ 電動車両（車両）
６１ バッテリマネジメントシステム
６２ 充電器
１０１，１０２，１０３，１０４，１０５ 状態診断システム
Ｄ２１ センサ情報（状態量）
Ｄ２５ 特性カーブ情報
Ｄ２６ 波形パラメータ
Ｄ２８ 波形近似モデル情報
Ｄ４１ タグ識別情報（識別情報）
Ｄ４３ バッテリ識別情報（識別情報）
Ｄ２５Ｃ 候補特性カーブ情報
Ｄ２６Ｃ 候補波形パラメータ

Claims (8)

1. 計測対象であるバッテリの状態量と、前記状態量から前記バッテリの特性カーブ情報を導出する波形近似モデル情報と、に基づいて前記特性カーブ情報を取得する波形近似モデル処理部と、
   前記特性カーブ情報における所定条件を満たす特徴点の位置に基づいて、前記バッテリの状態を推定するバッテリ状態推定部と、
   複数種類の前記バッテリに対応して各種類に応じた前記特性カーブ情報の波形形状を各々決定するための複数の波形パラメータを記憶するデータベースから、計測対象の前記バッテリに対応する前記波形パラメータを選択し、前記波形近似モデル情報に適用する波形パラメータ探索部と、を備え、
   前記特性カーブ情報は、前記バッテリの充電率に対する端子開放電圧の特性であり、
   前記特徴点は、前記特性カーブ情報に生じる変曲点であり、
   前記バッテリ状態推定部は、前記波形パラメータの候補値である候補波形パラメータに基づいて算出した候補特性カーブ情報が、広義単調増加条件を満たすことと、前記特性カーブ情報が所定の変曲点条件を満たすことと、を条件として、前記候補波形パラメータを前記波形パラメータに設定する
   ことを特徴とする状態診断装置。
2. 前記変曲点条件は、
   一つの前記特徴点は前記候補特性カーブ情報において前記充電率が５０％未満の領域に現れ、他の一つの前記特徴点は前記候補特性カーブ情報において前記充電率が５０％を超える領域に現れ、かつ、全ての前記特徴点の数は「２」または「３」である、という条件である
   ことを特徴とする請求項１に記載の状態診断装置。
3. バッテリの状態を判定する状態診断システムであって、
   前記バッテリの状態量に基づいて、前記バッテリの状態を推定する状態診断装置と、
   通信網を介して、前記状態診断装置との間でデータを送受信するサーバ装置と、を備え、
   前記サーバ装置は、前記バッテリを特定する識別情報または前記状態量と、前記バッテリの特性カーブ情報の波形形状を決定する波形パラメータと、を対応付けて記憶するデータベースを備え、
   波形近似モデル情報と、前記状態量と、に基づいて前記バッテリの前記特性カーブ情報を取得する波形近似モデル処理部と、
   前記特性カーブ情報における所定条件を満たす特徴点の位置に基づいて前記バッテリの状態を推定するバッテリ状態推定部と、
   前記特性カーブ情報の波形形状を決定する前記波形パラメータを前記データベースから読み出す波形パラメータ探索部と、を備え、
   前記特性カーブ情報は、前記バッテリの充電率に対する端子開放電圧の特性であり、
   前記特徴点は、前記特性カーブ情報に生じる変曲点であり、
   前記バッテリ状態推定部は、前記波形パラメータの候補値である候補波形パラメータに基づいて算出した候補特性カーブ情報が、広義単調増加条件を満たすことと、前記特性カーブ情報が所定の変曲点条件を満たすことと、を条件として、前記候補波形パラメータを前記波形パラメータに設定する
   ことを特徴とする状態診断システム。
4. 前記波形パラメータ探索部は、前記サーバ装置に設けられ、前記通信網を介して前記状態量を受信し、広義単調増加条件に基づいて波形パラメータを選択する機能を有する
   ことを特徴とする請求項３に記載の状態診断システム。
5. 前記波形近似モデル処理部および前記バッテリ状態推定部は、車両に搭載されたバッテリマネジメントシステムによって実現される
   ことを特徴とする請求項３または４に記載の状態診断システム。
6. 前記波形近似モデル処理部および前記バッテリ状態推定部は、バッテリを搭載した車両に接続される充電器によって実現される
   ことを特徴とする請求項３または４に記載の状態診断システム。
7. 前記波形近似モデル処理部および前記バッテリ状態推定部は、前記サーバ装置によって実現される
   ことを特徴とする請求項３または４に記載の状態診断システム。
8. 前記通信網に接続されたアプリケーションサーバ装置をさらに備え、
   前記波形近似モデル処理部および前記バッテリ状態推定部は、前記アプリケーションサーバ装置によって実現される
   ことを特徴とする請求項３に記載の状態診断システム。

Images