电力市场的输电阻塞管理

符号说明及名词解释
机组出力： 发电功率；
有功潮流： 输电功率和方向；
DC模型：  直流最优潮流（DC-OPF）；
阻塞成本（congestion cost）：由于线路容量限制，使得原先报价低成本的机组，不得不减少发电量，转由其它报价高的机组多发，由此导致的发电成本的增加；
相对安全裕度：即在应急情况下潮流绝对值可以超过限值的百分比的上限；

问题重述
    世界范围内的电力工业改革, 建立了各种形式的电力市场, 其核心目的是: 引入竞争, 提高效率, 降低电价, 改善服务, 持续发展。电力市场的各参与者以获得最大利润为直接目的, 而降低电价、改善服务又将使全社会受益。
随着我国电力市场改革的深入“出现了许多新的问题。电力市场中网络潮流的分布主要取决于电力交易的分布”而电力交易以利润最大化为导向"于是不可避免地出现系统中某些输电价格较低的线路或某些电能价格较低的发电机周围的线路往往承载着较重的负荷"增加了发生电阻塞的机会。所谓电阻塞就是假设某电网有若干台发电机组和若干条主要线路，每条线路上的有功潮流（输电功率和方向）取决于电网结构和各发电机组的出力。电网每条线路上的有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度（即在应急情况下潮流绝对值可以超过限值的百分比的上限）。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值，称为输电阻塞。同时电力已经作为商品市场化——“电力市场中系统各成员之间的关系是纯粹的经济利益关系”， 传统的以行政命令方式通过发电再调度来消除输电阻塞的方法在电力市场中已不再奏效，必须利用电价这一杠杆来对市场各成员进行协调。所以问题的实质就是建立依靠合理的交易与调度一体化的模式。电网公司在组织交易、调度和配送时，必须遵循电网“安全第一”的原则，同时要制订一个电力市场交易规则，按照购电费用最小的经济目标来运作。电网公司根据负荷预报和交易规则制订满足电网安全运行的调度计划——各发电机组的出力（发电功率）分配方案。
本文将利用问题中给出的数据和规则建立一套合理的机制希望能对我国电力市场的改革和发展添上自己的一块砖。

模型假设
基本假设：
1.遵循安全用电的原则，在超出 安全限值*（1+安全裕度）,调配无效时，要拉闸；
2.不考虑用户用电的优先权；
3.无白天和晚上用电价格的差异；
4.不考虑电力在传输过程中的损失；
5.无需考虑用电时段；
6.假设6条线路是相互连通的，但无需考虑电网的拓扑结构；
7.假设所有的电力调度是在瞬时完成的。

模型的建立与求解
问题一：
    电网有8台发电机组，6条主要线路，表1与表2的方案0给出了当前的机组的出力与有功潮流的关系，在现代的电力运输网络里是看不到机组的位置的电网公司负责电力的运输调度和配送，而机组发出的电力不再向特定的某条线路输送，而是由电网运营商制定的方案向整个电网传送，这样就很容易理解一台机组的出力发生变化就会影响到全部电网的有功潮流，表1与表2其余32种方案，发生的微调就引起了6 条电网的变化。
如图所示：
表1  各机组出力方案 （单位：兆瓦，记作MW）
方案/机组 1 2 3  4  5  6 7 8
0 120 73 180 80 125 125 81.1 90
1 133.02 73 180 80 125 125 81.1 90
2 129.63 73 180 80 125 125 81.1 90
3 158.77 73 180 80 125 125 81.1 90
4 145.32 73 180 80 125 125 81.1 90
5 120 78.596 180 80 125 125 81.1 90
6 120 75.45 180 80 125 125 81.1 90
7 120 90.487 180 80 125 125 81.1 90
8 120 83.848 180 80 125 125 81.1 90
9 120 73 231.39 80 125 125 81.1 90
10 120 73 198.48 80 125 125 81.1 90
11 120 73 212.64 80 125 125 81.1 90
12 120 73 190.55 80 125 125 81.1 90
13 120 73 180 75.857 125 125 81.1 90
14 120 73 180 65.958 125 125 81.1 90
15 120 73 180 87.258 125 125 81.1 90
16 120 73 180 97.824 125 125 81.1 90
17 120 73 180 80 150.71 125 81.1 90
18 120 73 180 80 141.58 125 81.1 90
19 120 73 180 80 132.37 125 81.1 90
20 120 73 180 80 156.93 125 81.1 90
21 120 73 180 80 125 138.88 81.1 90
22 120 73 180 80 125 131.21 81.1 90
23 120 73 180 80 125 141.71 81.1 90
24 120 73 180 80 125 149.29 81.1 90
25 120 73 180 80 125 125 60.582 90
26 120 73 180 80 125 125 70.962 90
27 120 73 180 80 125 125 64.854 90
28 120 73 180 80 125 125 75.529 90
29 120 73 180 80 125 125 81.1 104.84
30 120 73 180 80 125 125 81.1 111.22
31 120 73 180 80 125 125 81.1 98.092
32 120 73 180 80 125 125 81.1 120.44

表2  各线路的潮流值（各方案与表1相对应，单位：MW）
方案/线路 1 2 3 4 5 6
0 164.78 140.87 -144.25 119.09 135.44 157.69
1 165.81 140.13 -145.14 118.63 135.37 160.76
2 165.51 140.25 -144.92 118.7 135.33 159.98
3 167.93 138.71 -146.91 117.72 135.41 166.81
4 166.79 139.45 -145.92 118.13 135.41 163.64
5 164.94 141.5 -143.84 118.43 136.72 157.22
6 164.8 141.13 -144.07 118.82 136.02 157.5
7 165.59 143.03 -143.16 117.24 139.66 156.59
8 165.21 142.28 -143.49 117.96 137.98 156.96
9 167.43 140.82 -152.26 129.58 132.04 153.6
10 165.71 140.82 -147.08 122.85 134.21 156.23
11 166.45 140.82 -149.33 125.75 133.28 155.09
12 165.23 140.85 -145.82 121.16 134.75 156.77
13 164.23 140.73 -144.18 119.12 135.57 157.2
14 163.04 140.34 -144.03 119.31 135.97 156.31
15 165.54 141.1 -144.32 118.84 135.06 158.26
16 166.88 141.4 -144.34 118.67 134.67 159.28
17 164.07 143.03 -140.97 118.75 133.75 158.83
18 164.27 142.29 -142.15 118.85 134.27 158.37
19 164.57 141.44 -143.3 119 134.88 158.01
20 163.89 143.61 -140.25 118.64 133.28 159.12
21 166.35 139.29 -144.2 119.1 136.33 157.59
22 165.54 140.14 -144.19 119.09 135.81 157.67
23 166.75 138.95 -144.17 119.15 136.55 157.59
24 167.69 138.07 -144.14 119.19 137.11 157.65
25 162.21 141.21 -144.13 116.03 135.5 154.26
26 163.54 141 -144.16 117.56 135.44 155.93
27 162.7 141.14 -144.21 116.74 135.4 154.88
28 164.06 140.94 -144.18 118.24 135.4 156.68
29 164.66 142.27 -147.2 120.21 135.28 157.65
30 164.7 142.94 -148.45 120.68 135.16 157.63
31 164.67 141.56 -145.88 119.68 135.29 157.61
32 164.69 143.84 -150.34 121.34 135.12 157.64

由上图我们很容易看出对于机组的出力每一确定值，线路的有功潮流是一组随机变量，而可能遵循一定的分布而取值。也就是说，对于机组出力的值确定后，而线路的有功潮流不能随着确定[1]。
因此，我们可以选用回归分析的方法来确定有功潮流关于各发电机组的关系。我们利用Matlab的统计工具箱中提供的regress()函数来求解回归系数的最小二乘估计并建立回归的线性模型，误差估计及其置信区间，给出系数（复相关系数）统计量值[2]。
首先，建立m文件命名为New111（New111是线路1的有潮功流和各机组的出力建立的m文件，New112等类同）：
 我们给矩阵X的第一列加一列1，这样就在回归方程中出现了一个常数项，而我们从残差分析的结果中可以看出这种做法使得结果更接近于真实的情况。虽然New113和New114得出的残差图中出现了坏点，但从eint所得的值中很容易发现坏点的残差值New113的eint=-0.0850   -0.0056，非常的接近于0，这样坏点对整个问题的影响就可以忽略了。
这里只给出一个m文件的内容，余下的见附录。
New111：
y=[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.54 166.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 ];
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
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];
[b,bint,e,eint,stats]=regress(y,X,0.05)
rcoplot(e,eint)
回归残差图：
 

New111
y1=110.4775+0.0826x1+0.0478x2+0.0528x3+0.1199x4+-0.0257x5+0.1216x6+0.1220 x7-0.0015x8
b =110.4775 0.0826 0.0478   0.0528  0.1199  -0.0257  0.1216  0.1220  -0.0015

New112
y2=131.3521-0.0547x1+0.1275x2-0.0001x3+0.0332x4+0.0867x5-0.1127x6-0.0186x7+0.0985x8
b =131.3521  -0.0547  0.1275  -0.0001    0.0332  0.0867  -0.1127  -0.0186   0.0985

New113
y3=108.9928+0.0694x1-0.0620x2+0.1565x3+0.0099x4-0.1247x5-0.0024x6+0.0028x7+0.2012x8
b = 108.9928  0.0694  -0.0620  0.1565  0.0099  -0.1247 -0.0024  0.0028   0.2012

New114
y4=77.6116-0.0346x1-0.1028x2+0.2050x3-0.0209x4-0.0120x5+0.0057x6+0.1452x6+0.0763x8
b =77.6116  -0.0346  -0.1028   0.2050   -0.0209  -0.0120   0.0057    0.1452   0.0763

New115
y5=133.1334+0.0003x1+0.2428x2-0.0647x3-0.0412x4+-0.0655x5+0.0700x6-0.0039x7-0.0092x8
b =133.1334   0.0003    0.2428  -0.0647  -0.0412  -0.0655   0.0700  -0.0039  -0.0092

New116
y6=120.8481+0.2376x1-0.0607 x2 -0.0781x3+0.0929x4+0.0466x5-0.0003 x6+ 0.1664x7+0.0004x8
b =120.8481   0.237