潮流计算的基本算法及使用方法.doc

潮流计算的基本算法及使用方法一、 潮流计算的基本算法1. 牛顿拉夫逊法11 概述牛顿拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿拉夫逊法的核心。牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。12 一般概念对于非线性代数方程组即 (11)在待求量的某一个初始计算值附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组 (12)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 (13)将和相加，得到变量的第一次改进值。接着再从出发，重复上述计算过程。因此从一定的初值出发，应用牛顿法求解的迭代格式为 (14) (15)上两式中：是函数对于变量的一阶偏导数矩阵，即雅可比矩阵；为迭代次数。由式（14）和式子（15）可见，牛顿法的核心便是反复形成求解修正方程式。牛顿法当初始估计值和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。13 潮流计算的修正方程运用牛顿拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（节点）电压方程为从而得 进而有 （16）式（16）中，左边第一项为给定的节点注入功率，第二项为由节点电压求得的节点注入功率。他们二者之差就是节点功率的不平衡量。现在有待解决的问题就是各节点功率的不平衡量都趋近于零时，各节点电压应具有的价值。由此可见，如将式（16）作为牛顿拉夫逊中的非线性函数，其中节点电压就相当于变量。建立了这种对应关系，就可列出修正方程式，并迭代求解。但由于节点电压可有两种表示方式以直角做表或者极坐标表示，因而列出的迭代方程相应地也有两种，下面分别讨论。131 直角坐标表示的修正方程节点电压以直角坐标表示时，令、，且将导纳矩阵中元素表示为，则式（17）改变为 （17）再将实部和虚部分开，可得 （18）这就是直角坐标下的功率方程。可见，一个节点列出了有功和无功两个方程。对于节点（），给定量为节点注入功率，记为、，则由式（28）可得功率的不平衡量，作为非线性方程 （19）式中、分别表示第节点的有功功率的不平衡量和无功功率的不平衡量。对于节点（），给定量为节点注入有功功率及电压数值，记为、，因此，可以利用有功功率的不平衡量和电压的不平衡量表示出非线性方程，即有 （110）式中为电压的不平衡量。对于平衡节点（），因为电压数值及相位角给定，所以也确定，不需要参加迭代求节点电压。因此，对于个节点的系统只能列出个方程，其中有功功率方程个，无功功率方程个，电压方程个。将式（19）、式（110） 非线性方程联立，称为个节点系统的非线性方程组，且按泰勒级数在、（）展开，并略去高次项，得到以矩阵形式表示的修正方程如下 （111）上式中雅可比矩阵的各个元素则分别为 将（111）写成缩写形式 （112）对雅可比矩阵各元素可做如下讨论：当时，对于特定的，只有该特定点的和是变量，于是雅可比矩阵中各非对角元素表示为 当时，雅可比矩阵中各对角元素的表示式为由上述表达式可知，直角坐标的雅可比矩阵有以下特点：1) 雅可比矩阵是阶方阵，由于、等等，所以它是一个不对称的方阵。2) 雅可比矩阵中诸元素是节点电压的函数，在迭代过程中随电压的变化而不断地改变。3) 雅可比矩阵的非对角元素与节点导纳矩阵中对应的非对角元素有关，当中的为零时，雅可比矩阵中相应的、也都为零，因此，雅可比矩阵也是一个稀疏矩阵。132 极坐标表示的修正方程在牛顿拉夫逊计算中，选择功率方程作为非线性函数方程，把式中电压向量表示为极坐标形式则节点功率方程变为将上式分解成实部和虚部这就是功率方程的极坐标形式，由此可得到描述电力系统的非线性方程。对于节点，给定了 （113）对于节点，给定了、，而未知，式（113）中将失去作用，于是节点仅保留方程，以求得电压的相位角。（114）对于平衡节点，同样因为、已知，不参加迭代计算。将式（113）、式（114）联立，且按泰勒级数展开，并略去高次项后，得出矩阵形式的修正方程 （115）雅可比矩阵终，对节点，仍可写出两个方程的形式，但其中的元素以零元素代替，从而显示了雅可比矩阵的高度稀疏性。式中电压幅值的修正量采用的形式，并没有什么特殊意义，仅是为了雅可比矩阵中各元素具有相似的表达式。雅可比矩阵的各元素如下将式（115）写成缩写形式（116）以上得到了两种坐标系下的修正方程，这是牛顿拉夫逊潮流计算中需要反复迭代求解的基本方程式。2. 快速分解法21 概述快速分解法的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率的迭代分开来进行。快速分解法根据电力系统实际运行状态的物理特点，对牛顿-拉夫逊法潮流计算的数学模型进行合理的简化。22 基本公式在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化主要受电压相位的影响，无功功率的变化主要受母线电压幅值变化的影响。在修正方程式的系数矩阵中，偏导数和的数值相对于偏导数和是相当小的，作为简化的第一步，可以将方程式（21）中的子块和略去不计，即认为它们的元素都等于零。这样，阶的方程式便分解为一个阶和一个阶的方程式，即将式（21）简化为式（22）和式（23）。 (21) (22) (23)上述的简化大大地节省了计算机的内存和解题时间，但是矩阵和的元素都是节点电压幅值和相角差的函数，其数值在迭代过程中是不断变化的。因此，快速分解法潮流计算的第二个简化，也是最关键的一步简化就在于把系数矩阵和简化成在迭代过程中不变的常数对称矩阵。在一般情况下，线路两端电压的相角差是不大的（通常不超过）因此可以认为 ， （24）此外，与系统各节点无功功率相适应的导纳必远小于该节点自导纳的虚部，即 或 考虑到上面的关系，矩阵和的元素的表达式便被简化为 （i，j=1，2，n-1） （25） (i，j=1，2，m) （26） （27） （28）将式（27）和式（28）分别代入式（22）和（23），便得到：用和分别左乘以上两式便得简化了的修正方程式，可展开写成： （29） （210）式（29）和式（210）就是快速分解法潮流计算的修正方程式，其中系数矩阵都是由节点导纳矩阵的虚部构成，只是阶次不同，矩阵为阶，不含平衡节点对应的行和列，矩阵为阶，不含平衡节点和节点对应的行和列。 （211） （212）修正方程式（29）和（210）与功率误差方程式（211）和（212）构成了快速分解法迭代的基本计算公式。23 快速分解法的特点快速分解法与牛顿法潮流计算的主要差别表现在它们的修正方程上。快速分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅克比矩阵进行了有效的简化和改进，得到式（29）、式（210）所示的修正方程式。这两组方程式和牛顿法的修正方程相比主要有三个特点：a）快速分解法的修正方程式用两个阶线性方程组代替了一个阶线方程组。b）快速分解法的修正方程式中系数矩阵的所有元素在迭代过程中维持常数不变。c）快速分解法的修正方程式中系数矩阵是对称矩阵。这些特点在提高计算速度和减少内存方面的作用是很明显的：首先，因为修正方程式的系数矩阵是导纳矩阵的虚部，因此在迭代过程中不必像牛顿法那样每次都要重新计算雅克比矩阵，这样不仅减少了运算量，而且也大大简化了程序；其次，由于系数矩阵在迭代过程中维持不变，因此在求解修正方程式时，不必每次都对系数矩阵进行消去运算，只需要在进入迭代过程以前，将系数矩阵用三角分解形成因子表，然后反复利用因子表对不同的常数项或进行消去和回代运算，就可以迅速求得修正量，从而显著提高了迭代速度；第三，由于对称矩阵三角分解后，其上三角矩阵和下三角矩阵有非常简单的关系，所以在计算机中可以只存储上三角矩阵或下三角矩阵，从而也进一步节约了内存。快速分解法所采用的一系列简化假定只影响了修正方程的结构，也就是说只影响了迭代过程，但未影响最终结果。因为快速分解法和牛顿法都采用同样的数学模型，最后计算功率误差和判断收敛条件都是严格按照精确公式进行的，所以快速分解法和牛顿法一样都可以达到很高的精确度。为了改善快速分解法的收敛特性，修正方程的系数矩阵与一般并不简单的是电力系统导纳矩阵的虚部，下面讨论一下与的构成。与的阶数是不同的，为 阶，低于阶。因为式（210）不包含于节点有关的项，所以，如果系统有个节点，则应为阶。式（29）以有功功率误差为依据修正电压向量的角度，式（210）以无功功率误差依据修正电压幅值。为了加速收敛，使它们能够更有效地进行修正，可以考虑在中尽量去掉那些与有功功率及电压向量角度无关或影响较小的因素，而在中尽量去掉与无功功率及电压幅值影响较小的因素。所以，我们以电力系统导纳矩阵的虚部作为和时，可以在去掉充电电容和变压器变比的影响，在中去掉输电线路电阻对的影响。和的非对角元素和对角元素可分别按式（213）和（214）计算： （2.13） （2.14）式（213）中和分别为支路ij的电阻和感抗，式（214）中为节点接地支路的电纳。快速分解法改变了牛顿法迭代公式的结构，因此就改变了迭代过程的收敛性。牛顿法在迭代开始时收敛得较慢，当收敛到一定程度后，它的收敛速度非常之快，而快速分解法几乎是按同一速度收敛的，快速分解法每次迭代的计算量很小，因此快速分解法的计算速度比牛顿法有明显的提高。二、 潮流计算的使用方法1. 初始方式准备对任何潮流模拟操作计算，总是在某一个初始的运行方式上进行。这种初始方式可以是状态估计提供的实时运行方式，也可以是以往保存的历史运行方式。2. 调度操作模拟在准备好的初始潮流断面上，可以继续修改方式，模拟预想的潮流运行方式，再进行详细的潮流分析。模拟操作包括：1) 开关刀闸变位模拟2) 发电机功率调整3) 负荷功率设置4) 发电机分接头设置5) 线路停运、投入6) 变压器停运、投入7) 母线停运、投入8) 厂站停运、投入3. 运行参数维护潮流计算参数画面上可以设置算法、收敛判据、迭代次数、单/多平衡机等运行参数。平衡发电机是电气岛内的电压相角参考点，当采用“单平衡机”模式时，电网的不平衡功率（包括发电、负荷和网损）都将由设定平衡机吸收。当采用“多平衡机”模式时，电网的不平衡功率将由多台发电机负责平衡，多台发电机之间的不平衡功率分配方式包括容量、系数和平均三种方式。选择容量时将根据发电机的可调容量分配，选择系数时根据人工设置的系数按比例分配，选择平均时则平均分配不平衡功率。在分配过程中，确保发电机的出力在最大出力和最小出力范围内。发电机参数中可以设置发电机的调节特性，包括节点类型（平衡节点、PQ 节点、PV节点等），对于PV节点可以设定控制机端电压还是高压侧母线电压以及控制的目标电压值，对于按指定系数参与有功调节的机组可以设置比例系数。4. 计算结果分析潮流计算结束后，计算结果分析包括：1) 潮流计算状态2) 电气岛、迭代信息3) 潮流计算结果4) 设备越限和重载监视5) 运行信息5. 误差统计在潮流模拟计算完成后，如果现场很快发生了模拟的动作，可以从统计每个测点模拟计算值和实际量测值