JPH04372046A - Method and device for predicting demand amount - Google Patents

Method and device for predicting demand amount

Info

Publication number
JPH04372046A
JPH04372046A JP3148832A JP14883291A JPH04372046A JP H04372046 A JPH04372046 A JP H04372046A JP 3148832 A JP3148832 A JP 3148832A JP 14883291 A JP14883291 A JP 14883291A JP H04372046 A JPH04372046 A JP H04372046A
Authority
JP
Japan
Prior art keywords
function
demand
data
influencing
discontinuity
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
JP3148832A
Other languages
Japanese (ja)
Other versions
JP2681421B2 (en
Inventor
Takaharu Ishida
隆張 石田
Iwao Madori
岩男 真鳥
Shigeo Abe
阿部 重夫
Shigeru Tamura
滋 田村
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP14883291A priority Critical patent/JP2681421B2/en
Publication of JPH04372046A publication Critical patent/JPH04372046A/en
Application granted granted Critical
Publication of JP2681421B2 publication Critical patent/JP2681421B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

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

Abstract

PURPOSE:To improve the predicted accuracy of the demand. CONSTITUTION:The past recorded data of the demand amount with the influence factor is stored in a data base 2. On the basis of past record data extracted from the data base 2, a collation function I between the consecutive influence factor and the demand amount is obtained for each inconsecutive influence factor or for each combination. On the area other than the past record data, pseudo data is prepared by a pseudo data preparation means 11, and a function II with the pseudo data added is derived from a function derivation means 12 to be stored in function tables 4 and 13. A prediction execution means 7 reads out the function I or II corresponding to the predicted influence factor at the estimated objective point from a function table 4 by the prediction execution mean 7, predicting the amount of demand at the predicated objective point using the function. The addition of new past record data is desired for updating the functions I and II. In this case, when the updating is performed during the prescribed period, a function is desired to be newly obtained by initializing the functions I and II.

Description

【発明の詳細な説明】[Detailed description of the invention]

【0001】0001

【産業上の利用分野】本発明は、電力、用水、商品等の
将来の需要量予測方法及び装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and apparatus for predicting future demand for electricity, water, commercial products, etc.

【0002】0002

【従来の技術】電力や用水の需要量は天候の違い、曜日
の種類、社会的イベント、気温等の要因(以下、影響要
因と総称する。)によって大きく変動することが知られ
ている。同様に一定の商品についても、それらの影響要
因によって需要量が変動することが知られている。そし
て、それらの需要量変動に合わせて過不足なくかつ円滑
に電力、用水、商品等を供給するため、将来の需要量を
予測し、その予測に合わせて電力運用、給水運用、商品
生産・流通等の計画・管理を行うことが肝要である。 
 例えば、従来、電力の需要予測方法として、特開昭6
0−102822公報に記載されたものが提案されてい
る。これによれば、1日を複数時間帯に区分し、各時間
帯における電力需要量の実績値を、需要量の変動因子で
ある影響因子の態様と関連付けて収集し、収集したデー
タを需要量と影響因子の関係を対応させた実績データテ
ーブルの該当する要素に発生頻度として記憶しておき、
その実績データテーブルを用い、予測対象の時間帯及び
その時間帯における予測影響要因の態様に一致する要素
から需要量の発生頻度ベクトルを求め、これに基づいて
当該予測すべき時間帯における最大頻度の電力需要量を
予測値として決定するようにしている。
BACKGROUND OF THE INVENTION It is known that the amount of electricity and water demanded varies greatly depending on factors (hereinafter collectively referred to as influencing factors) such as differences in weather, types of days of the week, social events, and temperature. Similarly, it is known that the demand for certain products fluctuates depending on these influencing factors. In order to smoothly supply electricity, water, products, etc. in accordance with fluctuations in demand, we forecast future demand and adjust power operations, water supply operations, product production and distribution according to the forecast. It is important to plan and manage such matters.
For example, as a conventional power demand forecasting method,
The one described in Publication No. 0-102822 has been proposed. According to this, a day is divided into multiple time periods, the actual value of electricity demand in each time period is collected in association with the aspects of influencing factors that are fluctuation factors of demand, and the collected data is used to calculate demand The relationship between the relationship and the influencing factor is stored as the frequency of occurrence in the corresponding element of the actual data table.
Using the actual data table, find the demand occurrence frequency vector from the time period to be predicted and the elements that match the aspect of the prediction influence factor in that time period, and based on this, calculate the maximum frequency in the time period to be predicted. The power demand amount is determined as a predicted value.

【0003】また、特開昭60−106325公報に記
載されたものによれば、過去の日最高気温と電力需要量
の関係実績から、日最高気温に対する電力需要特性を最
小2乗法を用いた2次式で表現し、この関係式に予測対
象日の予想最高気温を代入することにより、予測対象日
の電力需要を予測するようにしている。
[0003] Furthermore, according to what is described in Japanese Patent Application Laid-Open No. 106325/1983, based on the past relationship between daily maximum temperature and power demand, the power demand characteristics with respect to daily maximum temperature are calculated using the least squares method. By expressing the following equation and substituting the expected maximum temperature on the prediction target day into this relational expression, the power demand on the prediction target day is predicted.

【0004】0004

【発明が解決しようとする課題】しかし、上記従来技術
によれば、次のような解決すべき課題がある。まず、特
開昭60−102822公報の技術によれば、実績デー
タテーブルを作成、更新する際に、実績データを累積的
に記憶するようにしていることから、過去のデータパタ
ーンの影響を比較的大きく受けるという問題がある。特
に、天候の急激な変化又は特異日の周辺における予測精
度が悪くなる傾向がある。
However, according to the above-mentioned prior art, there are the following problems to be solved. First, according to the technique disclosed in Japanese Patent Application Laid-open No. 60-102822, when creating and updating a performance data table, performance data is stored cumulatively, so the influence of past data patterns is comparatively reduced. There is a problem of receiving a lot of attention. In particular, prediction accuracy tends to deteriorate around sudden changes in weather or special days.

【0005】また、特開昭60−106325公報の技
術によれば、日最高気温と需要量の関係を表す2次式に
、予測対象日の予測日最高気温を単に代入して、予測対
象日の電力需要を予測するもので、過去の気温などの気
象変動による電力需要の変動を考慮していないことから
、予測精度に限界があるという問題がある。
[0005] Furthermore, according to the technique disclosed in Japanese Patent Application Laid-Open No. 106325/1983, the predicted maximum temperature on the forecasted day is simply substituted into the quadratic equation expressing the relationship between the daily maximum temperature and the amount of demand. The problem is that there is a limit to the accuracy of predictions because it does not take into account fluctuations in power demand due to weather changes such as past temperatures.

【0006】このような問題は、用水や商品の需要予測
の場合も同様である。
[0006] Such problems are similar in the case of forecasting the demand for water and products.

【0007】本発明の目的は、可能な限り最新の需要量
実績データを用い、かつ予測精度を向上できる需要量予
測方法及び装置を提供することにある。
[0007] An object of the present invention is to provide a demand forecasting method and apparatus that can improve forecast accuracy while using the latest demand record data as possible.

【0008】[0008]

【課題を解決するための手段】上記目的を達成するため
、本発明の需要量予測方法は、需要量の実績データを需
要量の増減に影響を与える複数の影響要因を付して収集
し、影響要因の変化と需要量の変化との関係が連続的な
連続性影響要因とその関係が不連続的な不連続性影響要
因とに別け、前記実績データの最近の一定期間の実績デ
ータを用い、不連続性影響要因ごとに又は不連続性影響
要因の組合せごとに前記連続性影響要因と需要量の相関
関数を求め、該関数の中から予測対象時点の予測影響要
因に対応する関数を用いて当該予測対象時点の需要量を
予測することにある。
[Means for Solving the Problems] In order to achieve the above object, the demand quantity forecasting method of the present invention collects actual demand quantity data with a plurality of influencing factors that influence increases and decreases in demand quantity, Separate into continuous influencing factors, where the relationship between changes in influencing factors and changes in demand is continuous, and discontinuous influencing factors, where the relationship is discontinuous, and use actual data for a recent certain period of the aforementioned actual data. , find a correlation function between the continuity influence factor and the demand quantity for each discontinuity influence factor or for each combination of discontinuity influence factors, and use the function corresponding to the prediction influence factor at the time of prediction from among the functions. The goal is to predict the amount of demand at the time of the forecast.

【0009】また、前記関数を求めるに際し、実績デー
タの最近の一定期間の実績データから外れる領域につい
て過去の実績データに基づいて予測対象時点の疑似デー
タを作成し、該疑似データと前記一定期間の実績データ
とを用いるようにすることが望ましい。
[0009] Furthermore, when determining the function, pseudo data at the time of prediction is created based on past performance data for areas that deviate from the performance data for a recent certain period of the performance data, and the pseudo data for the certain period are combined with the pseudo data for the certain period. It is desirable to use actual data.

【0010】また、それらを組み合わせ、最近の一定期
間の実績データを用いて第1の関数を求めるとともに、
前記疑似データと前記一定期間の実績データを用いて第
2の関数を求め、第1と第2の関数の中から予測対象時
点の予測影響要因に対応する関数を用いて当該予測対象
時点の需要量を予測することができる。
[0010]Also, by combining them and finding the first function using actual data over a recent certain period,
A second function is obtained using the pseudo data and the actual data for a certain period of time, and a function corresponding to the prediction influence factor at the time of the prediction target is used from among the first and second functions to calculate the demand at the time of the prediction target. The amount can be predicted.

【0011】さらに、予測結果に対応する需要量の新た
な実績データを前記一定期間の実績データに加えて前記
関数を更新することが望ましい。この場合、更新に係る
新たな実績データの期間が所定の追加期間に達したとき
前記関数を初期化し、該新たな実績データを含む実績デ
ータの最近の前記一定期間の実績データを用いて、前記
関数を新たに求めることが望ましい。
[0011] Furthermore, it is desirable to update the function by adding new performance data of the demand amount corresponding to the prediction result to the performance data for the certain period. In this case, when the period of new performance data related to the update reaches a predetermined additional period, the function is initialized, and the performance data of the recent certain period of performance data including the new performance data is used to It is desirable to find a new function.

【0012】前記関数を求める手段として、、重回帰モ
デル又はニューラルネットモデルを用いることができる
。その際、用いる実績データ等の信頼度に応じてモデル
の収束条件を変えたり、前記第2の関数を求める際に前
記実績データを用いる場合よりも前記疑似データを用い
る場合のモデルの収束条件を緩やかに設定することが望
ましい。
[0012] A multiple regression model or a neural network model can be used as a means for determining the function. At that time, the convergence conditions of the model may be changed depending on the reliability of the actual data used, or the convergence conditions of the model may be changed when using the pseudo data when calculating the second function rather than when using the actual data. It is desirable to set it loosely.

【0013】また、関数導出に用いる実績データを、実
際の実績データから変動しないベース分を除去したデー
タにすれば、関数の感度が高くなる。さらに、関数を導
出した後、該導出した関数から一定量外れる実績データ
を棄却し、残った実績データを用いて関数を再度導出し
たり、また導出した関数と実績データとを比較し、関数
からの実績データの外れの程度に応じて当該関数を補正
するようにすることが、予測精度向上の点で望ましい。
[0013] Furthermore, if the performance data used for deriving the function is data obtained by removing the base portion that does not fluctuate from the actual performance data, the sensitivity of the function will be increased. Furthermore, after deriving the function, we can discard the actual data that deviates from the derived function by a certain amount and derive the function again using the remaining actual data, or we can compare the derived function and the actual data and In order to improve prediction accuracy, it is desirable to correct the function according to the degree of deviation in the actual data.

【0014】本発明の需要予測装置は、上記の需要予測
方法を適用し、コンピュータ等の情報処理装置を用いて
構成することにより実現できる。
The demand forecasting device of the present invention can be realized by applying the above demand forecasting method and configuring it using an information processing device such as a computer.

【0015】すなわち、需要量の実績データが需要量の
増減に影響を与える複数の影響要因を付して格納される
データベースと、影響要因の変化と需要量の変化との関
係が不連続的な予め定められた不連続性影響要因ごと又
は不連続性影響要因の組合せごとに対応する最近の一定
期間の実績データを前記データベースから順次抽出し、
該抽出した実績データに基づいて不連続性影響要因ごと
に又は不連続性影響要因の組合せごとに、前記連続性影
響要因と需要量の相関関数を求める関数導出手段と、該
求められた関数が格納される関数テーブルと、予測対象
時点の予測影響要因を入力し該予測影響要因に対応する
関数を前記関数テーブルから読み出し、該読み出した関
数を用いて当該予測対象時点の需要量を予測する予測実
行手段とを含んで構成する。
[0015] In other words, there is a database in which actual data on demand is stored with a plurality of influencing factors that influence increases and decreases in demand, and a database in which the relationship between changes in influencing factors and changes in demand is discontinuous. Sequentially extracting performance data for a recent certain period corresponding to each predetermined discontinuity influencing factor or each combination of discontinuity influencing factors from the database,
a function deriving means for calculating a correlation function between the continuity influencing factor and the demand quantity for each discontinuity influencing factor or for each combination of discontinuity influencing factors based on the extracted performance data; A prediction that inputs a stored function table and a prediction influence factor at a prediction target time, reads a function corresponding to the prediction influence factor from the function table, and uses the read function to predict the demand amount at the prediction target time. and an execution means.

【0016】また、上記構成に、抽出した実績データか
ら外れる領域について過去の実績データに基づいて予測
対象時点の疑似データを作成する疑似データ作成手段を
加え、該疑似データを前記関数導出手段の入力データに
加えるようにする。
[0016]Furthermore, a pseudo data creation means is added to the above configuration for creating pseudo data at a prediction target point based on past performance data for a region that deviates from the extracted performance data, and the pseudo data is input to the function derivation means. Add it to the data.

【0017】[0017]

【作用】このように構成されることから、本発明によれ
ば次の作用により、本発明の目的が達成できる。
[Function] With this structure, according to the present invention, the object of the present invention can be achieved by the following function.

【0018】まず、不連続性影響要因ごとに又はその組
合せごとに、連続性影響要因と需要量の相関関数を求め
、該関数に基づいて予測対象時点の需要量を予測するこ
とから、曜日の違いや社会的なイベント等が行われる特
異日についても実績データに基づいた関数により予測で
きるので、予測精度が向上する。また、不連続性影響要
因をのぞいて、連続性影響因子についてのみの相関関数
を求めて予測することから、実測データの数が少なくて
も関数の信頼性を高くできるので、最近の一定期間の実
績データを用いることができ、過去の古い実績の影響を
必要最小限にできる。また、関数のモデルを重回帰モデ
ル又はニューラルネットを適用すれば、更に実績データ
の対象期間を短縮できる。
First, for each discontinuity influencing factor or for each combination thereof, a correlation function between the continuity influencing factor and the demand quantity is determined, and based on the function, the demand quantity at the time of the prediction target is predicted. Unique days on which differences or social events occur can also be predicted using functions based on past data, improving prediction accuracy. In addition, since predictions are made by calculating the correlation function only for continuity-affecting factors, excluding discontinuity-affecting factors, the reliability of the function can be increased even if the number of actual measurement data is small. Performance data can be used, and the influence of past performance can be minimized. Moreover, if a multiple regression model or a neural network is applied to the function model, the period covered by the performance data can be further shortened.

【0019】また、疑似データを作成するようにしたも
のによれば、実績データが少ない場合、又は予測対象時
点の連続性影響因子に対応する実績データが無い場合で
も、予測精度が向上する。
Furthermore, by creating pseudo data, the prediction accuracy is improved even when there is little actual data or when there is no actual data corresponding to the continuity influencing factor at the time of the prediction target.

【0020】また、新たな実績データを加えて関数を更
新すれば、一層予測精度が向上する。特に、更新に係る
新たな実績データの期間が所定の追加期間に達したとき
、関数を初期化して導出しなおせば、過去の古い実績デ
ータの影響を少なくできる。
[0020] Moreover, if new performance data is added and the function is updated, the prediction accuracy can be further improved. In particular, when the period of new performance data related to an update reaches a predetermined additional period, the influence of past old performance data can be reduced by initializing and re-deriving the function.

【0021】前記関数を求める手段として、重回帰モデ
ル又はニューラルネットモデルを用い、用いる実績デー
タ等の信頼度に応じてモデルの収束条件を変えたり、疑
似データを用いる場合のモデルの収束条件を緩やかに設
定することにより、妥当な予測を実現できる。
[0021] A multiple regression model or a neural network model is used as a means for determining the function, and the convergence conditions of the model are changed depending on the reliability of the actual data used, or the convergence conditions of the model are relaxed when pseudo data is used. By setting , reasonable predictions can be achieved.

【0022】[0022]

【実施例】以下、本発明を図示実施例を参照しながら説
明する。図1に、本発明を電力需要予測に適用してなる
一実施例の電力需要予測装置の構成図を示す。図示のよ
うに、実績データ入力手段1から、電力需要量の実績デ
ータがデータベース2に入力される。入力される実績デ
ータには、過去の蓄積されていた実績データの他に、運
用により得られた新たな実績データが含まれる。関数導
出手段(I)3はデータベース2内の実績データを適宜
読み出し、予め設定された手順により予測に用いる関数
を導出し、関数テーブル(I)4に格納するようになっ
ている。関数更新手段(I)5は基本的に関数導出手段
(I)3と同様の構成であるが、異なる点は、予測の実
行に応じて収集される新たな実績データをデータベース
2から読み出し、新たなデータを加味した関数を求めて
、関数テーブル内の関数を新たなものに置き換えるよう
になっている。
DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be explained below with reference to illustrated embodiments. FIG. 1 shows a configuration diagram of an electric power demand prediction device according to an embodiment in which the present invention is applied to electric power demand prediction. As shown in the figure, performance data on power demand is input into a database 2 from a performance data input means 1 . The input performance data includes new performance data obtained through operations in addition to past performance data accumulated. The function deriving means (I) 3 reads out actual performance data in the database 2 as appropriate, derives a function to be used for prediction according to a preset procedure, and stores it in a function table (I) 4. The function update means (I) 5 basically has the same configuration as the function derivation means (I) 3, but the difference is that it reads new performance data collected in accordance with the execution of prediction from the database 2, and It calculates a function that takes into account the data and replaces the function in the function table with a new one.

【0023】一方、予測の実行指令は予測指令入力手段
6から入力されるようになっている。この予測指令には
、予測対象日等の予測対象時点の指定、予測対象時点に
係る天候、曜日、気温条件等の影響要因等の予測条件が
含まれている。予測実行手段6は入力される予測条件に
該当する関数を関数テーブル(I)4から読み出し、そ
の関数を用いて予測対象時点の電力需要量を求め、予測
値として予測結果出力手段8に出力するようになってい
る。
On the other hand, a prediction execution command is inputted from prediction command input means 6. This prediction command includes designation of a prediction target time such as a prediction target date, and prediction conditions such as influencing factors such as weather, day of the week, and temperature conditions related to the prediction target time. The prediction execution means 6 reads a function corresponding to the input prediction condition from the function table (I) 4, uses the function to determine the amount of power demand at the time of the prediction target, and outputs it as a predicted value to the prediction result output means 8. It looks like this.

【0024】以上の構成が、本発明の基本構成であるが
、本実施例では更に予測の精度を向上させるため、疑似
データ作成手段11、関数導出手段(II)12、関数
テーブル(II)13及び関数更新手段(II)14が
付加されている。
The above configuration is the basic configuration of the present invention, but in this embodiment, in order to further improve prediction accuracy, pseudo data creation means 11, function derivation means (II) 12, and function table (II) 13 are provided. and function updating means (II) 14 are added.

【0025】次に、各部の詳細構成について説明する。 実績データ入力手段1から入力される実績データには、
電力需要量の増減に影響を与える複数の影響要因が付さ
れている。ここで、影響要因とは電力需要量の増減に影
響を及ぼす様々な因子をいい、これらの因子を同質の因
子群に分類した例を次に示す。また、影響因子群は、更
に影響要因の変化と需要量の変化との関係が連続的な連
続性影響要因と、その関係が不連続的な不連続性影響要
因とに別けられる。
Next, the detailed configuration of each part will be explained. The performance data input from the performance data input means 1 includes:
There are multiple influencing factors that affect the increase or decrease in electricity demand. Here, the influencing factors refer to various factors that influence increases and decreases in power demand, and an example of classifying these factors into homogeneous factor groups is shown below. In addition, the influencing factor group is further divided into continuous influencing factors, in which the relationship between changes in the influencing factor and changes in demand is continuous, and discontinuous influencing factors, in which the relationship is discontinuous.

【0026】(不連続性影響要因の例)・天候の種類(
晴、雨、曇、雪など) ・曜日の種類(平日、日祝日、土曜日、飛び石連休の中
日、工場休日など) ・一日の中の時刻 ・社会的なイベント(スポーツイベントなど)(連続性
影響要因の例) ・日最高気温 ・日最低気温 ・その他、体感温度、不快指数、湿度、日射量本発明で
は、予測の精度を上げるために、まず、不連続性の影響
要因ごと又はその組合せごとに場合分けし、その場合分
けの中で連続性の影響要因と電力需要量との相関を関数
化して予測するようにしている。
(Example of discontinuity influencing factors) Type of weather (
- Sunny, rainy, cloudy, snowy, etc.) - Type of day of the week (weekdays, Sundays and holidays, Saturdays, mid-day holidays, factory holidays, etc.) - Time of day - Social events (sporting events, etc.) (continuous (Example of factors that influence discontinuity) ・Daily maximum temperature, daily minimum temperature, other sensible temperature, discomfort index, humidity, solar radiation. Cases are divided into cases for each combination, and within each case, the correlation between the influence factor of continuity and the amount of power demand is converted into a function and predicted.

【0027】そこで、データベース2の構成は、実績デ
ータを図2に示す概念図のようにテーブル化して格納す
るようにしている。なお、以下、本実施例では、天候の
種類、曜日の種類、時刻、日最高気温、日最低気温を考
慮して、電力需要量を予測する場合について説明する。
Therefore, the database 2 is configured to store performance data in a tabular format as shown in the conceptual diagram shown in FIG. In the following, in this embodiment, a case will be described in which the amount of power demand is predicted in consideration of the type of weather, the type of day of the week, the time of day, the daily maximum temperature, and the daily minimum temperature.

【0028】関数導出手段(I)3は、不連続性影響要
因ごとに又は不連続性影響要因の組合せごとに、前記連
続性影響要因と需要量の相関関数を求めるようになって
おり、必要な実績データをデータベース2から抽出する
。本実施例では、天候の種類ごとに、かつ曜日の種類ご
とにデータを抽出し、これらの組合せのそれぞれについ
て、予測対象日の時刻に対する電力需要量と日最高気温
と日最低気温との関数を求める。なお、抽出対象の実績
データは関数導出時点から見て、最近の一定期間のもの
に制限するようにしている。
The function deriving means (I) 3 is configured to obtain a correlation function between the continuity influencing factor and the demand quantity for each discontinuity influencing factor or for each combination of discontinuity influencing factors, and performance data is extracted from the database 2. In this example, data is extracted for each type of weather and for each type of day of the week, and for each of these combinations, the functions of power demand, daily maximum temperature, and daily minimum temperature for the time of the prediction target day are calculated. demand. Note that the performance data to be extracted is limited to data for a recent certain period from the point of derivation of the function.

【0029】ここで、図3に電力需要量のある一日の変
化の例を示し、図4にある時刻における電力需要量と日
最高気温と日最低気温との関係の一例を示す。図3にお
いて、縦軸は電力需要量、横軸は時刻を示す。このよう
に変化する電力需要量を一時間ごとに予測する場合、関
数導出手段(I)3は各時刻ごとすなわち24個の関数
を導出することになる。また、年、日、分、秒ごとなど
の関数を求めるようにすることもできる。一方、図4に
示した、ある時刻における日最高気温と電力需要量の関
係曲線Hと、日最低気温と電力需要量の関係曲線Lから
明らかなように、気温と電力需要量との間にある種の非
線形な関係が観測できる。なお、図中縦軸は電力需要量
であり、横軸は気温を示す。このような傾向の関係があ
るのは、次の理由によるものと考えられる。すなわち、
気温が上昇する夏場はクーラーなどの冷房需要のために
、気温の上昇と共に電力需要量が急激に伸びる傾向にあ
る。一方、気温の変化が少ない初春から初夏、あるいは
初秋から初冬にかけた気温の変化がそれほど大きくない
時期は電力需要はほぼ一定である。また、冬季の気温が
低い季節は暖房需要が伸びるために、気温の低下と共に
電力需要が伸びる傾向にある。このような日最高気温、
日最低気温の関係を基に関数関係を導出し、不連続性影
響要因ごと又はその組合せごとに関数テーブル(I)4
に格納するようにしている。
FIG. 3 shows an example of a change in power demand over a certain day, and FIG. 4 shows an example of the relationship between power demand, daily maximum temperature, and daily minimum temperature at a certain time. In FIG. 3, the vertical axis indicates power demand and the horizontal axis indicates time. When predicting the power demand that changes in this way hourly, the function deriving means (I) 3 derives 24 functions for each time. It is also possible to calculate functions for each year, day, minute, second, etc. On the other hand, as is clear from the relationship curve H between daily maximum temperature and electricity demand at a certain time and the relationship curve L between daily minimum temperature and electricity demand at a certain time shown in Figure 4, there is a relationship between temperature and electricity demand. A kind of nonlinear relationship can be observed. Note that the vertical axis in the figure represents the power demand, and the horizontal axis represents the temperature. The existence of such a trend is considered to be due to the following reasons. That is,
In summer, when temperatures rise, electricity demand tends to increase rapidly due to the demand for air conditioners such as air conditioners. On the other hand, electricity demand remains almost constant during periods when temperature changes are small, such as from early spring to early summer, or from early autumn to early winter, when temperature changes are not large. In addition, demand for heating increases during the winter season when temperatures are low, so electricity demand tends to increase as the temperature decreases. The highest temperature on such a day,
A functional relationship is derived based on the relationship between daily minimum temperatures, and a function table (I) 4 is created for each discontinuity influencing factor or for each combination thereof.
I am trying to store it in .

【0030】関数導出手段(I)3の関数導出に係る計
算モデルとしては、周知の重回帰モデル又はニューラル
ネットを用いたモデルが適用できる。ここでは、ニュー
ラルネットを用いた具体例を図5に示して説明する。こ
れは、一般に対象が非線形である場合、関数関係を求め
る上でニューラルネットを用いる効用が大きいためであ
る。図5に示すように、ニューラルネットとは、ニュー
ロンと呼ばれる素子51を結合線52a,52b...
を用いて、入力層、1又は2以上の中間層、出力層から
なる階層状に連結し、人間の脳細胞が学習する動作をシ
ミューレーションするモデルのことである。この場合学
習を行うためのデータ、例えば入力である影響要因、例
えば日最高気温、日最低気温がそれぞれα(℃)、β(
℃)の時に、出力である電力需要量がγ(MW)である
というデータを教師データと呼ぶ。また、ニューラルネ
ットでの学習とは教師データである気温データが入力さ
れた時に、需要量である出力が対応した結果になるよう
に、結合線52a,52b,...に対する重み係数を
求めることである。重み係数を決定する方法は様々あり
、本実施例では以下にバックプロパゲーション法につい
て示すが、これに限る必要はない。
As a calculation model for function derivation by the function derivation means (I) 3, a well-known multiple regression model or a model using a neural network can be applied. Here, a specific example using a neural network will be described with reference to FIG. This is because, in general, when the object is nonlinear, the utility of using a neural network in finding functional relationships is great. As shown in FIG. 5, a neural network connects elements 51 called neurons to connecting lines 52a, 52b, . .. ..
This is a model that simulates the behavior learned by human brain cells by connecting an input layer, one or more intermediate layers, and an output layer in a hierarchical manner. In this case, data for learning, such as influence factors that are input, such as daily maximum temperature and daily minimum temperature, are α (℃) and β (
℃), the data that the power demand, which is the output, is γ (MW) is called teaching data. Learning in a neural network means that when temperature data, which is teacher data, is input, connection lines 52a, 52b, . .. .. The purpose is to find the weighting coefficient for . There are various methods for determining the weighting coefficients, and in this embodiment, a backpropagation method will be described below, but the method is not limited to this.

【0031】一例として、簡単のために3層のニューラ
ルネットに対するバックプロパゲーション法について説
明する。入力層ニューロンからの出力値をXk、中間層
ニューロン殻の出力値をXj、出力層ニューロンからの
出力値をXiとする。これらの中間層と出力層の出力値
Xi、Xjは以下の数式1で求める。
As an example, for the sake of simplicity, a backpropagation method for a three-layer neural network will be described. Let the output value from the input layer neuron be Xk, the output value of the intermediate layer neuron shell be Xj, and the output value from the output layer neuron be Xi. The output values Xi and Xj of these intermediate layer and output layer are determined by the following equation 1.

【0032】[0032]

【数1】Xi= f(ΣWij・Xj−θi)xj =
 f(ΣWjk・Xk−θj)ただし Wpq:結合線52a,52b,...に対する重み係
数 p,q:i,j,k i,j,k:1,2,・・・ θp:ニューロンが発火を始めるしきい値f  :シグ
モイド関数 f(n)=(1/(1+exp(−n))いま、ある教
師データを入力した場合、出力Xiが期待するXi0と
異なり、ずれていたとする。学習が終了することは以下
の不等式が成立することと等価である。
[Math. 1] Xi= f(ΣWij・Xj−θi)xj=
f(ΣWjk・Xk−θj) where Wpq: bond lines 52a, 52b, . .. .. Weighting coefficient p, q: i, j, k i, j, k: 1, 2,... θp: Threshold value f at which the neuron starts firing: Sigmoid function f(n)=(1/(1+exp( -n)) Now, suppose that when certain teacher data is input, the output Xi is different from the expected Xi0 and deviates. Completion of learning is equivalent to the following inequality being satisfied.

【0033】[0033]

【数2】Σ(Xi−Xi0)2<ε ただし Σ:総和を表す記号 ε:予め設定した誤差のしきい値 そこで、この数式2を満たすために、Wijを数式2の
左辺の誤差に比例して逆修正する。このようなニューラ
ルネットを用いて関数導出手段(I)3が構成されてい
る。
[Equation 2] Σ(Xi-Xi0)2<ε Where Σ: Symbol representing the sum ε: Preset error threshold Therefore, in order to satisfy Equation 2, Wij is proportional to the error on the left side of Equation 2. and reversely correct it. The function deriving means (I) 3 is configured using such a neural network.

【0034】ここで、上記のように構成される実施例に
ついて、基本構成部分の更に詳細な構成と動作を図6乃
至図8を用いて説明する。また、合わせて付加部分の詳
細構成を動作とともに説明する。関数導出手段(I)3
の処理手順は図6に示すようになっている。まず、ステ
ップ61において予測対象条件の内の不連続性影響要因
の組合せについての全ての関数を求めるため、組合せを
順次設定・変更する。ステップ62では設定・変更され
た不連続性影響要因の組合せに対応する関数を初期化す
る。すなわち意味のない関数にする。次に、ステップ6
3にて、設定された組合せに対応する実績データをデー
タテーブル2から抽出する。例えば、最近の一定期間に
おける日実績データのうちから、天候が「晴」、曜日が
「平日」の実績データを読み出す。この読み出したデー
タに基づいて、ステップ64にて各時刻ごとの日最高気
温又は日最低気温と電力需要量の相関関数をそれぞれ導
出する。このようにして導出された関数は、関数テーブ
ル(I)4に格納される。そして、ステップ65におい
て全ての組合せに係る関数の導出が終了したか否かを判
断し、終了していない場合はステップ61に戻って組合
せを変更し、全ての関数の導出を終了するまで、上記の
処理を繰り返す。
Now, regarding the embodiment configured as described above, the more detailed configuration and operation of the basic components will be explained using FIGS. 6 to 8. Additionally, the detailed configuration of the additional portion will be explained along with its operation. Function derivation means (I) 3
The processing procedure is shown in FIG. First, in step 61, in order to obtain all functions for combinations of discontinuity influencing factors among prediction target conditions, combinations are sequentially set and changed. In step 62, a function corresponding to the set/changed combination of discontinuity influencing factors is initialized. In other words, make it a meaningless function. Next, step 6
3, performance data corresponding to the set combination is extracted from the data table 2. For example, from among daily performance data for a recent certain period, performance data for which the weather is "sunny" and the day of the week is "weekdays" is read out. Based on this read data, in step 64, a correlation function between daily maximum temperature or daily minimum temperature and power demand is derived for each time. The functions derived in this way are stored in the function table (I) 4. Then, in step 65, it is determined whether or not the derivation of functions related to all combinations has been completed. If not, the process returns to step 61 to change the combination, and the above-mentioned procedure is continued until derivation of all functions is completed. Repeat the process.

【0035】ここで、上記ステップ64における処理を
具体的に説明する。いま、影響要因と電力需要量との相
関関数を数式3で表し、数式4の条件を満たすようにニ
ューラルネットの重み係数を設定して関数を導出する。
[0035] Here, the processing in step 64 will be specifically explained. Now, the correlation function between the influencing factors and the power demand is expressed by Equation 3, and the weighting coefficient of the neural network is set so as to satisfy the condition of Equation 4 to derive the function.

【0036】[0036]

【数3】 P(i,h)=g(X1,X2,...,Xn)[Math 3] P(i,h)=g(X1,X2,...,Xn)

【00
37】
00
37]

【数4】|P(i,h)−P’(i,h)|<εただし
、記号の意味は以下のとおりである。
[Equation 4] |P(i,h)-P'(i,h)|<ε However, the meanings of the symbols are as follows.

【0038】P(i,h):i日の時時刻hにおける電
力需要量(MW) P’(i,h):i日の時刻hにおけるニューラルネッ
トの出力 g:関数関係を示す記号 Xj:気候、社会的イベントなどを表す影響要因ただし
(j=1,2,...) ε:出力誤差のしきい値 ニューラルネットによる関数導出の信頼性を高めるには
、一般に行われているように、教師データによる学習を
繰返し行わせるのが望ましい。例えば、1つの日最高気
温と需要量の組合せについて学習を20回繰返し、更に
組合せを変えて全体として100回繰り返す。これは、
前の学習結果の影響が残っているので、繰返しによりそ
の影響を除去するのである。このようにして、連続性の
影響要因に係る需要量の予測関数が高い精度で得られる
。なお、不連続性影響要因を関数化しないのは、例えば
天候のように同じ晴であっても、曜日の種類によって需
要量の傾向がまったく違う場合があり、関数化の精度が
上がらないからである。
P(i,h): Power demand (MW) at time h on day i P'(i,h): Output of neural network at time h on day i g: Symbol indicating functional relationship Xj: Influencing factors representing climate, social events, etc. (j = 1, 2,...) ε: Output error threshold In order to increase the reliability of function derivation using a neural network, as is generally done, , it is desirable to repeatedly perform learning using teacher data. For example, learning is repeated 20 times for one combination of daily maximum temperature and demand amount, and then the combination is changed and repeated 100 times in total. this is,
Since the influence of previous learning results remains, this influence is removed through repetition. In this way, a demand quantity prediction function related to continuity influencing factors can be obtained with high accuracy. The reason why discontinuity influencing factors are not converted into functions is because, for example, even if the weather is the same, the trend in demand may be completely different depending on the type of day of the week, and the accuracy of function conversion cannot be improved. be.

【0039】図7に、予測実行手段7の処理手順を示す
。図示のように、ステップ71にて予測指令手段6から
予測対象の予測条件を入力する。予測条件として、例え
ば何月何日の電力需要量の日変化を予測するとし、その
日の予想天候は晴、曜日は平日、その日の予想日最高気
温と日最低気温等を入力する。次にステップ72で予測
対象時刻の初期値を設定する。次のステップ73にて、
後で説明するが、予測対象の条件が、実績データの存在
する領域内か領域外かの判断、つまり内挿か否かを判断
する。内挿の場合は、ステップ74に進んで関数テーブ
ル3から数式3に示した対応する関数のニューラルネッ
トを読みだし、その入力層に予測条件を入力し、その出
力を当該時刻の電力需要量の予測値とする。続いてステ
ップ75において、全ての時刻について予測を完了した
か否かを判断し、終了するまでステップ72に戻って上
記の予測処理を繰り返す。全ての予測が完了したときは
、ステップ76に進んで、予測結果を出力する。
FIG. 7 shows the processing procedure of the prediction execution means 7. As shown in the figure, in step 71, the prediction conditions of the prediction target are input from the prediction command means 6. As a prediction condition, for example, assume that the daily change in power demand is to be predicted on what month and day, and input the expected weather for that day as sunny, the day of the week as a weekday, the expected daily maximum temperature and daily minimum temperature for that day, etc. Next, in step 72, an initial value of the prediction target time is set. In the next step 73,
As will be explained later, it is determined whether the condition of the prediction target is within or outside the area in which actual performance data exists, that is, it is determined whether or not interpolation is required. In the case of interpolation, proceed to step 74, read out the neural network of the corresponding function shown in Equation 3 from function table 3, input the prediction condition to its input layer, and calculate the output of the power demand at the relevant time. Use as predicted value. Subsequently, in step 75, it is determined whether prediction has been completed for all times, and the process returns to step 72 to repeat the above prediction process until the prediction is completed. When all predictions are completed, the process proceeds to step 76 and the prediction results are output.

【0040】関数更新手段(I)5は、図8に示す手順
により、関数テーブル(I)4内の関数を更新して、新
たな実績データを反映する。上記の予測結果に対応する
時点の実績データが入力手段1を介してデータベース2
に随時入力される。関数更新手段(I)5は図8のステ
ップ81にて新たな実績データが入力されたか否かを判
断し、肯定のときはステップ82において、その新たな
データを加えて対応する関数を導出し、関数テーブル(
I)4を更新する。つまり、ニューラルネットの重み係
数を変更する。このステップ82の処理内容は、図6の
ステップ62,63,64,66,67に相当するので
、詳細は省略する。以下、新たな実績データが入力され
る度に、関数の更新が繰返しなされるが、この繰り返し
の期間が予め設定した期間(学習追加期間)を越えたか
否かを判断し(ステップ83)、肯定のときはステップ
84の進んで当該関数関係の初期化を行い、新たに初期
学習期間の影響要因と電力需要量の関数について、各時
点のニューラルネットを用いて求めて、関数を新たなも
のに設定しなおす。
The function updating means (I) 5 updates the functions in the function table (I) 4 to reflect new performance data according to the procedure shown in FIG. Actual data at the time corresponding to the above prediction result is sent to the database 2 via the input means 1.
is input from time to time. The function updating means (I) 5 determines whether new performance data has been input in step 81 of FIG. , function table (
I) Update 4. In other words, the weighting coefficients of the neural network are changed. The processing content of this step 82 corresponds to steps 62, 63, 64, 66, and 67 in FIG. 6, so the details will be omitted. Thereafter, each time new performance data is input, the function is repeatedly updated, but it is determined whether or not the period of this repetition exceeds a preset period (additional learning period) (step 83). In this case, proceed to step 84 to initialize the relevant functional relationship, and use the neural network at each point in time to newly find the influencing factors of the initial learning period and the function of power demand, and create a new function. Reset.

【0041】上記の手順により、電力需要予測を行った
結果の一例を次に説明する。 ・連続性影響要因の条件:日最高気温、日最低気温・不
連続性影響要因の条件:天候−晴、曜日−平日・初期学
習期間:1カ月 ・学習追加期間:1カ月 図9にある月の予測を行った結果を示す。○はニューラ
ルネットにおいて、教師データとなる予測実行1カ月前
の実績量であり、●は本需要予測方法により求めた予測
量、*は●に対応する実績量である。影響要因である気
温が○印の存在範囲中にある場合の予測量は実績量とほ
とんど一致し、精度のよい予測を行うことができている
(内挿のケース)。しかし、それ以外の範囲にある予測
量は一致の度合が非常に低く、精度が悪くなる(外挿の
ケース)。
An example of the results of power demand prediction performed using the above procedure will be explained next.・Continuity influencing factor conditions: daily maximum temperature, daily minimum temperature ・Discontinuity influencing factor conditions: Weather - sunny, day of the week - weekday ・Initial learning period: 1 month ・Additional learning period: 1 month Months in Figure 9 The results of the prediction are shown below. In the neural network, ○ is the actual amount one month before the prediction execution which becomes the training data, ● is the predicted amount obtained by this demand forecasting method, and * is the actual amount corresponding to ●. When the temperature, which is an influencing factor, is within the range marked with a circle, the predicted amount almost matches the actual amount, indicating that highly accurate predictions can be made (interpolation case). However, predicted quantities in other ranges have a very low degree of agreement and poor accuracy (extrapolation case).

【0042】これらのことに基づいて、需要予測の精度
評価を内挿のケースのみ、外挿のケースと内挿のケース
をひとまとめにしたものの2通りについて評価を行う。 需要予測の精度評価の基準は大きく2つに分類され、翌
日予測の場合次の2点である。
Based on these, the accuracy of demand forecasting is evaluated in two ways: only the interpolation case and the extrapolation case and the interpolation case together. Criteria for evaluating the accuracy of demand forecasting are broadly classified into two categories, and in the case of next-day forecasting, the following two points are used.

【0043】・日平均誤差を減少すること・日平均誤差
のばらつきを減少することここで、日平均誤差は次式5
の精度評価式で行う。
・Reducing the daily average error ・Reducing the variation in the daily average error Here, the daily average error is calculated by the following formula 5.
This is done using the accuracy evaluation formula.

【0044】[0044]

【数5】[Math 5]

【0045】ただし、    abs:絶対値を表す記
号Σ:総和を表す記号 i=1,2,...,24 上記式は、各時点iでの需要予測値と需要実績値との差
の絶対値を実績値で除することによって相対誤差を求め
、各時点毎の相対誤差の総和を時点数で除することによ
って日平均の誤差を求めるものである。
[0045] However, abs: symbol representing absolute value Σ: symbol representing sum i=1, 2, . .. .. , 24 The above formula calculates the relative error by dividing the absolute value of the difference between the predicted demand value and the actual demand value at each time point i by the actual value, and then divides the sum of the relative errors at each time point by the number of time points. By doing this, the daily average error is determined.

【0046】このような誤差の評価法としては、種々の
変形が考えられる。また、本実施例では相対誤差を用い
て評価を行ったが、予測値と需要値の絶対値を評価基準
に用いることも勿論可能である。
Various modifications can be considered as a method for evaluating such an error. Moreover, although the evaluation was performed using relative errors in this embodiment, it is of course possible to use the absolute values of the predicted value and the demand value as evaluation criteria.

【0047】誤差の一例として、最近1年5カ月分の実
績データに基づいた予測結果を、図10に示す。図に示
すように、日平均誤差の平均値が内挿のケースだけの場
合は2.18%であり、外挿のケースを含んだ場合は2
.84%と、外挿のケースを含んだ場合は内挿のケース
だけの場合よりも精度が悪くなることが観測できる。
As an example of the error, a prediction result based on the past 1 year and 5 months' worth of actual data is shown in FIG. As shown in the figure, the average value of the daily average error is 2.18% when only the interpolation case is involved, and 2.18% when the extrapolation case is included.
.. It can be observed that the accuracy is 84%, which is worse when including the extrapolation case than when only the interpolation case is included.

【0048】ここで、図1実施例の付加部分11〜14
の詳細について説明する。この付加部分は予測精度の向
上のために付加されたものであり、特に疑似データ作成
手段11を設けたことに特徴がある。この疑似データ作
成手段11を設けた理由について説明する。精度のよい
需要予測を行うために、外挿のケースのもとで需要予測
を行うことは望ましくないことを前記した。そこで、外
挿のケースを減少するために、もともと外挿のケースに
当たる予想時点の影響要因を、内挿のケースに変化させ
るため、疑似的な教師データを関数導出手段(II)1
2に入力するようにしたのである。疑似データの一例を
図11を用いて説明する。図11は春先のある月の実績
データをもとに、翌月におけるそれぞれの日々の電力需
要量を求める場合の例である。図中、黒で塗りつぶした
楕円はある月の実績データの範囲、ハッチングを施した
楕円は予測対象であるその翌月の実績データの範囲であ
る。また縦軸は電力需要量、横軸は気温を示す。学習デ
ータとなる前記実績データの存在する気温の範囲中には
予想実行月の実績データが存在しないので予測の実行に
際し、これは外挿のケースに当たり、前述したように予
測の精度が悪くなる。一般に、そのような外挿のケース
を解消するために、教師データをクリアせずに蓄積する
方法がある。しかし、この方法には以下にあげる2つの
欠点がある。
Here, the additional parts 11 to 14 of the embodiment in FIG.
The details will be explained below. This additional part was added to improve prediction accuracy, and is particularly characterized by the provision of pseudo data creation means 11. The reason for providing this pseudo data creation means 11 will be explained. As mentioned above, in order to perform accurate demand forecasting, it is not desirable to perform demand forecasting based on the extrapolation case. Therefore, in order to reduce the number of cases of extrapolation, in order to change the influencing factors at the predicted time point, which originally corresponds to the case of extrapolation, to the case of interpolation, pseudo training data is used as function derivation means (II) 1.
2. An example of pseudo data will be explained using FIG. 11. FIG. 11 is an example in which the daily power demand for the next month is calculated based on actual data for a certain month in early spring. In the figure, the black ellipse is the range of actual data for a certain month, and the hatched ellipse is the range of actual data for the next month, which is the target of prediction. Also, the vertical axis shows the amount of electricity demand, and the horizontal axis shows the temperature. Since there is no actual data for the month in which the prediction is executed within the temperature range in which the actual data serving as learning data exists, this corresponds to a case of extrapolation when executing the prediction, and as described above, the accuracy of the prediction deteriorates. Generally, in order to resolve such cases of extrapolation, there is a method of accumulating training data without clearing it. However, this method has the following two drawbacks.

【0049】1)関数関係の導出に多くの時間を要する
2)過去のデータの影響を大きく受け、予測精度が低下
する このような問題点を解決するために、本実施例では、疑
似的な白抜きの楕円の範囲で表されるデータを付け加え
、外挿のケースを減少する方法を導入している。疑似デ
ータの作成は、計算機等を用いて原則として自動的にお
こなう。しかし、必要に応じて、人間が持つノウハウを
オペレータが手動で入力することによる疑似データ生成
も可能である。
1) It takes a lot of time to derive the functional relationship. 2) The prediction accuracy is reduced due to the large influence of past data. In order to solve these problems, in this example, a pseudo A method is introduced to reduce the number of cases of extrapolation by adding data represented by the range of the white ellipse. In principle, pseudo data is created automatically using a computer, etc. However, if necessary, it is also possible to generate pseudo data by manually inputting human know-how by an operator.

【0050】疑似データ作成手段11における疑似デー
タの作成方法としては様々な方法が考えられる。本実施
例ではその中の代表的な一例を以下に示す。 いま、Tmin(i):予想月の前年同月のデータの日
最低気温 Tmax(i):予想月の前年同月のデータの日最低気
温P(i,h):予想月の前年同月のデータの時刻hに
おける電力需要量 とする。
Various methods can be considered for creating pseudo data in the pseudo data creating means 11. In this embodiment, a typical example thereof is shown below. Now, Tmin(i): Daily minimum temperature of data for the same month of the previous year of the forecast month Tmax(i): Daily minimum temperature of data of the same month of the previous year of the forecast month P(i, h): Time of data of the same month of the previous year of the forecast month Let it be the power demand at h.

【0051】まず、1年前の同時期の同時点のデータ(
Tmin(i),Tmax(i),P(i,h)))か
ら、(Tmin(i),Tmax(i),P”(i,h
))を求める。ただし、
First, data at the same time one year ago (
Tmin(i), Tmax(i), P(i, h))), then (Tmin(i), Tmax(i), P"(i, h
)). however,

【0052】[0052]

【数6】[Math 6]

【0053】Pmean(h):予想月の前月の時刻h
における電力需要量の平均値 P’mean(h):予想月の前年同月の時刻hに於け
る電力需要量の平均値 なる形式で疑似データを作成する。上式で、Pmean
(h)/P’mean(h)の代わりに、年負荷成長率
などを用いることも可能である。また、影響要因の一つ
である電気製品の売上高の前年比、又は社会的イベント
などの要因もこの部位に加味することが可能である。
Pmean(h): Time h of the previous month of the predicted month
P'mean(h): Pseudo data is created in the form of the average value of power demand at time h in the same month of the previous year of the forecast month. In the above formula, Pmean
Instead of (h)/P'mean(h), it is also possible to use the annual load growth rate or the like. In addition, factors such as the year-on-year change in sales of electrical products, which is one of the influencing factors, or social events can also be taken into consideration in this section.

【0054】また、関数導出手段(II)12で上記の
ような疑似データを用いて非線形関数を導出する場合、
疑似データの信頼性は100%でないから、本実施例で
はデータの信頼性に応じて関数の導出法を変更させるよ
うにしている。具体的には、ニューラルネットを関数の
モデルとして用いた場合は、ニューラルネットの重み係
数を緩い収束条件で求める一方、実績データが存在する
場合は、ニューラルネットの重み係数をきびしい収束条
件で求める方式とする。つまり、ニューラルネットの出
力値と、疑似データ出力値との間の誤差を、実績データ
を用いた場合よりも、高めに設定するような関数導出法
を用いる。すなわち、数式7の関数型に対し、数式8の
条件を満たすように、ニューラルネットの重み係数を更
新する。
Furthermore, when the function deriving means (II) 12 derives a nonlinear function using the above-mentioned pseudo data,
Since the reliability of pseudo data is not 100%, in this embodiment, the method for deriving the function is changed depending on the reliability of the data. Specifically, when a neural network is used as a function model, the weighting coefficients of the neural network are determined under loose convergence conditions, while when actual data exists, the weighting coefficients of the neural network are determined under strict convergence conditions. shall be. In other words, a function derivation method is used in which the error between the output value of the neural network and the pseudo data output value is set higher than when actual data is used. That is, the weighting coefficient of the neural network is updated so that the function type of Equation 7 satisfies the condition of Equation 8.

【0055】[0055]

【数7】 P”(i,h)=g(X1,X2,...,Xn)[Math 7] P”(i,h)=g(X1,X2,...,Xn)

【0
056】
0
056]

【数8】|P”(i,h)−P’(i,h)|<ε’た
だし、記号の意味は以下のとおりである。
[Equation 8] |P''(i, h)-P'(i, h)|<ε' However, the meanings of the symbols are as follows.

【0057】P”(i,h):i日の時時刻hにおける
電力需要量の疑似データ P’(i,h):i日の時刻hにおけるニューラルネッ
トの出力 g:関数関係を示す記号 Xj:気候、社会的イベントなどを表す影響要因ただし
(j=1,2,...) ε’:疑似データに対する出力誤差のしきい値(ε<ε
’) このようにして生成された疑似データを実績データに加
え、関数導出手段(II)12で関数IIを求める。こ
の関数導出手段(II)12の構成は関数導出手段(I
)3と同様であるから、説明を省略する。また、疑似デ
ータ作成と関数導出の処理手順は、図6のステップ66
,67に挿入されており、関数Iと関数IIは並行して
求められる。
P''(i,h): Pseudo data of power demand at time h on day i P'(i,h): Output of neural network at time h on day i g: Symbol Xj indicating functional relationship : Influencing factors representing climate, social events, etc., where (j=1, 2,...) ε': Threshold of output error for pseudo data (ε<ε
') The pseudo data thus generated is added to the actual data, and the function deriving means (II) 12 calculates the function II. The structure of the function deriving means (II) 12 is as follows:
) Since it is the same as 3, the explanation will be omitted. Furthermore, the processing procedure for creating pseudo data and deriving functions is shown in step 66 in FIG.
, 67, and Function I and Function II are calculated in parallel.

【0058】なお、図6の処理手順に代え、予測実行手
段(I)5における予測対象の条件が関数Iの基礎とな
る実績データの範囲から外れるか否かの判断結果(ステ
ップ73)を、疑似データ作成手段11に出力するよう
にし、これに応動させて疑似データを作成し、更に関数
導出手段(II)12を起動し関数IIを導出させるよ
うにしてもよい。
Note that instead of the processing procedure shown in FIG. 6, the judgment result (step 73) of whether or not the condition of the prediction target in the prediction execution means (I) 5 falls outside the range of the actual performance data that is the basis of the function I is The pseudo data may be output to the pseudo data creation means 11, pseudo data may be created in response to this, and the function derivation means (II) 12 may be activated to derive the function II.

【0059】また、図1、図6,7の実施例では、関数
導出手段(I)3,12と予測実行手段7の処理は直接
関連性がなく、それぞれ独立に処理を実行する場合を例
に説明した。しかし、本発明はこれに限られるものでは
なく、予測指令入力手段6から入力される予測指令及び
条件を関数導出手段(I)3,12に入力し、その予測
条件に合った実績データのみをデータベース2から抽出
するとともに、必要に応じて疑似データを作成するよう
にし、予測指令の入力に応動してかつ予測対象に関係す
る関数だけを導出するようにすることができる。これに
よれば、不必要な関数の導出を省略できる。
Furthermore, in the embodiments shown in FIGS. 1, 6, and 7, the processing of the function derivation means (I) 3, 12 and the prediction execution means 7 are not directly related, and the processing is executed independently. explained. However, the present invention is not limited to this, and the prediction commands and conditions inputted from the prediction command input means 6 are input to the function derivation means (I) 3, 12, and only the actual data matching the prediction conditions are input. In addition to extracting from the database 2, pseudo data can be created as needed, and only functions related to the prediction target can be derived in response to the input of the prediction command. According to this, derivation of unnecessary functions can be omitted.

【0060】以下、更に上記実施例の変形例を説明する
。図12を参照して、更に予測精度を向上させる実施例
について説明する。この例は、図12(a)に示すよう
に、電力需要量の日変化はある一定の最低需要電力量B
をベースとしてそれよりも低くならないことに鑑み、そ
のベース除去分Bを差し引いた需要量(図12(b))
を予測対象とすることにより、予測感度を高くし、これ
により予測精度を向上させようとするものである。すな
わち、前述した予測の実績データ(Tmin(i),T
max(i),P(i,h))を、(Tmin(i),
Tmax(i),P(i,h)−B)に変更して、関数
導出手段4,12を作動させるようにする。ここで、 Tmin(i):i番目の教師データの日最低気温Tm
ax(i):i番目の教師データの日最高気温P(i,
h)  :i番目の教師データの時刻hにおける電力需
要量 B  :ベース除去分 また、図13及び14を参照して、上記実施例で導出し
た関数を修正することにより、予測の精度を向上させる
実施例について説明する。図13は、関数の修正の概念
を便宜的に二次元平面で表現したものである。図13及
び14に示すように、関数導出手段(I)3,12にお
いて、関数を導出した結果に基づき、導出した関数と実
績データを比較し、実績データの中から需要傾向に合わ
ない不適なデータを予め設定された基準に基づいて検出
し(ステップ101,111)、それらの不適なデータ
を棄却した適正なデータのみを用いて再度関数を導出す
るようにしたものである(ステップ102,112)。 この方法は関数I,IIのいずれの導出にも適用できる
。なお、データ棄却の具体的方法としては、関数再導出
の際に、ニューラルネットの出力誤差のしきい値を緩や
かに調整することにより、実質的に棄却することができ
る。
Modifications of the above embodiment will be further explained below. With reference to FIG. 12, an example for further improving prediction accuracy will be described. In this example, as shown in FIG. 12(a), the daily change in power demand is a certain minimum power demand B
Considering that it will not be lower than the base value, the demand quantity after subtracting the base removal amount B (Figure 12(b))
By targeting the prediction target, the prediction sensitivity is increased and the prediction accuracy is thereby improved. That is, the actual prediction data (Tmin(i), T
max(i), P(i, h)), (Tmin(i),
Tmax(i), P(i,h)-B), and the function deriving means 4 and 12 are activated. Here, Tmin(i): Daily minimum temperature Tm of the i-th teacher data
ax(i): Daily maximum temperature P(i,
h): Electricity demand amount at time h of the i-th training data B: Base removal amount Also, with reference to FIGS. 13 and 14, the accuracy of prediction is improved by correcting the function derived in the above example. An example will be explained. FIG. 13 conveniently expresses the concept of function modification on a two-dimensional plane. As shown in FIGS. 13 and 14, the function derivation means (I) 3 and 12 compare the derived function with actual data based on the result of deriving the function, and identify inappropriate functions that do not match the demand trend from among the actual data. Data is detected based on preset criteria (steps 101 and 111), and a function is derived again using only proper data by rejecting the inappropriate data (steps 102 and 112). ). This method can be applied to the derivation of either function I or II. As a specific method for rejecting data, it is possible to substantially reject the data by gently adjusting the threshold value of the output error of the neural network when re-deriving the function.

【0061】図15及び16を参照して、更に予測精度
を向上させるために、関数を補正するようにした実施例
を説明する。図15は本実施例の概念を説明するもので
あり、次に定義する誤差Eを補正値として、導出した関
数を補正するようにしたものである(図16のステップ
103,104)。
Referring to FIGS. 15 and 16, an embodiment will be described in which the function is corrected in order to further improve the prediction accuracy. FIG. 15 explains the concept of this embodiment, in which the derived function is corrected using the error E defined next as a correction value (steps 103 and 104 in FIG. 16).

【0062】[0062]

【数9】[Math. 9]

【0063】ここで、g(X1,...Xn):導出し
た関数 X1,...,Xn:影響要因 Yi:電力需要量 m  :実績データの数 n  :影響要因の数 上式9で求めた補正値Eにより、図15に示すように、
補正前の関数gをEだけ平行移動して補正後の関数g’
を求める。これにより、予測精度を向上できる。
Here, g(X1,...Xn): derived functions X1,... .. .. , Xn: Influencing factor Yi: Power demand m: Number of actual data n: Number of influencing factors By the correction value E obtained by the above formula 9, as shown in FIG. 15,
The function g before correction is translated by E to obtain the function g' after correction.
seek. Thereby, prediction accuracy can be improved.

【0064】また、上記実施例では考慮していないが、
需要予測を行う際、特異日、例えば飛び石連休の中日、
お盆明け最初の平日などの予測は予測値と実績値との誤
差が大きいことが多い。この点に対処するために、過去
の特異日の予測値と実績値との誤差をデータベース2に
保存、または実績データ入力手段1を通じて入力し、こ
のデータを用いて関数導出手段(I)3,12で求めた
関数I,IIを補正することにより、精度のよい予測を
実行することが可能である。
[0064] Although not taken into consideration in the above embodiment,
When performing demand forecasting, special days, such as the middle of the Topping Stones holiday,
Forecasts for the first weekday after the Obon holidays often have large errors between predicted values and actual values. In order to deal with this point, the error between the predicted value and the actual value of the past singular day is stored in the database 2 or inputted through the actual data input means 1, and this data is used to calculate the function derivation means (I) 3, By correcting the functions I and II obtained in step 12, it is possible to perform accurate prediction.

【0065】以上、本発明を電力需要量の予測に適用し
た実施例に基づいて説明したが、前述したように、本発
明は電力に限らず、上水需要量やパンなどのような季節
や気温などの影響要因により需要量が増減する対象に適
用できることは明らかである。
[0065] The present invention has been described above based on an embodiment in which the present invention is applied to predicting the amount of electricity demand. However, as mentioned above, the present invention is not limited to electricity, but can also be applied to predictions of water demand, bread demand, etc. It is clear that this method can be applied to objects whose demand increases or decreases depending on influencing factors such as temperature.

【0066】[0066]

【発明の効果】以上説明したように、本発明によれば、
次の効果がある。不連続性影響要因ごとに又はその組合
せごとに、連続性影響要因と需要量の相関関数を求め、
該関数に基づいて予測対象時点の需要量を予測すること
から、曜日の違いや社会的なイベント等が行われる特異
日についても実績データに基づいた関数により予測でき
るので、予測精度が向上する。
[Effects of the Invention] As explained above, according to the present invention,
It has the following effects. For each discontinuity influencing factor or for each combination thereof, find the correlation function between the continuity influencing factor and the demand quantity,
Since the demand amount at the time of the prediction target is predicted based on the function, it is possible to predict differences in days of the week and special days when social events are held using the function based on actual data, thereby improving prediction accuracy.

【0067】また、連続性影響因子についての相関関数
を求めて予測することから、実測データの数が少なくて
も関数の信頼性が高いので、最近の一定期間の実績デー
タを用いることができ、過去の古い実績の影響を必要最
小限にできる。特に、関数のモデルを重回帰モデル又は
ニューラルネットを適用すれば、更に実績データの対象
期間を短縮できる。
[0067] Furthermore, since the correlation function for the continuity influencing factor is calculated and predicted, the reliability of the function is high even if the number of actual measurement data is small, so it is possible to use actual data for a recent certain period. The influence of past performance can be minimized to the necessary minimum. In particular, if a multiple regression model or a neural network is applied to the function model, the period covered by the performance data can be further shortened.

【0068】また、疑似データを作成するようにしたも
のによれば、実績データが少ない場合、又は予測対象時
点の連続性影響因子に対応する実績データが無い場合で
も、予測精度が向上する。
Furthermore, by creating pseudo data, the prediction accuracy is improved even when there is little actual data or when there is no actual data corresponding to the continuity influencing factor at the time of the prediction target.

【0069】また、新たな実績データを加えて関数を更
新するようにしたのものによれば、一層予測精度が向上
する。特に、更新に係る新たな実績データの期間が所定
の追加期間に達したとき、関数を初期化して導出しなお
せば、過去の古い実績データの影響を少なくできる。
Furthermore, if the function is updated by adding new performance data, the prediction accuracy is further improved. In particular, when the period of new performance data related to an update reaches a predetermined additional period, the influence of past old performance data can be reduced by initializing and re-deriving the function.

【0070】前記関数を求める手段として、重回帰モデ
ル又はニューラルネットモデルを用い、実績データ等の
信頼度に応じてモデルの収束条件を変えたり、疑似デー
タを用いる場合のモデルの収束条件を緩やかに設定する
ようにしたものによれば、妥当な予測を実現できる。
[0070] As a means of determining the function, a multiple regression model or a neural network model is used, and the convergence conditions of the model are changed depending on the reliability of actual data, etc., or the convergence conditions of the model are relaxed when pseudo data is used. Depending on the settings, reasonable predictions can be achieved.

【0071】また、実績データからベース除去分を除去
した実績値を用いたり、導出した関数から大きくはずれ
る実績データを棄却して再度関数を導出したり、導出し
た関数と実績データとのずれの程度に応じて関数を補正
するようにしたものによれば、更に予測精度を向上させ
ることができる。
[0071] In addition, it is possible to use actual values obtained by removing base removal from actual data, to reject actual data that deviates significantly from the derived function and to derive the function again, or to calculate the degree of deviation between the derived function and actual data. If the function is corrected in accordance with the above, the prediction accuracy can be further improved.

【図面の簡単な説明】[Brief explanation of drawings]

【図1】本発明の一実施例需要予測装置のブロック構成
図である。
FIG. 1 is a block diagram of a demand forecasting device according to an embodiment of the present invention.

【図2】データベースの一実施例のテーブル構成図であ
る。
FIG. 2 is a table configuration diagram of an embodiment of a database.

【図3】電力需要量の日変化の一例を示す図である。FIG. 3 is a diagram showing an example of daily changes in power demand.

【図4】気温と電力需要量の相関の一例を示す線図であ
る。
FIG. 4 is a diagram showing an example of the correlation between temperature and power demand.

【図5】関数導出モデルの一例のニューラルネットの構
成図である。
FIG. 5 is a configuration diagram of a neural network as an example of a function derivation model.

【図6】関数導出手段の処理手順のフローチャートであ
る。
FIG. 6 is a flowchart of the processing procedure of the function deriving means.

【図7】予測実行手段の処理手順のフローチャートであ
る。
FIG. 7 is a flowchart of the processing procedure of the prediction execution means.

【図8】関数更新手段の処理手順のフローチャートであ
る。
FIG. 8 is a flowchart of a processing procedure of a function updating means.

【図9】予測対象の影響要因が導出した関数の実績デー
タの範囲の内挿ケースと外挿ケースの場合の予測結果を
比較説明する図である。
FIG. 9 is a diagram for comparing and explaining the prediction results in the case of interpolation and extrapolation of the range of performance data of the function derived from the influencing factor of the prediction target.

【図10】内挿ケースと外挿ケースの予測結果の誤差を
比較した図である。
FIG. 10 is a diagram comparing errors in prediction results between the interpolation case and the extrapolation case.

【図11】疑似データの作成方法を説明するための図で
ある。
FIG. 11 is a diagram for explaining a method of creating pseudo data.

【図12】需要実績量からベース除去分を差し引いた値
を需要実績値として関数を導出する実施例の説明図であ
り、(a)はベース除去分を含んだ需要実績値、(b)
はベース除去分を差し引いた需要実績値の線図である。
FIG. 12 is an explanatory diagram of an example in which a function is derived using the actual demand value obtained by subtracting the base removal amount from the actual demand amount, (a) is the actual demand value including the base removal amount, (b)
is a diagram of the actual demand value after subtracting the base removal amount.

【図13】導出した関数から外れる実績データを棄却す
る実施例を説明するための線図である。
FIG. 13 is a diagram for explaining an example in which performance data that deviates from the derived function is rejected.

【図14】導出した関数から外れる実績データを棄却す
る実施例の処理手順の主要部フローチャートである。
FIG. 14 is a flowchart of the main part of the processing procedure of the embodiment for rejecting performance data that deviates from the derived function.

【図15】導出した関数から外れる実績データの程度に
応じて関数を補正する実施例を説明するための線図であ
る。
FIG. 15 is a diagram for explaining an embodiment in which a function is corrected according to the degree of performance data that deviates from the derived function.

【図16】導出した関数から外れる実績データの程度に
応じて関数を補正する実施例の処理手順の主要部フロー
チャートである。
FIG. 16 is a flowchart of the main part of the processing procedure of the embodiment for correcting a function according to the degree of performance data that deviates from the derived function.

【符号の説明】[Explanation of symbols]

1  実績データ入力手段 2  データベース 3  関数導出手段(I) 4  関数テーブル(I) 5  関数更新手段(I) 6  予測指令入力手段 7  予測実行手段 8  予測結果出力手段 11  疑似データ作成手段 12  関数導出手段(II) 13  関数テーブル(II) 14  関数更新手段(II) 1. Actual data input means 2 Database 3 Function derivation means (I) 4 Function table (I) 5 Function update means (I) 6 Prediction command input means 7 Prediction execution means 8 Prediction result output means 11 Pseudo data creation means 12 Function derivation means (II) 13 Function table (II) 14 Function update means (II)

Claims (23)

【特許請求の範囲】[Claims] 【請求項1】  需要量の実績データを需要量の増減に
影響を与える複数の影響要因を付して収集し、影響要因
の変化と需要量の変化との関係が連続的な連続性影響要
因とその関係が不連続的な不連続性影響要因とに別け、
前記実績データの最近の一定期間の実績データを用い、
不連続性影響要因ごとに又は不連続性影響要因の組合せ
ごとに前記連続性影響要因と需要量の相関関数を求め、
該関数の中から予測対象時点の予測影響要因に対応する
関数を用いて当該予測対象時点の需要量を予測すること
を含む需要量予測方法。
[Claim 1] Actual data on demand quantity is collected with a plurality of influencing factors that influence increases and decreases in demand quantity, and the relationship between changes in the influencing factors and changes in demand quantity is continuous influencing factors. and discontinuous influence factors whose relationship is discontinuous,
Using the performance data for a recent certain period of the performance data,
Determining the correlation function between the continuity influencing factor and the demand quantity for each discontinuity influencing factor or for each combination of discontinuity influencing factors,
A demand amount forecasting method comprising predicting a demand amount at a prediction target time using a function corresponding to a prediction influence factor at the prediction target time from among the functions.
【請求項2】  需要量の実績データを需要量の増減に
影響を与える複数の影響要因を付して収集し、影響要因
の変化と需要量の変化との関係が連続的な連続性影響要
因とその関係が不連続的な不連続性影響要因とに別け、
前記実績データの最近の一定期間の実績データから外れ
る領域について過去の実績データに基づいて予測対象時
点の疑似データを作成し、該疑似データと前記一定期間
の実績データとを用い、不連続性影響要因ごとに又は不
連続性影響要因の組合せごとに前記連続性影響要因と需
要量の関数を求め、該関数の中から予測対象時点の予測
影響要因に対応する関数を用いて当該予測対象時点の需
要量を予測することを含む需要量予測方法。
[Claim 2] Actual data on the amount of demand is collected with multiple influencing factors that influence increases and decreases in the amount of demand, and the relationship between changes in the influencing factors and changes in the amount of demand is a continuous influencing factor. and discontinuous influence factors whose relationship is discontinuous,
Create pseudo data at the time of prediction target based on past performance data for areas that deviate from the performance data for a recent certain period of the performance data, and use the pseudo data and the performance data for the certain period to calculate discontinuity effects. For each factor or for each combination of discontinuity influencing factors, a function of the continuity influencing factor and the demand quantity is determined, and from among the functions, the function corresponding to the predicted influencing factor at the forecasting point is used to calculate the function at the forecasting point. A demand forecasting method that involves forecasting demand.
【請求項3】  需要量の実績データを需要量の増減に
影響を与える複数の影響要因を付して収集し、影響要因
の変化と需要量の変化との関係が連続的な連続性影響要
因とその関係が不連続的な不連続性影響要因とに別け、
前記実績データの最近の一定期間の実績データを用い、
不連続性影響要因ごとに又は不連続性影響要因の組合せ
ごとに前記連続性影響要因と需要量の第1の関数を求め
るとともに、前記一定期間の実績データから外れる領域
について過去の実績データに基づいて予測対象時点の疑
似データを作成し、該疑似データと前記一定期間の実績
データを用い、不連続性影響要因ごとに又は不連続性影
響要因の組合せごとに前記連続性影響要因と需要量の第
2の関数を求め、第1と第2の関数の中から予測対象時
点の予測影響要因に対応する関数を用いて当該予測対象
時点の需要量を予測することを含む需要量予測方法。
[Claim 3] Actual data on the amount of demand is collected with a plurality of influencing factors that influence increases and decreases in the amount of demand, and the relationship between changes in the influencing factors and changes in the amount of demand is a continuous influencing factor. and discontinuous influence factors whose relationship is discontinuous,
Using the performance data for a recent certain period of the performance data,
Calculate the first function of the continuity influencing factor and demand quantity for each discontinuity influencing factor or for each combination of discontinuity influencing factors, and based on past actual data for areas that deviate from the actual data for the certain period of time. Create pseudo data at the time of the forecast target, and use the pseudo data and the actual data for the certain period to calculate the continuity impact factors and demand for each discontinuity impact factor or for each combination of discontinuity impact factors. A demand quantity forecasting method comprising calculating a second function and predicting a demand quantity at a prediction target time using a function corresponding to a prediction influence factor at the prediction target time from among the first and second functions.
【請求項4】  請求項1乃至3のいずれかにおいて、
前記予測需要量に対応する需要量の新たな実績データを
前記一定期間の実績データに加えて前記関数を更新する
ことを特徴とする需要量予測方法。
[Claim 4] In any one of claims 1 to 3,
A demand forecasting method, comprising: updating the function by adding new performance data of the demand corresponding to the predicted demand to the performance data for the certain period.
【請求項5】  請求項4において、前記更新に係る新
たな実績データの期間が所定の追加期間に達したとき前
記関数を初期化し、該新たな実績データを含む実績デー
タの最近の前記一定期間の実績データを用いて、前記関
数を新たに求めることを特徴とする需要量予測方法。
5. In claim 4, the function is initialized when the period of new performance data related to the update reaches a predetermined additional period, and the function is initialized when the period of new performance data related to the update reaches a predetermined additional period, and A method for predicting demand quantity, characterized in that the function is newly determined using actual data of.
【請求項6】  請求項1乃至3のいずれかにおいて、
前記関数導出に用いる実績データが実際の実績データか
らベース分を除去したデータであることを特徴とする需
要予測方法。
[Claim 6] In any one of claims 1 to 3,
A demand forecasting method characterized in that the performance data used for deriving the function is data obtained by removing a base portion from actual performance data.
【請求項7】  請求項1乃至3のいずれかにおいて、
前記関数を導出した後、該導出した関数から一定量外れ
る実績データを棄却し、残った実績データを用いて関数
を再度導出することを特徴とする需要予測方法。
[Claim 7] In any one of claims 1 to 3,
A demand forecasting method characterized in that, after deriving the function, performance data that deviates by a certain amount from the derived function is discarded, and the function is derived again using the remaining performance data.
【請求項8】  請求項1乃至3のいずれかにおいて、
前記関数を導出した後、該導出した関数と実績データと
を比較し、関数からの実績データの外れの程度に応じて
を当該関数を補正することを特徴とする需要量予測方法
[Claim 8] In any one of claims 1 to 3,
After deriving the function, the derived function is compared with actual data, and the function is corrected according to the degree of deviation of the actual data from the function.
【請求項9】  請求項1において、前記関数を、重回
帰モデル又はニューラルネットモデルを用いて求めるこ
とを特徴とする需要量予測方法。
9. The demand forecasting method according to claim 1, wherein the function is determined using a multiple regression model or a neural network model.
【請求項10】  請求項1乃至3において、前記関数
を、重回帰モデル又はニューラルネットモデルを用いて
求め、その際用いる実績データ等の信頼度に応じてモデ
ルの収束条件を変えることを特徴とする需要量予測方法
10. Claims 1 to 3, characterized in that the function is determined using a multiple regression model or a neural network model, and the convergence conditions of the model are changed depending on the reliability of actual data used at that time. demand forecasting method.
【請求項11】  請求項3において、前記第1と第2
の関数を重回帰モデル又はニューラルネットモデルを用
いて求め、前記第2の関数を求める際に前記実績データ
を用いる場合よりも前記疑似データを用いる場合のモデ
ルの収束条件を緩やかに設定することを特徴とする需要
量予測方法。
11. In claim 3, the first and second
The function of is determined using a multiple regression model or a neural network model, and when determining the second function, the convergence conditions of the model are set more gently when using the pseudo data than when using the actual data. Characteristic demand forecasting method.
【請求項12】  需要量の実績データが需要量の増減
に影響を与える複数の影響要因を付して格納されるデー
タベースと、影響要因の変化と需要量の変化との関係が
不連続的な予め定められた不連続性影響要因ごと又は不
連続性影響要因の組合せごとに対応する最近の一定期間
の実績データを前記データベースから順次抽出し、該抽
出した実績データに基づいて不連続性影響要因ごとに又
は不連続性影響要因の組合せごとに前記連続性影響要因
と需要量の相関関数を求める関数導出手段と、該求めら
れた関数が格納される関数テーブルと、予測対象時点の
予測影響要因を入力し該予測影響要因に対応する関数を
前記関数テーブルから読み出し、該読み出した関数を用
いて当該予測対象時点の需要量を予測する予測実行手段
とを含んでなる需要量予測装置。
[Claim 12] A database in which actual data on demand quantity is stored with a plurality of influencing factors that influence increases and decreases in demand quantity, and a database in which the relationship between changes in the influencing factors and changes in demand quantity is discontinuous. Actual data for a recent certain period corresponding to each predetermined discontinuity influencing factor or combination of discontinuity influencing factors is sequentially extracted from the database, and the discontinuity influencing factor is determined based on the extracted actual data. a function deriving means for determining a correlation function between the continuity influence factor and the demand quantity for each combination of discontinuity influence factors, a function table in which the determined function is stored, and a prediction influence factor at the time of prediction target A demand amount forecasting device comprising: a prediction execution unit that inputs a function that corresponds to the prediction influence factor, reads a function corresponding to the prediction influence factor from the function table, and predicts the demand amount at the prediction target time using the read function.
【請求項13】  需要量の実績データが需要量の増減
に影響を与える複数の影響要因を付して格納されるデー
タベースと、影響要因の変化と需要量の変化との関係が
不連続的な予め定められた不連続性影響要因ごと又は不
連続性影響要因の組合せごとに対応する最近の一定期間
の実績データを前記データベースから順次抽出し、該抽
出した実績データから外れる領域について過去の実績デ
ータに基づいて予測対象時点の疑似データを作成する疑
似データ作成手段と、該疑似データと前記抽出した実績
データとに基づいて不連続性影響要因ごとに又は不連続
性影響要因の組合せごとに、前記連続性影響要因と需要
量の相関関数を求める関数導出手段と、該求められた関
数が格納される関数テーブルと、予測対象時点の予測影
響要因を入力し該予測影響要因に対応する関数を前記関
数テーブルから読み出し、該読み出した関数を用いて当
該予測対象時点の需要量を予測する予測実行手段とを含
んでなる需要量予測装置。
[Claim 13] A database in which the actual data of the demand quantity is stored with a plurality of influencing factors that influence an increase or decrease in the demand quantity, and a database in which the relationship between changes in the influencing factors and changes in the demand quantity is discontinuous. Performance data for a recent certain period corresponding to each predetermined discontinuity influencing factor or combination of discontinuity influencing factors is sequentially extracted from the database, and past performance data is extracted for areas that deviate from the extracted performance data. pseudo data creation means for creating pseudo data at the prediction target point based on the above, and for each discontinuity influencing factor or for each combination of discontinuity influencing factors, based on the pseudo data and the extracted actual data; A function deriving means for obtaining a correlation function between a continuity influencing factor and a demand quantity, a function table in which the determined function is stored, and a function deriving means for inputting a predicted influencing factor at a prediction target time and a function corresponding to the predicted influencing factor. 1. A demand forecasting device comprising a prediction execution unit that reads out a function table and uses the read function to predict a demand at the time of the prediction target.
【請求項14】  需要量の実績データが需要量の増減
に影響を与える複数の影響要因を付して格納されるデー
タベースと、影響要因の変化と需要量の変化との関係が
不連続的な予め定められた不連続性影響要因ごと又は不
連続性影響要因の組合せごとに対応する最近の一定期間
の実績データを前記データベースから順次抽出し、該抽
出した実績データに基づいて不連続性影響要因ごとに又
は不連続性影響要因の組合せごとに、前記連続性影響要
因と需要量の第1の相関関数を求める第1の関数導出手
段と、前記抽出した実績データから外れる領域について
過去の実績データに基づいて予測対象時点の疑似データ
を作成する疑似データ作成手段と、該疑似データと前記
抽出した実績データとに基づいて不連続性影響要因ごと
に又は不連続性影響要因の組合せごとに、前記連続性影
響要因と需要量の第2の相関関数を求める第2の関数導
出手段と、該求められた第1と第2の関数が格納される
関数テーブルと、予測対象時点の予測影響要因を入力し
該予測影響要因に対応する前記第1又は第2関数を前記
関数テーブルから読み出し、該読み出した関数を用いて
当該予測対象時点の需要量を予測する予測実行手段とを
含んでなる需要量予測装置。
[Claim 14] A database in which the actual data on the demand quantity is stored with a plurality of influencing factors that influence increases and decreases in the demand quantity, and a database in which the relationship between changes in the influencing factors and changes in the demand quantity is discontinuous. Actual data for a recent certain period corresponding to each predetermined discontinuity influencing factor or combination of discontinuity influencing factors is sequentially extracted from the database, and the discontinuity influencing factor is determined based on the extracted actual data. a first function deriving means for calculating a first correlation function between the continuity influencing factor and the demand quantity for each case or for each combination of discontinuity influencing factors; pseudo data creation means for creating pseudo data at the prediction target point based on the above, and for each discontinuity influencing factor or for each combination of discontinuity influencing factors, based on the pseudo data and the extracted actual data; a second function deriving means for calculating a second correlation function between the continuity influencing factor and the demand quantity; a function table storing the determined first and second functions; and a prediction influencing factor at the prediction target time point. A forecasting execution means for reading out the first or second function inputted and corresponding to the forecasting influence factor from the function table, and predicting the demand at the forecasting target time using the read function. Prediction device.
【請求項15】  請求項12乃至14のいずれかにお
いて、前記関数導出手段は、前記予測需要量に対応する
需要量の新たな実績データを前記一定期間の実績データ
に加えて前記関数を更新することを特徴とする需要量予
測装置。
15. In any one of claims 12 to 14, the function deriving means updates the function by adding new actual data of the demand amount corresponding to the predicted demand amount to the actual data for the certain period. A demand forecasting device characterized by:
【請求項16】  請求項15において、前記関数導出
手段は、前記更新に係る新たな実績データの期間が所定
の追加期間に達したとき前記関数を初期化し、該新たな
実績データを含む実績データの最近の前記一定期間の実
績データを用いて、前記関数を新たに求めることを特徴
とする需要量予測装置。
16. In claim 15, the function deriving means initializes the function when the period of the new performance data related to the update reaches a predetermined additional period, and the function derivation means initializes the function when the period of the new performance data related to the update reaches a predetermined additional period, and generates the performance data including the new performance data. A demand forecasting device characterized in that the function is newly calculated using actual performance data for the recent certain period.
【請求項17】  請求項12乃至14のいずれかにお
いて、前記関数導出手段は、不連続性影響要因ごとに又
は不連続性影響要因の組合せごとに前記連続性影響要因
と需要量の相関関数を導出する重回帰モデル又はニュー
ラルネットモデルを含んでなることを特徴とする需要量
予測装置。
17. In any one of claims 12 to 14, the function deriving means calculates a correlation function between the continuity influencing factor and the demand quantity for each discontinuity influencing factor or for each combination of discontinuity influencing factors. A demand forecasting device comprising a multiple regression model or a neural network model.
【請求項18】  請求項12乃至14のいずれかにお
いて、前記関数導出手段は、それぞれ不連続性影響要因
ごとに又は不連続性影響要因の組合せごとに前記連続性
影響要因と需要量の相関関数を導出する重回帰モデル又
はニューラルネットモデルを含んでなり、用いる実績デ
ータ等の信頼度に応じて各モデルの収束条件を変えるこ
とを特徴とする需要量予測装置。
18. In any one of claims 12 to 14, the function deriving means calculates a correlation function between the continuity influencing factor and the demand quantity for each discontinuity influencing factor or for each combination of discontinuity influencing factors. 1. A demand forecasting device comprising a multiple regression model or a neural network model for deriving the following, and changing the convergence conditions of each model depending on the reliability of actual performance data, etc. used.
【請求項19】  請求項14において、前記第2の関
数導出手段は、前記実績データを用いる場合よりも前記
疑似データを用いる場合のモデルの収束条件が緩やかに
設定されてなることを特徴とする需要量予測装置。
19. In claim 14, the second function deriving means is characterized in that the convergence conditions of the model are set more gently when using the pseudo data than when using the actual data. Demand forecasting device.
【請求項20】  請求項1乃至11のいずれかにおい
て、前記需要量が電力、水、商品のいずれか一つである
ことを特徴とする需要予測方法。
20. The demand forecasting method according to claim 1, wherein the demand amount is one of electricity, water, and merchandise.
【請求項21】  請求項1乃至11のいずれかにおい
て、前記需要量が電力であり、前記不連続性影響要因が
天候の種類、曜日の種類、社会的なイベント等のいずれ
か一つを含み、前記連続性影響要因が日最高気温、日最
低気温、体感温度、不快指数、湿度、日射量等のいずれ
か一つを含むものであることを特徴とする需要予測方法
21. In any one of claims 1 to 11, the demand amount is electricity, and the discontinuity influencing factor includes any one of weather type, day of the week type, social event, etc. . A demand forecasting method, wherein the continuity influencing factor includes any one of daily maximum temperature, daily minimum temperature, sensible temperature, discomfort index, humidity, amount of solar radiation, etc.
【請求項22】  請求項12乃至19のいずれかにお
いて、前記需要量が電力、水、商品のいずれか一つであ
ることを特徴とする需要予測装置。
22. The demand forecasting device according to claim 12, wherein the demand amount is one of electricity, water, and merchandise.
【請求項23】  請求項12乃至19のいずれかにお
いて、前記需要量が電力であり、前記不連続性影響要因
が天候の種類、曜日の種類、社会的なイベント等のいず
れか一つを含み、前記連続性影響要因が日最高気温、日
最低気温、体感温度、不快指数、湿度、日射量等のいず
れか一つを含むものであることを特徴とする需要予測装
置。
23. In any one of claims 12 to 19, the demand amount is electricity, and the discontinuity influencing factor includes any one of a type of weather, a type of day of the week, a social event, etc. . A demand forecasting device, wherein the continuity influencing factor includes any one of daily maximum temperature, daily minimum temperature, sensible temperature, discomfort index, humidity, amount of solar radiation, etc.
JP14883291A 1991-06-20 1991-06-20 Demand forecasting method and device Expired - Lifetime JP2681421B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP14883291A JP2681421B2 (en) 1991-06-20 1991-06-20 Demand forecasting method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP14883291A JP2681421B2 (en) 1991-06-20 1991-06-20 Demand forecasting method and device

Publications (2)

Publication Number Publication Date
JPH04372046A true JPH04372046A (en) 1992-12-25
JP2681421B2 JP2681421B2 (en) 1997-11-26

Family

ID=15461732

Family Applications (1)

Application Number Title Priority Date Filing Date
JP14883291A Expired - Lifetime JP2681421B2 (en) 1991-06-20 1991-06-20 Demand forecasting method and device

Country Status (1)

Country Link
JP (1) JP2681421B2 (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004355460A (en) * 2003-05-30 2004-12-16 Sanki Eng Co Ltd Method for forecasting energy demand and method for managing energy
WO2008117392A1 (en) * 2007-03-26 2008-10-02 Vpec, Inc. Power system
JP2008259402A (en) * 2007-03-09 2008-10-23 Sanyo Electric Co Ltd Demand control system, demand controller, demand program, and demand control method
JP2009237832A (en) * 2008-03-26 2009-10-15 Tokyo Gas Co Ltd Variable prediction model construction method and variable prediction model construction system
JP2012044739A (en) * 2010-08-16 2012-03-01 Tokyo Electric Power Co Inc:The Load estimation method of power distribution section and system for controlling power distribution system
JP2012079229A (en) * 2010-10-05 2012-04-19 Hitachi Ltd Demand prediction device and water operation monitoring system
JP2014153180A (en) * 2013-02-08 2014-08-25 Mitsubishi Electric Corp Data processing device and program
WO2014141841A1 (en) * 2013-03-15 2014-09-18 株式会社 東芝 Power demand estimating device, method, and program, and demand restraint planning device
JP2015097059A (en) * 2013-11-15 2015-05-21 株式会社日立パワーソリューションズ Power demand prediction system and power demand prediction method
JP2016194845A (en) * 2015-04-01 2016-11-17 株式会社東芝 Accident occurrence prediction device and accident occurrence prediction method
WO2017222030A1 (en) * 2016-06-22 2017-12-28 真実 梶田 Event prediction device and event prediction method
JP2018010521A (en) * 2016-07-14 2018-01-18 新日鐵住金株式会社 Product state prediction device, product state control device, product state prediction method and program
JP2018124727A (en) * 2017-01-31 2018-08-09 株式会社東芝 Electric power demand prediction device
WO2020189585A1 (en) * 2019-03-15 2020-09-24 住友重機械工業株式会社 Monitoring device, display device, monitoring method and monitoring program
WO2022113201A1 (en) * 2020-11-25 2022-06-02 三菱電機ビルテクノサービス株式会社 Power consumption monitoring device and program
WO2023218593A1 (en) * 2022-05-12 2023-11-16 日本電信電話株式会社 Estimation device, estimation method, and estimation program

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5614366A (en) * 1979-07-13 1981-02-12 Toshiba Corp Forecasting unit
JPS5616261A (en) * 1979-07-17 1981-02-17 Toshiba Corp Forecasting method for water demand
JPS6053904A (en) * 1983-09-05 1985-03-28 Hitachi Ltd Ridge type light guide
JPS6237763A (en) * 1985-08-12 1987-02-18 Fujitsu Ltd Automatic calculation system for catching rate of capacity planning
JPH02189492A (en) * 1989-01-18 1990-07-25 Koito Ind Ltd Estimating method for freezing of road surface
JPH03134703A (en) * 1989-10-20 1991-06-07 Hitachi Ltd Demand pattern predicting method and plant operating method
JPH03135496A (en) * 1989-10-20 1991-06-10 Hitachi Ltd Method for controlling residual chlorine concentration in pipeline
JPH03202515A (en) * 1989-12-28 1991-09-04 Meidensha Corp Device for predicting demand of water delivery amount
JPH04195603A (en) * 1990-11-28 1992-07-15 Toshiba Corp Daily water distribution amount prediction device
JPH04277860A (en) * 1991-03-06 1992-10-02 Hitachi Ltd Event numerical correlation display system, event numerical correlation model preparing system and intension decision supporting system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5614366A (en) * 1979-07-13 1981-02-12 Toshiba Corp Forecasting unit
JPS5616261A (en) * 1979-07-17 1981-02-17 Toshiba Corp Forecasting method for water demand
JPS6053904A (en) * 1983-09-05 1985-03-28 Hitachi Ltd Ridge type light guide
JPS6237763A (en) * 1985-08-12 1987-02-18 Fujitsu Ltd Automatic calculation system for catching rate of capacity planning
JPH02189492A (en) * 1989-01-18 1990-07-25 Koito Ind Ltd Estimating method for freezing of road surface
JPH03134703A (en) * 1989-10-20 1991-06-07 Hitachi Ltd Demand pattern predicting method and plant operating method
JPH03135496A (en) * 1989-10-20 1991-06-10 Hitachi Ltd Method for controlling residual chlorine concentration in pipeline
JPH03202515A (en) * 1989-12-28 1991-09-04 Meidensha Corp Device for predicting demand of water delivery amount
JPH04195603A (en) * 1990-11-28 1992-07-15 Toshiba Corp Daily water distribution amount prediction device
JPH04277860A (en) * 1991-03-06 1992-10-02 Hitachi Ltd Event numerical correlation display system, event numerical correlation model preparing system and intension decision supporting system

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004355460A (en) * 2003-05-30 2004-12-16 Sanki Eng Co Ltd Method for forecasting energy demand and method for managing energy
JP2008259402A (en) * 2007-03-09 2008-10-23 Sanyo Electric Co Ltd Demand control system, demand controller, demand program, and demand control method
WO2008117392A1 (en) * 2007-03-26 2008-10-02 Vpec, Inc. Power system
US8401708B2 (en) 2007-03-26 2013-03-19 Vpec, Inc. Electric power system
JP2009237832A (en) * 2008-03-26 2009-10-15 Tokyo Gas Co Ltd Variable prediction model construction method and variable prediction model construction system
JP2012044739A (en) * 2010-08-16 2012-03-01 Tokyo Electric Power Co Inc:The Load estimation method of power distribution section and system for controlling power distribution system
JP2012079229A (en) * 2010-10-05 2012-04-19 Hitachi Ltd Demand prediction device and water operation monitoring system
JP2014153180A (en) * 2013-02-08 2014-08-25 Mitsubishi Electric Corp Data processing device and program
US10345770B2 (en) 2013-03-15 2019-07-09 Kabushiki Kaisha Toshiba Power demand estimating apparatus, method, program, and demand suppressing schedule planning apparatus
WO2014141841A1 (en) * 2013-03-15 2014-09-18 株式会社 東芝 Power demand estimating device, method, and program, and demand restraint planning device
JP2015097059A (en) * 2013-11-15 2015-05-21 株式会社日立パワーソリューションズ Power demand prediction system and power demand prediction method
JP2016194845A (en) * 2015-04-01 2016-11-17 株式会社東芝 Accident occurrence prediction device and accident occurrence prediction method
JPWO2017222030A1 (en) * 2016-06-22 2019-04-11 真実 梶田 Event prediction apparatus and event prediction method
WO2017222030A1 (en) * 2016-06-22 2017-12-28 真実 梶田 Event prediction device and event prediction method
JP2018010521A (en) * 2016-07-14 2018-01-18 新日鐵住金株式会社 Product state prediction device, product state control device, product state prediction method and program
JP2018124727A (en) * 2017-01-31 2018-08-09 株式会社東芝 Electric power demand prediction device
WO2020189585A1 (en) * 2019-03-15 2020-09-24 住友重機械工業株式会社 Monitoring device, display device, monitoring method and monitoring program
JPWO2020189585A1 (en) * 2019-03-15 2020-09-24
EP3940491A4 (en) * 2019-03-15 2022-05-11 Sumitomo Heavy Industries, Ltd. Monitoring device, display device, monitoring method and monitoring program
TWI791949B (en) * 2019-03-15 2023-02-11 日商住友重機械工業股份有限公司 Monitoring device, display device, monitoring method, and monitoring program
WO2022113201A1 (en) * 2020-11-25 2022-06-02 三菱電機ビルテクノサービス株式会社 Power consumption monitoring device and program
JPWO2022113201A1 (en) * 2020-11-25 2022-06-02
WO2023218593A1 (en) * 2022-05-12 2023-11-16 日本電信電話株式会社 Estimation device, estimation method, and estimation program

Also Published As

Publication number Publication date
JP2681421B2 (en) 1997-11-26

Similar Documents

Publication Publication Date Title
CN111563706A (en) Multivariable logistics freight volume prediction method based on LSTM network
Khotanzad et al. A neuro-fuzzy approach to short-term load forecasting in a price-sensitive environment
CN101383023B (en) Neural network short-term electric load prediction based on sample dynamic organization and temperature compensation
CN112288164B (en) Wind power combined prediction method considering spatial correlation and correcting numerical weather forecast
JPH04372046A (en) Method and device for predicting demand amount
CN111091196B (en) Passenger flow data determination method and device, computer equipment and storage medium
CN110443417A (en) Multi-model integrated load prediction method based on wavelet transformation
CN115270965A (en) Power distribution network line fault prediction method and device
CN103793887B (en) Short-term electric load on-line prediction method based on self-adaptive enhancement algorithm
CN112149879A (en) New energy medium-and-long-term electric quantity prediction method considering macroscopic volatility classification
CN113393057A (en) Wheat yield integrated prediction method based on deep fusion machine learning model
CN113361195B (en) Power distribution network line ultra-short term load prediction method based on LSTM and similar days
Capuno et al. Very Short‐Term Load Forecasting Using Hybrid Algebraic Prediction and Support Vector Regression
CN115186923A (en) Photovoltaic power generation power prediction method and device and electronic equipment
CN116526473A (en) Particle swarm optimization LSTM-based electrothermal load prediction method
CN111697560B (en) Method and system for predicting load of power system based on LSTM
CN115470862A (en) Dynamic self-adaptive load prediction model combination method
CN115496290A (en) Medium-and-long-term runoff time-varying probability prediction method based on &#39;input-structure-parameter&#39; full-factor hierarchical combination optimization
CN112465266A (en) Bus load prediction accuracy analysis method and device and computer equipment
CN114971090A (en) Electric heating load prediction method, system, equipment and medium
CN112036598A (en) Charging pile use information prediction method based on multi-information coupling
CN117930776A (en) Accurate steam supply pre-regulation method and system for white spirit production based on model prediction
JP3650304B2 (en) Demand forecasting device
CN115936791A (en) Short-term and medium-term price prediction method for electric power market
CN112446550B (en) Short-term building load probability density prediction method