一种改进的pmu最优配置整数线性规划方法
【技术领域】
[0001] 本发明属于电力系统状态估计技术领域，特别是涉及一种改进的PMU最优配置整 数线性规划方法。
【背景技术】
[0002] 近年来，我国大容量远距离输电系统的建设和大规模电网互联的实现，一方面使 得电力资源更大范围的优化配置成为可能，另一方面也使得电网复杂度不断增加，从而对 现代电力系统安全稳定问题提出了新的挑战，对电网安全监控技术提出了更高的要求。互 联电网的实时动态检测是安全稳定控制的基础，现有的SCADA/EMS系统只能够提供电网运 行的稳态信息，但已不能满足电力系统在线安全控制的现实要求。而基于GPS的相量测量装 置(PMU)的出现能够实现对电力系统进行同步测量与分析，因而得到了广泛的应用。PMU可 以直接提供在统一时标下的电网暂态过程中的电压、相位角、功率等关键参数的变化曲线， 为电网动态检测、状态估计、在线安全监控等提供方便。
[0003] 建立在PMU基础上的广域测量系统对电网的动态监控提供了有力支持，但由于PMU 的成本巨大以及对通讯设备要求较高等因素，因而目前对每个厂站都配置PMU并不现实。所 以如何实现PMU的最优配置是该领域的研究热点。
[0004] 目前，PMU最优配置问题主要分为两类，即全网可观条件下和网络不完全可观条件 下的PMU最优配置，其主要方法包括整数规划方法、随机算法、启发式算法以及带有启发规 则的随机算法。这些方法的主要缺点是对于零注入节点的处理方式不够严谨;单纯的启发 式算法虽然简单易行，但在网络规模增大时效果并不理想。
[0005] 由于电网节点数目众多，拓扑复杂，有限数量的PMU配置往往不能保证全网可观 测，所以对PMU分阶段安装方案的优化研究尤为必要。部分研究提出了两种常用的不可观测 深度的概念，但省略了形成配置方案的优化模型及诸如初始点选择等一些关键步骤。部分 研究详细、缜密地给出了加入不可观测深度分别为1和2的约束的配置模型，但在网络较大 时较为繁琐，难以实施。
【发明内容】
[0006] 为了解决上述问题，本发明的目的在于提供一种改进的PMU最优配置整数线性规 划方法。
[0007] 为了达到上述目的，本发明提供的改进的PMU最优配置整数线性规划方法包括按 顺序执行的下列步骤：
[0008] 步骤一、判断全网状态是否为完全可观的S01阶段:根据拓扑可观判断方法判断全 网状态是否为完全可观，如果判断结果为"是"，则采用全网可观的PMU配置方法，进入S06阶 段，否则采用全网不完全可观的PMU配置方法，进入S02阶段；
[0009] 步骤二、对是否选择定义1进行判断的S02阶段：如果判断结果为"是"，进入S03阶 段，否则进入S04阶段；
[0010] 步骤三、建立基于定义1的PMU配置模型的S03阶段；
[0011] 步骤四、建立基于定义2的PMU配置模型的S04阶段，然后进入S05阶段；
[0012] 步骤五、执行定义2下的后验算法S5阶段；
[0013] 步骤六、遍历全网各节点，判断其是否有零注入节点的S06阶段:若判断结果为"是 "，则下一步进入S08阶段，否则下一步进入S07阶段；
[0014] 步骤七、执行无零注入节点PMU配置算法的S07阶段；
[0015]步骤八、执行存在零注入节点的PMU配置算法的S08阶段：
[0016] 步骤九、对是否存在零注入节点的可观性进行判断的S09阶段：如果判断结果为 "是"，则接受相应的方案，否则进入S10阶段；
[0017] 步骤十、执行存在零注入节点不可观时的后验算法的S10阶段。
[0018] 在步骤三中，所述的建立基于定义1的PMU配置模型的方法是：
[0019] 定义1:节点i的不可观测深度ri(i)表示其达到最近的可观测节点所需的支路数 目；定义全网的不可观测深度m为单个节点不可观测深度的最大值，即：
[0020] qi=max{ri(i)} ,i = l ,2,··· ,n (16)
[0021] 式中η为系统的节点总数；
[0022] 基于矩阵组合的方法，m = 1时PMU配置模型为：
[0023]
[0024]
[0025] ru = 2时PMU配置模型模型为：
[0026]
[0027]
[0028]以上两式中，T为节点关联矩阵;X= [X1，X2，···，XN]T，其中XiG {〇, 1}，表示第i个节 点是否安装PMU，值为0表示未安装，值为1表示安装;Y = TX，Y中不为0的节点包含了 Xi和与Xi 相连的节点，即直接量测和间接量测所包含的节启
[0029] 在步骤四中，所述的建立基于定义2的PMU配置模型是：
[0030] 定义2:全网的不可观测深度112定义为全网中每个不可观测区域内连通的不可观 节点个数的最大值，BP:
[0031] ri2=max{size(Di)}，i = 1，2，…，m (21)
[0032] 式中，m表示不可观测区域的数量，Di表示第i个不可观测区域；
[0033] ri2=l时，PMU配置模型如下：
[0034]
(22)
[0035](23) Ι'? ~
^
[0036] 式中，C的每一行为全网中任意两个相连节点的组合，若节点i，j相连，则矩阵中必 存在某一行的第i和第j个元素值为1，其他元素值为〇;1^是维数与C的行数相同的全1向量。 [0037] 当ri2 = 2时，PMU配置模型如下：
[0038] (2.4)
[0039] (25)
[0040] 在步骤五中，所述的执行定义2下的后验算法的方法是：
[0041] 1)采用分支定界法求解S04阶段获得的模型，得到若干PMU配置方案，对每个方案 判断是否拓扑可观，若可观，则该方案可作为PMU配置的可选方案;若不可观，则说明存在T 形接线节点不可观，剔除对应配置方案后进行下一步；
[0042] 2)判断剩余方案数是否为零，若否，则输出这些配置方案;若是，则增加1个PMU作 为约束条件后转步骤1);
[0043] 后验算法执行完成后即得到了定义2下的最优PMU配置方案。
[0044] 在步骤七中，所述的执行无零注入节点PMU配置算法的方法是：
[0045] 建立无零注入节点PMU配置模型为：
[0046]
[0047]
[0048] 由于^为0或1，故这是一个整数优化问题，采用分支定界法或割平面法等针对整 数规划的算法进行求解;若存在多组最优解，则在最优解中寻找使系统冗余度最大的一组 解作为最终的PMU配置方案。
[0049] 在步骤八中，所述的执行存在零注入节点的PMU配置算法的方法是：
[0050] 建立存在零注入节点的PMU配置模型为：
[0051]
[0052]
[0053] 其中T，X，b的定义与前文相同，G为ηΧη矩阵，对角线元素均为1，若存在含有k个节 点的零注入区域，所有k个节点构成集合D，则gij = 1 /k，i，j e D，i矣j，其他元素均为0;
[0054] 采用分支定界法或割平面法求解该模型，得到初始的最优配置方案。
[0055] 在步骤十中，所述的执行存在零注入节点不可观时的后验算法的方法是：
[0056] 如果方案存在类似步骤九所提的不可观情况，则加入约束式(45)后重新求解存在 零注入节点的PMU配置模型，直到得到合格的方案为止，如果合格的最优方案数t>l，则选择 其中冗余度较大的方案作为最优配置方案；
[0057]
i=l
[0058] 本发明提供的改进的PMU最优配置整数线性规划方法针对全网可观和不完全可观 两种情况提出了新的PMU配置方案，使得在全网可观时能够简单、准确地处理存在零注入或 注入已知节点的情况，得到使PMU数量最少、冗余度最大的最优解;在全网不完全可观的条 件下，能够针对不同情况，在较短时间内得到满足不可观测深度要求的PMU数量最少的最优 解，为解决可观测系统和不可观测系统的PMU优化配置问题提供了新的解决方案。由于方案 直观简单，能够准确、快速得到最优解，且针对PMU配置各个阶段的多种情况，故易于在实际 操作中应用和推广。
【附图说明】
[0059] 图1为本发明提供的改进的PMU最优配置整数线性规划方法流程图。
[0060] 图2为零注入节点区域示意图。
【具体实施方式】
[0061] 下面结合附图和具体实施例对本发明提供的改进的PMU最优配置整数线性规划方 法进行详细说明。
[0062] 如图1所示，本发明提供的改进的PMU最优配置整数线性规划方法包括按顺序执行 的下列步骤：
[0063] 步骤一、判断全网状态是否为完全可观的S01阶段:根据拓扑可观判断方法判断全 网状态是否为完全可观，如果判断结果为"是"，则采用全网可观的PMU配置方法，进入S06阶 段，否则采用全网不完全可观的PMU配置方法，进入S02阶段；
[0064]步骤二、对是否选择定义1进行判断的S02阶段：如果判断结果为"是"，进入S03阶 段，否则进入S04阶段；
[0065] 步骤三、建立基于定义1的PMU配置模型的S03阶段：
[0066] 定义1:节点i的不可观测深度n(i)表示其达到最近的可观测节点所需的支路数 目；定义全网的不可观测深度m为单个节点不可观测深度的最大值，即：
[0067] qi=max{ri(i)} ,i = l ,2,··· ,n (31)
[0068] 式中η为系统的节点总数；
[0069] 基于矩阵组合的方法，ru = 1时PMU配置模型为：
[0072] ru = 2时PMU配置模型模型为：
[0070]
[0071]
[0073]
[0074]
[0075]以上两式中，T为节点关联矩阵;Χ= [χι，Χ2，···，xn]T，其中XiG {〇, 1}，表示第i个节 点是否安装PMU，值为0表示未安装，值为1表示安装;Y = TX，Y中不为0的节点包含了 Xi和与Xi 相连的节点，即直接量测和间接量测所包含的节点；6 = [U. ;
[0076]由于以上两种模型中的^为整数，所以仍可以利用整数线性规划的方法，如分支 定界法求解。值得注意的是，矩阵组合算法在全