电力系统最小负荷功率裕度的快速求取.docx

电力系统最小负荷功率裕度的快速求取王惠喜 胡林献哈尔滨工业大学摘要：从系统初始运行点出发，以雅可比矩阵零奇异值对应的左奇异向量为初始负荷增长方向，用改进连续潮流法搜索电压崩溃点，再以该崩溃点雅可比矩阵零奇异值对应的左奇异向量作为新的负荷增长方向重新计算，直到负荷增长方向与崩溃点处崩溃点曲面的法线方向重合为止，从而得出了一种最危险的负荷增长方式和最小负荷功率裕度的快速求解方法，IEEE14节点系统实例计算验证了该方法的可行性。关键词：连续潮流 最危险负荷增长方向 最近崩溃点 最小负荷功率裕度Fast Solving Method of the Minimal Load Margin in Power SystemWang Huixi Hu Linxian Abstract: In this article, from the initial operating point I searched the voltage collapse point with the improved continuation power flow method, using the left singular vector corresponding to the zero singular value of the Jacobi matrix as the initial load increasing direction, and then used the left singular vector of the voltage collapse point just gotten as the new load increasing direction to recalculate till the load increasing direction coincided with the normal direction of the collapse hood face at the new voltage collapse point. Accordingly I put forward a fast solving method for the worst load increasing direction and the minimal load margin, and validate the feasibility of this method with the calculation of the IEEE14 system in the end.Keywords: continuation power flow the worst load increasing direction the closest collapse point the minimal load margin 71引言近年来，随着工业生产的高度发展和经济、环境各方面的改变，现代电力系统发生了许多变化，规模越来越大，负荷越来越重，大容量发电机组、交流超高压线路以及直流输电和新型控制装置等得到普遍应用，但受输电线路成本及环境保护的限制，往往不能建成相应强大的输电网，负荷中心地区发输电设施的建设速度远小于负荷需求的增长速度，能源供应中心与负荷中心呈现出逐步分离的趋势，使电网结构相对脆弱，所有这些都给电力系统的安全运行带来严重威胁，其中之一便是由于系统不能维持受端电压水平而引起的电压失稳或电压崩溃事故，因此提供能够准确计算系统最危险的负荷增长方向及最小负荷功率裕度的方法具有重要的意义。为了可靠、快速地计算出系统临界点，进而计算功率裕度，人们研究了不少方法，如通过等效变换改变系统结构参数的重负荷导纳法；在常规潮流方程的基础上增加参数方程、求解增广潮流方程的连续潮流法1和崩溃点直接计算法2,3等。但这些方法主要研究给定负荷增长方向下的功率裕度求取，而在实际电力系统中，运行人员主要关心最危险的负荷增长方向和目前运行点与最近电压崩溃点的距离，以便进行预防控制。不同负荷增长方式对应不同的崩溃点，这些崩溃点构成了一个崩溃点曲面，该曲面与当前运行点最近的点即为最近崩溃点，这两点的连线方向也就是最危险的负荷增长方向，最近崩溃点的裕度即最小负荷功率裕度。理论证明，如果初始运行点指向崩溃点的方向和该崩溃点在崩溃点曲面的法线方向重合，那么该崩溃点就是最近崩溃点4。本文给出了崩溃点曲面上法线的求解方法，证明了电压崩溃点雅可比矩阵零奇异值对应的左奇异向量就是该崩溃点在崩溃点曲面上的法向量。在此基础上提出一种用改进连续潮流法快速求取电力系统最危险的负荷增长方式及最小负荷功率裕度的方法。该方法由系统初始运行点出发，按雅可比矩阵零奇异值对应的左奇异向量方向增长负荷，然后利用连续潮流法求取这个方向对应的电压崩溃点，再以该崩溃点雅可比矩阵零奇异值对应的左奇异向量作为新的负荷增长方向重新计算，直到负荷增长方向与崩溃点的法线方向重合为止。2最短路径算法原理电力系统的初始潮流方程为 （1）式中：；x是系统的状态向量，即各节点电压的幅值和角度；是系统的注入量，即负荷的有功、无功和发电机的有功。令，是原潮流方程的雅可比矩阵。当注入量沿某一方向增长至雅可比矩阵奇异，系统即到达电压稳定临界点，这时的参数向量为，状态向量为，即（，）奇异，而的集合组成N维空间的超平面S。在S上，。问题是给定一个系统初始运行点，如何在S超平面上找到参数向量，使得与之间的距离 (定义为系统电压稳定裕度) 达到极小。假设S为在邻域充分光滑的超平面，可得S面在的一个法向量，=W*T，式中，W*T是零奇异值对应的左奇异向量，将规格化成=1，就是使为最小的负荷增长方向。对任何系统，寻找从起始负荷点到奇异超平面S最短距离的一般过程为：1）在系统初始运行点假定一个初始负荷增长向量，且规格化为=1，也可由=求得。2）不断沿的方向增加，直到奇异，即确定，使得在S超平面上。3）设=，且=1。4）循环1）、2）、3）步，直到收敛于，这时=，是最危险的负荷增长方向向量，是最临近的电压崩溃点，即为最小负荷功率裕度。3 算法的数学推证该算法的推证可化为以下两部分：1）为S在点的法向量。系统的潮流方程可写为式（1），则有，当处于稳定临界点时雅可比矩阵奇异，rankspan我们知道组成S面上的切向量空间，由左奇异向量定义有W*TW*T=0这里是矩阵的零奇异值，即W*T，所以W*T为S面上的法向量。2）从沿在S面的法向量方向到S面的距离是到S面的极小距离（在本算法中已排除是极大值的可能性）。设距离为d, 则有T,由极值条件T，i=1,N-1即垂直于S面的切空间为S面的法向量时，到S面的距离极小，此时即为。因此，系统接近电压崩溃时，在W*T方向上离稳定边界最近，在此方向上的稳定裕度也最小，而按W*T的反方向控制系统可使系统以最快的速度远离崩溃点，达到最优控制的目的。4 初始方向的确定在系统正常运行时，雅可比矩阵J非奇异 ，将式（1）线性化，由奇异值分解法5可得，式中：为的奇异值；，分别为对应的左、右奇异向量。在重负荷情形下，雅可比矩阵接近奇异，其最小奇异值接近零，设为，反映了雅可比矩阵接近奇异的程度，所以系统响应由该最小奇异值和它对应的左、右奇异向量和决定，因此有T若令 则有 即当负荷按最小奇异值对应的左奇异向量方向增长时，状态变量变化最大，所以为初始点附近最危险的负荷增长模式，可把它作为从初始点到最近电压崩溃点的负荷增长方向。5 改进连续潮流法我们知道，在电压崩溃点处，弱节点的电压值一般在0.7 pu以下，为了能快速地接近电压稳定临界点，减少中间运行点的求解个数，本文以弱节点电压作为连续变量，进行大步长下降。由于是以弱节点电压作为连续变量，即使采用大步长一般也不会导致潮流解的发散。为了确定哪个节点属于弱节点，可在初始运行点进行一次以负荷变化因数为连续变量的潮流计算，电压变化量最大的点即为弱节点。大步长下降后，为获得精确的临界点，仍需以弱节点电压作为连续变量，采用小步长追踪。小步长可由下式确定：step=k1e其中， 为控制步长计算的变量，可自行定义，本文分别取0.08，15.0，0.000 2；为系统内幅值最小的节点电压幅值。6 计算步骤本文所提算法的步骤为：1）计算系统的初始运行点，再将所有的PV节点转化为PQ节点，计算时令这些节点的无功输出等于它们的输出上限；2）求解初始点处雅可比矩阵最小奇异值对应的左奇异向量，并将其规格化，作为初始负荷增长方向；3）用改进连续潮流法求出沿方向的临界稳定点及初始运行点到该点的距离；4）求所得临界点处雅可比矩阵最小奇异值对应的左奇异向量，并将其规格化后作为新的负荷增长方向；5）若（为事先设定的一个小数），转6），否则，令=，转3）；6）判断最后得到的临界点处各无功源节点的端电压是否低于指定的电压，若不低于，把端电压高的节点转回PV节点回到3）重新计算。7 算例本文以IEEE14节点系统为例进行验证。由于在重负荷下才关心最小负荷功率裕度，将所有节点的负荷都增加1.5倍进行计算，以=0.001作为收敛条件。本文分别用常规方法和所提出的方法进行计算，以作比较。7.1 常规方法任意选取一个负荷的初始增长方向，本文令负荷节点的有功、无功保持初始工作点时的功率因数和各节点间的比例不变同步增长，多余的负荷由各发电机根据当前出力多少，按比例分配，平衡机不限制。结果如表1所示。表1 常规方法求最小负荷功率裕度初始点与崩溃点距离中间点个数崩溃点10.168 318崩溃点20.142 418崩溃点30.142 218崩溃点40.1422187.2 本文方法用本文所提改进连续潮流法搜索崩溃点，结果如表2所示。表2 改进连续潮流法求最小负荷功率裕度初始点与崩溃点距离中间点个数崩溃点10.14357崩溃点20.14227崩溃点30.14227最近崩溃点处PV节点2、3、6、8的电压为0.917 8，0.799 3，0.702 4，0.780 6，均小于指定的电压，简化条件成立。求得初始方向与最危险负荷增长方向的角度差为=8.01°，可见初始方向的预估比较准确，切实可行。另外比较表1和2可知，常规法搜索了4个崩溃点，每个崩溃点经过了18个中间点，共计算了72次潮流；本文方法搜索了3个崩溃点，每个崩溃点经过了7个中间点，共计算了21次潮流，大大减少了计算量，提高了效率。8 结论本文证明了最危险的负荷增长方向即为最近电压崩溃点的法向量，提出的以雅可比矩阵最小奇异值对应的左奇异向量为初始负荷增长方向的方法提高了搜索最近崩溃点的效率，形成了一种最危险的负荷增长方式和最小负荷功率裕度的快速求解方法，IEEE14节点系统实例计算验证了本文方法的正确性。参考文献1 周双喜，冯治鸿，杨宁.大型电力系统PV曲线的求取J.电网技术，1996，20（8）：4-82 王成山，江伟，江晓东.一种新的电力系统鞍型分叉点计算方法J.中国电机工程学报，1999，19（8）：20-243 郭瑞鹏，韩祯祥.电压崩溃临界点计算的改进零特征根法J.中国电机工程学报，2000，20（5）：63-664 周双喜.电力系统最临近静态电压失稳点的确定.清华大学学报(自然科学版)，1997，37（S1）：88-905 周双喜，朱凌志，郭锡玖等.电力系统电压稳定性及其控制.中国电力出版社，2004：211-213