[发明专利]一种基于HACS算法求解电力系统无功优化的方法

权利要求书

1.一种基于HACS算法求解电力系统无功优化的方法，其特征在于，包括以下步骤：

1)、建立电力系统无功优化数学模型，即确定目标函数和约束条件，其中目标函数选择有功功率损耗最小作为目标函数，系统的约束条件包括等式约束和不等式约束两类，并获取电力系统的运行参数；

2)、种群初始化，初始种群通过一个矩阵表示，每一列代表一个鸟巢，鸟巢的位置即代表无功优化问题的潜在解，在约束条件范围内，在d维空间内随机产生n个鸟巢；

3)、确定包含罚函数的目标函数，通过潮流计算求得初始种群的适应度值；

4)、执行布谷鸟算法中的莱维飞行机制更新鸟巢位置；

5)、引入包含全局收敛引导变异算子的人工蜂群算法中雇佣蜂和跟随蜂的开采行为，更新全局最优值；

6)、通过布谷鸟算法中的发现概率P_a淘汰适应度值较差的个体；

7)、判断迭代次数是否满足kk_max，若满足，则返回到步骤4)重复执行步骤4)-步骤7)，否则循环结束，输出无功优化最优解；

所述步骤5)引入包含全局收敛引导变异算子的人工蜂群算法ABC中雇佣蜂和跟随蜂的开采行为，更新全局最优值，具体包括：

①全局收敛引导的变异算子

全局收敛引导的变异算子在进化初期采用均匀的搜索方式，在后期采用局部化的搜索方式，引导个体向最优个体飞行，其定义方式如下：

式中，为第i个鸟巢位置的第j维分量；为当前搜索中得到的最优鸟巢的第j维分量；γ随机取0或1；函数Δ(k,y)为[0,y]之间的一个值：

式中rand为[0,1]之间的随机数；k_max为最大迭代次数；

②变异时机的选择

通过评价最优鸟巢适应度值的变化率的方式确定变异时机，判定方式为：

式中，为第k次迭代时最优鸟巢的适应度值，Δ为阈值；

其思想为:本次迭代所求得的最佳个体适应度值与之前第t次迭代所取得的适应度值的差值的绝对值，再与本次迭代所得适应度值之比值，若其比值小于限定的某个阈值，则引入ABC算法中雇佣蜂和跟随蜂的开采行为进行变异；

③ABC算法中雇佣蜂和跟随蜂的开采行为根据开采方程式更新全局最优值：

式中，j∈[1,2,…,d],m＝[1,2,…,n]，n为种群规模，且m≠i；φ_i,j是[-1,1]之间的随机数，左边项为全局收敛引导变异结束后的第i个位置的第j个元素；等式右边第一个因子为全局收敛引导变异结束后未进行变异前的第i个位置的第j个元素；为全局收敛引导变异结束后未进行变异前的随机的第m个位置。

2.根据权利要求1所述的基于HACS算法求解电力系统无功优化的方法，其特征在于，所述步骤1)考虑电力系统的经济性，选择经典的有功功率损耗最小作为目标函数：

式中，P_loss为电网中的有功功率损耗(MW)；N_L为系统的支路数；G_k为节点i、j之间的第k条支路的电导；V_i、V_j分别为节点i和j的电压幅值；δ_i、δ_j分别为节点i和j的电压相角。

3.根据权利要求1所述的基于HACS算法求解电力系统无功优化的方法，其特征在于，所述等式约束来源于潮流计算中有功和无功功率的平衡：

式中，NB、N和N_PQ分别表示与节点i相连的节点数、总的节点数和PQ节点数；其中总的节点是平衡节点以外的其他节点的总和，P_Gi和P_Di分别表示第i个PV节点和第i个PQ节点的有功功率；Q_Gi和Q_Di分别表示第i个PV节点和第i个PQ节点的无功功率；G_ij和B_ij分别表示Y_ij的实部和虚部，其中Y_ij为节点导纳矩阵Y中第i行第j列的元素；δ_ij＝δ_i-δ_j为节点i和j之间的电压相角差。