电力场的输电阻塞管理

1电力市场的输电阻塞管理摘要：随着电力系统改革的进行和用电紧张的缓解，电力市场化将进入新一轮的发展。本文根据电力市场交易规则和阻塞管理原则，运用统计学、数学规划等工具建立了一个电力市场输电阻塞管理的数学模型，同时对于给定的题设条件，对模型进行了求解和讨论。首先，运用逐步回归的方法拟合了各线路潮流值关于机组出力的实验数据，得出了6个线路潮流值的经验回归公式。其次，主要解决本文的两个核心问题：阻塞费用的计算；阻塞时如何调整出力以满足最大的安全和负荷需求。为了便于数学上处理和计算，我们采用的阻塞费用()iU计算公式是最简单的线性形式，然后以81()iiU最小为目标函数建立一个规划模型，然后分别根据问题3——问题5的条件，分别计算求解。结果发现，在负荷需求为982.4MW时，基本上可以通过阻塞管理可以较安全地满足需求，出力分配为1150x，279x，3180x，499.5x，5125x，6140x，795x，8113.9x，但是负荷需求为1052.8MW时是无论如何都不可能满足的，所以必须拉闸限电。此时是可以看成一个双目标规划，即要安全性尽量高，出力和又要尽可能大。求解时，采用列举不同安全裕度的形式，得到一个相对较优的解：1153x，288x，3228x，499.5x，598x，6100.1x，7102.1x，8117x一条经验性的规律是：线路潮流上限是模型的最主要的约束，是电力运营的瓶颈。本文通过一定的合理的简化和假设，建立了一个较为简单的优化（规划）模型，并借助Matlab程序提供了简单的求解方法，最后给出了建议和评价。关键词：线形回归数学规划阻塞费用2一问题重述否能否是是否电网公司下一时段的负荷需求预报分配预案电力市场交易规则计算各线路潮流值，判断是否阻塞清算，完成一个时段交易在保证线路安全下再分配各机组出力判断能否满足负荷需求发电厂商判断能否满足负荷需求限电3上图给出了本问题的一个处理流程：电网公司根据各机组当前出力情况以及下一时段的负荷需求预报，发电厂商则根据市场交易原则，得出下一时段的各机组的出力分配预案，网方以此计算出各线路的有功潮流，判断是否会出现输电阻塞的情况。如果不出现，接受各机组的出力分配预案；否则，根据阻塞管理原则进行调整。根据阻塞管理原则，当改变出力方案时，就会出现序内容量不能出力的部分以及报价高于清算价的序外容量的部分，使得发电方产生损失，因此网方应该给予一定的补偿，这部分就是阻塞费用。现在考虑的电网有6条主要线路，厂方有8台发电机组，当前各机组的出力方案，各线路上的有功潮流以及围绕当前方案的一些实验数据给定。要求：设计一种简明、合理的阻塞费用计算规则，可以均衡的考虑序内容量不能出力的部分以及报价高于清算价的部分，并能尽可能的减少阻塞费用。以此为基础，加上有功潮流的近似式，当线路发生阻塞的时候，相应的对出力方案作出一些调整，消除输电阻塞。二基本假设以及相关符号说明1．出力连续性假设各机组的出力是连续的，受到机组爬坡速度的限制，不存在出力的突变。2．户负荷的需求假设对每时段的的负荷预报仅仅是该时段平均负荷的预报，在每个考察时段内实际负荷会在预报负荷的周围波动，但总体上不影响用户的用电，也不存在安全问题。3．户需求的满足假设满足用户的负荷需求只是在该时段的开始时刻达到预报负荷即可，而在该时段内出力会不断变化，以达到下一时段的用户负荷预报需求。4．清算价格的假设假设方案改变后，机组的同一时段内的清算价固定，不发生变化5．其他假设为了更好地考虑阻塞费用，考虑发电商(厂方)和电网公司(网方)两个实体，网方从该厂方购买电能，但该发电商可能有多个买家(电网公司)。符号说明：0ix1,2,3....8i第i台机组的当前出力值（即方案0）ix1,2,3....8i第i台机组的下一时段的出力方案0iiixxx1,2,3....8i第i台机组的出力值差额iv1,2,3....8i第i台机组的爬坡速率iy1,2,3....6i第i条线路的有功潮流值4Ai1,2,3....6i第i条线路的限值。iz1,2,3....6i第i条线路的限值的裕度p表示下一阶段的负荷预报三具体建模过程【1】由于给定了围绕当前各机组的出力以及各线路上的有功潮流，我们利用这32组数据，采取线性回归的的方法，对数据进行拟合。以各线路上的有功潮流为被解释变量，用各机组的出力情况作为解释变量，先用SPSS统计软件作多元函数的一次线性拟合，即0112233445566778857123468(,,,,,,,)xxxxxxxxyfxxxxxxxx得到各线路上有功潮流关于各发电机组出力的函数关系式，通过对该拟合作F检验，以及对拟合优度的比较，发现在0.05以及0.01的置信度上，拟合都是高度显著。由于采用的是逐步回归（Stepwise）方法，随着拟合的进行，在每一次取舍之后，拟合优度都是逐步增大的，因此可以认为该结果还是比较准确的反映了对实际数据的拟合情况。从另一方面考虑，由于各机组是独立发电，互不影响，而每条线路的有功潮流均机会均等的与8台机组的出力情况有关，因此，该回归模型中应该不存在交互效应，至于较高次数的回归，我们认为也没有太大的实际意义，因为从上述回归来看拟合已经是高度显著，在将次数提高结果也不会有太大的的改进，况且从数据的变化来看，1—4方案第一机组的出力变化，而其他机组出力不变，对应于1—4方案的各线路的有功潮流第五条线路的负荷几乎不变，而其他线路均有不同程度的变化，得到的关系式中也是如此，第五线路的潮流值与第一机组的出力无关，其他关系式也满足；而从实际情况来考虑，用一次多元函数来模拟也是比较合理，符合一定的实际情况，因此我们用一次多元函数来解释各线路上的有功潮流与各机组的出力情况的关系，以上拟合得经验回归函数如下：57123468(,,,,,,,)iiyfxxxxxxxx，即：51123476110.01640.08310.04880.05320.12000.02510.12240.1211yxxxxxxx52124678131.21890.05460.12790.03330.08690.11240.01890.0987yxxxxxxx5312348108.39220.07010.06040.15710.01010.12380.2021yxxxxxx554123476877.48170.03450.10240.20520.02080.01180.00600.14490.0766yxxxxxxxx55234768133.12490.24300.06470.04120.06540.07010.00400.0091yxxxxxxx5612347120.79300.23770.06040.07800.09290.04680.1662yxxxxxx由上述函数关系表达式可以看出，各线路上的有功潮流值基本在一个固定值的上下波动，由于各机组出力的权数很小，均小于1，所以该波动值一般不会很大，在考虑每条线路是否会发生阻塞时，可以先用（110，131，108，77，133，120）即函数表达式的常数项，与各线路的潮流限值相比较，便可以对各线路的阻塞情况有个大概的了解，以便以后的方案调整可以简化一点，其中如果潮流值为负值仅仅表示方向的不同而已，不涉及正负的关系。在上述关系式中，我们仅选了围绕0方案的32组数据作回归，把0方案的数据当作当前的运转情况，在以后的计算中会进一步用到该组数据。【2】这里我们给出一种最简单的阻塞费用的计算规则：1()||*4iiiUxxC这里C表示清算价格()iiUx表示对第i台机组的补偿。当0ix时，表示有序外容量，厂方以低于对应报价的清算价出力，厂方损失，故网方应对厂方补偿。当0ix时，表示有序内容量，厂方丧失了向其他网方输出电能赢利的机会(即“机会成本”)，因而网方同样要对厂方补偿。由于改变出力方案后，各机组的出力或升高或降低，即ix有正有负，所以这里我们取了绝对值，表明改变预案责任在于网方，所以无论如何都是网方赔偿。这里取1()||*4iiiUxxC，是因为它比较简单直观，易于计算。这样，我们得到如下的规划模型：【3】根据市场交易原则，这一阶段要得出各机组的预分配方案。由市场交易规则可知市场交易中心根据机组当前报价、当前出力以及出力改变速率，按段价从低到高选取各机组的段容量或其部分，直到其和等于等于预报负荷。我们将各机组的段价由低到高进行排序便可以得到相应的段容量的入选部分，直到段容量之和等于预报负荷6具体的算法为：设ijs为第i机组第j段的段容量，ijr第i机组第j段的段价1,2,3....8i,1,2,3....8j;将ijr按照由低到高排序，对应的每个ijr有自己的段容量，将段容量累加直到等于负荷需求为止，此时负荷P=11mnijijs（其中ijs可能取某容量的一部分，m,n为入选的段个数）；1niijjxs（1,2,3....8j）清算价为最后一个入选的ijs的ijr；P=982.4时，由以上的排序可以得到各机组各阶段的容量值，所以对各机组第一到第六段容量入选部分累加就可以得到各机组的出力分配预案，结果如下：表（一）机组/段1234567ix17005000300150230020815607931100400300018045551010101001005755150151501256950102001501407501551510009587002002003.4113.4此时的段价为g=302元/MWh由于有机组爬坡速度的限制，在一个交易时段内要考虑在当前的出力情况下机组能否攀升到预案所分配给各机组的出力，因此我们要作以下的判断：001515iiiiixvxxv（*）将各机组当前的出力及各机组的爬坡速率代入上述判断式中可以得到：187153x725888x3132288x460.599.5x598152x695155x769102.1x863117x可以看出第四机组的出力超出了其所能承受的最大值，因此还要对方案进行改进。我们将第四机组中超出的部分去掉，差额再按照市场交易规则，对剩下的段容量再重复以上的过程，因此可以得到重新分配后的出力情况，如下：表（二）机组/段1234567ix17005000300150230020815607931100400300018045551010109.5099.557551501515012569501020015014075015515103.409587002002003.9113.9以上结果即是分配预案。此时的清算价为303元/MWh【4】根据第一部分中得到的线路潮流值关于各机组出力的关系式以及第三部分中的预分配方案，可以得到各线路的潮流值，具体如下：表（三）1y2y3y4y5y6y173.35141.01-150.96120.912136.82168.51由以上数值与线路潮流限值以及起安全裕度相比较可以看出，虽然6条线路都在安全裕度以内，但第1、5、6条线路已经超过限值，已经产生了线路阻塞，所以要对原方案进行调整，调整的目标是对原预案作出适当的调整使各线路的潮流值尽量不要超过限值，即不产生阻塞，并使网方赔偿的阻塞费用最小，还应满足原负荷需求。我们只需以阻塞费用为目标函数作最优化的规划，由第二8部分的的阻塞函数可以得到以下规划：81810()s.t.982.4|-|*151,2,...,8y1,2,...,6iiiiiiiiiiMinUxxxxviAi由于88111()*||4iiiiiUxCx含又绝对值难以分析处理，我们将目标函数替换成821*()iiCx，这里我们认为821*()iiCx的最优解同时也是81*||iiCx的最优解。利用Matlab程序fn1.m(见附录)