[发明专利]基于模块度优化的双阶段无功电压分区方法及系统

权利要求书

1.一种基于模块度优化的双阶段无功电压分区方法，其特征在于，包括以下步骤：

步骤1，根据电力系统信息进行中心吸引聚类，本步骤为所述双阶段无功电压分区方法的第一阶段；

所述电力系统信息包括节点信息、支路信息和容量信息；根据电力系统信息以中心吸引聚类方法对电力系统无功电压的节点进行初始无功分区；所述初始无功分区包括邻接节点归并和电源节点归并；所述中心吸引聚类方法包括以下步骤：

步骤1.1，建立复杂网络模型；

根据电力系统的节点信息、支路信息和容量信息建立复杂网络模型，即复杂网络G(V,E)，其中V为根据节点信息建立节点的集合，E为根据支路信息建立边的集合；V中包括中心节点、邻接节点和电源节点；

步骤1.2，建立邻接矩阵；

邻接矩阵A_N×N中的矩阵元素为任意2个节点v_i和v_j的组合；如果v_i和v_j之间存在边，则v_i和v_j所代表的矩阵元素的值a_ij或a_ji为1/(2M)；如果v_i和v_j之间不存在边，则该矩阵元素的值为0；即：

其中，N为复杂网络G(V,E)中的节点的个数，M为复杂网络G(V,E)中边的数量总和；

步骤1.3，计算节点度；

将任意节点v_i与复杂网络G(V,E)内其他节点连接的边的个数设为v_i的节点度k_i；节点v_i的节点度k_i越大，则表示其在复杂网络G(V,E)中的重要性或影响力就越大；邻接矩阵二次幂A_N×N⁽²⁾的对角元素等于v_i的节点度k_i，即：

步骤1.4，选取中心节点；

计算复杂网络G(V,E)中各节点的中心性评价值，选取中心节点；

步骤1.5，进行邻接节点归并；

根据步骤1.2得出的矩阵元素的值a_ij和步骤1.2得出的节点度k_i计算得出任意节点v_i与任意节点v_j的模块度增量ΔQ，计算方法为：

ΔQ＝a_ij+a_ji-2a_ia_j＝2(a_ij-a_ia_j)；

a_i＝k_i/(2M)；

其中，a_i为任意节点v_i对应的辅助向量中元素；

根据ΔQ，对复杂网络G(V,E)进行邻接节点归并，形成各自以中心节点为中心的若干社团结构；

步骤1.6，进行电源节点归并；

对复杂网络G(V,E)中的电源节点进行电源节点归并，从而完成对复杂网络G(V,E)的初始无功分区，即将电源节点归并于各社团结构的中心节点；

步骤2，利用粒子群优化算法对初始种群的模块度Q进行优化，得到模块度Q的最优解和与其相对应最优无功分区；本步骤为所述双阶段无功电压分区方法的第二阶段；

步骤2.1，建立初始种群；

将通过步骤1的中心吸引聚类方法进行的初始无功分区和通过随机聚类方法进行的多个随机无功分区所得出的解集合并为初始种群；

步骤2.2，建立基于电气距离d_ij的权重网络模型D；

根据电力系统建立牛顿-拉夫逊潮流计算方程，即：

其中，J_Pθ，J_Qθ，J_PV和J_QV为雅可比矩阵，ΔP为任意节点v_i注入有功功率的变化量，ΔQ为任意节点v_i注入无功功率的变化量，Δθ为任意节点v_i的节点电压的相角变化，ΔV为任意节点v_i的节点电压的幅值变化；

在实际电力系统运行时，有功功率的变化量对节点电压的影响比较小，因此不考虑有功功率的变化量对节点电压的影响，则令ΔP＝0，将其代入上述牛顿-拉夫逊潮流计算方程，得：

其中，和为雅可比矩阵中无功功率与电压幅值有关的部分；

定义电力系统中任意节点v_i和任意节点v_j之间的电压/电压灵敏度α_ij来表征这两个节点间的紧密程度，即：

对α_ij进行对称性和正定性处理，得到电气距离d_ij，即：

d_ij＝d_ji＝-lg(a_ija_ji)；

其中，电气距离d_ij表征任意节点v_i和任意节点v_j之间的紧密程度；

建立基于电气距离d_ij的权重网络模型D，D的元素取值如下：

步骤2.3，计算初始种群的模块度Q；

将步骤2.1中建立的权重网络模型D作为网络模型中的连边权重，则有A_ij＝D_ij；将初始种群中的各个社团结构进行编号；提出衡量复杂网络中初始种群质量优劣的模块度Q，并计算初始种群中各解集的模块度Q，模块度Q定义为：

其中

A_ij为网络模型中节点i和节点j之间的连边权重；m为整个网络的连边权重；l_i为与节点i连接的连边权重之和；l_j为与节点j连接的连边权重之和；l_il_j/2m为网络内部节点随机相连时，节点i和节点j之间连边权重的期望值；c_i、c_j分别为节点i和节点j所在社团结构的编号；若节点i和节点j在同一社团结构，则δ(c_i,c_j)＝1，否则，δ(c_i,c_j)＝0；

步骤2.4，通过粒子群优化算法进行初始种群的模块度Q的优化求解；

将初始种群中各解集作为粒子群优化算法的粒子；在D维空间中，由N个粒子构成粒子群，则

X_i＝[x_i1,x_i2,...,x_iD](i＝1,2,...,N)，

V_i＝[v_i1,v_i2,…,v_iD](i＝1,2,…,N)；

其中，X_i表示第i个粒子的位置变量；x_id(d＝1,2,…,D)为节点的社团编号；V_i表示第i个粒子的速度变量；v_id(d＝1,2,…,D)为节点的社团的编号修正量；

待优化问题的潜在解由粒子位置v_id表示；粒子的优劣程度由待优化问题的目标函数模块度Q来评价；选出初始种群中模块度Q较大的解集作为初始粒子，各初始粒子根据粒子群算法迭代公式分别更新自身的速度与位置，粒子群算法迭代公式为：

其中，k为迭代次数，k＝1,2,3,…,K，其中K为最大迭代次数；和分别是粒子i在第k次迭代时第d维速度分量和位置分量；rand()为在[0,1]范围内变化的随机函数；c₁和c₂为加速常数，也为学习因子；ω为惯性权重；为异或操作，涉及的两个数相同时取值为0，不同时取值为1；C(rand(n_∈d))表示在邻接节点中随机选取一个，将其所在的社团编号传递给x_id；

迭代过程中，各粒子获取的目标函数值最大的位置P_best记为

P_best＝[p_best1,p_best2,…,p_bestD]；

初始种群的目标函数值最大的位置，即为初始种群的全局最优解G_best记为

G_best＝[g_best1,g_best2,…,g_bestD]；

步骤2.5，经过步骤2.4的优化求解，检测粒子群优化算法是否满足迭代次数，若不满足，则按照步骤2.4中粒子的速度与位置更新方式进行信息更新；若满足，则输出结果为初始种群的模块度Q的最优解，从而得到其相应的最优无功分区；

所述步骤1.4，选取中心节点；包括以下步骤：

步骤1.4.1，计算节点度中心性；

将节点度k_i进行归一化计算，得到节点度中心性C_D(v_i)，节点度中心性C_D(v_i)反应任意节点v_i与其连接的各节点中的中心程度，所述归一化计算方法为：

C_D(v_i)＝k_i/(N-1)；

其中，N为复杂网络G(V,E)中的节点的个数；

步骤1.4.2，计算节点介数和节点介数中心性；

复杂网络G(V,E)中任意不相邻的节点v_j和v_k之间的最短路径途经有可能经过节点v_l，经过节点v_i的最短途径越多，则表示节点v_i在复杂网络G(V,E)中的重要性或影响力就越大；这种重要性或影响力可用节点v_i的节点介数B_i来表示，节点介数B_i定义为：

其中，n_jk为任意不相邻的节点v_j和v_k之间的最短路径的个数；n_jk(i)为v_j和v_k之间的最短路径经过节点v_i的个数；N为复杂网络G(V,E)中的节点的个数；由此可见，节点介数B_i就是复杂网络G(V,E)中所有最短路径中经过节点v_i的数量比例；

将节点介数B_i进行归一化计算，得到节点介数中心性C_B(v_i)，即：

C_B(v_i)＝2B_i/[(N-1)(N-2)]；

步骤1.4.3，计算中心性评价值，选取中心节点；

将节点度中心性C_D(v_i)和节点介数中心性C_B(v_i)进行综合考虑，得出中心性评价值L_l，即：

L_i＝αC_D(v_i)+βC_B(v_i)

其中，α和β为权值，取值范围为[0,1],α+β＝1；

根据L_i的大小进行排序，选取L_i最大的若干节点作为中心节点；选取中心节点的数量与电网中电源的数量相同；将中心节点和电源节点以外的其他节点定义为邻接节点；

所述步骤1.5，所述邻接节点归并的方法为：在v_i为邻接节点、v_j为中心节点时，在v_i与各中心节点的模块度增量ΔQ中，选择最大的ΔQ所对应的中心节点，将v_i归并到此中心节点；

所述步骤1.6，所述电源节点归并的方法为：如果电源节点只与一个中心节点之间存在连接的边，则优先将此电源节点归并于该中心节点；如果电源节点与多个中心节点之间存在连接的边，且各社团结构不满足多负荷水平下的静态无功平衡要求，优先将电源节点归并于无功不足的社团结构的中心节点；如果如果电源节点与多个中心节点之间存在连接的边，且各社团结构皆满足负荷水平下的静态无功平衡要求，则将电源节点归并于与各电源节点短路阻抗距离最小的中心节点；

所述步骤2.4，所述粒子群算法迭代公式中，关于惯性权重ω，现采用一种基于粒子目标函数的自适应调节策略，即每个粒子根据当前目标函数来进行调节；ω值变大，则对应粒子的大范围寻优能力增强；ω值变小，则算法的局部搜索能力增强；因此采用如下搜索方式：

即对于任意粒子而言，当粒子位置远离G_best时，增大ω值；而当接近G_best时，减小ω值；

其中，ω_max为权重系数ω的上限，ω_min为权重系数ω的下限；Q_av为粒子群当前目标函数值的平均值；Q_max为当前G_best的目标函数值；Q_i为粒子当前的适应度值，也就是粒子当前的目标函数值。

2.采用权利要求1所述的基于模块度优化的双阶段无功电压分区方法的系统，包括：电网信息采集模块、初始分区模块、优化分区模块和分区显示模块；

所述电网信息采集模块包括PMU信息采集器和中继传输器；所述PMU信息采集器用于采集电力系统的节点信息、支路信息和容量信息，并通过所述中继传输器传输到所述初始分区模块；

所述初始分区模块用于对电力系统进行初始无功分区；所述初始分区模块包括复杂网络模型单元、中心节点筛选单元、邻接矩阵建立单元、邻接节点归并单元和电源节点归并单元；所述复杂网络模型单元用于建立复杂网络模型；所述邻接矩阵建立单元用于建立邻接矩阵；所述中心节点筛选单元用于计算节点度并选取中心节点；所述邻接节点归并单元用于进行邻接节点归并；所述电源节点归并单元用于进行电源节点归并；

所述优化分区模块用于利用粒子群优化算法得到最优无功分区；所述优化分区模块包括初始种群构建单元、权重网络模型建立单元、模块度计算单元和粒子群优化单元；所述初始种群构建 单元用于建立初始种群；所述权重网络模型建立单元用于建立基于电气距离d_ij的权重网络模型D；所述模块度计算单元用于计算初始种群的模块度Q；所述粒子群优化单元用于通过粒子群优化算法进行初始种群的模块度Q的优化求解，并得出最优无功分区；

所述分区显示模块包括信息接收单元、图像转换单元和分区显示单元；所述信息接收单元用于接收优化分区模块的最优无功分区的相关信息，并通过所述图像转换单元在分区显示单元上显示。