JP5606114B2 - 発電量予測装置、予測方法及び予測プログラム - Google Patents

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

Info

Publication number
JP5606114B2
JP5606114B2 JP2010065083A JP2010065083A JP5606114B2 JP 5606114 B2 JP5606114 B2 JP 5606114B2 JP 2010065083 A JP2010065083 A JP 2010065083A JP 2010065083 A JP2010065083 A JP 2010065083A JP 5606114 B2 JP5606114 B2 JP 5606114B2
Authority
JP
Japan
Prior art keywords
power generation
generation amount
data
predicted value
prediction
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.)
Active
Application number
JP2010065083A
Other languages
English (en)
Other versions
JP2011200040A (ja
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.)
Toshiba Corp
Original Assignee
Toshiba 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
Priority to JP2010065083A priority Critical patent/JP5606114B2/ja
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to MYPI2012700591A priority patent/MY156889A/en
Priority to TW100109378A priority patent/TWI533251B/zh
Priority to PCT/JP2011/001648 priority patent/WO2011114752A1/ja
Priority to CN201180013084.7A priority patent/CN102792542B/zh
Priority to EP11755942.7A priority patent/EP2549641A4/en
Priority to US13/583,792 priority patent/US20130013233A1/en
Publication of JP2011200040A publication Critical patent/JP2011200040A/ja
Application granted granted Critical
Publication of JP5606114B2 publication Critical patent/JP5606114B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/004Generation forecast, e.g. methods or systems for forecasting future energy generation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2260/00Function
    • F05B2260/82Forecasts
    • F05B2260/821Parameter estimation or prediction
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems 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/20Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Description

本発明は、例えば、風力発電や太陽光発電などの自然エネルギーを用いた発電機の発電量を予測するための発電量予測装置、予測方法及び予測プログラムに関する。
風力や太陽光などの自然エネルギーによる発電は、燃料費が不要で、温暖化ガスを排出しないという利点がある。しかし、一般的に、自然エネルギーは天候に左右されるため、その発電量は、不確実で予測が困難であるという欠点がある。これまでは、主に、気象予報の結果を用いて、自然エネルギーによる発電量を予測していた。この分野における先行技術としては、例えば、特許文献1がある。
特開2004−289918号公報
村田晃伸、山口浩、大谷謙仁、「広域的に多数台導入された太陽光発電に関する出力変動幅の推定法」、電学論B、127巻5号、2007年、pp.645−652
ところで、特許文献1における予測技術では、自然エネルギーによる発電量を気象情報から予測し、不足分は発電で補うとしている。しかし、かかる予測には、変動幅があることが考慮されていない。
一方、自然エネルギーを利用した発電機が複数あるときには、相互の変動が重なり合って、変動が小さくなることが知られている(非特許文献1参照)。従って、複数の自然エネルギーによる発電機の発電量の合計を予測するためには、相互の相関を考慮して予測値を合計する必要がある。
例えば、非特許文献1では、相関を考慮した合計出力の変動を、近似式で予測している。しかし、かかる手法では、誤差が大きいという問題があった。また、従来技術では、特定の発電所の合計出力を、相関を考慮して予測することも困難であった。
一方、需要量と発電量とのバランスが正確にとれていないと、周波数が所定値から変動するなど、電力系統が不安定になる。このため、自然エネルギーによる発電を導入した場合には、電力系統の管理を行う上で、特別の注意が必要となる。
ところが、周波数の変動は、発電量の変動と関係がある。そして、上述のように、風力発電や太陽光発電のような自然エネルギーによる発電は、発電機の発電量が不確実に変動する。このため、従来技術による予測方法では、この変動を考慮して、電力系統を安定して管理することが困難であった。
風力発電機の発電量は、風向と風速が決まれば、ほぼ予測することができる。例えば、図16は、ある年の8月に、東京のある地点の風速を1時間ごとに示したものである。また、図17は、同じ地点の風向を1時間ごとに示したものである。図17の風向を表示するにあたっては、図18のように方向を数値化してある。
これらの図16、図17を一見すると、風速や風向はランダムに変動しているように見える。しかし、よく見ると、各日で一定の傾向があることが分かる。例えば、この時期のこの場所での風速は、午前に小さく午後に大きくなり、風速が大きいほど変動も大きい。一方、風向は、風速が小さい場合には大きく変わるが、風速が大きな場合には一定になる傾向がある。
また、太陽光発電の発電量は、場所と発電機の設置方法と日射量が決まれば、ほぼ予測することができる。例えば、図19は、ある年の12月の東京のある地点での日射量を示している。この図19からも、日射量は全くランダムというわけではなく、ある曲線に近いデータが多く、午前中の日射量が小さい場合には午後の日射量も小さいといった特徴があることが分かる。
このように、風力データ、日射量データの両者に共通することは、どちらも全くランダムというわけではないが、完全に予測することは困難であるということである。従って、一定の誤差を考慮して予測するしかない。また、予測にあたっては、過去のデータを参考にして、一定の傾向を満たすように予測することが必要である。この傾向は、場所や時期によって決まる。また、午前と午後、あるいは時間ごとの関係も重要である。この時間ごとの関係を時刻相関と呼ぶ。
風力発電機でも太陽光発電機でも、発電量が気象条件に左右されるため、同じ場所にある発電機の発電量は似通ってくるが、遠くにある発電機の発電量は全く異なるパターンを示す。数km程度であれば、気象条件は似通っていると考えられるため、発電量のパターンには類似性が見られるはずである。これを位置相関と呼ぶ。これらの発電機は、仮に、隣り合う場所に設置されていたとしても、風速や日射量が正確に同じになるわけではないため、位置相関を考慮する必要がある。
以上に述べたように、風力発電機でも太陽光発電機でも1台の発電機の発電量を予測する場合には、時刻相関を考慮する必要がある。また、複数台の発電機の発電量を同時に予測する場合には、位置相関を考慮する必要がある。しかし、これまでの予測手法では、時刻相関や位置相関を考慮した予測を行うことは、容易ではなかった。
本発明は、上記のような従来技術の問題点を解決するために提案されたものであり、その目的は、発電量が不確実に変動する発電機の場合であっても、統計的相関を考慮して、発電量を予測することができる発電量予測装置、予測方法及び予測プログラムを提供することにある。
上記目的を達成するために、本発明は、発電機の発電量に関する予測値を求める発電量予測装置において、少なくとも一つの発電機の過去の発電量に関する過去データを記憶する過去データ記憶部と、前記過去データにおける異なる時刻の間の統計的相関若しくは異なる発電機の位置の間の統計的相関に基づいて、発電機の発電量に関する予測値を、発生確率を含む時系列データとして算出する予測値算出部と、を有し、前記予測値算出部は、前記過去データに基づいて、分散共分散行列を作成する分散共分散行列作成部と、前記分散共分散行列に基づいて、確率分布に従う乱数を発生させる乱数発生部と、を有することを特徴とする。
以上のような本発明では、発電量に関する過去データについて、異なる時刻の間の統計的相関を考慮して、予測値を算出するので、特定の発電機の発電量を、精度よく予測することができる。また、異なる発電機の位置の間の統計的相関を考慮して、予測値を算出するので、複数台の発電機の発電量を、精度よく予測することができる。
以上のような本発明によれば、発電量が不確実に変動する発電機の場合であっても、統計的相関を考慮して、発電量を予測可能な発電量予測装置、予測方法及び予測プログラムを提供することができる。
本発明の発電量予測装置の一実施形態を示すブロック図 図1の実施形態において、風速・風向データに基づいて風力発電の発電量を予測する処理を示すフローチャート 図1の実施形態において、発電量データに基づいて風力発電の発電量を予測する処理を示すフローチャート 太陽光発電の過去の類似日の発電量データの一例を示す説明図 図1の発電量データを重ね書いた一例を示す説明図 図1の発電量データに基づく太陽光発電の発電量予測値のシミュレーション結果を示す説明図 図6の相関係数行列の一例を示す説明図 風力発電の過去の風速データ(a)、風向データ(b)、風速・風向データを重ね書いたもの(c)を示す説明図 風速・風向を含む時系列の予測結果(a)、風速の予測結果(b)、風向の予測結果(c)を示す説明図 2箇所の太陽光発電所の過去の発電量データ(a)(b)、両者を含む時系列データ(c)を示す説明図 2箇所の太陽光発電所を含む発電量の予測値(a)、各々の太陽光発電所の発電量の予測値(b)(c)を示す説明図 図11(b)の発電出力の頻度分布(a)、確率密度分布(b)、図11(c)の発電出力の頻度分布(c)、確率密度分布(d)を示す説明図 図11(b)と図11(c)の発電量の合計を示す説明図 従来の予測モデルの一例を示す概念図 本発明を、従来の予測モデルに適用した一例を示す概念図 ある年の東京の8月の風速データの一例を示す説明図 ある年の東京の8月の風向データの一例を示す説明図 風向データの換算表の一例を示す説明図 ある年の東京の12月の日照量の一例を示す説明図
本発明を実施するための形態(以下、実施形態とする)を、図面を参照して説明する。
[1.実施形態の構成]
まず、本実施形態の構成を説明する。すなわち、図1に示すように、本実施形態の発電量予測装置1は、パーソナルコンピュータやサーバ装置のような汎用のコンピュータを、所定のプログラムで制御することによって実現できる。この場合のプログラムは、コンピュータのハードウェアを物理的に活用することで、以下に述べるような各部の処理を実現するものであり、かかる処理を行う方法、プログラム及びプログラムを記録した記録媒体は単独でも本発明の一態様である。
発電量予測装置1は、予測処理部10、記憶部20、入力部30及び出力部40等を有している。予測処理部10は、上記のようなプログラムの働きにより機能する類似日選択部11、時系列データ生成部12、予測値算出部13、発電量合計部14等を有している。
類似日選択部11は、予測対象日の気象予報データ等と、過去のデータ(風速・風向データ、発電量データ、日射量データ等)とに基づいて、予測対象日と類似した日を選択する処理部である。時系列データ生成部12は、選択された各類似日におけるデータを時系列で抽出する処理部である。
予測値算出部13は、各発電機の発電量の予測値を算出する処理部である。発電量合計部14は、各発電機の発電量の予測値に基づいて、合計発電量の予測値を算出する処理部である。
なお、予測値算出部13は、分散共分散行列作成部131、乱数発生部132、発電量計算部133を有している。分散共分散行列作成部131は、抽出された時系列のデータに基づいて、分散共分散行列を作成する処理部である。乱数発生部132は、作成された分散共分散行列に基づく乱数を発生することにより、予測値(風速・風向データ、発電量データ、日射量データ等)を作成する処理部である。発電量計算部133は、予測値(風速・風向データ、日射量データ等)に基づいて、各発電機の発電量の予測値を算出する処理部である。
記憶部20は、予測処理部10における処理に必要なデータ及び演算された予測値等、各種の情報を記憶する構成部である。これらの情報は、下記の入力部30から入力されるか、外部から通信ネットワークを介して入力され、記憶される。記憶部20は、典型的には、内蔵された若しくは外部接続された各種メモリ、ハードディスク、光ディスク等により構成できるが、現在又は将来において利用可能なあらゆる記憶媒体を利用可能である。
入力部30は、発電量予測処理装置1に必要な情報の入力、処理の選択や指示を入力する構成部である。この入力部30としては、例えば、キーボード、マウス、タッチパネル(表示装置に構成されたものを含む)等が考えられる。但し、現在又は将来において利用可能なあらゆる入力装置を含む。したがって、音声により入力する音声入力装置も含まれる。また、単純なスイッチ等であってもよい。
出力部40は、各種データ、演算された予測値等を、管理者、操作者等も含むユーザが認識可能となるように出力する手段である。この出力部40としては、例えば、表示装置、プリンタ等が考えられる。但し、現在又は将来において利用可能なあらゆる出力装置を含む。したがって、音声により出力する音声出力装置も含まれる。
[2.実施形態の作用]
[2−1.処理の概要]
本実施形態は、気象データ若しくは発電量データを、多次元正規分布を用いてモデル化し、分散共分散行列を用いて多次元正規分布に従う乱数を多数発生させることにより、発生確率が付与された複数の時系列データの組を生成することにより、予測値を算出するものである。
[2−2.風速・風向データに基づく予測]
次に、本実施形態の作用を説明する。まず、風速・風向データに基づいて、風力発電の発電機の発電量を予測する一例を、図2を参照して説明する。なお、図2は、風速・風向データを予測し、予測した風速・風向データから、発電量の予測値を計算する手順とデータの流れを示すフローチャートである。
まず、類似日選択部11が、予測対象日の気象予報データと、各発電機の位置における過去の風速・風向データ(D1−1)とに基づいて、予測対象日と類似した日を、N日分選択する(ステップ101)。過去のデータから、該当日と類似する日を選択する手法としては、特定のものには限定されない。
例えば、午前中が晴れで午後が曇りの日というように、ある情報が一致する日を選択するといった簡単な手法が考えられる。また、晴天指数や雲量など、ある情報について近似する日を選択するといった手法なども用いることができる。
時系列データ生成部12は、選択されたN日分の類似日における時系列データ(D1−2)として、風速と風向の2種のデータの組(K時点、N日分、M台分)を生成する(ステップ102)。
分散共分散行列作成部131は、これらの時系列データ(D1−2)から、(2×K×M)行(2×K×M)列の分散共分散行列を作成する(ステップ103)。次に、乱数発生部132が、分散共分散行列を用いて、(2×K×M)個の乱数の組をL組発生させることによって、各風力発電機の位置における風向と風速の2種の予測値(D1−3)が、L組のデータ(K時点、M台、L通り)として作成される(ステップ104)。なお、このような分散共分散行列の作成処理、乱数の発生処理の詳細については、後述する。
しかる後に、発電量計算部133が、風速および風向の予測値(D1−3)から、K時点における各発電機の発電量の予測値(D1−4)をL組(K時点、M台)生成する(ステップ105)。このL組のデータは元のN日分のデータと同じ時刻相関、位置相関を持ったデータであり、発電量の予測結果に相当する。このため、発電量合計部14が、各発電機の発電量の予測値(D1−4)を合計することにより、全体の発電量の予測値(D1−5)が求められる(ステップ106)
[2−3.日射量データに基づく予測]
本実施形態では、風向・風速データの代わりに、日射量データを用いて、太陽光発電の発電機の発電量を予測することもできる。その手順は、基本的には、上記と同様である。すなわち、気象予報データと過去の日射量データに基づく類似日の選択、日射量の時系列データの生成、分散共分散行列の作成、乱数の発生が行われる。
この場合、M台の発電機のK時点のデータを用いると、(K×M)行(K×M)列の分散共分散行列が得られ、最終的には(K×M)個の乱数の組がL組得られる。これが、日射量の予測値となり、この予測値に基づいて、発電量の予測値及び合計値が求められる。
[2−4.発電量に基づく予測]
さらに、上記の例では、風向・風速の予測値や日射量の予測値を、発電量の予測値へ変換(発電量計算部133による発電量の予測値の計算)する処理が行われていた。しかし、このような変換処理は、必ずしも必要ではない。すなわち、本実施形態においては、過去の発電量データに基づいて、発電量の予測値を求めることもできる。
このように発電量データを直接求める例を、図3のフローチャートを参照して説明する。まず、類似日選択部11が、予測対象日の気象予報データと、過去の発電量データ(D2−1)とに基づいて、同様な条件を持つ過去の日(例えば、上述の類似日と同じ手法)をN日選択する(ステップ201)。
時系列データ生成部12は、選択されたN日分の類似日における時系列データ(D2−2)として、各発電機の発電量データの組(K時点、N日分、M台分)を抽出する(ステップ202)。分散共分散行列作成部131は、これらの時系列データ(D2−2)を統計的に分析し、分散共分散行列を作成する(ステップ203)。
次に、乱数発生部132が、分散共分散行列を用いて、乱数の組をL組生成することによって、各発電機の発電量の予測値(D2−3)を求める(ステップ204)。さらに、発電量合計部14は、各発電機の発電量の予測値(D2−3)を合計することにより、全体の発電量の予測値(D2−4)を求める(ステップ205)。
このように、風向・風速の予測値や日射量の予測値からの変換を経ることなく、発電量の予測値(K時点、M台、L通り)と、その合計値(K時点、L通り)を得ることができる。なお、図3は、風力発電機の場合であるが、太陽光発電機の場合も同様に適用できる。
[2−5.分散共分散行列の作成と乱数の発生]
次に、分散共分散行列の作成と乱数の発生の手法を、太陽光発電機の発電量の予測を一例として説明する。ここでは、K時点を、1時間ごとの24点とし、1台の発電機(M=1)の発電量を、時刻相関を考慮して予測する手法を例にとって説明する。
すなわち、あらかじめ複数の日から各時刻ごとの発電量データを集めたデータ{xi}が、記憶部20に記憶されている。類似日選択部11は、例えば、分類アルゴリズム等を用いることにより、類似日を選択する。時系列データ生成部12は、選択された類似日の発電量を、1時間ごとの時刻iの発電量xiに分解し、これをx1からx24まで並べたベクトルである時系列データを作成する。
この時系列データから、分散共分散行列作成部131は、時刻間の発電量の関係を考慮するために、[数1]のような24×24の成分を持つ分散共分散行列を作成する。ここで、var(xi)は{xi}の分散であり、cov(xi,xj)は、{xi}と{xj}の共分散である。分散や共分散を計算するためには、発電量データを取得する日は最低でも2日以上必要である。発電量データを取得する日数は2日でもかまわないが、日数が多いほど予測精度は高くなる。
Figure 0005606114
乱数発生部132は、分散共分散行列に基づいて、24次元の乱数を多数発生し、発電量の予測値24を作成する。
この場合、各時刻の発電量が従う確率分布を定める必要がある。ここでは、24次元の正規分布を用いた場合を一例として説明する。すなわち、24次元のべクトルx=(x1,x2,x3,...,x24)が、24次元正規分布に従う場合、xが満たす確率密度関数は、xの平均値ベクトルμ=(μ1,μ2,μ3,...μ24) を用いて、[数2]のように表される。
Figure 0005606114
一日分の発電量データは24次元のベクトルxとみなすことができる。ここでは、これを発電量シナリオと呼ぶ。[数2]の確率分布に従う24個の乱数の組を多数発生させることにより、任意の数の発電量シナリオを作成することができる。この発電量シナリオは、過去の発電量データの統計的な性質を再現している。
すなわち、各時刻の変動幅は各時刻の分散によって再現され、時刻間の相関も共分散によって再現されている。しかも、過去の発電量データがわずか2日分しかなくても、例えば百万個の発電量シナリオを作成することもできる。
次に、特定の確率分布に従い、かつ、相互に与えられた相関を持つ乱数の組をL組だけ作成する手法を以下に説明する。ここでは、互いに相関を持つ24個の乱数の組をL組生成する場合を例にとって説明する。
まず、各時刻ごとに過去のデータの分布を求め、その分布に従う乱数をL個発生する。特定の分布に従う乱数を作成する最も簡単な方法の一例は、棄却法によるものである。これは、分布関数の定義域と値域によって定められる領域で一様な乱数を生成し、その中から分布関数の値を超えるものを棄却する方法である。
この方法は、あらゆる分布関数に用いることができるが、必ずしも効率はよくない。また、定義域が有限の領域でない場合にも、適用しにくい。もし、特定の分布の逆関数が存在すれば、一様乱数から容易に特定の分布に従う乱数を作成することができる。これを、24時間の各時刻において繰り返す。
以上のプロセスにより、L個の乱数が、24時間分作成される。この段階では、各時刻の乱数の組は、他の時刻の乱数の組と相関がなく、互いに独立である。もっとも、通常のアルゴリズムで生成される乱数は、疑似乱数であり、独立な乱数を作成したつもりでも、実際には相関がある可能性がある。従って、必要に応じて主成分分析により互いに独立な乱数の組に変換しておく。
この結果、L行24列(L×24)の成分を持つ乱数行列が得られる。この段階で各乱数は標準偏差1、平均0になるように変数変換により正規化しておく。このようにして、正規乱数行列Gが得られる。
次に、24×24の相関係数行列Rをコレスキー分解し、上三角行列Tと下三角行列Tに分解する。相関係数行列は分散共分散行列Σを規格化することで得られる。
[数3]
R=T
その後、正規乱数行列Gに上三角行列Tを右から掛けることで、相関を持った正規乱数行列が作成できる。
[数4]
G’=GT
これに、各時刻の標準偏差を対角線上におき、他の成分が0である24行24列(24×24)の標準偏差行列Sを右から掛け、平均値ベクトルμをL行並べたL行24列(L×24)の平均値行列Aを加える。これにより、最終的に求める平均と分散および相関を持った乱数行列G”が得られる。
[数5]
G”=G’S+A
この行列G”が予測結果そのものであり、G”のN個の各行ベクトル(1行24列)が、各予測結果に相当する。
なお、時刻分割数が大きい場合や、発電量の相関が強い場合などは、数値計算上の問題として桁落ちによる誤差が大きくなり、相関係数行列Rのコレスキー分解が困難な場合がある。このような場合には、必要に応じて固有値分解や特異値分解を用いる。なお、上記の作成手法において、1日を24時間に分割する代わりに、任意の時間で分割することもできる。つまり、分割の単位は自由である。
[3.具体的適用結果]
次に、本発明(本実施形態)による予測処理を、具体的な事例に適用した結果について説明する。
[3−1.太陽光発電の発電量予測]
図4は、過去の太陽光発電の発電量をプロットした例である。この場合、2時間ごと、7日分のデータが示されている。すなわち、K=12、N=7である。ここでは、類似日選択処理により、すでに同様な気象条件の日が選択されているものとする。同じ気象条件でも、わずかな気温の違いや雲の量の違いから発電量が変動する場合がある。図5は、図4のデータを重ね書きしたものである。
図6は、本発明の手法で生成した発電量シナリオである。この場合、乱数を用いて50通りの予測結果(L=50)が生成されている。但し、予測数は50に限る必要はない。この図6の発電量の予測結果は、図5のデータをよく再現していることが分かる。この図6で重要な点は、元のデータにおける時刻間の相関が、予測結果においても再現されていることである。図7は、本例の場合の相関係数行列である。これは、分散共分散行列を規格化したものである。
太陽光発電機が複数ある場合には、相互がランダムに変動していると、平滑化の効果が見られる。相互の発電量に相関があれば、平滑化の効果は小さくなる。本発明では、このような状況のシミュレーションも可能である。
[3−2.風力発電の発電量予測]
図8は、風力発電の発電量を予測するために、風速と風向を予測する場合の例である。図8において、(a)は過去の風速データの例であり、(b)は対応する風向データの例である。ここでは、図17と同様に、風向が数値化されている。
風向と風速は相互に関連があるために、同時に予測する必要がある。このため、図8(c)のように、対応するデータを形式的に並べて時系列データを作成する。風速、風向ともに24点のデータがあるため、この風力データは48時点からなる時系列となる。図8(c)には、このような時系列が5本重ねてプロットされている。
これらの5本の時系列から、各時点におけるデータの分散、異なる時点間の共分散を計算することで、48×48の分散共分散行列を作成することができる。そして、この分散共分散行列を用いて、48次元の正規分布に従う乱数を生成することで、48時点の長さを持つ時系列を任意の数だけ生成することができる。
図9(a)は、このようにして生成した時系列データである。この時系列の1本、1本が、風速と風向の両方の情報を持っており、各時点において元のデータと同じ分散を持ち、時点間の共分散も元のデータの情報を再現している。すなわち、例えば、13時の風速と18時の風向の相関などが考慮された予測結果になっている。図9(a)の時系列を、中央で分離することで、各時刻の風速と風向をプロットしたものが図9(b)と図9(c)である。
[3−3.複数の発電機の発電量予測]
次に、自然エネルギーによる発電機が複数ある場合の予測処理を、太陽光発電を一例として説明する。複数の太陽光発電機の発電量が、互いに相関がない場合、これらの発電量を独立に予測することができる。しかし、複数の太陽光発電機の発電量は気象データを通じて相互に相関があるのが普通である。複数の太陽光発電機の発電量を合計すると変動が小さくなることが知られている。この場合、相関が小さければ変動が平滑化される可能性が高いが、正の相関が強ければ変動は平滑化されない。
例えば、太陽光発電機がM台の場合の予測を行なう場合、時点数が24個の場合、(M×24)行(M×24)列の成分からなる分散共分散行列を作成し、(M×24)次元の乱数ベクトルを多数発生させることで、各太陽光発電機の発電量を予測することができる。
図10は、2箇所の太陽光発電所における発電量をプロットしたものである。図10(a)の太陽光データ1と図10(b)の太陽光データ2は、過去において1日24時点において、同時に測定されたデータであり、10日分が重ね書かれている。これらのデータは、それぞれ相関を持っている。
そこで、同じ日のデータを形式的に繋いで、48時点の時系列データとしてプロットしたものが図10(c)である。この図10(c)のデータから、前述の方法と同様な方法で、48行48列の分散共分散行列を作成し、乱数を発生させて、生成した時系列データが図11である。図11(a)は計算結果であり、50通りの時系列が表示されている。これを分離したものが図11(b)と図11(c)であり、それぞれ太陽光発電所1と太陽光発電所2の発電量の予測値を示している。
図11では、多数の時系列データが線としてプロットされているが、この場合には各時系列データの実現確率は等しい。従って、図11のある時刻、ある出力の点での時系列データの線の密度がその時刻におけるその出力の実現確率に比例することになる。
図12(a)は、図11(b)の太陽光発電1における12時の発電出力の頻度分布を示したものである。この場合、1000個の時系列データをもとに頻度が計算されている。また、図12(b)は、その確率密度分布を示している。この確率密度分布から予測結果の期待値や発生確率を計算することもできる。例えば、95%の確率で発電出力が4kW〜4.9kWの間になると予測することができる。同様に、図12(c)(d)は、図11(c)の太陽光発電所2における15時の出力の頻度分布と確率密度分布を示している。
図11(b)の発電量と図11(c)の発電量の変動は独立ではなく、太陽光発電所1においてR番目の時系列データが実現した場合には太陽光発電所2においてもR番目の時系列データが実現することを意味している。従って、太陽光発電所1と太陽光発電所2の発電量の合計を求めることは容易であり、対応する時系列データで時刻ごとに和を計算すればよい。
図13は、上記の2箇所の太陽光発電所の発電量の合計を示している。この図から、合計の発電量がどの程度変動するかが見て取れる。この変動量は従来の予測モデルにおける発電量の予測誤差を合計しただけでは得られない。この変動量が実際に系統に与える変動であり、需給制御や負荷周波数制御において考慮すべき変動量である。
負荷周波数制御においては、例えば周波数の変動幅は負荷や他の発電機の発電量が一定であれば、自然エネルギーによる発電量の変動幅に比例する。このときの比例定数は、負荷周波数定数である。従って、図1の予測処理部10に設けられた周波数演算部によって、この図13の変動幅を系統の全発電量で割り、負荷周波数定数の逆数を掛けることで、周波数の変動率を求めることができる。図13の例では太陽光発電機は2台であったが、発電機の数が多い場合にも同様な手法を用いることができる。
[4.実施形態の効果]
以上のような本実施形態によれば、気象データや発電量データの異なる時刻の間の発電量の統計的相関、あるいは位置が異なる発電機の発電量の同時刻の発電量の統計的相関を考慮して、各発電機の発電量を発生確率が付与された複数の時系列データの組として予測することができる。したがって、発電量が不確実に変動する発電機を持つ発電プラントや電力系統であっても、統計確率的な手法を用いることにより、発電量を正確に予測することができる。
また、多数台の発電機の発電量を同時に予測する場合であっても、時刻相関及び位置相関を考慮した予測が可能となる。さらに、発電機、発電プラントが複数ある場合に、それらの発電量を合計(特定のグループ単位でも可)して、全発電量を予測することが可能となる。さらに、多数の発電機の合計の発電量の予測結果から、系統の発電量の変動量を評価し、系統の周波数変動を評価することもできる。
また、風力発電の発電量予測をする場合、風向や風速を予測して発電量を評価することもできるし、発電量を直接に予測することもできる。太陽光発電の発電量の予測をする場合、日照量の予測に用いることもできるが、発電量を直接に予測することもできる。このため、発電量予測の過程で、風向や風速の予測値、日射量の予測値を得ることもできる。
なお、世の中で同時に起こる一連の出来事は一つのシナリオと考えることができる。同一のシナリオに属する太陽光発電シナリオ、風力発電シナリオ、負荷シナリオは、同時に起こる可能性が高いシナリオであり、これらから通常の発電所で発電することが必要な発電量を評価することもできる。また、これにともない、必要なコストや燃料の量なども変動を含めて評価することができる。従って、これらからコストの期待値なども評価可能である。
[5.他の実施形態]
本発明は、上記のような実施形態に限定されるものではない。例えば、本発明で用いられる各種の情報の具体的な内容、値は自由であり、特定の内容、数値には限定されない。請求項の「発電量に関する予測値」、「発電量に関する過去データ」は、気象データ、発電量データ、誤差データ(後述)等、発電量に関する情報を含む広い概念である。
また、上記の例では、多次元の正規分布に従う乱数を用いたが、確率分布は正規分布に限る必要はない。例えば、対数正規分布、指数分布、ワイブル分布、ベータ分布、その他の任意の分布を用いることができる。実際の発電量から分布形状を関数近似により作成することもできる。太陽光発電量のように、理論的に上限(快晴時の日射量による)と下限(=0)が決まっている場合には、例えば、ベータ分布などを用いることでフィッティングの精度を上げることができる。
つまり、気象データ若しくは発電量データを、上記のような分布を用いてモデル化し、分散共分散行列を用いて、かかる分布に従う乱数を多数発生させることにより、発生確率が付与された複数の時系列データの組を生成することにより、予測値を算出するものも、本発明である。
また、本発明は既存の予測モデルと組み合わせることで、精度を向上させることも可能である。ここで、図14は、従来の予測モデルを示した概念図である。太陽光発電に限らず、将来を予測するためのモデルは多数提案されている。これらは全て、現在までに得られている情報から、将来の値を予測するものである。
このために、回帰モデルやニューラルネットなど種々多数のモデルが用いられる。この場合、通常は過去の気象情報と発電量の関係をモデル化する。過去の実績値(実際の気象データ、発電量データ等)を用いる場合もあるが、普通は過去の予測値(気象データ、発電量データ等の予測値)は用いない。また、将来の予測値の誤差や変動を評価するためには、エラーバーや標準偏差などが用いられる。
一方、図15は、本発明を従来の予測方法に適用したモデルの一例である。この場合、典型的には、過去の予測値と過去の実績値のみが用いられる。過去の予測値を得るためには、過去において何らかの予測モデルAが用いられていたわけであるが、この予測モデルの方法に関しては、本発明では、何であってもかまわない(本発明の予測モデルでもよい)。
このような過去の実績値と過去の予測値の差から、過去の誤差データが得られる。この誤差データは、例えば、予測処理部10に設けられた誤差データ演算部により演算することができる。または、入力部30若しくは通信ネットワークを介して外部から入力され、記憶部20に記憶された誤差データを用いることもできる。
ここで、誤差データは、時系列データ生成部12によって、誤差データの時系列の組として生成される。例えば、過去の一日、あるいは一週間、あるいは一ヶ月など、一定周期ごとに計算され、データが存在する日数、あるいは週数、あるいは月数分だけ誤差データの時系列の組が得られる。この誤差データの時系列の組が、K個(K≧2)あるとする。
一方、将来の予測値は、上記の何らかの予測モデルAによって予測されるわけであるが、本発明では、上記の誤差データの時系列の組を用いて、将来の予測結果の誤差を評価することができる。すなわち、前述の手法によって、K個の誤差データの時系列の組から、分散共分散行列作成部131が分散共分散行列を作成し、乱数発生部132が誤差データの予測値に相当する時系列をL組生成する。これは、現在の予測モデルAを用いた場合に期待される誤差の予測値である。
従って、予測処理部10に設けられた補正部が、現在の予測モデルAによる予測結果に、この時系列を加える補正を行うことによって、実現する確率の高い予測結果が得られる。例えば、現在の予測モデルが実績値に比べて小さめの予測結果を出す傾向が強い場合、誤差の時系列データにおいてもマイナスになる傾向が見られる。予測結果の時系列データは、そのような特徴を反映したものになる。
以上のような態様では、過去の気象データあるいは発電量データから、自然エネルギーによる発電機など発電量が不確実な発電機の将来の発電量を誤差を含めて予測することが可能になる。いったん、発電量シナリオが作成されれば、従来と同様の技術により発電計画を作成することも容易であり、作成された発電計画シナリオを用いて、発電量や発電コスト、収益等の期待値を計算したり、最適な発電計画を選択したりすることも可能になる。
なお、かかる場合には、現状の予測モデルは何であってもかまわない。予測モデルの精度は高い方が望ましいが、精度が悪いモデルであっても、本発明を用いることで、精度が改善される。予測モデルは一般に時間が経つにつれて精度が悪くなるものであるが、本発明を用いると同じ予測モデルを長期間使用できるという利点もある。
また、本発明の対象となる時系列データは、太陽光発電、風力発電には限らない。地熱発電や波力発電、潮力発電、その他の自然エネルギーによる発電はもちろん、特定の地域の需要データなどを予測することもできる。もちろん、これらの誤差時系列でもかまわない。さらに、潮力発電など気象データからの予測が困難な場合であっても、過去の発電量の実績データからだけで、発電量の変動量の予測を行うことができる。
1…発電量予測装置
10…予測処理部
11…類似日選択部
12…時系列データ生成部
13…予測値算出部
14…発電量合計部
20…記憶部
30…入力部
40…出力部
131…分散共分散行列作成部
132…乱数発生部
133…発電量計算部

Claims (9)

  1. 発電機の発電量に関する予測値を求める発電量予測装置において、
    少なくとも一つの発電機の過去の発電量に関する過去データを記憶する過去データ記憶部と、
    前記過去データにおける異なる時刻の間の統計的相関若しくは異なる発電機の位置の間の統計的相関に基づいて、発電機の発電量に関する予測値を、発生確率を含む時系列データとして算出する予測値算出部と、
    を有し、
    前記予測値算出部は、
    前記過去データに基づいて、分散共分散行列を作成する分散共分散行列作成部と、
    前記分散共分散行列に基づいて、確率分布に従う乱数を発生させる乱数発生部と、
    を有することを特徴とする発電量予測装置。
  2. 前記過去データは、発電量データであることを特徴とする請求項1記載の発電量予測装置。
  3. 前記過去データは、発電量に影響を与える気象データであり、
    前記予測値算出部は、気象データの予測値を、発電量に換算する発電量計算部を有することを特徴とする請求項1記載の発電量予測装置。
  4. 複数の発電機における発電量の予測値に基づいて、複数を合計した発電量の予測値を算出する発電量合計部を有することを特徴とする請求項1記載の発電量予測装置。
  5. 前記合計した発電量の予測値に基づいて、周波数の変動率を演算する周波数演算部を有することを特徴とする請求項記載の発電量予測装置。
  6. 前記過去データは、予測値の誤差データであり、
    前記予測値は、誤差データの予測値であることを特徴とする請求項1記載の発電量予測装置。
  7. 気象データの予測値若しくは発電量データの予測値を、前記誤差データの予測値によって補正する補正部を有することを特徴とする請求項記載の発電量予測装置。
  8. コンピュータにより発電量を予測する発電量予測方法において、
    前記コンピュータは、過去データ記憶部と予測値算出部とを有し、
    前記過去データ記憶部が、少なくとも一つの発電機の過去の発電量に関する過去データを記憶する処理と、
    前記予測値算出部が、前記過去データにおける異なる時刻の間の統計的相関若しくは異なる発電機の位置の間の統計的相関に基づいて、発電機の発電量に関する予測値を、発生確率を含む時系列データとして算出する予測値算出処理と、
    を実行し、
    前記予測値算出部は、分散共分散行列作成部と乱数発生部とを有し、
    前記予測値算出処理を実行するにあたり、
    前記分散共分散行列作成部は、前記過去データに基づいて、分散共分散行列を作成する処理と、
    前記乱数発生部は、前記分散共分散行列に基づいて、確率分布に従う乱数を発生させる処理と、
    を実行することを特徴とする発電量予測方法。
  9. コンピュータにより実行されることにより、発電量を予測する発電量予測プログラムにおいて、
    前記コンピュータに、
    少なくとも一つの発電機の過去の発電量に関する過去データを記憶する処理と、
    前記過去データにおける異なる時刻の間の統計的相関若しくは異なる発電機の位置の間の統計的相関に基づいて、発電機の発電量に関する予測値を、発生確率を含む時系列データとして算出する予測値算出処理と、
    を実行させ、
    前記コンピュータに前記予測値算出処理を実行させるにあたり、
    前記過去データに基づいて、分散共分散行列を作成する処理と、
    前記分散共分散行列に基づいて、確率分布に従う乱数を発生させる処理と、
    を実行させることを特徴とする発電量予測プログラム。
JP2010065083A 2010-03-19 2010-03-19 発電量予測装置、予測方法及び予測プログラム Active JP5606114B2 (ja)

Priority Applications (7)

Application Number Priority Date Filing Date Title
JP2010065083A JP5606114B2 (ja) 2010-03-19 2010-03-19 発電量予測装置、予測方法及び予測プログラム
TW100109378A TWI533251B (zh) 2010-03-19 2011-03-18 發電量預測裝置及其方法以及記錄媒體
PCT/JP2011/001648 WO2011114752A1 (ja) 2010-03-19 2011-03-18 発電量予測装置とその方法及びプログラム
CN201180013084.7A CN102792542B (zh) 2010-03-19 2011-03-18 发电量预测装置和其方法及程序
MYPI2012700591A MY156889A (en) 2010-03-19 2011-03-18 Electric-power-generation level predicting apparatus, method and program
EP11755942.7A EP2549641A4 (en) 2010-03-19 2011-03-18 DEVICE FOR FORECASTING THE AMOUNT OF A GENERATED CURRENT, METHOD THEREFOR AND PROGRAM THEREFOR
US13/583,792 US20130013233A1 (en) 2010-03-19 2011-03-18 Electric-power-generation level predicting apparatus, method and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010065083A JP5606114B2 (ja) 2010-03-19 2010-03-19 発電量予測装置、予測方法及び予測プログラム

Publications (2)

Publication Number Publication Date
JP2011200040A JP2011200040A (ja) 2011-10-06
JP5606114B2 true JP5606114B2 (ja) 2014-10-15

Family

ID=44648862

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010065083A Active JP5606114B2 (ja) 2010-03-19 2010-03-19 発電量予測装置、予測方法及び予測プログラム

Country Status (7)

Country Link
US (1) US20130013233A1 (ja)
EP (1) EP2549641A4 (ja)
JP (1) JP5606114B2 (ja)
CN (1) CN102792542B (ja)
MY (1) MY156889A (ja)
TW (1) TWI533251B (ja)
WO (1) WO2011114752A1 (ja)

Families Citing this family (55)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5505191B2 (ja) * 2010-08-16 2014-05-28 東京電力株式会社 太陽光発電量予測方法および配電系統制御システム
JP5661594B2 (ja) * 2011-11-01 2015-01-28 日本電信電話株式会社 予測モデル構築装置、方法、及びプログラム、並びに発電量予測装置、及び方法
JP6178045B2 (ja) 2012-04-16 2017-08-09 株式会社東芝 エネルギー管理システム、エネルギー管理方法、プログラムおよびサーバ装置
JP5797599B2 (ja) * 2012-04-19 2015-10-21 株式会社日立パワーソリューションズ 発電量予測方法及びそのシステム
TWI509942B (zh) * 2012-06-21 2015-11-21 Au Optronics Corp 多電源供電系統之電力管理方法
JP5940922B2 (ja) * 2012-07-12 2016-06-29 トヨタ自動車株式会社 自然エネルギー量予測装置
JP2014054048A (ja) * 2012-09-06 2014-03-20 Toshiba Corp 電力予測方法、装置、及びプログラム
JP6054122B2 (ja) 2012-09-28 2016-12-27 株式会社東芝 鉄道電力管理装置
JP5946742B2 (ja) * 2012-10-16 2016-07-06 一般財団法人電力中央研究所 自然エネルギー型分散電源群の合計出力の変動推定方法、変動推定装置及び変動推定プログラム
TWI492158B (zh) * 2012-11-15 2015-07-11 Inventec Corp 負載預測方法與電子裝置
JP6009976B2 (ja) 2013-03-07 2016-10-19 株式会社東芝 エネルギー管理システム、エネルギー管理方法、プログラムおよびサーバ
JP2014217092A (ja) * 2013-04-22 2014-11-17 清水建設株式会社 発電電力予測装置、発電電力予測方法及びプログラム
US20150100386A1 (en) * 2013-10-04 2015-04-09 International Business Machines Corporation Digital framework for business model innovation
KR101575071B1 (ko) 2013-12-02 2015-12-07 두산중공업 주식회사 풍력 발전 단지의 발전량 제어 방법
JP6181589B2 (ja) * 2014-03-31 2017-08-16 株式会社Ihi 計測値解析装置及び計測値解析方法
KR101585961B1 (ko) 2014-04-23 2016-02-02 전자부품연구원 풍력-파력 복합 발전에서의 상태 감시 및 예측 시스템
JP6303909B2 (ja) * 2014-08-11 2018-04-04 富士通株式会社 計画策定方法、計画策定システム及び計画策定プログラム
JP2016053925A (ja) * 2014-09-04 2016-04-14 株式会社東芝 電力市場価格予測装置、電力市場価格予測方法およびプログラム
JP6400403B2 (ja) * 2014-09-12 2018-10-03 株式会社東芝 契約電力選択装置、契約電力選択方法、及び、契約電力選択プログラム
JP5989732B2 (ja) * 2014-10-01 2016-09-07 中国電力株式会社 発電量予測装置、発電量予測方法、及びプログラム
JP5937167B2 (ja) * 2014-10-01 2016-06-22 中国電力株式会社 需要予測装置、需要予測方法、及びプログラム
JP6543028B2 (ja) * 2014-11-05 2019-07-10 東芝エネルギーシステムズ株式会社 異常診断装置及び異常診断方法
JP6582418B2 (ja) * 2015-01-23 2019-10-02 中国電力株式会社 発電量予測装置、発電量予測方法、及びプログラム
JPWO2016121202A1 (ja) * 2015-01-30 2017-11-09 日本電気株式会社 予測装置、予測方法、及び、プログラム
CN107251356B (zh) * 2015-02-23 2020-09-11 学校法人冲绳科学技术大学院大学学园 确定可再生能源波动的预测误差的系统和方法
JP6600516B2 (ja) * 2015-09-14 2019-10-30 株式会社東芝 アグリゲーション管理装置及び方法
US9773303B2 (en) * 2015-11-04 2017-09-26 Institute of Nuclear Energy Research, Atomic Energy Council, Executive Yuan, R.O.C. Prediction system for short-term solar irradiance and method for operating the same
TWI642019B (zh) * 2015-11-17 2018-11-21 財團法人資訊工業策進會 發電系統效能估測方法
JP6427090B2 (ja) 2015-12-07 2018-11-21 株式会社日立製作所 発電量予測装置、発電量予測方法、系統安定化装置、並びに系統安定化方法
CN105515042A (zh) * 2015-12-08 2016-04-20 上海电力设计院有限公司 风光同场项目光伏接入容量优化算法
US10509868B2 (en) * 2016-01-29 2019-12-17 Fujitsu Limited Operating a solar power generating system
JP6702408B2 (ja) * 2016-03-08 2020-06-03 日本電気株式会社 電力制御装置、電力制御システム、電力制御方法、及び、プログラム
US10824119B2 (en) * 2016-03-29 2020-11-03 International Business Machines Corporation Intelligent energy switch
JP6742894B2 (ja) * 2016-06-09 2020-08-19 株式会社日立製作所 データ予測システムおよびデータ予測方法
WO2017212880A1 (ja) * 2016-06-09 2017-12-14 株式会社日立製作所 データ予測システムおよびデータ予測方法
JP6725406B2 (ja) * 2016-12-12 2020-07-15 株式会社東芝 発電計画策定装置、発電計画策定方法、および発電計画策定プログラム
US10418822B2 (en) * 2017-02-14 2019-09-17 Mitsubishi Electric Research Laboratories, Inc. Energy production and frequency regulation Co-optimization for power generation systems
TWI623890B (zh) * 2017-07-04 2018-05-11 台灣電力股份有限公司 由多個類神經網路組成的發電量預測系統與其方法
KR102395474B1 (ko) 2017-08-24 2022-05-09 삼성전자주식회사 반도체 소자의 특성 예측 방법 및 반도체 소자의 특성 예측 장치
TWI632516B (zh) * 2018-01-11 2018-08-11 國立臺灣師範大學 日射量比對系統與方法
CN110676868A (zh) * 2018-07-02 2020-01-10 赫普科技发展(北京)有限公司 一种风光气互补耦合发电系统及方法
CN109165777A (zh) * 2018-08-10 2019-01-08 广东电网有限责任公司 一种居民用电电量预测方法、装置、设备及存储介质
KR102534367B1 (ko) * 2018-10-08 2023-05-22 한국전력공사 전력 시계열 데이터의 전송처리 시스템 및 방법
CN109412200B (zh) * 2018-10-24 2021-12-10 安徽科达售电有限公司 一种可再生能源发电的接入控制系统
CN110008981B (zh) * 2019-01-11 2023-08-08 国网浙江省电力有限公司 基于矢量模糊矩阵模型的电水气用能替代可引导评价方法
CN111985673A (zh) * 2019-05-22 2020-11-24 三菱重工业株式会社 预测装置、预测系统、预测方法以及记录介质
US11261846B2 (en) 2019-11-01 2022-03-01 General Electric Company System and method for designing and operating a wind turbine power system based on statistical analysis of operational and/or grid data thereof
JP7436980B2 (ja) * 2020-01-22 2024-02-22 日本エア・リキード合同会社 液化装置
US11636390B2 (en) 2020-03-19 2023-04-25 International Business Machines Corporation Generating quantitatively assessed synthetic training data
JP7443173B2 (ja) 2020-07-08 2024-03-05 株式会社東芝 電力系統状態推定方法、電力系統状態推定プログラム、電力系統状態推定プログラムを記憶した記憶媒体、及び電力系統状態推定装置
CN112186761B (zh) * 2020-09-30 2022-03-01 山东大学 一种基于概率分布的风电功率场景生成方法及系统
JP2023023114A (ja) * 2021-08-04 2023-02-16 トヨタ自動車株式会社 電力管理システム、電力管理サーバ、および、電力管理方法
FR3135798A1 (fr) * 2022-05-19 2023-11-24 IFP Energies Nouvelles Procédé de prévision d’une puissance produite par au moins un panneau photovoltaïque
CN116123028A (zh) * 2022-12-20 2023-05-16 中国华能集团清洁能源技术研究院有限公司 风电场级mppt预测模型控制方法和装置
KR102642564B1 (ko) * 2023-06-30 2024-02-29 식스티헤르츠 주식회사 재생 에너지 발전기 이상 진단 방법 및 시스템

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3142489B2 (ja) * 1996-08-26 2001-03-07 沖電気工業株式会社 目標運動解析方法及び目標運動解析装置
US20030189420A1 (en) * 2001-09-28 2003-10-09 Hiroyuki Hashimoto Method for assisting in planning of power supply schedule
JP3994910B2 (ja) * 2003-05-08 2007-10-24 株式会社日立製作所 電力売買支援システム
US6966073B2 (en) * 2003-12-01 2005-11-22 Timothy Lee Sock
EP1744058A1 (en) * 2004-05-07 2007-01-17 Mitsubishi Denki Kabushiki Kaisha Wind power generation evaluation system and prediction control service system for wind power generator
JP4245583B2 (ja) * 2005-04-15 2009-03-25 日本電信電話株式会社 分散型エネルギーシステムの制御装置、制御方法、プログラム、および記録媒体
JP2007056686A (ja) * 2005-08-22 2007-03-08 Univ Of Ryukyus 風速予測に基づく風力発電機の数時間先発電電力予測装置
JP2007120362A (ja) * 2005-10-26 2007-05-17 Toshiba Corp 発電電力予測装置
US7716100B2 (en) * 2005-12-02 2010-05-11 Kuberre Systems, Inc. Methods and systems for computing platform
JP4886322B2 (ja) * 2006-02-28 2012-02-29 東北電力株式会社 風力発電出力予測方法、風力発電出力予測装置およびプログラム
JP2007281060A (ja) * 2006-04-04 2007-10-25 Nippon Telegr & Teleph Corp <Ntt> 太陽光発電システムの発電予測方法、装置、およびプログラム
US8401708B2 (en) * 2007-03-26 2013-03-19 Vpec, Inc. Electric power system
JP2009138523A (ja) * 2007-12-03 2009-06-25 Mitsubishi Electric Corp 風力発電出力予測方法
JP5020106B2 (ja) * 2008-01-11 2012-09-05 中国電力株式会社 風力発電量予測システム、方法及びプログラム
JP5190311B2 (ja) * 2008-07-09 2013-04-24 株式会社東芝 取引対象物に関する将来需要の予測シナリオ作成方法およびその作成装置
KR101075958B1 (ko) * 2009-11-16 2011-10-21 일진전기 주식회사 분산전원의 전력 거래 시스템 및 거래 방법

Also Published As

Publication number Publication date
JP2011200040A (ja) 2011-10-06
EP2549641A4 (en) 2016-11-30
MY156889A (en) 2016-04-15
CN102792542B (zh) 2015-01-07
US20130013233A1 (en) 2013-01-10
EP2549641A1 (en) 2013-01-23
CN102792542A (zh) 2012-11-21
WO2011114752A1 (ja) 2011-09-22
TWI533251B (zh) 2016-05-11
TW201205492A (en) 2012-02-01

Similar Documents

Publication Publication Date Title
JP5606114B2 (ja) 発電量予測装置、予測方法及び予測プログラム
Ding et al. A novel adaptive discrete grey model with time-varying parameters for long-term photovoltaic power generation forecasting
Behera et al. Solar photovoltaic power forecasting using optimized modified extreme learning machine technique
Zakaria et al. Uncertainty models for stochastic optimization in renewable energy applications
Mohan et al. A data-driven strategy for short-term electric load forecasting using dynamic mode decomposition model
Persson et al. Multi-site solar power forecasting using gradient boosted regression trees
Almeida et al. PV power forecast using a nonparametric PV model
Pousinho et al. A risk-averse optimization model for trading wind energy in a market environment under uncertainty
Catalão et al. Mixed-integer nonlinear approach for the optimal scheduling of a head-dependent hydro chain
Shayesteh et al. Multi-station equivalents for short-term hydropower scheduling
Muzumdar et al. Designing a robust and accurate model for consumer-centric short-term load forecasting in microgrid environment
Schinke-Nendza et al. Probabilistic forecasting of photovoltaic power supply—A hybrid approach using D-vine copulas to model spatial dependencies
Inage Development of an advection model for solar forecasting based on ground data first report: Development and verification of a fundamental model
Sharma et al. Short-term wind speed forecasting: Application of linear and non-linear time series models
Bashari et al. Forecasting electric load by aggregating meteorological and history-based deep learning modules
WO2022168357A1 (ja) 発電量管理システム及び発電量管理方法
Bødal et al. Capacity expansion planning with stochastic rolling horizon dispatch
Matsuo et al. Forecasting daily electric load by applying artificial neural network with fourier transformation and principal component analysis technique
Memarzadeh et al. Solar power generation forecasting by a new hybrid cascaded extreme learning method with maximum relevance interaction gain feature selection
Zhang et al. Long-term energy and peak power demand forecasting based on sequential-XGBoost
JP7346356B2 (ja) 予測装置および予測方法
Abubakar et al. Intelligent modeling and optimization of solar plant production integration in the smart grid using machine learning models
Pramanik et al. An ensemble-based approach for short-term load forecasting for buildings with high proportion of renewable energy sources
CN116316617A (zh) 多场站智能融合的新能源发电功率区域预测方法和系统
Kronqvist et al. Alternating mixed-integer programming and neural network training for approximating stochastic two-stage problems

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120808

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131203

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140131

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140729

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140826

R151 Written notification of patent or utility model registration

Ref document number: 5606114

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151