扣除背景影响量的运行三相电压互感器误差计算方法.pdf

《扣除背景影响量的运行三相电压互感器误差计算方法.pdf》由会员分享，可在线阅读，更多相关《扣除背景影响量的运行三相电压互感器误差计算方法.pdf（12页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010988824.X (22)申请日 2020.09.18 (71)申请人 国网重庆市电力公司营销服务中心 地址 401123 重庆市渝北区青枫北路20号 凤凰D座 申请人 国网重庆市电力公司电力科学研究 院 国家电网有限公司 (72)发明人 周峰冯凌程瑛颖田芝华 肖冀张家铭江金洋杜杰 周华勇谭时顺常仕亮 (74)专利代理机构 北京众合诚成知识产权代理 有限公司 11246 代理人 贺春林 (51)Int.Cl. G01R 35/02(2006.01) (54)发明名称 。

2、一种扣除背景影响量的运行三相电压互感 器误差计算方法 (57)摘要 一种扣除背景影响量的运行三相电压互感 器误差计算方法， 包括以下步骤： 采集两组运行 三相电压互感器的计量绕组三相信号， 获得两组 三相电压互感器采集的数据矩阵X1、 X2， 运行的 两组三相电压互感器为同一电压等级， 且运行的 两组三相电压互感器的三相信号为同相序三相 电压等比例变化信号； 对采集到的两组三相信号 数据矩阵X1、 X2进行数据处理， 将处理后的两组 三相电压相减得到差值数据X； 利用差值数据X开 展优化主元分析计算， 得到效应量的加权离均差 平方和Q的统计量； 利用Q统计量判定运行三相电 压互感器的异常； 本。

3、申请将原有对运行中电压互 感器误差计算差异从千分位提高到万分位， 准确 度提升一个数量级。 权利要求书2页 说明书5页 附图4页 CN 112068064 A 2020.12.11 CN 112068064 A 1.一种扣除背景影响量的运行三相电压互感器误差计算方法， 其特征在于， 具体步骤 如下： 1)采集两组运行三相电压互感器的计量绕组三相信号， 获得两组三相电压互感器采集 的数据矩阵X1、 X2， 运行的两组三相电压互感器为同一电压等级， 且运行的两组三相电压互 感器的三相信号为同相序三相电压等比例变化信号； 2)对采集到的两组三相信号数据矩阵X1、 X2进行数据处理， 将处理后的两组三。

4、相电压相 减得到差值数据X； 3)利用差值数据X开展优化主元分析计算， 得到效应量的加权离均差平方和Q的统计 量； 4)利用Q统计量判定三相电压互感器的运行状态。 2.如权利要求1所述的扣除背景影响量的运行三相电压互感器误差计算方法， 其特征 在于， 步骤1)中所述运行三相电压互感器为同一电压等级的具体条件是： 三相电压互感器一 次电压等级范围为351000kV， 准确度等级为0.2级或0.5级； 步骤1)中所述同相序三相电压等比例变化的具体条件是： 多组三相电压互感器的A相 电压、 B相电压和C相电压， 各自等比例变化； 步骤1)中所述采集两组运行三相电压互感器的计量绕组三相信号为模拟信号采。

5、样， 具 体方法为： 采用模拟采样电路板实现采样， 高速采集频率大于等于12.8kHz， 两组三相电压 互感器的数据矩阵X1、 X2； 其中： 其中：为第一组中A相电压互感器；为第一组中B相电压互感器；为第一组中C 相电压互感器；为第二组中A相电压互感器；为第二组中B相电压互感器；为第二 组中C相电压互感器； n为采集的样本数。 3.如权利要求1所述的扣除背景影响量的运行三相电压互感器误差计算方法， 其特征 在于， 对采集到的两组三相信号数据矩阵X1、 X2进行数据处理， 将处理后的两组三相电压相 减得到差值数据X的具体方法为： 差值数据X计算为： XX1-X2 (1)。 4.如权利要求1所述。

6、的扣除背景影响量的运行三相电压互感器误差计算方法， 其特征 在于， 步骤3)中利用差值数据X开展优化主元分析计算的具体步骤如下： 3-1)将差值数据矩阵X进行分解： 式中，为差值数据矩阵X的主元子空间模型， ETePeT为数据矩阵X的残差子空 间模型； T为主元得分矩阵， P为主元载荷矩阵， Te为残差得分矩阵， Pe为残差载荷矩阵； 权利要求书 1/2 页 2 CN 112068064 A 2 3-2)对差值数据矩阵X的协方差矩阵进行奇异值分解， 得到载荷矩阵P和Pe： RXTX/(N-1)PPePPeT (3) 式中， diag( 1， 2， . m)， 1 2. m为协方差矩阵R的特征值。

7、， PPe为对应 的特征向量组成的负荷向量； 3-3)计算统计量Q统计量: Q统计量为效应量的加权离均差平方和， 即加权平方和； Q统计量反映总离散度； 3-4)计算显著性水平 的Q统计量控制阀值QC。 5.如权利要求1所述的扣除背景影响量的运行三相电压互感器误差计算方法， 其特征 在于， 步骤4)中利用Q统计量判定三相电压互感器运行状态的具体方法如下： 若QQC， 则判定三相电压互感测量误差正常； 若QQC， 则判定三相电压互感中有互感器的测量误差异常。 权利要求书 2/2 页 3 CN 112068064 A 3 一种扣除背景影响量的运行三相电压互感器误差计算方法 技术领域 0001 本发。

8、明涉及电压互感器误差计算领域， 特别是一种扣除背景影响量的运行三相电 压互感器误差计算方法。 背景技术 0002 电力系统的稳定运行与安全可靠发展是重点民生问题。 统计数据显示， 2012年全 国电力消费总量达49726.64亿千瓦小时， 2015年全国电力消费总量达58019.91亿千瓦小 时， 涨幅达16.67， 为社会的全面发展做出了举足轻重的作用。 在电力系统迅速发展的同 时， 用户对电网的运行可靠性提出了越来越高的要求； 作为电力系统中关键的信号提供设 备， 互感器的作用是实现一次侧高压部分与二次设备可靠的电气隔离， 并在保证二次设备 和用电安全的情况下实现一次电压、 电流的准确测量。

9、， 以便提供电能计量、 状态监控和继电 保护的可靠依据， 保证电力系统的安全、 稳定和经济运行； 然而， 长期运行发现互感器在长 期运行中受到外界和内在双重不稳定因素影响， 在实际运行过程中会出现超差现象， 特别 是电容式电压互感器(CVT)。 根据多年的现场运行经验， 110kV及以上电压等级的互感器中， CVT的故障率约为电磁式电压互感器的五倍， 是电磁式电流互感器的10倍； 因此， 需要对运 行中电压互感器的准确性开展评估分析。 发明内容 0003 本发明的目的就是提供一种扣除背景影响量的运行三相电压互感器误差计算方 法， 他可以用于电压互感器误差的计算。 0004 本发明的目的是通过这。

10、样的技术方案实现的， 它包括以下步骤： 0005 1)采集两组运行三相电压互感器的计量绕组三相信号， 获得两组三相电压互感器 采集的数据矩阵X1、 X2， 运行的两组三相电压互感器为同一电压等级， 且运行的两组三相电 压互感器的三相信号为同相序三相电压等比例变化信号； 0006 2)对采集到的两组三相信号数据矩阵X1、 X2进行数据处理， 将处理后的两组三相电 压相减得到差值数据X； 0007 3)利用差值数据X开展优化主元分析计算， 得到效应量的加权离均差平方和Q的统 计量； 0008 4)利用Q统计量判定三相电压互感器的运行状态。 0009 进一步， 步骤1)中所述运行三相电压互感器为同一。

11、电压等级的具体条件是： 三相 电压互感器一次电压等级范围为351000kV， 准确度等级为0.2级或0.5级； 0010 步骤1)中所述同相序三相电压等比例变化的具体条件是： 多组三相电压互感器的 A相电压、 B相电压和C相电压， 各自等比例变化； 0011 步骤1)中所述采集两组运行三相电压互感器的计量绕组三相信号为模拟信号采 样， 具体方法为： 采用模拟采样电路板实现采样， 高速采集频率大于等于12.8kHz， 两组三相 电压互感器的数据矩阵X1、 X2； 其中： 说明书 1/5 页 4 CN 112068064 A 4 0012 0013其中：为第一组中A相电压互感器；为第一组中B相电压。

12、互感器；为第一 组中C相电压互感器；为第二组中A相电压互感器；为第二组中B相电压互感器；为 第二组中C相电压互感器； n为采集的样本数。 0014 进一步， 对采集到的两组三相信号数据矩阵X1、 X2进行数据处理， 将处理后的两组 三相电压相减得到差值数据X的具体方法为： 0015 差值数据X计算为： 0016 XX1-X2 (1)。 0017 进一步， 步骤3)中利用差值数据X开展优化主元分析计算的具体步骤如下： 0018 3-1)将差值数据矩阵X进行分解： 0019 0020式中，为差值数据矩阵X的主元子空间模型， ETePeT为数据矩阵X的残差 子空间模型。 T为主元得分矩阵， P为主元。

13、载荷矩阵， Te为残差得分矩阵， Pe为残差载荷矩阵； 0021 3-2)对差值数据矩阵X的协方差矩阵进行奇异值分解， 得到载荷矩阵P和Pe： 0022 RXTX/(N-1)PPePPeT (3) 0023 式中， diag( 1， 2， . m)， 1 2. m为协方差矩阵R的特征值， PPe为 对应的特征向量组成的负荷向量； 0024 3-3)计算统计量Q统计量: 0025 Q(XPePeT)(XPePeT)TXPePeTXT (4) 0026 Q统计量为效应量的加权离均差平方和， 即加权平方和； Q统计量反映总离散度； 0027 3-4)计算显著性水平 的Q统计量控制阀值QC。 0028。

14、 进一步， 步骤4)中利用Q统计量判定三相电压互感器运行状态的具体方法如下： 0029 若QQC， 则判定三相电压互感无测量误差异常； 0030 若QQC， 则判定三相电压互感中有互感器的测量误差异常。 0031 由于采用了上述技术方案， 本发明具有如下的优点： 0032 1、 本专利利用电源等比例变化时电压互感器的背景影响量相同原理， 设计了一种 扣除背景影响量的运行电压互感器误差计算方法， 将原有对运行中电压互感器误差计算差 异从千分位提高到万分位， 准确度提升一个数量级； 0033 2、 本专利将同电压等级同相序的两组或多组电压幅值相减， 再进行优化的主元分 析， 这种利用同电压同相序电。

15、压之间的差值开展优化主元分析计算， 不仅对电压互感器有 效， 对计算运行中电流互感器同样有效； 0034 3、 本专利对硬件要求优于0.05级， 本专利设计的软件算法优于0.01级， 被监测的 电压互感器计量绕组为0.2级， 硬件搭载软件后高于被监测电压互感器2个准确度等级， 分 析数据的不确定度不影响计算出的被监测电压互感器误差； 说明书 2/5 页 5 CN 112068064 A 5 0035 4、 本专利安装在在变电站控制室内电压互感器误差特性在线监测系统内， 实现对 运行中电压互感器误差的实时分析计算， 本算法对智能处理器要求较低， 普通的FPGA即可 搭载本算法， 可实现对电压互感。

16、器运行误差的边缘计算。 0036 5、 本专利采用数据并行分析原理， 将传统的线性主元分析计算需要采集一个月的 运行电压互感器数据量作为分析数据库支撑量缩短到一周数据可支撑开展误差分析， 大幅 提高计算效率。 0037 本发明的其他优点、 目标和特征在某种程度上将在随后的说明书中进行阐述， 并 且在某种程度上， 基于对下文的考察研究对本领域技术人员而言将是显而易见的， 或者可 以从本发明的实践中得到教导。 本发明的目标和其他优点可以通过下面的说明书和权利要 求书来实现和获得。 附图说明 0038 本发明的附图说明如下： 0039 图1为本发明的流程图。 0040 图2为本发明实验一中三相电压采。

17、集数据的统计图。 0041 图3为本发明实验一中三相电压Q统计值的统计图。 0042 图4为本发明实验二中A相电压采集数据部分数据减少了0.2的统计图。 0043 图5为本发明实验二中三相电压Q统计值的统计图。 0044 图6为本发明实验三中 母和母的A相电压波形幅值采集数据的统计图。 0045 图7为本发明实验三中A相电压Q统计值的统计图。 具体实施方式 0046 下面结合附图和实施例对本发明作进一步说明。 0047 一种扣除背景影响量的运行三相电压互感器误差计算方法， 如图1所示， 具体步骤 如下： 0048 1)采集两组运行三相电压互感器的计量绕组三相信号， 获得两组三相电压互感器 采集。

18、的数据矩阵X1、 X2， 运行的两组三相电压互感器为同一电压等级， 且运行的两组三相电 压互感器的三相信号为同相序三相电压等比例变化信号； 0049 2)对采集到的两组三相信号数据矩阵X1、 X2进行数据处理， 将处理后的两组三相电 压相减得到差值数据X； 0050 3)利用差值数据X开展优化主元分析计算， 得到效应量的加权离均差平方和Q的统 计量； 0051 4)利用Q统计量判定三相电压互感器的运行状态。 0052 步骤1)中所述运行三相电压互感器为同一电压等级的具体条件是： 三相电压互感 器一次电压等级范围为351000kV， 准确度等级为0.2级或0.5级； 0053 步骤1)中所述同相。

19、序三相电压等比例变化的具体条件是： 多组三相电压互感器的 A相电压、 B相电压和C相电压， 各自等比例变化； 0054 步骤1)中所述采集两组运行三相电压互感器的计量绕组三相信号为模拟信号采 样， 具体方法为： 采用模拟采样电路板实现采样， 高速采集频率大于等于12.8kHz， 两组三相 说明书 3/5 页 6 CN 112068064 A 6 电压互感器的数据矩阵X1、 X2； 其中： 0055 0056其中：为第一组中A相电压互感器；为第一组中B相电压互感器；为第一 组中C相电压互感器；为第二组中A相电压互感器；为第二组中B相电压互感器；为 第二组中C相电压互感器； n为采集的样本数。 0。

20、057 步骤2)中对采集到的两组三相信号数据矩阵X1、 X2进行数据处理， 将处理后的两组 三相电压相减得到差值数据X的具体方法为： 0058 差值数据X计算为： 0059 XX1-X2 (1)。 0060 步骤3)中利用差值数据X开展优化主元分析计算的具体步骤如下： 0061 3-1)将差值数据矩阵X进行分解： 0062 0063式中，为差值数据矩阵X的主元子空间模型， ETePeT为数据矩阵X的残差 子空间模型。 T为主元得分矩阵， P为主元载荷矩阵， Te为残差得分矩阵， Pe为残差载荷矩阵； 0064 3-2)对差值数据矩阵X的协方差矩阵进行奇异值分解， 得到载荷矩阵P和Pe： 006。

21、5 RXTX/(N-1)PPePPeT (3) 0066 式中， diag( 1， 2， . m)， 1 2. m为协方差矩阵R的特征值， PPe为 对应的特征向量组成的负荷向量； 0067 3-3)计算统计量Q统计量: 0068 Q(XPePeT)(XPePeT)TXPePeTXT (4) 0069 Q统计量为效应量的加权离均差平方和， 即加权平方和； Q统计量反映总离散度； 0070 3-4)计算显著性水平 的Q统计量控制阀值QC。 0071 步骤4)中利用Q统计量判定三相电压互感器运行状态的具体方法如下： 0072 若QQC， 则判定三相电压互感无测量误差异常； 0073 若QQC， 则。

22、判定三相电压互感中有互感器的测量误差异常。 0074 实验一： 0075 采集某220kV变电站线路 母三相运行中三相电压互感器采集的电压数据片段， 如 图2所示， 通过主元分析法求得三相电压的主元分析Q值， 如图3所示； 可见， 传统的主元分析 表示的是三相输出电压的综合变化情况。 0076 实验二： 0077 采集某220kV变电站线路 母三相运行中三相电压互感器采集的电压数据片段， 将 A相电压采集数据的部分数据减少0.2， 如图4所示， 通过主元分析得到的Q统计量， 如图5 所示； 由图5可知， 部分Q值下降， 可知若某相电压互感器的电源出现异动后会影响Q值计算 的准确度， 这就是背景。

23、误差产生的计算偏差， 同时证明传统的主元分析的准确性依赖于电 说明书 4/5 页 7 CN 112068064 A 7 源的稳定。 0078 实验三： 0079 某220kV变电站的 母和母并联运行， 及存在两组电压互感器来自于同一个电压 等级， 获取 母A相和母A相的电压互感器计量绕组输出电压信号， 如图6所示， 将两组电压 值标记为M1和M2， 采用扣除背景影响量的优化主元分析得到M1和M2的Q统计量， 如图7所示， 比较实验一中的图3和实验三中的图7可知， 传统主元分析的Q值变化范围在-0.015+ 0.035， 变化在千分位， 优化后的主元分析Q值变化范围在-0.0023+0.007，。

24、 变化在万分位。 之后再对Q值按原来的主元分析方式求解电压互感器的运行误差值， 证明了本专利对于运 行误差计算准确度上有较大提升。 0080 本领域内的技术人员应明白， 本申请的实施例可提供为方法、 系统、 或计算机程序 产品。 因此， 本申请可采用完全硬件实施例、 完全软件实施例、 或结合软件和硬件方面的实 施例的形式。 而且， 本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机 可用存储介质(包括但不限于磁盘存储器、 CD-ROM、 光学存储器等)上实施的计算机程序产 品的形式。 0081 本申请是参照根据本申请实施例的方法、 设备(系统)、 和计算机程序产品的流程 图和/或方框。

25、图来描述的。 应理解可由计算机程序指令实现流程图和/或方框图中的每一流 程和/或方框、 以及流程图和/或方框图中的流程和/或方框的结合。 可提供这些计算机程序 指令到通用计算机、 专用计算机、 嵌入式处理机或其他可编程数据处理设备的处理器以产 生一个机器， 使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实 现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。 0082 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特 定方式工作的计算机可读存储器中， 使得存储在该计算机可读存储器中的指令产生包括指 令装置的制造品， 该指令装置实现。

26、在流程图一个流程或多个流程和/或方框图一个方框或 多个方框中指定的功能。 0083 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上， 使得在计 算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理， 从而在计算机或 其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一 个方框或多个方框中指定的功能的步骤。 0084 最后应当说明的是： 以上实施例仅用以说明本发明的技术方案而非对其限制， 尽 管参照上述实施例对本发明进行了详细的说明， 所属领域的普通技术人员应当理解： 依然 可以对本发明的具体实施方式进行修改或者等同替换， 而未脱离本发明精神和范围的任何 修改或者等同替换， 其均应涵盖在本发明的权利要求保护范围之内。 说明书 5/5 页 8 CN 112068064 A 8 图1 说明书附图 1/4 页 9 CN 112068064 A 9 图2 图3 说明书附图 2/4 页 10 CN 112068064 A 10 图4 图5 说明书附图 3/4 页 11 CN 112068064 A 11 图6 图7 说明书附图 4/4 页 12 CN 112068064 A 12 。