基于改进粒子群算法的配电网重构优化方法技术

本发明专利技术公开了一种基于改进粒子群算法的配电网重构优化方法，包括：配电网络简化并进行粒子群编码；粒子群解码，计算每个粒子对应的适应度函数；随机初始化粒子群，对每一粒子随机赋予一个可行解；对每个粒子所对应的简化网络进行拓扑检查，直至全部粒子符合配电网拓扑要求；对每个粒子通过算法迭代进行位置更新，并进行简化网络拓扑检查；计算每个粒子更新位置后的适应度值：如果粒子i的适应值优于此前个体极值位置Pi对应的适应值，则将其位置赋值给Pi；如果当前极值位置Pi的适应值优于此前全局极值位置Pg的适应值，则将其位置赋值给Pg，直到达到预设最大迭代次数，Pg即为最优解。本发明专利技术方法具有快速、高效的全局寻优能力。

【技术实现步骤摘要】
基于改进粒子群算法的配电网重构优化方法
本专利技术涉及一种基于改进粒子群算法的配电网重构优化方法。
技术介绍
配电网系统具有闭环系统，开环运行的特点。配电网重构是优化配电系统运行的有效手段，它主要通过切换联络开关和分段开关的开合状态来改变网络拓扑结构，以改变网络中的潮流流动。在正常的运行条件下，根据运行情况进行开关操作以调整网络结构，一方面能够平衡负荷，消除过载，提高供电电压质量；另一方面，能够提高系统运行的经济性。配电网重构是在保证配电网呈现辐射状结构、满足馈线热容、节点压降和变压器容量等条件的前提下，使配电网某一指标或者多重指标最佳的配电网运行方式。由于配电网结构复杂，存在着大量的分段开关和联络开关，因此配电网重构是一个复杂的多目标高维数非线性混合优化问题。目前的求解算法大致可分为以下几类：1)传统数学优化算法：即直接利用现有的数学优化原理进行配电网重构。随维数的增加，该方法将面临“组合爆炸”问题，难以实际应用；2)启发式方法，主要有最优流模式法和支路交换法，由于结合了配电网重构问题的物理特性，计算速度有了很大提高，但重构结果与打开开关的顺序或网络中开关的初始状态有关，每一次优化计算只搜索了整个解空间的一部分子空间，缺乏数学意义上的全局最优性；3)人工智能算法，如模拟退火法、遗传算法、禁忌算法等，这些方法中的某些参数随着网络规模等因素的变化而变化，难以确定，并且存在收敛速度慢，容易出现大量不可行解，不成熟收敛等问题。
技术实现思路
本专利技术的目的在于克服现有技术中的不足，提供一种基于改进粒子群算法的配电网重构优化方法，能够提高重构优化的全局寻优能力，解决配电网重构中出现大量不可行解的技术问题。为达到上述目的，本专利技术所采用的技术方案是：基于改进粒子群算法的配电网重构优化方法，包括如下步骤：步骤一：配电网络简化并进行粒子群编码，具体步骤如下：S1-1：编码不考虑不在环路上的支路；S1-2：解环效果一样的支路合并成一个支路组，称之为广义支路；S1-3：用一个粒子维度代表网络中的一条广义支路，粒子的大小代表这条广义支路对应的原始网络中的支路编号；步骤二：粒子群解码，并通过配电网潮流计算求得每个粒子对应的适应度函数fitness；步骤三：设定改进粒子群算法参数，包括粒子群规模为M、学习因子c1和c2、最大迭代次数k；步骤四：随机初始化粒子群，对于种群中的每一粒子位置Xθ在搜索空间中随机赋予一个可行解，其中θ＝1，2…，M；步骤五：对每个粒子所对应的简化网络进行拓扑检查，直至全部粒子符合配电网拓扑要求，采用基于广度优先树搜索的配电网拓扑分析；步骤六：计算每个粒子的适应度值，设定每个粒子的当前位置为其个体极值位置pθ，取全局极值位置Pg为当前群体最优适应度值对应的位置；步骤七：对每个粒子通过算法迭代进行位置更新，并进行简化网络拓扑检查；所述的算法迭代中，第k+1次迭代时，粒子位置Xθ的更新公式为：其中，为点对点乘法；L(β)为莱维飞行的随机搜索路径其中，L(β)计算方法如下：式中：s为莱维飞行跳跃步长；参数β取值范围为0<β<2，参数μ、ν为正态分布随机数，服从下式的正态分布：式中标准差σμ、σv的取值为：步骤八：计算每个粒子更新位置后的适应度值：如果粒子θ的适应度值优于此前个体极值位置pθ对应的适应度值，则将其位置赋值给pθ；如果当前极值位置pθ的适应度值优于此前全局极值位置Pg的适应度值，则将其位置赋值给Pg；步骤九：进入下一次迭代过程，直到达到预设最大迭代次数k；步骤十：粒子解码并输出结果：Pg即为最优解，对其解码可得重构优化后的网络开关状态，最优目标函数值为fitness(Pg)。步骤五的具体操作步骤如下：S5-1：初始化节点信息列表Node，分别存放节点编号、属层和遍历；初始化支路信息列表Branch，分别存放支路首端节点编号、末端节点编号、遍历和开关状态；S5-2：将Node和Branch列表中的遍历置0、Node中的属层置1；S5-3：确定根节点，从根节点开始，由Branch列表逐个搜索与之相连的节点，即开关状态为1的节点，将搜索到的节点属层加1，将搜索过的节点和支路的遍历加1；S5-4：从属层为2的节点开始，逐个搜索与之相连的节点，将搜索到的节点属层加1，将搜索过的节点和支路的遍历加1；S5-5：重复第S5-4步，直到所有闭合支路遍历均为1时，完成了整个网络的拓扑分析；S5-6：检查所有节点遍历是否均为1，当所有节点遍历均为1时，表示满足辐射状拓扑结构要求；否则，表示不满足配电网拓扑要求。适应度函数fitness的公式如下：式中，ηnetloss为网损率；Pnetloss为网络损耗；Pload_i为负荷节点i的有功功率；N为全网负荷节点数；Rij为第ij支路电阻；Vij为第ij支路末端的电压幅值；kij表示该支路是否断开，如果断开则为0，闭合为1；Pij和Qij分别为流过第ij支路的有功和无功功率；i、j分别表示第ij支路的两个负荷节点。与现有技术相比，本专利技术所达到的有益效果是：采用改进粒子群算法进行配电网重构优化，具有快速、高效的全局寻优能力；对原始配电网络进行简化，基于广义之路的粒子编码方案，减小了粒子维度，同时能避免大部分不可行解得产生，通过简化网络的拓扑检查，完全避免了不可行解，提高了重构优化的计算速度；将莱维飞行矢量引入粒子群搜索，改善了粒子群算法易陷入局部最优的问题，保留了其高效性和鲁棒性，提高了粒子群算法的全局收敛能力；与传统基于遗传算法的配电网重构优化相比，采用粒子群算法能同时实现对离散变量和连续变量的编码，因此可扩展本专利技术，用于解决含分布式电源配电网的综合优化问题。附图说明图1为本专利技术的流程图。图2为本专利技术中改进粒子群算法流程图。图3为本专利技术中简化网络拓扑检查流程图。图4为PG&E69节点配电网络。图5为图4简化后网络。图6为传统粒子群算法和改进粒子群算法的重构优化收敛曲线对比图。具体实施方式下面结合附图对本专利技术作进一步描述。步骤一：配电网络简化并进行粒子群编码，具体方法如下：S1-1：编码不考虑不在环路上的支路；S1-2：解环效果一样的支路合并成一个支路组，称之为广义支路；S1-3：用一个粒子维度代表网络中的一条广义支路，粒子的大小代表这条广义支路对应的原始网络中的支路编号。步骤二：粒子群解码，并通过配电网潮流计算求得每个粒子对应的适应度函数fitness；步骤三：设定改进粒子群算法参数，包括粒子群规模为M、学习因子c1和c2、最大迭代次数k；步骤四：随机初始化粒子群，对于种群中的每一粒子位置Xθ在搜索空间中随机赋予一个可行解，其中θ＝1，2…，M；步骤五：对每个粒子所对应的简化网络进行拓扑检查，直至全部粒子符合配电网拓扑要求，采用基于广度优先树搜索的配电网拓扑分析，具体操作步骤如下：S5-1：初始化节点信息列表Node，分别存放节点编号、属层和遍历；初始化支路信息列表Branch，分别存放支路首端节点编号、末端节点编号、遍历和开关状态；S5-2：将Node和Branch列表中的遍历置0、Node中的属层置1；S5-3：确定根节点，从根节点开始，由Branch列表逐个搜索与之相连的节点，即开关状态为1的节点，将搜索到的节点属层加1，将搜索过的节点和支路本文档来自技高网...

【技术保护点】
基于改进粒子群算法的配电网重构优化方法，其特征在于：包括如下步骤：步骤一：配电网络简化并进行粒子群编码；步骤二：粒子群解码，并通过配电网潮流计算求得每个粒子对应的适应度函数fitness；步骤三：设定改进粒子群算法参数，包括粒子群规模为N、学习因子c1和c2、最大迭代次数k步骤四：随机初始化粒子群，对于种群中的每一粒子Xi在搜索空间中随机赋予一个可行解，其中i＝1，2…，N；步骤五：对每个粒子所对应的简化网络进行拓扑检查，直至全部粒子符合配电网拓扑要求，采用基于广度优先树搜索的配电网拓扑分析；步骤六：计算每个粒子的适应度值，设定每个粒子的当前位置为其个体极值位置Pi，取全局极值位置Pg为当前群体最优适应值对应的位置；步骤七：对每个粒子通过算法迭代进行位置更新，并进行简化网络拓扑检查；步骤八：计算每个粒子更新位置后的适应度值：如果粒子i的适应值优于此前个体极值位置Pi对应的适应值，则将其位置赋值给Pi；如果当前极值位置Pi的适应值优于此前全局极值位置Pg的适应值，则将其位置赋值给Pg；步骤八：进入下一次迭代过程，直到达到预设最大迭代次数k；步骤九：粒子解码并输出结果：Pg即为最优解，对其解码可得重构优化后的网络开关状态，最优目标函数值为fitness(Pg)。...

【技术特征摘要】
1.基于改进粒子群算法的配电网重构优化方法，其特征在于：包括如下步骤：步骤一：配电网络简化并进行粒子群编码，具体步骤如下：S1-1：编码不考虑不在环路上的支路；S1-2：解环效果一样的支路合并成一个支路组，称之为广义支路；S1-3：用一个粒子维度代表网络中的一条广义支路，粒子的大小代表这条广义支路对应的原始网络中的支路编号；步骤二：粒子群解码，并通过配电网潮流计算求得每个粒子对应的适应度函数fitness；步骤三：设定改进粒子群算法参数，包括粒子群规模为M、学习因子c1和c2、最大迭代次数k；步骤四：随机初始化粒子群，对于种群中的每一粒子位置Xθ在搜索空间中随机赋予一个可行解，其中θ＝1，2…，M；步骤五：对每个粒子所对应的简化网络进行拓扑检查，直至全部粒子符合配电网拓扑要求，采用基于广度优先树搜索的配电网拓扑分析；步骤六：计算每个粒子的适应度值，设定每个粒子的当前位置为其个体极值位置pθ，取全局极值位置Pg为当前群体最优适应度值对应的位置；步骤七：对每个粒子通过算法迭代进行位置更新，并进行简化网络拓扑检查；所述的算法迭代中，第k+1次迭代时，粒子位置Xθ的更新公式为：其中，为点对点乘法；L(β)为莱维飞行的随机搜索路径其中，L(β)计算方法如下：式中：s为莱维飞行跳跃步长；参数β取值范围为0<β<2，参数μ、ν为正态分布随机数，服从下式的正态分布：式中标准差σμ、σv的取值为：步骤八：计...

【专利技术属性】
技术研发人员：何小栋，张强，左智勇，
申请(专利权)人：南京工程学院，
类型：发明
国别省市：江苏;32