《电力系统分析》复习题1.分别列出下列潮流算法的迭代格式、收敛判据,并从收敛性、计算量和内存占用量比较其算法特点及适用范围。(1)直角坐标的N-R法;(2)极坐标的N-R法;(3)快速解耦潮流算法(P-Q分解法);(4)二阶潮流算法(保留非线性潮流算法);(5)最优乘子法。答:(1)极坐标N-R法:迭代格式:PHNQMLUU1kkkUUU,1kkk。收敛判据:P或Q牛顿潮流算法的特点1)其优点是收敛速度快,若初值较好,算法将具有平方收敛特性,一般迭代4~5次便可以收敛到非常精确的解,而且其迭代次数与所计算网络的规模基本无关。2)牛顿法也具有良好的收敛可靠性,对于对高斯-塞德尔法呈病态的系统,牛顿法均能可靠地敛。3)初值对牛顿法的收敛性影响很大。解决的办法可以先用高斯-塞德尔法迭代1~2次,以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次求得一个较好的角度初值,然后转入牛顿法迭代。(2)直角坐标N-R法:迭代格式:2PHNeQMLfRSU1kkkeee,1kkkfff收敛判据:P或Q特点同极坐标N-R(3)P-Q分解法:迭代格式:'PUB,''QUBU1kkkUUU,1kkk收敛判据:maxiiiPU且maxiiiQU特点:(1)用解两个阶数几乎减半的方程组(n-1阶和n-m-1阶)代替牛顿法的解一个(2n-m-2)阶方程组,显著地减少了内存需求量及计算量。(2)牛顿法每次迭代都要重新形成雅可比矩阵并进行三角分解,而P-Q分解法的系数矩阵B’和B’’是常数阵,因此只需形成一次并进行三角分解组成因子表,在迭代过程可以反复应用,显著缩短了每次迭代所需的时间。(3)雅可比矩阵J不对称,而B’和B’’都是对称阵,为此只要形成并贮存因子表的上三角或下三角部分,减少了三角分解的计算量并节约了内存。由于上述原因,P-Q分解法所需的内存量约为牛顿法的60%,而每次迭代所需时间约为牛顿法的1/5。PS:处理R/X大比值的两种方法:对大R/X比值支路的参数加以补偿(串补并补,串补有电压畸形问题,并补没有更好);对算法加以改进(BX方案有明显优势)。(4)保留非线性潮流算法迭代格式:1001kksxJxyxyyx收敛判据:1maxkkiiiiiyxyx特点:1)保留非线性快速潮流算法采用的是用初值(0)x计算而得的恒定雅可比矩阵,整个计算过程只需一次形成,可用三角分解构成因子表。所以每次迭代所需时间可以节省很多。2)两种算法x的含义不同。牛顿法的()kx是相对于上一次迭代所得到的迭代点()kx的修正量;而保留非线性快速潮流算法的()kx则是相对于始终不变的初始估计值(0)x的修正量。3)保留非线性快速潮流算法达到收敛所需的迭代次数比牛领法要多,但由于每次迭代所需的计算量比牛顿法节省很多,所以总的计算速度比牛顿法可提高很多。4)由于不具对称性质的雅可比矩阵经三角分解后,其上下三角元素都需要保存,和牛顿法的一种方案仅需保存上三角元素相比,此算法所需的矩阵存储量将比要牛顿法增加35%~40%。5)由于利用以初始值计算得到的恒定雅可比矩阵进行迭代,初始值的选择对保留非线性快速潮流算法的收敛特性有很大影响。保留非线性快速潮流算法比牛顿法优越,但与快速解耦法相比:计算速度稍慢,内存相差太大。(5)最优乘子法:迭代格式:100112kksxJxyxyyx收敛判据:1maxkkiiiiiyxyx特点:采用带有最优乘子k的牛顿潮流算法以后,潮流计算不会发散,即从算法上保证了计算过程的收敛性,从而有效地解决了病态潮流的计算问题。而通过k的具体数值,提供了在给定的运算条件下,潮流问题是否存在解的一个判断标志。2.五节点电力系统的节点类型和支路参数如下:节点类型:1PV节点;2、3、4PQ节点;5平衡节点。支路标幺参数(非标准变比在支路首端):1-22-33-44-52-4电阻R0.080.050.09电抗X0.20.30.20.30.35电纳1/2B0.60.40.75变比K1.050.975(1)请写出采用P-Q分解法进行潮流计算时,修正方程的B’和B”矩阵。解:5500511.193.332.86'03.333.33502.8657.86B9.443.112.68''3.116.824.712.684.719.75B(注意:电纳正为容性,与支路感性符号相反;非标准变比需要等效为型等效电路,线路和末端处理不一样)(2)请写出采用直角坐标N-R法进行潮流计算时的修正方程结构。(注:状态变量和运行变量用符号表示,雅可比矩阵中非零元素用“×”表示,零元素用“0”表示。)解:11211222233334444000000000000000000PeUfPeQfPeQfPeQf(3)如果考虑负荷节点的静态电压特性或非标准变比,雅可比矩阵元素有什么不同?答:考虑负荷节点的静态电压特性时,雅各比矩阵元素iiN、iiL将和N-R法的雅各比矩阵对应元素不一样。考虑非标准变比时,当支路i-j是用来调整节点j电压的变压器支路时,雅各比矩阵中与节点j的电压向量所对应的一列元素除了对角元素外,只有一组非零对角元素。3.什么是病态潮流?哪些潮流算法适用于求解病态潮流?答:满足下列某个或多个条件的系统的潮流问题叫做病态潮流:(1)节点间相位角差很大的重负荷系统;(2)包含有负电抗支路(如某些三绕组变压器或线路串联电容等)的系统;(3)具有较长的辐射形线路的系统;(4)长线路与短线路接在同一节点上,而且长短线路的长度比值又很大的系统。牛拉法(N-R)、P-Q分解法的BX方案、保留非线性法和最小化潮流等。4.最优潮流与传统经济调度的区别是什么?答:由于基于协调方程式的经典经济调度方法虽然具有方法简单,计算速度快,适宜于实时应用等优点,但协调方程式在处理节点电压越界及线路过负荷等安全约束的问题上却显得无能为力。以数学规划问题作为基本模式的最优潮流在约束条件的处理上具有很强的能力。最优潮流能够在模型中引入能表示成状态变量和控制变量函数的各种不等式约束,将电力系统对于经济性、安全性以及电能质量三方面的要求,完美地统一起来。5.最优潮流计算的简化梯度法计算步骤。不等约束条件的处理方法。启动输入原始数据给定(0)u0k进行潮流计算即求解得()()()0kkgx,u()kx根据不等式约束被违犯的情况决定惩罚项的组成()()(,)kkWux由式(1-216)计算()kλ根据式(1-217)计算简化梯度()kf()?kf决定修正步长因子()kc修正控制变量()()()()kkkkuucf对越界控制变量加以限制1kk输出计算结果停机是控制变量不等式约束:控制变量的不等式约束比较容易处理,若按照(1)()()kkkuuu对控制变量进行修正,如果得到的()ku使得任何一个(1)ku超过其限值maxiu或者miniu时,则该越界的控制变量就被强制定在响应的界上。函数不等式约束:采用罚函数处理法,即越界不等式约束以罚函数项的形式附加在原来的目标函数(,)fux上,从而构成一个新的目标函数(即惩罚函数),对这个新的目标函数按无约束求极值的方法求解,使得最终求得的节点在满足上述约束条件的前提下能使原来的目标函数达到最小。6.什么是电力系统的状态估计?电力系统状态估计的意义是什么?答:电力系统状态估计:对给定的系统结构及量测配置,在量测量有误差的情况下,通过计算得到可靠地并且位数最少的状态变量值----各母线上的电压相角与模值及各元件上的潮流。电力系统状态估计的意义:(1)提高数据精度,去除不良数据;(2)计算出难以测量的电气量,相当于补充了量测量。(3)状态估计为建立一个高质量的数据库提供数据信息,以便于进一步实现在线潮流、安全分析及经济调度等功能。7.运行状态估计必须具备什么基本条件?状态估计如何实现的?主要完成的功能有哪些?答:实现状态估计需要的条件:1.量测冗余度:量测冗余度是指量测量个数m与待估计的状态量个数n之间的比值m/n。系统冗余度越高,对状态估计采用一定的估计方法排除不良数据以及消除误差影响就越好。冗余量测的存在是状态估计可以实现提高数据精度的基础。2.分析系统可观性:当收集到的量测量通过量测方程能够覆盖所有母线的电压幅值和相角时,则通过状态估计可以得到这些值,称该系统是可观测的。状态估计如何实现:完成的主要功能:进行不良数据的检测与辨识,最终建立完整的电力系统模型;提供补充的测量量;用于负荷预测。8.状态估计与常规潮流计算的区别和联系。答:潮流计算方程式的数目等于未知数的数目。而状态估计的测量向量的维数一般大于未知状态向量的维数,即方程数的个数多于未知数的个数。其中,测量向量可以是节点电压、节点注入功率、线路潮流等测量量的任意组合。两者求解的数学方法也不同。潮流计算一般用牛顿-拉夫逊法求解个非线性方程组。而状态估计则是根据一定的估计准则,按估计理论的方法求解方程组状态估计中的“估计”不意味着不准确,相反,对于实际运行的系统来说,不能认为潮流计算是绝对准确的,而状态估计的值显然更准确。状态估计可认为是一种广义潮流,而常规潮流计算是一种狭义潮流,及状态估计中m=n的特例。9.某电力系统及其测量系统配置如图,图中:“×”表示支路潮流量测;“”表示节点注入功率量测;“U”表示节点电压量测。网络结构处理不良数据检测与辨识可观察性检验状态估计器负荷预计实时数据库有无10.UU134256(1)写出该系统的节点导纳矩阵Y、测量方程的雅可比矩阵H和信息矩阵A的结构(注:矩阵中非零元素用“×”表示,零元素用“0”表示)1515252534344343515152526262126PQ0000PQ0000PQ00000000PQ0000000PQ0000000YHPQ0000000PQ000000U00000000U00000P000::0000000000A00000000:(2)画出与信息矩阵A对应的网络结构。13425610、简明扼要地填写下表:物理量服从分布数学期望方差协方差状态估计误差正态分布011()THRH对角元11()THRH量测估计误差正态分布011()TTHHRHH对角元11()TTHHRHH残差正态分布011()TTWRRHHRHH对角元11()TTWRRHHRHH加权残差正态分布01WWRWR对角元1WWRWR标准化残差标准正态分布0111DWRD目标函数2分布K(测量冗余度)2K无11、状态估计的算法通常有那些?各种算法的优缺点是什么?答:状态估计算法主要有:加权最小二乘法,快速解耦状态估计,支路潮流状态估计法三种方法优缺点:1)加权最小二乘法:估计质量好,收敛性好,但内存需要大,计算时间长,适用于小型电力系统。2)快速分解法: