机组组合问题Word文件下载.docx
《机组组合问题Word文件下载.docx》由会员分享,可在线阅读,更多相关《机组组合问题Word文件下载.docx(24页珍藏版)》请在冰豆网上搜索。
![机组组合问题Word文件下载.docx](https://file1.bdocx.com/fileroot1/2023-1/25/4811112a-a45f-4975-abc1-a9bb74194063/4811112a-a45f-4975-abc1-a9bb741940631.gif)
对于问题2,在问题1的基础上进行了进一步的讨论,还需要满足机组的稳定运行出力约束,启动停运时的处理约束、最小停运及运行时间约束等,目标函数与问题1相同,仍然是使总费用最小,同样可建立优化模型进行求解。
对于问题3,用IEEE—118节点的电力系统对问题2的求解模型进行测试,由于该问题属于大规模机组组合问题,可使用遗传算法进行求解,从而得出最优机组组合计划。
三、模型的假设
1)一个小时之内发电机发电出力和电力负荷均保持不变。
2)输电线路不存在电能的损耗。
3)3母线系统空载成本和增量成本之和随发电出力呈分段线性增长。
4)除题目所给的约束外,不受其它无关因素的影响。
四、符号的约定
机组i在t时段的发电出力
t时段的系统备用要求
机组i在t时段的运行状态,
机组i在t时段的空载成本和增量成本之和
机组i的启动成本
负荷j在t时段的负荷量
机组i的最大出力
机组i的最小稳定运行出力
t时段母线
上的注入功率
t时段第k根输电线路上流过的电能
第k根输电线路最大传输容量
第k根输电线路第l条母线的线性转移因子
机组i的最大增出力
机组i的最大减出力
Y
总成本
机组i的最小运行时间
机组i的最小停运时间
m
负荷数
n
机组台数
T
所研究的计划周期
N
母线数
M
输电线路数
i=1,2,…,n;
j=1,2,…,m;
k=1,2,…,M;
l=0,1,…,N;
t=1,2,…,T
五、模型的建立与求解
5.13母线系统Ⅰ模型的建立及求解
5.1.1成本分析
机组组合是一个最优化问题,即在一定的约束条件下,确定某些可控制的合理取值,使选定的目标达到最优的问题[1]。
发电机组的成本包括启动成本、空载成本和增量成本三部分组成。
1)启动成本:
当机组从停运状态(不发电)变为运行状态(发电)时发生的成本。
2)空载成本:
只要机组处于运行状态,就会发生的成本。
3)增量成本:
与该机组发电量有关的成本。
发电机组的发电成本如下表1所示
表1机组成本(所有小时都相同)
机组
启动成本(¥)
空载成本(¥)
增量发电量(MWh)
增量成本(¥/MWh)
G1
350
100
100,100
10,14
G2
200
60,40
12,15
则该发电机组的空载成本和增量成本如图2所示。
G1:
当该机组的发电出力在[0,100]MW时范围时,每MWh的成本为10¥/MWh;
当该机组的发电出力在[100,200]MW时范围时,每MWh的成本为14¥/MWh。
G2:
当该机组的发电出力在[0,60]MW时范围时,每MWh的成本为12¥/MWh;
当该机组的发电出力在[60,100]MW时范围时,每MWh的成本为15¥/MWh。
图2空载成本和增量成本之和随该机组发电出力增长走势图
由上图可知,机组G1和机组G2在t(t=1,2,3,4)时段的的成本特性可表示为G1:
(1)
(2)
若用
表示机组i在t时段的运行状态,其中
表示运行,
表示停止,则
表示机组i在t-1时段的运行状态。
所以启动成本可表示为
(3)
其中
为机组i的启动成本。
由题意知,
总成本=启动成本+空载成本+增量成本
设所研究的计划周期为
,机组台数为
,则总成本Y可表示为
(4)
5.1.2约束条件分析
1)负载平衡约束:
任何小时,电力负荷之和必须等于发电机发出电力出力之和。
即
(5)
其中,
为负荷j在t时段的负荷量,m为负荷数。
该问题负荷L见下表2
表2Bus2负荷L
小时
1
2
3
4
L(MW)
130
170
140
2)系统备用约束:
发电机的备用容量为处于运行状态的发电机的最大发电能力减去其出力,其中处于停运状态的发电机的备用容量为0。
任何小时,发电机的容量之和必须大于系统备用要求。
即在t时段,
(6)
为机组i的最大出力,
为t时段的系统备用要求,如表3
表3系统备用要求
备用(MW)
20
30
50
40
3)输电线路传输容量约束:
线路传输的电能必须在它的传输容量范围内。
线路k上流过的电能为
,(7)
为t时段母线
(busl)上的注入功率,
为第k根输电线路第l条母线的线性转移因子,
为t时段第k根输电线路上流过的电能,且k=1表示线路0-1,k=2表示线路0-2,k=3表示线路1-2。
线性转移因子如表4。
表4线性转移因子表
bus0
bus1
bus2
line-1
-0.6667
-0.3333
line-2
line-3
0.3333
由于线路传输的电能不能大于其最大传输容量,又知线路上的电流是有方向的,故
(8)
所以
(9)
为第k根输电线路最大传输容量,见表5。
表5输电线路最大传输容量(所有小时都相同)
线路
线路容量(MW)
4)发电机组出力范围约束:
处于运行状态的发电机组i的发电出力必须小于其最大发电能力。
(10)
5)机组增降出力约束:
发电机组在增加发电出力时,增加出力的速度要小
于其最大增出力;
发电机组在减少发电出力时,减少出力的速度要小于其最大减出力。
(11)
为机组i的最大增出力,
为机组i的最大减出力。
机组G1和机组G2的物理特性及初始状态分别见表6和表7。
表6机组物理特性
所在母线
最大出力(MW)
最大增出力(MW/h)
最大减出力(MW/h)
Bus0
Bus1
60
表7机组初始状态
状态小时数
兆瓦(MW)
运行2小时
关机3小时
5.1.3模型的建立及求解
机组组合的目的是针对在指定的周期内,满足系统负荷、备用容量、发电出力、传输容量、增降出力等限制,优化确定各机组的启停机计划和优化分配其发电负荷,使发电总费用最小。
因此,要以机组的费用最小为依据建立相应的目标函数[3]。
故该3母线系统优化模型可建立如下:
(12)
用Lingo程序求解[3](见附录1),可得结果如表8
表83母线系统Ⅰ的最优机组组合计划
发电机组
状态
运行
停运
发电出力(MW)
备用容量(MW)
70
110
90
总成本(¥)
6680
5.23母线系统Ⅱ模型的建立及求解
在问题1的基础上,考虑发电机组的下列物理特性约束,重新制定最优机组组合计划。
1)发电机组的稳定出力范围约束:
处于运行状态的发电机组i的发电出力必须大于其最小稳定运行出力,即发电机组不能在[0,pimin]出力范围内稳定运行。
故
(13)
为机组i的最小稳定运行出力。
2)机组启动和停运时的出力约束:
当机组从停运状态变为运行状态时,机组在该小时的发电出力必须为其最小稳定运行出力
,即运行状态的第一个小时出力必须为
。
故当
时,
(14)
当机组从运行状态变为停运状态时,机组在该小时的发电出力必须为其最小稳定运行出力
(15)
3)机组最小运行时间约束:
由表1和表7可知,机组G1的启动成本远高于G2,且增量发电量一定时的增量成本又低于G2,又已知机组G1已在运行,故在时间较短时,应避免使G1的停运和再次启动。
从问题一中不难看出,当费用最小时机组G1在各时段的状态均为运行,且其发电出力基本满足出力约束,所以该问题中不必对其运行的时间作限制。
当机组2在t时段停止且在t-1时段运行时
(16)
为机组i的最小运行时间。
否则
(17)
(18)
4)机组最小停运时间约束:
当机组i在t时段运行且在t-1时段停止时
(19)
为机组i的最小停运时间。
(20)
(21)
机组G1和机组G2的物理特性(补充)见表9,
表9机组物理特性(补充)
最小出力(MW)
最小运行时间(h)
最小停运时间(h)
目标函数和其它约束条件同问题1,故该3母线系统优化模型建立如下:
(22)
用Lingo程序求解[3](见附录2),可得结果如表10
表103母线系统Ⅱ的最优机组组合计划
80
6820
5.3IEEE118系统模型的建立及求解
由于该问题中的机组规模较大,可采用二次函数对空载成本和增量成本曲线参数进行拟合,机组i在t时段的曲线方程为
,
分别为曲线方程的二次项系数、一次项系数及常数项。
同问题二,该IEEE118系统的数学模型如下:
(23)
可利用模拟退火算法对该模型进行测试,该算法的流程图如图3
图3遗传算法流程图
六、模型的评价
6.1模型的评价
6.1.1模型的优点
1)该问题对目标函数和约束条件分析透彻,采用的数学模型有成熟的理论基础,可信度较高。
2)利用Lingo编程的方法对模型严格求解,具有科学性。
3)本文建立的模型与实际联系紧密,模型求解简便,可适用于其它类似的问题,通用性、推广性较强。
6.1.2模型的缺点
1)该问题中把空载成本和增量成本之和随发电出力的变化看成线性增长,因此计算误差可能较大;
2)该问题对机组发电出力的安排过于理想化,与现实情况存在一定的差距。
参考文献
[1]姜启源等.数学模型.北京:
高等教育出版社,2003
[2]吴建国.数学建模案例精编.北京:
中国水利水电出版社,2005
[3]谢金星等.优化建模与LINDO/LINGO软件.北京:
清华大学出版社,2005
附录
附录1求解3母线系统Ⅰ的Lingo程序
min=x11*(1-x10)*z1+x12*(1-x11)*z1+x13*(1-x12)*z1+x14*(1-x13)*z1+x21*(1-x20)*z2+x22*(1-x21)*z2+x23*(1-x22)*z2+x24*(1-x23)*z2+x11*y11+x12*y12+x13*y13+x14*y14+x21*y21+x22*y22+x23*y23+x24*y24;
100=x11*p11+x21*p21;
130=x12*p12+x22*p22;
170=x13*p13+x23*p23;
140=x14*p14+x24*p24;
x11*(200-p11)+x21*(100-p21)>
20;
x12*(200-p12)+x21*(100-p22)>
30;
x13*(200-p13)+x23*(100-p23)>
50;
x14*(200-p14)+x24*(100-p24)>
40;
-0.6667*x21*p21+0.3333*100<
=200;
-0.6667*x21*p21+0.3333*100>
=-200;
-0.3333*x21*p21+0.6667*100<
=100;
-0.3333*x21*p21+0.6667*100>
=-100;
0.3333*x21*p21+0.3333*100<
0.3333*x21*p21+0.3333*100>
-0.6667*x22*p22+0.3333*130<
-0.6667*x22*p22+0.3333*130>
-0.3333*x22*p22+0.6667*130<
-0.3333*x22*p22+0.6667*130>
0.3333*x22*p22+0.3333*130<
0.3333*x22*p22+0.3333*130>
-0.6667*x23*p23+0.3333*170<
-0.6667*x23*p23+0.3333*170>
-0.3333*x23*p23+0.6667*170<
-0.3333*x23*p23+0.6667*170>
0.3333*x23*p23+0.3333*170<
0.3333*x23*p23+0.3333*170>
-0.6667*x24*p24+0.3333*140<
-0.6667*x24*p24+0.3333*140>
-0.3333*x24*p24+0.6667*140<
-0.3333*x24*p24+0.6667*140>
0.3333*x24*p24+0.3333*140<
0.3333*x24*p24+0.3333*140>
p11<
200;
p21<
100;
p12<
p22<
p13<
p23<
p14<
p24<
x11*p11-x10*p10>
=-50;
x11*p11-x10*p10<
=30;
x21*p21-x20*p20>
=-60;
x21*p21-x21*p20<
=40;
x12*p12-x11*p11>
x12*p12-x11*p11<
x22*p22-x21*p21>
x22*p22-x21*p21<
x13*p13-x12*p12>
x13*p13-x12*p12<
x23*p23-x22*p22>
x23*p23-x22*p22<
x14*p14-x13*p13>
x14*p14-x13*p13<
x24*p24-x23*p23>
x24*p24-x23*p23<
y11=(10*p11+100)*((p11#ge#(-0.1))#and#(100.1#ge#p11))+(14*p11-300)*
((p11#ge#100.1)#and#(200#ge#p11));
y12=(10*p12+100)*((p12#ge#(-0.1))#and#(100.1#ge#p12))+(14*p12-300)*
((p12#ge#100.1)#and#(200#ge#p12));
y13=(10*p13+100)*((p13#ge#(-0.1))#and#(100.1#ge#p13))+(14*p13-300)*
((p13#ge#100.1)#and#(200#ge#p13));
y14=(10*p14+100)*((p14#ge#(-0.1))#and#(100.1#ge#p14))+(14*p14-300)*
((p14#ge#100.1)#and#(200#ge#p14));
y21=(12*p21+200)*((p21#ge#(-0.1))#and#(60.1#ge#p21))+(15*p21+20)*
((p21#ge#60)#and#(100#ge#p21));
y22=(12*p22+200)*((p22#ge#(-0.1))#and#(60.1#ge#p22))+(15*p22+20)*
((p22#ge#60.1)#and#(100#ge#p22));
y23=(12*p23+200)*((p23#ge#(-0.1))#and#(60.1#ge#p23))+(15*p23+20)*
((p23#ge#60.1)#and#(100#ge#p23));
y24=(12*p24+200)*((p24#ge#(-0.1))#and#(60.1#ge#p24))+(15*p24+20)*
((p24#ge#60.1)#and#(100#ge#p24));
@bin(x11);
@bin(x12);
@bin(x13);
@bin(x14);
@bin(x21);
@bin(x22);
@bin(x23);
@bin(x24);
x10=1;
x20=0;
z1=350;
z2=100;
p10=100;
p20=0;
结果:
Localoptimalsolutionfoundatiteration:
1381
Objectivevalue:
6679.800
VariableValueReducedCost
X111.0000000.000000
X101.0000000.000000
Z1350.00000.000000
X121.0000000.000000
X131.0000000.000000
X141.0000000.000000
X210.0000000.000000
X200.0000000.000000
Z2100.00000.000000
X221.0000000.000000
X231.0000000.000000
X240.0000000.000000
Y111100.0000.000000
Y121100.9990.000000
Y131240.0000.000000
Y141660.0000.000000
Y211520.0000.000000
Y22558.80060.000000
Y23920.00000.000000
Y24204.60550.000000
P11100.00000.000000
P21100.00000.000000
P12100.09990.000000
P2229.900050.000000
P13110.00000.000000
P2360.000000.000000
P14140.00000.000000
P240.38378930