两机五节点网络潮流计算牛拉法设计说明书

1 两机五节点网络潮流计算牛拉法设计说明书 第一章 前言 流计算 潮流计算概述 潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供 依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据 。 具体表现在以下方面： (1)在电网规划阶段 ,通过潮流计算 ,合理规划电源容量及接入点 ,合理规划网架 ,选择无功补偿方案 ,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时 ,在预计负荷增长及新设备投运基础上 ,选择典型方式进行潮流计算 ,发现电网中薄弱环节 ,供调度员日常调度控制参考 ,并对规划、基建部门提出改进网架结构 ,加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算 ,用于 日运行方式的编制 ,指导发电厂开机方式 ,有功、无功调整方案及负荷调整方案 ,满足线路、变压器热稳定要求及电压质量要求。 2 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力 系统运行状态的实时监控中，则采用在线潮流计算。 此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。 流计算的要求 电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下： 1. 节点电压应满足 m i n m a x ( 1 , 2 , )i i U i n 式 (1 从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附 近。 点电压幅值必须按上述条件给定。因此，这一约束条件对 点而言。 2. 节点的有功功率和无功功率应满足 m i n m a xm i n m a xG i G i G iG i G i G Q 式 (1 3 点的有功功率和无功功率，以及 点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的 以及 应按上述条件进行检验。 3. 节点之间电压的相位差应满足 m a x| | | | | |i j i j i j 式 (1 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。 因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。 流计算的优势 电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算 。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用于电力系统稳态研究，安全估计或最优潮流等，实际电力系统的潮流计算主要采用牛顿 由 以轻松实现计算复杂的电力系统潮流分布。 介 含义是矩阵实验室（ 是美国 司于 1982推出的一套高性能的数值计算可视化软件，包括 程序、 态 4 系统仿真包和各种专业工具箱它集数值分析，矩阵计算，信号处理和图形显示于一体，构成一个方便的，界面友好的用户环境，具有极强大的计算功能和极高的编程效率，特别适合于科学计算、数值分析、系统仿真和信号处理等任务。 学设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率的解决工业计算问题，特别是 关于矩阵和矢量的计算。通过 言，可以用类似数学公式的方式来编写算法，大大降低了程序需要的难度别难并节省了时间，从而可把主要的经历集中在算法的构思而不是编程上。学习运用 算电力系统潮流分布是本次课程设计的重点，可以说，作为工科学生，会运用 解决工程问题已成为必须。到目前为止， 发展成为国际上最优秀的科技应用软件之一。它的功能十分强大，不仅仅可以实现计算潮流分布，还可以模拟仿真各式各样的数值系统，工程。这里将借助 完成用直角牛顿 分布。 言 1、 言相关指令 表 1理命令和函数 线帮助文件 入超文本说明 、 出 过 位函数和文件 行演示程序 制 搜索路径 5 管理变量和工作空间 出当前变量 出当前变量（长表） 磁盘文件中恢复变量 存工作空间变量 内存中清除变量和函数 理工作空间内存 阵的尺寸 量的长度 示矩阵或 与文件和操作系统有关的命令 变当前工作目录 录列表 除文件 取环境变量值 ! 执行 行 存 务 控制命令窗口 置命令行编辑 命令窗口 标置左上角 6 置输出格式 稿文件内使用的回显命令 命令窗口中控制分页输出 启动和退出 出 用 所执行的 启动 一般信息 统信息及 为 订购用户 服务程序的识别代号 说明书中未包含的新信息 本信息 操作符和特殊字符 + 加 减 * 矩阵乘法 .* 数组乘法 矩阵幂 . 数组幂 左除或反斜杠 / 右除或斜杠 7 ./ 数组除 量积 : 冒号 ( ) 圆括号 方括号 . 小数点 . 父目录 继续 , 逗号 ; 分号 % 注释 ! 感叹号 转置或引用 = 赋值 = = 相等 关系操作符 & 逻辑与 | 逻辑或 逻 辑非 辑异或 逻辑函数 查变量或函数是否存在 8 量的任一元为真，则其值为真 量的所有元为真，则其值为真 出非零元素的索引号 三角函数 弦 曲正弦 正弦 n ?解修正方程式，由 、 、 和 J 求各节点电压变量 、 ； i = 1 , 2 , ， n , i s() ) 2| | 、 | | ()ke m a x()a k k+1计算各节点电压新值：= + ， i = 1 ,2 , ,n= + ; i s( 1 )1 )()i i + 1求出 | | 、 | |()a x m a 线路功率sS ； ， 设节点电压 、 ； ， 置迭代次数对 节点：按式 计算 、对 节点：按式 计算 、置节点号按式 计算雅可比矩阵元素 、 、 、 、 、按式 计算雅可比矩阵元素 、 、 、 、 、雅可比矩阵 是否已全部形成，解修正方程式，由 、 、 和 求各节点电压变量 、 ； ， 解修正方程式，由 、 、 和 求各节点电压变量 、 ； ， 迭代是否收敛， 、 迭代是否收敛， 、 增大迭代次数， 计算各节点电压新值： ， 计算各节点电压新值： ， 增大节点号 求出 、 求出 、 计算平衡节点功率 和线路功率计算平衡节点功率 23 参数矩阵： 1,2,j;1,3,j;2,3,j;2,4,j;2,5,j;3,4,j;4,5,j; 源程序： =请输入节点总数为 :N=); nl=请输入支路数为： ); 请输入平衡节点为 :); pr=请输入迭代精度为 :); B1=请输入由支路参数形成的矩阵为 :); Y=); i=1:nl p=B1(i,1);q=B1(i,2); Y(p,q)=Y(p,q)B1(i,3); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./B1(i,3); Y(p,p)=Y(p,p)+1./B1(i,3); 导纳矩阵为 :Y=); %求出导纳矩阵 ;1;1;1;0;0;0;0;0;P=0;Q=0; 24 );k=0; G=);B=); m=1:N if m=; n=1:N pp=pp+ei(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)+fi(m)*(G(m,n)*fi(n)+B(m,n)*ei(n); qq=qq+fi(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)m)*(G(m,n)*fi(n)+B(m,n)*ei(n); P(m)=P(m) %有功功率的不平衡量 DQ(m)=Q(m) %无功功率的不平衡量 I(m)= (P(m)-Q(m)/ei(m)+fi(m)*i); H(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)+(m); N(m,m)=G(m,m)*ei(m)+B(m,m)*fi(m)+(m); J(m,m)=-G(m,m)*ei(m)-B(m,m)*fi(m)+(m); L(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)(m); *H(m,m);*N(m,m);*J(m,m);*L(m,m); 25 m=1:N if m=n=1:N if n=if n=m (m,n)=-B(m,n)*ei(m)+G(m,n)*fi(m); N(m,n)=G(m,n)*ei(m)+B(m,n)*fi(m); J(m,n)=-B(m,n)*fi(m)-G(m,n)*ei(m); L(m,n)=G(m,n)*fi(m)-B(m,n)*ei(m); *H(m,n);*N(m,n);*J(m,n);*L(m,n); 雅克比矩阵 ); J); 26 m=1:N if m=W(2*DP(m);*DQ(m); %形成不平衡量的矩阵 J %解修正方程式 pr=Y); 第 M 次修正方程的解 ); Y); n=1:N if n=fi(n)=fi(n) +* ei(n)=ei(n) +* %计算新值 u=ei+fi*i; 节点电压的第 C(k)次近似值 :); u); 各点的电压实部 位 :V)为 (节点号从小到大排列 ):); 各点的电压虚部 位 :V)为 (节点号从小到大排列 ):); k=k+1; 迭代次数 :); 27 k); m=1:N I(m)=Y(1,m)*u(m); 平衡节点的功率 ); S1=u(1)*)%计算平衡节点的功率 m=1:N n=1:N S(m,n)=u(m)*(u(m)u(n)*Y(m,n);%计算各支路功率 各支路功率 );) %结束 运行结果 导纳矩阵为 :Y= 0 0 0 0 0 28 雅克比矩阵 0 0 0 0 0 0 0 0 M 次修正方程的解 点电压的第 C(k)次近似值 : 29 点的电压实部 位 :V)为 (节点号从小到大排列 ): 点的电压虚部 位 :V)为 (节点号从小到大排列 ): 0 代次数 : 1 雅克比矩阵 30 0 0 0 0 0 0 0 0 M 次修正方程的解 点电压的第 C(k)次近似值 : 31 点的电压实部 位 :V)为 (节点号从小到大排列 ): 点的电压虚部 位 :V)为 (节点号从小到大排列 ): 0 代次数 : 2 雅克比矩阵 0 0 32 0 0 0 0 0 0 M 次修正方程的解 点电压的第 C(k)次近似值 : 33 点的电压实部 位 :V)为 (节点号从小到大排列 ): 点的电压虚部 位 :V)为 (节点号从小到大排列 ): 0 代次数 : 3 雅克比矩阵 0 0 34 0 0 0 0 0 0 M 次修正方程的解 点电压的第 C(k)次近似值 : 35 点的电压实部 位 :V)为 (节点号从小到大排列 ): 点的电压虚部 位 :V)为 (节点号从小到大排列 ): 0 代次数 : 4 平衡节点的功率 36 支路功率 0 0 0 0 0 0 0 0 0 0 参考文献 1、电力系统稳态分析，陈珩，中国电力出版社， 2007，第三版， 1272、电力系统分析课程实际设计与综合实验，祝书萍，中国电力出版社， 2007，第一版 ,1133、关于 谢 感谢刘景霞老师帮助我们学习电力系统稳态分析这本书并和杨培宏老师耐心地给我们辅导课程设计书！对我们以后的学习有了更大的帮住！由于以前没学过开始一直没