[发明专利]一种基于网络层次分析法的电网节点综合脆弱性指标计算方法

摘要

本发明具体涉及一种基于网络层次分析法的电网节点综合脆弱性指标计算方法，针对节点结构脆弱性和节点状态脆弱性不是独立的，两者是相互联系的，而一般的层次分析法在理论上要求两者相互独立，这在实际应用中会使脆弱性的评估结果不合理，为此本发明基于网络层次分析法，在专家决策方面考虑了节点结构脆弱性和节点状态脆弱性评估结果的相互关系，得到的权重更加精确，另一方面利用熵权法体现客观的权重数值，来进行修正网络层次分析法得到的电网节点结构脆弱性和节点状态脆弱性权重，从而计算得到电网各节点的综合脆弱性评估结果，并进行排序，以此判断电网中存在的薄弱环节。本发明可适应于大电网互联背景下电网节点综合脆弱性评估的需要。

主权项

1.一种基于网络层次分析法的电网节点综合脆弱性指标计算方法，其特征在于，包括步骤1、基于已知电网得到电网的拓扑结构和线路电流、节点电压和注入节点有功功率与无功功率的运行参数，为步骤2中得到节点结构脆弱性指标和状态脆弱性指标提供数据支持；步骤2、基于负荷转移系数和短路容量的方法分别得到电网节点结构脆弱性指标和节点状态脆弱性指标；步骤2.1、采用“电气介数”作为节点结构脆弱性评估模型，以负荷转移因子为基础，获得各节点注入电流变化对每条支路电流的影响，从而得到线路结构脆弱性指标其中node表示电网中所有节点的集合，A_ij,k表示节点k处电流变化对线路ij的影响，P_k、Q_k分别代表节点k当前运行方式下实际输出有功功率和无功功率；进而得到考虑有功功率和无功功率的节点结构脆弱性指标其中，是指与节点i连接的线路l_k的线路脆弱性指标，n表示与节点i连接的线路总数；将节点结构脆弱性指标ω_si进行归一化处理；步骤2.2、采用基于短路容量裕度的节点状态脆弱性指标，对负荷节点进行戴维南等值，通过计算比较负荷节点临界最小短路容量和实际短路容量得到负荷节点的状态脆弱性指标，电源节点和中间节点的状态脆弱性指标为0；并将节点状态脆弱性指标进行归一化处理；步骤3、深入分析影响电网中结构脆弱性和状态脆弱性的因素后，建立相应算例的网络结构分析模型，控制层为电网节点综合脆弱性，网络层为节点状态脆弱性和节点结构脆弱性；步骤4、在网络结构确定后，进行元素间两两比较构造判断矩阵；由于网络层次分析法中被比较的元素之间不是独立的，是有相互联系的，因此比较需要两种方法进行：步骤4.1、元素相互独立时，用直接优势度，即在一个准则下，两元素对该准则的重要程度进行比较；步骤4.2、元素相互联系时，用间接优势度，即两元素在给定的准则下对第三个元素的影响程度进行比较；控制层中有元素综合脆弱性指标，网络层中有元素组C₁,...,C_n，其中C_i有元素e_i1,...,e_in，以控制层元素脆弱性脆弱性指标为准则，以C_j中元素e_jl(l＝1,...,n)为次准则，元素组C_i中元素按其对e_jl的影响力大小进行间接优势度比较，使用1～9标度法构造判断矩阵，并由特征根法得到排序向量1～9标度法、判断矩阵；若元素i和元素j比较得到的评标标度为A_ij，则元素j和元素i相比的结果为1/A_ij；步骤5、一致性检验；判断矩阵的构造过程中要求判断有大体的一致性，以防止出现极端情况；因此在求取排序向量后要对其进行一致性检验，具体步骤如下：步骤5.1、一致性指标C.I.的计算其中，n为判断矩阵的阶数，λ_max为判断矩阵的特征根；步骤5.2、平均随机一致性指标R.I.的计算平均随机一致性指标是指经过多次重复随机判断矩阵特征值的计算后取算数平均数得到的；步骤5.3、一致性比例C.R.的计算当C.R.＜0.1时，则满足一致性检验；当C.R.＞0.1，代表指标权重不满足一致性要求，需要重新调整判断矩阵的元素取值步骤6、求取排序矩阵W_ij：其中列向量为C_i中元素对C_j中元素的影响程度排序向量；当没有影响时，W_ij＝0；最终可综合脆弱性指标控制层下的超矩阵W；步骤7、构造加权超矩阵；上述方法构造的控制层下的超矩阵是非负矩阵，但是W不是列归一化的，所以需要构造加权矩阵，使得超矩阵W为列归一化矩阵；以综合脆弱性指标为准则，对综合脆弱性指标下各组元素对准则C_j的重要性进行比较，得到在控制层下两两判断矩阵，与C_j无关的元素组的排序向量分量为零，由此可得加权矩阵A：其中a_ij∈[0,1]，且则可得加权超矩阵其中的元素为：步骤8、求取极限超矩阵；在网络层次分析法系统中不存在单个元素或最高层次来支配整体，所以我们需要关系两类排序：步骤8.1、相对排序，这是指所有元素对于元素j的影响作用排序；由于在网络层次分析法中元素之间的影响是交互的，比如A影响B，B影响C，C影响A；所有必须找到这种影响的极限情况，也就是极限相对排序，记为LIP；步骤8.2、绝对排序，这是指在给定所有元素在系统中初始重要性排序后，考虑系统元素间累计影响作用所得到的排序；这个排序是对整个系统而言的，不是对某个元素的；此时求取的极限绝对排序，记为LAP；记超矩阵W的k次幂为则一般地，考虑元素i对元素j的累计k步的影响为因此，当W_k在k→∞时如果极限存在，则有同样的，W^kν⁽⁰⁾表示超矩阵对初始绝对排序的累计k步的影响；当W^kν⁽⁰⁾在k→∞时如果极限存在，则有若极限绝对排序ν^∞与ν⁽⁰⁾无关，那么这种系统称为遍历的；如果存在整数k≥1，使得称j由i经k步可达；系统中的子系统C中的任意元素i和不再C中的任意元素j，使得则称C是封闭的；对于不可约矩阵的元素形成最小封闭子集，此时称这样的子集为不可约子集；在系统中存在两个及以上的封闭子集时，系统可进行分解；其中，如果W是非负不可约列随机矩阵，则存在的充分必要条件是W是素阵，素阵是非负不可约矩阵，且满足最大特征根(主特征根)唯一，重数等于1的条件；此时LIP和LAP均存在，且LIP为由LAP作为列组成的矩阵，即极限超矩阵；步骤9、求取脆弱性指标相应的权重；计算上述得到的极限超矩阵，可得极限排序权重，在此基础上可计算某一层次或某一元素的标准化权重；若各元素间极限排序权重为a_i,i＝1,...,n，且某一层次或元素组包含k(k＜n)个元素，则第k个元素的标准化权重ω_k为：最终得到结构脆弱性和状态脆弱性指标对应的的专家权重ξ₁、ξ₂，但是专家权重系数W＝[ξ₁,ξ₂]^T取决于结构脆弱性和状态脆弱性在脆弱性评估体系中的重要程度，需要依据实际算例进行具体分析；网络层次法求取的节点结构脆弱性权重和节点状态脆弱性权重具有较强的主观性，需要用熵权法进行一定的修正，使得权重更加合理准确；熵权法赋权是一种客观赋权方法，在脆弱性的评价中，通过对熵的计算确定权重，就是根据各项评价指标值的差异程度，确定各评价指标的权重；将电网的全部节点作为评价对象，同时建立两个评价指标：结构脆弱性、状态脆弱性；以期根据客观条件的变化来修正综合脆弱性中其相应权重；使用熵权法赋权主要包括以下3个步骤：步骤9.1、原始数据矩阵进行标准化设电网中有m个节点，即原始矩阵中m个评价对象；有n个评价指标得到的原始矩阵为：式中x_ij表示第i个评价节点在第j个脆弱性指标上的指标值；对于某项脆弱性指标j，指标值x_ij的差异越大，则该指标在综合脆弱性评估中所起的作用越大；若某项脆弱性指标的指标值全部相等，则该指标在综合脆弱性评估中几乎不起作用；步骤9.2、原始矩阵的标准化综合脆弱性评估问题可看成一个多目标决策问题，将综合脆弱性的单项指标看作多目标决策中的目标后，综合脆弱性评估就成为多目标决策问题；相对优属度借助模糊数学的概念以类似于隶属度的含义描述相对“优”的程度，并根据目标类型及特点进行确定；将2项脆弱性指标划分为成本型、效益型2种目标类型；1)成本型指标属性值越小越好的指标，其相对优属度为：式中，为指标x_j所有评估点实测值的最大值和最小值；2)效益型指标指属性值越大越好的指标，其相对优属度为：电力系统的结构脆弱性评估，可以科学辨识系统中的关键节点；节点结构脆弱性指标越大，该节点在网架结构中重要性越高，脆弱度越高，属于效益型节点；状态脆弱性反映了系统承受干扰的能力，节点的状态脆弱性越大，说明节点的运行状态越不稳定，脆弱度越高，属于效益型节点；对含有m个评估节点，n个评价指标进行评估，先根据各观测点的脆弱性形成目标评价矩阵X＝(x_ij)_m×n，根据每个指标的属性运用公式(14)、(15)将目标评价矩阵转换为相对优属度矩阵R＝(r_ij)_m×n；步骤9.3、定义熵和熵权在有m个评价节点，n个评价指标的评估问题中，第j个指标的熵定义为：式中，k＝1/lnm，当f_ij＝0时，令f_ijlnf_ij＝0；定义了第j个指标的熵之后，第j个指标的熵权定义为：式中，0≤τ_j≤1，指标的熵权反映该指标在各节点中的差异；熵权越大，表示该指标在各节点中的差异越大，越能为最终的决策提供更多的有用信息，熵权越小，表示该指标在各节点中的差异度越小，当该指标在所有节点中完全一样时，该指标的熵权为0；用熵权法来修正网络层次分析法得到的权重，使权重在体现专家的主观判断的同时还能包含数据本身信息；通过熵权法修正后的综合脆弱性指标权重，即熵权法求取的权重乘以网络层次法的专家权重：式中η_j为第j个指标的综合权重，τ_j为第j个指标下的熵权，ξ_j为由网络层析分析法确定的第j个指标下的权重；将得到的结构脆弱性和状态脆弱性权重分别乘以节点结构脆弱性评估结果和节点状态脆弱性评估结果并进行求和得到电网节点综合脆弱性指标，将上述指标进行归一化后进行排序，指标越大则说明该节点越脆弱，相反的，指标越小则说明该节点越不脆弱；节点综合脆弱性指标可有效全面的评估电网中的薄弱环节，为电网运行和规划人员进行指导。