(57)【要約】      （修正有）【課題】発電機の動揺変化が同じ傾向にあるものをグルーピングし、代表指定した発電機は詳細な安定度計算を実施し、同グループに属する他の発電機は、換算係数を用いて簡易計算する、電力系統解析装置およびその方法を提供する。【解決手段】電力系統解析装置は、複数の発電機の相差角の変化傾向に基づいてグループ分けを行うグルーピング作成部１０２と、グループにおける代表発電機を決定する代表発電機決定部１０３と、グループ及び代表発電機の情報に基づいて解析対象の系統モデルのデータを作成する系統モデル作成部１０５と、代表以外の発電機の相差角の変化量を求めるための換算係数を算出する換算係数算出部と、系統モデル及び換算係数に基づいて系統事故を想定した場合の系統の安定度を計算する安定度計算部１０６とを備える。【選択図】図１
2017118614.tif 000002
【特許請求の範囲】
【請求項１】
複数の発電機の相差角の変化傾向に基づいてグループ分けを行うグルーピング作成部と、
前記グループにおける代表発電機を決定する代表発電機決定部と、
前記グループ及び前記代表発電機の情報に基づいて解析対象の系統モデルのデータを作成する系統モデル作成部と、
前記代表以外の発電機の相差角の変化量を求めるための換算係数を作成する換算係数作成部と、
前記系統モデル及び前記換算係数に基づいて系統事故を想定した場合の系統の安定度を計算する安定度計算部と、を備えることを特徴とする電力系統解析装置。
【請求項２】
請求項１に記載の電力系統解析装置において、
前記安定度計算部は、安定度計算した代表発電機の相差角の変化量及び前記代表以外の発電機の換算係数を用いて、前記代表以外の発電機の相差角の変化量を計算することを特徴とする電力系統解析装置。
【請求項３】
請求項１又は請求項２の何れかに記載の電力系統解析装置において、
前記グルーピング作成部は、発電機毎の容量、単位慣性定数、又は発電機間の距離に応じて、複数の発電機をグルーピングすることを特徴とする電力系統解析装置。
【請求項４】
請求項１に記載の電力系統解析装置において、
前記代表発電機決定部は、発電機毎の容量、単位慣性定数、又は定格出力に基づいて代表の発電機を決定することを特徴とする電力系統解析装置。
【請求項５】
請求項１に記載の電力系統解析装置において、
前記安定度計算部は、前記代表発電機決定部で決定した代表発電機の動揺を詳細な時間領域シミュレーションすることを特徴とする電力系統解析装置。
【請求項６】
請求項１乃至請求項５の何れかに記載の電力系統解析装置において、
前記安定度計算部は、前記代表発電機の相差角の変化量を求めることを特徴とする電力系統解析装置。
【請求項７】
請求項１乃至請求項６の何れかに記載の電力系統解析装置において、
前記安定度計算部は、前記代表発電機の相差角の変化量と前記代表以外の発電機の換算係数が線形関係であることに基づいて、前記代表以外の発電機の相差角変化量を求めることを特徴とする電力系統解析装置。
【請求項８】
請求項１乃至請求項７の何れかに記載の電力系統解析装置において、
前記安定度計算部は、前記代表発電機の相差角の変化量と前記代表以外の発電機の換算係数が非線形関係であることに基づいて、前記代表以外の発電機の相差角変化量を求めることを特徴とする電力系統解析装置。
【請求項９】
請求項１に記載の電力系統解析装置において、
前記換算係数作成部は、前記換算係数発電機の容量、単位慣性定数、又は定格出力に基づいて前記代表以外の発電機の換算係数を求めることを特徴とする電力系統解析装置。
【請求項１０】
請求項１に記載の電力系統解析装置において、
前記安定度計算部は、発電機の相差角を計算する途中で、前記グループを修正することを特徴とする電力系統解析装置。
【請求項１１】
請求項１乃至請求項１０の何れかに記載の電力系統解析装置において、
前記安定度計算部は、発電機の相差角を計算する途中に、所定時間前の時刻における発電機の相差角の変化量の大きさに基づいて、前記グループを修正することを特徴とする電力系統解析装置。
【請求項１２】
請求項１に記載の電力系統解析装置において、
前記安定度計算部は、発電機の相差角を計算する途中に、同グループの代表発電機を修正することを特徴とする電力系統解析装置。
【請求項１３】
請求項１乃至請求項１２の何れかに記載の電力系統解析装置において、
前記安定度計算部は、発電機の相差角を計算する途中に、所定時間前の時刻における発電機の相差角の変化量の大きさに基づいて、同グループの代表発電機を修正することを特徴とする電力系統解析装置。
【請求項１４】
請求項１に記載の電力系統解析装置は、
前記安定度計算部で計算した系統の安定度に基づいて、系統の安定性に与える影響を評価する安定度評価部を更に備えることを特徴とする電力系統解析装置。
【請求項１５】
複数の発電機の相差角の変化傾向に基づいてグループ分けを行い、前記グループにおける代表発電機を決定し、前記グループ及び前記代表発電機の情報に基づいて解析対象の系統モデルのデータを作成し、前記代表以外の発電機の相差角の変化量を求めるための換算係数と前記系統モデルとに基づいて系統事故を想定した場合の系統の安定度を計算することを特徴とする電力系統解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統の安定度解析を行う電力系統解析装置及び方法に関する。
【背景技術】
【０００２】
電力会社は、電力系統の安定性を確保するため、事前に想定事故ケースを解析して系統が不安定にならないように安定度解析装置や安定化システムを用いて系統内の発電機や機器を制御する。解析対象の系統規模拡大や分散電源の導入量増加に伴い、安定度解析の計算時間が非常に長くかかってしまう。そのため、解析時間を短縮するための高速化手法の一つとして、電力系統を縮約する技術がある。本技術分野の背景技術として、特開２００７−１８１３８７号公報（特許文献１）がある。この公報には、想定事故毎に安定度シミュレーションを繰り返し、等価な２機系モデルを作成して安定度を判定する、という内容が記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００７−１８１３８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１には、電力系統の過渡安定度判定方法およびその装置が記載されている。しかし、系統を縮約することで安定度解析を高速化できるが、縮約後には下位系統の詳細な情報が失われてしまう。 
【０００５】
そこで、本発明は、発電機の動揺変化が同じ傾向にあるものをグルーピングし、各発電機に関する情報に基づいて代表指定した発電機は詳細な安定度計算を実施して発電機相差角の変化量を求め、同じグループに属する他の発電機は、その変化量と発電機毎に与えられた換算係数を用いて動揺変化を簡易計算することができる、電力系統解析装置およびその方法を提供する。
【課題を解決するための手段】
【０００６】
上記課題を解決する為に本発明に係る電力系統解析装置は、複数の発電機の相差角の変化傾向に基づいてグループ分けを行うグルーピング作成部と、前記グループにおける代表発電機を決定する代表発電機決定部と、前記グループ及び前記代表発電機の情報に基づいて解析対象の系統モデルのデータを作成する系統モデル作成部と、前記代表以外の発電機の相差角の変化量を求めるための換算係数を作成する換算係数作成部と、前記系統モデル及び前記換算係数に基づいて系統事故を想定した場合の系統の安定度を計算する安定度計算部と、を備えることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、電力系統の安定度解析において下位系統の情報を失わずに、解析時間を短縮できる。上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
【図面の簡単な説明】
【０００８】
【図１】電力系統解析装置の構成図の例である。
【図２】電力系統解析装置と電力系統の関係図の例である。
【図３】発電機のグルーピングに関するデータの一例を示す図である。
【図４】発電機をグループ分けしたときの構成図の例である。
【図５】発電機情報記憶部の一例を示す図である。
【図６】代表の発電機決定部の処理の全体を示すフローチャートの例である。
【図７】安定度計算部の処理の全体を示すフローチャートの例である。
【図８】安定度計算部の処理で代表以外の発電機に関して簡易動揺計算の処理を示すフローチャートの例である。
【図９】発電機毎の換算係数の一例を示す図である。
【図１０】実施例２における電力系統解析装置の構成図の例である。
【図１１】換算係数算出部の処理の全体を示すフローチャートの例である。
【図１２】実施例３における電力系統解析装置の構成図の例である。
【図１３】実施例３における安定度計算部の処理の全体を示すフローチャートの例である。
【発明を実施するための形態】
【０００９】
以下、本発明の実施例を、図面を用いて説明する。尚、下記はあくまでも実施の例に過ぎず、下記具体的内容に発明自体が限定されることを意図する趣旨ではない。
【実施例１】
【００１０】
図１は、本発明の一実施形態に係る電力系統解析装置の構成図の例である。
【００１１】
図１に示す通り、電力系統解析装置１００では、データ入力部１０１、グルーピング作成部１０２、代表の発電機決定部１０３、発電機情報記憶部１０４、系統モデル作成部１０５、安定度計算部１０６、換算係数入力部１０７、安定度評価部１０８、出力部１０９を有する。
【００１２】
データ入力部１０１は、解析対象の系統に関する情報や、発電機に関する情報、発電機のグルーピング作成に必要な情報などを外部入力する。例えば、系統の送電線インピーダンスや発電機の単位慣性定数、発電機容量などである。発電機のグルーピングに必要な情報として、発電機毎に指定する値を直接入力しても良いし、グルーピングの指標となる発電機容量や単位慣性定数、定格出力などの値を入力しても良い。
【００１３】
グルーピング作成部１０２は、発電機のうち、相差角の動揺変化が同じ傾向のものを同一グループとなるように、複数のグループを作成する。例えば、データ入力部１０１で入力されたグループ分けのテーブル情報に基づいてグルーピングを実施する。この入力データは、グルーピングされた情報でも良いし、グルーピングに必要な指標でも良い。グルーピングの指標は、例えば、発電機間の距離や発電機容量などである。
【００１４】
代表の発電機決定部１０３は、グルーピング作成部１０２で作成されたグループ情報を基に、各グループの代表となる発電機を決定する。発電機に関する情報は発電機情報記憶部１０４に記憶しておくことで、繰り返し同じ電力系統を解析する場合にデータ入力を簡略化し、解析時間を短縮することができる。なお、代表の発電機決定方法の一例は後述する。
【００１５】
発電機情報記憶部１０４は、代表の発電機を決定するために必要な発電機毎の情報を記憶する。発電機の情報は、例えば、発電機容量、単位慣性定数、定格出力、制御系のパラメータなどである。
【００１６】
系統モデル作成部１０５は、データ入力部１０１と、グルーピング作成部１０２と、代表の発電機決定部の各種情報を用いて、解析対象の電力系統モデルのデータを整理し、電力系統モデルのデータ一式を作成する。電力系統モデルのデータは、例えば、発電機に関する情報、負荷に関する情報、送電線に関する情報などである。
【００１７】
安定度計算部１０６は、系統モデル作成部１０５で作成した電力系統モデル情報を用いて、系統事故を想定した場合の電力系統の安定性への影響を時間領域シミュレーションで計算する。なお、時間領域シミュレーションの計算方法の一例は、後述する。
【００１８】
換算係数入力部１０７は、安定度計算部１０６で代表以外の発電機の相差角変化量を簡易計算するときに用いるための換算係数を外部から入力する。データの一例は後述する。
【００１９】
安定度評価部１０８は、安定度計算部１０６で得られた電力系統の安定性への影響を評価する。
【００２０】
出力部１０９は、安定度計算および安定度評価の結果を外部に出力する。出力は、パソコンのディスプレイに出力しても良いし、紙に印刷しても良い。
【００２１】
図２は、電力系統解析装置と電力系統の関係を示した図の例である。
【００２２】
図２に示す通り、電力系統解析装置１００は、電力系統安定化制御装置２００と通信ネットワーク２０１を介して電力系統２０２や発電機２０３−２０８と接続する。発電された電力は、電力系統２０２を介して負荷２０９が消費する。電力系統解析装置１００は、安定度計算を実施し、その結果を評価する。電力系統安定化制御装置は、その評価結果を基づいて、電力系統の安定性を保つように各発電機や電力系統に設置されている機器に制御指令を送ることにより、電力系統の安定性を確保する。
【００２３】
電力系統安定化制御装置２００は、電力系統解析装置の結果に基づいて、電力系統の安定性を保つように各発電機や電力系統に設置されている機器に制御指令を送る。
【００２４】
通信ネットワーク２０１は、電力系統２０２や発電機２０３−２０８と、電力系統安定化制御装置２００を接続する。電力系統２０２の系統状態に関する情報や、発電機２０３−２０８に関する情報が、電力系統安定化制御装置２００と電力系統解析装置１００に送られる。一方、電力系統安定化制御装置２００は、通信ネットワーク２０１を介して、制御指令を電力系統２０２や発電機２０３−２０８に送信する。
【００２５】
電力系統２０２は、発電機と負荷を接続するネットワークである。例えば、送電線や変圧器などで構成されている。
【００２６】
発電機２０３−２０８は、電力系統に電力を供給する。発電機は６台に限るものではなく、１台でも良いし１０台以上でも良い。
【００２７】
負荷２０９は、電力系統２０２を介して発電機２０３−２０８で発電された電力を受け取り、消費する。負荷は一つに限るものではなく、二つでも良いし、それ以上でも良い。
【００２８】
図３は、発電機のグルーピングに関するデータの一例を示す。例えば、データ入力部１０１から発電機毎のグループＩＤを入力すると、図３に示すようなテーブルがグルーピング生成部１０２で作られる。グループＩＤは、データ入力部から直接、値を入力しても良いし、データ入力部に入力された情報で例えば発電機の単位慣性定数や発電機間の距離などから決定しても良い。
【００２９】
図４は、発電機をグループ分けしたときの構成図の例である。既に説明した図２に示された同一の符号を付した構成と、同一の機能を有する部分については、説明を省略する。例えば、図３に示したグルーピング結果の一例に基づいて、発電機グループを２つのグループに分けた場合の構成図を示している。これは、発電機グループを２つに限定するものではなく、３つでも良いし、それ以上でも良い。図４に示した発電機グループ毎に、代表の発電機決定部１０３で代表となる発電機を選定する。
【００３０】
図５は、発電機情報記憶部１０４に記憶されるデータテーブルを示す図の例である。各発電機に関する情報として、発電機容量や単位慣性定数、定格出力などが格納される。このとき、数値の単位は実際の値でも良いし、規格化された値でも良い。
【００３１】
図６は、代表の発電機決定部１０３の処理の全体を示すフローチャートの例である。グルーピング作成部１０２で作成された発電機のグループ毎に代表の発電機を選定する。選定には、データ入力部１０１で入力された発電機に関する情報を基に選定する。このとき、選定の指標には、発電機容量でも良いし、発電機の単位慣性定数でも良い。ここでは、発電機容量を指標として、発電機容量が最も大きなものが代表となるように選定する場合について説明する。
【００３２】
ステップＳ１０１では、発電機情報である発電機ＩＤやグルーピング作成部１０２で付与された発電機グループＩＤ、発電機容量などの情報を入力する。
【００３３】
ステップＳ１０２では、発電機グループＩＤ毎に代表を選定するため、発電機グループＩＤであるＧＩＤを初期化する。
【００３４】
ステップＳ１０３では、データ番号ｎを初期化する。以下、データ番号が終了するまでステップＳ１０９までの処理を繰り返し実行する。
【００３５】
ステップＳ１０４では、発電機容量の最大値Ｃ＿ｍａｘ（ＧＩＤ）を初期化する。ここでは、発電機容量が最大のものを代表の発電機とするため、Ｃ＿ｍａｘ（ＧＩＤ）を０と設定するように初期化する。
【００３６】
ステップＳ１０５では、発電機グループＩＤが同一のものか判定する。ＧＩＤが同一であれば、ステップＳ１０６に進む。もし、ＧＩＤが同一でなければ、Ｓ１０８に進んで次のデータを参照するように処理する。
【００３７】
ステップＳ１０６では、発電機容量が最大値かどうかを判定する。発電機容量が最大である場合には、ステップＳ１０７に進む。もし、発電機容量が最大でなければ、ステップＳ１０８に進んで次のデータを参照するように処理する。
【００３８】
ステップＳ１０７では、発電機容量をＣ＿ｍａｘ（ＧＩＤ）に記憶する。ステップＳ１０５とステップＳ１０６の判定を通過した場合は、同一の発電機グループで発電機容量が最大であるため、発電機容量をＣ＿ｍａｘ（ＧＩＤ）に記憶する。
【００３９】
ステップＳ１０８では、データ番号を一つ増やし、次のデータを参照するように処理を進める。
【００４０】
ステップＳ１０９では、データ番号が終了したかどうか判定する。もしデータ番号が終了していなければ、ステップＳ１０４に戻って処理を繰り返す。もしデータ番号が終了している場合は、同一の発電機グループでの代表選定が終了したことになる。そのため、次の発電機グループＩＤで代表を選定する処理に移行する。
【００４１】
ステップＳ１１０では、次の発電機グループの代表を選定するためにグループＩＤを進める。
【００４２】
ステップＳ１１１では、全ての発電機グループで処理を実施したか判定する。もし全ての発電機グループで選定が終了していなければ、ステップＳ１０３に戻って処理を繰り返す。もし全ての発電機グループで代表を選定していれば、処理は全て完了となる。
【００４３】
以上のように発電機グループ毎に代表を決定し、系統モデル作成部１０５で解析する電力系統のデータを作成する。
【００４４】
図７は、安定度計算部１０６の処理において、代表発電機の時間領域シミュレーションの処理を示すフローチャートの例である。安定度計算部１０６では、発電機グループ内の代表の発電機について、想定事故時の時間領域シミュレーションを実施する。これにより、発電機の相差角の時間変化が得られる。さらに、このとき得られた相差角の変化量と、発電機毎に与えられた換算係数に基づいて、同一の発電機グループ内に属する他の発電機の相差角の動揺を簡易的に求める。図７は前者の時間領域シミュレーション全体の処理を示したものである。後者の簡易動揺計算の処理は、図８に示し、後述する。
【００４５】
ステップＳ２０１は、初期潮流状態を算出する。初期潮流計算は、ニュートンラフソン法を用いても良いし、直流法を用いても良い。
【００４６】
ステップＳ２０２は、想定する事故条件を入力する。例えば、２回線送電線＃１で三相地絡事故が発生し、事故から０．１秒後に送電線＃１の１回線を開放することによって事故を除去する。その後事故から０．２秒後に再閉路して系統を元の状態に戻す、といった事故を想定する。そのとき、事故除去・再閉路のタイムスケールや、系統状態の変化に伴う送電線の線路インピーダンス情報を入力する。
【００４７】
ステップＳ２０３は、時間領域シミュレーションの開始時刻を設定する。以下、ステップＳ２０９までの処理を繰り返すことで時間領域シミュレーションを実行する。
【００４８】
ステップＳ２０４は、系統状態の変更を判定する。ステップＳ２０２で入力した事故条件に基づいて、事故除去や再閉路に伴って系統状態が変化する場合、ステップＳ２０５に進んで系統状態を変更する。もし、次の時刻でも系統状態が変化しない場合には、ステップＳ２０６に進むように処理する。
【００４９】
ステップＳ２０５は、ステップＳ２０４で系統状態の変更の必要があると判定された場合に、系統状態を変更する。例えば、線路インピーダンスの値を変更する。
【００５０】
ステップＳ２０６は、系統計算を実行する。系統状態や発電機端子電圧の状態に応じて系統の線路潮流やノード電圧が変化する。そのため、潮流計算を実施して系統の状態を求める。このとき、潮流計算はニュートンラフソン法でも良いし、直流法でも良い。ただし、計算が収束しない場合には数値解法を変更して解を求める。
【００５１】
ステップＳ２０７は、積分計算を実施する。時刻が進んで系統状態が変化することで、発電機の動揺がどのように変化するのか、発電機の動揺方程式を解いて相差角の変化量を求める。このとき、数値解法はルンゲクッタ法でも良いし、台形法でも良い。
【００５２】
ステップＳ２０８は、ステップＳ２０７で求めた代表発電機の相差角変化量と換算係数を用いて、代表以外の発電機の動揺を簡易的に算出する。処理の詳細は、図８で後述する。
【００５３】
ステップＳ２０９は、時間刻みΔｔだけ時間を進める。
【００５４】
ステップＳ２１０は、時間領域シミュレーションの終了時刻かどうか判定する。もしシミュレーションの終了時刻ではない場合は、ステップ２０４に戻って処理を繰り返す。もしシミュレーションの終了時刻であれば、処理を終了する。
【００５５】
以上の処理を実施することで、発電機グループ内の代表の発電機について、想定事故時の時間領域シミュレーションを実施することができる。
【００５６】
図８は、安定度計算部１０６の処理で代表以外の発電機に関して簡易動揺計算の処理を示すフローチャートの例である。
【００５７】
ステップＳ３０１は、ステップＳ２０７で算出した代表発電機の相差角変化量を入力する。
【００５８】
ステップＳ３０２は、発電機の換算係数を入力する。換算係数は、発電機毎の違いを示す指標であり、例えば発電機容量や発電機の単位慣性定数、定格出力などの比がそれに当たる。
【００５９】
以上の２つの情報から、代表発電機以外の発電機の動揺を、以下に示す処理で求める。
【００６０】
ステップＳ３０３は、発電機ＩＤを初期化する。これは、全ての発電機で計算するためである。
【００６１】
ステップＳ３０４は、代表発電機か否か判定する。代表発電機であれば計算の必要が無いため、ステップ３０６に進む。もし代表発電機でなければ、ステップＳ３０５に進む。
【００６２】
ステップＳ３０５は、代表発電機以外の相差角の動揺を求める。例えば、簡易動揺計算の一例を示すとするならば、代表発電機の相差角変化量をΔδ、代表発電機以外の相差角変化量をΔδ’、換算係数をＫとすると、次のように計算することができる。
【００６３】
（数１）
Δδ’＝Ｋ×Δδ
【００６４】
ただし、換算係数は比例関係に限らず、線形関係でも良いし、係数の二乗に比例してもよいし、Δδの二乗に比例しても良い。
【００６５】
ステップＳ３０６は、発電機ＩＤを次のものに進める。
【００６６】
ステップＳ３０７は、全ての発電機の計算ができたか判定する。もし全ての発電機の計算ができていなければ、ステップＳ３０４に戻って処理を繰り返す。もし全ての発電機の計算が終了していたら、処理を終える。
【００６７】
以上のようにして代表発電機以外の相差角変化量の計算を終えたら、ステップ２０９に戻り、次の時刻のシミュレーション処理を進める。
【００６８】
図９は、発電機毎の換算係数の一例を示す図である。換算係数は、発電機毎に入力される値である。
【００６９】
以上のようにして、処理を実施することにより、電力系統の安定度解析において下位系統の情報を失わずに、解析時間を短縮することができる。
【実施例２】
【００７０】
本実施例では、換算係数算出部３０１で発電機の換算係数を自動的に求めて安定度計算を実施する電力系統解析装置３００の例を説明する。
【００７１】
図１０は、実施例２における電力系統解析装置の構成図の例である。
【００７２】
図１０に示す通り、電力系統解析装置３００では、データ入力部１０１、グルーピング作成部１０２、代表の発電機決定部１０３、発電機情報記憶部１０４、系統モデル作成部１０５、安定度計算部１０６、安定度評価部１０８、出力部１０９に加えて、換算係数算出部３０１、換算係数記憶部３０２を有する。既に説明した図１に示された同一の符号を付した構成と、同一の機能を有する部分については、説明を省略する。
【００７３】
換算係数算出部３０１は、発電機情報記憶部１０３に記憶された発電機に関する情報を用いて換算係数を算出する。算出方法の例については、図１１を用いて後述する。算出した換算係数は、安定度計算部１０６に出力する。
【００７４】
換算係数記憶部３０２は、換算係数算出部３０１で計算した換算係数を格納する。また、換算係数を換算係数算出部３０１で算出する場合に、過去の計算結果を参照してもよい。
【００７５】
図１１は、換算係数算出部の処理の全体を示すフローチャートの例である。
【００７６】
ステップＳ４０１では、発電機に関する情報を発電機情報記憶部１０４から入力する。発電機に関する情報は、例えば、発電機容量や単位慣性定数、定格出力などである。
【００７７】
ステップＳ４０２では、データ番号ｎを初期化する。以下、ステップＳ４０５までの処理を繰り返すことによって、全てのデータについて換算係数を算出する。
【００７８】
ステップＳ４０３では、発電機の情報に基づいて換算係数を算出する。算出する方法は、例えば、発電機容量を用いて算出する方法を説明すると、以下の式を用いて算出することができる。
【００７９】
（数２）
Ｋ＝Ｃａｐ（ｎ）／Ｃａｐ（１）
【００８０】
ここで、Ｃａｐ（ｎ）は、データ番号ｎの発電機容量、Ｃａｐ（１）は代表発電機の発電機容量を示している。このように、例えば、換算係数は発電機容量の比を用いて算出することができる。これは、発電機容量だけでも良いし、単位慣性定数だけでも良い。また、発電機容量と単位慣性定数を組み合わせても良い。
【００８１】
ステップＳ４０４では、データ番号を進める。
【００８２】
ステップＳ４０５では、データ番号が終了したかどうかを判定する。もしデータ番号が終了していなければ、ステップＳ４０３に戻って処理を進める。もしデータ番号が終了した場合は、全ての換算係数の算出が完了したため、処理を終了する。
【実施例３】
【００８３】
本実施例では、安定度計算部１０６で安定度を解析する場合に、発電機のグルーピングや換算係数を計算途中で修正することができる電力系統解析装置４００の例を説明する。
【００８４】
図１２は、実施例３における電力系統解析装置の構成図の例である。
【００８５】
図１２に示す通り、電力系統解析装置４００では、データ入力部１０１、グルーピング作成部１０２、代表の発電機決定部１０３、発電機情報記憶部１０４、系統モデル作成部１０５、安定度評価部１０８、出力部１０９、換算係数算出部３０１、換算係数記憶部３０２に加えて、安定度計算部４０１を有する。既に説明した図１および図１０に示された同一の符号を付した構成と、同一の機能を有する部分については、説明を省略する。
【００８６】
安定度計算部４０１では、系統モデル作成部１０５で作成した電力系統モデル情報を用いて、系統事故を想定した場合の電力系統の安定性への影響を時間領域シミュレーションで計算する。このとき、時間領域シミュレーションにおいて、発電機のグルーピングや換算係数を計算途中で修正するため、グルーピング作成部１０２や換算係数算出部３０１に修正を指示することができる。なお、時間領域シミュレーションの計算方法の一例は、以下で説明する。
【００８７】
図１３は、実施例３における安定度計算部の処理の全体を示すフローチャートの例である。
【００８８】
ステップＳ５０１は、初期潮流状態を算出する。初期潮流計算は、ニュートンラフソン法を用いても良いし、直流法を用いても良い。
【００８９】
ステップＳ５０２は、想定する事故条件を入力する。例えば、２回線送電線＃１で三相地絡事故が発生し、事故から０．１秒後に送電線＃１の１回線を開放することによって事故を除去する。その後事故から０．２秒後に再閉路して系統を元の状態に戻す、といった事故を想定する。そのとき、事故除去・再閉路のタイムスケールや、系統状態の変化に伴う送電線の線路インピーダンス情報を入力する。
【００９０】
ステップＳ５０３は、時間領域シミュレーションの開始時刻を設定する。以下、ステップＳ５１５までの処理を繰り返すことで時間領域シミュレーションを実行する。
【００９１】
ステップＳ５０４は、発電機のグルーピンの修正が必要かどうか判定する。判定には、例えば、一つ前の時刻における発電機相差角の変化量の大きさを用いることが考えられる。一つ前の時刻における発電機相差角の変化量が一定値を超えた場合、動揺変化が大きな状態であるため、発電機毎の相差角の変化の傾向を判別しやすい。そのため、上記の場合にグルーピングを見直すため、ステップＳ５０５の処理に進んでグルーピング作成部１０２で処理を実施する。もしグルーピング修正が必要ない場合は、ステップＳ５０７に進む。
【００９２】
ステップＳ５０５は、グルーピング作成部１０２で発電機相差角の動揺変化が同じ傾向のものを同一グループとなるように、発電機のグループを作成する。例えば、発電機間の距離や発電機容量などを指標として、グルーピングを実施する。この指標となるデータは、データ入力部１０１で入力する。
【００９３】
ステップＳ５０６は、ステップＳ５０５で修正した発電機グループ毎に、代表の発電機を選定する。選定方法は、図６に示したように、データ入力部１０１で入力された発電機に関する情報を基に選定する。このとき、選定の指標には、発電機容量でも良いし、発電機の単位慣性定数でも良い。
【００９４】
ステップＳ５０７は、換算係数の修正が必要かどうか判定する。判定には、例えば、一つ前の時刻における発電機相差角の変化量の大きさを用いることが考えられる。一つ前の時刻における発電機相差角の変化量が一定値を超えた場合、動揺変化が大きな状態であるため、発電機毎の違いを区別しやすい。そのため、上記の場合に換算係数を見直すため、ステップＳ５０８の処理に進んで換算係数算出部３０１で処理を実施する。もし換算係数の修正が必要ない場合は、ステップＳ５０９に進む。
【００９５】
ステップＳ５０８は、換算係数を修正する。換算係数の算出方法は、例えば、換算係数算出部３０１において、図１１に示したような処理を実施して、計算する。
ステップＳ５０９は、系統状態の変更を判定する。ステップＳ５０２で入力した事故条件に基づいて、事故除去や再閉路に伴って系統状態が変化する場合、ステップＳ５１０に進んで系統状態を変更する。もし、次の時刻でも系統状態が変化しない場合には、ステップＳ５１１に進むように処理する。
【００９６】
ステップＳ５１０は、ステップＳ５０９で系統状態の変更の必要があると判定された場合に、系統状態を変更する。例えば、線路インピーダンスの値を変更する。
【００９７】
ステップＳ５１１は、系統計算を実行する。系統状態や発電機端子電圧の状態に応じて系統の線路潮流やノード電圧が変化する。そのため、潮流計算を実施して系統の状態を求める。このとき、潮流計算はニュートンラフソン法でも良いし、直流法でも良い。ただし、計算が収束しない場合には数値解法を変更して解を求める。
【００９８】
ステップＳ５１２は、積分計算を実施する。時刻が進んで系統状態が変化することで、発電機の動揺がどのように変化するのか、発電機の動揺方程式を解いて相差角の変化量を求める。このとき、数値解法はルンゲクッタ法でも良いし、台形法でも良い。
【００９９】
ステップＳ５１３は、ステップＳ５１２で求めた代表発電機の相差角変化量と換算係数を用いて、代表以外の発電機の動揺を簡易的に算出する。処理の詳細は、図８で説明したとおりである。
【０１００】
ステップＳ５１４は、時間刻みΔｔだけ時間を進める。
【０１０１】
ステップＳ５１５は、時間領域シミュレーションの終了時刻かどうか判定する。もしシミュレーションの終了時刻ではない場合は、ステップ５０４に戻って処理を繰り返す。もしシミュレーションの終了時刻であれば、処理を終了する。
【０１０２】
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。 
【０１０３】
また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、ＳＳＤ（Solid State Drive）等の記録装置、または、ＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に置くことができる。 
【０１０４】
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。
【符号の説明】
【０１０５】
１００  電力系統解析装置
１０１  データ入力部
１０２  グルーピング作成部
１０３  代表の発電機決定部
１０４  発電機情報記憶部
１０５  系統モデル作成部
１０６  安定度計算部
１０７  換算係数入力部
１０８  安定度評価部
１０９  出力部
２００  電力系統安定化制御装置
３０１  換算係数算出部
３０２  換算係数記憶部
【図１】
2017118614.tif 000003
【図２】
2017118614.tif 000004
【図３】
2017118614.tif 000005
【図４】
2017118614.tif 000006
【図５】
2017118614.tif 000007
【図６】
2017118614.tif 000008
【図７】
2017118614.tif 000009
【図８】
2017118614.tif 000010
【図９】
2017118614.tif 000011
【図１０】
2017118614.tif 000012
【図１１】
2017118614.tif 000013
【図１２】
2017118614.tif 000014
【図１３】
2017118614.tif 000015