전력계통 DAE 비반복 해법: 미분변환을 이용한 전력 흐름·동적 해석 혁신

본 논문은 전력계통의 미분대수방정식(DAE) 해석을 위한 새로운 비반복 알고리즘을 제안한다. 기존의 전력계통 시뮬레이션은 미분 방정식은 수치 적분(예: 암시적 사다리법)으로, 대수 방정식은 뉴턴‑라프슨 등 반복적 방법으로 각각 해결한다. 이러한 이중 반복 구조는 시간 스텝이 작아야 할 필요성과 비선형 네트워크 방정식의 수렴 어려움으로 인해 계산 부하가 크게 증가한다. 이를 극복하고자 저자는 ‘미분변환(Differential Transformation, DT)’이라는 수학적 도구를 도입한다. DT는 시간 함수 x(t)를 무한 파워 시리즈 ΣX(k)t^k 로 전개하고, 각 계수 X(k)를 직접 구할 수 있는 변환 규칙을 제공한다. 고차 도함수를 직접 계산할 필요가 없으며, 합성·곱·역함수 등에 대한 벡터·행렬 형태의 규칙을 제시해 전력계통 모델에 바로 적용할 수 있다. 논문은 먼저 전력계통 DAE를 상태 변수 x, 버스 전압 v, 전류 주입 i, 네트워크 어드미턴스 Ybus 로 구성된 상태‑공간 형태(2) 로 제시한다. DT를 적용하면 원래의 비선형 DAE는 계수들 X(k), V(k), F(k), I(k) 사이의 새로운 방정식 집합(4) 으로 변환된다. 여기서 중요한 점은 (4a)에서 X(k)만이 좌변에 나타나고, RHS는 이전 차수까지의 계수만 포함하므로 X(k)는 명시적 식(5) 로 바로 구할 수 있다는 것이다. 반면 (4b)에서는 V(k)가 좌·우변 모두에 등장해 직접 해석이 어려워 보이지만, 저자는 전류 주입 I(k)와 버스 전압 V(k) 사이에 선형 관계식(6)을 증명한다. 이 선형식은 행렬 A와 B가 이전 차수의 V(·)만을 포함하므로, V(k)를 식(7) 로 명시적으로 계산할 수 있다. 따라서 전체 알고리즘은 초기 조건 X(0), V(0)만 알면, 차례대로 X(1), V(1), X(2), V(2)… 를 순차적으로 구해 나가는 구조가 된다. 이 과정에서 반복적 수치 해석이 전혀 필요 없으며, 각 차수마다 행렬 연산과 컨볼루션만 수행한다. 논문은 이러한 이론적 기반을 바탕으로 실제 전력계통에 적용하기 위해 다음과 같은 절차를 제시한다. ① 전력계통 모델(발전기, 제어기, ZIP 부하 등)을 DT 규칙에 맞게 전개한다. ② 전류 주입식(9‑12)과 부하식(13‑20)을 벡터화된 변환 규칙(8) 로 변환한다. ③ 행렬 A·B 를 이전 차수 계수로부터 계산하고, 식(7) 로 V(k)를 구한다. ④ 식(5) 로 X(k)를 구하고, 이를 이용해 다음 차수로 진행한다. 이 과정을 원하는 시간 구간까지 반복하면, 시간에 대한 파워 시리즈 형태의 해를 얻는다. 파워 시리즈는 원하는 정밀도에 따라 차수를 선택할 수 있어, 실시간 요구사항에 맞게 조정 가능하다. 성능 검증을 위해 폴란드 2383버스 시스템에 적용하였다. 테스트에서는 기존의 암시적 사다리법(Trapezoidal)과 명시적 수정 오일러·뉴턴‑라프슨 혼합 스킴을 비교했다. 동일한 시뮬레이션 정확도(전압·주파수 오차 <10⁻⁴) 하에, DT 기반 비반복 알고리즘은 평균 5~10배 빠른 실행 시간을 기록했으며, 특히 큰 스텝(Δt=0.05 s)에서도 수렴 문제가 발생하지 않았다. 메모리 사용량도 비슷하거나 약간 낮았으며, 전압·전류 프로파일은 기존 방법과 거의 일치했다. 이러한 결과는 비선형 네트워크 방정식의 반복을 완전히 제거함으로써 계산 복잡도가 크게 감소했음을 보여준다. 마지막으로 논문은 향후 연구 방향으로 (1) 고차 비선형 요소(예: 비선형 제어기, 고장 모델)와의 통합, (2) 실시간 디지털 트윈 및 온라인 감시 시스템에의 적용, (3) 병렬/GPU 가속을 통한 대규모 시스템 확장성을 제시한다. 전반적으로 미분변환을 이용한 비반복 DAE 해법은 전력계통 시뮬레이션의 속도와 안정성을 동시에 개선할 수 있는 강력한 도구로 평가된다.