(57)【特許請求の範囲】
【請求項１】
電力系統の線路インピーダンス及び前記電力系統の力率改善用コンデンサを線形回路素子ＲＬＣよりなる解析モデルで模擬して電力系統の状態微分方程式を導き、
前記状態微分方程式の係数行列の固有値で定まる時定数及び固有周波数を持つＲＬＣ回路による等価回路を導き、
前記等価回路の力率改善用コンデンサ電流の高次高調波電流源の電流に対する電流増幅率に基づいて電力系統における高次高調波の共振周波数特性を推定することを特徴とする高次高調波共振周波数特性推定方法。
【請求項２】
表皮効果による前記線路インピーダンスの周波数特性を線形素子ＲＬのＲＬ並列回路よりなる解析モデルで追加模擬して前記状態微分方程式を導くことを特徴とする請求項１記載の高次高調波共振周波数特性推定方法。
【請求項３】
前記解析モデルは、前記ＲＬ並列回路を直列多段接続し、模擬可能な周波数領域を拡大したことを特徴とする請求項２記載の高次高調波共振周波数特性推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統における高次高調波の共振周波数特性を推定する高次高調波共振周波数特性推定方法に関する。
【背景技術】
【０００２】
近年、パワーエレクトロニクスによる制御には、高速スイッチングによるＰＷＭ方式を採用する機器が多い。例えば、風力発電や燃料電池などにおいては直流電力を発電し、パワーエレクトロニクス機器によるＰＷＭ方式で直流電力を交流電力に変換して電力系統に供給する。
【０００３】
パワーエレクトロニクス機器は、ＰＷＭ方式による交直変換の際に周波数ｆ０の高次高調波電流を発生し電力系統に流出する。一方、電力系統には、力率改善用コンデンサが設置されており、パワーエレクトロニクス機器が発生する周波数ｆ０と近隣の電力系統の力率改善用コンデンサによる固有周波数ｆｎとが近接している場合には、共振によりパワーエレクトロニクス機器が発生する周波数ｆ０の電流が増幅され、力率改善用コンデンサの過電流を誘起し電磁騒音や過熱など電力品質問題が生ずる。
【０００４】
このため、パワーエレクトロニクス機器の系統連系時には、電力品質問題が生ずるか否かを事前に把握することが重要となってきている。なお、ここでの高次高調波は次数が１３次を越え、上限を１０ｋＨｚとする周波数領域としている。
【０００５】
高次高調波共振を解析的に検討する手法としては、高次高調波電流を発生するパワーエレクトロニクス機器の連系点に等価な電流源を置き、周波数応答法で力率改善用コンデンサの電流の増幅特性を得る方法がある。また、高調波発生源の特定方法として、配電系統を数値的に模擬した模擬回路を設定し、この模擬回路の固有値とその固有ベクトルを算出し、固有値の複素数部からその固有モードの共振周波数を、それぞれの固有ベクトルの各要素の絶対値からその固有モードにおける各ノードの感度を算出し、一方、配電系統に並列に接続されている力率改善用コンデンサの内の一つの充電波形に基づいて周波数分析を行って実際の発生モードでの共振周波数とその強度を求め、これらを比較して転流振動を発生させる整流器負荷を特定し、高調波発生源と見なすようにしたものがある（例えば、特許文献１参照）。
【特許文献１】特開平４−６７７２７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかし、特許文献１のものでは、模擬回路の固有値及びその固有ベクトルから求めた共振周波数及び感度と、配電系統で実際に発生している発生モードでの共振周波数及びその強度を求めて、転流振動を発生させる整流器負荷を特定して高調波発生源を特定するものであり、新たにパワーエレクトロニクス機器を系統に連系しようとした場合の高次高調波の共振周波数特性を把握することはできない。従って、新たにパワーエレクトロニクス機器を系統に連系しようとする際に、その連系に問題がないか否かの判定に用いることができない。
【０００７】
一方、新たに系統連系しようとするパワーエレクトロニクス機器の連系点に等価な電流源を置き、周波数応答法で力率改善用コンデンサの電流の増幅特性を得る方法では、多数の周波数に対する計算例が必要となり煩雑である。
【０００８】
また、線路インピーダンスＲＬは表皮効果などで周波数により変化する。高次高調波の領域では、特に抵抗分Ｒが急激に増大する。この抵抗分Ｒの増大を無視すると精度の良い周波数応答を得ることができず実用性に乏しくなるので、これを考慮した解析手法が望まれている。
【０００９】
本発明の目的は、新たにパワーエレクトロニクス機器を系統に連系しようとした場合に高次高調波の共振周波数特性を精度よく把握できる高次高調波共振周波数特性推定方法を提供することである。
【課題を解決するための手段】
【００１０】
請求項１の発明に係わる高次高調波共振周波数特性推定方法は、電力系統の線路インピーダンス及び前記電力系統の力率改善用コンデンサを線形回路素子ＲＬＣよりなる解析モデルで模擬して電力系統の状態微分方程式を導き、前記状態微分方程式の係数行列の固有値で定まる時定数及び固有周波数を持つＲＬＣ回路による等価回路を導き、前記等価回路の力率改善用コンデンサ電流の高次高調波電流源の電流に対する電流増幅率に基づいて電力系統における高次高調波の共振周波数特性を推定することを特徴とする。
【００１１】
請求項２の発明に係わる高次高調波共振周波数特性推定方法は、請求項１の発明において、表皮効果による前記線路インピーダンスの周波数特性を線形素子ＲＬのＲＬ並列回路よりなる解析モデルで追加模擬して前記状態微分方程式を導くことを特徴とする。
【００１２】
請求項３の発明に係わる高次高調波共振周波数特性推定方法は、請求項２の発明において、前記解析モデルは、前記ＲＬ並列回路を直列多段接続し、模擬可能な周波数領域を拡大したことを特徴とする。
【発明の効果】
【００１３】
本発明によれば、電力系統の状態微分方程式の係数行列の固有値で定まる時定数及び固有周波数を持つＲＬＣ回路による等価回路を導き、その等価回路の力率改善用コンデンサ電流の高次高調波電流源の電流に対する電流増幅率に基づいて電力系統における高次高調波の共振周波数特性を推定するので、高次高調波電流源であるパワーエレクトロニクス機器を新たに系統に連系しようとする際に、その連系に問題がないか否かの判定を容易に行うことができる。
【００１４】
また、表皮効果による線路インピーダンスの周波数特性を線形素子ＲＬのＲＬ並列回路からなる解析モデルで追加模擬して電力系統の状態微分方程式を導くので、表皮効果を考慮に入れた線路インピーダンスの高次高調波の共振周波数特性を精度よく推定できる。さらに、ＲＬ並列回路を直列多段接続した解析モデルを用いることにより、模擬可能な周波数領域を拡大できる。
【発明を実施するための最良の形態】
【００１５】
まず、本発明に至るまでの検討事項について説明する。電力系統の高次高調波共振を検討するあたっては、電力系統の解析モデルを作成し、その解析モデルを用いて固有値及びその固有ベクトルを求めて共振周波数及びその感度を得ることはできる。しかし、解析モデルの固有値及びその固有ベクトルでは、振動が発生する周波数及び振動の強弱は分かるが、新たにパワーエレクトロニクス機器を系統に連系した場合に、そのパワーエレクトロニクス機器が発生する高次高調波電流が力率改善用コンデンサによりどの程度電流増幅されるかどうかを把握することはできない。
【００１６】
そこで、力率改善用コンデンサ電流の高次高調波電流源の電流に対する電流増幅率を求めるあたり、本発明では、解析モデルの固有値は電力系統では一般に複素数であり、その固有値の実数部は時定数、虚数部は固有周波数となることに着目し、固有値で定まる時定数及び固有周波数を持つＲＬＣ回路による等価回路を導くこととした。そして、導き出した等価回路に高次高調波電流源を接続して高次高調波電流を変化させ、等価回路での力率改善用コンデンサ電流の高次高調波電流源の電流に対する電流増幅率を求め、その電流増幅率に基づいて電力系統における高次高調波の共振周波数特性を推定するようにした。
【００１７】
図１は本発明の実施の形態に係わる高次高調波共振周波数特性推定方法の内容を示すフローチャートである。まず、電力系統の解析モデルの固有値を求めるために電力系統の状態微分方程式を求める（Ｓ１）。以下、電力系統の一例として高圧配電系統について説明する。また、電力系統の解析モデルを簡潔化するため、次のような仮定を設定する。
【００１８】
（１）線路インピーダンスは三相平衡であるとする。
【００１９】
（２）線路キャパシタンスは力率改善用コンデンサＣに比べ十分小であるので無視する。
【００２０】
これにより、解析モデルは、図２に示すように、作用インダクタンスＬ、抵抗Ｒ、力率改善用コンデンサＣなどの定数からなる単相回路となる。また、解析時の変数は、上位系や線路インピーダンスＲＬのブランチ１１ａ〜１１ｎのブランチ電流Ｉ（Ｉ１、Ｉ２…Ｉｎ）と、力率改善用コンデンサＣ（Ｃ１、Ｃ２…Ｃｎ）の接続点におけるブランチ１２ａ〜１２ｎのノード電圧Ｖ（Ｖ１、Ｖ２…Ｖｎ）とを採用する。また、電源は高次高調波の等価電流源Ｊのみとする。以上の条件下で電力系統の状態微分方程式を求める。
【００２１】
図２に示す力率改善用コンデンサＣ（Ｃ１、Ｃ２…Ｃｎ）の接続点であるノード１２ａ〜１２ｎでの電流バランス（ノード電流平衡式）を求め、ブランチ１１ａ〜１１ｎでの電圧バランス（ブランチ電圧平衡式）を求める。ノード電流平衡式は下記（１）式で示され、ブランチ電圧平衡式は下記（２）式で示される。
【数１】
0004845133.tif 000002
【数２】
0004845133.tif 000003
【００２２】
（１）、（２）式により、電力系統の状態微分方程式は（３）式で示される。
【数３】
0004845133.tif 000004
【００２３】
次に、（３）式で示される状態微分方程式の係数行列の固有値λを求める（Ｓ２）。この（３）式の係数行列の固有値λは、電力系統では一般に複素数となり、よく知られるように、時定数Ｔ０、固有周波数ｆ０としたとき、実数部及び虚数部は（４）式で示すような関係となる。
【数４】
0004845133.tif 000005
【００２４】
この固有値λは一般に複数個求まるので、固有値λごとに時定数Ｔ０及び固有振動数ｆ０を求める（Ｓ３）。
【００２５】
次に、時定数Ｔ０及び固有周波数ｆ０を持つＲＬＣ回路による等価回路を求める（Ｓ４）。この等価回路は固有値λの個数だけ求め、各固有値λの固有モードでの現象を個々のモード別に考察することで、その物理的評価が容易かつ明確となる。
【００２６】
図３は、時定数Ｔ０及び固有周波数ｆ０を持つＲＬＣ回路による等価回路の回路図である。等価回路は、線路インピーダンスＲ０、Ｌ０、力率改善用コンデンサＣのＲＬＣ回路であり、力率改善用コンデンサＣに高次高調波の電流源Ｊλが並列配置されている。この等価回路の抵抗Ｒ０、作用インダクタンスＬ０、力率改善用コンデンサＣと、時定数Ｔ０、固有周波数ω０（ω０＝２πｆ０）との関係は、次の（５）式のようになる。
【数５】
0004845133.tif 000006
【００２７】
そして、得られた等価回路の力率改善用コンデンサＣの電流ｉｃを求め、この電流ｉｃの電流源Ｊλに対する増幅率Ｑ（ω）を求める（Ｓ５）。この増幅率Ｑ（ω）の特性により、共振周波数特性を推定する（Ｓ６）。
【００２８】
すなわち、図３に示す等価回路の方程式は下記の（６）式で示される。
【数６】
0004845133.tif 000007
【００２９】
この（６）式より増幅率Ｑ（ω）は（７）式で示される。
【数７】
0004845133.tif 000008
【００３０】
また、（７）式内のＬ０Ｃ、ＣＲ０は（５）式から（８）式で示される。
【数８】
0004845133.tif 000009
【００３１】
従って、（７）式のＬ０Ｃ、ＣＲ０に（８）式を代入すると、各固有モードの増幅率Ｑ（ω）は、角周波数ωと、時定数Ｔ０と固有角周波数ω０との簡潔な式で示されるので、各固有モードの増幅率Ｑ（ω）は計算で容易に求まる。
【００３２】
このように、電力系統の状態微分方程式の係数行列の固有値λで定まる時定数Ｔ０及び固有周波数ｆ０を持つＲＬＣ回路による等価回路を導き、その等価回路の力率改善用コンデンサ電流ｉｃの高次高調波電流源Ｊλの電流に対する増幅率Ｑ（ω）に基づいて電力系統における高次高調波の共振周波数特性を推定するので、高次高調波電流源であるパワーエレクトロニクス機器を新たに系統に連系しようとする際に、その連系に問題がないか否かの判定を容易に行うことができる。
【００３３】
次に、線路インピーダンスＲＬの表皮効果の考慮について説明する。線路インピーダンスＲＬは図４に示すように周波数特性を持っており、特に抵抗Ｒは周波数ｆの上昇により急激に増大する一方、作用インダクタンスＬは僅かながら低下する。共振による電流増幅は（Ｌ／Ｒ）に比例するため、これを無視しては実用的な精度確保ができない。そこで、次に抵抗Ｒ及び作用インダクタンスＬの周波数ｆ依存性の模擬方法について検討する。
【００３４】
周波数ｆの依存性を考慮する場合のポイントは、固有モード法の利点を生かせる線形の状態微分方程式となる模擬法が必須のことにある。関数の近似法に多く用いるラグランジェ補間など多項式近似では、方程式は非線形となり固有値などは得られない。
【００３５】
そこで、状態微分方程式を線形に保つ模擬方法として、図５に示す線形素子Ｒｐ、ＬｐのＲＬ並列回路を考える。このＲＬ並列回路の等価抵抗Ｒｅ、等価インダクタンスＬｅは、よく知られるように次の（９）式で与えられる。
【数９】
0004845133.tif 000010
【００３６】
（９）式より周波数ｆとともに等価抵抗Ｒｅは増大し、等価インダクタンスＬｅは減少するので、線路インピーダンスＲＬの対周波数特性は模擬可能となる。これより図６に示すように、模擬回路は、ベース分を抵抗Ｒｓ及び作用インダクタンスＬｓの直列ＲＬ回路で、対周波数特性を等価抵抗Ｒｅ及び等価インダクタンスＬｓのＲＬ並列回路で分担する回路となる。
【００３７】
一方、等価抵抗Ｒｅは図７のように変曲点ｆｍを境にして増大率は低下するので、模擬すべき周波数領域が広い場合には、変曲点ｆｍより高い周波数領域の近似精度は低下する。変曲点ｆｍは（１０）式で示される。
【数１０】
0004845133.tif 000011
【００３８】
また、等価抵抗Ｒｅは（１１）式で示される。
【数１１】
0004845133.tif 000012
【００３９】
変曲点ｆｍを高く設定すれば、（１１）式より等価抵抗Ｒｅは、増大開始の周波数ｆが高くなり低い周波数領域の近似精度が低下する。図８は１つのＲＬ並列回路を用いた抵抗Ｒの対周波数特性を示す特性図である。点線は抵抗Ｒ（Ｒ＝Ｒｓ＋Ｒｅ）の対周波数特性、実線は抵抗Ｒ（Ｒ＝Ｒｓ）の対周波数特性である。点線の抵抗Ｒ（Ｒ＝Ｒｓ＋Ｒｅ）の対周波数特性から分かるように、模擬可能な周波数ｆは３ｋＨｚまで達していない。
【００４０】
これにより、１つのＲＬ並列回路で近似可能な周波数領域は限定されるのが分かる。そこで、変曲点ｆｍがより高い第２のＲＬ並列回路を追加し直列接続することで、より高い周波数領域を模擬する。
【００４１】
図９は対周波数特性を模擬する２つのＲＬ並列回路を直列接続した模擬回路図である。この場合の、１０ｋＨｚまでの周波数領域で、最大誤差は一般に約２％程度であり、十分な実用性を有している。さらに、精度向上や模擬周波数領域を拡大する場合には、ＲＬ並列回路の直列段数を増すことで対応可能である。
【００４２】
次に、ＲＬ並列回路の追加により、周波数特性を模擬した場合の回路方程式が線形であることを、２つのＲＬ並列回路を直列接続した場合を例に取り説明する。ＲＬ並列回路の追加による状態微分方程式への影響はインダクタンスと抵抗の係数行列に現れる。
【００４３】
いま、図９において、ブランチｎの両端ノードをｊ、ｋとし、電流Ｉｎはｊからｋに流れるとする。そして、第１のＲＬ並列回路の抵抗Ｒｐ１を流れる電流をＩｐ１，ｎ、第２のＲＬ並列回路の抵抗Ｒｐ２を流れる電流をＩｐ２，ｎ、とする。この場合、電圧平衡式は（１２）式で示される。
【数１２】
0004845133.tif 000013
【００４４】
これを全ブランチ体で表記すれば、（１３）式で示される。
【数１３】
0004845133.tif 000014
【００４５】
これより、全体の回路方程式は次の（１４）式のように、固有値計算可能な線形の状態微分方程式になり、ＲＬ並列回路が２つ以上に増えた場合でも、一般性を失わない。
【数１４】
0004845133.tif 000015
【００４６】
前述の（３）式の場合と同様に、（１４）式で示される状態微分方程式の係数行列の固有値λを求め、固有値λごとに時定数Ｔ０及び固有振動数ｆ０を求め、時定数Ｔ０及び固有周波数ｆ０を持つＲＬＣ回路による等価回路を求める。そして、得られた等価回路の力率改善用コンデンサＣの電流ｉｃを求め、この電流ｉｃの電流源Ｊλに対する増幅率Ｑ（ω）を求め、この増幅率Ｑ（ω）の特性により共振周波数特性を推定する（Ｓ６）。
【００４７】
このように、線路インピーダンスＲＬの表皮効果による周波数特性を線形素子ＲＬのＲＬ並列回路よりなる解析モデルで追加模擬して電力系統の状態微分方程式を導くので、線路インピーダンスＲＬの表皮効果を考慮に入れた高次高調波の共振周波数特性を推定できる。さらに、線形素子ＲＬのＲＬ並列回路を直列多段接続した解析モデルを用いることにより、模擬可能な周波数領域を拡大できる。
【図面の簡単な説明】
【００４８】
【図１】本発明の実施の形態に係わる高次高調波共振周波数特性推定方法の内容を示すフローチャート。
【図２】本発明の実施の形態における電力系統の解析モデルの一例を示す模擬図。
【図３】本発明の実施の形態における電力系統の解析モデルの時定数及び固有周波数を持つＲＬＣ回路による等価回路の回路図。
【図４】電力系統の線路インピーダンスＲＬの周波数特性図。
【図５】電力系統の線路インピーダンスＲＬの対周波数特性を模擬するＲＬ並列回路の回路図。
【図６】対周波数特性を模擬するＲＬ並列回路を備えた電力系統の線路インピーダンスＲＬの模擬回路。
【図７】対周波数特性を模擬するＲＬ並列回路の等価抵抗の周波数特性図。
【図８】１つのＲＬ並列回路を用いた抵抗Ｒの対周波数特性を示す特性図。
【図９】対周波数特性を模擬する２つのＲＬ並列回路を直列接続した模擬回路図。
【符号の説明】
【００４９】
１１…ブランチ、１２…ノード
【図１】
0004845133.tif 000016
【図２】
0004845133.tif 000017
【図３】
0004845133.tif 000018
【図４】
0004845133.tif 000019
【図５】
0004845133.tif 000020
【図６】
0004845133.tif 000021
【図７】
0004845133.tif 000022
【図８】
0004845133.tif 000023
【図９】
0004845133.tif 000024