윤치명
(Chi-Myeong Yun)
1iD
안태풍
(Tae-Pung An)
2
이지혜
(Jihye Lee)
2
민명환
(Myung-Hwan Min)
†iD
-
(Dept. of Transportation Engineering, Korea University of Science and Technology, Korea)
-
(ENTEC Electric & Electronic Co., Ltd., Korea.)
Copyright © The Korean Institute of Electrical Engineers(KIEE)
Key words
Regenerative energy, DC railway, Planning, Train Performance Simulation, Load flow
1. 서 론
전 세계적으로 기후변화에 대응하기 위해 탄소 배출 감소를 목적으로 많은 연구와 노력이 이루어지고 있다. 탄소 배출 감소를 위해 효율적인 에너지 소모가
필연적이며, 철도시스템 분야에서도 에너지 효율을 높이기 위해 회생제동에 의해 발생되는 에너지 활용을 위한 많은 연구가 수행되고 있다(1,2). 도시철도의 경우 역 간의 거리가 짧고, 차량의 역행과 제동이 빈번하여 회생제동에 의해 높은 에너지가 발생되는 특징이 있다. 차량 제동 시 짧은
시간 동안 높은 에너지가 발생되어 적절한 회생 에너지 활용이 필요하다(3).
도시철도에서 에너지 효율을 높이는 방법으로 차량 운행 시격 최적화, 기계적 제동과 회생제동을 혼합한 최적 제동을 통한 에너지 효율 증대 등 다양한
연구가 진행되고 있고, 에너지저장장치(ESS, Energy Storage System), 회생 인버터 그리고 회생 저항기를 이용해 회생 에너지를 활용하는
연구가 진행되고 있다. 회생 에너지 활용 설비 각각에 대한 효율 향상 또는 설치 위치 최적화를 통해 회생 에너지를 최대한 활용하여 운영요금 절감 등을
목표로 하고 있다(4-6). 하지만, 회생 에너지 활용설비 설계 시 구체적인 설계방안 및 인접 변전소 설치 시 영향 등이 고려되지 않고 있다.
전체적인 계통의 조류해석을 통한 회생 에너지 활용률을 고려하여 회생 에너지 활용설비들을 설치하여야 한다. 본 논문에서는 TPS(Train Performance
Simulation)를 이용하여 차량 한 대의 운행 프로파일과 필요 전력량을 산출한 후 DCPS(DC Power Simulation)를 이용해 운행
시격에 따른 조류해석을 수행한다. 조류해석 후 최저 소거법 방식을 이용하여 최저 회생 에너지 발생 변전소를 순차적으로 제거하여 설계하는 방식에 대해
소개한다. 제안하는 방법론의 진위성을 검증하기 위해 실제 운영 중인 도시철도 A 노선 데이터를 활용하여 이를 검증하였다.
2. TPS를 이용한 직류전기철도의 조류해석 방안
2.1 Train Performance Simulation
회생 에너지 활용설비 설치 설계 방법 소개에 앞서 직류 급전계통조류해석 방안에 대한 분석이 필요하다. 도시철도 계통 특성상 같은 노선 내에서도 선로의
구배나 커브에 의해 발생되는 회생 에너지 발생량이 상이하므로 급전 계통 해석을 통한 분석이 필요하다. 직류 급전 계통 조류해석은 두 단계로 나누어져
수행된다. Train Performance Simulation(TPS)를 이용해 차량의 속도 프로파일과 차량 운행에 필요한 구간별 필요 전력량 및
회생 발생량을 산출하고, DC Power Simulation(DCPS)를 이용해 설정된 운행 시격에 따라 전체 계통에 대한 조류해석 결과와 전압 및
회생 전력 발생량을 도출할 수 있다. 전체 시스템 급전계통 해석에 대한 순서도는 그림 1과 같다(7).
그림 1 직류 급전계통 조류해석 방법
Fig. 1 Flowchart of load flow analysis for DC railway power system
TPS에서는 전기철도 선로의 각종 정수를 기준으로 해당 선로의 차량 운행에 필요한 주행 동력과 차량의 운행 패턴을 출력한다. 차량 한 대가 이동하면서
노선의 각 위치에 따른 속도와 필요 전력을 계산하고, 결괏값은 DCPS의 입력값으로 사용한다. TPS의 입력값으로 열차의 제원과 구배와 커브를 포함한
계통에 대한 정보, 그리고 속도폭, dwell 시간을 포함한 운전 데이터가 연산에 필요하다. TPS는 주행 동력을 계산하기 위해 뉴턴의 운동법칙을
기반하여 도출한다(8).
여기서 $F_{ETF}$는 유효 견인력(Effective Traction Force), $m_{dyn}$은 유효질량(Dynamic Mass), $F_{MTF}$는
차량 전동기 견인력(Motor Traction Force), $R_{T}$는 열차의 총저항이다. 열차 총저항은 주행저항($R_{R}$), 커브저항($R_{C}$),
그리고 구배저항($R_{G}$)으로 구성된다. 유효견인력이 계산되면 열차의 속도에 따라 기계적 Power를 계산할 수 있다. 전기적 Power는 기계적
Power에 인버터, 모터 그리고 기어의 전체 효율을 곱하거나 나누어 견인력(Traction Power,$P_{t}$)과 제동력(Braking Power,$P_{b}$)을
구할 수 있다.
기계적 Power가 0보다 작은 경우 제동력, 0보다 큰 경우 역행력으로 구분한다. 여기서, 견인력 효율은 $T_{eff}$, 제동력 효율은 $B_{eff}$로
나타낸다. 차량 자체에서 공조기나 조명설비 등에 필요한 전력은 일정하므로 상수로 고정하고 보조 동력($P_{aux}$)으로 정의한다. 역행력, 제동력
그리고 보조 동력까지 모두 고려한 차량에 필요한 Power($P_{tra\in}$)는 식 (6)과 같다(9).
2.2 DC Power Simulation
TPS에서 도출된 노선상 열차의 위치에 따라 필요한 Power가 계산되면 여러 대의 차량을 Headway에 맞춰 운행할 때 DCPS를 통해 조류해석이
가능해진다. 일반적인 전력 시스템과는 다르게 전기철도의 급전 계통은 부하의 특성에서 차이가 있다. 부하가 시공간적으로 변하고, 회생제동 시 전원 소스로
동작한다는 특성이 존재한다. 조류해석에 앞서 직류 급전 계통을 Norton 등가회로로 표현하여 Nodal Analysis를 통해 해석할 수 있다(10).
차량 부하 특성상 비선형성을 가지고 있어 노드 방정식을 세우면 차량 노드의 전류로 인해 선형 식이 아닌 비선형 식으로 표현된다. 이러한 이유로 iteration을
이용한 반복법을 통해 차량 부하에 대해 가정치 전류를 설정 후 전압을 구하는 방식을 사용하여 전압이 오차 이내 들어올 경우 iteration을 끝내는
방식을 적용하고 있다(11).
3. 전압 안정화 및 회생 에너지 활용률을 고려한 최적워치 선정방법
3.1 회생 에너지와 전압의 상관관계 분석
회생 에너지 활용설비들은 Threshold Voltage($V_{th}$)를 기준으로 동작을 수행한다. 즉, 각 변전소에서 발생되는 전압 발생빈도,
최대 전압과 같이 설비들의 동작과 관련되므로 고려해야 한다. 특히, 계통의 안정성을 고려하여 집전전압의 국제 허용 규정치 기준인 IEC 60850의
최고 영구 전압 및 최고 비영구 전압 기준을 넘지 않아야 한다. 그에 따라, 변전소에서 발생되는 전압과 회생 에너지 발생량 간의 관계가 우선적으로
검토되어야 한다.
전압 발생빈도, 최대 전압 그리고 회생 에너지 발생량 차이에 대한 검토가 필요하다. 국내 실제 운영노선인 A 도시철도 노선의 데이터를 활용하여 TPS
및 DCPS로 분석하였다. A 노선은 750V 공칭 전압을 사용하며, 12개의 변전소가 설치되어 있고, 회생 발생량, 최대 전압 그리고 IEC 60850
기준 영구 전압 초과 빈도수에 따른 변전소 순위는 표 1과 그림 2와 같다.
표 1 구분에 따른 변전소 순위
Table 1 Substation ranking by classification
Substation Number
|
Regenerative energy
|
Maximum voltage(V)
|
Voltage counts
|
1
|
10
|
12
|
10
|
2
|
8
|
11
|
6
|
3
|
9
|
10
|
6
|
4
|
7
|
9
|
6
|
5
|
12
|
8
|
2
|
6
|
11
|
2
|
10
|
7
|
6
|
1
|
10
|
8
|
5
|
5
|
9
|
9
|
3
|
6
|
2
|
10
|
4
|
7
|
2
|
11
|
2
|
3
|
2
|
12
|
1
|
4
|
1
|
그림 2 분야에 따른 변전소별 발생량 a) 회생 에너지 b) 최대전압 c) 최고영구전압 초과 빈도수
Fig. 2 Amount generated by substation according to field a) Regenerative energy b)
Maximum voltage c) Frequency of exceeding maximum permanent voltage
그림 2와 표 1과 같이 1시간의 시뮬레이션 결과 전압 또는 회생 에너지 발생량에 따라 순위에 차이는 있지만, 상관관계가 존재한다. 최대전압, 전압 빈도수, 그리고
회생 에너지 발생량에 따른 상관관계를 비교한 결과는 표 2와 같다.
표 2 변수별 상관관계 분석결과
Table 2 Correlation analysis result by variable
$V_{\max}$and $F_{V}$
|
$V_{\max}$and $P_{re\ge n}$
|
$F_{V}$and $P_{re\ge n}$
|
0.0962
|
0.4171
|
0.5984
|
전압 빈도수와 회생 에너지 발생량 사이 높은 상관관계가 존재하고, 적절한 회생 에너지 활용에 따라 전압 안정화를 도모할 수 있다. 즉, 회생 에너지를
활용하기 위한 설비들의 배치가 계통 안정성과 운영요금 절감에 필수적이다. 따라서 본 논문에서는 회생 에너지 활용량에 따른 최저 소거법(Minimum
Elimination Method, MEM) 방식을 소개한다.
3.2 최저 소거법을 이용한 회생 에너지 활용설비 설계
계통에 대한 조류해석을 TPS 및 DCPS를 통해 수행 후 계통 내 변전소 회생 에너지 발생량과 순위를 도출하여 설치될 설비의 용량을 계산할 수 있다.
실제 사용 가능한 회생 에너지의 순시 최댓값 $P_{av,\:\max}$은 식 (7)과 같이 계산된다.
여기서, $P_{r,\:\max}$는 계통에서 발생한 회생 에너지 순시 값의 최대치이며, $P_{loss}$는 선로손실 및 차량 자체에 필요한 공조설비
등의 보조 동력을 합한 값이다. 경제적인 요건 등을 고려하여 설치 가능한 회생 설비 수 $N$이 정해지면, 가장 회생 에너지 활용이 낮은 변전소를
제거 후 다시 조류해석을 수행한다. 순차적으로 $N$개의 설비가 남을 때까지 순차적으로 제거하여 설치 위치를 선정할 수 있다. 표 1처럼 회생 발생량 순위(Maximum Selection Method, MSM)로 6개의 변전소를 선정할 경우, 7~12번 변전소에 설치된다. 이처럼
모두 인접한 변전소로 설치가 될 경우, 전류 분배법칙에 의해 회생 전류 $i_{r}$는 양옆 변전소로 식 (8),(9)와 같이 분배된다.
여기서, $i_{sn}$은 n번 변전소로 흐르는 전류이며, $R_{sn}$은 차량 위치부터 n번 변전소까지의 전차선 저항을 나타낸다. 회생 에너지
활용설비가 연속으로 몰려있게 될 경우, 회생 전류의 전류 분배가 인접한 두 변전소로 흐르므로 설비의 용량을 넘는 회생 에너지가 발생할 수 있다. 변전소가
한쪽으로 몰리는 현상을 배제하기 위해 가운데를 기준으로 나누어 각 구간별 순위를 도출하여 선택하는 Division Selection(DS) 방법이
존재한다. MEM을 이용하게 되면 최저 변전소부터 제거해나가며 전류 분배가 상대적으로 멀어지는 변전소로 전류의 흐름이 상대적으로 적게 발생된다. 인접한
변전소들은 회생 에너지 활용량이 높아지게 되고, 상대적으로 낮은 다른 위치의 변전소를 소거해나감으로써 모두 인접한 위치의 변전소 설계를 배제할 수
있다. 각 방법에 대한 개략도와 MEM에 대한 순서도는 그림 3과 4와 같다.
그림 3 방법에 따른 회생 에너지 활용설비 설계방식 a) Maximum Selection Method b) Division Selection Method,
c) Minimum Elimination Method
Fig. 3 Design method of regenerative energy utilization facility according to method
그림 4 최대 회생 에너지 활용 선정 방법 순서도
Fig. 4 Flowchart of selection method for maximum regenerative energy utilization
4. 시뮬레이션 및 결과
4.1 DCPS를 이용한 최저 소거법 시뮬레이션
제안하는 방법에 대한 검증을 위해 국내 실제 운영되고 있는 도시철도 노선 A의 데이터를 활용하였다. A 노선은 27개의 정거장과 12개의 변전소로
구성되어 있으며, DC 750V의 공칭 전압을 사용한다. 자세한 TPS와 DCPS에 사용된 파라미터는 표 3에 나타냈다. 시뮬레이션 결과, 그림 3a와 같이 5번 변전소에서 가장 낮은 회생 에너지 발생량을 보였고, 순차적으로 N=6이 될 때까지 최저 소거법을 이용하여 소거한 결과는 표 4와 같다.
표 3 TPS 및 DCPS 시뮬레이션 파라미터 입력값
Table 3 Simulation parameters for TPS and DCPS
Simulation Parameters
|
Values
|
Dwell Time [s]
|
20
|
Auxiliary power of train [kW]
|
96
|
Headway [s]
|
360
|
Train mass [ton]
|
75.78
|
Motion resistance [kN]
|
18.2966+1.2666$v$+0.09462$v^{2}$
|
Maximum train speed [km/s]
|
80
|
Train speed width [km/s]
|
5
|
No-load voltage [V]
|
804
|
Simulation Time [s]
|
3,600
|
N
|
6
|
표 4 최저 소거법을 이용한 시뮬레이션 결과
Table 4 Simulation result of Minimum Elimination method
Sub. Number
|
Iteration 1
|
Iteration 2
|
Iteration 3
|
$P_{re\ge n}$
|
Rank
|
$P_{re\ge n}$
|
Rank
|
$P_{re\ge n}$
|
Rank
|
1
|
20.22
|
11
|
0
|
11
|
0
|
10
|
2
|
20.473
|
10
|
24.197
|
2
|
24.963
|
3
|
3
|
20.698
|
8
|
23.3
|
4
|
24.212
|
6
|
4
|
21.416
|
4
|
23.184
|
5
|
24.345
|
5
|
5
|
0
|
12
|
0
|
11
|
0
|
10
|
6
|
20.534
|
9
|
21.525
|
8
|
23.272
|
9
|
7
|
20.993
|
6
|
21.678
|
7
|
23.98
|
7
|
8
|
20.762
|
7
|
21.313
|
10
|
0
|
10
|
9
|
21.435
|
3
|
21.903
|
6
|
24.361
|
4
|
10
|
21.057
|
5
|
21.451
|
9
|
23.515
|
8
|
11
|
23.706
|
2
|
24.043
|
3
|
25.805
|
2
|
12
|
25.843
|
1
|
26.162
|
1
|
27.831
|
1
|
Sub. Number
|
Iteration 4
|
Iteration 5
|
Iteration 6
|
$P_{re\ge n}$
|
Rank
|
$P_{re\ge n}$
|
Rank
|
$P_{re\ge n}$
|
Rank
|
1
|
0
|
9
|
0
|
8
|
0
|
7
|
2
|
26.7
|
5
|
27.687
|
6
|
33.292
|
2
|
3
|
26.25
|
7
|
27.432
|
7
|
0
|
7
|
4
|
26.892
|
3
|
28.386
|
5
|
33.05
|
3
|
5
|
0
|
9
|
0
|
8
|
0
|
7
|
6
|
0
|
9
|
0
|
8
|
0
|
7
|
7
|
26.784
|
4
|
29.541
|
4
|
31.952
|
6
|
8
|
0
|
9
|
0
|
8
|
0
|
7
|
9
|
26.377
|
6
|
30.22
|
3
|
32.031
|
5
|
10
|
25.2
|
8
|
0
|
8
|
0
|
7
|
11
|
27.235
|
2
|
31.33
|
2
|
32.722
|
4
|
12
|
29.184
|
1
|
33.044
|
1
|
34.352
|
1
|
4.2 Maximum Selection 방식과 Division Selection 방식과의 비교
DS 방식(7~12번), DS 방식(2,3,4,9,11,12) 그리고 최저 소거법을 이용한 방법(2,4,7,9,11,12번)의 회생 에너지 발생량을
비교한 결과는 표 5와 같다. MS 방식은 약 176.5 kWh, DS 방식은 195.6 kWh, 그리고 ME 방식은 약 197.5 kWh로 기존 MS 방식 대비 약
11.9%의 에너지 활용률을 높일 수 있었다.
표 5 방식별 변전소 회생에너지 발생량
Table 5 Substation regenerative energy generation by methods
Sub. no.
|
MS Method
(kWh)
|
DS Method
(kWh)
|
ME Method
(kWh)
|
1
|
0
|
0
|
0
|
2
|
0
|
29.7
|
33.3
|
3
|
0
|
29.8
|
0
|
4
|
0
|
31.4
|
33.1
|
5
|
0
|
0
|
0
|
6
|
0
|
0
|
0
|
7
|
31.9
|
0
|
32.0
|
8
|
29.3
|
0
|
0
|
9
|
28.7
|
34.3
|
32.0
|
10
|
27.1
|
0
|
0
|
11
|
28.8
|
34.5
|
32.7
|
12
|
30.7
|
35.9
|
34.4
|
Total
|
176.5
|
195.6
|
197.5
|
5. 결 론
본 논문에서는 회생 에너지 활용률 향상을 위해 최저 소거법을 이용한 직류 전기철도 회생 에너지 활용설비 설계방법에 대해 소개하였다. TPS를 이용하여
분석하려는 노선에 대한 차량 운행 프로파일과 필요 차량 전력값을 도출한 후, DCPS를 통해 운행시격에 맞추어 차량을 운행하며 전체 계통에 대한 조류해석
방법을 먼저 설명하였다. 계통의 안정성과 운영비용 저감을 목적으로 최대 전압, IEC 60850 기준 영구 전압 초과 빈도수, 그리고 회생 에너지
발생량에 대해 상관관계를 도출한 결과, 빈도수와 회생 에너지 발생량 간 0.5984의 높은 상관관계를 보였다. 회생 에너지의 시의적절한 활용을 통한
계통 안정성 향상을 도모하기 위해 MEM을 소개하였다. 최대 회생 에너지 발생량을 기준으로 설치하는 방법의 단점을 개선하고자 제일 낮은 회생 에너지
활용 변전소를 1개씩 소거해가며 조류해석을 TPS와 DCPS를 통해 반복적으로 iteration을 수행하였다. 그 결과, 기존 방식 대비 11.9%의
에너지 활용률을 증가시킬 수 있었다.
Acknowledgements
본 연구는 국토교통부/국토교통과학기술진흥원의 지원으로 수행되었음(과제번호 RS-2022-00143874).
References
H. Jung, 2022, An Optimal Charging and Discharging Scheduling Algorithm of Energy
Storage System to Save Electricity Pricing Using Reinforcement Learning in Urban Railway
System. J. Electr. Eng. Technol, 17, pp. 727-735
C.-M. Yun, G.-J. Cho, H. Kim, H. A. Jung, 2022, Study on the Train Brake Position-Based
Control Method for Regenerative Inverters, Energies, Vol. 15
H.J. Kim, H. Jung, Y.J. Ko, E.S. Chae, H.J. Kim, I.S. Hwang, J.H. Heo, J.Y. Park,
2021, Cooperative Operation Schedules of En-ergy Storage System and Demand Response
Resources Considering Urban Railway Load Characteristic under a Time-of-Use Tariff,
J. Electr. Eng. Technol, Vol. 16, pp. 1273-1284
H. Cho, J. Kim, H. Jung, H. Kim, 2022, Simultaneous DC Railway Power System Analysis
Method Using Model-Based TPS, Appl. Sci, 12, 6929
S. Su, X. Li, T. Tang, 2013, A subway train timetable optimization approach based
on energy-efficient operation strategy, IEEE Trans. Intell. Transp. Syst, Vol. 14,
pp. 883-893
S. Lu, P. Weston, S. Hillmansen, H.B. Gooi, C. Roberts, 2014, Increasing the regenerative
braking energy for railway vehicles, IEEE Trans. Intell. Transp. Syst, Vol. 15, pp.
2506-2515
Y. Cai, M.R. Irving, S.H. Case, 1995, Iterative techniques for the solution of complex
DC-rail-traction systems including regener-ative braking, IEEE Proc.-Gener. Transm.
Distrib, Vol. 142, pp. 445-452
H.S. Lee, 2010, Development of a Novel Powerflow Algorithm for Energy Storage Optimization
in DC Electric Railway Systems, Ph.D. Thesis, Korea University, Seoul, Korea
S.G. Chung, S.J. Lee, 2003, The computer algorithm for DC traction power supply system
analysis including regenerative brak-ing vehicles, Trans. Korean Inst. Electr. Eng,
Vol. 52, pp. 639-646
J. Włodzimierz, S. Adam, 2018, The multi-criteria optimization method for implementation
of a regenerative inverter in a 3kV DC traction system, Electr. Power Syst. Res, Vol.
161, pp. 61-73
A. Gonzalez-Gil, R. Palacin, P. Batty, 2013, Sustainable urban rail systems: Strategies
and technologies for optimal management of regenerative braking energy, Energy Convers.
Manag, Vol. 75, pp. 374-388
저자소개
He received a B.S. degrees in 2019, from the College of Electric and Electrical Enginnering,
Seoul, Hongik University.
He reveived a M.S. degree from the College of Electrical and Computer Engineering,
Sungkyunkwan University, Suwon, South Korea.
At present, he is enrolled in the doctor’s program in the Department of Transportation
Engineering, Korea University of Science and Technology, Republic of Korea.
His research interests include regenerative energy efficiency enhancement and traction
power system analysis.
He received BS degree in electric engineering from Sungkyunkwan University.
He is currently a senior research engineer with ENTEC Electric & Electronic Co., Ltd.
His research interests are protection, control, communication monitoring system for
power electric system.
She received MS degree in electric engineering from Seoul National Univ. of Science
& Technology.
She is currently a research engineer with ENTEC Electric & Electronic Co., Ltd.
Her research interests are power system protection engineering and its digital processing
algorithms.
He received MS degree in electric engineering from Sungkyunkwan University.
He is currently a senior research engineer with ENTEC Electric & Electronic Co., Ltd.
His research interests are protection algorithms for power electric system and its
analysis.