2022年潮流计算的基本算法及使用方法 .pdf

.潮流计算的基本算法及使用方法一、潮流计算的基本算法1.牛顿拉夫逊法11 概述牛顿拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿拉夫逊法的核心。牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为 0，幅值为 1)启动即在此邻域内。12 一般概念对于非线性代数方程组0 xf即0,21nixxxfni,2,1 (11)在待求量x的某一个初始计算值0 x附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组0000 xxfxf (12)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量0100 xfxfx (13)将0 x和0 x相加，得到变量的第一次改进值1x。接着再从1x出发，重复上述计算过程。因此从一定的初值0 x出发，应用牛顿法求解的迭代格式为kkkxfxxf (14)kkkxxx1 (15)上两式中：xf是函数xf对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 12 页 -.次数。由式（1 4）和式子（15）可见，牛顿法的核心便是反复形成求解修正方程式。牛顿法当初始估计值0 x和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。13 潮流计算的修正方程运用牛顿拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（i节点）电压方程为njiijijUSUY1从而得njjijiiUYUS1进而有01jnjijiiiUYUjQP（16）式（16）中，左边第一项为给定的节点注入功率，第二项为由节点电压求得的节点注入功率。他们二者之差就是节点功率的不平衡量。现在有待解决的问题就是各节点功率的不平衡量都趋近于零时，各节点电压应具有的价值。由此可见，如将式（1 6）作为牛顿拉夫逊中的非线性函数0XF，其中节点电压就相当于变量X。建立了这种对应关系，就可列出修正方程式，并迭代求解。但由于节点电压可有两种表示方式以直角做表或者极坐标表示，因而列出的迭代方程相应地也有两种，下面分别讨论。13 1 直角坐标表示的修正方程节点电压以直角坐标表示时，令iiijfeU、jjjjfeU，且将导纳矩阵中元素表示为ijijijjBGY，则式（17）改变为01njjjijijiiiijfejBGjfejQP（17）再将实部和虚部分开，可得0011njjijjijijijjijiinjjijjijijijjijiieBfGefBeGfQeBfGffBeGeP（18）这就是直角坐标下的功率方程。可见，一个节点列出了有功和无功两个方程。对于PQ节点（1,21mi，），给定量为节点注入功率，记为iP、iQ，则由式（28）可得功率的不平衡量，作为非线性方程名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 12 页 -.njjijjijijijjijiiinjjijjijijijjijiiieBfGefBeGfQQeBfGffBeGePP11（19）式中iP、iQ分别表示第i节点的有功功率的不平衡量和无功功率的不平衡量。对于PV节点（nmmi,2,1），给定量为节点注入有功功率及电压数值，记为iP、iU，因此，可以利用有功功率的不平衡量和电压的不平衡量表示出非线性方程，即有22221iiiinjjijjijijijjijiiifeUUeBfGffBeGePP（110）式中iU为电压的不平衡量。对于平衡节点（mi），因为电压数值及相位角给定，所以SsSjfeU也确定，不需要参加迭代求节点电压。因此，对于n个节点的系统只能列出12 n个方程，其中有功功率方程1n个，无功功率方程1m个，电压方程mn个。将式（19）、式（110）非线性方程联立，称为n个节点系统的非线性方程组，且按泰勒级数在0if、0ie（mini,2,1）展开，并略去高次项，得到以矩阵形式表示的修正方程如下nnppnnnnnpnpnnnnnnnnnpnpnnnnpnpnpppppppppnpnppppppppnnppnnppnnppnnppnnppefefefefSRSRSRSRNHNHNHNHSRSRSRSRNHNHNHNHLJLJLJLJNHNHNHNHLJLJLJLJNHNHNHNHUPUPQPQP22112211221122112211222222222121222222222121111112121111111112121111222211（111）上式中雅可比矩阵的各个元素则分别为jiijfPHjiijePN名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 12 页 -.jiijfQJjiijeQLjiijfUR2jiijeUS2将（111）写成缩写形式efefSRLJNHUQPJ2（112）对雅可比矩阵各元素可做如下讨论：当ij时，对于特定的j，只有该特定点的if和ie是变量，于是雅可比矩阵中各非对角元素表示为iijiijjiijfGeBfPHiijiijjiijfBeGePNiijiijjiijeGfBfQJiijiijjiijeBfGeQL02jiijfUR02jiijeUS当ij时，雅可比矩阵中各对角元素的表示式为iiiiiinjjijjijjiijeBfGeBfGfPH1iiiiiinjjijjijjiijfBeGfBeGePN1fBeGfBeGfQJiiiiinjjijjijjiij1iiiiiinjjijjijjiijeBfGeBfGeQL1ijiijffUR22ijiijeeUS22由上述表达式可知，直角坐标的雅可比矩阵有以下特点：名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 12 页 -.1)雅可比矩阵是12 n阶方阵，由于jiijHH、jiijNN等等，所以它是一个不对称的方阵。2)雅可比矩阵中诸元素是节点电压的函数，在迭代过程中随电压的变化而不断地改变。3)雅可比矩阵的非对角元素与节点导纳矩阵BY中对应的非对角元素有关，当BY中的ijY为零时，雅可比矩阵中相应的ijH、ijN、ijJ、ijL也都为零，因此，雅可比矩阵也是一个稀疏矩阵。13 2 极坐标表示的修正方程在牛顿拉夫逊计算中，选择功率方程njjijiiiUYUjQP10作为非线性函数方程，把式中电压向量表示为极坐标形式iiijiijUeUUisincosjjjjjjjUeUUjsincos则节点功率方程变为0sincossincos1jjjnjijijiiiiijUjBGjUjQP将上式分解成实部和虚部0sincos1njijijijijjiiBGUUP0cossin1njijijijijjiiBGUUQ这就是功率方程的极坐标形式，由此可得到描述电力系统的非线性方程。对于PQ节点，给定了njijijijijjiiinjijijijijjiiiBGUUQQBGUUPP11cossinsincos121mi、（113）对于PV节点，给定了iP、iU，而iQ未知，式（113）中iQ将失去作用，于是PV节点仅保留iP方程，以求得电压的相位角。njijijijijjiiiBGUUPP1sincosnmmi、21名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 12 页 -.（114）对于平衡节点，同样因为sU、s已知，不参加迭代计算。将式（1 13）、式（1 14）联立，且按泰勒级数展开，并略去高次项后，得出矩阵形式的修正方程npnnnpnnnnpnppppppnpnpnpnpnpUUUUHHNHNHHHNHNHLJLJLJNHNHNHLJLJLJHHNHNHPPQPQP2221112211221122212121212221212121111212111111121211112211（115）雅可比矩阵终，对PV节点，仍可写出两个方程的形式，但其中的元素以零元素代替，从而显示了雅可比矩阵的高度稀疏性。式中电压幅值的修正量采用UU的形式，并没有什么特殊意义，仅是为了雅可比矩阵中各元素具有相似的表达式。雅可比矩阵的各元素如下ijijijijjijiijBGUUPHcossinnijjijijijijjiiiiiBGUUPH1cossinijijijijjijjiijBGUUUUPNsincosiiinijjijijijijjiiiiiiGUBGUUUUPN212sincosijijijijjijiijBGUUQJsincosnijjijijijijjiiiiiBGUUQJ1sincosijijijijjijiijBGUUQLsincosiiinijjijijijijjijiiiiBUBGUUUQL212cossin名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 12 页 -.将式（1 15）写成缩写形式UULJNHQP（116）以上得到了两种坐标系下的修正方程，这是牛顿拉夫逊潮流计算中需要反复迭代求解的基本方程式。2.快速分解法21 概述快速分解法的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率的迭代分开来进行。快速分解法根据电力系统实际运行状态的物理特点，对牛顿-拉夫逊法潮流计算的数学模型进行合理的简化。22 基本公式在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化主要受电压相位的影响，无功功率的变化主要受母线电压幅值变化的影响。在修正方程式的系数矩阵中，偏导数Q和VP的数值相对于偏导数VQ和P是相当小的，作为简化的第一步，可以将方程式（21）中的子块N和K略去不计，即认为它们的元素都等于零。这样，mn1阶的方程式便分解为一个1n阶和一个m阶的方程式，即将式（21）简化为式（22）和式（23）。VVLKNHQP1D (21)HP (22)VLVQ1D (23)上述的简化大大地节省了计算机的内存和解题时间，但是矩阵H和L的元素都是节点电压幅值和相角差的函数，其数值在迭代过程中是不断变化的。因此，快速分解法潮流计算的第二个简化，也是最关键的一步简化就在于把系数矩阵H和L简化成在迭代过程中不变的常数对称矩阵。在一般情况下，线路两端电压的相角差是不大的（通常不超过1020）因此可以认为1cosij，ijijijBG sin（24）此外，与系统各节点无功功率相适应的导纳LDiB必远小于该节点自导纳的虚部，即名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 12 页 -.iiiiLDiBVQB2或iiiiBVQ2考虑到上面的关系，矩阵H和L的元素的表达式便被简化为ijjiijBVVH（i，j=1，2，n-1）（25）ijjiijBVVL (i，j=1，2，m)（26）11,1122,1111,1111,222222121211,1121211111nnnnnnnnnnnnVBVVBVVBVVBVVBVVBVVBVVBVVBVH（27）mmmmmmmmmmmmVBVVBVVBVVBVVBVVBVVBVVBVVBV221122222212121121211111L（28）将式（27）和式（28）分别代入式（22）和（23），便得到：D1D1VBVPVBVQD2用11DV和12DV分别左乘以上两式便得简化了的修正方程式，可展开写成：1122111,12,11,11,222211,11211112211nnnnnnnnnnVVVBBBBBBBBBVPVPVP（29）mmmmmmmmmVVVBBBBBBBBBVQVQVQ212122221112112211（210）式（29）和式（210）就是快速分解法潮流计算的修正方程式，其中系数矩阵都是由节点导纳矩阵的虚部构成，只是阶次不同，矩阵B为1n阶，不含平衡节点对应的行和名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 12 页 -.列，矩阵B为m阶，不含平衡节点和PV节点对应的行和列。njijijijijjiisiisiBGVVPPPP1sincos（211）njijijijijjiisiisiBGVVQQQQ1)cossin(（212）修正方程式（29）和（210）与功率误差方程式（211）和（212）构成了快速分解法迭代的基本计算公式。23 快速分解法的特点快速分解法与牛顿法潮流计算的主要差别表现在它们的修正方程上。快速分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅克比矩阵进行了有效的简化和改进，得到式（29）、式（210）所示的修正方程式。这两组方程式和牛顿法的修正方程相比主要有三个特点：a）快速分解法的修正方程式用两个n阶线性方程组代替了一个n2阶线方程组。b）快速分解法的修正方程式中系数矩阵的所有元素在迭代过程中维持常数不变。c）快速分解法的修正方程式中系数矩阵是对称矩阵。这些特点在提高计算速度和减少内存方面的作用是很明显的：首先，因为修正方程式的系数矩阵是导纳矩阵的虚部，因此在迭代过程中不必像牛顿法那样每次都要重新计算雅克比矩阵，这样不仅减少了运算量，而且也大大简化了程序；其次，由于系数矩阵在迭代过程中维持不变，因此在求解修正方程式时，不必每次都对系数矩阵进行消去运算，只需要在进入迭代过程以前，将系数矩阵用三角分解形成因子表，然后反复利用因子表对不同的常数项VP/或VQ/进行消去和回代运算，就可以迅速求得修正量，从而显著提高了迭代速度；第三，由于对称矩阵三角分解后，其上三角矩阵和下三角矩阵有非常简单的关系，所以在计算机中可以只存储上三角矩阵或下三角矩阵，从而也进一步节约了内存。快速分解法所采用的一系列简化假定只影响了修正方程的结构，也就是说只影响了迭代过程，但未影响最终结果。因为快速分解法和牛顿法都采用同样的数学模型，最后计算功率误差和判断收敛条件都是严格按照精确公式进行的，所以快速分解法和牛顿法一样都可以达到很高的精确度。为了改善快速分解法的收敛特性，修正方程的系数矩阵B与B一般并不简单的是电力系统导纳矩阵的虚部，下面讨论一下B与B的构成。B与B的阶数是不同的，B为1n阶，B低于1n阶。因为式（210）不包含名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页，共 12 页 -.于PV节点有关的项，所以，如果系统有r个PV节点，则B应为1rn阶。式（29）以有功功率误差为依据修正电压向量的角度，式（210）以无功功率误差依据修正电压幅值。为了加速收敛，使它们能够更有效地进行修正，可以考虑在B中尽量去掉那些与有功功率及电压向量角度无关或影响较小的因素，而在B中尽量去掉与无功功率及电压幅值影响较小的因素。所以，我们以电力系统导纳矩阵的虚部作为B和B时，可以在B去掉充电电容和变压器变比的影响，在B中去掉输电线路电阻对B的影响。B和B的非对角元素和对角元素可分别按式（213）和（2 14）计算：22ijijijijxrxBijijijijiixrxB22（2.13）ijxB1ijiijbxB01（2.14）式（213）中ijr和ijx分别为支路ij的电阻和感抗，式（214）中0ib为节点i接地支路的电纳。快速分解法改变了牛顿法迭代公式的结构，因此就改变了迭代过程的收敛性。牛顿法在迭代开始时收敛得较慢，当收敛到一定程度后，它的收敛速度非常之快，而快速分解法几乎是按同一速度收敛的，快速分解法每次迭代的计算量很小，因此快速分解法的计算速度比牛顿法有明显的提高。二、潮流计算的使用方法1.初始方式准备对任何潮流模拟操作计算，总是在某一个初始的运行方式上进行。这种初始方式可以是状态估计提供的实时运行方式，也可以是以往保存的历史运行方式。2.调度操作模拟在准备好的初始潮流断面上，可以继续修改方式，模拟预想的潮流运行方式，再进行详细的潮流分析。模拟操作包括：1)开关刀闸变位模拟2)发电机功率调整3)负荷功率设置4)发电机分接头设置5)线路停运、投入名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 12 页 -.6)变压器停运、投入7)母线停运、投入8)厂站停运、投入3.运行参数维护潮流计算参数画面上可以设置算法、收敛判据、迭代次数、单/多平衡机等运行参数。平衡发电机是电气岛内的电压相角参考点，当采用“单平衡机”模式时，电网的不平衡功率（包括发电、负荷和网损）都将由设定平衡机吸收。当采用“多平衡机”模式时，电网的不平衡功率将由多台发电机负责平衡，多台发电机之间的不平衡功率分配方式包括容量、系数和平均三种方式。选择容量时将根据发电机的可调容量分配，选择系数时根据人工设置的系数按比例分配，选择平均时则平均分配不平衡功率。在分配过程中，确保发电机的出力在最大出力和最小出力范围内。发电机参数中可以设置发电机的调节特性，包括节点类型（平衡节点、PQ 节点、PV 节点等），对于 PV 节点可以设定控制机端电压还是高压侧母线电压以及控制的目标电压值，对于按指定系数参与有功调节的机组可以设置比例系数。4.计算结果分析潮流计算结束后，计算结果分析包括：1)潮流计算状态2)电气岛、迭代信息3)潮流计算结果4)设备越限和重载监视5)运行信息5.误差统计在潮流模拟计算完成后，如果现场很快发生了模拟的动作，可以从统计每个测点模拟计算值和实际量测值相比的误差，并统计出全网平均误差，统计方法遵循实用化考核细则。在表格中全部列出所有测点的SCADA 量测值、潮流模拟计算值、考核基准值以及测点误差等内容。如果只关心部分厂站的误差情况，局部误差统计中选择好需要关心的厂站，局部误差统计程序会过滤掉没有选中的厂站，只列出选中厂站的所有设备的误差统计情况。如果想要保存当前误差统计的运行断面，可以保存误差断面，将当前的断面以及误差统计结果一起保存起来。误差统计历史记录包含所有保存的误差断面，内容包括统计时间、平名师资料总结-精品资料欢迎下载-名师精心整理-第 11 页，共 12 页 -.均误差、执行用户、值班主机、断面名称以及操作信息等内容。名师资料总结-精品资料欢迎下载-名师精心整理-第 12 页，共 12 页 -