电力市场输电阻塞管理模型--数学建模全国赛优秀论文1内部有源程序

电力市场输电阻塞管理模型(轩辕杨杰整理)摘要本文通过设计合理的阻塞费用计算规则，建立了电力市场的输电阻塞管理模型。首先通过对各机组出力方案实验数据的分析，用最小二乘法进行拟合，得到了各线路上有功潮流关于各发电机组出力的近似表达式。按照电力市场规则，确定各机组的出力分配预案。如果执行该预案会发生输电阻塞，则调整方案，并对引起的部分序内容量和序外容量的收益损失，设计了阻塞费用计算规则。我们引入危险因子反映输电线路的安全性，根据安全且经济的原则，把输电阻塞管理问题归结为：以求解阻塞费用和危险因子最小值为目标的双目标规划问题。采用“两步走”的策略，把双目标规划转化为两次单目标规划：首先以危险因子为目标函数，得到其最小值；然后以其最小值为约束，找出使阻塞管理费用最小的机组出力分配方案。当预报负荷为982.4MW时，分配预案的清算价为303元/MWh，购电成本为74416.8元，此时发生输电阻塞，经过调整后可以消除，阻塞费用为3264元。当预报负荷为1052.8MW时，分配预案的清算价为356元/MWh，购电成本为93699.2元，此时发生输电阻塞，经过调整后可以使用线路的安全裕度输电，阻塞费用为1437.5元。进一步，我们分析了各线路的潮流限值调整对最大负荷的影响，并据此给电网公司提出了建议。最后，我们对模型进行了评价，并提出了模型的改进方案。1一、问题的重述我国电力系统的市场化改革正在积极、稳步地进行，随着用电紧张的缓解，电力市场化将进入新一轮的发展，这给有关产业和研究部门带来了可预期的机遇和挑战。电网公司在组织电力的交易、调度和配送时，必须遵循电网“安全第一”的原则，同时按照购电费用最小的经济目标，制订如下电力市场交易规则：1、以15分钟为一个时段组织交易，每台机组在当前时段开始时刻前给出下一个时段的报价。各机组将可用出力由低到高分成至多10段报价，每个段的长度称为段容量，每个段容量报一个段价，段价按段序数单调不减。2、在当前时段内，市场交易-调度中心根据下一个时段的负荷预报、每台机组的报价、当前出力和出力改变速率，按段价从低到高选取各机组的段容量或其部分，直到它们之和等于预报的负荷，这时每个机组被选入的段容量或其部分之和形成该时段该机组的出力分配预案。最后一个被选入的段价称为该时段的清算价，该时段全部机组的所有出力均按清算价结算。电网上的每条线路上有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值，称为输电阻塞。当发生输电阻塞时，需要按照以下原则进行调整：1、调整各机组出力分配方案使得输电阻塞消除；2、如果1做不到，可以使用线路的安全裕度输电，以避免拉闸限电，但要使每条线路上潮流的绝对值超过限值的百分比尽量小；3、如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度，则必须在用电侧拉闸限电。调整分配预案后，一些通过竞价取得发电权的发电容量不能出力；而一些在竞价中未取得发电权的发电容量要在低于对应报价的清算价上出力。因此，发电商和网方将产生经济利益冲突。网方应该为因输电阻塞而不能执行初始交易结果付出代价，网方在结算时应该适当地给发电商以经济补偿，由此引起的费用称之为阻塞费用。网方在电网安全运行的保证下应当同时考虑尽量减少阻塞费用。现在需要完成的工作如下：1、某电网有8台发电机组，6条主要线路，附件1中表1和表2的方案0给出了各机组的当前出力和各线路上对应的有功潮流值，方案1~32给出了围绕方案0的一些实验数据，试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。2、设计一种简明、合理的阻塞费用计算规则，除考虑电力市场规则外，还需注意：在输电阻塞发生时公平地对待序内容量不能出力的部分和报价高于清算价的序外容量出力的部分。3、假设下一个时段预报的负荷需求是982.4MW，附件1中的表3、表4和表5分别给出了各机组的段容量、段价和爬坡速率的数据，试按照电力市场规则给出下一个时段各机组的出力分配预案。4、按照表6给出的潮流限值，检查得到的出力分配预案是否会引起输电阻塞，并在发生输电阻塞时，根据安全且经济的原则，调整各机组出力分配方案，并给出与该方案相应的阻塞费用。5、假设下一个时段预报的负荷需求是1052.8MW，重复3~4的工作。2二、问题的分析市场交易-调度中心在一个时段内的工作流程如图1所示。首先根据电力市场交易规则及负荷预报需求确定下一时段各机组出力的分配预案，再通过计算各线路潮流值判断是否会出现输电阻塞。若出现，则按输电阻塞管理原则对预案进行调整。图1市场交易-调度中心工作流程图是当前时段开始计算各机组当前出力值确定下一时段分配方案计算各线路有功潮流出现输电阻塞是否可消除是否可用安全裕度输电计算阻塞费用确定分配方案当前时段结束进入下一时段拉闸限电否是否是否3根据功率的叠加原理，各线路上有功潮流应为各发电机组出力的线性组合，考虑对所有实验数据采用最小二乘法进行线性拟合，从而得到各线路有功潮流关于各发电机组出力的近似表达式。得到分配预案后，代入近似表达式便可计算各线路上的潮流值。为保证电网的安全，每条线路潮流的绝对值超过潮流限值的百分比应尽量小。若使各线路中潮流超出的百分比中最大的值尽量小，就可保证所有线路上潮流超出的百分比较小，即电网相对较为安全。在电网安全运行的保证下应当同时考虑尽量减少阻塞费用。阻塞费用分为两个部分：一是对序内容量不能出力部分的补偿；二是对报价高于清算价的序外容量出力部分的补偿。以每个机组各自的报价作为其边际成本，则该机组单位出力的绝对盈利为清算价与报价的差值，因此，补偿的主要目的是解决由于方案调整导致的获利变化的问题。该阻塞管理问题归结为在一定约束条件下的最优化问题。优化目标为使潮流超出现值的百分比尽量小，同时尽可能减少阻塞费用。三、基本假设1、机组当前出力是对机组在当前时段结束时刻实际出力的预测值；2、每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻；3、机组在单位时间内能增加或减少的出力相同，出力值为爬坡速率；4、各个发电机组出力相互独立，即出力不受其他机组影响。四、定义符号说明1、名词解释电力市场：电力的买方和卖方相互作用以决定其电价和电量的过程（见[1]第4页）；边际成本：在一定的时期内，增加一个单位产量所需支付的成本；序内容量：在电力市场中通过竞价取得发电权的发电容量；序外容量：在竞价中未取得发电权的发电容量；爬坡速率：机组在单位时间内能增加或减少的出力值；最终报价：进行结算时，机组分配到的出力对应的报价。2、符号说明ix：第i个机组的出力值；单位：兆瓦，记作MWix：调整后第i个机组的出力值；单位：MWiv：第i个机组的爬坡速率；单位：MW/分钟jl：第j条线路的有功潮流值；单位：MWjA：第j条线路的初始潮流值；单位：MWjL：第j条线路的潮流限值；单位：MWja：第j条线路的潮流的安全裕度；ip：分配预案中第i个机组的最终报价；单位：元/MWhip：调整方案后第i个机组的最终报价；单位：元/MWhif：对第i个序外容量的补偿；单位：元if：对第i个序内容量的补偿；单位：元4X：负荷预报；单位：MWP：清算价；单位：元/MWhT：时段长，T为15分钟；五、模型的建立1、建模前的准备1）有功潮流近似表达式的确定每条线路上的有功潮流取决于电网结构和各发电机组的出力，问题所研究的电网有8台发电机组，6条主要线路，附件1中的表1和表2的方案0给出了各机组的当前出力和各线路上对应的潮流值，方案1~32给出了围绕方案0的一些实验数据。根据功率的叠加原理，我们认为各线路上有功潮流应为各发电机组出力的线性组合，随机抽取几组方案进行检验，得到以下图形：125130135140145150155160165.5166166.5167167.516874767880828486889092141141.5142142.5143143.5线路1受机组1的影响线路2受机组2的影响65707580859095100-144.35-144.3-144.25-144.2-144.15-144.1-144.05-144606264666870727476135135.1135.2135.3135.4135.5135.6135.7135.8135.9136线路3受机组4的影响线路5受机组7的影响图2对实验方案的分析从图形中我们发现，有功潮流受到各机组的影响近似成线性关系，因此假设有功潮流关于各个机组出力的函数关系式为jiijijAxkl81其中jl表示第j条线路上的潮流值，jik表示第j条线路受第i台机组影响的比例系数，ix表示第i台机组的出力，jA表示第j条线路对应的初始潮流值。5对应每一条线路，根据表1表2中的数据可列出关于未知数jik（i=1,2,…,8）的32个方程的超定方程组，在Matlab下编程求解方程组（源程序见附件2），得到结果如下：0466.00655.00120.01247.00867.00257.00929.00412.00209.00099.00332.01199.00781.00647.02050.01565.00001.00528.00607.02428.01028.00620.01275.00478.02376.00003.00346.00694.00547.00826.054321654321xxxxxllllll+8481.1201334.1336116.779928.1083521.1314775.1100004.00092.00763.02012.00985.00015.01664.00039.01452.00028.00186.0122.00003.00700.00057.00024.01127.01216.0876xxx2）阻塞费用计算规则的设计当改变根据电力市场交易规则得到的各机组出力分配预案时，一些通过竞价取得发电权的发电容量（称序内容量）不能出力；而一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的清算价上出力。以机组的最终报价作为其边际成本，则该机组单位出力的绝对盈利为清算价与报价的差值，因此，补偿的主要目的是解决由于方案调整导致获利变化的问题。我们设计的阻塞费用计算规则如下：①对于序内容量：由于方案的调整，使得一些机组的出力值减少，减少部分的获利值消失。为解决这部分冲突，网方赔偿该机组应得的获利值，有调整量调整前报价）（清算价补偿费用即)()(iiiixxpPf②对于序外容量：方案调整后，一些机组由于出力增加，其边际成本（报价）也随之增加，但由于清算价保持不变，机组不得不在低于其报价的清算价上出力，导致了获利损失。因此，网方对调整的出力部分造成的损失应给予补偿，有调整量清算价）（调整后报价补偿费用即)()(iiiixxPpf总的阻塞费用即为81)(iiifff2、约束条件的讨论1）爬坡速率由假设1，在当前时段，市场交易-调度中心预测出各机组结束时刻的实际出力，即当前出力值，由于爬坡速率的约束，当前出力在时段长15分钟内改变的值有限，有6Tvxxiii0min，Tvxxiii0max其中，对于第i个机组，0ix为当前的出力值，minix为其下一时段出力值的下限，maxix为其下一时段出力值的上限，iv为爬坡速率，T为时段长。因此，下一时段的出力值须满足：],[maxminiiixxx2）线路