电力系统稳态分析摘要电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗。
所以，电力系统潮流计算是进行电力系统故障计算，继电保护整定，安全分析的必要工具。
本文介绍了基于MATLAB软件的牛顿-拉夫逊法和P-Q分解法潮流计算的程序，该程序用于计算中小型电力网络的潮流。
在本文中，采用的是一个5节点的算例进行分析，并对仿真结果进行比较，算例的结果验证了程序的正确性和迭代法的有效性。
关键词：电力系统潮流计算；MATLAB；牛顿-拉夫逊法；P-Q分解法；目次1 绪论 (1)1.1背景及意义 (1)1.2相关理论 (1)1.3本文的主要工作 (2)2 潮流计算的基本理论 (3)2.1节点的分类 (3)2.2基本功率方程式（极坐标下） (3)2.3本章小结 (4)3 潮流计算的两种算法 (5)3.1牛顿—拉夫逊算法 (5)3.2PQ分解算法 (10)3.3本章小结 (15)4 算例 (16)4.1系统模型 (16)4.2结果分析 (16)4.3本章小结 (19)结论 (20)参考文献 (21)附录 (22)1 绪论1.1背景及意义电力系统稳态分析是研究电力系统运行和规划方案最重要和最基本的手段。
电力系统稳态分析根据给定的发电运行方式和系统接线方式来确定系统的稳态运行状态，其中潮流计算针对电力系统的各种正常的运行方式进行稳态分析。
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。
通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。
电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都离不开迭代。
潮流计算方法的改进过程中，经历了高斯-赛德尔迭代法、阻抗法、分块阻抗法、牛顿-拉夫逊法、改进牛顿法、P-Q分解法等。
现在比较常用的方法就是牛顿-拉夫逊法和P-Q分解法。
对潮流计算的要求可以归纳为下面几点：（1）计算方法的可靠性或收敛性；（2）对计算机内存量的要求；（3）计算速度；（4）计算的方便性和灵活性。
1.2 相关理论所谓潮流计算，就是已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、各支路电流与功率及网损。
对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。
对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。
潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。
在运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以校验运行可靠性。
在电力系统调度运行的多个领域都涉及到电网潮流计算。
潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。
1.3 本文的主要工作本文介绍了电力系统潮流计算方法中的牛顿-拉夫逊法和PQ快速分解法的相关知识及其基本原理，并用MATLAB编写程序，最后通过一个5节点的算例来验证该程序的正确性，并对两种算法的结果进行了分析，对比了两种算法。
2 潮流计算的基本理论2.1 节点的分类用一般的电路理论求解网络方程，目的是给出电压源（或电流源）研究网络内的电流（或电压）分布，作为基础的方程式，一般用线性代数方程式表示。
然而在电力系统中，给出发电机或负荷连接母线上电压或电流（都是向量）的情况是很少的，一般是给出发电机母线上发电机的有功功率和母线电压的幅值，给出负荷母线上负荷消耗的有功功率和无功功率。
主要目的是由这些已知量去求电力系统内的各种电气量。
所以，根据电力系统中各节点性质的不同，很自然的把节点分成三种类型。
（1）PQ节点对这一类节点，事先给定的是节点有功功率和无功功率（P、Q），待求的未知量是节点电压向量(V， )，所以叫“PQ节点”。
通常变电所母线都是PQ 节点，当某些发电机的输出功率P、Q给定时，也作为PQ节点。
在潮流计算中，系统大部分节点属于PQ节点。
(2)PV节点这类节点给出的参数是该节点的有功功率P及电压幅值V，待求量为该节点的无功功率Q及电压向量的相角θ。
这类节点在运行中往往要有一定可调节的无功电源，用于维持给定的电压值。
通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线作PV节点处理。
(3)平衡节点在潮流计算中，这类节点一般只设一个。
对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。
也就是说，对平衡节点给定的运行参数是V和θ，因此又称为Vθ节点，而待求量是该节点的P，Q，整个系统的功率平衡由这一节点承担。
关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂（或发电机），有时也可能按其他原则选择，例如，为提高计算的收敛性，可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。
以上三类节点4个运行参数P、Q、V、θ中，已知量都是两个，待求量也是两个，只是类型不同而已。
2.2 基本功率方程式（极坐标下）在潮流计算中任何复杂的电力系统都可以归结为以下元件（参数）组成：发电机（注入电流或功率）；负荷（注入负的电流或功率）；输电线支路（电阻、电抗）；变压器支路（电阻、电抗、变比）；母线上对地支路（阻抗和导纳）；线路上的对地支路（一般为线路充电电容导纳）。
必须指出，如果仅研究稳态情况下的潮流而不涉及暂态过程的计算，则不需要发电机和负荷的阻抗参数，只需要给出发电机和负荷的注入功率和电流，并且规定发电机和负荷的注入功率和电流取正，而负荷取负。
在潮流计算中，节点功率可表示为：^^j i j ij i i i V Y V jQ P ∑∈•=+ （i=1,2,...n) （2.1）若把电压表示为极坐标的形式，即i j i i e V V θ=•（2.2）将导纳矩阵中元素表示为ij ij ij jB G Y += （2.3）这样，我们可以得到：∑∈+-=+ij ij ij ij ij j i i i j jB G V V jQ P )sin )(cos (θθ (i=1,2,...n) （2.4）按实部和虚部展开，得到(cos sin )(sin cos )i i j ij ij ij ij j iii j ij ij ij ij j i P V V G B Q V V G B θθθθ∈∈⎧=+⎪⎨=-⎪⎩∑∑ （2.5） 上式就是功率的极坐标方程式。
该方程组在牛顿法和P-Q 分解法中起到了重要作用。
2.3 本章小结在本章里主要介绍了电力系统的潮流计算的基本的理论。
首先对电力系统中三种节点进行了详细的阐述；其次，介绍了在极坐标的情况下的电力系统的潮流计算的功率方程，为下文的潮流计算分析打下基础。
3 潮流计算的两种算法 3.1 牛顿—拉夫逊算法牛顿-拉夫逊算法产生于50 年代末期，是一种实用且有竞争力的电力系统潮流计算方法，求解非线性方程式的典型方法。
在稀疏矩阵技巧和高斯消去法被应用以后，其真正的价值才体现出来。
该方法有较好的收敛性，迭代次数少，在电力系统潮流计算中也得到应用。
目前，牛顿法潮流计算是最为广泛、效果最好的一种潮流计算方法。
3.1.1 基本原理该方法把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程，通常称为逐次线性化过程，经过一步步的迭代，得到最终的结果。
设有非线性方程式：0)(=x f （3.1） 设(0)x 为该方程式的初值。
(0)(0)x x x =-∆ （3.2）(0)x ∆为初值(0)x 的修正量。
将式（3.2）代入式（3.1），可以得到(0)(0)()0f x x -∆= （3.3）按泰勒级数展开。
因为(0)x ∆很小，(0)x ∆二次以及二次以上的各项均可以略去，可以得到简化的方程式：(0)'(0)(0)()()0f x f x x -∆= （3.4）上式是对于变量修正量(0)x ∆的线性方程式，即修正方程式，解得：(0)(0)'(0)()()f x xf x ∆= （3.5） 上式求得的)0(x ∆不是方程真正的解，我们需要进行反复的迭代，一步步的趋近方程式的解，得到最逼近的值。
这样反复下去，就构成了不断求解非线性方程式的逐次线性化过程。
第t 次迭代时的参数方程为0)(')()()()(=∆-t t t x x f x f （3.6）当0)()(→t x f 时，就满足了原方程式（3.1），因而)(t x 就成为该方程的解。
式中)(')(t x f 是函数0)(=x f 在)(t x 点的一次导数，也就是曲线在)(t x 点的斜率，如图3.1所示，()'()()t t tg f x α= （3.7）修正量)(t x ∆则是由)(t x 点的切线与横轴的交点来确定，由图3.1可以直观的看出牛顿法的求解过程。
yxx(1)t x +∆(1)t x +()t x∆()t α(1)()t f x +()y f x =()()t f x o图3.1 牛顿法的几何解释现在把牛顿法推广到多变量非线性方程组的情况。
设有变量12,,,n x x x L 的非线性联立方程组：11221212(,,,)0(,,,)0(,,,)0n n n n f x x x f x x x f x x x =⎫⎪=⎪⎬⎪⎪=⎭L L M L (3.8)给定各变量初值)0()0(2)0(1,,,n x x x Λ，假设)0()0(2)0(1,,,n x x x ∆∆∆Λ为其修正量，并使其满足⎪⎪⎭⎪⎪⎬⎫=∆-∆-∆-=∆-∆-∆-=∆-∆-∆-0),,,(0),,,(0),,,()0()0()0(2)0(2)0(1)0(1)0()0()0(2)0(2)0(1)0(12)0()0()0(2)0(2)0(1)0(11n n n n n n nx x x x x x f x x x x x x f x x x x x x f ΛMΛΛ （3.9）对以上n 个方程式分别按泰勒级数展开，当忽略)0()0(2)0(1,,,n x x x ∆∆∆Λ所组成的二次项和高次项时，可以得到⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎬⎫=⎥⎥⎦⎤⎢⎢⎣⎡∆∂∂+∆∂∂+∆∂∂-=⎥⎥⎦⎤⎢⎢⎣⎡∆∂∂+∆∂∂+∆∂∂-=⎥⎥⎦⎤⎢⎢⎣⎡∆∂∂+∆∂∂+∆∂∂-0),,,(0),,,(0),,,()0(0)0(202)0(101)0()0(2)0(1)0(02)0(2022)0(1012)0()0(2)0(12)0(01)0(2021)0(1011)0()0(2)0(11n n n n n n n n n n n n nx x f x x f x x f x x x f x x f x x f x x f x x x f x x f x x f x x f x x x f ΛΛM ΛΛΛΛ （3.10） 式中：ii x f ∂∂为函数),,,(21n i x x x f Λ对自变量j x 的偏导数在点（)0()0(2)0(1,,,n x x x ∆∆∆Λ）处的值。