(57)【要約】【課題】通常稼働している蓄電システムから得られる情報から、蓄電装置を評価する蓄電池評価装置、蓄電池評価方法及びコンピュータプログラムを提供する。【解決手段】本実施形態の蓄電池評価装置は、充電又は放電を指示する充放電指令値に応じて充放電する蓄電装置の稼働中の計測データに基づき、前記蓄電装置の充電量が閾値条件を満たす複数の時刻を検出し、検出した前記複数の時刻に基づき、複数の時間区間を選出する時間区間選出部と、前記時間区間に属する前記充放電指令値に基づき、前記時間区間の充放電効率を算出する効率算出部と、前記充放電効率の集合を複数のクラスタに分割し、前記充放電指令値及び前記計測データの少なくとも一方に基づき、前記クラスタの特徴情報を生成するクラスタリング部と、前記クラスタに属する前記充放電効率の推移データと、前記クラスタの特徴情報とを表示するための出力データを生成する出力部と、を備える。【選択図】図１
2020109720.tif 000002
【特許請求の範囲】
【請求項１】
充電又は放電を指示する充放電指令値に応じて充放電する蓄電装置の稼働中の計測データに基づき、前記蓄電装置の充電量が閾値条件を満たす複数の時刻を検出し、検出した前記複数の時刻に基づき、複数の時間区間を選出する時間区間選出部と、
前記時間区間に属する前記充放電指令値に基づき、前記時間区間の充放電効率を算出する効率算出部と、
前記充放電効率の集合を複数のクラスタに分割し、前記充放電指令値及び前記計測データの少なくとも一方に基づき、前記クラスタの特徴情報を生成するクラスタリング部と、
前記クラスタに属する前記充放電効率の推移データと、前記クラスタの特徴情報とを表示するための出力データを生成する出力部と、
を備えた蓄電池評価装置。
【請求項２】
前記出力データを表示する表示部
を備えた請求項１に記載の蓄電池評価装置。
【請求項３】
前記閾値条件は、前記充電量が閾値に達すること、又は前記閾値に対して一定の範囲に入ることである
請求項１又は２に記載の蓄電池評価装置。
【請求項４】
前記閾値は、前記蓄電装置の充電量が前記蓄電装置の定格容量の所定割合である
請求項３に記載の蓄電池評価装置。
【請求項５】
前記時間区間選出部は、前記計測データに基づき前記充電量の分布を生成し、前記分布に基づき前記閾値条件を決定する
請求項１〜４のいずれか一項に記載の蓄電池評価装置。
【請求項６】
前記推移データは、前記充放電効率を時刻と関連付けたデータを含む
請求項１〜５のいずれか一項に記載の蓄電池評価装置。
【請求項７】
前記推移データは、前記充放電効率を前記蓄電装置の累積放電量又は累積充電量と関連づけたデータを含む
請求項１〜６のいずれか一項に記載の蓄電池評価装置。
【請求項８】
前記クラスタリング部は、前記クラスタを近似する関数を計算し、
前記推移データは、前記関数のグラフを含む
請求項１〜７のいずれか一項に記載の蓄電池評価装置。
【請求項９】
前記クラスタリング部は、前記時間区間に属する前記充放電指令値及び前記計測データの少なくとも一方の標準偏差を計算し、前記標準偏差の統計情報を、前記クラスタの特徴情報として生成する
請求項１〜８のいずれか一項に記載の蓄電池評価装置。
【請求項１０】
前記クラスタリング部は、前記クラスタの特徴情報として、前記クラスタを近似する直線の傾き情報を計算する
請求項１〜９のいずれか一項に記載の蓄電池評価装置。
【請求項１１】
前記時間区間選出部は、前記時刻間の時間差が所定値以上である区間を前記時間区間として選出する
請求項１〜１０のいずれか一項に記載の蓄電池評価装置。
【請求項１２】
前記時刻間の時間差が所定値以上である区間は、前記計測データのサンプル数が前記所定値以上の区間である
請求項１１に記載の蓄電池評価装置。
【請求項１３】
前記時間区間選出部は、前記充放電指令値が前記充電及び前記放電のいずれも行わないことを指示している時間長さの合計が所定値以下である前記時刻間の区間を、前記時間区間として選出する
請求項１〜１２のいずれか一項に記載の蓄電池評価装置。
【請求項１４】
前記時間区間選出部は、前記時刻間の充電量の平均値と、閾値との差が所定値以下である前記時刻間の区間を、前記時間区間として選出する
請求項１〜１３のいずれか一項に記載の蓄電池評価装置。
【請求項１５】
前記閾値条件は、前記充電量が第１閾値又は第２閾値に達することであり、
前記第１閾値および前記第２閾値のそれぞれについて前記時間区間選出部、前記効率算出部、前記クラスタリング部、前記出力部の処理を行い、
前記第１閾値に基づく前記クラスタの特徴情報は前記第１閾値を特定する情報を含み、前記第２閾値に基づく前記クラスタの特徴情報は前記第２閾値を特定する情報を含む
請求項１〜１４のいずれか一項に記載の蓄電池評価装置。
【請求項１６】
前記効率算出部は、前記充放電効率を、前記充放電効率が前記算出された前記時間区間に属する計測データが示す温度に応じて補正する
請求項１〜１５のいずれか一項に記載の蓄電池評価装置。
【請求項１７】
充電又は放電を指示する充放電指令値に応じて充放電する蓄電装置の稼働中の計測データに基づき、前記蓄電装置の充電量が閾値条件を満たす複数の時刻を選出するステップと、
検出した前記複数の時刻に基づき、複数の時間区間を選出するステップと、
前記時間区間に属する前記充放電指令値に基づき、前記時間区間の充放電効率を算出するステップと、
前記充放電効率の集合を複数のクラスタに分割し、前記充放電指令値及び前記計測データの少なくとも一方に基づき、前記クラスタの特徴情報を生成するステップと、
前記クラスタに属する前記充放電効率の推移データと、前記クラスタの特徴情報とを表示するための出力データを生成するステップと、
を備えた蓄電池評価方法。
【請求項１８】
充電又は放電を指示する充放電指令値に応じて充放電する蓄電装置と、
蓄電池評価装置とを備え、
前記蓄電池評価装置は、
前記蓄電装置の稼働中の計測データに基づき、前記蓄電装置の充電量が閾値条件を満たす複数の時刻を検出し、検出した前記複数の時刻に基づき、複数の時間区間を選出する時間区間選出部と、
前記時間区間に属する前記充放電指令値に基づき、前記時間区間の充放電効率を算出する効率算出部と、
前記充放電効率の集合を複数のクラスタに分割し、前記クラスタに属する前記計測データに基づき、前記クラスタの特徴情報を生成するクラスタリング部と、
前記クラスタに属する前記充放電効率の推移データと、前記クラスタの特徴情報とを表示するための出力データを生成する出力部と、
を備えた蓄電池評価システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、蓄電池評価装置、蓄電池評価方法及び蓄電池評価システムに関する。
【背景技術】
【０００２】
蓄電池の劣化は、蓄電池の長期に渡る利用や充放電の繰り返しにより進む。蓄電池の劣化が進むと、蓄電池に定格電力量を充電できなくなる。したがって、蓄電池の劣化具合を随時把握しておかないと、予定量の充放電が行えず、システム全体の効率が低下する可能性がある。
【０００３】
蓄電池の劣化状態を評価する方法として、時間をかけて完放電と満充電を行い、この間の測定結果に基づき、蓄電池の劣化状態を評価する方法がある。また、EV（Electric Vehicle）や携帯電話のような連続充電を行うときに測定を行い、測定結果に基づき蓄電池の劣化状態を評価する方法がある。これらの方法は、特殊な充放電を実行している時、またはオフライン時（蓄電池の非運転時）を前提としているものであった。
【０００４】
一方、蓄電池の通常運転中の充放電電圧から評価を行う電圧標準偏差法もある。しかしながら、この方法は、事前に充放電電圧標準偏差と電池健全度とを対応づけるためのデータ収集が必要である。このため、蓄電池の診断を開始するまでに時間を要するという課題がある。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第４６６８０１５号公報
【特許文献２】特許第５０１７０８４号公報
【特許文献３】特許第６１３４４３８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明の実施形態は、通常稼働している蓄電システムから得られる情報を用いて蓄電装置を評価する蓄電池評価装置、蓄電池評価方法及び蓄電池評価システムを提供する。
【課題を解決するための手段】
【０００７】
本実施形態の蓄電池評価装置は、充電又は放電を指示する充放電指令値に応じて充放電する蓄電装置の稼働中の計測データに基づき、前記蓄電装置の充電量が閾値条件を満たす複数の時刻を検出し、検出した前記複数の時刻に基づき、複数の時間区間を選出する時間区間選出部と、前記時間区間に属する前記充放電指令値に基づき、前記時間区間の充放電効率を算出する効率算出部と、前記充放電効率の集合を複数のクラスタに分割し、前記クラスタに属する前記計測データに基づき、前記クラスタの特徴情報を生成するクラスタリング部と、前記クラスタに属する前記充放電効率の推移データと、前記クラスタの特徴情報とを表示するための出力データを生成する出力部と、を備える。
【図面の簡単な説明】
【０００８】
【図１】本実施形態に係る蓄電池評価システムの基本構成を示すブロック図。
【図２】蓄電池の構成例を示す図。
【図３】電池モジュールの構成の一例を示す図。
【図４】データ記憶部に記憶された充放電履歴データの例を示す図。
【図５】線形補間により充電量が閾値θになった時刻を特定する例を示す図。
【図６】線形補間により充電量が閾値θになった時刻を特定する例を示す図。
【図７】充放電効率を計算する具体例を示す図。
【図８】充放電履歴データの他の例を示す図。
【図９】図８の充放電履歴データを、横軸を時刻、縦軸を充電量とする座標系に展開してグラフ表示した例を示す図。
【図１０】データを、時刻及び充放電効率を軸とする座標系にプロットした例を示す図。
【図１１】クラスタリングの例を模式的に示す図。
【図１２】長い測定期間で取得したデータをプロットした例を示す図。
【図１３】長い測定期間で取得したデータをプロットした例を示す他の図。
【図１４】閾値θとしてθ１とθ２を使った場合に表示される画面の例を模式的に示す図。
【図１５】本実施形態に係る蓄電池評価装置の動作の一例のフローチャート。
【図１６】図１２の各点の充放電効率を代表温度に応じて補正した例を示す図。
【図１７】図１３の各点の充放電効率を代表温度に応じて補正した例を示す図。
【図１８】探索開始時刻から探索終了時刻までの間の充電量のヒストグラムを模式的に示す図。
【図１９】本発明の実施形態に係る蓄電池評価装置のハードウェア構成例を示す図。
【発明を実施するための形態】
【０００９】
以下、図面を参照しながら、本発明の実施形態について説明する。
【００１０】
図１は、本実施形態に係る蓄電池評価システムの基本構成を示すブロック図である。蓄電池評価システムは、評価対象となる蓄電池を含む蓄電システム（ＥＳＳ：Ｅｎｅｒｇｙ  Ｓｔｏｒａｇｅ  Ｓｙｓｔｅｍ）１０１と、蓄電池を評価する蓄電池評価装置２０１とを備える。蓄電池は二次電池とも呼ばれる。以下では蓄電池に呼び方を統一する。
【００１１】
蓄電システム１０１は、通信ネットワークを介して、ＳＣＡＤＡ３０１に接続されている。ＳＣＡＤＡ３０１（Ｓｕｐｅｒｖｉｓｏｒｙ  Ｃｏｎｔｒｏｌ  Ａｎｄ  Ｄａｔａ  Ａｃｑｕｉｓｉｔｉｏｎ）は、地域内に存在する様々な蓄電システム（ＥＳＳ）を１つの大型ＥＳＳに見立て、個々のＥＳＳに、時間に応じて充電指令値及び放電指令値の少なくとも一方（以下、充放電指令値）を送る。充放電指令値には、充放電の実行時刻が付加されていてもよい。あらかじめ時刻同期がされている場合などは、充放電指令値に実行時刻を付加しない構成も可能である。ＳＣＡＤＡ３０１は、電力系統４０１に対する個々のＥＳＳの充電及び放電の少なくとも一方（以下、充放電）を制御する。ＳＣＡＤＡ３０１は、一例として、電力会社の中央給電指令所など、上位のエネルギー・マネジメント・システムからの指令、または下位の配電側の個々のエネルギー管理システム（Ｅｎｅｒｇｙ  Ｍａｎａｇｅｍｅｎｔ  Ｓｙｓｔｅｍ）からの指令、またはこれらの両方等に基づき、個々のＥＳＳの充放電の制御を行う。
【００１２】
蓄電システム１０１は、ＳＣＡＤＡ３０１からの充放電指令値に応じて、電力系統４０１との間で蓄電池の充放電を行う機能を有する。蓄電システム１０１は、制御部１１１と、蓄電池１１２と、交直変換機１１３とを備える。交直変換機１１３は、電力系統４０１に接続されている。交直変換機１１３と電力系統４０１との間に、変圧器が存在してもよい。
【００１３】
蓄電システム１０１の制御部１１１は、ＳＣＡＤＡ３０１から充放電指令値を受信し、当該充放電指令値および蓄電池１１２の充電量に基づき、交直変換機１１３に対する実行可能な充放電指令値を生成する。ＳＣＡＤＡ３０１の充放電指令値と区別するため、便宜上、制御部１１１が生成した充放電指令値を、充放電命令と称する場合がある。制御部１１１は、生成した充放電命令に基づき交直変換機１１３を制御して、蓄電池１１２の充放電を行う。
【００１４】
交直変換機１１３は、電力系統４０１と蓄電池１１２間を接続し、電力系統４０１側の交流電力と蓄電池１１２側の直流電力とを、双方向に変換する機能を有する。交直変換機１１３は、単一のＡＣ／ＤＣコンバータを含む構成でも良いし、ＡＣ／ＤＣコンバータ、ＤＣ／ＤＣコンバータ、ＡＣ／ＡＣコンバータのうちの２種類以上のコンバータを任意に接続した構成でもよい。交直変換機１１３は、制御部１１１から提供される充放電命令に応じて、蓄電池１１２に放電指示又は充電指示を出力する。これにより、蓄電池１１２の充放電を行う。
【００１５】
蓄電池１１２は、電気エネルギーを蓄積（充電）および放電可能な蓄電池である。蓄電池１１２に蓄積されている電力量を充電量又は電荷量と呼ぶ。
【００１６】
蓄電池１１２は、１つ以上の電池盤を備え、各電池盤は、一例として、１つ以上の電池モジュールと、１つのＢＭＵ（Ｂａｔｔｅｒｙ  Ｍａｎａｇｅｍｅｎｔ  Ｕｎｉｔ：電池監視部）とを備える。各電池モジュールは、複数の単位電池（セル）を備える。
【００１７】
蓄電池１１２は、交直変換機１１３から放電指示を受信した場合、セル群に蓄電している電力を交直変換機１１３に放電する。蓄電池１１２は、交直変換機１１３から充電指示を受信した場合、電力系統４０１から供給される電力をセル群に充電する。セル、電池モジュール、電池盤および蓄電池１１２は、それぞれいずれも、電気エネルギーを内部に蓄積する蓄電装置の一形態である。
【００１８】
以下、蓄電池１１２の構成例について詳細に説明する。
【００１９】
図２は、蓄電池１１２の構成例を示す図である。蓄電池１１２は、複数の電池盤１１が並列に接続された電池アレイを含む。各電池盤１１において、複数の電池モジュール１２が直列に接続されている。各電池盤１１は、ＢＭＵ１３を備えている。複数の電池モジュール１２が並列に接続されてもよいし、複数の電池モジュール１２が直並列に接続されてもよい。また、複数の電池盤が、直列または直並列に接続されていてもよい。ＢＭＵ１３は、蓄電池評価装置２０１と情報を送受信する通信部を備えていてもよい。当該通信部は、電池盤１１の内部に配置されていても、電池盤１１の外部に配置されていてもよい。
【００２０】
図３は、各電池モジュール１２の構成の一例を示す図である。電池モジュール１２は、直並列に接続された複数のセル１４を備える。複数のセル１４が直列に接続された構成、並列に接続された構成、および直列と並列を組み合わせた構成も可能である。電池モジュール１２がＣＭＵ（Ｃｅｌｌ  Ｍｏｎｉｔｏｒｉｎｇ  Ｕｎｉｔ：セル監視部）を備えていてもよい。各セル１４は、充放電が可能な単位電池である。例えば、リチウムイオン電池、リチウムイオンポリマー電池、鉛蓄電池、ニッケルカドミウム電池、ニッケル水素電池などが挙げられる。
【００２１】
各セル１４に対して、電圧、電流、温度等のパラメータを計測する計測部（図示せず）が配置されている。同様に、各電池モジュール１２に対して、電池モジュールの電圧、電流、温度等のパラメータを計測する計測部（図示せず）が配置されている。また、各電池盤１１に対して、電池盤の電圧、電流、温度等のパラメータを計測する計測部（図示せず）が配置されている。また電池アレイに対して、電池アレイの電圧、電流、温度等のパラメータを計測する計測部（図示せず）が配置されている。ここでは、セル、電池モジュール、電池盤、および電池アレイのすべてに対して、電圧、電流、温度等を計測する計測部が配置されているとしたが、これらのうちの一部の種類にのみ配置されていてもよい。また、すべてのセルではなく、診断対象となる一部のセルにのみ計測部が配置されてもよい。このことは、セル以外の蓄電装置（電池モジュール、電池盤）についても同様である。また、各計測部は、電圧、電流、温度以外に、湿度など他のパラメータを計測してもよい。なお、計測部のサンプリング頻度は、充放電指令値の受信頻度により高くても、同じでもよい。
【００２２】
蓄電池は、蓄電池に関する計測データを蓄電池評価装置２０１に送出する。計測データは、計測部により計測されたパラメータ（電圧、電流、温度等）および計測時刻を含む。計測時刻は、蓄電池１１２内又は蓄電システム１０１に時計を配置し、パラメータの取得時に当該時計の時刻を計測時刻として取得してもよい。なお、蓄電池評価装置２０１に時計を設け、蓄電池評価装置２０１が蓄電システム１０１から計測データを取得した時に、取得した計測データに、時計の時刻を関連づけてもよい。この場合、計測データに計測時刻を含めなくてもよい。また、本実施形態では、蓄電池１１２を評価対象とする例を示すが、ここでいう蓄電池１１２の評価とは、蓄電池１１２に含まれる電池セルすべての集合体を評価することであり、蓄電池に含まれる計測部や、セル監視部、コントローラを評価するのとは異なる。また、本実施形態は、セル単体から、セルが複数集まった集合が階層構造まで、計測値がある任意の階層に対して評価を行うことができる。本実施形態に係る蓄電装置の評価は、セル、電池モジュール（実際には電池モジュールに含まれるセルの集合体）、電池盤（実際には電池盤に含まれるセルの集合体）など、計測値のある電池セル集合単位を評価することに相当する。
【００２３】
蓄電池評価装置２０１は、データ取得部２１１、データ記憶部２１２、時間区間選出部２１３、効率算出部２１４、クラスタリング部２１５、出力部２１６及び表示部２１７を備える。
【００２４】
データ取得部２１１は、蓄電システム１０１から蓄電池の計測データを取得する。また、データ取得部２１１は、蓄電システム１０１の制御部１１１から充放電指令値（充放電の実行時刻が付加されているものとする）取得する。ここでいう充放電指令値は、制御部１１１がＳＣＡＤＡの充放電指令値と蓄電池１１２の充電量に基づき生成した充放電指令値（前述した充放電命令）に対応する。データ取得部２１１は、計測データ及び充放電指令値の取得を一定時間毎に行ってもよい。データ取得部２１１は、取得した計測データ及び充放電指令値をデータ記憶部２１２に格納する。
【００２５】
データ記憶部２１２は、データ取得部２１１により取得された計測データ及び充放電指令値を含むデータを充放電履歴データとして内部に格納する。
【００２６】
図４に、データ記憶部２１２に記憶された充放電履歴データの例を示す。時刻ｔ１〜ｔｎまで、充放電指令値と、蓄電池の温度及び蓄電池の充電量が格納されている。時刻ｔ１〜ｔｎはサンプル時刻に対応し、本実施形態では充放電指令値の実行時刻及びパラメータの計測時刻は一致するとする。
【００２７】
Ｐ１,Ｐ２、・・・,Ｐｎは充放電指令値である。Ｐ１、Ｐ２、・・・、Ｐｎは符号付の数値である。一例として、正値は放電、負値は充電の指示を表すとするが、逆でもよい。
【００２８】
Ｑ１,Ｑ２、・・・Ｑｎは蓄電池１１２に蓄積されている電荷量（充電量）を表す。蓄電池１１２を放電すると、蓄電池１１２の充電量は減り、蓄電池１１２を充電すると蓄電池の充電量は増える。電力量を定格容量、すなわち劣化前の最大の充電量で除した割合をＳＯＣ（State of Charge）と呼ぶことにする。充電量をＳＯＣで表してもよい。
【００２９】
Ｔ１,Ｔ２、・・・Ｔｎは、蓄電池１１２の温度を表す。
【００３０】
例えば、時刻ｔ２のデータは、充放電指令値Ｐ２、充電量Ｑ２及び温度Ｔ２を含む。これは、充放電指令値Ｐ２の充放電が時刻ｔ２で実行され、この結果、蓄電池１１２の充電量はＱ２になり、またこのときの蓄電池１１２の温度がＴ２であることを意味する。他の時刻のデータも同様にして解釈される。
【００３１】
時間区間選出部２１３は、データ記憶部２１２に格納されている充放電履歴データに基づき、充電量が閾値条件を満たす時刻を探索する。探索は、一例として、予め与えられた探索開始時刻から探索終了時刻までの時間範囲で行う。閾値条件は、一例として、充電量が閾値θに達することである。充電量が閾値θに達するとは、充電量が閾値θより大きい値から閾値θになる場合、充電量が閾値θより低い値から閾値θになる場合のいずれも含む。その他の例として、閾値を基準として一定の範囲（例えばθ−Ｌ以上かつθ＋Ｍ以下の範囲）に充電量が入ること、つまり、一定の範囲外から一定の範囲内に充電量が入ることでもよい。閾値に基づく限り、ここで述べた以外の条件でもよい。
【００３２】
以下では、閾値条件を、充電量が閾値θに達する場合を想定して説明を行う。この場合、充電量が閾値条件を満たす時刻は、一例として、充電量が閾値θになった時刻である。充電量が閾値θになった時刻がサンプル時刻に一致する場合は、サンプル時刻が、充電量が閾値θになった時刻である。充電量が閾値θになった時刻がサンプル時刻間の場合は、線形補間により閾値θになった時刻を算出することも可能である。または、充電量が閾値θになった直後（又は直前）のサンプル時刻を、充電量が閾値θになった時刻としてもよい。
【００３３】
以下、図５及び図６を用いて、線形補間により充電量が閾値θになった時刻を特定する例を示す。
【００３４】
図５に、時刻ｔ（ｉ）の充電量が閾値θよりも小さく、その次の時刻ｔ（ｉ＋１）での充電量が閾値θよりも大きい場合の充放電履歴データの例を示す。充電量が閾値θになった時刻は、
ｔ＝（１−α）ｔ（ｉ）＋αｔ（ｉ＋１）として算出できる。
また、α＝（θ−Ｑ（ｉ））／（Ｑ（ｉ＋１）−Ｑ（ｉ））である。
【００３５】
図６に、横軸を時刻、縦軸を充電量とした座標系に、時刻ｔ（ｉ）とｔ（ｉ＋１）のデータをマッピングし、マッピングされた座標を通過する直線のグラフを記載した例を示す。この直線のグラフと、閾値θの直線とが交わる点の時刻が、充電量が閾値θになった時刻である。
【００３６】
充電量が閾値θになった時刻における充放電指令値についても、線形補間により、
（１−α）Ｐ（ｉ）＋αＰ（ｉ＋１）として算出できる。
【００３７】
同様に、充電量が閾値θになった時刻における温度についても、線形補間により、
（１−α）Ｔ（ｉ）＋αＴ（ｉ＋１）として算出できる。
【００３８】
充電量が閾値θになった直後のサンプル時刻を、充電量が閾値θになった時刻とする場合、図５の例では、時刻ｔ（ｉ＋１）が、充電量が閾値θになった時刻である。また、このときの充放電指令値及び温度は、Ｐ（ｉ＋１）及びＴ（ｉ＋１）である。
【００３９】
閾値θの値は、予め与えられていてもよいし、予め与えられた閾値の値を調整したものでもよい。また、予め与えられた複数の閾値候補の中から選択することで、使用する閾値θを決定してもよい。
【００４０】
時間区間選出部２１３は、探索により検出した時刻群を時間順に並べる。時間区間選出部２１３は、抽出条件を満たす時刻間の区間を、時間区間として選出する。時間区間は互いに重複しない（同じ区間を共有しない）とするが、重複を許容してもよい。
【００４１】
ここで、抽出条件について具体的に説明する。第１時刻と第２時刻で、充電量Ｑが閾値条件を満たす時刻として特定されているとする。
【００４２】
（第１の例）第１時刻と第２時刻間の時間差が、一定時間Δｔa以上である場合に、抽出条件が満たされる。この場合、第１時刻と第２時刻間の区間を、時間区間として選出する。時間差が一定時間Δｔa未満の場合は、抽出条件は満たされない。この場合、当該時間区間を選出しない。時間差としてデータサンプル数を用いてもよい。この場合、データサンプル数が所定値ｎ以上である場合に、抽出条件が満たされ、所定値ｎ未満の場合に、抽出条件が満たされない。
【００４３】
（第２の例）第１時刻と第２時刻間の充放電指令値が０である時間帯の長さが一定時間Δｔb以下である場合に抽出条件が満たされる。充放電指令値が０の場合としては、ＳＣＡＤＡ３０１からの充放電指令値が０であることにより制御部１１１が生成する充放電指令値（充放電命令値）も０であることがある。また、ＳＣＡＤＡ３０１からの充放電指令値が０ではないが、蓄電池１１２が満充電状態でこれ以上充電できない場合の充電指令、あるいは蓄電池１１２が完放電状態でこれ以上放電できない場合の放電指令に対して、制御部１１１が０の充放電指令値を生成する場合もある。このように充放電指令値が０である時間帯の長さが一定時間Δｔｂ以下の場合に、第１時刻と第２時刻間の区間を、時間区間として選出する。一定時間Δｔｂ以上の場合は、抽出条件は満たされない。この場合、当該時間区間を選出しない。
【００４４】
（第３の例）第１時刻と第２時刻間の充電量の平均値μＱと、閾値θとの差が所定値ｈ以下である場合に抽出条件が満たされる。この場合、第１時刻と第２時刻間の区間を時間区間として選出する。所定値ｈ未満の場合は、抽出条件は満たされない。この場合、当該時間区間を選出しない。
【００４５】
（第４の例）
第１〜第３の例の抽出条件を２つ以上組み合わせ、組み合わせた抽出条件が全て満たされた場合にのみ、第１時刻と第２時刻間の区間を時間区間として選出する。組み合わせた抽出条件のうちの少なくとも１つが満たされない場合は、第１時刻と第２時刻間の区間を選出しない。
【００４６】
効率算出部２１４は、時間区間選出部２１３が選出した各時間区間について、放電指令値の総和と、充電指令値の総和とを計算する。効率算出部２１４は、放電指令値の総和と、充電指令値の総和との比により、充放電効率を評価指標として計算する。
【００４７】
図７を用いて、充放電効率を計算する具体例を示す。図７は、ある充放電履歴データの一部を示す。充放電指令値の符号から、充放電指令値Ｐ１，Ｐ２，Ｐ４，Ｐ５，Ｐ６が放電指示であり、Ｐ３，Ｐ７，Ｐ８が充電指示である。充電量が閾値条件を満たした時刻としてｔ２，ｔ５，ｔ７が特定されている。これらのうち時刻ｔ２とｔ７間の区間が抽出条件を満たすため、時間区間Ａとして選出されている。時刻ｔ５とｔ７間の区間は抽出条件を満たさないため、時間区間として選出されていない。
【００４８】
この場合、時間区間Ａにおいて放電指令値により指示された放電量の総和は、
｜Ｐ２｜（ｔ２−ｔ１）＋｜Ｐ４｜（ｔ４−ｔ３）＋｜Ｐ５｜（ｔ５−ｔ４）
である。
時間区間Ａにおいて充電指令値により指示された充電量の総和は、
｜Ｐ３｜（ｔ３−ｔ２）＋｜Ｐ６｜（ｔ６−ｔ５）＋｜Ｐ７｜（ｔ７−ｔ６）である。
よって、時間区間Ａの充放電効率ｒＡは、以下の式で算出される。
充放電効率ｒＡ＝
（｜Ｐ２｜（ｔ２−ｔ１）＋｜Ｐ４｜（ｔ４−ｔ３）＋｜Ｐ５｜（ｔ５−ｔ４））／（｜Ｐ３｜（ｔ３−ｔ２）＋｜Ｐ６｜（ｔ６−ｔ５）＋｜Ｐ７｜（ｔ７−ｔ６））
【００４９】
充放電効率ｒＡが大きいほど、蓄電池の使用効率が高いことを意味する。
【００５０】
図８に、充放電履歴データの他の例を示す。時刻ｔ１〜ｔ８までのデータは、図７と同じある。図８の充放電履歴データでは、図７のデータに加えて、時刻ｔ９〜ｔ２２のデータが含まれる。ｔ１３〜ｔ１５では充放電が指示されていない（充放電指令値が０）。充電量が閾値条件を満たす時刻として、時刻ｔ２，ｔ５，ｔ７に加えて、時刻ｔ１１，ｔ１６，ｔ２１が検出されている。これらの時刻に基づき、時間区間Ａ（ｔ２〜ｔ７）、時間区間Ｂ（ｔ７〜ｔ１１）、時間区間Ｃ（ｔ１６〜ｔ２１）が、抽出条件を満たすとして選出されている。
【００５１】
時間区間Ａについては、前述した計算方法で充放電効率ｒＡを算出する。時間区間Ｂ及び時間区間Ｃについても時間区間Ａと同様の計算方法で、充放電効率ｒＢ及び充放電効率ｒＣを算出する。
【００５２】
図９に、図８の充放電履歴データを、横軸を時刻、縦軸を充電量とする座標系に展開してグラフを表示した例を示す。グラフが閾値θと交わる時刻が、充電量が閾値条件を満たした時刻に対応する。図８の充放電履歴データに対して特定された時間区間Ａ、Ｂ、Ｃが記載されている。
【００５３】
出力部２１６は、効率算出部２１４から得られる充放電効率を、時間区間を代表する時刻（代表時刻）と関連づけたデータを生成する。代表時刻は、一例として、充放電効率を算出する元となった時間区間の終了時刻である。当該データは、例えば（代表時刻、充放電効率）として表される。出力部２１６は、当該データを、時刻と充放電効率を軸とする座標系にプロットする。表示部２１７は、座標系にプロットされたデータの集合、すなわち充放電効率の推移データを画面に表示する。なお、当該データを出力部２１６が生成する代わりに、後述するクラスタリング部２１５が生成して、生成したデータを後述のクラスタ識別情報又はクラスタ特徴情報を付加して出力部２１６に提供してもよい。
【００５４】
図８の充放電履歴データの場合、出力部２１６は、時間区間Ａ、Ｂ、Ｃに対して計算された充放電効率ｒＡ，ｒＢ，ｒＣを、それぞれの代表時刻ｔＡ，ｔＢ，ｔＣと関連づけることにより、データ（ｔＡ，ｒＡ），（ｔＢ，ｒＢ），（ｔＣ，ｒＣ）を生成する。出力部２１６は、これらのデータを、時刻及び充放電効率を軸とする座標系にプロットする。
【００５５】
図１０に、これらのデータを、時刻及び充放電効率を軸とする座標系にプロットした例を示す。ｔＡ，ｔＢ，ｔＣは、充放電効率ｒＡ，ｒＢ，ｒＣを算出した時間区間の終了時刻である。すなわち、ｔＡ，ｔＢ，ｔＣは、それぞれｔ７，ｔ１１，ｔ２１である。
【００５６】
充放電効率を代表時刻ではなく、蓄電池１１２を使用し始めてから当該時間区間の終了時刻までの累積放電量又は累積充電量と関連付けることにより、データを生成してもよい。この場合、当該データを、充放電効率と、蓄電池の使用開始時からの累積放電量又は累積充電量をそれぞれ軸とする座標系にプロットすればよい。
【００５７】
クラスタリング部２１５は、上記の座標系にプロットされたデータの集合をクラスタリングして、データの集合を複数のクラスタに分割する。クラスタリングの手法は、一例として、一般的に知られた階層的方法、又は凝縮法を用いることができる。階層的方法は、点群をトップダウンに分割していく方法、凝縮法は、予め任意の方法で細分化したクラスタを統合していく方法である。あるいは、特開２００６−３５０７３０号公報に開示されたクラスタリング方法を用いることもできる。この公報に記載のクラスタリング方法では人間の直感に近いクラスタリングを可能としている。本実施形態で用いるクラスタリングの手法は特定のものに限らず、任意の手法を用いることができる。
【００５８】
図１１は、クラスタリングの例を模式的に示す図である。この例では横軸は時刻であるが、横軸が累積放電量又は累積充電量の場合も可能である。図示の座標系に、１３個のデータがプロットされている。クラスタリング部２１５がこのデータの集合に対してクラスタリングを行う。この結果、２つのクラスタＣ１，Ｃ２が生成されている。クラスタリング部２１５は、クラスタＣ１，Ｃ２毎に、クラスタを近似する関数（近似関数）を計算する。各クラスタについて計算された近似関数のグラフが表示されている。この例では、１次関数（直線）によりクラスタを近似しているが、より高次の関数によりクラスタを近似してもよい。この場合は、近似関数は曲線（曲線のグラフ）になり得る。各クラスタに属するデータ、又は各クラスタを近似する関数のグラフは、各クラスタの充放電効率の推移データに対応する。クラスタＣ１，Ｃ２を近似する一次関数の傾きはそれぞれａ１，ａ２である。これらの傾きの情報が表示部２１７に表示されている。また、クラスタＣ１，Ｃ２を識別する情報（クラスタ識別情報）として、各クラスタを囲む枠を表示してもよいし、クラスタＣ１，Ｃ２に属するデータの点の色又は背景色をクラスタごとに変えてもよい。
【００５９】
クラスタリング部２１５は、各クラスタの特徴情報（クラスタ特徴情報）を生成する。クラスタ特徴情報は、一例として、クラスタに属する各点（各データ）の情報から算出される。クラスタリング部２１５は、生成したクラスタ特徴情報を、プロットしたデータの集合と同じ画面又は別の画面に表示する。各クラスタには、各クラスタ内で共通した蓄電池の使用条件の下で取得されたデータ（例えば電池温度、充放電指令値又はＳＯＣなどの要因が類似しているデータ）が属しているといえる。つまり、クラスタ特徴情報は、当該クラスタに属しているデータが得られたときに共通している蓄電池１１２の使用条件を代表的に表しているといえる。前述した傾きの情報及びクラスタを識別する情報（クラスタ識別情報）は、クラスタ特徴情報の一例であるが、その他のクラスタ特徴情報の生成例を示す。
【００６０】
第１の例として、クラスタに属する各データについて、充放電指令値の標準偏差を計算する。つまり各データが算出される元となった時間区間に含まれる充放電指令値の標準偏差を計算する。そして、これらの標準偏差の統計情報（平均、中央値、最大値、最小値等）を、当該クラスタのクラスタ特徴情報とする。クラスタＣ１のクラスタ特徴情報として、充放電指令値の標準偏差の平均値、クラスタＣ２のクラスタ特徴情報として、充放電指令値の標準偏差の平均値がある。充放電指令値の標準偏差の代わりに、充放電指令値の平均値又は分散等、他の値を計算し、計算した値の統計情報を、クラスタ特徴情報としてもよい。図１１ではクラスタ数は２であるが、クラスタ数は３以上でもよい。
【００６１】
第２の例として、クラスタに属する各データについて、各データが算出される元となった時間区間に含まれる温度の標準偏差を計算する。そして、これらの標準偏差の統計情報（平均、中央値、最大値、最小値等）を、当該クラスタのクラスタ特徴情報とする。温度の標準偏差の代わりに、温度の平均値又は分散等、他の値を計算し、計算した値の統計情報を、クラスタ特徴情報としてもよい。
【００６２】
第３の例として、クラスタに属する各データについて、各データが算出される元となった時間区間に含まれる充電量（電荷量）の標準偏差を計算する。そして、これらの標準偏差の統計情報（平均、中央値、最大値、最小値等）を、当該クラスタのクラスタ特徴情報とする。充電量の標準偏差の代わりに、充電量の平均値又は分散等、他の値を計算し、計算した値の統計情報を、クラスタ特徴情報としてもよい。
【００６３】
第４の例として、クラスタに属する各データについて、各データが算出される元となった時間区間の区間長（時間長）の標準偏差を計算する。そして、これらの標準偏差の統計情報（平均、中央値、最大値、最小値等）を、当該クラスタのクラスタ特徴情報とする。区間長の標準偏差の代わりに、区間長の平均値又は分散等、他の値を計算し、計算した値の統計情報を、クラスタ特徴情報としてもよい。
【００６４】
上に記載したクラスタ特徴情報は一例に過ぎず、充放電指示値、充電量、温度、時刻等に基づく情報であれば、他の情報でもよい。
【００６５】
蓄電池評価装置２０１のオペレータは、表示部２１７に表示された画面を参照することで、充放電効率の良い蓄電池の使用条件を把握できる。図１１の例では、クラスタＣ１の充放電効率は、クラスタＣ２に比べて高いことが分かる。またクラスタＣ１の近似直線の傾きａ１が、クラスタＣ２の近似直線の傾きａ２よりも緩やかであることから、クラスタＣ１の方が劣化の速度も緩やかであると判断できる。そこで、オペレータは、蓄電池１１２の使用条件として、クラスタ特徴情報又はこれに近似した条件に従って充放電を行うことで、充放電の効率を高めることができ、かつ蓄電池１１２の劣化の進行も抑制できると判断できる。
【００６６】
このような蓄電池の使用条件を特定する情報を、オペレータの指示により蓄電池評価装置２０１又は別のユーザデバイスから電力系統４０１側のサーバに送信してもよい。サーバ側ではこの情報に基づき、充放電効率が高くなるように充放電指令値を生成してもよい。一例として、蓄電池評価装置２０１からサーバにクラスタＣ１のクラスタ特徴情報を送信し、サーバ側では、クラスタ特徴情報が満たされるように充放電指令値を生成してもよい。
【００６７】
図１２は、図１０又は図１１の例よりも長い測定期間に基づくデータをプロットした例を示す。横軸は時刻、縦軸は充放電効率である。
【００６８】
図１３は、図１０又は図１１の例よりも長い測定期間に基づくデータをプロットした例を示す。図１１と異なり、横軸は累積放電量、縦軸は充放電効率である。
【００６９】
図１２又は図１３では、図１０又は図１１よりも多数のデータがプロットされている。いずれの場合も、図１１を用いて説明したのと同様の方法によってデータのクラスタリングを行うことにより、複数のクラスタを生成する。そして、生成された各クラスタについて、クラスタを近似する関数の計算、及び、クラスタ特徴情報の生成を行う。
【００７０】
上述した実施形態では閾値θは１つであったが、閾値θの個数は、複数でもよい。複数の閾値θを用いた場合、閾値θごとに、データのプロット及びクラスタリング等を行う。この場合、クラスタ特徴情報の一部として、閾値θの値を含めて、表示部２１７に表示してもよい。これにより、オペレータは、閾値θと充放電効率との関係の傾向も把握することができる。
【００７１】
図１４に、閾値θとして２つの閾値θ１とθ２を使った場合に表示される画面の例を模式的に示す。閾値θ１を用いた場合にクラスタＣ１，Ｃ２が生成され、閾値θ２を用いた場合にクラスタＣ３，Ｃ４が生成されている。なお、図１４では近似直線の傾き等、その他の情報の表示は省略している。
【００７２】
図１５は、本実施形態に係る蓄電池評価装置２０１の動作の一例のフローチャートである。
【００７３】
ステップＳ１において、データ取得部２１１が充放電指令値及び計測データ（電池温度、温、計測時刻等）を取得し、取得し充放電指令値及び計測データを充放電履歴データとしてデータ記憶部２１２に格納する。
【００７４】
ステップＳ２において、時間区間選出部２１３が、充放電履歴データにおいて閾値条件を満たす時刻を検出する。時間区間選出部２１３が、検出した時刻に基づき、抽出条件を満たす時刻間の区間を時間区間として選出する。
【００７５】
ステップＳ３において、効率算出部２１４は、選出した時間区間に対して充放電効率を計算する。
【００７６】
ステップＳ４において、クラスタリング部２１５は、充放電効率と、時間区間の代表時刻（あるいは累積放電量、累積充電量）とを含むデータをクラスタリングして、複数のクラスタを生成する。クラスタリング部２１５は、各クラスタを近似する関数の計算、及び、クラスタ特徴情報の生成を行う。
【００７７】
ステップＳ５において、出力部２１６は、クラスタリングされたデータを座標系にマッピングし、マッピングされたデータを表示部２１７に出力する。また、クラスタ特徴情報、近似関数のグラフ、各クラスタを近似する直線の傾きの情報等も、表示部２１７に出力する。表示部２１７は、出力部２１６から入力された各種の情報を画面に表示する。なお、座標系へのデータのマッピングと、データのクラスタリングの順序はどちらが先に行われてもよい。
【００７８】
本実施形態によれば、通常稼働している蓄電システムから受信した計測データ及び充放電指令値に基づき、充電量が閾値条件を満たす時刻を特定し、抽出条件を満たす時刻間を時間区間として選出する。選出した時間区間ごとに充放電効率を評価指標として算出し、充放電効率をクラスタリングする。各クラスタの特徴情報を比較することで、充放電効率の良い蓄電池の使用条件（例えば電池温度、充放電指令値又はＳＯＣなどの条件）を評価できる。また、通常稼働している蓄電システムから得られる計測データを利用するため、蓄電システムの評価のために、蓄電システムをオフラインにする必要はなく、特殊な充放電パターンの充放電を実行させる必要もない。充放電効率の変化を時間、累積放電量又は累積充電量を軸とする座標系にプロットすることで、蓄電池の健全性の時間変化を可視化できる。
【００７９】
（変形例１）
充放電効率を、時間区間を代表する温度（代表温度）により補正してもよい。代表温度は、例えば、時間区間の平均温度、最高温度、最低温度等である。一般に蓄電池の内部抵抗は、電池温度が高くなると低くなる傾向がある。このため、稼働時の電池温度が高いと、温度が低い場合に比べて、算出される充放電効率が高くなることがある。そこで、予めある基準温度を決めておき、基準温度より低い場合は、充放電効率を高くなるように充放電効率を補正し、基準温度より高い場合は、充放電効率が低くなるように充放電効率を補正する。このために、事前に、温度と充放電効率との関係を試験又はシミュレーション等により取得しておく。そして、この関係を利用して、基準温度の充放電効率を基準として、充放電効率を正規化（補正）する。例えば、温度及び充放電効率を入力変数とし、補正された温度を出力変数とする式を作成しておく。この式に代表温度及び充放電効率を入力することで、充放電効率を補正する。
【００８０】
図１６に、図１２の各点の充放電効率を代表温度に応じて補正し、補正した充放電効率をマッピングした例を示す。また、図１７に、図１３の各点の充放電効率を代表温度に応じて補正し、補正した充放電効率をマッピングした例を示す。
【００８１】
充放電効率を補正した後の処理は、上述した実施形態と同様である。すなわち、補正した充放電指令値を代表時刻と関連づけたデータを生成し、データを座標系にマッピング及びクラスタリング等する。またクラスタを近似する関数を計算する。また、クラスタ特徴情報の生成を行う。
【００８２】
（変形例２）
閾値θを決定する方法を示す。一例として、閾値θは、蓄電池の定格容量の所定割合の値とする。例えば、蓄電池の定格容量の１／２を閾値θとする。この場合、充電量をＳＯＣで表すと、閾値θは５０％となる。
【００８３】
他の方法として、充電量の分布を作成し、作成した分布に基づき、閾値θを決定してもよい。
【００８４】
図１８は、探索開始時刻から探索終了時刻までの間の充電量Ｑのヒストグラムを模式的に示す。仮に閾値がヒストグラムの上限値より大きいθｚである場合、充電量がこの閾値θｚに達することはないから、閾値条件を満たす時刻を検出できない。従って、時間区間の選出もできない。そこで、ヒストグラムの上限と下限の間に閾値を設定する。例えば、図示のθｘ又はθｙ又はこれらの両方を設定する。ヒストグラムの中心又はその近くに閾値を設定することで、多くの時間区間を選出することが期待される。これにより、より多くの時間区間について充放電効率を計算でき、クラスタリングのためのサンプル数を増やすことができる。複数の閾値を設定した場合、閾値ごとに、閾値条件を満たす時刻の検出、時間区間の選出、充放電効率の計算、マッピング、クラスタリング、クラスタ特徴情報の生成等を行う。
【００８５】
図１９は、本発明の実施形態に係る蓄電池評価装置のハードウェア構成例を示す。図１９のハードウェア構成はコンピュータ１５０として構成される。コンピュータ１５０は、ＣＰＵ１５１、入力インタフェース１５２、表示装置１５３、通信装置１５４、主記憶装置１５５、外部記憶装置１５６を備え、これらはバス１５７により相互に通信可能に接続される。
【００８６】
入力インタフェース１５２は、蓄電池で計測された計測データを配線等を介して取得する。通信装置１５４は、無線または有線の通信手段を含み、蓄電システム１０１と有線または無線の通信を行う。通信装置１５４を介して、計測データを取得してもよい。入力インタフェース１５２及び通信装置１５４は、それぞれ別個の集積回路等の回路で構成されていてもよいし、単一の集積回路等の回路で構成されてもよい。
【００８７】
外部記憶装置１５６は、例えば、ＨＤＤ、ＳＳＤ、メモリ装置、ＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＡＭ、ＤＶＤ−Ｒ等の記憶媒体等を含む。外部記憶装置１５６は、蓄電池評価装置の各処理部の機能を、プロセッサであるＣＰＵ１５１に実行させるためのプログラムを記憶している。また、データ記憶部２１２も、外部記憶装置１５６に含まれる。ここでは、外部記憶装置１５６を１つのみ示しているが、複数存在しても構わない。
【００８８】
主記憶装置１５５は、ＣＰＵ１５１による制御の下で、外部記憶装置１５６に記憶された制御プログラムを展開し、当該プログラムの実行時に必要なデータ、当該プログラムの実行により生じたデータ等を記憶する。主記憶装置１５５は、例えば揮発性メモリ（ＤＲＡＭ、ＳＲＡＭ等）または不揮発性メモリ（ＮＡＮＤフラッシュメモリ、ＭＲＡＭ等）など、任意のメモリまたは記憶部を含む。主記憶装置１５５に展開された制御プログラムがＣＰＵ１５１により実行されることで、蓄電池評価装置２０１の各処理部の機能が実行される。
【００８９】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００９０】
１０１：蓄電システム
１１１：制御部
１１２：蓄電池
１１３：交直変換機
２０１：蓄電池評価装置
２１１：データ取得部
２１２：データ記憶部
２１３：時間区間選出部
２１４：効率算出部
２１５：クラスタリング部
２１６：出力部
２１７：表示部
３０１：ＳＣＡＤＡ
４０１：電力系統
１５１：ＣＰＵ
１５２：入力インタフェース
１５３：表示装置
１５４：通信装置
１５５：主記憶装置
１５６：外部記憶装置
１５７：バス
【図１】
2020109720.tif 000003
【図２】
2020109720.tif 000004
【図３】
2020109720.tif 000005
【図４】
2020109720.tif 000006
【図５】
2020109720.tif 000007
【図６】
2020109720.tif 000008
【図７】
2020109720.tif 000009
【図８】
2020109720.tif 000010
【図９】
2020109720.tif 000011
【図１０】
2020109720.tif 000012
【図１１】
2020109720.tif 000013
【図１２】
2020109720.tif 000014
【図１３】
2020109720.tif 000015
【図１４】
2020109720.tif 000016
【図１５】
2020109720.tif 000017
【図１６】
2020109720.tif 000018
【図１７】
2020109720.tif 000019
【図１８】
2020109720.tif 000020
【図１９】
2020109720.tif 000021