matlab电力系统分析报告课程设计.doc

课 程 设 计课程名称： 电力系统分析 设计题目：基于Matlab计算程序的电力系统运营分析学 院： 电力工程学院 专 业： 电气工程自动化 年 级： 2023 学生姓名： 胡健（） 指导教师： 单节杉、宋琪、谢静 日 期： 2023年1月12日 教 务 处 制目录前 言 1第一章 参数计算 2 一、目的电网接线图 2 二、电网模型的建立 3第二章 潮流计算 7 一系统参数的设立 7 二 程序的调试 8 三、对运营结果的分析 17第三章 短路故障的分析计算 19 一、三相短路 19 二、不对称短路 22 三、由上面表对运营结果的分析及在短路中的一些问题 33心得体会 35参考文献 36前 言电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运营状态的计算。潮流计算的目的是求取电力系统在给定运营状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足规定，功率的分布和分派是否合理以及功率损耗等。对现有电力系统的运营和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中也许发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的重要因素是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运营时，要采用适当的措施减少短路故障的发生概率。短路计算可认为设备的选择提供原始数据。第一章 参数计算一、 目的电网接线图系统参数表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/30920.0471.788-9LGJ-240/301220.0471.78说明：线路零序电抗为正序电抗3倍。表2. 变压器参数表线路编号变压器型号变压器变比（kV）短路电压百分数(%)2-7SSPL-22023024232.5%/2010.433-9SSPL-12023024232.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）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、系统参数的计算（1）线路参数计算公式如下： 各条线路参数的结果：4-5： 4-6： 5-7： 6-9： 7-8： 8-9： （2）变压器参数的计算：（3）发电机参数的计算：（暂态分析时，只用到发电机的暂态电抗来代替另一方面暂态电抗，故只求出暂态电抗）（4）负荷节点的计算 2系统等值电路图的绘制根据以上计算结果，得到系统等值电路图如下： 第二章 潮流计算一系统参数的设立设计中规定所有结点电压不得低于1.0p.u.，也不得高于1.05p.u.，若电压不符合该条件，可采用下面的方法进行调压：（1） 改变发电机的机端电压（2） 改变变压器的变比（即改变分接头）（3） 改变发电机的出力（4） 在电压不符合规定的结点处增长无功补偿调压方式应属于逆调压。结点的分类：根据电力系统中各结点性质的不同，将结点分为三类：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.00001。由支路参数形成的矩阵:B1=1.0000 4.0000 0.0555i 0 1.0000 0; 2.0000 7.0000 0.0550i 0 1.0000 0;3.0000 9.0000 0.0605i 0 1.0000 0;4.0000 5.0000 0.0130+0.1107i 0.1154i 1.0000 0; 4.0000 6.0000 0.0229+0.1240i 0.1067i 1.0000 0;5.0000 7.0000 0.0144+0.0727i 0.0681i 1.0000 0;6.0000 9.0000 0.0240+0.1041i 0.1098i 1.0000 0;7.0000 8.0000 0.0087+0.0744i 0.0775i 1.0000 0;8.0000 9.0000 0.0120+0.1025i 0.1068i 1.0000 0;%支路参数矩阵由各节点参数形成的矩阵:B2=2.2023+1.2395i 0 1.0500 1.0000 0 1.0000;1.8000+1.1155i 0 1.0000 1.0000 0 3.0000;1.0000+0.6197i 0 1.0000 1.0000 0 3.0000;0 0 1.0000 0 0 2.0000;0 1.2500+0.5000i 1.0000 0 0 2.0000;0 0.9000+0.3000i 1.0000 0 0 2.0000;0 0 1.0000 0 0 2.0000;0 1.0000+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;二 程序的调试1. 未调试前，原始参数运营结果如下：输入“var=1”时将选用牛顿-拉夫逊法来进行潮流计算；输入“var=其他数字时将选用PQ分解法来进行潮流计算”var=2迭代次数 8每次没有到达精度规定的有功功率个数为 5 8 8 8 7 4 2 0每次没有达成精度规定的无功功率个数为 6 6 6 6 6 5 2 0各节点的电压标幺值E为（节点号从小到大排）: Columns 1 through 7 1.0000 0.9892 + 0.1467i 0.9964 + 0.0852i 0.9813 - 0.0217i 0.9464 - 0.0355i 0.9577 - 0.0466i 0.9769 + 0.0448i Columns 8 through 9 0.9631 - 0.0074i 0.9855 + 0.0236i各节点的电压U大小(节点号从小到大排)为: 1.0000 1.0000 1.0000 0.9816 0.9470 0.9588 0.9779 0.9631 0.9858低于电压下限的节点为: 5低于下限的电压值为： 0.9470各节点的电压相角O（节点号从小到大排）: 0 8.4374 4.8883 -1.2657 -2.1475 -2.7877 2.6269 -0.4381 1.3696各节点的功率S（节点号从小到大排）为: Columns 1 through 7 0.3907 + 0.3365i 1.8000 + 0.4934i 1.0000 + 0.2662i -0.0000 + 0.0001i -1.2500 - 0.5000i -0.9000 - 0.3000i 0.0000 - 0.0000i Columns 8 through 9 -1.0000 - 0.3500i 0.0000 - 0.0000i各支路的首段功率Si(顺序同您输入B1时同样）为: 0.3907 + 0.3365i 1.8000 + 0.4934i 1.0000 + 0.2662i 0.1630 + 0.2325i 0.2276 + 0.0894i -1.0884 - 0.1728i -0.6741 - 0.1194i 0.6922 + 0.0944i -0.3123 - 0.2212i各条支路的末端功率Sj（顺序同您输入B1时同样）为: -0.3907 - 0.3217i -1.8000 - 0.3018i -1.0000 - 0.2023i -0.1616 - 0.3272i -0.2259 - 0.1806i 1.1078 + 0.2073i 0.6861 + 0.0676i -0.6877 - 0.1288i 0.3139 + 0.1338i各条支路的功率损耗DS（顺序同您输入B1时同样）为: 0 + 0.0148i 0 + 0.1916i -0.0000 + 0.0648i 0.0015 - 0.0948i 0.0017 - 0.0912i 0.0193 + 0.0346i 0.0120 - 0.0518i 0.0045 - 0.0344i 0.0016 - 0.0874i以下是每次迭代后各节点的电压值(如图所示)（用PQ分解法及N-L法运营的结果同样，故只给出N-L法运营的结果） 由运营结果可知，节点4、5、6、7、8、9电压均不满足规定。故需进行调试，以期各结点电压均满足规定。一下分别用N-L法及PQ分解法进行潮流计算。2.采用NL法进行潮流的计算和分析。1）第一次调试将2号发电机电压初值由1.0000升高至1.0200，则修改结果如下：运营结果如下： 如图所示：节点5、6、7、8、9都不满足要在1.00001.0500范围内的规定，改变发电机机端电压的影响不大。2）第二次调试 将1号变压器变比初值由1.0000改至1.0250，则修改结果如下： 运营结果如下： 如图所示：节点5、6、7、8、9的值，虽然较第一次调试仍然不满足要在1.00001.0500范围内的规定，但电压有所升高。 将1号变压器变比由1.0250改至1.0500，则修改结果如下： 运营结果如下： 如图所示：节点4、5、7、8的值，不满足要在1.00001.0500范围内的规定，且节点4的值超过了范围上限1.0500。则应当对与节点4临近的变压器T1的变比进行再次调节；同时调节与节点5、7、8临近的变压器T2的变比。3）第三次调试将1号变压器变比由1.0500重新改至1.0250，同时把2号变压器变比由1.0000改至1.0250，则修改结果如下：运营结果如下：如图所示：节点5、6、8、的值，不满足要在1.00001.0500范围内的规定。4. 第四次调试将2号变压器变比由1.0250改至1.0500，则修改结果如下：运营结果如下：（满足条件）NL法 程序运营结果如下：3.采用PQ分解法法进行潮流的计算和分析。1）第一次调试将2号发电机电压初值由1.0000升高至1.0200，则修改结果如下：运营结果如下：如图所示：节点4、5、6、7、8、9都不满足要在1.00001.0500范围内的规定。2）第二次调试 将1号变压器变比初值由1.0000改至1.0250，则修改结果如下： 运营结果如下： 如图所示：节点4、5、6、7、8、9的值，虽然较第一次调试仍然不满足要在1.00001.0500范围内的规定，但电压有所升高。 则，将1号变压器变比由1.0250改至1.0500，则修改结果如下： 运营结果如下： 如图所示：节点5、6、7、8、9的值，不满足要在1.00001.0500范围内的规定。3）第三次调试 将2号变压器变比初值由1.0000改至1.0250，则修改结果如下： 运营结果如下： 如图所示：节点5、6、8、9的值，仍然不满足要在1.00001.0500范围内的规定，但电压有所升高。 则，将2号变压器变比由1.0250改至1.0500，则修改结果如下： 运营结果如下： 如图所示：节点5、6、8的值，仍然不满足要在1.00001.0500范围内的规定。4）第四次调试将3号变压器变比由1.0000改至1.0250，则修改结果如下：运营结果如下：如图所示：只有节点5的值，不满足要在1.00001.0500范围内的规定，应当对与节点5临近的变压器T2的变比进行再次调节。5）第五次调试将2号变压器变比由1.0500改至1.0750，则修改结果如下：运营结果如下：（满足条件）PQ分解法 程序运营结果如下：三、对运营结果的分析：1、为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的？哪些可调？哪些不可调？ 答：潮流计算时功率方程是非线性，多元的具有多解。初始条件给定后得到的结果不一定能满足约束条件规定，要进行调整初值后才干满足。其约束条件有：，。负荷的PQ量为扰动变量，发电机的PQ为控制变量，各节点的V为状态变量。扰动变量是不可控变量，因而也是不可调节的，状态变量是控制变量的函数，因而状态变量和控制变量是可以调节的。所以，计算机对某网络初次进行潮流计算时往往是要调潮流的。2、潮流控制的重要手段有哪些？ 答：潮流控制的重要手段有：（1）改变发电机的机端电压（2）改变变压器的变比（即改变分接头）（3）改变发电机的出力（4）在电压不符合规定的节点处增长无功补偿3、 牛顿拉夫逊法与PQ分解法有哪些联系？有哪些区别？两者的计算性能如何？答：（1）联系：它们采用相同的数学模型和收敛判据。当电路的电抗远大于电阻时，可以简化牛顿拉夫逊极坐标的修正方程的系数矩阵得到PQ分解法，且简化后并未改变节点功率平衡方程和收敛判据，因而不会减少计算结果的精度。（2）区别：P-Q分解法的修正方程结构和牛顿拉夫逊的结构不同。pq分解法由于雅可比矩阵常数化，计算过程中减少了很大的计算量，并且有功和电压幅值，无功和电压相角的完全割裂也大大的对矩阵降维数，减少了一半的计算量，但是他雅克比矩阵常数化是经验值，丧失了一部分稳定收敛的特性，并且当支路电阻与电抗比值较大的时候收敛性也特别差，甚至不收敛（3）P-Q法按几何级数收敛，牛顿拉夫逊法按平方收敛。PQ分解法把节点功率表达为电压向量的极坐标方程式，抓住重要矛盾，把有功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开进行。它密切地结合了电力系统的固有特点，无论是内存占用量还是计算速度方面都比牛顿-拉夫逊法有了较大的改善。4、 选取PQ分解法的数据来分析减少网损的方法：支路未调整前：调整后：支路首端功率支路末端功率支路功率损耗支路首端功率支路末端功率支路功率损耗1-40.3907 + 0.3365i-0.3907 - 0.3217i0+ 0.0148i0.3866+0.3017i-0.3866-0.2883i0+ 0.0133i2-71.8000 + 0.4934i-1.8000 - 0.3018i0 + 0.1916i1.8000+0.6750i-1.8000-0.4717i0 + 0.2033i4-51.0000 + 0.2662i-1.0000 - 0.2023i-0.0000 + 0.0648i1.0000-0.0424i-1.0000+0.0182i0 + 0.0606i3-90.1630 + 0.2325i-0.1616 - 0.3272i0.0015 - 0.0948i0.1553+0.1536i-0.1545-0.2664i0.0009-0.1128i4-60.2276 + 0.0894i-0.2259 - 0.1806i0.0017 - 0.0912i0.2314+0.1342i-0.2294-0.2345i0.0019-0.1003i5-7-1.0884 - 0.1728i1.1078 + 0.2073i0.0193 + 0.0346i-1.0955-0.2336i1.1131+ 0.2509i0.0176+0.0173i6-9-0.6741 - 0.1194i0.6861 + 0.0676i0.0120 - 0.0518i-0.6706-0.0654i0.6813- 0.0013i0.0107-0.0667i7-80.6922 + 0.0944i-0.6877 - 0.1288i0.0045 - 0.0344i0.6869+0.2209i-0.6825-0.2657i0.0043-0.0449i8-9-0.3123 - 0.2212i0.3139 + 0.1338i0.0016 - 0.0874i-0.3175-0.0842i0.3187- 0.0170i0.0012-0.1012i（1）提高机端电压电压和节点电压一定可以使有功损耗减少，但是对于无功损耗来说为正的是可以减少的，为负的则是提高了；（2）此外适当提高负荷的功率因数、改变电力网的运营方式，对原有电网进行技术改造都可以减少网损。5、 发电机节点的注入无功为负值说明了什么？答：由于线路无功潮流最有也许的流向由电压的幅值大小决定：由幅值高的节点流向幅值低的节点。由此看出发电机的电压小于节点电压而无功功率的方向是从高电压到低电压，所以发电机的注入无功为负值。6、负荷功率因数对系统潮流有什么影响？答：负荷功率因数减少，无功功率就会增大，其输电线路的总电流就会相应增大，从而会导致电压损耗的升高，从而会改变无功功率潮流的大小，严重时甚至会改变方向；反之亦然。三 绘制潮流分布图第三章 短路故障的分析计算一、三相短路1、三相短路电流的计算原理:运用结点阻抗矩阵和导纳矩阵都可以计算短路电流，其算法有所不同。运用结点阻抗阵时，只要形成了阻抗阵，计算网络中任意一点的对称短路电流和网络中电流、电压的分布非常方便，计算工作量小，但是，形成阻抗阵的工作量大，网络变化时的修改也比较麻烦，并且结点阻抗矩阵是满阵，需要计算机存储量较大。运用结点导纳矩阵计算短路电流，实质是先用它计算出与指定的短路点有关的结点阻抗矩阵的相应列的元素，然后运用公式（为连接节点i和j的支路阻抗）进行短路电流的计算。然而，导纳阵是对称、稀疏阵，极易形成，且网络结构变化时也易于修改。程序运营环节及对变量的解释： （1） 请输入短路点的数目NF（2） 请输入节点数n（3） 请输入支路数nl（4） 输入各支路参数矩阵B 矩阵B的每行是由下列参数构成的： 某支路的首端号P； 末端号Q，且PQ； 支路的阻抗（R+jX）; 支路的对地电纳； 支路的变比K； 折算到哪一侧的标志（假如支路的首端P处在高压侧请输入“1”，否则请输入“0”）。（5） 输入由短路点号，短路点阻抗形成的矩阵D（6） 请输入由各节点的初电压标幺值形成的列矩阵V0（7） 形成节点阻抗矩阵Z（8） 求短路点电流（9） 求网络中各节点的电压（10）求网络中各支路电流参数矩阵B = 1.0000 4.0000 0 + 0.0682i 0 1.0000 0 2.0000 7.0000 0 + 0.0824i 0 1.0000 0 3.0000 9.0000 0 + 0.1452i 0 1.0000 0 4.0000 5.0000 0.0119 + 0.1013i 1.0000 1.0000 0 4.0000 6.0000 0.0210 + 0.1134i 2.0000 1.0000 0 5.0000 7.0000 0.0131 + 0.0665i 3.0000 1.0000 0 6.0000 9.0000 0.0219 + 0.0953i 4.0000 1.0000 0 7.0000 8.0000 0.0080 + 0.0681i 4.0000 1.0000 0 8.0000 9.0000 0.0110 + 0.0938i 5.0000 1.0000 0 ; 2、 三相短路的等值网络图：5节点、6节点发生三相短路的程序结果：1、 5节点发生三相短路： 1）短路点的电流：标幺值： 0.4481 - 6.1435i有名值（kA）： 0.1125 - 1.5421i 2）各节点的电压：电压标幺电压有名值（kV）1118.37521213115.7541.0327237.4985 -2.2204e-016 -2.7756e-017i061.0049232.73770.0328 - 0.0000i6.693+ 0.0000i80.0094 - 0.0000i3542+ 0.0000i90.0166 - 0.0000i8.142+ 0.0000i 3）各支路电流：电流标幺电流有名值（kA）140 - 0.2416i0 - 0.0606i270 -11.7656i0 - 2.9534i390 - 6.7755i0 - 1.7008i451.1813 -10.0557i0.2965 - 2.5242i460.0439 - 0.2370i0.0110 - 0.0595i57-0.0935 + 0.4748i-0.0235 + 0.1192i692.2636 - 9.8502i0.5682 - 2.4726i780.0398 - 0.3389i0.0100 - 0.0851i89-0.0089 + 0.0757i-0.0022 + 0.0190i2、6节点发生三相短路 1）短路点的电流：标幺值： 0.6773 - 5.6526i有名值（kA）：0.1700 - 1.4189i 2）各节点的电压：电压标幺电压有名值（kV）1118.37521213115.7540.7485 - 0.0375i166.75-9.026i51.0077230.22960071.0405236.99281.0171233.77291.0243238.28 3 )各支路电流：电流标幺电流有名值（kA）140.5851 - 5.1245i0.1469 - 1.2864i270 - 0.3895i0 - 0.0978i390 - 0.0047i0 - 0.0012i45-0.7794 + 3.1135i-0.1956 + 0.7815i460.7212 - 5.8897i0.1810 - 1.4784i57-0.0935 + 0.4748i-0.0235 + 0.1192i69-2.3460 +10.2090i-0.5889 + 2.5627i780.0398 - 0.3389i0.0100 - 0.0851i89-0.0089 + 0.0757i-0.0022 + 0.0190i3、对运营结果的分析:由上面的结果可知，当发生三相短路时，其短路点的电压为0，短路点的电流为最大，这与理论都是相符的。二、不对称短路1、短路计算原理：简朴不对称故障（涉及横向和纵向故障）与对称故障的计算环节是一致的，一方面算出故障口的电流，接着算出网络中各个结点的电压，由结点电压即可拟定支路电流，所不同的是，要分别按三个序进行。（1）系统三序等值网络图如下：不对称短路计算的三序等值网络图（将发电机与变压器支路合并为一条支路）：1）系统的正序等值网络图：2）系统的负序等值网络图：3） 系统的零序等值网络图：2、5节点和6节点发生不对称短路的程序结果（1）5节点发生不对称短路的程序结果：1）故障点和各支路各序电流标幺值：正序电流标幺值负序电流标幺值零序电流标幺值单相接地短路短路点(5号点)0.0795 - 1.4956i0.0795 - 1.4956i0.0795 - 1.4956i两相短路0.2240 - 3.0717i-0.2240 + 3.0717i0两相接地短路0.2486 - 3.6605i-0.1994 + 2.4830i-0.0492 + 1.1776i 单相接地短路140 - 0.2416i0027-0.0545 - 3.1550i-0.0545 - 2.7656i0.1976 -10.5853i39-0.0325 - 1.6507i-0.0325 - 1.6460i0.1885 -10.1025i450.2618 - 2.6343i0.2332 - 2.3909i0.0983 - 1.6993i460.0439 - 0.2370i0057-0.0935 + 0.4748i00690.4685 - 2.2593i0.5129 - 2.4527i 0.1716 - 1.7984i780.0398 - 0.3389i0079-0.0089 + 0.0757i00两相短路140 - 0.2416i00270 - 6.0776i0 + 5.6881i0390 - 3.3901i0 + 3.3854i0450.6049 - 5.1496i-0.5763 + 4.9061i0460.0439 - 0.2370i0057-0.0935 + 0.4748i00691.1096 - 4.8284i-1.1540 + 5.0218i0780.0398 - 0.3389i0079-0.0089 + 0.0757i00两相接地短路140 - 0.2416i0027-0.0338 - 7.1654i-0.0338 + 4.6003i-0.0612 + 8.3310i39-0.0201 - 4.0376i-0.0201 + 2.7380i-0.0584 + 7.9510i450.6860 - 6.0913i-0.4952 + 3.9644i-0.0622 + 1.3380i460.0439 - 0.2370i0057 -0.0935 + 0.4748i00691.3005 - 5.7957i-0.9631 + 4.0545i-0.1190 + 1.4167i780.0398 - 0.3389i0079-0.0089 + 0.0757i002）故障点和各支路各序电压标么值：节点单相接地故障两相短路故障两相接地故障各节点的正序电压标幺值为11.0000 1.00001.000021.0000 1.00001.050031.0000 1.00001.048041.0327 1.03271.032750.7777 - 0.0026i0.5038 + 0.0000i0.4164 - 0.0015i61.0049 1.00491.004970.8105 - 0.0026i0.5366 + 0.0000i0.4492 - 0.0015i80.7871 - 0.0026i0.5132 + 0.0000i0.4258 - 0.0015i90.7943