基于功率灵敏度和经济补偿最小化的线路过负荷紧急控制方法
2017-05-24樊世通
徐 岩 ,郅 静 ,樊世通
(1.华北电力大学 新能源电力系统国家重点实验室,河北 保定 071003;2.国网江苏省电力公司检修分公司苏州运维分部,江苏 苏州 215000;3.国网江苏省电力公司苏州供电公司,江苏 苏州 215000)
0 引言
电力系统中,线路因发生故障被切除时会引发潮流转移,这很容易导致其他线路过载,从而引起过负荷后备保护动作切除线路,易导致连锁跳闸,甚至大停电事故[1]。针对因潮流转移导致过载的线路,快速闭锁其后备保护,并采取调整发电机出力甚至切负荷等紧急控制措施消除线路过载,可有效防止线路连锁跳闸的发生。
现有紧急控制方法有灵敏度法[2-4]和优化规划法[5-6]。灵敏度法根据电网各节点对过载线路的控制灵敏度逐个选择最佳控制节点,每一轮控制都选择当前控制灵敏度最大的控制节点组,该方法容易实现调整量最小的目标,但无法考虑紧急控制中的经济因素。优化规划法根据各节点对过载线路的控制情况及电网各线路的功率限制约束构造数学模型,利用最优化方法求解该模型得到控制方案,有较好的全局性和经济性,但有时计算量较大,同时,在电力市场环境下,独立系统调度员ISO(Independent System Operator)在采取紧急控制措施时,需要对被迫调整出力的发电机和被迫切负荷的用户提供一定的经济补偿[7],因此,制定能够消除线路过载且经济补偿最小的紧急控制方案具有实用价值。综上所述,本文采用优化规划法,并利用相应方法减少计算量。
已有的紧急控制方法[2-5,8-10]在计算节点对过载线路的控制灵敏度时,大都采用直流潮流法中的发电转移分布因子GSDF(Generation Shift Distribution Factor)。GSDF的计算依赖于平衡机的选择,令平衡机对所有线路的灵敏度为0[3],不使平衡机参与紧急控制,易导致无法得到最优控制方案,同时,当平衡机不同时,得到的紧急控制方案也可能不同[11]。文献[12]提出了功率灵敏度的概念,计算出包括平衡机在内的电网各节点对过载线路的控制情况,证明功率灵敏度在紧急控制中应用效果好,同时,在计算节点与线路间的功率灵敏度时,将局部灵敏度计算法和稀疏向量技术结合,大幅减少了计算量和内存量。
本文提出一种基于功率灵敏度和经济补偿最小化的线路过负荷紧急控制方法。利用功率灵敏度矩阵反映各节点注入功率变量对线路功率变量的影响,以经济补偿最小化为目标,利用非线性优化得到满足线路过负荷消除约束、节点可调量约束、系统功率平衡约束及正常线路冗余量约束的紧急控制方案。该方法可计算平衡机对电网各线路功率的灵敏度,令平衡机参与到紧急控制中,提高了电网保持稳定运行的能力;同时,以经济补偿最小为目标,使得到的紧急控制方案更具有实用价值。仿真算例验证了本文方法的正确性和优越性。
1 功率灵敏度
文献[12]中详细介绍了功率灵敏度矩阵的推导过程。当电网发生断线故障时,根据广域测量系统提供的实测信息刷新电网拓扑结构和节点电压相量数据,可以快速计算此时的功率灵敏度矩阵,得到各节点注入功率变量对支路功率变量的控制情况。
电网中支路Bk的功率变量与节点Ni的注入功率变量之间的功率灵敏度βk-i为:
其中,λk-i为支路Bk的电流相量与节点Ni的注入电流相量之间的电流相关度系数;Uk,B和 φk,B分别为支路 Bk的首端电压模值和相角;Ui,N和 φi,N分别为节点Ni的电压模值和相角。
功率灵敏度矩阵结合网络拓扑参数和实际电网运行情况衡量节点注入功率对线路潮流的控制情况,其计算过程不依赖于平衡机的选择,能准确计算平衡机对各线路功率的灵敏度,令平衡机参与到紧急控制中。
设电网中发电机节点数为m,负荷节点数为f,当采取调节发电机出力或者切负荷等紧急控制措施时,支路 Bk的功率变量 ΔPk,B如式(2)所示。
其中,ΔPi,G为发电机节点 NGi的发电调整量;ΔPj,L为负荷节点NLj的切负荷量。
2 最优化求解紧急控制方案
2.1 最优化目标
为保证供电可靠性,紧急控制中应尽量依靠调整发电机出力消除过载,避免切负荷。当仅通过调整发电机出力无法得到满足系统安全运行的控制方案时,需要进一步考虑采取相应的切负荷措施[8]。同时,在消除线路过载的紧急控制中,ISO需要以一定的补偿价格对被迫调整出力的发电机和被迫切负荷的用户进行经济补偿[7]。因此,在尽量避免切负荷的情况下,使经济补偿最小化的最优化目标如式(3)所示。
其中,ρi为发电机节点NGi的单位功率变化补偿价格;λj为负荷节点NLj的切除单位负荷的补偿价格;为尽量避免切负荷,取M为一个很大的正数,保证在调整发电机出力无法消除线路过载时才采取切负荷措施,本文取M=103。
式(3)所示的最优化目标能够有效保证所制定的紧急控制方案在尽量避免切负荷的情况下,使经济补偿最小化。
2.2 最优化约束
2.2.1 线路过负荷消除约束
对于过载线路集合A中任一条过载线路Ba,紧急控制方案都可消除其过载,即满足:
其中,Pa,B为紧急控制前线路 Ba的传输功率,Pa,Bmax为线路Ba的最大传输功率;考虑到功率灵敏度及调整量计算误差等因素,定义可靠系数η<1,本文取η=0.9,实际电网中根据情况适当取值。
2.2.2 节点可调量及系统功率平衡约束
(1)节点分类。
根据各节点对所有过载线路的综合控制作用,定义节点Ns的综合灵敏度如下:
其中,fA-s为节点Ns的综合灵敏度;ξa为线路Ba的过载率,线路过载程度越大,其在综合灵敏度计算公式中的重要性权重越大。
根据节点的综合灵敏度对节点进行分类:
a.若fA-s>0且节点Ns为可调整出力的发电机节点,说明减少该节点的发电机出力可以降低所有过载线路的整体过载程度,因此,将该节点归入减出力节点集合;
b.若fA-s<0且节点Ns为可调整出力的发电机或可切负荷,说明增加该节点的发电机出力或减少该节点的负荷可以降低所有过载线路的整体过载程度,因此,将该节点归入加出力节点集合。
(2)节点可调量及系统功率平衡约束。
减出力节点集合中发电机节点NGi的出力调整量 ΔPi,G需满足:
其中,Pi,G为发电机节点 NGi的当前发电量;Pi,Gmin为发电机节点NGi的最小发电量。
加出力节点集合中,ΔPi,G需满足:
其中,Pi,Gmax为发电机节点NGi的最大发电量。
当平衡机为加出力节点时,考虑平衡机在电网运行中需要承担的调频及平衡电网微小功率差额的任务[13-14],计算其可增出力时要预留出一定裕度,即对于平衡机Gp,需要满足的约束是:
其中,Py为平衡机预留裕度,本文中取Py=10 MW,实际电网运行中可根据情况适当取值。
加出力节点集合中负荷节点NLj的切负荷量ΔPj,L不大于该节点的总负荷量 Pj,L,即满足:
同时,为保证系统功率平衡,所有发电机节点的出力调整量和负荷节点的切负荷量需满足:
2.2.3 正常线路冗余量约束
紧急控制中要保证系统不出现新的线路过载,对于正常线路Bl,需满足如下约束:
其中,Pl,B为紧急控制前线路 Bl的传输功率;Pl,Bmax为线路Bl的最大传输功率。
本文将支路热稳定功率极限作为其最大传输功率,实际电网运行中依情况而定[15]。
利用功率灵敏度矩阵反映各节点注入功率变量对线路功率变量的影响,以经济补偿最小化为目标,利用非线性优化可得到满足线路过负荷消除约束、节点可调量约束、系统功率平衡约束及正常线路冗余量约束的紧急控制方案。
2.3 计算量分析
制定紧急控制方案时只需计算发电机节点和负荷节点对电网线路的功率灵敏度,利用局部灵敏度计算法[12]和稀疏向量技术,可以大幅减少计算量和内存量。同时,正常线路冗余量约束保证了正常线路在紧急控制过程中不会出现过载,有效避免了潮流校验,提高了紧急控制速度。
3 仿真算例
在IEEE39节点系统中对本文方法进行验证,其系统接线图如图1所示。系统各发电机有功出力情况及单位出力变化补偿价格如表1所示,该系统中节点31是平衡机节点,计算该发电机的可增出力量时需留出10 MW的裕度,同时,令各发电机的最小发电量为0 MW,实际电网运行中依情况而定。因篇幅有限,仅示出该系统部分负荷节点的负荷量及单位切负荷补偿价格,如表2所示。
图1 IEEE 39节点系统Fig.1 IEEE 39-bus system
表1 各发电机有功出力情况及单位出力变化补偿价格Table1 Power output and compensation per output-changing for generators
表2 部分负荷节点的负荷量及单位切负荷补偿价格Table2 Load and compensation per load-shedding for some load nodes
将本文方法定义为方法1,将以功率灵敏度为基础以系统总调整量最小为优化目标的方法定义为方法2,将以GSDF为基础以系统总调整量最小为优化目标的方法(即文献[5]中的方法)定义为方法3。在MATLAB中利用优化规划法编程实现3种方法,当潮流转移导致系统中线路过载时,分别利用这3种方法制定相应的紧急控制方案,在PSASP中进行潮流计算验证各紧急控制方案的可行性,比较各方案的总调整量及经济补偿,验证本文所提方法的优越性。
3.1 仿真1
IEEE39节点系统中,线路B2-3和B26-27因故障断开,导致线路B2-1和B25-2过载,其线路功率分别为738.91 MW和505.26 MW,相应的过载率分别为1.23和1.01。发电机节点及部分负荷节点对过载线路的功率灵敏度如表3所示。
分别利用3种方法制定紧急控制方案,各紧急控制方案如表4所示。利用PSASP软件验证各紧急控制方案的可行性,结果表明3种方法均只需调节发电机出力即可消除线路过载,使系统恢复安全运行状态,各方案的总调节量和经济补偿如表5所示。
3.2 仿真2
IEEE 39节点系统中,线路B11-6和B14-4因故障断开, 导致线路 B16-17、B3-4、B17-18、B13-14、B14-15和 B10-13过载,过载线路功率如表6所示。
发电机节点及部分负荷节点对过载线路的功率灵敏度如表7所示。
分别利用3种方法制定紧急控制方案,方法1、2均只需调节发电机出力即可消除线路过载,方法3需要切负荷才能消除所有线路过载。各紧急控制方案如表8所示(负荷节点的调整量即为切负荷量),各方案的总调节量和经济补偿如表9所示。
表3 线路B2-3和B26-27断开时部分节点对过载线路的功率灵敏度Table3 Power sensitivity of partial nodes to overload line when Line B2-3and B26-27are cut off
表4 仿真1的紧急控制方案Table4 Emergency control scheme for Simulation 1
表5 仿真1中各方法的总调节量和经济补偿Table5 Total power regulation and economic compensation for different strategies of Simulation 1
表6 线路B11-6和B14-4断开时的过载线路功率Table6 Power of overload lines when Line B11-6and B14-4are cut off
表7 线路B11-6和B14-4断开时部分节点对过载线路的功率灵敏度Table7 Power sensitivity of partial nodes to overload line when Line B11-6and B14-4are cut off
表8 仿真2的紧急控制方案Table8 Emergency control scheme for Simulation 2
表9 仿真2中各方法的总调节量和经济补偿Table9 Total power regulation and economic compensation for different strategies of Simulation 2
a.由表4和表9可知,方法1虽然在总调节量上略多于方法2,但在总经济补偿代价上有较大优势,因此,以经济补偿最小化为目标得到的紧急控制方案有更大的实用价值。
对比方法2和方法3可知,方法2在总调节量和总经济补偿上都少于方法3。结合表4和表7中这2种方法的紧急控制方案可知,平衡机在紧急控制中发挥了重要作用。因此,利用功率灵敏度矩阵衡量节点注入功率变量对线路功率变量的控制情况,计算平衡机对各线路功率的灵敏度,令平衡机参与到紧急控制中,可以有效减少消除线路过载所需的总调节量和总经济补偿,得到更好的紧急控制方案。和方法3相比,方法1在总控制量和总经济补偿上都有很大优势,再次验证了本文方法的优越性。
b.仿真1中,3种方法均只需调节发电机出力即可消除线路过载,使系统恢复安全运行状态;仿真2中,方法1和方法2均只需调整发电机出力即可消除线路过载,而方法3除了调整发电机出力之外,还需要切负荷才能消除线路过载,证明了使平衡机参与紧急控制的优势。
需要说明的是,仿真2中方法3所得的紧急控制方案中,在系统加出力发电机并未完全用尽其可增出力量的情况下就采取了切负荷措施,这是因为为保证紧急控制中不出现新的线路过载,系统正常线路冗余量限制了部分发电机的加出力量。
c.本文仿真中所采用的补偿价格数据,是在文献[7]的基础上进行修改后得到的,实际电网中依情况而定,仿真所用数据不影响对本文方法优越性的验证。
d.式(3)中M的取值并无规定,只要能够体现负荷节点优先级低于发电机节点的合适取值均可。
4 结论
本文提出一种基于功率灵敏度和经济补偿最小化的线路过负荷紧急控制方法。利用功率灵敏度矩阵反映各节点注入功率变量对线路功率变量的影响,计算平衡机对各线路功率的灵敏度,令平衡机参与到紧急控制中。以经济补偿最小化为目标,利用非线性优化得到满足线路过负荷消除约束、节点可调量约束、系统功率平衡约束及正常线路冗余量约束的紧急控制方案。仿真算例验证了本文方法的正确性和优越性。
本文方法虽有一定的有效性,但是在考虑紧急控制的经济代价时只考虑了经济补偿最小化,如何综合考虑事故风险及紧急控制过程中其他经济代价等因素来全面衡量经济损失是下一步的研究方向。
[1]任建文,魏俊姣,谷雨峰.基于多目标粒子群优化算法的连锁跳闸预防控制[J].电力自动化设备,2016,36(7):53-59.REN Jianwen,WEI Junjiao,GU Yufeng.Preventive control based on multi-objective particle swarm optimization algorithm for cascading trips[J].Electric Power Automation Equipment,2016,36(7):53-59.
[2]程临燕,张保会,郝治国,等.基于综合灵敏度分析的快速控制算法研究[J].电力自动化设备,2009,29(4):46-49.CHENG Linyan,ZHANG Baohui,HAO Zhiguo,et al.Fast control algorithm based on integrative sensitivity analysis[J].Electric Power Automation Equipment,2009,29(4):46-49.
[3]邓佑满,黎辉,张伯明,等.电力系统有功安全校正策略的反向等量配对调整法[J].电力系统自动化,1999,23(18):5-8.DENG Youman,LI Hui,ZHANG Boming,et al.Adjustment of equal and opposite quantities in pairs for strategy of active power security correction of power systems[J].Automation of Electric Power Systems,1999,23(18):5-8.
[4]李刚,王增平,任建文,等.基于图论分区与改进BFS算法搜索安全约束集的防联锁过载控制策略[J].电工技术学报,2012,27(11):219-229.LI Gang,WANG Zengping,REN Jianwen,et al.A control strategy to prevent cascading overload with security constraint set searched by graph theory partition and improved BFS algorithm[J].Transactions of China Electrotechnical Society,2012,27(11):219-229.
[5]程临燕,郝治国,张保会,等.基于内点法消除输电断面过载的实时控制算法[J].电力系统自动化,2011,35(17):51-55.CHENG Linyan,HAO Zhiguo,ZHANG Baohui,et al.Fast elimination of overload in transmission line section based on simplified primal-dual interior point method [J].Automation of Electric Power Systems,2011,35(17):51-55.
[6]赵晋泉,江晓东,张伯明.一种基于连续线性规划技术的在线静态安全校正算法[J].电网技术,2005,29(5):25-30.ZHAO Jinquan,CHIANG Hsiaodong,ZHANG Boming.A successive linear programming based on-line static security corrective control approach[J].Power System Technology,2005,29(5):25-30.
[7]孔祥玉,房大中,钟德成.电力市场稳定安全约束下的经济代价最小数学模型[J].电力系统自动化,2004,28(6):25-29.KONG Xiangyu,FANG Dazhong,ZHONG Decheng.Mathematic models for minimum cost in security management of power market[J].Automation of Electric Power Systems,2004,28(6):25-29.
[8]姜臻,苗世洪,刘沛,等.一种基于粒子群优化算法的转移潮流控制策略[J].电力系统自动化,2010,34(18):16-20.JIANG Zhen,MIAO Shihong,LIU Pei,et al.A particle swarm optimization based power flow transferring control strategy [J].Automation of Electric Power Systems,2010,34(18):16-20.
[9]曲正伟,王京波,张坤,等.考虑不确定性成本的含风电场群电力系统短期优化调度[J].电力自动化设备,2016,36(4):137-144.QU Zhengwei,WANG Jingbo,ZHANG Kun,et al.Short-term optimal dispatch considering uncertainty cost for power system with wind farms[J].Electric Power Automation Equipment,2016,36(4):137-144.
[10]袁晓丹,张会强.多支路开断潮流转移识别及防连锁过载策略研究[J].现代电力,2014,31(5):74-79.YUAN Xiaodan,ZHANG Huiqiang.Study on power flow transferring identification of multi-line tripping and the strategy to prevent cascading overload strategy [J].Modern Electric Power,2014,31(5):74-79.
[11]沈瑜,夏清,康重庆.发电联合转移分布因子及快速静态安全校核算法[J].电力系统自动化,2003,27(18):13-17.SHEN Yu,XIA Qing,KANG Chongqing.Generation Jointly Shift Distribution Factor(GJSDF) and fast security analysis algorithm in economic dispatch[J].Automation of Electric Power Systems,2003,27(18):13-17.
[12]徐岩,郅静.基于功率灵敏度的线路过载划区域紧急控制策略[J].电工技术学报,2015,30(15):60-72.XU Yan,ZHI Jing.A zone-divided emergency control strategy for overload lines based on power sensitivity[J].Transactions of China Electrotechnical Society,2015,30(15):60-72.
[13]陈珩.电力系统稳态分析[M].北京:中国电力出版社,2007:201-219.
[14]任建文,李莎,严敏敏,等.基于潮流跟踪算法的线路过负荷紧急控制策略[J].电网技术,2013,37(2):392-397.REN Jianwen,LI Sha,YAN Minmin,et al.Emergency control strategy for line overload based on power flow tracing algorithm[J].Power System Technology,2013,37(2):392-397.
[15]杨文辉.预防连锁跳闸的关键线路后备保护与紧急控制策略研究[D].北京:华北电力大学,2012.YANG Wenhui.Study on backup protection and control strategy for critincal lines against cascading trips[D].Beijing:North China Electric Power University,2012.