(57)【要約】【課題】取得可能な発電システムの情報が限られている場合であっても、より詳細な発電量を推定可能にする。【解決手段】推定システムは、対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とに基づいて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する変換パラメータ推定手段１０１を備える。【選択図】図１７
2018007461.tif 000002
【特許請求の範囲】
【請求項１】
対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、前記対象発電システムの設置位置または前記対象地域における、前記実績発電量を算出した期間中の前記第１の時間単位の日射量の系列である日射量列とに基づいて、前記第２の時間単位中の前記第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する変換パラメータ推定手段を備えた
ことを特徴とする推定システム。
【請求項２】
予め、変換情報として、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報を記憶する変換情報記憶手段を備え、
変換パラメータ推定手段は、前記変換情報によって特定される変換係数列の候補または制約条件を用いて、変換パラメータを推定する
請求項１に記載の推定システム。
【請求項３】
変換情報記憶手段は、変換情報として、複数の変換係数列の候補の情報を記憶する変換情報記憶手段を備え、
変換パラメータ推定手段は、前記変換情報によって示される変換係数列の候補を用いて仮の積算発電量を算出した結果に基づいて変換係数列を決定し、決定された変換係数列または該変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定する
請求項１または請求項２に記載の推定システム。
【請求項４】
第１の時間単位の日射量は、直達日射量と散乱日射量とに分かれており、
変換情報は、前記直達日射量に適用する第１の変換係数列の候補の情報と、前記散乱日射量に適用する第２の変換係数列の候補の情報とを含み、
変換パラメータ推定手段は、前記変換情報によって示される第１の変換係数列の候補と第２の変換係数列の候補とを少なくとも用いて仮の積算発電量を算出した結果に基づいて、１組の変換係数列を決定し、決定された変換係数列の組または該組の各変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定する
請求項３に記載の推定システム。
【請求項５】
変換情報は、複数の方位角および複数の傾斜角の組み合わせにおける水平面日射量と傾斜面日射量との比率が少なくとも表現された変換係数列の候補の情報を含む
請求項３または請求項４に記載の推定システム。
【請求項６】
第２の時間単位は、太陽赤緯の周期性における一周期以上であり、
変換情報は、太陽赤緯と同じの周期性をもった変換係数列の候補の情報を含む
請求項３から請求項５のうちのいずれか１項に記載の推定システム。
【請求項７】
変換情報は、一定値の変換係数列の候補の情報を含む
請求項３から請求項６のうちのいずれか１項に記載の推定システム。
【請求項８】
変換情報は、変換係数列の各候補の生成に用いた発電システムに関するパラメータの情報を含む
請求項３から請求項７のうちのいずれか１項に記載の推定システム。
【請求項９】
変換パラメータ推定手段は、変換情報によって示される変換係数列の候補の中から選択した２以上の変換係数列の候補を掛け合わせて新たな変換係数列の候補を生成し、生成された新たな変換候補列の候補を含む変換係数列の候補を用いて仮の積算発電量を算出する
請求項３から請求項８のうちのいずれか１項に記載の推定システム。
【請求項１０】
変換情報は、変換係数列の候補の生成に用いた発電システムに関するパラメータの種類別に変換係数列の候補が登録される複数の変換係数列候補テーブルを含み、
変換パラメータ推定手段は、変換係数列候補テーブルの各々から１または複数の変換係数列の候補を選択し、それらを掛け合わせて新たな変換係数列の候補を生成する
請求項９に記載の推定システム。
【請求項１１】
変換情報記憶手段は、変換情報として、所定の数理計画モデルに与える制約条件の情報を記憶し、
変換パラメータ推定手段は、前記変換情報によって示される制約条件を用いて前記所定の数理計画モデルを解くことにより変換係数列を決定し、決定された変換係数列または該変換係数列を決定した際の制約条件が示す発電システムに関するパラメータを、変換パラメータとして推定する
請求項２に記載の推定システム。
【請求項１２】
変換パラメータと、対象発電システムの設置位置または対象地域における任意の期間中の第１の時間単位の日射量とに基づいて、前記任意の期間中の前記第１の時間単位の発電量を推定する発電量推定手段を備えた
請求項１から請求項１１のうちのいずれか１項に記載の推定システム。
【請求項１３】
位置情報と、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、前記第１の時間単位の発電量を求める推定期間に関する期間情報とを少なくとも入力する入力手段と、
前記位置情報で示される位置における前記実績発電量を算出した期間中の前記第１の時間単位の日射量の系列である日射量列、および前記位置情報で示される位置における前記推定期間中の前記第１の時間単位の日射量を推定する日射量推定手段と、
前記実績発電量と、前記実績発電量を算出した期間中の前記日射量列とに基づいて、前記第２の時間単位中の前記第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する変換パラメータ推定手段と、
前記変換パラメータと、前記推定期間中の前記第１の時間単位の日射量とに基づいて、前記推定期間中の前記第１の時間単位の発電量を推定する発電量推定手段と、
推定結果を出力する出力手段とを備えた
ことを特徴とする推定システム。
【請求項１４】
変換パラメータ推定手段は、新たな実績発電量と、前記実績発電量を算出した期間中の前記第１の時間単位の日射量の系列である日射量列とが入力される度に、変換パラメータを推定する
請求項１から請求項１３のうちのいずれか１項に記載の推定システム。
【請求項１５】
変換パラメータ推定手段は、気象情報が変化する所定の期間ごとに、変換パラメータを推定する
請求項１から請求項１３のうちのいずれか１項に記載の推定システム。
【請求項１６】
変換パラメータ推定手段は、変換パラメータにより推定された期間中の第１の時間単位の発電量の積算値と、該期間中の実績発電量との差から算出される予測精度が所定の閾値を下回った場合に、前記変換パラメータを推定しなおす
請求項１から請求項１３のうちのいずれか１項に記載の推定システム。
【請求項１７】
対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、前記対象発電システムの設置位置または前記対象地域における、前記実績発電量を算出した期間中の前記第１の時間単位の日射量の系列である日射量列とが入力されると、予め記憶されている変換情報であって、前記第２の時間単位中の前記第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の前記第１の時間単位の発電量の積算値である積算発電量が前記実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報である変換情報によって示される変換係数列の候補または制約条件を用いて、前記第２の時間単位中の前記第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する
ことを特徴とする推定方法。
【請求項１８】
コンピュータに、
対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、前記対象発電システムの設置位置または前記対象地域における、前記実績発電量を算出した期間中の前記第１の時間単位の日射量の系列である日射量列とが入力されると、予め記憶されている変換情報であって、前記第２の時間単位中の前記第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の前記第１の時間単位の発電量の積算値である積算発電量が前記実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報である変換情報によって示される変換係数列の候補または制約条件を用いて、前記第２の時間単位中の前記第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する処理
を実行させるための推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、発電量を推定するための推定システム、推定方法および推定プログラムに関する。
【背景技術】
【０００２】
太陽光発電などの再生可能エネルギーは気候や設置場所の影響を受けやすく、出力量（例えば、発電量）は時間によって大きく変動する。
【０００３】
また、近年、家庭用のソーラーパネルなど小規模の発電システムが増えている。このような状況に鑑みて、電力会社などが管理する大規模発電システムだけでなく、これら小規模の発電システムを含めて、発電量を詳細に把握したいとの要望がある。しかし、小規模の発電システムを含む発電システムの発電効率や設置条件等を全て管理するのは容易ではない。
【０００４】
発電量の算出技術に関して、例えば特許文献１、２がある。例えば、特許文献１には、時間帯別最大発電量と時間帯別日射強度との対応関係を抽出し、該対応関係に基づいて、日射強度予報値からその時間帯の太陽光発電量を予測する太陽光発電量予測装置が記載されている。
【０００５】
また、例えば、特許文献２には、日射強度を発電量に変換する発電係数を算出し、算出された発電係数を用いて日射強度の予測値から発電量の予測値を算出する予測システムが記載されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００７−１７３６５７号公報
【特許文献２】特開２０１５−５６４１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、特許文献１に記載の太陽光発電予測装置は、最大発電量と日射強度とを対応づける際に、日射強度の実績値と同じ時間単位の太陽光発電の実績値が取得されることを前提としている。また、特許文献２に記載の予測システムは、日射強度の予測値から発電量に変換するための発電係数を最適化するために、やはり日射量実績値と同じ時間単位の発電量実績値を必要としている。
【０００８】
上述したように、全ての発電システムから所望する時間単位の発電量の実績値を取得できるとは限らない。例えば、発電システムから比較的短時間単位の発電量を漏れなく得ようとした場合、その発電システムに、該時間単位の発電量を計測するセンサと、計測結果をサーバ等に通知する通知手段とが必要となる。しかし、小規模発電システムを含む発電システムの全てに対して、そのようなセンサと通知手段とを設けるのは非常に困難である。
【０００９】
なお、特許文献２に記載の方法によれば、比較的短時間単位の発電量が得られなくても、発電システムの発電効率や設置条件などを与えることにより、日射強度の予測値から発電量に変換できる。しかし、上記の発電量の実績値と同様、小規模発電システムを含む発電システムの全ての発電特性や設置条件等の情報を管理するのは非常に困難である。
【００１０】
そこで、本発明は、取得可能な発電システムの情報が限られている場合であっても、より詳細な発電量を推定可能な推定システム、推定方法および推定プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明による推定システムは、対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とに基づいて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する変換パラメータ推定手段を備えたことを特徴とする。
【００１２】
また、本発明による推定システムは、位置情報と、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、第１の時間単位の発電量を求める推定期間に関する期間情報とを少なくとも入力する入力手段と、位置情報で示される位置における実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列、および位置情報で示される位置における推定期間中の第１の時間単位の日射量を推定する日射量推定手段と、実績発電量と、実績発電量を算出した期間中の日射量列とに基づいて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する変換パラメータ推定手段と、変換パラメータと、推定期間中の第１の時間単位の日射量とに基づいて、推定期間中の第１の時間単位の発電量を推定する発電量推定手段と、推定結果を出力する出力手段とを備えたことを特徴とする。
【００１３】
また、本発明による推定方法は、対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とが入力されると、予め記憶されている変換情報であって、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報である変換情報によって示される変換係数列の候補または制約条件を用いて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定することを特徴とする。
【００１４】
また、本発明による推定プログラムは、コンピュータに、対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とが入力されると、予め記憶されている変換情報であって、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報である変換情報によって示される変換係数列の候補または制約条件を用いて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する処理を実行させることを特徴とする。
【発明の効果】
【００１５】
本発明によれば、取得可能な発電システムの情報が限られている場合であっても、より詳細な発電量を推定できる。
【図面の簡単な説明】
【００１６】
【図１】第１の実施形態の推定システムの例を示す構成図である。
【図２】第１の時間単位と第２の時間単位の関係を示す説明図である。
【図３】第１の実施形態の推定システムの動作の一例を示すフローチャートである。
【図４】変換パラメータの決定方法の一例を示す説明図である。
【図５】変換パラメータの決定方法の他の例を示す説明図である。
【図６】変換係数列候補テーブルの例を示す説明図である。
【図７】変換係数列候補の例を示す説明図である。
【図８】変換係数列候補テーブルの他の例を示す説明図である。
【図９】変換係数列候補テーブルの他の例を示す説明図である。
【図１０】第２の実施形態の推定システムの例を示すブロック図である。
【図１１】第２の実施形態の推定システムの動作の一例を示すフローチャートである。
【図１２】第３の実施形態の推定システムの例を示すブロック図である。
【図１３】第３の実施形態の推定システムの動作の一例を示すフローチャートである。
【図１４】第４の実施形態の推定システムの例を示すブロック図である。
【図１５】第４の実施形態の推定システムの動作の一例を示すフローチャートである。
【図１６】本発明の各実施形態にかかるコンピュータの構成例を示す概略ブロック図である。
【図１７】本発明の概要を示すブロック図である。
【図１８】本発明による推定システムの他の構成を示すブロック図である。
【図１９】本発明による推定システムの他の構成を示すブロック図である。
【発明を実施するための形態】
【００１７】
実施形態１．
以下、本発明の実施形態を図面を参照して説明する。図１は、第１の実施形態の推定システムの例を示す構成図である。図１に示す推定システム１は、日射量列取得手段１１と、実績発電量取得手段１２と、変換情報記憶手段１３と、変換パラメータ推定手段１４とを備える。
【００１８】
日射量列取得手段１１は、対象発電システムの設置位置における所定の評価期間中の第１の時間単位の日射量の系列である日射量列を取得する。取得される日射量の種別は特に問わない。日射量は、例えば、平均全天日射量であってもよいし、直達日射量と散乱日射量とに分かれていてもよい。
【００１９】
実績発電量取得手段１２は、対象発電システムの上記の評価期間における発電量の実績値である実績発電量を取得する。
【００２０】
ここで、評価期間は、実績発電量が取得可能な過去の期間であればよい。また、第１の時間単位は、評価期間の長さである第２の時間単位よりも短い時間単位とする。第１の時間単位は、例えば、推定したい発電量の時間間隔が設定されてもよい。なお、第２の時間単位は、２以上の第１の時間単位を含む。以下、実績発電量を、第２の時間単位の発電量と表現する場合がある。
【００２１】
図２は、第１の時間単位と第２の時間単位の関係を示す説明図である。図２（ａ）〜図２（ｂ）に示すように、評価期間は、第２の時間単位の時間長をもつ時系列上の所定の期間であればよい。また、図２（ａ）は、第２の時間単位が、一定の第１の時間単位を含む例である。なお、図２（ｂ）に示すように、第１の時間単位は、一定でなくてもよい。すなわち、第２の時間単位は、異なる長さの第１の時間単位を含んでいてもよい。これは、例えば、日中の日射量を３０分単位で取得して、夜間の日射量をそれよりも大きい時間単位で取得する場合等に相当する。また、第２の時間単位の大きさは固定でなくてもよい。これは、例えば、第２の時間単位を１ヶ月とした場合に、ある月では３０日とされ、ある月では３１日とされる場合等に相当する。なお、図２において、ｔは、第２の時間単位に含まれる第１の時間単位を識別する変数を表している。
【００２２】
また、対象発電システムの設置位置における日射量は、例えば、対象発電システムが実際に受光した日射量（実績値）であってもよいし、対象発電システムが設置されている所定の地域における平均的な日射量であってもよい。
【００２３】
日射量の取得方法としては、対象発電システムの設置位置もしくはその近傍に設けられた日射計を用いる方法や、該設置位置を含む所定の地域の気象に関する予報値を含むＧＰＶ（Grid Point Value）データから推定する方法や、該設置位置を含む所定の地域の上空を撮影した衛星画像から推定する方法などが挙げられる。
【００２４】
実績発電量の取得方法としては、例えば、売電契約を行って、対象発電システムで計測される値を取得する方法であってもよい。なお、上記の例はあくまで一例であって、日射量および実績発電量の取得方法はこの限りではない。
【００２５】
変換情報記憶手段１３は、日射量列に含まれる各日射量を発電量に変換するための変換パラメータを決定するための変換情報を予め記憶する。ここで、変換パラメータは、例えば、日射量列に含まれる各日射量を発電量に変換する変換式に用いられる変換係数の系列である変換係数列であってもよいし、そのような変換係数列を決定するための所定の数理計画モデルに与える制約条件であってもよい。また、変換パラメータは、変換係数列の生成に用いた対象発電システムに関するパラメータや制約条件が意味する対象発電システムに関するパラメータ、例えば対象発電システムの発電特性（効率）や発電容量や設置条件（例えば、傾斜角や方位角や太陽方向に対する遮蔽物の有無）等またはその時間変化を示す情報であってもよい。
【００２６】
変換情報は、例えば、変換係数列の候補の情報であってもよいし、制約条件の情報であってもよい。変換係数列の候補の情報としては、候補そのものであってもよいし、候補とともにその適用方法（優先順位や組み合わせルールなど）を示す情報が含まれていてもよい。また、制約条件の情報としては、制約条件とされる対象発電システムのパラメータ（発電効率や設置条件）の候補や、値の範囲や、値の精度や、値の時間変化の情報であってもよいし、それらの情報とともにその適用方法を示す情報が含まれていてもよい。
【００２７】
変換パラメータ推定手段１４は、評価期間中の日射量列および実績発電量と、変換情報記憶手段１３に記憶されている変換情報とに基づいて、変換パラメータを決定する。変換パラメータ推定手段１４は、例えば、日射量列に含まれる各日射量を発電量に変換する変換式に、変換情報によって示される変換係数列の候補を適用して仮の積算発電量を算出し、その結果に基づいて変換係数列を決定することにより、変換パラメータを推定してもよい。そのような場合において、変換パラメータ推定手段１４は、決定された変換係数列または該変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定してもよい。
【００２８】
また、変換パラメータ推定手段１４は、例えば、前記変換情報によって示される制約条件を用いて、所定の数理計画モデルを解いて得られた変換係数列または該変換係数列を決定した際の制約条件が示す発電システムに関するパラメータを、変換パラメータとして推定してもよい。変換パラメータ推定手段１４は、より具体的には、変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルを解くことによって、変換係数列を決定してもよい。そのような場合において、変換パラメータ推定手段１４は、決定された変換係数列または該変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定してもよい。
【００２９】
次に、本実施形態の動作について説明する。図３は、本実施形態の推定システムの動作の一例を示すフローチャートである。図３に示す例では、まず、日射量列取得手段１１および実績発電量取得手段１２が、評価期間中の日射量列および実績発電量を取得する（ステップＳ１０１）。
【００３０】
次に、変換パラメータ推定手段１４が、ステップＳ１０１で取得された評価期間中の日射量列および実績発電量と、変換情報記憶手段１３に記憶されている変換情報とに基づいて、変換パラメータを決定する（ステップＳ１０２）。
【００３１】
変換パラメータ推定手段１４は、決定した変換パラメータを外部の装置に出力してもよいし、変換情報記憶手段１３に記憶させてもよい。
【００３２】
次に、変換情報および変換パラメータの決定方法についてより詳細に説明する。以下でも、第２の時間単位に含まれる第１の時間単位を識別する変数として、時間ｔ＝１，・・・，ｎを用いる。なお、ｎは第２の時間単位に含まれる第１の時間単位数（最大値）を表す。なお、最大値としているのは、第２の時間単位の時間長が変動する場合を考慮した為である。
【００３３】
［第１の例］
第１の例は、変換情報として、日射量列が含む日射量の各々に適用させる変換係数を含む変換係数列の候補である変換係数列候補が複数登録された変換係数列候補テーブルが与えられる例である。
【００３４】
図４は、本例における変換パラメータの決定方法の一例を示す説明図である。本例の変換パラメータ推定手段１４は、変換係数列候補テーブルから変換係数列候補を１つ選択し、変換式に、取得された日射量列に含まれる日射量と、選択した変換係数列候補に含まれる該日射量の時間ｔに対応する変換係数とを順次適用して、仮の積算発電量を算出する。ここで、仮の積算発電量は、変換式を用いて算出された第１の時間単位の発電量の積算値である。変換パラメータ推定手段１４は、変換係数列候補テーブルの中から、取得された評価期間中の実績発電量に最も近い仮の積算発電量が得られた変換係数列を決定して、変換パラメータを推定する。
【００３５】
本例では、日射量列に含まれる日射量の最大数と同じ数の変換係数を含む変換係数列候補が変換係数列候補テーブルに登録されているものとする。図４に示す例では、取得された日射量列を、Ｉ＝｛ｉ１，・・・，ｉｎ｝で表し、変換係数列候補テーブル内の任意の変換係数列候補を、Ｃｊ＝｛ｃｊ１，・・・，ｃｊｎ｝で表している。ここで、ｊは変換係数列候補テーブル内の変換係数列候補を識別する変数である。
【００３６】
変換パラメータ推定手段１４は、例えば、変換係数列候補テーブルから変換係数列候補を１つずつ選び、選んだ変換係数列候補Ｃｊと日射量列Ｉとを用いて、仮の積算発電量ＰＭｊを、以下の式（１）を用いて算出する。式（１）において、ｃｊｔは、変換係数列候補Ｃｊに含まれる時間ｔに対応する変換係数を表す。また、ｉｔは、日射量列Ｉに含まれるち時間ｔに対応する日射量を表す。
【００３７】
ＰＭｊ＝Σｔ（ｃｊｔ×ｉｔ）  ・・・（１）
【００３８】
変換パラメータ推定手段１４は、全ての変換係数列候補Ｃｊについて仮の積算発電量を求め、得られた仮の積算発電量と、取得した実績発電量（第２の時間単位の発電量の実績値）との差が最も小さくなる変換係数列候補を、対象発電システムにおける変換係数列に決定すればよい。
【００３９】
変換係数列候補テーブルには、例えば、対象発電システムの発電効率や、発電容量や、対象発電システムの設置角度による日射量減衰率や、対象発電システムの周辺環境による日射量減衰率など、日射量から発電量に変換する際に影響を与える発電システムに関するパラメータの値やその値の時間変化が様々なパターンで表現された変換係数列候補が含まれていることが好ましい。
【００４０】
このような変換係数列候補を予め用意しておくことにより、対象発電システムの発電特性や設置条件等が分からなくても、それらの影響が反映された変換パラメータを評価期間中の日射量列と実績発電量とから推定することができる。変換パラメータが推定できれば、対象発電システムの任意の期間における第１の時間単位の日射量から当該期間における第１の時間単位の発電量を推定することができる。
【００４１】
例えば、変換パラメータとして、変換係数列Ｃ１＝｛ｃ１１，・・・，ｃ１ｎ｝が推定されたとする。その場合、任意の期間上に設定される第１の時間単位の発電量ｐは、該期間の第１の時間単位の日射量ｉと、変換係数列に含まれるその期間が対応している第２の時間単位中の第１の時間単位を識別する時間ｔに対応した変換係数ｃ１ｔとを用いて算出できる（式（２）参照）。
【００４２】
ｐ＝ｃ１ｔ×ｉ  ・・・（２）
【００４３】
時間ｔは、求めたい期間が対応する第２の時間単位中における第１の時間単位である。例えば、第２の時間単位が１ヶ月、第１の時間単位が１時間、第２の時間単位の起算日が各月の１日０時であったとする。そのような場合、例えば、求めたい期間がある月の１日の０時〜１時までの１時間であれば、ｔ＝１とし、求めたい期間がある月の１日の１時〜２時までの１時間であれば、ｔ＝２とすればよい。
【００４４】
［第２の例］
第２の例は、変換情報として、２種類の変換係数列候補テーブルが与えられる例である。図５は、本例における変換パラメータの決定方法の一例を示す説明図である。本例は、日射量が直達日射量と散乱日射量とに分かれて取得される場合の例である。そのような場合、直達日射量と散乱日射量とに適用させる変換係数列候補が別々に登録された変換係数列候補テーブルが予め用意されてもよい。図５に示す例では、直達日射量列に適用させる第１の変換係数列候補が登録された第１の変換係数列候補テーブルと、散乱日射量列に適用させる第１の変換係数列候補が登録された第２の変換係数列候補テーブルとが用意されている。なお、図示省略してるが、いずれにも適用可能な変換係数列候補が登録される第３の変換係数列候補テーブルがさらに用意されてもよい。
【００４５】
本例では、変換パラメータ推定手段１４は、第１の変換係数列の候補と第２の変換係数列の候補とを少なくとも用いて仮の積算発電量を算出し、その結果に基づいて、１組の変換係数列を決定してもよい。そのような場合、変換パラメータ推定手段１４は、決定された変換係数列の組または該組の各変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定してもよい。
【００４６】
例えば、変換パラメータ推定手段１４は、日射量列に含まれる各時間ｔの直達日射量に、第１の変換係数列候補テーブルから選択した第１の変換係数列候補に含まれる時間ｔに対応した変換係数を掛けたものと、日射量列に含まれる各時間ｔの散乱日射量に、第２の変換係数列候補テーブルから選択した第２の変換係数列候補に含まれる時間ｔに対応した変換係数を掛けたものとの総和をとり、仮の積算発電量としてもよい。
【００４７】
例えば、第１の変換件数列候補テーブルから選択された第１の変換係数列候補がＣａｊ＝｛ｃａｊ１，・・・，ｃａｊｎ｝であり、第２の変換件数列候補テーブルから選択された第２の変換係数列候補がＣｂｋ＝｛ｃｂｋ１，・・・，ｃｂｋｎ｝であったとする。ここで、ｊ，ｋは、その変換係数列候補テーブル内の変換係数列候補を識別する変数である。また、取得された日射量列がＩ＝｛｛ｄｉ１，ｓｉ１｝，・・・，｛ｄｉｎ，ｓｉｎ｝｝であったとする。以下、日射量列のうち、直達日射量ｄｉｊのみからなる直達日射量列をＤＩ＝｛ｄｉ１，・・・，ｄｉｎ｝、散乱日射量ｓｉｊのみからなる散乱日射量列をＳＩ＝｛ｓｉ１，・・・，ｓｉｎ｝で表す。このような場合において、変換パラメータ推定手段１４は、例えば、仮の積算発電量ＰＭｊｋを、以下の式（３）を用いて算出してもよい。なお、式（３）において、ｃａｊｔ，ｃｂｋｔは、変換係数列候補Ｃａｊ，Ｃｂｋに含まれる変換係数のうち時間ｔに対応する変換係数を表す。また、ｄｉｔ，ｓｉｔは、時間ｔに対応する直達日射量，散乱日射量を表す。
【００４８】
ＰＭｊｋ＝Σｔ（（ｃａｊｔ×ｄｉｔ）＋（ｃｂｋｔ×ｓｉｔ））  ・・・（３）
【００４９】
変換パラメータ推定手段１４は、変換係数列候補Ｃａｊ，Ｃｂｋの全ての組み合わせについて仮の積算発電量を求め、得られた仮の積算発電量と、取得した実績発電量との差が最も小さくなる変換係数列候補の組を決定してもよい。
【００５０】
第１の変換係数列候補テーブルには、例えば、対象発電システムの発電特性や、発電容量や、対象発電システムの設置条件による直達日射量減衰率や、周辺環境による直達日射量減衰率など、直達日射量から発電量に変換する際に影響を与える発電システムに関するパラメータの値やその値の時間変化が様々なパターンで表現された変換係数列候補が含まれていることが好ましい。
【００５１】
また、第２の変換係数列候補テーブルには、例えば、対象発電システムの発電特性や、発電容量や、対象発電システムの設置条件による散乱日射量減衰率や、周辺環境による散乱日射量減衰率など、散乱日射量から発電量に変換する際に影響を与える発電システムに関するパラメータの値やその値の時間変化が様々なパターンで表現された変換係数列候補が含まれていることが好ましい。
【００５２】
なお、別途いずれにも適用可能な変換係数列候補がある場合には、第１の変換係数列候補と第２の変化係数列候補のそれぞれにその変換係数列候補を掛け合わせたものを、新たな第１の変換係数列候補と第２の変化係数列候補としてもよい。なお、変換係数列候補の掛け合わせにより新たな変換係数候補を生成する方法は、後述の第４の例にも示されている。
【００５３】
これにより、対象発電システムの発電特性や設置条件等が分からなくても、それらの影響が反映された変換パラメータを評価期間中の日射量列（直達日射量列と散乱日射量列）と実績発電量とから推定することができる。なお、直達日射量と散乱日射量とに適用する変換係数列を別々に推定することにより、例えば、設置角度や周囲の環境による影の影響といった設置条件による日射量減衰率など、直達日射量と散乱日射量とで値が異なる場合にも対応できるので、より精度よく変換することができる。
【００５４】
例えば、選択された変換係数列の組がＣａ１＝｛ｃａ１１，・・・，ｃａ１ｎ｝とＣｂ２＝｛ｃｂ２１，・・・，ｃｂ２ｎ｝であったとする。その場合、任意の期間上に設定される第１の時間単位の発電量ｐは、該期間の第１の時間単位の直達日射量ｄｉおよび散乱日射量ｓｉと、各変換係数列に含まれるその期間が対応している第２の時間単位中の第１の時間単位を識別する時間ｔに対応した変換係数ｃａ１ｔおよびｃｂ２ｔとを用いて算出できる（式（４）参照）。
【００５５】
ｐ＝（ｃａ１ｔ×ｄｉ）＋（ｃｂ２ｔ×ｓｉ）  ・・・（４）
【００５６】
［第３の例］
第３の例は、変換情報として、複数の方位角および複数の傾斜角の組み合わせにおける水平面日射量と傾斜面日射量との比率が少なくとも表現された変換係数列の候補が登録された変換係数列候補テーブルが与えられる例である。なお、本例では、第１の時間単位の日射量として、水平面日射量が取得されることを想定している。
【００５７】
図６は、本例の変換係数列候補テーブルの例を示す説明図である。本例の変換係数列候補テーブルには、発電システムがある方位角およびある傾斜角をとった場合の、該発電システムにおける水平面日射量と傾斜面日射量の比率が少なくとも表現された変換係数からなる変換系数列候補が登録される。なお、図６に示すように、各変換係数列候補に対応づけて、その変換係数列候補の生成に用いたパラメータの情報（本例では、方位角および傾斜角）が付与されていてもよい。
【００５８】
図中において、θは発電システムの方位角を表し、φは発電システムの傾斜角を表している。このとき、各変換係数列候補に含まれる変換係数は、その方位角および傾斜角での水面日射量と傾斜面日射量の比率を表す、０〜１の値であってもよい。なお、変換係数は、水平面日射量と傾斜面日射量の比率を表す値に、他の比率（例えば、発電効率や、日射量減衰率）を表す値を掛け合わせたものであってもよい。そのような場合、その変換係数列候補の生成に用いた方位角および傾斜角の情報に加えて、上記の他の比率の計算に用いた情報（発電効率や、時間帯別の遮蔽物の有無等）が付与されていてもよい。
【００５９】
水面日射量と傾斜面日射量の比率は、例えば、発電システムが対象とされた方位角および傾斜角で設置されていると仮定して、各時間帯（ｔ＝１，・・・，ｔ）において晴天だった場合の該発電システムの単位面積当たりで受光される日射量予測値（傾斜面予測値）の、水平面日射量予測値に対する減衰率、またはそれらを任意の定数で割ったものであってもよい。
【００６０】
本例の変換パラメータの決定方法は、第１の例および第２の例と同様でよい。
【００６１】
なお、第２の例のように、取得される日射量が直達日射量と散乱日射量とに分かれている場合には、第２の例と同様に、水平面直達日射量に適用する第１の変換係数列候補テーブルと、水平面散乱日射量に適用する第２の変換係数列候補テーブルの２種類を用意すればよい。
【００６２】
これにより、対象発電システムの発電特性や設置条件等が分からなくても、それらの影響が反映された変換パラメータを評価期間中の日射量列と実績発電量とから推定することができる。なお、発電システムにおける水平面日射量と傾斜面日射量の比率のみが表現された変換係数からなる変換係数列候補の場合、該変換係数列候補は、水平面日射量を傾斜面日射量に変換する用途に使用できる。そのような場合には、後述の第４の例に示すように、他の変換係数列候補と掛け合わせて適用すればよい。それにより、様々な条件による影響を考慮することができる。
【００６３】
［第４の例］
第４の例は、変換パラメータ推定手段１４が、１つまたは複数の変換係数列候補テーブルから選択された少なくとも２以上の変換係数列候補を掛け合わせたものを、新たな変換係数列候補として用いる例である。
【００６４】
本例では、変換情報として、発電システムのパラメータの種類別に変換係数列候補テーブルが用意されていてもよいし、１つの変換係数列候補テーブルに、２以上のパラメータの様々な組み合わせに対応した変換係数列候補が含まれていてもよい。
【００６５】
図７は、本例の変換係数列候補テーブルの例および変換係数列候補の選択例を示す説明図である。図７（ａ）に示すように、変換パラメータ推定手段１４は、複数の変換係数列候補テーブルからそれぞれ１つの変換係数列候補を選択し、それらを掛け合わせて新たな変換係数列候補としてもよい。なお、変換係数列候補テーブルが３以上ある場合は、そのうちの少なくとも２以上の変換係数列候補テーブルからそれぞれ１つの変換係数列候補を選択し、それらを掛け合わせて新たな変換係数列候補としてもよい。
【００６６】
掛け合わせる変換係数列候補の組み合わせは、特に限定されないが、各変換係数列候補テーブルに対して予め定められたルールに従って、掛け合わせる変換係数列候補の組み合わせを選択してもよい。
【００６７】
例えば、変換件数列候補テーブルＡから選択された変換係数列候補がＡｊ＝｛ａｊ１，・・・，ａｊｎ｝であり、変換件数列候補テーブルＢから選択された変換係数列候補がＢｋ＝｛ｂｋ１，・・・，ｂｋｎ｝であったとする。ここで、ｊ，ｋは、それぞれ変換係数列候補テーブル内の変換係数列候補を識別する変数である。また、取得された日射量列をＩ＝｛ｉ１，・・・，ｉｎ｝とする。このような場合、変換パラメータ推定手段１４は、例えば、仮の積算発電量ＰＭａｊｂｋを、以下の式（５）を用いて算出する。式（５）において、ａｊｔ，ｂｋｔは、変換係数列候補Ａｊ，Ｂｋに含まれる変換係数のうち時間ｔに対応する変換係数を表す。
【００６８】
ＰＭａｊｂｋ＝Σｔ（（ａｊｔｂｋｔ）×ｉｔ）  ・・・（５）
【００６９】
また、図７（ｂ）に示すように、変換パラメータ推定手段１４は、１つの変換係数列候補テーブルから複数の変換係数列候補を選択し、それらを掛け合わせて新たな変換係数列候補とすることも可能である。例えば、太陽方向の遮蔽物の有無に対応した変換係数列候補テーブルを考える。当該変換係数列候補テーブルには、時間帯別に、太陽方向の遮蔽物の有無による日射量減衰率が表現された変換係数からなる変換係数列候補が登録される。各変換係数列候補は、例えば、発電システムから見てある１つの方位に１つの遮蔽物があると仮定したときの日射量減衰率の時間変化が表された変換係数列候補であってもよい。時間帯別の太陽方向の遮蔽物の有無のパターンを全て表現しようとすると登録数が非常に多くなるが、例えば、上記のように条件を限定して方位別に変換係数列候補を用意する方法であれば、登録数を削減できる。変換パラメータ推定手段１４は、登録されている変換係数列候補を複数組み合わせれば、２以上の遮蔽物に対応した変化係数列候補を新たに生成できる。
【００７０】
例えば、１つの変換係数列候補テーブルから選択された変換係数列候補が、Ｃ１＝｛ｃ１１，・・・，ｃ１ｎ｝とＣ２＝｛ｃ２１，・・・，ｃ２ｎ｝であったとする。取得された日射量列をＩ＝｛ｉ１，・・・，ｉｎ｝とすると、変換パラメータ推定手段１４は、例えば、仮の積算発電量ＰＭｃ１ｃ２を、以下の式（６）を用いて算出する。式（６）において、ｃ１ｔ，ｃ２ｔは、変換係数列候補Ｃ１，Ｃ２に含まれる変換係数のうち時間ｔに対応する変換係数を表す。なお、１つの変換係数列候補テーブルから選択される変換係数列候補の数は２つに限定されない。
【００７１】
ＰＭｃ１ｃ２＝Σｔ（（ｃ１ｔｃ２ｔ）×ｉｔ）  ・・・（６）
【００７２】
なお、図７（ａ）に示す例と、図７（ｂ）に示す例とを組み合わせることも可能である。すなわち、あるテーブルからは１つの変換係数列候補を選択し、他のあるテーブルからは２以上の変換係数列候補を選択して、それらを掛け合わせて新たな変換係数列候補としてもよい。そのような場合、各変換係数列候補テーブルに、選択可能な変換係数列候補の数や、同じ変換係数列候補を重複して選択可能かといったルールを定めておいてもよい。
【００７３】
［第５の例］
第５の例は、変換係数列候補が、２４時間や１日といった太陽赤緯と同様の周期性をもつ値で表現される例である。なお、本例は、第２の時間単位が太陽赤緯の周期性における一周期以上であることを前提としている。本例では、変換情報として、そのような周期性をもつ値で表現された変換係数列候補が登録された変換係数列候補テーブルが予め用意される。なお、各変換係数列候補は、第２の時間単位中の第１の時間単位数に相当するｎ個の変換係数を含んでいてもよいし、ｎ個より少ない一周期分の変換係数のみを含んでいてもよい。後者の場合には、変換パラメータ推定手段１４は、日射量に変換係数列を適用する際に、該変換係数列が、登録されている変換係数の繰り返し構造であるとして、周期毎に変換係数列の候補を展開して適用する。
【００７４】
図８は、本例の変換係数列候補テーブルの例を示す説明図である。なお、図８（ａ）は、ｎ個の変換係数を含む変換係数列候補が登録された変換係数列候補テーブルの例であり、図８（ｂ）は、一周期分の変換係数のみを含む変換係数列候補が登録された変換係数列候補テーブルの例である。なお、図８に示す例では、各変換係数が０〜１の値をとる例が示されているが、変換係数は０〜１の値に限られない。
【００７５】
このような周期性を利用することにより、テーブル内における変換係数列候補の組み合わせを限定することができる。例えば、一ヶ月分の日射量に対して、第３の例のような様々な発電システムの方位角および傾斜角に対応した日射量減衰率に相当する値を変換係数として含む変換係数列候補を用意することを考える。ここで、第１の時間単位は３０分、第２の時間単位は１ヶ月とする。すると、ｎ＝３１（日）×２４（時間）／０．５＝１４８８である。また、方位角のパターン数をｐθ、傾斜角のパターン数をｐφとする。なお、両者を合わせたパターン数はｐθ×ｐφである。
【００７６】
例えば、周期性を考慮しないで変換係数列候補を生成した場合、ｎ×ｐθ×ｐφの変数を独立に扱う必要がある。日周期性があることを考慮すると、変換係数列候補に含まれる変換係数の数が１日分（ｎ’＝２４（時間）／０．５＝４８）で済むことから、変換係数列候補テーブルにおいて考慮する変数はｎ’×ｐθ×ｐφ個で済む。
【００７７】
このように周期性を制約条件とすることにより、変換係数列候補テーブルにおいて考慮する変数の数を減らすことができるので、テーブル内を探索するのが容易になる。また、周期性を制約条件とすることにより、第１の時間単位の各日射量の計測（推定誤差）に対して過剰に変換係数をフィッティングしてしまう過学習を防ぐことができるので、決定される変換係数列の精度向上が期待できる。なお、本例においても他の変換係数列候補との掛け合わせは有効である。以下、特に断りがない場合、本例以外のいずれも変換係数列候補も他の変換係数候補との掛け合わせも有効とする。
【００７８】
［第６の例］
第６の例は、変換係数列候補が一定値の変換係数の系列である例である。図９は、本例の変換係数列候補テーブルの例を示す説明図である。本例の変換係数列候補テーブルには、図９（ａ）に示すように、すべて一定値の変換係数の系列である変換係数列候補が登録されている。なお、１つの変換係数列候補内における変換係数は一定値であればよく、必ずしも０〜１の範囲でなくてもよい。例えば、発電容量に対応する変換係数の場合には、１以上の値もとり得る。
【００７９】
また、図９（ｂ）に示すように、各変換係数列候補は、一定値とされる１つの値を表した１つの変換係数のみを含んでいてもよい。この場合、変換パラメータ推定手段１４は、各変換係数列候補を、登録されている１つの変換係数の繰り返し構造とみなせばよい。
【００８０】
また、発電容量に応じた変換を行うスケール調整に関して、予めスケール調整用に変換係数列候補を用意せずに、例えば、取得された日射量列および実績発電量に応じて、（実績発電量／日射量列の積算値）±αの範囲内に収まるような一定値を持った変換係数列候補を新たに生成してもよい。ここで、αは特に限定されないが、例えば、（実績発電量／日射量列の積算値）としてもよい。
【００８１】
一定値の変換係数からなる変換係数列の候補を用意することにより、変換係数列候補テーブルにおいて考慮する変数の数を減らすことができるので、テーブル内を探索するのが容易になる。
【００８２】
［第７の例］
第７の例は、変換情報として、変換係数列を決定するための所定の数理計画モデルに与える制約条件の情報が与えられる例である。
【００８３】
本例の変換パラメータ推定手段１４は、日射量列に含まれる各日射量を発電量に変換するための変換係数列を、一般的な二次計画問題を解くことによって算出する。二次計画問題では、無数に考えられる変換係数列候補すなわち変換係数の組み合わせに対して、制約条件によって取り得る候補（変換係数の組み合わせ）の限定を行い、その中から最もよいものを選び出すことができる。
【００８４】
以下、二次計画問題の一例を説明する。例えば、発電システムの周囲に鏡などの反射物や、太陽以外の光源がないと仮定すると、日射量は、太陽からの直接日射と晴天からの一様散乱日射の和によって表されると考えられる。
【００８５】
今、ある時間帯における対象発電システムの設置位置における水平面直達日射量をＤとし、その時間帯における対象発電システムの方位角および傾斜角による減衰率をζＤＳとする。すると、周辺に遮蔽物が何もないと仮定したときのその時間帯における対象発電システムの設置位置、方位角および傾斜角における傾斜面直達日射量は、ζＤＳＤで表される。また、その時間帯における該地点の太陽方向の遮蔽物による減衰率をζＤＩとする。すると、対象発電システムが実際に受光する直達日射量である実傾斜面直達日射量は、ζＤＩζＤＳＤで表される。
【００８６】
また、ある時間帯における対象発電システムが設置された地点に対する水平面散乱日射量をＳとし、その時間帯における対象発電システムの方位角および傾斜角による減衰率をζＳＳとする。すると、周辺に遮蔽物が何もないと仮定したときのその時間帯における対象発電システムの設置地位、方位角および傾斜角における傾斜面散乱日射量は、ζＳＳＳで表される。また、その時間帯における該地点の太陽方向の遮蔽物による減衰率をζＳＩとする。すると、対象発電システムが実際に受光する散乱日射量である実傾斜面散乱日射量は、ζＳＩζＳＳＳで表される。
【００８７】
また、ηを対象発電システムの発電効率とすると、単位時間あたりの対象発電システムの発電量ｐは、以下の式（７）で表される。
【００８８】
ｐ＝η（ζＤＩζＤＳＤ＋ζＳＩζＳＳＳ）  ・・・（７）
【００８９】
上記の仮定によれば、定性的に、ζＤＳ、ζＤＩおよびζＳＳは周期性があり、ζＳＩおよびηは一定という関係が成り立つ。したがって、式（７）を展開して整理すると、以下の式（８）のように、発電量ｐは、周期性の変換係数と一定値の変換係数の組み合わせで表すことができる。ただし、Ｃｃｏｎｓｔ１＝η、Ｃｃｙｃｌｅ１＝ζＤＩζＤＳ、Ｃｃｏｎｓｔ２＝ηζＳＩ、Ｃｃｙｃｌｅ２＝ζＳＳである。
【００９０】
ｐ＝Ｃｃｏｎｓｔ１Ｃｃｙｃｌｅ１Ｄ＋Ｃｃｏｎｓｔ２Ｃｃｙｃｌｅ２Ｄ  ・・・（８）
【００９１】
各時間ｔで式（８）が成り立つとすると、その時間積算も同様に成り立つ。したがって、以下の式（９）が成り立つ。
【００９２】
【数１】
2018007461.tif 000003
【００９３】
ここで、各時間ｔの直達日射量Ｄｔおよび散乱日射量Ｓｔと、各時間ｔの発電量の積算値ＰＭとが既知であれば、以下の式（１０）で表される非線形問題を解くことにより、Ｃｃｏｎｓｔ１、Ｃｃｙｃｌｅ１、Ｃｃｏｎｓｔ２、Ｃｃｙｃｌｅ２を得ることができる。
【００９４】
【数２】
2018007461.tif 000004
【００９５】
ところで、Ｃｔｃｏｎｓｔ１とＣｔｃｏｎｓｔ２とは、周期整数列と読み替えることができ、すると、以下の式（１１）で表される一般的な二次計画問題を解くことにより、変換係数列Ｃ’ｔｃｙｃｌｅ１およびＣ’ｔｃｙｃｌｅ２が得られる。なお、Ｃ’ｔｃｙｃｌｅ１＝Ｃｔｃｏｎｓｔ１Ｃｔｃｙｃｌｅ１、Ｃ’ｔｃｙｃｌｅ２＝Ｃｔｃｏｎｓｔ２Ｃｔｃｙｃｌｅ２である。
【００９６】
【数３】
2018007461.tif 000005
【００９７】
本例の変換情報は、例えば、上記の二次計画問題に与える制約条件の候補を示す情報であってもよい。より具体的には、対象発電システムのパラメータ（発電効率や設置条件）の値やその時間変化やそれに基づいて生成される変換係数列候補（ζＤＳ、ζＤＩ、ζＳＳ、ζＳＩおよびηの一周期分の値の組み合わせの候補）を示す情報であってもよい。また、変換情報は、制約条件の適用方法（例えば、各パラメータの値の範囲、精度、周期性、各組み合わせルール、優先順位等）を示す情報を含んでいてもよい。例えば、上記の各例の変換係数列候補の生成に用いた発電システムに関するパラメータの候補を、制約条件として与えることも可能である。
【００９８】
なお、制約条件の適用方法が固定である場合や、１つの対象システムに対して十分な学習データ（日射量列と実績発電量の組）が用意される場合など、制約条件を外部から与える必要がない場合には、変換情報を省略することも可能である。
【００９９】
例えば、変換パラメータ推定手段１４は、複数の評価期間についてそれぞれ二次計画問題を解き、その結果を機械学習して最終的な変換パラメータを推定してもよい。なお、変換パラメータ推定手段１４は、変換係数列候補を用いて算出される仮の積算発電量の算出結果から変換係数列を決定する場合も同様に、複数の評価期間についての変換係数列の決定結果を機械学習して最終的な変換パラメータを推定してもよい。
【０１００】
以上のように、本実施形態によれば、第１の時間単位の発電量が取得できない場合であっても、第１の時間単位の日射量の系列と実績発電量とが取得されれば、第１の時間単位の日射量から第１の時間単位の発電量に変換するための変換パラメータを得ることができる。
【０１０１】
比較例として、上記の特許文献２に記載の方法では、以下のように発電係数を算出した上で、日射量と発電係数とをもとに発電予測を行っている。なお、該文献の実施形態によれば、例えば１０分ごとの発電係数を算出している。
発電係数（ｔ）＝発電量ピーク値（ｔ）÷日射量強度ピーク値（ｔ）
【０１０２】
しかし、そのような発電係数の求め方は、発電量と日射量とが同じ時間単位で取得する必要があり、発電量と日射量とが異なる時間単位で取得されている場合には用いることができない。一般的に、日射量は気象に関する情報であることから短周期で取得することが可能であるが、発電量に関する情報は、発電システムごとにセンサを設けるなどの特別の手間をかけなければ、月に一度の電力量計にて計測される情報しか取得することができない。
【０１０３】
上記のような理由から、相対的に短い周期で取得できる日射量と、相対的に長い周期でしか取得できない発電量とに基づいて変換係数を求めようとすると、特許文献２に記載されたような単純に発電量と日射量で割るだけでは求めることができないため、本実施形態のように、テーブルなどにより係数列を準備する必要がある。
【０１０４】
実施形態２．
次に、本発明の第２の実施形態について説明する。図１０は、第２の実施形態の推定システムの例を示すブロック図である。図１０に示す推定システム２は、図１に示す第１の推定システムの構成に加えて、さらに発電量推定手段１５を備えている。
【０１０５】
以下、上記の実施形態と同じ手段については同じ符号を付し、説明を省略する。
【０１０６】
発電量推定手段１５は、変換パラメータ推定手段１４が推定した変換パラメータと、日射量列取得手段１１が取得した日射量列とに基づいて、該日射量列を取得した期間（評価期間）中の第１の時間単位の発電量を推定する。発電量推定手段１５は、変換パラメータによって特定される変換係数列を用いて、評価期間中の第１の時間単位の日射量の各々を発電量に変換することにより、評価期間中の第１の時間単位の発電量を推定する。
【０１０７】
発電量推定手段１５は、例えば、上記の式（２）または式（４）を用いて、評価期間中の第１の時間単位の日射量の各々に対して、変換パラメータによって特定される変換係数列のうち当該日射量に対応する時間ｔの変換係数を掛け合わせた値を、評価期間中の該第１の時間単位の発電量としてもよい。
【０１０８】
図１１は、本実施形態の推定システムの動作の一例を示すフローチャートである。なお、第１の実施形態の動作と同じ部分については同じ符号を付し、説明を省略する。
【０１０９】
図１１に示す例では、変換パラメータが決定されると、発電量推定手段１５が、決定された変換パラメータと、取得された日射量列とに基づいて、該日射量列を取得した評価期間中の第１の時間単位の発電量を推定する（ステップＳ２０１）。
【０１１０】
以上のように、本実施形態によれば、日射量列および実績発電量を得た過去の期間における第１の時間単位の発電量を得ることができる。
【０１１１】
実施形態３．
次に、本発明の第３の実施形態について説明する。図１２は、第３の実施形態の推定システムの例を示すブロック図である。図１３に示す推定システム３は、図１０に示す第２の推定システムの構成に加えて、さらに予測日射量取得手段１６を備えている。
【０１１２】
予測日射量取得手段１６は、ある対象発電システムについて、第１の時間単位の発電量を予測したい将来の期間（以下、推定期間という）が指定されると、対象発電システムの設置位置における該推定期間中の第１の時間単位の日射量を取得する。予測日射量取得手段１６は、例えば、指定された推定期間中の、対象発電システムの設置位置を含む所定の地域の気象に関する予報値を含むＧＰＶデータから推定期間中の第１の時間単位の日射量を推定してもよい。また、予測日射量取得手段１６は、他のシステムがＧＰＶデータ等から推定した、対象発電システムの設置位置を含む地域の推定期間中の第１の時間単位の日射量を取得してもよい。なお、日射量の取得方法は特に問わない。
【０１１３】
また、本実施形態の発電量推定手段１５は、変換パラメータ推定手段１４が推定した変換パラメータと、予測日射量取得手段１６が取得した推定期間中の第１の時間単位の日射量とに基づいて、推定期間中の第１の時間単位の発電量を推定する。
【０１１４】
発電量推定手段１５は、例えば、上記の式（２）または式（４）を用いて、推定期間中の第１の時間単位の日射量の各々に対して、変換パラメータによって特定される変換係数列のうち当該日射量に対応する時間ｔの変換係数を掛け合わせた値を、推定期間中の該第１の時間単位の発電量としてもよい。
【０１１５】
本実施形態で用いる変換パラメータは、対象発電システムについて既に決定ずみの変換パラメータがあればそれをそのまま用いてもよいし、推定期間が指定されたタイミングで過去の任意の評価期間中の日射量列と実績発電量とから新たに決定されたものであってもよい。なお、新たに決定する場合においても、評価期間は特に限定されないが、直近の過去の第２の時間単位の期間を評価期間としてもよいし、推定期間と最も気象情報（季節や気温や太陽の位置など）が近い過去の第２の時間単位の期間やを評価期間としてもよい。
【０１１６】
図１３は、本実施形態の推定システムの動作の一例を示すフローチャートである。なお、上記の実施形態の動作と同じ部分については同じ符号を付し、説明を省略する。
【０１１７】
図１３に示す例では、既に対象発電システムについて変換パラメータが決定されているものとする。推定期間が指定されると、予測日射量取得手段１６は、対象発電システムの設置位置における該推定期間中の第１の時間単位の日射量を取得する（ステップＳ３０１）。
【０１１８】
次に、発電量推定手段１５が、対象発電システムの変換パラメータと、取得された該推定期間中の第１の時間単位の日射量とに基づいて、推定期間中の第１の時間単位の発電量を推定する（ステップＳ３０２）。
【０１１９】
以上のように、本実施形態によれば、日射量の予測値が得られる期間であれば、将来のの任意の期間における第１の時間単位の発電量を得ることができる。
【０１２０】
実施形態４．
次に、本発明の第４の実施形態について説明する。図１４は、第４の実施形態の推定システムの例を示すブロック図である。図１４に示す推定システム４は、対象発電システムの位置情報と、実績発電量と、発電量を求めたい期間と時間単位に関する期間情報とを入力すると、期間情報で指定した期間中の指定した時間単位での発電量の系列である発電量列を出力するシステムである。
【０１２１】
本例において、実績発電量が算出された期間は、上記の評価期間とされる。なお、実績発電量とともに、その実績発電量が算出された期間（評価期間）の情報が入力されてもよい。また、期間情報で指定される期間および時間単位は、上記の推定期間および第１の時間単位とされる。また、位置情報で指定される位置が、対象発電システムの位置とされる。なお、既に説明したように、第１の時間単位は一定でなくてもよい。推定期間には、任意の期間が指定可能である。推定期間は、例えば、過去の期間であってもよいし、将来の期間であってもよいし、過去と将来の両方を含んだ期間であってもよい。一方、評価期間は、日射量の周期性（太陽赤緯と同様の周期性）を考慮して、該周期性の少なくとも一周期分に相当する長さが含まれていることが好ましい。また、位置情報は、例えば、緯度および経度の情報であってもよい。
【０１２２】
本実施形態の推定システムは、第２の実施形態の推定システム２の機能と第３の実施形態の推定システム３の機能とを有する。図１４に示す推定システム４は、変換情報記憶手段１３と、変換パラメータ推定手段１４と、発電量推定手段１５と、日射量推定手段１７と、入力手段１８と、出力手段１９とを含む。
【０１２３】
日射量推定手段１７は、指定された期間中の指定された位置における、指定された時間単位の日射量を推定する。なお、日射量推定手段１７は、上記の日射量列取得手段１１および予測日射量取得手段１６の機能をまとめた手段に相当する。
【０１２４】
日射量推定手段１７は、例えば、過去の期間の日射量であれば、指定された期間中に指定された位置を含む地域の上空を所定の時間間隔で撮影した衛星画像を取得して、該衛星画像から得られる雲量等に基づいて、指定された時間単位の日射量を推定してもよい。また、日射量推定手段１７は、例えば、指定された期間中の指定された位置または該位置近傍の地点における日射量を所定の時間間隔で計測した計測データを取得して、該計測データから指定された時間単位の日射量を推定してもよい。また、日射量推定手段１７は、例えば、指定された期間中の、指定された位置を含む地域上空の気象に関する予報値を含むＧＰＶデータを取得して、該ＧＰＶデータから得られる雲量等に基づいて、指定された時間単位の日射量を推定してもよい。また、日射量推定手段１７は、例えば、将来の期間の日射量であれば、指定された期間中の、指定された位置を含む地域上空の気象に関する予報値を含むＧＰＶデータを取得して、該ＧＰＶデータから得られる雲量等に基づいて、指定された時間単位の日射量を推定してもよい。
【０１２５】
なお、衛星画像や計測データやＧＰＶデータを得る時間間隔は、必ずしも第１の時間単位と同じでなくてもよい。そのような場合には、例えば、それらのデータから得られた日射強度等の情報を、日射量の時間変化を模式化した日射カーブにフィッティングさせて、所望の時間単位の日射量を得ればよい。また、日射量推定手段１７が推定する日射量は、全天日射量であってもよいし、直達日射量と散乱日射量とに分かれていてもよい。
【０１２６】
また、入力手段１８は、上記の実施形態の実績発電量取得手段１２に代わって、対象発電システムの実績発電量を取得する手段であるとともに、発電量の推定条件とされる情報を取得する手段である。入力手段１８は、例えば、対象発電システムの位置情報と、実績発電量と、期間情報とを入力する。入力手段１８は、これらをユーザに入力させるための入力インタフェースを有していてもよい。なお、実績発電量は２以上の評価期間中の実績発電量であってもよい。また、後述するように、設定済みの変換パラメータを用いる場合には、実績発電量の入力は省略可能である。
【０１２７】
また、入力手段１８は、発電量の推定条件とされる情報（位置情報と期間情報）とが入力されると、評価期間中の対象発電システムの設置位置における第１の時間単位の日射量を、日射量推定手段１７に推定させる。ここでは、評価期間が複数の第１の時間単位を含むことから、複数の日射量すなわち評価期間中の日射量列が推定される。また、入力手段１８は、推定期間中の対象発電システムの位置における第１の時間単位の日射量を、日射量推定手段１７に推定させる。
【０１２８】
また、出力手段１９は、発電量推定手段１５によって推定された推定期間中の第１の時間単位の発電量を出力する。
【０１２９】
本実施形態の変換パラメータ推定手段１４は、入力手段１８より実績発電量が入力されると、入力された実績発電量と、日射量推定手段１７によって推定された評価期間中の第１の時間単位の日射量の系列である日射量列と、変換情報記憶手段１３に記憶されている変換情報とに基づいて、変換パラメータを決定する。
【０１３０】
また、発電量推定手段１５は、日射量推定手段１７によって推定された推定期間中の第１の時間単位の日射量と、変換パラメータ推定手段１４が推定した変換パラメータとに基づいて、推定期間中の第１の時間単位の発電量を推定する。
【０１３１】
図１５は、本実施形態の推定システムの動作の一例を示すフローチャートである。図１５に示す例では、まず、入力手段１８が、位置情報、実績発電量および期間情報を入力する（ステップＳ４０１）。既に説明したように、設定済みの変換パラメータを用いる場合には、実績発電量の入力は省略可能である。
【０１３２】
次に、日射量推定手段１７が、評価期間中の対象発電システムの位置における日射量列（第１の時間単位の日射量の系列）を推定する（ステップＳ４０２）。なお、実績発電量が複数入力された場合には、日射量推定手段１７は、各実績発電量が算出された評価期間中の対象発電システムの位置における日射量列を推定する。
【０１３３】
次に、変換パラメータ推定手段１４は、入力手段１８から出力される実績発電量と、推定された該実績発電量の評価期間中の日射量列と、変換情報記憶手段１３に記憶されている変換情報とに基づいて、変換パラメータを決定する（ステップＳ４０３）。なお、実績発電量と日射量列の組が複数入力された場合には、変換パラメータ推定手段１４は、各組につき変換パラメータを決定した上で、最終的な変換パラメータを決定する。ステップＳ４０３の動作は、第１〜第３の実施形態のステップＳ１０２の動作と同様でよい。また、設定済みの変換パラメータを用いる場合には、ステップＳ４０３の動作は省略される。
【０１３４】
次に、日射量推定手段１７が、推定期間中の対象発電システムの位置における第１の時間単位の日射量を推定する（ステップＳ４０４）。ここで、推定期間が複数の第１の時間単位を含む場合、日射量推定手段１７は、複数の第１の時間単位の日射量すなわち推定期間中の日射量列を推定する。
【０１３５】
次に、発電量推定手段１５が、推定された推定期間中の第１の時間単位の日射量と、変換パラメータ推定手段１４が決定した変換パラメータとに基づいて、推定期間中の第１の時間単位の発電量を推定する（ステップＳ４０５）。ここで、推定期間が複数の第１の時間単位を含む場合、発電量推定手段１５は、複数の第１の時間単位の発電量すなわち推定期間中の発電量列を推定する。なお、発電量推定手段１５は、推定期間中の第１の時間単位の日射量の各々に対して、変換パラメータによって特定される変換係数列のうち該日射量の時間ｔに対応した変換係数を適用して、対応する時間ｔの発電量を得ればよい。
【０１３６】
最後に、出力手段１９が、推定結果である推定期間中の第１の時間単位の発電量を出力する（ステップＳ４０６）。
【０１３７】
以上のように、本実施形態によれば、対象発電システムの任意の期間の任意の時間単位の発電量を、その期間および時間単位の情報と、対象発電システムの位置情報と、実績発電量とを指定するだけで得ることができる。さらに、変換パラメータの設定後は、発電量の推定条件の指定を行うだけで、対象発電システムの任意の期間の任意の時間単位の発電量を得ることができる。
【０１３８】
なお、上記の各実施形態では、変換パラメータ推定手段１４が、第１の時間単位の日射量と実績発電量とが入力されたタイミングで変換パラメータを推定し、その結果またはその結果に基づく学習結果等を、最終的な変換パラメータに設定する例を示したが、変換パラメータ推定手段１４は、変換パラメータの設定を、例えば、以下のようなタイミングで行ってもよい。
【０１３９】
例えば、変換パラメータ推定手段１４は、一度設定した後は、同じ変換パラメータを利用しつづけてもよい。また、例えば、変換パラメータ推定手段１４は、第１の時間単位の発電量の推定や予測を行う度に、対応する評価期間中の第１の時間単位の日射量と実績発電量とから変換パラメータを推定し、再設定してもよい。そのような場合、変換パラメータ推定手段１４は、例えば、新たに、実績発電量と該実績発電量を算出した期間中の第１の時間単位の日射量とが入力される度に、変換パラメータを推定し、再設定してもよい。
【０１４０】
また、例えば、変換パラメータ推定手段１４は、季節の変化などにより気温や太陽の位置など気象情報が変化する所定の期間（例えば、１週間、１ヶ月、３カ月等）ごとに、変換パラメータを推定しなおしてもよい。例えば、システムが、周期的に直近の評価期間中の第１の時間単位の日射量と実績発電量とを取得する場合に、変換パラメータ推定手段１４は、そのうちの気象情報が変化する所定の期間ごとに、現在設定されている変換パラメータをクリアして、新たに取得された情報から変換パラメータを推定しなおしてもよい。
【０１４１】
また、例えば、変換パラメータ推定手段１４は、現在設定されている変換パラメータにより推定したある推定期間中の発電量（積算発電量）とその後得られた該推定期間中の実際の発電量の実績値（実績発電量）との差から算出される予測精度が一定の閾値を下回った場合に、変換パラメータを再設定するようにしてもよい。例えば、変換パラメータ推定手段１４は、推定期間が明けたタイミングで予測精度を評価し、その結果予測精度が一定の閾値を下回っていた場合に、現在設定されている変換パラメータをクリアして、その後に実績発電量と日射量列が入力されたタイミング（次回の予測依頼時や、一定周期後等）で変換パラメータを推定しなおしてもよい。なお、変換パラメータ推定手段１４は、変換パラメータの再設定を決定した場合に、即座に実績発電量と日射量列とを要求して、推定しなおすことも可能である。予測精度に応じて変換パラメータを再設定する方法によれば、計算コストを抑えつつ、例えば新しいビルの建築などの周囲状況の変化によって日照量に変化が生じた場合やパネルの経年的な劣化が生じた場合であっても、精度の高い変換パラメータを得ることができる。
【０１４２】
なお、上記の変換パラメータの設定・クリア・再設定は、そのまま変換情報に含まれる変換係数列の候補や発電システムに関するパラメータや制約条件の決定・決定状態のクリア・再決定と読み替えてもよい。
【０１４３】
また、上記の各実施形態では、発電システムごとに、変換パラメータや第１の時間単位の発電量を推定する例を示したが、変換パラメータおよび第１の時間単位の発電量の推定は、例えば、所定の地域単位で行うことも可能である。
【０１４４】
そのような場合、例えば、対象地域で得られる評価期間中の実績発電量と、評価期間中の対象地域における平均日射量の系列とを基に、変換パラメータを推定すればよい。変換パラメータの推定方法は、発電システム単位の場合と同様である。地域単位で変換パラメータを推定すれば、該地域内の発電システムの数やそれらの仕様が明らかでない場合にも、その地域の日射量の情報と実績発電量とから、実績発電量を得た時間単位より短い時間単位の発電量の情報を容易に得ることができる。なお、変換パラメータの推定対象とする単位が異なる場合においても、変換パラメータの設定タイミングの例として上記の例が挙げられる。
【０１４５】
次に、本発明の各実施形態にかかるコンピュータの構成例を示す。図１６は、本発明の各実施形態にかかるコンピュータの構成例を示す概略ブロック図である。コンピュータ１０００は、ＣＰＵ１００１と、主記憶装置１００２と、補助記憶装置１００３と、インタフェース１００４と、ディスプレイ装置１００５とを備える。
【０１４６】
上述の各実施形態の推定システムは、コンピュータ１０００に実装されてもよい。その場合、推定システムの動作は、プログラムの形式で補助記憶装置１００３に記憶されていてもよい。ＣＰＵ１００１は、プログラムを補助記憶装置１００３から読み出して主記憶装置１００２に展開し、そのプログラムに従って各実施形態における所定の処理を実施する。
【０１４７】
補助記憶装置１００３は、一時的でない有形の媒体の一例である。一時的でない有形の媒体の他の例として、インタフェース１００４を介して接続される磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等が挙げられる。また、このプログラムが通信回線によってコンピュータ１０００に配信される場合、配信を受けたコンピュータは１０００がそのプログラムを主記憶装置１００２に展開し、各実施形態における所定の処理を実行してもよい。
【０１４８】
また、プログラムは、各実施形態における所定の処理の一部を実現するためのものであってもよい。さらに、プログラムは、補助記憶装置１００３に既に記憶されている他のプログラムとの組み合わせで各実施形態における所定の処理を実現する差分プログラムであってもよい。
【０１４９】
また、実施形態における処理内容によっては、コンピュータ１０００の一部の要素は省略可能である。例えば、推定システムがユーザに情報を提示しないのであれば、ディスプレイ装置１００５は省略可能である。また、図１６には図示省略しているが、実施形態における処理内容によっては、コンピュータ１０００は、入力デバイスを備えていてもよい。例えば、第４の実施形態の推定システム４の場合に、位置情報や実績発電量や期間情報を入力するための入力デバイスを備えていてもよい。
【０１５０】
また、各装置の各構成要素の一部または全部は、汎用または専用の回路（Circuitry）、プロセッサ等やこれらの組み合わせによって実施される。これらは単一のチップによって構成されてもよいし、バスを介して接続される複数のチップによって構成されてもよい。また、各装置の各構成要素の一部又は全部は、上述した回路等とプログラムとの組み合わせによって実現されてもよい。
【０１５１】
各装置の各構成要素の一部又は全部が複数の情報処理装置や回路等により実現される場合には、複数の情報処理装置や回路等は、集中配置されてもよいし、分散配置されてもよい。例えば、情報処理装置や回路等は、クライアントアンドサーバシステム、クラウドコンピューティングシステム等、各々が通信ネットワークを介して接続される形態として実現されてもよい。
【０１５２】
次に、本発明の概要について説明する。図１７は、本発明の概要を示すブロック図である。図１７に示すように、本発明の推定システムは、変換パラメータ推定手段１０１を備えている。
【０１５３】
変換パラメータ推定手段１０１は、対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とに基づいて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する。
【０１５４】
このように、比較的長期間の発電量の実績値と、該実績値を得た期間中の比較的短期間の日射量の系列とを基に、変換係数列を推定することにより、対象発電システムの特性や設置条件が不明であっても比較的短時間の日射量から発電量に変換できるので、より詳細な情報を得ることができる。
【０１５５】
なお、本発明の推定システムは、例えば、図１８や図１９に示すような構成であってもよい。図１８および図１９は、本発明による推定システムの他の構成を示すブロック図である。
【０１５６】
すなわち、推定システムは、予め、変換情報として、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報を記憶する変換情報記憶手段１０２を備え、変換パラメータ推定手段１０１は、変換情報によって特定される変換係数列の候補または制約条件を用いて、変換パラメータを推定してもよい。
【０１５７】
また、変換情報記憶手段１０２は、変換情報として、複数の変換係数列の候補の情報を記憶し、変換パラメータ推定手段１０３は、変換情報によって示される変換係数列の候補を用いて仮の積算発電量を算出した結果に基づいて変換係数列を決定し、決定された変換係数列または該変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定する
【０１５８】
また、第１の時間単位の日射量は、直達日射量と散乱日射量とに分かれており、変換情報は、直達日射量に適用する第１の変換係数列の候補の情報と、散乱日射量に適用する第２の変換係数列の候補の情報とを含み、変換パラメータ推定手段１０１は、変換情報によって示される第１の変換係数列の候補と第２の変換係数列の候補とを少なくとも用いて仮の積算発電量を算出した結果に基づいて、１組の変換係数列を決定し、決定された変換係数列の組または該組の各変換係数列の生成に用いた発電システムに関するパラメータを、変換パラメータとして推定してもよい。
【０１５９】
また、変換情報は、複数の方位角および複数の傾斜角の組み合わせにおける水平面日射量と傾斜面日射量との比率が少なくとも表現された変換係数列の候補の情報を含んでいてもよい。
【０１６０】
また、第２の時間単位は、太陽赤緯の周期性における一周期以上であり、変換情報は、太陽赤緯と同じの周期性をもった変換係数列の候補の情報を含んでいてもよい。
【０１６１】
また、変換情報は、一定値の変換係数列の候補の情報を含んでいてもよい。
【０１６２】
また、変換情報は、変換係数列の各候補の生成に用いた発電システムに関するパラメータの情報を含んでいてもよい。
【０１６３】
また、変換パラメータ推定手段１０１は、変換情報によって示される変換係数列の候補の中から選択した２以上の変換係数列の候補を掛け合わせて新たな変換係数列の候補を生成し、生成された新たな変換候補列の候補を含む変換係数列の候補を用いて仮の積算発電量を算出してもよい。
【０１６４】
また、変換情報は、変換係数列の候補の生成に用いた発電システムに関するパラメータの種類別に変換係数列の候補が登録される複数の変換係数列候補テーブルを含み、変換パラメータ推定手段１０１は、変換係数列候補テーブルの各々から１または複数の変換係数列の候補を選択し、それらを掛け合わせて新たな変換係数列の候補を生成してもよい。
【０１６５】
また、変換情報記憶手段は、変換情報として、所定の数理計画モデルに与える制約条件の情報を記憶し、変換パラメータ推定手段１０１は、変換情報によって示される制約条件を用いて所定の数理計画モデルを解くことにより変換係数列を決定し、決定された変換係数列または該変換係数列を決定した際の制約条件が示す発電システムに関するパラメータを、変換パラメータとして推定してもよい。
【０１６６】
また、推定システムは、変換パラメータと、対象発電システムの設置位置または対象地域における任意の期間中の第１の時間単位の日射量とに基づいて、任意の期間中の第１の時間単位の発電量を推定する発電量推定手段１０３を備えていてもよい。
【０１６７】
なお、上記の各実施形態は以下の付記のようにも記載できる。
【０１６８】
（付記１）位置情報と、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、第１の時間単位の発電量を求める推定期間に関する期間情報とを少なくとも入力する入力手段と、位置情報で示される位置における実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列、および位置情報で示される位置における推定期間中の第１の時間単位の日射量を推定する日射量推定手段と、実績発電量と、実績発電量を算出した期間中の日射量列とに基づいて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する変換パラメータ推定手段と、変換パラメータと、推定期間中の第１の時間単位の日射量とに基づいて、推定期間中の第１の時間単位の発電量を推定する発電量推定手段と、推定結果を出力する出力手段とを備えたことを特徴とする推定システム。
【０１６９】
（付記２）対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とが入力されると、予め記憶されている変換情報であって、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報である変換情報によって示される変換係数列の候補または制約条件を用いて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定することを特徴とする推定方法。
【０１７０】
（付記３）位置情報と、実績発電量と、第１の時間単位の発電量を求める推定期間に関する期間情報とが少なくとも入力されると、位置情報で示される位置における実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列および位置情報で示される位置における推定期間中の第１の時間単位の日射量を取得し、実績発電量と、実績発電量を算出した期間中の日射量列とに基づいて、変換パラメータを推定し、変換パラメータと、推定期間中の第１の時間単位の日射量とに基づいて、推定期間中の第１の時間単位の発電量を推定する付記２に記載の推定方法。
【０１７１】
（付記４）コンピュータに、対象発電システムまたは対象地域の、第１の時間単位よりも長い第２の時間単位の発電量の実績値である実績発電量と、対象発電システムの設置位置または対象地域における、実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列とが入力されると、予め記憶されている変換情報であって、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換する変換式に用いられる変換係数の系列である変換係数列の候補の情報、または変換後の第１の時間単位の発電量の積算値である積算発電量が実績発電量に最も近くなる変換係数列を解く所定の数理計画モデルに与える制約条件の情報である変換情報によって示される変換係数列の候補または制約条件を用いて、第２の時間単位中の第１の時間単位の日射量の各々を発電量に変換するための変換パラメータを推定する処理を実行させるための推定プログラム。
【０１７２】
（付記５）コンピュータに、位置情報と、実績発電量と、第１の時間単位の発電量を求める推定期間に関する期間情報とを少なくとも入力する処理、位置情報で示される位置における実績発電量を算出した期間中の第１の時間単位の日射量の系列である日射量列、および位置情報で示される位置における推定期間中の第１の時間単位の日射量を取得する処理、および推定された変換パラメータと、推定期間中の第１の時間単位の日射量とに基づいて、推定期間中の第１の時間単位の発電量を推定する処理を実行させる付記４に記載の推定プログラム。
【０１７３】
以上、本実施形態および実施例を参照して本願発明を説明したが、本願発明は上記実施形態および実施例に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
【産業上の利用可能性】
【０１７４】
本発明は、対象発電システムの比較的短時間の発電量だけでなく、対象発電システムのパラメータや、１つ以上の発電システムを含む所定の地域の比較的短時間の発電量などを得る用途に好適に適用可能である。
【符号の説明】
【０１７５】
１、２、３、４  推定システム
１１  日射量列取得手段
１２  実績発電量取得手段
１３  変換情報記憶手段
１４  変換パラメータ推定手段
１５  発電量推定手段
１６  予測日射量取得手段
１７  日射量推定手段
１８  入力手段
１９  出力手段
１０１  変換パラメータ推定手段
１０２  変換情報記憶手段
１０３  発電量推定手段
１０００  コンピュータ
１００１  ＣＰＵ
１００２  主記憶装置
１００３  補助記憶装置
１００４  インタフェース
１００５  ディスプレイ装置
【図１】
2018007461.tif 000006
【図２】
2018007461.tif 000007
【図３】
2018007461.tif 000008
【図４】
2018007461.tif 000009
【図５】
2018007461.tif 000010
【図６】
2018007461.tif 000011
【図７】
2018007461.tif 000012
【図８】
2018007461.tif 000013
【図９】
2018007461.tif 000014
【図１０】
2018007461.tif 000015
【図１１】
2018007461.tif 000016
【図１２】
2018007461.tif 000017
【図１３】
2018007461.tif 000018
【図１４】
2018007461.tif 000019
【図１５】
2018007461.tif 000020
【図１６】
2018007461.tif 000021
【図１７】
2018007461.tif 000022
【図１８】
2018007461.tif 000023
【図１９】
2018007461.tif 000024