(57)【要約】【課題】設備のパラメータ設定値に誤りがあった場合、その誤りデータを検出し、推定精度を向上する。【解決手段】少なくとも電力系統を構成する各設備１１０の有効電力、無効電力、および端子電圧を含む計測値、並びに送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切情報、および変圧器のタップ位置情報を含む設備のパラメータを含む状態量に基づいて、系統を構成する母線あるいは送電線における状態の推定値を算出する電力系統の状態推定方法において、前記各設備の時系列データから得られる設備特性係数マトリクスを用いて各計測値の大きさのばらつきを分析して誤りデータを検出する。【選択図】図１
2008253076.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも電力系統を構成する各設備の有効電力、無効電力、および端子電圧を含む計測値、並びに送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切情報、および変圧器のタップ位置情報を含む設備のパラメータを含む状態量に基づいて、系統を構成する母線あるいは送電線における状態の推定値を算出する電力系統の状態推定方法において、
前記各設備の時系列データから得られる設備特性係数マトリクスを用いて各計測値の大きさのばらつきを分析して誤りデータを検出することを特徴とする電力系統の状態推定方法。
【請求項２】
少なくとも電力系統を構成する各設備の有効電力、無効電力、および端子電圧を含む計測値、並びに送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切情報、および変圧器のタップ位置情報を含む設備のパラメータを含む状態量に基づいて、系統を構成する母線あるいは送電線における状態の推定値を算出する電力系統の状態推定方法において、
前記各設備の時系列データから得られる設備特性係数マトリクスを用いて各計測値の大きさのばらつきを分析して誤りデータを検出するとともに検出した誤りデータを除去した時系列データを用いて状態値の推定値を再計算することにより推定精度を向上することを特徴とする電力系統の状態推定方法。
【請求項３】
少なくとも電力系統を構成する各設備の有効電力、無効電力、および端子電圧を含む計測値、並びに送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切情報、および変圧器のタップ位置情報を含む設備のパラメータを含む状態量に基づいて、系統を構成する母線あるいは送電線における状態の推定値を算出する電力系統の状態推定装置において、
前記状態の計測値および各設備のパラメータを収集して記憶する記憶部と、該記憶部に記憶されたデータを用いて前記状態の推定値を演算し、該演算の結果と前記計測値の偏差をもとに誤りデータを検出し、検出した誤りデータを除去したデータを用いて推定値を再度演算する誤りデータ判定部を備え、前記状態の推定値を再計算することにより推定精度を向上することを特徴とする電力系統の状態推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統の状態推定方法に係り、特に、設備のパラメータ設定値に誤りがあった場合、その誤りデータを検出し、推定精度をより向上することのできる電力系統の状態推定技術に関する。
【背景技術】
【０００２】
電力系統において、電力潮流を解析する際には、互いに独立変数である、例えば各母線の電圧とその位相角を「状態」として解析する。電力系統における情報収集は、数秒の間隔で、大量のデータを計測することになるため、収集情報にはある程度の計測誤差あるいは不等時性が含まれてくる。このため、これらの測定データをそのまま使用したのでは、キルヒホッフの法則を満足できないことがあるため、状態推定法が考えられた。
【０００３】
状態推定法によれば、特許文献１に述べられているように、電力系統の、ある設備において、最小限の計測値が分かれば、キルヒホッフの法則により、それに関連する値を知ることができる。しかし計測値が正確でなく、ある程度の測定誤差があるときには、単に前記最小限の測定値のみならず、数多くの測定値を利用した方が、より良い状態推定を行うことができる。このため、従来の電力系統状態推定法では、１つの計測ポイントの推定値の精度を良くする為に、その計測ポイントの周辺設備のデータを集めて、これに最小２乗法あるいは重み付けなどの経時要素を伴わない断面的な計算手法を適用している。
【特許文献１】特開平１１−１１３１７４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
電力系統に変更が加えられ、その変更が監視制御システムのデータベースに適切に反映されなかった場合、あるいは系統を構成する機器、電子装置等の誤操作あるいは誤測定により設備のパラメータ設定値に誤りがある場合、すなわち、インピーダンス値あるいはタップ比に誤りがある場合に、従来の状態推定法を適用すると、計算に用いる周辺の設備からの有効注入電力測定値、無効注入電力測定値、電圧測定値の整合性を整えるため、周りの設備の推定値が調整される。このため、パラメータデータの誤りがある設備を検出することができない。
【０００５】
本発明は、このような問題点に鑑みてなされたもので、設備のパラメータ設定値に誤りがあった場合、その誤りデータを検出し、推定精度をより向上することのできる電力系統の状態推定技術を提供する。
【課題を解決するための手段】
【０００６】
本発明は上記課題を解決するため、次のような手段を採用した。
【０００７】
少なくとも電力系統を構成する各設備の有効電力、無効電力、および端子電圧を含む計測値、並びに送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切情報、および変圧器のタップ位置情報を含む設備のパラメータを含む状態量に基づいて、系統を構成する母線あるいは送電線における状態の推定値を算出する電力系統の状態推定方法において、前記各設備の時系列データから得られる設備特性係数マトリクスを用いて各計測値の大きさのばらつきを分析して誤りデータを検出する。
【発明の効果】
【０００８】
本発明は、以上の構成を備えるため、設備のパラメータ設定値に誤りがあった場合、その誤りデータを検出し、推定精度をより向上することができる。
【発明を実施するための最良の形態】
【０００９】
以下、最良の実施形態を添付図面を参照しながら説明する。図１は、本発明の実施形態にかかる電力系統状態推定装置を説明する図である。図１において、１は電力系統状態推定装置、１ないしＮは電力系統を構成する設備機器に設定したサイトであり、それぞれ設備１１０、並びに前記設備の状態（電圧、電流、位相等）およびパラメータ（線路のインピーダンス、変圧器のタップ位置等）を計測する計測器１０９を備える。１０８は電力系統状態推定装置１０１と各サイト１ないしＮを接続する広域ネットワークである。
【００１０】
１０６は通信部、１０７は記憶部、１０２は情報収集部である。情報収集部１０２は通信部を介して収集したデータを記憶部１０７に格納する。１０３は潮流計算部であり、記憶部に格納されたデータをもとに潮流を計算する。１０４は誤りデータ判定部であり、記憶部に格納されたデータをもとに系統状態の推定値を算出し、算出した推定値と計測値（あるいは予め計測して格納しておいた設定値）の偏差（あるいはばらつき）が所定値以上であるとき前記計測値を誤りータとして検出する。１０５は出力部であり、誤りデータ判定部が判定した誤りデータを表示装置に出力する。
【００１１】
設備１ないしＮの状態を表すデータ、すなわち、有効電力、無効電力、および端子電圧を含む計測値、送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切状態、および変圧器のタップ位置を表すデータは計測器１０９によって計測され、得られたデータは、広域ネットワーク１０８を介して電力系統状態推定装置１０１に送信される。
【００１２】
電力系統状態推定装置１０１は、通信部１０６を介して前記データを受け取り、記憶部１０７に格納する。潮流計算部１０３は、前記受信したデータと、例えば予め測定し、記憶部に格納してある送電線のインピーダンス、変圧器のインピーダンス等の設備パラメータをもとに電力の潮流を計算する。また、誤りデータ判定部１０４は、前記受信したデータと、例えば予め測定し記憶部に格納してある送電線のインピーダンス、変圧器のインピーダンス等の設備パラメータをもとに系統を構成する母線あるいは送電線における状態の推定値を算出する。
【００１３】
図２は、誤りデータ判定部の処理を説明する図である。まず、電力系統から系統を構成する各設備の有効電力、無効電力、および端子電圧等を含む計測値、並びに送電線のインピーダンス、変圧器のインピーダンス、開閉器の入・切情報、および変圧器のタップ位置情報等を含む設備のパラメータを含む情報を収集する（ステップＳＴ１１）。次に、計算量を低減するため、系統の解析対象部から見て遠方の系統を等価的に縮約表現する（ステップＳＴ１２）。次に、縮約された電力系統データを用いて、電力系統の状態を表す関数行列式（ヤコビアン）を作成する（ステップＳＴ１３）。次に、作成した行列式を用いて、電力系統の状態を推定演算する（ステップＳＴ１４）。
【００１４】
推定演算結果は、バッドデータ判定部に伝送される。バッドデータ判定部は演算結果と計測値とを比較し、計測値と演算結果の差が予め設定した閾値以上である場合に、前記計測値を異常データ（バッドデータ）として除去する（ステップＳＴ１５）。また、演算結果はパラメータ判定部に伝送される。パラメータ誤り判定部は図３に示す方法で取得した推定値をヤコビアン作成部に戻し、ヤコビアン作成部は再度推定演算を行う。また、前記計測値と演算結果の差が予め設定した閾値未満である場合には異常なしと判定し、データを出力部に出力する（ステップＳＴ１７）。
【００１５】
図３は、パラメータ誤り判定部の処理を説明する図である。なお、ここでは公知のカルマンフィルタを用いた処理（ＳＴ２３〜ＳＴ２７）によりパラメータの推定値を演算している。
【００１６】
まず、前記ステップＳＴ１４における演算結果である状態推定結果を受け取り（ステップＳＴ２１）、系統機器の接続状態（接続されているか否か）を確認する。接続状態に異常がある場合には再度確認を行う（ステップＳＴ２２）。
【００１７】
接続状態が正常であれば、予め設定した計測値（計測値の推定値）ｘ０（数式および図面中ではアッパーラインで示しているが明細書中ではアンダーラインで示す）および初期誤差値Ｍ０を初期値として与える（ステップＳＴ３）。次に、与えられた計測値ｘ０および初期誤差値Ｍ０をもとに設備特性係数マトリクスＰ０，Ｇ０，Ｚ０を求める（ステップＳＴ２４）。次に、前記求めた設備特性係数マトリクスＰ０，Ｇ０，Ｚ０を用いて、より精度の高い計測値（推定値）ｘ１および誤差値Ｍ１を求める（ステップＳＴ２５）。次に前記計測値（推定値）ｘ１および誤差値Ｍ１をもとにより精度の高い設備特性マトリクスＰ１，Ｇ１，Ｚ１を求める（ステップＳＴ２４）。次に、前記求めた設備特性係数マトリクスＰ１，Ｇ１，Ｚ１を用いて、より精度の高い計測値（推定値）ｘ２および誤差値Ｍ２を求める（ステップＳＴ２５）。この計算（ステップＳＴ２３，２４，２５，２７）を繰り返すことにより、より精度の高い推定値を得ることができる。
【００１８】
次に、推定値（例えばｘ１）を出力し、出力した推定値と計測値（実際値）を比較し、偏差が閾値以上の場合にはパラメータ異常と判定し、前記推定値を計測値として出力する（ステップＳＴ２６，２８，２９）。
【００１９】
推定値を求める計算では、時刻ｋにおいては設備の真の状態値はｘｋ、推定値はｘｋとする。真値ｘｋと推定値ｘｋの差の二乗の総和を最小にするような係数行列を求める。
【数１】
2008253076.tif 000003
【数２】
2008253076.tif 000004
【数３】
2008253076.tif 000005
【００２０】
ここで、ｙｋは観測値、Ｇｋ−１,Ｑｍは係数行列、Ｚは補正項である。
【００２１】
ここで求めるべきＱｍは無限個あるが、特性係数マトリックスによって現在の推定値は１つ前の推定値に既知の定行列をかけ、現在の観測値で決まる補正項を加えることによって求められる。
【００２２】
特性係数マトリックスは変動方程式と観測方程式からなる。変動方程式はある設備の状態を式で表したものである。例えば、設備状態を考えると未来時刻の設備状態は過去のデータから変動方程式を通して推測できる。観測方程式は、観測値と関連する別の数値である。
【００２３】
ここで、変動方程式は次式で表すことができる。
【数４】
2008253076.tif 000006
【００２４】
ここで、ｘｋは状態変数、Ａｋは状態遷移行列、Ｂｋは変動行列、ｕｋは入力である。
【００２５】
また、 観測方程式は次式で表すことができる。
【数５】
2008253076.tif 000007
【００２６】
ここで、ｙｋは観測値、Ｃｋは観測行列、ｗｋはノイズである。
【００２７】
これらの式を用いて推定値ｘｋを求める。
【数６】
2008253076.tif 000008
【数７】
2008253076.tif 000009
【００２８】
ｘｋは、計測したい設備の状態値と仮定する。Ｍｋはその誤差である。Ａ，Ｂ，Ｃは変動方程式と観測方程式の係数行列で、Ｕ，Ｗは本発明で用いた設備のデータ変動率とデータそのものの標準偏差である。Ａ，Ｂ，Ｃ，Ｕ，Ｗはすでに分かっているものである。Ｐｋはｘｋの共分散である。ここで、ｘ０とＭ０に初期値を与えると、Ｍ０よりＰ０が得られ、Ｐ０よりＧ０が得られ、Ｐ０とｙ０よりＺが得られ、Ｇ０とＺよりｘ１が得られる。次にＰ０よりＭ１が得られる。
【００２９】
次に図３の処理の処理（カルマンフィルタを用いた処理（ステップＳＴ２３〜ＳＴ２７））を説明する。
【００３０】
計測値ｘ０と予め決めた初期誤差値Ｍ０を与える（ステップＳＴ２３）。
【数８】
2008253076.tif 000010
【００３１】
これによって設備特性係数マトリックスＰ０とＧ０とＺ０が得られる（ステップＳＴ２４）。
【数９】
2008253076.tif 000011
【００３２】
次に、Ｐ０とＧ０とＺ０によってｘ１とＭ１が得られる（ステップＳＴ２５）。
【数１０】
2008253076.tif 000012
【００３３】
次に、ｘ１とＭ１を用いて（ステップＳＴ２７）より正確な設備特性係数マトリックスＰ１とＧ１とＺ１を得る。
【数１１】
2008253076.tif 000013
【００３４】
Ｐ１とＧ１とＺ１によってより精度のいいｘ２とＭ２が得られる（ステップＳＴ２３〜２７を繰り返す）。
【００３５】
このように順次計算していくと、ｘｋが求められるたびに係数行列Ｇｋ−１と補正項Ｚｋ−１が更新されるので、精度が徐々に上昇する。
【００３６】
以上説明したように、本実施形態によれば、予測と更新の手順が反復して、計測値と推定値の補正項が常に更新される。このためパラメータデータに誤りがある設備を検出し、推定値を精度よくすることができる。
【図面の簡単な説明】
【００３７】
【図１】実施形態にかかる電力系統状態推定装置を説明する図である。
【図２】誤りデータ判定部の処理を説明する図である。
【図３】パラメータ誤り判定部の処理を説明する図である。
【符号の説明】
【００３８】
１０１  電力系統状態推定装置
１０２  情報収集部
１０３  潮流計算部    
１０４  誤りデータ判定部
１０５  出力部
１０６  通信部
１０７  記憶部
１０８  広域ネットワーク
１０９  計測器
１１０  設備
【図１】
2008253076.tif 000014
【図２】
2008253076.tif 000015
【図３】
2008253076.tif 000016