特許 7201253 劣化状態判定装置および劣化状態判定方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-12-26

(45)【発行日】2023-01-10

(54)【発明の名称】劣化状態判定装置および劣化状態判定方法

(51)【国際特許分類】

   H02J 7/00 20060101AFI20221227BHJP

   H01M 10/48 20060101ALI20221227BHJP

   H01M 10/42 20060101ALI20221227BHJP

   G01R 31/367 20190101ALI20221227BHJP

   G01R 31/389 20190101ALI20221227BHJP

   G01R 31/392 20190101ALI20221227BHJP

【ＦＩ】

H02J7/00 Q

H01M10/48 P

H01M10/42 P

H02J7/00 A

G01R31/367

G01R31/389

G01R31/392

【請求項の数】 6

(21)【出願番号】P 2020125817

(22)【出願日】2020-07-22

(62)【分割の表示】P 2019569847の分割

【原出願日】2018-10-12

(65)【公開番号】P2020191779

(43)【公開日】2020-11-26

【審査請求日】2021-08-04

(31)【優先権主張番号】P 2018061133

(32)【優先日】2018-03-28

(33)【優先権主張国・地域又は機関】JP

(73)【特許権者】

【識別番号】595098011

【氏名又は名称】東洋システム株式会社

(74)【代理人】

【識別番号】110000800

【氏名又は名称】特許業務法人創成国際特許事務所

(72)【発明者】

【氏名】宗像 一郎

(72)【発明者】

【氏名】庄司 秀樹

【審査官】佐藤 卓馬

(56)【参考文献】

【文献】特開２００５－０３７２３０（ＪＰ，Ａ）

【文献】特開２００７－０５７３７９（ＪＰ，Ａ）

【文献】特開２００９－０７１９８６（ＪＰ，Ａ）

【文献】特開２０１７－２２７５３５（ＪＰ，Ａ）

【文献】国際公開第２０１４／１８５１６３（ＷＯ，Ａ１）

【文献】特開２０１０－０６０３８４（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｈ０２Ｊ ７／００

Ｈ０１Ｍ １０／４８

Ｈ０１Ｍ １０／４２

Ｇ０１Ｒ ３１／３６７

Ｇ０１Ｒ ３１／３８９

Ｇ０１Ｒ ３１／３９２

(57)【特許請求の範囲】

【請求項1】

二次電池の一の特性を表わす第１特性パラメーターとしての前記二次電池の電流（Ｉ）の測定値（Ｉ（ｋ））に基づき、指定アルゴリズムにしたがって第１基礎指標値（ｆ（Ｉ（ｋ）））を取得し、判定対象時点（ｋ）を含む判定対象期間にわたる前記第１基礎指標値（ｆ）の値の累積値（Σ_kｆ（Ｉ（ｋ）））または積分値を前記二次電池の今回状態を表わす第１指標値（Ｆ₁）として取得する機能と、
前記二次電池の電流（Ｉ）および前記二次電池の他の特性を表わす第２特性パラメーターとしての前記二次電池の温度（Ｔ）のそれぞれの前記判定対象時点（ｋ）における測定値としての今回測定値（Ｉ（ｋ），Ｔ（ｋ））に基づき、前記二次電池と同一規格の参照二次電池の初期特性を表わす初期特性モデルにしたがって、前記初期特性モデルの出力として電流（Ｉ）の前記判定対象時点（ｋ）よりも過去の初期時点（０）における値を前記第１特性パラメーター（Ｉ）の初期特性推定値（Ｉ（０←ｋ））として取得し、 前記第１特性パラメーター（Ｉ）の初期特性推定値（Ｉ（０←ｋ））に基づき、前記指定アルゴリズムにしたがって第２基礎指標値（ｆ（Ｉ（０←ｋ）））を計算し、前記第２基礎指標値（ｆ（Ｉ（０←ｋ）））の前記判定対象期間にわたる累積値（Σ_kｆ（Ｉ（０←ｋ）））または積分値を前記二次電池の初期状態を表わす第２指標値（Ｆ₂）として取得する機能と、
前記第１指標値（Ｆ₁）および前記第２指標値（Ｆ₂）の偏差が大きいほど、前記二次電池の初期状態を基準とした劣化度が高くなるように当該劣化度（Ｄ）を導出する機能と、を備えていることを特徴とする劣化状態判定装置。

【請求項2】

請求項１に記載の劣化状態判定装置において、
前記劣化度（Ｄ）に応じた劣化診断情報Ｉｎｆ（Ｄ（ｉ））を生成する機能と、
対象機器を構成する入力インターフェースに前記劣化診断情報Ｉｎｆ（Ｄ（ｉ））を送信または入力することにより、前記対象機器を構成する出力インターフェースに前記劣化度（Ｄ）に関する情報を出力させる機能と、をさらに備えていることを特徴とする劣化状態判定装置。

【請求項3】

請求項２に記載の劣化状態判定装置において、
前記二次電池が搭載されている前記対象機器を構成する演算処理装置により特性パラメーターセンサーの出力信号に基づいて測定された前記第１特性パラメーターおよび前記第２特性パラメーターの測定値を前記出力インターフェースから受信するまたは受け付ける機能をさらに備えていることを特徴とする劣化状態判定装置。

【請求項4】

請求項２または３に記載の劣化状態判定装置と、前記対象機器と、により構成されている劣化状態判定システム。

【請求項5】

前記二次電池が搭載されている対象機器に対して劣化判定用ソフトウェアを送信することにより、前記対象機器を構成する演算処理装置に対して請求項１～３のうちいずれか１項に記載の劣化状態判定装置としての機能を付与することを特徴とするソフトウェアサーバー。

【請求項6】

二次電池の一の特性を表わす第１特性パラメーターとしての前記二次電池の電流（Ｉ）の測定値（Ｉ（ｋ））に基づき、指定アルゴリズムにしたがって第１基礎指標値（ｆ（Ｉ（ｋ）））を取得し、判定対象時点（ｋ）を含む判定対象期間にわたる前記第１基礎指標値（ｆ）の値の累積値（Σ_kｆ（Ｉ（ｋ）））または積分値を前記二次電池の今回状態を表わす第１指標値（Ｆ₁）として取得する工程と、
前記二次電池の電流（Ｉ）および前記二次電池の他の特性を表わす第２特性パラメーターとしての前記二次電池の温度（Ｔ）のそれぞれの前記判定対象時点（ｋ）における測定値としての今回測定値（Ｉ（ｋ），Ｔ（ｋ））に基づき、前記二次電池と同一規格の参照二次電池の初期特性を表わす初期特性モデルにしたがって、前記初期特性モデルの出力として電流（Ｉ）の前記判定対象時点（ｋ）よりも過去の初期時点（０）における値を前記第１特性パラメーター（Ｉ）の初期特性推定値（Ｉ（０←ｋ））として取得し、 前記第１特性パラメーター（Ｉ）の初期特性推定値（Ｉ（０←ｋ））に基づき、前記指定アルゴリズムにしたがって第２基礎指標値（ｆ（Ｉ（０←ｋ）））を計算し、前記第２基礎指標値（ｆ（Ｉ（０←ｋ）））の前記判定対象期間にわたる累積値（Σ_kｆ（Ｉ（０←ｋ）））または積分値を前記二次電池の初期状態を表わす第２指標値（Ｆ₂）として取得する工程と、
前記第１指標値（Ｆ₁）および前記第２指標値（Ｆ₂）の偏差が大きいほど、前記二次電池の初期状態を基準とした劣化度が高くなるように当該劣化度（Ｄ）を導出する工程と、を含んでいることを特徴とする劣化状態判定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、リチウムイオンバッテリー等の二次電池の劣化状態を判定するシステム等に関する。

【背景技術】

【0002】

二次電池の電流値、電圧値および周囲温度の今回測定結果に基づき、二次電池の電気的等価回路モデルにおけるモデルパラメーターの値を推定することにより、当該二次電池の劣化を診断する技術的手法が提案されている（特許文献１参照）。

【先行技術文献】

【特許文献】

【0003】

【文献】特開２０１７－０１６９９１号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかし、当該先行技術によれば、モデルパラメーターの推定値の時系列的な変化量に基づき、二次電池の劣化が診断されるため、二次電池が初期状態からどの程度劣化しているかを判定するためには、当該二次電池それ自体の初期電流値等の測定が前提とされている。よって、対象となる二次電池それ自体の電流値等の初期測定結果が存在しない場合、当該二次電池の劣化状態の判定が困難になる。

【0005】

そこで、本発明は、対象となる二次電池それ自体の特性パラメーターの初期測定結果が存在しない場合でも、当該二次電池の劣化状態を判定しうる装置等を提供することを目的とする。

【課題を解決するための手段】

【0006】

本発明に係る劣化状態判定装置は、
二次電池の一の特性を表わす第１特性パラメーターとしての前記二次電池の電流（Ｉ）の測定値（Ｉ（ｋ））に基づき、指定アルゴリズムにしたがって第１基礎指標値（ｆ（Ｉ（ｋ）））を取得し、判定対象時点（ｋ）を含む判定対象期間にわたる前記第１基礎指標値（ｆ）の値の累積値（Σ _k ｆ（Ｉ（ｋ）））または積分値を前記二次電池の今回状態を表わす第１指標値（Ｆ ₁ ）として取得する機能と、
前記二次電池の電流（Ｉ）および前記二次電池の他の特性を表わす第２特性パラメーターとしての前記二次電池の温度（Ｔ）のそれぞれの前記判定対象時点（ｋ）における測定値としての今回測定値（Ｉ（ｋ），Ｔ（ｋ））に基づき、前記二次電池と同一規格の参照二次電池の初期特性を表わす初期特性モデルにしたがって、前記初期特性モデルの出力として電流（Ｉ）の前記判定対象時点（ｋ）よりも過去の初期時点（０）における値を前記第１特性パラメーター（Ｉ）の初期特性推定値（Ｉ（０←ｋ））として取得し、 前記第１特性パラメーター（Ｉ）の初期特性推定値（Ｉ（０←ｋ））に基づき、前記指定アルゴリズムにしたがって第２基礎指標値（ｆ（Ｉ（０←ｋ）））を計算し、前記第２基礎指標値（ｆ（Ｉ（０←ｋ）））の前記判定対象期間にわたる累積値（Σ _k ｆ（Ｉ（０←ｋ）））または積分値を前記二次電池の初期状態を表わす第２指標値（Ｆ ₂ ）として取得する機能と、
前記第１指標値（Ｆ ₁ ）および前記第２指標値（Ｆ ₂ ）の偏差が大きいほど、前記二次電池の初期状態を基準とした劣化度が高くなるように当該劣化度（Ｄ）を導出する機能と、を備えている。

【発明の効果】

【0007】

本発明に係る劣化状態判定装置によれば、対象の第１特性パラメーターおよび第２特性パラメーターのそれぞれの「今回測定値（判定対象時点における測定値）」に基づき、当該対象の特性を表わす初期特性モデルにしたがって第１特性パラメーターの「初期特性推定値（初期状態としての推定値）」が計算される。さらに、第１特性パラメーターの測定値の今回時系列（判定対象時点を含む判定対象期間における測定値の時系列）に基づいて「第１指標値」が計算され、第１特性パラメーターの初期特性推定値の今回時系列（判定対象時点を含む判定対象期間における初期特性推定値の時系列）に基づいて「第２指標値」が計算される。そして、「第１指標値」および「第２指標値」に基づき、対象の劣化状態が判定される。

【0008】

このため、劣化状態の判定対象そのものの特性パラメーターの初期または過去の測定結果が存在しない場合でも、当該対象の劣化状態が判定されうる。

【図面の簡単な説明】

【0009】

【図1】本発明の一実施形態としての劣化状態判定装置の構成に関する説明図。

【図2】劣化状態判定装置の機能を表わすブロックダイヤグラム。

【図3】バッテリーの電気的等価回路モデルの一例に関する説明図。

【図4】本発明の一実施形態としての劣化状態判定方法に関する説明図。

【図5】第１指標値および第２指標値の計算のための積分区間に関する説明図。

【図6】劣化診断情報の出力例に関する説明図。

【図7】本発明に係る劣化状態判定方法のタイムチャート。

【図8】従来技術に係る劣化状態判定方法のタイムチャート。

【発明を実施するための形態】

【0010】

（劣化状態判定装置の構成）
図１に示されている本発明の一実施形態としての劣化状態判定装置は、対象機器２００とネットワークを介して通信可能な一または複数のサーバーにより構成されている。劣化状態判定装置１００は、対象機器２００に電源として搭載されている二次電池２２０の劣化状態を判定する。

【0011】

劣化状態判定装置１００は、入力要素１０２と、出力要素１０４と、第１演算処理要素１１０と、第２演算処理要素１２０と、第３演算処理要素１３０と、を備えている。

【0012】

入力要素１０２は、二次電池２２０が搭載されている対象機器２００から、特性パラメーターの測定値を受信する。出力要素１０４は、二次電池２２０の劣化状態の判定結果またはこれに応じて生成された劣化診断情報を対象機器２２０に対して送信する。

【0013】

第１演算処理要素１１０、第２演算処理要素１２０および第３演算処理要素１３０のそれぞれは、プロセッサー（演算処理装置）、メモリー（記憶装置）およびＩ／Ｏ回路等により構成されている。当該メモリーまたはこれとは別個の記憶装置には、二次電池２２０の特性を表わす特性パラメーターなどの様々なデータのほか、プログラム（ソフトウェア）が記憶保持されている。二次電池２２０またはこれが搭載されている対象機器２００の種類を識別するための複数の識別子のそれぞれと、複数のモデルのそれぞれとが対応付けられてメモリーに記憶保持されている。プロセッサーがメモリーから必要なプログラムおよびデータを読み取り、当該データに基づき、当該プログラムにしたがった演算処理を実行することにより、各演算処理要素１１０、１２０、１３０に割り当てられたタスクが実行される。

【0014】

第１演算処理要素１１０は、二次電池２２０の特性を表わす複数の特性パラメーターｐ＝（ｐ₁，‥ｐ_n）のそれぞれの今回測定値ｐ（ｋ）を初期特性モデルに入力することにより、複数の特性パラメーターｐのうち第１特性パラメーターｐ₁の初期特性推定値ｐ₁（０←ｋ）を初期特性モデルの出力として計算する。「初期特性モデル」は、規格または種類の別に応じた二次電池の初期特性を表わす。初期特性モデルは、関数Ｇにより関係式（１１０）にしたがって表わされる。

【0015】

ｐ₁（０←ｋ）＝Ｇ（ｐ（ｋ）） ‥（１１０）。

【0016】

関数Ｇは、複数の特性パラメーターｐの今回測定値ｐ（ｋ）に応じた第１初期特性パラメーターｑ₁（ｐ（ｋ））の従変数としての第２初期特性パラメーターｑ₂（ｑ₁（ｐ（ｋ）））と、第１特性パラメーターｐ₁とは別の第２特性パラメーターｐ₂の今回測定値ｐ₂（ｋ）と、を主変数とする多変数関数Ｇ₁として関係式（１１１）により定義されていてもよい。

【0017】

Ｇ＝Ｇ₁（ｑ₂（ｑ₁（ｐ（ｋ））），ｐ₂（ｋ）） ‥（１１１）。

【0018】

関数Ｇは、第２初期特性パラメーターｑ₂（ｑ₁（ｐ（ｋ）））と、第２特性パラメーターｐ₂の今回測定値ｐ₂（ｋ）と、に加えて、第２特性パラメーターｐ₂の測定値ｐ₂（ｊ）が０になった時点ｔ＝ｊにおける第１特性パラメーターｐ₁の測定値ｐ₁（ｊ）を主変数とする多変数関数Ｇ₂として関係式（１１２）により定義されていてもよい。

【0019】

Ｇ＝Ｇ₂（ｑ₂（ｑ₁（ｐ（ｋ））），ｐ₁（ｊ），ｐ₂（ｋ）） ‥（１１２）。

【0020】

第２演算処理要素１２０は、二次電池２２０の複数の特性パラメーターｐの測定値ｐ（ｋ）の時系列Ｐ（ｉ）＝｛ｐ（ｉ）｜ｉ＝ｋ，ｋ＋１，‥｝に基づき、多変数関数ｆ（ｐ）にしたがって算出される値ｆ（ｐ（ｋ））＝ｆ（ｐ₁（ｋ），ｐ₂（ｋ），‥）の累積値または時間積分値を第１指標値Ｆ₁として計算する（関係式（１２１）参照）。

【0021】

Ｆ₁＝Σ_kｆ（ｐ（ｋ）） ‥（１２１）。

【0022】

第２演算処理要素１２０は、第１特性パラメーターｐ₁の初期特性推定値ｐ₁（０←ｋ）およびその他の特性パラメーターｐ_u（ｕ＝２，３‥）の測定値ｐ_u（ｋ）に基づき、同一の多変数関数ｆ（ｐ）にしたがって算出される値ｆ（ｐ₁（０←ｋ），ｐ₂（ｋ），‥）の累積値または時間積分値を第２指標値Ｆ₂として計算する（関係式（１２２）参照）。

【0023】

Ｆ₂＝Σ_kｆ（Ｖ（０←ｋ），Ｉ（ｋ），Ｔ（ｋ）） ‥（１２２）。

【0024】

第３演算処理要素１３０は、第１指標値Ｆ₁（ｉ）および第２指標値Ｆ₂（ｉ）に基づき、関係式（１３０）にしたがって二次電池２２０の劣化度Ｄ（ｉ）を計算する。

【0025】

Ｄ（ｉ）＝｛Ｆ₂（ｉ）－Ｆ₁（ｉ）｝／Ｆ₂（ｉ） ‥（１３０）。

【0026】

対象機器２００は、入力インターフェース２０２と、出力インターフェース２０４と、制御装置２１０と、二次電池２２０と、センサー群２２２と、を備えている。対象機器２００は、パソコン、携帯電話（スマートフォン）、家電製品または電動自転車等の移動体など、二次電池２２０を電源とするあらゆる機器を含んでいる。

【0027】

制御装置２１０は、プロセッサー（演算処理装置）、メモリー（記憶装置）およびＩ／Ｏ回路等により構成されている。当該メモリーまたはこれとは別個の記憶装置には、特性パラメーターの測定値の時系列などの様々なデータが記憶保持される。制御装置２１０は、二次電池２２０から供給電力に応じて作動し、通電状態において対象機器２００の動作を制御する。対象機器２００の動作には、当該対象機器２００を構成するアクチュエータ（電動式アクチュエータなど）の動作が含まれる。制御装置２１０を構成するプロセッサーがメモリーから必要なプログラムおよびデータを読み取り、当該データに基づき、当該プログラムにしたがった演算処理を実行することにより、割り当てられたタスクを実行する。

【0028】

二次電池２２０は、例えばリチウムイオンバッテリーであり、ニッケル・カドミウム電池等のその他の二次電池であってもよい。センサー群２２２は、二次電池２２０の特性パラメーターのほか、対象機器２００の制御に必要なパラメーターの値を測定する。センサー群２２２は、例えば二次電池２２０の電圧、電流および温度のそれぞれに応じた信号を出力する電圧センサー、電流センサーおよび温度センサーにより構成されている。

【0029】

劣化状態判定装置１００は対象機器２００に搭載されていてもよい。この場合、ソフトウェアサーバー（図示略）が、対象機器２００が備えている制御装置２１０を構成する演算処理装置に対して劣化判定用ソフトウェアを送信することにより、当該演算処理装置に対して劣化状態判定装置１００としての機能を付与してもよい。

【0030】

（各演算処理要素の構成）
図２には、各演算処理要素１１０、１２０、１３０の機能を表わすブロックダイヤグラムが示されている。

【0031】

第１演算処理要素１１０は、第１初期特性モデルパラメーター記憶部１１２、第２初期特性モデルパラメーター初期特性推定値計算部１１４、および第１特性パラメーター初期特性推定値計算部１１６としての機能を備えている。第２演算処理要素１２０は、第１指標値計算部１２１および第２指標値計算部１２２としての機能を備えている。第３演算処理要素１３０は、劣化度計算部としての機能を備えている。

【0032】

第１初期特性モデルパラメーター記憶部１１２は、参照二次電池と同一規格または同一種類の任意の二次電池の初期特性を表わす第１初期特性モデルパラメーターｑ₁（ｐ）を記憶保持する。第１初期特性モデルパラメーターｑ₁（ｐ）は、二次電池の規格または種類を識別するための複数の識別子ＩＤおよび複数の特性パラメーターｐのさまざまな測定値ｐに対応する複数の値を有する。

【0033】

本実施形態では、二次電池２２０の端子間電圧Ｖ、電流Ｉおよび温度Ｔが、特性パラメーターｐ＝（ｐ₁，ｐ₂，ｐ₃）として測定される。特性パラメーターｐ＝（ｐ₁，ｐ₂，ｐ₃）は、（Ｖ，Ｔ，Ｉ）、（Ｉ，Ｖ，Ｔ）、（Ｉ，Ｔ，Ｖ）、（Ｔ，Ｉ，Ｖ）または（Ｔ，Ｖ，Ｉ）であってもよい。特性パラメーターｐは２つ（ｐ₁，ｐ₂）であってもよく、４つ以上（ｐ₁，‥ｐ_N）（４≦Ｎ）であってもよい。

【0034】

例えば、図３に示されている電気的等価回路モデルに、第１初期特性モデルパラメーターｑ₁が適用されることにより初期特性モデルが定義される。この電気的等価回路は、起電力Ｖ₀の内部電源および抵抗値ｒの内部抵抗により構成されている。二次電池２２０の電気的等価回路モデルの電気特性は、二次電池２２０の端子間電圧Ｖ、電流Ｉ、内部電源の起電力Ｖ₀および内部抵抗の抵抗値ｒに基づき、関係式（２１０）にしたがって定義される。

【0035】

Ｖ＝Ｖ₀－Ｉ・ｒ ‥（２１０）。

【0036】

規格が相違する複数の参照二次電池の電圧Ｖ（第１特性パラメーターｐ₁）、電流Ｉ（第２特性パラメーターｐ₂）および温度Ｔ（第３特性パラメーターｐ₃）と、初期状態における当該二次電池の内部抵抗の抵抗値ｒとの関係は、基準電圧Ｖ_a、基準電流Ｉ_aおよび基準温度Ｔ_aに基づき、関係式（２１２）により近似的に表わされる。

【0037】

ｒ（Ｖ，Ｉ，Ｔ）＝ｒ（Ｖ_a，Ｉ_a，Ｔ_a）＋（∂ｒ／∂Ｖ）（Ｖ－Ｖ_a）
＋（∂ｒ／∂Ｉ）（Ｉ－Ｉ_a）
＋（∂ｒ／∂Ｔ）（Ｔ－Ｔ₀） ‥（２１２）。

【0038】

基準電流Ｉ_a＝０である場合、基準電圧Ｖ_aは内部電源の起電力Ｖ₀に相当するので、関係式（２１２）は関係式（２１４）により表わされる。

【0039】

ｒ（Ｖ，Ｉ，Ｔ）＝ｒ（Ｖ₀，０，Ｔ_a）＋（∂ｒ／∂Ｖ）（Ｖ－Ｖ₀）
＋（∂ｒ／∂Ｉ）Ｉ
＋（∂ｒ／∂Ｔ）（Ｔ－Ｔ_a） ‥（２１４）。

【0040】

初期状態における二次電池の電圧Ｖ、電流Ｉおよび温度Ｔが測定され、当該測定結果に基づき、さまざまなｐ＝（Ｖ，Ｉ，Ｔ）の組み合わせごとに、偏微分係数（∂ｒ／∂ｐ_s）（ｓ＝１，２，３）が算定される。このようにして任意のｐに対して算定された偏微分係数（∂ｒ／∂ｐ_s）が、第１初期特性モデルパラメーターｑ₁（ｐ）として第１初期特性モデルパラメーター記憶部１１２に記憶保持される。離散的な偏微分係数（∂ｒ／∂ｐ_s）の算定結果に基づき、ｐ＝（Ｖ，Ｉ，Ｔ）に対する（∂ｒ／∂ｐ_s）が連続関数により近似表現され、当該連続関数が第１初期特性モデルパラメーターｑ₁（ｐ）として第１初期特性モデルパラメーター記憶部１１２に記憶保持されてもよい。

【0041】

内部抵抗の過渡応答特性が勘案された形で電気的等価回路モデルが構築されていてもよい。例えば、図３右上部分に示されているように、内部抵抗が、抵抗（抵抗値Ｒ_s）および３つの複合抵抗（インピーダンスＺ₁、Ｚ₂、Ｚ₃）が直列接続されることにより構成されていてもよい。第１複合抵抗は、抵抗（抵抗値Ｒ₀）およびインダクタ（インダクタンスＬ₀）の並列回路により構成されている。第２複合抵抗は、抵抗（抵抗値Ｒ₁）およびキャパシタ（キャパシタンスＣ₁）の並列回路により構成されている。第２複合抵抗は、一対の抵抗（抵抗値Ｒ₁）および抵抗（インピーダンスＺ_W）の直列回路とキャパシタ（キャパシタンスＣ₂）との並列回路により構成されている。内部抵抗の抵抗値ｒは、関係式（２１６）にしたがって定義される。

【0042】

ｒ＝Ｒ_s＋Ｚ₁（Ｌ₀，Ｒ₀）＋Ｚ₂（Ｃ₁，Ｒ₁）＋Ｚ₃（Ｃ₂，Ｒ₂＋Ｚ_W） ‥（２１６）。

【0043】

この場合、Ｒ_s、Ｚ₁、Ｚ₂およびＺ₃のそれぞれが、式（１２）と同様の形で近似されることにより、任意のｐ_{_m}＝（Ｖ_{_m}，Ｉ_{_m}，Ｔ_{_m}）に対して算定された偏微分係数（∂Ｚ／∂ｐ_i）（Ｚ＝Ｒ_s、Ｚ₁、Ｚ₂、Ｚ₃）が、第１初期特性モデルパラメーターｑ₁（Ｖ，Ｉ，Ｔ）として第１初期特性モデルパラメーター記憶部１１２に記憶保持されてもよい。

【0044】

第２初期特性モデルパラメーター初期特性推定値計算部１１４は、二次電池２２０の識別子ＩＤに加えて複数の特性パラメーターＶ，ＩおよびＴのそれぞれの今回測定値（Ｖ（ｋ），Ｉ（ｋ），Ｔ（ｋ））に対応する第１初期特性モデルパラメーターｑ₁を第１初期特性モデルパラメーター記憶部１１２から読み取る。第２初期特性モデルパラメーター初期特性推定値計算部１１４は、二次電池２２０の複数の特性パラメーターの今回測定値ｐ（ｋ）＝（Ｖ（ｋ），Ｉ（ｋ），Ｔ（ｋ））と、当該今回測定値ｐ（ｋ）に応じた第１初期特性モデルパラメーターｑ₁（ｐ（ｋ））と、に基づき、関係式（２１８）にしたがって、内部抵抗の抵抗値ｒを第２初期特性モデルパラメーターｑ₂の今回初期特性推定値ｑ₂（０←ｋ）として計算する。

【0045】

ｒ（０←ｋ）＝ｒ（Ｖ₀，０，Ｔ₀）＋（∂ｒ／∂Ｖ）（Ｖ（ｋ）－Ｖ₀）
＋（∂ｒ／∂Ｉ）Ｉ（ｋ）
＋（∂ｒ／∂Ｔ）（Ｔ（ｋ）－Ｔ₀） ‥（２１８）。

【0046】

第１特性パラメーター初期特性推定値計算部１１６は、二次電池２２０の電流測定値Ｉ（ｔ）が０であった最後の時点ｔ＝ｊにおける電圧Ｖの測定値Ｖ（ｊ）である基準測定値Ｖ₀（ｊ）と、電圧Ｖの今回測定値Ｖ（ｋ）と、内部抵抗の今回初期特性推定値ｒ（０←ｋ）と、に基づき、関係式（２２０）にしたがって電圧Ｖの今回初期特性推定値Ｖ（０←ｋ）を計算する（関係式（１１２）参照）。

【0047】

Ｖ（０←ｋ）＝Ｖ₀（ｊ）－Ｉ（ｋ）・ｒ（０←ｋ） ‥（２２０）。

【0048】

二次電池２２０の電圧Ｖの測定値Ｖ（ｊ）が、第１特性モデルパラメーター初期特性推定値計算部１１６を構成する記憶部に記憶保持される。その後、二次電池２２０の電流測定値Ｉ（ｔ）が０になった場合、その時点における二次電池２２０の電圧Ｖの測定値Ｖ（ｊ）が、初期特性モデルにおける内部電源の今回起電力Ｖ₀（ｋ）として当該記憶部に上書き保存される。

【0049】

第１演算処理要素１１０は、二次電池２２０の特性を表わす複数の特性パラメーターｐ＝（ｐ₁，‥ｐ_N）のそれぞれの今回測定値ｐ（ｋ）＝（ｐ₁（ｋ），‥ｐ_N（ｋ））を、初期特性モデルを表わす関数Ｇに入力することにより、当該初期特性モデルの出力Ｇ（ｐ（ｋ））として複数の特性パラメーターｐのうち第１特性パラメーターｐ₁の今回初期特性推定値ｐ₁（０←ｋ）を計算する。当該関数Ｇは、例えば、関係式（１４）の右辺第２項に関係式（１２２）を代入することにより得られる関数である。

【0050】

第１指標値計算部１２１は、二次電池２２０の特性パラメーターｐの測定値ｐ（ｋ）＝（Ｖ（ｋ），Ｉ（ｋ），Ｔ（ｋ））に基づき、多変数関数ｆ（ｐ）にしたがって算出される値ｆ（Ｖ（ｋ），Ｉ（ｋ），Ｔ（ｋ））の累積値または時間積分値を第１指標値Ｆ₁として計算する（関係式（１２１）参照）。第２指標値計算部１２２は、電圧Ｖの初期特性推定値Ｖ（０←ｋ）ならびに電流Ｉの測定値Ｉ（ｋ）および温度Ｔの測定値Ｔ（ｋ）に基づき、同一の多変数関数ｆ（Ｖ，Ｉ，Ｔ）にしたがって算出される値ｆ（Ｖ（０←ｋ），Ｉ（ｋ），Ｔ（ｋ））の累積値または時間積分値を第２指標値Ｆ₂として計算する（関係式（１２２）参照）。

【0051】

図５には、初期状態における二次電池２２０の電圧Ｖ（第１特性パラメーターｐ₁（ｔ））の時間変化態様が実線で示され、劣化状態における二次電池２２０の電圧Ｖの時間変化態様が破線で示されている。図５に示されているように任意の期間が累積区間または積分区間［ｔ₁，ｔ₂］として採用されてもよい。電流Ｉが０になった最後の時点からの経過期間が所定期間以内であるような区間が当該累積区間として採用されてもよい。

【0052】

当該関数ｆに代えて、少なくとも電圧Ｖ（第１特性パラメーターｐ₁）を主変数とする関数ｆ₁（Ｖ）、ｆ₂（Ｖ，Ｉ）およびｆ₃（Ｖ，Ｔ）のうち少なくとも１つにしたがって、関係式（１２１）（１２２）と同様に当該関数（従変数）の値が累積または積分されることにより第１指標値Ｆ₁および第２指標値Ｆ₂が計算されてもよい。

【0053】

第３演算処理要素１３０（劣化度計算部）は、第１指標値Ｆ₁（ｉ）および第２指標値Ｆ₂（ｉ）に基づき、関係式（１３０）にしたがって二次電池２２０の劣化度Ｄ（ｉ）を計算する。

【0054】

（劣化状態判定方法）
前記構成の劣化状態判定装置により実行される対象の劣化状態判定方法について説明する。

【0055】

対象機器２００において、通電状態の制御装置２１０により第１条件が満たされているか否かが判定される（図４／ＳＴＥＰ２０２）。「第１条件」としては、対象機器２００において、入力インターフェース２０２を通じて二次電池２２０の劣化状態判定の要請があったこと、所定のアプリケーションソフトウェアが対象機器２００において起動されたこと、二次電池２２０の特性パラメーターの測定値が第１の変化態様を示したことなどの条件が採用される。

【0056】

第１条件が満たされていないと判定された場合（図４／ＳＴＥＰ２０２‥ＮＯ）、第１条件の充足性判定処理が再び実行される（図４／ＳＴＥＰ２０２）。第１条件の充足性判定処理（図４／ＳＴＥＰ２０２）は省略されてもよい。

【0057】

第１条件が満たされていると判定された場合（図４／ＳＴＥＰ２０２‥ＹＥＳ）、センサー群２２２の出力信号に基づき、二次電池２２０の特性を表わす複数の特性パラメーターｐの測定値ｐ（ｋ）＝（Ｖ（ｋ），Ｉ（ｋ），Ｔ（ｋ））が取得される（図４／ＳＴＥＰ２０４）。「ｋ」はサンプリング周期に応じた離散的な時刻を表わす指数である。二次電池２２０の電圧Ｖが第１特性パラメーターｐ₁として測定される。二次電池２２０の電流Ｉ（充電電流および放電電流を含む。）が第２特性パラメーターｐ₂として測定される。二次電池２２０の温度Ｔ（ハウジングの周囲温度または表面温度）が第３特性パラメーターｐ₃として測定される。この測定の過程で、電流Ｉが０になった最後の時点ｔ＝ｊにおける電圧Ｖの測定値Ｖ（ｊ）が、基準測定値Ｖ₀（ｊ）としてメモリーに記憶保持される。

【0058】

続いて、制御装置２１０により第２条件が満たされているか否かが判定される（図４／ＳＴＥＰ２０６）。「第２条件」としては、最後に第１条件が満たされたと判定された第１時点から所定時間が経過した第２時点に至ったこと、特性パラメーターｐ（ｋ）の測定結果を表わすデータの第１時点からの累計値が閾値に達したこと、センサー群２２２により測定される二次電池２２０の特性パラメーターの測定値が第２の変化態様を示したことなどの条件が採用される。

【0059】

第２条件が満たされていないと判定された場合（図４／ＳＴＥＰ２０６‥ＮＯ）、指数ｋが「１」だけ増やされ（図４／ＳＴＥＰ２０８）、そのうえで各特性パラメーターｐ_s（ｓ＝１，２，３）の測定値ｐ_s（ｋ）が取得される（図４／ＳＴＥＰ２０４）。この際、特性パラメーターｐ_sの測定値ｐ_s（ｋ）の今回時系列Ｐ_s（ｉ）＝｛ｐ_s（ｉ）｜ｉ＝ｋ，ｋ＋１，‥｝が累積的にメモリーに記憶保持される。

【0060】

第２条件が満たされていると判定された場合（図４／ＳＴＥＰ２０６‥ＹＥＳ）、特性パラメーターｐ_sの測定値ｐ_s（ｋ）の今回時系列Ｐ_S（ｉ）が、出力インターフェース１０２を構成する送信装置により対象機器２００から劣化状態判定装置１００に対して送信される（図４／ＳＴＥＰ２１０）。この際、二次電池２２０の規格または種類を識別するための識別子ＩＤ、および、電圧Ｖの基準測定値Ｖ₀（ｊ）も対象機器２００から劣化状態判定装置１００に対して送信される。当該送信に応じて指数ｉが「１」だけ増やされ（図４／ＳＴＥＰ２１２）、その上で第１条件の充足性判定以降の処理が再び実行される（図４／ＳＴＥＰ２０２→ＳＴＥＰ２０４→‥ＳＴＥＰ２１０参照）。

【0061】

第２条件の充足性判定処理（図４／ＳＴＥＰ２０６）が省略され、特性パラメーターｐ_sの今回測定値ｐ_s（ｋ）が、識別子ＩＤおよび電圧Ｖの基準測定値Ｖ₀（ｊ）とともに、対象機器２００から劣化状態判定装置１００に対して逐次的に送信されてもよい。

【0062】

劣化状態判定装置１００において、入力要素２０１が、特性パラメーターｐ_s（ｋ）の測定値の今回時系列Ｐ_s（ｉ）、識別子ＩＤおよび電圧Ｖの基準測定値Ｖ（ｊ）を受信する（図４／ＳＴＥＰ１０２）。

【0063】

第１演算処理要素１１０（第２初期特性モデルパラメーター計算部１１４）が、識別子ＩＤおよび特性パラメーターｐの測定値ｐ（ｋ）に対応する第１初期特性モデルパラメーターｑ₁（ｐ（ｋ））をメモリー（図２／第１初期特性モデルパラメーター記憶部１１２）から読み取る（図４／ＳＴＥＰ１０４）。第１初期特性モデルパラメーターｑ₁は、前記関係式（１２）における偏微分係数（∂ｒ／∂ｐ_s）である。

【0064】

第１演算処理要素１１０（第２初期特性モデルパラメーター計算部１１４）が、二次電池２２０の特性パラメーターｐの今回測定値ｐ（ｋ）と、第１初期特性モデルパラメーターｑ₁（ｐ（ｋ））と、に基づき、関係式（１２２）にしたがって、初期状態モデルにおける二次電池の内部抵抗の抵抗値ｒの今回初期特性推定値ｒ（０←ｋ）を第２初期特性モデルパラメーターｑ₂（０←ｋ）として計算する（図４／ＳＴＥＰ１０６）。

【0065】

第１演算処理要素１１０（第１特性パラメーター初期特性推定値計算部１１６）が、電圧Ｖの基準測定値Ｖ₀（ｊ）、電流Ｉの今回測定値Ｉ（ｋ）および初期モデルにおける内部抵抗の抵抗値ｒの今回初期特性推定値ｒ（０←ｋ）に基づき、関係式（２２０）にしたがって電圧Ｖの今回初期特性推定値Ｖ（０←ｋ）を計算する（図４／ＳＴＥＰ１０８）。

【0066】

第２演算処理要素１２０（第１指標値計算部１２１）が、少なくとも第１特性パラメーターｐ₁としての電圧Ｖの測定値Ｖ（ｋ）の今回時系列Ｖ（ｉ）に基づき、例えば関係式（１２１）にしたがって第１指標値Ｆ₁（ｉ）を計算する（図４／ＳＴＥＰ１１０）。第２演算処理要素１２０（第２指標値計算部１２２）が、少なくとも第１特性パラメーターｐ₁としての電圧Ｖの測定値Ｖ（ｋ）の今回時系列Ｖ（ｉ）に基づき、例えば関係式（１２２）にしたがって第２指標値Ｆ₂（ｉ）を計算する（図４／ＳＴＥＰ１１０）。

【0067】

第３演算処理要素１３０（劣化度計算部）が、第１指標値Ｆ₁（ｉ）および第２指標値Ｆ₂（ｉ）に基づき、例えば関係式（１３０）にしたがって二次電池２２０の劣化度Ｄ（ｉ）を計算する（図４／ＳＴＥＰ１１２）。

【0068】

第３演算処理要素１３０が、二次電池２２０の劣化度Ｄ（ｉ）に応じた劣化診断情報Ｉｎｆ（Ｄ（ｉ））を生成する（図４／ＳＴＥＰ１１４）。出力要素１０２を構成する送信装置により診断情報Ｉｎｆ（Ｄ（ｉ））が劣化状態判定装置１００から対象機器２００に対して送信される（図４／ＳＴＥＰ１１６）。

【0069】

対象機器２００において、入力インターフェース２０１を構成する受信装置が劣化診断情報Ｉｎｆ（Ｄ（ｉ））を受信する（図４／ＳＴＥＰ２２２）。出力インターフェース２０２を構成するディスプレイ装置に、劣化診断情報Ｉｎｆ（Ｄ（ｉ））が出力表示される（図４／ＳＴＥＰ２２４）。これにより、例えば図６に示されているように、二次電池２２０の劣化度Ｄ（ｉ）を示すグラフ表示のほか、「バッテリーの劣化度は３０％です。あと１５０日で交換することをおススメします。」といった劣化度Ｄ（ｉ）に応じた対処方法などに関するメッセージがディスプレイ装置に表示される。

【0070】

（効果）
本発明の劣化状態判定装置および劣化状態判定方法によれば、劣化状態の判定対象そのものである二次電池２２０の初期時点・初期期間における特性パラメーターｐの測定を省略し、図５および図７に示されているように任意の始点および終点を有する今回期間ｔ＝ｔ₁～ｔ₂における特性パラメーターｐの測定結果に基づいて二次電池２２０の劣化状態が判定される。このため、二次電池２２０の劣化状態判定に要する期間の短縮を図ることができる。

【0071】

これに対して、前記先行技術によれば、二次電池２２０の劣化状態を判定するためには、図８に示されているように一定期間ｔ＝ｔ_x1～ｔ_x2における二次電池２２０の特性パラメーターｐ＝（Ｖ，Ｉ，Ｔ）の測定値の今回時系列のみならず、今回以前（例えば初期状態）の一定期間ｔ＝ｔ₀₁～ｔ₀₂における二次電池２２０の特性パラメーターｐ＝（Ｖ，Ｉ，Ｔ）の測定値の過去時系列を用いることが必要になる。

【符号の説明】

【0072】

１００‥劣化状態判定装置、１０２‥入力要素、１０４‥出力要素、１１０‥第１演算処理要素、１１２‥第１初期特性モデルパラメーター記憶部、１１４‥第２初期特性モデルパラメーター初期特性推定値計算部、１１６‥第１特性パラメーター初期特性推定値計算部、１２０‥第２演算処理要素、１２１‥第１指標値計算部、１２２‥第２指標値計算部、１３０‥第３演算処理要素（劣化度計算部）、２００‥対象機器、２０２‥入力インターフェース、２０４‥出力インターフェース、２１０‥制御装置、２２０‥二次電池、２２２‥センサー群。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】