一种考虑零注入约束的双线性WLAV状态估计方法技术

本发明专利技术公开了一种考虑零注入约束的双线性WLAV状态估计方法，其特征包括：1获得电力系统的网络参数；2得到节点导纳矩阵；3从电力系统中获取遥测数据；4计算各量测量的权重；5阵计算雅克比矩阵；6用等式约束来描述零注入约束计及到双线性WLAV状态估计。本发明专利技术有效改善零注入节点的估计精度，从而能得到更精确的状态估计结果。

【技术实现步骤摘要】
一种考虑零注入约束的双线性WLAV状态估计方法
本专利技术涉及一种考虑零注入约束的双线性WLAV状态估计方法，属于电力系统状态估计领域。
技术介绍
电力系统状态估计是能量管理系统中的核心基础模块,准确的状态估计结果能为现代电网的实时分析、控制提供保障。实际电力系统中存在许多零注入节点，理论上零注入节点的注入功率严格为零，然而实际状态估计中这些节点的注入功率往往不为0，导致这些节点上的状态估计结果不能满足潮流方程，影响后续的安全评估结果，因此如何保证状态估计结果更好地符合零注入约束值得更深入的研究。目前状态估计中处理零注入约束的方法主要有大权重法和拉格朗日乘子法。大权重法给精确的零注入功率量测赋以大权重参与状态估计，形式简洁，也不用对现有程序加以修改，但这种方法会严重劣化信息矩阵条件数，严重时甚至造成程序的不收敛。拉格朗日乘子法可极大地提高求解的稳定性和精确性，但由于实际电力系统中零注入节点的比例很高，拉格朗日乘子法将导致信息矩阵规模显著增大，影响状态估计的效率。加权最小绝对值(WeightedLeastAbsoluteValue,WLAV)是一种研究广泛的抗差状态估计方法，计算精度高，但非线性内点法的计算效率难以满足工程需求。变量代换方法，将状态估计问题分解为两步线性化过程，提出了双线性WLAV抗差估计，使运算效率得到显著提升。但由于该方法没有考虑零注入约束，应用于实际系统时将产生较大的零注入节点功率失配量，影响电力系统的安全评估。
技术实现思路
本专利技术所要解决的技术问题是针对双线性WLAV状态估计难以满足对零注入节点功率的平衡要求这一缺陷，提供一种考虑零注入约束的双线性WLAV状态估计方法，以期能有效改善零注入节点的估计精度，从而能得到更精确的状态估计结果。本专利技术为实现上述目的，采用如下技术方案：本专利技术一种考虑零注入约束的双线性WLAV电力系统状态估计方法的特点是按如下步骤进行：步骤1、获得电力系统的网络参数，包括：输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗；步骤2、初始化，包括：对状态量x设置初值，所述状态量x包括节点电压的幅值V和相角θ；根据所述电力系统的网络参数得到节点导纳矩阵；步骤3、从所述电力系统中获取遥测数据z，包括：电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率；所述遥测数据z即为电力系统状态估计的量测量；步骤4、根据量测仪表的精度，计算各量测量的权重；步骤5、取中间变量y，令：式(1)中，n为系统的节点个数；Ｖi，Ｖj分别为节点i、节点j的电压幅值，θij为连接节点i和节点j的支路两端相角差，θij＝θi-θj，θi、θj分别为节点i、节点j的电压相角；步骤6、根据所述节点导纳矩阵计算雅克比矩阵C；步骤7、对所述雅克比矩阵C进行排序，建立如式(2)所示的排序后的方程：式(2)中，zB和zN分别表示零注入量测和非零注入量测，yB和yN分别表示将中间变量y的元素进行排序后的与零注入节点有直接电气联系的部分以及与零注入节点没有直接电气联系的部分，CBB、CBN、CNB和CNN分别表示排序后的雅克比矩阵C的子矩阵；步骤8、令式(2)中的零注入量测zB为零，则得到式(3)：步骤9、求解式(4)得到中间变量y，并将所得y中元素的顺序恢复到式(1)中的顺序，式(4)中，w为权重矩阵，是对角元素为各量测量权重的对角矩阵，l，u为松弛变量，且l,u≥0；步骤10、对中间变量y进行非线性变换，获得变换后的中间变量y′，如式(5)所示，步骤11、根据所述变换后的中间变量y′，利用内点法求得中间变量x′＝[lnVθ]T；步骤12、利用式(6)得到所述状态量x，从而实现电力系统状态估计：与现有技术相比，本专利技术的有益效果体现在：考虑到零注入节点的功率失配量对电力系统安全评估有着重要影响，本专利技术基于零注入节点注入功率为0的特性，用等式约束来描述零注入约束，计及到双线性WLAV状态估计中，该方法既能保留双线性WLAV状态估计原有的计算效率高的优点，与传统不考虑零注入约束的双线性WLAV相比，又能减小零注入节点的功率失配量，从而提高了状态估计结果的精度。附图说明图1为本专利技术的流程图；图2为本专利技术采用的线路π形等值电路图。具体实施方式下面结合附图对专利技术的技术方案进行详细说明：本实施例中，如图1所示，一种考虑零注入约束的双线性WLAV电力系统状态估计方法是按如下步骤进行：步骤1、获得电力系统的网络参数，包括：输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗；步骤2、初始化，包括：对状态量x设置初值，状态量x包括节点电压的幅值V和相角θ；根据电力系统的网络参数得到节点导纳矩阵；步骤3、从电力系统中获取遥测数据z，包括：电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率；遥测数据z即为电力系统状态估计的量测量；电力系统状态估计属于非线性优化问题，其量测方程可表示为：z＝h(x)+ε(1)式(1)中，h(x)为非线性函数；ε为量测误差向量。步骤4、根据量测仪表的精度，计算各量测量的权重；步骤5、取中间变量y，令：式(1)中，n为系统的节点个数；Vi，Vj分别为节点i、节点j的电压幅值，θij为连接节点i和节点j的支路两端相角差，θij＝θi-θj，θi、θj分别为节点i、节点j的电压相角；步骤6、根据节点导纳矩阵计算雅克比矩阵C；结合电力系统的特性，通过引入中间变量可将非线性方程式(1)分解为：z＝Cy+ε(3)y′＝f(y)(4)y′＝Dx′+εy′(5)x＝t(x′)(6)式(3)-式(6)中，y，y′，x′为中间变量；f(y)和t(x′)为非线性函数；C和D为常系数矩阵；εy′为中间变量y′的误差向量。状态估计量测量包括节点电压量测、节点功率量测、支路有功和无功量测，在极坐标系下，其计算方程为：式(7)中，Pi，Qi分别为节点i的注入有功及无功功率；Pij，Qij分别为支路ij上的有功和无功功率；Gij，Bij分别为节点i和节点j之间的互电导和互电纳；Gii，Bii分别为节点i的自电导和自电纳；gij，bij，yc分别为图2所示的支路ij的π形等效电路的串联电导、串联电纳和对地电纳。由上式可得雅克比矩阵C的各元素值如表1所示。表1雅克比矩阵C步骤7、对雅克比矩阵C进行排序，建立如式(8)所示的排序后的方程：式(8)中，zB和zN分别表示零注入量测和非零注入量测，yB和yN分别表示将中间变量y的元素进行排序后的与零注入节点有直接电气联系的部分以及与零注入节点没有直接电气联系的部分，CBB、CBN、CNB和CNN分别表示排序后的雅克比矩阵C的子矩阵；步骤8、令式(8)中的零注入量测zB为零，则得到式(9)：步骤9、求解式(10)得到中间变量y，并将所得y中元素的顺序恢复到式(2)中的顺序，式(10)中，w为权重矩阵，是对角元素为各量测量权重的对角矩阵，l，u为松弛变量，且l,u≥0；构造出式(10)的拉格朗日函数：式(11)中，α，β，η为拉格朗日乘子。由KKT(Karush-本文档来自技高网...

【技术保护点】
1.一种考虑零注入约束的双线性WLAV电力系统状态估计方法，其特征是按如下步骤进行：步骤1、获得电力系统的网络参数，包括：输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗；步骤2、初始化，包括：对状态量x设置初值，所述状态量x包括节点电压的幅值V和相角θ；根据所述电力系统的网络参数得到节点导纳矩阵；步骤3、从所述电力系统中获取遥测数据z，包括：电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率；所述遥测数据z即为电力系统状态估计的量测量；步骤4、根据量测仪表的精度，计算各量测量的权重；步骤5、取中间变量y，令：

【技术特征摘要】
1.一种考虑零注入约束的双线性WLAV电力系统状态估计方法，其特征是按如下步骤进行：步骤1、获得电力系统的网络参数，包括：输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗；步骤2、初始化，包括：对状态量x设置初值，所述状态量x包括节点电压的幅值V和相角θ；根据所述电力系统的网络参数得到节点导纳矩阵；步骤3、从所述电力系统中获取遥测数据z，包括：电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率；所述遥测数据z即为电力系统状态估计的量测量；步骤4、根据量测仪表的精度，计算各量测量的权重；步骤5、取中间变量y，令：式(1)中，n为系统的节点个数；Vi，Vj分别为节点i、节点j的电压幅值，θij为连接节点i和节点j的支路两端相角差，θij＝θi-θj，θi、θj分别为节点i、节点j的电压相角；步骤...

【专利技术属性】
技术研发人员：杨友情，黄启建，吴常胜，李进，赵葆雯，徐健，钱瑛，李朝东，邓晓艳，王为臻，
申请(专利权)人：国家电网公司，国网安徽省电力公司池州供电公司，
类型：发明
国别省市：北京,11