潮流计算问题

1、潮流计算的定义（课后题）各种潮流计算模型和算法的特点、适用范围以及相互之间的区别和联系（课后题）影响潮流收敛性的因素，以及如何改善潮流计算的收敛性（课后题）通过功率方程说明为什么潮流计算的数学模型是非线性的应该采用什么样的数学方法 求解（03A、05A）电力系统的潮流计算有哪些常规算法有哪些扩展算法（05B）潮流计算的目的是什么其数学模型是什么有何特点（06B）简要说明潮流计算的概念、模型及计算方法。（07B）高斯赛德尔迭代法和牛顿拉夫逊迭代法是常规的潮流计算方法，请介绍一下最优潮流（OPF）算法的原理及其应用。（04电科院）潮流计算的目的：常规潮流计算的目的是在已知电力网络参数和各节点的注入

2、量的条件下,求解各节点电 压。目的1：在电网规划阶段，通过潮流计算，合理规划电源容量和接入点，合理规划网架，选择无 功补偿方案，满足规划水平年的大小方式下潮流交换控制、调峰、调相、调压的要求。在编制年运行方式，在预计复合增长及新设备投运基础上，选择典型方式进行潮流计算， 发现电网中的薄弱环节，供调度人员异常调度控制参考，并对规划、基建部门提出改进 网架结构，加快基建进度的建议。正常检修及特殊运行方式下的潮流计算，用于日常运行方式的编制，指导发电厂开机方 式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求及电压质量要 求。预想事故、设备退出运行对静态安全分析的影响及做出预想的运行方

3、式调整方案。目的2：检查电力系统各元件是否过负荷；检查电力系统各节点的电压是否满足电压质量的要求；根据对各种运行方式的潮流分布计算，可以正确的选择系统接线方式，合理调整负荷，以保证电力系统安全、可靠地运行，向用户供给高质量的电能;根据功率分布，可以选择电力系统的电气设备和导线截面积，可以为电力系统继电保护整定计算提供必要的数据等；为电力系统扩建和规划提供依据；为调压计算、经济运行计算、短路计算等提供必要的数据。数学模型：数学模型为：潮流计算所用的电力网络系由变压器、输电线路、电容器、电 抗器等静止线性元件所构成，并用集中参数表示的串联或并联等值支路来模拟。普遍采用节 点法，I = YU来建立潮

4、流计算的数学模型。在实际工程中，节点注入量不是电流，而是节点功率，因此节点电压方程要进行修改：I = Pi - jQi ,(i = 1,2,.,n)，进一步得到*UiPi - jQi =祝Y U* ,(i = 1,2,., n)，上式为电压的非线性隐函数，无法直接求解，必须通Uj=i j ji过一定的算法求近似解。这是潮流计算问题最基本的方程式，是一个以节点电压U为变量 的非线性代数方程组，采用节点功率作为节点注入量是造成方程组呈非线性的根本原因。对 于每个节点，要确定其运行状态，需要四个变量，P、Q、U、0，个节点共4n个运行变 量需要确定，如果将节点电压方程式的实部和虚部拆开，形成2n个实

5、数方程，在潮流计算 前，必须先确定2n个变量作为已知量。这样潮流方程就可解。根据节点电压表示方式的不 同(Ui = ei + jfi，U. = Ui)，可以得到直角坐标系和极坐标系下的潮流方程。直角坐标系下功率方程：P - e (G e - B f )- f (G f + B e )= 0(i = 1,2,，n)i i ij j ij j i ij j ij j j=1j=1Q - f (G e - B f )+ e (G f + B e )= 0i i ij j ij j i ij j ij jj=1j=1极坐标系下功率方程：UG cos(5 -5 ) + B sin(5 -5 ) = 0j

6、 ij i j ij i jP - U (i = 1,2,，n)iiQi + Ui j=1j=1UB cos(5 -5 )-G sin(5 -5 ) = 0 j L ij i j ij i j 常规算法有：高斯-赛德尔法、牛顿-拉夫逊法、快速解耦法扩展算法有：保留非线性潮流算法、最小化潮流算法、最优潮流、直流潮流法、随机潮流法、 三相潮流高斯-赛德尔法：高斯-赛德尔法的迭代格式为：U(5= -!P - jQi-ILy U 俄)Ki = 2,3,., n)i Y* (k)ij jii Uj=iij色收敛判据为：max J WUJ优点：原理简单，程序设计容易。导纳矩阵是一个对称且高度稀疏的矩阵，因

7、此占用内存非 常节省。缺点：收敛速度很慢，算法收敛所需的迭代次数与所计算网络的节点数目有密切关系，在系 统病态的情况下，收敛困难。重负荷节点；负电抗支路；较长辐射型线路；长短线路接在同一节点上；且长短线路的比值很大；牛顿-拉夫逊法：该算法实际上是非线性方程或非线性方程组的多次线性逼近。| f( X (k )Ax (k) = - f (x( k)牛顿法的迭代格式为：fx (k+1) = x (k) + Ax (k)修正方程有极坐标形式和直角坐标形式：-AP -=-H“1PA0_和AP -AQ HMN_LAeAQML 1 AU / URSAf一 一|_ _AU 2_修正方程的特点：在PV节点所占比

8、例不大时，两者方程的数目都接近2(n-1)。雅可比矩阵的元素都是节点电压的函数，每次迭代，雅可比矩阵都需要重新形成。3）按节点序号顺序而构成的分块雅可比矩阵将和节点导纳矩阵具有同样的稀疏结构， 是一个高度稀疏的矩阵。4）雅可比矩阵不是对称矩阵。牛顿法的核心就是反复形成并求解修正方程式。提高牛顿法性能：采用稀疏技术，排零存储，排零运算。求解过程边形成、边消元、边存储。节点编号优化，采用半动态法。（静态法：按节点静态连接支路的多少顺序编号；半动态法： 按节点动态连接支路的多少顺序编号；动态法：按节点动态增加支路的多少顺序编号;） 牛顿法的性能和特点：1）平方收敛，开始时收敛比较慢，在几次迭代后，收

9、敛得非常快，其迭代次数和系统 的规模关系不大，如果程序设计良好，每次迭代的计算量仅与节点数成正比。2）对初值很敏感，有时需要其他算法为其提供初值。如果初值选择不当，可能根本不 收敛或收敛到一个无法运行的解点上。3）对函数的平滑性敏感，所处理的函数越接近线性，收敛性越好，为改善功率方程的 非线性，实用中可以通过限制修正量的幅度来达到目的。但幅度不能太小。4）对以节点导纳矩阵为基础的G-S法呈病态的系统，N-L法一般都能可靠收敛。快速解耦法：N-L法的J阵在每次迭代的过程中都要发生变化，需要重新形成和求解，这占 据了 N-L法的大部分计算时间，这也是N-L法速度不能提高的原因。N-L法可以简化成为

10、定 雅可比矩阵法，如果固定的迭代矩阵构造得当，定雅可比矩阵法可以收敛，但只有线性收敛 速度。APH N 蜀一AQM LJ/U 弟一步假设:由于R X 有功无功解耦国希JaU7U nw = L .AU / U 第二步假设：一般线路两端电压相角差 七较小（一般1020度），且G|j Bq 有： cos5 * 1 G sin6 B cos6 得到：H =UUB , L =UUB ijij ij ij ijij i J ij ij i J ijNP H I 0 区AP = H M第三步假设：H =-Q -U2B , L = Q -U2B为正常情况下节点i的注入无功功 iii i ii ii i i i

11、i率；此时其他节点未接地：U2Bii为除i节点外其他节点接地时，由节点i注入的无功功率;所以Q. U2B，得： H =-U2B，L =-U2B。修正方程缩写为:P/U = -B UA5Q/U = -B”U继续简化：1）形成B时略去那些主要影响无功功率和电压幅值，而对有功功率及电压相角关系 很少的因素。这包括输电线路的充电电容以及变压器非标准变比。2）为了减少迭代过程中无功功率及电压幅值对有功迭代的影响，将式上式1右端的电 压均置为标幺值。3）形成B时，略去串联元件的电阻。AP/U=-BAS最终表达式为：aq/u=-bau算法特点：（等斜率法，从平方收敛退化为线性收敛，所以迭代次数比牛顿法多）1

12、）用两个阶数几乎减半的方程组代替原方程组，显著减少了内存量和计算量；2）迭代矩阵为常数阵，只需形成求解一次，大大缩短每次迭代所需时间；3）迭代矩阵对称，可上（下）三角存储，减少内存量和计算量；基于以上原因，该算法内存需要量为N-L法的60%，每次迭代所需时间为N-L法的1/5。4）线性收敛，收敛次数多于N-L法，但总的计算速度任能大幅度提高。5）对R/X过大的病态条件以及线路特别重载的情况下，可能不收敛，一般适用于110kV 及以上的电网。6）由于算法的精确程度取决于，P-Q分解法的近似处理只影响计算过程，并不影响 结果的精度。关于元件大R/X比值病态问题，采用补偿法或者对算法加以改进：在于对

13、B和B”元件 电阻的取舍问题：若在B中不计串联元件电阻，仅用其电抗值X，而在B”中仍用精确的电纳值B,或者在8中忽略串联元件电阻而用精确的电纳值B,而在B中采用元件电抗值X, 分别称其为XB方案和BX方案。保留非线性潮流计算法：潮流问题其实是求解一个不含变量一次项的二次方程组，泰勒 级数只要取三项就能够得到一个没有截断误差的精确展开式。y, (x) = y,(x(o) + I (0)j=j在初值了(o)附近展开，可得到如下没有截断误差的精确展开式：Ax +室工IAx Axx=x(0)j 2!dx dx x=x(0)j kj= k Tj kAx Ax1Ax Ax，迭代格式为:写成矩阵形式:y $

14、 = y(x(0) + JAx + H2Ax AxAx (k) Ax( k)1Ax (k) Ax( k)2，但H的计算非常复杂和耗时，研究表明有简nAx(k+i) = J-iy$ - y(x(0) - H2Ax( k) Ax( k)n便的方法进行计算。将x,写成x= x(0)+ Ax涅于是xx = (x(0) + Ax)(x(0) +Ax ) = x(0)x(0) + x(0) Ax + x(0) Ax + Ax Axi jiijjijijjiy $ = ax (0) x(0)1x (0) x(0)2：+Ax(0) Ax x (0) Ax2：+ A1 (0) Ax x(0)2+ AAx AxA

15、x Ax2x (0) x(0) nx (0) AxnAx x (0)nAx Axn通过类比，可将泰勒展开式改写为y$ = y(x(0) + JAx + y(Ax)，迭代式改写为Ax(k+1) = J-iy$ -y(x(0)-y(Ax)，收敛判据max y (Ax(k+i)-y (Ax(k)1 2nc=c ,c ,L,c t =-y(Ax)12 n原方程组可简写为：f(x) = a +pb +p2c带入目标函数F (x) = f (x)2 = (a +财+3c )2 =0( p)，极值条件 TOC o 1-5 h z ii iid(p)_ dd pd pi=1i=1:+pb + p2C )2 i

16、ii2L(a + pb + p2c )(b + 2pc ) = 0i ii ii -i=1采用带有最优乘子的牛顿潮流算法后，潮流计算永远不会发散，即从算法上保证了计算 过程的收敛性，从而有效地解决了病态潮流的计算问题。最优潮流：基本潮流可归结为针对一定的扰动变量，根据给定的控制变量，求出相应的 状态变量，从而确定系统的一个运行状态。但基本潮流不能解决以下问题：当系统的状态变 量超出了它们的运行条件限制时，没有简便的手段使其恢复正常；当系统安全运行的方式很 多时，无法得到其中最经济的一种。最优潮流定义：就是当系统的结构参数及负荷情况给定时，通过控制变量的优选，所找 到的能满足所有指定的约束条件，

17、并使系统的某一个性能指标或目标函数达到最优时的潮流 分布。1）最优潮流与基本潮流的不同点：2）基本潮流的控制变量是给定的，而最优潮流中的控制变量通过优选得到；3）最优潮流除了要满足潮流等式约束外，还必须满足大量的不等式约束条件；4）基本潮流计算是求解非线性代数方程组，而最优潮流是一个非线性规划问题；5）基本潮流仅仅完成计算功能，而最优潮流可以根据实际需要自动优选控制变量，具 有指导系统进行优化调整的决策功能。最优潮流与传统经济调度的区别：传统经济调度只对有功进行优化，虽然考虑了线损修 正，也只考虑了有功功率引起的线损优化，同时传统经济调度一般不考虑母线电压的约束, 对安全约束一般也难以考虑。最

18、优潮流除了对有功和耗量进行优化外，还对无功及网损进行 了优化。此外，最优潮流还考虑了母线电压的约束及线路潮流的安全约束。最优潮流的控制变量有：1）除平衡节点外，其它发电机的有功出力2）所有发电机节点及其具有可调无功补偿设备节点的电压模值（或无功出力）。3）分接头可调变压器的变比状态变量1）除平衡节点外，其它所有节点的电压相角。2）除所有发电机节点及其具有可调无功补偿设备节点之外，其它所有节点的电压模值。最优潮流的目标函数：（简化表示为：f = f （u, X）1）全系统的发电燃料总耗量（或总费用）2）有功网损等式约束：潮流的基本方程g （u, x） = 0不等式约束：1）有功电源出力上下限约束

19、2）可调无功电源出力上下限约束3）分接头可调变压器变比调整范围4）节点电压模值上下限约束5）元件中通过的最大电流或视在功率约束6）线路通过的最大潮流约束7）线路两端节点电压最优潮流的数学模型：min f （u, x）、us.t g（u, x） = 0 最优潮流计算是一个典型的有约束非线性规划问题h（u, x） 0采用不同的目标函数并选择不同的控制变量，再和相应的约束条件相结合，就可以构成 不同应用目的的最优潮流问题：1）目标函数采用发电燃料耗量（或费用）最小，以除去平衡节点以外的所有有功电源 出力及所有可调无功电源出力（或用相应的节点电压），还有带负荷调压变压器的 变比作为控制变量，则就是对有

20、功及无功进行综合优化的通常泛称的最优潮流问 题。2）若目标函数同1），仅以有功电源出力作为控制变量而将无功电源出力（或相应的 节点电压模值）固定，则就成为有功最优潮流。3）若目标函数采用系统的有功网损最小，将各有功电源固定而以可调无功电源出力（或相应节点电压模值）及调压变压器变比作为控制变量，则就成为无功优化潮流。以上三种是目前用得最多的最优潮流问题。最优潮流的几种算法:1）简化梯度算法（特点：收敛缓慢；罚函数处理不等式约束会使收敛性变差）2）牛顿算法（具有二阶收敛性，速度快；对初值敏感；要求函数二阶连续可微；海森 矩阵及其逆阵计算量大）3）解耦最优潮流算法（分解为有功优化和无功优化两个优化子

21、问题，是大系统协调优 化问题）直流潮流法：以上的潮流计算，都属于精确的潮流计算。直流潮流法是所有潮流算法中 最快的。用于系统规划、在线安全分析等对计算速度要求较高的场合。对于输电线ij，支路有功潮流P = V2g 皿（g cos。+ b sin0 ）ij i ij i j ij ij ij ij00j很小,忽略电且j 令cos0 = 1,sin0 =0.-0 . , V = V = 1,忽略支路对地电容,p = = , P =P =Zb(00 ) = B0+&0 ,ijX1j匚 ij 七 j11 iij jjjn sB =Z ,BB =- , P = B e,这是一个线性方程组，可以一次直接求

22、解饿ii x ij x0j i lJlJ到结果，因而计算速度非常快。在可不计支路的无功潮流后，一条交流网络的支路就可以看成是一条直流支路，两端电 压值为0i和0j，支路电阻等于支路电抗X，支路电流值为相应的有功功率P。随机潮流：以上都是确定性潮流计算，随机潮流就是把潮流计算的已知量都作为随机变 量来处理的一种潮流计算方法。三相潮潮流：以上都是针对三相对称系统而言，超高压输电线各相间存在有不对称的耦 合，用对称分量法进行分析已经失去了各需网络相互独立的特点，所以研究三相潮流，目前 较多采用abc相坐标系统而不用120对称分量坐标系。总结牛拉法仍然不失为最基本、最重要的一种算法，它是其他一些派生算法的基础，其快速 的收敛特性和良好的收敛可靠性,使它在单纯的潮流计算以及在优化、稳定等程序的应用中,