《电力系统分析潮流计算课程序的设计和MATLAB程序的设计说明》由会员分享，可在线阅读，更多相关《电力系统分析潮流计算课程序的设计和MATLAB程序的设计说明（23页珍藏版）》请在金锄头文库上搜索。
1、 电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算 学 院 电气工程学院 专业班级 学生姓名 学 号 班内序号 指导教师 房大中 提交日期 2015年05月04日 . .目录一、程序设计目的2二、程序设计要求4三、13节点配网潮流计算43.1主要流程4第一步的前推公式如下1-1-1-5：4第二步的回代公式如下1-61-9：43.2配网前推后代潮流计算的原理53.3配网前推后代潮流计算迭代过程53.3计算原理6四、计算框图流程7五、确定前推回代支路次序7六、前推回代计算输入文件8主程序：8输入文件清单：9计算结果：10数据分析：10七、配电网潮流计算的要点11八、自我总结11九、参考
2、文献12附录一 MATLAB的简介12. .一、程序设计目的开式网络潮流计算：配电网的结构特点呈辐射状,在正常运行时是开环的；配电网的潮流计算采用的方法是前推回代法,本程序利用前推回代法的基本原理、收敛性。在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压的要求。 在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。 正常检修及特殊运行方式下的潮流计算,用于日运行方式的编
3、制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求。及电压质量要求。 预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 图1 13节点配电网结构图表1 系统支路参数支路RXB/2 S123.3673.6850.0232.3562.5410.0341.1451.280.0454.5245.040.0260.8561.140.0672.7452.9650.0283.7434.2510.0892.2372.7560.03104.144.6960.03111.3281.7630.011122.4362.8660.04133.5213.9660.
4、0表2 系统负荷参数节点编号节点类型节点初始电kVPiMVAQiMVA1根节点10.4002中间节点10.00.03420.03013中间节点10.00.06930.06424中间节点10.00.08450.07635叶节点10.00.02950.02616中间节点10.00.04740.04097叶节点10.00.11760.09578中间节点10.00.09460.08579叶节点10.00.09160.085910叶节点10.00.02710.022911中间节点10.00.06960.064312叶节点10.00.06760.057913叶节点10.00.02980.0242二、程序设
5、计要求1. 看懂前推回代法计算程序；2. 报告叙述计算原理及计算流程；3. 绘制计算流程框图；4. 确定前推回代支路次序广度优先,或深度优先并编写前推回代计算输入文件,然后进行潮流计算；5. 整理专利计算结果；6. 总结配电网潮流计算的要点；三、设计内容1.根据电力系统网络推导电力网络数学模型,写出节点导纳矩阵；2.赋予各节点电压变量直角坐标系形式初值后,求解不平衡量；3.形成雅可比矩阵；4.求解修正量后,重新修改初值,从2开始重新循环计算；5.求解的电压变量达到所要求的精度时,再计算各支路功率分布、功率损耗和节点 6.上机编程调试；连调；7.计算分析给定系统潮流分析并与手工计算结果作比较分析
6、。8.准备计算机演示答辩,书写该课程设计说明书必须计算机打印四、13节点配网潮流计算4.1牛拉法的原理及其基本方程牛顿迭代法Newtons method又称为牛顿-拉夫逊拉弗森方法Newton-Raphson method,它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。设r是f = 0的根,选取x0作为r初始近似值,过点x0,f做曲线y = f的切线L,L的方程为y = f f,求出L与x轴交点的横坐标 x1 = x0-f/f,称x1为r的一次近似值。过点x1,f做曲线y = f的切线,并求该切线与x轴的横坐标 x2 = x1-f/f,称x2为r的二次近似值。重复以上过程,
7、得r的近似值序列,其中x=xfx/fx,称为r的n+1次近似值,上式称为牛顿迭代公式。解非线性方程f=0的牛顿法是把非线性方程线性化的一种近似方法。把f在x0点附近展开成泰勒级数 f = f+f+2*f/2! + 取其线性部分,作为非线性方程f = 0的近似方程,即泰勒展开的前两项,则有f+f=f=0 设f0则其解为x1=x0f/f 这样,得到牛顿法的一个迭代序列：x=xfx/fx。4.2 PQ分解法的原理及其基本方程PQ分解法的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻值,则系统母线电压幅值的微小变化对母线有功功率的改变影响很小。同样,母线电压相角的少许改变 ,也不会引起
8、母线有功功率的明显改变 ,因此,节点功率方程在用极坐标形式表示时,它的修正方程式可简化为 2-1这就是把2n-1阶的线性方程组变成了两个n-1阶的线性方程组,将P和Q分开来进行迭代计算,因而大大地减少了计算工作量。但是H、L在迭代过程中仍然在不断的变化而且又都是不对称矩阵。对牛顿法的进一步简化,即把式2-1中的系数矩阵简化为在迭代过程中不变的对称矩阵。在一般情况下,线路两端的电压相角 是不大的不超过1020,因此,可以认为 2-2此外,与系统各节点无功功率相应的导纳远小于该节点自导纳的虚部,即因而 2-3考虑到以上关系,式2-1的系数矩阵中的个元素可以表示为 而系数矩阵H和L则可以分别写成：
9、将式2-4和2-5代人式中,得到用 和 分别左乘以上两式,便得 这就是简化了的修正方程式,它们也可以展开写成 在这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部,因而系数矩阵是对称矩阵,且在迭代过程中保持不变。这就大大减少了计算工作量。用极坐标表示的节点功率增量为 式2-8、2-9和2-10构成了PQ分解法迭代过程的基本方程式。4.3配网前推后代潮流计算的原理前推回代法在配电网潮流计算中简单实用,所有的数据都是以矢量形式存储,因此节省了大量的计算机内存,对于任何种类的配电网只要有合理的 R/X 值,此方法均可保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下,算法的收敛阶数
10、越高,算法的稳定性越差,前推回代法的收敛阶数为一阶,因此它也具有较好的稳定性。比较而言,前推回代法充分利用了网络呈辐射状的结构特点,数据处理简单,计算效率高,具有较好的收敛性,被公认是求解辐射状配电网潮流问题的最佳算法之一。4.4配网前推后代潮流计算迭代过程配网潮流前推回代潮流算法第n+1步的迭代公式如下：节点i的前推计算公式为由图1可知节点i+1的注入有功功率和无功功率分别为：NBNB-1P= PL+ LPQ= QL+ LQ其中：i=1,2,。NB-1,NB为节点数,PL为第i节点所带负荷有功功率；QL为第i节点所带负荷无功功率；LP为第i条线段上的有功功率损耗；LQ为第i条线段上的无功功率
11、损耗。收敛条件以前后两次迭代的电压偏差作为迭代收敛条件,式表明节点电压幅值最大偏差小于设定阈值,即认为迭代收敛,则迭代过程结束。max|Vk-Vk-1|.在开始计算的时候先把各节点上的LP、LQ分别设为零,利用式计算出P、Q,分别代入式计算各节点电压幅值,作为初始化的节点电压值.然后把此电压值代入公式来修正线路的损耗LP、LQ,用经过修正的线路损耗求新的节点注入功率,再根据节点注入功率求新的节点电压,根据前后两次迭代的电压偏差是否小于设定阈值判断是否收敛,如果满足,则输出结果,如果不满足,必须重复上述过程直到满足收敛条件为止。4.5计算原理第一步,从与叶节点联系的支路开始,该支路的末端功率即等于叶节点功率,利用这个功率和对应的节点电压计算支路功率损耗,求得支路的首端功率。当以某节点为始节点的各支路都计算完毕后,便想象将这些支路都拆去,是的该节点成为新的叶节点,其节点功率等于原有的负荷功率与该节点为始节点的各支路首端功率之和。于是计算便可延续下去,直到全部支路计算完毕。第二步,利用第一步所得的支路首端功率和本步骤刚算出的本支路始节点的电压对电源点为已知电压,从电源点开始逐条支路进行计算,求得各支路终节点的电压。在迭代计算开始之前,要先处理好支路的计算顺序问题。介绍两种确定之路计算顺序的方法。第一种方法是,按与叶节点联接的之路顺序
