电力系统潮流计算课程设计论文

1摘要：潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压，包括电压的幅值和相角，各元件流过的功率，整个系统的功率损耗等一系列系统中的潮流数据。潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。将N-R法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性，稀疏性及节点编号顺序优划等技巧，使N-R法在收敛性，占用内存，计算速度等方面的优点都超过了阻抗法。同时本文也应用了电力系统潮流计算仿真软件DDTR与利用程序计算的结果进行比较，使计算的结果更加准确。利用成形的程序对系统中出现的各种情况，例如负荷的变化以及线路上所发生的变化进行计算，并对母线上不满足范围的数据进行调控，使得系统处于一个较稳定的状态。关键词牛顿-拉夫逊法（Newton-Raphson）潮流计算与DDRT仿真潮流计算结果分析与调控2一、问题重述（即课程设计任务书）1、系统图如图1所示：2、发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为（300MW），母线3为机压母线，机压母线上装机容量为（100MW），最大负荷和最小负荷分别为50MW和20MW；发电厂二总装机容量为（200MW）。3、变电所资料：(一)变电所1、2、3、4低压母线的电压等级分别为：35KV10KV35KV10KV(二)变电所的负荷分别为：50MW50MW40MW60MW变电所1变电所2变电所3变电所435kV母线10kV母线35kV母线10kV母线一次侧电压220kV一次侧电压220kV母线1母线3母线2线路长为100km线路长为50km线路长为90km线路长为100km2*QFQ-50-22*QFS-50-2TQN-100-22*TQN-100-2。。。。。。。。。。。。。线路长为80km。。。。。。。。。。。。。电厂一电厂二3(三)每个变电所的功率因数均为cosφ=0.85；(四)变电所1和变电所3分别配有两台容量为75MVA的变压器，短路损耗414KW，短路电压（%）=16.7；变电所2和变电所4分别配有两台容量为63MVA的变压器，短路损耗为245KW，短路电压（%）=10.5；4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为17.0，单位长度的电抗为0.402，单位长度的电纳为S-610*2.78。通过以上给出数据得到该系统的等值电路，并通过计算比较与分析得到系统的潮流分布，再由不同的的负荷变化调控系统的潮流分布。二、问题分析（一）牛顿—拉夫逊法概要首先对一般的牛顿—拉夫逊法作一简单的说明。已知一个变量X函数为：0)(Xf到此方程时，由适当的近似值)0(X出发，根据：,......)2,1()()()()()()1(nXfXfXXnnnn反复进行计算，当)(nX满足适当的收敛条件就是上面方程的根。这样的方法就是所谓的牛顿—拉夫逊法。这一方法还可以做下面的解释，设第n次迭代得到的解语真值之差，即)(nX的误差为时，则：0)()(nXf把)()(nXf在)(nX附近对用泰勒级数展开0......)(!2)()()()(2)()()(nnnnXfXfXfXf上式省略去2以后部分0)()()()(nnXfXf4)(nX的误差可以近似由上式计算出来。)()()()(nnXfXf比较两式，可以看出牛顿—拉夫逊法的休整量和)(nX的误差的一次项相等。用同样的方法考虑，给出n个变量的n个方程：0),,,(0),,,(0),,,(21212211nnnnXXXfXXXfXXXf对其近似解1X得修正量1X可以通过解下边的方程来确定：nnnnnnnnnnnXXXxfxfxfxfxfxfxfxfxfXXXfXXXfXXXf2121222121211121212211),,,(),,,(),,,(式中等号右边的矩阵nnxf都是对于nXXX,,,21的值。这一矩阵称为雅可比（JACOBI）矩阵。按上述得到的修正向量nXXX,,,21后，得到如下关系nnnXXX这比nXXX,,,21更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反复计算满足1112121111,,,maxnnnnnnnnXXXXXX为预先规定的小正数，1nnX是第n次迭代nX的近似值。（二）牛顿法的框图及求解过程1、用牛顿法计算潮流时，有以下的步骤：（1）给这各节点电压初始值)0()0(,fe；（2）将以上电压初始值代入公式，求修正方程的常数项向量5)0(2)0()0()(,,VQP；（3）将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。（4）解修正方程式)0()0(,fe；（5）修正各节点电压)0()0()1(eee，)0()0()1(fff；（6）将)1(e，)1(f在带入方程式，求出)1(2)1()1()(,,VQP；（7）检验是否收敛，即)()(,maxkikiQP（8）如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回（2）进行下次迭代计算，直到收敛为止。2、程序框图如下6三、问题求解（一）等值电路的计算过程电压是衡量电力系统电能质量的标准之一。电压过高或过低，都将对人身及其用电设备产生重大的影响。保证用户的电压接近额定值是电力系统调度的基本任务之一。当系统的电压偏离允许值时，电力系统必须应用电压调节技术调节系统电压的大小，使其维持在允许值范围内。本文经过手算形成了等值电路图，并编写好了程序得出节点电压标幺值，使其满足所要求的调整范围。我们首先对给定的程序输入部分作了简要的分析，程序开始需要我们确定输入节点数、支路数、平衡母线号、支路参数矩阵、节点参数矩阵。（1）为了保证整个系统潮流计算的完整性，我们把凡具有母线及发电机处均选作节点，这样，共有10条母线以及4部发电机组即选定14个节点，我们确定电厂一母线3上的发电机作为平衡节点，节点号为①，其它机组作为PV节点，节点号为②③④，其余节点均为PQ节点，节点号见等值电路图。（2）确定完节点及编号后，各条支路也相应确定了，网络中总计有13条支路，我们对各支路参数进行了计算。根据所给实际电路图和题中的已知条件，有以下公式计算各输电线路的阻抗和对地支路电容的标幺值和变压器的阻抗标幺值。选择电压基准值为220BUKV和功率基准值100BSMVA所以484BBBUZS7（二）等值电路及参数：利用DDRT软件可得计算所需系统图如图3所示：图3（下端发电机组编号依次为①②③④）（三）潮流计算1、潮流分布根据任务书中的条件，四个变电所的总负荷为200MW,两个发电厂的总装机容量为500MVA(未包括母线3上的损耗)，那么就需要调整发电机的发出功率以免发电机发出负有功造成浪费以及经济方面的不合理。利用MATLAB软件，8分析不同的发出功率，以上面的条件为依据，经过比较与选取，可得一组较好的发出功率：①号发电机的发出功率50MW100KM②号的发出功率为50MW③号的发出功率为50MW④号发电机的发出功率为80MW。将计算所需数据再重新输入到潮流计算程序中，得到潮流计算结果，程序及结果如附录1所示。同时将所得结果与DDRT潮流模拟结果进行比较，发现两组计算结果相差为1%左右，可知此时计算的计算结果较准确。下图经MATLAB计算数据所得的四组数据的结果坐标图：0510150.9511.051.11.15节点号电压标幺值051015-10-505节点号电压角度1234567891011121314-0.4-0.200.20.4节点注入有功12345678910111213-0.200.20.40.6支路首端无功9下图为经DDRT软件仿真得到的潮流分布节点号1234567891011121314电压相01.22001.99743.2813-0.9710.9299-1.968-1.439-1.313-0.872-5.482-3.995-4.040-3.95010角72747938111节点的相角最大，说明在11节点处的无功功率较其它节点大。在该潮流中，各条输电线上的线损都很大，可以满足电能输送的经济性，若考虑线损范围并进行有功调控还需系统传输合理范围。各节点电压均在范围之内，特别是四个变电所低压侧电压均在任务书要求范围之内。2、系统4个变电所的负荷同时以2%的比例增大对于负荷的改变，在潮流计算中重新输入数据进行计算得到了在新的负荷下的系统潮流分布（程序相同，详细结果见附录）。如下给出了结果的坐标图以及仿真图0510150.9511.051.11.15节点号电压标幺值051015-10-505节点号电压角度1234567891011121314-0.4-0.200.20.4节点注入有功12345678910111213-0.200.20.40.6支路首端无功11节点号1234567891011121314电压相角01.04281.82073.0282-1.15030.6753-2.1645-1.6273-1.5335-1.1177-5.7623-4.2407-4.324-4.2672当负荷增大时，相角相应也增大。同时整个潮流也发生了相应的变化。123、4个变电所的负荷同时以2%的比例下降上图为负荷端以2%的比例下降时得到的系统潮流分布。在该潮流中，变电所3的电压为36.0108KV,超过了35~36KV的范围，所以需要对节点13的电压进行调控，即调整变压器的变比。在前面的负荷条件下变压器的变比均为1，现在将变比调为1.025，在调整后得到新的潮流分布，与用程序得到的数据比较基本一致。变电所编号123413调整前35.518910.35136.010810.2457调整后35.517910.350735.07710.2453单位：KV下图为调整后的新的潮流分布。144、1号和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升。上图给出的是利用软件的仿真图，利用程序得到的结果如下：节点编号1234567891011121314电压值231231231232218.2114220.4426214.1849215.8758216.2921226.337036.022235.00636.66736.5743电压01.232.013.33-0.90.97-1.9-1.4-1.2-0.8-5.3-4.0-4.0-3.815相角854045369633538582627063986953715、节点5与节点11之间的双回输电线路改为单回。利用仿真软件得到的潮流图为：当将节点5与节点11之间的双回线改为单回线后，节点11上的电压降低。需要进行无功补偿。补偿前后的结果为：变电所编号1234补偿前电压34.611310.33435.908510.2165补偿后电压35.035310.390535.954410.226316补偿后的潮流分布图：（6）依次断开50,90,80,100KMKMKMKM的双回线的一条变为单回线，进行潮流计算，并依照所给出的变电所侧的电压限制条件进检验。同时利用MATLAB运用程序进行潮流计算得到下表的潮流电压数据。而后给出几种情况利用仿真软件得到的系统潮流分图1750KM线路双回线变单回线：程序结果：节点号1234567891011121314电压值231231231231218.1555220.4243214.1270213.3927216.2475216.3036