JPS61273130A - 電力系統の状態推定方式 - Google Patents
電力系統の状態推定方式Info
- Publication number
- JPS61273130A JPS61273130A JP60111837A JP11183785A JPS61273130A JP S61273130 A JPS61273130 A JP S61273130A JP 60111837 A JP60111837 A JP 60111837A JP 11183785 A JP11183785 A JP 11183785A JP S61273130 A JPS61273130 A JP S61273130A
- Authority
- JP
- Japan
- Prior art keywords
- power system
- estimated
- value
- measurement
- measurement error
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
この発明は電力系統の状態推定方式に関するものである
。
。
第１図は、例えばリライアプル パット データ プロ
セシング フォー リアル−タイム ステート エステ
イメイショｙ（Ｒｅｌｉａｂｌｅ Ｂａｄ Ｄａｔａ
Ｐｒｏｃｅｓｓｉｎｇ ｆｏｒＲｅａｌ −Ｔｉ
ｍｅ ５ｔａｔｅ Ｅｓｔｉｍａｔｉｏｎ）ＩＥＥＥ
ＰＡ８−１０２、ム５．ＰＰ１１２６−１１３９に示
された従来の状態推定方式のフローチャートである。図
において５Ｔ−１は電力系統の状態量ベクトルＸ（母線
の電圧絶対値および電圧位相角で構成されるベクトル）
、ならびに推定残差ベクトル（実計測値ベクトルＺと推
定される計測値ベクトルｈ（３）との差）に対する重み
行列Ｗの初期値を与えるＸ、Ｗの初期設定ブロック、５
Ｔ−２は推定残差ベクトルのＷによる重み付け二乗和（
評価関数Ｊ（３）】を最小とする電力系統の状態量の推
定値Ｘを求めるブロック、５Ｔ−３は推定残差ベクトル
の各要素の大きさよシ実計測値の中に、計測誤差が大き
いと疑われる計測値が存在するか否かを判定するブロッ
クで疑わしい計測値が存在しない場合には処理を終了し
、存在する場合には、５Ｔ−４で示すブロックへ進む。
セシング フォー リアル−タイム ステート エステ
イメイショｙ（Ｒｅｌｉａｂｌｅ Ｂａｄ Ｄａｔａ
Ｐｒｏｃｅｓｓｉｎｇ ｆｏｒＲｅａｌ −Ｔｉ
ｍｅ ５ｔａｔｅ Ｅｓｔｉｍａｔｉｏｎ）ＩＥＥＥ
ＰＡ８−１０２、ム５．ＰＰ１１２６−１１３９に示
された従来の状態推定方式のフローチャートである。図
において５Ｔ−１は電力系統の状態量ベクトルＸ（母線
の電圧絶対値および電圧位相角で構成されるベクトル）
、ならびに推定残差ベクトル（実計測値ベクトルＺと推
定される計測値ベクトルｈ（３）との差）に対する重み
行列Ｗの初期値を与えるＸ、Ｗの初期設定ブロック、５
Ｔ−２は推定残差ベクトルのＷによる重み付け二乗和（
評価関数Ｊ（３）】を最小とする電力系統の状態量の推
定値Ｘを求めるブロック、５Ｔ−３は推定残差ベクトル
の各要素の大きさよシ実計測値の中に、計測誤差が大き
いと疑われる計測値が存在するか否かを判定するブロッ
クで疑わしい計測値が存在しない場合には処理を終了し
、存在する場合には、５Ｔ−４で示すブロックへ進む。
８Ｔ−４は、疑わしい計測値の中から最も計測誤差が大
きいと疑わしい計測値を選択し、当該計測値に対応する
Ｗの列をすてる（零とする）ブロックで、この変更され
た重み行列Ｗによシ再び８Ｔ−２で示すブロックにて電
力系統の状態量ベクトルＸが推定される。
きいと疑わしい計測値を選択し、当該計測値に対応する
Ｗの列をすてる（零とする）ブロックで、この変更され
た重み行列Ｗによシ再び８Ｔ−２で示すブロックにて電
力系統の状態量ベクトルＸが推定される。
従来の電力系統の状態推定方式においては、計測誤差が
大きいと疑われる計測値が存在する場合には第１図の５
Ｔ−２で示す評価関数Ｊ（Ｘ）を特徴とする請求めるブ
ロックの処理を複数回行う必要があったため処理時間が
長くなるという問題点があった。
大きいと疑われる計測値が存在する場合には第１図の５
Ｔ−２で示す評価関数Ｊ（Ｘ）を特徴とする請求めるブ
ロックの処理を複数回行う必要があったため処理時間が
長くなるという問題点があった。
この発明は上記のような問題点を解消するためになされ
たもので、第１図の５Ｔ−２の処理を１同案行するのみ
で計測誤差の大きい計測値を検出し、処理時間を大巾に
短縮できるようにすることを目的とする。
たもので、第１図の５Ｔ−２の処理を１同案行するのみ
で計測誤差の大きい計測値を検出し、処理時間を大巾に
短縮できるようにすることを目的とする。
この発明に係る電力系統の状態推定方式は１推定残差を
用いて、重み付け最小二乗法によシ１実計測値に含まれ
る計測誤差を推定し、計測誤差の大きい不良計測データ
を検出するとともに、推定された計測誤差を用いて、電
力系統の状態量に関する推定値を補正しよシ良い電力系
統の状態量に関する推定値を得ることによシ処理を複数
回行うことを不要とし、処理の高速化を達成したもので
ある。
用いて、重み付け最小二乗法によシ１実計測値に含まれ
る計測誤差を推定し、計測誤差の大きい不良計測データ
を検出するとともに、推定された計測誤差を用いて、電
力系統の状態量に関する推定値を補正しよシ良い電力系
統の状態量に関する推定値を得ることによシ処理を複数
回行うことを不要とし、処理の高速化を達成したもので
ある。
この発明における状態推定方式は計測誤差が大きいと疑
われる計測値を選択し、当該計測値の推定残差に対する
感度行列を先ず求め、続いて、評価関数の最小値を求め
、計測誤差の推定値で補正してから状態の最適の推定値
を求める。
われる計測値を選択し、当該計測値の推定残差に対する
感度行列を先ず求め、続いて、評価関数の最小値を求め
、計測誤差の推定値で補正してから状態の最適の推定値
を求める。
以下、この発明の一実施例を図について説明する。
図中、第２図と同一の部分は同一の符号をもって図示し
た第２図において、５Ｔ−４０は計測誤差が大きいと疑
われる計測値を選択し、当該計測値の推定残差に対する
感度行列Ｓを作成するブロック、５Ｔ−５０は、推定残
差几を用いて重み付け最小二乗法によシ、実計測値に含
まれる計測誤差ΔＺを推定するつまシ、評価関数Ｊｏ（
Δ２）を最小とするΔ２を求めるブロック、５Ｔ−６０
は、推定されたΔ２によシ、電力系統の状態量に関すツ
クである。
た第２図において、５Ｔ−４０は計測誤差が大きいと疑
われる計測値を選択し、当該計測値の推定残差に対する
感度行列Ｓを作成するブロック、５Ｔ−５０は、推定残
差几を用いて重み付け最小二乗法によシ、実計測値に含
まれる計測誤差ΔＺを推定するつまシ、評価関数Ｊｏ（
Δ２）を最小とするΔ２を求めるブロック、５Ｔ−６０
は、推定されたΔ２によシ、電力系統の状態量に関すツ
クである。
以下、本発明の主要部分に関し詳細な説明を行う。
（１）感度行倒Ｓ
まず、計測誤差が大きいと疑われる計測値の、推定残差
に対する感度行列Ｓは、下記（１）式の行列の、計測誤
差が大きいと疑われる計測値に対応する列のみをぬき出
したものである。
に対する感度行列Ｓは、下記（１）式の行列の、計測誤
差が大きいと疑われる計測値に対応する列のみをぬき出
したものである。
Ｉ −Ｈ（Ｈ’ Ｗ Ｈ）−’ ｆ（’ｗ ＋ｊ
＋＋＋＋＋＋＋＋ （１）ただし、Ｉ：、行×ｍ列の単
位行列 ｍは計測値の数 θＮ ＨＩｄ ｍ行×ｎ列の行列 ｎは電力系統の状態量の数 Ｗ：＠行×ｍ列の重み行列 （２）評価関数Ｊｏ（Δ２） また、状態推定の本来の評価関数Ｊ□Ｑは、最適な推定
値であるＸに対しても、零とはならないことが一般的で
ある。これは、計測値に誤差が含まれるため、次の（２
）式で示す推定残差Ｒが０とはならないためである。
＋＋＋＋＋＋＋＋ （１）ただし、Ｉ：、行×ｍ列の単
位行列 ｍは計測値の数 θＮ ＨＩｄ ｍ行×ｎ列の行列 ｎは電力系統の状態量の数 Ｗ：＠行×ｍ列の重み行列 （２）評価関数Ｊｏ（Δ２） また、状態推定の本来の評価関数Ｊ□Ｑは、最適な推定
値であるＸに対しても、零とはならないことが一般的で
ある。これは、計測値に誤差が含まれるため、次の（２
）式で示す推定残差Ｒが０とはならないためである。
几＝ｚ−ｈ（ｘ） ・・・・・・・
・・ （２）本発明では、前記（２）式の几が０となら
ないのは、主として計測誤差が大きいと疑われる小数の
計測値の持つ計測誤差Δ２に起因していると考え、新し
い評価関数Ｊｏ（Δ２）を（３）式を用いて導入する。
・・ （２）本発明では、前記（２）式の几が０となら
ないのは、主として計測誤差が大きいと疑われる小数の
計測値の持つ計測誤差Δ２に起因していると考え、新し
い評価関数Ｊｏ（Δ２）を（３）式を用いて導入する。
（３）式（Ｑ Ｊｏ （Δ２）は、本来の評価関数を計
測誤差Δ２の関数と見たものであシ１本発明ではＪｏ（
Δりを最小とするΔ子を計測誤差の推定値とする。なお
、ＪＯ（Δ多）に最小値を与えるΔ２は（４）式によシ
容易に求めることができる。
測誤差Δ２の関数と見たものであシ１本発明ではＪｏ（
Δりを最小とするΔ子を計測誤差の推定値とする。なお
、ＪＯ（Δ多）に最小値を与えるΔ２は（４）式によシ
容易に求めることができる。
（３）Ｘの補正
計測値に含まれる計測誤差の推定値Δ多が求まったこと
によシ最適の推定値Ｘを計測誤差の行列を導入し線形補
正を（５）式を用いて行う。
によシ最適の推定値Ｘを計測誤差の行列を導入し線形補
正を（５）式を用いて行う。
なお、上記の説明では、電力系統の状態推定を例にとっ
て説明したが、本発明は重み付け最小二乗法によシ状態
推定を行うあらゆる分野に適用可能である。
て説明したが、本発明は重み付け最小二乗法によシ状態
推定を行うあらゆる分野に適用可能である。
以上のように１この発明によれば、電力系統の状態推定
処理の中で大きな時間をしめる状態推定値Ｘを求める演
算処理を、１同案行するのみでよいように演算による推
定方式を与えたので、処理時間が大幅に短縮できる効果
がある。
処理の中で大きな時間をしめる状態推定値Ｘを求める演
算処理を、１同案行するのみでよいように演算による推
定方式を与えたので、処理時間が大幅に短縮できる効果
がある。
また、この発明によれば、複数の不良計測データがある
場合にも、計測誤差が大きいと疑われる計測値を複数一
括して扱い計測誤差を推定するので、明確な不良データ
の検出が可能となる効果がある。
場合にも、計測誤差が大きいと疑われる計測値を複数一
括して扱い計測誤差を推定するので、明確な不良データ
の検出が可能となる効果がある。
第１図は本発明の一実施例による電力系統の状態推定方
式を示すフローチャート、第２図は従来の電力系統の状
態推定方式のフローチャートである。 図において、 ５Ｔ−４０：感度行列Ｓを求めるブロック５Ｔ−５０：
評価関数を最小とする計測誤差の推定値を求めるブロッ
ク ５Ｔ−６０：最適の推定値を計測誤差の推定値によシ補
正し、よシ良い推定値を求め るブロック 手続補正書（自発）
式を示すフローチャート、第２図は従来の電力系統の状
態推定方式のフローチャートである。 図において、 ５Ｔ−４０：感度行列Ｓを求めるブロック５Ｔ−５０：
評価関数を最小とする計測誤差の推定値を求めるブロッ
ク ５Ｔ−６０：最適の推定値を計測誤差の推定値によシ補
正し、よシ良い推定値を求め るブロック 手続補正書（自発）
Claims (3)
- （１）推定残差ベクトルによる重み付け最小二乗法を用
いる電力系統の状態推定方式において、前記電力系統の
状態量の実計測値及び推定される計測値との推定残差を
用い、前記重みと同一の重みによる重み付け最小二乗法
により実計測に含まれる計測誤差を感度行列を用いて線
形補正し、状態推定するようにしたことを特徴とする電
力系統の状態推定方式。 - （２）前記状態推定方式の評価関数式を計測誤差の関数
として捕え、該評価関数を最小とする計測誤差の推定値
を求めるようにしたことを特徴とする特許請求の範囲第
１項記載の電力系統の状態推定方式。 - （３）前記実計測値に含まれる計測誤差の推定値を用い
電力系統の状態量に関する推定値を補正するようにした
ことを特徴とする特許請求の範囲第１項記載の電力系統
の状態推定方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60111837A JPS61273130A (ja) | 1985-05-24 | 1985-05-24 | 電力系統の状態推定方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60111837A JPS61273130A (ja) | 1985-05-24 | 1985-05-24 | 電力系統の状態推定方式 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS61273130A true JPS61273130A (ja) | 1986-12-03 |
Family
ID=14571412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60111837A Pending JPS61273130A (ja) | 1985-05-24 | 1985-05-24 | 電力系統の状態推定方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS61273130A (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6373831A (ja) * | 1986-09-16 | 1988-04-04 | 株式会社日立製作所 | 電力系統の状態推定方法 |
| JP2011024286A (ja) * | 2009-07-13 | 2011-02-03 | Mitsubishi Electric Corp | 電力系統監視制御システム及び制御方法 |
| JP2013208051A (ja) * | 2012-03-27 | 2013-10-07 | Mitsubishi Electric Corp | 電力系統のバスの状態を求める方法及びシステム |
| JP2013219902A (ja) * | 2012-04-06 | 2013-10-24 | Fuji Electric Co Ltd | 電力系統の状態推定方法 |
| CN105375464A (zh) * | 2015-10-30 | 2016-03-02 | 国家电网公司 | 一种计及智能可控设备的暂态安全稳定紧急控制方法 |
| JP2017034739A (ja) * | 2015-07-28 | 2017-02-09 | 株式会社ダイヘン | 系統インピーダンス推定装置、インバータ装置、および、系統インピーダンス推定方法 |
-
1985
- 1985-05-24 JP JP60111837A patent/JPS61273130A/ja active Pending
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6373831A (ja) * | 1986-09-16 | 1988-04-04 | 株式会社日立製作所 | 電力系統の状態推定方法 |
| JP2011024286A (ja) * | 2009-07-13 | 2011-02-03 | Mitsubishi Electric Corp | 電力系統監視制御システム及び制御方法 |
| JP2013208051A (ja) * | 2012-03-27 | 2013-10-07 | Mitsubishi Electric Corp | 電力系統のバスの状態を求める方法及びシステム |
| JP2013219902A (ja) * | 2012-04-06 | 2013-10-24 | Fuji Electric Co Ltd | 電力系統の状態推定方法 |
| JP2017034739A (ja) * | 2015-07-28 | 2017-02-09 | 株式会社ダイヘン | 系統インピーダンス推定装置、インバータ装置、および、系統インピーダンス推定方法 |
| CN105375464A (zh) * | 2015-10-30 | 2016-03-02 | 国家电网公司 | 一种计及智能可控设备的暂态安全稳定紧急控制方法 |
| CN105375464B (zh) * | 2015-10-30 | 2018-06-05 | 国家电网公司 | 一种计及智能可控设备的暂态安全稳定紧急控制方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0078592A2 (en) | Method of and apparatus for automatically compensating for variations in output response characteristics of sensors and the like | |
| EP2982994A1 (en) | Interference compensating single point detecting current sensor for a multiplex busbar | |
| US4520321A (en) | Phase difference detector with phase inversion of one input signal when phase difference is small | |
| JPS61273130A (ja) | 電力系統の状態推定方式 | |
| CN108593203A (zh) | 基于8051的压力传感器校准算法 | |
| CN103605904B (zh) | 基于误差估算的自补偿电力系统幅值算法 | |
| JP3676486B2 (ja) | 地絡検出装置 | |
| JPS62218813A (ja) | 圧力検出装置 | |
| CN114061524A (zh) | 一种钢卷轮廓测量方法及装置 | |
| EP3705902A1 (en) | Method of determining an absolute angle of a magnetic field | |
| Hirokawa et al. | Isotachophoretic determination of mobility and pKa by means of computer simulation: II. Evaluation of mo and pKa of 65 anions | |
| WO2022030148A1 (ja) | 測定装置及びその測定方法並びに測定プログラム | |
| JPH11142465A (ja) | 地絡故障点検出方法 | |
| JPH04109174A (ja) | インピーダンス測定装置 | |
| JPS628721B2 (ja) | ||
| JP2002296209A (ja) | 複合溶液の濃度計測システム | |
| JP3265681B2 (ja) | 三相交流測定における補正方法 | |
| JPS61231837A (ja) | 電力系統の状態推定方式 | |
| JP3318649B2 (ja) | 測定データ処理装置および方法 | |
| JP3099542B2 (ja) | 交流入力振巾のデジタル推定演算方法 | |
| JP7056363B2 (ja) | 電池システム | |
| JPH0744781B2 (ja) | 電力系統の状態推定に用いる重みの調整方式 | |
| CN115752456A (zh) | 一种惯组产品温补参数校正方法及系统 | |
| JPS60146143A (ja) | 湿度計 | |
| JPH0572000B2 (ja) |