13论文高海浪王芝佳王虎Word格式.docx
- 文档编号:21645245
- 上传时间:2023-01-31
- 格式:DOCX
- 页数:17
- 大小:159.10KB
13论文高海浪王芝佳王虎Word格式.docx
《13论文高海浪王芝佳王虎Word格式.docx》由会员分享,可在线阅读,更多相关《13论文高海浪王芝佳王虎Word格式.docx(17页珍藏版)》请在冰豆网上搜索。
模型1虽然考虑到成本付出问题,但并没有给出成本应该如何确定以及考虑序内容量的分段。
针对这两个问题又提出模型2,在模型2中考虑了不同段价和发电变化量的跨度范围,并通过成本和利润的转化,在一定程度上解决了模型1中的成本问题。
在模型评价中我们还提出了用贡献因子分配费用的办法。
第三问根据各机组的段容量、段价,考虑爬坡速率和预报的负荷需求的约束,设计线性规划数学模型,利用Matlab求得当下一个时段预报的负荷需求是982.4MW时,各机组的出力分配预案如下:
机组1
机组2
机组3
机组4
机组5
机组6
机组7
机组8
150.0000
79.0000
180.0000
99.5000
125.0000
140.0000
95.0000
113.9000
清算价为303元/MWh。
第四问在第三问的基础上考虑安全限值和相对安全裕度,研究如何制订既安全又经济的调度计划,使各个线路潮流值在一定的相对安全裕度内。
首要目标是使阻塞费用最少。
使用线路的安全裕度输电,以避免拉闸限电时,,在制定目标函数时将相对安全裕度利用率(
)考虑在内,即
,在无法调整潮流值使之在安全裕度内时首要目标是使拉闸限电时减少的电量最少以满足用户需求。
当下一个时段预报的负荷需求是982.4MW时,可以通过调整各机组出力分配方案使得输电阻塞消除,调整后的方案为:
150.6896
88.0000
228.0000
80.0276
152.0000
96.6827
70.0000
117.0000
此时的阻塞费用S=11800.14
同样对于第五问,我们得出的分配预案如下:
81.0000
218.2000
135.0000
102.1000
虽然此时不用拉闸,但有线路潮流值超出潮流限值,调整后的方案如下:
153.0000
120.0718
95.2282
此时的阻塞费用S=2013.74。
一、问题重述
电力从生产到使用的四大环节——发电、输电、配电和用电是瞬间完成的。
我国电力市场初期是发电侧电力市场,采取交易与调度一体化的模式。
电网公司在组织交易、调度和配送时,必须遵循电网“安全第一”的原则,同时要制订一个电力市场交易规则,按照购电费用最小的经济目标来运作。
市场交易-调度中心根据负荷预报和交易规则制订满足电网安全运行的调度计划――各发电机组的出力(发电功率)分配方案;
在执行调度计划的过程中,还需实时调度承担AGC(自动发电控制)辅助服务的机组出力,以跟踪电网中实时变化的负荷。
设某电网有若干台发电机组和若干条主要线路,每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力。
电网每条线路上的有功潮流的绝对值有一安全限值,限值还具有一定的相对安全裕度(即在应急情况下潮流绝对值可以超过限值的百分比的上限)。
如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值,称为输电阻塞。
当发生输电阻塞时,需要研究如何制订既安全又经济的调度计划。
你需要做的工作如下:
1.某电网有8台发电机组,6条主要线路,表1和表2中的方案0给出了各机组的当前出力和各线路上对应的有功潮流值,方案1~32给出了围绕方案0的一些实验数据,试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。
2.设计一种简明、合理的阻塞费用计算规则,除考虑上述电力市场规则外,还需注意:
在输电阻塞发生时公平地对待序内容量不能出力的部分和报价高于清算价的序外容量出力的部分。
3.假设下一个时段预报的负荷需求是982.4MW,表3、表4和表5分别给出了各机组的段容量、段价和爬坡速率的数据,试按照电力市场规则给出下一个时段各机组的出力分配预案。
4.按照表6给出的潮流限值,检查得到的出力分配预案是否会引起输电阻塞,并在发生输电阻塞时,根据安全且经济的原则,调整各机组出力分配方案,并给出与该方案相应的阻塞费用。
5.假设下一个时段预报的负荷需求是1052.8MW,重复3~4的工作。
二、模型假设
a)每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻。
b)机组当前出力是对机组在当前时段结束时刻实际出力的预测值。
c)数据基本有效
d)电网公司在组织交易、调度和配送时,必须遵循电网“安全第一”的原则
e)电网每条线路上的有功潮流的绝对值有一安全限值
三、符号说明
第i台机组当前出力值
为第i台机组第j段实际用量
预测第i台机组第j段实际用量
第i台机组第j段的出力系数
第i台机组第j段的段容量
第i台机组第j段报价
清算价
预报的负荷需求
实际交易的负荷
第i台机组爬坡速率
第m条线路的有功潮流的限值
第m条线路相对安全裕度
第m条线路相对安全裕度利用率
第m条线路的有功潮流
阻塞费用
四、问题分析
这道题是在“安全第一”的原则,同时制订一个电力市场交易规则,按照购电费用最小的经济目标来运作,其实质是多目标规划问题。
在分析第一问时,先用散点图描述各线路上有功潮流关于各发电机组出力的关系,从图中可以看出这两者之前具有很强的线性关系,于是我们建立多元线性回归模型,并解出方程,从而能基本反映回归关系。
在分析第二问时,应注意到序内容量不能出力的部分实际上并没有付出发电成本,因此网方不应该赔偿全部费用。
而对于序外容量这里采用的是以差价计算出费用,从而制定了第一个模型。
并且在第一个模型的基础上设计了第二个模型,以解决模型1中关于分段和比例K值的问题。
对于第三问,通过机组爬坡速率和预报的负荷需求的约束,以购电费用最小为主要目标,求得各机组的出力分配预案。
这里采取的是节点电价法[2]。
关于第四问和第五问,由输电阻塞管理原则我们知道要分为三种情况加以讨论建模。
首先是调整各机组出力分配方案使得输电阻塞消除,这对应的是第一个模型。
当无法消除输电阻塞时,还可以使用线路的安全裕度输电,以避免拉闸限电(强制减少负荷需求),目标函数是各条线路对相对安全裕度的利用率中的最大值为最小。
如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度,则必须在用电侧拉闸限电。
这个情况与前两个模型都不一样,主要区别在于首要目标是使强制减少的电量最小。
因为此时要在安全的基础上尽量满足用户需求。
五、模型的建立和求解
(一)问题一:
为了分析Cm与Xi的关系,现利用所给出的表一和表二的数据分别作出Cm对Xi的散点图,这里给出C1对X1与X3的散点图如下:
如图所示我们可以知道Xi与Cm之间有较强的线性关系。
这里考虑常数项问题并建立两个模型:
(1).
(2).
经过编程计算,我们得出
(1)的相关系数小于
(2),说明
(1)的拟合效果不如
(2),由此我们可得出模型1.1如下:
,其中
是常数,
是待确定的系数。
通过Matlab得到有功潮流关于各发电机组出力的近似表达式的系数如下:
110.2965
0.0828
0.0483
0.0530
0.1199
-0.0254
0.1220
0.1216
-0.0012
131.2289
-0.0546
0.1279
-0.0000
0.0333
0.0868
-0.1124
-0.0189
0.0987
-108.8732
-0.0695
0.0616
-0.1566
-0.0099
0.1245
0.0021
-0.0025
-0.2014
77.4817
-0.0345
-0.1024
0.2052
-0.0208
-0.0118
0.0060
0.1449
0.0765
132.9745
0.0005
0.2433
-0.0646
-0.0411
-0.0652
0.0703
-0.0043
-0.0089
120.6633
0.2378
-0.0602
-0.0779
0.0930
0.0469
0.0001
0.1659
0.0007
得到的相关系数R2如下:
1
2
3
4
5
6
0.999509753303878
0.999602436696781
0.999871391216843
0.999887631467390
0.999587794434228
0.999835313607957
相关系数明显大于0.8,说明拟合合程度很高,因此模型2能够比较好地描述有功潮流与各发电机组出力的回归关系。
(二)问题二:
1)模型2.1
当改变根据电力市场交易规则得到的各机组出力分配预案时,一些通过竞价取得发电权的发电容量(称序内容量)不能出力;
而一些在竞价中未取得发电权的发电容量(称序外容量)要在低于对应报价的清算价上出力。
因此,发电商和网方将产生经济利益冲突。
网方应该为因输电阻塞而不能执行初始交易结果付出代价,网方在结算时应该适当地给发电商以经济补偿。
前者我们在计算时只需要序内容量不能出力的部分用清算价计价后的费用计入,因为这部分是发电厂与电网公司之前约定购售的数量,发电商本来应该得到的利润。
但是因为这部分实际上并没有出力,所以实际发电商并没有付出这么多成本,损失只为该费用的一部分,我们用K(K>0且K<1)表示这部分所占比例。
对于序外容量要在低于对应报价的清算价上出力的部分,发电商的损失是由实际报价和清算价之间差价引起的,所以我们以差价计算出费用,具体公式如下:
因此,总的费用为
模型2.1虽然考虑到成本付出问题,当并没有给出K应该如何确定以及具体考虑序内容量不能出力的部分是处在各机组的段容量的哪一段上,此时有两种情况:
a)不能出力的部分位于同一段上
b)不能出力的部分分为两段。
这两种情况显然计价时不能等同的,针对成本和分段问题,我们设计了模型2.2
2)模型2.2:
对于序内容量部分,如果处在两段以上则计费也应分开计费。
我们以第二机组为例子说明。
假设负荷需求是982.4MW时,求得的预案分配中机组2的发电量应该是D1(79),现在由于输电阻塞调整后机组2的实际发电量是D2(65),此时发电的变化量处在第五段和第六段之内,所以应该分为两部分计价。
对于序外容量,现在假设机组2的发电量调整为D3(89),此时增加的发电量分别处在第七段和第十段之内,所以应该分为两部分计价。
所以,这种方法考虑了不同段价和发电变化量的跨度范围,而且在处理序内容量损失费用时,是通过清算价与段价之间差价来衡量成本,这在一定程度上解决了模型2.1中K值的问题。
之所以这样计费是因为把发电商给出的段价作为发电商的成本,而这种算法描述了发电商实际上应得的最少利润,这样发电商能得到最少应得补偿而网方也不会有无谓损失。
这样计算所得的费用更加符合常理,准确度也有一定提高。
设计的模型2.2如下
注:
为
矩阵中的元素,所得到的矩阵元素之和就是要求的阻塞费用。
(三)问题三~五
1.各机组的出力分配预案模型的建立
因为制订的电力市场交易规则,必须按照购电费用最小的经济目标来运作。
这就意味着网方按电力市场付费原则所付给发电商的的购电费用为最小。
此费用为:
由于机组爬坡速率的约束,可能导致选取它的某个段容量的部分。
而为了使得各机组计划出力之和等于预报的负荷需求,清算价对应的段容量可能只选取部分,这些构成了约束条件如下
s.t.:
此时为了更加分析和编程方便,我们采用了
作为第i台机组第j段的出力系数以表示方案中第i台机组第j段用量,即
。
因此,为算出出下一个时段各机组的出力分配预案,我们建立以下模型3.1
2.调整模型的建立
给出处理方案后,因为每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力,电网每条线路上的有功潮流的绝对值有一安全限值,限值还具有一定的相对安全裕度。
因此,可能发生输电阻塞,这就需要研究如何制订既安全又经济的调度计划,并给出阻塞费用。
此时我们选择第二个阻塞费用模型,即
矩阵中的元素
基于以上原则,设计的模型如下:
(1).调整各机组出力分配方案使得输电阻塞消除,模型4.1如下
(2).如果
(1)做不到,可以使用线路的安全裕度输电,以避免拉闸限电(强制减少负荷需求),则在制定目标函数时将相对安全裕度的限制条件考虑在内。
这时有两种算法。
第一种是假设所有线路对相对安全裕度的利用率(即
)表示都会对整个电网的安全产生影响,此时的目标函数应该是
,
若此时我们用木桶定理的思维考虑相对安全裕度的限制,则目标函数应该是各条线路对相对安全裕度的利用率中的最大值为最小,即
.,所以最终模型4.2为
另外,若目标函数是使每条线路上潮流的绝对值超过限值的百分比尽量小,则应为
,这种模型本质上和目标函数是各条线路对相对安全裕度的利用率中的最大值为最小的模型相同,但是后者能够更加明显清楚地说明每条线路上的实际潮流与安全临界值之间的关系,所以我们选择的是模型4.2。
(3).如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度,则必须在用电侧拉闸限电.拉闸限电应该在符合安全的基础上尽量满足用户需求,所以被强制减少的电量应该达到最少的目标。
因此,建立的模型4.3如下:
(四)问题的解决
下一个时段预报的负荷需求是982.4MW时,通过模型3.1编程得到各个机组每个序段的出力系数为元素的矩阵e为:
机组\段
7
8
9
10
1.000
0.000
0.950
0.195
由此得出各机组出力分配方案为:
150.000
79.000
180.000
99.500
125.000
140.000
95.000
113.900
清算价为303元/MWh
代入有功潮流关于各发电机组出力的近似表达式可知各线路的潮流值为(这里给出的潮流值是矢量,有方向):
线路1
线路2
线路3
线路4
线路5
线路6
173.3164
141.0129
150.9312
120.9198
136.8368
168.5309
而由题目得,各线路的潮流限值(单位:
MW)和相对安全裕度
线路
限值
165
150
160
155
132
162
安全裕度
13%
18%
9%
11%
15%
14%
最大值
186.45
177
174.4
172.5
151.8
184.68
比较可知,线路1,5,6超过限值,需要调整各机组出力分配方案。
通过模型调整后各机组出力分配方案为
各线路的潮流值为
165.0000
149.4714
155.0407
126.2649
132.0000
159.7219
此时各线路的潮流值均在限制之内,不用利用相对安全裕度。
通过模型2.2得到此时的阻塞费用S=11800.14
假设下一个时段预报的负荷需求是1052.8MW,各机组出力分配方案为:
81.000
218.200
135.000
102.100
117.000
由有功潮流关于各发电机组出力的近似表达式得到的有功潮流为:
177.2615
141.1833
156.1667
129.7598
134.8507
167.0859
比较可知,线路1,5,6超过限值,但没有超过安全裕度限制范围,此时的安全裕度利用率为
0.5716
0.1440
0.2242
可看出线路1的安全裕度利用率过高,不符合“安全第一”的原则,我们利用模型4.2,只考虑相对安全裕度利用率时,调整后机组出力分配方案为:
各线路的潮流值为:
173.4476
146.8859
155.4084
129.5750
132.7380
166.2697
此时的安全裕度利用率为:
0.3938
0.0373
0.1883
通过模型2.2得到此时的阻塞费用S=2013.74
六、模型评价
在第一问中,我们得到的多元线性回归模型虽然有一定的通用性,但同时又具有很大的随机性,模型不够精确。
在第二问中,虽然第二个模型比较合理地解决了阻塞费用的问题,但是这个模型却是站在发电商的角度来看待的,因为此时计算而得到的费用是通过计算发电商应得而未得利润和不应付而付出损失实现的。
我们在这里重新站在网方角度考虑分配阻塞费用[4]。
具体思路如下:
网方所重视的是序内容量A和序外容量B分别对输电阻塞消除所做的贡献,因此通过模型2得到的总阻塞费用后,还应该利用第一问求出来的有功潮流关于各发电机组出力的近似表达式,通过微分可得各发电机组对调整时有功潮流变化量的影响程度,分别是
进而求得序内容量和序外容量分别对输电阻塞消除所做的贡献因子大小。
分别为:
再通过
和
分配模型2所得的总费用S,这才是站在网方角度的赔偿费用合理归属。
在问题5的解决中,如果在目标函数中将相对安全利用率的取值固定,考虑阻塞费用最小,在计算过程中,我们发现目标函数中相对安全利用率的取值的限制与阻塞费用之间的关系如下
相对安全裕度利用率越低则安全度越高,相对安全裕度利用率越低阻塞费用越高。
综合考虑安全和经济因素,注重安全的网方可以选择相对安全裕度利用率低的,而注重经济的网方可以选择相对安全裕度利用率高的,我们的建议是采用适中的相对安全裕度利用率。
参考文献
[1]姜启源,谢金星,叶俊.数学模型(第三版)[M].北京:
高等教育出版社,2004年.
[2]杨洪明,段献忠,何仰
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 13 论文 海浪 王芝佳王虎