一种配电网节点导纳矩阵的辨识方法技术

本发明专利技术公开了一种配电网节点导纳矩阵的辨识方法，基于总线注入模型，对于不存在隐藏节点的电网系统，从对应于不同稳态的复电压和电流测量的序列中唯一地辨识导纳矩阵；对于存在隐藏节点的星状电网系统，则采用基于图分解的图理论方法以及最大集合的搜索和组合辨识出退化的导纳矩阵；对于存在隐藏节点的网格电网系统，则采用基于低秩和稀疏矩阵分解的方法辨识出退化的导纳矩阵；所辨识出的导纳矩阵可应用于电网故障诊断。

Identification method of node admittance matrix in distribution network

The invention discloses an identification method of distribution network node admittance matrix, bus injection based on power system does not exist hidden nodes, only to identify the sequence admittance matrix from complex voltage and current measurements corresponding to different steady state; for the existence of hidden nodes of the star grid system, using graph theory based on the graph decomposition and the largest collection of search and identify the combination admittance matrix degradation; for grid grid system of hidden nodes, the identification method of low rank and sparse matrix decomposition based on admittance matrix degradation; admittance matrix identifying can be applied to power system fault diagnosis.

【技术实现步骤摘要】
一种配电网节点导纳矩阵的辨识方法
本专利技术属于电力系统安全分析
，更具体地，涉及一种配电网节点导纳矩阵的辨识方法。
技术介绍
近年来，电力行业已经发生了深刻的变化，分布式能源广泛被采用，用户在新兴能源市场的积极参与，量测、通讯及控制基础设施迅速部署的驱动，使配电网在能观性和可控性达到了前所未有的高度。这些变化为系统操作者提供了充分的机会来提高电力系统的效率和稳定性，尽管通过新颖的控制和优化技术增加了不确定性。这种技术需要知道实时的网络拓扑结构，逆潮流问题，就是从许多节点测得的电压电流幅值和相位角来推导节点导纳矩阵的网络拓扑结构。当系统完全可观时，导纳矩阵可由一列不同的稳态量测来唯一地确定，而当系统中包含一些隐藏节点时，可以确定一个退化导纳矩阵；逆潮流问题关注的是节点导纳矩阵的估计，可以从同步电压和电流幅值、相位角量测来获得网络拓扑结构。
技术实现思路
针对现有技术的以上缺陷或改进需求，本专利技术提供了一种配电网节点导纳矩阵的辨识方法，其目的在于基于总线注入模型(BIM)，提供一种高效的辨识导纳矩阵的方法。为实现上述目的，按照本专利技术的一个方面，提供了一种配电网节点导纳矩阵的辨识方法，包括如下步骤：(1)当配电网所有节点都可测时，采用等价的无约束的最小二乘法获取节点导纳矩阵Y；其中，所有节点可测是指所有节点的复电压和电流量均可通过测量得到；(2)当配电网中有隐藏节点时，对于星状网络，则把星状网络分解为一个树和多个团的并集后采用图论的方法获取导纳矩阵；对于网格网络，将通过将导纳矩阵分解为一个稀疏矩阵和低秩矩阵，通过优化问题来获得导纳矩阵。优选地，上述配电网节点导纳矩阵的辨识方法，其步骤(1)包括如下子步骤：(1.1)建立节点注入模型：其中，Ni为与节点i直接相连的节点集；N是指电网网络的节点数，i是指第i个节点的编号，j是指第j个节点的编号，k∈{1,...,K}，K是指第K时刻，H是指矩阵的共轭转置，si(k)是指第i个节点的功率,yij是指第i个节点与第j个节点之间支路的导纳，vi(k)是指第i个节点k时刻的电压,vj(k)是指第j个节点k时刻的电压；(1.2)在所有时间指标下，将上述节点注入模型转换为向量形式；对于节点i：其中，VK是指k时刻的电压，YiN是指第i个节点与其他节点的互导纳，VN(K)是指第N个节点K时刻的电压，IiK,H是指第i个节点的注入电流共轭转置；(1.3)根据上述向量将导纳矩阵辨识问题转换为如下优化问题；min0,s.t:VKY＝IK,H,Y∈SN其中，是指N个节点在指时刻K的注入电流的共轭转置，SN是指对称矩阵；Y是指节点导纳矩阵，VK是指K时刻的电压；(1.4)将上述优化问题转换成约束最小二乘问题，建立目标函数与约束：其中，F是指F范数，s.t:是指subjectto，是指受约束于；(1.5)将向量算子应用到上述约束最小二乘问题的目标函数与约束，建立约束l2问题：其中，是指克罗内克乘积；(1.6)通过将复对称矩阵映射到复向量，将vec(Y)＝Qf(Y)映射到向量化的导纳矩阵，将上述约束l2问题可转化为如下无约束的l2优化问题：其中，Qf(Y)＝vec(Y)，Q是指映射矩阵，M是矩阵乘积的简化表达，f(Y)＝[Y21Y31...YN1Y32Y42...YNN-1]；(1.7)将上述无约束的l2优化问题转换为：其中，获得导纳矩阵，导纳矩阵的实部为Yτ，导纳矩阵的虚部为Yi。优选地，上述步骤(1.6)中，复对称矩阵映射到复向量的方法具体为：令f:为一个从复对称矩阵到复向量的映射，则f(Y)＝[Y21Y31...YN1Y32Y42...YNN-1]，f是任何一个Y∈SN的双射；其中，SN是指对称矩阵。优选地，上述配电网节点导纳矩阵的辨识方法，其步骤(2)中，辨识网格网络的退导纳矩阵的方法，包括如下子步骤：(I)根据下式将退化的导纳矩阵分解为一个稀疏矩阵A和低秩矩阵B：(II)将分解式子作如下变换：(III)将辨识到导纳矩阵的问题转换为下列凸优化问题：其中A是指稀疏矩阵，B是指低秩矩阵；由此，将导纳矩阵Y划分为四个子矩阵，使得Y22∈Ch×h仅对应于隐藏节点的互导纳：其中，C是指复数域，h是指矩阵维数，I1是指可测得的电流，Y11是指非隐藏节点的互导纳，Y21是指隐藏节点和非隐藏节点的互导纳，Y12是指隐藏节点和非隐藏节点的互导纳，Y22是指隐藏节点的互导纳，V1是指将通过计算得到的非隐藏节点的电压，V2是指将通过计算得到的隐藏节点的电压；根据可测量节点的电压和电流时间序列，获得退化的导纳矩阵总体而言，通过本专利技术所构思的以上技术方案与现有技术相比，能够取得下列有益效果：(1)本专利技术提供的配电网节点导纳矩阵的辨识方法，当系统不存在隐藏节点时，从对应于不同稳态的复电压和电流量测序列中唯一确定导纳矩阵；在电力系统运行分析中，往往要计算不同接线方式下的运行状态，当网络接线改变事，节点导纳矩阵也要做相应的修改，本专利技术中所提出的方法实施简便，计算复杂度低，精度高；(2)本专利技术提供的配电网节点导纳矩阵的辨识方法，当电网中存在隐藏节点时，对于星状网络，采用基于图分解的图理论方法以及最大集合的搜索和组合，来辨识出退化导纳矩阵；对于网格网络，基于低秩和稀疏矩阵分解辨识出退化导纳矩阵；根据退化导纳矩阵可得到真实导纳矩阵的部分信息，可进一步用于对电网拓扑的结构变化做出推测。附图说明图1是实施例提供的导纳矩阵辨识方法的流程示意图；图2是实施例所基于的IEEE14母线系统示意图；图3是在没有隐藏节点情况下，采用实施例提供的导纳矩阵辨识方法的故障诊断的结果示意图；图4是在含有隐藏节点情况下，采用实施例提供的导纳矩阵辨识方法的故障诊断的结果示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。电网系统中，变压器、电容、开关的操作和一些故障均可能改变节点间的导纳，甚至可以改变网络的拓扑；导纳矩阵可以反映这些变化；本专利技术提供了一种配电网节点导纳矩阵的辨识方法，包括对无隐藏节点电网和有隐藏节点电网的导纳矩阵辨识；这种方法所需要的数据量小，准确度高，可持续的监测电网系统的变化。根据采集到的电压、电流测量值估计导纳矩阵式的变化；根据导纳矩阵的变化可确定电网事件发生的时间；根据电网系统发生变化前后系统导纳矩阵差的内在稀疏性来定位电网事件所在的位置。实施例中，预设以下条件：给定一个具有一组可量测节点的电力系统，有其中，分别为时刻T前、后的系统退化导纳矩阵；实施例提供的配电网节点导纳矩阵的辨识方法，其流程如图1所示；具体包括如下步骤：(1)为发现导纳矩阵的变化，利用k时刻已知导纳矩阵和电压向量来估计k时刻的注入电流；通过比较k时刻的注入电流估计值和测量值，来计算预测误差其中，I(k)是指k时刻注入电流测量值，是指k时刻注入电流的估计值，V(k)是指k时刻的电压；实施例中，当导纳矩阵没有变化，则序列e(·)为一白噪声序列；当预测误差||e(k)||超过了预设的阈值γ＝10-3，则判定导纳矩阵发生了本文档来自技高网...

【技术保护点】
一种配电网节点导纳矩阵的辨识方法，其特征在于，包括如下步骤：(1)对于所有节点都可测的配电网，采用等价的无约束的最小二乘法获取节点导纳矩阵Y；(2)对应存在有隐藏节点的配电网中，对于星状配电网络，则把星状网络分解为一个树和多个团的并集后采用图论的方法获取导纳矩阵；对于网格配电网络，将通过将导纳矩阵分解为一个稀疏矩阵和低秩矩阵，通过优化问题来获得导纳矩阵。

【技术特征摘要】
1.一种配电网节点导纳矩阵的辨识方法，其特征在于，包括如下步骤：(1)对于所有节点都可测的配电网，采用等价的无约束的最小二乘法获取节点导纳矩阵Y；(2)对应存在有隐藏节点的配电网中，对于星状配电网络，则把星状网络分解为一个树和多个团的并集后采用图论的方法获取导纳矩阵；对于网格配电网络，将通过将导纳矩阵分解为一个稀疏矩阵和低秩矩阵，通过优化问题来获得导纳矩阵。2.如权利要求1所述的辨识方法，其特征在于，所述步骤(1)包括如下子步骤：(1.1)建立节点注入模型：其中，Ni为与节点i直接相连的节点集；N是指电网网络的节点数，i是指第i个节点的编号，j是指第j个节点的编号，k∈{1,...,K}，K是指第K时刻，si(k)是指第i个节点的功率,yij是指第i个节点与第j个节点之间支路的导纳，H是指矩阵的共轭转置，vi(k)是指第i个节点k时刻的电压,vj(k)是指第j个节点k时刻的电压；(1.2)在所有时间指标下，将所述节点注入模型转换为向量形式；对于节点i：其中，VK是指k时刻的电压，YiN是指第i个节点与其他节点的互导纳，VN(K)是指第N个节点K时刻的电压，IiK,H是指第i个节点的注入电流共轭转置；(1.3)根据向量形式将导纳矩阵辨识问题转换为如下优化问题；min0,s.t:VKY＝IK,H,Y∈SN其中，是指N个节点在指时刻K的注入电流的共轭转置，SN是指对称矩阵；Y是指节点导纳矩阵；(1.4)将所述优化问题转换成约束最小二乘问题，建立如下目标函数与约束：其中；F是指F范数；(1.5)将向量算子应用到所述约束最小二乘问题的目标函数与约束中，建立约束l2问题：其中，是指克罗内克乘积；(1....

【专利技术属性】
技术研发人员：袁烨，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：湖北,42