jacobhan.me
Power Systems Simulation

Power Flow · TSA · EMT — 전력망 시뮬레이션 세 계층 완전 비교

한 순간의 평형점, 수 초 동안의 흔들림, 그리고 파형 그 자체. 같은 전력망을 세 개의 시간 척도에서 푸는 세 가지 도구 이야기.

전력망을 컴퓨터로 모의하는 방법은 하나가 아니다. 무엇을 알고 싶은가에 따라 보는 시간의 단위가 달라지고, 시간의 단위가 달라지면 푸는 방정식의 종류 자체가 달라진다. 조류계산(Power Flow), 과도안정도 해석(Transient Stability Analysis, TSA), 전자기 과도 해석(Electromagnetic Transient, EMT)은 정밀도와 계산 비용을 맞바꾸는 세 개의 서로 다른 선택지다. 이 글은 세 도구가 각각 무엇을 가정하고 무엇을 버리는지, 그리고 왜 최근 들어 가장 무거운 EMT 해석이 변두리에서 표준의 자리로 밀려 올라오고 있는지를 정리한다.

0왜 도구가 셋이나 필요한가

전력 계통에서 벌어지는 현상은 시간 폭이 극단적으로 넓다. 송전선에 떨어진 낙뢰가 만드는 과전압은 100만 분의 1초(마이크로초) 단위에서 진행되고, 발전기 회전자가 서로 밀고 당기며 만드는 동요는 수 초에 걸쳐 나타나며, 수요와 공급을 맞추는 경제급전은 수 분에서 수 시간을 다룬다. 같은 전력망이라도 마이크로초의 세계와 시간 단위의 세계는 지배하는 물리 법칙의 우선순위가 다르다.

한쪽 끝을 모두 담을 수 있을 만큼 정밀한 모형은 그만큼 무겁다. 모든 도선의 인덕턴스와 커패시턴스를 일일이 미분방정식으로 따라가면 파형 하나하나가 보이지만, 그 정밀도로 수도권 전체 계통의 한나절을 모의하는 것은 현실적인 계산 시간 안에 끝나지 않는다. 반대로 충분히 가벼운 모형은 빠르고 큰 계통을 다룰 수 있지만, 그 대가로 빠른 현상을 통째로 지워 버린다.

그래서 도구가 셋이다. 정밀도(detail)와 규모·속도(scale·speed)는 한쪽을 얻으면 다른 쪽을 내주는 교환 관계에 있고, 세 도구는 이 교환선 위의 서로 다른 세 지점을 차지한다. 조류계산은 시간 자체를 버려 가장 빠르고, EMT는 파형까지 붙들어 가장 무거우며, 과도안정도 해석은 그 사이에 선다.

쉽게 풀어보면

지도를 떠올리면 편하다. 전국 고속도로망의 큰 흐름만 보고 싶으면 축척이 큰 전국 지도 한 장이면 된다. 한 동네의 골목 구조를 알고 싶으면 동네 지도를 펴야 하고, 특정 교차로에서 신호가 어떻게 바뀌는지를 초 단위로 보려면 그 교차로만 확대한 그림이 필요하다. 세 지도는 모두 같은 나라를 그리지만 담는 시간과 공간의 단위가 다르다.

조류계산은 전국 지도, 과도안정도 해석은 동네 지도, EMT는 교차로 확대도에 가깝다. 한 장으로 전부를 대신할 수 없기 때문에, 알고 싶은 것이 무엇이냐에 따라 펼치는 지도를 바꾸는 것이다.

1시간 척도의 스펙트럼

세 도구의 자리를 가장 직관적으로 보여 주는 것은 시간축이다. 가로축을 로그 척도로 펼쳐 100만 분의 1초부터 몇 시간까지를 한 줄에 놓으면, 각 물리 현상이 어느 구간에서 일어나는지가 한눈에 들어온다.

전력계통 현상의 시간 척도와 해석 도구마이크로초의 낙뢰부터 시간 단위 경제급전까지의 시간축 위에 EMT·과도안정도·조류계산의 적용 범위를 표시1µs0.1ms10ms1s100s2.8h물리 현상낙뢰개폐·스위칭고조파동기 이하 진동전기기계 동요주파수·장기시간영역 해석 도구의 적용 범위EMT (전자기 과도)TSA (과도안정도)Power Flow (조류계산): 시간을 적분하지 않는다 — 한 순간의 정상상태 평형점만 계산
전력 계통 현상의 시간 척도와 해석 도구의 적용 범위. 위쪽 회색 막대는 물리 현상이 일어나는 시간 구간을, 아래쪽 색 막대는 EMT과도안정도 해석이 다루는 범위를 나타낸다. 조류계산은 시간을 적분하지 않으므로 시간축 위에 올라가지 않고 별도 박스로 표시했다.

낙뢰와 개폐 서지, 고조파처럼 마이크로초에서 밀리초 사이의 빠른 현상은 EMT의 영역이다. 발전기 회전자가 서로 동기를 유지하느냐를 다투는 전기기계적 동요는 대략 0.1초에서 수십 초 구간에서 일어나고, 이는 과도안정도 해석의 본령이다. 두 도구의 적용 범위는 1초 안팎에서 겹친다. 인버터 기반 자원이 늘면서 새로 문제가 된 동기 이하 진동(sub-synchronous oscillation)은 공교롭게도 바로 이 겹치는 구간에 걸쳐 있는데, 이 점이 뒤에서 다룰 핵심이다.

조류계산이 이 그림에서 시간축 아래 별도 박스로 빠져 있는 이유는 분명하다. 조류계산은 시간을 따라가지 않는다. 어느 한 순간 계통이 정착해 있을 평형 상태를 단번에 푸는 일이라, 처음도 끝도 없이 "지금 이대로라면 어떻게 균형이 잡히는가"라는 한 장면만 계산한다.

쉽게 풀어보면

사진과 동영상의 차이를 생각하면 된다. 조류계산은 셔터를 한 번 눌러 찍은 정지 사진이다. 시간이 흐르지 않으니 "그 순간 모든 것이 어떻게 놓여 있었나"만 담긴다.

과도안정도 해석과 EMT는 동영상이다. 다만 같은 동영상이라도 초당 프레임 수가 다르다. 과도안정도 해석은 몇 밀리초마다 한 장면을 찍어 수 초짜리 영상을 만들고, EMT는 100만 분의 1초마다 한 장면을 찍어 아주 짧지만 훨씬 매끄러운 영상을 만든다. 프레임이 촘촘할수록 빠른 움직임이 또렷이 보이지만, 같은 길이를 찍으려면 저장 공간과 시간이 그만큼 더 든다.

2Power Flow — 정지한 평형

조류계산 조류계산은 전력 계통 해석의 가장 기본이 되는 정상상태(steady-state) 계산이다. 발전과 부하가 어딘가에서 균형을 이루고 있다고 보고, 그 균형점에서 각 모선(bus, 송전망의 접속점)의 전압 크기와 위상, 그리고 각 선로에 흐르는 유효전력과 무효전력을 구한다. 시간 미분이 등장하지 않는 대수방정식(algebraic equation)만으로 이루어져 있어, 세 도구 중 압도적으로 가볍다는 표현 대신 — 가장 적은 계산으로 답을 낸다.

조류계산은 교류 파형을 매 순간 추적하는 대신, 정현파 한 주기를 하나의 복소수(페이저, phasor)로 압축해 다룬다. 크기와 위상만 남기고 시간 변화는 지워 버리는 것이다. 또한 삼상이 완벽히 평형을 이룬다고 가정해 한 상만 대표로 계산하며(정상분, positive-sequence), 계통 전체가 하나의 주파수에서 안정적으로 진동한다고 전제한다.

모선은 역할에 따라 세 종류로 나뉜다. 전압과 위상의 기준이 되는 슬랙 모선(slack bus), 발전기가 붙어 전압 크기와 유효전력이 정해지는 PV 모선, 부하가 붙어 유효·무효전력이 정해지는 PQ 모선이다. 미지수가 비선형으로 얽혀 있어 한 번에 풀리지 않으므로, 추정값에서 출발해 반복적으로 답에 다가가는 수치 해법을 쓴다. 가우스–자이델법(Gauss–Seidel)이 가장 단순하고, 뉴턴–랩슨법(Newton–Raphson)은 답 근처에서 오차가 제곱으로 줄어드는 빠른 수렴을 보여 대규모 계통에서 표준으로 쓰이며, 고속 분리해법(Fast Decoupled)은 전력망의 물리적 특성을 이용해 계산을 더 줄인다.

이 해법들이 공통으로 만족시키려는 것은 각 모선에서 들어오고 나가는 전력이 같아야 한다는 전력 평형 조건이다. 이를 식으로 쓰면 다음과 같다.

전력 조류 방정식 — 정상상태 전력 평형
유효전력 주입
Pi = |Vi| k |Vk| ( Gik cos θik + Bik sin θik )
무효전력 주입
Qi = |Vi| k |Vk| ( Gik sin θik Bik cos θik )
뉴턴–랩슨 보정 — 추정값을 반복 갱신
[ΔPΔQ] = J [ΔθΔ|V|]
P·Q는 모선 i의 유효·무효전력 주입, |V|는 전압 크기, θik = θiθk(두 모선의 위상차), Gik + jBik는 어드미턴스 행렬의 성분이다. 합(∑)은 모선 i에 이어진 모든 모선 k에 대해 더한다. J는 자코비안(편미분으로 이루어진 행렬)으로, 추정값을 이 식으로 보정해 전력 평형에 수렴시킨다. 시간 미분이 전혀 없다는 점이 정상상태 계산임을 드러낸다.

조류계산이 답하는 질문은 이런 것들이다. 이 송전선에 과부하가 걸리지는 않는가, 어느 모선의 전압이 규정 범위를 벗어나는가, 발전기 한 대나 선로 하나가 갑자기 빠져도 나머지가 견디는가(N-1 상정고장 해석). 계통 계획, 경제급전, 대량의 시나리오를 빠르게 걸러 내는 일이 모두 조류계산 위에서 이루어진다.

쉽게 풀어보면

여러 개의 물탱크가 굵기가 제각각인 파이프로 연결되어 있다고 하자. 어떤 탱크에는 펌프가 물을 밀어 넣고(발전), 어떤 탱크에서는 물을 퍼낸다(부하). 충분한 시간이 지나면 물은 더 이상 출렁이지 않고 각 탱크의 수위와 파이프의 물 흐름이 일정한 값에 정착한다.

조류계산은 바로 그 "출렁임이 멎은 뒤의 최종 수위와 흐름"만 계산한다. 물이 어떻게 출렁이다 잦아들었는지는 묻지 않는다. 그래서 빠르지만, 펌프가 갑자기 꺼졌을 때 물탱크가 잠시 어떻게 요동치는지는 알려 주지 못한다. 그 요동을 보려면 다음 도구가 필요하다.

3TSA — 수 초 동안의 흔들림

과도안정도 과도안정도 해석은 큰 사고가 난 직후 계통이 새로운 안정 상태로 무사히 옮겨 가는지, 아니면 발전기들이 서로 동기를 잃고 흩어지는지를 따진다. 송전선 단락 같은 외란이 들어온 뒤 수 초 동안 발전기 회전자의 각도와 계통 주파수가 어떻게 움직이는지를 시간에 따라 추적한다.

이 해석의 뼈대는 미분-대수방정식(Differential-Algebraic Equation, DAE)이다. 발전기의 회전 운동과 제어기의 동특성은 시간 미분을 포함하는 미분방정식으로, 송전망 자체는 매 시각의 전력 균형을 맞추는 대수방정식으로 표현된다. 회전자의 운동을 지배하는 핵심은 동요방정식(swing equation)으로, 회전자에 들어오는 기계적 입력과 전기적으로 빠져나가는 출력의 차이가 회전자를 가속하거나 감속한다는 관계다.

과도안정도 해석의 지배 방정식
미분 방정식 — 발전기·제어기 동특성
dx/dt = f(x, y)
대수 방정식 — 송전망 전력 균형
0 = g(x, y)
동요방정식 — 회전자 운동
2Hωs d2δdt2 = Pm Pe
전기 출력 (단순 모델)
Pe = |E′| |V|X sin δ
x는 상태변수(회전자 각 δ·각속도, 제어기 내부 상태), y는 대수변수(모선 전압·위상 등)이다. H는 관성정수, ωs는 동기 각속도, Pm은 기계 입력, Pe는 전기 출력이다. 단순 모델에서 E′는 발전기 내부 기전력, V는 단자전압, X는 그 사이 리액턴스이며, 회전자 각 δ가 벌어질수록 sin δ를 통해 전기 출력이 달라진다. 미분식과 대수식이 한데 묶여 푸는 점이 미분-대수방정식의 핵심이다.

여기에 발전기 단자전압을 조절하는 자동전압조정기(Automatic Voltage Regulator, AVR)와 여자 시스템, 회전 속도를 다스리는 조속기(governor)의 동특성이 함께 묶인다. 다만 송전망의 전압과 전류는 여전히 실효값(Root Mean Square, RMS)을 담은 페이저로, 정상분으로, 기본 주파수만으로 다룬다. 회전자처럼 느린 것은 미분으로 따라가되, 송전망처럼 빠른 것은 매 시각 즉시 균형이 잡힌다고 보는 것이다. 시간 간격은 대략 1~10밀리초, 모의 구간은 수 초에서 수십 초 정도다.

과도안정도가 드러내는 것: 회전자 각의 운명고장 제거 후 발전기 회전자 각이 감쇠하며 회복하면 안정, 계속 증가해 동기에서 이탈하면 불안정012390°180°270°동기 이탈 임계(예시)안정: 감쇠 후 회복불안정: 동기 이탈고장 제거 후 시간 (초)회전자 각 δ (도)
고장 제거 후 발전기 회전자 각 δ의 거동. 안정한 경우(녹색)에는 몇 차례 흔들린 뒤 새로운 평형각으로 감쇠해 정착하지만, 불안정한 경우(붉은색)에는 각이 계속 벌어져 동기 이탈에 이른다.

이 그림이 과도안정도 해석의 본질을 압축한다. 같은 사고라도 계통 조건에 따라 회전자 각은 출렁이다 가라앉기도 하고, 임계선을 넘어 끝없이 벌어지기도 한다. 후자가 되면 발전기는 동기를 잃고, 보호계전기가 이를 분리하면서 정전으로 번질 수 있다. 그 갈림길을 미리 계산해 두는 것이 이 해석의 목적이다. 대표적인 도구로는 PSS/E, PSLF, PowerFactory의 RMS 모드, TSAT 등이 있다.

쉽게 풀어보면

손가락 끝에 긴 막대를 세워 균형을 잡고 있다고 하자. 누가 막대를 툭 쳤을 때, 약하게 쳤다면 막대가 몇 번 휘청이다 다시 똑바로 선다. 그러나 세게 쳐서 어떤 각도를 넘어 기울면, 그때부터는 무엇을 해도 막대가 점점 더 기울어 결국 쓰러진다. 그 "넘으면 돌이킬 수 없는 각도"가 안정과 불안정을 가르는 임계점이다.

과도안정도 해석은 사고라는 충격을 받은 발전기들이 휘청인 뒤 제자리로 돌아오는지, 아니면 임계각을 넘어 쓰러지는지를 수 초 동안 지켜보는 일이다. 다만 막대 표면의 미세한 진동까지 보는 것은 아니다. 막대가 전체적으로 기우는 큰 움직임에 집중한다.

4EMT — 파형 그 자체

전자기 과도 전자기 과도 해석은 페이저라는 압축을 풀고 전압과 전류의 순시값(instantaneous value), 즉 시간에 따라 출렁이는 정현파 그 자체를 직접 계산한다. 한 주기를 복소수 하나로 뭉뚱그리지 않고, 파형의 모든 점을 따라간다.

이것이 가능하려면 계통의 모든 인덕턴스와 커패시턴스를 빠짐없이 미분방정식으로 모형화해야 한다. 도선과 변압기의 인덕턴스에 걸리는 전압은 전류의 시간 변화율에 비례하고, 커패시턴스에 흐르는 전류는 전압의 시간 변화율에 비례한다. 바로 이 미분 관계 때문에 전류와 전압은 한순간에 도약하지 못하고 연속적으로만 변할 수 있는데, 조류계산이나 과도안정도 해석이 사고를 순간적인 계단으로 처리하는 것과 달리 EMT는 그 연속적인 변화 과정 전체를 그려 낸다.

EMT는 삼상을 각각 따로 다루므로 불평형 상태를 그대로 표현하고, 고조파와 차단기 개폐에 따른 과도 현상, 직류 성분이 얹힌 비대칭 고장 전류까지 담아낸다. 계산의 토대는 1969년 도멜(Hermann Dommel)이 제시한 방법으로, 사다리꼴 적분(trapezoidal integration)으로 미분방정식을 차분화하고 각 소자를 등가 회로(companion model)로 바꾼 뒤 절점 해석(nodal analysis)으로 매 시각의 전압을 푸는 방식이다. 송전선은 신호가 선로를 타고 전파되는 시간까지 반영하는 분포정수 모형이나 주파수 의존 모형으로 다룬다.

도멜 이산화 — 컴패니언 모델과 절점 방정식
인덕터의 사다리꼴 적분
iL(t) = Δt2L vL(t) + Ihist
이력 전류원 — 과거 값으로 결정
Ihist = iL(t−Δt) + Δt2L vL(t−Δt)
절점 방정식 — 매 시간스텝마다 푼다
G v(t) = i(t)
Δt는 시간 간격, L은 인덕턴스이다. 인덕터는 한 개의 저항과 한 개의 전류원으로 이루어진 등가 회로(컴패니언 모델)로 바뀌는데, 그 전류원 Ihist는 바로 직전 시각의 전압·전류로 정해진다. 커패시터도 같은 방식으로 바뀐다. G는 이 등가 저항들을 모은 컨덕턴스 행렬, v(t)는 절점 전압 벡터, i(t)는 이력 전류원을 포함한 주입 전류 벡터다. 결국 매 시간스텝에서 이 선형 연립식 한 번을 풀어 순시 전압을 얻고, 이를 수없이 반복해 파형을 그린다.

대가는 계산량이다. 시간 간격이 대략 1~50마이크로초로 과도안정도 해석보다 천 배가량 촘촘하다. 그만큼 같은 길이를 모의해도 단계 수가 폭증하고, 다룰 수 있는 구간은 보통 수 밀리초에서 길어야 수 초에 그친다. 대표 도구로는 PSCAD/EMTDC, EMTP, ATP가 있고, 실시간으로 돌려 실제 제어기 하드웨어와 연결해 시험하는 RTDS, OPAL-RT 같은 실시간 모의 장비도 이 범주에 속한다.

쉽게 풀어보면

앞 절에서 과도안정도 해석은 막대가 전체적으로 기우는 큰 움직임을 본다고 했다. EMT는 거기에 더해 막대 표면을 타고 흐르는 미세한 떨림 하나하나까지 초고속 카메라로 잡아내는 일에 가깝다.

음악으로 비유하면, 조류계산이 "이 곡은 라장조다"라고 조성만 말하고 과도안정도 해석이 멜로디의 큰 흐름을 그린다면, EMT는 공기를 진동시키는 음파의 파형 그 자체를 기록한다. 가장 충실하지만, 단 몇 초의 소리를 그렇게 기록하려 해도 엄청난 양의 데이터가 필요하다. 충실함과 무게는 같이 간다.

5한눈에 보는 핵심 비교

세 도구의 차이는 같은 사고를 각각 어떻게 그려 내는지를 나란히 놓고 보면 가장 선명하다. 송전선에 단락 고장이 한 번 일어나고 보호장치가 이를 제거하는 동일한 사건을, 세 해석은 전혀 다른 그림으로 내놓는다.

같은 사고를 세 해석이 보는 방식한 번의 고장에 대해 조류계산은 두 평형점, 과도안정도는 RMS 크기의 포락선, EMT는 순시 파형을 보여준다Power Flow두 평형점 (점)TSARMS 포락선 (선)EMT순시 파형사고 전사고 후|V|, 사이 구간은 풀지 않음|V| 포락선, 창=0–1초v(t) 순시값, 창=0–100ms (확대)점(상태) → 포락선(크기) → 파형(순시값): 같은 사건, 점점 더 정밀한 시선
하나의 고장을 세 해석이 보는 방식. 조류계산은 고장 전과 후의 두 평형 전압을 점으로 찍을 뿐 그 사이를 풀지 않는다. 과도안정도 해석은 전압 크기(RMS)가 떨어졌다 회복하는 포락선을 그린다. EMT는 같은 구간을 확대해 순시 전압 파형이 직류 성분을 머금고 출렁이는 모습까지 보여 준다.

조류계산은 사고 직전과 직후의 두 정상상태를 따로 풀어 점으로 비교할 뿐, 그 사이에 무슨 일이 있었는지는 묻지 않는다. 과도안정도 해석은 전압의 크기가 어떻게 가라앉았다 돌아오는지 그 외피(포락선)를 그리지만, 파형 한 주기 안의 떨림은 매끈하게 뭉개 버린다. EMT만이 같은 고장을 100분의 1초 단위로 확대해, 전류와 전압이 한순간에 뛰지 못하고 직류 성분을 끌며 연속적으로 변하는 실제 파형을 드러낸다.

항목Power Flow
조류계산
TSA
과도안정도
EMT
전자기 과도
다루는 양정상상태 페이저 (전압 크기·위상)RMS 페이저의 시간 거동순시값 파형 그 자체
방정식대수방정식만미분-대수방정식 (DAE)모든 인덕턴스·커패시턴스의 미분방정식
시간없음 (한 순간)1~10 ms 간격 · 수 초~수십 초1~50 µs 간격 · ms~수 초
상(phase) 표현정상분 (평형 가정)정상분 (평형 가정)삼상 개별 (불평형 가능)
주파수단일 기본 주파수기본 주파수 중심광대역 (고조파 포함)
사고 표현전·후 정적 비교순간 계단 변화연속 파형 변화
대표 도구뉴턴–랩슨 솔버, PSS/E 조류PSS/E, TSAT, PowerFactory RMSPSCAD/EMTDC, EMTP, RTDS
전형적 용도계통계획, 경제급전, N-1 선별과도·주파수 안정도낙뢰·개폐 과도, HVDC·인버터 상세, 고조파
상대 계산비용가장 낮음중간가장 높음
쉽게 풀어보면

같은 교통사고를 세 사람이 다르게 기록한다고 해보자. 첫 번째 사람은 사고 직전과 직후의 차량 위치만 사진 두 장으로 남긴다(조류계산). 두 번째 사람은 사고 순간부터 몇 초 동안 차들이 어떻게 밀리고 멈췄는지 그 궤적을 스케치한다(과도안정도). 세 번째 사람은 충돌 순간을 초고속 카메라로 찍어 범퍼가 찌그러지는 0.01초의 변형까지 담는다(EMT).

세 기록은 모두 같은 사고를 다루지만, 알려 주는 것이 다르다. 보험금 정산에는 첫 번째면 충분하고, 사고 경위 분석에는 두 번째가 필요하며, 충돌 안전성 설계에는 세 번째가 있어야 한다. 표의 모든 항목은 결국 "얼마나 확대해서 보느냐"라는 한 가지 선택의 결과다.

6어떤 도구를 언제 — 그리고 왜 EMT가 떠오르는가

오랫동안 세 도구의 역할 분담은 명확했다. 일상적인 계통 계획과 운영은 조류계산이, 큰 외란에 대한 안정도 판정은 과도안정도 해석이 맡고, EMT는 낙뢰와 개폐 과도, 고압직류송전(High Voltage Direct Current, HVDC)이나 유연송전시스템(Flexible AC Transmission System, FACTS)처럼 특수한 전력전자 설비를 다룰 때만 꺼내는 전문 도구였다. EMT는 무겁고 느렸기에, 꼭 필요한 곳에만 제한적으로 썼다.

이 구도를 흔든 것이 인버터 기반 자원(Inverter-Based Resource, IBR)의 급증이다. 태양광, 풍력, 배터리는 회전하는 동기기가 아니라 전력전자 인버터를 거쳐 계통에 연결된다. 인버터는 물리적 관성이 거의 없고, 그 거동은 순시 전압을 읽어 한 주기보다 훨씬 짧은 시간 안에 반응하는 제어 소프트웨어가 결정한다. 무엇을 출력할지가 회전체의 관성이 아니라 코드로 정해진다는 뜻이다.

페이저 모형이 놓치는 것

과도안정도 해석을 포함한 페이저 기반 모형은 "송전망은 매 순간 즉시 균형을 잡고, 느린 동특성만 시간에 따라 추적하면 된다"는 시간 척도 분리(time-scale separation)를 전제로 한다. 발전이 거대한 회전 기계로 이루어지던 시절에는 이 전제가 잘 들어맞았다. 그러나 인버터의 제어 루프가 한 주기 안쪽의 시간 척도에서 동작하면, 빠른 것과 느린 것을 분리한다는 전제 자체가 무너진다. 페이저 모형은 사고를 전압 크기의 순간적인 계단 변화로 처리하지만, 정작 인버터가 실제로 읽는 것은 그 순간순간의 파형 위 전압이다. 모형이 보는 것과 인버터가 보는 것이 어긋나는 것이다.

그 어긋남은 이미 현실의 진동으로 나타났다. 텍사스 계통운영기관(Electric Reliability Council of Texas, ERCOT) 관할 계통에서는 약 4헤르츠 대역의 전압 진동이, 캐나다 하이드로원(Hydro One) 계통에서는 약 3.5헤르츠 대역의 유효·무효전력 진동이 관측되었다. 과거 직렬 보상 송전선에서 주로 보이던 동기 이하 제어 상호작용(sub-synchronous control interaction)이, 이제는 직렬 보상이 없는 약한 계통에서도 인버터 제어기 사이의 상호작용으로 나타나기 시작한 것이다. 이런 현상은 RMS 기반 모형으로는 예측되지 않는 경우가 많다.

규제가 EMT를 요구하기 시작했다

북미전력신뢰도공사(North American Electric Reliability Corporation, NERC)는 2025년 5월, 대규모 인버터 자원이 예기치 않게 동시 탈락하는 사고들이 잇따르고 기존 RMS 모형이 이를 사전에 잡아내지 못한 점을 들어 높은 수준의 신뢰도 경보를 발령했다. 여기에는 발전사업자와 송전계획 주체들이 자신들의 EMT 모형화 절차를 보고하도록 하는 요구가 포함되었다. EMT 검토가 권장 사항에서 점차 의무에 가까운 절차로 옮겨 가고 있음을 보여 주는 신호다. 인버터 자원의 계통 연계를 다루는 국제 표준인 IEEE 2800-2022와 분산자원 연계 표준인 IEEE 1547-2018도 이러한 흐름 위에 놓여 있다.

그렇다고 페이저 모형이 폐기되는 것은 아니다. 수천 개의 상정고장을 빠르게 걸러 내고, 설비 용량을 계획하고, 넓은 범위의 계통 연계를 개략적으로 검토하는 일에서 조류계산과 과도안정도 해석은 여전히 가장 효율적인 도구다. 달라진 것은 경계다. 인버터 비중이 높은 계통, 고압직류송전 연계, 빠른 전력전자 제어와 고조파, 불평형 고장, 한 주기 안쪽의 상호작용이 문제가 되는 지점에서는 이제 EMT가 선택이 아니라 기본이 되어 가고 있다.

쉽게 풀어보면

과거의 발전기는 무겁고 둔한 대형 트럭과 같았다. 도로 사정이 조금 바뀌어도 관성 때문에 천천히 반응하니, 큰 흐름만 봐도 거동을 예측할 수 있었다. 반면 인버터는 운전자의 손끝에 즉각 반응하는 가벼운 경주용 차에 가깝다. 노면의 미세한 요철 하나에도 즉시 핸들을 트니, 큰 흐름만 보는 지도로는 그 빠른 반응을 따라갈 수 없다.

도로에 둔한 트럭만 다닐 때는 동네 지도(과도안정도 해석)로 충분했지만, 민첩한 경주차가 많아진 지금은 그 차들이 노면과 어떻게 실시간으로 주고받는지를 보는 확대도(EMT)가 필요해졌다. 규제 기관이 EMT 검토를 요구하기 시작한 것은, 도로 위 차종이 바뀌었으니 들여다보는 지도도 바꾸라는 요구인 셈이다.

7세 도구는 경쟁자가 아니라 포개진 층이다

지금까지 셋을 나란히 비교했지만, 더 정확한 그림은 서로 포개진 동심원이다. EMT가 모형화하는 물리는 과도안정도 해석이 다루는 것을 품고, 과도안정도 해석은 다시 조류계산을 품는다. 안쪽으로 갈수록 가정이 늘어 가벼워지고, 바깥으로 갈수록 가정을 벗고 무거워진다.

각 해석이 풀고 포착하는 동특성의 범위정상상태 대수해석을 안쪽에 두고 전기기계 동특성, 전자기 동특성 순으로 포함 범위가 넓어지는 동심 구조EMT · 전자기 과도순시값 · 모든 L·C 미분방정식 · 3상 · 고조파 · 스위칭TSA · 과도안정도회전자·제어 미분 + 매 시점 망(대수) · RMS 위상Power Flow · 조류계산대수방정식 · 정상상태안쪽일수록 단순·고속, 바깥쪽일수록 상세·저속 — 바깥 해석은 안쪽이 보는 현상을 모두 포함한다
세 해석의 모형화 범위는 경쟁 관계가 아니라 포함 관계에 가깝다. 바깥의 EMT는 순시 파형과 모든 소자의 동특성을, 가운데 과도안정도는 회전자·제어기 동특성과 RMS 거동을, 안쪽 조류계산은 정상상태 평형만을 다룬다.

이 포함 관계가 실무에서 가장 중요한 절충안을 낳는다. 거대한 계통 전체를 EMT로 푸는 것은 여전히 비현실적이지만, 정작 빠른 현상이 문제가 되는 곳은 인버터 군집이나 고압직류송전 연계점처럼 계통의 일부다. 그래서 관심 영역만 EMT로 상세히 모형화하고, 나머지 광역 계통은 RMS로 가볍게 두어 경계에서 둘을 주고받게 하는 연성 해석(co-simulation)이 쓰인다. 필요한 곳에만 EMT의 정밀도를 사고, 나머지는 RMS의 속도를 누리는 전략이다.

그 중간 어디쯤에는 동기 회전 좌표계(dq 좌표계) 기반의 해석도 있다. 송전망의 동특성을 미분방정식으로 살려 두되 파형의 모든 점을 좇는 대신 회전 좌표계 위에서 다루어, 순시값 EMT보다는 가볍고 RMS보다는 더 많은 빠른 동특성을 담는다. 시간 척도 분리가 깨지지만 완전한 파형 해석까지는 과한 경우에 적합하다.

여기에 계산 자체를 빠르게 만드는 흐름이 더해진다. 다수의 연산을 병렬로 처리하는 그래픽처리장치(Graphics Processing Unit, GPU) 가속과, 모형을 실제 시간과 같은 속도로 돌려 실물 제어기 하드웨어를 연결해 시험하는 실시간 모의(하드웨어 인 더 루프, Hardware-in-the-Loop)가 그것이다. 이런 기술은 EMT 수준의 정밀도를 점점 더 큰 계통과 더 긴 구간으로 확장하고 있다. 정밀도와 규모의 교환선 자체를 바깥으로 밀어내는 셈이다.

쉽게 풀어보면

디지털 지도 앱을 떠올리면 된다. 전국을 한눈에 보다가 한 도시로, 다시 한 골목으로 손가락을 벌려 확대해 들어가도 그것은 모두 같은 하나의 지도다. 축척이 바뀔 뿐, 큰 그림과 골목이 서로 다른 지도인 것은 아니다. 세 해석도 같은 전력망을 서로 다른 배율로 본 것이다.

연성 해석은 이 비유를 한 걸음 더 민다. 지도 전체는 멀리서 보되, 사고가 잦은 교차로 한 곳만 확대해 끼워 넣는 것과 같다. 전부를 확대하면 너무 무거우니, 꼭 봐야 할 곳만 골라 정밀하게 들여다보는 것이다.

8맺으며

조류계산과 과도안정도 해석, EMT의 차이는 결국 "시간을 얼마나 잘게 쪼개 보느냐"라는 한 가지 선택으로 수렴한다. 시간을 통째로 버리면 가장 빠르고 큰 계통을 다룰 수 있고, 마이크로초까지 쪼개면 파형의 진실에 가장 가까워지지만 다룰 수 있는 범위가 좁아진다. 세 도구는 이 교환 관계 위의 세 좌표이며, 어느 하나가 다른 것을 대체하지 않는다.

다만 그 좌표 위에서 무게중심은 분명히 움직이고 있다. 회전 기계가 지배하던 계통에서는 페이저 모형만으로 대부분의 질문에 답할 수 있었지만, 코드로 움직이는 인버터가 계통을 채워 가면서 파형 단위의 진실을 요구하는 질문이 늘고 있다. 가장 무거웠던 EMT가 변두리에서 표준의 자리로 올라오는 흐름, 그리고 그 무게를 연성 해석과 가속 기술로 덜어 내려는 노력은 같은 변화의 양면이다.