潮流计算的相关问题2010

1§4.5牛顿-拉夫逊法计算潮流有关问题一、初值、收敛性和多值解1.初值：初值选择不好，比较大，破坏了牛顿法的基础，不收敛。选择的原则。2.收敛性：牛顿-拉夫逊法具有平方收敛特性，高斯-塞德尔法、PQ分解法为一阶收敛特性。X23.多值解对于非线性方程组，解的可能性有：•有实际意义的解•有解，但在实际中无意义（PV节点或平衡节点的无功功率超过允许值，平衡节点的有功功率超过允许值；节点的电压过高或过低）对策：调整运行参数，PV节点、PQ节点相互转化•无解，或无实数解给定的网络结构和运行方式不合理；PV节点数目过少对策：调整运行方式，增加PV节点问题很复杂，至今尚未很好解决3二、稀疏矩阵技术1.稀疏矩阵表示法节点导纳矩阵：高度稀疏的N阶复数对称方阵。因此记录矩阵的下三角。用数组表示数组1：记录矩阵对角元素的数值；数组2：记录矩阵非对角元素的数值（按列存储）；数组3：记录矩阵非对角元素的行号；数组4：记录矩阵非对角元素的按行排的位置数；4非对角元素用指针表示，一个指针用结构表示：行号；列号；幅值；角度；指针（指向下一个非零元素）。对角元素用一个一维数组表示。雅可比矩阵：高度稀疏的2N阶实数方阵，其形式对称但数值不对称。其稀疏程度与节点导纳矩阵相同，可根据节点导纳矩阵形成。52.高斯消去法求解牛顿-拉夫逊法潮流计算的修正方程，可以采用矩阵求逆的方法。但是由于潮流计算的雅可比矩阵通常是一个高度稀疏的矩阵，其逆阵则是一个满矩阵，因此用求逆的方法会增加额外的存储单元和计算工作量。而用高斯消去法则可以保持方程组原有的稀疏性，可以大大减少计算所需的内存和时间。63.节点的优化编号静态优化法：按静态联结支路数的多少编号。统计好网络中各节点联结的支路数后，按联结支路数的多少，由少到多，顺序编号。半动态优化法：按动态联结支路数的多少编号。先只编一个联结支路数最小的节点号，并立即将其消去；再编消去第一个节点后联结支路数最小的节点号，再立即将其消去……依此类推。动态优化法：按动态增加支路数的多少编号。不首先进行节点编号，而是寻找消去后出现的新支路数最少的节点，并为其编号，且立即将其消去；然后再寻找第二个消去后出现的新支路数最少的节点并为其编号，再立即将其消去……依此类推。7三、直流潮流计算直流法的特点：简单、计算工作量小、没有收敛性问题，易于快速地处理投入或断开线路等操作。广泛应用于电力系统规划、静态安全分析以及牛顿-拉夫逊法潮流的初值计算等需要大量计算或运行条件不十分理想的场合。直流法的适用范围：110KV以上的超高压线路。直流法的经常处理的问题：处理开断问题，例如，在电力系统规划和电力系统静态安全分析时，需要进行一种所谓N-1校核计算，即对于某一种运行方式要逐一开断系统中的线路或变压器，检查是否存在支路过载情况。8直流法计算潮流的过程电力网中每条支路i-j中通过的有功功率为：根据电力系统的实际条件可做如下假设：1.实际电力系统中输电线路（或变压器）的电阻远小于其电抗，对地电导可忽略不计2.在正常运行时线路两端相位差很少超过20°3.节点电压值的偏移很少超过10%，且对有功功率分布影响不大****2ReRecossinijijijiiijiijijijijijijPUIUyUUUgUUgB9用式子表示：从而可得：各节点的注入功率为与该节点相连各支路功率之和：1.0,1/2.sin,cos13.1ijijijijijijijgbxUU/ijijijijijPbx(cossin)iijijijijijjiPUUGB10令B表示正常运行时电力网节点导纳矩阵的负数，则所有节点注入功率可用矩阵表示为：解方程求出各节点的相角后，可利用前面的式子求出各支路的有功潮流。1PBBP1()11(,)niijijijiijjijjjijijijiiijiiijjijiijijjijiPBBBBBBBBXX11§4.6潮流计算其它相关问题12一、潮流计算的发展历史Gauss法Newton法FDLF法计及非线性法最优潮流法含直流或FACTS元件的潮流Gauss法1、1956年，基于导纳矩阵的简单迭代法参考文献：WardJB，HaleHW．DigitalComputerApplicationsSolutionofPowerFlowPr-oblems．AIEETrans，1956，75，III：398~404该法特点：原理简单、内存需求较少、算法收敛性差2、1963年，基于阻抗矩阵的的算法参考文献：BrownHE，etal．PowerFlowSolutionbyImpedanceMatrixIterativmethod．IEEETransonPowerApparatusandSystems，1963，PAS-82：1~10特点：收敛性好、内存占用量大大增加（限制解题规模）1967年，Newton法参考文献：TinneyWF，HartCE．PowerFlowSolutionbyNewton’sMethod．IEEETransonPowerApparatusandSystems，Nov1967，PAS-86：1449～14601974年，FDLF法参考文献：StottB，AlsacO．FastDecoupledLoadFlow．IEEETransonPowerApparatusandSystems，May/June1974，PAS-93（3）：859～8691、1978年，保留非线性的快速潮流算法参考文献：IwamotoS，TamuraY．AFastLoadFlowMethodRetainingNonlinearity．IEEETrans．PAS．1978．97（5）：1586～15992、1982年，包括二阶项的快速潮流算法参考文献：RaoPSNagendra，RaoKSPrakasa，NandaJ．AnExactFastLoadFlowMethodIncludingSecondOrderTermsinRectangularCoordinates．IEEETrans．PAS．1982．101（9）：3261～3268最优潮流法1、1962年，最优潮流数学模型参考文献：JCarpentier．Contributional’etudeduDispatchingEconomique．Bull．Soc．Fr．Elec．1962．88（10）：1577～15812、1968年，最优潮流的简化梯度法参考文献：DommelHW，TinneyWF．OptimalPowerFlowSolutions．IEEETrans．PAS．1968．87（10）：1866～18763、1984年，最优潮流计算的牛顿算法参考文献：SunDI，etal．OptimalPowerFlowbyNewtonApproach．IEEETrans．PAS．1984．103（10）：2864～2880含直流和FACTS元件的潮流计算1、1976年，交直流潮流计算参考文献：BraunagelDA，KraftLA，WhysongJL．InclusionofDCConverterandTransmisstionEquationsDirectlyinaNewtonPowerFlow．IEEETrans．PAS．1976．95（1）：76～882、1992年，含Facts元件的潮流计算参考文献：GNTaranto，LMVGPinto，MVFPereira．Repres-EntationofFACTSDevicesinPowerFlowEconomicDispatch．IEEETrans．OnPowerSystem，1992，7（1）：572～57613二、特殊性质的潮流计算1、直流潮流这是一种近似算法，不计支路无功潮流，计算速度是所有潮流算法中最快的。应用场合：电力系统规划设计、实时安全分析的预想事故筛选等2、随机潮流这是一种把潮流计算的已知量和待求量都作为随机变量来处理的一种潮流计算方法，也叫概率潮流。计算结果具有概率统计特性（如期望值、方差、概率分布函数等）。3、三相潮流针对三相不对称的系统，已知量和待求量是单相潮流的三倍，建立三相潮流计算模型后，其计算方法类似单相潮流。14二、特殊性质的潮流计算6、谐波潮流谐波潮流计算考虑非线性元件对系统电能质量的影响，除了计算常规潮流计算中的基波潮流外，还要计算高次谐波。4、动态潮流动态潮流是计算系统存在不平衡功率情况下的稳态潮流，这种潮流计算中V节点和平衡节点不是一个概念，V节点只有一个，但是平衡节点有多个，不平衡功率在多台发电机中分配，还可以考虑系统功率不平衡时的频率调节效应。5、开断潮流开断潮流研究的开断包括：输电线路（变压器）开断、发电机开断和负荷开断。15二、特殊性质的潮流计算8、最优潮流在网络结构和参数及系统负荷给定的条件下，确定系统的控制变量u，使得描述系统运行效益的某一给定的目标函数取最小值。7、约束潮流除了满足常规潮流方程外，还要满足更多的约束条件，如节点电压、支路潮流minmaxiiiUUUminmaxijijijPPPminmaxijijijQQQmin()..()0()Objfxsthxggxg16三、潮流计算软件介绍1、国际上几种电力系统分析计算软件包17三、潮流计算软件介绍2、国内用得较多的几种潮流计算软件简介(1)BPA潮流计算程序简介：美国帮涅维尔电力局（BPA，BonnevillePowerAdministr-ation）开发，被中国电力科学院引进吸收，从1984年开始在中国得到推广应用。程序提供两种潮流计算方法：P_Q分解法和牛顿法(2)PSASP潮流计算程序简介：中国电力科学院开发。程序提供五种潮流计算方法：P_Q分解法、牛顿法(功率式)、最佳乘子法、牛顿法（电流式）、P_Q分解法转牛顿法(电流式)(3)PSS/E潮流计算程序简介：美国PTI开发，70年代推向市场，目前已有40个国家200多家公司应用该程序。提供5种潮流计算方法：牛顿法、解耦牛顿法、快速牛顿法、高斯－塞德尔法、改进的高斯－塞德尔法18四、潮流计算实例说明：1）采用西安交通大学发电教研室开发的GraphEdit1.1版2）采用IEEE22节点系统作为算例191、IEEE22节点电网接线202、潮流计算条件设置计算方法：牛顿法初始电压：Vx=1.0Vy=0.0计算精度：0.0001电压限值：Vmax＝1.2Vmin=0.8213、潮流计算过程节点号电压幅值电压相角注入有功注入无功11.0000.0001.5001.00021.0000.0003.0000.00031.0000.0002.8000.00041.0000.0003.0001.80051.0000.0001.5000.00061.0000.0000.0000.00071.0000.000-2.000-0.80081.0000.000-2.000-4.00091.0000.0000.0000.000101.0000.0000.0000.000111.0000.000-2.000-1.000121.0000.0000.0000.000131.0000.0000.0000.000141.0000.0000.0000.000151.0000.0000.0000.000161.0000.0000.0000.000171.0000.0000.0000.000181.0000.000-2.000-1.500191.0000.000-0.500-0.500201.0000.000-2.500-4.800211.0000.000-1.000-1.500221.0000.0000.0000.000(1)迭代前的初值列表（优化编号后）22(2)迭代前雅可比矩阵和迭代后雅可比矩阵观察比较：雅可比矩阵元素