(57)【要約】【課題】電力需要の予測精度の向上を図ることができる情報処理システム、および情報処理プログラムを提供することである。【解決手段】実施形態の情報処理システムは、個別予測モデル導出部と、グループ生成部とを持つ。個別予測モデル導出部は、複数の需要家の過去の個別の電力需要を示す個別需要データと、気象データとに基づき、前記需要家ごとの電力需要の個別予測モデルを導出する。グループ生成部は、前記個別予測モデル導出部によって導出された前記複数の需要家の前記個別予測モデルの相互の類似度に基づき前記複数の需要家を複数のグループに分ける。【選択図】図１
2018113817.tif 000002
【特許請求の範囲】
【請求項１】
複数の需要家の過去の個別の電力需要を示す個別需要データと、気象データとに基づき、前記需要家ごとの電力需要の個別予測モデルを導出する個別予測モデル導出部と、
前記個別予測モデル導出部によって導出された前記複数の需要家の前記個別予測モデルの相互の類似度に基づき前記複数の需要家を複数のグループに分けるグループ生成部と、
を備える情報処理システム。
【請求項２】
前記個別予測モデル導出部は、回帰モデルに基づき前記個別予測モデルを導出し、
前記グループ生成部は、前記複数の需要家の前記個別予測モデルに対して変数選択を実施し、前記個別予測モデルに含まれる変数の数を減らした状態で前記類似度を導出する
請求項１に記載の情報処理システム。
【請求項３】
前記個別予測モデル導出部は、回帰モデルに基づき前記個別予測モデルを導出し、
前記グループ生成部は、少なくとも前記個別予測モデルに含まれる変数の種類に基づき、前記類似度を導出する
請求項１または請求項２に記載の情報処理システム。
【請求項４】
前記グループ生成部は、前記変数の種類に加えて前記変数の重みに基づき、前記類似度を導出する
請求項３に記載の情報処理システム。
【請求項５】
前記個別予測モデル導出部は、回帰モデルに基づき前記個別予測モデルを導出し、
前記グループ生成部は、前記複数の需要家の前記個別予測モデルに含まれる変数の相互相関に基づき、前記類似度を導出する
請求項１から請求項４のいずれか１項に記載の情報処理システム。
【請求項６】
前記グループ生成部は、前記複数の需要家の前記個別予測モデルに含まれる変数の間に時間差成分がある場合、前記時間差成分にも基づき前記類似度を導出する
請求項３から請求項５のいずれか１項に記載の情報処理システム。
【請求項７】
前記グループ生成部により分けられたグループごとの前記個別需要データに関するデータと、気象データとに基づき、前記グループごとの電力需要の予測モデルを導出するグループ予測モデル導出部と、
前記グループ予測モデル導出部により導出された予測モデルに基づき、前記グループごとの電力需要を予測するグループ需要予測部と、
前記グループ需要予測部により予測された前記グループごとの電力需要に基づき、前記複数の需要家の全体需要を導出するグループデータ合成部と、
請求項１から請求項６のいずれか１項に記載の情報処理システム。
【請求項８】
前記個別予測モデル導出部は、回帰モデルに基づき前記個別予測モデルを導出し、
前記グループ予測モデル導出部は、前記個別予測モデルに含まれる係数によって重み付けされた前記需要家ごとの前記個別需要データが前記グループごとに纏められたデータに基づき、前記グループごとの電力需要の予測モデルを導出する
請求項７に記載の情報処理システム。
【請求項９】
異なるグループ数でグループ分けされた結果に基づいて予測された全体需要予測と電力需要の実測値との予測誤差を、分けられたグループの数ごとに導出する予測誤差導出部と、
前記予測誤差導出部により導出された複数の予測誤差に基づいて、予め決められた評価基準に従って予測誤差を評価し、評価結果に基づいてグループ設定値を決定するグループ設定値決定部と、
をさらに備える請求項７または請求項８に記載の情報処理システム。
【請求項１０】
コンピュータに、
複数の需要家の過去の個別の電力需要を示す個別需要データと、気象データとに基づき、前記需要家ごとの電力需要の個別予測モデルを導出させ、
導出された前記複数の需要家の前記個別予測モデルの相互の類似度に基づき前記複数の需要家を複数のグループに分けさせる、
情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、情報処理システム、および情報処理プログラムに関する。
【背景技術】
【０００２】
近年、電力小売が全面自由化され、多くの新電力会社（以後、単に新電力と記す）が事業に参入している。新電力は、電力を供給している需要家の総需要を予測し、広域機関に提出する必要がある。予測と実需要が乖離した場合、乖離度に応じたペナルティが発生し、新電力の収益を悪化させる。このため、需要予測の高精度化が要望されている。
【０００３】
スマートメーターなどの設置により、新電力顧客の需要カーブは個別に計測可能となった。これにより、総需要を予測する方法として、様々な戦略をとることが可能となった。第一の方法として、需要家の需要カーブを合計し、需要を予測するという方法がある。第二の方法として、需要家個別に需要を予測し、その合計の需要カーブを生成するという方法がある。第一の方法は、合計した需要カーブは個別の需要家のノイズが打ち消されて滑らかな形状になるという特徴がある。第二の方法は、需要家ごとの特性を個々の需要家の予測モデルに反映できるという特徴がある。両者の予測性能はケースバイケースである。ところで、電力需要の予測は、さらなる精度の向上が期待されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平１０−１１７４３７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明が解決しようとする課題は、予測精度の向上を図ることができる情報処理システム、および情報処理プログラムを提供することである。
【課題を解決するための手段】
【０００６】
実施形態の情報処理システムは、個別予測モデル導出部と、グループ生成部とを持つ。個別予測モデル導出部は、複数の需要家の過去の個別の電力需要を示す個別需要データと、気象データとに基づき、前記需要家ごとの電力需要の個別予測モデルを導出する。グループ生成部は、前記個別予測モデル導出部によって導出された前記複数の需要家の前記個別予測モデルの相互の類似度に基づき前記複数の需要家を複数のグループに分ける。
【図面の簡単な説明】
【０００７】
【図１】第１実施形態の情報処理システム１の構成図。
【図２】第１実施形態の個別需要データの一例を示す図。
【図３】第１実施形態のグループ分類部３２の動作の流れの一例を示すフローチャート。
【図４】第１実施形態の類似度テーブルの一例を示す図。
【図５】第２実施形態の情報処理システム１Ａの構成図。
【図６】第２実施形態の全体需要予測部３３の動作の流れの一例を示すフローチャート。
【図７】第３実施形態の類似度テーブルの他の例を示す図。
【図８】第４実施形態の情報処理システム１Ｂの構成図。
【図９】第４実施形態の処理部３Ｂの動作の流れの一例を示すフローチャート。
【図１０】第４実施形態におけるグループ数と予測誤差との関係を示すグラフ。
【図１１】コンピュータ５の構成図。
【発明を実施するための形態】
【０００８】
以下、実施形態の情報処理システムおよび情報処理プログラムを、図面を参照して説明する。なお、本願でいう「ＸＸに基づく」とは、少なくともＸＸに基づくことを意味し、ＸＸに加えて別の要素に基づく場合も含む。また「ＸＸに基づく」とは、ＸＸを直接に用いる場合に限定されず、ＸＸに対して演算や加工が行われたものに基づく場合も含む。「ＸＸ」は、任意の要素（例えば、任意の指標や物理量、その他の情報）である。
【０００９】
（第１実施形態）
図１は、実施形態の情報処理システム１の構成図である。情報処理システム１は、例えば、通信部２と、処理部３と、記憶部４とを備える。なお、情報処理システム１の各構成部は、１つの装置内に存在するものであってもよく、複数の装置に分かれて存在する構成であってもよい。
【００１０】
通信部２は、ネットワークを介して、各需要家のスマートメーターや気象データ等を保持するサーバ、またはそれらから情報を取得したサーバ等と通信する。通信部２は、スマートメーターやサーバなどから送信された情報を受信し、受信した情報を処理部３に出力する。また、通信部２は、処理部３の指示に従って、スマートメーターなどに情報を送信してもよい。
【００１１】
処理部３は、例えば、データ管理部３１と、グループ分類部３２とを備える。処理部３の各機能部は、例えばＣＰＵ（Central Processing Unit）等のプロセッサがプログラムメモリに格納されたプログラムを実行することにより実現される。また、処理部３の機能部のうち一部または全部は、ＬＳＩ（Large Scale Integration）、ＡＳＩＣ（Application Specific Integrated Circuit）、またはＦＰＧＡ（Field-Programmable Gate Array）等のハードウェアにより実現されてもよい。
【００１２】
データ管理部３１は、例えば、スマートメーターから受信した個別需要データ、及び、サーバから受信した気象データ等を、記憶部４に格納する。個別需要データとは、各需要家の過去の消費電力カーブの集合である。個別需要データは、例えば、スマートメーターにおいて計測された電力需要と、計測された日時を示す日時情報とを含む。気象データとは、過去の天候に関する情報であって、例えば、電力需要の計測日時における気温、湿度、降水量、風速等を含む。なお、気象データは、これに限られず、予測される将来の天候に関する情報であってもよい。
【００１３】
グループ分類部３２は、需要家ごとの個別需要データと気象データとに基づいて個別予測モデルを導出し、導出した個別予測モデルの相互の類似度に基づき複数の需要家を複数のグループに分類する。次いで、グループ分類部３２は、分類したグループごとに、各グループに分類された個別需要データを合成して、一つ以上のグループ需要データを生成する。ここで、グループ需要データとは、同一のグループに分類された複数の個別需要データを時間ごとに加算して生成された仮想的な消費電力カーブである。グループ分類部３２は、これらの処理を実行するため、個別予測モデル導出部３２１と、グループ生成部３２２と、個別データ合成部３２３とを備える。
【００１４】
個別予測モデル導出部３２１は、個別需要データと、気象データとに基づいて、需要家ごとの電力需要を予測する個別予測モデルを導出し、グループ生成部３２２に出力する。個別予測モデル導出部３２１は、例えば、回帰モデルに基づき、被説明変数（目的変数）に電力需要、説明変数に過去の電力需要に関する変数および気象に関数する変数を用いて、個別予測モデルを導出する。なお、説明変数はこれらに限られず、他の変数を用いてもよい。
【００１５】
グループ生成部３２２は、個別予測モデル導出部３２１によって導出された複数の需要家の個別予測モデルの相互の類似度を導出し、導出した類似度に基づき、複数の需要家を複数のグループに分ける。類似度とは、個別予測モデル同士がどれだけ似ているかを示す数値であってもよく、どれだけ似ていないかを示す数値（つまり、距離）であってもよい。グループ生成部３２２は、導出した類似度を記憶部４内のテーブルに書き込む。
【００１６】
個別データ合成部３２３は、グループ生成部３２２によって同一のグループに分けられた一つ以上の個別需要データを合成してグループ需要データを生成する。個別データ合成部３２３は、生成したグループ需要データを記憶部４に格納してもよく、通信部２を用いて外部装置に送信してもよい。
【００１７】
記憶部４は、気象データＤＢ４１と、個別需要データＤＢ４２と、類似度テーブル４３と、グループ需要データＤＢ４４とが記憶されている。記憶部４は、例えば、ＨＤＤ（Hard Disc Drive）、フラッシュメモリ、ＥＥＰＲＯＭ（Electrically Erasable Programmable Read Only Memory）、ＲＯＭ（Read Only Memory）、またはＲＡＭ（Random Access Memory）等により実現される。また、記憶部４には、上記の各種情報の他、アプリケーションプログラムなどの各種プログラム、処理部３による処理結果などが記憶される。
【００１８】
気象データＤＢ４１は、例えば、過去の天候に関する情報であって、気温、湿度、降水量、風速等と日時情報とが対応付けられた情報である。なお、気象データＤＢ４１は、これに限られず、予測される将来の天候に関する情報（予報値）を含んでもよい。
【００１９】
個別需要データＤＢ４２は、各需要家の個別需要データと、日時情報と、需要家を示す情報とが対応付けられた情報である。
【００２０】
類似度テーブル４３は、一の需要家の個別予測モデルに対する他の需要家の個別予測モデルの類似度を対応付けた情報である。
【００２１】
グループ需要データＤＢ４４は、処理部３により導出されたグループ需要データと、日時情報と、需要家を示す情報とが対応付けられた情報である。
【００２２】
次に、図２を参照して、個別需要データの一例について説明する。図２は、個別需要データの一例を示す図である。図２に示すグラフで、横軸は時間、縦軸は電力需要である。個別需要データｄ１（Ａ）〜ｄ１（Ｄ）は、需要家Ａ〜Ｄの過去の個別需要データである。個別需要データは、３０分を１単位で電力需要を集計したデータであって、各時刻における電力需要を結んだ線で示される。図示の例では、需要家Ａと需要家Ｃの電力波形において三日おきに同じ周期が現れ、需要家Ｂと需要家Ｄの電力波形において二日おきに同じ周期が現れている。なお、図示は省略するが、気象データも、個別需要データと同様に、３０分を１単位で、気温や湿度、降水量や風力等の気象データが用意されている。
【００２３】
次に、図３を参照して、グループ分類部３２による処理の一例について説明する。図３は、グループ分類部３２の動作の流れの一例を示すフローチャートである。ここでは、個別需要データｄ１（Ａ）〜ｄ１（Ｄ）を用いた処理の具体例について説明する。
【００２４】
（第一のステップ）
まず、個別予測モデル導出部３２１は、記憶部４から個別需要データｄ１（Ａ）〜ｄ１（Ｄ）と、気象データとを読み出し、読み出した情報に基づいて、それぞれの個別需要に対応する個別予測モデルを導出する（ステップＳ１）。ここでは、以下の通りの個別需要予測モデルｍ１（Ａ）〜ｍ１（Ｄ）が導出されたものとする。ここで、個別予測モデル導出部３２１は、回帰モデルに基づき、被説明変数に電力需要、説明変数に過去の電力需要に関する変数および気象に関数する変数を用いて、個別予測モデルを導出する。
【００２５】
ｍ１（Ａ）：demandA(t)=α1 * demandA(t-48*3) +α2 * temp(t) +α3…式（１）
ｍ１（Ｂ）：demandB(t)=β1 * demandB(t-48*2) +β2 * temp(t) +β3…式（２）
ｍ１（Ｃ）：demandC(t)=γ1 * demandC(t-48*3) +γ2 * temp(t) +γ3…式（３）
ｍ１（Ｄ）：demandD(t)=δ1 * demandD(t-48*2) +δ2 * temp(t) +δ3…式（４）
【００２６】
また、個別予測モデル導出部３２１は、以下の通り、各係数を計算したものとする。
α1=0.9, α2=0.2, α3=0.1 
β1=0.8, β2=0.2, β3=2.1
γ1=0.7, γ2=0.4, γ3=3.1
δ1=0.2, δ2=0.7, δ3=5.1
【００２７】
式（１）〜式（４）において、demand(t)は、将来のある日時tの電力需要である。temp(t)は、将来のある日時tの気温を表す。また、demand(t-48*3)は、時刻tより48*3*30分前の電力需要、すなわち、三日前の同時刻の電力需要を表すものとする。
【００２８】
（第二のステップ）
次に、グループ生成部３２２は、各需要家の類似度を導出する（ステップＳ２）。類似度は、さまざまなものが利用可能である。例えば、グループ生成部３２２は、個別予測モデルに含まれる変数の種類に基づき、類似度を導出する。ここでは、個別予測モデルに含まれる同一の変数の個数を類似度とする。式（１）〜（４）では、(t-48*3)とtemp(t)とが共通である需要家Ａと需要家Ｃが類似度２であり、(t-48*2)とtemp(t)とが共通である需要家Ｂと需要家Ｄが類似度２であり、temp(t)だけが共通である需要家Ａと需要家Ｂまたは需要家Ｄ、および需要家Ｃと需要家Ｂまたは需要家Ｄは、類似度１である。このように導出された類似度を図４に示す。図４は、類似度テーブルの一例を示す図である。
【００２９】
グループ生成部３２２は、導出した類似度に基づいて、類似している需要家同士を同一のグループに分類し、分類結果を個別データ合成部３２３に出力する（ステップＳ３）。ここで、グループ生成部３２２は、各需要家について類似度が最も大きい需要家を同一グループに分類する。図４に示す例では、需要家Ａについて類似度が最も大きい需要家はＣであり、需要家Ｂについて類似度が最も大きい需要家はＤである。よって、グループ生成部３２２は、需要家Ａと需要家Ｃとが同一のグループであり、需要家Ｂと需要家Ｄとが同一のグループであると判定し、判定結果を個別データ合成部３２３に出力する。なお、グループ生成部３２２は、これに限られず、類似度が閾値以上となる需要家同士（または類似度がある閾値範囲に含まれる需要家同士）を同一のグループに分類してもよく、類似度が大きい方から所定数の需要家同士を同一のグループに分類してもよい。前者の場合、閾値は、需要家の総数や類似度の中の最大値等に基づいて決定されてもよく、分類するグループ数に応じて決定されてもよい。後者の場合、所定数は、類似度の分布に応じて決定されてもよく、需要家の総数や分類するグループ数に応じて決定されてもよい。
【００３０】
（第三のステップ）
次いで、個別データ合成部３２３は、グループ生成部３２２から入力された情報に基づいて、個別需要データｄ１（Ａ）とｄ１（Ｃ）、個別需要データｄ１（Ｂ）とｄ１（Ｄ）とをそれぞれ加算し、グループ需要データｄ２（ＡＣ）とｄ２（ＢＤ）を生成する（ステップＳ４）。以下に、生成されたグループ需要データの数式を示す。
グループ需要データｄ２（ＡＣ）＝ｄ１（Ａ）＋ｄ１（Ｃ）…式（５）
グループ需要データｄ２（ＢＤ）＝ｄ１（Ｂ）＋ｄ１（Ｄ）…式（６）
【００３１】
上述の通り、各需要家に対応する個別予測モデルの類似度に基づいてグループ分けをすることにより、電力需要の傾向が類似する需要家同士をまとめることができる。このようにして分けられたグループごとに、電力需要を予測することで、予測精度を向上することができる。なお、この理由は、第２の実施形態のなかで詳しく述べる。
【００３２】
なお、グループ分類部３２による処理は、上述の内容に限られない。上述の例では、個別予測モデルの説明変数を２つとしたが、これに限られず、２以上の多数であってもよい。この場合、グループ生成部３２２は、説明変数が多数のまま類似度を導出してもよく、あるいは、複数の需要家の個別予測モデルに対して変数選択を実施し、個別予測モデルに含まれる変数の数を減らした状態で類似度を導出してもよい。
【００３３】
また、グループ生成部３２２は、個別予測モデルに含まれる変数の重みに基づき、類似度を導出してもよい。すなわち、類似度の選定基準は、共通変数に限られず、共通変数に加えて係数（重み）まで含めて類似度を算出してもよい。例えば、グル―プ生成部３２２は、需要家をグループに分ける際、類似度が一定基準を満たさない場合は、同一のグループと判断しない方法を採用してもよい。これにより、共通変数の係数の類似性を考慮することができる。
【００３４】
例えば、共通変数の係数であるα1とγ1、β1とδ1において、α1とγ1の差は0.2であり、β1とδ1の差は0.6である。ここで、一定基準を、変数の差分が第１閾値（例えば、０．５）以下であるとする。この場合、グル―プ生成部３２２は、差分が第１閾値より大きい共通変数は共通とみなさいため、需要家Bと需要家Dの類似度を1と導出する。よって、グル―プ生成部３２２は、需要家Ａと需要家Ｃとが一のグループであり、需要家Ｂが一のグループであり、需要家Ｄが一のグループであると判定する。
【００３５】
また、グル―プ生成部３２２は、類似度を算出する際、共通変数の差分に基づいて類似度を決定してもよい。例えば、グル―プ生成部３２２は、共通変数の差分がない場合に類似度＝１とし、共通変数の差分が第２閾値（例えば、０，５）以下である場合に類似度＝０．５とし、共通変数の差分が第２閾値より大きい場合に類似度＝０としてもよい。
【００３６】
また、第一のステップにおいて、個別予測モデル導出部３２１は、線形回帰モデルを用いて予測モデルを導出したが、変数の数を絞り込むために、AICのようなモデル選定指標を使ってAIC最小のモデルを導出してもよい。あるいは、個別予測モデル導出部３２１は、LASSO回帰のように、変数選択と回帰モデル決定を同時に実行可能な予測モデルを導出してもよい。
【００３７】
（第２実施形態）
図５は、第２実施形態の情報処理システム１Ａの構成図である。情報処理システム１Ａは、処理部３Ａと記憶部４Ａとを備える点で、情報処理システム１と異なる。処理部３Ａは、主に、全体需要予測部３３を備える点で処理部３と異なる。記憶部４Ａは、主に、グループ需要データＤＢ４４がなく、全体需要予測結果ＤＢ４５を備える点で記憶部４と異なる。
【００３８】
全体需要予測部３３は、複数のグループ需要データに基づいて、グループ分けした全ての需要家が存在する対象地域における電力需要の予測値である全体需要予測結果を生成する。全体需要予測部３３は、これらの処理を実行するため、グループ予測モデル導出部３３１と、グループ需要予測部３３２と、グループデータ合成部３３３とを備える。
【００３９】
グループ予測モデル導出部３３１は、記憶部４から、グループ需要データと気象データとを読み出し、読み出した情報に基づいてグループ需要予測モデルを導出する。グループ需要予測モデルとは、グループごとの電力需要を予測するモデルであって、回帰モデルを利用した予測モデルである。グループ予測モデル導出部３３１は、被説明変数に電力需要、説明変数に過去の電力需要に関する変数および気象に関数する変数を用いて、グループ需要予測モデルを導出する。なお、説明変数はこれらに限られず、他の変数を用いてもよい。なお、グループ予測モデル導出部３３１は、グループ需要データのようにグループごとに個別需要データが合成されたデータに限られず、合成される前の個別需要データを取得し、取得した情報をグループごとに合成して用いてもよい。
【００４０】
グループ需要予測部３３２は、グループ予測モデル導出部３３１により導出されたグループ需要予測モデルを用いて、グループ需要データと気象データからグループごとの需要予測結果を計算する。ここで、気象データは、予測する将来の日時における予報の情報である。
【００４１】
グループデータ合成部３３３は、グループごとの需要予測結果を時間ごとに合計して、全体需要予測結果を計算する。
【００４２】
記憶部４Ａは、気象データＤＢ４１と、個別需要データＤＢ４２と、類似度テーブル４３と、全体需要予測結果４５とを備える。全体需要予測結果４５は、処理部３Ａにより導出された全体需要予測結果に、日時情報と、対象地域を示す情報とを対応付けた情報である。
【００４３】
次に、図６を参照して、全体需要予測部３３による処理の一例について説明する。全体需要予測部３３の動作の流れの一例を示すフローチャートである。ここでは、グループ分類部の個別データ合成部３２３から、全体需要予測部３３に、例えばグループ需要データｄ２（ＡＣ）とｄ２（ＢＤ）が入力されるものとする。
【００４４】
（第四のステップ）
まず、グループ予測モデル導出部３３１は、入力するグループ需要データｄ２（ＡＣ）とｄ２（ＢＤ）と、記憶部４から読み出した気象データとに基づいて、各グループに対応するグループ予測モデルを導出する（ステップＳ５）。ここでは、以下に示すグループ需要予測モデルｍ２（Ａ）〜ｍ２（Ｄ）が生成されたものとする。
【００４５】
ｍ２（ＡＣ）：demandAC(t)=ε1 * demandAC(t-48*3) +ε2 * temp(t) +ε3…式（７）
ｍ２（ＢＤ）：demandBD(t)=ζ1 * demandBD(t-48*2) +ζ2 * temp(t) +ζ3…式（８）
【００４６】
また、グループ予測モデル導出部３３１は、以下に示す通り、各係数を計算したものとする。
ε1=0.6, ε2=0.3, ε3=1.5
ζ1=0.5, ζ2=0.5, ζ3=4.1
【００４７】
（第五のステップ）
次に、グループ需要予測部３３２は、グループ予測モデル導出部３３１により導出されたグループ需要予測モデルｍ２（ＡＣ）とｍ２（ＢＤ）とを利用して、設定された予測日の需要予測を行い、需要予測結果をグループデータ合成部３３３に出力する（ステップＳ６）。グループ需要予測部３３２は、例えば、翌日の２４時間分の需要予測結果であるdemandACとdemandBDとを導出する。
【００４８】
（第六のステップ）
最後に、グループデータ合成部３３３は、グループごとに導出されたグループごとの需要予測結果demandACとdemandBDとを加算して、加算結果である全体需要予測結果を、記憶部４に保存する。
【００４９】
上述の通り、個別需要データを類似度に基づいてグループに分けた後、グループごとに電力需要を予測し、予測結果に基づいて全体需要予測結果を得ることにより、以下に説明する第１の方法および第２の方法に比べて、精度の高い総需要予測を実現可能なシステムを得ることができる。これは、以下に説明する問題を解消することができるからである。つまり、個別需要データごとに予測結果を生成し、これらを合計することで全体需要予測結果を生成する第１の方法では、特に個別需要データの数が少ない場合、個別需要データの偏りで、個別需要予測モデルが真のモデルから乖離する可能性があり、その合計値も真の総需要を予測するモデルから乖離する可能性がある。さらに、ノイズが大きくなる可能性もある。また、個別需要データを同一時刻で合計して、総需要データを生成し、このデータから予測結果を生成する第２の方法では、総需要データを生成した時点で、個別需要の特性が失われ、真の総需要を予測するモデルから乖離する可能性がある。一方で、第２実施形態に係る情報処理システム１Ａによれば、第１の方法における個別需要データの偏りの問題と、第２の方法におけるモデルに個別需要の特性を持たせる問題とを解消することにより、精度の高い総需要予測を実現することができる。
【００５０】
なお、第四のステップでは、第一のステップと同様に、線形回帰モデルを用いて予測モデルを導出したが、モデルの予測精度向上のために、AICを使った変数選択や、LASSO回帰を使った変数選択と回帰モデルの同時決定を実施してもよい。
【００５１】
（第３実施形態）
第３実施形態は、第２実施形態と同様の構成であるため、構成についての詳細な説明は省略するが、グループ分類部３２と全体需要予測部３３における処理の一部が以下の点で異なる。
【００５２】
第一のステップにおいて、個別予測モデル導出部３２１は、第２実施形態と同様に式（１）〜式（４）に示す個別需要予測モデルｍ１（Ａ）〜ｍ１（Ｄ）を導出する。なお、係数は以下の通りである。
α1=0.9, α2=0.2, α3=0.1, 
β1=0.8, β2=0.2, β3=2.1,
γ1=0.7, γ2=0.4, γ3=3.1,
δ1=0.2, δ2=0.7, δ3=5.1,
【００５３】
第二のステップは、第２実施形態と同様であり、グループ生成部３２２は、例えば、需要家Ａと需要家Ｃとが同一のグループであり、需要家Ｂと需要家Ｄとが同一のグループであると判定したとする。
【００５４】
第三のステップにおいて、個別データ合成部３２３は、上述の式（５），（６）に示した通り、同じグループに分類された個別需要データを単純に合計するグループ需要データｄ２（ＡＣ）とｄ２（ＢＤ）に加え、以下に式（９），（１０）に示すグループ需要データｄ２´（ＡＣ）とｄ２´（ＢＤ）を導出する。つまり、個別データ合成部３２３は、個別予測モデル導出部３２１により導出された個別予測モデル内の係数を用いて個別需要データを重み付けした値を合成したグループ需要データｄ２´（ＡＣ）とｄ２´（ＢＤ）も導出する。
グループ需要データｄ２´（ＡＣ）＝α1 *ｄ１（Ａ）＋γ1 *ｄ１（Ｃ）…式（９）
グループ需要データｄ２´（ＢＤ）＝β1 *ｄ１（Ｂ）＋δ1 *ｄ１（Ｄ）…式（１０）
【００５５】
第四のステップにおいて、グループ予測モデル導出部３３１は、個別予測モデルに含まれる係数によって重み付けされた需要家ごとの個別需要データがグループごとに纏められたデータに基づき、グループごとの電力需要の予測モデルを導出する。つまり、グループ予測モデル導出部３３１は、入力するグループ需要データｄ２（ＡＣ），ｄ２（ＢＤ），ｄ２´（ＡＣ），ｄ２´（ＢＤ）と、記憶部４から読み出した気象データとに基づいて、各グループに対応するグループ予測モデルを導出する。ここでは、以下の通りのグループ需要予測モデルｍ２（Ａ）〜ｍ２（Ｄ）が生成されたものとする。
【００５６】
ｍ２（ＡＣ）：demandAC(t)=ε1 * demandAC’(t-48*3) +ε2 * temp(t) +ε3…式（１１）
ｍ２（ＢＤ）：demandBD(t)=ζ1 * demandBD’(t-48*2) +ζ2 * temp(t) +ζ3…式（１２）
【００５７】
また、グループ予測モデル導出部３３１は、以下の通り、各係数を計算したものとする。
ε1=1.1, ε2=0.3, ε3=1.5,
ζ1=0.9, ζ2=0.5, ζ3=4.1,
【００５８】
なお、式（１１），（１２）の説明変数のうち、電力需要に関連するものは、重みを加えた個別需要データであるdemandAC’(t)およびdemandBD’(t)により定義されている。なお、demandAC’(t)は、式（９）に示すグループ需要データｄ２´（ＡＣ）であり、demandBD’(t)は、式（１０）に示すグループ需要データｄ２´（ＢＤ）である。一方、式（７），（８）に示すグループ需要データｄ２（ＡＣ），ｄ２（ＢＤ）は、グループ予測モデル導出部３３１において、係数を導出する際に利用される。
【００５９】
第五〜六のステップは、第２実施形態と同様である。つまり、グループ需要予測部３３２は、グループ予測モデル導出部３３１により導出されたグループ需要予測モデルｍ２（ＡＣ）とｍ２（ＢＤ）とを利用して、設定された予測日の需要予測を行い、需要予測結果をグループデータ合成部３３３に出力する。グループデータ合成部３３３は、グループごとに導出されたグループごとの需要予測結果demandACとdemandBDとを加算して、加算結果である全体需要予測結果を、記憶部４に保存する。
【００６０】
上述の通り、個別需要データに重み付けをして得られたグループ需要データを用いてグループ需要予測モデルを導出することにより、総需要の予測精度の向上が期待できる。すなわち、個別予測モデルを生成した際に、需要予測に過去の需要が影響を与える度合いが回帰モデルの係数（重み）として与えられている。ここで、合計需要を単純に合計してグループ需要予測モデルを生成した場合、個別の需要家の過去の需要の影響度が同程度に影響を与えることになる。一方で、本実施形態のように、個別需要データに重み付けをしてグループ需要予測モデルを生成することで、個別の需要家の過去の需要の影響度を、もともとの個別需要予測モデルの算出で計算された量と一致させることができるためである。これにより、個別の需要家の過去の需要の影響度をグループ需要予測モデルに反映させることができ、総需要の予測精度の向上を期待することができる。
【００６１】
（変形例）
次に、第３実施形態の変形例について説明する。ここで、グループ生成部３２２は、複数の需要家の個別予測モデルに含まれる変数の間に時間差成分がある場合、時間差成分にも基づき類似度を導出する。以下、具体的に説明する。
【００６２】
第一のステップにおいて、個別予測モデル導出部３２１は、式（１）〜（４）に替えて、下の式（１３）〜（１６）に示す個別予測モデルｍ２（Ａ）〜ｍ２（Ｄ）を導出したものとする。
【００６３】
ｍ２（Ａ）：demandA(t)=α1 * demandA(t-48*3) +α2 * temp(t) +α3…式（１３）
ｍ２（Ｂ）：demandB(t)=β1 * demandB(t-48*2) +β2 * temp(t) +β3…式（１４）
ｍ２（Ｃ）：demandC(t)=γ1 * demandC(t-48*3+1) +γ2 * temp(t) +γ3…式（１５）
ｍ２（Ｄ）：demandD(t)=δ1 * demandD(t-48*2+1) +δ2 * temp(t) +δ3…式（１６）
【００６４】
また、個別予測モデル導出部３２１は、以下の通り、各係数を計算したものとする。
α1=0.9, α2=0.2, α3=0.1
β1=0.8, β2=0.2, β3=2.1
γ1=0.7, γ2=0.4, γ3=3.1
δ1=0.2, δ2=0.7, δ3=5.1
【００６５】
式（１５），（１６）に示す通り、需要家Ｃと需要家Ｄの個別予測モデルにおいて電力需要の変数の時間が１単位（つまり３０分）加算されている点で、上述の個別予測モデルと異なる。しかし、グループ生成部３２２は、時間が１単位ずれた電力需要の変数と時間の加算のない同じ種類の電力需要の変数との相互相関が高ければ、これらの変数を同一とみなす。つまり、グループ生成部３２２は、変数における時間周期にずれがあるが、そのずれ幅が閾値以下である場合、これらの変数を同一とみなすことができる。変数における時間周期にずれがある場合とは、例えば、ｔ−（４８×２−１）に対して、ｔ−４８×２は、１単位ずれていることを表している。
【００６６】
ここでは、個別予測モデルに含まれる同一の変数の個数を類似度とする。ただし、変数が完全に同一でない場合であっても、相関係数が類似度に加算される。例えば、需要家Ａと需要家Ｃとの類似度について説明する。グループ生成部３２２は、需要家Ａの過去の電力需要に関する変数［ｔ−４８×３］と需要家Ｃの過去の電力需要に関する変数［ｔ−（４８×３−１）］の相関係数を、０．９と算出する。また、需要家Ａと需要家Ｃとは、気象に関する変数［temp(t)］が一致している。よって、グループ生成部３２２は、需要家Ａと需要家Ｃとの類似度を、０．９＋１＝１．９と算出する。このように導出された類似度を図４に示す。図４は、類似度テーブルの一例を示す図である。このように導出された類似度を図７に示す。図７は、類似度テーブルの他の例を示す図である。
【００６７】
これに限られず、グループ生成部３２２は、系列の差を時刻差で重み付けしてもよい。例えば、グループ生成部３２２は、類似度を重みｗとして計算し、同じ時刻の系列の重みｗを1、時間が1単位（３０分）ずれるごとに１から0.1減算した値を重みｗと計算してもよい。
【００６８】
また、第三のステップにおいてグループ需要データを算出する際、個別データ合成部３２３は、同時刻の個別需要データを合成するのではなく、時刻のずれ、および回帰モデルの係数に基づいて、類似度の高い時間の個別需要データ同士を合成して、グループ需要データを導出してもよい。
【００６９】
（第４実施形態）
図８は、第４実施形態の情報処理システム１Ｂの構成図である。情報処理システム１Ｂは、処理部３Ｂと記憶部４Ｂとを備える点で、情報処理システム１Ａと異なる。処理部３Ｂは、予測結果比較部３４を備える点で処理部３Ａと異なる。記憶部４Ｂは、予測誤差ＤＢ４６と、グループ設定値４７とを備える点で記憶部４Ａと異なる。
【００７０】
予測結果比較部３４は、異なるグループ数でそれぞれ全体需要予測結果を算出し、それら全体需要予測結果と電力需要の実測値との予測誤差を算出し、算出された予測誤差に基づいて、グループ生成部３２２において設定されるグループ数（以下、グループ設定値と記す）を決定する。グループ設定値とは、上述した通りグループ需要データや全体需要予測結果を導出する際に、グループ生成部３２２による処理において、需要家を分類するグループ数として予め設定される値のことである。予測結果比較部３４は、これらの処理を実行するため、予測誤差導出部３４１と、グループ設定値決定部３４２とを備える。
【００７１】
予測誤差導出部３４１は、グループ分類部３２に対して、異なるグループ数でそれぞれグループ需要データを導出するよう指示するとともに、全体需要予測部３３に対して、グループ分類部３２において導出されたグループ需要データに基づいて、異なるグループ数のそれぞれに対応する全体需要予測結果を導出するよう指示する。ここで、予測誤差導出部３４１は、予測日時を過去の日時（例えば、昨日の２４時間）であることも合わせて指示する。これらの指示を、以下、グループ設定指示と記す。「異なるグループ数」において、最も大きいグループ数は需要家の数と同数であり、最も小さいグループ数は１である。グループ数が最も大きい数の場合、全ての需要家が個別に需要予測され、グループ数が最も小さい数の場合、すべての需要家が一つのグループにまとめられる。
【００７２】
グループ設定指示には、異なるグループ数の数だけ全体需要予測結果を導出する処理（以下、予測処理と記す）を繰り返すよう制御する指示も含まれる。また、グループ設定指示には、グループ数を大きい方から小さい方に変化させながら予測処理を繰りかえすよう制御する指示が含まれていてもよい。ただし、以下の説明では、予測誤差導出部３４１が異なるグループ数を決定し、グループ生成部３２２に指示する例について説明する。予測誤差導出部３４１は、例えば、１回目の予測処理において需要家の数をグループ数とし、２回目以降の予測処理において前回の予測処理で用いたグループ数を２で割った数をグループ数とする。
【００７３】
また、予測誤差導出部３４１は、全体需要予測部３３により導出された全体需要予測結果と、過去の個別需要データとの比較結果である予測誤差を、異なるグループ数のそれぞれについて導出し、導出した予測誤差をグループ設定値決定部３４２に出力する。予測誤差は、例えば、複数の需要家の過去の個別需要データの合計値と全体需要予測結果との差分である。
【００７４】
グループ設定値決定部３４２は、複数の予測誤差に基づいて、予め決められた評価基準に従って予測誤差を評価し、評価結果に基づいてグループ設定値（グループ数）を決定する。グループ設定値決定部３４２は、決定したグループ設定値を記憶部４に格納する。評価基準は、例えば、予測誤差が最小となる全体需要予測結果が最適であることが決められている。これに限られず、評価基準は、予測誤差が閾値以下となる全体需要予測結果が所定の条件を満たすと決められていてもよい。また、グループ設定値決定部３４２は、最適または所定の条件を満たすと判定された全体需要予測結果を導出した際に用いたグループ数をグループ設定値に決定してもよい。グループ設定値決定部３４２は、最適または所定の条件を満たすと判定した全体需要予測結果を記憶部４に格納してもよく、通信部２を介して送信してもよい。
【００７５】
グループ生成部３２２は、記憶部４を参照し、グループ設定値決定部３４２により決定されたグループ設定値が示す数のグループに複数の需要家を分類する。
【００７６】
予測誤差ＤＢ４６は、予測誤差導出部３４１により導出された予測誤差に、予測誤差が導出された際に用いられたグループ数が対応付けられた情報である。
【００７７】
グループ設定値４７は、グループ設定値決定部３４２により決定されたグループ設定値に、対象とした複数の需要家を含む対象地域を示す情報が対応付けられた情報である。
【００７８】
次に、図９を参照して、第４実施形態に係る処理の一例について説明する。図９は、第４実施形態に係る処理の流れの一例を示すフローチャートである。
【００７９】
まず、予測誤差導出部３４１が、グループ分類部３２と全体需要予測部３３とに対して、グループ設定指示を出力する（ステップＳ１１）。これにより、個別予測モデル導出部３２１は、記憶部４から複数の需要家の過去の個別需要データと、グループ設定指示により指定された過去の気象データとを読み出し、読み出した情報に基づいて、それぞれの個別需要に対応する個別予測モデルを導出する（ステップＳ１２）。
【００８０】
次いで、グループ生成部３２２は、各需要家の類似度を導出し、導出した類似度を類似度テーブル４３に書き込む（ステップＳ１３）。そして、グループ生成部３２２は、グループ設定指令に基づくグループ数で、各需要家を分類する（ステップＳ１４）。１回目の予測処理では、グループ数＝最大値であるため、各需要家は需要家の数に分類される。つまり、一つのグループに属する需要家の数は一つである。このように、グループ数＝最大値である場合、グループ生成部３２２は、分類という処理をしないものであってもよい。個別データ合成部３２３は、グループ生成部３２２による処理結果に基づいて、グループ需要データを生成し、グループ予測モデル導出部３３１に出力する（ステップＳ１５）。１回目の予測処理では、各需要家の個別需要データがグループ需要データとなる。
【００８１】
グループ予測モデル導出部３３１は、入力するグループ需要データと、記憶部４から読み出したグループ設定指示により指定された過去の気象データとに基づいて、各グループに対応するグループ予測モデルを導出し、グループ需要予測部３３２に出力する（ステップＳ１６）。グループ需要予測部３３２は、グループ予測モデル導出部３３１により導出されたグループ需要予測モデルを利用して、グループ設定指示により指定された複数の日時の需要予測を行い、需要予測結果をグループデータ合成部３３３に出力する（ステップＳ１７）。グループデータ合成部３３３は、グループごとに導出された各グループの需要予測結果を加算して、加算結果である全体需要予測結果を、予測誤差導出部３４１に出力する（ステップＳ１８）。そして、予測誤差導出部３４１は、全体需要予測部３３により導出された全体需要予測結果と、グループ設定指示により指定された日時の複数の需要家の個別需要データの合計値との比較結果である予測誤差を導出し、導出した予測誤差にグループ数を対応付けて、予測誤差ＤＢ４６に書き込む（ステップＳ１９）。
【００８２】
次いで、予測誤差導出部３４１は、グループ数＝１であるか否かを判定する（ステップＳ２０）。グループ数＝１でない場合、予測誤差導出部３４１は、前回のグループ数に基づいてグループ数を変更し、変更後のグループ数をグループ生成部３２２に出力する（ステップＳ２０）。ここで、予測誤差導出部３４１は、グループ数を前回のグループ数よりも少ないグループ数に変更する。例えば、予測誤差導出部３４１は、前回のグループ数の半分の数にグループ数を変更してもよい。そして、ステップＳ１４に戻って、グループ生成部３２２は、予測誤差導出部３４１により変更されたグループ数で、各需要家を分類する。
【００８３】
一方、ステップＳ２０において、グループ数＝１である場合、予測誤差導出部３４１は、その旨をグループ設定値決定部３４２に出力する。グループ設定値決定部３４２は、予測誤差ＤＢ４６から、複数の予測誤差を読み出し、予め決められた評価基準に従って予測誤差を評価し、評価結果に基づいてグループ設定値を決定する（ステップＳ２２）。
【００８４】
需要家Ａ〜Ｄを例に説明すると、１回目の予測処理において、需要家Ａ，Ｂ，Ｃ，Ｄの個別需要データに基づく全体需要予測とその予測誤差Ｅ１とが導出される。２回目の予測処理において、需要家Ａと需要家Ｃ、需要家Ｂと需要家Ｄとがグループに分けられたとする。そして、各グループのグループ需要データに基づく全体需要予測とその予測誤差Ｅ２とが導出される。次いで、３回目の予測処理では、需要家Ａ〜Ｄが同じグループに分けられる。そして、全てが一つのグループとするグループ需要データに基づく全体需要予測とその予測誤差Ｅ３とが導出される。このようにして導出された予測誤差Ｅ１〜Ｅ３の関係を図１０に示す。図１０は、グループ数と予測誤差との関係を示すグラフである。そして、グループ設定値決定部３４２は、図１０に示す予測誤差Ｅ１〜Ｅ３を比較して、予測誤差が最小となるグループ数「２」をグループ設定値に決定する。
【００８５】
上述の通り、異なるグループ数の複数の予測処理を行い、各予測処理における予測誤差同士を比較することにより、予測精度の高いグループ数を求めることができる。よって、以降の予測処理において、予測精度の高いグループ数で電力需要を予測することにより、予測誤差を削減し、予測精度を高めることができる。
【００８６】
なお、上述の実施形態における処理部３，３Ａ，３Ｂは、図１１を参照して説明するコンピュータ５上のソフトウェアとして構成されてもよい。図１１は、ＰＣの構成図の一例である。コンピュータ５は、ＣＰＵ５１、メモリ５２、ハードディスク５３、インタフェース（ＩＦ）５４、および、グラフィックインタフェース（ＧＵＩ）５５を備え、それぞれがバス５６を介して接続されている。処理部３，３Ａ，３Ｂの機能を実現するプログラムは、ハードディスク５３上に格納され、実行時にメモリ５２上へと展開された後、手順に従って実行される。中間データであるグループ需要データは、メモリ５２上で管理されてもよく、ハードディスク５３上に格納されてもよい。また、処理部３，３Ａ，３Ｂは、これらのデータをＧＵＩ５５を通じで運用者に可視的に出力してもよい。
【００８７】
以上説明した少なくともひとつの実施形態によれば、個別予測モデル導出部によって導出された複数の需要家の個別予測モデルの相互の類似度に基づき複数の需要家を複数のグループに分けるグループ生成部を持つことにより、電力需要の予測精度の向上を図ることができる。
【００８８】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００８９】
１…情報処理システム、２…通信部、３…処理部、４…記憶部、５…ＰＣ、３１…データ管理部、３２…グループ分類部、３３…全体需要予測部、３４…予測結果比較部、４１…気象データＤＢ、４２…個別需要データＤＢ、４３…類似度テーブル、４４…グループ需要データＤＢ、４５…全体需要予測結果ＤＢ、４６…予測誤差ＤＢ、４７…グループ設定値、５１…ＣＰＵ、５２…メモリ、５３…ハードディスク、５４…インタフェース、５５…ＧＵＩ、５６…バス、３２１…個別予測モデル導出部、３２２…グループ生成部、３２３…個別データ合成部、３３１…グループ予測モデル導出部、３３２…グループ需要予測部、３３３…グループデータ合成部、３４１…予測誤差導出部、３４２…グループ設定値決定部
【図１】
2018113817.tif 000003
【図２】
2018113817.tif 000004
【図３】
2018113817.tif 000005
【図４】
2018113817.tif 000006
【図５】
2018113817.tif 000007
【図６】
2018113817.tif 000008
【図７】
2018113817.tif 000009
【図８】
2018113817.tif 000010
【図９】
2018113817.tif 000011
【図１０】
2018113817.tif 000012
【図１１】
2018113817.tif 000013