电力系统潮流计算课程设计(共32页).doc

《电力系统潮流计算课程设计(共32页).doc》由会员分享，可在线阅读，更多相关《电力系统潮流计算课程设计(共32页).doc（32页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上一、问题重述1.1 课程设计要求1、在读懂程序的基础上画出潮流计算基本流程图2、通过输入数据，进行潮流计算输出结果3、对不同的负荷变化，分析潮流分布，写出分析说明。4、对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控制的参数变化。5、打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。1.2 课程设计题目1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。变电所1变电所2变电所3变电所435kV母线10kV母线35kV母线10kV母线一次侧电压220kV一次侧电压220kV母线1母线3母线2线路长为100km线路长为80

2、km线路长为60km线路长为80km线路长为100km2*QFQ-50-22*QFS-50-22*TQN-100-22*TQN-100-2。线路长为80km。电厂一电厂二2、发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为（400MW），母线3为机压母线，机压母线上装机容量为（100MW），最大负荷和最小负荷分别为50MW和30MW；发电厂二总装机容量为（200MW）。3、变电所资料：（1）变电所1、2、3、4低压母线的电压等级分别为：10KV 35KV 10KV 35KV（2）变电所的负荷分别为：50MW 40MW 50MW 60MW （3）每个变电所的功率因数均为cos=0.8

3、5；（4）变电所2和变电所4分别配有两台容量为75MVA的变压器，短路损耗414KW，短路电压（%）=16.7；变电所1和变电所3分别配有两台容量为63MVA的变压器，短路损耗为245KW，短路电压（%）=10.5；4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为。1.3 课程设计基本内容1.对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。2.输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。3.跟随变电所负荷按一定比例发生变化，进行潮流计算分析。1)4

4、个变电所的负荷同时以2%的比例增大；2)4个变电所的负荷同时以2%的比例下降3)1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升；4.在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电压的调整。（变电所低压母线电压10KV要求调整范围在9.5-10.5之间；电压35KV要求调整范围在35-36之间）5.轮流断开环网一回线，分析潮流的分布。6.利用DDRTS软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比较。7.最终形成课程设计成品说明书。二、问题分析2.1 节点设置及分类 根据系统图可知此系统为两端供电网路，将母线1，2设

5、为节点1，10，将变电所1、2、3、4的高低压侧分别设为节点2、3、4、5、6、7、8、9。并且，将节点1设为平衡节点，将节点10设为PV节点，其余节点设为PQ节点。2.2 参数求取设定基准值，所以根据题目原始资料，计算发电厂、变压器及线路的参数。（1）运用下列公式计算变压器参数： （2）计算线路参数 （3）变电所负荷分别为： 变电所1 =50+j30.987 变电所2 =40+j27.79 变电所3 =50+j30.987 变电所4 =60+j37.18将参数整理，见下表：首端号末端号阻抗有名值阻抗标幺值电纳有名值电纳标幺值128.5+j20.10.0176+j0.0415j0.j0.269

6、11413.6+j32.160.0218+j0.0664j0.j0.10761613.6+j32.160.0218+j0.0664j0.j0.1076231.495+j40.3350.0031+j0.083300451.78+j53.8850.037+j0.1113004610.2+j24.120.0211+j0.0498j0.j0.0807671.495+j40.3350.0031+j0.083300686.8+j16.080.0140+j0.0332j0.j0.2153891.78+j53.8850.0037+j0.1113008108.5+j20.10.0176+j0.0415j0.j0

7、.2691（4）计算变压器分接头变比变压器有5个抽头，电压调节范围为2*2.5%， 对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比以备调压时选用对变电所低压母线为35KV时，非标准变比与10KV时相同。2.3 计算方法利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用PSAT进行潮流分析，将两者进行比较。2.4 牛顿拉夫逊法1、牛顿拉夫逊法概要首先对一般的牛顿拉夫逊法作一简单的说明。已知一个变量X函数为： 到此方程时，由适当的近似值出发，根据：反复进

8、行计算，当满足适当的收敛条件就是上面方程的根。这样的方法就是所谓的牛顿拉夫逊法。这一方法还可以做下面的解释，设第次迭代得到的解语真值之差，即的误差为时，则：把在附近对用泰勒级数展开上式省略去以后部分的误差可以近似由上式计算出来。比较两式，可以看出牛顿拉夫逊法的休整量和的误差的一次项相等。用同样的方法考虑，给出个变量的个方程：对其近似解得修正量可以通过解下边的方程来确定：式中等号右边的矩阵都是对于的值。这一矩阵称为雅可比（JACOBI）矩阵。按上述得到的修正向量后，得到如下关系这比更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反复计算满足为预先规定的小正数，是第n次迭代的近似值。2、牛

9、顿法的框图及求解过程（1）用牛顿法计算潮流时，有以下的步骤：给这各节点电压初始值；将以上电压初始值代入公式，求修正方程的常数项向量 ；将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。解修正方程式；修正各节点电压，；将，在带入方程式，求出；检验是否收敛，即如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回进行下次迭代计算，直到收敛为止。（2）程序框图如下：启动输入原始数据形成节点导纳矩阵分解各节点初始电压的实部和虚部迭代次数K=0求PQ节点的,，求PV节点的，置节点号i=0雅克比矩阵是否形成，in求得雅克比矩阵各元素增大节点号i=i+1把雅克

10、比矩阵单位化求解修正方程，得，求解最大修正量，是否收敛 回带各电压新值，K=K+1计算输出电压大小及相角，节电功率及支路损耗停止否是否是三、问题求解3.1 等值电路的计算电压是衡量电力系统电能质量的标准之一。电压过高或过低，都将对人身及其用电设备产生重大的影响。保证用户的电压接近额定值是电力系统调度的基本任务之一。当系统的电压偏离允许值时，电力系统必须应用电压调节技术调节系统电压的大小，使其维持在允许值范围内。本文经过手算形成了等值电路图，并编写好了程序得出节点电压标幺值，使其满足所要求的调整范围。我们首先对给定的程序输入部分作了简要的分析，程序开始需要我们确定输入节点数、支路数、平衡母线号、

11、支路参数矩阵、节点参数矩阵。（1）为了保证整个系统潮流计算的完整性，我们把凡具有母线及发电机处均选作节点，这样，可以确定有10个节点，节点号见等值电路图。（2）确定完节点及编号后，各条支路也相应确定了，网络中总计有13条支路，我们对各支路参数进行了计算。根据所给实际电路图和题中的已知条件，有以下公式计算各输电线路的阻抗和对地支路电容的标幺值和变压器的阻抗标幺值。选择电压基准值为和功率基准值所以3.2 依据题目要求及原始资料画出系统等值电路图：j 0.S=50+j30.987S=40+j27.79S=50+j30.987S=40+j24.8发电厂一发电厂二1:k11:k21:k31:k41102

12、48635798.5+j20.113.6+j32.161.495+j40.3351.78+j53.8851.495+j40.3351.78+j53.88510.2+j24.1213.6+j32.166.8+j16.088.5+j20.1j0.j0.j 0.4j0.j 0.=0U=231VP=200MWU=231V图1 系统等值电路利用psat软件可得计算所需系统图如图2所示。图23.3 潮流计算3.3.1 B1、B2矩阵的形成：根据所求参数，以及B1矩阵的含义，列写B1矩阵如下:B1=1 2 8.5+20.1i 0.i 1 0; 1 4 13.6+32.16i 0.i 1 0; 1 6 13.

13、6+32.16i 0.i 1 0; 2 3 1.495+40.335i 0 1 1; 4 5 1.78+53.885i 0 1 1; 4 6 10.2+24.12i 0.i 1 0; 6 7 1.495+40.335i 0 1 1; 6 8 6.8+16.08i 0.i 1 0; 8 9 1.78+53.885i 0 1 1; 8 10 8.5+20.1i 0.i 1 0;由各个变电所负荷功率可以计算出总功率为200MW，而发电厂一、二的总装机容量分别为400 MW和200 MW。令发电厂二的功率为200 MW，为了减小线路上的损耗，令发电机的电压为额定电压的1.05倍。根据前面叙述的节点分类

14、，形成B2矩阵如下：B2=0 0 231 231 0 1; 0 0 220 0 0 2; 0 50+30.987i 220 0 0 2; 0 0 220 0 0 2; 0 40+27.79i 220 0 0 2; 0 0 220 0 0 2; 0 50+30.987i 220 0 0 2; 0 0 220 0 0 2; 0 60+37.18i 220 0 0 2; 200 0 231 231 0 3;3.3.2 编写程序并运行从结果中可以观察到各个节点电压标幺值分别为:节点12345678910电压1.05001.03351.00551.03241.01901.02871.00061.0283

15、0.98171.0500折算到有名值以后，可以观察到节点3、5、7、9的电压都与题目给定的范围相比偏高，因此调节变压器分接头和发电厂电压，通过对系统进行多次调整，最终得到合理结果，调整结果如下: 调节方法：调节方法电厂1电压电厂2电压分接头1分接头2分接头3分接头4未调整1.051.051111调整后1.04091.04091.051.0251.0251.025 调节结果： 节点号12345678910未调整1.05001.03351.00551.03241.01901.02871.00061.02830.98171.0500调整后1.04091.02400.94541.01590.95621

16、.01610.96201.01700.94351.0409由上表观察到，进行最后一次调节时节点3，5，7，9的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,2)支路0.-28.0836i0.-27.5332i(1,4)支路0.-11.1554i0.-10.3807i(1,6)支路0.95644-9.3679i1.1139-8.7545i(2,3)支路0.10571+2.8521i0.11959+3.2266i(4,5)支路0.+1.8257i0.+2.8887i(4,6)支路0.51037-7.3677i0.45013-7

17、.2695i(6,7)支路0.10676+2.8803i0.11549+3.116i(6,8)支路2.74909-16.2733i2.75176-15.7397i(8,9)支路0.19012+5.7553i0.20583+6.2309i(8,10)支路6.44415-13.8237i6.52583-13.0626i总损耗11.8159-72.7582i12.2915-67.478i由此可以看出，随着变压器分接头的调高有功损耗逐渐增大，符合实际情况，而最终调节使得电压在规定范围内时，有功损耗增加了0.4756MW，相对与实际损耗12.29MW小很多，可以认为是合理的。具体的潮流分布见下表各条支路

18、的首端功率Si各条支路的末端功率SjS(1,2)=50.608+6.68043iS(2,1)=-50.1196-34.2136iS(1,4)=-2.65728+34.5734iS(4,1)=3.0825-44.9541iS(1,6)=-35.6592+49.1553iS(6,1)=36.7731-57.9099iS(2,3)=50.1196+34.2136iS(3,2)=-50-30.987iS(4,5)=40.0954+30.6787iS(5,4)=-40-27.79iS(4,6)=-43.1779+14.2754iS(6,4)=43.628-21.5448iS(6,7)=50.1155+3

19、4.103iS(7,6)=-50-30.987iS(6,8)=-130.5166+45.35175iS(8,6)=133.2683-61.09145iS(8,9)=60.2058+43.4109iS(9,8)=-60-37.18iS(8,10)=-193.4742+17.6806iS(10,8)=200-30.74322i各节点的功率S为(节点号从小到大排列)：0.1229 + 0.9041i 0.0000 - 0.0000i -0.5000 - 0.3099i 0.0000 - 0.0000i -0.4000 - 0.2779i -0.0000 - 0.0000i -0.5000 - 0.3

20、099i -0.0000 + 0.0000i -0.6000 - 0.3718i 2.0000 - 0.3074i得到的各节点的电压标幺值，电压角度，节点注入无功，支路首端无功如图3所示：图33.3.3 4个变电所负荷同时以2%的比例增大由于变电所负荷以2%的比例增大，所以各个变电所负荷变为：变电所1 =51+j31.6067变电所2 =40.8+j28.3458变电所3 =51+j31.6067变电所4 =61.2+j37.9236相应地， B2矩阵变为：B2=0 0 231 231 0 1; 0 0 220 0 0 2; 0 51+31.6067i 220 0 0 2; 0 0 220 0

21、 0 2; 0 40.8+28.3458i 220 0 0 2; 0 0 220 0 0 2; 0 51+31.6067i 220 0 0 2; 0 0 220 0 0 2; 0 61.2+37.9236i 220 0 0 2; 200 0 231 231 0 3;将程序中的B2矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表：调节方法：调节方法电厂1电压电厂2电压分接头1分接头2分接头3分接头4调整前1.051.051111调整后1.04091.04091.051.0251.0251.025调节结果： 节点电压12345678910调整前1.05001.033

22、01.00441.02500.99071.02520.99641.02610.97841.0500调整后1.04091.02350.94421.01520.95481.01550.96071.01630.94181.0409由上表观察到，进行最后一次调节时节点3，5，7，9的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,2)支路0.-28.0178i0.-27.4654i(1,4)支路0.-10.6092i0.-10.3521i(1,6)支路1.057-9.0914i1.0928-8.7972i(2,3)支路0.11022

23、+2.9737i0.12472+3.365i(4,5)支路0.+2.793i0.+3.0144i(4,6)支路0.43301-7.4596i0.44047-7.2813i(6,7)支路0.11201+3.022i0.12048+3.2506i(6,8)支路2.6589-16.3602i2.70477-15.8222i(8,9)支路0.19916+6.029i0.21494+6.5066i(8,10)支路6.41198-13.8387i6.51806-13.0635i总损耗11.9856-70.5592i12.2612-66.6451i由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题

24、目要求，可认为是合理的。具体的潮流分布如下：各条支路的首端功率Si各条支路的末端功率SjS(1,2)=51.622+6.89876iS(2,1)=-51.1223-34.9067iS(1,4)=-1.44053+34.6728iS(4,1)=1.8613-45.2634iS(1,6)=-34.1421+48.9161iS(6,1)=35.206-57.9894iS(2,3)=51.1223+34.9067iS(3,2)=-51-31.6067iS(4,5)=40.9028+31.458iS(5,4)=-40.8-28.3458iS(4,6)=-42.7641+13.8053iS(6,4)=43

25、.1965-21.2633iS(6,7)=51.1181+34.7927iS(7,6)=-51-31.6067iS(6,8)=-129.5206+44.46004iS(8,6)=132.1795-60.8131iS(8,9)=61.4105+44.2948iS(9,8)=-61.2-37.9236iS(8,10)=-193.5899+16.51826iS(10,8)=200-30.35715i各节点的功率S为(节点号从小到大排列)： 0.1626+0.9202i -0.0000+0.0000i -0.5100-0.3161i 0.0000+0.0000i -0.4080-0.2835i -0.

26、0000+0.0000i -0.5100-0.3161i -0.0000+0.0000i -0.6120-0.3792i 2.0000-0.2919i得到的图像如下：图43.3.4 4个变电所负荷同时以2%的比例下降由于变电所负荷以2%的比例增大，所以各个变电所负荷变为：变电所1 =49+j30.3673变电所2 =39.2+j27.2342变电所3 =49+j30.3673变电所4 =58.8+j36.4364相应地， B2矩阵变为：B2=0 0 231 231 0 1; 0 0 220 0 0 2; 0 49+30.3673i 220 0 0 2; 0 0 220 0 0 2; 0 39.

27、2+27.2342i 220 0 0 2; 0 0 220 0 0 2; 0 49+30.3673i 220 0 0 2; 0 0 220 0 0 2; 0 58.8+36.4364i 220 0 0 2; 200 0 231 231 0 3;将程序中的B2矩阵替换后，运行并进行调节，对应的调节方法和结果见下表：调节方法：调节方法电厂1电压电厂2电压分接头1分接头2分接头3分接头4调整前1.051.051111调整后1.04091.04091.051.0251.0251.025调节结果见下表：节点号12345678910调整前1.05001.03401.00661.02620.99341.02

28、650.99891.02730.98171.0500调整后1.04091.02450.94651.01660.95761.01680.96331.01760.94531.0409由上表观察到，进行最后一次调节时节点3，5，7，9的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,2)支路0.-28.148i0.-27.5994i(1,4)支路0.-10.6572i0.-10.4065i(1,6)支路1.1018-8.9996i1.1366-8.7081i(2,3)支路0.1013+2.7332i0.11458+3.0914i(

29、4,5)支路0.+2.5705i0.+2.766i(4,6)支路0.45176-7.4365i0.45987-7.2573i(6,7)支路0.10287+2.7756i0.11062+2.9845i(6,8)支路2.75146-16.1972i2.79906-15.6563i(8,9)支路0.18259+5.5273i0.19695+5.9622i(8,10)支路6.42888-13.833i6.53433-13.06i总损耗12.0587-71.6649i12.3269+67.8835i具体的潮流分布如下表各条支路的首端功率Si各条支路的末端功率SjS(1,2)=49.5807+5.8593

30、4iS(2,1)=-49.1146-33.4587iS(1,4)=-3.95599+34.0901iS(4,1)=4.37337-44.4966iS(1,6)=-37.2978+48.8636iS(6,1)=38.4344-57.5717iS(2,3)=49.1146+33.4587iS(3,2)=-49-30.3673iS(4,5)=39.2914+30.0002iS(5,4)=-39.2-27.2342iS(4,6)=-43.6647+14.4964iS(6,4)=44.1246-21.7537iS(6,7)=49.1106+33.3518iS(7,6)=-49-30.3673iS(6,8

31、)=-131.6697+45.97358iS(8,6)=134.4687-61.62989iS(8,9)=58.997+42.3986iS(9,8)=-58.8-36.4364iS(8,10)=-193.4657+19.23132iS(10,8)=200-32.29128i各节点的功率S为(节点号从小到大排列)： 0.0833+0.8881i 0.0000+0.0000i -0.4900-0.3037i -0.0000-0.0000i -0.3920-0.2723i 0.0000-0.0000i -0.4900-0.3037i -0.0000-0.0000i -0.5880-0.3644i 2

32、.0000-0.3229i得到的各节点的电压标幺值，电压角度，节点注入无功，支路首端无功情况见图7。图53.3.5 1、4号变电所负荷同时以2%的比例下降，2、3号变电所负荷同时以2%的比例增加当1、4号变电所负荷同时以2%的比例下降，2、3号变电所负荷同时以2%的比例增加，各个变电所负荷变为：变电所1 =49+j30.3673变电所2 =40.8+j28.3458变电所3 =51+j31.6067变电所4 =58.8+j36.4364相应地， B2矩阵变为：B2=0 0 231 231 0 1; 0 0 220 0 0 2; 0 49+30.3673i 220 0 0 2; 0 0 220

33、0 0 2; 0 40.8+28.3458i 220 0 0 2; 0 0 220 0 0 2; 0 51+31.6067i 220 0 0 2; 0 0 220 0 0 2; 0 58.8+36.4364i 220 0 0 2; 200 0 231 231 0 3;将程序中的B2矩阵替换后，运行并进行调节，对应的调节方法和结果见下表：调节方法：调节方法电厂1电压电厂2电压分接头1分接头2分接头3分接头4调整前1.051.051111调整后1.04091.04091.051.0251.0251.025调节结果见下表：节点号12345678910调整前1.05001.03401.00661.02

34、510.99081.02560.99671.02680.98121.0500调整后1.04091.02450.94651.01550.95501.01590.96111.01710.94471.0409由上表观察到，进行最后一次调节后节点3，5，7，9的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,2)支路0.-28.148i0.-27.5994i(1,4)支路0.-10.6052i0.-10.352i(1,6)支路1.0853-9.0286i1.1204-8.736i(2,3)支路0.1013+2.7332i0.1145

35、8+3.0914i(4,5)支路0.+2.7991i0.+3.013i(4,6)支路0.44707-7.4308i0.45507-7.2518i(6,7)支路0.11192+3.0197i0.12038+3.2479i(6,8)支路2.73794-16.1975i2.78523-15.6571i(8,9)支路0.18279+5.5334i0.19717+5.969i(8,10)支路6.42155-13.8362i6.52731-13.0623i总损耗12.0503+71.1609i12.3211+67.3373i 由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理

36、的。具体的潮流分布如下表各条支路的首端功率Si各条支路的末端功率SjS(1,2)=49.5807+5.85934iS(2,1)=-49.1146-33.4587iS(1,4)=-2.18717+35.0781iS(4,1)=2.62244-45.4301iS(1,6)=-35.4724+49.5049iS(6,1)=36.5929-58.2409iS(2,3)=49.1146+33.4587iS(3,2)=-49-30.3673iS(4,5)=40.8995+31.3588iS(5,4)=-40.8-28.3458iS(4,6)=-43.522+14.0713iS(6,4)=43.977-21

37、.3231iS(6,7)=51.1204+34.8546iS(7,6)=-51-31.6067iS(6,8)=-131.6903+44.70935iS(8,6)=134.4755-60.36641iS(8,9)=58.9972+42.4054iS(9,8)=-58.8-36.4364iS(8,10)=-193.4727+17.96105iS(10,8)=200-31.02333i各节点的功率S为(节点号从小到大排列)： 0.1192+0.9044i 0.0000+0.0000i -0.4900-0.3037i -0.0000-0.0000i -0.4080-0.2835i -0.0000-0.

38、0000i -0.5100-0.3161i -0.0000-0.0000i -0.5880-0.3644i 2.0000-0.3102i得到的各节点的电压标幺值，电压角度，节点注入无功，支路首端无功情况见图6。图63.3.6 断开1、4支路时的潮流分析当断开1、4支路时，导致支路数变为9和B1矩阵改变，其他参数不变。用MATLAB软件进行潮流计算进行调整。B1矩阵变为B1=1 2 8.5+20.1i 0.i 1 0; 1 6 13.6+32.16i 0.i 1 0; 2 3 1.495+40.335i 0 1.05 1; 4 5 1.78+53.885i 0 1.05 1; 4 6 10.2+

39、24.12i 0.i 1 0; 6 7 1.495+40.335i 0 1.025 1; 6 8 6.8+16.08i 0.i 1 0; 8 9 1.78+53.885i 0 1.025 1; 8 10 8.5+20.1i 0.i 1 0;节点号12345678910调整前1.05001.03390.95570.99460.92031.01270.95791.01960.94731.0500调整后1.03641.01980.94190.98220.93201.00050.94551.00850.93591.0409由上表观察到，进行最后一次调节后节点3，5，7，9的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗