《潮流计算作业A4要点.pdf》由会员分享,可在线阅读,更多相关《潮流计算作业A4要点.pdf(15页珍藏版)》请在三一文库上搜索。
1、电力系统潮流计算综述 学院:电气工程学院 专业:电力系统及其自动化 学号: s13080802015 姓名:张雪 - 1 - 摘 要 电力系统潮流计算是电力系统分析中最基本的一项计算。本文对电力系统潮 流计算进行了综述。 首先简单回顾了潮流计算的发展历史,对当前基于计算机的 各种潮流算法的原理及其优缺点,作了简要介绍和比较, 并介绍了它们采用的一 些特别技术及程序设计技巧; 接着简要分析了三种新型的潮流计算方法的计算原 理及优缺点, 它们分别是基于人工智能的潮流计算方法、基于L1范数和现代内点 理论的电力系统潮流计算方法、基于符号分析的潮流计算方法等。除此之外还介 绍了配电系统潮流计算算法。
2、关键词 :电力系统;潮流计算;综述;新型潮流计算方法;配电系统 1 概述 电力系统潮流计算是研究电力系统稳态运行的一项基本运算。它根据给定系 统的网络结构及运行条件来确定整个系统的运行状态:主要是各节点电压(幅值 和相角 ),网络中功率分布及功率损耗等。它既是对电力系统规划设计和运行方 式的合理性、 可靠性及经济性进行定量分析的依据,又是电力系统静态和暂态稳 定计算的基础。 潮流计算经历了一个由手工,利用交、直流计算台到应用数字电 子计算机的发展过程。现在的潮流算法都以计算机的应用为前提。1956年ward 等人编制成实用的计算机潮流计算程序,标志着电子计算机开始在电力系统潮流 计算中应用。基
3、于导纳矩阵的高斯 塞德尔法是电力系统中最早得到应用的潮流 计算方法。因它对病态条件(所谓具有病态条件的系统是指:重负荷系统;包含 有负电抗支路的系统; 具有较长辐射型线路的系统; 长线路与短线路接在同一节 点,且其长度比值又很大的系统;或平衡节点位于网络远端的系统)特别敏感, 又发展了基于阻抗阵的高斯塞德尔法, 但此法中阻抗阵是满阵占大量内存,而 限制了其应用。 1961年VanNes等人提出用牛顿法求解系统潮流问题,经后人的 不断改进,而得到广泛应用并出现了多种变型以满足不同的需要,如快速解耦法、 直流法、保留非线性算法等。 同时,60年代初开始出现运用非线性规划的最优潮 流算法。 60年代
4、末 Dom-8mel和Tinney提出最优潮流的简化梯度法,70年代有人 提出海森矩阵法, 80年代SunDl提出最优潮流牛顿算法,还可把解耦技术应用于 最优潮流 ,从而形成解耦型最优潮流牛顿算法,还可把解祸技术应用于最优潮流, 从而形成解耦型最优潮流牛顿算法。随着直流输电技术的发展, 交直流联合电力 系统的潮流计算方法相应出现。另外,其它各种潮流算法如最小化潮流算法、随 机潮流算法等也不断涌现。 至于用于特殊用途的潮流算法如谐波潮流、适于低压 配电网的潮流算法也得到了较快的发展。 潮流算法多种多样,但一般要满足四个基本要求:(i)可靠收敛; (ii) 计算速 - 2 - 度快;(iii) 使
5、用方便灵活; (iv)内存占用量少。它们也是对潮流算法进行评价的主 要依据。 在潮流计算中, 给定的量应该是负荷吸收的功率、发电机发出的功率或者发 电机的电压。这样,按照给定量种类的不同,可以将节点分为以下三类 1 : (1)PQ节点。给定节点的注入有功功率P和注入无功功率 Q。这类节点对应于 实际系统中纯负荷节点 (如变电所母线 )、 有功和无功功率都给定的发电机节点(包 括节点上带有负荷 ),以及联络节点 (注入有功和无功功率都等于零)。这类节点占 系统中的绝大多数,它们的节点电压有效值和相位未知。 (2)PV节点。给定节点的注入有功功率P和节点电压有效值 U,待求量是节点 的注入无功功率
6、 Q和电压的相位 。这类节点通常为发电机节点,其有功功率给 定而且具有比较大无功容量, 它们能依靠自动电压调节器的作用使母线电压保持 给定值。有时将一些装有无功补偿设备的变电站母线也处理为PV节点。 (3)平衡节点。在潮流计算中,必须设置一个平衡节点,其电压有效值为给 定值,电压相位为 =0,即系统中其它各点的电压相位都以它为参考;而注入的 有功功率和无功功率都是待求量。实际上,由于所有的PQ节点和 PV节点的注入 有功功率都已经给定, 而网络中的总有功功率损耗是未知的,因此平衡节点的注 入有功功率必须平衡全系统的有功功率和有功损耗而不能加以给定。 需要注意的是以上介绍的节点分类只是一般的原则
7、,而不是一成不变的。 2 潮流计算主要方法与评价 2.1 潮流计算问题的数学模型 电力系统潮流的基本方程为 2 : * 1 n ii ijj j i PjQ Y U U (i=1,2,3 n) (1) 或 * 1 n ji iij j j PjQ UZ U (i=1,2,3 n) (2) 其中, ij Y , ij Z 分别为节点导纳矩阵和节点阻抗矩阵的相应元素,n为系统节点数。 这就是潮流计算问题最基本的方程式, 是一个以节点电压 U 为变量的非线性 代数方程组。 由此可见,采用节点功率作为节点注入量是造成方程组呈非线性的 根本原因。由于方程组为非线性的,因此必须采用数值计算方法,通过迭代来
8、求 解。根据在计算中对这个方程组的不同应用和处理,就形成了不同的潮流算法。 对于电力系统中的每个节点,要确定其运行状态,需要四个变量:有功注入 P、无功注入 Q、电压模值 U及电压相角 。n个节点总共有 4n个运行变量要确定。 再观察式 (1)和式(2),总共包括 n个复数方程式,如果将实部与虚部分开,则形成 - 3 - 2n个变量作为已知量而预先给以指定。也即对每个节点, 要给定其两个变量的值 作为已知条件,而另两个变量作为待求量。 按照电力系统的实际运行条件,根据预先给定的变量的不同, 电力系统中的 节点又可分为 PQ节点、PV节点及 V 节点或平衡节点三种类型。 对应于这些节点, 分别对
9、其注入有功、 无功功率,有功功率及电压模值以及电压模值和相角加以指 定;并且对平衡节点来说,其电压相角一般作为系统电压相角的基准(即=0)。 交流电力系统中的复数电压变量可以用两种坐标形式来表示 i j ii UU e(3) 或 iii Uejf(4) 而复数导纳为 ijijij YGjB(5) 将(3)、 式(4)以及式(5)代入以导纳矩阵为基础的式(1),并将实部与虚部分开, 可得到以下两种形式的潮流方程。 潮流方程的直角坐标形式为: (G eB f )() iiijiijjiijjijj j ij i PefG fB e(6) ( i1, 2 , n (G eB f )() iiijii
10、jjiijjijj j ij i QfeG fB e(7) (i1,2,n) 潮流方程的极坐标形式为: U ( G c o sB si n iiji ji ji ji j j i PU(8) ( i1 , 2 , n U (G sinB cos) iijijijijij j i QU(9) ( i1 , 2 , n 以上各式中,ji表示 号后的标号为 j节点必须直接和节点 i相联,并包括 j=i 的情况。这两种形式的潮流方程统称为节点功率方程,是牛顿-拉夫逊等潮流 算法所采用的主要数学模型。 对于以上潮流方程中的有关运行变量,还可以按其性质的不同加以分类,这 对于进行例如灵敏度分析以及最优潮流
11、的研究等都是比较方便的。 每个节点的注入功率是该节点的电源输入功率 Gt P 、 Gt Q 和负荷需求功率 Li P 、 Li Q 的代数和。负荷需求的功率取决于用户,是无法控制的,所以称之为不可控 - 4 - 变量或扰动变量。 而某个电源所发的有功、 无功功率则是可以由运行人员控制或 改变的变量,是自变量或称为控制变量。至于各个节点的电压模值或相角,则属 于随着控制变量的改变而变化的因变量或状态变量;当系统中各个节点的电压模 值及相角都知道以后,则整个系统的运行状态也就完全确定了。若以p、u、x分 别表示扰动变量、 控制变量、状态变量,则潮流方程可以用更简洁的方式表示为: f(x,u,p)=
12、0 (10) 根据式 (10),潮流计算的含义就是针对某个扰动变量p,根据给定的控制变 量u,求出相应的状态变量 x。 电力系统的潮流计算需要求解一组非线性代数方程。目前求解非线性代数方 程一般采用的是迭代方法, 而应用电子数字计算机进行迭代计算可以得到非常精 确的结果。常用的潮流算法有牛顿拉夫逊法、快速解耦法(PQ分解法)、直 流潮流法、极小化潮流算法、最优潮流算法、保留非线性法、交直流潮流法等。 2.2 牛顿-拉夫逊法 牛顿拉夫逊法简称牛顿法, 是求解非线性代数方程的一种有效且收敛速度 快的迭代计算方法, 而形成雅可比矩阵和求解修正方程式是牛顿法潮流计算中的 主体。牛顿 拉夫逊法将潮流方程
13、 f(x)=0用泰勒级数展开,并略去二阶及以上高 阶项,然后求解。 其实质是逐次线性化, 求解过程的核心是反复形成并求解修正 方程。其迭代格式为: (k)(k)(k) (k 1)(k)(k) (X)X(X)ff XXX (11) 式中: (X)f是(X)f对于变量 X的一阶偏导数矩阵; k为迭代次数。 各种形式牛顿法的共同优点是:(1)收敛速度快,具有平方收敛特性,其迭 代次数与系统规模基本无关;(2)能求解大部分有病态条件的问题;(3)利用了保 持稀疏性技术,所需内存适中。它是六十年代以来,广泛应用的方法。但具有以 下缺点: (1)由于雅可比矩阵的维数约为节点总数的两倍而且在迭代过程中不断
14、改变,因此在大规模电力系统中应用牛顿法计算潮流比较费时;(2)编程复杂;(3) 需要良好的初值 (可由高斯 -赛德尔法给出 ),否则不收敛或收敛到无法运行的解 上;(4)对重病态条件可能不收敛; 而快速分解法则是通过不断简化迭代过程中变化的矩阵来进行潮流计算,一 般来说,快速分解法所需要的迭代次数比牛顿法多,但每次迭代的计算工作量远 小于牛顿法, 因此总的来说迭代求解过程所需要的时间要少得多。直流潮流算法 是一种十分近似的方法, 它主要用于系统中有功功率分布的近似估算。极小化潮 流算法是将功率方程式的求解问题转化为一个求函数的极小值问题,然后应用数 - 5 - 学规划方法进行求解, 极小化潮流
15、算法的主要缺点是所需要的计算机内存计算时 间比常规牛顿法更多。 根据f(x) 的表达式不同,牛顿法又分功率偏差型算法和电流偏差型算法。根 据复电压变量采用的坐标形式不同,牛顿法又有直角坐标形式、 极坐标形式和混 合形式。 2.2.1 功率偏差型算法 令 i j ii UU e,可得极坐标形式修正方程式为: PH N QM LU U (12) 令 iii Uejf,可得直角坐标形式修正方程式为: 2 PHN e QML f R S U (13) 其特点是每次迭代, 雅可比矩阵都需要重新形成。雅可比矩阵虽非对称,但 分块雅可比矩阵为一高度稀疏阵。程序设计中采用的技巧有:(i)稀疏存储 (ii) 按
16、 行消去,采用边形成、边消元、边存储的方式(iii) 节点编号优化 (最优顺序消去 )。 2.2.2 电流偏差型算法 修正方程式由复电流偏差量构成,以直角坐标形式表示。当网络中只有PQ 节点和所有节点都为零注入的情况下,修正方程式为: r m IB Gf IGBe (14) 式中: r I 、 m I 分别代表节点 i电流偏差的实部和虚部。其优点是对存在许 多无负荷或轻负荷节点的网络,潮流计算有较好的收敛性, 可用于暂态稳定计算 中。 2.3 改进牛顿法 为了改进牛顿法在内存占用量及计算速度方面的不足,1974年提出快速解耦 法(又称P一Q分解法 ),是较成功的一种算法; 为了改进牛顿法在处理
17、病态条件时 的缺陷,提高算法的收敛性能,70年代后期开始采用一种保留非线性潮流算法, 将泰勒级数的高阶项也考虑进来,效果较好。 2.3.1 快速解耦法 - 6 - 它是密切结合电力系统固有特点,对牛顿法改进后得到的一种方法。原理是 根据系统有功主要决定于电压相角的变化,而无功主要决定于电压模值的变化这 一特性,并进行合理假设:(i)线路两端的相角差不大,且 ijij GB,即认为 cos1 ij ;sin ijijij GB ;(ii) 与节点无功功率对应的导纳 2 ii Q U远小于节点的 自导纳 ii B ,即 2 iiii QU B ,最后得修正方程式: P UB Q UBU (15)
18、式中: B、 B是由节点导纳阵的虚部构成的常数对称矩阵,但作了下述修 改:(i)在形成 B时,略去那些主要影响无功和电压模值的因素包括输电线的充电 电容及变压器非标准变比; (ii)在形成 B时,不计串联元件的电阻。 这种形成 B及 B的方案,又称为 XB方案。 此法特点包括: (i)以一个 (n一l)阶和一个 (n一m一1)阶方程代替了牛顿法的 (2n一m一2)阶方程,减少了内存需量及计算量。这里n为系统节点数; m为PV节 点数。(ii)用常数阵替代了随迭代过程变化的雅可比阵J,缩短了每次迭代的时间。 (iii) 用对称阵 B、 B替换了不对称阵 J,减少了三角分解的计算量并节约了内存。
19、因而,它具有简单、快速、内存节省且收敛可靠的优点,是广泛应用于在线处理 计算的方法,并已成为当前国内外最优先使用的算法。存在的问题是R/X比值过 大及因线路特别重载致使两节点间相角差特大时,收敛特性变坏或不收敛。 对大R/X比值病态问题的解决: (i)参数补偿,又分为串联补偿法和并联补偿 法两种; (ii) 算法改进,与上述传统的XB方案相反,在形成 B时采用精确的导纳 矩阵虚部,而在形成 B时只计串联元件的电抗值,并采用严格的,PQU交 替迭代方案。 2.3.2 保留非线性潮流算法 此类方法是将潮流方程组用泰勒级数展开,并保留高阶项。 (l)带二阶项的直角坐标形式牛顿算法 因直角坐标形式的潮
20、流方程为一个二次方程组,所以泰勒级数只取前三项即 为精确展开式。设潮流方程为(x)Y(x)Y0 s f。则由其泰勒级数展开式可得 1 2001 (X )JXH 2 s n XX XX YY XX (16) 式 中 : 0 J为 初 始 运 行 条 件 下 的 雅 可 比 矩 阵 ; H 为 一 常 数 矩 阵 ; - 7 - (0) 12 , T n xxxxxx为修正量向量。 程序设计技巧是泰勒展开式的第三项可写成和第一项相同的函数表达式仅 变量不同,以X代X。即有 (0) (x)JXY(X) s YY(17) 由式(17)得迭代公式: (k 1)1(0)(k) (X)YY(x) s XJY
21、(18) 此法主要特点是: (i)采用由初值 (0) X计算得的恒定雅可比阵,因而计算速度 较牛顿法快,但仍慢于快速解耦法;(ii) 处理病态的能力提高; (iii) 内存需量较牛 顿法大。 (2)带二阶项的直角坐标形式快速潮流算法 NagendraRao 等运用了两个技巧对计算进行简化:(i)改造导纳阵的对角元; (ii) 所有节点电压初值取为平衡节点电压,最后将雅可比矩阵化为常数对称矩阵。 对二阶项的计算,也充分利用了前一次的几个迭代量,避免了繁复计算。 此法主要特点是: (i)所需内存因雅可比阵的对称性大为减少;(ii)计算速度 比牛顿法快 40%50%,接近快速解耦法; (iii) 受
22、病态条件影响小,比快速解耦法 有更好的收敛可靠性。 2.4 最小化潮流算法 最小化潮流算法把潮流计算问题归结为求由潮流方程构成的一个目标函数 的最小值间题, 一般要采用数学规划方法或最小化技术。这种算法能有效的解决 病态电力系统的潮流计算间题,并已得到了广泛应用, 特别是可用于电力系统电 压稳定问题的研究。 采用最小化算法的优点是:(i)潮流计算永不会分散,(ii) 提供了给定运行条 件下,潮流问题是否有解的判断标志。 2.5 特殊性质的潮流算法 2.5.1 直流潮流 它对系统作了以下简化假定: ijij gb, ij 数值很小, ij UU ,略去线路 电阻及所有对地支路,并不计支路的无功潮
23、流。因而得其数学模型: 0 PB。 它的优点是计算速度最快;缺点是计算精度不高,且无法计算无功Q。 2.5.2 三相潮流 它是处理不对称电力系统的潮流计算方法。目前采用的方法主要有两种: (i) 相分量法。直接采用 abc相坐标系统,各元件以相参数表示,各已知或待求量均 - 8 - 较对称系统以三倍数增加。(ii) 序分量法 (或对称分量法 )。采用对称分量坐标,将 系统各量分为正、负、零序分量,并对系统中不对称元件的序分量之间的耦合, 通过加电流源补偿的方法使之解耦。 2.5.3 交直流联合电力系统的潮流计算 交直流联合电力系统的潮流计算是根据交流系统各节点给定的负荷和发电 情况,结合直流系
24、统指定的控制方式,通过计算来确定整个系统的运行状态它和 纯交流电力系统相比,有以下特点:(i)增加直流电力系统变量,与交流电力系统 变量通过换流站中交直流换流器建立联系;(ii)换流器一方面实现了交直流电力 系统间的有功功率传递,另一方面又从系统中吸取无功;(iij) 直流系统的运行须 对各个换流器的运行控制方式加以指定,直流系统的状态量是给定的直流控制量 和换流器交流端电压的函数。 主要有联合求解法和交替求解法两种计算方法,前者是将交流系统潮流方程 组和直流系统的方程组联立起来,统一求解出交流及直流系统中所有未知变量。 后者则将交流系统潮流方程组和直流系统的方程组分开来求解,求解直流系统方
25、程组时各换流站的交流母线电压由交流系统潮流的解算结果提供;而在进行交流 系统潮流方程组的解算时, 将每个换流站处理成接在相应交流节点上的一个等效 的有功、无功负荷,其数值则取自直流系统潮流的解算结果。这样交替迭代计算, 直到收敛。 2.5.4 随机潮流 把潮流计算的已知量和待求量都作为随机变量来处理,最后求得各节点电压 及支路潮流等的概率统计特性。此法最早是用直流模型, 后发展为线性化的交流 模型及采用最小二乘法并保留非线性的交流模型。其突出优点是通过一次计算就 提供了系统运行和规划的全面信息。 除上述之外, 还有其它一些用途不同的特殊潮流问题,如谐波潮流、动态潮 流等,在此不一一列举。 2.
26、6 潮流概念的推广 2.6.1 状态估计 实质是一种广义潮流计算。 一般潮流计算时, 已知量和方程式数等于未知量 数。而在状态估计中,已知量和方程式数大于待求未知量数,利用冗余变量,在 实际测量系统有偏差的情况下获得表征系统实际运行状态的状态量。主要方法有 最小二乘估计法、支路潮流状态估计法、递推状态估计法等。 - 9 - 2.6.2 最优潮流 所谓最优潮流, 就是当系统的结构参数及负荷情况给定时,通过控制变量的 优选,所找到的能满足所有指定的约束条件,并使系统的某一性能指标或目标函 数达到最优时的潮流分布.最优潮流计算是电力系统优化规划与运行的基础,它 将成为能量管理系统 (EMS)中的核心
27、应用软件之一。其数学模型可表示为: min(u,x) s.t. g(u,x)0 h(u, x)0 f (19) 选用不同的目标函数的控制变量,加上相应的约束条件, 就构成不同应用目 的的最优潮流问题。最优潮流的求解方法主要有: (1)最优潮流的简化梯度算法 此法采用了简化梯度, 并应用拉格朗日乘子和罚函数将等式和不等式约束加 在目标函数中, 从而把有约束问题变为无约束问题。优点是原理简单, 设计简便。 缺点是迭代点向最优点接近时走的是曲折路线,罚因子的选择比较困难。 (2)最优潮流的牛顿算法 对最优潮流问题: m i n( x ) s . t . g ( x )0 h ( x )0 f (20
28、) 先不考虑不等式约束,构造拉格朗日函数:(x, )f(x)g(x) T L,定义向 量, T Zx,则应用海森矩阵法求最优解点 * Z的迭代方程为: Wz=-d,式中: W,d分别为 L对于Z的海森矩阵及梯度向量。本方法的关键是充分开发并在迭代 过程中保持 W矩阵的高度稀疏性, 另外在求解时采用特殊的稀疏技巧。对不等式 约束的处理有两种方法:(i)罚函数法; (ii) 不等式约束化为等式方程法。 (3)解耦最优潮流 把最优潮流的整体最优化问题分解为有功优化和无功优化两个子优化问题。 它有一个特别的优点是容许根据两个子优化问题各自的特性而采用不同的求解 方法。 3 几种新型的潮流计算方法 3.
29、1 潮流计算的人工智能方法 近年来,人工智能作为一种新兴的方法,越来越广泛的应用到电力系统潮流 计算中。该方法不像传统方法那样依赖于精确的数学模型,这种方法只能基于对 自然界和人类本身活动的有效类比而获得启示。具有代表性的有遗传法、 模拟退 - 10 - 火法、粒子群优化算法等。 遗传算法是 20世纪80年代出现的新型优化算法, 近年来迅速发展, 它的机理 源于自然界中生物进化的选择和遗传,通过选择 (Selection)、杂交(Crossover)和变 异(Mutation)等核心操作,实现“优胜劣汰” 。遗传算法优点是具有很好的全局寻 优能力,优化结果普遍比传统优化方法好。缺点是计算量比较
30、大,计算时间长。 模拟退火算法是基于热力学原理建立的随机搜索算法,也可以视为一种进化 优化方法, 是一种有效的通用启发式随机搜索方法。算法思想来源于固体退火原 理:将固体加温至充分高温, 再让其徐徐冷却, 加温时固体内部粒子随温升变为 无序状态,内能增大,而徐徐冷却时粒子渐趋有序,在每个温度都达到平衡态, 最后在常温时达到基态,内能减为最小。其算法原理比较简单,只是对常规的迭 代寻优算法进行一点修正,允许以一定的概率接受比前次稍差的解作为当前解。 粒子群优化算法源自对群鸟捕食行为的研究,本质上属于迭代的随机搜索算 法,具有并行处理特征,鲁棒性好,易于实现。该算法原理上可以以较大的概括 找到优化
31、问题的全局最优解, 计算效率较高, 已成功地应用于求解电力系统中各 种复杂的优化问题。 人工智能方法的优点是: (1)与导数无关性。工程上很多优化问题的目标函数是不可导的,若采取前 一类方法只能对其进行假设和近似,这显然影响到解的真实性。 若采取非导数优 化方法,则不需要知道函数的导数信息,只依赖于对目标函数的重复求值运算; (2)随机性,容易跳出局部极值点,它们是一类全局优化算法,特别适用于 非线性大规模问题以及问题的解空间分布不规则的情况; (3)内在并行性,它的操作对象是一组可行解,而非单个可行解,搜索轨道 有多条,而非单条, 这种内在的可并行处理性大大提高了处理复杂优化问题的速 度,对
32、其内在并行性的开发可在一定程度上克服其性能上的不足。 其缺点是:表现不稳定,算法在同一问题的不同实例计算中会有不同的效果, 造成计算结果的可信度不高;按概率进行操作,不能保证百分之百获得最优解, 通常得到的解是与最优解很接近的次最优解;算法中的某些控制参数需要凭经验 人为地给出,需要一定量的试验或专家经验。 3. 2 基于 L1范数和现代内点理论的电力系统潮流计算方法 一般潮流计算采用迭代的计算方法。然而,这些直接迭代求解的方法有一个 共同的缺点:病态潮流计算问题。在一些病态电力系统的计算中,算法常常出现 振荡和不收敛的现象。 针对上述情况, 研究人员提出了基于非线性规划模型的算 法。该类算法
33、在数学上可表示为求一个由潮流方程构成的目标函数最小值问题。 在给定运行条件下,若潮流问题有解,则目标值为零;若潮流向题无解,则目标 - 11 - 值为一不为零的正值。因此,即使是在病态系统的情况下,计算过程不会发散。 国内专家学者对解决此问题也进行了许多有益的探讨。 提出了一种基于内点非线性规划的潮流计算模型和算法。基于L1范数的计算 原理,潮流方程的求解可以转化为求解一个新的非线性规划模型L1LF,并结合 现代内点算法来进行求解。和过去的模型相比,该模型非常的简洁、直观,易于 编程。与现代内点算法相结合的求解过程表现出了良好的收敛性和快速性,计算 结果准确、 可靠,计算各种病态系统均可良好的
34、收敛,基于L1范数的数学规划模 型将传统电力系统潮流的直接迭代求解转化为对一简单规划问题的求解后,对系 统运行中各部分的控制可更加简便。增加适当的不等式约束和相关控制变,即可 获得近似于最优潮流的计算模型,可方便的进行潮流计算中的调整。 3. 3 电力系统潮流计算的符号分析方法 随着电力系统规模的扩大, 电力系统的实时计算问题显得日益重要,但长期 以来受算法的计算效率所限, 潮流计算的速度难以得到实质性的突破。根据电力 网络在实际运行中的特点, 结合网络图论理论提出了运用符号分析方法求解电力 网络潮流的新思路,有望克服传统数值计算方法在收敛性、冗余项对消、计算机 有效字长效应等方面的不足。 基
35、于符号分析方法的潮流计算方法通过建立电力网络的拓扑模型生成拓扑 网络的全部树和 2-树,应用网络的 k-树树支导纳乘积对电力网络的节点电压方程 进行拓扑求解,进而得出所求变量 (即各节点电压 )的符号表达式 (即关于元件参数 符号的显式表达式 )。这种方法避免了求解非线性方程,不必进行行列式的展开 和代数余子式的计算, 而且不需要写出行列式和代数余子式,克服了传统数值计 算的不足。同时它还带来一个附加的好处,即在构造函数时自然地产生并行处, 以及由它的拓扑性质带来的电力网络运行方式改变后计算的灵活性。这些特点将 在电力系统的在线计算、静态安全分析等领域发挥明显优势。另外,传统的潮流 计算方法都
36、是纯 “ 数值计算 ” ,利用这些方法计算出来的结果是数字而不是函数, 它们的特点是逐点进行完整的数值计算,因此不可避免地存在收敛性问题、冗余 项对消问题、计算机有效字长效应问题和相近数值求差时发生的浮点运算误差问 题。 基于符号分析方法的潮流计算方法在电力系统在线静态安全分析、短路计 算、灵敏度计算等领域中也可推广使用。 4 配电系统潮流计算算法的研究 4.1 线性规划 - 12 - 线性规划是数学规划的一个重要分支,它在理论和实践上都比较成熟, 因而, 在二十世纪十年代以前,线性规划发展很快,在配电系统经济运行、 水库调度以 及物资合理调运等方面,都得到了应用。 线性规划法是在一组线性约束
37、条件下,寻找线性目标函数的最大值或最小值 的优化方法。 对于配电系统潮流计算问题, 线性规划方法一般将非线性方程和约 束使用泰勒级数近似线性化处理,或将目标函数分段线性化。 它用非负变量的线 性化形式来处理问题的目标函数和约束条件,线性规划解配电系统优化问题,是 将问题的目标函数和约束条件线性化,并把注意力集中在顶点, 有步骤地在顶点 中寻优,从而保证了最优值的唯一性。这是一个很重要的特性。 一般线性规划问题用矩阵表示如下: min 0 CX st AXb X (21) 其中A是m n矩阵, C是n维行向量, b是n维行列向量。 线性规划求解配电系统潮流计算又分为单纯形法和对偶单纯形法、混合整
38、数 规划法。 (1) 单纯形法和对偶单纯形法 求解线性规划问题最基本的方法是单纯形法和对偶单纯形法。单纯形法的基 本思路是从一个基本可行解出发,求一个使目标函数有所改善的基本可行解;通 过不断改进基本可行解,力图达到最优基本可行解。它是G. B. Dantzig 在1947年 提出来的,后来许多学者又对其进行了改进,如修正单纯形法。 修正单纯形法的 基本思路是给定初始基本可行解后,通过修改旧基的逆来获得新基的逆,进而完 成单纯形法的其他运算。 对偶单纯形的基本思路是从原问题的一个对偶可行的基 本解出发, 求改进的对偶可行的基本解, 当得到的对偶可行的基本解是原问题的 可行解时,就达到最优解。在
39、使用单纯形法和对偶单纯形法求解线性规划问题时, 每次迭代都要把整个表格重新计算一遍,不必要地计算了许多与迭代过程无关的 数据,从而使计算机的存贮量大,计算量也大。文献3研究了线性规划对偶单 纯形法的改进。 通过改进原始的单纯形法思想,建立了标准型线性规划对偶单纯 形法的一种改进算法。 与原对偶单纯形法相比, 改进算法的存贮量和计算量大大 减少。文献 4给出了一种可以避免退化情况的单纯型方法和对偶单纯型方法。 (2) 混合整数规划法 除目标函数和约束函数是线性函数外,决策变量部分是整数值, 则称此类问 题为混合整数规划。 在配电系统中既存在像发电机输出功率、节点电压等连续变 量,又存在像变压器变
40、比、 可调电容等离散变量, 因而潮流计算也属于混合整数 - 13 - 规划问题。这类方法求解的难点在于离散变量的处理。 对于潮流计算问题, 线性规划方法一般利用泰勒级数对非线性方程和约束条 件做近似线性化处理,或将目标函数分段线性化。线性化以后,求解可以用改进 的单纯形法或对偶线性规划法。这对于严格的凸函数优化线性规划方法很有效 果,但对有功无功耦合的目标函数优化,尤其是对以网损最小化为目标的优化效 果不好,加之在潮流计算问题中, 要考虑的等式约束方程, 即每个节点的有功和 无功功率注入平衡方程是典型的非线性方程,因此在耦合的潮流计算问题中较少 使用线性规划法求解。 但由于有功潮流可以以很好的
41、精度线性化,而配电系统经 济调度主要对发电厂有功进行分配,因此线性规划方法能够在安全约束经济调度 中广泛地应用。 4.2 非线性规划 一般的非线性规划问题可描述为满足非线性约束条件的非线性函数的最小 值问题,非线性规划是配电系统最优运行最早使用的一类最优化方法,因为它所 描述的结构与电网络的物理模型结构很相似。该方法解法较多, 很多在实际应用 中已用于解决实时在线和离线运行等问题。 非线性规划法是在等式约束条件或不等式约束条件下优化目标函数,其中等 式约束、不等式约束或目标函数至少有一个为非线性函数。它的一般形式可以写 成: min f(x) st AXb Exe (22) 其中f(x)为目标
42、函数, A为m n矩阵,E为l n矩阵,b和e分别为 m维和l维列向 量。在非线性规划方法中,最突出的几种方法是简化梯度法、牛顿法和内点法。 二次规划法是非线性规划方法中的特殊情况。 4.3 二次规划 二次规划是一种特定形式的非线性规划,其目标函数是二次的, 约束是线性 的。相对于非线性规划来说, 二次规划的形式比较简单, 但也可大致地反映配电 系统的物理特性,并且其Hessian矩阵是常数矩阵,一阶偏导数矩阵是线性的, 这对于解潮流计算是很有利的条件。此外,二次规划还可以转化为线性规划问题 来解算。这都使问题得以简化。 二次规划法与非线性规划还是有许多相似之处:精度比较高, 但对大型系统 的
43、收敛性比较差。而且在许多地方,二次规划法还不如非线性规划。 - 14 - 5 结束语 本文简要列举了几种常用及特殊的潮流计算方法,但到目前为止,牛顿-拉 夫逊仍然不失为最基本、最重要的一种算法。它是其它一些派生算法的基础,其 快速的收敛性和良好的收敛可靠性,使它在单纯的潮流计算以及在优化、稳定等 程序的应用中,继续占有重要的地位。其它的算法各有特点, 但算法之间都互相 补充。虽然经过三十多年的发展,潮流算法已比较成熟。但是,仍存在不少尚待 解决的问题,例如: (i)各种牛顿法潮流算法,对于某些病态条件可能导致不收敛 的问题; (ii) 潮流计算的多解现象及其机理;(iii) 在重负荷情况下,邻
44、近多根与电 压不稳定问题的关联。 由于最优潮流能够把系统的安全性和经济性融为一体,并 能够提供用于提高系统安全经济性能的决策依据;将成为现代能量管理系统的核 心应用软件之一。当前无论在实践上还是在理论上,均有许多问题需待解决,特 别是如何快速求解成千上万个变量的大规模非线性规划问题。 同时,本文还列举了几种新型的电力系统潮流算法,对比常用的牛顿-拉夫 逊法,快速解耦法等在收敛性、实用性、扩展性等方面都有了较大改进。但无论 在理论上还是在实践上,仍然还有许多问题需要解决。 随着各种新型算法的出现, 能否将新型算法和传统算法有机结合、 取长补短, 使得配电系统潮流计算进一步发展和完善,这将是今后研
45、究和探讨的立足点。 参考文献 1夏道止 .电力系统分析 .2版.中国电力出版社 ,2004.84-115 2诸骏伟 .电力系统分析 .上册.北京:水利电力出版社 ,1995.1-73 3赵晋泉 ,侯志俭 ,吴际舜 .改进潮流计算牛顿算法有效性的对策研究.中国电机 工程学报 ,1999.19(12),70-75. 4郝玉国 ,刘广一 ,于尔镫 .一种基于 Karmarkar内点法的潮流计算算法.中国电机 工程学报 ,1996,16(6),410 412. 5王守相 ,刘玉田 .电力系统潮流计算研究现状.山东电力技术 ,1996(5):8-12 6秦景 ,牛卢璐 ,路一平 .电力系统潮流计算的最新进展.河北建筑工程学院学 报,26(4):61-64. 7钟世霞 ,袁荣湘 ,内点法在电力系统中的应用评述.高电压技术 , 2005, 31 8阳育德 ,韦化.基于L1范数和现代内点理论的电力系统潮流计算.广西大学学报 , 2007, 32: 6 9罗日成 ,李卫国 .电力系统潮流计算的符号分析方法.电网技术 , 2005, 29: 10
链接地址:https://www.31doc.com/p-5210229.html