JP2002204527A - 電力伝送ネットワークの安定性を評価するための方法及び装置 - Google Patents

電力伝送ネットワークの安定性を評価するための方法及び装置

Info

Publication number
JP2002204527A
JP2002204527A JP2001365506A JP2001365506A JP2002204527A JP 2002204527 A JP2002204527 A JP 2002204527A JP 2001365506 A JP2001365506 A JP 2001365506A JP 2001365506 A JP2001365506 A JP 2001365506A JP 2002204527 A JP2002204527 A JP 2002204527A
Authority
JP
Japan
Prior art keywords
current
impedance
stability
pair
measurements
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.)
Pending
Application number
JP2001365506A
Other languages
English (en)
Inventor
Christian Pinzon
ピンツォン クリスティアン
Joachim Bertsch
ベルトシュ ヨアヒム
Christian Rehtanz
レータンツ クリスティアン
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.)
ABB Power Automation AG
Original Assignee
ABB Power Automation AG
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 ABB Power Automation AG filed Critical ABB Power Automation AG
Publication of JP2002204527A publication Critical patent/JP2002204527A/ja
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/242Arrangements for preventing or reducing oscillations of power in networks using phasor measuring units [PMU]
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • 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/22Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units

Abstract

(57)【要約】 【課題】 現在の安定性の限度の値を連続的に計算及び
出力すると共にその品質の測定値を供給する。 【解決手段】 各電圧及び電流フェーザーを有する第1
及び第2測定点(P1,P2)を備えた少なくとも一対
の測定値が処理され、Theveninインピーダンス(Zt)
及び現在の安定性の限度(dS(k))の値が計算さ
れ、妥当性インジケータ(v)が、第1及び第2測定点
の間に差があるかどうか、及び対応する推定された第1
及び第2負荷インピーダンス(Za)の間に差があるか
どうかによって計算される。少なくとも一対の測定値の
すべてに関連するすべての妥当性インジケータ(v)か
ら、品質インジケータ(q)が、Thevenin値(Zt)及
び現在の安定性の限度(dS(k))に関連させて計算
される。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、大規模電力伝送ネ
ットワークに関し、より詳細には、それぞれ請求項1,
10及び11の前文による、電力伝送ネットワークの安
定性を評価するための方法、電力ネットワークの安定性
推定装置及びコンピュータプログラム製品に関する。
【0002】
【従来の技術】電力伝送及び分配システム又はネットワ
ークは、地理的に分割された領域を接続すると共に電線
路間の接続を切換えるための高電圧連絡線を備えてい
る。幾つかの変電所を有するネットワークでの電力生成
及び負荷流れは中央エネルギ管理システムにより制御さ
れている。電力生成及び負荷流れの制御での重要な論点
はネットワークを安定に保つこと、すなわち、電圧崩壊
及び振動を避けることである。電圧の限界に基づいて、
ネットワークの安定性を評価する方法は、電力工業コン
ピュータアプリケーション(PICA)、1997年5
月12〜16日、IEEEの K.Vuらによる論文「電圧の安
定性の限度を推定する局部測定値の使用」、及び199
9年6月の電力システムコンピュータ会議(PSCC)
での「電圧の安定性プレディクター(VIP)及びその
適用」に開示されている。両方の記事はここにインコー
ポレイテッドバイリファレンスされている。これらの記
事は、電圧崩壊への接近を推論するため、局部的に電流
及び電圧を測定する「電圧の不安定性プレディクター」
(VIP)を説明している。VIPの概念は図1に示さ
れている。電力システムの一部分は電源として扱われ、
他の部分は負荷として扱われる。電源は、電源電圧Et
及びThevenin又は電源インピーダンスZtとそのTheven
in等価物1により表わされる。負荷は明らかな負荷イン
ピーダンスZaにより表わされている。Theveninインピ
ーダンスZt及び明らかな負荷インピーダンスZaの両
方はVIP装置2による電流及び電圧測定値から推定さ
れる。これらのインピーダンスの関係は、安定性の限度
又は電力の限度により表わされ、電力システム又はネッ
トワークがどの程度崩壊に近いかを示す。
【0003】第1電圧フェーザーV1及び第1電流フェ
ーザーI1を備えた第1測定点を与えて、第1負荷イン
ピーダンスZa1は として計算される。
【0004】第2電圧フェーザーV2及び第2電流フェ
ーザーI2を備えた第2測定点を与えて、Theveninイン
ピーダンスZtは として計算される。
【0005】インピーダンスと同様にフェーザーはベク
トル量であることに注意しなさい。実際の困難さはThev
eninインピーダンスZtを正確に推定することにある。
特に、第1及び第2測定点による明らかな負荷インピー
ダンスZaが変化しない場合、その後、Theveninインピ
ーダンスZtのための値は無効である。
【0006】米国特許5,631,569には、ネットワークの
パラメータの仮の測定値のためのシステムが示されてい
る。TheveninインピーダンスZtの大きさは電力ネット
ワークの固有波動を利用することにより推定され、測定
電圧及び電流を変化させる。RMS(自乗平均)電圧及
びRMS電流は周期的に測定される。各サイクルのた
め、スカラー見掛け負荷インピーダンスが計算される。
異なる負荷インピーダンスの2つのサイクルが認識され
た場合、Theveninインピーダンスが計算されるのみであ
る。好適な実施例では、負荷インピーダンスの10%の
差は、Theveninインピーダンスを計算するために有用で
あると考えられる。Theveninインピーダンスの連続値が
格納され、統計的データは連続値のために生成されると
共に維持される。標準偏差が1シグマ以下になると、値
の平均が表示される。
【0007】しかし、VIPは、時々の表示のための単
一値だけでなく、安定性の限度の発生及びネットワーク
又は変電所制御システムでの使用に適しているThevenin
インピーダンスZtのための推定値の連続的発生をも要
求する。
【0008】
【発明が解決しようとする課題】そのため、本発明の目
的は、電力伝送ネットワークの安定性を評価するための
方法、ネットワーク又は変電所制御システムでの使用に
適した安定性の限度の連続的発生を提供する電力ネット
ワーク安定性推定装置及びコンピュータプログラム製品
を作り出すことである。
【0009】
【課題を解決するための手段】これらの目的は、それぞ
れ請求項1及び10による、電力伝送ネットワークの安
定性を評価するための方法及び電力ネットワーク安定性
推定装置、及び請求項11による、コンピュータプログ
ラム製品により達成される。
【0010】電力伝送ネットワークの安定性を評価する
ための発明の方法及び装置において、第1及び第2測定
点を含む少なくとも一対の測定値であって、それぞれの
側定点が電圧及び電流フェーザーを含む測定値が処理さ
れ、Theveninインピーダンス及び安定性の限度値が計算
され、妥当性インジケータは、第1及び第2測定点間に
所定の差があるかどうか、そしてこれが事実の場合に
は、対応する第1及び第2負荷インピーダンス間に所定
の差があるかどうかによっても計算される。少なくとも
一対の測定値のすべてと関連するすべての妥当性インジ
ケータから、品質インジケータはThevenin値及び安定性
の限度値に関連させて計算される。
【0011】本発明は、安定性の限度を連続的に計算す
ると共に出力し、一定のその品質を供給することを許容
する。本発明の好適な実施例では、安定性の限度値及び
関連の品質インジケータはネットワーク及び変電所制御
システムに伝送される。例えば、これは制御システムを
ネットワークでの問題に自動的に反応させ、関連の品質
インジケータqsが安定性の限度が高品質である、すな
わち、所定の閾値より高いことを示したときのみ、安定
性の限度により表示される。
【0012】本発明の好適な変形では、妥当性インジケ
ータvは関連の一対の測定値が有効で、不変の又は無効
のTheveninインピーダンスを供給するかどうかを表示す
る。
【0013】本発明の好適な実施例では、装置は品質イ
ンジケータqを表示する手段、及び又は品質インジケー
タq及び安定性の限度値を少なくとも電力ネットワーク
の一部で制御する制御システムに伝送する手段を備えて
いる。
【0014】本発明によるコンピュータプログラム製品
はコンピュータ読取り可能な媒体を備え、そこには、前
記プログラムが負荷された時に作るコンピュータプログ
ラムコード手段を有し、コンピュータは発明の方法によ
る電力伝送ネットワークの安定性を評価するための方法
を実行する。
【0015】さらなる好適な実施例は従属請求項から明
らかである。
【0016】
【発明の実施の形態】本発明の主題は、添付図面に示さ
れた好適な実施例に関連させて、以下の本文において、
より詳細に説明する。
【0017】図面で使用される参照符号及びそれらの意
味は参照符号リストの概要形式で列挙されている。原則
として、図面において同一部分は同一参照符号で与えら
れている。
【0018】図1はネットワークの安定性を評価するた
めの概念構造を概略的に示しており、電力システムの一
部分は電源として取扱われ、別の部分は負荷として扱わ
れる。電源は、Thevenin又は電源インピーダンスとのそ
のThevenin等価物1により表わされている。負荷は見掛
け負荷インピーダンスZaにより表わされる。VIV装
置2は、例えば、変電所のベイレベルのフィーダー又は
変電所電線路に沿った分岐点にあるフェーザー測定装置
でフェーザーデータを決定する。電圧フェーザーは、例
えば、フィーダー又は電線路の電圧を表わし、電流フェ
ーザーはフィーダー又は電線路を流れる電流を表わして
いる。
【0019】フェーザーデータはフェーザーを表わし、
極数であってもよく、その絶対値は品質の実数の大きさ
又はRMS値のいずれか、及びゼロ時間での位相角に対
する位相偏角に対応する。代わりに、フェーザーは実及
び虚部を有する複素数であってもよく、フェーザーは矩
形又は指数表記法を使用してもよい。フェーザーは、電
圧、電流、電力又は相導体又は電子回路と関連する電力
又はエネルギーのような量を表わすために使用されても
よい。対照すると、電力ネットワークで使用される伝統
的な検知装置は、通常、電圧、電流等のRMS値のよう
な、スカラーの平均表示だけを測定する。
【0020】幾つかのVIPの適用において、フェーザ
ーデータは大きい地理学領域上、すなわち、数千キロメ
ートル上に分配されるフェーザー測定装置から収集され
る。これらの異なる電源からのフェーザーデータが結合
して分析される適用のため、それらは共通の位相標準に
関連していなければならない。換言すれば、異なるフェ
ーザー測定装置は所定の精度内に同期されるローカルク
ロックを有していなければならない。フェーザー測定装
置のそのような同期は、好ましくは、公知な時間の分配
システム、例えば、衛星航法システム(GPS)で達成
される。通常の実行では、フェーザーデータ9は、少な
くとも200又は100又は好ましくは40ミリ秒毎
に、好ましくは1ミリ秒以下の時間分割で、決定され
る。本発明の好適な実施例では、時間分割は10ミリ秒
以下であり、これは0.2度の位相誤差に相当する。各
測定値は同期したローカルクロックから得たタイムスタ
ンプと関連する。そのため、フェーザーデータはタイム
スタンプデータを備えている。
【0021】図2は電力ネットワークの所定点での電圧
及び電流のプロットを示し、VIP装置が配置されてい
る。これは連絡線又は負荷のフィーダーにあってもよ
い。負荷がより長くなると、電流は増加し、電圧は降下
する。電流曲線に対する実際の電圧(V/I曲線)は引
出し線3により示されている。実際の性質は全体のネッ
トワークにより影響を受けるが、VIP装置によりロー
カル測定値から推定される。VIP計算はこの曲線の推
定を供給し、それは破線4により示されている。特に、
曲線の勾配はVIP装置により計算されたTheveninイン
ピーダンスZtの大きさ|Zt|に対応する。この推定
曲線を供給して、MVAにおける最大電力は、現在のネ
ットワーク状態の下、前記フィーダーを通って伝送され
てもよい。最大電力は電流と電圧の積が最大となる推定
V/I曲線上の点に対応する。電力の限度は、現在の動
作点5にとって、最大電力と現在伝送されている電力と
の差である。現在の時間の電流及び電圧の大きさをそれ
ぞれI(k)及びV(k)としよう。それから、前記フ
ィーダーに対応する現在の電力の限界又は現在の安定性
の限度dS(k)は、 である。
【0022】要約すると、安定性の限度の計算方法は、
オンライン電流及び電圧フェーザー測定値に基いて、ど
の程度より多くの電力が連絡線を通って送られ、又、ネ
ットワークが崩壊する前にどの程度のより多くの電力が
変電所により引き出されるかの連続的な適応評価を与え
る。これは、重要な状態を示すテルテール破壊事象なし
に、ネットワーク状態がゆっくり不安定状態に移動する
状態において、特に、有利である。
【0023】本発明の概念は、2次元のインピーダンス
平面への投影によって説明される。しかし、本発明は、
同一の概念が電圧、電流及びインピーダンスのためのフ
ェーザーを表わすベクトル値を含んでいる。上記安定性
の限度dSの計算においてのみ、スカラー値、すなわ
ち、対応する値の大きさが使用される。
【0024】図3は、電流軸I及び電圧軸Vにより張ら
れたインピーダンス平面の第1測定点P1及び第2測定
点P2を示す図である。第1測定点P1として同一の見
掛け負荷インピーダンスZaとなるすべての測定点は第
1測定点P1を通る第1インピーダンス及び電流と電圧
軸の交差部にある。所定の最大差以下だけ前記見掛け負
荷インピーダンスZaと異なるすべての測定値は、軸の
交差部を通過し、例えば、見掛けの負荷インピーダンス
Za±10%に対応する第2インピーダンス電線路7及
び第3インピーダンス電線路8により定義される領域に
ある。所定の最大ベクトル差以下だけ第1測定点P1と
異なるすべての測定点は、第1測定点P1に中心がある
円9にある。前記差を計算するため、電圧及び電流値
は、例えば、名目上の動作値に関して、測定される。第
1測定点P1及び第2測定点P2の両方を通過する電源
インピーダンス電線路10は、2つの測定点P1,P2
から計算されたTheveninインピーダンスZtの値に対応
する勾配を有している。第1測定点P1の第1インピー
ダンス電線路6にある測定点は、第1測定点P1に関連
して、電源インピーダンス電線路10の勾配についての
情報を供給しない。
【0025】図4はインピーダンス平面の一組の連続側
定点P(k−n)・・・P(k−2),P(k−1),
P(k)を示す図である。本発明によれば、如何なる対
の点はTheveninインピーダンスZtを決定する試みで使
用されてもよい。本発明は、好ましくは多数の対を試験
し、どの対が有効値を供給し、一定のその品質と共に合
成TheveninインピーダンスZtを計算するかを決定す
る。
【0026】本発明の方法の機能は以下の通りである。
VIP装置2は、電力システムの所定点及び複数の瞬間
で、複数の測定点を繰返し測定し、それぞれは電圧フェ
ーザー及び電流フェーザーからなっている。少なくとも
一対の測定値が決定され、第1測定点P1及び第2測定
点P2を備え、第1測定点P1は第1電圧フェーザーV
1及び第1電流フェーザーI1を備え、第2測定点P2
は第2電圧フェーザーV2及び第2電流フェーザーI2
を備えている。
【0027】本発明の第1の好適な実施例では、各対の
測定値は一対の連続する測定点に対応する。すなわち、
一連の測定点P(k−n)・・・P(k−2),P(k
−1),P(k)を与えると、第1対の側定値は{P
(k−n),P(k−n+1}となり、第2対の測定値
は、{P(k−1),P(k)}まで、{P(k−n+
1),P(k−n+2}等・・・となる。
【0028】本発明の第2の好適な実施例では、測定点
の所定数が測定されると共に格納される。複数対の測定
値はその組のすべての他の部材でその組の測定点の各部
材を対にすることにより決定される。すなわち、一連の
n+1の測定点P(k−n)・・・P(k−2),P
(k−1),P(k)を与えると、対は、m=k−n
で、{P(m),P(m+1)},{P(m),P(m+
2)},{P(m),P(m+3)},・・・{P(m),
P(k)},{P(m+1),P(m+2)},{P(m+
1),P(m+3)},{P(m+1),P(m+
4)},・・・{P(m+1),P(k)},{P(m+
2),P(m+3)},{P(m+2),P(m+
4)},{P(m+2),P(m+5)},・・・{P(m
+2),P(k)},・・・{P(k−3),P(k−
2)},{P(k−3),P(k−1)},{P(k−
3),P(k)},{P(k−2),P(k−1)},{P
(k−2),P(k)},{P(k−1),P(k)}と
なる。古い及び新しい測定点からなる各順序対の測定値
は、以下のステップの繰返し実行における第1測定点P
1及び第2測定点P2に対応する。各対の測定値のた
め、第1及び第2測定点間の変化又は差を示す第1ブー
ル値が計算される。第1ブール値は、第1及び第2測定
点(P1,P2)間に最小差がある場合のみ、真とな
る。
【0029】本発明の好適な実施例では、第1ブール値
は、第1及び第2電圧フェーザー間の電圧差と第1及び
第2電流フェーザー間の電流差の少なくとも1つが所定
の各電圧又は電流差の閾値を超えるかどうかを表示す
る。図3を参照すると、これは、第2測定点P2が軸に
平行な矩形(図示せず)の外側にあり、第1測定点P1
の回りに中心があるかどうかの試験に対応する。本発明
の好適な実施例では、電圧閾値は動作電圧の約1又は
1.5%以下であり、電流閾値は動作電流の約1〜2%
である。
【0030】本発明の好適な別の実施例では、第1ブー
ル値は、第1及び第2測定点間のベクトル差が所定の距
離の閾値を超えたかどうかを表示する。図3を参照する
と、これは、第2測定点P2が第1測定点P1の回りに
中心がある円9又は楕円(図示せず)の外側にあるかど
うかの試験に対応する。
【0031】インピーダンスがベクトル値であるから、
上記説明において言及した円、楕円及び矩形は高次元の
等価物に相当する。
【0032】第1ブール値が偽の場合、その後、第1及
び第2測定点P1,P2間に最小差はなく、すなわち、
第2測定点P2は前記矩形又は楕円内にあり、その後、
測定点P1,P2は本質的に等しいと考えられ、対応す
るカウンタが増加する。各対の測定値が一対の連続測定
点に対応する実施例では、最後の有効なTheveninインピ
ーダンスZtが計算され、そのため、制御システムによ
り表示され、又は制御システムにより使用されているか
ら、この同等はネットワークの状態が本質的に変わらな
かったことを示す。
【0033】第1及び第2測定点P1,P2間に最小差
がある場合、その後、第1測定点から第1負荷インピー
ダンス が計算され、第2測定点から第2負荷インピーダンス が計算され、そして第2ブール値が計算され、第1及び
第2負荷インピーダンス間の差が所定のインピーダンス
の閾値を超えるかどうかを表示する。前記差は、好まし
くは、負荷インピーダンスの1つによって表わされ、例
えば、第2負荷インピーダンスが±5%以上、第1負荷
インピーダンスと異なることで、表わされる。図3を参
照すると、これは第2測定点P2が第2インピーダンス
電線路7と第3インピーダンス電線路8の間の領域の外
側にある状態に対応する。
【0034】第1及び第2負荷インピーダンス間の差が
所定のインピーダンスの閾値を超えない場合、その後、
ネットワークの状態は変化しないが、対の測定値はThev
eninインピーダンスZtの有効な評価を与えるために使
用することはできず、対応するカウンタが増加される。
【0035】第1及び第2負荷インピーダンス間の差が
所定のインピーダンス閾値を超える場合、その後、ネッ
トワークの状態は変化し、有効なTheveninインピーダン
スZtは対の測定値から計算可能であり、対応するカウ
ンタが増加される。TheveninインピーダンスZtは として計算される。今までに説明したステップの結果と
して、各対の測定値は第1及び任意に第2ブーン値に対
応する妥当性インジケータvと関連している。妥当性イ
ンジケータvは関連の対の測定値が有効で、不変の又は
無効のTheveninインピーダンスを与えるかどうかを表示
する。複数の測定点及び関連の対の測定値は、不変で無
効及び有効なTheveninインピーダンスZtに対応する3
つのカウンタの値に関連している。
【0036】例えば、10秒の間に200ms毎に起こ
る一連の測定値は51測定点となる。本発明の第1の好
適な実施例による処理、次々の側定点、すなわち、お互
いに続く測定点だけを対にすることは一組の50対の測
定点を与える。本発明の第2の好適な実施例により、す
べての測定点をお互いに対にすることは一組の50+4
9+48+47+・・・+3+2+1=1275対を与
える。
【0037】今までに説明した方法は、組の対の測定値
のどのくらいの対が不変で無効及び有効なTheveninイン
ピーダンスZtとなるかを決定する。例えば、50対の
測定値のうち、10対が不変値に相当し、5対が無効と
なり、35対が有効となるかもしれず、その組に含まれ
た情報の十分な品質を表示する。他方、40対の不変、
5対の無効及び5対の安定なものは対応するTheveninイ
ンピーダンスZt及び安定性の限度dS値が低品質であ
ることを表示する。
【0038】対の測定値の組のために決定されたTheven
inインピーダンスZtの合成値を得るため、前記Theven
inインピーダンスZtは、例えば、それらを平均化する
ことにより、統計的に濾過され又は処理される。
【0039】さらなるステップでは、対の測定値の組の
すべての対と関連するすべての妥当性インジケータか
ら、品質インジケータqは、Theveninインピーダンス値
及び現在の安定性の限度dS(k)値と関連させて計算
される。例えば、一組のNtot対での有効なTheveninイ
ンピーダンスZtの比較的高い数Nvalは高品質を示
す。本発明の好適な実施例では、品質インジケータqは
また、TheveninインピーダンスZt値及び又は有効なTh
eveninインピーダンスZtの絶対値の一定の分配を組込
む。因子w1,w2,w3、前記分配の標準偏差を重み
付けると、好適な品質インジケータqは と定義される。
【0040】合成TheveninインピーダンスZt及び現在
の測定点P(k)から、ネットワークの安定性を特徴づ
ける現在の安定性の限度dS(k)が計算される。本発
明の好適な実施例では、現時点の安定性の限度dS
(k)は として計算され、V(k)は現在の電圧の大きさで、I
(k)は現在の電流の大きさで、共に現在の測定点P
(k)に相当する。
【0041】安定性の限度dSはそれぞれの有効なThev
eninインピーダンスZtのために形成され、その後、他
の方法で平均化又は濾過されてもよく、又は、安定性の
限度dSは、1つだけの平均化又は濾過されたThevenin
インピーダンスZtのため、本発明の本質を変えること
なく計算されてもよい。
【0042】本発明の好適な実施例では、合成のTheven
inインピーダンスZt及び又は関連の安定性の限度dS
は、数値形式又は図表のいずれかで、例えば、時間上の
軌跡又はトレンドグラフとしてオペレータに表示され
る。品質インジケータqは、例えば、それに応じて表示
数又は時間の履歴を色付けすることにより、図表で表わ
される。別の実施例では、TheveninインピーダンスZt
のトレンドグラフは、有効、無効又は不変のTheveninイ
ンピーダンスZtを供給する対の測定値の相対数の少な
くとも1つのトレンドグラフと共に表示される。
【0043】本発明の方法及び装置は連続的に更新され
る安定性の限度dS及び関連の品質インジケータqを有
するTheveninインピーダンスZtを与え、それらは電力
ネットワークの動作で使用可能である。本発明の好適な
実施例では、これらの値はタイムスタンプ又は時間情報
と関連している。これは電力システムに及ぶ値の統合を
可能にする。
【0044】品質インジケータqにより供給されるさら
なる情報は、以下の例に示されているようにオペレータ
又は自動処理のいずれかによる区別された使用を可能に
する。
【0045】一組の測定値のすべての対の半分以上が、
本質的に等しい電流及び電圧を有する場合、その後、電
力システムは動的変化を示さず、新しい値が計算される
必要がある。TheveninインピーダンスZt又は安定性の
限度dSの現存の有効値は有効のままである。
【0046】一組の測定値のすべての対の半分以上が、
本質的に異なる電流及び電圧を有する場合、その後、電
力システムは動的変化を示すが、大きな負荷が変化する
ことなく、TheveninインピーダンスZt又は安定性の限
度dSが有益な特徴のみを有すべきである。
【0047】一組の測定値のすべての対の半分以上が、
本質的に異なる電流及び電圧を有する場合、半分以上は
本質的に異なる負荷インピーダンスを有し、その後、電
力システムは大きな負荷の変化を伴う動的変化を示す。
TheveninインピーダンスZt又は安定性の限度dSは、
負荷制限のような保護処置で決定するために使用可能で
ある。
【0048】本発明のさらなる好適な実施例では、品質
インジケータ及び現在の安定性の限度dS(k)値によ
り、制御動作が自動的に開始される。これをするため、
TheveninインピーダンスZt及び関連の安定性の限度d
Sは電力ネットワークの制御システム又はエネルギ管理
システム、又は電気変電所の制御システムに伝送され
る。負荷制限のような一定の重要な制御動作は、ネット
ワークの安定性が危険な場合に必要となるので、対応す
る制御プログラムは、安定性の限度dSが所定の最小品
質を有する場合のみ、所定の負荷の制限に影響を与え
る。少ない重要な負荷のため、もっと低品質が要求され
てもよい。
【0049】図5は本発明の方法のフロー図である。ス
テップ51では、タイムスタンプされた電圧及び電流値
が測定される。ステップ52では、これらの値はバッフ
ァに格納され、以下のステップが動作する一組の側定点
を供給する。ステップ53では、一対の測定値が選択さ
れ、2つの測定値の少なくとも電圧又は電流値が異なる
場合に作られるかどうかをチェックする。これがそうで
ない場合、すなわち、電圧と電流の両方の値が変化しな
い場合、その後、ステップ54で、Theveninインピーダ
ンスZtが不変の測定値対の数のためのカウンタが増加
される。そうでなければ、ステップ55では、測定値対
に対応する負荷インピーダンスが異なるかどうかがチェ
ックされる。その後、そうでない場合には、ステップ5
6において、無効なTheveninインピーダンスZtの測定
値対の数のためのカウンタが増加される。その後、そう
である場合には、ステップ57において、Theveninイン
ピーダンスZtが計算されると共に格納され、ステップ
58において、有効なTheveninインピーダンスZtの測
定値対の数のためのカウンタが増加される。それぞれの
増加させるステップ54,56,58の後、ステップ5
9はすべての測定値対が処理されたかどうかを試験す
る。これがそうでない場合、その後、ステップ53で処
理は再開される。そうでなければ、ステップ60におい
て、品質インジケータ及び最良のTheveninインピーダン
スの推定は格納されたTheveninインピーダンスZt及び
カウンタの値から計算される。カウンタはゼロにリセッ
トされ、ステップ61において、最良のTheveninインピ
ーダンスの推定及び品質インジケータは、関連の時間情
報と共に表示され、及び又は制御システムに伝送され
る。
【図面の簡単な説明】
【図1】最先端技術によりネットワークの安定性を評価
するための概念構造を概略的に示している。
【図2】電力ネットワークの選択点での電流及び電圧を
示す図である。
【図3】インピーダンス平面の2つの測定点を示す図で
ある。
【図4】インピーダンス平面の一組の測定点を示す図で
ある。
【図5】本発明の方法のフロー図である。
【符号の説明】
1 Thevenin等価物 2 電圧安定性プレディクタVIP 3 実際のV/I曲線 4 推定V/I曲線 5 動作点 6 第1インピーダンス電線路 7 第2インピーダンス電線路 8 第3インピーダンス電線路 9 円 10 電源インピーダンス電線路 P(k) 現在の測定点 P(k−n)・・・P(k−2),P(k−1),P
(k)測定点 P1 第1測定点 P2 第2測定点 V1 第1電圧フェーザー V2 第2電圧フェーザー V(k) 現在の電圧の大きさ I(k) 現在の電流の大きさ I1 第1電流フェーザー I2 第2電流フェーザー v 妥当性インジケータ q 品質インジケータ dS 安定性の限度 dS(k) 現在の安定性の限度 Za 見掛け負荷インピーダンス Zt Theveninインピーダンス、電源インピーダ
ンス 51 測定及び時間タグ電圧及び電流 52 バッファへの測定値の格納 53 差のための測定値対の試験 54 Theveninインピーダンスが不変の測定値対
の増加数 55 差のための負荷インピーダンス対の試験 56 無効なTheveninインピーダンスの測定値対
の増加数 57 Theveninインピーダンスの計算 58 有効なTheveninインピーダンスの測定値対
の増加数 59 すべての測定値対が処理されたかどうかの
試験 60 品質インジケータ及び最良のTheveninイン
ピーダンスの推定 61 結果の表示及び又は移動
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ヨアヒム ベルトシュ スイス ツェーハー−5405 バーデン−デ ットヴィル ゼゲルホフシュトラッセ 36 エー (72)発明者 クリスティアン レータンツ スイス ツェーハー−5405 バーデン−デ ットヴィル デットヴィラーシュトラッセ 6 Fターム(参考) 5G066 AA03 AE07 AE09

Claims (11)

    【特許請求の範囲】
  1. 【請求項1】 電力伝送ネットワークの安定性を評価す
    るための方法であって、前記ネットワークは、複数の変
    電所、母線及び電線路を備え、 a)電力ネットワークの所定点及び複数の瞬間に、それ
    ぞれ電圧フェーザーと電流フェーザーとからなる複数の
    測定点を測定し、 b)第1測定点(P1)及び第2測定点(P2)を備え
    た少なくとも一対の測定値を決定し、前記第1測定点
    (P1)は第1電圧フェーザーV1と第1電流点I1と
    を備え、前記第2測定点(P2)は第2電圧フェーザー
    V2と第2電流フェーザーI2とを備え、 c)TheveninインピーダンスZtを として計算し、 d)前記TheveninインピーダンスZt及び現在の測定点
    (P(k))からネットワークの安定性を特徴づける現
    在の安定性の限界dS(k)を計算する、ステップを繰
    返すことを含む方法であって、その方法が、少なくとも
    一対の測定値のそれぞれのため、 e)前記第1測定点(P1)と前記第2測定点(P2)
    との間に最小の差が存在する場合にのみ、真となる第1
    ブール値を計算し、 f)前記第1ブール値が真の場合にのみ、前記第1測定
    点(P1)から第1負荷インピーダンス を計算し、 g)前記第1ブール値が真の場合のみ、前記第2測定点
    (P2)から第2負荷インピーダンス を計算し、 h)前記第1ブール値が真の場合のみ、前記第1負荷イ
    ンピーダンスZa1と前記第2負荷インピーダンスZa
    2との差が所定のインピーダンスの閾値を超えたかどう
    かを示す第2ブール値を計算し、 i)前記第1ブール値から、及び前記第1ブール値が真
    の場合には、前記第2ブール値からも、一対の測定値と
    関連する妥当性インジケータ(v)を計算し、そして、
    前記少なくとも一対の測定値のすべてのため、 j)それぞれが別の対の測定値と関連する複数のTheven
    inインピーダンスZtから、合成のTheveninインピーダ
    ンスZtを計算し、 k)それぞれ少なくとも一対の測定値のすべてと関連す
    るすべての妥当性インジケータ(v)から、合成のThev
    eninインピーダンス値及び現在の安定性の限界dS
    (k)と関連する品質インジケータ(q)を計算する、
    ことを含んでいることを特徴とする電力伝送ネットワー
    クの安定性を評価するための方法。
  2. 【請求項2】 現時点での前記安定性の限界dS(k)
    が、 として計算され、V(k)は現在の電圧の大きさであ
    り、I(k)は現在の電流の大きさであり、共に現在の
    測定点(P(k))に対応している請求項1に記載の電
    力伝送ネットワークの安定性を評価するための方法。
  3. 【請求項3】 前記第1及び第2電圧フェーザー間の電
    圧差と、前記第1及び第2電流フェーザー間の電流差の
    少なくとも一方が所定の各閾値を超えた場合にのみ、前
    記第1ブール値が真となる請求項1に記載の電力伝送ネ
    ットワークの安定性を評価するための方法。
  4. 【請求項4】 前記第1及び第2測定点(P1,P2)
    の間のベクトル差が所定の距離の閾値を超えた場合にの
    み、前記第1ブール値が真となる請求項1に記載の電力
    伝送ネットワークの安定性を評価するための方法。
  5. 【請求項5】 前記方法のステップの繰返しの実行が連
    続測定点に適用される請求項1に記載の電力伝送ネット
    ワークの安定性を評価するための方法。
  6. 【請求項6】 ステップa)において、一組の測定点が
    測定され、ステップb)において、複数対の測定値は前
    記組の測定点の各メンバをすべての他のメンバと繰返し
    対にすることにより決定される請求項1に記載の電力伝
    送ネットワークの安定性を評価するための方法。
  7. 【請求項7】 前記関連の対の測定値が、有効、不変又
    は無効のTheveninインピーダンスを与えるかどうかを前
    記妥当性インジケータ(v)が表示する請求項1に記載
    の電力伝送ネットワークの安定性を評価するための方
    法。
  8. 【請求項8】 前記品質インジケータ(q)の視覚表示
    が表示装置に表示される請求項1に記載の電力伝送ネッ
    トワークの安定性を評価するための方法。
  9. 【請求項9】 前記品質インジケータ(q)及び前記現
    在の安定性の限度dS(k)により、制御動作が開始さ
    れる請求項1に記載の電力伝送ネットワークの安定性を
    評価するための方法。
  10. 【請求項10】 電力伝送ネットワークの安定性を推定
    する装置であって、前記ネットワークは、複数の変電
    所、母線及び電線路を備え、前記装置は、電力ネットワ
    ークの所定点及び複数の瞬間に、それぞれ電圧フェーザ
    ーと電流フェーザーとからなる複数の測定点を測定する
    ように構成可能であり、 a)第1測定点P1及び第2測定点P2を備えた少なく
    とも一対の測定値を決定する手段であって、前記第1測
    定点(P1)は第1電圧フェーザーV1と第1電流点I
    1とを備え、前記第2測定点P2は第2電圧フェーザー
    V2と第2電流フェーザーI2とを備えている手段と、 b)TheveninインピーダンスZtを として計算する手段と、 c)前記TheveninインピーダンスZt及び現在の測定点
    (P(k))からネットワークの安定性を特徴づける現
    在の安定性の限界dS(k)を計算する手段と、を備え
    た装置であって、その装置が、 e)少なくとも一対の測定値のそれぞれのため、前記第
    1測定点(P1)と前記第2測定点(P2)との間に最
    小の差が存在する場合にのみ、真となる第1ブール値を
    計算する手段と、 f)前記第1ブール値が真の場合にのみ、前記少なくと
    も一対の測定値のそれぞれのため、前記第1測定点から
    第1負荷インピーダンス を計算する手段と、 g)前記第1ブール値が真の場合のみ、前記少なくとも
    一対の測定値のそれぞれのため、前記第2測定点から第
    2負荷インピーダンス を計算する手段と、 h)前記第1ブール値が真の場合のみ、前記少なくとも
    一対の測定値のそれぞれのため、前記第1負荷インピー
    ダンスZa1と前記第2負荷インピーダンスZa2との
    差が所定のインピーダンスの閾値を超えたかどうかを示
    す第2ブール値を計算する手段と、 i)前記第1ブール値から、及び前記第1ブール値が真
    の場合には、前記第2ブール値からも、一対の測定値と
    関連する妥当性インジケータ(v)を計算する手段とを
    備え、さらに、その装置は、前記少なくとも一対の測定
    値のすべてのため、 j)それぞれが別の対の測定値と関連する複数のTheven
    inインピーダンスZtから、合成のTheveninインピーダ
    ンスZtを計算し、 k)前記少なくとも一対の測定値のすべてとそれぞれ関
    連するすべての妥当性インジケータ(v)から、合成の
    Theveninインピーダンス値及び現在の安定性の限界dS
    (k)と関連する品質インジケータ(q)を計算する、
    手段を備えていることを特徴とする電力伝送ネットワー
    クの安定性を評価する装置。
  11. 【請求項11】 コンピュータで読取り可能な媒体を備
    えたコンピュータプログラム製品であって、前記プログ
    ラムがロードされた時に、作成するコンピュータプログ
    ラムコード手段を有し、前記コンピュータは請求項1に
    よる電力伝送ネットワークの安定性を評価するための方
    法を実行するコンピュータプログラム製品。
JP2001365506A 2000-12-01 2001-11-30 電力伝送ネットワークの安定性を評価するための方法及び装置 Pending JP2002204527A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP00811144A EP1211775B1 (en) 2000-12-01 2000-12-01 Method and device for assessing the stability of an electric power transmission network
EP00811144.5 2000-12-01

Publications (1)

Publication Number Publication Date
JP2002204527A true JP2002204527A (ja) 2002-07-19

Family

ID=8175058

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001365506A Pending JP2002204527A (ja) 2000-12-01 2001-11-30 電力伝送ネットワークの安定性を評価するための方法及び装置

Country Status (5)

Country Link
US (1) US6690175B2 (ja)
EP (1) EP1211775B1 (ja)
JP (1) JP2002204527A (ja)
AT (1) ATE271274T1 (ja)
DE (1) DE60012193D1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015190434A1 (ja) * 2014-06-12 2015-12-17 国立研究開発法人産業技術総合研究所 配電線路のインピーダンス推定装置及び推定方法

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2400580A1 (en) * 2002-09-03 2004-03-03 Sureshchandra B. Patel Systems of advanced super decoupled load-flow computation for electrical power system
EP1416603B1 (en) 2002-11-01 2015-01-21 ABB Research Ltd. Protection of an electric power transmission network
WO2005015366A2 (en) 2003-08-08 2005-02-17 Electric Power Group, Llc Real-time performance monitoring and management system
JP4038484B2 (ja) * 2004-01-13 2008-01-23 三菱電機株式会社 同期フェーザ測定装置
KR100576511B1 (ko) * 2004-05-11 2006-05-03 한국전기연구원 시계열 데이터를 이용한 전력계통 실시간 전압안정도위험지수 계산 시스템 및 방법
US7366623B2 (en) * 2005-05-26 2008-04-29 Texas Instruments Incorporated Method and apparatus for characterizing a load on a data line
CA2585820A1 (fr) * 2007-04-18 2008-10-18 Hydro-Quebec Localisation de defaut sur un reseau electrique par mesures de tension distribuees
GB0712749D0 (en) * 2007-07-02 2007-08-08 Areva T & D Uk Ltd Method of determining voltage stability margin for load shedding within an electrical power system
US7462159B1 (en) 2007-07-24 2008-12-09 Tal Shlomovitz Knee-ankle-foot orthotic device
CN100523840C (zh) * 2007-09-17 2009-08-05 重庆大学 通过辨识电网薄弱环节实时识别电网电压稳定性的方法
MX2010002165A (es) * 2007-09-28 2010-06-02 Schweitzer Engineering Lab Inc Sistemas y metodos para deteccion de variacion de energia y perdida de sincronia utilizando datos marcados en tiempo.
WO2009076769A1 (en) * 2007-12-17 2009-06-25 The University Of British Columbia A method and system for protecting an electrical power transmission network
KR100987167B1 (ko) * 2008-10-01 2010-10-11 한국전력공사 다기 facts 기기의 제어방법 및 장치
US8554385B2 (en) * 2009-09-11 2013-10-08 Schweitzer Engineering Laboratories Inc. Systems and methods for monitoring and controlling electrical system stability
EP2715907B1 (en) * 2011-05-30 2016-10-05 Danmarks Tekniske Universitet Assessment of power systems
WO2013036855A2 (en) * 2011-09-07 2013-03-14 Quanta Associates, L.P. Real-time monitoring of electric power system voltage stability margins
RU2490767C2 (ru) * 2011-11-03 2013-08-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Братский государственный университет" Способ согласования трехпроводной линии электропередачи с электрической нагрузкой
US8773829B2 (en) 2012-10-31 2014-07-08 General Electric Company Method and system for power swing detection in a generator
US9876352B2 (en) * 2013-02-26 2018-01-23 Washington State University Voltage stability monitoring in power systems
RU2547271C2 (ru) * 2013-03-26 2015-04-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Братский государственный университет" Способ согласования несимметричной трехфазной трехпроходной линии электропередачи с электрической нагрузкой
RU2551126C2 (ru) * 2013-07-23 2015-05-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Братский государственный университет" Способ учета стрелы провеса проводов трехфазной трехпроводной линии электропередачи при ее согласовании с электрической нагрузкой
GB2521414B (en) 2013-12-19 2016-01-27 Univ Cape Town Optimal currents for power injection or extraction in a power network
RU2557663C1 (ru) * 2014-01-10 2015-07-27 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Братский государственный университет" Способ учета стрелы провеса линейных проводов трехфазной линии электропередачи при ее согласовании с электрической нагрузкой
RU2558659C1 (ru) * 2014-02-20 2015-08-10 Виктор Александрович Козлов Способ учета стрелы провеса каждого линейного провода трехфазной трехпроводной линии электропередачи при ее согласовании с нагрузкой
RU2561410C1 (ru) * 2014-06-10 2015-08-27 Виктор Александрович Козлов Способ определения места несанкционированного подключения нагрузки неизвестной мощности к линии электрической передачи
US10263426B2 (en) * 2014-10-31 2019-04-16 Hitachi, Ltd. System stabilizing control device and method
RU2577557C1 (ru) * 2015-02-02 2016-03-20 Виктор Александрович Козлов Способ определения места несанкционированного подключения нагрузки к линии электропередачи
DE102016102593B4 (de) 2016-02-15 2017-08-31 Maschinenfabrik Reinhausen Gmbh Verfahren zum Steuern eines Regeltransformators und elektrische Anlage zum Koppeln zweier Wechselstromnetze
US20180269681A1 (en) * 2017-03-14 2018-09-20 V & R Energy Systems Research, Inc. Power systems and related voltage stability methods
US11309702B2 (en) * 2018-02-27 2022-04-19 Hitachi Energy Switzerland Ag Systems and methods for islanding protection

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2184091A1 (en) * 1994-02-28 1995-09-08 Damir Novosel One-terminal data fault location system
US5514969A (en) 1994-05-24 1996-05-07 Reliable Power Meters, Inc. Impedance measurement in a high-voltage power system
US6249719B1 (en) * 1998-05-15 2001-06-19 Abb Power T&D Company, Inc. Applications and methods for voltage instability predictor (VIP)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015190434A1 (ja) * 2014-06-12 2015-12-17 国立研究開発法人産業技術総合研究所 配電線路のインピーダンス推定装置及び推定方法
JP2016001968A (ja) * 2014-06-12 2016-01-07 国立研究開発法人産業技術総合研究所 配電線路のインピーダンス推定装置及び推定方法

Also Published As

Publication number Publication date
ATE271274T1 (de) 2004-07-15
DE60012193D1 (de) 2004-08-19
EP1211775A1 (en) 2002-06-05
US20020097055A1 (en) 2002-07-25
EP1211775B1 (en) 2004-07-14
US6690175B2 (en) 2004-02-10

Similar Documents

Publication Publication Date Title
JP2002204527A (ja) 電力伝送ネットワークの安定性を評価するための方法及び装置
US9178386B2 (en) State-matrix-independent dynamic process estimation method in real-time for weakly observable measurement nodes without PMU
CN100431232C (zh) 电力输送网状态估计
US7494271B2 (en) Power-line sag calculation by way of power-system state estimation
US8812256B2 (en) System and apparatus for measuring the accuracy of a backup time source
Hamzah et al. A new approach to locate the voltage sag source using real current component
US7200500B2 (en) Determining parameters of an equivalent circuit representing a transmission section of an electrical network
EP2770600B1 (en) Method and system for determining power consumption of loads
AU2002302818B2 (en) Method of monitoring a high voltage grid power system
US11784489B2 (en) Supply and demand adjustment monitoring device of power grid and supply and demand adjustment monitoring method for power grid
EP1783499A1 (en) Method and apparatus for verifying the accuracy of state estimation calculations
Geisler et al. State estimation based external network solution for on-line security analysis
US20020123849A1 (en) Method and device for assessing the stability of an electric power transmission network
Tang et al. Online coherency analysis of synchronous generators in a power system
Baran Challenges in state estimation on distribution systems
EP3264114A1 (en) Fault point locating device and method, electric power system monitoring system, and facility planning support system
US10175304B2 (en) Power-system short-circuit capacity monitoring method and system thereof
CN109901022A (zh) 基于同步量测数据的配电网区域定位方法
US7138924B2 (en) Disturbance direction detection in a power monitoring system
CN109061262A (zh) 窃电监控方法及装置
Muñoz-Arango et al. A power-based fault direction estimation method for active distribution networks
RU2785216C1 (ru) Способ анализа качества электрической энергии в трехфазной системе промышленного электроснабжения
JPH11252795A (ja) 電力系統状態決定装置
Babu et al. A-star algorithm based on admissible searching for strategically placing PMU considering redundancy and cost/benefit analysis
JP2001177991A (ja) 電力系統の負荷推定装置