基于Matlab计算程序的电力系统运行分析.docx

1、基于Matlab计算程序的电力系统运行分析 课 程 设 计课程名称： 电力系统分析 设计题目：基于Matlab计算程序的电力系统运行分析学 院： 电力工程学院 专 业： 电气工程及其自动化 年 级： 2011级 学生姓名： 指导教师： 、 、 日 期： 教 务 处 制目录 前 言1第1章 参数计算 2 第一节 原始数据 2 第二节 电网模型的建立 3第二章 潮流计算 6 第一节 系统参数的设置 6 第二节 程序的调试 7第三节 网损和潮流分析 13第三章 短路故障的分析计算15第一节 序网分析15 第二节 短路计算18心得体会35参考文献36 前 言电力系统潮流计算是电力系统分析中的一种最基本

2、的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。

3、其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。第1章 参数计算第一节 原始数据一、 目标电网接线图图11二、系统参数表1. 线路参数表线路编号线路型号线路长度（km）线路电阻/km线路正序电抗/km线路容纳之半S/km4-5LGJ-240/301130.0470.41.784-6LGJ-120/701200.0741.475-7LGJ-120/251650.0791.606-9LGJ-95/551660.0921.807-8LGJ-240/309

4、20.0471.788-9LGJ-240/301220.0471.78说明：线路零序电抗为正序电抗3倍。表2. 变压器参数表线路编号变压器型号变压器变比（kV）短路电压百分数(%)2-7SSPL-22000024232.5%/2010.433-9SSPL-12000024232.5%/155.811-4SSPL-24000024232.5%/17.511.42说明：变压器零序电抗与正序电抗相等，且均为/Y0接法。 表3. 发电机参数表1发电机额定功率MW额定电压kV额定功率因数120016.50.852180180.85310013.80.85表4. 发电机参数表2发电机母线名（S）（S）（S

5、）1147.2800.320.130.210.218.962212.8001.930.261.870.436.000.535336.0201.510.211.450.298.590.60表5. 负荷数据表节点号有功负荷（MW）无功负荷（MVA）51355061003088035第二节 电网模型的建立本设计中，采用精确计算算法，选取=100MVA， =220KV，将所有支路的参数都折算到220KV电压等级侧；一、计算公式 （1）线路参数的计算： （2）变压器参数的计算：（3）发电机参数的计算：（暂态分析时，只用到发电机的暂态电抗来代替其次暂态电抗，故只求出暂态电抗）（4）负荷节点的计算二、使用编

6、程的结果如下：线路电抗X= 0.0934 0.0992 0.1364 0.1372 0.0760 0.1008线路电阻R= 0.0109 0.0183 0.0269 0.0316 0.0089 0.0118线路正序、负序阻抗Z1=Z2= Columns 1 through 5 0.0109 + 0.0937i 0.0183 + 0.0992i 0.0269 + 0.1364i 0.0316 + 0.1372i 0.0089 + 0.0760i Column 6 0.0118 + 0.1008i 线路容纳之半B/2= Columns 1 through 5 0 + 0.0974i 0 + 0.0

7、854i 0 + 0.1280i 0 + 0.1450i 0 + 0.0793i Column 6 0 + 0.1050i线路零序阻抗Z0: Columns 1 through 5 0.0109 + 0.2802i 0.0183 + 0.2975i 0.0269 + 0.4091i 0.0316 + 0.4112i 0.0089 + 0.2281i Column 6 0.0118 + 0.3024i 变压器正、负序电抗值Xt： 0 + 0.0576i 0 + 0.0574i 0 + 0.0586i变压器零序电抗值Xt0： 0 + 0.0576i 0 + 0.0574i 0 + 0.0586i发

8、电机暂态电抗Xdg= 0 + 0.0476i 0 + 0.0803i 0 + 0.1103i负荷节点阻抗Zl= 0.6514 - 0.2413i 0.9174- 0.2752i 1.049 - 0.4590i以上结果均采用程序计算。三、等值电路模型该等值电路图全部数值以标幺值表示，具体图12所示：图12 第二章 潮流计算第一节 系统参数的设置设计中要求所有结点电压不得低于1.0p.u.，也不得高于1.05p.u.，若电压不符合该条件，可采取下面的方法进行调压：（1） 改变发电机的机端电压（2） 改变变压器的变比（即改变分接头）（3） 改变发电机的出力（4） 在电压不符合要求的结点处增加无功补偿

9、调压方式应属于逆调压。结点的分类：根据电力系统中各结点性质的不同，将结点分为三类：PQ结点、PV结点和平衡结点，在潮流计算中，大部分结点属于PQ结点，小部分结点属于PV结点，一般只设一个平衡结点。对于平衡结点，给定其电压的幅值和相位，整个系统的功率平衡由这一点承担。本设计中，选1号节点为平衡节点；2、3号节点为P、U节点；4、5、6、7、8、9号结点为P、Q节点。设计中，节点数:n=9，支路数:nl=9，平衡母线节点号:isb=1，误差精度:pr=0.0001。由支路参数形成的矩阵:B1= 1 4 0.0576i 0 1.0000 0; 2 7 0.0574i 0 1.0000 0; 3 9

10、0.0586i 0 1.0000 0; 4 5 0.0109+0.0934i 2*0.0974i 1.0000 0; 4 6 0.0183+0.0992i 2*0.0854i 1.0000 0; 5 7 0.0269+0.1364i 2*0.1280i 1.0000 0; 6 9 0.0316+0.1372i 2*0.1450i 1.0000 0; 7 8 0.0089+0.0760i 2*0.0793i 1.0000 0; 8 9 0.0118+0.1008i 2*0.1050i 1.0000 0; ;%支路参数矩阵由各节点参数形成的矩阵:B2=0 0 1.0000 1.0000 0 1.0

11、000;1.8000 0 1.0000 1.0000 0 3.0000;1.0000 0 1.0000 1.0000 0 3.0000;0 0 1.0000 0 0 2.0000;0 1.3500+0.5000i 1.0000 0 0 2.0000;0 0.0000+0.3000i 1.0000 0 0 2.0000;0 0 1.0000 0 0 2.0000;0 0.8000+0.3500i 1.0000 0 0 2.0000;0 0 1.0000 0 0 2.0000;%节点参数矩阵由节点号及其对地阻抗形成的矩阵:X=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0

12、;第二节 程序的调试一、未调试前，原始参数运行结果如下：注意：输入“var=1”时将选用牛顿-拉夫逊法来进行潮流计算；输入“var=其他数字时将选用PQ分解法来进行潮流计算”var=1迭代次数 4没有达到精确要求的个数 14 16 8 0各节点的实际电压标么值E为(节点号从小到大排列):Columns 1 through 4 1.0000 0.9759 + 0.2183i 0.9906 + 0.1368i 0.9840 - 0.0240iColumns 5 through 8 0.9498 - 0.0429i 0.9689 - 0.0388i 0.9863 + 0.1148i 0.9814 +

13、 0.0641iColumn 9 0.9969 + 0.0785i各节点的电压大小V为(节点号从小到大排列):Columns 1 through 7 1.0000 1.0000 1.0000 0.9843 0.9508 0.9697 0.9929Columns 8 through 9 0.9835 1.0000各节点的电压相角O为(节点号从小到大排列):Columns 1 through 7 0 12.6101 7.8603 -1.3973 -2.5839 -2.2909 6.6372Columns 8 through 9 3.7377 4.5008各节点的功率S为(节点号从小到大排列):Co

14、lumns 1 through 4 0.4167 + 0.2781i 1.8000 + 0.2175i 1.0000 + 0.0298i 0.0000 - 0.0000iColumns 5 through 8 -1.3500 - 0.5000i -1.0000 - 0.3000i -0.0000 + 0.0000i -0.8000 - 0.3500iColumn 9 0.0000 - 0.0000i各条支路的首端功率Si为(顺序同您输入B1时一样): 0.4167 + 0.2781i 1.8000 + 0.2175i 1.0000 + 0.0298i 0.2456 + 0.2320i 0.17

15、11 + 0.0316i -1.1063 - 0.1017i -0.8297 - 0.1097i 0.6573 - 0.0157i -0.1467 - 0.2444i各条支路的末端功率Sj为(顺序同您输入B1时一样): -0.4167 - 0.2636i -1.8000 - 0.0288i -1.0000 + 0.0289i -0.2437 - 0.3983i -0.1703 - 0.1903i 1.1427 + 0.0445i 0.8528 - 0.0711i -0.6533 - 0.1056i 0.1472 + 0.0422i各条支路的功率损耗DS为(顺序同您输入B1时一样): 0 + 0

16、.0145i 0.0000 + 0.1887i 0.0000 + 0.0587i 0.0019 - 0.1663i 0.0008 - 0.1587i 0.0364 - 0.0572i 0.0232 - 0.1808i 0.0039 - 0.1213i 0.0005 - 0.2022i以下是每次迭代后各节点的电压值(如图所示) 图21 由运行结果可知，节点4、5、6、7、8电压均不满足要求。故需进行调试，以期各结点电压均满足要求。下面用N-L法进行潮流计算。二、采用NL法进行潮流的计算和分析（1）第一次调试将变压器1的变比初始值由1.0000改为1.0250，则修改结果如下：B1=1.0000

17、4.0000 0.0576i 0 1.0250 0; 2.0000 7.0000 0.0574i 0 1.0000 0;运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0022 0.9638 0.9823 0.9968 0.9875 1.0038 节点5、6、7、8都不满足要在1.00001.0500范围内的要求，改变变压器1的变比不满足要求。（2）第二次调试将变压器2的变比初始值由1.0000改为1.0250，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0250 0; 2.0000 7.0000 0.057

18、4i 0 1.0250 0; 运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0066 0.9745 0.9871 1.0147 1.0001 1.0085节点5、6的值，虽然第二次调试仍然不满足要在1.00001.0500范围内的要求，但电压有所升高。（3）第三次调试将变压器3的变比初始值由1.0000改为1.0250，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0250 0; 2.0000 7.0000 0.0574i 0 1.0250 0; 3.0000 9.0000 0.0586i 0 1.0250

19、0; 运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0109 0.9795 0.9977 1.0196 1.0109 1.0265 节点5、6的值，虽然第三次调试仍然不满足要在1.00001.0500范围内的要求，但电压有所升高。（4）第四次调试将变压器1的变比初始值由1.0250改为1.0500，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0500 0; 2.0000 7.0000 0.0574i 0 1.0250 0; 3.0000 9.0000 0.0586i 0 1.0250 0; 运行结果如下：各

20、节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0287 0.9924 1.0103 1.0236 1.0150 1.0305 节点5的值，不满足要在1.00001.0500范围内的要求，且变压器T1的变比以调至最大。应调节与节点5临近的变压器T2的变比。（5）第五次调试将变压器2的变比初始值由1.0250改为1.0500，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0500 0; 2.0000 7.0000 0.0574i 0 1.0500 0; 3.0000 9.0000 0.0586i 0 1.0250 0;运行结果如

21、下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0332 1.0031 1.0152 1.0413 1.0276 1.0353所有节点的值，都满足在1.00001.0500范围内的要求。NL法 程序运行结果如下：注意：输入“var=1”时将选用牛顿-拉夫逊法来进行潮流计算；输入“var=其他数字时将选用PQ分解法来进行潮流计算”var=1迭代次数 4没有达到精确要求的个数 14 16 10 0各节点的实际电压标么值E为(节点号从小到大排列):Columns 1 through 4 1.0000 0.9783 + 0.2074i 0.9914 + 0

22、.1307i 1.0329 - 0.0248iColumns 5 through 8 1.0022 - 0.0423i 1.0145 - 0.0384i 1.0356 + 0.1087i 1.0258 + 0.0609iColumn 9 1.0325 + 0.0755i各节点的电压大小V为(节点号从小到大排列):Columns 1 through 7 1.0000 1.0000 1.0000 1.0332 1.0031 1.0152 1.0413Columns 8 through 9 1.0276 1.0353各节点的电压相角O为(节点号从小到大排列):Columns 1 through 7

23、0 11.9693 7.5084 -1.3780 -2.4143 -2.1662 5.9892Columns 8 through 9 3.3961 4.1824各节点的功率S为(节点号从小到大排列):Columns 1 through 4 0.4108 + 0.2829i 1.8000 + 0.2384i 1.0000 - 0.1423i -0.0000 + 0.0000iColumns 5 through 8 -1.3500 - 0.5000i -1.0000 - 0.3000i -0.0000+0.0000i -0.8000 - 0.3500iColumn 9 -0.0000 - 0.00

24、00i各条支路的首端功率Si为(顺序同您输入B1时一样): 0.4108 + 0.2829i 1.8000 + 0.2384i 1.0000 - 0.1423i 0.2366 + 0.2034i 0.1743 + 0.0651i -1.1150 - 0.1078i -0.8267 - 0.0608i 0.6518 + 0.0403i -0.1519 - 0.1709i各条支路的末端功率Sj为(顺序同您输入B1时一样): -0.4108 - 0.2685i -1.8000 - 0.0491i -1.0000 + 0.2021i -0.2350 - 0.3922i -0.1733 - 0.2392

25、i 1.1482 + 0.0088i 0.8479 - 0.1521i -0.6481 - 0.1791i 0.1521 - 0.0500i各条支路的功率损耗DS为(顺序同您输入B1时一样): -0.0000 + 0.0143i -0.0000 + 0.1892i 0.0000 + 0.0598i 0.0015 - 0.1888i 0.0009 - 0.1741i 0.0332 - 0.0990i 0.0212 - 0.2128i 0.0036 - 0.1388i 0.0003 - 0.2209i以下是每次迭代后各节点的电压值(如图所示)图22三、对运行结果的分析：（1）、为什么在用计算机对某

26、网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的？哪些可调？哪些不可调？ 答：潮流计算时功率方程是非线性，多元的具有多解。初始条件给定后得到的结果不一定能满足约束条件要求，要进行调整初值后才能满足。其约束条件有：，。负荷的PQ量为扰动变量，发电机的PQ为控制变量，各节点的V为状态变量。扰动变量是不可控变量，因而也是不可调节的，状态变量是控制变量的函数，因而状态变量和控制变量是可以调节的。所以，计算机对某网络初次进行潮流计算时往往是要调潮流的。（2）、牛顿拉夫逊法与PQ分解法有哪些联系？有哪些区别？二者的计算性能如何？答：在

27、潮流计算中，牛顿拉夫逊法的实质 是一种逐步线性化的方法，它有很好的收敛性，但要求有合适的初值。 PQ分解法来源于N-L法极坐标形式，又充分考虑在交流高压电网中，变压器和线路的电抗远远大于电阻，且系统中母线有功功率的变化主要受电压相位的影响，而无功功率的变化主要受母线电压幅值的影响，在此基础之上对H、L阵进行简化。由于这种简化只涉及修正方程的系数矩阵，并为改变结点功率平衡方程和收敛判据，从NL法和PQ分解法程序运行结果可以看出不会降低计算结果的精度。由图22和图23比较可知，PQ分解法与N-L法相比大大提高了计算速度，但迭代次数却增多；PQ分解法为恒速率收敛，比起N-L法，其收敛速度变慢。第三节

28、 网损和潮流分析一 、选取NL法的数据来分析网损表21调整前后网络损耗比较节点编号1-42-73-94-54-65-76-97-88-9调整前0.0145i0.1887i0.0587i0.0019 - 0.1663i0.0008 - 0.1587i0.0364- 0.0572i0.0232 - 0.1808i0.0039 - 0.1213i0.0005 - 0.2022i牛拉法调整后0.0143i0.1892i0.0598i0.0015 - 0.1888i0.0009 - 0.1741i0.0332 - 0.0990i 0.0212 - 0.2128i0.0036 - 0.1388i0.0003 - 0.2209i调节前有功损耗：P=（0+0+0+0.0019+0.0008+0.03