第四章 电力系统潮流的计算机算法

第一节电力网络的数学模型第二节等值变压器模型及其应用第三节节点导纳矩阵的形成和修改第四节功率方程和变量及节点分类第五节高斯－塞德尔法潮流计算第六节牛顿－拉夫逊法潮流计算第七节P－Q分解法潮流计算第四章电力系统潮流的计算机算法机算法的步骤：1、建模2、确定解算方法3、制定计算流程4、编程电力网络的数学模型第一节电力网络的数学模型电力网络的数学模型是将网络有关参数和变量及其相互关系归纳起来，组成可以反映网络性能的数学方程式组，即一种数学描述。1、节点电压方程2、回路电流方程3、割集电压方程1、节点导纳矩阵表示2、节点阻抗矩阵表示1、节点导纳矩阵表示的节点电压方程BBBUYI的节点数为网络中除参考节点外阶节点导纳矩阵：节点电压的列向量量：节点注入电流的列向nn:YUBBnIBnnnnnnnnUUUYYYYYYYYYIII.............................2121222211121121元素－互导纳节点导纳矩阵的非对角素－自导纳节点导纳矩阵的对角元:Y:Yijii互导纳自导纳:Y:Yijii.ii1iYyYiijijii导纳之和相连支路的也等于与节点向网络中注入的电流，经节点，，其他节点全部接地时＝施加单位电压等于在节点＝iiiiiUIUI.jij1iYyyYiijiijjij支路元件导纳的负值之间所连、也等于与节点向网络中注入的电流，经节点，，其他节点全部接地时＝施加单位电压等于在节点＝－＝－iijUIUI结论：1、节点导纳矩阵为对称矩阵2、节点I、j之间没有支路相连时，3、节点导纳矩阵是稀疏矩阵0yyjiij＝例1求节点导纳矩阵75.425.225.755.255.7jjjjjjjjjYB网络中的参数均以电抗标么制给定，试求电力网络的节点导纳矩阵。123j4j2j0.5j0.4j0.22、节点阻抗矩阵表示的节点电压方程BBBUYI的节点数为网络中除参考节点外阶节点阻抗矩阵：节点电压的列向量量：节点注入电流的列向nn:ZUBBnIBnnnnnnnnUUUIIIZZZZZZZZZ.............................2121212222111211阻抗矩阵的互阻抗阻抗矩阵的自阻抗:Z:ZijiiBBBUIZ互阻抗自阻抗:Z:Zijii上的电压值节点，其他节点全部开路时＝注入单位电流等于在节点i1IiZUIUZiiiiiiii的电压。节点，，其他节点全部开路时＝注入单位电流等于在节点＝j1IiZZZijjijijiijUIU结论：1、节点阻抗矩阵为对称矩阵2、所有节点电压都不为零，互阻抗皆为非零元素3、节点导纳矩阵是满矩阵，没有非零元素返回第二节等值变压器模型及其应用不论采用有名制或标么制，凡涉及多电压级网络的计算，都必须将网络中所有参数和变量归算至同一电压级。实际上在计算中有些变压器的实际变比不等于变压器两侧所选电压基准值之比，这些等值电路模型并不能体现变压器实际具有的电压变换功能。以下将介绍另一种可等值地体现变压器电压变换功能的模型，这种模型可体现电压变换，在多电压级网络计算中采用这种变压器模型后，就可不必进行参数和变量的归算。一.变压器为非标准变比时的修正步骤一：从一个未作电压级归算的简单网络入手。令变压器的导纳或励磁支路和线路的导纳支路都可略去，设变压器两侧线路的阻抗都未经归算，即分别为高、低压侧或I，II侧线路的实际阻抗．变压器本身的阻抗则归在低压侧；设变压器的变比为k，其值为高、低压绕组电压之比。''ZZI、二.等值变压器模型步骤二：如在变压器阻抗ZT的左侧串联一变比为K的理想变压器如图示，就如同将变压器及其低压侧线路的阻抗都归算至高压侧，或将高压侧线路的阻抗归算至低压侧，从而实际获得将所有参数和变量都归算在同一侧的等值网络。Z1U1k:1ZTU2Z2I1U1/kI21I2IkII/21步骤三：计算等值模型kSS212111IUIU＝为流出理想变压器的功率＝为流入理想变压器的功率～～～～＝21SSk2111IUIUTZIUkU221/22021021121111kkYykkYkZkykYkZyyTTTTT根据双端口原理：TTTTZUkZUIkZUkZUI2122211Z1U1k:1ZTU2Z2I1U1/kI21I2I221201212212112101)()(UyyUyIUyUyyIZ1YT/kZ2YT(1-k)/k2YT(k-1)/kZ1kZTZ2k2ZT/(1-k)kZT/(k-1)等值三绕组变压器模型p77三.等值变压器模型的应用1、采用有名值，线路参数都未经归算，变压器参数则归算在低压侧。线路阻抗为Z`I、Z`II23210NkTSUPRNkTSUUX100(%)2UI、UII:分别为与变压器高、低压绕组实际匝数相对应的电压。实际变比（理想变压器的变比）k=UI/UII2、采用有名值，线路参数都已按选定变比UIN/UIIN归算至高压侧。23210NkTSUPRNkTSUUX100(%)2UIN、UIIN:分别为归算参数时任选的高低压电压。理想变压器的变比线路阻抗为ZI＝Z`I；ZII＝Z`II2NINUU2NINUU2NINUU)()(*UUUUKININ非标准变比k*3、采用标么值，线路和变压器参数都已按选定的基准电压UIB、UIIB折算为标么值。232*10NkTSUPRNkTSUUX100(%)2*UIB、UIIB:分别为折算参数时任选的变压器高低压侧的基准电压。理想变压器的变比线路阻抗标么值为)()(*UUUUKIBIB2'*IBBIIUSZZ2'*BBUSZZ2BBUS2BBUS返回第三节节点导纳矩阵的形成和修改一.节点导纳矩阵的形成nnnnnnYYYYYYYYY.....................212222111211(1)节点导纳矩阵的阶数等于除参考结点以外的节点数。(2)节点导纳矩阵是稀疏矩阵，非对角非零元素的个数等于对应节点所连的不接地支路数。(3)对角元素（自导纳）等于相应节点所连支路的导纳之和。(4)非对角元素（互导纳）等于两节点间支路导纳的负值。(5)节点导纳矩阵是对称方阵，只需求上三角或是下三角元素。BBBUYI标准变比：在采用有名值时，是指归算参数时所取的变比。采用标么值时，是指折算参数时所取各基准电压之比。非标准变比：统指与这两种情况相对应的理想变压器之比K(6)对于网络中的变压器，采用以导纳表示的等值电路。实际上，往往考虑接入非标准变比的变压器支路I、j时，对原来的节点导纳矩阵的修正。①增加非零非对角元素（ij之间互导纳）②节点I的自导纳，增加一个改变量为③节点j的自导纳，增加一个改变量为*jiijYYKYTTTTiiYkYkkYY***11Z1YT/kZ2YT(k-1)/kYT(1-k)/k2ij*2**11kYkkYYTTjj二.节点导纳矩阵的修改(2)在原有节点ij之间增加一条支路(1)从原来的网络中引出一条新的支路，同时增加一个新的节点。ijZijijijjiijijijjjijijiiz1yYYz1yYz1yY＝－－＝导纳矩阵增加一阶ijZijijijjiijijijjjijijiiz1yYYz1yYz1yY＝－－＝(3)在原有节点ij之间切除一条阻抗为Zij的支路ij－Zijijijjiijijijjjiiz1yYYz1yYY＝＝(4)原有节点ij之间阻抗由Zij变为Zij’ij－ZijZ’ijijijijijjiijijijijijjjiiz'1z1y'yYYz1z'1yy'YY＝＝(4)原有节点ij之间变压器的变比由K*变为K*’时。ij-ZTZTK’*:1K*:1T2*2*T2*2*jij0jjT**T**T**i0iji0iiT**T**jiijZ1)K1K'1(Y)K1K'1(yyYY)K1K'1(YK1KYK'1K'y0yyYZ1)K1K'1(Y)K1K'1(YY＝返回Z1YT/kZ2YT(k-1)/kYT(1-k)/k2第四节功率方程和变量及节点分类一.功率方程（自学）BBBIUYBS*BBBBUSUY～～121U2U~1LS~2LS~1GS~2GS的等值电源功率、分别是母线、2122~211~1GGGGGGjQPSjQPS的等值负荷功率、分别是母线、2122~211~1LLLLLLjQPSjQPS的注入功率、分别是母线、2121~~1~2~~1~1LLSSSSSSGG的注入电流、分别是母线、21222111LGLGIIIIII～～121U2U~1LS~2LS~1GS~2GSy12y20y10～～121U2U~1S~2Sy12y20y101I2I2*212122221*12121111USUUIUSUUIYYYY*1*212*2*2222*2*121*1*1111UUUUSUUUUSYYYY～～])cos[(cos])cos[(cos])sin[(sin])sin[(sin122122222212121111122121222212121111msLGmsLGmsLGmsLGsUUyUyQQQsUUyUyQQQUUyUyPPPUUyUyPPPmsmsmsms是电力系统网络参数、、、其中msmsyy)90(122112)90(212012102211YYYYmsjmjseyyeyyyyy令212211UUjjeUeU二.变量的分类])cos[(cos])cos[(cos])sin[(sin])sin[(sin122122222212121111122121222212121111msLGmsLGmsLGmsLGsUUyUyQQQsUUyUyQQQUUyUyPPPUUyUyPPPmsmsmsms12个变量：负荷的有功和无功——电源发出的有功和无功——母线电压和相位角——2211LLLLQPQP、、、2211GGGGQPQP、、、2121、、、UU不可控变量（扰动变量）d可控变量uU受QG控制，受PG控制－状态变量x对于有n个节点的电力系统（参考节点除外）：1、扰动变量d：2n个2、控制变量u：2n个共6n个3、状态变量x：2n个2n个扰动变量是已知的，给定2（n－1）个控制变量，给定2个状态变量，要求确定2（n－1）个状态变量。已知：4n个变量，待求：2n个变量约束条件：0Qmax.min.max.min.＝＝没有电源的节点：GiGiGiGiGiGiGiGiPQQQPPP控制变量状态变量maxmax.min.jijiiiiUUU三.节点的分类分类给定量待求节点举例PQ节点给定有功、无功的发电厂母线、无电源的变电所母线PV节点有一定无功储备的发电厂、有一定无功电源的变电所母线平衡节点担任调整系统频率的发电厂母线GiLiPP、GiGiLiLiQPQP、、、iiU、LSLsQP、iLiUQ、iiQ、ssU、GsGsQP、其中平衡