一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法技术方案

本发明专利技术公开了一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法，该方法包括(1)首先获取用户智能电表电压时间序列及配变TTU电压时间序列数据，对获取的数据进行降噪、特征提取等预处理；(2)数据预处理后，利用马尔可夫随机场对配电系统网络结构建模，建立描述配电网中各节点之间相关关系的联合概率分布；(3)通过最大似然估计法求解联合概率分布，得到配电网节点间的相关性矩阵，根据节点相关性矩阵，建立配电系统网络结构的邻接矩阵，完成配电系统网络结构的识别，本发明专利技术利用历史电表电压数据等已知数据对低压配电网络结构展开识别，不需要新增设备，减少了投资，具有较高的识别准确率。

A network topology identification method of low voltage distribution system based on Markov random field

【技术实现步骤摘要】
一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法
本专利技术涉及电力系统配电
，尤其涉及一种基于用户智能电表历史数据，利用马尔可夫随机场对低压配电网拓扑进行识别的方法。
技术介绍
目前，国内外的中高压配电系统已经具备了完备的配电管理系统，通过数据采集与监视控制系统和地理信息系统获取配电网拓扑结构、运行状态、负荷状态等一系列数据，实现状态评估；并以此建立各类配电网分析和决策算法，包括处理配电网中负荷和可再生能源不确定性的配电网随机优化调度方法、以故障恢复或者三相不平衡治理为目标的配电网重构方法等。相比之下，低压配电系统面临有限信息环境，缺乏有效的系统建模与状态评估方法，难以展开配电系统潮流计算以及其他配电管理高级功能。尤其是，低压配电系统的拓扑结构建模，是建立低压配电系统最优潮流的基础，也是治理低压配电系统三相不平衡、网络损耗、消纳可再生能源、接纳电动汽车充电负荷、提升用户用电可靠性等等一系列问题的基础。因此，识别配电网拓扑是实现低压配电系统可视化和管理的必备条件。考虑到低压配电系统中能够获取到的运行数据非常有限，如何有效利用智能电表数据以及配变运行数据，实现配电系统拓扑识别，建立配电系统基本运行模型，是当前发展低压配电系统的首要目的。
技术实现思路
本部分的目的在于概述本专利技术的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和专利技术名称中可能会做些简化或省略以避免使本部分、说明书摘要和专利技术名称的目的模糊，而这种简化或省略不能用于限制本专利技术的范围。鉴于上述现有基于马尔可夫随机场的低压配电系统网络拓扑识别方法存在的问题，提出了本专利技术。因此，本专利技术目的是提供一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法，其利用历史电表电压数据等已知数据对低压配电网络结构展开识别，不需要新增设备，减少了投资，具有较高的识别准确率。为解决上述技术问题，本专利技术提供如下技术方案：此种基于马尔可夫随机场的低压配电系统网络拓扑识别方法包括如下步骤：(1)获取用户智能电表电压时间序列数据及配变TTU电压时间序列数据；(2)对步骤(1)中获取的电压数据进行一阶差分、分箱处理、独热编码处理；(3)利用马尔可夫随机场对配电系统网络结构建模，建立描述配电系统网络结构中各节点之间相关关系的联合概率分布；(4)通过最大似然估计法求解联合概率分布，得到配电系统网络结构节点间的相关性矩阵；(5)根据节点相关性矩阵，建立配电系统网络结构的邻接矩阵。作为本专利技术所述基于马尔可夫随机场的低压配电系统网络拓扑识别方法的一种优选方案，其中：所述步骤(1)中，数据采集的要求样本数在500个时间截面以上。作为本专利技术所述基于马尔可夫随机场的低压配电系统网络拓扑识别方法的一种优选方案，其中：步骤(2)的数据预处理具体包括：步骤2.1：对步骤(1)中获取的电压数据在时序上做一阶前向差分，其计算表达式为：Vi_Diff＝Vi(t)-Vi(t-1)其中Vi(t)表示节点i在t时刻的电压，Vi(t-1)表示节点i在t-1时刻的电压，计算结果存储在Vi_Diff中，所述Vi_Diff表示节点i的电压前向差分数据；步骤2.2：在完成步骤2.1的一阶差分之后，对得到的电压数据Vi_Diff进行分箱处理，将离散的数据变换为0,1,2，…，N这样的整数编码，分箱处理过程如下：其中，Vmax_Diff和Vmin_Diff分别指一阶前向差分后一个时间截面内各节点电压的最大值和最小值，[]为取整符号，计算结果存储在Bin_Vi中，Bin_Vi表示节点i的电压前向差分数据得到的分箱数据；步骤2.3：在上述两步数据预处理之后，最后对数据Bin_Vi进行独热编码，采用N位状态寄存器来对N个状态进行编码，每个状态都有独立的寄存器位，并且在任意时候只有一位有效，第n个状态对应的n位为1，其他位为零，将编码后的i节点的信息存储在Xi中，作为节点i的输入变量。作为本专利技术所述基于马尔可夫随机场的低压配电系统网络拓扑识别方法的一种优选方案，其中：所述步骤2.2中，将电压数据分成了0到20之间的21种不同的状态，Bin_Vi＝n，n＝0，1，2，…，20。作为本专利技术所述基于马尔可夫随机场的低压配电系统网络拓扑识别方法的一种优选方案，其中：步骤(3)的马尔可夫随机场建模具体包括：步骤3.1：将配电系统网络结构定义为一个四元组M＝(X,E,Φ,Ψ)，其中(X，E)是无向的随机变量图像；X为图的顶点集，表示配电系统中的节点，E为图的边集，表示连接配电系统节点的电力线路；然后分别建立描述配电系统节点与连线关系的势函数Φ、Ψ，来定量刻画配电网节点和连线的变量特征，描述配电系统节点与连线关系的势函数分别定义为：Φ(Xs)＝exp{E(Xs)}Ψ(Xst)＝exp{E(Xs,Xt)}式中，Xs和Xt分别是代表节点x和节点s的随机变量，Xst表示连接节点s和节点t的边缘所对应的随机变量，E称为能量函数，参数化形式如下所示：E(Xs)＝Vsds(Xs)式中ds称为状态特征，是定义在节点上的特征函数，只考虑当前节点，bst称为转移特征，是定义在边缘上的特征函数，Vs是配电系统中节点s的权重，Wst是连接节点s、t线路的权重；步骤3.2：在完成对单个节点和对相连节点的建模后，接下来建立能够描述整个配电系统中所有节点，即整个拓扑网络相关关系的联合概率分布，如下：其中，Z为分区函数，是一个规范化因子、归一化函数，其作用是保证P(Y)构成一个联合概率分布，定义为所有可能赋值的总和。作为本专利技术所述基于马尔可夫随机场的低压配电系统网络拓扑识别方法的一种优选方案，其中：步骤(4)的配电系统网络结构联合概率分布的求解具体包括：步骤4.1：利用最大似然估计法求解联合概率分布，如下:其中，ll(Θ)＝logP(Y)是P(Y)的对数似然形式；步骤4.2，根据上步求解得到的节点相关概率，使用梯度下降算法，根据步骤(1)中获取的电压数据，求解目标函数得到边权重参数W；步骤4.3，将权重参数Wijst根据合成最后得到表征配电网中各个节点连接关系的节点相关性矩阵K，其中，K是节点相关性矩阵，Kst是K中元素，Wijst是连接节点s、t线路的权重，e是自然常数。作为本专利技术所述基于马尔可夫随机场的低压配电系统网络拓扑识别方法的一种优选方案，其中：步骤(5)将节点相关性矩阵K中相关性最大的节点相连，即将最大相关性所在位置元素置1，其他置0，最终得到低压配电系统网络结构的邻接矩阵。本专利技术的有益效果：传统的配电网结构识别方法需要知道部分配电网络结构而本专利技术能够根据已知的历史电压数据直接生成低压配电网的拓扑结构，本专利技术利用历史电表电压数据等已知数据对低压配电网络结构展开识别，不需要新增设备，减少了投资，具有较高的识别准确率。附图说明...

【技术保护点】
1.一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法，其特征在于：包括如下步骤：/n(1)获取用户智能电表电压时间序列数据及配变TTU电压时间序列数据；/n(2)对步骤(1)中获取的电压数据进行一阶差分、分箱处理、独热编码处理；/n(3)利用马尔可夫随机场对配电系统网络结构建模，建立描述配电系统网络结构中各节点之间相关关系的联合概率分布；/n(4)通过最大似然估计法求解联合概率分布，得到配电系统网络结构节点间的相关性矩阵；/n(5)根据节点相关性矩阵，建立配电系统网络结构的邻接矩阵。/n

【技术特征摘要】
1.一种基于马尔可夫随机场的低压配电系统网络拓扑识别方法，其特征在于：包括如下步骤：
(1)获取用户智能电表电压时间序列数据及配变TTU电压时间序列数据；
(2)对步骤(1)中获取的电压数据进行一阶差分、分箱处理、独热编码处理；
(3)利用马尔可夫随机场对配电系统网络结构建模，建立描述配电系统网络结构中各节点之间相关关系的联合概率分布；
(4)通过最大似然估计法求解联合概率分布，得到配电系统网络结构节点间的相关性矩阵；
(5)根据节点相关性矩阵，建立配电系统网络结构的邻接矩阵。


2.如权利要求1所述的基于马尔可夫随机场的低压配电系统网络拓扑识别方法，其特征在于：所述步骤(1)中，数据采集的要求样本数在500个时间截面以上。


3.如权利要求1或2所述的基于马尔可夫随机场的低压配电系统网络拓扑识别方法，其特征在于：步骤(2)的数据预处理具体包括：
步骤2.1：对步骤(1)中获取的电压数据在时序上做一阶前向差分，其计算表达式为：
Vi_Diff＝Vi(t)-Vi(t-1)
其中Vi(t)表示节点i在t时刻的电压，Vi(t-1)表示节点i在t-1时刻的电压，计算结果存储在Vi_Diff中，所述Vi_Diff表示节点i的电压前向差分数据；
步骤2.2：在完成步骤2.1的一阶差分之后，对得到的电压数据Vi_Diff进行分箱处理，将离散的数据变换为0,1,2，…，N这样的整数编码，分箱处理过程如下：



其中，Vmax_Diff和Vmin_Diff分别指一阶前向差分后一个时间截面内各节点电压的最大值和最小值，[]为取整符号，计算结果存储在Bin_Vi中，Bin_Vi表示节点i的电压前向差分数据得到的分箱数据；
步骤2.3：在上述两步数据预处理之后，最后对数据Bin_Vi进行独热编码，采用N位状态寄存器来对N个状态进行编码，每个状态都有独立的寄存器位，并且在任意时候只有一位有效，第n个状态对应的n位为1，其他位为零，将编码后的i节点的信息存储在Xi中，作为节点i的输入变量。


4.如权利要求3所述的基于马尔可夫随机场的低压配电系统网络拓扑识别方法，其特征在于：所述步骤2.2中，将电压数据分成了0到20之间的21种不同的状态，Bin_Vi＝n，n＝0，1，2，…，20。


5.如权利要求3或4所述的基于马尔可夫随机场的...

【专利技术属性】
技术研发人员：赵健，李梁，徐明昕，张秋石，王小宇，边晓燕，
申请(专利权)人：上海电力大学，
类型：发明
国别省市：上海;31