(57)【特許請求の範囲】
【請求項１】
水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水電力量の予測値を算出する予測システムであって、
前記水力発電所の貯水池へ流入する水の流入量の実績値と、前記水力発電所における発電電力量の実績値と、溢水電力量の実績値と、を、日時を示す情報と関連付けて記憶する第１記憶部と、
前記第１記憶部に記憶されている、前記流入量の実績値と、前記発電電力量の実績値と、前記溢水電力量の実績値と、に基づいて、流入量、発電電力量及び溢水電力量の関係を表す第１の回帰式を算出する第１分析部と、
予測対象日時における流入量の予測値と、発電電力量の予測値と、前記第１の回帰式と、に基づいて、前記予測対象日時における前記溢水電力量の予測値を算出する第１溢水量予測部と、
を有することを特徴とする予測システム。
【請求項２】
前記日時を示す情報は、季節を表す季節情報と関連づけて記憶されており、
前記第１分析部は、前記第１記憶部に記憶されている、季節ごとの前記流入量の実績値と、前記発電電力量の実績値と、前記溢水電力量の実績値と、に基づいて、各季節の前記第１の回帰式を算出し、
前記第１溢水量予測部は、前記予測対象日時の季節に係る前記第１の回帰式を用いて、前記予測対象日時における前記溢水電力量の予測値を算出する
ことを特徴とする請求項１に記載の予測システム。
【請求項３】
降水量の実績値を、日時を示す情報と関連付けて記憶する第２記憶部と、
前記降水量の実績値と、前記流入量の実績値に基づいて、降水量と、流入量の関係を表す第２の回帰式を算出し、予測対象日時における前記降水量の予測値と、前記第２の回帰式に基づいて、前記流入量の予測値を算出する流入量予測部を有する
ことを特徴とする請求項１又は２に記載の予測システム。
【請求項４】
気温の実績値を、日時を示す情報と関連付けて記憶する第３記憶部と、
前記気温の実績値と、前記発電電力量の実績値と、に基づいて、前記予測対象日時に係る気温と発電電力量の関係を表す第３の回帰式を算出し、前記気温の予測値と、前記第３の回帰式に基づいて、前記流入量の予測値を算出する発電電力量予測部を有する
ことを特徴とする請求項１乃至３いずれか一項に記載の予測システム。
【請求項５】
前記第１分析部は、前記第１記憶部に記憶されている、前記流入量の実績値と、前記発電電力量の実績値と、前記溢水電力量の実績値と、に基づいて、最小二乗法により、流入量、発電電力量及び溢水電力量の関係を表す前記第１の回帰式を算出する
ことを特徴とする請求項１乃至４いずれか一項に記載の予測システム。
【請求項６】
前記第１溢水量予測部は、溢水電力量の予測値と、発電電力量の予測値とを加算することにより、前記水力発電所における平均可能電力量を算出する
ことを特徴とする請求項１乃至５いずれか一項に記載の予測システム。
【請求項７】
水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水量の予測値を算出する予測システムであって、
前記水力発電所の貯水池へ流入する水の流入量の実績値と、前記水力発電所における発電電力量の実績値と、前記溢水量の実績値と、を、日時を示す情報と関連付けて記憶する第４記憶部と、
前記第４記憶部に記憶されている、前記流入量の実績値と、前記発電電力量の実績値と、前記溢水量の実績値と、に基づいて、流入量、発電電力量及び溢水量の関係を表す第４の回帰式を算出する第２分析部と、
予測対象日時における流入量の予測値と、発電電力量の予測値と、前記第４の回帰式と、に基づいて、前記予測対象日時における前記溢水量の予測値を算出する第２溢水量予測部と、
を有することを特徴とする予測システム。
【請求項８】
水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水電力量の予測値を算出する予測システムの動作方法であって、
前記予測システムは、前記水力発電所の貯水池へ流入する水の流入量の実績値と、前記水力発電所における発電電力量の実績値と、前記溢水電力量の実績値と、を日時を示す情報と関連付けて記憶する第１記憶部を備え、
前記第１記憶部に記憶されている、前記流入量の実績値と、前記発電電力量の実績値と、前記溢水電力量の実績値と、に基づいて、流入量、発電電力量及び溢水電力量の関係を表す第１の回帰式を算出する第１の工程と、
予測対象日時における流入量の予測値と、発電電力量の予測値と、前記第１の回帰式と、に基づいて、前記予測対象日時における前記溢水電力量の予測値を算出する第２の工程と、
を実行することを特徴とする予測システムの動作方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、水力発電所における溢水量、溢水電力量、平均可能電力量を予測する予測システム、予測システムの動作方法に関する。
【背景技術】
【０００２】
水力発電所は、貯水池に、電気を水のかたちで蓄えることが可能であり、需要の変動への対応が極めて容易であることから、急激な需要の変動、ピーク需要への対応供給力として活用されている。また、水力発電所は、石油・石炭等の燃料を用いないことから、クリーンエネルギーとして期待されている。
【０００３】
近年、このような水力発電所を有効利用するため、水力発電所への水の流入量を予測し、最適な発電計画を算出する試みがなされている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−２８２４３１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、流込式の水力発電所等は、河川の水を利用して発電することから、河川の水の増減に多大な影響を与えてしまう。そして、河川の下流へ放水する水の量は、水力発電所において発電に使用した水の量のみならず、設備の作業、事故、その他の原因により水力発電所から放流される溢水量も影響を与える場合がある。
【０００６】
一方、非常時等における電力供給能力を正確に把握する観点等から、溢水による溢水電力量、又は平均可能電力量（発電電力量に溢水電力量を加えた値）を予測することが求められている。
【０００７】
そこで、本発明は、水力発電所における溢水量、溢水電力量、及び平均可能電力量を予測する予測システム、予測システムの動作方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
前述した課題を解決する主たる本発明は、水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水電力量の予測値を算出する予測システムであって、前記水力発電所の貯水池へ流入する水の流入量の実績値と、前記水力発電所における発電電力量の実績値と、溢水電力量の実績値と、を、日時を示す情報と関連付けて記憶する第１記憶部と、前記第１記憶部に記憶されている、前記流入量の実績値と、前記発電電力量の実績値と、前記溢水電力量の実績値と、に基づいて、流入量、発電電力量及び溢水電力量の関係を表す第１の回帰式を算出する第１分析部と、予測対象日時における流入量の予測値と、発電電力量の予測値と、前記第１の回帰式と、に基づいて、前記予測対象日時における前記溢水電力量の予測値を算出する第１溢水量予測部と、を有することを特徴とする予測システムと、を有することを特徴とする予測システムである。
【０００９】
本発明の他の特徴については、添付図面及び本明細書の記載により明らかとなる。
【発明の効果】
【００１０】
本発明によれば、高い精度で、溢水量、溢水電力量、及び平均可能電力量を予測することができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の第１実施形態における水力発電所の構成を示す図である。
【図２】本発明の第１実施形態における予測システムを示す図である。
【図３Ａ】本発明の第１実施形態における予測装置のハードウェア構成を示す図である。
【図３Ｂ】本発明の第１実施形態における管理装置のハードウェア構成を示す図である。
【図４Ａ】本発明の第１実施形態における過去データのデータテーブルを示す図である。
【図４Ｂ】本発明の第１実施形態における過去データのデータテーブルを示す図である。
【図４Ｃ】本発明の第１実施形態における予測データのデータテーブルを示す図である。
【図４Ｄ】本発明の第１実施形態における予測データのデータテーブルを示す図である。
【図５】本発明の第１実施形態における予測装置の機能構成を示す図である。
【図６】本発明の第１実施形態における予測システムの動作のフローを示す図である。
【図７Ａ】水力発電所への水の流入量の時間ごとの変化を表すグラフを示す図である。
【図７Ｂ】水力発電所への水の流入量の時間ごとの変化を表すグラフを示す図である。
【図７Ｃ】水力発電所への水の流入量の時間ごとの変化を表すグラフを示す図である。
【図８】本発明の第２実施形態における確率分布を示すイメージ図である。
【発明を実施するための形態】
【００１２】
本明細書および添付図面の記載により、少なくとも以下の事項が明らかとなる。
【００１３】
＜第１実施形態＞
＝＝＝水力発電所の構成について＝＝＝
図１は、本実施形態における水力発電所の一例を説明するための図である。本実施形態における水力発電所は、貯水池１、取水口２、発電機３、放水口４、堰５を含んで構成される流込式の水力発電所である。
【００１４】
貯水池１は、河川や山林から水の流入を受ける（以下、貯水池１に流入した水の量を「流入量」という）。そして、貯水池１は、堰５により流入した水を貯め、所定量の水を取水口２より発電機３に放水する。発電機３は、取水口２より取水された水を受けて発電を行う。そして、発電機３は発電に使用した水（以下、発電機３が使用した水の量を「取水量」という）を放水口４に放水する。その後、放水口４に貯められた水は、一定量に達すると、河川に放水される構成となっている。
【００１５】
また、貯水池１に流入した水は、その他、河川維持のために堰５から放水されたり、排砂等のため、堰５から定期的に放水されることによっても、増減する（以下、それぞれの放水量を「責任放流量」、「ゲート放流量」という）。
【００１６】
このとき、取水口２より取水された発電機３が発電に使用する予定であった水の一部は、設備の作業、事故、その他の原因により、発電に利用されずに溢水する。溢水は、設備の作業により放水口４に放水する場合、水位の関係から地上に放水される場合、取水口２より取水されたにも関わらず、ゴミ等の影響により予期せず発電に使用されずに放水口４に放水する場合等、貯水池１に残留しなくなることを意味する（以下、溢水した水の量を「溢水量」という）。すなわち、取水口２より取水された水は、発電機３が発電に使用した取水量となる場合と、逸失により溢水量となる場合がある。
【００１７】
以上より、貯水池１に流入した流入量は、式（１）として表せる。
【００１８】
【数１】
0005908506.tif 000002
本実施形態では、溢水量は、発電に利用する取水口付近のゴミの量と関係を有するという理解に基づいて、ゴミの量に起因する説明変数を用いた重回帰モデルにより、溢水量を予測するものである。具体的には、貯水池１への流入量に応じて、ゴミが運ばれてくる量も増えるため溢水量は増加する点、取水口２に取水された際にゴミの影響により発電に使用されずに溢水となることから、発電電力量（取水量）に応じて、溢水量は増加する点に着目し、流入量、発電電力量（取水量）を、ゴミの量に起因する説明変数として用いる。すなわち、本実施形態では、当該回帰モデルから算出される溢水量予測式と、流入量の予測値と、発電電力量の予測値により、溢水量の予測値を算出する。
【００１９】
尚、発電機３が取水口２から取水する部分におけるゴミは、生活用品、木材、土砂等である。ゴミは、取水口２から取水された水とともに発電機３に取り込まれ、結果として発電機３の発電効率を悪化させ、取水した水量に対する発電電力量を得られなくする。また、その他、取り込まれたゴミは、発電機３のタービンと取水口２間に使用水量調整用の貯水槽にたまり、貯水槽から水を溢水させる原因となる。
【００２０】
＝＝＝予測システムの構成について＝＝＝
以下、図２〜図５を参照して、本実施形態における予測システムの構成について説明する。
【００２１】
図２に、本実施形態における溢水量の予測を実現する予測システムの一例を示す。本実施形態に係る予測システムは、予測装置１００と管理装置２００から構成される。両装置は、ＬＡＮ接続等による通信網３００を利用して、データの送受信を行う。
【００２２】
予測装置１００は、使用者が操作を行うコンピュータである。管理装置２００は、過去の発電電力量の実績値、過去の溢水電力量の実績値、過去の水力発電所への水の流入量の実績値等を記憶し、予測装置１００からのリクエストに応じてデータを送信するコンピュータである。本実施形態に係る予測装置１００は、管理装置２００から、予測対象となる水力発電所における過去の発電電力量の実績値、過去の溢水電力量の実績値、過去の流入量の実績値、予測対象期間における発電電力量の予測値、予測対象期間における流入量の予測値等を取得することで、予測対象期間における溢水電力量の予測値を算出する構成としている。
【００２３】
図３Ａに、本実施形態の予測装置１００のハードウェア構成を示す。
【００２４】
予測装置１００は、制御手段１００Ａ、記憶手段１００Ｂ、通信手段１００Ｃ、入力手段１００Ｄ、表示手段１００Ｅを有している。
【００２５】
制御手段１００Ａは、ＣＰＵ等であり、バス等を介して、記憶手段１００Ｂ、通信手段１００Ｃ、入力手段１００Ｄ、表示手段１００Ｅと接続されている。そして、制御手段１００Ａは、記憶手段１００Ｂに記憶されたコンピュータプログラムに基づいて、記憶手段１００Ｂ、通信手段１００Ｃ、入力手段１００Ｄ、表示手段１００Ｅとデータ通信を行うとともに、それらの動作を制御する。
【００２６】
記憶手段１００Ｂは、揮発性メモリー（ＲＡＭ）、不揮発性メモリー（フラッシュメモリー、磁気ディスク）等からなる。そして、記憶手段１００Ｂには、予測装置１００を制御するためのコンピュータプログラム等が記憶されている。尚、記憶手段１００Ｂは、後述する溢水量を予測する回帰モデル、各機能部で計算された中間データ、最終データ、取得した分析対象データ等を記憶する記憶部も有している（図示せず）。
【００２７】
通信手段１００Ｃは、通信コントローラ等であり、有線や無線によるＬＡＮ接続による通信網３００等を利用して、管理装置２００とデータの送受信を行う。
【００２８】
入力手段１００Ｄは、スイッチ、タッチパネル等であり、予測装置１００に対する使用者の操作指示を受付ける。
【００２９】
表示手段１００Ｅは、各種の情報を表示する液晶ディスプレイ等であり、後述する本システムにより算出した溢水電力量等を表示する。
【００３０】
図３Ｂに、本実施形態の管理装置２００のハードウェア構成を示す。
【００３１】
管理装置２００は、制御手段２００Ａ、記憶手段２００Ｂ、通信手段２００Ｃを有している。
【００３２】
制御手段２００Ａは、ＣＰＵ等であり、バス等を介して、記憶手段２００Ｂ、通信手段２００Ｃと接続されている。そして、制御手段２００Ａは、記憶手段２００Ｂに記憶されたコンピュータプログラムに基づいて、通信手段２００Ｃ、記憶手段２００Ｂとデータ通信を行うとともに、それらの動作を制御する。
【００３３】
記憶手段２００Ｂは、揮発性メモリー（ＲＡＭ）、不揮発性メモリー（フラッシュメモリー、磁気ディスク）等からなる。そして、記憶手段２００Ｂには、管理装置２００を制御するためのコンピュータプログラムと、後述する過去のデータ（以下、「過去データ」という）、後述する予測データ（以下、「予測データ」という）等が記憶されている。尚、過去データは、発電電力量、溢水電力量の過去の実績値に関するデータ（データテーブルＭ１）、及び水力発電所への流入量の過去の実績値に関するデータ（データテーブルＭ２）から構成される。また、予測データは、予測対象期間の予測される発電電力量（データテーブルＭ３）、予測対象期間の水力発電所への流入量に関するデータ（データテーブルＭ４）から構成される。
【００３４】
通信手段２００Ｃは、通信コントローラ等であり、有線や無線によるＬＡＮ接続による通信網３００等を利用して、予測装置１００とデータの送受信を行う。
【００３５】
図４Ａ、図４Ｂに、本実施形態の管理装置２００の記憶手段２００Ｂに記憶された過去データの一例を示す。
【００３６】
図４Ａは、水力発電所における過去の発電電力量の実績値、溢水電力量の実績値に関するデータテーブルＭ１である。このデータテーブルＭ１には、１時間ごとに、当該日時と対応させて発電電力量の実績値、及び溢水電力量の実績値が記憶されている。また、データテーブルＭ１には、当該日時と対応させて季節を表す季節情報が記憶されている。
【００３７】
ここで、発電電力量の実績値は、発電した発電電力量から算出されるものである。また、取水口２から発電機３に取り入れられた水の量を流量計により測定した値の実績値と、発電機３の設計上の定められた電水比（単位流量あたりの発電機出力）より算出されるものであってもよい。
【００３８】
また、溢水電力量の実績値は、溢水量の実績値と、発電機３の電水比（単位流量あたりの発電機出力）より算出される。尚、本実施形態では、溢水量は、上記式（１）から算出された値を用いている。具体的には、水力発電所における発電電力量と、発電機３の発電機３の設計上の定められた電水比（単位流量あたりの発電機出力）から取水量を算出し、流入量、責任放流量、ゲート放流量、貯水地の貯水量の変動量の測定結果より、式（１）に基づいて、溢水量を算出している。尚、溢水量の大部分が水力発電所に設置された貯水槽からの溢水である場合には、溢水量は、式（１）から算出する方法に代えて、貯水槽からの溢水量を流量計等で測定する方法であってもよい。
【００３９】
また、季節情報は、例えば、春、夏、秋、冬等の４区分で記憶されており、季節に依存したゴミの量の変化を反映させることを可能としている。
【００４０】
図４Ｂは、水力発電所における過去の貯水地１への水の流入量に関するデータテーブルＭ２である。このデータテーブルＭ２には、１時間ごとに、当該日時と対応させて水力発電所への水の流入量が記憶されている。
【００４１】
ここで、流入量は、例えば、河川を管理する自治体などが提供する河川の流量の測定結果である。また、その他、貯水池１における水深の変化から測定した値であってもよい。
【００４２】
図４Ｃ、図４Ｄに、本実施形態の管理装置２００の記憶手段２００Ｂに記憶された予測データの一例を示す。
【００４３】
図４Ｃは、予測対象となる水力発電所における発電電力量の予測値に関するデータテーブルＭ３である。このデータテーブルＭ３には、１時間ごとに、予測対象の日時と対応させて発電電力量の予測値が記憶されている。発電電力量の予測値は、例えば、電力需要量の予測値等から算出された、当該水力発電所における発電予定の電力量である。
【００４４】
図４Ｄは、予測対象となる水力発電所における、予測対象の日時の流入量の予測値に関するデータテーブルＭ４である。このデータテーブルＭ４には、１時間ごとに、予測対象の日時と対応させて流入量の予測値が記憶されている。流入量の予測値は、例えば、気象情報等から算出された、当該水力発電所に流入する水の量である。
【００４５】
尚、上記図４Ａ、図４Ｂ、図４Ｃ、図４Ｄでは、１時間ごとのデータとしたが、それらのデータを常時測定できる場合には、各時間との関係を対応付けたデータとして記憶されたものであってもよい。その場合、分析部１０２、溢水量予測部１０３は、流入量に起因するゴミの量をより正確に反映させるため、流入量の過去データ、及び予測データを、予測対象の日時よりも５分程度前にずらす等の処理を行ってもよい。
【００４６】
図５に、本実施形態の予測装置１００の機能構成の一例を示す。
【００４７】
予測装置１００は、記憶手段１００Ｂに記憶されたコンピュータプログラム、及び上記したハードウェア構成（１００Ａ〜１００Ｅ）により、以下に説明する取得部１０１、分析部１０２、溢水量予測部１０３、提示部１０４の機能を実現する。
【００４８】
取得部１０１は、通信手段１００Ｃにより管理装置２００と通信して、予測対象となる水力発電所に関する過去データ及び予測データを取得する。
【００４９】
分析部１０２は、発電電力量、溢水電力量及び流入量に関する過去データ（データテーブルＭ１、データテーブルＭ２）に基づいて、発電電力量、溢水電力量及び流入量の関係式を算出する。
【００５０】
溢水量予測部１０３は、分析部１０２により算出した発電電力量、溢水電力量及び流入量の関係式と、発電電力量に関する予測データ（データテーブルＭ３）と、流入量の予測データ（データテーブルＭ４）に基づいて、溢水電力量の予測値を算出する。
【００５１】
提示部１０４は、溢水量予測部１０４により算出された溢水電力量の予測値を、例えば、テキストデータ等として予測装置１００の表示手段１００Ｅに表示する。
【００５２】
尚、これらの機能の詳細は、後述する。
【００５３】
＝＝＝予測システムの動作について＝＝＝
次に、予測システムの動作の一例について説明する。
【００５４】
図６に、本実施形態のフローチャートを示す。
【００５５】
本実施形態では、予測装置１００の使用者が、水力発電所における予測対象日時の平均可能電力量を予測する場合について説明する。
【００５６】
（Ｓ１）は、使用者が、予測装置１００に対して予測対象の日時（例えば、２０１４／２／１の７時〜８時）を入力する工程である。
【００５７】
（Ｓ２）は、予測装置１００の取得部１０１が、管理装置２００に対して、過去データのうち、予測対象の水力発電所に関する、現在の季節と同一の季節情報の日の流入量の実績値、発電電力量の実績値、溢水電力量の実績値を要求する工程である。
【００５８】
尚、過去データとは、過去の所定日時における発電電力量の実績値、溢水電力量の実績値（データテーブルＭ１）、及び過去の所定日時に測定された流入量の測定結果（データテーブルＭ２）である。
【００５９】
ここで、取得すべきデータを現在の季節と同一の季節情報の日に限定しているのは、河川の上流から流れてくる水に含まれるゴミ、特に枯葉や木材の量は、季節によるところが大きいためである。したがって、季節ごとに溢水量予測式を算出することで、より精度の高い予測式とすることができる。尚、季節情報をデータテーブルＭ１に記憶していない場合には、月日に応じて季節を判断してもよい。
【００６０】
（Ｓ３）は、管理装置２００が当該要求を受けて、過去データの実測値に係るデータテーブルＭ１、Ｍ２から要求されたデータを取得し、予測装置１００に対して当該データを送信する工程である。
【００６１】
（Ｓ４）は、予測装置１００の分析部１０２が、取得した過去データに基づいて、溢水量予測式を算出する工程である。具体的には、予測装置１００は、取得した過去データに基づいて、発電電力量、溢水電力量、流入量を日時ごとに対応付けた、複数のデータからなる分析対象データテーブルを生成する。そして、予測装置１００は、当該分析対象データテーブルのデータに基づいて、回帰分析を実行し、発電電力量、溢水電力量、流入量の関係を示す溢水量予測式を算出する。尚、この算出処理の詳細は、後述する。
【００６２】
（Ｓ５）は、予測装置１００の取得部１０１が、管理装置２００に対して、予測データのうち、予測対象の日時（例えば、２０１４／２／１の７時〜８時）に関する当該水力発電所への流入量の予測値、及び発電電力量の予測値を要求する工程である。
【００６３】
（Ｓ６）は、管理装置２００が、当該要求を受けて、予測データのデータテーブルＭ３、Ｍ４から、予測対象の日時（例えば、２０１４／２／１の７時〜８時）に関する水力発電所への流入量、及び水力発電所における発電電力量の予測値を取得し、予測装置１００に対して、当該データを送信する工程である。
【００６４】
（Ｓ７）は、予測装置１００の溢水量予測部１０３が、（Ｓ６）で取得した、流入量の予測値、発電電力量の予測値、及び溢水量予測式に基づいて、予測対象の日時（例えば、２０１４／２／１の７時〜８時）の溢水電力量を算出する工程である。
【００６５】
尚、この工程で取得する流入量の予測値、発電電力量の予測値は、周知の流入量予測手法、発電電力量予測手法で算出された予測値であってよい。当該流入量の予測手法、発電電力量の予測手法の一例については、後述する。
【００６６】
（Ｓ８）は、予測装置１００の溢水量予測部１０３が、（Ｓ７）で算出した溢水電力量と、発電電力量の予測値とを加算することにより、当該水力発電所における平均可能電力量を算出する工程である。また、このとき、発電機３の電水比（単位流量あたりの発電機出力）により電力換算を行って、溢水電力量の予測値から、溢水量を算出してもよい。
【００６７】
（Ｓ９）は、予測装置１００の提示部１０４が、算出した１時間単位の溢水量、溢水電力量、平均可能電力量に対して所定の画像処理を施して、予測装置１００の使用者が認識できる態様で提示する工程である。提示部１０４は、例えば、１時間単位の溢水量、溢水電力量、平均可能電力量をテキストデータとして表示手段１００Ｅに表示する。
【００６８】
本実施形態の予測システムは、以上のようなプロセスによって、予測対象の日時の溢水量、溢水電力量、平均可能電力量の予測値を算出している。
【００６９】
尚、上記実施形態では、予測対象の日時として、例えば、２０１４／２／１の７時〜８時として１時間の溢水量を予測する場合について説明したが、時間帯ごとの溢水量を予測して、１日の推移を予測するものであってもよい。
【００７０】
また、上記実施形態では、予測装置１００の取得部１０１は、過去データと予測データとを別個のプロセスにより取得したが、過去データを取得する際、あわせて予測データを取得する態様としてもよい。
【００７１】
また、上記実施形態では、本予測システムは、予測対象の日時が設定されるに応じて、溢水量予測式を算出する工程（Ｓ４）を行うとしたが、（季節ごとの）当該溢水量予測式を予め生成しておき、予測対象の日時が設定されるに応じて、対応する季節の溢水量予測式を用いて、溢水量の予測値を算出してもよい。
【００７２】
＝＝＝溢水量予測式の算出について＝＝＝
次に、分析部１０２が、溢水量予測式を算出する工程（Ｓ４）について詳細に説明する。
【００７３】
本実施形態では、溢水量（溢水電力量）は、発電に利用する取水口付近のゴミの量と関係を有するという理解に基づいて、ゴミの量に起因する、水力発電所への水の流入量、及び取水量（発電電力量）を説明変数として用いる。
【００７４】
すなわち、本実施形態に係る予測装置１００は、溢水量予測式を算出するため、所定日時の溢水電力量Ｙｔを目的変数とし、少なくとも水力発電所への流入量Ｒｔ、水力発電所における発電電力量Ｆｔを説明変数とする回帰モデルである式（２）を用いる（式（２）では誤差項Ｅｔは省略して記載している）。
【００７５】
【数２】
0005908506.tif 000003
（ただし、α0は切片、α1、α２は回帰係数、Ｙは溢水電力量、Ｒは流入量、Ｆは発電電力量を表す。また、各変数の末尾のｔは該当する日時を表す。）
ここで、分析部１０２は、過去の発電電力量、溢水電力量、流入量に関する、日時ごとに対応付けられた、複数のデータからなる分析対象データテーブルに基づいて、回帰分析（例えば、最小二乗法）により、上記式（２）の回帰モデルの切片α0、回帰係数α1、α２を算出し、溢水量予測式を算出する。
【００７６】
尚、回帰分析で利用したサンプルと溢水量予測式から算出される決定係数Ｒ２値は、０．７程度と算出され、本回帰モデルにより、ばらつきが小さく、高い精度で予測することができる。尚、決定係数Ｒ２は、回帰モデルからの誤差Ｅｉの大きさを表す指標であり、決定係数Ｒ２が、０．５よりも大きな値となるとき、当該回帰モデルはばらつきが少なく精度が高いと判断できる。
【００７７】
上記実施形態では、回帰モデル（２）は、所定日時の溢水電力量Ｙを目的変数とし、少なくとも水力発電所への流入量Ｒ、水力発電所における発電電力量Ｆを説明変数とした。しかしながら、流入量Ｒに関するデータ、発電電力量Ｆに関するデータに基づいて、溢水電力量Ｙに関するデータを算出するものであれば、種々にバリエーション変更することができる。例えば、目的変数を溢水電力量Ｙに代えて溢水量としてもよいし、平均可能電力量としてもよい。この場合、例えば、回帰モデルは、式（３）として表すことができる。
【００７８】
【数３】
0005908506.tif 000004
（ただし、δ0は切片、δ1、δ２は回帰係数、Ｙ’は溢水量、Ｒは流入量、Ｆは発電電力量を表す。また、各変数の末尾のｔは該当する日時を表す。）
したがって、管理装置２００の記憶手段２００Ｂに、溢水量、発電電力量の実測値が記憶されている場合（図示せず）は、分析部１０２＊は、式（３）を用いて溢水量予測式を算出する。そして、溢水量予測部１０３＊は、当該溢水量予測式（３）を用いて、溢水量を算出する。また、回帰モデル式（２）において、流入量、発電電力量以外にも、説明変数として、上流で測定された土砂の量に関する情報等を追加してもよい。また、サンプルの分散を安定化させるため、分散安定化変換を行った回帰モデルを用いたり、１次回帰モデルに代えて高次回帰モデルとしてもよい。
【００７９】
また、上記実施形態では、溢水量予測式を算出する際、管理装置２００から、予測対象日の日時と同じ季節の過去データを取得する構成とした。しかし、これに代えて、予め式（２）の回帰モデルにダミー変数として季節に係る変数を追加し、分析部１０２が回帰分析を行う際に、他の季節情報が除外される方法であってもよい。この場合、同季節であるときには１、それ以外は０となるように設定すればよい。
【００８０】
＝＝＝流入量の予測について＝＝＝
水力発電所への水の流入量の予測値は、管理装置２００が備える流入量予測部２０１（図示せず）が算出する。以下、流入量予測部２０１が行う水力発電所への水の流入量の予測値の算出工程の一例を示す。
【００８１】
本実施形態の流入量予測部２０１は、水力発電所の上流における降水情報、降雪情報、気温、前の時間帯の流入量に基づいて、当該時間帯における水力発電所への水の流入量が決定するという理解に基づく流入量予測モデルを用いる。
【００８２】
図７Ａ、図７Ｂ、図７Ｃに、水力発電所への水の流入量の時間ごとの変化を表すグラフを示す。図７Ａ、図７Ｂ、図７Ｃは、それぞれ降水や降雪がない場合、降水があった場合、降雪があった場合の流入量の時間ごとの変化を表している。
【００８３】
図７Ａに示すように、降水や降雪がない場合にも、例えば、山林などからの滲出によって、所定の流入量が存在する。そして、その流入量は、降水や降雪がない場合、時間に応じて次第に減少し、一定量に収束する（以下、「均衡流入量」という）。一方、降水があった場合、図７Ｂに示すように、流入量は一時的に増加し、降水が止むと図７Ａと同様の均衡流入量に向けて減少する。そして、降水が降雪であった場合、図７Ｃに示すように、融雪するに応じて河川へ流入するため、流入量の増加速度及び減少速度は、通常の降水の場合と比して遅いものとなる。すなわち、降雪に基づく流入量は、降水とは別に融雪に基づく算出式が必要となる。
【００８４】
以上の理解に基づいて、本実施形態の流入量予測部２０１は、１時間毎に測定された流入量Ｒｔ、降水量Ｐ１ｔ、降雪量Ｐ２ｔ、融雪量Ｍｔ、積雪量Ｄｔ、気温Ｔｔに基づいて、以下の回帰モデル（４）、（５）、（６）、（７）により流入量予測式を算出する。
【００８５】
まず水力発電所への流入量Ｒｔは、均衡流入量μへ逓減する前の時間（本実施形態では１時間前）の流入量Ｒｔ-1からの変化量（μ-Ｒｔ-1）、降水による流入量の変化Ｐ１ｔ、融雪による流入量の変化Ｍｔを説明変数とする流入量モデル（４）を用いて表す。
【００８６】
尚、以下の回帰モデルでは、誤差項Ｅｔは省略して記載している。
【００８７】
【数４】
0005908506.tif 000005
（β1、β２、β3は回帰係数、Ｒは流入量、Ｐ１は降水量の一部、Ｍは融雪量、μは均衡流入量を表す。また、各変数の末尾のｔは該当する日時を表す。）
ここで、流入量モデル（４）の各説明変数は、以下のような構成となる。
【００８８】
まず当該時間帯の降水量Ｐｔは、降水量モデル（５）を用いて表せる。
【００８９】
【数５】
0005908506.tif 000006
（Ｐは降水量、Ｐ１は流入する降水量、Ｐ２は降雪となる降雪量を表す。）
降水量モデル（５）は、降水量Ｐｔが、当該日時の降水量を表し、そのときの気温Ｔｔが０℃以上である場合はすぐに流入するＰ１ｔとなるのに対し、０℃以下である場合は降雪となりＰｔは降雪量Ｐ２ｔとなることを表している。
【００９０】
そして、降雪となった場合、降雪量Ｐ２ｔは積雪した後、融雪し、河川へ流入するから、まず積雪量Ｄとして積雪量モデル（６）を用いて表せる。
【００９１】
【数６】
0005908506.tif 000007
（Ｄは積雪量、Ｐ２は降雪量、Ｍは融雪量を表す。）
積雪量モデル（６）は、積雪量Ｄｔが、前の時間帯の積雪量Ｄｔ-1に当該時間帯の積雪量を加算して、そのときの融雪量を減算して表せることを意味している。
【００９２】
そして、気温が上昇した場合、及び降水があった場合、積雪量Ｄｔ-1の一部が融雪して、融雪量Ｍｔとなるから、融雪量Ｍｔは、融雪量モデル（７）を用いて表せる。
【００９３】
【数７】
0005908506.tif 000008
（γ1、γ２は回帰係数、Ｔ０は融雪基準温度、Ｔは気温、Ｄは積雪量、Ｐ１は降水量、Ｍは融雪量を表す。）
融雪量モデル（７）は、積雪量Ｄｔ-1が残存してない場合には融雪量Ｍｔは０となり、気温Ｔｔが融雪基準温度Ｔ０より低い場合には、気温による融雪はないことを意味している。
【００９４】
以上より、水力発電所への流入量Ｒｔは、均衡流入量μへ逓減する前の時間（本実施形態では１時間前）の流入量Ｒｔ-1からの変化量（μ-Ｒｔ-1）、降水による流入量の変化Ｐ１ｔ、融雪による流入量の変化Ｍｔとの関係式（４）で表せる。
【００９５】
尚、管理装置２００は、均衡流入量μ、流入量Ｒｔ、降水量Ｐｔ、融雪量Ｍｔ、積雪量Ｄｔ、気温Ｔｔ等、過去の気象情報、流入情報を記憶手段２００Ｂに記憶している。また、管理装置２００は、気象庁等から発表された予測対象日における、水力発電所の上流の降水量の予測値と予測気温（例えば、ＧＰＶデータ）を記憶手段２００Ｂに記憶している（図示せず）。尚、河川への流入が遅い地域であれば、現在の雨量観測所での測定結果（例えば、アメダス雨量データ）を降水量の予測値と予測気温として、流入量の予測値を算出してもよい。
【００９６】
そして、管理装置２００の流入量予測部２０１は、過去の気象情報、流入情報に基づいて、回帰分析、例えば、最小二乗法により、それぞれの回帰係数β1、β２、β3、γ1、γ２を算出する。
【００９７】
これによって、流入量予測部２０１は、流入量モデル（４）、降水量モデル（５）、積雪量モデル（６）、融雪量モデル（７）からなる、流入量予測式を算出することができる。そして、流入量予測部２０１は、予測対象日の降水量の予測値と、予測気温より、例えば、１時間単位で流入量Ｒの予測値を算出する。
【００９８】
尚、上記実施形態では、流入量モデル（４）、降水量モデル（５）、積雪量モデル（６）、融雪量モデル（７）からなる、流入量予測式を算出したが、本実施形態で用いる回帰モデルは、上記に限る必要はない。例えば、温かい地域においては、降雪は考慮する必要がないため、流入量は降水量のみを説明変数とする式（４）のみの回帰モデルとすることもできる。この場合、回帰モデルは、式（８）として表すことができる。
【００９９】
【数８】
0005908506.tif 000009
他方、より精度を向上させるためには、当該河川に流入する地域をメッシュ状に複数に区分して、それぞれの箇所からの流入速度等を勘案して、流入量を算出してもよい。また、必要に応じて、説明変数として天気に関する情報等を追加してもよい。また、サンプルの分散を安定化させるため、分散安定化変換を行った回帰モデルを用いてもよい。
【０１００】
＝＝＝発電電力量の予測について＝＝＝
水力発電所における発電電力量の予測値は、管理装置２００が備える発電電力量予測部２０２（図示せず）が算出する。以下、発電電力量予測部２０２が行う水発電電力量予測部２０２の予測値の算出工程の一例を示す。
【０１０１】
本実施形態の発電電力量予測部２０２は、１日の中での時間帯、当該時間帯の気温に基づいて、水力発電所における発電電力量が決定するという理解に基づいて、発電電力量予測モデル（９）を用いる。
【０１０２】
【数９】
0005908506.tif 000010
（δ0ｊは切片、δ1ｊは回帰係数、Ｆは発電電力量、Ｔは気温、｜Ｔｔ―１８｜は冷暖房が使用される基準温度１８℃との差の絶対値を表す。また、変数の末尾のｔは該当する日時を表し、切片δ0ｊ、回帰係数δ1ｊの末尾のｊは１日の中での時間帯を表す。）
発電電力量予測モデル（９）は、１日の２４時間を１時間単位に区切るとともに、発電電力量を冷暖房需要量による説明変数｜Ｔｔ―１８｜により表すことを意味している。発電電力量は電力需要量に合致するように算出されるが、電力需要量は、１日の中での時間帯に依存する。例えば、電力需要量の大部分を占める工場が稼働する時間帯や、施設における冷暖房は、夜間には小さくなるが、昼間は大きくなる。一方、冷暖房による電力の需要量は、気温に依存する。
【０１０３】
したがって、本実施形態の発電電力量予測部２０２は、１日の２４時間を１時間単位に区切るとともに、冷暖房による電力需要量の説明変数｜Ｔｔ―１８｜を用いて、各時間帯に対してδ0ｊ、δ1ｊをそれぞれ算出する。尚、管理装置２００は、過去の発電電力量Ｆｔ、気温Ｔｔ等を記憶手段２００Ｂに記憶している（図示せず）。そして、管理装置２００の発電電力量予測部２０２は、これらに基づいて、回帰分析、例えば、最小二乗法により、各時刻に対して、切片δ0ｊ、回帰係数δ1ｊを算出する。
【０１０４】
これによって、発電電力量予測部２０２は、発電電力量予測式を算出することができる。尚、本実施形態における時間帯の意義は、１日の時間の推移の中での所定の期間を意味し、時間帯による電力需要量の変化を表せる範囲を意味する。したがって、変動が小さければ１時間単位に代えて２時間単位であってもよい。また、必ずしも時分で表された期間に限る必要はない。
【０１０５】
尚、上記実施形態では、発電電力量予測モデル（９）を用いたが、本実施形態で用いる回帰モデルは、上記に限る必要はない。例えば、工場立地帯では、冷暖房による電力需要量の変動幅は、施設稼働のために使用する電力量と比較して小さいから、説明変数｜Ｔｔ―１８｜は省略してもよい。
【０１０６】
他方、より精度を向上させるためには、電力需要量に時刻による変動が異ならせるため、地域をメッシュ状に複数に区分して、それぞれの箇所からの電力需要量を勘案して、総合した電力需要量を算出してもよい。また、電力供給者が、水力発電所以外にも原子力発電所や火力発電所等、複数の発電設備を有している場合には、それらの稼働に伴うコストを勘案して、必要とされる電力需要量のうちの一部を水力発電所における発電電力量として算出してもよい。また、必要に応じて、説明変数として天気に関する情報等を追加してもよい。また、サンプルの分散を安定化させるため、分散安定化変換を行った回帰モデルを用いてもよい。
【０１０７】
以上、本実施形態によれば、発電電力量の予測値と、流入量の予測値とに基づいて、溢水量、溢水電力量、又は平均可能電力量を算出することができる。これによって、例えば、水力発電所の下流の河川の水の増減を予測することができ、非常時等における電力供給能力を正確に予測することができる。
【０１０８】
＜第２実施形態＞
本実施形態の予測装置１００は、水力発電所における溢水電力量の予測値を、確率分布として算出する点で、第１実施形態と異なっている。以下、本実施形態の態様について説明する。尚、本実施形態では、分析部１０２’、溢水量予測部１０３’の構成のみ第１実施形態と異なっているため、その他の構成については説明を省略する。
【０１０９】
溢水電力量の予測値は、回帰式より算出されるものであるから、一定のばらつきを含む値となっている。そこで、本実施形態では、溢水電力量の予測値を確率分布として算出し、予測システムの使用者が、溢水電力量の実現し得る値が、予測値からどの程度外れるか等、リスクを把握できる態様としている。
【０１１０】
尚、本実施形態における確率分布（以下、「確率分布データ」という）とは、実現し得る値の予測値からのばらつきを示すものであり、例えば、予測対象日（例えば、２０１３／２／１の７時）に９５％の確率で実現し得る溢水電力量の幅（信頼区間）を意味する。また、確率分布データは、予測値に関する確率密度関数を表すデータや、実現し得る確率と予測値の幅の対応関係を示すデータ、予測値に関する分散係数等であってもよい。
【０１１１】
本実施形態では、図６の（Ｓ４）の工程において分析部１０２’が、式（２）に関して、ばらつきを示す係数を算出し（Ｓ４’）、（Ｓ７）の工程において溢水量予測部１０３’が、溢水予測式と当該ばらつきを示す係数に基づいて、確率分布データを算出する（Ｓ７’）点で、第１実施形態と異なっている。
【０１１２】
（Ｓ４’）では、式（１０）に基づいて、ばらつきを示す係数を算出する。
【０１１３】
溢水量のばらつきを示す係数は、例えば、式（２）の回帰モデルにおける誤差項Ｅｉ（式（２）では省略して記載している）の標準偏差σである。誤差項の標準偏差σは、サンプルが、平均的に回帰式からどれだけ散らばっているかを表す指標の１つである。
【０１１４】
誤差項の標準偏差σは、式（１０）の平方根で表すことができる。
【０１１５】
【数１０】
0005908506.tif 000011
（Ｙ＾ｔは、回帰分析においてサンプルとした各日時の流入量Ｒｔ、発電電力量Ｆｔについて、回帰式（２）より算出したＹｔ値を表す。ｎは、サンプル数を表す。）
本実施形態では、各実測値の分布が、以下の正規分布に従うとして、溢水電力量Ｙの確率分布は、説明変数の値によらず、一律に誤差項の標準偏差σ２を母分散とする正規分布とみなす。
【０１１６】
【数１１】
0005908506.tif 000012
この場合、上記溢水予測式、及び標準偏差σから、溢水電力量Ｙの確率分布を表すことができる。例えば、溢水量Ｙの確率分布を、確率密度関数ｆ（Ｙ’）で表す場合、式（１２）のように表される。
【０１１７】
【数１２】
0005908506.tif 000013
（Ｙ’は算出された溢水予測式における溢水電力量の予測値Ｙからずれた値を表す）
また、溢水電力量Ｙｔの確率分布データは、信頼区間として算出することもできる。
【０１１８】
尚、信頼区間とは、一定確率の範囲内で現実に起こり得る数値範囲を意味する。図８は、溢水量Ｙｔの確率密度関数ｆ（Ｙ’）に基づいて算出された信頼区間のイメージ図である。本実施形態の場合、６８．２％の信頼区間は、溢水量の予測値ＹからσだけずれたＹ±σの間の値となる。このとき、確率に応じた信頼区間は、同様に確率密度関数ｆ（Ｙ’）にしたがって、標準偏差σに対して信頼度に応じた係数を乗ずることによって、算出することができる。例えば、信頼度７０％の場合、標準偏差σに対して係数１．２８を乗ずればよいし、信頼度９０％の場合、標準偏差σに対して係数１．５３を乗ずればよい。
【０１１９】
尚、上記実施形態では、溢水電力量の確率分布は、説明変数（流入量Ｒ、発電電力量Ｆ）の値によらず、一律に誤差項の標準偏差σ２を母分散とする正規分布とみなした。しかし、溢水電力量の信頼区間は、説明変数の値（流入量Ｒ、発電電力量Ｆの値）によって分布関数を異ならせる、t分布に基づく信頼区間として算出してもよい。その場合、溢水電力量のばらつきに関する係数は、自由度ｎ−３で割った不偏推定量として、誤差項の標準偏差σ’は式（１３）で表せる。
【０１２０】
【数１３】
0005908506.tif 000014
このとき、回帰式（２）の分布が、以下の正規分布に従うと推定して、ｔ値より、溢水電力量の予測値Ｙの信頼区間を算出することができる。
【０１２１】
【数１４】
0005908506.tif 000015
（各観測点の平均値が回帰式（２）上の値となり、分散係数が誤差項の標準偏差σ’に基づいて決定されることを表す。尚、Ｄ２は流入量Ｒ、発電電力量Ｆにより定まるマハラノビスの距離の２乗を表す）
この場合の溢水電力量の予測値Ｙの信頼区間は、流入量Ｒ、発電電力量Ｆ、誤差項の標準偏差σ’に基づいて、式（１５）によって算出することができる。
【０１２２】
【数１５】
0005908506.tif 000016
以上のように本実施形態によれば、溢水電力量の予測値を確率分布として算出することができる。これによって、予測システムの使用者が、溢水電力量の実現し得る値が、予測値からどの程度外れるか等、リスクを把握できる態様としている。
【０１２３】
また、上記実施形態では、溢水電力量の確率分布を示すデータとして、他の分布関数を利用したものであってもよく、χ２分布やＦ分布に基づくデータであってもよい。
【０１２４】
また、上記実施形態では、流入量、発電電力量の予測値のばらつきは、考慮しない態様とした。しかしながら、式（４）、式（９）より算出されるそれらの誤差項に関しても、誤差項の標準偏差を算出し、溢水電力量の確率分布を算出する際に考慮に入れる態様としてもよい。その場合、周知の誤差伝搬の法則により、標準偏差を統合すればよい。
【０１２５】
尚、上記の計算式は、統計学上周知の計算方法のため、詳細な説明は省略する。
【０１２６】
＜その他実施形態＞
上記実施形態では、回帰分析の手法として、最小二乗法を示した。しかし、最小二乗法に代えて、モーメント法や最丈法を利用してもよい。
【０１２７】
また、上記実施形態では、予測装置１００が、分析部１０２、溢水量予測部１０３を機能部として有する構成とした。しかし、これらの機能部又はその一部は、他の装置に分散されていてもよい。同様に、各記憶手段に記憶されたデータの記憶領域は、任意の場所でよい。例えば、予測装置１００に集約されていてもよいし、複数のコンピュータから構成されるクラウドシステム上に分散して記憶される構成であってもよい。
【０１２８】
＝＝＝結言＝＝＝
以上より、上記各実施形態は、次のように記載できる。
【０１２９】
上記各実施形態は、水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水電力量の予測値（Ｙ）を算出する予測システムであって、水力発電所の貯水池へ流入する水の流入量の実績値（Ｒｔ）と、水力発電所における発電電力量の実績値（Ｆｔ）と、溢水電力量の実績値（Ｙｔ）と、を、日時を示す情報と関連付けて記憶する第１記憶部（データテーブルＭ１、Ｍ２に対応する）と、第１記憶部に記憶されている、流入量の実績値と、発電電力量の実績値と、溢水電力量の実績値と、に基づいて、流入量、発電電力量及び溢水電力量の関係を表す第１の回帰式（式（２）に対応する）を算出する第１分析部（１０２、１０２’に対応）と、予測対象日時における流入量の予測値（Ｒ）（データテーブルＭ４に対応する）と、発電電力量の予測値（Ｆ）（データテーブルＭ３に対応する）と、回帰式と、に基づいて、予測対象日時における溢水電力量の予測値を算出する第１溢水量予測部（１０３、１０３’に対応する）と、を有することを特徴とする予測システムを開示するものである。
【０１３０】
これによって、高い精度で、溢水量、溢水電力量、及び平均可能電力量を予測することができる。
【０１３１】
ここで、日時を示す情報は、季節を表す季節情報と関連づけて記憶されており、第１分析部は、第１記憶部に記憶されている、季節ごとの流入量の実績値と、発電電力量の実績値と、溢水電力量の実績値と、季節情報に基づいて、予測対象日時の季節に関して係る、季節ごとに、各季節の流入量と、発電電力量と、溢水電力量との関係を表す第１の回帰式を算出し、第１溢水量予測部は、予測対象日時の季節に係る第１の回帰式を用いて、予測対象日時における溢水電力量の予測値を算出するものであってもよい。
【０１３２】
これによって、季節によるゴミの量の違いを反映させ、より高い精度で、溢水量、溢水電力量、及び平均可能電力量を予測することができる。
【０１３３】
ここで、予測システムは、降水量の実績値を、日時を示す情報と関連付けて記憶する第２記憶部（記憶手段２００Ｂに記憶されている（図示せず））と、降水量の実績値と、流入量の実績値に基づいて、降水量と、流入量の関係を表す第２の回帰式を算出し、予測対象日時における降水量の予測値と、第２の回帰式（式（８）に対応する）に基づいて、流入量の予測値を算出する流入量予測部（２０１）を有していてもよい。
【０１３４】
ここで、予測システムは、気温の実績値を、日時を示す情報と関連付けて記憶する第３記憶部（記憶手段２００Ｂに記憶されている（図示せず））と、気温の実績値と、発電電力量の実績値と、に基づいて、予測対象の時間帯に係る気温と発電電力量の関係を表す第３の回帰式（式（９）に対応する）を算出し、予測対象日時における時間帯と、気温の予測値と、第３の回帰式に基づいて、流入量の予測値を算出する発電電力量予測部（２０２）を有していてもよい。
【０１３５】
ここで、第１分析部は、第１記憶部に記憶されている、流入量の実績値と、発電電力量の実績値と、溢水電力量の実績値と、に基づいて、最小二乗法により、流入量、発電電力量及び溢水電力量の関係を表す回帰式を算出してもよい。
【０１３６】
ここで、第１溢水量予測部は、溢水電力量の予測値と、発電電力量の予測値とを加算することにより、水力発電所における平均可能電力量を算出してもよい。
【０１３７】
また、上記各実施形態は、水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水量の予測値を算出する予測システムであって、水力発電所の貯水池へ流入する水の流入量の実績値と、水力発電所における発電電力量の実績値と、溢水量の実績値と、を、日時を示す情報と関連付けて記憶する第４記憶部（記憶手段２００Ｂに記憶されている（図示せず））と、第４記憶部に記憶されている、流入量の実績値と、発電電力量の実績値と、溢水量の実績値と、に基づいて、流入量、発電電力量及び溢水量の関係を表す第４の回帰式（式（３）に対応する）を算出する第２分析部（１０２＊に対応する）と、予測対象日時における流入量の予測値と、発電電力量の予測値と、第４の回帰式と、に基づいて、予測対象日時における溢水量の予測値を算出する第２溢水量予測部（１０３＊に対応する）と、を有することを特徴とする予測システムを開示するものである。
【０１３８】
これによって、高い精度で、溢水量、溢水電力量、及び平均可能電力量を予測することができる。
【０１３９】
また、上記各実施形態は、水力発電所において発電用に取水される水のうち、溢水によって逸失する溢水電力量の予測値を算出する予測システムの動作方法であって、前記予測システムは、前記水力発電所の貯水池へ流入する水の流入量の実績値と、前記水力発電所における発電電力量の実績値と、前記溢水電力量の実績値と、を日時を示す情報と関連付けて記憶する第１記憶部を備え、前記第１記憶部に記憶されている、前記流入量の実績値と、前記発電電力量の実績値と、前記溢水電力量の実績値と、に基づいて、流入量、発電電力量及び溢水電力量の関係を表す第１の回帰式を算出する第１の工程と、予測対象日時における流入量の予測値と、発電電力量の予測値と、前記第１の回帰式と、に基づいて、前記予測対象日時における前記溢水電力量の予測値を算出する第２の工程と、を実行することを特徴とする予測システムの動作方法を開示するものである。
【０１４０】
これによって、高い精度で、溢水量、溢水電力量、及び平均可能電力量を予測することができる。
【０１４１】
以上、本発明の具体例を詳細に説明したが、これらは例示にすぎず、請求の範囲を限定するものではない。請求の範囲に記載の技術には、以上に例示した具体例を様々に変形、変更したものが含まれる。
【符号の説明】
【０１４２】
１００  予測装置
１０１  取得部
１０２  分析部
１０３  溢水量予測部
１０４  提示部
２００  管理装置
２０１  流入量予測部
２０２  発電電力量予測部
３００  通信網
Ｙ  溢水電力量
Ｒ  流入量
Ｆ  発電電力量
Ｐ  降水量
Ｔ  気温
【図１】
0005908506.tif 000017
【図２】
0005908506.tif 000018
【図３Ａ】
0005908506.tif 000019
【図３Ｂ】
0005908506.tif 000020
【図４Ａ】
0005908506.tif 000021
【図４Ｂ】
0005908506.tif 000022
【図４Ｃ】
0005908506.tif 000023
【図４Ｄ】
0005908506.tif 000024
【図５】
0005908506.tif 000025
【図６】
0005908506.tif 000026
【図７Ａ】
0005908506.tif 000027
【図７Ｂ】
0005908506.tif 000028
【図７Ｃ】
0005908506.tif 000029
【図８】
0005908506.tif 000030