(57)【要約】【課題】状態推定計算が収束する可能性を高めることができるようにした電力系統の状態推定装置を提供すること。【解決手段】電力系統の状態推定装置１は、電力系統に設けられる各センサの計測値のうち、電力系統についての状態推定計算が収束すると予測される所定の計測値の組合せを選択する計測値選択部１０２と、計測値選択部の選択した所定の計測値の組合せと、電力系統の系統モデルとに基づいて、電力系統の状態を推定する状態推定計算部１０１と、を備える。【選択図】図１
2018098986.tif 000002
【特許請求の範囲】
【請求項１】
電力系統の状態を推定する状態推定装置であって、
前記電力系統に設けられる各センサの計測値のうち、前記電力系統についての状態推定計算が収束すると予測される所定の計測値の組合せを選択する計測値選択部と、
前記計測値選択部の選択した前記所定の計測値の組合せと、前記電力系統の系統モデルとに基づいて、前記電力系統の状態を推定する状態推定計算部と、
を備える電力系統の状態推定装置。
【請求項２】
前記計測値選択部は、前記各センサの計測値の組合せを変更しながら前記状態推定計算部での計算を行わせることにより、前記所定の計測値の組合せを、前記状態推定計算部の計算結果が収束する計測値の組合せとして選択する、
請求項１に記載の電力系統の状態推定装置。
【請求項３】
前記各センサは、センサ属性に応じた優先度が設定されており、
前記計測値選択部は、前記各センサの中から前記優先度に基づいて複数のセンサを抽出し、前記抽出した複数のセンサの計測値の組合せを変更しながら前記状態推定計算部での計算を行うことにより、前記状態推定計算部の計算結果が収束する計測値の組合せを前記所定の計測値の組合せとして選択する、
請求項２に記載の電力系統の状態推定装置。
【請求項４】
前記優先度は、精度の高いセンサほど高優先度となるように設定されており、
前記計測値選択部は、前記各センサの中から前記優先度の高いセンサを抽出する、
請求項３に記載の電力系統の状態推定装置。
【請求項５】
前記計測値選択部は、前記各センサのうち優先度が高い高精度センサの計測値は前記所定の計測値の組合せの一部として全て選択し、前記各センサのうち前記高精度センサ以外のセンサの計測値は、その精度に応じて前記所定の計測値の組合せの他の一部として部分的に選択する、
請求項４に記載の電力系統の状態推定装置。
【請求項６】
前記計測値選択部は、前記状態推定計算部の計算した推定値と計測値との乖離を評価指標として算出し、前記評価指標に基づいて、前記各計測値の中から前記所定の計測値の組合せを選択する、
請求項１〜５のいずれか一項に記載の電力系統の状態推定装置。
【請求項７】
前記計測値選択部による選択結果を出力する出力部を備える、
請求項６に記載の電力系統の状態推定装置。
【請求項８】
電力系統の状態を計算機を用いて推定する状態推定方法であって、
前記計算機は、
前記電力系統に設けられる各センサの計測値のうち、前記電力系統についての状態推定計算が収束すると予測される所定の計測値の組合せを選択し、
前記選択した所定の計測値の組合せと、前記電力系統の系統モデルとに基づいて、前記電力系統の状態を推定する、
電力系統の状態推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統の状態推定装置および状態推定方法に関する。
【背景技術】
【０００２】
電力系統の潮流状態（有効電力、無効電力、電圧、電圧位相）を把握することは、電力系統の監視制御に有効であるとともに、潮流計算等に使用する解析モデルを構築するためにも有用である。
【０００３】
現状の電力系統に近い状態を解析モデルで再現することで、今後起こりうる系統状態を予測したり、系統事故に備えた安定化対策を実施したりできる。一般に、電力系統に設置された計測センサで電気量（有効電力、無効電力、電圧値、電圧位相、電流値など）の計測値を測定し、それら計測値を用いて状態推定計算を行うことで、電力系統の状態を把握することができる。以下、電力系統を系統と略記する場合がある。
【０００４】
系統の状態を把握するためには、系統の電力方程式を解くために十分な数の計測値を用意し、可観測状態とすることが重要である。しかし、系統の構成や潮流条件などによっては、状態推定計算が収束しない場合が生じうる。例えば、線路の抵抗成分が大きい区間が存在する場合、無効電力潮流が大きい区間が存在する場合、隣接する計測値の例えば電圧値が大きく異なる場合などで、状態推定計算が収束しない可能性がある。
【０００５】
仮に、状態推定計算が収束した場合であっても、計測値と推定値とに大きい乖離（残差）が発生し、得られた推定結果（潮流状態）の信頼性が低下するという課題もある。
【０００６】
非特許文献１には、計測値と推定値との間に大きい乖離が生じる場合、乖離の生じた計測値を除外して、状態推定計算を再度実施する手法が示されている。
【０００７】
特許文献１には、状態推定の方程式を数学的に解くことができない場合でも、推定できる範囲の状態を推定する手法が示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０００−３０８２５８号公報
【非特許文献】
【０００９】
【非特許文献１】「電力系統の状態推定における不良データ検出法」，電気学会論文誌Ｂ，Ｖｏｌ．９８，Ｎｏ．１０（１９８０）
【発明の概要】
【発明が解決しようとする課題】
【００１０】
特許文献１に記載の従来技術には、部分的な系統の状態推定結果を計算することは可能であるが、元の系統の状態量をすべて求めることはできない。
【００１１】
非特許文献１に記載の従来技術では、状態推定が収束しない場合に、非収束となる原因の計測値がいずれであるかを判別できず、状態推定結果を得ることができない。
【００１２】
本発明は、上記課題に鑑みてなされたもので、その目的は、状態推定計算が収束する可能性を高めることができるようにした電力系統の状態推定装置および状態推定方法を提供することにある。
【課題を解決するための手段】
【００１３】
上記課題を解決すべく、本発明に従う電力系統の状態推定装置は、電力系統に設けられる各センサの計測値のうち、電力系統についての状態推定計算が収束すると予測される所定の計測値の組合せを選択する計測値選択部と、計測値選択部の選択した所定の計測値の組合せと、電力系統の系統モデルとに基づいて、電力系統の状態を推定する状態推定計算部と、を備える。
【発明の効果】
【００１４】
本発明によれば、各計測値の中から状態推定計算が収束すると予測される所定の計測値の組合せを選択して電力系統の状態を推定するため、状態推定計算の収束する可能性を高めて、使い勝手を向上することができる。
【図面の簡単な説明】
【００１５】
【図１】電力系統の状態推定装置の機能ブロック図である。
【図２】電力系統の解析モデルの例である。
【図３】電力系統の状態推定装置の構成図である。
【図４】状態推定の全体処理を示すフローチャートである。
【図５】状態推定計算処理のフローチャートである。
【図６】計測値のボトルネック指標を計算する処理のフローチャートである。
【図７】計測値の組合せを変更することで最適な組合せを探索する処理を示すフローチャートである。
【図８】計測値の組合せを探索する例を示すテーブルである。
【発明を実施するための形態】
【００１６】
以下、図面に基づいて、本発明の実施の形態を説明する。本実施形態は、例えば、電力系統の監視制御システム、潮流計算の入力データ作成システム、各種系統計算の初期状態作成システムの機能として活用できる。また、得られた状態推定の計算結果を初期値とする安定度解析システムまたは最適潮流計算システムへ本実施形態を用いることで、例えば、オンラインの信頼度監視システム、安定度監視システム、対策検討システム、中央給電指令所システムなどでも活用することができる。
【００１７】
本実施形態では、後述するように、電力系統の状態推定計算の入力データとして用いる計測値について、状態推定計算が収束し、より確からしい推定結果を得るための、計測値の組合せを求める。
【００１８】
このため本実施形態の電力系統の状態を推定する装置は、系統モデルデータベースＤＢ１と、計測値データベースＤＢ２と、入力データ補正データベースＤＢ３と、結果データベースＤＢ４とを持ち、さらに、状態推定計算部１０１と、入力データ計測値組合せ探索部１０２と、ボトルネック指標計算部１０３と、計算結果記録部１０４とを有する。
【００１９】
本実施形態では、状態推定に用いる計測値の組合せ、または除外する計測値の組合せを自動的に求めることができる。これにより、本実施形態では、状態推定計算が収束しない事態の発生を回避できる可能性を高めることができる。さらに、本実施形態では、入力データの作成に必要となる労力および時間を短縮可能である。さらに、本実施形態では、計測値の組合せを最適にできるため、計測値と推定値の差である推定残差をより小さくすることができ、より確からしい推定結果が得られる。さらに、本実施形態では、状態推定計算が収束しない原因となる不良データを判別することもできるため、センサの故障などの不具合を把握することもできる。
【実施例１】
【００２０】
図１〜図８を用いて実施例を説明する。以下に述べる実施例は、一つの例であって、本発明は実施例の構成に限定されない。
【００２１】
図１を用いて、電力系統の状態推定装置１の処理機能構成について説明する。状態推定装置１は、図３で後述するように、計算機を用いて構成することができる。状態推定装置１は、例えば、状態推定計算部１０１と、入力データ計測値組合せ探索部１０２と、ボトルネック指標計算部１０３と、計算結果記録部１０４と、系統モデルデータベースＤＢ１と、計測値データベースＤＢ２と、入力データ補正データベースＤＢ３と、結果データベースＤＢ４とを含む。
【００２２】
状態推定計算部１０１は、電力系統の状態推定計算手法として確立された一般的なアルゴリズムを用いて、電力系統の状態を推定する。入力データ計測値組合せ探索部１０２は、「計測値選択部」の例である。以下、組合せ探索部１０２と略記する場合がある。ボトルネック指標計算部１０３は、後述するように、「評価指標」の例であるボトルネック指標ΔＦｍｎを計算する。計算結果記録部１０４は、状態推定装置１での計算結果を結果ＤＢ４へ出力する。それら計算結果の一部または全部を、ユーザの指示に応じて、ディスプレイに表示させてもよい。
【００２３】
なお、図１に示す処理機能１０１〜１０４およびデータベースＤＢ１〜ＤＢ４間の接続関係、接続方向の矢印の向きは一例であり、図示の構成に限定されない。
【００２４】
図２を用いて、電力系統の状態を推定する手法の概要を説明する。図２は、電力系統の解析モデルの一例を示す。
【００２５】
系統解析モデルが表す電力系統は、例えば、ノード（母線）２０１と、ノード２０１に接続される発電機２０２と、各ノード２０１を接続する送電線路２０３と、ノード２０１に接続される負荷２０４とを含んで構成される。なお、点線の四角で囲んだ数字は、便宜上設定したノード番号を示す。
【００２６】
各ノード付近の線路には、実系統では該当場所にセンサ２０５，２０５Ａが設置されている。センサ２０５では、例えば線路の有効電力、無効電力、電圧、電流が測定されており、それら計測値のデータは図示せぬ監視制御装置を介して収集されている。
【００２７】
センサ２０５Ａは、通常のセンサ２０５よりも高精度なセンサである。高精度のセンサ２０５Ａには、例えばPhase Measurement Unit（ＰＭＵ）のような、高精度かつ電圧位相情報も計測できるような計測センサを用いることができる。
【００２８】
図２に示すような電力系統の解析モデルを構築して、潮流を計算するためには、線路の接続構成や線路インピーダンスのような系統構成情報に加えて、発電機２０２の発電量（有効電力、無効電力）や端子電圧、負荷２０４の有効・無効電力量が必要である。
【００２９】
現状の系統に近い状態の解析モデルを作成するためには、系統の計測センサ２０５および２０５Ａの計測情報を元に、前述の電気量を推定する必要がある。ここで、各センサ２０５，２０５Ａの計測値には、ＣＴ（Current Transformer：変流器）やＰＴ（Potential Transformer：計器用変圧器）の誤算、デジタルサンプリング、実効値演算等に起因する誤差が、通常含まれている。したがって、系統の電力方程式を解くために最低限必要となる計測情報よりも多くのセンサからの計測情報から、各計測値と計算値の二乗誤差が最小となる潮流状態を求める状態推定計算が行われる。
【００３０】
状態推定計算では、系統の構成や潮流条件によっては、状態推定計算が収束しない場合が生じうる。計算が収束した場合でも、計測値と推定値に大きい乖離（残差）が発生し、得られた推定結果（潮流状態）の信頼性が低下する場合がある。状態推定計算を収束させるためには、多くの計測値の中から、状態推定計算が収束し、かつより確からしい推定結果を得るための、状態推定計算の入力値の組合せを求めることが重要となる。以下、適切な入力値（計測値）の組合せを得るための構成例を説明する。
【００３１】
図３は、電力系統の状態推定装置１の構成例を示す。状態推定装置１は、計算機システムで構成されている。状態推定装置１は、例えば、表示装置１１、キーボードやマウス等の入力装置１２、ＣＰＵ１３、通信装置１４、メモリ１５、および各種データベースＤＢ１〜ＤＢ５がバス線１６に接続されている。
【００３２】
ＣＰＵ（Central Processing Unit ）１３は、所定のコンピュータプログラムを実行することで、後述する所定の各機能を実現する。ＣＰＵ１３は、表示装置１１に表示すべき画像データの生成を指示したり、各種データベース内のデータ検索等を行ったりする。
【００３３】
メモリ１５は、各種データを記憶する記憶装置である。図示するメモリ１５は、例えば、ＲＡＭ（Random access memory）のみならず、ハードディスクドライブやフラッシュメモリデバイスなどの補助記憶装置を含んでもよい。
【００３４】
メモリ１５は、系統構成や系統定数などの系統モデルデータ、現在や過去の計測値データ、状態推定に使用する入力データの補正データ、演算結果データなどを格納する。入力データの補正データ（入力データ補正データとも呼ぶ）は、例えば、使用優先度や計測精度良否などの計測値の属性を示すデータである。演算結果データは、例えば、状態推定計算の結果、使用した計測値の組合せ、ボトルネックとなる計測値の指標などを含む。
【００３５】
ＣＰＵ１３は、メモリ１５に格納されたデータに基づいて必要な画像データを生成し、生成した画像データを表示装置１１に表示させる。
【００３６】
表示装置１１は、状態推定装置１がシステム管理者などのユーザに対して、情報を提供するための装置である。表示装置１１は、例えば、ディスプレイから構成される。状態推定装置１は、ディスプレイに加えて、プリンタや音声出力装置などを備えてもよい。
【００３７】
入力装置１２は、ユーザが状態推定装置１にデータなどを入力するための装置である。入力装置１２は、例えば、キーボード、マウス、タッチパネル、音声認識装置などから構成される。
【００３８】
通信装置１４は、状態推定装置１が通信網を介して外部計算機や外部装置（いずれも不図示）と通信するための通信インターフェース回路である。
【００３９】
ユーザとの情報交換に用いる操作用端末を状態推定装置１と別に設けてもよい。ユーザは、操作用端末を用いて状態推定装置１にアクセスすることで、状態推定装置１と情報を交換することもできる。
【００４０】
電力系統の状態推定装置１には、大きく分けて５つのデータベースＤＢ１〜ＤＢ５が搭載されている。
【００４１】
系統モデルデータベースＤＢ１には、線路２０３や発電機２０２などの電力系統を構成する各設備に関するデータが記憶されている。線路に関するデータとしては、例えば、抵抗、リアクタンス、対地静電容量がある。発電機に関するデータとしては、例えば、容量、過渡リアクタンスなどがある。これら系統構成に関するデータは、電力系統の状態推定計算を行う場合に必要とされる。
【００４２】
計測値データベースＤＢ２には、過去の計測値のデータと現在の計測値のデータとが記憶されている。
【００４３】
入力データ補正データベースＤＢ３には、計測値の属性などの、状態推定のために入力するデータ（入力データ）に関連するデータが記憶されている。計測値の属性には、例えば、その計測値を使用する際に考慮する優先度、計測精度の良否などがある。
【００４４】
結果データベースＤＢ４には、電力系統の状態推定演算結果である系統潮流状態、系統解析モデルデータ、計測値補正値、計測値毎のボトルネック指標値、計測値を取得するための優先順位となるポーリング順位情報が記憶されている。結果データベースＤＢ４に格納されたデータにより、状態推定計算結果を確認することができるとともに、潮流計算等の入力データを準備することが可能となる。
【００４５】
プログラムデータベースＤＢ５は、状態推定プログラムＰＲ１、計算値ボトルネック指標計算プログラムＰＲ２、計測値組合せ探索プログラムＰＲ３を格納する。これらのコンピュータプログラムは、必要に応じてＣＰＵ１３に読み出されて実行される。状態推定プログラムＰＲ１が実行されると、図１で述べた状態推定計算部１０１が実現する。計算値ボトルネック指標計算プログラムＰＲ２が実行されると、計測値ボトルネック指標計算部１０３が実現される。計測値組合せ探索プログラムＰＲ３が実行されると、計測値組合せ探索部１０２が実現される。
【００４６】
図４のフローチャートを用いて、状態推定計算の入力値の組合せを求める処理を含む、電力系統の状態推定処理アルゴリズムの一例を説明する。
【００４７】
状態推定装置１は、まずステップＳ１１において、系統条件を設定する。このステップＳ１１は、組合せ探索部１０２で実施される。
【００４８】
ステップＳ１２では、状態推定装置１は、計測値組合せの初期値を設定する。例えば、状態推定装置１は、計測値の最大未使用数の設定や未使用項目の設定を、各計測値の属性データを元に設定する。状態推定装置１は、例えば、ＰＭＵのような高精度なセンサは状態推定計算に必ず用いるように設定し、その他のセンサは可観測性が保たれる範囲（上限値以内）を未使用とできるように設定する。このように、センサの属性に応じて、状態推定計算に用いるセンサの優先順位を複数設けておくことで、後述のセンサ組合せパターン数を低減し、計算時間を短縮することが可能となる。
【００４９】
ステップＳ１３において、状態推定装置１は、ステップＳ１１およびＳ１２で用意したデータを元に状態推定計算を実施し、計算過程で得られる推定残差（計測値と推定値の差ε）を含めて、計算結果を一時的に保存する。この処理は、状態推定計算部１０１で実施される。
【００５０】
ステップＳ１４において、状態推定装置１は、状態推定に用いる計測値の組合せを探索するための初期設定として、カウンタｋを１に初期化する。
【００５１】
ステップＳ１５では、状態推定装置１は、ステップＳ１３で計算された残差εを規定値と比較する。残差εが規定値以下であれば（Ｓ１５：ＹＥＳ）、十分有効な解が得られたと判定して本処理を正常終了する。残差εが規定値以上であれば（Ｓ１５：ＮＯ）、ステップＳ１６へ進む。
【００５２】
ステップＳ１６では、状態推定装置１は、処理の繰り返し回数をカウントするために、カウンタｋを１だけインクリメントする。
【００５３】
ステップＳ１７において、状態推定装置１は、計測値のボトルネック指標△Ｆｘを計算する。このボトルネック指標は、後述のように、計測地点ｍ毎、かつ計測項目ｎ毎に、計算される。
【００５４】
ステップＳ１８において、状態推定装置１は、状態推定計算に用いる計測値の組合せを変更する処理を実施する。この組合せ変更処理は、例えば後述のようにタブサーチ手法のようなヒューリスティック手法を用いてることで行うことができる。または、系統規模が小さい場合は、計測値の組合せを総探索することで解を求めてもよい。
【００５５】
ステップＳ１９では、状態推定装置１は、ステップＳ１８で求めた計測値の組合せに従って状態推定計算を実施し、推定残差を含めて計算結果を一時的に保存する。
【００５６】
ステップＳ２０では、状態推定装置１は、計測値組合せパターンの処理についてのカウント値であるカウンタｋが上限値ｋｍａｘ以上になったか判定する。カウンタｋが上限値ｋｍａｘ以上になった場合（Ｓ２０：ＹＥＳ）、本処理は正常終了する。カウンタｋが上限値ｋｍａｘ未満の場合（Ｓ２０：ＮＯ）、ステップＳ１５に戻り、計測値組合せ探索を繰り返す。処理を終了する場合、状態推定装置１は、計算された解のうち最も良い解を結果データベースＤＢ４へ保存する。この処理は、計算結果記録部１０４で行われる。
【００５７】
図５のフローチャートを用いて、状態推定計算の一例を説明する。ステップＳ３１では、状態推定装置１は、式１で表される観測値Ｚ、系統状態値ｆ（Ｘ）、系統誤差ｅの関係から、式２の値が最小となるように、状態値Ｘを決定する。
【００５８】
Ｚ＝ｆ（Ｘ）＋ｅ・・・（式１）
【００５９】
｛Ｚ−Ｆ（Ｘ）｝２  →  最小・・・（式２）
【００６０】
ここで、Ｚ：観測値（Ｐ、Ｑ、Ｖ）、Ｘ：系統状態（Ｖ、θなど）、Ｆ（Ｘ）：状態値  （回路の接続状態およびインピーダンスから決まる回路方程式）である。
【００６１】
ステップＳ３２は、状態推定装置１は、残差εを計算し、誤りデータを削除する。状態推定装置１は、例えば、残差が大きい計測値を除去したり、代わりの擬似計測値（例えば過去の計測値データ）に置き換えたりして、誤りデータを除去する。
【００６２】
ステップ３３は、計算された残差εを規定値と比較する。残差εが規定値以下であれば本処理を正常終了し、残差εが規定値以上であれば（Ｓ３３：ＮＯ）、ステップＳ３１へ戻る。
【００６３】
なお、このような電力系統の状態推定計算は、非特許文献１にも示されているように確立された計算手法であり、一般的なアルゴリズムを用いることで計算可能である。
【００６４】
図６のフローチャートを用いて、図１のボトルネック指標計算部１０３や図４のステップＳ１７の一例を示す処理概要を説明する。
【００６５】
ステップＳ４１において、状態推定装置１は、状態推定計算を実施し、計算過程で得られる推定残差を含めて計算結果を一時的に保存する。
【００６６】
ステップＳ４２において、状態推定装置１は、探索カウンタである計測点ｍおよびその計測項目ｎの値を初期化する（それぞれ１に設定する）。
【００６７】
ステップＳ４３において、状態推定装置１は、計測点の計測項目を状態推定の入力データとして使用するか否かを決定する。具体的には、状態推定装置１は、計測点ｍの項目ｎが「使用」の場合は「未使用」に変更し、また計測点ｍの項目ｎが「未使用」の場合は「使用」に変更する。
【００６８】
ステップＳ４４では、状態推定装置１は、ステップＳ４３で設定した計測項目の組合せを用いて、Ｓ４１と同様に状態推定計算を実施し、計算過程で得られる推定残差を含めて計算結果を一時的に保存する。
【００６９】
ステップＳ４５において、状態推定装置１は、計測値ボトルネック指標△Ｆｍｎを式３に従って計算する。
【００７０】
△Ｆｍｎ ＝ ε０ − εｍｎ・・・（式３）
【００７１】
ステップＳ４６では、状態推定装置１は、Ｓ４３で実施された計測点ｍの項目ｎの「使用」または「未使用」の設定変更を元に戻す。
【００７２】
ステップＳ４７では、状態推定装置１は、計測点毎の計測項目カウンタｎが上限ｎｍａｘに達したか判定する。上限に達しない場合（Ｓ４７：ＮＯ）、Ｓ４３に戻る。
【００７３】
ステップＳ４８では、状態推定装置１は、計測点カウンタｍが上限ｍｍａｘに達したか判定する。上限に達しない場合（Ｓ４８：ＮＯ）、Ｓ４３へ戻る。状態推定装置１は、すべての計測点についてボトルネック指標を計算すると、本処理を正常終了する。
【００７４】
図７のフローチャートを用いて、図１の組合せ探索部１０２や図４のステップＳ１８の計測値組合せ変更（タブサーチ） の処理概要を説明する。
【００７５】
ステップＳ５１では、状態推定装置１は、タブサーチの探索回数を計測するカウンタｋを１に初期化する。
【００７６】
ステップＳ５２において、状態推定装置１は、計測値ボトルネック指標△Ｆｍｎを計算する。この処理は、前述した図６の処理フローで実施される。
【００７７】
ステップＳ５３では、状態推定装置１は、評価指標Ｆ（＝Σ△Ｆｍｎ）の合計値が最大となるように、タブリストに無い計測値（計測点ｍの項目ｎ）の使用・未使用設定を変更する。
【００７８】
ステップＳ５４において、状態推定装置１はタブリストを更新する。
【００７９】
ステップＳ５５では、状態推定装置１は、状態推定計算を実施し、計算過程で得られる推定残差を含めて計算結果を一時的に保存する。
【００８０】
ステップＳ５６では、状態推定装置１は、ステップＳ５５で計算された残差εをこれまでの探索結果と比較する。状態推定装置１は、最新の残差εがこれまで探索された状態推定結果の最小値よりも小さければ（Ｓ５６：ＹＥＳ）、ステップＳ５７へ進み、残差εおよび計測値を最適組合せとして記録し、それ以外の場合（Ｓ５６：ＮＯ）、ステップＳ５８へ進む。
【００８１】
ステップＳ５８では、状態推定装置１は、探索回数のカウンタｋが上限ｋｍａｘに達したか判定する。上限に達しない場合（Ｓ５８：ＮＯ）、Ｓ５１へ戻る。カウンタｋが上限ｋｍａｘに達した場合（Ｓ５８：ＹＥＳ）、本処理は正常終了する。
【００８２】
図８に示すテーブルを用いて、提案手法の有効性を模式例によって説明する。図示のテーブルは、各センサの各計測項目（例えば、有効電力Ｐ、無効電力Ｑ、電圧Ｖなど）のボトルネック指標の数値と、評価指標残差の総計Σεとを表す。
【００８３】
ボトルネック指標の数値は、「＋」の場合はセンサの使用・未使用を変更した方が残差が小さくなり（結果が改善される）、「−」の場合は変更した方が残差が大きくなる（結果が悪化する）ことを示す。数値の変わりに、「×」が記載されたところは、状態推定が収束しないためにボトルネック指標が計算できないことを示す。 黒く塗られた場所は、その繰り返し回数の状態推定計算に該当する計測値を使用していないことを示す。
【００８４】
各行は、計測点・計測項目の組合せを変更しながら、状態推定計算を繰り返す、探索ステップを表している。今回の例では、この探索を１７回繰り返した例を示している。
【００８５】
計測点・計測項目の組合せを変更する対象は、通常精度のセンサから選択し、高精度のセンサ（ＰＭＵなど）は対象外としている（常時利用するように設定している）。
【００８６】
各繰り返し回数において、ボトルネック指標の数値が＋側に最も大きい項目の、使用・未使用の設定を変更することで、計測項目の組合せを変更している。ここで、タブリスト長は３としている。すなわち、一旦使用・未使用の設定が変更された計測項目は、３回探索が繰り返されるまでは設定変更が禁止される（タブサーチの処理手順）。
【００８７】
この事例では、繰り返し回数９回目が、状態推定評価指標である残差合計値が最も小さく、最良解であることがわかる。また、探索回数によっては状態推定が収束しない回もある（「非収束」と表示されている回）。
【００８８】
このように、状態推定が収束しないケースがある場合でも、演算に用いる計測点・計測項目の組合せを変更することで、収束する解を得る可能性を高めるとともに、残差のより少ない解（もっともらしい解）を求めることが可能となる。
【００８９】
このように構成される本実施例によれば、複数の計測値の中から状態推定計算が収束すると予測される計測値の組合せを選択して、電力系統の状態を推定するため、状態推定装置１の信頼性や使い勝手を向上させることができる。
【００９０】
本実施例によれば、計測値の組合せを最適化することで残差を小さくでき、より信頼性の高い推定が可能となる。
【００９１】
本実施例によれば、残差の大きい不良データを検出できるため、残差からセンサの不具合を遠隔診断することもできる。
【００９２】
本実施例によれば、高精度なセンサ２０５Ａの計測値はすべて使用し、それ以外のセンサ２０５の計測値は精度に応じた所定範囲内で使用するため、効率的にかつ高い信頼性で電力系統の状態を推定することができる。
【００９３】
本実施例では、例えば図８に示すテーブルに含まれる情報の全部または一部を表示装置１１に表示させることもできるため、ユーザは、状態推定計算の詳細を容易に確認することができ、使い勝手が向上する。
【００９４】
なお、本発明は、上述した実施形態に限定されない。当業者であれば、本発明の範囲内で、種々の追加や変更等を行うことができる。上述の実施形態において、添付図面に図示した構成例に限定されない。本発明の目的を達成する範囲内で、実施形態の構成や処理方法は適宜変更することが可能である。
【００９５】
また、本発明の各構成要素は、任意に取捨選択することができ、取捨選択した構成を具備する発明も本発明に含まれる。さらに特許請求の範囲に記載された構成は、特許請求の範囲で明示している組合せ以外にも組み合わせることができる。
【符号の説明】
【００９６】
１：電力系統の状態推定装置、１０１：状態推定計算部、１０２：計測値組合せ探索部、１０３：計測値ボトルネック指標計算部
【図１】
2018098986.tif 000003
【図２】
2018098986.tif 000004
【図３】
2018098986.tif 000005
【図４】
2018098986.tif 000006
【図５】
2018098986.tif 000007
【図６】
2018098986.tif 000008
【図７】
2018098986.tif 000009
【図８】
2018098986.tif 000010