CN109444791B - 一种电容式电压互感器误差状态评估方法和系统 - Google Patents

一种电容式电压互感器误差状态评估方法和系统 Download PDF

Info

Publication number
CN109444791B
CN109444791B CN201811332075.4A CN201811332075A CN109444791B CN 109444791 B CN109444791 B CN 109444791B CN 201811332075 A CN201811332075 A CN 201811332075A CN 109444791 B CN109444791 B CN 109444791B
Authority
CN
China
Prior art keywords
matrix
data
phase
voltage transformer
data matrix
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
CN201811332075.4A
Other languages
English (en)
Other versions
CN109444791A (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.)
WUHAN GLORY ROAD INTELLIGENT TECHNOLOGY Co.,Ltd.
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201811332075.4A priority Critical patent/CN109444791B/zh
Publication of CN109444791A publication Critical patent/CN109444791A/zh
Application granted granted Critical
Publication of CN109444791B publication Critical patent/CN109444791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/02Testing or calibrating of apparatus covered by the other groups of this subclass of auxiliary devices, e.g. of instrument transformers according to prescribed transformation ratio, phase angle, or wattage rating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/005Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
    • G01R35/007Standards or reference devices, e.g. voltage or resistance standards, "golden references"

Abstract

本发明公开了一种电容式电压互感器误差状态在线监测方法及系统,包括:在同步信号单元的同步触发下,采集校准投运的三相电容式电压互感器的测量数据,利用主元分析的方法建立表征电容式电压互感器误差状态的统计特征参数和标准正常运行状态的统计特征量的阈值;采集运行状态下的三相电容式电压互感器的测量数据,根据的统计特征参数,利用可自适应的主元分析的方法表征三相电容式电压互感器误差状态的统计特征量,并与统计阈值比较判断电容式电压互感器的误差状态是否正常;当测量数据的统计特征量超越其统计阈值时利用变量贡献率法对电容式电压互感器的误差异常进行诊断。本发明能够在不停电条件下实现电容式电压互感器误差状态的实时在线评估。

Description

一种电容式电压互感器误差状态评估方法和系统
技术领域
本发明属于电力技术领域,更具体地,涉及一种电容式电压互感器(CapacitorVoltage Transformer,CVT)误差状态评估方法及系统。
背景技术
测量用互感器是变电站内重要的一次主设备,是电能计量装置中的重要组成部分。其中,电容式电压互感器因其良好的绝缘性能和经济性而被广泛用于110kV及以上电压等级的电力系统中。相比于传统的电磁式电压互感器,电容式电压互感器的结构更为复杂,运行过程中的误差稳定性不高,极易出现超差现象,影响到电量交易的公平与准确,严重情况下甚至会威胁电网的安全运行。作为测量设备,计量误差的长期稳定性是衡量电容式电压互感器运行性能最重要的参数之一,需要对运行中电容式电压互感器的计量误差状态进行评估和风险预警。
目前在运电容式电压互感器计量误差状态评估所采用的方式一般是定期离线停电检修的形式,国家标准规定的检修时间为4年。但由于高压输变电线路停电检修的困难,大量电容式电压互感器处于超检定周期运行状态,对电网的安全、稳定、经济运行造成了一定的隐患。传统的定期离线检修的方式已不适合计量误差稳定性较差的电容式电压互感器,需要在不停电的条件下开展电容式电压互感器计量误差的状态评估。
作为测量设备,计量性能是CVT状态评估的最重要内容。状态评估在高压主设备如变压器、断路器上已经有了较为广泛的应用,而在运行互感器上的应用还属于空白,主要原因是:为了评估在运互感器最为关注的误差性能,需要提供高精度的标准器对被评估互感器进行校准,而高精度的标准器由于运行环境的制约,难以投入高压环境中,即使在线投入,也只能短时校准,这也造成了在运互感器的误差状态评估研究基本是空白的现状。
发明内容
针对现有技术的缺陷,本发明的目的在于解决传统的定期离线检修的方式已不适合计量误差稳定性较差的电容式电压互感器,以及评估在运互感器需要采用高精度的标准器,而高精度的标准器由于运行环境的制约,难以投入高压环境中,即使在线投入,也只能短时校准的技术问题。
为实现上述目的,第一方面,本发明提供一种电容式电压互感器计量误差状态评估方法,包括以下步骤:
(1.1)采集校准投运后的三相电容式电压互感器的二次输出信号,构成标准幅值数据矩阵和标准相位数据矩阵:
(1.2)利用主元分析法分解得到步骤(1.1)中的幅值数据矩阵与其主元模型值和其残差模型值之间的关系,以及得到步骤(1.1)中的相位数据矩阵与其主元模型值和其残差模型值之间的关系,利用步骤(1.1)中的幅值数据矩阵的协方差矩阵求得其残差模型值,利用步骤(1.1)中的相位数据矩阵的协方差矩阵求得其残差模型值,根据步骤(1.1)中的幅值数据矩阵和其残差模型值确定计量误差状态评估的第一评估标准量,根据步骤(1.1)中的相位数据矩阵和其残差模型值确定计量误差状态评估的第二评估标准量;
(1.3)采集在线运行的三相电容式电压互感器的实时测量数据,得到对应的在线运行三相电容式电压互感器的幅值数据矩阵和相位数据矩阵,分别根据正常运行状态下建立的幅值数据矩阵和相位数据矩阵的残差模型确定三相电容式电压互感器测量数据的第一统计量和第二统计量;
(1.4)将第一统计量和第二统计量与第一评估标准量和第二评估标准量进行比较,若第一统计量小于或等于第一评估标准量且第二统计量小于或等于第二评估标准量,则此时三相电容式电压互感器的计量误差状态正常,利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量;若第一统计量大于第一评估标准量,表明三相电容式电压互感器的幅值误差状态出现异常,或若第二统计量大于第二评估标准量,表明三相电容式电压互感器的相位误差状态出现异常,则停止更新残差模型、第一评估标准量以及第二评估标准量;
(1.5)当三相电容式电压互感器的幅值计量误差状态存在异常,或三相电容式电压互感器的相位计量误差状态存在异常时,通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器。
在一个可选的实例中,幅值数据矩阵X1和相位数据矩阵X2如下:
其中,uAi、θAi分别指A相电容式电压互感器第i个采集样本的幅值数据和相位数据,uBi、θBi分别指B相电容式互感器第i个采集样本的幅值数据和相位数据,uCi、θCi指C相电容式互感器第i个采集样本的幅值数据和相位数据,1≤i≤n,n为采集的样本数。
在一个可选的实例中,第一统计量、第二统计量、第一评估标准量和第二评估标准量通过以下步骤确定:
(2.1)将数据矩阵X1和X2按照下式进行分解:
式中,为幅值数据矩阵X1的主元模型值,E1为幅值数据矩阵X1的残差模型值;为相位数据矩阵X2的主元模型值,E2为相位数据矩阵X2的残差模型值,T1为幅值数据矩阵X1的主元得分矩阵,P1 T为幅值数据矩阵X1的主元载荷矩阵,T2为相位数据矩阵X2的主元得分矩阵,P2 T为相位数据矩阵X2的主元载荷矩阵,Te1为幅值数据矩阵X1的残差得分矩阵,Pe1 T为幅值数据矩阵X1的残差载荷矩阵,Te2为相位数据矩阵X2的残差得分矩阵,Pe2 T为相位数据矩阵X2的残差载荷矩阵;
通过对数据矩阵的协方差矩阵R进行奇异值分解进行计算数据矩阵的主元模型和残差模型:
式中,Λ1=diag(λ11,λ12,λ13),λ11>λ12>λ13为幅值数据矩阵X1的协方差矩阵R1的特征值,[P1 Pe1]为特征值λ11,λ12,λ13对应的特征向量组成的特征矩阵;Λ2=diag(λ21,λ22,λ23),λ21>λ22>λ23为相位数据矩阵X2的协方差矩阵R2的特征值,[P2 Pe2]为特征值λ21,λ22,λ23对应的特征向量组成的特征矩阵;
(2.2)量化评估标准量:
当所述幅值数据矩阵X1为在线运行三相电容式电压互感器的幅值数据矩阵,当所述相位数据矩阵X2为在线运行三相电容式电压互感器的相位数据矩阵时,获得正常运行状态下测量数据的Q统计量的统计控制限,其中,Q统计量的计算方法为:
其中,Q1为第一统计量,Q2为第二统计量;
当所述幅值数据矩阵X1为标准幅值数据矩阵,当所述相位数据矩阵X2为标准相位数据矩阵时,在显著性水平为α时Q统计量的控制限为:
式中, Cα是正态分布在检测水平为α下的临界值;Qc1表示第一评估标准量,Qc2表示相第二评估标准量。
在一个可选的实例中,步骤(1.4)中利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量,包括以下步骤:
(4.1)矩阵I过渡至矩阵II
从移动窗中抛弃最陈旧的数据矩阵数据矩阵I变成数据矩阵II,此时新的数据矩阵的均值相量和各变量的标准方差分别为:
式中 为数据矩阵的第i列,和bk(i)分别为向量和bk的第i个元素,nk+1为第k+1次数据更新,m为变量个数,表示第i个变量的标准差,L表示数据矩阵的长度;
对数据矩阵II进行标准化处理可得:
式中 为矩阵I标准化处理后的值,为矩阵II标准化处理后的值,数据矩阵I的协方差矩阵Rk与数据矩阵II的协方差矩阵之间的关系为:
式中∑k表示数据矩阵I的标准差矩阵,xk表示第k次数据采样值;
用初始数据矩阵的变量方差来标准化新的采样数据,此时数据矩阵II的协方差矩阵递归更新为:
(4.2)矩阵II过渡至矩阵III
在k+1时刻获得新的数据矩阵此时的数据矩阵II将转变为数据矩阵III:与步骤1类似,数据矩阵III的均值向量bk+1和协方差矩阵Rk+1的递归计算为:
其中将递归更新后的数据矩阵II的协方差矩阵代入数据矩阵III的协方差矩阵Rk+1,得到Rk+1的递归表达式为:
得到新的数据矩阵的协方差矩阵后,根据步骤(2.1)中描述对协方差矩阵进行奇异值分解,求取对应的主元子空间、残差子空间以及统计控制控制限,实现对三相电容式电压互感器计量误差状态评估标准量的递归更新。
在一个可选的实例中,步骤(1.5)中通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器,具体方法为:
其中,X1m、X2m分别为m相电容式电压互感器的幅值数据和相位数据,分别为m相电容式电压互感器的幅值重构数据和相位重构数据,Q1m为m相电容式电压互感器幅值数据对Q统计量的贡献率,Q2m为m相电容式电压互感器相位数据对Q统计量的贡献率,Q1m中的最大值对应相的互感器的幅值数据发生异常,Q2m中的最大值对应相的互感器的相位数据发生异常。
其中,m相可以为A相、B相或C相。
另一方面,本发明提供一种电容式电压互感器计量误差状态评估系统,包括:采集单元、评估标准量确定单元、统计量确定单元、误差状态判断及参数更新单元以及误差定位单元;
采集单元,用于采集校准投运后的三相电容式电压互感器的二次输出信号,构成标准幅值数据矩阵和标准相位数据矩阵:
评估标准量确定单元,用于利用主元分析法分解得到标准幅值数据矩阵与其主元模型值和其残差模型值之间的关系,以及得到标准相位数据矩阵与其主元模型值和其残差模型值之间的关系,利用标准幅值数据矩阵的协方差矩阵求得其残差模型值,利用标准相位数据矩阵的协方差矩阵求得其残差模型值,根据标准幅值数据矩阵和其残差模型值确定计量误差状态评估的第一评估标准量,根据标准相位数据矩阵和其残差模型值确定计量误差状态评估的第二评估标准量;
统计量确定单元,用于采集在线运行的三相电容式电压互感器的实时测量数据,得到对应的在线运行三相电容式电压互感器的幅值数据矩阵和相位数据矩阵,分别根据正常运行状态下建立的幅值数据矩阵和相位数据矩阵的残差模型确定三相电容式电压互感器测量数据的第一统计量和第二统计量;
误差状态判断及参数更新单元,用于将第一统计量和第二统计量与第一评估标准量和第二评估标准量进行比较,若第一统计量小于或等于第一评估标准量且第二统计量小于或等于第二评估标准量,则此时三相电容式电压互感器的计量误差状态正常,利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量;若第一统计量大于第一评估标准量,表明三相电容式电压互感器的幅值误差状态出现异常,或若第二统计量大于第二评估标准量,表明三相电容式电压互感器的相位误差状态出现异常,则停止更新残差模型、第一评估标准量以及第二评估标准量;
误差定位单元,用于当三相电容式电压互感器的幅值计量误差状态存在异常,或三相电容式电压互感器的相位计量误差状态存在异常时,通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器。
在一个可选的实例中,幅值数据矩阵X1和相位数据矩阵X2如下:
其中,uAi、θAi分别指A相电容式电压互感器第i个采集样本的幅值数据和相位数据,uBi、θBi分别指B相电容式互感器第i个采集样本的幅值数据和相位数据,uCi、θCi指C相电容式互感器第i个采集样本的幅值数据和相位数据,1≤i≤n,n为采集的样本数。
在一个可选的实例中,所述第一统计量、第二统计量、第一评估标准量和第二评估标准量通过以下步骤确定:
将数据矩阵X1和X2按照下式进行分解:
式中,为幅值数据矩阵X1的主元模型值,E1为幅值数据矩阵X1的残差模型值;为相位数据矩阵X2的主元模型值,E2为相位数据矩阵X2的残差模型值,T1为幅值数据矩阵X1的主元得分矩阵,P1 T为幅值数据矩阵X1的主元载荷矩阵,T2为相位数据矩阵X2的主元得分矩阵,P2 T为相位数据矩阵X2的主元载荷矩阵,Te1为幅值数据矩阵X1的残差得分矩阵,Pe1 T为幅值数据矩阵X1的残差载荷矩阵,Te2为相位数据矩阵X2的残差得分矩阵,Pe2 T为相位数据矩阵X2的残差载荷矩阵;
通过对数据矩阵的协方差矩阵R进行奇异值分解进行计算数据矩阵的主元模型和残差模型:
式中,Λ1=diag(λ11,λ12,λ13),λ11>λ12>λ13为幅值数据矩阵X1的协方差矩阵R1的特征值,[P1 Pe1]为特征值λ11,λ12,λ13对应的特征向量组成的特征矩阵;Λ2=diag(λ21,λ22,λ23),λ21>λ22>λ23为相位数据矩阵X2的协方差矩阵R2的特征值,[P2 Pe2]为特征值λ21,λ22,λ23对应的特征向量组成的特征矩阵;
量化评估标准量:
当所述幅值数据矩阵X1为在线运行三相电容式电压互感器的幅值数据矩阵,当所述相位数据矩阵X2为在线运行三相电容式电压互感器的相位数据矩阵时,获得正常运行状态下测量数据的Q统计量的统计控制限,其中,Q统计量的计算方法为:
其中,Q1为第一统计量,Q2为第二统计量;
当所述幅值数据矩阵X1为标准幅值数据矩阵,当所述相位数据矩阵X2为标准相位数据矩阵时,在显著性水平为α时Q统计量的控制限为:
式中, Cα是正态分布在检测水平为α下的临界值;Qc1表示第一评估标准量,Qc2表示相第二评估标准量。
在一个可选的实例中,误差状态判断及参数更新单元通过以下步骤更新残差模型、第一评估标准量以及第二评估标准量:
矩阵I过渡至矩阵II
从移动窗中抛弃最陈旧的数据矩阵数据矩阵I变成数据矩阵II,此时新的数据矩阵的均值相量和各变量的标准方差分别为:
式中 为数据矩阵的第i列,和bk(i)分别为向量和bk的第i个元素,nk+1为第k+1次数据更新,m为变量个数,表示第i个变量的标准差,L表示数据矩阵的长度;
对数据矩阵II进行标准化处理可得:
式中 为矩阵I标准化处理后的值,为矩阵II标准化处理后的值,数据矩阵I的协方差矩阵Rk与数据矩阵II的协方差矩阵之间的关系为:
式中∑k表示数据矩阵I的标准差矩阵,xk表示第k次数据采样值;
用初始数据矩阵的变量方差来标准化新的采样数据,此时数据矩阵II的协方差矩阵递归更新为:
矩阵II过渡至矩阵III
在k+1时刻获得新的数据矩阵此时的数据矩阵II将转变为数据矩阵III:与步骤1类似,数据矩阵III的均值向量bk+1和协方差矩阵Rk+1的递归计算为:
其中将递归更新后的数据矩阵II的协方差矩阵代入数据矩阵III的协方差矩阵Rk+1,得到Rk+1的递归表达式为:
得到新的数据矩阵的协方差矩阵后,根据步骤(2.1)中描述对协方差矩阵进行奇异值分解,求取对应的主元子空间、残差子空间以及统计控制控制限,实现对三相电容式电压互感器计量误差状态评估标准量的递归更新。
在一个可选的实例中,误差定位单元通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器,具体包括:
其中,X1m、X2m分别为m相电容式电压互感器的幅值数据和相位数据,分别为m相电容式电压互感器的幅值重构数据和相位重构数据,Q1m为m相电容式电压互感器幅值数据对Q统计量的贡献率,Q2m为m相电容式电压互感器相位数据对Q统计量的贡献率,Q1m中的最大值对应相的互感器的幅值数据发生异常,Q2m中的最大值对应相的互感器的相位数据发生异常。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)本发明提供的电容式电压互感器误差状态在线监测方法及系统,基于设备监测数据对设备状态进行实时连续采集,并进行计量误差状态的分析,同传统误差状态评估方法相比由于评估过程不需要停电,可以有效减少评估成本,具有较高的经济性。
(2)本发明提供的电压互感器误差状态在线监测方法及系统,基于设备监测数据对设备误差状态进行在线评估,同传统的离线评估方法相比可有效减少电网的非故障停电时间,提高电力系统的运行稳定性和可靠性。
附图说明
图1是本发明实施例提供的电容式电压互感器计量误差状态在线评估方法的流程示意图;
图2是本发明实施例提供的电容式电压互感器自适应评估模型的更新策略示意图;
图3是本发明实施例提供的电容式电压互感器计量误差状态在线评估系统示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提出了一种基于自适应主元分析的电容式电压互感器计量误差状态评估方法和系统,其目的在于实现在不停电条件下对在线运行的电容式电压互感器的计量误差状态进行在线评估。CVT的作用是准确测量电网节点的一次电压信号。由于受节点负荷变化和运行环境的影响,电网节点的一次电压信号是一个非平稳随机信号。一般情况下CVT计量误差的变化小于一次物理电压的随机波动,难以将CVT的计量误差变化从电网自身状态随机变化的影响中剥离。
本发明的总体思路是:通过分析在运CVT测量数据信息特征以及电力系统一次电气物理量的波动特点,利用数学分析的方法将电网的一次信息波动与CVT自身运行异常造成的计量偏差相互分离,通过测量数据提取表征CVT运行状态的特征统计量,为测量数据的准确性校验提供标准参量,实现了脱离高级标准电压互感器条件下的CVT误差性能在线准确监测、评估与预警。
本发明提供的基于自适应主元分析的电容式电压互感器计量误差状态评估方法和系统,其评估对象可以为110kV及以上电压等级变电站一个间隔内的三相电容式电压互感器,其思路在于:根据输变电系统三相对称运行的电气物理特性,利用主元分析的方法对三相电容式电压互感器的幅值测量数据和相位测量数据进行相关性分析,建立计量误差状态评估和异常诊断方法,实现不停电条件下对在线运行的电容式电压互感器的计量误差状态进行评估;可有效减少变电站的非故障停电时间;及时判断电容式电压互感器的运行状态,为电容式电压互感器的检修提供数据支撑,保障电网的安全稳定运行。以下结合附图和实施例来具体阐述。
本发明实施例提供的电容式电压互感器计量误差状态评估方法和系统,其流程如图1所示,包括如下步骤:
(1)利用高精度同步信号采集装置,在同步信号单元的触发作用下对三相电容式电压互感器的二次输出信号进行模数转换,计算电容式电压互感器的幅值数据x1=[uAt uBtuCt],相位数据x2=[θAt θBt θCt];uAt、θAt指A相电容式电压互感器的的幅值数据和相位数据,uBt、θBt指B相电容式互感器的的幅值数据和相位数据,uCt、θCt指C相电容式互感器的的幅值数据和相位数据。
(2)幅值数据矩阵X1和相位数据矩阵X2如下所示:
其中,uAn、θAn指A相电容式电压互感器的幅值数据和相位数据,uBn、θBn指B相电容式互感器的幅值数据和相位数据,uCn、θCn指C相电容式互感器的幅值数据和相位数据,n为采集的样本数;
(3)利用主元分析的方法根据上述数据矩阵X1和X2建立三相电容式电压互感器测量数据的残差模型,并计算正常运行状态Q统计量的统计控制限,量化评估标准量,具体包括如下子步骤:
(2.1)将数据矩阵X1和X2按照下式进行分解:
式中,为幅值数据矩阵X1的主元模型值,E1为幅值数据矩阵X1的残差模型值;为幅值数据矩阵X2的主元模型值,E2为幅值数据矩阵X2的残差模型值。
数据矩阵的主元模型和残差模型可以通过对数据矩阵的协方差矩阵R进行奇异值分解进行计算:
式中,Λ1=diag(λ11,λ12,λ13),λ11>λ12>λ13为协方差矩阵R1的特征值,[P1 Pe1]为特征值λ11,λ12,λ13对应的特征向量组成的特征矩阵;Λ2=diag(λ21,λ22,λ23),λ21>λ22>λ23为协方差矩阵R2的特征值,[P2 Pe2]为特征值λ21,λ22,λ23对应的特征向量组成的特征矩阵;
(2.2)量化评估标准量,即获得正常运行状态下测量数据的Q统计量的统计控制限Qα,其中,Q统计量的计算方法为:
在显著性水平为α时Q统计量的控制限为:
式中, Cα是正态分布在检测水平为α下的临界值;Qc1表示第一评估标准量,Qc2表示相第二评估标准量。
(4)基于移动窗的原理更新三相电容式电压互感器的残差模型和统计控制限,使得电容式电压互感器误差状态评估模型对过程的时变特性具备自适应能力,减少误判率。
幅值数据矩阵X1和相位数据矩阵X2具有相同的更新方法,图2是本发明实施例提供的电容式电压互感器自适应评估模型的更新策略示意图,如图2所示,具体步骤为:
步骤1:矩阵I过渡至矩阵II
从移动窗中抛弃最陈旧的数据块数据矩阵I变成数据矩阵II,此时新的数据矩阵的均值相量和各变量的标准方差分别为:
式中 为数据矩阵的第i列,和bk(i)分别为向量和bk的第i个元素。
对数据矩阵II进行标准化处理可得:
式中数据矩阵I的协方差矩阵Rk与数据矩阵II的协方差矩阵之间的关系为:
可以用初始数据矩阵的变量方差来标准化新的采样数据,此时数据矩阵II的协方差矩阵可以递归更新为:
步骤2:矩阵II过渡至矩阵III
在k+1时刻获得新的数据块此时的数据矩阵II将转变为数据矩阵III:与步骤1类似,数据矩阵III的均值向量bk+1和协方差矩阵Rk+1的递归计算为:
其中将式(3-6)代入式(3-8),可以得到Rk+1的递归表达式为:
当得到新的数据矩阵的协方差矩阵后,可根据步骤(3)的描述对协方差矩阵进行奇异值分解,求取对应的主元子空间、残差子空间以及统计控制控制限Qα,实现电容式电压互感器计量误差状态评估标准量的递归更新。
(5)采集三相电容式电压互感器的实时测量数据矩阵,根据步骤(4)计算的残差模型,计算三相电容式电压互感器的实时Q统计量Q1、Q2,并与统计控制Qc1、Qc2进行比较,判定三相电容式电压互感器是否存在互感器的计量误差状态发生异常,若运行正常,则重复步骤(4),更新评估模型。
(6)根据步骤(5),当电容式电压互感器的计量误差状态出现异常时,通过计算统计量的贡献率进行异常互感器的辨识,具体方法为:
其中,X1m、X2m分别为m相电容式电压互感器的幅值数据和相位数据,分别为m电容式电压互感器的幅值重构数据和相位重构数据。其中,m相可以为A相、B相或C相。根据三相电容式电压互感器的测量数据对异常Q统计的贡献率来判断具体哪一个电容式电压互感器的误差状态发生异常,贡献率最大的那一相电容式电压互感器发生异常。可根据识别出的具体某一相电容式电压互感器的误差状态故障来指导相关工作人员进行电容式电压互感器的检修工作。
本发明实施例提供的电容式电压互感器误差状态在线监测系统,其监测对象为智能变电站中一个间隔内的三相电容式式电压互感器,其系统构成如图3所示,包括:采集单元、评估标准量确定单元、统计量确定单元、误差状态判断及参数更新单元以及误差定位单元;
采集单元,用于采集校准投运后的三相电容式电压互感器的二次输出信号,构成标准幅值数据矩阵和标准相位数据矩阵:
评估标准量确定单元,用于利用主元分析法分解得到标准的幅值数据矩阵与其主元模型值和其残差模型值之间的关系,以及得到标准相位数据矩阵与其主元模型值和其残差模型值之间的关系,利用标准幅值数据矩阵的协方差矩阵求得其残差模型值,利用标准相位数据矩阵的协方差矩阵求得其残差模型值,根据标准幅值数据矩阵和其残差模型值确定计量误差状态评估的第一评估标准量,根据标准相位数据矩阵和其残差模型值确定计量误差状态评估的第二评估标准量;
统计量确定单元,用于采集在线运行的三相电容式电压互感器的实时测量数据,得到对应的在线运行三相电容式电压互感器的幅值数据矩阵和相位数据矩阵,分别根据正常运行状态下建立的幅值数据矩阵和相位数据矩阵的残差模型确定三相电容式电压互感器测量数据的第一统计量和第二统计量;
误差状态判断及参数更新单元,用于将第一统计量和第二统计量与第一评估标准量和第二评估标准量进行比较,若第一统计量小于或等于第一评估标准量且第二统计量小于或等于第二评估标准量,则此时三相电容式电压互感器的计量误差状态正常,利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量;若第一统计量大于第一评估标准量,表明三相电容式电压互感器的幅值误差状态出现异常,或若第二统计量大于第二评估标准量,表明三相电容式电压互感器的相位误差状态出现异常,则停止更新残差模型、第一评估标准量以及第二评估标准量;
误差定位单元,用于当三相电容式电压互感器的幅值计量误差状态存在异常,或三相电容式电压互感器的相位计量误差状态存在异常时,通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器。
具体各个单元的功能可参照前述方法实施例中的介绍,在此不做赘述。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种电容式电压互感器计量误差状态评估方法,其特征在于,包括以下步骤:
(1.1)采集校准投运后的三相电容式电压互感器的二次输出信号,构成标准幅值数据矩阵和标准相位数据矩阵:
(1.2)利用主元分析法分解得到步骤(1.1)中的幅值数据矩阵与其主元模型值和其残差模型值之间的关系,以及得到步骤(1.1)中的相位数据矩阵与其主元模型值和其残差模型值之间的关系,利用步骤(1.1)中的幅值数据矩阵的协方差矩阵求得其残差模型值,利用步骤(1.1)中的相位数据矩阵的协方差矩阵求得其残差模型值,根据步骤(1.1)中的幅值数据矩阵和其残差模型值确定计量误差状态评估的第一评估标准量,根据步骤(1.1)中的相位数据矩阵和其残差模型值确定计量误差状态评估的第二评估标准量;
(1.3)采集在线运行的三相电容式电压互感器的实时测量数据,得到对应的在线运行三相电容式电压互感器的幅值数据矩阵和相位数据矩阵,分别根据正常运行状态下建立的幅值数据矩阵和相位数据矩阵的残差模型确定三相电容式电压互感器测量数据的第一统计量和第二统计量;
(1.4)将第一统计量和第二统计量与第一评估标准量和第二评估标准量进行比较,若第一统计量小于或等于第一评估标准量且第二统计量小于或等于第二评估标准量,则此时三相电容式电压互感器的计量误差状态正常,利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量;若第一统计量大于第一评估标准量,表明三相电容式电压互感器的幅值误差状态出现异常,或若第二统计量大于第二评估标准量,表明三相电容式电压互感器的相位误差状态出现异常,则停止更新残差模型、第一评估标准量以及第二评估标准量;
(1.5)当三相电容式电压互感器的幅值计量误差状态存在异常,或三相电容式电压互感器的相位计量误差状态存在异常时,通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器;
幅值数据矩阵X1和相位数据矩阵X2如下:
其中,uAi、θAi分别指A相电容式电压互感器第i个采集样本的幅值数据和相位数据,uBi、θBi分别指B相电容式互感器第i个采集样本的幅值数据和相位数据,uCi、θCi指C相电容式互感器第i个采集样本的幅值数据和相位数据,1≤i≤n,n为采集的样本数;
第一统计量、第二统计量、第一评估标准量和第二评估标准量通过以下步骤确定:
(2.1)将数据矩阵X1和X2按照下式进行分解:
式中,为幅值数据矩阵X1的主元模型值,E1为幅值数据矩阵X1的残差模型值;为相位数据矩阵X2的主元模型值,E2为相位数据矩阵X2的残差模型值,T1为幅值数据矩阵X1的主元得分矩阵,P1 T为幅值数据矩阵X1的主元载荷矩阵,T2为相位数据矩阵X2的主元得分矩阵,P2 T为相位数据矩阵X2的主元载荷矩阵,Te1为幅值数据矩阵X1的残差得分矩阵,Pe1 T为幅值数据矩阵X1的残差载荷矩阵,Te2为相位数据矩阵X2的残差得分矩阵,Pe2 T为相位数据矩阵X2的残差载荷矩阵;
通过对数据矩阵的协方差矩阵R进行奇异值分解进行计算数据矩阵的主元模型和残差模型:
式中,Λ1=diag(λ11,λ12,λ13),λ11>λ12>λ13为幅值数据矩阵X1的协方差矩阵R1的特征值,[P1 Pe1]为特征值λ11,λ12,λ13对应的特征向量组成的特征矩阵;Λ2=diag(λ21,λ22,λ23),λ21>λ22>λ23为相位数据矩阵X2的协方差矩阵R2的特征值,[P2 Pe2]为特征值λ21,λ22,λ23对应的特征向量组成的特征矩阵;
(2.2)量化评估标准量:
当所述幅值数据矩阵X1为在线运行三相电容式电压互感器的幅值数据矩阵,当所述相位数据矩阵X2为在线运行三相电容式电压互感器的相位数据矩阵时,获得正常运行状态下测量数据的Q统计量的统计控制限,其中,Q统计量的计算方法为:
其中,Q1为第一统计量,Q2为第二统计量;
当所述幅值数据矩阵X1为标准幅值数据矩阵,当所述相位数据矩阵X2为标准相位数据矩阵时,在显著性水平为α时Q统计量的控制限为:
式中, Cα是正态分布在检测水平为α下的临界值; Qc1表示第一评估标准量,Qc2表示相第二评估标准量;
步骤(1.4)中利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量,包括以下步骤:
(4.1)矩阵Ⅰ过渡至矩阵Ⅱ
从移动窗中抛弃最陈旧的数据矩阵数据矩阵Ⅰ变成数据矩阵Ⅱ,此时新的数据矩阵的均值相量和各变量的标准方差分别为:
式中 为数据矩阵的第i列,和bk(i)分别为向量和bk的第i个元素,nk+1为第k+1次数据更新,m为变量个数,表示第i个变量的标准差,L表示数据矩阵的长度;
对数据矩阵Ⅱ进行标准化处理可得:
式中 为矩阵Ⅰ标准化处理后的值,为矩阵Ⅱ标准化处理后的值,数据矩阵Ⅰ的协方差矩阵Rk与数据矩阵Ⅱ的协方差矩阵之间的关系为:
式中∑k表示数据矩阵Ⅰ的标准差矩阵,xk表示第k次数据采样值;
用初始数据矩阵的变量方差来标准化新的采样数据,此时数据矩阵Ⅱ的协方差矩阵递归更新为:
(4.2)矩阵Ⅱ过渡至矩阵Ⅲ
在k+1时刻获得新的数据矩阵此时的数据矩阵Ⅱ将转变为数据矩阵Ⅲ:与步骤1类似,数据矩阵Ⅲ的均值向量bk+1和协方差矩阵Rk+1的递归计算为:
其中将递归更新后的数据矩阵Ⅱ的协方差矩阵代入数据矩阵Ⅲ的协方差矩阵Rk+1,得到Rk+1的递归表达式为:
得到新的数据矩阵的协方差矩阵后,根据步骤(2.1)中描述对协方差矩阵进行奇异值分解,求取对应的主元子空间、残差子空间以及统计控制控制限,实现对三相电容式电压互感器计量误差状态评估标准量的递归更新。
2.根据权利要求1所述的电容式电压互感器计量误差状态评估方法,其特征在于,步骤(1.5)中通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器,具体方法为:
其中,X1m、X2m分别为m相电容式电压互感器的幅值数据和相位数据,分别为m相电容式电压互感器的幅值重构数据和相位重构数据,Q1m为m相电容式电压互感器幅值数据对Q统计量的贡献率,Q2m为m相电容式电压互感器相位数据对Q统计量的贡献率,Q1m中的最大值对应相的互感器的幅值数据发生异常,Q2m中的最大值对应相的互感器的相位数据发生异常。
3.一种电容式电压互感器计量误差状态评估系统,其特征在于,包括:采集单元、评估标准量确定单元、统计量确定单元、误差状态判断及参数更新单元以及误差定位单元;
采集单元,用于采集校准投运后的三相电容式电压互感器的二次输出信号,构成标准幅值数据矩阵和标准相位数据矩阵:
评估标准量确定单元,用于利用主元分析法分解得到标准幅值数据矩阵与其主元模型值和其残差模型值之间的关系,以及得到标准相位数据矩阵与其主元模型值和其残差模型值之间的关系,利用标准幅值数据矩阵的协方差矩阵求得其残差模型值,利用标准相位数据矩阵的协方差矩阵求得其残差模型值,根据标准幅值数据矩阵和其残差模型值确定计量误差状态评估的第一评估标准量,根据标准相位数据矩阵和其残差模型值确定计量误差状态评估的第二评估标准量;
统计量确定单元,用于采集在线运行的三相电容式电压互感器的实时测量数据,得到对应的在线运行三相电容式电压互感器的幅值数据矩阵和相位数据矩阵,分别根据正常运行状态下建立的幅值数据矩阵和相位数据矩阵的残差模型确定三相电容式电压互感器测量数据的第一统计量和第二统计量;
误差状态判断及参数更新单元,用于将第一统计量和第二统计量与第一评估标准量和第二评估标准量进行比较,若第一统计量小于或等于第一评估标准量且第二统计量小于或等于第二评估标准量,则此时三相电容式电压互感器的计量误差状态正常,利用当前测量数据更新残差模型、第一评估标准量以及第二评估标准量;若第一统计量大于第一评估标准量,表明三相电容式电压互感器的幅值误差状态出现异常,或若第二统计量大于第二评估标准量,表明三相电容式电压互感器的相位误差状态出现异常,则停止更新残差模型、第一评估标准量以及第二评估标准量;
误差定位单元,用于当三相电容式电压互感器的幅值计量误差状态存在异常,或三相电容式电压互感器的相位计量误差状态存在异常时,通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器;
幅值数据矩阵X1和相位数据矩阵X2如下:
其中,uAi、θAi分别指A相电容式电压互感器第i个采集样本的幅值数据和相位数据,uBi、θBi分别指B相电容式互感器第i个采集样本的幅值数据和相位数据,uCi、θCi指C相电容式互感器第i个采集样本的幅值数据和相位数据,1≤i≤n,n为采集的样本数;
所述第一统计量、第二统计量、第一评估标准量和第二评估标准量通过以下步骤确定:
(2.1)将数据矩阵X1和X2按照下式进行分解:
式中,为幅值数据矩阵X1的主元模型值,E1为幅值数据矩阵X1的残差模型值;为相位数据矩阵X2的主元模型值,E2为相位数据矩阵X2的残差模型值,T1为幅值数据矩阵X1的主元得分矩阵,P1 T为幅值数据矩阵X1的主元载荷矩阵,T2为相位数据矩阵X2的主元得分矩阵,P2 T为相位数据矩阵X2的主元载荷矩阵,Te1为幅值数据矩阵X1的残差得分矩阵,Pe1 T为幅值数据矩阵X1的残差载荷矩阵,Te2为相位数据矩阵X2的残差得分矩阵,Pe2 T为相位数据矩阵X2的残差载荷矩阵;
通过对数据矩阵的协方差矩阵R进行奇异值分解进行计算数据矩阵的主元模型和残差模型:
式中,Λ1=diag(λ11,λ12,λ13),λ11>λ12>λ13为幅值数据矩阵X1的协方差矩阵R1的特征值,[P1 Pe1]为特征值λ11,λ12,λ13对应的特征向量组成的特征矩阵;Λ2=diag(λ21,λ22,λ23),λ21>λ22>λ23为相位数据矩阵X2的协方差矩阵R2的特征值,[P2 Pe2]为特征值λ21,λ22,λ23对应的特征向量组成的特征矩阵;
(2.2)量化评估标准量:
当所述幅值数据矩阵X1为在线运行三相电容式电压互感器的幅值数据矩阵,当所述相位数据矩阵X2为在线运行三相电容式电压互感器的相位数据矩阵时,获得正常运行状态下测量数据的Q统计量的统计控制限,其中,Q统计量的计算方法为:
其中,Q1为第一统计量,Q2为第二统计量;
当所述幅值数据矩阵X1为标准幅值数据矩阵,当所述相位数据矩阵X2为标准相位数据矩阵时,在显著性水平为α时Q统计量的控制限为:
式中, Cα是正态分布在检测水平为α下的临界值;Qc1表示第一评估标准量,Qc2表示相第二评估标准量;
误差状态判断及参数更新单元通过以下步骤更新残差模型、第一评估标准量以及第二评估标准量:
(4.1)矩阵Ⅰ过渡至矩阵Ⅱ
从移动窗中抛弃最陈旧的数据矩阵数据矩阵Ⅰ变成数据矩阵Ⅱ,此时新的数据矩阵的均值相量和各变量的标准方差分别为:
式中 为数据矩阵的第i列,和bk(i)分别为向量和bk的第i个元素,nk+1为第k+1次数据更新,m为变量个数,表示第i个变量的标准差,L表示数据矩阵的长度;
对数据矩阵Ⅱ进行标准化处理可得:
式中 为矩阵Ⅰ标准化处理后的值,为矩阵Ⅱ标准化处理后的值,数据矩阵Ⅰ的协方差矩阵Rk与数据矩阵Ⅱ的协方差矩阵之间的关系为:
式中∑k表示数据矩阵Ⅰ的标准差矩阵,xk表示第k次数据采样值;
用初始数据矩阵的变量方差来标准化新的采样数据,此时数据矩阵Ⅱ的协方差矩阵递归更新为:
(4.2)矩阵Ⅱ过渡至矩阵Ⅲ
在k+1时刻获得新的数据矩阵此时的数据矩阵Ⅱ将转变为数据矩阵Ⅲ:与步骤1类似,数据矩阵Ⅲ的均值向量bk+1和协方差矩阵Rk+1的递归计算为:
其中将递归更新后的数据矩阵Ⅱ的协方差矩阵代入数据矩阵Ⅲ的协方差矩阵Rk+1,得到Rk+1的递归表达式为:
得到新的数据矩阵的协方差矩阵后,根据(2.1)中描述对协方差矩阵进行奇异值分解,求取对应的主元子空间、残差子空间以及统计控制控制限,实现对三相电容式电压互感器计量误差状态评估标准量的递归更新。
4.根据权利要求3所述的电容式电压互感器计量误差状态评估系统,其特征在于,误差定位单元通过三相电容式电压互感器中各相互感器的误差贡献率确定异常互感器,具体包括:
其中,X1m、X2m分别为m相电容式电压互感器的幅值数据和相位数据,分别为m相电容式电压互感器的幅值重构数据和相位重构数据,Q1m为m相电容式电压互感器幅值数据对Q统计量的贡献率,Q2m为m相电容式电压互感器相位数据对Q统计量的贡献率,Q1m中的最大值对应相的互感器的幅值数据发生异常,Q2m中的最大值对应相的互感器的相位数据发生异常。
CN201811332075.4A 2018-11-09 2018-11-09 一种电容式电压互感器误差状态评估方法和系统 Active CN109444791B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811332075.4A CN109444791B (zh) 2018-11-09 2018-11-09 一种电容式电压互感器误差状态评估方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811332075.4A CN109444791B (zh) 2018-11-09 2018-11-09 一种电容式电压互感器误差状态评估方法和系统

Publications (2)

Publication Number Publication Date
CN109444791A CN109444791A (zh) 2019-03-08
CN109444791B true CN109444791B (zh) 2020-09-18

Family

ID=65551333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811332075.4A Active CN109444791B (zh) 2018-11-09 2018-11-09 一种电容式电压互感器误差状态评估方法和系统

Country Status (1)

Country Link
CN (1) CN109444791B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109902735A (zh) * 2019-02-22 2019-06-18 武汉格蓝若智能技术有限公司 一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法
CN110082698A (zh) * 2019-04-16 2019-08-02 贵州电网有限责任公司 一种电容式电压互感器综合运行状态评估仿真系统
CN110082699B (zh) * 2019-05-10 2021-03-19 国网天津市电力公司电力科学研究院 一种低压台区智能电能表运行误差计算方法及其系统
CN111507579B (zh) * 2020-03-26 2022-04-15 国网四川省电力公司电力科学研究院 一种基于聚类法的全类别电压互感器运行性能评估方法
CN112285634B (zh) * 2020-12-18 2021-03-23 华中科技大学 一种双母线结构中高压电压互感器测量误差异常识别方法
CN113050017A (zh) * 2021-03-02 2021-06-29 合肥工业大学 一种电子式互感器误差状态智能监测及故障诊断系统
CN113239132B (zh) * 2021-07-12 2021-09-24 武汉格蓝若智能技术有限公司 一种电压互感器的超差在线辨识方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104598552A (zh) * 2014-12-31 2015-05-06 大连钜正科技有限公司 一种支持增量式更新的大数据特征学习的方法
CN106771645B (zh) * 2016-12-12 2019-06-25 国网四川省电力公司电力科学研究院 电容式电压互感器介损及电容量在线监测方法及监测系统
CN107255792A (zh) * 2017-06-20 2017-10-17 华中科技大学 一种电子式电压互感器误差在线监测方法及系统
CN107944747A (zh) * 2017-12-11 2018-04-20 国网江苏省电力公司宿迁供电公司 一种基于改进主成分分析的低压台区评价方法

Also Published As

Publication number Publication date
CN109444791A (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
CN109444791B (zh) 一种电容式电压互感器误差状态评估方法和系统
CN104483575B (zh) 用于非侵入式电力监测的自适应负荷事件检测方法
CN107016236B (zh) 基于非线性量测方程的电网假数据注入攻击检测方法
CN111474510B (zh) 一种非平稳输出的电压互感器的误差评估方法及系统
CN204064509U (zh) 一种高压开关柜温度在线监测系统
CN112285634B (zh) 一种双母线结构中高压电压互感器测量误差异常识别方法
CN112098915B (zh) 双母线分段接线下多台电压互感器继发性误差的评估方法
CN109239546A (zh) 一种变压器绝缘寿命预测与可靠性评估方法
CN108549047A (zh) 一种基于加速退化试验的电能表寿命评估方法
CN111507579B (zh) 一种基于聚类法的全类别电压互感器运行性能评估方法
CN111025222A (zh) 基于递推主元分析的电子式电压互感器状态监测方法
CN113050017A (zh) 一种电子式互感器误差状态智能监测及故障诊断系统
CN112485747A (zh) 一种电容式电压互感器误差状态评估方法
CN109406923B (zh) 一种列车供电系统接地漏电故障预测方法及装置
Li et al. Real time evaluation algorithm for measurement performance of substation voltage transformer based on artificial neural network
CN109901022A (zh) 基于同步量测数据的配电网区域定位方法
CN106524896A (zh) 基于电路阻抗法的变压器绕组变形在线监测方法
CN103575982B (zh) 一种实时频率稳定度分析方法
Feng et al. State Assessment of Measuring Error and Analysis of Risk Early Warning Technology for Operating CVT
Xu et al. Real-time evaluation of error for electronic voltage transformers based on autocorrelation and large-dimensional random matrices
He et al. A Method to Predict the Trends of Operating Performance of Capacitive Voltage Transformers
Sujie et al. Research on the Model of CVT Operation Performance Based on the PHM System
CN110579323B (zh) 一种用于高压断路器绝缘气体的机器人测值识别算法
Ma et al. Battery Resistance Compensation in Access System
Bai et al. Application of Transformer Online Monitoring Data in State Evaluation

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211122

Address after: Room 1803-1805, building 2-07, guanggu.core center, 303 Guanggu Avenue, Donghu New Technology Development Zone, Wuhan City, Hubei Province, 430000

Patentee after: WUHAN GLORY ROAD INTELLIGENT TECHNOLOGY Co.,Ltd.

Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037

Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY