JPH07114534B2 - 保護継電装置 - Google Patents
保護継電装置Info
- Publication number
- JPH07114534B2 JPH07114534B2 JP63064217A JP6421788A JPH07114534B2 JP H07114534 B2 JPH07114534 B2 JP H07114534B2 JP 63064217 A JP63064217 A JP 63064217A JP 6421788 A JP6421788 A JP 6421788A JP H07114534 B2 JPH07114534 B2 JP H07114534B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- sampling
- substation
- fourier coefficient
- calculated
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/25—Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Emergency Protection Circuit Devices (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 この発明は電力設備の保護継電装置に関するものであ
る。
る。
時刻ｔの関数として変化する標本ｘ（ｔ）が時間周期Ｔ
で周期的であり、この標本ｘ（ｔ）の一周期分の記録が
等間隔に時間ｈだけ離れたＮ点でサンプリングされ、こ
のサンプリング値から任意の時刻ｔにおけるｘ（ｔ）を
算出する方法については、一般に標本化定理として知ら
れている。
で周期的であり、この標本ｘ（ｔ）の一周期分の記録が
等間隔に時間ｈだけ離れたＮ点でサンプリングされ、こ
のサンプリング値から任意の時刻ｔにおけるｘ（ｔ）を
算出する方法については、一般に標本化定理として知ら
れている。
これは式（１）と式（２）に示すものであり、例えば文
献「裳華房出版、佐藤洋著、情報理論、昭和50年８月15
日、P56〜P95」等に等価な式が示されている。
献「裳華房出版、佐藤洋著、情報理論、昭和50年８月15
日、P56〜P95」等に等価な式が示されている。
ただしｔは時刻で、その起点はｎ′＝Ｎサンプリング時
点がｔ＝Ｎ×ｈとなる様にとったもの。
点がｔ＝Ｎ×ｈとなる様にとったもの。
ｎ′はサンプリングの時系列を表わす整数であるが、対
象とする１周期分、すなわちＮ点のデータに対して古い
ものからｎ′＝１、２、…、Ｎとしたもの。
象とする１周期分、すなわちＮ点のデータに対して古い
ものからｎ′＝１、２、…、Ｎとしたもの。
ｑは被サンプリング量の基本周波数の倍調波数を表わす
整数。
整数。
Ｎは１周期Ｔ内のサンプリング点数で、Ｎ＝T/h（ここ
でｈはＴの整数分の１とする）。
でｈはＴの整数分の１とする）。
ω＝２π/T Ａ（ｑ）、Ｂ（ｑ）はＮ点のサンプリングデータに対し
てのｑ倍調波のフーリエ係数である。
てのｑ倍調波のフーリエ係数である。
ここで、保護継電装置のように、サンプリングデータが
１点得られる毎に、最新のＮ点のデータからフーリエ係
数を算出することが望まれる場合には、サンプリングデ
ータが１点得られる毎に対象とするＮ点に対してｎ′＝
１、２、…、Ｎと対応させて式（１）を演算し、最新の
サンプリング時点がｔ＝Ｎ×ｈとなる様に時刻の起点を
考えて式（２）の計算を行い、これによって任意の時刻
ｔにおけるｘ（ｔ）を算出することができる。なお、
Ｎ、ｈの具体数値例としては、50Hzの電力系統であれ
ば、Ｔ＝1/50秒で、Ｎ＝12、ｈ＝1/600秒とする等が考
えられる。
１点得られる毎に、最新のＮ点のデータからフーリエ係
数を算出することが望まれる場合には、サンプリングデ
ータが１点得られる毎に対象とするＮ点に対してｎ′＝
１、２、…、Ｎと対応させて式（１）を演算し、最新の
サンプリング時点がｔ＝Ｎ×ｈとなる様に時刻の起点を
考えて式（２）の計算を行い、これによって任意の時刻
ｔにおけるｘ（ｔ）を算出することができる。なお、
Ｎ、ｈの具体数値例としては、50Hzの電力系統であれ
ば、Ｔ＝1/50秒で、Ｎ＝12、ｈ＝1/600秒とする等が考
えられる。
保護継電装置に式（１）、（２）を応用した例は特願昭
62−276084号にあり、そこでは複数端子における同一時
刻の電流値を求める為のサンプリング同期ずれ時間補正
に用いた実施例が示されている。
62−276084号にあり、そこでは複数端子における同一時
刻の電流値を求める為のサンプリング同期ずれ時間補正
に用いた実施例が示されている。
第５図はフーリエ係数を算出して保護演算を行う従来の
保護継電装置を示すハードウェアブロック図で、被保護
電力設備として送電線を対象としている。
保護継電装置を示すハードウェアブロック図で、被保護
電力設備として送電線を対象としている。
図中、（11）はＡ変電所とＢ変電所を結ぶ送電線、（1
2）はA/D変換部、（７）は高速マイクロプロセッサ、
（14A）と（14B）は光電気変換部、（15）は情報伝送
部、（16）は変復調部、（17）は搬送端局装置、（18）
は無線装置、（19）は全２重回線、（20）はI/O部、（2
1）はしゃ断器である。
2）はA/D変換部、（７）は高速マイクロプロセッサ、
（14A）と（14B）は光電気変換部、（15）は情報伝送
部、（16）は変復調部、（17）は搬送端局装置、（18）
は無線装置、（19）は全２重回線、（20）はI/O部、（2
1）はしゃ断器である。
また、第６図は高速マイクロプロセッサ（７）にて行な
う処理のプログラムフローを示したもので、Ａ変電所で
のサンプリングとＢ変電所でのサンプリングの時間差を
得る処理ステップS31、処理の起動要因判別処理ステッ
プS32、Ｂ変電所でのサンプリング時刻と同一時刻にお
けるＡ変電所での電流値を後述の従来の式に基いて算出
する処理ステップS63、電流値の和が０に近いかどうか
を判定する処理ステップS34、及びトリップ指令を出す
処理ステップS35を有している。
う処理のプログラムフローを示したもので、Ａ変電所で
のサンプリングとＢ変電所でのサンプリングの時間差を
得る処理ステップS31、処理の起動要因判別処理ステッ
プS32、Ｂ変電所でのサンプリング時刻と同一時刻にお
けるＡ変電所での電流値を後述の従来の式に基いて算出
する処理ステップS63、電流値の和が０に近いかどうか
を判定する処理ステップS34、及びトリップ指令を出す
処理ステップS35を有している。
さらに、第７図は処理ステップS63を更に詳細に記した
プログラムフローで、従来の方法によるフーリエ係数の
算出を行う処理ステップS63Aとフーリエ係数をもとに所
望時間の電流値を算出する処理ステップS63Bとを有して
いる。
プログラムフローで、従来の方法によるフーリエ係数の
算出を行う処理ステップS63Aとフーリエ係数をもとに所
望時間の電流値を算出する処理ステップS63Bとを有して
いる。
次に動作について説明する。以下、Ａ変電所側での動作
を中心に説明するが、Ｂ変電所側でも同様の動作であ
る。送電線（11）の電気量である電流をA/D変換部（1
2）によってサンプルホールド後ディジタル値に変換
し、高速マイクロプロセッサ（７）に入力する。一方、
A/D変換部（12）の出力するディジタル値は光電気変換
部（14A）、光電気変換部（14B）、情報伝送部（15）、
変復調部（16）、搬送端局装置（17）、無線装置（18）
及び全２重回線（19）を経由してＢ変電所に送信される
ようになされ、他方Ｂ変電所でサンプリングされた電流
のディジタル値は同様に全２重回線（19）を経由してＡ
変電所に送信され、これを受けてＡ変電所では、このデ
ィジタル値を無線装置（18）、搬送端局装置（17）、変
復調部（16）、情報伝送部（15）、光電気変換部（14
B）及び光電気変換部（14A）を経由してマイクロプロセ
ッサ（７）に入力するようになされている。
を中心に説明するが、Ｂ変電所側でも同様の動作であ
る。送電線（11）の電気量である電流をA/D変換部（1
2）によってサンプルホールド後ディジタル値に変換
し、高速マイクロプロセッサ（７）に入力する。一方、
A/D変換部（12）の出力するディジタル値は光電気変換
部（14A）、光電気変換部（14B）、情報伝送部（15）、
変復調部（16）、搬送端局装置（17）、無線装置（18）
及び全２重回線（19）を経由してＢ変電所に送信される
ようになされ、他方Ｂ変電所でサンプリングされた電流
のディジタル値は同様に全２重回線（19）を経由してＡ
変電所に送信され、これを受けてＡ変電所では、このデ
ィジタル値を無線装置（18）、搬送端局装置（17）、変
復調部（16）、情報伝送部（15）、光電気変換部（14
B）及び光電気変換部（14A）を経由してマイクロプロセ
ッサ（７）に入力するようになされている。
しかして、Ａ変電所におけるサンプリングとＢ変電所に
おけるサンプリングは相手変電所と非同期な自走クロッ
クを基準にして一定時間周期ｈで行われる。
おけるサンプリングは相手変電所と非同期な自走クロッ
クを基準にして一定時間周期ｈで行われる。
高速マイクロプロセッサ（７）では第６図のフローに示
す処理を行う。この処理はＡ変電所でのA/D変換完了時
とＢ変電所からの電流値受信時に行われるものである。
まず処理ステップS31によってＡ変電所でのサンプリン
グとＢ変電所でのサンプリングの時間差を得る。これは
例えばサンプリング周期ｈよりも充分に短い時間周期毎
にカウントアップされるカウンタ値を、Ａ変電所でのA/
D変換完了時とＢ変電所からの電流値受信時にとり込
み、更に、Ｂ変電所からＡ変電所への電流値を伝送する
伝送時間を考慮することによって算出できるものであ
る。Ｂ変電所からＡ変電所への伝送時間は電送路が不変
である限り一定時間であり、例えば理論的に算出してお
くか、又は、一度だけ測定しておき、常時その値を用い
ることで対応できる。
す処理を行う。この処理はＡ変電所でのA/D変換完了時
とＢ変電所からの電流値受信時に行われるものである。
まず処理ステップS31によってＡ変電所でのサンプリン
グとＢ変電所でのサンプリングの時間差を得る。これは
例えばサンプリング周期ｈよりも充分に短い時間周期毎
にカウントアップされるカウンタ値を、Ａ変電所でのA/
D変換完了時とＢ変電所からの電流値受信時にとり込
み、更に、Ｂ変電所からＡ変電所への電流値を伝送する
伝送時間を考慮することによって算出できるものであ
る。Ｂ変電所からＡ変電所への伝送時間は電送路が不変
である限り一定時間であり、例えば理論的に算出してお
くか、又は、一度だけ測定しておき、常時その値を用い
ることで対応できる。
次に、処理ステップS32によってＢ変電所からの電流値
受信時には処理ステップS63に分岐し、そうでない時に
は処理を終了し次の処理開始条件成立を待つ。
受信時には処理ステップS63に分岐し、そうでない時に
は処理を終了し次の処理開始条件成立を待つ。
処理ステップS63ではＢ変電所から着信した電流値をサ
ンプリングしたのと同一時刻におけるＡ変電所での電流
値を算出するが、これはフーリエ係数を算出する処理ス
テップS63Aと、フーリエ係数をもとに所望の時刻の電流
値を算出する処理ステップS63Bとに分けられる。
ンプリングしたのと同一時刻におけるＡ変電所での電流
値を算出するが、これはフーリエ係数を算出する処理ス
テップS63Aと、フーリエ係数をもとに所望の時刻の電流
値を算出する処理ステップS63Bとに分けられる。
処理ステップS63Aでは前述した式（１）の演算を行ない
フーリエ係数Ａ（ｑ）とＢ（ｑ）を算出する。また、処
理ステップS63Bでは、式（２）の演算を行う。
フーリエ係数Ａ（ｑ）とＢ（ｑ）を算出する。また、処
理ステップS63Bでは、式（２）の演算を行う。
ただし、Ｂ変電所から着信した電流値をサンプリングし
たのと同一時刻ｔにおけるＡ変電所での電流値をｘ
（ｔ）とし、時刻ｔの起点はｎ′＝Ｎのサンプリング時
点がｔ＝Ｎ×ｈとなる様にとったものとする。
たのと同一時刻ｔにおけるＡ変電所での電流値をｘ
（ｔ）とし、時刻ｔの起点はｎ′＝Ｎのサンプリング時
点がｔ＝Ｎ×ｈとなる様にとったものとする。
時刻ｔはＡ変電所でのサンプリング時点を基準とした時
にＢ変電所でのサンプリング時点を相対的に示すもの
で、処理ステップS31によって算出した値を基をもとに
決定できる。
にＢ変電所でのサンプリング時点を相対的に示すもの
で、処理ステップS31によって算出した値を基をもとに
決定できる。
次に処理ステップS34によってＢ変電所でサンプリング
された電流値と、処理ステップS63によって算出された
同一時刻のＡ変電所の電流値との和が所定の範囲内で０
に近いかどうかを判定し、０に近くない場合には処理ス
テップS35によってトリップ指令を出す。このトリップ
指令はI/O部（20）を経由してしゃ断器（21）に伝えら
れ、しゃ断器（21）を開く。
された電流値と、処理ステップS63によって算出された
同一時刻のＡ変電所の電流値との和が所定の範囲内で０
に近いかどうかを判定し、０に近くない場合には処理ス
テップS35によってトリップ指令を出す。このトリップ
指令はI/O部（20）を経由してしゃ断器（21）に伝えら
れ、しゃ断器（21）を開く。
なお、上式は、被サンプリング量ｘ（ｔ）が時間周期Ｎ
・ｈで周期的な時に任意のｔに対して成立する。従って
例えば50Hzの交流電流をｈ＝1/600秒毎にサンプリング
する場合に適用するときには、電流は1/50秒で周期的で
あるから、Ｎ＝（1/50）／（1/600）＝12とすればよ
い。ただし、被サンプリング量ｘ（ｔ）には1/2h Hz以
上の周波数成分がないと実用上みなせる場合に上式が成
立する。
・ｈで周期的な時に任意のｔに対して成立する。従って
例えば50Hzの交流電流をｈ＝1/600秒毎にサンプリング
する場合に適用するときには、電流は1/50秒で周期的で
あるから、Ｎ＝（1/50）／（1/600）＝12とすればよ
い。ただし、被サンプリング量ｘ（ｔ）には1/2h Hz以
上の周波数成分がないと実用上みなせる場合に上式が成
立する。
従来の保護継電装置は以上のように、多くのサンプリン
グ値のサンメーションの形で表わされる式（１）に示さ
れるフーリエ係数を算出し、これに基いてさらに電気量
ｘ（ｔ）を求めるようになされていたので、サンプリン
グデータが１点得られる毎に多くの乗算と加算を必要と
していた。
グ値のサンメーションの形で表わされる式（１）に示さ
れるフーリエ係数を算出し、これに基いてさらに電気量
ｘ（ｔ）を求めるようになされていたので、サンプリン
グデータが１点得られる毎に多くの乗算と加算を必要と
していた。
この発明は上記のような問題点を解決するためになされ
たもので、フーリエ係数の算出に要するデジタル演算構
成の負荷を低減し、その分、安価で高機能な保護継電装
置を得ることを目的とする。
たもので、フーリエ係数の算出に要するデジタル演算構
成の負荷を低減し、その分、安価で高機能な保護継電装
置を得ることを目的とする。
この発明に係る保護継電装置は、複数の電気量検出端子
における周期的な各電気量を一定時間周期ｈでサンプリ
ングし、そのサンプリング値をディジタル符号に変換し
た値Ｘ（ｎ×ｈ）をもとに算出したフーリエ係数と端子
間のサンプリング同期ずれ時間情報を用いて同一時刻に
おける各端子の電気量を算出し、電力設備の保護のため
の演算を行う保護継電装置において、上記フーリエ係数
Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）をサンプリングデータＸ
（ｎ×ｈ）が１点得られる毎に下記の漸化式に基いて算
出すると共に、該漸化式によるフーリエ係数Ａ（ｑ、
ｎ）とＢ（ｑ、ｎ）を用いて任意の時刻での複数の電気
量検出端子における周期的な電機量Ｘ（ｔ）を下記の算
出式に基いて算出することを特徴とするものである。
における周期的な各電気量を一定時間周期ｈでサンプリ
ングし、そのサンプリング値をディジタル符号に変換し
た値Ｘ（ｎ×ｈ）をもとに算出したフーリエ係数と端子
間のサンプリング同期ずれ時間情報を用いて同一時刻に
おける各端子の電気量を算出し、電力設備の保護のため
の演算を行う保護継電装置において、上記フーリエ係数
Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）をサンプリングデータＸ
（ｎ×ｈ）が１点得られる毎に下記の漸化式に基いて算
出すると共に、該漸化式によるフーリエ係数Ａ（ｑ、
ｎ）とＢ（ｑ、ｎ）を用いて任意の時刻での複数の電気
量検出端子における周期的な電機量Ｘ（ｔ）を下記の算
出式に基いて算出することを特徴とするものである。
ここで、漸化式は次式で示される。
Ａ（ｑ、ｎ）＝Ａ（ｑ、ｎ−１）−Ｘ（（ｎ−Ｎ）×
ｈ） ×2/N×cos（ｑ×ω×ｍ（ｎ−Ｎ）/N）） ＋Ｘ（ｎ×ｈ）×2/N×cos（ｑ×ω×ｍ（ｎ）/N） Ｂ（ｑ、ｎ）＝Ｂ（ｑ、ｎ−１）−Ｘ（（ｎ−Ｎ）×
ｈ）× 2/N×sin（ｑ×ω×ｍ（ｎ−Ｎ）/N）） ＋Ｘ（ｎ×ｈ）×2/N×sin（ｑ×ω×ｍ（ｎ）/N） 但し、（Ｉ） ｔは時刻、 （II） ｎはサンプリングの時系列を表す整数 （ｎ＝１、２、３、・・・） （III） ｈはサンプリング周期 （IV） Ｎは１周忌Ｔ内のサンプリング点数で、Ｎ＝T/
h （ここで、ｈはＴの整数分の１とする） （Ｖ） ω＝２π/T （VI） Ｖ（（ｎ−Ｎ）×ｈ）はサンプリングした時刻
（ｎ×ｈ）における電気量 Ｘ（ｎ×ｈ）に対してサンプリング数Ｎ前の計測電気量 （VII） Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）はサンプリング
データがｎ点まで得られた時点でのｑ倍調波のフーリエ
係数 なお、ｎ＝０の時、Ａ（ｑ、０）＝Ｂ（ｑ、０）＝０ （VIII） ｑは波サンプリング量の基本周波数の倍調波
数を表す整数 （IX） ｍ（ｎ）＝｛（ｎ−１）をＮで除算した時の余
り＋１｝ （Ｘ） Ｘ（−Ｎ＋１）×ｈ）＝Ｘ（（−Ｎ＋２）×
ｈ）＝・・・＝Ｘ（−ｈ）＝Ｘ（０）＝０ また、算出式は次式で示される。
ｈ） ×2/N×cos（ｑ×ω×ｍ（ｎ−Ｎ）/N）） ＋Ｘ（ｎ×ｈ）×2/N×cos（ｑ×ω×ｍ（ｎ）/N） Ｂ（ｑ、ｎ）＝Ｂ（ｑ、ｎ−１）−Ｘ（（ｎ−Ｎ）×
ｈ）× 2/N×sin（ｑ×ω×ｍ（ｎ−Ｎ）/N）） ＋Ｘ（ｎ×ｈ）×2/N×sin（ｑ×ω×ｍ（ｎ）/N） 但し、（Ｉ） ｔは時刻、 （II） ｎはサンプリングの時系列を表す整数 （ｎ＝１、２、３、・・・） （III） ｈはサンプリング周期 （IV） Ｎは１周忌Ｔ内のサンプリング点数で、Ｎ＝T/
h （ここで、ｈはＴの整数分の１とする） （Ｖ） ω＝２π/T （VI） Ｖ（（ｎ−Ｎ）×ｈ）はサンプリングした時刻
（ｎ×ｈ）における電気量 Ｘ（ｎ×ｈ）に対してサンプリング数Ｎ前の計測電気量 （VII） Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）はサンプリング
データがｎ点まで得られた時点でのｑ倍調波のフーリエ
係数 なお、ｎ＝０の時、Ａ（ｑ、０）＝Ｂ（ｑ、０）＝０ （VIII） ｑは波サンプリング量の基本周波数の倍調波
数を表す整数 （IX） ｍ（ｎ）＝｛（ｎ−１）をＮで除算した時の余
り＋１｝ （Ｘ） Ｘ（−Ｎ＋１）×ｈ）＝Ｘ（（−Ｎ＋２）×
ｈ）＝・・・＝Ｘ（−ｈ）＝Ｘ（０）＝０ また、算出式は次式で示される。
〔作用〕 この発明による保護継電装置においては、フーリエ係数
を算出する時の時刻の起点のとり方を改め、サンプリン
グデータが１点得られた時に、１つ前のデータが得られ
た時に算出したフーリエ係数に若干の演算を行うだけで
最新のフーリエ係数が得られる漸化式によるフーリエ係
数の算出を採用することにより、デジタル演算構成の負
荷を低減させ、安価で高機能な保護継電装置を得ること
ができる。
を算出する時の時刻の起点のとり方を改め、サンプリン
グデータが１点得られた時に、１つ前のデータが得られ
た時に算出したフーリエ係数に若干の演算を行うだけで
最新のフーリエ係数が得られる漸化式によるフーリエ係
数の算出を採用することにより、デジタル演算構成の負
荷を低減させ、安価で高機能な保護継電装置を得ること
ができる。
以下、この発明の一実施例について説明する。
なお、前述した標本化定理におけ標本ｘ（ｔ）を、以下
電気量Ｘ（ｔ）と表記し取り扱う。
電気量Ｘ（ｔ）と表記し取り扱う。
この発明においては、複数の電気量検出端子において時
間周期Ｔで周期的な電気量Ｘ（ｔ）を一定時間周期ｈで
サンプリングし、そのサンプリング値をディジタル符号
に変換した値Ｘ（ｎ×ｈ）が１点得られる毎に次式
（３）によりフーリエ係数を算出する。
間周期Ｔで周期的な電気量Ｘ（ｔ）を一定時間周期ｈで
サンプリングし、そのサンプリング値をディジタル符号
に変換した値Ｘ（ｎ×ｈ）が１点得られる毎に次式
（３）によりフーリエ係数を算出する。
ただし、ｎはサンプリングの時系列を表わす整数（ｎ＝
１、２、３、…）、ｑは被サンプリング量の基本周波数
の倍調波数を表わす整数 ｍ（ｎ）＝｛（ｎ−１）をＮで除算した時の余り＋１｝ Ｘ（（−Ｎ＋１）×ｈ）＝Ｘ（（−Ｎ＋２）×ｈ） ＝…（＝Ｘ（−ｈ）＝Ｘ（０）＝Ｏ Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）はサンプリングデータがｎ
点まで得られた時点でのｑ倍調波のフーリエ係数 Ａ（ｑ、０）＝Ｂ（ｑ、０）＝０ である。
１、２、３、…）、ｑは被サンプリング量の基本周波数
の倍調波数を表わす整数 ｍ（ｎ）＝｛（ｎ−１）をＮで除算した時の余り＋１｝ Ｘ（（−Ｎ＋１）×ｈ）＝Ｘ（（−Ｎ＋２）×ｈ） ＝…（＝Ｘ（−ｈ）＝Ｘ（０）＝Ｏ Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）はサンプリングデータがｎ
点まで得られた時点でのｑ倍調波のフーリエ係数 Ａ（ｑ、０）＝Ｂ（ｑ、０）＝０ である。
係数Ａ（ｑ、ｎ）とＢ（ｑ、ｎ）は、ｎ番目までのサン
プリングが得られた時点で算出されるものであることを
明確にする為にｎを含めて記した。
プリングが得られた時点で算出されるものであることを
明確にする為にｎを含めて記した。
このＡ（ｑ、ｎ）とＢ（ｑ、ｎ）は「ｎ＝（Ｎの整数
倍）でｔ＝０」と考えて時刻の起点をとった場合の最新
のＮ点のデータから算出した式（１）のフーリエ係数に
相当する。これは、以下の理由による。
倍）でｔ＝０」と考えて時刻の起点をとった場合の最新
のＮ点のデータから算出した式（１）のフーリエ係数に
相当する。これは、以下の理由による。
すなわち、Ｘ（ｔ）、cos（ｑ×ω×ｔ）、sin（ｑ×ω
×ｔ）は時間周期Ｔで周期的であるから、第１図に示す
ように、αを適当な整数としてｎ＝α×Ｎの時点をｔ＝
０と考えて、そこからの１周期分のＮ点をｎ′＝１、
２、３、…、Ｎとして式（１）により算出した係数Ａ
（ｑ）とＢ（ｑ）を用いて式（２）より任意の時刻にお
けるＸ（ｔ）の算出が可能となる。ただし、先に述べた
ように時刻の起点はｎ＝α×Ｎの時にｔ＝０とする。こ
こで、第１図に示すように、最新のデータ点Ｎ点と過去
のデータ点Ｎ点との時間差（サンプリング周期の整数
倍）を示すものとして位置付ける値βを適当な整数とす
ると、最新のＮ点はｎ＝α×Ｎ＋β＋１からｎ＝α×Ｎ
＋β＋Ｎ番目と表わせ、かつ、Ｘ（ｔ）、cos（ｑ×ω
×ｔ）、sin（ｑ×ω×ｔ）は時間周期Ｔで周期的であ
るから であり、最新のＮ点に対してＸ（ｎ×ｈ）×2/N×cos
（ｑ×ω×ｍ（ｎ）/N）又はＸ（ｎ×ｈ）×2/N×sin
（ｑ×ω×ｍ（ｎ）/N）を合計したものはｎ＝α×Ｎの
時点をｔ＝０と考え、最新のＮ点をその周期性に従って
ｔ＝ｈ、２×ｈ、３×ｈ、Ｎ×ｈに得られたものと想定
して式（１）によって算出したＡ（ｑ）、Ｂ（ｑ）に一
致するのである。
×ｔ）は時間周期Ｔで周期的であるから、第１図に示す
ように、αを適当な整数としてｎ＝α×Ｎの時点をｔ＝
０と考えて、そこからの１周期分のＮ点をｎ′＝１、
２、３、…、Ｎとして式（１）により算出した係数Ａ
（ｑ）とＢ（ｑ）を用いて式（２）より任意の時刻にお
けるＸ（ｔ）の算出が可能となる。ただし、先に述べた
ように時刻の起点はｎ＝α×Ｎの時にｔ＝０とする。こ
こで、第１図に示すように、最新のデータ点Ｎ点と過去
のデータ点Ｎ点との時間差（サンプリング周期の整数
倍）を示すものとして位置付ける値βを適当な整数とす
ると、最新のＮ点はｎ＝α×Ｎ＋β＋１からｎ＝α×Ｎ
＋β＋Ｎ番目と表わせ、かつ、Ｘ（ｔ）、cos（ｑ×ω
×ｔ）、sin（ｑ×ω×ｔ）は時間周期Ｔで周期的であ
るから であり、最新のＮ点に対してＸ（ｎ×ｈ）×2/N×cos
（ｑ×ω×ｍ（ｎ）/N）又はＸ（ｎ×ｈ）×2/N×sin
（ｑ×ω×ｍ（ｎ）/N）を合計したものはｎ＝α×Ｎの
時点をｔ＝０と考え、最新のＮ点をその周期性に従って
ｔ＝ｈ、２×ｈ、３×ｈ、Ｎ×ｈに得られたものと想定
して式（１）によって算出したＡ（ｑ）、Ｂ（ｑ）に一
致するのである。
ここで、最新のＮ点のデータに対してＸ（ｎ×ｈ）×2/
N×cos（ｑ×ω×ｍ（ｎ）/N）又はＸ（ｎ×ｈ）×2/N
×sin（ｑ×ω×ｍ（ｎ）/N）を合計するということは
データが１点得られる毎に式（３）の漸化式によって計
算するのと同じことである。
N×cos（ｑ×ω×ｍ（ｎ）/N）又はＸ（ｎ×ｈ）×2/N
×sin（ｑ×ω×ｍ（ｎ）/N）を合計するということは
データが１点得られる毎に式（３）の漸化式によって計
算するのと同じことである。
このような時刻の起点をとれば、式（３）によるフーリ
エ係数Ａ（ｑ、ｎ）とＢ（ｑ、ｎ）を用いて式（４）に
より任意の時刻のｘ（ｔ）を算出できる。
エ係数Ａ（ｑ、ｎ）とＢ（ｑ、ｎ）を用いて式（４）に
より任意の時刻のｘ（ｔ）を算出できる。
ただし、サンプリングデータはＮ点以上得られているも
のとする。
のとする。
以上の様に、複数端子間でサンプリングが同期していな
くても、この発明により、複数端子で同一時刻における
電気量を高速に算出でき、すなわち、サンプリング時刻
だけの電気量から任意の時刻での電気量を算出できるも
のであるから、複数端子で電気量を任意のタイミングで
勝手にサンプリングしていても、各端子でのサンプリン
グデータから各端子同一時刻における電気量を得ること
ができる。例えば各端子から同一時刻における流出電流
の和が０であるかどうかを判定して事故検出を行う保護
アルゴリズムの実現が容易となる。
くても、この発明により、複数端子で同一時刻における
電気量を高速に算出でき、すなわち、サンプリング時刻
だけの電気量から任意の時刻での電気量を算出できるも
のであるから、複数端子で電気量を任意のタイミングで
勝手にサンプリングしていても、各端子でのサンプリン
グデータから各端子同一時刻における電気量を得ること
ができる。例えば各端子から同一時刻における流出電流
の和が０であるかどうかを判定して事故検出を行う保護
アルゴリズムの実現が容易となる。
ここで、１つのサンプリングデータが得られる毎にフー
リエ係数を算出する為の演算量は、従来の式（１）によ
るものよりも、この発明の式（３）によるものの方が明
らかに少ないものとなっている。
リエ係数を算出する為の演算量は、従来の式（１）によ
るものよりも、この発明の式（３）によるものの方が明
らかに少ないものとなっている。
また、ｎ＝１、２、３、…に対してｍ（ｎ）は１、２、
３、…、N,1、２、３、…、N,1、２、…、N,…の様に１
からＮをくり返すので、デジタル計算機においては、ｍ
の初期値を１としておいて、以降、サンプリングデータ
が１点得られる毎に１を加算し、Ｎを越えたら１に戻す
というプログラムでｎからｍを得る機能を容易に実現で
きる。
３、…、N,1、２、３、…、N,1、２、…、N,…の様に１
からＮをくり返すので、デジタル計算機においては、ｍ
の初期値を１としておいて、以降、サンプリングデータ
が１点得られる毎に１を加算し、Ｎを越えたら１に戻す
というプログラムでｎからｍを得る機能を容易に実現で
きる。
更に、ｑとｍ（ｎ）に対応する2/n×cos（ｑ×ω×ｍ
（ｎ）/N）や2/n×sin（ｑ×ω×ｍ（ｎ）/N）は高々qm
×Ｎ個しかないので、事前にこれらを計算してメモリに
格納しておけば、サンプリングデータが得られる毎に計
算する必要はない。
（ｎ）/N）や2/n×sin（ｑ×ω×ｍ（ｎ）/N）は高々qm
×Ｎ個しかないので、事前にこれらを計算してメモリに
格納しておけば、サンプリングデータが得られる毎に計
算する必要はない。
第２図はこの発明によるフーリエ係数演算を行いこれに
基づく保護演算を行う保護継電装置のハードウェアブロ
ック図で、第５図と同一部分は同一符号を用いて示して
いる。図中、（６）は低速マイクロプロセッサである。
基づく保護演算を行う保護継電装置のハードウェアブロ
ック図で、第５図と同一部分は同一符号を用いて示して
いる。図中、（６）は低速マイクロプロセッサである。
そして、第３図はこの発明における低速マイクロプロセ
ッサ（６）にて行なう処理のプログラムフローを示した
もので、Ｂ変電所でのサンプリング時刻と同一時刻にお
けるＡ変電所での電流値をこの発明の式に基いて算出す
る処理ステップS33を有している。また、第４図は処理
ステップS33を更に詳細に記したプログラムフローで、
漸化式によるフーリエ係数の算出を行う処理ステップS3
3Aとフーリエ係数をもとに所望時刻の電流値を算出する
処理ステップS33Bとを有している。
ッサ（６）にて行なう処理のプログラムフローを示した
もので、Ｂ変電所でのサンプリング時刻と同一時刻にお
けるＡ変電所での電流値をこの発明の式に基いて算出す
る処理ステップS33を有している。また、第４図は処理
ステップS33を更に詳細に記したプログラムフローで、
漸化式によるフーリエ係数の算出を行う処理ステップS3
3Aとフーリエ係数をもとに所望時刻の電流値を算出する
処理ステップS33Bとを有している。
次に動作について説明する。以下、Ａ変電所側での動作
を中心に説明するが、Ｂ変電所側でも同様の動作であ
る。送電線（11）の電気量である電流をA/D変換部（1
2）によってサンプルホールド後ディジタル値に変換
し、低速マイクロプロセッサ（６）に入力する。一方、
A/D変換部（12）の出力するディジタル値は光電気変換
部（14A）、光電気変換部（14B）、情報伝送部（15）、
変復調部（16）、搬送端局装置（17）、無線装置（18）
及び全２重回線（19）を経由してＢ変電所に送信される
ようになされ、他方Ｂ変電所でサンプリングされた電流
のディジタル値は同様に全２重回線（19）を経由してＡ
変電所に送信され、これを受けてＡ変電所では、このデ
ィジタル値を無線装置（18）、搬送端局装置（17）、変
復調部（16）、情報伝送部（15）、光電気変換部（14
B）及び光電気変換部（14A）を経由して低速マイクロプ
ロセッサ（６）に入力するようになされている。
を中心に説明するが、Ｂ変電所側でも同様の動作であ
る。送電線（11）の電気量である電流をA/D変換部（1
2）によってサンプルホールド後ディジタル値に変換
し、低速マイクロプロセッサ（６）に入力する。一方、
A/D変換部（12）の出力するディジタル値は光電気変換
部（14A）、光電気変換部（14B）、情報伝送部（15）、
変復調部（16）、搬送端局装置（17）、無線装置（18）
及び全２重回線（19）を経由してＢ変電所に送信される
ようになされ、他方Ｂ変電所でサンプリングされた電流
のディジタル値は同様に全２重回線（19）を経由してＡ
変電所に送信され、これを受けてＡ変電所では、このデ
ィジタル値を無線装置（18）、搬送端局装置（17）、変
復調部（16）、情報伝送部（15）、光電気変換部（14
B）及び光電気変換部（14A）を経由して低速マイクロプ
ロセッサ（６）に入力するようになされている。
しかして、Ａ変電所におけるサンプリングとＢ変電所に
おけるサンプリングは相手変電所と非同期な自走クロッ
クを基準にして一定時間周期ｈで行われる。
おけるサンプリングは相手変電所と非同期な自走クロッ
クを基準にして一定時間周期ｈで行われる。
低速マイクロプロセッサ（６）では第３図のフローに示
す処理を行う。この処理はＡ変電所でのA/D変換完了時
とＢ変電所からの電流値受信時に行われるものである。
まず処理ステップS31によってＡ変電所でのサンプリン
グとＢ変電所でのサンプリングの時間差を得る。これは
例えばサンプリング周期ｈよりも充分に短い時間周期毎
にカウントアップされるカウンタ値を、Ａ変電所でのA/
D変換完了時とＢ変電所からの電流値受信時にとり込
み、更に、Ｂ変電所からＡ変電所へ電流値を伝送する伝
送時間を考慮することによって算出できるものである。
Ｂ変電所からＡ変電所への伝送時間は電送路が不変であ
る限り一定時間であり、例えば理論的に算出しておく
か、又は、一度だけ測定しておき、常時その値を用いる
ことで対応できる。
す処理を行う。この処理はＡ変電所でのA/D変換完了時
とＢ変電所からの電流値受信時に行われるものである。
まず処理ステップS31によってＡ変電所でのサンプリン
グとＢ変電所でのサンプリングの時間差を得る。これは
例えばサンプリング周期ｈよりも充分に短い時間周期毎
にカウントアップされるカウンタ値を、Ａ変電所でのA/
D変換完了時とＢ変電所からの電流値受信時にとり込
み、更に、Ｂ変電所からＡ変電所へ電流値を伝送する伝
送時間を考慮することによって算出できるものである。
Ｂ変電所からＡ変電所への伝送時間は電送路が不変であ
る限り一定時間であり、例えば理論的に算出しておく
か、又は、一度だけ測定しておき、常時その値を用いる
ことで対応できる。
次に、処理ステップS32によってＢ変電所からの電流値
受信時には処理ステップS33に分岐し、そうでない時に
は処理を終了し次の処理開始条件成立を待つ。
受信時には処理ステップS33に分岐し、そうでない時に
は処理を終了し次の処理開始条件成立を待つ。
処理ステップS33ではＢ変電所から着信した電流値をサ
ンプリングしたのと同一時刻におけるＡ変電所での電流
値を算出するが、これは漸化式によってフーリエ係数を
算出する処理ステップS33Aとフーリエ係数をもとに所望
の時刻の電流値を算出する処理ステップS33Bとに分けら
れる。
ンプリングしたのと同一時刻におけるＡ変電所での電流
値を算出するが、これは漸化式によってフーリエ係数を
算出する処理ステップS33Aとフーリエ係数をもとに所望
の時刻の電流値を算出する処理ステップS33Bとに分けら
れる。
処理ステップS33Aでは式（３）によるフーリエ係数Ａ
（ｑ、ｎ）、Ｂ（ｑ、ｎ）の算出を行う。
（ｑ、ｎ）、Ｂ（ｑ、ｎ）の算出を行う。
また、処理ステップS33Bでは式（４）による電気量Ｘ
（ｔ）の演算を行う。
（ｔ）の演算を行う。
ただし、Ｘ（ｔ）はＢ変電所から着信した電流値をサン
プリングしたのと同一時刻ｔにおけるＡ変電所での電流
値であり、時刻の起点は、前述の様に「ｎ＝（Ｎの整数
倍）でｔ＝０」と考えてとったものであり、時刻ｔは処
理ステップS31によって算出した値をもとに決定でき
る。
プリングしたのと同一時刻ｔにおけるＡ変電所での電流
値であり、時刻の起点は、前述の様に「ｎ＝（Ｎの整数
倍）でｔ＝０」と考えてとったものであり、時刻ｔは処
理ステップS31によって算出した値をもとに決定でき
る。
次に処理ステップS34によってＢ変電所でサンプリング
された電流値と、処理ステップS33によって算出された
同一時刻のＡ変電所の電流値との和が所定の範囲内で０
に近いかどうかを判定し、０に近くない場合には処理ス
テップS35によってトリップ指令を出す。このトリップ
指令はI/O部（20）を経由してしゃ断器（21）に伝えら
れ、しゃ断器（21）を開く。
された電流値と、処理ステップS33によって算出された
同一時刻のＡ変電所の電流値との和が所定の範囲内で０
に近いかどうかを判定し、０に近くない場合には処理ス
テップS35によってトリップ指令を出す。このトリップ
指令はI/O部（20）を経由してしゃ断器（21）に伝えら
れ、しゃ断器（21）を開く。
なお、上式は、被サンプリング量Ｘ（ｔ）が時間周期Ｎ
・ｈで周期的な時に任意のｔに対して成立する。従って
例えば50Hzの交流電流をｈ＝1/600秒毎にサンプリング
する場合に適用するときには、電流は1/50秒で周期的で
あるから、Ｎ＝（1/50）／（1/600）＝12とすればよ
い。ただし、被サンプリング量Ｘ（ｔ）には1/2h Hz以
上の周波数成分がないと実用上みなせる場合に上式が成
立する。
・ｈで周期的な時に任意のｔに対して成立する。従って
例えば50Hzの交流電流をｈ＝1/600秒毎にサンプリング
する場合に適用するときには、電流は1/50秒で周期的で
あるから、Ｎ＝（1/50）／（1/600）＝12とすればよ
い。ただし、被サンプリング量Ｘ（ｔ）には1/2h Hz以
上の周波数成分がないと実用上みなせる場合に上式が成
立する。
なお、この発明は、各端子のサンプリング値自体を端子
間で通信する他、各端子において算出したフーリエ係数
を端子間で通信しても同様の効果を発揮する。ただし、
この場合、各端子でどのように時刻の起点をとっている
かを知る情報が必要なので、例えば、フーリエ係数Ａ
（q.n）、Ｂ（ｑ、ｎ）に加えてｍ（ｎ）を通信するこ
とが必要である。
間で通信する他、各端子において算出したフーリエ係数
を端子間で通信しても同様の効果を発揮する。ただし、
この場合、各端子でどのように時刻の起点をとっている
かを知る情報が必要なので、例えば、フーリエ係数Ａ
（q.n）、Ｂ（ｑ、ｎ）に加えてｍ（ｎ）を通信するこ
とが必要である。
フーリエ係数は先に示したｑの全ての値に対して算出し
なくても、保護演算のために必要な周波数成分に対する
ものだけでもよい。例えば50Hzの電力系統にｈ＝1/600
秒Ｎ＝12にて適用する場合、その基本周波数50Hz成分の
みが必要ならｑ＝１に対してのみ計算すればよい。
なくても、保護演算のために必要な周波数成分に対する
ものだけでもよい。例えば50Hzの電力系統にｈ＝1/600
秒Ｎ＝12にて適用する場合、その基本周波数50Hz成分の
みが必要ならｑ＝１に対してのみ計算すればよい。
更に、サンプリングした電気量とは意味の異なる量をフ
ーリエ係数から求め、これによ事故判定を行う保護継電
装置においても同様の効果を発揮する。例えばフーリエ
係数のＡ（ｑ、ｎ）を実部、Ｂ（ｑ、ｎ）を虚部にとっ
た複素ベクトルは、その大きさが被サンプリング量のｑ
倍調波成分の大きさで、偏角が被サンプリング量のｑ倍
調波成分の位相を表わすので、これをもとに事故判別を
行う保護継電装置においても同様の効果を発揮する。
ーリエ係数から求め、これによ事故判定を行う保護継電
装置においても同様の効果を発揮する。例えばフーリエ
係数のＡ（ｑ、ｎ）を実部、Ｂ（ｑ、ｎ）を虚部にとっ
た複素ベクトルは、その大きさが被サンプリング量のｑ
倍調波成分の大きさで、偏角が被サンプリング量のｑ倍
調波成分の位相を表わすので、これをもとに事故判別を
行う保護継電装置においても同様の効果を発揮する。
又、複数の端子で非同期でサンプリングしていて各端子
で上記の複素ベクトルを得ている場合には、サンプリン
グの周期ずれ時間に比例する角度だけ複素平面上でベク
トルを回転させることによって、同一時刻に各端子がサ
ンプリグしたのと同じ複素ベクトルを各端子について得
ることが出来るので、これらのベクトルをもとに事故判
別を行う保護継電装置に対してもこの発明は同様の効果
を発揮する。
で上記の複素ベクトルを得ている場合には、サンプリン
グの周期ずれ時間に比例する角度だけ複素平面上でベク
トルを回転させることによって、同一時刻に各端子がサ
ンプリグしたのと同じ複素ベクトルを各端子について得
ることが出来るので、これらのベクトルをもとに事故判
別を行う保護継電装置に対してもこの発明は同様の効果
を発揮する。
この発明においては、ｍ（ｎ）を１からＮまでとした
が、これを０からＮ−１としても若干の式の変更を行う
だけで同様の効果を実現できる。
が、これを０からＮ−１としても若干の式の変更を行う
だけで同様の効果を実現できる。
以上の様に、この発明によれば、フーリエ係数を算出す
る時の時刻の起点のとり方を改め、サンプリングデータ
が１点得られた時に、１つ前のデータが得られた時に算
出したフーリエ係数に若干の演算を加えるだけで最新の
フーリエ係数が得られる漸化式による演算を実現し、こ
れを採用して保護継電装置を構成したので、デジタル演
算構成の負荷を低減させ、安価で高機能な保護継電装置
を得ることができる。
る時の時刻の起点のとり方を改め、サンプリングデータ
が１点得られた時に、１つ前のデータが得られた時に算
出したフーリエ係数に若干の演算を加えるだけで最新の
フーリエ係数が得られる漸化式による演算を実現し、こ
れを採用して保護継電装置を構成したので、デジタル演
算構成の負荷を低減させ、安価で高機能な保護継電装置
を得ることができる。
第１図はこの発明におけるフーリエ係数の算出における
時刻起点のとり方の説明図、第２図はこの発明の一実施
例を示す保護継電装置のハードウェアブロック図、第３
図はこの発明によるプログラムフロー図、第４図は第３
図中の処理ステップS33を詳細に示すプログラムフロー
図、第５図は従来の保護継電装置のハードウェアブロッ
ク図、第６図は従来の方法によるプログラムフロー図、
第７図は第６図中の処理ステップS63を詳細に示すプロ
グラムフロー図である。 図中、（６）は低速マイクロプロセッサ（電気量算出手
段）、（11）は送電線、（12）はA/D変換部、S31はＡ変
電所でのサンプリングとＢ変電所でのサンプリングの時
間差を得る処理ステップ、S32は処理の起動要因判別処
理ステップ、S33はＢ変電所でのサンプリング時刻と同
一時刻におけるＡ変電所での電流値をこの発明による式
に基いて算出する処理ステップ、S33Aは漸化式によるフ
ーリエ係数の算出処理ステップ、S33Bはこの発明による
時刻起点のとり方による所望時刻の電流値算出処理ステ
ップである。 なお、各図中、同一符号は同一又は相当部分を示す。
時刻起点のとり方の説明図、第２図はこの発明の一実施
例を示す保護継電装置のハードウェアブロック図、第３
図はこの発明によるプログラムフロー図、第４図は第３
図中の処理ステップS33を詳細に示すプログラムフロー
図、第５図は従来の保護継電装置のハードウェアブロッ
ク図、第６図は従来の方法によるプログラムフロー図、
第７図は第６図中の処理ステップS63を詳細に示すプロ
グラムフロー図である。 図中、（６）は低速マイクロプロセッサ（電気量算出手
段）、（11）は送電線、（12）はA/D変換部、S31はＡ変
電所でのサンプリングとＢ変電所でのサンプリングの時
間差を得る処理ステップ、S32は処理の起動要因判別処
理ステップ、S33はＢ変電所でのサンプリング時刻と同
一時刻におけるＡ変電所での電流値をこの発明による式
に基いて算出する処理ステップ、S33Aは漸化式によるフ
ーリエ係数の算出処理ステップ、S33Bはこの発明による
時刻起点のとり方による所望時刻の電流値算出処理ステ
ップである。 なお、各図中、同一符号は同一又は相当部分を示す。
Claims (1)
- 【請求項１】複数の電気量検出端子における周期的な各
電気量を一定時間周期ｈでサンプリングし、そのサンプ
リング値をディジタル符号に変換した値Ｘ（ｎ×ｈ）を
もとに算出したフーリエ係数と端子間のサンプリング同
期ずれ時間情報を用いて同一時刻における各端子の電気
量を算出し、電力設備の保護のための演算を行う保護継
電装置において、上記フーリエ係数Ａ（ｑ、ｎ）、Ｂ
（ｑ、ｎ）をサンプリングデータＸ（ｎ×ｈ）が１点得
られる毎に下記の漸化式（１）に基いて算出すると共
に、該漸化式によるフーリエ係数Ａ（ｑ、ｎ）とＢ
（ｑ、ｎ）を用いて任意の時刻での複数の電気量検出端
子における周期的な電気量Ｘ（ｔ）を下記の算出式
（２）に基いて算出することを特徴とする保護継電装
置。 ここで、漸化式（１）は次式で示される。 Ａ（ｑ、ｎ）＝Ａ（ｑ、ｎ−１）−Ｘ（（ｎ−Ｎ）×
ｈ） ×2/N×cos（ｑ×ω×ｍ（ｎ−Ｎ）/N）） ＋Ｘ（ｎ×ｈ）×2/N×cos（ｑ×ω×ｍ（ｎ）/N） Ｂ（ｑ、ｎ）＝Ｂ（ｑ、ｎ−１）−Ｘ（（ｎ−Ｎ）×
ｈ） ×2/N×sin（ｑ×ω×ｍ（ｎ−Ｎ）/N）） ＋Ｘ（ｎ×ｈ）×2/N×sin（ｑ×ω×ｍ（ｎ）/N） 但し、（Ｉ） ｔは時刻、 （II） ｎはサンプリングの時系列を表す整数 （ｎ＝１、２、３、・・・） （III） ｈはサンプリング周期 （IV） Ｎは１周期Ｔ内のサンプリング点数で、Ｎ＝T/
h （ここで、ｈはＴの整数分の１とする） （Ｖ） ω＝２π/T （VI） Ｖ（（ｎ−Ｎ）×ｈ）はサンプリングした時刻
（ｎ×ｈ）における電気量 Ｘ（ｎ×ｈ）に対してサンプリング数Ｎ前の計測電気量 （VII） Ａ（ｑ、ｎ）、Ｂ（ｑ、ｎ）はサンプリング
データがｎ点まで得られた時点でのｑ倍調波のフーリエ
係数 なお、ｎ＝０の時、Ａ（ｑ、０）＝Ｂ（ｑ、０）＝０ （VIII） ｑは被サンプリング量の基本周波数の倍調波
数を表す整数 （IX） ｍ（ｎ）＝｛（ｎ−１）をＮで除算した時の余
り＋１｝ （Ｘ） Ｘ（−Ｎ＋１）×ｈ）＝Ｘ（（−Ｎ＋２）×
ｈ）＝・・・＝Ｘ（−ｈ）＝Ｘ（０）＝０ また、算出式（２）は次式で示される。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63064217A JPH07114534B2 (ja) | 1988-03-16 | 1988-03-16 | 保護継電装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63064217A JPH07114534B2 (ja) | 1988-03-16 | 1988-03-16 | 保護継電装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH01238416A JPH01238416A (ja) | 1989-09-22 |
| JPH07114534B2 true JPH07114534B2 (ja) | 1995-12-06 |
Family
ID=13251705
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63064217A Expired - Fee Related JPH07114534B2 (ja) | 1988-03-16 | 1988-03-16 | 保護継電装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07114534B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CH683721A5 (de) * | 1990-05-03 | 1994-04-29 | Landis & Gyr Business Support | Verfahren zur Ermittlung von Schätzwerten der Momentanwerte von Parametern mindestens eines sinusförmigen Signals mit konstanter und vorbekannter Frequenz. |
| CN116048376A (zh) * | 2022-12-30 | 2023-05-02 | 四方继保（武汉）软件有限公司 | 一种消除冒大数导致误动的方法及系统 |
-
1988
- 1988-03-16 JP JP63064217A patent/JPH07114534B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH01238416A (ja) | 1989-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | A precise calculation of power system frequency and phasor | |
| US6519537B1 (en) | Apparatus providing on-line indication of frequency of an AC electric power system | |
| US5671112A (en) | Digital integrator V/Hz relay for generator and transformer over-excitation protection | |
| CA1239664A (en) | Measuring circuit and method for power distribution equipment | |
| CN107782966A (zh) | 确定交流信号的频率 | |
| US4073009A (en) | Apparatus for calculating amplitude values of sinusoidal waves | |
| JPH06105266B2 (ja) | デジタル保護継電器 | |
| JPH04230869A (ja) | サイン形状信号のパラメータの瞬時値の推定値を求める方法及び装置 | |
| EP0371192B1 (en) | Electric quantity detecting method | |
| JPH07114534B2 (ja) | 保護継電装置 | |
| US4212046A (en) | Distance relaying systems | |
| Fakruddin et al. | Application of Haar functions for transmission line and transformer differential protection | |
| Kamwa et al. | A Kalman filter‐based technique for combined digital estimation of voltage flicker and phasor in power distribution systems | |
| Stojanovic et al. | Table based algorithm for inverse-time overcurrent relay | |
| US4073008A (en) | Apparatus for calculating amplitude values | |
| JP3047036B2 (ja) | 電力測定装置 | |
| EP0367563B1 (en) | Detector of quantity of electricity | |
| EP0801745B1 (en) | Measuring method for determining the amplitude and phase of the fundamental tone of an alternating voltage | |
| US11906557B2 (en) | Processing sinewave signals of variable frequency in a device with fixed processing rates | |
| RU2848935C1 (ru) | Измеритель симметричных составляющих многофазной сети переменного тока | |
| JPH08126190A (ja) | 電力系統のディジタル保護方法 | |
| Phadke | Precise synchronization of phasor measurements in electric power systems | |
| JPH0726976B2 (ja) | 電気量検出器 | |
| JPH0435968B2 (ja) | ||
| JP2508175B2 (ja) | ディジタルリレ―の動作判定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |