三机九节点潮流暂态MATLAB仿真.docVIP

电力系统仿真 课程报告 院 系： 自动化学院 专 业：电力系统及其自动化 学 号： 姓 名： 时 间： 研究对象 三机九节点系统 图1.1 WSCC-9系统模型 图1.1是一个三机九节点的系统阻抗图，给出的阻抗参数都是以MVA为基准的标幺值。包括三台发电机，双绕组变压器，）三个。将对该系统动态过程进行相应的仿真分析。参数 参数 节点类型，节点分为已知参数如下表： .1 节点已知参数幅值角度有功无功有功无功.040 0 0.7160 0.2705 0 0 2 PV 1.025 1.6300 0.0665 0 0 3 PV 1.025 0.8500 -0.1086 0 0 4 PQ 0 0 5 PQ 1.2500 0.5000 6 PQ 0.9000 0.3000 7 PQ 0 0 8 PQ 1．0000 0.3500 9 PQ 0 0 上表中发电机有功无功负荷的有功无功均为以MVA为基准时的标幺值。 参数 .2 支路参数 首节点 末节点 电阻 电抗 电纳一半.0100 0.0850 0.0880 4 6 0.0170 0.0920 0.0790 5 7 0.0320 0.1610 0.1530 6 9 0.0390 0.1700 0.1790 7 8 0.0085 0.0720 0.0745 8 9 0.0119 0.1008 0.1045 1 4 0.0000 0.0576 0.0000 2 7 0.0000 0.0625 0.0000 3 9 0.0000 0.0586 0.0000 上表中所有的均为，对于变压器支路最后三行三台参数，已经变压器等效电抗并直接在表格中给出。 参数 发电机，采用二阶经典模型系统如下假设： 的机械功率保持恒定； 阻尼效应 负荷采用恒阻抗模型。 .3 发电机参数 1 0.1460 0.0608 23.64 2 0.8958 0.1198 6.40 3 1.3125 0.1813 3.01 以上阻抗参数标幺值表示，额定转速下存储的能量（）以MVA为基准的标幺值。 求解的动态过程 运行时，s时刻在线路– 7末端靠近母线发生三相短路故障，故障持续周期（.083s）s – 0.083s时的动态过程然后– 7线路，.083s – 2s时的系统动态过程[1]。 仿真过程 流程 .1 仿真过程总流程 流程 计算采用——拉夫逊，的良好的收敛性，其解算过程理论依据如下 极坐标形式的功率方程 将其改写为增量的形式 有功功率： 无功功率： 以上两式经泰勒级数展开可线性化为： 上式可简写为： 进一步改写为： 其中为雅可比矩阵在每次潮流计算后通过式修正量，根据对x的值进行修正，并再次潮流，精度达到给定的精度限制值。 .2 潮流子程序流程图 计算得到节点的电压幅值和角度，得到每个节点的功率和无功功率，是进一步分析动态过程基础。 发电机初态 动态分析时，应将发电机和负荷用相应的模型等效本文等效为二阶模型将负荷等效为恒阻抗负荷。 .3 发电机等效模型 中具有三台发电机，因此会引入三个内节点节点导纳矩阵Y将增广到[2]。可表示为： 中： 和中： 中： 将网络等效在发电机内节点构成网络中，可以得到一个的节点导纳矩阵，矩阵可由以下关系解： 于故障和故障后的降阶导纳矩阵，利用式计算出对应的降阶节点导纳矩阵对于故障中的节点矩阵，的增广矩阵中去掉故障母线所在的那一行和那一列，利用下式计算降阶矩阵。 初态中包含的电压幅值和功角，根据图.3，可以求解发电机初始状态。 列写发电机动态方程 采用经典模型，其动态方程为： 发电机共有六个状态量定义 得到： 状态 龙格——库塔积分方法，求解对发电机状态方程积分，求解发电机的状态。认为故障过程中，发电机内电压的幅值是不变的，只有功角改变发电机的机械功率是不变的其值等于故障前的发电机电磁功率。 程序分析 图3.1 潮流计算结果矩阵 .2 降阶矩阵 发电机初态 .3 发电机初态 过程 .083秒时刻的发电机 图3.4 0.083s时的发电机状态 时刻的发电机状态 .5 2s时刻的发电机状态 和转速的变化曲线 图3.6 功角变化曲线 图3.7 功角差变化曲线 图3.8 转速变化曲线 分析 三台发电机、九母线的系统，在– 2s的时间内，给出了摇摆方程的数值积分。图.6给出了功角变化曲线，图.7给出了功角差的变化曲线。从中可以看出