CN108896863B - 一种频率响应绕组变形分析用的线性相关系数计算方法 - Google Patents

一种频率响应绕组变形分析用的线性相关系数计算方法 Download PDF

Info

Publication number
CN108896863B
CN108896863B CN201810503405.5A CN201810503405A CN108896863B CN 108896863 B CN108896863 B CN 108896863B CN 201810503405 A CN201810503405 A CN 201810503405A CN 108896863 B CN108896863 B CN 108896863B
Authority
CN
China
Prior art keywords
phase
frequency response
winding
samples
data
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
CN201810503405.5A
Other languages
English (en)
Other versions
CN108896863A (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
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Liaoning 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, Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201810503405.5A priority Critical patent/CN108896863B/zh
Publication of CN108896863A publication Critical patent/CN108896863A/zh
Application granted granted Critical
Publication of CN108896863B publication Critical patent/CN108896863B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/72—Testing of electric windings

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)

Abstract

本发明涉及电力变压器绕组变形检测技术领域，尤其涉及一种频率响应绕组变形分析用的线性相关系数计算方法。本发明是根据被诊断变压器两次不同时间测量的变压器频率响应数据，通过计算逐相纵向的线性相关系数、计算横向或相间线性相关系数的方法判断绕组变形情况；或者根据被诊断变压器本次测量的频率响应数据计算横向或相间线性相关系数；换算变形判断需要的相关系数，依据设定的判断标准评估绕组变形程度。本发明根据频率响应法获得的曲线数据，可以快速、准确的计算判断变形程度的相关系数，结合曲线中峰值变化能够判断出绕组的故障类型及程度，为评估变压器的健康状况提供依据。

Description

一种频率响应绕组变形分析用的线性相关系数计算方法
技术领域
本发明涉及电力变压器绕组变形检测技术领域，尤其涉及一种频率响应绕组变形分析用的线性相关系数计算方法。
背景技术
电力变压器是电力***中非常重要和昂贵的设备，若变压器发生故障将影响整个电力***运行的稳定性和可靠性。据统计，大约25％以上的变压器故障由绕组变形故障引起，即绕组发生变形后，会给大型电力变压带来严重的故障隐患，若变压器再次遭受短路电流冲击，往往会引起严重的停电事故，造成重大的社会经济损失。因此变压器绕组变形的诊断对保证变压器的安全稳定运行及电网的供电可靠性有重大意义。
目前，国内外先后提出的用于检测变压器绕组变形的方法主要有低电压阻抗法、电容法、频率响应法及其结合多种检测手段的综合法等。上述的方法在实际应用中各有其优缺点，但在时间检测中有一定成功的经验。
频率响应分析法最早由加拿大的Dick提出，20世纪90年代初我国开始讲频率响应法检测绕组变形情况，标准DL/T 911 2004针对频率响应法检测绕组变形做出了详细的测试要求，但是给出的计算相关系数计算分析方法有一定缺陷。
发明内容
针对上述相关系数计算分析方法中存在的问题，本发明提出了一种频率响应绕组变形分析用的线性相关系数计算方法，其目的是为了提供一种计算快速、准确，适应性强，具有较高实际应用价值的频率响应分析用的线性相关系数计算方法。
为了实现上述发明目，本发明是通过以下技术方案来实现的：
一种频率响应绕组变形分析用的线性相关系数计算方法，是根据被诊断变压器两次不同时间测量的变压器频率响应数据，通过计算逐相纵向的线性相关系数、计算横向或相间线性相关系数的方法判断绕组变形情况；或者根据被诊断变压器本次测量的频率响应数据计算横向或相间线性相关系数；换算变形判断需要的相关系数，依据设定的判断标准评估绕组变形程度。
所述两次不同时间测量的变压器频率响应数据，包括：
设被诊断变压器的相数为m；
两次不同时间测量的变压器频率响应数据分为上次测量的频率响应数据及本次测量数据；
上次测量的频率响应数据：
如果m＝3，上次第x个绕组A相的频率响应数据为X1xA,第x个绕组B相的频率响应数据为X1xB,第x个绕组C相的频率响应数据为X1xC；
如果m＝1，上次第x个绕组的频率响应数据为X1x；
本次次测量的频率响应数据：
如果m＝3，本次第x个绕组A相的频率响应数据为X2xA,第x个绕组B相的频率响应数据为X2xB,第x个绕组C相的频率响应数据为X2xC；
如果m＝1，本次第x个绕组的频率响应数据为X2x。
所述线性相关系数的计算方法包括：
第1：纵向比较的线性相关系数计算；
以三相变压器高压绕组A相表示：
计算两个样本的均值E(X1xA)及E(X2xA)
Figure BDA0001670592570000021
式(1)中，k为样本的个数，X1xA(i)为上次测量数据中第i个频率响应值，i为样本的序号，E(X1xA)为样本X1xA的均值，X1xA为高压绕组A相上次测量数据，X2xA为高压绕组A相本次次测量数据；
Figure BDA0001670592570000022
式(2)中，k为样本的个数，X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为样本X2xA的均值，X2xA为高压绕组A相本次次测量数据，i为样本的序号；
计算两个样本的均方差D(X1xA)及D(X2xA):
Figure BDA0001670592570000023
式(3)中，k为样本的个数，i为样本的序号,X1xA(i)为上次测量数据中第i个频率响应值，E(X1xA)为上次测量A相数据的均值；
Figure BDA0001670592570000024
式(4)中，k为样本的个数，i为样本的序号,X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量B相数据的均值；
计算两个样本的协方差
Figure BDA0001670592570000025
Figure BDA0001670592570000031
式(5)中，k为样本的个数，i为样本的序号,X1xA(i)为上次测量数据中第i个频率响应值，E(X1xA)为上次测量A相数据的均值，X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量B相数据的均值；
高压侧A相频率响应的线性相关系数
Figure BDA0001670592570000038
Figure BDA0001670592570000032
式(6)中，两个样本的均方差分别为D(X1xA)及D(X2xA)，计算两个样本的协方差为
Figure BDA0001670592570000039
同理可以求得高压侧B、C相纵向频率响应的线性相关系数ρ纵高B、ρ纵高C，可以求得其他绕组纵向频率响应的线性相关系数；
第2：横向比较的线性相关系数计算；
以三相变压器高压绕组A、B相横向比较为例表示：
计算两个样本的均值E(X2xA)及E(X2xB)
Figure BDA0001670592570000033
式(1)中，k为样本的个数，i为样本的序号,X2xA为高压绕组A相本次次测量数据，X2xA(i)为本次测量A相数据中第i个频率响应值；
Figure BDA0001670592570000034
式(2)中，k为样本的个数，i为样本的序号，X2xB为高压绕组B相本次次测量数据,X2xB(i)为本次测量B相数据中第i个频率响应值；
计算两个样本的均方差D(X2xA)及D(X2xB)；
Figure BDA0001670592570000035
式(3)中，k为样本的个数，i为样本的序号，X2xA为高压绕组A相本次次测量数据，X2xA(i)为本次测量A相数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值；
Figure BDA0001670592570000036
式(4)中，k为样本的个数，i为样本的序号，X2xB为高压绕组B相本次次测量数据,X2xB(i)为本次测量B相数据中第i个频率响应值，E(X2xB)为本次测量B相数据的均值；
计算两个样本的协方差Cov横高AB
Figure BDA0001670592570000037
式(5)中，k为样本的个数，i为样本的序号，X2xA(i)为本次测量A相数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值，X2xB(i)为本次测量B相数据中第i个频率响应值，E(X2xB)为本次测量B相数据的均值；
高压侧AB相频率响应的线性相关系数ρ横高AB
Figure BDA0001670592570000041
式(6)中，D(X2xA)及D(X2xB)为两个样本的均方差，
Figure BDA0001670592570000043
为两个样本的协方差；
同理可以求得高压侧BC、CA横向频率响应的线性相关系数ρ横高BC、ρ横高CA，可以求得其他绕组纵向频率响应的线性相关系数。
所述变形判断需要的相关系数的计算包括：
不论是三相变压器还是单相变压器，无论是纵向比较还是横向比较，最终变形相关系数的计算方法相同；
变形相关系数为V：
Figure BDA0001670592570000042
通过本方法计算变压器各侧不同绕组高、中、低频段的纵向变形相关系数；计算三相绕组两相间高、中、低频段的变形相关系数；低频段为1kHz～100kHz，中频段为100kHz～600kHz，低频段为600kHz～1000kHz。
所述变形程度的判断标准包括：
对于纵向比较绕组变形程度及横向绕组变形程度判断采用相同的标准，
(1)绕组变形程度：严重变形；
变形相关系数(V)：V低频<0.6；
(2)绕组变形程度：明显变形；
变形相关系数(V)：0.6≤V低频<1.0或V中频<0.6；
(3)绕组变形程度：轻微变形；
变形相关系数(V)：1.0≤V低频<2.0或0.6≤V中频<1.0；
(4)绕组变形程度：无变形；
变形相关系数(V)：V低频≥2.0和V中频≥1.0和V高频≥2.0；
其中：
V低频为曲线在低频段(1kHz～100kHz)内的变形相关系数；
V中频为曲线在中频段(100kHz～600kHz)内的变形相关系数；
V高频为曲线在高频段(600kHz～1000kHz)内的变形相关系数。
所述变形判断需要的相关系数的计算方法还包括：
V1＝a-n|ρ|-n+…+a-2|ρ|-2+a-1|ρ|-1+a0+a1|ρ|+a2|ρ|2+…+an|ρ|n
上式中：变形相关系数为V1，a-n，a-(n-1)，…，a-1，a1，…，a(n-1)，an不全为0的实数，ρ为两条曲线的线性相关系数，n为趋近无穷大的正整数；
V2＝logb(a-n|ρ|-n+…+a-2|ρ|-2+a-1|ρ|-1+a0+a1|ρ|+a2|ρ|2+…+an|ρ|n)
上式中：：变形相关系数为V2，a-n，a-(n-1)，…，a-1，a1，…，a(n-1)，an不全为0的实数，b为正实数，ρ为两条曲线的线性相关系数，n为趋近无穷大的正整数；
低频段为V低频，满足ω1≤V低频≤ω2；
中频段为V中频，满足ω2<V中频≤ω3；
高频段为V高频，满足ω3<V高频≤ω4；
上式中0<ω1<ω2<ω3<ω4≤2000kHz；
上式中ω表示响应绕组变形试验外施正弦激励源的角频率。
本发明具有如下优点及有益效果：
根据频率响应法获得的曲线数据，可以快速、准确的计算判断变形程度的相关系数，结合曲线中峰值变化能够判断出绕组的故障类型及程度，为评估变压器的健康状况提供依据。
下面结合附图和具体实施例，对本发明作进一步详细的说明，但不受本实施例所限。
附图说明
图1是本发明变压器绕组变形相关系数计算及判断流程图。
具体实施方式
实施例1：
本发明是一种频率响应绕组变形分析用的线性相关系数计算方法，根据被诊断变压器两次不同时间测量的变压器频率响应数据，通过计算逐相纵向的线性相关系数、计算横向或相间线性相关系数的方法判断绕组变形情况；或者根据被诊断变压器本次测量的频率响应数据计算横向或相间线性相关系数；换算变形判断需要的相关系数，依据设定的判断标准评估绕组变形程度。具体计算及判断方法如图1所示，下面详细介绍。
被诊断变压器的相数为m。
一、两次不同时间测量的变压器频率响应数据。
两次不同时间测量的变压器频率响应数据分为上次测量的频率响应数据及本次测量数据。
上次测量的频率响应数据：
如果m＝3，上次第x个绕组A相的频率响应数据为X1xA,第x个绕组B相的频率响应数据为X1xB,第x个绕组C相的频率响应数据为X1xC。
如果m＝1，上次第x个绕组的频率响应数据为X1x。
本次次测量的频率响应数据：
如果m＝3，本次第x个绕组A相的频率响应数据为X2xA,第x个绕组B相的频率响应数据为X2xB,第x个绕组C相的频率响应数据为X2xC。
如果m＝1，本次第x个绕组的频率响应数据为X2x。
二、线性相关系数的计算方法。
1、纵向比较的线性相关系数计算。
以三相变压器高压绕组A相为例表示：
计算两个样本的均值E(X1xA)及E(X2xA)
Figure BDA0001670592570000061
式(1)中，k为样本的个数，X1xA(i)为上次测量数据中第i个频率响应值，i为样本的序号，E(X1xA)为样本X1xA的均值，X1xA为高压绕组A相上次测量数据，X2xA为高压绕组A相本次次测量数据。
Figure BDA0001670592570000062
式(2)中，k为样本的个数，X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为样本X2xA的均值，X2xA为高压绕组A相本次次测量数据，i为样本的序号。
计算两个样本的均方差D(X1xA)及D(X2xA):
Figure BDA0001670592570000063
式(3)中，k为样本的个数，i为样本的序号,X1xA(i)为上次测量数据中第i个频率响应值，E(X1xA)为上次测量A相数据的均值。
Figure BDA0001670592570000064
式(4)中，k为样本的个数，i为样本的序号,X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量B相数据的均值。
计算两个样本的协方差
Figure BDA0001670592570000065
Figure BDA0001670592570000071
式(5)中，k为样本的个数，i为样本的序号,X1xA(i)为上次测量数据中第i个频率响应值，E(X1xA)为上次测量A相数据的均值，X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量B相数据的均值。
高压侧A相频率响应的线性相关系数
Figure BDA0001670592570000078
Figure BDA0001670592570000072
式(6)中，两个样本的均方差分别为D(X1xA)及D(X2xA)，计算两个样本的协方差为
Figure BDA0001670592570000079
同理可以求得高压侧B、C相纵向频率响应的线性相关系数ρ纵高B、ρ纵高C，可以求得其他绕组纵向频率响应的线性相关系数。
2、横向比较的线性相关系数计算。
以三相变压器高压绕组A、B相横向比较为例表示：
计算两个样本的均值E(X2xA)及E(X2xB)
Figure BDA0001670592570000073
式(1)中，k为样本的个数，i为样本的序号,X2xA为高压绕组A相本次次测量数据，X2xA(i)为本次测量A相数据中第i个频率响应值。
Figure BDA0001670592570000074
式(2)中，k为样本的个数，i为样本的序号，X2xB为高压绕组B相本次次测量数据,X2xB(i)为本次测量B相数据中第i个频率响应值。
计算两个样本的均方差D(X2xA)及D(X2xB)
Figure BDA0001670592570000075
式(3)中，k为样本的个数，i为样本的序号，X2xA为高压绕组A相本次次测量数据，X2xA(i)为本次测量A相数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值。
Figure BDA0001670592570000076
式(4)中，k为样本的个数，i为样本的序号，X2xB为高压绕组B相本次次测量数据,X2xB(i)为本次测量B相数据中第i个频率响应值，E(X2xB)为本次测量B相数据的均值。
计算两个样本的协方差Cov横高AB
Figure BDA0001670592570000077
式(5)中，k为样本的个数，i为样本的序号，X2xA(i)为本次测量A相数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值，X2xB(i)为本次测量B相数据中第i个频率响应值，E(X2xB)为本次测量B相数据的均值。
高压侧AB相频率响应的线性相关系数ρ横高AB
Figure BDA0001670592570000081
式(6)中，D(X2xA)及D(X2xB)为两个样本的均方差，
Figure BDA0001670592570000083
为两个样本的协方差。
同理可以求得高压侧BC、CA横向频率响应的线性相关系数ρ横高BC、ρ横高CA，可以求得其他绕组纵向频率响应的线性相关系数。
三、变形相关系数计算。
不论是三相变压器还是单相变压器，无论是纵向比较还是横向比较，最终变形相关系数的计算方法相同。
变形相关系数为V：
Figure BDA0001670592570000082
通过本方法可以计算变压器各侧不同绕组高、中、低频段的纵向变形相关系数；可以计算三相绕组两相间高、中、低频段的变形相关系数。低频段为1kHz～100kHz，中频段为100kHz～600kHz，低频段为600kHz～1000kHz。
四、变形程度判断标准。
对于纵向比较绕组变形程度及横向绕组变形程度判断采用相同的标准，具体标准如表1所示。
本发明的特定实施例已对本发明的内容做出了详尽的说明，但不局限本实施例，本领域技术人员根据本发明的启示所做的任何显而易见的改动，都属于本发明权利保护的范围。
实施例2：
在实施1的前两个步骤的基础上本发明中所述另两种计算变形相关系数计算方法，是：
V1＝a-n|ρ|-n+…+a-2|ρ|-2+a-1|ρ|-1+a0+a1|ρ|+a2|ρ|2+…+an|ρ|n
上式中：变形相关系数为V1，a-n，a-(n-1)，…，a-1，a1，…，a(n-1)，an不全为0的实数，ρ为两条曲线的线性相关系数，n为趋近无穷大的正整数。
V2＝logb(a-n|ρ|-n+…+a-2|ρ|-2+a-1|ρ|-1+a0+a1|ρ|+a2|ρ|2+…+an|ρ|n)
上式中：：变形相关系数为V2，a-n，a-(n-1)，…，a-1，a1，…，a(n-1)，an不全为0的实数，b为正实数，ρ为两条曲线的线性相关系数，n为趋近无穷大的正整数。
低频段为V低频，满足ω1≤V低频≤ω2；
中频段为V中频，满足ω2<V中频≤ω3；
高频段为V高频，满足ω3<V高频≤ω4。
上式中0<ω1<ω2<ω3<ω4≤2000kHz。
上式中ω表示响应绕组变形试验外施正弦激励源的角频率。
表1：纵向比较绕组变形程度及横向绕组变形程度判断标准表。
Figure BDA0001670592570000091

Claims (1)

1.一种频率响应绕组变形分析用的线性相关系数计算方法，其特征是：根据被诊断变压器两次不同时间测量的变压器频率响应数据，通过计算逐相纵向的线性相关系数的方法判断绕组变形情况；或者根据被诊断变压器本次测量的频率响应数据计算横向或相间线性相关系数；换算变形判断需要的相关系数，依据设定的判断标准评估绕组变形程度；
所述两次不同时间测量的变压器频率响应数据，包括：
设被诊断变压器的相数为m；
两次不同时间测量的变压器频率响应数据分为上次测量的频率响应数据及本次测量数据；
上次测量的频率响应数据：
如果m＝3，上次第x个绕组A相的频率响应数据为X1xA,第x个绕组B相的频率响应数据为X1xB,第x个绕组C相的频率响应数据为X1xC；
如果m＝1，上次第x个绕组的频率响应数据为X1x；
本次测量的频率响应数据：
如果m＝3，本次第x个绕组A相的频率响应数据为X2xA,第x个绕组B相的频率响应数据为X2xB,第x个绕组C相的频率响应数据为X2xC；
如果m＝1，本次第x个绕组的频率响应数据为X2x：
所述线性相关系数的计算方法包括：
第1：纵向比较的线性相关系数计算；
以三相变压器高压绕组A相表示：
计算两个样本的均值E(X1xA)及E(X2xA)
Figure FDA0002809517450000011
式(1)中，k为样本的个数，X1xA(i)为上次测量数据中第i个频率响应值，i为样本的序号，E(X1xA)为上次测量A相数据的均值，X1xA为上次第x个绕组A相的频率响应数据，X2xA为本次第x个绕组A相的频率响应数据；
Figure FDA0002809517450000012
式(2)中，k为样本的个数，X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值，X2xA为本次第x个绕组A相的频率响应数据，i为样本的序号；
计算两个样本的均方差D(X1xA)及D(X2xA):
Figure FDA0002809517450000021
式(3)中，k为样本的个数，i为样本的序号,X1xA(i)为上次测量数据中第i个频率响应值，E(X1xA)为上次测量A相数据的均值；
Figure FDA0002809517450000022
式(4)中，k为样本的个数，i为样本的序号,X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值；
两个样本的协方差Cov纵高A
Figure FDA0002809517450000023
式(5)中，k为样本的个数，i为样本的序号,X1xA(i)为上次测量数据中第i个频率响应值，E(X1xA)为上次测量A相数据的均值，X2xA(i)为本次测量数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值；
高压侧A相频率响应的线性相关系数ρ纵高A
Figure FDA0002809517450000024
式(6)中，两个样本的均方差分别为D(X1xA)及D(X2xA)，两个样本的协方差为Cov纵高A；
同理求得高压侧B、C相纵向频率响应的线性相关系数ρ纵高B、ρ纵高C；
第2：横向比较的线性相关系数计算；
以三相变压器高压绕组A、B相横向比较为例表示：
计算两个样本的均值E(X2xA)及E(X2xB)
Figure FDA0002809517450000025
式(1)中，k为样本的个数，i为样本的序号,X2xA为本次第x个绕组A相的频率响应数据，X2xA(i)为本次测量A相数据中第i个频率响应值；
Figure FDA0002809517450000026
式(2)中，k为样本的个数，i为样本的序号，X2xB为本次第x个绕组B相的频率响应数据,X2xB(i)为本次测量B相数据中第i个频率响应值；
计算两个样本的均方差D(X2xA)及D(X2xB)；
Figure FDA0002809517450000027
式(3)中，k为样本的个数，i为样本的序号，X2xA为本次第x个绕组A相的频率响应数据，X2xA(i)为本次测量A相数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值；
Figure FDA0002809517450000031
式(4)中，k为样本的个数，i为样本的序号，X2xB为本次第x个绕组B相的频率响应数据,X2xB(i)为本次测量B相数据中第i个频率响应值，E(X2xB)为本次测量B相数据的均值；
计算两个样本的协方差Cov横高AB
Figure FDA0002809517450000032
式(5)中，k为样本的个数，i为样本的序号，X2xA(i)为本次测量A相数据中第i个频率响应值，E(X2xA)为本次测量A相数据的均值，X2xB(i)为本次测量B相数据中第i个频率响应值，E(X2xB)为本次测量B相数据的均值；
高压侧AB相频率响应的线性相关系数ρ横高AB
Figure FDA0002809517450000033
式(6)中，D(X2xA)及D(X2xB)为两个样本的均方差，Cov纵高AB为两个样本的协方差；
同理求得高压侧BC、CA横向频率响应的线性相关系数ρ横高BC、ρ横高CA；
所述变形判断需要的相关系数的计算包括：
不论是三相变压器还是单相变压器，无论是纵向比较还是横向比较，最终变形相关系数的计算方法相同；
变形相关系数为V：
Figure FDA0002809517450000034
通过本方法计算变压器各侧不同绕组高、中、低频段的纵向变形相关系数；计算三相绕组两相间高、中、低频段的变形相关系数；低频段为1kHz～100kHz，中频段为100kHz～600kHz，高频段为600kHz～1000kHz；
所述变形程度的判断标准包括：
对于纵向比较绕组变形程度及横向绕组变形程度判断采用相同的标准，
(1)绕组变形程度：严重变形；
变形相关系数(V)：V低频＜0.6；
(2)绕组变形程度：明显变形；
变形相关系数(V)：0.6≤V低频＜1.0或V中频＜0.6；
(3)绕组变形程度：轻微变形；
变形相关系数(V)：1.0≤V低频＜2.0或0.6≤V中频＜1.0；
(4)绕组变形程度：无变形；
变形相关系数(V)：V低频≥2.0和V中频≥1.0和V高频≥0.6；
其中：
V低频为曲线在低频段(1kHz～100kHz)内的变形相关系数；
V中频为曲线在中频段(100kHz～600kHz)内的变形相关系数；
V高频为曲线在高频段(600kHz～1000kHz)内的变形相关系数；
所述变形判断需要的相关系数的计算方法还包括：
V1＝a-n|ρ|-n+…+a-2|ρ|-2+a-1|ρ|-1+a0+a1|ρ|+a2|ρ|2+…+an|ρ|n
上式中：变形相关系数为V1，a-n，a-(n-1)，…，a-1，a1，…，a(n-1)，an不全为0的实数，ρ为两条曲线的线性相关系数，n为趋近无穷大的正整数；
V2＝logb(a-n|ρ|-n+…+a-2|ρ|-2+a-1|ρ|-1+a0+a1|ρ|+a2|ρ|2+…+an|ρ|n)
上式中：变形相关系数为V2，a-n，a-(n-1)，…，a-1，a1，…，a(n-1)，an不全为0的实数，b为正实数，ρ为两条曲线的线性相关系数，n为趋近无穷大的正整数；
低频段为V低频，满足ω1≤V低频≤ω2；
中频段为V中频，满足ω2＜V中频≤ω3；
高频段为V高频，满足ω3＜V高频≤ω4；
上式中0＜ω1＜ω2＜ω3＜ω4≤2000kHz；
上式中ω表示响应绕组变形试验外施正弦激励源的角频率。
CN201810503405.5A 2018-05-23 2018-05-23 一种频率响应绕组变形分析用的线性相关系数计算方法 Active CN108896863B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810503405.5A CN108896863B (zh) 2018-05-23 2018-05-23 一种频率响应绕组变形分析用的线性相关系数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810503405.5A CN108896863B (zh) 2018-05-23 2018-05-23 一种频率响应绕组变形分析用的线性相关系数计算方法

Publications (2)

Publication Number Publication Date
CN108896863A CN108896863A (zh) 2018-11-27
CN108896863B true CN108896863B (zh) 2021-02-23

Family

ID=64343196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810503405.5A Active CN108896863B (zh) 2018-05-23 2018-05-23 一种频率响应绕组变形分析用的线性相关系数计算方法

Country Status (1)

Country Link
CN (1) CN108896863B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110161381B (zh) * 2019-04-29 2021-04-13 云南电网有限责任公司电力科学研究院 一种基于振荡波的变压器套管绝缘受潮状态评估方法
CN110361611B (zh) * 2019-06-23 2020-09-01 西南交通大学 变压器绕组径向变形下频率响应测试平台及其测评方法
CN110609198A (zh) * 2019-09-11 2019-12-24 广东电网有限责任公司 一种电缆线路的运维监测方法及相关设备
CN111856343B (zh) * 2019-12-27 2021-06-29 广东电网有限责任公司电力科学研究院 一种变压器绕组变形分析方法、装置、设备和存储介质
CN111983524B (zh) * 2020-08-26 2021-06-08 西南交通大学 一种基于振荡波时频变换的变压器绕组故障评估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106338237A (zh) * 2015-07-10 2017-01-18 华北电力大学（保定） 基于频响阻抗法的变压器绕组变形检测方法
CN106679558A (zh) * 2016-12-19 2017-05-17 国网内蒙古东部电力有限公司检修分公司 一种变压器绕组变形的评价方法
CN107478949A (zh) * 2017-08-07 2017-12-15 中国电力科学研究院 一种变压器绕组变形在线诊断方法及***
JP2018042358A (ja) * 2016-09-06 2018-03-15 株式会社三社電機製作所 試験用電源装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6565858B2 (ja) * 2016-10-11 2019-08-28 Tdk株式会社 ワイヤレス電力伝送装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106338237A (zh) * 2015-07-10 2017-01-18 华北电力大学（保定） 基于频响阻抗法的变压器绕组变形检测方法
JP2018042358A (ja) * 2016-09-06 2018-03-15 株式会社三社電機製作所 試験用電源装置
CN106679558A (zh) * 2016-12-19 2017-05-17 国网内蒙古东部电力有限公司检修分公司 一种变压器绕组变形的评价方法
CN107478949A (zh) * 2017-08-07 2017-12-15 中国电力科学研究院 一种变压器绕组变形在线诊断方法及***

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
电力变压器绕组变形的频率响应分析法;中华人民共和国国家发展和改革委员会;《中华人民共和国电力行业标准》;20041214;3-6 *

Also Published As

Publication number Publication date
CN108896863A (zh) 2018-11-27

Similar Documents

Publication Publication Date Title
CN108896863B (zh) 一种频率响应绕组变形分析用的线性相关系数计算方法
WO2020258835A1 (zh) 超导电缆pplp绝缘老化程度的检测方法及***
Liu et al. A study of the sweep frequency impedance method and its application in the detection of internal winding short circuit faults in power transformers
CN110361686B (zh) 基于多参数的电容式电压互感器故障检测方法
CN105093132B (zh) 一种大功率整流器开路故障快速在线诊断方法
Jacob et al. On-line monitoring of partial discharges in a HVDC station environment
CN102890226B (zh) 电力***xlpe电缆水树老化状态测试***
CN103576061A (zh) 一种用于变压器放电性故障诊断的方法
CN105093054B (zh) 一种大功率整流器开关管直通故障快速在线诊断方法
CN102735959A (zh) 超高压线路电力变压器不拆线试验方法
CN110988525B (zh) 一种短时间内重复短路对电力变压器冲击的评估方法
CN102087316B (zh) 一种v/x接线牵引变压器短路电抗在线监测的方法
CN112666231A (zh) 一种换流变压器固体绝缘含水量的测试方法
CN115267640A (zh) 一种电流互感器异常校验方法、***、设备和介质
CN116224210A (zh) 一种交直流共存环境下电流互感器校验***及方法
CN105429187B (zh) 一种多直流馈入评估方法
CN114509649A (zh) 一种线圈类设备匝间绝缘缺陷诊断方法及***
CN112595941B (zh) 一种换流变压器绕组层间绝缘监测方法及***
Jacob et al. Online partial discharge measurement of a high-voltage direct current converter wall-bushing
Zhang et al. Leakage inductance variation based monitoring of transformer winding deformation
CN104977488A (zh) 一种基于差动电流梯度角近似熵的变压器励磁涌流识别方法
CN108732450A (zh) 大型变压器卷铁心片间短路的检测***及方法
CN111413640B (zh) 高阻抗变压器的差动保护接线检测方法及装置
Soltani et al. Compensation of the effects of electrical sensors in measuring partial discharge signals
CN113156276A (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