JP2018018329A - Power amount prediction program, power amount prediction device and power amount prediction method - Google Patents

Power amount prediction program, power amount prediction device and power amount prediction method Download PDF

Info

Publication number
JP2018018329A
JP2018018329A JP2016148561A JP2016148561A JP2018018329A JP 2018018329 A JP2018018329 A JP 2018018329A JP 2016148561 A JP2016148561 A JP 2016148561A JP 2016148561 A JP2016148561 A JP 2016148561A JP 2018018329 A JP2018018329 A JP 2018018329A
Authority
JP
Japan
Prior art keywords
prediction
solar radiation
amount
radiation amount
power generation
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
JP2016148561A
Other languages
Japanese (ja)
Other versions
JP6772625B2 (en
Inventor
博之 樋口
Hiroyuki Higuchi
博之 樋口
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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2016148561A priority Critical patent/JP6772625B2/en
Publication of JP2018018329A publication Critical patent/JP2018018329A/en
Application granted granted Critical
Publication of JP6772625B2 publication Critical patent/JP6772625B2/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
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A30/00Adapting or protecting infrastructure or their operation
    • 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)

Abstract

PROBLEM TO BE SOLVED: To determine a prediction error of a power generation amount.SOLUTION: A power amount prediction device 100 generates a boundary value 143 on the basis of a prediction amount of solar radiation in weather prediction/actual past data 141 and an actual amount of solar radiation. The power amount prediction device 100 estimates if a prediction error of the amount of solar radiation swings to an increase direction or a decrease direction by comparing a clearness index of an object date with the boundary value 143. The power amount prediction device 100 specifies a prediction error of the amount of solar radiation, which satisfies a guarantee probability in the estimated direction, and determines a prediction error of a power generation amount.SELECTED DRAWING: Figure 1

Description

本発明は、電力量予測プログラム等に関する。   The present invention relates to a power amount prediction program and the like.

太陽光発電(PV:Photovoltaics)やバッテリによる電力供給が可能となった結果、PVの発電量およびバッテリの放電量に応じて、電力供給に関わる電力コストを制御することが、需要者・供給者の双方で求められている。   As a result of the fact that photovoltaic power generation (PV: Photovoltaics) and battery power supply became possible, it is possible to control the power cost related to power supply according to the amount of PV power generation and the amount of battery discharge. Both are required.

ここで、天候により日射量が変化し、日射量の変化によりPVの発電量が変化するため、電力の供給の予測を行う際は、天候を考慮した予測を行う。例えば、従来技術として、過去の気象情報に基づき予測される日射量(あるいは日射量の予測誤差)の発生確率が正規分布に従うと仮定し、所定の保障確率を満たすような日射量を特定することで、発電量を見積もる技術がある。   Here, the amount of solar radiation changes depending on the weather, and the amount of PV power generated changes depending on the amount of solar radiation. Therefore, when the power supply is predicted, the prediction is performed in consideration of the weather. For example, as a conventional technique, assuming that the probability of occurrence of solar radiation predicted based on past weather information (or prediction error of solar radiation) follows a normal distribution, specify the amount of solar radiation that satisfies a predetermined guarantee probability There is a technology for estimating the amount of power generation.

図12は、従来技術を説明するための図である。図12の横軸は日射量を示し、縦軸は発生確率を示す。例えば、保障確率を90%とすると、従来技術では、日射量をxと仮定して、PVの発電量を見積もる。 FIG. 12 is a diagram for explaining the prior art. The horizontal axis in FIG. 12 indicates the amount of solar radiation, and the vertical axis indicates the occurrence probability. For example, if the guarantee probability to 90%, in the prior art, the amount of solar radiation on the assumption that x 0, estimate the power generation amount of PV.

特開2006−304402号公報JP 2006-304402 A

しかしながら、上述した従来技術では、発電量の予測誤差を求めることができないという問題がある。   However, the above-described conventional technology has a problem that it is impossible to obtain a power generation amount prediction error.

上述した従来技術では、日射量の発生確率の分布として、正規分布を仮定しているが、実際には、日射量の発生確率の分布は、正規分布とはならない。このため、同様に、日射量の予測誤差の分布も、正規分布とはならない。   In the above-described prior art, a normal distribution is assumed as the distribution of the probability of occurrence of solar radiation, but in reality, the distribution of the probability of occurrence of solar radiation is not a normal distribution. For this reason, similarly, the distribution of the prediction error of the solar radiation amount is not a normal distribution.

図13は、従来技術の問題を説明するための図である。図13において、横軸は晴天指数を示し、縦軸は発生確率を示す。晴天指数は、水平面日射量を大気外日射量で割ることで得られる値である。図13の丸印は、日照なしの場合の晴天指数と発生確率との関係を示すものであり、曇りや雨の場合における晴天指数と発生確率との関係に対応する。ダイヤ印は、日照ありの場合の晴天指数と発生確率との関係を示すものであり、晴れの場合における晴天指数と発生確率との関係に対応する。丸印の分布もダイヤ印の分布も左右対称の正規分布とは大きく異なる。四角印は、日照なしの場合および日照ありの場合を合わせた全体の晴天指数と発生確率との関係を示す。   FIG. 13 is a diagram for explaining a problem of the conventional technique. In FIG. 13, the horizontal axis indicates the clear sky index, and the vertical axis indicates the occurrence probability. The clear sky index is a value obtained by dividing the horizontal solar radiation amount by the atmospheric solar radiation amount. The circles in FIG. 13 indicate the relationship between the clear sky index and the occurrence probability when there is no sunshine, and correspond to the relationship between the clear sky index and the occurrence probability when it is cloudy or rainy. The diamond mark indicates the relationship between the clear sky index and the occurrence probability when there is sunshine, and corresponds to the relationship between the clear sky index and the occurrence probability when the sky is clear. The distribution of circles and the distribution of diamonds are very different from the symmetrical normal distribution. The square marks indicate the relationship between the overall clear sky index and the probability of occurrence when there is no sunshine and when there is sunshine.

図13に示すように、曇りや雨の場合における晴天指数は、0付近の発生確率が大きくなり、晴れの場合における晴天指数は、大気外日射量付近の発生確率が大きくなる。   As shown in FIG. 13, the clear sky index in the case of cloudy or rain has a high probability of occurrence near 0, and the clear sky index in the case of clear weather has a high probability of occurrence near the solar radiation amount outside the atmosphere.

ここで、曇りや雨の予想が外れ、晴れとなった場合には、想定していた日射量が増加することはあっても、減少することは少ないため、曇りや雨の日射量の予測誤差は、減少方向の振れ幅がほぼ0であるのに対して、増加方向の振れ幅は一定量ある。   Here, if the forecast of cloudy or rain is off and clear, the expected amount of solar radiation may increase, but it will rarely decrease. In contrast, the fluctuation width in the decreasing direction is almost zero, while the fluctuation width in the increasing direction is a certain amount.

一方、晴れの予想が外れ、曇りや雨になった場合には、想定していた日射量が減少することはあっても、増加することはない。このため、晴れの日射量の予測誤差は、増加方向の振れ幅がほぼ0であるのに対して、減少方向の振れ幅は一定量ある。   On the other hand, if the forecast of clear weather is not met and it becomes cloudy or rainy, the expected amount of solar radiation may decrease but not increase. For this reason, the prediction error of the sunny solar radiation amount has a constant amount of fluctuation in the decreasing direction, while the fluctuation amount in the increasing direction is almost zero.

従来技術では、日射量の予測誤差の分布を正規分布と仮定している。このため、例えば、晴れの日射量を予測する場合に、考慮しなくても良い増加方向の予測誤差を考慮してしまう場合や、曇りや雨の日射量を予測する場合に、考慮しなくてもよい減少方向の予測誤差を考慮してしまう場合があった。   In the prior art, the distribution of solar radiation prediction errors is assumed to be a normal distribution. For this reason, for example, when predicting sunny solar radiation, it is not necessary to consider when forecasting errors in an increasing direction that do not need to be considered, or when predicting cloudy or rainy solar radiation. In some cases, the prediction error in the decreasing direction was taken into consideration.

1つの側面では、本発明は、発電量の予測誤差を求めることができる電力量予測プログラム、電力量予測装置および電力量予測方法を提供することを目的とする。   In one aspect, an object of the present invention is to provide a power amount prediction program, a power amount prediction apparatus, and a power amount prediction method capable of obtaining a power generation amount prediction error.

第1の案では、コンピュータに下記の処理を実行させる。コンピュータは、予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、次の処理を行う。コンピュータは、実測日射量に対する予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成する。コンピュータは、予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する。   In the first plan, the computer executes the following processing. The computer performs the following processing based on the predicted solar radiation amount predicted at the prediction target point and date and the actually measured solar radiation amount at the prediction target point and date and time. The computer generates prediction error information that specifies whether the prediction error of the predicted solar radiation amount with respect to the actually measured solar radiation amount is in an increasing direction or a decreasing direction for each weather prediction. The computer calculates a power generation amount generated by the solar power generation device and a power generation amount prediction error based on the prediction error information, the solar radiation amount prediction information, and the solar power generation device information.

発電量の予測誤差を求めることができる。   The prediction error of the power generation amount can be obtained.

図1は、本実施例に係る電力量予測装置の構成を示す機能ブロック図である。FIG. 1 is a functional block diagram illustrating a configuration of the power amount prediction apparatus according to the present embodiment. 図2は、天候予測・実績過去データのデータ構造の一例を示す図である。FIG. 2 is a diagram illustrating an example of a data structure of weather prediction / actual past data. 図3は、天候予測データのデータ構造の一例を示す図である。FIG. 3 is a diagram illustrating an example of a data structure of weather prediction data. 図4は、第1度数分布データの一例を示す図である。FIG. 4 is a diagram illustrating an example of the first frequency distribution data. 図5は、第2度数分布データの一例を示す図である。FIG. 5 is a diagram illustrating an example of the second frequency distribution data. 図6は、限界日射量算出部の処理を説明するための図(1)である。FIG. 6 is a diagram (1) for explaining the process of the limit solar radiation amount calculation unit. 図7は、限界日射量算出部の処理を説明するための図(2)である。FIG. 7 is a diagram (2) for explaining the processing of the limit solar radiation amount calculation unit. 図8は、充放電計画最適化部の処理を説明するための図である。FIG. 8 is a diagram for explaining the processing of the charge / discharge plan optimization unit. 図9は、電力量予測装置が境界値および度数分布データを生成する処理手順を示すフローチャートである。FIG. 9 is a flowchart illustrating a processing procedure in which the power amount prediction apparatus generates boundary values and frequency distribution data. 図10は、電力量予測装置が充放電計画を算出する処理手順を示すフローチャートである。FIG. 10 is a flowchart illustrating a processing procedure in which the power amount prediction apparatus calculates a charge / discharge plan. 図11は、電力量予測プログラムを実行するハードウェア構成の一例を示す図である。FIG. 11 is a diagram illustrating an example of a hardware configuration that executes the power amount prediction program. 図12は、従来技術を説明するための図である。FIG. 12 is a diagram for explaining the prior art. 図13は、従来技術の問題を説明するための図である。FIG. 13 is a diagram for explaining a problem of the conventional technique.

以下に、本願の開示する電力量予測プログラム、電力量予測装置および電力量予測方法の実施例を図面に基づいて詳細に説明する。なお、この実施例によりこの発明が限定されるものではない。   Hereinafter, embodiments of a power amount prediction program, a power amount prediction apparatus, and a power amount prediction method disclosed in the present application will be described in detail with reference to the drawings. Note that the present invention is not limited to the embodiments.

図1は、本実施例に係る電力量予測装置の構成を示す機能ブロック図である。図1に示すように、この電力量予測装置100は、通信部110と、入力部120と、表示部130と、記憶部140と、制御部150とを有する。   FIG. 1 is a functional block diagram illustrating a configuration of the power amount prediction apparatus according to the present embodiment. As illustrated in FIG. 1, the power amount prediction apparatus 100 includes a communication unit 110, an input unit 120, a display unit 130, a storage unit 140, and a control unit 150.

通信部110は、ネットワークを介して、外部のサーバ等とデータ通信を実行する処理部である。例えば、電力量予測装置100は、通信部110を用いて、外部のサーバと通信を行い、後述する天候予測・実績過去データ141、天候予測データ144等を取得しても良い。   The communication unit 110 is a processing unit that performs data communication with an external server or the like via a network. For example, the power amount prediction apparatus 100 may communicate with an external server using the communication unit 110 to acquire weather prediction / actual past data 141, weather prediction data 144, and the like, which will be described later.

入力部120は、各種の情報を電力量予測装置100に入力するための入力装置である。入力部120は、キーボードやマウス、タッチパネル等に対応する。   The input unit 120 is an input device for inputting various types of information to the power amount prediction device 100. The input unit 120 corresponds to a keyboard, a mouse, a touch panel, or the like.

表示部130は、制御部150から出力される情報を表示する表示装置である。表示部130は、液晶ディスプレイやタッチパネル等に対応する。   The display unit 130 is a display device that displays information output from the control unit 150. The display unit 130 corresponds to a liquid crystal display, a touch panel, or the like.

記憶部140は、天候予測・実績過去データ141と、度数分布データ142と、境界値143と、天候予測データ144と、発電機器データ145を有する。記憶部140は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ(Flash Memory)などの半導体メモリ素子、またはハードディスク、光ディスクなどの記憶装置に対応する。   The storage unit 140 includes weather prediction / result past data 141, frequency distribution data 142, boundary values 143, weather prediction data 144, and power generation equipment data 145. The storage unit 140 corresponds to, for example, a semiconductor memory device such as a random access memory (RAM), a read only memory (ROM), or a flash memory, or a storage device such as a hard disk or an optical disk.

天候予測・実績過去データ141は、過去に予測した過去の日時に対する天候予測および過去の天候の実績に関する情報を保持するテーブルである。図2は、天候予測・実績過去データのデータ構造の一例を示す図である。図2に示すように、この天候予測・実績過去データ141は、日時と、経度・緯度と、天気と、予測日射量と、実測日射量とを対応付ける。   The weather prediction / actual past data 141 is a table that holds information on the weather forecast for the past date and time predicted in the past and the past weather results. FIG. 2 is a diagram illustrating an example of a data structure of weather prediction / actual past data. As shown in FIG. 2, the weather prediction / actual past data 141 associates date / time, longitude / latitude, weather, predicted solar radiation amount, and actual solar radiation amount.

天候予測・実績過去データ141の日時は、現在日時よりも過去の日時を示す。経度・緯度は、日射量の予測の対象となった位置、日射量の観測を行う位置に対応する。天気は、例えば、晴れ、曇り、雨のいずれかを示すものである。予測日射量は、該当する日時、経度・緯度において、予測された日射量である。実測日射量は、該当する日時、経度・緯度において、実際に測定された日射量の実測値である。   The date and time of the weather prediction / result past data 141 indicates a date and time that is earlier than the current date and time. The longitude and latitude correspond to the position where the solar radiation amount is predicted and the position where the solar radiation amount is observed. The weather indicates, for example, sunny, cloudy or rainy. The predicted amount of solar radiation is the amount of solar radiation predicted at the relevant date and time, longitude and latitude. The actually measured amount of solar radiation is an actually measured value of the amount of solar radiation actually measured at the corresponding date and time, longitude and latitude.

度数分布データ142は、晴天指数誤差と頻度との関係を示す情報である。度数分布データ142に関する具体的な説明は後述する。   The frequency distribution data 142 is information indicating the relationship between the clear sky index error and the frequency. A specific description regarding the frequency distribution data 142 will be described later.

境界値143は、予測誤差の方向が日射量の増加する方向になる晴天指数と、予測誤差の方向が日射量の減少する方向になる晴天誤差との境界となる晴天指数を示すものである。   The boundary value 143 indicates a clear sky index that is a boundary between a clear sky index in which the direction of the prediction error is a direction in which the amount of solar radiation increases and a clear sky error in which the direction of the prediction error is a direction in which the amount of solar radiation decreases.

天候予測データ144は、現在日時よりも先の日時において、予測される天候に関する情報を保持するテーブルである。図3は、天候予測データのデータ構造の一例を示す図である。図3に示すように、この天候予測データ144は、日時と、経度・緯度と、天気と、予測日射量とを対応付ける。   The weather prediction data 144 is a table that holds information related to the weather that is predicted at a date earlier than the current date. FIG. 3 is a diagram illustrating an example of a data structure of weather prediction data. As shown in FIG. 3, the weather prediction data 144 associates date and time, longitude / latitude, weather, and predicted amount of solar radiation.

天候予測データ144の日時は、現在日時よりも先の日時を示す。経度・緯度は、日射量の予測の対象となる位置を示す。天気は、予測される天気を示す。予測日射量は、該当する日時、経度・緯度において、予測される日射量である。   The date / time of the weather prediction data 144 indicates a date / time ahead of the current date / time. The longitude / latitude indicates the position for which the amount of solar radiation is to be predicted. The weather indicates the predicted weather. The predicted amount of solar radiation is the amount of solar radiation predicted at the relevant date and time, longitude and latitude.

発電機器データ145は、日射量と、この日射量によって太陽光発電機器により発電される発電量を定義する情報である。   The power generation device data 145 is information that defines the amount of solar radiation and the amount of power generated by the solar power generation device based on the amount of solar radiation.

制御部150は、収集部151と、生成部152と、算出部153とを有する。制御部150は、例えば、ASIC(Application Specific Integrated Circuit)や、FPGA(Field Programmable Gate Array)などの集積装置に対応する。また、制御部150は、例えば、CPUやMPU(Micro Processing Unit)等の電子回路に対応する。   The control unit 150 includes a collection unit 151, a generation unit 152, and a calculation unit 153. The control unit 150 corresponds to an integrated device such as an application specific integrated circuit (ASIC) or a field programmable gate array (FPGA). Moreover, the control part 150 respond | corresponds to electronic circuits, such as CPU and MPU (Micro Processing Unit), for example.

収集部151は、ネットワークを介して、外部のサーバ等にアクセスし、天候予測・実績過去データ141および天候予測データ144を収集する処理部である。なお、利用者は、入力部120を操作して、天候予測・実績過去データ141、天候予測データ144を、電力量予測装置100に入力しても良い。この場合には、収集部151は、入力部120から、天候予測・実績過去データ141、天候予測データ144を収集する。収集部151は、記憶部140に、天候予測・実績過去データ141、天候予測データ144を格納する。   The collection unit 151 is a processing unit that accesses an external server or the like via a network and collects the weather prediction / actual past data 141 and the weather prediction data 144. The user may input the weather prediction / actual past data 141 and the weather prediction data 144 to the power amount prediction apparatus 100 by operating the input unit 120. In this case, the collection unit 151 collects the weather prediction / actual past data 141 and the weather prediction data 144 from the input unit 120. The collection unit 151 stores the weather prediction / result past data 141 and the weather prediction data 144 in the storage unit 140.

生成部152は、天候予測・実績過去データ144に格納された予測日射量および実測日射量を基にして、境界値143および度数分布データ142を生成する処理部である。以下において、生成部152が、境界値143を生成する処理と、度数分布データ142を生成する処理について説明する。   The generation unit 152 is a processing unit that generates the boundary value 143 and the frequency distribution data 142 based on the predicted solar radiation amount and the actual solar radiation amount stored in the weather prediction / actual past data 144. Hereinafter, a process in which the generation unit 152 generates the boundary value 143 and a process of generating the frequency distribution data 142 will be described.

生成部152が、境界値143を生成する処理の一例について説明する。生成部152は、天候予測・実績過去データ144の未選択のレコードを選択し、選択したレコードから、予測日射量を取得し、第1晴天指数を算出する。第1晴天指数は、式(1)により算出される。式(1)において、大気外日射量は、大気圏外まで達する日射量である。例えば、生成部152は、外部のサーバ等から、大気外日射量に関する情報を予め取得しておく。   An example of processing in which the generation unit 152 generates the boundary value 143 will be described. The generation unit 152 selects an unselected record of the weather prediction / actual past data 144, acquires the predicted solar radiation amount from the selected record, and calculates the first clear sky index. The first clear sky index is calculated by equation (1). In equation (1), the solar radiation amount outside the atmosphere is the amount of solar radiation reaching outside the atmosphere. For example, the generation unit 152 acquires in advance information on the amount of solar radiation outside the atmosphere from an external server or the like.

第1晴天指数=予測日射量/大気外日射量・・・(1)   1st clear sky index = predicted amount of solar radiation / amount of solar radiation outside the atmosphere (1)

生成部152は、選択したレコードに含まれる実測日射量と、予測日射量とを比較して、実測日射量に対する予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを特定する。生成部152は、特定した予測誤差の方向と、式(1)により算出した第1晴天指数とを対応付ける。生成部152は、未選択のレコードを選択し、予測誤差の方向と第1晴天指数とを対応付ける処理を繰り返し実行し、予測誤差の方向が切り替わる第1晴天指数を、境界値143として特定する。生成部152は、境界値143を、記憶部140に格納する。   The generation unit 152 compares the measured solar radiation amount included in the selected record with the predicted solar radiation amount, and determines whether the prediction error of the predicted solar radiation amount with respect to the actual solar radiation amount is increasing or decreasing. To do. The generation unit 152 associates the identified direction of the prediction error with the first clear sky index calculated by Expression (1). The generation unit 152 selects an unselected record, repeatedly executes the process of associating the direction of the prediction error with the first clear sky index, and specifies the first clear sky index at which the direction of the prediction error is switched as the boundary value 143. The generation unit 152 stores the boundary value 143 in the storage unit 140.

なお、上記処理では、経度・緯度を固定せずに、境界値143を求める例を示したが、経度・緯度を固定して、経度・緯度毎に境界値143を求めても良い。この場合には、生成部152は、ある同一の経度・緯度に対応する、実測日射量と、予測日射量とを比較して、増加方向となるのか、減少方向となるのかを特定する。   In the above processing, the example in which the boundary value 143 is obtained without fixing the longitude / latitude is shown. However, the boundary value 143 may be obtained for each longitude / latitude by fixing the longitude / latitude. In this case, the generation unit 152 compares the measured solar radiation amount and the predicted solar radiation amount corresponding to a certain same longitude and latitude, and specifies whether the direction is an increasing direction or a decreasing direction.

続いて、生成部152が、度数分布データ142を生成する処理の一例について説明する。生成部152は、天候予測・実績過去データ144の未選択のレコードを選択し、選択したレコードから、予測日射量および実測日射量を選択し、選択した情報を基にして、晴天指数誤差を算出する。晴天指数誤差は、式(2)により算出される。式(2)に含まれる第2晴天指数は、式(3)により算出される。   Next, an example of processing in which the generation unit 152 generates the frequency distribution data 142 will be described. The generation unit 152 selects an unselected record of the weather prediction / actual past data 144, selects a predicted solar radiation amount and an actual solar radiation amount from the selected record, and calculates a clear sky index error based on the selected information. To do. The clear sky index error is calculated by equation (2). The second clear sky index included in Equation (2) is calculated by Equation (3).

晴天指数誤差=第2晴天指数−第1晴天指数・・・(2)   Clear weather index error = 2nd clear sky index-1st clear sky index (2)

第2晴天指数=実測日射量/大気外日射量・・・(3)   2nd clear sky index = actual solar radiation amount / outside air solar radiation amount (3)

生成部152は、未選択のレコードを選択し、晴天指数誤差を算出する処理を繰り返し実行し、晴天指数誤差と、この晴天指数誤差の出現頻度とを対応付けた度数分布データ142を生成する。生成部152は、晴天指数誤差がマイナスとなる場合の、第1度数分布データ142aと、晴天指数誤差がプラスとなる場合の、第2度数分布データ142bとを生成する。度数分布データ142は、第1度数分布データ142aと、第2度数分布データ142bとを含むデータとなる。   The generation unit 152 selects an unselected record, repeatedly executes a process of calculating a clear sky index error, and generates frequency distribution data 142 in which the clear sky index error is associated with the appearance frequency of the clear sky index error. The generation unit 152 generates the first frequency distribution data 142a when the clear sky index error is negative, and the second frequency distribution data 142b when the clear sky index error is positive. The frequency distribution data 142 is data including first frequency distribution data 142a and second frequency distribution data 142b.

図4は、第1度数分布データの一例を示す図である。図4に示すように、第1度数分布データ142aの縦軸は出現頻度を示し、横軸は晴天指数誤差に対応する。第1度数分布データ142aの横軸は、右側から左側に向かって、晴天指数予測誤差が減少方向に変化する。   FIG. 4 is a diagram illustrating an example of the first frequency distribution data. As shown in FIG. 4, the vertical axis of the first frequency distribution data 142a indicates the appearance frequency, and the horizontal axis corresponds to the clear sky index error. On the horizontal axis of the first frequency distribution data 142a, the clear sky index prediction error changes in a decreasing direction from the right side to the left side.

図5は、第2度数分布データの一例を示す図である。第2度数分布データ142bの縦軸は出現頻度を示し、横軸は晴天指数誤差に対応する。第2度数分布データ142bの横軸は、左側から右側に向かって、晴天指数予測誤差が増加方向に変化する。   FIG. 5 is a diagram illustrating an example of the second frequency distribution data. The vertical axis of the second frequency distribution data 142b indicates the appearance frequency, and the horizontal axis corresponds to the clear sky index error. On the horizontal axis of the second frequency distribution data 142b, the clear sky index prediction error changes in the increasing direction from the left to the right.

なお、上記処理では、経度・緯度を固定せずに、第1度数分布データ142aおよび第2度数分布データ142bを求める例を示したが、経度・緯度を固定して、経度・緯度毎に第1度数分布データ142aおよび第2度数分布データ142bを求めてもよい。   In the above process, the first frequency distribution data 142a and the second frequency distribution data 142b are obtained without fixing the longitude / latitude. However, the longitude / latitude is fixed and the first frequency distribution data 142a and the second frequency distribution data 142b are fixed. The first frequency distribution data 142a and the second frequency distribution data 142b may be obtained.

図1の説明に戻る。算出部153は、記憶部140に格納された各情報141〜145を基にして、太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する処理部である。算出部153は、誤差方向判定部154と、限界日射量算出部155と、電力コストリスク方向推定部156と、充放電計画最適化部157とを有する。   Returning to the description of FIG. The calculation unit 153 is a processing unit that calculates a power generation amount generated by the photovoltaic power generation apparatus and a prediction error of the power generation amount based on the pieces of information 141 to 145 stored in the storage unit 140. The calculation unit 153 includes an error direction determination unit 154, a limit solar radiation amount calculation unit 155, a power cost risk direction estimation unit 156, and a charge / discharge plan optimization unit 157.

誤差方向判定部154は、日射量の予測誤差の方向が、増加方向となるのか、減少方向となるのかを判定する処理部である。誤差方向判定部154は、判定結果を限界日射量算出部155および電力コストリスク方向推定部156に出力する。   The error direction determination unit 154 is a processing unit that determines whether the direction of the prediction error of the solar radiation amount is an increasing direction or a decreasing direction. The error direction determination unit 154 outputs the determination result to the limit solar radiation amount calculation unit 155 and the power cost risk direction estimation unit 156.

誤差方向判定部154は、予測の対象となる日時および経度・緯度を特定し、特定した日時および経度・緯度の予測日射量を、天候予測データ144から取得する。以下の説明では、予測の対象となる日時および経度・緯度を、適宜、対象日時と表記する。誤差方向判定部154は、天候予測テーブル144から取得した予測日射量と、式(1)とを基にして、対象日時の第1晴天指数を算出する。   The error direction determination unit 154 specifies the date / time and the longitude / latitude to be predicted, and acquires the predicted date / time and the predicted solar radiation amount of the longitude / latitude from the weather prediction data 144. In the following description, the date / time and longitude / latitude to be predicted are appropriately described as the target date / time. The error direction determination unit 154 calculates the first clear sky index of the target date and time based on the predicted solar radiation amount acquired from the weather prediction table 144 and the formula (1).

誤差方向判定部154は、第1晴天指数と、境界値143とを比較し、第1晴天指数が境界値143よりも大きい場合には、対象日時の日射量の予測誤差の方向が減少方向になると判定する。一方、誤差方向判定部154は、第1晴天指数と、境界値143とを比較し、第1晴天指数が境界値143よりも大きくない場合には、対象日時の日射量の予測誤差の方向が増加方向になると判定する。   The error direction determination unit 154 compares the first clear sky index and the boundary value 143. If the first clear sky index is larger than the boundary value 143, the direction of the prediction error of the solar radiation amount at the target date and time is in the decreasing direction. It is determined that On the other hand, the error direction determination unit 154 compares the first clear sky index and the boundary value 143. If the first clear sky index is not larger than the boundary value 143, the direction of the prediction error of the solar radiation amount at the target date and time is It is determined that the direction is increasing.

限界日射量算出部155は、限界日射量を算出する処理部である。限界日射量算出部155は、限界日射量を算出する前に、誤差方向判定部154により判定された日射量の予測誤差の方向に基づき、予測誤差方向に累積確率が保障確率と等しくなる晴天指数誤差yを算出する。なお、保障確率は、利用者によって予め指示されているものとする。 The limit solar radiation amount calculation unit 155 is a processing unit that calculates the limit solar radiation amount. The marginal solar radiation amount calculation unit 155 calculates a clear sky index based on the direction of the prediction error of the solar radiation amount determined by the error direction determination unit 154 before the critical solar radiation amount is calculated, so that the cumulative probability is equal to the guaranteed probability in the prediction error direction. to calculate the error y 0. It is assumed that the guarantee probability is instructed in advance by the user.

図6および図7は、限界日射量算出部の処理を説明するための図である。日射量の予測誤差の方向が、減少方向であると判定された場合の、限界日射量算出部155の処理について説明する。図6に示すように、限界日射量算出部155は、第1度数分布データ142a上に、晴天指数誤差yを設定する。限界日射量算出部155は、式(4)に基づいて算出される累積確率が、保障確率と等しくなるように、晴天指数誤差yを調整する。式(4)において、第1領域は、縦軸10a、横軸10b、線10cによって囲まれる領域である。第2領域は、縦軸10a、横軸10b、線10c、線10dによって囲まれる領域である。 6 and 7 are diagrams for explaining the processing of the limit solar radiation amount calculation unit. The process of the limit solar radiation amount calculation unit 155 when it is determined that the direction of the prediction error of the solar radiation amount is the decreasing direction will be described. As shown in FIG. 6, the limit solar radiation amount calculation unit 155, on the first frequency distribution data 142a, it sets the sunny index error y 0. The limit solar radiation amount calculation unit 155 adjusts the clear sky index error y 0 so that the cumulative probability calculated based on the equation (4) is equal to the guarantee probability. In Expression (4), the first region is a region surrounded by the vertical axis 10a, the horizontal axis 10b, and the line 10c. The second region is a region surrounded by the vertical axis 10a, the horizontal axis 10b, the line 10c, and the line 10d.

累積確率=第2領域の面積/第1領域の面積×100(%)・・・(4)   Cumulative probability = area of second region / area of first region × 100 (%) (4)

日射量の予測誤差の方向が、増加方向であると判定された場合の、限界日射量算出部155の処理について説明する。図7に示すように、限界日射量算出部155は、第2度数分布データ142b上に、晴天指数誤差yを設定する。限界日射量算出部155は、式(5)に基づいて算出される累積確率が、保障確率と等しくなるように、晴天指数誤差yを調整する。式(5)において、第3領域は、縦軸20a、横軸20b、線20cによって囲まれる領域である。第4領域は、縦軸20a、横軸20b、線20c、線20dによって囲まれる領域である。 The process of the limit solar radiation amount calculation unit 155 when it is determined that the direction of the prediction error of the solar radiation amount is the increasing direction will be described. As shown in FIG. 7, the limit solar radiation amount calculation unit 155, on the second frequency distribution data 142b, it sets the sunny index error y 0. The limit solar radiation amount calculation unit 155 adjusts the clear sky index error y 0 so that the cumulative probability calculated based on the equation (5) becomes equal to the guarantee probability. In Expression (5), the third region is a region surrounded by the vertical axis 20a, the horizontal axis 20b, and the line 20c. The fourth region is a region surrounded by the vertical axis 20a, the horizontal axis 20b, the line 20c, and the line 20d.

累積確率=第4領域の面積/第3領域の面積×100(%)・・・(5)   Cumulative probability = Area of the fourth region / Area of the third region × 100 (%) (5)

限界日射量算出部155は、累積確率が保障確率と等しくなる晴天指数誤差yを求めた後に、式(6)に基づいて、第3晴天指数を算出する。第3晴天指数は、累積確率が保障確率と等しくなる晴天指数に対応するものである。 The limit solar radiation amount calculation unit 155 calculates the third clear sky index based on the equation (6) after obtaining the clear sky index error y 0 in which the cumulative probability is equal to the guaranteed probability. The third clear sky index corresponds to a clear sky index whose cumulative probability is equal to the guaranteed probability.

第3晴天指数=対象日の第1晴天指数+晴天指数誤差y・・・(6) 3rd clear sky index = 1st clear sky index + clear sky index error y 0 (6)

限界日射量算出部155は、式(7)に基づいて、限界日射量を算出する。限界日射量は、保障確率を満たす予想日射量である。   The limit solar radiation amount calculation unit 155 calculates the limit solar radiation amount based on the equation (7). The marginal solar radiation amount is an expected solar radiation amount that satisfies the guarantee probability.

限界日射量=第3晴天指数×大気外日射量・・・(7)   Marginal solar radiation amount = 3rd clear sky index x solar radiation amount outside the atmosphere ... (7)

限界日射量算出部155は、経度・緯度を固定して、日時を変更しつつ、上記処理を繰り返し実行することで、ある経度・緯度に関する日時毎の限界日射量を算出する。限界日射量算出部155は、ある経度・緯度に関する日時毎の限界日射量の情報を、充放電計画最適化部157に出力する。例えば、ある経度・緯度は、太陽光発電機器が設置された位置に対応する。   The limit solar radiation amount calculation unit 155 calculates the limit solar radiation amount for each date and time related to a certain longitude and latitude by repeatedly executing the above processing while fixing the longitude and latitude and changing the date and time. The limit solar radiation amount calculation unit 155 outputs information on the limit solar radiation amount for each date and time related to a certain longitude / latitude to the charge / discharge plan optimization unit 157. For example, a certain longitude / latitude corresponds to a position where a photovoltaic power generation device is installed.

電力コストリスク方向推定部156は、誤差方向判定部154の判定結果を基にして、電力コストが増える方向か減る方向かを推定する処理部である。電力コストリスク方向推定部156は、推定結果を、充放電計画最適化部157に出力する。   The power cost risk direction estimation unit 156 is a processing unit that estimates whether the power cost increases or decreases based on the determination result of the error direction determination unit 154. The power cost risk direction estimation unit 156 outputs the estimation result to the charge / discharge plan optimization unit 157.

電力コストリスク方向推定部156は、日射量の予測誤差の方向が増加方向となる場合には、電力コストが減る方向であると推定する。日射量の予測誤差の方向が増加方向となる場合には、太陽光発電機器による発電量が増えることを意味し、発電量が増えると、電力会社から買う電力量が減少するため、電力コストが減る。   The power cost risk direction estimation unit 156 estimates that the power cost is in a decreasing direction when the direction of the solar radiation amount prediction error is an increasing direction. When the amount of solar radiation forecast error increases, it means that the amount of power generated by solar power generation equipment increases.If the amount of power generation increases, the amount of power purchased from the power company decreases, so the power cost decrease.

一方、電力コストリスク方向推定部156は、日射量の予測誤差の方向が減少方向となる場合には、電力コストが増える方向であると推定する。日射量の予測誤差の方向が減少方向となる場合には、太陽光発電機器による発電量が減ることを意味し、発電量が減ると、電力会社から買う電力量が増加するため、電力コストが増える。   On the other hand, the power cost risk direction estimation unit 156 estimates that the power cost increases when the direction of the prediction error of the solar radiation amount decreases. When the amount of solar radiation forecast error is decreasing, this means that the amount of power generated by solar power generation equipment will decrease.If the amount of power generation decreases, the amount of power purchased from the power company will increase. Increase.

充放電計画最適化部157は、限界日射量と、発電機器データ145とを基にして、対象日時における、太陽光発電機器が発電する発電量を算出する。なお、上記のように、限界日射量は、日射量の予測誤差を加味した日射量であるため、充放電計画最適化部157が算出する発電量は、発電量の予測誤差を加味したものとなる。充放電計画最適化部157は、各対象日時において、太陽光発電機器が発電する発電量をそれぞれ算出する。   The charge / discharge plan optimization unit 157 calculates the power generation amount generated by the solar power generation device at the target date and time based on the limit solar radiation amount and the power generation device data 145. Note that, as described above, the limit solar radiation amount is the solar radiation amount that takes into account the prediction error of the solar radiation amount. Therefore, the power generation amount calculated by the charge / discharge plan optimization unit 157 includes the prediction error of the power generation amount. Become. The charge / discharge plan optimization unit 157 calculates the amount of power generated by the solar power generation device at each target date and time.

充放電計画最適化部157は、各対象日時において、電力需要量と発電量とを用いてピーク電力削減が最大となる蓄電池放電量を、線形計画法を用いて算出する。例えば、ピーク電力削減が最大となる各時刻の蓄電池充放電量が、充放電計画となる。   The charge / discharge plan optimization unit 157 calculates, using linear programming, the storage battery discharge amount that maximizes the peak power reduction using the power demand amount and the power generation amount at each target date and time. For example, the charge / discharge amount of the storage battery at each time when the peak power reduction is the maximum is the charge / discharge plan.

図8は、充放電計画最適化部の処理を説明するための図である。図8の横軸は時刻に対応し、縦軸は電力に対応する。図8に示す例では、d1からd2を差し引いた値が、削減されるピーク電力に対応する。   FIG. 8 is a diagram for explaining the processing of the charge / discharge plan optimization unit. The horizontal axis in FIG. 8 corresponds to time, and the vertical axis corresponds to power. In the example shown in FIG. 8, the value obtained by subtracting d2 from d1 corresponds to the peak power to be reduced.

図8において、Dは、各時刻(対象日時)における需要電力に対応する。Mは、各時刻における太陽光発電機器における発電量に対応する。Xは、各時刻における電池充放電量に対応する。蓄電池容量をA、蓄電池初期充電量をBとする。このとき、ピーク電力削減が最大となる充放電計画を、式(8)の線形計画問題の解を解くことで求める。D、A、Bの情報は、予め指定されているものとする。   In FIG. 8, D corresponds to the power demand at each time (target date and time). M corresponds to the power generation amount in the solar power generation device at each time. X corresponds to the battery charge / discharge amount at each time. The storage battery capacity is A, and the storage battery initial charge is B. At this time, the charge / discharge plan that maximizes the peak power reduction is obtained by solving the linear programming problem of Equation (8). Information on D, A, and B is specified in advance.

Figure 2018018329
Figure 2018018329

式(8)に含まれるZは、需要電力を発電量と蓄電池放電とで削減した電力値の上界であり、線形計画問題の解minZが、最小のピーク電力となる。問題後半のA、Bに関する制約条件は、蓄電池の充電量が容量以下にならないことを意味している。   Z included in the equation (8) is an upper limit of the power value obtained by reducing the demand power by the power generation amount and the battery discharge, and the solution minZ of the linear programming problem is the minimum peak power. The constraint conditions related to A and B in the latter half of the problem mean that the charge amount of the storage battery does not fall below the capacity.

充放電計画最適化部157は、ピーク電力削減が最大となる充放電計画を算出する。そして、充放電計画最適化部157は、充放電計画の算出結果と、電力コストリスク方向推定部156の推定結果とを対応付けて、出力する。充放電計画最適化部157は、表示部130に出力しても良いし、ネットワークを介して、他の端末装置に出力しても良い。   The charge / discharge plan optimization unit 157 calculates a charge / discharge plan that maximizes the peak power reduction. Then, the charge / discharge plan optimization unit 157 associates and outputs the calculation result of the charge / discharge plan and the estimation result of the power cost risk direction estimation unit 156. The charge / discharge plan optimization unit 157 may output to the display unit 130 or may output to another terminal device via a network.

例えば、電力コストリスク方向推定部156により、電力コストが減る方向であると推定されている場合には、算出した充放電計画に対応する電力コストは、電力コストの上限であることを意味する。   For example, when the power cost risk direction estimation unit 156 estimates that the power cost is decreasing, the power cost corresponding to the calculated charge / discharge plan is the upper limit of the power cost.

一方、電力コストリスク方向推定部156により、電力コストが増える方向であると推定されている場合には、算出した充放電計画に対応する電力コストは、電力コストの下限であることを意味する。   On the other hand, when the power cost risk direction estimation unit 156 estimates that the power cost is increasing, the power cost corresponding to the calculated charge / discharge plan is the lower limit of the power cost.

次に、本実施例に係る電力量予測装置の処理手順の一例について説明する。図9は、電力量予測装置が境界値および度数分布データを生成する処理手順を示すフローチャートである。図9に示すように、電力量予測装置100の生成部152は、天候予測・実績過去データ141を参照し、過去の予測日射量と大気外日射量とを基にして、複数の第1晴天指数を算出する(ステップS101)。   Next, an example of a processing procedure of the power amount prediction apparatus according to the present embodiment will be described. FIG. 9 is a flowchart illustrating a processing procedure in which the power amount prediction apparatus generates boundary values and frequency distribution data. As illustrated in FIG. 9, the generation unit 152 of the power amount prediction apparatus 100 refers to the weather prediction / result past data 141, and based on the past predicted solar radiation amount and the outdoor solar radiation amount, An index is calculated (step S101).

生成部152は、過去の実測日射量と大気外日射量とを基にして、複数の第2晴天指数を算出する(ステップS102)。生成部152は、過去の同一日時における第1晴天指数と第2晴天指数との関係から、境界値143を算出する(ステップS103)。   The generation unit 152 calculates a plurality of second clear sky indexes based on the past actually measured solar radiation amount and the atmospheric solar radiation amount (step S102). The generation unit 152 calculates the boundary value 143 from the relationship between the first clear sky index and the second clear sky index at the same date and time in the past (step S103).

生成部152は、過去の同一日時における第1晴天指数および第2晴天指数の差分に基づき、晴天指数誤差を算出する(ステップS104)。生成部152は、晴天指数誤差の度数分布データ142を生成する(ステップS105)。   The generation unit 152 calculates a clear sky index error based on the difference between the first clear sky index and the second clear sky index at the same date and time in the past (step S104). The generating unit 152 generates the frequency distribution data 142 of the clear sky index error (step S105).

図10は、電力量予測装置が充放電計画を算出する処理手順を示すフローチャートである。図10に示すように、電力量予測装置100の算出部153は、対象日時の予測日射量と大気外日射量とを基にして、対象日時における第1晴天指数を算出する(ステップS201)。算出部153は、第1晴天指数が境界値143よりも大きいか否かを判定する(ステップS202)。   FIG. 10 is a flowchart illustrating a processing procedure in which the power amount prediction apparatus calculates a charge / discharge plan. As illustrated in FIG. 10, the calculation unit 153 of the power amount prediction apparatus 100 calculates the first clear sky index at the target date and time based on the predicted solar radiation amount and the atmospheric solar radiation amount at the target date and time (step S201). The calculation unit 153 determines whether or not the first clear sky index is larger than the boundary value 143 (step S202).

算出部153は、第1晴天指数が境界値143よりも大きくない場合には(ステップS202,No)、誤差方向を「日射量の増える方向」と推定し(ステップS203)、ステップS205に移行する。一方、第1晴天指数が境界値143よりも大きい場合には(ステップS202,Yes)、誤差方向を「日射量の減る方向」と推定し(ステップS204)、ステップS205に移行する。   When the first clear sky index is not larger than the boundary value 143 (No at Step S202), the calculation unit 153 estimates the error direction as “a direction in which the amount of solar radiation increases” (Step S203), and proceeds to Step S205. . On the other hand, when the first clear sky index is larger than the boundary value 143 (step S202, Yes), the error direction is estimated as “a direction in which the amount of solar radiation decreases” (step S204), and the process proceeds to step S205.

算出部153は、推定した誤差方向に応じて度数分布データ142を選択し、累積確率が保障確率と等しくなる晴天指数誤差yを特定する(ステップS205)。例えば、ステップS205において、誤差方向が減少方向である場合には、算出部153は、第1度数分布データ142aを選択して、晴天指数誤差yを特定する。誤差方向が増加方向である場合には、算出部153は、第2度数分布データ142bを選択して、晴天指数誤差yを特定する。 The calculation unit 153 selects the frequency distribution data 142 in accordance with the estimated error direction, and specifies the clear sky index error y 0 in which the cumulative probability is equal to the guaranteed probability (step S205). For example, in step S205, when the error direction is the decreasing direction, the calculating unit 153 selects the first frequency distribution data 142a, identifies the sunny index error y 0. If the error direction is the increasing direction, the calculating unit 153 selects the second frequency distribution data 142b, identifies a sunny index error y 0.

算出部153は、対象日時の第1晴天指数に晴天指数誤差yを加算することで、第3晴天指数を算出する(ステップS206)。算出部153は、第3晴天指数と、大気外日射量とを乗算することで、限界日射量を算出する(ステップS207)。算出部153は、推定した誤差方向を基にして、電力コストが増える方向か減る方向かを推定する(ステップS208)。算出部153は、充放電計画を生成する(ステップS209)。 Calculating unit 153, by adding the clear sky index error y 0 in the first sunny index of the target date, and calculates a third sunny index (step S206). The calculating unit 153 calculates the limit solar radiation amount by multiplying the third clear sky index and the atmospheric solar radiation amount (step S207). The calculation unit 153 estimates whether the power cost increases or decreases based on the estimated error direction (step S208). The calculation unit 153 generates a charge / discharge plan (step S209).

次に、本実施例に係る電力量予測装置100の効果について説明する。電力量予測装置100は、天候予測・実績過去データ141の予測日射量と実測日射量とを基にして、境界値143を生成する。そして、電力量予測装置100は、対象日時の晴天指数と境界値143との比較により、日射量の予測誤差が増加方向に振れるのか減少方向に振れるのかを推定する。電力量予測装置100は、推定した方向により保障確率を満たす日射量の予測誤差を特定し、発電量を求めることで、正規分布による誤差分布の仮定を用いないで、発電量の予測誤差を求めることができる。   Next, the effect of the power amount prediction apparatus 100 according to the present embodiment will be described. The power amount prediction apparatus 100 generates the boundary value 143 based on the predicted solar radiation amount and the actual solar radiation amount of the weather prediction / actual past data 141. And the electric energy prediction apparatus 100 estimates whether the prediction error of the amount of solar radiation fluctuates in an increasing direction or a decreasing direction by comparing the clear sky index of the target date and time with the boundary value 143. The power amount prediction apparatus 100 determines a prediction error of the amount of solar radiation that satisfies the guarantee probability based on the estimated direction, and calculates a power generation amount, thereby obtaining a power generation amount prediction error without using an error distribution assumption based on a normal distribution. be able to.

電力量予測装置100は、天候予測・実績過去データ141を基にして、予測誤差の方向と第1晴天指数とを対応付ける処理を繰り返し実行し、予測誤差の方向が切り替わる第1晴天指数を、境界値143として特定する。このため、最適な境界値によって、予測誤差の方向を精度良く推定することができる。   The power amount prediction apparatus 100 repeatedly executes a process of associating the direction of the prediction error with the first clear sky index based on the weather prediction / actual past data 141, and the first clear sky index in which the direction of the prediction error is switched to the boundary. It is specified as a value 143. For this reason, it is possible to accurately estimate the direction of the prediction error using the optimum boundary value.

なお、従来技術では、保障確率を満たす充放電計画を生成する場合に、通常の予測誤差分布を用いると、予測誤差が増加する方向と減少する方向の双方について、充放電計画を生成するため、計算コストが高くなっていた。これに対して、電力量予測装置100では、予測誤差の方向を一方に決定し、決定した方向について、充放電計画を生成すればよいので、計算コストを従来技術のものと比較して、削減することができる。   In the prior art, when generating a charge / discharge plan that satisfies the guarantee probability, if a normal prediction error distribution is used, the charge / discharge plan is generated for both the direction in which the prediction error increases and the direction in which the prediction error decreases. The calculation cost was high. On the other hand, in the electric energy prediction apparatus 100, the direction of the prediction error is determined to one side, and a charge / discharge plan may be generated for the determined direction. can do.

ところで、上記実施例では、対象日時の晴天指数と境界値143との比較により、日射量の予測誤差が増加方向に振れるのか減少方向に振れるのかを推定していたが、対象日時の天気に基づき、推定を行っても良い。例えば、電力量予測装置100の誤差方向判定部154は、対象日時の天気が「晴れ」の場合には、対象日時の日射量の予測誤差の方向が「減少方向」になると判定する。一方、誤差方向判定部154は、対象日時の天気が「曇り・雨」の場合には、対象日時の日射量の予測誤差の方向が「増加方向」になると判定する。このように、境界値と比較する代わりに、天気に基づいて予測誤差の方向を算出することで、処理を簡略化して、処理負荷を軽減することができる。   By the way, in the said Example, it was estimated by the comparison with the clear sky index | exponent of the object date and time, and the boundary value 143, whether the prediction error of the solar radiation amount fluctuated in the increase direction or the decrease direction, but based on the weather of the object date The estimation may be performed. For example, the error direction determination unit 154 of the power amount prediction apparatus 100 determines that the direction of the prediction error of the amount of solar radiation at the target date / time is “decreasing direction” when the weather at the target date / time is “sunny”. On the other hand, when the weather at the target date / time is “cloudy / rainy”, the error direction determination unit 154 determines that the direction of the prediction error of the amount of solar radiation at the target date / time is “increase direction”. In this way, by calculating the direction of the prediction error based on the weather instead of comparing with the boundary value, the processing can be simplified and the processing load can be reduced.

次に、上記実施例に示した電力量予測装置100と同様の機能を実現する電力量予測プログラムを実行するハードウェア構成の一例について説明する。図11は、電力量予測プログラムを実行するハードウェア構成の一例を示す図である。   Next, an example of a hardware configuration that executes a power amount prediction program that realizes the same function as that of the power amount prediction apparatus 100 described in the above embodiment will be described. FIG. 11 is a diagram illustrating an example of a hardware configuration that executes the power amount prediction program.

図11に示すように、コンピュータ200は、各種演算処理を実行するCPU201と、ユーザからのデータの入力を受け付ける入力装置202と、ディスプレイ203とを有する。また、コンピュータ200は、記憶媒体からプログラム等を読み取る読み取り装置204と、ネットワークを介して他のコンピュータとの間でデータの授受を行うインタフェース装置205とを有する。また、コンピュータ200は、各種情報を一時記憶するRAM206と、ハードディスク装置207とを有する。そして、各装置201〜207は、バス208に接続される。   As illustrated in FIG. 11, the computer 200 includes a CPU 201 that executes various arithmetic processes, an input device 202 that receives input of data from a user, and a display 203. The computer 200 also includes a reading device 204 that reads programs and the like from a storage medium, and an interface device 205 that exchanges data with other computers via a network. The computer 200 also includes a RAM 206 that temporarily stores various information and a hard disk device 207. The devices 201 to 207 are connected to the bus 208.

ハードディスク装置207は、生成プログラム207aおよび算出プログラム207bを有する。CPU201は、生成プログラム207aおよび算出プログラム207bを読み出してRAM206に展開する。   The hard disk device 207 has a generation program 207a and a calculation program 207b. The CPU 201 reads the generation program 207 a and the calculation program 207 b and expands them in the RAM 206.

生成プログラム207aは、生成プロセス206aとして機能する。算出プログラム207bは、算出プロセス206bとして機能する。   The generation program 207a functions as a generation process 206a. The calculation program 207b functions as a calculation process 206b.

例えば、生成プロセス206aの処理は、生成部152の処理に対応する。算出プロセス206bの処理は算出部153の処理に対応する。   For example, the process of the generation process 206 a corresponds to the process of the generation unit 152. The process of the calculation process 206b corresponds to the process of the calculation unit 153.

なお、生成プログラム207a、算出プログラム207bについては、必ずしも最初からハードディスク装置207に記憶させておかなくても良い。例えば、コンピュータ200に挿入されるフレキシブルディスク(FD)、CD−ROM、DVDディスク、光磁気ディスク、ICカードなどの「可搬用の物理媒体」に各プログラムを記憶させておく。そして、コンピュータ200が各プログラム207a,207bを読み出して実行するようにしてもよい。   Note that the generation program 207a and the calculation program 207b are not necessarily stored in the hard disk device 207 from the beginning. For example, each program is stored in a “portable physical medium” such as a flexible disk (FD), a CD-ROM, a DVD disk, a magneto-optical disk, and an IC card inserted into the computer 200. Then, the computer 200 may read and execute the programs 207a and 207b.

以上の各実施例を含む実施形態に関し、さらに以下の付記を開示する。   The following supplementary notes are further disclosed with respect to the embodiments including the above examples.

(付記1)コンピュータに、
予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、前記実測日射量に対する前記予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成し、
前記予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、前記太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する
処理を実行させることを特徴とする電力量予測プログラム。
(Supplementary note 1)
Based on the predicted solar radiation amount predicted at the prediction target point and date and the actual solar radiation amount actually measured at the prediction target point and date and time, the prediction error of the predicted solar radiation amount with respect to the actual solar radiation amount increases. Generate prediction error information that identifies whether the direction is the direction or the direction of decrease for each weather forecast,
Based on the prediction error information, the solar radiation amount prediction information, and the solar power generation device information, a process of calculating the power generation amount generated by the solar power generation device and the power generation amount prediction error is executed. A power amount prediction program characterized by that.

(付記2)前記天候予測は、晴れ、曇り、雨の区分を含む天気予測であることを特徴とする付記1に記載の電力量予測プログラム。 (Supplementary note 2) The electric energy prediction program according to supplementary note 1, wherein the weather prediction is a weather prediction including clear, cloudy, and rainy categories.

(付記3)前記天候予測は、日射量に関連した指数であることを特徴とする付記1に記載の電力量予測プログラム。 (Supplementary note 3) The electric energy prediction program according to supplementary note 1, wherein the weather prediction is an index related to the amount of solar radiation.

(付記4)前記天候予測は、前記予測日射量を大気外日射量で除算した晴天指数であることを特徴とする付記1に記載の電力量予測プログラム。 (Supplementary note 4) The electric energy prediction program according to supplementary note 1, wherein the weather prediction is a clear sky index obtained by dividing the predicted solar radiation amount by an outdoor solar radiation amount.

(付記5)前記算出する処理は、前記日射量の予測情報を基にして、前記予測誤差情報に含まれる増加方向の予測誤差、または、前記予測誤差情報に含まれる減少方向の予測誤差を利用して、前記発電量の予測誤差を算出することを特徴とする付記1に記載の電力量予測プログラム。 (Supplementary Note 5) The calculation process uses an increase direction prediction error included in the prediction error information or a decrease direction prediction error included in the prediction error information based on the solar radiation amount prediction information. The power amount prediction program according to appendix 1, wherein a prediction error of the power generation amount is calculated.

(付記6)予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、前記実測日射量に対する前記予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成する生成部と、
前記予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、前記太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する算出部と
を有することを特徴とする電力量予測装置。
(Supplementary note 6) Prediction of the predicted solar radiation amount with respect to the actual solar radiation amount based on the predicted solar radiation amount predicted at the prediction target point and date and the actual solar radiation amount actually measured at the prediction target point and date and time A generation unit that generates prediction error information that specifies whether the error is in an increasing direction or a decreasing direction for each weather prediction;
Based on the prediction error information, solar radiation amount prediction information, and solar power generation device information, a power generation amount generated by the solar power generation device and a calculation unit that calculates a prediction error of the power generation amount An apparatus for predicting electric energy, comprising:

(付記7)前記天候予測は、晴れ、曇り、雨の区分を含む天気予測であることを特徴とする付記6に記載の電力量予測装置。 (Supplementary note 7) The electric energy prediction apparatus according to supplementary note 6, wherein the weather prediction is a weather prediction including clear, cloudy, and rainy categories.

(付記8)前記天候予測は、日射量に関連した指数であることを特徴とする付記6に記載の電力量予測装置。 (Supplementary note 8) The electric energy prediction device according to supplementary note 6, wherein the weather prediction is an index related to the amount of solar radiation.

(付記9)前記天候予測は、前記予測日射量を大気外日射量で除算した晴天指数であることを特徴とする付記6に記載の電力量予測装置。 (Supplementary note 9) The electric energy prediction apparatus according to supplementary note 6, wherein the weather prediction is a clear sky index obtained by dividing the predicted solar radiation amount by the atmospheric solar radiation amount.

(付記10)前記算出部は、前記日射量の予測情報を基にして、前記予測誤差情報に含まれる増加方向の予測誤差、または、前記予測誤差情報に含まれる減少方向の予測誤差を利用して、前記発電量の予測誤差を算出することを特徴とする付記6に記載の電力量予測装置。 (Supplementary Note 10) The calculation unit uses an increase direction prediction error included in the prediction error information or a decrease direction prediction error included in the prediction error information based on the solar radiation amount prediction information. The power amount prediction apparatus according to appendix 6, wherein a prediction error of the power generation amount is calculated.

(付記11)コンピュータが実行する電力量予測方法であって、
予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、前記実測日射量に対する前記予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成し、
前記予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、前記太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する
処理を実行することを特徴とする電力量予測方法。
(Additional remark 11) It is the electric energy prediction method which a computer performs, Comprising:
Based on the predicted solar radiation amount predicted at the prediction target point and date and the actual solar radiation amount actually measured at the prediction target point and date and time, the prediction error of the predicted solar radiation amount with respect to the actual solar radiation amount increases. Generate prediction error information that identifies whether the direction is the direction or the direction of decrease for each weather forecast,
Based on the prediction error information, the solar radiation amount prediction information, and the solar power generation device information, a process of calculating the power generation amount generated by the solar power generation device and the power generation amount prediction error is executed. An electric energy prediction method characterized by the above.

(付記12)前記天候予測は、晴れ、曇り、雨の区分を含む天気予測であることを特徴とする付記11に記載の電力量予測方法。 (Supplementary note 12) The electric energy prediction method according to supplementary note 11, wherein the weather prediction is a weather prediction including clear, cloudy, and rainy categories.

(付記13)前記天候予測は、日射量に関連した指数であることを特徴とする付記11に記載の電力量予測方法。 (Supplementary note 13) The electric energy prediction method according to supplementary note 11, wherein the weather prediction is an index related to the amount of solar radiation.

(付記14)前記天候予測は、前記予測日射量を大気外日射量で除算した晴天指数であることを特徴とする付記11に記載の電力量予測方法。 (Supplementary note 14) The electric energy prediction method according to supplementary note 11, wherein the weather prediction is a clear sky index obtained by dividing the predicted solar radiation amount by an outdoor solar radiation amount.

(付記15)前記算出する処理は、前記日射量の予測情報を基にして、前記予測誤差情報に含まれる増加方向の予測誤差、または、前記予測誤差情報に含まれる減少方向の予測誤差を利用して、前記発電量の予測誤差を算出することを特徴とする付記11に記載の電力量予測方法。 (Supplementary Note 15) The calculation process uses an increase direction prediction error included in the prediction error information or a decrease direction prediction error included in the prediction error information based on the solar radiation amount prediction information. The power amount prediction method according to appendix 11, wherein a prediction error of the power generation amount is calculated.

100 電力量予測装置
110 通信部
120 入力部
130 表示部
140 記憶部
150 制御部
151 収集部
152 生成部
153 算出部
DESCRIPTION OF SYMBOLS 100 Electric power amount prediction apparatus 110 Communication part 120 Input part 130 Display part 140 Storage part 150 Control part 151 Collection part 152 Generation part 153 Calculation part

Claims (7)

コンピュータに、
予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、前記実測日射量に対する前記予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成し、
前記予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、前記太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する
処理を実行させることを特徴とする電力量予測プログラム。
On the computer,
Based on the predicted solar radiation amount predicted at the prediction target point and date and the actual solar radiation amount actually measured at the prediction target point and date and time, the prediction error of the predicted solar radiation amount with respect to the actual solar radiation amount increases. Generate prediction error information that identifies whether the direction is the direction or the direction of decrease for each weather forecast,
Based on the prediction error information, the solar radiation amount prediction information, and the solar power generation device information, a process of calculating the power generation amount generated by the solar power generation device and the power generation amount prediction error is executed. A power amount prediction program characterized by that.
前記天候予測は、晴れ、曇り、雨の区分を含む天気予測であることを特徴とする請求項1に記載の電力量予測プログラム。   The electric energy prediction program according to claim 1, wherein the weather prediction is a weather prediction including sunny, cloudy, and rainy categories. 前記天候予測は、日射量に関連した指数であることを特徴とする請求項1に記載の電力量予測プログラム。   The electric energy prediction program according to claim 1, wherein the weather prediction is an index related to solar radiation. 前記天候予測は、前記予測日射量を大気外日射量で除算した晴天指数であることを特徴とする請求項1に記載の電力量予測プログラム。   The electric energy prediction program according to claim 1, wherein the weather prediction is a clear sky index obtained by dividing the predicted solar radiation amount by an outdoor solar radiation amount. 前記算出する処理は、前記日射量の予測情報を基にして、前記予測誤差情報に含まれる増加方向の予測誤差、または、前記予測誤差情報に含まれる減少方向の予測誤差を利用して、前記発電量の予測誤差を算出することを特徴とする請求項1に記載の電力量予測プログラム。   The process of calculating, based on the prediction information of the solar radiation amount, using the prediction error in the increasing direction included in the prediction error information, or the prediction error in the decreasing direction included in the prediction error information, The power amount prediction program according to claim 1, wherein a power generation amount prediction error is calculated. 予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、前記実測日射量に対する前記予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成する生成部と、
前記予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、前記太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する算出部と
を有することを特徴とする電力量予測装置。
Based on the predicted solar radiation amount predicted at the prediction target point and date and the actual solar radiation amount actually measured at the prediction target point and date and time, the prediction error of the predicted solar radiation amount with respect to the actual solar radiation amount increases. A generation unit that generates prediction error information specifying for each weather prediction whether it is a direction or a decreasing direction;
Based on the prediction error information, solar radiation amount prediction information, and solar power generation device information, a power generation amount generated by the solar power generation device and a calculation unit that calculates a prediction error of the power generation amount An apparatus for predicting electric energy, comprising:
コンピュータが実行する電力量予測方法であって、
予測対象の地点および日時において予測された予測日射量と、前記予測対象の地点および日時において実測された実測日射量とを基にして、前記実測日射量に対する前記予測日射量の予測誤差が、増加方向となるのか、減少方向となるのかを天候予測毎に特定した予測誤差情報を生成し、
前記予測誤差情報と、日射量の予測情報と、太陽光発電機器の情報とを基にして、前記太陽光発電機器が発電する発電量と、発電量の予測誤差とを算出する
処理を実行することを特徴とする電力量予測方法。
A method for predicting the amount of power executed by a computer,
Based on the predicted solar radiation amount predicted at the prediction target point and date and the actual solar radiation amount actually measured at the prediction target point and date and time, the prediction error of the predicted solar radiation amount with respect to the actual solar radiation amount increases. Generate prediction error information that identifies whether the direction is the direction or the direction of decrease for each weather forecast,
Based on the prediction error information, the solar radiation amount prediction information, and the solar power generation device information, a process of calculating the power generation amount generated by the solar power generation device and the power generation amount prediction error is executed. An electric energy prediction method characterized by the above.
JP2016148561A 2016-07-28 2016-07-28 Electric energy prediction program, electric energy prediction device and electric energy prediction method Active JP6772625B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016148561A JP6772625B2 (en) 2016-07-28 2016-07-28 Electric energy prediction program, electric energy prediction device and electric energy prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016148561A JP6772625B2 (en) 2016-07-28 2016-07-28 Electric energy prediction program, electric energy prediction device and electric energy prediction method

Publications (2)

Publication Number Publication Date
JP2018018329A true JP2018018329A (en) 2018-02-01
JP6772625B2 JP6772625B2 (en) 2020-10-21

Family

ID=61076876

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016148561A Active JP6772625B2 (en) 2016-07-28 2016-07-28 Electric energy prediction program, electric energy prediction device and electric energy prediction method

Country Status (1)

Country Link
JP (1) JP6772625B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108808730A (en) * 2018-06-12 2018-11-13 国网山东省电力公司聊城供电公司 Consider the distribution network system reserve capacity for load variation in power computational methods and system of photovoltaic time space distribution
CN109345046A (en) * 2018-12-05 2019-02-15 国能日新科技股份有限公司 The prediction technique and device of photovoltaic power
DE102019101184A1 (en) 2018-02-05 2019-08-08 Sumitomo Heavy Industries, Ltd. Fault diagnosis system
KR20220157852A (en) * 2021-05-21 2022-11-29 성균관대학교산학협력단 Method and apparatus for determining generated power for economic load dispatch

Citations (4)

* 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
JP2013253805A (en) * 2012-06-05 2013-12-19 Fuji Electric Co Ltd Information processing device, and control method and program for the same
JP2016038348A (en) * 2014-08-11 2016-03-22 富士通株式会社 Plan formulation method, plan formulation system and plan formulation program
JP2016111822A (en) * 2014-12-05 2016-06-20 中国電力株式会社 Prediction device

Patent Citations (4)

* 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
JP2013253805A (en) * 2012-06-05 2013-12-19 Fuji Electric Co Ltd Information processing device, and control method and program for the same
JP2016038348A (en) * 2014-08-11 2016-03-22 富士通株式会社 Plan formulation method, plan formulation system and plan formulation program
JP2016111822A (en) * 2014-12-05 2016-06-20 中国電力株式会社 Prediction device

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102019101184A1 (en) 2018-02-05 2019-08-08 Sumitomo Heavy Industries, Ltd. Fault diagnosis system
CN108808730A (en) * 2018-06-12 2018-11-13 国网山东省电力公司聊城供电公司 Consider the distribution network system reserve capacity for load variation in power computational methods and system of photovoltaic time space distribution
CN109345046A (en) * 2018-12-05 2019-02-15 国能日新科技股份有限公司 The prediction technique and device of photovoltaic power
CN109345046B (en) * 2018-12-05 2020-07-07 国能日新科技股份有限公司 Photovoltaic power prediction method and device
KR20220157852A (en) * 2021-05-21 2022-11-29 성균관대학교산학협력단 Method and apparatus for determining generated power for economic load dispatch
KR102591679B1 (en) * 2021-05-21 2023-10-19 성균관대학교산학협력단 Method and apparatus for determining generated power for economic load dispatch

Also Published As

Publication number Publication date
JP6772625B2 (en) 2020-10-21

Similar Documents

Publication Publication Date Title
JP5940922B2 (en) Renewable energy prediction device
JP2018018329A (en) Power amount prediction program, power amount prediction device and power amount prediction method
US9958921B2 (en) Power management to change power limits based on device skin temperature
US20160266629A1 (en) Changing power limits based on device state
US20150013748A1 (en) Maximum power point tracking (mppt)
TW201520784A (en) Renewable energy power generation prediction system and method and power allocation system
JP2017200360A (en) Photovoltaic power generation amount prediction system, photovoltaic power generation amount prediction method, and program
JP5734936B2 (en) Power generation amount prediction apparatus and method
JP2015104171A (en) Power demand prediction device, power demand prediction method, and power demand prediction program
JP2012010508A (en) Power generation amount prediction device, power generation amount prediction system, power generation amount prediction method, and power generation amount prediction program
JP2017010111A (en) Device, method, and program for assisting with noise data rejection
WO2014117659A1 (en) Predicting a generated power of a photovoltaic power generation device
JP6884318B2 (en) Meteorological prediction system, meteorological prediction method, meteorological prediction program
CN114221353B (en) Power grid regulation method, device, equipment, storage medium and computer program product
KR20190024162A (en) Method and apparatus for estimating generated photovoltaic power
JP2017011780A (en) Power generation amount prediction control device, power generation amount prediction method and power generation amount prediction program
JP5447295B2 (en) Power generation amount estimation system
JP6893323B2 (en) Power generation equipment information estimation system and power generation equipment information estimation method
AU2021238365A1 (en) Systems and methods for enhanced reactive power management in a hybrid environment
JP2019091335A (en) Power generation amount prediction device, power generation amount prediction system, power generation amount prediction method, and power generation amount prediction program
KR20150074462A (en) Method for predict a generating energy of the light of the sun
US20200321781A1 (en) Method, server, program, and terminal for generating information relating to electrical output of solar generator
JP6452854B2 (en) Power generation amount calculation device, power generation amount calculation system, power generation amount calculation method, and program
CN104281740A (en) Method for acquiring satellite surface potential on basis of non-uniform grid division
KR20200112790A (en) Method and apparatus for estimating generated photovoltaic power

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190409

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200331

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200528

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200914

R150 Certificate of patent or registration of utility model

Ref document number: 6772625

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150