JP2011142790A - Prediction system of solar power generation amount and prediction method of solar power generation amount - Google Patents

Prediction system of solar power generation amount and prediction method of solar power generation amount Download PDF

Info

Publication number
JP2011142790A
JP2011142790A JP2010003377A JP2010003377A JP2011142790A JP 2011142790 A JP2011142790 A JP 2011142790A JP 2010003377 A JP2010003377 A JP 2010003377A JP 2010003377 A JP2010003377 A JP 2010003377A JP 2011142790 A JP2011142790 A JP 2011142790A
Authority
JP
Japan
Prior art keywords
amount
power generation
value
solar radiation
photovoltaic power
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
JP2010003377A
Other languages
Japanese (ja)
Other versions
JP5493886B2 (en
Inventor
Tatsuya Iizaka
達也 飯坂
Yasuhisa Kanazawa
康久 金澤
Takeshi Ono
健 大野
Shinsuke Nii
真介 仁井
Naotaka Miyamura
尚孝 宮村
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.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Co 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 Fuji Electric Co Ltd filed Critical Fuji Electric Co Ltd
Priority to JP2010003377A priority Critical patent/JP5493886B2/en
Publication of JP2011142790A publication Critical patent/JP2011142790A/en
Application granted granted Critical
Publication of JP5493886B2 publication Critical patent/JP5493886B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To attain stabilization of a power system by suppressing increase in the number of capacitors and the standby amount of generators, even when a photovoltaic power generator is increased. <P>SOLUTION: In the prediction system of solar power generation amount which predicts the amount of photovoltaic power generation of photovoltaic power generators dispersed and installed in the power system, a data input means collects each pieces of data, including the predicted values and result values of atmospheric temperature and the amount of solar radiation in each region and the electrical energy of the power system. A parameter-estimating means utilizes the atmospheric temperature, the result value regarding the amount of solar radiation, and the electrical energy at substantially the same time, regarding a plurality of different days in each collected data, and computes regression coefficients, by using the atmospheric temperature and the results value regarding the amount of solar radiation as predicted variables and the electrical energy as an object variable in each or the plurality of regions. A solar power generation amount predicting means computes the predictor of the amount of solar power generation at each region from the regression coefficient at each region computed by the parameter estimating means and the predicted values of the amount of solar radiation. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、太陽光発電装置の発電量を予測する技術に係り、特に電力系統における需給計画および系統安定化制御に用いられる太陽光発電量予測システムおよび太陽光発電量予測方法に関する。   The present invention relates to a technology for predicting a power generation amount of a solar power generation apparatus, and more particularly to a solar power generation amount prediction system and a solar power generation amount prediction method used for supply and demand planning and system stabilization control in an electric power system.

電力系統では、周波数や電圧の変動を一定の範囲に抑えるために様々な制御が行われている。例えば、電力会社の管轄する系統全体において、発電量に対して需要が増減すると周波数の変化として現れるため、系統安定化用の発電機の出力を制御している。また一部の地域で発電量と需要のバランスが崩れると電圧の変化として現れるため、変圧器のタップの変更・無効電力の増減などにより電圧調整を行っている。   In the electric power system, various controls are performed in order to suppress fluctuations in frequency and voltage within a certain range. For example, in the entire system under the jurisdiction of the power company, when the demand increases or decreases with respect to the amount of power generation, it appears as a change in frequency, so the output of the generator for system stabilization is controlled. In some areas, if the balance between power generation and demand is lost, it will appear as a change in voltage. Therefore, the voltage is adjusted by changing the tap of the transformer or increasing or decreasing the reactive power.

一般に電力制御機器を停止状態から制御可能な状態にするには時間を要する。たとえば発電機はその種類によっては停止状態から待機状態にするまで1日かかる場合がある。そのために、電力会社は、気象や暦の情報から電力需要量を予測し、電力需要量に対して過不足が無いように発電機を待機させて制御が必要なときには待機状態から稼動状態へ移行させて制御をしている。   Generally, it takes time to change the power control device from a stopped state to a controllable state. For example, depending on the type of the generator, it may take one day from the stop state to the standby state. To that end, the power company predicts the power demand from weather and calendar information, makes the generator stand by so that there is no excess or deficiency in the power demand, and shifts from the standby state to the operating state when control is required. Let me control.

一方、近年、太陽光発電に代表される自然エネルギーの導入量が増加している。太陽光は、日射量により発電量が変化するので、他の発電機のように発電量を制御することができない。太陽光発電による電力が系統へ供給されると、これまでは需要の変化に追従して制御していたものを、太陽光などの制御不能な発電量まで考慮して、需要と発電量のバランスを取るように制御しなければならなくなった。このため、太陽光発電による電力の増加に伴い、太陽光発電による発電量の予測精度を向上させ、太陽光発電量の変動を十分に制御できるだけの制御機器(2次電池(NAS),分散電源(燃料電池,ガスタービン),SVC,SVRなど)を待機させることが重要になってくる。   On the other hand, in recent years, the amount of natural energy introduced as represented by solar power generation has increased. Since the power generation amount of sunlight changes depending on the amount of solar radiation, the power generation amount cannot be controlled like other generators. When power from solar power generation is supplied to the grid, the balance between demand and power generation, taking into account uncontrollable power generation, such as solar power, that was controlled according to changes in demand. Had to control to take. For this reason, with the increase in the power generated by solar power generation, the control accuracy (secondary battery (NAS), distributed power supply) that can improve the prediction accuracy of the power generation amount by solar power generation and can sufficiently control the fluctuation of the solar power generation amount. It is important to wait (fuel cell, gas turbine), SVC, SVR, etc.).

この太陽光発電量を予測する基本式として、非特許文献1には下記(1)式が示されている。
太陽光発電量=設備容量×日射量×係数 ・・・(1)
ここで、係数とは、太陽光パネルの汚れ、パネル温度の変化、インバータのロスなど様々な要因による発電量の低下要因からなる係数である。
As a basic formula for predicting the amount of photovoltaic power generation, Non-Patent Document 1 shows the following formula (1).
Photovoltaic power generation amount = Facility capacity x Solar radiation x Coefficient (1)
Here, the coefficient is a coefficient that is composed of factors that decrease the amount of power generation due to various factors such as contamination of the solar panel, changes in the panel temperature, and loss of the inverter.

しかしながら、(1)式を用いて、ある地域・地区の発電量を予測する場合、「日射量」は気象予報から得ることができるものの、「設備容量」と「係数」については分からない場合が多いので簡単に予測することができない。太陽光発電装置は各家庭に備えられていることが一般的だからである。各家庭の契約内容を調査することで、設備容量をある程度推計することも可能であるが、係数まで把握することは困難である。係数はパネルの設置向き、ビルなどの障害物の有無など様々なものに影響されるからである。   However, when using formula (1) to predict the amount of power generation in a certain region / district, the “irradiation amount” can be obtained from the weather forecast, but “equipment capacity” and “coefficient” may not be known. There are so many that it cannot be predicted easily. This is because a solar power generator is generally provided in each home. It is possible to estimate the installed capacity to some extent by investigating the contract contents of each household, but it is difficult to grasp the coefficient. This is because the coefficient is influenced by various factors such as the installation direction of the panel and the presence or absence of obstacles such as buildings.

このため、特許文献1,2では、(1)式を用いずに設備容量や係数が分からなくても太陽光発電量を予測する技術として、晴れ、曇、雨など典型的な気象条件と、過去の発電量の実績値から、該当日の発電量を予測する技術が開示されている。具体的には、過去の単位時間あたりの発電量と気象状態を制御装置単位で蓄積して、そのデータをもとに予測するというものである。   For this reason, in Patent Documents 1 and 2, as a technique for predicting the amount of photovoltaic power generation without knowing the equipment capacity and coefficient without using the equation (1), typical weather conditions such as sunny, cloudy, and rain, A technique for predicting the power generation amount on the corresponding day from the past actual power generation value is disclosed. Specifically, the amount of power generation and weather conditions per unit time in the past are accumulated in units of control devices, and prediction is based on the data.

しかしながら、現状、電力会社が把握できる情報は、変電所などの送り出し側(送電端)における、分単位の電力量(すなわち、本来の電力需要量から各家庭の太陽光発電量を減算した値)などに限られているため、電力会社が特許文献1の発明を実現するには、各家庭の発電量等の情報を分単位で収集しなければならず、多大な設備コストがかかる。   However, currently, the information that can be grasped by the power company is the amount of power in minutes on the sending side (transmission end) of substations (ie, the value obtained by subtracting the amount of solar power generated from each household from the actual amount of power demand) Therefore, in order for the electric power company to realize the invention of Patent Document 1, it is necessary to collect information such as the amount of power generation of each household in units of minutes, which requires a large facility cost.

また、実際に予測値と実績値では少なからず誤差が生ずるので、特許文献2の技術では、太陽光発電装置が普及し、その発電量が多くなると誤差を吸収するために蓄電装置の数や発電機の待機量を多くする必要がある。   In addition, since an error occurs in the predicted value and the actual value in practice, in the technique of Patent Document 2, the number of power storage devices and power generation are reduced in order to absorb the error when the amount of power generation increases. It is necessary to increase the standby amount of the machine.

特開2009−60704号公報JP 2009-60704 A 特開2006−33908号公報JP 2006-33908 A

日本工業規格JIS C 8907:2005(太陽光発電システムの発電電力量推定方法)Japanese Industrial Standard JIS C 8907: 2005 (Method for estimating the amount of power generated by a photovoltaic power generation system)

本発明は、上述のかかる事情に鑑みてなされたものであり、太陽光発電装置による発電量を効率的に予測し精度の高い需給計画を可能にすると共に、各家庭の太陽光発電装置が増加した場合でも、蓄電装置の数や発電機の待機量の増加を抑制して電力系統の安定化を実現することのできる太陽光発電量予測システムおよび太陽光発電量予測方法を提供することを目的とする。   The present invention has been made in view of the above-described circumstances, and it is possible to efficiently predict the amount of power generated by a solar power generation device and enable a highly accurate supply and demand plan, and increase the number of solar power generation devices in each household. The purpose of the present invention is to provide a photovoltaic power generation amount prediction system and a photovoltaic power generation amount prediction method capable of realizing stabilization of the power system by suppressing an increase in the number of power storage devices and generator standby amounts. And

上記目的を達成するため、本発明に係る太陽光発電量予測システムは、電力系統に分散設置された太陽光発電装置の太陽光発電量を予測する太陽光発電量予測システムであって、地域ごとの気温、日射量の予報値と実績値および電力系統の電力量の各データを収集するデータ入力手段と、前記収集した各データのうち、異なる複数の日について略同一時刻の気温と日射量実績値と電力量とを用いて、気温と日射量実績値を説明変数とし、電力量を目的変数としたときの回帰係数(後述する係数a)を一または複数の地域ごとに算出するパラメータ推定手段と、パラメータ推定手段で算出した地域ごとの回帰係数、および日射量予報値から太陽光発電量の予測値を地域ごとに算出する太陽光発電量予測手段と、を備えたことを特徴とする。   In order to achieve the above object, a photovoltaic power generation amount prediction system according to the present invention is a photovoltaic power generation amount prediction system that predicts the amount of photovoltaic power generation of photovoltaic power generation devices distributed and installed in an electric power system. Data input means for collecting each data of predicted temperature and actual value of solar radiation amount, solar radiation amount and electric energy of power system, and actual temperature and solar radiation amount results at substantially the same time for a plurality of different days among the collected data Parameter estimation means for calculating a regression coefficient (coefficient a to be described later) for one or a plurality of regions, using the value and the electric energy as an explanatory variable with the actual temperature and the amount of solar radiation as an explanatory variable and the electric energy as an objective variable And a photovoltaic power generation amount prediction means for calculating a predicted value of the photovoltaic power generation amount for each region from the regression coefficient for each region calculated by the parameter estimation means and the solar radiation amount forecast value.

本発明では、略同一時刻の気温と日射量実績値と電力量の各データを蓄積し、これら過去のデータを用いて、回帰係数(係数a)を演算する。そして、この回帰係数と日射量予報値から太陽光発電量を予測する。   In the present invention, temperature, actual solar radiation value and electric energy at substantially the same time are accumulated, and the regression coefficient (coefficient a) is calculated using these past data. Then, the photovoltaic power generation amount is predicted from the regression coefficient and the solar radiation amount forecast value.

好ましくは、電力需要が日射量変化の影響を受けにくい正午付近のデータを用いるようにするのが良い。   Preferably, it is preferable to use data around noon where the power demand is less susceptible to changes in solar radiation.

また、本発明に係る太陽光発電量予測システムは、電力系統に分散設置された太陽光発電装置の太陽光発電量を予測する太陽光発電量予測システムであって、地域ごとの日射量の予報値と実績値および電力系統の電力量の各データを収集するデータ入力手段と、前記収集した各データを用いて、日射量実績値の単位時間の変化量を説明変数とし、電力量の単位時間の変化量を目的変数としたときの回帰係数を一または複数の地域ごとに算出するパラメータ推定手段と、パラメータ推定手段で算出した地域ごとの回帰係数、および日射量予報値から太陽光発電量の予測値を地域ごとに算出する太陽光発電量予測手段と、を備えたことを特徴とする。   The solar power generation amount prediction system according to the present invention is a solar power generation amount prediction system for predicting the solar power generation amount of the solar power generation devices distributed and installed in the power system, and forecasts the solar radiation amount for each region. Data input means for collecting each value data, actual value, and electric energy of the power system, and using the collected data, the amount of change in unit time of solar radiation actual value is an explanatory variable, and unit time of electric energy The parameter estimation means for calculating the regression coefficient for each region or regions with the amount of change in the objective variable as the objective variable, the regression coefficient for each region calculated by the parameter estimation means, and the solar radiation forecast value And a photovoltaic power generation amount prediction means for calculating a predicted value for each region.

本発明では、日射量と電力量の変化分から回帰係数(係数a)を演算して、この回帰係数を用いて太陽光発電量を予測する。なお、単位時間は、例えば1分単位など任意の時間で良いが、系統安定化用の制御装置の応答時間に合わせるようにするのが好ましい。   In the present invention, a regression coefficient (coefficient a) is calculated from the amount of solar radiation and the amount of power, and the photovoltaic power generation amount is predicted using this regression coefficient. The unit time may be an arbitrary time such as a minute unit, but it is preferable to match the response time of the control device for stabilizing the system.

さらに、本発明に係る太陽光発電量予測システムの太陽光発電量予測手段は、回帰係数の信頼区間を演算すると共にデータ入力手段によって収集した日射量の予報値と実績値から日射量予報値誤差の信頼区間を演算し、当該回帰係数の信頼区間および当該日射量予報値誤差の信頼区間をもとに太陽光発電量の上限値と下限値とを含む予測値を算出することを特徴とする。   Further, the photovoltaic power generation amount prediction means of the photovoltaic power generation amount prediction system according to the present invention calculates the confidence interval of the regression coefficient, and calculates the solar radiation amount forecast value error from the predicted value and the actual value of the solar radiation amount collected by the data input means. And calculating a prediction value including an upper limit value and a lower limit value of the amount of photovoltaic power generation based on the confidence interval of the regression coefficient and the confidence interval of the solar radiation amount forecast value error. .

本発明では、予測値を信頼区間の上下限値として算出することにより、電力の需給計画や制御装置の出力制御の精度を向上させることができる。   In the present invention, by calculating the predicted value as the upper and lower limit values of the confidence interval, it is possible to improve the accuracy of power supply and demand planning and output control of the control device.

また、本発明に係る太陽光発電量予測方法は、電力系統に分散設置された太陽光発電装置の発電量の実績値を計測せずに太陽光発電量を予測する太陽光発電量予測方法であって、地域ごとの日射量の予報値と実績値および電力系統の送電端の電力量の各データを収集するデータ入力段階と、前記収集した各データを用いて、日射量と電力量の単位時間の変化量の実績値から日射量を説明変数とし、電力量を目的変数としたときの回帰係数を地域ごとに算出するパラメータ推定段階と、任意の確率の前記回帰係数の信頼区間および日射量予報値の信頼区間を算出すると共に、当該算出した信頼区間の値、パラメータ推定段階で算出した地域ごとの回帰係数、および日射量予報値から太陽光発電量の上限値と下限値とを含む予測値を地域ごとに演算する太陽光発電量予測段階と、地域ごとに当該地域の予測値を当該地域を管轄する系統安定化用の制御装置へ送信する予測値送信段階と、を備えたことを特徴とする。   Moreover, the photovoltaic power generation amount prediction method according to the present invention is a photovoltaic power generation amount prediction method for predicting the amount of photovoltaic power generation without measuring the actual value of the amount of power generation of the photovoltaic power generation devices distributed and installed in the power system. A data input stage for collecting the forecasted and actual values of solar radiation amount for each region and the power amount data at the transmission end of the power system, and using the collected data, the unit of solar radiation amount and power amount A parameter estimation stage that calculates the regression coefficient for each region using the amount of solar radiation as the explanatory variable and the amount of power as the objective variable from the actual value of the amount of change in time, the confidence interval of the regression coefficient of arbitrary probability, and the amount of solar radiation A prediction interval that includes the confidence interval of the forecast value, and includes the calculated confidence interval value, the regression coefficient for each region calculated in the parameter estimation stage, and the upper and lower limits of photovoltaic power generation from the solar radiation forecast value Show values by region A photovoltaic power generation amount prediction step of, characterized by comprising prediction value transmission step of transmitting the predicted value of the region for each region to the controller for system stabilization having jurisdiction over the area, the.

また、本発明に係る太陽光発電量予測方法は、太陽光発電量予測段階で演算した太陽光発電量の予測値に基づいて、制御装置の起動または停止の判定、並びに、制御装置の出力を決定することを特徴とする。   Further, the photovoltaic power generation amount prediction method according to the present invention is based on the predicted value of the photovoltaic power generation amount calculated in the photovoltaic power generation amount prediction stage, and determines whether to start or stop the control device, and outputs the control device. It is characterized by determining.

本発明では、予測した上下限値によって制御装置をの起動・停止の判定や出力を決定するので、制御装置の効率的な運転が可能になる。   In the present invention, since the determination of the start / stop of the control device and the output are determined based on the predicted upper and lower limit values, the control device can be efficiently operated.

本発明では、太陽光発電量予測の基本式のうち、設備容量と係数を一体として扱って地域ごとに太陽光発電量を予測することができるので、電力会社は各家庭に備えられた太陽光発電装置の発電量や係数、設備容量を取得する必要がなく、作業効率の向上、設備コストの削減を図ることができる。また、地域ごとに太陽光発電量を上下限の幅を持って予測することにより、制御装置を効率的に運転することができる。   In the present invention, among the basic formulas for predicting the amount of photovoltaic power generation, it is possible to predict the amount of photovoltaic power generation for each region by treating the installed capacity and coefficient as a unit. It is not necessary to acquire the power generation amount, coefficient, and equipment capacity of the power generation device, and work efficiency can be improved and equipment costs can be reduced. In addition, the control device can be operated efficiently by predicting the amount of photovoltaic power generation with a range of upper and lower limits for each region.

本発明の一実施の形態による電力系統を含む太陽光発電量予測システムの構成図である。1 is a configuration diagram of a photovoltaic power generation amount prediction system including an electric power system according to an embodiment of the present invention. 図1の太陽光発電量予測システムの処理概要を示すフローチャートである。It is a flowchart which shows the process outline | summary of the solar energy generation amount prediction system of FIG. 本発明の実施の形態による太陽光発電量の推定の考え方の説明図である。It is explanatory drawing of the view of the estimation of the amount of photovoltaic power generation by embodiment of this invention. 図1のパラメータ推定手段の係数aの算出に用いるデータの選別方法の概念図である。It is a conceptual diagram of the selection method of the data used for calculation of the coefficient a of the parameter estimation means of FIG. 本発明の実施の形態による方法1のパラメータ推定処理に用いられる入力データファイル31の一例である。It is an example of the input data file 31 used for the parameter estimation process of the method 1 by embodiment of this invention. 気温、天候、電力量の関係を示す説明図である。It is explanatory drawing which shows the relationship between temperature, the weather, and electric energy. 本発明の実施の形態による方法2のパラメータ推定処理に用いられる入力データファイル31の一例である。It is an example of the input data file 31 used for the parameter estimation process of the method 2 by embodiment of this invention. 日射量の誤差分布(図8(a))と係数aの誤差分布(図8(b))の説明図である。It is explanatory drawing of the error distribution (FIG.8 (a)) of the solar radiation amount, and the error distribution (FIG.8 (b)) of the coefficient a. 日射量誤差の信頼区間の説明図であり、図9(a)は天気と日射量誤差との関係の説明図、図9(b)は日ごとの日射量誤差から信頼区間の上下限値を求める方法の説明図である。FIG. 9A is an explanatory diagram of the confidence interval of the solar radiation error, FIG. 9A is an explanatory diagram of the relationship between the weather and the solar radiation error, and FIG. 9B is the upper and lower limits of the confidence interval from the daily solar radiation error. It is explanatory drawing of the method to obtain | require. 太陽光発電量予測処理の係数aの信頼区間の説明図であり、図10(a)は、パラメータ推定処理の結果求めた係数aを含むデータ例であり、図10(b)は太陽光発電量予測処理の係数aの信頼区間の上下限値の求め方の説明図である。It is explanatory drawing of the confidence interval of the coefficient a of photovoltaic power generation amount prediction processing, Fig.10 (a) is a data example containing the coefficient a calculated | required as a result of parameter estimation processing, FIG.10 (b) is photovoltaic power generation. It is explanatory drawing of how to obtain | require the upper and lower limit value of the confidence interval of the coefficient a of quantity prediction processing. 図1のパラメータファイル32に保存される係数aの信頼区間のデータの一例である。It is an example of the data of the confidence interval of the coefficient a preserve | saved at the parameter file 32 of FIG. 図1のパラメータファイル32に保存される日射量予測値の信頼区間のデータの一例である。It is an example of the data of the confidence interval of the solar radiation amount estimated value preserve | saved at the parameter file 32 of FIG. 太陽光発電量の上下限値の予測値から需給計画および制御装置に対する制御量の求め方の一例の説明図である。It is explanatory drawing of an example of how to obtain | require the control amount with respect to a supply-and-demand plan and a control apparatus from the predicted value of the upper and lower limit value of photovoltaic power generation. 図1の出力データファイル33の一例である。It is an example of the output data file 33 of FIG. 図1の太陽光発電量予測手段23の日ごとの太陽光発電量予測値(上下限値)算出ルーチンの処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of the daily photovoltaic power generation amount prediction value (upper / lower limit value) calculation routine of the solar power generation amount prediction means of FIG. 図1の太陽光発電量予測手段23の毎時の太陽光発電量予測値算出ルーチンの処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of the hourly photovoltaic power generation amount predicted value calculation routine of the photovoltaic power generation amount prediction means of FIG.

以下、本発明の実施の形態を説明する。図1は本発明の実施の形態による太陽光発電量予測システムの構成と同システムが適用される電力系統の概要を表した図である。   Embodiments of the present invention will be described below. FIG. 1 is a diagram showing a configuration of a photovoltaic power generation amount prediction system according to an embodiment of the present invention and an outline of a power system to which the system is applied.

この図において、発電所5で発電された電力を電力系統(変電設備を含む。)を通して電力需要家である各家庭3へ供給している。また、各家庭3には、太陽光発電装置4が設置されており、普段は自宅の消費電力の一部または全部を賄う一方、余剰電力は電力系統へ供給する。このように、各家庭は、電力需要家であると共に電力供給者でもある。なお、本実施の形態では、一般の家庭に太陽光発電装置4が設けられた場合を例に説明するが、工場などの事業者であっても良い。   In this figure, the electric power generated at the power plant 5 is supplied to each household 3 that is an electric power consumer through an electric power system (including substation equipment). Each home 3 is provided with a solar power generation device 4 that usually covers part or all of the power consumption of the home while supplying surplus power to the power system. Thus, each household is a power consumer as well as a power supplier. In the present embodiment, a case where the solar power generation device 4 is provided in a general household will be described as an example, but an operator such as a factory may be used.

また、電力系統は、地域ごとに区分され、各地域には、電力安定化のための制御装置2、温度計7、日射量計8が設けられている。制御装置2は、周波数・電圧安定化用の発電機あるいは2次電池(NAS)や燃料電池であり、電力系統と接続されている。   The electric power system is divided into regions, and a control device 2, a thermometer 7, and a solar radiation meter 8 for power stabilization are provided in each region. The control device 2 is a generator for frequency / voltage stabilization, a secondary battery (NAS), or a fuel cell, and is connected to the power system.

さらに、温度計7、日射量計8も地域ごとに設置され、それぞれ気温、日射量を測定している。なお、 制御装置2が設置される地域区分と温度計7、日射量計8が設置される地域区分は同じでも良いし、異なっていても良い。   Furthermore, the thermometer 7 and the solar radiation meter 8 are also installed for each region, and measure the temperature and the solar radiation amount, respectively. In addition, the area division in which the control apparatus 2 is installed and the area division in which the thermometer 7 and the solar radiation meter 8 are installed may be the same or different.

太陽光発電量予測システム1には、電力系統の変電所の送電端に設けられた電力量計6から電力量データを収集し、また、各地域の温度計7、日射量計8から気温データ、日射量データを収集している。なお、電力量計6は電力量を直接収集する機器に限らず、電力量を演算可能な電気データを収集する機器も含む。例えば、電力系統に備えられたPT,CTを介して収集された電気データによって電力量を演算することができるので、PT,CTも電力量計に含まれる。   The photovoltaic power generation amount prediction system 1 collects power amount data from a watt hour meter 6 provided at the transmission end of a substation in the power system, and also records temperature data from a thermometer 7 and a solar radiation amount meter 8 in each region. , Collecting solar radiation data. The watt-hour meter 6 is not limited to a device that directly collects the amount of power, but also includes a device that collects electrical data that can calculate the amount of power. For example, since the electric energy can be calculated from the electrical data collected via PT and CT provided in the electric power system, PT and CT are also included in the watt-hour meter.

また、気象業者のサーバ9からは、気象予報や気温、日射量の予報値を逐次入力している。これらのデータの収集は、電力量計6、温度計7、日射量計8に通信機能を設け、データ通信によって収集しても良いし、専用のデータ収集装置(図示せず)を介して収集するようにしても良い。   In addition, weather forecasts, temperature values, and solar radiation forecast values are sequentially input from the weather trader's server 9. These data can be collected by providing a communication function to the watt-hour meter 6, the thermometer 7, and the solar radiation meter 8 and collecting the data by data communication or via a dedicated data collection device (not shown). You may make it do.

太陽光発電量予測システム1は、電力量計6、温度計7、日射量計8、気象業者のサーバ9とデータ通信を行う通信装置11、データの演算処理を実行する演算処理部12、データを記憶する記憶部13から構成されている。記憶部13と演算処理部12は、それぞれサーバ構成にすることも可能である。     The photovoltaic power generation amount prediction system 1 includes a watt-hour meter 6, a thermometer 7, a solar radiation meter 8, a communication device 11 that performs data communication with a server 9 of a weather trader, an arithmetic processing unit 12 that executes data arithmetic processing, and data It is comprised from the memory | storage part 13 which memorize | stores. The storage unit 13 and the arithmetic processing unit 12 can each be configured as a server.

また、演算処理部12は、通信装置11から渡された受信データを記憶部13に保存するデータ入力手段21、受信した各種データをもとに太陽光発電量予測のためのパラメータを演算するパラメータ推定手段22、演算したパラメータを用いて太陽光発電量の予測値を演算する太陽光発電量予測手段23、予測値を制御装置2へ送信する予測値送信手段24を備えている。各手段21〜24は、プログラムとしてCPUによって実行可能な機能である。   Further, the arithmetic processing unit 12 is a data input means 21 for storing the received data passed from the communication device 11 in the storage unit 13, and a parameter for calculating a parameter for photovoltaic power generation amount prediction based on the received various data. An estimation unit 22, a photovoltaic power generation amount prediction unit 23 that calculates a predicted value of the photovoltaic power generation amount using the calculated parameters, and a predicted value transmission unit 24 that transmits the prediction value to the control device 2 are provided. Each means 21 to 24 is a function that can be executed by the CPU as a program.

なお、気温データ、日射量データは、各地域の温度計7、日射量計8から直接通信回線を介して入力しても良いし、気象業者から入力するようにしても良い。   The temperature data and solar radiation data may be input directly from the local thermometer 7 and solar radiation meter 8 via a communication line, or may be input from a weather company.

次に、上記の構成を有する太陽光発電量予測システム1の動作を説明する。
太陽光発電量予測システム1は、電力系統に送り出されている電力を電力量計6を介して分・秒単位で計測している。電力量計6から送られてくる電力量データは、通信装置11で受信され、演算処理部12のデータ入力手段21によって記憶部13の入力データファイル31に保存される。
Next, the operation of the photovoltaic power generation amount prediction system 1 having the above configuration will be described.
The photovoltaic power generation amount prediction system 1 measures the electric power sent to the electric power system in units of minutes and seconds via the watt-hour meter 6. The power amount data sent from the watt hour meter 6 is received by the communication device 11 and stored in the input data file 31 of the storage unit 13 by the data input means 21 of the arithmetic processing unit 12.

ここで、データ入力手段21を介して収集される電力量は、下記(2)式のように、各家庭で消費されている電力需要量から、各家庭の太陽光発電装置による発電量(以下、「太陽光発電量」という。)を差し引いたものである。
電力量=電力需要量−太陽光発電量 ・・・(2)
Here, the amount of power collected via the data input means 21 is the amount of power generated by the solar power generation device of each home (hereinafter referred to as the amount of power consumed by each home) as shown in the following equation (2). , Referred to as “photovoltaic power generation amount”).
Electricity amount = Electricity demand amount-Solar power generation amount (2)

一方、各地域の温度計7、日射量計8からも同様に分・秒単位で気温、日射量を計測し、計測したデータは通信装置11で受信され、データ入力手段21によって記憶部13の入力データファイル31に保存される。   On the other hand, the temperature and solar radiation amount are similarly measured in units of minutes and seconds from the thermometer 7 and solar radiation meter 8 in each region, and the measured data is received by the communication device 11 and is stored in the storage unit 13 by the data input means 21. It is stored in the input data file 31.

なお、日射量データについては、日射量計8から入力することに代えて、ある一家庭における太陽光発電装置4の係数と設備容量を太陽光発電量予測システム1の記憶部13に保存しておき、その家庭の太陽光発電量を分や秒単位に計測して、(1)式から日射量を計算するようにしても良い。   In addition, about solar radiation amount data, it replaces with inputting from the solar radiation meter 8, and the coefficient and installation capacity of the solar power generation device 4 in one certain household are preserve | saved in the memory | storage part 13 of the solar power generation amount prediction system 1. Alternatively, the solar power generation amount of the household may be measured in minutes or seconds, and the solar radiation amount may be calculated from equation (1).

太陽光発電量予測システム1は、記憶部13の入力データファイル31に蓄積されているデータを用いて、太陽光発電量の予測を行う。その処理手順は、図2に示すように、パラメータ推定手段22によって実行される太陽光発電量予測を行うためのパラメータを推定する処理(以下、「パラメータ推定処理」という。)と、太陽光発電量予測手段23によって実行される太陽光発電量を予測する処理(以下、「太陽光発電量予測処理」という。)からなる。   The photovoltaic power generation amount prediction system 1 uses the data stored in the input data file 31 of the storage unit 13 to predict the photovoltaic power generation amount. As shown in FIG. 2, the processing procedure includes a process for estimating a solar power generation amount prediction executed by the parameter estimation unit 22 (hereinafter referred to as “parameter estimation process”), and solar power generation. It consists of the process (henceforth "solar power generation prediction process") which estimates the amount of photovoltaic power generation performed by the amount prediction means 23.

なお、気象予報値は、気象業者のサーバ9から配信されるが、ここでは気象予報値として日射量の予報値が配信されるものとする。   In addition, although a weather forecast value is delivered from the server 9 of a weather company, the forecast value of solar radiation amount shall be delivered as a weather forecast value here.

<ステップS1:パラメータ推定処理>
次に、太陽光発電量予測システム1のパラメータ推定手段22の処理内容について説明する。
<Step S1: Parameter estimation processing>
Next, the processing content of the parameter estimation means 22 of the photovoltaic power generation amount prediction system 1 will be described.

パラメータ推定手段22は、オペレータの起動要求や事前に設定した任意のタイミングによって起動され、まず、一定期間蓄積された日射量データを用いて日射量データが瞬間的に変化している時間帯を切り出す(ステップS1a)。   The parameter estimation unit 22 is activated by an operator activation request or an arbitrary timing set in advance. First, the solar radiation amount data accumulated for a certain period is used to cut out a time zone in which the solar radiation amount data is instantaneously changing. (Step S1a).

切り出す条件は、単位時間当たりの日射量が一定値以上変化した場合、ある一定比率変化した場合、あるいは、ある一定周波数成分以上などがある。この切り出し条件でデータを抽出する。   The cutting conditions include a case where the amount of solar radiation per unit time changes by a certain value or more, a certain ratio changes, or a certain frequency component or more. Data is extracted under this cutting condition.

パラメータ推定手段22は、次に、ステップS1aで切り出された時間帯における日射量の変化と太陽光発電量の変化の関係を演算する(ステップS1b)。このとき、太陽光発電の特性式(1)を利用しようとすると、設備容量と係数を求めなければならない。しかし、電力会社が家庭ごとの太陽光発電量データを収集するには多大なコストがかかるため現実的ではない。このため、本実施の形態では、太陽光発電量は日射量の変化に対して極めて速く追従する一方、電力需要量は日射量が急変してもすぐには変化しないという物理的特性を利用する。すなわち、図3(c)のグラフで示されるトータルの電力需要量は、図3(b)で示される日中時々刻々変化する太陽光発電量と図3(a)で示される電力会社で発電される電力量(すなわち電力量計6で計測した電力量)の和であると捉え、図3(a)の計測値をもとに図3(b)の太陽光発電量を推測するというものである。なお、日射量が小さくなると暗くなるため、電灯が点けられ電力需要量が変化する場合があるが、多くの場合電灯の点滅は人間が操作するため瞬時に変化するものではなく、その操作による電力変化もある時間範囲の中で平均化される。   Next, the parameter estimation means 22 calculates the relationship between the change in the amount of solar radiation and the change in the amount of photovoltaic power generation in the time period cut out in step S1a (step S1b). At this time, if the characteristic formula (1) of photovoltaic power generation is to be used, the equipment capacity and coefficient must be obtained. However, it is not realistic for an electric power company to collect the amount of photovoltaic power generation data for each household because of the great cost. For this reason, in this embodiment, the amount of photovoltaic power generation follows the change in the amount of solar radiation very quickly, while the demand for power does not change immediately even if the amount of solar radiation changes suddenly. . That is, the total power demand shown in the graph of FIG. 3 (c) is the amount of photovoltaic power generation that changes momentarily during the day shown in FIG. 3 (b) and generated by the electric power company shown in FIG. 3 (a). It is assumed that it is the sum of the amount of electric power (that is, the amount of electric power measured by the electric energy meter 6), and the amount of photovoltaic power generation in FIG. 3 (b) is estimated based on the measured value of FIG. 3 (a) It is. In addition, since it becomes dark when the amount of solar radiation becomes small, there is a case where the light is turned on and the power demand amount changes, but in many cases the blinking of the light does not change instantaneously because it is operated by a human, and the power generated by that operation Changes are also averaged over a time range.

パラメータの推定方法としては、後述するように日射量と電力量の変化分を測定することによって推定する方法(方法1)と所定の時刻(たとえば正午)の電力量と気温、日射量とを測定することによって推定する方法(方法2)がある。方法2は、方法1に対して比較的大規模な電力系統に対して効果がある。   As a parameter estimation method, as will be described later, a method (method 1) for estimating by measuring changes in the amount of solar radiation and the amount of power, and the amount of power, temperature, and amount of solar radiation at a predetermined time (for example, noon) are measured. There is a method (method 2) of estimating by doing. Method 2 is more effective than method 1 for a relatively large power system.

(方法1.日射量と電力量の変化分の測定による推定)
具体的には、ステップS1aで切り出された時間のデータから最小二乗法で、次の(3)式の係数aを求める。
△太陽光発電量=係数a×△日射量 ・・・・(3)
ここで、係数a=「設備容量」×「係数」。また、△太陽光発電量は、太陽光発電量について前回収集したデータに対する変化量(以下、「変化分」という。)、△日射量も同様に前回サンプリング値に対する変化分を意味する。
(Method 1. Estimation by measuring changes in solar radiation and electric energy)
Specifically, the coefficient a of the following equation (3) is obtained from the time data cut out in step S1a by the least square method.
△ Solar power generation amount = Coefficient a × △ Solar radiation amount (3)
Here, coefficient a = “equipment capacity” × “coefficient”. Further, Δ solar power generation amount means a change amount (hereinafter referred to as “change amount”) of the solar power generation amount with respect to previously collected data, and Δ solar radiation amount also means a change amount with respect to the previous sampling value.

上記(3)式の係数aは、(1)式の「設備容量」と「係数」の積に相当するものである。従来のように「設備容量」と「係数」とを別々に把握する必要はない。   The coefficient a in the equation (3) corresponds to the product of “equipment capacity” and “coefficient” in the equation (1). There is no need to separately grasp “equipment capacity” and “coefficient” as in the past.

図4は、係数aの算出に用いるデータの選別方法の概念図である。この図において、図4(a1)は、電力会社による発電量(即ち、電力量計6で計測された電力量)の一日の時間変化のグラフ、図4(a2)は、図4(a1)と同じ日の日射量の時間変化のグラフである。通常1分ごとにサンプリングするがサンプリング間隔はこれに限らない。複数日についてデータを採取する。   FIG. 4 is a conceptual diagram of a data selection method used for calculating the coefficient a. In this figure, FIG. 4 (a1) is a graph of a daily time change in the amount of power generated by an electric power company (that is, the amount of power measured by the watt hour meter 6), and FIG. 4 (a2) is a graph of FIG. 4 (a1). It is a graph of the time change of the solar radiation amount on the same day. Usually, sampling is performed every minute, but the sampling interval is not limited to this. Collect data for multiple days.

次に、これらのデータから、図4(b1)、図4(b2)のように、それぞれ電力量、日射量について、変化分を計算する。ここで、Δ電力量=−Δ太陽光発電量の関係にある。そして、日射量の変化分が予め定めた一定値以上になったときの太陽光発電量の変化分を抽出する。図4(b2)の例では、日射量変化分が閾値C1以上になった時点の太陽光発電量変化分と、日射量変化が閾値C2以下になったときの太陽光発電量変化分を図4(b1)のデータから抽出する。   Next, from these data, as shown in FIG. 4 (b1) and FIG. 4 (b2), changes are calculated for the electric energy and the solar radiation amount, respectively. Here, there is a relationship of Δ power amount = −Δ solar power generation amount. Then, the amount of change in the amount of solar power generated when the amount of change in the amount of solar radiation exceeds a predetermined value is extracted. In the example of FIG. 4 (b2), the amount of change in solar power generation when the amount of change in solar radiation becomes greater than or equal to the threshold C1 and the amount of change in solar power generation when the change in solar radiation becomes less than or equal to the threshold C2 are shown. Extracted from the data of 4 (b1).

この抽出処理により、図4(c)に示すように、日射量変化分を横軸、発電量変化分を縦軸にとったときに、抽出されたデータは、第一象限と第三象限にある程度密集して出現するので、最小二乗法等の手法により、直線近似式の係数aを算出する。これにより、日射量変化と相関のある太陽光発電の変化量を処理負荷を抑えて簡便に求めることができる。   With this extraction process, as shown in FIG. 4C, when the amount of solar radiation change is taken on the horizontal axis and the amount of power generation change is taken on the vertical axis, the extracted data is in the first and third quadrants. Since they appear densely to some extent, the coefficient a of the linear approximation formula is calculated by a method such as the least square method. Thereby, the variation | change_quantity of the photovoltaic power generation correlated with the solar radiation amount change can be calculated | required simply, suppressing processing load.

なお、系統範囲が比較的大きい場合は、日射量計を複数設置する必要があり、日射量計を設置した予め定めた一定の領域ごとに、次式のように右辺を複数項にして多変量解析等の手法により、地域ごとの係数a1,a2,・・・を算出する。
△太陽光発電量=係数a1×△日射量1+係数a2×△日射量2+・・
・・・(3a)
ここで、係数ai、△日射量iは、それぞれ地域区分i(iは自然数)の係数a、日射量変化分を意味する。また、△太陽光発電量は、太陽光発電量の変化分を意味する。
In addition, when the system range is relatively large, it is necessary to install multiple solar radiation meters.For each of the predetermined areas where the solar radiation meters are installed, multiple variables with the right side as multiple terms The coefficients a1, a2,... For each region are calculated by a method such as analysis.
△ Solar power generation = coefficient a1 x △ solar radiation 1 + coefficient a2 x △ solar radiation 2 +
... (3a)
Here, the coefficient a i and the Δ solar radiation amount i mean the coefficient a and the amount of solar radiation change in the region division i (i is a natural number), respectively. Moreover, Δ solar power generation amount means a change in solar power generation amount.

パラメータ推定手段22は、算出した係数aiを記憶部13のパラメータファイル32に保存する。以下の説明において、特に指定の無い限り、係数aは、任意の地域区分iの係数aiを意味するものとする。   The parameter estimation unit 22 stores the calculated coefficient a i in the parameter file 32 of the storage unit 13. In the following description, unless otherwise specified, the coefficient a means a coefficient ai of an arbitrary regional section i.

次に、簡単な数値例で上記のパラメータ推定処理を説明する。図5は、1分ごとに入力したデータを蓄積した入力データファイル31の例を示す。この図において、電力量は、電力会社で計測した電力系統への送り出し側の電力量である。日射量は、該当地域に設置されている日射量である。ここで、Δ電力量は1分前の電力量からの変化分を意味する。なお、Δ太陽光発電量(すなわち、1分前の太陽光発電量からの変化分)は、上述したように概ね次式であるとして捉えることができる。
Δ太陽光発電量 = −Δ電力量
したがって、電力変化分のデータを収集することによって、太陽光発電量の変化量を把握することができる。
Next, the parameter estimation process will be described with a simple numerical example. FIG. 5 shows an example of an input data file 31 in which data input every minute is accumulated. In this figure, the electric energy is the electric energy on the sending side to the electric power system measured by the electric power company. The amount of solar radiation is the amount of solar radiation installed in the area. Here, the Δ electric energy means a change from the electric energy one minute ago. Note that the Δ solar power generation amount (that is, the change from the solar power generation amount one minute ago) can be generally regarded as the following equation as described above.
ΔPhotovoltaic power generation amount = −ΔPower amount Therefore, the amount of change in the solar power generation amount can be grasped by collecting data corresponding to the power change.

本実施の形態では、Δ日射量が10W/m以上のものだけを抽出するとすれば、図5の例では、12:02、および12:03の行のデータが抽出されることになる。もちろん、Δ日射量が10以上変化した他の時間のデータがあれば、それも同時に抽出する。
そして、抽出したデータから最小二乗法により係数aを求める。図5のデータでは、係数aは、約1になる。
In the present embodiment, if only those having a solar radiation amount of 10 W / m 2 or more are extracted, in the example of FIG. 5, the data of the rows of 12:02 and 12:03 are extracted. Of course, if there is data of other time when the Δ solar radiation amount has changed by 10 or more, it is also extracted at the same time.
Then, the coefficient a is obtained from the extracted data by the least square method. In the data of FIG. 5, the coefficient a is about 1.

なお、(1)式の「係数」や(3)式の「係数a」は、太陽光パネルの表面温度にも依存し、その温度は気温から推測することが可能であるが、このためには、個々の家庭に設置されている太陽光発電装置の特性を調査する必要がある。しかしながら、次の方法2で説明する手法によれば、この調査をすることなく気温の要素を考慮した係数aを算出することができる。   The “coefficient” in equation (1) and “coefficient a” in equation (3) also depend on the surface temperature of the solar panel, and the temperature can be estimated from the air temperature. Therefore, it is necessary to investigate the characteristics of the photovoltaic power generation devices installed in individual households. However, according to the method described in the following method 2, the coefficient a considering the temperature element can be calculated without conducting this investigation.

(方法2.所定時刻の太陽光発電の発電量による推定)
太陽光発電装置による時間帯ごとの発電量は、電力需要量が気温と天候に依存する性質を利用して推定する。図6は、天気をパラメータにしたときの、気温と電力量(電力量計6で計測される電力量)との関係を表した図である。気温ごとに比較したときに、晴れのときは、太陽光発電量が多いため、電力会社から見た電力量は小さくなる。一方、雨のときは太陽光発電量がほとんどないため、電力量が多く計測される。天気は日射量に影響するため、電力量、気温、日射量は下記の(4)式の関係になる。そして、同じ時間帯(例えば12時)のデータを用いて、次式から係数aを求めることができる。
電力量=f(気温)−係数a×日射量・・・(4)
ここで、f(気温)は、通常は2次式であるが、気温が高い領域、低い領域など範囲を区切って1次式にしても良い。また、気温だけでなく湿度、暦情報など他の因子を含んでもよい。用いるデータは、分・秒単位のデータでなく、1時間値の方がばらつきの影響が少なく望ましい。
(Method 2. Estimation based on the amount of photovoltaic power generation at a given time)
The amount of power generated by the solar power generation device for each time period is estimated using the property that the power demand depends on the temperature and the weather. FIG. 6 is a diagram showing the relationship between the temperature and the amount of power (the amount of power measured by the power meter 6) when the weather is used as a parameter. When compared for each temperature, when it is clear, the amount of photovoltaic power generation is large, so the amount of power viewed from the power company is small. On the other hand, when there is rain, there is almost no photovoltaic power generation, so a large amount of power is measured. Since the weather affects the amount of solar radiation, the amount of power, temperature, and amount of solar radiation are in the relationship of the following equation (4). And the coefficient a can be calculated | required from following Formula using the data of the same time slot | zone (for example, 12:00).
Electric energy = f (temperature)-coefficient a x solar radiation amount (4)
Here, f (temperature) is usually a quadratic expression, but it may be a linear expression by dividing a range such as a region where the temperature is high or a region where the temperature is low. In addition to temperature, other factors such as humidity and calendar information may be included. The data to be used is not data in minutes / seconds, but one-hour values are preferable because they are less affected by variations.

ちなみに、(4)式をより正確に記述すると、以下の式になる。
電力量=f(気温)−係数a×温度補正係数b(温度)×日射量
温度補正係数b=1−(気温ー基準温度)×係数c/100
ここで、基準温度,係数cは、太陽光パネルの種類により変化する係数であり、気温が基準温度以下の場合は補正しないものであるが、温度補正係数bの変化は,計測誤差と比較し小さいので、(4)式を用いれば実用上問題はない。
Incidentally, when the expression (4) is described more accurately, the following expression is obtained.
Electricity amount = f (air temperature) −coefficient a × temperature correction coefficient b (temperature) × irradiation amount temperature correction coefficient b = 1− (air temperature−reference temperature) × coefficient c / 100
Here, the reference temperature and the coefficient c are coefficients that vary depending on the type of solar panel, and are not corrected when the temperature is lower than the reference temperature. However, the change in the temperature correction coefficient b is compared with the measurement error. Since it is small, there is no practical problem if the equation (4) is used.

次に、(4)式を適用した具体例を簡単な数値例で説明する。図7は、正午から1分間の電力量その他の入力データである。日射量、気温を説明変数とし、電力量を目的変数とした重回帰式を構築する。図7のデータを用いると、上記(4)式は次のようになる。
電力量=(2×気温+100)−1×日射量 ・・・(5)
すなわち、係数aは1となる。
Next, a specific example to which the expression (4) is applied will be described with a simple numerical example. FIG. 7 shows the amount of power and other input data for 1 minute from noon. A multiple regression equation is constructed with the amount of solar radiation and temperature as explanatory variables and the amount of electricity as an objective variable. Using the data of FIG. 7, the above equation (4) becomes as follows.
Electricity amount = (2 × temperature + 100) −1 × irradiation amount (5)
That is, the coefficient a is 1.

この数値例は、簡単のため、一次式で構築したが、気温の項を2次式にしても良いし、湿度等他の気象データを用いてもよい。   Although this numerical example is constructed by a linear expression for simplicity, the term of temperature may be a quadratic expression, or other weather data such as humidity may be used.

なお、日射量の変化は太陽光発電量だけでなく、照明のための電力需要(照明需要)の変化から前述の電力需要量にも影響を与える。つまり係数aは、厳密には太陽光発電量だけでなく、照明需要の変化を合算した係数になっている。しかし、一般的に12時前後は太陽光発電量が多いこと、照明需要は、12時前後は少ないことから、上記(5)式の日射量の係数1は、ほぼ太陽光発電量とみなすことができる。したがって、上記(5)式は、係数aを近似しているとみなすことができる。   Note that the change in the amount of solar radiation affects not only the amount of photovoltaic power generation but also the above-described amount of power demand from the change in the power demand for lighting (lighting demand). In other words, the coefficient a is strictly a coefficient obtained by adding not only the amount of photovoltaic power generation but also the change in lighting demand. However, since the amount of photovoltaic power generation is generally large around 12:00 and the demand for lighting is small around 12:00, the coefficient of solar radiation 1 in the above equation (5) is regarded as almost the amount of photovoltaic power generation. Can do. Therefore, the above equation (5) can be regarded as approximating the coefficient a.

なお、実用的には、上記の(3a)式に対応させて、次の(4a)式によって地域ごとに係数aを求めるのが良い。
電力量=f(気温1,気温2, ・・)−係数a1×日射量1−係数a2×日射量2−・・
・・・(4a)
ここで、気温i,係数ai,日射量iは、それぞれ地域区分i(iは自然数)の気温、係数a、日射量を示す。
In practice, it is preferable to obtain the coefficient a for each region by the following equation (4a) corresponding to the above equation (3a).
Electric energy = f (temperature 1, temperature 2, ...)-coefficient a1 x solar radiation 1-coefficient a2 x solar radiation 2-
... (4a)
Here, the temperature i, the coefficient ai, and the amount of solar radiation i indicate the temperature, the coefficient a, and the amount of solar radiation of the region division i (i is a natural number), respectively.

<太陽光発電量予測処理>
次に、太陽光発電量予測手段23の処理内容について説明する。
<Solar power generation forecasting process>
Next, the processing content of the photovoltaic power generation amount prediction means 23 will be described.

太陽光発電量予測手段23は、次の(6)式をもとに地域ごとの太陽光発電量を演算する。
太陽光発電量i=係数ai×日射量i ・・・(6)
ここで、iは、地域区分の識別番号に対応している。
The photovoltaic power generation amount prediction means 23 calculates the photovoltaic power generation amount for each region based on the following equation (6).
Photovoltaic power generation amount i = coefficient ai × irradiation amount i (6)
Here, i corresponds to the identification number of the region segment.

なお、日射量は、気象業者から配信される値を用いる。この係数aiは、逐次予測結果と実績値のずれから係数を修正するようにすると良い。これにより、地域区分iの領域内で、家庭の太陽光発電装置が導入された場合でも、逐次係数aiを補正して精度の高い太陽光発電量の予測をすることができる。   In addition, the value delivered from a weather trader is used for the amount of solar radiation. The coefficient ai may be corrected from the difference between the sequential prediction result and the actual value. Thereby, even when a household solar power generation device is introduced in the region of the regional section i, it is possible to correct the sequential coefficient ai and predict the solar power generation amount with high accuracy.

次に、予測値の上下限値の求め方について説明する。
太陽光発電量予測の誤差要因は、「気象予報の誤差」と「予測モデルの誤差」(即ち、係数aの誤差)の2つに分類できる。そこで、この2つの誤差要因から、太陽光発電量予測の上下限値を求める。
Next, how to obtain the upper and lower limit values of the predicted value will be described.
The error factors in the prediction of the amount of photovoltaic power generation can be classified into two categories: “meteorological forecast error” and “prediction model error” (that is, error of coefficient a). Therefore, the upper and lower limit values of the photovoltaic power generation amount prediction are obtained from these two error factors.

最初に気象予報誤差について説明する。過去において蓄積された気象予報(日射量)と気象実績値から、その誤差(予報値−実績値)の分布を求める。日射量の誤差の分布は、図8(a)のような形状のグラフになり正規分布にならない。その理由は、晴れと予報した場合は、日射量がそれ以上になることはなく、かならず雨方向にしか誤差が発生しないためである(データ例は図9(a)参照)。   First, the weather forecast error will be described. The distribution of the error (forecast value−actual value) is obtained from the weather forecast (irradiation amount) and the actual weather value accumulated in the past. The distribution of errors in the amount of solar radiation is a graph having a shape as shown in FIG. The reason is that when it is predicted to be sunny, the amount of solar radiation does not exceed that, and an error always occurs only in the rain direction (see FIG. 9A for an example of data).

図9(a)において、複数の日の一定時刻における日射量の実績値をu、日射量の予報値をvとしたとき、日射量誤差(v−u)を算出し、この日射量誤差データを蓄積していく。   In FIG. 9 (a), when the actual value of the solar radiation amount at a fixed time on a plurality of days is u and the predicted value of the solar radiation amount is v, the solar radiation amount error (v−u) is calculated, and this solar radiation amount error data Will accumulate.

次に、この蓄積した日射量誤差データから、例えば±95%の信頼区間に入る誤差範囲を天候(晴れ、曇、雨)ごとに求める。95%信頼区間とは、誤差が下側に外れる2.5%、上側に外れる2.5%を除いた範囲のことである。もちろん、99%など他の信頼区間を求めても良い。   Next, an error range that falls within a confidence interval of, for example, ± 95% is determined for each weather (sunny, cloudy, rainy) from the accumulated solar radiation amount error data. The 95% confidence interval is a range excluding 2.5% in which the error deviates downward and 2.5% in which the error deviates upward. Of course, other confidence intervals such as 99% may be obtained.

もし、200日分の晴れデータがある場合、図9(b)に示すように日射量誤差を大きさの順番に並べて、上側5日分、下側5日分を除いた190日分が日射量誤差の95%信頼区間内のデータとなる。図9(b)の例では、日射量誤差の95%信頼区間の上側のデータは13であり、日射量誤差の95%信頼区間の下側のデータは1である。   If there is sunny data for 200 days, as shown in Fig. 9 (b), the solar radiation errors are arranged in order of magnitude, and the solar radiation for 190 days excluding the upper 5 days and the lower 5 days is solar radiation. The data is within the 95% confidence interval of the quantity error. In the example of FIG. 9B, the data on the upper side of the 95% confidence interval of the solar radiation amount error is 13, and the data on the lower side of the 95% confidence interval of the solar radiation amount error is 1.

次に、予測モデルの誤差を求める。予測モデルの誤差とは、係数aの信頼区間を意味する。上記のパラメータ推定処理の方法1または方法2で求めた係数aは、多量のサンプルデータから求めた平均的な値である。各サンプルデータごとに求めた係数aを同様に分布を求めると図8(b)に示す形状になる。信頼区間は、前述と同様に求めることができる(数値例は図10(a)参照)。   Next, the error of the prediction model is obtained. The error of the prediction model means a confidence interval of the coefficient a. The coefficient a obtained by the above-described parameter estimation processing method 1 or method 2 is an average value obtained from a large amount of sample data. When the distribution of the coefficient a obtained for each sample data is obtained in the same manner, the shape shown in FIG. The confidence interval can be obtained in the same manner as described above (see FIG. 10A for a numerical example).

図10(a)において、Δ太陽光発電量は電力量の前回値から今回値を差し引いた値、Δ日射量は日射量の今回値から前回値を差し引いた値である。そして、サンプリング時刻ごとにΔ太陽光発電量/Δ日射量、すなわち係数aが算出されている。   In FIG. 10A, Δ solar power generation amount is a value obtained by subtracting the current value from the previous value of the electric energy, and Δ solar radiation amount is a value obtained by subtracting the previous value from the current value of the solar radiation amount. Then, Δ solar power generation amount / Δ solar radiation amount, that is, coefficient a is calculated for each sampling time.

係数aの簡易的な演算方法として、もし、係数aを200個のデータから求めるとした場合、図10(b)に示すように、係数aを大きさの順番に並べて、上側5個分、下側5個分を除いた190個分が係数aの95%信頼区間内のデータとなる。図10(b)の例では、係数aの95%信頼区間の上側のデータは1.15、下側のデータは0.84である。図10(a)(b)は上記方法1によって求める係数aについて示したものであるが、方法2によって求める係数aについても同様の処理で信頼区間を求めることができる。
なお、係数aが正規分布である場合は、次の演算式によって算出することができる。
As a simple calculation method of the coefficient a, if the coefficient a is obtained from 200 pieces of data, as shown in FIG. 10 (b), the coefficients a are arranged in the order of magnitude, The data for 190 pieces excluding the lower five pieces is the data within the 95% confidence interval of the coefficient a. In the example of FIG. 10B, the upper data of the 95% confidence interval of the coefficient a is 1.15, and the lower data is 0.84. FIGS. 10A and 10B show the coefficient a obtained by the method 1, but the confidence interval can be obtained by the same process for the coefficient a obtained by the method 2.
When the coefficient a is a normal distribution, it can be calculated by the following arithmetic expression.

係数aの95%信頼区間の上側の信頼区間=係数aの平均値+係数aの標準偏差×2
係数aの95%信頼区間の下側の信頼区間=係数aの平均値-係数aの標準偏差×2
Confidence interval above 95% confidence interval of coefficient a = average value of coefficient a + standard deviation of coefficient a x 2
Lower confidence interval of 95% confidence interval of coefficient a = average value of coefficient a-standard deviation of coefficient a x 2

演算した係数aの信頼区間、日射量予報値の信頼区間は、記憶部13のパラメータファイル32に保存される。図11に係数aの信頼区間のデータ例、図12に日射量予報値の信頼区間のデータ例を示す。
上記より、太陽光発電量の予測値の上下限値とは次式より求められる。
The calculated confidence interval of the coefficient a and the confidence interval of the solar radiation amount forecast value are stored in the parameter file 32 of the storage unit 13. FIG. 11 shows a data example of the confidence interval of the coefficient a, and FIG. 12 shows a data example of the confidence interval of the solar radiation amount forecast value.
From the above, the upper and lower limit values of the predicted value of the amount of photovoltaic power generation can be obtained from the following equation.

太陽光発電量上限値
=(設備容量×係数)の上側の信頼区間×日射量予報値の上側の信頼区間
=係数aiの上側の信頼区間×(地域区分iの日射量予報値−日射量誤差の上側の信頼区間)
Upper limit of photovoltaic power generation amount = confidence interval on the upper side of (equipment capacity x coefficient) x confidence interval on the upper side of forecasted solar radiation value = confidence interval on the upper side of coefficient ai x (predicted solar radiation amount for region division i-solar radiation error) Upper confidence interval)

太陽光発電量下限値
=(設備容量×係数)の下側の信頼区間×日射量予報値の下側の信頼区間
=係数aiの下側の信頼区間×(地域区分iの日射量予報値−日射量誤差の下側の信頼区間)
Lower limit of photovoltaic power generation amount = lower confidence interval of (equipment capacity x coefficient) x lower confidence interval of forecasted solar radiation value = lower confidence interval of coefficient ai x (irradiance forecast value of region division i- Lower confidence interval for solar radiation error)

演算した太陽光発電量の上下限値は、記憶部13の出力データファイル33に保存される。   The upper and lower limit values of the calculated photovoltaic power generation amount are stored in the output data file 33 of the storage unit 13.

この太陽光発電量の予測値の上限値と下限値をもとに、少なくともこの上下限値の幅に予め定めた所定のマージンを含む発電量の変動を吸収できるように、運用する制御装置2を待機状態にする。たとえば、運用日の前日に、翌日(運用日)の太陽光発電量の上下限の予測値を算出して、この上下限値の幅の電力量を供給できるように停止状態にある制御装置2を待機状態にする。   Based on the upper limit value and the lower limit value of the predicted value of the photovoltaic power generation amount, the control device 2 that is operated so as to absorb the fluctuation of the power generation amount including a predetermined margin at least in the range of the upper and lower limit values. To the standby state. For example, on the day before the operation day, the predicted value of the upper and lower limit of the photovoltaic power generation amount of the next day (operation day) is calculated, and the control device 2 in the stopped state so that the amount of power in the range of the upper and lower limit value can be supplied. To the standby state.

また、太陽光発電量下限値に基づいて、効率は高いが出力変動が遅い火力発電機(ベースロード運転用の発電機)の発電量を決定する。具体的には、図13(a)において、太陽光発電量予測値の下限値よりも低い領域(B)の発電量を、太陽光発電を考慮しないベースロード運転用の発電機の需給計画値から差し引いく。また、太陽光発電量下限値から上限値の範囲(A)は、制御装置2、すなわち効率は高くはないが、出力変動が早い小型の発電機(分散電源)、燃料電池や2次電池(NAS)、水力発電機などの制御装置2の待機量として決定する。このように、太陽光発電量の上限値、下限値を求めることにより、全体としてエネルギー効率が高い発電計画が可能となる。   Further, based on the lower limit value of the photovoltaic power generation amount, the power generation amount of a thermal power generator (generator for base load operation) having high efficiency but slow output fluctuation is determined. Specifically, in FIG. 13 (a), the power supply amount in the region (B) lower than the lower limit value of the predicted amount of photovoltaic power generation is the planned supply and demand value of the generator for base load operation that does not consider solar power generation. Deduct from. Further, the range (A) from the lower limit value of the photovoltaic power generation amount to the upper limit value is the control device 2, that is, a small generator (distributed power source), fuel cell, or secondary battery (with low output fluctuation, although the efficiency is not high. NAS) and the standby amount of the control device 2 such as a hydroelectric generator. Thus, by obtaining the upper limit value and the lower limit value of the amount of photovoltaic power generation, a power generation plan with high energy efficiency as a whole becomes possible.

また、当日の運用状態では、図13(b)に示すように時刻tの太陽光発電量予測値の上限値y(t)からその時刻の日射量をもとに(6)式によって算定した実績値z(t)を差し引いた値を記憶部13の出力データファイル33に保存し、制御装置2の目標値として設定する。このとき(6)式で用いられる係数aiに上限と下限が存在する場合は、その中央値あるいは平均値を用いるようにする。この目標出力の演算は、制御装置が管轄する地域単位で行う。演算結果は、図14に示す出力データファイル33に保存される。この出力データファイル33は、地域区分ごとに作成され、データ入力手段21によって、当日の天気予報データ(晴れ、曇り、雨等)と時刻ごとの日射量予報値が予め入力されている。   Moreover, in the operational state of the day, as shown in FIG. 13 (b), the upper limit value y (t) of the predicted amount of photovoltaic power generation at time t was calculated from equation (6) based on the amount of solar radiation at that time. A value obtained by subtracting the actual value z (t) is stored in the output data file 33 of the storage unit 13 and set as a target value of the control device 2. At this time, when the upper limit and the lower limit exist in the coefficient ai used in the equation (6), the median value or the average value is used. The calculation of the target output is performed in units of regions that the control device has jurisdiction over. The calculation result is stored in the output data file 33 shown in FIG. This output data file 33 is created for each region, and the weather input data for the day (sunny, cloudy, rain, etc.) and the solar radiation amount forecast value for each time are input in advance by the data input means 21.

以下、上述した太陽光発電量予測手段23の処理手順を図15、図16を用いて詳述する。
太陽光発電量予測手段23の日ごとの太陽光発電量予測値(上下限値)算出ルーチンは、当日の日の出前または前日の一定時刻に起動されると、まず、パラメータファイル32にアクセスして、当日の天気に該当する日射量誤差の信頼区間データ(下限値、上限値)を抽出する(S103)。そして、出力データファイル33に保存されている各時刻の日射量予報値から日射量誤差信頼区間の上限値を減算し、減算結果を出力データファイル33の日射量予測値の下限値として保存する(S104)。同様に各時刻の日射量予報値から日射量誤差信頼区間の下限値を差し引いた値を出力データファイル33の日射量予測値の上限値として保存する(S105)。なお、天気予報データが一日のうち途中で変わる場合があるが、もし、天気予報データが変わる場合は、その切り替わり時刻から新たな天気予報データに対応する日射量誤差信頼区間データをパラメータファイル32から抽出して、このデータを用いて日射量予測値の上下限値を計算する。
Hereinafter, the processing procedure of the above-mentioned photovoltaic power generation amount prediction means 23 will be described in detail with reference to FIGS. 15 and 16.
When the solar power generation amount prediction value (upper / lower limit value) calculation routine of the solar power generation amount prediction means 23 is started at a certain time before sunrise or the previous day of the day, first, the parameter file 32 is accessed. Then, confidence interval data (lower limit value, upper limit value) of the solar radiation amount error corresponding to the weather of the day is extracted (S103). Then, the upper limit value of the solar radiation amount error confidence interval is subtracted from the solar radiation amount predicted value at each time stored in the output data file 33, and the subtraction result is stored as the lower limit value of the solar radiation amount predicted value in the output data file 33 ( S104). Similarly, a value obtained by subtracting the lower limit value of the solar radiation amount error confidence interval from the predicted solar radiation amount at each time is stored as the upper limit value of the predicted solar radiation amount in the output data file 33 (S105). Note that the weather forecast data may change during the day. If the weather forecast data changes, the parameter file 32 indicates the solar radiation amount error confidence interval data corresponding to the new weather forecast data from the switching time. And the upper and lower limit values of the predicted solar radiation amount are calculated using this data.

次に、太陽光発電量予測手段23は、日射量予測値の下限値にパラメータファイル32に保存されているその地域の係数aの信頼区間の下限値を掛け合わせて、太陽光発電量予測値の下限値を算出して、出力データファイル33に保存する(S106)。太陽光発電量予測手段23は、同様に、日射量予測値の上限値にパラメータファイル32に保存されているその地域の係数aの信頼区間の上限値を掛け合わせて、太陽光発電量予測値の上限値を算出して、出力データファイル33に保存する(S107)。以上の処理を全ての時刻、地域について実行する。   Next, the photovoltaic power generation amount predicting means 23 multiplies the lower limit value of the predicted solar radiation amount by the lower limit value of the confidence interval of the coefficient a stored in the parameter file 32, thereby calculating the predicted photovoltaic power generation amount. Is stored in the output data file 33 (S106). Similarly, the photovoltaic power generation amount predicting unit 23 multiplies the upper limit value of the predicted solar radiation amount by the upper limit value of the confidence interval of the coefficient a stored in the parameter file 32 to obtain the predicted photovoltaic power generation amount value. Is stored in the output data file 33 (S107). The above processing is executed for all times and regions.

太陽光発電量予測手段23は、上記の処理を終了すると、太陽光発電量予測値の地域ごと時刻ごとの下限値を予測値送信手段24を介して図示しない中央給電指令所へ送信する(S108)。中央給電指令所では、各地域を管轄するベースロード運転用の発電機の需給計画値(太陽光発電量を考慮しない需給計画値)から受信した太陽光発電量予測値の時刻ごとの下限値を減算して、太陽光発電量を考慮した需給計画値を算出する。この需給計画値にしたがってベースロード運転用の発電機が運転される。   When the above-described processing is completed, the photovoltaic power generation amount predicting unit 23 transmits the lower limit value for each region of the photovoltaic power generation amount predicted value for each region to the central power supply command station (not shown) via the predicted value transmitting unit 24 (S108). ). At the Central Power Supply Command Center, the lower limit value for each hour of the predicted amount of photovoltaic power generation received from the planned supply and demand values for generators for base load operation in each region (the demand and supply plan values not considering the amount of photovoltaic power generation) Subtract and calculate the supply and demand plan value considering the amount of photovoltaic power generation. The generator for base load operation is operated according to this supply and demand plan value.

その後、当日の日中、日射量実績値がデータ入力手段21を介して時々刻々入力され、入力されたデータは、出力データファイル33の日射量実績値として保存されていく。   Thereafter, during the day of the day, the actual solar radiation amount value is input from time to time via the data input means 21, and the input data is stored as the actual solar radiation amount value in the output data file 33.

太陽光発電量予測手段23の毎時の太陽光発電量予測値算出ルーチンは、各サンプリング時刻の日射量実績値の入力完了ごとに起動され、入力された日射量実績値にパラメータファイル32の係数aの信頼区間の中央値を掛け合わせて、当日の日射量実績値に基づく太陽光発電量(z)を算出して出力データファイル33に格納する(S202)。そして、太陽光発電量予測手段23は、太陽光発電量予測値の上限値(y)と上記太陽光発電量(z)の差分(y−z)を計算して、この差分値を出力データファイル33に格納する(S203)。   The hourly solar power generation amount prediction value calculation routine of the solar power generation amount prediction means 23 is started every time the input of the solar radiation amount actual value at each sampling time is completed, and the coefficient a of the parameter file 32 is added to the input solar radiation amount actual value. The solar power generation amount (z) based on the actual solar radiation amount on that day is calculated by multiplying the median values of the confidence intervals and stored in the output data file 33 (S202). And the photovoltaic power generation amount prediction means 23 calculates the difference (yz) between the upper limit value (y) of the predicted photovoltaic power generation amount and the photovoltaic power generation amount (z), and outputs this difference value as output data. The file 33 is stored (S203).

この差分(y−z)は、予測値送信手段24を介して、その地域を管轄する制御装置2側へ送信される(S205)。制御装置2では、受信した差分値に基づいて待機容量を確保する。   This difference (yz) is transmitted to the control device 2 side having jurisdiction over the area via the predicted value transmitting means 24 (S205). The control device 2 secures standby capacity based on the received difference value.

また、太陽光発電量予測手段23は、次の時刻の太陽光発電量の予測値と現在の太陽光発電量の実績値との差分を計算し(S204)、この計算結果を出力指令値として予測値送信手段24を介してその地域を管轄する制御装置2へ送信する(S205)。制御装置2では、この出力指令値に基づいて出力を変化させる。   Moreover, the photovoltaic power generation amount prediction means 23 calculates the difference between the predicted value of the photovoltaic power generation amount at the next time and the actual value of the current photovoltaic power generation amount (S204), and uses this calculation result as an output command value. It transmits to the control apparatus 2 which has jurisdiction over the area via the predicted value transmission means 24 (S205). The control device 2 changes the output based on this output command value.

本実施の形態によれば、地域ごとに太陽光発電量を上下限の幅を持って予測することができるので、制御装置を効率的に運転することができる。また、各家庭に備えられた太陽光発電装置による電力供給量や係数、設備容量を取得する必要がないので作業効率の向上、設備コストの削減を図ることができる。   According to the present embodiment, the amount of photovoltaic power generation can be predicted for each region with the upper and lower limits, so that the control device can be operated efficiently. Moreover, since it is not necessary to acquire the power supply amount, coefficient, and equipment capacity of the solar power generation device provided in each home, it is possible to improve work efficiency and reduce equipment costs.

本発明は、上述の実施の形態に限定されること無く、その要旨を逸脱しない範囲で種々変形して実現することができる。たとえば、太陽光発電量予測システムは、太陽光発電のみでなく風力発電等の他の電力調整値と合わせて需給計画や制御装置への目標出力を演算するようにしてもよい。また、太陽光発電量予測システムは太陽光発電量予測処理を実行する専用の装置として実現することもできるが、たとえば中央給電指令所の需給計画演算機能の一部として実現するなど、電力系統制御の他の機能と一体として実現することも可能である。   The present invention is not limited to the above-described embodiment, and can be implemented with various modifications without departing from the scope of the invention. For example, the photovoltaic power generation amount prediction system may calculate the supply and demand plan and the target output to the control device in combination with not only photovoltaic power generation but also other power adjustment values such as wind power generation. In addition, the photovoltaic power generation amount prediction system can be realized as a dedicated device for executing the photovoltaic power generation amount prediction processing, but for example, it is realized as a part of the supply and demand plan calculation function of the central power supply command station. It is also possible to realize it as a unit with other functions.

1 太陽光発電量予測システム
2 制御装置
3 家庭(電力需要家)
4 太陽光発電装置
5 発電所
6 電力量計
7 温度計
8 日射量計
9 気象業者のサーバ
11 通信装置
12 演算処理部
13 記憶部
21 データ入力手段
22 パラメータ推定手段
23 太陽光発電量予測手段
24 予測値送信手段
31 入力データファイル
32 パラメータファイル
33 出力データファイル
1 Photovoltaic power generation prediction system 2 Control device 3 Home (electric power consumer)
DESCRIPTION OF SYMBOLS 4 Solar power generation device 5 Power plant 6 Electricity meter 7 Thermometer 8 Solar radiation meter 9 Meteorological server 11 Communication device 12 Arithmetic processing unit 13 Storage unit 21 Data input means 22 Parameter estimation means 23 Solar power generation amount prediction means 24 Predicted value transmission means 31 Input data file 32 Parameter file 33 Output data file

Claims (6)

電力系統に分散設置された太陽光発電装置の太陽光発電量を予測する太陽光発電量予測システムであって、
地域ごとの気温、日射量の予報値と実績値および電力系統の電力量の各データを収集するデータ入力手段と、
前記収集した各データのうち、異なる複数の日について略同一時刻の気温と日射量実績値と電力量とを用いて、気温と日射量実績値を説明変数とし、電力量を目的変数としたときの回帰係数を一または複数の地域ごとに算出するパラメータ推定手段と、
前記パラメータ推定手段で算出した地域ごとの回帰係数、および日射量予報値から太陽光発電量の予測値を地域ごとに算出する太陽光発電量予測手段と、
を備えたことを特徴とする太陽光発電量予測システム。
A photovoltaic power generation amount prediction system for predicting the amount of photovoltaic power generation of photovoltaic power generation devices distributed in a power system,
A data input means for collecting each data of temperature, forecast value and actual value of solar radiation amount and electric power amount of electric power system for each region,
Of the collected data, when the temperature, the solar radiation actual value, and the electric energy at substantially the same time for a plurality of different days are used as the explanatory variables and the electric energy is the target variable Parameter estimation means for calculating the regression coefficient of each region or regions,
A solar power generation amount prediction means for calculating a prediction value of the solar power generation amount for each region from the regression coefficient for each region calculated by the parameter estimation means and the solar radiation amount forecast value;
A photovoltaic power generation amount prediction system characterized by comprising:
電力系統に分散設置された太陽光発電装置の太陽光発電量を予測する太陽光発電量予測システムであって、
地域ごとの日射量の予報値と実績値および電力系統の電力量の各データを収集するデータ入力手段と、
前記収集した各データを用いて、日射量実績値の単位時間の変化量を説明変数とし、電力量の単位時間の変化量を目的変数としたときの回帰係数を一または複数の地域ごとに算出するパラメータ推定手段と、
前記パラメータ推定手段で算出した地域ごとの回帰係数、および日射量予報値から太陽光発電量の予測値を地域ごとに算出する太陽光発電量予測手段と、
を備えたことを特徴とする太陽光発電量予測システム。
A photovoltaic power generation amount prediction system for predicting the amount of photovoltaic power generation of photovoltaic power generation devices distributed in a power system,
A data input means for collecting the forecast value and actual value of solar radiation amount for each region and the power amount data of the power system;
Using each collected data, the regression coefficient is calculated for each region or regions where the amount of change in unit time of solar radiation actual value is used as an explanatory variable and the amount of change in unit time of power is used as a target variable. Parameter estimation means for
A solar power generation amount prediction means for calculating a prediction value of the solar power generation amount for each region from the regression coefficient for each region calculated by the parameter estimation means and the solar radiation amount forecast value;
A photovoltaic power generation amount prediction system characterized by comprising:
前記太陽光発電量予測手段は、前記回帰係数の信頼区間を演算すると共に前記データ入力手段によって収集した日射量の予報値と実績値から日射量予報値誤差の信頼区間を演算し、当該回帰係数の信頼区間および当該日射量予報値誤差の信頼区間をもとに太陽光発電量の上限値と下限値を含む予測値を算出することを特徴とする請求項1または2に記載の太陽光発電量予測システム。   The photovoltaic power generation amount prediction means calculates a confidence interval of the regression coefficient and calculates a confidence interval of a solar radiation forecast value error from the forecast value and actual value of the solar radiation amount collected by the data input means, and the regression coefficient The solar power generation according to claim 1 or 2, wherein a predicted value including an upper limit value and a lower limit value of the solar power generation amount is calculated based on the confidence interval of the solar power generation amount and the confidence interval of the solar radiation amount forecast value error. Quantity prediction system. 地域ごとに前記太陽光発電量予測手段によって算出された太陽光発電量の予測値を、該地域を管轄する系統安定化用の制御装置へ送信する予測値送信手段を備えたことを特徴とする請求項1乃至3のいずれか一項に記載の太陽光発電量予測システム。   It is provided with a predicted value transmission means for transmitting the predicted value of the photovoltaic power generation amount calculated by the photovoltaic power generation amount prediction means for each region to a control device for system stabilization having jurisdiction over the region. The photovoltaic power generation amount prediction system according to any one of claims 1 to 3. 前記太陽光発電量予測手段によって演算された太陽光発電量の予測値に基づいて、前記制御装置の出力の目標値が算出されることを特徴とする請求項4に記載の太陽光発電量予測システム。   5. The photovoltaic power generation amount prediction according to claim 4, wherein a target value of the output of the control device is calculated based on a predicted value of the photovoltaic power generation amount calculated by the photovoltaic power generation amount prediction unit. system. 電力系統に分散設置された太陽光発電装置の発電量の実績値を計測せずに太陽光発電量を予測する太陽光発電量予測方法であって、
地域ごとの日射量の予報値と実績値および電力系統の送電端の電力量の各データを収集するデータ入力段階と、
前記収集した各データを用いて、日射量と電力量の単位時間の変化量の実績値から日射量を説明変数とし、電力量を目的変数としたときの回帰係数を地域ごとに算出するパラメータ推定段階と、
任意の確率の前記回帰係数の信頼区間および日射量予報値の信頼区間を算出すると共に、当該算出した信頼区間の値、前記パラメータ推定段階で算出した地域ごとの回帰係数、および日射量予報値から太陽光発電量の上限値と下限値とを含む予測値を地域ごとに算出する太陽光発電量予測段階と、
地域ごとに当該地域の予測値を当該地域を管轄する系統安定化用の制御装置へ送信する予測値送信段階と、を含み、
前記太陽光発電量予測段階で算出した前記予測値に基づいて、制御装置の起動または停止の判定、並びに、制御装置の出力を決定することを特徴とする太陽光発電量予測方法。
A photovoltaic power generation amount prediction method for predicting the amount of photovoltaic power generation without measuring the actual value of the amount of power generation of the photovoltaic power generation devices distributed in the power system,
A data input stage that collects forecast and actual values of solar radiation for each region and data on the power consumption at the transmission end of the power system,
Using each collected data, parameter estimation to calculate the regression coefficient for each region using the amount of solar radiation as the explanatory variable and the amount of power as the objective variable from the actual values of the amount of solar radiation and the amount of change in unit time Stages,
Calculate the confidence interval of the regression coefficient and the prediction value of the solar radiation amount with an arbitrary probability, and also calculate the confidence interval value, the regression coefficient for each region calculated in the parameter estimation step, and the solar radiation amount prediction value. A photovoltaic power generation prediction stage for calculating a predicted value including an upper limit value and a lower limit value of the solar power generation amount for each region;
A predicted value transmission stage for transmitting the predicted value of the region for each region to a control device for system stabilization that has jurisdiction over the region,
A method for predicting a photovoltaic power generation amount, comprising determining whether the control device is activated or stopped and determining an output of the control device based on the predicted value calculated in the photovoltaic power generation amount prediction step.
JP2010003377A 2010-01-08 2010-01-08 Photovoltaic power generation amount prediction system and solar power generation amount prediction method Expired - Fee Related JP5493886B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010003377A JP5493886B2 (en) 2010-01-08 2010-01-08 Photovoltaic power generation amount prediction system and solar power generation amount prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010003377A JP5493886B2 (en) 2010-01-08 2010-01-08 Photovoltaic power generation amount prediction system and solar power generation amount prediction method

Publications (2)

Publication Number Publication Date
JP2011142790A true JP2011142790A (en) 2011-07-21
JP5493886B2 JP5493886B2 (en) 2014-05-14

Family

ID=44458224

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010003377A Expired - Fee Related JP5493886B2 (en) 2010-01-08 2010-01-08 Photovoltaic power generation amount prediction system and solar power generation amount prediction method

Country Status (1)

Country Link
JP (1) JP5493886B2 (en)

Cited By (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012044740A (en) * 2010-08-16 2012-03-01 Tokyo Electric Power Co Inc:The Estimation method for amount of solar power generation and distribution system control system
JP2013029923A (en) * 2011-07-27 2013-02-07 Chugoku Electric Power Co Inc:The Photovoltaic power generation amount grasping system, load estimation device employing the same and load adjustment device
JP2013046494A (en) * 2011-08-24 2013-03-04 Denso Corp Charge control system for electric vehicle
JP2013066278A (en) * 2011-09-16 2013-04-11 Tokyo Gas Co Ltd Power system supporting system
JP2013192401A (en) * 2012-03-14 2013-09-26 Toshiba Corp Power demand control apparatus
JP2013208050A (en) * 2012-03-29 2013-10-07 Mitsubishi Electric Corp Method for predicting output of photovoltaic power generation device
WO2013157135A1 (en) * 2012-04-20 2013-10-24 株式会社日立製作所 Power monitoring device and method
JP2013232147A (en) * 2012-05-01 2013-11-14 Tokyo Electric Power Co Inc:The Generated power estimation device, generated power estimation method, and generated power estimation program
JP2013253805A (en) * 2012-06-05 2013-12-19 Fuji Electric Co Ltd Information processing device, and control method and program for the same
WO2014002274A1 (en) * 2012-06-29 2014-01-03 株式会社日立製作所 Fossil-fuel power generating equipment, renewable energy power generating plant, and control method for same
WO2014203431A1 (en) * 2013-06-21 2014-12-24 株式会社 東芝 Prediction system, prediction device, and prediction method
JP2016038348A (en) * 2014-08-11 2016-03-22 富士通株式会社 Plan formulation method, plan formulation system and plan formulation program
JP2016073154A (en) * 2014-10-01 2016-05-09 中国電力株式会社 Power prediction device, power prediction method, and program
JP2016073155A (en) * 2014-10-01 2016-05-09 中国電力株式会社 Demand prediction device, demand prediction method, and program
JP2016073156A (en) * 2014-10-01 2016-05-09 中国電力株式会社 Power generation amount prediction device, power generation amount prediction method, and program
JP2016152644A (en) * 2015-02-16 2016-08-22 北陸電力株式会社 Estimation method for photovoltaic power generation output, and device
JP2016163511A (en) * 2015-03-05 2016-09-05 中国電力株式会社 Power demand amount prediction system, power demand amount prediction method, and program
JP2016197948A (en) * 2015-04-03 2016-11-24 富士電機株式会社 Power estimation device and power estimation method
JPWO2014129045A1 (en) * 2013-02-19 2017-02-02 日本電気株式会社 Power flow control system and power flow control method
JP2017060224A (en) * 2015-09-14 2017-03-23 株式会社東芝 Aggregation management device and aggregation management method
JP2017084360A (en) * 2015-10-23 2017-05-18 富士通株式会社 Method for operating photovoltaic power generation system
JP2017135964A (en) * 2016-01-29 2017-08-03 富士通株式会社 Operation of photovoltaic power generation system
JP2018007370A (en) * 2016-06-29 2018-01-11 関西電力株式会社 Photovoltaic power generation output estimation device, and photovoltaic power generation output estimation method
JP2018007461A (en) * 2016-07-05 2018-01-11 日本電気株式会社 Estimation system, estimation method, and estimation program
JP2018011494A (en) * 2016-07-15 2018-01-18 関西電力株式会社 Device and method for estimating photovoltaic power generation output
JP2018018329A (en) * 2016-07-28 2018-02-01 富士通株式会社 Power amount prediction program, power amount prediction device and power amount prediction method
JP2018067155A (en) * 2016-10-19 2018-04-26 株式会社東芝 Power generation plan formulation device, power generation plan formulation method, and power generation plan formulation program
JP2018207668A (en) * 2017-06-05 2018-12-27 中国電力株式会社 Power generation amount prediction device and power generation amount prediction method
JP2020028195A (en) * 2018-08-13 2020-02-20 一般財団法人電力中央研究所 Photovoltaic power generation system distribution estimation device, photovoltaic power generation system distribution estimation method, photovoltaic power generation system distribution estimation program, and photovoltaic power generation output prediction device
JP2020072552A (en) * 2018-10-31 2020-05-07 株式会社日立パワーソリューションズ Power generation system, control device, and setting method of power generation target value
CN112761896A (en) * 2020-09-24 2021-05-07 国网内蒙古东部电力有限公司 Calculation method and device for improving power generation amount prediction accuracy of wind power station and computer equipment
JP7043661B1 (en) 2021-05-31 2022-03-29 株式会社東芝 Information processing equipment, information processing methods, programs and power generation systems
WO2022168357A1 (en) * 2021-02-08 2022-08-11 株式会社日立製作所 System for managing power generation amount and method for managing power generation amount
CN115906476A (en) * 2022-11-18 2023-04-04 国网湖北省电力有限公司经济技术研究院 Mountain land photovoltaic power generation capacity calculation method
CN117040137A (en) * 2023-10-09 2023-11-10 国网山东省电力公司聊城供电公司 Ring main unit temperature rise early warning method, system, terminal and medium based on multi-source data

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000333367A (en) * 1999-05-19 2000-11-30 Toshiba Corp Thermoelectric consumption predicting system
JP2004032989A (en) * 2002-04-30 2004-01-29 Atsushi Iga Storage battery capacity of system combined with storage battery in solar cell, merit calculation method, and storage battery charge discharge employment method
JP2004082408A (en) * 2002-08-23 2004-03-18 Ricoh Co Ltd Image forming apparatus
JP2004289918A (en) * 2003-03-20 2004-10-14 Fujitsu Ltd Power supply method
JP2006033908A (en) * 2004-07-12 2006-02-02 Nippon Telegr & Teleph Corp <Ntt> Method, device and program for estimating amount of power generation of solar light generating system
JP2009077498A (en) * 2007-09-19 2009-04-09 Osaka Gas Co Ltd Load reduction plan decision system and demand control system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000333367A (en) * 1999-05-19 2000-11-30 Toshiba Corp Thermoelectric consumption predicting system
JP2004032989A (en) * 2002-04-30 2004-01-29 Atsushi Iga Storage battery capacity of system combined with storage battery in solar cell, merit calculation method, and storage battery charge discharge employment method
JP2004082408A (en) * 2002-08-23 2004-03-18 Ricoh Co Ltd Image forming apparatus
JP2004289918A (en) * 2003-03-20 2004-10-14 Fujitsu Ltd Power supply method
JP2006033908A (en) * 2004-07-12 2006-02-02 Nippon Telegr & Teleph Corp <Ntt> Method, device and program for estimating amount of power generation of solar light generating system
JP2009077498A (en) * 2007-09-19 2009-04-09 Osaka Gas Co Ltd Load reduction plan decision system and demand control system

Cited By (47)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012044740A (en) * 2010-08-16 2012-03-01 Tokyo Electric Power Co Inc:The Estimation method for amount of solar power generation and distribution system control system
JP2013029923A (en) * 2011-07-27 2013-02-07 Chugoku Electric Power Co Inc:The Photovoltaic power generation amount grasping system, load estimation device employing the same and load adjustment device
JP2013046494A (en) * 2011-08-24 2013-03-04 Denso Corp Charge control system for electric vehicle
JP2013066278A (en) * 2011-09-16 2013-04-11 Tokyo Gas Co Ltd Power system supporting system
JP2013192401A (en) * 2012-03-14 2013-09-26 Toshiba Corp Power demand control apparatus
JP2013208050A (en) * 2012-03-29 2013-10-07 Mitsubishi Electric Corp Method for predicting output of photovoltaic power generation device
WO2013157135A1 (en) * 2012-04-20 2013-10-24 株式会社日立製作所 Power monitoring device and method
JPWO2013157135A1 (en) * 2012-04-20 2015-12-21 株式会社日立製作所 Power monitoring apparatus and power monitoring method
US20150073737A1 (en) * 2012-04-20 2015-03-12 Tatsuki Inuzuka Power monitoring apparatus and power monitoring method
JP2013232147A (en) * 2012-05-01 2013-11-14 Tokyo Electric Power Co Inc:The Generated power estimation device, generated power estimation method, and generated power estimation program
JP2013253805A (en) * 2012-06-05 2013-12-19 Fuji Electric Co Ltd Information processing device, and control method and program for the same
WO2014002274A1 (en) * 2012-06-29 2014-01-03 株式会社日立製作所 Fossil-fuel power generating equipment, renewable energy power generating plant, and control method for same
JPWO2014129045A1 (en) * 2013-02-19 2017-02-02 日本電気株式会社 Power flow control system and power flow control method
US10069302B2 (en) 2013-02-19 2018-09-04 Nec Corporation Power flow control system and power flow control method
JP2015005641A (en) * 2013-06-21 2015-01-08 株式会社東芝 Predicting system, predicting apparatus and predicting method
CN105308853A (en) * 2013-06-21 2016-02-03 株式会社东芝 Prediction system, prediction device, and prediction method
US20140379311A1 (en) * 2013-06-21 2014-12-25 Kabushiki Kaisha Toshiba Estimation system, estimation apparatus, and estimation method
WO2014203431A1 (en) * 2013-06-21 2014-12-24 株式会社 東芝 Prediction system, prediction device, and prediction method
JP2016038348A (en) * 2014-08-11 2016-03-22 富士通株式会社 Plan formulation method, plan formulation system and plan formulation program
JP2016073154A (en) * 2014-10-01 2016-05-09 中国電力株式会社 Power prediction device, power prediction method, and program
JP2016073155A (en) * 2014-10-01 2016-05-09 中国電力株式会社 Demand prediction device, demand prediction method, and program
JP2016073156A (en) * 2014-10-01 2016-05-09 中国電力株式会社 Power generation amount prediction device, power generation amount prediction method, and program
JP2016152644A (en) * 2015-02-16 2016-08-22 北陸電力株式会社 Estimation method for photovoltaic power generation output, and device
JP2016163511A (en) * 2015-03-05 2016-09-05 中国電力株式会社 Power demand amount prediction system, power demand amount prediction method, and program
JP2016197948A (en) * 2015-04-03 2016-11-24 富士電機株式会社 Power estimation device and power estimation method
JP2017060224A (en) * 2015-09-14 2017-03-23 株式会社東芝 Aggregation management device and aggregation management method
JP2017084360A (en) * 2015-10-23 2017-05-18 富士通株式会社 Method for operating photovoltaic power generation system
JP2017135964A (en) * 2016-01-29 2017-08-03 富士通株式会社 Operation of photovoltaic power generation system
JP2018007370A (en) * 2016-06-29 2018-01-11 関西電力株式会社 Photovoltaic power generation output estimation device, and photovoltaic power generation output estimation method
JP2018007461A (en) * 2016-07-05 2018-01-11 日本電気株式会社 Estimation system, estimation method, and estimation program
JP2018011494A (en) * 2016-07-15 2018-01-18 関西電力株式会社 Device and method for estimating photovoltaic power generation output
JP2018018329A (en) * 2016-07-28 2018-02-01 富士通株式会社 Power amount prediction program, power amount prediction device and power amount prediction method
JP2018067155A (en) * 2016-10-19 2018-04-26 株式会社東芝 Power generation plan formulation device, power generation plan formulation method, and power generation plan formulation program
JP2018207668A (en) * 2017-06-05 2018-12-27 中国電力株式会社 Power generation amount prediction device and power generation amount prediction method
JP7257117B2 (en) 2018-08-13 2023-04-13 一般財団法人電力中央研究所 Photovoltaic power generation system distribution estimation device, photovoltaic power generation system distribution estimation method, and photovoltaic power generation output prediction device
JP2020028195A (en) * 2018-08-13 2020-02-20 一般財団法人電力中央研究所 Photovoltaic power generation system distribution estimation device, photovoltaic power generation system distribution estimation method, photovoltaic power generation system distribution estimation program, and photovoltaic power generation output prediction device
JP2020072552A (en) * 2018-10-31 2020-05-07 株式会社日立パワーソリューションズ Power generation system, control device, and setting method of power generation target value
CN112761896A (en) * 2020-09-24 2021-05-07 国网内蒙古东部电力有限公司 Calculation method and device for improving power generation amount prediction accuracy of wind power station and computer equipment
CN112761896B (en) * 2020-09-24 2024-05-14 国网内蒙古东部电力有限公司 Calculation method, device and computer equipment for improving prediction accuracy of power generation amount of wind power station
WO2022168357A1 (en) * 2021-02-08 2022-08-11 株式会社日立製作所 System for managing power generation amount and method for managing power generation amount
JP7410897B2 (en) 2021-02-08 2024-01-10 株式会社日立製作所 Power generation management system and power generation management method
JP7043661B1 (en) 2021-05-31 2022-03-29 株式会社東芝 Information processing equipment, information processing methods, programs and power generation systems
JP2022183783A (en) * 2021-05-31 2022-12-13 株式会社東芝 Information processing device, information processing method, program, and power generation system
CN115906476A (en) * 2022-11-18 2023-04-04 国网湖北省电力有限公司经济技术研究院 Mountain land photovoltaic power generation capacity calculation method
CN115906476B (en) * 2022-11-18 2023-09-01 国网湖北省电力有限公司经济技术研究院 Mountain photovoltaic power generation capacity calculation method
CN117040137A (en) * 2023-10-09 2023-11-10 国网山东省电力公司聊城供电公司 Ring main unit temperature rise early warning method, system, terminal and medium based on multi-source data
CN117040137B (en) * 2023-10-09 2024-05-07 国网山东省电力公司聊城供电公司 Ring main unit temperature rise early warning method, system, terminal and medium based on multi-source data

Also Published As

Publication number Publication date
JP5493886B2 (en) 2014-05-14

Similar Documents

Publication Publication Date Title
JP5493886B2 (en) Photovoltaic power generation amount prediction system and solar power generation amount prediction method
EP2472694B1 (en) Methods and systems for controlling a power conversion system
JP5806132B2 (en) Power generation amount prediction device, power generation amount prediction correction method, and natural energy power generation system
CN108808743B (en) Communication network-based multi-energy micro-grid energy prediction and optimal scheduling method
US10298056B2 (en) Power control system, power control method, and recording medium
JP6543145B2 (en) Peak power prediction device, power management system and peak power prediction method
CA3032428A1 (en) Electrical energy storage system with battery resistance estimation
US9638545B2 (en) Power management apparatus, power management system and power management method
JP2014098601A (en) Insolation amount calculation method and power supply determining method
JP2010193605A (en) Load estimating method of power distribution section and power distribution system control method
US11018523B2 (en) Utility grid, intermittent energy management system
JP5962770B2 (en) Power system control system and power system control method
KR102396712B1 (en) Energy management system and method for minimizing the power purchase cost of microgrid using the same
JP5565193B2 (en) Load prediction method for distribution section and distribution system control system
JP2018160990A (en) Load management apparatus and method
JP2012043857A (en) Bright and clear time determination device and photovoltaic power generation amount prediction system
US9851734B2 (en) Alert presentation apparatus and alert presentation method
Schulte et al. Forecast-based charging strategy to prolong the lifetime of lithium-ion batteries in standalone PV battery systems in Sub-Saharan Africa
JP6833303B1 (en) Power generation forecaster
Schweighofer et al. Optimized selection of component models for photovoltaic and energy storage system simulations
Hokmabad et al. Optimizing size and economic feasibility assessment of photovoltaic and energy storage setup in residential applications
JP5893507B2 (en) Alarm presenting device and alarm presenting method
JP5882135B2 (en) Demand information presentation apparatus and demand information presentation method
Folesa et al. Implementation of PV Self-Consumption with Ramp-Rate Control Algorithm using a Day-to-Day Forecast battery charging, with a Vanadium Redox Flow Battery.
Xingyu et al. Uncertainty Analysis for Day Ahead Power Reserve Quantification in an Urban Microgrid Including PV Generators 2

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20110422

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20121114

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20131121

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131126

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131226

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140217

R150 Certificate of patent or registration of utility model

Ref document number: 5493886

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees