课程设计电力系统潮流计算.doc

电力系统潮流计算课程设计任务书一 .题目原始资料1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。变电所1变电所2变电所3变电所435kV母线10kV母线35kV母线10kV母线一次侧电压220kV一次侧电压220kV母线1母线3母线2线路长为70km线路长为90km线路长为90km线路长为100km线路长为80km2*QFQ-50-22*QFS-50-2TQN-100-22*TQN-100-2。线路长为70km。电厂一电厂二线路长为80km2、发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为40MW和20MW；发电厂二总装机容量为（ 200MW ）。3、变电所资料：(一) 变电所1、2、3、4低压母线的电压等级分别为：10kV 10kV 35kV 35kV(二) 变电所的负荷分别为： （4）50MW 50MW 60MW 70MW （三）每个变电所的功率因数均为cos=0.85；（四）变电所3和变电所4分别配有两台容量为75MVA的变压器，短路损耗414kW，短路电压（%）=16.7；变电所1和变电所2分别配有两台容量为63MVA的变压器，短路损耗为245kW，短路电压（%）=10.5；4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为。二、 课程设计基本内容：1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。1) 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. 最终形成课程设计成品说明书。三、课程设计成品基本要求：1. 在读懂程序的基础上画出潮流计算基本流程图2. 通过输入数据，进行潮流计算输出结果3. 对不同的负荷变化，分析潮流分布，写出分析说明。4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控制的参数变化。5. 打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。摘要本文运用MATLAB软件进行潮流计算，对给定题目进行分析计算，在应用DDRTS软件，构建系统图进行仿真，最终进行合理的系统潮流。潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种，当系统相对简单，节点数相对较少时，我们可以通过手算获得各参数，但当系统网络较复杂时，想要通过手算各参数显然是不可能的，这就为计算机算法的出现提出了要求。牛顿拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，利用MATLAB编制计算程序。 最后，应用DDRTS软件，构建系统图，对给定负荷重新进行分析，潮流计算后的结果也能满足相应的参数要求。关键词：牛顿-拉夫逊法 潮流计算 MATLAB DDRTS一、 牛顿拉夫逊潮流计算法的步骤 用牛顿法计算潮流时，有以下的步骤：（1） 给这各节点电压初始值；（2） 将以上电压初始值代入公式，求修正方程的常数项向量；（3） 将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。（4） 解修正方程式；（5） 修正各节点电压，；（6） 将，在带入方程式，求出；（7） 检验是否收敛，即如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功二、 求解方法利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分布，根据题目的不同要求对参数进行调整，通过调节变压器变比和无功补偿功率求解出合理的潮流分布，最后利用DDRTS进行潮流分析，将两者进行比较。 题目求解一、 程序对应的流程图 二、依据题目要求及原始资料画出等值电路图三、变电所负荷为题目所给数据进行求解 1、线路参数：支路首端号支路末端号线路阻抗 对地电纳920.28+0.66i0.1076i950.003+0.083i0910.025+0.058i0.0942260.003+0.083i0210.032+0.057i0.121i340.025+0.058i0.0942i370.0037+0.111i0310.018+0.042i0.269i3100.028+0.066i0.1076i480.0037+0.111i04100.032+0.075i0.i121 2、形成B1、B2矩阵 根据所求参数以及B1矩阵的含义，列写B1矩阵如下：B1=9 2 0.28+0.066i 0.1076i 1 0； 9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.018+0.042i 0.269i 1 0; 3 10 0.28+0.066i 0.1076i 1 0; 4 8 0.0037+0.111i 0 1.05 0; 4 10 0.032+0.075i 0.121i 1 0; 根据发电厂、变电所的容量及B2矩阵的含义，列写B2矩阵：B2=3 0 1.05 1.05 0 1; 0 0 1 0 0 2; 0 0 1 0 0 2; 0 0 1 0 0 2; 0 0.50+0.31i 1 0 0 2; 0 0.50+0.31i 1 0 0 2; 0 0.60+0.37i 1 0 0 2; 0 0.70+0.434i 1 0 0 2; 0 0 1 0 0 2; 2 0 1.05 1.05 0 3; 3、变压器低压侧电压标幺值的变化范围 四个变电所对应的节点号依次是5、6、7、8 .变电所号1234节点号5678电压下限0.950.9511电压上限1.051.051.02851.0285 4、利用MATLAB程序运行的结果未调整时各节点的标幺值如下节点12345678910电压1.05001.01781.00340.99651.04571.03891.00360.98611.02411.0500由调压范围可知变电所4不满足情况。 调节方法：通过改变变电所变压器变比或无功功率补偿来调整电压。调节方法变电所1无功补偿变电所2无功补偿变电所3无功补偿变电所4无功补偿变电所1变比变电所2变比变电所3变比变电所4变比未调整00001.051.051.051.05调整后00001.051.051.051.077 调整后的B1矩阵、数据及图：B1=9 2 0.28+0.066i 0.1076i 1 0； 9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.018+0.042i 0.269i 1 0; 3 10 0.28+0.066i 0.1076i 1 0; 4 8 0.0037+0.111i 0 1.077 0; 4 10 0.032+0.075i 0.121i 1 0;B2矩阵保持不变。调整后各节点的标幺值如下节点12345678910电压1.05001.01781.00340.99651.04571.03891.00361.01151.02411.0500图像如下： 各支路的功率损耗如下表：功率损耗未调整调整后（9,2）支路0.0176302-11.2117i0.0176302-11.2117i（9,5）支路0.10818+2.8964i0.10818+2.8964i（9,1）支路0.78049-8.3217i0.78049-8.3217i(2,6) 支路0.1096+2.9343i0.1096+2.9343i(2,1) 支路0.86259-10.916i0.86259-10.916i(3,4) 支路2.3582-3.9485i2.3582-3.9485i(3,7) 支路0.20124+6.0373i0.20124+6.0373i(3,1) 支路3.41261-20.4083i3.41261-20.4083i(3,10) 支路11.5125-8.63479i11.5125-8.63479i(4,8) 支路0.28457+8.5371i0.28457+8.5371i(4,10) 支路7.8174+5.6441i7.8174+5.6441i 各节点的功率S为（节点号从小到大排列） 0.5746 + 1.6014i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.5000 - 0.3100i -0.5000 - 0.3100i -0.6000 - 0.3700i -0.7000 - 0.4340i 0.0000 - 0.0000i 2.0000 - 0.5513i调整后的具体潮流分布各条支路的首端功率Si各条支路的末端功率SjS(9,2)=2.4652-6.3679iS(2,9)=-2.4476-4.8438iS(9,5)=50.1082+33.8964iS(5,9)=-50-31iS(9,1)=-52.5734-27.5285iS(1,9)=53.3539+19.2068iS(2,6)=50.1096+33.9343iS(6,2)=-50-31iS(2,1)=-47.662-29.0905iS(1,2)=48.5246+18.1745iS(3,4)=-83.0508+46.2485iS(4,3)=85.409-50.1971iS(3,7)=60.2012+43.0373iS(7,3)=-60-37iS(3,1)=47.82618-143.1649iS(1,3)=-44.41357+122.7566iS(3,10)=-24.9766+53.8791iS(10,3)=36.4891-62.5139iS(4,8)=70.2846+51.9371iS(8,4)=-70-43.4iS(4,10)=-155.6935-1.740071iS(10,4)=163.5109+7.384212i四、4个变电所的负荷同时以2%的比例增大；各变电所变化后功率的标幺值：变电所1: 0.51+0.3162i变电所2: 0.51+0.3162i变电所3: 0.612+0.3774i变电所4: 0.714+0.4427i B1矩阵不变，B2矩阵变为：B2=3 0 1.05 1.05 0 1;0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 0.51+0.3162i 1 0 0 2;0 0.51+0.3162i 1 0 0 2;0 0.612+0.3774i 1 0 0 2;0 0.714+0.4427i 1 0 0 2;0 0 1 0 0 2; 2 0 1.05 1.05 0 3;调节方法：调节方法变电所1无功补偿变电所2无功补偿变电所3无功补偿变电所4无功补偿变电所1变比变电所2变比变电所3变比变电所4变比未调整00001.051.051.051.05调整后00001.051.051.051.077调整后的B1矩阵、数据及图：B1=9 2 0.28+0.066i 0.1076i 1 0； 9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.018+0.042i 0.269i 1 0; 3 10 0.28+0.066i 0.1076i 1 0; 4 8 0.0037+0.111i 0 1.077 0;4 10 0.032+0.075i 0.121i 1 0;B2矩阵不变调整前后各节点的标幺值如下：节点电压12345678910未调整1.05001.01691.00280.99561.04431.03731.00180.98371.02341.0500调整后1.05001.01691.00280.99561.04431.03731.00181.00901.02341.0500由调压范围可知变电所4不满足情况。各支路的功率损耗如下表：功率损耗未调整调整后（9,2）支路0.0185012-11.1937i0.0185012-11.1937i（9,5）支路0.11285+3.0215i0.11285+3.0215i（9,1）支路0.81706-8.2299i0.81706-8.2299i(2,6) 支路0.11438+3.0624i0.11438+3.0624i(2,1) 支路0.903181-10.8098i0.903181-10.8098i(3,4) 支路2.3002-4.0685i2.3002-4.0685i(3,7) 支路0.21012+6.3036i0.21012+6.3036i(3,1) 支路3.39598-20.4307i3.39598-20.4307i(3,10) 支路11.3924-8.6565i11.3924-8.6565i(4,8) 支路0.29754+8.9263i0.29754+8.9263i(4,10) 支路7.8237+5.6702i7.8237+5.6702i各节点的功率S为(节点号从小到大排列)： 0.6199 + 1.6224i -0.0000 + 0.0000i 0.0000 - 0.0000i -0.0000 - 0.0000i -0.5100 - 0.3162i -0.5100 - 0.3162i -0.6120 - 0.3774i -0.7140 - 0.4427i 0.0000 - 0.0000i 2.0000 - 0.5340i调整后的具体潮流分布各条支路的首端功率Si各条支路的末端功率SjS(9,2)=2.5293-6.3574iS(2,9)=-2.5108-4.8363iS(9,5)=51.1129+34.6415iS(5,9)=-51-31.62iS(9,1)=-53.6421-28.2841iS(1,9)=54.4592+20.0542iS(2,6)=51.1144+34.6824iS(6,2)=-51-31.62iS(2,1)=-48.6036-29.8461iS(1,2)=49.5068+19.0363iS(3,4)=-81.6404+46.1307iS(4,3)=83.9406-50.1992iS(3,7)=61.4101+44.0436iS(7,3)=-61.2-37.74iS(3,1)=45.37611-143.5805iS(1,3)=-41.98012+123.1499iS(3,10)=-25.1458+53.4062iS(10,3)=36.5382-62.0627iS(4,8)=71.6975+53.1963iS(8,4)=-71.4-44.27iS(4,10)=-155.6381-2.997048iS(10,4)=163.4618+8.667289i图像如下：五、4个变电所的负荷同时以2%的比例减小；各变电所变化后功率的标幺值：变电所1: 0.49+0.3038i变电所2: 0.49+0.3038i变电所3: 0.588+0.3626i变电所4: 0.686+0.4253i B1矩阵不变，B2矩阵变为：B2=3 0 1.05 1.05 0 1;0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 0.49+0.3038i 1 0 0 2;0 0.49+0.3038i 1 0 0 2;0 0.588+0.3626i 1 0 0 2;0 0.686+0.4253i 1 0 0 2;0 0 1 0 0 2; 2 0 1.05 1.05 0 3;未调整时各节点的标幺值如下节点12345678910电压1.05001.01871.00400.99741.04711.04051.00540.98851.02481.0500由调压范围可知变电所4不满足情况。调节方法：调节方法变电所1无功补偿变电所2无功补偿变电所3无功补偿变电所4无功补偿变电所1变比变电所2变比变电所3变比变电所4变比未调整00001.051.051.051.05调整后00001.051.051.051.077调整后的B1矩阵、数据及图：B1=9 2 0.28+0.066i 0.1076i 1 0； 9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.018+0.042i 0.269i 1 0; 3 10 0.28+0.066i 0.1076i 1 0; 4 8 0.0037+0.111i 0 1.077 0; 4 10 0.032+0.075i 0.121i 1 0;B2矩阵不变调整后各节点的标幺值如下节点12345678910电压1.05001.01871.00400.99741.04711.04051.00541.01391.02481.0500图像如下：各支路的功率损耗如下表：功率损耗未调整调整后（9,2）支路0.0167858-11.2297i0.0167858-11.2297i（9,5）支路0.10362+2.7742i0.10362+2.7742i（9,1）支路0.7449-8.4112i0.7449-8.4112i(2,6) 支路0.10493+2.8095i0.10493+2.8095i(2,1) 支路0.823097-11.0197i0.823097-11.0197i(3,4) 支路2.417-3.8266i2.417-3.8266i(3,7) 支路0.19259+5.7778i0.19259+5.7778i(3,1) 支路3.43175-20.38i3.43175-20.38i(3,10) 支路11.633-8.61292i11.633-8.61292i(4,8) 支路0.27196+8.1589i0.27196+8.1589i(4,10) 支路7.812+5.6203i7.812+5.6203i调整后的具体潮流分布各条支路的首端功率Si各条支路的末端功率SjS(9,2)=2.4013-6.3781iS(2,9)=-2.3845-4.8515iS(9,5)=49.1036+33.1542iS(5,9)=-49-30.38iS(9,1)=-51.5049-26.7761iS(1,9)=52.2498+18.3649iS(2,6)=49.1049+33.1895iS(6,2)=-49-30.38iS(2,1)=-46.7204-28.3379iS(1,2)=47.5435+17.3182iS(3,4)=-84.4577+46.3697iS(4,3)=86.87467-50.19633iS(3,7)=58.9926+42.0378iS(7,3)=-58.8-36.26iS(3,1)=50.27346-142.7569iS(1,3)=-46.84171+122.3768iS(3,10)=-24.8084+54.3494iS(10,3)=36.4414-62.9623iS(4,8)=68.872+50.6889iS(8,4)=-68.6-42.53iS(4,10)=-155.7466-0.492552iS(10,4)=163.5586+6.112814i各节点的功率S为(节点号从小到大排列)： 0.5295 + 1.5806i -0.0000 + 0.0000i 0.0000 - 0.0000i -0.0000 - 0.0000i -0.4900 - 0.3038i -0.4900 - 0.3038i -0.5880 - 0.3626i -0.6860 - 0.4253i -0.0000 - 0.0000i 2.0000 - 0.5685i六、1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升；各变电所变化后功率的标幺值：变电所1: 0.49+0.3038i变电所2: 0.51+0.3162i变电所3: 0.612+0.3774i变电所4: 0.686+0.4253i B1矩阵不变，B2矩阵变为：B2=3 0 1.05 1.05 0 1;0 0 1 0 0 2;0 0 1 0 0 2;0 0 1 0 0 2;0 0.49+0.3038i 1 0 0 2;0 0.51+0.31662i 1 0 0 2;0 0.612+0.3774i 1 0 0 2;0 0.686+0.4253i 1 0 0 2;0 0 1 0 0 2; 2 0 1.05 1.05 0 3;未调整时各节点的标幺值如下节点12345678910电压1.05001.01711.00320.99701.04691.03751.00220.98801.02471.0500由调压范围可知变电所4不满足情况。调节方法：调节方法变电所1无功补偿变电所2无功补偿变电所3无功补偿变电所4无功补偿变电所1变比变电所2变比变电所3变比变电所4变比未调整00001.051.051.051.05调整后00001.051.051.051.077调整后的B1矩阵、数据及图：B1=9 2 0.28+0.066i 0.1076i 1 0； 9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.018+0.042i 0.269i 1 0; 3 10 0.28+0.066i 0.1076i 1 0; 4 8 0.0037+0.111i 0 1.077 0; 4 10 0.032+0.075i 0.121i 1 0;B2矩阵不变调整后各节点的标幺值如下节点12345678910电压1.05001.01711.00320.99701.04691.03751.00221.01341.02471.0500图像如下：各支路的功率损耗如下表：功率损耗未调整调整后（9,2）支路0.0245492-11.2081i0.0245492-11.2081i（9,5）支路0.10365+2.7751i0.10365+2.7751i（9,1）支路0.75875-8.3776i0.75875-8.3776i(2,6) 支路0.11434+3.0614i0.11434+3.0614i(2,1) 支路0.890437-10.8415i0.890437-10.8415i(3,4) 支路2.3949-3.8655i2.3949-3.8655i(3,7) 支路0.20994+6.2983i0.20994+6.2983i(3,1) 支路3.4421-20.3333i3.4421-20.3333i(3,10) 支路11.5865-8.61485i11.5865-8.61485i(4,8) 支路0.27225+8.1675i0.27225+8.1675i(4,10) 支路7.7977+5.5923i7.7977+5.5923i调整后的具体潮流分布各条支路的首端功率Si各条支路的末端功率SjS(9,2)=2.948-6.3661iS(2,9)=-2.9234-4.8419iS(9,5)=49.1036+33.1551iS(5,9)=-49-30.38iS(9,1)=-52.0516-26.789iS(1,9)=52.8103+18.4113iS(2,6)=51.1143+34.6814iS(6,2)=-51-31.62iS(2,1)=-48.1909-29.8395iS(1,2)=49.0814+18.998iS(3,4)=-84.2877+45.6256iS(4,3)=86.6826-49.4911iS(3,7)=61.4099+44.0383iS(7,3)=-61.2-37.74iS(3,1)=47.93872-143.718iS(1,3)=-44.49662+123.3847iS(3,10)=-25.061+54.0541iS(10,3)=36.6474-62.6689iS(4,8)=68.8723+50.6975iS(8,4)=-68.6-42.53iS(4,10)=-155.5549-1.206402iS(10,4)=163.3526+6.798688i各节点的功率S为(节点号从小到大排列)： 0.5740 + 1.6079i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.0000 - 0.0000i -0.4900 - 0.3038i -0.5100 - 0.3162i -0.6120 - 0.3774i -0.6860 - 0.4253i 0.0000 - 0.0000i 2.0000 - 0.5587i七、断开9、2节点之间一回线时的潮流计算B1=9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.018+0.042i 0.269i 1 0; 3 10 0.28+0.066i 0.1076i 1 0; 4 8 0.0037+0.111i 0 1.05 0;4 10 0.032+0.075i 0.121i 1 0;B1矩阵的支路数变为10，B2矩阵不变。B2=3 0 1.05 1.05 0 1; 0 0 1 0 0 2; 0 0 1 0 0 2; 0 0 1 0 0 2; 0 0.50+0.31i 1 0 0 2; 0 0.50+0.31i 1 0 0 2; 0 0.60+0.37i 1 0 0 2; 0 0.70+0.434i 1 0 0 2; 0 0 1 0 0 2; 2 0 1.05 1.05 0 3;未调整时各节点的标幺值如下节点12345678910电压1.05001.01331.00340.99651.04241.03391.00360.98611.02101.0500由调压范围可知变电所4不满足情况。调节方法：调节方法变电所1无功补偿变电所2无功补偿变电所3无功补偿变电所4无功补偿变电所1变比变电所2变比变电所3变比变电所4变比未调整00001.051.051.051.05调整后00001.051.051.051.077调整后的B1矩阵、数据及图：B1=9 5 0.0031+0.083i 0 1.05 0; 9 1 0.025+0.058i 0.0942i 1 0; 2 6 0.0031+0.083i 0 1.05 0; 2 1 0.032+0.075i 0.121i 1 0; 3 4 0.025+0.058i 0.0942i 1 0; 3 7 0.0037+0.111i 0 1.05 0; 3 1 0.01