JP2681421B2 - Demand forecasting method and device - Google Patents

Demand forecasting method and device

Info

Publication number
JP2681421B2
JP2681421B2 JP14883291A JP14883291A JP2681421B2 JP 2681421 B2 JP2681421 B2 JP 2681421B2 JP 14883291 A JP14883291 A JP 14883291A JP 14883291 A JP14883291 A JP 14883291A JP 2681421 B2 JP2681421 B2 JP 2681421B2
Authority
JP
Japan
Prior art keywords
function
data
power demand
prediction
actual
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.)
Expired - Lifetime
Application number
JP14883291A
Other languages
Japanese (ja)
Other versions
JPH04372046A (en
Inventor
隆張 石田
岩男 真鳥
重夫 阿部
田村  滋
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

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

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は、電力の将来の需要量を
予測する電力需要量予測装置に関する。
BACKGROUND OF THE INVENTION The present invention is designed to meet future demand for electric power.
The present invention relates to a power demand forecasting device for forecasting .

【0002】[0002]

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

【0003】また、特開昭60−106325公報に記
載されたものによれば、過去の日最高気温と電力需要量
の関係実績から、日最高気温に対する電力需要特性を最
小2乗法を用いた2次式で表現し、この関係式に予測対
象日の予想最高気温を代入することにより、予測対象日
の電力需要を予測するようにしている。
According to the one disclosed in Japanese Patent Laid-Open No. S60-106325, the power demand characteristic with respect to the daily maximum temperature is calculated by using the least squares method based on the past record of the maximum daily temperature and the amount of power demand. It is expressed by the following formula, and by substituting the predicted maximum temperature of the prediction target day into this relational expression, the power demand is predicted.

【0004】[0004]

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

【0005】また、特開昭60−106325公報の技
術によれば、日最高気温と需要量の関係を表す2次式
に、予測対象日の予測日最高気温を単に代入して、予測
対象日の電力需要を予測するもので、過去の気温などの
気象変動による電力需要の変動を考慮していないことか
ら、予測精度に限界があるという問題がある。
Further, according to the technique disclosed in Japanese Patent Laid-Open No. S60-106325, the predicted maximum daily temperature of the predicted target day is simply substituted into the quadratic expression representing the relationship between the maximum daily temperature and the demand amount, and the predicted target day is calculated. There is a problem in that there is a limit to the accuracy of the prediction because it predicts the electric power demand of the above and does not consider the fluctuation of the electric power demand due to the weather change such as the past temperature.

【0006】このような問題は、用水や商品の需要予測
の場合も同様である。
[0006] Such a problem also applies to the demand forecast of water and products.

【0007】本発明の目的は、可能な限り最新の需要量
実績データを用い、かつ予測精度を向上できる電力需要
量予測装置を提供することにある。
An object of the present invention, using the latest demand result data as possible, and power demand that can improve the prediction accuracy
It is to provide a quantity predicting device .

【0008】[0008]

【課題を解決するための手段】上記目的を達成するた
め、本発明の電力需要量予測装置は、電力需要量の実績
データが電力需要量の増減に影響を与える複数の影響要
因を付して前記コンピュータの記憶手段に格納してなる
データベースと、それら複数の影響要因を影響要因の変
化に対して電力需要量が連続的に増減する連続性影響要
因と不連続的に増減する不連続性影響要因とに別けて、
予め定められた不連続性影響要因ごと又は不連続性影響
要因の組合せごとに対応する最近の一定期間の実績デー
タを前記データベースから順次読み出し、該読み出した
実績データに基づいて不連続性影響要因ごとに又は不連
続性影響要因の組合せごとに、前記連続性影響要因と電
力需要量の第1の相関関数を求める第1の関数導出手段
と、前記読み出した実績データから外れる領域について
過去の実績データに基づいて予測対象時点の疑似データ
を作成する疑似データ作成手段と、該疑似データと前記
読み出した実績データとに基づいて不連続性影響要因ご
とに又は不連続性影響要因の組合せごとに、前記連続性
影響要因と電力需要量の第2の相関関数を求める第2の
関数導出手段と、前記第1と第2の関数導出手段により
求められた第1と第2の関数を前記記憶手段に格納して
なる関数テーブルと、前記コンピュータの入力手段から
入力される予測対象時点の予測影響要因に対応する前記
第1又は第2関数を前記関数テーブルから読み出し、該
読み出した関数を用いて当該予測対象時点の電力需要量
を予測する予測実行手段とをコンピュータを用いて形成
したことにある。
In order to achieve the above-mentioned object, the power demand forecasting apparatus of the present invention is a power demand track record.
There are several impacts that the data will have
It is stored in the storage means of the computer
The database and the multiple influencing factors
Continuity impact that power demand continuously increases and decreases
Factors and discontinuity influencing factors that increase and decrease discontinuously,
Each predetermined discontinuity influence factor or discontinuity influence
Actual data for a certain recent period that corresponds to each combination of factors
Data sequentially read from the database and read
Discontinuity based on actual data
For each combination of continuity influencing factors,
First function deriving means for obtaining a first correlation function of the power demand amount
And about the area deviating from the read actual data
Pseudo data at the time of prediction based on past actual data
Means for creating a pseudo data, the pseudo data and
Factors affecting discontinuity based on the read actual data
And the continuity for each combination of discontinuity influencing factors.
The second correlation function of the influential factors and the power demand
The function deriving means and the first and second function deriving means
Storing the obtained first and second functions in the storage means
From the function table and the input means of the computer
Corresponding to the prediction influencing factors of the input prediction target time
Reading the first or second function from the function table,
Electric power demand at the time of the target of prediction using the read function
Prediction execution means for predicting
I did it.

【0009】[0009]

【0010】[0010]

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

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

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

【0014】[0014]

【0015】[0015]

【0016】[0016]

【0017】[0017]

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

【0018】まず、不連続性影響要因ごとに又はその組
合せごとに、連続性影響要因と需要量の相関関数を求
め、該関数に基づいて予測対象時点の需要量を予測する
ことから、曜日の違いや社会的なイベント等が行われる
特異日についても実績データに基づいた関数により予測
できるので、予測精度が向上する。また、不連続性影響
要因をのぞいて、連続性影響因子についてのみの相関関
数を求めて予測することから、実測データの数が少なく
ても関数の信頼性を高くできるので、最近の一定期間の
実績データを用いることができ、過去の古い実績の影響
を必要最小限にできる。また、関数のモデルを重回帰モ
デル又はニューラルネットを適用すれば、更に実績デー
タの対象期間を短縮できる。
First, the correlation function between the continuity influencing factors and the demand amount is obtained for each discontinuity influencing factor or each combination thereof, and the demand amount at the prediction target time is predicted based on the function, so that the day of the week Since the function based on the actual data can be used to predict even a singular day when a difference or social event occurs, the prediction accuracy is improved. Also, since the correlation function for only the continuity influencing factors is obtained and predicted excluding the discontinuity influencing factors, the reliability of the function can be increased even if the number of actually measured data is small, so The actual data can be used, and the influence of past old actual results can be minimized. Further, if a multiple regression model or a neural network is applied to the function model, the target period of the actual data can be further shortened.

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

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

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

【0022】[0022]

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

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

【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 the accuracy of prediction, the pseudo data creating means 11, the function deriving means (II) 12, and the function table (II) 13 are provided. And function updating means (II) 14
Is added.

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

【0026】(不連続性影響要因の例) ・天候の種類(晴、雨、曇、雪など) ・曜日の種類(平日、日祝日、土曜日、飛び石連休の中
日、工場休日など) ・一日の中の時刻 ・社会的なイベント(スポーツイベントなど) (連続性影響要因の例) ・日最高気温 ・日最低気温 ・その他、体感温度、不快指数、湿度、日射量 本発明では、予測の精度を上げるために、まず、不連続
性の影響要因ごと又はその組合せごとに場合分けし、そ
の場合分けの中で連続性の影響要因と電力需要量との相
関を関数化して予測するようにしている。
(Examples of factors affecting discontinuity) -Type of weather (clear, rain, cloudy, snow, etc.)-Type of day of the week (weekdays, holidays, Saturday, medium day of stepping stone holidays, factory holidays, etc.)-One Time of day ・ Social events (sports events, etc.) (Examples of factors affecting continuity) ・ Maximum daily temperature ・ Minimum daily temperature ・ Others: sensible temperature, discomfort index, humidity, solar radiation amount In order to improve the accuracy, first, the cases are classified according to the discontinuity influencing factors or combinations thereof, and the correlation between the continuity influencing factors and the electric power demand is made into a function in the case classification and the prediction is performed. ing.

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

【0028】関数導出手段(I)3は、不連続性影響要
因ごとに又は不連続性影響要因の組合せごとに、前記連
続性影響要因と需要量の相関関数を求めるようになって
おり、必要な実績データをデータベース2から抽出す
る。本実施例では、天候の種類ごとに、かつ曜日の種類
ごとにデータを抽出し、これらの組合せのそれぞれにつ
いて、予測対象日の時刻に対する電力需要量と日最高気
温と日最低気温との関数を求める。なお、抽出対象の実
績データは関数導出時点から見て、最近の一定期間のも
のに制限するようにしている。
The function deriving means (I) 3 is adapted to obtain a correlation function between the continuity influencing factors and the demand amount for each discontinuity influencing factor or each combination of the discontinuity influencing factors, and is required. The actual 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, a function of the power demand amount and the daily maximum temperature and the minimum daily temperature with respect to the time of the prediction target day is calculated. Ask. It should be noted that the actual data to be extracted is limited to the data for a certain fixed period from the time of deriving the function.

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

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

【0031】一例として、簡単のために3層のニューラ
ルネットに対するバックプロパゲーション法について説
明する。入力層ニューロンからの出力値をXk、中間層
ニューロン殻の出力値をXj、出力層ニューロンからの
出力値をXiとする。これらの中間層と出力層の出力値
Xi、Xjは以下の数式1で求める。
As an example, a backpropagation method for a three-layer neural network will be described for simplicity. The output value from the input layer neuron is Xk, the output value of the intermediate layer neuron shell is Xj, and the output value from the output layer neuron is Xi. The output values Xi and Xj of the intermediate layer and the output layer are obtained by the following mathematical expression 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と異なり、ずれていたとする。学習が終了す
ることは以下の不等式が成立することと等価である。
## EQU1 ## Xi = f (ΣWij · Xj−θi) xj = f (ΣWjk · Xk−θj) where Wpq: coupling lines 52a, 52b ,. . . Weighting coefficient for p, q: i, j, k i, j, k: 1, 2, ... θp: threshold value at which neuron starts firing f: sigmoid function f (n) = (1 / (1 + exp ( -N)) Now, when a certain teacher data is input, it is assumed that the output Xi is different from the expected Xi 0. The end of learning is equivalent to the following inequality.

【0033】[0033]

【数2】Σ(Xi−Xi02<ε ただし Σ:総和を表す記号 ε:予め設定した誤差のしきい値 そこで、この数式2を満たすために、Wijを数式2の
左辺の誤差に比例して逆修正する。このようなニューラ
ルネットを用いて関数導出手段(I)3が構成されてい
る。
Σ (Xi−Xi 0 ) 2 <ε where Σ: symbol representing the sum ε: preset error threshold Therefore, in order to satisfy Expression 2, Wij is set to the error on the left side of Expression 2. Inversely correct in proportion. 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, with regard to the embodiment configured as described above, the more detailed structure and operation of the basic components will be described with reference to FIGS. 6 to 8. Also, the detailed configuration of the additional portion will be described together with the operation. Function deriving means (I) 3
The processing procedure of is as shown in FIG. First, in step 61, in order to obtain all the functions for the combinations of discontinuity influencing factors in the conditions to be predicted, the combinations are sequentially set / changed. In step 62, the function corresponding to the combination of the set / changed discontinuity affecting factors is initialized. That is, it makes a meaningless function. Next, step 6
At 3, the actual data corresponding to the set combination is extracted from the data table 2. For example, from the daily actual data in the recent fixed period, the actual data in which the weather is “fine” and the day of the week is “weekday” is read. Based on the read data, in step 64, the correlation function between the daily maximum temperature or the minimum daily temperature and the electric power demand for each time is derived. The function thus derived is stored in the function table (I) 4. Then, in step 65, it is judged whether or not the derivation of the functions related to all the combinations is completed, and if the derivation is not completed, the process returns to step 61 to change the combinations, and until the derivation of all the functions is completed, The process of is repeated.

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

【0036】[0036]

【数3】 P(i,h)=g(X1,X2,...,Xn)## EQU00003 ## P (i, h) = g (X1, X2, ..., Xn)

【0037】[0037]

【数4】|P(i,h)−P’(i,h)|<ε ただし、記号の意味は以下のとおりである。## EQU4 ## | 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): Electric power demand (MW) at time h on day i P '(i, h): Output of the neural network at time h on day i g: Symbol indicating a functional relationship Xj: Influencing factors representing climate, social events, etc. (j = 1, 2, ...) ε: Threshold of output error To increase the reliability of the function derivation by the neural network, as is generally done. , It is desirable to repeat the learning with teacher data. For example, the learning is repeated 20 times for one combination of the daily maximum temperature and the demand amount, and the combination is further changed and repeated 100 times as a whole. This is because the influence of the previous learning result remains, and the influence is removed by repetition. In this way, the demand quantity prediction function relating to the factors that influence continuity can be obtained with high accuracy. The reason why the discontinuity influencing factor is not made into a function is that even if the weather is the same, for example, the tendency of the demand amount may be completely different depending on the type of day of the week, and the accuracy of the function cannot be improved. is there.

【0039】図7に、予測実行手段7の処理手順を示
す。図示のように、ステップ71にて予測指令手段6か
ら予測対象の予測条件を入力する。予測条件として、例
えば何月何日の電力需要量の日変化を予測するとし、そ
の日の予想天候は晴、曜日は平日、その日の予想日最高
気温と日最低気温等を入力する。次にステップ72で予
測対象時刻の初期値を設定する。次のステップ73に
て、後で説明するが、予測対象の条件が、実績データの
存在する領域内か領域外かの判断、つまり内挿か否かを
判断する。内挿の場合は、ステップ74に進んで関数テ
ーブル3から数式3に示した対応する関数のニューラル
ネットを読みだし、その入力層に予測条件を入力し、そ
の出力を当該時刻の電力需要量の予測値とする。続いて
ステップ75において、全ての時刻について予測を完了
したか否かを判断し、終了するまでステップ72に戻っ
て上記の予測処理を繰り返す。全ての予測が完了したと
きは、ステップ76に進んで、予測結果を出力する。
FIG. 7 shows the processing procedure of the prediction executing means 7. As shown, in step 71, the prediction condition of the prediction target is input from the prediction command means 6. As the prediction condition, for example, it is assumed that the day-to-day change in the power demand on what day is predicted, the expected weather of the day is fine, the weekday is weekday, the expected maximum daily temperature of the day and the minimum daily temperature are input. Next, at step 72, the initial value of the prediction target time is set. In the next step 73, as will be described later, it is determined whether the condition to be predicted is inside or outside the region where the actual data exists, that is, whether or not the condition is interpolation. In the case of interpolation, the process proceeds to step 74, the neural network of the corresponding function shown in Formula 3 is read from the function table 3, the prediction condition is input to the input layer, and the output is the power demand amount at that time. Use as a predicted value. Subsequently, in step 75, it is determined whether or not the prediction has been completed for all times, and the process returns to step 72 until the end, and the above-described prediction process is repeated. When all predictions have been completed, the routine proceeds to step 76, where the prediction result is 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 function in the function table (I) 4 by the procedure shown in FIG. 8 to reflect the new performance data. The actual data at the time corresponding to the above prediction result is stored in the database 2 via the input means 1.
Will be input at any time. The function updating means (I) 5 judges in step 81 of FIG. 8 whether or not new record data has been input. When the result is affirmative, in step 82, the new data is added to derive a corresponding function. , Function table (I) 4 is updated. That is, the weighting factor of the neural network is changed. The processing content of this step 82 is shown in FIG.
Since it corresponds to steps 62, 63, 64, 66, and 67 of, the details are omitted. Hereinafter, each time new performance data is input, the function is repeatedly updated, but it is judged whether or not the repetition period exceeds a preset period (learning addition period) (step 83), and the result is affirmative. If so, the process proceeds to step 84 to initialize the function relation, and newly finds the function of the influence factor of the initial learning period and the power demand amount by using the neural net at each time point, and renews the function. Reset it.

【0041】上記の手順により、電力需要予測を行った
結果の一例を次に説明する。 ・連続性影響要因の条件:日最高気温、日最低気温 ・不連続性影響要因の条件:天候−晴、曜日−平日 ・初期学習期間:1カ月 ・学習追加期間:1カ月 図9にある月の予測を行った結果を示す。○はニューラ
ルネットにおいて、教師データとなる予測実行1カ月前
の実績量であり、●は本需要予測方法により求めた予測
量、*は●に対応する実績量である。影響要因である気
温が○印の存在範囲中にある場合の予測量は実績量とほ
とんど一致し、精度のよい予測を行うことができている
(内挿のケース)。しかし、それ以外の範囲にある予測
量は一致の度合が非常に低く、精度が悪くなる(外挿の
ケース)。
An example of the result of power demand forecasting according to the above procedure will be described below.・ Conditions affecting the continuity: maximum daily temperature, minimum daily temperature ・ Conditions affecting the discontinuity: weather-clear, day-weekday ・ Initial learning period: 1 month ・ Additional learning period: 1 month The results of the prediction are shown. In the neural network, ○ is the actual amount one month before the prediction execution, which is the teacher 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 in the range of existence of the circle, the predicted amount is almost the same as the actual amount, and accurate prediction can be performed (interpolation case). However, the prediction amount in the other range has a very low degree of coincidence, resulting in poor accuracy (case of extrapolation).

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

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

【0044】[0044]

【数5】 (Equation 5)

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

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

【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 of actual data is shown in FIG. As shown in the figure, the average value of the daily mean error is 2.18% when only the interpolation case is included, and 2.84% when the extrapolation case is included, including the extrapolation case. It can be observed that the case is less accurate than the case of interpolation alone.

【0048】ここで、図1実施例の付加部分11〜14
の詳細について説明する。この付加部分は予測精度の向
上のために付加されたものであり、特に疑似データ作成
手段11を設けたことに特徴がある。この疑似データ作
成手段11を設けた理由について説明する。精度のよい
需要予測を行うために、外挿のケースのもとで需要予測
を行うことは望ましくないことを前記した。そこで、外
挿のケースを減少するために、もともと外挿のケースに
当たる予想時点の影響要因を、内挿のケースに変化させ
るため、疑似的な教師データを関数導出手段(II)1
2に入力するようにしたのである。疑似データの一例を
図11を用いて説明する。図11は春先のある月の実績
データをもとに、翌月におけるそれぞれの日々の電力需
要量を求める場合の例である。図中、黒で塗りつぶした
楕円はある月の実績データの範囲、ハッチングを施した
楕円は予測対象であるその翌月の実績データの範囲であ
る。また縦軸は電力需要量、横軸は気温を示す。学習デ
ータとなる前記実績データの存在する気温の範囲中には
予想実行月の実績データが存在しないので予測の実行に
際し、これは外挿のケースに当たり、前述したように予
測の精度が悪くなる。一般に、そのような外挿のケース
を解消するために、教師データをクリアせずに蓄積する
方法がある。しかし、この方法には以下にあげる2つの
欠点がある。
Here, the additional portions 11 to 14 of the embodiment shown in FIG.
Will be described in detail. This added portion is added to improve the prediction accuracy, and is characterized in that the pseudo data creating means 11 is provided. The reason for providing the pseudo data creating means 11 will be described. It was mentioned above that it is not desirable to make a demand forecast under the extrapolation case in order to make an accurate demand forecast. Therefore, in order to reduce the extrapolation case, in order to change the influencing factor at the time of prediction, which originally corresponds to the extrapolation case, to the extrapolation case, the pseudo teacher data is used as the function deriving means (II) 1
I input it in 2. An example of the pseudo data will be described with reference to FIG. FIG. 11 is an example of a case where each daily power demand in the next month is calculated based on the actual data of a certain month in early spring. In the figure, the ellipse filled with black is the range of the actual data of a certain month, and the ellipse hatched is the range of the actual data of the next month which is the prediction target. The vertical axis shows the power demand and the horizontal axis shows the temperature. Since the actual data of the predicted execution month does not exist in the temperature range in which the actual data serving as the learning data exists, this is an extrapolation case when executing the prediction, and the accuracy of the prediction deteriorates as described above. Generally, in order to eliminate such an extrapolation case, there is a method of accumulating teacher 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 influence of past data is great, and the prediction accuracy is lowered. The data represented by the range of an open ellipse is added, and the method of reducing the case of extrapolation is introduced. In principle, pseudo data is created automatically using a computer. However, it is also possible to generate pseudo data by the operator manually inputting the know-how possessed by humans, if necessary.

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

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

【0052】[0052]

【数6】 (Equation 6)

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

【0054】また、関数導出手段(II)12で上記の
ような疑似データを用いて非線形関数を導出する場合、
疑似データの信頼性は100%でないから、本実施例で
はデータの信頼性に応じて関数の導出法を変更させるよ
うにしている。具体的には、ニューラルネットを関数の
モデルとして用いた場合は、ニューラルネットの重み係
数を緩い収束条件で求める一方、実績データが存在する
場合は、ニューラルネットの重み係数をきびしい収束条
件で求める方式とする。つまり、ニューラルネットの出
力値と、疑似データ出力値との間の誤差を、実績データ
を用いた場合よりも、高めに設定するような関数導出法
を用いる。すなわち、数式7の関数型に対し、数式8の
条件を満たすように、ニューラルネットの重み係数を更
新する。
When the function deriving means (II) 12 derives a non-linear function using the above pseudo data,
Since the reliability of the pseudo data is not 100%, the method of deriving the function is changed according to the reliability of the data in this embodiment. Specifically, when a neural network is used as a model of a function, the weighting coefficient of the neural network is obtained under a gentle convergence condition, while when actual data exists, the weighting coefficient of the neural network is obtained under a severe convergence condition. And That is, 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 the actual data is used. That is, the weight coefficient of the neural network is updated so that the condition of Expression 8 is satisfied with respect to the function type of Expression 7.

【0055】[0055]

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

【0056】[0056]

【数8】|P”(i,h)−P’(i,h)|<ε’ ただし、記号の意味は以下のとおりである。## EQU8 ## | P "(i, h) -P '(i, h) | <[epsilon]' 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 electric 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 : Influential factors representing climate, social events, etc. (j = 1, 2, ...) ε ': Output error threshold for pseudo data (ε <
ε ′) The pseudo data generated in this way is added to the actual data, and the function deriving means (II) 12 obtains the function II. The structure of the function deriving means (II) 12 is the same as that of the function deriving means (I) 3, and the description thereof will be omitted. Further, the processing procedure of pseudo data creation and function derivation is inserted in steps 66 and 67 of FIG. 6, and the function I and the function II are obtained in parallel.

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

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

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

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

【0062】[0062]

【数9】 (Equation 9)

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

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

【0065】以上、本発明を電力需要量の予測に適用し
た実施例に基づいて説明したが、前述したように、本発
明は電力に限らず、上水需要量やパンなどのような季節
や気温などの影響要因により需要量が増減する対象に適
用できることは明らかである。
Although the present invention has been described based on the embodiment in which it is applied to the prediction of electric power demand, as described above, the present invention is not limited to electric power, and the demand for water supply, the season such as bread, and the like. It is obvious that it can be applied to the target where the demand amount changes due to the influence factors such as temperature.

【0066】[0066]

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

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

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

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

【0070】前記関数を求める手段として、重回帰モデ
ル又はニューラルネットモデルを用い、実績データ等の
信頼度に応じてモデルの収束条件を変えたり、疑似デー
タを用いる場合のモデルの収束条件を緩やかに設定する
ようにしたものによれば、妥当な予測を実現できる。
As a means for obtaining the function, a multiple regression model or a neural network model is used, the convergence condition of the model is changed according to the reliability of the actual data, or the convergence condition of the model when using pseudo data is loosened. According to what is set, a reasonable prediction can be realized.

【0071】また、実績データからベース除去分を除去
した実績値を用いたり、導出した関数から大きくはずれ
る実績データを棄却して再度関数を導出したり、導出し
た関数と実績データとのずれの程度に応じて関数を補正
するようにしたものによれば、更に予測精度を向上させ
ることができる。
Further, the actual value obtained by removing the base removed portion from the actual data is used, the actual data that greatly deviates from the derived function is rejected, and the function is again derived, or the degree of deviation between the derived function and the actual data According to the one in which the function is corrected according to, the prediction accuracy can be further improved.

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

【図1】本発明の一実施例需要予測装置のブロック構成
図である。
FIG. 1 is a block diagram of a demand forecasting apparatus 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 illustrating an example of a daily change in power demand.

【図4】気温と電力需要量の相関の一例を示す線図であ
る。
FIG. 4 is a diagram showing an example of a correlation between air 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 a processing procedure of a function deriving unit.

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

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

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

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

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

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

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

【図16】導出した関数から外れる実績データの程度に
応じて関数を補正する実施例の処理手順の主要部フロー
チャートである。
FIG. 16 is a main part flowchart of a processing procedure of an embodiment in which a function is corrected according to the degree of actual 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 updating means (II)

───────────────────────────────────────────────────── フロントページの続き (72)発明者 田村 滋 茨城県日立市久慈町4026番地 株式会社 日立製作所 日立研究所内 (56)参考文献 特公 昭60−53904(JP,B2) ─────────────────────────────────────────────────── --- Continuation of the front page (72) Inventor Shigeru Tamura 4026 Kuji-machi, Hitachi City, Ibaraki Prefecture Hitachi Research Laboratory, Hitachi Ltd. (56) References Japanese Patent Publication Sho 60-53904 (JP, B2)

Claims (6)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 電力需要量の実績データが電力需要量の
増減に影響を与える複数の影響要因を付して前記コンピ
ュータの記憶手段に格納してなるデータベースと、 それら複数の影響要因を影響要因の変化に対して電力
要量が連続的に増減する連続性影響要因と不連続的に増
減する不連続性影響要因とに別け、予め定められた不
連続性影響要因ごと又は不連続性影響要因の組合せごと
に対応する最近の一定期間の実績データを前記データベ
ースから順次読み出し、該読み出した実績データに基づ
いて不連続性影響要因ごとに又は不連続性影響要因の組
合せごとに、前記連続性影響要因と電力需要量の第1の
相関関数を求める第1の関数導出手段と、 前記読み出した実績データから外れる領域について過去
の実績データに基づいて予測対象時点の疑似データを作
成する疑似データ作成手段と、 該疑似データと前記読み出した実績データとに基づいて
不連続性影響要因ごとに又は不連続性影響要因の組合せ
ごとに、前記連続性影響要因と電力需要量の第2の相関
関数を求める第2の関数導出手段と、前記第1と第2の関数導出手段により 求められた第1と
第2の関数を前記記憶手段に格納してなる関数テーブル
と、前記コンピュータの入力手段から入力される 予測対象時
点の予測影響要因に対応する前記第1又は第2関数を前
記関数テーブルから読み出し、該読み出した関数を用い
て当該予測対象時点の電力需要量を予測する予測実行手
段とをコンピュータを用いて形成してなる電力需要量予
測装置。
1. A performance data of the power demand is denoted by the plurality of influencing factors affecting the increase or decrease of the power demand said competent
A database stored in the storage means of the computer, and the plurality of influencing factors continuously increase or decrease the power demand amount with respect to changes in the influencing factors Separating the discontinuous impact factors, sequentially reads the performance data of the recent predetermined period corresponding to each combination of pre-per-determined discontinuities influence factors or discontinuities influence factors from said database, said read-out performance for each combination of or discontinuities influencing factors for each discontinuity influence factors on the basis of the data, the first function deriving means for obtaining a first correlation function of the continuity effect factor and the power demand, said read the basis of the random data generating means for generating a pseudo data of the prediction target time based on historical data for the area out of the actual data, to the actual data read said the該疑similar data For each combination of or discontinuities influencing factors for each discontinuity contributors, a second function deriving means for obtaining the second correlation function of the continuity effect factor and the power demand, the first and second a function table in which the first and obtained comprising storing the second function in the memory means by the function deriving unit, the first that corresponds to the predicted impact factors of the prediction target time input from the input means of the computer Alternatively, a power demand forecasting apparatus that uses a computer to form a second function that reads out the second function from the function table and uses the read-out function to forecast the power demand at the time of the prediction target.
【請求項2】 請求項において、前記第1と第2の
数導出手段は、予電力需要量に対応する電力需要量の
新たな実績データを前記一定期間の実績データに加えて
前記関数を更新することを特徴とする電力需要量予測装
置。
2. A method according to claim 1, wherein the first and second function <br/> number deriving means, prediction performance data of the predetermined period of new performance data of the power demand corresponding to the power demand power demand prediction apparatus and updates the function in addition to.
【請求項3】 請求項において、前記第1と第2の
数導出手段は、前記更新に係る新たな実績データの期間
が所定の追加期間に達したとき前記関数を初期化し、該
新たな実績データを含む実績データの最近の前記一定期
間の実績データを用いて、前記関数を新たに求めること
を特徴とする電力需要量予測装置。
3. The method according to claim 2 , wherein the first and second function derivation means initialize the function when a period of new record data related to the update reaches a predetermined additional period. The power demand forecasting apparatus is characterized in that the function is newly obtained by using the actual result data including the new actual result data for the latest fixed period.
【請求項4】 請求項1乃至3のいずれかにおいて、前
記第1と第2の関数導出手段は、前記関数の導出に用い
る実績データとして、実際の実績データからベース分を
除去したデータを用いることを特徴とする電力需要量予
測装置。
4. The method according to claim 1, wherein
The first and second function deriving means are used for deriving the function.
As the actual performance data,
Power demand forecast characterized by using the removed data
Measuring device.
【請求項5】 請求項1乃至3のいずれかにおいて、前
記第1と第2の関数導出手段は、前記関数を導出した
後、該導出した関数から一定量外れる実績データを棄却
し、残った実績データを用いて関数を再導出することを
特徴とする電力需要量予測装置。
5. The method according to any one of claims 1 to 3, wherein
The first and second function deriving means derives the function.
Then, reject the actual data that deviates from the derived function by a certain amount.
And re-derive the function using the remaining actual data.
A characteristic power demand forecasting device.
【請求項6】 請求項1乃至3のいずれかにおいて、前
記第1と第2の関数導出手段は、前記関数を導出した
後、該導出した関数と実績データとを比較し、関数から
の実積データの外れの程度に応じて、当該関数を補正す
ることを特徴とする電力需要量予測装置。
6. The method according to claim 1, wherein
The first and second function deriving means derives the function.
Then, compare the derived function with the actual data,
The function is corrected according to the degree of deviation of the actual product data of
A power demand forecasting device, characterized in that
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 JPH04372046A (en) 1992-12-25
JP2681421B2 true 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)

Families Citing this family (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
JP5288782B2 (en) * 2007-03-09 2013-09-11 三洋電機株式会社 Demand control system, demand controller, demand program, and demand control method
EP2136450A4 (en) * 2007-03-26 2013-06-19 Vpec Inc Power system
JP5119022B2 (en) * 2008-03-26 2013-01-16 東京瓦斯株式会社 Variable prediction model construction method and variable prediction model construction system
JP5565193B2 (en) * 2010-08-16 2014-08-06 東京電力株式会社 Load prediction method for distribution section and distribution system control system
JP5561643B2 (en) * 2010-10-05 2014-07-30 株式会社日立製作所 Demand forecasting device and water operation monitoring system
JP6150545B2 (en) * 2013-02-08 2017-06-21 三菱電機株式会社 Data processing apparatus and program
JP2014180187A (en) 2013-03-15 2014-09-25 Toshiba Corp Power demand prediction device, method and program and demand suppression scheduling device
JP6230055B2 (en) * 2013-11-15 2017-11-15 株式会社日立パワーソリューションズ Electric power demand prediction apparatus and electric power demand prediction method
JP6542008B2 (en) * 2015-04-01 2019-07-10 株式会社東芝 Accident occurrence forecasting device and accident occurrence forecasting method
JP6935928B2 (en) * 2016-06-22 2021-09-15 株式会社Singular Perturbations Event prediction device and event prediction method
JP6834209B2 (en) * 2016-07-14 2021-02-24 日本製鉄株式会社 Product status prediction device, product status control device, product status prediction method and program
JP2018124727A (en) * 2017-01-31 2018-08-09 株式会社東芝 Electric power demand prediction device
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
WO2023218593A1 (en) * 2022-05-12 2023-11-16 日本電信電話株式会社 Estimation device, estimation method, and estimation program

Family Cites Families (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
JPS592945B2 (en) * 1979-07-17 1984-01-21 株式会社東芝 Water demand forecasting method
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
JPH03135496A (en) * 1989-10-20 1991-06-10 Hitachi Ltd Method for controlling residual chlorine concentration in pipeline
JPH03134703A (en) * 1989-10-20 1991-06-07 Hitachi Ltd Demand pattern predicting method and plant operating method
JP2712688B2 (en) * 1989-12-28 1998-02-16 株式会社明電舎 Demand forecasting device for water distribution
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

Also Published As

Publication number Publication date
JPH04372046A (en) 1992-12-25

Similar Documents

Publication Publication Date Title
JP2681421B2 (en) Demand forecasting method and device
USRE42440E1 (en) Robust modeling
CN112149879B (en) New energy medium-and-long-term electric quantity prediction method considering macroscopic volatility classification
Ye et al. An ensemble method for short-term wind power prediction considering error correction strategy
JP3141164B2 (en) Power demand forecasting method and device
JP6481468B2 (en) Predictive distribution estimation system, predictive distribution estimation method, and predictive distribution estimation program
Capuno et al. Very short-term load forecasting using hybrid algebraic prediction and support vector regression
CN108204944A (en) The Buried Pipeline rate prediction method of LSSVM based on APSO optimizations
CN112884236B (en) Short-term load prediction method and system based on VDM decomposition and LSTM improvement
Dekkers et al. Optimum linear selection indexes for multiple generation objectives with non-linear profit functions
CN114362175A (en) Wind power prediction method and system based on depth certainty strategy gradient algorithm
CN116526473A (en) Particle swarm optimization LSTM-based electrothermal load prediction method
CN113537582A (en) Photovoltaic power ultra-short-term prediction method based on short-wave radiation correction
CN114065649A (en) Short-term prediction method and prediction system for top layer oil temperature of distribution transformer
CN115049115A (en) RDPG wind speed correction method considering NWP wind speed transverse and longitudinal errors
CN114118634A (en) Soil moisture prediction method and device
CN114330814A (en) Short-term load prediction method based on VMD decomposition and improved double-layer BILSTM network
CN111461378A (en) Power grid load prediction method and device
CN111697560B (en) Method and system for predicting load of power system based on LSTM
CN115496290A (en) Medium-and-long-term runoff time-varying probability prediction method based on &#39;input-structure-parameter&#39; full-factor hierarchical combination optimization
CN111080000A (en) Ultra-short term bus load prediction method based on PSR-DBN
CN116205508A (en) Distributed photovoltaic power generation abnormality diagnosis method and system
CN114881506A (en) Heat supply demand load assessment method and system based on room temperature and IBA-LSTM
CN115293406A (en) Photovoltaic power generation power prediction method based on Catboost and Radam-LSTM
CN116796635A (en) Short-term load prediction method based on CNN-CS-ILSTM model