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

기술 요약: Teager-Kaiser 에너지 연산자를 이용한 전력 계통의 순시 주파수 추정

문제 정의
순시 주파수(Instantaneous Frequency, IF)의 추정은 전력 계통의 비정상 현상을 특성화하는 데 매우 중요하다. 최근의 해석들은 전력 계통의 주파수를 복소 주파수 운동학 및 전압 궤적의 기하학적 구조와 연결 짓고 있으나, 기존의 추정 방식들은 흔히 Teager-Kaiser 에너지 연산자(Teager-Kaiser Energy Operator, TKEO)에 기반한 에너지 분리 알고리즘(Energy Separation Algorithm, ESA)에 의존한다. 이러한 전통적인 추정기들은 협대역 신호와 "느린 포락선(slow-envelope)" 조건(진폭 및 주파수 변조율이 반송파 주파수에 비해 작은 경우)을 가정한다. 전력 계통, 특히 과도 상태나 전압 크기가 크게 변하는 상황에서는 이러한 가정들이 빈번하게 위배된다. 결과적으로, 기존의 추정기들은 전압 크기 변화와 관련된 동적 항들을 무시하기 때문에 편향된 주파수 추정치를 생성할 수 있다.

방법론
저자들은 TKEO를 복소 신호로 확장(CTKEO)하고 동적 신호 항등식(dynamic-signal identity)을 적용하여 도출한 새로운 국소 IF 추정기를 제안한다. 방법론은 다음과 같이 진행된다:

1. 복소 TKEO 정식화: 표준 실수 값 TKEO를 대칭 이선형 형식(symmetric bilinear form)을 사용하여 복소 신호로 일반화한다. 복소 신호를 $\bar{x}$라 할 때, 연산자는 $\Psi_C(\bar{x}) = |\bar{x}'|^2 - \text{Re}(\bar{x}\bar{x}''^*)$로 정의된다. 이 정식화는 데카르트 성분에 대해 가법적(additive)이며 상수 위상 변화에 불변하므로, 협대역 가정을 필요로 하지 않는다.
2. 정확한 항등식 유도: 복소 신호를 크기 $|\bar{x}|$와 위상 $\phi$로 표현하고, 복소 주파수 $\bar{\eta} = \rho + j\omega$ (여기서 $\rho$는 크기 변화율, $\omega$는 각주파수)의 정의를 활용하여, 저자들은 다음과 같은 정확한 연속 시간 항등식을 유도한다:
   $$\Psi_C(\bar{x}) = |\bar{x}|^2 (2\omega^2 - \rho')$$
   이 식은 회전 에너지 성분과 크기 역학($\rho'$)과 관련된 병진(translation) 항을 명시적으로 분리한다.
3. 제안된 추정기:
   • $\omega_{TV}$ (시변형, Time-Varying): 항등식을 $\omega$에 대해 풀면, 크기 역학 보정 항($\rho'$ 및 $\rho^2$)을 유지하는 추정기를 얻는다. 이는 전압 크기가 급격히 변하는 조건에서 편향이 교정된 추정치를 제공한다.
   • $\omega_{TI}$ (시불변형, Time-Invariant): 크기가 일정하거나 느리게 변한다고 가정하여 보정 항을 무시하는 단순화된 근사치를 유도한다. 이는 $\omega_{TI} = \sqrt{\Psi_C(\bar{x}) / (2|\bar{x}|^2)}$가 되며, 계산이 더 간단하고 노이즈에 덜 민감하지만, 상당한 크기 과도 상태에서는 편향을 유발할 수 있다.
4. 기하학적 해석: 제안된 추정기는 복소 주파수 운동학과 일치하며, 비틀림(torsion)이 0이라고 가정할 때 $(v_\alpha, v_\beta)$ 평면에서의 전압 궤적 곡률과 관련된 기하학적 해석을 허용한다.

주요 기여

• 새로운 표현식: 본 논문은 포락선 곡률 항을 명시적으로 포함하는 CTKEO 기반의 IF 표현식을 도입한다. 이는 크기 역학을 무시함으로써 발생하는 편향을 교정한다.
• 편향 교정: 고전적인 ESA와 달리, 제안된 정식화($\omega_{TV}$)는 "느린 포랍선" 근사에 의존하지 않으므로, 연속 시간에서 매끄럽고 0이 아닌 2차원 전압 궤적에 대해 이론적으로 정확하다.
• 기하학적 일관성: 제안된 추정기는 복소 주파수 운동학과 일치하며, 명시적인 위상 언래핑(phase unwrapping) 없이도 주파수에 대한 기하학적 해석을 제공한다.
• 실용적 근사: 크기 변화가 느린 시나리오를 위해 계산 효율적인 근사치($\omega_{TI}$)를 제공하여, 정확도와 노이즈 민감도 사이의 절충안을 제시한다.

결과
제안된 방법들은 두 가지 사례 연구를 통해 검증되었다:

1. EMT 시뮬레이션 (IEEE 39-Bus 시스템): 노이즈가 존재하는 균형 상태에서의 3상 고장 시뮬레이션 결과, 제안된 추정기들($\omega_{TV}$ 및 $\omega_{TI}$)은 기하학적 참조 주파수($\omega_v$)와 밀접하게 일치하였다. 불균형 시나리오에서는 제안된 추정기들이 참조값에 비해 낮은 진동을 보였으며, 이는 불균형으로 인해 발생하는 곡률 왜곡에 덜 영향을 받음을 시사한다.
2. 현장 측정 (PV 발전소): 스페인의 PV 발전소에서 얻은 실제 전압 데이터에 대한 테스트 결과, 두 추정기 모두 기하학적 참조의 지배적인 경향을 따랐다. 높은 노이즈와 미분 민감도( $\omega_{TV}$는 2계 도함수를 요구함)가 발생하는 구간에서 차이가 나타났으나, 결과적으로 $\omega_{TI}$가 실제 운전 조건 하의 본래의 시변 추정치와 실질적으로 매우 유사함을 확인하였으며, 이는 크기 변화의 역학이 본질적으로 느리기 때문이다.

진폭 변조 신호에 대한 정량적 분석(Table I)에 따르면, 시변 추정기($\omega_{TV}$)는 시불변 추정기($\omega_{TI}$)에 비해, 특히 낮은 샘플링 주파수에서 제곱평균제곱오차(RMSE)와 평균 편향 오차(MBE)를 유의미하게 감소시켰다.

의의 및 주장
저자들은 본 연구의 주요 의의가 전압 크기 역학을 명시적으로 고려하는 IF 표현식을 유도함으로써, 단순화된 추정기들이 가진 편향을 교정했다는 점에 있다고 주장한다. 본 논문은 제안된 정식화가 비선형 에너지 해석과 기하학적 해석을 통합한다고 단언한다. 본 연구는 모든 맥락에서 기존 추정기를 대체하려는 것이 아니라, $\omega_{TV}$가 복소 주파수의 각 성분에 대한 정밀한 국소 실현을 제공하는 반면, $\omega_{TI}$는 크기 역학이 무시할 만한 수준일 때 신뢰할 수 있는 더 단순한 근사치로서 역할을 한다는 점을 강조한다. 본 연구는 향후 3차원 궤적(비제로 비틀림)으로의 확장 및 다른 주파수 추정 기법들과의 광범위한 검증을 위한 이론적 토대를 마련한다.