PQ分解法在潮流分析中应用Word文件下载.docx
《PQ分解法在潮流分析中应用Word文件下载.docx》由会员分享，可在线阅读，更多相关《PQ分解法在潮流分析中应用Word文件下载.docx（10页珍藏版）》请在冰豆网上搜索。
1or<
1tr.tapatbusnl0YujCfmUCw
linedata=[120.01920.05750.026401
130.04520.18520.020401
240.05700.17370.018401
340.01320.03790.004201
250.04720.19830.020901
260.05810.17630.018701
460.01190.04140.004501
570.04600.11600.010201
670.02670.08200.008501
680.01200.04200.004501
690.00.20800.00.978
6100.556000.969
9110.208001
9100.110001
4120.256000.932
12130.140001
lfybus%formthebusadmittancematrixeUts8ZQVRd
lfDecouple%LoadflowsolutionbydecouplemethodsQsAEJkW5T
busout%PrintsthepowerflowsolutiononthescreenGMsIasNXkA
lineflow%ComputesanddisplaysthelineflowandlossesTIrRGchYzg
以上是网络支路中的数据，通过首、末端节点来表示，包括了支路的电阻、电抗、0.5的电纳和变比标幺值，其中标幺值为1时表示线路，大于1时表示降压变压器，小于1时表示升压变压器。
其后的一系列是对函数的声明，lfDecouple是形成节点导纳矩阵程序，lfDecouple是PQ分解法程序，busout是输出母线的潮流结果程序，算出各个节点的电压的幅值和相位，lineflow是计算各条支路潮流分布和网损的程序。
7EqZcWLZNX
节点导纳矩阵程序lfbus：
%ThisprogramobtainsthBusAdmittanceMatrixforpowerflowsolutionlzq7IGf02E
%Copyright（c>
1998byH.Saadat
j=sqrt（-1>
。
i=sqrt（-1>
nl=linedata（:
1>
nr=linedata（:
2>
R=linedata（:
3>
zvpgeqJ1hk
X=linedata（:
4>
Bc=j*linedata（:
5>
a=linedata（:
6>
NrpoJac3v1
nbr=length（linedata（:
>
nbus=max（max（nl>
max（nr>
1nowfTG4KI
Z=R+j*X。
y=ones（nbr,1>
./Z。
%branchadmittancefjnFLDa5Zo
forn=1:
nbr
ifa（n>
<
=0a（n>
=1。
elseend
Ybus=zeros（nbus,nbus>
%initializeYbustozerotfnNhnE6e5
%formationoftheoffdiagonalelements
fork=1:
nbr。
Ybus（nl（k>
nr（k>
=Ybus（nl（k>
-y（k>
/a（k>
Ybus（nr（k>
nl（k>
end
end
%formationofthediagonalelements
forn=1:
nbus
fork=1:
ifnl（k>
==n
Ybus（n,n>
=Ybus（n,n>
+y（k>
/（a（k>
^2>
+Bc（k>
elseifnr（k>
+Bc（k>
else,end
clearPgg
输出母线的潮流结果程序busout：
%ThisprogramprintsthepowerflowsolutioninatabulatedformHbmVN777sL
%onthescreen.
%Copyright（C>
1998byH.Saadat.
%clc
disp（tech>
fprintf（'
MaximumPowerMismatch=%g\n'
maxerror>
V7l4jRB8Hs
No.ofIterations=%g\n\n'
iter>
83lcPA59W9
head=['
BusVoltageAngle------Load---------Generation---Injected'
mZkklkzaaP
'
No.Mag.DegreeMWMvarMWMvarMvar'
AVktR43bpw
]。
ORjBnOwcEd
disp（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.3f\n'
Qsh（n>
uEh0U1Yfmh
\n'
Total'
IAg9qLsgBX
Pdt>
Qdt>
Pgt>
Qgt>
%9.3f\n\n'
Qsht>
WwghWvVhPE
计算各条支路潮流分布和网损的程序lineflow：
%ThisprogramisusedinconjunctionwithlfgaussorlfNewtonasfpsfpi4k
%forthecomputationoflineflowandlinelosses.ooeyYZTjj1
%
1998H.Saadat
SLT=0。
\n'
LineFlowandLosses\n\n'
BkeGuInkxI
--Line--Poweratbus&
lineflow--Lineloss--Transformer\n'
PgdO0sRlMo
fromtoMWMvarMVAMWMvartap\n'
3cdXwckm15
busprt=0。
forL=1:
ifbusprt==0
%6g'
P（n>
*basemva>
h8c52WOngM
%9.3f'
Q（n>
%9.3f\n'
abs（S（n>
v4bdyGious
busprt=1。
ifnl（L>
==nk=nr（L>
In=（V（n>
-a（L>
*V（k>
*y（L>
/a（L>
^2+Bc（L>
^2*V（n>
J0bm4qMpJ9
Ik=（V（k>
-V（n>
+Bc（L>
Snk=V（n>
*conj（In>
*basemva。
Skn=V（k>
*conj（Ik>
SL=Snk+Skn。
SLT=SLT+SL。
elseifnr（L>
==nk=nl（L>
-V（k>
*V（n>
^2*V（k>
XVauA9grYP
==n|nr（L>
%12g'
k>
real（Snk>
imag（Snk>
bR9C6TJscw
abs（Snk>
real（SL>
==n&
a（L>
~=1
imag（SL>
a（L>
pN9LBDdtrd
else,fprintf（'
SLT=SLT/2。
Totalloss'
DJ8T7nHuGT
real（SLT>
imag（SLT>
QF81D7bvUA
clearIkInSLSLTSknSnk
2、程序运行
在matlab中运行以上程序，得到各个节点的电压幅值和相位，各条支路的潮流和网损，即潮流分析计算的结果：
4B7a9QFw9h
<
1）改变负荷侧无功补偿容量
根据实验要求，分析改变无功补偿容量后对潮流和网损的影响，于是我们增大负荷节点10处的无功补偿容量，将19Mvar改为38Mvar，对应数据如下：
ix6iFA8xoX
basemva=100。
wt6qbkCyDE
Kp5zH46zRk
%BusBusVoltageAngle---Load-----------Generator-----StaticMvarYl4HdOAA61
%NocodeMag.DegreeMWMvarMWMvarQminQmax+Qc/-Qlch4PJx4BlI
busdata=[111.060.00.00.00.00.0000qd3YfhxCzo
221.0430.021.7012.740.00.0-40500E836L11DO5
301.00.02.41.20.00.0000S42ehLvE3M
401.060.07.61.60.00.0000501nNvZFis
521.010.094.219.00.00.0-40400jW1viftGw9
601.00.00.00.00.00.0000xS0DOYWHLP
701.00.022.810.90.00.0000LOZMkIqI0w
821.010.030.030.00.00.0-30400ZKZUQsUJed
901.00.00.00.00.00.0000dGY2mcoKtT
1001.00.05.82.00.00.0-62438rCYbSWRLIA
1121.0820.00.00.00.00.0000FyXjoFlMWh
1201.0011.27.500000TuWrUpPObX
1321.071000.000-62407qWAq9jPqE
140106.21.600000llVIWTNQFk
150108.22.500000yhUQsDgRT1
160103.51.800000MdUZYnKS8I
170109.05.80000009T7t6eTno
180103.20.900000e5TfZQIUB5
190109.53.400000s1SovAcVQM
200102.20.700000GXRw1kFW5s
2101017.511.200000UTREx49Xj9
2201000.0000008PQN3NDYyP
230103.21.600000mLPVzx7ZNw
240108.76.700004.3AHP35hB02d
2501000.000000NDOcB141gT
260103.52.3000001zOk7Ly2vA
2701000.000000fuNsDv23Kh
2801000.000000tqMB9ew4YX
290102.40.900000HmMJFY05dE
3001010.61.900000]。
ViLRaIt6sk
运行上述程序，得到潮流分析计算的结果：
Totalloss17.580MW21.688MVar<
无功补偿改为两倍即38MVar时）
Totalloss17.598MW22.245MVar<
无功补偿为19MVar时）
由上面潮流计算的结果看出，增大负荷无功补偿的容量，系统各节点的潮流均发生了变化，系统总的网损，不论是有功功率损耗，还是无功功率损耗，都减小了，即提高了线路输送功率的能力。
9eK0GsX7H1
（2）改变回路数
在实际电力系统中有可能出现回路数变化的现象，通过改变支路中的数据，可以减小系统的回路数，其中采用的是将部分支路的电阻修改为无穷大inf，也相当于相应断路器断开的情况，其数据如下:
naK8ccr8VI
%BusbusRX1/2B=1forlinesB6JgIVV9ao
1tr.tapatbusnlP2IpeFpap5
26inf0.17630.018701
46inf0.04140.004501
1214inf.255901
1215.0662.130401
1216.0945.198701
1415.2210.199701
1617inf.192301
1518.1073.218501
1819.0639.129201
1920.0387.068001
1020.0167.209001
1017.0324.084501
1021.0348.074901
1022.0727.149901
2122.0116.023601
1523.0227.202001
2224.1150.179001
2324.1320.270001
2425.1885.329201
2526.0255.380001
2527.1093.208701
28270.396000.968
2729.2198.415301
2730.3202.602701
2930.2399.453301
828.0636.20000.02141
628.0169.05990.0651]。
lfybus%formthebusadmittancematrix3YIxKpScDM
lfgauss%LoadflowsolutionbyGauss-SeidelmethodgUHFg9mdSs
busout%PrintsthepowerflowsolutiononthescreenuQHOMTQe79
lineflow%ComputesanddisplaysthelineflowandlossesIMGWiDkflP
在matlab中运行修改后的程序，得到潮流计算的结果:
Totalloss37.990MW108.014MVar（无功补偿为19Mvar>
Totalloss17.598MW22.245MVar<
由以上的潮流计算结果看出，减少电力系统的回路数，系统的潮流发生了变化，系统总网损发生了很大的变化，总的有功损耗减小，可是总的无功损耗确增大了，即电力系统的稳定性变差了。
WHF4OmOgAw
4、实验分析
在matlab中运行以上程序，得到各个节点的电压幅值和相位，各条支路的潮流和网损。
此程序的运行结果和实验程序给出的结果是一致的。
说明程序无误，但在精度上有微小差异，这主要是和导纳矩阵的精度以及现实精度有关。
通过改变支路中的数据，可以减小系统的回路数，其中采用的是将部分支路的电阻修改为无穷大，也相当于相应断路器断开。
aDFdk6hhPd
5、实验心得
本程序分模块进行，先是排序，再是求导纳矩阵，然后求雅阁比，再进行迭代运算。
程序本身很简洁明了，运行的时候只要在matlab中输入main就行了，然后打开bus和line所在的m文件中，结果就自动存在rustl文件中。
通过编写牛顿拉夫逊法matlab潮流计算程序复习了潮流计算知识，也实现了计算机算法。
ozElQQLi4T
电压是电能质量的重要指标之一，网损是电力企业的一项重要综合性技术经济指标。
长期以来电力系统网络损耗问题比较突出，而无功补偿是降低线损的有效手段。
随着电力系统负荷的增加，对无功功率的需求也日益增加。
在电网中的适当位置装设无功补偿装置成为满足电网无功需求的必要手段。
提高电压质量、降损节能是供电企业当前的工作重点之一。
对电网给予合理的无功补偿，能够有效地维持系统的电压水平，降低有功网损，提高网络输送容量满足负荷的无功需求。
节点负荷
升级会员 