電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算する複数の電力ミスマッチ計算部と、前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から1つを最良の選択解として選択して出力する選択部と、を有することを特徴とする潮流計算装置。前記複数の電力ミスマッチ計算部は、前記複数の解に対して、各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出し、前記選択部は、前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの最大値のうちの最小値に対応する解を前記選択解として選択して出力することを特徴とする請求項1に記載の潮流計算装置。前記複数の解と、前記複数の解に対する前記電力ミスマッチの最大値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする請求項2に記載の潮流計算装置。前記複数の電力ミスマッチ計算部は、前記複数の解に対して、各ノードにおける電力ミスマッチの平均値をそれぞれ算出し、前記選択部は、前記複数の解のうち、前記複数の解に対する前記電力ミスマッチの平均値のうちの最小値に対応する解を前記選択解として選択して出力することを特徴とする請求項1に記載の潮流計算装置。前記複数の解と、前記複数の解に対する前記電力ミスマッチの平均値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする請求項4に記載の潮流計算装置。電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部を有し、前記複数の解を出力することを特徴とする潮流計算装置。前記複数の解に対して、各ノードにおける電力ミスマッチを計算し、計算された各ノードにおける電力ミスマッチのうちの最大値をそれぞれ抽出する複数の電力ミスマッチ計算部をさらに有し、前記複数の解と、前記複数の解に対する前記電力ミスマッチの最大値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする請求項6に記載の潮流計算装置。前記複数の解に対して、各ノードにおける電力ミスマッチを計算し、計算された各ノードにおける電力ミスマッチの平均値をそれぞれ算出する複数の電力ミスマッチ計算部をさらに有し、前記複数の解と、前記複数の解に対する前記電力ミスマッチの平均値とを対応付けて表示装置に表示させ、電力系統運用者に提示することを特徴とする請求項6に記載の潮流計算装置。電力系統の状態を示す系統情報を格納する記憶装置を備えるコンピュータに、前記系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める潮流計算処理と、前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算する電力ミスマッチ計算処理と、前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から1つを最良の選択解として選択して出力する選択処理と、を実行させることを特徴とする潮流計算プログラム。電力系統の状態を示す系統情報を格納する記憶装置を備えるコンピュータに、前記系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める潮流計算処理と、前記複数の解を出力する出力処理と、を実行させることを特徴とする潮流計算プログラム。電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求め、前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算し、前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から1つを最良の選択解として選択して出力することを特徴とする潮流計算方法。電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求め、前記複数の解を出力することを特徴とする潮流計算方法。潮流計算をニュートンラフソン法で行う場合、対象となる系統によっては解が収束しない場合や、低め解と呼ばれる運用上不安定となってしまう解を算出する場合がある。このような場合、従来では、初期値を試行錯誤で変化させて解を収束させたり、解探索時の更新値Δxに小さな係数例えば0.1を乗じて解を収束させたりしていた。しかしながら、いずれの場合もノウハウや試行錯誤が必要であり、試行錯誤におけるパラメータ調整など、電力系統の運用者に負荷がかかることとなる。前述した課題を解決する主たる本発明は、電力系統の状態を示す系統情報に基づいて、互いに異なる潮流計算手法を用いて潮流計算を行い、少なくとも各ノードにおける電圧および位相角を含む複数の解をそれぞれ求める複数の潮流計算部と、前記複数の解に対して、各ノードにおける電力ミスマッチをそれぞれ計算する複数の電力ミスマッチ計算部と、前記複数の解に対する前記電力ミスマッチに基づいて、前記複数の解から1つを最良の選択解として選択して出力する選択部と、を有することを特徴とする潮流計算装置である。本発明の他の特徴については、添付図面及び本明細書の記載により明らかとなる。本発明によれば、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく最良の解を出力することが可能となり、運用者の負荷を軽減することができる。本明細書および添付図面の記載により、少なくとも以下の事項が明らかとなる。===潮流計算装置および潮流計算システムの構成===以下、図1および図2を参照して、本発明の一実施形態における潮流計算装置および潮流計算システムの構成について説明する。図2は、本実施形態における潮流計算システム全体の構成を示している。図2に示されている潮流計算システムは、キーボード101などの入力装置と、ディスプレイ102などの表示装置と、系統情報データベース103および選択解データベース104を格納する記憶装置を備えた、PCなどのコンピュータ100で構成されている。コンピュータ100は、電力系統の運用者によって操作され、キーボード101から入力信号INとして入力される運用者の指示に応じて、系統情報データベース103に格納されている系統情報に基づいて潮流計算を行う。そして、コンピュータ100は、選択解Sselkを選択解データベース104に格納するとともに、解S1kないしS5kと、それらに対する電力ミスマッチの最大値Mmax1ないしMmax5とを対応付けてディスプレイ102に表示させ、運用者に提示する。すなわち、コンピュータ100は、図1に示す潮流計算装置1aの機能を実現する。図1に示されている潮流計算装置1aは、潮流計算部11ないし15、電力ミスマッチ計算部21ないし25、および選択部41を含んで構成されている。潮流計算部11ないし15は、それぞれ、系統情報データベース103を参照して、解S1kないしS5kを出力している。また、電力ミスマッチ計算部21ないし25には、それぞれ解S1kないしS5kが入力され、電力ミスマッチ計算部21ないし25からは、電力ミスマッチの最大値Mmax1ないしMmax5が出力されている。そして、選択部41には、解S1kないしS5k、および電力ミスマッチの最大値Mmax1ないしMmax5が入力され、選択部41は、それらとともに選択解Sselkを出力している。===潮流計算装置の動作===以下、図3を適宜参照して、本実施形態における潮流計算装置の動作について説明する。前述したように、潮流計算装置1aの機能は、コンピュータ100によって実現することができる。例えば、系統情報データベース103を格納する記憶装置を備えるコンピュータに潮流計算プログラムを実行させることによって、潮流計算部11ないし15、電力ミスマッチ計算部21ないし25、および選択部41に相当する機能を実現することができる。図3は、これら各部に相当する機能をコンピュータに実現させるためのプログラムの動作を示している。まず、運用者の指示などに応じてプログラムの処理が開始されS11、記憶装置に格納された系統情報データベース103から、電力系統の状態を示す系統情報を取得するS12。次に、系統情報に基づいて、複数の潮流計算手法を用いて潮流計算を行い、それぞれ、少なくとも各ノードkにおける電圧および位相角を含む解S1kないしS5kを求めるS13。本実施形態では、一例として、ニュートンラフソン法による解S1k、直流法によるによる解S2k、フロー直流法による解S3k、フロー交流法による解S4k、およびファーストデカップルド法による解S5kを計算するものとする。以上のように、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求めることができる。したがって、S12およびS13の処理は、潮流計算部11ないし15に相当する機能潮流計算処理を実現する。次に、解S1kないしS5kに対して、それぞれ、各ノードkにおける電力ミスマッチM1kないしM5kを計算するS14。より具体的には、複数の潮流計算手法によって得られた解電圧および位相角において、式1を解くことによって各ノードkの有効電力P_calkおよび無効電力Q_calkを計算し、それぞれの解に対する電力ミスマッチを計算する。PQ指定ノードでは、各ノードkに対して有効電力と無効電力とを事前に指定するため、電力ミスマッチは、以下の式9のように表される。ここで、P_setkは指定した有効電力、Q_setkは指定した無効電力である。一方、PV指定ノードでは、電力ミスマッチを以下の式10のように計算する。次に、電力ミスマッチM1kないしM5kから、それぞれの最大値Mmax1ないしMmax5を抽出するS15。以上のように、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチのうちの最大値Mmax1ないしMmax5をそれぞれ抽出することができる。したがって、S14およびS15の処理は、電力ミスマッチ計算部21ないし25に相当する機能電力ミスマッチ計算処理を実現する。なお、電力ミスマッチは、各手法による潮流計算の精度の評価指標となり、その最大値は、各手法において最も精度が低いノードにおける潮流計算の精度を示すこととなる。次に、電力ミスマッチの最大値Mmax1ないしMmax5のうちの最小値を選択し、解S1kないしS5kのうち、当該最小値に対応する解を選択解Sselkとして選択するS16。例えば、Mmax1が最小である場合、解S1kが選択解Sselkとして選択される。そして、選択解Sselkとともに、解S1kないしS5k、および電力ミスマッチの最大値Mmax1ないしMmax5を出力しS17、処理を終了するS18。したがって、S16およびS17の処理は、選択部41に相当する機能選択処理を実現する。このようにして、潮流計算装置1aは、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチのうちの最大値Mmax1ないしMmax5をそれぞれ抽出し、それらの最小値に対応する解を選択解Sselkとして選択して出力する。前述したように、電力ミスマッチの最大値は、各潮流計算手法において最も精度が低いノードにおける潮流計算の精度を示すため、何れかの手法で解が収束しない場合であっても、その最小値に対応する選択解Sselkを最良の解として安定的に出力することができる。なお、潮流計算装置1aから出力された選択解Sselkは、選択解データベース104に格納される。また、潮流計算装置1aは、解S1kないしS5kと、それらに対する電力ミスマッチの最大値Mmax1ないしMmax5とを対応付けてディスプレイ102に表示させ、運用者に提示する。したがって、何れかの潮流計算手法で低め解を算出した場合であっても、運用者が当該低め解などを除外して適切な解を選択することができる。===潮流計算装置の他の構成例===上記実施形態では、選択部41選択処理において、電力ミスマッチの最大値Mmax1ないしMmax5のうちの最小値に対応する選択解Sselkを選択しているが、これに限定されるものではない。例えば図4に示すように、選択部41を備えない潮流計算装置1bから、解S1kないしS5k、および電力ミスマッチの最大値Mmax1ないしMmax5を出力し、それらを対応付けてディスプレイ102に表示させる構成としてもよい。この場合、図5に示すように、S26の処理によってディスプレイ102に表示された電力ミスマッチの最大値Mmax1ないしMmax5から、運用者が手作業で最小値を選択し、解S1kないしS5kのうち、当該最小値に対応する解を選択解Sselkとして選択するS27。なお、図5において、S22ないしS25の処理は、図3におけるS12ないしS15の処理と同様である。また、例えば図6に示すように、電力ミスマッチ計算部21ないし25に代えて、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチの平均値Mave1ないしMave5をそれぞれ算出する電力ミスマッチ計算部31ないし35を備える構成としてもよい。当該潮流計算装置1cでは、図7に示すように、S35の処理によって算出された電力ミスマッチの平均値Mave1ないしMave5から最小値を選択し、解S1kないしS5kのうち、当該最小値に対応する解を選択解Sselkとして選択するS36。そして、選択解Sselkとともに、解S1kないしS5k、および電力ミスマッチの平均値Mave1ないしMave5を出力しS37、処理を終了するS38。なお、図7において、S32ないしS34の処理は、図3におけるS12ないしS14の処理と同様である。前述したように、潮流計算装置1aにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求め、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチをそれぞれ計算し、それらに基づいて選択解Sselkを選択することによって、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、各手法による潮流計算の精度に基づく最良の解を安定的に出力することが可能となり、運用者の負荷を軽減することができる。また、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチのうちの最大値Mmax1ないしMmax5をそれぞれ抽出し、それらの最小値に対応する解を選択解Sselkとして選択することによって、何れかの手法で解が収束しない場合であっても、各潮流計算手法において最も精度が低いノードにおける潮流計算の精度に基づいて、選択解Sselkを最良の解として安定的に出力することができる。さらに、解S1kないしS5kと、それらに対する電力ミスマッチの最大値Mmax1ないしMmax5とを対応付けてディスプレイ102に表示させ、運用者に提示することによって、何れかの潮流計算手法で低め解を算出した場合であっても、運用者が当該低め解などを除外して適切な解を選択することができる。また、潮流計算装置1cにおいて、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチの平均値Mave1ないしMave5をそれぞれ算出し、それらの最小値に対応する解を選択解Sselkとして選択することによって、何れかの手法で解が収束しない場合であっても、各潮流計算手法における平均的な潮流計算の精度に基づいて、選択解Sselkを最良の解として安定的に出力することができる。さらに、解S1kないしS5kと、それらに対する電力ミスマッチの平均値Mave1ないしMave5とを対応付けてディスプレイ102に表示させ、運用者に提示することによって、何れかの潮流計算手法で低め解を算出した場合であっても、運用者が当該低め解などを除外して適切な解を選択することができる。また、前述したように、潮流計算装置1bにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求めて出力することによって、ニュートンラフソン法で解が収束しない場合や、低め解を算出した場合であっても、それ以外の手法を評価し、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、運用者が最良の解を選択することが可能となり、運用者の負荷を軽減することができる。また、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチのうちの最大値Mmax1ないしMmax5をそれぞれ抽出し、それらを解S1kないしS5kと対応付けてディスプレイ102に表示させ、運用者に提示することによって、何れかの手法で解が収束しない場合や、低め解を算出した場合であっても、運用者が低め解などを適切に除外しつつ、各潮流計算手法において最も精度が低いノードにおける潮流計算の精度に基づいて、選択解Sselkを最良の解として選択することができる。また、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチの平均値Mave1ないしMave5をそれぞれ算出し、それらを解S1kないしS5kと対応付けてディスプレイ102に表示させ、運用者に提示することによって、何れかの手法で解が収束しない場合や、低め解を算出した場合であっても、運用者が低め解などを適切に除外しつつ、各潮流計算手法における平均的な潮流計算の精度に基づいて、選択解Sselkを最良の解として選択することができる。また、潮流計算装置1aの機能をコンピュータに実現させるためのプログラムにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求め、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチをそれぞれ計算し、それらに基づいて選択解Sselkを選択することによって、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、各手法による潮流計算の精度に基づく最良の解を安定的に出力することが可能となり、運用者の負荷を軽減することができる。また、潮流計算装置1bの機能をコンピュータに実現させるためのプログラムにおいて、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求めて出力することによって、ニュートンラフソン法で解が収束しない場合や、低め解を算出した場合であっても、それ以外の手法を評価し、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、運用者が最良の解を選択することが可能となり、運用者の負荷を軽減することができる。また、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求め、解S1kないしS5kに対して、各ノードkにおける電力ミスマッチをそれぞれ計算し、それらに基づいて選択解Sselkを選択することによって、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、各手法による潮流計算の精度に基づく最良の解を安定的に出力することが可能となり、運用者の負荷を軽減することができる。また、系統情報に基づいて、互いに異なる潮流計算手法による解S1kないしS5kをそれぞれ求めて出力することによって、ニュートンラフソン法で解が収束しない場合や、低め解を算出した場合であっても、それ以外の手法を評価し、解を収束させるために試行錯誤でのパラメータ調整などを行うことなく、運用者が最良の解を選択することが可能となり、運用者の負荷を軽減することができる。なお、上記実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、その趣旨を逸脱することなく、変更、改良され得るとともに、本発明にはその等価物も含まれる。