三机九节点电力系统仿真matlab

电力系统仿真作业------------ 三机九节点电力系统三机九节点电力系统 暂态仿真暂态仿真 学院能源与动力工程学院学院能源与动力工程学院 专业电力系统及其自动化专业电力系统及其自动化 学号学号 姓名于永生姓名于永生 导师导师 授课教师授课教师 目录目录 一、一、 概述概述.1 二、二、 课程主要任务课程主要任务 .1 1. 1. 2. 2. 3. 3. 4. 4. 5. 5. 6. 6. 系统数据系统数据 1 潮流计算潮流计算 2 负荷等效和支路简化负荷等效和支路简化 .4 求解电磁功率求解电磁功率 .5 求解运动方程求解运动方程 .5 程序清单程序清单 7 1. 2. 3. 4. 主程序.7 极坐标转换成直角坐标函数pol2rectV,del16 直角坐标转换成极坐标函数rect2polZ16 求解微分方程所用的得到微分量的函数Gen_fwt,X,Y_Gen,E,Pm0,Tj.16 三、三、 课程总结及心得体会课程总结及心得体会 .16 四、四、 参考文献参考文献 .17 一、一、概述概述 在动态稳定分析中，系统由线性化的微分方程组和代数方程组描写，并用经典的或现 代的线性系统理论来进行稳定分析，分析可以在时域或频域进行。当用计算机和现代线性 系统理论分析时，常把系统线性化的微分方程组和代数方程组消去代数变量，化为状态方 程形式，并广泛采用特征分析进行稳定分析。 电力系统是由不同类型的发电机组、多种电力负荷、不同电压等级的电力网络等组成 的十分庞大复杂的动力学系统。其暂态过渡过程不仅包括电磁方面的过渡过程，而且还有 机电方面的过渡过程。 由此可见， 电力系统的数学模型是一个强非线性的高维状态方程组。 在动态稳定仿真中使用简单的电力系统模型，发电机用三阶模型表示。 二、二、课程主要任务课程主要任务 本次课程主要应用 P. M. Anderson and A. A. Fouad编写的Power System Control and Stability 一书中所引用的 Western System Coordinated Council WSCC三机九节点系统模型。 1.1. 系统数据系统数据 其中，节点数据如下 节点数据 节点电压电压发电机发电机负荷负荷节点 号幅值相角有功无功有功无功类型1PQ 2PV 3平 衡 N[11.0400.71640.2705003 21.025 01.630.0665002 31.025 00.85-0.1086002 41000001 510001.250.51 610000.90.31 71000001 8100010.351 91000001]; 其中，支路数据如下 线路数据 首端末端电阻电抗电纳（1/2变压器非标准变比 L[450.010.0850.0881 460.0170.0920.0791 570.0320.1610.1531 690.0390.170.1791 780.00850.0720.07451 890.01190.10080.10541 1400.057601 2700.062501 3900.058601]; 发电机数据如下 发电机母线XdXdTd0XqXqTq0’TjXf Ge[110.14600.06088.960.09690.0969047.280.0576 220.89580.11986.000.86450.19690.53512.800.0625 331.31250.18138.591.25780.25000.6006.020.0585]; 系统电路结构拓扑图如下 图 1 WSCC 3 机 9 节点系统 （所有参数以 100MVA 为基准值的标幺值） 2.2. 潮流计算潮流计算 首先进行潮流计算，采用牛顿拉夫逊迭代法，电力系统潮流计算是电力系统运行和规 划中最基本和最经常的计算，其任务是在已知某些运行参数的情况下，计算出系统中全部 的运行参数，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平 衡节点除外） ，可以根据网络结构形成节点导纳矩阵， 然后由节点导纳矩阵和网络拓扑结构 列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计 算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述 功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未 知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这 样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差 方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值， 将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡 量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿拉夫逊算法修正方程 W -JΔV 其中 W 是节点不平衡量向量，包括有功，无功，电压；J 是雅克比矩阵；ΔV 是节点 电压修正量。 令 V i  e i  jf i ;Y ij  G ij  jB ij 则直角坐标形式的功率不平衡量方程为 PQ 节点 ， P i  P is e i G ijej  B ij f j  f j G ij f j  B ij e j  0 j1j1 nn nn Q i  Q is  f i G ijej  B ij f j  e j G ij f j  B ijej  0 j1j1 PV 节点 P i  P is e i G ijej  B ij f j  f j G ij f j  B ij e j  0 j1j1 nn V i 2V is 2V i 2V is 2e i 2 f i 2 极坐标形式的功率不平衡量方程 P i  P is V i V j G ij cos ij  B ij sin ij  0 j1 n n Q i Q is V i V j G ij sin ij  B ij cos ij  0 j1 雅可比矩阵 J 各元素的表达式 H J  M N  L 当 j≠i 时 H ij  P i B ijei G ij f i f j N ij  P i G ijei  B ij f i e j Q i B ij f i G ijei f j M ij  L ij  Q i G ij f i  B ijei e j 当 ji 时 H ii  P i B iiei G ii f i b i f i P i G iiei  B ii f i a i e i Q i G iiei  B ii f i a i f i N ii  M ii  L ii  Q i B iiei G ii f