(57)【特許請求の範囲】
【請求項１】
電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、
前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算する複数の電力ミスマッチ計算部と、
前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から１つを最良の選択解として選択して出力する選択部と、
を有し、
前記複数の電力ミスマッチ計算部は、前記複数の解に対して、各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出し、
前記選択部は、前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの最大値のうちの最小値に対応する解を前記選択解として選択して出力する
ことを特徴とする潮流計算装置。
【請求項２】
前記複数の解と、前記複数の解に対する前記電力ミスマッチの最大値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする請求項１に記載の潮流計算装置。
【請求項３】
電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、
前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算する複数の電力ミスマッチ計算部と、
前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から１つを最良の選択解として選択して出力する選択部と、
を有し、
前記複数の電力ミスマッチ計算部は、前記複数の解に対して、各ノードにおける電力ミスマッチの平均値をそれぞれ算出し、
前記選択部は、前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの平均値のうちの最小値に対応する解を前記選択解として選択して出力することを特徴とする潮流計算装置。
【請求項４】
前記複数の解と、前記複数の解に対する前記電力ミスマッチの平均値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする請求項３に記載の潮流計算装置。
【請求項５】
電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、
前記複数の解に対して、各ノードにおける電力ミスマッチを計算し、計算された各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出する複数の電力ミスマッチ計算部と、
を有し、
前記複数の解と、前記複数の解に対する前記電力ミスマッチの最大値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする潮流計算装置。
【請求項６】
電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、
前記複数の解に対して、各ノードにおける電力ミスマッチを計算し、計算された各ノードにおける電力ミスマッチの平均値をそれぞれ算出する複数の電力ミスマッチ計算部と、
を有し、
前記複数の解と、前記複数の解に対する前記電力ミスマッチの平均値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする潮流計算装置。
【請求項７】
電力系統の状態を示す系統情報を格納する記憶装置を備えるコンピュータに、
前記系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める潮流計算処理と、
前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算すると共に、計算された各ノードにおける前記電力ミスマッチのうちの最大値をそれぞれ抽出する電力ミスマッチ計算処理と、
前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの最大値のうちの最小値に対応する解を最良の選択解として選択して出力する選択処理と、
を実行させることを特徴とする潮流計算プログラム。
【請求項８】
電力系統の状態を示す系統情報を格納する記憶装置を備えるコンピュータに、
前記系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める潮流計算処理と、
前記複数の解に対して、各ノードにおける電力ミスマッチを計算し、計算された各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出する複数の電力ミスマッチ計算処理と、
前記複数の解と、前記複数の解に対する前記電力ミスマッチの最大値とを対応付けて表示装置に表示させ、電力系統運用者に提示する出力処理と、
を実行させることを特徴とする潮流計算プログラム。
【請求項９】
電力系統の状態を示す系統情報を格納する潮流計算装置による潮流計算方法であって、
前記潮流計算装置が、前記系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求め、
前記潮流計算装置が、前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算すると共に、計算された各ノードにおける前記電力ミスマッチのうちの最大値をそれぞれ抽出し、
前記潮流計算装置が、前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの最大値のうちの最小値に対応する解を最良の選択解として選択して出力する
ことを特徴とする潮流計算方法。
【請求項１０】
電力系統の状態を示す系統情報を格納する潮流計算装置による潮流計算方法であって、
前記潮流計算装置が、前記系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求め、
前記潮流計算装置が、前記複数の解に対して、各ノードにおける電力ミスマッチを計算し、計算された各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出し、
前記潮流計算装置が、前記複数の解と、前記複数の解に対する前記電力ミスマッチの最大値とを対応付けて表示装置に表示させ、電力系統運用者に提示する
ことを特徴とする潮流計算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、潮流計算装置、潮流計算プログラム、および潮流計算方法に関する。
【背景技術】
【０００２】
電力系統では、発電所で発電された電力が様々な負荷地点に流れている。この「電力の流れ」を「潮流」と呼び、電力系統の状態を示す系統情報（系統データ）に基づいて、各送電線や変圧器の潮流や各ノード（母線）の電圧および位相角を計算するのが潮流計算である。
【０００３】
潮流計算の際に入力される系統情報には、例えば系統構成、設備定数、および運用条件に関する情報が含まれる。系統構成情報としては、発電機、変圧器、および送電線負荷の接続状態や、系統図などが入力される。また、設備定数情報としては、各設備・機器のインピーダンスやアドミタンス、変圧器タップ情報などが入力される。また、運用条件情報としては、発電機や負荷の（有効）電力および無効電力や、基準電圧などが入力される。一方、潮流計算の結果として出力される解には、少なくとも各ノードにおける電圧および位相角が含まれ、また、それらと系統情報とから求まる各送電線や変圧器の（有効）電力および無効電力潮流が含まれ得る。
【０００４】
潮流計算は、電力方程式を解くことによって解が求まる。電力方程式は以下の式（１）で示される。  
0005830370.tif 000002
ここで、Ｉｋはノードｋへ流入する電流、Ｙｋｊはノードｋからｊに対するアドミタンス（Ｙ：アドミタンス行列）、Ｖｊはノードｊに対する電圧、Ｓｋはノードｋの電力、＊は共役複素数を示す。なお、潮流計算には、各ノードに対して事前に指定する条件によって、ＰＱ指定とＰＶ指定とがあり、ＰＱ指定またはＰＶ指定により系統の電圧や潮流が決定される。ＰＱ指定の場合、式（１）において、各ノードｋに対して、有効電力Ｐｋと無効電力Ｑｋとを事前に指定し、電圧値｜Ｖｋ｜と位相角δｋとを計算する。一方、ＰＶ指定の場合、各ノードｋに対して、有効電力Ｐｋと電圧値｜Ｖｋ｜とを事前に指定し、無効電力Ｑｋと位相角δｋとを計算する。
【０００５】
ところで、潮流計算手法としては、いくつかの手法が知られている。例えば、特許文献１では、ニュートン・ラフソン（Newton-Raphson）法による潮流計算が開示されている。また、例えば、特許文献２では、直流法による潮流計算が開示されている。また、例えば、特許文献３では、フロー直流法（簡略直流計算法）による潮流計算が開示されている。また、例えば、特許文献４では、フロー交流法による潮流計算が開示されている。また、例えば、特許文献５では、ファースト・デカップルド（Fast Decoupled Load Flow）法による潮流計算が開示されている。
【０００６】
式（１）の電力方程式は非線形方程式であり、これを解くには、未知数に適当な初期値を設定し、電力方程式から近似解を求め、さらにこれを電力方程式に代入して、より精度の高い近似解を求めるような反復計算が主に用いられる。反復計算の手法としては、上記のニュートン・ラフソン法のほか、ガウス・ザイデル（Gauss-Seidel）法も知られているが、一般的にニュートン・ラフソン法の方が丸め誤差が小さく、収束回数が少ないとされ、潮流計算手法としては、ニュートン・ラフソン法が主流となっている。以下にニュートン・ラフソン法の概略を示す。
【０００７】
極座標形式で示された電力方程式の有効電力と無効電力は、電力方程式から、以下の式（２）で示される非線形代数方程式となる。  
0005830370.tif 000003
ここで、ｋはノード番号、δは位相角、｜Ｖ｜は電圧値、Ｐは有効電力、Ｑは無効電力、ｆは非線形代数方程式である。
【０００８】
また、式（２）において、左辺をｙ１〜ｙｎとし、変数をｘ１〜ｘｎとおき、以下の式（３）のように簡易的に示す。  
0005830370.tif 000004
また、式（３）の左辺に対して、初期値をｘ＝ｘ（０）として１次の項までのテイラー展開を行う。  
0005830370.tif 000005
また、ベクトル表示を用いて、式（４）を以下の式（５）のように表す。  
0005830370.tif 000006
ここで、Ｊで表された偏微分導関数の行列は、ヤコビアン行列と呼ばれる。
【０００９】
また、式（５）から、変数ｘに関する式に整理する。  
0005830370.tif 000007
さらに、式（６）をｘ（０）から開始し、反復的にｘ（ｉ）まで計算を行う。  
0005830370.tif 000008
そして、上記計算を繰り返し行うことで、ｙ＝ｆ（ｘ）となる解を計算する。
【００１０】
式（７）を電力方程式に置き換えると、以下の式（８）となる。  
0005830370.tif 000009
ここで、左辺は位相角および電圧値の変化量、右辺はヤコビアン行列の逆行列と電力ミスマッチ行列との積である。そして、式（８）を解くことで、解である各ノードの位相角および電圧値を計算することができる。
【００１１】
このようにして、ニュートン・ラフソン法を用いて潮流計算を行い、各ノードにおける電圧および位相角や、さらに、それらから求まる各送電線や変圧器の有効電力および無効電力潮流などを含む解を求めることができる。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開昭５６−１５９６８１号公報
【特許文献２】特開平１−１３８９２６号公報
【特許文献３】特開平３−７０２３号公報
【特許文献４】特開昭５５−１４４７３２号公報
【特許文献５】特開昭６１−１９９４３０号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
潮流計算をニュートン・ラフソン法で行う場合、対象となる系統によっては解が収束しない場合や、「低め解」と呼ばれる運用上不安定となってしまう解を算出する場合がある。このような場合、従来では、初期値を試行錯誤で変化させて解を収束させたり、解探索時の更新値Δｘに小さな係数（例えば０．１）を乗じて解を収束させたりしていた。  
しかしながら、いずれの場合もノウハウや試行錯誤が必要であり、試行錯誤におけるパラメータ調整など、電力系統の運用者に負荷がかかることとなる。
【課題を解決するための手段】
【００１４】
前述した課題を解決する主たる本発明は、電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算する複数の電力ミスマッチ計算部と、前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から１つを最良の選択解として選択して出力する選択部と、を有し、前記複数の電力ミスマッチ計算部は、前記複数の解に対して、各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出し、前記選択部は、前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの最大値のうちの最小値に対応する解を前記選択解として選択して出力することを特徴とする潮流計算装置である。
【００１５】
本発明の他の特徴については、添付図面及び本明細書の記載により明らかとなる。
【発明の効果】
【００１６】
本発明によれば、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく最良の解を出力することが可能となり、運用者の負荷を軽減することができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の一実施形態における潮流計算装置の構成を示すブロック図である。
【図２】本発明の一実施形態における潮流計算システム全体の構成を示す模式図である。
【図３】本発明の一実施形態における潮流計算装置の各部の機能をコンピュータに実現させるためのプログラムの動作を説明するフローチャートである。
【図４】潮流計算装置の他の構成例を示すブロック図である。
【図５】図４に示した潮流計算装置の各部の機能をコンピュータに実現させるためのプログラムの動作を説明するフローチャートである。
【図６】潮流計算装置のさらに他の構成例を示すブロック図である。
【図７】図６に示した潮流計算装置の各部の機能をコンピュータに実現させるためのプログラムの動作を説明するフローチャートである。
【発明を実施するための形態】
【００１８】
本明細書および添付図面の記載により、少なくとも以下の事項が明らかとなる。
【００１９】
＝＝＝潮流計算装置および潮流計算システムの構成＝＝＝
以下、図１および図２を参照して、本発明の一実施形態における潮流計算装置および潮流計算システムの構成について説明する。
【００２０】
図２は、本実施形態における潮流計算システム全体の構成を示している。図２に示されている潮流計算システムは、キーボード１０１などの入力装置と、ディスプレイ１０２などの表示装置と、系統情報データベース１０３および選択解データベース１０４を格納する記憶装置を備えた、ＰＣなどのコンピュータ１００で構成されている。
【００２１】
コンピュータ１００は、電力系統の運用者によって操作され、キーボード１０１から入力信号ＩＮとして入力される運用者の指示に応じて、系統情報データベース１０３に格納されている系統情報に基づいて潮流計算を行う。そして、コンピュータ１００は、選択解Ｓｓｅｌ（ｋ）を選択解データベース１０４に格納するとともに、解Ｓ１（ｋ）ないしＳ５（ｋ）と、それらに対する電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５とを対応付けてディスプレイ１０２に表示させ、運用者に提示する。すなわち、コンピュータ１００は、図１に示す潮流計算装置１ａの機能を実現する。
【００２２】
図１に示されている潮流計算装置１ａは、潮流計算部１１ないし１５、電力ミスマッチ計算部２１ないし２５、および選択部４１を含んで構成されている。
【００２３】
潮流計算部１１ないし１５は、それぞれ、系統情報データベース１０３を参照して、解Ｓ１（ｋ）ないしＳ５（ｋ）を出力している。また、電力ミスマッチ計算部２１ないし２５には、それぞれ解Ｓ１（ｋ）ないしＳ５（ｋ）が入力され、電力ミスマッチ計算部２１ないし２５からは、電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５が出力されている。そして、選択部４１には、解Ｓ１（ｋ）ないしＳ５（ｋ）、および電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５が入力され、選択部４１は、それらとともに選択解Ｓｓｅｌ（ｋ）を出力している。
【００２４】
＝＝＝潮流計算装置の動作＝＝＝
以下、図３を適宜参照して、本実施形態における潮流計算装置の動作について説明する。
【００２５】
前述したように、潮流計算装置１ａの機能は、コンピュータ１００によって実現することができる。例えば、系統情報データベース１０３を格納する記憶装置を備えるコンピュータに潮流計算プログラムを実行させることによって、潮流計算部１１ないし１５、電力ミスマッチ計算部２１ないし２５、および選択部４１に相当する機能を実現することができる。図３は、これら各部に相当する機能をコンピュータに実現させるためのプログラムの動作を示している。
【００２６】
まず、運用者の指示などに応じてプログラムの処理が開始され（Ｓ１１）、記憶装置に格納された系統情報データベース１０３から、電力系統の状態を示す系統情報を取得する（Ｓ１２）。
【００２７】
次に、系統情報に基づいて、複数の潮流計算手法を用いて潮流計算を行い、それぞれ、少なくとも各ノードｋにおける電圧および位相角を含む解Ｓ１（ｋ）ないしＳ５（ｋ）を求める（Ｓ１３）。本実施形態では、一例として、ニュートン・ラフソン法による解Ｓ１（ｋ）、直流法によるによる解Ｓ２（ｋ）、フロー直流法による解Ｓ３（ｋ）、フロー交流法による解Ｓ４（ｋ）、およびファースト・デカップルド法による解Ｓ５（ｋ）を計算するものとする。
【００２８】
以上のように、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求めることができる。したがって、Ｓ１２およびＳ１３の処理は、潮流計算部１１ないし１５に相当する機能（潮流計算処理）を実現する。
【００２９】
次に、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、それぞれ、各ノードｋにおける電力ミスマッチＭ１（ｋ）ないしＭ５（ｋ）を計算する（Ｓ１４）。より具体的には、複数の潮流計算手法によって得られた解（電圧および位相角）において、式（１）を解くことによって各ノードｋの有効電力Ｐ＿ｃａｌ（ｋ）および無効電力Ｑ＿ｃａｌ（ｋ）を計算し、それぞれの解に対する電力ミスマッチを計算する。
【００３０】
ＰＱ指定ノードでは、各ノードｋに対して有効電力と無効電力とを事前に指定するため、電力ミスマッチは、以下の式（９）のように表される。  
0005830370.tif 000010
ここで、Ｐ＿ｓｅｔ（ｋ）は指定した有効電力、Ｑ＿ｓｅｔ（ｋ）は指定した無効電力である。
【００３１】
一方、ＰＶ指定ノードでは、電力ミスマッチを以下の式（１０）のように計算する。  
0005830370.tif 000011
次に、電力ミスマッチＭ１（ｋ）ないしＭ５（ｋ）から、それぞれの最大値Ｍｍａｘ１ないしＭｍａｘ５を抽出する（Ｓ１５）。
【００３２】
以上のように、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチのうちの最大値Ｍｍａｘ１ないしＭｍａｘ５をそれぞれ抽出することができる。したがって、Ｓ１４およびＳ１５の処理は、電力ミスマッチ計算部２１ないし２５に相当する機能（電力ミスマッチ計算処理）を実現する。なお、電力ミスマッチは、各手法による潮流計算の精度の評価指標となり、その最大値は、各手法において最も精度が低いノードにおける潮流計算の精度を示すこととなる。
【００３３】
次に、電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５のうちの最小値を選択し、解Ｓ１（ｋ）ないしＳ５（ｋ）のうち、当該最小値に対応する解を選択解Ｓｓｅｌ（ｋ）として選択する（Ｓ１６）。例えば、Ｍｍａｘ１が最小である場合、解Ｓ１（ｋ）が選択解Ｓｓｅｌ（ｋ）として選択される。そして、選択解Ｓｓｅｌ（ｋ）とともに、解Ｓ１（ｋ）ないしＳ５（ｋ）、および電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５を出力し（Ｓ１７）、処理を終了する（Ｓ１８）。したがって、Ｓ１６およびＳ１７の処理は、選択部４１に相当する機能（選択処理）を実現する。
【００３４】
このようにして、潮流計算装置１ａは、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチのうちの最大値Ｍｍａｘ１ないしＭｍａｘ５をそれぞれ抽出し、それらの最小値に対応する解を選択解Ｓｓｅｌ（ｋ）として選択して出力する。前述したように、電力ミスマッチの最大値は、各潮流計算手法において最も精度が低いノードにおける潮流計算の精度を示すため、何れかの手法で解が収束しない場合であっても、その最小値に対応する選択解Ｓｓｅｌ（ｋ）を最良の解として安定的に出力することができる。なお、潮流計算装置１ａから出力された選択解Ｓｓｅｌ（ｋ）は、選択解データベース１０４に格納される。
【００３５】
また、潮流計算装置１ａは、解Ｓ１（ｋ）ないしＳ５（ｋ）と、それらに対する電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５とを対応付けてディスプレイ１０２に表示させ、運用者に提示する。したがって、何れかの潮流計算手法で低め解を算出した場合であっても、運用者が当該低め解などを除外して適切な解を選択することができる。
【００３６】
＝＝＝潮流計算装置の他の構成例＝＝＝
上記実施形態では、選択部４１（選択処理）において、電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５のうちの最小値に対応する選択解Ｓｓｅｌ（ｋ）を選択しているが、これに限定されるものではない。
【００３７】
例えば図４に示すように、選択部４１を備えない潮流計算装置１ｂから、解Ｓ１（ｋ）ないしＳ５（ｋ）、および電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５を出力し、それらを対応付けてディスプレイ１０２に表示させる構成としてもよい。この場合、図５に示すように、Ｓ２６の処理によってディスプレイ１０２に表示された電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５から、運用者が手作業で最小値を選択し、解Ｓ１（ｋ）ないしＳ５（ｋ）のうち、当該最小値に対応する解を選択解Ｓｓｅｌ（ｋ）として選択する（Ｓ２７）。なお、図５において、Ｓ２２ないしＳ２５の処理は、図３におけるＳ１２ないしＳ１５の処理と同様である。
【００３８】
また、例えば図６に示すように、電力ミスマッチ計算部２１ないし２５に代えて、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチの平均値Ｍａｖｅ１ないしＭａｖｅ５をそれぞれ算出する電力ミスマッチ計算部３１ないし３５を備える構成としてもよい。当該潮流計算装置１ｃでは、図７に示すように、Ｓ３５の処理によって算出された電力ミスマッチの平均値Ｍａｖｅ１ないしＭａｖｅ５から最小値を選択し、解Ｓ１（ｋ）ないしＳ５（ｋ）のうち、当該最小値に対応する解を選択解Ｓｓｅｌ（ｋ）として選択する（Ｓ３６）。そして、選択解Ｓｓｅｌ（ｋ）とともに、解Ｓ１（ｋ）ないしＳ５（ｋ）、および電力ミスマッチの平均値Ｍａｖｅ１ないしＭａｖｅ５を出力し（Ｓ３７）、処理を終了する（Ｓ３８）。なお、図７において、Ｓ３２ないしＳ３４の処理は、図３におけるＳ１２ないしＳ１４の処理と同様である。
【００３９】
前述したように、潮流計算装置１ａにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求め、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチをそれぞれ計算し、それらに基づいて選択解Ｓｓｅｌ（ｋ）を選択することによって、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、各手法による潮流計算の精度に基づく最良の解を安定的に出力することが可能となり、運用者の負荷を軽減することができる。
【００４０】
また、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチのうちの最大値Ｍｍａｘ１ないしＭｍａｘ５をそれぞれ抽出し、それらの最小値に対応する解を選択解Ｓｓｅｌ（ｋ）として選択することによって、何れかの手法で解が収束しない場合であっても、各潮流計算手法において最も精度が低いノードにおける潮流計算の精度に基づいて、選択解Ｓｓｅｌ（ｋ）を最良の解として安定的に出力することができる。さらに、解Ｓ１（ｋ）ないしＳ５（ｋ）と、それらに対する電力ミスマッチの最大値Ｍｍａｘ１ないしＭｍａｘ５とを対応付けてディスプレイ１０２に表示させ、運用者に提示することによって、何れかの潮流計算手法で低め解を算出した場合であっても、運用者が当該低め解などを除外して適切な解を選択することができる。
【００４１】
また、潮流計算装置１ｃにおいて、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチの平均値Ｍａｖｅ１ないしＭａｖｅ５をそれぞれ算出し、それらの最小値に対応する解を選択解Ｓｓｅｌ（ｋ）として選択することによって、何れかの手法で解が収束しない場合であっても、各潮流計算手法における平均的な潮流計算の精度に基づいて、選択解Ｓｓｅｌ（ｋ）を最良の解として安定的に出力することができる。さらに、解Ｓ１（ｋ）ないしＳ５（ｋ）と、それらに対する電力ミスマッチの平均値Ｍａｖｅ１ないしＭａｖｅ５とを対応付けてディスプレイ１０２に表示させ、運用者に提示することによって、何れかの潮流計算手法で低め解を算出した場合であっても、運用者が当該低め解などを除外して適切な解を選択することができる。
【００４２】
また、前述したように、潮流計算装置１ｂにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求めて出力することによって、ニュートン・ラフソン法で解が収束しない場合や、低め解を算出した場合であっても、それ以外の手法を評価し、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、運用者が最良の解を選択することが可能となり、運用者の負荷を軽減することができる。
【００４３】
また、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチのうちの最大値Ｍｍａｘ１ないしＭｍａｘ５をそれぞれ抽出し、それらを解Ｓ１（ｋ）ないしＳ５（ｋ）と対応付けてディスプレイ１０２に表示させ、運用者に提示することによって、何れかの手法で解が収束しない場合や、低め解を算出した場合であっても、運用者が低め解などを適切に除外しつつ、各潮流計算手法において最も精度が低いノードにおける潮流計算の精度に基づいて、選択解Ｓｓｅｌ（ｋ）を最良の解として選択することができる。
【００４４】
また、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチの平均値Ｍａｖｅ１ないしＭａｖｅ５をそれぞれ算出し、それらを解Ｓ１（ｋ）ないしＳ５（ｋ）と対応付けてディスプレイ１０２に表示させ、運用者に提示することによって、何れかの手法で解が収束しない場合や、低め解を算出した場合であっても、運用者が低め解などを適切に除外しつつ、各潮流計算手法における平均的な潮流計算の精度に基づいて、選択解Ｓｓｅｌ（ｋ）を最良の解として選択することができる。
【００４５】
また、潮流計算装置１ａの機能をコンピュータに実現させるためのプログラムにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求め、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチをそれぞれ計算し、それらに基づいて選択解Ｓｓｅｌ（ｋ）を選択することによって、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、各手法による潮流計算の精度に基づく最良の解を安定的に出力することが可能となり、運用者の負荷を軽減することができる。
【００４６】
また、潮流計算装置１ｂの機能をコンピュータに実現させるためのプログラムにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求めて出力することによって、ニュートン・ラフソン法で解が収束しない場合や、低め解を算出した場合であっても、それ以外の手法を評価し、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、運用者が最良の解を選択することが可能となり、運用者の負荷を軽減することができる。
【００４７】
また、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求め、解Ｓ１（ｋ）ないしＳ５（ｋ）に対して、各ノードｋにおける電力ミスマッチをそれぞれ計算し、それらに基づいて選択解Ｓｓｅｌ（ｋ）を選択することによって、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、各手法による潮流計算の精度に基づく最良の解を安定的に出力することが可能となり、運用者の負荷を軽減することができる。
【００４８】
また、系統情報に基づいて、互いに異なる潮流計算手法による解Ｓ１（ｋ）ないしＳ５（ｋ）をそれぞれ求めて出力することによって、ニュートン・ラフソン法で解が収束しない場合や、低め解を算出した場合であっても、それ以外の手法を評価し、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、運用者が最良の解を選択することが可能となり、運用者の負荷を軽減することができる。
【００４９】
なお、上記実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、その趣旨を逸脱することなく、変更、改良され得るとともに、本発明にはその等価物も含まれる。
【符号の説明】
【００５０】
１ａ〜１ｃ  潮流計算装置
１１〜１５  潮流計算部
２１〜２５、３１〜３５  電力ミスマッチ計算部
４１        選択部
１００      コンピュータ
１０１      キーボード
１０２      ディスプレイ
１０３      系統情報データベース
１０４      選択解データベース
【図１】
0005830370.tif 000012
【図２】
0005830370.tif 000013
【図３】
0005830370.tif 000014
【図４】
0005830370.tif 000015
【図５】
0005830370.tif 000016
【図６】
0005830370.tif 000017
【図７】
0005830370.tif 000018