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

研究与开发年第期6电力系统故障录波数据分析邵玉槐许三宜何海祥丁周方(太原理工大学电气与动力工程学院，太原030024摘要电力系统故障录波数据是电力系统故障分析和保护动作判据的重要依据。本文提出了据电力系统故障录波数据完善了频率分析、谐波分析、故障定位的数学分析方法。采用java编程语言完成部分过程的编制工作。同时针对目前双端测距存在的伪根问题，提出了一种新的求解过程。关键词:电力故系统故障分析;故障录波数据;双端测距PowerSystemFaultRecorderDataAnalysisShaoYuhuaiXuSanyiHeHaixiangDingZhoufang(CollegeofElectricalandPowerEngineering,TaiyuanUniversityofTechnology,Taiyuan030024AbstractThepowersystemfaultrecorderdataprovidestheimportantbasisforfaultanalysisandprotectiveoperatingcriterion.Thepaperimprovedfrequencymeasurementmathematicalanalysisalgorithmandharmonicanalysismathematicalanalysisalgorithmaswellasfaultlocationmathematicalanalysisalgorithmbyuseofthosedata.Usingjavaprogramminglanguageasdevelopmenttoolsandaccomplishsomefunction.Atthesametime,thepaperproposesanewsolvingprocessaimingatfalserootsintwo-terminalfaultlocation.Keywords:powersystemfaultanalysis;faultrecorderdata;two-terminalfaultlocation1引言电力系统故障录波系统是电力系统发生故障及振荡时能自动记录的一种系统或一种装置。近年来，不同类型的故障录波器已在电力系统中得到广泛应用，所记录的各种故障录波数据为电力系统故障分析及各种保护动作行为的分析和评价提供了数据来源和依据。目前，电网调度端已能通过专用网或电话网将电网故障录波数据集中到一起，但如何有效管理和利用这些信息进行必要的故障分析、保护动作行为评价及故障测距等并没有统一的标准［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的强大功能可以生成了电压和电流波形图。具体的单通道波形生成程序为：山西省自然科学基金资助项目(6320097220001102电号技术研究与开发年第期图1JFreeChartchartA=ChartFactory.createXYLineChart("",//charttitle"X",//xaxislabel"Y",//yaxislabeldataset,//dataPlotOrientation.VERTICAL,true,//includelegendtrue,//tooltipsfalse//urls;importorg.jfree.data.general.DatasetChangeEvent;importorg.jfree.data.xy.AbstractXYDataset;importorg.jfree.data.xy.XYDataset;(2频率计算采样电压信号为一正弦电压，假设采样期间电压幅值保持不变，系统频率也不迅速改变，系统电压信号采样可用下式表示(sin(2nvtVft=设每周波的采样点数为N，每0.1s的采样点数为n,当电压信号以T为时间间隔采样时，第，1,2,3kkkk+++的采样值可以表示为sin(2nkvVft0=+1sin[2n(]kvVftT0+=++2sin[2n(2]kvVftT0+=++3sin[2n(3]kvVftT0+=++由上式可以推导出：2312cos(3n14sin(ncos(nkkkKvvfTfTvvfT++++==+(1为了消除120kkvv+++=的影响，利用等比定理若0acebdf==>厕aceabdfb++=++经过n-3次计算的214sin(nfT按上式取绝对值和进行补偿得到321312114sin(nnkkknkkkvvfTvv+=++=+=+ZZ(2今33131kkknkkkvvxvv+=++=+=+¥z(3根据采样值可以计算出每隔0.1s的x值由式(2和式(3可以推导出121narcsin[(1]2fTx=(4注:f=50Hz,T为采样周期，1TfN=根据式(4求出的x为x0,将121arcsin[(1]2x在x0处进行泰勒展开得21200033200022000523200011arcsin[(1]n(2223(23(1{[(23223(23](1}(2xfTxxxxxxxxxxxxxxxx=++++++(5可得ce5nfT=式([2](3全周傅氏变换算法这是在电力系统中应用很广泛的一种较好的算法，尤其是在电力系统微机保护提取基波分量时占有重要的地位。当电力系统发生故障时往往产生较大的衰减,非周期分量全周傅氏差分算法则可以消除这种误差。为了降低直流分量的影响采用了差分变换，即用采样值之差(1(xnxn+代替(xn。假定被采样信号是一个周期性的时间函数，除基波外还含有不衰减的直流分量和各次谐波，可表示为n，其为自然数，代表谐波次数。na和nb则分别代表各次谐波的余弦项和正弦项的振幅。1a和1b分别代表基波分量的余弦项和正弦项的振幅。结合全周傅式算法和差分采用矩形法可求得:112(sin(2n/NnkkkaxxknNNE2009727研究与开发年第期10NnkkkbxxknNN+==Z222/2/2[/2sin][/2cos]nsnsXanTbnTroro=+narctan((122nnsabfnT=(4双端测距利用双端数据的测距算法，在原理上完全不受故障过渡电阻(阻抗的大小、性质及双端系统阻抗的影响,可以保证测距的精度，但其主要缺点是需要通道传递对端信息。双端数据同步和消除测距方程伪根也是保证测距精度必不可少的手段。双端测距方法大多利用双端电压和电流量,列出从两端至故障点的输电厂线路的分布参数电压方程。对图1所示电路，可列出如下的电路方程：(((fFMFsMFcsssssUUchxIZshxyy=(6d(s((fFNFsNFcss((ssjUUchlxIZshlxd8YY=+(7(1(1(1NPMPIMPcs1UUchlIZshlYY=+(8(1(1(1MPNP1MP1csUIshlIchlZ丫丫=+(9式中,d8为两端数据的不同步相角差,s=0,1,2为序量标号广F”代表故障后,“P”代表故障前，((1/2cs(/ssllZZY=,“N”代表故障前。“f”代表故障点。图1单相接地故障由式(8或式(9得(1(1MP1MPcs1d(1NP(1(1argUchlIZshlUUshlIchlZI丫丫8丫丫=+=(10由式(6和式(7得dd(((MFNFsNFcss1(((sMFNFsNFcss1sjsssjssUeUchlIZshlxthIeUchlIZchl88yyyyy+=+(11对非对称故障，在式(1和(2中取s=1,2得d(1(1(1NF1NFcl1((jUchxIZshxUchlxIZshlxe8YYYY=+(12d(2(2MF1MFc11(2(2NF1NFc11((jUchxIZshxUchlxIZshlxe8YYYY=+(13将式(12和式(13相除，消去双端数据不同步相角差dje8(1(1MF1MFc11MFIMFc11(1(1NF1NFc11(2(2NF1NFc11(((((UchxIZshxDxUchxIZshxUchlxIZshlxUchlxIZshlx丫丫丫丫丫丫丫丫=++(14对于对称故障，可根据故障前的双端电压、电流列出如下方程d(1(1MP1MPc11(1(1NP1NPc11((]jUchxIZshxUchlxIZshlxe8YYYY=+(15将式(12和式(16相除，消去双端数据不同步相角差dje8(1(1MF1MFc112(1(1MP1MPc11(1(1NF1(1(1NP1NPc11(((((UchxIZshxDxUchxIZshxUchlxIZshlxUchlxIZshlxyyyyyyyy=++(16[3]其中，式(16适用于线路的全部短路类故障。式(14和式(16的求解过程繁琐，求出的结果往往是复数，存在的误差较大，为了剔除这些伪根，提出以下求解方法。牛顿法实质上是一种线性化方法,其基本思想是将非线性方法(Ofx=逐步归结为某种线性方程求解。它的优点是收敛快。令1(Dx和2(Dx为0,式(14和式(16变为(fx的函数，将函数(fx在点kx展开，有'((((kkkfxfxfxxx-+(17则1'((kkkkfxxxfx+=(0，1，2k=则下山法能保证函数值稳定下降，防止迭代发散，即具有单调性：1((kkfxfx+<(18结合二者的优点，令1'((kkkkfxxxfx+=(19与前一步的近似值kx适当加权平均作为新的改进值11(1kkkxxx久久++=+(202009728研究与开发

电区技7ft文件⑥囹形⑥分析㈤羞口曲TMCDL^_.止拯彼形图文件⑥囹形⑥M（y）分析㈤羞口曲TMCDm（aL^_.止拯彼形图年第期其中,01X<<,X称为下山因子。贝谥1'((kkkkfxxxfxX+=[4](0,1,2k=(21

(fx从X=1开始,逐次将X减半进行试算，直到能使下降条件(18成立为止。为了消除不对称故障时(14和(16的解出来的距离差距很大的取舍问题,同时根据模值平方相等建立2(((MFMFsMFcss2((NFsNFcss(((sssssfxUUchxIZshxUchlxIZshlx丫丫丫丫=+(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.⑵范岩，何勇.转子磁场定向速度自适应估算感应电机控制系统的研究[J].电气传动,2006,36(9:6-8.[3]李德华.交流调速控制系统[M].北京:电子工业出版社,2003.[4]MakoufA;DialloD;BenbouzidH,etal;Apracticalschemeforinductionmotorspeedsensorlessfield-orientedcontrol[J].EnergyConversion,IEEETransactionsonVolume19,Issue1,March2004:230-231.邱阿瑞，尹雁，王光辉，李永东.基于DSP的无速度传感器异步电机矢量控制系统[J].清华大学学报,2001,41(3:21-24.[6]黄永安马路刘慧敏.Matlab7.0/Simulink6.0建模仿真开发与高级工程应用[M].北京:清华大学出版社,2005.作者简介许琤(1981-,男，硕士，浙江舟山，研究方向:电力电子与电力传动。韩幸军(1972-,男,汉族,工程师，浙江舟山,研究方向:电机控制及风力发电技术。万山明，吴芳，黄声华.基于高频电压信号注入的永磁同步电机转子初始位置估计[J].中国电机工程学报,2008,28(33:82-86.陈荣.基于增量式光电编码盘的永磁同步电机转子位置初始定位[J].电机与控制应用,2007,34(3:32-34[7]徐永向，邹继明，赵猛,姚郁.基于单