目标函数:给定恢复时间内发电量最大
约束条件:稳态频率约束+惯量响应约束
风电场内有27台风机,共分为5组。
每个操作时步有两种可能:辅机启动、风机启动;
每组辅机或风机一起启动,且某组辅机启动后下一时步风机启动。
第一时步预估后,设置启动功率。
warning('off');
clear
clc
yalmip;
%新能源场站内节点上的辅机及新能源机组信息
%第1列:节点号 第2列:辅机功率 第3列:额定容量 第4列:功频静特性系数(机组) 第5列:功频静特性系数(辅机) 第6列:惯性时间常数
Ebus=[
1,18000,2000000,25,1,10;
2,18000,2000000,25,1,10;
3,18000,2000000,25,1,10;
4,18000,2000000,25,1,10;
5,18000,2000000,25,1,10;
6,18000,2000000,25,1,10;
7,18000,2000000,25,1,10;
8,18000,2000000,25,1,10;
9,18000,2000000,25,1,10;
10,18000,2000000,25,1,10;
11,18000,2000000,25,1,10;
12,18000,2000000,25,1,10;
13,18000,2000000,25,1,10;
14,18000,2000000,25,1,10;
15,18000,2000000,25,1,10;
16,18000,2000000,25,1,10;
17,18000,2000000,25,1,10;
18,18000,2000000,25,1,10;
19,18000,2000000,25,1,10;
20,18000,2000000,25,1,10;
21,18000,2000000,25,1,10;
22,18000,2000000,25,1,10;
23,18000,2000000,25,1,10;
24,18000,2000000,25,1,10;
25,18000,2000000,25,1,10;
26,18000,2000000,25,1,10;
27,18000,2000000,25,1,10;
];
ND=length(Ebus);%有27个节点
%新能源场站内线路信息
%第1列:线路号 第2列:首节点 第3列:末节点 第4列:rij 第5列:xij
Ebranch=[1,0,1;
2,1,2;
3,2,3;
4,3,4;
5,4,5;
6,0,6;
7,6,7;
8,7,8;
9,8,9;
10,9,10;
11,0,11;
12,11,12;
13,12,13;
14,13,14;
15,14,15;
16,0,16;
17,16,17;
18,17,18;
19,18,19;
20,19,20;
21,20,21;
22,0,22;
23,22,23;
24,23,24;
25,24,25;
26,25,26;
27,26,27;
];
NC=length(find(0==Ebranch(:,2)));%有5组
%% 额定常数
NT=2*NC;%最大总时步为2倍分组数
fmax=0.2/50;%最大频率偏差;0.2/50P.U.
R=0.5/50;%频率变化量最大值 0.5HZ/s
SG=2500000;%外部带电域或储能的额定容量 2.5MW
KG=50;%外部带电域或储能的功频静特性系数 50p.u.
HG=5;%外部带电域或储能的惯性时间常数 5s
PL=Ebus(:,2);
SW=Ebus(:,3);
KW=Ebus(:,4);
KL=Ebus(:,5);
HW=Ebus(:,6);
%% 决策变量
K&