区域电力市场双边交易阻塞管理实用计算方法_张永平

合集下载

电力市场的输电阻塞管理答案审批稿

电力市场的输电阻塞管理答案审批稿

电力市场的输电阻塞管理答案YKK standardization office【 YKK5AB- YKK08- YKK2C- YKK18】电力市场的输电阻塞管理摘要本文讨论电力市场中输电阻塞管理问题。

通过合理的假设简化,根据交易规则和阻塞管理原则建立数学模型。

对问题一,本文利用多元线性回归得到各线路有功潮流关于各发电机组出力的函数,利用matlab 软件对回归方程进行了残差分析,得到各路线有功潮流关于各发电机组出力的近似表达式。

对问题二,按照购电费用小的经济原则,本着公平的原则建立了两个阻塞费用计算规则对两种不同的情形进行分析,引进正、负偏差变量法,根据双方满意原则,得到合理的阻塞费用计算规划模型。

对问题三,依据电网公司遵循的最小化原则,以购电费用最小为目标,以满足负荷需求为约束条件,建立最优化模型,通过MATLAB 程序编程,得到各组分配最优预案,得到各组分配预案分别为:150、79、18、、125、14、95、(单位:兆瓦),这时清算价格为MWh /303元。

对问题四:本文根据安全且经济的原则 ,以及阻塞费用的计算规则,将各线路潮流值与潮流限值进行对比,调整各机组分配方案,建立了阻塞费用关于各机组出力的非线性模型,利用Lingo 软件求解,得到了相应的阻塞费用。

对问题五:将下一个时段预报的负荷需求改为MW 8.1052,将各线路潮流值与潮流限值比较,运用Lingo 软件求解,得到结果:不考虑潮流限值时,各组分配预案分别为:150、81、、、135、150、、117(单位:兆瓦),这时交易价格为MWh /356元。

考虑潮流限值时,得到各组分配预案分别为:153、88、228、、152、155、、117(单位:兆瓦),这时由于调整需要阻塞费用为元。

关键词:多元线性回归模型目标规划残差分析非线性规划模型.一、问题重述在电力市场化过程中,我国电力市场采取交易与调度一体化的模式。

从使电力从生产到使用的四大环节瞬间即可完成。

电力市场双边交易的阻塞管理

电力市场双边交易的阻塞管理

ma a e n ly n i o tn oei lcrct r e . I h sp p r h s fs n iii a e n DC p we lw s n g me tpa sa mp ra t l n ee tiiyma k t n t i a e ,te u e o e st t b s d o o rfo i r vy
算 中应用直流法的基础上 .提 出 了应 用灵敏度 的概 念解决 输 电阻塞问题 ,将 削减量 与所 占支路的容 量的 比例 结
合起采・做到公 正、佘平 的阻塞管理 。最后给 出了一个 5节点 系统的算俐 ,证明该方 法的有效性 。
关键词 :电力市场 ;双边 交易;潮流计算直流法;灵敏 度 系数 ;mg 管理 -
电 力市 场 双 边 交 易 的 阻塞 管理
李嘉龙 ,姚常青
(.广东省电力调度 中心 .广州 50 ( ;2 1 161 .天津大唐 国际盘 山发电有限责任公 司设备部 ,天津 3 10 ) 0 O9 0

要 :双边交 易是电力市场 中的一种 重要 交 易形式 ,其 相应 的 阻塞管理 也具 有重要 的意 叉。鍪此 .在潮流 计
P we nrt nC t F ni 1 90 hn ) o r Ge eai o ・L d - i j 3) 0 .C ia o a n 1
Ab ta t Bi t rl rd i n m p ra t r d mo e n o r ak t a d o r s o d n ta s sin o g sin src : l e a ta e s i o tn ta e a a d i p we m r e , n c rep n ig r n mi o c n e t s o

电力市场的输电阻塞优化管理

电力市场的输电阻塞优化管理

• 101•本文提出新的阻塞费用计算方法来解决目前电力市场出现的输电阻塞现象,利用实验数据建立各路线上有功潮流的多元线性回归模型,各机组出力分配模型和出力方案调整模型。

各模型均具有较强的通用性,能针对不同的实验数据、预测负荷需求,给出各机组出力分配预案的算法和出现输电阻塞后重新调整预案的模型。

当下,电力市场化将进入新一轮的发展,用电紧张缓解的市场形势也使有关产业及研究部门面临着新的挑战。

电力从生产到使用需要经过发电、输电、配电和用电共计四个步骤,并且这四个步骤是瞬间完成的。

电网公司在规划设计、施工安装时,必须遵守“安全第一”的原则。

考虑到线路走廊、电源点的分布、经济条件和社会基础的约束,电网要制订合理的规划,即实现购电费用最小的经济目标。

通常来说,输电阻塞是指制定的各机组出力分配方案使某条线路上有功潮流的绝对值超出其限值。

当出现这种现象时,电网亟需研究出既安全又经济的出力方案调整模型。

1 电网阻塞费用计算当出现输电阻塞,需要改变原有的各机组出力分配预案,一部分通过竞价获得发电权的发电容量(称序内容量)不能出力;而另一部分在竞价中未获得发电权的发电容量(称序外容量)要在低于对应报价的清算价上出力。

因此,发电商和网方将产生经济利益冲突,由此产生的费用称为阻塞费用。

网方应在电网安全运行的前提下同时尽可能减少阻塞费用。

阻塞费用包括两部分:(1)对不能出力或者少出力的序内容量的补偿;(2)对在低于对应报价的清算价上出力的序外容量的补偿,分别称之为欠发电补偿和过发电补偿。

1.1 序内容量不能出力的补偿设机组属于欠发电补偿情况,则其未出力的序内容量为 x i -x i ' ( x i ' ≤ x i )。

相对于方案调整前,机组少赚的钱为(x i -x i ' )×c 。

那么这些序内容量不能出力或者少出力导致发电商应得的利益损失的,为了公平合理,按照原价补偿,即欠发电补偿为:(1)其中x i 为第i 台发电机组的出力,x i '为第i 台发电机组在方案调整后的出力,c 为预案清算价。

电力市场的输电阻塞管理

电力市场的输电阻塞管理

电力市场的输电阻塞管理摘要输电线路的革新主要取决于输电和配电的工业化。

本文通过数据分析得到发电机组的负荷的改变量与各条线路上的潮流值的改变近似的成线性关系,利用题目中给定的数据计算出它们近似关系。

电网中心每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力。

电网每条线路上的有功潮流的绝对值有一安全限值,限值还具有一定的相对安全裕度(即在应急情况下潮流绝对值可以超过限值的百分比的上限)。

如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值,就造成了输电阻塞。

根据电力市场交易规则,有一种调配方案需要没有发电权的机组出力,电网公司在这种方案中需要支付阻塞费用,如何使阻塞费用达到最小,并且让竞争双方觉得公平,是要解决的中心问题,再就是如何运用我们所制定的阻塞费用规则。

任务一:通过对表1和表2的分析,看出了各个组的变化呈现出一定的规律,相对于方案0,这八个发电机组中的一个组的出力值变化四次,其余组保持出力值不变。

当一个组发生变化时,就会导致六条线路的有功潮流发生某些变化。

,总的负可以判断出是否发生了输电阻塞。

任务五:此问的分析方法同第三问和第四问相同二、模型假设1.第j台机组对第i条线路的改变量的影响是线性的2.假设竞价双方不存在市场权力的滥用,交易中心通过对竞争强度的控制,避免竞价电量被少数电厂所垄断的局面发生,同时也可以通过对竞争强度系数的适当调整,避免“平均主义”的情况发生。

3.为了使阻塞费用达到最小,则所选取得发电机组在那些没有发电权的机组中的报价应该是最小的。

4.假设爬破速率为均匀的。

即此过程中发电机组没有机器故障,输电网络也不存在故障。

x i0——根据电力市场交易规则确定的各机组出力改变量的分配计划y——根据电力市场交易规则确定的各条线路潮流改变量的分配计划i0——为阻塞费用的参照系数四、模型建立阻塞费用的计算规则:以偏差率的大小作为衡量阻塞费用的标准。

偏差率定义:实际改变量x i与计划改变量x i 0之差除以计划出力m x i i +0的比值的绝对值,补偿系数为m x x x i i i i +-00β,阻塞费用为m x x x ii i i +-00β()x x i i 0- 负荷档次划分规则:当各条线路上的有功潮流均达到限值时,各机组相对于0方案负荷最大变化总量为43.82MW ,此时各发电机组所能承受的最大无输电阻塞五、模型求解及结果分析对任务一的求解:由每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力,因此有功潮流与电机组的出力成一的关系可以设他们的关系式用f 来表示,为机组的出力量与线路的潮流量存在着线形的关系,因此令:根据全微分公式: 由假设知:ji x f ∂∂为常数,所以 根据泰勒公式的展开式得到:当第一台机组发生一次变化时,根据第一条潮流的变化,就可以算出11a ,因为第一台机组在其他组不变的情况下变了四次,所以就得到了四个11a ,为了ij a 形j 台发电 只考虑安全因素时,以线路中超过限值的最大百分比最小为目标,利用LINGO 软件进行求解:目标函数:min[max(s r y iii - 约束条件:t v x i i ≤仅考虑安全的最优结果为:只考虑经济因素时,以阻塞费用最小为目标。

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

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

电力市场的输电阻塞管理数学建模电力市场输电阻塞管理数学建模随着电力市场的发展和扩张,电力输电阻塞问题逐渐凸显出来。

输电阻塞是指电力系统中由于电网容量不足或者输电线路过载等原因导致的输电能力不足,无法满足电力需求的现象。

为了有效管理输电阻塞,避免电力市场的混乱和电力供需失衡,需要进行数学建模来优化输电系统。

输电阻塞管理的数学模型主要包括输电网络模型、输电能力模型和输电阻塞约束模型。

首先是输电网络模型。

电力系统可以看作是一个复杂的网络结构,输电线路和变电站等组成了这个网络。

输电网络模型将电力系统抽象为一个图,其中节点表示变电站,边表示输电线路。

每条边上有一个输电能力的限制,表示该线路所能承载的最大电流。

通过建立输电网络模型,可以清晰地描述电力系统的拓扑结构和电力流动情况。

接下来是输电能力模型。

输电线路的输电能力受到多种因素的影响,包括线路材料、长度、温度等。

输电能力模型通过考虑这些因素,计算出每条输电线路的实际输电能力。

这个模型可以帮助管理者了解每条线路的实际输电能力,为后续的优化决策提供依据。

最后是输电阻塞约束模型。

输电阻塞约束模型是整个数学建模的核心。

它将输电阻塞问题转化为一个优化问题,通过调整输电网络中各个节点的电压和相应的功率分配,使得电力系统达到最优状态,同时满足输电能力的限制条件。

这个模型可以通过数学优化方法求解,得到最佳的输电方案。

在实际应用中,输电阻塞管理数学建模可以应用于电力市场的日常运行和应急管理中。

通过建立准确的数学模型,可以对电力系统进行实时监测和调度,及时发现潜在的输电阻塞问题,并采取相应的措施进行调整和优化。

同时,在突发情况下,可以利用数学建模预测和分析输电阻塞的发生可能性,提前做好应急准备工作,保障电力市场的平稳运行。

电力市场输电阻塞管理数学建模是一项重要的工作,它可以帮助管理者更好地了解电力系统的输电能力和运行状态,预测和避免输电阻塞问题的发生,保障电力市场的稳定供应。

如何解析电力市场的输电阻塞管理

如何解析电力市场的输电阻塞管理

电力市场的输电阻塞管理摘要本模型根据下一时段的负荷预报、每台机组的报价、当前出力以及交易规则建立了一个分配预案模型。

当执行分配预案时,如果有线路出现输电阻塞,遵循电网“安全第一”的原则和输电阻塞管理原则建立了输电阻塞管理的优化模型。

就题目所给的五个问题本文建立了四个模型。

模型1通过一个物理学的公式,用差值的计算方法求出了有功潮流关于各机组出力的近似表达式,巧妙的避开了因改变量相对较小,误差不可忽略时而使拟和精度很低的情况,并通过一个数学公式给出了阻塞费用的计算规则。

在模型3中引入了一个定理,使得只要计算较少的负荷差的优化选取即可求得最优方案,同时根据爬坡速率的限制给出了可以选取的范围,大大减少了运算量。

对于模型4建立了一个优化模型,求解阻塞费用时,根据安全且经济的原则,建立了以阻塞费用、电网系统的稳定性和安全性为目标的多目标函数,但通过将电网系统的稳定性和安全性的约束转换为“损失的钱数”的约束,将问题最终转化为求各项“钱数”之和最少的单目标优化问题,有效的降低了问题求解的难度,并且模型具有很高的推广性。

针对MW的负荷需求,通过调整各机组的出力分配方案可以使输电阻塞消除,通过matlab编程求出最小的阻塞费用为元,但对于MW的负荷需求,无论怎样调整各机组的出力分配方案都不能消除输电阻塞,此时要使用线路的安全裕度输电,同样用matlab编程求出“钱数”之和最小为元。

一、问题重述电力从生产到使用的四大环节——发电、输电、配电和用电是瞬间完成的。

我国电力市场初期是发电侧电力市场,采取交易与调度一体化的模式。

电网公司在组织交易、调度和配送时,必须遵循电网“安全第一”的原则,同时要制订一个电力市场交易规则,按照购电费用最小的经济目标来运作。

设某电网有若干台发电机组和若干条主要线路,每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力。

电网每条线路上的有功潮流的绝对值有一安全限值,限值还具有一定的相对安全裕度。

电力市场的输电阻塞管理

电力市场的输电阻塞管理

电力市场的输电阻塞管理 This manuscript was revised by the office on December 10, 2020.电力市场的输电阻塞管理摘要输电线路的革新主要取决于输电和配电的工业化。

本文通过数据分析得到发电机组的负荷的改变量与各条线路上的潮流值的改变近似的成线性关系,利用题目中给定的数据计算出它们近似关系。

根据输电阻塞的限值和相对安全裕度,我们利用LINGO软件分别算出在输电阻塞限值条件和相对安全裕度下的各机组的最大负荷总量。

根据这两个值我们将负荷分为三个档次:阻塞可消除档次、安全限度内有阻塞档次和拉闸限电档次。

在安全限度有阻塞档次内,我们以安全和经济作为两个目标,先分别考虑这两个目标,然后将他们综合起来考虑,使这两个目标最优化。

我们还对输电阻塞的费用作了简单的化假设,用偏差率作为衡量阻塞费用的一个标准。

得到了如下的结果:以及清算值为356。

一、问题重述与分析电网中心每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力。

电网每条线路上的有功潮流的绝对值有一安全限值,限值还具有一定的相对安全裕度(即在应急情况下潮流绝对值可以超过限值的百分比的上限)。

如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值,就造成了输电阻塞。

根据电力市场交易规则,有一种调配方案需要没有发电权的机组出力,电网公司在这种方案中需要支付阻塞费用,如何使阻塞费用达到最小,并且让竞争双方觉得公平,是要解决的中心问题,再就是如何运用我们所制定的阻塞费用规则。

任务一:通过对表1和表2 的分析,看出了各个组的变化呈现出一定的规律,相对于方案0,这八个发电机组中的一个组的出力值变化四次,其余组保持出力值不变。

当一个组发生变化时,就会导致六条线路的有功潮流发生某些变化。

任务二:因为电网公司要满足安全又经济的原则,我们采用三步法,第一步只考虑安全问题,第二步只考虑经济问题,第三步是将安全和经济结合起来考虑制定了这种阻塞费用规则。

电力市场输电阻塞管理模型

电力市场输电阻塞管理模型

电力市场输电阻塞管理模型杨双红刘刚...本文对公平开放市场条件下,独市电网的输电阻塞管理问题做了模型研究。

首先,在局部线性化假设下,利用多元线性回归求取线路潮流分布与机组出力分配之间的近似公式。

本文对带有常数项和没有常数项的两种线性回归模型分别做了回归分析和细致的假设检验。

并山电力系统分析的背景知识,阐明了电网潮流分布与机组出力只有统计规律性,带有常数项的回归模型更合理。

根据阻塞调整产生的影响。

本文设计了“按损失成比例补偿”和“按市场规则确定费用”两种阻塞费用计算规则,并做了详细地比较讨论。

根据电力市场交易规则,兼顾计算的时间效率,本文利用递归策略给出了简单易行的出力分配预案计算方法及其流程图,在机组数不多时简单的手工计算很容易求得分配顶案。

对阻塞调整问题,本文按电网“安全第,兼顾经济”的原则,提出分阶段(共分四个阶段:阻塞检查、调整预案、裕度输电、拉闸限电)按步骤规划的计算流程,并对各个规划阶段分别建立了数学模型:阻塞检查为判断一组不等式;调整预案是求解以阻塞费用最小为目标的规划问题;裕度输电规划先以裕度占用率最小为目标,再在裕度占用率不增的条件下以阻塞费用最小为目标做规划;拉闸限电规划则是在保证电网最低安全水平的条件下,以总出力最大为目标做规划。

化简后,各阶段的规划模型,除调整预案模型是线性约束条件下阶梯函数族的最大最小规划外,其余阶段规划模型均为线性规划。

出于计算效率的考虑,结合题目特点,本文发现以Huffman树作为决策树时,阻塞管理问题的规划流程具有最高计算效率,此时通过对几条简单的规则的判定即可确定应该进行哪个阶段的规划,从而不必一步步按部就班地进行。

本文还对Huffmnan决策树流程的一些技术细节及改进节点定义的最优流程做了讨论。

另外,本文从广义函数角度对阶梯函数的数学分析性能及优化解法做了讨论,并给出了求解以阶梯函数为目标的优化问题的求解建议及两种简单易行的启发式算法,并在附录中,给出了些典型算法的流程图。

2004年全国数学建模 B题 电力市场的阻塞管理

2004年全国数学建模 B题  电力市场的阻塞管理

电力市场的输电阻塞管理(CUMCM 2004B)B10020614 孙朋B10020609 潘姝亚B10020701 王琪2013-08-10 NUPT****************摘要本题是关于电力市场的输电管理问题,首先要求附件中的数据确定发电机组出力与输电线路潮流值的关系,然后根据电力市场交易调度规则建立适当的模型,求解给定负荷的发电机组出力方案,并检测方案是否会引起输电阻塞。

当发生阻塞时考虑不同因素,在安全且经济的原则下,给出调整后的方案。

问题一:要求结合附件表1和表2的数据,确定8台机组的出力和6条线路的潮流值的关系式。

首先,我们利用散点图观察二者的定性关系,发现线路潮流值与机组出力值大致成线性关系。

用SPSS软件计算各机组出力的Pearson相关系数,结果表明各机组出力之间线性无关。

因此尝试对输电线路与机组出力进行多元线性回归分析,用MATLAB编程求出回归系数。

最后,利用统计检验回归方程(见表3),结果显示回归方程显著说明可行。

问题二:要求根据电力市场规则设计一种阻塞费用计算规则。

阻塞费用包括补偿不能出力的序内容量和补偿在低于对应报价的清算价上出力的序外容量。

发电商的损失与其减少的出力量、序内容量的报价及清算价格有关。

出于简化问题的考虑,我们认为网方给予发电商阻塞补偿即为输电阻塞所导致发电商的损失。

对于某一机组有确定的出力时,其各个段容的出力是容易确定的。

结合调整后的方案,分段容计算阻塞费用,然后求和,最后对各个机组的阻塞费用进行求和即可得到阻塞费用。

问题三:要求在已知下一个时段预报的负荷需求为982.4MW的情况下,按照电力市场规则,同时结合各机组段容量、段价和爬坡速率等因素,给出下一个时段各机组的出力分配预案。

这是一个最优化问题,目标函数为最小购电费用。

电力市场规则规定确定的负荷有一个统一的价格即清算价,同时考虑各机组爬坡速率的限制等因素确定约束条件。

利用LINGO软件编程求解,最后得到分配方案如表6所示。

电力市场的输电阻塞管理论文

电力市场的输电阻塞管理论文

论文2电力市场的输电阻塞管理参赛臥员:、、幫要:本文根据电力市场的交易规此就目前我国电力系筑中各个发电机组的出力分配硕案和各线路的有功潮渝冋题进行了深人分析,并对产生输电皿塞的分配预案进行了调整, 得到了较好的岀力分配方案。

1•根据各机组岀力和各线路潮流的关系建立了一个多元线性回旧模塑(见模塑一),利用所给实验数振呆用最小二乘法6I归,得到每条线路上的潮流值关于各发电机组出力的的ifi lfl表这式,并对毎一个表这式进11 Til差分ffi, II岀各表达式的复相关系数,可以看岀我们的回归效果显著,说明我们的模型是可靠、合理的。

2. 我们果用pool模式下的输电81塞贾用廿算方法,公平对待序内序外两fttiift.aif 岀了一种简明、合理的阻塞费用廿算规!M:第一、果用序外多发电量按18发电报价廿算; 第二、序内少发电量按涓算价与发电报价之差价给算。

并建立了一个合理的计算皿塞费用模型。

3. °在下一时目预报负荷需求为982.4MW tf]条件下,根据市场规则,以最小即电费用为目标、以机组的段容量,爬玻速率作为约束条件,采用动态规划算法建立了一f单目标规划模型,通过数学軟件M ATLAB编程给出各机组的岀力分配预案,各台机组的出力分别为(MW): 150、79、180、99.5、125、140、95、113.9。

总费用为:C. =74416.8 元。

清算价为:g = 3O3S/MWh4•通过对预案分析it算可得,第一、五、X线路岀现输电阻塞现象,根稠安全且经济的原则,利用排序算法进行了调整,得到了消除输电皿塞的分配方案,分别是:152.5947,8&0,228.0,81.7434,152.0,9&0,65.0592,117.0。

其清算价怡Q =303 元/MWh 阻塞费用:C = 12026.0元;总费用为:C, =86140.8元5.同理对下一时段硕报负荷需求为1052.8MW的条件下,重复步骤3、4的工作,借到分配预案为(MW):152.0,81.0,218.2,99.5,135.0,150.0,102.1,117.0.总费用为:Q =93699.2元。

基于改进粒子群优化算法的阻塞管理

基于改进粒子群优化算法的阻塞管理

2 模型
在 PO O L模式下 , 台竞 价机组提交非减阶梯状 每
的报 价 曲线 。运 行 人 员 根 据 电 力 系 统 运 行 的约 束 条 件 , 化 系统 中发 电机 的 出力 , 总 发 电成 本最 小。其 优 使
目标 函数如 下 :
效用取代购 电费用最小作为 区域 阻塞管理 的 目标 函 数 , 出了全新 的基 于风 险 的 区 域 阻 塞管 理 算 法 。文 提
SU n Big
( an g i n g o e u p & ti tnE g er gC . Ld N nig 30 2 C i ) N n i a i w r p l Uiz i ni e n o ,t, a n 0 2 , h a n jn n P S y la o n i n5 n
维普资讯
<电气开关》 20 . o5 (0 8 N . )
2 9
文章 编号 :0 4— 8 X( 0 8 O 0 2 o 10 2 9 20 )5— 0 9一 3
基 于 改 进 粒 子群 优 化 算 法 的阻 塞 管 理
苏冰
( 南宁建 宁供 用电工程有 限责任公 司 , 西 南宁 5o2 ) 广 302
摘 要 : 电力 市场环境 下, 在 阻塞 管理是 确 保 电力 系统安全 经 济运行 的 重要 途 径 。针 对 电力 系统 阻塞管理 问题 的 求解 , 把混沌寻优思想引入到粒子群优化算法 中, 出一种改进粒子群算法, 提 解决 了基本粒子群算法易陷入局部
极小点的缺陷, 从而提高了收敛精度。算例结果表明了提 出的算法相对于其他算法的有效性与可行性。 关键词 : 电力系统; 阻塞管理; 改进粒子群算法 中图分 类号 :M 1 T 71 文献 标 识码 : B Co g sin M a a e e tBa e n I r v d P r il wa ms Op i z t n n e to n g m n s d o mp o e a tce S r t mi a i o

电力市场输电阻塞管理问题的规划模型求解[1]

电力市场输电阻塞管理问题的规划模型求解[1]

第35卷第5期2005年5月数学的实践与认识M A TH EM A T I CS I N PRA CT I CE AND TH EO R YV o l 135 N o 15 M ay, 2005 建 模电力市场输电阻塞管理问题的规划模型求解陈以平(湖北民族学院理学院,湖北恩施 445000)摘要: 针对CUM CM 2004B 电力市场的输电阻塞管理问题,建立目标规划模型,详细给出利用M A TLAB优化工具箱函数linp rog 及fgoalattain 求解模型的方法,指出正确使用数学软件在建模活动中的重要性.关键词: 输电阻塞管理;目标规划模型;M A TLAB 优化工具箱;linp rog ;fgoalattain1 引 言收稿日期:2004211208基金项目:湖北省教育厅优秀中青年项目(Q 200529001) 在用数学建模的方法解决实际问题时,往往需要通过计算机解决大量的计算,尤其在数学建模竞赛中,由于时间短,任务重,参赛选手是否具备熟练应用计算机高级语言编程求解的能力,从某种意义上说,是决定能否成功参赛的关键.M A TLAB 以其简洁的语言、方便的操作以及开放的体系,在数值计算、数据分析、图形显示、应用编程等方面得到了广泛的接受和应用,因此也就成为了数学建模教学和赛前培训的重要内容.本文以全国大学生数学建模竞赛(CUM C M )2004B 电力市场的输电阻塞管理问题为例,详细讨论如何利用M A TLAB 优化工具箱相关函数解决规划问题.原题目详见[1],略.本文所述表1~6与赛题表格序号对应.2 记号说明b 1ij :第i 个方案第j 台机组出力值(对应表1中方案1~32),记b 1=(b 1ij )32×8;l i :第i 台机组当前时段出力值(对应表1中方案0),记L =(l i )T1×8;b 2ij :第i 个方案第j 条线路潮流值(对应表2中方案1~32),记b 2=(b 2ij )32×6;b 3ij :第i 台机组第j 个序段的段容量值(对应表3),记b 3=(b 3ij )8×10;b 4ij :第i 台机组第j 个序段的段价(对应表4),记b 4=(b 4ij )8×10;b 5i :第i 台机组的爬坡速率(对应表5),记b 5=(b 5i )T 1×8;b 6ij :对应表6,b 61j 为第j 条线路的潮流限值,b 62j 为第j 条线路的相对安全裕度,记b 6=(b 6ij )2×6;b ij :第i 台机组第j 个序段的出力系数(段容量使用度因子)(0Φb ij Φ1),记b =(b ij )8×10;Y b :下一时段的负荷预报值;x i :第i 台机组的出力;y i :第i 条线路上的有功潮流.3 模型的建立根据市场交易规则及文献[2—6]相关知识,建立模型,过程分析略.311 电网中各线路有功潮流关于发电机组出力的数学模型y i =a i 0+∑8j =1aijx j , i =1,2,…,6(1) 记a =(a i 0)1×6,A =(a ij )6×8,A 为第j 台机组对第i 条线路出力的系数矩阵,a 为常数向量.注1 理论与实验研究表明电网中各主要线路有功潮流值与各机组出力之间具有近似线性叠加关系[3],据此建立模型(1),该模型实际上是一个多元线性回归模型.312 下一时段电网中各机组出力分配预案的规划模型m in∑8i =1∑10j =1(b ijb 3ijb 4ij)◊购电费用最小s .t -15b 5i Φ∑10j =1b ijb3ij-l i Φ15b 5i , i =1,2,…,8◊机组爬坡速率约束 ∑8i =1∑10j =1b ij b 3ij =Y b◊系统负荷平衡约束 0Φb ij Φ1, i =1,2,…,8,j =1,2,…,10(2) 注2 按照市场交易-调度中心在当前时段内要完成的具体操作过程,利用模型(2)求出下一时段电网中各机组每个序段的出力系数b ij ,再计算出各机组出力分配预案,然后利用模型(1)计算当执行各机组出力分配预案时电网各主要线路上的有功潮流,据此判断是否会出现输电阻塞.注3 各机组出力计算公式为:∑10j =1b ijb 3ij,其中i 表示机组序号.313 调整电网中各机组出力分配方案的规划模型m in∑8i =1∑10j =1(b ijb 3ijb4ij)◊购电费用最小s .t -15b 5i Φ∑10j =1b ijb3ij-l i Φ15b 5i ,i =1,2,…,8◊机组爬坡速率约束 ∑8j =1a ij∑10k =1bjkb 3jk +a i 0Φb 61i , i =1,2,…,6◊线路潮流限值约束 ∑8i =1∑10j =1b ij b 3ij =Y b◊系统负荷平衡约束 b i 1b 3i 1=b 3i 1, i =1,2,…,8◊最小技术出力约束 0Φb ij Φ1, i =1,2,…,8(3)2数 学 的 实 践 与 认 识35卷m in ∑8i =1∑10j =1(b ijb3ijb 4ij )◊购电费用最小m in∑8j =1a ij∑10k =1bjkb 3jk +a i 0-b 61ib 61i , i =1,…,6◊各线路潮流超限值比例最小s .t -15b 5i Φ∑10j =1b ijb3ij-l i Φ15b 5i , i =1,2,…,8◊机组爬坡速率约束∑8j =1a ij∑10k =1bjkb 3jk +a i 0<b 61i +b 61i b 62i , i =1,…,6◊相对安全裕度约束 ∑8i =1∑10j =1b ij b 3ij =Y b◊系统负荷平衡约束 b i 1b 3i 1=b 3i 1, i =1,2,…,8◊最小技术出力约束 0Φb ij Φ1, i =1,2,…,8,j =1,2,…,10(4) 注4 模型(3)、(4)是根据输电阻塞管理原则,在模型(2)的基础上,通过增加约束条件及目标函数得到的(最小技术出力约束可改为机组出力上下限约束,这里不考虑阻塞费用作目标函数).首先利用模型(3)计算,若无可行解(说明不能消除输电阻塞),则用模型(4),若模型(4)无可行解,则只能通过用电侧拉闸限电解决输电阻塞问题.拉闸限电量的计算只需将模型(4)中的系统负荷平衡约束条件修改为一个目标函数,该目标函数为各线路出力之和的相反数最小,然后利用该模型求出不出现输电阻塞的各机组最大出力,再用负荷预报值减各机组最大出力之和即得.以上所列的几个规划模型含有80个变量,近103个约束条件,手工计算不可能,只能借助计算机求解.虽然在处理线性、非线性规划等问题时,采用SA S 及L ingo 等软件要方便些,但在数学建模的教学和赛前培训中,多数介绍的是M A TLAB 软件的使用.事实上,利用M A TLAB 优化工具箱相关函数解决规划问题亦十分方便,不仅编程效率高,计算功能强,而且使用简便.4 模型的求解这里,我们以M A TLAB 611为运行环境[7],详细介绍以上几个模型的求解.411 原始数据的录入利用M A TLAB 强大的数据导入功能,先将问题中已知数据表格1~6转化为矩阵变量数据格式(3.m at )的文件,保存在磁盘中,以备调用.按前述记号,文件名依次为:b 1、L 、b 2、b 3、b 4、b 5、b 6.具体做法是:利用从h ttp : m c m .edu .cn m c m 04 P rob lem s 2004c .h tm 下载的w o rd 文档,分别选定有关表格中的数据,通过复制、粘贴,分别建立这些数据的文本格式(3.tx t )的文件,然后打开m atlab 611,从文件下拉菜单中单击I m po rt D ata 打开包含前面建立的文件的窗口,选定拟导入的某个数据文本文件,打开I m po rt W izard 窗口,按提示完成操作,即得选定文本文件中的数据对应的矩阵变量数据格式(3.m at )的文件,打开w o rk sp ace ,选定相应文件存盘.依此操作,即将已知数据表格1~6全部转化为矩阵变量数据格式(3.m at )的文件保存.35期陈以平:电力市场输电阻塞管理问题的规划模型求解注5 按以上操作方法,大约3分钟时间,即可完成原始数据的录入,生成矩阵变量,不仅数据准确,而且节约了大量时间.在建模软件教学中,有针对性地加强培训十分必要.412 优化工具箱中求解目标规划模型的函数l i nprog 、fgoa la tta i n 用法介绍[7]4.2.1 l i nprog 函数功能:求解单目标线性规划问题数学模型:m in xf T xs .t A x Φb A eq x =beq lb Φx Φub式中f ,x ,b ,beq ,lb 和ub 为向量,A 和A eq 为矩阵.常用的调用格式:[x ,f va l ,ex itf lag ]=linp rog (f ,A ,b ,A eq ,beq ,lb ,ub ) 其中x 记录返回的最小值点(即所求的解),f va l 记录解x 处的目标函数最小值,ex itf lag 记录程序退出条件,取值为1,表示运行成功,取值为0,表示未达到给定精度便已超出最大迭代次数,取值为-1,表示无可行解或解无界.41212 fgoa la tta i n 函数功能:求解多目标达到问题数学模型:m in x ,ΚΧF (x )-w eig h t ΧΦg oa lc (x )Φ0ceq (x )=0A x Φb A eq x =beq lb Φx Φub式中,x ,w eig h t ,g oa l ,b ,beq ,lb 和ub 为向量,A 和A eq 为矩阵,c (x ),ceq (x )和F (x )为函数,返回向量.F (x ),c (x )和ceq (x )可以是非线性函数.注意F (x )为目标函数向量.常用的调用格式:[x ,f va l ,a tta inf actor ,ex itf lag ]=f g oa la tta in (f un ,x 0,g oa l ,w eig h t ,A ,b ,A eq ,beq ,lb ,ub )其中x 记录返回的最小值点(即所求的解),f va l 记录解x 处的目标函数最小值;a tta inf actor 记录解x 处的目标达到因子,若a tta inf actor 为负,则目标已经溢出,若a tta inf actor 为正则目标个数还未达到;ex itf lag 记录程序退出条件.f g oa la tta in 函数调用的详细说明参见[7].413 模型的向量、矩阵表示41311 补充记号说明(为方便,部分记号使用M A TLAB 语言说明)x :矩阵b 中的元素按行顺序连接生成的80维列向量.B :各机组的段容量(表3)与其对应的段价(表4)相乘构成的8×10阶矩阵(即b 3.3b 4)中的元素按行顺序连接生成的80维行向量.C :矩阵b 3中的元素按行顺序连接生成的80维行向量.令e =z eros (1,10),用D 1表示如下矩阵:4数 学 的 实 践 与 认 识35卷D 1=[b 3(1,:),e ,e ,e ,e ,e ,e ,e ;e ,b 3(2,:),e ,e ,e ,e ,e ,e ;e ,e ,b 3(3,:),e ,e ,e ,e ,e ;e ,e ,e ,b 3(4,:),e ,e ,e ,e ;e ,e ,e ,e ,b 3(5,:),e ,e ,e ;e ,e ,e ,e ,e ,b 3(6,:),e ,e ;e ,e ,e ,e ,e ,e ,b 3(7,:),e ;e ,e ,e ,e ,e ,e ,e ,b 3(8,:)] 记D =[D 1;-D 1];E =[L +153b 5;-L +153b 5];G =[b 6(1,:)T -a ;b 6(1,:)T +a ];q =b 3(:,1);H =[b 6(1,:)T +b 6(1,:)T .3b 6(2,:)T -a ;b 6(1,:)T +b 6(1,:)T .3b 6(2,:)T+a ].I i (i =1,2,…,6):矩阵A 的第i 行第j 个元素(j =1,2,…,8)与矩阵b 3的第j 行元素(j =1,2,…,8)逐个相乘所得新矩阵按行顺序连接生成的80维行向量,即I i =[a i 1b 311,a i 1b 312,…,a i 1b 3110,a i 2b 321,a i 2b 322,…,a i 2b 3210,…,a i 8b 381,…,a i 8b 3810]J :以I i 为第i 行生成的6×80维矩阵,记N =[J ;-J ].Q :表示第i 行第103i -9列元素为b 3i 1(i =1,2,…,8),其它元素为0的8×80维矩阵.41312 模型的向量、矩阵表示利用前述有关记号,模型(2)、(3)、(4)可分别改写为:m in B x s .t D x ΦE C x =Y b 0Φx Φ1; m in B xs .t D x ΦE N x ΦG C x =Y b Q x =q 0Φx Φ1; m in B xm in ( I i x +a i 0 -b 61i )b 61i , i =1,2,…,6s .t D x ΦE N x ΦH C x =Y b Q x =q 0Φx Φ1 注6 能否应用向量和矩阵正确表示模型是用M A TLAB 优化工具箱相关函数解决规划问题的关键.414 程序实现4.4.1 模型(1)的程序实现可利用函数regress 求解模型(1),假设以95%为显著性水平.先编制M 文件m odel -1.m 如下:%m odel -1.mclear ,load b 1,load b 2,b 11=[ones (32,1)b 1];fo r i =1:6[b ,b in t ,r ,rin t ,stats ]=regress (b 2(:,i ),b 11);stats 1(i ,:)=stats ;a (i )=b (1);A (i ,:)=b (2:9)′;enda =a ′;a ,A ,fo rm at long ,stats 1,fo rm at sho rt ,save A 1a A 运行m odel -1求得模型(1)中的回归系数及其检验结果,见下表.55期陈以平:电力市场输电阻塞管理问题的规划模型求解6数 学 的 实 践 与 认 识35卷表1 模型(1)中的回归系数及其检验结果线路a ij(i=1,…,6;j=1,…,8)a i0R2F p序号10.0828,0.0483,0.0530,0.1199,-0.0254,0.1220,0.1216,-0.0012110.29650.99955861.502-0.0546,0.1279,-0.0000,0.0333,0.0868,-0.1124,-0.0189,0.0987131.22890.99967228.70 3-0.0695,0.0616,-0.1566,-0.0099,0.1245,0.0021,-0.0025,-0.2014-108.87320.9999223520 4-0.0345,-0.1024,0.2052,-0.0208,-0.0118,0.0060,0.1449,0.076577.48170.999925583050.0005,0.2433,-0.0646,-0.0411,-0.0652,0.0703,-0.0043,-0.0089132.97450.99966971.8060.2378,-0.0602,-0.0779,0.0930,0.0469,0.0001,0.1659,0.0007120.66330.9998174550注7 从上表1中检验结果可看出:相关系数平方值均近似于110,说明拟合程度非常好;对多元线性回归作方差分析,F值均大于置信限F0.05(8,32-8-1),回归显著,线性相关密切;显著性概率p均为0,小于0105,因此拒绝零假设,回归方程中至少有一个系数不为零,回归方程有意义.另外,对残差r的正态分布检验,接受残差r服从正态分布的假设.因此,用多元线性回归模型拟合得非常好.注意程序中用save A1a A将a及A存盘,以备调用.41412 模型(2)的程序实现模型(2)属单目标线性优化问题,我们用函数linp rog求解.首先编制子程序m odel-21.m以备后面模型求解调用.%m odel-21.mload b3,load b4,load b5,load b6,load L,load A1B1=b3.3b4;B=[B1(1,:),B1(2,:),B1(3,:),B1(4,:),B1(5,:),B1(6,:), B1(7,:),B1(8,:)]′;save B BC=[b3(1,:),b3(2,:),b3(3,:),b3(4,:),b3(5,:),b3(6,:),b3(7,:),b3(8,:)];save C Ce=zero s(1,10);D1=[b3(1,:),e,e,e,e,e,e,e;e,b3(2,:),e,e,e,e,e,e;e,e,b3(3,:),e,e,e,e,e;e,e, e,b3(4,:),e,e,e,e;e,e,e,e,b3(5,:),e,e,e;e,e,e,e,e,b3(6,:),e,e;e,e,e,e,e,e,b3(7,:),e;e,e,e,e,e,e,e,b3(8,:)];D=[D1;-D1];E=[L+153b5;-L+153b5];lb=zero s(80,1);ub=ones(80,1);接下来编制模型(2)的求解程序m odel-22.m%m odel-22.m,求出力分配预案使用clear,m odel-21,Yb=inp u t(‘请输入下一时段负荷预报值’)[x,fval,ex itflag]=linp rog(B,D,E,C,Yb,lb,ub);M=x′;disp(‘以各机组每个序段的出力系数为元素的矩阵b为’)b=[M([1:10]);M([11:20]);M([21:30]);M([31:40]);M([41:50]);M([51: 60]);M([61:70]);M([71:80])]fval,ex itflag若下一时段负荷预报值98214MW,运行程序m odel-22得:µm odel -22请输入下一时段负荷预报值982.4Yb = 982.4000O p ti m izati on term inated successfu lly .以各机组每个序段的出力系数为元素的矩阵b 为b = 1.000001.0000001.0000000.00001.000001.00001.00001.00001.00000.0000000.00001.000001.000001.000000.00000.000000.00001.00001.00001.00001.00001.00000.95000.000000.00001.00001.00001.000001.00001.000000.00000.00000.00001.00001.0001.000001.00000.00000.000000.00001.00001.00001.00001.00001.00000.00000.00000.00000.00000.00001.000001.000001.000000.19500.00000.00000.0000fval = -2.4665e +005ex itflag = 1由于ex itflag =1表示迭代收敛,得到最优解.显然,第4台机组的第6个段容量由于爬坡速率约束,只选取了其中的95%,最后一个被选取的段容量为第8台机组的第7个段容量,且只选取了其中的1915%.由此得该时段的清算价为:303元.根据机组出力计算公式为:∑10j =1b ijb3ij(其中i 表示机组序号)及模型(1),编制子程序m odel -23.m :%m odel -23.mFA 1=zero s (1,80);fo r i =1:80 FA 1(i )=C (i )3M (i );endFA =zero s (1,8);fo r i =1:10:80j =(i -1)10+1;FA (j )=sum (FA 1(i :i +9));enddisp (‘机组出力分配方案为’),FA disp (‘执行此出力分配方案时电网中各线路的有功潮流值为’),CL =(A 3FA ′+a )′然后在程序m odel -22后面增加语句m odel -23,并存盘,再次运行m odel -22,即可求得当下一时段负荷预报值为98214MW 时各机组的出力分配预案为:150,79,180,9915,125,75期陈以平:电力市场输电阻塞管理问题的规划模型求解8数 学 的 实 践 与 认 识35卷140,95,113.9;当执行此分配预案时电网中各线路的有功潮流值为:17313164,14110129, -15019312,12019198,13618368,16815309.对照题中表6给出的各线路潮流限值及相对安全裕度,显然除第2、3、4条线路外,其余线路潮流绝对值均超出限值,出现输电阻塞,但仍然在相对安全裕度范围内.下一步通过模型(3)调整出力分配方案,看是否能在潮流限值范围内消除输电阻塞.仿此可讨论负荷预报值为1052.8MW时的情况.4.4.3 模型(3)、(4)的程序实现模型(3)亦属单目标线性优化问题,首先编制子程序m odel-31.m以备后面模型求解调用.%m odel-31.mJ=zero s(6,80);I=zero s(8,10);fo r f=1:6 fo r i=1:8 fo r j=1:10 I(i,j)=A(f,i)3b3(i,j); end end J(f,:)=[I(1,:),I(2,:),I(3,:),I(4,:),I(5,:),I(6,:),I(7,:),I(8,:)];endN=[J;-J];G=[b6(1,:)′-a;b6(1,:)′+a];Q=zero s(8,80);fo r i=1:8 Q(i,103i-9)=b3(i,1);end接下来编制模型(3)的求解程序m odel-32.m%m odel-32.m,对应输电阻塞管理原则1clear,m odel-21,m odel-31,Yb=inp u t(‘请输入下一时段负荷预报值’)F2=[D;N];G1=[E;G];C1=[C;Q];Y1=[Yb;b3(:,1)];[x,fval,ex itflag]=linp rog(B,F2,G1,C1,Y1,lb,ub);M=x′;disp(‘以各机组每个序段的出力系数为元素的矩阵b为’)b=[M([1:10]);M([11:20]);M([21:30]);M([31:40]);M([41:50]);M([51: 60]);M([61:70]);M([71:80])]fval,ex itflag,m odel-23取下一时段负荷预报值98214MW,运行程序m odel-32,有ex itflag=1,表示迭代收敛,得到最优解,说明应用模型(3)(对应输电阻塞管理原则1)调整机组出力方案后,能在潮流限值范围内消除输电阻塞,具体求解结果运行程序m odel-32即得,这里限于篇幅略.但当取下一时段负荷预报值105218MW,运行程序m odel-32,ex itflag=-1,无可行解,仍出现输电阻塞.下面利用模型(4)判断是否能消除输电阻塞.模型(4)属多目标达到问题,需用函数fgoalattain求解,其求解程序m odel-4.m如下:%m odel-4.m,对应输电阻塞管理原则2clear ,m odel -21,m odel -31,Yb =inp u t (‘请输入下一时段负荷预报值’)H =[b 6(1,:)′+b 6(1,:)′.3b 6(2,:)′-a ;b 6(1,:)′+b 6(1,:)′.3b 6(2,:)′+a ];F 2=[D ;N ];G 1=[E ;H ];C 1=[C ;Q ];Y 1=[Yb ;b 3(:,1)];goal =[-11230.131 230.181 230.091 230.111 230.151 230.14];w eigh t =[11230.131 230.181 230.091 230.111 230.151 230.14];x 0=ones (80,1);[x ,fval ,attainfacto r ,ex itflag ,ou tp u t ]=fgoalattain (@op tq -4,x 0,goal ,w eigh t ,F 2,G 1,C 1,Y 1,lb ,ub );fo r i =1:2 x 0=x ;[x ,fval ,attainfacto r ,ex itflag ,ou tp u t ]=fgoalattain (@op tq -4,x 0,goal ,w eigh t ,F 2,G 1,C 1,Y 1,lb ,ub );endM =x ′;disp (‘以各机组每个序段的出力系数为元素的矩阵b 为’)b =[M ([1:10]);M ([11:20]);M ([21:30]);M ([31:40]);M ([41:50]);M ([51:60]);M ([61:70]);M ([71:80])]fval ,attainfacto r ,ex itflag ,ou tp u t ,m odel -23其中目标函数op tq -4的M 文件为:%op tq -4.mfuncti on f =op tq -4(x )load B ,load J ,load b 6,load A 1f (1)=B ′3x ;f (2)=(ab s (J (1,:)3x +a (1))-b 6(1,1)) b 6(1,1);f (3)=(ab s (J (2,:)3x +a (2))-b 6(1,2)) b 6(1,2);f (4)=(ab s (J (3,:)3x +a (3))-b 6(1,3)) b 6(1,3);f (5)=(ab s (J (4,:)3x +a (4))-b 6(1,4)) b 6(1,4);f (6)=(ab s (J (5,:)3x +a (5))-b 6(1,5)) b 6(1,5);f (7)=(ab s (J (6,:)3x +a (6))-b 6(1,6)) b 6(1,6);当取下一时段负荷预报值1052.8MW ,运行程序m odel -4,得ex itflag =1,表示过程成功收敛于解x 解处.说明应用模型(4)(对应输电阻塞管理原则2)调整机组出力方案后,能在安全裕度范围内输电.调整后的各机组的出力分配方案为:153,88,228,9915,152,12116785,9316215,117;当执行此分配方案时电网中各线路的有功潮流值为:17314483,14617357,-1551401,12913518,13218578,16610033;利用表6给出的各线路潮流限值及相对安全裕度,可求得当执行此分配方案时各线路有功潮流绝对值与线路极限负荷之差为-1310017,-3012643,-1819990,-4216982,-1819422,-1816767.显然此时输电阻塞消除.注8 若当下一时段给出的某负荷预报值,导致模型(4)无可行解,则可判断此时无法消除输电阻塞,只能在用电侧拉闸限电.根据前面介绍的拉闸限电量计算模型(见注4),相应修改目标函数op tq -4和程序m odel -4.m ,不难求出拉闸限电量.另外,由于模型(4)属多目标95期陈以平:电力市场输电阻塞管理问题的规划模型求解01数 学 的 实 践 与 认 识35卷规划模型,目前虽有一些较成熟的求解算法,但对大型问题一般计算速度慢,且会出现未达到给定精度便已超出最大迭代次数的情况(此时ex itflag=0),另外特别是对一些实际上根本无可行解的情况,更费时(如:我们在建立模型(1)时若不考虑常数项,相应修改模型(4)的目标函数及安全裕度约束条件,再利用该模型编程计算,若取负荷预报值为105218,则有ex itflag=0,且很费时,而事实上,这时模型(4)无可行解),针对这种情况,可通过减少目标函数,借助辅助模型来解决问题.参考文献:[1] Hom e Page of CUM CM.赞助单位2本站.h ttp: m c m c m04 P roblem s2004c.h tm[2] 柯进,管霖.电力市场下的输电阻塞管理技术[J].电力系统自动化,2002,26(14):20—24.[3] 王锡凡,王秀丽,陈皓勇.电力市场基础[M].西安:西安交通大学出版社,2003.[4] 袁辉等.电力市场中爬坡约束问题的一种实用算法[J].电力系统自动化,2001,25(12):17—19.[5] 高旭矗,伍永刚.考虑机组爬坡率的火电系统日发电计划编制[J].华中电力,2002,15(2):1—2.[6] 谷源盛.运筹学[M].重庆:重庆大学出版社,2001.[7] 苏金明,阮沈勇.M A TLAB611实用指南(下册)[M].北京:电子工业出版社,2002.The Solution to the Programm i ng M odel of Tran s m ission Congestion M anagem en t i n Power M arketCH EN Y i2p ing(Co llege of Science,H ubei In stitu te fo r N ati onalities,En sh i445000,Ch ina)Abstract: T he p rogramm ing models are con structed concern ing the m anagem en t oftran s m issi on congesti on in pow er m arket in CUM C M2004.A nd an app roach is p ropo sed toso lve the model in tail by u sing the functi on linp rog and fgoalattain in M A TLAB op ti m izati ontoo lbox.In additi on,the exact app licati on of m athem atical softw are is a key in them athem atical modeling.Keywords: tran s m issi on congesti on m anagem en t;ob jective p rogramm ing model;M A TLABop ti m izati on too lbox;linp rog;fgoalattain。

电力市场的输电阻塞管理

电力市场的输电阻塞管理

电力市场的输电阻塞管理摘要输电线路的革新主要取决于输电和配电的工业化。

本文通过数据分析得到发电机组的负荷的改变量与各条线路上的潮流值的改变近似的成线性关系,利用题目中给定的数据计算出它们近似关系。

根据输电阻塞的限值和相对安全裕度,我们利用LINGO软件分别算出在输电阻塞限值条件和相对安全裕度下的各机组的最大负荷总量。

根据这两个值我们将负荷分为三个档次:阻塞可消除档次、安全限度内有阻塞档次和拉闸限电档次。

在安全限度有阻塞档次内,我们以安全和经济作为两个目标,先分别考虑这两个目标,然后将他们综合起来考虑,使这两个目标最优化。

我们还对输电阻塞的费用作了简单的化假设,用偏差率作为衡量阻塞费用的一个标准。

得到了如下的结果:以及清算值为356。

一、问题重述与分析电网中心每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力。

电网每条线路上的有功潮流的绝对值有一安全限值,限值还具有一定的相对安全裕度(即在应急情况下潮流绝对值可以超过限值的百分比的上限)。

如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值,就造成了输电阻塞。

根据电力市场交易规则,有一种调配方案需要没有发电权的机组出力,电网公司在这种方案中需要支付阻塞费用,如何使阻塞费用达到最小,并且让竞争双方觉得公平,是要解决的中心问题,再就是如何运用我们所制定的阻塞费用规则。

任务一:通过对表1和表2 的分析,看出了各个组的变化呈现出一定的规律,相对于方案0,这八个发电机组中的一个组的出力值变化四次,其余组保持出力值不变。

当一个组发生变化时,就会导致六条线路的有功潮流发生某些变化。

任务二:因为电网公司要满足安全又经济的原则,我们采用三步法,第一步只考虑安全问题,第二步只考虑经济问题,第三步是将安全和经济结合起来考虑制定了这种阻塞费用规则。

任务三:当前时段各机组负荷量总和为KW1.4.982,874,而下个时段的负荷需求为KW总的负荷量改变了KW108,因为负荷量发生了改变,所以就得改变分配方3.案,并尽可能使得购买费用最小。

电力市场双边交易的阻塞管理

电力市场双边交易的阻塞管理

电力市场双边交易的阻塞管理
李嘉龙;姚常青
【期刊名称】《广东电力》
【年(卷),期】2006(19)11
【摘要】双边交易是电力市场中的一种重要交易形式,其相应的阻塞管理也具有重要的意义.鉴此,在潮流计算中应用直流法的基础上,提出了应用灵敏度的概念解决输电阻塞问题,将削减量与所占支路的容量的比例结合起来,做到公正、公平的阻塞管理.最后给出了一个5节点系统的算例,证明该方法的有效性.
【总页数】4页(P6-9)
【作者】李嘉龙;姚常青
【作者单位】广东省电力调度中心,广州,510600;天津大唐国际盘山发电有限责任公司设备部,天津,301900
【正文语种】中文
【中图分类】TM732;F407.61
【相关文献】
1.电力市场双边交易模式下应用边际网损系数法的网损分摊 [J], 徐耀良;夏磊;杨晓红;杨宁
2.基于可中断负荷拍卖模型电力双边交易阻塞管理 [J], 张森林;王钦;张尧;陈皓勇;荆朝霞;屈少青
3.基于可中断负荷拍卖模型电力双边交易阻塞管理 [J], 张森林;王钦;张尧;陈皓勇;荆朝霞;屈少青
4.基于可中断负荷拍卖模型的电力双边交易阻塞管理 [J], 张森林;王钦;张尧;陈皓勇;荆朝霞;屈少青
5.双边交易模式下考虑静态电压安全性的阻塞管理 [J], 王成山;魏炜
因版权原因,仅展示原文概要,查看原文内容请购买。

电力市场输电阻塞管理

电力市场输电阻塞管理
电力市场输电阻塞管理 摘要
本文拟合了各线路有功潮流关于各机组出力的近似表达式, 并设计阻塞费用 计算规则。 建立了不考虑潮流限值时建立下一时段出力分配模型和考虑潮流限值 时模型修改与最色费用计算。 首先根据功率的叠加性, 我们得出各线路上有功潮流与各发电机组出力有一 定的线性关系,因此建立多元线性回归模型拟合数据,运用 Matlab 求解得出带 系数的模型拟合度更高。 然后建立阻塞费用的计算规则。 除了需考虑电力市场规则外还考虑了序内容 量不能出力的部分和报价高于清算价的序外容量出力的部分, 我们同时考虑了发 电商实际叫价与清算价的差值和机组实际出力与预定出力的差值, 得到一个较公 平的补偿规则。 对于任务三, 我们在建立预定配电方案时先不考虑各个机组的爬坡速率, 只 对各机组的段容量按段价从大到小进行排序得到附录 1. 表 7,根据负荷预报值 得出分配方案, 再考虑爬坡速率的影响对比限值修改方案,最后得出各路线出力 值为: 机组 1 2 3 4 5 6 7 8 出力值 150 79 180 99.5 125 140 95 113.9 对于任务四, 我们先将上表数据带入模型一中得到初始分配预案,因考虑安全问 题,再与潮流限值进行比对得出是否发生输电阻塞。当发生阻塞时,我们根据建 立的阻塞费用计算模型提出目标函数, 将其转化为线性规划问题, 提出约束条件, 运用 LINGO 求解得出优化方案:
符号 T
符号说明 时段长(为 15 分钟) 第 i 个机组的出力值 调整后第 i 个机组的出力值 第 i 个机组的爬坡速率 第 j 条路线的有功潮流值 第 j 条路线的有功潮流限值 第 j 条路线的潮流安全裕度 分配预案中第 i 个机组的最终报价 调整方案后第 i 个机组的最终报价 清算价 对第 i 个序外容量的补偿 对第 i 个序内容量的补偿 负荷预报 MW MW MW/分钟 MW MW % 元/MWh 元/MWh 元/MWh 元 元 MW

电力市场联营-双边交易模式下的阻塞费用分摊

电力市场联营-双边交易模式下的阻塞费用分摊

电力市场联营-双边交易模式下的阻塞费用分摊
刘路登;张粒子;杨静
【期刊名称】《电网技术》
【年(卷),期】2006(30)15
【摘要】针对联营-双边交易模式电力市场的特点和我国处于电力市场改革初期的具体情况,提出了一种新的阻塞费用分摊方法,即首先按照约束分配系数将系统阻塞费用分摊到阻塞线路,然后利用潮流追踪方法得出各市场成员对阻塞线路的利用程度,最后将线路阻塞费用分摊给全部市场成员。

5节点系统算例结果表明,该方法能够正确有效地分摊市场成员的阻塞费用。

【总页数】4页(P50-53)
【关键词】电力市场;联营-双边交易模式;阻塞费用;约束分配系数;潮流追踪
【作者】刘路登;张粒子;杨静
【作者单位】华北电力大学电力市场研究所
【正文语种】中文
【中图分类】TM744
【相关文献】
1.联营与双边混合交易模式下结算盈余的分摊方法 [J], 帅勇;彭建春;龚演平
2.电力市场双边交易模式下的网损分摊方法 [J], 王秋梅;刘梓洪
3.联营-双边混合交易模式下的阻塞费用分摊 [J], 杨右虹;杨宗权;潘永贺;张义金
4.基于交流潮流的双边交易模式下阻塞费用分摊 [J], 苏健;彭建春;江辉;廖菁
5.双边交易模式下基于Aumann-Shapley值的阻塞费用分摊方法研究 [J], 杨洪明;段献忠
因版权原因,仅展示原文概要,查看原文内容请购买。

一种适用于我国电力市场的输电定价方法

一种适用于我国电力市场的输电定价方法

一种适用于我国电力市场的输电定价方法张瑞友;韩水;张近朱;汪定伟【期刊名称】《中国电机工程学报》【年(卷),期】2008(28)4【摘要】为使英国的输电定价模式更好地应用于我国电力市场,对其进行深入研究。

英国的电网使用费中不考虑机组报价,这主要适用于双边交易模式,而我国以集中交易模式为主。

英国的电网使用费中假设线路传输容量为无穷大,这适用于成熟期的电网,而我国电网扩建任务较重,经常有阻塞发生。

因此,在电网使用费的计算中引入了机组的报价和线路的传输容量,而且体现出阻塞线路和非阻塞线路的差别。

理论分析及最后的应用算例表明:改进后的方法保留了原来方法的主要优点,而且适用于我国当前集中交易为主的市场模式。

【总页数】5页(P78-82)【关键词】电力市场;输电定价;阻塞管理;集中交易;双边交易;电力交易新模式【作者】张瑞友;韩水;张近朱;汪定伟【作者单位】东北大学信息科学与工程学院,辽宁省沈阳市110004;国家电力监管委员会东北监管局,辽宁省沈阳市110006;东北电力科学研究院,辽宁省沈阳市110006【正文语种】中文【中图分类】TM731;F123【相关文献】1.资产一体化、市场均衡与资产定价--一种适用于企业资产定价的资本资产定价模型(CAPM)扩展形式 [J], 王永海;范明2.基于差别定价的电力市场输电服务定价研究 [J], 吴军;涂光瑜;罗毅;罗朝春3.发展中国家电力市场的定价之三——网嵌入的边际成本输电定价新方法 [J],4.在电力市场环境下网嵌入的边际成本输电定价新方法 [J], 言茂松;辛洁晴5.发展中国家电力市场的定价之三——网嵌入的输电定价方法 [J], 辛洁晴;言茂松因版权原因,仅展示原文概要,查看原文内容请购买。

输电阻塞管理的新方法述评(一):——基于潮流的可交易输电权

输电阻塞管理的新方法述评(一):——基于潮流的可交易输电权

输电阻塞管理的新方法述评(一):——基于潮流的可交易输电

方军;张永平;魏萍;祁达才;焦连伟;刁勤华;陈寿孙;倪以信;卢强
【期刊名称】《电网技术》
【年(卷),期】2001(25)7
【摘要】开放的电力市场需要输电网开放和有效的管理。

引入可交易输电权、利用市场机制解决网络阻塞问题是阻塞管理的重要趋势和研究热点。

该文分析了基于潮流的输电权的基本思想、特点、实质及其实施过程 ;三节点示例的分析表明输电权的引入可有效地实现市场价格 ,发现、缓解阻塞问题 ,同时具有规避价格风险的功能 ;对引入输电权所带来的问题进行了讨论。

【总页数】5页(P4-8)
【关键词】输电网;阻塞管理;潮流;输电权;电力系统;电力市场
【作者】方军;张永平;魏萍;祁达才;焦连伟;刁勤华;陈寿孙;倪以信;卢强
【作者单位】清华大学电机系;香港大学电机电子工程学系
【正文语种】中文
【中图分类】TM73;F407.61
【相关文献】
1.金融输电权和输电期权在输电阻塞管理中应用 [J], 郭金;谭忠富
2.基于分布式最优潮流算法的跨区输电阻塞管理研究 [J], 刘科研;盛万兴;李运华
3.输电阻塞管理的新方法述评(二):金融性输电权及与FGR之比较 [J], 张永平;方军;魏萍;祁达才;焦连伟;刁勤华;陈寿孙;倪以信;卢强
4.基于潮流影响系数的输电阻塞管理 [J], 李春燕;俞集辉;张谦
5.基于潮流输电权的阻塞管理分析 [J], 史江凌;俞集辉;张谦;彭光金
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

区域电力市场双边交易阻塞管理实用计算方法张永平1,焦连伟1,倪以信2,文福拴2,陈寿孙1,吴复立2(1.清华大学电机系,北京市100084;2.香港大学电机电子工程学系,香港)摘要:为了实现更大范围的资源优化配置,网络互联与区域间的功率交换日益频繁,并由此产生了跨大区的输电阻塞管理问题。

针对这一问题提出了跨大区电力市场双边交易的阻塞管理数学模型和实用计算方法。

在优化目标函数中,引入了合同的“权重因子”,以反映经济上的重要性或政策要求等因素,从而区别对待各合同,并使合同的调整量最小;在约束中,同时考虑了单条联络线的热稳定或动稳定约束和联络线割集的动稳定约束。

导出了组对组(点集对点集)合同内发电功率对线路功率的转移分布因子,从而将约束转化为发电机功率增量向量的线性函数,并通过阻塞优化校正算法得到最终的可行交易集合。

改进的IEEE 30节点算例证明了该算法的正确性和有效性。

关键词:电力市场;跨大区阻塞管理;双边合同;分布因子中图分类号:TM73;F123.9收稿日期:2002-12-06;修回日期:2003-02-24。

国家重点基础研究专项经费资助项目(G 1998020305);香港政府大学研究基金(RGC )资助项目。

0 引言作为保证系统安全的重要措施,输电阻塞管理从物理意义上讲并不是一个新的课题。

在传统的电力系统运行中,调度员通过“安全校正”,即直接调度发电资源或切除负荷来消除系统可能出现的过负荷。

但在市场环境下,由于信息和决策呈现出分散化趋势,同时电能交易方式呈现多样化,传统的安全校正措施难以完全适应新的运行环境。

因此,寻求适用于市场环境下多种交易模式的阻塞调度方法成为亟待解决的重要课题。

国内外学者对各种交易模式下的输电阻塞管理进行了研究[1~3]。

为了在更大范围内实现资源的优化配置,网络互联与区域间的功率交换日益频繁。

由于各区域之间联络线容量限制的存在,由此产生了区域间的阻塞管理问题。

文献[4~7]提出了基于分布式OPF 的跨大区阻塞管理方法,由于算法的复杂性和信息交换要求高等缺点,使得其在实际应用中存在一定限制。

文献[8]提出利用基于潮流的输电权(FGR )这一经济手段来解决阻塞问题,但当存在数量巨大的输电“瓶颈”(flow gate )时,这一方法的实施无疑会非常复杂。

文献[9]分析了北欧电力市场的“市场分裂”阻塞管理方法。

北欧电力系统是辐射状网络,因此这一方法是否适用于环网,需要进一步研究。

随着我国大区联网、西电东送等工程的实施,对于大区间输电阻塞问题的研究将具有重要的现实意义。

当存在大区间电能交易时,整个电力市场将呈现分层结构。

为了便于与区域电力市场及省级电力市场协调,一般区域间的电能交易将采用远期双边交易为主、现货交易为辅的模式。

针对这一模式,本文提出了跨大区电力市场双边交易的阻塞管理模型和实用计算方法。

在此模型中,某一区域的大发电机可以和其他区域的大负荷直接签订双边合同,也可以以区域电网公司为代表签订区域间的双边合同,通过本文提出的阻塞调度算法得到最终的可行交易集合。

改进的IEEE 30节点算例证明了本文算法的正确性和有效性。

1 阻塞管理数学模型1.1 多级电力市场的协调如前所述,当区域间存在电能交易时,将引发多级电力市场的协调问题。

针对分布式OPF 计算复杂以及信息交换量大的缺点,本文提出利用交易时序来解决市场间协调的方案。

其大体思路为:区域间的电能交易以远期双边交易为主,其时间跨度分别为年、月、星期;区域内也设立远期合约市场,但在时序上落后于区域间双边交易;区域间可设立现货市场(为进一步减少信息交换量,也可不设立),在时序上落后于区域内远期双边市场,但领先于区域内现货市场。

因此,市场交易时序依次为区域间远期双边、区域内远期双边、区域间现货、区域内现货。

时序上落后的交易均在已确认交易的基础上进行。

通过时序的先后进行多级电力市场的协调,能够有18第27卷 第14期2003年7月25日 电 力 系 统 自 动 力Automation of Electric Power Systems Vol .27 N o .14July 25,2003效地减少区域间的信息交换量,并且简单易行。

本文重点研究区域间的双边电能交易,对区域内的双边交易可采用同样的方法处理。

区域间双边电能交易在远期合同市场(forw ard contract market )进行(为简单起见,本文统称其为期货交易,这与通常的商品期货概念不同)。

本文的研究重点为区域间双边交易的安全校核和阻塞调度,该研究基于直流潮流假定,以充分利用其线性特征及可叠加性,必要时可进一步进行交流潮流校验,以及基于交流灵敏度的小量调整。

1.2 阻塞管理总体计算流程对于时间跨度为年、月、星期的双边期货交易,其可行性校核均采用相同的算法。

所不同的是,下一时间跨度交易可使用的输电容量为上一时间跨度交易校核结束后的剩余容量。

对于每一个时间跨度,相应的阻塞分析和调度计算流程如图1所示,其核心部分为阻塞调度算法模块,以消除线路潮流或联络线割集潮流的越界。

图1 双边合同校核流程Fig .1 Flo w chart of bilateral contract feasibility check1.3 阻塞调度优化问题的数学模型1.3.1 合同的表示本文考虑的区域间双边期货合同有2种:点对点合同及组对组(点集对点集)合同。

点对点合同规定了某区域的某台发电机i 向另一区域某负荷j 的供电功率。

点对点合同可表示为:τ=[0 … P ij … -P ij … 0]T (i ) (j )式中:交易τ为N F 维向量;N F 为参加区域间期货交易的发电机和负荷总数;P ij 表示从发电机i 到负荷j 的功率流动。

组对组(点集对点集)合同一般为各区域系统调度员(SO )代表各自区域内发电机或负荷与其他区域发电机或负荷签订的双边合同。

它规定了从某区域的若干台发电机向另一区域若干个负荷的供电功率。

组对组合同可表示为:τ=[0 … P m … -P n … 0]T式中:非零元素表示相应节点有功率注入(正值)或析出(负值)。

非零元素代数和应为0(忽略网损)。

点对点合同可以看做是组对组合同的特殊形式。

区域间的所有期货交易(均为双边)可以表示为交易矩阵T :T =[τ1 τ2 … τt ](1)式中:t 为发电机和负荷之间的总合同数。

参与区域间期货交易的发电机和负荷节点的总注入功率矢量可表示为:P =TU (2)式中:U =[1 … 1]T,为t 维列向量。

此总功率注入矢量用于潮流计算。

1.3.2 优化目标函数考虑到不同的双边交易重要性和经济性不同,本文引入权重因子作为交易重要性和经济性的指标,权重因子由SO 根据市场规则给定。

当权重因子为1时,各交易平等对待;当某交易有优先级时(如愿意支付更高的输电费用,或根据政策需要优先考虑),权重因子取较大的值。

可行性校验的目标函数为计及权重系数后的交易调整量最小:min∑i ∈{G P}w i ΔP 2i+∑Kk =1w k∑i ∈{G k}ΔP i ,k2(3)式中:w i 为点对点交易i 的支付因子;w k 为组对组交易k 的支付因子;ΔP i 为点对点交易调整量(发电节点i 的调整量);ΔP i ,k 为对于第k 个组对组交易,发电节点i 的功率调整量;{G P }为所有点对点交易的发电机集合;{G k }为第k 个组对组交易的发电机集合;K 为所有组对组交易的个数。

1.3.3 约束条件优化目标函数相应的约束条件如下:a .点对点交易平衡约束ΔP i =ΔD j i ∈{G P },j ∈{D P }(4)式中:ΔD j 为点对点交易调整量(负荷节点j 的调整量);{D P }为所有点对点交易的负荷集合;节点i 和j 属于同一个交易。

b .组对组交易平衡约束∑i ∈{G k}ΔP i ,k=∑j ∈{D k}ΔD j ,ki ∈{G k },j ∈{D k }, k ∈K(5)19·电力市场专栏· 张永平等 区域电力市场双边交易阻塞管理实用计算方法式中:ΔD j ,k 为第k 个组对组交易负荷节点j 的调整量;{D k }为第k 个组对组交易的负荷集合。

c .发电机出力约束P i ,m in -P G i ,0≤ΔP G i ≤P i ,max -P G i ,0 i ∈{G }(6)式中:P i ,m ax 和P i ,min 分别为发电机i 出力上、下限;P G i ,0为发电机i 初始出力;ΔP G i 为发电机i 出力调整量;{G }为所有发电机集合。

d .负荷变化约束D j ,mi n -D j ,0≤ΔD j ≤D j ,max -D j ,0 j ∈{D }(7)式中:D j ,max 和D j ,m in 分别为负荷j 调整范围的上、下限;D j ,0为负荷j 初始量;ΔD j 为负荷j 调整量;{D }为所有负荷集合。

e .线路潮流约束-P l ,ma x -P l ,0≤ΔP l ≤P l ,ma x -P l ,0 l ∈{L }(8)式中:P l ,max 为支路l 的潮流限制,可以取热稳定和动稳定约束中值小者;P l ,0为支路l 的初始潮流;ΔP l 为支路l 的潮流变化量;{L }为所有线路集合。

f .联络线割集潮流约束-P T ,max -∑l ∈TP l ,0≤∑l ∈TΔP l ≤P T ,m ax -∑l ∈TP l ,0T ∈{T c ut }(9)式中:P T ,m ax 为割集T 的潮流限制,可取为动稳定约束;∑l ∈TP l ,0为割集T 的初始潮流和;∑l ∈TΔP l 为割集T的潮流变化量;{T c ut }为所有联络线割集集合。

1.3.4 组对组合同中负荷的调整对于组对组合同,当发电机进行调整时,必须按照某种原则对合同中的负荷进行调整。

本文按照等比例原则进行调整,即:D j ,k =D j ,k ,0∑i ∈{G k}P i ,k∑i ∈{G k}P i ,k ,0 j ∈{D k }, k ∈K(10)式中:P i ,k ,0和D j ,k ,0分别为对于第k 个组对组交易,发电节点i 的初始功率和负荷节点j 的初始功率。

2 约束的计算及处理2.1 分布因子的定义及计算本文利用直流潮流假定和分布因子的概念,充分利用了其线性特征和可叠加性以加快计算速度。

a .点对点交易的分布因子对于点对点交易,定义功率传输分布因子(PTDF )为某节点对(m ,n )之间的电能传输发生单位改变时,某支路i -j 的潮流变化量,即:F PTDF ,ij ,mn =ΔP ij ΔP mn =X im -X jm -X in +X jnx i j(11)式中:ΔP mn 为节点m 和节点n 之间的电能交易变化量;ΔP ij 为支路i -j 的潮流变化量;x i j 为支路i -j 的电抗;对于直流潮流方程B θ=P ,令X =B -1,X im 为X 矩阵第i 行第m 列元素。

相关文档
最新文档