电力系统故障录波数据分析.docx

电力系统故障录波数据分析.docx

《电力系统故障录波数据分析.docx》由会员分享，可在线阅读，更多相关《电力系统故障录波数据分析.docx（44页珍藏版）》请在冰豆网上搜索。

电力系统故障录波数据分析

。

研究与开发

年第期

6

电力系统故障录波数据分析

邵玉槐

许三宜

何海祥

丁周方

（太原理工大学电气与动力工程学院,太原030024

摘要电力系统故障录波数据是电力系统故障分析和保护动作判据的重要依据。

本文提出了据电力系统故障录波数据完善了频率分析、谐波分析、故障定位的数学分析方法。

采用java编程语言完成部分过程的编制工作。

同时针对目前双端测距存在的伪根问题,提出了一种新的求解过程。

关键词:

电力故系统故障分析;故障录波数据;双端测距

PowerSystemFaultRecorderDataAnalysis

ShaoYuhuaiXuSanyiHeHaixiangDingZhoufang

（CollegeofElectricalandPowerEngineering,TaiyuanUniversityofTechnology,Taiyuan030024AbstractThepowersystemfaultrecorderdataprovidestheimportantbasisforfaultanalysisandprotectiveoperatingcriterion.Thepaperimprovedfrequencymeasurementmathematicalanalysisalgorithmandharmonicanalysismathematical

精选资料，欢迎下载

。

analysisalgorithmaswellasfaultlocationmathematicalanalysisalgorithmbyuseofthosedata.Usingjavaprogramminglanguageasdevelopmenttoolsandaccomplishsomefunction.Atthesametime,thepaperproposesanewsolvingprocessaimingatfalserootsintwo-terminalfaultlocation.

Keywords:

powersystemfaultanalysis;faultrecorderdata;two-terminalfaultlocation

1引言

电力系统故障录波系统是电力系统发生故障及振荡时能自动记录的一种系统或一种装置。

近年来,不同类型的故障录波器已在电力系统中得到广泛应用,所

记录的各种故障录波数据为电力系统故障分析及各种保护动作行为的分析和评价提供了数据来源和依据。

目前,电网调度端已能通过专用网或电话网将电网故障录波数据集中到一起,

但如何有效管理和利用这些信息进行必要的故障分析、保护动作行为评价及故障测距等并没有统一的标准[1]。

2系统总体设计

java的最大优势就是跨平台,通俗地说可以用于各种操作系统,本系统是以java为平台开发的基于IEEE标准的COMTRADE数据格式的面向对象的可视化程序,下面简单说一下设计思路:

（1数据采用的格式

目前故障录波器基本上采用IEEE的COMTRADE标准。

每个COMTRADE记录都有一组4个与其相关的文件,其中CFG和DAT文件有严格的格式,用于存储通道数据和相关解释信息;HDR没有固定格式。

COMTRADE文件遵循固定的记录格

精选资料，欢迎下载

。

式,这使得编写程序读取数据成为可能。

文件以一行为单位记录录波信息,每行种又以逗号隔开各类信息或数据。

（2图形用户界面的显示

在本系统中,无论是波形分析（如波形的横向、纵向放大缩小;波形的瞬时值还

是故障分析（如谐波分析、序分量分析以及故障测距计算等结果,都要通过视图显

示到屏幕上,实现信息从机器到人的传递,因此,设计一个直观、友好的GUI

（GraphicalUserInterface,即图形用户界面对程序编写的基本要求[1]。

（3基本过程思路如图1。

3算法及程序

（1波形再现

现场采集的数据是一系列离散的点,由COMTRADE文件可以

计算出每个时刻对应的电压和电流瞬时值,利用jfreechart的强大功能可以生

成了电压和电流波形图。

具体的单通道波形生成程序为:

山西省自然科学基金资助项目（63

20097220001102

研究与开发

精选资料，欢迎下载

。

年第

期

图1

JFreeChartchartA=ChartFactory.createXYLineChart（

"",//charttitle"X",//xaxislabel"Y",//yaxislabeldataset,//data

PlotOrientation.VERTICAL,true,//includelegendtrue,//tooltipsfalse//urls;

import

org.jfree.data.general.DatasetChangeEvent;

importorg.jfree.data.xy.AbstractXYDataset;import

org.jfree.data.xy.XYDataset;（2频率计算

精选资料，欢迎下载

。

采样电压信号为一正弦电压,假设采样期间电压幅值保持不变,系统频率也不

迅速改变,系统电压信号采样可用下式表示

（sin（2π

vtVft=设每周波的采样点数为N,每0.1s的采样点数

为n,当电压信号以T为时间间隔采样时,第,1,2,3kkkk+++的

采样值可以表示为

sin（2π

kvVftθ=+1sin[2π（]kvVftTθ+=++2sin[2π（2]

kvVftTθ+=++3sin[2π（3]kvVftTθ+=++由上式可以推导出:

2312

cos（3π14sin（π

cos（π

kkkKvvfTfTvvfT++++=

=+（1

为了消除120kkvv+++=的影响,利用等比定理若0aceb

d

f

==>,则acea

bdf

精选资料，欢迎下载

。

b

++=++经过n-3次计算的214sin（πfT按上式取绝对值和进行补偿得到

3

32

1

3

1

2

1

14sin（πnk

kknkkkv

vfTv

v+=++=+=+∑∑（2

令

3

3

13

1

精选资料，欢迎下载

。

2

1

n

k

kknkkkv

vxv

v+=++=+=+∑∑（3

根据采样值可以计算出每隔0.1s的x值

由式（2和式（3可以推导出

121

πarcsin[（1]

2

fTx=（4

注:

f=50Hz,T为采样周期,1

TfN=根据式（4求出的x为x0,将1

21

arcsin[（1]

2

x在x0处进行泰勒展开得

精选资料，欢迎下载

。

2

2

2

2

1

2

00

0332

0002

2

00

0523

2

0001

1

arcsin[（1

]π（

2223

（23

精选资料，欢迎下载

。

（1

1

（

{[（23

2

23（23]（1}（2

xfT

xxxxxxxx

xxxxxxx

x=++++++（5

可得ce5πfT

=式（

[2]

（3全周傅氏变换算法

这是在电力系统中应用很广泛的一种较好的算法,尤其是在电力系统微机保护

提取基波分量时占有重要的地位。

当电力系统发生故障时往往产生较大的衰减,非周期分量全周傅氏差分算法则可以消除这种误差。

为了降低直流分量的影响采用了差分变换,即用采样值之差（1（xnxn+

代替（xn。

精选资料，欢迎下载

。

假定被采样信号是一个周期性的时间函数,除基波外还含有不衰减的直流分量和各次谐波,可表示为n,其为自然数,代表谐波次数。

na和nb则分别代表各次谐波的余弦项和正弦项的振幅。

1a和1b分别代表基波分量的余弦项和正弦项的振幅。

结合全周傅式算法和差分采用矩形法可求得:

1

1

2（sin（2π/

Nnkkkax

xknNN

+==

∑2009727

研究与开发

年第

期

精选资料，欢迎下载

。

1

10

2（cos（2π/

NnkkkbxxknNN+==

∑2

2

2

/2/2[/2sin][/2cos]

nsnsXanTbnTωω=+π

arctan（

（122

nnsabfnT=（4双端测距

利用双端数据的测距算法,在原理上完全不受故障过渡电阻（阻抗的大小、性质及双端系统阻抗的影响,可以保证测距的精度,但其主要缺点是需要通道传递对端信息。

双端数据同步和消除测距方程伪根也是保证测距精度必不可少的手段。

双端测距方法大多利用双端电压和电流量,列出从两端至故障点的输电厂线路的分布参数电压方程。

对图1所示电路,可列出如下的电路方程:

（（（fFMFsMFcsssssUUchxIZshxγγ=（6

d

（s（（

精选资料，欢迎下载

。

fFNFs

NFcss

（（ssjUUchlxIZshlxdδ

γγ=+（7（1（1（1

NPMP1MPcs1UUchlIZshl

γγ=+（8（1（1

（1MPNP

1MP1

cs

UI

shlIchlZγγ=+（9

式中,d

δ为两端数据的不同步相角差,s=0,1,2为序量标号,“F”代表故障后,“P”

代表故障前,

（（1/2cs（/ssllZZY=,“N”代表故障前。

“f”代表故障

点。

图1单相接地故障由式（8或式（9得

精选资料，欢迎下载

。

（1（1MP1MPcs1d（1NP

（1（1

MP1MP1cs

（1NP

arg

Uchl

IZshl

U

UshlIchlZIγγδγγ=+=

（10

由式（6和式（7得

dd

（（（

MFNFsNFcss

1（

（（

sMFNFsNFcss1sjsssjssUeUchlIZshlxthIeUchlIZchl

δδγγγγγ+=+（11

精选资料，欢迎下载

。

对非对称故障,在式（1和（2中取s=1,2得

d

（1（1

MF1MFc11

（1

（1NF1NFc1

1（（jUchxIZshxUchlxIZshlxe

δ

γγγγ=+（12

d

（2（2MF1MFc11（2（2

NF1NFc11（（jUchxIZshx

U

chlxI

Zshl

xe

δ

γγγγ=+（13

将式（12和式（13相除,消去双端数据不

精选资料，欢迎下载

。

同步相角差d

jeδ

（1（1MF1MFc11

1（2

（2

MF1MFc11

（1（1NF1NFc11

（2（2NF1NFc11（（（（（UchxIZshxDxUchxIZshxUchl

xIZshlxUchlxIZshlxγγγγγγγγ=++（14

对于对称故障,可根据故障前的双端电压、电

流列出如下方程

d

（1（1MP1MPc11（1（1

NP1NPc11（（]jUchxIZshxUchl

xIZshl

xeδ

γγγγ=+（15

将式（12和式（16相除,消去双端数据不

同步相角差d

精选资料，欢迎下载

。

jeδ

（1（1

MF1

MFc112（1（1MP1MPc11（1（1NF1

NFc11

（1（1

NP

1NP

c11（（（（（

UchxIZshxDxUchx

IZshxUchlxIZshlxU

chlxIZshl

xγγγγγγγγ=++（16[3]

其中,式（16适用于线路的全部短路类故障。

式（14和式（16的求解过程繁琐,

求出的结果往往是复数,存在的误差较大,为了剔除这些伪根,提出以下求解方

法。

牛顿法实质上是一种线性化方法,其基本思想是将非线性方法（0fx=逐步

归结为某种线性方程求解。

它的优点是收敛快。

令1（Dx和2（Dx为0,式

（14和式（16变为（fx的函数,将函数（fx在点kx展开,有

精选资料，欢迎下载

。

'（（（（kkkfxfxfxxx

≈+（17

则1'（（kkk

kfxxxfx+=（0,1,2k=

则下山法能保证函数值稳定下降,防止迭代发

散,即具有单调性:

1（（

kkfxfx+<（18结合二者的优点,令

1'（（

kkk

kfxxxfx+=（19

与前一步的近似值kx适当加权平均作为新的改进值

11（1kkkxxxλλ++=+（20

2009728

研究与开发

精选资料，欢迎下载

。

精选资料，欢迎下载

。

年第期

其中,01λ<≤,λ称为下山因子。

则式

1'（（

k

kk

kfxxxfxλ+=[4]

（0,1,2k=

（21

精选资料，欢迎下载

。

（fx从λ=1开始,逐次将λ减半进行试算,直到能使下降条件（18成立为

止。

为了消除不对称故障时（14和（16的解出来的距离差距很大的取舍问题,同时

根据模值平方相等建立

2

（（（

MFMFsMFcss2

（

（

NFs

NFcss

（（（

sssssfxUUchxIZshx

UchlxIZshlxγγγγ=+（22

作为附加判据。

选择令（fx为最小值的解。

4试验结果

针对现场采集的一组COMTRADE数据的分

析结果如下:

精选资料，欢迎下载

。

1,-45833,-592,944,-336,3,-567,981,-398,2,4,-1,4,3,-2,-45000,-

764,873,-93,2,-755,922,-152,3,5,-1,3,3,-3,-44166,-885,743,155,3,-

891,802,104,3,5,-1,2,3,-4,-43333,-947,561,394,3,-968,625,353,2,5,0,2,4,-

单通道波形结果如图2。

叠加通道波形结果如图3。

频率分析结果如图4

。

图2

单通道波形

图3

叠加通道波形

图4频率分析

5结论

电力系统自动化的重要程度日益提高,为电网调度自动化赋予了更重要的意义。

本文设计的基于COMTRADE故障录波信息的电力故障分析软件是应太原南瑞继保电力公司开发的一套软件,立足于实际条件,可实现波形再现、频率分析、谐波分析、

故障测距分析等功能。

目前,该软件已经完成了现场采集数据的部分析功能。

参考文献

[1]李媛,刘涤尘,杜新伟,王和春,叶菁.电力故障波形再现及分析系统的开发.电网技术,2006,30（5:

106-109.

（下转第33页

2009729

精选资料，欢迎下载

。

PE电力电子

年第期

33

图7

三相定子电流输出波形

7结论

仿真结果表明,转子磁链估算的波形基本保持不变,说明本文提出的改进定子

磁链观测器具有参数不敏感性,用该方法估算的定子磁链完全可以满足转速估算精度的要求。

从波形可以看出当给定转矩发生阶跃变化时,各参数的动态调节时间

较短,说明该系统具有较好的动态特性和较强的鲁棒性。

但是在实际系统中,电

精选资料，欢迎下载

。

机参数变化会影响系统性能,所以实现时还需要加上参数辨识和误差校正环节来

提高系统抗参数变化和干扰的鲁棒性[5-6]。

参考文献

[1]

李永东.交流电机数字控制系统[M].北京:

机械工业

出版社,2002.[2]范岩,何勇.转子磁场定向速度自适应估算感应电机控

制系统的研究[J].电气传动,2006,36（9:

6-8.[3]李德华.交流调速控制系统

[M].北京:

电子工业出版社,2003.[4]

MakoufA;DialloD;BenbouzidH,etal;Apracticalschemeforinductionmotorspeedsensorlessfield-orientedcontrol[J].EnergyConversion,IEEETransactionsonVolume19,Issue1,March2004:

230-231.

[5]邱阿瑞,尹雁,王光辉,李永东.基于DSP的无速度传感器异步电机矢量控制系统[J].清华大学学报,2001,41（3:

21-24.[6]

黄永安马路刘慧敏.Matlab7.0/Simulink6.0建模仿真开发与高级工程应

用[M].北京:

清华大学出版社,2005.

作者简介

许琤（1981-,男,硕士,浙江舟山,研究方向:

电力电子与电力传动。

韩幸军（1972-,男,汉族,工程师,浙江舟山,研究方向:

电机控制及风力发电技术。

（上接第25页

精选资料，欢迎下载

。

[5]万山明,吴芳,黄声华.基于高频电压信号注入的永磁同步电机转子初始位置估计[J].中国电机工程学报,2008,28（33:

82-86.

[6]陈荣.基于增量式光电编码盘的永磁同步电机转子位置初始定位[J].电

机与控制应用,2007,34（3:

32-34[7]

徐永向,邹继明,赵猛,姚郁.基于单霍尔传感器的PMSM位置检

测方法[J].哈尔滨工业大学学报,2008,40（11:

1766-1768.

[8]周卫平,吴正国,唐劲松,刘大明.SVPWM的等效算法及SVPWM与

SPWM的本质联系[J].中国电机工程学报,2006,26（2:

133-137.

[9]陈伯时.交流电机变频调速讲座.第四讲,电压空间矢量脉宽调制

（SVPWM控制[J].电力电子,2007（3.

[10]姜淑忠.空间矢量PWM算法的理解[J].伺服控制,2006（1:

38-41.

[11]祝琴,钟祥微,章鸿.SVPWM原理及其Simulink仿真

[J].电机电器技术,2005（4:

49-51.

[12]王兴,杨振强.电压空间矢量的原理及其在DSP上的实现[J].机电工程技术,2005（3:

16-18.

作者简介

许碧华（1973-,女,本科,工程师,从事纺织机器的电气设计工作。

（上接第29页