电力系统课程设计用c语言进行复杂电力系统的潮流计算

电力系统课程设计报告班级学号姓名2013年3月7日用C语言进行复杂电力系统的潮流计算电力系统的潮流计算是对电力系统分析的最基本步骤也是最重要的步骤，是指在一定的系统结构和运行条件下，确定系统运行状态的计算，也即是对各母线（节点）电压，各元件（支路）传输电线或功率的计算。通过计算出的节点电压和功率分布用以检查系统各元件是否过负荷，各点电压是否合理，以及功率损耗等。即使对于一个简单的电力系统，潮流计算也不是一件简单就可以完成的事，其运算量很大，因此如果对于一个大的、复杂的电网来说的话，由于其节点多，分支杂，其计算量可想而知，人工对其计算也更是难上加难了。特别是在现实生活中，遇到一个电力系统不会像我们期望的那样可以知道它的首端电压和首端功率或者是末端电压和末端功率，而是只知道它的首端电压和末端功率，更是使计算变的头疼万分。为了使计算变的简单，我们就可以利用计算机，用C语言编程来实现牛顿拉夫逊（NEWTONRAPHSON）迭代法，最终实现对电力系统潮流的计算。一用牛顿拉夫逊迭代法进行电力系统潮流计算的相关概念1节点导纳矩阵如图所示的电力网络,将节点I和J的电压用I和J表示,它们之间的支路导纳表示为UYIJ,那么有基尔霍夫电流定律可知注入接点I的电流I设流入节点的电流为正等于离开节点I的电流之和,因此有JIJIII1100JINIJIJNIJIJIUIIY12YUINIJJNIJI00如令则可将12改写为INIJJYY0IJIJYI1,2,N13NIJIJIUYI1上式也可以写为IYU14其中Y为节点导纳矩阵，也称为稀疏的对称矩阵，它是NN阶方阵。对角元YII称为自导纳，它等于与该节点I直接相连的所有支路导纳总和；非对角元YIJ（IJ）称为互导纳或转移导纳，它等于连结节点I，J支路导纳的负数，且有YIJYJI，当节点I，J之间没有支路直接相连时，YIJYJI0。电力系统的分析计算中，往往要作不同的运行方式下的潮流计算，如果系统发生变化，如投切一条线路或一台变压器，由于改变了一条支路的状态或参数只影响该支路两端节点的自导纳和他们之间的互导纳，因而对每一种运行方式不必重新形成导纳矩阵，只需对原有导纳矩阵作相应的修改即可。2潮流计算的功率方程在实际的电力系统中，已知的条件往往不是节点的注入电流而是负荷和发电机功率，且这些功率一般不随节点的电压变化而变化，而节点的电流则是随电压的变化而变化的，因此在已知节点导纳矩阵的情况下，必须用已知的节点功率来替代未知的节点注入电流，才能求出节点电压，每一个节点的注入功率方程式为IIIUGIIGIJQPSLIILIJQPS（15）LIGII（16）IIIQ17NJJIIUYI1节点注入电流用功率和电压表示为（18）IILIGILIGIIIUJQPJPSI功率方程可以表示为（19）NJJIIUYQ13节点分类对于有N个节点的电力网络，可以列出N个功率方程，由图可知一个节点有四个变量注入有功功率PI，注入无功功率QI，节点电压幅值UI和相角。N个节点有4N个变量,但IG只有2N个关系式,所以为了使潮流有确定解,必须给定其中2N个变量。根据给定节点变量的不同，可以有以下三种类型的节点（1）PQ节点给定注入功率PI，QI，即已知PGI，PLI，QGI，QLI，待求UI，I。例如降压变电所母线（负荷节点），固定出力的发电厂母线。（2）PV节点给定了注入有功功率PI（PGI，PLI），UI和QLI，待求QGI（QI），I。例如有一定无功电源的降压变电所母线，有一定储备的发电厂母线。（3）平衡节点给定了UI，I和PLI，QLI，待求PGI，QGI，即PI，QI，用来平衡全电网的功率，通常在一个独立的电力系统中只设一个平衡节点。4牛顿拉夫逊迭代法牛顿拉夫逊迭代法将解非线性方程组的过程转化为反复求与之相对应的线性方程的求解过程。对于一个N维非线性方程组N1,2,3,NINIYXF,21假定其初值为X10,X20,XN0,也即其近似解,它与真值之间的误差为也即各变量与真解之间的修正量。0021,NX将这N个方程式都在初值的附近展开成TAYLOR级数且忽略二次项及高次项,则可得修正方程,I1,2,N110INIYXXF0010N021FFX,X将修正方程写成矩阵形式111其中令J,称之为雅可比JACOBI方阵。0N01N0N101XFFFF它的第I行，第J列交点的元素为第I个函数对第J个变量XJ的偏导数在点,21NIXF（X10,X20,XN0）的值，所以方程组是线性方程，可用于求出，0021,N从而得到新的近似解，（112）001IIIXX于是得到一般迭代式01N01NN1010N0211XFFFFX,X,NFYF（113）1N1NN11KNK211XFFFFX,X,KKKKKKNFYF于是得到近似解（114）1KIKIKIX迭代一直进行到MAX|YIFIX10,X20,XN0|K1KKM四用编程方法求解实际问题如图所示的一个电力网络，如何来求解它的潮流分布呢已知为定值，其余四个节点都是PQ节点，且给定的注入功率分别061JU为106,054,154,2032JSJSJSJS程序清单如下INCLUDEINCLUDEFLOATDIVREB1,B2,B3,B4FLOATB1,B2,B3,B4FLOATA1RA1RB1B3B2B4/B3B3B4B4RETURNA1RFLOATDIVIMB1,B2,B3,B4FLOATB1,B2,B3,B4FLOATA1IA1IB2B3B1B4/B3B3B4B4RETURNA1IFLOATMULREB1,B2,B3,B4FLOATB1,B2,B3,B4FLOATA2RA2RB1B3B2B4RETURNA2RFLOATMULIMB1,B2,B3,B4FLOATB1,B2,B3,B4FLOATA2IA2IB2B3B1B4RETURNA2IFLOATMAXFLOATA,INTNINTIFLOATMAXFORI0IAI1MAXAIAIAI1AI1MAXRETURNMAXMAININTI,J,K,N,KMFLOATEPS,SUMPI1,SUMPI2,SUMQI1,SUMQI2,MAX,SUMIR,SUMII,I1R,I1IFLOATPI05,QI05,DETPI5,DETQI5,IIR05,III05,J088,DETSI8,DETUI8,U88,L88,Y8,UI18,H44,N44,J44,L44,EI15,FI15STATICFLOATYBR556250,5000,1250,0,0,5000,10834,1667,1667,2500,1250,1667,12917,10000,0,0,1667,10000,12917,1250,0,2500,0,1250,3750STATICFLOATYBI5518750,15000,3750,0,0,15000,32500,5000,5000,7500,3750,5000,38750,30000,0,0,5000,30000,38750,3750,0,7500,0,3750,11250FLOATEI05106,10,10,10,10FLOATFI050,0,0,0,0FLOATPI50,02,045,04,06FLOATQI50,02,015,005,01K0KM6EPS000001DOK1PRINTF“NOWSTARTN“PRINTF“THEDTIMESN“,KFORI1I0IDETUIIYIFORJI1JEPSPRINTF“ALLDODTIMESN“,KSUMIR0SUMII0FORI0I5II1RMULREYBR0I,YBI0I,EI0I,FI0II1IMULIMYBR0I,YBI0I,EI0I,FI0ISUMIRI1RSUMIII1IPI0MULREEI00,FI00,SUMIR,SUMIIQI0MULIMEI00,FI00,SUMIR,SUMIIPRINTF“S1FJFN“,PI0,QI0EI10EI00FI10FI00FORI0I5IPRINTF“UDFFN“,I1,SQRTEI1IEI1IFI1IFI1I,ATANFI1I/EI1I180/314159最终求出平衡节点的功率和线路的功率以及各个节点的电压为247098161JIJSR06U396024712853104U924715利用这中方法进行潮流的计算，其收敛的速度是显而易见的，如上面这到题，在经过三次迭代以后已经达到我们要求的水平，但必须建立在有一个比较优的初值的基础上。实验总结通过这次为期两周的课程设计，让我知道了潮流计算的基本步骤和方法，要让电力系统稳定运行，必须进行精密的计算和设计，明白了潮流计算对于电力系统的重要性，尤其是对牛顿拉夫逊法进行潮流计算有了深刻的理解。在用牛顿拉夫逊法进行潮流计算的过程中，用到了导纳矩阵，雅可比方阵等数学方法，用建立的数学模型结合C语言进行编程，这就要求掌握C语言的编程方法，还要熟悉VISUALC编程软件，复习以前学习过的知识，查阅相关的资料，完成课程设计。