(57)【特許請求の範囲】
【請求項１】
分散電源が接続された系統に系統基本波の非整数倍周波数の次数間高調波を注入し、
系統の注入次数の次数間高調波の計測に基づき、系統の前記注入次数の次数間高調波のインピーダンス又はアドミタンスを検出し、
系統の前記注入次数の次数間高調波の検出インピーダンス又は検出アドミタンスの変化から、系統停止時の前記分散電源の単独運転を検出する分散電源の単独運転検出方法であって、
系統に前記注入次数の次数間高調波を単相注入し、
系統各相の前記注入次数の次数間高調波の計測された電圧，電流それぞれを対称成分に分解して前記電圧，前記電流それぞれの正相分，逆相分を検出し、
前記系統基本波の正相分の変化量が逆相分の変化量以下であれば、
前記検出インピーダンス又は前記検出アドミタンスとして、前記注入次数の次数間高調波の電圧，電流それぞれの前記正相分から、系統の前記注入次数の次数間高調波の正相インピーダンス又は正相アドミタンスを検出し、
前記系統基本波の正相分の変化量が逆相分の変化量より大きければ、
前記検出インピーダンス又は前記検出アドミタンスとして、前記注入次数の次数間高調波の電圧，電流それぞれの前記逆相分から、系統の前記注入次数の次数間高調波の逆相インピーダンス又は逆相アドミタンスを検出する
ことを特徴とする分散電源の単独運転検出方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、系統に接続された分散電源の系統停止時の単独運転を検出する分散電源の単独運転検出方法に関する。
【０００２】
【従来の技術】
従来、電力系統の事故停電等による停止時に需要家のこの種分散電源の単独運転を検出してその運転を停止するため、本出願人は、例えば特開平１０−２４８１６８号公報（Ｈ０２Ｊ  ３／３８）等に記載されているように、電力系統に系統基本波の非整数倍周波数の次数間高調波（中間次数調波）の電流を注入し、系統の電圧，電流の計測信号から注入次数の次数間高調波の電圧，電流を検出し、この電圧，電流から、系統の注入次数の次数間高調波のインピーダンス又はアドミタンスを検出し、これらの変化（変動）から、系統停止時の分散電源の単独運転を検出する発明を、既に出願している。
【０００３】
【発明が解決しようとする課題】
前記従来の単独運転検出方法においては、３相系統に適用する場合、前記既出願の明細書，図面等にも記載されているように、注入電源（中間調波電流注入装置）を３相電源とし、系統各相に注入次数の次数間高調波の３相電流を注入する必要がある。
【０００４】
したがって、次数間高調波の注入電源として、高価かつ大型の３相電源を要する問題点がある。
【０００５】
そこで、注入電源を単相電源にしてコストダウン及び小型化を図り、系統の２相間に次数間高調波を単相注入して系統停止時の分散電源の単独運転を検出することが考えられるが、この単相注入による具体的な単独運転の検出手法は発明されておらず、単相注入した場合には系統のバンクトランス等での相順変化による注入相と検出相ののずれ等が問題となる。
【０００６】
さらに、注入次数の次数間高調波の電圧，電流が、系統の負荷変動等による基本波の変動等の影響を受け、誤検出が生じるおそれもある。
【０００７】
そして、単独運転の誤検出を極力防止し、系統正常時に分散電源を系統から誤って切離さないようにすることが望まれる。
【０００８】
本発明は、次数間高調波を単相注入して系統の注入次数の次数間高調波のインピーダンス又はアドミタンスの変化から系統停止時の分散電源の単独運転を検出する際に、系統の相順変化を気にすることなく、しかも、基本波の変動に基づく系統の変化に応じた適切な検出方法により、誤検出を防止して分散電源の単独運転を確実に検出することを課題とする。
【０００９】
【課題を解決するための手段】
前記の課題を解決するために、本発明の分散電源の単独運転検出方法は、分散電源が接続された系統に系統基本波の非整数倍周波数の次数間高調波を注入し、
系統の注入次数の次数間高調波の計測に基づき、系統の前記注入次数の次数間高調波のインピーダンス又はアドミタンスを検出し、
系統の注入次数の次数間高調波の検出インピーダンス又は検出アドミタンスの変化から、系統停止時の分散電源の単独運転を検出する分散電源の単独運転検出方法であって、
系統に注入次数の次数間高調波を単相注入し、
系統各相の注入次数の次数間高調波の計測された電圧，電流それぞれを対称成分に分解して前記電圧，前記電流それぞれの正相分，逆相分を検出し、
系統基本波の正相分の変化量が逆相分の変化量以下であれば、
検出インピーダンス又は検出アドミタンスとして、注入次数の次数間高調波の電圧，電流それぞれの正相分から、系統の注入次数の次数間高調波の正相インピーダンス又は正相アドミタンスを検出し、
系統基本波の正相分の変化量が逆相分の変化量より大きければ、
検出インピーダンス又は検出アドミタンスとして、注入次数の次数間高調波の電圧，電流それぞれの前記逆相分から、系統の前記注入次数の次数間高調波の逆相インピーダンス又は逆相アドミタンスを検出する。
【００１０】
したがって、例えば３相系統であっても次数間高調波は単相注入すればよく、注入電源として３相電源より安価かつ小型の単相電源を用いることができる。
【００１１】
つぎに、単相注入に基づく系統各相の注入次数の次数間高調波の計測された電圧，電流それぞれの正相分又は逆相分が検出される。
【００１２】
ところで、系統の負荷変動等により、系統基本波の正相分の変化量が逆相分の変化量以下になるときは、系統に注入された次数間高調波の正相分が、逆相分より系統基本波の変動の影響を受けにくく、逆に、系統基本波の正相分の変化量が逆相分の変化量より大きくなるとき、換言すれば、系統基本波の逆相分の変化量が正相分の変化量より小さくなるときは、系統に注入された次数間高調波の逆相分が正相分より系統基本波の変動を受けにくいことが判明した。
【００１３】
そして、系統基本波の正相分の変化量が逆相分の変化量以下になる系統状態であれば、系統の注入次数の次数間高調波の検出インピーダンス又は検出アドミタンスとして、正相インピーダンス又は正相アドミタンスが検出され、この正相インピーダンス又は正相アドミタンスの変化から、系統基本波の変動の影響を極力防止して系統停止時の分散電源の単独運転が検出される。
【００１４】
また、系統基本波の正相分の変化量が逆相分の変化量より大きくなる系統状態であれば、系統の注入次数の次数間高調波の検出インピーダンス又は検出アドミタンスとして、逆相インピーダンス又は逆相アドミタンスが検出され、この逆相インピーダンス又は逆相アドミタンスの変化から、系統基本波の変動の影響を極力防止して系統停止時の分散電源の単独運転が検出される。
【００１５】
したがって、コストダウン及び小型化を図った次数間高調波の単相注入により、系統の相順変化等を気にすることなく、しかも、基本波の変動に基づく系統の変化に応じた検出方法で、系統基本波の変動の影響を極力防止して、系統停止時の分散電源の単独運転が確実に検出される。
【００１６】
【発明の実施の形態】
本発明の実施の１形態につき、図１〜図５を参照して説明する。
図１は電力系統の１例である３相配電系統の単線結線図であり、上位系統１に変電所２の１又は複数の変圧器３の１次側が接続され、各変圧器３の２次側から遮断器４を介して１又は複数の下位系統５が分枝状に引出される。
【００１７】
これらの系統５は、分散電源６を有する需要家設備７，分散型電源６が設けられていない一般需要家設備８等の複数の需要家設備が接続される。
【００１８】
そして、需要家設備７は他の需要家設備と同様、系統５に引込線９の遮断器１０を介して負荷母線１１が接続され、この負荷母線１１に各負荷フィーダ１２の変圧器１３を介してそれぞれの負荷が接続される。
【００１９】
また、負荷母線１１に遮断器１４が接続され、この遮断器１４に解列用の開閉器１５を介して分散電源６が接続されるとともに、次数間高調波の単相の電流注入装置１６が接続される。
【００２０】
この電流注入装置１６は次数間高調波の単相の注入電流を出力するインバータ等の安価かつ小型の単相の電源部１７，この電源部１７と負荷母線１１との間に設けられた単相注入用の変圧器１８により形成される。
【００２１】
また、引込線９の遮断器１０より負荷母線１１側に、それぞれ３相の変圧器１９及び変流器２０が設けられ、それらの各相の電圧，電流の計測信号が系統停止検出処理装置２１のサンプル・ホールド回路２２に供給される。
【００２２】
このサンプル・ホールド回路２２はタイミング指令部２３の一定周波数のサンプリング指令のタイミング信号により、電圧，電流の計測信号をサンプル・ホールドし、その出力が後段のＡ／Ｄ変換回路２４によりデジタル信号に変換されて電圧，電流のサンプリングデータとなり、この電圧，電流のサンプリングデータが演算処理部２５に供給される。
【００２３】
そして、演算処理部２５はマイクロコンピュータ等で形成され、そのソフトウェア処理により、両サンプリングデータの公知のフーリエ変換のデジタルフィルタ演算を実行して電流注入装置１６から系統５に注入された注入次数の次数間高調波を検出し、その変化から、遮断器４が開放する系統停止時の分散電源６の単独運転を監視して検出する。
【００２４】
さらに、この単独運転を検出すると、演算処理部２５は開閉器１５に解列の指令を供給し、開閉器１５を開放して分散電源６を系統５から切離す。
【００２５】
ところで、系統正常時は、遮断器４，１０，１４及び開閉器１５がいずれも閉成され、上位系統１の電力が系統５に給電され、系統５は電力供給状態にある。
【００２６】
このとき、分散電源６は系統５に連系して運転され、その出力は自設備７内で消費されるとともに余剰分が引込線９を介して系統５に出力される。
【００２７】
一方、前記の変圧器１９，変流器２０は、図１の受電点Ｐの各相（３相）の電圧，受電点Ｐの各相の引込線９の電流を常時計測する。
【００２８】
そして、演算処理部２５がタイミング指令部２３の一定周波数のタイミング信号に同期して電源部１７に周期的に起動指令を出力し、この指令に基づき、電源部１７が前記タイミング信号に同期した１又は複数の周波数（チャンネル）の次数間高調波の単相の電流を形成し、これらの単相の注入電流が変圧器１８，負荷母線１１，引込線９を介して受電点Ｐから系統５の適当な２相間に注入される。
【００２９】
この注入に基づき、変圧器１９，変流器２０の各相の計測信号には、注入次数の次数間高調波の計測電圧，計測電流が含まれる。
【００３０】
つぎに、演算処理部２５による単独運転の検出処理について説明する。
まず、系統５の停止により遮断器４が開放されると、電圧，電流の計測点である受電点Ｐからみた系統上流側のインピーダンスいわゆる短絡インピーダンスから開放インピーダンスに増大変化し、受電点Ｐからみた系統５の注入次数の次数間高調波についての検出インピーダンス又は検出アドミタンス（以下検出インピーダンス等という）が変化するため、この検出インピーダンス等の変化から、系統５の停止時に、分散電源６の単独運転を検出できる。
【００３１】
そこで、演算処理部２５は、系統５の注入次数の次数間高調波のインピーダンス又はアドミタンスを検出し、その変化から系統停止時の分散電源６の単独運転を監視して検出する。
【００３２】
ところで、系統に単相注入される電流は、基本波を含むｎ次高調波（ｎ：１以上の整数）とｎ＋１次高調波との間の本来系統には存在しない非整数次数Ｍ（ｎ＜Ｍ＜ｎ＋１）の次数間高調波の電流であり、検出結果が系統に存在するコンデンサ設備の入，切等の影響を受けないようにするため、例えば、２＜Ｍ＜３の２次高調波と３次高調波との間の次数間高調波の電流である。
【００３３】
具体的には、次数間高調波の注入チャンネルの周波数単位間隔であるベース周波数をｆ０，系統基本波周波数（６０Ｈｚ又は５０Ｈｚ）をｆ１とし、注入次数Ｍの次数間高調波の周波数（注入周波数）をｆＭとすると、ｆ１＝６０Ｈｚの系統において、例えば、ｆ０＝２Ｈｚに設定し、１２０（＝２・６０）Ｈｚ＜ｆＭ（＝１２０＋２・Ａ）＜１８０（＝３・６０）Ｈｚ，（１≦Ａ≦ｋ（ｋ＝ｆ１／ｆ０＝６０／２＝３０））の単相電流を受電点Ｐから系統５に注入する。
【００３４】
このとき、Ａ＝１５の注入次数Ｍ＝２．５（＝ｎ＋（Ａ／３０））の次数間高調波の電流であれば、その注入周波数ｆＭは１５０Ｈｚである。
【００３５】
つぎに、受電点Ｐからみた系統上流側は、一般に、系統の線路インピーダンス等の抵抗分が極めて小さいため、図２の単結線の等価回路図に示すように、系統電源Ｓと受電点Ｐとの間に、それぞれ誘導性リアクタンスの電源インピーダンスＺＳ，配電線路インピーダンスＺＬを直列接続した回路で表わされる。
【００３６】
この図２の等価回路図において、インピーダンスＺＳ，ＺＬのリアクタンスをＸＳ，ＸＬとし、受電点Ｐからみた系統上流側のインピーダンスをＺＰとすれば、ＺＳ＝ｊＸＳ，ＺＬ＝ｊＸＬであり、ＺＰ＝ＲＰ＋ｊＸＰ≒ｊ（ＸＳ＋ＸＬ），（ＲＰ：抵抗分，ＸＰ：リアクタンス分）となる。
【００３７】
また、アドミタンスＹＰ（＝１／ＺＰ）は、ＹＰ＝ＧＰ＋ｊ（−ＢＰ）≒−ｊ（１／（ＸＳ＋ＸＬ）），（ＧＰ：コンダクタンス分，（−ＢＰ）：サセプタンス分）となる。
【００３８】
一方、受電点Ｐからみた系統下流側（負荷側）のインピーダンスは、上流側（電源側）より十分に大きい。
【００３９】
そして、受電点Ｐの注入次数Ｍの次数間高調波の電圧（注入電圧）をＶｉＭ  ，注入電流をＩｉＭとし、注入次数Ｍの次数間高調波についてのインピーダンスＺＰ  ，アドミタンスＹＰをＺＰＭ，ＹＰＭとすれば、他の需要家等の影響がない理想的な場合、２ＰＭ，ＹＰＭは、それぞれＺＰＭ＝ＶｉＭ／ＩｉＭ，ＹＰＭ＝ＩｉＭ／ＶｉＭから求まり、実用上は、インピーダンスＺＰＭとして、そのリアクタンス分ＸＰＭを、ＶｉＭ／ＩｉＭの虚数部Ｉｍａｇｅ（ＶＭ／ＩＭ）から求めればよく、アドミタンスＹＰＭとして、そのサセプタンス（−ＢＰＭ）を、ＩｉＭ／ＶｉＭの（−虚数部Ｉｍａｇｅ（ＩｉＭ／ＶｉＭ））から求めれば十分である。
【００４０】
このとき、系統５の３相を相順のａ，ｂ，ｃとし、電源部１７からｂ相，ｃ相間に次数間高調波の単相電流を注入する場合、注入次数Ｍの次数間高調波についての受電点Ｐからみた系統上流側の３相等価回路は、図３のΔ形の３相負荷回路２６で表わされる。
【００４１】
そのため、次数間高調波の単相電流を注入すると、系統上流側には、その次数間高調波の３相電流が流れる。
【００４２】
そして、注入次数Ｍの次数間高調波についての負荷回路２６のΔ結線されたアドミタンス（インピーダンス）をＹａｂ（Ｚａｂ），Ｙｂｃ（Ｚｂｃ），Ｙｃａ（Ｚｃａ）とし、各相の線電圧をＶａ，Ｖｂ，Ｖｃ，線電流をＩａ，Ｉｂ，Ｉｃとすると、それぞれ対称成分に分解して得られる対称座標法での零相分Ｙ０，Ｖ０，Ｉ０，正相分Ｙ１，Ｖ１，Ｉ１，逆相分Ｙ２，Ｖ２，Ｉ２は、つぎのようにして求めることができる。
【００４３】
なお、電圧Ｖａ〜Ｖｃ，電流Ｉａ〜Ｉｃ等は、実数部と虚数部とを有するベクトル値である。
【００４４】
そして、電圧Ｖａ〜Ｖｃはつぎの数１の式で表され、式中のαはα＝ｅｘｐ（ｊ・（２π／３））の定数である。
【００４５】
【数１】
0003601517.tif 000002
【００４６】
そのため、各対称成分Ｖ０〜Ｖ２は、電圧Ｖａ〜Ｖｃを成分Ｖ０〜Ｖ２に分解するつぎの数２の式の演算から求めることができる。
【００４７】
【数２】
0003601517.tif 000003
【００４８】
同様に、各対称成分Ｉ０〜Ｉ２は、電流Ｉａ〜Ｉｃを成分Ｉ０〜Ｉ２に分解するつぎの数３の式の演算から求めることができる。
【００４９】
【数３】
0003601517.tif 000004
【００５０】
さらに、対称成分Ｖ０〜Ｖ２，Ｉ０〜Ｉ２が求まれば、対称成分Ｙ０〜Ｙ２は、つぎの数４の式の演算から求めることができる。
【００５１】
【数４】
0003601517.tif 000005
【００５２】
そして、配電系統のような非接地系統にあっては、零相分Ｉ０  が０になることから、数４の式の演算により、受電点Ｐからみた注入次数Ｍの次数間高調波のアドミタンスＹｐＭとして、正相アドミタンスＹ１＝Ｉ１／Ｖ１，逆相アドミタンスＹ２＝Ｉ２／Ｖ２が得られる。
【００５３】
このとき、正相分Ｙ１，逆相分Ｙ２は系統停止によって共に変化し、この変化を監視すれば、系統停止時の分散電源６の単独運転を検出することができる。
【００５４】
ところで、アドミタンスＹ１，Ｙ２はいずれもサセプタンス分が大きく、実用上は、アドミタンスＹ１，Ｙ２として、それぞれの虚数部―Ｉｍａｇｅ（Ｉ１／Ｖ１），−Ｉｍａｇｅ（Ｉ２／Ｖ２）を求めればよい。
【００５５】
つぎに、系統５の各需要家の負荷設備の入，切等により、そのインピーダンス（アドミタンス）が変化して系統基本波が変動し、系統特性が変わるときは、注入次数の次数間高調波に対してそのような負荷設備等が例えば図２に示す変動発生源ＮＧとして作用する。
【００５６】
そして、図２に示すように、受電点Ｐから系統５に注入される注入次数Ｍの次数間高調波の電流ＩｉＭ  に対して、変動発生源ＮＧから等価的に系統５に注入される注入次数Ｍの次数間高調波の電流をＩｇＭとする。
【００５７】
この場合、変動時の受電点Ｐの注入電圧ＶｉＭ  は、変動発生源ＮＧがないため、つぎの数５の式の演算から求まる。
【００５８】
【数５】
ＶｉＭ＝ＺＰＭ・ＩｉＭ
【００５９】
そして、変動前のサセプタンス分を（−ＢＰＭ）とすれば、この変動前のサセプタンス（−ＢＰＭ）は、つぎの数６の式の演算から求まる（Ｉｍａｇｅ：虚数部の関数）。
【００６０】
【数６】
−ＢＰＭ＝−Ｉｍａｇｅ（ＩｉＭ／ＶｉＭ）
【００６１】
一方、変動により電流ＩｇＭが注入されると、この電流ＩｇＭもインピーダンスＺｓ  ，ＺＬを流れるため、電流ＩｇＭによって、つぎの数７の式で示される注入次数Ｍの次数間高調波の電圧ＶｇＭが発生する。
【００６２】
【数７】
ＶｇＭ＝ＺＰＭ・ＩｇＭ
【００６３】
そして、インピーダンスＺＰＭに電圧ＶｉＭ，ＶｇＭの直列電圧が印加された状態になり、このとき注入電流はＩｉＭであるから、変動後のサセプタンス分を（−ＢＰＭ’）とすれば、変動後のサセプタンス（−ＢＰＭ’）は、つぎの数８の式の演算から求まる。
【００６４】
【数８】
−ＢＰＭ’＝−Ｉｍａｇｅ（ＩｉＭ／（ＶｉＭ＋ＶｇＭ））
【００６５】
そして、数６の式の絶対値（大きさ）｜ＩｉＭ／ＶｉＭ｜と数８の式の絶対値（大きさ）｜ＩｉＭ／（ＶｉＭ＋ＶｇＭ）｜とを比較すると、数８の式の分母（ＶｉＭ＋ＶｇＭ）が数６の式の分母ＶｇＭより大きくなることから、変動後のサセプタンス（−ＢＰＭ’）は変動前のサセプタンス（−ＢＰＭ）より小さくなり、系統停止時の値に近づく。
【００６６】
そして、系統基本波の変動により、受電点Ｐからみた系統上流側の注入次数Ｍの次数間高調波についてのサセプタンス分が（−ＢＰＭ）から（−ＢＰＭ’）に減少変化すると、系統停止でないにもかかわらず、分散電源６の単独運転を誤検出する。
【００６７】
ところで、この誤検出は、前記したように、注入次数Ｍの次数間高調波の電流が系統基本波の変動の影響を受けてＩｉＭからＩｉＭ＋ＩｇＭに変化することで発生する。
【００６８】
一方、単相注入の場合、前記数２〜数４の式から、次数間高調波の電流、電圧の大きさがほぼ等しくなることが、次のようにしてわかる。
【００６９】
例えば、ａ相、ｂ相間に単相注入する場合は、ａ相、ｂ相の線電流は位相が反転しており、ｃ相の線電流は零である。この時、数３の式にこの条件をあてはめると、正相電流と逆相電流は、位相は異なるが、大きさは等しくなる。
【００７０】
また、電圧に関しては、数４より、正相アドミタンスと逆相アドミタンスが等しければ、正相電圧と逆相電圧の大きさは等しくなる事がわかる。なぜなら、単相注入の場合、正相電流と逆相電流の大きさは等しいからである。また、受電点Ｐからみた系統の３相負荷回路は、概ね３相間でバランスしていると考えてよいため、アドミタンス（インピーダンス）の正相と逆相はほぼ等しい。
【００７１】
以上より、単相注入の場合は次数間高調波の電圧、電流の大きさがほぼ等しくなるといえる。
【００７２】
そして、変圧器１９，変流器２０の３相の計測信号のフーリエ変換で系統５の注入次数Ｍの次数間高調波の電圧、電流を検出（抽出）した場合、その検出における系統基本波の影響は、具体的には、正相分と逆相分とに分けて考案する必要がある。
【００７３】
つぎに、次数間高調波抽出における系統基本波の変動の影響について、具体的に説明する。
まず、以下の説明においては、前記のようにベース周波数，系統基本波周波数（商用基本波周波数）をｆ０，ｆ１とし、注入次数Ｍの代わりに、ベース周波数ｆ０を基準（基本波）にした整数次数ｍ＝ｆＭ／ｆ０を使用し、周波数ｆＭ＝ｆｍとする。
【００７４】
（イ）系統基本波の変動がない場合
この場合、周波数ｆ０，ｆ１，ｆＭの各周波数ω０，ω１，ωｍは、ω０＝２πｆ０，ω１＝２πｆ１＝ｋω０，ωｍ＝２πｆｍ＝ｍω０，（ｍ，ｋ：ｍ＞ｋの整数）になる。
【００７５】
そして、次数ｍの次数間高調波の電流ＩｉＭの注入により、変流器２０の計測信号に基づく時刻ｔの３相の計測電流（線電流）を、ベクトル値の電流Ｉａ（ｔ），Ｉｂ（ｔ），Ｉｃ（ｔ）とすると、これらの電流Ｉａ（ｔ）〜Ｉｃ（ｔ）  は、系統基本波周波数に次数ｍの次数間高調波が重畳しているため、つぎの数９の３式で表すことができ、それぞれの右辺第１項が系統基本波の成分であり、右辺第２項が次数間高調波の成分である。
【００７６】
【数９】
Ｉａ（ｔ）＝Ｉ１ａ・ｓｉｎ（ω１ｔ＋φ１ａ）＋  Ｉｍａ・ｓｉｎ（ωｍｔ＋φｍａ）
Ｉｂ（ｔ）＝Ｉ１ｂ・ｓｉｎ（ω１ｔ＋φ１ｂ）＋  Ｉｍｂ・ｓｉｎ（ωｍｔ＋φｍｂ）
Ｉｃ（ｔ）＝Ｉ１ｃ・ｓｉｎ（ω１ｔ＋φ１ｃ）＋  Ｉｍｃ・ｓｉｎ（ωｍｔ＋φｍｃ）
【００７７】
なお、数９の各式中のＩ１ａ，Ｉ１ｂ，Ｉ１ｃ，Ｉｍａ，Ｉｍｂ，Ｉｍｃは最大値であり、φ１ａ，φ１ｂ，φ１ｃ，φｍａ，φｍｂ，φｍｃは初期位相である。
【００７８】
そして、この９式の電流Ｉａ（ｔ）〜Ｉｃ（ｔ）  から注入次数Ｍの次数間高調波を検出（抽出）する場合、一般につぎのフーリエ変換が用いられる。
【００７９】
いま、ベクトル値の電流Ｉ（ｔ）を周期信号とし、その基本の角周波数をω０とすると、電流Ｉ（ｔ）に含まれるｎ次高調波の成分（高調波成分）は、その正弦成分をＩｎＳ，余弦成分をＩｎＣとして、つぎの数１０の２式のフーリエ変換式から算出される。
【００８０】
【数１０】
0003601517.tif 000006
なお、式中のωｎ  はｎ次高調波の角周波数、Ｓ字状の記号は積分記号である。
【００８１】
そして、数１０の２式に基づくｎ次高調波の成分をＩ（ｎ）とするとこの高調波成分Ｉ（ｎ）は、つぎの数１１の複素数表現の式で表される。
【００８２】
【数１１】
0003601517.tif 000007
【００８３】
ところで、電流Ｉ（ｔ）＝Ｉ・ｓｉｎ（ωｔ＋φ）のときには、電流Ｉ（ｔ）がつぎの数１２のベクトル式で表される。
【００８４】
【数１２】
Ｉ（ｔ）＝Ｉ・ｓｉｎ（ωｔ＋φ）＝（１／２ｊ）・［Ｉ・ｅｘｐ（ｊ（ωｔ＋φ））−Ｉ・ｅｘｐ（−ｊ（ωｔ＋φ））］
【００８５】
そのため、つぎの数１３の式が成立する。
【００８６】
【数１３】
0003601517.tif 000008
【００８７】
そして、数９の式の電流Ｉａ（ｔ）から次数ｍの次数間高調波を検出（抽出）する場合は、つぎの数１４の式の演算を行えばよい。
【００８８】
【数１４】
0003601517.tif 000009
【００８９】
この数１４の式の演算を実際に行うと、数１３の式の関係を用いることにより、つぎの数１５の式が得られる。
【００９０】
【数１５】
0003601517.tif 000010
【００９１】
さらに、数１３の式において、ω１−ωｍ＝（ｋ−ｍ）ω０，ω１＋ωｍ＝（ｋ＋ｍ）ω０，ωｍ−ωｍ＝０，ωｍ＋ωｍ＝２ｍω０であり、ｋ−ｍ，ｋ＋ｍ，２ｍが整数になるため、数１５の式において、つぎの数１６の４式が成立する。
【００９２】
【数１６】
0003601517.tif 000011
【００９３】
したがって、数１５の式の演算結果は、Ｉｍａ・ｅｘｐ（ｊφｍａ）となる。
そして、数９の式の電流Ｉｂ（ｔ），Ｉｃ（ｔ）についても同様の演算を行うことにより、電流Ｉａ（ｔ），Ｉｂ（ｔ），Ｉｃ（ｔ）の次数ｍの次数間高調波の成分Ｉａ（ｍ），Ｉｂ（ｍ），Ｉｃ（ｍ）は、系統基本波の成分に変動がない場合、つぎの数１７の３式のフーリエ変換の演算により、系統基本波の影響を受けることなく検出することができる。
【００９４】
【数１７】
0003601517.tif 000012
【００９５】
（ロ）系統基本波が変動した場合
つぎに、時刻ｔ０  に負荷変動等が発生して系統基本波が変動したとすると、ｔ０前（ｔ＜ｔ０）には前記数９の３式で表される各相の電流Ｉａ（ｔ）〜Ｉｃ（ｔ）が、ｔ０以後（ｔ≧ｔ０）は、つぎの数１８の３式で表されるようになる。
【００９６】
【数１８】
Ｉａ（ｔ）＝Ｉ１ａ’・ｓｉｎ（ω１ｔ＋φ１ａ’）＋Ｉｍａ・ｓｉｎ（ωｍｔ＋φｍａ）
Ｉｂ（ｔ）＝Ｉ１ｂ’・ｓｉｎ（ω１ｔ＋φ１ｂ’）＋Ｉｍｂ・ｓｉｎ（ωｍｔ＋φｍｂ）
Ｉｃ（ｔ）＝Ｉ１ｃ’・ｓｉｎ（ω１ｔ＋φ１ｃ’）＋Ｉｍｃ・ｓｉｎ（ωｍｔ＋φｍｃ）
【００９７】
なお、数１８の各式中の右辺第１項が変動後の系統基本波の成分であり、Ｉ１ａ’，Ｉ１ｂ’，Ｉ１ｃ’は変動後の最大振幅、φ１ａ’，φ１ｂ’，φ１ｃ’は変動後の初期位相である。
また、数１８の各式の右辺第２項の次数間高調波の成分は変化しない。
【００９８】
そして、数１８の各式で表される電流Ｉａ（ｔ）〜Ｉｃ（ｔ）につき、（ｔ−２π／ω０）≦ｔ０≦ｔの範囲で前記数１７の各式のフーリエ変換を施して電流Ｉａ（ｍ）〜Ｉｃ（ｍ）を求める。
このとき、例えばａ相の電流Ｉａ（ｍ）は、つぎの数１９の式で表される。
【００９９】
【数１９】
0003601517.tif 000013
【０１００】
さらに、この数１９の式に数１２の式の関係を適用すると、ω１＝ｋω０，ωｍ＝ｍω０，ｋ−ｍ，ｋ＋ｍが整数であることから、数１９の式の最右辺の第１項につき、つぎの数２０の式が成立する。
【０１０１】
【数２０】
0003601517.tif 000014
【０１０２】
また、数１９の式の最右辺の第２項につき、つぎの数２１の式が成立する。
【０１０３】
【数２１】
0003601517.tif 000015
【０１０４】
そのため、数１９の式の電流Ｉａ（ｍ）は、つぎの数２２の式に示すようになる。
【０１０５】
【数２２】
0003601517.tif 000016
【０１０６】
この数２２の式において、Ｉ１ａ’・ｅｘｐ（ｊφ１ａ’）−Ｉ１ａ・ｅｘｐ（ｊφ１ａ）は系統基本波の電流成分の変化による差分ベクトルであり、これをΔＩａとし、その共役量を＊を付けてΔＩａ＊とすると、つぎの数２３の２式が成立する。
【０１０７】
【数２３】
ΔＩａ＝Ｉ１ａ’・ｅｘｐ（ｊφ１ａ’）−Ｉ１ａ・ｅｘｐ（ｊφ１ａ）
ΔＩａ＊＝Ｉ１ａ’・ｅｘｐ（−ｊφ１ａ’）−Ｉ１ａ・ｅｘｐ（−ｊφ１ａ）
【０１０８】
そして、ｂ相の電流Ｉｂ（ｍ），ｃ相の電流Ｉｃ（ｍ）についても同様であることから、ｔ−２π／ω０≦ｔ０≦ｔにおける、次数ｍの次数間高調波成分の検出量（抽出量）は、つぎの数２４の３式に示すようになる。
【０１０９】
【数２４】
0003601517.tif 000017
【０１１０】
なお、数２４の各式中の＊を付したものは共役量を示し、ζ，ξは、つぎの数２５の２式それぞれに示す成分である。
【０１１１】
【数２５】
0003601517.tif 000018
【０１１２】
そして、数２４の式の各相の電流Ｉａ（ｍ），Ｉｂ（ｍ），Ｉｃ（ｍ）から、次数ｍの次数間高調波を対称成分に分解したときの正相分（正相検出量）Ｉ１（ｍ），逆相分（逆相検出量）Ｉ２（ｍ）を求めると、一般に３相電流Ｉａ，Ｉｂ，Ｉｃに基づく正相分Ｉ１，逆相分Ｉ２がＩ１＝（１／３）（Ｉａ＋αＩｂ＋α２ＩＣ），Ｉ２＝（１／３）（Ｉａ＋α２Ｉｂ＋αＩＣ），（但し、α＝ｅｘｐ（ｊ２／３π），α２＝ｅｘｐ（ｊ４／３π））になることから、つぎの数２６の２式に示すようになる。
【０１１３】
【数２６】
Ｉ１（ｍ）＝ζ・（ΔＩ１１）−ξ・（ΔＩ１２）＊＋Ｉｍ１
Ｉ２（ｍ）＝ζ・（ΔＩ１２）−ξ・（ΔＩ１１）＊＋Ｉｍ２
【０１１４】
この数２６の２式において、（ΔＩ１１）＊，（ΔＩ１２）＊は（ΔＩ１１），（ΔＩ１２）の共役量を示し、右辺の成分ζ，ξを含む第１，第２項が系統基本波の変動の影響により生じる誤差成分ΔＩｍ１，ΔＩｍ２であり、右辺の第３項が本来の注入成分である。
【０１１５】
そして、成分ζ，ξは、数２５の２式の演算により、つぎの数２７の２式に示すようになる。
【０１１６】
【数２７】
0003601517.tif 000019
【０１１７】
この数２７の２式に基づき、成分ζ，ξの大きさ｜ζ｜，｜ξ｜は、つぎの数２８の２式に示すようになる。
【０１１８】
【数２８】
｜ζ｜≦｜（１／（ｋ−ｍ）π）・ｓｉｎ［（ｋ−ｍ）ω０（ｔ−ｔ０）／２］｜≦（１／（ｍ−ｋ）π）
｜ξ｜≦｜（１／（ｋ＋ｍ）π）・ｓｉｎ［（ｋ＋ｍ）ω０（ｔ−ｔ０）／２］｜≦（１／（ｍ＋ｋ）π）
【０１１９】
したがって、前記数２６の２式で示される正相分Ｉ１（ｍ），逆相分Ｉ２（ｍ）は系統基本波の変動の影響を受け、その変動による誤差成分の大きさ｜ΔＩｍ１｜，｜ΔＩｍ２｜
は、前記数２８の２式に基づき、つぎの数２９の２式それぞれで示される。
【０１２０】
【数２９】
0003601517.tif 000020
【０１２１】
この数２９の２式の最右辺に示すように、誤差成分｜ΔＩｍ１｜，｜ΔＩｍ２｜は、誤差成分｜ΔＩｍ１｜の式中の分子のｋ（｜ΔＩ１１｜―｜ΔＩ１２｜）が、誤差成分｜ΔＩｍ２｜では、ｋ（｜ΔＩ１２｜―｜ΔＩ１１｜）になっている点で異なる。
【０１２２】
そして、前記数２９の２式中の｜ΔＩ１１｜，｜ΔＩ１２｜は、｜ΔＩ１１｜が系統基本波の変動により生じたその正相分の変化量（変動量）であり、｜ΔＩ１２｜が系統基本波の変動により生じたその逆相分の変化量（変動量）である。
【０１２３】
また、前記数２９の２式のｍ，ｋは、ｆ０＝２Ｈｚ，ｆ１＝６０Ｈｚ，ｆｍ＝１５０Ｈｚであれば、ｍ＝１５０／２＝７５，ｋ＝６０／２＝３０であり、このとき、ｍ−ｋ＝４５，ｍ＋ｋ＝１０５になる。
【０１２４】
以上より、系統５の負荷変動等で系統基本波が変動した場合は、数２４の３式で示される次数ｍの次数間高調波の電流Ｉａ（ｍ），Ｉｂ（ｍ），Ｉｃ（ｍ）が系統基本波の変動の影響を受け、それらを対称成分に分解すると、数２６の２式で示される正相分Ｉ１（ｍ），逆相分Ｉ２（ｍ）は、単相注入に基づく本来の成分Ｉｍ１，Ｉｍ２だけでなく、数２９に示す系統基本波成分の変動に基づくそれぞれの誤差成分ΔＩｍ１，ΔＩｍ２それぞれを含むことが判明した。
【０１２５】
さらに、誤差成分ΔＩｍ１，ΔＩｍ２により、例えば、変動後のサセプタンス（−ＢＰｍ’）（＝−ＢＰＭ）’）が変動前のサセプタンス（−ＢＰｍ）（＝−ＢＰＭ）から動し、このときサセプタンスＢｉｍ’の変動の大きさがその正相分と逆相分とで異なることも判明した。
【０１２６】
（ハ）系統基本波の影響
まず、前記数２６の２式で示される正相分Ｉ１（ｍ），逆相分Ｉ２（ｍ）において、本来の注入成分の正相量Ｉｍ１と逆相量Ｉｍ２とは、単相注入であるため、その大きさ｜Ｉｍ１｜，｜Ｉｍ２｜は等しく、｜Ｉｍ１｜＝｜Ｉｍ２｜である。
【０１２７】
つぎに、誤差成分ΔＩｍ１，ΔＩｍ２の大きさ｜ΔＩｍ１｜，｜ΔＩｍ２｜は、前記数２９の２式で示され、系統基本波の正相分の変化量｜ΔＩ１１｜と逆相分の変化量｜ΔＩ１２｜との大小関係によって異なる。
【０１２８】
そして、正相分Ｉ１（ｍ）と逆相分Ｉ２（ｍ）の系統基本波の変動の影響は、変化量ΔＩ１１，ΔＩ１２の大小関係によって変化し、つぎの（ｉ），（ｉｉ）のようになる。
【０１２９】
（ｉ）｜ΔＩ１１｜≦｜ΔＩ１２｜のときは、正相分Ｉ１（ｍ）の誤差成分ΔＩｍ１が逆相分Ｉ２（ｍ）の誤差成分ΔＩｍ２より小さく、正相分Ｉ１（ｍ）の方が系統基本波の変動の影響を受けにくい。
（ｉｉ）  ｜ΔＩ１１｜＞｜ΔＩ１２｜のときは、逆相分Ｉ２（ｍ）の誤差成分ΔＩｍ２が正相分Ｉ１（ｍ）の誤差成分ΔＩｍ１より小さく、逆相分Ｉ２（ｍ）の方が系統基本波の変動の影響を受けにくい。
【０１３０】
したがって、｜ΔＩ１１｜≦｜ΔＩ１２｜であれば、正相分Ｉ１（ｍ）を用いて次数ｍの次数間高調波についての検出インピーダンス又は検出アドミタンスの正相分を検出・監視し、その変化から系統停止時の分散電源６の単独運転を検出し、｜ΔＩ１１｜＞｜ΔＩ１２｜であれば、逆相分Ｉ２（ｍ）を用いて次数ｍの次数間高調波についての検出インピーダンス又は検出アドミタンスの正相分を検出・監視し、その変化から系統停止時の分散電源６の単独運転を検出すれば、系統基本波の変動の影響が極めて少なくなる。
【０１３１】
そのため、検出アドミタンスの変化から分散電源６の単独運転を検出する場合、系統基本波の正相変化量ΔＩ１１と逆相変化量ΔＩ１２の大小を比較し、｜ΔＩ１１｜≦｜ΔＩ１２｜であれば、検出アドミタンスＹＰｍ（＝ＹＰＭ）を正相アドミタンスＹ１とし、その変化から分散電源６の運転停止を検出し、｜ΔＩ１１｜＞｜ΔＩ１２｜であれば、検出アドミタンスＹＰｍを逆相アドミタンスＹ２とし、その変化から分散電源６の運転停止を検出する。
【０１３２】
具体的には、系統５について、事前に正相変化量ΔＩ１１と逆相変化量ΔＩ１２を測定し、その大小の比較結果から、検出アドミタンスＹＰｍを、正相アドミタンスＹ１又は逆相アドミタンスＹ２に設定する。
【０１３３】
そして、｜ΔＩ１１｜≦｜ΔＩ１２｜であって、検出アドミタンスＹＰｍが正相アドミタンスＹ１に設定される場合、演算処理部２５は、まず、変圧器１９，変換器２０の各相の計測信号に含まれた注入次数ｍの次数間高調波の計測電圧，計測電流それぞれを抽出するため、前記のフーリエ変換のデジタルフィルタ演算として、例えば回帰形ＤＦＴ演算を実行する。
【０１３４】
このとき、過去ＮサンプリングのデータでＤＦＴ演算するため、各相のサンプリングされた計測信号の電圧，電流を、電圧Ｖ（ｑ），電流Ｉ（ｑ），（ｑ：０，１，２，…，Ｎ−１）とすると、直前の変化結果を利用するつぎの数３０の２式の演算から、各相の注入次数Ｍの次数間高調波の計測電圧Ｖｍ（ｑ），計測電流ｍ（ｑ）を求める。
【０１３５】
【数３０】
Ｖｍ（ｑ）＝（２／Ｎ）・（Ｖｍ（ｑ−１）−Ｖ（ｑ−Ｎ）＋Ｖ（ｑ））・ｘ−１
Ｉｍ（ｑ）＝（２／Ｎ）・（Ｉｍ（ｑ−１）−Ｉ（ｑ−Ｎ）＋Ｉ（ｑ））・ｘ−１
【０１３６】
なお、データ数Ｎは、例えば系統基本波１サイクル当り６４サンプリングで３２波長分とすると、（６４×３２＝）４０９６個である。
また、数５の式中のｘはｘ＝ｅｘｐ（−ｊ２πｍ／Ｎ）のパラメータである。
【０１３７】
つぎに、前記数５の式の演算により抽出された系統各相の次数ｍの次数間高調波の電圧，電流を、Ｖａｍ（ｑ），Ｖｂｍ（ｑ），Ｖｃｍ（ｑ），Ｉａｍ（ｑ），Ｉｂｍ（ｑ），Ｉｃｍ（ｑ）とすると、正相分Ｙ１，Ｖ１，Ｉ１をＹｍ１（ｑ），Ｖｍ１（ｑ），Ｉｍ１（ｑ）として、つぎの数３１の２式の正相演算により、電圧Ｖｍ（ｑ），電流Ｉｍ（ｑ）それぞれの正相分Ｖｍ１（ｑ），Ｉｍ１（ｑ）を求める。
【０１３８】
【数３１】
Ｖｍ１（ｑ）＝（Ｖａｍ（ｑ）＋α・Ｖｂｍ（ｑ）＋α２・Ｖｃｍ（ｑ））／３
Ｉｍ１（ｑ）＝（Ｉａｍ（ｑ）＋α・Ｉｂｍ（ｑ）＋α２・Ｉｃｍ（ｑ））／３
【０１３９】
さらに、正相アドミタンスＹｍ１（ｑ）として、この形態では、正相分Ｖｍ１（ｑ），Ｉｍ１（ｑ）に基づくつぎの数３２の式の演算から、そのサセプタンス分（−Ｂｍ１（ｑ））を求める。
【０１４０】
【数３２】
（−Ｂｍ１（ｑ））＝−Ｉｍａｇｅ（Ｉｍ１（ｑ）／Ｖｍ１（ｑ））
【０１４１】
そして、このサセプタンス分（−Ｂｍ１（ｑ））の算出をくり返すとともに、算出したサセプタンス分（−Ｂｍ１（ｑ））と設定された判定値とを比較して、サセプタンス分（−Ｂｍ１（ｑ））の変化を監視する。
【０１４２】
すなわち、算出したサセプタンス分（−Ｂｍ１（ｑ））が、設定された継続整定値の期間連続して判定値より小さくなり、受電点Ｐからみた系統上流側のインピーダンスが開放インピーダンスに増大変化するか否かを監視し、この増大変化から系統停止による分散電源６の単独運転を検出する。
【０１４３】
この正相変化による単独運転の検出処理は、具体的には、次数ｍを次数Ｍとして、演算処理部２５が例えば図４に示す正相検出処理のフローチャートの各ステップＳ１〜Ｓ６を実行することで実現される。
【０１４４】
一方、｜ΔＩ１１｜＞｜ΔＩ１２｜であって、検出アドミタンスＹＰｍ（＝ＹＰＭ）が逆相アドミタンスＹ２に設定される場合は、演算処理部２５により、前記の数３１の正相演算の代わりにつぎの数３３の逆相演算を実行し、電圧Ｖｍ（ｑ），電流Ｉｍ（ｑ）それぞれの逆相分Ｖｍ２（ｑ），Ｉｍ２（ｑ）を求める。
【０１４５】
【数３３】
Ｖｍ２（ｑ）＝（Ｖａｍ（ｑ）＋α２・Ｖｂｍ（ｑ）＋α・Ｖｃｍ（ｑ）／３
Ｉｍ２（ｑ）＝（Ｉａｍ（ｑ）＋α２・Ｉｂｍ（ｑ）＋α・Ｉｃｍ（ｑ）／３
【０１４６】
さらに、逆相分Ｖｍ２（ｑ），Ｉｍ２（ｑ）に基づくつぎの数３４の式の演算から、逆相アドミタンスＹＩ２（ｑ）として、そのサセプタンス分（−Ｂｍ２（ｑ））を求める。
【０１４７】
【数３４】
（−Ｂｍ２（ｑ））＝−Ｉｍａｇｅ（Ｉｍ２（ｑ）／Ｖｍ２（ｑ））
【０１４８】
そして、サセプタンス分（−Ｂｍ２（ｑ））の代わりにサセプタンス分（−Ｂｍ２（ｑ））の変化を監視して分散電源６の単独運転を検出する。
【０１４９】
なお、この場合は、演算処理部２３が図中の検出処理のフローチャートの代わりに図５の逆相検出処理のフローチャートの各ステップＳ１’〜Ｓ６’を実行する。
【０１５０】
したがって、系統の負荷変動の傾向等に応じた適切な検出方法により、系統基本波の変動の影響を極力受けないようにして系統停止時の分散電源６の単独運転を確実に検出することができ、誤検出による分散電源６の系統５からの誤った切離しを防止することがてきる。
【０１５１】
そして、注入次数Ｍの次数間高調波の電流を受電点Ｐから系統５の例えばｂ相，ｃ相間に単相注入すればよいため、注入電源（電源部１７）等が３相注入の場合より極めて安価かつ小型になる。
【０１５２】
また、系統の注入次数Ｍの次数間高調波についての検出アドミタンスとして、その対称成分である正相又は逆相のアドミタンスのサセプタンス分を算出し、その変化から単独運転を検出したため、系統の相順変化等の影響を受けることなく、系統停止側の分散電源６の単独運転を検出することができ、系統の相順変化等を気にする必要がない。
【０１５３】
なお、実系統においては、例えば、単相負荷が多く、３相間の負荷変動が不規則であれば、｜ΔＩ１１｜≦｜ΔＩ１２｜，｜ΔＩ１１｜＞｜ΔＩ１２｜のいずれにもなり得ることから、この場合正相分Ｉ１（ｍ），逆相分Ｉ２（ｍ）のいずれを用いて検出してもよいが、通常は、扱う量として一般的な正相分とすることが好ましい。
【０１５４】
一方、３相負荷の同時開閉等により、３相間の負荷変動がある程度規則的であれば、｜ΔＩ１１｜＞｜ΔＩ１２｜になるため、逆相分Ｉ２（ｍ）を用いて系統基本波の変動の影響を受けないようにする。
【０１５５】
そして、３相の負荷が同時開閉されて３相同時に増，減変化するときは、前記のｆ０＝２Ｈｚ，ｆ１＝６０Ｈｚ，ｋ＝６０／２＝３０とすると、｜ΔＩ１２｜＝０の条件下、｜ΔＩｍ１｜≦（１／π）・（｜ΔＩ１１｜／（ｍ−３０）），｜ΔＩｍ２｜≦（１／π）・（｜ΔＩ１１｜／（ｍ＋３０））となり、｜ΔＩｍ２｜＜｜ΔＩｍ１｜になることから、逆相分Ｉ２（ｍ）を用いた方が系統基本波の変動の影響は少ないことは明らかである。
【０１５６】
ところで、検出アドミタンスＹＰｍとして、正相アドミタンスＹ１（＝Ｙｍ１（ｑ）），逆相アドミタンスＹ２（＝Ｙｍ２（ｑ））の実数部のコンダクタンス及び虚数部のサセプタンスの両方を求め、正相アドミタンスＹ１（＝Ｙｍ１（ｑ）），逆相アドミタンスＹ２（＝Ｙｍ２（ｑ））のベクトル変化から分散電源６の単独運転を検出するようにしてもよい。
【０１５７】
また、検出アドミタンスＹＰｍの逆数の検出インピーダンスＺＰｍの変化から、系統停止側の分散電源６の単独運転を検出する場合は、正相アドミタンスＹ１（＝Ｙｍ１（ｑ）），逆相アドミタンスＹ２（＝Ｙｍ２（ｑ））の代わりに、｜ΔＩ１１｜≦｜ΔＩ１２｜であれば、正相インピーダンスＺ１（＝Ｚｍ１（ｑ））を、Ｚｍ１（ｑ）＝Ｖｍ１（ｑ）／Ｉｍ１（ｑ）の演算から算出し、｜ΔＩ１１｜＞｜ΔＩ１２｜であれば、逆相インピーダンスＺ２（＝Ｚｍ２（ｑ））を、Ｚｍ２（ｑ）＝Ｖｍ２（ｑ）／Ｉｍ２（ｑ）の演算から算出する。
【０１５８】
さらに、サセプタンス分Ｂｍ１（ｑ），Ｂｍ２（ｑ）に対応する正相インピーダンスＺｍ１（ｑ）のリアクタンス分Ｘｍ１（ｑ），逆相インピーダンスＺｍ２（ｑ）のリアクタンス分Ｘｍ２を、つぎの数３５の式から算出し、リアクタンス分Ｘｍ１（ｑ），Ｘｍ２（ｑ）の変化を監視して検出すればよい。
【０１５９】
【数３５】
Ｘｍ１（ｑ）＝Ｉｍａｇｅ（Ｖｍ１（ｑ）／Ｉｍ１（ｑ））
Ｘｍ２（ｑ）＝Ｉｍａｇｅ（Ｖｍ２（ｑ）／Ｉｍ２（ｑ））
【０１６０】
なお、リアクタンス分Ｘｍ１（ｑ），Ｘｍ２（ｑ）の変化でなく、正相インピーダンスＺｍ１（ｑ），逆相インピーダンスＺｍ２（ｑ）のベクトル変化を監視して検出してもよい。
【０１６１】
つぎに、変化量ΔＩ１１，ΔＩ１２の大小に応じて検出アドミタンスと、正相アドミタンスＹ１又は逆相アドミタンスＹ２に設定する場合、前記のように運用前の事前の変化量ΔＩ１１，ΔＩ１２の設定に基づいて、検出アドミタンスを正相アドミタンスＹ１又は逆相アドミタンスＹ２に設定する代わりに、例えば、自動的に又は手動操作に基づき、演算処理部２５により、周期的に又は任意の時点で変流器２０の計測信号に含まれた系統基本波の電流を検出（抽出）し、検出した系統基本波の電流を対称成分に分解して、その正相分，逆相分それぞれの変化量ΔＩ１１，ΔＩ１２の大小比較により、｜ΔＩ１１｜≦｜ΔＩ１２｜，｜ΔＩ１１｜＞｜ΔＩ１２｜を判別して検出アドミタンスを正相アドミタンスＹ１又は逆相アドミタンスＹ２に変更自在に設定してもよい。
【０１６２】
この場合は、｜ΔＩ１１｜≦｜ΔＩ１２｜，｜ΔＩ１１｜＞｜ΔＩ１２｜の判別後、演算処理部２５により、図４の検出処理プログラムと図５の検出処理プログラムとを選択的に実行する。
【０１６３】
また、系統の基本波電流の逆相分の変動がその正相分の変動より小さく、注入次数の次数間高調波の逆相分が正相分より系統基本波の変動の影響を受けにくいときは、単相注入に基づく系統各相の次数間高調波の電圧，電流それぞれの逆相分に基づき、系統の注入次数の次数間高調波についての逆相インピーダンス又は逆相アドミタンスが、検出インピーダンス又は検出アドミタンスとして算出されて検出され、その変化から、系統基本波の変動の影響を受けることなく、系統停止時の分散電源の単独運転を検出することができる。
【０１６４】
そして、単相注入に基づく系統各相の注入次数の次数間高調波の電圧，電流を対称成分に分解してそれぞれの正相分，逆相分を検出するため、系統の相順変化等を気にすることなく、検出インピーダンスとしての注入次数の次数間高調波の正相インピーダンス，正相アドミタンス又は逆相インピーダンス，逆相アドミタンスを検出することができる。
【０１６５】
なお、変化量Δ｜Ｉ１１｜，｜ΔＩ１２｜の大小に応じて、検出インピーダンスを正相インピーダンスＺ１又は逆相インピーダンスＺ２に設定する場合も、同様にして変更自在に設定してもよいのは勿論である。
【０１６６】
そして、本発明は、３相系統以上の多相の種々の電力系統の分散電源の単独運転検出に適用することができ、その際、次数間高調波の単相電源は系統の任意の２相間に注入すればよい。
【０１６７】
【発明の効果】
本発明は、以下に記載する効果を奏する。
まず、例えば３相系統であっても次数間高調波を単相注入すればよく、注入電源として、３相電源より安価かつ小型の単相電源を用いることができる。
【０１６８】
つぎに、系統５の基本波電流の正相分の変動がその逆相分の変動以下で注入次数の次数間高調波の正相分が系統基本波電流の変動の影響を受けにくいときは、単相注入に基づく系統各相の注入次数の次数間高調波の電圧，電流それぞれの正相分に基づき、系統５の注入次数の次数間高調波についての正相インピーダンス又は正相アドミタンスが、検出インピーダンス又は検出アドミタンスとして算出されて検出され、その変化から、系統５の基本波の変動の影響を受けることなく、系統停止時の分散電源６の単独運転を検出することができる。
【０１６９】
したがって、コストダウン及び小型化を図った次数間高調波の単相注入により、系統５の相順変化等を気にすることなく、系統基本波の変動の影響を極力排除して、系統５の注入次数の次数間高調波のインピーダンス又はアドミタンスの変化から、系統停止時の分散電源６の単独運転を確実に検出する。
【図面の簡単な説明】
【図１】本発明の実施の１形態の単線結線図である。
【図２】図１の次数間高調波についての等価回路図である。
【図３】図２の単相注入時の３相結線の等価回路図である。
【図４】図１の正相検出処理のフローチャートである。
【図５】図１の逆相検出処理のフローチャートである。
【符号の説明】
１，５  系統
６  分散電源
１６  電流注入装置
２１  系統停止検出処理装置
【図１】
0003601517.tif 000021
【図２】
0003601517.tif 000022
【図３】
0003601517.tif 000023
【図４】
0003601517.tif 000024
【図５】
0003601517.tif 000025