东北电力大学电力系统潮流计算课程设计报告书

目录一、设计任务 (1)1.1 课程设计要求 (1)1.2 课程设计题目 (1)1.3 课程设计基本容 (2)二、问题分析 (3)2.1 节点设置及分类 (3)2.2 参数求取 (3)2.3 计算方法 (4)三、问题求解 (7)3.1 等值电路的计算 (7)3.2画出系统等值电路图： (7)3.3 潮流计算 (8)四、误差分析 (29)五、心得体会及总结 (38)附录: (39)参考文献 (39)程序 (39)电力系统潮流计算课程设计一、设计任务1.1 课程设计要求1、在读懂程序的基础上画出潮流计算基本流程图2、通过输入数据，进行潮流计算输出结果3、对不同的负荷变化，分析潮流分布，写出分析说明。

4、对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控制的参数变化。

5、打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。

1.2 课程设计题目系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。

变电所1 变电所2母线电厂一电厂二发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和20MW ；发电厂二总装机容量为（ 200MW ）。

变电所资料：(一)变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV (二)变电所的负荷分别为：60MW 40MW 70MW 50MW (三)每个变电所的功率因数均为cos φ=0.85；(四)变电所1和变电所3分别配有两台容量为75MVA 的变压器，短路损耗414KW ，短路电压（%）=16.7；变电所2和变电所4分别配有两台容量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5；输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。

1.3 课程设计基本容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. 最终形成课程设计成品说明书。

二、问题分析2.1 节点设置及分类2.2 参数求取设定基准值MVA S B 100=，KV U B 220=，所以4842==BB B U S Z Ω根据题目原始资料，计算发电厂、变压器及线路的参数。

（1）运用下列公式计算变压器参数：S UP RN N kT221000=S UUXNNKT1002%=（2）计算线路参数L jx r jX R Z )(+=+= Y = jB / 2 （3）变电所负荷分别为：变电所1 60+j37.184 MVA 变电所2 40+j24.790 MVA 变电所3 70+j43.382 MVA 变电所4 50+j30.987 MVA （4）计算变压器分接头变比变压器有5个抽头，电压调节围为U N ±2*2.5%， U N 对应分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比k *以备调压选用955.02201110%)51(2201*=⨯⨯+⨯=k932.02201110%)5.21(2202*=⨯⨯+⨯=k909.022011102203*=⨯⨯=k886.02201110%)5.21(2204*=⨯⨯-⨯=k864.02201110%)51(2205*=⨯⨯-⨯=k对变电所低压母线为35KV 时，非标准变比与10KV 时相同。

2.3 计算方法利用牛顿拉夫逊法进行求解，用MATLAB 软件编程，可以求解系统潮流分 布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用PSAT 进行潮流分析，将两者进行比较。

1、牛顿—拉夫逊法概要首先对一般的牛顿拉夫逊法作简单说明。

已知一个变量X 函数为：0)(=X f 到此方程时，由适当的近似值)0(X出发，根据：,......)2,1()()()()()()1(='-=+n X f X f XXn n n n 反复进行计算，当)(n X 满足适当的收敛条件就是上面方程的根。

这样的方法就是所谓的牛顿—拉夫逊法。

这一方法还可以做下面的解释，设第n 次迭代得到的解语真值之差，即)(n X 的误差为ε时，则：0)()(=+εn X f把)()(ε+n Xf 在)(n X附近对ε用泰勒级数展开0......)(!2)()()()(2)()()(=+''+'+=+n n n n X f Xf Xf Xf εεε上式省略去2ε以后部分 0)()()()(≈'+n n X f X f ε )(n X的误差可以近似由上式计算出来。

)()()()(n n X f X f '-≈ε比较两式，可以看出牛顿—拉夫逊法的休整量和)(n X 的误差的一次项相等。

用同样的方法考虑，给出n 个变量的n 个方程：⎪⎪⎩⎪⎪⎨⎧===0),,,(0),,,(0),,,(21212211n n n n X X X f X X X f X X X f ΛΛΛΛΛΛK Λ 对其近似解1X '得修正量1X '∆可以通过解下边的方程来确定： ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂-=⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡'''''''''n n n n n nn n n n n X X X x f x f xf x f x f x f x f x f x f X X X f X X X f X X X f M ΛM KΛM M ΛΛ2121222121211121212211),,,(),,,(),,,( 式中等号右边的矩阵nnx f ∂∂都是对于n X X X ''',,,21Λ的值,称为雅可比（JACOBI ）矩阵。

按上述得到的修正向量n X X X '∆'∆'∆,,,21Λ后，得到关系式:n n n X X X ∆+'='' 这比n X X X ''',,,21Λ更接近真实值。

这一步在收敛到希望的值以前重复进行，一般要反复计算满足{}ε<---++++++1112121111,,,max n nn nn n n n X X X X X X Λε为预先规定的小正数，1+n n X 是第n 次迭代n X 的近似值。

2、牛顿法的框图及求解过程（1）用牛顿法计算潮流时，有以下的步骤：①给这各节点电压初始值)0()0(,fe ；②将以上电压初始值代入公式，求修正方程的常数项向量)0(2)0()0()(,,V Q P ∆∆∆；③将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。

④解修正方程式)0()0(,f e ∆∆；⑤修正各节点电压)0()0()1(e e e ∆+=，)0()0()1(f f f ∆+=； ⑥将)1(e ，)1(f 在带入方程式，求出)1(2)1()1()(,,V Q P ∆∆∆；⑦检验是否收敛，即{}ε<∆∆)()(,max k ik iQ P如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。

如果不收敛，转回②进行下次迭代计算，直到收敛为止。

（2三、问题求解3.1 等值电路的计算电压是衡量电力系统电能质量的标准之一。

电压过高或过低，都将对人身及其用电设备产生重大的影响。

保证用户的电压接近额定值是电力系统调度的基本任务之一。

当系统的电压偏离允许值时，电力系统必须应用电压调节技术调节系统电压的大小，使其维持在允许值围。

本文经过手算形成了等值电路图，并编写好了程序得出节点电压标幺值，使其满足所要求的调整围。

我们首先对给定的程序输入部分作了简要的分析，程序开始需要我们确定输入节点数、支路数、平衡母线号、支路参数矩阵、节点参数矩阵。

（1）为了保证整个系统潮流计算的完整性，我们把凡具有母线及发电机处均选作节点，这样，可以确定有10个节点，节点号见等值电路图。

（2）确定完节点及编号后，各条支路也相应确定了，我们对网络中各支路参数进行了计算。

根据所给实际电路图和题中的已知条件，有以下公式计算各输电线路的阻抗和对地支路电容的标幺值和变压器的阻抗标幺值。

选择电压基准值为220B U KV =和功率基准值100B S MVA = 所以484BB BU Z S ==Ω 3.2画出系统等值电路图：( 注1:图中数据均为有名值,且已归算到220kV 侧 )( 注2:节点1处发电厂机压母线存在20-50MW 负荷 )3.3 潮流计算3.3.1 B1、B2矩阵的形成：B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳5、支路的变比；6、支路首端处于K侧为1，1侧为0B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量6、节点分类标号根据所求参数，以及B1矩阵的含义，列写B1矩阵如下:B1=[1 2 8.5+20.1i 0.000556i 1 0;1 4 4.25+10.05i 0.000278i 1 0;1 6 7.65+18.09i 0.0005004i 1 0;2 3 1.78+53.885i 0 0.909 1;4 5 1.495+40.335i 0 0.909 1;6 7 1.78+53.885i 0 0.909 1;6 8 6.8+16.008i 0.0004448i 1 0;8 9 1.495+40.335i 0 0.909 1;8 10 8.6+20.1i 0.000556i 1 0;]由各个变电所负荷功率可以计算出总功率为220MW，而发电厂一、二的总装机容量分别为300 MW和200 MW。