%% Please email questions, comments, and significant results to: Robert.C.Green@gmail.com. Thanks! %% %% (c) 2009-2010 Zhu Wang, Robert Green %% %% ALL RIGHTS RESERVED WORLDWIDE %% %% THIS PROGRAM IS FREEWARE. IT MAY BE COPIED AND %% DISTRIBUTED WITHOUT LIMITATION AS LONG AS THIS %% COPYRIGHT NOTICE AND REFERENCE %% INFORMATION BELOW ARE INCLUDED WITHOUT MODIFICATION, %% AND AS LONG AS NO FEE OR COMPENSATION IS CHARGED, %% INCLUDING "TIE-IN" OR "BUNDLING" FEES CHARGED FOR %% OTHER PRODUCTS. %% This code implements the IEEE-RTS79 model for Matpower %% References concerning this code include: %% 1 - R. Green, L. Wang, Z. Wang, M. Alam, and C. Singh, “Power System Reliability Assessment Using Intelligent State Space Pruning Techniques: A Comparative Study” Submitted to 2010 Conference on Power System Technology, Hangzou China, October 2010. %% 2 - R. Green, L. Wang, M. Alam, and C. Singh, “State space pruning for Reliability Evaluation using Binary Particle Swarm Optimization,” Submitted to Hawaii International Conference on System Sciences,University of Hawaii at Manoa, January 2011. %% 3 - R. Green, L. Wang, and C. Singh, “State space pruning for power system reliability evaluation using genetic algorithms,” IEEE Power & Energy Society General Meeting 2010, Minneapolis, MN, July 2010. %% 4 - R. Green, Z. Wang, L. Wang, M. Alam, and C. Singh, “Evaluation of loss of load probability for power systems using intelligent search based state space pruning,” The 11th International Conference on Probabilistic Methods Applied to Power Systems, Singapore, June 2010 function [baseMVA, bus, gen, branch, areas, gencost] = case79 %% system MVA base baseMVA = 100; %% bus data % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin bus = [ 1 2 108 22 0 0 11 1 0 138 1 1.1 0.9; 2 2 97 20 0 0 11 1 0 138 2 1.1 0.9; 3 4 180 37 0 0 11 1 0 138 1 1.1 0.9; 4 4 74 15 0 0 11 1 0 138 1 1.1 0.9; 5 4 71 14 0 0 11 1 0 138 1 1.1 0.9; 6 4 136 28 0 1 11 1 0 138 2 1.1 0.9; 7 2 125 25 0 0 11 1 0 138 2 1.1 0.9; 8 4 171 35 0 0 11 1 0 138 2 1.1 0.9; 9 4 175 36 0 0 11 1 0 138 3 1.1 0.9; 10 4 195 40 0 0 11 1 0 138 3 1.1 0.9; 11 4 0 0 0 0 11 1 0 230 3 1.1 0.9; 12 4 0 0 0 0 11 1 0 230 3 1.1 0.9; 13 3 265 54 0 0 12 1 0 230 4 1.1 0.9; 14 2 194 39 0 0 12 1 0 230 6 1.1 0.9; 15 2 317 64 0 0 12 1 0 230 6 1.1 0.9; 16 2 100 20 0 0 12 1 0 230 6 1.1 0.9; 17 4 0 0 0 0 12 1 0 230 7 1.1 0.9; 18 2 333 68 0 0 12 1 0 230 7 1.1 0.9; 19 4 181 37 0 0 12 1 0 230 5 1.1 0.9; 20 4 128 26 0 0 12 1 0 230 5 1.1 0.9; 21 2 0 0 0 0 12 1 0 230 7 1.1 0.9; 22 2 0 0 0 0 12 1 0 230 7 1.1 0.9; 23 2 0 0 0 0 12 1 0 230 5 1.1 0.9; 24 4 0 0 0 0 12 1 0 230 6 1.1 0.9; ]; %% generator data % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin gen = [ 1 20 0 10 0 1.035 24 1 20 0; 1 20 0 10 0 1.035 24 1 20 0; 1 76 14.1 30 -25 1.035 89 1 76 0; 1 76 14.1 30 -25 1.035 89 1 76 0; 2 20 0 10 0 1.035 24 1 20 0; 2 20 0 10 0 1.035 24 1 20 0; 2 76 7 30 -25 1.035 89 1 76 0; 2 76 7 30 -25 1.035 89 1 76 0; 7 100 17.2 60 0 1.025 118 1 100 0; 7 100 17.2 60 0 1.025 118 1 100 0; 7 100 17.2 60 0 1.025 118 1 100 0; 13 197 40.7 80 0 1.020 232 1 197 0; 13 197 40.7 80 0 1.020 232 1 197 0; 13 197 40.7 80 0 1.020 232 1 197 0; 15 12 0 6 0 1.014 14 1 12 0; 15 12 0 6 0 1.014 14 1 12 0; 15 12 0 6 0 1.014 14 1 12 0; 15 12 0 6 0 1.014 14 1 12 0; 15 12 0 6 0 1.014 14 1 12 0; 15 155 0.05 80 -50 1.014 182 1 155 0; 16 155 25.22 80 -50 1.017 182 1 155 0; 18 400 137.4 200 -50 1.050 471 1 400 0; 21 400 108.2 200 -50 1.050 471 1 400 0; 22 50 -4.96 16 -10 1.050 53 1 50 0; 22 50 -4.96 16 -10 1.050 53 1 50 0; 22 50 -4.96 16 -10 1.050 53 1 50 0; 22 50 -4.96 16 -10 1.050 53 1 50 0; 22 50 -4.96 16 -10 1.050 53 1 50 0; 22 50 -4.96 16 -10 1.050 53 1 50 0; 23 155 31.79 80 -50 1.050 182 1 155 0; 23 155 31.79 80 -50 1.050 182 1 155 0; 23 350 71.78 150 -25 1.050 412 1 350 0; ]; %% branch data % fbus tbus r x b rateA rateB rateC ratio angle status branch = [ 1 2 0.0026 0.0139 0.4611 193 200 200 0 0 1; 1 3 0.0546 0.2112 0.0572 208 220 220 0 0 1; 1 5 0.0218 0.0845 0.0229 208 220 220 0 0 1; 2 4 0.0328 0.1267 0.0343 208 220 220 0 0 1; 2 6 0.0497 0.1920 0.0520 208 220 220 0 0 1; 3 9 0.0308 0.1190 0.0322 208 220 220 0 0 1; 3 24 0.0023 0.0839 0.000 510 600 600 1.015 0 1; 4 9 0.0268 0.1037 0.0281 208 220 220 0 0 1; 5 10 0.0228 0.0883 0.0239 208 220 220 0 0 1; 6 10 0.0139 0.0605 2.459 193 200 200 0 0 1; 7 8 0.0159 0.0614 0.0166 208 220 220 0 0 1; 8 9 0.0427 0.1651 0.0447 208 220 220 0 0 1; 8 10 0.0427 0.1651 0.0447 208 220 220 0 0 1; 9 11 0.0023 0.0839 0.000 510 600 600 1.03 0 1; 9 12 0.0023 0.0839 0.000 510 600 600 1.03 0 1; 10 11 0.0023 0.0839 0.000 510 600 600 1.015 0 1; 10 12 0.0023 0.0839 0.000 510 600 600 1.015 0 1; 11 13 0.0061 0.0476 0.0999 600 625 625 0 0 1; 11 14 0.0054 0.0418 0.0879 600 625 625 0 0 1; 12 13 0.0061 0.0476 0.0999 600 625 625 0 0 1; 12 23 0.0124 0.0966 0.2030 600 625 625 0 0 1; 13 23 0.0111 0.0865 0.1818 600 625 625 0 0 1; 14 16 0.0050 0.0389 0.0818 600 625 625 0 0 1; 15 16 0.0022 0.0173 0.0364 600 625 625 0 0 1; 15 21 0.0063 0.0490 0.1030 600 625 625 0 0 1; 15 21 0.0063 0.0490 0.1030 600 625 625 0 0 1; 15 24 0.0067 0.0519 0.1091 600 625 625 0 0 1; 16 17 0.0033 0.0259 0.0545 600 625 625 0 0 1; 16 19 0.0030 0.0231 0.0485 600 625 625 0 0 1; 17 18 0.0018 0.0144 0.0303 600 625 625 0 0 1; 17 22 0.0135 0.1053 0.2212 600 625 625 0 0 1; 18 21 0.0033 0.0259 0.0545 600 625 625 0 0 1; 18 21 0.0033 0.0259 0.0545 600 625 625 0 0 1; 19 20 0.0051 0.0396 0.0833 600 625 625 0 0 1; 19 20 0.0051 0.0396 0.0833 600 625 625 0 0 1; 20 23 0.0028 0.0216 0.0455 600 625 625 0 0 1; 20 23 0.0028 0.0216 0.0455 600 625 625 0 0 1; 21 22 0.0087 0.0678 0.1424 600 625 625 0 0 1; ]; %% area data areas = [ 1 1; 12 13; ]; %% generator cost data % 1 startup shutdown n x1 y1 ... xn yn % 2 startup shutdown n c(n-1) ... c0 gencost = [ 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; 2 1 0 3 1 1 1; ];