이질적 동조 발전기의 정확한 차원축소 모델링

본 논문은 전력망에서 동조 발전기 집합의 집합적 주파수 동역학을 정확히 모델링하고, 이를 효율적으로 차원축소하는 새로운 방법론을 제시한다. 먼저, 각 발전기의 선형화된 전송함수 \(g_i(s)\) 와 네트워크 라플라시안 \(L\) 을 이용해 전체 시스템의 전송 행렬 \(T(s)=(I+\operatorname{diag}\{g_i(s)\}L/s)^{-1}\operatorname{diag}\{g_i(s)\}\) 을 정의한다. 라플라시안의 두 번째 고유값인 대수 연결성 \(\lambda_2(L)\) 이 충분히 클 경우, 모든 버스의 주파수 응답이 동일한 형태 \(\hat w(s)=\hat g(s)\sum_{i=1}^n u_i(s)\) 로 수렴한다는 정리 1을 증명한다. 여기서 \(\hat g(s)=\big(\sum_{i=1}^n g_i^{-1}(s)\big)^{-1}\) 는 “동조 동역학”이라 불리며, 개별 발전기의 역함수 합의 역함수 형태로 나타난다. 동조 동역학 \(\hat g(s)\) 는 발전기 모델에 따라 다양한 형태를 갖는다. 1차 모델(동기 발전기 \(g_i(s)=1/(m_i s+d_i)\) 또는 인버터 \(g_i(s)=k_{P,i}/(\tau_{P,i}s+1)\))의 경우, \(\hat g(s)\)는 동일한 1차 형태로 간단히 합산된 관성·감쇠(또는 드롭·필터) 파라미터로 표현된다. 그러나 터빈 드롭 제어가 포함된 2차 모델 \(g_i(s)=1/(m_i s+d_i+r_i^{-1}/(\tau_i s+1))\) 에서는 각 터빈의 시간 상수 \(\tau_i\) 가 서로 다르면 \(\hat g(s)\)는 \(\sum_{i} r_i^{-1}\tau_i s+1\) 와 같이 서로 다른 극점을 가진 고차 전송함수가 된다. 이 경우 단순히 관성·감쇠만 합산한 저차 모델은 실제 동조 동역학을 제대로 재현하지 못한다. 이를 해결하기 위해 저자들은 주파수 가중 균형 절단(Frequency‑Weighted Balanced Truncation, FWBT) 기법을 도입한다. FWBT는 특정 주파수 구간(주로 저주파)에 높은 가중치를 부여해 해당 구간에서 시스템의 핵심 동적을 보존하면서 차원을 축소한다. 두 가지 축소 전략을 제시한다. 첫 번째 전략은 터빈 집합 \(\hat g_t(s)=\sum_{i=1}^n r_i^{-1}\tau_i s+1\) 에 직접 FWBT를 적용한다. k‑1 차의 축소 모델 \(\tilde g_{t,k-1}(s)\) 를 얻고, 이를 전체 모델에 결합해 \(\tilde g_{tb,k}(s)=1/(\hat m s+\hat d+\tilde g_{t,k-1}(s))\) 를 만든다. k=2이면 2차(단일 발전기) 모델이 되며, DC 이득 불일치가 존재하지만 적절한 가중치 선택으로 최소화한다. k=3이면 3차 모델(두 개의 1차 터빈 병렬)으로, 실제 시뮬레이션에서 원본 고차 \(\hat g(s)\)와 거의 동일한 응답을 제공한다. 두 번째 전략은 전체 고차 \(\hat g(s)\) 에 직접 FWBT를 적용하는 것으로, 차수 k 를 자유롭게 선택해 저주파 응답을 정확히 맞출 수 있다. 이 방법은 터빈 동역학뿐 아니라 관성·감쇠·드롭 전체를 동시에 축소하므로, 보다 유연한 모델링이 가능하다. 논문은 아이슬란드 전력망(중간 수준 연결성) 사례를 통해 제안된 이론을 검증한다. 시뮬레이션 결과, 개별 발전기의 고주파 진동은 존재하지만, 코이 주파수(전체 관성 가중 평균)와 \(\hat g(s)\) 의 응답은 매우 일치한다. 2차·3차 축소 모델을 적용했을 때, 원본 고차 모델과 거의 동일한 스텝 응답을 얻었으며, 특히 3차 모델은 거의 완벽한 근사를 제공한다. 결론적으로, 본 연구는 (1) 대수 연결성에 기반한 동조 동역학의 엄밀한 정의, (2) 이질 터빈 동역학을 고차 전송함수로 정확히 표현, (3) 주파수 가중 균형 절단을 이용한 단계적 차원축소 프레임워크를 제시함으로써 기존 관성·감쇠 합산 방식보다 높은 정확도와 해석 가능성을 제공한다는 점에서 전력 시스템 동역학 모델링에 중요한 기여를 한다.