4复杂的潮流计算(8学时)教程.ppt

补充说明： 1.负荷需求的功率取决于用户，是无法控制的。 2.电源点发出的有功和无功功率是可以由运行人员控制或改变的。 3.各个节点的电压幅值或相角随着控制变量的改变而变化。 补充说明： 1.一般在潮流计算中只设一个平衡节点s，即它的电压幅值为给定值，相位为零，即该节点的电压方向作为参考方向。平衡节点的待求量是节点注入功率 ，实际上，整个系统的功率平衡是由这一节点完成的 。 2.潮流计算中有大量的这一类节点，即它的节点功率给定，待求量为Ui， θi。只要可以给定功率的发电厂或负荷节点，均可处理为PQ节 。 3.在实际电力系统中PV节点对应于有一定无功功率储备的发电厂电压母线和有无功补偿设备的变电站母线。这类节点的给定量为有功 Pi 和电压幅值 Vi ，待求量为电压相角θi 和节点无功Qi 。 补充说明： 1.很明显可以看出来，不管是直角坐标形式还是极坐标形式，潮流方程都是一组非线性代数方程。求解非线性方程的最基本方法就是迭代法。 2.20世纪50年代，计算机开始应用于电力系统。最早的潮流计算机方法是以节点导纳矩阵为基础的逐次逼近法（高斯－塞德尔迭代导纳法），后来为解决导纳法的收敛性较差的问题，又出现了以阻抗矩阵为基础的逐次逼近法（高斯－塞德尔迭代的阻抗法）； 到20世纪60年代，针对阻抗法占用大量计算机内存的问题又提出了分块阻抗法以及牛顿－拉夫逊法（N－R法）。N－R法是数学上解非线性方程式的有效方法，具有良好的收敛性，并且，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使N－R法在收敛性、占用内存、计算速度方面的优点都超过了阻抗法，因而得到了广泛应用。 1974年，B.Stott和O.Alsac提出了快速分解法（也称为PQ分解法或PQ解耦法），通过对网络模型进行简化，计算速度大大超过了N－R法，不但能应用于离线潮流计算，而且也能应用于在线潮流计算。对于110kV及以上电压等级的系统，架空线的 r x 时，这种简化方法还是比较合理的。 *-207 4节点系统 *-207 4节点系统的数据文件： 4，4，1.05，0.000001 0 1，1，2，0.1，0.4，0.01528 2，1，4，0.12，0.5，0.01920 3，2，4，0.08，0.4，0.01413 0 1，1，3，0，0.3，0.98 0 0 1，1，0，0，0.3，0.18 2，2，0，0，0.55，0.13 3，3，0.5，0，0，0 0 1，3，1.1，0，0 0 *-207 清零：U，I，G，B，P，Q 读入节点数，平衡节点，平衡节点电压，计算精度 读线路参数 GIJ=r/(r×r+x×x)，BIJ=-x/(r×r+x×x) G(j,i)=G(i,j)=-GIJ，B(j,i)=B(i,j)=-BIJ G(i,i)+= GIJ，B(i,i)+= BIJ +b/2 G(j,j)+= GIJ，B(j,j)+= BIJ +b/2 读变压器参数 GIJ=r/(r×r+x×x)，BIJ=-x/(r×r+x×x) G(j,i)=G(i,j)=-GIJ/k0，B(j,i)=B(i,j)=-BIJ/k0，G(i,i)+= GIJ/k0/k0，B(i,i)+= BIJ/k0/k0，G(j,j)+= GIJ，B(j,j)+= BIJ 读接地支路参数 G(i,i)+=gi，B(i,i)+=bi 转换成G+jB=Y∠α 形成节点导纳矩阵流程图 *-207 第一次消元 第k次消元 第n次消元 消元过程 回代过程 高斯消去法流程图 *-207 利用LU（Crout）分解求解修正方程 对于r=2,3,…n计算： 计算U的第r行元素： 计算L的第r列元素： *-207 二. 计算实例 说明： 1）采用中国版BPA潮流程序2.1版 2）采用IEEE22节点系统作为算例 *-207 1、IEEE22节点电网接线 *-207 2、潮流计算条件设置 发电机节点 有功出力P 无功出力Q B2 6.00 3.20 B3 3.10 0.50 B4 1.60 0.70 B5 4.30 3.34 B6 4.00 0.32 负荷节点 有功负荷P 无功负荷Q B8 2.87 1.44 B9 3.76 2.21 B16 5.0 2.9 B18 3.50 2.60 B19 0.86 0.66 B20 0.72 0.47 B21 0.70 0.50 计算方法：牛顿法 初始电压：Vx = 1.0 Vy = 0.0 计算精度：0.0001 电压限值：Vmax＝1.2 Vmin = 0.8 *-207 3、潮流计算过程 节点号