JP2002040067A - 次数間高調波検出方法 - Google Patents

次数間高調波検出方法

Info

Publication number
JP2002040067A
JP2002040067A JP2000226365A JP2000226365A JP2002040067A JP 2002040067 A JP2002040067 A JP 2002040067A JP 2000226365 A JP2000226365 A JP 2000226365A JP 2000226365 A JP2000226365 A JP 2000226365A JP 2002040067 A JP2002040067 A JP 2002040067A
Authority
JP
Japan
Prior art keywords
order
inter
injection
harmonic
frequency
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.)
Granted
Application number
JP2000226365A
Other languages
English (en)
Other versions
JP3788212B2 (ja
Inventor
Tokuo Emura
徳男 江村
Giko Haneda
儀宏 羽田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nissin Electric Co Ltd
Original Assignee
Nissin Electric Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nissin Electric Co Ltd filed Critical Nissin Electric Co Ltd
Priority to JP2000226365A priority Critical patent/JP3788212B2/ja
Priority to US09/803,945 priority patent/US6545885B2/en
Publication of JP2002040067A publication Critical patent/JP2002040067A/ja
Application granted granted Critical
Publication of JP3788212B2 publication Critical patent/JP3788212B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Abstract

(57)【要約】 【課題】 電力系統の電圧又は/及び電流の計測信号を
固定周波数サンプリング方式でサンプリングし、フーリ
エ変換のフィルタ演算により電力系統に注入した次数間
高調波を抽出する際、系統基本波及び高調波の影響を演
算により極力少なくする。 【解決手段】 固定周波数サンプリング方式でサンプリ
ングした電力系統5のサンプリングデータ結果のフーリ
エ変換のフィルタ演算により注入次数の次数間高調波及
び注入次数の上,下の非注入次数の次数間高調波の成分
を抽出し、両非注入次数の次数間高調波の抽出成分の平
均から注入次数の次数間高調波の抽出成分に含まれた系
統基本波及び整数倍周波数の高調波の系統成分に基づく
誤差を補間演算し、注入次数の次数間高調波の成分から
誤差を減算して除去し、注入次数の次数間高調波を検出
する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、いわゆる分散型電
源の単独運転検出等に用いられる次数間高調波検出方法
に関する。
【0002】
【従来の技術】従来、電力系統の事故停電等による停止
時に需要家の分散型電源の単独運転を検出してその運転
を停止するため、本出願人は、特開平10−24816
8号公報(H02J 3/38)等に記載されているよ
うに、電力系統に系統基本波に同期したその非整数倍の
周波数の次数間高調波(中間次数調波)を注入し、系統
の電圧又は/及び電流の計測信号から次数間高調波を検
出し、この検出結果に基づき系統の停止を検出する発明
を、既に出願している。
【0003】この場合、次数間高調波が系統に本来存在
しない周波数であることから、例えば系統電源の0.1
%程度の少量の次数間高調波の供給(注入)で系統の停
止を検出できる利点がある。
【0004】
【発明が解決しようとする課題】前記公報に記載されて
いるように系統電圧にPLL同期したサンプリング周波
数のタイミング指令を形成し、この指令に基づいて系統
の電圧,電流の計測信号をサンプリングすると、そのサ
ンプリングデータを公知のフーリエ変換によるデジタル
フィルタ処理により、系統基本波及びその整数倍周波数
の高調波,すなわち系統成分の影響を受けることなく、
注入次数の次数間高調波を精度よく抽出して検出し得る
が、複雑なPLL回路等を要し、安価かつ簡単に検出す
ることができない欠点がある。
【0005】一方、PLL同期をとるためのPLL回路
を省き、電力系統の電圧(以下系統電圧という)とは非
同期にサンプリングのタイミング指令を形成し、この指
令に基づいて計測信号をサンプリングすると、PLL回
路を用いる場合より安価かつ簡単な構成になるが、フー
リエ変換の抽出に系統電圧の誤差が生じ、精度の高い検
出が困難で、場合によっては検出できなくなる問題点が
ある。
【0006】そして、前記の分散型電源の単独運転の検
出の場合だけでなく、次数間高調波を系統に注入し、系
統の電圧又は/及び電流の計測信号をサンプリングして
計測信号に含まれた次数間高調波をフーリエ変換の手法
で抽出し、系統の高調波特性を計測する場合等にも、固
定周波数サンプリング方式を採用すると同様の問題点が
生じる。
【0007】本発明は、電力系統の電圧又は/及び電流
の計測信号を周波数一定の固定周波数サンプリング方式
でサンプリングし、そのサンプリングデータのフーリエ
変換によるフィルタ演算により、系統成分の影響を極力
少なくして注入次数の次数間高調波を精度よく抽出して
検出することを課題とする。
【0008】
【課題を解決するための手段】前記の課題を解決するた
めに、本発明の次数間高調波検出方法は、電力系統の電
圧又は/及び電流をサンプリング周波数一定の固定周波
数サンプリング方式でサンプリングし、サンプリングデ
ータのフーリエ変換によるフィルタ演算により注入次数
の次数間高調波及び注入次数の上,下の非注入次数の次
数間高調波の成分を抽出し、両非注入次数の次数間高調
波の抽出成分の平均から注入次数の次数間高調波の抽出
成分に含まれた系統電圧に基づく誤差を補間演算し、注
入次数の次数間高調波の抽出成分から誤差を減算して除
去し、注入次数の次数間高調波を検出する。
【0009】したがって、電力系統に注入された次数間
高調波を、PLL回路等を用いることなく、サンプリン
グ周波数を固定して安価かつ容易にサンプリングして抽
出することができる。
【0010】そして、サンプリングデータのフーリエ変
換で抽出された注入次数の次数間高調波につき、その系
統周波数とサンプリング周波数との非同期に起因した系
統電圧に基づく誤差が注入次数の上,下の非注入次数の
次数間高調波の成分にも含まれ、その次数(周波数)の
近い範囲では、誤差が線形変化することから、それらの
成分の平均値から注入次数での前記誤差が求められて注
入次数の抽出成分から減算除去される。
【0011】そのため、固定周波数サンプリング方式で
精度よく注入次数の次数間高調波を検出することができ
る。
【0012】そして、注入次数の次数間高調波が、k次
高調波とk+1次高調波(kは2以上の整数)との間の
次数間高調波であり、次数間高調波の注入周波数を、k
f±mf0 ,(fは系統基本波周波数,mは1,2,
…,n−1の整数(nは2以上の整数,f0は次数間高
調波の注入間隔を示し、f0=f/n)から決定するこ
とが実用的で好ましい。
【0013】
【発明の実施の形態】本発明の実施の1形態につき、図
1〜図3を参照して説明する。図1は電力系統の1例の
単線結線図であり、上位系統1に変電所2の1又は複数
の変圧器3の1次側が接続され、各変圧器3の2次側か
ら遮断器4を介して1又は複数の下位の電力系統5が分
枝状に引出される。
【0014】これらの電力系統5は、分散型電源6を有
する需要家設備7,分散型電源6が設けられていない一
般需要家設備8等の複数の需要家設備が接続される。
【0015】そして、本発明が適用される需要家設備7
は他の需要家設備と同様、電力系統5に引込線9の遮断
器10を介して負荷母線11が接続され、この負荷母線
11に各負荷フィーダ12の変圧器13を介してそれぞ
れの負荷が接続される。
【0016】また、負荷母線11に遮断器14が接続さ
れ、この遮断器14に解列用の開閉器15を介して分散
型電源6が接続されるとともに、次数間高調波の電流注
入装置16が接続される。
【0017】この電流注入装置16は次数間高調波の注
入電流を出力するインバータ等の電源部17,この電源
部17と負荷母線11との間に設けられた注入用の変圧
器18により形成される。
【0018】また、引込線9の遮断器10より負荷母線
11側に、受電点変圧器19及び受電点変流器20が設
けられ、それらの電圧,電流の計測信号が系統停止検出
処理装置21のサンプル・ホールド回路22に供給され
る。
【0019】このサンプル・ホールド回路22は水晶発
振器等の定周波発生器で形成されたタイミング指令部2
3の一定周波数のサンプリング指令のタイミング信号に
より、固定周波数サンプリング方式で系統電圧,系統電
流に非同期に両計測信号をサンプル・ホールドし、その
出力が後段のA/D変換回路24によりデジタル信号に
変換され、サンプリングデータとなる。
【0020】そして、A/D変換回路24によりデジタ
ル信号に変換された電圧,電流のサンプリングデータが
演算処理部25に供給され、この処理部25はマイクロ
コンピュータ等で形成され、そのソフトウェア処理によ
り、両サンプリングデータを公知のフーリエ変換のデジ
タルフィルタ演算により電流注入装置16から電力系統
5に注入された注入次数の次数間高調波を抽出して検出
し、その変化から遮断器4が解放する系統停止,すなわ
ち分散型電源6の単独運転を監視して検出する。
【0021】さらに、この単独運転の検出時は、演算処
理部25から開閉器15に解列の指令が供給されて開閉
器15が開放され、分散型電源6が電力系統5から切離
される。
【0022】ところで、系統正常時は、遮断器4,1
0,14及び開閉器15がいずれも閉成され、上位系統
1の電力が電力系統5に給電され、電力系統5が電力供
給状態にある。
【0023】このとき、分散型電源6は電力系統5に連
系して運転され、その出力は自設備7内で消費されると
ともに余剰分が引込線9を介して電力系統5に出力され
る。
【0024】さらに、需要家設備7においては、受電点
変圧器19,受電点変流器20により図1の受電点Aの
電圧及び受電点Aを出入する電流が常時計測される。
【0025】また、演算処理部25がタイミング指令部
23の一定周波数のタイミング信号に同期して電源部1
7に周期的に起動指令を出力し、この指令に基づき、電
源部17が前記タイミング信号に同期した1又は複数の
周波数の次数間高調波の電流を形成し、これらの系統基
本波に非同期の注入電流が変圧器18,負荷母線11,
引込線9を介して需要家設備7の受電点Aから電力系統
5に注入される。
【0026】この注入に基づき、受電点変圧器19,受
電点変流器20の計測信号に次数間高調波の電圧,電流
が含まれる。
【0027】そして、これらの計測信号はタイミング指
令部23のタイミング信号に基づき、前記したようにサ
ンプル・ホールド回路22の固定周波数サンプリング方
式のサンプリングにより、電力系統5の系統電圧,系統
電流に非同期にサンプル・ホールドされる。
【0028】さらに、電圧,電流のホールド出力がA/
D変換回路24によりデジタル信号に変換され、これら
のデジタル信号に変換されたサンプリングデータが演算
処理部25で処理され、注入された次数間高調波の電
圧,電流が抽出されて検出される。
【0029】そして、電力系統5の停止により遮断器4
が開放されると、受電点Aからみた電力系統5の電源イ
ンピーダンスが短絡インピーダンスから開放インピーダ
ンスに変わるため、受電点Aからみた電力系統5の注入
次数の次数間高調波についてのインピーダンス又はアド
ミタンス(以下インピーダンス等という)の変化から電
力系統5の停止,すなわち分散型電源6の単独運転を検
出できる。
【0030】また、演算の簡素化等を図る場合、検出し
た次数間高調波の電圧又は/及び電流の変化から分散型
電源6の単独運転を検出してもよい。
【0031】そこで、演算処理部25は、注入次数の次
数間高調波の電圧又は/及び電流の変化或いは電力系統
5の注入次数の次数間高調波のインピーダンス又はアド
ミタンスの変化から、系統停止の発生に基づく分散型電
源6の単独運転を監視して検出する。
【0032】このとき、注入装置16の次数間高調波の
電流注入及びサンプル・ホールド回路22のサンプリン
グが、タイミング指令部23のタイミング信号により、
系統電圧,系統電流に非同期に行われるため、注入次数
の次数間高調波のサンプリングデータにこの非同期に起
因する誤差が含まれ、フーリエ変換で抽出される注入次
数の次数間高調波の抽出成分に前記誤差の影響が含まれ
る。
【0033】そこで、演算処理部25は、つぎのように
して前記誤差を補正演算して低減,除去する。
【0034】まず、電力系統5に注入されてサンプリン
グされる次数間高調波について説明する。
【0035】電力系統5の60Hz又は50Hzの系統基本
波の周波数(商用電源の規定の周波数)をf,次数間高
調波の注入周波数をfinj ,その注入数を定める周波数
fの分割数をn,この分割数nに基づく注入周波数fin
j の間隔をf0 とすると、f 0=f/nであり、角周波
数ω0はω0=2πf0=2πf/nである。そして、f
=60Hz,n=30であれば、f0 は60/30=2Hz
である。
【0036】また、次数間高調波はk次高調波とk+1
次高調波の間にあり(k=2,3,…)、それらをk次
高調波に近いものから順のmチャンネル(m=1,2,
3,…,(n−1))とすると、k.m次の次数間高調
波の周波数finj ,その角周波数ωinjはつぎの数1,
数2の2式で表される。
【0037】
【数1】
【0038】
【数2】
【0039】このとき、k.m次の次数間高調波をFin
j(t)とすると(tは時間)、それの振幅を1として、
又、それの初期位相をθmとして、この次数間高調波F
inj(t)はつぎの数3の式で表すことができる。
【0040】
【数3】
【0041】一方、電力系統5の系統基本波及びその整
数倍周波数の高調波,すなわち系統成分については、そ
の周波数変動Δを考慮すると、その基本波からl次高調
波までの既存の系統成分F(t) はl次高調波の振幅をV
l,初期位相をφlとして、つぎの数4の式で表わされ
る。
【0042】
【数4】
【0043】そのため、電力系統5に次数間高調波の電
流又は電圧が注入されたときに、変圧器19,変流器2
0の計測信号,すなわちサンプル・ホールド回路22の
電圧,電流に含まれる入力信号Fin(t) は、数3,数4
の式を加算合成したつぎの数5の式に示すようになる。
【0044】
【数5】
【0045】そして、この数5の式の右辺の第1項(F
(t))と第2項(Finj(t))とは、F(t) が周波数変動
Δを含むのに対して、Finj(t)がそれを含まないことか
ら、非同期である。
【0046】つぎに、フーリエ変換(フーリエ積分)に
よる次数間高調波の抽出演算式と系統成分の誤差につい
て、フーリエ変換の積分原理式から説明する。
【0047】まず、入力信号Fin(t) に含まれる数3の
式で示されるk.m次の次数間高調波Finj(t)の余弦成
分(cos成分) Cm,正弦成分(sin成分) Smはつぎ
の数6,数7の公知のフーリエ変換(フーリエ積分)で
抽出することができる。なお、これらの数6,数7の式
が系統周波数とは非同期の周波数固定サンプリング方式
でのデジタルフィルタ演算の原理式となるものである。
【0048】
【数6】
【0049】
【数7】
【0050】そして、数6,数7の式のFin(t) に数5
の式を代入すると、つぎの数8,数9の式が得られる。
【0051】
【数8】
【0052】
【数9】
【0053】この数8,数9の式の右辺の第1項(積分
項)がそれぞれ系統成分の影響に基づく誤差ΔC,ΔS
の項である。
【0054】そして、誤差ΔC,ΔSに基づく抽出成分
全体の誤差ΔS+jΔCは、つぎの数10の式から求ま
る。
【0055】
【数10】
【0056】この数10の式の最右辺の第1項,第2項
の時間積分の結果は、つぎの数11,数12に示すよう
になる。
【0057】
【数11】
【0058】
【数12】
【0059】そして、数10の式に数11,数12の式
を代入して整理すると、誤差ΔS,ΔCは、つぎの数1
3,数14の式に示すようになる。
【0060】
【数13】
【0061】
【数14】
【0062】したがって、数5の式の入力信号Fin(t)
からフーリエ変換のフィルタ演算で注入次数(k.m
次)の次数間高調波を抽出すると、その余弦成分Cm,
正弦成分Smは、つぎの数15,数16の式に示すよう
に、本来の余弦成分(cos成分),正弦成分(sin成分)
と、系統成分(系統基本波及び高調波)に基づく誤差成
分ΔC,ΔSとの和になる。
【0063】
【数15】
【0064】
【数16】
【0065】つぎに数15,数16の式の誤差成分Δ
C,ΔSの振幅の具体的な値について説明する。
【0066】まず、誤差ΔCの振幅GCmはつぎの数1
7の式から求まり、誤差ΔSの振幅GSmはつぎの数1
8の式から求まる。
【0067】
【数17】
【0068】
【数18】
【0069】そして、n=30,k=1,2,3,4,
周波数変動Δ=−5%〜+5%,l=1,2,3,4,
5,6,7,m=1,2,…,29としてシミュレーシ
ョンにより検証したところ、つぎの表1に示すように、
例えばn=30,k=3,Δ=−5%,l=1,m=1
のときに、GCm=−0.00776,GSm=−0.
00243になることがわかった。
【0070】
【表1】
【0071】なお、表1のCm,SmはゲインGCm,
GSmを示す。
【0072】そして、前記のシミュレーションによる
と、例えばk=3,m=15の3.15次の次数間高調
波の周波数は60Hz系統で60Hz×3+2Hz×15=1
80+30=210Hzであり、このとき、系統の周波数
変動Δが−2.5%であれば、基本波(l=1)に基づ
く誤差(残留誤差)の振幅はGCm=0.00465
(=0.465%),GSm=0.00129(=0.
129%)になり、周波数変動Δが−5%であれば、基
本波に基づく誤差ΔC,ΔSの振幅はGCm=−0.0
0655(=−0.655%),GSm=−0.001
78(=−0.178%)になる。一方、基本波を10
0%とすると、次数間高調波の注入量は0.1%程度で
ある。
【0073】したがって、前記の各誤差ΔC,ΔSの振
幅(ゲイン)GCm,GSmは、次数間高調波の注入量
と同じオーダ(パーセント)の大きさで抽出したcos 成
分、sin 成分に含まれ、この場合、誤差ΔC,ΔSが大
きく、数15,数16の抽出した余弦成分Cm,正弦成
分Smそのままでは、これらから本来の次数間高調波の
成分cosθm,sinθmを分離して精度よく検出すること
はできない。
【0074】そして、つぎの数19,数20に示す数1
5,数16の系統成分による誤差ΔC,ΔSは振幅GC
m,GSmがk,mに依存して変化し、しかも、mの近
接値間では、mに対してほぼ線形に変化する。
【0075】
【数19】
【0076】
【数20】
【0077】具体的には、k=3,Δ=−5%の条件
下、l=1,3におけるm=5,6,7の振幅GC5
GC6,GC7及びGS5,GS6,GS7は、表1より、
つぎのようになる。 l=1: GC5=-0.00736,GC6=-0.00727,GC7
=-0.00718 GS5=-0.00221,GS6=-0.00216,GS7=-0.00211 l=3: GC5=0.03527 , GC6=0.03207 ,GC7
=0.02942 GS5=0.03174 ,GS6=0.02856 ,GS7=0.02593
【0078】そして、振幅GC6,GS6については、つ
ぎの数21〜数24に示すように、その上,下の中間次
数調波の誤差成分ΔC,ΔSのゲインGC5とGC7,G
5とGS7の平均値から補間演算して求めることができ
る。
【0079】
【数21】
【0080】
【数22】
【0081】
【数23】
【0082】
【数24】
【0083】すなわち、m−1,m,m+1のk.m−
1次,k.m次,k.m+1次の隣接3チャンネルの次
数間高調波間において、系統成分の誤差(ΔC,ΔS)
の振幅(GCm-1,GSmm-1),(GCm,GS
m),(GCm+1,GSmm+1)は、図2の破線に示す
ようにその周波数にほぼ線形変化し、中間のk.m次の
系統成分の誤差ΔC,ΔSのゲインGCm,GSmは、
その上,下のk.m−1次,k.m+1次の系統成分の
誤差ΔC,ΔSのゲインGCm-1とGCm+1,GSm-1
とGSmm+1の平均値を求めるつぎの数25,数26の
補間演算からゲインGCm’,GSm’として求めるこ
とができる。なお、図2は振幅GCm-1 ,GCm,G
Sm+1を示し、振幅GSm-1,GSm,GSm+1につい
ても同様である。
【0084】
【数25】
【0085】
【数26】
【0086】そして、抽出した成分Cm,Smに含まれ
た誤差ΔC,ΔSの振幅GCm,GSmと、数25,数
26の補間演算から求まる振幅GCm’,GSm’につ
き、つぎの数27,数28の式からそれらの差の振幅G
Cm”,GSm”をシミュレーションして演算したとこ
ろ、つぎの表2の結果が得られた。
【0087】
【数27】
【0088】
【数28】
【0089】
【表2】
【0090】なお、表2のCm’,Sm’が振幅GC
m”,GSm”であり、この表2から明らかなように、
振幅GCm”,GSm”はほぼ0になり、前記の補間演
算から極めて精度よく振幅GCm,GSmが求まること
が確かめられた。
【0091】そこで、k.m次の次数間高調波の電流の
みが注入され、その上,下のk.m−1次、k.m+1
次の次数間高調波が非注入の場合、演算処理部25はフ
ーリエ変換のフィルタ演算によりk.m−1次, k.m
+1次の成分Cm-1 ,Cm+ 1,Sm-1,Sm+1を抽出
し、成分Cm-1,Cm+1の平均値,成分Sm-1,Sm
+1の平均値をk.m次の誤差成分ΔCm,ΔSmとし
て補間演算し、これらの平均値を同時に抽出したk.m
次の成分Cm,Smから減算し、成分Cm,Smに含ま
れた数15,数16の誤差成分ΔC,ΔSを除去して成
分Cm,Smを補正し、本来の次数間高調波の余弦成分
(cosθm),正弦成分(sinθm)を精度よく抽出して
検出する。
【0092】この場合、複雑なPLL回路を用いること
なく安価かつ簡単な固定周波数サンプリング方式で計測
信号のサンプリングを行うことができる。
【0093】そして、サンプリングデータのフーリエ変
換で得られた注入次数(k.m次)の次数間高調波の抽
出成分に、その上,下のk.m−1次、k.m+1次の
非注入次数の抽出成分の平均値に基づく簡単な補正演算
を施して、系統成分の影響をほとんど受けることなく、
注入次数の次数間高調波の電圧,電流を精度よく抽出し
て検出することができる。
【0094】なお、次数間高調波の注入周波数finj
は、finj=kf−mf0 から決定してもよい。
【0095】また、単独運転の誤検出防止等を図るた
め、複数次数の次数間高調波を注入し、その検出結果か
ら多数決の原理で単独運転を検出する場合等は、各次数
間高調波は、数1の式にしたがって、k次高調波とk+
1次高調波との間にf0 の間隔でその注入周波数(注入
次数)を決定することが好ましい。
【0096】この場合、決定した注入周波数に基づき、
例えばm=1,3,5,7,…又はm=2,4,6,…
の1つおきの注入周波数の次数間高調波を実際に注入
し、それぞれの上,下隣りに非注入次数の空きの次数間
高調波が存在するようにし、注入次数の次数間高調波に
ついて前記実施の形態と同様にして個別に系統成分の補
正を施してもよいが、つぎに説明するように一括して補
正を施すようにしてもよい。
【0097】すなわち、数1の式にしたがって例えばm
=m−1,mの連続する複数次数の次数間高調波を決定
して電力系統5に注入し、このとき、図3に示すように
それらの上,下のm=m−2,m+1の非注入の次数間
高調波の振幅GCm-2,GCm+1(GSm-2,GS
+1)もゲインGCm-1,GCm(GSm-1,GSm)
と同様にほぼ線形変化することから、振幅GCm-2とG
Cm+1(GSm-2とGSm +1)のいわゆる重み付けの平
均値により、振幅GCm-1,GCm(GSm-1,GS
m)を補間演算する。
【0098】具体的には、つぎの数31〜数32の式の
重み付けの平均により振幅GCm,GSm,GCm+1
GSm+1をk.m−1次,k.m次の誤差成分ΔC
-1,ΔSm-1,ΔCm,ΔSmを補間演算し、成分
(Cm-1,Sm-1,Cm,Sm)を一括して補正すれば
よい。
【0099】
【数29】
【0100】
【数30】
【0101】
【数31】
【0102】
【数32】
【0103】ところで、次数間高調波の周波数を決定す
る場合、k次,k+1次の高調波の影響を極力少なくす
るため、極力それらの中間の周波数に決定することが好
ましく、具体的にはk=2,f=60Hzの場合、次数間
高調波の注入周波数は、120Hz(=2×60Hz)と1
80Hz(=3×60Hz)との中間の150Hzに近い周波
数に決定することが好ましい。
【0104】そして、本発明は分散型電源の単独運転の
検出における次数間高調波検出だけでなく、例えば次数
間高調波を用いて電力系統の特性を検出する際の次数間
高調波の検出にも適用することができるのは勿論であ
る。
【0105】
【発明の効果】本発明は、以下に記載する効果を奏す
る。まず、請求項1の場合は、電力系統5に注入した次
数間高調波をPLL回路等を用いることなく、一定サン
プリング周波数の固定周波数サンプリング方式で安価か
つ簡単にサンプリングすることができる。
【0106】さらに、サンプリングデータのフーリエ変
換で抽出された次数間高調波の系統基本波及び高調波の
影響に基づく誤差が簡単な平均値の補間演算と減算によ
り除去され、系統周波数とサンプリング周波数との非同
期に起因した誤差を補正して精度の高い次数間高調波の
検出を行うことができる。
【0107】したがって、固定周波数サンプリング方式
で精度の高い次数間高調波の検出を行うことができ、分
散型電源の単独運転検出等に好適な次数間高調波検出方
法を提供することができる。
【0108】また、請求項2の場合、注入する次数間高
調波の周波数を簡単な式の演算から決定して請求項1と
同様の効果を得ることができ、極めて実用的である。
【図面の簡単な説明】
【図1】本発明の実施の1形態の単線結線図である。
【図2】図1の演算処理装置の補間演算の1例を説明す
る波形図である。
【図3】図1の演算処理装置の補間演算の他の例を説明
する波形図である。
【符号の説明】
5 電力系統 16 注入装置 22 サンプル・ホールド回路 23 タイミング指令部 24 A/D変換部 25 演算処理部

Claims (2)

    【特許請求の範囲】
  1. 【請求項1】 電力系統の電圧又は/及び電流の計測信
    号をサンプリングし、サンプリングデータのフーリエ変
    換により、前記電力系統に注入された系統基本波の非整
    数倍周波数の注入次数の次数間高調波を抽出して検出す
    る次数間高調波検出方法において、 前記電力系統の電圧又は/及び電流をサンプリング周波
    数一定の固定周波数サンプリング方式でサンプリング
    し、 サンプリングデータのフーリエ変換によるフィルタ演算
    により前記注入次数の次数間高調波及び前記注入次数の
    上,下の非注入次数の次数間高調波の成分を抽出し、 前記両非注入次数の次数間高調波の抽出成分の平均から
    前記注入次数の次数間高調波の抽出成分に含まれた前記
    系統基本波及び整数倍周波数の高調波の系統成分に基づ
    く誤差を補間演算し、 前記注入次数の次数間高調波の抽出成分から前記誤差を
    減算して除去し、 前記注入次数の次数間高調波を検出することを特徴とす
    る次数間高調波検出方法。
  2. 【請求項2】 注入次数の次数間高調波がk次高調波と
    k+1次高調波(kは2以上の整数)との間の次数間高
    調波であり、前記次数間高調波の注入周波数を、kf±
    mf0 ,(fは系統基本波周波数,mは1,2,…,n
    −1の整数nは2以上の整数,f0は次数間高調波の注
    入間隔を示し、f0=f/n)から決定することを特徴
    とする請求項1記載の次数間高調波検出方法。
JP2000226365A 2000-03-13 2000-07-27 次数間高調波検出方法 Expired - Lifetime JP3788212B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2000226365A JP3788212B2 (ja) 2000-07-27 2000-07-27 次数間高調波検出方法
US09/803,945 US6545885B2 (en) 2000-03-13 2001-03-13 Isolated operation prevention device for distributed power supply and interharmonic detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000226365A JP3788212B2 (ja) 2000-07-27 2000-07-27 次数間高調波検出方法

Publications (2)

Publication Number Publication Date
JP2002040067A true JP2002040067A (ja) 2002-02-06
JP3788212B2 JP3788212B2 (ja) 2006-06-21

Family

ID=18719999

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000226365A Expired - Lifetime JP3788212B2 (ja) 2000-03-13 2000-07-27 次数間高調波検出方法

Country Status (1)

Country Link
JP (1) JP3788212B2 (ja)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006266684A (ja) * 2005-03-22 2006-10-05 Meidensha Corp 高調波検出装置
JP2010136548A (ja) * 2008-12-05 2010-06-17 Kansai Electric Power Co Inc:The 分散電源の単独運転検出装置
CN102520245A (zh) * 2011-09-28 2012-06-27 天津大学 基于三次样条插值波形重构的微网谐波及间谐波分析方法
CN103424622A (zh) * 2013-06-21 2013-12-04 国家电网公司 一种分数谐波源的检测和判断方法
CN104181390A (zh) * 2014-08-18 2014-12-03 信元瑞电气有限公司 一种基于三相四线制零序电流分离的谐波检测方法
CN105891632A (zh) * 2016-04-07 2016-08-24 中国矿业大学 一种基于二阶低通陷波器的同步信号检测方法
CN106872777A (zh) * 2017-02-15 2017-06-20 国网江苏省电力公司苏州供电公司 一种谐波和间谐波分离分析方法
CN107064634A (zh) * 2017-04-18 2017-08-18 辽宁大学 电力系统谐波的检测方法
CN111190042A (zh) * 2019-12-30 2020-05-22 中国电力科学研究院有限公司 一种智能传感终端及对电网信号全频段测量的方法
CN111257637A (zh) * 2020-02-21 2020-06-09 深圳市中电电力技术股份有限公司 超高次谐波测量方法、系统
CN112667962A (zh) * 2020-11-11 2021-04-16 福州大学 一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法
CN113315126A (zh) * 2021-05-31 2021-08-27 华中科技大学 一种有源电力滤波器指定次谐波抑制二次采样方法及系统
CN113517724A (zh) * 2021-04-09 2021-10-19 国网山西省电力公司电力科学研究院 一种交直流混合微电网直流侧电压纹波的抑制方法
CN114172209A (zh) * 2021-11-29 2022-03-11 国网福建省电力有限公司电力科学研究院 基于滑动窗滤波和容量自适应分配的冲击功率检测方法
US11378604B2 (en) * 2018-02-02 2022-07-05 Sankyo Seisakusho Co. Method and device for detecting positional change amount due to movement of moving body

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006266684A (ja) * 2005-03-22 2006-10-05 Meidensha Corp 高調波検出装置
JP2010136548A (ja) * 2008-12-05 2010-06-17 Kansai Electric Power Co Inc:The 分散電源の単独運転検出装置
CN102520245A (zh) * 2011-09-28 2012-06-27 天津大学 基于三次样条插值波形重构的微网谐波及间谐波分析方法
CN103424622A (zh) * 2013-06-21 2013-12-04 国家电网公司 一种分数谐波源的检测和判断方法
CN104181390A (zh) * 2014-08-18 2014-12-03 信元瑞电气有限公司 一种基于三相四线制零序电流分离的谐波检测方法
CN105891632A (zh) * 2016-04-07 2016-08-24 中国矿业大学 一种基于二阶低通陷波器的同步信号检测方法
CN105891632B (zh) * 2016-04-07 2018-11-02 中国矿业大学 一种基于二阶低通陷波器的同步信号检测方法
CN106872777A (zh) * 2017-02-15 2017-06-20 国网江苏省电力公司苏州供电公司 一种谐波和间谐波分离分析方法
CN107064634A (zh) * 2017-04-18 2017-08-18 辽宁大学 电力系统谐波的检测方法
CN107064634B (zh) * 2017-04-18 2019-07-02 辽宁大学 电力系统谐波的检测方法
US11378604B2 (en) * 2018-02-02 2022-07-05 Sankyo Seisakusho Co. Method and device for detecting positional change amount due to movement of moving body
CN111190042A (zh) * 2019-12-30 2020-05-22 中国电力科学研究院有限公司 一种智能传感终端及对电网信号全频段测量的方法
CN111190042B (zh) * 2019-12-30 2022-06-21 中国电力科学研究院有限公司 一种智能传感终端及对电网信号全频段测量的方法
CN111257637B (zh) * 2020-02-21 2022-05-20 深圳市中电电力技术股份有限公司 超高次谐波测量方法、系统
CN111257637A (zh) * 2020-02-21 2020-06-09 深圳市中电电力技术股份有限公司 超高次谐波测量方法、系统
CN112667962A (zh) * 2020-11-11 2021-04-16 福州大学 一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法
CN112667962B (zh) * 2020-11-11 2023-02-14 福州大学 一种电网谐波潮流相量矩阵计算方法
CN113517724A (zh) * 2021-04-09 2021-10-19 国网山西省电力公司电力科学研究院 一种交直流混合微电网直流侧电压纹波的抑制方法
CN113517724B (zh) * 2021-04-09 2022-11-25 国网山西省电力公司电力科学研究院 一种交直流混合微电网直流侧电压纹波的抑制方法
CN113315126A (zh) * 2021-05-31 2021-08-27 华中科技大学 一种有源电力滤波器指定次谐波抑制二次采样方法及系统
CN113315126B (zh) * 2021-05-31 2022-07-05 华中科技大学 一种有源电力滤波器指定次谐波抑制二次采样方法及系统
CN114172209A (zh) * 2021-11-29 2022-03-11 国网福建省电力有限公司电力科学研究院 基于滑动窗滤波和容量自适应分配的冲击功率检测方法
CN114172209B (zh) * 2021-11-29 2023-08-01 国网福建省电力有限公司电力科学研究院 基于滑动窗滤波和容量自适应分配的冲击功率检测方法

Also Published As

Publication number Publication date
JP3788212B2 (ja) 2006-06-21

Similar Documents

Publication Publication Date Title
US6545885B2 (en) Isolated operation prevention device for distributed power supply and interharmonic detection method
JP3788212B2 (ja) 次数間高調波検出方法
Moore et al. A new numeric technique for high-speed evaluation of power system frequency
Chakir et al. Extended C37. 118.1 PMU algorithms for joint tracking of fundamental and harmonic phasors in stressed power systems and microgrids
US7180300B2 (en) System and method of locating ground fault in electrical power distribution system
US7328114B2 (en) Methods and systems for measuring a rate of change of frequency
Verma et al. An improved hybrid prefiltered open-loop algorithm for three-phase grid synchronization
US4261038A (en) Protection of electrical power supply systems
US6326796B1 (en) Harmonic measuring method and a current injection device for harmonic measurement
US5832413A (en) Generator protection system and method for phasor estimation and frequency tracking during frequency ramping
US9423434B2 (en) Method and system for measuring electrical quantity in electrical network
US10613127B2 (en) Determining the frequency of an alternating signal
JP2012163384A (ja) 短絡容量計測システム
JP3367371B2 (ja) 分散型電源の単独運転防止装置
Elsahwi et al. A digital frequency adaptive synchronization unit for on-and off-grid systems
US11043803B2 (en) Reference signal generating method for distance and directional protection elements
Grando et al. A method for synchronized harmonic phasor measurement based on hardware‐enhanced Fast Fourier Transform
JP5518274B1 (ja) 保護リレー装置
JP4316103B2 (ja) 故障点標定装置
Zhang et al. An innovative timestamp-based convolution integral method in synchrophasor estimation within digital relays
JP3896051B2 (ja) 分散型電源の単独運転検出方法
JP2943355B2 (ja) 交流電圧の差電圧演算方法
Kusic et al. Measurement of power system phase differences by means of GPS timing
JP2005304227A (ja) 分散電源の単独運転検出装置
JP2003244851A (ja) 分散電源の単独運転検出方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040421

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20051219

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060307

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060320

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 3788212

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100407

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100407

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110407

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120407

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130407

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140407

Year of fee payment: 8

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term