CN110275129B - 一种确定高压电能计量装置合成误差的方法和系统 - Google Patents

一种确定高压电能计量装置合成误差的方法和系统 Download PDF

Info

Publication number
CN110275129B
CN110275129B CN201910526765.1A CN201910526765A CN110275129B CN 110275129 B CN110275129 B CN 110275129B CN 201910526765 A CN201910526765 A CN 201910526765A CN 110275129 B CN110275129 B CN 110275129B
Authority
CN
China
Prior art keywords
transformer
error
voltage
electric energy
difference
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910526765.1A
Other languages
English (en)
Other versions
CN110275129A (zh
Inventor
李贺龙
于海波
林繁涛
袁晓蕾
赵进全
王春雨
刘佳
王兴媛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Hubei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910526765.1A priority Critical patent/CN110275129B/zh
Publication of CN110275129A publication Critical patent/CN110275129A/zh
Application granted granted Critical
Publication of CN110275129B publication Critical patent/CN110275129B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/04Testing or calibrating of apparatus covered by the other groups of this subclass of instruments for measuring time integral of power or current

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明提供一种确定高压电能计量装置合成误差的方法和系统。所述方法和系统将高压电能计量装置等效为多输入单输出系统,通过建立互感器误差的多参数降维模型,将互感器历史检定数据进行预处理后,采用支持向量机回归算法求解多参数降维模型,确定互感器实际运行工况下的比差和角差,再与在线监测获得的二次回路误差,电能表误差进行合成,最后确定高压电能计量装置的合成误差。所述方法和系统实现了对难以直接在线监测的电压互感器和电流互感器的动态计量误差的准确评估,实时性好、精度高,再结合在线监测的二次回路、电能表的误差进行电能计量装置的合成误差计算,从而有效地评价实际运行工况下电能计量装置的整体动态误差水平。

Description

一种确定高压电能计量装置合成误差的方法和系统
技术领域
本发明涉及电能计量领域,并且更具体地,涉及一种确定高压电能计量装置合成误差的方法和系统。
背景技术
电能计量是电网经济核算的依据,其准确度直接关系到供、受电双方的经济效益。电能计量装置的误差由电压互感器误差、电流互感器误差、二次回路压降误差以及电能表误差四部分综合而成,只有明确综合误差才能合理评价电能计量的准确程度,由于误差来源因素的多样性和内在关联性,电能计量整体误差呈现动态变化的特征。传统的电能计量装置误差合成方法是直接对电压互感器、电流互感器、二次回路和电能表的误差进行综合计算,但计量装置各部分的误差与实际的电压、电流、二次负荷等相关,传统方法虽然计算简单但无法评价电能计量装置实际运行工况下的综合误差;因此如何确定一种实时性强、精度高的电能计量装置误差合成方法和系统就成为一个亟需解决的问题。
发明内容
为了解决现有技术中电能计量装置合成误差方法无法评价实际运行工况下的综合误差的技术问题,本发明提供一种确定高压电能计量装置合成误差的方法,所述方法包括:
基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数,其中,所述互感器包括电压互感器和电流互感器;
根据互感器实际运行工况下的测量数据,基于所述最优回归函数确定互感器比差和角差;
基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差。
进一步地,所述方法采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数之前还包括将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程,其中:
所述互感器的测量误差ε表达式为:
ε=f(fPTPT,fCTCT)
取互感器一次电压、电流和二次负荷作为互感器误差的影响因子,则电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT的表达式为:
Figure BDA0002098440710000021
式中,V(t)为电压互感器一次电压,I(t)为电流互感器一次电流,
S2U(t)为电压互感器二次负荷,S2I(t)为电流互感器二次负荷;
所述确定互感器测量误差的多参数降维方程为:
a=Xb
式中,a=(fPTPT,fCTCT),b=(V(t),I(t),S2U(t),S2I(t)),X为所述多参数方程的解。
进一步地,所述方法在将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程之前还包括对互感器的历史检定数据进行预处理,其中:
二次负荷与互感器误差呈线性关系,不同负荷下的误差按照一阶线性插值法扩充,假设第i点和第j点的二次负荷为S2(i)和S2(j),数据量增加至N倍,则两个负荷点之间的插值序列为:
Figure BDA0002098440710000031
一次电压或电流与互感器误差的关系与励磁特性有关,故不同电压或电流下的互感器测量误差按照分段三次Hermite插值进行数据扩充,假设第i点和第j点一次电压为V1(i)和V1(j),两个节点的一阶导数为fi′和fj′,则采用三次Hermite多项式H3(x)作为插值函数:
H3(x)=V1(i)αi(x)+V1(j)αj(x)+fi′βi(x)+fj′βj(x)
式中,
Figure BDA0002098440710000032
Figure BDA0002098440710000033
进一步地,所述基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数包括:
设含有l个训练样本的训练集为{(xi,yi),i=1,2,…l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,
Figure BDA0002098440710000034
yi∈R为对应的输出值,设在高维特征空间中建立的线性回归函数为
f(x)=wΦ(x)+b
其中,Φ(x)为非线性映射函数,w表示权系数相量,b表示阈值;
通过训练集找到w和b的最优解,使得|f(x)-wΦ(x)-b|≤ε,则可根据回归函数得到任意输入下的输出值;
定义ε线性不敏感损失函数,其表达式为:
Figure BDA0002098440710000035
其中,f(x)为回归函数返回的预测值;y为对应的真实值,即表示若f(x)与y之间的差别小于等于ε,则损失等于0;
引进松弛变量
Figure BDA0002098440710000041
则通过训练集寻找w和b的最优解的表达式为:
Figure BDA0002098440710000042
其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大;ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小;
引入Largrange函数,将寻找w和b的最优解的表达式转换为对偶格式,所述对偶格式的表达式为:
Figure BDA0002098440710000043
其中,K(xi,xj)=Φ(xi)Φ(xj)为核函数;
当求解对偶格式的表达式得到的最优解为α=[α12,…,αl],
Figure BDA0002098440710000044
则有
Figure BDA0002098440710000045
Figure BDA0002098440710000046
其中,Nnsv为支持向量个数;
根据确定的w*和b*可知多参数降维方程的最优回归函数表达式为:
Figure BDA0002098440710000047
进一步地,所述基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差,其计算公式为:
Figure BDA0002098440710000051
式中,εall为高压电能计量装置的合成误差,fIa、fIb、fIc、δIa、δIb、δIc为三相电流互感器比差和角差,fUa、fUb、fUc、δUa、δUb、δUc为三相电压互感器比差和角差,fa、fb、fc为二次回路压降误差,εe为电能表误差,
Figure BDA0002098440710000052
为负载功率因数角。
根据本发明的另一方面,本发明提供一种确定高压电能计量装置合成误差的系统,所述系统包括:
数据采集单元,其用于采集高压输电线路的互感器测量数据;
最优函数单元,其用于基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数,其中,所述互感器包括电压互感器和电流互感器;
互感器误差单元,其用于根据互感器实际运行工况下的测量数据,基于所述最优回归函数确定互感器比差和角差;
合成误差单元,其用于基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差。
进一步地,所述系统还包括方程设置单元,其用于将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程,其中:
所述互感器的测量误差ε表达式为:
ε=f(fPTPT,fCTCT)
取互感器一次电压、电流和二次负荷作为互感器误差的影响因子,则电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT的表达式为:
Figure BDA0002098440710000061
式中,V(t)为电压互感器一次电压,I(t)为电流互感器一次电流,S2U(t)为电压互感器二次负荷,S2I(t)为电流互感器二次负荷;
所述确定互感器测量误差的多参数降维方程为:
a=Xb
式中,a=(fPTPT,fCTCT),b=(V(t),I(t),S2U(t),S2I(t)),X为所述多参数方程的解。
进一步地,所述系统还包括数据预处理单元,其用于对互感器的历史检定数据进行预处理,其中:
二次负荷与互感器误差呈线性关系,不同负荷下的误差按照一阶线性插值法扩充,假设第i点和第j点的二次负荷为S2(i)和S2(j),数据量增加至N倍,则两个负荷点之间的插值序列为:
Figure BDA0002098440710000062
一次电压或电流与互感器误差的关系与励磁特性有关,故不同电压或电流下的互感器测量误差按照分段三次Hermite插值进行数据扩充,假设第i点和第j点一次电压为V1(i)和V1(j),两个节点的一阶导数为fi′和fj′,则采用三次Hermite多项式H3(x)作为插值函数:
H3(x)=V1(i)αi(x)+V1(j)αj(x)+fi′βi(x)+fj′βj(x)
式中,
Figure BDA0002098440710000063
Figure BDA0002098440710000064
进一步地,所述最优函数单元基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数包括:
设含有l个训练样本的训练集为{(xi,yi),i=1,2,…l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,
Figure BDA0002098440710000071
yi∈R为对应的输出值,设在高维特征空间中建立的线性回归函数为
f(x)=wΦ(x)+b
其中,Φ(x)为非线性映射函数,w表示权系数相量,b表示阈值,通过训练集找到w和b的最优解,使得|f(x)-wΦ(x)-b|≤ε,则可根据回归函数得到任意输入下的输出值;
定义ε线性不敏感损失函数,其表达式为:
Figure BDA0002098440710000072
其中,f(x)为回归函数返回的预测值;y为对应的真实值,即表示若f(x)与y之间的差别小于等于ε,则损失等于0;
引进松弛变量
Figure BDA0002098440710000073
则通过训练集寻找w和b的最优解的表达式为:
Figure BDA0002098440710000074
其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大;ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小;
引入Largrange函数,将寻找w和b的最优解的表达式转换为对偶格式,所述对偶格式的表达式为:
Figure BDA0002098440710000075
其中,K(xi,xj)=Φ(xi)Φ(xj)为核函数;
当求解对偶格式的表达式得到的最优解为α=[α12,…,αl],
Figure BDA0002098440710000081
则有
Figure BDA0002098440710000082
Figure BDA0002098440710000083
其中,Nnsv为支持向量个数;
根据确定的w*和b*可知多参数降维方程的最优回归函数表达式为:
Figure BDA0002098440710000084
进一步地,所述合成误差单元基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差,其计算公式为:
Figure BDA0002098440710000085
式中,εall为高压电能计量装置的合成误差,fIa、fIb、fIc、δIa、δIb、δIc为三相电流互感器比差和角差,fUa、fUb、fUc、δUa、δUb、δUc为三相电压互感器比差和角差,fa、fb、fc为二次回路压降误差,εe为电能表误差,
Figure BDA0002098440710000086
为负载功率因数角。
本发明技术方案提供的确定高压电能计量装置合成误差的方法和系统将高压电能计量装置等效为多输入单输出系统,通过建立互感器误差的多参数降维模型,将互感器历史检定数据进行预处理后,采用支持向量机回归算法求解多参数降维模型,确定互感器实际运行工况下的比差和角差,再与在线监测获得的二次回路误差,电能表误差进行合成,最后确定高压电能计量装置的合成误差。所述确定高压电能计量装置合成误差的方法和系统通过检测实际运行工况下的二次电压、电流、二次负荷得到电能计量装置实际运行工况下的动态计量误差,实现了对难以直接在线监测的电压互感器和电流互感器的动态计量误差的准确评估,再结合在线监测的二次回路、电能表的误差进行电能计量装置的合成误差计算,从而有效地评价实际运行工况下电能计量装置的整体动态误差水平,可靠评价电能计量的准确度,该方法实时性好、精度高而且无需建立复杂模型。而且本发明通过对历史检定数据进行回归计算,提高了模型的抗扰性,且算法无需大量样本数据,利用小样本即可获得全局最优解。
附图说明
通过参考下面的附图,可以更为完整地理解本发明的示例性实施方式:
图1为根据本发明优选实施方式的确定高压电能计量装置合成误差的方法的流程图;
图2为根据本发明优选实施方式的确定高压电能计量装置合成误差的系统的结构示意图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
图1为根据本发明优选实施方式的确定高压电能计量装置合成误差的方法的流程图。如图1所示,本优选实施方式所述的确定高压电能计量装置合成误差的方法100从步骤101开始。
在步骤101,将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程。
在本优选实施方式中,所述互感器的测量误差ε表达式为:
ε=f(fPTPT,fCTCT)
互感器的比差和角差与电压、电流、二次负荷、频率、温度等均有关联,为简化模型,本优选实施方式只选取电压、电流、二次负荷作为互感器误差的影响因素,则电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT的表达式为:
Figure BDA0002098440710000101
式中,V(t)为电压互感器一次电压,I(t)为电流互感器一次电流,S2U(t)为电压互感器二次负荷,S2I(t)为电流互感器二次负荷;
所述确定互感器测量误差的多参数降维方程为:
a=Xb
式中,a=(fPTPT,fCTCT),b=(V(t),I(t),S2U(t),S2I(t)),X为所述多参数方程的解。
在步骤102,对互感器的历史检定数据进行预处理。在本优选实施方式中,由于只选取电压、电流、二次负荷作为互感器误差的影响因素,因此历史检定数据的预处理主要是针对一次电压、一次电流和二次负荷。
二次负荷与互感器误差呈线性关系,不同负荷下的误差按照一阶线性插值法扩充,假设第i点和第j点的二次负荷为S2(i)和S2(j),数据量增加至N倍,则两个负荷点之间的插值序列为:
Figure BDA0002098440710000102
一次电压或电流与互感器误差的关系与励磁特性有关,故不同电压或电流下的互感器测量误差按照分段三次Hermite插值进行数据扩充,假设第i点和第j点一次电压为V1(i)和V1(j),两个节点的一阶导数为fi′和fj′,则采用三次Hermite多项式H3(x)作为插值函数:
H3(x)=V1(i)αi(x)+V1(j)αj(x)+fi′βi(x)+fj′βj(x)
式中,
Figure BDA0002098440710000111
Figure BDA0002098440710000112
在步骤103,基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数,其中,所述互感器包括电压互感器和电流互感器。支持向量机回归算法对多参数降维方程求解的核心是对预处理的数据进行回归拟合分析时,寻找一个最优函数,使得所有的数据样本离该最优函数的误差最小。
在步骤104,根据互感器实际运行工况下的测量数据,基于所述最优回归函数确定互感器比差和角差。在本发明中,通过确定的最优回归函数,将互感器实际运行工况下的一次电压、一次电流和二次负荷代入,即可得到动态的互感器计量误差。
在步骤105,基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差。
优选地,所述基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数包括:
设含有l个训练样本的训练集为{(xi,yi),i=1,2,…l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,
Figure BDA0002098440710000113
yi∈R为对应的输出值,设在高维特征空间中建立的线性回归函数为
f(x)=wΦ(x)+b
其中,Φ(x)为非线性映射函数,w表示权系数相量,b表示阈值;
通过训练集找到w和b的最优解,使得|f(x)-wΦ(x)-b|≤ε,则可根据回归函数得到任意输入下的输出值;
定义ε线性不敏感损失函数,其表达式为:
Figure BDA0002098440710000121
其中,f(x)为回归函数返回的预测值;y为对应的真实值,即表示若f(x)与y之间的差别小于等于ε,则损失等于0;
引进松弛变量
Figure BDA0002098440710000122
则通过训练集寻找w和b的最优解的表达式为:
Figure BDA0002098440710000123
其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大;ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小;
引入Largrange函数,将寻找w和b的最优解的表达式转换为对偶格式,所述对偶格式的表达式为:
Figure BDA0002098440710000124
其中,K(xi,xj)=Φ(xi)Φ(xj)为核函数;
当求解对偶格式的表达式得到的最优解为α=[α12,…,αl],
Figure BDA0002098440710000125
则有
Figure BDA0002098440710000126
Figure BDA0002098440710000127
其中,Nnsv为支持向量个数;
根据确定的w*和b*可知多参数降维方程的最优回归函数表达式为:
Figure BDA0002098440710000131
优选地,所述基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差,其计算公式为:
Figure BDA0002098440710000132
式中,εall为高压电能计量装置的合成误差,fIa、fIb、fIc、δIa、δIb、δIc为三相电流互感器比差和角差,fUa、fUb、fUc、δUa、δUb、δUc为三相电压互感器比差和角差,fa、fb、fc为二次回路压降误差,εe为电能表误差,
Figure BDA0002098440710000133
为负载功率因数角。
在Matlab软件中搭建电能计量全系统模型,对本发明进行了仿真验证,得到电能计量装置合成误差的绝对误差不超过0.0015%,可见本方法计算精确。所述仿真验证的具体步骤为:
步骤1:在Matlab/Simulink中建立500kV输电线路及电能计量全系统的仿真模型,为模拟现场实际检定情况,对电压互感器80%、100%、105%电压点,100%,25%,5%二次负荷点的误差进行测量,对电流互感器1%、5%、20%、100%、120%电流点,100%,25%,5%二次负荷点的误差进行测量,作为历史检定数据。
步骤2:对原始数据预处理、归一化后进行SVR求解多参数降维方程,SVR核函数选择RBF核函数,利用交叉验证方法寻找最佳参数c(惩罚因子)和参数g(RBF核函数中的方差),然后利用最佳的参数训练模型。
步骤3:互感器误差在线评估:为验证SVR回归的准确性,随机设置几组不同的一次电压和二次负荷进行测试,得到电压互感器的比差和角差,其测试结果与真实结果比较接近。电流互感器的测试结果同电压互感器相似。
步骤4:电能计量装置误差合成:模拟实际输电线路电压波动、电流波动和互感器二次负荷波动的情况,对电能计量装置的误差进行合成计算,得到合成误差,所述合成误差与实际误差相比,绝对误差不超过0.0015%,相对误差不超过0.31%。
图2为根据本发明优选实施方式的确定高压电能计量装置合成误差的系统的结构示意图。在本优选实施方式中,所述确定高压电能计量装置合成误差的系统200包括:
方程设置单元201,其用于将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程,其中:
所述互感器的测量误差ε表达式为:
ε=f(fPTPT,fCTCT)
取互感器一次电压、电流和二次负荷作为互感器误差的影响因子,则电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT的表达式为:
Figure BDA0002098440710000141
式中,V(t)为电压互感器一次电压,I(t)为电流互感器一次电流,S2U(t)为电压互感器二次负荷,S2I(t)为电流互感器二次负荷;
所述确定互感器测量误差的多参数降维方程为:
a=Xb
式中,a=(fPTPT,fCTCT),b=(V(t),I(t),S2U(t),S2I(t)),X为所述多参数方程的解。
数据预处理单元202,其用于对互感器的历史检定数据进行预处理,其中:
二次负荷与互感器误差呈线性关系,不同负荷下的误差按照一阶线性插值法扩充,假设第i点和第j点的二次负荷为S2(i)和S2(j),数据量增加至N倍,则两个负荷点之间的插值序列为:
Figure BDA0002098440710000151
一次电压或电流与互感器误差的关系与励磁特性有关,故不同电压或电流下的互感器测量误差按照分段三次Hermite插值进行数据扩充,假设第i点和第j点一次电压为V1(i)和V1(j),两个节点的一阶导数为fi′和fj′,则采用三次Hermite多项式H3(x)作为插值函数:
H3(x)=V1(i)αi(x)+V1(j)αj(x)+fi′βi(x)+fj′βj(x)
式中,
Figure BDA0002098440710000152
Figure BDA0002098440710000153
数据采集单元203,其用于采集高压输电线路的互感器测量数据;
最优函数单元204,其用于基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数,其中,所述互感器包括电压互感器和电流互感器;
互感器误差单元205,其用于根据互感器实际运行工况下的测量数据,基于所述最优回归函数确定互感器比差和角差;
合成误差单元206,其用于基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差。
优选地,所述最优函数单元204基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数包括:
设含有l个训练样本的训练集为{(xi,yi),i=1,2,…l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,
Figure BDA0002098440710000154
yi∈R为对应的输出值,设在高维特征空间中建立的线性回归函数为
f(x)=wΦ(x)+b
其中,Φ(x)为非线性映射函数,w表示权系数相量,b表示阈值,通过训练集找到w和b的最优解,使得|f(x)-wΦ(x)-b|≤ε,则可根据回归函数得到任意输入下的输出值;
定义ε线性不敏感损失函数,其表达式为:
Figure BDA0002098440710000161
其中,f(x)为回归函数返回的预测值;y为对应的真实值,即表示若f(x)与y之间的差别小于等于ε,则损失等于0;
引进松弛变量
Figure BDA0002098440710000162
则通过训练集寻找w和b的最优解的表达式为:
Figure BDA0002098440710000163
其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大;ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小;
引入Largrange函数,将寻找w和b的最优解的表达式转换为对偶格式,所述对偶格式的表达式为:
Figure BDA0002098440710000164
其中,K(xi,xj)=Φ(xi)Φ(xj)为核函数;
当求解对偶格式的表达式得到的最优解为α=[α12,…,αl],
Figure BDA0002098440710000165
则有
Figure BDA0002098440710000171
Figure BDA0002098440710000172
其中,Nnsv为支持向量个数;
根据确定的w*和b*可知多参数降维方程的最优回归函数表达式为:
Figure BDA0002098440710000173
优选地,所述合成误差单元206基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差,其计算公式为:
Figure BDA0002098440710000174
式中,εall为高压电能计量装置的合成误差,fIa、fIb、fIc、δIa、δIb、δIc为三相电流互感器比差和角差,fUa、fUb、fUc、δUa、δUb、δUc为三相电压互感器比差和角差,fa、fb、fc为二次回路压降误差,εe为电能表误差,
Figure BDA0002098440710000175
为负载功率因数角。
已经通过参考少量实施方式描述了本发明。然而,本领域技术人员所公知的,正如附带的专利权利要求所限定的,除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
通常地,在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释,除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例,除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行,除非明确地说明。

Claims (6)

1.一种确定高压电能计量装置合成误差的方法,其特征在于,所述方法包括:
将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程,其中:
所述互感器的测量误差ε表达式为:
ε=f(fPTPT,fCTCT)
取互感器一次电压、电流和二次负荷作为互感器误差的影响因子,则电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT的表达式为:
Figure FDA0003804419500000011
式中,V(t)为电压互感器一次电压,I(t)为电流互感器一次电流,S2U(t)为电压互感器二次负荷,S2I(t)为电流互感器二次负荷;
确定互感器测量误差的多参数降维方程为:
a=Xb
式中,a=(fPTPT,fCTCT),b=(V(t),I(t),S2U(t),S2I(t)),X为所述多参数降维方程的解;
基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数,其中,所述互感器包括电压互感器和电流互感器;
根据互感器实际运行工况下的测量数据,基于所述最优回归函数确定互感器比差和角差;
基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差,其计算公式为:
Figure FDA0003804419500000021
式中,εall为高压电能计量装置的合成误差,fIa、fIb、fIc、fIa、fIb、δIc为三相电流互感器比差和角差,fUa、fUb、fUc、δUa、δUb、δUc为三相电压互感器比差和角差,fa、fb、fc为二次回路压降比差,δa、δb、δc为二次回路压降角差,εe为电能表误差,
Figure FDA0003804419500000022
为负载功率因数角。
2.根据权利要求1所述的方法,其特征在于,所述方法在将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程之前还包括对互感器的历史检定数据进行预处理,其中:
二次负荷与互感器误差呈线性关系,不同负荷下的误差按照一阶线性插值法扩充,假设第i点和第j点的二次负荷为S2(i)和S2(j),数据量增加至N倍,则两个负荷点之间的插值序列为:
Figure FDA0003804419500000023
一次电压或电流与互感器误差的关系与励磁特性有关,故不同电压或电流下的互感器测量误差按照分段三次Hermite插值进行数据扩充,假设第i点和第j点一次电压为V1(i)和V1(j),两个节点的一阶导数为fi′和fj′,则采用三次Hermite多项式H3(x)作为插值函数:
H3(x)=V1(i)αi(x)+V1(j)αj(x)+fi′βi(x)+fj′βj(x)
式中,
Figure FDA0003804419500000024
Figure FDA0003804419500000025
3.根据权利要求1所述的方法,其特征在于,所述基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数包括:
设含有l个训练样本的训练集为{(xi,yi),i=1,2,…l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,
Figure FDA0003804419500000031
yi∈R为对应的输出值,设在高维特征空间中建立的线性回归函数为
f(x)=wΦ(x)+b
其中,Φ(x)为非线性映射函数,w表示权系数相量,b表示阈值;
通过训练集找到w和b的最优解,使得|f(x)-wΦ(x)-b|≤ε,则可根据回归函数得到任意输入下的输出值;
定义ε线性不敏感损失函数,其表达式为:
Figure FDA0003804419500000032
其中,f(x)为回归函数返回的预测值;y为对应的真实值,即表示若f(x)与y之间的差别小于等于ε,则损失等于0;
引进松弛变量ξi,
Figure FDA0003804419500000033
则通过训练集寻找w和b的最优解的表达式为:
Figure FDA0003804419500000034
其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大;ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小;
引入Largrange函数,将寻找w和b的最优解的表达式转换为对偶格式,所述对偶格式的表达式为:
Figure FDA0003804419500000041
其中,K(xi,xj)=Φ(xi)Φ(xj)为核函数;
当求解对偶格式的表达式得到的最优解为α=[α12,…,αl],
Figure FDA0003804419500000042
则有
Figure FDA0003804419500000043
Figure FDA0003804419500000044
其中,Nnsv为支持向量个数;
根据确定的w*和b*可知多参数降维方程的最优回归函数表达式为:
Figure FDA0003804419500000045
4.一种确定高压电能计量装置合成误差的系统,其特征在于,所述系统包括:
数据采集单元,其用于采集高压输电线路的互感器测量数据;
方程设置单元,其用于将互感器的测量误差ε等价为电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT,并基于所述电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT建立确定互感器测量误差的多参数降维方程,其中:
所述互感器的测量误差ε表达式为:
ε=f(fPTPT,fCTCT)
取互感器一次电压、电流和二次负荷作为互感器误差的影响因子,则电压互感器比差fPT和角差δPT,以及电流互感器比差fCT和角差δCT的表达式为:
Figure FDA0003804419500000051
式中,V(t)为电压互感器一次电压,I(t)为电流互感器一次电流,S2U(t)为电压互感器二次负荷,S2I(t)为电流互感器二次负荷;
确定互感器测量误差的多参数降维方程为:
a=Xb
式中,a=(fPTPT,fCTCT),b=(V(t),I(t),S2U(t),S2I(t)),X为所述多参数降维方程的解;
最优函数单元,其用于基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定预先设置的多参数降维方程的最优回归函数,其中,所述互感器包括电压互感器和电流互感器;
互感器误差单元,其用于根据互感器实际运行工况下的测量数据,基于所述最优回归函数确定互感器比差和角差;
合成误差单元,其用于基于计算确定的互感器比差和角差,并结合采集的二次回路压降误差和电能表误差确定高压电能计量装置的合成误差,其计算公式为:
Figure FDA0003804419500000052
式中,εall为高压电能计量装置的合成误差,fIa、fIb、fIc、δIa、δIb、δIc为三相电流互感器比差和角差,fUa、fUb、fUc、δUa、δUb、δUc为三相电压互感器比差和角差,fa、fb、fc为二次回路压降比差,δa、δb、δc为二次回路压降角差,εe为电能表误差,
Figure FDA0003804419500000053
为负载功率因数角。
5.根据权利要求4所述的系统,其特征在于,所述系统还包括数据预处理单元,其用于对互感器的历史检定数据进行预处理,其中:
二次负荷与互感器误差呈线性关系,不同负荷下的误差按照一阶线性插值法扩充,假设第i点和第j点的二次负荷为S2(i)和S2(j),数据量增加至N倍,则两个负荷点之间的插值序列为:
Figure FDA0003804419500000061
一次电压或电流与互感器误差的关系与励磁特性有关,故不同电压或电流下的互感器测量误差按照分段三次Hermite插值进行数据扩充,假设第i点和第j点一次电压为V1(i)和V1(j),两个节点的一阶导数为fi′和fj′,则采用三次Hermite多项式H3(x)作为插值函数:
H3(x)=V1(i)αi(x)+V1(j)αj(x)+fi′βi(x)+fj′βj(x)
式中,
Figure FDA0003804419500000062
Figure FDA0003804419500000063
6.根据权利要求4所述的系统,其特征在于,所述最优函数单元基于经过预处理的高压输电线路互感器历史检定数据,采用支持向量机回归算法确定所述预先设置的多参数降维方程的最优回归函数包括:
设含有l个训练样本的训练集为{(xi,yi),i=1,2,…l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,
Figure FDA0003804419500000064
yi∈R为对应的输出值,设在高维特征空间中建立的线性回归函数为
f(x)=wΦ(x)+b
其中,Φ(x)为非线性映射函数,w表示权系数相量,b表示阈值,通过训练集找到w和b的最优解,使得|f(x)-wΦ(x)-b|≤ε,则可根据回归函数得到任意输入下的输出值;
定义ε线性不敏感损失函数,其表达式为:
Figure FDA0003804419500000065
其中,f(x)为回归函数返回的预测值;y为对应的真实值,即表示若f(x)与y之间的差别小于等于ε,则损失等于0;
引进松弛变量ξi,
Figure FDA0003804419500000071
则通过训练集寻找w和b的最优解的表达式为:
Figure FDA0003804419500000072
其中,C为惩罚因子,C越大表示对训练误差大于ε的样本惩罚越大;ε规定了回归函数的误差要求,ε越小表示回归函数的误差越小;
引入Largrange函数,将寻找w和b的最优解的表达式转换为对偶格式,所述对偶格式的表达式为:
Figure FDA0003804419500000073
其中,K(xi,xj)=Φ(xi)Φ(xj)为核函数;
当求解对偶格式的表达式得到的最优解为α=[α12,…,αl],
Figure FDA0003804419500000074
则有
Figure FDA0003804419500000075
Figure FDA0003804419500000076
其中,Nnsv为支持向量个数;
根据确定的w*和b*可知多参数降维方程的最优回归函数表达式为:
Figure FDA0003804419500000081
CN201910526765.1A 2019-06-18 2019-06-18 一种确定高压电能计量装置合成误差的方法和系统 Active CN110275129B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910526765.1A CN110275129B (zh) 2019-06-18 2019-06-18 一种确定高压电能计量装置合成误差的方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910526765.1A CN110275129B (zh) 2019-06-18 2019-06-18 一种确定高压电能计量装置合成误差的方法和系统

Publications (2)

Publication Number Publication Date
CN110275129A CN110275129A (zh) 2019-09-24
CN110275129B true CN110275129B (zh) 2022-10-04

Family

ID=67961033

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910526765.1A Active CN110275129B (zh) 2019-06-18 2019-06-18 一种确定高压电能计量装置合成误差的方法和系统

Country Status (1)

Country Link
CN (1) CN110275129B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113820632A (zh) * 2020-06-18 2021-12-21 杭州普安科技有限公司 一种无源电力二次回路接线检查系统及方法
CN113156359B (zh) * 2021-04-16 2024-01-26 中国电力科学研究院有限公司 一种用于确定电容式电压互感器计量误差的方法及系统
CN113945881B (zh) * 2021-07-27 2023-09-22 国网浙江省电力有限公司营销服务中心 计量表箱内电能表运行误差监测装置
CN114290938B (zh) * 2021-12-29 2024-01-30 国网新疆电力有限公司经济技术研究院 直流充电桩电流误差补偿方法、系统、设备及可读存储介质
CN115524657B (zh) * 2022-11-24 2023-03-10 中国电力科学研究院有限公司 一种电能表的计量误差分析方法、计量装置和自校准方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101074987A (zh) * 2007-07-20 2007-11-21 太原市优特奥科电子科技有限公司 高压电能计量装置综合误差实时在线监测方法及监测设备
CN103886486A (zh) * 2014-03-21 2014-06-25 吉首大学 一种基于支持向量机svm的电子商务推荐方法
CN106446940A (zh) * 2016-09-13 2017-02-22 大连理工大学 一种基于支持向量机的超级电容器电容值退化趋势的预测方法
CN108805217A (zh) * 2018-06-20 2018-11-13 山东大学 一种基于支持向量机的锂离子电池健康状态估计方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2817086B1 (fr) * 2000-11-17 2003-01-03 Alstom Procede d'equilibrage statistique de mesures de transformateurs de courant pour une protection differentielle de jeu de barres

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101074987A (zh) * 2007-07-20 2007-11-21 太原市优特奥科电子科技有限公司 高压电能计量装置综合误差实时在线监测方法及监测设备
CN103886486A (zh) * 2014-03-21 2014-06-25 吉首大学 一种基于支持向量机svm的电子商务推荐方法
CN106446940A (zh) * 2016-09-13 2017-02-22 大连理工大学 一种基于支持向量机的超级电容器电容值退化趋势的预测方法
CN108805217A (zh) * 2018-06-20 2018-11-13 山东大学 一种基于支持向量机的锂离子电池健康状态估计方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于CA-FSVR方法的电流互感器误差试验参数预测;陈宏 等;《三峡大学学报(自然科学版)》;20161231;第38卷(第6期);92-95 *

Also Published As

Publication number Publication date
CN110275129A (zh) 2019-09-24

Similar Documents

Publication Publication Date Title
CN110275129B (zh) 一种确定高压电能计量装置合成误差的方法和系统
CN111474510A (zh) 一种非平稳输出的电压互感器的误差评估方法及系统
CN106980910B (zh) 中长期电力负荷测算系统及方法
CN105044649B (zh) 一种电流互感器误差水平动态检测方法及系统
CN110031788A (zh) 一种空心线圈电流互感器误差环境相关性分析方法
EP2024755B1 (en) A method for determining the linear electrical response of a transformer, generator or electrical motor
Zhang et al. Detection of the ratio error drift in CVT considering AVC
Stenbakken et al. Developing linear error models for analog devices
Li et al. Self-detecting the measurement error of electronic voltage transformer based on principal component analysis-wavelet packet decomposition
CN101107533A (zh) 表征电气部件线性性能的方法和设备
CN115062681A (zh) 在电网中确定相连接
CN116736212B (zh) 一种互感器运行误差在线监测及计量方法
Zhang et al. Research on data-driven self-diagnosis for measurement errors in capacitor voltage Transformers
CN103474993B (zh) 基于映射弹性势能的电网有功承载能力定量分析指标
CN117686102A (zh) 变压器异常温度主动预警的方法及装置
Melo et al. Voltage stability assessment using modal analysis based on power systems state estimation
Asprou et al. The effect of instrument transformer accuracy class on the WLS state estimator accuracy
CN112098066A (zh) 基于门控循环单元的高压并联电抗器故障诊断方法及系统
Mazadi et al. Instantaneous voltage estimation for assessment and monitoring of flicker indices in power systems
CN115951292A (zh) 一种电子式电压互感器误差状态在线评估方法及装置
CN116226620A (zh) 一种基于svm和lstm的次同步振荡模态辨识方法
CN116125358A (zh) 电容式电压互感器介损测量误差修正方法、装置
CN112557819A (zh) 一种主动配电网故障识别方法、系统及计算机可读储存介质
Okon et al. Influence of UPFC device on power system state estimation
CN111864741A (zh) 一种线路参数误差对功率分布影响的定量分析方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant