P-Q分解法潮流计算

合集下载

(完整word版)PQ分解法计算潮流

一、PQ 分解法的原理P —Q 分解法是牛顿-拉夫逊法潮流计算的一种简化方法。

P-Q 分解法利用了电力系统的一些特有的运行特性,对牛顿-拉夫逊法做了简化，以改进和提高计算速度。

的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻，则系统母线电压幅值的微小变化对用功功率的改变影响很小。

同样,母线电压相角的的改变对无功功率的影响较小.因此，节点功率方程在用极坐标形式表示时。

它的修正方程式可简化为：00P H Q L U U θ∆∆⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥∆∆⎣⎦⎣⎦⎣⎦将P 、Q 分开来迭代计算，因此大大地减少了计算工作量.但是H 、L 在迭代过程中仍将不断变化，而且又都是不对称矩阵。

对牛顿法的进一步简化。

为把上式中的系数矩阵简化成迭代过程中不变的对称矩阵。

在一般情况下线路两端的电压相角ij θ是不大的,因此可以认为：cos 1sin ij ij ijijG B θθ≈2ii ii Q U B考虑到上述关系,可以得到：ij i ij j ij i ij jH U B U L U B U ==节点的功率增量为：11(cos sin )(sin cos )ni is i j ij ij ij ij j ni is i j ij ij ij ij j P P U U G B Q Q U U G B θθθθ==∆=-+∆=--∑∑P —Q 分解法的特点：以一个n-1阶和一个n —m —1阶线性方程组代替原有的2n —m —1阶线性方程组；修正方程的系数矩阵B'和B”为对称常数矩阵，且在迭代过程中保持不变；P —Q 分解法具有线性收敛特性,与牛顿—拉夫逊法相比，当收敛到同样的精度时需要的迭代次数较多。

二、程序说明1.数据说明Branch1。

txt：支路参数矩阵第1列为支路的首端编号;第2列为支路的末端编号（首端编号小于末端编号）；第3列为之路的阻抗;第4为支路的对地容抗；第5列为支路的变比；第6列为折算到那一侧的标志Branch2。

P-Q分解法潮流计算解读

P-Q分解法的特点和性能分析
(1) 用一个n-1阶和一个m阶的线性方程组代替了牛顿法的n-1+m阶线性方程组，显著地减少了内存需求量及计算量。
(2)系数矩阵B’和B’’为常数矩阵。因此，不必像牛顿法那样每次迭代都要形成雅可比矩阵并进行三角分解，只需要在进入迭代过程以前一次形成雅可比矩阵并进行三角分解形成因子表，然后反复利用因子表对不同的常数项△P/V或△Q/V进行消去回代运算，就可以迅速求得修正量，从而显著提高了迭代速度。
在B'中尽量去掉那些对有功功率及电压相角影响较小的因素，如略去变压器非标准电压比和输电线路充电电容的影响；在B"中尽量去掉那些对无功功率及电压幅值影响较小的因素，如略去输电线路电阻的影响
即B’的非对角和对角元素分别按下式计算：
B”的非对角和对角元素分别按下式计算：
其中rij和xij分别为支路的电阻和感抗，bi0为节点i 的接地支路的电纳。（BX法）
由图2-3可以看出，牛顿法在开始时收敛得比较慢，当收敛到一定程度后，它的收敛速度就非常快，而P-Q分解法几乎是按同一速度收敛的。如果给出的收敛条件小于图中A点相应的误差，那么P-Q 分解法所需要的迭代次数要比牛顿法多几次。可以粗略地认为P-Q分解法的选代次数与精度的要求之间存在着线性关系。
(3)系数矩阵B’和B’’是对称矩阵。因此，只需要形成并贮存因子表的上三角或下三角部分，这样又减少了三角分解的计算量并节约了内存。
P-Q分解法的收敛特性
P-Q分解法所采取的一系列简化假定只影响了修正方程式的结构，也就是说只影响了迭代过程，并不影响最终结果。因为P-Q'分解法和牛顿法都采用相同的数学模型式，最后计算功率误差和判断收敛条件都是严格按照精确公式进行的，所以 P-Q分解法和牛顿法一样可以达到很高的精度。

P-Q分解法潮流计算方法改进综述

P-Q分解法潮流计算方法改进综述摘要：本文介绍了P-Q分解法潮流计算方法的数学模型，简化假设及特点，总结了P-Q分解法在低压配电网络中，随着支路R/X比值的增大所带来的迭代次数增大和不收敛性的解决方法，及该方法在不同假设条件下收敛性，并提出了自己的见解。

关键词： P-Q分解法；收敛性；大R/X比支路1 潮流计算的数学模型P-Q分解法又称为快速解耦法，是基于牛顿-拉夫逊法的改进，其基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，把有功功率误差作为修正电压向量角度的依据，把无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开进行【1】。

对一个有 n 个节点的系统，假定第1个为平衡节点，第 2~m+1号节点为PQ节点，第m+2~n号节点为PV节点，则对于每一个PQ或PV节点，都可以在极坐标形式下写出一个有功功率的不平衡方程式：这些假设密切地结合了电力系统的某些固有特点，作为电力系统潮流计算广泛使用的一种算法，P-Q分解法无论是内存占用量还是计算速度方面都比牛顿-拉夫逊法有了较大的改进，主要反映在以下三点：① 在修正方程式中，B’和B’’二者的阶数不同。

B’为n-1 阶，B ‘’为m阶方阵，简化了牛顿法的一个n+m-1的方程组，显著减少了方程组的求解难度，相应地也提高了计算速度。

②用常系数矩阵B’和B’’代替了变系数雅可比矩阵，而且系数矩阵的元素在迭代过程中保持不变。

系数矩阵的元素是由导纳矩阵元素的虚部构成的，可以在进行迭代过程以前，对系数矩阵形成因子表，然后反复利用因子表对不同的常数项△P/V 或△Q/V进行前代和回代运算，就可以迅速求得电压修正量，从而提高了迭代速度，大大地缩短了每次迭代所需的时间【2】。

③用对称的B’和B’’代替了不对称的雅可比矩阵，因此只需要存储因子表的上三角部分，这样减少了三角分解的计算量和内存【2】。

3 P-Q分解法的收敛性改进在各种文献中，都有对P-Q分解法从不同方面提出了讨论和改进，有些是对硬件的改进，如使用并行算法和相应的并行软件来替代原来的串行处理，有些是对算法程序做出了改进，方法众多，不在此累述。

第四节PQ分解法潮流计算

第四节 PQ 分解法潮流计算一、PQ 分解法的基本方程式60年代以来N —R 法曾经是潮流计算中应用比较普遍的方法，但随着网络规模的扩大（从计算几十个节点增加到几百个甚至上千个节点）以及计算机从离线计算向在线计算的发展，N —R 法在内存需要量及计算速度方面越来越不适应要求。

70年代中期出现的快速分解法比较成功的解决了上述问题，使潮流计算在N —R 法的基础上向前迈进了一大步，成为取代N —R 法的算法之一。

快速分解法（又称P —Q 分解法）是从简化牛顿法极坐标形式计算潮流程序的基础上提出来的。

它的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻值，则系统母线电压副值的微小变化V ∆对母线有功功率的改变P ∆影响很小。

同样，母线电压相角的少许改变θ∆，也不会引起母线无功功率的明显改变Q ∆。

因此，节点功率方程在用极坐标形式表示时，它的修正方程式可简化为：⎥⎦⎤⎢⎣⎡∆∆⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∆∆V V L H Q P /00θ （4—19）这就是把2（n —1）阶的线性方程组变成了两个n —1阶的线性方程组，将P 和Q 分开来进行迭代计算，因而大大地减少了计算工作量。

但是，H ，L 在迭代过程中仍然在不断的变化，而且又都是不对称的矩阵。

对牛顿法的进一步简化（也是最关键的一步），即把（4—19）中的系数矩阵简化为在迭代过程中不变的对称矩阵。

在一般情况下，线路两端电压的相角ij θ是不大的（不超过10○～20○）。

因此，可以认为：⎭⎬⎫<<≈ij ij ij ij B G θθsin 1cos （4—20）此外，与系统各节点无功功率相应的导纳B LDi 远远小于该节点自导纳的虚部，即 ii iiLDi B V Q B <<=2 因而 ii i i B V Q 2<< （4—21）考虑到以上关系，式（4—19）的系数矩阵中的各元素可表示为： ij j i ij B V V H = （i,j=1,2,………，n-1）（4—22）ij j i ij B V V L = （i,j=1,2,……………，m ） (4—23)而系数矩阵H 和L 则可以分别写成：⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=------------11,1122,1111,1111,222222121211,1121211111n n n n n n n n n n n n V B V V B V V B V V B V V B V V B V V B V V B V V B V H =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------1211,12,11,11,222211,11211121n n n n n n n n V V V B B B B B B B B B V V V =11D D BV V （4—24）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=m mm m m m m m m m m m m V B V V B V V B V V B V V B V V B V V B V V B V V B V L 22122222212121121211111 =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡m mm m m m m m V V V B B B B B B B B B V V V2121222211121121=22''D D V B V （4—25）将（4—24）和（4—25）式代入（4—19）中，得到[][][][][]θ∆'-=∆11D D V B V P[][][][]V B V Q D ∆-=∆''2用[]11-D V 和[]12-D V 分别左乘以上两式便得：[][][][][]θ∆-=∆-111'D D V B P V （4—26）[][][][]V B Q V D ∆-=∆-''12 （4—27）这就是简化了的修正方程式，它们也可展开写成：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆----------1122111,12,11,11,222211,11211112211n n n n n n n n n n V V V B B B B B B B B B V P V P V P θθθ（4—28）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆m mm m m m m m mV V V B B B B B B B B B V Q V Q V Q 212122221112112211 （4—29）在这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部，因而系数矩阵是对称矩阵，且在迭代过程中保持不变。

基于MATLAB软件的P-Q分解法潮流计算

基于MATLAB软件的P-Q分解法潮流计算摘要电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗。

在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性，可靠性和经济性。

所以，电力系统潮流计算是进行电力系统故障计算，继电保护整定，安全分析的必要工具。

随着电力系统网络的急剧扩大和不断复杂，运用手算进行潮流计算已经不现实。

但是，伴随着计算机技术的飞速发展，基于计算机的潮流计算也就应运而生了。

这样，通过潮流计算，实现对系统的分析成为可能。

本文结合潮流计算的三个基本要求，紧跟该领域的发展，介绍了基于MATLAB软件P-Q分解法潮流计算的程序，该程序用于粗略的计算中小型电力网络的潮流，实现对其的分析。

本文所设计的程序，在计算中，所用的算法通俗易懂并对以往的主流算法做了一些改进，提高了计算速度。

同时，该程序采用了GUI人机对话，将Excel表格、TXT文档与MATLAB程序紧密联系起来，使输入输出界面更加人性化。

关键词：电力系统潮流计算；P-Q分解法；MATLAB软件Power flow calculation of P-Q mode basedonMATLAB softwareAbstractPower flow calculation is one of the important calculations which are to study the operation of power system steady state analysis.It isbased on the given operating conditions and system wiring to identify the various parts of the power system operating state: the buses voltage, the stream components power, system power loss. both power system planning design and operation of existing power system mode of study are need to use the power flow calculation to quantitatively compare the program or run mode power supply reasonable, reliability and economy. Therefore, the power flow calculation is an essential tool for a calculation of power system faults, protection setting, security analysis. withthe rapid expansion of power system network and continuing to be more complex, using hand calculation for flow calculation has been unrealistic.But ,with the celerity development in computer technology, computer-based power flow calculation has also emerged. In this way, It is possible to analysis power system through the power flow calculation.Based on the three basic requirements of power flow calculation and followed by the development of the field, This paper introduces the PQ mode power flow calculationprocedure based on MATLAB software .It is used for a rough calculation of the small and medium power network to achieve its analysis. The algorithm used in the procedure mentioned in this paper is more easy to understand and made some improvements to enhance the computing speed rather than the past. At the same time, the program uses the GUI man-machine dialogue.So Excel table, TXT documents isclosely linked with the MATLAB program to make the input and output interfaces morehumanity.Keywords:power flow calculation；P-Q decomposition mode；MATLAB software目录摘要IAbstract II第1章绪论11.1 课题背景11.2 电力系统潮流计算11.2.1 电力系统潮流计算简介11.2.2 电力系统潮流计算的基本要求21.3 潮流计算的意义及其发展31.4 本次毕业设计主要工作4第2章潮流计算的原理及具体算法过程62.1 电力网络的数学模型62.1.1 电力网络的基本方程62.1.2 导纳矩阵的形成72.1.3 电力网络中几种特殊的数学模型82.2 电力系统潮流计算112.2.1 电力系统潮流计算数学模型112.2.2 电力系统节点分类122.2.3 潮流计算的约束条件132.3 牛顿-拉夫逊法求解潮流计算132.3.1 牛顿-拉夫逊法原理132.3.2 P-Q分解法潮流计算15第3章基于MATLAB软件 P-Q法潮流计算203.1 P-Q分解法程序框图203.2 计算步骤及实现各部分功能的程序213.2.1 原始数据的输入213.2.2 导纳矩阵及，形成233.2.3 计算不平衡功率ΔPi及修正相角Δθi253.2.4 计算不平衡功率ΔQi及修正相电压ΔVi263.2.5 程序运行结果的输出27第4章算例验证与分析284.1 算例说明及分析284.1.1 算例说明284.1.2 算例分析284.2 根据算例输入相应节点线路参数28 4.3 算例运行结果29结论32致谢33参考文献34附录Ａ36附录Ｂ46附录Ｃ63第1章绪论1.1课题背景电力是衡量一个国家经济发展的主要指标，也是反映人民生活水平的重要标志，它已成为现代工农业生产、交通运输以及城乡生活等许多方面不可或缺的能源和动力。

电力系统潮流分析与计算设计(P Q分解法)

电力系统潮流分析与计算设计（P Q分解法）电力系统潮流分析与计算设计（p-q分解法）摘要潮流排序就是研究电力系统的一种最基本和最重要的排序。

最初，电力系统潮流排序就是通过人工手算的，后来为了适应环境电力系统日益发展的须要，使用了交流排序台。

随着电子数字计算机的发生，1956年ward等人基本建设了实际可取的计算机潮流排序程序。

这样，就为日趋繁杂的大规模电力系统提供更多了极其有力的排序手段。

经过几十年的时间，电力系统潮流排序已经发展得十分明朗。

潮流排序就是研究电力系统稳态运转情况的一种排序,就是根据取值的运转条件及系统接线情况确认整个电力系统各个部分的运转状态，例如各母线的电压、各元件中穿过的功率、系统的功率损耗等等。

电力系统潮流排序就是排序系统动态平衡和静态平衡的基础。

在电力系统规划设计和现有电力系统运转方式的研究中，都须要利用电力系统潮流排序去定量的比较供电方案或运转方式的合理性、可靠性和经济性。

电力系统潮流计算分为离线计算和在线计算，离线计算主要用于系统规划设计、安排系统的运行方式，在线计算则用于运行中系统的实时监测和实时控制。

两种计算的原理在本质上是相同的。

实际电力系统的潮流技术主要使用pq水解法。

1974年，由scottb.在文献(@)中首次提出pq分解法，也叫快速解耦法（fastdecoupledloadflow，简写为fdlf）。

本设计就是使用pq水解法排序电力系统潮流的。

关键词：电力系统潮流排序pq水解法第一章概论1.1详述电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它是根据给定的运行条件及系统接线情况确定整个电力系统各个部分的运行状态，如各母线的电压、各元件中流过的功率、系统的功率损耗等等。

电力系统潮流计算是计算系统动态稳定和静态稳定的基础。

在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性。

PQ分解法

4 P-Q 分解法潮流计算 4.1P-Q 分解法的基本原理P-Q 分解法是从简化一极坐标表示的牛顿-拉夫逊法潮流修正方程基础上派生出来的，考虑到了电力系统本身的特点。

牛顿法潮流计算的核心是求解修正方程式。

当节点功率方程式采用极坐标系统时，修正方程式为[∆P ∆Q ]=[H N J L ][∆δ∆U/U ] (4.1)将其展开为{∆P =H∆δ+N(∆U/U)∆Q =J∆δ+L(∆U/U)(4.2) 对修正方程式的第一步简化是：计及电力网络中各元件的电抗远大于电阻，以致各节点电压相位角的改变主要影响各元件中的有功功率及各节点的注入有功功率；各节点电压大小的改变主要影响元件中的无功功率以及各节点的注入无功功率；式（4.2）中子阵N 及J 中各元素的数值相对很小，因此可以略去，从而将式（4.2）简化为 {∆P =H∆δ∆Q =L(∆U/U)（4.3）但是，H 、L 中的元素是电压的函数，在每次迭代中都要重新形成上述H 、L 矩阵，并且又都是不对称矩阵，仍然相当麻烦。

对修正方程式的第二步简化是：由于有对状态变量δi 的约束条件|δi −δj |<|δi −δj |max，即线路两端电压的相角差是不大的，再计及G ij ≪B ij ，可以认为cos δij ≈1 G ij sin δij ≪B ij 于是，H ij 和L ij 的表达式H ij =∂∆P i∂δj=−U i U j (G ij sin δij −B ij cos δij ) i ≠j L ij =∂∆Q i∂U j U j=−U i U j (G ij sin δij −B ij cos δij ) i ≠j 可简化为H ij =U i U j B ij L ij =U i U j B ij (4.4) 再由式H ii =U i 2B ii +Q i (当i =j ，sin δij ≈0，cos δij ≈1时) （4.5） L ii =∂∆Q i ∂U iU i =−U i ∑U j (G ij sin δij −B ij cos δij )+2U i 2B ii =U i 2B ii −Q i j=nj=1j≠i（4.6）按自导纳的定义，上两式中的U i 2B ii 项应为各元件电抗远大于电阻的前提下除节点i 外其他节点都接地时由节点i 注入的无功功率。

P-Q 分解法潮流计算程序设计

z=[0 0.10+0.4*i 0.3*i 0.12+0.5*i;0 0 0 0.08+0.4*i;0 0 0 0;0 0 00];%分别定义所需矩阵的初值，方便编程所需变量%n=length(z);k=zeros(n,n);Y=zeros(n,n);yd=zeros(n,n);y=zeros(n,n);%输入原始数据如下%yd=[0 0.01528*i 0 0.01920*i;0 01528*i 0 0 0.01413*i;0 0 0 0;0.01920*i 0.01413*i 0 0];%其中yd(i,j)表明i,j结点间与i结点接地电阻%%y(i,j)表明i,j结点间正常联接电阻%k(1,3)=1.1;%数据处理%for m=1:nfor j=1:nifz(m,j)~=0y(m,j)=1/z(m,j);y(j,m)=y(m,j);endendfor m=1:nfor j=1:nifk(m,j)~=0y(m,j)=k(m,j)/z(m,j);y(j,m)=y(m,j);yd(m,j)=(k(m,j)-1)*k(m,j)/z(m,j); yd(j,m)=(1-k(m,j))/z(m,j);endendendfor m=1:nfor j=1:nifm==jY(m,j)=sum(y(m,:))+sum(yd(m,:)); elseY(m,j)=-y(m,j);Y(j,m)=Y(m,j);endendA=[-0.3 -0.55 0.5 0;-0.18 -0.13 0 0;1 1 1.1 1.05;0 0 0 0]; G=real(Y);B=imag(Y);%修正方程式的系数矩阵就是导纳矩阵的虚部%B1=B([1,2,3],[1,2,3]);B2=B([1,2,],[1,2,]);for k1=0:100for m=1:(n-1)sum=0;forj=1:nh=A(m,3)*A(j,3)*(G(m,j)*cos(2*pi/360*(A(m,4)-A(j,4)))+B(m,j )*sin(2*pi/360*(A(m,4)-A(j,4))));sum=sum+h;endDP(m,1)=A(m,1)-sum;e n dV1=A([1,2,3],[3]);a=DP./V1;a=inv(-B1)*180/pi*a;DS=V1.\a;A([1,2,3],[4])=A([1,2,3],[4])+DS;for m=1:2sum=0;for j=1:nw=A(m,3)*A(j,3)*(G(m,j)*sin(2*pi/360*(A(m,4)-A(j,4)))-B(m,j )*cos(2*pi/360*(A(m,4)-A(j,4))));sum=sum+w;endDQ(m,1)=A(m,2)-sum;e n dV2=A([1,2],[3]);b=DQ./V2;b=inv(-B2)*b;V2=V2+b;A([1,2],[3])=A([1,2],[3])+b;ifmax(max(abs(DP)),max(abs(DQ)))<0.00001break;endendsum=0;sum1=0;sum2=0;for j=1:nx=A(4,3)*A(j,3)*(G(4,j)*cos(2*pi/360*(A(4,4)-A(j,4)))+B(4,j)* sin(2*pi/360*(A(4,4)-A(j,4))));sum=sum+x;c=A(4,3)*A(j,3)*(G(4,j)*sin(2*pi/360*(A(4,4)-A(j,4)))-B(4,j)* cos(2*pi/360*(A(4,4)-A(j,4))));sum1=sum1+c;d=A(3,3)*A(j,3)*(G(3,j)*sin(2*pi/360*(A(3,4)-A(j,4)))-B(3,j)* cos(2*pi/360*(A(3,4)-A(j,4))));sum2=sum2+d;endA(4,1)=sum;A(4,2)=sum1;A(3,2)=sum2;for i=1:4P(i)=A(i,1);Q(i)=A(i,2);V(i)=A(i,3);S(i)=A(i,4);endY disp('节点1 节点2 节点3 节点4');P disp('节点1 节点2 节点3 节点4');Q disp('节点1 节点2 节点3 节点4'); V disp('节点1 节点2 节点3 节点4'); S disp('节点1 节点2 节点3 节点4'); disp('P Q V S');disp(A);。

PQ分解法潮流计算

否
是
K01=0?
K01=0,t=t+1 否
K01=1
ERM(1)< ε& ERM(,功率误差的数值。EMP，寄存器迭代过程中
最大功率误差。 K01是0时为有功功率，K01是1时为无功功率。
Hij ViVj (Gij sin ij Bij cosij ) Nij ViVj (Gij cosij Bij sin ij ) M ij ViVj (Gij cosij Bij sin ij ) Lij ViVj (Gij sin ij Bij cosij )
化简为
Hij ViVj Bij (i, j 1, 2, , n 1) Lij ViVj Bij (i, j 1, 2, , m)
将上式代入可得到
P H
Q L(V / V )
在实际的P-Q分解法中，两个修正方程的系数矩阵并不相同，一般可以写为
H VBV L VBV 式中:V是由各节点电压幅值组成的对角阵。由于PV节点的存在， B’及B”的阶数不同，分别为n-1阶和m阶。（m<n-1）
由图2-3可以看出，牛顿法在开始时收敛得比较慢，当收敛到一定程度后，它的收敛速度就非常快，而P-Q分解法几乎是按同一速度收敛的。如果给出的收敛条件小于图中A点相应的误差，那么P-Q 分解法所需要的迭代次数要比牛顿法多几次。可以粗略地认为P-Q分解法的选代次数与精度的要求之间存在着线性关系。
P-Q分解法的修正方程式为
P / V B
Q / V BV
通过这一步简化，修正方程式中的系数矩阵B'和 B"由节点导纳矩阵的虚部构成，从而是常数对称矩阵。其区别只是阶数不同，矩阵B'为n -1阶，不含平衡节点对应的行和列，矩阵B"为m阶，不含平衡节点和PV节点所对应的行和列。但在实际 P-Q分解法程序中，为了提高收敛速度，对B'与 B"的构成作了下面一些修改:

P-Q分解法潮流计算

化简为
H ij VV i j Bij Lij VV i j Bij (i, j 1, 2,, n 1) (i, j 1, 2,, m)
将上式代入
可得到
Байду номын сангаас
P H Q L(V / V )
在实际的P-Q分解法中，两个修正方程的系数矩阵并不相同，一般可以写为
H ij VV i j (Gij sin ij Bij cos ij ) Nij VV i j (Gij cos ij Bij sin ij ) M ij VV i j (Gij cos ij Bij sin ij ) Lij VV i j (Gij sin ij Bij cos ij )
(1)一般情况下，线路两端电压的相角差不大(不超过10°~20°)，因此可以认为
cos ij 1,
Gij sin ij Bij
(2)与系统各节点无功功率相对应的导纳通常远小于该节点自导纳的虚部，即
B Li Qi 2 Bii Vi
Qi Vi2Bii
考虑到上述关系，略去相关项可将系数矩阵
P-Q分解法潮流计算
P-Q分解法潮流计算
PQ分解法是由极坐标形式的牛顿法演化而来，以有功功率作为修正电压向量角度的依据，以无功功率作为修正电压幅值的依据，把有功功率和无功功率迭代分开进行。
一、P-Q分解法的基本原理
极坐标形式的牛顿潮流算法的修正方程为
P H Q M
P-Q分解法改变了牛顿法迭代公式的结构，就改变了迭代过程的收敛特性。事实上，依一个不变的系数矩阵进行非线性方程组的迭代求解，在数学上属于“等斜率法”，其选代过程是按几何级数收敛的，若画在对数坐标系上，这种收敛特性基本上接近一条直线。而牛顿法是按平方收敛的，在对数坐标纸上基本上是一条抛物线，如图2-3所示。

第四章 PQ分解法潮流计算

j¹i
极坐标形式的修正方程式（重新排序）
（n1)
é DP1 ù é H11 H12 L H1p H1n M N11 N12 L
ê ê
D P2
ú ú
ê ê
H
21
H 22
L H2p
H2n
M
N21 N22 L
ù é Dd1 ù
úê úê
Dd 2
ú ú
êMú ê
L
M
L úê M ú
ê ê
D Pp
ú ú
ê ê
N22 M H2 p L22 M J 2 p
H2n J2n
ú ú ú
ê ê ê
Dd 2 DU2 /U
2
ú ú ú
êMú ê
L
M
L
úê
ú
ê ê
L
ú ú
=
ê ê
L
L
L
L
M
L
LL
L
L úú
ê ê
L
ú ú
2(nm )
PV节点
ê ê
DPp
ú ú
ê ê
H
p1
N p1
H p2
N p2
M H pp
雅可比矩ê 阵：ú ê
网
自学内容
n 稀疏技术 n 压缩存储 n 节点编号优化 n 高斯消去法 n 因子表法
\
ì í î
Pi Qi
= =
eiai + fiai -
f i bi eibi
极坐标形式的功率方程（电压用相量表示）：
å å Pi - jQi = UiÐ -qi (Gij + jBij )U jÐq j = Ui (Gij + jBij )U j (cosqij - j sinqij )

电力系统分析P-Q分解法潮流计算

专业课程设计报告P-Q分解法潮流计算系别电气工程系专业班级09级电气4班学生XX 钟剑帆学号8指导教师房大中提交日期2011年11月12日目录P-Q分解法潮流计算1一、原理分析1二、程序流程1三、设计内容23.1 程序设计23.2.1 输入数据73.2.2 输出数据83.2.3 结果分析11四、心得总结12一、原理分析从潮流计算的基本方程出发，采用PQ分解法并通过建立矩阵的修正方程来依次迭代，逐步逼近真值来计算电网的电压和功率分布。

二、程序流程三、设计内容3.1 程序设计主函数Sbase_MVA=100.fid=fopen('Nodedata.txt');N=textscan(fid, '%s %u %d %f %f %f %f %f %f')fclose(fid);busnumber=size(N{1},1)for i=1:busnumberBus(i).name=N{1}(i);Bus(i).type=N{2}(i);Bus(i).no=i;Bus(i).Base_KV=N{3}(i);Bus(i).PG=N{4}(i);Bus(i).QG=N{5}(i);Bus(i).PL=N{6}(i);Bus(i).QL=N{7}(i);Bus(i).pb=N{8}(i);Bus(i).V=1.0;Bus(i).angle=0;endfid=fopen('Aclinedata.txt');A=textscan(fid, '%s %s %f %f %f %f')fclose(fid);aclinenumber=size(A{1},1)for i=1:aclinenumberAcline(i).fbname=A{1}(i);Acline(i).tbname=A{2}(i);Acline(i).Base_KV=A{3}(i);Acline(i).R=A{4}(i);Acline(i).X=A{5}(i);Acline(i).hB=A{6}(i);for k=1:busnumberif strcmp(Acline(i).fbname, Bus(k).name)Acline(i).fbno=Bus(k).no;endif strcmp(Acline(i).tbname, Bus(k).name)Acline(i).tbno=Bus(k).no;endendendfid=fopen('Transdata.txt');T=textscan(fid, '%s %f %f %s %f %f %f %f')fclose(fid);tansnumber=size(T{1},1)for i=1:tansnumberTrans(i).fbname=T{1}(i);Trans(i).fbBase_KV=T{2}(i);Trans(i).fbrated_KV=T{3}(i);Trans(i).tbname=T{4}(i);Trans(i).tbBase_KV=T{5}(i);Trans(i).tbrated_KV=T{6}(i);Trans(i).R=T{7}(i);Trans(i).X=T{8}(i);for k=1:busnumberif strcmp(Trans(i).fbname, Bus(k).name)Trans(i).fbno=Bus(k).no;endif strcmp(Trans(i).tbname, Bus(k).name)Trans(i).tbno=Bus(k).no;endendTrans(i).k=Trans(i).tbrated_KV*Trans(i).fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbB ase_KV;tempx=Trans(i).fbrated_KV^2/Trans(i).fbBase_KV^2;Trans(i).X=tempx*Trans(i).X;Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)% for Y=G+ matrix[G,B,B2]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber) ;%B:=B';B2:=B"dlmwrite('Gmatrix.txt', G, 'delimiter', '\t','precision', 6);dlmwrite('Bmatrix.txt', B, 'delimiter', '\t','precision', 6);GBB2pause[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber);JPiJP=-inv(JP)JQiJQ=-inv(JQ)pause%maxiteration=0for i=1:busnumberNodeV(i)=Bus(i).V;Nodea(i)=Bus(i).angle;VX(i)=Bus(i).V*cos(Bus(i).angle);VY(i)=Bus(i).V*sin(Bus(i).angle);dQGQL(i)=Bus(i).QG-Bus(i).QL;dPGPL(i)=Bus(i).PG-Bus(i).PL;endNodeV=NodeV'Nodea=Nodea'%VX=VX'%VY=VY'dQGQL=dQGQL'dPGPL=dPGPL'pause%for nointer=1:10maxdP=1.;maxdQ=1.;epsilon=0.000001;noiteration=0;while (maxdP>epsilon)&(maxdP>epsilon)[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGP L,B,G,busnumber);deltaP;deltaQ;maxdP;maxdQ;da=iJP*deltaP;dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1;if noiteration>20breakendendfor i=1:busnumberBus(i).V=NodeV(i);NodeV(i)=NodeV(i)*Bus(i).Base_KV;Bus(i).angle=Nodea(i);Nodea(i)=Nodea(i)*180/pi;endnoiterationNodea=Nodea'NodeV=NodeV'Clear子函数%生成G、B矩阵function[G,B,X]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber) Y=zeros(busnumber);X=zeros(busnumber);for i=1:busnumberY(i,i)=Y(i,i)+Bus(i).pb*j;endfor i=1:aclinenumberf=Acline(i).fbno;t=Acline(i).tbno;Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(t,t)=Y(t,t)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(f,t)=Y(f,t)-1/(Acline(i).R+Acline(i).X*j);Y(t,f)=Y(t,f)-1/(Acline(i).R+Acline(i).X*j);X(f,f)=X(f,f)-1/Acline(i).X;X(t,t)=X(t,t)-1/Acline(i).X;X(f,t)=1/Acline(i).X;X(t,f)=1/Acline(i).X;endfor i=1:tansnumberf=Trans(i).fbno;t=Trans(i).tbno;Y(f,f)=Y(f,f)+1/(Trans(i).R+Trans(i).X*j);Y(t,t)=Y(t,t)+1/(Trans(i).R+Trans(i).X*j)/Trans(i).k^2;Y(f,t)=Y(f,t)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;Y(t,f)=Y(t,f)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;X(f,f)=X(f,f)-1/Trans(i).X;X(t,t)=X(t,t)-1/Trans(i).X;X(f,t)=1/Trans(i).X;X(t,f)=1/Trans(i).X;endG=real(Y);B=imag(Y);end%生成JP、JQ矩阵function [JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber)JP=B;JQ=B2;for i=1:busnumberif Bus(i).type==1for k=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;JP(i,k)=0.;JP(k,i)=0.;endJQ(i,i)=1.;JP(i,i)=1.;endif Bus(i).type==3for k=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;endJQ(i,i)=1.;endendend%计算偏节点PQ差量function[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G ,busnumber)deltaQ=dQGQL;deltaP=dPGPL;maxdP=0.;maxdQ=0.;for i=1:busnumberif Bus(i).type==1deltaQ(i)=0.;deltaP(i)=0.;endif Bus(i).type==3deltaQ(i)=0.;%y1=0;%y2=0;y3=0;for k=1:busnumberif (B(i,k)~=0|G(i,k)~=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k)));endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;endif Bus(i).type==2%y1=0;%y2=0;y3=0;y4=0;for k=1:busnumberif (B(i,k)~=0|G(i,k)~=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k))); y4=y4+NodeV(k)*(G(i,k)*sin(Nodea(i)-Nodea(k))-B(i,k)*cos(Nodea(i)-Nodea(k)));endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;deltaQ(i)=deltaQ(i)-y4*NodeV(i);%deltaQ2(i)=(deltaQ2(i)-(y1*VY(i)-y2*VX(i)))/Bus(i).V;endif maxdP<abs(deltaP(i));maxdP=abs(deltaP(i));endif maxdQ<abs(deltaQ(i));maxdQ=abs(deltaQ(i));enddeltaP(i)=deltaP(i)/NodeV(i);deltaQ(i)=deltaQ(i)/NodeV(i);endend3.2.1 输入数据节点数据（Nodedata.txt）bus1 1 18. 0. 0. 0. 0. 0.bus2 3 18. 1.63 0. 0. 0. 0.bus3 3 18. 0.85 0. 0. 0. 0.bus4 2 230. 0. 0. 0. 0. 0.bus5 2 230. 0. 0. 1.25 0.5 0.bus6 2 230. 0. 0. 0.9 0.3 0.bus7 2 230. 0. 0. 0. 0. 0.bus8 2 230. 0. 0. 1.0 0.35 0.bus9 2 230. 0. 0. 0. 0. 0.支路数据（Aclinedata.txt）bus4 bus5 230. 0.01 0.085 0.088bus4 bus6 230. 0.017 0.092 0.079bus5 bus7 230. 0.032 0.161 0.153bus6 bus9 230. 0.039 0.17 0.179bus7 bus8 230. 0.0085 0.072 0.0745bus8 bus9 230. 0.0119 0.1008 0.1045变压器数据（Transdata.txt）bus1 18.0 18.0 bus4 230. 230. 0.0 0.0576bus2 18.0 18.0 bus7 230. 230. 0.0 0.0625bus3 18.0 18.0 bus9 230. 230. 0.0 0.05863.2.2 输出数据Sbase_MVA =100N ={9x1 cell} [9x1 uint32] [9x1 int32] [9x1 double] [9x1 double] [9x1 double] [9x1 double] [9x1 double] [8x1 double]busnumber =9A ={6x1 cell} {6x1 cell} [6x1 double] [6x1 double] [6x1 double] [6x1 double]aclinenumber =6T ={3x1 cell} [3x1 double] [3x1 double] {3x1 cell} [3x1 double] [3x1 double] [3x1 double] [3x1 double]tansnumber =3NodeV =111111111 Nodea =0 dQGQL =-0.5000-0.3000-0.3500dPGPL =1.63000.8500-1.2500-0.9000-1.0000noiteration =9Nodea =0 9.6687 4.7711 -2.4066 -4.3499 -4.0173 3.7991 0.6215 1.9256NodeV =18 18 18 227 220 224 229 227 2313.2.3 结果分析busnumber=9，aclinenumber=6，tansnumber=3表示此次计算为9节点，6支路，3变压器的网络。

基于MATLAB软件的P_Q分解法潮流计算

在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性，可靠性和经济性。

所以，电力系统潮流计算是进行电力系统故障计算，继电保护整定，安全分析的必要工具。

随着电力系统网络的急剧扩大和不断复杂，运用手算进行潮流计算已经不现实。

但是，伴随着计算机技术的飞速发展，基于计算机的潮流计算也就应运而生了。

这样，通过潮流计算，实现对系统的分析成为可能。

本文所设计的程序，在计算中，所用的算法通俗易懂并对以往的主流算法做了一些改进，提高了计算速度。

同时，该程序采用了GUI人机对话，将Excel表格、TXT文档与MATLAB程序紧密联系起来，使输入输出界面更加人性化。

关键词：电力系统潮流计算；P-Q分解法；MATLAB软件Power flow calculation of P-Q mode basedon MATLAB softwareAbstractPower flow calculation is one of the important calculations which are to study the operation of power system steady state analysis. It is based on the given operating conditions and system wiring to identify the various parts of the power system operating state: the buses voltage, the stream components power, system power loss. both power system planning design and operation of existing power system mode of study are need to use the power flow calculation to quantitatively compare the program or run mode power supply reasonable, reliability and economy. Therefore, the power flow calculation is an essential tool for a calculation of power system faults, protection setting, security analysis. with the rapid expansion of power system network and continuing to be more complex, using hand calculation for flow calculation has been unrealistic. But ,with the celerity development in computer technology, computer-based power flow calculation has also emerged. In this way, It is possible to analysis power system through the power flow calculation.Based on the three basic requirements of power flow calculation and followed by the development of the field, This paper introduces the PQ mode power flow calculation procedure based on MATLAB software .It is used for a rough calculation of the small and medium power network to achieve its analysis. The algorithm used in the procedure mentioned in this paper is more easy to understand and made some improvements to enhance the computingspeed rather than the past. At the same time, the program uses the GUI man-machine dialogue. So Excel table, TXT documents is closely linked with the MATLAB program to make the input and output interfaces more humanity.Keywords:power flow calculation；P-Q decomposition mode；MA TLAB software目录摘要 (I)Abstract (II)第1章绪论 (6)1.1课题背景 (6)1.2电力系统潮流计算 (6)1.2.1电力系统潮流计算简介 (6)1.2.2电力系统潮流计算的基本要求 (7)1.3潮流计算的意义及其发展 (8)1.4本次毕业设计主要工作 (9)第2章潮流计算的原理及具体算法过程 (11)2.1电力网络的数学模型 (11)2.1.1电力网络的基本方程 (11)2.1.2导纳矩阵的形成 (12)2.1.3电力网络中几种特殊的数学模型 (13)2.2电力系统潮流计算 (16)2.2.1电力系统潮流计算数学模型 (16)2.2.2电力系统节点分类 (16)2.2.3潮流计算的约束条件 (17)2.3牛顿-拉夫逊法求解潮流计算 (18)2.3.1牛顿-拉夫逊法原理 (18)2.3.2 P-Q分解法潮流计算 (19)第3章基于MATLAB软件 P-Q法潮流计算 (24)3.1 P-Q分解法程序框图 (24)3.2计算步骤及实现各部分功能的程序 (25)3.2.1原始数据的输入 (25)3.2.2导纳矩阵及B'，B''形成 (27)3.2.3计算不平衡功率ΔP i及修正相角Δθi (29)3.2.4计算不平衡功率ΔQ i及修正相电压ΔV i (30)3.2.5程序运行结果的输出 (31)第4章算例验证与分析 (32)4.1算例说明及分析 (32)4.1.1算例说明 (32)4.1.2算例分析 (32)4.2根据算例输入相应节点线路参数 (32)4.3算例运行结果 (33)结论 (36)致 (37)参考文献 (38)附录Ａ (40)附录Ｂ (50)附录Ｃ (67)第1章绪论1.1课题背景电力是衡量一个国家经济发展的主要指标，也是反映人民生活水平的重要标志，它已成为现代工农业生产、交通运输以及城乡生活等许多方面不可或缺的能源和动力。

电力系统分析P-Q分解法潮流计算终稿

《电力系统分析》课程设计报告题目：复杂电力系统潮流分析程序设计所在学院电气工程学院专业班级12级电气工程及其自动化一班学生姓名陈剑秋学生学号************指导教师房大中提交日期 2014年 12 月 16日摘要电力系统潮流计算是电力系统分析课程基本计算的核心部分之一。

它既有本身的独立意义，又是电力系统规划设计，运行分析和理论研究的基础。

电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算，是电力系统计算分析中的一种最基本的计算。

潮流计算是电力系统的各种计算的基础，同时它又是研究电力系统的一项重要分析功能，是进行故障计算，继电保护鉴定，安全分析的工具。

电力系统潮流计算是计算系统动态稳定和静态稳定的基础。

在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性本课程设计采用P-Q分解法对三机九节点系统进行潮流计算。

计算结果的得出是通过MATLAB软件编程实现。

关键词：电力系统潮流计算，P-Q分解法，三机九节点系统AbstractPower flow calculation is the core part of the power system analysis course basic computing. It has both the independent meaning itself, and is the power system planning and design, the basic operation analysis and theoretical study.Power flow calculation is the calculation of the steady state operation of normal and fault conditions of complex power system, is one of the most basic calculation of power system calculation and analysis of. Power flow calculation is the basis of all kinds of power system calculation, and it is also an important function of power system analysis study, is the fault calculation, relay protection and identification, safety analysis tool. Power flow calculation is the basis of calculation system of dynamic and static stabilization. In the research of operation mode of power system planning and the existing power system, the rationality of power supply, reliability and economy of the program or run mode comparison requires the use of power flow calculation to quantitativeThe curriculum design using P-Q decomposition method for power flow calculation of the three machine nine bus system. It is concluded that the calculation result is achieved by MATLAB software programming.Keywords:power system power flow calculation, the P-Q decomposition method, three machine nine node system目录摘要 (I)Abstract (II)一、绪论 (1)1.1本课题的目的和意义 (1)1.2国内外发展现状 (1)1.2.1高斯-赛德尔迭代法 (2)1.2.2牛顿-拉夫逊法和P-Q分解法 (2)1.2.3基于MATLAB 的电力系统潮流计算发展前景 (3)二、设计目的 (4)三、设计要求和设计指标 (4)四、设计内容 (4)4.1 选题内容 (4)4.2 基础资料 (5)4.3算法原理 (6)4.3.1 节点导纳矩阵形成的计算机方法 (6)4.3.2 电力网络的潮流方程 (9)4.3.3 电力网络极坐标形式的潮流方程 (10)4.3.4 极坐标牛顿潮流算法的修正方程雅和可比矩阵 (10)4.3.5 PQ分解法潮流计算 (11)五、程序设计 (13)5.1 主函数 (13)5.2 子函数 (16)六、输入与输出结果 (21)6.1 输入数据 (21)6.2 输出结果 (22)6.3 当bus4-6的发生故障被切断后，系统的运行情况 (25)6.4 两机五节点网络 (30)七、结果分析 (34)7.1 三机九节点系统正常运行结果分析 (34)7.2 三机九节点系统故障切除后结果分析 (35)结论 (35)参考文献 (36)一、绪论1.1本课题的目的和意义电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。

新版本P-Q分解法

j 1
n
U i (Gij jBij )U j (cos i cos j sin i sin j j cos i sin j j cos j sin i )
j 1 n
n
U i (Gij jBij )U j (cos ij j sin ij )
对f(x)按泰勒级数展开并略去高阶项
(0) ( 0) f1 ( x1( 0) , x2 ,..., xn )
f1 f f x1 1 x2 ... 1 xn y1 x1 0 x2 0 xn 0 f1 f f x1 1 x2 ... 1 xn x1 0 x2 0 xn 0
n个节点的网络，有（n-m-1）个PQ节点，有 m个PV节点（1个平衡节点）。
P1 H11 Q J 1 11 P2 H 21 Q 2 = J 21 P p H p1 Pn H n1
标幺值的电压等级归算：将网络中的有名值归算到同一电压等级— —基本级，然后除以基本级的基准值；将基本级的基准值归算到各电压等级，然后用各电压下的有名值除以该电压等级下的基准值。

四、潮流计算算法
1、潮流方程节点导纳法，节点电压与节点电流之间的关系为：
. . I YU
(1)
展开得：
三、标幺制
标幺制计算结果清晰、便于迅速判断计算结果的正确性、可大量简化计算。
有名值（、S、kV、kA、MVA等）标幺值基准值（与相应有名值单位相同）
基准值的选择：

阻抗、导纳的基准值为每相阻抗、导纳；电压、电流的基准值为线电压、线电流；功率基准值为三相功率。

PQ分解法潮流计算实验

xxxx实验报告学生姓名：学号：专业班级：实验类型：□验证□综合■设计□创新实验日期： 2010.10.16 实验成绩：一、实验项目名称P-Q分解法潮流计算实验二、实验目的与要求：目的：电力系统分析的潮流计算是电力系统分析的一个重要的部分。

通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。

电力系统潮流的计算和分析是电力系统运行和规划工作的基础。

潮流计算是指对电力系统正常运行状况的分析和计算。

通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、P-Q分解法、直流潮流法，以及由高斯-塞德尔法、牛顿-拉夫逊法演变的各种潮流计算方法。

本实验采用P-Q分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。

通过实验教学加深学生对电力系统潮流计算原理的理解和计算，初步学会运用计算机知识解决电力系统的问题，掌握潮流计算的过程及其特点。

熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。

要求：编制调试电力系统潮流计算的计算机程序。

程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。

三、主要仪器设备及耗材每组计算机1台、相关计算软件1套四、实验内容：1.理论分析：P-Q分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。

牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，修正方程式为：⎥⎦⎤⎢⎣⎡∆∆⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∆∆V V L N J H Q P /δ 或展开为：VV L J Q VV N H P //∆⋅+∆⋅=∆∆⋅+∆⋅=∆δδ (4)以上方程式是从数学上推倒出来的，并没有考虑电力系统这个具体对象的特点。

PQ分解潮流算法简介

P( k ) H ( k )Q(k)M(k
)
N ( k ) θ( k )
L( k )
U
(
k
)
U
(
k
)
P( k ) P(θ( k ) ,U( k ) )
Q( k ) Q(θ( k ) ,U( k ) )
θ( k ) 1( k ) 2( k )
( k ) n1
T
U( k )
PQ分解潮流算法简介
前言
潮流计算的内容:
根据给定的电网结构、发电计划及负荷分布情况，求出整个电网的运行状态。（运行状态：节点母线的电压、相角。再由状态变量计算线路输送的有功和无功功率。）
潮流计算的意义:
（1）潮流计算，对于系统运行方式的分析，对电网规划阶段中设计方案的确定都是必不可少的。为判别这些运行方式及规划设计方案的合理性、安全性、可靠性及经济性提供了定量分析的依据。
Ui( k
)U
( j
k
)
(
Gij
sin
( ij
k
)
Bij
cos
( ij
k
)
)
Ui(
k
)U
( j
k
)
Bij
H ( k ) ii
Pi
i
U U( k )
[U
( i
k
)
]2
Bii
U
( i
k
)
U
( j
k
)
(
Gij
ji
sin
( ij
k
)
Bij
cos
( ij
k
)
)
θ θ( k )