JP2022183783A - Information processing device, information processing method, program, and power generation system - Google Patents

Information processing device, information processing method, program, and power generation system Download PDF

Info

Publication number
JP2022183783A
JP2022183783A JP2021091268A JP2021091268A JP2022183783A JP 2022183783 A JP2022183783 A JP 2022183783A JP 2021091268 A JP2021091268 A JP 2021091268A JP 2021091268 A JP2021091268 A JP 2021091268A JP 2022183783 A JP2022183783 A JP 2022183783A
Authority
JP
Japan
Prior art keywords
power
target
time period
information processing
generated
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
Application number
JP2021091268A
Other languages
Japanese (ja)
Other versions
JP7043661B1 (en
Inventor
博正 進
Hiromasa Shin
満 柿元
Mitsuru Kakimoto
慶明 志賀
Yoshiaki Shiga
健史 大場
Takeshi Oba
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
Toshiba Energy Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Energy Systems and Solutions 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 Toshiba Corp, Toshiba Energy Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2021091268A priority Critical patent/JP7043661B1/en
Application granted granted Critical
Publication of JP7043661B1 publication Critical patent/JP7043661B1/en
Publication of JP2022183783A publication Critical patent/JP2022183783A/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Photovoltaic Devices (AREA)

Abstract

To predict an amount of power generated by a power plant.SOLUTION: An information processing device includes a position calculation unit, a power amount calculation unit, and an output unit. The position calculation unit calculates a target position of a first time zone result value of an amount of power generated by one or more power plants on a target day with respect to a first sample group including a plurality of first time zone result values of an amount of power generated by the one or more power plants on a plurality of reference days. The power amount calculation unit calculates a power amount corresponding to the target position, in a second sample group including a plurality of prediction target time zone result values of an amount of power generated by the one or more power plants on the plurality of reference days. The output unit outputs the power amount corresponding to the target position, as a predicted power amount which is predicted to be generated by the one or more power plants in the prediction target time zone.SELECTED DRAWING: Figure 1

Description

本発明の実施形態は、情報処理装置、情報処理方法、プログラムおよび発電システムに関する。 TECHNICAL FIELD Embodiments of the present invention relate to an information processing device, an information processing method, a program, and a power generation system.

太陽光発電により発電された電力の買取方法が、固定価格の買取制度(FIT)から市場価格連動の買取制度(FIP)に移行する。市場価格連動の買取制度においては、発電量の計画値と実績値とを一致させるために、例えば、電力供給の1時間前に電力取引を行う1時間前市場の活用が検討されている。このため、市場価格連動の買取制度においては、比較的に近い将来に発電する電力量を予測する短時間予測技術の高精度化が求められている。 The method of purchasing electricity generated by photovoltaic power generation will shift from a fixed-price purchase system (FIT) to a market-price-linked purchase system (FIP). In the purchase system linked to the market price, in order to match the planned value and the actual value of the generated power, for example, utilization of the one-hour-ahead market, in which power is traded one hour before the power supply, is being considered. For this reason, in the purchase system linked to the market price, there is a demand for high-precision short-term prediction technology for predicting the amount of power to be generated in the relatively near future.

従来、太陽光発電における短時間予測技術として、発電実績値のみを利用する方法、並びに、発電実績値に加えて気象衛星および気象予報を利用する方法が存在する。 Conventionally, as short-term forecasting techniques in photovoltaic power generation, there are a method of using only actual power generation values and a method of using weather satellites and weather forecasts in addition to actual power generation values.

発電実績値のみを利用する方法として、快晴指数の持続予測方式が知られている。しかし、快晴指数の持続予測方式は、日射強度と発電量との相互変換が複雑であり、誤差も大きい。 As a method of using only actual power generation values, a continuous forecasting method of a clear weather index is known. However, the continuous forecasting method of the clear weather index is complicated in mutual conversion between solar radiation intensity and power generation amount, and the error is large.

発電実績値に加えて気象衛星および気象予報を利用する方法として、気象の類似日を検索する方式が知られている。しかし、この方式は、外部から取得するデータが増加することに伴いデータ処理が複雑となり、この結果、データがリアルタイムで取得できない等の不具合が生じるリスクも増加する。 As a method of using weather satellites and weather forecasts in addition to actual power generation values, a method of searching for days with similar weather is known. However, this method complicates data processing as the amount of data acquired from the outside increases, and as a result, the risk of problems such as failure to acquire data in real time also increases.

特開2020-108188号公報Japanese Patent Application Laid-Open No. 2020-108188 特開2013-84736号公報JP 2013-84736 A 特許第6785971号公報Japanese Patent No. 6785971

Stuart Coles, “An Introduction to Statistical Modeling of Extreme Values”, Springer, 2001Stuart Coles, “An Introduction to Statistical Modeling of Extreme Values,” Springer, 2001 Sokol Dervishi and Ardeshir Mahdavi, "Computing diffuse fraction of global horizontal solar radiation: A model comparison", in Solar Energy. 2012 Jun; 86(6): 1796-1802.Sokol Dervishi and Ardeshir Mahdavi, "Computing diffuse fraction of global horizontal solar radiation: A model comparison", in Solar Energy. 2012 Jun; 86(6): 1796-1802.

本発明が解決しようとする課題は、発電所により発電される電力量を、少ない種類のデータを用いて簡易に精度良く予測する情報処理装置、情報処理方法、プログラムおよび発電システムを提供することにある。 The problem to be solved by the present invention is to provide an information processing device, an information processing method, a program, and a power generation system that easily and accurately predict the amount of electric power generated by a power plant using a small number of types of data. be.

実施形態に係る情報処理装置は、位置算出部と、電力量算出部と、出力部と、を備える。前記位置算出部は、1または複数の発電所により発電された電力量の複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群に対する、前記1または複数の発電所により発電された電力量の対象日における前記第1時間帯の実績値の対象位置を算出する。前記電力量算出部は、前記1または複数の発電所により発電された電力量の前記複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群における、前記対象位置に対応する電力量を算出する。前記出力部は、前記対象位置に対応する電力量を、前記予測対象時間帯において前記1または複数の発電所により発電されると予測される予測電力量として出力する。 An information processing apparatus according to an embodiment includes a position calculator, a power amount calculator, and an output unit. The position calculation unit performs power generation by the one or more power plants with respect to a first sample group including a plurality of actual values of the amount of electric power generated by the one or more power plants in a first time period on a plurality of reference days. The target position of the actual value for the first time period on the target day of the electric energy calculated is calculated. The electric energy calculation unit corresponds to the target position in a second sample group including a plurality of actual values of the electric energy generated by the one or more power plants in the prediction target time period on the plurality of reference days. Calculate the amount of power. The output unit outputs the electric energy corresponding to the target position as a predicted electric energy expected to be generated by the one or more power plants in the prediction target time period.

第1実施形態に係る予測装置の機能構成を示す図。The figure which shows the functional structure of the prediction apparatus which concerns on 1st Embodiment. 予測装置の処理の流れを示すフローチャート。A flow chart showing a flow of processing of a prediction device. 発電所により発電された電力量の時間帯毎の実績値を示す図。FIG. 4 is a diagram showing actual values of electric energy generated by a power plant for each time zone; 予測対象時間帯を示す図。The figure which shows a prediction object time zone. 対象日の第1時間帯の実績値を示す図。The figure which shows the actual value of the 1st time slot|zone of an object day. 複数の参照日の第1時間帯の実績値を示す図。The figure which shows the actual value of the 1st time slot|zone of several reference days. 複数の参照日の予測対象時間帯の実績値を示す図。The figure which shows the actual value of the prediction object time zone of several reference days. 実績値の変化の一例を示す図。The figure which shows an example of the change of an actual value. 第2実施形態に係る予測装置の機能構成を示す図。The figure which shows the functional structure of the prediction apparatus which concerns on 2nd Embodiment. 発電所モデルのブロック構成を示す図。The figure which shows the block structure of a power plant model. 直散分離モデルであるErbs、OrgillおよびReindlの数値例を示す図。Numerical examples of Erbs, Orgill, and Reindl, which are linear separation models. Orgillを用いた場合の直散分離モデルの演算例を示す図。The figure which shows the calculation example of the direct separation model at the time of using Orgill. 太陽方向の数値例を示す図。The figure which shows the example of a numerical value of a sun direction. 傾斜面日射強度の数値例を示す図。The figure which shows the example of a numerical value of an inclined surface solar radiation intensity. 太陽方向の演算例を示す図。The figure which shows the calculation example of a sun direction. モジュール温度と周囲の気温との差を表す図。A diagram showing the difference between the module temperature and the ambient air temperature. 負荷率に対するインバータの変換効率を表す図。The figure showing the conversion efficiency of the inverter with respect to a load factor. 発電所モデルの処理の流れを示すフローチャート。4 is a flowchart showing the flow of processing of a power plant model; 予測装置のハードウェア構成を示す図。The figure which shows the hardware constitutions of a prediction apparatus.

以下、図面を参照しながら実施形態に係る予測装置10について説明する。予測装置10は、1または複数の発電所により発電された過去の電力量に基づき、1または複数の発電所により発電される予測対象時間帯における電力量を算出する。 A prediction device 10 according to an embodiment will be described below with reference to the drawings. The prediction device 10 calculates the amount of power generated by one or more power plants in a prediction target time period based on past power amounts generated by one or more power plants.

(第1実施形態)
図1は、第1実施形態に係る予測装置10の機能構成を示す図である。予測装置10は、例えば単体のコンピュータまたはサーバ装置等のコンピュータである。予測装置10は、1台のコンピュータであってもよいし、クラウドシステムのように複数台のコンピュータにより構成されていてもよい。コンピュータは、所定のプログラムを実行することにより、予測装置10として機能する。
(First embodiment)
FIG. 1 is a diagram showing the functional configuration of a prediction device 10 according to the first embodiment. The prediction device 10 is, for example, a single computer or a computer such as a server device. The prediction device 10 may be one computer, or may be composed of a plurality of computers like a cloud system. The computer functions as the prediction device 10 by executing a given program.

予測装置10は、受信部22と、記憶部24と、入力部26と、対象サンプル取得部28と、第1サンプル群取得部30と、第2サンプル群取得部32と、実績値補正部34と、位置算出部36と、電力量算出部38と、出力部40とを備える。 The prediction device 10 includes a reception unit 22, a storage unit 24, an input unit 26, a target sample acquisition unit 28, a first sample group acquisition unit 30, a second sample group acquisition unit 32, and a performance value correction unit 34. , a position calculator 36 , a power amount calculator 38 , and an output unit 40 .

受信部22は、発電所により発電された電力量の時間帯毎の実績値を受信する。発電所は、本実施形態においては、太陽光発電所である。なお、発電所は、水力発電所、風力発電所または地熱発電所等の他の種類の発電所であってもよい。また、発電所は、電力会社等の事業者により設置された設備であってもよいし、家庭等に設置された自家発電装置であってもよい。 The receiving unit 22 receives the actual value of the amount of electric power generated by the power plant for each time period. The power plant is a photovoltaic power plant in this embodiment. It should be noted that the power plants may be other types of power plants such as hydroelectric power plants, wind power plants or geothermal power plants. Also, the power plant may be a facility installed by a business operator such as an electric power company, or may be a private power generator installed in a home or the like.

また、受信部22は、発電する電力を同一の発電モデルにより推定可能な複数の発電所により発電された電力量の時間帯毎の実績値を受信してもよい。例えば、受信部22は、太陽光から得られるエネルギーがほぼ同一であり、気象条件も同一となる、所在地が近く同一の設備を用いた複数の太陽光発電所により発電された電力量の時間帯毎の実績値を受信してもよい。 In addition, the receiving unit 22 may receive actual values for each time period of the amount of power generated by a plurality of power plants that can estimate the power to be generated by the same power generation model. For example, the receiving unit 22 detects the amount of power generated by a plurality of photovoltaic power plants that are located close to each other and use the same equipment, in which the energy obtained from sunlight is almost the same and the weather conditions are the same. You may receive the performance value for each.

記憶部24は、受信部22により受信された時間帯毎の実績値を記憶する。例えば、記憶部24は、データベースにより管理され、年月日、時間帯および実績値をセットで記憶する。記憶部24は、年月日および時間帯を指定することにより、指定された年月日および時間帯の実績値を検索して出力する。記憶部24は、予測装置10内に実現されてもよいし、予測装置10の外部のデータベースサーバ等に実現されてもよい。 The storage unit 24 stores the performance value for each time period received by the reception unit 22 . For example, the storage unit 24 is managed by a database and stores a set of dates, time zones, and performance values. By designating the year, month, day and time period, the storage unit 24 retrieves and outputs the actual values for the designated year, month, day and time period. The storage unit 24 may be implemented within the prediction device 10 or may be implemented as a database server or the like external to the prediction device 10 .

入力部26は、予測日における予測対象時間帯の入力を受け付ける。入力部26は、例えばユーザにより操作入力された予測対象時間帯を受け付けてもよいし、他の装置から受信した予測対象時間帯を受け付けてもよいし、定期的に自動入力された予測対象時間帯を受け付けてもよい。入力部26は、受け付けた予測対象時間帯を、第2サンプル群取得部32および実績値補正部34に与える。 The input unit 26 receives an input of the prediction target time period on the prediction date. The input unit 26 may accept, for example, a prediction target time period input by an operation by a user, may accept a prediction target time period received from another device, or may receive a prediction target time period automatically input periodically. Obi may be accepted. The input unit 26 provides the received prediction target time period to the second sample group acquisition unit 32 and the actual value correction unit 34 .

対象サンプル取得部28は、記憶部24により記憶されている時間帯毎の実績値のうちの、現在時刻より前の時間帯であって、実質的な発電がされた1つの時間帯における実績値を取得する。例えば、対象サンプル取得部28は、記憶部24により記憶されている時間帯毎の実績値のうちの、実質的な発電がされた直近の時間帯における実績値を取得する。続いて、対象サンプル取得部28は、取得した実績値の時間帯を第1時間帯として特定する。また、対象サンプル取得部28は、取得した実績値の発電日を対象日として特定する。対象サンプル取得部28は、特定した第1時間帯を第1サンプル群取得部30および実績値補正部34に与える。また、対象サンプル取得部28は、取得した実績値を、対象日における第1時間帯の実績値として、位置算出部36に与える。 The target sample acquiring unit 28 obtains actual values in one time period before the current time, among the actual values for each time period stored in the storage unit 24, during which substantial power generation was performed. to get For example, the target sample acquiring unit 28 acquires the actual value in the most recent time period in which substantial power generation was performed among the actual values for each time period stored in the storage unit 24 . Subsequently, the target sample acquiring unit 28 specifies the time zone of the acquired performance values as the first time zone. In addition, the target sample acquisition unit 28 identifies the power generation date of the acquired actual value as the target date. The target sample acquisition unit 28 provides the specified first time period to the first sample group acquisition unit 30 and the performance value correction unit 34 . In addition, the target sample acquisition unit 28 gives the acquired performance value to the position calculation unit 36 as the performance value for the first time slot on the target day.

第1サンプル群取得部30は、複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群を、記憶部24から取得する。複数の参照日は、予測日より前の予め定められた規則の下で選択される複数の日である。第1サンプル群取得部30は、第1サンプル群を実績値補正部34に与える。 The first sample group acquisition unit 30 acquires from the storage unit 24 a first sample group including a plurality of actual values in a first time period on a plurality of reference dates. The multiple reference dates are multiple days selected under a predetermined rule prior to the forecast date. The first sample group acquisition unit 30 provides the first sample group to the performance value correction unit 34 .

第2サンプル群取得部32は、複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群を、記憶部24から取得する。第2サンプル群取得部32が実績値を取得する複数の参照日は、第1サンプル群取得部30が取得した第1時間帯の実績値と同一である。第2サンプル群取得部32は、第2サンプル群を実績値補正部34に与える。 The second sample group acquisition unit 32 acquires from the storage unit 24 a second sample group including a plurality of actual values for prediction target time periods on a plurality of reference dates. The plurality of reference dates from which the second sample group acquisition unit 32 acquires the actual values are the same as the actual values for the first time period acquired by the first sample group acquisition unit 30 . The second sample group acquisition unit 32 provides the second sample group to the performance value correction unit 34 .

実績値補正部34は、第1サンプル群に含まれる複数の実績値のそれぞれを、予め定められた発電所モデルに基づき算出された対象日計算値と、発電所モデルに基づき算出された参照日計算値との比率により補正する。発電所モデルは、発電所により発電される電力、または、発電所により発電される電力に比例する値を算出するモデルである。また、対象日計算値は、発電所モデルに基づき算出された、対象日における対応する時間帯において発電所により発電される対象日電力または、対象日電力に比例する値である。また、参照日計算値は、発電所モデルに基づき算出された、複数の参照日のうちの対応する参照日における対応する時間帯において発電所により発電される参照日電力または参照日電力に比例する値である。 The actual value correction unit 34 converts each of the plurality of actual values included in the first sample group into a target date calculated value calculated based on a predetermined power plant model and a reference date calculated based on the power plant model. Corrected by the ratio to the calculated value. The power plant model is a model that calculates the power generated by the power plant or a value proportional to the power generated by the power plant. The target day calculated value is the target day power generated by the power plant in the corresponding time period on the target day, or a value proportional to the target day power, which is calculated based on the power plant model. In addition, the reference day calculation value is proportional to the reference day power or the reference day power generated by the power plant in the corresponding time period on the corresponding reference day of the plurality of reference days calculated based on the power plant model. value.

また、実績値補正部34は、第2サンプル群に含まれる複数の実績値のそれぞれを、発電所モデルに基づき算出された対象日計算値と、発電所モデルに基づき算出された参照日計算値との比率により補正する。実績値補正部34は、補正された第1サンプル群を、位置算出部36に与える。また、実績値補正部34は、補正された第2サンプル群を、電力量算出部38に与える。 Further, the actual value correcting unit 34 converts each of the plurality of actual values included in the second sample group into the target date calculated value calculated based on the power plant model and the reference date calculated value calculated based on the power plant model. Corrected by the ratio of The actual value correction unit 34 gives the corrected first sample group to the position calculation unit 36 . In addition, the actual value correction unit 34 supplies the corrected second sample group to the power amount calculation unit 38 .

位置算出部36は、発電所により発電された電力量の複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群に対する、発電所により発電された電力量の対象日における第1時間帯の実績値の対象位置を算出する。本実施形態において、対象位置は、第1サンプル群に含まれる複数の実績値を昇順または降順のうちの何れかの第1順に並べたデータ列内における、対象日における第1時間帯の実績値に対応する電力量の順位である。 The position calculation unit 36 calculates the amount of electric power generated by the power plant on the target day for the first sample group including a plurality of actual values of the amount of electric power generated in the first time period on a plurality of reference days. Calculate the target position of the actual value for one hour. In this embodiment, the target position is the actual value of the first time period on the target date in the data string in which the multiple actual values included in the first sample group are arranged in the first order of either ascending order or descending order. is the order of power consumption corresponding to .

なお、順位は、自然数に限らず、実数値で表されてもよい。また、順位は、0%から100%の間に規格化したパーセンタイル順位であってもよい。位置算出部36は、算出した対象位置を電力量算出部38に与える。 Note that the ranking may be expressed not only by natural numbers but also by real numbers. Also, the rank may be a percentile rank normalized between 0% and 100%. The position calculator 36 provides the calculated target position to the power amount calculator 38 .

電力量算出部38は、発電所により発電された電力量の複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群における、対象位置に対応する電力量を算出する。ここで、対象位置に対応する電力量は、第2サンプル群に含まれる複数の実績値を第1順に並べたデータ列内における、対象位置に対応する電力量である。 The power amount calculation unit 38 calculates the power amount corresponding to the target position in the second sample group including a plurality of actual values of the power amount generated by the power plant in the prediction target time period on a plurality of reference days. Here, the power amount corresponding to the target position is the power amount corresponding to the target position in the data string in which the multiple actual values included in the second sample group are arranged in the first order.

本実施形態においては、電力量算出部38は、対象位置に対応する電力量として、第2サンプル群に含まれる複数の実績値を第1順に並べたデータ列内における、対象位置(順位)に対応する電力量を算出する。電力量算出部38は、算出した電力量を出力部40に与える。 In the present embodiment, the power amount calculation unit 38, as the power amount corresponding to the target position, in the data string in which the plurality of actual values included in the second sample group are arranged in the first order, at the target position (order) Calculate the corresponding power consumption. The power amount calculator 38 provides the calculated power amount to the output unit 40 .

出力部40は、電力量算出部38から受け取った対象位置に対応する電力量を、発電される予測対象時間帯の予測電力量として出力する。 The output unit 40 outputs the power amount corresponding to the target position received from the power amount calculation unit 38 as the predicted power amount for the prediction target time zone to be generated.

図2は、予測装置10の処理の流れを示すフローチャートである。予測装置10の処理の流れを、図2に示すフローチャート、および、図3から図7を参照しながら説明する。 FIG. 2 is a flowchart showing the processing flow of the prediction device 10. As shown in FIG. The flow of processing of the prediction device 10 will be described with reference to the flowchart shown in FIG. 2 and FIGS. 3 to 7. FIG.

まず、記憶部24は、図3に示すように、発電所により発電された電力量の時間帯毎の実績値を記憶している。本実施形態において、記憶部24は、30分毎の電力量を記憶する。より詳しくは、本実施形態において、記憶部24は、各時間における00分00秒から29分59秒までおよび30分00秒から59分59秒までの電力量を記憶する。なお、時間帯は、30分毎に限らず、5分毎、10分毎または60分毎等、どのような時間単位であってもよい。 First, as shown in FIG. 3, the storage unit 24 stores the actual value of the amount of electric power generated by the power plant for each time period. In the present embodiment, the storage unit 24 stores power consumption for every 30 minutes. More specifically, in the present embodiment, the storage unit 24 stores the electric energy from 00 minutes 00 seconds to 29 minutes 59 seconds and from 30 minutes 00 seconds to 59 minutes 59 seconds at each time. Note that the time period is not limited to every 30 minutes, and may be any time unit such as every 5 minutes, 10 minutes, or 60 minutes.

S11において、入力部26は、予測日における予測対象時間帯の入力を受け付ける。予測対象時間帯は、発電所により発電された電力量を予測する対象となる時間帯である。入力部26は、現在時刻より未来の時間帯を、予測日における予測対象時間帯として受け付ける。予測対象時間帯は、実績値を取得する時間単位に対応する。例えば、入力部26は、図4に示すような現在時刻において、図4に示すような予測日における予測対象時間帯を受け付ける。 In S<b>11 , the input unit 26 receives an input of the prediction target time period on the prediction date. The prediction target time period is the time period for which the electric energy generated by the power plant is to be predicted. The input unit 26 accepts a future time period from the current time as a prediction target time period on the prediction date. The prediction target time zone corresponds to the time unit for obtaining the actual value. For example, the input unit 26 receives the prediction target time period on the prediction date as shown in FIG. 4 at the current time as shown in FIG.

続いて、S12において、対象サンプル取得部28は、記憶部24により記憶されている時間帯毎の実績値のうちの、現在時刻より前の時間帯であって、実質的な発電がされた1つの時間帯における実績値を取得する。本実施形態において、対象サンプル取得部28は、記憶部24により記憶されている時間帯毎の実績値のうちの、実質的な発電がされた最新の時間帯における実績値を取得する。 Subsequently, in S12, the target sample acquisition unit 28 selects the actual value for each time period stored in the storage unit 24 for the time period prior to the current time, when the power generation was substantially performed. Get the actual value for one time period. In the present embodiment, the target sample acquiring unit 28 acquires the actual value in the latest time period in which substantial power generation was performed among the actual values for each time period stored in the storage unit 24 .

実質的な発電がされた時間帯とは、実績値が0または所定値以下である時間帯ではなく、実績値が0または所定値より大きい時間帯をいう。例えば、現在時刻が13時10分であり、予測日の12時30分から12時59分に所定値より大きい電力量で発電がされている場合、対象サンプル取得部28は、図5に示すように、予測日の12時30分から12時59分の実績値を取得する。また、例えば、現在時刻が8時00分であり、発電所による発電が17時00分から翌朝の8時59分まで停止されている場合、対象サンプル取得部28は、予測日の前日の16時30分から16時59分の時間帯の実績値を取得する。 The time period during which substantial power generation is performed does not mean the time period in which the actual value is 0 or less than or equal to the predetermined value, but the time period in which the actual value is 0 or greater than the predetermined value. For example, if the current time is 13:10, and power is generated from 12:30 to 12:59 on the predicted date with an amount of power that is greater than the predetermined value, the target sample acquisition unit 28 acquires Then, the actual values from 12:30 to 12:59 on the forecast day are acquired. Further, for example, if the current time is 8:00 and power generation by the power plant is stopped from 17:00 to 8:59 the next morning, the target sample acquisition unit 28 Acquire the actual value for the time period from 30 minutes to 16:59.

続いて、S13において、対象サンプル取得部28は、取得した実績値の時間帯を第1時間帯として特定する。また、対象サンプル取得部28は、取得した実績値の発電日を対象日として特定する。例えば、予測日の当日の実績値を取得した場合、対象サンプル取得部28は、予測日の当日を、対象日として特定する。また、例えば、予測日の前日の実績値を取得した場合、対象サンプル取得部28は、予測日の前日を、対象日として特定する。 Subsequently, in S<b>13 , the target sample acquisition unit 28 specifies the time zone of the acquired performance values as the first time zone. In addition, the target sample acquisition unit 28 identifies the power generation date of the acquired actual value as the target date. For example, when the actual value of the predicted date is obtained, the target sample obtaining unit 28 identifies the predicted date as the target date. Further, for example, when the actual values for the day before the prediction date are obtained, the target sample obtaining unit 28 identifies the day before the prediction date as the target date.

続いて、S14において、第1サンプル群取得部30は、複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群を、記憶部24から取得する。複数の参照日は、予測日より前の予め定められた規則の下で選択される複数の日である。本実施形態において、複数の参照日は、予測日の1日前から30日前までの期間に含まれる30個の日である。例えば、予測日が2021年5月25日である場合、第1サンプル群取得部30は、図6に示すような、2021年4月25日から2021年5月24日までの30個の日を複数の参照日として特定し、これらの日の第1時間帯の実績値を取得する。なお、複数の参照日は、このような期間に限られず、他の規則の下で選択された複数の日であってもよい。 Subsequently, in S<b>14 , the first sample group acquisition unit 30 acquires from the storage unit 24 a first sample group including a plurality of actual values for the first time period on a plurality of reference dates. The multiple reference dates are multiple days selected under a predetermined rule prior to the forecast date. In this embodiment, the multiple reference dates are 30 days included in the period from 1 day to 30 days before the predicted date. For example, if the prediction date is May 25, 2021, the first sample group acquisition unit 30 collects 30 days from April 25, 2021 to May 24, 2021, as shown in FIG. are identified as a plurality of reference dates, and the actual values for the first time period on these days are obtained. Note that the multiple reference dates are not limited to such a period, and may be multiple days selected under other rules.

続いて、S15において、第2サンプル群取得部32は、複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群を、記憶部24から取得する。第2サンプル群取得部32が実績値を取得する複数の参照日は、第1サンプル群取得部30が取得した第1時間帯の実績値と同一である。例えば、第1サンプル群取得部30が2021年4月25日から2021年5月24日までの30個の日の第1時間帯の実績値を取得した場合には、第2サンプル群取得部32は、図7に示すように、2021年4月25日から2021年5月24日までの30個の日の予測対象時間帯の実績値を取得する。 Subsequently, in S<b>15 , the second sample group acquisition unit 32 acquires from the storage unit 24 a second sample group including a plurality of actual values for the prediction target time period on a plurality of reference dates. The plurality of reference dates from which the second sample group acquisition unit 32 acquires the actual values are the same as the actual values for the first time period acquired by the first sample group acquisition unit 30 . For example, when the first sample group acquisition unit 30 acquires actual values for the first time period of 30 days from April 25, 2021 to May 24, 2021, the second sample group acquisition unit As shown in FIG. 7, 32 acquires the actual values of the prediction target time period for 30 days from April 25, 2021 to May 24, 2021.

続いて、S16において、実績値補正部34は、第1サンプル群に含まれる複数の実績値のそれぞれを、予め定められた発電所モデルに基づき算出された対象日計算値と、発電所モデルに基づき算出された参照日計算値との比率により補正する。より具体的には、実績値補正部34は、第1サンプル群に含まれる複数の実績値のそれぞれを、対象日における第1時間帯の対象日計算値と、対応する参照日における第1時間帯の参照日計算値との比率を乗じることにより補正する。 Subsequently, in S16, the performance value correction unit 34 converts each of the plurality of performance values included in the first sample group into the target day calculation value calculated based on the predetermined power plant model and the power plant model. Corrected by the ratio with the reference date calculated value calculated based on More specifically, the actual value correcting unit 34 combines each of the multiple actual values included in the first sample group with the target day calculated value for the first time slot on the target day and the first hour on the corresponding reference date. It is corrected by multiplying the ratio with the reference date calculation value of the band.

例えば、実績値補正部34は、式(1)のように補正をする。 For example, the performance value correcting unit 34 makes corrections as shown in Equation (1).

補正後の第1時間帯の実績値=補正前の第1時間帯の実績値×(対象日における第1時間帯の対象日計算値÷対応する参照日における第1時間帯の参照日計算値)…(1) Actual value for the 1st time period after correction = Actual value for the 1st time period before correction x (Calculated value for the 1st time period on the target date / Calculated value for the 1st time period on the corresponding reference date ) (1)

ここで、発電所モデルは、発電所により発電される電力、または、発電所により発電される電力に比例する値を算出するモデルである。また、対象日計算値は、発電所モデルに基づき算出された、対象日における対応する時間帯において発電所により発電される対象日電力、または、対象日電力に比例する値である。また、参照日計算値は、発電所モデルに基づき算出された、複数の参照日のうちの対応する参照日における対応する時間帯において発電所により発電される参照日電力、または、参照日電力に比例する値である。 Here, the power plant model is a model that calculates the power generated by the power plant or a value proportional to the power generated by the power plant. The target day calculated value is the target day power generated by the power plant in the corresponding time period on the target day, or a value proportional to the target day power, which is calculated based on the power plant model. In addition, the reference day calculated value is the reference day power generated by the power plant in the corresponding time period on the corresponding reference day of the plurality of reference days calculated based on the power plant model, or the reference day power It is a proportional value.

発電所が太陽光発電所である場合、発電所モデルは、一例として、設備情報および気象情報に基づき、対象日または参照日における指定された時刻において発電所により発電される電力を、対象日計算値または参照日計算値として算出する。 If the power plant is a solar power plant, for example, the power plant model calculates the power generated by the power plant at a specified time on the target date or reference date based on facility information and weather information. Calculated as a value or reference date calculation value.

設備情報は、レイアウト角、モジュール容量、および、インバータ容量を含む。レイアウト角は、太陽光発電所に設置される太陽電池モジュールの法線方向である。レイアウト角は、例えば、太陽電池モジュールの方位および傾斜角度により表されてもよい。モジュール容量は、太陽光発電所に設置される太陽電池モジュールの最大出力である。インバータ容量は、太陽光発電所の直流交流変換器の最大出力である。 Equipment information includes layout angle, module capacity, and inverter capacity. The layout angle is the normal direction of the solar cell modules installed in the solar power plant. The layout angle may be represented, for example, by the orientation and tilt angle of the solar cell module. Module capacity is the maximum output of a solar cell module installed in a solar power plant. Inverter capacity is the maximum output of a DC to AC converter in a photovoltaic power plant.

設備情報は、日毎および時間帯毎に変化が無く、対象日計算値と参照日計算値との比率に影響を与えない。従って、発電所モデルは、レイアウト角、モジュール容量、および、インバータ容量を予め定められた規格化された値を用いてもよい。この場合、発電所モデルは、太陽電池モジュールの方位を南向き(0°)とし、太陽電池モジュールの傾斜角度を30°~40°程度としてもよい。また、発電所モデルは、太陽電池モジュールの容量を1.0~2.0程度とし、インバータ容量を1.0程度としてもよい。なお、発電所の設置位置の緯度が高い場合、発電所モデルは、太陽電池モジュールの傾斜角度を30°~40°より大きくしてもよい。 The equipment information does not change for each day and for each time zone, and does not affect the ratio between the target date calculation value and the reference date calculation value. Therefore, the power plant model may use standardized values that predetermine layout angles, module capacities, and inverter capacities. In this case, in the power plant model, the orientation of the solar cell module may be south (0°) and the tilt angle of the solar cell module may be approximately 30° to 40°. In the power plant model, the capacity of the solar cell module may be approximately 1.0 to 2.0, and the inverter capacity may be approximately 1.0. If the latitude of the installation position of the power plant is high, the power plant model may have the inclination angle of the solar cell module larger than 30° to 40°.

気象情報は、例えば、全天水平面日射強度、大気上端水平面日射強度、太陽方向、地表アルベド、気温および風速等を含む。 The weather information includes, for example, all-sky horizontal solar radiation intensity, upper atmospheric horizontal solar radiation intensity, sun direction, ground surface albedo, temperature, wind speed, and the like.

全天水平面日射強度は、水平面に対する日射強度である。大気上端水平面日射強度は、大気上端での水平面の日射強度である。太陽方向は、発電所の設置位置に対する太陽の方向の立体角であり、発電所の設置位置に対する太陽の方位および高度により表される。地表アルベドは、地表面が太陽光を反射する割合である。 The global solar radiation intensity is the solar radiation intensity for the horizontal plane. The horizontal surface solar radiation intensity at the top of the atmosphere is the horizontal surface solar radiation intensity at the top of the atmosphere. The sun direction is the solid angle of the direction of the sun with respect to the installation position of the power plant, represented by the azimuth and altitude of the sun with respect to the installation position of the power plant. The surface albedo is the rate at which the surface of the earth reflects sunlight.

実績値補正部34は、太陽方向に代えて、発電所の所在地の緯度および経度および時刻を取得し、発電所の所在地の緯度および経度および時刻に基づき太陽方向を算出してもよい。また、実績値補正部34は、全天水平面日射強度、大気上端水平面日射強度および地表アルベドを、気象データベースまたは気象提供サービスサーバ等から取得することができる。また、実績値補正部34は、気温および風速を、発電所に設けられた計測機器から取得することができる。また、実績値補正部34は、気温および風速を、気象データベースまたは気象提供サービスサーバから取得してもよい。 The actual value correction unit 34 may acquire the latitude, longitude and time of the location of the power plant instead of the solar direction, and calculate the solar direction based on the latitude, longitude and time of the location of the power plant. Further, the actual value correcting unit 34 can acquire the all-sky horizontal solar radiation intensity, the upper atmospheric horizontal solar radiation intensity, and the surface albedo from a weather database, a weather providing service server, or the like. In addition, the actual value correction unit 34 can acquire the air temperature and the wind speed from measuring equipment provided at the power plant. Also, the actual value correction unit 34 may acquire the temperature and wind speed from a weather database or a weather providing service server.

なお、発電所モデルは、地表アルベドを取得できない場合、地表アルベドとして予め設定された値を用いて、発電所により発電される電力を算出してもよい。また、発電所モデルは、全天水平面日射強度または大気上端水平面日射強度を取得できない場合、全天水平面日射強度または大気上端水平面日射強度として、発電所が設けられた地域の平均的な値等の予め設定された値を用いて、発電所により発電される電力を算出してもよい。また、発電所モデルは、気温または風速を取得できない場合、過去の平均値を用いて、発電所により発電される電力を算出してもよい。 Note that the power plant model may calculate the power generated by the power plant using a preset value for the ground surface albedo when the ground surface albedo cannot be obtained. In addition, when the power plant model cannot obtain the global solar radiation intensity or the upper horizontal solar radiation intensity, the average value of the area where the power plant is installed is used as the global solar radiation intensity or the upper horizontal solar radiation intensity. A preset value may be used to calculate the power generated by the power plant. The power plant model may also calculate the power generated by the power plant using past average values when temperature or wind speed cannot be obtained.

また、太陽光発電所に設けられた太陽電池モジュールに対する日射強度を、傾斜面日射強度という。傾斜面日射強度は、太陽光発電所により発電される電力に対して、ほぼ比例する。従って、発電所モデルは、太陽光発電所により発電される電力に比例する値として、傾斜面日射強度を算出してもよい。この場合、発電所モデルは、対象日または参照日における指定された時刻における傾斜面日射強度を、対象日計算値または参照日計算値として算出する。 Moreover, the solar radiation intensity with respect to the solar cell module provided in the solar power plant is called the inclined plane solar radiation intensity. Inclined solar radiation intensity is approximately proportional to the power generated by the photovoltaic power plant. Therefore, the power plant model may calculate the slope solar radiation intensity as a value proportional to the power generated by the photovoltaic power plant. In this case, the power plant model calculates the slope solar radiation intensity at the specified time on the target date or reference date as the target date calculation value or reference date calculation value.

傾斜面日射強度を算出する場合、太陽電池モジュールは、設備情報、並びに、気象情報のうちの全天水平面日射強度、大気上端水平面日射強度、太陽方向および地表アルベドを用いる。この場合においても、発電所モデルは、地表アルベドを取得できない場合、地表アルベドとして予め設定された値を用いて、傾斜面日射強度を算出してもよい。また、発電所モデルは、全天水平面日射強度または大気上端水平面日射強度を取得できない場合、全天水平面日射強度または大気上端水平面日射強度として、発電所が設けられた地域の平均的な値等の予め設定された値を用いて、傾斜面日射強度を算出してもよい。 When calculating the sloped surface solar radiation intensity, the solar cell module uses facility information and the all-sky horizontal solar radiation intensity, top-of-atmosphere horizontal surface solar radiation intensity, solar direction, and ground surface albedo among the weather information. Even in this case, if the ground surface albedo cannot be acquired, the power plant model may calculate the inclined plane solar radiation intensity using a value preset as the ground surface albedo. In addition, when the power plant model cannot obtain the global solar radiation intensity or the upper horizontal solar radiation intensity, the average value of the area where the power plant is installed is used as the global solar radiation intensity or the upper horizontal solar radiation intensity. A preset value may be used to calculate the slanted surface solar radiation intensity.

なお、発電所モデルの一例については、図10~図18を参照して詳細を後述する。 An example of the power plant model will be described later in detail with reference to FIGS. 10 to 18. FIG.

続いて、S17において、実績値補正部34は、第2サンプル群に含まれる複数の実績値のそれぞれを、発電所モデルに基づき算出された対象日計算値と、発電所モデルに基づき算出された参照日計算値との比率により補正する。より具体的には、実績値補正部34は、第2サンプル群に含まれる複数の実績値のそれぞれを、対象日における予測対象時間帯の対象日計算値と、対応する参照日における予測対象時間帯の参照日計算値との比率を乗じることにより補正する。 Subsequently, in S17, the actual value correction unit 34 sets each of the plurality of actual values included in the second sample group to the target day calculated value calculated based on the power plant model and the calculated value calculated based on the power plant model. Corrected by the ratio with reference date calculation value. More specifically, the actual value correcting unit 34 sets each of the plurality of actual values included in the second sample group to the target day calculated value for the forecast target time period on the target day and the forecast target time on the corresponding reference date. It is corrected by multiplying the ratio with the reference date calculation value of the band.

例えば、実績値補正部34は、式(2)のように補正をする。 For example, the performance value correcting unit 34 makes corrections as shown in Equation (2).

補正後の予測対象時間帯の実績値=補正前の予測対象時間帯の実績値×(対象日における予測対象時間帯の対象日計算値÷対応する参照日における予測対象時間帯の参照日計算値)…(2) Actual value of the forecast target time period after correction = Actual value of the forecast target time period before correction x (calculated value of target day of forecast time period on target day / calculated value of reference date of target time period of corresponding reference date ) (2)

なお、発電所モデルは、S16において第1サンプル群に含まれる複数の実績値を補正したモデルと同一である。 The power plant model is the same as the model obtained by correcting a plurality of actual values included in the first sample group in S16.

続いて、S18において、位置算出部36は、補正された第1サンプル群に対する、発電所により発電された電力量の対象日における第1時間帯の実績値の対象位置を算出する。 Subsequently, in S18, the position calculation unit 36 calculates the target position of the actual value of the amount of power generated by the power plant in the first time period on the target day for the corrected first sample group.

ここで、対象位置は、補正された第1サンプル群に含まれる複数の実績値を昇順または降順のうちの何れかの第1順に並べたデータ列内における、対象日における第1時間帯の実績値に対応する電力量の順位である。本実施形態において、位置算出部36は、対象位置として、第1サンプル群に含まれる複数の実績値を昇順に並べたデータ列内における、対象日における第1時間帯の実績値に対応する電力量の順位を算出する。 Here, the target position is the performance of the first time period on the target day in the data string in which the multiple performance values included in the corrected first sample group are arranged in the first order of either ascending order or descending order. It is the order of the power amount corresponding to the value. In the present embodiment, the position calculation unit 36 calculates, as the target position, the power corresponding to the actual value of the first time period on the target day in the data string in which the multiple actual values included in the first sample group are arranged in ascending order. Calculate quantity ranking.

例えば、第1サンプル群がN個の実績値(Nは2以上の整数)を含む場合、位置算出部36は、式(3)のように、N個の実績値を並べたデータ列を生成する。 For example, when the first sample group includes N actual values (N is an integer equal to or greater than 2), the position calculation unit 36 generates a data string in which the N actual values are arranged as shown in Equation (3). do.

Y[1]≦Y[2]≦Y[3]≦…≦Y[N-2]≦Y[N-1]≦Y[N]…(3) Y[1]≤Y[2]≤Y[3]≤...≤Y[N-2]≤Y[N-1]≤Y[N] (3)

ここで、Y[・]は、第1サンプル群に含まれる実績値を表す。Y[・]のカッコの中は、順位を表す。 Here, Y[·] represents the actual value included in the first sample group. The brackets of Y[·] indicate the order.

位置算出部36は、例えば、このように並べたデータ列と、対象日における第1時間帯の実績値とを比較して、対象位置である順位を算出する。 The position calculation unit 36 compares, for example, the data string arranged in this way with the actual values for the first time slot on the target day, and calculates the order, which is the target position.

一例として、位置算出部36は、データ列内における、対象日における第1時間帯の実績値に最も近い値を検出する。そして、位置算出部36は、検出した値の順位を対象位置としてもよい。また、データ列内に、対象日における第1時間帯の実績値と同一の値が複数個含まれる場合には、位置算出部36は、複数の同一の値の中点の順位を、対象位置としてもよい。この場合、位置算出部36は、実数値により表される順位を、対象位置とする。 As an example, the position calculation unit 36 detects the closest value to the actual value for the first time slot on the target day in the data string. Then, the position calculation unit 36 may use the order of the detected values as the target position. In addition, when the data string contains a plurality of values that are the same as the actual value of the first time period on the target date, the position calculation unit 36 determines the rank of the midpoint of the plurality of the same values as the target position. may be In this case, the position calculation unit 36 uses the order represented by the real number as the target position.

他の一例として、位置算出部36は、データ列に含まれる複数の値のうちの、対象日における第1時間帯の実績値以下の最大の順位の値と、対象日における第1時間帯の実績値以上の最小の順位の値とを取得する。そして、位置算出部36は、取得した2つの値に基づき直線補間して得られる順位を、対象位置としてもよい。 As another example, the position calculation unit 36 calculates, among the plurality of values included in the data string, the highest ranking value that is equal to or lower than the actual value for the first time period on the target day, Get the lowest rank value equal to or greater than the actual value. Then, the position calculation unit 36 may use the order obtained by performing linear interpolation based on the two acquired values as the target position.

具体的には、対象日における第1時間帯の実績値をXとし、第1サンプル群に含まれる対象日における第1時間帯の実績値以下の最大の順位の値をY[n]とし、第1サンプル群に含まれる対象日における第1時間帯の実績値以上の最小の順位の値をY[n+1]とする。この場合、X、Y[n]およびY[n+1]の大小関係は、式(4)のようになる。
Y[n]≦X≦Y[n+1]…(4)
Specifically, let X be the actual value of the first time slot on the target day, let Y[n] be the highest ranking value below the actual value of the first time slot on the target day included in the first sample group, Let Y[n+1] be the lowest ranking value equal to or higher than the actual value for the first time period on the target day included in the first sample group. In this case, the magnitude relationship among X, Y[n] and Y[n+1] is given by Equation (4).
Y[n]≤X≤Y[n+1] (4)

そして、位置算出部36は、式(5)を演算することにより、対象日における第1時間帯の実績値に対応する順位であるnrを算出する。
nr=n×{(Y[n+1]-X)/(Y[n+1]-Y[n])}+(n+1)×{X-(Y[n])/(Y[n+1]-Y[n])}…(5)
Then, the position calculation unit 36 calculates nr, which is the ranking corresponding to the actual value of the first time period on the target day, by calculating the expression (5).
nr=n×{(Y[n+1]−X)/(Y[n+1]−Y[n])}+(n+1)×{X−(Y[n])/(Y[n+1]−Y[n ])}…(5)

なお、順位は、0%から100%の間に規格化したパーセンタイル順位であってもよい。 Note that the rank may be a percentile rank normalized between 0% and 100%.

また、さらに、他の一例として、位置算出部36は、最新の時間帯以前の2以上の所定数の時間帯のそれぞれを第1時間帯とし、2以上の所定数の時間帯のそれぞれについて対象位置を算出してもよい。また、位置算出部36は、過去の予測電力量の予測処理において、既に他の時間帯の対象位置を算出している場合には、過去の過去の予測電力量の予測処理で用いた対象位置を読み出してもよい。そして、位置算出部36は、算出した所定数の対象位置の平均値を算出し、算出した所定数の対象位置の平均値を、対象日における第1時間帯の実績値の対象位置として出力する。これにより、位置算出部36は、直近の時間帯における実績値が急激に変化した場合であっても、精度良く予測電力量を算出させることができる。位置算出部36は、例えば、直近の時間帯の実績値が、直近の時間帯のさらに1つ前の時間帯から所定値以上変動した場合、このような平均値を、対象日における第1時間帯の実績値の対象位置として出力してもよい。 Further, as another example, the position calculation unit 36 sets each of two or more predetermined number of time periods before the latest time period as the first time period, and sets two or more predetermined number of time periods as the first time period. A position may be calculated. In addition, when the position calculation unit 36 has already calculated the target position for another time period in the past prediction processing of the predicted power amount, the position calculation unit 36 may be read. Then, the position calculation unit 36 calculates the average value of the calculated predetermined number of target positions, and outputs the calculated average value of the predetermined number of target positions as the target position of the actual value for the first time period on the target date. . As a result, the position calculation unit 36 can calculate the predicted power amount with high accuracy even when the actual value in the most recent time zone changes abruptly. For example, when the actual value in the most recent time period has changed by a predetermined value or more from the time period immediately before the most recent time period, the position calculation unit 36 calculates such an average value as the first hour of the target day. You may output as a target position of the actual value of a band.

続いて、S19において、電力量算出部38は、補正された第2サンプル群における、対象位置に対応する電力量を算出する。ここで、対象位置に対応する電力量は、第2サンプル群に含まれる複数の実績値を第1順に並べたデータ列内における、対象位置に対応する電力量である。 Subsequently, in S19, the power amount calculator 38 calculates the power amount corresponding to the target position in the corrected second sample group. Here, the power amount corresponding to the target position is the power amount corresponding to the target position in the data string in which the multiple actual values included in the second sample group are arranged in the first order.

本実施形態においては、電力量算出部38は、対象位置に対応する電力量として、第2サンプル群に含まれる複数の実績値を昇順に並べたデータ列内における、対象位置(順位)に対応する電力量を算出する。 In this embodiment, the power amount calculator 38 corresponds to the target position (order) in the data string in which the actual values included in the second sample group are arranged in ascending order as the power amount corresponding to the target position. Calculate the amount of power to be used.

例えば、第2サンプル群がN個の実績値を含む場合、電力量算出部38は、式(6)のように、N個の実績値を並べたデータ列を生成する。 For example, when the second sample group includes N actual values, the power amount calculation unit 38 generates a data string in which the N actual values are arranged, as in Equation (6).

Z[1]≦Z[2]≦Z[3]≦…≦Z[N-2]≦Z[N-1]≦Z[N]…(6) Z[1]≤Z[2]≤Z[3]≤...≤Z[N-2]≤Z[N-1]≤Z[N] (6)

ここで、Z[・]は、第2サンプル群に含まれる実績値を表す。Z[・]のカッコの中は、順位を表す。 Here, Z[·] represents the actual value included in the second sample group. The brackets of Z[·] indicate the order.

電力量算出部38は、例えば、このように並べたデータ列と、対象位置(順位)とを比較して、対象位置(順位)に対応する電力量を算出する。 For example, the power amount calculation unit 38 compares the data string arranged in this way with the target position (ranking), and calculates the power amount corresponding to the target position (ranking).

一例として、電力量算出部38は、データ列内における、対象位置に最も近い順位の実績値を検出する。そして、電力量算出部38は、検出した実績値を対象位置(順位)に対応する電力量としてもよい。 As an example, the power amount calculation unit 38 detects the actual value of the order closest to the target position in the data string. Then, the power amount calculation unit 38 may use the detected actual value as the power amount corresponding to the target position (order).

他の一例として、電力量算出部38は、データ列に含まれる複数の値のうちの、対象位置(順位)以下の最大の順位の実績値と、対象位置(順位)以上の最小の順位の実績値とを取得する。そして、電力量算出部38は、取得した2つの実績値に基づき直線補間して得られる電力量を、対象位置に対応する電力量としてもよい。 As another example, the power amount calculation unit 38 calculates the actual value of the highest ranking below the target position (ranking) and the lowest ranking above the target position (ranking) among the plurality of values included in the data string. Acquire achievement values. Then, the power amount calculation unit 38 may use the power amount obtained by performing linear interpolation based on the two obtained actual values as the power amount corresponding to the target position.

具体的には、対象位置(順位)をnrとし、第2サンプル群に含まれる対象位置(順位)以下の最大の順位の実績値をZ[n]とし、第2サンプル群に含まれる対象位置(順位)以上の最小の順位の実績値をZ[n+1]とする。この場合、nr、nおよび(n+1)の大小関係は、式(7)のようになる。
n≦nr≦(n+1)…(7)
Specifically, let nr be the target position (rank), let Z[n] be the actual value of the highest rank below the target position (rank) included in the second sample group, and let Z[n] be the target position included in the second sample group. Let Z[n+1] be the actual value of the lowest ranking equal to or higher than (ranking). In this case, the magnitude relation among nr, n and (n+1) is given by equation (7).
n≦nr≦(n+1) (7)

そして、電力量算出部38は、式(8)を演算することにより、対象位置(順位)に対応する電力量であるX[nr]を算出する。
X[nr]=Z[n]×{(n+1)-nr}+Z[n+1]×{nr-n}…(8)
Then, the power amount calculation unit 38 calculates X[nr], which is the power amount corresponding to the target position (order), by calculating Expression (8).
X[nr]=Z[n]×{(n+1)−nr}+Z[n+1]×{nr−n} (8)

なお、S19においても、順位は、0%から100%の間に規格化したパーセンタイル順位であってもよい。 Also in S19, the rank may be a percentile rank normalized between 0% and 100%.

続いて、S20において、出力部40は、S19で算出した対象位置(順位)に対応する電力量を、予測日における予測対象時間帯において発電所により発電されると予測される予測電力量として出力する。 Subsequently, in S20, the output unit 40 outputs the electric power amount corresponding to the target position (order) calculated in S19 as the predicted electric power amount that is predicted to be generated by the power plant during the prediction target time period on the prediction date. do.

予測装置10は、以上のS11からS20までの処理を、予測対象時間帯が入力される毎に繰り返す。 The prediction device 10 repeats the above processing from S11 to S20 each time a prediction target time period is input.

なお、S18において、位置算出部36は、対象日における第1時間帯が、現在時刻から予め設定された時間以上前である場合、予め設定された位置(順位)を、対象位置としてもよい。 In S18, the position calculation unit 36 may set a preset position (order) as the target position when the first time period on the target date is earlier than the preset time from the current time.

例えば、発電所の1日における発電期間が9時00分から16時59分までであり、現在時刻が8時00分である場合、対象日における第1時間帯は、前日の16時30分から17時59分となる。このような場合、対象日における第1時間帯は、現在時刻から予め設定された時間以上前となる。また、予め設定された時間以上の通信不良等により、受信部22が実績値を受信することができない場合、対象日における第1時間帯は、現在時刻から予め設定された時間以上前となる。 For example, if the power generation period of the power plant in one day is from 9:00 to 16:59 and the current time is 8:00, the first time slot on the target day is from 16:30 to 17:00 on the previous day. It will be 59 minutes. In such a case, the first time period on the target date is earlier than the preset time from the current time. Also, if the receiving unit 22 cannot receive the actual value due to a communication failure for a preset time or longer, the first time period on the target day is set ahead of the current time by a preset time or longer.

このような場合、対象日における第1時間帯の実績値と、予測対象時間帯の予想電力量との間の相関性が小さくなり、予測装置10は、精度良く予想電力量を算出することができない可能性がある。従って、対象日における第1時間帯が現在時刻から予め設定された時間以上前である場合、予測装置10は、予め設定された位置(順位)を、対象位置とすることにより、電力量の予測精度を高めることができる。 In such a case, the correlation between the actual value for the first time slot on the target day and the expected power consumption for the prediction target time slot decreases, and the prediction device 10 can accurately calculate the expected power consumption. may not be possible. Therefore, when the first time period on the target date is earlier than the preset time from the current time, the prediction device 10 predicts the power amount by setting the preset position (order) as the target position. Accuracy can be improved.

また、S18において、位置算出部36は、対象日における第1時間帯が、現在時刻から予め設定された時間以上前である場合、第1サンプル群に含まれる複数の実績値における最も頻度が高い位置(順位)を、対象位置としてもよい。このようにしても、予測装置10は、精度良く予想電力量を算出することができない可能性がある場合において、電力量の予測精度を高めることができる。 In addition, in S18, if the first time period on the target date is earlier than the preset time from the current time, the position calculation unit 36 determines that the actual values included in the first sample group have the highest frequency. The position (order) may be the target position. Even in this way, the prediction device 10 can improve the prediction accuracy of the power amount when there is a possibility that the predicted power amount cannot be calculated with high accuracy.

図8は、複数の参照日のそれぞれの実績値の変化、および、対象日の実績値の変化の一例を示す図である。 FIG. 8 is a diagram showing an example of changes in performance values for a plurality of reference dates and changes in performance values on a target date.

太陽光発電所により発電された電力量は、日毎に異なるが、例えば直近の30日間等の近い季節での1日における変化の傾向は近似する。従って、例えば直近の30日間等の複数の参照日における、第1時間帯の電力量を昇順に並べた順位と、予測対象時間帯の電力量を昇順に並べた順位とは、近似する。 The amount of power generated by a photovoltaic power plant varies from day to day, but the trend of change during a day in close seasons, such as the last 30 days, is similar. Therefore, the ranking of the power amounts in the first time period arranged in ascending order on a plurality of reference days such as the most recent 30 days approximates the order of the power amounts in the prediction target time period arranged in ascending order.

そこで、予測装置10は、複数の参照日における第1時間帯の複数の実績値に対する、対象日における第1時間帯の実績値の対象位置(順位)を算出し、複数の参照日における予測対象時間帯の複数の実績値における、対象位置(順位)に対応する電力量を算出する。そして、予測装置10は、対象位置に対応する電力量を、予測対象時間帯において発電されると予測される予測電力量として出力する。 Therefore, the prediction device 10 calculates the target position (order) of the actual values of the first time period on the target date with respect to the multiple actual values of the first time period on the multiple reference dates, Calculate the amount of electric power corresponding to the target position (rank) in a plurality of actual values for the time period. Then, the prediction device 10 outputs the power amount corresponding to the target position as the predicted power amount that is predicted to be generated in the prediction target time zone.

このような予測装置10によれば、発電所により発電された電力量の実績値から、予測電力量を算出するので、少ない種類のデータを用いて簡易に精度良く予測電力量を算出することができる。 According to such a prediction device 10, the predicted amount of power is calculated from the actual value of the amount of power generated by the power plant, so it is possible to easily and accurately calculate the predicted amount of power using a small number of types of data. can.

また、予測装置10は、発電所モデルを用いて、複数の参照日における第1時間帯の複数の実績値および複数の参照日における予測対象時間帯の複数の実績値を補正する。これにより、予測装置10は、太陽方向等を含む気象情報に基づき実績値を補正することができるので、より精度の良い予測電力量を出力することができる。 In addition, the prediction device 10 corrects a plurality of actual values for the first time period on a plurality of reference dates and a plurality of actual values for the prediction target time period on a plurality of reference dates using the power plant model. As a result, the prediction device 10 can correct the actual value based on the weather information including the direction of the sun and the like, so that the prediction device 10 can output a more accurate predicted power amount.

(第2実施形態)
つぎに、第2実施形態に係る予測装置10について説明する。なお、第2実施形態の説明では、それまでの実施形態おいて説明した構成要素と略同一の構成および機能を有する構成要素については、同一の符号を付けて、相違点を除き詳細な説明を省略する。
(Second embodiment)
Next, the prediction device 10 according to the second embodiment will be described. In the description of the second embodiment, components having substantially the same configurations and functions as those of the components described in the previous embodiments are denoted by the same reference numerals, and detailed descriptions will be given except for the differences. omitted.

図9は、第2実施形態に係る予測装置10の機能構成を示す図である。第2実施形態に係る予測装置10は、第1実施形態と比較して、位置補正部52をさらに備える。 FIG. 9 is a diagram showing the functional configuration of the prediction device 10 according to the second embodiment. The prediction device 10 according to the second embodiment further includes a position corrector 52 as compared with the first embodiment.

位置補正部52は、気象情報等の付加情報に基づき、位置算出部36により算出された対象位置(順位)を補正して、電力量算出部38に与える。例えば、位置補正部52は、対象日における第1時間帯における位置に対する、予測対象時間帯における位置の予想変化率に基づき、対象位置(順位)を補正する。 The position correction unit 52 corrects the target position (order) calculated by the position calculation unit 36 based on additional information such as weather information, and supplies the corrected target position to the power amount calculation unit 38 . For example, the position correction unit 52 corrects the target position (ranking) based on the predicted change rate of the position in the prediction target time zone with respect to the position in the first time zone on the target day.

位置補正部52は、例えば、対象日における第1時間帯における気象情報と、予測日における予測対象時間帯における予想気象情報とを受け取り、予想変化率を算出してもよい。例えば、位置補正部52は、対象日における第1時間帯は晴天であるが、天候が急激に変化して、予測対象時間帯において雨天となるような場合、対象位置(順位)を低くするように補正する。また、例えば、位置補正部52は、対象日における第1時間帯は雨天であるが、天候が急激に変化して、予測対象時間帯において晴天となるような場合、対象位置(順位)を高くするように補正する。 For example, the position correction unit 52 may receive weather information for the first time slot on the target date and forecast weather information for the forecast target time slot on the forecast date, and calculate the forecast change rate. For example, the position correction unit 52 lowers the target position (rank) when the weather is fine in the first time slot on the target day, but the weather suddenly changes and it becomes rainy in the prediction target time slot. corrected to In addition, for example, the position correction unit 52 increases the target position (rank) when the weather is rainy in the first time period on the target day, but the weather suddenly changes and becomes sunny in the prediction target time period. Correct to

このような第2実施形態に係る予測装置10は、第1時間帯の気象情報と予測対象時間帯との気象情報との違い等により、対象位置(順位)を補正するので、より精度の良い予測電力量を出力することができる。 Since the prediction device 10 according to the second embodiment corrects the target position (order) based on the difference between the weather information of the first time period and the weather information of the prediction target time period, etc., the accuracy is improved. Predicted power consumption can be output.

(発電所モデル)
つぎに、太陽光発電所の発電所モデルについて説明する。
(Power plant model)
Next, the power plant model of the photovoltaic power plant will be explained.

図10は、太陽光発電所の発電所モデルのブロック構成を示す図である。発電所モデルは、設備情報として、レイアウト角、モジュール容量およびインバータ容量を受け取る。さらに、発電所モデルは、気象情報として、対象時刻における、全天水平面日射強度、大気上端水平面日射強度、太陽方向、および、地表アルベドを受け取る。そして、発電所モデルは、対象時刻において、太陽光発電所により発電される電力(kW)を算出する。 FIG. 10 is a diagram showing a block configuration of a power plant model of a photovoltaic power plant. The power plant model receives layout angles, module capacities and inverter capacities as facility information. Furthermore, the power plant model receives, as weather information, the global solar radiation intensity, the top horizontal solar radiation intensity, the direction of the sun, and the ground surface albedo at the target time. Then, the power plant model calculates the power (kW) generated by the solar power plant at the target time.

なお、対象時刻は、電力を推定する時刻である。本実施形態においては、対象時刻は、第1時間帯の中の何れかの時刻、または、予測対象時間帯の中の何れかの時刻となる。例えば、対象時刻は、第1時間帯または予測対象時間帯の中心時刻であってもよい。 Note that the target time is the time at which power is estimated. In the present embodiment, the target time is any time within the first time period or any time within the prediction target time period. For example, the target time may be the center time of the first time period or the prediction target time period.

発電所モデルは、日射変換部112と、温度補正部114と、電力変換部116とを含む。 The power plant model includes a solar conversion unit 112 , a temperature correction unit 114 and a power conversion unit 116 .

日射変換部112は、レイアウト角、全天水平面日射強度、大気上端水平面日射強度、太陽方向および地表アルベドを受け取る。日射変換部112は、これらの情報に基づき、直達散乱分離処理および傾斜面変換処理を実行して、傾斜面日射強度を算出する。 The solar radiation converter 112 receives the layout angle, the global solar radiation intensity, the upper atmospheric horizontal solar radiation intensity, the solar direction, and the ground surface albedo. Based on these pieces of information, the solar radiation conversion unit 112 performs direct scattering separation processing and inclined surface conversion processing to calculate the inclined surface solar radiation intensity.

温度補正部114は、気温、風速、および、日射変換部112により算出された傾斜面日射強度を受け取る。温度補正部114は、これらの情報に基づき、温度補正係数を算出する。 The temperature correction unit 114 receives the air temperature, the wind speed, and the inclined surface solar radiation intensity calculated by the solar radiation conversion unit 112 . The temperature correction unit 114 calculates a temperature correction coefficient based on these pieces of information.

電力変換部116は、モジュール容量、インバータ容量、日射変換部112により算出された傾斜面日射強度、および、温度補正部114により算出された温度補正係数を受け取る。電力変換部116は、これらの情報に基づき、太陽光発電所により発電される電力(kW)を算出する。 The power conversion unit 116 receives the module capacity, the inverter capacity, the inclined surface solar radiation intensity calculated by the solar conversion unit 112 , and the temperature correction coefficient calculated by the temperature correction unit 114 . The power converter 116 calculates the power (kW) generated by the solar power plant based on these pieces of information.

(直達散乱分離処理)
発電所モデルは、太陽電池モジュールの傾斜面に対する日射強度である傾斜面日射強度を取り扱う。気象予測および気象観測で取り扱う日射強度は、水平面に対する日射強度である全天水平面日射強度である。このため、発電所モデルは、全天水平面日射強度を傾斜面日射強度に変換する。
(Direct scattering separation treatment)
The power plant model deals with the sloped surface solar radiation intensity, which is the solar radiation intensity with respect to the sloped surface of the solar cell module. The solar radiation intensity handled in weather prediction and weather observation is the all-sky horizontal solar radiation intensity, which is the solar radiation intensity with respect to the horizontal plane. For this reason, the power plant model converts the global solar radiation intensity into the inclined surface solar radiation intensity.

全天水平面日射強度を傾斜面日射強度に変換するために、まず、日射変換部112は、直散分離モデルを用いて、全天水平面日射強度を、直達光成分と、散乱光成分とに分離する直達散乱分離処理を実行する。全天水平面日射強度をG、直達光成分をG、散乱光成分をGと表した場合、G、GおよびGの関係は、式(11)のように表される。

Figure 2022183783000002
In order to convert the all-sky solar radiation intensity into the inclined-plane solar radiation intensity, first, the solar radiation conversion unit 112 separates the all-sky horizontal solar radiation intensity into a direct light component and a scattered light component using a direct scattering model. Executes the direct scattering separation process. When the global plane solar radiation intensity is G, the direct light component is G b , and the scattered light component is G d , the relationship between G, G b and G d is expressed as in Equation (11).
Figure 2022183783000002

直散分離モデルは、全天水平面日射強度における散乱光成分の割合を算出する。直散分離モデルは、晴天時は散乱光成分が小さく、すなわち直達光成分が大きく、且つ、曇天時は散乱光成分が大きい、すなわち直達光成分が小さい、という定性的な傾向を表現する経験式である。 The direct dispersion model calculates the ratio of the scattered light component in the global solar radiation intensity. The direct scattering separation model is an empirical formula that expresses the qualitative tendency that the scattered light component is small in fine weather, that is, the direct light component is large, and the scattered light component is large in cloudy weather, that is, the direct light component is small. is.

全天水平面日射強度は、雲量以外に、太陽高度に依存する。直散分離モデルは、太陽高度の影響を取り除くため、全天水平面日射強度を、大気上端における水平面日射強度である大気上端水平面日射強度により規格化した快晴指数を用いる。快晴指数をk、大気上端水平面日射強度をGと表した場合、快晴指数は、式(12)のように表される。kは、0以上、1以下の実数となる。

Figure 2022183783000003
The global solar radiation intensity depends on the solar altitude as well as the cloud cover. In order to eliminate the influence of the sun's altitude, the direct-scattering separation model uses a clearness index obtained by normalizing the global solar radiation intensity by the horizontal surface solar radiation intensity at the top of the atmosphere. When the clearness index is represented by k, and the horizontal surface solar radiation intensity at the upper end of the atmosphere is represented by G0 , the clearness index is expressed as in Equation (12). k is a real number of 0 or more and 1 or less.
Figure 2022183783000003

直散分離モデルは、式(13)により表される。

Figure 2022183783000004
A direct dispersion model is represented by Formula (13).
Figure 2022183783000004

(k)は、快晴指数に基づき全天水平面日射強度における散乱光成分の割合を算出する関数であり、kが小さい場合にはkが大きく、kが大きい場合にはkが小さいという傾向を表す関数である。 f d (k) is a function for calculating the ratio of the scattered light component in the all-sky horizontal solar radiation intensity based on the clearness index. When k is small, k d is large, and when k is large, k d is small. It is a function that expresses the tendency of

また、直達光成分は、式(14)により表される。

Figure 2022183783000005
Also, the direct light component is represented by Equation (14).
Figure 2022183783000005

以上から、日射変換部112は、直達散乱分離処理において、式(15-1)、式(15-2)および式(15-3)の演算を実行して、直達光成分および散乱光成分を算出する。

Figure 2022183783000006
From the above, in the direct scattering separation process, the solar radiation conversion unit 112 performs the calculations of formulas (15-1), (15-2), and (15-3) to convert the direct light component and the scattered light component into calculate.
Figure 2022183783000006

図11は、直散分離モデルであるErbs、OrgillおよびReindlの数値例を示す図である。Erbs、OrgillおよびReindlの何れの直散分離モデルも近似した数値を算出する。従って、日射変換部112は、何れの直散分離モデルを用いても近似した計算結果を出力することができる。 FIG. 11 is a diagram showing numerical examples of Erbs, Orgill, and Reindl, which are linear separation models. Any of the Erbs, Orgill and Reindl direct segregation models will give approximate numerical values. Therefore, the solar radiation conversion unit 112 can output an approximated calculation result using any direct-scattering separation model.

図12は、Orgillを用いた場合の直散分離モデルの演算例を示す図である。日射変換部112は、例えば、Orgillの直散分離モデルを用いる場合、図12に示すような演算を実行する。 FIG. 12 is a diagram showing an example of computation of a direct-discrete model using Orgill. For example, when Orgill's direct-scattering separation model is used, the solar radiation conversion unit 112 executes calculations as shown in FIG. 12 .

(傾斜面変換処理)
傾斜面日射強度における直達光成分および散乱光成分の依存度は、太陽電池モジュールのレイアウト角と太陽方向との関係によって変化する。式(16)は、太陽電池モジュールのレイアウト角と太陽方向との角度関係を考慮した、傾斜面日射強度を表す関係式である。

Figure 2022183783000007
(Tilted surface conversion processing)
The degree of dependence of the direct light component and the scattered light component on the slanted surface solar radiation intensity changes depending on the relationship between the layout angle of the solar cell module and the direction of the sun. Equation (16) is a relational expression representing the inclined surface solar radiation intensity, taking into consideration the angular relationship between the layout angle of the solar cell module and the direction of the sun.
Figure 2022183783000007

Sは、傾斜面日射強度を表す。n^は、太陽電池モジュールの法線方向の立体角(レイアウト角)を表す単位ベクトルである。s^は、太陽方向の立体角を表す単位ベクトルである。z^は、鉛直方向の立体角を表す単位ベクトルである。なお、n^、s^およびz^は、数式においてはアルファベットの直上にハット記号が設けられている。また、ρは、地表アルベドを表し、地表面の反射率である。 S represents the slanted surface solar radiation intensity. n̂ is a unit vector representing a solid angle (layout angle) in the normal direction of the solar cell module. s^ is a unit vector representing the solid angle in the direction of the sun. z^ is a unit vector representing a solid angle in the vertical direction. Note that n̂, ŝ, and ẑ are provided with a hat symbol immediately above the alphabet in the mathematical expressions. Also, ρ represents the ground surface albedo, which is the reflectance of the ground surface.

式(16)の右辺の第1項は、太陽電池モジュールの傾斜面に対する直達光成分を表す。式(16)の右辺の第2項は、太陽電池モジュールの傾斜面に対する散乱光成分を表す。式(16)の右辺の第3項は、太陽電池モジュールの傾斜面に対する反射光成分を表す。 The first term on the right side of Equation (16) represents the direct light component with respect to the inclined surface of the solar cell module. The second term on the right side of Equation (16) represents the scattered light component for the inclined surface of the solar cell module. The third term on the right side of Equation (16) represents the reflected light component with respect to the inclined surface of the solar cell module.

傾斜面に対する直達光成分は、指向性のある日射成分であり、太陽電池モジュールの法線方向(n^)が、太陽方向(s^)一致する場合に最大となる。つまり、水平面に対する直達光成分(G)を、鉛直方向に対する太陽方向の立体角(s^,z^)で規格化した値(G/(s^,z^))が直達光成分の最大値となる。従って、式(16)の右辺の第1項は、直達光成分の最大値に、太陽電池モジュールの法線方向に対する太陽方向の立体角(s^,n^)を乗算した値となる。 The direct light component to the inclined surface is a directional solar radiation component, and becomes maximum when the normal direction (n̂) of the solar cell module coincides with the sun direction (ŝ). That is, the value (G b /(s^, z^)) obtained by normalizing the direct light component (G b ) with respect to the horizontal plane by the solid angle (s^, z^) in the direction of the sun with respect to the vertical direction is the direct light component. maximum value. Therefore, the first term on the right side of equation (16) is a value obtained by multiplying the maximum value of the direct light component by the solid angle (ŝ, n̂) in the direction of the sun with respect to the normal direction of the solar cell module.

傾斜面に対する散乱光成分は、指向性のない日射成分であり、天球の各方向(r^)から同じ強度で太陽電池モジュールに入射する成分(r^,n^)を、立体角平均した立体角に比例する。従って、傾斜面に対する散乱光成分は、式(17)に表される値に比例する。

Figure 2022183783000008
The scattered light component on the inclined surface is a solar radiation component without directivity, and is a solid angle average of the components (r^, n^) incident on the solar cell module with the same intensity from each direction (r^) of the celestial sphere. proportional to the angle. Therefore, the scattered light component for the inclined plane is proportional to the value expressed in equation (17).
Figure 2022183783000008

なお、式(17)において、<・>r^は、天球の各方向(r^)により立体角平均した値を表す。従って、式(16)の右辺の第2項は、式(17)の右辺にGを乗算した値となる。 In equation (17), <·> represents a solid angle average value for each direction (r̂) of the celestial sphere. Therefore, the second term on the right side of equation (16) is the value obtained by multiplying the right side of equation (17) by Gd .

傾斜面に対する反射光成分は、指向性のない日射成分であり、地表からの反射光をρGの散乱光の強度の散乱光と解釈すれば散乱光成分と同様である。従って、式(16)の右辺の第3項は、式(17)の右辺にρGを乗算した値となる。 The reflected light component on the inclined surface is a solar radiation component with no directivity, and is similar to the scattered light component if the reflected light from the ground surface is interpreted as the scattered light with the intensity of the scattered light of ρGd . Therefore, the third term on the right side of equation (16) is the value obtained by multiplying the right side of equation (17) by ρGd .

以上から、日射変換部112は、太陽光発電所の緯度および経度、時刻およびレイアウト角に基づき、n^、s^およびz^を算出する。そして、日射変換部112は、式(16)の関係式に基づき、傾斜面日射強度を算出する。 From the above, the solar radiation converter 112 calculates n^, s^, and z^ based on the latitude and longitude, time, and layout angle of the solar power plant. Then, the solar radiation conversion unit 112 calculates the inclined surface solar radiation intensity based on the relational expression (16).

図13は、太陽方向の数値例を示す図である。図14は、図13に示すように太陽方向が変化する場合における傾斜面日射強度の数値例を示す図である。 FIG. 13 is a diagram showing numerical examples of the sun direction. FIG. 14 is a diagram showing numerical examples of the inclined surface solar radiation intensity when the direction of the sun changes as shown in FIG.

図13および図14の例においては、kを0.8とし、発電所の緯度および経度を35.00N,135.00Eとしている。また、図13および図14の例においては、南から西までの方位を0°から90°までで表し、水平から天頂までの高度を0°から90°までで表している。図14においては、太陽電池モジュールの傾斜角度を0°、30°および60°で変化させ、太陽電池モジュールの方位を30°としている。従って、図14の例においては、傾斜角度が大きい場合、早朝に直達光成分が届かず、傾斜面日射強度は、小さくなる。 In the examples of FIGS. 13 and 14, k is 0.8 and the latitude and longitude of the power station are 35.00N and 135.00E. In the examples of FIGS. 13 and 14, the azimuth from south to west is represented by 0° to 90°, and the altitude from the horizontal to the zenith is represented by 0° to 90°. In FIG. 14, the tilt angle of the solar cell module is changed between 0°, 30° and 60°, and the azimuth of the solar cell module is 30°. Therefore, in the example of FIG. 14, when the inclination angle is large, the direct light component does not reach early in the morning, and the inclined surface solar radiation intensity becomes small.

図15は、太陽方向の演算例を示す図である。太陽方向は、緯度、経度および時刻に依存する関数により算出することが可能である。発電所モデルは、太陽光発電所の緯度、経度および対象時刻に基づき、図15に示されるような関数を用いて、太陽方向を算出してもよい。 FIG. 15 is a diagram showing an example of calculation of the direction of the sun. Sun direction can be calculated by a function that depends on latitude, longitude and time. The power plant model may calculate the direction of the sun using a function as shown in FIG. 15 based on the latitude, longitude and target time of the solar power plant.

(温度補正係数の算出)
太陽光発電所の効率は、主に、太陽電池モジュールの変換効率に依存する。一般に、太陽電池モジュールは、モジュール温度が低いほど変換効率が高くなる。モジュール温度は、周囲の気温、風速および傾斜面日射強度の影響を受ける。モジュール温度は、傾斜面日射強度が弱いと周囲の気温に近くなるが、傾斜面日射強度が強いと周囲の気温より高くなる。また、モジュール温度は、風速が大きいと日射による温度上昇が抑えられる。
(Calculation of temperature correction coefficient)
The efficiency of a photovoltaic power plant mainly depends on the conversion efficiency of the photovoltaic modules. In general, the lower the module temperature, the higher the conversion efficiency of a solar cell module. Module temperature is affected by ambient air temperature, wind speed, and sloping surface insolation intensity. The module temperature becomes close to the surrounding air temperature when the slanting surface solar radiation intensity is weak, but becomes higher than the surrounding air temperature when the slanting surface solar radiation intensity is strong. Also, the temperature rise of the module due to solar radiation can be suppressed when the wind speed is high.

傾斜面日射強度をS、傾斜面日射強度の吸収率をa、熱伝達効率をh、周囲の気温をT、モジュール温度をTとした場合、太陽電池モジュールから逃げる熱流束であるqは、式(18-1)のように表される。従って、モジュール温度であるTは、式(18-2)のように表される。

Figure 2022183783000009
If S is the solar radiation intensity on the inclined surface, a is the absorption rate of the solar radiation intensity on the inclined surface, h is the heat transfer efficiency, T is the ambient air temperature, and Tm is the module temperature, the heat flux q that escapes from the solar cell module is It is expressed as in formula (18-1). Therefore, Tm , which is the module temperature, is expressed as in Equation (18-2).
Figure 2022183783000009

熱伝達効率であるhは、一般に風速に応じて大きくなるが、流れの状態および設置条件に依存するので、容易に計算することができない。そこで、温度補正部114は、モジュール温度を、傾斜面日射強度、気温および風速を含む式(19)に示す経験式に基づき、算出してもよい。なお、式(19)において、Vは、風速である。また、式(19)において、AおよびBは、定数である。例えば、式(19)において、Aは約50であり、Bは約0.4である。

Figure 2022183783000010
The heat transfer efficiency, h, generally increases with wind speed, but cannot be easily calculated because it depends on flow conditions and installation conditions. Therefore, the temperature correction unit 114 may calculate the module temperature based on the empirical formula shown in formula (19) including the sloping surface solar radiation intensity, the air temperature, and the wind speed. In addition, in Formula (19), V is a wind speed. Also, in equation (19), A and B are constants. For example, in equation (19), A is approximately 50 and B is approximately 0.4.
Figure 2022183783000010

図16は、式(19)を用いて算出されたモジュール温度と周囲の気温との差を表す図である。温度補正部114は、例えば、図16に示す数値に基づき、傾斜面日射強度および風速から、モジュール温度を算出する。 FIG. 16 is a diagram showing the difference between the module temperature calculated using equation (19) and the ambient air temperature. The temperature correction unit 114 calculates the module temperature from the inclined surface solar radiation intensity and the wind speed based on the numerical values shown in FIG. 16, for example.

また、太陽電池モジュールの出力電流-出力電圧特性は、主に、短絡電流と開放電圧とにより特徴付けられる。短絡電流は、傾斜面日射強度に応じて増加する。しかし、開放電圧は、モジュール温度の上昇に伴い減少する。太陽電池モジュールの変換効率は、発生電力を傾斜面日射強度で割り算した値であり、温度上昇に伴い減少する。太陽電池モジュールの温度補正係数は、モジュール温度が25.0℃で規格化された変換効率である。 Also, the output current-output voltage characteristics of a solar cell module are mainly characterized by short-circuit current and open-circuit voltage. The short-circuit current increases according to the slanted surface solar radiation intensity. However, the open circuit voltage decreases as the module temperature increases. The conversion efficiency of a solar cell module is a value obtained by dividing the generated electric power by the slanted surface solar radiation intensity, and decreases as the temperature rises. The temperature correction coefficient of the solar cell module is the conversion efficiency normalized at a module temperature of 25.0°C.

温度補正係数とモジュール温度との関係は、温度係数を用いた経験式により近似することができる。温度補正係数をKPVとした場合、温度補正係数は、式(20)により表される。なお、αは、-0.4/100程度の値となる。

Figure 2022183783000011
The relationship between the temperature correction coefficient and the module temperature can be approximated by an empirical formula using the temperature coefficient. When the temperature correction coefficient is KPV , the temperature correction coefficient is represented by Equation (20). Note that α is a value of about -0.4/100.
Figure 2022183783000011

以上から、温度補正部114は、モジュール温度を算出した後、式(20)に基づき温度補正係数を算出する。例えば、温度補正係数は、モジュール温度が55℃の場合、発電効率が12%程度低下させる数値となる。 From the above, after calculating the module temperature, the temperature correction unit 114 calculates the temperature correction coefficient based on Equation (20). For example, when the module temperature is 55° C., the temperature correction coefficient is a numerical value that reduces the power generation efficiency by about 12%.

(電力の算出処理)
太陽光発電所の主要パラメータは、モジュール容量およびインバータ容量ある。インバータ容量に対するモジュール容量の比率は、過積載率と呼ばれており、稼働率を高めるために、1.0~2.0程度の値に設定される。この場合、太陽電池モジュールの直流出力電力は、式(21-1)により表され、インバータの交流出力電力は、式(21-2)により表される。

Figure 2022183783000012
(Power calculation process)
The main parameters of a photovoltaic power plant are module capacity and inverter capacity. A ratio of the module capacity to the inverter capacity is called an overload rate, and is set to a value of about 1.0 to 2.0 in order to increase the operating rate. In this case, the DC output power of the solar cell module is expressed by equation (21-1), and the AC output power of the inverter is expressed by equation (21-2).
Figure 2022183783000012

DCは、太陽電池モジュールの直流出力電力である。YACは、インバータの交流出力電力である。CPVは、モジュール容量である。CPCSは、インバータ容量である。KPVは、インバータ変換効率である。ηは、過積載率である。min(x,y)は、xおよびyのうち小さい方の値を選択する関数である。 Y DC is the DC output power of the solar module. Y AC is the AC output power of the inverter. C PV is the module capacitance. C PCS is the inverter capacity. K PV is the inverter conversion efficiency. η is the overload rate. min(x,y) is a function that selects the smaller of x and y.

なお、モジュール容量は、通常、モジュール温度が25℃および傾斜面日射強度が1.0kW/mの標準条件の下の出力電力と定義されている。従って、太陽電池モジュールの直流出力電力は、傾斜面日射強度に比例する簡単な計算式で表される。 It should be noted that the module capacity is usually defined as the output power under standard conditions of module temperature of 25° C. and slant surface irradiance of 1.0 kW/m 2 . Therefore, the DC output power of the solar cell module can be expressed by a simple formula proportional to the slanted surface solar radiation intensity.

太陽光発電所の最終出力電力は、インバータの交流出力電力である。インバータの交流出力電力は、太陽電池モジュールの直流出力電力を、インバータ容量で制限した値となる。 The final output power of the photovoltaic power plant is the AC output power of the inverter. The AC output power of the inverter is a value obtained by limiting the DC output power of the solar cell module by the inverter capacity.

図17は、負荷率に対する、ンバータの変換効率を表す図である。発電所モデルは、インバータの変換効率を、近似的に定数である1.0とみなしてもよい。しかし、発電所モデルは、経験的パラメータを含む関数を用いて、インバータ変換効率を算出してもよい。例えば、発電所モデルは、式(22-1)および式(22-2)の関数を用いて、インバータ変換効率を算出してもよい。

Figure 2022183783000013
FIG. 17 is a diagram showing the conversion efficiency of the inverter with respect to the load factor. The power plant model may consider the conversion efficiency of the inverter to be approximately a constant 1.0. However, the power plant model may use a function containing empirical parameters to calculate the inverter conversion efficiency. For example, the power plant model may calculate the inverter conversion efficiency using the functions of equations (22-1) and (22-2).
Figure 2022183783000013

xは、負荷率を表す。式(22-2)において、A、BおよびCは、経験的パラメータを表す。Aは、0.05程度である。Bは、1.0程度である。Cは、0.02程度である。式(22-2)の関数は、負荷率が0.63(x=√(C/A))の場合に、最大値が0/97(KPCS(x)=1/(B+2√AC))となる。発電所モデルは、式(22-2)の関数を用いてインバータ変換効率を算出することにより、低負荷時における電力の推定精度を良くすることができる。 x represents the load factor. In equation (22-2), A, B and C represent empirical parameters. A is about 0.05. B is about 1.0. C is about 0.02. The function of formula (22-2) has a maximum value of 0/97 (K PCS (x)=1/(B+2√AC) when the load factor is 0.63 (x=√(C/A)). ). The power plant model can improve the accuracy of estimating power at low load by calculating the inverter conversion efficiency using the function of equation (22-2).

(発電所モデルの処理フロー)
図18は、発電所モデルの処理の流れを示すフローチャートである。以上の説明をまとめると、発電所モデルは、図18に示す流れで処理を実行する。
(Processing flow of power plant model)
FIG. 18 is a flow chart showing the process flow of the power plant model. To summarize the above description, the power plant model executes processing according to the flow shown in FIG.

まず、S51において、発電所モデルは、以下の式(23)に示す設備情報を取得する。

Figure 2022183783000014
First, in S51, the power plant model acquires equipment information shown in the following equation (23).
Figure 2022183783000014

なお、発電所モデルは、レイアウト角、モジュール容量、および、インバータ容量を予め定められた規格化された値を取得してもよい。 The power plant model may obtain standardized values that are predetermined for the layout angle, module capacity, and inverter capacity.

続いて、S52において、発電所モデルは、以下の式(24)に示す気象情報を取得する。

Figure 2022183783000015
Subsequently, in S52, the power plant model acquires the weather information shown in Equation (24) below.
Figure 2022183783000015

なお、発電所モデルは、地表アルベドを取得できない場合、地表アルベドとして予め設定された値を取得してもよい。また、発電所モデルは、全天水平面日射強度または大気上端水平面日射強度を取得できない場合、全天水平面日射強度または大気上端水平面日射強度として、発電所が設けられた地域の平均的な値等の予め設定された値を取得してもよい。また、発電所モデルは、気温または風速を取得できない場合、過去の平均値を用いて、発電所により発電される電力を算出してもよい。 Note that the power plant model may acquire a preset value as the ground surface albedo when the ground surface albedo cannot be acquired. In addition, when the power plant model cannot obtain the global solar radiation intensity or the upper horizontal solar radiation intensity, the average value of the area where the power plant is installed is used as the global solar radiation intensity or the upper horizontal solar radiation intensity. A preset value may be acquired. The power plant model may also calculate the power generated by the power plant using past average values when temperature or wind speed cannot be obtained.

続いて、S53において、発電所モデルは、以下の式(25-1)、式(25-2)、式(25-3)および式(25-4)を演算して、直達散乱分離処理を実行する。

Figure 2022183783000016
Subsequently, in S53, the power plant model calculates the following equations (25-1), (25-2), (25-3) and (25-4) to perform direct scattering separation processing. Run.
Figure 2022183783000016

続いて、S54において、発電所モデルは、以下の式(26-1)、式(26-2)、式(26-3)および式(26-4)を演算して、傾斜面変換処理を実行する。これにより、発電所モデルは、傾斜面日射強度(S)を算出することができる。

Figure 2022183783000017
Subsequently, in S54, the power plant model calculates the following equations (26-1), (26-2), (26-3) and (26-4) to perform the slope conversion process. Run. Thereby, the power plant model can calculate the inclined surface solar radiation intensity (S).
Figure 2022183783000017

続いて、S55において、発電所モデルは、以下の式(27-1)および式(27-2)を演算して、温度補正係数を算出する。

Figure 2022183783000018
Subsequently, in S55, the power plant model calculates the temperature correction coefficient by calculating the following equations (27-1) and (27-2).
Figure 2022183783000018

続いて、S56において、発電所モデルは、以下の式(28-1)および式(28-2)を演算して、太陽光発電所から出力される電力を算出する。

Figure 2022183783000019
Subsequently, in S56, the power plant model computes the following equations (28-1) and (28-2) to calculate the power output from the solar power plant.
Figure 2022183783000019

そして、S57において、発電所モデルは、式(29)に示す電力を、対象日計算値または参照日計算値として出力する。

Figure 2022183783000020
Then, in S57, the power plant model outputs the power shown in Equation (29) as the target day calculated value or the reference day calculated value.
Figure 2022183783000020

なお、発電所モデルは、S51からS54までの処理を実行し、傾斜面日射強度(S)を、対象日計算値または参照日計算値として出力してもよい。この場合、発電所モデルは、S52において気温および風速を取得しなくてよい。 Note that the power plant model may execute the processes from S51 to S54 and output the inclined plane solar radiation intensity (S) as the calculated value for the target day or the calculated value for the reference day. In this case, the power plant model does not need to acquire temperature and wind speed in S52.

以上のように、発電所モデルは、第1サンプル群および第2サンプル群に含まれる複数の実績値を補正するための計算値を算出することができる。 As described above, the power plant model can calculate calculated values for correcting a plurality of actual values included in the first sample group and the second sample group.

(ハードウェア構成)
図19は、実施形態に係る予測装置10のハードウェア構成の一例を示す図である。本実施形態に係る予測装置10は、例えば図19に示すようなハードウェア構成の情報処理装置により実現される。予測装置10は、CPU(Central Processing Unit)201と、RAM(Random Access Memory)202と、ROM(Read Only Memory)203と、操作入力装置204と、記憶装置206と、通信装置207とを備える。そして、これらの各部は、バスにより接続される。
(Hardware configuration)
FIG. 19 is a diagram illustrating an example of the hardware configuration of the prediction device 10 according to the embodiment. The prediction device 10 according to this embodiment is realized by, for example, an information processing device having a hardware configuration as shown in FIG. 19 . The prediction device 10 includes a CPU (Central Processing Unit) 201 , a RAM (Random Access Memory) 202 , a ROM (Read Only Memory) 203 , an operation input device 204 , a storage device 206 and a communication device 207 . These units are connected by a bus.

CPU201は、プログラムに従って演算処理および制御処理等を実行するプロセッサである。CPU201は、RAM202の所定領域を作業領域として、ROM203および記憶装置206等に記憶されたプログラムとの協働により各種処理を実行する。 The CPU 201 is a processor that executes arithmetic processing, control processing, and the like according to programs. The CPU 201 uses a predetermined area of the RAM 202 as a work area and executes various processes in cooperation with programs stored in the ROM 203 and the storage device 206 or the like.

RAM202は、SDRAM(Synchronous Dynamic Random Access Memory)等のメモリである。RAM202は、CPU201の作業領域として機能する。ROM203は、プログラムおよび各種情報を書き換え不可能に記憶するメモリである。 A RAM 202 is a memory such as SDRAM (Synchronous Dynamic Random Access Memory). A RAM 202 functions as a work area for the CPU 201 . The ROM 203 is a memory that non-rewritably stores programs and various information.

操作入力装置204は、マウスおよびキーボード等の入力デバイスである。操作入力装置204は、ユーザから操作入力された情報を指示信号として受け付け、指示信号をCPU201に出力する。 The operation input device 204 is an input device such as a mouse and keyboard. The operation input device 204 receives information input by a user as an instruction signal, and outputs the instruction signal to the CPU 201 .

記憶装置206は、フラッシュメモリ等の半導体による記憶媒体、または、磁気的若しくは光学的に記録可能な記憶媒体等にデータを書き込みおよび読み出しをする装置である。記憶装置206は、CPU201からの制御に応じて、記憶媒体にデータの書き込みおよび読み出しをする。通信装置207は、CPU201からの制御に応じて外部の機器とネットワークを介して通信する。 The storage device 206 is a device that writes data to and reads data from a semiconductor storage medium such as a flash memory, or a magnetically or optically recordable storage medium. The storage device 206 writes data to and reads data from a storage medium under the control of the CPU 201 . A communication device 207 communicates with an external device via a network under the control of the CPU 201 .

本実施形態の予測装置10で実行されるプログラムは、受信モジュールと、入力モジュールと、対象サンプル取得モジュールと、第1サンプル群取得モジュールと、第2サンプル群取得モジュールと、実績値補正モジュールと、位置算出モジュールと、電力量算出モジュールと、出力モジュールとを備える。さらに、プログラムは、位置補正モジュールを含んでもよい。このプログラムは、CPU201(プロセッサ)によりRAM202上に展開して実行されることにより、情報処理装置を受信部22、入力部26、対象サンプル取得部28、第1サンプル群取得部30、第2サンプル群取得部32、実績値補正部34、位置算出部36、電力量算出部38、および、出力部40として機能させる。さらに、このプログラムは、情報処理装置を位置補正部52として機能させてもよい。なお、予測装置10は、受信部22、入力部26、対象サンプル取得部28、第1サンプル群取得部30、第2サンプル群取得部32、実績値補正部34、位置算出部36、電力量算出部38、出力部40および位置補正部52の少なくとも一部をハードウェア回路(例えば半導体集積回路)により実現した構成であってもよい。 The program executed by the prediction device 10 of the present embodiment includes a reception module, an input module, a target sample acquisition module, a first sample group acquisition module, a second sample group acquisition module, a performance value correction module, It comprises a position calculation module, a power calculation module and an output module. Additionally, the program may include a position correction module. This program is expanded on the RAM 202 and executed by the CPU 201 (processor), thereby operating the information processing apparatus as a receiving unit 22, an input unit 26, a target sample acquiring unit 28, a first sample group acquiring unit 30, a second sample, and so on. It functions as a group acquisition unit 32 , a performance value correction unit 34 , a position calculation unit 36 , a power amount calculation unit 38 , and an output unit 40 . Furthermore, this program may cause the information processing device to function as the position corrector 52 . The prediction device 10 includes a receiving unit 22, an input unit 26, a target sample acquiring unit 28, a first sample group acquiring unit 30, a second sample group acquiring unit 32, a performance value correcting unit 34, a position calculating unit 36, an electric energy At least part of the calculation unit 38, the output unit 40, and the position correction unit 52 may be realized by a hardware circuit (for example, a semiconductor integrated circuit).

また、本実施形態の予測装置10で実行されるプログラムは、コンピュータにインストール可能な形式または実行可能な形式のファイルで、CD-ROM、フレキシブルディスク、CD-R、DVD(Digital Versatile Disk)等のコンピュータで読み取り可能な記録媒体に記録されて提供される。 In addition, the program executed by the prediction device 10 of the present embodiment is a file in a format that can be installed in a computer or in a format that can be executed, and is stored on a CD-ROM, a flexible disk, a CD-R, a DVD (Digital Versatile Disk), etc. It is recorded on a computer-readable recording medium and provided.

また、本実施形態の予測装置10で実行されるプログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。また、本実施形態の予測装置10で実行されるプログラムをインターネット等のネットワーク経由で提供または配布するように構成してもよい。また、予測装置10で実行されるプログラムを、ROM203等に予め組み込んで提供するように構成してもよい。 Alternatively, the program executed by the prediction device 10 of this embodiment may be stored in a computer connected to a network such as the Internet, and may be provided by being downloaded via the network. Also, the program executed by the prediction device 10 of this embodiment may be configured to be provided or distributed via a network such as the Internet. Alternatively, the program to be executed by the prediction device 10 may be configured so as to be pre-installed in the ROM 203 or the like and provided.

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。 While several embodiments of the invention have been described, these embodiments have been presented by way of example and are not intended to limit the scope of the invention. These novel embodiments can be implemented in various other forms, and various omissions, replacements, and modifications can be made without departing from the scope of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are included in the scope of the invention described in the claims and equivalents thereof.

10 予測装置
22 受信部
24 記憶部
26 入力部
28 対象サンプル取得部
30 第1サンプル群取得部
32 第2サンプル群取得部
34 実績値補正部
36 位置算出部
38 電力量算出部
40 出力部
52 位置補正部
10 prediction device 22 reception unit 24 storage unit 26 input unit 28 target sample acquisition unit 30 first sample group acquisition unit 32 second sample group acquisition unit 34 performance value correction unit 36 position calculation unit 38 power amount calculation unit 40 output unit 52 position corrector

Sokol Dervishi and Ardeshir Mahdavi, "COMPARISON OF MODELS FOR THE DERIVATION OF DIFFUSE FRACTION OF GLOBAL IRRADIANCE DATA FOR VIENNA, AUSTRIA", Proceedings of Building Simulation 2011: 12th Conference of International Building Performance Simulation Association, Sydney, 14-16 November,P.765-771Sokol Dervishi and Ardeshir Mahdavi, "COMPARISON OF MODELS FOR THE DERIVATION OF DIFFUSE FRACTION OF GLOBAL IRRADIANCE DATA FOR VIENNA, AUSTRIA", Proceedings of Building Simulation 2011: 12th Conference of International Building Performance Simulation Association, Sydney, 14-16 November,P. 765-771

Claims (20)

1または複数の発電所により発電された電力量の複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群に対する、前記1または複数の発電所により発電された電力量の対象日における前記第1時間帯の実績値の対象位置を算出する位置算出部と、
前記1または複数の発電所により発電された電力量の前記複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群における、前記対象位置に対応する電力量を算出する電力量算出部と、
前記対象位置に対応する電力量を、前記予測対象時間帯において前記1または複数の発電所により発電されると予測される予測電力量として出力する出力部と、
を備える情報処理装置。
subjecting the amount of electricity generated by the one or more power plants to a first sample set comprising a plurality of actual values for a first time period on multiple reference days of the amount of electricity generated by the one or more power plants. a position calculation unit that calculates a target position for the actual value of the first time period in a day;
Electric energy for calculating the electric energy corresponding to the target position in a second sample group including a plurality of actual values of the electric energy generated by the one or more power plants in the prediction target time period on the plurality of reference days. a calculation unit;
an output unit that outputs the amount of power corresponding to the target position as a predicted amount of power that is predicted to be generated by the one or more power plants in the prediction target time period;
Information processing device.
前記対象位置は、前記第1サンプル群に含まれる複数の実績値を昇順または降順のうちの何れかの第1順に並べたデータ列内における、前記対象日における前記第1時間帯の実績値に対応する電力量の順位であり、
前記対象位置に対応する電力量は、前記第2サンプル群に含まれる複数の実績値を前記第1順に並べたデータ列内における、前記対象位置に対応する電力量である
請求項1に記載の情報処理装置。
The target position is the actual value of the first time period on the target date in a data string in which the multiple actual values included in the first sample group are arranged in the first order of either ascending order or descending order. is the order of the corresponding electric energy,
The power amount corresponding to the target position is the power amount corresponding to the target position in a data string in which a plurality of actual values included in the second sample group are arranged in the first order. Information processing equipment.
前記第1サンプル群に含まれる複数の実績値のそれぞれおよび前記第2サンプル群に含まれる複数の実績値のそれぞれを、予め定められた発電所モデルに基づき算出された対象日計算値と、前記発電所モデルに基づき算出された参照日計算値との比率により補正する実績値補正部をさらに備え、
前記発電所モデルは、前記1または複数の発電所により発電される電力、または、前記1または複数の発電所により発電される電力に比例する値を算出するモデルであり、
前記対象日計算値は、前記発電所モデルに基づき算出された、前記対象日における対応する時間帯において前記1または複数の発電所により発電される対象日電力または前記対象日電力に比例する値であり、
前記参照日計算値は、前記発電所モデルに基づき算出された、前記複数の参照日のうちの対応する参照日における対応する時間帯において前記1または複数の発電所により発電される参照日電力または前記参照日電力に比例する値である
請求項1または2に記載の情報処理装置。
Each of the plurality of actual values included in the first sample group and each of the plurality of actual values included in the second sample group are calculated based on a predetermined power plant model, and the Further includes an actual value correction unit that corrects according to the ratio with the reference date calculation value calculated based on the power plant model,
The power plant model is a model that calculates the power generated by the one or more power plants or a value proportional to the power generated by the one or more power plants,
The target day calculated value is the target day power generated by the one or more power plants in the corresponding time period on the target day or a value proportional to the target day power calculated based on the power plant model. can be,
The reference day calculated value is the reference day power generated by the one or more power plants in the corresponding time period on the corresponding reference date of the plurality of reference days, calculated based on the power plant model, or The information processing apparatus according to claim 1 or 2, wherein the value is proportional to the reference daily power.
前記実績値補正部は、
前記第1サンプル群に含まれる複数の実績値のそれぞれを、前記対象日における前記第1時間帯の前記対象日計算値と、対応する前記参照日における前記第1時間帯の前記参照日計算値との比率を乗じることにより補正し、
前記第2サンプル群に含まれる複数の実績値のそれぞれを、前記対象日における前記第1時間帯の前記対象日計算値と、対応する前記参照日における前記第1時間帯の前記参照日計算値との比率を乗じることにより補正する
請求項3に記載の情報処理装置。
The actual value correction unit,
Each of the plurality of actual values included in the first sample group is calculated as the target date calculated value for the first time period on the target date and the reference date calculated value for the first time period on the corresponding reference date. corrected by multiplying by the ratio of
Each of the plurality of actual values included in the second sample group is calculated as the target date calculated value for the first time period on the target date and the reference date calculated value for the first time period on the corresponding reference date. The information processing apparatus according to claim 3, wherein the correction is performed by multiplying a ratio of .
前記1または複数の発電所のそれぞれは、同一の前記発電所モデルを用いて発電する電力を推定可能な太陽光発電所であり、
前記発電所モデルは、太陽電池モジュールの法線方向を示すレイアウト角、全天水平面日射強度、大気上端水平面日射強度、太陽方向および地表アルベドに基づき、前記太陽電池モジュールの傾斜面に対する日射強度である傾斜面日射強度を算出する
請求項4に記載の情報処理装置。
Each of the one or more power plants is a solar power plant capable of estimating the power generated using the same power plant model,
The power plant model is the solar radiation intensity with respect to the inclined surface of the solar cell module based on the layout angle indicating the normal direction of the solar cell module, the solar radiation intensity on the all-sky horizontal plane, the solar radiation intensity on the top horizontal plane of the atmosphere, the solar direction, and the ground surface albedo. The information processing apparatus according to claim 4, wherein the inclined plane solar radiation intensity is calculated.
前記発電所モデルは、前記地表アルベドを取得できない場合、前記地表アルベドとして予め設定された値を用いて前記傾斜面日射強度を算出する
請求項5に記載の情報処理装置。
The information processing apparatus according to claim 5, wherein, when the ground surface albedo cannot be acquired, the power plant model calculates the inclined plane solar radiation intensity using a value preset as the ground surface albedo.
前記発電所モデルは、前記全天水平面日射強度または前記大気上端水平面日射強度を取得できない場合、前記全天水平面日射強度または前記大気上端水平面日射強度として予め設定された値を用いて前記傾斜面日射強度を算出する
請求項5または6に記載の情報処理装置。
The power plant model, when the global solar radiation intensity or the upper atmosphere horizontal surface solar radiation intensity cannot be obtained, uses a preset value as the global solar radiation intensity or the upper atmosphere horizontal surface solar radiation intensity to calculate the inclined surface solar radiation intensity. The information processing apparatus according to claim 5 or 6, wherein the intensity is calculated.
前記1または複数の発電所のそれぞれは、同一の前記発電所モデルを用いて発電する電力を推定可能な太陽光発電所であり、
前記発電所モデルは、太陽電池モジュールの法線方向を示すレイアウト角、全天水平面日射強度、大気上端水平面日射強度、太陽方向、地表アルベド、気温および風速に基づき、前記1または複数の発電所により発電される電力を算出する
請求項4に記載の情報処理装置。
Each of the one or more power plants is a solar power plant capable of estimating the power generated using the same power plant model,
The power plant model is based on the layout angle indicating the normal direction of the solar cell module, the all-sky solar radiation intensity, the top horizontal solar radiation intensity, the solar direction, the ground surface albedo, the temperature and the wind speed, and the one or more power plants The information processing apparatus according to claim 4, wherein the power to be generated is calculated.
前記発電所モデルは、前記地表アルベドを取得できない場合、前記地表アルベドとして予め設定された値を用いて、前記1または複数の発電所により発電される電力を算出する
請求項8に記載の情報処理装置。
9. The information processing according to claim 8, wherein when the ground surface albedo cannot be obtained, the power plant model calculates the electric power generated by the one or more power plants using a preset value as the ground surface albedo. Device.
前記発電所モデルは、前記全天水平面日射強度または前記大気上端水平面日射強度を取得できない場合、前記全天水平面日射強度または前記大気上端水平面日射強度として予め設定された値を用いて前記1または複数の発電所により発電される電力を算出する
請求項8または9に記載の情報処理装置。
When the global solar radiation intensity or the upper atmospheric horizontal solar radiation intensity cannot be obtained, the power plant model uses a preset value as the global horizontal solar radiation intensity or the upper atmospheric horizontal solar radiation intensity to generate the one or more The information processing apparatus according to claim 8 or 9, which calculates the electric power generated by the power plant.
前記発電所モデルは、前記気温または前記風速を取得できない場合、過去の平均値を用いて前記1または複数の発電所により発電される電力を算出する
請求項8から10の何れか1項に記載の情報処理装置。
11. The power plant model according to any one of claims 8 to 10, wherein, when the air temperature or the wind speed cannot be obtained, the power generated by the one or more power plants is calculated using past average values. information processing equipment.
前記第1時間帯は、前記1または複数の発電所により電力が発電された最新の時間帯であり、
前記予測対象時間帯は、前記対象日における前記第1時間帯よりも後の時間帯である
請求項2から11の何れか1項に記載の情報処理装置。
The first time period is the most recent time period in which power was generated by the one or more power plants;
The information processing apparatus according to any one of claims 2 to 11, wherein the prediction target time period is a time period later than the first time period on the target date.
前記位置算出部は、前記対象日における前記第1時間帯が、現在時刻から予め設定された時間以上前である場合、予め設定された値を前記対象位置とする
請求項12に記載の情報処理装置。
13. The information processing according to claim 12, wherein, when the first time period on the target date is earlier than a preset time from the current time, the position calculation unit sets a preset value as the target position. Device.
前記対象日における前記第1時間帯が、現在時刻から予め設定された時間以上前である場合、前記位置算出部は、前記第1サンプル群に含まれる複数の実績値における最も頻度が高い位置を、前記対象位置とする
請求項12に記載の情報処理装置。
When the first time period on the target date is earlier than the preset time from the current time, the position calculation unit calculates the position with the highest frequency among the plurality of actual values included in the first sample group. , is the target position.
前記位置算出部は、前記1または複数の発電所により電力が発電された最新の時間帯以前の2以上の所定数の時間帯のそれぞれを前記第1時間帯として、所定数の前記対象位置を算出し、算出した所定数の前記対象位置の平均値を、前記対象日における前記第1時間帯の実績値の対象位置として出力する
請求項2から11の何れか1項に記載の情報処理装置。
The position calculation unit calculates a predetermined number of target positions, each of a predetermined number of two or more time periods before the latest time period in which power was generated by the one or more power plants as the first time period. 12. The information processing apparatus according to any one of claims 2 to 11, wherein the average value of the calculated predetermined number of target positions is output as the target position of the actual value of the first time period on the target date. .
前記対象日における前記第1時間帯における位置に対する、前記予測対象時間帯における位置の予想変化率に基づき、前記対象位置を補正する位置補正部をさらに備える
請求項2から15の何れか1項に記載の情報処理装置。
16. The method according to any one of claims 2 to 15, further comprising a position correction unit that corrects the target position based on a predicted change rate of the position in the prediction target time zone with respect to the position in the first time zone on the target date. The information processing device described.
前記位置補正部は、前記対象日における前記第1時間帯における気象情報と、予測日における前記予測対象時間帯における予想気象情報とに基づき、前記予想変化率を算出する
請求項16に記載の情報処理装置。
17. The information according to claim 16, wherein the position correction unit calculates the predicted rate of change based on weather information for the first time slot on the target date and predicted weather information for the prediction target time slot on the prediction date. processing equipment.
情報処理装置により実行される情報処理方法であって、
前記情報処理装置が、1または複数の発電所により発電された電力量の複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群に対する、前記1または複数の発電所により発電された電力量の対象日における前記第1時間帯の実績値の対象位置を算出し、
前記情報処理装置が、記1または複数の発電所により発電された電力量の前記複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群における、前記対象位置に対応する電力量を算出し、
前記情報処理装置が、前記対象位置に対応する電力量を、前記予測対象時間帯において前記1または複数の発電所により発電されると予測される予測電力量として出力する
情報処理方法。
An information processing method executed by an information processing device,
The information processing device performs power generation by the one or more power plants for a first sample group including a plurality of actual values of the amount of electric power generated by the one or more power plants in a first time period on a plurality of reference days. Calculate the target position of the actual value of the first time period on the target day of the electric energy,
The power corresponding to the target position in a second sample group including a plurality of actual values of the amount of power generated by the one or more power plants in the prediction target time period on the plurality of reference dates. calculate the amount of
The information processing method, wherein the information processing device outputs the electric energy corresponding to the target position as a predicted electric energy expected to be generated by the one or more power plants in the prediction target time period.
情報処理装置を、電力量を予測する予測装置として機能させるためのプログラムであって、
前記情報処理装置を、
1または複数の発電所により発電された電力量の複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群に対する、前記1または複数の発電所により発電された電力量の対象日における前記第1時間帯の実績値の対象位置を算出する位置算出部と、
前記1または複数の発電所により発電された電力量の前記複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群における、前記対象位置に対応する電力量を算出する電力量算出部と、
前記対象位置に対応する電力量を、前記予測対象時間帯において前記1または複数の発電所により発電されると予測される予測電力量として出力する出力部と、
して機能させるプログラム。
A program for causing an information processing device to function as a prediction device that predicts power consumption,
the information processing device,
subjecting the amount of electricity generated by the one or more power plants to a first sample set comprising a plurality of actual values for a first time period on multiple reference days of the amount of electricity generated by the one or more power plants. a position calculation unit that calculates a target position for the actual value of the first time period in a day;
Electric energy for calculating the electric energy corresponding to the target position in a second sample group including a plurality of actual values of the electric energy generated by the one or more power plants in the prediction target time period on the plurality of reference days. a calculation unit;
an output unit that outputs the amount of power corresponding to the target position as a predicted amount of power that is predicted to be generated by the one or more power plants in the prediction target time period;
A program that works as
1または複数の発電所と、
前記1または複数の発電所により発電される電力量を予測する情報処理装置と、
を備え、
前記情報処理装置は、
1または複数の発電所により発電された電力量の複数の参照日における第1時間帯の複数の実績値を含む第1サンプル群に対する、前記1または複数の発電所により発電された電力量の対象日における前記第1時間帯の実績値の対象位置を算出する位置算出部と、
前記1または複数の発電所により発電された電力量の前記複数の参照日における予測対象時間帯の複数の実績値を含む第2サンプル群における、前記対象位置に対応する電力量を算出する電力量算出部と、
前記対象位置に対応する電力量を、前記予測対象時間帯において前記1または複数の発電所により発電されると予測される予測電力量として出力する出力部と、
を有する発電システム。
one or more power plants;
an information processing device that predicts the amount of power generated by the one or more power plants;
with
The information processing device is
subjecting the amount of electricity generated by the one or more power plants to a first sample set comprising a plurality of actual values for a first time period on multiple reference days of the amount of electricity generated by the one or more power plants. a position calculation unit that calculates a target position for the actual value of the first time period in a day;
Electric energy for calculating the electric energy corresponding to the target position in a second sample group including a plurality of actual values of the electric energy generated by the one or more power plants in the prediction target time period on the plurality of reference days. a calculation unit;
an output unit that outputs the amount of power corresponding to the target position as a predicted amount of power that is predicted to be generated by the one or more power plants in the prediction target time period;
A power generation system with
JP2021091268A 2021-05-31 2021-05-31 Information processing equipment, information processing methods, programs and power generation systems Active JP7043661B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021091268A JP7043661B1 (en) 2021-05-31 2021-05-31 Information processing equipment, information processing methods, programs and power generation systems

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021091268A JP7043661B1 (en) 2021-05-31 2021-05-31 Information processing equipment, information processing methods, programs and power generation systems

Publications (2)

Publication Number Publication Date
JP7043661B1 JP7043661B1 (en) 2022-03-29
JP2022183783A true JP2022183783A (en) 2022-12-13

Family

ID=81215099

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021091268A Active JP7043661B1 (en) 2021-05-31 2021-05-31 Information processing equipment, information processing methods, programs and power generation systems

Country Status (1)

Country Link
JP (1) JP7043661B1 (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011142790A (en) * 2010-01-08 2011-07-21 Fuji Electric Co Ltd Prediction system of solar power generation amount and prediction method of solar power generation amount
JP2021039739A (en) * 2019-08-28 2021-03-11 富士電機株式会社 Prediction device, prediction program, and prediction method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011142790A (en) * 2010-01-08 2011-07-21 Fuji Electric Co Ltd Prediction system of solar power generation amount and prediction method of solar power generation amount
JP2021039739A (en) * 2019-08-28 2021-03-11 富士電機株式会社 Prediction device, prediction program, and prediction method

Also Published As

Publication number Publication date
JP7043661B1 (en) 2022-03-29

Similar Documents

Publication Publication Date Title
Zhang et al. A solar time based analog ensemble method for regional solar power forecasting
Mayer et al. Techno-economic optimization of grid-connected, ground-mounted photovoltaic power plants by genetic algorithm based on a comprehensive mathematical model
Massidda et al. Use of multilinear adaptive regression splines and numerical weather prediction to forecast the power output of a PV plant in Borkum, Germany
Draxl et al. The wind integration national dataset (wind) toolkit
Carneiro et al. Review on photovoltaic power and solar resource forecasting: current status and trends
Abuella et al. Solar power forecasting using support vector regression
Constantinescu et al. A computational framework for uncertainty quantification and stochastic optimization in unit commitment with wind power generation
Junior et al. Regional forecasts and smoothing effect of photovoltaic power generation in Japan: An approach with principal component analysis
Kato Prediction of photovoltaic power generation output and network operation
Mayer et al. Pairing ensemble numerical weather prediction with ensemble physical model chain for probabilistic photovoltaic power forecasting
Lorenz et al. Overview of irradiance and photovoltaic power prediction
CN110472762A (en) A kind of correction method and system of gridding precipitation
Alomari et al. A predictive model for solar photovoltaic power using the Levenberg-Marquardt and Bayesian regularization algorithms and real-time weather data
Camargo et al. Electricity self-sufficiency of single-family houses in Germany and the Czech Republic
US10255393B2 (en) Optimally placing photovoltaic arrays to maximize value of energy production based on peak power production, local solar radiation, weather, electricity market prices and rate structures
CN115049171B (en) Photovoltaic power prediction method and system based on feature migration
CN107204615A (en) A kind of method and system for realizing power prediction
CN113256033B (en) Power station site selection and volume fixing method, device, system and medium based on mode complementation
Trigo-González et al. Photovoltaic power electricity generation nowcasting combining sky camera images and learning supervised algorithms in the Southern Spain
Haupt Short-range forecasting for energy
Eyring et al. High-resolution electricity generation model demonstrates suitability of high-altitude floating solar power
CN116260141B (en) Reconstruction method and system for power of photovoltaic power station and reconstruction terminal
Holmgren et al. A comparison of pv power forecasts using pvlib-python
Pierro et al. Machine learning-based PV power forecasting methods for electrical grid management and energy trading
JP2022183783A (en) Information processing device, information processing method, program, and power generation system

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211119

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211119

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20211119

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: 20220215

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220316

R150 Certificate of patent or registration of utility model

Ref document number: 7043661

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150