电力系统潮流计算1-概念方程及计算方法.ppt

* 电力系统潮流计算（1）概念、方程及算法 华北电力大学电气与电子工程学院 孙英云 手机Email: sunyy@ncepu.edu.cn 办公室：教五 C204 * 问题 什么是潮流计算？ 什么是潮流？ 什么是计算？ 为什么要进行潮流计算？ 原因：电力系统状态不可直接测量 潮流计算结果和电力系统运行状态之间关系 电力系统运行状态有什么用？ 如何进行潮流计算？ * 潮流计算发展简史 史前时代 手算、交流模拟台 50年代Y矩阵法（Gauss迭代法） 内存需求量小，收敛性差； 60年代初Z矩阵法 收敛性好，内存占用大； 60年代Newton－Raphson法； Tinney稀疏矩阵技术、节点优化编号； 1974年B Stott 提出快速分解法（Fast Decoupled Load Flow）； * 简单电力系统等值电路（实例） 发电机 输电线路 配电线路 降压变压器 负荷 降压变压器 升压变压器 G T1 T2 T3 L1 L2 K2ZT2 Z210 Z220 ZL2 YL2/2 YL2/2 K3ZT3 Z310 Z320 ZL1 YL1/2 YL1/2 PD+jQD K1ZT1 Z110 Z120 G * 电力系统稳态数学模型 发电机 出力可调，机端电压可控：PV或平衡节点 P=const、U=const P=const、Q=const 电力网络 节点导纳阵（Y） 负荷 恒功率模型（PQ节点） P=const，Q=const * 潮流计算数学模型节点功率平衡方程 电力网络—电路网络 节点电压方程 节点功率平衡方程： 将其代入可得： 即： 所有节点的功率平衡方程 问题：公式里的功率是什么功率？ 问题：公式里的电压和电流分别是 什么电压和电流？ * 直角坐标功率平衡方程 如果将节点电压用直角坐标表示，即令 则有： * 极坐标功率平衡方程 如果将节点电压用极坐标表示，即令 则有： * 从节点功率平衡方程到潮流方程——节点类型的划分 对于电力系统来讲，每个节点有四个运行变量（电压×2，功率×2），两个功率平衡方程（有功、无功） 负荷节点 负荷由需求决定，一般不可控，PQ节点 发电机节点 发电机励磁控制电压不变，PV给定，PV节点 考虑系统网损 电压、相角给定，平衡节点 * 从节点功率平衡方程到潮流方程——节点类型的划分 一个N个节点的电力网络，若选第N个节点为平衡节点，则剩下n（n=N-1）中有r个节点是PV节点，则PQ节点个数为n-r个。 已知量为：平衡节点的电压；除平衡节点外所有节点的有功注入量；PQ节点的无功注入量；PV节点的电压辐值 直角坐标下和极坐标下有不同的处理方法 * 直角坐标下潮流方程 直角坐标下待求变量 直角坐标下功率方程 * 直角坐标下潮流方程 直角坐标潮流方程的已知量和待求量？ * 极坐标潮流方程 极坐标潮流方程的已知量和待求量？ * 潮流方程的解法 潮流方程是一组高维非线性方程组 所有能用于求解非线性方程组的方法都可以用于求解潮流方程 Gauss法（简单迭代法） Newton法（包括其变形算法） 割线法 拟牛顿法 …… * 以Gauss法为基础的潮流方程解法 待求方程 高斯迭代法 当矩阵的谱半径小于1时收敛，谱半径越小，收敛性越好 * 以如下非线性方程为例进行说明 写成gauss法形式为？ 如果取初值为 X(1)=0.75 X(2)=0.8125 X(3)=0… X(100)= 0.9906925 * 基于节点导纳矩阵的高斯迭代法（P176） 令 则有 * 高斯法的讨论 高斯法可分为基于节点导纳阵的高斯法和基于阻抗阵的高斯法两种 高斯法的改进 高斯-赛德尔法 高斯法的PV节点处理较为困难 具体可参见 Kusic G L. Computer-aided power systems analysis. Prentice Hall, 1986 * 牛顿-拉夫逊法潮流计算 牛顿法的历史 牛顿法基本原理 对于非线性方程 给定初值 用Talor级数展开，有： 忽略高阶项，则有 * 牛顿-拉夫逊法潮流计算 牛顿法的几何意义 * 以如下非线性方程为例进行说明 写成牛顿法形式为？ 如果取初值为 X(1)=0.75 X(2)=0.875 X(3)=0.9375 X(4)=0.96875 X(5)=0.984375 X(6)=0.9921875 … X（20）=0.9999999 * 牛顿-拉夫逊法潮流计算 牛顿法计算流程 1 初始化，形成节点导纳阵，给出初值 2