Enhancing Power Grid Sustainability with Environment-Friendly and Resilient Operation Strategies
LEI Shunbo
- Dept. of Electrical & Electronic Eng.
The University of Hong Kong
April 25, 2018
with Environment-Friendly and Resilient Operation Strategies LEI - - PowerPoint PPT Presentation
Enhancing Power Grid Sustainability with Environment-Friendly and Resilient Operation Strategies LEI Shunbo Dept. of Electrical & Electronic Eng. The University of Hong Kong April 25, 2018 Outline Background and motivation
April 25, 2018
Sustainable power grids Environmental sustainability Economic sustainability Social sustainability Environment-friendly
Resilient operation RES integration Emission issue …… ……
Coal Intensive Power Systems,” IEEE Trans. Power Syst., vol. 28, no. 1, pp. 236-245, Feb. 2013.
Power Syst., vol. 21, no. 1, pp. 341-347, Feb. 2006. …
2
b b it i it i it i
x y z W i n d d i r e c t i
Δh x y z W i n d d i r e c t i
(a) (b)
H
r
2 2 2 2 2 2 2
1 exp 2 2
y it ijt y z x y x z x
d E H C w I I d I d I d
Considers meteorological conditions and the system’s geographical distribution, to estimate GLAPC due to emissions of generations.
Consider population density & background pollution Minimize people’s exposure to air pollutants
min( ) CC BDC APDC
1 1
( )
G
N T i it i it t i
CC su SU sd SD
2 1 1
[ ( ) ]
G
N T b b i it i it i it t i
BDC P P O
2 1 1 1
[ ( ) ]
G D
N N T b b j ij i it i it i it t j i
APDC K c P b P a O
2 j j j
s NP BP Commitment cost: start-up/shut-down of generators Base-case dispatch cost: real power outputs of generators
Uncertainty set of wind power outputs. Robust feasibility check: ensure adequate spinning reserve, to avoid environmentally unfriendly operations of more polluting coal-fired units in real time.
1) Set the number of iterations K=0. Choose a tolerance δ(>0) for the robust feasibility check. 2) Solve the MP to update the current optimal solution (x*, y1
*, y2,1 *,…, y2,K *).
3) K=K+1. Solve the SP with xK = x*, to obtain uK. 4) If objective_of_SP ≤ δ, return (x*, y1
*) and terminate. Otherwise, go to step 2).
A mixed-integer quadratic programing problem with a robust feasibility check
2 4 6 8 10 12 50 60 70 AUCC (k$) PGLAPC (µg/m3) 28 28.5 29 PGLAPC without wind AGLAPC without wind AGLAPC (µg/m3) PGLAPC with wind AGLAPC with wind
UC cost k$ Total emission Ton Peak GLAPC µg/m3 Average GLAPC µg/m3 Total exposure ·105µg/m3 Case 1 224.58 100.25 70.32 29.37 150.92 Case 2 200.19 93.43 70.32 29.22 150.42 Case 3 213.38 65.00 64.79 28.29 148.09 Case 4 209.79 67.90 56.07 28.10 145.81 Bus No. Back-ground PM2.5 µg/m3 BPj NPj Peak GLAPC µg/m3 Average GLAPC µg/m3 Case No. Case No. 1 2 3 4 1 2 3 4 4 22 3 9203 22.00 22.00 22.00 22.00 20.00 22.00 22.00 22.00 5 50 5 1463 70.32 70.32 64.79 56.07 59.97 59.15 56.04 52.53 7 22 3 320 22.00 22.00 22.00 22.00 22.00 22.00 22.00 22.00 9 22 3 5680 22.00 22.00 22.00 22.00 22.00 22.00 22.00 22.00 10 30 4 1733 30.02 30.02 30.02 30.01 30.01 30.01 30.00 30.00 11 36 5 674 53.28 53.28 53.28 53.28 46.14 46.01 46.13 46.75 12 23 3 1174 32.36 32.36 31.99 31.46 28.60 28.34 27.71 27.61 13 18 3 2600 21.02 21.02 20.97 20.90 20.57 20.37 19.73 19.97 14 13 2 2869 13.13 13.13 13.13 13.13 13.06 13.06 13.02 13.06
Case 1: w/o wind power, w/o APDC Case 2: w/ wind power, w/o APDC Case 3: w/ wind power, w/o APDC, w/ total emissions limit Case 4: w/ wind power, w/ APDC
Emission issues need to be explicitly considered, so as to utilize wind power’s benefits in
air pollution control.
Only limiting total emissions is less cost effective.
Again: 1) need to explicitly consider emission issues; and 2)
effective. Observation : wind power makes a system more cost effective in and more capable of air pollution control.
UC cost ·107 $ Total emission Ton Peak GLAPC µg/m3 Average GLAPC µg/m3 Total exposure ·108µg/m3 Case 1 1.6254 26640.70 74.36 58.88 133.59 Case 2 1.4596 24183.33 74.36 55.43 124.21 Case 3 1.5290 11500.00 74.36 33.60 66.54 Case 4 1.5158 19228.48 24.36 14.52 16.93
23.0 29.5 36.0 42.5 49.0 55.5 62.0 68.5 74.0 1 2 3 4 5 6 7 8 9 1013.0 19.5 26.0 32.5 39.0 45.5 52.0 58.5 PGLAPC without wind AGLAPC without wind PGLAPC with wind AGLAPC with wind AUCC (105$) PGLAPC (µg/m3) AGLAPC (µg/m3)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
Control center Substation
Communica- tion server Feeder terminal unit Remote-controlled switch Feeder breaker Remote terminal unit Communication network Workstation Database and server
Optimization model formulation Single/multiple objectives ANM strategies coordination Solution algorithm Hybrid Benders/C&CG Strong duality v.s. KKT Uncertainty consideration Stochastic characteristics Time-series issues …… …… ……
Obj.: minimize DG curtailments Number of selected critical switches
Switch-type-dependent operation constraints
Intro-day DSDR constraints
Innermost min: seek the optimal intro-day DSDR strategy to achieve MDGC, using selected switches, with given loads and DG outputs. Mid-level max: obtain loads and DG outputs that maximize MDGC. Outermost min: find critical switches that minimize the maximum MDGC.
Uncertainty set of loads Uncertainty set of DG outputs Modified obj., oriented by the worst case
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 DG DG DG DG DG DG
20 40 60 80 100 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 % Time Period Wind Solar 8.58 6.36 5.32 5.21 4.98 4.53 4.53 8.58 6.88 5.82 5.55 5.71 5.03 4.53 4 5 6 7 8 9 10 2 4 6 8 10 37 W-MDGC (MWh) Value of Π𝑆𝐷𝑇 Robust Model Deterministic Model
Robust Model Deterministic Model 2 12-13*, 14-15 12-13*, 16-17 4 6-26*, 27-28*, 12-13, 16-17 6-26*, 27-28*, 9-15, 17-18 6 6-26*, 7-8*, 16-17*, 27-28*, 12-13, 14-15 6-26*, 7-8*, 16-17*, 27-28*, 9-15, 20-21 8 20-21*, 27-28*, 2-3, 6-26, 12-13, 14-15, 16-17, 21-8 20-21*, 27-28*, 5-6, 6-7, 9-15, 17-18, 26-27, 29-30 10 6-26*, 8-9*, 9-15*, 12-13*, 27-28*, 2-19, 14-15, 16-17, 20-21, 21-8 6-26*, 8-9*, 9-15*, 12-13*, 27-28*, 2-3, 7-8, 11-12, 17-18, 24-25
Open Switches 2 9-15, 11-12, 14-15, 20-21, 27-28 4 6-26, 9-15, 11-12, 16-17, 20-21 6 7-8, 9-15, 11-12, 20-21, 27-28 8 2-3, 6-26, 9-15, 11-12, 16-17 10 6-26, 8-21, 9-15, 11-12, 20-21 Time Period Switching Actions 2 Hour 4 Open: 12-13. Close: 14-15. Hour 14 Open: 14-15. Close: 12-13. 4 Hour 4 Open: 12-13. Close: 16-17. Hour 6 Open: 16-17. Close: 12-13. Hour 22 Open: 27-28. Close: 6-26. Hour 23 Open: 6-26. Close:27-28 6 Hour 14 Open: 6-26, 14-15. Close: 7-8, 27-28. Hour 15 Open: 7-8, 27-28. Close: 6-26, 14-15. Hour 18 Open: 16-17. Close: 7-8. Hour 21 Open: 6-26. Close: 27-28. Hour 22 Open: 7-8, 27-28. Close: 6-26, 16-17. 8 Hour 4 Open: 21-8, 27-28. Close: 6-26, 16-17. Hour 6 Open: 14-15. Close: 2-3. Hour 14 Open: 6-26. Close: 27-28. Hour 15 Open: 16-17. Close: 14-15. Hour 17 Open: 20-21. Close: 21-8. Hour 21 Open: 2-3. Close: 20-21. 10 Hour 5 Open: 12-13. Close: 20-21. Hour 14 Open: 14-15. Close: 12-13. Hour 15 Open: 12-13, 16-17. Close: 14-15, 9-15. Hour 16 Open: 8-9. Close: 16-17. Hour 17 Open: 16-17. Close: 12-13. Hour 21 Open: 20-21. Close: 8-9. Hour 23 Open: 9-15. Close: 16-17.
Observed Outages to the Bulk Electric System (Source: Energy Information Administration) U.S. 2014 Billion-Dollar Weather and Climate Disasters (Source: National Oceanic and Atmospheric Administration)
Event Date Region/Division/State Customers Affected Superstorm Sandy October 2012 Northeast 8,100,000 Derecho July 2012 Middle Atlantic 4,200,000 Early season snow October 2011 New England 3,000,000 Tropical Storm Irene August 2011 Middle Atlantic 3,200,000 Wildfires July 2012 California, Colorado 2,000,000 Windstorm November 2011 Southern California 400,000 2011-2012 Natural Disasters and Reported Customers Affected by Power Outages
The GridWise Alliance, “Improving electric grid reliability and resilience: Lessons learned from Superstorm Sandy and other extreme events,” 2013 [Online].
3 2
[2] NERC, “Hurricane Sandy Event Analysis Report,” 2014 [Online].
Figure 1: Damaged high voltage transmission line Figure 2: Damaged substation equipment Figure 3: Damaged distribution system equipment Figure 4: Crews clearing miles of road to access and repair a line
1 4
2 3 4 5 1 6 7 8 2 3 4 5 1 6 7 8 2 3 4 5 1 6 7 8 2 3 4 5 1 6 7 8 2 3 4 5 1 6 7 8 2 3 4 5 1 6 7 8 (a) (b) (c) (d) (e) (f) Closed branch Open branch Load Breaker Closed switch Open switch
Global optimum & more efficient computation
C.-S. Chen, C.-H Lin, H.-J. Chuang, C.-S. Li, M.-Y. Huang, and C.-W. Huang, “Optimal placement of line switches for distribution automation systems using immune algorithm,” IEEE Trans. Power Syst., vol. 21,
sectionalizing switches in distribution networks using a fuzzy multiobjective approach,” IEEE Trans. Power Del., vol. 24, no. 1, pp. 268–276, Jan. 2009.
particle swarm optimization algorithm,” IEEE Trans. Power Del., vol. 23, no. 1, pp. 271-279, Jan. 2008.
by genetic algorithm,” Elect. Power Syst. Res., vol. 31, no. 2, pp. 97–102, Nov. 1994. ……
placement strategy in distribution systems,” IEEE Trans. Power Del., vol. 27, no. 1, pp. 362-370, Jan. 2012.
automation considering earth fault events,” IEEE Trans. Smart Grid, vol. 6, no. 2, pp.1010-1018, Mar. 2015.
Obj.: max restored loads & min switching actions Only allow RCSs to be utilized Reconfiguration constraints
Obj.: max reduction of customer interruption cost, & min maintenance and operation cost. Loads that can be restored: constrained by the number and locations of RCSs. Obj.: max reduction of system average interruption duration index, & min M&OC. Number of customers that can be restored: constrained by the number and locations of RCSs. Loads that can be restored: constrained by the number and locations of RCSs. Restorable loads using RCSs: equal to restorable loads using all switches (manual/remote-controlled). Obj.: min the number of RCSs to be installed.
Integer variables: from to Computational complexity: from to
Nonlinear New integer variables
424 in total (compared to 2^37 possible combinations of switching actions)
A small number of RCSs can reduce distribution system CIC effectively. The solution of this method is robust to the changes of CIC parameters.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 1 2 8 9 15 56 57 58 59 60 50 54 117 118 119 121 52 53 55 116 36 51 33 49 70 114 115 122 47 48 71 113 35 34 28 46 109 111 112 123 27 45 69 110 30 29 43 108 106 44 68 105 107 26 25 38 42 103 39 41 67 104 24 23 77 37 62 72 64 63 65 73 21 20 74 66 75 76 7 19 99 18 17 101 97 5 95 93 94 88 87 90 6 102 98 14 16 100 96 4 92 89 86 91 61 78 79 80 81 82 83 84 85 31 32 3 12 11 10 13 22 120 40
Fault location 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Number of PCRSs 2 83 74 11 9 8 4 5 5 4 2 2 1 1 Fault location 17 19 20 21 23 24 26 27 28 29 30 31 32
1 10 2 1 15 1 36 33 31 14 39 29 1
RCIC ($/yr) Optimal RCS locations 25% 4 2467.31 7257.55 8-9, 12-22, 18-33, 25-29 50% 5 3084.14 15571.38 8-9, 12-22, 18-33, 25-29, 30-31 100% 5 3084.14 31142.77 8-9, 12-22, 18-33, 25-29, 30-31 200% 5 3084.14 62285.53 8-9, 12-22, 18-33, 25-29, 30-31 400% 7 4317.80 124727.83 6-7, 8-21, 11-12, 12-22, 18-33, 25-29, 31-32
With the number of RCSs increasing, upgrading
marginal benefit in terms of SAIDI reduction.
Different objectives lead to quite different optimal RCS allocation schemes. With a specific objective, a proper model should be developed and applied (or, compare results of different models before finalize the decision).
5 10 15 20 25 30 35 40 45 50 45 40 35 30 25 20 RCS cost (102$/yr) RSAIDI (min/yr)
Allocation model Allocation scheme RCIC minus RCS cost ($/yr) RSAIDI (min/yr) Restoration capability CIC-oriented 8-9, 12-22, 18-33, 25-29, 30-31 28058.63 46.89 24/27 SAIDI-oriented 6-7, 8-21, 18-33, 25-29 27351.96 45.96 18/27 RL-oriented 9-10, 9-15, 12-22, 14-15, 18-33, 25-29, 28-29, 30-31 26249.58 47.14 27/27
FEMA, “Mitigation assessment team report: Hurricane Sandy in New Jersey and New York,” 2013 [Online].
A MEG supplying power for a data center after Hurricane Sandy (Manhattan, NY)
The GridWise Alliance, “Improving electric grid reliability and resilience: Lessons learned from Superstorm Sandy and
Case of Hurricane Sandy D. Barrett, “Few big FEMA generators humming,” 2012.
Before Sandy struck (Oct. 29) 400 industrial-size generators prepared. 3 days afterwards (Nov. 2) Only 50 of them working, while millions still without electricity. Connecticut: 1/51 MEGs working, 300,000 consumers without power. New York State: 7/13 MEGs working.
2011 Int. Conf. on Electr. Inf. and Control Eng., Apr. 2011, pp. 5419-5422.
telecommunications network," in Proc. 5th Int. Telecom. Energy Conf., Oct. 1983, pp. 144-150.
electric network," in Proc. 20th Int. Conf. and Exhib. on Electr. Distrib., Jun. 2009, pp. 1-4.
[Online] www.duquesnelight.com
and difficult to forecast system conditions accurately (statistical/simulation-based models).
S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 29 21 22 23 24 25 26 27
Modified 27-node system from: S. Civanlar, J. J. Grainger, H. Yin, and S.
(Pre-installed microgrids V.S. dynamically formed microgrids) (Pre-installed distributed generators V.S. dynamically dispatched MEGs)
Major decisions MEGs allocation: allocate each MEG to which node. Microgrids formation:
Natural disaster struck
Pre-position Real-time allocation Timeline
Weather/conditions forecast and monitor
Distribution system & road network damages monitor & assessment Week/days
Days/day/hours
Minutes/hours/day/days afterwards
S1 S2 S3 S4 100 102 106 108 112 114 150 116 118 120 122 110 124 126 156 154 152 132 188 190 158 164 130 128 148 146 144 142 134 160 136 140 162 138 299 275 225 201 202 203 227 213 212 205 242 244 228 229 209 230 231 232 208 233 234 237 238 211 241 240 210 236 235 206 220 204 214 218 207 224 222 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 S1 Substations
Distribution nodes Candidate distribution nodes for MEG connection DS1 DS2 DS3 DS4 Normally closed lines Normally closed lines
S1 S2 S3 S4 100 102 106 108 112 114 150 116 118 120 122 110 124 126 156 154 152 132 188 190 158 164 130 128 148 146 144 142 134 160 136 140 162 138 299 275 225 201 202 203 227 213 212 205 242 244 228 229 209 230 231 232 208 233 234 237 238 211 241 240 210 236 235 206 220 204 214 218 207 224 222 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 S1 Substations
Distribution nodes Candidate distribution nodes for MEG connection DS1 DS2 DS3 DS4 Normally closed lines Normally closed lines
Geographic information of road network & distribution systems Distribution systems
Boundaries
Candidate distribution nodes for MEG connection Staging locations
DS1 DS4 DS2 DS3
Road edges Intersections
S1 S2 S3 116 132 136 203 204 234 306 311 403 413 419
( lg( )) O E V V
min{ ( , ): , }
n n n n n
f
x y x Λ y Ω min{ ( , ): , , }
n n n n n n n n n n
f
x y x Λ y Ω A x h
Boundaries
Candidate distribution nodes for MEG connection Staging locations
DS1 DS4 DS2 DS3
Road edges Intersections
S1 S2 S3 116 132 136 203 204 234 306 311 403 413 419
S1 S2 S3 S4 100 102 106 108 112 114 150 118 120 122 110 124 126 156 154 152 132 188 190 158 164 130 128 148 146 144 142 134 160 140 162 138 299 275 225 201 202 203 227 213 212 205 242 244 228 229 209 230 231 232 208 233 234 237 238 211 241 240 210 236 235 206 220 214 218 207 224 222 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 401 402 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 S1 Substations
Distribution nodes Candidate distribution nodes for MEG connection DS1 DS2 DS3 DS4 Normally closed lines Normally closed lines
MG 204 MG 413 MG 403 403 204 136 116
Minimize sum of (load_priorities*load_size*outage_time) Minimize sum of unserved loads MEG Capacity utilization rate Times of CUR≥85% Capacity utilization rate Times of CUR≥85% 1 73.03 % 221/500 74.58 % 228/500 2 64.13 % 97/500 80.87 % 252/500 3 75.93 % 242/500 87.28 % 323/500 4 63.65 % 149/500 71.36 % 221/500 5 60.41 % 126/500 72.14 % 179/500 Average 67.43 % 167/500 77.25 % 240.6/500
—— Average amount of loads restored by MEGs Average travel time to assigned locations Non-optimal pre-position w/ traffic consideration 3866.8 kW 28.26 minutes Optimal pre-position w/ traffic consideration 3910.9 kW 21.58 minutes w/o traffic consideration 4082.5 kW 35.05 minutes