JP2022084398A - Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program - Google Patents

Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program Download PDF

Info

Publication number
JP2022084398A
JP2022084398A JP2020196254A JP2020196254A JP2022084398A JP 2022084398 A JP2022084398 A JP 2022084398A JP 2020196254 A JP2020196254 A JP 2020196254A JP 2020196254 A JP2020196254 A JP 2020196254A JP 2022084398 A JP2022084398 A JP 2022084398A
Authority
JP
Japan
Prior art keywords
power generation
photovoltaic power
generation output
coefficient
estimation
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
JP2020196254A
Other languages
Japanese (ja)
Other versions
JP7450522B2 (en
JP2022084398A5 (en
Inventor
弘 奥村
Hiroshi Okumura
一浩 安並
Kazuhiro Yasunami
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2020196254A priority Critical patent/JP7450522B2/en
Publication of JP2022084398A publication Critical patent/JP2022084398A/en
Publication of JP2022084398A5 publication Critical patent/JP2022084398A5/ja
Application granted granted Critical
Publication of JP7450522B2 publication Critical patent/JP7450522B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

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

Abstract

To provide a solar power generation output estimation device capable of estimating the solar power generation output of a solar power generation facility without newly installing many solar radiation meters or measuring instruments having sufficiently high time resolution.SOLUTION: In a solar power generation output estimation system 10, a solar power generation output estimation device 100 includes: an estimation unit 150 configured to estimate a second solar power generation output that is a power generation output of a second solar power generation facility installed within a predetermined distance from a first solar power generation facility, using a first solar power generation output that is a power generation output of the first solar power generation facility and a residual demand that is apparent power consumption of a consumer who has installed the second solar power generation facility; and an estimation accuracy verification unit 140 configured to determine whether or not estimation accuracy of the second solar power generation output satisfies a specified accuracy using the power generation output of a third solar power generation facility installed within a predetermined distance from the first solar power generation facility and a residual demand that is apparent power consumption of a consumer who has installed the third solar power generation facility.SELECTED DRAWING: Figure 2

Description

本開示は、太陽光発電設備の発電出力を推定する太陽光発電出力推定装置、太陽光発電出力推定方法および太陽光発電出力推定プログラムに関する。 The present disclosure relates to a photovoltaic power generation output estimation device for estimating the power output of a photovoltaic power generation facility, a photovoltaic power generation output estimation method, and a photovoltaic power generation output estimation program.

近年、再生可能エネルギーの利用を拡大することの重要性が高まっており、太陽光発電設備などの分散型電源を設置し、送電系統または配電系統(以下、これらを電力系統という)に電力を供給する需要家が増えてきている。一方、電力系統を運用する電力会社などでは、一部の太陽光発電設備の発電出力は把握しているものの、多くの太陽光発電設備の発電出力は把握できていない。このため、各需要家が設置した太陽光発電設備の発電出力を推定する装置が提案されている(例えば、特許文献1参照)。なお、上記の発電出力を把握できる太陽光発電設備は、全量買取制度の契約(以下、全量買取契約)の対象となっている太陽光発電設備であり、発電出力を把握できない太陽光発電設備は、余剰電力買取制度の契約(以下、余剰買取契約)の対象となっている太陽光発電設備である。これは、全量買取契約の対象となっている太陽光発電設備を保有する需要家には、太陽光発電設備の発電量を計測するスマートメータと、需要家内の負荷の消費電力量を計測するスマートメータが個別に設置されているのに対し、余剰買取契約の対象となっている太陽光発電設備を保有する需要家には、太陽発電設備の発電量と負荷の消費電力量との合算値を計測するスマートメータしか設置されていないためである。 In recent years, the importance of expanding the use of renewable energy has increased, and distributed power sources such as photovoltaic power generation facilities have been installed to supply power to power transmission systems or distribution systems (hereinafter referred to as power systems). The number of consumers is increasing. On the other hand, although electric power companies that operate electric power systems know the power generation output of some photovoltaic power generation facilities, they cannot grasp the power generation output of many photovoltaic power generation facilities. Therefore, a device for estimating the power generation output of the photovoltaic power generation equipment installed by each consumer has been proposed (see, for example, Patent Document 1). The above-mentioned photovoltaic power generation equipment that can grasp the power generation output is the photovoltaic power generation equipment that is subject to the contract of the total purchase system (hereinafter referred to as the total purchase contract), and the photovoltaic power generation equipment that cannot grasp the power generation output is , It is a photovoltaic power generation facility that is subject to the surplus power purchase system contract (hereinafter referred to as the surplus purchase contract). This is a smart meter that measures the amount of power generated by the photovoltaic power generation equipment and a smart that measures the power consumption of the load inside the customer for consumers who own the photovoltaic power generation equipment that is the subject of the total purchase contract. While the meters are installed individually, for consumers who own the photovoltaic power generation equipment that is the subject of the surplus purchase contract, the total value of the power generation amount of the solar power generation equipment and the power consumption of the load is calculated. This is because only a smart meter for measuring is installed.

特許文献1には、電力系統に接続された負荷に電力を供給する太陽光発電設備の所定時点での発電出力である太陽光発電出力を推定する太陽光発電出力推定装置であって、所定時点以前の期間であって、太陽の南中高度が所定時点での当該南中高度から所定範囲内の期間における、所定地点での日射強度と電力系統の有効電力とを用いて、太陽光発電出力を推定する発電出力推定部を備えることを特徴とする、太陽光発電出力推定装置が開示されている。この太陽光発電出力推定装置は、変電所周辺に設置した日射計などで測定した日射強度と、変電所などに設置した計測器で測定した有効電力とを用いて、電力系統に連系した太陽光発電設備の発電出力を推定する。 Patent Document 1 is a photovoltaic power generation output estimation device that estimates a photovoltaic power generation output, which is a photovoltaic power generation output at a predetermined time point of a photovoltaic power generation facility that supplies power to a load connected to a power system, and is a photovoltaic power generation output estimation device at a predetermined time point. In the previous period, when the south-middle altitude of the sun is within a predetermined range from the south-middle altitude at a predetermined point, the photovoltaic power output is used by using the solar radiation intensity at a predetermined point and the active power of the power system. A photovoltaic power generation output estimation device is disclosed, which comprises a power generation output estimation unit for estimating. This photovoltaic power generation output estimation device uses the solar radiation intensity measured by a solar radiation meter installed around the substation and the active power measured by a measuring instrument installed in the substation, etc., to connect the sun to the power system. Estimate the power output of the photovoltaic power generation facility.

特開2016-139270号公報Japanese Unexamined Patent Publication No. 2016-139270

電力会社などが多くの太陽光発電設備の発電出力を把握できていない状況下において、太陽光発電設備が大量に導入された場合、電力系統の運用には様々な問題が生じる。電力会社などでは、自身が所管する変電所などで、太陽光発電設備の発電出力を加味した見かけ上の電力系統の負荷(以下、「見かけ上の負荷」という)の消費電力を計測器などで計測できているものの、電力系統に連系されている全ての太陽光発電設備の発電出力合計値が未知であるため、実際の負荷(以下、「実負荷」という)の消費電力を正確に把握することができない。このため、太陽光発電設備が大量に導入された場合、重回帰分析などを利用して予測を行っている実負荷の消費電力の予測誤差が大きくなるという需給制御上のリスクが増大したり、系統事故後の復旧操作に支障が生じるという系統制御上のリスクが増大したりする。 If a large number of photovoltaic power generation facilities are introduced in a situation where electric power companies and the like cannot grasp the power generation output of many photovoltaic power generation facilities, various problems will occur in the operation of the power system. At electric power companies, etc., at substations under their jurisdiction, the power consumption of the apparent power system load (hereinafter referred to as "apparent load") that takes into account the power output of the photovoltaic power generation equipment is measured with a measuring instrument or the like. Although it can be measured, the total power output of all photovoltaic power generation facilities connected to the power system is unknown, so the power consumption of the actual load (hereinafter referred to as "actual load") can be accurately grasped. Can not do it. For this reason, when a large number of photovoltaic power generation facilities are introduced, there is an increased risk in supply and demand control that the prediction error of the actual load power consumption, which is predicted using multiple regression analysis, becomes large. There is an increased risk in system control that the recovery operation after a system accident will be hindered.

このため、太陽光発電設備の発電出力を正確に推定する必要がある。しかし、特許文献1の太陽光発電出力推定装置によれば、日射強度を測定するために多くの日射計を設置する必要性と、十分高い時間分解能を有する、すなわち十分高いサンプリング周期で計測と記録が可能な計測器を設置する必要性が生じる。また、変電所などの単位で推定を行うため、対象となる電力系統に連系された太陽光発電設備の導入量が実負荷と比較して小さくなるほど、推定精度が低下するという問題もある。それゆえ、電力系統に連系された太陽光発電設備の発電出力を高精度に推定するためには、技術的にもコスト的にも問題がある。 Therefore, it is necessary to accurately estimate the power generation output of the photovoltaic power generation equipment. However, according to the photovoltaic power generation output estimation device of Patent Document 1, it is necessary to install many pyranometers in order to measure the solar radiation intensity, and the measurement and recording have a sufficiently high time resolution, that is, a sufficiently high sampling cycle. It becomes necessary to install a measuring instrument that can be used. Further, since the estimation is performed in units of substations and the like, there is also a problem that the estimation accuracy decreases as the installed amount of the photovoltaic power generation equipment connected to the target power system becomes smaller than the actual load. Therefore, in order to estimate the power generation output of the photovoltaic power generation equipment connected to the power system with high accuracy, there are technical and cost problems.

本開示は、上記に鑑みてなされたものであって、多くの日射計および十分高い時間分解能をもった計測器を新たに設置することなく、太陽光発電設備の太陽光発電出力を推定することができる太陽光発電出力推定装置を得ることを目的とする。 This disclosure is made in view of the above, and estimates the photovoltaic output of a photovoltaic power generation facility without newly installing many pyranometers and measuring instruments with sufficiently high time resolution. The purpose is to obtain a photovoltaic power generation output estimation device that can be used.

上述した課題を解決し、目的を達成するために、本開示にかかる太陽光発電出力推定装置は、第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、第1太陽光発電設備の発電出力である第1太陽光発電出力と第2太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて推定する推定部、を備える。太陽光発電出力推定装置は、さらに、第1太陽光発電設備から予め定められた距離内に設置された第3太陽光発電設備の発電出力と第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定する推定精度検証部と、を備える。 In order to solve the above-mentioned problems and achieve the object, the photovoltaic power generation output estimation device according to the present disclosure is a second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility. The second photovoltaic power generation output, which is the power generation output, is the apparent power consumption of the consumer who has installed the first photovoltaic power generation output, which is the power generation output of the first photovoltaic power generation facility, and the second photovoltaic power generation facility. It is equipped with an estimation unit that estimates using the residual demand. The photovoltaic power generation output estimation device is further used for the power output of the third photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility and for the consumer in which the third photovoltaic power generation facility is installed. It is provided with an estimation accuracy verification unit for determining whether or not the estimation accuracy of the second photovoltaic power generation output satisfies the predetermined accuracy by using the residual demand which is the apparent power consumption.

本開示によれば、多くの日射計および十分高い時間分解能をもった計測器を新たに設置することなく、太陽光発電設備の太陽光発電出力を推定することができるという効果を奏する。 According to the present disclosure, it is possible to estimate the photovoltaic power generation output of a photovoltaic power generation facility without newly installing many pyranometers and measuring instruments having sufficiently high time resolution.

実施の形態1にかかる太陽光発電出力推定システムの構成例を示す図The figure which shows the configuration example of the photovoltaic power generation output estimation system which concerns on Embodiment 1. 実施の形態1の太陽光発電出力推定装置の機能構成例を示すブロック図Block diagram showing a functional configuration example of the photovoltaic power generation output estimation device of the first embodiment 実施の形態1の太陽光発電出力推定装置を実現するコンピュータシステムの構成例を示す図The figure which shows the configuration example of the computer system which realizes the photovoltaic power generation output estimation apparatus of Embodiment 1. 実施の形態1の太陽光発電出力推定装置による係数αの推定処理手順の一例を示すフローチャートA flowchart showing an example of an estimation processing procedure of the coefficient α by the photovoltaic power generation output estimation device of the first embodiment. 実施の形態1の場所情報データにより示される各需要家の太陽光発電設備の地理的位置を模式的に示す図The figure which shows typically the geographical position of the photovoltaic power generation facility of each consumer shown by the place information data of Embodiment 1. 実施の形態1の記憶部に格納される算出データ内のデータの一例を示す図The figure which shows an example of the data in the calculated data stored in the storage part of Embodiment 1. 実施の形態1の記憶部に格納される、太陽光発電設備の組み合わせごとの係数αの一例を示す図The figure which shows an example of the coefficient α for each combination of the photovoltaic power generation facilities stored in the storage part of Embodiment 1. 実施の形態1の第1共分散算出部による第1共分散の算出処理手順の一例を示すフローチャートA flowchart showing an example of a procedure for calculating the first covariance by the first covariance calculation unit of the first embodiment. 実施の形態1の第2共分散算出部による第2共分散の算出処理手順の一例を示すフローチャートA flowchart showing an example of a procedure for calculating the second covariance by the second covariance calculation unit of the first embodiment. 実施の形態1の係数αの算出処理手順の一例を示すフローチャートA flowchart showing an example of the calculation processing procedure of the coefficient α of the first embodiment. 実施の形態1の太陽光発電出力の推定処理手順の一例を示すフローチャートA flowchart showing an example of the procedure for estimating the photovoltaic power generation output of the first embodiment. 複数の需要家の太陽光発電設備を推定対象とする場合の第1太陽光発電設備および第2太陽光発電設備の一例を示す図The figure which shows an example of the 1st solar power generation facility and the 2nd solar power generation facility when the solar power generation facility of a plurality of consumers is estimated. 複数の需要家の太陽光発電設備を推定対象とする場合の係数αの一例を示す図A diagram showing an example of the coefficient α when the estimation target is the photovoltaic power generation equipment of multiple consumers. 実施の形態1の発電出力の推定方法の検証結果の一例を示す図The figure which shows an example of the verification result of the power generation output estimation method of Embodiment 1. 推定精度の検証を行わない場合の実施の形態1の発電出力の推定方法の検証結果の一例を示す図The figure which shows an example of the verification result of the estimation method of the power generation output of Embodiment 1 when the estimation accuracy is not verified. 図14に検証結果が示された日の前日の検証結果の一例を示す図FIG. 14 is a diagram showing an example of the verification result on the day before the day when the verification result is shown. 図15に検証結果が示された日の前日の検証結果の一例を示す図FIG. 15 is a diagram showing an example of the verification result on the day before the day when the verification result is shown. |λ-1|と第2太陽光発電出力の推定誤差との散布図の一例を示す図A diagram showing an example of a scatter diagram of | λ-1 | and the estimation error of the second photovoltaic power generation output. 発電出力の推定の順序の一例を示す模式図Schematic diagram showing an example of the order of estimation of power generation output 実施の形態2の係数αの算出処理手順の一例を示すフローチャートA flowchart showing an example of the calculation processing procedure of the coefficient α of the second embodiment. 実施の形態3の太陽光発電出力推定装置による係数αの推定処理手順の一例を示すフローチャートA flowchart showing an example of an estimation processing procedure of the coefficient α by the photovoltaic power generation output estimation device of the third embodiment.

以下に、実施の形態にかかる太陽光発電出力推定装置、太陽光発電出力推定方法および太陽光発電出力推定プログラムを図面に基づいて詳細に説明する。 Hereinafter, the photovoltaic power generation output estimation device, the photovoltaic power generation output estimation method, and the photovoltaic power generation output estimation program according to the embodiment will be described in detail based on the drawings.

実施の形態1.
図1は、実施の形態1にかかる太陽光発電出力推定システムの構成例を示す図である。本実施の形態の太陽光発電出力推定システム10は、太陽光発電出力推定装置100と、複数の需要家の下に設置された各種の設備とが、通信ネットワーク50を介して接続された構成である。図1において実線は電力の流れを表し、破線は情報の流れを表している。
Embodiment 1.
FIG. 1 is a diagram showing a configuration example of the photovoltaic power generation output estimation system according to the first embodiment. The photovoltaic power generation output estimation system 10 of the present embodiment has a configuration in which the photovoltaic power generation output estimation device 100 and various facilities installed under a plurality of consumers are connected via a communication network 50. be. In FIG. 1, the solid line represents the flow of electric power, and the broken line represents the flow of information.

図1では、複数の需要家として需要家200-1、需要家200-2、需要家200-3、需要家200-4、・・・、需要家200-m、・・・、需要家200-nを図示している。需要家200-1~200-nのそれぞれには、負荷201および太陽光発電設備(図では、PV(PhotoVoltaic)と略す)202が接続されている。なお、図1に示した例では、nは6以上の整数であり、mは5以上n未満の整数であるが、n,mの値は、図1に示した例に限定されない。また、各需要家200-1~200-nの負荷201は、電力を消費する1つ以上の設備を示しており、一般には、負荷201を構成する具体的な設備は全需要家200-1~200-nで同一ではない。太陽光発電設備202についても、同様に、一般には、全需要家200-1~200-nで同一ではない。 In FIG. 1, as a plurality of consumers, consumer 200-1, consumer 200-2, consumer 200-3, consumer 200-4, ..., consumer 200-m, ..., consumer 200 -N is illustrated. A load 201 and a photovoltaic power generation facility (abbreviated as PV (PhotoVoltaic) in the figure) 202 are connected to each of the consumers 200-1 to 200-n. In the example shown in FIG. 1, n is an integer of 6 or more and m is an integer of 5 or more and less than n, but the values of n and m are not limited to the example shown in FIG. Further, the load 201 of each consumer 200-1 to 200-n indicates one or more facilities that consume electric power, and in general, the specific equipment constituting the load 201 is all consumers 200-1. ~ 200-n is not the same. Similarly, the photovoltaic power generation facility 202 is not the same for all consumers 200-1 to 200-n.

需要家200-1に設けられている太陽光発電設備202は、全量買取契約対象の太陽光発電設備である。需要家200-2~200-nに設けられている太陽光発電設備202は、余剰買取契約対象の太陽光発電設備である。 The photovoltaic power generation facility 202 provided in the consumer 200-1 is a photovoltaic power generation facility subject to a purchase contract in its entirety. The photovoltaic power generation equipment 202 installed in the consumers 200-2 to 200-n is a photovoltaic power generation facility subject to a surplus purchase contract.

太陽光発電設備202は、例えば、需要家200-1~200-nが有する建築物の屋根などに設置されている太陽光発電設備であり、小規模な太陽光発電設備のみならず、いわゆるメガソーラーなどの大規模太陽光発電所を含む。太陽光発電設備202は、電力系統20に連系されており、発電した電力を電力系統20に供給する。また、電力系統20には、図1に示すように、各需要家200-1~200-nの負荷201が接続されている。負荷201は電力系統20から電力供給を受け、電力を消費する。 The photovoltaic power generation facility 202 is, for example, a photovoltaic power generation facility installed on the roof of a building owned by consumers 200-1 to 200-n, and is not only a small-scale solar power generation facility but also a so-called mega. Includes large-scale solar power plants such as solar. The photovoltaic power generation facility 202 is connected to the electric power system 20, and supplies the generated electric power to the electric power system 20. Further, as shown in FIG. 1, a load 201 of each consumer 200-1 to 200-n is connected to the power system 20. The load 201 receives power from the power system 20 and consumes power.

負荷201および太陽光発電設備202が設けられている需要家200-1~200-nは、電力系統20からの電力供給を受けるとともに太陽光発電設備202によって発電された電力を電力系統20に供給することになる。このため、電力系統20から見た各需要家200-1~200-nの見かけ上の消費電力は、負荷201の実際の消費電力と太陽光発電設備202による発電出力との両方に依存することになる。なお、負荷201による消費電力を太陽光発電設備202による発電出力が上回る場合には、電力系統20からみると当該需要家の発電出力が電力系統20に供給されることになるため、見かけ上の消費電力は負の値になる。以下では、電力系統20から見た、各需要家200-1~200-nの見かけ上の消費電力を、残余需要と呼ぶ。 Consumers 200-1 to 200-n provided with the load 201 and the photovoltaic power generation facility 202 receive power from the power system 20 and supply the power generated by the photovoltaic power generation facility 202 to the power system 20. Will be done. Therefore, the apparent power consumption of each consumer 200-1 to 200-n as seen from the power system 20 depends on both the actual power consumption of the load 201 and the power generation output of the photovoltaic power generation facility 202. become. When the power generation output of the photovoltaic power generation facility 202 exceeds the power consumption of the load 201, the power generation output of the consumer is supplied to the power system 20 from the viewpoint of the power system 20. The power consumption becomes a negative value. In the following, the apparent power consumption of each consumer 200-1 to 200-n as seen from the power system 20 will be referred to as residual demand.

余剰買取契約を締結している需要家200-2~200-nには、スマートメータ(図では、SM(Smart Meter)と略す)205が設置されている。iを2からnまでの任意の整数とすると、需要家200-iのスマートメータ205は、時刻をtとするとき、需要家200-iの太陽光発電設備202の発電出力PPVi(t)と負荷201による消費電力PLi(t)とが合算された電力を計測する。すなわち、スマートメータ205が計測する電力は上述した残余需要であり、需要家200-2~200-nごとに、負荷201の消費電力から太陽光発電設備202の発電出力が差し引かれた量である。スマートメータ205が計測した残余需要の計測データは、通信網25と通信ネットワーク50を介して、太陽光発電出力推定装置100に送られる。 A smart meter (abbreviated as SM (Smart Meter) in the figure) 205 is installed in the consumers 200-2 to 200-n who have concluded a surplus purchase contract. Assuming that i is an arbitrary integer from 2 to n, the smart meter 205 of the consumer 200-i has a power output PPVi (t) of the photovoltaic power generation facility 202 of the consumer 200-i when the time is t. And the power consumption PLi (t) by the load 201 are added together to measure the total power. That is, the electric power measured by the smart meter 205 is the residual demand described above, which is the amount obtained by subtracting the power generation output of the photovoltaic power generation facility 202 from the power consumption of the load 201 for each consumer 200-2 to 200-n. .. The measurement data of the residual demand measured by the smart meter 205 is sent to the photovoltaic power generation output estimation device 100 via the communication network 25 and the communication network 50.

例えば、需要家200-2では、スマートメータ205が、太陽光発電設備202の発電出力PPV2(t)と負荷201による消費電力PL2(t)との合算値を残余需要として計測および記録する。そして、スマートメータ205は、需要家200-2の残余需要の計測データを、通信網25と通信ネットワーク50とを介して、太陽光発電出力推定装置100に送信する。このようにして太陽光発電出力推定装置100は、余剰買取契約対象の太陽光発電設備202を保有する需要家200-2~200-nに設置されたスマートメータ205の計測値である残余需要の計測データを取得する。 For example, in the consumer 200-2, the smart meter 205 measures and records the total value of the power generation output PPV2 (t) of the photovoltaic power generation facility 202 and the power consumption PL2 (t) of the load 201 as the residual demand. .. Then, the smart meter 205 transmits the measurement data of the residual demand of the consumer 200-2 to the photovoltaic power generation output estimation device 100 via the communication network 25 and the communication network 50. In this way, the photovoltaic power generation output estimation device 100 is the residual demand, which is the measured value of the smart meter 205 installed in the consumers 200-2 to 200-n who own the photovoltaic power generation equipment 202 subject to the surplus purchase contract. Acquire measurement data.

また、需要家200-3には、需要家200-3に設けられた太陽光発電設備202の発電出力PPV3(t)を計測する計測器206が設けられている。計測器206は、太陽光発電設備202を制御するパワーコンディショナなどの装置の一部であってもよいし、個別に設けられた計測器であってもよい。計測器206は、計測した発電出力PPV3(t)の計測データを、通信網25と通信ネットワーク50を介して、太陽光発電出力推定装置100に送る。なお、ここでは、需要家200-3に計測器206が設けられる例を説明するが、需要家200-2~200-nのうち計測器206が設けられる需要家は0軒であってもよいし、また、複数であってもよい。 Further, the consumer 200-3 is provided with a measuring instrument 206 for measuring the power generation output PPV3 (t) of the photovoltaic power generation facility 202 provided in the consumer 200-3. The measuring instrument 206 may be a part of a device such as a power conditioner for controlling the photovoltaic power generation facility 202, or may be an individually provided measuring instrument. The measuring instrument 206 sends the measured data of the measured power generation output P PV3 (t) to the photovoltaic power generation output estimation device 100 via the communication network 25 and the communication network 50. Here, an example in which the measuring instrument 206 is provided in the consumer 200-3 will be described, but among the consumers 200-2 to 200-n, the number of consumers in which the measuring instrument 206 is provided may be 0. However, there may be more than one.

全量買取契約を締結している需要家200-1には、スマートメータ203とスマートメータ204とが設置されている。スマートメータ203は、需要家200-1の負荷201の消費電力PL1(t)を計測および記録する。スマートメータ204は、需要家200-1の太陽光発電設備202の発電出力PPV1(t)を計測および記録する。スマートメータ204が計測した需要家200-1の太陽光発電設備202の発電出力の計測データは、通信網25と通信ネットワーク50とを介して太陽光発電出力推定装置100に送られる。このようにして太陽光発電出力推定装置100は、需要家200-1の太陽光発電設備202の発電出力PPV1(t)の計測データを取得する。スマートメータ203の計測データも、通信網25と通信ネットワーク50とを介して太陽光発電出力推定装置100に送られる。なお、スマートメータ203とスマートメータ204の少なくともどちらか一方の装置、あるいは、それらの装置の代わりに図示しない装置が、スマートメータ203の計測データとスマートメータ204の計測データとの合算である見かけ上の消費電力、つまり、残余需要を計測(あるいは、算出)し、計測(あるいは、算出)された残余需要を太陽光発電出力推定装置100に送信してもよい。 A smart meter 203 and a smart meter 204 are installed in the customer 200-1 who has concluded a purchase contract for all the products. The smart meter 203 measures and records the power consumption PL1 (t) of the load 201 of the consumer 200-1. The smart meter 204 measures and records the power output P PV1 (t) of the photovoltaic power generation facility 202 of the consumer 200-1. The measurement data of the power generation output of the photovoltaic power generation facility 202 of the consumer 200-1 measured by the smart meter 204 is sent to the photovoltaic power generation output estimation device 100 via the communication network 25 and the communication network 50. In this way, the photovoltaic power generation output estimation device 100 acquires the measurement data of the power generation output PPV1 (t) of the photovoltaic power generation facility 202 of the consumer 200-1. The measurement data of the smart meter 203 is also sent to the photovoltaic power generation output estimation device 100 via the communication network 25 and the communication network 50. It should be noted that at least one of the smart meter 203 and the smart meter 204, or a device (not shown) instead of those devices, is apparently the sum of the measurement data of the smart meter 203 and the measurement data of the smart meter 204. The power consumption, that is, the residual demand may be measured (or calculated), and the measured (or calculated) residual demand may be transmitted to the photovoltaic power generation output estimation device 100.

以上のように、需要家200-1~200-nのうち、需要家200-1,200-3については、スマートメータ204または計測器206により太陽光発電設備202の発電出力の計測データを得ることができる。スマートメータ204も計測器の一部であると考えることができるため、需要家200-1および需要家200-3はいずれも、太陽光発電設備202の発電出力の計測データを取得できる需要家である。一方、需要家200-1~200-nのうち、スマートメータ204および計測器206のいずれも設けられていない需要家では、消費電力と合算された残余需要は計測されているものの、太陽光発電設備202の発電出力は直接計測されていないため正確に発電出力を求めることが困難である。本実施の形態では、太陽光発電出力推定装置100が、需要家200-1および需要家200-3のうち少なくとも一方を用いて、スマートメータ204および計測器206のいずれも設けられていない需要家の発電出力を正確に推定する方法について説明する。以下、発電出力の推定において基準として用いられる需要家200-1,200-3の太陽光発電設備を第1太陽光発電設備とも呼び、発電出力の推定対象となる太陽光発電設備を第2太陽光発電設備とも呼ぶ。 As described above, among the consumers 200-1 to 200-n, for the consumers 200-1 and 200-3, the measurement data of the power generation output of the photovoltaic power generation facility 202 is obtained by the smart meter 204 or the measuring instrument 206. be able to. Since the smart meter 204 can also be considered to be a part of the measuring instrument, both the consumer 200-1 and the consumer 200-3 are consumers who can acquire the measurement data of the power generation output of the photovoltaic power generation facility 202. be. On the other hand, among the consumers 200-1 to 200-n, the consumers who do not have either the smart meter 204 or the measuring instrument 206 are measured for the residual demand combined with the power consumption, but the photovoltaic power generation is performed. Since the power generation output of the equipment 202 is not directly measured, it is difficult to accurately obtain the power generation output. In the present embodiment, the photovoltaic power generation output estimation device 100 uses at least one of the consumer 200-1 and the consumer 200-3, and the smart meter 204 and the measuring instrument 206 are not provided in the consumer. The method of accurately estimating the power generation output of the above will be described. Hereinafter, the photovoltaic power generation equipment of consumers 200-1, 200-3 used as a reference in the estimation of the power generation output is also referred to as the first solar power generation equipment, and the solar power generation equipment for which the power generation output is estimated is the second solar power generation equipment. Also called a photovoltaic power generation facility.

太陽光発電出力推定装置100は、第1太陽光発電設備の発電出力である第1太陽光発電出力を利用して、第2太陽光発電設備の予め定められた時点(以下、「推定時点」と称する)の太陽光発電出力である第2太陽光発電出力を推定する。推定時点には、現在の時点のみならず、過去または未来の時点が含まれる。太陽光発電出力推定装置100は、パーソナルコンピュータなどの汎用のコンピュータシステムがプログラムを実行することによって実現されても良いし、専用のコンピュータシステムによって実現されても良い。太陽光発電出力推定装置100の詳細な構成については、後述する。なお、第1太陽光発電設備と第2太陽光発電設備との間の距離は、予め定められた値以下である。予め定められた値は、例えば数kmであるがこれに限定されない。 The photovoltaic power generation output estimation device 100 utilizes the first photovoltaic power generation output, which is the power generation output of the first photovoltaic power generation facility, at a predetermined time point of the second photovoltaic power generation facility (hereinafter, “estimated time point””. The second photovoltaic power generation output, which is the photovoltaic power generation output of (referred to as), is estimated. Estimated time points include not only current time points but also past or future time points. The photovoltaic power generation output estimation device 100 may be realized by a general-purpose computer system such as a personal computer executing a program, or may be realized by a dedicated computer system. The detailed configuration of the photovoltaic power generation output estimation device 100 will be described later. The distance between the first photovoltaic power generation facility and the second photovoltaic power generation facility is not more than a predetermined value. The predetermined value is, for example, several km, but is not limited to this.

次に、太陽光発電出力推定装置100の機能について説明する。太陽光発電出力推定装置100は、基準となる第1太陽光発電出力と推定対象の需要家の残余需要とを用いて、推定対象の需要家の太陽光発電設備202の発電出力である第2太陽光発電出力を推定する。第1太陽光発電出力は、第1需要家に設けられる第1太陽光発電設備の発電出力であり、第2太陽光発電出力は、第2需要家に設けられる第2太陽光発電設備の発電出力である。なお、太陽光発電設備202の発電出力と負荷201の消費電力には、有効電力と無効電力とが存在する。しかし、一般的な太陽光発電設備、特に家庭用の太陽光発電設備は力率一定で運転されているため、有効電力を推定することで無効電力も容易に推定することが可能である。このため、以下に記載する、太陽光発電設備202の発電出力と負荷201の消費電力は、有効電力を指すものとする。 Next, the function of the photovoltaic power generation output estimation device 100 will be described. The photovoltaic power generation output estimation device 100 uses the reference first photovoltaic power generation output and the residual demand of the consumer to be estimated, and is the second power output of the photovoltaic power generation facility 202 of the consumer to be estimated. Estimate the PV output. The first photovoltaic power generation output is the power generation output of the first photovoltaic power generation facility installed in the first consumer, and the second photovoltaic power generation output is the power generation of the second photovoltaic power generation facility installed in the second consumer. The output. The power output of the photovoltaic power generation facility 202 and the power consumption of the load 201 include active power and ineffective power. However, since general photovoltaic power generation equipment, especially household photovoltaic power generation equipment, is operated at a constant power factor, it is possible to easily estimate the ineffective power by estimating the active power. Therefore, the power generation output of the photovoltaic power generation facility 202 and the power consumption of the load 201 described below refer to the active power.

ここで、太陽光発電出力推定装置100における太陽光発電設備202の発電出力の推定方法の概要について説明する。ここでは、発電出力の推定対象として需要家200-2の太陽光発電設備202を例に挙げて説明し、また推定において基準として用いられる太陽光発電出力として需要家200-1の太陽光発電設備202の発電出力を例にあげて説明する。すなわち、この例において、第1太陽光発電出力は需要家200-1の太陽光発電設備202の発電出力であり、推定対象の第2太陽光発電出力は、需要家200-2の太陽光発電設備202の発電出力である。なお、推定対象が、需要家200-4~200-nのそれぞれの太陽光発電設備202の発電出力である場合、すなわち第2太陽光発電出力が需要家200-4~200-nのそれぞれの太陽光発電設備202の発電出力である場合も、残余需要としてそれぞれの需要家に対応する値を用いることで、同様に太陽光発電設備202の発電出力を推定することができる。また、基準として用いられる太陽光発電出力である第1太陽光発電出力として、需要家200-1の太陽光発電設備202の発電出力の代わりに、需要家200-3の太陽光発電設備202の発電出力、すなわち計測器206により計測される発電出力を用いる場合も、同様に、第2太陽光発電出力を推定することができる。 Here, the outline of the method of estimating the power generation output of the photovoltaic power generation facility 202 in the photovoltaic power generation output estimation device 100 will be described. Here, the photovoltaic power generation facility 202 of the consumer 200-2 will be described as an example of the estimation target of the power generation output, and the photovoltaic power generation facility of the consumer 200-1 will be used as the reference in the estimation. The power generation output of 202 will be described as an example. That is, in this example, the first photovoltaic power generation output is the photovoltaic power generation output of the photovoltaic power generation facility 202 of the consumer 200-1, and the second photovoltaic power generation output to be estimated is the photovoltaic power generation of the consumer 200-2. It is the power generation output of the equipment 202. In addition, when the estimation target is the power generation output of each of the solar power generation facilities 202 of the consumers 200-4 to 200-n, that is, the second solar power generation output is each of the consumers 200-4 to 200-n. Even in the case of the power output of the photovoltaic power generation facility 202, the power output of the photovoltaic power generation facility 202 can be similarly estimated by using a value corresponding to each consumer as the residual demand. In addition, as the first photovoltaic power generation output, which is the photovoltaic power generation output used as a reference, instead of the photovoltaic power generation output of the photovoltaic power generation facility 202 of the consumer 200-1, the photovoltaic power generation facility 202 of the consumer 200-3 Similarly, when the power generation output, that is, the power generation output measured by the measuring instrument 206 is used, the second photovoltaic power generation output can be estimated.

第2需要家である需要家200-2の残余需要をP(t)、需要家200-2の負荷201の消費電力をPL2(t)、第2太陽光発電出力である需要家200-2の太陽光発電設備202の発電出力をPPV2(t)とする。また、第1需要家である需要家200-1の太陽光発電設備202の発電出力、すなわち第1太陽光発電出力をPPV1(t)とする。このとき、以下の式(1)のように、P(t)に関する式が成り立ち、さらに、式(2)のように、PPV2(t)に関する式を仮定できる。 The residual demand of the consumer 200-2, which is the second consumer, is P 2 (t), the power consumption of the load 201 of the consumer 200-2 is PL 2 (t), and the consumer 200, which is the output of the second photovoltaic power generation. Let the power generation output of the photovoltaic power generation facility 202 of -2 be P PV2 (t). Further, the power generation output of the photovoltaic power generation facility 202 of the consumer 200-1, which is the first consumer, that is, the first photovoltaic power generation output is PPV1 (t). At this time, the equation relating to P 2 (t) holds as in the following equation (1), and further, the equation relating to P PV2 (t) can be assumed as in the equation (2).

Figure 2022084398000002
Figure 2022084398000002

Figure 2022084398000003
Figure 2022084398000003

ここで、αは第1太陽光発電出力を第2太陽光発電出力に変換する係数である。ε12(t)は時間的外乱である。τは遅延時間、すなわち、第1太陽光発電設備の設置地点と第2太陽光発電設備の設置地点との間を日射変動が伝播する時間である。第1太陽光発電設備の設置地点と第2太陽光発電設備の設置地点との間を日射変動が伝播する時間、すなわち第1太陽光発電設備と第2太陽光発電設備との間を日射変動が伝播する時間を、第1遅延時間とも呼ぶ。 Here, α is a coefficient for converting the first photovoltaic power generation output into the second photovoltaic power generation output. ε 12 (t) is a temporal disturbance. τ s is the delay time, that is, the time during which the solar radiation fluctuation propagates between the installation point of the first photovoltaic power generation facility and the installation point of the second photovoltaic power generation facility. The time during which the solar radiation fluctuation propagates between the installation point of the first photovoltaic power generation facility and the installation point of the second photovoltaic power generation facility, that is, the solar radiation fluctuation between the first photovoltaic power generation facility and the second photovoltaic power generation facility. The time for propagation is also called the first delay time.

そして、PL2(t)とε12(t)の変動が、それぞれ、PPV1(t)の変動と相関がないと仮定し、上記の式(1)、式(2)において、時間的な定常性が成立すると仮定すると、上記の式(1)と式(2)から、以下の式(3)が導かれる。 Then, assuming that the fluctuations of PL 2 (t) and ε 12 (t) do not correlate with the fluctuations of P PV1 (t), respectively, in the above equations (1) and (2), temporally. Assuming that the steady state is established, the following equation (3) is derived from the above equations (1) and (2).

Figure 2022084398000004
Figure 2022084398000004

ここで、Cov[]は2つの量の共分散関数を計算する演算子を意味しており、上記の式(3)におけるCov[PPV1(t),PPV1(t+τ+τ)]は、τ=-τのときに最大値をとる。このため、第1太陽光発電出力と第2需要家の残余需要との時系列データの共分散関数Cov[PPV1(t),P(t+τ)]が最小値(符号が「-」かつ絶対値が最大値)をとるときのタイムラグに-1を乗じた値が、遅延時間τである。 Here, Cov t [] means an operator for calculating a covariance function of two quantities, and Cov t [P PV1 (t), P PV1 (t + τ + τ s )] in the above equation (3) is , Τ = -τ s takes the maximum value. Therefore, the covariance function Cov t [P PV1 (t), P2 ( t + τ)] of the time series data of the first photovoltaic power generation output and the residual demand of the second consumer is the minimum value (the sign is “-”. And the value obtained by multiplying the time lag when taking (the absolute value is the maximum value) by -1 is the delay time τ s .

また、上記の式(3)にτ=0を代入すると、以下の式(4)が得られる。 Further, by substituting τ = 0 into the above equation (3), the following equation (4) is obtained.

Figure 2022084398000005
Figure 2022084398000005

そして、上記の式(4)を以下の式(5)に変形することで、係数αの推定値を得ることができる。 Then, by transforming the above equation (4) into the following equation (5), an estimated value of the coefficient α can be obtained.

Figure 2022084398000006
Figure 2022084398000006

Cov[PPV1(t),P(t)]は、第1太陽光発電出力PPV1(t)と残余需要P(t)との共分散であり、この共分散を第1共分散と呼ぶ。Cov[PPV1(t),PPV1(t+τ)]は、遅延時間τずれた2つの第1太陽光発電出力であるPPV1(t)とPPV1(t+τ)との自己共分散であり、この自己共分散を第2共分散と呼ぶ。そして、αは、第1共分散Cov[PPV1(t),P(t)]を第2共分散Cov[PPV1(t),PPV1(t+τ)]で除して-1を乗じて得られる係数である。 Cov t [P PV1 (t), P2 ( t)] is a covariance between the first photovoltaic power generation output P PV1 (t) and the residual demand P2 ( t), and this covariance is the first covariance. Called variance. Cov t [P PV1 (t), P PV1 (t + τ s )] is the autocovariance of P PV1 (t) and P PV1 (t + τ s ), which are the two first photovoltaic outputs with a delay time of τ s . It is a variance, and this autocovariance is called a second covariance. Then, α is obtained by dividing the first covariance Cov t [P PV1 (t), P2 (t)] by the second covariance Cov t [P PV1 (t), P PV1 (t + τ s )]. It is a coefficient obtained by multiplying by 1.

また、ε12(t+τ)は微小であり無視できると仮定すれば、上記の式(2)から、第2太陽光発電出力に関する近似式である以下の式(6)を得ることができる。 Further, assuming that ε 12 (t + τ) is minute and can be ignored, the following equation (6), which is an approximate equation for the second photovoltaic power generation output, can be obtained from the above equation (2).

Figure 2022084398000007
Figure 2022084398000007

本実施の形態の太陽光発電出力推定装置100は、以上で述べた推定方法を用いて第2太陽光発電出力を推定する。すなわち、太陽光発電出力推定装置100は、上述した第1共分散、第2共分散を算出し、算出した第1共分散、第2共分散を用いて式(5)により、係数αの推定値を求める。係数αを求める処理は、第2太陽光発電出力の推定対象の時刻である推定時点より前に行われていればよい。ただし、推定時点が過去の場合は、推定時点の未来にあたる推定時点から現在までに行うことでもよい。上述したように、係数αの算出には、各計測データの時系列データを用いるため、ある程度の期間の計測データが必要である。この期間は、推定時点以前であって、かつ推定時点から時間的に離れすぎていない範囲で設定される。例えば、この期間は、第2太陽光発電出力を推定する推定時点から1,2週間程度遡った期間までの間の任意の期間である。例えば、係数αの推定値を求める処理は、第2太陽光発電出力を推定する日の前日などに行われるが、係数αの推定値を求める処理を行うタイミングはこの例に限定されない。そして、太陽光発電出力推定装置100は、係数αの推定値に、推定時点tから遅延時間τずれた時点における第1太陽光発電出力PPV1(t+τ)を乗じることで、第2太陽光発電出力PPV2(t)の推定値を算出することができる。 The photovoltaic power generation output estimation device 100 of the present embodiment estimates the second photovoltaic power generation output by using the estimation method described above. That is, the photovoltaic power generation output estimation device 100 calculates the first covariance and the second covariance described above, and estimates the coefficient α by the equation (5) using the calculated first covariance and second covariance. Find the value. The process for obtaining the coefficient α may be performed before the estimated time point, which is the time of the estimation target of the second photovoltaic power generation output. However, if the estimated time point is in the past, it may be performed from the estimated time point, which is the future of the estimated time point, to the present. As described above, since the time series data of each measurement data is used for the calculation of the coefficient α, the measurement data for a certain period is required. This period is set before the estimated time point and within a range not too far from the estimated time point. For example, this period is an arbitrary period from the estimation time at which the second photovoltaic power generation output is estimated to the period retroactive by about 1 or 2 weeks. For example, the process of obtaining the estimated value of the coefficient α is performed on the day before the day when the second photovoltaic power generation output is estimated, but the timing of performing the process of obtaining the estimated value of the coefficient α is not limited to this example. Then, the photovoltaic power generation output estimation device 100 multiplies the estimated value of the coefficient α by the first photovoltaic power generation output P PV1 (t + τ s ) at the time when the delay time τ s deviates from the estimation time point t, so that the second sun The estimated value of the photovoltaic power generation output P PV2 (t) can be calculated.

なお、上記の算出方法では、PL2(t)とPPV1(t)の相関が無い、すなわち、Cov[PPV1(t),PL2(t)]は、極めて小さく0に近似できると仮定したが、この仮定が成り立たない場合がある。例えば、日射強度が弱く第1太陽光発電出力が小さい値の場合や、日射強度の変動が少なく第1太陽光発電出力の変動が小さい場合には、Cov[PPV1(t),PPV1(t+τ)]が小さくなるため、相対的にCov[PPV1(t),PL2(t)]の値を0と近似して無視することができなくなり、上記の過程が成り立たず、このような場合、係数αの推定精度が低下することから第2太陽光発電出力PPV2(t)の推定精度も低下する。したがって、上記仮定が成り立たない条件の期間は、係数αの推定に適さない期間すなわち推定精度が低下する期間であると言える。 In the above calculation method, there is no correlation between PL2 (t) and PPV1 (t), that is, Cov t [P PV1 (t), PL2 (t)] is extremely small and can be approximated to 0. I made an assumption, but sometimes this assumption does not hold. For example, when the solar radiation intensity is weak and the fluctuation of the first photovoltaic power generation output is small, or when the fluctuation of the solar radiation intensity is small and the fluctuation of the first photovoltaic power generation output is small, Covert [P PV1 ( t ), P PV1 Since (t + τ s )] becomes smaller, the value of Cov t [P PV1 (t), PL2 (t)] cannot be ignored by approximating it to 0, and the above process does not hold. In such a case, the estimation accuracy of the coefficient α is lowered, so that the estimation accuracy of the second photovoltaic power generation output PPV2 (t) is also lowered. Therefore, it can be said that the period under the condition that the above assumption does not hold is a period that is not suitable for the estimation of the coefficient α, that is, a period in which the estimation accuracy is lowered.

本実施の形態では、係数αの推定に適さない期間の計測データを用いて推定されたαの推定値が第2太陽光発電出力PPV2(t)の推定に用いられることを避けるため、推定値に適している期間であるか否かを判定し、係数αの推定に適さない期間である場合には、当該期間の計測データを用いて算出されたαの推定値を採用しない。後述するように、例えば、係数αの初期値が予め定められ、係数αを求める処理が行われると、当該処理により推定された係数αの値で更新されていくが、係数αの推定に適さない期間であると判定された場合は、係数αを求める処理が行われないまたは係数αを求める処理により算出された値での更新が行われない。 In the present embodiment, the estimation value of α estimated by using the measurement data of the period unsuitable for the estimation of the coefficient α is estimated in order to avoid being used for the estimation of the second photovoltaic power generation output PPV2 (t). It is determined whether or not the period is suitable for the value, and if the period is not suitable for estimating the coefficient α, the estimated value of α calculated using the measurement data of the period is not adopted. As will be described later, for example, when the initial value of the coefficient α is predetermined and a process for obtaining the coefficient α is performed, the value of the coefficient α estimated by the process is updated, but it is suitable for estimating the coefficient α. If it is determined that there is no period, the process for obtaining the coefficient α is not performed or the update with the value calculated by the process for obtaining the coefficient α is not performed.

具体的には、太陽光発電出力推定装置100は、上記式(5)に示した係数αと同様の算出方法で、第1太陽光発電出力を推定するための係数である自己係数λを算出する。自己係数λは推定精度の検証に用いる検証係数である。詳細には、式(5)の分子の第1共分散として、第1太陽光発電出力PPV1(t)と第2需要家の残余需要P(t)との共分散の代わりに、第1太陽光発電出力PPV1(t)と第1需要家の残余需要P(t)との共分散を用い、式(5)の分母における遅延時間τを除いた以下の式(7)により、自己係数λを算出する。なお、本実施の形態では、このような、第1太陽光発電出力PPV1(t)と第1需要家の残余需要P(t)との共分散についても、第1共分散と呼ぶ。すなわち、第1共分散は、第1太陽光発電出力と残余需要との共分散であり、この残余需要は、第1需要家以外の需要家の残余需要の場合もあり第1需要家の残余需要の場合もある。また、式(5)における分母には遅延時間τが含まれ、式(7)における分母には遅延時間τが含まれないが、遅延時間τを0とすれば両者は同じとなるため、以下では、式(7)における分母も、式(5)における分母と同様に第2共分散と呼ぶ。 Specifically, the photovoltaic power generation output estimation device 100 calculates the self-coefficient λ, which is a coefficient for estimating the first photovoltaic power generation output, by the same calculation method as the coefficient α shown in the above equation (5). do. The self-coefficient λ is a verification coefficient used for verification of estimation accuracy. Specifically, as the first covariance of the molecule of the formula (5), instead of the covariance of the first photovoltaic output P PV1 (t) and the residual demand P2 (t) of the second consumer, the first covariance is performed. 1 Using the covariance of the photovoltaic output P PV1 (t) and the residual demand P1 (t) of the first consumer, the following equation (7) excluding the delay time τ s in the denominator of the equation (5) To calculate the self-factor λ. In the present embodiment, such a covariance between the first photovoltaic power generation output P PV1 (t) and the residual demand P1 (t) of the first consumer is also referred to as the first covariance. That is, the first covariance is the covariance between the output of the first photovoltaic power generation and the residual demand, and this residual demand may be the residual demand of consumers other than the first consumer, and the residual of the first consumer. It may be demand. Further, the denominator in the equation (5) includes the delay time τ s , and the denominator in the equation (7) does not include the delay time τ s , but if the delay time τ s is set to 0, both are the same. Therefore, in the following, the denominator in the equation (7) is also referred to as the second covariance like the denominator in the equation (5).

Figure 2022084398000008
Figure 2022084398000008

このように、自己係数λは、第1太陽光発電出力の計測データから第1太陽光発電出力自身を推定するための変換係数に相当する。理想的にはλは1である。実際にはλは各種の誤差等により1にはならないがλが1に近いほど誤差は少ないと考えられる。したがって、λと1の差の絶対値である|λ-1|が0に近いと、λの誤差は少ないと考えられ、λと同様の手法により推定される係数αについても、|λ-1|が0に近いと推定誤差は少ないと考えられる。このため、本実施の形態では、|λ-1|を用いて、係数αの推定に適さない期間であるか否かを判定する。このように、本実施の形態では、係数αの推定に適さない期間であるか否かの判定を行い、係数αの推定に適さない期間と判定された場合、当該期間の計測データを用いた係数αを採用しないため、第2太陽光発電出力PPV2(t)の推定精度の低下を防ぐことができる。太陽光発電出力推定装置100の処理の詳細については後述する。 As described above, the self-coefficient λ corresponds to the conversion coefficient for estimating the first photovoltaic power generation output itself from the measurement data of the first photovoltaic power generation output. Ideally λ is 1. Actually, λ does not become 1 due to various errors and the like, but it is considered that the closer λ is to 1, the smaller the error. Therefore, when | λ-1 |, which is the absolute value of the difference between λ and 1, is close to 0, the error of λ is considered to be small, and the coefficient α estimated by the same method as λ is also | λ-1. When | is close to 0, the estimation error is considered to be small. Therefore, in the present embodiment, | λ-1 | is used to determine whether or not the period is not suitable for estimating the coefficient α. As described above, in the present embodiment, it is determined whether or not the period is not suitable for the estimation of the coefficient α, and if it is determined that the period is not suitable for the estimation of the coefficient α, the measurement data of the period is used. Since the coefficient α is not adopted, it is possible to prevent a decrease in the estimation accuracy of the second photovoltaic power generation output PPV2 (t). The details of the processing of the photovoltaic power generation output estimation device 100 will be described later.

次に、太陽光発電出力推定装置100の機能構成例について説明する。図2は、本実施の形態の太陽光発電出力推定装置100の機能構成例を示すブロック図である。太陽光発電出力推定装置100は、太陽光発電設備202の推定時点の発電出力である第2太陽光発電出力を推定する装置である。図2に示すように、太陽光発電出力推定装置100は、データ取得部110と、第1共分散算出部120と、第2共分散算出部130と、推定精度検証部140と、推定部150と、表示部160と、入力受付部170と、記憶部180とを備える。 Next, an example of the functional configuration of the photovoltaic power generation output estimation device 100 will be described. FIG. 2 is a block diagram showing a functional configuration example of the photovoltaic power generation output estimation device 100 of the present embodiment. The photovoltaic power generation output estimation device 100 is a device that estimates the second photovoltaic power generation output, which is the power generation output at the time of estimation of the photovoltaic power generation facility 202. As shown in FIG. 2, the photovoltaic power generation output estimation device 100 includes a data acquisition unit 110, a first covariance calculation unit 120, a second covariance calculation unit 130, an estimation accuracy verification unit 140, and an estimation unit 150. A display unit 160, an input reception unit 170, and a storage unit 180 are provided.

データ取得部110は、通信網25と通信ネットワーク50とを介して、スマートメータ203~205および計測器206から計測データを受信し、受信した計測データを記憶部180に格納する。詳細には、スマートメータ204および計測器206から受け取った計測データは太陽光発電出力データ182として記憶部180に格納され、スマートメータ205から受け取った計測データは残余需要データ183として記憶部180に格納される。これらの計測データは、計測値ごとに時刻と対応づけられている。なお、記憶部180に格納されるこれらの計測データの各計測値の時間間隔は、例えば1秒、10秒、1分、30分、または1時間などであるが、特にこれらには限定されず、ユーザによって適切な数値が定められる。また、これらの計測データは、需要家200-1~200-nの識別情報または太陽光発電設備202の識別情報とともに格納される。また、データ取得部110は、需要家200-1のスマートメータ203の計測データとスマートメータ204の計測データとを受信した場合は、これらの計測データを用いて需要家200-1の残余需要を算出して残余需要データ183として記憶部180に格納する。 The data acquisition unit 110 receives measurement data from the smart meters 203 to 205 and the measuring instrument 206 via the communication network 25 and the communication network 50, and stores the received measurement data in the storage unit 180. Specifically, the measurement data received from the smart meter 204 and the measuring instrument 206 is stored in the storage unit 180 as the photovoltaic power generation output data 182, and the measurement data received from the smart meter 205 is stored in the storage unit 180 as the residual demand data 183. Will be done. These measurement data are associated with the time for each measured value. The time interval of each measured value of these measurement data stored in the storage unit 180 is, for example, 1 second, 10 seconds, 1 minute, 30 minutes, or 1 hour, but is not particularly limited thereto. , The user determines the appropriate number. Further, these measurement data are stored together with the identification information of the consumers 200-1 to 200-n or the identification information of the photovoltaic power generation facility 202. Further, when the data acquisition unit 110 receives the measurement data of the smart meter 203 of the consumer 200-1 and the measurement data of the smart meter 204, the data acquisition unit 110 uses these measurement data to obtain the residual demand of the consumer 200-1. It is calculated and stored in the storage unit 180 as residual demand data 183.

なお、ここでは、太陽光発電出力推定装置100が、通信網25と通信ネットワーク50とを介してスマートメータ203~205および計測器206から計測データを取得する例を示すが、スマートメータ203~205の計測データは、いわゆる集約装置またはヘッドエンドシステムなどにより収集され、収集した装置から太陽光発電出力推定装置100へ送信されてもよい。すなわち、通信網25および通信ネットワーク50内に、計測データを収集する装置が存在し、当該装置から太陽光発電出力推定装置100が計測データを取得してもよい。また、太陽光発電出力推定装置100とスマートメータ203~205および計測器206との間の通信ルートは図1および図2に示した例に限定されず、太陽光発電出力推定装置100が計測器206を受信する通信ルートと、太陽光発電出力推定装置100がスマートメータ203~205の計測データを受信する通信ルートとは異なっていてもよい。また、太陽光発電出力推定装置100とスマートメータ203~205および計測器206との間の通信ネットワークは、通信網25と通信ネットワーク50とに分離されていなくてもよい。また、太陽光発電出力推定装置100は、複数の通信ルートでスマートメータ203~205の計測データを受信してもよい。 Here, an example is shown in which the photovoltaic power generation output estimation device 100 acquires measurement data from the smart meters 203 to 205 and the measuring instrument 206 via the communication network 25 and the communication network 50, but the smart meters 203 to 205 are shown. The measurement data of the above may be collected by a so-called aggregation device, a head-end system, or the like, and may be transmitted from the collected device to the photovoltaic power generation output estimation device 100. That is, there may be a device for collecting measurement data in the communication network 25 and the communication network 50, and the photovoltaic power generation output estimation device 100 may acquire the measurement data from the device. Further, the communication route between the photovoltaic power generation output estimation device 100 and the smart meters 203 to 205 and the measuring instrument 206 is not limited to the examples shown in FIGS. 1 and 2, and the photovoltaic power generation output estimation device 100 is the measuring instrument. The communication route for receiving 206 may be different from the communication route for the photovoltaic power generation output estimation device 100 to receive the measurement data of the smart meters 203 to 205. Further, the communication network between the photovoltaic power generation output estimation device 100 and the smart meters 203 to 205 and the measuring instrument 206 does not have to be separated into the communication network 25 and the communication network 50. Further, the photovoltaic power generation output estimation device 100 may receive the measurement data of the smart meters 203 to 205 through a plurality of communication routes.

記憶部180は、上述した太陽光発電出力データ182および残余需要データ183を格納するとともに、太陽光発電出力推定装置100において算出された各種のデータを算出データ185として記憶する。また、記憶部180には、例えば、ユーザにより入力受付部170を介して入力された場所情報データ181、係数初期値データ184が格納される。なお、場所情報データ181、係数初期値データ184は、入力受付部170を介して入力される代わりに他の装置から送信され、データ取得部110によって取得されて記憶部180に格納されてもよい。場所情報データ181は、太陽光発電設備202を保有する需要家200-1~200-nの住所または緯度経度などの場所情報のデータである。 The storage unit 180 stores the above-mentioned photovoltaic power generation output data 182 and residual demand data 183, and stores various data calculated by the photovoltaic power generation output estimation device 100 as calculated data 185. Further, the storage unit 180 stores, for example, the location information data 181 and the coefficient initial value data 184 input by the user via the input reception unit 170. The location information data 181 and the coefficient initial value data 184 may be transmitted from another device instead of being input via the input receiving unit 170, acquired by the data acquisition unit 110, and stored in the storage unit 180. .. The location information data 181 is data on location information such as the address or latitude / longitude of consumers 200-1 to 200-n who own the photovoltaic power generation facility 202.

第1共分散算出部120は、記憶部180の場所情報データ181を参照して、発電出力の推定対象の需要家である第2需要家との距離が、予め定められた距離以内であって推定の際に基準となる太陽光発電設備202を有する第1需要家を選択する。第1需要家は、上述したように、全量買取契約の対象の太陽光発電設備202を有する需要家200-1、または計測器206により太陽光発電設備202の発電出力が計測されている需要家200-3である。第1共分散算出部120は、第1需要家の太陽光発電設備202の発電出力である第1太陽光発電出力と残余需要との共分散である第1共分散を算出し、算出した第1共分散を記憶部180に算出データ185として格納する。第1共分散は、上述したとおり、第1需要家以外の需要家すなわち発電出力の推定対象の需要家の残余需要の場合もあり、第1需要家の残余需要の場合もある。発電出力の推定対象の需要家は、上述したように例えば需要家200-2であるが、需要家200-4~200-nであってもよい。需要家200-2,200-4~200-nの全てを順次推定対象として、需要家200-2,200-4~200-nごとに推定が行われてもよいし、その合計値が必要な場合はその合計値を一括で推定することにしてもよい。 The first covariance calculation unit 120 refers to the location information data 181 of the storage unit 180, and the distance from the second consumer, who is the consumer whose power generation output is estimated, is within a predetermined distance. The first consumer who has the photovoltaic power generation facility 202 as a reference at the time of estimation is selected. As described above, the first consumer is the consumer 200-1 having the photovoltaic power generation facility 202 subject to the total purchase contract, or the consumer whose power generation output of the photovoltaic power generation facility 202 is measured by the measuring instrument 206. It is 200-3. The first covariance calculation unit 120 calculates and calculates the first covariance, which is the covariance between the first photovoltaic power generation output, which is the power output of the photovoltaic power generation facility 202 of the first consumer, and the residual demand. 1 The covariance is stored in the storage unit 180 as calculated data 185. As described above, the first covariance may be the residual demand of a consumer other than the first consumer, that is, the consumer whose power generation output is estimated, or the residual demand of the first consumer. As described above, the consumer whose power generation output is estimated is, for example, the consumer 200-2, but may be the consumer 200-4 to 200-n. All of the consumers 200-2, 200-4 to 200-n may be sequentially estimated, and the estimation may be performed for each of the consumers 200-2, 200-4 to 200-n, and the total value is required. If this is the case, the total value may be estimated in a lump sum.

具体的には、第1共分散算出部120は、まず、上述した自己係数λの算出のために、第1期間における第1太陽光発電出力を記憶部180の太陽光発電出力データ182から抽出して読み出し、第1需要家の残余需要を記憶部180の残余需要データ183から抽出して読み出す。第1共分散算出部120は、読み出したデータを用いて第1太陽光発電出力と第1需要家の残余需要との共分散である第1共分散を算出して、記憶部180に算出データ185として格納する。第1期間は、推定時点から1,2週間程度遡った期間内の、例えば、6時間から10時間程度の期間であり、好ましくは8時間である。第1期間は、日射強度が強く、かつ日射強度の変動が大きい期間であるのが好ましい。ただし、本実施の形態では、第1期間が係数αの算出に適した期間であるか否かを判定するため、第1期間を日射強度によらずに選択したとしても、日射強度が弱い期間、日射強度の変動が小さい期間などが第1期間として設定されても、当該第1期間に対応する係数αは用いられないため、係数αの推定精度の低下を抑制することができる。このため、第1期間の設定時に、日射強度が考慮されていなくてもよい。 Specifically, the first covariance calculation unit 120 first extracts the first photovoltaic power generation output in the first period from the photovoltaic power generation output data 182 of the storage unit 180 in order to calculate the self-factor λ described above. Then, the residual demand of the first consumer is extracted from the residual demand data 183 of the storage unit 180 and read out. The first covariance calculation unit 120 calculates the first covariance, which is the covariance between the first photovoltaic power generation output and the residual demand of the first consumer, using the read data, and the calculated data is stored in the storage unit 180. Store as 185. The first period is, for example, a period of about 6 to 10 hours, preferably about 8 hours, within a period of about 1 to 2 weeks from the estimation time. The first period is preferably a period in which the solar radiation intensity is strong and the fluctuation of the solar radiation intensity is large. However, in the present embodiment, in order to determine whether or not the first period is a period suitable for calculating the coefficient α, even if the first period is selected regardless of the solar radiation intensity, the solar radiation intensity is weak. Even if a period in which the fluctuation of the solar radiation intensity is small is set as the first period, the coefficient α corresponding to the first period is not used, so that the deterioration of the estimation accuracy of the coefficient α can be suppressed. Therefore, it is not necessary to consider the solar radiation intensity when setting the first period.

また、第1共分散算出部120は、推定精度検証部140から係数αの算出のための第1共分散の算出を指示されると、第1期間の第2需要家の残余需要を記憶部180の残余需要データ183から抽出して読み出す。第1共分散算出部120は、読み出し済の第1太陽光発電出力と読み出した第2需要家の残余需要とを用いて、第1期間における第1太陽光発電出力と第2需要家の残余需要との共分散である第1共分散を算出して、記憶部180に算出データ185として格納する。さらに、第1共分散算出部120は、読み出し済の第1太陽光発電出力と、同じく読み出し済の第2需要家の残余需要とを用いて、第1期間における遅延時間τを算出し、記憶部180に算出データ185として格納する。なお、遅延時間は算出せずに、0と近似することにしてもよい。 Further, when the first covariance calculation unit 120 is instructed by the estimation accuracy verification unit 140 to calculate the first covariance for calculating the coefficient α, the first covariance calculation unit 120 stores the residual demand of the second consumer in the first period. It is extracted and read from the residual demand data 183 of 180. The first covariance calculation unit 120 uses the read-out first photovoltaic power generation output and the read-out residual demand of the second consumer to generate the first photovoltaic power generation output and the residual demand of the second consumer in the first period. The first covariance, which is the covariance with the demand, is calculated and stored in the storage unit 180 as the calculated data 185. Further, the first covariance calculation unit 120 calculates the delay time τ s in the first period by using the read-out first photovoltaic power generation output and the read-out residual demand of the second consumer. It is stored as calculated data 185 in the storage unit 180. The delay time may not be calculated and may be approximated to 0.

第2共分散算出部130は、自己係数λの算出のための第2共分散の算出では、第1期間における第1太陽光発電出力の時系列データを記憶部180の算出データ185から抽出して読み出し、読み出した第1太陽光発電出力の時系列データの分散を、第2共分散として算出する。 In the calculation of the second covariance for calculating the self-covariance λ, the second covariance calculation unit 130 extracts the time series data of the first solar power generation output in the first period from the calculation data 185 of the storage unit 180. The variance of the time-series data of the first solar power generation output read out is calculated as the second covariance.

また、第2共分散算出部130は、係数αの算出のための第2共分散の算出処理では、遅延時間τを記憶部180の算出データ185から抽出して読み出し、第1期間における第1太陽光発電出力と、第1期間より遅延時間τずれた第2期間における第1太陽光発電出力とを記憶部180の太陽光発電出力データ182から抽出して読み出す。そして、第2共分散算出部130は、遅延時間τずれた2つの第1太陽光発電出力の自己共分散である第2共分散を算出し、算出した第2共分散を記憶部180に算出データ185として格納する。詳細には、第2共分散算出部130は、記憶部180の太陽光発電出力データ182から、第1期間の第1太陽光発電出力の時系列データと、第1期間より遅延時間τずれた第2期間の第1太陽光発電出力の時系列データとを読み出す。そして、第2共分散算出部130は、第1期間における第1太陽光発電出力の時系列データと第2期間の第1太陽光発電出力の時系列データとの自己共分散である第2共分散を算出し、算出した第2共分散を記憶部180に算出データ185として格納する。なお、係数αの算出のための処理であって遅延時間τを0と近似する場合は、自己係数λの算出のための第2共分散の算出処理と同様に、第2共分散算出部130は、第1期間における第1太陽光発電出力の時系列データの分散を、第2共分散として算出する。 Further, in the calculation process of the second covariance for calculating the coefficient α, the second covariance calculation unit 130 extracts and reads the delay time τ s from the calculation data 185 of the storage unit 180, and reads the delay time τ s from the calculation data 185 of the storage unit 180. 1 The photovoltaic power generation output and the first photovoltaic power generation output in the second period deviated by the delay time τ s from the first period are extracted from the photovoltaic power generation output data 182 of the storage unit 180 and read out. Then, the second covariance calculation unit 130 calculates the second covariance, which is the autocovariance of the two first photovoltaic power generation outputs deviated by the delay time τ s , and stores the calculated second covariance in the storage unit 180. It is stored as calculated data 185. Specifically, the second co-dispersion calculation unit 130 deviates from the photovoltaic power generation output data 182 of the storage unit 180 with the time-series data of the first photovoltaic power generation output in the first period and the delay time τ s from the first period. The time-series data of the first photovoltaic power generation output in the second period is read out. Then, the second covariance calculation unit 130 is a second covariance which is an autocovariance between the time series data of the first solar power generation output in the first period and the time series data of the first solar power generation output in the second period. The variance is calculated, and the calculated second covariance is stored in the storage unit 180 as the calculated data 185. When the delay time τ s is approximated to 0 in the process for calculating the coefficient α, the second covariance calculation unit is the same as the process for calculating the self-variance λ. Reference numeral 130 calculates the variance of the time-series data of the first solar power generation output in the first period as the second covariance.

推定精度検証部140は、第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定する。詳細には、推定精度検証部140は、記憶部180に格納されている自己係数λの算出のための第1共分散および第2共分散を用いて、上述した自己係数λを算出し、自己係数λを用いて第1期間が係数αの算出に適した期間であるか否かを判定する。推定精度検証部140は、第1期間が係数αの算出に適した期間であると判定した場合、第1共分散算出部120へ係数αの算出のための第1共分散の算出を指示する。 The estimation accuracy verification unit 140 determines whether or not the estimation accuracy of the second photovoltaic power generation output satisfies a predetermined accuracy. Specifically, the estimation accuracy verification unit 140 calculates the self-coefficient λ described above by using the first covariance and the second covariance for calculating the self-coefficient λ stored in the storage unit 180, and self. Using the coefficient λ, it is determined whether or not the first period is a period suitable for calculating the coefficient α. When the estimation accuracy verification unit 140 determines that the first period is a period suitable for calculating the coefficient α, the estimation accuracy verification unit 140 instructs the first covariance calculation unit 120 to calculate the first covariance for calculating the coefficient α. ..

推定部150は、第1共分散および第2共分散を用いて第1太陽光発電出力と第2太陽光発電出力との比の推定値である係数αを算出し、係数αを推定時点から遅延時間ずれた第1太陽光発電出力に乗じることで第2太陽光発電出力を推定する。詳細には、推定部150は、記憶部180に算出データ185として格納されている係数αの算出のための第1共分散および第2共分散を読み出し、読み出した第1共分散および第2共分散を用いて発電出力の推定対象の需要家ごとの係数αを上記式(5)により推定する。推定部150は、算出データ185として格納されている係数αの値を推定値で更新する。また、推定部150は、推定時点から遅延時間τずれた時点の第1太陽光発電出力を記憶部180の太陽光発電出力データ182から抽出して読み出し、第1太陽光発電出力に係数αを乗じることで、推定対象の需要家の太陽光発電出力である第2太陽光発電出力を推定し、推定結果を算出データ185として記憶部180に格納する。なお、係数αは1,2週間などの短期間では変化しにくいものの、遅延時間τは風向または風速の変化によって短期間で変化しやすい。従って、この第2太陽光発電出力を推定する際の遅延時間τは、気象データなどを考慮してリアルタイムで変更されたものが用いられるのが好ましい。 The estimation unit 150 calculates a coefficient α, which is an estimated value of the ratio between the first photovoltaic power generation output and the second photovoltaic power generation output, using the first co-dispersion and the second co-dispersion, and calculates the coefficient α from the time of estimation. The second photovoltaic power generation output is estimated by multiplying the first photovoltaic power generation output with a delay time deviation. Specifically, the estimation unit 150 reads out the first covariance and the second covariance for calculating the coefficient α stored as the calculation data 185 in the storage unit 180, and reads out the first covariance and the second covariance. Using the variance, the coefficient α for each consumer whose power output is estimated is estimated by the above equation (5). The estimation unit 150 updates the value of the coefficient α stored as the calculated data 185 with the estimated value. Further, the estimation unit 150 extracts and reads out the first photovoltaic power generation output at the time when the delay time τ s deviates from the estimation time from the photovoltaic power generation output data 182 of the storage unit 180, and has a coefficient α for the first photovoltaic power generation output. By multiplying by, the second photovoltaic power generation output, which is the photovoltaic power generation output of the consumer to be estimated, is estimated, and the estimation result is stored in the storage unit 180 as calculation data 185. Although the coefficient α is unlikely to change in a short period such as 1 or 2 weeks, the delay time τ s tends to change in a short period due to a change in wind direction or speed. Therefore, it is preferable that the delay time τ s when estimating the second photovoltaic power generation output is changed in real time in consideration of meteorological data and the like.

表示部160は、記憶部180に格納された各種データを表示する。入力受付部170は、ユーザからの入力を受け付ける。 The display unit 160 displays various data stored in the storage unit 180. The input receiving unit 170 receives input from the user.

ここで、太陽光発電出力推定装置100のハードウェア構成について説明する。本実施の形態の太陽光発電出力推定装置100は、コンピュータシステム上で、太陽光発電出力推定装置100における処理が記述されたプログラムである太陽光発電出力推定プログラムが実行されることにより、コンピュータシステムが太陽光発電出力推定装置100として機能する。図3は、本実施の形態の太陽光発電出力推定装置100を実現するコンピュータシステムの構成例を示す図である。図3に示すように、このコンピュータシステムは、制御部101と入力部102と記憶部103と表示部104と通信部105と出力部106とを備え、これらはシステムバス107を介して接続されている。 Here, the hardware configuration of the photovoltaic power generation output estimation device 100 will be described. The photovoltaic power generation output estimation device 100 of the present embodiment is a computer system by executing a photovoltaic power generation output estimation program, which is a program in which the processing in the photovoltaic power generation output estimation device 100 is described, on the computer system. Functions as the photovoltaic power generation output estimation device 100. FIG. 3 is a diagram showing a configuration example of a computer system that realizes the photovoltaic power generation output estimation device 100 of the present embodiment. As shown in FIG. 3, this computer system includes a control unit 101, an input unit 102, a storage unit 103, a display unit 104, a communication unit 105, and an output unit 106, which are connected via a system bus 107. There is.

図3において、制御部101は、例えば、CPU(Central Processing Unit)等のプロセッサであり、本実施の形態の太陽光発電出力推定装置100における処理が記述された太陽光発電出力推定プログラムを実行する。入力部102は、たとえばキーボード、マウスなどで構成され、コンピュータシステムの使用者が、各種情報の入力を行うために使用する。記憶部103は、RAM(Random Access Memory),ROM(Read Only Memory)などの各種メモリおよびハードディスクなどのストレージデバイスを含み、上記制御部101が実行すべきプログラム、処理の過程で得られた必要なデータ、などを記憶する。また、記憶部103は、プログラムの一時的な記憶領域としても使用される。表示部104は、ディスプレイ、LCD(液晶表示パネル)などで構成され、コンピュータシステムの使用者に対して各種画面を表示する。通信部105は、通信処理を実施する受信機および送信機である。出力部106は、プリンタなどである。なお、図3は、一例であり、コンピュータシステムの構成は図3の例に限定されない。 In FIG. 3, the control unit 101 is, for example, a processor such as a CPU (Central Processing Unit), and executes a photovoltaic power generation output estimation program in which processing in the photovoltaic power generation output estimation device 100 of the present embodiment is described. .. The input unit 102 is composed of, for example, a keyboard, a mouse, or the like, and is used by a user of a computer system to input various information. The storage unit 103 includes various memories such as RAM (Random Access Memory) and ROM (Read Only Memory) and a storage device such as a hard disk, and is a necessary program obtained in the process of a program to be executed by the control unit 101. Memorize data, etc. The storage unit 103 is also used as a temporary storage area for the program. The display unit 104 is composed of a display, an LCD (liquid crystal display panel), and the like, and displays various screens to a user of a computer system. The communication unit 105 is a receiver and a transmitter that performs communication processing. The output unit 106 is a printer or the like. Note that FIG. 3 is an example, and the configuration of the computer system is not limited to the example of FIG.

ここで、本実施の形態の太陽光発電出力推定プログラムが実行可能な状態になるまでのコンピュータシステムの動作例について説明する。上述した構成をとるコンピュータシステムには、たとえば、図示しないCD(Compact Disc)-ROMドライブまたはDVD(Digital Versatile Disc)-ROMドライブにセットされたCD-ROMまたはDVD-ROMから、太陽光発電出力推定プログラムが記憶部103にインストールされる。そして、太陽光発電出力推定プログラムの実行時に、記憶部103から読み出された太陽光発電出力推定プログラムが記憶部103に格納される。この状態で、制御部101は、記憶部103に格納されたプログラムに従って、本実施の形態の太陽光発電出力推定装置100としての処理を実行する。 Here, an operation example of the computer system until the photovoltaic power generation output estimation program of the present embodiment becomes executable will be described. In a computer system having the above-described configuration, for example, solar power output estimation is performed from a CD-ROM or DVD-ROM set in a CD (Compact Disc) -ROM drive or DVD (Digital Versatile Disc) -ROM drive (not shown). The program is installed in the storage unit 103. Then, when the photovoltaic power generation output estimation program is executed, the photovoltaic power generation output estimation program read from the storage unit 103 is stored in the storage unit 103. In this state, the control unit 101 executes the process as the photovoltaic power generation output estimation device 100 of the present embodiment according to the program stored in the storage unit 103.

なお、上記の説明においては、CD-ROMまたはDVD-ROMを記録媒体として、太陽光発電出力推定装置100における処理を記述したプログラムを提供しているが、これに限らず、コンピュータシステムの構成、提供するプログラムの容量などに応じて、たとえば、通信部105を経由してインターネットなどの伝送媒体により提供されたプログラムを用いることとしてもよい。 In the above description, a program describing the processing in the solar power generation output estimation device 100 using a CD-ROM or a DVD-ROM as a recording medium is provided, but the present invention is not limited to this, and the configuration of the computer system, Depending on the capacity of the program to be provided, for example, a program provided by a transmission medium such as the Internet via the communication unit 105 may be used.

図2に示した第1共分散算出部120、第2共分散算出部130、推定精度検証部140および推定部150は、図3に示した記憶部103に記憶された太陽光発電出力推定プログラムが図3に示した制御部101により実行されることにより実現される。図2に示した記憶部180は、図3に示した記憶部103の一部である。図2に示したデータ取得部110は、図3に示した通信部105および制御部101により実現される。図2に示した表示部160は、図3に示した表示部104により実現され、図2に示した入力受付部170は、図3に示した入力部102により実現される。太陽光発電出力推定装置100は複数のコンピュータシステムにより実現されてもよい。 The first covariance calculation unit 120, the second covariance calculation unit 130, the estimation accuracy verification unit 140, and the estimation unit 150 shown in FIG. 2 are photovoltaic power generation output estimation programs stored in the storage unit 103 shown in FIG. Is realized by being executed by the control unit 101 shown in FIG. The storage unit 180 shown in FIG. 2 is a part of the storage unit 103 shown in FIG. The data acquisition unit 110 shown in FIG. 2 is realized by the communication unit 105 and the control unit 101 shown in FIG. The display unit 160 shown in FIG. 2 is realized by the display unit 104 shown in FIG. 3, and the input reception unit 170 shown in FIG. 2 is realized by the input unit 102 shown in FIG. The photovoltaic power generation output estimation device 100 may be realized by a plurality of computer systems.

例えば、本実施の形態の太陽光発電出力推定プログラムは、第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、第1太陽光発電設備の発電出力である第1太陽光発電出力と第2太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて推定するステップをコンピュータに実行させる。さらに、本実施の形態の太陽光発電出力推定プログラムは、第1太陽光発電設備の発電出力と第1太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定するステップ(判定ステップ)、をコンピュータに実行させる。 For example, the photovoltaic power generation output estimation program of the present embodiment calculates the second photovoltaic power generation output, which is the photovoltaic power generation output of the second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility. , A computer that estimates the steps using the first photovoltaic power generation output, which is the power generation output of the first photovoltaic power generation facility, and the residual demand, which is the apparent power consumption of the consumer in which the second photovoltaic power generation facility is installed. To execute. Further, the photovoltaic power generation output estimation program of the present embodiment uses the power output of the first photovoltaic power generation facility and the residual demand which is the apparent power consumption of the consumer in which the first photovoltaic power generation facility is installed. The second step (determination step) of determining whether or not the estimation accuracy of the photovoltaic power generation output satisfies the predetermined accuracy is executed by the computer.

次に、本実施の形態の太陽光発電出力推定装置100の動作の詳細について説明する。図4は、本実施の形態の太陽光発電出力推定装置100による係数αの推定処理手順の一例を示すフローチャートである。以下、図4のフローに沿って、太陽光発電出力推定装置100が第2太陽光発電出力を推定する処理について説明する。まず、第1共分散算出部120は、係数αの推定対象の太陽光発電設備202を設定し、発電出力を用いる太陽光発電設備を選択する(ステップS1)。係数αの推定対象の太陽光発電設備202が複数存在する場合は、複数の太陽光発電設備202の中から1つの太陽光発電設備202を選択する。なお、2回目以降のステップS1では、複数の太陽光発電設備202の中から、ステップS2~ステップS6の処理が行われていない太陽光発電設備202を、係数αの推定対象の太陽光発電設備202として設定する。なお、このステップS1は、推定部150により行われて、結果が第1共分散算出部120に通知されてもよい。 Next, the details of the operation of the photovoltaic power generation output estimation device 100 of the present embodiment will be described. FIG. 4 is a flowchart showing an example of an estimation processing procedure for the coefficient α by the photovoltaic power generation output estimation device 100 of the present embodiment. Hereinafter, the process of estimating the second photovoltaic power generation output by the photovoltaic power generation output estimation device 100 will be described along with the flow of FIG. First, the first covariance calculation unit 120 sets the photovoltaic power generation facility 202 to be estimated with the coefficient α, and selects the photovoltaic power generation facility using the power generation output (step S1). When there are a plurality of photovoltaic power generation facilities 202 for which the coefficient α is estimated, one photovoltaic power generation facility 202 is selected from the plurality of photovoltaic power generation facilities 202. In the second and subsequent steps S1, among the plurality of photovoltaic power generation facilities 202, the photovoltaic power generation facility 202 that has not been processed in steps S2 to S6 is selected as the photovoltaic power generation facility for which the coefficient α is estimated. Set as 202. Note that this step S1 may be performed by the estimation unit 150, and the result may be notified to the first covariance calculation unit 120.

ステップS1では、具体的には、第1共分散算出部120は、記憶部180の場所情報データ181を参照して、発電出力の推定対象の需要家である第2需要家との距離が、予め定められた距離以内であって推定の際に基準となる太陽光発電設備202を有する第1需要家を選択する。 In step S1, specifically, the first covariance calculation unit 120 refers to the location information data 181 of the storage unit 180, and the distance from the second consumer who is the consumer whose power generation output is estimated is determined. A first consumer who is within a predetermined distance and has a photovoltaic power generation facility 202 as a reference for estimation is selected.

図5は、本実施の形態の場所情報データ181により示される各需要家の太陽光発電設備202の地理的位置を模式的に示す図である。図5では、各需要家の太陽光発電設備202の地理的位置が当該需要家の地理的位置と一致すると仮定している。図5では、需要家200-i(i=1,2,…,8)の太陽光発電設備202を、PViとして示している。例えば、PV1は、需要家200-1の太陽光発電設備202であり、PV2は、需要家200-2の太陽光発電設備202である。実線の枠で示したPV1,PV3は、それぞれ需要家200-1,200-3の太陽光発電設備202であり、発電出力が計測されている。破線の枠で示したPV2などは、それぞれ需要家200-2など発電出力が計測されていない太陽光発電設備202である。 FIG. 5 is a diagram schematically showing the geographical position of the photovoltaic power generation facility 202 of each consumer shown by the location information data 181 of the present embodiment. In FIG. 5, it is assumed that the geographical position of the photovoltaic power generation facility 202 of each consumer coincides with the geographical position of the consumer. In FIG. 5, the photovoltaic power generation equipment 202 of the consumer 200-i (i = 1, 2, ..., 8) is shown as PVi. For example, PV1 is the photovoltaic power generation facility 202 of the consumer 200-1, and PV2 is the photovoltaic power generation facility 202 of the consumer 200-2. PV1 and PV3 shown in the solid line frame are the photovoltaic power generation facilities 202 of the consumers 200-1 and 200-3, respectively, and the power generation output is measured. PV2 and the like shown by the broken line frame are photovoltaic power generation facilities 202 such as the consumer 200-2 whose power generation output is not measured.

図5の範囲301は、需要家200-1の太陽光発電設備202であるPV1を中心とした上記予め定められた距離を半径とする円であり、範囲301内の需要家の太陽光発電設備202の発電出力の推定に、PV1の発電出力を用いることができる。同様に、図5の範囲302は、需要家200-3の太陽光発電設備202であるPV3を中心とした上記予め定められた距離を半径とする円であり、範囲302内の需要家の太陽光発電設備202の発電出力の推定に、PV3の発電出力を用いることができる。したがって、例えば、需要家200-2の太陽光発電設備202であるPV2が発電出力の推定対象として設定されている場合、PV2は範囲301内でありかつ範囲302内であるため、PV1に対応する需要家200-1とPV3に対応する需要家200-3との両方が第1需要家の候補となる。このように、複数の第1需要家の候補が存在する場合には、推定対象のPV2に近い方の需要家である需要家200-2が第1需要家として選択される。あるいは、複数の第1需要家の候補の計測データを合算して使うことにしてもよく、その場合は、複数の需要家が第1需要家として選択される。 The range 301 of FIG. 5 is a circle whose radius is the above-mentioned predetermined distance centered on PV1 which is the solar power generation facility 202 of the consumer 200-1, and the solar power generation facility of the consumer within the range 301. The power generation output of PV1 can be used to estimate the power generation output of 202. Similarly, the range 302 of FIG. 5 is a circle whose radius is the above-mentioned predetermined distance centered on PV3, which is the photovoltaic power generation facility 202 of the consumer 200-3, and the sun of the consumer within the range 302. The power generation output of PV3 can be used to estimate the power generation output of the photovoltaic power generation facility 202. Therefore, for example, when PV2, which is the photovoltaic power generation facility 202 of the consumer 200-2, is set as the estimation target of the power generation output, PV2 corresponds to PV1 because it is within the range 301 and within the range 302. Both the consumer 200-1 and the consumer 200-3 corresponding to PV3 are candidates for the first consumer. In this way, when there are a plurality of candidates for the first consumer, the consumer 200-2, which is the consumer closer to the PV2 to be estimated, is selected as the first consumer. Alternatively, the measurement data of a plurality of candidates for the first consumer may be combined and used, in which case the plurality of consumers are selected as the first consumer.

図4の説明に戻る。ステップS1の後、第1共分散算出部120は、第1期間を設定する(ステップS2)。詳細には、第1共分散算出部120は、記憶部180に格納されている太陽光発電出力データ182内の第1需要家の太陽光発電設備202の発電出力である第1太陽光発電出力の時系列データを参照し、日射強度が強くかつ日射強度の変動が大きい期間を探索する。例えば、第1共分散算出部120は、第2太陽光発電出力を推定する日の前日のうち、第1太陽光発電出力が大きくかつ第1太陽光発電出力の変動が大きい6時間から10時間程度の期間を探索する。そして、第1共分散算出部120は、探索した期間を第1期間として設定する。なお、第1共分散算出部120は、上記の期間(例えば、第2太陽光発電出力を推定する日の前日)の中から、第1期間を探索するのではなく、上記の期間自体を探索してから、当該期間の中から第1期間を探索することにしてもよい。なお、上述したとおり、この第1太陽光発電出力が大きくかつ第1太陽光発電出力の変動が大きい6時間から10時間程度の期間の探索を行わずに任意の方法で第1期間を設定してもよい。 Returning to the description of FIG. After step S1, the first covariance calculation unit 120 sets the first period (step S2). Specifically, the first covariance calculation unit 120 is the first photovoltaic power generation output which is the power generation output of the photovoltaic power generation facility 202 of the first consumer in the photovoltaic power generation output data 182 stored in the storage unit 180. Refer to the time series data of, and search for the period when the solar intensity is strong and the fluctuation of the solar intensity is large. For example, the first co-dispersion calculation unit 120 has 6 to 10 hours on the day before the day when the second photovoltaic power generation output is estimated, in which the first photovoltaic power generation output is large and the fluctuation of the first photovoltaic power generation output is large. Explore a period of time. Then, the first covariance calculation unit 120 sets the searched period as the first period. The first covariance calculation unit 120 does not search for the first period from the above period (for example, the day before the day when the second photovoltaic power generation output is estimated), but searches for the above period itself. Then, the first period may be searched from the period. As described above, the first period is set by an arbitrary method without searching for a period of about 6 to 10 hours in which the first photovoltaic power generation output is large and the fluctuation of the first photovoltaic power generation output is large. You may.

次に、第1共分散算出部120は、自己係数λを算出するための第1共分散を算出する(ステップS3)。第1共分散算出部120は、算出した第1共分散を記憶部180に算出データ185として格納する。第1共分散の算出の詳細については後述する。 Next, the first covariance calculation unit 120 calculates the first covariance for calculating the self-coefficient λ (step S3). The first covariance calculation unit 120 stores the calculated first covariance in the storage unit 180 as calculated data 185. The details of the calculation of the first covariance will be described later.

次に、第2共分散算出部130は、自己係数λを算出するための第2共分散を算出する(ステップS4)。第2共分散算出部130は、算出した第2共分散を記憶部180に算出データ185として格納する。第2共分散の算出の詳細については後述する。 Next, the second covariance calculation unit 130 calculates the second covariance for calculating the self-coefficient λ (step S4). The second covariance calculation unit 130 stores the calculated second covariance in the storage unit 180 as calculated data 185. The details of the calculation of the second covariance will be described later.

次に、推定精度検証部140は、自己係数λを算出する(ステップS5)。具体的には、推定精度検証部140は、記憶部180に格納されている自己係数λを算出するための第1共分散および第2共分散を用いて、上述した式(7)により自己係数λを算出する。 Next, the estimation accuracy verification unit 140 calculates the self-coefficient λ (step S5). Specifically, the estimation accuracy verification unit 140 uses the first covariance and the second covariance for calculating the self-coefficient λ stored in the storage unit 180, and the self-coefficient according to the above equation (7). Calculate λ.

次に、推定精度検証部140は、|λ-1|が閾値u以下であるか否かを判定する(ステップS6)。|λ-1|が閾値u以下である場合(ステップS6 Yes)、係数αの算出が行われ(ステップS7)、係数αの推定を終了するか否かの判定が行われる(ステップS8)。係数αの算出処理の詳細については後述する。ステップS8の処理は第1共分散算出部120または推定部150により行われる。ステップS8の係数αの推定を終了するか否かの判定は、係数αの推定対象の全ての太陽光発電設備202に関してステップS2~ステップS6の処理が行われたか否かの判定である。係数αの推定対象の全ての太陽光発電設備202に関してステップS2~ステップS6の処理が行われた場合、ステップS8では、Yesと判定され、係数αの推定処理は終了する。係数αの推定対象の太陽光発電設備202のうちステップS2~ステップS6の処理が行われていない太陽光発電設備202がある場合、ステップS8ではNoと判定され、ステップS1からの処理が繰り返される。 Next, the estimation accuracy verification unit 140 determines whether or not | λ-1 | is equal to or less than the threshold value u (step S6). When | λ-1 | is equal to or less than the threshold value u (step S6 Yes), the coefficient α is calculated (step S7), and it is determined whether or not to end the estimation of the coefficient α (step S8). The details of the calculation process of the coefficient α will be described later. The process of step S8 is performed by the first covariance calculation unit 120 or the estimation unit 150. The determination of whether or not to end the estimation of the coefficient α in step S8 is the determination of whether or not the processes of steps S2 to S6 have been performed for all the photovoltaic power generation facilities 202 to be estimated by the coefficient α. When the processes of steps S2 to S6 are performed for all the photovoltaic power generation facilities 202 for which the coefficient α is estimated, it is determined as Yes in step S8, and the coefficient α estimation process ends. If there is a photovoltaic power generation facility 202 that has not been processed in steps S2 to S6 among the photovoltaic power generation facilities 202 for which the coefficient α is estimated, it is determined as No in step S8, and the processing from step S1 is repeated. ..

一方、|λ-1|が閾値uを超える場合(ステップS6 No)、ステップS7は実施されず、処理はステップS8へ進む。上述したように、自己係数λは、第1需要家の残余需要を用いて、係数αと同様の手法で算出された値であり、理想的には、第1太陽光発電出力にλを乗じると第1太陽光発電出力に一致する。したがって、自己係数λが1から離れるにつれて自己係数λの誤差が大きくなることになる。自己係数λは係数αと同様の手法により算出されているため、自己係数λの誤差が大きいときには係数αの誤差も大きいと想定でき、自己係数λの算出に用いた第1期間は、係数αの算出に適していない期間であると推定される。このため、本実施の形態では、|λ-1|が閾値uを超える場合に、係数αの算出を行わないことで係数αの推定精度の低下を抑制することができる。ここでは、|λ-1|が閾値uを超える場合に、係数αの算出を行わないようにしたが、ステップS6の前に係数αの算出は行った上で、ステップS6でNoと判定された場合に、係数αの更新を行わないようにしてもよい。 On the other hand, when | λ-1 | exceeds the threshold value u (step S6 No), step S7 is not executed and the process proceeds to step S8. As described above, the self-coefficient λ is a value calculated by the same method as the coefficient α using the residual demand of the first consumer, and ideally, the first photovoltaic power generation output is multiplied by λ. And the first photovoltaic power generation output. Therefore, as the self-coefficient λ deviates from 1, the error of the self-coefficient λ increases. Since the self-coefficient λ is calculated by the same method as the coefficient α, it can be assumed that the error of the coefficient α is also large when the error of the self-coefficient λ is large. It is estimated that the period is not suitable for the calculation of. Therefore, in the present embodiment, when | λ-1 | exceeds the threshold value u, it is possible to suppress a decrease in the estimation accuracy of the coefficient α by not calculating the coefficient α. Here, when | λ-1 | exceeds the threshold value u, the coefficient α is not calculated, but after the coefficient α is calculated before step S6, it is determined as No in step S6. In that case, the coefficient α may not be updated.

なお、推定精度検証部140は、係数αの推定を行わない場合、すなわちステップS6でNoの場合、推定対象の第2需要家の推定がその時点までに一度も行われていないときには、記憶部180の算出データ185内の係数αとして、係数初期値データ184として格納されている値を格納する。なお、係数初期値データ184は、第2需要家ごとの係数αの初期値である。係数αの初期値は、第2需要家ごとに対応する第1需要家が予め定められている場合には、第2需要家の太陽光発電設備202と第1需要家の太陽光発電設備202との定格容量の比であってもよい。第2需要家ごとに対応する第1需要家が予め定められていない場合、第1需要家の候補となる全ての需要家ごとに、第2需要家の太陽光発電設備202と第1需要家の候補の太陽光発電設備202との定格容量の比が初期値として格納されていてもよい。 In addition, when the estimation accuracy verification unit 140 does not estimate the coefficient α, that is, when No in step S6, when the estimation of the second consumer to be estimated has not been performed even once by that time, the storage unit 140 As the coefficient α in the calculated data 185 of 180, the value stored as the coefficient initial value data 184 is stored. The coefficient initial value data 184 is an initial value of the coefficient α for each second consumer. The initial value of the coefficient α is the solar power generation facility 202 of the second consumer and the photovoltaic power generation facility 202 of the first consumer when the first consumer corresponding to each second consumer is predetermined. It may be the ratio of the rated capacity to. If the first consumer corresponding to each second consumer is not predetermined, the second consumer's photovoltaic power generation facility 202 and the first consumer are used for all the customers who are candidates for the first consumer. The ratio of the rated capacity to the candidate photovoltaic power generation facility 202 may be stored as an initial value.

図6は、本実施の形態の記憶部180に格納される算出データ185内のデータの一例を示す図である。図6に示すように、記憶部180には、算出データ185として、第1共分散、第2共分散、係数α、遅延時間、発電出力推定値が格納される。算出データ185内の第1共分散は、ステップS3で算出された第1共分散と、ステップS7の係数αの推定の過程で算出される第1共分散とを含む。算出データ185内の第2共分散は、ステップS4で算出された第2共分散と、ステップS7の係数αの推定の過程で算出される第2共分散とを含む。第1共分散、第2共分散については、一時的に記憶されていればよいため、発電出力推定対象の各需要家に対応する係数αの推定が終了するたびに消去されてもよい。係数αは、ステップS7の係数αの推定処理が行われると、推定値で更新される。係数αおよび遅延時間は、第2需要家と第1需要家との組み合わせごと、すなわち第2太陽光発電設備と第1太陽光発電設備との組み合わせごとに異なるため、この組み合わせごとに算出データ185として記憶部180に格納される。 FIG. 6 is a diagram showing an example of data in the calculated data 185 stored in the storage unit 180 of the present embodiment. As shown in FIG. 6, the storage unit 180 stores the first covariance, the second covariance, the coefficient α, the delay time, and the estimated power generation output as the calculated data 185. The first covariance in the calculated data 185 includes the first covariance calculated in step S3 and the first covariance calculated in the process of estimating the coefficient α in step S7. The second covariance in the calculated data 185 includes the second covariance calculated in step S4 and the second covariance calculated in the process of estimating the coefficient α in step S7. Since the first covariance and the second covariance need only be temporarily stored, they may be deleted each time the estimation of the coefficient α corresponding to each consumer whose power generation output is estimated is completed. The coefficient α is updated with the estimated value when the estimation process of the coefficient α in step S7 is performed. Since the coefficient α and the delay time are different for each combination of the second consumer and the first consumer, that is, for each combination of the second photovoltaic power generation facility and the first photovoltaic power generation facility, the calculated data 185 for each combination. Is stored in the storage unit 180.

図7は、本実施の形態の記憶部180に格納される、太陽光発電設備202の組み合わせごとの係数αの一例を示す図である。図7に示すように、ステップS7の係数αの推定が行われると、発電出力推定対象の太陽光発電設備202である第2太陽光発電設備と、推定に使用する太陽光発電設備202である第1太陽光発電設備との組み合わせごとに記憶部180の算出データ185として係数αの推定値が格納される。例えば、需要家200-2の太陽光発電設備202であるPV2は、需要家200-1の太陽光発電設備202であるPV1の発電出力を用いて発電出力の推定が行われ、推定に用いられる係数αの推定値はαである。上述したように、係数αの推定が一度も行われておらず、上述したステップS6でNoと判定された場合には、算出データ185の係数αには、係数初期値データ184内の対応する初期値が格納される。 FIG. 7 is a diagram showing an example of the coefficient α for each combination of the photovoltaic power generation equipment 202 stored in the storage unit 180 of the present embodiment. As shown in FIG. 7, when the coefficient α in step S7 is estimated, the second photovoltaic power generation facility 202, which is the photovoltaic power generation facility 202 for which the power generation output is estimated, and the photovoltaic power generation facility 202 used for the estimation. The estimated value of the coefficient α is stored as the calculated data 185 of the storage unit 180 for each combination with the first photovoltaic power generation facility. For example, PV2, which is the photovoltaic power generation facility 202 of the consumer 200-2, is estimated to have a power generation output using the power generation output of PV1, which is the photovoltaic power generation facility 202 of the consumer 200-1, and is used for the estimation. The estimated value of the coefficient α is α 1 . As described above, when the coefficient α has never been estimated and is determined to be No in step S6 described above, the coefficient α of the calculated data 185 corresponds to the coefficient initial value data 184. The initial value is stored.

図4のステップS3の処理である第1共分散の算出処理について説明する。図8は、本実施の形態の第1共分散算出部120による第1共分散の算出処理手順の一例を示すフローチャートである。なお、ステップS3およびステップS4は、自己係数λを算出するための処理であるが、上述した式(5)に示したように係数αの算出にも第1共分散および第2共分散が用いられる。図8では、自己係数λの算出のための第1共分散を算出する処理と係数αの算出のための第1共分散を算出する処理とで共用されるフローチャートを示す。同様に、後述する図9においても、自己係数λの算出のための第2共分散を算出する処理と係数αの算出のための第2共分散を算出する処理とで共用されるフローチャートを示す。なお、上述した式(5)および式(7)に示したように、係数αの算出には遅延時間が用いられるのに対し、自己係数λの算出には遅延時間が用いられない。このため、図8および図9に示す処理では、遅延時間を考慮するか否かの判断が行われる。 The calculation process of the first covariance, which is the process of step S3 of FIG. 4, will be described. FIG. 8 is a flowchart showing an example of the calculation processing procedure of the first covariance by the first covariance calculation unit 120 of the present embodiment. Although steps S3 and S4 are processes for calculating the self-coefficient λ, the first covariance and the second covariance are also used for calculating the coefficient α as shown in the above equation (5). Be done. FIG. 8 shows a flowchart shared by the process of calculating the first covariance for calculating the self-coefficient λ and the process of calculating the first covariance for calculating the coefficient α. Similarly, FIG. 9, which will be described later, also shows a flowchart shared by the process of calculating the second covariance for calculating the self-coefficient λ and the process of calculating the second covariance for calculating the coefficient α. .. As shown in the above equations (5) and (7), the delay time is used to calculate the coefficient α, whereas the delay time is not used to calculate the self-coefficient λ. Therefore, in the processes shown in FIGS. 8 and 9, it is determined whether or not to consider the delay time.

図8に示すように、まず、第1共分散算出部120は、第1期間での第1太陽光発電出力の時系列データと残余需要の時系列データとを取得する(ステップS11)。具体的には、第1共分散算出部120は、記憶部180の太陽光発電出力データ182のうち第1期間における第1太陽光発電出力を読み出すことで、第1期間での第1太陽光発電出力の時系列データを取得する。また、第1共分散算出部120は、自己係数λの算出のための第1共分散を算出するときには、記憶部180の残余需要データ183のうち第1期間における第1需要家の残余需要を読み出すことで、第1期間での残余需要の時系列データを取得する。また、係数αの算出のための第1共分散を算出するときには、記憶部180の残余需要データ183のうち第1期間における第2需要家の残余需要を読み出すことで、第1期間での残余需要の時系列データを取得する。なお、第1期間での第1太陽光発電出力の時系列データは、自己係数λの算出のための処理と係数αの算出のため処理との両方で用いられるため、第1共分散算出部120は、自己係数λの算出のために読み出した第1期間での第1太陽光発電出力の時系列データを、係数αの算出のため処理が行われるまで保持することで、係数αの算出のためのステップS11において、第1太陽光発電電力を記憶部180の太陽光発電出力データ182内から抽出して読み出す処理を省略してもよい。 As shown in FIG. 8, first, the first covariance calculation unit 120 acquires the time-series data of the first photovoltaic power generation output in the first period and the time-series data of the residual demand (step S11). Specifically, the first co-dispersion calculation unit 120 reads out the first photovoltaic power generation output in the first period from the photovoltaic power generation output data 182 of the storage unit 180, so that the first solar power in the first period can be read. Acquire time-series data of power generation output. Further, when the first covariance calculation unit 120 calculates the first covariance for calculating the self-covariance λ, the first covariance calculation unit 120 determines the residual demand of the first consumer in the first period of the residual demand data 183 of the storage unit 180. By reading, the time series data of the residual demand in the first period is acquired. Further, when calculating the first covariance for calculating the coefficient α, the residual demand in the first period is read out from the residual demand data 183 of the storage unit 180 by reading the residual demand of the second consumer in the first period. Get time series data of demand. Since the time-series data of the first photovoltaic power generation output in the first period is used in both the processing for calculating the self-factor λ and the processing for calculating the coefficient α, the first covariance calculation unit. 120 calculates the coefficient α by holding the time-series data of the first photovoltaic power generation output in the first period read for the calculation of the self-factor λ until the processing is performed for the calculation of the coefficient α. In step S11 for the above, the process of extracting and reading the first photovoltaic power generation power from the photovoltaic power generation output data 182 of the storage unit 180 may be omitted.

なお、太陽光発電出力データ182は、上述したようにスマートメータ204および計測器206により計測された計測データであるが、これらの計測データの代わりに推定値が用いられてもよい。例えば、太陽光発電出力推定装置100が、本実施の形態の手法または他の手法を用いて太陽光発電出力を推定し、推定した太陽光発電出力を太陽光発電出力データ182として記憶しておいてもよい。または、ユーザから、太陽光発電出力の推定値または実測値が入力され、この入力された推定値または実測値が太陽光発電出力データ182として記憶されてもよい。 The photovoltaic power generation output data 182 is the measurement data measured by the smart meter 204 and the measuring instrument 206 as described above, but an estimated value may be used instead of these measurement data. For example, the photovoltaic power generation output estimation device 100 estimates the photovoltaic power generation output by using the method of the present embodiment or another method, and stores the estimated photovoltaic power generation output as the photovoltaic power generation output data 182. You may. Alternatively, the user may input an estimated value or an actually measured value of the photovoltaic power generation output, and the input estimated value or the actually measured value may be stored as the photovoltaic power generation output data 182.

第1共分散算出部120は、第1太陽光発電出力の時系列データと残余需要の時系列データとの共分散を第1共分散として算出する(ステップS12)。第1共分散算出部120は、算出した第1共分散を、記憶部180に算出データ185として書き込む。 The first covariance calculation unit 120 calculates the covariance between the time-series data of the first photovoltaic power generation output and the time-series data of the residual demand as the first covariance (step S12). The first covariance calculation unit 120 writes the calculated first covariance in the storage unit 180 as calculated data 185.

さらに、第1共分散算出部120は、遅延時間を算出するか否かを判断し(ステップS13)、遅延時間を算出すると判断した場合(ステップS13 Yes)、読み出し済の第1太陽光発電出力と、同じく読み出し済の第2需要家の残余需要とを用いて、第1期間における遅延時間τを算出し(ステップS14)、記憶部180に算出データ185として格納し、第1共分散の算出処理を、終了する。ステップS13では、第1共分散算出部120は、自己係数λの算出のための処理すなわち図4に示したステップS3の処理である場合は遅延時間を算出しないと判断し、係数αの算出のための処理すなわち図4に示したステップS7における第1共分散の算出処理である場合は遅延時間を算出すると判断する。なお、係数αの算出のための処理である場合も、遅延時間は算出せずに、0と近似することにしてもよい。 Further, the first covariance calculation unit 120 determines whether or not to calculate the delay time (step S13), and when it is determined to calculate the delay time (step S13 Yes), the read-out first solar power generation output. And the residual demand of the second consumer who has also been read out, the delay time τ s in the first period is calculated (step S14), stored in the storage unit 180 as the calculated data 185, and the first covariance is obtained. The calculation process is terminated. In step S13, the first covariance calculation unit 120 determines that the delay time is not calculated in the case of the process for calculating the self-covariance λ, that is, the process in step S3 shown in FIG. 4, and calculates the coefficient α. In the case of the process for calculating the first covariance in step S7 shown in FIG. 4, it is determined that the delay time is calculated. Even in the case of processing for calculating the coefficient α, the delay time may not be calculated and may be approximated to 0.

ステップS14では、詳細には、第1共分散算出部120は、発電出力の推定対象である第2太陽光発電設備の設置地点と、第1太陽光発電設備の設置地点との間の遅延時間を算出し、記憶部180に算出データ185として格納する。この遅延時間は、上述したように第1太陽光発電設備の設置地点と第2太陽光発電設備の設置地点との間を日射変動が伝播する時間である。第1共分散算出部120は、例えば、第1太陽光発電出力と第2需要家の残余需要との時系列データの共分散関数Cov[PPV1(t),P(t+τ)]が最小値(符号が「-」かつ絶対値が最大値)をとるときのタイムラグに-1を乗じて遅延時間τを算出する。 In step S14, in detail, the first covariance calculation unit 120 has a delay time between the installation point of the second photovoltaic power generation facility, which is the estimation target of the power generation output, and the installation point of the first photovoltaic power generation facility. Is calculated and stored in the storage unit 180 as the calculated data 185. This delay time is the time for the solar radiation fluctuation to propagate between the installation point of the first photovoltaic power generation facility and the installation point of the second photovoltaic power generation facility as described above. The first covariance calculation unit 120 has, for example, a covariance function Cov t [P PV1 (t), P2 ( t + τ)] of time-series data of the first solar power output and the residual demand of the second consumer. The delay time τ s is calculated by multiplying the time lag when taking the minimum value (the sign is “−” and the absolute value is the maximum value) by -1.

ステップS13において遅延時間を算出しないと判断した場合(ステップS13 No)、第1共分散算出部120は、第1共分散の算出処理を、終了する。 When it is determined in step S13 that the delay time is not calculated (step S13 No), the first covariance calculation unit 120 ends the calculation process of the first covariance.

図4のステップS4の処理である第2共分散の算出処理について説明する。図9は、本実施の形態の第2共分散算出部130による第2共分散の算出処理手順の一例を示すフローチャートである。 The second covariance calculation process, which is the process of step S4 of FIG. 4, will be described. FIG. 9 is a flowchart showing an example of the calculation processing procedure of the second covariance by the second covariance calculation unit 130 of the present embodiment.

第2共分散算出部130は、遅延時間を使用するか否かを判断し(ステップS21)、遅延時間を使用しないと判断した場合(ステップS21 No)、遅延時間τを0とし、ステップS22を実施せずに、後述するステップS23以降の処理を実施する。遅延時間を使用すると判断した場合(ステップS21 Yes)、第2共分散算出部130は、記憶部180に格納されている算出データ185内の遅延時間τを読み出すことで遅延時間τを取得する(ステップS22)。なお、ステップS21では、第2共分散算出部130は、自己係数λの算出のための処理すなわち図4に示したステップS4の処理である場合は遅延時間を使用しないと判断し、係数αの算出のための処理すなわち図4に示したステップS7における第2共分散の算出処理である場合は遅延時間を使用すると判断する。 The second covariance calculation unit 130 determines whether or not to use the delay time (step S21), and if it determines that the delay time is not used (step S21 No), sets the delay time τ s to 0 and steps S22. Is not performed, but the processes after step S23, which will be described later, are carried out. When it is determined to use the delay time (step S21 Yes), the second covariance calculation unit 130 acquires the delay time τ s by reading the delay time τ s in the calculated data 185 stored in the storage unit 180. (Step S22). In step S21, the second covariance calculation unit 130 determines that the delay time is not used in the case of the process for calculating the self-covariance λ, that is, the process in step S4 shown in FIG. 4, and determines that the delay time is not used, and the coefficient α In the case of the processing for calculation, that is, the calculation processing of the second covariance in step S7 shown in FIG. 4, it is determined that the delay time is used.

次に、第2共分散算出部130は、遅延時間τずれた2つの第1太陽光発電出力の時系列データを取得する(ステップS23)。つまり、第2共分散算出部130は、第1期間における第1太陽光発電出力の時系列データと、第1期間より遅延時間τずれた第2期間における第1太陽光発電出力の時系列データとを、記憶部180の太陽光発電出力データ182から抽出して読み出して取得する。遅延時間τが0の場合は、遅延時間τずれた2つの第1太陽光発電出力は同一となるため、第2共分散算出部130は、記憶部180の太陽光発電出力データ182から第1期間の第1太陽光発電出力を読み出せばよい。なお、第1期間での第1太陽光発電出力の時系列データは、自己係数λの算出のための処理と係数αの算出のため処理との両方で用いられるため、第2共分散算出部130は、自己係数λの算出のために読み出した第1期間での第1太陽光発電出力の時系列データを、係数αの算出のため処理が行われるまで保持することで、係数αの算出のためのステップS23において、第1期間における第1太陽光発電電力を記憶部180の太陽光発電出力データ182内から抽出して読み出す処理を省略してもよい。 Next, the second covariance calculation unit 130 acquires time-series data of the two first photovoltaic power generation outputs deviated by the delay time τ s (step S23). That is, the second covariance calculation unit 130 uses the time series data of the first photovoltaic power generation output in the first period and the time series of the first photovoltaic power generation output in the second period, which is delayed by τ s from the first period. The data is extracted from the photovoltaic power generation output data 182 of the storage unit 180, read out, and acquired. When the delay time τ s is 0, the two first photovoltaic power generation outputs deviated by the delay time τ s are the same, so that the second covariance calculation unit 130 is based on the photovoltaic power generation output data 182 of the storage unit 180. The first photovoltaic power generation output of the first period may be read out. Since the time-series data of the first photovoltaic power generation output in the first period is used in both the processing for calculating the self-factor λ and the processing for calculating the coefficient α, the second covariance calculation unit. The 130 calculates the coefficient α by holding the time-series data of the first photovoltaic power generation output in the first period read for the calculation of the self-factor λ until the processing is performed for the calculation of the coefficient α. In step S23 for the above, the process of extracting and reading the first photovoltaic power generation power in the first period from the photovoltaic power generation output data 182 of the storage unit 180 may be omitted.

次に、第2共分散算出部130は、取得したデータを用いて、第2共分散を算出する(ステップS24)。詳細には、第2共分散算出部130は、遅延時間τずれた2つの第1太陽光発電出力の時系列データの自己共分散である第2共分散を算出する。具体的には、第2共分散算出部130は、第1期間における第1太陽光発電出力の時系列データと、第2期間における第1太陽光発電出力の時系列データとの自己共分散を第2共分散として算出する。なお、遅延時間τが0の場合は、第2共分散算出部130は、第1期間における第1太陽光発電出力の時系列データの分散を第2共分散として算出する。第2共分散算出部130は、算出した第2共分散を算出データ185に書き込み、第2共分散の算出処理を終了する。 Next, the second covariance calculation unit 130 calculates the second covariance using the acquired data (step S24). Specifically, the second covariance calculation unit 130 calculates the second covariance, which is the autocovariance of the time-series data of the two first photovoltaic power generation outputs deviated by the delay time τ s . Specifically, the second covariance calculation unit 130 performs autocovariance between the time-series data of the first photovoltaic power generation output in the first period and the time-series data of the first photovoltaic power generation output in the second period. Calculated as the second covariance. When the delay time τ s is 0, the second covariance calculation unit 130 calculates the variance of the time series data of the first photovoltaic power generation output in the first period as the second covariance. The second covariance calculation unit 130 writes the calculated second covariance in the calculation data 185, and ends the calculation process of the second covariance.

図4のステップS7の処理である係数αの算出処理について説明する。図10は、本実施の形態の係数αの算出処理手順の一例を示すフローチャートである。まず、第1共分散算出部120は、係数αを算出するための第1共分散を算出する(ステップS31)。詳細には、推定精度検証部140が、図4のステップS6でYesと判定すると、第1共分散算出部120へ、係数αの算出のための第1共分散の算出を指示し、第1共分散算出部120が、図8に示した処理により、第1共分散および遅延時間τを算出し、記憶部180に算出データ185として格納する。なお、係数αの算出のための第1共分散を算出するときには、第1共分散算出部120は、図8に示したステップS11において、記憶部180の残余需要データ183のうち第1期間における第2需要家の残余需要を読み出すことで、第1期間での残余需要の時系列データを取得する。これにより、第1太陽光発電出力の時系列データと第2需要家の残余需要の時系列データとの共分散が第1共分散として算出される。なお、上述したように、遅延時間τを0と近似し遅延時間τの算出を省いてもよい。 The calculation process of the coefficient α, which is the process of step S7 in FIG. 4, will be described. FIG. 10 is a flowchart showing an example of the calculation processing procedure of the coefficient α of the present embodiment. First, the first covariance calculation unit 120 calculates the first covariance for calculating the coefficient α (step S31). Specifically, when the estimation accuracy verification unit 140 determines Yes in step S6 of FIG. 4, it instructs the first covariance calculation unit 120 to calculate the first covariance for calculating the coefficient α, and first. The covariance calculation unit 120 calculates the first covariance and the delay time τ s by the process shown in FIG. 8, and stores the first covariance and the delay time τ s in the storage unit 180 as the calculated data 185. When calculating the first covariance for calculating the coefficient α, the first covariance calculation unit 120 in step S11 shown in FIG. 8 in the first period of the residual demand data 183 of the storage unit 180. By reading the residual demand of the second consumer, the time series data of the residual demand in the first period is acquired. As a result, the covariance between the time-series data of the first photovoltaic power generation output and the time-series data of the residual demand of the second consumer is calculated as the first covariance. As described above, the delay time τ s may be approximated to 0 and the calculation of the delay time τ s may be omitted.

次に、第2共分散算出部130は、係数αを算出するための第2共分散を算出する(ステップS32)。詳細には、第2共分散算出部130は、図9に示した手順により、第1期間における第1太陽光発電出力の時系列データと、第1期間より遅延時間τずれた第2期間における第1太陽光発電出力の時系列データとの自己共分散を第2共分散として算出し、記憶部180に算出データ185として格納する。 Next, the second covariance calculation unit 130 calculates the second covariance for calculating the coefficient α (step S32). Specifically, the second covariance calculation unit 130 uses the procedure shown in FIG. 9 to obtain the time-series data of the first solar power generation output in the first period and the second period in which the delay time τ s deviates from the first period. The autocovariance with the time-series data of the first solar power generation output in the above is calculated as the second covariance, and is stored in the storage unit 180 as the calculated data 185.

次に、推定部150は、係数αを算出して更新する(ステップS33)。詳細には、推定部150は、記憶部180の算出データ185内の係数αの算出のための第1共分散および第2共分散を読み出し、読み出した第1共分散を第2共分散で除して-1を乗じることで、係数αを算出する。そして、推定部150は、記憶部180の算出データ185内の係数αを、算出した値で更新する。以上により、係数αの推定処理が終了する。 Next, the estimation unit 150 calculates and updates the coefficient α (step S33). Specifically, the estimation unit 150 reads out the first covariance and the second covariance for calculating the coefficient α in the calculated data 185 of the storage unit 180, and divides the read first covariance by the second covariance. Then multiply by -1 to calculate the coefficient α. Then, the estimation unit 150 updates the coefficient α in the calculated data 185 of the storage unit 180 with the calculated value. As a result, the estimation process of the coefficient α is completed.

図4および図8~図10を用いて説明した処理により、係数αの推定値が算出される。また、図4に示したように、自己係数λを用いて係数αの推定に適していないと判断された場合には、係数αの算出は行われない。なお、以上述べた例では、自己係数λを用いて係数αの推定に適していないと判断された場合には、係数αの算出を行わないようにしたが、自己係数λの算出とともに、推定部150が係数αについても算出して一時的に保持しておき、自己係数λを用いて係数αの推定に適していると判断された場合に、記憶部180の算出データ185内の係数αを算出した値で更新するようにしてもよい。この場合、推定精度検証部140は、自己係数λを用いて係数αの推定に適していないと判断すると係数αの更新を行わないように推定部150に指示し、推定部150はこの指示を受けると記憶部180の算出データ185内の係数αの更新は行わない。 The estimated value of the coefficient α is calculated by the process described with reference to FIGS. 4 and 8 to 10. Further, as shown in FIG. 4, when it is determined that the self-coefficient λ is not suitable for estimating the coefficient α, the coefficient α is not calculated. In the above-mentioned example, when it is judged that the self-coefficient λ is not suitable for the estimation of the coefficient α, the calculation of the coefficient α is not performed, but it is estimated together with the calculation of the self-coefficient λ. When the unit 150 also calculates and temporarily holds the coefficient α and determines that it is suitable for estimating the coefficient α using the self-coefficient λ, the coefficient α in the calculated data 185 of the storage unit 180 is determined. May be updated with the calculated value. In this case, the estimation accuracy verification unit 140 instructs the estimation unit 150 not to update the coefficient α when it is determined that the coefficient α is not suitable for estimation using the self-coefficient λ, and the estimation unit 150 gives this instruction. Upon receiving it, the coefficient α in the calculated data 185 of the storage unit 180 is not updated.

次に、係数αを用いた太陽光発電出力の推定処理について説明する。図11は、本実施の形態の太陽光発電出力の推定処理手順の一例を示すフローチャートである。推定部150は、推定対象の太陽光発電設備である第2太陽光発電設備を決めると、推定時点における第2太陽光発電設備の発電出力の推定のために図11に示す処理を実施する。推定時点は、上述した第1期間内であってもよいし、第1期間と異なっていてもよい。なお、第2太陽光発電設備に対応する第1太陽光発電設備は、係数αの推定処理で選択されており、第2太陽光発電設備および第1太陽光発電設備の組み合わせごとに係数αが算出されているので、太陽光発電出力の推定処理においても、第2太陽光発電設備および第1太陽光発電設備の組み合わせは係数αの算出において用いられた組み合わせを用いる。例えば、図7に示すような組み合わせで係数αが算出されている場合、需要家200-2の太陽光発電設備202の発電出力の推定では、需要家200-1の太陽光発電設備202が第1太陽光発電設備として用いられる。 Next, the estimation process of the photovoltaic power generation output using the coefficient α will be described. FIG. 11 is a flowchart showing an example of the estimation processing procedure of the photovoltaic power generation output of the present embodiment. When the estimation unit 150 determines the second photovoltaic power generation facility, which is the photovoltaic power generation facility to be estimated, the estimation unit 150 performs the process shown in FIG. 11 for estimating the power generation output of the second photovoltaic power generation facility at the time of estimation. The estimated time point may be within the first period described above, or may be different from the first period. The first photovoltaic power generation facility corresponding to the second photovoltaic power generation facility is selected by the estimation process of the coefficient α, and the coefficient α is set for each combination of the second photovoltaic power generation facility and the first photovoltaic power generation facility. Since it has been calculated, the combination used in the calculation of the coefficient α is used as the combination of the second photovoltaic power generation facility and the first photovoltaic power generation facility also in the estimation process of the photovoltaic power generation output. For example, when the coefficient α is calculated by the combination as shown in FIG. 7, in the estimation of the power generation output of the photovoltaic power generation facility 202 of the consumer 200-2, the photovoltaic power generation facility 202 of the consumer 200-1 is the first. 1 Used as a solar power generation facility.

図11に示すように、推定部150は、まず、遅延時間を算出する(ステップS41)。詳細には、推定部150は、推定時点における第1太陽光発電設備の設置地点と第2太陽光発電設備の設置地点との遅延時間を算出する。推定部150は、記憶部180に算出データ185として格納されている遅延時間を推定時点における遅延時間として利用してもよいし、推定時点における気象データなどを反映して遅延時間を更新することにしてもよい。例えば、推定部150は、気象データに含まれる風向風速データから雲が通過する時間を算出して算出した値を用いて遅延時間を更新してもよいし、衛星画像等を用いて遅延時間を更新してもよいし、これら以外の方法により気象データを用いた更新を行ってもよい。また、遅延時間を0と近似することにしてもよい。 As shown in FIG. 11, the estimation unit 150 first calculates the delay time (step S41). Specifically, the estimation unit 150 calculates the delay time between the installation point of the first photovoltaic power generation facility and the installation point of the second photovoltaic power generation facility at the time of estimation. The estimation unit 150 may use the delay time stored as the calculated data 185 in the storage unit 180 as the delay time at the estimation time, or update the delay time by reflecting the weather data at the estimation time. You may. For example, the estimation unit 150 may update the delay time using a value calculated by calculating the time for the clouds to pass from the wind direction and wind velocity data included in the meteorological data, or may update the delay time using a satellite image or the like. It may be updated, or it may be updated using meteorological data by a method other than these. Further, the delay time may be approximated to 0.

次に、推定部150は、推定時点から遅延時間ずれた時刻の太陽光発電出力計測データを取得する(ステップS42)。詳細には、推定部150は、記憶部180の太陽光発電出力データ182内の第1太陽光発電設備に対応するデータのうち推定時点から遅延時間ずれた時刻のデータを取得する。なお、ここでは、太陽光発電出力データ182が計測データであるとして説明するが、上述したように、計測データのかわりに推定値が用いられてもよい。 Next, the estimation unit 150 acquires the photovoltaic power generation output measurement data at a time deviated from the estimation time by the delay time (step S42). Specifically, the estimation unit 150 acquires the data corresponding to the first photovoltaic power generation facility in the photovoltaic power generation output data 182 of the storage unit 180 at a time delayed from the estimation time. Although the photovoltaic power generation output data 182 will be described here as the measurement data, as described above, an estimated value may be used instead of the measurement data.

次に、推定部150は、係数αを取得する(ステップS43)。詳細には、推定部150は、記憶部180に算出データ185として格納されている係数αを読み出すことで、係数αを取得する。 Next, the estimation unit 150 acquires the coefficient α (step S43). Specifically, the estimation unit 150 acquires the coefficient α by reading the coefficient α stored as the calculated data 185 in the storage unit 180.

次に、推定部150は、係数αと太陽光発電出力計測データを用いて推定対象の太陽光発電設備の発電出力を推定し(ステップS44)、処理を終了する。詳細には、ステップS44では、第1太陽光発電出力であるステップS42で取得した太陽光発電出力計測データを用いて上記式(6)により第2太陽光発電出力を推定し、推定した結果である発電出力推定値を記憶部180の算出データ185として格納する。 Next, the estimation unit 150 estimates the power generation output of the photovoltaic power generation facility to be estimated using the coefficient α and the photovoltaic power generation output measurement data (step S44), and ends the process. Specifically, in step S44, the second photovoltaic power generation output is estimated by the above equation (6) using the photovoltaic power generation output measurement data acquired in step S42, which is the first photovoltaic power generation output, and the estimated result is used. A certain power generation output estimated value is stored as the calculated data 185 of the storage unit 180.

本実施の形態の太陽光発電出力推定装置100は、以上の処理により、太陽光発電出力を取得可能な太陽光発電設備が少なくとも1つあれば、当該太陽光発電設備の発電出力を用いることで、日射計および高い時間分解能を有する計測器を設置することなく、他の所望の太陽光発電設備の太陽光発電出力を個別に推定することが可能である。さらに、本実施の形態では、上述したように、自己係数λを用いて係数αの検証を行うことにより係数αの推定精度の低下を防ぐことができるため、第2太陽光発電出力を高精度に推定することができる。 The photovoltaic power generation output estimation device 100 of the present embodiment can use the photovoltaic power generation output of the photovoltaic power generation facility if there is at least one photovoltaic power generation facility capable of acquiring the photovoltaic power generation output by the above processing. It is possible to estimate the PV output of other desired PV equipment individually without installing a solarometer and a measuring instrument with high time resolution. Further, in the present embodiment, as described above, by verifying the coefficient α using the self-coefficient λ, it is possible to prevent the estimation accuracy of the coefficient α from deteriorating, so that the second photovoltaic power generation output is highly accurate. Can be estimated to.

なお、1つの第2需要家の第2太陽光発電設備の発電出力を1つの第1需要家の第1太陽光発電設備の発電出力を用いて推定する例を説明したが、複数の第2太陽光発電設備の発電出力の総量(合算値)を1つの第1需要家の第1太陽光発電設備の発電出力を用いて推定してもよい。また、複数の第2太陽光発電設備の発電出力の総量(合算値)を複数の第1需要家の第1太陽光発電設備の発電出力の総量(合算値)を用いて推定してもよい。例えば、柱上トランス単位などの配電区間単位、配電線単位、送電線単位などのように、複数の需要家の太陽光発電設備の発電出力の総量を推定する必要がある場合、発電出力の計測が行われていない需要家の発電出力の総量を、同様に推定することができる。 Although an example of estimating the power generation output of the second photovoltaic power generation facility of one second consumer using the power generation output of the first photovoltaic power generation facility of one first consumer has been described, a plurality of second ones have been described. The total amount (total value) of the power generation output of the photovoltaic power generation facility may be estimated using the power generation output of the first photovoltaic power generation facility of one first consumer. Further, the total amount (total value) of the power generation output of the plurality of second photovoltaic power generation facilities may be estimated using the total amount (total value) of the power generation output of the first photovoltaic power generation facility of the plurality of first consumers. .. For example, when it is necessary to estimate the total amount of power generation output of photovoltaic power generation equipment of multiple consumers, such as distribution section unit such as pillar transformer unit, distribution line unit, transmission line unit, etc., power generation output measurement The total amount of power generation output of consumers who have not done so can be estimated in the same way.

図12は、複数の需要家の太陽光発電設備を推定対象とする場合の第1太陽光発電設備および第2太陽光発電設備の一例を示す図である。図12に示した例では、図5に示した例と同様に、需要家200-i(i=1,2,…,8)の太陽光発電設備202を、PViとして示している。実線の枠で示したPV1,PV3は、発電出力が計測されており、破線の枠で示したPV2などは、発電出力が計測されていない太陽光発電設備202である。図12に示した例では、PV2,4~8は発電出力の推定対象の第2太陽光発電設備であり、PV1,3は発電出力が計測される第1太陽光発電設備である。この場合、PV2,4~8の発電出力の総量を推定対象の第2太陽光発電出力とし、PV1,3の発電出力の総量を第1太陽光発電出力として、1つの第2需要家の第2太陽光発電設備の発電出力を1つの第1需要家の第1太陽光発電設備の発電出力を用いて推定する例と同様に発電出力の推定を行うことができる。この場合も、PV1~PV8が予め定められた距離を半径とする範囲303内に存在すれば、第1太陽光発電設備と第2太陽光発電設備との間の距離は予め定められた距離以下となる。 FIG. 12 is a diagram showing an example of a first solar power generation facility and a second solar power generation facility when the solar power generation facilities of a plurality of consumers are targeted for estimation. In the example shown in FIG. 12, the photovoltaic power generation facility 202 of the consumer 200-i (i = 1, 2, ..., 8) is shown as PVi, as in the example shown in FIG. PV1 and PV3 shown in the solid line frame have the power generation output measured, and PV2 and the like shown in the broken line frame are the photovoltaic power generation facilities 202 in which the power generation output has not been measured. In the example shown in FIG. 12, PV2, 4 to 8 are the second photovoltaic power generation facilities whose power generation output is estimated, and PV1 and PV3 are the first photovoltaic power generation facilities whose power generation output is measured. In this case, the total amount of power generation output of PV2, 4 to 8 is defined as the second photovoltaic power generation output to be estimated, and the total amount of power generation output of PV1 and PV1 is defined as the first photovoltaic power generation output. (2) The power generation output can be estimated in the same manner as in the example of estimating the power generation output of the photovoltaic power generation facility using the power generation output of the first photovoltaic power generation facility of one first consumer. In this case as well, if PV1 to PV8 are within the range 303 whose radius is a predetermined distance, the distance between the first photovoltaic power generation facility and the second photovoltaic power generation facility is equal to or less than the predetermined distance. Will be.

図13は、複数の需要家の太陽光発電設備を推定対象とする場合の係数αの一例を示す図である。図13に示したように、複数の第2太陽光発電設備の発電出力を複数の第1需要家の第1太陽光発電設備の発電出力を用いて推定する場合、複数の太陽光発電設備の組み合わせごとに係数αが算出される。なお、このような場合には、遅延時間は、第1太陽光発電設備の設置地点の重心的位置と第2太陽光発電設備の設置地点の重心的位置との間を日射変動が伝播する時間となる。なお、遅延時間を0としてもよい。図13では、配電区間ごとに係数αを定めて配電区間ごとに発電出力を推定する例を示しているが、需要家を統合する単位はこれに限定されない。 FIG. 13 is a diagram showing an example of the coefficient α when the solar power generation equipment of a plurality of consumers is targeted for estimation. As shown in FIG. 13, when the power generation output of the plurality of second photovoltaic power generation facilities is estimated by using the power generation output of the first photovoltaic power generation facility of the plurality of first consumers, the power generation output of the plurality of photovoltaic power generation facilities is estimated. The coefficient α is calculated for each combination. In such a case, the delay time is the time for the solar radiation fluctuation to propagate between the position of the center of gravity of the installation point of the first photovoltaic power generation facility and the position of the center of gravity of the installation point of the second photovoltaic power generation facility. Will be. The delay time may be set to 0. FIG. 13 shows an example in which a coefficient α is set for each distribution section and the power generation output is estimated for each distribution section, but the unit for integrating consumers is not limited to this.

次に、太陽光発電出力推定装置100が奏する効果の検証結果を説明する。図14は、本実施の形態の発電出力の推定方法の検証結果の一例を示す図である。図14において、横軸は時刻を示しており、縦軸は第2太陽光発電出力(有効電力)を示している。破線は、本実施の形態で述べた推定方法により推定された第2太陽光発電出力の推定値を示し、実線は、実測値を示している。図14に示した例では、複数軒の全量買取契約の対象である需要家の太陽光発電設備を第1太陽光発電設備とし、複数軒の余剰買取契約の対象である需要家の太陽光発電設備を第2太陽光発電設備として、第2太陽光発電設備の発電出力の総量を推定した結果を示している。つまり、図14に示した例では、第1太陽光発電設備のスマートメータ204による計測データを第1太陽光発電出力として用い、第2太陽光発電設備を有する第2需要家のスマートメータ205の計測データを残余需要として用いて、第2太陽光発電設備の発電出力の総量を推定している。また、図14に示した例では、過去のある日の8時から16時までの30単位の16点をそれぞれ推定時点として、当該推定時点の日の前日に係数αを求め、推定時点の当日の第1太陽光発電出力の計測データを用いて第2太陽光発電出力を推定した。この検証では、第2太陽光発電設備の発電出力を別途特別に計測することで図14に求めた実測値を算出している。図14からわかるように、本実施の形態の発電出力の推定方法により、第2太陽光発電出力を高精度に推定できることが確認できた。 Next, the verification result of the effect of the photovoltaic power generation output estimation device 100 will be described. FIG. 14 is a diagram showing an example of verification results of the power generation output estimation method of the present embodiment. In FIG. 14, the horizontal axis indicates the time, and the vertical axis indicates the second photovoltaic power generation output (active power). The broken line shows the estimated value of the second photovoltaic power generation output estimated by the estimation method described in the present embodiment, and the solid line shows the measured value. In the example shown in FIG. 14, the solar power generation equipment of the consumer who is the target of the total purchase contract of the plurality of houses is set as the first solar power generation equipment, and the solar power generation of the consumer who is the target of the surplus purchase contract of the plurality of houses. The result of estimating the total amount of the power generation output of the 2nd photovoltaic power generation facility is shown with the facility as the 2nd photovoltaic power generation facility. That is, in the example shown in FIG. 14, the measurement data by the smart meter 204 of the first photovoltaic power generation facility is used as the first photovoltaic power generation output, and the smart meter 205 of the second consumer having the second photovoltaic power generation facility Using the measured data as the residual demand, the total amount of power generation output of the second photovoltaic power generation facility is estimated. Further, in the example shown in FIG. 14, the coefficient α is obtained on the day before the day of the estimation time, with 16 points of 30 units from 8:00 to 16:00 on a certain day in the past as the estimation time, and the day of the estimation time. The second photovoltaic power generation output was estimated using the measurement data of the first photovoltaic power generation output of. In this verification, the measured value obtained in FIG. 14 is calculated by separately specially measuring the power generation output of the second photovoltaic power generation facility. As can be seen from FIG. 14, it was confirmed that the second photovoltaic power generation output can be estimated with high accuracy by the power generation output estimation method of the present embodiment.

図15は、推定精度の検証を行わない場合の本実施の形態の発電出力の推定方法の検証結果の一例を示す図である。図15において、図14と同様に、破線は、本実施の形態で述べた推定方法により推定された第2太陽光発電出力の推定値を示し、実線は、実測値を示している。図15に示した例では、図14に示した例と同様の検証を別の日に行っている。ただし、図15に示した例では、自己係数λを用いた推定精度の検証を行っていない。このように、自己係数λを用いた推定精度の検証を行わない場合、第2太陽光発電出力の推定精度が低下する場合がある。 FIG. 15 is a diagram showing an example of verification results of the power generation output estimation method of the present embodiment when the estimation accuracy is not verified. In FIG. 15, similarly to FIG. 14, the broken line indicates the estimated value of the second photovoltaic power generation output estimated by the estimation method described in the present embodiment, and the solid line indicates the actually measured value. In the example shown in FIG. 15, the same verification as in the example shown in FIG. 14 is performed on another day. However, in the example shown in FIG. 15, the estimation accuracy using the self-coefficient λ is not verified. As described above, if the estimation accuracy using the self-coefficient λ is not verified, the estimation accuracy of the second photovoltaic power generation output may decrease.

ここで、推定精度の低下の要因について説明する。図16は、図14に検証結果が示された日の前日の検証結果の一例を示す図であり、図17は、図15に検証結果が示された日の前日の検証結果の一例を示す図である。図14および図15に結果を示した検証では、上述したように前日に係数αの推定を行っている。図16と図17を比較するとわかるように、図17に対応する日は、図16に対応する日に比べ、日射強度が非常に弱い。このため、図15に示した推定値を求めるために使用された係数αの算出に用いられる第1太陽光発電出力、すなわち図17に対応する日の第1太陽光発電出力も非常に小さくなる。したがって、Cov[PPV1(t),PPV1(t+τ)]が小さくなるため、相対的に、Cov[PPV1(t),PL2(t)]の値を0と近似して無視することができず、上述した式(5)の導出時に前提とされた仮定が成り立たたなくなり、係数αの推定精度が低下する。これにより、図15に示した例では、第2太陽光発電出力の推定精度が低下することになる。本実施の形態では、自己係数λを用いて、|λ-1|が閾値uを超える場合には、係数αの推定値を採用しないようにしているため、図15に示すような推定精度の低下を防ぐことができる。 Here, the factors that reduce the estimation accuracy will be described. FIG. 16 is a diagram showing an example of the verification result on the day before the day when the verification result is shown in FIG. 14, and FIG. 17 shows an example of the verification result on the day before the day when the verification result is shown in FIG. It is a figure. In the verification showing the results in FIGS. 14 and 15, the coefficient α is estimated the day before as described above. As can be seen by comparing FIGS. 16 and 17, the day corresponding to FIG. 17 has a much weaker illuminance intensity than the day corresponding to FIG. Therefore, the first photovoltaic power generation output used for calculating the coefficient α used to obtain the estimated value shown in FIG. 15, that is, the first photovoltaic power generation output on the day corresponding to FIG. 17 is also very small. .. Therefore, since the Cov t [P PV1 (t), P PV1 (t + τ s )] becomes small, the value of the Cov t [P PV1 (t), PL2 (t)] is relatively close to 0. It cannot be ignored, and the assumptions made when deriving the above-mentioned equation (5) do not hold, and the estimation accuracy of the coefficient α decreases. As a result, in the example shown in FIG. 15, the estimation accuracy of the second photovoltaic power generation output is lowered. In the present embodiment, the self-coefficient λ is used, and when | λ-1 | exceeds the threshold value u, the estimated value of the coefficient α is not adopted. Therefore, the estimation accuracy as shown in FIG. 15 is achieved. It is possible to prevent the decrease.

図18は、|λ-1|と第2太陽光発電出力の推定誤差との散布図の一例を示す図である。図18において、横軸は|λ-1|を示しており、縦軸は、第2太陽光発電出力の推定値の2乗平均平方根誤差(RMSE:Root Mean Squared Error)を示している。図18に示した散布図は、図14および図15に示した例と同様の検証を複数の日に関して実施し、複数の日に関してそれぞれ推定値の実測値に対するRMSEと|λ-1|を算出して得られた結果を示している。誤差が大きく推定精度が低い日は、図18に示すように|λ-1|も大きくなってなり、推定精度が高くなる日は|λ-1|が0に近い値になる。したがって、|λ-1|が閾値uを超える場合に、係数αの算出を行わないことで推定精度の低下を防ぐことができる。図18に示した例では、閾値uを0.5に設定している。閾値uは、第2太陽光発電出力の推定精度が所望の精度すなわち定められた精度を満たすように決定される。これにより、図4に示したステップS6の判定でYesと判定された場合に、第2太陽光発電出力の推定精度が所望の精度を満たすようになる。 FIG. 18 is a diagram showing an example of a scatter diagram of | λ-1 | and the estimation error of the second photovoltaic power generation output. In FIG. 18, the horizontal axis indicates | λ-1 |, and the vertical axis indicates the root-mean squared error (RMSE) of the estimated value of the second photovoltaic power generation output. In the scatter plot shown in FIG. 18, the same verification as in the examples shown in FIGS. 14 and 15 was performed for a plurality of days, and RMSE and | λ-1 | for the measured values of the estimated values were calculated for each of the plurality of days. The results obtained are shown. On days when the error is large and the estimation accuracy is low, | λ-1 | also becomes large as shown in FIG. 18, and on days when the estimation accuracy is high, | λ-1 | becomes a value close to 0. Therefore, when | λ-1 | exceeds the threshold value u, it is possible to prevent a decrease in estimation accuracy by not calculating the coefficient α. In the example shown in FIG. 18, the threshold value u is set to 0.5. The threshold value u is determined so that the estimation accuracy of the second photovoltaic power generation output satisfies a desired accuracy, that is, a predetermined accuracy. As a result, when the determination in step S6 shown in FIG. 4 determines Yes, the estimation accuracy of the second photovoltaic power generation output satisfies the desired accuracy.

閾値uの決定方法に特に制約はないが、例えば、推定精度検証部140が、検証により得られたRMSEと|λ-1|とを記憶部180に格納しておき、表示部160が、これらを散布図として表示する。すなわち、表示部16が、図18に例示した散布図を表示する。なお、このときのRMSEは、|λ-1|を用いた推定精度の検証を行わず、自己係数λの値によらずに全ての場合で係数αを更新して発電出力を推定して得られた結果とする。そして、ユーザがこの散布図を確認しながら、推定誤差であるRMSEが大きいグループと小さいグループとを適切に分離できるように閾値uを決定し、決定した閾値uを、入力受付部170を介して入力する。例えば、図18に示した閾値uの直線を入力受付部170の一例であるマウスなどにより左右に動かすことでユーザが閾値uを調整する。そして、例えば、図18に例示した散布図に確定ボタンを表示しておき、ユーザが、閾値uを示す直線を所望の位置に移動させて確定ボタンを押下することで閾値uが設定されるようにしてもよい。なお、閾値uの決定に利用する需要家は、太陽光発電出力を計測している第1需要家にしてもよいし、第2需要家に仮設的に計測器を設置して、太陽光発電出力を計測することでもよい。さらには、第2需要家に設置された太陽光発電出力の推定値を利用してもよい。ユーザによる閾値uの設定方法はこの例に限定されない。また、検証結果の代わりにシミュレーション結果などが用いられてもよい。 The method for determining the threshold value u is not particularly limited. For example, the estimation accuracy verification unit 140 stores the RMSE and | λ-1 | obtained by the verification in the storage unit 180, and the display unit 160 stores them. Is displayed as a scatter plot. That is, the display unit 16 displays the scatter plot illustrated in FIG. The RMSE at this time is obtained by updating the coefficient α and estimating the power generation output in all cases regardless of the value of the self-coefficient λ without verifying the estimation accuracy using | λ-1 |. The result is. Then, while the user confirms this scatter plot, the threshold value u is determined so that the group having a large RMSE and the group having a small estimation error can be appropriately separated, and the determined threshold value u is set via the input reception unit 170. input. For example, the user adjusts the threshold value u by moving the straight line of the threshold value u shown in FIG. 18 left and right with a mouse or the like, which is an example of the input receiving unit 170. Then, for example, a confirmation button is displayed on the scatter diagram illustrated in FIG. 18, and the user moves the straight line indicating the threshold value u to a desired position and presses the confirmation button so that the threshold value u is set. You may do it. The consumer used to determine the threshold value u may be the first consumer measuring the photovoltaic power generation output, or a measuring instrument may be temporarily installed in the second consumer to generate photovoltaic power generation. It may be possible to measure the output. Furthermore, the estimated value of the photovoltaic power generation output installed in the second consumer may be used. The method of setting the threshold value u by the user is not limited to this example. Further, a simulation result or the like may be used instead of the verification result.

閾値uは、機械学習などの計算により、太陽光発電出力推定装置100または他の装置により算出されてもよい。例えば、推定精度検証部140が、検証により得られたRMSEと|λ-1|とを用いて、機械学習により閾値uを決定してもよい。すなわち、推定精度検証部140は、第2太陽光発電出力の実測値に対する、推定部150により推定された第2太陽光発電出力の推定値の誤差(RMSE)を複数取得し、誤差と誤差に対応する自己係数λとで構成されるデータセットを複数用いて、|λ-1|が閾値u未満である場合の誤差が定められた値以下となるように機械学習により閾値を算出してもよい。機械学習としては、最大エントロピー法、決定木などを用いてもよいし、クラスタリングを行う手法により結果を分類することで閾値uを求めてもよく、どのような手法を用いてもよい。表示部160は、上述した散布図、すなわち第2太陽光発電出力の実測値に対する推定部150により推定された第2太陽光発電出力の推定値の誤差との関係を示す散布図を表示してもよい。さらに、表示部160は、この散布図に、計算により算出された閾値uを重畳させて例えば図18のように表示してもよい。これにより、ユーザが、閾値uが正しく設定されていかを確認することができる。また、ユーザが、計算により算出された閾値uを変更できるようにしてもよい。例えば、ユーザが閾値uを設定するときと同様に、計算により算出された閾値uを示す直線を左右にずらすことで、閾値uを変更できるようにしてもよい。 The threshold value u may be calculated by the photovoltaic power generation output estimation device 100 or another device by calculation such as machine learning. For example, the estimation accuracy verification unit 140 may determine the threshold value u by machine learning using the RMSE obtained by the verification and | λ-1 |. That is, the estimation accuracy verification unit 140 acquires a plurality of errors (RMSE) of the estimated value of the second photovoltaic power generation output estimated by the estimation unit 150 with respect to the measured value of the second photovoltaic power generation output, and uses the error and the error. Even if the threshold is calculated by machine learning so that the error when | λ-1 | is less than the threshold u is equal to or less than the specified value using a plurality of data sets composed of the corresponding self-factors λ. good. As machine learning, the maximum entropy method, a decision tree, or the like may be used, or the threshold value u may be obtained by classifying the results by a clustering method, and any method may be used. The display unit 160 displays the above-mentioned scatter diagram, that is, a scatter diagram showing the relationship with the error of the estimated value of the second photovoltaic power generation output estimated by the estimation unit 150 with respect to the measured value of the second photovoltaic power generation output. May be good. Further, the display unit 160 may superimpose the threshold value u calculated by calculation on the scatter plot and display it as shown in FIG. 18, for example. This allows the user to confirm whether the threshold value u is set correctly. Further, the user may be able to change the threshold value u calculated by calculation. For example, the threshold value u may be changed by shifting the straight line indicating the threshold value u calculated by calculation to the left or right, as in the case of setting the threshold value u by the user.

上記の例では、自己係数λを用いて推定精度の検証を行ったが、第1太陽光発電設備以外の太陽光発電設備202の残余需要を用いて係数αに相当する係数λ´を計算し、係数λ´を自己係数λの代わりに用いて係数αの推定に適しているか否かを判定してもよい。係数λ´も推定精度の検証に用いられる検証係数である。例えば、図1に示した需要家200-1の太陽光発電設備202が第1太陽光発電設備であり、図1に示した需要家200-2の太陽光発電設備202が第2太陽光発電設備であるとする。このとき、需要家200-3の太陽光発電設備202も計測器206により計測されており、需要家200-3の太陽光発電設備202の発電出力も記憶部180に、太陽光発電出力データ182として格納されている。需要家200-3の太陽光発電設備202を第3太陽光発電設備とし、需要家200-3を第3需要家とし、第3需要家の見かけ上の消費電力である残余需要をP(t)とすると、係数λ´は、以下の式(8)により算出することができる。 In the above example, the estimation accuracy was verified using the self-coefficient λ, but the coefficient λ'corresponding to the coefficient α was calculated using the residual demand of the photovoltaic power generation facility 202 other than the first photovoltaic power generation facility 202. , The coefficient λ'may be used instead of the self-coefficient λ to determine whether or not it is suitable for estimating the coefficient α. The coefficient λ'is also a verification coefficient used for verification of estimation accuracy. For example, the solar power generation facility 202 of the consumer 200-1 shown in FIG. 1 is the first solar power generation facility, and the solar power generation facility 202 of the consumer 200-2 shown in FIG. 1 is the second solar power generation facility. It is assumed that it is a facility. At this time, the photovoltaic power generation facility 202 of the consumer 200-3 is also measured by the measuring instrument 206, and the power generation output of the photovoltaic power generation facility 202 of the consumer 200-3 is also stored in the storage unit 180, and the photovoltaic power generation output data 182. It is stored as. The photovoltaic power generation facility 202 of the consumer 200-3 is the third photovoltaic power generation facility, the consumer 200-3 is the third consumer, and the residual demand, which is the apparent power consumption of the third consumer, is P 3 ( If t), the coefficient λ'can be calculated by the following equation (8).

Figure 2022084398000009
Figure 2022084398000009

係数λ´は、第1太陽光発電出力を第3太陽光発電設備の発電出力である第3太陽光発電出力に変換するための係数であるため、第1太陽光発電出力に係数λ´を乗じると、第3太陽光発電出力の推定値が得られる。一方、第3太陽光発電出力は、需要家200-3において計測器206により計測されており、理想的には、この計測データと推定値とが一致することになる。したがって、係数λ´は、理想的には第3太陽光発電出力(計測データ)であるPPV3(t)を、第1太陽光発電出力PPV1(t)で割った値になる。このため、図4のステップS6では、第3太陽光発電出力を第1太陽光発電出力で除した値を減じた値を検証係数から減じた結果の絶対値である|λ´-PPV3(t)/PPV1(t)|が閾値u以下であるか否かを判定することで、自己係数λを用いた場合と同様に、推定精度が低下している場合の係数αの採用を防ぐことができる。PPV3(t)/PPV1(t)は、第1期間において複数得られるため平均値を用いてもよいし中央値を用いてもよい。また、以下の式(9)で計算できるPPV3(t)とPPV1(t)の回帰係数ρを利用してもよい。 Since the coefficient λ'is a coefficient for converting the first photovoltaic power generation output into the third photovoltaic power generation output, which is the power generation output of the third photovoltaic power generation facility, the coefficient λ'is added to the first photovoltaic power generation output. Multiplying gives an estimate of the third PV output. On the other hand, the third photovoltaic power generation output is measured by the measuring instrument 206 in the consumer 200-3, and ideally, the measured data and the estimated value match. Therefore, the coefficient λ'is ideally a value obtained by dividing P PV3 (t), which is the third photovoltaic power generation output (measurement data), by the first photovoltaic power generation output P PV1 (t). Therefore, in step S6 of FIG. 4, it is an absolute value of the result of subtracting the value obtained by subtracting the value obtained by dividing the third photovoltaic power generation output by the first photovoltaic power generation output from the verification coefficient | λ'-P PV3 ( By determining whether or not t) / P PV1 (t) | is equal to or less than the threshold value u, it is possible to prevent the adoption of the coefficient α when the estimation accuracy is low, as in the case of using the self-coefficient λ. be able to. Since a plurality of P PV3 (t) / P PV1 (t) are obtained in the first period, the average value may be used or the median value may be used. Further, the regression coefficients ρ of P PV3 (t) and P PV1 (t) that can be calculated by the following equation (9) may be used.

Figure 2022084398000010
Figure 2022084398000010

このように、第3太陽光発電設備の発電出力である第3太陽光発電出力を用いて第2太陽光発電出力の推定精度が高いか否かを判定することができる。第2太陽光発電設備および第3太陽光発電設備は、第1太陽光発電設備と予め定められた距離内に設置されている。 In this way, it is possible to determine whether or not the estimation accuracy of the second photovoltaic power generation output is high by using the third photovoltaic power generation output, which is the power generation output of the third photovoltaic power generation facility. The second photovoltaic power generation facility and the third photovoltaic power generation facility are installed within a predetermined distance from the first photovoltaic power generation facility.

すなわち、本実施の形態の推定部150は、第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、第1太陽光発電設備の発電出力である第1太陽光発電出力と第2太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて推定する。そして、推定精度検証部140は、第1太陽光発電設備から予め定められた距離内に設置された第3太陽光発電設備の発電出力と第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定する。 That is, the estimation unit 150 of the present embodiment first sets the second photovoltaic power generation output, which is the power output of the second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility. It is estimated using the first photovoltaic power generation output, which is the power generation output of the photovoltaic power generation facility, and the residual demand, which is the apparent power consumption of the consumer in which the second photovoltaic power generation facility is installed. Then, the estimation accuracy verification unit 140 sees the power output of the third photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility and the apparentness of the consumer in which the third photovoltaic power generation facility is installed. It is determined whether or not the estimation accuracy of the second photovoltaic power generation output satisfies the predetermined accuracy by using the residual demand which is the above power consumption.

詳細には、第1共分散算出部120は、検証係数である係数λ´の算出のために、第1期間における、第1太陽光発電出力の時系列データと第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要の時系列データとの共分散である第3共分散を算出する。このとき、係数αの算出時と同様に、遅延時間も算出される。この遅延時間は、第1太陽光発電設備の設置地点と第3太陽光発電設備の設置地点との間を日射変動が伝播する時間であり、第2遅延時間とも呼ぶ。第2共分散算出部130は、第2遅延時間を用いて、係数αの算出時と同様に、第1期間における第1太陽光発電出力の時系列データと第1期間から第2遅延時間だけずれた第3期間の第1太陽光発電出力の時系列データとの自己共分散である第4共分散を算出する。第3共分散および第4共分散は、第3太陽光発電設備を第2太陽光発電設備としたときの係数αの算出における第1共分散および第2共分散にそれぞれ相当する。したがって、第3共分散を第1共分散ということもでき、第4共分散を第2共分散ということもできる。推定精度検証部140は、第3共分散および第4共分散を用いて第1太陽光発電出力と第3太陽光発電出力との比の推定値に相当する検証係数である係数λ´を算出し、係数λ´と、第3太陽光発電出力の第1太陽光発電出力に対する比、あるいは、回帰係数とを用いて、第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定する。そして、推定部150は、推定精度検証部140によって第2太陽光発電出力の推定精度が定められた精度を満たすと判定された場合は係数αを用いた第2太陽光発電出力の推定を行い、推定精度検証部140によって第2太陽光発電出力の推定精度が定められた精度を満たさないと判定された場合は、係数の初期値または過去に算出された係数αを用いて第2太陽光発電出力の推定を行う。第3太陽光発電設備は、第1太陽光発電設備であってもよいし第1太陽光発電設備とは異なる太陽光発電設備であってもよい。 Specifically, in the first covariance calculation unit 120, the time series data of the first photovoltaic power generation output and the third photovoltaic power generation facility in the first period are installed for the calculation of the coefficient λ'which is the verification coefficient. The third covariance, which is the covariance with the time-series data of the residual demand, which is the apparent power consumption of the consumer, is calculated. At this time, the delay time is also calculated in the same manner as when the coefficient α is calculated. This delay time is the time for the solar radiation fluctuation to propagate between the installation point of the first photovoltaic power generation facility and the installation point of the third photovoltaic power generation facility, and is also called the second delay time. Using the second delay time, the second covariance calculation unit 130 uses only the time-series data of the first solar power output in the first period and the second delay time from the first period, as in the case of calculating the coefficient α. The fourth covariance, which is the autocovariance with the time-series data of the first solar power output in the shifted third period, is calculated. The third covariance and the fourth covariance correspond to the first covariance and the second covariance in the calculation of the coefficient α when the third photovoltaic power generation facility is used as the second photovoltaic power generation facility, respectively. Therefore, the third covariance can be referred to as the first covariance, and the fourth covariance can be referred to as the second covariance. The estimation accuracy verification unit 140 calculates a coefficient λ'which is a verification coefficient corresponding to an estimated value of the ratio between the first photovoltaic power generation output and the third solar power generation output using the third co-dispersion and the fourth co-dispersion. Then, whether or not the estimation accuracy of the second photovoltaic power generation output satisfies the predetermined accuracy by using the coefficient λ'and the ratio of the third photovoltaic power generation output to the first photovoltaic power generation output or the regression coefficient. Is determined. Then, when the estimation accuracy verification unit 140 determines that the estimation accuracy of the second photovoltaic power generation output satisfies the determined accuracy, the estimation unit 150 estimates the second photovoltaic power generation output using the coefficient α. If it is determined by the estimation accuracy verification unit 140 that the estimation accuracy of the second photovoltaic power generation output does not satisfy the specified accuracy, the initial value of the coefficient or the coefficient α calculated in the past is used for the second photovoltaic power generation. Estimate the power output. The third photovoltaic power generation facility may be the first photovoltaic power generation facility or may be a photovoltaic power generation facility different from the first photovoltaic power generation facility.

図4を用いて説明した自己係数λを用いる例は、第3太陽光発電設備が第1太陽光発電設備である場合ということになり、自己係数λは、係数λ´の一例である。ここでは、第1太陽光発電設備が需要家200-1の太陽光発電設備202であり、第3太陽光発電設備が需要家200-3の太陽光発電設備202である例を説明したが、この逆に、第1太陽光発電設備が需要家200-3の太陽光発電設備202であり、第3太陽光発電設備が需要家200-1の太陽光発電設備202であってもよい。なお、第3太陽光発電設備と第1太陽光発電設備が同一の場合、第2遅延時間は0である。また、第3太陽光発電設備が第1太陽光発電設備と異なる場合に、第2遅延時間を0と近似してもよい。また、第1太陽光発電設備が需要家200-1の太陽光発電設備202である場合に、第1太陽光発電設備と予め定められた距離内に他の全量買取契約の対象の太陽光発電設備がある場合には、当該太陽光発電設備を第3太陽光発電設備としてもよい。さらに、複数の需要家の第1太陽光発電出力の総量を用いて第2太陽光発電出力を推定する場合、1つの第3需要家の第3太陽光発電出力および残余需要を用いて推定精度の検証を行ってもよいし、複数の第3需要家の第3太陽光発電出力の総量および残余需要の総量を用いて推定精度の検証を行ってもよい。 The example using the self-coefficient λ described with reference to FIG. 4 is the case where the third photovoltaic power generation facility is the first photovoltaic power generation facility, and the self-coefficient λ is an example of the coefficient λ ′. Here, an example has been described in which the first photovoltaic power generation facility is the photovoltaic power generation facility 202 of the consumer 200-1 and the third photovoltaic power generation facility is the photovoltaic power generation facility 202 of the consumer 200-3. On the contrary, the first photovoltaic power generation facility may be the photovoltaic power generation facility 202 of the consumer 200-3, and the third photovoltaic power generation facility may be the photovoltaic power generation facility 202 of the consumer 200-1. When the third photovoltaic power generation facility and the first photovoltaic power generation facility are the same, the second delay time is 0. Further, when the third photovoltaic power generation facility is different from the first photovoltaic power generation facility, the second delay time may be approximated to 0. Further, when the first photovoltaic power generation facility is the photovoltaic power generation facility 202 of the consumer 200-1, the photovoltaic power generation subject to another total purchase contract within a predetermined distance from the first photovoltaic power generation facility. If there is equipment, the solar power generation equipment may be used as a third solar power generation equipment. Furthermore, when estimating the second photovoltaic power generation output using the total amount of the first photovoltaic power generation output of a plurality of consumers, the estimation accuracy is estimated using the third photovoltaic power generation output and the residual demand of one third consumer. Or, the estimation accuracy may be verified using the total amount of the third photovoltaic power generation output of the plurality of third consumers and the total amount of the residual demand.

このように、図4を用いて説明した例は、第3太陽光発電設備が第1太陽光発電設備である場合すなわち第3太陽光発電設備と第1太陽光発電設備とが同一の場合であり、本実施の形態の太陽光発電出力推定プログラムがコンピュータに実行させる上述した判定ステップは、第3太陽光発電設備の発電出力と第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定するステップ、の一例である。 As described above, the example described with reference to FIG. 4 is when the third photovoltaic power generation facility is the first photovoltaic power generation facility, that is, when the third photovoltaic power generation facility and the first photovoltaic power generation facility are the same. Therefore, the above-mentioned determination step that the photovoltaic power generation output estimation program of the present embodiment causes the computer to execute is apparently the power output of the third photovoltaic power generation facility and the consumer in which the third photovoltaic power generation facility is installed. This is an example of a step of determining whether or not the estimation accuracy of the second photovoltaic power generation output satisfies the predetermined accuracy by using the residual demand which is the power consumption.

上述した例では、計測データである第1太陽光発電出力を用いて第2太陽光発電出力を推定するようにしたが、第2太陽光発電出力の推定後、推定された第2太陽光発電出力を第1太陽光発電出力として用いて、さらに別の太陽光発電設備の発電出力を推定してもよい。これを繰り返すことで、例えば、推定対象の太陽光発電設備202から予め定められた距離内に発電出力が計測されている太陽光発電設備202が存在しない場合であっても、推定対象の太陽光発電設備202の発電出力を推定することができる。 In the above example, the second photovoltaic power generation output is estimated using the first photovoltaic power generation output, which is the measurement data. However, after the estimation of the second photovoltaic power generation output, the estimated second photovoltaic power generation is generated. The output may be used as the first photovoltaic power generation output to estimate the power generation output of yet another photovoltaic power generation facility. By repeating this, for example, even if there is no photovoltaic power generation facility 202 whose power generation output is measured within a predetermined distance from the photovoltaic power generation facility 202 to be estimated, the solar power to be estimated is to be estimated. The power output of the power generation facility 202 can be estimated.

図19は、発電出力の推定の順序の一例を示す模式図である。図19では、需要家200-i(i=1,2,…,10)の太陽光発電設備202を、PViとして示している。実線の枠で示したPV1,PV3は、それぞれ需要家200-1,200-3の太陽光発電設備202であり、発電出力が計測されている。破線の枠で示したPV2などは、それぞれ需要家200-2など発電出力が計測されていない太陽光発電設備202である。範囲301,302は、図5に示した例と同様である。図19に示した黒矢印は発電出力の計測データを用いた推定を示し、白矢印は他の太陽光発電設備202の発電出力の推定結果を用いた発電出力の推定を示す。図5に示すように、PV2,PV5,PV6は、PV1の発電出力の計測データを用いて発電出力の推定が行われ、PV4,PV7,PV8は、PV3の発電出力の計測データを用いて発電出力の推定が行われる。その後、PV6の発電出力の推定値を第1太陽光発電出力として用いて、第2太陽光発電設備であるPV8の発電出力の推定が行われ、PV4の発電出力の推定値を第1太陽光発電出力として用いて、第2太陽光発電設備であるPV10の発電出力の推定が行われる。さらに、PV8の発電出力の推定値を第1太陽光発電出力として用いて、第2太陽光発電設備であるPV9の発電出力の推定が行われる。 FIG. 19 is a schematic diagram showing an example of the order of estimation of power generation output. In FIG. 19, the photovoltaic power generation equipment 202 of the consumer 200-i (i = 1, 2, ..., 10) is shown as PVi. PV1 and PV3 shown in the solid line frame are the photovoltaic power generation facilities 202 of the consumers 200-1 and 200-3, respectively, and the power generation output is measured. PV2 and the like shown by the broken line frame are photovoltaic power generation facilities 202 such as the consumer 200-2 whose power generation output is not measured. The ranges 301 and 302 are the same as the example shown in FIG. The black arrow shown in FIG. 19 indicates the estimation using the measurement data of the power generation output, and the white arrow indicates the estimation of the power generation output using the estimation result of the power generation output of the other photovoltaic power generation equipment 202. As shown in FIG. 5, PV2, PV5, and PV6 are estimated for power generation output using the measurement data of the power generation output of PV1, and PV4, PV7, and PV8 generate power using the measurement data of the power generation output of PV3. The output is estimated. After that, the estimated value of the power generation output of PV6 is used as the first photovoltaic power generation output, the power generation output of PV8, which is the second photovoltaic power generation facility, is estimated, and the estimated value of the power generation output of PV4 is used as the first photovoltaic power generation. The power generation output of PV10, which is the second photovoltaic power generation facility, is estimated by using it as the power generation output. Further, the estimated value of the power generation output of PV8 is used as the first solar power generation output to estimate the power generation output of PV9, which is the second solar power generation facility.

PV8,PV9,PV10は、PV1から予め定められた距離内に位置せず、PV3から予め定められた距離内に位置しないが、推定値を用いた発電出力の推定が順次行われることで、PV8,PV9,PV10の発電出力の推定を行うことができる。 PV8, PV9, and PV10 are not located within a predetermined distance from PV1 and are not located within a predetermined distance from PV3. However, by sequentially estimating the power generation output using the estimated value, PV8 , PV9, PV10 can be estimated.

また、スマートメータ203,204,205の計測データは、例えば30分の積算値などの積算値であってもよいし、瞬時値であってもよい。例えば、係数αの算出処理では、積算値が用いられ、推定時点における係数αを用いた第2太陽光発電出力の推定時に第1太陽光発電出力の瞬時値が用いられてもよい。 Further, the measurement data of the smart meters 203, 204, 205 may be an integrated value such as an integrated value for 30 minutes, or may be an instantaneous value. For example, in the calculation process of the coefficient α, the integrated value is used, and the instantaneous value of the first photovoltaic power generation output may be used when estimating the second photovoltaic power generation output using the coefficient α at the time of estimation.

以上に説明したように、実施の形態1の太陽光発電出力推定装置100は、第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、第1太陽光発電設備の発電出力である第1太陽光発電出力と第2太陽光発電設備が設置された需要家の残余需要とを用いて推定する推定部150を備える。従って、実施の形態1の太陽光発電出力推定装置100によれば、太陽光発電出力を取得可能な太陽光発電設備が少なくとも1つあれば、当該太陽光発電設備を第1太陽光発電設備とすることで、日射計および高い時間分解能を有する計測器を設置することなく、他の所望の太陽光発電設備の太陽光発電出力を個別に推定することが可能である。 As described above, the photovoltaic power generation output estimation device 100 of the first embodiment is the power output of the second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility. 2 The estimation unit 150 that estimates the photovoltaic power generation output using the first photovoltaic power generation output, which is the power output of the first photovoltaic power generation facility, and the residual demand of the consumer in which the second photovoltaic power generation facility is installed. Be prepared. Therefore, according to the photovoltaic power generation output estimation device 100 of the first embodiment, if there is at least one photovoltaic power generation facility capable of acquiring the photovoltaic power generation output, the photovoltaic power generation facility is referred to as the first photovoltaic power generation facility. By doing so, it is possible to individually estimate the photovoltaic power output of other desired photovoltaic power generation equipment without installing a photovoltaic meter and a measuring instrument having high time resolution.

さらに、太陽光発電出力推定装置100は、第1太陽光発電設備から予め定められた距離内に設置された第3太陽光発電設備の発電出力と第3太陽光発電設備が設置された需要家の残余需要とを用いて第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定する推定精度検証部140を備える。このため、推定精度が定められた精度を満たさないと判定された場合に、対策を講じることができる。この対策の一例は、算出された係数αを採用せず初期値または前回算出された係数αを用いることであるが、対策はこれに限定されない。例えば、推定精度が定められた精度を満たさないと判定された場合に別の推定方法により第2太陽光発電出力を推定するなどであってもよい。推定精度が定められた精度を満たさないと判定された場合に、対策を行うことで、第2太陽光発電出力の推定精度の低下を防ぐことができる。 Further, the photovoltaic power generation output estimation device 100 includes the power output of the third photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility and the consumer in which the third photovoltaic power generation facility is installed. It is provided with an estimation accuracy verification unit 140 for determining whether or not the estimation accuracy of the second photovoltaic power generation output satisfies a predetermined accuracy by using the residual demand of. Therefore, if it is determined that the estimation accuracy does not satisfy the predetermined accuracy, countermeasures can be taken. An example of this measure is to use the initial value or the previously calculated coefficient α without adopting the calculated coefficient α, but the measure is not limited to this. For example, when it is determined that the estimation accuracy does not satisfy the predetermined accuracy, the second photovoltaic power generation output may be estimated by another estimation method. When it is determined that the estimation accuracy does not satisfy the predetermined accuracy, it is possible to prevent the estimation accuracy of the second photovoltaic power generation output from deteriorating by taking measures.

実施の形態2.
次に、実施の形態2にかかる太陽光発電出力推定装置100について説明する。本実施の形態の太陽光発電出力推定装置100の機能構成およびハードウェア構成は、実施の形態1の太陽光発電出力推定装置100と同様の構成である。ただし、第1共分散算出部120および推定精度検証部140の動作が一部実施の形態1と異なる。実施の形態1と同様の機能を有する構成要素は実施の形態1と同一の符号を付して説明し、実施の形態1と重複する説明を省略する。以下、実施の形態1と異なる点を主に説明する。
Embodiment 2.
Next, the photovoltaic power generation output estimation device 100 according to the second embodiment will be described. The functional configuration and hardware configuration of the photovoltaic power generation output estimation device 100 of the present embodiment are the same as those of the photovoltaic power generation output estimation device 100 of the first embodiment. However, the operations of the first covariance calculation unit 120 and the estimation accuracy verification unit 140 are partially different from those of the first embodiment. The components having the same functions as those of the first embodiment will be described with the same reference numerals as those of the first embodiment, and the description overlapping with the first embodiment will be omitted. Hereinafter, the differences from the first embodiment will be mainly described.

実施の形態2にかかる太陽光発電出力推定装置100は、実施の形態1と同様に図4に示した手順で係数αの推定処理を実施するが、実施の形態2では、ステップS7の係数αの算出処理において、αが正の値であるかの判定が追加される。 The photovoltaic power generation output estimation device 100 according to the second embodiment carries out the estimation process of the coefficient α by the procedure shown in FIG. 4 in the same manner as in the first embodiment, but in the second embodiment, the coefficient α in step S7. In the calculation process of, the determination of whether α is a positive value is added.

図20は、実施の形態2の係数αの算出処理手順の一例を示すフローチャートである。図20に示したステップS31~ステップS32は、実施の形態1と同様である。ステップS33aでは、実施の形態1とステップS33と同様に係数αの算出が行われる(ステップS33a)。この時点で係数αの更新は行わない。ステップS33aの係数αの算出の後に、推定部150は係数αが正の値であるか否かを判定する(ステップS34)。係数αは、第2太陽光発電出力と第1太陽光発電出力との比であるため、正の値になるはずである。係数αが正の値にならない場合には、係数αの推定精度が低下していると考えられる。このため、係数αが正の値でない場合(ステップS34 No)、推定部150は、係数αの更新を行わずに、処理を終了する。係数αが正の値である場合(ステップS34 Yes)、推定部150は、記憶部180の算出データ185内の係数αを、算出した係数αで更新し(ステップS35)、処理を終了する。 FIG. 20 is a flowchart showing an example of the calculation processing procedure of the coefficient α of the second embodiment. Steps S31 to S32 shown in FIG. 20 are the same as those in the first embodiment. In step S33a, the coefficient α is calculated in the same manner as in the first embodiment and step S33 (step S33a). At this point, the coefficient α is not updated. After calculating the coefficient α in step S33a, the estimation unit 150 determines whether or not the coefficient α is a positive value (step S34). Since the coefficient α is the ratio of the second photovoltaic power generation output to the first photovoltaic power generation output, it should be a positive value. If the coefficient α does not become a positive value, it is considered that the estimation accuracy of the coefficient α is lowered. Therefore, when the coefficient α is not a positive value (step S34 No.), the estimation unit 150 ends the process without updating the coefficient α. When the coefficient α is a positive value (step S34 Yes), the estimation unit 150 updates the coefficient α in the calculated data 185 of the storage unit 180 with the calculated coefficient α (step S35), and ends the process.

また、上記の例では、ステップS34で係数αが0より大きいか否かを判定しているが、この代わりに、ステップS34において、係数αが0より大きくかつ定格容量比×k以下であるかを判定してもよい。定格容量比は、第1太陽光発電設備の定格容量に対する第2太陽光発電設備の定格容量の比である。kは、1より大きく、例えば1.5程度とすることができるが、検証などによって適宜定められればよくこの値に限定されない。係数αは、第1太陽光発電出力に乗算することで第2太陽光発電出力を求めるための変換係数であることから、係数αが定格比と大きく異なる場合には、係数αの精度が低下している可能性がある。このため、ステップS34において、係数αが0より大きくかつ定格容量比×k以下となるかを判定し、係数αが0より大きくかつ定格容量比×k以下となる場合にステップS35を実施することで、第2太陽光発電出力の推定精度の劣化を防ぐことができる。また、ステップS34で、係数αが定格容量比×k以上でありかつ定格容量比×k以下であるかを判定するようにしてもよい。kは1より小さく、kは1より大きい。kは、例えば、1.5程度であり、kは例えば0.5程度であるが、k,kは、検証などによって適宜定められていればよくこれらの値に限定されない。 Further, in the above example, whether or not the coefficient α is larger than 0 is determined in step S34, but instead, in step S34, is the coefficient α larger than 0 and the rated capacity ratio × k or less? May be determined. The rated capacity ratio is the ratio of the rated capacity of the second photovoltaic power generation facility to the rated capacity of the first photovoltaic power generation facility. k can be larger than 1, for example, about 1.5, but is not limited to this value as long as it is appropriately determined by verification or the like. Since the coefficient α is a conversion coefficient for obtaining the second photovoltaic power generation output by multiplying the first photovoltaic power generation output, the accuracy of the coefficient α decreases when the coefficient α is significantly different from the rated ratio. It may be. Therefore, in step S34, it is determined whether the coefficient α is larger than 0 and the rated capacity ratio × k or less, and when the coefficient α is larger than 0 and the rated capacity ratio × k or less, step S35 is performed. Therefore, it is possible to prevent deterioration of the estimation accuracy of the second photovoltaic power generation output. Further, in step S34, it may be determined whether the coefficient α is the rated capacity ratio × k 1 or more and the rated capacity ratio × k 2 or less. k 1 is less than 1 and k 2 is greater than 1. For example, k 2 is about 1.5 and k 1 is about 0.5, but k 1 and k 2 are not limited to these values as long as they are appropriately determined by verification or the like.

すなわち、推定部150は、係数αが定格比に第1値であるkを乗算した値より小さい場合、および係数αが定格比に第2値であるkを乗算した値より大きい場合に、係数の初期値または過去に算出された係数αを用いて第2太陽光発電出力の推定を行ってもよい。 That is, when the coefficient α is smaller than the value obtained by multiplying the rated ratio by the first value k 1 , and when the coefficient α is larger than the value obtained by multiplying the rated ratio by the second value k 2 . , The second solar power output may be estimated using the initial value of the coefficient or the coefficient α calculated in the past.

以上述べた以外の本実施の形態の動作は実施の形態1と同様である。本実施の形態では、実施の形態1と同様に、自己係数λを用いた推定精度の検証を行うとともに、係数αの値に基づいた推定精度の検証も実施する。これにより、実施の形態1と同様の効果が得られるとともに、さらなる推定精度の向上を図ることができる。 The operation of the present embodiment other than that described above is the same as that of the first embodiment. In the present embodiment, as in the first embodiment, the estimation accuracy is verified using the self-coefficient λ, and the estimation accuracy is also verified based on the value of the coefficient α. As a result, the same effect as that of the first embodiment can be obtained, and the estimation accuracy can be further improved.

実施の形態3.
次に、実施の形態3にかかる太陽光発電出力推定装置100について説明する。本実施の形態の太陽光発電出力推定装置100の機能構成およびハードウェア構成は、実施の形態1の太陽光発電出力推定装置100と同様の構成である。ただし、推定部150の動作が一部実施の形態1と異なる。実施の形態1と同様の機能を有する構成要素は実施の形態1と同一の符号を付して説明し、実施の形態1と重複する説明を省略する。以下、実施の形態1と異なる点を主に説明する。
Embodiment 3.
Next, the photovoltaic power generation output estimation device 100 according to the third embodiment will be described. The functional configuration and hardware configuration of the photovoltaic power generation output estimation device 100 of the present embodiment are the same as those of the photovoltaic power generation output estimation device 100 of the first embodiment. However, the operation of the estimation unit 150 is partially different from that of the first embodiment. The components having the same functions as those of the first embodiment will be described with the same reference numerals as those of the first embodiment, and the description overlapping with the first embodiment will be omitted. Hereinafter, the differences from the first embodiment will be mainly described.

図21は、本実施の形態の太陽光発電出力推定装置100による係数αの推定処理手順の一例を示すフローチャートである。ステップS1は、実施の形態1の図4に示したステップS1と同様である。ステップS1の後、第1共分散算出部120は、複数の第1期間を設定する(ステップS2a)。その後、複数の第1期間のうちの1つの期間に関して、実施の形態1と同様に、ステップS3~ステップS6が行われる。ステップS7aでは、実施の形態1と同様に係数αの算出が行われるが、算出された係数αは推定部150または記憶部180に一時的に記憶され、記憶部180の算出データ内の係数αの更新はこの時点では行われない。ステップS7aの後、第1共分散算出部120は、全ての第1期間の係数αの算出が行われたか否かを判断し(ステップS9)、係数αの算出が行われていない第1期間がある場合(ステップS9 No)には、係数αの算出が行われていない第1期間のうちの1つを選択する。以降、選択された第1期間に関してステップS3からの処理が繰り返される。 FIG. 21 is a flowchart showing an example of an estimation processing procedure for the coefficient α by the photovoltaic power generation output estimation device 100 of the present embodiment. Step S1 is the same as step S1 shown in FIG. 4 of the first embodiment. After step S1, the first covariance calculation unit 120 sets a plurality of first periods (step S2a). Then, as in the first embodiment, steps S3 to S6 are performed for one of the plurality of first periods. In step S7a, the coefficient α is calculated as in the first embodiment, but the calculated coefficient α is temporarily stored in the estimation unit 150 or the storage unit 180, and the coefficient α in the calculated data of the storage unit 180 is stored. Is not updated at this point. After step S7a, the first covariance calculation unit 120 determines whether or not the coefficients α of all the first periods have been calculated (step S9), and the first period in which the coefficients α have not been calculated. If there is (step S9 No), one of the first periods in which the coefficient α is not calculated is selected. After that, the process from step S3 is repeated for the selected first period.

全ての第1期間の係数αの算出が行われた場合(ステップS9 Yes)、第1共分散算出部120は推定部150に代表値の算出を指示し、推定部150は、係数αの代表値を算出し、記憶されている係数αを算出した値で更新する(ステップS10)。詳細には、推定部150は、一時的に保持されている第1期間ごとの係数αを用いて、係数αの代表値を算出する。係数αの代表値は、算出された複数の係数αの中央値であることが望ましいが、平均値または最頻値のような統計量であってもよい。ステップS8は、実施の形態1と同様である。 When all the coefficients α of the first period are calculated (step S9 Yes), the first covariance calculation unit 120 instructs the estimation unit 150 to calculate the representative value, and the estimation unit 150 is the representative of the coefficient α. A value is calculated, and the stored coefficient α is updated with the calculated value (step S10). Specifically, the estimation unit 150 calculates a representative value of the coefficient α by using the coefficient α for each first period that is temporarily held. The representative value of the coefficient α is preferably the median value of the plurality of calculated coefficients α, but may be a statistic such as an average value or a mode value. Step S8 is the same as that of the first embodiment.

なお、ステップS7aの前にステップS6において、第1期間が係数αの算出に適しているかすなわち推定精度は所望の精度を満たすかの判定が行われている。第1期間が係数αの算出に適していないと判定されると、係数αが算出されないため、この場合、算出される係数αの数は、ステップS2aで設定された第1期間の数より少なくなる。例えば、ステップS2aで10個の第1期間が設定され、そのうち1つの期間で係数αの算出に適していないという判定が行われたとする。この場合、推定部150は、ステップS10では、9つの係数αを保持していることになるため、9つの係数αを用いて代表値を決定する。このように複数の第1期間を設定して、第1期間が係数αの算出に適している場合に対応する係数αを算出し、算出した係数αを用いて代表値を算出し、算出した代表値を採用するようにした。このため、稀に発生する係数αの外れ値、すなわち推定誤差が大きい係数αの影響により発生する、第2太陽光発電出力の推定精度の低下を軽減することができる。 In addition, before step S7a, in step S6, it is determined whether the first period is suitable for the calculation of the coefficient α, that is, whether the estimation accuracy satisfies the desired accuracy. If it is determined that the first period is not suitable for calculating the coefficient α, the coefficient α is not calculated. Therefore, in this case, the number of calculated coefficients α is smaller than the number of the first period set in step S2a. Become. For example, it is assumed that 10 first periods are set in step S2a, and it is determined that one of the periods is not suitable for calculating the coefficient α. In this case, since the estimation unit 150 holds the nine coefficients α in step S10, the estimation unit 150 determines the representative value using the nine coefficients α. In this way, a plurality of first periods are set, the coefficient α corresponding to the case where the first period is suitable for the calculation of the coefficient α is calculated, and the representative value is calculated and calculated using the calculated coefficient α. The representative value is adopted. Therefore, it is possible to reduce the decrease in the estimation accuracy of the second photovoltaic power generation output, which is caused by the influence of the coefficient α, which rarely occurs, that is, the deviation value of the coefficient α, that is, the estimation error is large.

このように、第1期間は複数設定されてもよく少なくとも1つ設定されればよい。すなわち、第1共分散算出部120は、第2太陽光発電出力の推定時点より前の少なくとも1つの第1期間における、第1太陽光発電出力の時系列データと第2需要家の残余需要の時系列データとの共分散である第1共分散を算出すればよい。第2共分散算出部130、推定部150、推定精度検証部140は、各第1期間に関して実施の形態1に述べた動作を実施する。第1期間が複数設定されるときには、推定部150は、|λ-1|が閾値u未満であると判定された第1期間に対応する係数αを用いて係数αの代表値を算出し、代表値を用いて第2太陽光発電出力を推定する。 In this way, a plurality of first periods may be set, and at least one may be set. That is, the first covariance calculation unit 120 determines the time series data of the first photovoltaic power generation output and the residual demand of the second consumer in at least one first period before the estimation time of the second photovoltaic power generation output. The first covariance, which is the covariance with the time series data, may be calculated. The second covariance calculation unit 130, the estimation unit 150, and the estimation accuracy verification unit 140 carry out the operation described in the first embodiment for each first period. When a plurality of first periods are set, the estimation unit 150 calculates a representative value of the coefficient α using the coefficient α corresponding to the first period in which | λ-1 | is determined to be less than the threshold value u. The second photovoltaic output is estimated using the representative value.

なお、実施の形態2と本実施の形態を組み合わせて、さらに、係数αが正の値でない場合に係数αを更新しないようにしてもよい。この場合、推定部150は、係数αが正の値の場合のみ係数αを保持し、ステップS10で保持している係数αの代表値を求めてもよいし、係数αが負の値であっても一時的に保持しておき、ステップS10において正の値の係数αを用いて代表値を求めてもよい。すなわち、推定部150は、|λ-1|が閾値u未満であると判定された第1期間に対応する係数αのうち正の値の係数αを用いて係数αの代表値を算出し、代表値を用いて第2太陽光発電出力を推定することになる。同様に、係数αが0より大きくかつ定格容量比×k以下であるかの判定と組み合わせてもよいし、係数αが定格容量比×k以上でありかつ定格容量比×k以下であるかの判定と組み合わせてもよい。また、図21では、自己係数λを使用して係数αを採用するか否かを判定する例を示したが、実施の形態1で述べたように係数λ´を用いて係数αを採用するか否かを判定してもよい。 It should be noted that the second embodiment and the present embodiment may be combined so that the coefficient α is not updated when the coefficient α is not a positive value. In this case, the estimation unit 150 may hold the coefficient α only when the coefficient α is a positive value, and may obtain a representative value of the coefficient α held in step S10, or the coefficient α is a negative value. However, it may be temporarily held and the representative value may be obtained by using the coefficient α of the positive value in step S10. That is, the estimation unit 150 calculates a representative value of the coefficient α by using the coefficient α of a positive value among the coefficients α corresponding to the first period in which | λ-1 | is determined to be less than the threshold u. The second photovoltaic power generation output will be estimated using the representative value. Similarly, it may be combined with the determination of whether the coefficient α is larger than 0 and the rated capacity ratio × k or less, or the coefficient α is the rated capacity ratio × k 1 or more and the rated capacity ratio × k 2 or less. It may be combined with the determination. Further, in FIG. 21, an example of determining whether or not to adopt the coefficient α by using the self-coefficient λ is shown, but as described in the first embodiment, the coefficient α is adopted by using the coefficient λ'. It may be determined whether or not.

以上の実施の形態に示した構成は、一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、実施の形態同士を組み合わせることも可能であるし、要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。 The configuration shown in the above embodiments is an example, and can be combined with another known technique, can be combined with each other, and does not deviate from the gist. It is also possible to omit or change a part of the configuration.

10 太陽光発電出力推定システム、100 太陽光発電出力推定装置、110 データ取得部、120 第1共分散算出部、130 第2共分散算出部、140 推定精度検証部、150 推定部、160 表示部、170 入力受付部、180 記憶部、200-1~200-4,200-m,200-n 需要家、201 負荷、202 太陽光発電設備、203~205 スマートメータ、206 計測器。 10 Photovoltaic power generation output estimation system, 100 Photovoltaic power generation output estimation device, 110 Data acquisition unit, 120 1st covariance calculation unit, 130 2nd covariance calculation unit, 140 estimation accuracy verification unit, 150 estimation unit, 160 display unit , 170 Input reception unit, 180 storage unit, 200-1 to 200-4, 200-m, 200-n consumer, 201 load, 202 photovoltaic power generation equipment, 203-205 smart meter, 206 measuring instrument.

Claims (18)

第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、前記第1太陽光発電設備の発電出力である第1太陽光発電出力と前記第2太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて推定する推定部と、
前記第1太陽光発電設備から前記予め定められた距離内に設置された第3太陽光発電設備の発電出力と前記第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて前記第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定する推定精度検証部と、
を備えることを特徴とする太陽光発電出力推定装置。
The second photovoltaic power generation output, which is the power generation output of the second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility, is the power generation output of the first photovoltaic power generation facility. An estimation unit that estimates using the photovoltaic power generation output and the residual demand, which is the apparent power consumption of the consumer in which the second photovoltaic power generation facility is installed,
It is the power generation output of the third photovoltaic power generation facility installed within the predetermined distance from the first photovoltaic power generation facility and the apparent power consumption of the consumer in which the third photovoltaic power generation facility is installed. An estimation accuracy verification unit that determines whether or not the estimation accuracy of the second photovoltaic power generation output satisfies the specified accuracy using the residual demand, and
A photovoltaic power generation output estimation device characterized by being equipped with.
前記第2太陽光発電出力の推定時点より前の少なくとも1つの第1期間における、前記第1太陽光発電出力の時系列データと前記第2太陽光発電設備が設置された需要家の残余需要の時系列データとの共分散である第1共分散を算出する第1共分散算出部と、
前記第1期間における前記第1太陽光発電出力の時系列データと前記第1期間から第1遅延時間だけずれた第2期間の前記第1太陽光発電出力の時系列データとの自己共分散である第2共分散を算出する第2共分散算出部と、
を備え、
前記第1遅延時間は、第1太陽光発電設備と前記第2太陽光発電設備との間で日射変動が伝搬する時間であり、
前記推定部は、前記第1共分散および前記第2共分散を用いて前記第1太陽光発電出力と前記第2太陽光発電出力との比の推定値である係数を算出し、前記係数を前記推定時点から遅延時間ずれた前記第1太陽光発電出力に乗じることで前記第2太陽光発電出力を推定し、
前記推定精度検証部は、前記第1期間ごとに前記第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定することを特徴とする請求項1に記載の太陽光発電出力推定装置。
Time-series data of the first photovoltaic power generation output and residual demand of the consumer in which the second photovoltaic power generation facility is installed in at least one first period before the estimation time of the second photovoltaic power generation output. The first covariance calculation unit that calculates the first covariance, which is the covariance with the time series data,
By self-covariance between the time-series data of the first photovoltaic power generation output in the first period and the time-series data of the first photovoltaic power generation output in the second period deviated by the first delay time from the first period. A second covariance calculation unit that calculates a certain second covariance,
Equipped with
The first delay time is the time for the solar radiation fluctuation to propagate between the first photovoltaic power generation facility and the second photovoltaic power generation facility.
The estimation unit calculates a coefficient which is an estimated value of the ratio between the first photovoltaic power generation output and the second photovoltaic power generation output by using the first co-dispersion and the second co-dispersion, and calculates the coefficient. The second photovoltaic power generation output is estimated by multiplying the first photovoltaic power generation output whose delay time deviates from the estimation time point.
The photovoltaic power generation according to claim 1, wherein the estimation accuracy verification unit determines whether or not the estimation accuracy of the second solar power generation output satisfies a predetermined accuracy for each first period. Output estimator.
前記第1共分散算出部は、前記第1期間における、前記第1太陽光発電出力の時系列データと前記第3太陽光発電設備が設置された需要家の残余需要の時系列データとの共分散である第3共分散を算出し、
前記第2共分散算出部は、前記第1期間における前記第1太陽光発電出力の時系列データと前記第1期間から第2遅延時間だけずれた第3期間の前記第1太陽光発電出力の時系列データとの自己共分散である第4共分散を算出し、
前記第2遅延時間は、第1太陽光発電設備と前記第3太陽光発電設備との間で日射変動が伝搬する時間であり、
前記推定精度検証部は、前記第3共分散および前記第4共分散を用いて前記第1太陽光発電出力と前記第3太陽光発電設備の発電出力である第3太陽光発電出力との比の推定値である検証係数を算出し、前記検証係数と、前記第3太陽光発電出力の前記第1太陽光発電出力に対する比または回帰係数とを用いて、前記第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定し、
前記推定部は、前記推定精度検証部によって前記第2太陽光発電出力の推定精度が定められた精度を満たすと判定された場合は前記係数を用いた前記第2太陽光発電出力の推定を行い、前記推定精度検証部によって前記第2太陽光発電出力の推定精度が定められた精度を満たさないと判定された場合は、前記係数の初期値または過去に算出された前記係数を用いて前記第2太陽光発電出力の推定を行うことを特徴とする請求項2に記載の太陽光発電出力推定装置。
The first covariance calculation unit is a combination of the time-series data of the first photovoltaic power generation output and the time-series data of the residual demand of the consumer in which the third photovoltaic power generation facility is installed in the first period. Calculate the third covariance, which is the variance,
The second covariance calculation unit is the time-series data of the first photovoltaic power generation output in the first period and the first photovoltaic power generation output in the third period deviated from the first period by the second delay time. Calculate the fourth covariance, which is the autocovariance with the time series data,
The second delay time is the time for the solar radiation fluctuation to propagate between the first photovoltaic power generation facility and the third photovoltaic power generation facility.
The estimation accuracy verification unit uses the third co-dispersion and the fourth co-dispersion to compare the first photovoltaic power generation output with the third photovoltaic power generation output, which is the power generation output of the third photovoltaic power generation facility. The verification coefficient, which is an estimated value of, is calculated, and the ratio or regression coefficient of the third photovoltaic power generation output to the first photovoltaic power generation output is used to estimate the second photovoltaic power generation output. Judging whether the accuracy meets the specified accuracy,
When the estimation unit determines that the estimation accuracy of the second photovoltaic power generation output satisfies the predetermined accuracy by the estimation accuracy verification unit, the estimation unit estimates the second photovoltaic power generation output using the coefficient. If the estimation accuracy verification unit determines that the estimation accuracy of the second photovoltaic power generation output does not satisfy the determined accuracy, the first value of the coefficient or the coefficient calculated in the past is used. 2. The photovoltaic power generation output estimation device according to claim 2, wherein the photovoltaic power generation output is estimated.
前記係数は、前記第1共分散を前記第2共分散で除して-1を乗じた値であり、
前記検証係数は、前記第3共分散を前記第4共分散で除して-1を乗じた値であることを特徴とする請求項3に記載の太陽光発電出力推定装置。
The coefficient is a value obtained by dividing the first covariance by the second covariance and multiplying by -1.
The photovoltaic power generation output estimation device according to claim 3, wherein the verification coefficient is a value obtained by dividing the third covariance by the fourth covariance and multiplying by -1.
前記推定部は、前記係数が正の値でない場合に、前記係数の初期値または過去に算出された前記係数を用いて前記第2太陽光発電出力の推定を行うことを特徴とする請求項3または4に記載の太陽光発電出力推定装置。 3. The estimation unit is characterized in that, when the coefficient is not a positive value, the second photovoltaic power generation output is estimated using the initial value of the coefficient or the coefficient calculated in the past. Or the photovoltaic power generation output estimation device according to 4. 前記推定部は、前記係数が定格比に第1値を乗算した値より小さい場合、および前記係数が前記定格比に第2値を乗算した値より大きい場合に、前記係数の初期値または過去に算出された前記係数を用いて前記第2太陽光発電出力の推定を行い、
前記定格比は、前記第1太陽光発電設備の定格容量に対する前記第2太陽光発電設備の定格容量の比であり、前記第1値は0以上かつ1未満であり、前記第2値は1より大きいことを特徴とする請求項3または4に記載の太陽光発電出力推定装置。
The estimation unit determines the initial value of the coefficient or in the past when the coefficient is smaller than the value obtained by multiplying the rated ratio by the first value, and when the coefficient is larger than the value obtained by multiplying the rated ratio by the second value. Using the calculated coefficient, the second photovoltaic power generation output is estimated.
The rated ratio is the ratio of the rated capacity of the second photovoltaic power generation facility to the rated capacity of the first photovoltaic power generation facility, the first value is 0 or more and less than 1, and the second value is 1. The photovoltaic power generation output estimation device according to claim 3 or 4, characterized in that it is larger.
前記推定精度検証部は、前記第3太陽光発電出力を第1太陽光発電出力で除した値を減じた値を前記検証係数から減じた結果の絶対値が閾値以下である場合に前記第2太陽光発電出力の推定精度が定められた精度を満たすと判定し、前記絶対値が閾値以下である場合に前記第2太陽光発電出力の推定精度が定められた精度を満たさないと判定することを特徴とする請求項4に記載の太陽光発電出力推定装置。 The second estimation accuracy verification unit is the second when the absolute value of the result of subtracting the value obtained by subtracting the value obtained by dividing the third photovoltaic power generation output by the first photovoltaic power generation output from the verification coefficient is equal to or less than the threshold value. It is determined that the estimation accuracy of the photovoltaic power generation output satisfies the defined accuracy, and when the absolute value is equal to or less than the threshold value, it is determined that the estimation accuracy of the second photovoltaic power generation output does not satisfy the defined accuracy. The photovoltaic power generation output estimation device according to claim 4, wherein the photovoltaic power generation output is estimated. 前記第1期間は複数設定され、
前記推定部は、前記絶対値が閾値以下であると判定された前記第1期間に対応する前記係数を用いて前記係数の代表値を算出し、前記代表値を用いて前記第2太陽光発電出力を推定することを特徴とする請求項7に記載の太陽光発電出力推定装置。
A plurality of the first periods are set,
The estimation unit calculates a representative value of the coefficient using the coefficient corresponding to the first period in which the absolute value is determined to be equal to or less than the threshold value, and uses the representative value to generate the second photovoltaic power generation. The photovoltaic power generation output estimation device according to claim 7, wherein the output is estimated.
前記推定部は、前記絶対値が閾値以下であると判定された前記第1期間に対応する前記係数のうち正の値の前記係数を用いて前記係数の代表値を算出し、前記代表値を用いて前記第2太陽光発電出力を推定することを特徴とする請求項8に記載の太陽光発電出力推定装置。 The estimation unit calculates a representative value of the coefficient by using the coefficient having a positive value among the coefficients corresponding to the first period in which the absolute value is determined to be equal to or less than the threshold value, and obtains the representative value. The photovoltaic power generation output estimation device according to claim 8, wherein the second photovoltaic power generation output is estimated by using the device. 前記推定部は、前記絶対値が閾値以下であると判定された前記第1期間に対応する前記係数のうち、定格比に第1値を乗算した値より大きくかつ前記定格比に第2値を乗算した値以下の前記係数を用いて前記係数の代表値を算出し、前記代表値を用いて前記第2太陽光発電出力を推定し、
前記定格比は、前記第1太陽光発電設備の定格容量に対する前記第2太陽光発電設備の定格容量の比であり、前記第1値は0以上かつ1未満であり、前記第2値は1より大きいであることを特徴とする請求項8に記載の太陽光発電出力推定装置。
The estimation unit sets the second value to the rated ratio, which is larger than the value obtained by multiplying the rated ratio by the first value among the coefficients corresponding to the first period in which the absolute value is determined to be equal to or less than the threshold value. A representative value of the coefficient is calculated using the coefficient equal to or less than the multiplied value, and the second solar power output is estimated using the representative value.
The rated ratio is the ratio of the rated capacity of the second photovoltaic power generation facility to the rated capacity of the first photovoltaic power generation facility, the first value is 0 or more and less than 1, and the second value is 1. The photovoltaic power generation output estimation device according to claim 8, which is characterized by being larger.
前記推定精度検証部は、前記第2太陽光発電出力の実測値に対する、前記推定部により推定された前記第2太陽光発電出力の推定値の誤差を複数取得し、前記誤差と前記誤差に対応する前記検証係数とで構成されるデータセットを複数用いて、前記絶対値が前記閾値以下である場合の前記誤差が定められた値以下となるように機械学習により前記閾値を算出することを特徴とする請求項7から10のいずれか1つに記載の太陽光発電出力推定装置。 The estimation accuracy verification unit acquires a plurality of errors of the estimated value of the second photovoltaic power generation output estimated by the estimation unit with respect to the measured value of the second photovoltaic power generation output, and corresponds to the error and the error. It is characterized in that the threshold value is calculated by machine learning so that the error when the absolute value is equal to or less than the threshold value is equal to or less than a predetermined value by using a plurality of data sets composed of the verification coefficient. The photovoltaic power generation output estimation device according to any one of claims 7 to 10. 前記検証係数と、前記第2太陽光発電出力の実測値に対する前記推定部により推定された前記第2太陽光発電出力の推定値の誤差との関係を示す散布図を表示する表示部、
を備えることを特徴とする請求項7から11のいずれか1つに記載の太陽光発電出力推定装置。
A display unit that displays a scatter diagram showing the relationship between the verification coefficient and the error of the estimated value of the second photovoltaic power generation output estimated by the estimation unit with respect to the measured value of the second photovoltaic power generation output.
The photovoltaic power generation output estimation device according to any one of claims 7 to 11, wherein the solar power generation output estimation device is provided.
前記表示部は、前記散布図に重畳して前記閾値を表示することを特徴とする請求項12に記載の太陽光発電出力推定装置。 The photovoltaic power generation output estimation device according to claim 12, wherein the display unit displays the threshold value superimposed on the scatter diagram. 前記1遅延時間および前記第2遅延時間を0と近似することを特徴とする請求項3から13のいずれか1つに記載の太陽光発電出力推定装置。 The photovoltaic power generation output estimation device according to any one of claims 3 to 13, wherein the one delay time and the second delay time are approximated to 0. 前記第3太陽光発電設備は前記第1太陽光発電設備であることを特徴とする請求項1から14のいずれか1つに記載の太陽光発電出力推定装置。 The photovoltaic power generation output estimation device according to any one of claims 1 to 14, wherein the third photovoltaic power generation facility is the first photovoltaic power generation facility. 前記第3太陽光発電設備は前記第1太陽光発電設備とは異なる太陽光発電設備であることを特徴とする請求項1から14のいずれか1つに記載の太陽光発電出力推定装置。 The photovoltaic power generation output estimation device according to any one of claims 1 to 14, wherein the third photovoltaic power generation facility is a photovoltaic power generation facility different from the first photovoltaic power generation facility. 太陽光発電出力推定装置が、第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、前記第1太陽光発電設備の発電出力である第1太陽光発電出力と前記第2太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて推定し、
前記太陽光発電出力推定装置が、前記第1太陽光発電設備から前記予め定められた距離内に設置された第3太陽光発電設備の発電出力と前記第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて前記第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定することを特徴とする太陽光発電出力推定方法。
The photovoltaic power generation output estimation device converts the second photovoltaic power generation output, which is the power generation output of the second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility, into the first photovoltaic power generation. Estimated using the first photovoltaic power generation output, which is the power generation output of the facility, and the residual demand, which is the apparent power consumption of the consumer in which the second photovoltaic power generation facility is installed.
The demand that the photovoltaic power generation output estimation device installs the power output of the third photovoltaic power generation facility installed within the predetermined distance from the first photovoltaic power generation facility and the third photovoltaic power generation facility. A photovoltaic power generation output estimation method, characterized in that it is determined whether or not the estimation accuracy of the second photovoltaic power generation output satisfies a predetermined accuracy by using the residual demand which is the apparent power consumption of the house.
第1太陽光発電設備から予め定められた距離内に設置された第2太陽光発電設備の発電出力である第2太陽光発電出力を、前記第1太陽光発電設備の発電出力である第1太陽光発電出力と前記第2太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて推定するステップと、
第1太陽光発電設備から前記予め定められた距離内に設置された第3太陽光発電設備の発電出力と前記第3太陽光発電設備が設置された需要家のみかけ上の消費電力である残余需要とを用いて前記第2太陽光発電出力の推定精度が定められた精度を満たすか否かを判定するステップと、
をコンピュータに実行させることを特徴とする太陽光発電出力推定プログラム。
The second photovoltaic power generation output, which is the power generation output of the second photovoltaic power generation facility installed within a predetermined distance from the first photovoltaic power generation facility, is the power generation output of the first photovoltaic power generation facility. A step of estimating using the photovoltaic power generation output and the residual demand, which is the apparent power consumption of the consumer in which the second photovoltaic power generation facility is installed, and
The power output of the third photovoltaic power generation facility installed within the predetermined distance from the first photovoltaic power generation facility and the residual power consumption that is the apparent power consumption of the consumer in which the third photovoltaic power generation facility is installed. A step of determining whether or not the estimation accuracy of the second photovoltaic power generation output satisfies the defined accuracy by using the demand and
A photovoltaic power generation output estimation program characterized by having a computer execute.
JP2020196254A 2020-11-26 2020-11-26 Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program Active JP7450522B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020196254A JP7450522B2 (en) 2020-11-26 2020-11-26 Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020196254A JP7450522B2 (en) 2020-11-26 2020-11-26 Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program

Publications (3)

Publication Number Publication Date
JP2022084398A true JP2022084398A (en) 2022-06-07
JP2022084398A5 JP2022084398A5 (en) 2022-11-24
JP7450522B2 JP7450522B2 (en) 2024-03-15

Family

ID=81868190

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020196254A Active JP7450522B2 (en) 2020-11-26 2020-11-26 Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program

Country Status (1)

Country Link
JP (1) JP7450522B2 (en)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5547999B2 (en) 2010-03-24 2014-07-16 大阪瓦斯株式会社 Photovoltaic power generation prediction device
KR101808978B1 (en) 2014-03-31 2017-12-13 텐서 컨설팅 가부시키가이샤 Power generation system analysis device and method
JP6564264B2 (en) 2014-09-16 2019-08-21 積水化学工業株式会社 Power management system and power management method
JP6489846B2 (en) 2015-01-27 2019-03-27 関西電力株式会社 Photovoltaic power output estimation device, solar power output estimation system, and solar power output estimation method
KR101957197B1 (en) 2015-02-23 2019-03-12 엘에스산전 주식회사 Photovoltaic system
JP6782573B2 (en) 2016-07-15 2020-11-11 関西電力株式会社 Photovoltaic power generation output estimation device and photovoltaic power generation output estimation method
CN110059974B (en) 2019-04-24 2021-08-13 杭州国电电力科技发展有限公司 Photovoltaic socket type selection method based on photovoltaic optimized access

Also Published As

Publication number Publication date
JP7450522B2 (en) 2024-03-15

Similar Documents

Publication Publication Date Title
US9575475B2 (en) Systems and methods for generating an energy usage model for a building
EP3829055B1 (en) Method for calculating power generation amount reflecting farm environment and detailed technology for predicting solar farm weather, and real-time solar power generation amount prediction system using the method
JP5743881B2 (en) Power management system, power management method, customer terminal, and power management apparatus
Lorenz et al. Short term forecasting of solar irradiance by combining satellite data and numerical weather predictions
CN102687089A (en) MPPT controller,solar battery control device, solar power generation system,mppt control program,and control method for mppt controller
US20150168465A1 (en) Method and apparatus for electric power system distribution state estimations
Marchesoni-Acland et al. Intra-day solar irradiation forecast using RLS filters and satellite images
Lonij et al. Forecasts of PV power output using power measurements of 80 residential PV installs
Bashir et al. Solar-TK: A data-driven toolkit for solar PV performance modeling and forecasting
Mayer et al. Pairing ensemble numerical weather prediction with ensemble physical model chain for probabilistic photovoltaic power forecasting
JP2015138912A (en) Photovoltaic power generation amount prediction system and weather forecast system
Bosisio et al. Improving DTR assessment by means of PCA applied to wind data
Engerer et al. Himawari-8 enabled real-time distributed PV simulations for distribution networks
Cheng Hybrid solar irradiance now-casting by fusing Kalman filter and regressor
CN112780486A (en) System and method for controlling a wind power converter
Zelikman et al. Short-term solar irradiance forecasting using calibrated probabilistic models
JP6489846B2 (en) Photovoltaic power output estimation device, solar power output estimation system, and solar power output estimation method
JP6429523B2 (en) Photovoltaic power generation forecast system
JP6512745B2 (en) Photovoltaic system, terminal, analysis processing method, and program
JP2022084398A (en) Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program
Chaabene Measurements based dynamic climate observer
JP7390951B2 (en) Solar power generation output estimation device, solar power generation output estimation method, and solar power generation output estimation program
JP6763736B2 (en) Signal time synchronization device and signal time synchronization method
JP2024060369A (en) Photovoltaic power generation output estimation device, grid control system, supply and demand control system, facility formation support system, photovoltaic power generation output estimation method, and photovoltaic power generation output estimation program
EP3885860B1 (en) Method, computer-implemented tool and power plant control device for detecting power production degradation of solar power plants

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221115

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221115

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231026

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231031

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231225

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240305

R150 Certificate of patent or registration of utility model

Ref document number: 7450522

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150