(57)【要約】【課題】状態推定演算結果の異常を検出することができる電力系統監視制御システムを得ることを目的とする。【解決手段】電力系統監視制御システム１００の監視制御装置１は、観測値収集装置１１が収集した観測値に基づいた観測状態の要素値とその推定値とを評価関数にて評価した結果に基づいて電力系統１０の最も確からしい状態を推定する状態推定処理部２と、状態推定処理部２が推定した推定状態を判定する異常検出部３と、異常検出部３の判定結果に基づいて電力系統１０を制御する制御装置８に推定状態を出力する情報伝送部５とを有する。異常検出部３は、評価関数に基づいた判定式にて演算した演算結果を判定基準と比較し、異常であるか否かを判定する。情報伝送部５は、正常と判定された場合に推定状態を制御装置８に出力し、異常と判定された場合に推定状態を制御装置８に出力しない。【選択図】図２
2011024286.tif 000002
【特許請求の範囲】
【請求項１】
電力系統を観測する観測値収集装置と、前記観測値収集装置が収集した観測値を監視する監視制御装置を備えた電力系統監視制御システムであって、
前記監視制御装置は、前記観測値に基づいた観測状態の要素値とその推定値とを評価関数にて評価した結果に基づいて前記電力系統の最も確からしい状態を推定する状態推定処理部と、前記状態推定処理部が推定した推定状態を判定する異常検出部と、前記異常検出部の判定結果に基づいて前記電力系統を制御する制御装置に前記推定状態を出力する情報伝送部とを有し、
前記異常検出部は、前記評価関数に基づいた判定式にて演算した演算結果を判定基準と比較し、前記判定式の演算結果が前記判定基準を超えた場合に異常であると判定し、
前記情報伝送部は、前記異常検出部で正常と判定された場合に前記推定状態を前記制御装置に出力し、前記異常検出部で異常と判定された場合に前記推定状態を前記制御装置に出力しないことを特徴とした電力系統監視制御システム。
【請求項２】
前記判定式は前記評価関数であり、前記判定基準は所定の閾値であり、前記判定式の演算結果が前記閾値よりも大きい場合に異常であると判定することを特徴とした請求項１記載の電力系統監視制御システム。
【請求項３】
前記判定式は前記評価関数を前記電力系統の総需要で除したものであり、前記判定基準は所定の閾値であり、前記判定式の演算結果が前記閾値よりも大きい場合に異常であると判定することを特徴とした請求項１記載の電力系統監視制御システム。
【請求項４】
前記判定式は前記評価関数における今回の演算結果と以前の演算結果との変化を求めるものであり、前記判定基準は所定の上限閾値及び下限閾値の範囲であり、前記判定式の演算結果が前記上限閾値よりも大きい場合または前記下限閾値よりも小さい場合に異常であると判定することを特徴とした請求項１記載の電力系統監視制御システム。
【請求項５】
前記監視制御装置は、前記状態推定処理部にて推定状態が推定された際の前記評価関数の演算結果である評価値を記憶する評価値保存部と、前記評価値保存部に記憶された複数の前記評価値に基づいて新判定基準を生成し、前記判定基準を前記新判定基準に変更する判定基準生成部を有することを特徴とした請求項１乃至４のいずれか１項に記載の電力系統監視制御システム。
【請求項６】
電力系統の観測値に基づいて前記電力系統を監視制御する電力系統監視制御システムの制御方法であって、
前記観測値に基づいた観測状態の要素値とその推定値とを評価関数にて評価した結果に基づいて前記電力系統の最も確からしい状態を推定する推定手順と、前記評価関数に基づいた判定式にて演算した演算結果を判定基準と比較し、前記判定式の演算結果が前記判定基準を超えた場合に異常であると判定する異常検出手順と、前記異常検出手順で正常と判定した場合に前記電力系統を制御する制御装置に前記推定手順で推定した推定状態を出力する出力手順と、前記異常検出手順で異常と判定した場合に前記制御装置に前記推定手順で推定した推定状態を出力しない出力中止手順とを有することを特徴とした電力系統監視制御システムの制御方法。
【請求項７】
前記判定式は前記評価関数であり、前記判定基準は所定の閾値であり、前記判定式の演算結果が前記閾値よりも大きい場合に異常であると判定することを特徴とした請求項６記載の電力系統監視制御システムの制御方法。
【請求項８】
前記判定式は前記評価関数を前記電力系統の総需要で除したものであり、前記判定基準
は所定の閾値であり、前記判定式の演算結果が前記閾値よりも大きい場合に異常であると判定することを特徴とした請求項６記載の電力系統監視制御システムの制御方法。
【請求項９】
前記判定式は前記評価関数における今回の演算結果と以前の演算結果との変化を求めるものであり、前記判定基準は所定の上限閾値及び下限閾値の範囲であり、前記判定式の演算結果が前記上限閾値よりも大きい場合または前記下限閾値よりも小さい場合に異常であると判定することを特徴とした請求項６記載の電力系統監視制御システムの制御方法。
【請求項１０】
前記推定手順にて推定状態が推定された際の前記評価関数の演算結果である評価値を記憶する評価値保存手順と、前記評価値保存手順に記憶された複数の前記評価値に基づいて新判定基準を生成し、前記判定基準を前記新判定基準に変更する判定基準生成手順とを有することを特徴とした請求項６乃至９のいずれか１項に記載の電力系統監視制御システムの制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、電力系統の監視制御システムに関するものである。
【背景技術】
【０００２】
従来の電力系統監視制御システムは、監視制御対象の電力系統のノード（母線）やブランチ（送電線や変圧器）における有効電力や無効電力などの観測値を収集して、当該電力系統の最も確からしい状態を求めている。最も確からしい状態を求める手法として、状態推定方法と潮流計算法がある。特許文献１の電力系統監視制御システムは、状態推定方法で得られた推定値を潮流計算法の指定値として用いて、最も確からしい状態を求めていた。
【０００３】
算出された最も確からしい電力系統の状態は、電力系統監視制御システムに実装されている需給自動制御や電圧自動制御の初期値として利用されており、状態推定結果が異常な値となった場合は、電力系統の運用上の安全性、経済性が損なわれることになり、大きな問題となる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−９２６８５号公報（図２）
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来の電力系統監視制御システムは、状態推定方法と潮流計算法を組み合わせたとしても、最初に求める状態推定方法で数値演算においてかろうじて収束した場合は、推定された電力系統の状態が監視制御対象の電力系統の状態と不整合になり、需給自動制御や電力自動制御に影響を及ぼすという問題点があった。
【０００６】
この発明は、状態推定演算結果の異常を検出することができる電力系統監視制御システムを得ることを目的とする。
【課題を解決するための手段】
【０００７】
電力系統監視制御システムの監視制御装置は、観測値収集装置が収集した観測値に基づいた観測状態の要素値とその推定値とを評価関数にて評価した結果に基づいて電力系統の最も確からしい状態を推定する状態推定処理部と、状態推定処理部が推定した推定状態を判定する異常検出部と、異常検出部の判定結果に基づいて電力系統を制御する制御装置に推定状態を出力する情報伝送部とを有する。異常検出部は、評価関数に基づいた判定式にて演算した演算結果を判定基準と比較し、異常であるか否かを判定する。情報伝送部は、正常と判定された場合に推定状態を制御装置に出力し、異常と判定された場合に推定状態を制御装置に出力しない。
【発明の効果】
【０００８】
この発明に係る電力系統監視制御システムは、状態推定演算結果の異常を検出することができ、電力系統への運用上の影響を事前に防止することができる。
【図面の簡単な説明】
【０００９】
【図１】この発明の実施の形態１における電力系統監視制御システムを示す図である。
【図２】図１の監視制御装置を示す図である。
【図３】図１の監視制御装置のフローチャートである。
【図４】図１の電力系統の総需要と評価値を示す図である。
【図５】図１の電力系統の評価値の判定条件を説明する図である。
【図６】この発明の実施の形態２における電力系統の評価値の判定条件を説明する図である。
【図７】この発明の実施の形態３における電力系統の評価値の判定条件を説明する図である。
【図８】実施の形態３における監視制御装置の他のフローチャートである。
【図９】この発明の実施の形態４における監視制御装置を示す図である。
【発明を実施するための形態】
【００１０】
実施の形態１．
図１はこの発明の実施の形態１における電力系統監視制御システムを示す図であり、図２は監視制御装置を示す図である。電力系統監視制御システム１００は複数の観測値収集装置１１ａ、１１ｂ、１１ｃ、監視制御装置１、ネットワーク１３を有している。観測値収集装置１１ａ、１１ｂ、１１ｃは監視制御対象である電力系統１０のノード（母線）やブランチ（送電線や変圧器）における観測値を収集する。監視制御装置１は状態推定処理部２、異常検出部３、警報生成部４、情報伝送部５を有する。制御装置８は、電力系統１０の各種の機器である制御対象装置１２に対して需給自動制御や電圧自動制御などの制御を行う。
【００１１】
状態推定処理部２は複数の観測値収集装置１１から伝送された観測値から電力系統１０の最も確からしい状態を求める。観測値は電圧Ｖ、電流Ｉ、有効電力Ｐ、無効電力Ｑなどである。電力系統１０の状態は観測箇所のデータの集合であり、例えば電圧Ｖや位相θの集合である。位相θは有効電力Ｐに基づいて求めることができる。異常検出部３は状態推定処理部２が演算した推定値が異常かどうかを判定する。情報伝送部５は、異常検出部３で推定値が異常でない（正常である）と判定された場合に制御装置８に状態推定結果である推定値を伝送する。警報生成部４は、異常検出部３で推定値が異常であると判定された場合に警報を生成する。警報はブザーの鳴動や表示装置の点灯や点滅などである。
【００１２】
状態推定処理部２が行う電力系統１０の状態推定演算について説明する。状態推定演算は評価関数Ｊ（ｘ）を最小化するように演算を行う。評価関数Ｊ（ｘ）は式（１）の通りである。
Ｊ（ｘ）＝（１／２）×[ｚ−ｈ（ｘ）]ＴＷ[ｚ−ｈ（ｘ）]      ・・・（１）
ここで、観測値ベクトルｚは電力系統１０の複数の観測値の１次元ベクトルであり、推定値ベクトルｈ（ｘ）は電力系統１０の複数の観測箇所の推定値を要素とする１次元ベクトルである。係数行列Ｗは各計測値に対する重み係数を要素とした正方行列である。
【００１３】
観測値ベクトルｚは、観測値収集装置１１が測定し収集したデータである電圧Ｖ、電流Ｉ、有効電力Ｐ、無効電力Ｑや位相θを用いることができる。ここでは電圧Ｖと位相θを例として説明する。ｘは時間と関連付けされた数字であり、それ故に観測値ベクトルｚは、特定の時間ｔ１におけるｎ個の観測箇所の電圧Ｖ及び位相θを要素としている。具体的には、２ｎ行１列の行列形式で表され、１行からｎ行にはＶ１からＶｎが配置され、ｎ＋１行から２ｎ行にはθ１からθｎが配置される。ｈ（ｘ）は観測値ベクトルｚと同じ形式に配置された電圧Ｖ及び位相θの推定値を要素としている。[ｚ−ｈ（ｘ）]Ｔは[ｚ−ｈ（ｘ）]の転置行列であり、評価関数Ｊ（ｘ）は観測値ベクトルｚの各要素とこれに対応する推定値ベクトルｈ（ｘ）の各要素との差、すなわち誤差の２乗を重み付けした和となっている。なお、この場合における係数行列Ｗは、ｎ行ｎ列の行列である。
【００１４】
図３は電力系統１０の総需要と評価値を示す図である。図３（ａ）は電力系統１０における電気使用量の総和あるは供給電力量の総和である総需要Ｄの時間推移であり、図３（ｂ）は評価関数Ｊ（ｘ）の演算結果である評価値Ｊが最小であると推定した確定評価値Ｊｏの時間推移である。図３（ａ）において、縦軸は総需要Ｄであり、横軸は時間Ｔである。図３（ｂ）において、縦軸は確定評価値Ｊｏであり、横軸は時間Ｔである。図４は確定評価値Ｊｏと総需要Ｄの時間を細分化した拡大図である。横軸は時間Ｔであり、左側の縦軸は総需要Ｄであり、右側の縦軸は確定評価値Ｊｏである。三角印を有する線は総需要Ｄの時間推移を示し、丸印を有する線は確定評価値Ｊｏの時間推移を示している。状態推定演算は、一般的にある一定の周期（１分周期や３分周期など）毎に実行している。三角印や丸印は、例えば３分毎の値である。時間ｔ２において確定評価値Ｊｏ（以降Ｊｅとする）は他の値に比べて著しく異なる数値を示している。評価値Ｊｅは異常値である。
【００１５】
図４に示したＡ１は異常検出部３にて確定評価値Ｊｏが正常か異常かを判定する閾値である。判定基準である閾値Ａ１は、状態推定処理部２で演算した推定値が異常であったときの評価値Ｊｅと推定値が正常であったときの確定評価値Ｊｏを統計的に処理し、異常と判断すべき値を算出して求めた値である。図４において、閾値Ａ１を１．０にすることで、評価値Ｊｅを異常値であると判定することができる。
【００１６】
次に監視制御装置１の動作を図５のフローチャートを用いて説明する。図５は監視制御装置１のフローチャートである。監視制御装置１は複数の観測値収集装置１１から送られてくる観測値を受け取り、所定の数の要素数に達した場合に、観測値ベクトルｚを生成する。状態推定処理部２にて観測値ベクトルｚに対して推定値ベクトルｈ（ｘ）の初期値により評価関数Ｊ（ｘ）を演算する。推定値ベクトルｈ（ｘ）を順次変更しながら評価関数Ｊ（ｘ）が最小になるように演算を繰り返し、最も確からしい推定値ベクトルｈ（ｘ）を求める（ステップＳ１０１）。ステップＳ１０１の状態推定演算が収束したかどうかを判定する（ステップＳ１０２）。
【００１７】
ステップＳ１０２で収束しなかった場合は、終了する。ステップＳ１０２で収束したと判定した場合は、最も確からしいと推定した確定推定値ベクトルｈｏに対する評価関数Ｊ（ｘ）の結果である確定評価値Ｊｏを異常検出部３で判定する（ステップＳ１０３）。なお、ステップＳ１０３における判定式は評価関数Ｊ（ｘ）である。異常検出部３において、確定評価値Ｊｏが評価値Ｊの閾値Ａ１より大きい（Ｊｏ＞Ａ１）場合に異常であると判定する。異常ではなく正常と判定した場合には、状態推定結果である確定推定値ベクトルｈｏを、情報伝送部５を介して制御装置８に出力する（ステップＳ１０４）。
【００１８】
ステップＳ１０３にて、異常であると判定した場合には、状態推定結果である確定推定値ベクトルｈｏを、情報伝送部５を介して制御装置８に出力することを中止する（ステップＳ１０５）。その後、警報生成部４において、警報鳴動や警報表示等の警報を生成する（ステップＳ１０６）。
【００１９】
以上のように実施の形態１の電力系統監視制御システム１００は状態推定演算結果の異常を、状態推定演算の度にその直後において検出することができ、異常と判定した場合には制御装置８に状態推定演算結果を出力しないので、電力系統１０への運用上の影響を事前に防止することができる。電力系統監視制御システム１００は、電力系統１０の観測値の異常（例えば観測値収集装置１１の故障）により計算が収束しない場合も制御装置８に状態推定演算結果を出力しないので、電力系統１０への運用上の影響を事前に防止することができる。
【００２０】
従来は、数値演算においてかろうじて収束した場合は、その状態推定演算結果が正常であるか異常であるか判定できなかったので、その後の処理を行ってしまい、推定された電
力系統の状態が監視制御対象の電力系統の状態と不整合になり、需給自動制御や電力自動制御に影響を及ぼしていた。また、フェールセーフとして状態推定演算結果の収束の程度によって電力系統１０の監視制御を停止する場合には、オペレータ（運用者）によって判断されて再開するまでに多大な時間を要し、早急に電力系統１０を制御すべき事象に適切に対応できない問題があった。
【００２１】
しかしながら、実施の形態１の電力系統監視制御システム１００は、異常と判定された状態推定演算結果を制御装置８に状態推定演算結果を出力しないだけではなく、次の観測値を受け取り、状態推定演算を行うので、従来とは異なり、早急に電力系統１０を制御すべき事象に適切に対応することができる。
【００２２】
以上のように実施の形態１の電力系統監視制御システム１００によれば、監視制御装置１は、観測値収集装置１１が収集した観測値に基づいた観測状態ｚの要素値とその推定値とを評価関数Ｊ（ｘ）にて評価した結果に基づいて電力系統１０の最も確からしい状態ｈｏを推定する状態推定処理部２と、状態推定処理部２が推定した推定状態ｈｏを判定する異常検出部３と、異常検出部３の判定結果に基づいて電力系統１０を制御する制御装置８に推定状態ｈｏを出力する情報伝送部５とを有し、異常検出部３は、評価関数Ｊ（ｘ）に基づいた判定式Ｊ（ｘ）にて演算した演算結果を判定基準Ａ１と比較し、異常であるか否かを判定し、情報伝送部５は、正常と判定された場合に推定状態ｈｏを制御装置８に出力し、異常と判定された場合に推定状態ｈｏを制御装置８に出力しないので、状態推定演算結果の異常を検出することができ、電力系統への運用上の影響を事前に防止することができる。
【００２３】
実施の形態１の電力系統監視制御システム１００の制御方法によれば、観測値に基づいた観測状態ｚの要素値とその推定値とを評価関数Ｊ（ｘ）にて評価した結果に基づいて電力系統１０の最も確からしい状態ｈｏを推定する推定手順と、評価関数Ｊ（ｘ）に基づいた判定式Ｊ（ｘ）にて演算した演算結果を判定基準Ａ１と比較し、判定式Ｊ（ｘ）の演算結果が判定基準Ａ１を超えた場合に異常であると判定する異常検出手順と、異常検出手順で正常と判定した場合に電力系統１０を制御する制御装置８に推定手順で推定した推定状態ｈｏを出力する出力手順と、異常検出手順で異常と判定した場合に制御装置に推定手順で推定した推定状態ｈｏを出力しない出力中止手順とを有するので、状態推定演算結果の異常を検出することができ、電力系統への運用上の影響を事前に防止することができる。
【００２４】
なお、判定における閾値として１つを用いる場合で説明したが、時間帯を区切りその時間帯毎の閾値を用いてもよい。この場合は時間帯毎に最適な閾値を使用することができ、高精度に状態推定演算結果の異常を検出することができる。
【００２５】
実施の形態２．
実施の形態１では、異常検出部３における状態推定演算結果の判定を評価関数Ｊ（ｘ）の確定評価値Ｊｏが評価値Ｊの閾値Ａ１より大きい（Ｊｏ＞Ａ１）場合に異常であると判定していたが、判定式Ｊ（ｘ）／Ｄ（ｘ）の結果を判定することで精度の高い判定を行うことができる。判定式Ｊ（ｘ）／Ｄ（ｘ）は評価関数Ｊ（ｘ）を総需要Ｄで平準化したものである。以下に説明する。
【００２６】
図６はこの発明の実施の形態２における電力系統の評価値の判定条件を説明する図であり、確定評価値Ｊｏを総需要Ｄで除したＪｏ／Ｄの時間推移を示したものである。縦軸はＪｏ／Ｄであり、横軸は時間Ｔである。図６は図４と同一の時間におけるデータを用いたものである。
【００２７】
図４に示したように確定評価値Ｊｏは異常値を除くと総需要Ｄに比例している。これは
評価関数値が主として観測値収集装置１１のそれぞれの観測値と推定値の誤差の総和であるため、観測値の絶対量が多い場合、すなわち電力系統の総需要値が多い場合は、相対的に評価関数値も大きくなるためである。
【００２８】
例えば、電気の使用量が多い日中の確定評価値Ｊｏは、電気の使用量が少ない夜間の確定評価値Ｊｏよりも大きくなる。そのため、実施の形態１の方法では、閾値Ａ１の値が固定である場合は、閾値Ａ１は夜間における判定としては余裕が大きすぎて、電力系統を高精度に制御する上で改善する余地があった。
【００２９】
状態推定処理部２における状態推定演算は評価関数Ｊ（ｘ）を最小化するように演算を行う。異常検出部３の異常判定においては、評価関数Ｊ（ｘ）の演算結果とそのときの総需要Ｄ（ｘ）を受け、判定式Ｊ（ｘ）／Ｄ（ｘ）の結果を判定することで、総需要Ｄの変動に影響されること無く、高精度に異常を検出できる。図６に示したＡ２は異常検出部３にて確定評価値Ｊｏが正常か異常かを判定する閾値である。閾値Ａ２は、状態推定処理部２で演算した確定評価値Ｊｏが異常であったときの判定値Ｊｅ／Ｄｅと推定値が正常であったときの確定評価値Ｊｏ／Ｄを統計的に処理し、異常と判断すべき値を算出して求めた値である。図６において、判定基準である閾値Ａ２を１．０ｅ−４にすることで、判定値Ｊｏ／Ｄを異常値であると判定することができる。なお、Ｄｅは確定評価値Ｊｏが異常であったときの総需要である。
【００３０】
実施の形態２の電力系統監視制御システム１００は、総需要Ｄの変動による影響を極力排除し、実施の形態１に比べて高精度に異常を検出できる。したがって、実施の形態１に比べて電力系統１０を高精度に制御することができる。
【００３１】
なお、状態推定処理部２における状態推定演算は評価関数Ｊ（ｘ）を最小化するように演算を行う場合で説明したが、状態推定処理部２においてＪ（ｘ）／Ｄ（ｘ）を評価関数として使用し、Ｊ（ｘ）／Ｄ（ｘ）を最小化するように演算しても構わない。
【００３２】
実施の形態３．
実施の形態１では、異常検出部３における状態推定演算結果の判定を評価関数Ｊ（ｘ）の確定評価値Ｊｏが評価値Ｊの閾値Ａ１より大きい（Ｊｏ＞Ａ１）場合に異常であると判定していたが、確定評価値Ｊｏの時系列を用いても異常であると判定することができる。以下に説明する。
【００３３】
図７はこの発明の実施の形態３における電力系統の評価値の判定条件を説明する図であり、確定評価値Ｊｏの時系列の比Ｊｏ（−１）／Ｊｏの時間推移を示したものである。確定評価値Ｊｏ（−１）は確定評価値Ｊｏの一つまえのものである。縦軸はＪｏ（−１）／Ｊｏであり、横軸は時間Ｔである。図７は図４と同一の時間におけるデータを用いたものである。Ａ３及びＡ４は異常検出部３にて確定評価値Ｊｏが正常か異常かを判定する閾値である。異常検出部３は前回の確定評価値Ｊｏ（−１）を記憶する記憶部を有する。
【００３４】
状態推定演算は、一般的にある一定の周期（１分周期や３分周期など）毎に実行しているので、この程度の間隔であれば、総需要Ｄが大きく変化することはないため、前回の演算で推定した確定評価値Ｊｏと比較することにより異常を検出することができる。確定評価値Ｊｏの時系列の比Ｊｏ（−１）／Ｊｏを判定式として異常判定を行う。異常判定は、判定基準である下限閾値Ａ３と上限閾値Ａ４を用いて式（２）により判定する。
Ａ３≦Ｊｏ（−１）／Ｊｏ≦Ａ４      ・・・（２）
確定評価値Ｊｏの時系列の比Ｊｏ（−１）／Ｊｏが式（２）を満足する場合には、今回の演算結果である確定評価値Ｊｏは異常でない（正常である）と判定する。確定評価値Ｊｏの時系列の比Ｊｏ（−１）／Ｊｏが式（２）を満足しない場合には、今回の演算結果であ
る確定評価値Ｊｏは異常であると判定する。
【００３５】
状態推定処理部２における状態推定演算は評価関数Ｊ（ｘ）を最小化するように演算を行う。異常検出部３の異常判定においては、評価関数Ｊ（ｘ）の演算結果を受け、記憶していた前回の評価関数Ｊ（ｘ）の演算結果の比、すなわちＪｏ（−１）／Ｊｏを判定することで、総需要Ｄの変動に影響されること無く、高精度に異常を検出できる。閾値Ａ３及びＡ４は、状態推定処理部２で演算した確定評価値Ｊｏが異常であったときの判定値Ｊｅ（−１）／Ｊｅと推定値が正常であったときの確定評価値Ｊｏ（−１）／Ｊｏを統計的に処理し、異常と判断すべき値を算出して求めた値である。図６において、閾値Ａ３を０．９にし、閾値Ａ４を１．１にすることで、判定値Ｊｏ（−１）／Ｊｏを異常値であると判定することができる。
【００３６】
実施の形態３の電力系統監視制御システム１００は、総需要Ｄの変動による影響を極力排除し、実施の形態１に比べて高精度に異常を検出できる。したがって、実施の形態１に比べて電力系統１０を高精度に制御することができる。
【００３７】
ところで、電力系統１０の運用において影響が大きい事象として、停電等を伴う系統事故もある。停電を伴う大規模な事故が発生した場合は、電気の使用量が一時的に大きく減少し、暫く減少しているため、結果として状態推定演算の評価関数値も大きく減少し、暫くの間減少する。
【００３８】
したがって、状態推定演算の異常検出に前回の確定評価値Ｊｏ（−１）を用いることにより、状態推定演算結果の異常の検出だけでなく電力系統事故の際にも異常として検出することができる。異常と判定した場合には制御装置８に状態推定演算結果を出力しないので、電力系統１０への運用上の影響を事前に防止することができる。
【００３９】
なお、状態推定演算結果の異常と電力系統事故を区別して制御するには、異常検出部３の記憶部は異常と判定した確定評価値Ｊｏ、すなわち評価値Ｊｅを記憶し、Ｊｏ（−１）／Ｊｏに代えてＪｅ／Ｊｏを判定式にして事故かどうかを判定するようにする。監視制御装置１は図８のフローチャートのように動作する。図８は監視制御装置の他のフローチャートであり、図５のフローチャートとはステップＳ１０７、Ｓ１０８及びＳ１０９が追加された点で異なる。ステップＳ１０３にて異常検出部３において、確定評価値Ｊｏが異常と判定された場合はステップＳ１０７に進み、事故かどうかを判定する。ステップＳ１０３で異常と判定された評価値Ｊｅとの比、すなわちＪｅ／Ｊｏを判定式として判定する。
Ａ３≦Ｊｅ／Ｊｏ≦Ａ４      ・・・（３）
Ｊｅ／Ｊｏが式（３）を満足する場合には、今回の演算結果である確定評価値Ｊｏは事故でないと判定する。Ｊｅ／Ｊｏが式（３）を満足しない場合には、今回の演算結果である確定評価値Ｊｏは事故であると判定する。なお、ステップＳ１０３で異常と判定された直後は事故ではないと判定される。
【００４０】
ステップＳ１０７にて、事故ではないと判定した場合にはステップＳ１０５に進む。ステップＳ１０７にて事故であると判定した場合には、監視制御装置１は情報伝送部５を介して制御対象装置１２に自動制御の一時停止を指令する。
【００４１】
ステップＳ１０３で異常と判定された後に、ステップＳ１０３にて正常と判定された場合はステップＳ１０９へ進む。ステップＳ１０９では過去ｍ回前、例えば１回前にステップＳ１０３で異常と判定されたかを判定する。正しくない場合はステップＳ１０４に進む。ステップＳ１０９の判定が正しい場合はステップＳ１０７に進む。
【００４２】
事故が発生した場合には状態推定演算の評価関数値が暫く大きく減少する。したがって
複数回ステップＳ１０７にて判定を行うことで事故であると判定することができる。
【００４３】
制御対象装置１２は自動制御の一時停止の指令を受けると手動制御に切替える。このようにすることで電力系統１０の制御対象装置１２を一時停止させ、手動制御に切替えることができる。
【００４４】
電力系統１０の制御対象装置１２は、コスト抑制が図られており、限定された自動制御機能の動作をしている場合には、大規模な事故が発生した場合など緊急時は、オペレータ（運用者）による手動制御を優先すべきである。この場合には、図８のフローチャートのように制御することで、電力系統１０の制御対象装置１２を一時停止させ、手動制御に切替えることができる。
【００４５】
なお、実施の形態３では状態推定処理部２における状態推定演算は評価関数Ｊ（ｘ）を最小化し、確定評価値Ｊｏを異常検出部３に出力する場合で説明したが、状態推定処理部２において前回の確定評価値Ｊｏ（−１）を記憶しておき、Ｊｏ（−１）／Ｊｏを異常検出部３に出力しても構わない。この場合、状態推定処理部２は前回の確定評価値Ｊｏ（−１）を記憶する記憶部を有する。異常検出部３には前回の確定評価値Ｊｏ（−１）を記憶する記憶部は必要ない。
【００４６】
また、確定評価値Ｊｏの時系列の比Ｊｏ（−１）／Ｊｏで判定する場合で説明したが、確定評価値Ｊｏの時系列の変化を判定するには、他の式を用いても構わない。例えば以下のような式Ｆを用いることができる。
Ｆ＝（Ｊｏ−Ｊｏ（−１））／Ｊｏ      ・・・（４）
Ｆ＝（Ｊｏ−Ｊｏ（−１））／Ｊｏ（−１）      ・・・（５）
Ｆ＝ＡＢＳ（Ｊｏ−Ｊｏ（−１））      ・・・（６）
ここで、ＡＢＳは絶対値を求める関数である。なお、異常を判定する判定式が式（４）乃至（６）のいずれかが用いられた場合は、ステップＳ１０７での判定式は異常を判定する判定式におけるＪｏ（−１）をＪｅに代えた式を用いる。また、ステップＳ１０３の判定式とステップＳ１０７の判定式が異なっても構わない。事故判定する場合の閾値を、状態推定演算結果の異常判定する場合の閾値を異なるものにしても構わない。
【００４７】
実施の形態４．
実施の形態４では評価関数Ｊ（ｘ）の閾値あるいは評価関数Ｊ（ｘ）を含む他の評価関数の閾値を計算し、フィードバックするようにした。以下に説明する。図９はこの発明の実施の形態４における監視制御装置を示す図である。図１の監視制御装置とは評価値保存部６及び判定基準生成部７を有する点で異なる。以下実施の形態１の確定評価値Ｊｏと閾値Ａ１で判定する場合で説明する。
【００４８】
評価値保存部６は評価値保存手順を実行し、定基準生成部７は判定基準生成手順を実行する。評価値保存部６は評価関数Ｊ（ｘ）の推定した結果である確定評価値Ｊｏを記憶する。判定基準生成部７は確定評価値Ｊｏを統計的に処理し、推定値が正常であったときの平均値と標準偏差σに基づいて閾値Ａ１と比較を行う。例えば３σの上限値と閾値Ａ１を比較し、余裕が大きい場合には閾値Ａ１の値を小さい値に変更する。３σの上限値と閾値Ａ１を比較し、余裕が小さい場合には閾値Ａ１の値を大きい値に変更する。閾値Ａ１の値を大きい値に変更する場合は、推定値が異常であったときの平均値と標準偏差σとも比較して、推定値が異常の３σの下限値に近づきすぎないようにする。
【００４９】
実施の形態４の電力系統監視制御システム１００は、評価関数Ｊ（ｘ）の演算結果を記憶し、統計処理を行って異常判定に用いる閾値を変更することができる。したがって異常判定に用いる閾値を適切な値に更新するので、電力系統監視制御システム１００を安定し
て運用することができる。すなわち、閾値Ａ１が厳しすぎた場合に自動制御用のデータとして確定評価値Ｊｏが伝送されないために生じる急激な変更制御を防止し、緩やかな変更制御にすることができる。また、閾値Ａ１が緩すぎた場合に自動制御用のデータとして確定評価値Ｊｏが伝送されたために生じる不適正な変更制御を防止できる。
【００５０】
なお、実施の形態４の内容は実施の形態２乃至４にも適用できる。
【００５１】
なお、実施の形態１乃至４では、観測値ベクトルｚとしてノードやブランチ毎の値を用いる場合で説明したが、ノードやブランチに複数のセンサを設けて複数の観測値を収集する場合には、複数の観測値の平均を用いて当該ノードやブランチの観測値とすることで確定推定値ベクトルｈｏの精度を向上させることができる。観測値収集装置１１で処理して平均値を監視制御装置１に送信してもよいし、監視制御装置１で平均値の処理をしてもよい。
【００５２】
実施の形態２では、評価関数Ｊ（ｘ）を平準化するのに総需要Ｄを用いた例で説明したが、推定対象のノードやブランチの状態数Ｎで評価関数Ｊ（ｘ）を平準化することもできる。総需要Ｄが大きくなると必然的に推定対象のノードやブランチの状態数Ｎは大きくなるので、総需要Ｄは状態数Ｎにある程度比例する。したがって、Ｊ（ｘ）／Ｎを閾値で判定しても実施の形態２の効果を奏する。また、推定対象のノードやブランチの観測値の数Ｍで評価関数Ｊ（ｘ）を平準化することもできる。ノードやブランチの観測値の数Ｍはノードやブランチの状態数Ｎと等しい又は多い。総需要Ｄが大きくなるに従って観測値の数Ｍは大きくなるので、総需要Ｄは観測値の数Ｍにある程度比例する。したがってＪ（ｘ）／Ｍを閾値で判定しても実施の形態２の効果を奏する。
【産業上の利用可能性】
【００５３】
この発明に係る電力系統監視制御システムは、電力系統の監視及び制御を好適に実行することができる。
【符号の説明】
【００５４】
８      監視制御装置                  ２      状態推定処理部
３      異常検出部                    ５      情報伝送部
６      評価値保存部                  ７      判定基準生成部
８      制御装置                      １０    電力系統
１１    観測値収集装置                １１ａ  観測値収集装置
１１ｂ  観測値収集装置                １１ｃ  観測値収集装置
１００  電力系統監視制御システム      ｚ      観測値ベクトル
ｈｏ    確定推定値ベクトル            Ｊ（ｘ）  評価関数
Ｊ（ｘ）／Ｄ（ｘ）  判定式            Ｊｏ（−１）／Ｊｏ  判定式
Ａ１    閾値                          Ａ２    閾値
Ａ３    下限閾値                      Ａ４    上限閾値
【図１】
2011024286.tif 000003
【図２】
2011024286.tif 000004
【図３】
2011024286.tif 000005
【図４】
2011024286.tif 000006
【図５】
2011024286.tif 000007
【図６】
2011024286.tif 000008
【図７】
2011024286.tif 000009
【図８】
2011024286.tif 000010
【図９】
2011024286.tif 000011