利用PSASP进行潮流计算

目录1 设计目的 (1)2 关于 PSASP (1)2.1 软件简介 (1)2.2 PSASP软件的体系结构 (1)3 关于牛顿—拉夫逊算法 (2)3.1 牛顿—拉夫逊算法简介 (2)3.2牛顿—拉夫逊法计算潮流 (3)4 九节点电力系统的单线图及元件数据 (4)4.1单线图 (4)4.2 元件数据 (5)5 潮流计算结果 (8)6 结论 (11)7 参考文献 (12)1 设计目的电力系统潮流计算就是对复杂电力系统正常和故障条件下稳态运行状态的计算。

潮流计算的目标是求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统个元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗。

通过电力系统分析仿真软件PSASP7.0对任务书中所给出的9节点系统进行潮流计算，并导出结果。

2 关于PSASP2.1 软件简介PSASP是电力系统分析综合程序（Power System Analysis Software Package的简称。

PSASP 是一套由电科院开发的具有我国自主知识产权，便捷高效、高度集成的开放软件。

它基于电网基础数据库、固定模型库以及用户自定义模型库的支持，可以进行电力系统的各种分析计算，例如：潮流计算、短路电流计算、网损分析、静态安全分析等。

并可以输出各种计算结果。

2.2 PSASP软件的体系结构PSASP的体系分为三层，第一层为公用数据和模型的资源库；第二层为基于资源库的应用程序包；第三层为计算结果库和分析工具。

在使用PSASP时，用户首先利用电网基础数据库、模型库、用户程序库讲模型数据输入到PSASP中；然后使用应用程序包对输入的模型数据进行潮流计算、网损分析等计算；最后将计算结果输出到结果库并可以将结果用报表、图形、曲线等形式输出出来。

PSASP软件的特点1、有公用数据库作支持，可以共用基础数据。

2、有固定模型库和用户自定义模型库作支持。

3、有文本和图形两种方式计算。

4、有多种形式的结果分析输出。

5、有多种常用软件接口。

6、有能力计算大规模的交直流混合电力系统。

3 关于牛顿—拉夫逊算法3.1 牛顿—拉夫逊算法简介牛顿—拉夫逊法（Newton-Raphson 法）是求解非线性代数方程的计算方法中一种有效的计算方法。

在使用这个此算法进行迭代的过程中，非线性问题通线性化而得到逐步化简。

牛顿—拉夫逊法早在20世纪50年代末就已经被应用于求解电力系统潮流问题。

设有非线性方程组()()()⎪⎪⎪⎭⎪⎪⎪⎬⎫===y x x x fy x x x f y x x x f n n n n n ,,,,,,,,,2122121211 （3-1）设其的初值()()()x x xn 00201,,, ，它与真解的误差为x 1∆，x 2∆，…，x n ∆，则式3-1可写成()()()⎪⎪⎪⎭⎪⎪⎪⎬⎫=∆+∆+∆+=∆+∆+∆+=∆+∆+∆+y x x x x x x f y x x x x x x f y x x x x x x f n n n n n n n n )0(2)0(21)0(12)0(2)0(21)0(121)0(2)0(21)0(11,,,,,,,,,（3-2） 再将上式2-1用泰勒级数展开可得()()()()()()()()()()()()⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎬⎫=∆∂∂++∆∂∂+∆∂∂+=∆∂∂++∆∂∂+∆∂∂+=∆∂∂++∆∂∂+∆∂∂+n n n n n n n n n n n y x x f x x f x x f x x x f y x x f x x f x x f x x x f y x x f x x f x x f x x x f 1221011100,,20,10210221011100,,20,10211221011100,,20,101.........................（3-3）取第一式为例()()y x x f x x f x x f x x x f x x x x x x f n n n n n 112010)0()0(2)0(11)0(2)0(21)0(1112111,,,,,,=+∆++∆+∆+=∆+∆+∆+∂∂∂∂∂∂φ （3-4）式中：x f 110∂∂，x f 210∂∂，……，x f n ∂∂10分别表示以()()()x x x n 00201,,, 代入这些偏导数表示式时计算所得，φ1则是一包含x 1∆，x 2∆，…x n ∆的高次方与f 1的高阶偏导数乘积的函数。

如近似解x i )0(与精确解相差不大，则x i ∆的高次方可略去，从而φ1也可略去，由此可以得到一组线性方程组，常称为修正方程组。

它可以用矩阵的形式表示()()()⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡---∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂x x x x f x f x f x f x f x f x f x f x f x x x f y x x x f y x x x f y n n n n n n n n n n n n 21000000000)0()0(2)0(1)0()0(2)0(122)0()0(2)0(111212*********,,,,,,,,, （3-5）或者简写为 x J f ∆=∆ （3-6）式中：J 称为函数f i 的雅可比矩阵；x ∆为由x i ∆组成的列向量；f ∆则称不平衡量的列向量。

将x i )0(代入，可得f ∆、J 中的各元素。

然后用任何一种解线性代数方程的方法，可求得x i )0(∆，从而求得经第一次迭代后x i 的新值x x x i i i )0()0(∆+=。

再将求得的x i )1(代入，又可求得f ∆、J 中各元素的新值，从而解得x i )1(∆以及x x x i i i )1()1()2(∆+=。

如此循环不已，最后可获得对初始式子足够精确的解。

运用这种方法计算时，初值要选择比较接近他们的精确解，否则迭代过程可能不收敛3.2牛顿—拉夫逊法计算潮流形成了雅可比矩阵后并建立了修正方程式，运用牛顿—拉夫逊法计算潮流的核心问题已解决，已有可能列出基本计算步骤并编制流程图。

显然，其修正方程式有两中不同表示方式，但牛顿—拉夫逊法潮流计算的基本步骤却大体上一致，如下几步：1、形成节点导纳矩阵Y B ；2、设各节点电压的初值e i )0(，f i)0(；3、将各节点电压初始值代入式（2—9）中，求修正方程式中的不平衡量P i )0(∆,Q i)0(∆和U i)0(2∆； 4、将各节点电压的初值代入雅可比矩阵表达式中，求修正方程式的系数矩阵； 5、解修正方程式，求各节点电压的变化量，即修正量e i )0(∆、f i)0(∆；6、计算各节点电压的新值，即修正后的值e e e i ii )0()0(∆+= ； f f f iii )0()0(∆+=； (3-7)7、运用各节点电压新值顺次进行下一次迭代； 8、计算平衡点功率和线路功率。

其中平衡点功率为Q P U Y U S s s i n i is s s j +==*=*•∑1~ （3-8）线路功率为Q P I U S ij ij ij i ij j +==*•~ （3-9）Q P I U S ji ji ji j ji j +==*•~ （3-10）从而，线路上损耗的功率为Q P S S S ij ij ji ij ij j ∆+∆=+=∆~~~（3-11）4 九节点电力系统的单线图及元件数据4.1单线图由PSASP 软件可导出单线图，如图4-1所示，4.2 元件数据九节点电力系统中各元件数据如下面各表所示表4-1两绕组变压器数据报表数据组变压器的名称有效标记I侧节点名J侧节点名编号I侧连接方式J侧连接方式单位X1 Gm Bm Tk图4-15 潮流计算结果经潮流计算后以表格形式输出各表如下所示6 结论电力系统潮流计算就是对复杂电力系统正常和故障条件下稳态运行状态的计算。

潮流计算的目标是求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统个元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗。

该设计使用计算机软件PSASP并利用牛顿—拉夫逊算法对给出的9节点电力系统进行了潮流计算并导出了结果。

精确的潮流计算可以节约能源使能源充分的发挥作用，并使电力网处于稳定且平衡的状态，对电网的安全运行起到了重要作用。

在经过了这个学期中对《电力系统分析》的学习，我对电力系统各方面有了新的了解，也有了一定的能力去理解现实生活中复杂的电力系统是如何建立、如何运行的。

理论需要联系实践，这次课程设计给了我良好的契机去实践我所掌握的知识，并让我熟悉并掌握了使用PSASP软件进行电力系统稳态的潮流计算的方法。

7 参考文献[1]孟祥萍,高嬿.电力系统分析（第2版）[M].北京:高等教育出版社,2010.12[2]祝淑萍等.电力系统分析课程设计与综合实验[M].北京:中国电力出版社,2007[3]华智明,张瑞林.电力系统[M].重庆大学出版社,1997[4]徐政电力系统分析学习指导[M].北京机械工业出版社,2002。