JP3459158B2 - 電力系統シミュレーション方法 - Google Patents
電力系統シミュレーション方法Info
- Publication number
- JP3459158B2 JP3459158B2 JP33734496A JP33734496A JP3459158B2 JP 3459158 B2 JP3459158 B2 JP 3459158B2 JP 33734496 A JP33734496 A JP 33734496A JP 33734496 A JP33734496 A JP 33734496A JP 3459158 B2 JP3459158 B2 JP 3459158B2
- Authority
- JP
- Japan
- Prior art keywords
- power system
- electric
- analysis
- mechanical
- analysis process
- 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 - Fee Related
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Description
【０００１】
【発明の属する技術分野】この発明は、電力系統過渡現
象を電気系と機械系とに分離して正しい解析結果を得る
ことのできる電力系統シミュレーション方法に関するも
のである。
象を電気系と機械系とに分離して正しい解析結果を得る
ことのできる電力系統シミュレーション方法に関するも
のである。
【０００２】
【従来の技術】電力系統過渡現象を解析するための従来
一般に行われている手法は、「電力系統過渡解析論」
（関根泰次、オーム社、１９８４）に示されている。こ
の手法の特徴は電力系統の交流電気回路ベクトル変数を
実数に変換する。このため、過渡状態同期発電機交流電
気回路方程式のみで電気系の解が得られず、電力系統電
気系の過渡現象を解明できない。そこで電気系定数にい
くつかの補正設定を施し、電力系統電気系と電力系統機
械系の連立方程式で電気系および機械系の過渡現象を解
明する。
一般に行われている手法は、「電力系統過渡解析論」
（関根泰次、オーム社、１９８４）に示されている。こ
の手法の特徴は電力系統の交流電気回路ベクトル変数を
実数に変換する。このため、過渡状態同期発電機交流電
気回路方程式のみで電気系の解が得られず、電力系統電
気系の過渡現象を解明できない。そこで電気系定数にい
くつかの補正設定を施し、電力系統電気系と電力系統機
械系の連立方程式で電気系および機械系の過渡現象を解
明する。
【０００３】図１２は、この「電力系統過渡解析論」に
示されている電力系統過渡現象を解析するための従来の
手法をもとにした電力系統シミュレーション方法の構成
を示すフローチャートである。このフローチャートによ
れば、先ず、電力系統の初期潮流を計算し（ステップＳ
Ｔ１）、さらに電力系統の発電機諸量の初期値を計算す
る（ステップＳＴ２）。そして、電力系統の状態変化を
計算し（ステップＳＴ３）、電力系統各部分の電圧、電
流を計算し（ステップＳＴ４）、電力系統微分方程式を
計算し（ステップＳＴ５）、さらに前記計算した電力系
統微分方程式で次のステップを計算し（ステップＳＴ
６）、電力系統シミュレーションが終了するまでステッ
プＳＴ４からステップＳＴ６までの処理を繰り返す。
示されている電力系統過渡現象を解析するための従来の
手法をもとにした電力系統シミュレーション方法の構成
を示すフローチャートである。このフローチャートによ
れば、先ず、電力系統の初期潮流を計算し（ステップＳ
Ｔ１）、さらに電力系統の発電機諸量の初期値を計算す
る（ステップＳＴ２）。そして、電力系統の状態変化を
計算し（ステップＳＴ３）、電力系統各部分の電圧、電
流を計算し（ステップＳＴ４）、電力系統微分方程式を
計算し（ステップＳＴ５）、さらに前記計算した電力系
統微分方程式で次のステップを計算し（ステップＳＴ
６）、電力系統シミュレーションが終了するまでステッ
プＳＴ４からステップＳＴ６までの処理を繰り返す。
【０００４】
【発明が解決しようとする課題】従来の電力系統シミュ
レーション方法は以上のように構成されていたので、時
定数の異なる電力系統電気系と電力系統機械系とを分離
することができず、電力系統過渡現象を正しく解析する
ことが困難であり、電力系統の送電限界を求める場合な
どに障害となる課題があった。また、電力系統電気系解
析において電気回路の解析解が得られないため、電気系
の過渡現象が不明確な状態で電力系統シミュレーション
を行うため、電力系統の最適な高速制御を行うための正
しいシミュレーションが実現できない課題があった。
レーション方法は以上のように構成されていたので、時
定数の異なる電力系統電気系と電力系統機械系とを分離
することができず、電力系統過渡現象を正しく解析する
ことが困難であり、電力系統の送電限界を求める場合な
どに障害となる課題があった。また、電力系統電気系解
析において電気回路の解析解が得られないため、電気系
の過渡現象が不明確な状態で電力系統シミュレーション
を行うため、電力系統の最適な高速制御を行うための正
しいシミュレーションが実現できない課題があった。
【０００５】この発明は、上記のような課題を解決する
ためになされたもので、時定数の異なる電力系統電気系
と電力系統機械系とを分離し電力系統過渡現象を正しく
解析し、さらに電力系統の最適な高速制御を行うための
シミュレーションを実現する電力系統シミュレーション
方法を得ることを目的とする。
ためになされたもので、時定数の異なる電力系統電気系
と電力系統機械系とを分離し電力系統過渡現象を正しく
解析し、さらに電力系統の最適な高速制御を行うための
シミュレーションを実現する電力系統シミュレーション
方法を得ることを目的とする。
【０００６】
【課題を解決するための手段】請求項１記載の発明に係
る電力系統シミュレーション方法は、電力系統の電気系
過渡現象の解析にスパイラルベクトル理論を適用し、
る電力系統シミュレーション方法は、電力系統の電気系
過渡現象の解析にスパイラルベクトル理論を適用し、
【数３】 により、前記電力系統の電流一般解の解析を行う電力系
統電気系解析過程と、該電力系統電気系解析過程による
解析結果をもとに、
統電気系解析過程と、該電力系統電気系解析過程による
解析結果をもとに、
【数４】 により、前記電力系統の機械系の解析を行う電力系統機
械系解析過程と、該電力系統機械系解析過程による解析
結果をもとに、前記電力系統電気系解析過程で得られた
解析結果の各定数を修正して前記電力系統機械系解析過
程による解析結果を前記電力系統電気系解析過程へ反映
させる反映過程とを備えるようにしたものである。
械系解析過程と、該電力系統機械系解析過程による解析
結果をもとに、前記電力系統電気系解析過程で得られた
解析結果の各定数を修正して前記電力系統機械系解析過
程による解析結果を前記電力系統電気系解析過程へ反映
させる反映過程とを備えるようにしたものである。
【０００７】
請求項２記載の発明に係る電力系統シミュ
レーション方法は、シミュレートしようとする状況につ
いての電力系統の初期状態に応じた定数を設定する初期
条件設定過程を有し、該初期条件設定過程で前記定数が
設定された前記初期状態に応じた前記状況発生時の過渡
現象について前記電力系統の電気系の解析を行う電力系
統電気系解析過程と、該電力系統電気系解析過程による
前記状況発生時の解析結果をもとに前記状況発生時の前
記電力系統の機械系の解析を行う電力系統機械系解析過
程とを備えるようにしたものである。
レーション方法は、シミュレートしようとする状況につ
いての電力系統の初期状態に応じた定数を設定する初期
条件設定過程を有し、該初期条件設定過程で前記定数が
設定された前記初期状態に応じた前記状況発生時の過渡
現象について前記電力系統の電気系の解析を行う電力系
統電気系解析過程と、該電力系統電気系解析過程による
前記状況発生時の解析結果をもとに前記状況発生時の前
記電力系統の機械系の解析を行う電力系統機械系解析過
程とを備えるようにしたものである。
【０００８】
請求項３記載の発明に係る電力系統シミュ
レーション方法は、電力系統で発生する負荷変動の発生
条件初期状態に応じた定数を設定する初期条件設定過程
を有し、該初期条件設定過程で前記定数が設定された負
荷変動発生の際の条件である初期状態に応じた負荷変動
発生時の過渡現象について前記電力系統の電気系の解析
を行う電力系統電気系解析過程と、該電力系統電気系解
析過程による前記負荷変動発生時の解析結果をもとに前
記負荷変動発生時の前記電力系統の機械系の解析を行う
電力系統機械系解析過程とを備えるようにしたものであ
る。
レーション方法は、電力系統で発生する負荷変動の発生
条件初期状態に応じた定数を設定する初期条件設定過程
を有し、該初期条件設定過程で前記定数が設定された負
荷変動発生の際の条件である初期状態に応じた負荷変動
発生時の過渡現象について前記電力系統の電気系の解析
を行う電力系統電気系解析過程と、該電力系統電気系解
析過程による前記負荷変動発生時の解析結果をもとに前
記負荷変動発生時の前記電力系統の機械系の解析を行う
電力系統機械系解析過程とを備えるようにしたものであ
る。
【０００９】
請求項４記載の発明に係る電力系統シミュ
レーション方法は、負荷変動発生以降の過渡現象につい
ての電力系統の電気系の解析結果から得られた電力系統
同期機の発生出力を求める電力系統同期機発生出力演算
過程と、前記電力系統同期機発生出力演算過程で求めた
前記電力系統同期機の発生出力をもとに電力系統機械系
の状態変化量を求める状態変化量演算過程と、前記状態
変化量演算過程で求めた状態変化量に応じて電力系統電
気系解析過程で得られた解析結果の各定数を修正するこ
とで、電力系統機械系解析過程による解析結果を前記電
力系統電気系解析過程へ反映させる反映過程とを備える
ようにしたものである。
レーション方法は、負荷変動発生以降の過渡現象につい
ての電力系統の電気系の解析結果から得られた電力系統
同期機の発生出力を求める電力系統同期機発生出力演算
過程と、前記電力系統同期機発生出力演算過程で求めた
前記電力系統同期機の発生出力をもとに電力系統機械系
の状態変化量を求める状態変化量演算過程と、前記状態
変化量演算過程で求めた状態変化量に応じて電力系統電
気系解析過程で得られた解析結果の各定数を修正するこ
とで、電力系統機械系解析過程による解析結果を前記電
力系統電気系解析過程へ反映させる反映過程とを備える
ようにしたものである。
【００１０】
請求項５記載の発明に係る電力系統シミュ
レーション方法は、電力系統における事故の発生する条
件である初期状態に応じた定数を設定する初期条件設定
過程を有し、前記初期条件設定過程で前記定数が設定さ
れた初期状態に応じた事故発生時の過渡現象について前
記電力系統の電気系の解析を電力系統電気系解析過程で
行い、前記電力系統電気系解析過程による前記事故発生
時の解析結果をもとに前記事故発生時の前記電力系統の
機械系の解析を電力系統機械系解析過程で行うようにし
たものである。
レーション方法は、電力系統における事故の発生する条
件である初期状態に応じた定数を設定する初期条件設定
過程を有し、前記初期条件設定過程で前記定数が設定さ
れた初期状態に応じた事故発生時の過渡現象について前
記電力系統の電気系の解析を電力系統電気系解析過程で
行い、前記電力系統電気系解析過程による前記事故発生
時の解析結果をもとに前記事故発生時の前記電力系統の
機械系の解析を電力系統機械系解析過程で行うようにし
たものである。
【００１１】
請求項６記載の発明に係る電力系統シミュ
レーション方法は、事故発生以降の過渡現象についての
電力系統の電気系の解析結果から電力系統同期機の発生
出力を求める電力系統同期機発生出力演算過程と、該電
力系統同期機発生出力演算過程で求めた前記電力系統同
期機の発生出力をもとに電力系統機械系の状態変化量を
演算し求める状態変化量演算過程と、該状態変化量演算
過程で求めた状態変化量に応じて電力系統電気系解析過
程で得られた解析結果の各定数を修正することで、前記
電力系統機械系解析過程による解析結果を前記電力系統
電気系解析過程へ反映させる反映過程とを備えるように
したものである。
レーション方法は、事故発生以降の過渡現象についての
電力系統の電気系の解析結果から電力系統同期機の発生
出力を求める電力系統同期機発生出力演算過程と、該電
力系統同期機発生出力演算過程で求めた前記電力系統同
期機の発生出力をもとに電力系統機械系の状態変化量を
演算し求める状態変化量演算過程と、該状態変化量演算
過程で求めた状態変化量に応じて電力系統電気系解析過
程で得られた解析結果の各定数を修正することで、前記
電力系統機械系解析過程による解析結果を前記電力系統
電気系解析過程へ反映させる反映過程とを備えるように
したものである。
【００１２】
請求項７記載の発明に係る電力系統シミュ
レーション方法は、状態変化量演算過程で求めた電力系
統機械系の状態変化量を所定の閾値と比較する比較過程
と、前記状態変化量演算過程で求めた状態変化量が前記
所定の閾値を越えているとの比較結果が前記比較過程で
得られると電力系統電気系解析過程で得られた解析結果
の各定数を修正し、前記電力系統機械系解析過程による
解析結果を前記電力系統電気系解析過程へ反映過程が反
映させるようにしたものである。
レーション方法は、状態変化量演算過程で求めた電力系
統機械系の状態変化量を所定の閾値と比較する比較過程
と、前記状態変化量演算過程で求めた状態変化量が前記
所定の閾値を越えているとの比較結果が前記比較過程で
得られると電力系統電気系解析過程で得られた解析結果
の各定数を修正し、前記電力系統機械系解析過程による
解析結果を前記電力系統電気系解析過程へ反映過程が反
映させるようにしたものである。
【００１３】
請求項８記載の発明に係る電力系統シミュ
レーション方法は、状態変化量演算過程で求めた状態変
化量が所定の閾値を越えていないとの比較結果が前記比
較過程で得られると、電力系統同期機発生出力演算過程
での電力系統同期機の発生出力の演算を行い、この電力
系統同期機の発生出力をもとに前記状態変化量演算過程
で求めた状態変化量と前記所定の閾値との比較を前記比
較過程で行う再比較過程を備えるようにしたものであ
る。
レーション方法は、状態変化量演算過程で求めた状態変
化量が所定の閾値を越えていないとの比較結果が前記比
較過程で得られると、電力系統同期機発生出力演算過程
での電力系統同期機の発生出力の演算を行い、この電力
系統同期機の発生出力をもとに前記状態変化量演算過程
で求めた状態変化量と前記所定の閾値との比較を前記比
較過程で行う再比較過程を備えるようにしたものであ
る。
【００１４】
請求項９記載の発明に係る電力系統シミュ
レーション方法は、電力系統同期機発生出力演算過程で
求めた電力系統同期機の発生出力をもとに電力系統機械
系の位相角である状態変化量を電力系統機械系解析過程
の状態変化量演算過程が求めるようにしたものである。
レーション方法は、電力系統同期機発生出力演算過程で
求めた電力系統同期機の発生出力をもとに電力系統機械
系の位相角である状態変化量を電力系統機械系解析過程
の状態変化量演算過程が求めるようにしたものである。
【００１５】
請求項１０記載の発明に係る電力系統シミ
ュレーション方法は、事故発生以降および線路開放以降
の系統状態の急変の有無について判定する系統状態急変
判定過程を備えるようにしたものである。
ュレーション方法は、事故発生以降および線路開放以降
の系統状態の急変の有無について判定する系統状態急変
判定過程を備えるようにしたものである。
【００１６】
請求項１１記載の発明に係る電力系統シミ
ュレーション方法は、系統状態急変判定過程により系統
状態の急変があったことが判定されたときの電力系統電
気系の過渡現象について解析を行う系統状態急変時電気
系解析過程と、該系統状態急変時電気系解析過程の解析
結果をもとに前記電力系統が急変したときの前記電力系
統の機械系の解析を行う電力系統機械系解析過程とを備
えるようにしたものである。
ュレーション方法は、系統状態急変判定過程により系統
状態の急変があったことが判定されたときの電力系統電
気系の過渡現象について解析を行う系統状態急変時電気
系解析過程と、該系統状態急変時電気系解析過程の解析
結果をもとに前記電力系統が急変したときの前記電力系
統の機械系の解析を行う電力系統機械系解析過程とを備
えるようにしたものである。
【００１７】
請求項１２記載の発明に係る電力系統シミ
ュレーション方法は、初期条件設定過程で定数が設定さ
れた初期状態に応じた事故発生時の過渡現象について過
渡直流電流を含む電力系統の電気系の解析を行う電力系
統電気系解析過程を備えるようにしたものである。
ュレーション方法は、初期条件設定過程で定数が設定さ
れた初期状態に応じた事故発生時の過渡現象について過
渡直流電流を含む電力系統の電気系の解析を行う電力系
統電気系解析過程を備えるようにしたものである。
【００１８】
請求項１３記載の発明に係る電力系統シミ
ュレーション方法は、複数の電力系統同期機を有した電
力系統の電力系統電気系解析過程および電力系統機械系
解析過程における演算にそれぞれ並列演算ノードを設定
し、前記電力系統電気系解析過程による電力系統電気系
の解析および前記電力系統機械系解析過程による電力系
統機械系の解析を並列的に行うようにしたものである。
ュレーション方法は、複数の電力系統同期機を有した電
力系統の電力系統電気系解析過程および電力系統機械系
解析過程における演算にそれぞれ並列演算ノードを設定
し、前記電力系統電気系解析過程による電力系統電気系
の解析および前記電力系統機械系解析過程による電力系
統機械系の解析を並列的に行うようにしたものである。
【００１９】
【発明の実施の形態】以下、この発明の実施の一形態を
説明する。 実施の形態１． この実施の形態１の電力系統シミュレーション方法は、
同期発電機は、事故、負荷変動などに対し電気的出力に
より対応していることから、電力系統電気系の解析結果
を電力系統機械系に反映させ、前記電力系統機械系へ反
映させた前記同期発電機電気的出力などの反映結果（機
械系解析結果）をさらに前記電力系統電気系へフィード
バックするなどの手法を用い、前記電力系統電気系と前
記電力系統機械系とを別々に解析するものである。
説明する。 実施の形態１． この実施の形態１の電力系統シミュレーション方法は、
同期発電機は、事故、負荷変動などに対し電気的出力に
より対応していることから、電力系統電気系の解析結果
を電力系統機械系に反映させ、前記電力系統機械系へ反
映させた前記同期発電機電気的出力などの反映結果（機
械系解析結果）をさらに前記電力系統電気系へフィード
バックするなどの手法を用い、前記電力系統電気系と前
記電力系統機械系とを別々に解析するものである。
【００２０】
図１は、この実施の形態１の電力系統事故
発生時の電力系統シミュレーション方法の構成を示すフ
ローチャートである。このフローチャートによれば、先
ず電力系統の潮流計算を行い（ステップＳＴ１１）、電
力系統事故発生条件初期状態による定数を求める（ステ
ップＳＴ１２）（初期条件設定過程）。スパイラルベク
トル理論による電力系統電気系解析解を求め（ステップ
ＳＴ１３）（電力系統電気系解析過程）、さらに電力系
統それぞれの同期発電機の機械系微分方程式を解くこと
で同期機電気系の電気的出力より電力系統機械系の位相
角変化分を演算して求める（ステップＳＴ１４）（電力
系統機械系解析過程，電力系統同期機発生出力演算過
程，状態変化量演算過程）。この演算は数値積分法によ
り求める。そして、線路遮断器の操作により線路を開放
するなどして電力系統の事故をクリアした後に系統状態
急変が発生したか否かを判定し（ステップＳＴ１５）
（系統状態急変判定過程）、電力系統シミュレーション
の終了時間経過を判定する（ステップＳＴ１６）。
発生時の電力系統シミュレーション方法の構成を示すフ
ローチャートである。このフローチャートによれば、先
ず電力系統の潮流計算を行い（ステップＳＴ１１）、電
力系統事故発生条件初期状態による定数を求める（ステ
ップＳＴ１２）（初期条件設定過程）。スパイラルベク
トル理論による電力系統電気系解析解を求め（ステップ
ＳＴ１３）（電力系統電気系解析過程）、さらに電力系
統それぞれの同期発電機の機械系微分方程式を解くこと
で同期機電気系の電気的出力より電力系統機械系の位相
角変化分を演算して求める（ステップＳＴ１４）（電力
系統機械系解析過程，電力系統同期機発生出力演算過
程，状態変化量演算過程）。この演算は数値積分法によ
り求める。そして、線路遮断器の操作により線路を開放
するなどして電力系統の事故をクリアした後に系統状態
急変が発生したか否かを判定し（ステップＳＴ１５）
（系統状態急変判定過程）、電力系統シミュレーション
の終了時間経過を判定する（ステップＳＴ１６）。
【００２１】
ステップＳＴ１６における判定の結果、電
力系統シミュレーションの終了時間を経過していなけれ
ば電力系統機械系の同期発電機の位相角の変化とその位
相角の変化についての閾値とを比較する（ステップＳＴ
１７）（比較過程）。また、ステップＳＴ１６において
電力系統シミュレーションの終了時間を経過したときに
は電力系統シミュレーションを終了する。
力系統シミュレーションの終了時間を経過していなけれ
ば電力系統機械系の同期発電機の位相角の変化とその位
相角の変化についての閾値とを比較する（ステップＳＴ
１７）（比較過程）。また、ステップＳＴ１６において
電力系統シミュレーションの終了時間を経過したときに
は電力系統シミュレーションを終了する。
【００２２】
ステップＳＴ１７における閾値は、電力系
統シミュレーションを行う際の計算精度を追及する場合
には０．００１ラジアン以下、計算速度を追及する場合
には０．１ラジアン以上が好ましい。また、計算精度や
計算速度のバランスを考慮する場合には０．０１ラジア
ン程度が良いが、電力系統シミュレーションの目的によ
りこの閾値を変えるようにするのが望ましい。
統シミュレーションを行う際の計算精度を追及する場合
には０．００１ラジアン以下、計算速度を追及する場合
には０．１ラジアン以上が好ましい。また、計算精度や
計算速度のバランスを考慮する場合には０．０１ラジア
ン程度が良いが、電力系統シミュレーションの目的によ
りこの閾値を変えるようにするのが望ましい。
【００２３】
ステップＳＴ１７における判定の結果、前
記位相角の変化が前記閾値を越える場合には、スパイラ
ルベクトル理論により電力系統電気系の解析解を修正す
る（ステップＳＴ１９）（反映過程）。一方、ステップ
ＳＴ１７における判定の結果、前記位相角の変化が前記
閾値を越えない場合には、ステップＳＴ１４へ戻り、ス
テップＳＴ１４以降の処理を実行する（再比較過程）。
そして、これらステップＳＴ１４からステップＳＴ１７
までの処理を、電力系統シミュレーションの終了時間経
過まで繰り返す。以上の説明におけるスパイラルベクト
ルとは、図２に示すように複素平面上で反時計回りで減
衰しながら回転するベクトル
記位相角の変化が前記閾値を越える場合には、スパイラ
ルベクトル理論により電力系統電気系の解析解を修正す
る（ステップＳＴ１９）（反映過程）。一方、ステップ
ＳＴ１７における判定の結果、前記位相角の変化が前記
閾値を越えない場合には、ステップＳＴ１４へ戻り、ス
テップＳＴ１４以降の処理を実行する（再比較過程）。
そして、これらステップＳＴ１４からステップＳＴ１７
までの処理を、電力系統シミュレーションの終了時間経
過まで繰り返す。以上の説明におけるスパイラルベクト
ルとは、図２に示すように複素平面上で反時計回りで減
衰しながら回転するベクトル
【数５】 をいう。
【００２４】
次に、図３に示す電力系統図をもとに図１
のフローチャートに示す動作について説明する。図３に
おいて、Ｒ１は同期発電機の電機子抵抗、Ｒは線路抵
抗、Ｌｓは同期発電機の同期インダクタンス、Ｌは線路
インダクタンス、ＣＢ１，ＣＢ２，ＣＢ３，ＣＢ４は線
路遮断器、Ｆは事故点を示している。事故点Ｆにおける
地絡事故を想定し、一定時間後に線路の線路遮断器ＣＢ
１，ＣＢ２を動作させ、線路を開放するというシミュレ
ーションを行う。
のフローチャートに示す動作について説明する。図３に
おいて、Ｒ１は同期発電機の電機子抵抗、Ｒは線路抵
抗、Ｌｓは同期発電機の同期インダクタンス、Ｌは線路
インダクタンス、ＣＢ１，ＣＢ２，ＣＢ３，ＣＢ４は線
路遮断器、Ｆは事故点を示している。事故点Ｆにおける
地絡事故を想定し、一定時間後に線路の線路遮断器ＣＢ
１，ＣＢ２を動作させ、線路を開放するというシミュレ
ーションを行う。
【００２５】
同期発電機の内部電圧をｅ1 、同期発電機
の端子電圧をｖ1 、端子電流をｉ1（無限大母線電圧を
ｖb ）と仮定すると、次式（１）および次式（２）が成
立する。ここで、全ての変数がスパイラルベクトルであ
る。
の端子電圧をｖ1 、端子電流をｉ1（無限大母線電圧を
ｖb ）と仮定すると、次式（１）および次式（２）が成
立する。ここで、全ての変数がスパイラルベクトルであ
る。
【００２６】
ｅ1 ＝ｖ1 ＋ｉ1 Ｒ1 ＋Ｌs ｐｉ1 （１）
ｖ1 ＝ｖb ＋ｉ1 （Ｒ／２）＋（Ｌ／２）ｐｉ1 （２）
【００２７】
電力系統の定常状態においては、スパイラ
ルベクトルの大きさと回転速度は一定であるので、次式
（３）および次式（４）が成立する。
ルベクトルの大きさと回転速度は一定であるので、次式
（３）および次式（４）が成立する。
【００２８】
【数６】 ここで、Ｅ1 は電力系統同期機内部電圧定常ベクトル、
Ｖ1 は電力系統同期機端子電圧定常ベクトル、Ｖb は電
力系統無限大母線定常ベクトル、Ｉ1 は電力系統同期機
端子電流定常ベクトル、ωは電力系統定常角速度であ
り、式（３）および式（４）により図４に示すベクトル
図が成立する。
Ｖ1 は電力系統同期機端子電圧定常ベクトル、Ｖb は電
力系統無限大母線定常ベクトル、Ｉ1 は電力系統同期機
端子電流定常ベクトル、ωは電力系統定常角速度であ
り、式（３）および式（４）により図４に示すベクトル
図が成立する。
【００２９】
事故発生後および線路遮断器による線路の
開放で事故をクリアした後の電力系統の過渡状態では、
スパイラルベクトル理論による電流一般解ｉ1sは次式
（５）の通りである。
開放で事故をクリアした後の電力系統の過渡状態では、
スパイラルベクトル理論による電流一般解ｉ1sは次式
（５）の通りである。
【００３０】
【数７】
【００３１】ここで、λは減衰定数、Ｉ1sは電流定常
解、Ａは初期条件により決定する任意定数である。この
ようにして事故中（事故発生後および線路遮断器による
線路の開放で事故をクリアした後）における電力系統電
気系の過渡現象は解析解より解析することが可能であ
る。この式（５）は、事故発生後および線路遮断器によ
る線路の開放で事故をクリアした後に電力系統が変った
ときには相応の解析解となる（ステップＳＴ１８）。そ
して、式（５）と前記初期条件により次式（６）に示す
同期機電気的出力Ｐeを求めることができる。
解、Ａは初期条件により決定する任意定数である。この
ようにして事故中（事故発生後および線路遮断器による
線路の開放で事故をクリアした後）における電力系統電
気系の過渡現象は解析解より解析することが可能であ
る。この式（５）は、事故発生後および線路遮断器によ
る線路の開放で事故をクリアした後に電力系統が変った
ときには相応の解析解となる（ステップＳＴ１８）。そ
して、式（５）と前記初期条件により次式（６）に示す
同期機電気的出力Ｐeを求めることができる。
【００３２】
【数８】
【００３３】電力系統に複数の同期発電機が存在してい
るときには、このようにして求めた各同期発電機の同期
機電気的出力Ｐei（i は１，２，３・・・Ｎであり同期
発電機の数）を次式（７）へ代入することで、各同期発
電機の電力系統機械系の１ステップにおける位相角変化
分Δδi を求める（ステップＳＴ１４）。
るときには、このようにして求めた各同期発電機の同期
機電気的出力Ｐei（i は１，２，３・・・Ｎであり同期
発電機の数）を次式（７）へ代入することで、各同期発
電機の電力系統機械系の１ステップにおける位相角変化
分Δδi を求める（ステップＳＴ１４）。
【００３４】
【数９】
【００３５】ここで、δは同期発電機の位相角、Ｐm は
同期機機械的入力、Ｐe は同期機電気的出力、Ｍは同期
機慣性定数、ωは基本角速度である。そして、事故クリ
ア、負荷脱落などによる系統状態の急変があるか否かを
判定する（ステップＳＴ１５）。この判定の結果、系統
状態の急変があったときには新しい系統についての電力
系統電気系の解析的な解を再度スパイラルベクトル理論
により得て（ステップＳＴ１８）、ステップＳＴ１４へ
戻る。一方、ステップＳＴ１５において系統状態の急変
がないと判定されると、系統の構成は変わっていないの
で、ステップＳＴ１６へ進んでシミュレーション終了時
間が経過したか否かを判定する。この結果、シミュレー
ション終了時間が経過していなければ、シミュレーショ
ンを継続して電力系統機械系の各同期発電機の位相角の
変化Δδi がそれぞれの同期発電機のある一定の閾値δ
0i（i は１，２，３・・・Ｎであり同期発電機の数）を
越えているか否かを判定する（ステップＳＴ１７）。ス
テップＳＴ１６において、シミュレーション終了時間を
経過していれば電力系統シミュレーションを終了させ
る。
同期機機械的入力、Ｐe は同期機電気的出力、Ｍは同期
機慣性定数、ωは基本角速度である。そして、事故クリ
ア、負荷脱落などによる系統状態の急変があるか否かを
判定する（ステップＳＴ１５）。この判定の結果、系統
状態の急変があったときには新しい系統についての電力
系統電気系の解析的な解を再度スパイラルベクトル理論
により得て（ステップＳＴ１８）、ステップＳＴ１４へ
戻る。一方、ステップＳＴ１５において系統状態の急変
がないと判定されると、系統の構成は変わっていないの
で、ステップＳＴ１６へ進んでシミュレーション終了時
間が経過したか否かを判定する。この結果、シミュレー
ション終了時間が経過していなければ、シミュレーショ
ンを継続して電力系統機械系の各同期発電機の位相角の
変化Δδi がそれぞれの同期発電機のある一定の閾値δ
0i（i は１，２，３・・・Ｎであり同期発電機の数）を
越えているか否かを判定する（ステップＳＴ１７）。ス
テップＳＴ１６において、シミュレーション終了時間を
経過していれば電力系統シミュレーションを終了させ
る。
【００３６】
ステップＳＴ１７において、前記位相角の
変化Δδi が前記一定の閾値δ0iを越えている同期発電
機があれば、当該同期発電機の前記位相角の変化Δδi
を電力系統電気系に反映させ、電力系統電気系の修正処
理を行い、次式（５）で示される電力系統電気系の解析
解の各定数を修正し、系統状態の急変のない同じ系統の
当該同期発電機について新たな電気系の解析的な解をス
パイラルベクトル理論により得る。そして、機械系の位
相角の変化δi を零として（ステップＳＴ１９）、ステ
ップＳＴ１４へ戻る。一方、ステップＳＴ１７において
前記位相角の変化Δδi が前記一定の閾値δ0iを越えて
いる同期発電機がなければ、電力系統機械系の位相角の
変化分が小さく、現在の電力系統電気系の解析結果の精
度が充分高いとみなして、ステップＳＴ１９で行うよう
な電力系統電気系の修正処理を行わず、電力系統機械系
の位相角の変化分の演算処理であるステップＳＴ１４へ
戻る。
変化Δδi が前記一定の閾値δ0iを越えている同期発電
機があれば、当該同期発電機の前記位相角の変化Δδi
を電力系統電気系に反映させ、電力系統電気系の修正処
理を行い、次式（５）で示される電力系統電気系の解析
解の各定数を修正し、系統状態の急変のない同じ系統の
当該同期発電機について新たな電気系の解析的な解をス
パイラルベクトル理論により得る。そして、機械系の位
相角の変化δi を零として（ステップＳＴ１９）、ステ
ップＳＴ１４へ戻る。一方、ステップＳＴ１７において
前記位相角の変化Δδi が前記一定の閾値δ0iを越えて
いる同期発電機がなければ、電力系統機械系の位相角の
変化分が小さく、現在の電力系統電気系の解析結果の精
度が充分高いとみなして、ステップＳＴ１９で行うよう
な電力系統電気系の修正処理を行わず、電力系統機械系
の位相角の変化分の演算処理であるステップＳＴ１４へ
戻る。
【００３７】
以上説明したように、この実施の形態１に
よれば、スパイラルベクトルという現実に即した変数を
用い、電力系統機械系と時定数の異なる電力系統電気系
の過渡現象を解明し、その解明した電力系統電気系の過
渡現象をもとに機械系の過渡現象を解明し、さらにその
解明した機械系の過渡現象をもとに前記電力系統電気系
の解析解の各定数を修正して新しい解析解を得るという
手法を用いるため、電力系統過渡現象について現実の現
象に忠実な正しい解析結果が得られる効果がある。
よれば、スパイラルベクトルという現実に即した変数を
用い、電力系統機械系と時定数の異なる電力系統電気系
の過渡現象を解明し、その解明した電力系統電気系の過
渡現象をもとに機械系の過渡現象を解明し、さらにその
解明した機械系の過渡現象をもとに前記電力系統電気系
の解析解の各定数を修正して新しい解析解を得るという
手法を用いるため、電力系統過渡現象について現実の現
象に忠実な正しい解析結果が得られる効果がある。
【００３８】
実施の形態２．
この実施の形態２では、電力系統事故直後における過渡
直流電流を再現したシミュレートを行うことで、電力系
統事故直後におけるより現実的な電力系統シミュレーシ
ョンの可能な電力系統シミュレータを実現するための電
力系統シミュレーション方法について説明する。図５
は、この実施の形態２の電力系統過渡状態における動揺
モードを示す説明図である。この説明図に示されている
ように電力系統過渡状態においては全ての動揺モードを
大きく２種類に分類することが可能である。第１の動揺
モードは、電力系統電気系モードであり時定数は数ミリ
秒オーダである。第２の動揺モードは、電力系統機械系
モードであり時定数は数秒オーダである。
直流電流を再現したシミュレートを行うことで、電力系
統事故直後におけるより現実的な電力系統シミュレーシ
ョンの可能な電力系統シミュレータを実現するための電
力系統シミュレーション方法について説明する。図５
は、この実施の形態２の電力系統過渡状態における動揺
モードを示す説明図である。この説明図に示されている
ように電力系統過渡状態においては全ての動揺モードを
大きく２種類に分類することが可能である。第１の動揺
モードは、電力系統電気系モードであり時定数は数ミリ
秒オーダである。第２の動揺モードは、電力系統機械系
モードであり時定数は数秒オーダである。
【００３９】
図６は、これら電力系統電気系動揺モード
と電力系統機械系動揺モードのモード特性を示す説明図
である。この説明図に示すように、電力系統電気系動揺
モードは電力系統電気回路の抵抗、インダクタンス、同
期機の突極性などと関係しているが、電力系統同期機の
慣性定数、潮流量などとの関連は希薄である。一方、電
力系統機械系動揺モードは、電力系統同期機の慣性定
数、潮流量などと密接に関係している。なお、電力系統
電気系動揺モードに対しては磁束連続定理により過渡状
態と定常状態との電流の連続性を保つことが必要であ
り、また電力系統機械系動揺モードに対しては、エネル
ギー不変の観点から位相角の連続性が必要である。
と電力系統機械系動揺モードのモード特性を示す説明図
である。この説明図に示すように、電力系統電気系動揺
モードは電力系統電気回路の抵抗、インダクタンス、同
期機の突極性などと関係しているが、電力系統同期機の
慣性定数、潮流量などとの関連は希薄である。一方、電
力系統機械系動揺モードは、電力系統同期機の慣性定
数、潮流量などと密接に関係している。なお、電力系統
電気系動揺モードに対しては磁束連続定理により過渡状
態と定常状態との電流の連続性を保つことが必要であ
り、また電力系統機械系動揺モードに対しては、エネル
ギー不変の観点から位相角の連続性が必要である。
【００４０】
図５，図６に示す電力系統電気系動揺モー
ドと電力系統機械系動揺モードの比較において、両者間
では時定数の値が大きく異なっている。この特性を利用
して機械系変数の変化が無視できる時間領域については
この値を一定と仮定し、スパイラルベクトル理論に基づ
き電力系統電気系の解析解を求める。そして、この解を
用いて機械系の解析解を求める。さらに一定時間後の電
力系統機械系の変化を電力系統電気系に反映させ、新し
い解析的な解が求められる。図３に示す電力系統のＦ点
の事故直後において、電力系統電気系モードは次式
（８）により示される。
ドと電力系統機械系動揺モードの比較において、両者間
では時定数の値が大きく異なっている。この特性を利用
して機械系変数の変化が無視できる時間領域については
この値を一定と仮定し、スパイラルベクトル理論に基づ
き電力系統電気系の解析解を求める。そして、この解を
用いて機械系の解析解を求める。さらに一定時間後の電
力系統機械系の変化を電力系統電気系に反映させ、新し
い解析的な解が求められる。図３に示す電力系統のＦ点
の事故直後において、電力系統電気系モードは次式
（８）により示される。
【００４１】
【数１０】
【００４２】ここで、λは減衰定数、Ｉ1Sは電流定常
解、Ａは初期条件により決定する任意定数である。
解、Ａは初期条件により決定する任意定数である。
【数１１】 は電力系統電気系の過渡減衰直流であり、
【数１２】 は電力系統電気系の定常交流成分である。ある一定の時
間、例えば数１０ミリ秒において電力系統同期機の位相
角を一定と仮定し、（８）式より電力系統電気系モード
を解析的に求めることができる。一方、電力系統機械系
モードは次式（９）により表わされる。
間、例えば数１０ミリ秒において電力系統同期機の位相
角を一定と仮定し、（８）式より電力系統電気系モード
を解析的に求めることができる。一方、電力系統機械系
モードは次式（９）により表わされる。
【００４３】
【数１３】
【００４４】ここで、δは電力系統同期機位相角、Ｐm
は電力系統同期機機械的入力、ｉ1は電力系統同期機の
電流、Ｒ1 は電力系統同期機電機子抵抗、Ｍは電力系統
同期機慣性定数であり、数値積分で（９）式を解くこと
ができる。このようにして電力系統電気系モードと電力
系統機械系モードの分離を実現して解析を可能にしてい
る。
は電力系統同期機機械的入力、ｉ1は電力系統同期機の
電流、Ｒ1 は電力系統同期機電機子抵抗、Ｍは電力系統
同期機慣性定数であり、数値積分で（９）式を解くこと
ができる。このようにして電力系統電気系モードと電力
系統機械系モードの分離を実現して解析を可能にしてい
る。
【００４５】
従って、従来ではディジタル発電機モデル
を模擬する電力系統シミュレータにおいては、電力系統
電気系と電力系統機械系とを区別しないことにより事故
直後に発生する過渡直流電流を再現できず、実際の事故
発生時の現象を正確に再現できないのに対し、本実施の
形態の電力系統シミュレーション方法では事故直後に発
生する過渡直流電流を再現でき、精度の良い電力系統シ
ミュレータの実現に寄与できる効果がある。
を模擬する電力系統シミュレータにおいては、電力系統
電気系と電力系統機械系とを区別しないことにより事故
直後に発生する過渡直流電流を再現できず、実際の事故
発生時の現象を正確に再現できないのに対し、本実施の
形態の電力系統シミュレーション方法では事故直後に発
生する過渡直流電流を再現でき、精度の良い電力系統シ
ミュレータの実現に寄与できる効果がある。
【００４６】
実施の形態３．
図７は、図３に示した電力系統図において全ての電気抵
抗を無視したときの１機無限大母線系統図である。図７
において図３と同一または相当の部分については同一の
符号を付し説明を省略する。この系統図においてＦ点で
地絡事故を発生させ、一定の時間が経過した後、線路遮
断器ＣＢ１と線路遮断器ＣＢ２を動作させ地絡事故の発
生した線路を開放するシミュレーションを行う。同期機
の内部電圧をｅ1 、同期機の端子電圧をｖ1 、同期機の
電流をｉ1 （無限大母線電圧をｖb ）と仮定する。この
結果、次式（１０），（１１）が成立する。ここで全て
の変数がスパイラルベクトルである。
抗を無視したときの１機無限大母線系統図である。図７
において図３と同一または相当の部分については同一の
符号を付し説明を省略する。この系統図においてＦ点で
地絡事故を発生させ、一定の時間が経過した後、線路遮
断器ＣＢ１と線路遮断器ＣＢ２を動作させ地絡事故の発
生した線路を開放するシミュレーションを行う。同期機
の内部電圧をｅ1 、同期機の端子電圧をｖ1 、同期機の
電流をｉ1 （無限大母線電圧をｖb ）と仮定する。この
結果、次式（１０），（１１）が成立する。ここで全て
の変数がスパイラルベクトルである。
【００４７】
ｅ1 ＝ｖ1 ＋Ｌs ｐｉ1 （１０）
ｖ1 ＝ｖb ＋（Ｌ／２）ｐｉ1 （１１）
【００４８】
事故前の定常状態においては、スパイラル
ベクトルの大きさと回転速度は一定であるので、式（１
０），式（１１）は次式（１２），（１３）となる。図
８は事故前の定常状態における定常ベクトル図である。
ベクトルの大きさと回転速度は一定であるので、式（１
０），式（１１）は次式（１２），（１３）となる。図
８は事故前の定常状態における定常ベクトル図である。
【００４９】
【数１４】
【００５０】ここで、Ｅ1 は電力系統同期機内部電圧定
常ベクトル、Ｖ1 電力系統同期機端子電圧定常ベクト
ル、Ｖb は電力系統無限大母線定常ベクトル、Ｉ1 は電
力系統同期機端子電流定常ベクトル、ωは電力系統定常
角速度である。一方、電力系統内部電圧、電流、無限大
母線電圧の瞬時値は次式（１４），（１５），（１６）
により表わすことができる。
常ベクトル、Ｖ1 電力系統同期機端子電圧定常ベクト
ル、Ｖb は電力系統無限大母線定常ベクトル、Ｉ1 は電
力系統同期機端子電流定常ベクトル、ωは電力系統定常
角速度である。一方、電力系統内部電圧、電流、無限大
母線電圧の瞬時値は次式（１４），（１５），（１６）
により表わすことができる。
【００５１】
【数１５】
【００５２】また、電力系統同期機の有効電力Ｐｅおよ
び無効電力Ｑｅは次式（１７），（１８）により求める
ことができる。
び無効電力Ｑｅは次式（１７），（１８）により求める
ことができる。
【００５３】
【数１６】
【００５４】電力系統の過渡状態において、スパイラル
ベクトル理論より式（１０），（１１）の同期機電流一
般解ｉ1sは次式（１９）になる。
ベクトル理論より式（１０），（１１）の同期機電流一
般解ｉ1sは次式（１９）になる。
【００５５】
【数１７】
【００５６】ここで、Ｉ1sは電流定常解、Ａは初期条件
により決定する任意定数である。事故発生直後において
同期機電流定常解は次式（２０）の通りである。
により決定する任意定数である。事故発生直後において
同期機電流定常解は次式（２０）の通りである。
【００５７】
【数１８】
【００５８】この式（２０）と初期条件により式（１
９）の任意定数Ａの値を決定することができる。また、
内部電圧ｅ１の大きさは一定とする。一方、力学系の位
相角は式（７）により求められる。式（１９）から明ら
かなように、電力系統の電気抵抗を無視すれば電力系統
電気系の過渡現象は発生しない。この場合、電力系統の
過渡状態において電力系統機械系過渡現象のみ存在す
る。
９）の任意定数Ａの値を決定することができる。また、
内部電圧ｅ１の大きさは一定とする。一方、力学系の位
相角は式（７）により求められる。式（１９）から明ら
かなように、電力系統の電気抵抗を無視すれば電力系統
電気系の過渡現象は発生しない。この場合、電力系統の
過渡状態において電力系統機械系過渡現象のみ存在す
る。
【００５９】
以上のように、この実施の形態３によれ
ば、電力系統電気系の電気抵抗を無視することで、電力
系統電気系の過渡現象を無視することができ、電力系統
過渡状態において電力系統機械系の過渡現象に対応する
モードを算出することができるため、電力系統の解析や
安定化装置の設定の際の指針を与えることが可能となる
効果がある。
ば、電力系統電気系の電気抵抗を無視することで、電力
系統電気系の過渡現象を無視することができ、電力系統
過渡状態において電力系統機械系の過渡現象に対応する
モードを算出することができるため、電力系統の解析や
安定化装置の設定の際の指針を与えることが可能となる
効果がある。
【００６０】
実施の形態４．
この実施の形態４では、スパイラルベクトル理論を電力
系統シミュレーションに用いることで、電力系統電気系
過渡現象を解析的に解明できることを説明する。図３に
示した電力系統図においてＦ点で地絡事故を発生させ、
一定の時間が経過した後に線路遮断器ＣＢ１と線路遮断
器ＣＢ２を動作させ、線路を開放したときのシミュレー
ションを行ったときの、同期機の内部電圧ｅ1 、同期機
の端子電圧ｖ1 、同期機の電流をｉ1 （無限大母線電圧
をｖb ）と仮定すると、式（１）および式（２）が成立
し、また電力系統の定常状態においては式（３）および
式（４）が成立し、その定常ベクトル図は図４に示され
ることは前記実施の形態１で既に述べたところである。
そして、事故発生後および線路遮断器による線路の開放
で事故をクリアした後の電力系統の過渡状態では、スパ
イラルベクトル理論による電流一般解ｉ1sは式（５）の
通りであることも述べた。式（５）において
系統シミュレーションに用いることで、電力系統電気系
過渡現象を解析的に解明できることを説明する。図３に
示した電力系統図においてＦ点で地絡事故を発生させ、
一定の時間が経過した後に線路遮断器ＣＢ１と線路遮断
器ＣＢ２を動作させ、線路を開放したときのシミュレー
ションを行ったときの、同期機の内部電圧ｅ1 、同期機
の端子電圧ｖ1 、同期機の電流をｉ1 （無限大母線電圧
をｖb ）と仮定すると、式（１）および式（２）が成立
し、また電力系統の定常状態においては式（３）および
式（４）が成立し、その定常ベクトル図は図４に示され
ることは前記実施の形態１で既に述べたところである。
そして、事故発生後および線路遮断器による線路の開放
で事故をクリアした後の電力系統の過渡状態では、スパ
イラルベクトル理論による電流一般解ｉ1sは式（５）の
通りであることも述べた。式（５）において
【数１９】 は電力系統電気系の過渡直流減衰分、
【数２０】 は電力系統電気系の定常交流分である。図９の（ａ）に
過渡直流減衰分、同図（ｂ）に電力系統電気系の定常交
流分、同図（ｃ）にｉ1sの波形をそれぞれ示す。このよ
うに、スパイラルベクトル理論を適用することで得られ
た解により電力系統過渡状態における電力系統過渡現象
を解析的に分析することが可能になる。以上説明したよ
うに、スパイラルベクトル理論を用いることで電力系統
電気系過渡現象を解析することができ、電力系統の高速
安定化制御の実施を可能にする効果がある。
過渡直流減衰分、同図（ｂ）に電力系統電気系の定常交
流分、同図（ｃ）にｉ1sの波形をそれぞれ示す。このよ
うに、スパイラルベクトル理論を適用することで得られ
た解により電力系統過渡状態における電力系統過渡現象
を解析的に分析することが可能になる。以上説明したよ
うに、スパイラルベクトル理論を用いることで電力系統
電気系過渡現象を解析することができ、電力系統の高速
安定化制御の実施を可能にする効果がある。
【００６１】
実施の形態５．
この実施の形態５では、電力系統同期機至近端事故直後
における逆動揺現象が同期機電機子抵抗にあることを説
明する。図３に示す電力系統図のＦ点地絡事故直後の電
力系統機械系の運動方程式は、次式（２１）で表わすこ
とができる。
における逆動揺現象が同期機電機子抵抗にあることを説
明する。図３に示す電力系統図のＦ点地絡事故直後の電
力系統機械系の運動方程式は、次式（２１）で表わすこ
とができる。
【００６２】
【数２１】
【００６３】ここで、δは電力系統同期機位相角、Ｐm
は電力系統同期機機械的入力、ｉ1は電力系統同期機の
電流、Ｒ1 は電力系統同期機電機子抵抗、Ｍは電力系統
同期機慣性定数である。数十ミリ秒程度のスパンで見た
場合、この時間内での電力系統機械系の変数であるＰm
は一定とみなすことができる。一方、電力系統電気系の
変数である
は電力系統同期機機械的入力、ｉ1は電力系統同期機の
電流、Ｒ1 は電力系統同期機電機子抵抗、Ｍは電力系統
同期機慣性定数である。数十ミリ秒程度のスパンで見た
場合、この時間内での電力系統機械系の変数であるＰm
は一定とみなすことができる。一方、電力系統電気系の
変数である
【数２２】 の変動は大きい。その値はＰm より大きくなることも可
能であり、この結果、事故後ある一定時間内この現象を
続けると、電力系統同期機位相角δの値は事故直後にお
いて事故前の値よりも小さくなり、電力系統同期機位相
角δは安定しなくなって動揺する。これを電力系統同期
機逆動揺現象と言う。
能であり、この結果、事故後ある一定時間内この現象を
続けると、電力系統同期機位相角δの値は事故直後にお
いて事故前の値よりも小さくなり、電力系統同期機位相
角δは安定しなくなって動揺する。これを電力系統同期
機逆動揺現象と言う。
【００６４】
一方また、図７に示した電気抵抗を無視し
た電力系統図におけるＦ点地絡事故直後の電力系統機械
系の運動方程式も式（２１）で表わすことができる。し
かしながら、電力系統同期機の電機子抵抗を無視したた
め式（２１）の右辺の第２項は零である。従って、電力
系統同期機の電機子抵抗を無視できれば逆動揺現象は発
生しない。このことから、電力系統同期機至近端事故直
後における逆動揺現象の原因は同期機電機子抵抗にある
ことが明らかになる。
た電力系統図におけるＦ点地絡事故直後の電力系統機械
系の運動方程式も式（２１）で表わすことができる。し
かしながら、電力系統同期機の電機子抵抗を無視したた
め式（２１）の右辺の第２項は零である。従って、電力
系統同期機の電機子抵抗を無視できれば逆動揺現象は発
生しない。このことから、電力系統同期機至近端事故直
後における逆動揺現象の原因は同期機電機子抵抗にある
ことが明らかになる。
【００６５】
図１０は、電力系統同期機至近端事故直後
における位相角の変化を示す特性図である。δ0 は事故
発生前の電力系統同期機位相角、曲線Ａ−Ｂ−Ｃは電力
系統同期機の電機子抵抗の影響で逆動揺現象が発生した
ときの位相角の変化を示す特性曲線、曲線Ａ−Ｄ−Ｅは
電力系統同期機の電機子抵抗を無視したため逆動揺が発
生していないときの位相角の変化を示す特性曲線であ
る。
における位相角の変化を示す特性図である。δ0 は事故
発生前の電力系統同期機位相角、曲線Ａ−Ｂ−Ｃは電力
系統同期機の電機子抵抗の影響で逆動揺現象が発生した
ときの位相角の変化を示す特性曲線、曲線Ａ−Ｄ−Ｅは
電力系統同期機の電機子抵抗を無視したため逆動揺が発
生していないときの位相角の変化を示す特性曲線であ
る。
【００６６】
以上説明したように、この実施の形態５で
は、電力系統同期機至近端事故直後における逆動揺現象
の原因は同期機の電機子抵抗であることが判明し、これ
を電力系統同期機の設計や試験に際しての指針として用
いることが可能になる効果がある。
は、電力系統同期機至近端事故直後における逆動揺現象
の原因は同期機の電機子抵抗であることが判明し、これ
を電力系統同期機の設計や試験に際しての指針として用
いることが可能になる効果がある。
【００６７】
実施の形態６．
この実施の形態６では、電力系統電機系をスパイラルベ
クトル理論を適用して解析し、電力系統機械系を数値積
分により解析する解析手法により、並列演算処理を行う
ことが可能になって、リアルタイムな電力系統シミュレ
ーションや電力系統シミュレータの実現に寄与すること
を説明する。
クトル理論を適用して解析し、電力系統機械系を数値積
分により解析する解析手法により、並列演算処理を行う
ことが可能になって、リアルタイムな電力系統シミュレ
ーションや電力系統シミュレータの実現に寄与すること
を説明する。
【００６８】
図１１は、同期発電機Ａおよび同期発電機
Ｂを有した電力系統における電力系統電気系および電力
系統機械系の各演算部分に並列演算ノードを設定したと
きの説明図である。図において、１００は電力系統電気
系演算部分に設定された並列演算ノード、２００は電力
系統同期発電機Ａの機械系演算部分に設定された並列演
算ノード、３００は電力系統同期発電機Ｂの機械系演算
部分に設定された並列演算ノード、４００は電力系統誘
導機の機械系演算部分に設定された並列演算ノードを示
す。このような各並列演算ノードを設定して電力系統電
気系をスパイラルベクトル理論により解析し、電力系統
機械系は数値積分により解析する。この際に並列演算技
術を用いることでリアルタイムな電力系統シミュレーシ
ョン、電力系統シミュレータを実現できる。
Ｂを有した電力系統における電力系統電気系および電力
系統機械系の各演算部分に並列演算ノードを設定したと
きの説明図である。図において、１００は電力系統電気
系演算部分に設定された並列演算ノード、２００は電力
系統同期発電機Ａの機械系演算部分に設定された並列演
算ノード、３００は電力系統同期発電機Ｂの機械系演算
部分に設定された並列演算ノード、４００は電力系統誘
導機の機械系演算部分に設定された並列演算ノードを示
す。このような各並列演算ノードを設定して電力系統電
気系をスパイラルベクトル理論により解析し、電力系統
機械系は数値積分により解析する。この際に並列演算技
術を用いることでリアルタイムな電力系統シミュレーシ
ョン、電力系統シミュレータを実現できる。
【００６９】
【発明の効果】以上のように、請求項１記載の発明によ
れば、電力系統の電気系過渡現象の解析にスパイラルベ
クトル理論を適用し、
れば、電力系統の電気系過渡現象の解析にスパイラルベ
クトル理論を適用し、
【数２３】 により、前記電力系統の電流一般解の解析を電力系統電
気系解析過程により行い、前記電力系統電気系解析過程
による解析結果をもとに、
気系解析過程により行い、前記電力系統電気系解析過程
による解析結果をもとに、
【数２４】 により、前記電力系統の機械系の解析を電力系統機械系
解析過程で行い、この電力系統機械系解析過程による解
析結果を前記電力系統電気系解析過程で得られる解析結
果の各定数を修正して前記電力系統機械系解析過程によ
る解析結果を前記電力系統電気系解析過程へ反映過程で
反映させるように構成したので、スパイラルベクトル理
論の適用により現実に忠実な電力系統電気系の解析を行
うことができ、さらにこの現実に即した電力系統電気系
の解析結果により前記電気系の解析と同様に前記電力系
統の機械系の解析を現実に忠実に行うことができる効果
がある。また、反映過程で前記電力系統の機械系の解析
結果を電力系統電気系の解析に反映させるため前記電気
系と前記機械系とが互いに関連し合って電力系統で発生
する過渡現象を現実に忠実に模擬できる効果がある。ま
た、さらに電力系統電気系と電力系統機械系とを別々に
解析するため、電力系統シミュレーションを行う際の計
算機利用に対しての適合性が良い効果があり、さらに、
時定数の異なる前記電力系統の機械系と電気系の過渡現
象の解析およびそのシミュレーションを現実に忠実に行
うことができる効果がある。
解析過程で行い、この電力系統機械系解析過程による解
析結果を前記電力系統電気系解析過程で得られる解析結
果の各定数を修正して前記電力系統機械系解析過程によ
る解析結果を前記電力系統電気系解析過程へ反映過程で
反映させるように構成したので、スパイラルベクトル理
論の適用により現実に忠実な電力系統電気系の解析を行
うことができ、さらにこの現実に即した電力系統電気系
の解析結果により前記電気系の解析と同様に前記電力系
統の機械系の解析を現実に忠実に行うことができる効果
がある。また、反映過程で前記電力系統の機械系の解析
結果を電力系統電気系の解析に反映させるため前記電気
系と前記機械系とが互いに関連し合って電力系統で発生
する過渡現象を現実に忠実に模擬できる効果がある。ま
た、さらに電力系統電気系と電力系統機械系とを別々に
解析するため、電力系統シミュレーションを行う際の計
算機利用に対しての適合性が良い効果があり、さらに、
時定数の異なる前記電力系統の機械系と電気系の過渡現
象の解析およびそのシミュレーションを現実に忠実に行
うことができる効果がある。
【００７０】
請求項２記載の発明によれば、シミュレー
トしようとする状況についての電力系統の初期状態に応
じた定数を設定する初期条件設定過程を有し、該初期条
件設定過程で前記定数が設定された前記初期状態に応じ
た前記状況発生時の過渡現象について前記電力系統の電
気系の解析を行う電力系統電気系解析過程を備えるよう
に構成したので、前記初期状態に応じた前記状況につい
ての電力系統電気系の解析を行うことができる効果があ
る。
トしようとする状況についての電力系統の初期状態に応
じた定数を設定する初期条件設定過程を有し、該初期条
件設定過程で前記定数が設定された前記初期状態に応じ
た前記状況発生時の過渡現象について前記電力系統の電
気系の解析を行う電力系統電気系解析過程を備えるよう
に構成したので、前記初期状態に応じた前記状況につい
ての電力系統電気系の解析を行うことができる効果があ
る。
【００７１】
請求項３記載の発明によれば、電力系統で
発生する負荷変動の発生条件初期状態に応じた定数を設
定する初期条件設定過程を有し、該初期条件設定過程で
前記定数が設定された負荷変動発生の際の条件である初
期状態に応じた負荷変動発生時の過渡現象について前記
電力系統の電気系の解析を行う電力系統電気系解析過程
と、該電力系統電気系解析過程による前記負荷変動発生
時の解析結果をもとに前記負荷変動発生時の前記電力系
統の機械系の解析を行う電力系統機械系解析過程とを備
えるように構成したので、前記発生条件初期状態に応じ
た負荷変動時の状況についての電力系統電気系の解析を
行うことができる効果がある。
発生する負荷変動の発生条件初期状態に応じた定数を設
定する初期条件設定過程を有し、該初期条件設定過程で
前記定数が設定された負荷変動発生の際の条件である初
期状態に応じた負荷変動発生時の過渡現象について前記
電力系統の電気系の解析を行う電力系統電気系解析過程
と、該電力系統電気系解析過程による前記負荷変動発生
時の解析結果をもとに前記負荷変動発生時の前記電力系
統の機械系の解析を行う電力系統機械系解析過程とを備
えるように構成したので、前記発生条件初期状態に応じ
た負荷変動時の状況についての電力系統電気系の解析を
行うことができる効果がある。
【００７２】
請求項４記載の発明によれば、負荷変動発
生以降の過渡現象についての電力系統の電気系の解析結
果から得られた前記電力系統同期機の発生出力を求める
電力系統同期機発生出力演算過程と、該電力系統同期機
発生出力演算過程で求めた前記電力系統同期機の発生出
力をもとに電力系統機械系の状態変化量を求める状態変
化量演算過程と、該状態変化量演算過程で求めた状態変
化量に応じて電力系統電気系解析過程で得られた解析結
果の各定数を修正することで、電力系統機械系解析過程
による解析結果を前記電力系統電気系解析過程へ反映さ
せる反映過程とを備えるように構成したので、前記電力
系統同期機の発生出力と前記電力系統機械系の状態変化
量とをもとに、負荷変動発生以降の過渡現象について互
いに関連し合う電力系統電気系および電力系統機械系の
解析を行うことができる効果がある。
生以降の過渡現象についての電力系統の電気系の解析結
果から得られた前記電力系統同期機の発生出力を求める
電力系統同期機発生出力演算過程と、該電力系統同期機
発生出力演算過程で求めた前記電力系統同期機の発生出
力をもとに電力系統機械系の状態変化量を求める状態変
化量演算過程と、該状態変化量演算過程で求めた状態変
化量に応じて電力系統電気系解析過程で得られた解析結
果の各定数を修正することで、電力系統機械系解析過程
による解析結果を前記電力系統電気系解析過程へ反映さ
せる反映過程とを備えるように構成したので、前記電力
系統同期機の発生出力と前記電力系統機械系の状態変化
量とをもとに、負荷変動発生以降の過渡現象について互
いに関連し合う電力系統電気系および電力系統機械系の
解析を行うことができる効果がある。
【００７３】
請求項５記載の発明によれば、電力系統に
おける事故の発生する条件である初期状態に応じた定数
を設定する初期条件設定過程と、該初期条件設定過程で
前記定数が設定された初期状態に応じた事故発生時の過
渡現象について前記電力系統の電気系の解析を行う電力
系統電気系解析過程と、該電力系統電気系解析過程によ
る前記事故発生時の解析結果をもとに前記事故発生時の
前記電力系統の機械系の解析を行う電力系統機械系解析
過程とを備えるように構成したので、前記初期状態に応
じた前記電力系統の事故についての電力系統電気系の解
析を行うことができる効果がある。
おける事故の発生する条件である初期状態に応じた定数
を設定する初期条件設定過程と、該初期条件設定過程で
前記定数が設定された初期状態に応じた事故発生時の過
渡現象について前記電力系統の電気系の解析を行う電力
系統電気系解析過程と、該電力系統電気系解析過程によ
る前記事故発生時の解析結果をもとに前記事故発生時の
前記電力系統の機械系の解析を行う電力系統機械系解析
過程とを備えるように構成したので、前記初期状態に応
じた前記電力系統の事故についての電力系統電気系の解
析を行うことができる効果がある。
【００７４】
請求項６記載の発明によれば、事故発生以
降の過渡現象についての電力系統の電気系の解析結果か
ら電力系統同期機の発生出力を求める電力系統同期機発
生出力演算過程と、該電力系統同期機発生出力演算過程
で求めた前記電力系統同期機の発生出力をもとに電力系
統機械系の状態変化量を演算し求める状態変化量演算過
程と、該状態変化量演算過程で求めた状態変化量に応じ
て電力系統電気系解析過程で得られた解析結果の各定数
を修正することで、電力系統機械系解析過程による解析
結果を前記電力系統電気系解析過程へ反映させる反映過
程とを備えるように構成したので、初期状態に応じた前
記電力系統の事故について、前記電力系統同期機の発生
出力と前記電力系統機械系の状態変化量とをもとに、負
荷変動発生以降の過渡現象について互いに関連し合う電
力系統電気系および電力系統機械系の解析を行うことが
できる効果がある。
降の過渡現象についての電力系統の電気系の解析結果か
ら電力系統同期機の発生出力を求める電力系統同期機発
生出力演算過程と、該電力系統同期機発生出力演算過程
で求めた前記電力系統同期機の発生出力をもとに電力系
統機械系の状態変化量を演算し求める状態変化量演算過
程と、該状態変化量演算過程で求めた状態変化量に応じ
て電力系統電気系解析過程で得られた解析結果の各定数
を修正することで、電力系統機械系解析過程による解析
結果を前記電力系統電気系解析過程へ反映させる反映過
程とを備えるように構成したので、初期状態に応じた前
記電力系統の事故について、前記電力系統同期機の発生
出力と前記電力系統機械系の状態変化量とをもとに、負
荷変動発生以降の過渡現象について互いに関連し合う電
力系統電気系および電力系統機械系の解析を行うことが
できる効果がある。
【００７５】
請求項７記載の発明によれば、状態変化量
演算過程で求めた電力系統機械系の状態変化量を所定の
閾値と比較する比較過程と、前記状態変化量演算過程で
求めた状態変化量が前記所定の閾値を越えているとの比
較結果が前記比較過程で得られると電力系統電気系解析
過程で得られた解析結果の各定数を修正し、前記電力系
統機械系解析過程による解析結果を前記電力系統電気系
解析過程へ反映させる反映過程を備えるように構成した
ので、電力系統機械系の状態変化量が電力系統シミュレ
ーション結果に影響する量変化した場合に電力系統シミ
ュレーションの次のステップへ進むことになり、精度の
よい電力系統シミュレーションを行うことができる効果
がある。
演算過程で求めた電力系統機械系の状態変化量を所定の
閾値と比較する比較過程と、前記状態変化量演算過程で
求めた状態変化量が前記所定の閾値を越えているとの比
較結果が前記比較過程で得られると電力系統電気系解析
過程で得られた解析結果の各定数を修正し、前記電力系
統機械系解析過程による解析結果を前記電力系統電気系
解析過程へ反映させる反映過程を備えるように構成した
ので、電力系統機械系の状態変化量が電力系統シミュレ
ーション結果に影響する量変化した場合に電力系統シミ
ュレーションの次のステップへ進むことになり、精度の
よい電力系統シミュレーションを行うことができる効果
がある。
【００７６】
請求項８記載の発明によれば、状態変化量
演算過程で求めた状態変化量が所定の閾値を越えていな
いとの比較結果が前記比較過程で得られると、電力系統
同期機発生出力演算過程での電力系統同期機の発生出力
の演算を行い、この電力系統同期機の発生出力をもとに
前記状態変化量演算過程で求めた状態変化量と前記所定
の閾値との比較を前記比較過程で行う再比較過程を備え
るように構成したので、電力系統機械系の状態変化量が
電力系統シミュレーション結果に影響する量変化してい
ない場合には電力系統シミュレーションの次のステップ
へ進むことなく、前記状態変化量が電力系統シミュレー
ション結果に影響する量変化するまで待って電力系統シ
ミュレーションの次のステップへ進むことになり、精度
のよい電力系統シミュレーションを行うことができる効
果がある。
演算過程で求めた状態変化量が所定の閾値を越えていな
いとの比較結果が前記比較過程で得られると、電力系統
同期機発生出力演算過程での電力系統同期機の発生出力
の演算を行い、この電力系統同期機の発生出力をもとに
前記状態変化量演算過程で求めた状態変化量と前記所定
の閾値との比較を前記比較過程で行う再比較過程を備え
るように構成したので、電力系統機械系の状態変化量が
電力系統シミュレーション結果に影響する量変化してい
ない場合には電力系統シミュレーションの次のステップ
へ進むことなく、前記状態変化量が電力系統シミュレー
ション結果に影響する量変化するまで待って電力系統シ
ミュレーションの次のステップへ進むことになり、精度
のよい電力系統シミュレーションを行うことができる効
果がある。
【００７７】
請求項９記載の発明によれば、電力系統同
期機発生出力演算過程で求めた電力系統同期機の発生出
力をもとに電力系統機械系の位相角である状態変化量を
状態変化量演算過程が求めるように構成したので、初期
状態に応じた電力系統の負荷変動や事故について、前記
電力系統同期機の発生出力と前記電力系統機械系の状態
変化量とをもとに、負荷変動発生以降や事故発生以降の
過渡現象について、互いに関連し合う電力系統電気系お
よび電力系統機械系の解析を行うことができる効果があ
る。
期機発生出力演算過程で求めた電力系統同期機の発生出
力をもとに電力系統機械系の位相角である状態変化量を
状態変化量演算過程が求めるように構成したので、初期
状態に応じた電力系統の負荷変動や事故について、前記
電力系統同期機の発生出力と前記電力系統機械系の状態
変化量とをもとに、負荷変動発生以降や事故発生以降の
過渡現象について、互いに関連し合う電力系統電気系お
よび電力系統機械系の解析を行うことができる効果があ
る。
【００７８】
請求項１０記載の発明によれば、事故発生
以降および線路開放以降の系統状態の急変の有無につい
て判定する系統状態急変判定過程を備えるように構成し
たので、前記系統状態の急変の有無に応じた電気系と機
械系とが互いに関連し合って電力系統で発生する過渡現
象を現実に忠実に模擬できる効果がある。
以降および線路開放以降の系統状態の急変の有無につい
て判定する系統状態急変判定過程を備えるように構成し
たので、前記系統状態の急変の有無に応じた電気系と機
械系とが互いに関連し合って電力系統で発生する過渡現
象を現実に忠実に模擬できる効果がある。
【００７９】
請求項１１記載の発明によれば、系統状態
急変判定過程により系統状態の急変があったことが判定
されたときの電力系統電気系の過渡現象について解析を
行う系統状態急変時電気系解析過程と、該系統状態急変
時電気系解析過程の解析結果をもとに前記電力系統が急
変したときの前記電力系統の機械系の解析を行う電力系
統機械系解析過程とを備えるように構成したので、前記
系統状態の急変があったときの前記電気系と前記機械系
とが互いに関連し合って電力系統で発生する過渡現象を
現実に忠実に模擬できる効果がある。
急変判定過程により系統状態の急変があったことが判定
されたときの電力系統電気系の過渡現象について解析を
行う系統状態急変時電気系解析過程と、該系統状態急変
時電気系解析過程の解析結果をもとに前記電力系統が急
変したときの前記電力系統の機械系の解析を行う電力系
統機械系解析過程とを備えるように構成したので、前記
系統状態の急変があったときの前記電気系と前記機械系
とが互いに関連し合って電力系統で発生する過渡現象を
現実に忠実に模擬できる効果がある。
【００８０】
請求項１２記載の発明によれば、初期条件
設定過程で定数が設定された初期状態に応じた事故発生
時の過渡現象について過渡直流電流を含む電力系統の電
気系の解析を行う電力系統電気系解析過程を備えるよう
に構成したので、現実の現象に忠実な過渡直流電流を含
む電力系統シミュレーションを実現できる効果がある。
設定過程で定数が設定された初期状態に応じた事故発生
時の過渡現象について過渡直流電流を含む電力系統の電
気系の解析を行う電力系統電気系解析過程を備えるよう
に構成したので、現実の現象に忠実な過渡直流電流を含
む電力系統シミュレーションを実現できる効果がある。
【００８１】
請求項１３記載の発明によれば、複数の電
力系統同期機を有した電力系統の電力系統電気系解析過
程および電力系統機械系解析過程における演算にそれぞ
れ並列演算ノードを設定し、前記電力系統電気系解析過
程による電力系統電気系の解析および前記電力系統機械
系解析過程による電力系統機械系の解析を並列的に行な
うように構成したので、計算機利用、特にマルチタスク
処理を行うことのできる計算機を有効利用して電力系統
シミュレーションを行うことができる効果がある。
力系統同期機を有した電力系統の電力系統電気系解析過
程および電力系統機械系解析過程における演算にそれぞ
れ並列演算ノードを設定し、前記電力系統電気系解析過
程による電力系統電気系の解析および前記電力系統機械
系解析過程による電力系統機械系の解析を並列的に行な
うように構成したので、計算機利用、特にマルチタスク
処理を行うことのできる計算機を有効利用して電力系統
シミュレーションを行うことができる効果がある。
【図１】 この発明の実施の形態１における電力系統事
故発生時の電力系統シミュレーション方法の構成を示す
フローチャートである。
故発生時の電力系統シミュレーション方法の構成を示す
フローチャートである。
【図２】 この発明の実施の形態１による電力系統シミ
ュレーション方法におけるスパイラルベクトルを示す説
明図である。
ュレーション方法におけるスパイラルベクトルを示す説
明図である。
【図３】 この発明の実施の形態１における電力系統シ
ミュレーション方法が適用される電力系統図である。
ミュレーション方法が適用される電力系統図である。
【図４】 この発明の実施の形態１による電力系統シミ
ュレーション方法が適用される図３に示す電力系統の定
常状態における電力系統同期機内部電圧および電力系統
同期機端子電圧を示すベクトル図である。
ュレーション方法が適用される図３に示す電力系統の定
常状態における電力系統同期機内部電圧および電力系統
同期機端子電圧を示すベクトル図である。
【図５】 この発明の実施の形態２による電力系統過渡
状態における動揺モードを示す説明図である。
状態における動揺モードを示す説明図である。
【図６】 この発明の実施の形態２による電力系統過渡
状態における電力系統電気系動揺モードと電力系統機械
系動揺モードとの比較を示す説明図である。
状態における電力系統電気系動揺モードと電力系統機械
系動揺モードとの比較を示す説明図である。
【図７】 この発明の実施の形態３による電力系統シミ
ュレーション方法における電力系統の全ての電気抵抗を
無視したときの１機無限大母線系統図である。
ュレーション方法における電力系統の全ての電気抵抗を
無視したときの１機無限大母線系統図である。
【図８】 この発明の実施の形態３による事故前の定常
状態における電力系統同期機内部電圧および電力系統同
期機端子電圧を示すベクトル図である。
状態における電力系統同期機内部電圧および電力系統同
期機端子電圧を示すベクトル図である。
【図９】 この発明の実施の形態４における価値直流減
衰分、定常交流分および過渡電流を示す波形図である。
衰分、定常交流分および過渡電流を示す波形図である。
【図１０】 この発明の実施の形態５における電力系統
同期機至近端事故直後の電力系統同期機の位相角の変化
を示す特性図である。
同期機至近端事故直後の電力系統同期機の位相角の変化
を示す特性図である。
【図１１】 この発明の実施の形態６による同期発電機
Ａおよび同期発電機Ｂを有した電力系統における電力系
統電気系および電力系統機械系の各演算部分に並列演算
ノードを設定したときの説明図である。
Ａおよび同期発電機Ｂを有した電力系統における電力系
統電気系および電力系統機械系の各演算部分に並列演算
ノードを設定したときの説明図である。
【図１２】 従来の電力系統シミュレーション方法の構
成を示すフローチャートである。
成を示すフローチャートである。
ステップＳＴ１２ 初期条件設定過程、ステップＳＴ１
３ 電力系統電気系解析過程、ステップＳＴ１４ 電力
系統機械系解析過程，電力系統同期機発生出力演算過
程，状態変化量演算過程）、ステップＳＴ１５ 系統状
態急変判定過程、ステップＳＴ１７ 比較過程，再比較
過程、ステップＳＴ１８ 系統状態急変時電気系解析過
程、ステップＳＴ１９ 反映過程。
３ 電力系統電気系解析過程、ステップＳＴ１４ 電力
系統機械系解析過程，電力系統同期機発生出力演算過
程，状態変化量演算過程）、ステップＳＴ１５ 系統状
態急変判定過程、ステップＳＴ１７ 比較過程，再比較
過程、ステップＳＴ１８ 系統状態急変時電気系解析過
程、ステップＳＴ１９ 反映過程。
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平３−256528（ＪＰ，Ａ)
迎久雄、他、スパイラルベクトル理論
による同期機モデルの三相短絡事故解
析，電気学会論文誌Ｂ 電力・エネルギ
ー部門誌，日本，社団法人電気学会、
1995年 ５月20日，Ｖｏｌ．115−Ｂ,
ｐ．624−630
番場隆治、他，スパイラルベクトル理
論による同期機モデルを用いた系統解析
（その１），工学院大学研究報告，日
本，工学院大学，1993年 ４月30日，第
74号，ｐ．125−132
番場隆治、他，スパイラルベクトル理
論を用いた系統解析（その２），工学院
大学研究報告，日本，工学院大学，1995
年 ４月30日，第78号，ｐ．135−142
(58)調査した分野(Int.Cl.7，ＤＢ名)
H02J 3/00 - 5/00
ＪＩＣＳＴファイル（ＪＯＩＳ)
Claims (13)
- 【請求項１】 電力系統の電気系過渡現象の解析にスパ
イラルベクトル理論を適用し、【数１】 により、前記電力系統の電流一般解の解析を行う電力系
統電気系解析過程と、 該電力系統電気系解析過程による解析結果をもとに、 【数２】 により、前記電力系統の機械系の解析を行う電力系統機
械系解析過程と、 該電力系統機械系解析過程による解析結果をもとに、前
記電力系統電気系解析過程で得られた解析結果の各定数
を修正して、前記電力系統機械系解析過程による解析結
果を前記電力系統電気系解析過程へ反映させる反映過程
とを備えた電力系統シミュレーション方法。 - 【請求項２】 電力系統電気系解析過程では、 シミュレートしようとする状況についての電力系統の初
期状態に応じた定数を設定する初期条件設定過程を有
し、該初期条件設定過程で前記定数が設定された前記初
期状態に応じた状況発生時の過渡現象について前記電力
系統の電気系の解析を行い、 電力系統機械系解析過程では、前記電力系統電気系解析
過程による前記状況発生時の解析結果をもとに前記状況
発生時の前記電力系統の機械系の解析を行うことを特徴
とする請求項１記載の電力系統シミュレーション方法。 - 【請求項３】 電力系統電気系解析過程では、 電力系統で発生する負荷変動の発生条件初期状態に応じ
た定数を設定する初期条件設定過程を有し、該初期条件
設定過程で前記定数が設定された負荷変動発生の際の条
件である初期状態に応じた負荷変動発生時の過渡現象に
ついて前記電力系統の電気系の解析を行い、 電力系統機械系解析過程では、前記電力系統電気系解析
過程による前記負荷変動発生時の解析結果をもとに前記
負荷変動発生時の前記電力系統の機械系の解析を行うこ
とを特徴とする請求項１記載の電力系統シミュレーショ
ン方法。 - 【請求項４】 負荷変動発生以降の過渡現象についての
電力系統の電気系の解析結果から得られた電力系統同期
機の発生出力を求める電力系統同期機発生出力演算過程
と、 前記電力系統同期機発生出力演算過程で求めた前記電力
系統同期機の発生出力をもとに電力系統機械系の状態変
化量を求める状態変化量演算過程を有し、 反映過程では、前記状態変化量演算過程で求めた状態変
化量に応じて電力系統電気系解析過程で得られた解析結
果の各定数を修正することで、前記電力系統機械系解析
過程による解析結果を前記電力系統電気系解析過程へ反
映させることを特徴とする請求項３記載の電力系統シミ
ュレーション方法。 - 【請求項５】 電力系統における事故の発生する条件で
ある初期状態に応じた定数を設定する初期条件設定過程
を有し、電力系統電気系解析過程では、前記初期条件設
定過程で前記定数が設定された初期状態に応じた事故発
生時の過渡現象について前記電力系統の電気系の解析を
行い、 電力系統機械系解析過程では、前記電力系統電気系解析
過程による前記事故発生時の解析結果をもとに前記事故
発生時の前記電力系統の機械系の解析を行うことを特徴
とする請求項２記載の電力系統シミュレーション方法。 - 【請求項６】 事故発生以降の過渡現象についての電力
系統の電気系の解析結果から電力系統同期機の発生出力
を求める電力系統同期機発生出力演算過程を有し、 電力系統機械系解析過程は、前記電力系統同期機発生出
力演算過程で求めた前記電力系統同期機の発生出力をも
とに電力系統機械系の状態変化量を演算し求める状態変
化量演算過程を有し、 反映過程では、前記状態変化量演算過程で求めた状態変
化量に応じて電力系統電気系解析過程で得られた解析結
果の各定数を修正することで、前記電力系統機械系解析
過程による解析結果を前記電力系統電気系解析過程へ反
映させることを特徴とする請求項５記載の電力系統シミ
ュレーション方法。 - 【請求項７】 状態変化量演算過程で求めた電力系統機
械系の状態変化量を所定の閾値と比較する比較過程を備
え、反映過程では、前記状態変化量演算過程で求めた状
態変化量が前記所定の閾値を越えているとの比較結果が
前記比較過程で得られると電力系統電気系解析過程で得
られた解析結果の各定数を修正し、前記電力系統機械系
解析過程による解析結果を前記電力系統電気系解析過程
へ反映させることを特徴とする請求項４または請求項６
記載の電力系統シミュレーション方法。 - 【請求項８】 状態変化量演算過程で求めた状態変化量
が所定の閾値を越えていないとの比較結果が前記比較過
程で得られると、電力系統同期機発生出力演算過程での
電力系統同期機の発生出力の演算を行い、この電力系統
同期機の発生出力をもとに前記状態変化量演算過程で求
めた状態変化量と前記所定の閾値との比較を前記比較過
程で行う再比較過程を備えていることを特徴とする請求
項７記載の電力系統シミュレーション方法。 - 【請求項９】 電力系統機械系解析過程の状態変化量演
算過程は、電力系統同期機発生出力演算過程で求めた前
記電力系統同期機の発生出力をもとに電力系統機械系の
位相角である状態変化量を求めることを特徴とする請求
項４または請求項６から請求項８のうちのいずれか１項
記載の電力系統シミュレーション方法。 - 【請求項１０】 事故発生以降および線路開放以降の系
統状態の急変の有無について判定する系統状態急変判定
過程を備えていることを特徴とする請求項１から請求項
９のうちのいずれか１項記載の電力系統シミュレーショ
ン方法。 - 【請求項１１】 系統状態急変判定過程により系統状態
の急変があったことが判定されたときの電力系統電気系
の過渡現象について解析を行う系統状態急変時電気系解
析過程を有し、 電力系統機械系解析過程は、前記系統状態急変時電気系
解析過程の解析結果をもとに前記電力系統が急変したと
きの前記電力系統の機械系の解析を行うことを特徴とす
る請求項１０記載の電力系統シミュレーション方法。 - 【請求項１２】 電力系統電気系解析過程では、初期条
件設定過程で定数が設定された初期状態に応じた事故発
生時の過渡現象について過渡直流電流を含む電力系統の
電気系の解析を行うことを特徴とする請求項５から請求
項１１のうちのいずれか１項記載の電力系統シミュレー
ション方法。 - 【請求項１３】 複数の電力系統同期機を有した電力系
統の電力系統電気系解析過程および電力系統機械系解析
過程における演算にそれぞれ並列演算ノードを設定し、
前記電力系統電気系解析過程による電力系統電気系の解
析および前記電力系統機械系解析過程による電力系統機
械系の解析を並列的に行う請求項４、請求項６から請求
項１２のうちのいずれか１項記載の電力系統シミュレー
ション方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33734496A JP3459158B2 (ja) | 1996-12-17 | 1996-12-17 | 電力系統シミュレーション方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33734496A JP3459158B2 (ja) | 1996-12-17 | 1996-12-17 | 電力系統シミュレーション方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10174290A JPH10174290A (ja) | 1998-06-26 |
| JP3459158B2 true JP3459158B2 (ja) | 2003-10-20 |
Family
ID=18307752
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP33734496A Expired - Fee Related JP3459158B2 (ja) | 1996-12-17 | 1996-12-17 | 電力系統シミュレーション方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3459158B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1321490C (zh) * | 2003-09-28 | 2007-06-13 | 中国电力科学研究院 | 电力***数字仿真方法 |
| JP5078763B2 (ja) * | 2008-06-10 | 2012-11-21 | 三菱電機株式会社 | 送電線事故点標定装置および送電線事故点標定方法 |
| CN110556853B (zh) * | 2018-06-01 | 2023-02-03 | 中国电力科学研究院有限公司 | 一种为电磁暂态仿真提供初值的计算方法及*** |
-
1996
- 1996-12-17 JP JP33734496A patent/JP3459158B2/ja not_active Expired - Fee Related
Non-Patent Citations (3)
| Title |
|---|
| 番場隆治、他，スパイラルベクトル理論による同期機モデルを用いた系統解析（その１），工学院大学研究報告，日本，工学院大学，１９９３年 ４月３０日，第７４号，ｐ．１２５−１３２ |
| 番場隆治、他，スパイラルベクトル理論を用いた系統解析（その２），工学院大学研究報告，日本，工学院大学，１９９５年 ４月３０日，第７８号，ｐ．１３５−１４２ |
| 迎久雄、他、スパイラルベクトル理論による同期機モデルの三相短絡事故解析，電気学会論文誌Ｂ 電力・エネルギー部門誌，日本，社団法人電気学会、１９９５年 ５月２０日，Ｖｏｌ．１１５−Ｂ，ｐ．６２４−６３０ |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH10174290A (ja) | 1998-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Methods of interfacing rotating machine models in transient simulation programs | |
| Amrhein et al. | 3-D magnetic equivalent circuit framework for modeling electromechanical devices | |
| Fouad et al. | Power system transient stability analysis using the transient energy function method | |
| Tse et al. | Refinement of conventional PSS design in multimachine system by modal analysis | |
| Ojo et al. | An improved model of saturated induction machines | |
| Karaagac et al. | An efficient synchronous machine model for electromagnetic transients | |
| WO2018102720A1 (en) | System and method for a fast power network simulator | |
| Kundur et al. | Implementation of advanced generator models into power system stability programs | |
| JP3459158B2 (ja) | 電力系統シミュレーション方法 | |
| Lei et al. | Efficient Power System Transient Simulation Based on Frequency Response Optimized Integrators Considering Second Order Derivative | |
| JP3428885B2 (ja) | 同期機シミュレータ及び同期機シミュレーション方法 | |
| Ramanujam | Power system dynamics: analysis and simulation | |
| Kyesswa et al. | A Matlab-Based Simulation Tool For The Analysis Of Unsymmetrical Power System Transients In Large Networks. | |
| Kulasza et al. | Extending the frequency bandwidth of transient stability simulation using dynamic phasors | |
| Zhou et al. | Steady state analysis of synchronous generators by a coupled field-circuit method | |
| Sedaghati et al. | Transient stability study of power systems with high-order models based on hybridizing loop solving and vector computation | |
| El–Serafi et al. | Saturation representation in synchronous machine models | |
| Barber et al. | Simulation of a synchronous generator connected via a delta—star transformer | |
| Li et al. | Modeling and simulation of aircraft power supply system based on Dymola and Modelica | |
| Richards et al. | Limiting induction motor transient shaft torques following source discontinuities | |
| JPH07135738A (ja) | 電力系統の安定度判定方式 | |
| Wannous et al. | The development of the impedance measured by distance relay during near-to-generator short circuit | |
| Kohlrusz et al. | OrCAD PSpice implementation of a realistic three-phase PMSM model for diagnostic purposes | |
| Ntombela et al. | An investigation into the capabilities of MATLAB power system toolbox for small signal stability analysis in power systems | |
| Lange et al. | Circuit coupled simulation of a claw-pole alternator by a temporary linearization of the 3d-fe model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070808 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080808 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080808 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090808 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090808 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100808 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110808 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110808 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120808 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120808 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130808 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |