课件第4章计算机潮流20130323_ok.pptVIP

* * §4-5 牛顿-拉夫逊法潮流计算 * * ?N-R法是先将非线性方程组化成线性方程组，然后迭代求解； ?从修正方程式可见，此法是以切线逼近真解的方法，有较高的收敛速度 ；但对初值要求较高； 例如： 由上可见： §4-5 牛顿-拉夫逊法潮流计算 * * x y y=f(x) f(x(0)) X(2) X(0) X(1) X(0) X(1) X(2) N-R法要求初值接近于真解，否则可能不收敛 ?每次迭代J元素总是变化，从而计算量大； ? ΔX(0) §4-5 牛顿-拉夫逊法潮流计算 * * 潮流计算的方程组采用功率方程组，先求出其修正方程式： 两边取共轭： 按照电压的表达形式不同修正方程式有两种： §4-5 牛顿-拉夫逊法潮流计算 1. 节点电压以直角表示的修正方程式： * * 电压表示为： 导纳表示为： 带入功率方程： 功率方程： §4-5 牛顿-拉夫逊法潮流计算 ? 功率方程： * * 将实部、虚部分开： 每个节点列2个方程，n个节点的系统将有2n个方程； 假设：一个n个节点的系统， 所以，各节点所列的方程数和变量如下表： §4-5 牛顿-拉夫逊法潮流计算 * * 节点类型 已知条件 可列方程 未知量 总方程数 总变量数 PQ P、Q ΔPi=0 ΔQi=0 ei、fi 2m 2m PV P、V ΔPi=0 ei、fi 2(n-1-m) 2(n-1-m) Vδ ei、fi 不列方程 P、Q 总方程数和变量数相同：2m+2(n-1-m)=2(n-1)个； §4-5 牛顿-拉夫逊法潮流计算 * * 由于迭代的变量是ei、fi，平衡节点的e1、f1已知，不参加迭代；方程的数目确定后，可得N-R法的修正方程式： §4-5 牛顿-拉夫逊法潮流计算 ? N-R法的修正方程式： * * 功率电压不平衡量 雅克比矩阵 电压修正量 2m个PQ节点 2(n-m)-2个 PV节点 §4-5 牛顿-拉夫逊法潮流计算 * * ?不平衡量的计算：将迭代值直接带入功率方程和电压方程； ?雅克比矩阵元素计算：功率方程和电压方程求导后带入迭代值计算； 例如： 书上已经推出了全部的雅克比矩阵的元素公式，应该记住； 有元素计算公式 元素公式与电压有关，每次迭代元素就变化。 §4-5 牛顿-拉夫逊法潮流计算 说明： * * 可见：每迭代一次，雅克比矩阵元素就要重新计算一次，导致计算量大； 求得潮流计算的修正方程式后，可给出N-R法的计算步骤： ? 建立导纳矩阵[Y]; ? 设定电压初值： ? 带入电压、功率方程求出不平衡量： §4-5 牛顿-拉夫逊法潮流计算 ? 牛顿拉夫逊法的计算步骤： * * 注意：由于带入的初值并不是真值，算出的结果不为零 ? 带入雅克比矩阵[J]的元素公式，求雅克比矩阵[J]： ? 用修正方程式求出电压修正量： 得到第一次迭代值： ? 判断收敛： §4-5 牛顿-拉夫逊法潮流计算 * * ? 上式成立，迭代结束，并求出各支路功率，若不成立时： 注意：无功功率的越限处理： 对于PV节点，每次迭代要计算QP值，并用约束条件比较，当越限时，该点就转化为PQ节点了（这时，修正方程式的结构和元素公式要改变）。 (直角坐标) 1 §4-5 牛顿-拉夫逊法潮流计算 ? N-R法的计算机计算框图： * * 输入原始数据 形成[Y] 给定电压初始值 K=0 计算各节点不平衡量： 计算雅克比矩阵各元素 解修正方程式，得： 修正变量： K=K+1 计算平衡节点功率和电网的线路功率 输出结果 Y N 结束