CN110086168A - 一种基于多源数据融合的配电网运行拓扑辨识方法 - Google Patents

Abstract

本发明公开了一种基于多源数据融合的配电网运行拓扑辨识方法，在配电网安装DSCADA量测装置，在联络开关两端节点和度大于2的多支路节点安装μPMU量测装置；通过μPMU对其安装节点的电压相位进行量测，构建拓扑变化时刻辨识模型，确定拓扑变化时刻；基于拓扑变化前后的节点电压变化，构建拓扑变化趋势判据，确定运行状态发生变化的支路范围，枚举出重构后的所有可能的拓扑情形；将DSACAD和μPMU遥测数据进行融合，构建配电网状态估计模型，对所有可能的拓扑进行状态估计；构建拓扑相似度辨识模型，辨识出实际拓扑，提高了配电网拓扑辨识的可靠性。本发明在单联络开关动作和多联络开关动作情形下具有较好的适用性。

Description

一种基于多源数据融合的配电网运行拓扑辨识方法

技术领域

本发明属于电网运行拓扑技术领域，具体涉及一种基于多源数据融合的配电网运行拓扑辨识方法。

背景技术

随着智能电网建设的加快推进，配电网系统结构日趋复杂，运行拓扑变化更加频繁。配电网通过实时采集遥信数据，及时掌握网络拓扑变化情况，更新运行拓扑结构。然而配电网DSCADA的遥信数据在采集时常存在误报和未报的情形，且数据源较为单一，生成的网络拓扑的正确性具有较大的不确定性，对配电网运行拓扑的辨识可靠性不高。国内电网拓扑辨识研究大多是针对主网进行分析，对结构复杂的配电网拓扑辨识研究较少。而国外对配电网的拓扑辨识均是基于新型μPMU量测装置展开研究，并没有结合国内配电自动化安装情形，目前尚无适用于国内已安装DSCADA的配电网的拓扑辨识研究方法。

发明内容

本发明的针对现有技术中的不足，提供了一种基于多源数据融合的配电网运行拓扑辨识方法，该方法在配电网安装了配电网数据采集及监视控制系统(DSACAD)，并且在部分节点安装少量微型同步相量量测装置(μPMU)情形下，基于多源数据融合实现配电网运行拓扑辨识。

为实现上述目的，本发明采用以下技术方案：一种基于多源数据融合的配电网运行拓扑辨识方法，包括以下步骤：

在配电网安装DSCADA量测装置，在联络开关两端节点和度大于2的多支路节点安装μPMU量测装置；联络开关用来联络两个电源，或者两台变压器；

通过μPMU对其安装节点的电压相位进行量测，构建拓扑变化时刻辨识模型，确定拓扑变化时刻；

基于拓扑变化前后的节点电压变化，构建拓扑变化趋势判据，确定运行状态发生变化的支路范围，枚举出重构后的所有可能的拓扑情形；

将DSACAD和μPMU遥测数据进行融合，构建基于加权最小二乘法的配电网状态估计模型，对所有可能的拓扑进行状态估计；

由节点估计电压相位和μPMU实测电压相位，构建拓扑相似度辨识模型，辨识出实际拓扑。

为优化上述技术方案，采取的具体措施还包括：

所述拓扑变化时刻辨识模型是将配电网某一时刻运行拓扑的辨识问题，转换为确定网络重构时刻，基于原始网络辨识重构后拓扑的问题；

基于μPMU节点电压相位量测构建拓扑变化时刻辨识模型，模型如下：

式中：向量α(t)表示t时刻所有μPMU测量的节点电压相位值所构成的向量，向量的维度等于安装μPMU的节点数；向量α(t-i)表示各安装μPMU的节点在t-i时刻测得的电压相位值构成的向量；x是时间变量；delta(t)为以时间t为变量的1-范数函数；参数w是delta(t)函数的一个时间窗，w越大，t时刻发生网络运行方式的转换时，delta(t)与delta(t-1)值之间的差越大。

将安装μPMU的节点之间除联络开关支路以外的支路定义为一个支路组，将支路组内度为2的节点进行合并，简化配电网络；

同时，拓扑变化需要满足配电网重构的必要条件：

1)为了避免网络出现孤岛运行，在每个支路组中只能打开一个分段开关，即只断开一条支路；

2)当闭合了N个联络开关时，配电网中出现N个环路，只能断开不同环路下各支路组中的N条支路。

所述拓扑变化趋势判据如下：

判据1：当只有一个联络开关动作时，可能发生分段开关动作的支路包含在与该联络开关形成环网的所有支路组中，且只有一个分段开关动作；当有N个联络开关动作时，可能发生分段开关动作的支路包含在形成的N个环网的所有支路组中，且每个环网中各有一个分段开关动作；

判据2：某一支路的分段开关的状态发生变换，那么该支路的两端节点的电压幅值和相位，相比于上一时刻变化较大，推导出判据公式如下：

式中：节点i和j是安装μPMU量测装置的节点；θ_i(t)和θ_j(t)依次为t时刻的节点i和j的电压相位；θ_i(t-1)和θ_j(t-1)依次为t-1时刻的节点i和j的电压相位；U_i(t)和U_j(t)依次为t时刻的节点i和j的电压幅值；U_i(t-1)和U_j(t-1)依次为t-1时刻的节点i和j的电压幅值；U_N为额定电压；ε_θ为节点电压的相位波动约束值；ε_U为节点电压的幅值波动的约束值；

判据3：当某支路组内的某一分段开关状态发生变化时，该支路组两端的节点电压幅值之差的变化最大；

maxM＝||U_j(t)-U_i(t)|-|U_j(t-1)-U_i(t-1)||

式中：M表示t时刻和t-1时刻节点电压差的变化值。

通过拓扑变化趋势判据，确定可能发生开关变位的支路，从而形成可能的重构拓扑，其方法如下：

当网络中仅有1个联络开关动作时，由判据1和判据2初步确定可能动作的支路区间范围；由判据3选择μPMU量测节点间电压差变化最大的两个量测节点，进一步缩小可能存在开关动作的支路区间；

当网络中有N个联络开关动作时，由判据1和判据2确定所有可能动作的支路组后，在形成的每个环网中，由判据3分别选择支路组两端节点电压之差变化最大的支路组，确定为可能发生开关变位的支路区间；

由判据1～3确定可能的分段开关动作的支路范围后，通过枚举法罗列所有可能形成的拓扑情形；当网络闭合N个联络开关时，根据网络重构的必要条件，在N个属于不同环网的支路组内各选取一条支路，通过排列组合方式进行枚举，列出所有可能的运行拓扑。

所述基于加权最小二乘法的状态估计模型，其目标函数的矩阵表达式如下：

J(x)＝[z-h(x)]^TW[z-h(x)]

式中：W为权重矩阵，且W＝[w₁，w₂，…，w_m]^T，其中，w_n为第n个量测量z_n的权重值；J(x)为目标函数矩阵；z为量测向量，h(x)为量测函数向量；

量测误差的协方差对角矩阵R如下：

式中：为第n个量测量z_n的量测误差的方差；E(v)＝0，E(vv^T)＝R；v为量测误差向量，且v＝[v₁，v₂，…，v_m]^T，m为量测量数；

取量测权重矩阵W＝R^-1，那么

从而可以写成：

J(x)＝[z-h(x)]^TR^-1[z-h(x)]

其迭代式可以写成：

式中：k表示迭代次数；和依次表示第k+1次和第k次迭代计算得出的状态量估计值；表示第k次迭代计算中的状态修正量；为雅克比矩阵，由量测函数h(x)对x求导并令x取值得到；

采用收敛判据：

式中：门槛值ε取10^-4～10^-6，n是量测量个数变量。

所述加权最小二乘法的状态估计模型，其相应的量测方程有：

(1)注入功率量测

式中：m表示有网络中有m个节点；θ_ij＝θ_i-θ_j，θ_i和θ_j分别是节点i和节点j的电压相位；G_ij和B_ij为节点导纳矩阵中节点i和j之间元素的实部和虚部；和分别为节点i有功和无功注入量测；P_i和Q_i分别为节点i处的有功和无功量测；和为节点i处的注入功率量测误差，U_i、U_j是节点i和节点j的电压值；

(2)支路首端功率量测

式中：和分别为支路ij首端的有功和无功量测；P_ij是由节点i流向节点j的支路功率；Q_ij是由节点i流向j的无功功率；和分别为对应支路功率量测的量测误差；g_ij和b_ij分别为支路ij的电导和电纳。

(3)支路末端功率量测

式中：和分别为支路ij末端的有功和无功量测；P_ji表示由节点j流向节点i的功率；Q_ji是由节点i流向j的无功功率；

(4)支路电流幅值量测

式中：为支路ij电流幅值的平方；U_i和U_j分别为节点i和节点j的电压值；为电流幅值的量测误差；

根据状态变量和量测方程，求出配电网状态估计模型中对应的雅克比分块矩阵：

式中：矩阵和内的元素分别为节点注入功率以及支路功率对各节点电压的相位和幅值的偏导数；矩阵和内的元素分别为支路电流的幅值的平方对各节点电压的相位和幅值的偏导数。

所述拓扑相似度辨识模型根据内积原理，基于不同拓扑情形下的节点估计电压相位而构建：

式中：J为相似度值；α_i(t)表示在第i种拓扑下的节点电压相位估计向量；α(t)表示节点电压相位实测向量，且两个向量的维度等于μPMU安装数；M为所有可能的拓扑数目。

本发明的有益效果是：本发明在配电网安装DSCADA量测装置的情形下，在联络开关两端节点和度大于2的多支路节点安装少量的μPMU量测装置，通过μPMU对其安装节点的电压相位进行量测，构建拓扑变化时刻辨识模型，确定拓扑发生变化的时刻，进而构建拓扑变化趋势判据，确定运行状态发生变化的支路范围，枚举出重构后的所有可能拓扑情形。将μPMU和DSACAD的遥测数据进行融合，通过多源遥测数据的融合实现配电网运行拓扑辨识，提高了配电网拓扑辨识的可靠性。本发明还考虑了DSCADA和μPMU的量测误差对拓扑辨识准确度的影响，以及该算法在单联络开关动作和多联络开关动作情形下具有较好的适用性。

附图说明

图1是本发明的配电网运行拓扑的辨识方法流程图。

图2是本发明的配电网运行拓扑的辨识方法的具体流程图。

具体实施方式

现在结合附图对本发明作进一步详细的说明。

如附图1、2所示，本发明提供的一种基于多源数据融合的配电网运行拓扑辨识方法，包括以下步骤：

在配电网安装DSCADA量测装置，在联络开关两端节点和度大于2的多支路节点安装μPMU量测装置。

步骤一：通过μPMU量测装置对其安装节点的电压相位进行量测，构建拓扑变化时刻辨识模型，确定拓扑变化时刻t；

拓扑变化时刻辨识模型是将配电网某一时刻运行拓扑的辨识问题，转换为确定网络重构时刻，基于原始网络辨识重构后拓扑的问题，通过此模型确定拓扑发生变化的时刻t。

基于μPMU节点电压相位量测构建拓扑变化时刻辨识模型如下：

式中：向量α(t)表示t时刻所有μPMU测量的节点电压向量的相位值所构成的向量，向量的维度等于安装μPMU的节点数；向量α(t-i)表示各安装μPMU的节点在t-i时刻测得的电压相位值构成的向量；x是时间变量；读取t-1时刻的拓扑结构，然后在t-1时刻拓扑上，经过判别，确定可能发生开断变化的支路，组合形成t时刻可能的拓扑，再进行辨识和筛选；delta(t)为以时间t为变量的1-范数函数；参数w是delta(t)函数的一个时间窗，w越大，t时刻发生网络运行方式的转换时，delta(t)与delta(t-1)值之间的差越大。然而，若w的值设置的较大则会增大程序的计算量，因此，合理的w的值应该满足计算速度要求，且拓扑发生变化时，使得delta(t)值在t时刻产生陡增尖峰。

步骤二：为了简化配电网，需要对网络中的支路先进行分组。将安装μPMU的节点之间的支路(联络开关支路除外)定义为一个支路组，将支路组内度为2的节点进行合并，从而简化配电网络。读取t时刻联络开关处μPMU支路电流量测，确定闭合的联络开关数N。拓扑变化需要满足配电网重构的必要条件：

1)为了避免网络出现孤岛运行，在每个支路组中只能打开一个分段开关，即只断开一条支路；

2)当闭合了N个联络开关时，配电网中出现N个环路，为了构建辐射状运行的连通配电网络，只能断开不同环路下各支路组中的N条支路。

步骤三：确定拓扑发生变化的时刻后，通过μPMU对其安装节点的电压相位进行量测，基于拓扑变化前后的节点电压变化，构建拓扑变化趋势判据，确定运行状态可能发生变化的开关支路范围，枚举出重构后的所有可能的拓扑情形；

其中拓扑变化趋势判据如下：

判据1：当只有一个联络开关动作时，可能发生分段开关动作的支路包含在与该联络开关形成环网的所有支路组中，且只有一个分段开关动作；当有N个联络开关动作时，可能发生分段开关动作的支路包含在形成的N个环网的所有支路组中，且每个环网中各有一个分段开关动作；

由于联络开关两端安装μPMU装置，监测到支路电流时，可以确定联络开关支路为连通状态。根据判据1，可以确定保证配电网辐射状运行的所有可能动作的分段开关的支路组范围。

判据2：某一支路的分段开关的状态发生变换，那么该支路的两端节点的电压幅值和相位(电压量测值由幅值和相位值构成)，相比于上一时刻变化较大，推导出判据公式如下：

式中：节点i和j是安装μPMU量测装置的节点；θ_i(t)和θ_j(t)依次为t时刻的节点i和j的电压相位；θ_i(t-1)和θ_j(t-1)依次为t-1时刻的节点i和j的电压相位；U_i(t)和U_j(t)依次为t时刻的节点i和j的电压幅值；U_i(t-1)和U_j(t-1)依次为t-1时刻的节点i和j的电压幅值；U_N为额定电压；ε_θ为节点电压的相位波动约束值；ε_U为节点电压的幅值波动的约束值；优选地，ε_θ取0.1°，ε_U取4％，由于μPMU量测频率较高，当配电网稳定运行，不发生拓扑变化时，连续两次节点电压相位量测可以看做不变，电压相位波动在量测误差范围内。因此，ε_θ取0.1°；当节点电压相位波动大于约束值，说明该节点的供电路径发生了转换，节点处发生较大的潮流变化。国家标准书中关于电压波动的规定表明：在10kV及以下的配电网中，允许的供电电压波动范围在1.25％-4％之内。因此，ε_U取4％，当分段开关发生开关变位，网络发生重构时，该分段开关附近支路的潮流发生很大变化，附近节点的电压幅值波动较大，不在网络稳定运行时允许的波动范围之内。读取t时刻和t-1时刻下，支路组范围内各支路组两端的μPMU节点的电压幅值和相位量测。

判据3：当某支路组内的某一分段开关状态发生变化时，该支路组两端的节点电压幅值之差的变化最大；

maxM＝||U_j(t)-U_i(t)|-|U_j(t-1)-U_i(t-1)||

式中：M表示t时刻和t-1时刻节点电压差的变化值；U_i(t)和U_j(t)依次为节点i和j在t时刻的电压幅值；U_i(t-1)和U_j(t-1)依次为节点i和j在t-1时刻的电压幅值。

步骤四：通过拓扑变化趋势判据，确定可能发生开关变位的支路，从而形成可能的重构拓扑，筛选出判据1中满足判据2和判据3的支路组，其方法如下：

当网络中仅有1个联络开关动作时，由判据1和判据2初步确定可能动作的支路区间范围；由判据3选择μPMU量测节点间电压差变化最大的两个量测节点，进一步缩小可能存在开关动作的支路区间；

当网络中有N个联络开关动作时，由判据1和判据2确定所有可能动作的支路组后，在形成的每个环网中，由判据3分别选择支路组两端节点电压之差变化最大的支路组，确定为可能发生开关变位的支路区间；

由判据1～3确定可能的分段开关动作的支路范围后，通过枚举法罗列所有可能形成的拓扑情形；当网络闭合N个联络开关时，根据网络重构的必要条件，在N个属于不同环网的支路组内各选取一条支路，通过排列组合方式进行枚举，列出所有可能的运行拓扑。通过拓扑变化趋势判据，可以有效缩小需要列举的可能拓扑组合数目，从而提高辨识的运算速度。

步骤五：采集当前时刻t的遥测数据，基于DSCADA采集的节点注入功率，与μPMU装置采集的支路上的电流和功率信息进行融合，通过两种装置的遥测数据的融合，提高状态估计的冗余度和精度。考虑μPMU和DSCADA量测误差，通过WLS(加权最小二乘法)估计出所有可能的拓扑下的各节点电压相位估计向量。构建基于加权最小二乘法的状态估计模型，对所有可能的拓扑进行状态估计，其目标函数的矩阵表达式如下：

J(x)＝[z-h(x)]^TW[z-h(x)]

式中：W为权重矩阵，且W＝[w₁，w₂，…，w_m]^T，其中，w_n为第n个量测量z_n的权重值；J(x)为目标函数矩阵；z为量测向量，h(x)为量测函数向量；

量测误差的协方差对角矩阵R如下：

式中：为第n个量测量z_n的量测误差的方差；E(v)＝0，E(vv^T)＝R；v为量测误差向量，且v＝[v₁，v₂，…，v_m]^T，m为量测量数。

取量测权重矩阵^W＝R^-1，那么

从而可以写成：

J(x)＝[z-h(x)]^TR^-1[z-h(x)]

其迭代式可以写成：

式中：k表示迭代次数；和依次表示第k+1次和第k次迭代计算得出的状态量估计值；表示第k次迭代计算中的状态修正量；为雅克比矩阵，由量测函数h(x)对x求导并在取值求得。

采用工程中最常用的收敛判据：

式中：门槛值ε取10^-4～10^-6，n是量测量个数变量。

上述的基于加权最小二乘法的状态估计模型，其相应的量测方程有：

(1)注入功率量测

式中：m表示有网络中有m个节点；θ_ij＝θ_i-θ_j，θ_i和θ_j分别是节点i和节点j的电压相位；G_ij和B_ij为节点导纳矩阵中节点i和j之间元素的实部和虚部；和分别为节点i有功和无功注入量测；P_i和Q_i分别为节点i处的有功和无功量测；和为节点i处的注入功率量测误差，U_i、U_j是节点i和节点j的电压值。

(2)支路首端功率量测

式中：和分别为支路ij首端的有功和无功量测；P_ij是由节点i流向节点j的支路功率，由节点i流向节点j；Q_ij是由节点i流向j的无功功率；和分别为对应支路功率量测的量测误差；g_ij和b_ij分别为支路ij的电导和电纳；

(3)支路末端功率量测

式中：和分别为支路ij末端的有功和无功量测；P_ji表示由节点j流向节点i的功率；Q_ji是由节点i流向j的无功功率；

(4)支路电流幅值量测

式中：为支路ij电流幅值的平方；U_i和U_j分别为节点i和节点j的电压值；为电流幅值的量测误差；m表示网络中有m个节点，是下标i和j的取值范围；

根据状态变量和量测方程，求出配电网状态估计模型中对应的雅克比分块矩阵：

式中：矩阵和内的元素分别为节点注入功率以及支路功率对各节点电压的相位和幅值的偏导数；矩阵和内的元素分别为支路电流的幅值的平方对各节点电压的相位和幅值的偏导数。

步骤六：基于μPMU实际量测的节点电压与状态估计出的各可能拓扑下的节点估计电压，构建拓扑相似度辨识模型，由拓扑相似度辨识模型，计算出拓扑相似度J，辨识出实际运行的拓扑。该模型根据内积原理，基于不同拓扑情形下的节点估计电压相位而构建：

式中：J为相似度值；α_i(t)表示在第i种拓扑下的节点电压相位估计向量；α(t)表示节点电压相位实测向量，且两个向量的维度等于μPMU安装数；M为所有可能的拓扑数目。

将α(t)和α_i(t)进行单位化处理后，节点电压相位实测向量会和其中某拓扑的节点电压相位估计向量的夹角较为接近，其J值趋向于1，且和其他拓扑下的节点电压相位估计向量的J值小于1。因此，选取使J值最接近于1的拓扑，确定为实际运行拓扑。

本发明充分利用国内现有的配电自动化装置，在配电网安装DSCADA量测装置的情形下，结合国外新型μPMU，在联络开关两端节点和度大于2的多支路节点安装少量的μPMU量测装置，对提高辨识可靠性具有重要意义。首先，通过μPMU对其安装节点的电压相位进行量测，构建拓扑变化时刻辨识模型，确定拓扑发生变化的时刻。进而构建拓扑变化趋势判据，确定运行状态发生变化的支路范围，枚举出重构后的所有可能拓扑情形。将μPMU和DSACAD的遥测数据进行融合，对所有可能拓扑进行状态估计。由节点估计电压相位和μPMU实测电压相位，构建拓扑相似度辨识模型，辨识出实际运行的拓扑。并考虑DSCADA和μPMU的量测误差对拓扑辨识准确度的影响，以及该算法在单联络开关动作和多联络开关动作情形下的适用性，具有较好的适用性。

以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。

Claims (8)

1.一种基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，包括以下步骤：

在配电网安装DSCADA量测装置，在联络开关两端节点和度大于2的多支路节点安装μPMU量测装置；

通过μPMU对其安装节点的电压相位进行量测，构建拓扑变化时刻辨识模型，确定拓扑变化时刻；

基于拓扑变化前后的节点电压变化，构建拓扑变化趋势判据，确定运行状态发生变化的支路范围，枚举出重构后的所有可能的拓扑情形；

将DSACAD和μPMU遥测数据进行融合，构建基于加权最小二乘法的配电网状态估计模型，对所有可能的拓扑进行状态估计；

由节点估计电压相位和μPMU实测电压相位，构建拓扑相似度辨识模型，辨识出实际拓扑。

2.根据权利要求1所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，所述拓扑变化时刻辨识模型是将配电网某一时刻运行拓扑的辨识问题，转换为确定网络重构时刻，基于原始网络辨识重构后拓扑的问题；

基于μPMU节点电压相位量测构建拓扑变化时刻辨识模型，模型如下：

式中：向量a(t)表示t时刻所有μPMU测量的节点电压相位值所构成的向量，向量的维度等于安装μPMU的节点数；向量α(t-i)表示各安装μPMU的节点在t-i时刻测得的电压相位值构成的向量；x是时间变量；delta(t)为以时间t为变量的1-范数函数；参数w是delta(t)函数的一个时间窗，w越大，t时刻发生网络运行方式的转换时，delta(t)与delta(t-1)值之间的差越大。

3.根据权利要求2所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，将安装μPMU的节点之间除联络开关支路以外的支路定义为一个支路组，将支路组内度为2的节点进行合并，简化配电网络；

同时，拓扑变化需要满足配电网重构的必要条件：

1)为了避免网络出现孤岛运行，在每个支路组中只能打开一个分段开关，即只断开一条支路；

2)当闭合了N个联络开关时，配电网中出现N个环路，只能断开不同环路下各支路组中的N条支路。

4.根据权利要求3所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，所述拓扑变化趋势判据如下：

判据1：当只有一个联络开关动作时，可能发生分段开关动作的支路包含在与该联络开关形成环网的所有支路组中，且只有一个分段开关动作；当有N个联络开关动作时，可能发生分段开关动作的支路包含在形成的N个环网的所有支路组中，且每个环网中各有一个分段开关动作；

判据2：某一支路的分段开关的状态发生变换，那么该支路的两端节点的电压幅值和相位，相比于上一时刻变化较大，推导出判据公式如下：

式中：节点i和j是安装μPMU量测装置的节点；θ_i(t)和θ_j(t)依次为t时刻的节点i和j的电压相位；θ_i(t-1)和θ_j(t-1)依次为t-1时刻的节点i和j的电压相位；U_i(t)和U_j(t)依次为t时刻的节点i和j的电压幅值；U_i(t-1)和U_j(t-1)依次为t-1时刻的节点i和j的电压幅值；U_N为额定电压；ε_θ为节点电压的相位波动约束值；ε_U为节点电压的幅值波动的约束值；

判据3：当某支路组内的某一分段开关状态发生变化时，该支路组两端的节点电压幅值之差的变化最大；

maxM＝||U_j(t)-U_i(t)|-|U_j(t-1)-U_i(t-1)||

式中：M表示t时刻和t-1时刻节点电压差的变化值。

5.根据权利要求4所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，

通过拓扑变化趋势判据，确定可能发生开关变位的支路，从而形成可能的重构拓扑，其方法如下：

当网络中仅有1个联络开关动作时，由判据1和判据2初步确定可能动作的支路区间范围；由判据3选择μPMU量测节点间电压差变化最大的两个量测节点，进一步缩小可能存在开关动作的支路区间；

当网络中有N个联络开关动作时，由判据1和判据2确定所有可能动作的支路组后，在形成的每个环网中，由判据3分别选择支路组两端节点电压之差变化最大的支路组，确定为可能发生开关变位的支路区间；

由判据1～3确定可能的分段开关动作的支路范围后，通过枚举法罗列所有可能形成的拓扑情形；当网络闭合N个联络开关时，根据网络重构的必要条件，在N个属于不同环网的支路组内各选取一条支路，通过排列组合方式进行枚举，列出所有可能的运行拓扑。

6.根据权利要求1所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，所述基于加权最小二乘法的状态估计模型，其目标函数的矩阵表达式如下：

J(x)＝[z-h(x)]^TW[z-h(x)]

式中：W为权重矩阵，且W＝[_w1，_w2，…，w_n]^T，其中，w_n为第n个量测量z_n的权重值；J(x)为目标函数矩阵；z为量测向量，h(x)为量测函数向量；

量测误差的协方差对角矩阵R如下：

式中：为第n个量测量z_n的量测误差的方差；E(v)＝0，E(vv^T)＝R；v为量测误差向量，且v＝[v₁，v₂，…，v_m]^T，m为量测量数；

取量测权重矩阵W＝R^-1，那么

从而可以写成：

J(x)＝[z-h(x)]^TR^-1[z-h(x)]

其迭代式可以写成：

式中：k表示迭代次数；和依次表示第k+1次和第k次迭代计算得出的状态量估计值；表示第k次迭代计算中的状态修正量；为雅克比矩阵，由量测函数h(x)对x求导并令x取值得到；

采用收敛判据：

式中：门槛值ε取10^-4～10^-6，n是量测量个数变量。

7.根据权利要求6所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，所述加权最小二乘法的状态估计模型，其相应的量测方程有：

(1)注入功率量测

式中：m表示有网络中有m个节点；θ_ij＝θ_i-θ_j，θ_i和θ_j分别是节点i和节点j的电压相位；G_ij和B_ij为节点导纳矩阵中节点i和j之间元素的实部和虚部；P_i ^m和分别为节点i有功和无功注入量测；P_i和Q_i分别为节点i处的有功和无功量测；和为节点i处的注入功率量测误差，U_i、U_j是节点i和节点j的电压值；

(2)支路首端功率量测

式中：和分别为支路ij首端的有功和无功量测；P_ij是由节点i流向节点j的支路功率；Q_ij是由节点i流向j的无功功率；和分别为对应支路功率量测的量测误差；g_ij和b_ij分别为支路ij的电导和电纳；

(3)支路末端功率量测

式中：和分别为支路ij末端的有功和无功量测；P_ji表示由节点j流向节点i的功率；Q_ji是由节点i流向j的无功功率；

(4)支路电流幅值量测

式中：为支路ij电流幅值的平方；U_i和U_j分别为节点i和节点j的电压值；为电流幅值的量测误差；

根据状态变量和量测方程，求出配电网状态估计模型中对应的雅克比分块矩阵：

式中：矩阵和内的元素分别为节点注入功率以及支路功率对各节点电压的相位和幅值的偏导数；矩阵和内的元素分别为支路电流的幅值的平方对各节点电压的相位和幅值的偏导数。

8.根据权利要求1所述的基于多源数据融合的配电网运行拓扑辨识方法，其特征在于，所述拓扑相似度辨识模型根据内积原理，基于不同拓扑情形下的节点估计电压相位而构建：

式中：J为相似度值；α_i(t)表示在第i种拓扑下的节点电压相位估计向量；α(t)表示节点电压相位实测向量，且两个向量的维度等于μPMU安装数；M为所有可能的拓扑数目。