Instantaneous Frequency in Power Systems using the Teager-Kaiser Energy Operator

技術要約：テージェル・カイザー・エネルギー演算子を用いた電力系統における瞬時周波数推定

問題提起
瞬時周波数（IF）の推定は、電力系統における非定常現象を特徴付ける上で極めて重要である。近年の解釈では、電力系統の周波数は複素周波数キネマティクスおよび電圧軌道の幾何学に関連付けられているが、従来の推定法は、テージェル・カイザー・エネルギー演算子（TKEO）に基づくエネルギー分離アルゴリズム（ESA）に依存することが多い。これらの伝統的な推定法は、狭帯域信号および「スローエンベロープ（緩やかな包絡線）」条件（振幅変調および周波数変調の速度が搬送波周波数に対して小さい状態）を前提としている。しかし、電力系統においては、特に過渡現象時や電圧大きさが大きく変動する際に、これらの仮定が頻繁に崩れる。その結果、従来の推定法は電圧大きさの変化に伴う動的項を無視するため、偏った（バイアスのある）周波数推定値を生じさせる可能性がある。

手法
著者らは、TKEOを複素信号へと拡張（CTKEO）し、動的信号の恒等式を適用することで導出された、新しい局所的IF推定法を提案している。その手法は以下の通りである。

1. 複素TKEOの定式化： 実数値の標準的なTKEOを、対称双線型形式を用いて複素信号へと一般化する。複素信号 $\bar{x}$ に対して、演算子は $\Psi_C(\bar{x}) = |\bar{x}'|^2 - \text{Re}(\bar{x}\bar{x}''^*)$ と定義される。この定式化はデカルト成分に対して加算的であり、一定の位相シフトに対して不変であり、狭帯域の仮定を必要としない。
2. 厳密な恒等式の導出： 複素信号をその大きさ $|\bar{x}|$ と位相 $\phi$ を用いて表現し、複素周波数 $\bar{\eta} = \rho + j\omega$ （ここで $\rho$ は大きさの変化率、$\omega$ は角周波数）の定義を利用することで、著者らは厳密な連続時間恒等式を導出する：
   $$\Psi_C(\bar{x}) = |\bar{x}|^2 (2\omega^2 - \rho')$$
   この方程式は、回転エネルギー成分と、大きさのダイナミクスに関連する移動項（$\rho'$）を明示的に分離している。
3. 提案される推定器：
   • $\omega_{TV}$（時変型）： この恒等式を $\omega$ について解くことで、大きさのダイナミクスの補正項（$\rho'$ および $\rho^2$）を保持した推定器が得られる。これにより、電圧大きさが急速に変化する条件下でも、バイアスが補正された推定が可能となる。
   • $\omega_{TI}$（時不変型）： 大きさが一定である、あるいは緩やかに変化するという仮定に基づき、補正項を無視することで簡略化された近似が導かれる。これは $\omega_{TI} = \sqrt{\Psi_C(\bar{x}) / (2|\bar{x}|^2)}$ となり、計算がより単純でノイズに対する感度が低いが、電圧大きさの大きな過渡現象時にはバイアスが生じる可能性がある。
4. 幾何学的解釈： 提案された推定器は複素周波数キネマティクスと一致しており、捩れ（torsion）がゼロであると仮定した場合の $(v_\alpha, v_\beta)$ 平面における電圧軌道の曲率に関連した幾何学的解釈を許容する。

主な貢献

• 新しい表現式： 本論文は、エンベロープの曲率項を明示的に扱うCTKEOベースのIF表現を導入している。これにより、大きさのダイナミクスを無視することによる従来の推定法のバイアスを修正している。
• バイアス補正： 古典的なESAとは異なり、提案された定式化（$\omega_{TV}$）は「スローエンベロープ」近似に依存しないため、連続時間における滑らかで非ゼロの二次元電圧軌跡に対して理論的に厳密である。
• 幾何学的整合性： 提案された推定器は複素周波数キネマティクスと整合しており、明示的な位相アンラップ（unwrapping）を必要とせずに周波数の幾何学的解釈を提供する。
• 実用的な近似： 大きさの変化が緩やかなシナリオに対して、計算効率の高い近似（$\omega_{TI}$）を提供し、精度とノイズ感度のトレードオフを実現している。

結果
提案手法は、2つのケーススタディを通じて検証された：

1. EMTシミュレーション（IEEE 39バス系統）： 三相故障のシミュレーションにおいて、ノイズを含む平衡条件下では、提案された推定器（$\omega_{TV}$ および $\omega_{TI}$）は幾何学的参照周波数（$\omega_v$）と密接に一致した。不平衡シナリオにおいては、提案された推定器は参照値と比較して振動が小さく、不平衡によって導入される曲率の歪みの影響を受けにくいことが示された。
2. フィールド計測（太陽光発電所）： スペインの太陽光発電所から得られた実際の電圧データを用いたテストでは、両方の推定器が幾何学的参照の主要な傾向に従うことが示された。高ノイズや微分感度が高い区間（$\omega_{TV}$ は二階微分を必要とするため）では乖離が見られたものの、結果として、大きさの変化が本質的に緩やかであるため、簡略化された $\omega_{TI}$ が現実的な運用条件下において、フル・タイムバリアント（時変）推定値に実用的に近い値を示すことが確認された。

合成された振幅変調信号に関する定量分析（表I）では、サンプリング周波数が低い場合においても、時変型推定器（$\omega_{TV}$）が時不変型推定器（$\omega_{TI}$）と比較して、平方根平均二乗誤差（RMSE）および平均バイアス誤差（MBE）を有意に低減させることが示された。

意義および主張
著者らは、本研究の主要な意義は、電圧大きさのダイナミクスを明示的に考慮するIF表現を導出し、それによって簡略化された推定器に見られるバイアスを修正したことにあると主張している。本論文は、提案された定式化が非線形エネルギーと幾何学的解釈を統合するものであると述べている。これはあらゆる文脈ですべての既存の推定器に取って代わることを主張するものではなく、$\omega_{TV}$ が複素周波数の角成分の精密な局所実現を提供する一方で、$\omega_{TI}$ は大きさのダイナミクスが無視できる場合に信頼できる単純な近似として機能することを強調している。本研究は、三次元軌跡（非ゼロの捩れ）への拡張および他の周波数推定手法に対するより広範な検証のための理論的基礎を確立するものである。