CN106950470B - 一种基于大数据的变压器雷电冲击的故障诊断方法 - Google Patents
一种基于大数据的变压器雷电冲击的故障诊断方法 Download PDFInfo
- Publication number
- CN106950470B CN106950470B CN201710142739.XA CN201710142739A CN106950470B CN 106950470 B CN106950470 B CN 106950470B CN 201710142739 A CN201710142739 A CN 201710142739A CN 106950470 B CN106950470 B CN 106950470B
- Authority
- CN
- China
- Prior art keywords
- data
- transformer
- transfer function
- lightning impulse
- test
- 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
Links
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/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Relating To Insulation (AREA)
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
一种基于大数据的变压器雷电冲击的故障诊断方法，包括收集变压器雷电冲击耐压试验的历史数据，对试验测量的首端电压数据和中性点电流数据分别进行快速傅里叶变换。建立故障标准数据库，进行标准化处理得到传递函数数据。将新的雷电冲击耐压试验数据进行标准化处理传递函数数据，通过50％雷电冲击电压与全压下的传递函数数据对比，判断是否发生故障。选取故障数据库中相同型号变压器的相关数据，将新试验的传递函数与故障数据库中的传递函数数据求取相关系数ρ。按照排序的故障列表选取最接近的故障信息，设置偏差值η％，以故障对应位置为中心选取总匝数η％的绕组段为测试范围进行故障定位诊断。本发明对大量故障波形数据进行特征归类，可以确定故障大致范围，可大大缩减故障定位时间，使故障诊断更加快速有效。
Description
技术领域
本发明一种基于大数据的变压器雷电冲击的故障诊断方法，涉及变压器故障检测领域。
背景技术
变压器出厂前都需要进行雷电冲击耐压试验的考核，雷电冲击试验的故障诊断运用最多的是中性点电流法和传递函数法。主要通过50％电压与全压下的试验数据进行对比。对于较小的故障，通过首端电压和中性点电流法也难以判断，此时需要对电压和电流时域波形进行快速傅里叶变换得到各自的频谱，然后用电流信号的频谱除以电压信号的频谱得到传递函数值，通过传递函数波形可以判断较小的故障，但传递函数只与频率相关，不能反映故障时间，也无法确定故障在变压器中的位置。确定变压器存在故障后，可采用超高频局放信号或采用超声局放等故障诊断方法进行故障定位，但对于大型变压器，由于超高频或超声局放等方法每次针对性检测的范围较小，若不能确定故障的大致范围，需要经过多次排查才能进行故障定位，将耗费大量时间和人力物力。
发明内容
为解决上述技术问题，本发明提供一种基于大数据的变压器雷电冲击的故障诊断方法，对大量故障波形数据进行特征归类，虽然雷电冲击由于放电回路间隙击穿的随机性，每个波形都会有差异，但可以确定故障大致范围，可大大缩减故障定位时间，使故障诊断更加快速有效。
本发明采取的技术方案为：
一种基于大数据的变压器雷电冲击的故障诊断方法，包括以下步骤：
步骤1、收集变压器雷电冲击耐压试验的历史数据，按公式1对试验测量的首端电压数据和中性点电流数据分别进行快速傅里叶变换，式中x(n)为离散信号数据，如电压或是电流采样数据，X(k)为信号的快速傅里叶变换数据，为傅里叶变换的核函数，N为数据采样点数，e为自然对数的底数，k为各个采样数据的序号。然后将变换得到的电流信号频谱I(f)除以电压信号频谱U(f)，得到变压器的传递函数H(f)，传递函数也是一组个数为N的数据，每个数据包括频率和对应的幅值，选取对应频率段为0～1MHz的传递函数数据。
步骤2、建立故障标准数据库，包括：试验变压器的型号参数，包括变压器额定电压，额定容量，变压器绕组结构和匝数；雷电冲击耐压试验的记录，包括采样频率、采样点数、变压器绕组首端电压和中性点电流数据；按照步骤1计算得到的0～1MHz频段内传递函数数据；故障信息，即最后通过检测确定的故障类型与故障位置。
步骤3、将待诊断的雷电冲击耐压试验的首端电压和中性点电流数据分别进行快速傅里叶变换，按步骤1求取0～1MHz频率段内的传递函数数据，通过50％雷电冲击电压与全压下的传递函数数据对比，判断是否发生故障。
步骤4、待诊断的雷电冲击耐压试验数据通过传递函数对比判断存在故障，则选取故障数据库中相同型号变压器的相关数据，按公式2将新试验的传递函数与故障数据库中的传递函数数据求取相关系数ρ，式中，x(n)为新的故障传递函数数据，y(n)为故障数据库中的传递函数数据，xi和yi分别为x(n)和y(n)中的第i项数据，和分别为x(n)和y(n)的平均值。设置阈值T，当相关系数ρ≥T时，认为两者类似，将ρ≥T的故障类型按相关系数从大到小的顺序列表。
步骤5、按照排序的故障列表选取最接近的故障信息，设置偏差值η％，以故障对应位置为中心选取总匝数η％的绕组段为测试范围，利用超高频局放测试或超声局放测试进行故障定位诊断；若未检测处，则按列表依次选择下一个接近的故障信息，按上述方法选择测试范围进行测试，直到确定故障位置，完成诊断。
本发明一种基于大数据的变压器雷电冲击的故障诊断方法，利用变压器雷电冲击试验大量的历史数据进行标准化处理，求取雷电冲击下的变压器传递函数建立故障特征数据库，对于待诊断的雷电冲击试验数据，判断发生故障后，利用待诊断故障与数据库中故障数据的传递函数相关系数进行对比，推导新故障的大致类型与范围，然后利用超高频局放测试或超声局放测试进行更准确的故障定位，可减少故障诊断工作量，缩短故障定位时间。
附图说明
图1为本发明的流程图。
图2为雷电冲击下变压器绕组首端电压波形图。
图3为雷电冲击下变压器中性点电流波形图。
图4为50％雷电冲击电压与全压下的传递函数数据对比图。
具体实施方式
如图1所示，一种基于大数据的变压器雷电冲击的故障诊断方法，包括以下步骤：
步骤1)：收集变压器雷电冲击耐压试验的历史数据，例如试验中首端电压和中性点电流数据采样频率为100MHz，记录200μs时长的数据采样点数为20000个点，根据采样点数据绘制雷电冲击下变压器绕组首端电压波形如图1，绘制中性点电流波形如图2。按公式1对试验测量的首端电压数据和中性点电流数据分别进行快速傅里叶变换，式中x(n)为离散信号数据，如电压或是电流采样数据，X(k)为信号的快速傅里叶变换数据，为傅里叶变换的核函数，N为数据采样点数，e为自然对数的底数,k为各个采样数据的序号。然后将变换得到的电流信号频谱I(f)除以电压信号频谱U(f)，得到变压器的传递函数H(f)，传递函数也是一组个数为N的数据，每个数据包括频率和对应的幅值，选取对应频率段为0～1MHz的传递函数数据，如图4所示。
傅里叶变换过后得到的数据就是横坐标为频率，纵坐标为幅值或相角的数据，这里将电压、电流两个傅里叶变换后的数据数值相除，得到的仍是横坐标为频率，纵坐标为数值的图，因此选取对应频率的数值为传递函数数据。
由于雷电冲击波形波头很陡，采样频率较高，如100MHz，记录100μs就需要10000个点，一般都要记录更长，例如采样20000个点。
步骤2)：建立故障标准数据库，如表1所示，主要包括：试验变压器的型号参数，包括变压器额定电压，额定容量，变压器绕组结构和匝数；雷电冲击耐压试验的记录，包括采样频率、采样点数、变压器绕组首端电压和中性点电流数据；按照步骤1计算得到的0～1MHz频段内传递函数数据；故障信息，即最后通过检测确定的故障类型与故障位置。
步骤3)：将待诊断的雷电冲击耐压试验的首端电压和中性点电流数据分别进行快速傅里叶变换，按步骤1求取0～1MHz频率段内的传递函数数据，通过50％雷电冲击电压与全压下试验的传递函数数据对比，判断是否发生故障。
此处新的雷电冲击耐压试验数据指的是刚进行试验获得的待判断是否发生故障的试验数据。标准化处理就是将数据快速傅里叶变换。
步骤4)：待诊断的雷电冲击耐压试验数据通过传递函数对比判断存在故障，则选取故障数据库中相同型号变压器的相关数据，按公式2将新试验的传递函数与故障数据库中的传递函数数据求取相关系数ρ，式中，x(n)为新的故障传递函数数据，y(n)为故障数据库中的传递函数数据，xi和yi分别为x(n)和y(n)中的第i项数据，和分别为x(n)和y(n)的平均值。设置阈值T，当相关系数ρ≥T时，认为两者类似，将ρ≥T的故障类型按相关系数从大到小的顺序列表。
步骤5)：按照排序的故障列表选取最接近的故障信息，设置偏差值η％，以故障对应位置为中心选取总匝数η％的绕组段为测试范围，利用超高频局放测试或超声局放测试进行故障定位诊断。若未检测出，则按列表依次选择下一个接近的故障信息，按上述方法选择测试范围进行测试，直到确定故障位置，完成诊断。
具体算例：
本实例以一台10kV单相双绕组的试验变压器为例进行分析，此台变压器上设置了故障模拟抽头，可以模拟匝间短路，主要步骤如下：
1)、收集变压器雷电冲击耐压试验的历史数据，按公式1对试验测量的首端电压数据和中性点电流数据分别进行快速傅里叶变换，式中x(n)为离散信号数据，如电压或是电流采样数据，X(k)为信号的快速傅里叶变换数据，为傅里叶变换的核函数，N为数据采样点数，e为自然对数的底数。然后将变换得到的电流信号频谱I(f)除以电压信号频谱U(f)，得到变压器的传递函数H(f)，传递函数也是一组个数为N的数据，每个数据包括频率和对应的幅值，选取对应频率段为0～1MHz的传递函数数据。
2)、建立故障标准数据库，包括：试验变压器的型号参数，包括变压器额定电压，额定容量，变压器绕组结构和匝数；雷电冲击耐压试验的记录，包括采样频率、采样点数、变压器绕组首端电压和中性点电流数据；按照步骤1计算得到的0～1MHz频段内传递函数数据；故障信息，即最后通过检测确定的故障类型与故障位置。
3)、将待诊断的雷电冲击耐压试验的首端电压和中性点电流数据分别进行快速傅里叶变换，按步骤1求取0～1MHz频率段内的传递函数数据，通过50％雷电冲击电压与全压下的传递函数数据对比，如图2所述，两者有明显差别，判断发生了故障。
4)、通过对比判断存在故障，选取故障数据库中相同型号的10kV单相双绕组变压器的故障数据，按公式2将新试验的传递函数与故障数据库中的传递函数数据求取相关系数ρ，式中，x(n)为新的故障传递函数数据，y(n)为故障数据库中的传递函数数据，xi和yi分别为x(n)和y(n)中的第i项数据，和分别为x(n)和y(n)的平均值。设置阈值为0.9，当相关系数ρ≥0.9时，认为两者类似，将ρ≥T的故障类型按相关系数从大到小的顺序列表如表1。
表1
| 序号 | 相关系数 | 故障类型与故障位置 |
| 1 | 0.952 | 绝缘损坏造成的间隙放电故障，第151～158匝短路 |
| 2 | 0.948 | 绝缘损坏造成的间隙放电故障，第125～132匝短路 |
| 3 | 0.937 | 绝缘损坏造成的间隙放电故障，第159～166匝短路 |
| …… | …… | …… |
5)、按照排序的故障列表选取最接近的故障信息，设置偏差值5％，以故障对应位置为中心选取总匝数5％，总匝数为560匝，即选取28匝绕组，按序号1的故障位置选择141～168匝，利用超高频局放测试或超声局放测试仪对141～168匝绕组进行检测，诊断出故障位置为143～149匝的匝间故障。
Claims (1)
1.一种基于大数据的变压器雷电冲击的故障诊断方法，其特征在于包括以下步骤：
步骤1、收集变压器雷电冲击耐压试验的历史数据，按公式1对试验测量的首端电压数据和中性点电流数据分别进行快速傅里叶变换，式中x(n)为离散信号数据，如电压或是电流采样数据，X(k)为信号的快速傅里叶变换数据，为傅里叶变换的核函数，N为数据采样点数，e为自然对数的底数，k为各个采样数据的序号；然后将变换得到的电流信号频谱I(f)除以电压信号频谱U(f)，得到变压器的传递函数H(f)，传递函数也是一组个数为N的数据，每个数据包括频率和对应的幅值，选取对应频率段为0～1MHz的传递函数数据；
步骤2、建立故障标准数据库，包括：试验变压器的型号参数，包括变压器额定电压，额定容量，变压器绕组结构和匝数；雷电冲击耐压试验的记录，包括采样频率、采样点数、变压器绕组首端电压和中性点电流数据；按照步骤1计算得到的0～1MHz频段内传递函数数据；故障信息，即最后通过检测确定的故障类型与故障位置；
步骤3、将待诊断的雷电冲击耐压试验的首端电压和中性点电流数据分别进行快速傅里叶变换，按步骤1求取0～1MHz频率段内的传递函数数据，通过50％雷电冲击电压与全压下试验的传递函数数据对比，判断是否发生故障；
步骤4、待诊断的雷电冲击耐压试验数据通过传递函数对比判断存在故障，则选取故障数据库中相同型号变压器的相关数据，按公式2将新试验的传递函数与故障数据库中的传递函数数据求取相关系数ρ，式中，x(n)为新的故障传递函数数据，y(n)为故障数据库中的传递函数数据，xi和yi分别为x(n)和y(n)中的第i项数据，和分别为x(n)和y(n)的平均值；设置阈值T，当相关系数ρ≥T时，认为两者类似，将ρ≥T的故障类型按相关系数从大到小的顺序列表；
步骤5、按照排序的故障列表选取最接近的故障信息，设置偏差值η％，以故障对应位置为中心选取总匝数η％的绕组段为测试范围，利用超高频局放测试或超声局放测试进行故障定位诊断；若未检测出，则按列表依次选择下一个接近的故障信息，按上述方法选择测试范围进行测试，直到确定故障位置，完成诊断。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710142739.XA CN106950470B (zh) | 2017-03-10 | 2017-03-10 | 一种基于大数据的变压器雷电冲击的故障诊断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710142739.XA CN106950470B (zh) | 2017-03-10 | 2017-03-10 | 一种基于大数据的变压器雷电冲击的故障诊断方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106950470A CN106950470A (zh) | 2017-07-14 |
| CN106950470B true CN106950470B (zh) | 2019-03-29 |
Family
ID=59466843
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710142739.XA Active CN106950470B (zh) | 2017-03-10 | 2017-03-10 | 一种基于大数据的变压器雷电冲击的故障诊断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106950470B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108205662A (zh) * | 2017-12-29 | 2018-06-26 | 国网上海市电力公司 | 基于皮尔逊算法的输电线路雷电反击识别方法 |
| CN108267660B (zh) * | 2018-03-21 | 2020-02-14 | 广东电网有限责任公司电力科学研究院 | 一种变压器故障诊断方法及装置 |
| CN109870634A (zh) * | 2019-01-24 | 2019-06-11 | 深圳市新杰斯锐电子科技有限公司 | 一种变压器智能耐压检测系统 |
| CN110895322A (zh) * | 2019-10-31 | 2020-03-20 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | 一种固体绝缘电压互感器匝间短路诊断方法 |
| CN111583176B (zh) * | 2020-03-31 | 2023-09-15 | 北京全路通信信号研究设计院集团有限公司 | 一种基于图像的防雷模拟盘元件故障检测方法及系统 |
| CN112710929B (zh) * | 2020-12-16 | 2021-12-14 | 清华大学 | 一种电力电子变压器绕组层间绝缘监测方法及系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101251564B (zh) * | 2008-04-08 | 2010-06-09 | 昆明理工大学 | 一种利用可拓学与粗糙集理论相结合的电力变压器故障诊断方法 |
| CN103529337B (zh) * | 2013-10-30 | 2016-03-23 | 国家电网公司 | 设备故障与电气量信息间非线性相关关系的识别方法 |
| CN104198899B (zh) * | 2014-08-04 | 2017-01-25 | 西安交通大学 | 一种多放电源情况下变压器局部放电类型识别方法 |
| CN106338662A (zh) * | 2015-07-10 | 2017-01-18 | 华北电力大学（保定） | 基于数学形态学的变压器绕组变形诊断方法 |
| CN105807194A (zh) * | 2016-05-09 | 2016-07-27 | 三峡大学 | 一种基于lvq神经网络的雷电冲击下换流变压器绕组故障诊断方法 |
-
2017
- 2017-03-10 CN CN201710142739.XA patent/CN106950470B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106950470A (zh) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106950470B (zh) | 一种基于大数据的变压器雷电冲击的故障诊断方法 | |
| Samimi et al. | FRA interpretation using numerical indices: State-of-the-art | |
| Nurmanova et al. | A new transformer FRA measurement technique to reach smart interpretation for inter-disk faults | |
| CN105842596B (zh) | 一种高灵敏度电力电缆局部缺陷诊断方法 | |
| CN103792462B (zh) | 基于阻频曲线的电力变压器绕组匝间短路故障检测方法 | |
| CN109375060B (zh) | 一种配电网故障波形相似度计算方法 | |
| Bigdeli et al. | Detection of probability of occurrence, type and severity of faults in transformer using frequency response analysis based numerical indices | |
| CN106841949B (zh) | 三相异步交流电机定子绝缘在线监测方法和装置 | |
| WO2016065959A1 (zh) | 中性点不接地的10kV系统中铁磁谐振的诊断方法 | |
| CN108647479B (zh) | 一种避雷器故障暂态波形诊断方法及装置 | |
| CN109342883A (zh) | 一种用于电缆的局部老化故障检测定位方法 | |
| CN109521391B (zh) | 发电机电压互感器绕组匝间短路故障的检测装置及方法 | |
| CN104237713A (zh) | 基于离散小波变换的变压器绕组形变诊断方法 | |
| CN103454526A (zh) | 一种基于电力变压器绕组故障类型判定方法 | |
| CN110764027B (zh) | 基于频谱特征变化的电连接器发生间歇故障诊断方法 | |
| Nurmanova et al. | A new transformer fra test setup for advanced interpretation and winding short-circuit prediction | |
| CN102288881B (zh) | 变压器油纸绝缘尖刺放电缺陷严重程度诊断方法 | |
| CN110895322A (zh) | 一种固体绝缘电压互感器匝间短路诊断方法 | |
| Pant et al. | Detecting faults in power transformers using wavelet transform | |
| CN106526433B (zh) | 变压器电流波形特征试验方法和装置 | |
| CN110632481B (zh) | 一种中压电缆本体绝缘缺陷程度识别方法 | |
| CN110954818B (zh) | 一种中间继电器线圈软故障检测方法 | |
| Di Pasquale et al. | Frequency characterization of cast-resin transformers | |
| Župan et al. | Fleet Asset Management Opportunities Arising from Transient Monitoring of Power Transformers and Shunt Reactors | |
| CN113156275B (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 |