(57)【要約】  本発明は、機械学習による統計的手法を用いて自然エネルギー発電量、日射量又は風速を予測する技術において、予測の精度を高めることを課題とする。当該課題を解決するため、対象時刻よりもｍ時間前（ｍは２以上）から対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出部（１３）と、複数日分の特徴量に基づいて対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定部（第１の推定部（１４））と、を有する予測装置（１０）を提供する。
2016121202.tif 000002
【特許請求の範囲】
【請求項１】
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段と、
複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段と、を有する予測装置。
【請求項２】
前記推定手段は、
前記対象時刻までの気象データから抽出した特徴量を説明変数、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数として予測する予測式を用いて推定を行う請求項１に記載の予測装置。
【請求項３】
前記推定手段は、前記説明変数および前記目的変数の組み合わせからなる複数日分の教師データに基づいた予測式を用いて推定を行う請求項２に記載の予測装置。
【請求項４】
前記推定手段は、予測対象日と所定の属性が類似する複数日分の前記特徴量に基づいて推定を行う請求項１から３のいずれか１項に記載の予測装置。
【請求項５】
前記ｍの値は可変であり、前記ｍの値を設定するｍ値設定手段をさらに有する請求項１から４のいずれか１項に記載の予測装置。
【請求項６】
前記推定手段は、複数の地域の自然エネルギー発電量、日射量又は風速を推定し、
前記ｍ値設定手段は、地域毎に前記ｍの値を設定する請求項５に記載の予測装置。
【請求項７】
前記ｍ値設定手段は、予測対象日の属性に基づいて、前記ｍの値を設定する請求項５又は６に記載の予測装置。
【請求項８】
前記ｍ値設定手段は、ｍの値各々における推定の精度を出力する手段を備える請求項５から７のいずれか１項に記載の予測装置。
【請求項９】
前記推定手段は、予測対象日及び予測対象地点の少なくとも一方が特定された予測対象と、所定の属性が所定レベル以上類似する教師データに基づき生成された前記予測式を用いて推定を行う請求項１から８のいずれか１項に記載の予測装置。
【請求項１０】
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段と、
前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段と、を有する予測装置。
【請求項１１】
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段と、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段と、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段と、
前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段と、
を有する予測装置。
【請求項１２】
コンピュータが、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出工程と、
複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定工程と、
を実行する予測方法。
【請求項１３】
コンピュータを、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段、
複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、予測装置、予測方法、及び、プログラムに関し、より具体的には、自然エネルギー発電量、日射量及び／又は風速を予測する予測装置、予測方法、及び、プログラムに関する。
【背景技術】
【０００２】
特許文献１乃至３、及び、非特許文献１に、機械学習による統計的手法を用いて、気象データから太陽光発電の発電量や日射量を予測する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平９−２１５１９２号
【特許文献２】特許第３９８４６０４号
【特許文献３】特許第５３３９３１７号
【非特許文献】
【０００４】
【非特許文献１】Joao Gari da Silva Fonseca Junior, Takashi Oozeki, Takumi Takashima, and Kazuhiko Ogimoto, Analysis of the Use of Support Vector Regression and Neural Networks to Forecast Insolation for 25 Locations in Japan, Solar World Congress 2011 Proceedings, Germany, International Solar Energy Society, 2011, pp.4128-4135.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１乃至３、及び、非特許文献１に開示の技術の場合、予測の精度が十分でなかった。本発明は、機械学習による統計的手法を用いて自然エネルギー発電量、日射量及び／又は風速を予測する技術において、予測の精度を高めることを課題とする。
【課題を解決するための手段】
【０００６】
本発明によれば、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段と、
複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段と、を有する予測装置が提供される。【０００７】
また、本発明によれば、
コンピュータが、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出工程と、
複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定工程と、
を実行する予測方法が提供される。【０００８】
また、本発明によれば、
コンピュータを、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段、
複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段、
として機能させるためのプログラムが提供される。【０００９】
また、本発明によれば、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段と、
前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段と、を有する予測装置が提供される。【００１０】
また、本発明によれば、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段と、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段と、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段と、
前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段と、
を有する予測装置が提供される。
【発明の効果】
【００１１】
本発明によれば、機械学習による統計的手法を用いて太陽光発電の発電量や日射量を予測する技術において、予測の精度を高めることができる。
【図面の簡単な説明】
【００１２】
上述した目的、およびその他の目的、特徴および利点は、以下に述べる好適な実施の形態、およびそれに付随する以下の図面によってさらに明らかになる。【００１３】
【図１】本実施形態の装置のハードウエア構成の一例を概念的に示す図である。
【図２】本実施形態の予測装置の機能ブロック図の一例である。
【図３】本実施形態の予測式取得部の機能ブロック図の一例である。
【図４】本実施形態の予測装置が利用する過去データの一例を模式的に示す図である。
【図５】本実施形態の概要を説明するための図である。
【図６】本実施形態の予測装置の機能ブロック図の一例である。
【図７】本実施形態の予測装置が表示する情報の一例を示す図である。
【図８】本実施形態の予測装置が表示する情報の一例を示す図である。
【図９】本実施形態の予測装置が表示する情報の一例を示す図である。
【図１０】本実施形態の予測装置の機能ブロック図の一例である。
【図１１】本実施形態の予測装置の検証結果を示す図である。
【図１２】本実施形態の予測装置の機能ブロック図の一例である。
【発明を実施するための形態】
【００１４】
まず、本実施形態の装置のハードウエア構成の一例について説明する。本実施形態の装置が備える各部は、任意のコンピュータのＣＰＵ（Central Processing Unit）、メモリ、メモリにロードされるプログラム、そのプログラムを格納するハードディスク等の記憶ユニット（あらかじめ装置を出荷する段階から格納されているプログラムのほか、ＣＤ（Compact Disc）等の記憶媒体やインターネット上のサーバ等からダウンロードされたプログラムをも格納できる）、ネットワーク接続用インタフェイスを中心にハードウエアとソフトウエアの任意の組合せによって実現される。そして、その実現方法、装置にはいろいろな変形例があることは、当業者には理解されるところである。
【００１５】
図１は、本実施形態の装置のハードウエア構成の一例を概念的に示す図である。図示するように、本実施形態の装置は、例えば、バス１０Ａで相互に接続されるＣＰＵ１Ａ、ＲＡＭ（Random Access Memory）２Ａ、ＲＯＭ（Read Only Memory）３Ａ、表示制御部４Ａ、ディスプレイ５Ａ、操作受付部６Ａ、操作部７Ａ、通信部８Ａ、補助記憶装置９Ａ等を有する。なお、図示しないが、その他、外部機器と有線で接続される入出力インタフェイス、マイク、スピーカ等の他の要素を備えてもよい。
【００１６】
ＣＰＵ１Ａは各要素とともに装置のコンピュータ全体を制御する。ＲＯＭ３Ａは、コンピュータを動作させるためのプログラムや各種アプリケーションプログラム、それらのプログラムが動作する際に使用する各種設定データなどを記憶する領域を含む。ＲＡＭ２Ａは、プログラムが動作するための作業領域など一時的にデータを記憶する領域を含む。補助記憶装置９Ａは、例えばＨＤＤ（Hard Disc Drive）であり、大容量のデータを記憶可能である。
【００１７】
ディスプレイ５Ａは、例えば、表示装置（ＬＥＤ（Light Emitting Diode）表示器、液晶ディスプレイ、有機ＥＬ（Electro Luminescence）ディスプレイ等）である。ディスプレイ５Ａは、タッチパッドと一体になったタッチパネルディスプレイであってもよい。表示制御部４Ａは、ＶＲＡＭ（Video RAM）に記憶されたデータを読み出し、読み出したデータに対して所定の処理を施した後、ディスプレイ５Ａに送って各種画面表示を行う。操作受付部６Ａは、操作部７Ａを介して各種操作を受付ける。操作部７Ａは、操作キー、操作ボタン、スイッチ、ジョグダイヤル、タッチパネルディスプレイ、キーボードなどを含む。通信部８Ａは、有線及び／または無線で、インターネット、ＬＡＮ（Local Area Network）等のネットワークに接続し、他の電子機器と通信する。
【００１８】
以下、本実施の形態について説明する。なお、以下の実施形態の説明において利用する機能ブロック図は、ハードウエア単位の構成ではなく、機能単位のブロックを示している。これらの図においては、各装置は１つの機器により実現されるよう記載されているが、その実現手段はこれに限定されない。すなわち、物理的に分かれた構成であっても、論理的に分かれた構成であっても構わない。なお、同一の構成要素には同一の符号を付し、適宜説明を省略する。
【００１９】
＜第１の実施形態＞
本実施形態の予測装置１０は、対象時刻よりもｍ時間前（ｍは２以上）から当該対象時刻までの気象データから抽出した特徴量を説明変数とし、当該対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された予測式を利用して、予測対象日の当該対象時刻の自然エネルギー発電量、日射量又は風速を予測する。
【００２０】
自然エネルギー発電量とは、自然エネルギーを利用した発電の発電量を意味する。このような発電方式としては、太陽光を利用した発電、風力を利用した発電等が考えられる。以下、本実施形態を詳細に説明する。
【００２１】
図１２に、本実施形態の予測装置１０の機能ブロック図の一例を示す。図示するように、予測装置１０は、特徴量抽出部１３と、第１の推定部１４とを有する。特徴量抽出部１３は、対象時刻よりもｍ時間前（ｍは２以上）から当該対象時刻までの気象データから時系列の変動である特徴量を抽出する。第１の推定部１４は、複数日分の特徴量に基づいて対象時刻における自然エネルギー発電量、日射量又は風速を推定する。なお、第１の推定部１４は、対象時刻までの気象データから抽出した特徴量を説明変数、対象時刻の自然エネルギー発電量、日射量又は風速を目的変数として予測する予測式を用いて推定を行ってもよい。また、第１の推定部１４は、説明変数および目的変数の組み合わせからなる複数日分の教師データに基づいた予測式を用いて推定を行ってもよい。
【００２２】
図２に、本実施形態の予測装置１０の機能ブロック図の他の一例を示す。図示するように、予測装置１０は、予測式取得部１１と、気象データ取得部１２と、特徴量抽出部１３と、第１の推定部１４とを有する。以下、各部について説明する。
【００２３】
予測式取得部１１は、対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する。当該予測式は、対象時刻よりもｍ時間前から当該対象時刻までの気象データから抽出した特徴量を説明変数とし、当該対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成される。予測式取得部１１は、このような予測式を生成してもよいし、予測装置１０と通信可能になっている他の外部装置から、有線及び／又は無線での通信により取得してもよい。
【００２４】
予測式を生成する実施形態における予測式取得部１１の機能ブロック図の一例を、図３に示す。図示する予測式取得部１１は、過去データ記憶部２１と予測式生成部２２とを有する。なお、予測式取得部１１が外部装置から予測式を取得する場合、当該外部装置が過去データ記憶部２１及び予測式生成部２２を備えることとなる。
【００２５】
過去データ記憶部２１は、過去の日時毎に（所定の日の所定の時刻毎に）、気象データの実績値又は予測値（各時刻より前の所定のタイミングで発表された予測値）と、自然エネルギー発電量、日射量及び／又は風速の実績値と、これらの属性を示す属性値と、を対応付けた過去データを記憶している。過去データ記憶部２１は、複数日分（例：３０日分、６０日分、１年分、３年分等）の過去データを記憶する。
【００２６】
図４に、過去データ記憶部２１が記憶する過去データの一例を模式的に示す。図示する過去データは、日付と、時刻と、太陽光発電量と、日射量と、気象データと、属性データとが互いに対応付けられている。図示しないが、さらに、風速の実績値及び／又は風力発電の発電量の実績値が対応付けられていてもよい。
【００２７】
過去データは、所定の時間間隔で蓄積された複数のデータを含む。データの時間間隔は様々であり、５分毎、１５分毎、３０分毎、１時間毎等、任意に選択できる。なお、過去データは、さらに観測地毎に蓄積されてもよい。すなわち、過去データは、観測地毎に、所定の時間間隔で蓄積されていてもよい。
【００２８】
太陽光発電量及び日射量の欄には、対応付けられている日時を基準にして特定される所定の時間内における積算量の実績値が記載される。例えば、対応付けられている日時を中心としたＭ分間の積算量（Ｍは、例えば５、１５、３０、６０等）や、対応付けられている日時からＭ分後までの積算量等が考えられるが、これらに限定されない。実績データが観測地毎に蓄積されている場合、太陽光発電量及び日射量の欄には、各観測地における日射量の実績値、及び、各観測地に設置された太陽光発電装置で発電された太陽光発電量の実績値が記載される。
【００２９】
図示しないが、風力発電の欄を有する場合、同様に、対応付けられている日時を基準にして特定される所定の時間内における積算量の実績値が記載される。風速の欄には、対応付けられている日時における実績値、又は、対応付けられている日時を基準にして特定される所定の時間内における実績値の統計値（平均値、最大値、最頻値、中央値、最小値等）が記載される。
【００３０】
気象データの欄には、対応付けられている日時における実績値が記載される。なお、過去データの時間間隔と気象データのサンプリング間隔とが相違する等の理由により、過去データの日時ちょうどに測定された気象データがない場合、当該日時に最も近いタイミングで測定された気象データを採用してもよい。その他、気象データの欄には、対応付けられている日時を基準にして特定される所定の時間内における実績値の統計値（平均値、最大値、最頻値、中央値、最小値等）が記載されてもよい。また、気象データの欄には、実績値に代えて、対応付けられている時刻よりも前の所定のタイミングで発表された予測値が記載されてもよい。当該予測値は、前日等に発表される気象予報の値等が該当する。
【００３１】
気象データは、自然エネルギー発電量や日射量や風速に影響する項目のデータを少なくとも１つ含む。例えば、気象データは、気温、湿度、風向き、風速、降水量、天気、上層雲量、中層雲量、下層雲量、全雲量、地上気圧、海面気圧、日射量等の項目が考えられるが、これらに限定されない。実績データが観測地毎に蓄積されている場合、気象データの欄には、各観測地の実績値又は予測値が記載される。
【００３２】
属性データの欄には、各データの属性を示す値が記載される。属性データは、自然エネルギー発電量や日射量や風速に影響する項目のデータを少なくとも１つ含む。例えば、属性データは、観測地や、観測日の季節等が考えられるが、これらに限定されない。観測地は、市町村名で示されてもよいし、緯度・経度で示されてもよいし、その他であってもよい。
【００３３】
図３に戻り、予測式生成部２２は、過去データ記憶部２１が記憶する過去データに基づいて、対象時刻における自然エネルギー発電量、日射量又は風速を予測する予測式を生成する。具体的には、予測式生成部２２は、対象時刻よりもｍ時間前（ｍは２以上）から当該対象時刻までの気象データから抽出した特徴量を説明変数とし、当該対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で、予測式を生成する。
【００３４】
上記特徴量は、対象時刻よりもｍ時間前から対象時刻までの時間内における気象データの時系列な変動の特徴を表す。特徴量を抽出するアルゴリズムは様々である。例えば、当該時間内における所定の１つ又は複数の項目（気象データ）の値そのものを時系列に並べた１次元配列又は多次元配列を、特徴量としてもよい。又は、一方の軸に所定の項目（気象データ）の値をとり、他方の軸に時間をとったグラフ上にデータをプロットし、得られた波形の形状から、上記変動を表す任意の特徴量を抽出してもよい。また、複数の項目（気象データ）から当該手法（波形の形状）で特徴量を抽出し、これらを所定の項目順に並べた配列を特徴量としてもよい。
【００３５】
機械学習の手法には、重回帰、ニューラルネットワーク、サポートベクターマシン等のあらゆる手法を採用できる。
【００３６】
ｍの値の下限は２であり、好ましくは５、さらに好ましくは９である。以下の実施例で示すが、このようにすることで、自然エネルギー発電量、日射量又は風速の予測の精度を十分に高めることができる。ｍの値の上限は、例えば２０、好ましくは１３である。以下の実施例で示すように、ｍの値が所定値以下の場合、ｍの値が大きくなるほど予測の精度は高くなる。しかし、ｍの値が所定値を超えると、予測の精度はほぼ横ばいとなり、大きな変化は得られなくなる。ｍの上限をこのように定めることで、十分な予測の精度を実現しつつ、処理対象のデータ量を減らすことでコンピュータの処理負担を軽減できる。
【００３７】
予測式生成部２２は、互いに異なる複数の対象時刻各々に対応した複数の予測式を生成することができる。
【００３８】
気象データ取得部１２は、予測対象日の対象時刻までの気象データ（時系列データ）を取得する。気象データ取得部１２は、少なくとも、予測対象日の対象時刻よりもｍ時間前から当該対象時刻までの気象データを取得する。例えば、気象データ取得部１２は、有線及び／又は無線での通信により外部装置と通信し、当該気象データを取得することができる。気象データ取得部１２は、観測地毎に、当該気象データを取得してもよい。
【００３９】
気象データ取得部１２が取得する気象データは、実績値であってもよいし、予測値であってもよいし、これらが混在していてもよい。気象データ取得部１２が予測対象日の対象時刻よりもｍ時間前から当該対象時刻までの気象データを取得する時点で、その一部又は全部の実績値がいまだに公表されていない場合があり得る。全部の実績値が公表されていない場合、気象データ取得部１２は、予測対象日の対象時刻よりもｍ時間前から当該対象時刻までの気象データとして、予測値を取得する。一方、一部の実績値が公表され、一部の実績値が公表されていない場合、気象データ取得部１２は、公表されている実績値を取得するとともに、実績値が公表されていない時間帯においては予測値を取得してもよい。その他、一部の実績値が公表され、一部の実績値が公表されていない場合、気象データ取得部１２は、すべての時間帯において予測値を取得してもよい。
【００４０】
特徴量抽出部１３は、気象データ取得部１２が取得した気象データに基づいて、所定の処理を行う。具体的には、特徴量抽出部１３は、予測対象日の対象時刻よりもｍ時間前から対象時刻までの気象データから特徴量を抽出する。特徴量抽出部１３が抽出する特徴量は、予測式取得部１１が取得する予測式の生成において説明変数として用いられた特徴量と同種の特徴量である。
【００４１】
第１の推定部１４は、予測式取得部１１が取得した予測式と、特徴量抽出部１３が抽出した特徴量とに基づいて、予測対象日の対象時刻における自然エネルギー発電量、日射量又は風速を推定する。すなわち第１の推定部１４は、予測式取得部１１が取得した予測式に、特徴量抽出部１３が抽出した特徴量を入力することで、予測対象日の対象時刻における自然エネルギー発電量、日射量又は風速の推定値（出力）を得る。なお、第１の推定部１４は、日射量の推定値を得た場合、その後、日射量の推定値に変換係数を乗算するなどにより、太陽光発電量を算出してもよい。また、第１の推定部１４は、風速の推定値を得た場合、所定の式に入力し、風力発電量を算出してもよい。風力発電量は、ローター面積（予めユーザが指定）や、風速（推定値）の３乗に比例することが知られている。
【００４２】
ここで、図５に示す具体例を用いて、予測装置１０による処理の概念を説明する。例えば、予測対象日が２０１５年１月１日であり、対象時刻が１８時であり、ｍの値は１２であるとする。この場合、対象時刻よりもｍ時間前は６時となる。
【００４３】
図５には、気象データの一例として気温のデータを示している。当該例の場合、２０１５年１月１日（予測対象日）より前の任意の複数日（図の場合、ｐ日）分のデータを教師データとして、予測式が生成されている。具体的には、各日の６時から１８時までの気象データから抽出された特徴量が、説明変数となる。そして、各日の１８時における自然エネルギー発電量、日射量又は風速（図の場合、自然エネルギー発電量）が目的変数となる。予測式取得部１１は、このような説明変数及び目的変数の組み合わせからなる複数日分の教師データに基づいた機械学習により得られた予測式を取得する。当該予測式は、任意の日の１８時における自然エネルギー発電量、日射量又は風速を予測する式である。
【００４４】
そして、気象データ取得部１２は、少なくとも、２０１５年１月１日（予測対象日）の６時から１８時までの気象データを取得する。当該気象データは予測値であってもよいし、実績値と予測値とが混在していてもよい。混在の例としては、例えば、６時から１２時までは実績値であり、それ以降は予測値である等が考えられる。
【００４５】
特徴量抽出部１３は、気象データ取得部１２により取得された２０１５年１月１日（予測対象日）の６時から１８時までの気象データから、所定の特徴量を抽出する。当該特徴量は、２０１５年１月１日（予測対象日）の６時から１８時までの時間内における気象データの時系列な変動を表す。
【００４６】
第１の推定部１４は、上述の通り予測式取得部１１が取得した予測式と、上述の通り特徴量抽出部１３が抽出した特徴量とに基づいて、２０１５年１月１日（予測対象日）の１８時における自然エネルギー発電量、日射量又は風速を予測する。
【００４７】
対象時刻を変化させて上記処理を繰り返すことで、２０１５年１月１日（予測対象日）の終日の自然エネルギー発電量、日射量又は風速の予測が得られる。
【００４８】
次に、本実施形態の作用効果について説明する。本実施形態の予測装置１０は、対象時刻よりも所定時間（ｍ時間）前から対象時刻までの気象データの変動の特徴に基づいて、対象時刻における自然エネルギー発電量、日射量又は風速を推定する。以下の実施例に示すように、このような本実施形態によれば、自然エネルギー発電量、日射量又は風速の推定の精度を高めることができる。また、本実施形態の予測装置１０は、複数日分の教師データに基づいた機械学習で、予測式を生成することができる。このため、精度の高い予測式を生成することができる。
【００４９】
＜第２の実施形態＞
本実施形態は、測対象日及び予測対象地点の少なくとも一方が特定された予測対象と所定レベル以上類似する過去データを選択的に用いた機械学習により、推定式を生成する点で、第１の実施形態と異なる。以下、詳細に説明する。
【００５０】
本実施形態の機能ブロック図の一例は、第１の実施形態と同様、図２で示される。図示するように、本実施形態の予測装置１０は、予測式取得部１１と、気象データ取得部１２と、特徴量抽出部１３と、第１の推定部１４とを有する。以下、第１の実施形態と異なる点を説明する。
【００５１】
予測式取得部１１は、予測対象日及び予測対象地点の少なくとも一方が特定された予測対象と、所定の属性が所定レベル以上類似する教師データに基づき生成された予測式を取得する。以下、このような予測式を生成する処理について説明する。
【００５２】
まず、予測式生成部２２は、予測対象の属性値を取得する。上述の通り、予測対象は、予測対象日及び予測対象地点の少なくとも一方が特定されている。このような予測対象の属性値として、例えば、予測対象月、予測対象日の季節、予測対象日の気象データの予測値、予測対象地点等が取得できる。
【００５３】
その後、予測式生成部２２は、過去データ記憶部２１が記憶する過去データの中から、所定の属性が予測対象と所定レベル以上類似するデータを抽出する。例えば、予測対象地点（観測地）が一致する、又は、その差（距離）が所定値以下であるデータを抽出してもよい。その他、季節や月が一致するデータを抽出してもよい。その他、所定の時刻の所定の項目（気象データ）の値が一致、又は、その差が所定値以下のデータを抽出してもよい（予測対象の予測値と、過去データの実績値との比較）。その他、これらの条件を所定の論理式で組み合わせた条件を満たすデータを抽出してもよい。その他、類似度を算出する任意の手法を利用して類似度を算出し、類似度が所定レベル以上であるデータを抽出してもよい。
【００５４】
その後、予測式生成部２２は、抽出したデータを教師データとした機械学習により、予測式を生成する。
【００５５】
気象データ取得部１２は、予測対象の対象時刻までの気象データを取得する。特徴量抽出部１３は、当該気象データから特徴量を抽出する。第１の推定部１４は、当該特徴量と、予測式取得部１１が取得した予測式とに基づいて、予測対象の対象時刻における自然エネルギー発電量、日射量又は風速を推定する。
【００５６】
このような本実施形態によれば、予測装置１０は、予測対象の対象時刻における自然エネルギー発電量、日射量又は風速を推定において、所定の属性が予測対象に所定レベル以上類似した過去データを選択的に教師データとして用いて生成された予測式を利用することができる。
【００５７】
例えば、予測装置１０は、第１の観測地における自然エネルギー発電量、日射量又は風速を推定する場合、第１の観測地の過去データを選択的に教師データとして用いて生成された推定式に基づいて、自然エネルギー発電量、日射量又は風速を推定することができる。
【００５８】
また、予測装置１０は、１０月の任意の日における自然エネルギー発電量、日射量又は風速を推定する場合、１０月の過去データを選択的に教師データとして用いて生成された推定式に基づいて、自然エネルギー発電量、日射量又は風速を推定することができる。
【００５９】
また、予測装置１０は、予測気温（最高気温、最低気温等）がＭ℃となっている日（予測対象日）における自然エネルギー発電量、日射量又は風速を推定する場合、気温（最高気温、最低気温等の実績値）が上記予測気温と所定レベル以上類似する過去データを選択的に教師データとして用いて生成された推定式に基づいて、自然エネルギー発電量、日射量又は風速を推定することができる。
【００６０】
このような本実施形態の予測装置１０によれば、自然エネルギー発電量、日射量又は風速の推定の精度が向上する。
【００６１】
＜第３の実施形態＞
本実施形態の予測装置１０は、ｍの値が可変である点で、第１及び第２の実施形態と異なる。以下、詳細に説明する。
【００６２】
図６に、本実施形態の予測装置１０の機能ブロック図の一例を示す。図示するように、予測装置１０は、予測式取得部１１と、気象データ取得部１２と、特徴量抽出部１３と、第１の推定部１４と、ｍ値設定部１５とを有する。以下、第１及び第２の実施形態と異なる点を説明する。
【００６３】
ｍ値設定部１５は、ｍの値を設定する。例えば、ｍ値設定部１５は、過去データを用いた分析により、最適なｍの値を決定し、決定した値を設定してもよい。例えば、ｍ値設定部１５は、上記分析により、ｍの値各々における推定の精度を算出してもよい。そして、ｍ値設定部１５は、精度が最も高いｍの値を設定してもよい。その他、ｍ値設定部１５は、ユーザからｍの値を指定する入力を受付けてもよい。そして、ｍ値設定部１５は、受付けた値を設定してもよい。例えば、ｍ値設定部１５は、上記分析の結果（ｍの値各々における推定の精度）をユーザに向けて出力する手段と、ユーザからｍの値を指定する入力を受付ける手段とを備えてもよい。
【００６４】
予測式取得部１１は、ｍ値設定部１５により設定されたｍの値に基づいて生成された予測式を取得する。特徴量抽出部１３は、ｍ値設定部１５により設定されたｍの値に基づいて特徴量を抽出する。
【００６５】
ここで、ｍ値設定部１５が、過去データを用いた分析により、ｍの値各々における推定の精度を算出する処理の一例について説明する。ここでは、第１の対象時刻の推定に適したｍの値を決定する処理を説明する。
【００６６】
（１）まず、ｍ値設定部１５は、過去データ記憶部２１に記憶されている過去データの中から、予測式生成部２２による予測式の生成に利用されるデータ（以下、対象データ）を抽出する。
【００６７】
対象データは、例えば、予測対象と所定の属性が所定レベル以上類似するデータ（例：観測地が一致するデータ、季節が一致するデータ、予測対象日の月が一致するデータ、所定の項目の気象データが所定レベル以上類似するデータ等）であってもよいし、又は、予測対象日よりも所定日数前の日から予測対象日の前日までのデータであってもよい。
【００６８】
（２）次に、ｍ値設定部１５は、対象データに基づいて、複数のｍ値（例：１乃至１５）各々に対応した予測式（第１の対象時刻における予測を行う予測式）を生成する。
【００６９】
（３）その後、ｍ値設定部１５は、ｍ値毎に生成した予測式各々に、対象データの中の任意のサンプル日の特徴量（第１の対象時刻よりもｍ時間前から第１の対象時刻までの気象データから抽出した特徴量）を入力し、サンプル日の第１の対象時刻における自然エネルギー発電量、日射量又は風速の予測値を得る。
【００７０】
（４）その後、ｍ値設定部１５は、ｍ値毎に、サンプル日の第１の対象時刻における実績値と、上記（３）で算出した当該サンプル日の第１の対象時刻における予測値との差を算出する。
【００７１】
なお、任意のサンプル日を複数設定し、（３）及び（４）の処理をサンプル日ごとに行ってもよい。このようにすれば、ｍの値毎に複数の上記差が得られる。この場合、ｍ値設定部１５は、複数の差の統計値（例：平均値、最大値、最小値、最頻値、中央値等）を、各ｍ値の差の代表値としてもよい。
【００７２】
このようにして得られた差に基づいて、第１の対象時刻の推定における各ｍ値の推定の精度を評価できる。当該差が小さい程、予測の精度は高いことを意味する。例えば、ｍ値設定部１５は、当該差が最も小さいｍの値を設定してもよい。なお、ｍ値設定部１５は、対象時刻ごとに上記処理を行い、最適なｍ値を設定してもよい。
【００７３】
また、以下の実施例で示すが、本発明者らは、予測対象の属性（観測地点、季節、月、気象等）が異なると、予測精度を高める上で最適なｍの値が異なり得ることを見出した。
【００７４】
例えば、ある観測地点ではｍの値が１０の時に予測精度が最も高くなるが、他の観測地点ではｍの値が１２の時に予測精度が最も高くなる、という現象が起こり得る。同様に、季節、月、気象等に応じて、最適なｍの値が変わり得る。
【００７５】
本実施形態によれば、ｍ値設定部１５は、推定対象に応じた適切な対象データを選択し、観測地毎に（地域毎に）最適なｍの値を設定することができる。すなわち、観測地毎に最適化した推定式を利用することができる。また、ｍ値設定部１５は、予測対象日の属性（季節、月、気象等）に基づいて、予測対象日に最適なｍの値を設定することができる。すなわち、予測対象日毎に最適化した推定式を利用することができる。このような本実施形態によれば、自然エネルギー発電量、日射量又は風速の推定の精度が向上する。
【００７６】
また、推定した自然エネルギー発電量、日射量又は風速の用途によっては、ある程度の推定精度でよく、それよりも推定の処理速度を向上させることを希望する場合がある。ユーザがｍの値を指定できる本実施形態によれば、ユーザは、例えば、予測装置１０により提供された各ｍ値の推定の精度を考慮しながら、自己の用途に合ったｍの値を選択することができる。例えば、推定精度を重視する場合、ユーザは、処理速度が遅くなっても最適なｍの値（推定の精度を高めることができる値）を選択することができる。また、処理速度を重視する場合には、ある程度の推定の精度を得られる任意のｍの値を選択することができる。このように、本実施形態の予測装置１０によれば、ユーザフレンドリーな装置を実現できる。
【００７７】
＜第４の実施形態＞
本実施形態の予測装置１０は、所定の情報をユーザに提供する手段（情報出力部）を備える点で、第１乃至第３の実施形態と異なる。図７乃至図９に、本実施形態の情報出力手段が出力する情報の一例を示す。
【００７８】
図７に示す例では、設定されているパラメータを表示するエリア（パラメータ設定エリア）と、所定のメイン情報が表示されるメインエリア（図７の場合、入力変数Ｘｎの時間変化を示すグラフが表示されているエリア）と、メインエリアに表示する情報の選択内容を表示するエリア（画面切替エリア）と、が表示されている。
【００７９】
パラメータ設定エリアには、設定されている各種パラメータが表示されている。図の例の場合、対象地点（観測地）、対象日（予測対象日）、対象時刻、遡る時間の設定値（ｍの設定値）、入力変数（説明変数）の種類（気象データの１つ又は複数の項目）、学習日数（予測式の生成に利用した教師データの量）が示されている。
【００８０】
画面切替エリアには、メインエリアに表示する情報の選択内容が表示されている。当該エリアには、入力変数、予測値、実績値及びグラフ表示のパラメータが存在し、各々、Ｏｎ又はＯｆｆが対応付けられている。
【００８１】
入力変数及びグラフ表示がＯｎになり、予測値及び実績値がＯｆｆになっている図７の例の場合、メインエリアには、一方の軸に設定されている入力変数をとり、他方の軸に時刻を取ったグラフが表示されている。そして、当該グラフ上に、予測対象日の対象時刻（ｔ）よりもｍ時間前（ｔ−ｍ）から対象時刻（ｔ）までの気象データ（入力変数として設定されている項目の気象データ）から抽出した特徴量（入力変数）が表示されている。なお、設定されている入力変数の種類（気象データの項目）が複数ある場合、図示するようなグラフを並べて一覧表示してもよい。
【００８２】
図８に示す例の情報の骨格は、図７と同様である。図８に示す例の場合、画面切替エリアを観察すると、入力変数及び実績値がＯｎになり、予測値及びグラフ表示がＯｆｆになっている。当該例の場合、メインエリアには、予測式の生成に利用される教師データの値が一覧表示されている。図より、ｐ日分の教師データが表示され、説明変数（Ｘ１（ｔ）・・・）と、目的変数（対象時刻（ｔ）における実績値（自然エネルギー発電量、日射量又は風速））が表示されていることが分かる。
【００８３】
図９に示す例の情報の骨格は、図７及び図８と同様である。図９に示す例の場合、画面切替エリアを観察すると、予測値及びグラフ表示がＯｎになり、入力変数及び実績値がＯｆｆになっている。当該例の場合、メインエリアには、一方の軸に自然エネルギー発電量（実績値及び予測値）をとり、他方の軸に時刻を取ったグラフが表示されている。そして、当該グラフ上に、予測対象日の対象時刻（ｔ）までの自然エネルギー発電量（実績値及び予測値）の値が表示されている。例えば、グラフ上におけるすべての時刻において、第１の推定部１４により推定された予測値が表示されてもよい。その他、グラフ表示の時点までに自然エネルギー発電量の実績値が得られている時刻においては、実績値がプロットされてもよい。そして、実績値が得られていない時刻においては第１の推定部１４により推定された予測値が表示されてもよい。
【００８４】
第３の実施形態で説明したように、ｍ値設定部１５は、上記分析の結果（ｍの値各々における推定の精度）をユーザに向けて出力する手段と、ユーザからｍの値を指定する入力を受付ける手段とを備えることができる。例えば、ｍ値設定部１５は、図７乃至図９に示すような画面（例えばメインエリア）において、上記分析の結果を表示してもよい。そして、ｍ値設定部１５は、図７乃至図９に示すような画面（例えばパラメータ設定エリア）に、ｍの値を指定する入力を受付けるＧＵＩ（graphical user interface）部品を表示し、ｍの値の入力を受付けてもよい。
【００８５】
ｍ値設定部１５は、ｍの値を設定する。例えば、ｍ値設定部１５は、過去データを用いた分析により、最適なｍの値を決定し、決定した値を設定してもよい。例えば、ｍ値設定部１５は、上記分析により、ｍの値各々における推定の精度を算出してもよい。そして、ｍ値設定部１５は、精度が最も高いｍの値を設定してもよい。その他、ｍ値設定部１５は、ユーザからｍの値を指定する入力を受付けてもよい。そして、ｍ値設定部１５は、受付けた値を設定してもよい。例えば、ｍ値設定部１５は、上記分析の結果（ｍの値各々における推定の精度）をユーザに向けて出力する手段と、ユーザからｍの値を指定する入力を受付ける手段とを備えてもよい。
【００８６】
以上説明した本実施形態によれば、推定に用いる入力変数の詳細、推定式に利用された教師データの詳細、及び、推定結果を、所定の表示形式でユーザに向けて出力することができる。このような本実施形態によれば、ユーザは、推定結果のみならず、入力変数や教師データの詳細を確認することで、推定結果の妥当性を判断することができる。
【００８７】
＜第５の実施形態＞
本実施形態の予測装置３０は、予測対象日の対象時刻よりもｎ時間前（ｎは０より大）から対象時刻より所定時間（ｎより小さい時間）前までの実績データ（自然エネルギー発電量、日射量又は風速）に基づいた機械学習により、予測対象日の対象時刻における自然エネルギー発電量、日射量又は風速を推定する。ｎの値は可変である。以下、詳細に説明する。
【００８８】
図１０に、本実施形態の予測装置３０の機能ブロック図の一例を示す。図示するように、予測装置３０は、実績データ取得部３１と、第２の推定部３２と、ｎ値設定部３３とを有する。
【００８９】
実績データ取得部３１は、予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する。実績データ取得部３１は、少なくとも、予測対象日の対象時刻よりもｎ時間前（ｎは０より大）から、対象時刻よりも所定時間（ｎより小さい時間）前までの自然エネルギー発電量、日射量又は風速の実績データを取得する。
【００９０】
第２の推定部３２は、対象時刻よりもｎ時間前（ｎは０より大）から対象時刻より所定時間（ｎより小さい時間）前までの実績データ（自然エネルギー発電量、日射量又は風速）に基づいて、対象時刻における自然エネルギー発電量、日射量又は風速を推定する。推定には、例えば時系列分析のモデルを利用することができる。
【００９１】
ｎ値設定部３３は、ｎの値を設定する。例えば、ｎ値設定部３３は、過去データ記憶部２１が記憶する過去データを用いた分析により、ｎの値各々における推定（第２の推定部３２による推定）の精度を算出してもよい。そして、ｎ値設定部３３は、算出結果に基づいてｎの値を決定し、決定した値を設定してもよい。例えば、ｎ値設定部３３は、精度が最も高いｎの値を設定してもよい。その他、ｎ値設定部３３は、ユーザからｎの値を指定する入力を受付けてもよい。そして、ｎ値設定部３３は、受付けた値を設定してもよい。例えば、ｎ値設定部３３は、上記分析の結果（ｎの値各々における推定の精度）をユーザに向けて出力する手段と、ユーザからｎの値を指定する入力を受付ける手段とを備えてもよい。
【００９２】
ここで、ｎ値設定部３３が行う過去データを用いた分析の例を説明する。ここでは、第１の対象時刻の推定に適したｎの値を決定する処理を説明する。
【００９３】
（１）´まず、ｎ値設定部３３は、過去データ記憶部２１に記憶されている過去データの中から、所定のデータを抽出する。例えば、予測装置３０は、過去データ記憶部２１を有してもよい。または、予測装置３０と通信可能である外部装置が過去データ記憶部２１を有していてもよい。
【００９４】
ｎ値設定部３３は、例えば、予測対象日及び観測地の少なくとも一方が特定された予測対象と所定の属性が所定レベル以上類似するデータ（例：観測地が一致するデータ、季節が一致するデータ、予測対象日の月が一致するデータ、所定の項目の気象データが所定レベル以上類似するデータ等）を抽出してもよいし、又は、予測対象日よりも所定日数前の日から予測対象日の前日までのデータを抽出してもよい。
【００９５】
（２）´その後、ｎ値設定部３３は、抽出したデータを利用して、第１の対象時刻よりもｎ時間前（ｎは０より大）から第１の対象時刻より所定時間（ｎより小さい時間）前までの実績データ（自然エネルギー発電量、日射量又は風速）に基づき、第１の対象時刻の自然エネルギー発電量、日射量又は風速の予測を行う。ここでの予測には、第２の推定部３２が利用するアルゴリズムと同じものを利用する。
【００９６】
（３）´その後、ｎ値設定部３３は、算出した第１の対象時刻の予測値と、当該第１の対象時刻における実績値との差を算出する。なお、複数の日のデータ各々に基づき、各日の上記差を算出してもよい。そして、それらの統計値（例：平均値、最大値、最小値、最頻値、中央値等）を、の差の代表値を算出してもよい。
【００９７】
ｎ値設定部３３は、複数のｎの値各々に対して、上記（２）´及び（３）´の処理を行い、ｎの値毎に上記差を算出する。当該差に基づいて、各ｎ値の推定の精度を評価できる。当該差が小さい程、予測の精度は高いことを意味する。例えば、ｎ値設定部３３は、上記差が最も小さいｎの値を設定してもよい。
【００９８】
第２の推定部３２は、ｎ値設定部３３が設定したｎの値に基づいて、自然エネルギー発電量、日射量又は風速を推定する。
【００９９】
以上説明したように、本実施形態の予測装置３０は、予測対象日の対象時刻よりもｎ時間前（ｎは０より大）から対象時刻より所定時間（ｎより小さい時間）前までの実績データ（自然エネルギー発電量、日射量又は風速）に基づいた機械学習により、予測対象日の対象時刻における自然エネルギー発電量、日射量又は風速を推定する。
【０１００】
そして、ｎの値は、可変である。このような本実施形態の予測装置３０によれば、上記（１）´の処理で最適なデータを選択することで、例えば、観測地毎に最適なｎの値を決定したり、所定の属性（季節、月、気象等）毎に最適なｎの値を決定したりできる。このような本実施形態によれば、自然エネルギー発電量、日射量又は風速の推定の精度が向上する。
【０１０１】
＜実施例＞
以下の条件で、第１乃至第４の実施形態の予測装置１０の検証を行った。
【０１０２】
観測地：札幌、東京
予測対象日：２００８年６月から８月の各日
対象時刻：８時、９時、１０時、１１時、１２時、１３時、１４時、１５時、１６時及び１７時
ｍの値：０から１２各々
教師データ：予測対象日の直前６０日分のデータ
説明変数：対象時刻よりｍ時間前から対象時刻までの上層雲量、中層雲量、下層雲量、気温、湿度各々の１時間毎の値、および、対象時刻とその１時間前の大気外日射量の１時間毎の値
目的変数：対象時刻の日射量
推定式に入力する値：予測対象日の前日の１５時において発表されていた予測対象日の対象時刻よりｍ時間前から対象時刻までの気象データ（上記説明変数の項目）の１時間毎の予測値
予測実行時刻：予測対象日の前日の１８時に翌日の予測を行う
機械学習手法：サポートベクターマシン
【０１０３】
「ｍの値に応じた精度改善率」
まず、ＭＡＰＥ（Mean Absolute Percentage Error）を用いて、各ｍの値の予測誤差を算出した。ｘｉは、各対象時刻の日射量の実績値である。ｙｉは、上記条件で推定した各対象時刻の日射量の推定値である。ｎは、ｍの値各々に対応するサンプル数である。
【０１０４】
【数１】
2016121202.tif 000003
【０１０５】
そして、ｍ＝０の時のＭＡＰＥの値を基準値とし、基準値からｍ値各々のＭＡＰＥの値を引いた差分を、ｍ値各々の精度改善率とした。精度改善率が正の値の場合、ｍ＝０の場合に比べて精度が改善している。そして、その値が大きい程、精度が改善している。一方、精度改善率が負の値の場合、ｍ＝０の場合に比べて精度が悪化している。そして、その値が小さい程、精度が悪化している。
【０１０６】
図１１に、札幌及び東京各々での検証結果を示す。図に示す基準値は、ｍ＝０の時のＭＡＰＥの値を示す。図より、札幌及び東京いずれにおいても、ｍの値が１の場合より、ｍの値が２以上の場合の方が、精度改善率が高いことが分かる。また、札幌及び東京いずれにおいても、ｍの値が所定値以下の場合、ｍの値が大きくなるほど、精度改善率が大きくなる傾向が読み取れる。
【０１０７】
また、ｍの値が所定値を超えると精度改善率が飽和し、それ以上ｍの値が大きくなってもあまり変化しないことが分かる。そして、精度改善率が飽和するｍの値は、観測地毎に異なることが分かる。
【０１０８】
また、札幌で最も精度改善率が大きいｍの値は１２であり、東京で最も精度改善率が大きいｍの値は１０であることが分かる。すなわち、最適なｍの値は、観測地毎に異なることが分かる。
【０１０９】
以下、参考形態の例を付記する。
１．  対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段と、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段と、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段と、
前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段と、
を有する予測装置。
２．  １に記載の予測装置において、
前記ｍの値は可変であり、前記ｍの値を設定するｍ値設定手段をさらに有する予測装置。
３．  ２に記載の予測装置において、
前記第１の推定手段は、複数の地域の自然エネルギー発電量、日射量又は風速を推定し、
前記ｍ値設定手段は、地域毎に前記ｍの値を設定する予測装置。
４．  ２又は３に記載の予測装置において、
前記ｍ値設定手段は、前記予測対象日の属性に基づいて、前記ｍの値を設定する予測装置。
５．  １から４のいずれかに記載の予測装置において、
前記予測式取得手段は、前記予測対象日及び予測対象地点の少なくとも一方が特定された予測対象と、所定の属性が所定レベル以上類似する前記教師データに基づき生成された前記予測式を取得する予測装置。
６．  １から５のいずれかに記載の予測装置において、
前記特徴量は、前記対象時刻よりもｍ時間前から前記対象時刻までの時間内における気象データの変動の特徴を表す予測装置。
７．  コンピュータが、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得工程と、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得工程と、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出工程と、
前記予測式取得工程で取得した前記予測式と、前記特徴量抽出工程で抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定工程と、
を実行する予測方法。
８．  コンピュータを、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段、及び、
前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段、
として機能させるためのプログラム。
９．  予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する実績データ取得手段と、
前記対象時刻よりもｎ時間前（ｎは０より大）から前記対象時刻より前記所定時間前までの前記実績データに基づいて、前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第２の推定手段と、
前記ｎの値は可変であり、前記ｎの値を設定するｎ値設定手段と、
を有する予測装置。
１０．  コンピュータが、
予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する実績データ取得工程と、
前記対象時刻よりもｎ時間前（ｎは０より大）から前記対象時刻より前記所定時間前までの前記実績データに基づいて、前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第２の推定工程と、
前記ｎの値は可変であり、前記ｎの値を設定するｎ値設定工程と、
を実行する予測方法。
１１．  コンピュータを、
予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する実績データ取得手段、
前記対象時刻よりもｎ時間前（ｎは０より大）から前記対象時刻より前記所定時間前までの前記実績データに基づいて、前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第２の推定手段、及び、
前記ｎの値は可変であり、前記ｎの値を設定するｎ値設定手段、
として機能させるためのプログラム。
【０１１０】
この出願は、２０１５年１月３０日に出願された日本出願特願２０１５−０１７１０７号を基礎とする優先権を主張し、その開示の全てをここに取り込む。
【図１】
2016121202.tif 000004
【図２】
2016121202.tif 000005
【図３】
2016121202.tif 000006
【図４】
2016121202.tif 000007
【図５】
2016121202.tif 000008
【図６】
2016121202.tif 000009
【図７】
2016121202.tif 000010
【図８】
2016121202.tif 000011
【図９】
2016121202.tif 000012
【図１０】
2016121202.tif 000013
【図１１】
2016121202.tif 000014
【図１２】
2016121202.tif 000015
【国際調査報告】
2016121202.tif 000016
2016121202.tif 000017