JPWO2016121202A1 - 予測装置、予測方法、及び、プログラム - Google Patents

予測装置、予測方法、及び、プログラム Download PDF

Info

Publication number
JPWO2016121202A1
JPWO2016121202A1 JP2016571785A JP2016571785A JPWO2016121202A1 JP WO2016121202 A1 JPWO2016121202 A1 JP WO2016121202A1 JP 2016571785 A JP2016571785 A JP 2016571785A JP 2016571785 A JP2016571785 A JP 2016571785A JP WO2016121202 A1 JPWO2016121202 A1 JP WO2016121202A1
Authority
JP
Japan
Prior art keywords
prediction
value
target time
amount
target
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2016571785A
Other languages
English (en)
Inventor
鈴木　勝也
勝也 鈴木
康将 本間
康将 本間
耕治 工藤
耕治 工藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
NEC Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NEC Corp filed Critical NEC Corp
Publication of JPWO2016121202A1 publication Critical patent/JPWO2016121202A1/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00—Photometry, e.g. photographic exposure meter
    • G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R21/00—Arrangements for measuring electric power or power factor
    • G01R21/133—Arrangements for measuring electric power or power factor by using digital technique
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00—Machine learning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00—Computing arrangements based on specific mathematical models
    • G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00—Photometry, e.g. photographic exposure meter
    • G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G01J2001/4266—Photometry, e.g. photographic exposure meter using electric radiation detectors for measuring solar light
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00—Machine learning
    • G06N20/10—Machine learning using kernel methods, e.g. support vector machines [SVM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Environmental Sciences (AREA)
  • Ecology (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Atmospheric Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Artificial Intelligence (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本発明は、機械学習による統計的手法を用いて自然エネルギー発電量、日射量又は風速を予測する技術において、予測の精度を高めることを課題とする。当該課題を解決するため、対象時刻よりもｍ時間前（ｍは２以上）から対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出部（１３）と、複数日分の特徴量に基づいて対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定部（第１の推定部（１４））と、を有する予測装置（１０）を提供する。

Description

本発明は、予測装置、予測方法、及び、プログラムに関し、より具体的には、自然エネルギー発電量、日射量及び／又は風速を予測する予測装置、予測方法、及び、プログラムに関する。
特許文献１乃至３、及び、非特許文献１に、機械学習による統計的手法を用いて、気象データから太陽光発電の発電量や日射量を予測する技術が開示されている。
特開平９−２１５１９２号 特許第３９８４６０４号 特許第５３３９３１７号
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）を用いて、各ｍの値の予測誤差を算出した。ｘｉは、各対象時刻の日射量の実績値である。ｙｉは、上記条件で推定した各対象時刻の日射量の推定値である。ｎは、ｍの値各々に対応するサンプル数である。
Figure 2016121202
そして、ｍ＝０の時のＭＡＰＥの値を基準値とし、基準値からｍ値各々のＭＡＰＥの値を引いた差分を、ｍ値各々の精度改善率とした。精度改善率が正の値の場合、ｍ＝０の場合に比べて精度が改善している。そして、その値が大きい程、精度が改善している。一方、精度改善率が負の値の場合、ｍ＝０の場合に比べて精度が悪化している。そして、その値が小さい程、精度が悪化している。
図１１に、札幌及び東京各々での検証結果を示す。図に示す基準値は、ｍ＝０の時のＭＡＰＥの値を示す。図より、札幌及び東京いずれにおいても、ｍの値が１の場合より、ｍの値が２以上の場合の方が、精度改善率が高いことが分かる。また、札幌及び東京いずれにおいても、ｍの値が所定値以下の場合、ｍの値が大きくなるほど、精度改善率が大きくなる傾向が読み取れる。
また、ｍの値が所定値を超えると精度改善率が飽和し、それ以上ｍの値が大きくなってもあまり変化しないことが分かる。そして、精度改善率が飽和するｍの値は、観測地毎に異なることが分かる。
また、札幌で最も精度改善率が大きいｍの値は１２であり、東京で最も精度改善率が大きいｍの値は１０であることが分かる。すなわち、最適なｍの値は、観測地毎に異なることが分かる。
以下、参考形態の例を付記する。
１． 対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段と、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段と、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段と、
前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段と、
を有する予測装置。
２． １に記載の予測装置において、
前記ｍの値は可変であり、前記ｍの値を設定するｍ値設定手段をさらに有する予測装置。
３． ２に記載の予測装置において、
前記第１の推定手段は、複数の地域の自然エネルギー発電量、日射量又は風速を推定し、
前記ｍ値設定手段は、地域毎に前記ｍの値を設定する予測装置。
４． ２又は３に記載の予測装置において、
前記ｍ値設定手段は、前記予測対象日の属性に基づいて、前記ｍの値を設定する予測装置。
５． １から４のいずれかに記載の予測装置において、
前記予測式取得手段は、前記予測対象日及び予測対象地点の少なくとも一方が特定された予測対象と、所定の属性が所定レベル以上類似する前記教師データに基づき生成された前記予測式を取得する予測装置。
６． １から５のいずれかに記載の予測装置において、
前記特徴量は、前記対象時刻よりもｍ時間前から前記対象時刻までの時間内における気象データの変動の特徴を表す予測装置。
７． コンピュータが、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得工程と、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得工程と、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出工程と、
前記予測式取得工程で取得した前記予測式と、前記特徴量抽出工程で抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定工程と、
を実行する予測方法。
８． コンピュータを、
対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段、
予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段、
前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段、及び、
前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段、
として機能させるためのプログラム。
９． 予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する実績データ取得手段と、
前記対象時刻よりもｎ時間前（ｎは０より大）から前記対象時刻より前記所定時間前までの前記実績データに基づいて、前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第２の推定手段と、
前記ｎの値は可変であり、前記ｎの値を設定するｎ値設定手段と、
を有する予測装置。
１０． コンピュータが、
予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する実績データ取得工程と、
前記対象時刻よりもｎ時間前（ｎは０より大）から前記対象時刻より前記所定時間前までの前記実績データに基づいて、前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第２の推定工程と、
前記ｎの値は可変であり、前記ｎの値を設定するｎ値設定工程と、
を実行する予測方法。
１１． コンピュータを、
予測対象日の対象時刻より所定時間前までの自然エネルギー発電量、日射量又は風速の実績データを取得する実績データ取得手段、
前記対象時刻よりもｎ時間前（ｎは０より大）から前記対象時刻より前記所定時間前までの前記実績データに基づいて、前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第２の推定手段、及び、
前記ｎの値は可変であり、前記ｎの値を設定するｎ値設定手段、
として機能させるためのプログラム。
この出願は、２０１５年１月３０日に出願された日本出願特願２０１５−０１７１０７号を基礎とする優先権を主張し、その開示の全てをここに取り込む。

Claims (13)

  1. 対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段と、
    複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段と、を有する予測装置。
  2. 前記推定手段は、
    前記対象時刻までの気象データから抽出した特徴量を説明変数、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数として予測する予測式を用いて推定を行う請求項１に記載の予測装置。
  3. 前記推定手段は、前記説明変数および前記目的変数の組み合わせからなる複数日分の教師データに基づいた予測式を用いて推定を行う請求項２に記載の予測装置。
  4. 前記推定手段は、予測対象日と所定の属性が類似する複数日分の前記特徴量に基づいて推定を行う請求項１から３のいずれか１項に記載の予測装置。
  5. 前記ｍの値は可変であり、前記ｍの値を設定するｍ値設定手段をさらに有する請求項１から４のいずれか１項に記載の予測装置。
  6. 前記推定手段は、複数の地域の自然エネルギー発電量、日射量又は風速を推定し、
    前記ｍ値設定手段は、地域毎に前記ｍの値を設定する請求項５に記載の予測装置。
  7. 前記ｍ値設定手段は、予測対象日の属性に基づいて、前記ｍの値を設定する請求項５又は６に記載の予測装置。
  8. 前記ｍ値設定手段は、ｍの値各々における推定の精度を出力する手段を備える請求項５から７のいずれか１項に記載の予測装置。
  9. 前記推定手段は、予測対象日及び予測対象地点の少なくとも一方が特定された予測対象と、所定の属性が所定レベル以上類似する教師データに基づき生成された前記予測式を用いて推定を行う請求項１から８のいずれか１項に記載の予測装置。
  10. 対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段と、
    前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段と、を有する予測装置。
  11. 対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから抽出した特徴量を説明変数とし、前記対象時刻の自然エネルギー発電量、日射量又は風速を目的変数とし、複数日分の教師データに基づいた機械学習で生成された前記対象時刻の自然エネルギー発電量、日射量又は風速を予測する予測式を取得する予測式取得手段と、
    予測対象日の前記対象時刻までの気象データを取得する気象データ取得手段と、
    前記予測対象日の前記対象時刻よりもｍ時間前から前記対象時刻までの気象データから前記特徴量を抽出する特徴量抽出手段と、
    前記予測式取得手段が取得した前記予測式と、前記特徴量抽出手段が抽出した前記特徴量とに基づいて、前記予測対象日の前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する第１の推定手段と、
    を有する予測装置。
  12. コンピュータが、
    対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出工程と、
    複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定工程と、
    を実行する予測方法。
  13. コンピュータを、
    対象時刻よりもｍ時間前（ｍは２以上）から前記対象時刻までの気象データから時系列の変動である特徴量を抽出する特徴量抽出手段、
    複数日分の前記特徴量に基づいて前記対象時刻における自然エネルギー発電量、日射量又は風速を推定する推定手段、
    として機能させるためのプログラム。
JP2016571785A 2015-01-30 2015-11-18 予測装置、予測方法、及び、プログラム Pending JPWO2016121202A1 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2015017107 2015-01-30
JP2015017107 2015-01-30
PCT/JP2015/082353 WO2016121202A1 (ja) 2015-01-30 2015-11-18 予測装置、予測方法、及び、プログラム

Publications (1)

Publication Number Publication Date
JPWO2016121202A1 true JPWO2016121202A1 (ja) 2017-11-09

Family

ID=56542840

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016571785A Pending JPWO2016121202A1 (ja) 2015-01-30 2015-11-18 予測装置、予測方法、及び、プログラム

Country Status (3)

Country Link
US (1) US20170371073A1 (ja)
JP (1) JPWO2016121202A1 (ja)
WO (1) WO2016121202A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111582551A (zh) * 2020-04-15 2020-08-25 中南大学 风电场短期风速预测方法、系统及电子设备
CN115829140A (zh) * 2022-12-14 2023-03-21 广西电网有限责任公司 一种基于机器学习的风力发电场发电量预测方法及系统

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6628700B2 (ja) * 2016-08-05 2020-01-15 三菱電機株式会社 気象情報予測装置および電力需要予測装置
JP6848271B2 (ja) * 2016-09-01 2021-03-24 富士電機株式会社 風速予測装置、風速予測システム、風速予測方法及びプログラム
KR102407597B1 (ko) * 2017-12-14 2022-06-13 한국전기연구원 단조성을 반영한 일사량 예측 모델 생성 방법
JP2019157841A (ja) * 2018-03-16 2019-09-19 株式会社Ｌｉｘｉｌ 制御装置、及び制御方法
US10955584B2 (en) * 2018-04-25 2021-03-23 Microsoft Technology Licensing, Llc Predicting microclimate
CN109034478B (zh) * 2018-07-27 2021-09-24 中南大学 一种高速铁路沿线大风迭代竞争高精度预测方法
CN109886464B (zh) * 2019-01-20 2022-03-18 东北电力大学 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法
JP2020166622A (ja) * 2019-03-29 2020-10-08 三菱重工業株式会社 発電量予測装置、発電量予測方法、およびプログラム
JP7355548B2 (ja) * 2019-08-07 2023-10-03 三菱重工サーマルシステムズ株式会社 気象観測所選択装置、予測システム、空調システム、気象観測所選択方法およびプログラム
CN111091138A (zh) * 2019-11-14 2020-05-01 远景智能国际私人投资有限公司 辐照预报的处理方法、堆叠泛化模型的训练方法及装置
JP7409080B2 (ja) * 2019-12-27 2024-01-09 富士通株式会社 学習データ生成方法、学習データ生成プログラムおよび情報処理装置
WO2021190764A1 (en) * 2020-03-27 2021-09-30 Winji Ag Device, method and computer program for evaluating the operation of a power plant
US20230273344A1 (en) * 2020-07-31 2023-08-31 Mitsui Chemicals, Inc. Insolation correction method, insolation correction device, recording medium, model, model generating method, and model providing method
CN117494573B (zh) * 2023-11-16 2024-06-11 中山大学 一种风速预测方法、系统及电子设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003108909A (ja) * 2001-09-28 2003-04-11 Tohoku Electric Power Co Inc 短期予測システム
JP3984604B2 (ja) * 2004-07-12 2007-10-03 日本電信電話株式会社 太陽光発電システムの発電量予測方法、装置、およびプログラム
JP2007056686A (ja) * 2005-08-22 2007-03-08 Univ Of Ryukyus 風速予測に基づく風力発電機の数時間先発電電力予測装置
JP5606114B2 (ja) * 2010-03-19 2014-10-15 株式会社東芝 発電量予測装置、予測方法及び予測プログラム
US20140278163A1 (en) * 2013-03-15 2014-09-18 Gigawatt, Inc. Distributed solar power generation and monitoring system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111582551A (zh) * 2020-04-15 2020-08-25 中南大学 风电场短期风速预测方法、系统及电子设备
CN111582551B (zh) * 2020-04-15 2023-12-08 中南大学 风电场短期风速预测方法、系统及电子设备
CN115829140A (zh) * 2022-12-14 2023-03-21 广西电网有限责任公司 一种基于机器学习的风力发电场发电量预测方法及系统
CN115829140B (zh) * 2022-12-14 2023-07-25 广西电网有限责任公司 一种基于机器学习的风力发电场发电量预测方法及系统

Also Published As

Publication number Publication date
WO2016121202A1 (ja) 2016-08-04
US20170371073A1 (en) 2017-12-28

Similar Documents

Publication Publication Date Title
WO2016121202A1 (ja) 予測装置、予測方法、及び、プログラム
Solomon et al. Distinguishing the roles of natural and anthropogenically forced decadal climate variability: implications for prediction
Pi et al. A grey prediction approach to forecasting energy demand in China
Jury et al. Evaluation of CMIP5 models in the context of dynamical downscaling over Europe
JP5661594B2 (ja) 予測モデル構築装置、方法、及びプログラム、並びに発電量予測装置、及び方法
KR101827908B1 (ko) 해양 기상 분석 장치 및 방법
Corchado et al. A topology-preserving system for environmental models forecasting
JP2009225550A (ja) 電力の需要予測処理方法及び装置、並びに発電予測処理方法
Wehner et al. Evaluation of extreme sub-daily precipitation in high-resolution global climate model simulations
Potter et al. Visualization of uncertainty and ensemble data: Exploration of climate modeling and weather forecast data with integrated ViSUS-CDAT systems
JP2016224566A (ja) 予測装置、予測方法および予測プログラム
JP2017010111A (ja) ノイズデータ除去支援装置、方法、および、プログラム
JP2019087030A (ja) 予測モデル生成装置、予測モデル生成方法及び予測モデル生成プログラム
Brentan et al. Water demand time series generation for distribution network modeling and water demand forecasting
JP2021039739A (ja) 予測装置、予測プログラム、及び予測方法
Yen et al. Solar power prediction via support vector machine and random forest
Karozis et al. A deep learning approach for spatial error correction of numerical seasonal weather prediction simulation data
JP6494258B2 (ja) 予測システム、予測方法、および予測プログラム
Marshall et al. Visualizing and verifying probabilistic forecasts of the Madden‐Julian Oscillation
Plant A new modelling framework for statistical cumulus dynamics
Sun et al. Precipitation patterns and associated hydrological extremes in the Yangtze River basin, China, using TRMM/PR data and EOF analysis
Branstator et al. Is AMOC more predictable than North Atlantic heat content?
Durrant et al. Consensus forecasts of modeled wave parameters
Kim et al. Sampling error of observation impact statistics
Wang et al. Reconstruction of three-dimensional Pacific temperature with Argo and satellite observations

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181005

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190813

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20200303