电力系统潮流的计算机算法优秀PPT.ppt

《电力系统潮流的计算机算法优秀PPT.ppt》由会员分享，可在线阅读，更多相关《电力系统潮流的计算机算法优秀PPT.ppt（59页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、电力系统潮流的计算机算法你现在浏览的是第一页，共59页第四章第四章 复杂电力系统潮流的计算机算法复杂电力系统潮流的计算机算法u第二节第二节 潮流计算的节点功率方程和节点分类潮流计算的节点功率方程和节点分类u第三节第三节 潮流计算的牛顿潮流计算的牛顿-拉夫逊法拉夫逊法u第一节第一节 网络方程式网络方程式u第四节第四节 牛顿牛顿-拉夫逊潮流计算中的收敛性和稀疏技术拉夫逊潮流计算中的收敛性和稀疏技术u第五节第五节 其他潮流计算方法简介其他潮流计算方法简介你现在浏览的是第二页，共59页第一节第一节第一节第一节 网络方程式网络方程式网络方程式网络方程式网络方程2用节点导纳矩阵表示的网络方程式3用节点阻抗

2、矩阵表示的网络方程式1基本概念你现在浏览的是第三页，共59页网络方程：反映系统中电流与电压之间相互网络方程：反映系统中电流与电压之间相互关系的数学方程；关系的数学方程；如节点电压方程、回路电流方程、割集电压如节点电压方程、回路电流方程、割集电压方程等，一般来讲，由于系统的等值电路方程等，一般来讲，由于系统的等值电路中的接地支路较多，采用节点电压方程时中的接地支路较多，采用节点电压方程时的方程数比回路电流少，故在电力系统潮的方程数比回路电流少，故在电力系统潮流计算中大都采用节点电压方程。流计算中大都采用节点电压方程。基本概念基本概念返回返回你现在浏览的是第四页，共59页1 1、网络方程的形成、网

3、络方程的形成用节点导纳矩阵表示用节点导纳矩阵表示的网络方程式的网络方程式一般规定：外部向系统注入的功率为节点功率的一般规定：外部向系统注入的功率为节点功率的正方向，且在等值电路中，与节点注入功率相对正方向，且在等值电路中，与节点注入功率相对应的电流称为节点注入电流，它的规定正方向与应的电流称为节点注入电流，它的规定正方向与注入功率一致。注入功率一致。以书本图以书本图4-1为例形成节点导纳矩阵表示的网为例形成节点导纳矩阵表示的网络方程式。络方程式。你现在浏览的是第五页，共59页推导出一般情况简写简写U UY Y=.I.你现在浏览的是第六页，共59页自自导导纳纳互互导导纳纳节点节点i:加单位电压加

4、单位电压其余节点其余节点j:全部接地全部接地节点节点 i 注入网络电流注入网络电流Yii0节点节点i:加单位电压加单位电压其余节点其余节点j:全部接地全部接地由地流向节点由地流向节点j的电流的电流稀疏性：当稀疏性：当yij=0 时时Yij=0节点导纳距阵的特点：节点导纳距阵的特点：1、阶数、阶数2、对称性、对称性3、稀疏性、稀疏性2 2、节点导纳矩阵的物理意义和特点、节点导纳矩阵的物理意义和特点Yii节点节点i自导纳，等于与自导纳，等于与i相连所有支路导纳之和；相连所有支路导纳之和；Yij节点节点i，j间的自导纳，等于节点间的自导纳，等于节点i，j间支路导纳的负值间支路导纳的负值你现在浏览的是

5、第七页，共59页补充知识：导纳矩阵的修改补充知识：导纳矩阵的修改补充知识：导纳矩阵的修改补充知识：导纳矩阵的修改1)增加一节点增加一节点ijyij2)增加一条支路增加一条支路ijyij节点导纳矩阵增加一阶节点导纳矩阵增加一阶Yii=yij Yjj=yijYij =Yji =-yij 导纳矩阵的阶数不变导纳矩阵的阶数不变 Yii=Yjj=yij Yij=Yji=-yij 你现在浏览的是第八页，共59页3)切除一条支路切除一条支路4)修改一条支路的导纳值（修改一条支路的导纳值（yij 改变为改变为yij）ij-yij导纳矩阵的阶数不变导纳矩阵的阶数不变 Yii=Yjj=yij -yij Yij=Y

6、ji=yij -yij 导纳矩阵的阶数不变导纳矩阵的阶数不变相当于增加一导纳为（相当于增加一导纳为（-yij）的支路）的支路 Yii=Yjj=-yij Yij=Yji=yij ij-yijyij 你现在浏览的是第九页，共59页5)修改一条支路的变压器变比值（修改一条支路的变压器变比值（k*改变为改变为k*）ijyT/k*yT(k*-1)/k*yT(1-k*)/k*2 Yii=0 Yij=Yji=-(1/k*-1/k*)yTYjj=(1/k*2-1/k*2)yT返回返回你现在浏览的是第十页，共59页用节点阻抗矩阵表示用节点阻抗矩阵表示的网络方程式的网络方程式1 1、阻抗矩阵形式网络方程的形成、阻

7、抗矩阵形式网络方程的形成简写简写I IZ Z=.U.你现在浏览的是第十一页，共59页自自阻阻抗抗节点节点i:加单位电流加单位电流其余节点其余节点j:全部与外电路断开全部与外电路断开节点节点 i 的电压的电压节点阻抗距阵节点阻抗距阵的特点：的特点：1、N阶数、对阶数、对称性称性2、满阵、满阵3、不容易求得、不容易求得2 2、节点阻抗矩阵的物理意义和特点、节点阻抗矩阵的物理意义和特点互互阻阻抗抗节点节点i:加单位电流加单位电流其余节点其余节点j:全部与外电路断开全部与外电路断开节点节点 j 的电压的电压返回返回你现在浏览的是第十二页，共59页第二节第二节第二节第二节 潮流计算的节点功率方程和节潮流

8、计算的节点功率方程和节潮流计算的节点功率方程和节潮流计算的节点功率方程和节点分类点分类点分类点分类极坐标表示极坐标表示的节点功率方程的节点功率方程直角标表示直角标表示的节点功率方程的节点功率方程节点的分类节点的分类主要知识主要知识你现在浏览的是第十三页，共59页GG12y10y20y12极坐标表示极坐标表示的节点功率方程的节点功率方程简单系统的等值网络简单系统的等值网络你现在浏览的是第十四页，共59页第一步：第一步：第二步：第二步：第三步：第三步：相相位位差差决决定定潮潮流流分分布布返回返回你现在浏览的是第十五页，共59页直角标表示直角标表示的节点功率方程的节点功率方程节点电压相量用实部和虚部

9、表示节点电压相量用实部和虚部表示返回返回你现在浏览的是第十六页，共59页节点的分类节点的分类一般节点负荷节点：该节点上没有发电机而只有负荷；发电机节点：该节点上只有发电机而没有负荷；联络节点：该节点上既没有发电机，也没有负荷；你现在浏览的是第十七页，共59页约束条件约束条件实际电力系统运行要求：实际电力系统运行要求：电能质量约束条件：电能质量约束条件：电能质量约束条件：电能质量约束条件：U U U Uimin imin U U U Ui i U Uimaximax电压相角约束条件电压相角约束条件|ijij|=|=|=|=|i i-j j|ijmaxijmax,稳定运行的一个重要条件。稳定运行的

10、一个重要条件。有功、无功约束条件有功、无功约束条件 P Pimin imin P Pi i i i P Pimax imax Q Q Q Qimin imin Q Qi i Q Qimaximaximaximax你现在浏览的是第十八页，共59页潮流计算中的节点分类电力系统有电力系统有电力系统有电力系统有n n个节点，每个节点可能有个节点，每个节点可能有个节点，每个节点可能有个节点，每个节点可能有4 4个变量个变量个变量个变量，则共有则共有则共有则共有4n4n个个个个变量，而上述功率方程只有变量，而上述功率方程只有变量，而上述功率方程只有变量，而上述功率方程只有2n2n个，所以需要事先给定个，所

11、以需要事先给定个，所以需要事先给定个，所以需要事先给定2n2n个变量个变量个变量个变量的值。根据各个节点的已知量的不同，分成三类：的值。根据各个节点的已知量的不同，分成三类：的值。根据各个节点的已知量的不同，分成三类：的值。根据各个节点的已知量的不同，分成三类：PQPQ节点、节点、节点、节点、PV PV 节点、平衡节点。节点、平衡节点。节点、平衡节点。节点、平衡节点。PQPQ节点：给定的是注入有功功率节点：给定的是注入有功功率节点：给定的是注入有功功率节点：给定的是注入有功功率P P和注入无功功率和注入无功功率和注入无功功率和注入无功功率QQ；待求；待求；待求；待求量是节点电压有效值量是节点电

12、压有效值量是节点电压有效值量是节点电压有效值U U和电压的相位和电压的相位和电压的相位和电压的相位 。PVPV节点：节点：节点：节点：给定的是注入有功功率给定的是注入有功功率给定的是注入有功功率给定的是注入有功功率P P和节点电压有效值和节点电压有效值和节点电压有效值和节点电压有效值U U；待求；待求；待求；待求量是节点的注入无功功率和电压的相位量是节点的注入无功功率和电压的相位量是节点的注入无功功率和电压的相位量是节点的注入无功功率和电压的相位 。平衡节点：平衡节点：平衡节点：平衡节点：给定的是节点电压有效值给定的是节点电压有效值给定的是节点电压有效值给定的是节点电压有效值U U和电压的相位

13、和电压的相位和电压的相位和电压的相位 ；待；待；待；待求量是节点的注入无功功率求量是节点的注入无功功率求量是节点的注入无功功率求量是节点的注入无功功率QQ和注入有功功率和注入有功功率和注入有功功率和注入有功功率P P。也称为。也称为。也称为。也称为V-V-节点、松弛节点或电压参考节点节点、松弛节点或电压参考节点节点、松弛节点或电压参考节点节点、松弛节点或电压参考节点你现在浏览的是第十九页，共59页设置平衡节点的目的设置平衡节点的目的在结果未出来之前，网损是未知的，至少在结果未出来之前，网损是未知的，至少需要一个节点的功率不能给定，用来平衡全需要一个节点的功率不能给定，用来平衡全网功率。网功率。

14、电压计算需要参考节点。电压计算需要参考节点。返回返回你现在浏览的是第二十页，共59页第三节第三节第三节第三节 潮流计算的牛顿潮流计算的牛顿潮流计算的牛顿潮流计算的牛顿-拉夫逊法拉夫逊法拉夫逊法拉夫逊法一、N-R的原理和一般方法的原理和一般方法1.1.非线性方程的求解：非线性方程的求解：f(x)=0设：设：x x(0)(0)为的初始近似解，为的初始近似解，x x(0)(0)为与真实解的偏差为与真实解的偏差则：则：x=x(0)x(0)f(x(0)x(0)=0 按按TaylorTaylors s展开展开 f(x(0)x(0)=f(x(0)-f(x(0)x(0)+.+(-1)n fn(x(0)(x(0

15、)n/n!+.=0你现在浏览的是第二十一页，共59页 由于由于x x(0(0）较小，故忽略高次项后：较小，故忽略高次项后：f(xf(x(0)(0)x x(0)(0)=f(x)=f(x(0)(0)-f(x)-f(x(0)(0)x x(0)(0)=0=0 x x(0)(0)=f(x=f(x(0)(0)/f(x)/f(x(0)(0)x x(1)(1)=x=x(0)(0)-x x(0)(0)=x=x(0)(0)-f(x-f(x(0)(0)/f(x)/f(x(0)(0)k k次迭代时修正方程为：次迭代时修正方程为：f(x f(x(k)(k)-f(x)-f(x(k)(k)x x(k)(k)=0 =0 x

16、x(k)(k)=f(x=f(x(k)(k)/f(x)/f(x(k)(k)x x(k+1)(k+1)=x=x(k)(k)-f(x-f(x(k)(k)/f(x)/f(x(k)(k)你现在浏览的是第二十二页，共59页 结束迭代的条件（收敛）结束迭代的条件（收敛）：|f(x|f(x(k)(k)|)|1 1 或或|x x(k)(k)|RRijij，B BijijGGijij，ijij00。P P，Q UQ U），得出的一种简化形式。），得出的一种简化形式。图图形形解解释释你现在浏览的是第四十六页，共59页第五节第五节 P-Q P-Q分解法潮流计算分解法潮流计算二、二、P-QP-Q分解法的修正方程式分解法

17、的修正方程式重写极座标方程重写极座标方程(4-53)你现在浏览的是第四十七页，共59页第五节第五节 P-Q P-Q分解法潮流计算分解法潮流计算简写为简写为(4-54)进一步进一步(4-55)计及计及cosij1,G1,Gij ij sinij B Bijij你现在浏览的是第四十八页，共59页第五节第五节 P-Q P-Q分解法潮流计算分解法潮流计算(4-49a)(4-49b)(4-56a)(4-49c)(4-49d)(4-43b)(4-56b)你现在浏览的是第四十九页，共59页第五节第五节 P-Q P-Q分解法潮流计算分解法潮流计算(4-57)你现在浏览的是第五十页，共59页第五节第五节 P-Q

18、 P-Q分解法潮流计算分解法潮流计算(6-75)(4-58a)(4-58b)你现在浏览的是第五十一页，共59页P1/U1P2/U2Pn/UnB11B12B1nB21B22B2nBn1Bn2BnnU11U22Unn(4-59a)Q1/U1Q2/U2Qm/UmB11B12B21B22B2mBm1Bm2Bmm U1 U2 Um(4-59b)B1mP/U=BUQ/U=B U(4-60a)(4-60b)简写为：简写为：你现在浏览的是第五十二页，共59页P-QP-Q分解法的修正方程式的特点：分解法的修正方程式的特点：以一个以一个(n-1)(n-1)阶和一个阶和一个(m-1)(m-1)阶系数矩阵阶系数矩阵B

19、 B、B B 替代原有的替代原有的(n+m-2)(n+m-2)阶系数矩阵阶系数矩阵J J，提高了，提高了计算速度，降低了对存储容量的要求。计算速度，降低了对存储容量的要求。以迭代过程中不变的系数矩阵以迭代过程中不变的系数矩阵以迭代过程中不变的系数矩阵以迭代过程中不变的系数矩阵B B B B、B B 替代变化替代变化的系数矩阵的系数矩阵J J，显著地提高了计算速度。，显著地提高了计算速度。以对称的系数矩阵以对称的系数矩阵B B、B B B B替代不对称的系数替代不对称的系数矩阵矩阵J J，使求逆等运算量和所需的存储容量大为，使求逆等运算量和所需的存储容量大为减少。减少。你现在浏览的是第五十三页，

20、共59页 牛顿拉夫逊法牛顿拉夫逊法和和P PQ Q分解法的分解法的特性：特性：牛顿拉夫逊法牛顿拉夫逊法P PQ Q分解法分解法你现在浏览的是第五十四页，共59页三、三、P-QP-Q分解法的潮流计算的基本步骤分解法的潮流计算的基本步骤1)1)形成系数矩阵形成系数矩阵B B、B B ，并求其逆矩阵。，并求其逆矩阵。2)2)设各节点电压的初值设各节点电压的初值 I I(0)(0)(i=1,2,(i=1,2,n,i,n,i s)s)。U UI I(0)(0)(i=1,2,(i=1,2,m,i,m,i s)s)3)3)按式（按式（4 445a45a）计算有功不平衡量）计算有功不平衡量 P PI I(0)

21、(0)(i=1,2,(i=1,2,n,i,n,i s)s)。4)4)解修正方程式，求各节点电压相位的变量解修正方程式，求各节点电压相位的变量 I I(0)(0)(i=1,2,(i=1,2,n,i,n,i s)s)5)5)求各节点电压相位的新值求各节点电压相位的新值 I I(1)(1)=I I(0)(0)+I I(0)(0)(i=1,2,(i=1,2,n,i,n,i s)s)6)6)按式（按式（4 445a45a）计算无功不平衡量）计算无功不平衡量 Q QI I(0)(0)(i=1,2,(i=1,2,m,i,m,i s)s)。7)7)解修正方程式，求各节点电压幅值的变量解修正方程式，求各节点电压

22、幅值的变量 U UI I(0)(0)(i=1,2,(i=1,2,m,i,m,i s)s)8)8)求各节点电压幅值的新值求各节点电压幅值的新值U UI I(1)(1)=U UI I(0)(0)+U UI I(0)(0)(i=1,2,(i=1,2,m,i,m,i s)s)9)9)不收敛时，运用各节点电压的新值不收敛时，运用各节点电压的新值自第三步开始进入下一次迭代。自第三步开始进入下一次迭代。10)10)计算平衡节点功率和线路功率。计算平衡节点功率和线路功率。你现在浏览的是第五十五页，共59页见书上P175P-QP-Q分解的潮流分解的潮流 计算流程图计算流程图你现在浏览的是第五十六页，共59页作业

23、作业1 1、某系统的等值电路如图所示，试求：、某系统的等值电路如图所示，试求：1)1)写出节点导纳距阵。写出节点导纳距阵。2)2)如果节点如果节点2 2、4 4之间阻抗为之间阻抗为0 0，即，即Z Z2424=0=0，修改导纳距阵。，修改导纳距阵。2 2、潮流计算中，节点的分类及各自的已知和未知量是什么、潮流计算中，节点的分类及各自的已知和未知量是什么?1:1.1j0.1j0.2j0.25j0.13 3、试述、试述牛顿牛顿牛顿牛顿拉夫逊法潮流计算的基本步骤拉夫逊法潮流计算的基本步骤拉夫逊法潮流计算的基本步骤拉夫逊法潮流计算的基本步骤你现在浏览的是第五十七页，共59页参考书参考书 电力系统分析与设计电力系统分析与设计 (Power System Analysis and Design)(Power System Analysis and Design)（美）（美）J.Duncan Glover Mulukutla.S.SarmaJ.Duncan Glover Mulukutla.S.Sarma机械工业出版社机械工业出版社你现在浏览的是第五十八页，共59页ThanksThanks！第四章结束第四章结束你现在浏览的是第五十九页，共59页