(57)【特許請求の範囲】
【請求項１】
自然エネルギーを用いる発電設備を少なくとも一つ含む電力系統の管理を行う情報処理装置において、
過去の気象の推移を示す気象情報、および前記発電設備の過去の発電量の推移を示す発電量データを少なくとも記憶している記憶媒体としての情報保存部と、
前記情報保存部に記憶される過去の気象の推移を示す気象情報から、発電量の予測の対象となる日と気象の推移が類似する過去の日を少なくとも二日分選択する類似日選択手段と、
前記類似日選択手段により選択される日の発電量データを前記情報保存部から取得し、取得した前記発電量データを一定時間毎の時系列データとして生成する時系列データ生成手段と、
前記時系列データを用いて、前記類似日選択手段により選択される日の時刻間の発電量の相関を少なくとも示す数理モデルを生成する数理モデル生成手段と、
前記数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布を決定するとともに確率分布に従う乱数の組を複数生成することによって前記発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを複数生成するシナリオ生成手段と、
を具備することを特徴とする情報処理装置。
【請求項２】
自然エネルギーを用いる発電設備を少なくとも一つ含む電力系統の管理を行う情報処理装置において、
過去の気象の推移を示す気象情報、前記発電設備の過去の発電量の推移を示す発電量データ、および需要設備の過去の需要量の推移を示す需要量データを少なくとも記憶している情報保存部と、
前記情報保存部に記憶される過去の気象の推移を示す気象情報から、発電量および需要量の予測の対象となる日と気象の推移が類似する過去の日を少なくとも二日分選択する類似日選択手段と、
前記類似日選択手段により選択される日の発電量データおよび需要量データを前記情報保存部から取得し、取得した前記発電量データおよび前記需要量データを一定時間毎の時系列データとして生成する時系列データ生成手段と、
前記時系列データを用いて、前記類似日選択手段により選択される日の時刻間の発電量の相関を示すと共に時刻間の需要量の相関を示し、かつ、発電量と需要量との相関をも示す数理モデルを生成する数理モデル生成手段と、
前記数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布と、一定時間毎の需要量の確率分布とを決定するとともに、これらの確率分布にそれぞれ従う乱数の組を複数生成することによって、前記発電設備の一定時間毎の発電量の予測値を示す発電量シナリオおよび前記需要設備の一定時間毎の需要量の予測値を示す需要量シナリオを複数生成し、前記需要量シナリオから前記発電量シナリオを差し引くことによって前記発電設備以外の設備に求められる一定時間毎の必要発電量の予測値を示す必要発電量シナリオを複数生成するシナリオ生成手段と、
を具備することを特徴とする情報処理装置。
【請求項３】
前記数理モデル生成手段により生成される数理モデルは、自然エネルギーを用いる複数の発電設備のぞれぞれの発電量の相関をも示すものであることを特徴とする請求項１又は２に記載の情報処理装置。
【請求項４】
前記確率密度関数は、多次元正規分布を形成するものであることを特徴とする請求項１乃至３のいずれか１項に記載の情報処理装置。
【請求項５】
前記確率密度関数は、ベータ分布を形成するものであることを特徴とする請求項１乃至３のいずれか１項に記載の情報処理装置。
【請求項６】
前記シナリオ生成手段により生成される必要発電量シナリオを用いて、各発電設備の起動，停止，および出力増減のタイミングを一定期間に渡って定めた発電計画を生成する発電計画生成手段を更に具備することを特徴とする請求項２乃至４のいずれか１項に記載の情報処理装置。
【請求項７】
前記発電計画生成手段は、前記シナリオ生成手段により生成される必要発電量シナリオの平均値の変動に応じて電力系統に含まれる蓄電池を充放電させる制御を含む発電計画を生成することを特徴とする請求項６に記載の情報処理装置。
【請求項８】
前記シナリオ生成手段により生成される必要発電量シナリオに応じて電力系統に含まれる蓄電池の容量を決定する容量決定手段を更に具備することを特徴とする請求項６に記載の情報処理装置。
【請求項９】
前記発電計画生成手段により生成される発電計画に基づき、発電量の期待値、発電コストの期待値、および収益の期待値の少なくとも一つを求める期待値算出手段を更に具備することを特徴とする請求項６または７に記載の情報処理装置。
【請求項１０】
自然エネルギーを用いる発電設備を少なくとも一つ含む電力系統の管理を行うコンピュータに、
記憶媒体に記憶される過去の気象の推移を示す気象情報から、発電量の予測の対象となる日と気象の推移が類似する過去の日を少なくとも二日分選択する類似日選択機能と、
前記類似日選択機能により選択される日の発電量の推移を示す発電量データを前記記憶媒体から取得し、取得した前記発電量データを一定時間毎の時系列データとして生成する時系列データ生成機能と、
前記時系列データを用いて、前記類似日選択機能により選択される日の時刻間の発電量の相関を少なくとも示す数理モデルを生成する数理モデル生成機能と、
前記数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布を決定するとともに確率分布に従う乱数の組を複数生成することによって前記発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを複数生成するシナリオ生成機能と、
を実現させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自然エネルギー（又は再生可能エネルギーとも称す）を用いる発電設備（例えば、太陽光発電設備，風力発電設備等）を少なくとも一つ含む電力系統において、自然エネルギーを用いる発電設備の発電量を予測することが可能な情報処理装置に関する。
【背景技術】
【０００２】
自然エネルギーを用いる発電設備は、燃料費がかからず、且つ温室効果ガスの排出が少ないという利点があるが、一方で、気象状態に大きく依存するため、発電量が不確実であり、且つ予測が困難であるという難点がある。
【０００３】
関連する代表的な技術としては以下に示す、特許文献１〜特許文献３が挙げられる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−２８９９１８号公報
【特許文献２】特開２００７−２２８６７６号公報
【特許文献３】特開２００４−７２９００号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１では、自然エネルギーを用いる発電設備の発電量を、気象情報に基づいて予測し、不足分を他の発電設備の発電量により補うとしているが、予測に変動幅があることが考慮されていない。
【０００６】
また、発電量と需要量のバランスが保たれないと、周波数が変動するなど電力系統が不安定になる。このため、自然エネルギーを導入した場合には、電力系統の管理を行う上で特別の注意が必要である。
【０００７】
一つ又は複数の発電設備の起動、停止、および出力の増減のタイミングを一定期間に渡って定めたものを発電計画と称す。例えば、発電設備が一つであるときは、その発電設備は電力需要が存在する間は常に運転され、その出力は電力需要に合わせて調整すればよい。また、発電設備が複数あるときは、どの発電設備をいつ用いるかを決定する必要がある。通常は、発電コストの安い発電設備から優先的に使用する。また、火力発電等の発電設備を用いる場合、発電設備の運転中に燃料費がかかるだけでなく、発電設備を起動するときにも起動のためのコストがかかる。従って、発電設備が複数あるときには、どの発電設備をいつ用いるかを決定することは容易ではなく、複雑な数理計画問題が生じる。
【０００８】
大規模であれ、小規模であれ、複数の発電設備を有する電力系統を運用する場合、少なくとも前日には、翌日の発電計画を生成する必要がある。この場合、例えば、翌日の電力需要の予測結果に基づいて発電設備の最適な発電計画（又は、起動停止計画）を生成する。このためには、例えば、動的計画法、ラグランジュ緩和法、および２次計画法等の数理計画法が用いられる。このような分野における先行技術としては、特許文献２が挙げられる。
【０００９】
しかしながら、特許文献２は需要の変動は考慮されているが、個々の発電計画の生成においては、決まった需要変動に対して発電計画を生成している。このため、電力系統に自然エネルギーを用いる発電設備を有する場合、需給バランスを保つために必要な発電量が決定しないため、翌日の発電計画を生成することは困難である。
【００１０】
これに対する最も簡単な解決策としては、自然エネルギーを用いる発電設備を負の需要とみなすことが挙げられる。自然エネルギーを用いる発電設備が１台であり、需要が変動しないとみなせる場合には、この解決策でもよい。しかしながら、自然エネルギーを用いる発電設備が複数台ある又は需要変動がある場合には、全ての変動の組み合わせを考慮する必要があるため、この解決策を利用することは困難である。また、自然エネルギーを用いる発電設備の発電量が需要量を上回る場合にも、発電量を制御可能な発電設備の発電量が負となるため、この方法を利用することができない。
【００１１】
このように、自然エネルギーを用いる発電設備の発電量を需要変動も含めて予測し、最も適当な発電計画を生成することは従来の方法では困難であった。
【００１２】
なお、自然エネルギーを用いる発電設備を用いる場合の他の解決策としては、容量が十分に大きい蓄電池を利用し、その蓄電池の出力を制御することで、蓄電池の充放電量と自然エネルギーを用いる発電設備の発電量との合計量を一定に保つ方法が挙げられる。このような分野における、先行技術としては、特許文献３が挙げられる。
【００１３】
特許文献３では、他の電力系統からの電力の流出入（潮流）を用いており、自らの発電設備は用いていない。すなわち、必要に応じて自由に電力を調達し、送電することで需給バランスを保っているため、発電計画の生成は他者に任されている。
【００１４】
自然エネルギーを用いる発電設備の発電量は予測が困難であるため、前日の予測は諦めて、蓄電池と発電当日の制御技術とによって需要変動を最小にするという考え方もある。しかしながら、基本的には予測が困難であるとしても、可能な限り前日に発電量の予測を行い最善の発電計画を生成しておく必要がある。
【００１５】
また、従来の技術では、例えば、太陽光発電の発電量を予測する場合には、基本的には全天日射量を予測し、複雑な天文計算および幾何学的計算を行い、太陽光発電パネルの効率を評価することにより発電量を予測しなければならなかった。同様に、風力発電の発電量を予測する場合においても、気象予報の風況から発電量を予測しなければならなかった。
【００１６】
以上に示すように、自然エネルギーを用いる発電設備を含む電力系統において、前日等に発電量の予測を行う場合、誤差を考慮した予測を行うことは困難である。
【００１７】
本発明は上記実情に鑑みてなされたものであり、自然エネルギーを用いる発電設備からの発電量を容易に予測することが可能な情報処理装置を提供することを目的とする。
【課題を解決するための手段】
【００１８】
本発明の一態様による情報処理装置は、自然エネルギーを用いる発電設備を少なくとも一つ含む電力系統の管理を行う情報処理装置において、過去の気象の推移を示す気象情報、および前記発電設備の過去の発電量の推移を示す発電量データを少なくとも記憶している記憶媒体としての情報保存部と、前記情報保存部に記憶される過去の気象の推移を示す気象情報から、発電量の予測の対象となる日と気象の推移が類似する過去の日を少なくとも二日分選択する類似日選択手段と、前記類似日選択手段により選択される日の発電量データを前記情報保存部から取得し、取得した前記発電量データを一定時間毎の時系列データとして生成する時系列データ生成手段と、前記時系列データを用いて、前記類似日選択手段により選択される日の時刻間の発電量の相関を少なくとも示す数理モデルを生成する数理モデル生成手段と、前記数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布を決定するとともに確率分布に従う乱数の組を複数生成することによって前記発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを複数生成するシナリオ生成手段とを具備することを特徴とする。
【発明の効果】
【００１９】
本発明によれば、自然エネルギーを用いる発電設備からの発電量を容易に予測することが可能な情報処理装置を提供することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の各実施形態に共通する電力システムのハードウェア構成の一例を示す図。
【図２】本発明の第１の実施形態に係る情報処理装置の機能の一例を示す機能ブロック図。
【図３】分類アルゴリズムを用いて、過去の気象の推移を示す気象情報を天候に応じて分類する方法の一例を示す図。
【図４】ある年の１２時における東京の全天日射量の分布の一例を示すグラフ。
【図５】ある年の東京の全天日射量の時刻間の関係の一例を示すグラフ。
【図６】ある年の１２月における東京の一定時間毎の全天日射量の一例を示すグラフ。
【図７】ある年の１２月における東京で「快晴」又は「晴れ」の日のみの一定時間毎の全天日射量の一例を示すグラフ。
【図８】同実施形態に係る情報処理装置の具体的な動作の二つの例を示す図。
【図９】同実施形態に係る情報処理装置が生成する発電量シナリオの一例を示すグラフ。
【図１０】同実施形態に係る情報処理装置が生成する需要量シナリオの一例を示すグラフ。
【図１１】同実施形態に係る情報処理装置の機能の一例を示す機能ブロック図であり、風力発電設備の発電量シナリオを生成する一例を示す図。
【図１２】本発明の第２の実施形態に係る情報処理装置の構成の一例を示す機能ブロック図。
【図１３】同実施形態に係る情報処理装置の具体的な動作の一例を示す図。
【図１４】同実施形態に係る情報処理装置が生成する必要発電量シナリオの一例を示すグラフ。
【図１５】本発明の第３の実施形態に係る情報処理装置の構成の一例を示す機能ブロック図。
【発明を実施するための形態】
【００２１】
以下、図面を参照して、本発明の実施の形態について説明する。
【００２２】
（第１の実施形態）
先ず、図１〜図１１を参照して、本発明の第１の実施形態について説明する。
【００２３】
図１は、本発明の各実施形態に共通する電力システムのハードウェア構成の一例を示す図である。
【００２４】
本発明の各実施形態に共通する電力システムに適用される情報処理装置１は、ＣＰＵ２、メモリ３、入力装置４、および表示装置５を有し、電力系統１００に接続される各種発電設備、需要設備、蓄電池等と各種データやコマンドの通信を行える（各種発電設備から過去の発電量の推移を示す発電量データを取得したり、需要設備から過去の需要量の推移を示す需要量データを取得したり、発電計画に基づいて蓄電池に対して充放電を指示するコマンドを送信したりすることができる）ようになっている。各種発電設備は、自然エネルギーを用いる太陽光発電１１１〜１１Ｍ、風力発電１２１等の発電設備と、その他の発電設備１３１〜１３Ｎ（例えば、火力発電、原子力発電等）とに分類される。また、電力系統１００には、自然エネルギーを用いる発電設備等から生成される電力の充放電を行う蓄電池１４１や、電力を直流から交流に変換する電力変換装置１５１が設けられる。電力系統１００に供給された電力は、需要設備１６１へ供給される。
【００２５】
なお、図１では、自然エネルギーを用いる発電設備が複数存在する場合を示しているが、この第１の実施形態では、最初に技術の基本部分を理解しやすいものとするため、自然エネルギーを用いる発電設備が、例えば「太陽光発電設備」１台のみである場合を想定して説明する。
【００２６】
図２は、本発明の第１の実施形態に係る情報処理装置の機能の一例を示す機能ブロック図である。
【００２７】
本発明の第１の実施形態に係る情報処理装置１は、情報保存部１０，類似日選択部１１，時系列データ生成部１２，数理モデル生成部１３，およびシナリオ生成部１４を有する。
【００２８】
情報保存部１０は、過去の気象の推移を示す気象情報、および自然エネルギーを用いる発電設備の過去の発電量の推移を示す発電量データを少なくとも記憶している記憶媒体である。なお、情報保存部１０に記憶されるデータは前述のものに限らず、例えば、電力系統１００に接続される需要設備の需要量を示す需要量データ、電力系統１００に含まれる自然エネルギーを用いる発電設備以外の設備（ここでは、蓄電池等を含むものとする。）の発電量を示す必要発電量データ等が挙げられる。また、情報保存部１０に記憶される過去の気象の推移を示す気象情報は、例えば、図３に示すような分類アルゴリズムを用いて天候に応じて分類される。図３に示す分類アルゴリズムは、気象庁の発表等を用いて、頻度の高い天候パターンを１日の時間経過に沿って分類１〜分類３のように分類した後に、分類４のように出現頻度の高い天候パターンを５〜２０種類程度リストアップし、リストアップされた天候パターンのそれぞれにパターン番号を付すことで、気象情報を分類している。なお、分類１および分類３では、「快晴」、「晴れ」、「曇り」、「雨」、「雪」等の天候の分類を行い、分類２では「時々」、「後」、「一時」、「一」等の条件の分類を行う。また、分類４では、分類１〜分類３の結果を用いて、「快晴」、「晴れ」、「晴れ時々曇り」等の出現頻度の高い天候パターンの分類およびリストアップを行う。なお、図３に示した分類アルゴリズムは一例でありこれに限られるものでない。その他の分類アルゴリズムの具体的な例としては、例えば、クラスター分析、サポートベクターマシン等が挙げられる。
【００２９】
類似日選択部１１は、情報保存部１０に記憶される過去の気象の推移を示す気象情報から、発電量の予測の対象となる日と気象の推移が類似する過去の日（例えば、過去の一定の期間において、一日の天候パターンが最も類似する日、あるいは一日の日射量の変化が最も類似する日）（以下、「類似日」と称す）を少なくとも二日分選択する機能を有する。この類似日選択部１１は、例えば図３に示したパターン番号が付された天候パターンから、発電量の予測の対象となる日と同様の天候パターンの番号を選択することで類似日を選択する。なお、類似日の分類では、気象庁の発表する雲量を示す雲量データおよび快晴指数（＝表面日射量／大気外日射量）等により、さらに詳細な分類を行うことも可能である。また、類似日を選択する日数が多いほど後述する発電量の予測の精度は高くなる。
【００３０】
時系列データ生成部１２は、類似日選択部１１により選択される日の発電量データを情報保存部１０から取得し、取得した発電量データを一定時間毎（例えば、「１時間毎」、「１分毎」等）の時系列データとして生成する機能を有する。なお、時系列データ生成部１２の具体的な処理については後述する。
【００３１】
数理モデル生成部１３は、時系列データ生成部１２により生成された時系列データを用いて、類似日の時刻間の発電量の相関を少なくとも示す数理モデルを生成する機能を有する。なお、生成する数理モデルの具体的な例としては、後述の分散共分散行列等が挙げられる。
【００３２】
シナリオ生成部１４は、数理モデル生成部１３で生成された数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布を決定するとともに決定した確率分布に従う一定時間毎の乱数の組を複数生成することによって自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを複数生成する機能を有する。なお、確率密度関数を用いた具体的な計算式については後述する。
【００３３】
また、図２に示される情報処理装置１は、前述の機能の他に、太陽の位置等を示す天文データを記憶する情報保存部１５、太陽光パネル特性から太陽光発電量の理論上限を求める理論上限生成部１６、および確率分布を決定する確率密度関数のパラメータを生成するパラメータ生成部１７を含んでいる。なお、情報保存部１５、理論上限生成部１６、およびパラメータ生成部１７は同実施形態に係る情報処理装置１に含まれるとは限らない。また、情報保存部１５に記憶されるデータは前述のものに限られない。
【００３４】
次に、図４〜図７を参照して、全天日射量の変動と、太陽光発電設備の発電量の変動との関係について説明する。
【００３５】
図４は、「日付」を横軸にとり、「全天日射量」を縦軸にとった、ある年の１２時における東京の全天日射量の分布の一例を示すグラフである。
【００３６】
同図に示されるように、全天日射量には上限が存在し、上限から下限（０）までの間でランダムに分布している。このようにランダムな分布になる要因としては、例えば、雲量の変化などが挙げられる。
【００３７】
図５は、「９時の全天日射量」を横軸にとり、「１２時の全天日射量」を縦軸にとった、ある年の東京の時刻間の全天日射量の関係の一例を示すグラフである。
【００３８】
同図に示されるように、９時の全天日射量が増加するにつれて１２時の全天日射量も増加している。このため、時刻間の全天日射量には相関があるといえる。
【００３９】
図６は、「時刻」を横軸にとり、「全天日射量」を縦軸にとった、ある年の１２月における東京の一定時間毎の全天日射量の一例を示すグラフである。
【００４０】
図７は、「時刻」を横軸にとり、「全天日射量」を縦軸にとった、ある年の１２月における東京で「快晴」又は「晴れ」の日のみの一定時間毎の全天日射量の一例を示すグラフである。
【００４１】
図６および図７に示されるように、全天日射量には上限が存在し、日中（６時から１８時）に大きく変動している。このように全天日射量が日中に大きく変動する要因としては、例えば、一定時間毎の雲量の変化等が挙げられる。
【００４２】
図４〜図７より、全天日射量には上限が存在し、且つ時刻間の全天日射量には相関があることが示される。一方で、雲量の変化は、予測前日の段階で大まかな予測をすることは可能であるが、ある時刻におけるある場所の雲量を正確に予測することは困難であり、誤差を考慮した予測が必要となる。このため、誤差を考慮することが可能である確率密度関数を用いる予測を行うことが適当である。また、図７では「快晴」又は「晴れ」の日のみ示したが、その他の天候の日（例えば、「雨」、「曇り」等）の全天日射量の正確な観測は非常に困難である。そのため、同実施形態に係る情報処理装置１では、全天日射量を決定論的に予測するのではなく、全天日射量の変動も含めて予測する。
【００４３】
次に、図８を参照して、同実施形態に係る情報処理装置１の具体的な動作の二つの例を示す。
【００４４】
図８（ａ）は、情報処理装置１が、自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを生成するまでの動作の一例を示す図である。
【００４５】
先ず、類似日選択部１１が、情報保存部１０から類似日を少なくとも二日分選択する（ステップＳ１０１）。
【００４６】
次いで、時系列データ生成部１２が、選択された類似日の発電量データを情報保存部１０から取得し、取得した類似日の発電量データから一定時間毎の時系列データを生成する（ステップＳ１０２）。
【００４７】
次いで、数理モデル生成部１３が、生成した時系列データから時刻間の発電量の相関を示す数理モデルを生成する（ステップＳ１０３）。
【００４８】
最後に、シナリオ生成部１４が、生成した数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布を決定するとともに確率分布に従う乱数の組を複数生成することによって自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを複数生成する（ステップＳ１０４）。
【００４９】
また、図８（ｂ）は、図８（ａ）と同様に、自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを、全天日射量データから生成する動作の一例を示す図である。
【００５０】
先ず、類似日選択部１１が、情報保存部１０から類似日を少なくとも二日分選択する（ステップＳ２０１）。
【００５１】
次いで、時系列データ生成部１２が、選択された類似日の全天日射量データを情報保存部１０から取得し、取得した類似日の全天日射量データから一定時間毎の時系列データを生成する（ステップＳ２０２）。
【００５２】
次いで、数理モデル生成部１３が、生成した時系列データから時刻間の全天日射量の相関を示す数理モデルを生成する（ステップＳ２０３）。
【００５３】
次いで、シナリオ生成部１４が、生成した数理モデルを所定の確率密度関数に適用して一定時間毎の全天日射量の確率分布を決定するとともに確率分布に従う乱数の組を複数生成することによって一定時間毎の全天日射量の予測値を示す全天日射量シナリオを複数生成する（ステップＳ２０４）。
【００５４】
最後に、シナリオ生成部１４が、全天日射量シナリオについて全天日射量から発電量への変換を行うことにより、自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを複数生成する（ステップＳ２０５）。
【００５５】
次に、時系列データ生成部１２、数理モデル生成部１３、およびシナリオ生成部１４が行う処理の具体的な例について説明する。
【００５６】
時系列データ生成部１２は、類似日選択部１１で選択された類似日の発電量データを情報保存部１０から取得し、取得した発電量データを、例えば、１時間毎に纏めることで、類似日の１時間毎の発電量データを集めた発電量データ｛ｘｔ｝を生成する。また、時系列データ生成部１２は、生成した発電量データ｛ｘｔ｝をｘ１からｘ２４まで並べることで、１時間毎の時系列データを、Ｎ日分（２４時間×Ｎ日）生成する。なお、前述の「ｔ」は時刻を示す変数であり、前述の「Ｎ」は選択した類似日の日数を示す２以上の整数である。また、同実施形態における時系列データは、取得した発電量データを１時間毎のデータとして生成されているが、これに限られるものでなく、１分毎又は更に短い時間毎のデータとして生成されるようにしてもよい。
【００５７】
数理モデル生成部１３が生成する数理モデルは、例えば、次の（１）式のように生成される。
【数１】
0005537304.tif 000002
【００５８】
上述の（１）式に示されたｖａｒ（ｘｔ）は、時系列データ生成部１２から得た発電量データ｛ｘｔ｝の分散を表し、ｃｏｖ（ｘｔ，ｘｕ）は、発電量データ｛ｘｔ｝と発電量データ｛ｘｕ｝との共分散を表している。なお、各時刻の発電量の変動幅は各時刻の分散によって示され、時刻間の発電量の相関は共分散によって示されている。また、前述の「ｕ」は、時刻ｔとは別の時刻を示す変数であり、「｛ｘｕ｝」は、時刻ｕの発電量データを示している。数理モデル生成部１３において、このような分散および共分散を求めるためには、類似日を少なくとも二日分選択する必要があり、選択する類似日の日数が多いほど発電量の予測精度は高くなる。
【００５９】
シナリオ生成部１４が決定する確率分布は、例えば、次の（２）式に示すような確率密度関数を用いて決定される。
【数２】
0005537304.tif 000003
【００６０】
なお、（２）式に示すμはｘの平均値ベクトルμ＝（μ１，μ２，μ３，…，μ２４）を示し、Ｔは転置行列を示している。また、一日分の発電量データ｛ｘｔ｝は２４次元のベクトルとみなせるため、（２）式に従う２４個の乱数の組を多数生成することにより、任意の数のシナリオを生成することが可能である。このように生成されたシナリオは、類似日の統計的な性質を再現している。なお、乱数の組の生成には、例えば特願２００９−５７６８２号に記載された方法を利用してもよい。
【００６１】
シナリオ生成部１４で用いる確率分布は、上述の（２）式に示した多次元正規分布に限られず、例えば、対数正規分布、指数分布、ワイブル分布、ベータ分布等を用いることも可能である。また、確率分布は、情報保存部１０より取得した発電量データから関数近似により決定することも可能である。
【００６２】
以下に、シナリオ生成部１４が、多次元正規分布以外の分布を用いる具体的な例として、ベータ分布を用いる場合について説明する。
【００６３】
ベータ分布の確率密度関数を、次の（３）式に示す。
【数３】
0005537304.tif 000004
【００６４】
ベータ分布は、太陽光発電のように、発電量の下限が０であり、発電量の上限が太陽の位置等を示す天文学的データ、太陽と太陽光発電パネルとの位置関係、および発電効率によって決定されるような、上限と下限とが決定しているような場合に用いられ、フィッティングの精度を高めることが可能である。
【００６５】
なお、（３）式に示すＢ（ｐ，ｑ）はベータ関数，ａは分布の下限，ｂは分布の上限，ｐおよびｑは分布を決定するパラメータを示している。
【００６６】
なお、（３）式によるベータ分布を用いる場合、最初にベータ分布に従う乱数を複数生成させてから、その後に前述の分散共分散行列を適用することにより、複数の乱数に少なくとも時刻間の発電量の相関を持たせるようにしてもよい。
【００６７】
次に、図９および図１０を参照して、情報処理装置１が生成するシナリオについて説明する。
【００６８】
図９（ａ）は、「時刻」を横軸にとり、「発電量」を縦軸にとった、情報保存部１０に記憶される過去の発電量の推移の一例を示すグラフであり、過去４日分の発電量データが示されている。図９（ｂ）は、図９（ａ）と同様に、「時刻」を横軸にとり、「発電量」を縦軸にとった、同実施形態における情報処理装置１が生成する発電量シナリオの一例を示すグラフである。なお、図９（ｂ）では、シナリオ生成部１４で生成する乱数の組を用いて３０通りのシナリオを生成した結果を示しているが、生成されるシナリオの数はこれに限られるものでなく、選択する類似日が二日分しかなくとも、例えば１００万通りのシナリオを生成することが可能である。
【００６９】
この図９（ａ）および図９（ｂ）を比較すると、図９（ｂ）の発電量シナリオは、図９（ａ）に示した時刻間の発電量の相関を再現しているといえる。
【００７０】
また、図９では太陽光発電設備の一定時間毎の発電量の予測値を示す発電量シナリオについて説明したが、同実施形態に係る情報処理装置１は、図１０に示すように、電力系統１００に接続される需要設備の一定時間毎の需要量の予測値を示す需要量シナリオを生成することも可能である。
【００７１】
図１０（ａ）は、「時刻」を横軸にとり、「需要量」を縦軸にとった、情報保存部１０に記憶される過去の需要量の推移の一例を示すグラフである。図１０（ｂ）は、図１０（ａ）と同様に、「時刻」を横軸にとり、「需要量」を縦軸にとった、同実施形態における情報処理装置１が生成する需要量シナリオの一例を示すグラフである。
【００７２】
この図１０（ａ）および図１０（ｂ）を比較すると、図１０（ｂ）の需要量シナリオは、図１０（ａ）に示した時刻間の需要量の相関を再現しているといえる。
【００７３】
なお、同実施形態に係る情報処理装置１は、前述で示した太陽光発電設備に限らず、他の自然エネルギーを用いる発電設備に用いることも可能であり、例えば、図１１に示すように風力発電設備に用いることも可能である。
【００７４】
図１１は、同実施形態に係る情報処理装置１の機能の一例を示す機能ブロック図であり、風力発電設備の一定時間毎の発電量の予測値を示す発電量シナリオを生成する一例を示す。
【００７５】
図１１に示される風力発電設備の場合は、情報保存部１５に記憶される地域毎の風況を示す風況データから風速の理論上限を求めた後に、風速の理論上限と風車特性とから、理論上限生成部１６で風力発電設備の発電量の理論上限を求め、パラメータ生成部１７において、確率分布を決定する確率密度関数のパラメータを生成する。なお、図１１に示した類似日選択部１１、時系列データ生成部１２、数理モデル生成部１３、およびシナリオ生成部１４は、情報保存部１０に記憶される、風力発電設備による過去の発電量データ等を用いて、前述と同様に機能する。
【００７６】
この第１の実施形態によれば、自然エネルギーを用いる発電設備を含む電力系統において、時刻間の発電量の相関を少なくとも考慮した発電量シナリオを容易に生成することが可能となる。
【００７７】
（第２の実施形態）
次に、図１２〜図１４を参照して、本発明の第２の実施形態について説明する。
【００７８】
なお、この第２の実施形態においては、図２に示した第１の実施形態の構成と共通する部分には同一の符号を付し、重複する説明を省略する。以下では、第１の実施形態と異なる部分を中心に説明する。
【００７９】
この第２の実施形態では、電力系統１００に含まれる自然エネルギーを用いる発電設備以外の設備に求められる一定時間毎の必要発電量の予測値を示す必要発電量シナリオの生成と、発電設備の起動、停止、および出力増減のタイミングを一定期間に渡って定めた発電計画の生成とについて説明する。
【００８０】
図１２は、本発明の第２の実施形態に係る情報処理装置１の構成の一例を示す機能ブロック図である。
【００８１】
同実施形態に係る情報処理装置１は、第１の実施形態に示した構成に、発電計画生成部１８と、発電計画選択部１９と、充放電量決定部２０と、容量決定部２１と、期待値算出部２２とを更に設けた構成とする。
【００８２】
但し、同実施形態に係る情報保存部１０は、過去の気象の推移を示す気象情報、自然エネルギーを用いる発電設備の過去の発電量の推移を示す発電量データのほか、需要設備の過去の需要量の推移を示す需要量データを記憶している。
【００８３】
また、類似日選択部１１は、情報保存部１０に記憶される過去の気象の推移を示す気象情報から、発電量および需要量の予測の対象となる日と気象の推移が類似する過去の日を少なくとも二日分選択する。
【００８４】
また、時系列データ生成部１２は、類似日選択部１１により選択される日の発電量データおよび需要量データを情報保存部１０から取得し、取得した発電量データおよび需要量データを一定時間毎の時系列データとして生成する。
【００８５】
また、数理モデル生成部１３は、時系列データ生成部１２により生成された時系列データを用いて、類似日選択部１１により選択される日の時刻間の発電量の相関を示すと共に時刻間の需要量の相関を示し、かつ、発電量と需要量との相関をも示す数理モデルを生成する。
【００８６】
また、シナリオ生成部１４は、数理モデル生成部１３により生成された数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布と、一定時間毎の需要量の確率分布とを決定するとともに、これらの確率分布にそれぞれ従う乱数の組を複数生成することによって、自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオおよび需要設備の一定時間毎の需要量の予測値を示す需要量シナリオを複数生成し、需要量シナリオから発電量シナリオを差し引くことによって自然エネルギーを用いる発電設備以外の設備（ここでは、蓄電池等を含むものとする。）に求められる一定時間毎の必要発電量の予測値を示す必要発電量シナリオを複数生成する。
【００８７】
発電計画生成部１８は、シナリオ生成部１４により生成された必要発電量シナリオを用いて、発電設備の起動，停止，および出力増減のタイミングを一定期間に渡って定めた発電計画を生成する機能を有する。この発電計画生成部１８は、例えば、シナリオ生成部１４により生成された必要発電量シナリオの平均値の変動に応じて蓄電池１４１を充放電させる制御を含む発電計画を生成する。
【００８８】
発電計画選択部１９は、生成した複数の発電計画から最適な発電計画を選択する機能を有する。この発電計画選択部１９の機能は、発電計画生成部１８に含まれていてもよい。
【００８９】
なお、充放電量決定部２０，容量決定部２１，および期待値算出部２２については後述する。
【００９０】
次に、図１３を参照して、同実施形態に係る情報処理装置１の具体的な動作の一例を示す。
【００９１】
図１３は、情報処理装置１が時刻間の発電量の相関と、時刻間の需要量の相関とを示す必要発電量シナリオを生成し、発電計画を生成するまでの動作の一例を示す図である。
【００９２】
先ず、類似日選択部１１が、情報保存部１０から類似日を少なくとも二日分選択する（ステップＳ３０１）。
【００９３】
次いで、時系列データ生成部１２が、選択された類似日の発電量データと需要量データとを情報保存部１０から取得し、取得した類似日のデータから一定時間毎の時系列データを生成する（ステップＳ３０２）。例えば、時系列データ生成部１２は、類似日選択部で選択された類似日の発電量データと需要量データとを情報保存部１０から取得し、取得したデータを１時間毎に纏めることで、１時間毎の発電量データ｛ｘｔ｝と１時間毎の需要量データ｛ｙｔ｝とを生成する。また、時系列データ生成部１２は、生成した発電量データ｛ｘｔ｝と需要量データ｛ｙｔ｝とを、それぞれｘ１からｘ２４まで、ｙ１からｙ２４まで並べることで、１時間毎の４８個の時系列データを生成する。
【００９４】
次いで、数理モデル生成部１３が、生成した時系列データから時刻間の発電量の相関と、時刻間の需要量の相関とを示す数理モデルを生成する（ステップＳ３０３）。例えば、数理モデル生成部１３は、生成した４８個の時系列データから、時刻間の発電量の相関と、時刻間の需要量の相関とを示す４８×４８の分散共分散行列を生成する。
【００９５】
次いで、シナリオ生成部１４が、生成した数理モデルを所定の確率密度関数に適用して一定時間毎の発電量の確率分布と、一定時間毎の需要量の確率分布とを決定するとともに、これらの確率分布にそれぞれ従う乱数の組を複数生成することによって、自然エネルギーを用いる発電設備の一定時間毎の発電量の予測値を示す発電量シナリオおよび需要設備の一定時間毎の需要量の予測値を示す需要量シナリオを複数生成する（ステップＳ３０４）。例えば、シナリオ生成部１４は、生成した数理モデルを４８次元正規分布に適用することで、一定時間毎の発電量の確率分布と需要量の確率分布とを決定し、決定したそれぞれの確率分布に従う４８個の乱数の組を生成することによって、発電量シナリオおよび需要量シナリオを生成する。
【００９６】
次いで、シナリオ生成部１４が、生成した需要量シナリオから発電量シナリオを差し引くことによって自然エネルギーを用いる発電設備以外の設備に求められる一定時間毎の必要発電量の予測値を示す必要発電量シナリオを複数生成する（ステップＳ３０５）。例えば、シナリオ生成部１４が生成した４８個の乱数の組の内、最初の２４個を発電量シナリオの生成に使用し、残りの２４個を需要量シナリオの生成に使用することで、需要量シナリオと発電量シナリオとの差を求めて、必要発電量シナリオを生成する。
【００９７】
最後に、発電計画生成部１８（および発電計画選択部１９）が、シナリオ生成部１４により生成された必要発電量シナリオを用いて発電計画を生成する（ステップＳ３０６）。例えば、シナリオ生成部１４により生成された、発電量シナリオと、需要量シナリオと、必要発電量シナリオとから発電設備の起動，停止，および出力増減のタイミングを一定期間に渡って定めた発電計画を生成する。
【００９８】
なお、図１２は、自然エネルギーを用いる発電設備が１台である場合について示したが、これに限られるものでない。自然エネルギーを用いる発電設備が複数台ある場合は、例えば、数理モデル生成部１３により生成される数理モデルを、自然エネルギーを用いる複数台の発電設備のぞれぞれの発電量の相関を示すものとして作成することにより、自然エネルギーを用いる複数台の発電設備に係る発電量シナリオ、需要量シナリオ、必要発電量シナリオ、および発電計画を容易に生成することが可能である。
【００９９】
図１４は、「時刻」を横軸にとり、「必要発電量」を縦軸にとった、同実施形態に係る情報処理装置１が生成する必要発電量シナリオの一例を示すグラフである。
【０１００】
この図１４の必要発電量シナリオは、図１０（ｂ）の需要量シナリオから、図９（ｂ）の発電量シナリオを差し引いた結果を示している。
【０１０１】
同実施形態に係る情報処理装置１は、発電量と需要量とに大きな相関がない場合においても、従来のように全ての変動の組み合わせを考慮する必要がないため、容易に必要発電量シナリオを生成することが可能である。
【０１０２】
次に、同実施形態に係る情報処理装置１の充放電量決定部２０と、容量決定部２１と、期待値算出部２２とについて説明する。
【０１０３】
同実施形態に係る情報処理装置１が生成したシナリオおよび発電計画は、図１２に示した、充放電量決定部２０，容量決定部２１，および期待値算出部２２に用いることが可能である。
【０１０４】
充放電量決定部２０は、シナリオ生成部１４により生成されたシナリオのそれぞれの値を平均し、電力系統１００に含まれる蓄電池において充放電すべき充放電量を求める機能を有する。このように、充放電量決定部２０が、蓄電池の充放電量を求めることで、蓄電池の充電および放電を制御することが可能になり、複雑に変動する自然エネルギーを用いる発電設備の発電量の変動を平滑化することが可能である。
【０１０５】
容量決定部２１は、シナリオ生成部１４により生成されたシナリオのそれぞれの値に応じて、電力系統１００に含まれる蓄電池の容量を決定する機能を有する。このように、容量決定部２１が、蓄電池の適度な容量を決定することにより、発電量の変動を効率よく平滑化することが可能である。
【０１０６】
期待値算出部２２は、シナリオ生成部１４により生成された発電計画に基づいて、自然エネルギーを用いる発電設備の発電量の期待値，発電コストの期待値，および収益の期待値の少なくとも一つを求める機能を有する。
【０１０７】
また、発電計画選択部１９，充放電量決定部２０，容量決定部２１，および期待値算出部２２は、例えば、実行が困難である発電計画である場合に蓄電池を利用する等の事前対策を行うような非常の事態に備えたシミュレーションを行い、非常事態の発生時には当該シミュレーションに沿った処理を実施するようにしてもよい。
【０１０８】
この第２の実施形態によれば、自然エネルギーを用いる発電設備を含む電力系統において、発電量シナリオのみならず、需要量シナリオおよび必要発電量シナリオ、ならびに発電計画を容易に生成することが可能となる。
【０１０９】
（第３の実施形態）
次に、図１５を参照して、本発明の第３の実施形態について説明する。
【０１１０】
この第３の実施形態では、自然エネルギーを用いる発電設備が複数ある場合のシナリオおよび発電計画の生成について説明する。
【０１１１】
図１５は、本発明の第３の実施形態に係る情報処理装置１の構成の一例を示す機能ブロック図である。
【０１１２】
なお、同図では、自然エネルギーを用いる発電設備であるＭ台の太陽光発電設備および１台の風力発電設備の過去の発電量データ、ならびに需要設備の過去の需要量データが、情報処理装置１の情報保存部１０に記憶されている場合について例示している。
【０１１３】
ここで、自然エネルギーを用いる発電設備が複数存在する場合の、情報処理装置１がシナリオを生成する動作の一例について説明する。
【０１１４】
先ず、類似日選択部１１が、１日分の時系列データを生成するために、情報保存部１０から類似日を少なくとも二日分選択する。次いで、時系列データ生成部１２が、情報保存部１０から選択された類似日の発電量データおよび需要量データを取得し、例えば、１日分の太陽光発電設備の時系列データＫ個、１日分の風力発電設備の時系列データＬ個、および１日分の需要設備の時系列データ２４個を生成する。次いで、数理モデル生成部１３が、生成したそれぞれの時系列データから、時刻間の発電量の相関と時刻間の需要量の相関とを示す、（Ｋ×Ｍ＋Ｌ＋２４）×（Ｋ×Ｍ＋Ｌ＋２４）の成分を有する分散共分散行列を生成する。最後に、シナリオ生成部１４が、生成した数理モデルを所定の確率密度関数に適用することで、一定時間毎の発電量と需要量との分布を決定し、決定した分布に従う（Ｋ×Ｍ＋Ｌ＋２４）次元の一定時間毎の乱数の組を生成し、太陽光発電設備と、風力発電設備と、需要設備とのシナリオを生成する。
【０１１５】
なお、時系列データ生成部１２により生成された時系列データの数は、前述の数に限られない。
【０１１６】
また、前述に示すように生成した発電量シナリオおよび需要量シナリオは、同一条件の下で生成したシナリオであるため、同時に生じる高い可能性を有する。このため、図１２と同様の方法で、必要発電量シナリオを生成することが可能であり、更に、生成した必要発電量シナリオから発電計画を生成することが可能である。なお、電力系統に含まれる蓄電池は、揚水発電設備と同様に扱うことが可能であり、電力需要が高まった時に反覆利用される。
【０１１７】
この第３の実施形態によれば、自然エネルギーを用いる発電設備が複数ある場合において、各発電設備における時刻間の発電量、需要量、および必要発電量の相関を考慮するだけでなく、発電設備間の相関をも考慮した、発電量シナリオ、需要量シナリオ、および必要発電量シナリオ、ならびに発電計画を生成することが可能となる。
【０１１８】
上述した各実施形態で述べた各種の機能や処理手順は、コンピュータプログラムとして、コンピュータにより読み取り可能な記憶媒体（例えば磁気ディスク，光ディスク，半導体メモリ）に記憶させておき、必要に応じてそれをプロセッサにより読み出して実行するようにしてもよい。また、このようなコンピュータプログラムは、通信媒体を介してあるコンピュータから他のコンピュータに伝送することにより配布することも可能である。
【０１１９】
本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【符号の説明】
【０１２０】
１…情報処理装置、２…ＣＰＵ、３…メモリ、４…入力装置、５…表示装置、１０…情報保存部、１１…類似日選択部、１２…時系列データ生成部、１３…数理モデル生成部、１４…シナリオ生成部、１５…情報保存部、１６…理論上限生成部、１７…パラメータ生成部、１８…発電計画生成部、１９…発電計画選択部、２０…充放電量決定部、２１…容量決定部、２２…期待値算出部、１００…電力系統、１１１，１１Ｍ…太陽光発電設備、１２１…風力発電設備、１３１，１３Ｎ…発電設備、１４１…蓄電池、１５１…電力変換装置、１６１…需要設備。
【図１】
0005537304.tif 000005
【図２】
0005537304.tif 000006
【図３】
0005537304.tif 000007
【図４】
0005537304.tif 000008
【図５】
0005537304.tif 000009
【図６】
0005537304.tif 000010
【図７】
0005537304.tif 000011
【図８】
0005537304.tif 000012
【図９】
0005537304.tif 000013
【図１０】
0005537304.tif 000014
【図１１】
0005537304.tif 000015
【図１２】
0005537304.tif 000016
【図１３】
0005537304.tif 000017
【図１４】
0005537304.tif 000018
【図１５】
0005537304.tif 000019