(57)【要約】【課題】複雑な振動を伴う動揺波形や複数の動揺モードが干渉する中間領域の波形において極端な誤判別の生じない安定判別方法及びこれらの安定判別方法を組み合わせる事で、誤不動作はせず、誤動作をできる限り抑制する安定判別方法の提供。【解決手段】動揺波形前半部及び後半部における波形の振幅を算出し、その比較結果を用いて判別を行うことを特徴とする自動安定判別方法、動揺波形前半部及び後半部における波形の実効値を算出し、その比較結果を用いて判別を行うことを特徴とする自動安定判別方法及び時系列解析による動揺モード抽出結果を用いて判別を行うことを特徴とする自動安定判別方法のうち、任意の方法の判別結果の論理積を用いて安定判別を行う。【選択図】      図１
2004282882.tif 000002
【特許請求の範囲】
【請求項１】
複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺を一定の時間領域でシミュレーションした結果について、この電力系統の安定度をコンピュータにより判別する方法であって、
入力手段が前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力するステップと、
演算手段が前記時系列の前半部及び後半部それぞれにおける前記入力データの振幅を算出するステップと、
判別手段が前記前半部と後半部との振幅の比較結果を基準データとして電力系統の安定度の判別を行うステップから構成されることを特徴とする電力系統の安定判別方法。
【請求項２】
複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺を一定の時間領域でシミュレーションした結果について、この電力系統の安定度をコンピュータにより判別する方法であって、
入力手段が前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力するステップと、
演算手段が前記時系列の前半部及び後半部それぞれにおける入力データの波形の実効値を算出するステップと、
判別手段が前記前半部と後半部との前記実効値の比較結果を基準データとして電力系統の安定度の判別を行うステップから構成されることを特徴とする電力系統の安定判別方法。
【請求項３】
複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺を一定の時間領域でシミュレーションした結果について、この電力系統の安定度をコンピュータにより判別する方法であって、
入力手段が前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角を入力データとして時系列で入力するステップと、
演算手段が前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルの最適な次数とパラメータを設定すると共に、前記自己回帰モデルを用いて動揺モードの検出を行うステップと、
判別手段が前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルによる動揺モード検出結果を基準データとして安定度の判別を行うステップから構成されることを特徴とする電力系統の安定判別方法。
【請求項４】
前記演算手段が、前記動揺モード検出の際に、前記自己回帰モデルの同定に用いる前記入力データの時間領域の両端での値がその入力データの平均値に近づくように、前記時間領域の両端となる時刻を決定するステップと、前記演算手段が、前記動揺モード検出の際に、逆時系列による入力データから検出される動揺モードと照合するステップと、前記演算手段が、前記自己回帰モデルの最適な次数を設定する際に、次数が探索範囲の上限値となるとき、周期が重複する動揺モードがなくなるまで次数を低減するステップとから構成されることを特徴とする請求項３に記載された電力系統の安定判別方法。
【請求項５】
前記請求項１、請求項２、請求項３又は請求項４の電力系統の安定判別方法のうちの少なくとも二つの安定判別方法を併用し、正常終了手段が、併用された少なくとも二つの判別結果の論理積を判別結果とするステップを備えることを特徴とする電力系統の安定判別方法。
【請求項６】
複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺を一定の時間領域でシミュレーションした結果について、この電力系統の安定度を自動判別するために、コンピュータを、
前記シミュレーションにおける解析対象時間内の結果に脱調する発電機を検出する異常検出手段、
前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力し、前記時系列の前半部及び後半部それぞれにおける前記入力データの振幅を算出し、前記前半部と後半部との振幅の比較結果を基準データとして電力系統の安定度の判別を行う第一安定判別手段、
前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力し、前記時系列の前半部及び後半部それぞれにおける入力データの波形の実効値を算出し、前記前半部と後半部との前記実効値の比較結果を基準データとして電力系統の安定度の判別を行う第二安定判別手段、
前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角を入力データとして時系列で入力し、前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルの最適な次数とパラメータを設定し、前記自己回帰モデルを用いて動揺モードの検出を行い、前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルによる動揺モード検出結果を基準データとして安定度の判別を行う第三安定判別手段、
当該第三安定判別手段において、前記動揺モード検出の際に、前記自己回帰モデルの同定に用いる前記入力データの時間領域の両端での値がその入力データの平均値に近づくように、前記時間領域の両端となる時刻を決定し、前記動揺モード検出の際に、逆時系列による入力データから検出される動揺モードと照合し、前記自己回帰モデルの最適な次数を設定する際に、次数が探索範囲の上限値となるとき、周期が重複する動揺モードがなくなるまで次数を削減する第四安定判別手段、
前記第一安定判別手段及び第二安定判別手段並びに第三安定判別手段又は第四安定判別手段のいずれかの安定判別処理において安定と判別された場合に電力系統が安定であると判別し安定判別処理を終了する正常終了手段、
前記シミュレーションにおける一定の時間領域内に前記異常検出手段により脱調する発電機が有る場合や前記シミュレーションで前記変数を導く演算が収束せずに異常終了した場合、或いは、前記第一安定判別手段乃至第四安定判別手段のいずれかの安定判別処理において不安定と判断された場合に、電力系統が不安定であると判別し安定判別処理を終了する異常終了手段として機能させるための電力系統の安定判別用プログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、電力系統の自動安定判別方法及びそのプログラムに関する。
【０００２】
【従来の技術】
従来の電力系統の自動安定判別には、発電機の位相角差が閾値を超えたことを検出する方法や電力動揺の単調増加を検出する方法がある（例えば、特許文献１参照）。また、エネルギー関数を応用したものなどがある（例えば、特許文献２参照）。特許文献１に開示された方法では、「発電機が脱調しているか否かの判別に、発電機の位相角差に対する一定の閾値を用いる。また、各時間刻みの位相角差を基に波形を認識し、波形データの極大値、極小値より振幅を求め、各振幅ごとに発散傾向（不安定）か収束傾向（安定）か否かを判別することもできる。また、各発電機の位相角差の代わりに、各発電機の角速度や加速エネルギーを用いて脱調判別を行なうことも可能である。」との記載がある。
【０００３】
【特許文献１】
特開２００１−３５２６７８号公報
【特許文献２】
特開平１１−２９９０９７号公報
【０００４】
また、大規模・複雑化した基幹系統では、電力系統の動揺が事故除去後も継続し、やがて脱調にいたる中間領域での安定度が問題となっているが、この対策として、複数の動揺モードが混在した中間領域の系統動揺を判別するために、プローニー解析法を用いた方法（例えば、特許文献３参照）や、時系列の有する動揺モードを時系列解析モデルにより検出する方法も提案されている（例えば、特許文献４参照）。
【０００５】
【特許文献３】
特開２００１−３５２６７９号公報
【特許文献４】
特開２０００−３２６６７号公報
【０００６】
【発明が解決しようする課題】
しかしながら、上記特許文献１及び２の、従来技術によれば、図２に示すような複雑な振動を伴う動揺波形に対して正確な判別が行われない場合がある。これを回避するための例外処理を行うと、例外処理のための閾値の設定や場合分けなどが複雑になるという問題がある。
【０００７】
また、図３のような複数の動揺モードが干渉する中間領域の判別も困難である。この対策として上記特許文献３及び４の各方法が有効と考えられるが、特許文献３の方法はノイズを含んだ波形の判別が難しい。一方特許文献４の方法は、ノイズには強いが、この方法を単独で用いた場合は、発電機が系統に接続されている位置によっては、誤不動作が発生しやすいという問題がある。
【０００８】
ここで、誤不動作とは、実際には不安定となっている波形に対して安定の判別を行うことをいう。一方、実際には安定となっている波形に対して不安定の判別を行う場合を誤動作という。
【０００９】
上記のように、単一の安定判別方法によって、あらゆる波形の安定／不安定を確実に判別することは、困難である。また、このような安定判別を用いて系統運用する際には、誤動作は許容されたとしても、誤不動作は許容できない。
【００１０】
そこで、本発明は、複雑な振動を伴う動揺波形に対して、極端な誤判別の生じない安定判別方法を得ること、複数の動揺モードが干渉する中間領域の波形において極端な誤判別の生じない安定判別方法を得ること、これらの安定判別方法を組み合わせる事で、誤不動作はせず、誤動作をできる限り抑制する安定判別方法を得ることを課題とする。
【００１１】
【課題を解決するための手段】
上記課題を解決するために為された本発明による電力系統の安定判別方法は、複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺を一定の時間領域でシミュレーションした結果について、この電力系統の安定度をコンピュータにより判別する方法である。
【００１２】
そして、第一安定判別方法は、入力手段が前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力するステップと、演算手段が前記時系列の前半部及び後半部それぞれにおける前記入力データの振幅、すなわち最大値と最小値の差を算出するステップと、判別手段が前記前半部と後半部との振幅の比較結果を基準データとして電力系統の安定度の判別を行うステップから構成されることを特徴とする電力系統の安定判別方法である。
【００１３】
第二安定判別方法は、入力手段が前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力するステップと、演算手段が前記時系列の前半部及び後半部それぞれにおける入力データの波形の実効値を算出するステップと、判別手段が前記前半部と後半部との前記実効値の比較結果を基準データとして電力系統の安定度の判別を行うステップから構成されることを特徴とする電力系統の安定判別方法である。
【００１４】
第三安定判別方法は、入力手段が前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角を入力データとして時系列で入力するステップと、演算手段が前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルの最適な次数とパラメータを設定すると共に、前記自己回帰モデルを用いて動揺モードの検出を行うステップと、判別手段が前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルによる動揺モード検出結果を基準データとして安定度の判別を行うステップから構成されることを特徴とする電力系統の安定判別方法である。
【００１５】
第四安定判別方法は、前記第三安定判別方法のステップに加え、前記演算手段が、前記動揺モード検出の際に、前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角の時系列による入力データの時間領域の両端での値がその入力データの平均値に近づくように、前記時間領域の両端となる時刻を決定するステップと、前記演算手段が、前記動揺モード検出の際に、逆時系列による入力データから検出される動揺モードと照合するステップと、前記演算手段が、前記自己回帰モデルの最適な次数を設定する際に、次数が探索範囲の上限値となるとき、周期が重複した動揺モードがなくなるまで次数を低減するステップとから構成されることを特徴とする電力系統の安定判別方法である。
【００１６】
そして、上記各安定判別方法の利点を併せ持つ構成として、前記第一安定判別方法乃至第四安定判別方法のうちの少なくとも二つの安定判別方法を併用する電力系統の安定判別方法であって、正常終了手段が、併用された少なくとも二つの判別結果の論理積を判別結果とするステップを備えることを特徴とする電力系統の安定判別方法とする場合もある。
【００１７】
更に、上記課題を解決する為に為された本発明による電力系統の安定判別用プログラムとして、複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺を一定の時間領域でシミュレーションした結果について、この電力系統の安定度を自動判別するために、以下の安定判別用プログラムが挙げられる。
【００１８】
このプログラムは、コンピュータを、異常検出手段、第一判別手段、第二判別手段、第三判別手段、及び第四判別手段、並びに正常終了手段、及び異常終了手段として機能させるものである。
【００１９】
前記異常検出手段は、前記シミュレーションにおける解析対象時間内に脱調する発電機や、前記シミュレーションで前記変数を導く演算が収束しないことによる異常終了を検出するものである。
【００２０】
また、前記第一安定判別手段は、前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力し、前記時系列の前半部及び後半部のそれぞれにおける前記入力データの振幅を算出し、前記前半部と後半部との振幅の比較結果を基準データとして電力系統の安定度の判別を行うものである。
【００２１】
また、前記第二安定判別手段は、前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力し、前記時系列の前半部及び後半部それぞれにおける入力データの波形の実効値を算出し、前記前半部と後半部との前記実効値の比較結果を基準データとして電力系統の安定度の判別を行うものである。
【００２２】
また、前記第三安定判別手段は、前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角を入力データとして時系列で入力し、前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルの最適な次数とパラメータを設定し、前記自己回帰モデルを用いて動揺モードの検出を行し、前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルによる動揺モード検出結果を基準データとして安定度の判別を行うものである。
【００２３】
更に、前記第四安定判別手段は、当該第三安定判別手段において、前記動揺モード検出の際に、前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角の時系列による入力データの時間領域の両端での値がその入力データの平均値に近づくように、前記時間領域の両端となる時刻を決定し、前記動揺モード検出の際に、逆時系列による入力データから検出される動揺モードと照合し、前記自己回帰モデルの最適な次数を設定する際に、次数が探索範囲の上限値となるとき、周期が重複する動揺モードがなくなるまで次数を削減するものである。
【００２４】
一方、前記正常終了手段は、前記第一安定判別手段及び第二安定判別手段並びに第三安定判別手段又は第四安定判別手段のいずれかの安定判別処理において安定と判別された場合に電力系統が安定であると判別し安定判別処理を終了するものであり、前記異常終了手段は、前記シミュレーションにおける解析対象時間内に前記異常検出手段により脱調する発電機が有る場合や前記シミュレーションで前記変数を導く演算が収束せず異常終了した場合、或いは、前記第一安定判別手段乃至第四安定判別手段のいずれかの安定判別処理において不安定と判断された場合に、電力系統が不安定であると判別し安定判別処理を終了するものである。
【００２５】
【発明の実施の形態】
以下、本発明による電力系統の安定判別方法の実施形態を図面に基づき説明する。
本発明による電力系統の安定判別方法は、複数の発電機を含む電力系統において、想定した事故が発生した後の電力系統の動揺について、発電機の運動方程式や電力系統の回路方程式を収束計算することで、各発電機等の状態を表す変数を求め、これを時間刻みごとに逐次繰り返すことで、シミュレーションを行い、その結果に基づき電力系統が安定であるか否かを判別するものであって、ＣＰＵや記憶手段等、各種ハードウエア資源からなるコンピュータを、異常検出手段、第一判別手段、第二判別手段、第三判別手段、及び第四判別手段、並びに正常終了手段、及び異常終了手段として機能させる電力系統の安定判別用プログラムによって制御されたものである。
【００２６】
前記異常検出手段は、安定判別処理の前処理として、前記シミュレーションにおける解析対象時間内での入力データから、脱調する発電機の有無や、計算が収束せずに終了する異常終了の有無を検出するステップ（ステップ１）を行うものである。
【００２７】
当該前処理が正常に終了した場合は、続いてコンピュータを第一判別手段として機能させ（ステップ２）、当該第一判別手段による判別結果が安定の場合は、続いてコンピュータを第二判別手段として機能させ（ステップ３）、当該第二判別手段による判別結果が安定の場合は、続いてコンピュータを第三判別手段又は第四判別手段として機能させ（ステップ４）、当該第三判別手段又は第四判別手段による判別結果が安定の場合は、コンピュータを前記正常終了手段と機能させ総合判別結果が安定であると判別し一連の安定判別方法を終了する（ステップ５）。
【００２８】
また、当該前処理（ステップ１）で脱調と判定された場合や解析が異常終了した場合、乃至はステップ２〜４のいずれかで、不安定の判別結果が得られた場合は、コンピュータを前記異常終了手段として機能させ、総合判別結果として電力系統が不安定であると取り扱い当該一連の安定判別方法を終了する（ステップ６）。
【００２９】
以下、各安定判別方法についてさらに詳しく説明する。前記第一安定判別手段は、入力手段が前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力するステップと、演算手段が前記時系列の前半部及び後半部それぞれにおける前記入力データの振幅、すなわち最大値と最小値の差を算出するステップと、判別手段が前記前半部と後半部との振幅の比較結果を基準データとして電力系統の安定度の判別を行うステップとを行う第一安定判別方法を実施するものである。
【００３０】
具体的には、例えば、各発電機の状態を表す変数として発電機角速度偏差を用い、図４にプロットされたような発電機角速度偏差波形について、波形の前半部の振幅Ａ、及び波形の後半部の振幅ＢについてＢ／Ａを算出し、Ｂ／Ａがある閾値（例えば０．９５）を超えた場合に不安定と判断すると共に、それ以外の場合は、安定と判断する。
【００３１】
この方法は、設定を容易に行うことができ、また計算も単純であり高速に判別が可能であるという長所をもつ。一方、複数の動揺モードが干渉する中間領域の波形に対して、誤不動作を生ずる恐れがある。これを極力回避するためには、判別の閾値を小さくすることが考えられるが、誤動作が増加するため、あまり閾値を小さくすることはできない。また、急峻なスパイクのようなノイズがある場合は、誤判別の可能性が高くなるという欠点がある。
【００３２】
なお、この例では、発電機角速度偏差波形について判別を行ったが、発電機内部相差角波形や、発電機加速エネルギー波形を用いて判別を行うことも可能である。
【００３３】
第二安定判別手段は、入力手段が前記シミュレーションで得られた各発電機の状態を表す変数を入力データとして時系列で入力するステップと、演算手段が前記時系列の前半部及び後半部それぞれにおける入力データの波形の実効値を算出するステップと、判別手段が前記前半部と後半部との前記実効値の比較結果を基準データとして電力系統の安定度の判別を行うステップとを行う第二安定判別方法を実施するものである。
【００３４】
具体的には、例えば、各発電機の状態を表す変数として発電機角速度偏差を用い、図４にプロットされたような発電機角速度偏差波形について、実効値に相当する数値の減衰傾向を判別するものである。具体的には、波形の前半部  （ｉ＝１とする）と後半部（ｉ＝２とする）それぞれにおいて（１）式によりＷｉを算出し、Ｗ２／Ｗ１がある閾値（例えば０．９５）を超えた場合に不安定と判断すると共に、それ以外の場合は、安定と判断する。
【００３５】
【数１】
2004282882.tif 000003
ただし、Ｎ：発電機角速度偏差［ＰＵ］、Ｎａｖｇ：発電機角速度偏差時系列データの平均値［ＰＵ］、Ｎのサンプリング刻み時間はΔｔ［ｓｅｃ］である。
【００３６】
ｔｉ及びＴには、一定の定数を用いても一応判別が可能であるが、より誤判別を少なくするために、ｔｉとしては、波形がＮａｖｇを横切った瞬間、Ｔはｔｉから数周期後において波形がＮａｖｇを横切った瞬間となるように波形に応じて可変の値とすることも可能である。
【００３７】
この方法は、第一安定判別方法の（ステップ２）と比較すれば、若干計算が複雑になるが、比較的単純であり高速に判別が可能であるという長所をもつ。また、スパイク状のノイズにも強い。一方、複数の動揺モードが干渉する中間領域の波形に対しては、前記第一安定判別方法と同様、誤不動作を生ずる恐れがある。
【００３８】
なお、当該第二安定判別方法（ステップ３）においても、発電機内部相差角波形や、発電機加速エネルギー波形を用いて判別を行うことも可能である。
【００３９】
前記第三安定判別手段は、入力手段が前記シミュレーションで得られた送電線両端の位相差及び各発電機の内部位相角を入力データとして時系列で入力するステップと、演算手段が前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルの最適な次数とパラメータを設定するステップ、及び前記自己回帰モデルを用いて動揺モードの検出を行うステップと、判別手段が前記送電線両端の位相差及び各発電機の内部位相角の自己回帰モデルによる動揺モード検出結果を基準データとして安定度の判別を行うステップとを行う第三安定判別方法を実施するものである。
【００４０】
当該第三安定判別方法は、前記演算手段による時系列解析によるものである。時系列解析とは、動揺する前記送電線両端の位相差及び各発電機の内部位相角の振動波形を、後述する次数の自動決定手法により、ｍ次の自己回帰モデル（（２）式参照）で表現して、得られた差分方程式をＺ変換して得られる特性方程式（（３）式参照）の特性根ｐｉを求めることで、動揺モードを検出する方法である。
【００４１】
【数２】
2004282882.tif 000004
ただし、ｙは出力信号（系統動揺）の時系列データ、ｗは白色雑音の時系列データ、ａｉは自己回帰モデルのパラメータ、ｍは自己回帰モデルの次数である。
【００４２】
【数３】
2004282882.tif 000005
ただし、ｐｉは特性根（実根か共役複素根），Ｄ（ｚ）は特性方程式である。
【００４３】
時系列解析の結果、共役複素根の特性根が存在する場合、その特性根で表現されるモードは振動的に変化する性質がある。さらにその共役複素根の絶対値が１より大きい場合、減衰率が正の不安定なモードが１つ存在し、絶対値が１より小さい場合は、減衰率が負の安定なモードが１つ存在する。絶対値がちょうど１の場合は、減衰率が０であり、持続振動となるモードが１つ存在する。一方、実根の特性根が存在する場合，その特性根で表現されるモードは一方向に単調に変化する性質がある。即ち、特性根が共役複素根のときに、その実部から動揺モードの減衰率が求まり，虚部から動揺モードの周期が求まることとなる。
【００４４】
本発明では、系統状態や事故点などで大きく変化する動揺モードを精度よく検出するために、前記の如く自己回帰モデルの最適な次数とパラメータを自動設定する特徴を有する。この具体的な手順について以下に述べる。
【００４５】
図５に示すフローチャートのように、まず、前記入力手段としての機能により時系列の入力データを取り込む。次に、前記演算手段としての機能によりモデルの次数を探索範囲内で順次更新しながら、各次数のモデルにおけるパラメータおよび赤池のＦＰＥ（（４）式参照）を算出する。前記自己回帰モデルの次数を更新する探索範囲は、Ｎを時系列データ数として、１〜３√Ｎである。赤池のＦＰＥ（（４）式参照）を用いて同定精度を診断し，赤池のＦＰＥが最小となる次数の自己回帰モデルを用いて、動揺モードの解析を実施する。次に、検出された動揺モードを診断し、その変動量が時系列解析で得られる実根よりも大きく、且つその周期が、時系列の入力データのサンプリング時間の５倍から解析対象時間までの範囲にある動揺モードを有効な動揺モードと特定する。
【００４６】
【数４】
2004282882.tif 000006
ただし、ｎは時系列データの数、ｍは自己回帰モデルの次数、σ２は自己回帰モデル誤差の分散である。
【００４７】
即ち、前記ステップ４の処理の全容は図６のようになる。まずシミュレーションで得られた送電線両端の位相差について、前記演算手段としての機能による時系列解析により動揺モードの検出を行い、不安定モードを特定する。次に各発電機の内部位相角の時系列データについても動揺モードの検出を行い、不安定モードを特定する。そして、前記判別手段の機能により上記動揺モード検出結果を基準データとして送電線の不安定モードと一致する不安定モードが発電機に存在する場合は不安定と判別し、存在しない場合は安定と判別する。
【００４８】
また、送電線で減衰率が正の動揺モードが検出されない場合には安定と判別する代わりに、送電線で減衰率が正の動揺モードが検出されない場合、多数の発電機で検出される動揺モードの中に、減衰率が正で、送電線で検出される動揺モードと同一周期の動揺モードが存在すれば不安定と判別することで、誤不動作の抑制を図ることもできる。ただし、この論理を追加すると、誤動作が増加する傾向があるため、ステップ２乃至ステップ３の方法と組み合わせて判別を実施する場合は、この論理の追加は必要ない。
【００４９】
当該ステップ４において採用される処理としては、以上述べたような時系列解析法を用いる代わりに、上記第三安定判別方法とは以下の点で異なる第四安定判別方法を用いることもできる。
【００５０】
第四安定判別方法は、前記第三安定判別方法のステップに加え、前記演算手段が、前記自己回帰モデルにとって誤差の要因となる不要なデータを取り除くため、前記動揺モード検出の際に、前記自己回帰モデルの同定に用いるデータの時間領域の両端での値がそのデータの平均値に近づくように、前記時間領域の両端となる時刻を決定するステップを付加することができる。
【００５１】
また、前記演算手段が、発散傾向のデータの判別精度を向上させるため、前記動揺モード検出の際に、図７のフローチャートに示すように、正方向で特定された動揺モードに加えて、逆方向で特定される動揺モードと照合するステップを加えることができる。具体的には、図７に示すフローチャートのように、送電線両端の位相差の時系列データから極めて減衰の悪い動揺モード（以下、動揺モードＡと記す）が検出された場合、逆方向の時系列データを用いて動揺モードを検出する。前記動揺モードＡと同一周期の逆方向の動揺モード（以下、動揺モードＢと記す）が存在する場合は、この逆方向のデータから得られる動揺モードＢの減衰率と元の動揺モードＡの減衰率を比較し、同一の時間軸方向で見て、前者の減衰の方が悪ければ、前者の減衰率を用いて判別を行う。
【００５２】
更に、前記演算手段が、自己回帰モデルの最適な次数を設定する際に、赤池のＦＰＥによる基準値を評価関数として求められる次数が探索範囲の上限値となるとき、周期の重複した動揺モードが検出されて誤判別となることが多いため、その場合には、周期が重複する動揺モードがなくなるまで次数を低減するステップを加えることができる。
【００５３】
具体的には、図８に示すフローチャートのように、周期の重複した動揺モード（以下，重複モードと記す）の数Ｎを算出し、２Ｎだけ次数を低減して再度動揺モード検出を実施する。その後、重複モードの有無を確認し、重複モードが無ければその動揺モードを採用、重複モードがあれば再度重複モードの数を算出して、先の手順に戻り、重複モードがなくなるまで繰り返す。
【００５４】
最後にステップ５及び６では、ステップ２から４の判別結果を集約して最終判別結果を出力するが、当該例では、ステップ２から４のすべての判別結果が安定となった場合に限って、前記正常終了手段の機能により電力系統が安定であると判別され安定判別処理を終了するものであり、それ以外の場合、即ち、前記異常検出手段による前処理中に前記異常検出手段により脱調する発電機が有る場合や前記シミュレーションで前記変数を導く演算が収束せず異常終了した場合、或いは、前記第一安定判別手段乃至第三又は第四安定判別手段のいずれかの安定判別処理において不安定と判断された場合は、前記異常終了手段の機能により電力系統が不安定であると判別し安定判別処理を終了する。
【００５５】
これは、例えば、この判別結果によって、系統安定化を実施する場合に、誤動作によって安定化制御が若干増加してしまうことは止むを得ないが、誤不動作によって系統が崩壊してしまう事態は許容しがたいために、このような方法を採っているものである。しかしながら、逆に誤不動作が許容される場合は、ステップ２から４のすべての判別が不安定となった場合に限って、総合判別が不安定となり、それ以外の場合を安定とする方法とすることも可能である。
【００５６】
また、必ずしもステップ２から４の全てを実施する必要はなく、それぞれのステップを単独で用いて安定判別を実施すること、もしくは、任意の２つのステップを組み合わせて安定判別を実施することも、もちろん可能である。この場合は、必要に応じて閾値を厳しめの値とする等の対策により、誤不動作が生じないようにすることが可能である。
【００５７】
さらに、図１のステップ２から４の実施順序は、必ずしもこの通りである必要はなく、例えばステップ４の一部（送電線の動揺モードチェック）のみを先に実施してこの収束性が良い場合は、安定と判別し、それ以外の場合について、ステップ２、３を実施するという方法も可能である。
【００５８】
表１及び表２に、本発明による電力系統の安定判別方法を適用した結果の一例を示す。
表１は、図１における、ステップ１が終了した、総数２５１ケースのシミュレーション結果について図１における、ステップ２及び３を単独で適用した結果をまとめた例である。ステップ２、３では、誤不動作が少ない結果が得られる事がわかる。ただし、図３のような波形の判別は困難である。
【００５９】
【表１】
2004282882.tif 000007
【００６０】
表２は、図１における、ステップ１が終了した、総数２５１ケースのシミュレーション結果について図１における、ステップ４を単独で適用した結果をまとめた例である。ステップ４では、誤動作が少ない結果が得られる事がわかる。このため、ステップ２、３の判別と組み合わせて用いたとしても、総合判別において、誤動作が増加する恐れが少ない。
【００６１】
【表２】
2004282882.tif 000008
【００６２】
【発明の効果】
以上説明したように、前記第一安定判別方法及び第二安定判別方法によれば、複雑な振動を伴う動揺波形に対して、極端な誤判別の生じることがなく、かつ、閾値設定が容易であり、計算も比較的単純であり高速な安定判別が可能である。また、第三安定判別方法及び第四安定判別方法によれば、複数の動揺モードが干渉する中間領域の波形に対しても、極端な誤判別の生じない安定判別が可能である。さらに、それらを組み合わせる事で、誤不動作はせず、誤動作をできる限り抑制する電力系統の安定判別方法が得られる。
【図面の簡単な説明】
【図１】本発明による電力系統の安定判別方法の一例を示すフローチャートである。
【図２】複雑な振動を伴う動揺波形の例を示す図である。
【図３】複数の動揺モードが干渉する波形の例を示す図である。
【図４】図１中のステップ１の適用例の説明図である。
【図５】本発明による電力系統の安定判別方法における自己回帰モデルの次数を自動決定する手順を説明するフローチャートである。
【図６】図１中のステップ４の具体的な算出手順を説明するフローチャートである。
【図７】本発明による電力系統の安定判別方法における正方向と逆方向の時系列データから検出される動揺モードを照合する方法を説明するフローチャートである。
【図８】本発明による電力系統の安定判別方法における重複モードがなくなるまで次数を低減する手順を説明するフローチャートである。
【図１】
2004282882.tif 000009
【図２】
2004282882.tif 000010
【図３】
2004282882.tif 000011
【図４】
2004282882.tif 000012
【図５】
2004282882.tif 000013
【図６】
2004282882.tif 000014
【図７】
2004282882.tif 000015
【図８】
2004282882.tif 000016