本发明公开了一种数据驱动的配电网线变关系诊断方法、装置及系统，方法包括计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数；筛选出最大相关系数，作为各配电变压器的中级相关系数；基于虚拟代理思想修正所有中级相关系数，得到各配电变压器的高级相关系数；基于高级相关系数，计算出高级相关系数阈值，筛选出线变关系疑似异常配电变压器；修正线变关系疑似异常配电变压器修正到最有可能所属馈线，若修正后所属的馈线与原始所属的馈线不同，判定该疑似异常配电变压器为线变关系错误的配电变压器。本发明可用于识别任意规模的配电网线变关系错误，计算简单，原理明确，可帮助配网运行人员及时发现线变关系连接错误。
数据驱动的配电网线变关系诊断方法、装置及系统
技术领域
[0001] 本发明属于配电网线变关系诊断技术领域，具体涉及一种数据驱动的配电网线变关系诊断方法、装置及系统，尤其涉及一种中压数据驱动的配电网线变关系诊断方法、装置
及系统。
背景技术
[0002] 中压配电网线变关系由中压馈线、配电变压器负荷、联络开关和输电线路等相互连接组成，描述了中压馈线挂接配电变压器负荷的关系，由实际现场人员人工维护至能量
管理系统(Power Management System，PMS)。但是由于实际运行过程中，大量联络开关频繁
切换和更新配置，以及量测装置准确性低、人工操作主观性差等因素影响，使得当前配电网
线变关系难以与现场匹配一致。
[0003] 传统依赖直接考虑馈线首端与配电变压器负荷电压变化相似性能够识别部分线变关系错误现象，但传统方法未考虑2方面问题，一是配电变压器量测数据采样时间与馈线
量测数据采样时间不同步问题，本发明通过自适应时间窗滑动在量测数据“前移后推”过程
中寻找匹配最优的时间窗，并给出配电变压器量测数据时钟调整建议；二是如果线路较长、
配电变压器数量较多，将处于线路末端的配电变压器直接匹配首端母线电压变化的趋势并
不公平，本发明通过虚拟代理思想，从馈线首端出发逐步寻找优秀虚拟代理，形成代理链
条，克服线路过长引起的相关系数指标误差，精确实现线变关系异常诊断。
[0004] 正确的配电网线变关系是电力系统线损计算、规划设计和稳定运行的研究基础，没有一个及时准确的配电网线变关系，线损计算等重要工作将难以推进。因此，中压配电网
网络线变关系诊断方法是一个重要的研究课题，研究成果可以帮助运行人员及时发现线变
关系错误，进一步挖掘配网信息。
发明内容
[0005] 针对上述问题，本发明提出一种数据驱动的配电网线变关系诊断方法、装置及系统，可用于识别任意规模的中压配电网线变关系错误，该方法计算简单，原理明确，可帮助
配网运行人员及时发现线变关系连接错误。
[0006] 为了实现上述技术目的，达到上述技术效果，本发明通过以下技术方案实现：
[0007] 第一方面，本发明提供了一种数据驱动的配电网线变关系诊断方法，包括：
[0008] 基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数；
[0009] 筛选出各配电变压器与其原始所属馈线之间的最大相关系数，作为各配电变压器与其原始所属馈线之间的中级相关系数；
[0010] 基于虚拟代理思想，修正所有中级相关系数，得到各配电变压器与其原始所属馈线之间的高级相关系数；
[0011] 基于所述高级相关系数，计算出高级相关系数阈值，并筛选出线变关系疑似异常配电变压器；
[0012] 将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，若修正后所属的馈线与原始所属的馈线不同，则判定该疑似异常配电变压器为线变关系错误的配电变压
器。
[0013] 可选地，假设时间窗偏移数为m，所述不同时间窗内各配电变压器与其原始所属馈线之间的相关系数的计算方法包括：
[0014] 针对各配电变压器，分别取出其当前时间窗、以及超前和滞后m个时间窗的电压序列，计算出各时间窗内配电变压器与其原始所属馈线的相关系数，得到与各配电变压器对
应的线变关系指标集；
[0015] 所述中级相关系数的计算方法包括：
[0016] 筛选出各线变关系指标集中的最大值，作为各配电变压器与其原始所属馈线之间的中级相关系数。
[0017] 可选地，假设某配电变压器Tk连接在馈线f上，配电变压器Tk的电压量测数据为Vk＝[vk,1,vk,2,…,vk,i,…,vk,n]，馈线f的电压量测数据为Vf＝[vf,1,vf,2,…,vf,i,…,vf,n]，n
表示配电变压器与馈线电压量测数据采样点个数；
[0018] 所述配电变压器与其原始所属馈线的相关系数的计算方法包括：
[0019] 采用相关系数计算公式，计算出配电变压器与其原始所属的馈线之间的相关系数，所述相关系数计算公式为：
[0020]
[0021] 其中，Rk,f表示配电变压器Tk与其原始所属馈线f的相关系数； 表示配电变压器Tk的电压量测数据平均值； 表示馈线f的电压量测数据平均值。
[0022] 可选地，各配电变压器与其原始所属馈线之间的高级相关系数的计算方法包括：
[0023] 设共有g个配电变压器原始所属于馈线f，利用该g个配电变压器与其原始所属馈线之间的中级相关系数构建馈线f的中级相关系数集合R(2)；
[0024] 挑选出中级相关系数集合R(2)里中级相关系数大于设定虚拟代理阈值c的配电变压器，作为第一层虚拟代理，并将第一层虚拟代理中的各配电变压器与其原始所属馈线之
间的中级相关系数作为其与原始所属馈线之间的高级相关系数；
[0025] 对于除第一层虚拟代理以外的配电变压器，重复执行以下步骤，直至未找到下一层虚拟代理为止：
[0026] 以前一次计算出的第a层虚拟代理作为虚拟馈线，计算出剩余各配电变压器与第a层虚拟代理之间的中级相关系数，取出中级相关系数大于c的配电变压器作为第a+1层虚拟
代理，并将第a+1层虚拟代理中各配电变压器与第a层虚拟代理之间计算出的中级相关系
数，作为第a+1层虚拟代理中各配电变压器与原始所属馈线之间的高级相关系数。
[0027] 可选地，所述高级相关系数阈值的计算方法包括：
[0028] 将所有配电变压器与其原始所属馈线之间的高级相关系数组成指标集合；
[0029] 按照设定的区间步长将统计区间划分成若干个子区间，计算指标集合在各子区间的累计分布概率，形成累计分布概率集合；
[0030] 选择直线y＝Ax+B对所述累计分布概率集合的下降段拟合的平均误差最小时，拟合直线y＝Ax+B与x轴的交点作为高级相关系数阈值。
[0031] 可选地，记累计分布概率集合为I＝[I1,I2,…，Ib,…,It]，其中Ib表示指标集合中数据分布在第b个子区间的概率；
[0032] 所述选择直线y＝Ax+B对所述累计分布概率集合的下降段拟合的平均误差最小时，拟合直线y＝Ax+B与x轴的交点作为高级相关系数阈值，包括以下步骤：
[0033] 利用数据点(b,Ib)形成数据集S，b＝1,2,…,t；
[0034] 置初始待拟合数据规模为q，记录数据集S中的后q个数据点为数据集Sq；
[0035] 基于最小二乘法对数据集Sq作拟合，得到拟合曲线y＝Aqx+Bq；
[0036] 计算当x＝t‑q,t+1‑q,…,t时对应的拟合值y，并记录为数据集Sq*；
[0037] 按 计算选择q个数据点时拟合的平均误差Γq，并记入数据集Γ中；
[0038] 令q＝q+1，当q>t时，找到集合Γ中的最小值，其对应的直线y＝Aqx+Bq与x轴的交点即为最终选定的高级相关系数阈值。
[0039] 可选地，所述将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，包括以下步骤：
[0040] 利用经纬度信息，筛选出地理位置相近的馈线，形成目标馈线集合F＝[f1,f2,…,fl]，其中l表示地理位置相近的馈线条数；
[0041] 筛选出各馈线fd上线变关系正确的配电变压器和疑似异常配电变压器，d＝1,2,…,l；
[0042] 对于各馈线fd，将线变关系正确的配电变压器归为一类，并计算该类中配电变压器与馈线的相关系数平均值，作为该类的类别中心；
[0043] 计算各疑似异常配电变压器与各类别中心之间的欧式距离；
[0044] 基于欧式距离最小，依次确定出与各疑似异常配电变压器对应的类别，将疑似异常配电变压器修正到对应类中，更新对应类的聚类中心，完成将线变关系疑似异常配电变
压器修正到最有可能所属馈线中去。
[0045] 第二方面，本发明提供了一种数据驱动的配电网线变关系诊断装置，包括：
[0046] 计算单元，用于基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数；
[0047] 第一筛选单元，用于筛选出各配电变压器与其原始所属馈线之间的最大相关系数，作为各配电变压器与其原始所属馈线之间的中级相关系数；
[0048] 修正单元，用于基于虚拟代理思想，修正所有中级相关系数，得到各配电变压器与其原始所属馈线之间的高级相关系数；
[0049] 第二筛选单元，用于基于所述高级相关系数，计算出高级相关系数阈值，并筛选出线变关系疑似异常配电变压器；
[0050] 判断单元，用于将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，若修正后所属的馈线与原始所属的馈线不同，则判定该疑似异常配电变压器为线变关系错
误的配电变压器。
[0051] 可选地，假设时间窗偏移数为m，所述计算单元包括：
[0052] 提取模块，用于针对各配电变压器，分别取出其当前时间窗、以及超前和滞后m个时间窗的电压序列；
[0053] 计算模块，用于计算出各时间窗内配电变压器与其原始所属馈线的相关系数，得到与各配电变压器对应的线变关系指标集；
[0054] 所述第一筛选单元包括：
[0055] 筛选模块，用于筛选出各线变关系指标集中的最大值；
[0056] 赋值模块，用于将各线变关系指标集中的最大值作为各配电变压器与其原始所属馈线之间的中级相关系数。
[0057] 可选地，设共有g个配电变压器原始所属于馈线f，所述修正单元包括：
[0058] 构建模块，用于利用该g个配电变压器与其原始所属馈线之间的中级相关系数构建馈线f的中级相关系数集合R(2)；
[0059] 第一挑选模块，用于挑选出中级相关系数集合R(2)里中级相关系数大于设定虚拟代理阈值c的配电变压器，作为第一层虚拟代理，并将第一层虚拟代理中的各配电变压器与
其原始所属馈线之间的中级相关系数作为其与原始所属馈线之间的高级相关系数；
[0060] 第二挑选模块，用于对于除第一层虚拟代理以外的配电变压器，重复执行以下步骤，直至未找到下一层虚拟代理为止：
[0061] 以前一次计算出的第a层虚拟代理作为虚拟馈线，计算出剩余各配电变压器与第a层虚拟代理之间的中级相关系数，取出中级相关系数大于c的配电变压器作为第a+1层虚拟
代理，并将第a+1层虚拟代理中各配电变压器与第a层虚拟代理之间计算出的中级相关系
数，作为第a+1层虚拟代理中各配电变压器与原始所属馈线之间的高级相关系数。
[0062] 可选地，所述高级相关系数阈值的计算方法包括：
[0063] 将所有配电变压器与其原始所属馈线之间的高级相关系数组成指标集合；
[0064] 按照设定的区间步长将统计区间划分成若干个子区间，计算指标集合在各子区间的累计分布概率，形成累计分布概率集合；
[0065] 选择直线y＝Ax+B对所述累计分布概率集合的下降段拟合的平均误差最小时，拟合直线y＝Ax+B与x轴的交点作为高级相关系数阈值。
[0066] 可选地，所述将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，包括以下步骤：
[0067] 利用经纬度信息，筛选出地理位置相近的馈线，形成目标馈线集合F＝[f1,f2,…,fl]，其中l表示地理位置相近的馈线条数；
[0068] 筛选出各馈线fd上线变关系正确的配电变压器和疑似异常配电变压器，d＝1,2,…,l；
[0069] 对于各馈线fd，将线变关系正确的配电变压器归为一类，并计算该类中配电变压器与馈线的相关系数平均值，作为该类的类别中心；
[0070] 计算各疑似异常配电变压器与各类别中心之间的欧式距离；
[0071] 基于欧式距离最小，依次确定出与各疑似异常配电变压器对应的类别，将疑似异常配电变压器修正到对应类中，更新对应类的聚类中心，完成将线变关系疑似异常配电变
压器修正到最有可能所属馈线中去。
[0072] 第三方面，本发明提供了一种数据驱动的配电网线变关系诊断系统，包括：存储介质和处理器；
[0073] 所述存储介质用于存储指令；
[0074] 所述处理器用于根据所述指令进行操作以执行根据第一方面中任一项所述方法的步骤。
[0075] 与现有技术相比，本发明的有益效果：
[0076] (1)本发明考虑了配电变压器量测数据采样时间与馈线量测数据采样时间不同步问题，提出了基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所
属的馈线之间的相关系数，筛选出各配电变压器与其原始所属的馈线之间的最大相关系
数，作为各配电变压器与其原始所属的馈线之间的中级相关系数，中级相关系数对应的时
间窗即为各配电变压器的最优时间窗，根据最优时间窗给出各配电变压器采样时钟调整建
议。
[0077] (2)本发明考虑了线路较长、配电变压器数量较多时，处于线路末端的配电变压器直接匹配首端母线电压变化的趋势并不公平的问题，通过虚拟代理思想，从馈线首端出发
逐步寻找优秀虚拟代理，形成代理链条，克服线路过长引起的相关系数指标误差。
[0078] (3)本发明通过学习历史数据特征，通过概率分布估计方法计算评判线变关系是否异常的阈值，避免了人为指定。
[0079] (4)本发明基于类别修正方法，避免了直接聚类时初始化阶段的随机性，同时避免了直接聚类时若线变关系疑似异常配电变压器数量上占据主导地位时，聚类中心会向线变
关系疑似异常配电变压器方向偏移的问题。
[0080] (5)本发明计算简单，原理明确，可帮助配网运行人员及时发现线变关系错误，具有良好的应用前景。
附图说明
[0081] 为了使本发明的内容更容易被清楚地理解，下面根据具体实施例并结合附图，对本发明作进一步详细的说明，其中：
[0082] 图1为本发明的配电网线变关系诊断的流程图；
[0083] 图2为奇异值分解方法的原理图；
[0084] 图3为自适应时间窗滑动方法中时间窗滑动示意图。
具体实施方式
[0085] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于
限定本发明的保护范围。
[0086] 下面结合附图对本发明的应用原理作详细的描述。
[0087] 实施例1
[0088] 本发明实施例中提供了一种数据驱动的配电网线变关系诊断方法，包括以下步骤：
[0089] 基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数；
[0090] 筛选出各配电变压器与其原始所属馈线之间的最大相关系数，作为各配电变压器与其原始所属馈线之间的中级相关系数；
[0091] 基于虚拟代理思想，修正所有中级相关系数，得到各配电变压器与其原始所属馈线之间的高级相关系数；
[0092] 基于所述高级相关系数，计算出高级相关系数阈值，并筛选出线变关系疑似异常配电变压器；
[0093] 将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，若修正后所属的馈线与原始所属的馈线不同，则判定该疑似异常配电变压器为线变关系错误的配电变压
器。
[0094] 在本发明实施例的一种具体实施方式中，假设时间窗偏移数为m，所述不同时间窗内各配电变压器与其原始所属馈线之间的相关系数的计算方法包括：
[0095] 针对各配电变压器，分别取出其当前时间窗、以及超前和滞后m个时间窗的电压序列，计算出各时间窗内配电变压器与其原始所属馈线的相关系数，得到与各配电变压器对
应的线变关系指标集。
[0096] 所述中级相关系数的计算方法包括：
[0097] 筛选出各线变关系指标集中的最大值，作为各配电变压器与其原始所属馈线之间的中级相关系数。
[0098] 在具体实施过程中，所述基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数，可以采用以下具体步骤来实现：
[0099] 步骤1：针对配电变压器Tk，其电压量测数据为Vk＝[vk,1,vk,2,…,vk,i,…,vk,n]，假设时间窗偏移数为m，vk,‑m为vk,1前第m个电压数据，vk,n+m为vk,n后第m个电压数据；
[0100] 步骤2：取出配电变压器Tk超前1个时间窗的电压序列Vk(1)＝[vk,‑1,vk,1,…,vk,n‑1]T，(1)
计算与馈线f的相关系数Rk,f ；
[0101] 步骤3：取出配电变压器Tk滞后1个时间窗的电压序列Vk(1)＝[vk,2,vk,3,…,vk,n+1]T，(‑1)
计算与馈线f的相关系数Rk,f ；
[0102] 步骤4：重复上述步骤2和步骤3，直至计算完成超前和滞后m个时间窗对应的线变(‑m) (‑m+1) (m‑1) (m)
关系指标集[Rk,f ,Rk,f ,…,Rk,f,…,Rk,f ,Rk,f ]。
[0103] 所述中级相关系数可以采用以下步骤计算得到：
[0104] 步骤5：将步骤4中获得的指标集中最大值作为配电变压器Tk的中级相关系数，对应的时间窗即为配电变压器Tk最优时间窗，可以根据最优时间窗给出配电变压器采样时钟
调整建议。
[0105] 在本发明实施例的一种具体实施方式中，假设某配电变压器Tk连接在馈线f上，配电变压器Tk的电压量测数据为Vk＝[vk,1,vk,2,…,vk,i,…,vk,n]，馈线f的电压量测数据为Vf
＝[vf,1,vf,2,…,vf,i,…,vf,n]，n表示配电变压器与馈线电压量测数据采样点个数；
[0106] 所述配电变压器与其原始所属馈线的相关系数的计算方法包括：
[0107] 采用相关系数计算公式，计算出配电变压器与其原始所属的馈线之间的相关系数；所述相关系数计算公式为：
[0108]
[0109] 其中，Rk,f表示配电变压器Tk的相关系数，其值范围在[0,1]区间内； 表示配电变压器Tk的电压量测数据平均值； 表示馈线f的电压量测数据平均值。
[0110] 在本发明实施例的一种具体实施方式中，考虑由于线路供电长度不同引起配电变压器相关系数计算误差问题，本发明提出基于虚拟代理思想修正中级相关系数，最终计算
出配电变压器与其原始所属的馈线之间的高级相关系数；具体地，所述计算出各配电变压
器与其原始所属馈线之间的高级相关系数的计算方法包括：
[0111] 设共有g个配电变压器原始所属于馈线f，利用该g个配电变压器的中级相关系数构建馈线f的中级相关系数集合R(2)；
[0112] 挑选出中级相关系数集合R(2)里中级相关系数大于设定虚拟代理阈值c的配电变压器，作为第一层虚拟代理，并将第一层虚拟代理中的各配电变压器与其原始所属馈线之
间的中级相关系数作为其与原始所属馈线之间的高级相关系数；在具体应用过程中，所述
虚拟代理阈值c可以取0.95；
[0113] 对于除第一层虚拟代理以外的配电变压器，重复执行以下步骤，直至未找到下一层虚拟代理为止：
[0114] 以前一次计算出的第a层虚拟代理作为虚拟馈线，计算出剩余各配电变压器与第a层虚拟代理之间的中级相关系数，取出中级相关系数大于c的配电变压器作为第a+1层虚拟
代理，并将第a+1层虚拟代理中各配电变压器与第a层虚拟代理中各配电变压器之间计算出
的中级相关系数，作为第a+1层虚拟代理中各配电变压器与原始所属馈线之间的高级相关
系数；
[0115] 对于未被选中成为虚拟代理的各配电变压器，则将这些配电变压器与其原始所属馈线之间的中级相关系数作为其与原始所属馈线之间的高级相关系数。
[0116] 在实际应用过程中，所述计算出各配电变压器与其原始所属馈线之间的高级相关系数的计算方法包括，具体包括以下步骤：
[0117] 步骤A：挑选出中级相关系数集合R(2)里中级相关系数大于设定虚拟代理阈值c的配电变压器，作为第一层虚拟代理，并将第一层虚拟代理中的各配电变压器与其原始所属
馈线之间的中级相关系数作为其与原始所属馈线之间的高级相关系数；在具体应用过程
中，所述虚拟代理阈值c可以取0.95；
[0118] 步骤B：以第一层虚拟代理作为虚拟馈线，计算出剩余配电变压器与第一层虚拟代理之间的中级相关系数，取出中级相关系数大于c的配电变压器作为第2层虚拟代理；所述
剩余配电变压器与第一层虚拟代理之间的中级相关系数是通过针对各个剩余配电变压器，
分别进行以下操作获得的：基于自适应时间窗滑动方法，计算出不同时间窗内某剩余配电
变压器与第一层虚拟代理中各配电变压器之间的相关系数；筛选出该配电变压器与第一层
虚拟代理中各配电变压器之间的最大相关系数，作为该配电变压器与第一层虚拟代理之间
的中级相关系数；
[0119] 步骤C：将第2层虚拟代理中各配电变压器与第一层虚拟代理间计算出的中级相关系数，作为第2层虚拟代理中各配电变压器其与原始所属馈线之间的高级相关系数；
[0120] 步骤D：重复上述步骤B‑C直到未找到下一层虚拟代理为止。
[0121] 在本发明实施例的一种具体实施方式中，提出利用累计分布概率信息，基于最小二乘法计算高级相关系数阈值，所述高级相关系数阈值的计算方法包括：
[0122] 以历史日某地市为例，假设共有中压配电变压器h台，计算出所有配电变压器与其原始所属馈线之间的高级相关系数；
[0123] 将所有配电变压器与其原始所属馈线之间的高级相关系数组成指标集合；
[0124] 按照设定的区间步长将统计区间划分成若干个子区间，计算指标集合在各子区间的累计分布概率，形成累计分布概率集合；
[0125] 选择直线y＝Ax+B对所述累计分布概率集合的下降段拟合的平均误差最小时，拟合直线y＝Ax+B与x轴的交点作为高级相关系数阈值。
[0126] 在本发明实施例的一种具体实施方式中，记累计分布概率集合为I＝[I1,I2,…，Ib,…,It]，其中Ib表示指标集合中数据分布在第b个子区间的概率；
[0127] 所述选择直线y＝Ax+B对所述累计分布概率集合的下降段拟合的平均误差最小时，拟合直线y＝Ax+B与x轴的交点作为高级相关系数阈值，包括以下步骤：
[0128] 利用数据点(b,Ib)形成数据集S，b＝1,2,…,t；
[0129] 置初始待拟合数据规模为q，记录数据集S中的后q个数据点为数据集Sq；
[0130] 基于最小二乘法对数据集Sq作拟合，得到拟合曲线y＝Aqx+Bq；
[0131] 计算当x＝t‑q,t+1‑q,…,t时对应的拟合值y，并记录为数据集Sq*；
[0132] 按 计算选择q个数据点时拟合的平均误差Γq，并记入数据集Γ中；
[0133] 令q＝q+1，当q>t时，找到集合Γ中的最小值，其对应的直线y＝Aqx+Bq与x轴的交点即为最终选定的高级相关系数阈值。
[0134] 在具体实施过程中，可以采用以下具体步骤来实现：
[0135] 以历史日某地市为例，假设共有中压配电变压器h台，计算所有中压配电变压器与对应馈线的高级相关系数，组成指标集合H。
[0136] 以0.01为区间步长，将区间[0,1]划分成[0,0.01×b]，其中b＝1,2,…,100，计算指标集合H数据在各区间的累计分布概率，形成累计分布概率集合I＝[I1,I2,…,Ib,…,
I100]，其中Ib表示集合H数据中分布在区间[0,0.01×b]的概率；
[0137] 结合最小二乘法，使直线y＝Ax+B对累计分布概率集合I数据下降段拟合的平均误差最小时，此时拟合直线与x轴的交点作为高级相关系数阈值w。具体步骤如下：
[0138] 步骤1：计算出累计分布概率集合I。记录100个数据点(b,Ib)为数据集S；
[0139] 步骤2：计算拟合曲线。置初始待拟合数据规模q＝3，记录数据集S中的后q个数据点为数据集Sq。基于最小二乘法对集合Sq作拟合，得到拟合曲线y＝Aqx+Bq。计算当x＝100‑
*
q,101‑q,…,100时对应的拟合值y，并记录为数据集Sq；
[0140] 步骤3：计算拟合误差。按下式计算选择q个数据点时拟合的平均误差Γq，并记入数据集Γ中；
[0141]
[0142] 步骤4：令q＝q+1，当q>100时，转到步骤5，否则转到步骤2；
[0143] 步骤5：找到集合Γ中的最小值，其对应的直线y＝Aqx+Bq与x轴的交点即为最终选定的高级相关系数阈值w。
[0144] 所述筛选出线变关系疑似异常配电变压器，包括以下步骤：
[0145] 若配电变压器的高级相关系数大于阈值w，则认为此配电变压器线变关系正确，若配电变压器的高级相关系数小于阈值w，则认为此配电变压器线变关系疑似异常。
[0146] 在本发明实施例的一种具体实施方式中，所述将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，包括以下步骤：
[0147] 考虑到配电变压器至中压馈线的挂接点覆盖的地理范围有限，利用经纬度信息，筛选出地理位置相近的馈线，形成目标馈线集合F＝[f1,f2,…,fl]，其中l表示地理位置相
近的馈线条数；
[0148] 筛选出各馈线fd上线变关系正确的配电变压器和疑似异常配电变压器，d＝1,2,…,l；
[0149] 对于各馈线fd，将线变关系正确的配电变压器归为一类，并计算该类中配电变压器与馈线的相关系数平均值，作为该类的类别中心；
[0150] 计算各疑似异常配电变压器与各类别中心之间的欧式距离；
[0151] 基于欧式距离最小，依次确定出与各疑似异常配电变压器对应的类别，将此疑似异常配电变压器修正到该类中，更新该类的聚类中心，完成将线变关系疑似异常配电变压
器修正到最有可能所属馈线中去。
[0152] 实施例2
[0153] 基于实施例1，本发明实施例与实施例1的区别在于：
[0154] 如图2所示，所述基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数步骤之前，还包括：
[0155] 获取中压配电网中的配电变压器与馈线的相关数据：
[0156] 对获得的配电网中的配电变压器与馈线的相关数据进行去噪处理，获得去噪后的配电变压器与馈线的相关数据。
[0157] 在本发明实施例的一种具体实施方式中，所述获取中压配电网中的配电变压器与馈线的相关数据，可以通过以下步骤实现：
[0158] 在设备管理系统(该系统为现有技术中存在的系统)中选择待处理的配电网(比如中压配电网)，读取现有系统存储的“馈线‑配变负荷”的线变关系，导出一天每15min采样一
次的配变与馈线的电压数据，即数据长度n＝96，前述的每15min采样一次的采用频率可以
根据实际情况进行修改。
[0159] 在本发明实施例的一种具体实施方式中，所述对获得的配电网中的配电变压器与馈线的相关数据进行去噪处理，具体为：
[0160] 采用奇异值分解方法对获得的配电网中的配电变压器与馈线的相关数据进行去噪处理。
[0161] 所述采用奇异值分解方法对原始采集数据作去噪处理，其基本思想是将原始数据经奇异值分解方法分解为左奇异向量、奇异值和右奇异向量相乘得到形式，较小奇异值对
应了数据噪声，保留奇异值较大值并进行数据重构，重构后的数据即为去噪后的数据，其基
本原理图如图2所示。共有g个配电变压器原始所属于馈线f，设g个配变的电压数据构成的
T g×g n
g*n阶矩阵X＝[X1,X2,...,Xg]，在奇异值分解理论中，对于矩阵X，存在正交U∈R 和V∈R
×n
，使得：
[0162]
[0163] 式中：正交矩阵U＝[u1,u2,…,ug]的列向量为单位向量且相互正交，同时是矩阵T
XX的特征向量，称为左奇异向量；正交矩阵V＝[v1,v2,…,vn]的列向量也为单位向量且相
T
互正交，同时是矩阵XX的特征向量，称为右奇异向量；Λ1＝diag(λ1,λ2,…,λp)，其中p＝min
(g,n)。Λ1是对角矩阵，其对角元素为矩阵X的奇异值并按降序排列，即λ1≥λ2≥,…,≥λp。
认为较小奇异值对应的部分代表了能够忽略的噪声。这样，由上式，仅保留前q(q
主导地位的奇异值，矩阵X和某一配变数据Xi分别约简为：
[0164]
[0165]
[0166] 在忽略掉代表噪声的奇异值后并进行重构，即可得到去噪后的数据。
[0167] 实施例3
[0168] 本发明实施例中提供了一种数据驱动的配电网线变关系诊断装置，包括：
[0169] 计算单元，用于基于自适应时间窗滑动方法，计算出不同时间窗内各配电变压器与其原始所属馈线之间的相关系数；
[0170] 第一筛选单元，用于筛选出各配电变压器与其原始所属馈线之间的最大相关系数，作为各配电变压器与其原始所属馈线之间的中级相关系数；
[0171] 修正单元，用于基于虚拟代理思想，修正所有中级相关系数，得到各配电变压器与其原始所属馈线之间的高级相关系数；
[0172] 第二筛选单元，用于基于所述高级相关系数，计算出高级相关系数阈值，并筛选出线变关系疑似异常配电变压器；
[0173] 判断单元，用于将线变关系疑似异常配电变压器修正到最有可能所属馈线中去，若修正后所属的馈线与原始所属的馈线不同，则判定该疑似异常配电变压器为线变关系错
误的配电变压器。
[0174] 在本发明实施例的一种具体实施方式中，假设时间窗偏移数为m，所述计算单元包括：
[0175] 提取模块，用于针对各配电变压器，分别取出其当前时间窗、以及超前和滞后m个时间窗的电压序列；
[0176] 计算模块，用于计算出各时间窗内配电变压器与其原始所属馈线的相关系数，得到与各配电变压器对应的线变关系指标集；
[0177] 所述第一筛选单元包括：
[0178] 筛选模块，用于筛选出各线变关系指标集中的最大值；
[0179] 赋值模块，用于将各线变关系指标集中的最大值作为各配电变压器与其原始所属馈线之间的中级相关系数。
[0180] 在本发明实施例的一种具体实施方式中，设共有g个配电变压器原始所属于馈线f，所述修正单元包括：
[0181] 构建模块，用于利用该g个配电变压器与其原始所属馈线之间的中级相关系数构建馈线f的中级相关系数集合R(2)；
[0182] 第一挑选模块，用于挑选出中级相关系数集合R(2)里中级相关系数大于设定虚拟代理阈值c的配电变压器，作为第一层虚拟代理，并将第一层虚拟代理中的各配电变压器与
其原始所属馈线之间的中级相关系数作为其与原始所属馈线之间的高级相关系数；
[0183] 第二挑选模块，用于对于除第一层虚拟代理以外的配电变压器，重复执行以下步骤，直至未找到下一层虚拟代理为止：
[0184] 以前一次计算出的第a层虚拟代理作为虚拟馈线，计算出剩余各配电变压器与第a层虚拟代理之间的中级相关系数，取出中级相关系数大于c的配电变压器作为第a+1层虚拟
代理，并将第a+1层虚拟代理中各配电变压器与第a层虚拟代理之间计算出的中级相关系
数，作为第a+1层虚拟代理中各配电变压器与原始所属馈线之间的高级相关系数。
[0185] 在本发明实施例的一种具体实施方式中，所述高级相关系数阈值的计算方法包括：
[0186] 将所有配电变压器与其原始所属馈线之间的高级相关系数组成指标集合；
[0187] 按照设定的区间步长将统计区间划分成若干个子区间，计算指标集合在各子区间的累计分布概率，形成累计分布概率集合；
[0188] 选择直线y＝Ax+B对累计分布概率集合的下降段拟合的平均误差最小时，拟合直线y＝Ax+B与x轴的交点作为高级相关系数阈值。
[0189] 实施例4
[0190] 一种数据驱动的配电网线变关系诊断系统，包括：存储介质和处理器；
[0191] 所述存储介质用于存储指令；
[0192] 所述处理器用于根据所述指令进行操作以执行根据实施例1中任一项所述方法的步骤。
[0193] 本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实
施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机
可用存储介质(包括但不限于磁盘存储器、CD‑ROM、光学存储器等)上实施的计算机程序产
品的形式。
[0194] 本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流
程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序
指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产
生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实
现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0195] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指
令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或
多个方框中指定的功能。
[0196] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或
其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一
个方框或多个方框中指定的功能的步骤。
[0197] 以上结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员
在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多
形式，这些均属于本发明的保护之内。
[0198] 以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本
发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变
化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其
等效物界定。
