CN111428977A - 一种基于电压序列灰色关联度的离群配电变压器辨识方法 - Google Patents

Abstract

本发明公开了一种基于电压序列灰色关联度的离群配电变压器辨识方法，具体步骤有：1、将配电网中10kV馈线下所有配电变压器视为一个群集，并从电力营销部门用电采集系统导出群集所有配电变压器出口电压时间序列数据；2、对数据进行预处理，包括完整度计算、时标对齐、三相电压平衡归算；3、对预处理后的电压时间序列数据进行灰色关联度计算；4、绘制灰色关联度图谱；5、筛选离群配电变压器。本发明基于现有配电台区用电采集系统数据进行离群配电变压器辨识，利用线上数据分析方式替代现有人工实地巡线排查方式，大幅降低人工成本，能简单、实时、有效地解决10kV配电网拓扑关系混乱的问题。

Description

一种基于电压序列灰色关联度的离群配电变压器辨识方法

技术领域

本发明属于配电网技术领域，具体涉及一种基于电压序列灰色关联度的离群配电变压器辨识方法。

背景技术

我国配电网的特点是点多面广、改造频繁，导致其网络拓扑时常出现混乱。同时因均衡负荷、降低线损、可靠供电、经济运行等需求，需对配电网拓扑结构进行适时调整，导致配电网拓扑结构是动态的。实时、准确的网络拓扑是配电网开展网损计算、故障研判、停电管理等业务的基础。国家电网公司《泛在电力物联网建设大纲》明确要求实现“站-线-变-户”关系实时准确。因此实时准确的配电网拓扑是建设泛在电力物联网的内涵所在，是泛在电力物联网平台层建设的基础内容之一。

配电网拓扑关系主要包括变电站与10kV馈线、10kV馈线与配电变压器、配电变压器与负载用户以及用户相序四种。其中，10kV馈线与配电变压器拓扑关系是当前研究的(线变拓扑关系)难点问题，当前主要依靠人工巡线排查的方式进行线变拓扑关系辨识，人工巡线排查方式不仅实施难度大耗时费力，而且难以确保时效性。需耗费大量人力物力进行巡线且效果难保证时效性，尤其是对于城区配电网，大量敷设地下电缆，人工巡线排查实施难度较大。

当前解决上述问题，现有方法有诸如专利文件(CN107508297A)基于皮尔逊线性相关系数的线变关系校验方法，与之相比，本发明具有更高的辨识精度。再比如专利文件(CN107832927A)基于灰色关联分析法的10kV线路线变关系评价方法，是基于灰色关联对电能量及相关线损量进行关联分析，而本发明则是对配电变压器出口电压时间序列进行分析，两者存在本质差别。

发明内容

本发明提出了一种基于电压序列灰色关联度的离群配电变压器辨识方法，以解决上述背景技术中提出的问题。本发明的技术方案是这样实现的：

一种基于电压序列灰色关联度的离群配电变压器辨识方法，包含如下几个步骤：

步骤一：从用电信息采集系统中提取10kV馈线所挂n台配电变压器在时段T内的电压时间序列，

假设时段T为某日，电压时间序列的数据长度为m，则第i台配电变压器导出的三相出口电压序列U_iA、U_iB、U_iC分别为

U_iA＝{U_iA(k)|k＝1,2,3,…,m}

U_iB＝{U_iB(k)|k＝1,2,3,…,m}

U_iC＝{U_iC(k)|k＝1,2,3,…,m}

式中，k表示电压序列元素序号，其余各配电变压器出口数据均参照第i台配电变压器进行处理；

步骤二：对配电变压器三相出口电压时间序列数据进行数据预处理，具体包括计算完整度、时标对齐和三相电压平衡归算，

(1)计算电压时间序列的完整度是指校验导出的电压数据有效数据的占比，具体公式如下所示：

式中，a表示单相电压序列中有效元素的数量，当得到的完整度达不到计算要求，重新从系统中导出另一日数据，直至满足计算要求，

(2)将馈线下n台配电变压器电压序列中各元素的时标对齐一致，

(3)三相电压平衡归算，是指将已经偏移的中性点还原，以消除三相不平衡的影响，中性点偏移前、后，相电压因负荷不同发生改变，而线电压并未改变，因此根据三相电压的相量三角形的封闭性，可得出偏移前、后第i台配变的第k组三相电压满足以下方程组关系：

式中，

为A相和B相、B相和C相、C相和A相之间的相位角，U_iAB、U_iBC、U_iCA分别归算的A相和B相、B相和C相、C相和A相之间的线电压，三者幅值相等，通过求解上述方程组，得到三相平衡归算后的电压序列，

将第i台配电变压器全部同一时刻的三相电压数据带入上述方程组求解，即可得到三相平衡归算后的电压序列为：

U_i＝{U_i(k)|k＝1,2,3,…,m}

式中，U_i为三相平衡计算后的电压出口序列，其余n-1台配电变压器出口数据均参照第i台配电变压器进行处理：

U＝{U_i(k)|k＝1,2,3,…,m；i＝1,2,3…n}

式中，U为所有三相平衡计算后的电压出口序列集合；

步骤三：灰色关联度计算，令预处理后的第i台配电变压器数据序列为参考序列，预处理后的所有配电变压器出口电压序列数据为比较序列，第i台参考序列与第j台比较序列在第k个元素的相对差值计算如下：

式中，ρ称为分辨系数，ρ∈(0，1)，常取0.5，表示两层的最小差，

参考序列U_i与比较序列U_j灰色关联度系数为：

式中，M_ij是参考序列U_i(k)对比较序列U_j(k)的关联度系数，

根据上述方法分别求得一条10kV馈线下任意一台与n台配电变压器之间的灰色关联度值矩阵：

式中，M是n台配电变压器各自为参考序列的灰色关联度系数矩阵；

步骤四：将得到的灰色关联度矩阵M绘制成图谱，依据图谱差异特征筛选出离群配电变压器。

在本发明的基于电压序列灰色关联度的离群配电变压器辨识方法中，步骤一中的某时段T是某日、某周、某月的任意天数。

在本发明的基于电压序列灰色关联度的离群配电变压器辨识方法中，步骤一中用电信息采集系统中的每日电压时间序列数据长度m为24的正整数倍，具体依据用电信息采集系统采集时间步长决定。

实施本发明的这种基于电压序列灰色关联度的离群配电变压器辨识方法，具有以下有益效果：

本发明基于现有配电台区的电压时间序列数据进行离群配电变压器的识别和定位算法，利用线上数据分析方式替代现有人工实地巡线排查方式，不仅大幅降低人工工作量，节省成本，而且能简单、实时、有效地解决配电网点多面广拓扑关系混乱的问题，此外，本发明方法较现有皮尔逊相关系数方法就有更高的识别准确性。

附图说明

图1为本发明基于电压序列灰色关联度的离群配电变压器辨识方法的流程图；

图2为本发明基于电压序列灰色关联度的离群配电变压器辨识方法的灰色关联度图谱。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。

如图1和2所示的一种基于电压序列灰色关联度的离群配电变压器辨识方法，具体步骤有：步骤一、从用电信息采集系统中提取10kV馈线所挂n台配电变压器在时段T内的电压时间序列，某时段T是某日、某周、某月的任意天数，假设时段T为某日，电压时间序列的数据长度为m，m为24的正整数倍，可以是24、48或96等任意数据，具体依据用电信息采集系统采集时间步长决定，则第i台配电变压器导出的三相出口电压序列U_iA、U_iB、U_iC分别为：

U_iA＝{U_iA(k)|k＝1,2,3,…,m}

U_iB＝{U_iB(k)|k＝1,2,3,…,m}

U_iC＝{U_iC(k)|k＝1,2,3,…,m}

UiC＝{UiC(k)|k＝1，2，3，…，m}

式中，k表示电压序列元素序号，其余各配电变压器出口数据均参照第i台配电变压器进行处理；步骤二、对配电变压器三相出口电压时间序列数据进行数据预处理，具体包括计算完整度、时标对齐和三相电压平衡归算，(1)计算电压时间序列的完整度是指校验导出的电压数据有效数据的占比，具体公式如下所示：

式中，a表示单相电压序列中有效元素的数量，当得到的完整度达不到计算要求，重新从系统中导出另一日数据，直至满足计算要求，(2)将馈线下n台配电变压器电压序列中各元素的时标对齐一致，(3)三相电压平衡归算，是指将已经偏移的中性点还原，以消除三相不平衡的影响，中性点偏移前、后，相电压因负荷不同发生改变，而线电压并未改变，因此根据三相电压的相量三角形的封闭性，可得出偏移前、后第i台配变的第k组三相电压满足以下方程组关系：

式中，

为A相和B相、B相和C相、C相和A相之间的相位角，U_iAB、U_iBC、U_iCA分别归算的A相和B相、B相和C相、C相和A相之间的线电压，三者幅值相等，通过求解上述方程组，得到三相平衡归算后的电压序列，将第i台配电变压器全部同一时刻的三相电压数据带入上述方程组求解，即可得到三相平衡归算后的电压序列为U_i＝{U_i(k)|k＝1,2,3,...,m}式中，U_i为三相平衡计算后的电压出口序列，其余n-1台配电变压器出口数据均参照第i台配电变压器进行处理：U＝{U_i(k)|k＝1,2,3,…,m；i＝1,2,3…n}式中，U为所有三相平衡计算后的电压出口序列集合；步骤三、灰色关联度计算，令预处理后的第i台配电变压器数据序列为参考序列，预处理后的所有配电变压器出口电压序列数据为比较序列，第i台参考序列与第j台比较序列在第k个元素的相对差值计算如下：

式中，ρ称为分辨系数，ρ∈(0，1)，常取0.5，表示两层的最小差，参考序列U_i与比较序列U_j灰色关联度系数为：

式中，M_ij是参考序列U_i(k)对比较序列U_j(k)的关联度系数，根据上述方法分别求得一条10kV馈线下任意一台与n台配电变压器之间的灰色关联度值矩阵：

式中，M是n台配电变压器各自为参考序列的灰色关联度系数矩阵；步骤四、将得到的灰色关联度矩阵M绘制成图谱，依据图谱差异特征筛选出离群配电变压器。

本发明运用在一个具体的实际场景中。以某地的3条10kV馈线为例，进行离群变压器的辨识。三条线路下共挂有配电变压器47台，其中北环Ⅰ线有22台配电变压器、青宾线17台配电变压器、北清Ⅳ线有8台配电变压器，共分为三个群集，用电信息采集系统中对应的拓扑关系如表1所示。

表1

序号	10kV线路	配电变压器	序号	10kV线路	配电变压器
						1	10kV北环Ⅰ线	机械局01号箱变	25	10kV青宾线	10kV青宾线13B青山湖区检察院公变
2	10kV北环Ⅰ线	七中01号箱变	26	10kV青宾线	10kV青宾线10B阳明路电力大楼公变
						3	10kV北环Ⅰ线	10kV北环一线环城路北08号公变	27	10kV青宾线	10kV青宾线01B永外正街一号公变
4	10kV北环Ⅰ线	创业小区01号室变	28	10kV青宾线	10kV青宾线02B永外正街二号公变
						5	10kV北环Ⅰ线	10kV北环一线环城路北01号公变	29	10kV青宾线	10kV青宾线03B省车辆管理大楼公变
6	10kV北环Ⅰ线	10kV北环一线环城路北02号公变	30	10kV青宾线	10kV青宾线06B贤士一路三号公变
						7	10kV北环Ⅰ线	10kV北环一线洪城巷03号公变	31	10kV青宾线	10kV青宾线07B贤士一路四号公变
8	10kV北环Ⅰ线	10kV北环一线新赣南路03号公变	32	10kV青宾线	10kV青宾线08B贤士一路五号公变
						9	10kV北环Ⅰ线	10kV北环一线环城路北03号公变	33	10kV青宾线	10kV央央春天线01B央央春天小区二号配电室公变
10	10kV北环Ⅰ线	10kV北环一线环城路北04号公变	34	10kV青宾线	10kV央央春天线05B央央春天小区三号配电室公变
						11	10kV北环Ⅰ线	机械局02号箱变	35	10kV青宾线	10kV央央春天线02B央央春天小区二号配电室公变
12	10kV北环Ⅰ线	吉祥花园01号箱变	36	10kV青宾线	10kV央央春天线03B央央春天小区二号配电室公变
						13	10kV北环Ⅰ线	洪城巷01号箱变	37	10kV青宾线	10kV央央春天线04B央央春天小区二号配电室公变
14	10kV北环Ⅰ线	洪城巷02号箱变	38	10kV青宾线	10kV青宾线14B永外正街四号公变
						15	10kV北环Ⅰ线	环城路北05号箱变	39	10kV青宾线	10kV青宾线15号公变
16	10kV北环Ⅰ线	藕塘里01号公变	40	10kV北清Ⅳ线	10kV北清四线大公路东04号公变
						17	10kV北环Ⅰ线	吉祥花园03号箱变	41	10kV北清Ⅳ线	10kV北清四线大公路东03号公变
18	10kV北环Ⅰ线	环城路北09号箱变	42	10kV北清Ⅳ线	10kV北清四线大公路东02号公变
						19	10kV北环Ⅰ线	吉祥花园02号箱变	43	10kV北清Ⅳ线	10kV北清四线文清路08号公变
20	10kV北环Ⅰ线	环城路北#10公变	44	10kV北清Ⅳ线	10kV北清四线文清路07号公变
						21	10kV北环Ⅰ线	祥和2号公变	45	10kV北清Ⅳ线	10kV北清四线新赣南路04号公变
22	10kV北环Ⅰ线	祥和1号公变	46	10kV北清Ⅳ线	10kV北清四线新赣南路01号公变
						23	10kV青宾线	10kV青宾线11B永外正街三号公变	47	10kV北清Ⅳ线	10kV北清四线新赣南路02号公变
24	10kV青宾线	10kV青宾线12B永外正街与贤士横街公变

根据步骤一，在用电信息采集系统中导出3条10kV馈线下所有47台配电变压器在2019年8月1日的出口电压时间序列，由于篇幅限制，仅给出第1台配电变压器的出口电压序列数据，如表2所示。

表2

根据步骤二，对提取的电压数据对配电变压器三相出口电压时间序列数据进行数据预处理，具体包括计算完整度、时标对齐和三相电压平衡归算。

(1)所导出的配电变压器台电压序列数据长度m为47，其数据完整度计算如下：

三个配电变压器群集的电压时间序列均满足完整度满足要求。

(2)按照时间序列时标，将三个配电变压器群集电压序列中各元素的时标对齐一致。

(3)三相电压平衡归算后的电压序列，由于篇幅限制，仅给出第1台配电变压器的出口电压序列数据，如表2所示

表2

根据步骤三,根据上述方法分别求得三个10kV馈线配配电变压器群集的灰色关联度值矩阵M，如表3所示：

表3

续

根据步骤四，将灰色关联度矩阵M绘制成相应的灰色关联度图谱，如图2所示，通过观察灰色关联度图谱，可得出：

(1)9号、13号配电变压器与北环Ⅰ线群集中其他配电变压器相差较大，为该群集的离群配电变压器。

(2)23号配电变压器与青宾线群集中其他配电变压器相差较大，为该群集的离群配电变压器。

(3)43号配电变压器与北清Ⅳ线群集中其他配电变压器相差较大，为该群集的离群配电变压器。

实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于电压序列灰色关联度的离群配电变压器辨识方法，其特征在于，包含如下几个步骤：

步骤一：从用电信息采集系统中提取10kV馈线所挂n台配电变压器在时段T内的电压时间序列，

假设时段T为某日，电压时间序列的数据长度为m，则第i台配电变压器导出的三相出口电压序列

、

、

分别为

式中，k表示电压序列元素序号，其余各配电变压器出口数据均参照第i台配电变压器进行处理；

步骤二：对配电变压器三相出口电压时间序列数据进行数据预处理，具体包括计算完整度、时标对齐和三相电压平衡归算，

（1）计算电压时间序列的完整度是指校验导出的电压数据有效数据的占比，具体公式如下所示：

式中，a表示单相电压序列中有效元素的数量，当得到的完整度达不到计算要求，重新从系统中导出另一日数据，直至满足计算要求，

（2）将馈线下n台配电变压器电压序列中各元素的时标对齐一致，

（3）三相电压平衡归算，是指将已经偏移的中性点还原，以消除三相不平衡的影响，中性点偏移前、后，相电压因负荷不同发生改变，而线电压并未改变，因此根据三相电压的相量三角形的封闭性，可得出偏移前、后第i台配变的第k组三相电压满足以下方程组关系：

式中，

为A相和B相、B相和C相、C相和A相之间的相位角，U_iAB、U_iBC、U_iCA分别归算的A相和B相、B相和C相、C相和A相之间的线电压，三者幅值相等，通过求解上述方程组，得到三相平衡归算后的电压序列，

将第i台配电变压器全部同一时刻的三相电压数据带入上述方程组求解，即可得到三相平衡归算后的电压序列为：

式中，

为三相平衡计算后的电压出口序列，其余n-1台配电变压器出口数据均参照第i台配电变压器进行处理：

式中，U为所有三相平衡计算后的电压出口序列集合；

步骤三：灰色关联度计算，令预处理后的第i台配电变压器数据序列为参考序列，预处理后的所有配电变压器出口电压序列数据为比较序列，第i台参考序列与第j台比较序列在第k个元素的相对差值计算如下：

式中，ρ称为分辨系数，ρ∈（0，1），常取0.5，表示两层的最小差，

参考序列U_i与比较序列U_j灰色关联度系数为：

式中，M_ij是参考序列U_i(k)对比较序列U_j(k)的关联度系数，

根据上述方法分别求得一条10kV馈线下任意一台与n台配电变压器之间的灰色关联度值矩阵：

式中，

是n台配电变压器各自为参考序列的灰色关联度系数矩阵；

步骤四：将得到的灰色关联度矩阵

绘制成图谱，依据图谱差异特征筛选出离群配电变压器。

2.根据权利要求1所述的基于电压序列灰色关联度的离群配电变压器辨识方法，其特征在于，步骤一中的某时段T是某日、某周、某月的任意天数。

3.根据权利要求2所述的基于电压序列灰色关联度的离群配电变压器辨识方法，其特征在于，步骤一中用电信息采集系统中的每日电压时间序列数据长度m为24的正整数倍，具体依据用电信息采集系统采集时间步长决定。

Images