(57)【特許請求の範囲】
【請求項１】
複数の負荷を有する電力系統を、一つの模擬的な誘導電動機と一つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価方法であって、
データ入力手段が前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力過程と、
電気的入力計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の電気的入力を求める電気的入力計算過程と、
機械的出力計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の機械的出力を求める機械的出力計算過程と、
安定性評価手段が誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する安定性評価過程と、
結果出力手段が安定性評価の結果を出力する結果出力過程と、からなり、
前記安定性評価過程が、
非正常動作領域抽出手段がトルク−速度平面上において電気的入力の最大値を与える入力最大回転速度を抽出する非正常動作領域抽出過程と、
入出力値比較手段が０から入力最大回転速度までの速度領域において電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有するか否かを判断する入出力値比較過程と、
判定手段が電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有する場合に前記一負荷系統を不安定であると評価し、交点を有さない場合に安定であると評価する判定過程と、
からなることを特徴とする電力系統の電圧安定性評価方法。
【請求項２】
複数の負荷を有する電力系統を、一つの模擬的な誘導電動機と一つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価装置であって、
前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力手段と、
記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の電気的入力を求める電気的入力計算手段と、
記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の機械的出力を求める機械的出力計算手段と、
誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する安定性評価手段と、
安定性評価の結果を出力する結果出力手段と、を備え、
前記安定性評価手段が、
トルク−速度平面上において電気的入力の最大値を与える入力最大回転速度を抽出する非正常動作領域抽出手段と、
０から入力最大回転速度までの速度領域において電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有するか否かを判断する入出力値比較手段と、
電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有する場合に前記一負荷系統を不安定であると評価し、交点を有さない場合に安定であると評価する判定手段と、
からなることを特徴とする電力系統の電圧安定性評価装置。
【請求項３】
複数の負荷を有する電力系統を、一つの模擬的な誘導電動機と一つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価プログラムであって、
前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力ステップ（Ｓ１）と、
記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の電気的入力を求める電気的入力計算ステップ（Ｓ２）と、
記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の機械的出力を求める機械的出力計算ステップ（Ｓ３）と、
誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する安定性評価ステップ（Ｓ４）と、
安定性評価の結果を出力する結果出力ステップ（Ｓ５）と、をコンピュータに実行させるためのものであって、
前記安定性評価ステップ（Ｓ４）が、
トルク−速度平面上において電気的入力の最大値を与える入力最大回転速度を抽出する非正常動作領域抽出ステップ（Ｓ４ａ）と、
０から入力最大回転速度までの速度領域において電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有するか否かを判断する入出力値比較ステップ（Ｓ４ｂ）と、
電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有する場合に前記一負荷系統を不安定であると評価し、交点を有さない場合に安定であると評価する判定ステップ（Ｓ４ｃ）と、
からなることを特徴とする電力系統の電圧安定性評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統の電圧安定性を評価するのに好適な評価方法および装置ならびにプログラムに関する。
【背景技術】
【０００２】
近年、電力系統では無効電力の瞬動供給力である電源の遠隔化や長距離大電力送電に伴う系統無効電力損失の増加が進んでおり、系統運用に当たっては、電圧維持について従来にも増して高度な技術が必要となっている。このため、電力系統の電圧運用においては、系統電圧の安定性評価技術が重要となっている。
【０００３】
系統電圧の安定性を評価するためには、一般的にＰ−Ｖ特性が用いられる。例えば図９のような無限大母線一負荷系統を想定する。ここで、VGIは無限大母線電圧、VLIは負荷内部電圧、XGは無限大母線から変圧器２次までのリアクタンス、XLTは変圧器２次から負荷端子までのリアクタンス、XLZは定インピーダンス負荷の内部リアクタンス、RLZは定インピーダンス負荷の内部抵抗、Cは補償用コンデンサ、PAは負荷の有効電力をそれぞれ表す。負荷の有効電力PAと電圧VLIとの関係は、横軸をPA、縦軸をVLIとして図１０のようなＰ−Ｖカーブで表現される。Ｐ−Ｖカーブは、系統側から決まる受電端の有効電力と電圧の関係を示すものであり、系統の運転点は、負荷の有効電力と電圧の関係（負荷の有効電力の電圧特性）を示すカーブとＰ−Ｖカーブの交点となる。この交点が、図１０のＰ−Ｖカーブのノーズ（Ｄ点）より上側の範囲にある場合は電圧安定状態、下側の範囲にある場合は不安定状態であると判断される。
【０００４】
このＰ−Ｖカーブによる安定性評価方法では、系統負荷の中でも静的負荷のみを対象としており、変圧器タップや調相設備に起因する比較的時間レンジが長い（１０秒以上程度）電圧変化における電圧安定性の解析には有効であるが、実際の系統負荷には動的負荷（誘導電動機等）が含まれ、系統擾乱によって発生する数秒以内の電圧変化における動的負荷および系統電圧の安定性評価ができない。非特許文献１においても、Ｐ−Ｖカーブでは負荷の過渡的な動特性が引き起こす電圧不安定現象の取り扱いが難しいことが報告されている。
【非特許文献１】Sekine, Ohtsuki, 「Cascaded Voltage Collapse」, IEEE Transactions on Power Systems, Vol.5, no.1, p.250-256, Feb.1990
【０００５】
よって従来、負荷の過渡的な動特性を考慮して系統電圧の安定性を評価するためには、詳細な数値解析シミュレーションが必要であった。しかし、収束計算を行うため、電圧安定性の裕度を求める場合や系統状態が変更された場合は膨大な繰り返し計算が必要となり、簡単に評価することは困難であった。
【０００６】
そこで、シミュレーション手法を用いず、より簡単に負荷の過渡的な動特性を考慮した系統電圧の安定性を評価するためには、何らかの動的モデルを想定する必要がある。動的負荷の代表は誘導電動機であり、非特許文献２において、定リアクタンス背後の誘導電動機と定インピーダンスの並列モデルにより観測波形と良い一致が得られることが確認されている。さらに、非特許文献３において、動的負荷を考慮する場合の高精度な系統縮約法としてＹ結線法が提示されており、この手法によれば、負荷ブランチの導入により上位系統から配電用変電所もしくは特別高圧需要家に至る送電系統のインピーダンスを模擬することができる。
【非特許文献２】上田、駒見、「物理的イメージと観測データに基づく電力系統の動的負荷モデル」、平成17年電気学会Ｂ部門大会、2005年
【非特許文献３】山岸、上田、駒見、「動的負荷を含む電力系統の縮約方法の検討」、電気学会電力技術電力系統技術研究会資料、PE-05-41、PSE-05-48、2005年9月
【発明の開示】
【発明が解決しようとする課題】
【０００７】
本発明は、上記事情を鑑みたものであり、収束計算を必要とする数値解析シミュレーションを行うことなく、負荷の過渡的な動特性を考慮した系統電圧の安定性を評価する方法および装置ならびにプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明のうち請求項１の発明は、複数の負荷を有する電力系統を、一つの模擬的な誘導電動機と一つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価方法であって、データ入力手段が前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力過程と、電気的入力計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の電気的入力を求める電気的入力計算過程と、機械的出力計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の機械的出力を求める機械的出力計算過程と、安定性評価手段が誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する安定性評価過程と、結果出力手段が安定性評価の結果を出力する結果出力過程と、からなる。
【０００９】
そしてさらに、前記安定性評価過程が、非正常動作領域抽出手段がトルク−速度平面上において電気的入力の最大値を与える入力最大回転速度を抽出する非正常動作領域抽出過程と、入出力値比較手段が０から入力最大回転速度までの速度領域において電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有するか否かを判断する入出力値比較過程と、判定手段が電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有する場合に前記一負荷系統を不安定であると評価し、交点を有さない場合に安定であると評価する判定過程と、からなることを特徴とする。
【００１０】
また、請求項２の発明は、複数の負荷を有する電力系統を、一つの模擬的な誘導電動機と一つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価装置であって、前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力手段と、記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の電気的入力を求める電気的入力計算手段と、記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の機械的出力を求める機械的出力計算手段と、誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する安定性評価手段と、安定性評価の結果を出力する結果出力手段と、を備える。
【００１１】
そしてさらに、前記安定性評価手段が、トルク−速度平面上において電気的入力の最大値を与える入力最大回転速度を抽出する非正常動作領域抽出手段と、０から入力最大回転速度までの速度領域において電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有するか否かを判断する入出力値比較手段と、電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有する場合に前記一負荷系統を不安定であると評価し、交点を有さない場合に安定であると評価する判定手段と、からなることを特徴とする。
【００１２】
また、請求項３の発明は、複数の負荷を有する電力系統を、一つの模擬的な誘導電動機と一つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価プログラムであって、前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力ステップと、記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の電気的入力を求める電気的入力計算ステップと、記憶装置に保存された計算式を読み込み、前記数値データを代入して誘導電動機の機械的出力を求める機械的出力計算ステップと、誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する安定性評価ステップと、安定性評価の結果を出力する結果出力ステップと、をコンピュータに実行させる。
【００１３】
そしてさらに、前記安定性評価ステップが、トルク−速度平面上において電気的入力の最大値を与える入力最大回転速度を抽出する非正常動作領域抽出ステップと、０から入力最大回転速度までの速度領域において電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有するか否かを判断する入出力値比較ステップと、電気的入力と機械的出力のそれぞれのトルク−速度曲線が交点を有する場合に前記一負荷系統を不安定であると評価し、交点を有さない場合に安定であると評価する判定ステップと、からなることを特徴とする。
【発明の効果】
【００１４】
本発明によれば、動特性を有する負荷の安定性に着目し系統電圧の安定性を評価するため、詳細な数値解析シミュレーションによる収束計算を行うことなく、実系統で発生しうる電圧不安定現象をより正確に把握することが可能となる。さらに、電圧不安定に至る領域を予め定量的に捉えておくことが可能となり、系統運用におけるオンライン監視や予防制御、また系統運用における調相設備の検討等に有効となる。
【発明を実施するための最良の形態】
【００１５】
本発明は、複数の負荷を有する電力系統を、誘導電動機と定インピーダンスの並列からなる一負荷系統に縮約し、該一負荷系統についての電圧安定性を評価するものである。そこでまず、縮約された一負荷系統について説明する。負荷を誘導電動機と定インピーダンスの並列とみなしたとき、無限大母線一負荷系統は図１のように表される。ここで、VGIは無限大母線電圧、VCは変圧器２次側電圧、VLTは負荷端子電圧、XGは無限大母線から変圧器２次までのリアクタンス、XLTは変圧器２次から負荷端子までのリアクタンス、XLZは定インピーダンス負荷の内部リアクタンス、XLMは誘導電動機負荷の内部リアクタンス、RLZは定インピーダンス負荷の内部抵抗、RLMは誘導電動機負荷の内部抵抗、NTは変圧器のタップ比、Cは補償用コンデンサをそれぞれ表す。このように表された無限大母線一負荷系統は、縮約前の系統の特性を良く表しており、この無限大母線一負荷系統の電圧安定性を評価することにより、縮約前の系統の安定性を評価したものとみなすことができる。
【００１６】
次に、本発明の安定性評価方法の具体的な手順について説明する。第一に、誘導電動機の電気的入力を求める。図１の系統において、系統擾乱により短時間の間に大きく変化しうるのは誘導電動機の内部抵抗RLMだけで、他は一定とみなすことができる。RLMは、定格回転速度をω0、その時の内部抵抗をRLM0とすると、回転速度ωによって式（１）のように表される。
【数１】
0004584157.tif 000002
【００１７】
ここで、誘導電動機に供給される電気的入力PLMについて考える。系統を構成する要素の数値データである諸定数（VGI，XG，XLT，XLZ，XLM，RLZ，NT，C，RLM0）が与えられると、以下の式（２）から式（１１）が展開され、最終的に式（１２）が求まる。ただし、jは虚数単位を表す。
【００１８】
誘導電動機のインピーダンスZLMは、以下の式（２）で表される。
【数２】
0004584157.tif 000003
定インピーダンス負荷のインピーダンスZLZは、以下の式（３）で表される。
【数３】
0004584157.tif 000004
ZLM、ZLZ、XLTを合わせたインピーダンスZLは、以下の式（４）で表される。
【数４】
0004584157.tif 000005
変圧器２次側の系統のインピーダンスZLCは、以下の式（５）で表される。
【数５】
0004584157.tif 000006
系統全体のインピーダンスZは、以下の式（６）で表される。
【数６】
0004584157.tif 000007
無限大母線電流IGは、以下の式（７）で表される。
【数７】
0004584157.tif 000008
変圧器２次側電圧VCは、以下の式（８）で表される。
【数８】
0004584157.tif 000009
変圧器２次側電流ILは、以下の式（９）で表される。
【数９】
0004584157.tif 000010
負荷端子電圧VLTは、以下の式（１０）で表される。
【数１０】
0004584157.tif 000011
誘導電動機を流れる電流ILMは、以下の式（１１）で表される。
【数１１】
0004584157.tif 000012
誘導電動機に供給される電気的入力PLMは、以下の式（１２）で表される。
【数１２】
0004584157.tif 000013
式（１２）は誘導電動機の回転速度ωを含む式であるから、これより誘導電動機のトルク−速度曲線を描くことができる。このトルク−速度曲線は、系統および定インピーダンス負荷も考慮した等価電源からこの誘導電動機に供給される電気的入力ＰLMを意味する。
【００１９】
第二に、誘導電動機の機械的出力を求める。誘導電動機の機械的負荷TMは、機械の種類により異なるが、一般的に回転速度ωのn乗に比例し、式（１３）のように表される。ただし、PLM0は定格消費電力を表す。
【数１３】
0004584157.tif 000014
【００２０】
ここで、前述の式（１）から式（１３）によれば、定めるべき諸定数は、VGI，XG，XLT，XLZ，XLM，RLZ，NT，C，RLM0，ω0，nであるが、実際の電力系統に即した以下の考え方を導入し、これらの諸定数を簡略化する。
【００２１】
一般に負荷の内部抵抗と内部リアクタンスの比は一定値であり、ここではXLM/RLM=XLZ/RLZとする。また、負荷における誘導電動機の比率をmとすれば、ZLMおよびZLZは以下の式（１４）、式（１５）のように表される。ただし、RLおよびXLIはそれぞれ負荷全体の抵抗およびリアクタンスを表す。
【数１４】
0004584157.tif 000015
【数１５】
0004584157.tif 000016
RLは、皮相電力がVC2/(RL＋j(XLT＋XLI))で表されることから、その実部をとって、以下の式（１６）のように表される。ただし、PLは負荷の消費電力を表す。
【数１６】
0004584157.tif 000017
また、一般に補償用コンデンサCの値は変圧器２次側電圧VCの目標値を定めて設定されるため、ここではVCを定めてCの値を逆算することとする。Cは、VC点での上位系から流れ込む無効電力QG、下位系に流れ込む無効電力QLにより、以下の式（１７）で表される。
【数１７】
0004584157.tif 000018
QGは以下の式（１８）で表される。ただし、Iｍは複素数の虚部を、＊は複素共役を表す。
【数１８】
0004584157.tif 000019
式（１８）において、IGは以下の式（１９）から求められる。
【数１９】
0004584157.tif 000020
式（１９）において、θは以下の式（２０）から求められる。
【数２０】
0004584157.tif 000021
また、QLは以下の式（２１）で表される。
【数２１】
0004584157.tif 000022
式（２１）において、ILは以下の式（２２）で表される。
【数２２】
0004584157.tif 000023
さらに、PLM0は以下の式（２３）で表される。
【数２３】
0004584157.tif 000024
以上により、定めるべき諸定数は、VGI，XG，XLT，XLI，NT，VC，PL，m，ω0，nとなる。このようにすれば、誘導電動機と定インピーダンスのそれぞれの抵抗とリアクタンスを個別に導出する必要がなくなり、より簡単に安定性評価を行うことができる。
【００２２】
第三に、誘導電動機の電気的入力および機械的出力をトルク−速度平面上で比較することにより電圧安定性を評価する。電気的入力PLM（以下Ｐと表す）および機械的出力TM（以下Ｔと表す）のトルク−速度曲線は諸定数の値により変化するが、両曲線の関係という点では、代表的な例として、図２に示すように両曲線が三点で交わる場合と、図３に示すように両曲線が一点で交わる場合の二つの場合が挙げられる。なお、系統の諸定数は、VGI=1.1，XG=0.4，XLT=0.2，XLI=0.1，NT=1.0，VC=1.0，PL=1.0，ω0=0.98，n=2とし、図２においては誘導電動機比率を５０％、図３においては誘導電動機比率を４０％とした。
【００２３】
電気的入力Ｐおよび機械的出力Ｔを示すトルク−速度曲線が図２のようになる場合、二つの曲線の交点、すなわち、入出力が釣り合う平衡点は三個ある（回転速度の大きい側からそれぞれＡ点、Ｂ点、Ｃ点とする）。このうちＡ点およびＣ点の近傍では、回転速度が僅かに上がるとＴ＞Ｐとなるので減速し、回転速度が僅かに下がるとＴ＜Ｐとなるので加速する。従ってＡ点およびＣ点は安定平衡点である。一方、Ｂ点においては逆の挙動となるので、Ｂ点は不安定平衡点である。二個の安定平衡点のうち正常時の動作点は、電気的入力Ｐが最大トルクとなる回転速度（入力最大回転速度ωＭ）よりも速い回転速度領域（以下、正常動作領域とする）のＡ点である。Ｃ点は回転速度および内部抵抗が異常に低いため誘導電動機の有効電力消費が少なく、無効電力消費が非常に大きい状態で安定している。この状態はストール（失速）と呼ばれる状態である。
【００２４】
このような電気的入力Ｐと機械的出力Ｔの交点が三つとなる電力系統において系統擾乱が発生すると、正常時の動作点であるＡ点から系統擾乱のショックにより不安定平衡点であるＢ点を超えてＣ点に移ってしまう。誘導電動機には同期化力がないので電力動揺もなく、Ａ→Ｂ→Ｃの移動は非振動的であり、反動でＡ点に復帰することはない。そして、動的負荷が不安定性となると、負荷の内部抵抗は０に近づくため負荷近傍の系統電圧も０に近づく。すなわち、系統擾乱のショックにより低下した系統の電圧が復帰せず、負荷の不安定現象の結果として電圧不安定現象が発生することになる。
【００２５】
一方、電気的入力Ｐおよび機械的出力Ｔを示すトルク−速度曲線が図３のようになる場合、両曲線の交点は一つであり、図２の場合におけるＡ点にあたる正常動作領域の安定平衡点しか存在しないため、系統擾乱で異常に回転速度が低下しても速やかに元の正常な運転点に復帰できる。すなわち、系統擾乱のショックにより系統の電圧が低下しても、元の電圧に戻ることになる。
【００２６】
このように、誘導電動機の電気的入力Ｐと機械的出力Ｔのそれぞれのトルク−速度曲線を比較して、安定平衡点が正常動作領域のみにある場合、誘導電動機は安定運転状態であり、電力系統の電圧は安定であると評価できる。一方、安定平衡点が正常動作領域以外にある場合、誘導電動機は不安定運転状態であり、電力系統の電圧は不安定であると評価できる。
【００２７】
さらに、安定平衡点が正常動作領域のみにあり電力系統の電圧が安定であると判断される場合において、電気的入力Ｐと機械的出力Ｔの距離を把握することで電圧安定性の余裕を把握することができる。すなわち、電気的入力Ｐと機械的出力Ｔのそれぞれのトルク−速度曲線について、正常動作領域以外の領域での両曲線同士の距離が離れていればいるほど、電圧安定性に余裕があるといえる。
【００２８】
この評価結果を数値解析シミュレーションによる安定性評価と比較する。図１の系統について、同等の条件で誘導電動機比率が５０％の場合と４０％の場合を想定し、系統擾乱が発生したときの負荷端電圧の時間変化を数値解析シミュレーションにより求めたものが図４に示すグラフである。誘導電動機比率が５０％の場合、系統擾乱により変動した電圧が復帰せず、系統電圧が不安定状態になるのに対し、誘導電動機比率が４０％の場合、変動した電圧が元の電圧に復帰し、系統電圧は安定状態となることがわかる。この結果は、本発明の安定性評価方法による結果と一致している。
【００２９】
一方、同等の条件で従来の静的負荷モデルとＰ−Ｖカーブによる安定性評価を行うと、誘導電動機比率が５０％の場合のＰ−Ｖカーブは図５のようになる。なお、負荷の有効電力の電圧特性は、５０％の誘導電動機が定電力負荷であり、残りの５０％の定インピーダンス負荷が電圧の２乗に比例することから、静的負荷全体としては電圧の１乗に比例するものと設定すれば良い。このときのＰ−Ｖカーブと負荷特性カーブとの交点は、図５に示すように、Ｐ−Ｖカーブのノーズの上側にあるため、系統は安定と判断され、これは数値解析シミュレーションおよび本発明の方法による評価結果と異なる。
【００３０】
このように、従来の静的負荷モデルとＰ−Ｖカーブによる安定性評価方法は、系統の特性が十分に反映されず、本来不安定である系統を安定であると評価してしまう場合もある。それに対し、本発明の動的負荷モデルとトルク−速度曲線による安定性評価方法は、系統の特性を十分に反映し、数値解析シミュレーションと同等の結果を得ることができる。
【００３１】
実際に本発明の電力系統における安定性評価方法を実施するには、以下に示す安定性評価装置ならびにプログラムを利用する。
【００３２】
本発明の安定性評価プログラムは、キーボードやマウス等からなる入力装置、ディスプレイ等からなる出力装置、本プログラムの命令を順番に実行するＣＰＵ、本プログラムおよび本プログラムの実行に必要なデータおよび計算結果等を保存しておく記憶装置を構成要素とする標準的なコンピュータにより実行される。また、本発明の安定性評価装置は、前記プログラムが実行されるコンピュータにより構成される（図６）。
【００３３】
本プログラムをコンピュータに実行させた場合、コンピュータが各種の手段（データ入力手段、電気的入力計算手段、機械的出力計算手段、安定性評価手段、結果出力手段）として機能し、ＣＰＵからの指令によって図７に示すように、データ入力手段がデータ入力ステップＳ１を実行し、電気的入力計算手段が電気的入力計算ステップＳ２を実行し、機械的出力計算手段が機械的出力計算ステップＳ３を実行し、安定性評価手段が安定性評価ステップＳ４を実行し、結果出力手段が結果出力ステップＳ５を実行することで、電力系統の電圧安定性を評価する。なお、安定性評価手段は、さらに非正常動作領域抽出手段、入出力値比較手段および判定手段に細分され、それぞれ非正常動作領域抽出ステップＳ４ａ、入出力値比較ステップＳ４ｂおよび判定ステップＳ４ｃを実行する。
【００３４】
本プログラムを実行すると、まずデータ入力手段が機能して、データ入力ステップＳ１が実行される。データ入力ステップＳ１では、記憶装置に保存されているテーブル状の入出力フォーム１１を読み込んで、図８に示すような入出力フォーム１１を出力装置に表示する。入出力フォーム１１は、データを入力するための入力部１２と、結果を出力するための出力部１３と、トルク−速度曲線を表示するグラフ表示部１４と、から構成される。
【００３５】
入力部１２は、VGI：無限大母線電圧、XG：無限大母線から変圧器２次までのリアクタンス、XLT：変圧器２次から負荷端子までのリアクタンス、XLI：負荷リアクタンス、NT：変圧器のタップ比、VC：変圧器２次側電圧、PL：負荷の消費電力、m：誘導機比率、ω0：定格回転速度、n：機械的出力乗数、からなる見出しが横一列に記載され、その下にそれぞれの値の入力欄が設けられている。
【００３６】
出力部１３は、ω：回転速度、Ｐ：電気的入力、Ｔ：機械的出力、からなる見出しが横一列に記載され、ωの下には０から１までの値が任意の刻みで記載されており（図８の例では０．０１刻み）、ＰおよびＴの下には、それぞれ同じ行のωの値に対応する出力欄が設けられている。また、系統の安定性評価の結果を「安定」又は「不安定」と表示する評価欄（図８の例では「不安定」）が設けられている。
【００３７】
グラフ表示部１４は、縦軸に電気的入力Ｐおよび機械的出力Ｔをとり、横軸に回転速度ωをとったトルク−速度平面が設けられている。
【００３８】
使用者は、系統を構成する要素の数値データを図８に示すように各欄に入力装置から入力する。各入力データは記憶装置に保存されると共に、そのまま入力部１２に表示される。
【００３９】
次に、電気的入力計算手段が機能して、電気的入力計算ステップＳ２が実行される。電気的入力計算ステップＳ２では、入力された前記数値データおよびωの値を記憶装置から読み込み、式（１）、式（４）から式（１２）および式（１４）から式（２２）に代入して各ωの値に対する電気的入力Ｐの値を計算する。ただし、ω＝１に対しては式（１）が定義できず電気的入力Ｐを計算することができない。そこで、ω＝１に対する電気的入力Ｐの値については、当該欄に予め０を表示しておくと共に、記憶装置に保存しておくか、便宜的に、ω＝１の代わりに１に限りなく近い数（例えば、ω＝０．９９９９９）を用いれば良い。各計算結果は記憶装置に保存されると共に、出力部１３のωの値に対応する欄に表示される。
【００４０】
次に、機械的出力計算手段が機能して、機械的出力計算ステップＳ３が実行される。機械的出力計算ステップＳ３では、入力された前記数値データおよびωの値を記憶装置から読み込み、式（１３）および式（２３）に代入して、各ωの値に対する機械的出力Ｔの値を計算する。各計算結果は記憶装置に保存されると共に、出力部１３のωの値に対応する欄に表示される。
【００４１】
次に、安定性評価手段のうち非正常動作領域抽出手段が機能して、非正常動作領域抽出ステップＳ４ａが実行される。非正常動作領域抽出ステップＳ４ａは、各ωに対する電気的入力Ｐの値を記憶装置から読み込んで、Ｐが最大値となるときのωの値を抽出し、入力最大回転速度ωＭとする。ωＭの値は記憶装置に保存される。
【００４２】
次に、安定性評価手段のうち入出力値比較手段が機能して、入出力値比較ステップＳ４ｂが実行される。入出力値比較ステップＳ４ｂは、０からωＭまでの各ωに対する電気的入力Ｐおよび機械的出力Ｔの値を記憶装置から読み込んで、各ωについてＰとＴの値を比較する。ω＝０においてはＴ＝０であり必ずＴ＜Ｐであるから、Ｔ＞Ｐとなる点があるならば、両トルク−速度曲線は交わると判断される。交点を有するか否かの判断結果は記憶装置に保存される。
【００４３】
次に、安定性評価手段のうち判定手段が機能して、判定ステップＳ４ｃが実行される。判定ステップＳ４ｃは、電気的入力Ｐと機械的出力Ｔの両トルク−速度曲線が０からωＭまでの範囲に交点を有するか否かの判断結果を記憶装置から読み込み、交点を有する場合には、電力系統は不安定であると判定し、交点がない場合には、電力系統は安定であると判定する。不安定か安定かの判定結果は記憶装置に保存される。
【００４４】
次に、結果出力手段が機能して、結果出力ステップＳ５が実行される。結果出力ステップＳ５では、安定性評価の判定結果を記憶装置から読み込み、その結果に応じて、出力部１３の評価欄に「不安定」又は「安定」と表示する。さらに、各ωの値とそれに対応する電気的入力Ｐおよび機械的出力Ｔの値を記憶装置から読み込み、グラフ表示部１４のトルク−速度平面上に各点をプロットしたものを表示すると共に、記憶装置に保存する。この際、グラフ補間手段が機能してグラフ補間ステップが実行され、各点を線形補間やスプライン補間等の手法により補間してその曲線を表示するように本装置およびプログラムを構成すると、使用者がグラフをより見やすくなり望ましい。このグラフを表示することにより、安定平衡点を視覚的に確認することができ、また、正常動作領域以外の領域での両曲線同士の距離から、電圧安定性の余裕を確認することができる。
【００４５】
なお、上記各ステップにおいて、電気的入力計算ステップＳ２と、機械的出力計算ステップＳ３とは順序が入れ替わっても構わない。また、結果出力ステップＳ５のグラフ表示に関しては、安定性評価ステップＳ４の直前に実行されるものであっても構わない。
【００４６】
また、本発明の安定性評価プログラムは、専用のソフトウェアとして実行されるものであっても、汎用の表計算ソフトウェア等の上で実行されるものであっても良い。さらに、実際の系統を構成する要素の数値データを測定する計器と接続されたコンピュータ上で作動するプログラムであって、測定された数値データを読み込み、常時系統の安定性を評価するものであっても良い。
【図面の簡単な説明】
【００４７】
【図１】本発明の電圧安定性評価方法に用いられる動的負荷を考慮した無限大母線一負荷系統を示す図。
【図２】図１の系統において誘導電動機比率が５０％である場合のトルク−速度曲線を示すグラフ。
【図３】図１の系統において誘導電動機比率が４０％である場合のトルク−速度曲線を示すグラフ。
【図４】図１の系統において数値解析シミュレーションにより求めた電圧の時間変化を示すグラフ。
【図５】図１の系統において誘導電動機比率が５０％である場合に静的負荷のみを考慮して求められたＰ−Ｖカーブを示すグラフ。
【図６】電圧安定性評価装置の構成を示すブロック図。
【図７】電圧安定性評価プログラムの処理の流れを示すフローチャート。
【図８】電圧安定性評価装置の出力装置に表示される入出力フォームを示す図。
【図９】従来の電圧安定性評価方法に用いられた動的負荷を考慮しない無限大母線一負荷系統を示す図。
【図１０】従来の電圧安定性評価方法に用いられたＰ−Ｖカーブを示すグラフ。
【符号の説明】
【００４８】
Ｓ１  データ入力ステップ
Ｓ２  電気的入力計算ステップ
Ｓ３  機械的出力計算ステップ
Ｓ４  安定性評価ステップ
Ｓ４ａ  非正常動作領域抽出ステップ
Ｓ４ｂ  入出力値比較ステップ
Ｓ４ｃ  判定ステップ
Ｓ５  結果出力ステップ
【図１】
0004584157.tif 000025
【図２】
0004584157.tif 000026
【図３】
0004584157.tif 000027
【図４】
0004584157.tif 000028
【図５】
0004584157.tif 000029
【図６】
0004584157.tif 000030
【図７】
0004584157.tif 000031
【図８】
0004584157.tif 000032
【図９】
0004584157.tif 000033
【図１０】
0004584157.tif 000034