一种基于量子遗传算法的10kV配网无功优化方法
【专利摘要】本发明提供了一种基于量子遗传算法的10kV配网无功优化方法，首先构建有功网损最小的模型，对所述控制变量进行量子比特概率幅编码来构建模型解空间的个体，而不是直接采用控制变量数值本身进行编码，由于概率幅含义为全部可能状态的等概率相加，因此采用概率幅进行量子编码的长度，能够小于传统二进制编码的长度，从而提高算法的计算效率。并且本方法提出新的量子交叉的方法，进行量子旋转门更新的两个体的更新目标互相，能够改善种群多样性，避免个体收敛于自身的局部最优解的情况。因此本发明能够缩短运算周期，提高计算效率，有机会跳出局部最优解，进入新的搜索空间，使得算法最终收敛于全局最优解。
【专利说明】-种基于量子遗传算法的1〇kV配网无功优化方法
【技术领域】
[0001] 本发明涉及电力【技术领域】，尤其涉及一种基于量子遗传算法的l〇kV配网无功优 化方法。
【背景技术】
[0002] 配电网是直接向用电客户供电的电力网络，10KV配电线路的电压质量显得尤为重 要。作为影响电压的重要因素之一，解决好配电网络无功补偿的问题，对电网的安全性和降 损节能有着重要意义。
[0003] 配网无功补偿主要通过调节有载调压分接头和无功补偿容量投切控制器来实现， 目前，无功优化方法主要是智能方法，智能方法的代表之一是遗传算法但遗传算法的计算 效率低、优化速度慢、但参数选取不当时，会过早收敛于局部最优解。
[0004] 因此现在需要新型的遗传算法以缩短运算周期提高计算效率，有机会跳出局部最 优解，进入新的搜索空间，使得算法最终收敛于全局最优解。
【发明内容】
[0005] 本发明提供了一种基于量子遗传算法的10kV配网无功优化方法，本发明能够缩 短运算周期，提高计算效率，有机会跳出局部最优解，进入新的搜索空间，使得算法最终收 敛于全局最优解。
[0006] 为了实现上述目的本发明提供了以下技术特征：
[0007] -种基于量子遗传算法的10kV配网无功优化方法，包括：
[0008] 构建电能有功网损最小的模型，模型的约束条件为系统潮流约束，潮流约束中的 变量约束包括控制变量的约束条件和状态变量的约束条件；
[0009] 基于量子遗传算法对所述控制变量进行量子比特概率幅编码，构建模型的个体， 控制变量的参数不同，个体也随之不同；
[0010] 随机产生初始种群，将初始种群作为当前种群；
[0011] 对当前种群中的每个个体进行测量，获得每个个体的测量值，对每个个体进行适 应度计算获得适应值；获得种群中最优个体及最优适应值、局部最优个体及局部适应值； [0012] 判断是否达到终止进化条件；
[0013] 当达到终止进化条件则输出最优个体及最优适应值；
[0014] 当未达到终止进化条件，则以随机抽取当前种群中两个个体，及与两个个体对应 的局部最优解，将两个个体的局部最优解互换，利用量子交叉方法产生新个体，将新个体加 入当前种群获得新种群，将新种群作为当前种群，再对当前种群中的每个个体进行测量，直 到达到终止进化条件。
[0015] 优选的，基于量子遗传算法对所述控制变量进行量子比特概率幅编码，构建模型 的个体包括：
[0016] 采用k个量子比特来表示分布式电源发电机机端电压，k取决于分布式电源发电 机机端电压的量化精度；
[0017] 采用1个量子比特的有载可调变压器变比分接头档位，1取决于接头调控范围的 最大值；
[0018] 采用m个量子比特的无功补充装置投切组数，m取决于投切组数调控范围的最大 值；
[0019] 将k个量子比特的分布式电源发电机机端电压、1个量子比特的有载可调变压器 变比分接头档位和m个量子比特的无功补充装置投切组数，组成一个控制变量的个体。
[0020] 优选的，对每个个体进行适应度计算获得适应值包括：
[0021] 对每个个体进行测量，得到二进制编码，获得二进制个体；
[0022] 在所述约束条件下，将所述二进制个体解码为十进制个体；
[0023] 对种群每个十进制个体利用潮流计算方法进行计算，获得所有变量的值；
[0024] 利用变量值代入适应度函数，获得适应值。
[0025] 优选的，将两个个体的局部最优解互换，利用量子交叉方法产生新个体包括：
[0026] 将其中一个个体作为第一个体，另外一个作为第二个体，与第一个体对应的局部 最优解为第一最优目标，与第二个个体对应的局部最优解为第二最优目标；
[0027] 将第一最优目标作为第二个体的局部最优目标，对第二个体进行一次量子旋转门 更新；
[0028] 将第二最优目标作为第一个体的局部最优目标，对第一个体进行一次量子旋转门 更新；
[0029] 还原各自的局部最优目标，生成新个体。
[0030] 优选的，终止进化条件包括：
[0031] 当前种群达到预设代数、所有个体适应值均一致或预设数量代数的适应值没有得 到改善。
[0032] 优选的，还包括：
[0033] 当达到终止进化条件后，还输出所有变量的值。
[0034] 本发明提供了一种基于量子遗传算法的10kV配网无功优化方法，首先构建有功 网损最小的模型，对所述控制变量进行量子比特概率幅编码来构建模型解空间的个体，而 不是直接采用控制变量数值本身进行编码，由于概率幅含义为全部可能状态的等概率相 力口，因此采用概率幅进行量子编码的长度，能够小于传统二进制编码的长度，从而提高算法 的计算效率。
[0035] 并且本方法提出新的量子交叉的方法，进行量子旋转门更新的两个体的更新目标 互相，能够改善种群多样性，避免个体收敛于自身的局部最优解的情况。因此本发明能够缩 短运算周期，提高计算效率，有机会跳出局部最优解，进入新的搜索空间，使得算法最终收 敛于全局最优解。
【专利附图】
【附图说明】
[0036] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据 提供的附图获得其他的附图。
[0037] 图1为本发明实施例公开的一种基于量子遗传算法的10kV配网无功优化方法的 流程图；
[0038] 图2为本发明实施例公开的又一种基于量子遗传算法的10kV配网无功优化方法 的流程图；
[0039] 图3为本发明实施例公开的又一种基于量子遗传算法的10kV配网无功优化方法 的流程图；
[0040] 图4为本发明实施例公开的又一种基于量子遗传算法的10kV配网无功优化方法 的流程图。
【具体实施方式】
[0041] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0042] 在详细介绍本发明之前首先说明一下量子遗传算法，该算法是以量子计算的概念 和理论为基础，用量子位编码来表示个体，用量子旋转门更新完成进化搜索，具有种群规模 小、同时兼顾种群多样性和改善进化方向的能力、收敛速度快。本发明将量子遗传算法应用 于10kV配网无功优化，针对该问题，提出相应的编码方式、量子交叉方式，来提高算法的效 率和全局最优解的概率，下面介绍本发明的【具体实施方式】 ：
[0043] 如图1所示，本发明提供了一种基于量子遗传算法的10kV配网无功优化方法，包 括：
[0044] 步骤S101 :构建电能有功网损最小的模型，模型的约束条件为系统潮流约束，潮 流约束中的变量约束包括控制变量的约束条件和状态变量的约束条件；
[0045] 控制变量是随机产生的，约束条件是电网运行或者设备预先设定的，潮流约束，即 有功平衡，无功平衡，及相关变量的限制约束。
[0046] 步骤S102 :基于量子遗传算法对所述控制变量，进行量子比特概率幅编码构建模 型的个体，控制变量的参数不同，个体也随之不同；
[0047] 步骤S103 :随机产生初始种群，将初始种群作为当前种群。
[0048] 步骤S104:对当前种群中的每个个体进行测量，获得每个个体的测量值，对每个 个体进行适应度计算获得适应值；获得种群中最优个体及最优适应值、局部最优个体； [0049] 在求解适应值的过程中可以获得种群中的最优个体和大于最优个体对应的最优 适应值，依据多个局部最优个体及与其对应的局部适应值，将满足所有约束条件，或者与约 束条件相比出现越界行为，越界幅值最小的者为最优个体。
[0050] 步骤S105 :判断是否达到终止进化条件；
[0051] 本发明中，可以选取三种条件作为终止进化条件：
[0052] 第一种：当前种群达到预设代数；初始种群每更新一代，则增加一次代数，当当前 种群达到预设规定的代数后，则终止进化，理论上认为，更新预设代数后，便可出现最优化 个体，所以达到预设代数后，便停止进化，预设代数为人为依据现有技术经验得出的。
[0053] 第二种：所有个体适应值均一致；
[0054] 当所有的个体的适应值均一致，则表明此时各个个体之间的都是类似的，无最优 个体之分，所以此时便无需再进行更新和优化。
[0055] 第三种：预设数量代数的适应值没有得到改善，第三种为即便进行预设数量的更 新后，个体的适应值没有得到改善，还是不能够得到一个最优个体，此时则终止进化条件， 因为即便再更新，也不会有更优个体出现了，所以便停止进化。
[0056] 步骤S106 :当达到终止进化条件则输出最优个体及最优适应值；当达到终止进化 条件后，还会输出所有变量的值。即进行潮流计算的后获得的结果。
[0057] 步骤S107 :当未达到终止进化条件，则随机抽取当前种群中两个个体，及与两个 个体对应的局部最优解，将两个个体的局部最优解互换，利用量子交叉方法产生新个体，将 新个体加入当前种群获得新种群，将新种群作为当前种群，再进入步骤S104。
[0058] 当未达到终止条件，则要继续更新和进化当前种群，为了避免产生局部最优解，本 发明采用量子交叉算法对初始种群进行更新。优选的，将两个个体对应的局部最优个体进 行交换，产生新种群。
[0059] 本发明通过量子个体编码方式缩短传统二进制编码的长度，提高算法的计算效 率；将量子理论引入到遗传算法中改善种群的多样性，针对同一优化问题所采用的种群数 小于基本遗传算法；提出新的量子交叉的方法，将交叉概念融合到量子旋转门更新中，能够 解决10kV配网的无功优化方案中计算效率低、有可能收敛到局部最优解的问题，还能够得 到10kV配网无功优化配置及变压器分接头最佳位置。
[0060] 优选的，上述步骤S101中基于量子遗传算法对所述控制变量进行量子比特概率 幅编码构建控制变量的个体，如图2所示，包括：
[0061] 步骤S201 :采用k个量子比特来表示配网中分布式电源发电机机端电压，k取决于 分布式电源发电机机端电压的量化精度；
[0062] 将分布式电源发电机机端电压数值用量子比特编码，采用Uq表示，电压实数值可 转化为k位二进制码，每一位二进制码用一位量子比特来表示，k取决于量化精度。例如：k 值为2,其量化精度为l/2k = 0. 25 ;每个量子比特同时用概率幅α，β来表示，概率幅用来 表示个体的可能状态的概率。
[0063] 步骤S202 :采用1个量子比特的有载可调变压器变比分接头档位，1取决于接头调 控范围的最大值；
[0064] 可调分压器变比分接头的档位采用Tq表示，有多个调控范围，用来调节分压器， 使变压器工作在正常范围内，例如：调控范围有7档，则21〉= 7, 1 = 3, 一个量子比特代表 一位二进制码。
[0065] 步骤S203:采用m个量子比特的无功补充装置投切组数，m取决于投切组数调控范 围的最大值；
[0066] 无功补充装置投切组数采用Nq表示，编码方法同分接头档位量子编码。
[0067] 步骤S204 :将k个量子比特的分布式电源发电机机端电压Uq、1个量子比特的有 载可调变压器变比分接头档位Tq和m个量子比特的无功补充装置投切组数Nq，组成一个控 制变量的个体。
[0068] 将分布式电源发电机机端电压Uq、分接头档位Tq和投切组数Nq组合形成一个控 制变量个体，即[Uq Tq Nq]。
[0069] 本发明采用量子个体编码方式，相当于传统二进制编码的方式而言，缩短了长度， 方便了处理器进行运算。
[0070] 优选的，步骤S103中对每个个体进行适应度计算获得适应值，如图3所示，包括：
[0071] 步骤S301 :对每个个体进行测量得到二进制编码，获得二进制个体；
[0072] 步骤S302 :在所述约束条件下，将所述二进制个体解码为十进制个体；
[0073] 步骤S303 :对种群每个十进制个体利用潮流计算方法进行计算，获得所有变量的 值；
[0074] 步骤S304 :利用变量值构建适应度函数，求解适应度函数获得适应值。
[0075] 如图4所示，优选的，将两个个体的局部最优解互换，利用量子交叉方法产生新个 体包括：
[0076] 步骤S401 :将其中一个个体作为第一个体，另外一个作为第二个体，与第一个体 对应的局部最优解为第一最优目标，与第二个个体对应的局部最优解为第二最优目标；
[0077] 步骤S402 :将第一最优目标作为第二个体的局部最优目标，对第二个体进行一次 量子旋转门更新；
[0078] 步骤S403 :将第二最优目标作为第一个体的局部最优目标，对第一个体进行一次 量子旋转门更新；
[0079] 步骤S404 :还原各自的局部最优目标，生成新个体。
[0080] 量子交叉优化算法，即将第一最优个体，作为第二个体的最优目标，将第二最优个 体作为第一个体的最优目标，按新的目标进行更新，所以第一个体便不会再局限于自己的 第一最优个体，第二个体也不会局限于自己的第二最优个体，每次更新的时候，都会产生新 的个体进行更换，所以不会局限于局部最优解，而且还使种群具有多样性，解决了现有技术 中局部最优解的问题。
[0081] 下面介绍本发明的一种具体实施例，包括如下的步骤：
[0082] 1、采集电网运行参数，运行参数包括：电网拓扑、线路参数、变电站参数、各种运行 方式下的各节点负荷参数、发电机参数、节点电压允许运行范围、无功补偿设备参数、所有 控制变量的约束条件、状态变量的约束条件；
[0083] 2、构造由系统控制变量组成的个体，随机产生初始种群。
[0084] a、量子个体编码
[0085] 本发明提出一种新的量子个体编码，将各控制变量用概率幅来表示。
[0086]
【权利要求】
1. 一种基于量子遗传算法的lOkV配网无功优化方法，其特征在于，包括： 构建电能有功网损最小的模型，模型的约束条件为系统潮流约束，潮流约束中的变量 约束包括控制变量的约束条件和状态变量的约束条件； 基于量子遗传算法对所述控制变量进行量子比特概率幅编码，构建模型的个体，控制 变量的参数不同，个体也随之不同； 随机产生初始种群，将初始种群作为当前种群； 对当前种群中的每个个体进行测量，获得每个个体的测量值，对每个个体进行适应度 计算获得适应值；获得种群中最优个体及最优适应值、局部最优个体及局部适应值； 判断是否达到终止进化条件； 当达到终止进化条件则输出最优个体及最优适应值； 当未达到终止进化条件，则以随机抽取当前种群中两个个体，及与两个个体对应的局 部最优解，将两个个体的局部最优解互换，利用量子交叉方法产生新个体，将新个体加入当 前种群获得新种群，将新种群作为当前种群，再对当前种群中的每个个体进行测量，直到达 到终止进化条件。
2. 如权利要求1所述的方法，其特征在于，基于量子遗传算法对所述控制变量进行量 子比特概率幅编码，构建模型的个体包括： 采用k个量子比特来表示分布式电源发电机机端电压，k取决于分布式电源发电机机 端电压的量化精度； 采用1个量子比特的有载可调变压器变比分接头档位，1取决于接头调控范围的最大 值； 采用m个量子比特的无功补充装置投切组数，m取决于投切组数调控范围的最大值； 将k个量子比特的分布式电源发电机机端电压、1个量子比特的有载可调变压器变比 分接头档位和m个量子比特的无功补充装置投切组数，组成一个控制变量的个体。
3. 如权利要求1所述的方法，其特征在于，对每个个体进行适应度计算获得适应值包 括： 对每个个体进行测量，得到二进制编码，获得二进制个体； 在所述约束条件下，将所述二进制个体解码为十进制个体； 对种群每个十进制个体利用潮流计算方法进行计算，获得所有变量的值； 利用变量值代入适应度函数，获得适应值。
4. 如权利要求1所述的方法，其特征在于，将两个个体的局部最优解互换，利用量子交 叉方法产生新个体包括： 将其中一个个体作为第一个体，另外一个作为第二个体，与第一个体对应的局部最优 解为第一最优目标，与第二个个体对应的局部最优解为第二最优目标； 将第一最优目标作为第二个体的局部最优目标，对第二个体进行一次量子旋转门更 新； 将第二最优目标作为第一个体的局部最优目标，对第一个体进行一次量子旋转门更 新； 还原各自的局部最优目标，生成新个体。
5. 如权利要求1所述的方法，其特征在于，终止进化条件包括： 当前种群达到预设代数、所有个体适应值均一致或预设数量代数的适应值没有得到改 善。
6.如权利要求3所述的方法，其特征在于，还包括： 当达到终止进化条件后，还输出所有变量的值。
【文档编号】G06Q50/06GK104091214SQ201410347101
【公开日】2014年10月8日 申请日期:2014年7月21日 优先权日:2014年7月21日
【发明者】许家益, 胡振斌, 邵名声, 程金松, 汪宏华, 李敏, 吴哲, 朱兵 申请人:国家电网公司, 国网安徽省电力公司黄山供电公司