发布时间 : 星期六 文章5魏至诚、陈奔、卢洁玲：电力市场输电阻塞管理的数学模型--张传林更新完毕开始阅读
电力市场输电阻塞管理的数学模型
h1=@if(a#le#70,(-505*a),@if(a#le#120,124*a,@if(a#le#150,252*a,@if(a#le#190,489*a,489*a))));
h2=@if(b#le#30,(-560*b),@if(b#le#50,182*b,@if(b#le#58,203*b,@if(b#le#73,245*b,@if(b#le#79,300*b,@if(b#le#81,320*b,@if(b#le#89,495*b,495*b)))))));
h3=@if(c#le#110,(-610*c),@if(c#le#150,152*c,@if(c#le#180,233*c,@if(c#le#200,308*c,@if(c#le#240,356*c,@if(c#le#280,500*c,500*c))))));
h4=@if(d#le#55,(-500*d),@if(d#le#60,150*d,@if(d#le#70,170*d,@if(d#le#80,200*d,@if(d#le#90,255*d,@if(d#le#100,302*d,@if(d#le#115,325*d,@if(d#le#116,800*d,800*d))))))));
h5=@if(e#le#75,(-590*e),@if(e#le#80,0*e,@if(e#le#95,116*e,@if(e#le#110,188*e,@if(e#le#125,215*e,@if(e#le#135,310*e,@if(e#le#145,396*e,@if(e#le#155,510*e,510*e))))))));
h6=@if(f#le#95,(-607*f),@if(f#le#105,159*f,@if(f#le#125,173*f,@if(f#le#140,252*f,@if(f#le#150,305*f,@if(f#le#170,380*f,@if(f#le#180,520*f,520*f)))))));
h7=@if(g#le#50,(-500*g),@if(g#le#65,120*g,@if(g#le#70,180*g,@if(g#le#85,251*g,@if(g#le#95,260*g,@if(g#le#105,306*g,@if(g#le#110,315*g,@if(g#le#120,335*g,@if(g#le#123,348*g,@if(g#le#125,548*g,548*g))))))))));
h8=@if(h#le#70,(-800*h),@if(h#le#90,183*h,@if(h#le#110,253*h,@if(h#le#130,303*h,@if(h#le#140,318*h,@if(h#le#155,400*h,@if(h#le#160,800*h,800*h))))))); !各机组每小时的合理收入;
b1=@if(a#le#150,(303*150-p1)*0.4,(h1-p1)*0.6); b2=@if(b#le#79,(303*79-p2)*0.4,(h2-p2)*0.6); b3=@if(c#le#180,(303*180-p3)*0.4,(h3-p3)*0.6); b4=@if(d#le#99.5,(99.5*303-p4)*0.4,(h4-p4)*0.6); b5=@if(e#le#125,(303*125-p5)*0.4,(h5-p5)*0.6); b6=@if(f#le#140,(303*140-p6)*0.4,(h6-p6)*0.6); b7=@if(g#le#95,(303*95-p7)*0.4,(h7-p7)*0.6); b8=@if(h#le#113.9,(303*113.9-p8)*0.4,(h8-p8)*0.6);
!考虑补偿系数后各机组每小时的阻塞费用。对于序内容量，阻塞费用=0.4*（预案收入-实际收入）；
对于序外容量，阻塞费用=0.6*（合理收入-实际收入）;
0.0826*a+0.0478*b+0.0528*c+0.1199*d-0.0257*e+0.1216*f+0.122*g-0.0015*h+110.4775<=165;
-0.0547*a+0.1275*b-0.0001*c+0.0332*d+0.0867*e-0.1127*f-0.0186*g+0.0985*h+131.3521<=150;
-108.993-0.0694*a+0.062*b-0.1565*c-0.0099*d+0.1247*e+0.0024*f-0.0028*g-0.2012*h>=-160;
77.6116-0.0346*a-0.1028*b+0.205*c-0.0209*d-0.012*e+0.0057*f+0.1452*g+0.0763*h<=155; 133.1334+0.0003*a+0.2428*b-0.0647*c-0.0412*d-0.0655*e+0.07*f-0.0039*g-0.0092*h<=132; 120.8481+0.2376*a-0.0607*b-0.0781*c+0.0929*d+0.0466*e-0.0003*f+0.1664*g+0.0004*h<=162;
!保证各线路在安全裕度内工作;
a+b+c+d+e+f+g+h=982.4;!各机组出力之和为982.4MW;
21
电力市场输电阻塞管理的数学模型
a<=153; 87<=a; b<=88; 58<=b; 132<=c; c<=228; 60.5<=d; d<=99.5; 98<=e; e<=152; 95<=f; f<=155; 60.1<=g; g<=102.1; 63<=h; h<=117;
!现阶段各机组向上或向下爬坡后可达到的出力的范围; END
? 附录四：在安全裕度内调整的lingo实现程序
MODEL: min=k;
!是该最大值尽可能小;
k=@smax(y1,y2,y3,y4,y5,y6);
!求出各线路实际潮流值与其相应潮流限值的比例的最大值;
y1=@abs(0.0826*a+0.0478*b+0.0528*c+0.1199*d-0.0257*e+0.1216*f+0.122*g-0.0015*h+110.4775)/165;
y2=@abs(-0.0547*a+0.1275*b-0.0001*c+0.0332*d+0.0867*e-0.1127*f-0.0186*g+0.0985*h+131.3521)/150;
y3=@abs(-108.993-0.0694*a+0.062*b-0.1565*c-0.0099*d+0.1247*e+0.0024*f-0.0028*g-0.2012*h)/160;
y4=@abs(77.6116-0.0346*a-0.1028*b+0.205*c-0.0209*d-0.012*e+0.0057*f+0.1452*g+0.0763*h)/155;
y5=@abs(133.1334+0.0003*a+0.2428*b-0.0647*c-0.0412*d-0.0655*e+0.07*f-0.0039*g-0.0092*h)/132;
y6=@abs(120.8481+0.2376*a-0.0607*b-0.0781*c+0.0929*d+0.0466*e-0.0003*f+0.1664*g+0.0004*h)/162;
!求出各线路实际潮流值与其相应潮流限值的比例的最大值;
22
电力市场输电阻塞管理的数学模型
@abs(0.0826*a+0.0478*b+0.0528*c+0.1199*d-0.0257*e+0.1216*f+0.122*g-0.0015*h+110.4775)<=186.45;
@abs(-0.0547*a+0.1275*b-0.0001*c+0.0332*d+0.0867*e-0.1127*f-0.0186*g+0.0985*h+131.3521)<=177;
@abs(-108.993-0.0694*a+0.062*b-0.1565*c-0.0099*d+0.1247*e+0.0024*f-0.0028*g-0.2012*h)<=174.4;
@abs(77.6116-0.0346*a-0.1028*b+0.205*c-0.0209*d-0.012*e+0.0057*f+0.1452*g+0.0763*h)<=172.05;
@abs(133.1334+0.0003*a+0.2428*b-0.0647*c-0.0412*d-0.0655*e+0.07*f-0.0039*g-0.0092*h)<=151.8;
@abs(120.8481+0.2376*a-0.0607*b-0.0781*c+0.0929*d+0.0466*e-0.0003*f+0.1664*g+0.0004*h)<=184.68;
!保证各线路在安全裕度内工作;
a+b+c+d+e+f+g+h=1052.8; !各机组出力之和为1052.8MW; a<=153; 87<=a; b<=88; 58<=b; 132<=c; c<=228; 60.5<=d; d<=99.5; 98<=e; e<=152; 95<=f; f<=155; 60.1<=g; g<=102.1; 63<=h; h<=117;
!现阶段各机组向上或向下爬坡后可达到的出力的范围; END
? 附录五：在安全限值内调整的临界值的lingo实现程序
MODEL:
max=a+b+c+d+e+f+g+h; !最大出力;
@abs(0.0826*a+0.0478*b+0.0528*c+0.1199*d-0.0257*e+0.1216*f+0.122*g-0.0015*h+110.477
23
电力市场输电阻塞管理的数学模型
5)<=165;
@abs(-0.0547*a+0.1275*b-0.0001*c+0.0332*d+0.0867*e-0.1127*f-0.0186*g+0.0985*h+131.3521)<=150;
@abs(-108.993-0.0694*a+0.062*b-0.1565*c-0.0099*d+0.1247*e+0.0024*f-0.0028*g-0.2012*h)<=160;
@abs(77.6116-0.0346*a-0.1028*b+0.205*c-0.0209*d-0.012*e+0.0057*f+0.1452*g+0.0763*h)<=155;
@abs(133.1334+0.0003*a+0.2428*b-0.0647*c-0.0412*d-0.0655*e+0.07*f-0.0039*g-0.0092*h)<=132;
@abs(120.8481+0.2376*a-0.0607*b-0.0781*c+0.0929*d+0.0466*e-0.0003*f+0.1664*g+0.0004*h)<=162;
!保证各线路在限值内工作; a<=153; 87<=a; b<=88; 58<=b; 132<=c; c<=228; 60.5<=d; d<=99.5; 98<=e; e<=152; 95<=f; f<=155; 60.1<=g; g<=102.1; 63<=h; h<=117;
!现阶段各机组向上或向下爬坡后可达到的出力的范围; END
? 附录六：在安全裕度内调整的临界值的lingo实现程序
MODEL:
max=a+b+c+d+e+f+g+h; !最大出力;
@abs(0.0826*a+0.0478*b+0.0528*c+0.1199*d-0.0257*e+0.1216*f+0.122*g-0.0015*h+110.4775)<=186.45;
@abs(-0.0547*a+0.1275*b-0.0001*c+0.0332*d+0.0867*e-0.1127*f-0.0186*g+0.0985*h+131.3521)<=177;
@abs(-108.993-0.0694*a+0.062*b-0.1565*c-0.0099*d+0.1247*e+0.0024*f-0.0028*g-0.2012*h)
24