高等电力网络分析课件.pptVIP

高等电力网络分析

——牛顿-拉夫逊潮流计算仿真;电力系统潮流计算的意义和要求;给定电力系统的网络结构、参数和决定系统运行状况的边界条件，电力系统的稳态运行状态便随之确定。潮流计算就是要通过数值仿真的方法把电力系统的详细运行状态呈献给运行和规划人员，以便研究系统在给定条件下的稳态运行特点。;最初在计算机上实现的潮流计算方法是以导纳矩阵为基础的高斯迭代法。这种方法内存需求小，但收敛性差。后来发展了以阻抗矩阵为基础的算法。这种方法收敛性好，但内存占用量大大增加，限制了解题规模。牛顿-拉夫逊方法是解非线性代数方程组的一种基本方法，在潮流计算中也得到了应用。20世纪60年代中后期，稀疏矩阵技术和节点编号优化技术的提出使牛拉法的解题规模和计算效率进一步提高，至今仍是潮流计算中的广泛采用的优秀算法。;由于潮流计算在电力系统分析中所处的特殊地位和作用，对潮流计算的要求可以归纳为如下几点：

（1）算法的可靠性或收敛性

（2）计算速度和内存占用量

（3）计算的方便性和灵活性

潮流计算的结果意义重大，是电力系统运行、规划以及安全性、可靠性分析和优化的基础，也是各种电磁暂态和机电暂态分析的基础和出发点。;对于N个节点的电力网络（地作为参考节点不包括在内），如果网络结构和元件参数已知，则网络方程可以表示为：;电力系统的计算中，给定的运行变量是节点注入功率，不是节点注入电流，那么两者之间有如下关系：;上式就是潮流方程的复数形式，是N维的非线性复数代数方程组。将其展开，有：;式中;如果节点电压用极坐标表示，即令，代入式（1-3）中则有：;牛顿-拉夫逊法的求解步骤如下。在给定的初值

处作一阶泰勒展开：;用就得到x的新值。如果用k表示迭代次数，写成一般的表达式，有：;上式也可以写成下面的简单迭代法的计算格式;牛顿-拉夫逊法的极坐标形式;牛顿-拉夫逊法的极坐标形式;程序流程图;9;支路数据（标幺值，SB=100MVA）;变压器数据（标幺值）;IEEE14母线标准试验;IEEE14母线标准试验;IEEE14母线标准试验;IEEE14母线标???试验;IEEE14母线标准试验;IEEE14母线标准试验;求导纳矩阵：

form=1:14

t(m)=0;

forn=1:14

t(m)=t(m)+y(m,n);

end

ifm==4

t(m)=t(m)+0.190;

end

end

form=1:14

d(m)=0;

forn=1:14

ifm==n

;程序清单;计算雅克比矩阵各元素:

form=1:N1

forn=1:N1+1

h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));

n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));

j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));

L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));

end

ifm=9

H(m,m)=sum(h0)-u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));

N(m,m)=sum(n0)-2*u(m)^2*G(m,m)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));

J(m,m)=sum(j0)+u(m)*u(m)*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));

L(m,m)=sum(L0)+2*u(m)^2*B(m,m)+u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));

else

H(m,m)=sum(h0)-u(m)*u(m)*(G(m,m)*sin(delt(m)-delt(m))-B(m,m