考虑负荷不确定性的区间潮流计算方法 论文

1、考虑负荷不确定性的区间潮流计算方法考虑负荷不确定性的区间潮流计算方法 论文论文 关键字：迭代 结果 潮流 情况 得到 负荷 区间 方法 计算 考虑负荷不确定性的区间潮流计算方法 裴爱华 1, 刘明波 1, 张 弛 2(1.华南理工大学电力学院, 广州 510640; 2.广东省广电集团公司调度中心, 广州 510600) 摘 要: 潮流计算是进行电力系统分析的重要工具，而实际工程实践中因为负荷等量的不确定因素的存在，使得不确定潮流问题有其重要的实用性。将区间数学引入求解不确定负荷下的潮流问题，采用区间分析的方法进行潮流计算方程组的求解。以区间量为基础的潮流计算模型出发，采用 KrawczykM

2、oore 区间迭代法求解，推出了求解方法。并以 Ward & Hale 6 节点试验系统为算例，进行不确定负荷潮流计算。 关键词: 潮流计算; 不确定性负荷; 区间分析Power Flow Calculation with Unce rtain LoadPEI Aihua1, LIU Mingbo1, ctric Engineering,South China University of Technology, Guangzhou 510641, China; 2.Power Dispatching and Communication Center,Guang Dian Group L

3、td., Abstract: Power flow calculation is a fundamental tool for all kinds of studies in power systems. However, in practical system, some original data may be uncertain such a s the load. Hence, calculating power flow with uncertain loads i s important and practical. This paper presents a method to

4、solve the power flow problem with uncertain loads using interval arithmetic. Based on interval mathematics, an expression of the power flow model with interval variables, is propose d and the power flow calculation with flow result Key words: power flow; uncertain load; interval analysis 1 前言 处理负荷的不

5、确定性一直是潮流计算中的一个难点。考虑负荷不确定性的潮流计算目前主要有三类方法：1）随机潮流法,对于随机的信息，利用概率的方式来处理1，2；2）模糊数学法, 利用模糊数建立配电网潮流计算的模型，利用模糊隶属函数处理一部分不确定信息3；3）区间分析方法, 因为负荷预测的结果在一定范围内准确，则其各个量都为区间，得到的结果也必然在一定区间内，因而采用区间数学来求解非常适用。 区间分析方法最早由 Moore 于 1966 年提出3，迄今已有 30 多年的历史，已成为计算数学的一个重要分支。区间分析方法已在电力系统潮流计算中得到了应用4，5。 本文将区间分析方法应用于求解不确定负荷情况下的潮流问题，直

6、接用区间迭代方法求解潮流非线性方程组，得到节点电压幅值和相角的上下限。区间分析方法与随机潮流法相比，计算量小，与模糊数学法相比，不存在模糊隶属函数的人为给定问题，因而结果不受到任何人为影响。区间数学与区间算法 区间分析方法作为一种可能性方法，用区间数来描述参数的不确定性，并已建立一套完整的区间运算方法和规则3，6。 2.1 区间与区间数 对于给定的实数对。则闭有界数集合 X 就称为有界闭区间。其中，x 称为区间 X 的下端点，x 称为区间 X 的上端点。 若区间 X 的上、下端点相等，即=，定义 X=，为点区间。R 上所有闭界区间的集合，记为 I(R)。 论文考虑负荷不确定性的区间潮流计算方法

7、来自 若y，则定义 XY；若 X 的中点为 0，即=，则称为对称区间。 2.2 区间数的基本运算规则 2.3 非线性方程组的区间算法 引进映射 fD RnRn，则方程组可以简单地表示为 点迭代法的简单的迭代程序为 对于算子 的选择，Moore 首先给出的区间 Newton 算子为 其中：V(X)表示包含区间矩阵 F(X)的逆的区间矩阵，F为 f具有包含单调性的区间扩展，f为 J(X)，是 f(x)的雅可比矩阵。 因为这种方式每一步都要求 V(X)，即需要对区间矩阵求逆，因而这种方法并不是一种实用化 的方法。Krawczyk 为解决区间 Newton 法需要对区间求逆的缺陷，提出了改进的算子 即

8、为 Krawczyk 算子。其中：y 为函数 f(x)的零点，一般取初始迭代区间的中点；Y 为任意 n*n 阶非奇异矩阵；F(X)为 f(x)的雅可比矩阵的具有包含单调性的区间扩展。Krawczy k 算子中 Y 的选择可以有多种方法，若取 这种形式得到的 K 算子称为 Krawczyk 算子的 Moore 形式。3 区间潮流模型及算法 3.1 考虑负荷不确定性的区间潮流模型 设系统有 n 个节点，其中，节点 1，2，m 为 PQ 节点，节点(m+1)，(n-1)为 PV 节点，节点 n 为平衡节点。则潮流方程可写为 3.2 区间潮流算法 由计算模型式（12）和（13）得到的方程组为(n+m-

9、1)阶方程组，为了方便表述，将其表示为 首先求 Krawczyk-Moore 算子，给定初始区间为 其中：V 为电压幅值，为 m 维向量； 为电压相角，为(n-1)维向量；下标 max、分别表示上、下限。取中点 m(X)做为近似零点 x*，h(x)在此点的值为h(x*)。 方程组（14）中，h(x)的雅可比矩阵的具有包含单调性的区间扩展为 非奇异矩阵，否则必然会引起 T 当中 的元素异常。为了防止为奇异，可以在设置初始区间时尽量使上、下限不要关于 0 对称，这样就能避免该区间矩阵的元素的中点为 0。 这样，就得到 Krawczyk-Moore 算子为 即可得到方程组(15)(17)的所有的解。

10、 3.3 求解过程中实际问题的处理 3.3.1 初始区间的给定 电压幅值 V 由节点原有的上、下限值直接得到。除平衡节点外的其它节点电压相角 在角度90，90内。在实际计算中，可根据确定负荷的潮流计算的信息，将角度缩小到确定潮流解的附近区间，这样就能使迭代所需的次数减少，降低计算量。 3.3.2 终止条件的设定 在按照式（18）进行区间迭代的过程中，满足下列条件之一需要终止迭代： 1）迭代得到的区间量 X 的宽度足够窄，这时可认为迭代的点已经很接近方程组的零点，可以输出到结果集合中。 2）得到的区间量 X 代入到原潮流方程中，若其最大区间潮流偏差的区间宽度已经很小，则认为得到的区间已经足够达到

11、方程的零点，输出结果。 3）迭代仍然可以继续进行，但是继续迭代不会使区间再缩小，或者缩小得很慢，则认为零点在得到的区间内，此时应该再进行讨论考察是否需要对分区间。 3.3.3 对分区间的思想 当迭代进行到一定程度，X 还是宽度较大，但每次迭代不能有效缩小其宽度时，可以进行区间的对分。即选定其中的一维量，将其一分为二，这样原来的区间向量就分为两个，再分别进行迭代求解。区间对分的方法是区间数学求解过程中缩小区间宽度的有效方法。 3.3.4 多个结果的分析处理的问题 经过了区间对分，就会使得到的解向量的个数变多。在使用式（18）迭代的过程中，就将会出现空集，则认为该区间内没有零点，该区间自动去除。最

12、终输出的结果需要按照以下依据进行分析处理，满足要求才能作为潮流计算的合理解。分析处理以两个量作为依据： 1）区间量 X 代入到原潮流方程中，得到其最大区间潮流偏差的区间宽度。 2）区间量 X 的中点代入到原潮流方程中，得到其近似最大潮流偏差。 如果这两个量都很小，则该结果 X 可以输出为最后潮流计算的结果。结果分析与讨论 基于区间方法的潮流计算程序由 C+语言实现，在 Visual C+6.0 环境下编译。本文以 Ward & Hale 6 节点试验系统为例进行计算，采用负荷有微小不确定性时的区间潮流与确定负荷状况下的潮流结果进行比较，验证该算法的正确性。Ward & Hale

13、 6 节点试验系统的接线图如图 1 所示，包括 2 台发电机、2 台变压器、2 台并联电容器、4 个负荷点、7 条支路。本文确定潮流计算的程序参考文献7-9的思路设计。 4.1 考虑负荷不确定性的区间潮流计算结果与确定潮流计算结果的比较 图 1 给出了所计算系统的确定负荷状态，当负荷不确定时，我们不能确切地知道每个负荷点 的负荷大小，但是我们知道负荷在一定的范围内变化，本文所模拟的变化幅度为2.5%。该 系统在确定负荷情况下和在不确定负荷情况下潮流计算结果为表 1 和表 2 所示。 表 1 给出了在确定负荷情况下的电压幅值和相角的计算结果以及在不确定负荷情况下区间法 计算出的电压幅值和相角的区

14、间值。表 2 给出了分析出的网络有功标么值。平衡机的有功，在确定负荷情况下 PG1=0.937 932。在不确定负荷情况下 区间潮流算法迭代求解的时间，与终止条件当中的终止区间宽度有很大的关系，如果终止区间过窄，则在进行区间迭代的过程中，需要多次重复对区间进行对分，导致计算量大大增加。甚至可能到一定的步骤以后，式（18）的迭代总不能将区间缩小，中途不能分离出所有不可行的区间，导致最终出现无限的循环。在本文的计算中终止宽度设置为初始区间宽度的 1/3。 4.2 考虑负荷不确定性的区间潮流计算结果的验证与区间法进行潮流计算的意义 准确地验证区间潮流结果的做法可以采用确定潮流的仿真计算来实现，仿真的

15、方法有很多种，其本质都是把各节点的负荷功率的各个量的变化区间进行细分，然后把各种情况进行排列组合，分别计算出各种情况的确定潮流，从所有的情况中总结出网络潮流量的变化上限和下限。因为网络结构的复杂，各种因素的相互影响，这种仿真计算与仅仅针对所有负荷的上下界各作一次确定性潮流计算所得的结果是存在本质性差异的。显然，考虑负荷不确定性的区间潮流计算，由负荷的可能波动的区间出发，计算得到网络潮流的区间结果，与确定潮流计算存在着根本性的不同。 为了验证本文所提出的考虑负荷不确定性区间潮流计算的正确性，可以采用文献5提出的仿真计算方法。该算法的思路为：若计算的系统有 n 个负荷点，则首先分别产生 2n 个随

16、机数，然后利用这些随机数对负荷点的有功和无功量波动区间分别进行随机取值，这样就取出任意的一组负荷可能情况，进行确定潮流计算；然后再产生一组 2n 个随机数对负荷进行随机取值，再进行确定潮流计算，计算尽量多次，最后从所有的情况中总结出网络潮流量的变化上限和下限。这种仿真方法涵盖的负荷可能情况广泛，仿真符合实际情况，并且实现简便。 表 3 给出了在不确定负荷情况下，采用仿真计算分析网络潮流的网络电压幅值和相角最终结果以及区间潮流算法计算出的电压幅值和相角的区间值。可以看到，仿真计算与区间潮流计算 结果基本一致，很好地验证了区间潮流计算的正确性。结论 本文从以区间量为基础的潮流计算模型出发，采用 K

17、rawczykMoore 区间迭代法求解，得到节点的电压幅值、相角上下限，以及网络功率分布,比确定性的潮流算法更符合工程实际情况，因而具有广泛的实用前景。区间数学中所处理的变量是区间，即有界闭域，因而在处理一些问题时，若其变量在一固定闭区域内，则区间数学很适用。区间 Newton 法的优点是：具有全局收敛性，既能给出近似的解，又能给出解的误差估计。因而对于在一些含有区间量的问题上的建模和求解是非常适用的。本文的计算结果有效验证了本算法的正确性。但是区间运算有其复杂性，一般区间迭代程序的计算量远大于点迭代，花费机器时间增加，因此有必要在算法方面继续探讨。参考文献 1Dopazo J F,Klit

18、in O A,Sasson A M.Stochastic load flow methodJ. IEEE Trans,1975,94（5）:299-309 2Wang Z,Alvarado F L.Interval arithmetic in power flow analysisJ .IEEE Trans on Power Systems,1992,7(3):1341-1349 3Moore R E.Interval AnalysisM.Englewood Cliffs,New Jersey：，1966 4王守相,王成山,刘若沁.基于模糊区间算法的配电网潮流计算J.电力系 统自动化.2000,24（20）:19-22 5Biswarup Das.Radial distribution system power flow using interval arithmeticJ.International Jou