(57)【特許請求の範囲】
【請求項１】
現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定する判定部と、
前記所定の範囲から外れていると判定された場合に、現在時刻が属する日の後日における同一時刻に対応付けられた消費電力のデータを含む前記データ分布に基づいて、現在時刻に対応付けられた消費電力のデータが、前記所定の範囲から外れていると判定した回数を算出する外れ値回数算出部と、
算出した回数に基づいて、前記所定の範囲から外れているとの判定結果が継続する度合いを示す状態変化度を算出する状態変化度算出部と
を有することを特徴とするデータ処理装置。
【請求項２】
前記状態変化度を、前記消費電力のデータと時刻とが対応付けられた時系列データに対応付けて通知する通知部を更に有することを特徴とする請求項１に記載のデータ処理装置。
【請求項３】
前記通知部は、
前記時系列データの各時刻に対応する位置に配された直線領域において、前記状態変化度に応じて動作する識別子をプロットすることで、前記状態変化度を通知することを特徴とする請求項２に記載のデータ処理装置。
【請求項４】
前記通知部は、更に、
前記判定部により所定の範囲から外れていると判定されたことを通知することを特徴とする請求項２に記載のデータ処理装置。
【請求項５】
前記判定部は、更に、
現在時刻に対応付けられた消費電力の所定時間前からの変動率のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力の所定時間前からの変動率のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力の所定時間前からの変動率のデータが、所定の範囲から外れているか否かを判定し、
前記通知部は、
前記消費電力のデータが所定の範囲から外れていると判定されたこと、または、前記変動率のデータが所定の範囲から外れていると判定されたことを通知することを特徴とする請求項４に記載のデータ処理装置。
【請求項６】
前記判定部は、
前記データ分布において、統計的外れ値に該当するデータであるか否かを判定することで、前記所定の範囲から外れているか否かを判定することを特徴とする請求項１に記載のデータ処理装置。
【請求項７】
コンピュータに、
現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定し、
前記所定の範囲から外れていると判定された場合に、現在時刻が属する日の後日における同一時刻に対応付けられた消費電力のデータを含む前記データ分布に基づいて、現在時刻に対応付けられた消費電力のデータが、前記所定の範囲から外れていると判定した回数を算出し、
算出した回数に基づいて、前記所定の範囲から外れているとの判定結果が継続する度合いを示す状態変化度を算出する、
処理を実行させることを特徴とするデータ処理プログラム。
【請求項８】
現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定し、
前記所定の範囲から外れていると判定された場合に、現在時刻が属する日の後日における同一時刻に対応付けられた消費電力のデータを含む前記データ分布に基づいて、現在時刻に対応付けられた消費電力のデータが、前記所定の範囲から外れていると判定した回数を算出し、
算出した回数に基づいて、前記所定の範囲から外れているとの判定結果が継続する度合いを示す状態変化度を算出する、
処理をコンピュータが実行することを特徴とするデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、消費電力の時系列データを処理するデータ処理装置、データ処理プログラム及びデータ処理方法に関する。
【背景技術】
【０００２】
従来より、電力需要の日内変動を予測する予測装置が知られている。予測装置は、電力需要の日内変動を表す変動パターンを複数有しており、早朝時間帯における電力需要の時系列データを取得して各変動パターンとの比較を行い、類似する変動パターンを抽出することで、昼間の時間帯の電力需要の変動予測を行っている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−１１４６２９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、電力需要の日内変動は、例えば季節ごとに傾向があり、季節の変わり目等に、それまでとは異なる変動パターンが出現しはじめること（「状態変化」ともいう）がある。
【０００５】
従来の予測装置では、このような状態変化を、一定期間が経過した後に検出していたため、この期間の電力需要の変動予測において、予測精度を維持することが困難であった。一方で、状態変化の兆しを速やかに検出することができれば、予測精度を維持するための対策等を行うこともできる。
【０００６】
一つの側面では、状態変化の兆しを速やかに検出することを目的とする。
【課題を解決するための手段】
【０００７】
一態様によれば、データ処理装置は、現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定する判定部と、前記所定の範囲から外れていると判定された場合に、現在時刻が属する日の後日における同一時刻に対応付けられた消費電力のデータを含む前記データ分布に基づいて、現在時刻に対応付けられた消費電力のデータが、前記所定の範囲から外れていると判定した回数を算出する外れ値回数算出部と、算出した回数に基づいて、前記所定の範囲から外れているとの判定結果が継続する度合いを示す状態変化度を算出する状態変化度算出部とを有することを特徴とする。
【発明の効果】
【０００８】
状態変化の兆しを速やかに検出することが可能となる。
【図面の簡単な説明】
【０００９】
【図１】消費電力監視装置を備える監視システムの全体構成の一例を示す図である。
【図２】消費電力監視装置のハードウェア構成を示す図である。
【図３】消費電力の時系列データの一例を示す図である。
【図４】判定用データ生成部により生成される判定用データの一例を示す図である。
【図５】判定用データ生成部による判定用データ生成処理のフローチャートである。
【図６】判定用データ生成部により所定周期ごとに判定用データが生成される様子を示す図である。
【図７】消費電力判定部の詳細な機能構成の一例を示す図である。
【図８】変化発生有無判定部による変化発生有無判定処理のフローチャートである。
【図９】変化発生有無判定処理の具体例を示す図である。
【図１０】状態変化度算出部による状態変化度算出処理のフローチャートである。
【図１１】状態変化度算出処理の第１の具体例を示す図である。
【図１２】状態変化度算出処理の第２の具体例を示す図である。
【図１３】状態変化度算出処理の第３の具体例を示す図である。
【図１４】判定部による判定処理のフローチャートである。
【図１５】判定部による判定処理のフローチャートである。
【図１６】消費電力の変化が発生してからの経過日数に対する状態変化度の推移を示す図である。
【図１７】判定結果通知部により通知された判定結果の一例を示す図である。
【発明を実施するための形態】
【００１０】
以下、各実施形態について添付の図面を参照しながら説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複した説明を省く。
【００１１】
［第１の実施形態］
はじめに、本実施形態に係るデータ処理装置の一例である消費電力監視装置を備える監視システム１４０について説明する。図１は、消費電力監視装置を備える監視システムの全体構成の一例を示す図である。
【００１２】
図１に示すように、監視システム１４０は、消費電力監視装置１００と、表示装置１１０と、電力制御装置１２０とを備える。消費電力監視装置１００は表示装置１１０と通信可能に接続されている。また、消費電力監視装置１００は、電力制御装置１２０と通信可能に接続されている。
【００１３】
消費電力監視装置１００は、消費電力（例えば、１時間ごとの平均の消費電力（ｋＷ））と時刻とが対応付けられた"消費電力の時系列データ"を解析する。また、消費電力監視装置１００は、解析結果に基づいて季節の変わり目等における状態変化の兆しを検出するとともに、検出した兆しが状態変化の発生によるものであったか否か（変化の属性）を判定する。
【００１４】
なお、「状態変化の兆し」を検出するとは、季節の変わり目等において状態変化が発生する場合に、状態変化が発生したことを示す「消費電力の変化」を検出することをいう。ここで、「消費電力の変化」とは、消費電力の日内変動において、それまでの日内変動の変動パターンから外れた消費電力が出現することをいう。
【００１５】
消費電力監視装置１００には、時系列データ取得プログラム、判定用データ生成プログラム、消費電力判定プログラム、判定結果通知プログラムがインストールされている。消費電力監視装置１００は、これら各種プログラムが実行されることで、時系列データ取得部１０１、判定用データ生成部１０２、消費電力判定部１０３、判定結果通知部１０４として機能する。
【００１６】
時系列データ取得部１０１は、電力制御装置１２０から、消費電力の時系列データを順次取得する。消費電力の時系列データは、所定周期（例えば、１時間周期）で電力制御装置１２０から順次出力されるため、時系列データ取得部１０１では、所定周期で取得した時系列データを、複数日分、時系列データ格納部１１１に格納する。
【００１７】
判定用データ生成部１０２は、時系列データ格納部１１１から消費電力の時系列データを読み出す。また、判定用データ生成部１０２は、現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日（過去Ｎ日分）の消費電力のデータとを抽出することで、対象時刻データを生成し、対象時刻データ格納部１１２に格納する。
【００１８】
また、判定用データ生成部１０２は、生成した対象時刻データに基づいて、各日の対象時刻の所定時間前（例えば１時間前）の消費電力からの変動率を算出することで、第１の対象時刻変動率データを生成し、第１の対象時刻変動率データ格納部１１３に格納する。
【００１９】
また、判定用データ生成部１０２は、生成した対象時刻データに基づいて、各日の対象時刻の所定時間前（例えば２時間前）の消費電力からの変動率を算出することで、第２の対象時刻変動率データを生成し、第２の対象時刻変動率データ格納部１１４に格納する。
【００２０】
消費電力判定部１０３は、対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データそれぞれについて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行い、統計的外れ値の有無を判定する。なお、統計的外れ値の有無の判定とは、各データのデータ分布に対して、統計的処理を行った場合に所定の範囲から外れるデータの有無を判定することをいう。これにより、消費電力判定部１０３では、現在時刻において消費電力の変化が発生したか否かを判定する。
【００２１】
また、消費電力判定部１０３は、消費電力の変化が発生したと判定した場合に、当該消費電力の変化が、状態変化の発生によるものなのか否かを判定する。あるいは、状態変化の発生によるものではなく、そのとき限りの一時的な異常発生によるものなのか否かを判定する。つまり、消費電力判定部１０３は、消費電力の変化が発生したと判定した場合に、変化の属性についても判定する。
【００２２】
判定結果通知部１０４は、消費電力判定部１０３における判定結果を時系列データとともに、表示装置１１０に通知する。また、判定結果通知部１０４は、電力制御装置１２０にも判定結果を通知する。
【００２３】
表示装置１１０は、電力制御装置１２０を管理する電力管理担当者が、電力制御装置１２０により制御されている消費電力を確認するための装置である。表示装置１１０は、消費電力監視装置１００から通知された時系列データと判定結果とを表示する。これにより、電力管理担当者は、消費電力の変化が発生したことを認識することができる。また、消費電力の変化が、状態変化の発生によるものなのか、そのとき限りの一時的な異常発生によるものなのか、変化の属性を認識することができる。
【００２４】
電力制御装置１２０は、消費電力を制御する。また、制御中の消費電力を時刻と対応付けて、消費電力の時系列データとして、所定周期で消費電力監視装置１００に送信する。更に、電力制御装置１２０は、消費電力監視装置１００から通知された判定結果に基づいて、消費電力の日内変動の変動予測を行う。
【００２５】
次に、消費電力監視装置１００のハードウェア構成について説明する。図２は、消費電力監視装置１００のハードウェア構成を示す図である。図２に示すように、消費電力監視装置１００は、ＣＰＵ２０１、ＲＯＭ（Read Only Memory）２０２、ＲＡＭ（Random Access Memory）２０３を備える。また、消費電力監視装置１００は、補助記憶部２０４、操作部２０５、接続部２０６、通信部２０７、ドライブ部２０８を備える。なお、消費電力監視装置１００の各部は、バス２０９を介して相互に接続されている。
【００２６】
ＣＰＵ２０１は、補助記憶部２０４に格納された各種プログラム（例えば、時系列データ取得プログラム、判定用データ生成プログラム、消費電力判定プログラム、判定結果通知プログラム等）を実行するコンピュータである。
【００２７】
ＲＯＭ２０２は不揮発性メモリである。ＲＯＭ２０２は、補助記憶部２０４に格納された各種プログラムをＣＰＵ２０１が実行するために必要なプログラム、データ等を格納する主記憶部として機能する。具体的には、ＢＩＯＳ（Basic Input/Output System）やＥＦＩ（Extensible Firmware Interface）等のブートプログラム等を格納する主記憶部として機能する。
【００２８】
ＲＡＭ２０３は揮発性メモリであり、ＤＲＡＭ（Dynamic Random Access Memory）やＳＲＡＭ（Static Random Access Memory）等を含む。ＲＡＭ２０３は、補助記憶部２０４に格納された各種プログラムがＣＰＵ２０１によって実行される際に展開される、作業領域を提供する主記憶部として機能する。
【００２９】
補助記憶部２０４は、各種プログラムや、各種プログラムが実行されることで生成されるデータ等を格納するコンピュータ読み取り可能な媒体である。補助記憶部２０４の所定の領域は、例えば、時系列データ格納部１１１、対象時刻データ格納部１１２、第１の対象時刻変動率データ格納部１１３、第２の対象時刻変動率データ格納部１１４として機能する。
【００３０】
操作部２０５は、例えば、消費電力監視装置１００に対する各種操作を受け付ける。接続部２０６は、表示装置１１０と接続され、表示装置１１０に対して消費電力の時系列データ及び判定結果を通知する。
【００３１】
通信部２０７は電力制御装置１２０と通信し、電力制御装置１２０との間で、消費電力の時系列データ及び判定結果の送受信を行う。
【００３２】
ドライブ部２０８は記録媒体２１０をセットするためのデバイスである。ここでいう記録媒体２１０には、ＣＤ−ＲＯＭ、フレキシブルディスク、光磁気ディスク等のように情報を光学的、電気的あるいは磁気的に記録する媒体が含まれる。また、記録媒体２１０には、ＲＯＭ、フラッシュメモリ等のように情報を電気的に記録する半導体メモリ等が含まれていてもよい。
【００３３】
なお、本実施形態において、補助記憶部２０４に格納されている各種プログラムは、例えば、記録媒体２１０がドライブ部２０８にセットされ、該記録媒体２１０に記録された各種プログラムがドライブ部２０８により読み出されることで格納されるものとする。あるいは、通信部２０７を介してネットワークからダウンロードされることで格納されてもよい。
【００３４】
次に、電力制御装置１２０から送信され、時系列データ取得部１０１を介して時系列データ格納部１１１に格納される、消費電力の時系列データについて説明する。図３は消費電力を示す時系列データの一例を示す図である。
【００３５】
図３において、横軸は日付及び時刻を示しており、縦軸は消費電力を示している。なお、図３の例では、１時間ごとの平均の消費電力（ｋＷ）がプロットされている。
【００３６】
図３に示すように、消費電力は日単位で周期的に変動する。また、消費電力の日内変動は、日ごとに異なった変動パターンとなっており、それぞれの変動パターンを比較した場合に、差異が小さい変動パターンと大きい変動パターンとがある。例えば、図３の時点３０１、３０２等では、その前後で、消費電力の日内変動の変動パターンが大きく変化しており、差異が大きい。
【００３７】
次に、判定用データ生成部１０２により生成される判定用データ（対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データ）について説明する。図４は、判定用データ生成部により生成される判定用データの一例を示す図である。
【００３８】
このうち、図４（ａ）は、判定用データ生成部１０２により生成された対象時刻データの一例を示す図であり、横軸は日付及び時刻を示し、縦軸は消費電力を示している。図４（ａ）の対象時刻データは、消費電力の時系列データのうち、現在時刻＝"午前８時"として１９日分の時系列データを抽出することで生成されたものである。
【００３９】
具体的には、現在時刻が属する日の１日前の同一時刻"午前８時"における消費電力と、２日前の同一時刻"午前８時"における消費電力と、・・・、１８日前の同一時刻"午前８時"における消費電力とを抽出することで生成されたものである。
【００４０】
なお、図４（ａ）では、現在時刻＝"午前８時"として説明したが、現在時刻＝"午前９時"になれば、"午前９時"における１９日分の時系列データに基づいて、対象時刻データが生成される。
【００４１】
図４（ｂ）は、判定用データ生成部１０２により生成された第１の対象時刻変動率データの一例を示す図であり、横軸は日付及び時刻を示し、縦軸は消費電力の変動率Ｒ１を示している。第１の対象時刻変動率データは、対象時刻データに基づいて生成される。例えば、現在時刻＝"午前８時"の消費電力をＶａ、現在時刻の１時間前（＝"午前７時"）の消費電力をＶｂとすると、現在時刻における変動率Ｒ１は、下式のように算出される。
【００４２】
【数１】
0006558015.tif 000002
【００４３】
同様に、現在時刻が属する日の１日前の同一時刻"午前８時"の消費電力と、その１時間前の消費電力とに基づいて、現在時刻が属する日の１日前の同一時刻"午前８時"における変動率Ｒ１を算出する。以下、現在時刻が属する日の１８日前の変動率Ｒ１まで算出することで、図４（ｂ）の第１の対象時刻変動率データを得る。
【００４４】
図４（ｃ）は、判定用データ生成部１０２により生成された第２の対象時刻変動率データの一例を示す図であり、横軸は日付及び時刻を示し、縦軸は消費電力の変動率Ｒ２を示している。第２の対象時刻変動率データは、対象時刻データに基づいて生成される。例えば、現在時刻＝"午前８時"の消費電力をＶａ'、現在時刻の２時間前（＝"午前６時"）の消費電力をＶｂ'とすると、現在時刻における変動率Ｒ２は、下式のように算出される。
【００４５】
【数２】
0006558015.tif 000003
【００４６】
同様に、現在時刻が属する日の１日前の同一時刻"午前８時"の消費電力と、その２時間前の消費電力とに基づいて、現在時刻が属する日の１日前の同一時刻"午前８時"における変動率Ｒ２を算出する。以下、現在時刻が属する日の１８日前の変動率Ｒ２まで算出することで、図４（ｃ）の第２の対象時刻変動率データを得る。
【００４７】
次に、判定用データ生成部１０２による判定用データ生成処理の流れについて説明する。図５は、判定用データ生成部による判定用データ生成処理のフローチャートである。
【００４８】
ステップＳ５０１において、判定用データ生成部１０２は、時系列データ格納部１１１に格納されている、消費電力の時系列データを読み出す。なお、時系列データ格納部１１１には、既に現在時刻における消費電力のデータを含む時系列データが格納されているものとする。
【００４９】
ステップＳ５０２において、判定用データ生成部１０２は、読み出した時系列データから、現在時刻における消費電力のデータと、現在時刻が属する日の１日前から１８日前までの現在時刻と同一時刻における消費電力のデータとを抽出する。更に、抽出した消費電力のデータを、対象時刻データとして対象時刻データ格納部１１２に格納する。
【００５０】
ステップＳ５０３において、判定用データ生成部１０２は、現在時刻における変動率Ｒ１のデータ及び現在時刻が属する日の１日前から１８日前までの現在時刻と同一時刻における変動率Ｒ１のデータを算出する。更に、算出した変動率Ｒ１のデータを、第１の対象時刻変動率データとして第１の対象時刻変動率データ格納部１１３に格納する。
【００５１】
ステップＳ５０４において、判定用データ生成部１０２は、現在時刻における変動率Ｒ２のデータ及び現在時刻が属する日の１日前から１８日前までの現在時刻と同一時刻における変動率Ｒ２のデータを算出する。更に、算出した変動率Ｒ２のデータを、第２の対象時刻変動率データとして第２の対象時刻変動率データ格納部１１４に格納する。
【００５２】
ステップＳ５０５において、判定用データ生成部１０２は、ステップＳ５０１において時系列データを読み出してから、所定周期（例えば、１時間）が経過したか否かを判定する。ステップＳ５０５において、所定周期が経過していないと判定した場合には、所定周期が経過するまで待機する。
【００５３】
一方、ステップＳ５０５において、所定周期が経過したと判定した場合には、ステップＳ５０６に進み、判定用データ生成処理の終了指示が入力されていないか否かを判定し、入力していないと判定した場合には、ステップＳ５０１に戻る。
【００５４】
これにより、判定用データ生成処理の終了指示が入力されるまでは、所定周期ごとに、対象時刻データと、第１の対象時刻変動率データと、第２の対象時刻変動率データとが順次生成されていく。一方、判定用データ生成処理の終了指示が入力された場合には、判定用データ生成処理を終了する。
【００５５】
次に、判定用データ生成部１０２により所定周期ごとに生成される判定用データ（対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データ）について説明する。図６は、判定用データ生成部により所定周期ごとに判定用データが生成される様子を示した図である。
【００５６】
図６に示すように、例えば、現在時刻＝"午前８時"であった場合、判定用データ生成部１０２では、対象時刻データ６０１と、第１の対象時刻変動率データ６０２と、第２の対象時刻変動率データ６０３とを生成する。
【００５７】
また、所定周期（例えば、１時間）が経過し、現在時刻＝"午前９時"になると、判定用データ生成部１０２では、対象時刻データ６１１と、第１の対象時刻変動率データ６１２と、第２の対象時刻変動率データ６１３とを生成する。
【００５８】
更に、所定周期が経過し、現在時刻＝"午前１０時"になると、判定用データ生成部１０２では、対象時刻データ６２１と、第１の対象時刻変動率データ６２２と、第２の対象時刻変動率データ６２３とを生成する。
【００５９】
このように、判定用データ生成部１０２では、所定周期ごとに判定用データを生成していく。
【００６０】
次に、消費電力判定部１０３について説明する。図７は、消費電力判定部の詳細な機能構成の一例を示す図である。図７に示すように、消費電力判定部１０３は、データ検定部７０１、変化発生有無判定部７０２、外れ値回数算出部７０３、状態変化度算出部７０４、判定部７０５を有する。
【００６１】
データ検定部７０１は、所定周期ごとに、対象時刻データ格納部１１２に格納された対象時刻データのうち、現在時刻と同一時刻についての対象時刻データを読み出し、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行う。また、所定周期ごとに、第１の対象時刻変動率データ格納部１１３に格納された第１の対象時刻変動率データのうち、現在時刻と同一時刻についての第１の対象時刻変動率データを読み出し、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行う。更に、所定周期ごとに、第２の対象時刻変動率データ格納部１１４に格納された第２の対象時刻変動率データのうち、現在時刻と同一時刻についての第２の対象時刻変動率データを読み出し、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行う。
【００６２】
変化発生有無判定部７０２は、データ検定部７０１において行われた、対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データに対する検定の結果に基づいて、消費電力の変化が発生したか否かを判定する。具体的には、変化発生有無判定部７０２は、対象時刻データに対する検定の結果に基づいて、現在時刻における消費電力のデータが統計的外れ値に該当するか否かを判定し、統計的外れ値に該当すると判定した場合に、消費電力の変化が発生したと判定する。また、第１の対象時刻変動率データに対する検定の結果に基づいて、現在時刻における変動率Ｒ１のデータが統計的外れ値に該当するか否かを判定し、統計的外れ値に該当すると判定した場合に、消費電力の変化が発生したと判定する。更に、第２の対象時刻変動率データに対する検定の結果に基づいて、現在時刻における変動率Ｒ２のデータが統計的外れ値に該当するか否かを判定し、統計的外れ値に該当すると判定した場合に、消費電力の変化が発生したと判定する。
【００６３】
外れ値回数算出部７０３は、変化発生有無判定部７０２において検定が行われ消費電力の変化が発生したと判定された場合に、対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データそれぞれについて、統計的外れ値の回数を算出する。
【００６４】
具体的には、外れ値回数算出部７０３は、現在時刻における消費電力のデータが、現在時刻が属する日の後日（所定日数Ｍ−１日分）において、統計的外れ値であると判定された回数を算出する。例えば、所定日数Ｍ＝５の場合、現在時刻が属する日の４日後までの間に、現在時刻における消費電力のデータが統計的外れ値であると判定された回数（外れ値判定回数）を算出する。
【００６５】
同様に、外れ値回数算出部７０３は、現在時刻における変動率Ｒ１のデータが、現在時刻が属する日の後日（所定日数Ｍ−１日分）において、統計的外れ値であると判定された外れ値判定回数を算出する。更に、外れ値回数算出部７０３は、現在時刻における変動率Ｒ２のデータが、現在時刻が属する日の後日（所定日数Ｍ−１日分）において、統計的外れ値であると判定された外れ値判定回数を算出する。
【００６６】
状態変化度算出部７０４は、対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データそれぞれについて、外れ値回数算出部７０３において統計的外れ値であると判定された外れ値判定回数に基づいて、現在時刻についての状態変化度を算出する。
【００６７】
状態変化度とは、統計的外れ値であるとの判定結果が、以降、どの程度継続するかという継続性の度合いを、統計的外れ値であると判定された回数に基づいて数値化したものである。ここで、消費電力の変化が発生したと判定された現在時刻をＴ'、統計的外れ値であると判定された外れ値判定回数をＯＬ（Ｔ'）とおく。この場合、Ｋ日（Ｋ≦Ｍ−１、例えばＭ＝５）経過時の状態変化度（ｓｔ＿ｃｈｇ＿ｍｅａｓｕｒｅ（Ｔ'；Ｋ））は、下式により算出される。
【００６８】
【数３】
0006558015.tif 000004
状態変化度算出部７０４では、対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データそれぞれについて算出した、消費電力の変化が発生したと判定された現在時刻についての状態変化度を判定部７０５に通知する。
【００６９】
判定部７０５は、状態変化度算出部７０４により通知された状態変化度に基づいて、消費電力の変化が、状態変化の発生によるものなのか、そのとき限りの一時的な異常発生によるものなのかを判定する。つまり、判定部７０５は変化の属性を判定する。
【００７０】
ここで、状態変化度は０から１の範囲内で変化し、１に近いほど、状態変化の発生によるものであった可能性が高いと判定する。一方、状態変化度が０に近いほど、そのとき限りの一時的な異常発生によるものであった可能性が高いと判定する。
【００７１】
なお、消費電力判定部１０３により判定された判定結果は、判定結果通知部１０４により時系列データとともに表示装置１１０に送信される。ここで、判定結果には、消費電力の変化が発生したこと（現在時刻の消費電力のデータが統計的外れ値に該当すること）、及び、変化の属性が含まれる。また、判定結果には、消費電力の変化が発生したと判定された時刻についての状態変化度の値が含まれる。なお、判定結果は、電力制御装置１２０にも送信される。
【００７２】
次に、消費電力判定部１０３による変化発生有無判定処理の流れについて説明する。図８は、消費電力判定部による変化発生有無判定処理のフローチャートである。
【００７３】
ステップＳ８０１において、データ検定部７０１は、対象時刻データ格納部１１２に格納された対象時刻データのうち、現在時刻と同一時刻についての対象時刻データを読み出す。また、データ検定部７０１は、第１の対象時刻変動率データ格納部１１３に格納された第１の対象時刻変動率データのうち、現在時刻と同一時刻についての第１の対象時刻変動率データを読み出す。更に、データ検定部７０１は、第２の対象時刻変動率データ格納部１１４に格納された第２の対象時刻変動率データのうち、現在時刻と同一時刻についての第２の対象時刻変動率データを読み出す。
【００７４】
ステップＳ８０２において、データ検定部７０１は、ステップＳ８０１において読み出した対象時刻データに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行う。更に、変化発生有無判定部７０２は、検定結果に基づいて、現在時刻における消費電力のデータが統計的外れ値に該当するか否かを判定する。
【００７５】
ステップＳ８０３において、データ検定部７０１は、ステップＳ８０１において読み出した第１の対象時刻変動率データに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行う。更に、変化発生有無判定部７０２は、検定結果に基づいて、現在時刻における変動率Ｒ１が統計的外れ値に該当するか否かを判定する。
【００７６】
ステップＳ８０４において、データ検定部７０１は、ステップＳ８０１において読み出した第２の対象時刻変動率データに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行う。更に、変化発生有無判定部７０２は、検定結果に基づいて、現在時刻における変動率Ｒ２が統計的外れ値に該当するか否かを判定する。
【００７７】
ステップＳ８０５において、変化発生有無判定部７０２は、ステップＳ８０２からステップＳ８０４における判定の結果、統計的外れ値であると判定されたデータがあるか否かを判定する。ステップＳ８０５において、統計的外れ値であると判定されたデータがないと判定された場合には、ステップＳ８０７に進む。
【００７８】
一方、ステップＳ８０５において、統計的外れ値であると判定されたデータがあると判定された場合には、ステップＳ８０６に進む。ステップＳ８０６において、変化発生有無判定部７０２は、現在時刻において、消費電力の変化が発生したと判定する。
【００７９】
ステップＳ８０７において、データ検定部７０１は、ステップＳ８０１における判定用データの読み出し（対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データの読み出し）から所定周期が経過したか否かを判定する。ステップＳ８０７において、所定周期が経過していないと判定した場合には、所定周期が経過するまで待機する。
【００８０】
一方、ステップＳ８０７において、所定周期が経過したと判定した場合には、ステップＳ８０８に進み、変化発生有無判定処理の終了指示が入力されたか否かを判定し、終了指示が入力されていないと判定された場合には、ステップＳ８０１に戻る。この場合、所定周期が経過した後の現在時刻において、ステップＳ８０１からステップＳ８０７までの処理を行う。
【００８１】
一方、ステップＳ８０８において、変化発生有無判定処理の終了指示が入力されたと判定した場合には、変化発生有無判定処理を終了する。
【００８２】
次に、消費電力判定部１０３による変化発生有無判定処理の具体例について説明する。図９は、変化発生有無判定処理の具体例を示す図である。このうち、図９（ａ）は、対象時刻データ格納部１１２に格納された対象時刻データのうち、現在時刻と同一時刻についての対象時刻データを示している。図９（ａ）において、点線枠９０１は、現在時刻における消費電力のデータ及び過去１８日分の消費電力のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。図９（ａ）に示すように、現在時刻における消費電力のデータは、点線枠９０１の外側に位置しているため、統計的外れ値であると判定される。
【００８３】
図９（ｂ）は、第１の対象時刻変動率データ格納部１１３に格納された第１の対象時刻変動率データのうち、現在時刻と同一時刻についての第１の対象時刻変動率データを示している。図９（ｂ）において、点線枠９０２は、現在時刻における変動率Ｒ１のデータ及び過去１８日分の変動率Ｒ１のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。図９（ｂ）に示すように、現在時刻における変動率Ｒ１のデータは、点線枠９０２に含まれるため、統計的外れ値でないと判定される。
【００８４】
図９（ｃ）は、第２の対象時刻変動率データ格納部１１４に格納された第２の対象時刻変動率データのうち、現在時刻と同一時刻についての第２の対象時刻変動率データを示している。図９（ｃ）において、点線枠９０３は、現在時刻における変動率Ｒ２のデータ及び過去１８日分の変動率Ｒ２のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。図９（ｃ）に示すように、現在時刻における変動率Ｒ２のデータは、点線枠９０３の外側に位置しているため、統計的外れ値であると判定される。
【００８５】
このように、図９の例では、対象時刻データ及び第２の対象時刻変動率データにおいて、現在時刻におけるデータが統計的外れ値であると判定された。この結果、変化発生有無判定部７０２では、現在時刻において消費電力の変化が発生したと判定する。
【００８６】
次に、消費電力判定部１０３による状態変化度算出処理の流れについて説明する。図１０は、消費電力判定部による状態変化度算出処理のフローチャートである。なお、以降の説明において、"現在時刻"とは、変化発生有無判定部７０２において、消費電力の変化が発生したと判定された時刻（例えば"午前８時"）を指すものとする。
【００８７】
ステップＳ１００１において、状態変化度算出部７０４は、現在時刻において変化発生有無判定部７０２にて消費電力の変化が発生したと判定されたか否かを判定する。
【００８８】
消費電力の変化が発生したと判定されていなかった場合には、現在時刻についての状態変化度算出処理を終了する。一方、消費電力の変化が発生したと判定されていた場合には、ステップＳ１００２に進む。
【００８９】
ステップＳ１００２において、状態変化度算出部７０４は、経過日数Ｋに初期値である"１"を代入する。更に、ステップＳ１００３において、状態変化度算出部７０４は、現在時刻からＫ日が経過したか否かを判定する。
【００９０】
ステップＳ１００３において、Ｋ日が経過していないと判定された場合には、Ｋ日が経過したと判定されるまで待機する。一方、Ｋ日が経過したと判定された場合には、ステップＳ１００４に進む。
【００９１】
ステップＳ１００４において、状態変化度算出部７０４は、現在時刻が属する日の１日後からＫ日後までのそれぞれの対象時刻データに基づいて、現在時刻についての消費電力のデータが、統計的外れ値であると判定された外れ値判定回数を算出する。
【００９２】
ステップＳ１００５において、状態変化度算出部７０４は、ステップＳ１００４において算出された外れ値判定回数に基づいて、現在時刻についての状態変化度を算出する。これにより、現在時刻についてのＫ日経過時点での対象時刻データに基づく状態変化度が算出される。
【００９３】
ステップＳ１００６において、状態変化度算出部７０４は、現在時刻が属する日の１日後からＫ日後までのそれぞれの第１の対象時刻変動率データに基づいて、現在時刻についての変動率Ｒ１のデータが、統計的外れ値であると判定された外れ値判定回数を算出する。
【００９４】
ステップＳ１００７において、状態変化度算出部７０４は、ステップＳ１００６において算出された外れ値判定回数に基づいて、現在時刻についての状態変化度を算出する。これにより、現在時刻についてのＫ日経過時点での第１の対象時刻変動率データに基づく状態変化度が算出される。
【００９５】
ステップＳ１００８において、状態変化度算出部７０４は、現在時刻が属する日の１日後からＫ日後までのそれぞれの第２の対象時刻変動率データに基づいて、現在時刻についての変動率Ｒ２のデータが、統計的外れ値であると判定された外れ値判定回数を算出する。
【００９６】
ステップＳ１００９において、状態変化度算出部７０４は、ステップＳ１００８において算出された外れ値判定回数に基づいて、現在時刻についての状態変化度を算出する。これにより、現在時刻についてのＫ日経過時点での第２の対象時刻変動率データに基づく状態変化度が算出される。
【００９７】
ステップＳ１０１０において、状態変化度算出部７０４は、経過日数Ｋが所定の閾値（＝Ｍ−１、例えばＭ＝５）以下であるか否かを判定する。
【００９８】
ステップＳ１０１０において経過日数Ｋが所定の閾値以下であると判定された場合には、ステップＳ１０１１に進む。ステップＳ１０１１において、状態変化度算出部７０４は、経過日数Ｋをインクリメントした後、ステップＳ１００３に戻る。これにより、経過日数Ｋが所定の閾値を超えるまで、現在時刻における消費電力、変動率Ｒ１、変動率Ｒ２について、ステップＳ１００３からステップＳ１０１１までの処理が繰り返される。
【００９９】
一方、ステップＳ１０１０において、経過日数Ｋが所定の閾値を超えたと判定された場合には、状態変化度算出処理を終了する。
【０１００】
次に、消費電力判定部１０３による状態変化度算出処理の具体例について説明する。図１１から図１３は、状態変化度算出処理の具体例を示す図であり、図１１〜図１３は、いずれも現在時刻において、判定用データのいずれかで消費電力の変化が発生したと判定されたケースを示している。このうち、図１１は、現在時刻が属する日の後日（１日後）の同一時刻において変化発生有無判定処理が行われた結果、現在時刻における"消費電力の変化が発生したとの判定"が解消したケースを示している。一方、図１２は、現在時刻が属する日の後日（１日後）の同一時刻において変化発生有無判定処理が行われた結果、現在時刻における"消費電力の変化が発生したとの判定"が継続したケースを示している。更に、図１３は、現在時刻が属する日の後日（１日後）の同一時刻において変化発生有無判定処理が行われた結果、現在時刻における"消費電力の変化が発生したとの判定"が継続したケースを示している。なお、図１３のケースは、更に、１日後の同一時刻において"消費電力の変化が新たに発生したとの判定"がなされたケースを示している。以下、図１１から順に説明する。
【０１０１】
図１１（ａ）は、現在時刻が属する日の１日後の同一時刻において対象時刻データを読み出し、変化発生有無判定処理が行われた様子を示している。点線枠１１０１は、現在時刻が属する日の１日後の同一時刻における消費電力のデータ及び過去Ｎ日分（Ｎ＝１８）の消費電力のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。
【０１０２】
図１１（ａ）に示すように、現在時刻における消費電力のデータは、点線枠１１０１に含まれるため、１日後の同一時刻における変化発生有無判定処理では、統計的外れ値ではないと判定される。つまり、現在時刻における消費電力のデータは、現在時刻においては点線枠９０１の外側に位置していたため、統計的外れ値であると判定されたが、現在時刻が属する日の１日後の同一時刻においては、統計的外れ値であるとの判定は解消されたことになる。
【０１０３】
このため、現在時刻における消費電力のデータについての１日経過時点での外れ値判定回数は"１"となり、状態変化度算出部７０４では、当該外れ値判定回数を用いて、対象時刻データに基づく現在時刻についての１日経過時点での状態変化度を算出する。
【０１０４】
図１１（ｂ）は、現在時刻が属する日の１日後の同一時刻において第１の対象時刻変動率データを読み出し、変化発生有無判定処理が行われた様子を示している。点線枠１１０２は、現在時刻が属する日の１日後の同一時刻における変動率Ｒ１のデータ及び過去Ｎ日分（Ｎ＝１８）の変動率Ｒ１のデータについて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行い、統計的外れ値でないと判定された範囲を示している。
【０１０５】
図１１（ｂ）に示すように、現在時刻における変動率Ｒ１のデータは、点線枠１１０２に含まれるため、再び、統計的外れ値ではないと判定される。このため、第１の対象時刻変動率データに基づく現在時刻についての状態変化度は算出されない。
【０１０６】
図１１（ｃ）は、現在時刻が属する日の１日後の同一時刻において第２の対象時刻変動率データを読み出し、変化発生有無判定処理が行われた様子を示している。点線枠１１０３は、現在時刻が属する日の１日後の同一時刻における変動率Ｒ２のデータ及び過去Ｎ日分（Ｎ＝１８）の変動率Ｒ２のデータについて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行い、統計的外れ値でないと判定された範囲を示している。
【０１０７】
図１１（ｃ）に示すように、現在時刻における変動率Ｒ２のデータは、現在時刻においては点線枠９０３の外側に位置していたため、統計的外れ値であると判定された。しかし、現在時刻が属する日の１日後の同一時刻においては、統計的外れ値であるとの判定は解消されることになる。
【０１０８】
このため、現在時刻における変動率Ｒ２のデータについての１日経過時点での外れ値判定回数は"１"となる。状態変化度算出部７０４では、当該外れ値判定回数を用いて、第２の対象時刻変動率データに基づく現在時刻についての１日経過時点での状態変化度を算出する。
【０１０９】
図１２（ａ）は、現在時刻が属する日の１日後の同一時刻において対象時刻データを読み出し、変化発生有無判定処理が行われた他の様子を示している。点線枠１２０１は、現在時刻が属する日の１日後の同一時刻における消費電力のデータ及び過去Ｎ日分（Ｎ＝１８）の消費電力のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。
【０１１０】
図１２（ａ）に示すように、現在時刻における消費電力のデータは、点線枠１２０１の外側に位置するため、統計的外れ値であると判定される。つまり、現在時刻における対象時刻データは、現在時刻において点線枠９０１の外側に位置していたため、統計的外れ値であると判定された。そして、現在時刻が属する日の１日後の同一時刻においても依然として点線枠１２０１の外側に位置しているため、統計的外れ値であると判定され、統計的外れ値であるとの判定は継続されることになる。
【０１１１】
このため、現在時刻における消費電力のデータについての１日経過時点での外れ値判定回数は"２"となり、状態変化度算出部７０４では、当該外れ値判定回数を用いて、対象時刻データに基づく現在時刻についての１日経過時点での状態変化度を算出する。
【０１１２】
図１２（ｂ）は、現在時刻が属する日の１日後の同一時刻において第１の対象時刻変動率データを読み出し、変化発生有無判定処理が行われた他の様子を示している。点線枠１２０２は、現在時刻が属する日の１日後の同一時刻における変動率Ｒ１のデータ及び過去Ｎ日分（Ｎ＝１８）の変動率Ｒ１のデータについて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行い、統計的外れ値でないと判定された範囲を示している。
【０１１３】
図１２（ｂ）に示すように、現在時刻における変動率Ｒ１のデータは、点線枠１２０２に含まれるため、再び、統計的外れ値ではないと判定される。このため、第１の対象時刻変動率データに基づく現在時刻についての状態変化度は算出されない。
【０１１４】
図１２（ｃ）は、現在時刻が属する日の１日後の同一時刻において第２の対象時刻変動率データを読み出し、変化発生有無判定処理が行われた他の様子を示している。点線枠１２０３は、現在時刻が属する日の１日後の同一時刻における変動率Ｒ２のデータ及び過去Ｎ日分（Ｎ＝１８）の変動率Ｒ２のデータについて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行い、統計的外れ値でないと判定された範囲を示している。
【０１１５】
図１２（ｃ）に示すように、現在時刻における変動率Ｒ２のデータは、点線枠１２０３の外側に位置するため、統計的外れ値であると判定される。つまり、現在時刻における変動率Ｒ２のデータは、現在時刻においては点線枠９０３の外側に位置していたため、統計的外れ値であると判定された。そして、現在時刻が属する日の１日後の同一時刻においても依然として点線枠１２０３の外側に位置しているため、統計的外れ値であると判定され、統計的外れ値であるとの判定は継続されることになる。
【０１１６】
このため、現在時刻における変動率Ｒ２のデータについての１日経過時点での外れ値判定回数は"２"となる。状態変化度算出部７０４では、当該外れ値判定回数を用いて、第２の対象時刻変動率データに基づく現在時刻についての１日経過時点での状態変化度を算出する。
【０１１７】
図１３（ａ）は、現在時刻が属する日の１日後の同一時刻において対象時刻データを読み出し、変化発生有無判定処理が行われた他の様子を示している。点線枠１３０１は、現在時刻が属する日の１日後の同一時刻における消費電力のデータ及び過去Ｎ日分（Ｎ＝１８）の消費電力のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。
【０１１８】
図１３（ａ）に示すように、現在時刻における消費電力のデータは、点線枠１３０１の外側に位置するため、統計的外れ値であると判定される。なお、現在時刻における消費電力のデータは、現在時刻においても点線枠９０１の外側に位置していたため、現在時刻が属する日の１日後の同一時刻において、現在時刻についての外れ値判定回数は２回となる。このため、状態変化度算出部７０４では、当該外れ値判定回数を用いて、対象時刻データに基づく現在時刻についての１日経過時点での状態変化度を算出する。
【０１１９】
また、現在時刻が属する日の１日後の同一時刻における消費電力のデータは、点線枠１３０１の外側に位置しているため、現在時刻が属する日の１日後の同一時刻における消費電力のデータは、統計的外れ値であると判定される。なお、現在時刻が属する日の１日後の同一時刻において、現在時刻が属する日の１日後の同一時刻における消費電力のデータについての外れ値判定回数は１回となる。このため、状態変化度算出部７０４では、当該外れ値判定回数を用いて、対象データに基づく現在時刻が属する日の１日後の同一時刻についての状態変化度を算出する。
【０１２０】
図１３（ｂ）は、現在時刻が属する日の１日後の同一時刻において第１の対象時刻変動率データを読み出し、変化発生有無判定処理が行われた他の様子を示している。点線枠１３０２は、現在時刻が属する日の１日後の同一時刻における変動率Ｒ１のデータ及び過去Ｎ日分（Ｎ＝１８）の変動率Ｒ１のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。
【０１２１】
図１３（ｂ）に示すように、現在時刻における変動率Ｒ１のデータは、点線枠１３０２に含まれるため、統計的外れ値でないと判定される。このため、第１の対象時刻変動率データに基づく現在時刻についての状態変化度は算出されない。
【０１２２】
図１３（ｃ）は、現在時刻が属する日の１日後の同一時刻において第２の対象時刻変動率データを読み出し、変化発生有無判定処理が行われた他の様子を示している。点線枠１３０３は、現在時刻が属する日の１日後の同一時刻における変動率Ｒ２のデータ及び過去Ｎ日分（Ｎ＝１８）の変動率Ｒ２のデータに基づいて、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行った結果、統計的外れ値でないと判定された範囲を示している。
【０１２３】
図１３（ｃ）に示すように、現在時刻における変動率Ｒ２のデータは、点線枠１３０３の外側に位置するため、統計的外れ値であると判定される。なお、現在時刻における変動率Ｒ２のデータは、現在時刻においても点線枠９０３の外側に位置していたため、現在時刻が属する日の１日後の同一時刻において、現在時刻についての外れ値判定回数は２回となる。このため、状態変化度算出部７０４では、当該外れ値判定回数を用いて、第２の対象時刻変動率データに基づく現在時刻についての１日経過時点での状態変化度を算出する。
【０１２４】
また、現在時刻が属する日の１日後の同一時刻における変動率Ｒ２のデータは、点線枠１３０３の外側に位置しているため、現在時刻が属する日の１日後の同一時刻における変動率Ｒ２のデータは、統計的外れ値であると判定される。なお、現在時刻が属する日の１日後の同一時刻において、現在時刻が属する日の１日後の同一時刻における変動率Ｒ２のデータについての外れ値判定回数は１回となる。このため、状態変化度算出部７０４では、当該外れ値判定回数を用いて、第２の対象時刻変動率データに基づく現在時刻が属する日の１日後の同一時刻についての状態変化度を算出する。
【０１２５】
次に、消費電力判定部１０３による判定処理の流れについて説明する。図１４、図１５は、消費電力判定部による判定処理のフローチャートである。ステップＳ１４０１において、判定部７０５は、現在時刻において変化発生有無判定部７０２にて消費電力の変化が発生したと判定されたか否かを判定する。
【０１２６】
消費電力の変化が発生したと判定されていなかった場合には、現在時刻についての判定処理を終了する。一方、消費電力の変化が発生したと判定されていた場合には、ステップＳ１４０２に進む。
【０１２７】
ステップＳ１４０２において、経過日数Ｋに初期値である"１"を代入する。更に、ステップＳ１４０３において、判定部７０５は、現在時刻からＫ日が経過したか否かを判定する。
【０１２８】
ステップＳ１４０３において、Ｋ日が経過していないと判定された場合には、Ｋ日が経過したと判定されるまで待機する。一方、Ｋ日が経過したと判定された場合には、ステップＳ１４０４に進む。
【０１２９】
ステップＳ１４０４において、判定部７０５は、対象時刻データについて、現在時刻からＫ日経過した時点での、現在時刻についての状態変化度を状態変化度算出部７０４から取得する。
【０１３０】
ステップＳ１４０５において、判定部７０５は、ステップＳ１４０４において取得した状態変化度が所定の閾値以上であるか否かを判定する。ステップＳ１４０５において所定の閾値以上でないと判定された場合には、ステップＳ１４０６に進む。
【０１３１】
ステップＳ１４０６において、判定部７０５は、現在時刻における消費電力の変化は、そのとき限りの一時的な異常発生によるものであると判定する。
【０１３２】
一方、ステップＳ１４０５において所定の閾値以上であると判定された場合には、ステップＳ１４０７に進む。
【０１３３】
ステップＳ１４０７において、判定部７０５は、現在時刻における消費電力の変化は、状態変化の発生によるものであると判定する。
【０１３４】
図１５のステップＳ１５０１において、判定部７０５は、第１の対象時刻変動率データについて、現在時刻からＫ日経過した時点での、現在時刻についての状態変化度を状態変化度算出部７０４から取得する。
【０１３５】
ステップＳ１５０２において、判定部７０５は、ステップＳ１５０１において取得した状態変化度が所定の閾値以上であるか否かを判定する。ステップＳ１５０２において所定の閾値以上でないと判定された場合には、ステップＳ１５０３に進む。
【０１３６】
ステップＳ１５０３において、判定部７０５は、現在時刻における変動率Ｒ１の変化は、そのとき限りの一時的な異常発生によるものであると判定する。
【０１３７】
一方、ステップＳ１５０２において所定の閾値以上であると判定された場合には、ステップＳ１５０４に進む。
【０１３８】
ステップＳ１５０４において、判定部７０５は、現在時刻における変動率Ｒ１の変化は、状態変化の発生によるものであると判定する。
【０１３９】
ステップＳ１５０５において、判定部７０５は、第２の対象時刻変動率データについて、現在時刻からＫ日経過した時点での、現在時刻についての状態変化度を状態変化度算出部７０４から取得する。
【０１４０】
ステップＳ１５０６において、判定部７０５は、ステップＳ１５０５において取得した状態変化度が所定の閾値以上であるか否かを判定する。ステップＳ１５０６において所定の閾値以上でないと判定された場合には、ステップＳ１５０７に進む。
【０１４１】
ステップＳ１５０７において、判定部７０５は、現在時刻における変動率Ｒ２の変化は、そのとき限りの一時的な異常発生によるものであると判定する。
【０１４２】
一方、ステップＳ１５０６において所定の閾値以上であると判定された場合には、ステップＳ１５０８に進む。
【０１４３】
ステップＳ１５０８において、判定部７０５は、現在時刻における変動率Ｒ２の変化は、状態変化の発生によるものであると判定する。
【０１４４】
ステップＳ１５０９において、判定部７０５は、経過日数Ｋが所定の閾値（＝Ｍ−１、例えばＭ＝５）以下であるか否かを判定する。
【０１４５】
ステップＳ１５０９において経過日数Ｋが所定の閾値以下であると判定された場合には、ステップＳ１５１０に進む。ステップＳ１５１０において、判定部７０５は、経過日数Ｋをインクリメントした後、ステップＳ１４０３に戻る。これにより、経過日数Ｋが所定の閾値を超えるまで、現在時刻における消費電力、変動率Ｒ１、変動率Ｒ２について、ステップＳ１４０３からステップＳ１５０８の処理が繰り返される。
【０１４６】
一方、ステップＳ１５０９において、経過日数が所定の閾値を超えたと判定された場合には、判定処理を終了する。
【０１４７】
次に、消費電力の変化が発生してからの経過日数Ｋに対する状態変化度の推移について説明する。図１６は、消費電力の変化が発生してからの経過日数に対する状態変化度の推移を示す図である。
【０１４８】
図１６において、横軸は消費電力の変化が発生してからの経過日数を示しており、縦軸は状態変化度を示している。なお、ここでは、対象時刻データに基づいて算出された状態変化度の推移について説明する（第１の対象時刻変動率データ、第２の対象時刻変動率データについても同様であるため、ここでは説明は省略する）。
【０１４９】
図１６に示すように、消費電力の変化が発生した時刻で状態変化度算出部７０４により算出される状態変化度は０．５であり、その後、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定が行われるごとに変化していく。
【０１５０】
具体的には、消費電力の変化が発生した現在時刻に対して１日後の同一時刻において、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行うことで、再び、統計的外れ値であると判定された場合には、状態変化度は０．５より小さくなる。更に、２日後の同一時刻において、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行うことで、再び、統計的外れ値であると判定された場合には、状態変化度は更に小さくなる（図１６の黒四角点参照）。
【０１５１】
このように、消費電力の変化が発生した現在時刻における消費電力のデータが、統計的外れ値であると判定され続けると、状態変化度は０に近づいていく。この結果、消費電力の変化は、そのとき限りの一時的な異常発生によるものであったと判定することができる。
【０１５２】
一方、消費電力の変化が発生した現在時刻に対して１日後の同一時刻において、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行うことで、統計的外れ値でないと判定された場合には、状態変化度は０．５より大きくなる。更に、２日後の同一時刻において、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行うことで、統計的外れ値でないと判定された場合には、状態変化度は更に大きくなる（図１６の白四角点参照）。
【０１５３】
このように、消費電力の変化が発生した現在時刻における消費電力のデータについて、その後、統計的外れ値であるとの判定が解消されると、状態変化度は１に近づいていく。この結果、現在時刻における消費電力の変化は、状態変化の発生によるものであったと判定することができる。
【０１５４】
次に、判定結果通知部１０４により時系列データとともに通知され、表示装置１１０に表示される判定結果について説明する。図１７は、判定結果通知部により通知された判定結果の一例を示す図である。
【０１５５】
図１７に示すように、表示装置１１０には、判定結果通知部１０４により通知された時系列データが表示される。具体的には、予め定められた日数分の消費電力の時系列データが表示される。
【０１５６】
また、表示装置１１０には、判定結果通知部１０４により通知された判定結果が時系列データと対応付けて表示される。具体的には、時系列データの各時刻における状態変化度が、時系列データの各時刻に対応する位置に配された直線領域に、状態変化度の大きさに応じて動作する点（識別子）としてプロットされる。
【０１５７】
図１７において、時系列データに対応付けてプロットされた点は、状態変化度を示している。上述した通り、状態変化度は０から１．０の間の値となるため、各時刻における状態変化度は、図１７の右端に示された目盛に則した位置にプロットされる。
【０１５８】
図１６を用いて説明した通り、状態変化度は、消費電力の変化が発生した時点では０．５となるが、その後、日数が経過するにつれて、０または１に近づいていく。したがって、図１７に示すように、０から１の範囲内において状態変化度を示す点がプロットされることになる。
【０１５９】
このように、表示装置１１０には、消費電力の時系列データに加え、状態変化度が表示される。このため、電力制御装置１２０の電力管理担当者は、最新の時系列データが表示された際に、状態変化度を表す点がプロットされたか否かを確認することで、消費電力の変化が発生したか否かを直ちに認識することができる。
【０１６０】
また、状態変化度を表す点がプロットされた後、１に近づいていくのか０に近づいていくのかを監視することで、状態変化度を表す点がプロットされたのが、状態変化の発生によるものであったのか否かを判断することができる。
【０１６１】
図１７の拡大表示１７０１は、消費電力の変化が発生したと判定され、その後、日数が経過するにつれ、状態変化度が１．０になった例を示している。この場合、電力管理担当者は、状態変化度を表す点がプロットされた時点で、消費電力の変化が発生したことを認識する。その後、日数が経過することで、電力管理担当者は、当該消費電力の変化が、状態変化の発生によるものであったと認識することができる。
【０１６２】
一方、図１７の拡大表示１７０２は、消費電力の変化が発生したと判定され、その後、日数が経過するにつれ、状態変化度が０に近づいた例を示している。この場合、電力管理担当者は、状態変化度を表す点がプロットされた時点で、消費電力の変化が発生したことを認識する。その後、日数が経過することで、電力管理担当者は、当該消費電力の変化が、そのとき限りの一時的な異常発生によるものであったと認識することができる。
【０１６３】
以上の説明から明らかなように、本実施形態では、消費電力の時系列データに基づいて、現在時刻における消費電力のデータが、統計的外れ値であるか否かを判定し、統計的外れ値であると判定した場合に、消費電力の変化が発生したと判定する。
【０１６４】
このように、消費電力の日内変動の変動パターン全体に基づいて判定するのではなく、現在時刻における消費電力のデータに基づいて判定することで、本実施形態に係る消費電力監視装置によれば、状態変化の兆しをより早く検出することが可能になる。
【０１６５】
また、本実施形態に係る消費電力監視装置では、統計的外れ値であると判定された場合、その後、統計的外れ値であるとの判定がどの程度継続するのかを監視することで、統計的外れ値であるとの判定が、状態変化の発生によるものであったか否かを判定する。
【０１６６】
これにより、本実施形態に係る消費電力監視装置によれば、消費電力の変化が、状態変化の発生によるものであったのか、そのとき限りの一時的な異常発生によるものであったのか、変化の属性を判定することが可能となる。
【０１６７】
［第２の実施形態］
上記第１の実施形態では、判定用データ生成部１０２が、判定用データとして、対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データを生成した。しかしながら、消費電力の変化の発生有無の判定に用いる判定用データの組み合わせはこれに限定されない。例えば、このうちの１または２の判定用データを用いて消費電力の変化の発生有無を判定するようにしてもよい。あるいは、これら以外の判定用データを用いて消費電力の変化の発生有無を判定するようにしてもよい。
【０１６８】
また、上記第１の実施形態では、第１の対象時刻変動率データとして、１時間あたりの変動率を用いたが、第１の対象時刻変動率データは、１時間あたりの変動率に限定されない。同様に、上記第１の実施形態では、第２の対象時刻変動率データとして、２時間あたりの変動率を用いたが、第２の対象時刻変動率データは、２時間あたりの変動率に限定されない。
【０１６９】
また、上記第１の実施形態では、判定用データ生成部１０２が対象時刻データ、第１の対象時刻変動率データ、第２の対象時刻変動率データを生成するにあたり、過去１８日分のデータを用いた。しかしながら、判定用データの生成に用いる範囲はこれに限定されず、１８日より長くても短くてもよい。
【０１７０】
また、上記第１の実施形態では、状態変化度を算出する期間を、現在時刻から４日が経過するまでとしたが、状態変化度を算出する期間は、現在時刻から４日が経過するまでに限定されない。また、上記第１の実施形態では、状態変化度を判定する際の閾値について特に限定しなかったが、状態変化度を判定する際の閾値は、例えば、０．５に設定することができる。
【０１７１】
また、上記第１の実施形態では、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定を行うことで統計的外れ値であるか否かを判定した。しかしながら、統計的外れ値であるか否かの判定に用いる検定は、Ｓｍｉｒｎｏｖ−Ｇｒｕｂｂｓ検定に限定されず、他の検定方法であってもよい。
【０１７２】
また、上記第１の実施形態では、消費電力監視装置１００から表示装置１１０に送信された判定結果のうち、消費電力の変化が発生したことを、プロットにより表現することとした。また、状態変化度は、点がプロットされた位置により表現することとした。しかしながら、消費電力の変化が発生したことを示す情報の表示態様及び状態変化度の表示態様はこれに限定されない。また、判定結果のうち、変化の属性（状態変化の発生によるものか、そのとき限りの一時的な異常発生によるものなのか）についての表示態様は特に言及しなかったが、例えば、プロットする点の色または形を変更することで表現するようにしてもよい。
【０１７３】
また、上記第１の実施形態では、判定用データそれぞれにおいて、消費電力の変化が発生したと判定された場合のプロット方法について特に言及しなかったが、例えば、優先順位をつけてプロットするようにしてもよい。具体的には、判定用データのうち、対象時刻データと第２の対象時刻変動率データそれぞれにおいて、消費電力の変化が発生したと判定され、それぞれにおいて算出された状態変化度が異なっていたとする。この場合、例えば、対象時刻データにおいて算出された状態変化度を優先してプロットするようにしてもよい。
【０１７４】
また、上記第１の実施形態では、消費電力監視装置１００を表示装置１１０及び電力制御装置１２０の近傍に配置するものとして説明した。しかしながら、消費電力監視装置１００は、ネットワークを介して表示装置１１０及び電力制御装置１２０と接続し、表示装置１１０及び電力制御装置１２０に対して遠隔に配置してもよい。
【０１７５】
なお、開示の技術では、以下に記載する付記のような形態が考えられる。
（付記１）
現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定する判定部と、
前記判定部により所定の範囲から外れていると判定されたことを通知する通知部と
を有することを特徴とするデータ処理装置。
（付記２）
前記判定部は、更に、
現在時刻に対応付けられた消費電力の所定時間前からの変動率のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力の所定時間前からの変動率のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力の所定時間前からの変動率のデータが、所定の範囲から外れているか否かを判定し、
前記通知部は、
前記消費電力のデータが所定の範囲から外れていると判定されたこと、または、前記変動率のデータが所定の範囲から外れていると判定されたことを通知することを特徴とする付記１に記載のデータ処理装置。
（付記３）
現在時刻に対応付けられた消費電力の所定時間前からの変動率のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力の所定時間前からの変動率のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力の所定時間前からの変動率のデータが、所定の範囲から外れているか否かを判定する判定部と、
前記判定部により所定の範囲から外れていると判定されたことを通知する通知部と
を有することを特徴とするデータ処理装置。
（付記４）
現在時刻に対応付けられた消費電力に応じたデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力に応じたデータとを含むデータ分布において、現在時刻に対応付けられた消費電力に応じたデータが、前記所定の範囲から外れているか否かを判定する判定部と、
前記判定部により所定の範囲から外れていると判定されたことを通知する通知部と
を有することを特徴とするデータ処理装置。
（付記５）
前記判定部は、
前記データ分布において、統計的外れ値に該当するデータであるか否かを判定することで、前記所定の範囲から外れているか否か判定することを特徴とする付記１に記載のデータ処理装置。
（付記６）
前記判定部は、
前記所定の範囲から外れていると判定された場合に、現在時刻が属する日の後日における同一時刻に対応付けられた消費電力のデータを含む前記データ分布に基づいて、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れていると判定した回数を算出し、
算出した回数に基づいて、前記所定の範囲から外れているとの判定結果が継続する度合いを示す状態変化度を算出することを特徴とする付記５に記載のデータ処理装置。
（付記７）
前記通知部は、
前記状態変化度を、前記消費電力のデータと時刻とが対応付けられた時系列データに対応付けて通知することを特徴とする付記５に記載のデータ処理装置。
（付記８）
前記通知部は、
前記時系列データの各時刻に対応する位置に配された直線領域において、前記状態変化度に応じて動作する識別子をプロットすることで、前記状態変化度を通知することを特徴とする付記７に記載のデータ処理装置。
（付記９）
コンピュータに、
現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定し、
前記所定の範囲から外れていると判定されたことを通知する、
処理を実行させることを特徴とするデータ処理プログラム。
（付記１０）
現在時刻に対応付けられた消費電力のデータと、現在時刻と同一時刻に対応付けられた他の日の消費電力のデータとを含むデータ分布において、現在時刻に対応付けられた消費電力のデータが、所定の範囲から外れているか否かを判定し、
前記所定の範囲から外れていると判定されたことを通知する、
処理をコンピュータが実行することを特徴とするデータ処理方法。
【０１７６】
なお、上記実施形態に挙げた構成等に、その他の要素との組み合わせ等、ここで示した構成に本発明が限定されるものではない。これらの点に関しては、本発明の趣旨を逸脱しない範囲で変更することが可能であり、その応用形態に応じて適切に定めることができる。
【符号の説明】
【０１７７】
１００                  ：消費電力監視装置
１０１                  ：時系列データ取得部
１０２                  ：判定用データ生成部
１０３                  ：消費電力判定部
１０４                  ：判定結果通知部
１１０                  ：表示装置
１２０                  ：電力制御装置
１４０                  ：監視システム
６０１、６１１、６２１  ：対象時刻データ
６０２、６１２、６２２  ：第１の対象時刻変動率データ
６０３、６１３、６２３  ：第２の対象時刻変動率データ
７０１                  ：データ検定部
７０２                  ：変化発生有無判定部
７０３                  ：外れ値回数算出部
７０４                  ：状態変化度算出部
７０５                  ：判定部
【図１】
0006558015.tif 000005
【図２】
0006558015.tif 000006
【図３】
0006558015.tif 000007
【図４】
0006558015.tif 000008
【図５】
0006558015.tif 000009
【図６】
0006558015.tif 000010
【図７】
0006558015.tif 000011
【図８】
0006558015.tif 000012
【図９】
0006558015.tif 000013
【図１０】
0006558015.tif 000014
【図１１】
0006558015.tif 000015
【図１２】
0006558015.tif 000016
【図１３】
0006558015.tif 000017
【図１４】
0006558015.tif 000018
【図１５】
0006558015.tif 000019
【図１６】
0006558015.tif 000020
【図１７】
0006558015.tif 000021