-
Notifications
You must be signed in to change notification settings - Fork 13
/
Copy pathdataaustralian14gen.m
351 lines (327 loc) · 8.46 KB
/
dataaustralian14gen.m
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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
PSSMODEL=1;
%bus p(MW) q (MVar)
load_case4=[...
102 270 30
205 235 25
206 80 10
207 1130 120
208 125 15
211 1060 110
212 1000 110
215 290 30
216 1105 120
217 750 80
306 900 90
307 470 50
308 620 100
309 140 15
312 92 10
313 1625 165
314 180 20
405 730 75
406 540 55
407 0 0
408 110 10
409 190 20
410 390 40
411 420 45
412 922 100
504 180 20
507 640 65
508 490 50
509 122 15];
load_case4(20,2)=1;
load_mag=load_case4;
%bus MW(each unit) KV H Xa Xd Xq Xd' Tdo' Xd" Tdo" Xq' Tqo' Xq" Tqo" r
%1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
G_data=[...
101 333.3 15 3.6 0.14 1.1 0.65 0.25 8.5 0.25 0.05 0 0 0.25 0.2 0
201 666.7 20 3.2 0.2 1.8 1.75 0.3 8.5 0.21 0.04 0.7 2 0.21 0.25 0
202 555.6 20 2.8 0.17 2.2 2.1 0.3 4.5 0.2 0.04 0.5 1.5 0.21 0.06 0
203 555.6 20 2.6 0.2 2.3 1.7 0.3 5 0.25 0.03 0.4 2 0.25 0.25 0
204 666.7 20 3.2 0.2 1.8 1.75 0.3 8.5 0.21 0.04 0.7 2 0.21 0.25 0
301 666.7 20 2.8 0.2 2.7 1.5 0.3 7.5 0.25 0.04 0.85 0.85 0.25 0.12 0
302 444.4 20 3.5 0.15 2 1.8 0.25 7.5 0.2 0.04 0 0 0.2 0.25 0
401 444.4 20 2.6 0.2 2.3 1.7 0.3 5 0.25 0.03 0.4 2 0.25 0.25 0
402 333.3 20 3 0.2 1.9 1.8 0.3 6.5 0.26 0.035 0.55 1.4 0.26 0.04 0
403 444.4 20 2.6 0.2 2.3 1.7 0.3 5 0.25 0.03 0.4 2 0.25 0.25 0
404 333.3 20 4 0.18 2.2 1.4 0.32 9 0.24 0.04 0.75 1.4 0.24 0.13 0
501 333.3 20 3.5 0.15 2.2 1.7 0.3 7.5 0.24 0.025 0.8 1.5 0.24 0.1 0
502 250 15 4 0.2 2 1.5 0.3 7.5 0.22 0.04 0.8 3 0.22 0.2 0
503 166.7 15 7.5 0.15 2.3 2 0.25 5 0.17 0.022 0.35 1 0.17 0.035 0];
G_data(1,13)=1;
G_data(7,13)=1;
active_gen_case4=[2 4 3 3 4 6 2 3 2 3 3 2 3 1];
G_data(:,2)=G_data(:,2).*active_gen_case4';
gentestcase=[2:8,11,12,14];
genrem=setdiff(1:14,gentestcase);
Hinv=1/2*diag(G_data(gentestcase,4))^-1;
% 1 2 3 4 5 6 7 8 9 10 11 12
% bus Tr KA TA TB TC KE TE KF TF TB1(s) TC1(s)
Exc=[...
101 0 200 0.1 13.25 2.5 0 0 0 0 0 0
%201 0 300 0.01 0.7 0.35 0 0 0 0 0 0
201 0 400 0.02 1.12 0.5 0 0 0 0 0 0
202 0 400 0.02 0 0 1 1 0.029 1 0 0
203 0 300 0.01 0.7 0.35 0 0 0 0 0 0
%204 0 300 0.01 0.7 0.35 0 0 0 0 0 0
204 0 400 0.02 1.12 0.5 0 0 0 0 0 0
301 0 400 0.05 6.42 1.14 0 0 0 0 0 0
302 0 200 0.05 0 0 1 1.333 0.02 0.8 0 0
401 0 300 0.1 40 4 0 0 0 0 0 0
402 0.02 300 0.05 9.8 1.52 0 0 0 0 0 0
403 0 300 0.01 0.7 0.35 0 0 0 0 0 0
404 0 250 0.2 0.0232 0.136 0 0 0 0 0 0
501 0 1000 0.04 0 0 1 0.87 0.004 0.27 0 0
502 0 400 0.5 16 1.4 0 0 0 0 0.05 0.6
503 0 300 0.01 0.8 0.2 0 0 0 0 0 0
];
Exc(:,2)=0.02;
Exc(13,9)=.001;
Exc(13,10)=1;
% bus rating pg (MW) Heavy load case
% bus MVA(per unit) p(per unit) q(per unit) active units
p_case4=[...
101 333.3 .001 -97.4 2
201 666.7 540 -30.8 4
202 555.6 460 -2.5 3
203 555.6 470 9.4 3
204 666.7 399.3 -43.6 4
301 666.7 555 16.6 6
302 444.4 380 -9.3 2
401 444.4 320 -21.9 3
402 333.3 290 -2.4 2
403 444.4 320 14.2 3
404 333.3 217 -3.5 3
501 333.3 280 -52.5 2
502 250 180 -1.8 3
503 166.7 150 2.2 1];
p=p_case4;
p(:,3)=p(:,3).*p(:,5);
p(:,4)=p(:,4).*p(:,5);
pref=p(:,3)./(p(:,2).*p(:,5));
bus_v=[...
101 15
102 330
201 20
202 20
203 20
204 20
205 330
206 330
207 330
208 330
209 330
210 500
211 330
212 330
213 500
214 330
215 330
216 330
217 330
301 20
302 20
303 500
304 500
305 500
306 500
307 500
308 500
309 330
310 330
311 330
312 220
313 220
314 220
315 275
401 20
402 20
403 20
404 20
405 275
406 275
407 275
408 275
409 275
410 275
411 275
412 275
413 275
414 330
415 330
416 330
501 20
502 15
503 15
504 275
505 275
506 275
507 275
508 275
509 275];
% No from To v1(kv) v2(kv) rating x%(on rating) total number of units
%1 2 3 4 5 6 7 8
trans=[...
1 101 102 15 330 333.3 12 12
2 201 206 20 330 666.7 16 6
3 202 209 20 330 555.6 16 5
4 203 208 20 330 555.6 17 4
5 204 215 20 330 666.7 16 6
6 209 210 330 500 625 17 4
7 213 214 500 330 625 17 4
8 301 303 20 500 666.7 16 8
9 302 312 20 220 444.4 15 4
10 304 313 500 220 500 16 2
11 305 311 500 330 500 12 2
12 305 314 500 220 700 17 2
13 308 315 500 275 370 10 2
14 401 410 20 275 444.4 15 4
15 402 408 20 275 333.3 17 3
16 403 407 20 275 444.4 15 4
17 404 405 20 275 333.3 17 6
18 413 414 275 330 750 6 3
19 501 504 20 275 333.3 17 2
20 502 505 15 275 166.7 16 4
21 503 506 15 275 250 16.7 6];
active_trans_case4=[2,4,3,3,4,4,4,6,2,2,2,2,2,3,2,3,3,3,2,3,1];
trans(:,6)=trans(:,6).*active_trans_case4';
trans(:,7)=trans(:,7)/100;
%from to r X b(pu) v(kv) No
line=[...
102 217 0.0084 0.0667 0.817 330 1
102 217 0.0078 0.062 0.76 330 2
102 309 0.0045 0.0356 0.437 330 3
102 309 0.0109 0.0868 0.76 330 4
205 206 0.0096 0.076 0.931 330 5
205 416 0.0037 0.046 0.73 330 6
206 207 0.0045 0.0356 0.437 330 7
206 212 0.0066 0.0527 0.646 330 8
206 215 0.0066 0.0527 0.646 330 9
207 208 0.0018 0.014 0.171 330 10
207 209 0.0008 0.0062 0.076 330 11
208 211 0.0031 0.0248 0.304 330 12
209 212 0.0045 0.0356 0.437 330 13
210 213 0.001 0.0145 1.54 500 14
211 212 0.0014 0.0108 0.133 330 15
211 214 0.0019 0.0155 0.19 330 16
212 217 0.007 0.0558 0.684 330 17
214 216 0.001 0.0077 0.095 330 18
214 217 0.0049 0.0388 0.475 330 19
215 216 0.0051 0.0403 0.494 330 20
215 217 0.0072 0.0574 0.703 330 21
216 217 0.0051 0.0403 0.494 330 22
303 304 0.001 0.014 1.48 500 23
303 305 0.0011 0.016 1.7 500 24
304 305 0.0003 0.004 0.424 500 25
305 306 0.0002 0.003 0.32 500 26
305 307 0.0003 0.0045 0.447 500 27
306 307 0.0001 0.0012 0.127 500 28
307 308 0.0023 0.0325 3.445 500 29
309 310 0.009 0.0713 0.874 330 30
310 311 0 -0.0337 0 330 31
312 313 0.002 0.015 0.9 220 32
313 314 0.0005 0.005 0.52 220 33
315 509 0.007 0.05 0.19 275 34
405 406 0.0039 0.0475 0.381 275 35
405 408 0.0054 0.05 0.189 275 36
405 409 0.018 0.122 0.79 275 37
406 407 0.0006 0.0076 0.062 275 38
407 408 0.0042 0.0513 0.412 275 39
408 410 0.011 0.128 1.01 275 40
409 411 0.0103 0.0709 0.46 275 41
410 411 0.0043 0.0532 0.427 275 42
410 412 0.0043 0.0532 0.427 275 43
410 413 0.004 0.0494 0.4 275 44
411 412 0.0012 0.0152 0.122 275 45
414 415 0.002 0.025 0.39 330 46
415 416 0.0037 0.046 0.73 330 47
504 507 0.023 0.15 0.56 275 48
504 508 0.026 0.019 0.87 275 49
505 507 0.0008 0.0085 0.06 275 50
505 508 0.0025 0.028 0.17 275 51
506 507 0.0008 0.0085 0.06 275 52
506 508 0.003 0.028 0.14 275 53
507 508 0.002 0.019 0.09 275 54
507 509 0.03 0.22 0.9 275 55
];
zbase=(line(:,6).^2)/100;
line(:,3)=line(:,3).*zbase;
line(:,4)=line(:,4).*zbase/(2*pi*50);
line(:,5)=line(:,5)./zbase/(2*pi*50);
% No bus QL Qc
comp=[...
1 211 0 0
2 212 0 400
3 216 0 300
4 409 0 60
5 411 0 30
6 414 30 0
7 415 60 0
8 416 60 0
9 504 0 0
];
SVC=[...
1 205 1.045 -39.3
2 313 1.015 86.7
3 412 1 -52.2
4 507 1.01 -4
5 509 1.03 -109.3];
PSSMODEL=1;
G=1; FL=.1; KL=10; FI=1; KI=10; FH=13; KH=120;
Tr=.001;
Ka=200; Ta=.001;
Ke=1; Te=0;
Tb=0; Tc=0;
Kf=0; Tf=0;
Efmin=-5; Efmax=5;
alpha(1:15)=0*1e11*ones(15,1);
alpha(16:30)=0*1e8*ones(15,1);
nred=90;
%PLL parameters
PLL_Ki=2/500;
PLL_Kp=2/500;
PLL_tau=1/500;
visat=inf;%(40+2.9)*1e6;
Isel=eye(2*15);
Isel_forming=eye(3*15);
thetavsc0(1)=6.6166;
thetavsc0(2)=21.2703;
thetavsc0(3)=15.5222;
thetavsc0(4)=28.8703;
thetavsc0(5)=14.0249;
thetavsc0(6)=-6.3379;
thetavsc0(7)=4.7928;
thetavsc0(8)=3.6837;
thetavsc0(9)=-1.0986;
thetavsc0(10)=20.7;
thetavsc0(11)=10.5944;
thetavsc0(12)=1.8335;
thetavsc0(13)=-12.7055;
thetavsc0(14)=-12.4459;
thetavsc0(15)=-14.7643;
Kvsc(1,:)=1*ones(1,15);
Kvsc(2,:)=1*ones(1,15);
dist.inp=zeros(6,1);
dist.sc=ones(6,1);
dist.time=50;
c_fil=1e-5; % 1e-5
r_fil=5e-1;
l_fil=5e-3;
% Grid forming
Pinit=zeros(15,1);
Vbase(1)=15e3;
Vbase(2)=20e3;
Vbase(3)=20e3;
Vbase(4)=20e3;
Vbase(5)=20e3;
Vbase(6)=20e3;
Vbase(7)=15e3;
Vbase(8)=20e3;
Vbase(9)=20e3;
Vbase(10)=20e3;
Vbase(11)=20e3;
Vbase(12)=20e3;
Vbase(13)=20e3;
Vbase(14)=20e3;
Vbase(15)=15e3;