PQ分解法在潮流分析中应用Word文件下载.docx

1、 1 or 1998 by H. Saadatj=sqrt（-1。 i = sqrt（-1nl = linedata（:,1 nr = linedata（:,2 R = linedata（:,3zvpgeqJ1hkX = linedata（:,4 Bc = j*linedata（:,5 a = linedata（:, 6NrpoJac3v1nbr=length（linedata（: nbus = max（max（nl, max（nr1nowfTG4KIZ = R + j*X。 y= ones（nbr,1./Z。 %branch admittancefjnFLDa5Zofor n = 1:nbr

2、if a（n = 1。 else endYbus=zeros（nbus,nbus % initialize Ybus to zerotfnNhnE6e5 % formation of the off diagonal elementsfor k=1:nbr。 Ybus（nl（k,nr（k=Ybus（nl（k-y（k/a（k Ybus（nr（k,nl（k endend % formation of the diagonal elementsfor n=1:nbus for k=1: if nl（k=n Ybus（n,n = Ybus（n,n+y（k/（a（k2 + Bc（k elseif nr（

3、k +Bc（k else, endclear Pgg输出母线的潮流结果程序busout：% This program prints the power flow solution in a tabulated formHbmVN777sL% on the screen.% Copyright （C 1998 by H. Saadat.%clcdisp（techfprintf（ Maximum Power Mismatch = %g n, maxerrorV7l4jRB8Hs No. of Iterations = %g nn, iter83lcPA59W9head = Bus Voltage

4、Angle -Load- -Generation- InjectedmZkklkzaaP No. Mag. Degree MW Mvar MW Mvar Mvar AVktR43bpw。ORjBnOwcEddisp（head fprintf（ %5g, n, fprintf（ %7.3f, Vm（n, %8.3f, deltad（n %9.3f, Pd（n,2MiJTy0dTT, Qd（n, Pg（n,gIiSpiue7A %9.3f , Qg（n %8.3fn, Qsh（nuEh0U1Yfmh n Total IAg9qLsgBX, Pdt, Qdt, Pgt, Qgt %9.3fnn, Q

5、shtWwghWvVhPE计算各条支路潮流分布和网损的程序lineflow：% This program is used in conjunction with lfgauss or lf Newtonasfpsfpi4k% for the computation of line flow and line losses.ooeyYZTjj1% 1998 H. SaadatSLT = 0。n Line Flow and Losses nnBkeGuInkxI -Line- Power at bus & line flow -Line loss- TransformernPgdO0sRlMo f

6、rom to MW Mvar MVA MW Mvar tapn3cdXwckm15busprt = 0。 for L = 1: if busprt = 0%6g, P（n*basemvah8c52WOngM%9.3f, Q（n%9.3fn, abs（S（nv4bdyGious busprt = 1。 if nl（L=n k = nr（L In = （V（n - a（L*V（k*y（L/a（L2 + Bc（L2*V（nJ0bm4qMpJ9 Ik = （V（k - V（n + Bc（L Snk = V（n*conj（In*basemva。 Skn = V（k*conj（Ik SL = Snk +

7、Skn。 SLT = SLT + SL。 elseif nr（L=n k = nl（L - V（k*V（n2*V（kXVauA9grYP=n | nr（L%12g, k, real（Snk, imag（SnkbR9C6TJscw, abs（Snk, real（SL =n & a（L = 1, imag（SL, a（LpN9LBDdtrd else, fprintf（SLT = SLT/2。 Total loss DJ8T7nHuGT, real（SLT, imag（SLTQF81D7bvUAclear Ik In SL SLT Skn Snk2、程序运行 在matlab中运行以上程序，得到各个

8、节点的电压幅值和相位，各条支路的潮流和网损，即潮流分析计算的结果：4B7a9QFw9h1）改变负荷侧无功补偿容量 根据实验要求，分析改变无功补偿容量后对潮流和网损的影响，于是我们增大负荷节点10处的无功补偿容量，将19Mvar改为38Mvar，对应数据如下：ix6iFA8xoX basemva = 100。wt6qbkCyDEKp5zH46zRk% Bus Bus Voltage Angle -Load- -Generator- Static MvarYl4HdOAA61% No code Mag. Degree MW Mvar MW Mvar Qmin Qmax +Qc/-Qlch4PJx4

9、BlIbusdata=1 1 1.06 0.0 0.0 0.0 0.0 0.0 0 0 0qd3YfhxCzo 2 2 1.043 0.0 21.70 12.7 40.0 0.0 -40 50 0E836L11DO5 3 0 1.0 0.0 2.4 1.2 0.0 0.0 0 0 0S42ehLvE3M 4 0 1.06 0.0 7.6 1.6 0.0 0.0 0 0 0501nNvZFis 5 2 1.01 0.0 94.2 19.0 0.0 0.0 -40 40 0jW1viftGw9 6 0 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0xS0DOYWHLP 7 0 1.0

10、0.0 22.8 10.9 0.0 0.0 0 0 0LOZMkIqI0w 8 2 1.01 0.0 30.0 30.0 0.0 0.0 -30 40 0ZKZUQsUJed 9 0 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0dGY2mcoKtT 10 0 1.0 0.0 5.8 2.0 0.0 0.0 -6 24 38rCYbSWRLIA 11 2 1.082 0.0 0.0 0.0 0.0 0.0 0 0 0FyXjoFlMWh 12 0 1.0 0 11.2 7.5 0 0 0 0 0TuWrUpPObX 13 2 1.071 0 0 0.0 0 0 -6 24 07qW

11、Aq9jPqE 14 0 1 0 6.2 1.6 0 0 0 0 0llVIWTNQFk 15 0 1 0 8.2 2.5 0 0 0 0 0yhUQsDgRT1 16 0 1 0 3.5 1.8 0 0 0 0 0MdUZYnKS8I 17 0 1 0 9.0 5.8 0 0 0 0 009T7t6eTno 18 0 1 0 3.2 0.9 0 0 0 0 0e5TfZQIUB5 19 0 1 0 9.5 3.4 0 0 0 0 0s1SovAcVQM 20 0 1 0 2.2 0.7 0 0 0 0 0GXRw1kFW5s 21 0 1 0 17.5 11.2 0 0 0 0 0UTREx

12、49Xj9 22 0 1 0 0 0.0 0 0 0 0 08PQN3NDYyP 23 0 1 0 3.2 1.6 0 0 0 0 0mLPVzx7ZNw 24 0 1 0 8.7 6.7 0 0 0 0 4.3AHP35hB02d 25 0 1 0 0 0.0 0 0 0 0 0NDOcB141gT 26 0 1 0 3.5 2.3 0 0 0 0 01zOk7Ly2vA 27 0 1 0 0 0.0 0 0 0 0 0fuNsDv23Kh 28 0 1 0 0 0.0 0 0 0 0 0tqMB9ew4YX 29 0 1 0 2.4 0.9 0 0 0 0 0HmMJFY05dE 30 0

13、 1 0 10.6 1.9 0 0 0 0 0。ViLRaIt6sk运行上述程序，得到潮流分析计算的结果：Total loss 17.580MW21.688MVar无功补偿改为两倍即38MVar时）Totalloss17.598MW 22.245MVarTotalloss 17.598MW 22.245MVar 由以上的潮流计算结果看出，减少电力系统的回路数，系统的潮流发生了变化，系统总网损发生了很大的变化，总的有功损耗减小，可是总的无功损耗确增大了，即电力系统的稳定性变差了。WHF4OmOgAw4、实验分析在matlab中运行以上程序，得到各个节点的电压幅值和相位，各条支路的潮流和网损。此程

14、序的运行结果和实验程序给出的结果是一致的。说明程序无误，但在精度上有微小差异，这主要是和导纳矩阵的精度以及现实精度有关。通过改变支路中的数据，可以减小系统的回路数，其中采用的是将部分支路的电阻修改为无穷大，也相当于相应断路器断开。aDFdk6hhPd5、实验心得 本程序分模块进行，先是排序，再是求导纳矩阵，然后求雅阁比，再进行迭代运算。程序本身很简洁明了，运行的时候只要在matlab中输入main就行了，然后打开bus和line所在的m文件中，结果就自动存在rustl文件中。通过编写牛顿拉夫逊法matlab潮流计算程序复习了潮流计算知识，也实现了计算机算法。ozElQQLi4T电压是电能质量的重要指标之一，网损是电力企业的一项重要综合性技术经济指标。长期以来电力系统网络损耗问题比较突出，而无功补偿是降低线损的有效手段。随着电力系统负荷的增加，对无功功率的需求也日益增加。在电网中的适当位置装设无功补偿装置成为满足电网无功需求的必要手段。提高电压质量、降损节能是供电企业当前的工作重点之一。对电网给予合理的无功补偿，能够有效地维持系统的电压水平，降低有功网损，提高网络输送容量满足负荷的无功需求。节点负荷