(57)【要約】  内部抵抗（ｒ）について設定された第１の分布Ｐに基づき、評価関数（Ｅ）が最小となる充放電パターン（ｕ０：ｔｆ）を生成し、事後分布（Ｑ）を近似する第２の分布（Ｒｔ）を生成する充放電パターン生成部（１４）、蓄電池（９００）に対する測定結果から算出される、時刻ごとの第３の分布（Ｓｔ）を生成する蓄電池診断部（１６）、および第３の分布の分散（ＶＳｔ［ｒ］）が第２の分布の分散（ＶＲｔ［ｒ］）より大きな場合は、充放電の大きさ（ｕ）の絶対値を増加させ、小さな場合は、絶対値を減少させる充放電パターン調整部（１７）、を備える。
2021014482.tif 000002
【特許請求の範囲】
【請求項１】
充放電の大きさが時系列で変化する充放電パターンを用いて、蓄電池に充放電させた際の電圧降下から前記蓄電池の内部抵抗を診断する蓄電池診断装置であって、
確率分布を状態変数とする離散時間確率システムにおいて、前記内部抵抗の確率分布における事前分布として設定された第１の分布に基づき、前記内部抵抗の推定精度と前記充放電に伴うコストの和を示す評価関数が最小となるよう前記充放電パターンを生成し、かつ、前記充放電パターンによって得られる前記第１の分布の事後分布を近似する第２の分布を生成する充放電パターン生成部、
前記充放電パターンを用いた前記蓄電池に対する充放電での測定結果から推定される、前記内部抵抗の第３の分布を生成する蓄電池診断部、および
前記時系列における時刻ごとに、前記第３の分布の分散と前記第２の分布の分散とを比較し、前記第３の分布の分散が前記第２の分布の分散より大きな場合は、前記充放電の大きさの絶対値を増加させ、小さな場合は、前記絶対値を減少させる充放電パターン調整部、
を備えたことを特徴とする蓄電池診断装置。
【請求項２】
前記充放電パターン調整部は、前記第３の分布の分散を前記第２の分布の分散で除した商を底とする正の実数のべき乗を、前記充放電の大きさに乗じることを特徴とする請求項１に記載の蓄電池診断装置。
【請求項３】
前記第１の分布、前記第２の分布、および前記第３の分布は、対数正規分布、または対数ロジスティック分布であることを特徴とする請求項１または２に記載の蓄電池診断装置。
【請求項４】
前記蓄電池診断部は、前記測定結果が得られる前に推定または算出された前記第３の分布を事前分布として前記測定結果が得られた時の事後分布を計算し、前記事後分布との分布間距離が最小となる、前記事前分布と同一の分布族に属する分布を、新たな第３の分布として更新することを特徴とする請求項１から３のいずれか１項に記載の蓄電池診断装置。
【請求項５】
前記充放電パターン調整部には、
当該時刻における前記第２の分布と、前記蓄電池診断部によって更新された前記第３の分布から、対数写像の値である第１の接ベクトルを計算する対数変換部、
前記対数変換部によって計算された第１の接ベクトルに係数を乗じ、接空間における修正値である第２の接ベクトルを計算する乗算部、および
前記乗算部によって計算された第２の接ベクトルと、前記充放電パターン生成部によって生成された充放電パターンの指数写像の値を前記充放電パターンの修正値として計算する指数変換部、
が設けられていることを特徴とする請求項４に記載の蓄電池診断装置。
【請求項６】
前記分布間距離として、カルバック・ライブラー情報量を用いることを特徴とする請求項４に記載の蓄電池診断装置。
【請求項７】
蓄電池に充放電させた際の電圧降下から前記蓄電池の内部抵抗を診断するための、前記充放電の大きさが時系列で変化する充放電パターンを設定する方法であって、
確率分布を状態変数とする離散時間確率システムにおいて、前記内部抵抗の確率分布における事前分布として設定された第１の分布に基づき、前記内部抵抗の推定精度と前記充放電に伴うコストの和を示す評価関数が最小となるよう前記充放電パターンを生成する充放電パターン生成ステップ、
前記充放電パターンによって得られる前記第１の分布の事後分布を近似する前記第１の分布と同一の分布族に属する第２の分布を生成する第２の分布生成ステップ、
前記充放電パターンを用いた前記蓄電池に対する充放電での測定結果から推定される、前記内部抵抗の第３の分布を生成する第３の分布生成ステップ、および
前記時系列における時刻ごとに、前記第３の分布の分散と前記第２の分布の分散とを比較し、前記第３の分布の分散が前記第２の分布の分散より大きな場合は、前記充放電の大きさの絶対値を増加させ、小さな場合は、前記絶対値を減少させる充放電パターン調整ステップ、
を含むことを特徴とする充放電パターン設定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、蓄電池診断装置および充放電パターン設定方法に関するものである。
【背景技術】
【０００２】
蓄電池は、使用状況、環境等の影響を受けて性能が劣化するため、交換時期を判断するために、蓄電池の内部抵抗、容量等の健全度を診断することが求められる。その際、診断用の充放電に伴う電力ロスを低減するため、組電池のうち、内部抵抗が最も高い電池をサンプル電池として選定し、サンプル電池のみに充放電試験を行う蓄電池診断装置（例えば、特許文献１参照。）が提案されている。あるいは、本来、緩和時間を長く設けて端子開放電圧を測定すべきところ、鉄道車両の運行パターンの特徴を利用して、短時間の緩和時間で充電率を算出する鉄道車両用特電池制御システム（例えば、特許文献２参照。）が提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００７−２４６４０号公報（段落００１６〜００６９、図１〜図７）
【特許文献２】特開２０１２−７８０９５号公報（段落００２１〜００４０、図１〜図５）
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、健全度を所望の精度で推定するために、どのような測定パターン、あるいは運用パターンを用いればよいかについては、いずれの文献でも検討されておらず、特許文献２に開示された技術では、運用パターンが不規則な態様には適用すらできない。つまり、診断に伴うロスを抑制するとともに、必要な精度で蓄電池の健全性を診断することは困難であった。【０００５】
本願は、上記のような課題を解決するための技術を開示するものであり、診断に伴うロスを抑制し、必要な精度で蓄電池の健全性を診断できる蓄電池診断装置、あるいは充放電パターン設定方法を得ることを目的とする。
【課題を解決するための手段】
【０００６】
本願に開示される蓄電池診断装置は、充放電の大きさが時系列で変化する充放電パターンを用いて、蓄電池に充放電させた際の電圧降下から前記蓄電池の内部抵抗を診断する蓄電池診断装置であって、確率分布を状態変数とする離散時間確率システムにおいて、前記内部抵抗の確率分布における事前分布として設定された第１の分布に基づき、前記内部抵抗の推定精度と前記充放電に伴うコストの和を示す評価関数が最小となるよう前記充放電パターンを生成し、かつ、前記充放電パターンによって得られる前記第１の分布の事後分布を近似する第２の分布を生成する充放電パターン生成部、前記充放電パターンを用いた前記蓄電池に対する充放電での測定結果から推定される、前記内部抵抗の第３の分布を生成する蓄電池診断部、および前記時系列における時刻ごとに、前記第３の分布の分散と前記第２の分布の分散とを比較し、前記第３の分布の分散が前記第２の分布の分散より大きな場合は、前記充放電の大きさの絶対値を増加させ、小さな場合は、前記絶対値を減少させる充放電パターン調整部、を備えたことを特徴とする。【０００７】
本願に開示される充放電パターン設定方法は、蓄電池に充放電させた際の電圧降下から前記蓄電池の内部抵抗を診断するための、前記充放電の大きさが時系列で変化する充放電パターンを設定する方法であって、確率分布を状態変数とする離散時間確率システムにおいて、前記内部抵抗の確率分布における事前分布として設定された第１の分布に基づき、前記内部抵抗の推定精度と前記充放電に伴うコストの和を示す評価関数が最小となるよう前記充放電パターンを生成する充放電パターン生成ステップ、前記充放電パターンによって得られる前記第１の分布の事後分布を近似する前記第１の分布と同一の分布族に属する第２の分布を生成する第２の分布生成ステップ、前記充放電パターンを用いた前記蓄電池に対する充放電での測定結果から推定される、前記内部抵抗の第３の分布を生成する第３の分布生成ステップ、および前記時系列における時刻ごとに、前記第３の分布の分散と前記第２の分布の分散とを比較し、前記第３の分布の分散が前記第２の分布の分散より大きな場合は、前記充放電の大きさの絶対値を増加させ、小さな場合は、前記絶対値を減少させる充放電パターン調整ステップ、を含むことを特徴とする。
【発明の効果】
【０００８】
本願に開示される蓄電池診断装置あるいは充放電パターン設定方法によれば、診断に伴うロスを抑制し、必要な精度で蓄電池の健全性を診断できる。
【図面の簡単な説明】
【０００９】
【図１】実施の形態１にかかる蓄電池診断装置の構成を説明するためのブロック図である。
【図２】実施の形態１にかかる蓄電池診断装置の動作、あるいは充放電パターン設定方法を示すフローチャートである。
【図３】実施の形態１にかかる蓄電池診断装置の演算処理を実行する部分の構成例を示すブロック図である。
【図４】実施の形態２にかかる蓄電池診断装置の充放電パターン調整部の詳細を説明するためのブロック図である。
【図５】実施の形態２にかかる蓄電池診断装置の充放電パターン生成部が生成する疑似コードを示す図である。
【発明を実施するための形態】
【００１０】
実施の形態１．
図１〜図３は、実施の形態１にかかる蓄電池診断装置の構成と動作、あるいは充放電パターン設定方法について説明するためのものであり、図１は蓄電池診断装置の全体構成を示すブロック図、図２は蓄電池診断装置の充放電パターンを設定する動作を説明するためのフローチャート、図３は蓄電池診断装置の演算処理を実行する部分の構成例を示すブロック図である。
【００１１】
本実施の形態にかかる蓄電池診断装置の構成、あるいは充放電パターン設定方法の詳細についての説明の前に、蓄電池の診断で実施する充放電パターンの調整に関する本願特有の考え方について、内部抵抗の診断を例にとって説明する。内部抵抗ｒは、蓄電池を流れる電流をｕ、その時の電圧降下の観測値をｙとして、オーム（Ohm）の法則（ｙ＝ｒｕ＋ｗ）によって推定できる。ただし、ｗは電圧の観測ノイズである。
【００１２】
ここで、ノイズｗが分散σ２のガウスノイズであるとすると、内部抵抗ｒの推定精度の上界はクラメール・ラオ（Cramer-Rao）の不等式（式（１））によって評価できる。
V[r]≧(-∂2/(∂r2)Logｐ(y│r))-1＝σ2/u2  ・・・（１）
ただし、Ｖ［ｒ］は内部抵抗ｒの分散、ｐ（ｙ│ｒ)は、内部抵抗がrであるときの電圧降下ｙの条件付き分布の確率密度を表す。
【００１３】
これにより電流ｕの絶対値が大きければ分散Ｖ［ｒ］が小さい、つまり高精度の推定が可能であり、電流ｕの絶対値が小さければ内部抵抗ｒの推定精度は低くならざるを得ないことが示される。さて、電池に電流を流すということは、内部抵抗ｒによるジュール（Joule）損失ｒｕ２が発生することを意味する。従って、内部抵抗ｒの推定精度を高めるためには、大きな消費電力、すなわちコストを要することがわかる。
【００１４】
ここで、内部抵抗ｒの推定精度と、推定に要するコストのトレードオフをとることを考える。すなわち、αを所定の定数として、ｒｕ２とαＶ［ｒ］の和を最小化するような電流ｕを求めたい。ただし、今、内部抵抗ｒは未知であり、分散Ｖ［ｒ］は未知の観測値ｙに依存するため、これらの事前分布に関する期待値をとったもの（式（２））を最小にするような電流ｕを求めればよい。
【００１５】
【数１】
2021014482.tif 000003
ここで、分散は内部抵抗の推定値ハットｒの分散であり、期待値は電圧降下の観測値ｙに関する期待値であることに注意する。
【００１６】
式（２）のＳ１（ｕ）を詳細に解析する。まず抵抗の事前分布をＰとすると、ベイズ(Bayes)の定理より、推定値ハットｒの分布（以下では事後分布と呼ぶ）Ｑは式（３）で与えられる。
dQ(r│y)=(p(y│r)dP(r))/(∫p(y│r)dP(r))  ・・・（３）
ただし、ｄＰ、ｄＱはそれぞれ事前分布Ｐ、事後分布Ｑの確率密度である。
【００１７】
Ｓ１（ｕ）の最適化のためには、事後分布Ｑの分散（すなわちＶ［ハットｒ］のことであるが、明確さのため、以下ではＶＱ［ｒ］のように書く）を、観測値ｙおよび電流ｕの関数として表現する必要がある。しかし、分母の積分評価は一般に難しく、また事前分布Ｐと事後分布Ｑは異なる分布族に属するため、事後分布Ｑを新たな事前分布として推定を繰り返すことができないという課題がある。
【００１８】
そこで、本願の各実施の形態にかかる蓄電池診断装置では、事前分布Ｐと同一の分布族に属する分布Ｒを用いて、分布間距離を示す、式（４）で定義されるカルバック・ライブラー（Kullback-Leibler）情報量を最小化することで、事後分布Ｑを近似する。
【００１９】
【数２】
2021014482.tif 000004
ただし、ｌｎは、自然対数を表す。例えば、事前分布Ｐとして対数平均μ、対数分散ω２の対数正規分布をとると、式（５）のようになり、分布Ｒの対数平均をν、対数分散をη２とすると、式（６）を最小にするようなν、ηを求める。
【００２０】
【数３】
2021014482.tif 000005
ただし、ＥＲとは分布Ｒに関する期待値を意味する。ここで期待値をとる関数はｌｎｒとｒの２次式であるから、式群（７）に示す対数正規分布の性質を用いれば、ＤＫＬ（Ｒ||Ｑ）は、解析的に式（８）と表すことができる。
【００２１】
【数４】
2021014482.tif 000006
したがって、ＤＫＬ（Ｒ||Ｑ）を最小にするν、ηは、数値的に容易に求めることができる。さらに、分布Ｒの観測値ｙ、電流ｕに関する偏微分を求めることができる。
【００２２】
今、時刻tにおける推定の事前分布をＲｔ（ただし、Ｒ０＝Ｐとする）とおき、事前分布Ｒｔのもとでの事後分布をＱｔとおくと、時刻ｔ＋１における事前分布Ｒｔ＋１は、ＤＫＬ（Ｒ||Ｑ）を最小にする分布Ｒである。ここで、事前分布Ｒｔ＋１は、事前分布Ｒｔ、電流ｕｔ、および観測値ｙｔの関数であるから、これを、式（９）と表す。
Rt+1＝F(Rt,ut,yt)  ・・・（９）
これは、確率分布を状態変数とする離散時間確率システムを示す。
【００２３】
本願の各実施の形態にかかる蓄電池診断装置は、この離散時間確率システムにおいて、評価関数Ｅ（式（１０））を最小にするような充放電パターンｕ0：ｔｆを生成する。
【００２４】
【数５】
2021014482.tif 000007
ただし、ｕ0：ｔｆとは、ｕｔをｔ＝０からｔｆ−１までまとめたものである。関数Ｓ、Ｌは、例えば内部抵抗ｒの診断においては、式群（１１）のように定義される。
【００２５】
【数６】
2021014482.tif 000008
さらに、実際の診断実行時において、時刻（または、時系列上の時点）ｔにおける診断結果が電流パターン生成時の分布とずれたときに、時刻ｔ＋１における電流ｕを動的に決定するためのフィードバック係数βを、同時に生成する。
【００２６】
これにより、本願の各実施の形態にかかる蓄電池診断装置１は、診断に要するコストを考慮して必要十分な精度で診断を行うことができる。さらに、実施の形態１において例示するように、診断実行時において、ある時刻ｔでの診断結果に応じて次の時刻ｔ＋１における電流を動的に決定することが可能となる効果を奏する。
【００２７】
実施の形態１にかかる蓄電池診断装置１は、図１に示すように、診断のための蓄電池９００への充放電を制御する充放電制御部１８と、充放電中の計測値に基づいて蓄電池９００を診断する蓄電池診断部１６と、診断結果表示用の診断結果表示部１９を備えている。さらに、充放電パターンｕ0：ｔｆを生成する充放電パターン生成部１４と、生成された充放電パターンｕ0：ｔｆを調整する充放電パターン調整部１７、および演算に必要なデータを格納する事前分布格納部１３と推定分布格納部１５を備えている。以下、詳細に説明する。
【００２８】
事前分布格納部１３は、蓄電池９００の健全度に関する内部抵抗ｒの事前分布である第１の分布Ｐを格納する。事前分布格納部１３は、後述する計算機ないしマイクロコントローラのメモリとして実現されてもよい。あるいは変更の必要がなければ電子的な回路として実現されてもよい。また、推定分布格納部１５は、推定開始時点において、上述した第１の分布Ｐを格納し、後述する計算機ないしマイクロコントローラのメモリとして実現されてもよい。
【００２９】
＜充放電パターン生成部＞
充放電パターン生成部１４は、事前分布格納部１３に格納された内部抵抗ｒの事前分布である第１の分布Ｐから、評価関数Ｅが最小となる充放電パターンｕ0：ｔｆを生成する。そして、生成した充放電パターンｕ0：ｔｆによって得られる蓄電池９００の健全度に関する事後分布Ｑを近似する第２の分布の時系列Ｒ０：ｔｆを生成する。
【００３０】
さらに、後述する蓄電池診断部１６から出力された診断結果Ｓｔが、第２の分布の時系列Ｒ０：ｔｆ中の当該時刻ｔの分布Ｒｔとずれた場合に、次の時刻の充放電パターンｕｔ＋１をどのように修正すべきかを表すフィードバック係数βを生成する。そして、生成した充放電パターンｕ0：ｔｆ、第２の分布の時系列Ｒ０：ｔｆのうちの時刻ｔにおける分布Ｒｔ、およびフィードバック係数βを充放電パターン調整部１７に出力する。
【００３１】
＜蓄電池診断部＞
蓄電池診断部１６は、推定分布格納部１５に格納された分布と、計測された蓄電池９００の電流または電圧、あるいはその双方から、蓄電池９００の診断結果から得られた健全度に関する分布である第３の分布Ｓｔを推定する。推定した第３の分布Ｓｔは、充放電パターン調整部１７に出力されるとともに、推定分布格納部１５に出力され、推定分布格納部１５内の格納内容を更新する。
【００３２】
＜充放電パターン調整部＞
充放電パターン調整部１７は、健全度に関する第３の分布Ｓｔと、充放電パターンｕ0：ｔｆと、時刻ｔにおける分布Ｒｔと、フィードバック係数βを受信し、次の時刻ｔ＋１における充放電パターンｕ0：ｔｆの修正値νｔを計算する。充放電制御部１８は、充放電パターン調整部１７が計算した充放電パターンｕ0：ｔｆにおける次の時刻ｔ＋１の充放電の大きさｕｔ＋１の修正値νｔにもとづき、蓄電池９００の充放電を行う。もちろん、修正値νｔが生成されない場合は、充放電パターン調整部１７を介して受信した充放電パターンｕ0：ｔｆで設定されたｕｔ＋１にもとづき、蓄電池９００の充放電を行う。
【００３３】
ここで、充放電パターンは、上述した蓄電池９００の充放電電流ｕを時系列（ｔ＝０，１，・・・，ｔｆ−１）上で変化させたものであってよい。あるいは、蓄電池９００の充放電電力の時系列上の変化であってもよい。また、定電圧操作をする場合は、電圧値の時系列上の変化であってもよい。さらに、電流ｕまたは電力と、蓄電池９００の環境温度の組の時系列上の変化であってよい。
【００３４】
健全度とは、少なくとも蓄電池９００の内部抵抗ｒを対象とするが、蓄電池９００の等価回路に含まれる容量（capacitance）、蓄電池９００の蓄電容量（charge capacity）のいずれか、またはその組み合わせが含まれていてもよい。また、蓄電池９００の健全度に関する分布（第３の分布Ｓｔ）は、例えば半無限区間［０，∞）を台とする分布であってよい。さらに具体的には指数分布、対数正規分布、対数ロジスティック分布、あるいはそれらの混合分布、それらに従う確率変数の和の分布であってよい。
【００３５】
上記構成を前提に、図２のフローチャートを参考にして、動作について説明する。
はじめに、事前分布格納部１３は、第１の分布Ｐの母数を記憶する。母数とは、例えば第１の分布Ｐが対数正規分布の場合、対数平均と、対数分散の組、あるいはそれに可逆な変換をかけた値であってよい。
【００３６】
充放電パターン生成部１４は、評価関数ＥのＪ（ｕ0：ｔｆ)を最小にする充放電パターンｕ0：ｔｆを生成する（ステップＳ１００）。評価関数Ｅにおいて、Ｊ（ｕ0：ｔｆ)の定義は期待値を含んでいるが、この最小化は、例えば、式（１２）に示す分布に独立に従うｙ０：ｔｆを、疑似乱数を用いて生成し、確率的最急降下法を用いることで実行できる。
O(yt)＝∫p(yt│r)dP(r)  ・・・（１２）
【００３７】
さらに、式（９）より、充放電パターンｕ0：ｔｆ、およびＲ０（＝Ｐ）から、健全度に関する事後分布Ｑを近似する第２の分布の時系列Ｒ０：ｔｆを計算する（ステップＳ１１０）。一方、実際の診断時に得られた第３の分布Ｓｔが、第２の分布の時系列Ｒ０：ｔｆに含まれる当該時刻ｔの第２の分布Ｒｔとずれた場合に、どのように充放電パターンｕｔ＋１を修正すべきかを表すフィードバック係数βについても計算する。フィードバック係数βは、本実施の形態においては、例えば所定の定数であってよい。
【００３８】
充放電パターンｕ０：ｔｆ、第２の分布の時系列Ｒ０：ｔｆ、およびフィードバック係数βは充放電パターン調整部１７に出力される。ここで、図示しない充放電開始指令を受けると、充放電制御部１８は、充放電パターン調整部１７が保持している充放電パターンｕ0：ｔｆにしたがい、蓄電池９００に対して充放電試験を実施する（ステップＳ１２０）。
【００３９】
蓄電池診断部１６は、計測前の第３の分布Ｓｔを事前分布として測定結果が得られた時の事後分布を計算し、事後分布との分布間距離（ＤＫＬ（Ｓｔ＋１||Ｗｔ））が最小となる、事前分布と同一の分布族に属する分布を、新たな第３の分布として更新する。つまり、充放電パターン生成部１４と同様に、事後分布Ｑを近似することで推定を繰りかえす。具体的には、例えば計測された電圧降下ｙｔに関し、分布間距離ＤＫＬ（Ｓｔ＋１||Ｗｔ）を最小にするような、第３の分布Ｓｔ＋１を求める（ステップＳ１３０）。ただし、dWt(r│y)は、式（１３）のとおりである。
dWt(r│y)＝(p(y│r)dSt(r))/(∫p(y│r)dSt(r))  ・・・（１３）
【００４０】
推定分布格納部１５は、第３の分布Ｓｔの母数を記憶する。母数とは、例えば、第３の分布Ｓｔが対数正規分布の場合、対数平均と、対数分散の組、あるいはそれに可逆な変換をかけた値であってよい。
【００４１】
充放電パターン調整部１７は、例えば、抵抗推定の場合、ＶＳｔ［ｒ］／ＶＲｔ［ｒ］を計算し、フィードバック係数βを用いて、充放電パターンｕ0：ｔｆの当該時刻ｔの次の時刻ｔ＋１における電流ｕｔ＋１の修正値νｔを計算する。このとき、推定分散ＶＳｔ［ｒ］がＶＲｔ［ｒ］より小さければ（ステップＳ１４０で「＜１」）、修正値νｔの絶対値を減少させる（ステップＳ１６０）。一方、推定分散ＶＳｔ［ｒ］がＶＲｔ［ｒ］より大きければ（ステップＳ１４０で「＞１」）修正値νｔの絶対値を増大させる（ステップＳ１５０）。
【００４２】
その際、例えば、式（１４）のように修正値νｔを定義し、フィードバック係数βとして、正の実数を設定する。すると、上述したステップＳ１４０〜Ｓ１６０のように、場合分け工程を設けなくても、一つの工程で、推定分散ＶＳｔ［ｒ］／ＶＲｔ［ｒ］の大小に応じて、自動的に修正値νｔの増大と減少を切り替えることができる。
νt＝((VSt[r]／VRt[r]))βｕｔ＋１  ・・・（１４）
【００４３】
修正値νｔが計算されると、時刻ｔ＋１における充放電のため、ｕｔ＋１として、充放電制御部１８に出力される。さらに、時刻ｔ＋１における充放電での診断結果に基づく推定分散ＶＳｔ＋１［ｒ］とＶＲｔ＋１［ｒ］とを比較するため、修正値νｔは、充放電パターン生成部１４にも出力される。充放電パターン生成部１４では、修正値νｔにより、時刻ｔ＋１における第２の分布Ｒｔ＋１が更新され、更新された第２の分布Ｒｔ＋１が充放電パターン調整部１７に返される（ステップＳ１７０）。これらの動作が、ひとつの充放電パターンｕ0：ｔｆによる充放電が完了するまで（ｔ＝ｔｆ−１になるまで）繰り返される（ステップＳ１８０）。
【００４４】
なお、充放電パターン生成部１４、蓄電池診断部１６、充放電パターン調整部１７は、例えば計算機またはマイクロコントローラのプログラムとして実現されてよい。あるいは特化された集積回路として実現されてもよい。その場合、例えば、診断結果表示部１９と充放電制御部１８を除いた演算処理を行う部分（事前分布格納部１３〜充放電パターン調整部１７）は、図３に示すように、プロセッサ１０１と記憶装置１０２を備えた一つのハードウェア１０によって構成することが考えられる。
【００４５】
事前分布格納部１３と推定分布格納部１５を記憶装置１０２とすると、図示していないが、ランダムアクセスメモリ等の揮発性記憶装置と、フラッシュメモリ等の不揮発性の補助記憶装置とを具備する。また、フラッシュメモリの代わりにハードディスクの補助記憶装置を具備してもよい。プロセッサ１０１は、記憶装置１０２から入力されたプログラムを実行する。この場合、補助記憶装置から揮発性記憶装置を介してプロセッサ１０１にプログラムが入力される。また、プロセッサ１０１は、演算結果等のデータを記憶装置１０２の揮発性記憶装置に出力してもよいし、揮発性記憶装置を介して補助記憶装置にデータを保存してもよい。
【００４６】
実施の形態２．
本実施の形態２においては、分布Ｒがリーマン多様体とみなせる場合の充放電パターン調整部の構成、および充放電パターン設定方法について説明する。図４と図５は、本実施の形態２にかかる蓄電池診断装置および充放電パターン設定方法について説明するためのもので、図４は充放電パターン調整部の詳細を説明するためのブロック図、図５は充放電パターン生成部が生成する疑似コードを示す図である。なお、本実施の形態２にかかる蓄電池診断装置の充放電パターン調整部以外の構成、あるいは充放電パターンの調整以外の動作については、実施の形態１の説明に用いた図１と同様であるため、構成の詳細な説明は繰り返さない。
【００４７】
実施の形態２にかかる蓄電池診断装置１の充放電パターン調整部１７は、図４に示すように、対数変換部１７１と乗算部１７２と指数変換部１７３に細分化される。ただし、その説明の前に、本実施の形態２にかかる蓄電池診断装置１あるいは充放電パターン設定方法にて扱う、分布Ｒについて説明する。
【００４８】
分布Ｒは、その母数を局所座標とすることで、多様体とみなすことができる。さらに、指数分布、対数正規分布、あるいは対数ロジスティック分布のように局所座標がベクトル空間の部分集合であるような分布の場合、可微分多様体となる。さらに、その接空間の内積を例えば、フィッシャー（Fisher）計量とすることで、分布Ｒをリーマン多様体（Riemannian manifold）とみなすことができる。
【００４９】
実施の形態１で説明した評価関数Ｅにおいて、Ｊ（ｕ0：ｔｆ)がＪ（ｕ0：ｔｆ)＝Ｊ０（Ｒ０,ｕ0：ｔｆ)であれば、部分和をとることで、式群（１５）のように再帰的に展開できる。
Jt(Rt,ut:tf):＝E[L(Rt,ut)+Jt+1(Rt+1,ut+1:tf)│Rt]，
Jtf(Rtf):＝S(Rtf)                        ・・・（１５）
ここで、動的計画法を適用すると、Ｊ（ｕ0：ｔｆ)の最小化は、式群（１６）のように展開できる。
【００５０】
【数７】
2021014482.tif 000009
ここで、Ｇ（Ｒ,ｕ,μ,ｗ）:＝Ｆ（Ｒ,ｕ,ｅμｕ＋ｗ）とおき、写像Ｇと写像Ｖｔ＋１の合成写像Ｈｔを点（Ｒバー,ｕバー,μバー,０）のまわりで展開すると、式（１７）を得る。
【００５１】
【数８】
2021014482.tif 000010
ここで、ＥｘｐＲバーＸ、ＥｘｐｕバーＵは指数写像であり、∇ＸＨｔ等はＨｔの勾配、∇Ｘ，ＸＨｔ等はヘッシアン（Hessian）である。ここから、式（１８）のようになる。
【００５２】
【数９】
2021014482.tif 000011
ただし、ν（Ｒ）とは、分布Ｒの対数分散を表す。
さらに、式（１９）とおくと、点（Ｒバー，ｕバー）の周りでの、式（２０）に示す２次近似を得る。
【００５３】
【数１０】
2021014482.tif 000012
ただし、Γｉｊｋを分布のなす多様体の接続係数として、式群（２１）のようになる。
【００５４】
【数１１】
2021014482.tif 000013
このように、Ｑｔの勾配およびヘッシアンは、Ｖｔ＋１の勾配およびヘッシアンを用いて計算できる。また、Ｑｔ（Ｒ，ｕ）の最小を与えるｕは、近似的に式（２２）で与えられる。
【００５５】
【数１２】
2021014482.tif 000014
ただし、ＬｏｇＲバーＲは対数写像である。
【００５６】
さらに、ＶｔはＱｔの最小値であることから、Ｖｔの勾配とヘッシアンは、式群（２３）のようになる。
∇i´Vt＝∇i´Qt-∇kQt・∇klQt・∇i´lQt,
∇i´j´Vt＝∇i´j´Qt-∇ki´Qt・∇klQt・∇i´lQt  ・・・（２３）
これは、Ｑｔの勾配とヘッシアンから求められる。ここで、Ｖｔｆ＝Ｓに注意すると、Ｑｔの勾配とヘッシアンを時刻ｔｆから再帰的に求められることがわかる。
【００５７】
従って、任意の充放電パターンｕ０：ｔｆ（ｋ）があったとき、その充放電パターンによる分布Ｒ０：ｔｆ（ｋ）を計算し、充放電パターンｕ０：ｔｆ（ｋ）と得られたＲ０：ｔｆ（ｋ）からＱｔの勾配とヘッシアンを計算する。そして、それを用いて新たな充放電パターンｕ０：ｔｆ（ｋ＋１）を計算することで、評価関数ＥにおけるＪ（ｕ0：ｔｆ)を最小化する充放電パターンを求めることができる。
【００５８】
上述した事項を前提とし、充放電パターン生成部１４が生成する擬似コードについて図５を用いて説明する。１行目において、任意の初期充放電パターンｕ０：ｔｆ（０）を用意する。２行目から５行目にかけて、初期充放電パターンに対応する第２の分布の時系列を計算する。ただし、ＥＮはノイズｗに関する期待値を表す。８行目から１１行目では、Ｑｔの勾配とヘッシアンを再帰的に計算する。１２行目から１８行目にかけて、Ｑｔの勾配とヘッシアンをもとに充放電パターンｕ０：ｔｆ（ｋ）を更新する。以上をＱｔのＵに関する勾配が十分に小さくなるまで繰り返す。
【００５９】
アルゴリズムが終了したとき、充放電パターン生成部１４は充放電パターンｕ０：ｔｆ（ｋ）、そのときの第２の分布の時系列Ｒ０：ｔｆ（ｋ）、および式（２４）で表されるＫｔによるフィードバック係数の時系列Ｋ０：ｔｆを出力する。
Kt＝-(∇U,UQt(k))-1∇U,XQt(k)  ・・・（２４）
【００６０】
これらを踏まえ、図４に示す充放電パターン調整部１７について説明する。対数変換部１７１は、蓄電池診断部１６が実際の診断によって算出した第３の分布Ｓｔと、充放電パターン生成部１４が生成した第２の分布Ｒｔから、第１の接ベクトルＸ＝ＬｏｇＲｔＳｔを計算する。乗算部１７２は充放電パターン生成部１４が生成したフィードバック係数Ｋｔと、対数変換部１７１が計算した第１の接ベクトルの値Ｘから、第２の接ベクトルＵ＝Ｋｔ・Ｘを計算する。指数変換部１７３は、乗算部１７２が計算した第２の接ベクトルＵと、充放電パターン生成部１４が生成した充放電パターンｕｔから、指数写像νｔ＝ＥｘｐｕｔＵを計算する。
【００６１】
充放電制御部１８は、指数変換部１７３から出力された指数写像νｔを、次の充放電パターンｕ０：ｔｆ（ｋ＋１）の修正値として、蓄電池９００を充放電する。このように構成することで、蓄電池９００の健全度の推定精度と、推定に要するコストと重み付け和を最小にするように蓄電池９００を診断することができる。さらに実際の診断時において、蓄電池９００の健全度に関する分布が、計画時と異なる場合においても、重み付け和を最小にするように充放電パターンを修正することができる。
【００６２】
実施の形態３．
本実施の形態３においては、上記実施の形態２と同じ構成である。しかし、充放電パターンがユークリッド空間の元である場合を想定し、実施の形態２での説明よりも指数写像を容易に計算する例について説明する。
【００６３】
上述したように、充放電パターンｕｔがユークリッド空間の元であるとき、指数写像は、ＥｘｐｕＵ＝ｕ＋Ｕにより、容易に計算できる。一方、分布Ｒはユークリッド空間ではないため、その対数写像は容易に計算できるとは限らない。実際、対数正規分布においては、指数写像の解析的な式は存在するが、その逆関数である対数写像は反復演算によって求めることになる。
【００６４】
そこで、実施の形態３にかかる蓄電池診断装置１では、Ｒのなすリーマン多様体の接続を、通常のレヴィ・チヴィタ（Levi-Civita）接続に代えて、双対平坦接続（∇,∇＊）を用いる。このとき、関数ｆ（Ｒ）のヘッシアンは、接続∇＊の接続係数Δｉｊｋを用いて式（２５）となり、実施の形態２で示した式群（２１）のΓｉｊｋもまた、Δｉｊｋによって置き換えられる。
∇ijf＝∂ijf-Δijk∂kf  ・・・（２５）
【００６５】
このとき、対数写像が式（２６）のように容易に計算できる。
LogRS＝(ξ-ν,ζ2-η2)  ・・・（２６）
ここで、ξとζ２は、それぞれＳの対数平均と対数分散を表す。本実施の形態３のように、充放電パターンｕｔがユークリッド空間の元であることに着目することで、計算量の削減が達成される。
【００６６】
なお、本願は、様々な例示的な実施の形態および実施例が記載されているが、１つ、または複数の実施の形態に記載された様々な特徴、態様、および機能は特定の実施の形態の適用に限られるのではなく、単独で、または様々な組み合わせで実施の形態に適用可能である。従って、例示されていない無数の変形例が、本願明細書に開示される技術の範囲内において想定される。例えば、少なくとも１つの構成要素を変形する場合、追加する場合または省略する場合、さらには、少なくとも１つの構成要素を抽出し、他の実施の形態の構成要素と組み合わせる場合が含まれるものとする。
【００６７】
例えば、充放電パターンｕ０：ｔｆと第２の分布Ｒｔの双方を充放電パターン生成部１４で生成する例を示したが、これに限ることはなく、別々の部位で生成してもよい。また、フィードバック係数βについては、充放電パターン生成部１４で設定する例を挙げたが、これに限ることはなく、例えば、充放電パターン調整部１７で設定してもよい。さらには、修正値νｔ（充放電パターンｕｔ＋１）に対する第２の分布Ｒｔ＋１の更新は、充放電パターン調整部１７で行うようにしてもよい。
【００６８】
また、上述した蓄電池診断装置１あるいは充放電パターン設定方法において、ある時刻ｔにおける第３の分布の分散ＶＳｔ［ｒ］と第２の分布の分散ＶＳｔ［ｒ］との比較結果による調整対象として、以下のような対象が考えられる。例えば、実施の形態１で示すように、充放電パターンｕ０：ｔｆにおける次の時刻ｔ＋１の充放電の大きさｕｔ＋１を対象とすれば、動的な調整が可能となる。あるいは、実施の形態２で示すように、今回の充放電パターンｕ０：ｔｆ（ｋ）に対して、新たな充放電パターンｕ０：ｔｆ（ｋ＋１）における同じ時刻ｔの充放電の大きさｕｔ（ｋ＋１）を対象とすれば、次回の充放電において、系統的に好適な充放電が行える。しかしこれに限ることはなく、その時刻ｔに関連する充放電の大きさｕを調整対象とすればよい。
【００６９】
以上のように、各実施の形態にかかる蓄電池診断装置１によれば、充放電の大きさｕが時系列で変化する充放電パターンｕ０：ｔｆを用いて、蓄電池９００に充放電させた際の電圧降下ｙから蓄電池９００の内部抵抗ｒを診断する蓄電池診断装置１であって、確率分布を状態変数とする離散時間確率システムにおいて、内部抵抗ｒの確率分布における事前分布として設定された第１の分布Ｐに基づき、内部抵抗ｒの推定精度と充放電に伴うコストの和を示す評価関数Ｅが最小となるよう充放電パターンｕ０：ｔｆを生成し、かつ、充放電パターンｕ０：ｔｆによって得られる第１の分布Ｐの事後分布Ｑを近似する第２の分布Ｒｔを生成する充放電パターン生成部１４、充放電パターンｕ０：ｔｆを用いた蓄電池９００に対する充放電での測定結果から推定される、内部抵抗ｒの第３の分布Ｓｔを生成する蓄電池診断部１６、および時系列における時刻ｔごとに、第３の分布の分散ＶＳｔ［ｒ］と第２の分布の分散ＶＲｔ［ｒ］とを比較し、第３の分布の分散ＶＳｔ［ｒ］が第２の分布の分散ＶＲｔ［ｒ］より大きな場合は、充放電の大きさｕの絶対値を増加させ、小さな場合は、絶対値を減少させる充放電パターン調整部１７、を備えるように構成したので、診断に伴うロスを抑制し、必要な精度で蓄電池９００の健全性を診断することができる。
【００７０】
とくに、充放電パターン調整部１７は、第３の分布の分散ＶＳｔ［ｒ］を第２の分布の分散ＶＲｔ［ｒ］で除した商を底とする正の実数（フィードバック係数β）のべき乗（（(VSt[r]／VRt[r]))β）を、充放電パターンｕ０：ｔｆにおける次の時刻ｔ＋１の充放電の大きさｕｔ＋１に乗じる（式（１４））ようにすれば、第３の分布の分散ＶＳｔ［ｒ］と第２の分布の分散ＶＲｔ［ｒ］の比に応じて自動的に充放電の大きさｕｔ＋１の絶対値を増減させることができる。
【００７１】
また、第１の分布Ｐ、第２の分布Ｒｔ、および第３の分布Ｓｔは、対数正規分布、または対数ロジスティック分布であれば、可微分多様体となり、フィッシャー軽量で分布Ｒをリーマン多様体とみなすことができるので、Ｑｔの勾配とヘッシアンを時刻ｔｆから再帰的に求めることができる。
【００７２】
蓄電池診断部１６は、測定結果が得られる前に推定または算出された第３の分布Ｓｔを事前分布として測定結果が得られた時の事後分布を計算し、事後分布との分布間距離（ＤＫＬ（Ｓｔ＋１||Ｗｔ））が最小となる、事前分布Ｒと同一の分布族に属する分布を、新たな第３の分布として更新するように構成すれば、分布Ｓｔが対数正規分布の場合、事前分布として、対数平均と対数分散の組、あるいはそれに可逆な変換をかけた値を母数として、記憶しておけばよい。
【００７３】
充放電パターン調整部１７には、当該時刻ｔにおける第２の分布Ｒｔと、蓄電池診断部１６によって更新された第３の分布Ｓｔから、対数写像の値である第１の接ベクトルＸ＝ＬｏｇＲｔＳｔを計算する対数変換部１７１、対数変換部１７１によって計算された第１の接ベクトルＸに係数Ｋｔを乗じ、接空間における修正値である第２の接ベクトルＵを計算する乗算部１７２、および乗算部１７２によって計算された第２の接ベクトルＵと、充放電パターン生成部１４によって生成された充放電パターンｕ０：ｔｆの指数写像の値ＥｘｐｕｔＵを充放電パターンの修正値νｔとして計算する指数変換部１７３、が設けられているように構成すれば、測定結果で得られた第３の分布が計画時と異なっている場合においても、容易に評価関数Ｅを最小化し、診断に伴うロスを抑制し、必要な精度で蓄電池９００の健全性を診断することができる。
【００７４】
分布間距離として、カルバック・ライブラー情報量を用いるようにすれば、事前分布Ｐと事後分布Ｑが異なる分布族に属しており、推定の繰り返しが困難であっても、事前分布Ｐと同一の分布族に属するＲを用いて、事後分布Ｑを近似することで、推定の繰り返しが可能となる。
【００７５】
また、各実施の形態にかかる充放電パターン設定方法によれば、蓄電池９００に充放電させた際の電圧降下ｙから蓄電池９００の内部抵抗ｒを診断するための、充放電の大きさｕが時系列で変化する充放電パターンｕ０：ｔｆを設定する方法であって、確率分布を状態変数とする離散時間確率システムにおいて、内部抵抗ｒの確率分布における事前分布として設定された第１の分布Ｐに基づき、内部抵抗ｒの推定精度と充放電に伴うコストの和を示す評価関数Ｅが最小となるよう充放電パターンｕ０：ｔｆを生成する充放電パターン生成ステップ（ステップＳ１００）、充放電パターンｕ０：ｔｆによって得られる第１の分布Ｐの事後分布Ｑを近似する第１の分布Ｐと同一の分布族に属する第２の分布Ｒｔを生成する第２の分布生成ステップ（ステップＳ１１０）、充放電パターンｕ０：ｔｆを用いた蓄電池９００に対する充放電での測定結果から推定される、内部抵抗ｒの第３の分布Ｓｔを生成する第３の分布生成ステップ（ステップＳ１２０、ステップＳ１３０）、および時系列における時刻ｔごとに、第３の分布の分散ＶＳｔ［ｒ］と第２の分布の分散ＶＳｔ［ｒ］とを比較し、第３の分布の分散ＶＳｔ［ｒ］が第２の分布の分散ＶＳｔ［ｒ］より大きな場合は、充放電の大きさｕの絶対値を増加させ、小さな場合は、絶対値を減少させる充放電パターン調整ステップ（ステップＳ１４０〜Ｓ１６０）、を含むように構成したので、診断に伴うロスを抑制し、必要な精度で蓄電池９００の健全性を診断することができる。
【符号の説明】
【００７６】
１：蓄電池診断装置、  １３：事前分布格納部、  １４：充放電パターン生成部、  １５：推定分布格納部、  １６：蓄電池診断部、  １７：充放電パターン調整部、  １７１：対数変換部、  １７２：乗算部、  １７３：指数変換部、  ９００：蓄電池、  ＤＫＬ（Ｒ||Ｑ）：分布間距離、  ＤＫＬ（Ｓｔ＋１||Ｗｔ）：分布間距離、  Ｅ：評価関数、  Ｐ：事前分布（第１の分布）、  ｒ：内部抵抗、  Ｒｔ：第２の分布、  Ｓｔ：診断結果（第３の分布）、  ｕ：充放電の大きさ、  ｕ０：ｔｆ：充放電パターン、  ｙ：電圧降下。
【図１】
2021014482.tif 000015
【図２】
2021014482.tif 000016
【図３】
2021014482.tif 000017
【図４】
2021014482.tif 000018
【図５】
2021014482.tif 000019
【国際調査報告】
2021014482.tif 000020
2021014482.tif 000021
2021014482.tif 000022
2021014482.tif 000023