JPS6035894B2 - 脱調予測方法 - Google Patents

脱調予測方法

Info

Publication number
JPS6035894B2
JPS6035894B2 JP55018037A JP1803780A JPS6035894B2 JP S6035894 B2 JPS6035894 B2 JP S6035894B2 JP 55018037 A JP55018037 A JP 55018037A JP 1803780 A JP1803780 A JP 1803780A JP S6035894 B2 JPS6035894 B2 JP S6035894B2
Authority
JP
Japan
Prior art keywords
power system
generator
potential energy
equation
energy
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
Application number
JP55018037A
Other languages
English (en)
Other versions
JPS56115138A (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.)
Mitsubishi Electric Corp
Tokyo Electric Power Co Holdings Inc
Original Assignee
Tokyo Electric Power Co Inc
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tokyo Electric Power Co Inc, Mitsubishi Electric Corp filed Critical Tokyo Electric Power Co Inc
Priority to JP55018037A priority Critical patent/JPS6035894B2/ja
Publication of JPS56115138A publication Critical patent/JPS56115138A/ja
Publication of JPS6035894B2 publication Critical patent/JPS6035894B2/ja
Expired legal-status Critical Current

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Description

【発明の詳細な説明】 この発明は、電力系統がある規定事故に対して脱調つま
り不安定になるか否かを判定する手法に関するものであ
る。
電力系統のふるまいは、非線形の微分方程式により表現
される。
そのため従来は、これをステップバィステップで解くや
り方と、エネルギーを計算してこれがある値を越えるか
どうかで判定する手法の2通りがあった。まず従釆のシ
ミュレーション手法について説明する。
電力系統は、発電機、送電線、負荷より構成されており
、それらは以下のような非線形の微分方程式で表現され
る。Pi=2ViVi{Giicos(8i−8i)十
Bijsin((8i−8j)} ……(1ー
Qi=2ViVi{GiPin(8i−8j)−Bij
cos(ai−8j)} .・..・.(2
)PLi=k,十k2Vi+K3Vi2
……{3,QLi=k′+k2′Vj十k′3Vi
2 ・・・・・・{4}Pしi+Pi=0
・・・・・・【5}QLi+Q
i=0 ・・・・・・■Mia
i=PMi−Pi(i:発電機に対応するノード番号
)・・・・・・(7)第1図は電力系統の一例を示す構
成図で、1は電力系統をVi,aiはノードiの電圧ベ
クトルの大きさ及びアングルを、Pi,Qiはノードで
より出ていく有効電力、無効電力をPしi,Qしiは、
ノードiの負荷のとる有効電力、無効電力を示す。
‘1},‘2)式はノードiにつながっている送電線よ
り出ていく有効、無効電力を、【3’,‘4’式は負荷
の電圧に対する特性を示す。
{3’,■式のk,,k,′は定電力分をk2,k2′
は定電流分をk3,k3′は定インピーダンス分を示し
ている。
又(5),■式はこれらの和が零となることを示してい
る。(7}式は発電機の運動方程式を表わしMjは慣性
定数をPMiは機械入力を表わす。機械入力PMiは普
通過渡安定度を問題にする時は一定として扱う。又{G
ii+iBii}は電力系統のアドミタンス・マトリク
スを表わす。
このように電力系統のふるまいは、比較的時間の短い範
囲を問題とする過渡安定度の場合m〜{7}式で表現で
きるため、この方程式に初期値を与え、想定される事故
や操作シークェンスに従って各時刻の値を解いていき発
電機のアングルを見て同期がとれている。
つまり安定かどうかの判定をするのがシミュレーション
手法である。
一方エネルギー法というのは、負荷をも含めた系統のア
ドミタンスマトリクスを{GTii+由Tij}とした
時GTii(j≠i)を無視してを求める。
ここでois、(i=1〜n)は安定平衡点のアングル
を示す。なおnは発電機数を示す。そしてWc=W(o
,aic)として求めたWcと事故クリア直後のWの値
を比較してWSWcなら安定、W>Wcなら不安定と判
定する。
従来のシミュレーション手法は、事故後(又は事故後の
一連の操作のあと)の安定を判別するのにかなり長い間
(1)〜‘7}式の微分方程式を解かねばならず、又発
電機のアングルから安定・不安定を判定するのは、かな
り経験的な要素がありこれを計算機で判定させる指標を
得ることが困難であった。
一方エネルギー法は、負荷をアドミタンスマトリクスに
組み込むため、負荷の非線形特性が全く考慮できないと
いう欠点を有していた。この発明は上記のような従来の
ものの欠点を除去するためになされたもので、多機系を
等価的に2機におくことより、非線形負荷が考慮でき、
しかも演算時間の速い脱議予測のアルゴリズム及びそれ
を用いた安定判別装置を提供することを目的としている
以下この発明の一実施例について説明する。
第2図に予測手法のシークェンスを示す。まず‘2}の
部分で等価2機系におく。これは注目した脱調モードも
こおいて各発電機の位相が2グループに分かれるものと
し、一方をaグループ他方をbグループとすると、aグ
ループに属する発電機アングルをaa慣性定数をMa、
機械入力をPMaとして運動方程式はNねaa=PMa
−Pea ・・・・・・{9)但しとなる
同様にしてbに対してもMD 8b=PMb−Peb
......{1■となる。
次に とおき次の計算をおこなう。
W(81,82,81,82)=享Ma8≧十事M比0
隻−{PMa(8,一8,s)+PMb(82一8が} 十′aISPead81十′82Pebd82..・.
・・(12)但しaa,aa,8b,ab直後の値を用
いる。
次にステップ‘4}の部分でWc=W(0,0,8,c
,02c)を求める。
ここで0,c,82cは不安定平衡点の値である。そし
てステップ3で求めたWとステップ4で求めたWcを比
較しW>Wc →不安定 WミWc →安 定 と判定する。
この判定がステップ5でおこなわれる。以下本発明の原
理について説明する。
まず第3図のような多様系において脱調モ−ドがa群と
b群とに分かれるケースについて考える。
一般的に発電機を表わす運動方程式は、 Mjai=PMi−Pei(i=1〜n).・・.・・
(13) で表わされる。
今グループaに属する発電機をi=1〜m、残りをグル
ープbに属するとすると8i=81。
十aa(i=1〜m) ・・・・・・(1心8j=81
。十ab(i=m+1〜n).・・.・・(15) とおける。
但し8ioは発電機での初期アングルである。これを(
16)式に代入しとおくと NGaaニPM一Pea ・・・・・・(
19DM伍ab=PMb一Peb ・・・…
(20)となり等価的な2機系として表現できる。
この2機系について非線形負荷まで考慮した形のエネル
ギーを考える。
MT全Ma+Mb ……(21)8
09(NG8a+Mbab) ……(22)8・ea
a−8。
・・・・・・(23)02eob−8
……(2心とおくと非線形負荷まで考
慮したエネルギーはW(8・’82’8・,82)=鱗
州裏M比0室 −{PM1((0,一0,s)十PM2(82−a交)
}十′8,SPe,do,十′82SPe2d02
……(25)となる。
ここで0,s,02s,Pe,,Pe2はグループa,
bの安定平衡点および発電機出力をあらわす。又(25
)式の第1,2項は発電機の運動エネルギーを、第3,
4項は発電機回転子のポテンシャルエネルギーを、第5
,6項は電力系統に蓄えられるポテンシャルエネルギー
をあらわす。又(22),(23),(24)式より8
,=M伍(8a−8b2 82 =Ma(8b−8a
)MT MT….・・(26) となるのでこれらはaa−8bのみの関数となっている
第3図に示すように、電力系統6に蓄えられるポテンシ
ャルエネルギーは発電機グループa,bとの接続点の電
圧および相対位相aa−abが決まれば決定する。よっ
てPel,Pe2はたとえ電力系統に非線形負荷が入っ
ていてもoa−8bのみによって決まるので何通りかの
8a−abについて潮流計算を操り返すことによって第
4図のようにこれを求めることができ、(25)式は、
具体的に求めることができる。一方、不安定平衡点のエ
ネルギーは、グループaおよびbの発電機の加速度が等
しくなるという条件、すなわち8a=生と土筆 P単に
土塾=8b ..,...(27)Ma = Mbとな
るアングルを求めその点のポテンシャルエネルギーをW
cとすれば求まりと判定することが可能となる。
なお最後の判定の部分は、(25)式で求めたエネルギ
ーは保存されるので、(25)式のうち運動エネルギー
がすべてポテンシャルエネルギーにかわった時にもまだ
不安定平衡点に達していなければ系は安定となる、とい
うことにより導かれる。
なお上記実施例では、事故クリア後の脱調判定を行なっ
ているが、これは事故クリア後ある操作が加わった場合
の安定判別にも用いられる。又上記の手法は、負荷が非
線形でも用いることができる特色を有するが、線形でも
全く同様の予測手法が可能である。以上のように、この
発明によれば多機系を等価2機系で表現し、それに対し
て非線形負荷まで考慮したエネルギーを考えたので、非
線形負荷が入りしかも演算時間の速い脱調子側が行なえ
る効果がある。
又これを用いて系統の安定・不安定判定装置を構成する
ことも可能である。図面の簡単な説明第1図は、複数の
発電機が接続されている電力系統の一例を示す構成図、
第2図はこの発明に係る脱調子側方法の一実施例の説明
図、第3図はこの発明により、発電機を等価的な2機系
統とした場合の系統図、第4図はこの発明の方法で求め
られた電力系統のポテンシャルエネルギーを2機系銃の
夫々について示した説明図である。
図において、2,3,4,5はこの発明に係る脱調予測
方法の一実施例の各ステップ、6は電力系統を示す。
第1図 第2図 第3図 第4図

Claims (1)

    【特許請求の範囲】
  1. 1 電力系統に接続された複数の発電機について上記系
    統で生じた事故による脱調を予測する方法において、上
    記発電機を第1および第2のグループに分け、事故発生
    直後の各グループの発電機の位相、周波数、相対位相差
    から事故除去後の各グループの発電機の運動エネルギー
    、回転子のポテンシヤルエネルギーおよび電力系統のポ
    テンシヤルエネルギーの和を計算し、この値が予め求め
    た不安定平衡点のポテンシヤルエネルギーより大である
    ときは脱調、小であるときは脱調しないものと判別する
    ことを特徴とする脱調予測方法。
JP55018037A 1980-02-16 1980-02-16 脱調予測方法 Expired JPS6035894B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP55018037A JPS6035894B2 (ja) 1980-02-16 1980-02-16 脱調予測方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP55018037A JPS6035894B2 (ja) 1980-02-16 1980-02-16 脱調予測方法

Publications (2)

Publication Number Publication Date
JPS56115138A JPS56115138A (en) 1981-09-10
JPS6035894B2 true JPS6035894B2 (ja) 1985-08-17

Family

ID=11960469

Family Applications (1)

Application Number Title Priority Date Filing Date
JP55018037A Expired JPS6035894B2 (ja) 1980-02-16 1980-02-16 脱調予測方法

Country Status (1)

Country Link
JP (1) JPS6035894B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH041993U (ja) * 1990-04-23 1992-01-09

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH041993U (ja) * 1990-04-23 1992-01-09

Also Published As

Publication number Publication date
JPS56115138A (en) 1981-09-10

Similar Documents

Publication Publication Date Title
CN110070286B (zh) 电网多灾种耦合连锁故障分析方法及系统
Broadwater et al. Power flow analysis of unbalanced multiphase radial distribution systems
Rao et al. Nonlinear network reduction for distribution networks using the holomorphic embedding method
JP2017200275A (ja) 電力系統解析システム及び方法
CN120527962B (zh) 三相不平衡潮流线损分析方法、系统、设备及介质
Machowski et al. Damping of power swings by control of braking resistors
Pereira et al. Composite generation-transmission expansion planning
Ramirez-Gonzalez et al. Simplified fuzzy logic controller and its application as a power system stabilizer
Kasusky et al. Assessment of the SVC's effect on nonlinear instabilities and voltage collapse in electric power systems
CN112836165A (zh) 一种基于全纯嵌入的暂态稳定网络方程算法
Cao et al. Power system voltage small-disturbance stability studies based on the power flow equation
CN112234609B (zh) 基于电力系统随机动力学模型的非线性励磁控制方法
JP4116594B2 (ja) 電力系統の縮約モデル作成方法、装置およびプログラム
Murty et al. Design of decentralized variable structure stabilizers for multimachine power systems
CN114819482A (zh) 一种单机无穷大时变系统电压稳定性分析方法
Sakaguchi et al. Transient stability assessment for power system via Lur'e type Lyapunov function
Sharma et al. Transient stability assessment and synchronization of multimachine power system using kuramoto model
Pawellek et al. Simulation of dynamic power flows in the international grid control cooperation
CN116522043B (zh) 短路电流计算的方法及装置
Rajkumar et al. Nonlinear self-tuning control of a flexible AC transmission system
Padiyar et al. A novel structure preserving energy function for direct stability evaluation of power systems with known modes of instability
CN121643129B (zh) 一种用于多vsg并联系统的暂态稳定性评估方法及装置
JP2703091B2 (ja) 電力系統模擬装置
Laus et al. SIMO and MIMO PHIL methods for distributed generation in LV networks
Alkadhim Simulation of Multi-Machine Transient Stability