一种在线识别失准电能表的方法
阅读说明：本技术 一种在线识别失准电能表的方法 (Method for on-line identifying misalignment electric energy meter ) 是由 夏桃芳 林华 高琛 李建新 王雅平 陈前 詹世安 詹文 鄢盛腾 丁忠安 陈吴晓 于 2021-07-07 设计创作，主要内容包括：本发明提出一种在线识别失准电能表的方法,包括以下步骤；步骤S1：选取与待校验电能表所在表箱A最邻近的表箱B,结合用电信息采集系统电压、电流曲线与基尔霍夫电压、电流定律,递推得到与待校验电能表同相位的最邻近表箱节点电压与电流幅值；步骤S2：依据用电信息采集系统相位信息,确定同箱同相电能表数量；步骤S3：进行多元线性回归分析,得到待校验电能表运行误差；步骤S4：当待校验电能表运行误差大于允许阈值,则识别为失准电能表；本发明仅需要同表箱及邻近表箱电能表一天的电压、电流曲线数据即可准确定位待校验电能表运行误差,同时线性回归原理简单,计算高效,能够迅速及时定位失准电能表,具有较强的工程实用性。(The invention provides a method for identifying a misalignment electric energy meter on line, which comprises the following steps of; step S1: selecting a meter box B which is most adjacent to a meter box A where the electric energy meter to be verified is located, combining a voltage and current curve of an electricity utilization information acquisition system, kirchhoff voltage and current law, and recurrently obtaining a node voltage and a current amplitude value of the most adjacent meter box which is in the same phase as the electric energy meter to be verified; step S2: determining the number of the same-box and same-phase electric energy meters according to the phase information of the power utilization information acquisition system; step S3: performing multiple linear regression analysis to obtain the running error of the electric energy meter to be verified; step S4: when the running error of the electric energy meter to be verified is larger than an allowable threshold value, identifying the electric energy meter as a misalignment electric energy meter; the method can accurately position the running error of the electric energy meter to be verified only by voltage and current curve data of the electric energy meter in the same meter box and the adjacent meter box in one day, has simple linear regression principle and high calculation efficiency, can quickly and timely position the misaligned electric energy meter, and has stronger engineering practicability.)
技术领域
本发明涉及电网运维
技术领域
，尤其是一种在线识别失准电能表的方法。背景技术
电能的采集、计量、存储与传输主要依靠电能表，其准确性关系到电网企业和每个用户的切身利益。为保证计量公正性，目前电网企业主要通过拆回检定、周期轮换、定期抽检与现场试验等方式对计量失准装置进行管控，需要耗费较大的人力、物力与财力，但有时仍无法及时发现并更换失准表计。因此，基于电力数据实现电能表运行误差在线监测，是实现包括现场的计量装置全寿命周期管理的重要手段。目前，在线识别失准电能表主要有以下方法：
(1)基于统计学的方法：该方法考虑电能计量层级的总-分性，结合统计学分析方法进行数据拟合，并结合拟合系数和电能表运行误差相关性，定位失准电能表。但是目前常规方法是采用电量数据进行拟合，所需样本量的持续时长需几个月，无法及时发现失准表计。
(2)基于聚类的方法：该方法认为正常的同类用户在负荷曲线形态等指标上具有相似性，为了凸显失准用户常见的用电量趋势性下降、日负荷曲线异常以及报装容量利用率偏低等异常，基于聚类的方法多采用月/日用电量、电压/电流三相不平衡度和负荷变化率与波动率等指标进行聚类，但工程实际中大量用户用电行为具有明显的离散性，聚类识别异常用户的误报率往往还高于基于分类的方法。
(3)基于状态估计的方法：该方法利用线路各节点的电压、注入功率与线路电气参数等状态量具有强耦合性的特点辨识用户上报的异常电量数据，通过同一线路上用户电气量的强耦合性检测异常。但该方法对数据准确性要求高，配电线路电气参数的准确性等因素可能影响检测结果。
上述方法虽然能够用数据方法定位部分失准电能表，但是存在工程应用上的局限性，如准确性不高、数据依赖程度较高或所需样本时间跨度较长等。
发明内容
本发明提出一种在线识别失准电能表的方法，仅需要同表箱及邻近表箱电能表一天的电压、电流曲线数据即可准确定位待校验电能表运行误差，同时线性回归原理简单，计算高效，能够迅速及时定位失准电能表，具有较强的工程实用性。
本发明采用以下技术方案。
一种在线识别失准电能表的方法，用于对低压配电线路的表箱内的高误差电能表进行识别，所述方法包括以下步骤；
步骤S1：选取与待校验电能表所在表箱A最邻近的表箱B，结合用电信息采集系统电压、电流曲线与基尔霍夫电压、电流定律，递推得到与待校验电能表同相位的最邻近表箱节点电压与电流幅值；
步骤S2：依据用电信息采集系统相位信息，确定待校验电能表所在表箱的同箱同相电能表数量；
步骤S3：结合用电信息采集系统电压、电流曲线数据和同箱同相电能表数量，带入对应公式进行多元线性回归分析，得到待校验电能表运行误差；
步骤S4：当待校验电能表运行误差大于允许阈值，则识别为失准电能表。
所述步骤S3中，若表箱A内的电能表为两只并联且同相的电能表，则分析过程为同相别递推，具体如下；
设低压配电线路的电阻为R；则按基尔霍夫电压定律，电能表节点1和电能表节点2的电压电流以公式表述为
U12＝U1+R1I1＝U2+R2I2 公式一；
表箱A和表箱B的电压电流以公式表述为
U0＝U12+R12I12＝U3+R3I3 公式二；
式中，U1,U2分别为表箱A内电能表1和电能表2电压幅值；I1,I2为对应电流幅值；R1,R2为对应线路电阻值；U12、I12为电能表1和电能表2的公共上层节点电压幅值与电流幅值，在这里认为是表箱节点A的电压幅值与电流幅值；U3、I3为表箱A最邻近表箱B的节点电压幅值与电流幅值，可结合基尔霍夫电压、电流定律，由邻近表箱内同相电能表递推得到；R12、R3为表箱A与表箱B对应上层节点的电阻值；U0为上层节点电压幅值；
由公式二推导得公式三：U3＝U12+R12I12-R3I3 公式三；
把公式一表述的电能表节点1和电能表节点2的电压电流分别代入公式三，得到以下公式
U3＝U1+R1I1+R12I12-R3I3 公式四；
U3＝U2+R2I2+R12I12-R3I3 公式五；
由基尔霍夫电流定律可知，表箱A电流满足如下关系：
I12＝I1+I2 公式六；
把公式六代入公式四、公式五，得到以下公式
U3＝U1+(R1+R12)I1+R12I2-R3I3 公式七；
U3＝U2+R12I1+(R2+R12)I2-R3I3 公式八；
假设电能表1为失准电能表，电能表示数为I′1,该电流实际值为kI′1，k为修正误差，则公式七、公式八可转换为以下公式
U3＝U1+(R1+R12)kI’1+R12I2-R3I3 公式九；
U3＝U2+kR12I’1+(R2+R12)I2-R3I3 公式十；
通过电压、电流曲线对公式九、公式十进行线性回归，构建多元线性回归方程，以公式表述为
U3，i＝β1U1，i+β2I’1，i+β3I2，i-β4I3，i 公式十一；
U3，i＝β5U2，i+β6I’1，i+β7I2，i-β8I3，i 公式十二；
式中，i为时间序号1,2,3，……96，选取曲线数据样本。在多元线性拟合回归方程中，U3,i为回归方程的因变量样本，U1,i、U2,i、I′1,i、I2,i、I3,i为回归方程自变量样本，β为对应拟合系数；
结合公式九、公式十、公式十一、公式十二，电流I′1,i和I2,i回归系数β6和β3的关系式表述为
式中，k为电能表1的修正误差；电能表1的计算运行误差如下式所示：
按上述方法计算表箱A内其它电能表的误差；当式中计算结果ε为正数时，则电能表存在正误差，即电能表多计量电量；当计算结果ε为负数时，则电能表存在负误差，即电能表少计量电量；如果计算运行误差ε超过规定的运行误差阈值，则研判为失准电能表。
所述步骤S3中，若表箱A内的电能表为并联同相的电能表，且电能表总数大于两只，则分析过程为同相别递推，具体如下；
由基尔霍夫电压定律可知，电能表节点1、2和3的电压电流以公式表述如下
U总＝U1+R1I1＝U2+R2I2＝U3+R3I3 公式十五；
表箱A和表箱B的电压电流以公式表述为
U0＝U总+R总I总＝U41+R41I41 公式十六；
式中，U1,U2,U3分别为表箱A内电能表1、2和3电压幅值；I1,I2,I3为对应电流幅值；R1,R2,R3为对应线路电阻值；U总认为是表箱节点A的电压幅值；U41、I41为表箱A最邻近表箱B的节点电压幅值与电流幅值，可结合基尔霍夫电压、电流定律，由邻近表箱内同相电能表递推得到；R总、R41为表箱A与表箱B对应上层节点的电阻值；U0为上层节点电压幅值；
由公式推导得U41＝U总+R总I总-R41I41 公式十七；
将U总＝U2+R2I2代入公式十七，得
U41＝U2+R2I2+R总I总-R41I41 公式十八；
由基尔霍夫电流定律可得，表箱A电流满足如下关系：
I总＝I1+I2+...In 公式十九；
把公式十九代入公式十八，得
U41＝U2+(R2+R总)I2+R总I1+R总(I3+I4+...In)-R41I41公式二十；
假设电能表1为失准电能表，电能表示数为I′1,该电流实际值为kI′1，k为修正误差，则公式二十可转化为
U41＝U2+(R2+R总)I2+R总kI′1+R总(I3+I4+...In)-R41I41公式二十一；
通过电压、电流曲线对公式二十一进行线性回归，构建多元线性回归方程为
U41，i＝β11U2，i+β2I2，i+β1I′1，i+β3I3，i+β4I4，i+...βnIn，i-β41I41，i公式二十二；
式中，i为时间序号1,2,3，……96，选取曲线数据样本；在多元线性拟合回归方程中，U41,i为回归方程的因变量样本，U2,i、I′1,i、I2,i、I3,i、I4,i、…In,i、I41,i为回归方程自变量样本，β为对应拟合系数；
结合公式二十一、公式二十二，电流I′1,i和I3,i回归系数β3和β4的关系以公式表述为
电能表运行误差为
按上述方法计算表箱A内其它电能表的误差；当式中计算结果ε为正数时，则电能表存在正误差，即电能表多计量电量；当计算结果ε为负数时，则电能表存在负误差，即电能表少计量电量；如果计算运行误差ε超过规定的运行误差阈值，则研判为失准电能表。
在对公式九、公式十进行线性回归分析时，选取的曲线数据样本为一天96点的电压、电流曲线数据。
在对公式二十一进行线性回归分析时，选取的曲线数据样本为一天96点的电压、电流曲线数据。
所述表箱B是多表位表箱，或是离散的单表位表箱。
在步骤S1中结合营销系统档案信息，选取与待校验电能表所在表箱A最邻近的表箱B。
本发明结合基尔霍夫电压电流定律与多元线性回归分析，由下而上递推得到含拟合系数的多个电能表间电压与电流对应关系，并经公式推导确定待校验电能表运行误差与拟合系数相关关系，从而得到电能表运行误差准确表达式，精准识别失准电能表。
本发明所述方法的优点在于：
1.该方法仅需要结合一天电压、电流曲线数据即可进行失准电能表分析计算，所需样本时间跨度小，及时迅速定位失准电能表，实现计量装置含现场的全寿命周期管理。
2.多元线性拟合算法简单，计算速度快，具有较好的工程实用性。
3.为电能表失准更换提供全新的校验方法。
附图说明
下面结合附图和
具体实施方式
对本发明进一步详细的说明：
附图1是同表箱同相电能表总数为二时的电气拓扑结构示意图；
附图2是同表箱同相电能表总数大于二时的电气拓扑结构示意图；
附图3是本发明所述方法的流程示意图。
具体实施方式
如图所示，一种在线识别失准电能表的方法，用于对低压配电线路的表箱内的高误差电能表进行识别，所述方法包括以下步骤；
步骤S1：选取与待校验电能表所在表箱A最邻近的表箱B，结合用电信息采集系统电压、电流曲线与基尔霍夫电压、电流定律，递推得到与待校验电能表同相位的最邻近表箱节点电压与电流幅值；
步骤S2：依据用电信息采集系统相位信息，确定待校验电能表所在表箱的同箱同相电能表数量；
步骤S3：结合用电信息采集系统电压、电流曲线数据和同箱同相电能表数量，带入对应公式进行多元线性回归分析，得到待校验电能表运行误差；
步骤S4：当待校验电能表运行误差大于允许阈值，则识别为失准电能表。
所述步骤S3中，若表箱A内的电能表为两只并联且同相的电能表，则分析过程为同相别递推，具体如下；
设低压配电线路的电阻为R；则按基尔霍夫电压定律，电能表节点1和电能表节点2的电压电流以公式表述为
U12＝U1+R1I1＝U2+R2I2 公式一；
表箱A和表箱B的电压电流以公式表述为
U0＝U12+R12I12＝U3+R3I3 公式二；
式中，U1,U2分别为表箱A内电能表1和电能表2电压幅值；I1,I2为对应电流幅值；R1,R2为对应线路电阻值；U12、I12为电能表1和电能表2的公共上层节点电压幅值与电流幅值，在这里认为是表箱节点A的电压幅值与电流幅值；U3、I3为表箱A最邻近表箱B的节点电压幅值与电流幅值，可结合基尔霍夫电压、电流定律，由邻近表箱内同相电能表递推得到；R12、R3为表箱A与表箱B对应上层节点的电阻值；U0为上层节点电压幅值；
由公式二推导得公式三：U3＝U12+R12I12-R3I3 公式三；
把公式一表述的电能表节点1和电能表节点2的电压电流分别代入公式三，得到以下公式
U3＝U1+R1I1+R12I12-R3I3 公式四；
U3＝U2+R2I2+R12I12-R3I3 公式五；
由基尔霍夫电流定律可知，表箱A电流满足如下关系：
I12＝I1+I2 公式六；
把公式六代入公式四、公式五，得到以下公式
U3＝U1+(R1+R12)I1+R12I2-R3I3 公式七；
U3＝U2+R12I1+(R2+R12)I2-R3I3 公式八；
假设电能表1为失准电能表，电能表示数为I′1,该电流实际值为kI′1，k为修正误差，则公式七、公式八可转换为以下公式
U3＝U1+(R1+R12)kI’1+R12I2-R3I3 公式九；
U3＝U2+kR12I’1+(R2+R12)I2-R3I3 公式十；
通过电压、电流曲线对公式九、公式十进行线性回归，构建多元线性回归方程，以公式表述为
U3，i＝β1U1，i+β2I’1，i+β3I2，i-β4I3，i 公式十一；
U3，i＝β5U2，i+β6I’1，i+β7I2，i-β8I3，i 公式十二；
式中，i为时间序号1,2,3，……96，选取曲线数据样本。在多元线性拟合回归方程中，U3,i为回归方程的因变量样本，U1,i、U2,i、I′1,i、I2,i、I3,i为回归方程自变量样本，β为对应拟合系数；
结合公式九、公式十、公式十一、公式十二，电流I′1,i和I2,i回归系数β6和β3的关系式表述为
式中，k为电能表1的修正误差；电能表1的计算运行误差如下式所示：
按上述方法计算表箱A内其它电能表的误差；当式中计算结果ε为正数时，则电能表存在正误差，即电能表多计量电量；当计算结果ε为负数时，则电能表存在负误差，即电能表少计量电量；如果计算运行误差ε超过规定的运行误差阈值，则研判为失准电能表。
所述步骤S3中，若表箱A内的电能表为并联同相的电能表，且电能表总数大于两只，则分析过程为同相别递推，具体如下；
由基尔霍夫电压定律可知，电能表节点1、2和3的电压电流以公式表述如下
U总＝U1+R1I1＝U2+R2I2＝U3+R3I3 公式十五；
表箱A和表箱B的电压电流以公式表述为
U0＝U总+R总I总＝U41+R41I41 公式十六；
式中，U1,U2,U3分别为表箱A内电能表1、2和3电压幅值；I1,I2,I3为对应电流幅值；R1,R2,R3为对应线路电阻值；U总认为是表箱节点A的电压幅值；U41、I41为表箱A最邻近表箱B的节点电压幅值与电流幅值，可结合基尔霍夫电压、电流定律，由邻近表箱内同相电能表递推得到；R总、R41为表箱A与表箱B对应上层节点的电阻值；U0为上层节点电压幅值；
由公式推导得U41＝U总+R总I总-R41I41 公式十七；
将U总＝U2+R2I2代入公式十七，得
U41＝U2+R2I2+R总I总-R41I41 公式十八；
由基尔霍夫电流定律可得，表箱A电流满足如下关系：
I总＝I1+I2+...In 公式十九；
把公式十九代入公式十八，得
U41＝U2+(R2+R总)I2+R总I1+R总(I3+I4+...In)-R41I41公式二十；
假设电能表1为失准电能表，电能表示数为I′1,该电流实际值为kI′1，k为修正误差，则公式二十可转化为
U41＝U2+(R2+R总)I2+R总kI′1+R总(I3+I4+...In)-R41I41公式二十一；
通过电压、电流曲线对公式二十一进行线性回归，构建多元线性回归方程为
U41，i＝β11U2，i+β2I2，i+β1I′1，i+β3I3，i+β4I4，i+...βnIn，i-β41I41，i公式二十二；
式中，i为时间序号1,2,3，……96，选取曲线数据样本；在多元线性拟合回归方程中，U41,i为回归方程的因变量样本，U2,i、I′1,i、I2,i、I3,i、I4,i、…In,i、I41,i为回归方程自变量样本，β为对应拟合系数；
结合公式二十一、公式二十二，电流I′1,i和I3,i回归系数β3和β4的关系以公式表述为
电能表运行误差为
按上述方法计算表箱A内其它电能表的误差；当式中计算结果ε为正数时，则电能表存在正误差，即电能表多计量电量；当计算结果ε为负数时，则电能表存在负误差，即电能表少计量电量；如果计算运行误差ε超过规定的运行误差阈值，则研判为失准电能表。
在对公式九、公式十进行线性回归分析时，选取的曲线数据样本为一天96点的电压、电流曲线数据。
在对公式二十一进行线性回归分析时，选取的曲线数据样本为一天96点的电压、电流曲线数据。
所述表箱B是多表位表箱，或是离散的单表位表箱。
在步骤S1中结合营销系统档案信息，来选取与待校验电能表所在表箱A最邻近的表箱B。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：基于数字示波器的波形展示方法和数字示波器