(57)【要約】【課題】  Ｐ−Ｖカーブによって、電力系統におけるモータのストール現象を適切に説明できる電圧安定性評価方法を提供する。【解決手段】  電力系統を、１つのモータと１つの定インピーダンスが並列接続された負荷を有する一負荷系統に縮約して、電圧安定性を評価する電圧安定性評価方法であって、数値データの入力を受け付けるデータ入力過程と、負荷母線電圧を求める負荷母線電圧計算過程と、モータの電気入力トルク、モータの機械出力トルクおよび定インピーダンスの消費電力を求めるトルク・消費電力計算過程と、電力−電圧平面上において、モータの電気入力トルクと定インピーダンスの消費電力の和を示す曲線と、モータの機械出力トルクと定インピーダンスの消費電力の和を示す曲線が、不安定平衡点を有する場合に不安定、有しない場合に安定であると評価する安定性評価過程と、からなる。【選択図】      図５
2017051068.tif 000002
【特許請求の範囲】
【請求項１】
複数の負荷を有する電力系統を、１つの模擬的なモータと１つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価方法であって、
データ入力手段が前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力過程と、
負荷母線電圧計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入して負荷母線電圧を求める負荷母線電圧計算過程と、
トルク・消費電力計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入してモータの電気入力トルク、モータの機械出力トルクおよび定インピーダンスの消費電力を求めるトルク・消費電力計算過程と、
安定性評価手段が、電力−電圧平面上において、モータの電気入力トルクと定インピーダンスの消費電力の和を示す曲線と、モータの機械出力トルクと定インピーダンスの消費電力の和を示す曲線が、不安定平衡点を有する場合に前記一負荷系統を不安定であると評価し、不安定平衡点を有しない場合に安定であると評価する安定性評価過程と、
からなることを特徴とする電力系統の電圧安定性評価方法。
【請求項２】
複数の負荷を有する電力系統を、１つの模擬的なモータと１つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価プログラムであって、
前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力ステップ（Ｓ１）と、
記憶装置に保存された計算式を読み込み、前記数値データを代入して負荷母線電圧を求める負荷母線電圧計算ステップ（Ｓ２）と、
記憶装置に保存された計算式を読み込み、前記数値データを代入してモータの電気入力トルク、モータの機械出力トルクおよび定インピーダンスの消費電力を求めるトルク・消費電力計算ステップ（Ｓ３）と、
電力−電圧平面上において、モータの電気入力トルクと定インピーダンスの消費電力の和を示す曲線と、モータの機械出力トルクと定インピーダンスの消費電力の和を示す曲線が、不安定平衡点を有する場合に前記一負荷系統を不安定であると評価し、不安定平衡点を有しない場合に安定であると評価する安定性評価ステップ（Ｓ５）と、
をコンピュータに実行させることを特徴とする電力系統の電圧安定性評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統における、モータ（誘導電動機）を考慮した電圧安定性評価方法およびプログラムに関する。
【背景技術】
【０００２】
近年、電力系統の電圧安定性解析において、モータ負荷を考慮することの重要性が認識されるようになっており、そのために、モータのストール現象を説明できることが必要とされている。従来、特許文献１に示すように、モータのストール現象は、トルク−速度曲線（Ｔ−ωカーブ）によって適切に説明できることが知られている。すなわち、特許文献１の図２に示すＴ−ωカーブにおいて、Ａ点およびＣ点が安定平衡点、Ｂ点が不安定平衡点であるが、２つの安定平衡点のうち、正常時の動作点はＡ点であり、Ｃ点においては、回転速度および内部抵抗が異常に低いため、モータの有効電力消費が少なく、無効電力消費が非常に大きな状態で安定しており、モータがストールした状態である。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第４５８４１５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、一般に、電力系統の電圧安定性解析には、Ｐ−Ｖカーブが用いられているが、従来、Ｐ−Ｖカーブはモータのストール現象によって生じる電圧不安定現象の説明には向いていないとされていた。しかしながら、電圧安定性解析に広く用いられているＰ−Ｖカーブによって、モータのストール現象も説明することができれば、当業者にとってはより利便性が高い。
【０００５】
本発明は、上記事情を鑑みたものであり、Ｐ−Ｖカーブによって、電力系統におけるモータのストール現象を適切に説明できる電圧安定性評価方法およびプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明のうち請求項１の発明は、複数の負荷を有する電力系統を、１つの模擬的なモータと１つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価方法であって、データ入力手段が前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力過程と、負荷母線電圧計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入して負荷母線電圧を求める負荷母線電圧計算過程と、トルク・消費電力計算手段が記憶装置に保存された計算式を読み込み、前記数値データを代入してモータの電気入力トルク、モータの機械出力トルクおよび定インピーダンスの消費電力を求めるトルク・消費電力計算過程と、安定性評価手段が、電力−電圧平面上において、モータの電気入力トルクと定インピーダンスの消費電力の和を示す曲線と、モータの機械出力トルクと定インピーダンスの消費電力の和を示す曲線が、不安定平衡点を有する場合に前記一負荷系統を不安定であると評価し、不安定平衡点を有しない場合に安定であると評価する安定性評価過程と、からなることを特徴とする。
【０００７】
本発明のうち請求項２の発明は、複数の負荷を有する電力系統を、１つの模擬的なモータと１つの模擬的な定インピーダンスが並列接続された負荷を有する一負荷系統に縮約し、該一負荷系統の電圧安定性を縮約前の系統の電圧安定性とみなす電力系統の電圧安定性評価プログラムであって、前記一負荷系統を構成する要素の数値データの入力を受け付けるデータ入力ステップと、記憶装置に保存された計算式を読み込み、前記数値データを代入して負荷母線電圧を求める負荷母線電圧計算ステップと、記憶装置に保存された計算式を読み込み、前記数値データを代入してモータの電気入力トルク、モータの機械出力トルクおよび定インピーダンスの消費電力を求めるトルク・消費電力計算ステップと、電力−電圧平面上において、モータの電気入力トルクと定インピーダンスの消費電力の和を示す曲線と、モータの機械出力トルクと定インピーダンスの消費電力の和を示す曲線が、不安定平衡点を有する場合に前記一負荷系統を不安定であると評価し、不安定平衡点を有しない場合に安定であると評価する安定性評価ステップと、をコンピュータに実行させることを特徴とする。
【発明の効果】
【０００８】
本発明によれば、従来、電圧安定性解析に広く用いられている、電力−電圧平面上の電気入力トルクおよび機械出力トルクに相当する曲線（Ｐ−Ｖカーブ）によって、モータのストール現象も説明することができるので、当業者にとってはより利便性が高い。
【図面の簡単な説明】
【０００９】
【図１】本発明の電圧安定性評価に用いる系統モデルの模式図。
【図２】故障発生時の負荷母線電圧などの時間変化をシミュレーションした結果を示すグラフ。
【図３】Ｔ−ωカーブを示すグラフ。
【図４】モータのみについてのＰ−Ｖカーブを示すグラフ。
【図５】モータと定インピーダンス負荷を合わせたＰ−Ｖカーブを示すグラフ（不安定の場合）。
【図６】モータと定インピーダンス負荷を合わせたＰ−Ｖカーブを示すグラフ（安定の場合）。
【図７】電圧安定性評価プログラムの処理の流れを示すフローチャート。
【図８】電圧安定性評価装置の出力装置の表示例を示す図。
【発明を実施するための形態】
【００１０】
本発明は、複数の負荷を有する電力系統を、モータ１０１と定インピーダンス１０２が並列接続された１つの負荷を有する一負荷系統に縮約し、この一負荷系統についての電圧安定性を評価するものであって、図１に、その系統モデルを示す。無限大母線１０３から系統リアクタンス１０４を介して負荷母線１０５に送電されており、モータ１０１は一次抵抗と励磁アドミタンスを無視したすべりモデルであって、定リアクタンス１０１ａと可変コンダクダンス１０１ｂからなり、さらに調相用設備に相当するサセプタンス１０６を有する。
【００１１】
そしてまず、この系統モデルに基づいて、故障発生時の数値解析シミュレーションを行う。ここでは、系統リアクタンスの負荷側で三相地絡故障（故障継続時間１００ｍｓ）を与え、故障リアクタンスをＸｆ＝０，ｊ０．１とした２ケースでシミュレーションした。なお、系統定数およびパラメータは表１に示すとおりである。
【表１】
2017051068.tif 000003
【００１２】
このシミュレーションの結果を示したのが図２のグラフであり、故障リアクタンスがＸｆ＝ｊ０．１のケース（故障の規模が小さい場合）では、瞬低後、電圧（Ｖｒ）が回復しているが、Ｘｆ＝０のケース（故障の規模が大きい場合）では、瞬低後、電圧が回復しない。モータのストール現象によって、電圧不安定現象が生じているためである。
【００１３】
そして、本発明の電圧安定性評価方法によれば、上記のような数値解析シミュレーションによることなく、モータのストール現象を説明できる。以下において、その具体的な手順について説明する。まず、モータの電気入力トルクおよび機械出力トルクを求める式を導出する。図１において、ＶＳ、ＸＳ、Ｖｒ、ＰＺ、Ｐｍ、Ｘｍおよびモータ二次抵抗Ｒ２を与えると、以下の各値が求められる。
【００１４】
無限大母線と負荷母線の相差角δは、以下の式（１）で表される。
【数１】
2017051068.tif 000004
無限大母線における無効電力ＱＳは、以下の式（２）で表される。
【数２】
2017051068.tif 000005
負荷母線における無効電力Ｑｒは、以下の式（３）で表される。
【数３】
2017051068.tif 000006
モータによって消費される無効電力Ｑｍは、以下の式（４）で表される。
【数４】
2017051068.tif 000007
調相用設備から供給される無効電力ＱＺは、以下の式（５）で表される。
【数５】
2017051068.tif 000008
定インピーダンス負荷のコンダクタンスＧＺは、以下の式（６）で表される。
【数６】
2017051068.tif 000009
調相用設備のサセプタンスＢＺは、以下の式（７）で表される。
【数７】
2017051068.tif 000010
モータ消費電力初期値Ｐｍ０は、以下の式（８）で表される。ただし、ｓ０はすべりの初期値である（ｓ０＝１−ω０）。
【数８】
2017051068.tif 000011
さらに、式（８）より、式（９）が導かれる。ただし、２解のうち有効な方を解とする。
【数９】
2017051068.tif 000012
モータ内部抵抗初期値Ｒｍ０は、以下の式（１０）で表される。
【数１０】
2017051068.tif 000013
【００１５】
続いて、Ｐ−Ｖカーブを描くにあたり、負荷母線電圧Ｖｒ（ω）およびモータ消費電力Ｐｍ（ω）を導出する。（ω）はωの関数であることを示す。モータ内部抵抗Ｒｍ（ω）を式（１１）とおいて、モータのインピーダンスＺｍ（ω）［ベクトル量］は、式（１２）で表される。ただし、ｓはモータのすべりである。
【数１１】
2017051068.tif 000014
【数１２】
2017051068.tif 000015
モータ、定インピーダンス負荷、調相用設備込のアドミタンスＹｒ（ω）［ベクトル量］は、以下の式（１３）で表される。
【数１３】
2017051068.tif 000016
アドミタンスに系統リアクタンスを含んだインピーダンスＺＳ（ω）［ベクトル量］は、以下の式（１４）で表される。
【数１４】
2017051068.tif 000017
系統リアクタンスを流れる電流ＩＳ（ω）［ベクトル量］は、以下の式（１５）で表される。
【数１５】
2017051068.tif 000018
負荷母線電圧Ｖｒ（ω）［ベクトル量］およびその絶対値Ｖｒ（ω）は、以下の式（１６）および式（１７）で表される。
【数１６】
2017051068.tif 000019
【数１７】
2017051068.tif 000020
モータ消費電力Ｐｍ（ω）は、以下の式（１８）で表される。
【数１８】
2017051068.tif 000021
【００１６】
さらに、モータ消費電力Ｐｍ（ω）は、モータの電気入力トルクＴｅ（ω）に等しいことを示す。モータの二次入力電力Ｐ２（ω）は、二次電流をＩ２（ω）とすると、以下の式（１９）で表される。
【数１９】
2017051068.tif 000022
モータの二次軸出力電力Ｐ２Ｓ（ω）は、以下の式（２０）で表される。
【数２０】
2017051068.tif 000023
式（１９）および式（２０）より、電気入力トルクＴｅ（ω）、二次軸出力電力Ｐ２Ｓ（ω）および二次入力電力Ｐ２（ω）の関係は、以下の式（２１）で表される。ただし、ω＝１−ｓである。
【数２１】
2017051068.tif 000024
式（２１）より、電気入力トルクＴｅ（ω）は、二次入力電力Ｐ２（ω）に等しいことがわかる。さらに、一次抵抗と励磁リアクタンスを無視すると、モータ消費電力Ｐｍ（ω）は、二次入力電力Ｐ２（ω）に等しくなる。すなわち、以下の式（２２）が導出される。
【数２２】
2017051068.tif 000025
また、機械出力トルクＴｍ（ω）は、モータ回転速度初期値をω０として、以下の式（２３）で表される。
【数２３】
2017051068.tif 000026
【００１７】
以上の式（１）〜（１８）、（２２）、（２３）から、図３に示すように、トルク−速度曲線（Ｔ−ωカーブ）を描くことができる。このＴ−ωカーブによってモータのストール現象を説明できることは、前記のとおり、従来知られている。ここで、この図３のグラフ上において、負荷母線電圧Ｖｒ（ω）を併せて描く。ただし、見易さを考慮して、０．４倍して表示する。これによれば、回転速度ωが増大すると、負荷母線電圧Ｖｒ（ω）も増大することがわかる。したがって、回転速度ωの代わりに負荷母線電圧Ｖｒ（ω）を横軸にとって表現することも可能である。そして、縦軸（トルク）と横軸（電圧）を入れ替えれば、図４に示すように、Ｔ−ωカーブを等価変換したＰ−Ｖカーブが得られる。Ｔ−ωカーブはモータだけに着目したものであるから、図４のＰ−Ｖカーブもモータだけに着目したものである。
【００１８】
そこで、次にモータと定インピーダンス負荷を合わせたＰ−Ｖカーブを描く。定インピーダンス消費電力ＰＺ（ω）は、以下の式（２４）で表される。
【数２４】
2017051068.tif 000027
電気入力トルクＴｅ（ω）を表す式（１８）、（２２）と、機械出力トルクＴｍ（ω）を表す式（２３）に、式（２４）を加えて、それぞれＴｅ＋ＰＺ、Ｔｍ＋ＰＺとすることで、モータと定インピーダンス負荷を合わせたＰ−Ｖカーブを描くことができる。
【００１９】
図５および図６に示したのは、この手法により描いたＰ−Ｖカーブであり、図５は、系統定数およびパラメータを、前記の図２に示したシミュレーションと同じ条件にした場合、図６は、系統リアクタンスＸＳ＝０．５とした場合（その他の系統定数およびパラメータは図５の場合と同じ）である。図５においては、電気入力トルクを示す曲線（Ｔｅ＋ＰＺ）と機械出力トルクを示す曲線（Ｔｍ＋ＰＺ）が３点で交わっており、図６においては、電気入力トルクを示す曲線（Ｔｅ＋ＰＺ）と機械出力トルクを示す曲線（Ｔｍ＋ＰＺ）が１点で交わっている。両曲線の交点は、入力と出力が釣り合う平衡点であり、この平衡点に基づいて電力系統の電圧安定性を評価できる。
【００２０】
まず、図５に示すように、３個の平衡点を有する場合において、各平衡点を電圧が大きい側からＡ点、Ｂ点、Ｃ点とする。Ａ点とＢ点の間では、電気入力トルクＴｅ（＋ＰＺ）が機械出力トルクＴｍ（＋ＰＺ）を上回るので、モータは加速し、電圧は上昇する。一方、Ｂ点とＣ点の間では、電気入力トルクＴｅ（＋ＰＺ）が機械出力トルクＴｍ（＋ＰＺ）を下回るので、モータは減速し、電圧は低下する。よって、Ｂ点から電圧が少し上昇（すなわちモータが少し加速）すると、加速領域に入って電圧はさらに上昇し、Ｂ点から電圧が少し低下（すなわちモータが少し減速）すると、減速領域に入って電圧はさらに低下する。したがって、Ｂ点は不安定平衡点である。そして、Ａ点およびＣ点では逆の挙動になるので、Ａ点とＣ点は安定平衡点である。２つの安定平衡点のうち、Ａ点は、モータが定格回転速度付近で動作する正常動作点である。一方、Ｃ点は、回転速度および内部抵抗が異常に低いため、モータの有効電力消費が少なく、無効電力消費が非常に大きな状態で安定しており、モータがストールした状態である。
【００２１】
このような、電力−電圧平面上で電気入力トルクＴｅ（＋ＰＺ）と機械出力トルクＴｍ（＋ＰＺ）の交点が３つとなる（すなわち、２つの安定平衡点の間に１つの不安定平衡点を有する）電力系統において、正常動作中（Ａ点の状態）に故障が発生した場合、故障の規模が小さければ、動作点はＡ点に戻るが、故障の規模が大きいと、動作点はＡ点からＢ点を超えてＣ点に移ってしまい、Ａ点には戻らない。すなわち、故障により低下した系統の電圧が復帰せず、電圧不安定現象が発生することになる。よって、このように電力−電圧平面上の電気入力トルクと機械出力トルクを示す曲線（Ｐ−Ｖカーブ）が不安定平衡点を有する場合、変動によりその点を超えると元の動作点に戻らなくなるので、その電力系統は、モータのストール現象によって電圧不安定になる可能性があると評価される。
【００２２】
なお、この評価結果は、図２に示したシミュレーションの結果と一致する。すなわち、故障リアクタンスＸｆ＝ｊ０．１のケースが、故障の規模が小さい場合であり、このケースでは、瞬低後、電圧が回復している（動作点がＡ点に戻っている）。一方、Ｘｆ＝０のケースが、故障の規模が大きい場合であり、このケースでは、瞬低後、モータがストールしており、電圧が回復しない（動作点がＢ点を超えてＣ点に移っている）。
【００２３】
また、図６に示すように、１個の平衡点（Ａ点）を有する場合においては、その点が安定平衡点となり、不安定平衡点は存在しない。この場合、故障により電圧が変動しても、速やかに元の動作点に復帰する。よって、このように電力−電圧平面上の電気入力トルクと機械出力トルクを示す曲線（Ｐ−Ｖカーブ）が不安定平衡点を有しない場合、変動が生じても常に元の動作点に戻るので、その電力系統は安定であると評価される。
【００２４】
以上のように、本発明の電圧安定性評価方法によれば、電気入力トルクと機械出力トルクについて、Ｔ−ωカーブを等価変換したＰ−Ｖカーブを描き、それが不安定平衡点を有する場合、その電力系統は不安定であると評価され、不安定平衡点を有しない場合、その電力系統は安定であると評価される。そして、電圧安定性解析に広く用いられているＰ−Ｖカーブによって、モータのストール現象も説明することができるので、当業者にとってはより利便性が高い。
【００２５】
実際に本発明の電圧安定性評価方法を実施するには、以下に示す電圧安定性評価プログラムを利用する。本発明の電圧安定性評価プログラムは、キーボードやマウスなどからなる入力装置、ディスプレイなどからなる出力装置、本プログラムの命令を順番に実行するＣＰＵ、本プログラムならびに本プログラムの実行に必要なデータおよび計算結果などを保存する記憶装置を構成要素とする標準的なコンピュータにより実行される。本プログラムを実行することにより、コンピュータが電圧安定性評価装置として機能する。
【００２６】
本プログラムをコンピュータに実行させると、コンピュータが各種の手段（データ入力手段、負荷母線電圧計算手段、トルク・消費電力計算手段、グラフ表示手段、安定性評価手段）として機能する。そして、ＣＰＵからの指令によって、図７に示すように、データ入力手段がデータ入力ステップＳ１を実行し、負荷母線電圧計算手段が負荷母線電圧計算ステップＳ２を実行し、トルク・消費電力計算手段がトルク・消費電力計算ステップＳ３を実行し、グラフ表示手段がグラフ表示ステップＳ４を実行し、安定性評価手段が安定性評価ステップＳ５を実行する。
【００２７】
本プログラムを実行すると、まずデータ入力手段が機能して、データ入力ステップＳ１が実行される。データ入力ステップＳ１では、記憶装置に保存されているテーブル状の入出力フォームのデータを読み込み、図８に示すような入出力フォーム１１を出力装置に表示する。入出力フォーム１１は、データを入力するための入力部１２と、結果を出力するための出力部（定数表示部１３ａ、関数表示部１３ｂ、判定表示部１３ｃ）と、Ｔ−ωカーブおよびＰ−Ｖカーブを表示するグラフ表示部１４から構成される。
【００２８】
入力部１２は、ＶＳ：無限大母線電圧、ＸＳ：系統リアクタンス、Ｖｒ：負荷母線電圧、ＰＺ：定インピーダンス消費電力、Ｐｍ：モータ消費電力、Ｘｍ：モータリアクタンス、Ｒ２：モータ二次抵抗、からなる見出しが横一列に記載され、その下にそれぞれの値の入力欄が設けられている。
【００２９】
出力部のうち、定数表示部１３ａは、入力部１２の下側に設けられており、δ：無限大母線と負荷母線の相差角、ＱＳ：無限大母線における無効電力、Ｑｒ：負荷母線における無効電力、ＱＺ：調相用設備から供給される無効電力、Ｑｍ：モータによって消費される無効電力、ＧＺ：定インピーダンス負荷のコンダクタンス、ＢＺ：調相用設備のサセプタンス、ｓ０：モータのすべり初期値、Ｒｍ０：モータの内部抵抗初期値、からなる見出しが横一列に記載され、その下にそれぞれの値の出力欄が設けられている。また、関数表示部１３ｂは、ω：回転速度、Ｖｒ：負荷母線電圧、Ｔｅ：電気入力トルク、Ｔｍ：機械出力トルク、ＰＺ：定インピーダンス消費電力、Ｔｅ＋ＰＺ、Ｔｍ＋ＰＺ、Ｔｅ−Ｔｍ、その他、前記の計算過程で算出される値の見出しが横一列に記載され、ωの下には０＜ω＜１の値が所定の刻み（１の近傍は変化量が大きいので細かく刻んである）で記載されており、その他の見出しの下には、それぞれ同じ行のωの値に対応する出力欄が設けられている。さらに、判定表示部１３ｃは、入力部１２の右側に設けられており、電力系統の安定性評価の結果を「安定」または「不安定」と表示する（図８の例では「不安定」）。
【００３０】
グラフ表示部１４は、縦軸に電気入力トルクＴｅおよび機械出力トルクＴｍをとり、横軸に回転速度ωをとった、Ｔ−ωカーブを記載するグラフと、縦軸に負荷母線電圧Ｖｒをとり、横軸に負荷消費電力（Ｔｅ＋ＰＺ，Ｔｍ＋ＰＺ）をとったＰ−Ｖカーブを記載するグラフが設けられている。
【００３１】
使用者は、系統を構成する要素の数値データを、図８に示すように、各欄に入力装置から入力する。各入力データは記憶装置に保存されると共に、そのまま入力部１２に表示される。
【００３２】
そして、入力された前記数値データを記憶装置から読み込み、前記の各式に代入して、モータのすべり初期値ｓ０、モータの内部抵抗初期値Ｒｍ０、定インピーダンス負荷のコンダクタンスＧＺ、調相用設備のサセプタンスＢＺなどの各定数の値を計算する。各計算結果は記憶装置に保存されると共に、定数表示部１３ａの各欄に表示される。
【００３３】
次に、負荷母線電圧計算手段が機能して、負荷母線電圧計算ステップＳ２が実行される。負荷母線電圧計算ステップＳ２では、入力された前記数値データおよびωの値を記憶装置から読み込み、前記の各式に代入して各ωの値に対する負荷母線電圧Ｖｒの値を計算する。各計算結果は記憶装置に保存されると共に、関数表示部１３ｂのωの値に対応する欄に表示される。
【００３４】
次に、トルク・消費電力計算手段が機能して、トルク・消費電力計算ステップＳ３が実行される。トルク・消費電力計算ステップＳ３では、入力された前記数値データおよびωの値を記憶装置から読み込み、前記の各式に代入して各ωの値に対する電気入力トルクＴｅ、機械出力トルクＴｍおよび定インピーダンス消費電力ＰＺの値を計算する。各計算結果は記憶装置に保存されると共に、関数表示部１３ｂのωの値に対応する欄に表示される。
【００３５】
次に、グラフ表示手段が機能して、グラフ表示ステップＳ４が実行される。グラフ表示ステップＳ４では、各ωの値とそれに対応するＴｅ、Ｔｍの値を記憶装置から読み込み、グラフ表示部１４のトルク−速度平面上に各点をプロットしたものを表示すると共に、記憶装置に保存する。また、各ωの値に対応するＶｒ（ω）の値と、それに対応するＴｅ＋ＰＺ、Ｔｍ＋ＰＺの値を記憶装置から読み込み、グラフ表示部１４の電力−電圧平面上に各点をプロットしたものを表示すると共に、記憶装置に保存する。この際、グラフ補間手段が機能してグラフ補間ステップが実行され、各点を線形補間やスプライン補間などの手法により補間してその曲線を表示するように本プログラムを構成すると、使用者がグラフをより見やすくなり望ましい。
【００３６】
次に、安定性評価手段が機能して、安定性評価ステップＳ５が実行される。安定性評価ステップＳ５は、各ω（負荷母線電圧Ｖｒ）の値に対する電気入力トルクＴｅおよび機械出力トルクＴｍの値を記憶装置から読み込み、Ｔｅ−Ｔｍ（＝（Ｔｅ＋ＰＺ）−（Ｔｍ＋ＰＺ））の値を求めて、その正負を確認する。ここで、不安定平衡点とは、その点よりも高電圧の領域ではＴｅ＋ＰＺ＞Ｔｍ＋ＰＺとなり、その点よりも低電圧の領域ではＴｅ＋ＰＺ＜Ｔｍ＋ＰＺとなる交点であるから、高電圧側から走査してＴｅ−Ｔｍの値が正から負に転じる部分があれば、そこに不安定平衡点を有することになり、正から負に転じる部分がなければ、不安定平衡点を有しないことになる（Ｔｅ−Ｔｍの欄の隣にフラグの欄を設け、Ｔｅ−Ｔｍの値が正から負に転じる部分はフラグ欄を１とし、そうでない部分はフラグ欄を０として、フラグ欄の中に１があれば不安定平衡点を有し、１がなければ不安定平衡点を有しないと判断する）。そして、不安定平衡点を有する場合には、電力系統は不安定であると判定し、不安定平衡点を有しない場合には、電力系統は安定であると判定する。不安定か安定かの判定結果は、判定表示部１３ｃに表示されると共に、記憶装置に保存される。なお、使用者は、グラフ表示部１４に表示されたグラフにより、安定平衡点および不安定平衡点を視覚的に確認することができる。
【００３７】
なお、上記各ステップにおいて、負荷母線電圧計算ステップＳ２とトルク・消費電力計算ステップＳ３は、一部の順序が入れ替わったり、同時並行的に実行されたりしても構わない。また、グラフ表示ステップＳ４は、安定性評価ステップＳ５の後に実行されるものであっても構わない。さらに、安定性評価ステップＳ５において、不安定平衡点の条件を満たす交点が存在するか否かの判断方法については、上記の方法に限られず、どのようなものであってもよい。
【００３８】
また、本発明の電圧安定性評価プログラムは、専用のソフトウェアとして実行されるものであっても、汎用の表計算ソフトウェアなどの上で実行されるものであってもよい。さらに、実際の系統を構成する要素の数値データを測定する計器と接続されたコンピュータ上で作動するプログラムであって、測定された数値データを読み込み、常時系統の安定性を評価するものであってもよい。
【符号の説明】
【００３９】
Ｓ１  データ入力ステップ
Ｓ２  負荷母線電圧計算ステップ
Ｓ３  トルク・消費電力計算ステップ
Ｓ４  グラフ表示ステップ
Ｓ５  安定性評価ステップ
【図１】
2017051068.tif 000028
【図２】
2017051068.tif 000029
【図３】
2017051068.tif 000030
【図４】
2017051068.tif 000031
【図５】
2017051068.tif 000032
【図６】
2017051068.tif 000033
【図７】
2017051068.tif 000034
【図８】
2017051068.tif 000035