JP2018036196A - Wind speed prediction device, wind speed prediction system, wind speed predication method and program - Google Patents
Wind speed prediction device, wind speed prediction system, wind speed predication method and program Download PDFInfo
- Publication number
- JP2018036196A JP2018036196A JP2016170918A JP2016170918A JP2018036196A JP 2018036196 A JP2018036196 A JP 2018036196A JP 2016170918 A JP2016170918 A JP 2016170918A JP 2016170918 A JP2016170918 A JP 2016170918A JP 2018036196 A JP2018036196 A JP 2018036196A
- Authority
- JP
- Japan
- Prior art keywords
- wind speed
- prediction
- learning data
- evaluation index
- time
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims description 34
- 230000010365 information processing Effects 0.000 claims description 6
- 108090000672 Annexin A5 Proteins 0.000 claims 1
- 238000012545 processing Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 238000010248 power generation Methods 0.000 description 10
- 238000010276 construction Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
本発明は、風速予測装置、風速予測システム、風速予測方法及びプログラムに関する。 The present invention relates to a wind speed prediction apparatus, a wind speed prediction system, a wind speed prediction method, and a program.
計算機等が、風力発電において発電される電力を気象条件等に基づいて計算して予測する方法が知られている。 A method is known in which a computer or the like calculates and predicts electric power generated in wind power generation based on weather conditions or the like.
例えば、風力発電出力を精度良く予測する方法が知られている。具体的には、まず、風力発電出力予測装置が、過去の地形条件、気象条件及び実測値を用いて、複数の物理モデル及び複数の統計モデルについて集団学習し、複数の物理モデル及び複数の統計モデルを最適に組み合わせた最終モデルを決定する。次に、風力発電出力予測装置が、最終モデルを用いて、過去よりも後の特定の時刻における地形条件、気象条件及び実測値に基づいて、特定の時刻における対象地点の風力発電出力の最終予測値を算出する。このように、最適化された最終モデルを用いることで、風力発電出力予測装置が風力発電出力を精度良く予測できる方法が知られている(例えば、特許文献1)。 For example, a method for accurately predicting wind power generation output is known. Specifically, first, the wind power generation output prediction device performs collective learning on a plurality of physical models and a plurality of statistical models using past terrain conditions, weather conditions, and actually measured values, and then a plurality of physical models and a plurality of statistics. Determine the final model that best combines the models. Next, the wind power generation output prediction device uses the final model, and finally predicts the wind power generation output at the target point at a specific time based on the topographic conditions, weather conditions, and actual measurement values at a specific time after the past. Calculate the value. As described above, there is known a method in which the wind power generation output prediction apparatus can accurately predict the wind power generation output by using the optimized final model (for example, Patent Document 1).
他にも、計算機等が、風速の3乗に基づいて、発電される電力を計算して予測する方法等が知られている。 In addition, a method is known in which a computer or the like calculates and predicts the generated power based on the cube of the wind speed.
また、自然エネルギー量予測装置が、まず、自然エネルギー量の変化に関連する説明変数を取得して、確率分布を計算する。そして、自然エネルギー量予測装置が、自然エネルギー量の予測値の信頼区間を確率分布に基づいて計算する。このようにして、自然エネルギーから得られる単位時間当たりのエネルギー量の予測値を予測誤差の分布を含めて精度良く予測する方法が知られている(例えば、特許文献2)。 In addition, the natural energy amount prediction apparatus first acquires an explanatory variable related to a change in the natural energy amount and calculates a probability distribution. Then, the natural energy amount prediction device calculates the confidence interval of the predicted value of the natural energy amount based on the probability distribution. Thus, a method for accurately predicting the predicted value of the amount of energy per unit time obtained from natural energy including the distribution of prediction errors is known (for example, Patent Document 2).
しかしながら、従来の方法では、季節等による風速の変動特性等が考慮されていない。すなわち、風速は、季節等によって変動特性があるため、変動特性を考慮しないで風速の予測が行われると、予測誤差の信頼区間を精度良く予測できない場合がある。 However, the conventional method does not consider the fluctuation characteristics of the wind speed due to the season or the like. That is, since the wind speed has a variation characteristic depending on the season and the like, if the wind speed is predicted without considering the variation characteristic, the confidence interval of the prediction error may not be accurately predicted.
本発明の1つの側面は、このような問題に鑑みてなされたものであり、風速の予測において、予測誤差の信頼区間を精度良く予測することを目的とする。 One aspect of the present invention has been made in view of such problems, and an object thereof is to accurately predict a confidence interval of a prediction error in wind speed prediction.
上述した課題を解決し、目的を達成するため、本発明の一実施形態における、風速を予測する風速予測装置は、
現時点より過去の時点における風速を計測した実績値及び前記実績値を予測した予測値を示す学習データを入力する入力部と、
前記学習データに基づいて、前記現時点から所定の時間後となる予測時点における風速のモデルと、前記モデルに基づいて前記予測時点における予測誤差の信頼区間とを少なくとも計算する計算部と、
前記予測誤差の信頼区間の評価指標を計算する評価指標計算部と、
前記評価指標に基づいて、前記学習データのうち、予測に用いる期間データを決定する決定部と
を含む。
In order to solve the above-described problems and achieve the object, the wind speed prediction apparatus for predicting the wind speed in one embodiment of the present invention is
An input unit for inputting learning data indicating an actual value obtained by measuring the wind speed at a time point in the past from the present time and a predicted value obtained by predicting the actual value;
Based on the learning data, a calculation unit for calculating at least a wind speed model at a prediction time after a predetermined time from the current time, and a prediction error confidence interval at the prediction time based on the model;
An evaluation index calculation unit for calculating an evaluation index of the confidence interval of the prediction error;
And a determination unit that determines period data to be used for prediction among the learning data based on the evaluation index.
本発明によれば、風速の予測において、予測誤差の信頼区間を精度良く予測することができる。 According to the present invention, it is possible to accurately predict a confidence interval of a prediction error in wind speed prediction.
以下、本発明の実施形態を図面に基づいて説明する。
1.風速予測装置のハードウェア構成例
2.風速予測装置による全体処理例
3.風速予測装置の機能構成例
≪ 1. 風速予測装置のハードウェア構成例 ≫
風速予測装置は、例えば、PC(Personal Computer)等の情報処理装置である。以下、風速予測装置がPCである例で説明する。
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
1. 1. Hardware configuration example of wind
≪ 1. Example of hardware configuration of wind speed prediction device ≫
The wind speed prediction apparatus is an information processing apparatus such as a PC (Personal Computer). Hereinafter, an example in which the wind speed prediction apparatus is a PC will be described.
図1は、本発明の一実施形態における風速予測装置のハードウェア構成の一例を示すブロック図である。具体的には、PC10は、CPU(Central Processing Unit)101と、記憶装置102と、ネットワークI/F(interface)103と、入力I/F104と、出力I/F105とを有する。つまり、風速予測装置は、プログラムを実行して処理を行うコンピュータである。
FIG. 1 is a block diagram illustrating an example of a hardware configuration of a wind speed prediction apparatus according to an embodiment of the present invention. Specifically, the PC 10 includes a CPU (Central Processing Unit) 101, a
CPU101は、各種処理及び各種制御を実現するための演算と各種データの加工とを行う演算装置である。さらに、CPU101は、PC10が有するハードウェアを制御する制御装置である。
The
記憶装置102は、PC10が使うデータ、プログラム及び設定値等を記憶する。また、記憶装置102は、いわゆるメモリ(memory)等である。なお、記憶装置102は、ハードディスク(harddisk)等の補助記憶装置等を有してもよい。
The
ネットワークI/F103は、LAN(Local Area Network)等のネットワークを介して接続される装置と各種データ等を送受信する。例えば、ネットワークI/F103は、NIC(Network Interface Controller)及びLANケーブルを接続させるコネクタ等である。なお、ネットワークI/F103は、ネットワークを利用するI/Fに限られず、ケーブル、無線又はコネクタ等によって外部装置と送受信するI/Fであってもよい。
The network I /
入力I/F104は、PC10を使うユーザとのインタフェースである。具体的には、入力I/F104は、ユーザが行う各種操作を入力する。例えば、入力I/F104は、キーボード等の入力装置及び入力装置をPC10に接続させるコネクタ等によって構成される。
The input I /
出力I/F105は、PC10を使うユーザとのインタフェースである。具体的には、出力I/F105は、PC10が行う各種処理の処理結果等をユーザに出力する。例えば、出力I/F105は、ディスプレイ等の出力装置及び出力装置をPC10に接続させるコネクタ等である。
The output I / F 105 is an interface with a user who uses the PC 10. Specifically, the output I /
なお、PC10は、各ハードウェア資源による処理等を補助する補助装置を更に有する構成でもよい。また、PC10は、各種処理を並列、冗長又は分散して処理するため、装置を内部又は外部に更に有してもよい。さらに、PC10は、複数の情報処理装置で構成されてもよい。 Note that the PC 10 may further include an auxiliary device that assists the processing by each hardware resource. Further, the PC 10 may further include an apparatus inside or outside in order to perform various processes in parallel, redundantly, or distributedly. Further, the PC 10 may be composed of a plurality of information processing apparatuses.
≪ 2. 風速予測装置による全体処理例 ≫
図2は、本発明の一実施形態における風速予測装置による全体処理の一例を示すフローチャートである。例えば、PCは、図示するような処理を行う。
≪ 2. Example of overall processing by wind speed prediction device ≫
FIG. 2 is a flowchart showing an example of overall processing by the wind speed prediction apparatus according to the embodiment of the present invention. For example, the PC performs processing as illustrated.
≪ 学習データの入力例(ステップS101) ≫
ステップS101では、PCは、学習データDLを入力する。学習データDLには、風速の予測を行う現時点より過去の時点において、風速を計測した結果(以下「実績値」という。)が含まれる。また、学習データDLには、実績値を予測した結果(以下「予測値」という。)が含まれる。例えば、学習データDL、すなわち、実績値及び予測値は、以下のようなデータである。
≪ Example of learning data input (step S101) ≫
In step S101, the PC inputs learning data DL. The learning data DL includes a result of measuring the wind speed (hereinafter referred to as “actual value”) at a time point in the past from the present time when the wind speed is predicted. Further, the learning data DL includes a result of predicting the actual value (hereinafter referred to as “predicted value”). For example, the learning data DL, that is, the actual value and the predicted value are the following data.
図3は、本発明の一実施形態における学習データの一例を示す図である。図では、横軸を時間とし、縦軸を風速とする。また、この例では、予測を行う現在の時点を現時点TCとする。 FIG. 3 is a diagram illustrating an example of learning data according to an embodiment of the present invention. In the figure, the horizontal axis is time, and the vertical axis is wind speed. In this example, the current time point when the prediction is performed is a current time TC.
図示する学習データDLは、現時点TCより過去に、6時点において実績値が計測された場合の例のデータである。具体的には、図示する例では、第1時間T1、第2時間T2、第3時間T3、第4時間T4、第5時間T5及び第6時間T6の時点で、それぞれの実績値がそれぞれ計測されたとする。図では、各時点で計測されたそれぞれの実績値RESを点で示す。なお、実績値RESは、不定期に計測された値でもよい。 The learning data DL shown in the figure is data of an example in the case where the actual value is measured at six time points in the past from the current time point TC. Specifically, in the illustrated example, each actual value is measured at the time of the first time T1, the second time T2, the third time T3, the fourth time T4, the fifth time T5, and the sixth time T6. Suppose that In the figure, each actual value RES measured at each time point is indicated by a point. The actual value RES may be a value measured irregularly.
そして、図示するように、学習データDLには、実績値RESを予測した予測値CALが入力される。すなわち、各時点より更に過去の時点において、各時点の風速を予測した結果が予測値CALである。このように、学習データDLは、複数の時点において計測された実績値と、実績値を予測した予測値とがあらかじめ入力されたデータである。例えば、PCは、1年分の学習データDLを入力する。以下、学習データDLが1年分のデータである例で説明する。 As shown in the figure, a predicted value CAL obtained by predicting the actual value RES is input to the learning data DL. That is, the predicted value CAL is the result of predicting the wind speed at each time point in the past. As described above, the learning data DL is data in which the actual value measured at a plurality of time points and the predicted value obtained by predicting the actual value are input in advance. For example, the PC inputs learning data DL for one year. Hereinafter, an example in which the learning data DL is data for one year will be described.
≪ 学習データの構築例(ステップS102) ≫
図2に戻り、ステップS102では、PCは、学習データを構築する。具体的には、まず、PCは、学習データを分割する。なお、PCは、分割されたそれぞれのデータ(以下「期間データ」という。)が、均等に、所定の期間となるように学習データを分割する。例えば、所定の期間が「1か月」とあらかじめ設定される場合には、1年分の学習データは、1か月単位で期間データに分割され、各期間データは、1か月分の各時点における実績値及び予測値を示すデータとなる。そして、このように1か月単位で分割されると、1年分の学習データから、12個の期間データが分割によって生成される。以下、学習データは、「m」個の期間データに分割されるとする。なお、所定の期間は、例えば、1分単位、1時間単位、1日単位、1週間単位又は季節単位等でもよい。そして、所定の期間は、例えば、ユーザ操作等によってPCにあらかじめ設定される。
<< Example of construction of learning data (step S102) >>
Returning to FIG. 2, in step S102, the PC constructs learning data. Specifically, first, the PC divides the learning data. The PC divides the learning data so that each of the divided data (hereinafter referred to as “period data”) is evenly within a predetermined period. For example, when the predetermined period is set as “1 month” in advance, the learning data for one year is divided into period data in units of one month, and each period data is The data shows the actual value and the predicted value at the time. When divided in units of one month in this way, twelve period data are generated by division from learning data for one year. Hereinafter, it is assumed that the learning data is divided into “m” period data. The predetermined period may be, for example, 1 minute unit, 1 hour unit, 1 day unit, 1 week unit, or season unit. The predetermined period is preset in the PC by a user operation or the like, for example.
また、期間データを以下の説明では、下記(1)式のように示す。 In the following description, the period data is represented by the following formula (1).
次に、PCは、分割によって生成された複数の期間データのうち、1つ又は2つ以上の組み合わせを選択する。以下、現時点から近い順に、複数の期間データが、PCによって選択されるとする。 Next, the PC selects one or a combination of two or more of the plurality of period data generated by the division. Hereinafter, it is assumed that a plurality of pieces of period data are selected by the PC in the order from the current time point.
≪ 学習データに基づくモデルの計算例(ステップS103) ≫
図2に戻り、ステップS103では、PCは、学習データに基づいてモデルを計算する。すなわち、PCは、ステップS102によって分割された期間データのうち、1つ又は複数の期間データを用いて、モデルを計算する。
<< Example of model calculation based on learning data (step S103) >>
Returning to FIG. 2, in step S103, the PC calculates a model based on the learning data. That is, the PC calculates a model using one or a plurality of period data among the period data divided in step S102.
以下、予測を行う現時点を「i」とする。そして、PCは、「i」の現時点から、所定の時間後となる予測時点の風速を予測する。また、以下の説明では、所定の時間単位を「ステップ」という。また、予測時点は、「kステップ」先の時点という。具体的には、所定の時間を「1時間」単位とすると、「1ステップ」は、「1時間」となる。したがって、この例では、例えば、5時間後の予測時点は、「5ステップ」先の時点となる。 Hereinafter, the current time point when the prediction is performed is “i”. Then, the PC predicts the wind speed at the predicted time point after a predetermined time from the current time point of “i”. In the following description, the predetermined time unit is referred to as “step”. The predicted time point is referred to as “k step” ahead. Specifically, assuming that the predetermined time is in units of “1 hour”, “1 step” becomes “1 hour”. Therefore, in this example, for example, the predicted time point after 5 hours is the time point “5 steps” ahead.
そして、過去の「i」時点において、「kステップ」先の予測時点における予測値及び実績値に基づいて、予測時点における予測誤差は、下記(2)式のように計算される。 Then, at the past “i” time point, the prediction error at the prediction time point is calculated as in the following equation (2) based on the prediction value and the actual value at the prediction time point “k steps” ahead.
≪ 評価指標の計算例(ステップS104) ≫
図2に戻り、ステップS104では、PCは、評価指標をモデルごとに計算する。例えば、評価指標は、PICP(Prediction Interval Coverage Probability)である。また、PICPは、例えば、下記(5)式によって計算される値である。
<< Calculation example of evaluation index (step S104) >>
Returning to FIG. 2, in step S104, the PC calculates an evaluation index for each model. For example, the evaluation index is PICP (Prediction Interval Coverage Probability). Also, PICP is a value calculated by the following equation (5), for example.
PICPのより詳しい計算方法は、例えば、「Abbas Khosravi, Saeid Nahavandi and Doug Creighton,"Construction of Optimal Prediction Intervals for Load Forecasting Probrem",IEEE Transactions on Power Systems, Vol.25, Issue 3, pp.1496―1503(2010)」に記載されている方法等である。
More detailed calculation method of PICP is, for example, "Abbas Khosravi, Saeid Nahavandi and Doug Creighton," Construction of Optimal Prediction Intervals for Load Forecasting Probrem ", IEEE Transactions on Power Systems, Vol.25,
これらの推定値は、上記(5)式における「lj」及び「uj」となる。一方で、「1σ」となる標準偏差では、予測誤差の信頼区間に理論上、実績値が含まれる確率は、約68.27パーセントとなる。さらに、PICPの理論値は、標準偏差の整数倍ごとに、下記(6)式のようになる。 These estimated values are “l j ” and “u j ” in the above equation (5). On the other hand, with the standard deviation of “1σ”, the probability that the actual value is theoretically included in the confidence interval of the prediction error is about 68.27%. Further, the theoretical value of PICP is expressed by the following equation (6) for each integer multiple of the standard deviation.
≪ 「k=kmax」であるか否かの判断例(ステップS105) ≫
図2に戻り、ステップS105では、PCは、「k=kmax」であるか否かを判断する。そして、「k=kmax」でないとPCが判断すると(ステップS105でNO)、PCは、ステップS106に進む。一方で、「k=kmax」であるとPCが判断すると(ステップS105でYES)、PCは、ステップS107に進む。
<< Judgment example of whether or not “k = k max ” (step S105) >>
Returning to FIG. 2, in step S <b> 105, the PC determines whether or not “k = k max ”. If the PC determines that “k = k max ” is not satisfied (NO in step S105), the PC proceeds to step S106. On the other hand, when the PC determines that “k = k max ” (YES in step S105), the PC proceeds to step S107.
「kmax」は、予測を行うステップの最大値である。具体的には、「80」ステップ先まで予測を行う場合には、「kmax=80」となる。 “K max ” is the maximum value of the step for performing prediction. Specifically, when prediction is performed up to “80” steps ahead, “k max = 80”.
≪ 「k=k+1」の計算例(ステップS106) ≫
ステップS106では、PCは、「k=k+1」を計算する。すなわち、ステップS106によって、PCは、「k」を「+1」として、「k」が「kmax」となるまで、ステップS103乃至ステップS105を繰り返し行う。
<< Calculation Example of “k = k + 1” (Step S106) >>
In step S106, the PC calculates “k = k + 1”. That is, in step S106, the PC repeats steps S103 to S105 until “k” becomes “+1” and “k” becomes “k max ”.
≪ 「i=m」であるか否かの判断例(ステップS107) ≫
ステップS107では、PCは、「i=m」であるか否かを判断する。そして、「i=m」でないとPCが判断すると(ステップS107でNO)、PCは、ステップS108に進む。一方で、「i=m」であるとPCが判断すると(ステップS107でYES)、PCは、ステップS109に進む。
<< Judgment example of whether or not "i = m" (step S107) >>
In step S107, the PC determines whether or not “i = m”. If the PC determines that “i = m” is not satisfied (NO in step S107), the PC proceeds to step S108. On the other hand, if the PC determines that “i = m” (YES in step S107), the PC proceeds to step S109.
≪ 「i=i+1」の計算例(ステップS108) ≫
ステップS108では、PCは、「i=i+1」を計算する。すなわち、ステップS108によって、PCは、「i」を「+1」として、「i」が「m」となるまで、ステップS102乃至ステップS107を繰り返し行う。
<< Calculation Example of “i = i + 1” (Step S108) >>
In step S108, the PC calculates “i = i + 1”. That is, in step S108, the PC sets “i” to “+1” and repeats steps S102 to S107 until “i” becomes “m”.
以上のように、ステップS105乃至ステップS108によって、ステップごと、かつ、学習データごとに、モデル及び評価指標が計算される。そして、あらかじめ設定される判定条件となる「k=kmax」かつ「i=m」となると、モデル及び評価指標を計算する繰り返しが終了する。 As described above, the model and the evaluation index are calculated for each step and for each learning data in steps S105 to S108. Then, when “k = k max ” and “i = m”, which are determination conditions set in advance, are satisfied, the iteration of calculating the model and the evaluation index is completed.
以上のような処理を図示すると、例えば、以下のように示せる。 The above process can be illustrated as follows, for example.
図4は、本発明の一実施形態における全体処理例を示す模式図である。図では、横軸を時間とし、縦軸を「i」、すなわち、1回の図2に示すステップS102乃至ステップS104の処理結果とする。したがって、図では、縦軸において、上から(「i=1」からとする。)下に向かって(すなわち、「i=m」となるまで)繰り返しステップS102乃至ステップS104が行われるとする。 FIG. 4 is a schematic diagram showing an example of overall processing in an embodiment of the present invention. In the figure, the horizontal axis is time, and the vertical axis is “i”, that is, the processing results of steps S102 to S104 shown in FIG. Accordingly, in the figure, it is assumed that steps S102 to S104 are repeatedly performed from the top (from “i = 1”) toward the bottom (that is, until “i = m”) on the vertical axis.
例えば、図示するように、現時点TCより過去の実績値及び予測値を含む学習データDLがあらかじめPCに入力される(図2に示すステップS101)。そして、図示するように、学習データDLは、所定の期間ごとに分割され、分割によって、あらかじめ期間データDTが生成されるとする(図2に示すステップS102)。 For example, as shown in the figure, learning data DL including past actual values and predicted values from the current TC is input to the PC in advance (step S101 shown in FIG. 2). As illustrated, the learning data DL is divided every predetermined period, and the period data DT is generated in advance by the division (step S102 shown in FIG. 2).
図示する例では、まず、「i=1」において、PCは、第1期間TM1の期間データDTを学習データDLから選択し、第1期間TM1の期間データDTを用いて、第1モデルf1を計算する(図2に示すステップS103)。このようにすると、「i=1」において、PCは、第1モデルf1に基づいて、例えば、上記(5)式を用いて、第1評価指標EV1を計算する(図2に示すステップS104)。 In the illustrated example, first, in “i = 1”, the PC selects the period data DT for the first period TM1 from the learning data DL, and uses the period data DT for the first period TM1 to determine the first model f1. Calculation is performed (step S103 shown in FIG. 2). In this way, at “i = 1”, the PC calculates the first evaluation index EV1 based on the first model f1, for example, using the equation (5) (step S104 shown in FIG. 2). .
次に、図2に示すステップS108によって、PCは、「i=2」とする。そして、「i=2」において、PCは、「i=1」と同様に、まず、期間データDTを選択する。図示する例では、PCは、第1期間TM1及び第2期間TM2の期間データDTを学習データDLから選択する(図2に示すステップS102)。次に、PCは、第1期間TM1及び第2期間TM2の期間データDTを用いて、第2モデルf2を計算する(図2に示すステップS103)。このようにすると、「i=2」において、PCは、第2モデルf2に基づいて、例えば、上記(5)式を用いて、第2評価指標EV2を計算する(図2に示すステップS104)。 Next, in step S108 shown in FIG. 2, the PC sets “i = 2”. In “i = 2”, the PC first selects the period data DT as in “i = 1”. In the illustrated example, the PC selects the period data DT for the first period TM1 and the second period TM2 from the learning data DL (step S102 shown in FIG. 2). Next, the PC calculates the second model f2 using the period data DT of the first period TM1 and the second period TM2 (step S103 shown in FIG. 2). In this way, at “i = 2”, the PC calculates the second evaluation index EV2 based on the second model f2 using, for example, the above equation (5) (step S104 shown in FIG. 2). .
以下、図示するように、PCは、「i=m」となるまで繰り返し「i=1」等と同様の処理を行う。 Hereinafter, as illustrated, the PC repeatedly performs the same processing as “i = 1” or the like until “i = m”.
なお、PCは、予測時点ごとにそれぞれのモデルを計算、すなわち、予測時点ごとに、それぞれの予測に用いる期間データをそれぞれ決定するのが望ましい。具体的には、「k=1」のステップにおけるモデルを計算するのに用いられる期間データと、「k=2」のステップにおけるモデルを計算するのに用いられる期間データとは、別々に決定されるのが望ましい。したがって、ステップごとに、モデルの計算等は、独立して別個に行われるのが望ましい。具体的には、図2に図示するように、ステップS106で「k=k+1」となるたびに、ステップS103でモデルが計算されるのが望ましい。 It is desirable that the PC calculates each model for each prediction time point, that is, determines the period data used for each prediction for each prediction time point. Specifically, the period data used to calculate the model at the step “k = 1” and the period data used to calculate the model at the step “k = 2” are determined separately. Is desirable. Therefore, it is desirable that the calculation of the model and the like be performed independently for each step. Specifically, as illustrated in FIG. 2, it is desirable that the model is calculated in step S103 every time “k = k + 1” in step S106.
このように、ステップごとに、それぞれの期間データが選択されると、用いられる期間データが異なる場合があるため、計算によって生成されるモデルも、ステップごとにそれぞれ異なるモデルとなる場合がある。このようにして、予測時点ごとに最適な期間データが決定されると、PCは、それぞれの予測時点により適したモデルを計算することができる。 As described above, when each period data is selected for each step, the period data to be used may be different. Therefore, a model generated by calculation may be a different model for each step. In this way, when optimum period data is determined for each prediction time point, the PC can calculate a model more suitable for each prediction time point.
≪ 予測誤差の信頼区間に基づく期間データの決定例(ステップS109) ≫
図2に戻り、ステップS109では、PCは、まず、予測誤差の信頼区間を決定する。予測誤差の信頼区間は、例えば、以下のように決定される。
<< Example of determining period data based on confidence interval of prediction error (step S109) >>
Returning to FIG. 2, in step S109, the PC first determines a confidence interval of the prediction error. The confidence interval of the prediction error is determined as follows, for example.
図5は、本発明の一実施形態における予測誤差の信頼区間の一例を示す図である。例えば、モデルに基づいて、上記(4)式のように、モデルに基づいて標準偏差が計算されると、予測誤差の信頼区間の上限値と、下限値とが決定するため、予測誤差の信頼区間が定まる。 FIG. 5 is a diagram illustrating an example of a confidence interval for a prediction error in an embodiment of the present invention. For example, when the standard deviation is calculated based on the model based on the model as shown in the above equation (4), the upper limit value and the lower limit value of the prediction error confidence interval are determined. The section is determined.
図5(A)に示すように、例えば、標準偏差が「1σ」となる予測誤差の信頼区間には、約68.27パーセントの確率で実績値が含まれることになる。また、予測誤差の信頼区間は、例えば、図5(B)のようにも図示できる。図示するように、予測誤差の信頼区間は、上限値と下限値とによって定まる予測誤差の平均値を中心とした幅である。なお、予測値と、予測誤差の平均値は、同意である。 As shown in FIG. 5A, for example, the confidence interval of the prediction error where the standard deviation is “1σ” includes the actual value with a probability of about 68.27%. The confidence interval of the prediction error can also be illustrated as shown in FIG. 5B, for example. As shown in the figure, the confidence interval for the prediction error is a width centered on the average value of the prediction error determined by the upper limit value and the lower limit value. Note that the predicted value and the average value of the prediction error are agreement.
そして、図5(B)に示すように、現時点TCから「kステップ」先の予測時点において、予測値を計算するのに、決定された予測誤差の信頼区間となるように、例えば、以下のように学習データが選定される。 Then, as shown in FIG. 5B, at the prediction time point “k steps” ahead of the current TC, the prediction value is calculated so that the confidence interval of the determined prediction error is obtained. Learning data is selected as follows.
図6は、本発明の一実施形態における予測誤差の信頼区間に基づく学習データの選定例を示す図である。図では、縦軸をPICPとし、横軸をステップ「k」とする。図は、標準偏差が「1σ」となる予測誤差の信頼区間にする例である。 FIG. 6 is a diagram showing an example of selecting learning data based on the prediction error confidence interval in the embodiment of the present invention. In the figure, the vertical axis is PICP and the horizontal axis is step “k”. The figure shows an example of a prediction error confidence interval in which the standard deviation is “1σ”.
図示するように、学習データの期間データが、現時点TCから直近「15日」分、直近「30日」分、直近「45日」分、直近「60日」分及び直近「75日」分のデータである例で説明する。そして、この例は、各ステップ「k」において、各期間データに基づいて計算されるモデルに対するPICPの計算結果を示す。 As shown in the figure, the period data of the learning data is the latest “15 days”, the latest “30 days”, the latest “45 days”, the latest “60 days”, and the latest “75 days” from the current TC. An example of data will be described. This example shows the PICP calculation result for the model calculated based on each period data at each step “k”.
そして、上記(7)式に基づいて、各PICPの計算結果のうち、理論値に最も近い計算結果となる期間データが、ステップごとに選定される。すなわち、選定結果DRは、各ステップにおいて、最も理論値に近いPICPとなった期間データの組み合わせである。 Then, based on the above equation (7), among the calculation results of each PICP, the period data that is the calculation result closest to the theoretical value is selected for each step. That is, the selection result DR is a combination of period data having a PICP closest to the theoretical value at each step.
具体的には、ステップが「1」(すなわち、「k=1」である。)の場合には、「30日」分の期間データを用いると、PICPが理論値に最も近い値であるため、この例では、ステップ「1」では、「30日」分の期間データが選定される選定結果DRとなる。次に、ステップが「2」(すなわち、「k=2」である。)の場合には、「45日」分の期間データを用いると、PICPが理論値に最も近い値であるため、この例では、ステップ「2」では、「45日」分の期間データが選定される選定結果DRとなる。このようにして、各ステップにおいて、PICPに基づいて、いずれかの期間データが選定されて、選定結果DRが定まる。 Specifically, when the step is “1” (that is, “k = 1”), if period data for “30 days” is used, PICP is the closest value to the theoretical value. In this example, in step “1”, a selection result DR is selected in which period data for “30 days” is selected. Next, when the step is “2” (that is, “k = 2”), using the period data for “45 days”, the PICP is the closest value to the theoretical value. In the example, in the step “2”, the period data for “45 days” is selected as the selection result DR. In this way, in each step, any period data is selected based on the PICP, and the selection result DR is determined.
以上のように決定される期間データを用いて、各ステップでは、PCによって、風速が予測される。つまり、図示する例では、ステップ「1」では、学習データのうち、直近「30日」分の期間データが用いられて風速の予測が行われる。このように、PICP等の評価指標に基づいて決定された学習データが用いられると、PCは、精度良く予測誤差の信頼区間を予測することができる。 Using the period data determined as described above, the wind speed is predicted by the PC in each step. In other words, in the illustrated example, in step “1”, the wind speed is predicted using the latest “30 days” of the learning data. As described above, when learning data determined based on an evaluation index such as PICP is used, the PC can accurately predict a confidence interval of a prediction error.
なお、選定される期間データは、現時点から連続するデータでなくともよい。すなわち、例えば、PCは、現時点から直近の期間データを選択しないでもよい。例えば、直近「1日」分、すなわち、昨日の期間データが選定されず、直近「2日」から「16日」までの「15日」分等というように、PCは、直近の期間データが含まれない選定を行ってもよい。直近の期間データが選定されると、季節性の影響が大きくなる場合がある。そこで、季節性の影響を少なくするため、PCは、現時点から直近の期間データを選択しない選定を行ってもよい。 The selected period data may not be continuous data from the present time. That is, for example, the PC may not select data for the most recent period from the current time. For example, as for the latest “1 day”, that is, yesterday's period data is not selected, and the latest period data is “15 days” from the latest “2 days” to “16 days”, etc. Selections that are not included may be made. When the latest period data is selected, the influence of seasonality may increase. Therefore, in order to reduce the influence of seasonality, the PC may perform selection not to select data for the most recent period from the current time.
図7は、本発明の一実施形態における効果の一例を示す図である。まず、図示するように、PICPが計算され、PICPが最も理論値に近くなる期間データDTが選定される。次に、PCは、選定された期間データDTを用いて、風速を予測する。このようにすると、PCは、精度良く予測誤差の信頼区間を予測することができる。 FIG. 7 is a diagram illustrating an example of an effect in one embodiment of the present invention. First, as shown in the figure, PICP is calculated, and period data DT in which PICP is closest to the theoretical value is selected. Next, the PC predicts the wind speed using the selected period data DT. In this way, the PC can accurately predict the confidence interval of the prediction error.
なお、評価指標は、PICPに限られない。例えば、評価指標は、NMPIL(Normalized Mean Prediction Interval Length)等でもよい。PICPを用いると、PCは、分布の信頼度が評価できるのに対して、NMPILを用いると、PCは、分布の分離度を評価することができる。NMPILの計算方法等は、例えば、「Abbas Khosravi, Saeid Nahavandi and Doug Creighton,"Construction of Optimal Prediction Intervals for Load Forecasting Probrem",IEEE Transactions on Power Systems, Vol.25, Issue 3, pp.1496―1503(2010)」に記載されている通りである。
Note that the evaluation index is not limited to PICP. For example, the evaluation index may be NMPIL (Normalized Mean Prediction Interval Length) or the like. When PICP is used, the PC can evaluate the reliability of the distribution, whereas when NMPIL is used, the PC can evaluate the degree of separation of the distribution. Calculation method and the like of NMPIL is, for example, "Abbas Khosravi, Saeid Nahavandi and Doug Creighton," Construction of Optimal Prediction Intervals for Load Forecasting Probrem ", IEEE Transactions on Power Systems, Vol.25,
≪ 3. 風速予測装置の機能構成例 ≫
図8は、本発明の一実施形態における風速予測装置の機能構成の一例を示す機能ブロック図である。例えば、PC10は、図示するように、入力部FN1と、計算部FN2と、評価指標計算部FN3と、決定部FN4とを含む機能構成である。
≪ 3. Example of functional configuration of wind speed prediction device ≫
FIG. 8 is a functional block diagram illustrating an example of a functional configuration of the wind speed prediction apparatus according to the embodiment of the present invention. For example, the
入力部FN1は、実績値及び予測値を示す学習データDLを入力する。例えば、入力部FN1は、ネットワークI/F103(図1参照)又は入力I/F104(図1参照)等によって実現される。 The input unit FN1 inputs learning data DL indicating the actual value and the predicted value. For example, the input unit FN1 is realized by the network I / F 103 (see FIG. 1) or the input I / F 104 (see FIG. 1).
計算部FN2は、学習データDLに基づいて、まず、モデルを計算する。次に、計算部FN2は、モデルに基づいて、それぞれの予測時点におけるそれぞれの予測誤差の信頼区間を計算する。例えば、計算部FN2は、CPU101(図1参照)等によって実現される。 The calculation unit FN2 first calculates a model based on the learning data DL. Next, the calculation unit FN2 calculates a confidence interval of each prediction error at each prediction time based on the model. For example, the calculation unit FN2 is realized by the CPU 101 (see FIG. 1) or the like.
評価指標計算部FN3は、計算部FN2が計算する予測誤差の信頼区間のそれぞれの評価指標を計算する。例えば、評価指標計算部FN3は、CPU101(図1参照)等によって実現される。 The evaluation index calculation unit FN3 calculates each evaluation index of the confidence interval of the prediction error calculated by the calculation unit FN2. For example, the evaluation index calculation unit FN3 is realized by the CPU 101 (see FIG. 1) or the like.
決定部FN4は、評価指標計算部FN3が計算する評価指標に基づいて、学習データDLのうち、予測に用いる期間データを決定する。例えば、決定部FN4は、CPU101(図1参照)等によって実現される。 The determination unit FN4 determines period data used for prediction in the learning data DL based on the evaluation index calculated by the evaluation index calculation unit FN3. For example, the determination unit FN4 is realized by the CPU 101 (see FIG. 1) or the like.
≪ まとめ ≫
まず、PC10は、入力部FN1によって、あらかじめ計測した風速の実績値と、実績値を予測した予測値とを示す学習データDLを入力する。そして、学習データDLが入力されると、図4に示すように、PC10は、計算部FN2によって、モデルを計算することができる。次に、モデルが計算されると、PCは、モデルに基づいて、標準偏差又は標準偏差の整数倍を計算し、予測誤差の信頼区間を計算する。
≪ Summary ≫
First, the
続いて、PC10は、評価指標計算部FN3によって、例えば、モデルに基づいて、上記(5)式等により、評価指標を計算する。このように、評価指標が計算されると、PCは、例えば、上記(7)式に基づいて、決定部FN4によって、図6に示すように最も理論値に近いPICPとなる期間データを予測に用いるように決定できる。
Subsequently, the
このように、予測に用いる期間データを決定できると、PCは、風速の予測において、予測誤差の信頼区間を精度良く予測することができる。 Thus, if the period data used for prediction can be determined, the PC can accurately predict the confidence interval of the prediction error in the wind speed prediction.
精度が良い予測誤差の信頼区間であると、例えば、風力発電において蓄電池の容量等を決定する目安等にすることができる。具体的には、蓄電池は、風力発電によって発電された電力を蓄電するため、設置の際、発電量に応じて蓄電できる容量を決める必要がある。つまり、予測される最大風速等に応じて蓄電池の容量が決定される場合がある。そこで、予測誤差の信頼区間が精度良く予測されると、無駄に容量が大きい蓄電池が選ばれてしまうこと等を少なくすることができる。 If the confidence interval of the prediction error has good accuracy, it can be used as a guideline for determining the capacity of the storage battery in wind power generation, for example. Specifically, since the storage battery stores the electric power generated by wind power generation, it is necessary to determine the capacity that can be stored according to the amount of power generation at the time of installation. That is, the capacity of the storage battery may be determined according to the predicted maximum wind speed or the like. Therefore, when the confidence interval of the prediction error is predicted with high accuracy, it is possible to reduce a case where a storage battery having a large capacity is unnecessarily selected.
なお、本発明の一実施形態に係る各処理の全部又は一部は、アセンブラ等の低水準言語、C言語、Java(登録商標)言語等の高水準言語又はこれらを組み合わせて記述されるプログラムによって実現されてもよい。すなわち、プログラムは、情報処理装置又は情報処理システム等のコンピュータに風速予測方法に係る各手順を実行させるためのコンピュータプログラムである。 Note that all or part of each processing according to an embodiment of the present invention is performed by a low-level language such as an assembler, a high-level language such as C language, Java (registered trademark) language, or a program described in combination of these. It may be realized. That is, the program is a computer program for causing a computer such as an information processing apparatus or an information processing system to execute each procedure related to the wind speed prediction method.
さらに、本発明の一実施形態に係る各処理の全部又は一部は、1以上の情報処理装置を有する風速予測システムによって、処理の全部又は一部が並行、分散、冗長又はこれらの組み合わせで処理されてもよい。 Further, all or a part of each process according to an embodiment of the present invention is processed by a wind speed prediction system having one or more information processing devices, and all or a part of the process is processed in parallel, distributed, redundant, or a combination thereof. May be.
また、本発明の一実施形態に係る各処理は、図示した順序に限られない。例えば、各処理の一部又は全部は、異なる順序、並行、分散又は省略されて処理されてもよい。 Moreover, each process which concerns on one Embodiment of this invention is not restricted to the order shown in figure. For example, some or all of the processes may be processed in different orders, in parallel, distributed, or omitted.
以上、本発明の好ましい実施例について詳述したが、本発明は、上述の実施形態に限定されず、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形又は変更が可能である。 The preferred embodiments of the present invention have been described in detail above, but the present invention is not limited to the above-described embodiments, and various modifications or changes can be made within the scope of the gist of the present invention described in the claims. Is possible.
10 PC
FN1 入力部
FN2 計算部
FN3 評価指標計算部
FN4 決定部
DL 学習データ
DT 期間データ
DR 選定結果
CAL 予測値
RES 実績値
10 PC
FN1 Input unit FN2 Calculation unit FN3 Evaluation index calculation unit FN4 Determination unit DL Learning data DT Period data DR Selection result CAL Predicted value RES Actual value
Claims (7)
現時点より過去の時点における風速を計測した実績値及び前記実績値を予測した予測値を示す学習データを入力する入力部と、
前記学習データに基づいて、前記現時点から所定の時間後となる予測時点における風速のモデルと、前記モデルに基づいて前記予測時点における予測誤差の信頼区間とを少なくとも計算する計算部と、
前記予測誤差の信頼区間の評価指標を計算する評価指標計算部と、
前記評価指標に基づいて、前記学習データのうち、予測に用いる期間データを決定する決定部と
を含む風速予測装置。 A wind speed prediction device for predicting wind speed,
An input unit for inputting learning data indicating an actual value obtained by measuring the wind speed at a time point in the past from the present time and a predicted value obtained by predicting the actual value;
Based on the learning data, a calculation unit for calculating at least a wind speed model at a prediction time after a predetermined time from the current time, and a prediction error confidence interval at the prediction time based on the model;
An evaluation index calculation unit for calculating an evaluation index of the confidence interval of the prediction error;
A wind speed prediction apparatus including a determination unit that determines period data to be used for prediction out of the learning data based on the evaluation index.
現時点より過去の時点における風速を計測した実績値及び前記実績値を予測した予測値を示す学習データを入力する入力部と、
前記学習データに基づいて、前記現時点から所定の時間後となる予測時点における風速のモデルと、前記モデルに基づいて前記予測時点における予測誤差の信頼区間とを少なくとも計算する計算部と、
前記予測誤差の信頼区間の評価指標を計算する評価指標計算部と、
前記評価指標に基づいて、前記学習データのうち、予測に用いる期間データを決定する決定部と
を含む風速予測システム。 A wind speed prediction system having one or more information processing devices and predicting wind speed,
An input unit for inputting learning data indicating an actual value obtained by measuring the wind speed at a time point in the past from the present time and a predicted value obtained by predicting the actual value;
Based on the learning data, a calculation unit for calculating at least a wind speed model at a prediction time after a predetermined time from the current time, and a prediction error confidence interval at the prediction time based on the model;
An evaluation index calculation unit for calculating an evaluation index of the confidence interval of the prediction error;
A wind speed prediction system including a determination unit that determines period data used for prediction among the learning data based on the evaluation index.
前記風速予測装置が、現時点より過去の時点における風速を計測した実績値及び前記実績値を予測した予測値を示す学習データを入力する入力手順と、
前記風速予測装置が、前記学習データに基づいて、前記現時点から所定の時間後となる予測時点における風速のモデルと、前記モデルに基づいて前記予測時点における予測誤差の信頼区間とを少なくとも計算する計算手順と、
前記風速予測装置が、前記予測誤差の信頼区間の評価指標を計算する評価指標計算手順と、
前記風速予測装置が、前記評価指標に基づいて、前記学習データのうち、予測に用いる期間データを決定する決定手順と
を含む風速予測方法。 A wind speed prediction method performed by a wind speed prediction device that predicts wind speed,
Input procedure for inputting learning data indicating a predicted value obtained by predicting the actual value and the actual value obtained by measuring the actual wind speed from the present time by the wind speed prediction device;
The wind speed prediction device calculates at least a wind speed model at a prediction time after a predetermined time from the current time based on the learning data, and a prediction error confidence interval at the prediction time based on the model. Procedure and
The wind speed prediction apparatus calculates an evaluation index calculation procedure for calculating an evaluation index of a confidence interval of the prediction error;
A wind speed prediction method, wherein the wind speed prediction device includes a determination procedure for determining period data used for prediction among the learning data based on the evaluation index.
前記コンピュータが、現時点より過去の時点における風速を計測した実績値及び前記実績値を予測した予測値を示す学習データを入力する入力手順と、
前記コンピュータが、前記学習データに基づいて、前記現時点から所定の時間後となる予測時点における風速のモデルと、前記モデルに基づいて前記予測時点における予測誤差の信頼区間とを少なくとも計算する計算手順と、
前記コンピュータが、前記予測誤差の信頼区間の評価指標を計算する評価指標計算手順と、
前記コンピュータが、前記評価指標に基づいて、前記学習データのうち、予測に用いる期間データを決定する決定手順と
を実行させるためのプログラム。 A program for causing a computer for predicting wind speed to execute a wind speed prediction method,
An input procedure for inputting learning data indicating an actual value obtained by measuring the wind speed at a time point in the past from the present time and a predicted value obtained by predicting the actual value;
A calculation procedure in which the computer calculates at least a wind speed model at a prediction time point after a predetermined time from the current time based on the learning data, and a prediction error confidence interval at the prediction time point based on the model; ,
An evaluation index calculation procedure in which the computer calculates an evaluation index of a confidence interval of the prediction error;
The program for making the said computer perform the determination procedure which determines the period data used for prediction among the said learning data based on the said evaluation parameter | index.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016170918A JP6848271B2 (en) | 2016-09-01 | 2016-09-01 | Wind speed prediction device, wind speed prediction system, wind speed prediction method and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016170918A JP6848271B2 (en) | 2016-09-01 | 2016-09-01 | Wind speed prediction device, wind speed prediction system, wind speed prediction method and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018036196A true JP2018036196A (en) | 2018-03-08 |
JP6848271B2 JP6848271B2 (en) | 2021-03-24 |
Family
ID=61564831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016170918A Active JP6848271B2 (en) | 2016-09-01 | 2016-09-01 | Wind speed prediction device, wind speed prediction system, wind speed prediction method and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6848271B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583075A (en) * | 2018-11-26 | 2019-04-05 | 湖南科技大学 | Permanent magnet direct-drive wind-force machine military service quality evaluating method based on temperature parameter prediction |
CN110707744A (en) * | 2019-10-09 | 2020-01-17 | 广东电网有限责任公司电网规划研究中心 | Method and device for monitoring power generation state of wind power plant cluster and storage medium |
CN112132360A (en) * | 2020-09-29 | 2020-12-25 | 中交建冀交高速公路投资发展有限公司 | Wind speed short-time confidence interval prediction method and system based on LSTM and GMM |
JP7431377B1 (en) | 2023-08-29 | 2024-02-14 | 株式会社日立パワーソリューションズ | Wind condition observation equipment placement support device and wind condition observation equipment placement support method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6975925B1 (en) * | 2002-03-19 | 2005-12-13 | Windlynx Systems, B.V. | Forecasting an energy output of a wind farm |
JP2007056686A (en) * | 2005-08-22 | 2007-03-08 | Univ Of Ryukyus | Device for predicting power generated after several hours on basis of predicted wind velocity in wind power generator |
JP2007233639A (en) * | 2006-02-28 | 2007-09-13 | Tohoku Electric Power Co Inc | Wind power generation output prediction method, wind power generation output prediction unit, and program |
JP2014026625A (en) * | 2012-07-30 | 2014-02-06 | Yoshito Hirata | Prediction mechanism for power variation in recyclable energy based on real time prediction method for weather change |
CN105046374A (en) * | 2015-08-25 | 2015-11-11 | 华北电力大学 | Power interval predication method based on nucleus limit learning machine model |
WO2016121202A1 (en) * | 2015-01-30 | 2016-08-04 | 日本電気株式会社 | Prediction device, prediction method, and program |
-
2016
- 2016-09-01 JP JP2016170918A patent/JP6848271B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6975925B1 (en) * | 2002-03-19 | 2005-12-13 | Windlynx Systems, B.V. | Forecasting an energy output of a wind farm |
JP2007056686A (en) * | 2005-08-22 | 2007-03-08 | Univ Of Ryukyus | Device for predicting power generated after several hours on basis of predicted wind velocity in wind power generator |
JP2007233639A (en) * | 2006-02-28 | 2007-09-13 | Tohoku Electric Power Co Inc | Wind power generation output prediction method, wind power generation output prediction unit, and program |
JP2014026625A (en) * | 2012-07-30 | 2014-02-06 | Yoshito Hirata | Prediction mechanism for power variation in recyclable energy based on real time prediction method for weather change |
WO2016121202A1 (en) * | 2015-01-30 | 2016-08-04 | 日本電気株式会社 | Prediction device, prediction method, and program |
CN105046374A (en) * | 2015-08-25 | 2015-11-11 | 华北电力大学 | Power interval predication method based on nucleus limit learning machine model |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583075A (en) * | 2018-11-26 | 2019-04-05 | 湖南科技大学 | Permanent magnet direct-drive wind-force machine military service quality evaluating method based on temperature parameter prediction |
CN109583075B (en) * | 2018-11-26 | 2022-12-02 | 湖南科技大学 | Permanent magnet direct-drive wind turbine service quality evaluation method based on temperature parameter prediction |
CN110707744A (en) * | 2019-10-09 | 2020-01-17 | 广东电网有限责任公司电网规划研究中心 | Method and device for monitoring power generation state of wind power plant cluster and storage medium |
CN112132360A (en) * | 2020-09-29 | 2020-12-25 | 中交建冀交高速公路投资发展有限公司 | Wind speed short-time confidence interval prediction method and system based on LSTM and GMM |
JP7431377B1 (en) | 2023-08-29 | 2024-02-14 | 株式会社日立パワーソリューションズ | Wind condition observation equipment placement support device and wind condition observation equipment placement support method |
Also Published As
Publication number | Publication date |
---|---|
JP6848271B2 (en) | 2021-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20200063710A1 (en) | System and methods for hyper short-term wind power prediction using real-time wind parameter measurements | |
JP5888640B2 (en) | Photovoltaic power generation prediction apparatus, solar power generation prediction method, and solar power generation prediction program | |
KR101875329B1 (en) | Forecasting apparatus and method for power consumption | |
CN117039895B (en) | Wind power prediction method and system for energy storage auxiliary black start | |
JP6848271B2 (en) | Wind speed prediction device, wind speed prediction system, wind speed prediction method and program | |
JP2004086896A (en) | Method and system for constructing adaptive prediction model | |
CN112308281A (en) | Temperature information prediction method and device | |
JP7010674B2 (en) | Power Demand Forecasting Device, Power Demand Forecasting Method and Program | |
CN110474808A (en) | A kind of method for predicting and device | |
Langsari et al. | Optimizing effort and time parameters of COCOMO II estimation using fuzzy multi-objective PSO | |
CN110969306A (en) | Power distribution low-voltage distribution area load prediction method and device based on deep learning | |
CN106980874B (en) | A kind of multi-time Scales dimension data fusion method towards distribution big data | |
CN112396231A (en) | Modeling method and device for spatio-temporal data, electronic equipment and readable medium | |
Feng et al. | Scenario reduction for stochastic unit commitment with wind penetration | |
CN105844353A (en) | Aquatic product price prediction method and device | |
WO2015145978A1 (en) | Energy-amount estimation device, energy-amount estimation method, and recording medium | |
CN111144721A (en) | Power grid project demand evaluation model construction method and device and computing equipment | |
CN104010029A (en) | Distributed computing environment performance predicting method based on transverse and longitudinal information integration | |
US10998725B2 (en) | Electric power generation prediction method based on expected value calculation, electric power generation prediction system based on expected value calculation, and electric power generation prediction program product based on expected value calculation | |
Langsari et al. | Optimizing COCOMO II parameters using particle swarm method | |
CN111209930A (en) | Method and device for generating credit granting strategy and electronic equipment | |
JP6877956B2 (en) | Demand forecasting device and demand forecasting method | |
Hu et al. | Stochastic multi-attribute group decision-making for long-term comprehensive operation of cascade reservoirs under multiple uncertainties | |
WO2020166300A1 (en) | Measurement guide device, and simulation computation device used in same | |
CN118523200B (en) | Suspension cable type crossing frame sag self-adaption method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190809 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200806 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200818 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201015 |
|
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: 20210202 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210215 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6848271 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |