Topology and data of modified Barry Island test case.xls

	A	B	C	D	E	F	G	H	I	J	K	L	M	N	O	P	Q	R	S	T	U	V
1	%%----- Topo Data -----%%
2	%% system MVA base
3	mpc.baseMVA = 100;
4
5	%% bus data
6	%	bus_i	type	Pd	Qd	Gs	Bs	area	Vm	Va	baseKV	zone	Vmax	Vmin
7	mpc.bus = [ %% (Pd and Qd are specified in kW & kVAr here, and this value is the original value for reference, and load power is shown below)
8		1	3	0	0	0	0	1	1	0	12.66	1	1	1;
9		2	1	100	60	0	0	1	1	0	12.66	1	1.1	0.9;
10		3	1	90	40	0	0	1	1	0	12.66	1	1.1	0.9;
11		4	1	120	80	0	0	1	1	0	12.66	1	1.1	0.9;
12		5	1	60	30	0	0	1	1	0	12.66	1	1.1	0.9;
13		6	1	60	20	0	0	1	1	0	12.66	1	1.1	0.9;
14		7	1	200	100	0	0	1	1	0	12.66	1	1.1	0.9;
15		8	1	200	100	0	0	1	1	0	12.66	1	1.1	0.9;
16		9	1	60	20	0	0	1	1	0	12.66	1	1.1	0.9;
17		10	1	60	20	0	0	1	1	0	12.66	1	1.1	0.9;
18		11	1	45	30	0	0	1	1	0	12.66	1	1.1	0.9;
19		12	1	60	35	0	0	1	1	0	12.66	1	1.1	0.9;
20		13	1	60	35	0	0	1	1	0	12.66	1	1.1	0.9;
21		14	1	120	80	0	0	1	1	0	12.66	1	1.1	0.9;
22		15	1	60	10	0	0	1	1	0	12.66	1	1.1	0.9;
23		16	1	60	20	0	0	1	1	0	12.66	1	1.1	0.9;
24		17	1	60	20	0	0	1	1	0	12.66	1	1.1	0.9;
25		18	2	90	40	0	0	1	1	0	12.66	1	1.1	0.9;
26		19	1	90	40	0	0	1	1	0	12.66	1	1.1	0.9;
27		20	1	90	40	0	0	1	1	0	12.66	1	1.1	0.9;
28		21	1	90	40	0	0	1	1	0	12.66	1	1.1	0.9;
29		22	1	90	40	0	0	1	1	0	12.66	1	1.1	0.9;
30		23	1	90	50	0	0	1	1	0	12.66	1	1.1	0.9;
31		24	1	420	200	0	0	1	1	0	12.66	1	1.1	0.9;
32		25	1	420	200	0	0	1	1	0	12.66	1	1.1	0.9;
33		26	1	60	25	0	0	1	1	0	12.66	1	1.1	0.9;
34		27	1	60	25	0	0	1	1	0	12.66	1	1.1	0.9;
35		28	1	60	20	0	0	1	1	0	12.66	1	1.1	0.9;
36		29	1	120	70	0	0	1	1	0	12.66	1	1.1	0.9;
37		30	1	200	600	0	0	1	1	0	12.66	1	1.1	0.9;
38		31	1	150	70	0	0	1	1	0	12.66	1	1.1	0.9;
39		32	1	210	100	0	0	1	1	0	12.66	1	1.1	0.9;
40		33	2	60	40	0	0	1	1	0	12.66	1	1.1	0.9;
41	];
42
43	%% generator data(here, status=1 means this gen is a CHP unit, status=2 means this gen is a thermal unit)
44	%	bus	Pg	Qg	Qmax	Qmin	Vg	mBase	status	Pmax	Pmin	Pc1	Pc2	Qc1min	Qc1max	Qc2min	Qc2max	ramp_agc	ramp_10	ramp_30	ramp_q	apf
45	mpc.gen = [
46		1	0	0	100	-100	1.05	100	1	50	0	0	0	0	0	0	0	0	0	0	0	0;
47		18	0	0	100	-100	1.03	100	1	50	0	0	0	0	0	0	0	0	0	0	0	0;
48		33	0	0	100	-100	1.03	100	1	50	0	0	0	0	0	0	0	0	0	0	0	0;
49	];
50
51	%% branch data
52	%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status	angmin	angmax
53	mpc.branch = [ %% (r and x specified in ohms here)
54		1	2	0.0922	0.047	0	0	0	0	0	0	1	-360	360;
55		2	3	0.493	0.2511	0	0	0	0	0	0	1	-360	360;
56		3	4	0.366	0.1864	0	0	0	0	0	0	1	-360	360;
57		4	5	0.3811	0.1941	0	0	0	0	0	0	1	-360	360;
58		5	6	0.819	0.707	0	0	0	0	0	0	1	-360	360;
59		6	7	0.1872	0.6188	0	0	0	0	0	0	1	-360	360;
60		7	8	0.7114	0.2351	0	0	0	0	0	0	1	-360	360;
61		8	9	1.03	0.74	0	0	0	0	0	0	1	-360	360;
62		9	10	1.044	0.74	0	0	0	0	0	0	1	-360	360;
63		10	11	0.1966	0.065	0	0	0	0	0	0	1	-360	360;
64		11	12	0.3744	0.1238	0	0	0	0	0	0	1	-360	360;
65		12	13	1.468	1.155	0	0	0	0	0	0	1	-360	360;
66		13	14	0.5416	0.7129	0	0	0	0	0	0	1	-360	360;
67		14	15	0.591	0.526	0	0	0	0	0	0	1	-360	360;
68		15	16	0.7463	0.545	0	0	0	0	0	0	1	-360	360;
69		16	17	1.289	1.721	0	0	0	0	0	0	1	-360	360;
70		17	18	0.732	0.574	0	0	0	0	0	0	1	-360	360;
71		2	19	0.164	0.1565	0	0	0	0	0	0	1	-360	360;
72		19	20	1.5042	1.3554	0	0	0	0	0	0	1	-360	360;
73		20	21	0.4095	0.4784	0	0	0	0	0	0	1	-360	360;
74		21	22	0.7089	0.9373	0	0	0	0	0	0	1	-360	360;
75		3	23	0.4512	0.3083	0	0	0	0	0	0	1	-360	360;
76		23	24	0.898	0.7091	0	0	0	0	0	0	1	-360	360;
77		24	25	0.896	0.7011	0	0	0	0	0	0	1	-360	360;
78		6	26	0.203	0.1034	0	0	0	0	0	0	1	-360	360;
79		26	27	0.2842	0.1447	0	0	0	0	0	0	1	-360	360;
80		27	28	1.059	0.9337	0	0	0	0	0	0	1	-360	360;
81		28	29	0.8042	0.7006	0	0	0	0	0	0	1	-360	360;
82		29	30	0.5075	0.2585	0	0	0	0	0	0	1	-360	360;
83		30	31	0.9744	0.963	0	0	0	0	0	0	1	-360	360;
84		31	32	0.3105	0.3619	0	0	0	0	0	0	1	-360	360;
85		32	33	0.341	0.5302	0	0	0	0	0	0	1	-360	360;
86		21	8	2	2	0	0	0	0	0	0	0	-360	360;
87		9	15	2	2	0	0	0	0	0	0	0	-360	360;
88		12	22	2	2	0	0	0	0	0	0	0	-360	360;
89		18	33	0.5	0.5	0	0	0	0	0	0	0	-360	360;
90		25	29	0.5	0.5	0	0	0	0	0	0	0	-360	360;
91	];
92
93	%%----- OPF Data -----%%
94	%% generator cost data (only thermal unit)
95	%	1	startup	shutdown	n	x1	y1	...	xn	yn
96	%	2	startup	shutdown	n	c(n-1)	...	c0
97	mpc.gencost = [
98		18	0	0	0	38	0.172	0
99		33	0	0	0	36	0.186	0
100	];