C语言进行潮流计算

-1-电力系统课程设计C语言潮流计算学院：电气工程班级：电092班学号：0912002020学生姓名：闵凯2013.3.7-2-电力系统的潮流计算是对电力系统分析的最基本步骤也是最重要的步骤，是指在一定的系统结构和运行条件下，确定系统运行状态的计算，也即是对各母线（节点）电压，各元件（支路）传输电线或功率的计算。通过计算出的节点电压和功率分布用以检查系统各元件是否过负荷，各点电压是否合理，以及功率损耗等。即使对于一个简单的电力系统，潮流计算也不是一件简单就可以完成的事，其运算量很大，因此如果对于一个大的、复杂的电网来说的话，由于其节点多，分支杂，其计算量可想而知，人工对其计算也更是难上加难了。特别是在现实生活中，遇到一个电力系统不会像我们期望的那样可以知道它的首端电压和首端功率或者是末端电压和末端功率，而是只知道它的首端电压和末端功率，更是使计算变的头疼万分。为了使计算变的简单，我们就可以利用计算机，用C语言编程来实现牛顿-拉夫逊（Newton-Raphson）迭代法，最终实现对电力系统潮流的计算。一．用牛顿-拉夫逊迭代法进行电力系统潮流计算的相关概念1.节点导纳矩阵如图所示的电力网络,将节点i和j的电压用Ui和.Uj表示,它们之间的支路导纳表示为yij,那么有基尔霍夫电流定律可知注入接点I的电流.Ii(设流入节点的电流为正)等于离开节点I的电流之和,因此有jijIIiI)(.00jinijijnijijiUUIIy(1-1)UyyUInijijnijiji00(1-2)如令iinijijYy0ijijYy则可将(1-2)改写为:nijijijiUYI1I=1,2,…,n.(1-3)上式也可以写为:I=YU(1-4)-3-其中Y为节点导纳矩阵，也称为稀疏的对称矩阵，它是n×n阶方阵。对角元Yii称为自导纳，它等于与该节点I直接相连的所有支路导纳总和；非对角元Yij（i≠j）称为互导纳或转移导纳，它等于连结节点I，j支路导纳的负数，且有Yij=Yji，当节点I，j之间没有支路直接相连时，Yij=Yji=0。电力系统的分析计算中，往往要作不同的运行方式下的潮流计算，如果系统发生变化，如投切一条线路或一台变压器，由于改变了一条支路的状态或参数只影响该支路两端节点的自导纳和他们之间的互导纳，因而对每一种运行方式不必重新形成导纳矩阵，只需对原有导纳矩阵作相应的修改即可。2．潮流计算的功率方程在实际的电力系统中，已知的条件往往不是节点的注入电流而是负荷和发电机功率，且这些功率一般不随节点的电压变化而变化，而节点的电流则是随电压的变化而变化的，因此在已知节点导纳矩阵的情况下，必须用已知的节点功率来替代未知的节点注入电流，才能求出节点电压，每一个节点的注入功率方程式为：IiiiUUGiGiGijQPS~LiLiLijQPS~LiGiiPPP（1-5）LiGiiQQQ（1-6）njjijiUYI1(1-7)节点注入电流用功率和电压表示为：iiiiLiGiLiGiiiiUjQPUQQjPPUSIˆˆ)()(ˆˆ（1-8）功率方程可以表示为：njjijiiiUYUjQP1ˆ（1-9）3．节点分类对于有n个节点的电力网络，可以列出n个功率方程，由图可知一个节点有四个变量：注入有功功率Pi，注入无功功率Qi，节点电压幅值Ui和相角i。n个节点有4n个变量,但只有2n个关系式,所以为了使潮流有确定解,必须给定其中2n个变量。根据给定节点变量的不同，可以有以下三种类型的节点：（1）PQ节点：给定注入功率Pi，Qi，即已知PGi，PLi，QGi，QLi，待求Ui，δi。例如：降压变电所母线（负荷节点），固定出力的发电厂母线。G-4-（2）PV节点：给定了注入有功功率Pi（PGi，PLi），Ui和QLi，待求QGi（Qi），δi。例如：有一定无功电源的降压变电所母线，有一定储备的发电厂母线。（3）平衡节点：给定了Ui，δi和PLi，QLi，待求PGi，QGi，即Pi，Qi，用来平衡全电网的功率，通常在一个独立的电力系统中只设一个平衡节点。4．牛顿-拉夫逊迭代法牛顿-拉夫逊迭代法将解非线性方程组的过程转化为反复求与之相对应的线性方程的求解过程。对于一个n维非线性方程组：iniyxxxf),...,,(21n=1,2,3,…,n假定其初值为x1(0),x2(0),…,xn(0),也即其近似解,它与真值之间的误差为)0()0(2)0(1,...,,nxxx也即各变量与真解之间的修正量。将这n个方程式都在初值的附近展开成Taylor级数且忽略二次项及高次项,则可得修正方程iniyxxf)0(0nn)0(1011(0)n(0)2(0)1xf...xf)x,,x,x(,I=1,2,…n.(1-10)将修正方程写成矩阵形式:(1-11)其中令J=0nn01n0n1011xf...xf.........xf...xf,称之为雅可比(Jacobi)方阵。它的第I行，第j列交点的元素为第I个函数),...,(21nixxxf对第j个变量xj的偏导数在点（x1(0),x2(0),…,xn(0)）的值，所以方程组是线性方程，可用于求出)0()0(2)0(1,...,,nxxx，从而得到新的近似解，)0()0()1(iiixxx（1-12）于是得到一般迭代式：)0()0(10nn01n0n1011(0)n(0)2(0)1(0)n(0)2(0)111...xf...xf.........xf...xf)x,,x,x(:)x,,x,x(nnnxxfyfy-5-)()(1nn1nn111(k)n(k)2(k)1(k)n(k)2(k)111...xf...xf.........xf...xf)x,,x,x(:)x,,x,x(knkkkkknnxxfyfy（1-13）于是得到近似解：)()()1(kikikixxx（1-14）迭代一直进行到Max{|yi-fi(x1(0),x2(0),…,xn(0))|}ε或Max{|Δxi(k)|}ε为止。二．用牛顿-拉夫逊迭代法进行潮流计算设网络中除参考节点外有n个节点，其中1个平衡节点（并令第n个节点为平衡节点），m个PQ节点（为第1~m个节点），有n-m-1个PV节点。运用牛顿-拉夫逊法直接求解功率方程（1-9），并将Yij=Gij+jBij及...iiifeU代入得：iiiinjijijiijQPjfejBGjfe))((1将实部和虚部分开得：njijijjijijijjijiPeBfGffBeGe1)]()([(2-1)njijijjijiiijjijiQeBfGefBeGf1)]()([(2-2)此外由于系统中还有PV节点，所以还应补充一组方程：222iiiUfe(2-3)在式（2-3）中，ei、fi分别为迭代过程中求得的节点电压的实部和虚部，Pi为PQ节点和PV节点的注入有功功率，Qi为PQ节点的注入无功功率，Ui为PV节点的电压大小。由（2-1），（2-2），（2-3）三式所组成的方程组一共有2（n-1）个独立方程，其中，式（2-1）类型的有(n-1)个,包括除平衡节点以外所有节点的有功功率Pi的表达式;式(3-2)类型的有(m-1)个,包括所有PQ节点无功功率Qi的表达式,式(2-3)类型的有(n-m)个,包括所有PV节点电压Ui2的表达式.平衡节点s的功率和电压之所以不包括在这方程组内,是由于平衡节点的注入功率不可能事先给定,从而不可能列出Ps,Qs的表达式,而平衡节点的电压sssjfeU则不必求取.于是建立修正方程式如下:-6-nnppnnnnnpnpnnnnnnnnnpnpnnnnpnpnpppppppppnpnppppppppnnppnnppnnppnnppnnppefefefefSRSRSRSRNHNHNHNHSRSRSRSRNHNHNHNHLJLJLJLJNHNHNHNHLJLJLJLJNHNHNHNHUPUPQPQP22112211221122112211222222222121222222222121111112121111111112121111222211(2-4)式中的njjijjijijijjijiiieBfGffBeGePP1)]()([(2-5)njjijjijiiijjijiiieBfGefBeGfQQ1)]()([(2-6))(2222iiiifeUU(2-7)当j≠I时雅可比方阵的各个元素分别为：iijiijjiijfGeBfPH；iijiijjiijfBeGePN；ijjiijNfQJ；ijiiijHeQL；;021jijfUR021jijeUS当j=I时雅可比方阵的各个元素分别为：iiiiiiiiiiiibfGeBfPH;iiiiiiiiiiiafBeGePN;iiiiiiijiiiafBeGeQN;iiiiiiiiiiiibfGeBfPH;iiiiiffUR22;iiiiieeUS22其中：)()(1nijiijiijiiiiiiiifBfGfBeGa-7-nijjijjijjiijiiiieBfGfBeGb1)()(三．用C语言编程计算潮流的流程图否是否是开始始输入原始数据形成节点导纳矩阵YB设节点电压初值ei(0),fi(0)置迭代次数k=1计算ΔPi(k),ΔQi(k),ΔUi(k)2形成雅可比矩阵求各节点电压变量Δei(k),Δfi(k)求出|Δe(k)|max、|Δf(k)|max迭代是否收敛，|Δe(k)|max、|Δf(k)|max≤ε？计算平衡节点功率sS~和线路功率ijS~停止计算各节点电压新值：)()()1(kikikieee)()()1(kikikifff增大迭代次数，k-k+1kkm?-8-四．用编程方法求解实际问题如图所示的一个电力网络，G0.08+j.240.01+j0.030.06+j0.180.06+j0.180.04+j0.120.08+j0.240.02+j0.060.45+j0.150.4+j0.050.6+j0.1-(0.2+j0.2)123450G0.08+j.240.01+j0.030.06+j0.180.06+j0.180.04+j0.120.08+j0.240.02+j0.060.45+j0.150.4+j0.050.6+j0.1-(0.2+j0.2)123450已知：006.11jU为定值，其余四个节点都是PQ节点，且给定的注入功率分别为：10.060.0~,05.040.0~,15.045.0~,20.020.0~5432jSjSjSjS由上图可得相应的节点导纳矩阵Y=[6.25-18.75i-5+15i-1.25+3.75i00;-5+15i10.834-32.5i-1.667+5i-1.667+5i-2.5+7.5i;-1.25+3.75i-1.667+5i12.917-38.75i-10+30i0;0-1.667+5i-10+30i12.917-38.7