电力系统课程设计潮流计算

1、1 潮流计算概述 在电力系统运行和规划中，都需要研究电力系统稳定运行情况，确定电力系统的稳态运行状态。给定电力系统的网络结构、参数和决定电力系统运行状况的边界条件，确定电力系统运行的方法之一是朝流计算。从数学上说,朝流计算是要求解一组有潮流方程描述的非线性方程组。电力系统潮流计算是电力系统分析中最重要最基本的计算，是电力运行、规划以及安全性、可靠性分析和优化的基础，也是各种电磁暂态和机电暂态分析的基础和出发点。 潮流计算方法的发展是与人们所使用的计算工具的发展相联系的。早期，除了手动潮流外，人们用交流计算台通过物理模拟的方法来分析电力系统稳态运行状态。这种方法虽然直观，物理概念清楚，但受到系统

2、规模等因素的影响，分析大网络的超流会遇到困难。 作为研究电力系统稳定运行情况的一种基本电气计算 ，电力系统常规潮流计算的任务是根据给定的网络结构及运行条件，求出整个网络的运行状态，其中包括各母线的电压、网络中的功率分布以及功率损耗等。 潮流计算的结果，无论是对于现有系统运行方式的分析研究，还是对规划中供电方案的分析比较，都是必不可少的。它为判别这些运行方式计规划设计方案的合理性、安全可靠性及经济性提供了定量分析的依据。 此外，在进行电力系统静态及暂态稳定计算时，要利用潮流计算的结果作为其计算的基础；一些故障分析以及优化计算也需要有相应的潮流计算作配合；潮流计算往往成为上述计算程序的一个重要组成

3、部分。以上这些，主要是在系统规划设计及运行方式安排中的应用，属于离线计算范畴。 随着现代化的调度中心的建立，为了对电力系统进行实时安全监控，需要根据实时数据库所提供的信息。2 题目要求 在图1所示的简单电力系统中，系统中节点1、2、3为节点，节点4为节点，节点5为平衡节点，已给定，网络各元件参数的标幺值如表1所示，给定电压的初始值如表2所示，收敛系数。试求： 采用极坐标下的分解法计算图1网络的潮流分布。图1 简单电力系统表2.1 网络各元件参数的标幺值支路电阻电抗输电线路变压器变比k110.04120.1450.5810.021150.0820.427220.05230.1040.5180.0

4、18340.0310.2480.95350.1630.754 节点i123451.0+j0.01.0+j0.01.0+j0.01.0+j0.01.05+j0.0 表2.2 各节点电压（初值）标幺值参数 3 p-q分解计算方法3.1 形成节点导纳矩阵 自导纳的形成。对节点i其自导纳yii是节点i以外的所有节点都接地时节点i对地的总导纳。显然，yii应等于与节点i相接的各支路导纳之和，即式中，yi0为节点i与零电位节点之间的支路导纳；yij为节点i与节点j之间的支路导纳。互导纳的形成。对节点i与节点k之间的互导纳是节点i、k之间的支路导纳的负值，即不难理解。若节点i和k没有支路直接相连时，便有yi

5、k=0。含变压器支路的处理。若节点p、q间接有变压器，如下图所示，则可作出其型等值电路为：图1 变压器型等值电路则p、q的自导纳和节点间的互导纳分别为3.2 计算不平衡功率p、q并形成修正方程式对每一个pq节点或每一个pv节点都根据下列公式计算出有功功率增量p而对于每一个pq节点还可以根据下面的公式计算出无功功率增量q在有功功率增量和无功功率增量不满足如下约束条件时利用pq分解法则可以形成如下修正方程3.3 利用因子表法求解修正方程在电网计算中经常遇到这样的问题，对方程组需要反复多次求解，而每次求解仅改变常数项f，系数矩阵保持不变。按照一般的高斯消去法，对每一改变的常数项，形成包括常数项及系数

6、矩阵在内的增广矩阵，然后消去回代求出其解。可以看出，每次对增广矩阵中a矩阵元素的消元都是重复的，为了避免这种重复，我们把对相同的系数矩阵重复进行的消去与对不同的常数项进行的消去分开进行，因此对系数矩阵的消去只需进行一次，并在消去的过程中将对常数项进行消去运算的运算因子保存下来，形成所谓因子表，这就是因子表法。因为因子表记录了高斯消去法对常数项进行消去的全部信息，利用它便可对不同常数项进行消去，形成上三角矩阵，最后求出全部未知数。在使用pq分解法时，其系数矩阵是在迭代过程中保持不变的，所以为了节省内存和缩短运算时间我们采取了因子表法。同时由于电网的节点导纳矩阵矩阵是稀疏阵和对称阵，于是我们可以采

7、取只保存系数矩阵的上三角阵来使运算更为简化。若线性方程组一般形式如下：其中称为系数矩阵，称为未知数向量， 称为常数项向量。将矩阵a的元素进行如下处理：得到因子表其中 ；再利用因子表进行前代过程，求出每次迭代后的常数项。其前代公式是：求得向量；再由因子表与前代得到的向量f，得到方程组求解出此方程即可得到线性方程组的解向量。3.4 多次迭代最终求得v和以及全线路功率利用上面所介绍的方法求解修正方程组 可以求得和。再利用求得每次迭代后的结果。多次迭代当其满足约束条件和时，迭代结束。迭代结束后即可得到各节点的v和，再根据v、来计算pv节点的无功功率q和平衡节点的功率以及网络中的功率分布。pv节点及平衡

8、节点无功功率计算公式为：平衡节点有功功率计算公式为：以下图所标示的正方向，输电线路功率的计算公式如下：图2 支路功率计算对其进行实部虚部进行分解可得p、q计算公式为：4 matlab简介目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法的要求：计算速度快；内存需要少；计算结果有良好的可靠性和可信性；适应性好，亦即能处理变压器变比调整；系统元件的不同描述和与其它程序配合的能力强、简单。matlab是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图

9、等方面也具有强大的功能。matlab程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。matlab与c语言和fortran语言相比更容易被掌握。通过m语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。另外，matlab提供了一种特殊的工具：工具箱（toolboxes）.这些工具箱主要包括：信号处理（signal processing）、控制系统（control systems）、神经网络（neural networks）、模糊逻辑

10、(fuzzy logic)、小波(wavelets)和模拟(simulation)等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。matlab设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。矩阵的运算矩阵是matlab数据存储的基本单元，而矩阵的运算是matlab语言的核心，在matlab语言系统中几乎一切运算均是以对矩阵的操作为基础的。矩阵的基本数学运算包括矩阵的四则运算、与常数的运算、逆运算、行列式运算、秩运算、特征值运算等基本函数

11、运算，这里进行简单介绍。四则运算矩阵的加、减、乘运算符分别为“+，*” ，用法与数字运算几乎相同，但计算时要满足其数学要求 在matlab中矩阵的除法有两种形式：左除“”和右除“/”。在传统的matlab算法中，右除是先计算矩阵的逆再相乘，而左除则不需要计算逆矩阵直接进行除运算。通常右除要快一点，但左除可避免被除矩阵的奇异性所带来的麻烦。与常数的运算 常数与矩阵的运算即是同该矩阵的每一元素进行运算。但需注意进行数除时，常数通常只能做除数。基本函数运算矩阵的函数运算是矩阵运算中最实用的部分，常用的主要有以下几个：det(a) 求矩阵a的行列式eig(a) 求矩阵a的特征值inv(a)或a (-1

12、) 求矩阵a的逆矩阵rank(a) 求矩阵a的秩trace(a) 求矩阵a的迹（对角线元素之和）我们在进行工程计算时常常遇到矩阵对应元素之间的运算。这种运算不同于前面讲的数学运算，为有所区别，我们称之为数组运算。基本数学运算数组的加、减与矩阵的加、减运算完全相同。而乘除法运算有相当大的区别，数组的乘除法是指两同维数组对应元素之间的乘除法，它们的运算符为“.*”和“./”或“.”。前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系”的规定，数组与常数之间的除法运算没有任何限制。另外，矩阵的数组运算中还有幂运算（运算符为 . ）、指数运算（exp）、对数运算（log）、和开方

13、运算（sqrt）等。有了“对应元素”的规定，数组的运算实质上就是针对数组内部的每个元素进行的。矩阵的幂运算与数组的幂运算有很大的区别。逻辑关系运算 逻辑运算是matlab中数组运算所特有的一种运算形式，也是几乎所有的高级语言普遍适用的一种运算。 5 程序及说明 本程序采用matlab程序设计语言进行程序的设计。matlab程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。matlab与c语言和fortran语言相比更容易被掌握。通过m语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了

14、时间，从而可把主要的精力集中在算法的构思而不是编程上。matlab设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。5.1程序流程图开始读取节点数据读取线路数据并形成y矩阵形成矩阵b及b并求取因子表设定pq节点电压初值及各节点电压相角的初值迭代计算计算平衡节点功率及全部线路功率结束图5.1 程序主流程图kp=1，kq=1kp=0且kq=0kq=1count+kp=0计算p，p/nkp=0 ?求+=kp=0且kq=0kq=0计算q、q/nkq=0？求vv+=v结束开始图5.2 迭代部

15、分流程图5.2具体程序%电力系统pq分解法潮流计算disp(电力系统极坐标下的pq分解法潮流计算:);clearn=input(请输入结点数：n=);n1=input(请输入pv结点数：n1=);n2=input(请输入pq结点数：n2=);isb=input(请输入平衡结点：isb=);pr=input(请输入精确度：pr=);k=input(请输入变比矩阵看:k=);c=input(请输入支路阻抗矩阵：c=);y=input(请输入支路导纳矩阵：y=);u=input(请输入结点电压矩阵：u=);s=input(请输入各结点的功率：s=);z=zeros(1,n);n=zeros(n2,n

16、2+n1);l=zeros(n1+n2,n2);qt1=zeros(1,n1+n2);for m=1:n for r=1:n c(m,m)=c(m,m)+y(m,r); if k(m,r)=0 c(m,m)=c(m,m)+1/(k(m,r)*c(m,r)/(k(m,r)-1); c(r,r)=c(r,r)+1/(k(m,r)2*c(m,r)/(1-k(m,r); c(m,r)=c(m,r)*k(m,r); c(r,m)=c(m,r); endendendfor m=1:n for r=1:n if m=r z(m)=z(m)+1/c(m,r); end endendfor m=1:n for

17、r=1:n if m=r y(m,m)=c(m,m)+z(m); else y(m,r)=-1/c(m,r); end endenddisp(结点导纳矩阵:);disp(y);disp(迭代中关于b的矩阵:);g=real(y);b=imag(y);o=angle(u);u1=abs(u);k=0;pr=1;p=real(s);q=imag(s);while prpr for m=1:n2 ud(m)=u1(m); end for m=1:n1+n2 for r=1:n pt(r)=u1(m)*u1(r)*(g(m,r)*cos(o(m)-o(r)+b(m,r)*sin(o(m)-o(r);

18、end pt1(m)=sum(pt); pp(m)=p(m)-pt1(m); pp1(k+1,m)=pp(m); end for m=1:n2 for r=1:n qt(r)=u1(m)*u1(r)*(g(m,r)*sin(o(m)-o(r)-b(m,r)*cos(o(m)-o(r); end qt1(m)=sum(qt); qq(m)=q(m)-qt1(m); qq1(k+1,m)=qq(m); end pr1=max(abs(pp); pr2=max(abs(qq); pr=max(pr1,pr2); for m=1:n1+n2 for r=1:n1+n2 b1(m,r)=b(m,r);

19、end end for m=1:n2 for r=1:n2 b2(m,r)=b(m,r); end end jj=b1 l;n b2; disp(jj); for m=1:n1+n2 pp2(m)=pp(m)/u(m); end for m=1:n2 qq2(m)=qq(m)/u(m); end pq=pp2;qq2; da=-inv(jj)*pq; da1=da; for m=1:n1+n2 oo(m)=da1(m)/u(m); end for m=n:n1+n2+n2 uu(m-n1-n2)=da1(m); end for m=1:n1+n2 o(m)=o(m)+oo(m); end fo

20、r m=1:n2 u1(m)=u1(m)+uu(m); end for m=1:n1+n2 o(k+1,m)=180/pi*o(m); end for m=1:n2 u(k+1,m)=u1(m); end k=k+1;endfor m=1:n b(m)=u1(m)*cos(o(m); c(m)=u1(m)*sin(o(m);endu=b+i*c;for r=1:n ph1(r)=u(isb)*conj(y(isb,r)*conj(u(r);endph=sum(ph1);for m=1:n for r=1:n if m=r c1(m,r)=1/c(m,r); else c1(m,m)=c(m,m

21、); end endendfor m=1:n for r=1:n if (c(m,r)=inf)&(m=r) ss(m,r)=u1(m)2*conj(c1(m,m)+u(m)*(conj(u(m)-conj(u(r)*conj(c1(m,r); end endenddisp(迭代中的p：);disp(pp1);disp(迭代中的q：);disp(qq1);disp(迭代中相角:);disp(o);disp(迭代中电压的模:);disp(u);disp(平衡结点的功率:);disp(ph);disp(全部线路功率分布:);disp(ss);5.3 程序计算结果结点导纳矩阵:0.8381 - 3.

22、8179i -0.4044 + 1.6203i 0 0 -0.4337 + 2.2586i-0.4044 + 1.6203i 0.7769 - 3.4370i -0.3726 + 1.8557i 0 0 0 -0.3726 + 1.8557i 1.1608 - 7.0930i -0.5224 + 4.1792i -0.2739 + 1.2670i0 0 -0.5224 + 4.1792i 0.5499 - 4.3991i 0 -0.4337 + 2.2586i 0 -0.2739 + 1.2670i 0 0.7077 - 3.5257i迭代中关于b的矩阵: -3.8179 1.6203 0 0

23、 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.

24、6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0

25、 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370

26、1.8557 0 0 0 0 0 1.8557 -7.0930 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930 4.1792 0 0 0 0 0 4.1792 -4.3991 0 0 0 0 0 0 0 -3.8179 1.6203 0 0 0 0 0 1.6203 -3.4370 1.8557 0 0 0 0 0 1.8557 -7.0930迭代中的p： -0.1983 -0.1800 -0.2482 0.3225 0.0104 0.0183 -0.0040 0.0055 0.0037 -0.0004 0.0039 -0.0012 -0.0000