配电网无功功率优化研究

1、配电网无功功率优化研究摘要配电网的无功功率的有效优化与合理控制既能提高电力系统运行时 的电压质量，也能有效减少网损，节约能源，是保证电力系统安全经济运 行的重要措施，对电网调度和规划具有重要的指导意义。无功优化的核心问题主要集中在数学模型和优化算法两方面，其中数 学模型问题是根据解决问题的重点不同来选取不同的目标函数；而优化算 法的研究则大量集中在提高计算速度、改善收敛性能上。本文选取有功网 损最小作为数学模型的冃标函数，数学模型的约束条件有各节点的注入有 功、无功功率的等式约束和各节点电压、发电机输出无功功率、可调变压 器变比、并联补偿电容量、发电机机端电压均在各自的上下限之内的不等 式约束

2、，优化方法采用遗传算法。设计和编制了牛顿拉夫逊直角坐标 matlab潮流计算程序以及遗传算法无功优化的matlab潮流计算程序。通过 ieee30节点系统的算例分析，得出基于遗传算法的无功优化能有效降低系 统网损、提高电压水平，验证了该算法在解决多变量、非线性、不连续、 多约束问题时的独特优势，并指出了该算法的不足之处以及如何改善。关键词：牛顿拉夫逊法，无功优化，遗传算法research of reactive power optimization distribution networkabstractreactive power with reas on able optimization

3、 and control of powersystem can not only improve the stability of power system, but also effectively reduce network losses and save energy. it ensures the safety and economic operation of power systems and improve the voltage quality. it is important for planning departments on grid reactive power s

4、cheduling reactive power optimization focuses on mathematical models and optimization algorithms. the mathematical model is selected depending on the focus of problem-solving. optimization algorithm is concentrated in improving the calculation speed and improve the convergence performance this paper

5、 selects the active power loss minimum objective function as a mathematical model, the constraints of mathematical model are each node of the injected active and reactive power equality constraint and the node voltage and reactive power of generator output, adjustable transformer ratio, parallel cap

6、acitanee compensation, the generator terminal voltage within the respective upper and lower limits of the in equality con strai nts, optimizati on method using genetic algorithms. design cartesian coordinate newt on raphson power flow calculation method and genetic algorithm matlab calculate the rea

7、ctive power optimization procedures through a numerical example of the ieee 30 node system, we can draw reactive power optimization based on genetic algorithm can effectively reduce system loss and improve voltage level and verify the algorithm have unique advantages to solve multivariable, nonlinea

8、r, discontinuous, multi-constraint problem. key words: newton raphson method; reactive power optimization; genetic algorithm第一章绪1.1引x-.1.2配电网的特1.3无功优化的基本概31.4无功优化的数学模.31.5现状和发展趋3第二章基于牛顿拉夫逊法的潮流计42.1电力系统潮流计算方法概42.2节点导纳矩念型势算述阵2.3牛顿拉夫逊法的计2.4牛顿拉夫逊发的基本流程错误！未定义书签。第三章电力系统无功优化问题及其遗传算法优化求解133.1无功优化问题描述及其模型133

9、.2遗传算法的理论基础143.3遗传算法基本原理及操作过程153.3.1适应度函数定标153.3.2初始解的形成153.3.3遗传操作163.4基于遗传算法的无功优化与电压控制实现的步骤173.2遗传算法的流程图20第四章算例分析204.1ieee14节点系统214.2ieee14节点系统算例分析错误!未定义书签。4.3ieee30节点系统214.4ieee30节点系统算例分析错误!未定义书签。第五章结论与展望27参考文献33附录a遗传算法无功优化matlab程序34第一章绪论1.1引言电能是现今社会最主要的能源，人们工作生活中都离不开电能。随着 社会的不断发展，电能的重要性显著增加。提供安全

10、、可靠、稳定、环保 的电能是现今电力系统发展的首要目标。最优潮流被提岀以后就一直用于电力系统的经济和安全运行及规划。最优潮流是指当系统的结构参数和负 荷情况都已给定时，调节可利用的控制变量（如发电机输出功率、可调变压 器抽头等）来找到能满足所有运行约束的，并使系统的某一性能指标（如发 电成本或网络损耗）达到最优值下的潮流分布。这一大系统非线性规划问题, 通常分为两个子问题：调节发电机的有功出力以减少发电费用；调节pv 节点和平衡节点的电压及可调变压器的分接头位置以改善电压分布及减 少系统的有功网损，后者即为无功优化问题。电力系统无功优化控制是指在满足各种电力系统运行条件的约束下， 对系统进行尽

11、量少的无功补偿，使电力系统中的各个节点电压得到最大限 度的改善，系统的有功网损降低，达到提高电力系统运行稳定性与经济性 的冃的。它涉及选择无功补偿装置地点、确定无功补偿容量、调节变压器 分接头和发电机机端电压的配合等，是一个动态、多目标、多约束的非线 性规划问题，也是电力系统分析中的一个难题。无功功率的最优分布包括无功功率电源的最优分布和无功功率负荷的最优补偿两个方面。电力系统的无功优化和电压控制是相互作用的，合 理的无功潮流分布是维持电压稳定的前提。无功功率的流动将在电网中产 生压降，造成电力系统节点电压偏移。当节点处的无功功率过剩时，往往 意味着电压的升高，相反，当节点处的无功功率不足时，

12、常常会使电压水 平降低。电力系统无功优化与控制是保证电力系统安全经济运行、提高电压质 量的重要措施，对指导调度人员安全运行和计划部门进行电网规划具有重 要意义。电力系统无功优化与控制不仅能改善电压质量，提高电力系统运 行的稳定性，更能有效的减少网损，节约能源。因此研究无功优化与控制 问题具有重要意义。11.2配电网特点配电网具有以下显著特点：(1) 闭环设计'开环运行,一般呈辐射状分布;适合于独立进化优化计算。(2) 节点和线路都较多，接线复杂;有的变电站出线可达到二十多条。要 求算法计算速度快，能应用于大规模系统。因此，配电网无功优化从数学模型的建立到优化算法和优化方式的选 择都应适

13、应配电网特点。1.3无功优化的基本概念电力系统无功优化是指在电力系统有功负荷、有功电源及有功潮流分 布已经给定的情况下，以发电机端电压幅值、无功补偿电源容量和可调变压 器分接头位置作为控制变量，而以发电机无功岀力、负荷节点电压幅值和支 路输送功率作为状态变量，应用优化技术和人工智能技术，在满足电力系统 无功负荷的需求下，谋求合理的无功补偿点和最佳补偿容量，使电力系统安 全、经济地向用户供电。配电网自动化水平的不断提高，为实现无功优化控制提供了条件，也使 英成为当前迫切需要研究解决的问题。因此，本文将主要研究配电网的无功 运行优化问题。1.4无功优化的数学模型电力系统无功优化问题的数学模型包括目

14、标函数、功率方程约束、变量约束。无功优化的目标函数根据具体需要有很多种，从技术指标方面或经济 指标方面看各有侧重。常见的有(1) 电压质量最好；(2) 全网有功网损最小；(3) 电网新增加无功补偿容量最小；(4) 系统总的费用最少；(5) 控制变量变化次数最小等。本文主要研究运行时的无功优化问题，宜有功网损最小作为数学模型 的目标函数。1.5现状和发展趋势在无功优化问题这一研究领域内，已有多种解决方法，例如：线性规 划、非线性规2划、混合整数规划、灵敏度分析、遗传算法等。这些方法都有各自的 优越性，也冇一定程度的局限性。线性规划是比较成熟的，它速度快、收敛性好、算法稳定，但在处理 无功规划优化

15、时需要将目标函数和约束函数线性化，要求优化问题可微, 对离散性问题缺乏指导性；若迭代步长选取不合适，可能会引发振荡或收 敛缓慢。非线性规划能直接处理非线性的目标函数和约朿函数，但非线性 规划还没有一个成熟的算法，现有算法存在计算量大、收敛性差、稳定性 不好等问题。尽管基于灵敏度和梯度法的数学优化方法能用来解决电力系 统的无功优化问题，但与线性规划法同样要求假设控制量是连续的，而但通常只能求得局部最优解。混合整数规划可以较好地处理离散性整数问题, 但在实际中由于操作复杂而得不到推广应用。1967 年 j. d. bagley 首次提出了遗传算法(genetic algorithm,简称 ga)

16、的概念。1975年左右美国密执安大学教授john h. holland等研究出了具有 开创意义的遗传算法理论和方法。在研究遗传算法的专家学者中，d.e. goldberg的贡献最为突出。他不但建立并完善了整个ga体系，而且将其 应用到优化、搜索及机器学习等领域，为ga的发展拓展了天地。遗传算法把自然界中基于自然遗传和自然选择的机制引入到数学理 论中来，提出了一种全新的寻优算法。它是利用目标函数本身的信息建立 寻优方向，因此不要求函数的连续性和可导性，有能力在一个复杂的、多 极值点、具有不确定性的空间中寻找全局最优解遗传算法随着计算机技 术的高速发展已经引起越来越多的注意，并已经应用于求解许多领

17、域中的 难题。在许多情况下，遗传算法表现得优于传统的优化算法。近年来，遗传算法在搜索与最优化问题方面已取得较大的进展。在电力 系统技术中，这一应用已经覆盖了负荷预测,电力系统设计与规划，电力系统 的进度安排与调度,单位投入和其它电力系统控制问题。遗传算法在寻求电 力系统问题全局最优解方面是强有力的工具，并被广泛应用于最优化及数 学问题上。在电力系统研究中，遗传算法具有随机搜索、灵活高效、稳定、 多目标处理和对复杂因素进行处理等优点。3第二章基于牛顿一拉夫逊法的电力系统潮流计算2.1电力系统潮流计算方法概述电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂 电力系统止常和故障条件下稳态

18、运行状态的计算。潮流计算的目标是求取 电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统 各元件是否过负荷、各点电压是否满足要求，功率的分布和分配是否合理 以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规 划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。随着用数字计算机解电力系统潮流问题的开始，许多计算方法被应用 如：阻抗法，pq节点法，牛顿一拉夫逊法等。因为牛顿一拉夫逊法是数 学中解决非线性方程式的典型方法，有较好的收敛性，用matlab仿真时在 内存要求、速度方面都较好。所以，本文采用牛顿一拉夫逊法计算电力系 统的初始网损。2.2节点导纳矩阵节

19、点导纳矩阵对角线元素yii(i=l,2, -n)为自导纳。节点i的自导纳yii 在数值上就等于与该节点直接连接的所有支路导纳的总和。节点导纳矩阵 非对角线元素yij (i=l,2, n；冋2,n； ihj)为互导纳。节点i、j之间 的互导纳yij数值上就等于连接节点i、j的支路导纳的负值。显然，yij=yjio 假如两节点不直接相连，也不计两支路之间，比如两相邻电力线路之间的 互感时，yij二yji=o节点电压方程ib?ybubo(2-4)注：ib节点注入电流列向量；ub节点电压的列向量；yb nxn阶节点导纳矩阵。它可展开为4?i1?y11?i?y?2?21?i3?=?y31?i?y ?n?

20、 ?nly12y22y32?yn2y13y23y33?yn3?yln?ul?u?y2n?2?y3n?u3?(2-5)? ?ynn?un?根据定义求取节点导纳矩阵时，仅需注意以下几点：(1) 节点导纳矩阵是方阵，一般也是对称矩阵。这是由网络的互易特性 所决定的。通常情况下取大地编号为零，作为参考节点。(2) 节点导纳矩阵是稀疏矩阵，每行非零非对角元素数就等于该行所对 应节点所连接的不接地支路数。节点导纳矩阵的对角线元素等于各该节点所连接导纳的总和。因此, 与没有接地支路的节点对应的行中，对角元为非对角元之和的负值。(4)网络中的变压器，设变压器两侧线路的阻抗都未经归算，即分别为i (高压测)、i

21、i (低压侧)线路的实际阻抗，变压器本身的阻抗归在低 压侧，设变压器变比为k(高、低压绕组电压之比)。2.3牛顿-拉夫逊法潮流计算方法牛顿-拉夫逊法法是常用的解非线性方程组的方法，也是当前广泛采用的计算潮流的方法，其原理如下。设有非线性方程组?fl(xl,x2,.zxn)?yl?f(x,x,.,x)?y?212n2? ?fn(xl,x2/.,xn)?y n(2?6)(0)(0)其近似解为xl(0), x2，xn,与精确解分别相差?xl, ?x2,，?xn, 则下式成立(0)?fl(xl(0)?xl,x2?x2/./x(0)?xn)?yl?(0)(0)(0)?f2(xl?xl,x2?x2/./x

22、n?xn)?y2? ?f(x(0)?x,x(0)?x,./x(0)?x)?y12nn2n?nl(2?7)上式中的任何一式都可按泰勒级数展开。以下则以第一式为例子加以说明,5(0)(0)(0)fl(xl(0)?xl,x2(0)?x2,.,x(0)?x)?f(x,x,./x)?nlnl2 n?fl?f?f?xl?l?x2?.?l?xn?l?yl?xl0?x20?xn0式中：?fl?f?f(o)(o),1,，1分别表示以xl(o), x2,，xn代入这些偏导数表示式?xl0?x20?xn0(0)(0)的计算所得，？1则是一包含xl(0), x2,，xn的高次方与fl的高阶 偏导数乘积的函数。如近似解

23、?xi(0)与精确解相差不大，贝ij?xi的高次方可略略去，从而?1 也可略去。由此可得?fl?f?f(0)(0)(0)?xl?l?x2?.?l?xn?yl?fl(xl,x2/./xn)?xl0?x20?xn0 ?f2?f2?f2(0)(0)(0)?f(x,x,./x)?x?x?.?xn?y2?21122n?x?x?x?1020n0 ? ?fn?fn?fn(0)?fn(xl(0),x2/.,x(0)?x?x?.?xn?ynl2 n?x?x?x?1020n0?(2?8)这是一组线性方程或线性化了的方程组，成称为修正方程组。它可改写为如下的矩阵方程：?fl?fl?f.1?x?x?xn(0)(0)?

24、yl?fl(xl(0),x2,./xn)?1020?f?f2?f(0)(0)(0).2?y2?f2(xl,x2,./xn)?2?xl0?x20?xn ?y?f(x(0),x(0)/.,x(0)?n 2n?nl?f?f?fn?n n?xl0?x20?xn?xl?x?2?0?.?x?n?0?(2?9)或简写为：?f?j?x其中：j称为函数fi的雅克比矩阵；？x为由?xi组成的列向量；？f则称 不平衡量的列向量。将xi(o)代入，可得?f、j屮的各元素。然后运用一种解 线性代数方程的方法，可求的?xi(o)从而求得经第一次迭代后的xi的新值 xi(l)?xi(0)?xi(0)o再将求得的xi代入,乂

25、可得?f、j中的各元素的新值, 从而解得?xi以及xi(2)?xi(l)?xi(l)o如此循环不己，最后可获得此非线性 方程式组足够精确的解。运用这种方法计算时，xi的初值要选择地比较接 近它们的精确解，否则迭代过程可能不收敛。在这系统中，给定一对状态变量us、?s,只要求确定(nl)对状态变量 ui、?io 给6定的?s通常就赋以零值。这实际上就相当于取节点s的电压向量为参 考轴。给定的us 般可取标幺值1.0左右，以使系统中各节点的电压水平 在额定值附近。这样，原则上可从2n个方程式中解出2n个未知量。但是，这个解还 应满足一些约束条件，这些约束条件是保证系统正常运行必不可少的。对控制变量

26、的约束条件是：pgimin?pgi?pgimax； qgimin?qgi?qgimax对无电源的节点，约束条件则为：pgi=o； qgi二0这些pgimin> pgimax> qgimin、qgimax的确定一方面要参照发电机 的运行极限，另一方面还耍计及动力机械所受到的约束。对状态变量ui的约朿条件则是：uimin?ugi?uimax对有些状态变量?i还有如下的约束条件：?i?j?i?jmax这条件主要是保证系统运行的稳态性所要求的。由于扰动变量pli、 qli不可控，对它们没有约束。计算电力系统潮流时，若运用牛顿拉夫逊法可直接用以求解功率方 程。ui?yijuj?pi?jqij

27、?l?(2-10)?e?jf待入式(2-10),并将实数部分和虚数部分分别列出：将yij?gij?jbij,uiiipi?ei(gijej?bijfj)?fi(gijfj?bijej)?j?nqi?(ijfj?bijej?fi(gijej?bijfj)?eig?)j?n(?2bl 1)此外，由于电力系统中还有电压幅值给定的pv节点，还应补充一组 方程注：ei2?fi2?ui2(2?lcl )ei和fi分别为迭代过程中求得的节点电压实部和虚部；7pi为pq节点和pv节点的注入冇功功率；qi为pq节点的注入无功功率；ui为pv节点的电压幅值。牛顿型潮流计算的核心问题是修正方程式的建立和求解。为说明

28、这一 修正方程式的建立过程，先对网络中各类节点的编号作如下约定：(1)网络中共有n个节点，编号为1, 2,，n,其中包含一个平衡节 点，编号为s； (2)网络中有个pq节点，编号为1, 2,，m,包含 编号为s的平衡节点；(3)网络屮有n-m个pv节点，编号为m+1, m+2,， n。由式(2-lla)> (2llb)、(2-llc)所组成的方程式组中共有2(n-l)个独立方 程式。其中，式(2-lla)类型的有(nl)个，包括除平衡节点外所有节点有功 功率pi的表达式，即i=l, 2,，n, ihs；式(2-llb)类型的有(ml)个， 包括所有pq节点无功功率qi的表达式，即曰,2,

29、m %；式(2llc) 类型的有(n-l)-(m-l)=n-m个，包括所有节点pv节点电压ui2的表达式，即 i=m+l, m+2, ?, no平衡节点s的功率和电压之所以不包括在这方程组 内，是由于平衡节点的注入功率不可能事先给定。综上所述：就可以建立类似式(29)的修正方程式如式(2-12)0?p?1?hll?ql?j?ll?p2?h21?q?j?2?21?pp?hpl ?r?2?pl?u?p?p?hnl?n?r2?nl?u?n?n11l11n21l21?nplsplnnlsnlh12j12h22j22?hp2rp2hn2rn2n12l12n22l22?np2sp2nn2sn2hlpjlp

30、h2pj2p?hpprpphnprnpnlpllpn2pl2p?nppsppnnpsnphlnjlnh2nj2n?hpnrpnhnnrn nnln?fl?e?lln?l?n2n?f2?l2n?e2?npn?fp?spn?ep?nnn?fn?snn?en?(2?12)式中的?pi、？qi、？ui2分别为注入功率和节点电压平方的不平衡量。由式(2口)可见，它们分别为?pi?pi?ei(gijej?bijfj)?fi(gijfj?bijej)?(2?13a)?qi?qi?(ijfj?bijej?fi(gijej?bijfj)?eig?)j?n)(?2bl3?ui2?ui2?(ei2?fi2)(2?1

31、3c)8式中雅克比矩阵的各个元素则分别为?p?p?hij?i;nij?i?fj?ej?qi?qi?j?;l?ijij?f?ejj?22?u?uii?rij?;sij?f?ej?j?(2?14)为求取这些偏导数，可将pi、qi、ui2分别展开如下:pi?ei(giiej?biifi)?fi(giifi?biiei)?qi?fi(giiei?biifi)?ei(giifi?biiei)?j?n?ei(gijej?bijfj)?fi(gijfj?bijej)?(2?15a) ?fi(gijej?bijfj)?ei(gijfj?bijej)?(2?15b)j?nei2?fi2?ui2(2?15c)jh

32、i时，由于对特定的j,只有该特定节点的fj和ej是变量，由式(2-14)>(2-15)?p?p?hij=i=-bijei+gijfi; nij=i=gijei+bijfi?fj?ej?qi?q?=-bijfi-gijei=-nij; lij=i=gijfi-bijei=hij?jij=?fj?ej?22?u?uii?rij=o; sij=o?f?e?jj?(2?16a)冃时，为使这些偏导数的表示式更简洁，先引入节点注入电流的表示 式如下?+?yu?l?i 二 yiiuiijjj=lj?ij=n=?giiei-biifi?+?(gijej-bijfj)+j?giifi+biiei?+?(g

33、ijfj+bijej)冋?ij 二 nj 二 n i =au+jbn然后由式(2-14).式(2j5)和上式可得9j=n?pi二biiei+2giifi+biiei+? (gijfj+bijej) =-biiei+giifi+bii?hii=?fj=li?j?i?j=n?pi?nii=?e=2giiei-biifi+biifi+?(gijej-bijfj)=giiei+biifi+aiij=li?j?i?j=n?qi?=-2biifi+giiei-giiei+?(gijej-bijfj)=-giiei-biifi+aii?jii=?fj=li?j?ij=n ?qi?(gijfj+bijej)

34、=-biiei+giifi-bii?lii=?e=giifi-giifi-2biiei-?j=li?j?i?22?r=?ui=2f; s=?ui=2eiiii?ii?f?en ?(2?16b)由式(216a)可见，如yij?gij?jbij?o,即节点i, j之间无直接关系，这些 元素都?hij等于零。从而，如将雅克比矩阵分块，而将每个2x2阶子阵？?jijnij?hij、?lij?rijnij?sij?作分块矩阵的元素时，分块雅克比矩阵和节点导纳矩阵yb将有相同 的结构。但前与后者不同,前者因hij?hji、nij?nji、jij?jji、lij?lji不是对称 矩阵。 2.4潮流计算的基本

35、流程形成了雅克比矩阵并建立了修正方程式，运用牛顿拉夫逊法计算潮流的核心问题已解决，下面列出基木计算步骤并编制流程图。形成节点导纳矩阵yb。设置各节点电压初始值ei(o)> fi(0)o把各节点电压初始值代到式(2-13a)-(2-13c)求修正方程式中的不平衡量?pi(o)、?qi(o)以及?ui(0)2o(4) 把各节点电压初始值代到式(216a)式(216b),求修正方程式的雅克比矩阵中的(0)(0)(0)(0)(0)srhij、nij、jij、l(0)、ijijo ij(0)(0)(0)(5) 把不平衡量?pi(0)、?qi(0)以及?ui(0)2,雅克比矩阵中的hij、nij、j

36、ij、l(0)ij、(0)(0)rij、sij代到式(212),求各节点电压的变化量，即修正量?ei(0)、?fi(0)o(6) 得出各节点电压的新值,ei(l)?ei(0)?ei(0), fi(l)?fi(0)?fi(0)o(刀检查修正量中的最大值，与给定的允许误差作比较，判断是否收敛,如果不收敛，10则以各节点电压新值作为初值自第3步重新进行下一次迭代，否则进 行下一步。(8)计算平衡节点功率和线路功率。其屮，平衡节点功率为s?i?n?s?us?ysiui?ps?qsi?l线路功率为s?ij?uilij?uiuiyio?(ui?uj)yij?pij?jqijs?ji?ujlji?ujujy

37、jo?(uj?ui)yji?pji?jqji线路上的损耗功率为?sij?sij?sji?pij?j?qij(9)输岀结果，完毕。(2?1刀(2?18a)(2?18b) (2?19)112.5牛顿拉夫逊法潮流计算流程图12第三章电力系统无功优化问题及英遗传算法优化求解3.1无功优化问题描述及其模型在给定负荷和无功补偿装置地点的基础上，以有功网损最小为目标函 数，主要考虑了变压器分接头位置、并联电容补偿容量和发电机机端电压 的控制作用。这一类无功优化问题的数学模型如下：目标函数minf?pl?l?(?vi?qi)2?2?()2vimax?viminqimax?qimin(3?l)注：pl为有功网损

38、;?l?(?vi)2为对各节点电压越限的惩罚函数项；vimax?vimin?qi)2为对发电机无功功率越限的惩罚函数项；qimax?qiminvi?vimaxvimin?vi?vimaxvi?vimin?2?(?vi?vimax?vi?0?v?v?imini?qi?qimax?qi?o?q?qi?iminqi?qimaxqimin?qi?qimax qi?qimin约束条件等式约束条件：?pi?vi?v(jgicjo?s?ibj?j?h?sij?bij?qi?vi?vj(gijco?j?h?注：sini?js?iijn)ij)(3?2)pi、qi、vi表示节点i处注入的有功功率、无功功率和节点

39、电压；gij、bij ?ij表示节点i、j之间的电导、电纳和电压之间的相差角；13不等式约束条件：min?vimin?vi?vimax?qgimin?qgi?cigimax?timin?ti?timax?c?c?ciimax?i?vgimin?vgi?vgimax(3?3)注：vimin> vimax表示节点i电压的下、上限；qgmin、qgmax表示发电机输出无功功率的下、上限;timin、timax表示变压器变比的下、上限；cimin、cimax表示并联补偿电容量的下、上限；vgimin> vgimax表示发电机机端电压的下、上限；3.2遗传算法的理论基础遗传算法是建立在达尔文

40、的生物进化论和孟德尔的遗传学说基础上 的算法。牛物体可以通过遗传和变异来适应于外界环境。遗传算法借鉴了 进化原理、遗传原理、随机统计理论。在求解过程中，先形成一个初始群 体，然后一代一代地寻找问题的最优个体，直至满足收敛判据或达到预先设 定的迭代次数才停止。在进化论中认为，每一物种在不断的进化过程中都是越来越适应环境, 物种的每个个体的基本特征被后代所继承，但后代又不完全同于父代这些 新变化，如果适应环境，则被保留下来，否则，将被淘汰。在遗传学说 中认为，每个基因有特殊的位置并控制某个特殊的性质。杂交和突变可能 产生对环境适应性强的后代，进而保存适应值高的基因结构。遗传算法是一种基于自然选择和

41、群体遗传机理的搜索算法，它模拟了 自然选择和白然遗传过程中发生的繁殖、朵交和突变现象。求解问题时， 问题的一个可能解被编码成一个个体，若干个个体构成了群体。在遗传算 法开始时，随机地产生一些个体，在一定约朿条件下，根据目标函数转化 为适应度函数，对每个个体计算出一个适应度值。根据适应度值，选择个 体用来复制下一代。适应度函数值高的较适应度函数值低的有较大的繁殖 能力和机会，选择操作体现了 “适者生存”原理。选择其中相对优秀的个 体进行交叉、变异等遗传操作，产生下一代，重复上述过程，逐步朝着更遗传算法可以看作是一个由可行解组成的群体逐代进化的过程。3.3遗传算法基本原理及操作过程3.3.1适应度

42、函数定标遗传算法在进化搜索中基本上不利用外部信息，仅以适应度函数为依 据，利用群体中每个个体的适应度值来进行搜索，所以选取适应度函数非 常关键。适应度函数是由目标函数转化而成的。对目标函数值域的某种映 射变换称为适应度函数定标。用遗传算法求解优化问题时，要求适应度函 数是无约束的单目标函数。因此将目标模型中的约束条件以惩罚项的形式 加入冃标函数，其它的约束条件则在形成染色体时加以考虑，这样原单冃 标多约束的优化模型就转化为无约束条件的优化模型。由于遗传算法利用 仿真的轮盘赌来寻优，因此只可以求极大值。若优化问题是求极小值的问 题，需对目标函数进行改造变求极小值为求极大值。3.3.2初始解的形成

43、采用遗传算法求解问题时，首先确定问题的变量和目标函数，然后对 变量编码。在遗传算法中，问题的解是用数字串来表示的，而且遗传操作 算子也是直接对数字串进行操作的。编码方式可分为二进制和十进制。具 体使用哪种编码方式，要根据实际的优化问题来确定。二进制编码方式有如下优点：与计算机码制一致，码串的每一位，只有 1和0两个码值，操作简单；表示的范围广，如l位码串最多可表示2l个不 同的变量；适合于表示离散变量。但也存在如下缺点：对于大规模的优化 问题，如果用二进制表示其变量，同时乂要确保解具有一定的精度，数字串位数就会很多，这会使计算量很大，计算用时增多，同时占用了很大的 计算机内存；用二进制来表示变

44、量时，需耍对参数进行编码和译码，用以转 换为十进制，造成了转换误差。采用十进制编码的优点：数字串会比采用二进制表示的数字串短得多, 计算量也会减少，计算用时也会降低；在优化过程中不需对参数进行编码 和译码，也就不存在解的精度问题。遗传算法的计算：?t?s.maxf(x) xmin?x?xmax 15式中x?xlx2?xn, n为控制变量个数。初始化的过程是在所给定的控制变量域xmin,xmax中，随机选取一个 变量群体ixl(0)/x2(0),?/xp(0),其中每个个体x(0),包含n个控制变量，每个控 制变量又可以用j位二进制码串来表示，二进制码串中的每一位码也是随机选取 而来。j取值越大

45、精度越高，但过高的位数将影响到计算的速度，因此应 该同时兼顾两者。每个个体代表了系统的一种状态。p为该群体的规模数。 为了保证在生产下一代的过程中成对选择，p为偶数。3.3.3遗传操作遗传操作是参考生物基因的操作，冃的是根据个体的适应度值对其进 行一定的操作。遗传操作可使解逐代地优化，逼近最优解。遗传操作有三个遗传算子：选择、交叉、变异。1）选择选择建立在适应度评估的基础上。适应度越大的个体，被选择的可能 性就越大。选择岀来的个体放入配对库中作为交叉和变异的前提。轮盘赌方法是遗传算法中最常用的选择方法。设种群的规模为n,个 体的适应度为j?nfj0即概率p反映了个体i的适应度值在整个群体fi,

46、则个体i被 选择的概率pi?fi/j?i?l适应度值总和中所占的比例。个体的适应度值越大，被选中的概率就 越高，被选中的个体被放入配对库中。2）交叉杂交被称为遗传算法中的特色操作，是遗传算法中寻找最优个体的最 主要手段，也是遗传算法区别于其他优化方法的主要标志。交叉是指在配 对库中随机选择两个个体作为双亲，将个体的部分结构加以变换从而生成 新个体。遗传算法的搜索能力通过交叉操作得以飞跃性的提高。交叉操作是按照一定的概率pc,在配对库中随机地选取两个个体进行 的。pc般取0.6到0.9之间。交叉算子有几种分别是：一点交叉：在个体数字串中随机选定一个交叉点，两个个体在该点前 或后进行互换，产生新的

47、个体。例如：16图3.1 点交叉图两点交叉：随机地产生两个交叉点。多点交叉：是两点交叉的推广。一致交叉：通过设置屏蔽字来决定父代个体中相应的基因如何遗传给 子代。当屏蔽字位为1时,父代的两个个体相应位交换生成子代的相应位； 如果屏蔽字位为0,则父代的两个个体的相应位直接复制给子代的相应位。例如：图3.2 一致交叉图3）变异变异就是以很小的概率pm随机地改变个体的某些基因的值。对于交 叉后产生的个体的每一个基因值，产生一个0到1之间的随机数x,若 x<pm,就进行变异操作。在二进制编码中，变异算子随机地将某个基因 值取反,即“0”变成“1”,或“t变为“0”。变异与选择、交叉算了结

48、合在一起，就能避免上两种两算了引起的某 些基因信息的永久性丢失，保证了遗传算法的有效性；同时使得遗传算法 保持群体的多样性，以防止出现未成熟收敛。变异操作是一种防止算法早 熟的措施。在变异操作中，变异率不能取的太大，如果pm&gt；0.5,遗传算法就 退化为随机搜索，而遗传算法的一些重要的数学特性和搜索能力也就不复 存在了。3.4基于遗传算法的无功优化与电压控制实现的步骤（1）确定无功优化数学模型及目标函数采用3.1节所述模型及冃标函数。（2）确定适应度函数因为遗传算法只可求得极大值，而无功优化是求极小值问题，因而需对目标函数进17行改造，变求极小值为求极大值。如下式：max(l/f)

49、所以，适应度函数为：fit?l/?(pl(x)?l?(3) 形成初始解 vi?vilim2q?q)?2?(iilim)2vimax?vilimqimax?qilim输入潮流计算数据，主要包括支路数据(首末节点号、支路电阻、支路 电抗、变压器变比或节点电容等)和节点数据(节点电压幅值及相角、输入 有功无功和负荷有功及无功功率);设置运行参数，如种群规模n,交叉率和变 异率的初始值;进行初始潮流计算,并计算网损值，随机产牛初始化种群；4) 进入计算潮流的循环对每个个体进行潮流计算，将结果代入适应度函数，得到对应的适应 度函数值。循(k)环n次，就可以得到n个适应度值，从中选出最大的适应度函数值 f

50、itmaxo(5) 选择将得到的n个适应度值从大到小排序，用最靠前的1/4代替最后的1/4, 重新形成n个个体。这样做，既选择了比较优良的个体，又保证了种群的 多样性。(6) 交叉进入交叉循环，这里交叉率pc设定为0.9。随机生成一个01之间的 数，如果这个数小于pc,就进行交叉，否则，维持不变。循环n/2次，不 重复交叉。交叉公式为：?xi?(l?a)xi?axj 实型变量交叉公式：? ?x?j?axi?(l?a)xj ?xi?round(l?a)xi?axj整型变量交叉公式：? ?xj?roundaxi?(l?a)xj 注：xi、xj是要交叉的两个体；xi?、x?j是新产生的两个体；a是随机生成的一个0-1之