电力市场输电阻塞管理的数学模型

?G?983.4829MW??983.4829MW?G?1094.600MW?G?1094.600MW?在安全限值内调整消除输电阻塞在安全裕度内进行调整电侧拉闸限电

五． 模型的评价

5.1 模型的优缺点

优点：

? 思路清晰，方法简明易懂。本模型运用合理的假设，将相当复杂的电

力市场输电管理、出力分配，运用线性规划的方法，在有限的规划式子里面得到了较为全面的解答。

? 在设计阻塞费用计算方法时，对序内容量不能出力与序外容量出力部

分的阻塞费用以不同方案计算，并考虑其单位发电成本的变动，计算方法较简明、公平。

? 当阻塞费用计算方法改变时，模型仍然适用。

? 本模型求解有功潮流与机组出力的相关函数时，比较了各种方案，同

时了解其物理意义，得出的函数回归效果显著。 缺点：

? 本模型未给出将有功潮流与各机组出力视为线形关系的适用范围。 ? 计算阻塞费用时未考虑机组出力调整幅度。

5.2 模型的改进

此模型可以在阻塞费用的计算方法上进行改进，即使各机组计算阻塞费用时补偿系数随出力调整幅度增大而增大，更加公平、合理。

参考文献

[1] 王文波，数学建模及其基础知识详解，武汉大学出版社，2006年5月 [2] 西北工业大学数学建模指导委员会，数学建模简明教程，高等教育出版社，2008年9月

[3] 耿建 王锡凡 陈皓勇 丁晓莺，分段竞价电力市场的运营方式、规则和模型，西安交通大学学报第37卷第6期，2003年6月。

[4] 肖鸣 李林川 孙伟 茅波，分段电力市场的出清算法实现，电力系统自动化第27卷第23期，2003年12月10日。

[5] 黄光东 李平 史川北 霍振中 詹浩凯，输电线路阻塞费用的管理，数学的实践与认识第35卷第6期，2005年6月。

17

电力市场输电阻塞管理的数学模型

附录

? 附录一：实际收入Hj的8个分段函数

???505X10?X1?70H?124X170

电力市场输电阻塞管理的数学模型

??500X7?120X7??180X7??251X7?H7??260X7?306X7??315X7??335X7??348X70?X7?5050?X7?6565?X7?7570?X7?8585?X7?9595?X7?105105?X7?110110?X7?120120?X7?123??800X8?183X8??253X8? H8??303X8?318X8??400X8??800X80?X8?7070?X8?9090?X8?110110?X8?130 130?X8?140140?X8?155155?X8?160? 附录二：各线路上有功潮流关于各发电机组出力表达式的matlab实现程序

x=[1 120 73 180 80 125 125 81.1 90 1 133.02 73 180 80 125 125 81.1 90 1 129.63 73 180 80 125 125 81.1 90 1 158.77 73 180 80 125 125 81.1 90 1 145.32 73 180 80 125 125 81.1 90 1 120 78.596 180 80 125 125 81.1 90 1 120 75.45 180 80 125 125 81.1 90 1 120 90.487 180 80 125 125 81.1 90 1 120 83.848 180 80 125 125 81.1 90 1 120 73 231.39 80 125 125 81.1 90 1 120 73 198.48 80 125 125 81.1 90 1 120 73 212.64 80 125 125 81.1 90 1 120 73 190.55 80 125 125 81.1 90 1 120 73 180 75.857 125 125 81.1 90 1 120 73 180 65.958 125 125 81.1 90 1 120 73 180 87.258 125 125 81.1 90 1 120 73 180 97.824 125 125 81.1 90 1 120 73 180 80 150.71 125 81.1 90 1 120 73 180 80 141.58 125 81.1 90 1 120 73 180 80 132.37 125 81.1 90 1 120 73 180 80 156.93 125 81.1 90 1 120 73 180 80 125 138.88 81.1 90 1 120 73 180 80 125 131.21 81.1 90 1 120 73 180 80 125 141.71 81.1 90 1 120 73 180 80 125 149.29 81.1 90 1 120 73 180 80 125 125 60.582 90 1 120 73 180 80 125 125 70.962 90 1 120 73 180 80 125 125 64.854 90 1 120 73 180 80 125 125 75.529 90

19

电力市场输电阻塞管理的数学模型

1 120 73 180 80 125 125 81.1 104.84 1 120 73 180 80 125 125 81.1 111.22 1 120 73 180 80 125 125 81.1 98.092 1 120 73 180 80 125 125 81.1 120.44 ] ;

y1=[164.78 165.81 165.51 167.93 166.79 164.94 164.8 165.59 165.21 167.43 165.71 166.45 165.23 164.23 163.04 165.54166.88 164.07 164.27 164.57 163.89 166.35 165.54 166.75 167.69 162.21 163.54 162.7 164.06 164.66 164.7 164.67 164.69 ];

y2=[140.87 140.13 140.25 138.71 139.45 141.5 141.13 143.03 142.28 140.82 140.82 140.82 140.85 140.73 140.34 141.1 141.4 143.03 142.29 141.44 143.61 139.29 140.14 138.95 138.07 141.21141 141.14 140.94 142.27 142.94141.56 143.84 ];

y3=[-144.25 -145.14 -144.92 -146.91 -145.92 -143.84 -144.07 -143.16 -143.49 -152.26 -147.08 -149.33 145.82 -144.18 -144.03 -144.32 -144.34 -140.97 -142.15 -143.3 -140.25 -144.2 -144.19 -144.17 -144.14 -144.13 -144.16 -144.21 -144.18 -147.2 -148.45 -145.88 -150.34 ]; y4=[119.09 118.63 118.7 117.72 118.13 118.43 118.82 117.24 117.96 129.58 122.85 125.75 121.16 119.12 119.31 118.84 118.67 118.75 118.85 119 118.64 119.1 119.09 119.15 119.19 116.03 117.56 116.74 118.24 120.21 120.68 119.68 121.34 ];

y5=[135.44 135.37 135.33 135.41 135.41 136.72 136.02 139.66 137.98 132.04 134.21 133.28 134.75 135.57 135.97 135.06 134.67 133.75 134.27 134.88 133.28 136.33 135.81 136.55 137.11 135.5 135.44 135.4 135.4 135.28 135.16 135.29 135.12 ];

y6=[157.69 160.76 159.98 166.81 163.64 157.22 157.5 156.59 156.96 153.6 156.23 155.09 156.77 157.2 156.31 158.26 159.28 158.83 158.37 158.01 159.12 157.59 157.67 157.59 157.65 154.26 155.93 154.88 156.68 157.65 157.63 157.61 157.64 ]; [b bint r rint stats]=regress(y1’,x); [b bint r rint stats]=regress(y2’,x); [b bint r rint stats]=regress(y3’,x); [b bint r rint stats]=regress(y4’,x); [b bint r rint stats]=regress(y5’,x); *b bint r rint stats+=regress(y6’,x);

? 附录三：在安全限值内调整的lingo实现程序

MODEL:

min=(b1+b2+b3+b4+b5+b6+b7+b8)*0.25;!使阻塞费用最小;

p1=a*303; p2=b*303; p3=c*303; p4=d*303; p5=e*303; p6=f*303; p7=g*303;

p8=h*303;!各机组每小时实际收入;

20