JP2002058163A - 電力系統の周波数特性同定方法 - Google Patents
電力系統の周波数特性同定方法Info
- Publication number
- JP2002058163A JP2002058163A JP2000238795A JP2000238795A JP2002058163A JP 2002058163 A JP2002058163 A JP 2002058163A JP 2000238795 A JP2000238795 A JP 2000238795A JP 2000238795 A JP2000238795 A JP 2000238795A JP 2002058163 A JP2002058163 A JP 2002058163A
- Authority
- JP
- Japan
- Prior art keywords
- power system
- identification
- value
- fourier
- frequency characteristic
- 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
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
(57)【要約】
【課題】離散値系に好適な電力系統の周波数特性の同定
方法を提供する。 【解決手段】同定信号生成回路２１からの同定用信号の
指令値（ｉC *）をｘ［ｎ］とし、電力系統２の電流（ｉ
Sh）をｙ［ｎ］とし、このｘ［ｎ］とｙ［ｎ］とをそれ
ぞれフーリエ変換し、さらに、これらの相互相関関数の
フーリエ変換した値ＤＦＴ［ｒxy［ｍ］］に含まれるｋ
次高調波の振幅を、指令値（ｉc *）のｋ次の振幅の２乗
で割った値をｋ次のゲイン−周波数特性の同定値とす
る。
方法を提供する。 【解決手段】同定信号生成回路２１からの同定用信号の
指令値（ｉC *）をｘ［ｎ］とし、電力系統２の電流（ｉ
Sh）をｙ［ｎ］とし、このｘ［ｎ］とｙ［ｎ］とをそれ
ぞれフーリエ変換し、さらに、これらの相互相関関数の
フーリエ変換した値ＤＦＴ［ｒxy［ｍ］］に含まれるｋ
次高調波の振幅を、指令値（ｉc *）のｋ次の振幅の２乗
で割った値をｋ次のゲイン−周波数特性の同定値とす
る。
Description
【０００１】
【発明の属する技術分野】この発明は、電力系統と連系
運転される半導体電力変換装置などを介して、該電力系
統の周波数特性の同定方法に関する。
運転される半導体電力変換装置などを介して、該電力系
統の周波数特性の同定方法に関する。
【０００２】
【従来の技術】交流電動機駆動用インバータに代表され
る半導体電力変換装置の普及に伴って、電力系統の高調
波歪みの増大など電力環境の悪化が問題となり、この悪
化を補正するために該電力系統と連系運転される半導体
電力変換装置として、例えばアクティブフィルタなどを
設置することが行われている。このアクティブフィルタ
において、上述の電力環境の悪化をより良く補正するた
めには電力系統の周波数特性を同定する必要がある。
る半導体電力変換装置の普及に伴って、電力系統の高調
波歪みの増大など電力環境の悪化が問題となり、この悪
化を補正するために該電力系統と連系運転される半導体
電力変換装置として、例えばアクティブフィルタなどを
設置することが行われている。このアクティブフィルタ
において、上述の電力環境の悪化をより良く補正するた
めには電力系統の周波数特性を同定する必要がある。
【０００３】この同定方法として、周知のスペクトル密
度解析方法では、発振周波数の異なる複数の正弦波信号
を合成してなる同定用信号〔ｘ（ｔ）〕を電力系統に注
入し、この注入による電力系統の応答〔ｙ（ｔ）〕を求
め、前記〔ｘ（ｔ）〕と〔ｙ（ｔ）〕との間の相互相関
関数〔Ｒxy（τ）〕をフーリエ変換して得られる相互ス
ペクトル密度を、前記〔ｘ（ｔ）〕の自己相関関数〔Ｒ
xx（τ）〕をフーリエ変換して得られる自己スペクトル
密度で除算演算した値をこの電力系統の周波数特性の同
定値としていた。
度解析方法では、発振周波数の異なる複数の正弦波信号
を合成してなる同定用信号〔ｘ（ｔ）〕を電力系統に注
入し、この注入による電力系統の応答〔ｙ（ｔ）〕を求
め、前記〔ｘ（ｔ）〕と〔ｙ（ｔ）〕との間の相互相関
関数〔Ｒxy（τ）〕をフーリエ変換して得られる相互ス
ペクトル密度を、前記〔ｘ（ｔ）〕の自己相関関数〔Ｒ
xx（τ）〕をフーリエ変換して得られる自己スペクトル
密度で除算演算した値をこの電力系統の周波数特性の同
定値としていた。
【０００４】上述のスペクトル密度解析方法は、ノイズ
に比較的強いため、微弱な同定用信号で前記周波数特性
を同定でき、また、フーリエ変換を用いて複数の前記発
振周波数点における特性を一度に求められるので短時間
に前記周波数特性を同定できるが、前記電力系統の負荷
として高調波電流発生源を有する負荷が接続されている
場合には、この高調波電流発生源の高調波電流が外乱と
なり、このスペクトル密度解析方法で得られる前記同定
値には前記外乱による誤差が含まれている。
に比較的強いため、微弱な同定用信号で前記周波数特性
を同定でき、また、フーリエ変換を用いて複数の前記発
振周波数点における特性を一度に求められるので短時間
に前記周波数特性を同定できるが、前記電力系統の負荷
として高調波電流発生源を有する負荷が接続されている
場合には、この高調波電流発生源の高調波電流が外乱と
なり、このスペクトル密度解析方法で得られる前記同定
値には前記外乱による誤差が含まれている。
【０００５】上述の周知のスペクトル密度解析方法にお
ける外乱誤差を解消するためには、同一出願人による平
成１１年特許願第２２９７５１号に記載の同定方法があ
り、この同定方法では、発振周波数の異なる複数の正弦
波信号を合成してなる同定用信号〔ｘ（ｔ）〕を用い、
この同定信号を構成する個々の正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ：複数）ずつ変位させて
電力系統に注入し、この注入による電力系統の前記ｎ個
の応答〔ｙ（ｔ）〕を求め、前記それぞれの〔ｘ
（ｔ）〕と〔ｙ（ｔ）〕との間の相互相関関数〔Ｒ
xy（τ）〕のベクトル的加算平均値を求めることにより
先述の高調波電流成分を除去し、この加算平均値をフー
リエ変換して得られる相互スペクトル密度を、前記〔ｘ
（ｔ）〕の自己相関関数〔Ｒxx（τ）〕をフーリエ変換
して得られる自己スペクトル密度で除算演算した値をこ
の電力系統の周波数特性の同定値としている。
ける外乱誤差を解消するためには、同一出願人による平
成１１年特許願第２２９７５１号に記載の同定方法があ
り、この同定方法では、発振周波数の異なる複数の正弦
波信号を合成してなる同定用信号〔ｘ（ｔ）〕を用い、
この同定信号を構成する個々の正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ：複数）ずつ変位させて
電力系統に注入し、この注入による電力系統の前記ｎ個
の応答〔ｙ（ｔ）〕を求め、前記それぞれの〔ｘ
（ｔ）〕と〔ｙ（ｔ）〕との間の相互相関関数〔Ｒ
xy（τ）〕のベクトル的加算平均値を求めることにより
先述の高調波電流成分を除去し、この加算平均値をフー
リエ変換して得られる相互スペクトル密度を、前記〔ｘ
（ｔ）〕の自己相関関数〔Ｒxx（τ）〕をフーリエ変換
して得られる自己スペクトル密度で除算演算した値をこ
の電力系統の周波数特性の同定値としている。
【０００６】
【発明が解決しようとする課題】上述の同一出願人によ
る平成１１年特許願第２２９７５１号に記載の同定方法
では前記相関関数，フーリエ変換など比較的計算量が多
い演算を行う必要があり、その結果、電力系統と連系運
転され、デジタル制御，演算を行う半導体電力変換装置
において、この半導体電力変換装置を流用して該電力系
統の周波数特性を同定するときには、その演算時間が長
くなるという難点があった。
る平成１１年特許願第２２９７５１号に記載の同定方法
では前記相関関数，フーリエ変換など比較的計算量が多
い演算を行う必要があり、その結果、電力系統と連系運
転され、デジタル制御，演算を行う半導体電力変換装置
において、この半導体電力変換装置を流用して該電力系
統の周波数特性を同定するときには、その演算時間が長
くなるという難点があった。
【０００７】この発明の目的は、上記難点を解消した電
力系統の周波数特性同定方法を提供することにある。
力系統の周波数特性同定方法を提供することにある。
【０００８】
【課題を解決するための手段】この第１の発明は、発振
周波数が単一の正弦波信号からなる同定用信号又は発振
周波数が異なる複数の正弦波信号を合成してなる同定用
信号を生成し、この同定用信号を電力系統に注入したと
きの該電力系統の応答を求め、求めた応答及び前記同定
用信号を個別にフーリエ変換し、さらに、これらのフー
リエ変換した値の相互相関関数に対するフーリエ変換を
行い、この相互相関関数のフーリエ変換した値から導出
される前記発振周波数毎の振幅を、この発振周波数に対
応する前記正弦波信号の振幅の２乗で除算演算した値を
前記周波数特性の同定値とする電力系統の周波数特性同
定方法を行う。
周波数が単一の正弦波信号からなる同定用信号又は発振
周波数が異なる複数の正弦波信号を合成してなる同定用
信号を生成し、この同定用信号を電力系統に注入したと
きの該電力系統の応答を求め、求めた応答及び前記同定
用信号を個別にフーリエ変換し、さらに、これらのフー
リエ変換した値の相互相関関数に対するフーリエ変換を
行い、この相互相関関数のフーリエ変換した値から導出
される前記発振周波数毎の振幅を、この発振周波数に対
応する前記正弦波信号の振幅の２乗で除算演算した値を
前記周波数特性の同定値とする電力系統の周波数特性同
定方法を行う。
【０００９】第２の発明は、発振周波数が単一の正弦波
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を電力系統に注入したときの該電力系統の応答
を求め、求めた応答及び前記同定用信号を個別にフーリ
エ変換し、さらに、これらのフーリエ変換した値の相互
相関関数に対するフーリエ変換を行い、前記同定用信号
をフーリエ変換した値の自己相関関数に対するフーリエ
変換を行い、前記相互相関関数のフーリエ変換した値か
ら導出される前記発振周波数毎の振幅を前記自己相関関
数のフーリエ変換した値から導出される該発振周波数の
振幅で除算演算した値を、前記周波数特性の同定値とす
る電力系統の周波数特性同定方法を行う。
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を電力系統に注入したときの該電力系統の応答
を求め、求めた応答及び前記同定用信号を個別にフーリ
エ変換し、さらに、これらのフーリエ変換した値の相互
相関関数に対するフーリエ変換を行い、前記同定用信号
をフーリエ変換した値の自己相関関数に対するフーリエ
変換を行い、前記相互相関関数のフーリエ変換した値か
ら導出される前記発振周波数毎の振幅を前記自己相関関
数のフーリエ変換した値から導出される該発振周波数の
振幅で除算演算した値を、前記周波数特性の同定値とす
る電力系統の周波数特性同定方法を行う。
【００１０】第３の発明は、発振周波数が単一の正弦波
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を構成するそれぞれの正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位させつ
つ、電力系統に注入したときの該電力系統の応答を求
め、前記ｎ個の求めた応答及び同定用信号を個別にフー
リエ変換し、さらに、これらのフーリエ変換した値の相
互相関関数に対するフーリエ変換のベクトル的加算平均
値を求め、このベクトル的加算平均値から導出される前
記発振周波数毎の振幅をこの発振周波数に対応する前記
正弦波信号の振幅の２乗で除算演算した値を、前記周波
数特性の同定値とする電力系統の周波数特性同定方法を
行う。
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を構成するそれぞれの正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位させつ
つ、電力系統に注入したときの該電力系統の応答を求
め、前記ｎ個の求めた応答及び同定用信号を個別にフー
リエ変換し、さらに、これらのフーリエ変換した値の相
互相関関数に対するフーリエ変換のベクトル的加算平均
値を求め、このベクトル的加算平均値から導出される前
記発振周波数毎の振幅をこの発振周波数に対応する前記
正弦波信号の振幅の２乗で除算演算した値を、前記周波
数特性の同定値とする電力系統の周波数特性同定方法を
行う。
【００１１】第４の発明は、発振周波数が単一の正弦波
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を構成するそれぞれの正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位させつ
つ、電力系統に注入したときの該電力系統の応答を求
め、前記ｎ個の求めた応答及び同定用信号を個別にフー
リエ変換し、さらに、これらのフーリエ変換した値の相
互相関関数に対するフーリエ変換のベクトル的加算平均
値を求め、前記ｎ個の同定用信号をフーリエ変換した値
の自己相関関数に対するフーリエ変換のベクトル的加算
平均値を求め、前記相互相関関数のフーリエ変換のベク
トル的加算平均値から導出される前記発振周波数毎の振
幅を、前記自己相関関数のフーリエ変換した値のうちの
１つ、若しくわ前記自己相関関数のフーリエ変換した値
のベクトル的加算平均値導出される該発振周波数の振幅
で除算演算した値を前記周波数特性の同定値とする電力
系統の周波数特性同定方法を行う。
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を構成するそれぞれの正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位させつ
つ、電力系統に注入したときの該電力系統の応答を求
め、前記ｎ個の求めた応答及び同定用信号を個別にフー
リエ変換し、さらに、これらのフーリエ変換した値の相
互相関関数に対するフーリエ変換のベクトル的加算平均
値を求め、前記ｎ個の同定用信号をフーリエ変換した値
の自己相関関数に対するフーリエ変換のベクトル的加算
平均値を求め、前記相互相関関数のフーリエ変換のベク
トル的加算平均値から導出される前記発振周波数毎の振
幅を、前記自己相関関数のフーリエ変換した値のうちの
１つ、若しくわ前記自己相関関数のフーリエ変換した値
のベクトル的加算平均値導出される該発振周波数の振幅
で除算演算した値を前記周波数特性の同定値とする電力
系統の周波数特性同定方法を行う。
【００１２】第５の発明は、発振周波数が単一の正弦波
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を電力系統に注入したときの該電力系統の応答
を求め、求めた応答及び前記同定用信号を個別にフーリ
エ変換し、さらに、これらのフーリエ変換した値の相互
相関関数に対するフーリエ変換を行い、この相互相関関
数をフーリエ変換した値から導出される前記発振周波数
毎の位相量を前記周波数特性の同定値とする電力系統の
周波数特性同定方法を行う。
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を電力系統に注入したときの該電力系統の応答
を求め、求めた応答及び前記同定用信号を個別にフーリ
エ変換し、さらに、これらのフーリエ変換した値の相互
相関関数に対するフーリエ変換を行い、この相互相関関
数をフーリエ変換した値から導出される前記発振周波数
毎の位相量を前記周波数特性の同定値とする電力系統の
周波数特性同定方法を行う。
【００１３】第６の発明は、発振周波数が単一の正弦波
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を構成するそれぞれの正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位させつ
つ、電力系統に注入したときの該電力系統の応答を求
め、前記ｎ個の求めた応答及び同定用信号を個別にフー
リエ変換し、さらに、これらのフーリエ変換した値の相
互相関関数に対するフーリエ変換のベクトル的加算平均
値を求め、このベクトル的加算平均値から導出される前
記発振周波数毎の位相量を前記周波数特性の同定値とす
る電力系統の周波数特性同定方法を行う。
信号からなる同定用信号又は発振周波数が異なる複数の
正弦波信号を合成してなる同定用信号を生成し、この同
定用信号を構成するそれぞれの正弦波信号の位相を前回
より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位させつ
つ、電力系統に注入したときの該電力系統の応答を求
め、前記ｎ個の求めた応答及び同定用信号を個別にフー
リエ変換し、さらに、これらのフーリエ変換した値の相
互相関関数に対するフーリエ変換のベクトル的加算平均
値を求め、このベクトル的加算平均値から導出される前
記発振周波数毎の位相量を前記周波数特性の同定値とす
る電力系統の周波数特性同定方法を行う。
【００１４】第７の発明は、前記第１〜第６の発明の電
力系統の周波数特性同定方法において、前記同定用信号
は、電力系統と連系運転する半導体電力変換装置に付加
的に指令される同定用信号指令値により該半導体電力変
換装置で生成することを特徴とする。
力系統の周波数特性同定方法において、前記同定用信号
は、電力系統と連系運転する半導体電力変換装置に付加
的に指令される同定用信号指令値により該半導体電力変
換装置で生成することを特徴とする。
【００１５】第８の発明は、前記第７の発明の電力系統
の周波数特性同定方法において、前記同定用信号を前記
同定用信号指令値とし、前記応答を上位の電力系統と前
記電力系統との間の電流又は電圧として前記周波数特性
の同定値を演算することを特徴とする。
の周波数特性同定方法において、前記同定用信号を前記
同定用信号指令値とし、前記応答を上位の電力系統と前
記電力系統との間の電流又は電圧として前記周波数特性
の同定値を演算することを特徴とする。
【００１６】第９の発明は、前記第７の発明の電力系統
の周波数特性同定方法において、前記同定用信号を前記
半導体電力変換装置の出力電流又は出力電圧とし、前記
応答を上位の電力系統と前記電力系統との間の電流又は
電圧として前記周波数特性の同定値を演算することを特
徴とする。
の周波数特性同定方法において、前記同定用信号を前記
半導体電力変換装置の出力電流又は出力電圧とし、前記
応答を上位の電力系統と前記電力系統との間の電流又は
電圧として前記周波数特性の同定値を演算することを特
徴とする。
【００１７】第１０の発明は、前記第１〜第９の発明の
電力系統の周波数特性同定方法において、前記同定用信
号を構成する正弦波信号の発振周波数を、前記電力系統
の基本波周波数に比例した値としたことを特徴とする。
電力系統の周波数特性同定方法において、前記同定用信
号を構成する正弦波信号の発振周波数を、前記電力系統
の基本波周波数に比例した値としたことを特徴とする。
【００１８】この発明によれば、相互相関関数のフーリ
エ変換した値を、後述の如く、同定用信号のフーリエ変
換と前記応答のフーリエ変換との積から直接に得るよう
にしているため、前記周波数特性の同定のための演算時
間を短縮することとなり、また、デジタル演算の際に好
適な高速フーリエ変換のアルゴリズムを用いることがで
きるので、さらに演算時間を短縮することができる。
エ変換した値を、後述の如く、同定用信号のフーリエ変
換と前記応答のフーリエ変換との積から直接に得るよう
にしているため、前記周波数特性の同定のための演算時
間を短縮することとなり、また、デジタル演算の際に好
適な高速フーリエ変換のアルゴリズムを用いることがで
きるので、さらに演算時間を短縮することができる。
【００１９】また、第３，第４，第６の発明によれば、
前記電力系統の負荷として高調波電流発生源を有する負
荷が接続されている場合には、この高調波電流発生源の
高調波電流に起因する前記同定値の誤差を除去すること
ができる。
前記電力系統の負荷として高調波電流発生源を有する負
荷が接続されている場合には、この高調波電流発生源の
高調波電流に起因する前記同定値の誤差を除去すること
ができる。
【００２０】
【発明の実施の形態】図１は、この発明の電力系統の周
波数特性同定方法の第１の実施の形態を示す構成図であ
り、１は上位の電力系統、１ａは前記電力系統１を一次
側とした変圧器、２は変圧器１ａの二次側から給電さ
れ、この発明の対象とする電力系統、３は前記電力系統
２の負荷、４は電力系統２から負荷３へ流れる電流（ｉ
Sh）を検出するＣＴ、１１は補正値生成回路、１２は切
替スイッチ、１３は電流調節器、１４はインバータ、１
５はインバータ１４が出力する電流（ｉC ）を検出する
ＣＴ、２１は同定用信号生成回路、２２は同定値演算回
路を示す。
波数特性同定方法の第１の実施の形態を示す構成図であ
り、１は上位の電力系統、１ａは前記電力系統１を一次
側とした変圧器、２は変圧器１ａの二次側から給電さ
れ、この発明の対象とする電力系統、３は前記電力系統
２の負荷、４は電力系統２から負荷３へ流れる電流（ｉ
Sh）を検出するＣＴ、１１は補正値生成回路、１２は切
替スイッチ、１３は電流調節器、１４はインバータ、１
５はインバータ１４が出力する電流（ｉC ）を検出する
ＣＴ、２１は同定用信号生成回路、２２は同定値演算回
路を示す。
【００２１】図１において、補正値生成回路１１と、紙
面下側の接点に閉路したときの切替スイッチ１２と、電
流調節器１３と、インバータ１４とでアクティブフィル
タを構成し、補正値生成回路１１では前記ＣＴ４の検出
値から、このアクティブフィルタが電力系統２に対して
補正すべき電流補正値を生成し、電流調節器１３では前
記電流補正値と前記ＣＴ１５の検出値との偏差を零にす
る調節演算を行い、この演算結果をインバータ１４への
出力電圧指令値とすることにより、インバータ１４から
電力系統２へ補正すべき電流を注入することができる。
面下側の接点に閉路したときの切替スイッチ１２と、電
流調節器１３と、インバータ１４とでアクティブフィル
タを構成し、補正値生成回路１１では前記ＣＴ４の検出
値から、このアクティブフィルタが電力系統２に対して
補正すべき電流補正値を生成し、電流調節器１３では前
記電流補正値と前記ＣＴ１５の検出値との偏差を零にす
る調節演算を行い、この演算結果をインバータ１４への
出力電圧指令値とすることにより、インバータ１４から
電力系統２へ補正すべき電流を注入することができる。
【００２２】次に、切替スイッチ１２の接点を紙面上側
に閉路させ、電力系統２の周波数特性を同定するための
動作について説明する。
に閉路させ、電力系統２の周波数特性を同定するための
動作について説明する。
【００２３】図２は図１に示した同定信号生成回路２１
の動作例を説明する波形図であり、紙面最上段から順
に、電力系統２の基本波周波数（紙面では、基本周期と
称している）に対応した５次高調波，７次高調波，１１
次高調波の正弦波信号と、これらを合成した同定用信号
の指令値（ｉC *）とを示している。
の動作例を説明する波形図であり、紙面最上段から順
に、電力系統２の基本波周波数（紙面では、基本周期と
称している）に対応した５次高調波，７次高調波，１１
次高調波の正弦波信号と、これらを合成した同定用信号
の指令値（ｉC *）とを示している。
【００２４】電流調節器１３では前記指令値（ｉC *）と
前記ＣＴ１５の検出値との偏差を零にする調節演算を行
い、この演算結果をインバータ１４への出力電圧指令値
とすることにより、インバータ１４から電力系統２の周
波数特性を同定するための電流を注入することができ
る。
前記ＣＴ１５の検出値との偏差を零にする調節演算を行
い、この演算結果をインバータ１４への出力電圧指令値
とすることにより、インバータ１４から電力系統２の周
波数特性を同定するための電流を注入することができ
る。
【００２５】図３はこの発明の第１の実施例を示し、図
１に示した同定値演算回路２２の動作を説明するフロー
チャートである。
１に示した同定値演算回路２２の動作を説明するフロー
チャートである。
【００２６】すなわち図３において、先ず、同定信号生
成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、係数の取り方などの違いにより複数の定義
があるが、例えば下記数１の式に示す如く、このｘ
［ｎ］を離散値系でフーリエ変換したＸ［ｋ］を求め
る。
成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、係数の取り方などの違いにより複数の定義
があるが、例えば下記数１の式に示す如く、このｘ
［ｎ］を離散値系でフーリエ変換したＸ［ｋ］を求め
る。
【００２７】
【数１】
【００２８】ここで、Ｎはデータ点数である。
【００２９】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、係数の取り方などにより複数の定義がある
が、例えば下記数２の式に示す如く、このｙ［ｎ］を離
散値系でフーリエ変換したＹ［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、係数の取り方などにより複数の定義がある
が、例えば下記数２の式に示す如く、このｙ［ｎ］を離
散値系でフーリエ変換したＹ［ｋ］を求める。
【００３０】
【数２】
【００３１】ここで、Ｎはデータ点数である。
【００３２】次に、上記数１及び数２のフーリエ変換し
た値の相互相関関数のフーリエ変換したＤＦＴ［ｒ
xy［ｍ］］は、下記数３の式で表される。
た値の相互相関関数のフーリエ変換したＤＦＴ［ｒ
xy［ｍ］］は、下記数３の式で表される。
【００３３】
【数３】
【００３４】ここで、Ｘ*［ｋ］はＸ［ｋ］の複素共役
である。
である。
【００３５】このＤＦＴ［ｒxy［ｍ］］に含まれる先述
の５次高調波，７次高調波，１１次高調波それぞれの振
幅を、同定信号生成回路２１が出力する５次高調波，７
次高調波，１１次高調波それぞれの正弦波信号の振幅の
２乗で個々に除算演算したときの商は、前記高調波周波
数における電力系統２のゲイン−周波数特性になる。
の５次高調波，７次高調波，１１次高調波それぞれの振
幅を、同定信号生成回路２１が出力する５次高調波，７
次高調波，１１次高調波それぞれの正弦波信号の振幅の
２乗で個々に除算演算したときの商は、前記高調波周波
数における電力系統２のゲイン−周波数特性になる。
【００３６】すなわち、一般的には、前記ｘ［ｎ］のｋ
次の振幅がＰ［ｋ］ときに、前記ＤＦＴ［ｒxy［ｍ］］
のｋ次の振幅がＲ［ｋ］だったとすると、求めるｋ次の
振幅比はＲ［ｋ］／Ｐ2［ｋ］となる。
次の振幅がＰ［ｋ］ときに、前記ＤＦＴ［ｒxy［ｍ］］
のｋ次の振幅がＲ［ｋ］だったとすると、求めるｋ次の
振幅比はＲ［ｋ］／Ｐ2［ｋ］となる。
【００３７】例えば、前述のアクティブフィルタの補正
値生成回路１１において、単に共振点の周波数が解れば
良いときには前記振幅比の値でよいが、この振幅比の絶
対値が必要な場合にはフーリエ変換の定義式によって
は、さらに補正係数αが必要となり、このときの振幅比
はαＲ［ｋ］／Ｐ2［ｋ］となる。なお、前記αはフー
リエ変換の定義式と、データの点数（Ｎ）とから一意に
決まる定数である。
値生成回路１１において、単に共振点の周波数が解れば
良いときには前記振幅比の値でよいが、この振幅比の絶
対値が必要な場合にはフーリエ変換の定義式によって
は、さらに補正係数αが必要となり、このときの振幅比
はαＲ［ｋ］／Ｐ2［ｋ］となる。なお、前記αはフー
リエ変換の定義式と、データの点数（Ｎ）とから一意に
決まる定数である。
【００３８】図４は、この発明の電力系統の周波数特性
同定方法の第２の実施の形態を示す構成図であり、図１
に示した構成と同一機能を有するものには同一符号を付
して、その説明を省略する。
同定方法の第２の実施の形態を示す構成図であり、図１
に示した構成と同一機能を有するものには同一符号を付
して、その説明を省略する。
【００３９】すなわち図４に示した構成が図１に示した
構成と異なる点は、同定値演算回路２２に代えて、同定
値演算回路２３又は同定値演算回路２４のいずれかを備
えていることである。
構成と異なる点は、同定値演算回路２２に代えて、同定
値演算回路２３又は同定値演算回路２４のいずれかを備
えていることである。
【００４０】図５はこの発明の第２の実施例を示し、図
４に示した同定値演算回路２３の動作を説明するフロー
チャートである。
４に示した同定値演算回路２３の動作を説明するフロー
チャートである。
【００４１】すなわち図５において、先ず、同定信号生
成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記数１の式に示した如く、このｘ［ｎ］
を離散値系でフーリエ変換したＸ［ｋ］を求める。
成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記数１の式に示した如く、このｘ［ｎ］
を離散値系でフーリエ変換したＸ［ｋ］を求める。
【００４２】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
【００４３】次に、前記Ｘ［ｋ］の自己相関関数のフー
リエ変換した値ＤＦＴ［ｒxx［ｍ］］として、｜Ｘ
［ｋ］｜2 を求める。
リエ変換した値ＤＦＴ［ｒxx［ｍ］］として、｜Ｘ
［ｋ］｜2 を求める。
【００４４】また、前記ｘ［ｎ］とｙ［ｎ］との相互相
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
【００４５】このＤＦＴ［ｒxy［ｍ］］に含まれる先述
の５次高調波，７次高調波，１１次高調波それぞれの振
幅を、前記｜Ｘ［ｋ］｜2 に含まれる５次高調波，７次
高調波，１１次高調波それぞれの振幅で個々に除算演算
したときの商は、前記高調波周波数における電力系統２
のゲイン−周波数特性になる。
の５次高調波，７次高調波，１１次高調波それぞれの振
幅を、前記｜Ｘ［ｋ］｜2 に含まれる５次高調波，７次
高調波，１１次高調波それぞれの振幅で個々に除算演算
したときの商は、前記高調波周波数における電力系統２
のゲイン−周波数特性になる。
【００４６】すなわち、一般的には、前記ＤＦＴ［ｒxx
［ｍ］］のｋ次の振幅がＱ［ｋ］ときに、前記ＤＦＴ
［ｒxy［ｍ］］のｋ次の振幅がＲ［ｋ］だったとする
と、求めるｋ次の振幅比はＲ［ｋ］／Ｑ［ｋ］となる。
［ｍ］］のｋ次の振幅がＱ［ｋ］ときに、前記ＤＦＴ
［ｒxy［ｍ］］のｋ次の振幅がＲ［ｋ］だったとする
と、求めるｋ次の振幅比はＲ［ｋ］／Ｑ［ｋ］となる。
【００４７】図６はこの発明の第３の実施例を示し、図
４に示した同定値演算回路２４の動作を説明するフロー
チャートである。
４に示した同定値演算回路２４の動作を説明するフロー
チャートである。
【００４８】すなわち図６において、先ず、同定信号生
成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記数１の式に示した如く、このｘ［ｎ］
を離散値系でフーリエ変換したＸ［ｋ］を求める。
成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記数１の式に示した如く、このｘ［ｎ］
を離散値系でフーリエ変換したＸ［ｋ］を求める。
【００４９】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
【００５０】次に、前記ｘ［ｎ］とｙ［ｎ］との相互相
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
【００５１】このＤＦＴ［ｒxy［ｍ］］における各次の
位相が、電力系統２の位相−周波数特性となる。
位相が、電力系統２の位相−周波数特性となる。
【００５２】図７は、この発明の電力系統の周波数特性
同定方法の第３の実施の形態を示す構成図であり、図１
に示した構成と同一機能を有するものには同一符号を付
して、その説明を省略する。
同定方法の第３の実施の形態を示す構成図であり、図１
に示した構成と同一機能を有するものには同一符号を付
して、その説明を省略する。
【００５３】すなわち図７に示した構成が図１に示した
構成と異なる点は、同定値演算回路２２に代えて、同定
値演算回路２５又は同定値演算回路２６のいずれかを備
えていることである。
構成と異なる点は、同定値演算回路２２に代えて、同定
値演算回路２５又は同定値演算回路２６のいずれかを備
えていることである。
【００５４】図８はこの発明の第４の実施例を示し、図
７に示した同定値演算回路２５の動作を説明するフロー
チャートである。
７に示した同定値演算回路２５の動作を説明するフロー
チャートである。
【００５５】すなわち図８において、先ず、インバータ
１４から注入され電流（ｉC ）のＣＴ１５を介した検出
値をｘ［ｎ］とし、前記数１の式に示した如く、このｘ
［ｎ］を離散値系でフーリエ変換したＸ［ｋ］を求め
る。
１４から注入され電流（ｉC ）のＣＴ１５を介した検出
値をｘ［ｎ］とし、前記数１の式に示した如く、このｘ
［ｎ］を離散値系でフーリエ変換したＸ［ｋ］を求め
る。
【００５６】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
【００５７】次に、前記Ｘ［ｋ］の自己相関関数のフー
リエ変換した値ＤＦＴ［ｒxx［ｍ］］として、｜Ｘ
［ｋ］｜2 を求める。
リエ変換した値ＤＦＴ［ｒxx［ｍ］］として、｜Ｘ
［ｋ］｜2 を求める。
【００５８】また、前記Ｘ［ｋ］とＹ［ｋ］との相互相
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
【００５９】このＤＦＴ［ｒxy［ｍ］］に含まれる先述
の５次高調波，７次高調波，１１次高調波それぞれの振
幅を、前記｜Ｘ［ｋ］｜2 に含まれる５次高調波，７次
高調波，１１次高調波それぞれの振幅で個々に除算演算
したときの商は、前記高調波周波数における電力系統２
のゲイン−周波数特性になる。
の５次高調波，７次高調波，１１次高調波それぞれの振
幅を、前記｜Ｘ［ｋ］｜2 に含まれる５次高調波，７次
高調波，１１次高調波それぞれの振幅で個々に除算演算
したときの商は、前記高調波周波数における電力系統２
のゲイン−周波数特性になる。
【００６０】すなわち、一般的には、前記ＤＦＴ［ｒxx
［ｍ］］のｋ次の振幅がＱ［ｋ］ときに、前記ＤＦＴ
［ｒxy［ｍ］］のｋ次の振幅がＲ［ｋ］だったとする
と、求めるｋ次の振幅比はＲ［ｋ］／Ｑ［ｋ］となる。
［ｍ］］のｋ次の振幅がＱ［ｋ］ときに、前記ＤＦＴ
［ｒxy［ｍ］］のｋ次の振幅がＲ［ｋ］だったとする
と、求めるｋ次の振幅比はＲ［ｋ］／Ｑ［ｋ］となる。
【００６１】図９はこの発明の第５の実施例を示し、図
７に示した同定値演算回路２６の動作を説明するフロー
チャートである。
７に示した同定値演算回路２６の動作を説明するフロー
チャートである。
【００６２】すなわち図９において、先ず、インバータ
１４から注入された電流（ｉC ）のＣＴ１５を介した検
出値をｘ［ｎ］とし、前記数１の式に示した如く、この
ｘ［ｎ］を離散値系でフーリエ変換したＸ［ｋ］を求め
る。
１４から注入された電流（ｉC ）のＣＴ１５を介した検
出値をｘ［ｎ］とし、前記数１の式に示した如く、この
ｘ［ｎ］を離散値系でフーリエ変換したＸ［ｋ］を求め
る。
【００６３】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記数２の式に示した如く、このｙ［ｎ］
を離散値系でフーリエ変換したＹ［ｋ］を求める。
【００６４】次に、前記ｘ［ｎ］とｙ［ｎ］との相互相
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
関関数のフーリエ変換した値として、前記数３の式に示
した如く、ＤＦＴ［ｒxy［ｍ］］を求める。
【００６５】このＤＦＴ［ｒxy［ｍ］］における各次の
位相が、電力系統２の位相−周波数特性となる。
位相が、電力系統２の位相−周波数特性となる。
【００６６】図１０は、この発明の電力系統の周波数特
性同定方法の第４の実施の形態を示す構成図であり、図
１に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
性同定方法の第４の実施の形態を示す構成図であり、図
１に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
【００６７】すなわち図１０に示した構成が図１に示し
た構成と異なる点は、同定用信号生成回路２１に代えて
同定用信号生成回路２１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
た構成と異なる点は、同定用信号生成回路２１に代えて
同定用信号生成回路２１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
【００６８】このＰＬＬ（フェーズ・ロックド・ルー
プ）回路２７と補助変圧器５とにより、同定用信号生成
回路２１ａが出力する同定信号の指令値（ｉC *）の位相
を電力系統２の基本波の位相に同期させているが、電力
系統２の周波数特性を同定させるための動作は、図１に
示した実施の形態における動作と同じである。
プ）回路２７と補助変圧器５とにより、同定用信号生成
回路２１ａが出力する同定信号の指令値（ｉC *）の位相
を電力系統２の基本波の位相に同期させているが、電力
系統２の周波数特性を同定させるための動作は、図１に
示した実施の形態における動作と同じである。
【００６９】図１１は、この発明の電力系統の周波数特
性同定方法の第５の実施の形態を示す構成図であり、図
４に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
性同定方法の第５の実施の形態を示す構成図であり、図
４に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
【００７０】すなわち図１１に示した構成が図４に示し
た構成と異なる点は、同定用信号生成回路２１に代えて
同定用信号生成回路２１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
た構成と異なる点は、同定用信号生成回路２１に代えて
同定用信号生成回路２１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
【００７１】このＰＬＬ回路２７と補助変圧器５とによ
り、同定用信号生成回路２１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図４に示した実施の形態における動作と同
じである。
り、同定用信号生成回路２１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図４に示した実施の形態における動作と同
じである。
【００７２】図１２は、この発明の電力系統の周波数特
性同定方法の第６の実施の形態を示す構成図であり、図
６に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
性同定方法の第６の実施の形態を示す構成図であり、図
６に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
【００７３】すなわち図１２に示した構成が図６に示し
た構成と異なる点は、同定用信号生成回路２１に代えて
同定用信号生成回路２１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
た構成と異なる点は、同定用信号生成回路２１に代えて
同定用信号生成回路２１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
【００７４】このＰＬＬ回路２７と補助変圧器５とによ
り、同定用信号生成回路２１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図６に示した実施の形態における動作と同
じである。
り、同定用信号生成回路２１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図６に示した実施の形態における動作と同
じである。
【００７５】図１３は、この発明の電力系統の周波数特
性同定方法の第７の実施の形態を示す構成図であり、図
１に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
性同定方法の第７の実施の形態を示す構成図であり、図
１に示した構成と同一機能を有するものには同一符号を
付して、その説明を省略する。
【００７６】すなわち図１３に示した構成が図１に示し
た構成と異なる点は、同定信号生成回路２１と同定値演
算回路２２とに代えて、同定信号生成回路３１と同定値
演算回路３２とを備えていることである。
た構成と異なる点は、同定信号生成回路２１と同定値演
算回路２２とに代えて、同定信号生成回路３１と同定値
演算回路３２とを備えていることである。
【００７７】図１４は図１３に示した同定信号生成回路
３１の動作例を説明する波形図であり、この同定信号生
成回路３１では電力系統２の基本波周波数（紙面では、
基本周期と称している）に対応した５次高調波の正弦波
信号を２π／４〔ｒａｄ〕ずつ位相を違えて、例えば、
紙面上段から順に計４回発生させる機能を有し、これを
同定用信号の指令値（ｉC *）としている。
３１の動作例を説明する波形図であり、この同定信号生
成回路３１では電力系統２の基本波周波数（紙面では、
基本周期と称している）に対応した５次高調波の正弦波
信号を２π／４〔ｒａｄ〕ずつ位相を違えて、例えば、
紙面上段から順に計４回発生させる機能を有し、これを
同定用信号の指令値（ｉC *）としている。
【００７８】図１５はこの発明の第６の実施例を示し、
図１３に示した同定値演算回路３２の動作を説明するフ
ローチャートである。
図１３に示した同定値演算回路３２の動作を説明するフ
ローチャートである。
【００７９】すなわち図１５において、先ず、同定信号
生成回路３１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、図１６のタイミングチャートに示す如く、
図１４に示した正弦波信号を順にパターン１〜パターン
４として前記同定用信号を出力し、その都度サンプリン
グを行い、このときのｘ［ｎ］に対して前記数１の式に
示した如く、離散値系でのフーリエ変換したＸ［ｋ］を
求める。
生成回路３１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、図１６のタイミングチャートに示す如く、
図１４に示した正弦波信号を順にパターン１〜パターン
４として前記同定用信号を出力し、その都度サンプリン
グを行い、このときのｘ［ｎ］に対して前記数１の式に
示した如く、離散値系でのフーリエ変換したＸ［ｋ］を
求める。
【００８０】またインバータ１４から注入された前記指
令値（ｉC *）それぞれに対応する電流（ｉC ）を含みつ
つ、負荷３へ流れる電流（ｉSh）のＣＴ４を介し、前記
サンプリング毎の検出値をｙ［ｎ］とし、前記数２の式
に示した如く、このｙ［ｎ］を離散値系でフーリエ変換
したＹ［ｋ］を求める。
令値（ｉC *）それぞれに対応する電流（ｉC ）を含みつ
つ、負荷３へ流れる電流（ｉSh）のＣＴ４を介し、前記
サンプリング毎の検出値をｙ［ｎ］とし、前記数２の式
に示した如く、このｙ［ｎ］を離散値系でフーリエ変換
したＹ［ｋ］を求める。
【００８１】次に、前記サンプリング毎のＸ［ｋ］，Ｙ
［ｋ］のフーリエ変換した値の相互相関関数のフーリエ
変換した値として、前記数３の式に示した如き、Ｘ
*［ｋ］Ｙ［ｋ］を求める。
［ｋ］のフーリエ変換した値の相互相関関数のフーリエ
変換した値として、前記数３の式に示した如き、Ｘ
*［ｋ］Ｙ［ｋ］を求める。
【００８２】さらに、このサンプリング毎のＸ*［ｋ］
Ｙ［ｋ］をベクトル的な加算平均値をＤＦＴ［ｒ
xy［ｍ］］とし、これに含まれる先述の５次高調波の振
幅を、同定信号生成回路３１が出力する５次高調波の正
弦波信号の振幅の２乗で除算演算したときの商は、電力
系統２の５次高調波のゲイン−周波数特性になり、ここ
で、図１７のベクトル図に示す如く、同定用信号の応答
成分（図示の太字ベクトル）ではない誤差成分（図示の
細字ベクトル）は２π／４ずつ異なった位相（図示のφ
i1〜φi4）で現われるため、前記ベクトル的な加算平均
値ＤＦＴ［ｒxy［ｍ］］では消去される。
Ｙ［ｋ］をベクトル的な加算平均値をＤＦＴ［ｒ
xy［ｍ］］とし、これに含まれる先述の５次高調波の振
幅を、同定信号生成回路３１が出力する５次高調波の正
弦波信号の振幅の２乗で除算演算したときの商は、電力
系統２の５次高調波のゲイン−周波数特性になり、ここ
で、図１７のベクトル図に示す如く、同定用信号の応答
成分（図示の太字ベクトル）ではない誤差成分（図示の
細字ベクトル）は２π／４ずつ異なった位相（図示のφ
i1〜φi4）で現われるため、前記ベクトル的な加算平均
値ＤＦＴ［ｒxy［ｍ］］では消去される。
【００８３】なお、上述の同定信号生成回路３１では５
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的なゲイン−周波数特性を、高調波電
流発生源を有する負荷３に対して、この高調波電流成分
を除去した形で求めることができる。
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的なゲイン−周波数特性を、高調波電
流発生源を有する負荷３に対して、この高調波電流成分
を除去した形で求めることができる。
【００８４】図１８は、この発明の電力系統の周波数特
性同定方法の第８の実施の形態を示す構成図であり、図
１３に示した構成と同一機能を有するものには同一符号
を付して、その説明を省略する。
性同定方法の第８の実施の形態を示す構成図であり、図
１３に示した構成と同一機能を有するものには同一符号
を付して、その説明を省略する。
【００８５】すなわち図１８に示した構成が図１３に示
した構成と異なる点は、同定値演算回路３２に代えて、
同定値演算回路３３又は同定値演算回路３４のいずれか
を備えていることである。
した構成と異なる点は、同定値演算回路３２に代えて、
同定値演算回路３３又は同定値演算回路３４のいずれか
を備えていることである。
【００８６】図１９はこの発明の第７の実施例を示し、
図１８に示した同定値演算回路３３の動作を説明するフ
ローチャートである。
図１８に示した同定値演算回路３３の動作を説明するフ
ローチャートである。
【００８７】すなわち図１９において、先ず、同定信号
生成回路３１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記サンプリング毎の前記数１の式に示し
た如き、このｘ［ｎ］を離散値系でフーリエ変換したＸ
［ｋ］を求める。
生成回路３１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記サンプリング毎の前記数１の式に示し
た如き、このｘ［ｎ］を離散値系でフーリエ変換したＸ
［ｋ］を求める。
【００８８】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
【００８９】次に、前記Ｘ［ｋ］の自己相関関数のフー
リエ変換した値の前記サンプリング毎の｜Ｘ［ｋ］｜2
を求める。
リエ変換した値の前記サンプリング毎の｜Ｘ［ｋ］｜2
を求める。
【００９０】また、前記Ｘ［ｋ］とＹ［ｋ］との相互相
関関数のフーリエ変換した値の前記サンプリンク毎のＤ
ＦＴ［ｒxy［ｍ］］を求める。
関関数のフーリエ変換した値の前記サンプリンク毎のＤ
ＦＴ［ｒxy［ｍ］］を求める。
【００９１】さらに、このサンプリング毎のＸ*［ｋ］
Ｙ［ｋ］をベクトル的な加算平均値をＤＦＴ［ｒ
xy［ｍ］］とし、これに含まれる先述の５次高調波の振
幅を、前記サンプリング毎の前記｜Ｘ［ｋ］｜2 のベク
トル的な加算平均値ＤＦＴ［ｒxX［ｍ］］に含まれる５
次高調波の振幅で除算演算したときの商は、前記５次高
調波における電力系統２のゲイン−周波数特性になる。
Ｙ［ｋ］をベクトル的な加算平均値をＤＦＴ［ｒ
xy［ｍ］］とし、これに含まれる先述の５次高調波の振
幅を、前記サンプリング毎の前記｜Ｘ［ｋ］｜2 のベク
トル的な加算平均値ＤＦＴ［ｒxX［ｍ］］に含まれる５
次高調波の振幅で除算演算したときの商は、前記５次高
調波における電力系統２のゲイン−周波数特性になる。
【００９２】なお、上述の同定信号生成回路３１では５
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的なゲイン−周波数特性を、高調波電
流発生源を有する負荷３に対して、この高調波電流成分
を除去した形で求めることができる。
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的なゲイン−周波数特性を、高調波電
流発生源を有する負荷３に対して、この高調波電流成分
を除去した形で求めることができる。
【００９３】図２０はこの発明の第８の実施例を示し、
図１８に示した同定値演算回路３４の動作を説明するフ
ローチャートである。
図１８に示した同定値演算回路３４の動作を説明するフ
ローチャートである。
【００９４】すなわち図２０において、先ず、同定信号
生成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記サンプリング毎の前記数１の式に示し
た如く、このｘ［ｎ］を離散値系でフーリエ変換したＸ
［ｋ］を求める。
生成回路２１からの同定用信号の指令値（ｉC *）をｘ
［ｎ］とし、前記サンプリング毎の前記数１の式に示し
た如く、このｘ［ｎ］を離散値系でフーリエ変換したＸ
［ｋ］を求める。
【００９５】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
【００９６】次に、前記サンプリング毎に得られた前記
Ｘ［ｋ］とＹ［ｋ］との相互相関関数のフーリエ変換し
た値として、前記数３の式に示した如き、Ｘ*［ｋ］Ｙ
［ｋ］を求める。
Ｘ［ｋ］とＹ［ｋ］との相互相関関数のフーリエ変換し
た値として、前記数３の式に示した如き、Ｘ*［ｋ］Ｙ
［ｋ］を求める。
【００９７】さらに、前記サンプリング毎に得られたＸ
*［ｋ］Ｙ［ｋ］のベクトル的加算平均値をＤＦＴ［ｒ
xy［ｍ］］とし、このＤＦＴ［ｒxy［ｍ］］における位
相が、電力系統２の５次高調波の位相−周波数特性とな
る。
*［ｋ］Ｙ［ｋ］のベクトル的加算平均値をＤＦＴ［ｒ
xy［ｍ］］とし、このＤＦＴ［ｒxy［ｍ］］における位
相が、電力系統２の５次高調波の位相−周波数特性とな
る。
【００９８】なお、上述の同定信号生成回路３１では５
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的な位相−周波数特性を、高調波電流
発生源を有する負荷３に対して、この高調波電流成分を
除去した形で求めることができる。
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的な位相−周波数特性を、高調波電流
発生源を有する負荷３に対して、この高調波電流成分を
除去した形で求めることができる。
【００９９】図２１は、この発明の電力系統の周波数特
性同定方法の第９の実施の形態を示す構成図であり、図
１３に示した構成と同一機能を有するものには同一符号
を付して、その説明を省略する。
性同定方法の第９の実施の形態を示す構成図であり、図
１３に示した構成と同一機能を有するものには同一符号
を付して、その説明を省略する。
【０１００】すなわち図２１に示した構成が図１３に示
した構成と異なる点は、同定値演算回路３２に代えて、
同定値演算回路３５又は同定値演算回路３６のいずれか
を備えていることである。
した構成と異なる点は、同定値演算回路３２に代えて、
同定値演算回路３５又は同定値演算回路３６のいずれか
を備えていることである。
【０１０１】図２２はこの発明の第９の実施例を示し、
図２１に示した同定値演算回路３５の動作を説明するフ
ローチャートである。
図２１に示した同定値演算回路３５の動作を説明するフ
ローチャートである。
【０１０２】すなわち図２２において、先ず、インバー
タ１４から注入された電流（ｉC ）のＣＴ１５を介した
検出値をｘ［ｎ］とし、前記サンプリング毎の前記数１
の式に示した如き、このｘ［ｎ］を離散値系でフーリエ
変換したＸ［ｋ］を求める。
タ１４から注入された電流（ｉC ）のＣＴ１５を介した
検出値をｘ［ｎ］とし、前記サンプリング毎の前記数１
の式に示した如き、このｘ［ｎ］を離散値系でフーリエ
変換したＸ［ｋ］を求める。
【０１０３】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
【０１０４】次に、前記サンプリング毎に得られる前記
Ｘ［ｋ］の自己相関関数のフーリエ変換した値として、
｜Ｘ［ｋ］｜2 を求める。
Ｘ［ｋ］の自己相関関数のフーリエ変換した値として、
｜Ｘ［ｋ］｜2 を求める。
【０１０５】また、前記サンプリング毎に得られる前記
Ｘ［ｋ］とＹ［ｋ］との相互相関関数のフーリエ変換し
た値として、Ｘ*［ｋ］Ｙ［ｋ］を求める。
Ｘ［ｋ］とＹ［ｋ］との相互相関関数のフーリエ変換し
た値として、Ｘ*［ｋ］Ｙ［ｋ］を求める。
【０１０６】このＸ*［ｋ］Ｙ［ｋ］のベクトル的な加
算平均値としてのＤＦＴ［ｒxy［ｍ］］に含まれる先述
の５次高調波の振幅を、前記｜Ｘ［ｋ］｜2 のベクトル
的な加算平均値としてのＤＦＴ［ｒxx［ｍ］］含まれる
５次高調波の振幅で除算演算したときの商は、５次高調
波における電力系統２のゲイン−周波数特性になる。
算平均値としてのＤＦＴ［ｒxy［ｍ］］に含まれる先述
の５次高調波の振幅を、前記｜Ｘ［ｋ］｜2 のベクトル
的な加算平均値としてのＤＦＴ［ｒxx［ｍ］］含まれる
５次高調波の振幅で除算演算したときの商は、５次高調
波における電力系統２のゲイン−周波数特性になる。
【０１０７】なお、上述の同定信号生成回路３１では５
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的なゲイン−周波数特性を、高調波電
流発生源を有する負荷３に対して、この高調波電流成分
を除去した形で求めることができる。
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的なゲイン−周波数特性を、高調波電
流発生源を有する負荷３に対して、この高調波電流成分
を除去した形で求めることができる。
【０１０８】図２３はこの発明の第１０の実施例を示
し、図２１に示した同定値演算回路３６の動作を説明す
るフローチャートである。
し、図２１に示した同定値演算回路３６の動作を説明す
るフローチャートである。
【０１０９】すなわち図２３において、先ず、インバー
タ１４から注入された電流（ｉC ）のＣＴ１５を介した
検出値をｘ［ｎ］とし、前記サンプリング毎の前記数１
の式に示した如き、このｘ［ｎ］を離散値系でフーリエ
変換したＸ［ｋ］を求める。
タ１４から注入された電流（ｉC ）のＣＴ１５を介した
検出値をｘ［ｎ］とし、前記サンプリング毎の前記数１
の式に示した如き、このｘ［ｎ］を離散値系でフーリエ
変換したＸ［ｋ］を求める。
【０１１０】またインバータ１４から注入された前記指
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
令値（ｉC *）に対応する電流（ｉC）を含みつつ、負荷
３へ流れる電流（ｉSh）のＣＴ４を介した検出値をｙ
［ｎ］とし、前記サンプリング毎の前記数２の式に示し
た如き、このｙ［ｎ］を離散値系でフーリエ変換したＹ
［ｋ］を求める。
【０１１１】次に、前記サンプリング毎の前記Ｘ［ｋ］
とＹ［ｋ］との相互相関関数のフーリエ変換した値とし
て、Ｘ*［ｋ］Ｙ［ｋ］を求める。
とＹ［ｋ］との相互相関関数のフーリエ変換した値とし
て、Ｘ*［ｋ］Ｙ［ｋ］を求める。
【０１１２】このＸ*［ｋ］Ｙ［ｋ］のベクトル的な加
算平均値としてのＤＦＴ［ｒxy［ｍ］］における５次高
調波の位相が、この電力系統２の５次高調波の位相−周
波数特性となる。
算平均値としてのＤＦＴ［ｒxy［ｍ］］における５次高
調波の位相が、この電力系統２の５次高調波の位相−周
波数特性となる。
【０１１３】なお、上述の同定信号生成回路３１では５
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的な位相−周波数特性を、高調波電流
発生源を有する負荷３に対して、この高調波電流成分を
除去した形で求めることができる。
次高調波を発生する例で示したが、電力系統２の基本波
に比例した発振周波数の正弦波信号とすることにより、
電力系統２の全体的な位相−周波数特性を、高調波電流
発生源を有する負荷３に対して、この高調波電流成分を
除去した形で求めることができる。
【０１１４】図２４は、この発明の電力系統の周波数特
性同定方法の第１０の実施の形態を示す構成図であり、
図１３に示した構成と同一機能を有するものには同一符
号を付して、その説明を省略する。
性同定方法の第１０の実施の形態を示す構成図であり、
図１３に示した構成と同一機能を有するものには同一符
号を付して、その説明を省略する。
【０１１５】すなわち図２４に示した構成が図１３示し
た構成と異なる点は、同定用信号生成回路３１に代えて
同定用信号生成回路３１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
た構成と異なる点は、同定用信号生成回路３１に代えて
同定用信号生成回路３１ａを備え、さらに、補助変圧器
５とＰＬＬ回路２７とが付加されている。
【０１１６】このＰＬＬ（フェーズ・ロックド・ルー
プ）回路２７と補助変圧器５とにより、同定用信号生成
回路３１ａが出力する同定信号の指令値（ｉC *）の位相
を電力系統２の基本波の位相に同期させているが、電力
系統２の周波数特性を同定させるための動作は、図１３
に示した実施の形態における動作と同じである。
プ）回路２７と補助変圧器５とにより、同定用信号生成
回路３１ａが出力する同定信号の指令値（ｉC *）の位相
を電力系統２の基本波の位相に同期させているが、電力
系統２の周波数特性を同定させるための動作は、図１３
に示した実施の形態における動作と同じである。
【０１１７】図２５は、この発明の電力系統の周波数特
性同定方法の第１１の実施の形態を示す構成図であり、
図１８に示した構成と同一機能を有するものには同一符
号を付して、その説明を省略する。
性同定方法の第１１の実施の形態を示す構成図であり、
図１８に示した構成と同一機能を有するものには同一符
号を付して、その説明を省略する。
【０１１８】すなわち図２６に示した構成が図１８に示
した構成と異なる点は、同定用信号生成回路３１に代え
て同定用信号生成回路３１ａを備え、さらに、補助変圧
器５とＰＬＬ回路２７とが付加されている。
した構成と異なる点は、同定用信号生成回路３１に代え
て同定用信号生成回路３１ａを備え、さらに、補助変圧
器５とＰＬＬ回路２７とが付加されている。
【０１１９】このＰＬＬ回路２７と補助変圧器５とによ
り、同定用信号生成回路３１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図１８に示した実施の形態における動作と
同じである。
り、同定用信号生成回路３１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図１８に示した実施の形態における動作と
同じである。
【０１２０】図２６は、この発明の電力系統の周波数特
性同定方法の第１２の実施の形態を示す構成図であり、
図２１に示した構成と同一機能を有するものには同一符
号を付して、その説明を省略する。
性同定方法の第１２の実施の形態を示す構成図であり、
図２１に示した構成と同一機能を有するものには同一符
号を付して、その説明を省略する。
【０１２１】すなわち図２６に示した構成が図２１に示
した構成と異なる点は、同定用信号生成回路３１に代え
て同定用信号生成回路３１ａを備え、さらに、補助変圧
器５とＰＬＬ回路２７とが付加されている。
した構成と異なる点は、同定用信号生成回路３１に代え
て同定用信号生成回路３１ａを備え、さらに、補助変圧
器５とＰＬＬ回路２７とが付加されている。
【０１２２】このＰＬＬ回路２７と補助変圧器５とによ
り、同定用信号生成回路３１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図２１に示した実施の形態における動作と
同じである。
り、同定用信号生成回路３１ａが出力する同定信号の指
令値（ｉC *）の位相を電力系統２の基本波の位相に同期
させているが、電力系統２の周波数特性を同定させるた
めの動作は、図２１に示した実施の形態における動作と
同じである。
【０１２３】なお、上述の第１〜第１２の実施の形態で
は、アクティブフィルタを流用し、それぞれの電流に基
づく電力系統の周波数特性同定方法について述べたが、
このアクティブフィルタを用いて、それぞれの電圧に基
づく電力系統の周波数特性同定方法を行うこともでき、
また、同定信号を電力系統に注入する手段として、電力
系統と連系運転する他のインバータなどを流用すること
もできる。
は、アクティブフィルタを流用し、それぞれの電流に基
づく電力系統の周波数特性同定方法について述べたが、
このアクティブフィルタを用いて、それぞれの電圧に基
づく電力系統の周波数特性同定方法を行うこともでき、
また、同定信号を電力系統に注入する手段として、電力
系統と連系運転する他のインバータなどを流用すること
もできる。
【０１２４】
【発明の効果】この発明によれば、電力系統の周波数特
性を同定するために、相関関数のフーリエ変換を入力さ
れるｘ［ｎ］，ｙ［ｎ］のフーリエ変換から直接演算す
ることにより、計算量を減らすことができる。また、こ
のフーリエ変換に高速フーリエ変換のアルゴリズムを用
いることにより、さらに計算量を減らすことができる。
性を同定するために、相関関数のフーリエ変換を入力さ
れるｘ［ｎ］，ｙ［ｎ］のフーリエ変換から直接演算す
ることにより、計算量を減らすことができる。また、こ
のフーリエ変換に高速フーリエ変換のアルゴリズムを用
いることにより、さらに計算量を減らすことができる。
【図１】この発明の第１の実施の形態を示す構成図
【図２】図１の動作を説明する波形図
【図３】この発明の第１の実施例を示すフローチャート
【図４】この発明の第２の実施の形態を示す構成図
【図５】この発明の第２の実施例を示すフローチャート
【図６】この発明の第３の実施例を示すフローチャート
【図７】この発明の第３の実施の形態を示す構成図
【図８】この発明の第４の実施例を示すフローチャート
【図９】この発明の第５の実施例を示すフローチャート
【図１０】この発明の第４の実施の形態を示す構成図
【図１１】この発明の第５の実施の形態を示す構成図
【図１２】この発明の第６の実施の形態を示す構成図
【図１３】この発明の第７の実施の形態を示す構成図
【図１４】図１３の動作を説明する波形図
【図１５】この発明の第６の実施例を示すフローチャー
ト
ト
【図１６】図１５の動作を説明するタイミングチャート
【図１７】図１５の動作を説明するベクトル図
【図１８】この発明の第８の実施の形態を示す構成図
【図１９】この発明の第７の実施例を示すフローチャー
ト
ト
【図２０】この発明の第８の実施例を示すフローチャー
ト
ト
【図２１】この発明の第９の実施の形態を示す構成図
【図２２】この発明の第９の実施例を示すフローチャー
ト
ト
【図２３】この発明の第１０の実施例を示すフローチャ
ート
ート
【図２４】この発明の第１０の実施の形態を示す構成図
【図２５】この発明の第１１の実施の形態を示す構成図
【図２６】この発明の第１２の実施の形態を示す構成図
１…電力系統、１ａ…変圧器、２…電力系統、３…負
荷、４…ＣＴ、５…補助変圧器、１１…補正値生成回
路、１２…切替スイッチ、１３…電流調節器、１４…イ
ンバータ、１５…ＣＴ、２１，２１ａ…同定信号生成回
路、２２〜２６…同定値演算回路、２７…ＰＬＬ回路、
３１，３１ａ…同定信号生成回路、３２〜３６…同定値
演算回路。
荷、４…ＣＴ、５…補助変圧器、１１…補正値生成回
路、１２…切替スイッチ、１３…電流調節器、１４…イ
ンバータ、１５…ＣＴ、２１，２１ａ…同定信号生成回
路、２２〜２６…同定値演算回路、２７…ＰＬＬ回路、
３１，３１ａ…同定信号生成回路、３２〜３６…同定値
演算回路。
───────────────────────────────────────────────────── フロントページの続き Ｆターム(参考） 2G028 AA05 BE06 BF01 CG20 DH05 GL06 MS01 2G036 AA28 BA01 CA12 5G066 EA03 EA10
Claims (10)
- 【請求項１】 発振周波数が単一の正弦波信号からなる
同定用信号又は発振周波数が異なる複数の正弦波信号を
合成してなる同定用信号を生成し、 この同定用信号を電力系統に注入したときの該電力系統
の応答を求め、 求めた応答及び前記同定用信号を個別にフーリエ変換
し、さらに、これらのフーリエ変換した値の相互相関関
数に対するフーリエ変換を行い、 この相互相関関数のフーリエ変換した値から導出される
前記発振周波数毎の振幅を、この発振周波数に対応する
前記正弦波信号の振幅の２乗で除算演算した値を前記周
波数特性の同定値としたことを特徴とする電力系統の周
波数特性同定方法。 - 【請求項２】 発振周波数が単一の正弦波信号からなる
同定用信号又は発振周波数が異なる複数の正弦波信号を
合成してなる同定用信号を生成し、 この同定用信号を電力系統に注入したときの該電力系統
の応答を求め、 求めた応答及び前記同定用信号を個別にフーリエ変換
し、さらに、これらのフーリエ変換した値の相互相関関
数に対するフーリエ変換を行い、 前記同定用信号をフーリエ変換した値の自己相関関数に
対するフーリエ変換を行い、 前記相互相関関数のフーリエ変換した値から導出される
前記発振周波数毎の振幅を前記自己相関関数のフーリエ
変換した値から導出される該発振周波数の振幅で除算演
算した値を、前記周波数特性の同定値としたことを特徴
とする電力系統の周波数特性同定方法。 - 【請求項３】 発振周波数が単一の正弦波信号からなる
同定用信号又は発振周波数が異なる複数の正弦波信号を
合成してなる同定用信号を生成し、 この同定用信号を構成するそれぞれの正弦波信号の位相
を前回より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位
させつつ、電力系統に注入したときの該電力系統の応答
を求め、 前記ｎ個の求めた応答及び同定用信号を個別にフーリエ
変換し、さらに、これらのフーリエ変換した値の相互相
関関数に対するフーリエ変換のベクトル的加算平均値を
求め、 このベクトル的加算平均値から導出される前記発振周波
数毎の振幅をこの発振周波数に対応する前記正弦波信号
の振幅の２乗で除算演算した値を、前記周波数特性の同
定値としたことを特徴とする電力系統の周波数特性同定
方法。 - 【請求項４】 発振周波数が単一の正弦波信号からなる
同定用信号又は発振周波数が異なる複数の正弦波信号を
合成してなる同定用信号を生成し、 この同定用信号を構成するそれぞれの正弦波信号の位相
を前回より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位
させつつ、電力系統に注入したときの該電力系統の応答
を求め、 前記ｎ個の求めた応答及び同定用信号を個別にフーリエ
変換し、さらに、これらのフーリエ変換した値の相互相
関関数に対するフーリエ変換のベクトル的加算平均値を
求め、 前記ｎ個の同定用信号をフーリエ変換した値の自己相関
関数に対するフーリエ変換のベクトル的加算平均値を求
め、 前記相互相関関数のフーリエ変換のベクトル的加算平均
値から導出される前記発振周波数毎の振幅を、前記自己
相関関数のフーリエ変換した値のうちの１つ、若しくわ
前記自己相関関数のフーリエ変換した値のベクトル的加
算平均値から導出される該発振周波数の振幅で除算演算
した値を前記周波数特性の同定値としたことを特徴とす
る電力系統の周波数特性同定方法。 - 【請求項５】 発振周波数が単一の正弦波信号からなる
同定用信号又は発振周波数が異なる複数の正弦波信号を
合成してなる同定用信号を生成し、 この同定用信号を電力系統に注入したときの該電力系統
の応答を求め、 求めた応答及び前記同定用信号を個別にフーリエ変換
し、さらに、これらのフーリエ変換した値の相互相関関
数に対するフーリエ変換を行い、 この相互相関関数をフーリエ変換した値から導出される
前記発振周波数毎の位相量を前記周波数特性の同定値と
したことを特徴とする電力系統の周波数特性同定方法。 - 【請求項６】 発振周波数が単一の正弦波信号からなる
同定用信号又は発振周波数が異なる複数の正弦波信号を
合成してなる同定用信号を生成し、 この同定用信号を構成するそれぞれの正弦波信号の位相
を前回より２π／ｎ〔ｒａｄ〕，（ｎ＝複数）ずつ変位
させつつ、電力系統に注入したときの該電力系統の応答
を求め、 前記ｎ個の求めた応答及び同定用信号を個別にフーリエ
変換し、さらに、これらのフーリエ変換した値の相互相
関関数に対するフーリエ変換のベクトル的加算平均値を
求め、 このベクトル的加算平均値から導出される前記発振周波
数毎の位相量を前記周波数特性の同定値としたことを特
徴とする電力系統の周波数特性同定方法。 - 【請求項７】 請求項１乃至６のいずれかに記載の電力
系統の周波数特性同定方法において、 前記同定用信号は、電力系統と連系運転する半導体電力
変換装置に付加的に指令される同定用信号指令値により
該半導体電力変換装置で生成することを特徴とする電力
系統の周波数特性同定方法。 - 【請求項８】 請求項７に記載の電力系統の周波数特性
同定方法において、 前記同定用信号を前記同定用信号指令値とし、 前記応答を上位の電力系統と前記電力系統との間の電流
又は電圧として前記周波数特性の同定値を演算すること
を特徴とする電力系統の周波数特性同定方法。 - 【請求項９】 請求項７に記載の電力系統の周波数特性
同定方法において、 前記同定用信号を前記半導体電力変換装置の出力電流又
は出力電圧とし、 前記応答を上位の電力系統と前記電力系統との間の電流
又は電圧として前記周波数特性の同定値を演算すること
を特徴とする電力系統の周波数特性同定方法。 - 【請求項１０】 請求項１乃至９のいずれかに記載の電
力系統の周波数特性同定方法において、 前記同定用信号を構成する正弦波信号の発振周波数を、
前記電力系統の基本波周波数に比例した値としたことを
特徴とする電力系統の周波数特性同定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000238795A JP2002058163A (ja) | 2000-08-07 | 2000-08-07 | 電力系統の周波数特性同定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000238795A JP2002058163A (ja) | 2000-08-07 | 2000-08-07 | 電力系統の周波数特性同定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2002058163A true JP2002058163A (ja) | 2002-02-22 |
Family
ID=18730442
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000238795A Pending JP2002058163A (ja) | 2000-08-07 | 2000-08-07 | 電力系統の周波数特性同定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2002058163A (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007119855A1 (ja) * | 2006-04-13 | 2007-10-25 | Hitachi, Ltd. | 電力変換装置およびその制御方法 |
| US7848897B2 (en) | 2008-01-30 | 2010-12-07 | Southern Company Services, Inc. | Dynamic real-time power system monitoring |
| CN112067927A (zh) * | 2020-09-09 | 2020-12-11 | 国网冀北电力有限公司电力科学研究院 | 中高频振荡检测方法及装置 |
| JP2021022861A (ja) * | 2019-07-29 | 2021-02-18 | アルパイン株式会社 | 適応同定システム、適応同定装置、及び適応同定方法 |
-
2000
- 2000-08-07 JP JP2000238795A patent/JP2002058163A/ja active Pending
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007119855A1 (ja) * | 2006-04-13 | 2007-10-25 | Hitachi, Ltd. | 電力変換装置およびその制御方法 |
| US7848897B2 (en) | 2008-01-30 | 2010-12-07 | Southern Company Services, Inc. | Dynamic real-time power system monitoring |
| JP2021022861A (ja) * | 2019-07-29 | 2021-02-18 | アルパイン株式会社 | 適応同定システム、適応同定装置、及び適応同定方法 |
| JP7252086B2 (ja) | 2019-07-29 | 2023-04-04 | アルパイン株式会社 | 適応同定システム、適応同定装置、及び適応同定方法 |
| CN112067927A (zh) * | 2020-09-09 | 2020-12-11 | 国网冀北电力有限公司电力科学研究院 | 中高频振荡检测方法及装置 |
| CN112067927B (zh) * | 2020-09-09 | 2023-08-25 | 国网冀北电力有限公司电力科学研究院 | 中高频振荡检测方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Freijedo et al. | A generic open-loop algorithm for three-phase grid voltage/current synchronization with particular reference to phase, frequency, and amplitude estimation | |
| Lee et al. | AC voltage and current sensorless control of three-phase PWM rectifiers | |
| US7177165B2 (en) | System and method for unbalanced independent AC phase voltage control of a 3-phase, 4-wire output DC/AC inverter | |
| Cirrincione et al. | Sensorless control of induction machines by a new neural algorithm: The TLS EXIN neuron | |
| US20120206945A1 (en) | Controller for a power converter and method of operating the same | |
| EP1741178B1 (en) | Power converter apparatus and methods using a phase reference derived from a dc bus voltage | |
| CN114665512B (zh) | 一种并网逆变器无交流电压传感器预测控制装置及方法 | |
| CN106546851A (zh) | 一种mmc换流阀运行试验电路的稳定控制方法及装置 | |
| Patil et al. | Modified SOGI based shunt active power filter to tackle various grid voltage abnormalities | |
| Guo et al. | A virtual-flux state observer-based inductance identification method for model predictive control of grid-tied inverters with a finite phase angle set-based PLL | |
| JP2002058163A (ja) | 電力系統の周波数特性同定方法 | |
| Hu et al. | Predictive current control of grid-connected voltage source converters during network unbalance | |
| JP3528475B2 (ja) | 電力用アクティブフィルタ | |
| Ebadpour et al. | EKF digital twinning of induction motor drives for the metaverse | |
| JP4055992B2 (ja) | インバータの電流検出装置 | |
| Suresh | Discrete switching strategy for direct torque control of induction motor using adaptive pulse width modulator | |
| JP3676133B2 (ja) | 電力系統の周波数特性同定方法 | |
| JP2644750B2 (ja) | 電圧形インバータの制御方法 | |
| Hamasaki et al. | A novel method for active filter current regulation using deadbeat control | |
| JP4298040B2 (ja) | 電源装置及びこれを用いた磁気共鳴イメージング装置 | |
| Chatterjee et al. | Current Zero-Crossing Shift During Compensation of Dead-Time Distortion in PWM VSI and Its Corrective Action in Current Sensor-Less Compensation | |
| JP2774246B2 (ja) | 電流形コンバータの制御装置 | |
| Hu et al. | Improved predictive current control of grid-connected DC-AC converters under unbalanced grid voltage conditions | |
| JPH11252796A (ja) | 電力系統高調波リアルタイムシミュレータ | |
| KR100635717B1 (ko) | 전력변환설비용 위상 동기 루프 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060215 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20071127 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071213 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080515 |