电力系统-潮流计算技术报告.ppt

* 电力系统分析 稳态潮流计算详解 电气学院 吴俊勇 求解单变量非线性方程组的牛顿－拉夫逊法 单变量非线性方程 假定一个近似解 和误差 泰勒级数展开 忽略高次项 （修正方程） 得到一次近似解 一般迭代形式 收敛判据 求解多变量非线性代数方程组的牛顿－拉夫逊法 多变量非线性代数方程组 选定初始值，泰勒级数展开，忽略高次项 写成矩阵形式 得到一次近似解 一般形式 或 （修正方程） 收敛判据 节点电压方程 复功率方程 统一潮流方程 直角坐标下的牛顿－拉夫逊潮流计算 节点电压 节点导纳矩阵元素 带入统一潮流方程，并按实部和虚部展开 假设第1至m号节点为PQ节点，第i节点的给定有功和无功为 和 ，则节点方程 假设第m＋1至n－1号节点为PV节点，第i节点的给定有功和电压为 和 ， 则节点方程 第n节点为平衡节点，电压已给定，不参与迭代。 因此，共有2（n－1）个方程， 2（n－1）个待求变量。 修正方程 m个PQ节点 n－m－1个PV节点 J：雅克比矩阵 修正方程也可写成方块矩阵的形式 节点导纳矩阵 雅克比矩阵 像把伞！ 对PQ节点 对PV节点 PQ 节点 PV 节点 维数：2(n-1)X2(n-1) 对角线元素， ： 非对角线元素， ： 雅克比矩阵的特点： 各元素都是节点电压的函数，在迭代过程中不断改变； 雅克比矩阵的非对角子块Jij中的元素只用到了节点导纳矩阵中非对角元素Yij,如果Yij＝0，则Jij＝0。可见，雅克比矩阵与节点导纳矩阵一样，具有稀疏性； 节点导纳矩阵的对角线元素不为零，所以雅克比矩阵的对角子块的各元素也不为零； 节点导纳矩阵具有对称性，但雅克比矩阵的元素和子块都不具有对称性； 直角坐标下潮流计算的流程 ①输入电网络的原始数据，用追加支路法形成节点导纳矩阵； ②选定PQ节点、PV节点和平衡节点，给定各节点电压的初值 ，迭代 次数k＝0； ③根据给定的节点电压，计算各类节点的不平衡量 ； ④校验收敛条件：若 ，则转入⑨； ⑤计算雅克比矩阵各子块的元素，形成修正方程； ⑥用高斯消去法或三角分解法求解修正方程，得到各节点电压的 修正量 ； ⑦修正各节点电压： ； ⑧迭代次数k=k+1,返回第③步，继续迭代； ⑨迭代结束后，计算网络中的各支路功率分布、平衡节点的功率和网损。 输电线路功率的计算公式 推荐一种电网络原始数据的格式： 对一个6节点、6支路、3台发电机、3个负荷的电网络： network.dat： #节点数 6 #支路数 6 #各支路参数：起点编号，终点编号，电阻，电抗，电导，电纳 1,2,0.001,0.034,0.0024,0.00005 …… #各节点参数：节点编号，类型，注入有功，注入无功，电压幅值，电压相位 #类型：1＝PQ节点，2＝PV节点，0＝平衡节点 1,0,0.0,0.0,1.0,0.0 2,1,50.0,35.0,0.0,0.0 3,2,30.0,0.0,1.02,0.0 4,1,-18.0,-16.0,0.0,0.0 5,1,-30.0,-26.0,0.0,0.0 6,1,-50.0,-45.0,0.0,0.0 极坐标下的潮流计算 节点坐标用极坐标表示 带入统一潮流方程，按实部和虚部展开： 假设1至m为PQ节点，m+1至n-1为PV节点，n为平衡节点。 每一个PQ节点和每一个PV节点都可以列写有功功率的不平衡方程： 每一个PV节点都可以列写无功功率的不平衡方程： 待求变量为所有节点的电压相角： 和PQ节点的电压幅值： 共有n-1+m个待求变量，共有n-1+m个方程。 写成分块矩阵形式的修正方程为 其中， 雅克比矩阵的子块为： 雅克比矩阵的对角线元素， 雅克比矩阵的非对角线元素， 极坐标下的潮流计算的流程 （与直角坐标下的潮流计算的流程类似，在此省略。） PQ分解法潮流计算 前提是极坐标下的修正方程 简化理由：在高压交流输电系统中线路电抗要远大于电阻。 因此：有功功率的变化主要取决于节点电压相位的变化， 无功功率的变化主要取决于节点电压幅值的变化。 在实际计算中，雅克比矩阵的非对角线元素K和N可以忽略不计，这样上式的n-1+m个方程就变成了一个n-1阶方程和一个m阶方程。 即：有功功率