雷电冲击电压波形K因子零相位数字滤波器设计
2021-03-16
来源:爱问旅游网
第31卷第2期 2012年4月 电工电能新技术 Advanced Technology of Electrical Engineering and Energy Vo1.31,No.2 Apr.2012 雷电冲击电压波形 因子零相位数字滤波器设计 司文荣 ,傅晨钊 ,黄 华 ,金 珩 ,李彦明 (1.上海市电力公司电力科学研究院,上海200437; 2.电力设备电气绝缘国家重点实验室(西安交通大学),陕西西安710049) 摘要:针对最新版IEC60060—1《高电压试验技术,第一部分:一般定义及试验要求》中叠加过冲和 振荡的标准雷电冲击波形的数字处理程序,其引入与频率相关的试验电压因数Ii}( ,介绍一种设 计等价于尼(. 对应转移函数的 因子零相位IIR数字滤波器的方法。对零相位数字滤波原理进 行了分析与仿真实验。基于IEC60060—1中定义的试验电压因数k(.厂),给出了设计K因子零相位 IIR数字滤波器的推导过程,以及可用于开发冲击电压测量软件的时域滤波差分方程。最后对最 新版IEC61083-2((高压冲击试验测量用仪器和软件,第二部分:对软件的要求》附带的TDG 2.04中 几种典型雷电冲击电压波形进行了滤波处理与参数提取,验证了该方法的可行性和准确性。为研 制冲击电压测量系统提供了一种实用的数字滤波技术。 关键词:试验电压因数 ( ;零相位;数字滤波;雷电冲击;参数提取 中图分类号:TM835 文献标识码:A 文章编号:1003-3076(2012)02-0011-04 1 引言 电力系统在实际运行中不可避免地会遭受雷电 冲击,为了研究电力设备在遭受雷电过电压时的绝 缘性能,许多电气设备在型式试验、出厂试验或大修 后利用冲击试验发生器进行雷电冲击电压试验。雷 后的数据;④对剩余数据用双指数函数模型拟合: “d(£)=A[exp(一(t—D)/8)一 exp(一(t—D)/c)] (1) 式中,t是时间;u (t)是双指数电压函数;A,B,C 和D是拟合参数。⑤将U (t)中点D至点 (£)中减去 之间 的数据替换为u (t),得基准曲线U (£);⑥从 (t)以获得剩余曲线U (t)= 电冲击电压的破坏作用不仅决定于幅值,还与波前 陡度有关。为了保证多次试验结果的重复性和各实 验室间试验结果的可比性,对波形类型及波形参数 定义应有明确规定,为此国际电工委员会(IEC)和 国家标准(GB)规定了标准雷电冲击全波和截波的 (t)一 (t);⑦构建一个数字滤波器,其转移函 数为试验电压因子函数Jc( ,其与频率相关(见图 1): ( =1/(1+ ) (2) 具体特征参数…。 最新出版的IEC60060-1给出的叠加过冲和振 式中,口=2.2,f为频率,单位为MHz。⑧对剩余曲 线U ( )使用数字滤波器获得滤波后的曲线 荡的标准雷电冲击参数计算程序如下 引:①从输 入电压为零所记录的开始部分计算电压的平均值, U (t);⑨将曲线U ,(t)叠加到U (t)上获得试验电 压曲线 .(t),从而进行参数提取(包括试验电压值 求取记录曲线的基准水平;②从记录曲线 (£)中 去掉基准水平偏置,并用该曲线 (t)进行后续步 u 、波前时间 .、半峰值时间 和相对过冲幅值 。 与以往旧版本IEC60060-1标准中给出的冲击 电压参数计算程序比较可得出,其最大的改动在于 骤,求取其极值U ;③求取波前小于0.2U 电压值 和波尾大于0.4U 电压值的两个采样点 :(最后 一对剩余曲线 (f)进行了 因子滤波。因此针对 式(2)定义的试验电压因数Jj}( ,如何设计一个准 个点)和U0(第一个点),删除点£, 前和点 . . 收稿日期:2011-03-21 基金项目:国家自然科学基金资助项目(50377034) 作者简介:司文荣(1981-),男,江苏籍,工程师,博士,从事高电压试验与设备诊断技术研究; 金珩(1961一),男,山东籍,首席工程师,长期从事电力设备检测、故障诊断及状态评估技术研究。 l2 电工电能新技术 第31卷 } i i;l f 。 i -…一 ”^}¨ m~ ・¨・ …-・{…・ I_l 越 搽 圜 出 }整 。‘ I暮…. --¨●.I}^¨|一 f *…“÷…●” 1…・ …鹭 图1试验电压因数 ( Fig.1 ( factor characteristic 确且可以用于冲击电压波形参数提取的数字滤波 器,或给工程技术人员提供开发冲击电压测量软件 的时域滤波差分方程,是目前需要解决的关键问题。 笔者则介绍使用零相位数字滤波器技术,基于定义 的试验电压因数 ( ,给出如何设计用于冲击电压 参数提取用的K因子零相位数字滤波器的方法。 2零相位数字滤波器 2.1方法与原理 零相位数字滤波器的实现方法可简述如下 : ①先将输入序列按顺序滤波(forward filter);②再将 所得结果逆转后反向通过滤波器(reverse filter);③ 最后将所得结果逆转后输出(reverse output),即得 到精确无相位失真的输出序列。其滤波原理可由如 下的时域和频域描述解释。零相位数字滤波的时域 描述: Y (n)= (n){l:h(n) (3) Y (n)=Y (N一1一rt) (4) Y (n)=Y (rt);l:h(n) (5) Yd(n)=Y。(^ 一1一n) (6) 式中, (n)表示输人序列; (n)为所用数字滤波器 冲激响应序列;Y (n)为零相位滤波的最终输出结 果。上述时域描述对应的频域描述: (e ): (e )¥ (e ) (7) Yb(e )=e一 ‘ 。 y日(e一 ) (8) (e )=Yb(e ) H(e ) (9) yd(n)=e ‘ ’l,c(e一扣) (10) 则可得: yd(e )=I H(e )I X(e ) (11) 由式(11)可见,输出 (e )与输入X(e )之 间为零相位关系,但幅频特性关系为1日(e‘j )I 。 2.2仿真实验 为验证前面的论述,笔者做了如下对比实验。 采用5阶契比雪夫(II)低通滤波器 ,数字截止频 率为100Hz,50dB衰减。对采样率1000Hz下获得 的叠加随机噪声的50Hz正弦波采样获得的波形, 进行零相位数字滤波,结果如图2所示。该结果表 明原始信号与滤波后信号在相位上为零误差。 图2零相位数字滤波器仿真实验 Fig.2 Simulation test of zero—phase digital filtering 3 因子零相位数字滤波器设计 为求取试验电压因数,必须建立其幅频响应等 于试验电压因数函数的数字滤波器。 由式(2)定义的幅频特性及∞=2 ,定义: K( )=一20log1o(1+口 /4,rr ) (12) 考虑到零相位滤波器的工作方式及式(11),滤 波先正向后反向,因此滤波器的衰减仅为所需的一 半,式(12)可写成: K( )/2=一10logl0(1+ao) /4"n" ) (13) 由低通滤波器,在连续域定义G(s): G(s)=1/(s/ +1) (14) 式中,s为拉普拉斯算子(Laplace operator);∞ 为截 止频率。 于是,在截止频率处滤波器的幅值可以定义为: K(∞)/2=一10log10(2) (15) 式(13)和式(15)可得定义的截止频率为 ∞ =21T√1/a (16) 利用双线性变换: s= T 力 = L¨ 式中, 为采样间隔;F=tan(历 1'/2)。 笙 塑 ————————————~司文荣,等: :二:::::!! 重皇 童皇 ,望 ! 位数字滤波器设计 l3 把式(17)代入式(14)可得: ’ f南 1 一 I-F)(18) 式(8)可写成: 一, ‘ ) (南 式中, _ta (订 )/(】一 一-) ( 9) 二 图3 LI—M2参数提取结果 南 ¨ )+ , 1) Fig・3 Par ̄lnleters evaluati。n f0r LI.M2 Jh-Jl ̄i T g ̄lOOMS/s采样获取的冲击电压波形=lO.n 以得 雨平 …由用于数字滤波): 下的差 … —:~ _。一 々^ 且拔 ~-1/\f ~i一 } 。、一 … …。。。。 ~ ÷---一一{一 ~ 1 ~ {— ~I--y‘ 。 。02074433774260(x(i)+ ( : -.. 其中甘由,y(一1):0 ,、.o_95851132514810y( 一1) 1))+ 一(2… {一一一一{~ ‘1) c i i…。。 ;三 一_ J … a原始波形 b基准水平 c记录曲线 d基准曲线 e剩余曲线 f剩余滤波曲线 g试验电压 ’ 、一。一 』一f …一 C十 警 取试验电压因数的数字篓 篓 .与最新版IEC60060.1附录c(资 滤波 出的相同 …… 刚 蒜 陆 图4 u-M5参数提取结果 n昏4 Parameters evaluati。n f0r LI—M5 耋芏妻 譬波形冬项参数的数字计算软件可 兰 软 要 编制 。为检 j;I IEC q ̄~田 —结果的准确性. :ex ̄-ir;ig ,曼 c6 3-2《高压冲击试验测量 茇 软件的要求》。 萎 软件,称作试 :兰毒 。 。 1996 年颁布了图5 LI—M6参数提取结果 Fjg‘5 Parameters evaluati。n f0r LIM6 —一 曼 有多种冲击 操作波论 竺拟实际测量; ,;其中有光嘉 茗暑 荡的;还有 娄 . r; : 薹 苎 标准还规定手 这些冲击波形已 允许的偏差范围。 ≯ … 0一一---0 …~ 一. ., a原始被形 ’。。。 一 ~ …~~~-、 一--一~0… .0 … 一b基准承平 c记录曲线 ■一- ~0.~… d基准曲线 e剩余曲线 ...____ 、~ 曲线 曼篓 篓据,耋 苎 被检的数字计 荐 数在标准 一 、一I、 -…一 .~~I… ~ 图6 LI—M7参数提取结果 匮测量中使用 软件才是合格的;一~ 伍 田否 喜 . Fig・6 Paraine£er8 evalua f0r uM7 .'D 0 4 型实测雷电冲击波形参数提取 TDG2曼 三!尊波形曲线(等等)的定义原始波形、基准水 乏 苎 中4个 妻 给予过多解释见文献[ 蚕 (因子零相位滤波的结果。其中 曩 (软件中按式( 1 平 ~ ’尥土 、骨 。 表1对4个典型波形的参数提取结果与标准给 14 电工电能新技术 第31卷 出的参考值进行了对比。表中参数 、 和 的 误差计算定义为: i=(P 一P )/P x 100% (22) 式中,P。为参数的计算值;P 为参数的标准值。 而参数口 的误差计算定义为: 2=P 一P (23) 表1 雷电冲击电压波形参数计算对比 Tab.1 Comparison of calculated parameters for lightning impulse 从表1数据可以看出对剩余曲线进行K因子 零相位滤波后建立的试验电压曲线,再求取的波形 参数,其与标准给出的误差均在可接受的范围内。 5 结论 本文给出了用于最新版IEC60060-1《高电压试 验技术,第一部分:一般定义及试验要求》中叠加过 冲和振荡的标准雷电冲击电压的数字处理程序,其 引入与频率相关的试验电压因数k( 对应转移函 数的K因子零相位IIR数字滤波器的设计方法及其 推导过程,并给出了可用于开发冲击电压测量软件 的时域滤波差分方程。 对最新版IEC61083-2(高压冲击试验测量用仪 器和软件,第二部分:对软件的要求》附带的TDG 2.04中几种典型雷电冲击电压波形进行了滤波处 理与参数提取,参数值均在最新标准规定的允许误 差范围内。该试验结果表明该 因子零相位数字 滤波器可以用于研制新型冲击电压测量系统。 参考文献(References): [1]郭弘,司文荣,李彦明,等(Guo Hong,Si Wenrong,Li Yanming,et a1.).带过冲和振荡雷电冲击波形的参数 提取研究(Parameter extraction of lightning impulse with overshoot or oscillation)[J].高压电器(High Voltage Apparatus),2010,46(I):76-79. [2]IEC 60060-1 Ed.3.0,High・voltage test techniques Part 1:General deifnitions and test requirements,2008[s]. [3]Paul L Lewin,Trung N Tran,David J Swafield.Zero— phase filtering for lightning impulse evaluation:a K—factor iflter for the revision of IEC 60060.1 and-2『J].IEEE Trans.on Power Delivery,2008,23(1):3—12. [4]纪跃波,秦树人,汤宝平(Ji Yuebo,Qin Shuren,Tang Baoping).零相位数字滤波器(Digital filtering with ze— ro phase error)[J].重庆大学学报(Journal of Chongqing Univ.),2000,23(6):4-7. [5]郑重,谈克雄,王猛,等(Zheng Zhong,Tan Kexiong, Wang Meng,et a1.)基于脉冲波形时域特征的局部放 电识别(Partial discharge recognition based on time do— main features of pulse waveform)[J].电工电能新技术 (Adv.Tech.of Elec.Eng.&Energy),2001,20(1): 20-24. [6]司文荣,李军浩,李彦明,等(Si Wenrong,Li Junhao, Li Yanming,et a1.).基于宽带检测的局放脉冲波形快 速特征提取技术(Fast feature extraction technique for PD pulse shape based on wideband detection)[J].电工 电能新技术(Adv.Tech.of Elec.Eng.&Energy), 2oo8,27(2):11.15. [7]DL/T 992—2006,冲击电压测量实施细则(Detailed im— plementation guide for impulse voltage measurement) [s]. 『8]IEC 61083.2 Ed.2.0,Instruments and software used for measurement in high・voltage impulse tests Part 2:Re- quirement ofr software,2008[S]. (下转第46页,cont.on P.46) 电工电能新技术 第31卷 Research on modulation method of five・leg inverter fed dual induction motor drive system LIU Jia。,MEI Yang ,WANG Li.na。SUN Kai。 ,(1.School of Automation Science and Elec.Eng.,Beihang University,Beijing 100191,China; 2.Power Electronics&Motor Driving Eng.Center,North China University of Technology。 Beijing 100144,China;3.State Key Lab of Power Systems,Dept.of Elec.Eng., Tsinghua University,Beijing 100084,China) Abstract:In electric trains,electric vehicles and some industry applications,two or more motors are needed to op- crate in paralle1.Therefore。multiple inverters are needed in these drive systems.In order to reduce the system cost,a dual induction motor drive system fed by a five-leg inverter with independent control is investigated in this paper.An improved 7-section modulation method,which is based on the traditional Space Vector PWM(SVPWM) technology and which divides the control period into two parts,is proposed to control the five—leg inverter.This method is suitable to be implemented by DSP,which is its main advantage.Through simulation and experiments, the feasibility and effectiveness of the proposed method are proved. Key words:five・leg inverter;modulation method;induction motor;SVPWM (上接第l4页。cont.from P.14) Design of zero-phase digital iflter with K factor used in parameters evaluation for lightning impulse sI Wen.rong ,FU Chen.zhao ,HUANG Hua ,JIN Heng ,LI Yan.ming (1.Electric Power Research Institute,SMEPC,Shanghai 200437,China; 2.State Key Lab of Elec.Insulation and Power Equipment(Xi’an Jiaotong University), Xi’an 710049,China) Abstract:According to the procedures for calculation of parameters of standard lightning impulses with superim- posed overshoot or oscillations in IEC 60060—1 Ed.3.0:High—voltage test techniques Part 1:General definitions and test requirements,and the introduced ̄equency dependent test voltage factor k( ,a way to design of zero— phase digital filter corresponding to the ( with IIR method is given in this paper.Principle of the zero-phase dig— ital filter is shown with simulation test using signal processing to sine wave.The difference equation with the same  ̄equency characteristic as Ij}( defined in IEC 60060-1 Ed.3.0,used to develop the software for measurement in impulse test is given.At last,the method is validated by processing several typical lightning impulses in TDG2.04 of IEC 61083-2 Ed.2.0:Instruments and software used for measurement in.high-voltage impulse tests Part 2:re- quirements for software,which shows its feasibility and veracity.This provides an effective digital filtering technolo— gY to develop impulse test system. Key words:test factor k( ;zero—phase;digital filtering;lightning impulse;parameters evaluation