电力市场分配问题

基于多元线性回归模型解决电力分配问题

【摘要】

对问题一，在题目给出的表1和表2中我们很容易看出对于机组的出力每一确定值，线路的有功潮流是一组随机变量，而可能遵循一定的分布而取值。也就是说，对于机组出力的值确定后，而线路的有功潮流不能随着确定，对于此题使用多元线性回归的知识来解决，多元线性回归的模型是,所以可以将本题所给出表1和表2中各发电机组的出力数据与各条线路的有功潮流导入进表格中，使用表格中的数据分析工具，求出每条线路多元线性回归方程的系数与常数，列出测试值与比较值使用残差分析比较，观察多元线性回归表达式准确度，准确度高之后可以列出各线路上有功潮流关于各发电机组出力的近似表达式。

针对问题二，当执行各机组的分配预案时，如果主要线路上的有功潮流超过了限值，网方会重新对分配预案进行调整，改变初始交易规则，产生阻塞费用，考虑的公平性，设计了序内容量的阻塞费用计算规则为：，序内容量的阻塞费用计算规则为：，详细的设计思路见正文。

针对问题三，在给出下一时段预报的负荷需求982.4，根据各机组的段容量、段价和爬坡速率的数据，按照电力市场交易规则，给出下一时段各机组的出力分配预案。首先按照市场交易规则，按段价从低到高选取各个机组的段容量，直到选取的段容量之和等于负荷需求982.4，通过方案0和各个机组的爬坡速率算出下一时刻出力的最大值和最小值，让选取每个机组的段容量之和小于或等于下一时刻的出力的最大值，最终得出当负荷需求为982.4时，各机组的出力分配预案为：机组一是150、机组二是79、机组三是180、机组四是99.5、机组五是125、机组六是140、机组七是95、机组八是113.9，清算价是303元，购电费用为74416.8元，详细的计算方法见正文。

针对问题五，当下一时段的预报负荷需求为1052.8，与问题三计算各个机组的分配预案的方法类似，得出清算价是356元，购电费用为93699.2元，各机组的分配预案为：机组一是150、机组二是81、机组三是218.2、机组四是99.5、机组五是135、机组六是150、机组七是102.1、机组八是117。

关键词：软件 多元线性回归 lingo程序 excel表格

 

 

一、问题重述

1引言

我国电力系统的市场化改革正在积极、稳步地进行。2003年3月国家电力监管委员会成立，2003年6月该委员会发文列出了组建东北区域电力市场和进行华东区域电力市场试点的时间表，标志着电力市场化改革已经进入实质性阶段。可以了解，随着我国用电紧张的缓解，电力市场化将进入新一轮的发展，这给有关产业和研究部门带来了可预期的机遇和挑战。

市场交易-调度中心根据负荷预报和交易规则制订满足电网安全运行的调度计划――各发电机组的出力（发电功率）分配方案；在执行调度计划的过程中，还需实时调度承担AGC（自动发电控制）辅助服务的机组出力，以跟踪电网中实时变化的负荷。设某电网有若干台发电机组和若干条主要线路，每条线路上的有功潮流（输电功率和方向）取决于电网结构和各发电机组的出力。

电网每条线路上的有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度（即在应急情况下潮流绝对值可以超过限值的百分比的上限）。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值，称为输电阻塞。当发生输电阻塞时，需要研究如何制订既安全又经济的调度计划。了解给出的电力市场交易规则、市场交易-调度中心在当前时段内要完成的具体操作过程、输电阻塞管理原则完成以下题目。

2问题提出

（1）某电网有8台发电机组，6条主要线路，表1和表2中的方案0给出了各机组的当前出力和各线路上对应的有功潮流值，方案给出了围绕方案0的一些实验数据，试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。

（2）设计一种简明、合理的阻塞费用计算规则，除考虑上述电力市场规则外，还需注意：在输电阻塞发生时公平地对待序内容量不能出力的部分和报价高于清算价的序外容量出力的部分。

   （3）假设下一个时段预报的负荷需求是982.4，表3、表4和表5分别给出了各机组的段容量、段价和爬坡速率的数据，试按照电力市场规则给出下一个时段各机组的出力分配预案。

（4）按照表6给出的潮流限值，检查得到的出力分配预案是否会引起输电阻塞，并在发生输电阻塞时，根据安全且经济的原则，调整各机组出力分配方案，并给出与该方案相应的阻塞费用。

（5）假设下一个时段预报的负荷需求是1052.8，重复的工作。

二、问题分析

1针对问题一进行分析

由上表我们很容易看出对于机组的出力每一确定值，线路的有功潮流是一组随机变量，而可能遵循一定的分布而取值。也就是说，对于机组出力的值确定后，而线路的有功潮流不能随着确定，对于本题就可以使用多元线性回归的知识来解决。

所以可以将本题所给出表1和表2中各发电机组的出力数据与各条线路的有功潮流导入进表格中，使用表格中的数据分析工具，建立多元线性回归表达式并运用工具中的回归分析求出相关信息，能够列出各线路上有功潮流关于各发电机组出力的近似表达式。

2针对问题二进行分析

问题2是要在公平对待序内容量和序外容量的基础上，设计简明、合理的阻塞费用计算规则。

电网每条线路上的有功潮流的绝对值具有一安全限值，限值还具有一定的相对安全裕度（即在应急情况下潮流绝对值可以超过限值的百分比的上限）。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出了限值，就称之为输电阻塞。

对不能出力的序内容量，本应得到的利润无法得到，所以损失费用就为清算价与减少出力值的乘积，总损失就是所有不能出力的段容量的损失之和。  

对序外容量，其损失费用就是报价与清算价的差额与容量的乘积，总损失就是所有低于清算价出力的段容量的损失之和。

对于本题就需要对序内容量和序外容量，分别赔偿其所有损失。

3针对问题三进行分析

假设下一个时段预报的负荷需求是982.4，需要求出按照电力市场规则给出下一个时段各机组的出力分配预案，并使购电花费的价钱最少。

根据电力市场规则可以知道下个时段预报的负荷等于每台机组的报价加上当前机组的出力再加上当前机组出力的改变速率。

不需要考虑到线路输电安全的问题。由当前各个机组出力值和各机组的爬坡速度，可以得到各机组下一时段出力值的范围。在不超出范围的前提下，按段价从低到高选取各级组的段容量，直到机组各段容量或其部分之和等于预报负荷需求，最后一个被选入的段价，即为清算价。由此可以知道能够建立优化模型，使用软件求解，得到分配预案。

4针对问题四进行分析

此问题建立在问题1，2，3的基础上，是根据给定的潮流值，判断问题3得到的分配预案是否会引起输电阻塞。将问题3得到的预案中的数据代入由问题1得出的有功潮流的计算表达式，将结果与题中所给表6给定的潮流限值比较，就可以判断此时是否发生输电阻塞。如果有阻塞，就给出调整后的方案，并计算阻塞费用，可能需要对不同的情况建立不同的线性优化模型，列出会发生的所以情况：

（1）当输电阻塞能够消除时，优化的目标是阻塞费用最小，对模型使用软件求解。

（2）当输电阻塞不能够消除时，但是可以使每条线路的潮流绝对值超过限值的百分比小于相对安全裕度时，以每条线路上潮流的绝对值超过限值的百分比最小作为优化目标，以确保各个线路相对安全，避免拉闸限电，对模型使用软件求解。

（3）当无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度，必须在用电侧拉闸限电时，以限制的电量最小为优化目标。对模型使用软件求解。

5针对问题五进行分析

此题相当于将预报负荷需求量增大到1052.8重复3，4问题流程。

三、符号说明

符号	代表含义
	第条线路的有潮功流
	第台发电机组的出力
	系数
	发电机组
	线路
	清算价
	第个机组第个时段的报价
	第个机组第个时段的段容量
	第个机组的爬坡速率
	方案0第个机组的出力
	第个机组第个时段的实际出力
	双下标”0-1“变量

四、模型假设

1假设本题给出的数据准确无误；

2假设发电机组出力不发生损耗；

3假设不考虑白天与晚上消耗电量的差异；

4假设发电机组持续发电不断电；

5假设阻塞时阻塞费用一定存在；

6假设每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻；

7假设每台机组单位时间内能增加或减少的出力相同；

8假设线路输电时绝对安全；

9各机组爬坡时速率稳定不变。

五、模型的建立与求解

1求出各线路的有功潮流与各发电机组出力的近似表达式

1.1多元线性回归模型

表1与表2中的每一条线路与各发电机组存在联系， 电网有8台发电机组，6条主要线路，表1与表2的方案0给出了当前的机组的出力与有功潮流的关系，在现代的电力运输网络里是看不到机组的位置的电网公司负责电力的运输调度和配送，而机组发出的电力不再向特定的某条线路输送，而是由电网运营商制定的方案向整个电网传送，一台机组的出力发生变化就会影响到全部电网的有功潮流，表1与表2其余32种方案，发生的微调就引起了6 条电网的变化。

由此可以知道，对于此题使用多元线性回归来列出每条线路与各发电机组出力的方程，了解到多元线性回归的模型为：

1.2 表格求系数

根据题目中所给出表一表二的数据，可以将各机组出力与各线路的潮流值所有数据导入进表格中，使用表格中的数据分析求得。

第一条线路的到的所有数据：

=110.4775，=0.0826，=0.0478，=0.0528，=0.1199，=-0.0257，=0.1216，=0.1220，=-0.0015。

第二条线路的到的所有数据：

=131.3521，=-0.0547，=0.1275，=-0.0001，=0.0332，=0.0867，=-0.1127，=-0.0186，=0.0985。

第三条线路的到的所有数据：

=-108.993，=0.0694，=-0.062，=-0.1565，=-0.0099，=0.1247，=0.0024，=-0.0028，=-0.2012。

第四条线路的到的所有数据：

=77.6116，=-0.0346，=-0.1028，=0.205，=-0.0209，=-0.012，=0.0027，=0.1452，=0.0763。

第五条线路的到的所有数据：

=133.1334，=0.0003，=0.2428，=-0.0647，=-0.0412，=-0.0655，=0.07，=-0.0039，=-0.0092。

第六条线路的到的所有数据：

=120.8481，=0.2376，=-0.0607，=-0.0781，=0.0929，=-0.0466，=-0.0002，=0.1664，=-0.0004。

1.3系数矩阵与常数矩阵

根据以上使用表格求出的系数数据可以列出系数矩阵，与常数矩阵。

系数矩阵：

常数矩阵:

1.4对残差图进行分析

对于多元线性回归方程是否准确可以通过分析此时段预测值与实际值的残差来判断。

图1 线路6残差图

由上残差图与附件中其余时段的残差图可以知道，预测值与实际值相差不大，没有极度偏移的残差点，说明对于本题建立的每条线路与各发电机组出力的表达式准确度高，误差值小。

1.5列出近似表达式：

2设计阻塞费用计算规则

2.1序内容量不能出力的部分

对不能出力的序内容量，本应得到的利润无法得到，所以损失费用就为清算价与报价的差额与容量的乘积，总损失就是所有不能出力的段容量的损失之和，设需要减少的出力为，为清算价，即对序内容量需要支付的阻塞费用为：

  

2.2序外容量出力的部分

对序外容量，其损失费用就是报价与清算价的差额与容量的乘积，总损失就是所有低于清算价出力的段容量的损失之和，设需要增加的出力为，对应的报价为，即对序外容量需要支付的阻塞费用为：

3建立购电费用最小优化模型求出各机组的出力分配方案

3.1目标函数的确定

题目给出了下一个时段预报的负荷需求值，负荷需求值为，根据表3、表4和表5分别给出的各机组的段容量、段价和爬坡速率的数据，按照电力市场规则给出下一个时段各机组的出力分配预案。

了解到下个时段预报的负荷等于每台机组的报价加上当前机组的出力再加上当前机组出力的改变速率，因为要使各机组分配最优化而且花费的价钱也要最低，花费的价钱等于负荷需求值乘以清算价，所以可以使花费费用为目标，列出以下目标函数：

3.2条件的约束

（1）实际出力值与爬坡速率之间的约束

方案0为当前时刻各机组的出力方案，根据方案0和各机组的爬坡速率约束，算出下一时刻各机组的最大出力值和最小出力值，即下一时刻所有机组的实际出力值不能够大于方案0的各机组的出力值加上爬坡正值，小于方案0的各机组的出力值加上爬坡负值，即具体约束如下：

（2）段容量的约束

根据电力市场交易规则知道，每台机组在当前时段开始时刻前给出下一个时段的报价与段容量，所以在不同时段的段容量是根据上一时段的段容量所决定的为了使得各机组计划出力之和等于预报的负荷需求，清算价对应的段容量可能只选取部分，当时，不选择第组第段的段容量，当时，各机组的实际出力小于或等于下一时段的段容量，即具体约束如下：

（3）清算价的约束

在电力市场交易规则可以知道最后一个被选入的段价（最高段价）称为该时段的清算价，该时段全部机组的所有出力均按清算价结算，因为各机组将可用出力由低到高分成至多10段报价，每个段的长度称为段容量，每个段容量报一个价（称为段价），段价按段序数单调不减，所以清算价相当于机组报价的最大值，即具体约束如下：

（4）预报的负荷需求值的约束

由题目可以知道给出下一时刻预报的负荷需求值等于982.4，并且负荷需求值等于所有机组实际出力值的总和，即具体约束如下：

（5）“0-1”变量的约束

根据电力市场交易规则知道，每台机组在当前时段开始时刻前给出下一个时段的报价，若是在该时段没有选择的机组和选择机组的段容量为0的机组，就在下一时段的报价为0，在这一情况下设置一个“0-1”变量来决定出下一时段的报价是否为0，即具体约束如下：

3.3模型的建立

目标：

 

3.4各机组出力分配预案

由当前各个机组出力值和各机组的爬坡速度，可以得到各机组下一时段出力值的范围，从题目中给出的注释内容可以知道各机组计划出力之和等于预报的负荷需求。

了解到电力市场交易规则以15分钟为一个时段组织交易，每台机组在当前时段开始时刻前给出会下一个时段的报价。各机组将可用出力由低到高分成至多10段报价，每个段的长度称为段容量，每个段容量报一个价（称为段价），段价按段序数单调不减，因为在时段一时刻报价数最低，预估时段一中所有发电机组都会被选出，根据花费价钱最低还要使预报负荷达标，所以我们根据题目中的表4选择在时段内报价最少的发电机组，根据各机组的段容量与各机组的爬坡速率可以求出各机组在不同时段的出力值，并且使各个发电机组出力值的总和等于预报出的负荷需求。

列出以下表格：

表1 分配完成的机组下一时段段容量(单位：兆瓦）

机组\段	1	2	3	4	5	6	7
1	70	0	50	0	0	30	0
2	30	0	20	8	15	6	0
3	110	0	40	0	30	0	0
4	55	5	10	10	10	9.5	0
5	75	5	15	0	15	15	0
6	95	0	10	20	0	15	0
7	50	15	5	15	10	0	0
8	70	0	20	0	20	0	3.9

各机组的出力分配预案为：机组一是150、机组二是79、机组三是180、机组四是99.5、机组五是125、机组六是140、机组七是95、机组八是113.9，上述表格可以知道在不同时段分配时分配的机组是不同的，在时段一分配到的机组为1至8所有机组；在时段二分配到的机组为机组4、机组5与机组7；在时段三分配到的机组为1至8所有机组；在时段四分配到的机组为机组2、机组4、机组6与机组7；在时段五分配到的机组为机组2、机组3、机组4、机组5、机组7与机组8；在时段六分配到的机组为机组1、机组2、机组4、机组5与机组6；在时段7分配到的机组为机组8；最后三个时段不分配机组。并且以上总值等于预报的电荷需求值符合题意。

3.5下一时刻分配预案算法的具体描述

Step1：首先按照各机组的段价从低到高选取各机组的段容量或其部分，直到它们之和等于预报的负荷。

Step2：跟据当前时刻各机组的出力值和各个机组的爬坡速率计算出下一时刻各机组出力的最大值和最小值。

Step3：在已入选的段容量中对每一个机组的段容量求和，观察段容量的和是否超出下一时刻的最大出力值，若超出下一时刻的最大出力值，则只取该机组已入选的段容量部分让段容量等于该机组在下一时刻的最大出力值。

Step4：不足下一时刻负荷需求的部分再继续从各机组段价中从低到高选择，已达到下一时刻的最大出力值的机组，不在选取该机组的段容量，直到满足下一时刻的负荷需求。

3.6具体结果

通过以上这种算法，得出清算价为303元，购电费用为74416.8元，和下一时刻的各机组的出力分配预案为：机组一是150、机组二是79、机组三是180、机组四是99.5、机组五是125、机组六是140、机组七是95、机组八是113.9。

4调整各机组出力分配方案，给出与该方案相应的阻塞费用

4.1判断是否阻塞

将问题三给出的各机组分配预案代进问题算出的有功潮流与各机组的近似表达式算出有功潮流，再根据题目中的表6给出的各线路的潮流限值进行比较，判断是否超过限值。

5预报的负荷需求改变求出对应的各机组的出力分配方案

六、模型的优点与改进

1 模型的优点

（1）能够有效的运用题目中的数据；

（2）对于列出的各线路上有功潮流关于各发电机组出力的近似表达式误差值小；

（3）模型造成的空间复杂度与时间复杂度低；

（4）有理有据，令人服气

附件

1 软件代码

问题三最小购电费用模型的Lingo代码

sets:

aa/1..8/:v,w,c;

bb/1..10/:;

cc(aa,bb):p,b,x,f;

endsets

data:

v=2.2,1,3.2,1.3,1.8,2,1.4,1.8;

w=120,73,180,80,125,125,81.1,90;

p=@ole('D:\math\·¢µç.xlsx','pp');

b=@ole('D:\math\·¢µç.xlsx','bb');

enddata

min=@max(cc(i,j):p(i,j)*f(i,j));

@for(cc(i,j):x(i,j)<=b(i,j)*f(i,j));

@for(aa(i):@sum(bb(j):x(i,j))<=w(i)+v(i)*15);

@for(aa(i):@sum(bb(j):x(i,j))>=w(i)-v(i)*15);

@sum(cc(i,j):x(i,j))=982.4;

@for(cc(i,j):@bin(f(i,j)));

@for(aa(i):c(i)=@sum(bb(j):x(i,j)));

2 各时段有潮功流与各机组出力值近似表达式残差值

图2 线路1残差图

图3 线路2残差图

图4 线路3残差图

图5 线路4残差图

图6 线路5残差图