JPH0678459A - Demand control unit - Google Patents

Demand control unit

Info

Publication number
JPH0678459A
JPH0678459A JP4102438A JP10243892A JPH0678459A JP H0678459 A JPH0678459 A JP H0678459A JP 4102438 A JP4102438 A JP 4102438A JP 10243892 A JP10243892 A JP 10243892A JP H0678459 A JPH0678459 A JP H0678459A
Authority
JP
Japan
Prior art keywords
demand
time
value
sampling
feature parameter
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
JP4102438A
Other languages
Japanese (ja)
Other versions
JP2671245B2 (en
Inventor
Yuji Sawaguchi
勇治 沢口
Kunio Satodate
邦男 里館
Jun Abe
純 阿部
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.)
Osaki Electric Co Ltd
Original Assignee
Osaki Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Osaki Electric Co Ltd filed Critical Osaki Electric Co Ltd
Priority to JP4102438A priority Critical patent/JP2671245B2/en
Publication of JPH0678459A publication Critical patent/JPH0678459A/en
Application granted granted Critical
Publication of JP2671245B2 publication Critical patent/JP2671245B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

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

Abstract

PURPOSE:To improve a flow-up characteristic and stability in demand control, by carrying out fuzzy inference with regard to a demand characteristic parameter for indicating a historical changing trend of working power or a stationary degree, and calculating a single predictive value from temporary predictive values calculated by a plurality of temporary prediction calculating means. CONSTITUTION:In a fuzzy characteristic operation unit 7, fuzzy inference is carried out with regard to a demand characteristic parameter generated from a demand characteristic parameter extracting unit 2 for indicating a historical changing trend of service power or a stationary degree so that weighting values for each predictive value from predicting units a3, b4, and c5 are determined on the basis of a given fuzzy rule by referring a membership function in a demand characteristic membership function unit 6. Then, a weight adding step is carried out by a weight adding unit 8 so that the determined weight values correspond to each predictive value from the predicting units 3 to 5. A center-of- gravity operation unit 9 generates a final power demand predictive value from a given formula on the basis of the predictive value from the predicting units 3 to 5 together with each corresponding weighting value.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、デマンド時限における
使用電力が契約電力を越えないように制御するデマンド
コントロ−ル装置の改良に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an improvement of a demand control device for controlling power consumption during demand time so as not to exceed contracted power.

【0002】[0002]

【従来の技術】従来、デマンドコントロ−ル装置におい
ては、使用電力の予測演算式として、例えば下記の
(1)式の様な計算式を用い、これによって算出される
値をデマンド時限終了時における使用電力の予測値とし
ていた。
2. Description of the Related Art Conventionally, in a demand control device, a calculation formula such as, for example, the following formula (1) is used as a predictive calculation formula of power consumption, and the value calculated by this formula is used at the end of the demand time period. It was used as the predicted value of power consumption.

【0003】 R=P+(Δp/Δt)×Tn …………(1) ここで、 R :デマンド時限終了時点における使用電力の予測
値 P :デマンド時限における使用電力の現在値 Δp :サンプリング時間における使用電力(使用電力
増加分) Δt :サンプリング時間 Tn :デマンド時限の残り時間 である。実際上はΔpは一定の電力値換算重みをもった
パルスの形で入力されるので、離散的な量である。
R = P + (Δp / Δt) × T n (1) where R: predicted value of power consumption at the end of demand time period P: current value of power consumption at demand time period Δp: sampling time Power consumption (increase in power consumption) Δt: Sampling time T n : Demand time remaining time. Practically, Δp is a discrete quantity because it is input in the form of a pulse having a constant power value conversion weight.

【0004】上記(1)式より明らかなように、デマン
ド時限の残り時間Tn が多い場合においては、サンプリ
ング時間Δtが同一であっても、その間の使用電力Δp
の僅かな変化が予測値Rの大きな変動となってあらわれ
てくることになる。この変動幅は、残り時間Tn が多い
場合、或は、サンプリング時間Δtが短い場合に顕著に
なる。
As is clear from the above formula (1), when the remaining time T n of the demand time period is large, even if the sampling time Δt is the same, the power consumption Δp during that period is
A slight change in the value will appear as a large change in the predicted value R. This fluctuation range becomes significant when the remaining time T n is large or when the sampling time Δt is short.

【0005】デマンド制御の追従性の面、例えば負荷が
急激に重くなった場合には速やかにこれを予測値に反映
させる、という観点からは、サンプリング時間Δtが短
い方が望ましく、また、デマンド制御の安定性、例えば
負荷の短時間の変動や負荷の突入電流による継続的でな
い電力変動によってデマンド遮断が起こるのを避ける、
という観点からは、サンプリング時間Δtが長い方が望
ましい、ということになる。
From the viewpoint of the followability of the demand control, for example, when the load suddenly becomes heavy, this is promptly reflected in the predicted value, it is desirable that the sampling time Δt is short. Stability, e.g. avoiding demand interruption due to short-term fluctuations in the load or inconsistent power fluctuations due to load inrush current,
From this point of view, it is desirable that the sampling time Δt is long.

【0006】従来の予測方式の一つとして、例えばデマ
ンド時限のうち、残り時間Tn が多いとみなせる範囲に
おいては、残り時間Tn が少ないとみなせる範囲に比べ
てサンプリング時間Δtを長くとり、上記の条件を満た
そうと試みるものもあった。すなわち、サンプリング時
間Δtを長くすることによって、該サンプリング時間Δ
tより短い継続時間を持った電力変動を平均化し、予測
値の変動を抑制しようとしたものである。
As one of the conventional prediction methods, for example, in the demand time period, in the range where the remaining time T n can be considered to be large, the sampling time Δt is set to be longer than the range where the remaining time T n is considered to be small, and Some tried to satisfy the conditions of. That is, by increasing the sampling time Δt, the sampling time Δt
This is an attempt to suppress fluctuations in the predicted value by averaging power fluctuations having a duration shorter than t.

【0007】[0007]

【発明が解決しようとする課題】上記の予測値の変動を
抑制する、という目的からすれば、サンプリング時間Δ
tは長ければ長い程効果的であるが、平均化の区間が長
ければ、例えば大容量の負荷が立ち上がった場合におい
ては予測値に実際に影響があらわれるのに時間遅れを伴
うことになり、負荷の大きさによってはデマンド制御
(負荷遮断・復帰)のタイミングを失する可能性もあ
る。
For the purpose of suppressing the fluctuation of the above-mentioned predicted value, the sampling time Δ
The longer t is, the more effective it is. However, if the averaging section is long, for example, when a large-capacity load rises, the predicted value will actually be affected, but there will be a time delay. Depending on the size of the, the timing of demand control (load cutoff / return) may be lost.

【0008】例えば、図11は比較的長くとったサンプ
リング時間Δt内においてパルス(電力パルス)がほぼ
等間隔で入力される場合であり、図12は同じく比較的
長くとったサンプリング時間Δt内において著しくパル
スの到着間隔にかたよりが存在する場合である。前者は
負荷が一定の定常運転中であり、後者は軽負荷から重負
荷への立ち上がりであることが想定される。
For example, FIG. 11 shows a case where pulses (power pulses) are input at substantially equal intervals within a relatively long sampling time Δt, and FIG. 12 remarkably occurs within a relatively long sampling time Δt. This is the case where there is a bias in the arrival intervals of the pulses. It is assumed that the former is in a steady operation with a constant load, and the latter is a rise from a light load to a heavy load.

【0009】以上から従来の予測演算式(1)におい
て、残り時間Tn が多い場合であっても、サンプリング
時間Δtを単純に長くしたのでは、この両者(図11と
図12の状態)の区別が困難であり、前述した様にデマ
ンド制御の遅れにつながる恐れがあった。
From the above, in the conventional prediction calculation formula (1), even if the remaining time T n is large, if the sampling time Δt is simply lengthened, both of them (states of FIGS. 11 and 12) The distinction was difficult, and there was a risk of delaying demand control as described above.

【0010】本発明の目的は、上記の点に鑑み、デマン
ド制御の追従性及び安定性を向上させることのできるデ
マンドコントロ−ル装置を提供することである。
In view of the above points, an object of the present invention is to provide a demand control device capable of improving the followability and stability of demand control.

【0011】[0011]

【課題を解決するための手段】本発明は、デマンド時限
終了時点における使用電力を仮に予測する複数の仮予測
演算手段と、過去のサンプリング時間におけるサンプリ
ング値から、現在までの使用電力の変化傾向や定常であ
る度合いの需要特徴パラメ−タを算出する需要特徴パラ
メ−タ算出手段と、該手段にて算出された需要特徴パラ
メ−タにつき、予め特定の需要パタ−ンに対応させたメ
ンバシップ関数を参照し、所定のファジィル−ルに基づ
いて前記複数の仮予測演算手段からの仮予測値それぞれ
への重み付け量を決定し、これにより前記各仮予測値の
重み付けを行って最終的な単一の予測値を算出するファ
ジィ演算手段とを備えている。
SUMMARY OF THE INVENTION The present invention provides a plurality of provisional prediction calculation means for tentatively predicting power consumption at the end of a demand time period, and a change tendency of power consumption up to the present from sampling values at past sampling times. A demand feature parameter calculating means for calculating a demand feature parameter having a steady degree, and a membership function for the demand feature parameter calculated by the means, which corresponds to a specific demand pattern in advance. , A weighting amount for each of the temporary prediction values from the plurality of temporary prediction calculation means is determined based on a predetermined fuzzy rule, and thereby the temporary prediction values are weighted to obtain a final single value. And a fuzzy calculation means for calculating a predicted value of

【0012】[0012]

【作用】ファジィ演算手段は、現在までの使用電力の変
化傾向や定常である度合いを示す需要特徴パラメ−タに
つき、予め特定の需要パタ−ンに対応させたメンバシッ
プ関数を参照して、所定のファジィル−ルに基づいて複
数の仮予測演算手段からの仮予測値それぞれへの重み付
け量を決定し、この各重み付け量によって前記各仮予測
値に重み付けを行い、つまりファジィ推論を行い、最終
的な単一の予測値を算出している。
The fuzzy calculating means refers to a membership function corresponding to a specific demand pattern in advance for a demand feature parameter indicating a change tendency or a steady state of the power consumption up to the present, and determines the predetermined value. Based on the fuzzy rule, a weighting amount for each of the temporary prediction values from the plurality of temporary prediction calculation means is determined, and each temporary prediction value is weighted by each weighting amount, that is, fuzzy inference is performed, and finally A single predictive value is calculated.

【0013】[0013]

【実施例】以下、本発明を図示の実施例に基づいて詳細
に説明する。
DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be described in detail below based on the illustrated embodiments.

【0014】図1はデマンドコントロ−ル装置11のブ
ロック図を示し、デマンドコントロ−ル装置11は、図
示せぬ負荷の使用電力を計量し、計量値に比例したパル
スを発信する発信器付電力量計1からパルスを受信す
る。このパルスは、独立に予測演算をなす予測器3,
4,5(以下、予測器a,b,cと記す)、及び需要特
徴パラメ−タ抽出器2にそれぞれ入力される。
FIG. 1 is a block diagram of the demand control device 11. The demand control device 11 measures electric power used by a load (not shown) and outputs electric power with a transmitter for transmitting a pulse proportional to the measured value. Receive pulses from meter 1. This pulse is used by the predictor 3, which independently performs prediction calculation.
4, 5 (hereinafter referred to as predictors a, b, c) and a demand feature parameter extractor 2 respectively.

【0015】前記予測器a〜cは、例えば前述の予測演
算式(1)におけるサンプリング時間Δtが互いに異な
るか、或は移動平均電力値を用いる等の特徴によって互
いに他と区別される。
The predictors a to c are distinguished from each other by the characteristics such as different sampling times Δt in the above-mentioned prediction calculation formula (1) or using a moving average power value.

【0016】需要特徴パラメ−タ抽出器2は、端的に言
えば、使用電力が増加しつつあるのか、或は減少しつつ
あるのかを定量的に判別する為のものである。具体的に
は、図2で示すΔt1 〜Δt3 について、例えば次の様
な計算を行い、その判別を行う。
In short, the demand feature parameter extractor 2 is for quantitatively determining whether the power consumption is increasing or decreasing. Specifically, for Δt 1 to Δt 3 shown in FIG. 2, for example, the following calculation is performed and the determination is performed.

【0017】t0 をデマンド時限開始から現時点までの
経過時間、この経過時間t0 の間に使用された電力をP
0 とすると、t0 までの平均使用電力PA0は PA0=(P0 /t0 ) で表される。さらに、図2中の区間A,区間B,区間C
それぞれにおける使用電力をP1 ,P2 ,P3 とする
と、これら各区間毎における平均使用電力PA1,PA2
A3は PA1=(P1 /Δt1 ) PA2=(P2 /Δt2 ) PA3=(P3 /Δt3 ) と表せる。
The time t 0 is the elapsed time from the start of the demand time period to the present time, and the electric power used during this elapsed time t 0 is P
When it is set to 0 , the average power consumption P A0 up to t 0 is represented by P A0 = (P 0 / t 0 ). Furthermore, section A, section B, and section C in FIG.
Letting P 1 , P 2 , and P 3 be the respective power consumptions, the average power consumptions P A1 , P A2 ,
P A3 can be expressed as P A1 = (P 1 / Δt 1 ) P A2 = (P 2 / Δt 2 ) P A3 = (P 3 / Δt 3 ).

【0018】このとき、デマンド時刻t1 についての使
用電力の相対変動量δP1 を次式で定義する。
At this time, the relative variation amount δP 1 of the electric power used at the demand time t 1 is defined by the following equation.

【0019】δP1 =(PA1−PA0)/PA0 同様に、デマンド時刻t2 ,t3 についての使用電力の
相対変動量δP2 ,δP3 を次式で定義する。
ΔP 1 = (P A1 −P A0 ) / P A0 Similarly, the relative fluctuation amounts δP 2 and δP 3 of the used power at the demand times t 2 and t 3 are defined by the following equations.

【0020】 δP2 =(PA2−PA0)/PA0 δP3 =(PA3−PA0)/PA0 さて、この需要特徴パラメ−タ抽出間隔ΔT(サンプリ
ング時間Δtとは無関係に定められている)、すなわち
デマンド時刻t0 からt3 の間の負荷の変動について考
えると、例えば典型的な負荷変動は次のようになる。
ΔP 2 = (P A2 −P A0 ) / P A0 δP 3 = (P A3 −P A0 ) / P A0 Now, the demand feature parameter extraction interval ΔT (determined regardless of the sampling time Δt , I.e., load fluctuations between the demand times t 0 and t 3 , typical load fluctuations are as follows.

【0021】δP1 <δP2 <δP3 ……負荷は漸増し
つつある。
ΔP 1 <δP 2 <δP 3 ... The load is gradually increasing.

【0022】δP1 =δP2 =δP3 ……負荷はほぼ定
常である。
ΔP 1 = δP 2 = δP 3 ... The load is almost steady.

【0023】δP1 >δP2 >δP3 ……負荷は漸減し
つつある。
ΔP 1 > δP 2 > δP 3 ... The load is gradually decreasing.

【0024】図3にこれらの場合の相対変動量δPの代
表線を示す。
FIG. 3 shows a representative line of the relative variation amount δP in these cases.

【0025】例えば、相対変動量δPが常に代表線Aよ
り上側にあるならば、負荷はかなり急激に立ち上がって
おり、代表線Bの付近にあるならば、負荷はほぼ定常で
あり、代表線Cより下側にあるならば、負荷はかなり急
激に立ち下がっていると考えられる。
For example, if the relative variation amount δP is always above the representative line A, the load rises fairly rapidly, and if it is near the representative line B, the load is almost steady, and the representative line C. If it is below, the load is likely to fall fairly quickly.

【0026】一般的な場合を考えると、代表線A〜Cで
示されるような直線状にはならないので、その場合の取
り扱いが必要である。
Considering a general case, the straight lines as shown by the representative lines A to C are not formed, and therefore handling in that case is necessary.

【0027】本実施例では、図1のファジィ特徴演算器
7が、需要特徴メンバシップ関数部6に保持されている
メンバシップ関数を参照することによって、相対変動量
δPの変化が直線状でない場合への対応を行う。
In this embodiment, when the fuzzy feature calculator 7 of FIG. 1 refers to the membership function held in the demand feature membership function unit 6, the relative variation amount δP is not linear. Take action.

【0028】需要特徴メンバシップ関数部6は、例えば
図4(a) 〜(c) で表される様なメンバシップ関数を保持
しており、図4(a) はデマンド時刻t1 における負荷の
増減を判定するためのメンバシップ関数であり、図4
(b),(c) はそれぞれデマンド時刻t2 ,t3 についての
増減判定メンバシップ関数である。
The demand feature membership function section 6 holds a membership function as shown in, for example, FIGS. 4 (a) to 4 (c). FIG. 4 (a) shows the load at the demand time t 1 . This is a membership function for determining increase / decrease, and FIG.
(b) and (c) are increase / decrease determination membership functions for demand times t 2 and t 3 , respectively.

【0029】デマンド時刻t1 における代表線への一致
度は、例えば図5に示されるように、横軸上に相対変動
量δP1 をとった時の、それぞれのメンバシップ関数と
の交点を読み取ることによって計算できる。相対変動量
δP2 ,δP3 についても同様な計算を行うものとす
る。このようにして、例えば代表線Aに対して平均電力
の増減傾向の一致度がCfAt1 ,CfAt2 ,CfAt3
してそれぞれ求められる。この平均値、すなわち (CfAt1 +CfAt2 +CfAt3 )/3 をもって代表線Aに対する平均一致度とする。また、そ
れぞれの一致度CfBt1〜CfBt3 ,CfCt1 〜CfCt3
の、平均一致度からの偏差の自乗平均、すなわち分散
を計算しておく。
The degree of coincidence with the representative line at the demand time t 1 is, for example, as shown in FIG. 5, the intersection with each membership function when the relative variation amount δP 1 is taken on the horizontal axis is read. It can be calculated by Similar calculations are performed for the relative fluctuation amounts δP 2 and δP 3 . Thus, for example degree of coincidence of increasing or decreasing trend of the average power with respect to the representative line A are obtained respectively as Cf At1, Cf At2, Cf At3 . This average value, that is, (Cf At1 + Cf At2 + Cf At3 ) / 3 is taken as the average degree of coincidence with the representative line A. In addition, each of the matching degree Cf Bt1 ~Cf Bt3, Cf Ct1 ~Cf Ct3
The root mean square of the deviation from the average degree of coincidence, that is, the variance is calculated.

【0030】ファジィ特徴演算器7は、この算出した平
均一致度と分散から、図6に示されるようなメンバシッ
プ関数を用いると共に、例えば図7のようなル−ル(ル
−ル1,2……)を用いてファジィ推論を行い、各予測
器a〜cに対する予測値の重み付け量を決定する。
The fuzzy feature calculator 7 uses the membership function as shown in FIG. 6 from the calculated average degree of coincidence and variance, and, for example, a rule as shown in FIG. ..) is used to determine the weighting amount of the prediction value for each of the predictors a to c.

【0031】図7のル−ルが例えば図3の代表線Aに対
応するものだとして、具体的な計算手順について説明す
る。なお、代表線A及びB,Cは直線に限るものではな
く、また、3本以上であってもよいものとする。
A specific calculation procedure will be described assuming that the rule of FIG. 7 corresponds to the representative line A of FIG. 3, for example. The representative lines A, B, and C are not limited to straight lines, and may be three or more.

【0032】この場合、ル−ル中での平均一致度とは、
代表線Aに対する平均一致度のことであり、分散につい
ても同様である。
In this case, the average degree of coincidence in the rule is
It is the average degree of coincidence with the representative line A, and the same applies to the variance.

【0033】図7の、以下に示す「ル−ル1」について
説明する。
"Rule 1" shown in FIG. 7 will be described.

【0034】『if(“平均一致度が高く”and“分
散が小さい”)then(“予測器cの重みを「1」と
する”なおかつ“予測器aの重みを「0」とする”』上
記の“平均一致度が高く”に対しては、図8に示すよう
に、現在の平均一致度と、“高い”に対応するメンバシ
ップ関数の交点の値を計算することによってその所属度
rが得られる。また、“分散が小さい”に対しては、図
9に示すように、現在の分散と、「0」に対応するメン
バシップ関数の交点を計算することによって所属度sが
得られる。
“If (“ the average degree of coincidence is high ”and“ variance is small ”) then (“ the weight of the predictor c is “1” ”and“ the weight of the predictor a is “0” ””] For the above-mentioned “high average degree of coincidence”, as shown in FIG. 8, by calculating the current average degree of coincidence and the value of the intersection of the membership functions corresponding to “high”, the degree of membership r For “small variance”, the membership level s is obtained by calculating the intersection of the current variance and the membership function corresponding to “0”, as shown in FIG. .

【0035】「ル−ル1」のif内では、and演算が
行われる。つまり、例えば所属度rと所属度sのうち、
小さい方の値をとるべく演算(Min演算)が行われる
(小さい方の値をとるのは、後述する予測器の重み決定
時に、負荷の上下変動が大きい時はサンプリング時間の
短い予測器の重みを小さくしたい為である)。この値を
wとする。
An and operation is performed within the "rule 1" if. That is, for example, of the belonging degree r and the belonging degree s,
The calculation (Min calculation) is performed to take the smaller value. (The smaller value is used when determining the weight of the predictor to be described later, and when the vertical fluctuation of the load is large, the weight of the predictor having a short sampling time is used. To make it smaller). Let this value be w.

【0036】「ル−ル1」のthen以下では、“予測
器cの重みを「1」とする”とある。すなわち、この
「ル−ル1」での予測器cに対する重み出力値は、「1
×w=wC1」として計算する。また、“予測器aの重み
を「0」とする”ことから、予測器aについても同様に
して、「0×w=wa1」を計算してこれを予測器aに対
する重み出力値とする。なお、デマンド時限内での時刻
経過とともに変化する関数fC (t) を用いて、予測器c
に対する重み出力値を「fC (t) ×w=wC1」として演
算しても良い。
Below "the rule 1" is "the weight of the predictor c is" 1 "", that is, the weight output value for the predictor c in the "rule 1" is "1
Xw = w C1 ”. Further, since “the weight of the predictor a is“ 0 ””, “0 × w = w a1 ” is similarly calculated for the predictor a, and this is set as the weight output value for the predictor a. . The predictor c is calculated using the function f C (t) that changes with the passage of time within the demand time limit.
The weighted output value with respect to may be calculated as “f C (t) × w = w C1 ”.

【0037】一般に、ル−ルは複数個あるので、各ル−
ル全てにわたって上記演算を行うことにより、例えば予
測器cに対する重みとして、wC1,wC2,wC3,……w
cnを得る。このうち、最も大きい値を予測器cに対する
重みwcとして確定する(Max演算)。同様に、他の
予測器a,bに対してもそれぞれに対応する重みが一つ
確定する。
Generally, since there are a plurality of rules, each rule is
By performing the above-mentioned calculation over all the parameters, for example, as weights for the predictor c, w C1 , w C2 , w C3 , ... W
get cn Of these, the largest value is determined as the weight wc for the predictor c (Max operation). Similarly, one weight corresponding to each of the other predictors a and b is determined.

【0038】図1の重み付加器8は、それぞれの予測器
a〜cから出力される予測値に対して上記演算の結果確
定した重みを対応させるものである。
The weight adder 8 of FIG. 1 associates the predicted values output from the predictors a to c with the weights determined as a result of the above calculation.

【0039】重心演算器9は、予測器a〜cから出力さ
れたそれぞれの予測値(以下、Ra,Rb,Rcと記
す)とそれぞれに対応した各重み値(以下、wa,w
b,wcと記す)とから、例えば次式を用いて単一の予
測値Rを出力する機能を有するものである。
The center-of-gravity computing unit 9 outputs the respective predicted values (hereinafter, referred to as Ra, Rb, Rc) output from the predictors a to c and the respective weight values (hereinafter, wa and w).
b), and a function of outputting a single predicted value R using the following equation, for example.

【0040】 R=(Ra・wa+Rb・wb+Rc・wc)/(wa+wb+wc) 以上の動作は、デマンド時刻計数器10によって発せら
れる適切なタイミングを用いて行われる。
R = (Ra.wa + Rb.wb + Rc.wc) / (wa + wb + wc) The above operation is performed using appropriate timing issued by the demand time counter 10.

【0041】なお、ファジィ推論に用いるメンバシップ
関数はデマンド時刻の経過によって変化してもよいし、
平均使用電力の絶対値によって変化してもよい。又、メ
ンバシップ関数は三角形状または台形状であることを要
件としないし、対称的である必要もない。
The membership function used for fuzzy inference may change with the passage of demand time,
It may change depending on the absolute value of the average power consumption. Also, the membership function need not be triangular or trapezoidal, and need not be symmetric.

【0042】また、需要特徴パラメ−タ抽出間隔ΔTは
定値であることを要しないし、図10で示すように、デ
マンド時限に対して時間経過とともに移動(ΔT1 →Δ
2→ΔT3 )しても差し支えない。又、特徴デ−タサ
ンプリング区間の長さが互いに整数比をなす必要もない
し、区間の数も3個に限定されない。更に又、各区間の
サンプリング開始が同時である必要もない。更に又、予
測器a〜cのサンプリング動作と需要特徴パラメ−タ抽
出器2のサンプリング動作は同期している必要はない。
Further, the demand feature parameter extraction interval ΔT does not need to be a constant value, and as shown in FIG. 10, it moves with the passage of time with respect to the demand period (ΔT 1 → Δ
There is no problem even if T 2 → ΔT 3 ). Further, it is not necessary for the lengths of the characteristic data sampling sections to form an integer ratio with each other, and the number of sections is not limited to three. Furthermore, it is not necessary to start sampling in each section at the same time. Furthermore, the sampling operations of the predictors a to c and the sampling operation of the demand feature parameter extractor 2 do not have to be synchronized.

【0043】更に、需要特徴パラメ−タとして本実施例
では分散を用いているが、これは分散の平方根、すなわ
ち、標準偏差であってもよく、標準偏差を平均一致度で
割った、いわゆる変動係数であっても良いし、分散では
なく偏差の平均であっても良い。
Further, although the variance is used as the demand feature parameter in this embodiment, it may be the square root of the variance, that is, the standard deviation, or the so-called fluctuation obtained by dividing the standard deviation by the average degree of coincidence. It may be a coefficient or an average of deviations instead of variance.

【0044】ここで、ル−ルの意味付けについて説明す
る。
Here, the meaning of rules will be described.

【0045】予測器aは、例えば従来の予測演算式
(1)でサンプリング時間Δtを「7分」としたもので
あり、予測器bは、例えばサンプリング時間を「3分」
としたものであり、予測器cは、例えばサンプリング時
間を「1分」としたものである。したがって、例えば予
測器aは負荷が定常的である場合に適し、予測器cは負
荷が急峻な変化を示すときに適す。
The predictor a has a sampling time Δt of "7 minutes" in the conventional prediction calculation formula (1), and the predictor b has a sampling time of "3 minutes", for example.
The predictor c has, for example, a sampling time of “1 minute”. Therefore, for example, the predictor a is suitable when the load is steady, and the predictor c is suitable when the load shows a sharp change.

【0046】さて、代表線Aに対する平均一致度が高い
ということは、平均電力は概ね上昇傾向にあることを指
示しており、なおかつその分散が小さいと言うことは、
上昇度合いのばらつきが少ない、すなわち一時的な負荷
変動ではなく、確実に負荷が上昇していることを示すも
のと考えられる。
The fact that the average degree of coincidence with the representative line A is high indicates that the average power is on the rise, and the fact that the variance is small means that
This is considered to indicate that there is little variation in the degree of increase, that is, it is not a temporary change in load but that the load is increasing reliably.

【0047】このような状況ではサンプリング時間Δt
の短い予測方式が適することは常識であり、つまりはサ
ンプリング時間Δtの短い、予測器cの重みを「1」に
近いものと設定して差し支えないものである(「ル−ル
1」参照)。
In such a situation, the sampling time Δt
It is common sense that the short prediction method is suitable, that is, the weight of the predictor c having a short sampling time Δt can be set to be close to "1" (see "Rule 1"). .

【0048】逆に、負荷が短時間に変動しているとき、
或は、需要特徴パラメ−タ抽出間隔ΔTの程度の時間内
で負荷が定常的であると判断される場合には、例えばサ
ンプリング時間Δtの長い予測器aの重みを「1」に近
付けることによって、予測値の安定度を向上させられ
る。
On the contrary, when the load fluctuates in a short time,
Alternatively, when it is determined that the load is stationary within the time of about the demand feature parameter extraction interval ΔT, for example, by making the weight of the predictor a having a long sampling time Δt close to “1”. , The stability of the predicted value can be improved.

【0049】このような動作を単純な2値論理的な判断
で行うと、予測器aや予測器cの予測方式毎に予測値は
当然異なるため、負荷の変動が連続的であっても、予測
値は連続にならない不具合が生じるが、本実施例ではフ
ァジィ推論により、予測器の重み付けが連続的に変化す
る構成をとることにより、予測値の変化も連続的である
ように動作させることが可能である。
When such an operation is performed by a simple binary logical judgment, the prediction value naturally differs for each prediction method of the predictor a and the predictor c, so that even if the load fluctuation is continuous, Although the predicted value does not become continuous, the present embodiment adopts a configuration in which the weighting of the predictor is continuously changed by fuzzy inference, so that the predicted value can be continuously changed. It is possible.

【0050】以上の様に、本実施例では、サンプリング
時間Δtの短い予測値が出力されるのは、平均電力が連
続的に上昇している、という客観的かつ確実性の高い条
件下であり、このことはデマンド制御の負荷遮断・復帰
の判断上において従来の予測方式よりも本実施例による
予測方式が有利であることを示している。
As described above, in the present embodiment, the short predicted value of the sampling time Δt is output under the objective and highly reliable condition that the average power continuously increases. This indicates that the prediction method according to the present embodiment is more advantageous than the conventional prediction method in determining the load cutoff / recovery of the demand control.

【0051】本実施例によれば、使用電力を予測する場
合に、従来の予測演算式で例えば単純にデマンド時限の
残り時間によってサンプリング時間を変更していた点を
改め、サンプリング時間を長く設定している場合であっ
ても、電力値の変動が過去の経験値に照らして、“十分
に急峻であるらしい”と推定できる使用状況にある場合
には、使用電力を長時間平均化することなく直ちに予測
値に反映させるようにしている為、デマンド制御の追従
性を向上させることができ、また、使用電力が“ほぼ一
定であるらしい”と推定できる使用状況にある場合に
は、サンプリング時間を長くとった予測演算式を用いる
ようにしている為、パルス到着間隔の揺らぎによる予測
値の変動を抑制でき、デマンド制御の安定性を向上させ
ることができる。
According to the present embodiment, when predicting the power consumption, the point that the sampling time is simply changed by the remaining time of the demand time limit by the conventional prediction calculation formula is amended, and the sampling time is set to be long. Even if it is, the fluctuations in the power value are estimated to be “sufficiently steep” in the light of past experience values. Since it is immediately reflected in the predicted value, the followability of demand control can be improved, and if the power consumption is estimated to be "almost constant," the sampling time can be changed. Since the long prediction calculation formula is used, it is possible to suppress the fluctuation of the prediction value due to the fluctuation of the pulse arrival interval and improve the stability of the demand control.

【0052】(変形例)本実施例では、デマンド時限の
残り時間、発信器付電力量計からのパルス数より得られ
る、サンプリング時間における使用電力であるサンプリ
ング値及びデマンド時限内において該デマンド時限開始
時から現在時までの使用電力の現在値から、デマンド時
限終了時点における使用電力を仮予測する様にしている
が、これに限定されるものではなく、例えば冷凍庫が負
荷であった場合、温度と内部に入れられる容量によって
も仮予測することができる。また、使用電力の現在値
は、デマンド時限開始時から現在時までの値であるとし
ているが、デマンド時限内において停電したような場合
には、必要に応じて前回のデマンド時限において使用さ
れた電力値を適宜用いるようにしても良い。
(Modification) In the present embodiment, the demand time is started within the demand time, the sampling value which is the electric power used in the sampling time obtained from the remaining time of the demand time, the number of pulses from the watt hour meter with a transmitter, and the demand time. From the current value of the power used from the time to the current time, the power used at the end of the demand time period is provisionally predicted, but the present invention is not limited to this.For example, when the freezer is a load, the temperature and It can be tentatively predicted by the capacity put inside. In addition, the current value of power consumption is assumed to be the value from the start of the demand time period to the current time, but if there is a power outage within the demand time period, the power used in the previous demand time period may be changed as necessary. You may make it use a value suitably.

【0053】[0053]

【発明の効果】以上説明したように、本発明によれば、
現在までの使用電力の変化傾向や定常である度合いを示
す需要特徴パラメ−タにつき、予め特定の需要パタ−ン
に対応させたメンバシップ関数を参照して、所定のファ
ジィル−ルに基づいて複数の仮予測演算手段からの仮予
測値それぞれへの重み付け量を決定し、この各重み付け
量によって前記各仮予測値に重み付けを行い、つまりフ
ァジィ推論を行い、最終的な単一の予測値を算出してい
る。
As described above, according to the present invention,
Regarding the demand feature parameters that show the tendency of change in power consumption up to the present and the degree of steady state, a plurality of them are based on a predetermined fuzzy rule by referring to the membership function that corresponds to a specific demand pattern in advance. Determines a weighting amount for each temporary prediction value from the temporary prediction calculation means, weights each temporary prediction value by each weighting amount, that is, performs fuzzy inference, and calculates a final single predicted value. is doing.

【0054】よって、デマンド制御の追従性及び安定性
を向上させることが可能となる。
Therefore, it is possible to improve the followability and stability of demand control.

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

【図1】本発明の一実施例におけるデマンドコントロ−
ル装置を示すブロック図である。
FIG. 1 is a demand controller according to an embodiment of the present invention.
FIG. 3 is a block diagram showing a mobile device.

【図2】本発明の一実施例において需要特徴パラメ−タ
を抽出するためのサンプリング間隔の一例を示す図であ
る。
FIG. 2 is a diagram showing an example of sampling intervals for extracting demand feature parameters in an embodiment of the present invention.

【図3】本発明の一実施例において需要特徴パラメ−タ
抽出するための負荷変動の代表線の一例を示す図であ
る。
FIG. 3 is a diagram showing an example of a representative line of load fluctuation for extracting demand feature parameters in one embodiment of the present invention.

【図4】図3の負荷変動の代表線に対応したメンバシッ
プ関数の一例を示す図である。
FIG. 4 is a diagram showing an example of a membership function corresponding to a representative line of load fluctuations in FIG.

【図5】本発明の一実施例において具体的に平均電力の
相対変動量が得られた場合の代表線に対する一致度の一
例を示す図である。
FIG. 5 is a diagram showing an example of the degree of coincidence with the representative line when the relative fluctuation amount of the average power is specifically obtained in the embodiment of the present invention.

【図6】本発明の一実施例において平均一致度と分散を
メンバシップ関数とした場合及び重み付け関数が定数で
ない場合の一例を示す図である。
FIG. 6 is a diagram showing an example of a case where the average degree of coincidence and variance are membership functions and a case where the weighting function is not a constant in one embodiment of the present invention.

【図7】本発明の一実施例におけるファジィ推論ル−ル
の一例を示す図である。
FIG. 7 is a diagram showing an example of a fuzzy reasoning rule according to an embodiment of the present invention.

【図8】本発明の一実施例において平均一致度のメンバ
シップ関数から具体的な所属度を計算するための一例を
示す図である。
FIG. 8 is a diagram showing an example for calculating a specific degree of belonging from a membership function of the average degree of coincidence in one embodiment of the present invention.

【図9】本発明の一実施例において分散のメンバシップ
関数から具体的な所属度を計算するための一例を示す図
である。
FIG. 9 is a diagram showing an example for calculating a specific degree of membership from a membership function of distribution according to an embodiment of the present invention.

【図10】図2の需要特徴パラメ−タ抽出間隔がデマン
ド時刻の経過とともに移動していく場合の一例を示す図
である。
10 is a diagram showing an example of a case where the demand feature parameter extraction interval of FIG. 2 moves with the passage of demand time.

【図11】負荷変動が殆ど無いと予想される場合におけ
る比較的長いサンプリング時間内での電力パルスの到着
間隔を示す図である。
FIG. 11 is a diagram showing arrival intervals of power pulses within a relatively long sampling time when it is expected that there is almost no load fluctuation.

【図12】負荷変動が大きいと予想される場合における
比較的長いサンプリング時間内での電力パルスの到着間
隔を示す図である。
FIG. 12 is a diagram showing arrival intervals of power pulses within a relatively long sampling time when it is expected that load fluctuations are large.

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

1 発信器付電力量計 2 需要特徴パラメ−タ抽出器 3 予測器(a) 4 予測器(b) 5 予測器(c) 6 需要特徴メンバシップ関数部 7 ファジィ特徴演算器 8 重み付加器 9 重心演算器 10 デマンド時刻係数器 11 デマンドコントロ−ル装置 1 Energy Meter with Transmitter 2 Demand Feature Parameter Extractor 3 Predictor (a) 4 Predictor (b) 5 Predictor (c) 6 Demand Feature Membership Function Part 7 Fuzzy Feature Calculator 8 Weighting Unit 9 Center of gravity calculator 10 Demand time coefficient device 11 Demand control device

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 デマンド時限の残り時間を指示する残時
間計数手段と、デマンド時限終了時点における使用電力
を仮に予測する複数の仮予測演算手段と、過去のサンプ
リング時間におけるサンプリング値から、現在までの使
用電力の変化傾向や定常である度合いの需要特徴パラメ
−タを算出する需要特徴パラメ−タ算出手段と、該需要
特徴パラメ−タ算出手段にて算出された需要特徴パラメ
−タにつき、予め特定の需要パタ−ンに対応させたメン
バシップ関数を参照し、所定のファジィル−ルに基づい
て前記複数の仮予測演算手段からの仮予測値それぞれへ
の重み付け量を決定し、これにより前記各仮予測値の重
み付けを行って最終的な単一の予測値を算出するファジ
ィ演算手段とを備えたデマンドコントロ−ル装置。
1. A remaining time counting means for instructing a remaining time of a demand time period, a plurality of temporary prediction calculation means for temporarily predicting electric power used at the end time of the demand time period, and a sampling value from a past sampling time to a present time. A demand feature parameter calculating means for calculating a change tendency of power consumption and a demand feature parameter of a steady degree, and a demand feature parameter calculated by the demand feature parameter calculating means are specified in advance. The membership function corresponding to the demand pattern is referred to, and the weighting amount for each of the temporary prediction values from the plurality of temporary prediction calculation means is determined based on a predetermined fuzzy rule. A demand control device comprising: a fuzzy calculation means for weighting the prediction values to calculate a final single prediction value.
【請求項2】 複数の仮予測演算手段は、残時間計数手
段より指示される残り時間、発信器付電力量計からのパ
ルス数より得られる、サンプリング時間における使用電
力であるサンプリング値及びデマンド時限内において該
デマンド時限開始時から現在時までの使用電力の現在値
から、デマンド時限終了時点における使用電力を仮に予
測する手段であることを特徴とする請求項1記載のデマ
ンドコントロ−ル装置。
2. The plurality of provisional prediction calculation means are a remaining time instructed by the remaining time counting means, a sampling value which is electric power used at a sampling time obtained from the number of pulses from a watt hour meter with a transmitter, and a demand time limit. 2. The demand control device according to claim 1, wherein the demand control device is means for tentatively predicting the power consumption at the end of the demand time period from the current value of the power consumption from the start of the demand time period to the current time.
JP4102438A 1992-03-30 1992-03-30 Demand control device Expired - Lifetime JP2671245B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP4102438A JP2671245B2 (en) 1992-03-30 1992-03-30 Demand control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP4102438A JP2671245B2 (en) 1992-03-30 1992-03-30 Demand control device

Publications (2)

Publication Number Publication Date
JPH0678459A true JPH0678459A (en) 1994-03-18
JP2671245B2 JP2671245B2 (en) 1997-10-29

Family

ID=14327476

Family Applications (1)

Application Number Title Priority Date Filing Date
JP4102438A Expired - Lifetime JP2671245B2 (en) 1992-03-30 1992-03-30 Demand control device

Country Status (1)

Country Link
JP (1) JP2671245B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7343360B1 (en) * 1998-05-13 2008-03-11 Siemens Power Transmission & Distribution, Inc. Exchange, scheduling and control system for electrical power

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5543934A (en) * 1978-09-21 1980-03-28 Tokyo Electric Power Co Maximum demand power monitor
JPH03128624A (en) * 1989-10-12 1991-05-31 Omron Corp Power supply controller

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5543934A (en) * 1978-09-21 1980-03-28 Tokyo Electric Power Co Maximum demand power monitor
JPH03128624A (en) * 1989-10-12 1991-05-31 Omron Corp Power supply controller

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7343360B1 (en) * 1998-05-13 2008-03-11 Siemens Power Transmission & Distribution, Inc. Exchange, scheduling and control system for electrical power

Also Published As

Publication number Publication date
JP2671245B2 (en) 1997-10-29

Similar Documents

Publication Publication Date Title
CN107578134B (en) Yellow river upstream step reservoir flood control scheduling method considering early warning
WO2012127595A1 (en) Power-leveling control apparatus, power leveling power-storage apparatus, power-leveling control method, and leveling program
JP2008109813A (en) Demand controller and power consumption system
CN104554212A (en) Control method and system of vehicle brakes
KR19980024536A (en) Drainage control device
CN115940224A (en) Charge-discharge control method, device, equipment and medium of energy storage system
CN106875681B (en) A kind of traffic congestion tendency prediction technique and system
CN114275487A (en) Belt speed control method, device and equipment for belt conveyor and storage medium
CN106969465B (en) Central air-conditioning in office building system magnetic suspension cold water host control method and device
Henze et al. The impact of forecasting uncertainty on the performance of a predictive optimal controller for thermal energy storage systems
JP2671245B2 (en) Demand control device
CN102930178A (en) Method for detecting abnormal load data
JPWO2019150586A1 (en) Centralized voltage controller and centralized voltage control system
JP2002165362A (en) Forecasting method of electric energy consumption and control method
CN114021902B (en) Reservoir dispatching method for dynamic planning dimension reduction based on dynamic rope collection and discrete mechanism
CN108292558A (en) Method and system for the cooling system for controlling power equipment
JP7156023B2 (en) Continuous casting operation support device, continuous casting operation support method, and program
CN102156822A (en) Pedestrian traffic data assembly multi-step forecasting method
JPH05265513A (en) Method for controlling, number of pumps in operation at place of drainage pump and device therefor
CN104407256B (en) A kind of Electrical Power System Dynamic state updating method for handling measuring section time delay
CN108984487B (en) Real-time sliding fitting method for water level flow relation
JPH0861244A (en) Rainwater drainage pump station control system
CN116518068A (en) Control method of transmission intermediate shaft brake, vehicle and storage medium
CN108637230A (en) A kind of copper disc casting control method and device
CN115328222A (en) Intelligent control system and method for molten iron packaging in front of blast furnace and computer readable storage medium

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080711

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090711

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100711

Year of fee payment: 13

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110711

Year of fee payment: 14

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120711

Year of fee payment: 15

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120711

Year of fee payment: 15