非光滑雷电冲击波形参数的求取.docxVIP

非光滑雷电冲击波形参数的求取 0 可做平均曲线的试验测量 标准波形参数是根据光滑波形定义的。高电压试验中获得的波形往往在高电压测试中重叠噪声。因此，对于iec60-1中的非光滑波形参数，请参考以下内容。“对于几个测量器，当激发电压的最大值处有振动或激发电压，激发频率为0.5hz（lif）或激发持续时间为1s（liso）时，应进行平均曲线。在测量过程中，取互相重叠的频率s（lis）或s（ligo）的平均值作为试验电压，波前和波后的时间也应通过平均曲线计算。关于平均曲线的上述规定是以图例的方式说明的, 给用计算机自动求取波形参数带来的困难。90年代以来出现了多种获取平均曲线的算法。其中主要有拟合法、滤波法及小波变换方法。IEC及CIGRE90年代初为此成立专项工作。1996年IEC正式发布了IEC1083-2。该标准并不限制各实验室所用的求取波形参数的算法及软件, 而是给出了若干典型冲击波形 (以软盘形式) , 要求使用者处理相应波形 (参数设置, 如垂直分辨率、采样速率、内部噪声等, 应与使用者所用硬件相对应) 时所得波形参数应在标准所列范围内。 本文在文的基础上, 根据IEC1083-2, 改进并完善了双指数函数拟合算法。 1 算法的原理 1.1 采样间隔的拟合 在双指数模型下, 时间离散的雷电全波和操作波的标准形式为: A[e?αt(n?n0)?e?βt(n?n0)](1)A[e-αt(n-n0)-e-βt(n-n0)](1) 式中,α、β分别波尾、波前时间常数的倒数;A为幅值系数;t为采样间隔;n0为拟合后波形的原点。 对采样数据, 求得电压基值后, 经整理, 用最小二乘原理, 通过两步拟合, 可以分别求得参数A、α和β、n0, 得到拟合波形F[n]: {F[n]=0,F[n]=A[e?αt(n?n0)?e?βt(n?n0)],nn0n≥n0(2){F[n]=0,nn0F[n]=A[e-αt(n-n0)-e-βt(n-n0)],n≥n0(2) 解式 (2) 可得到其波形参数。 (1) um电压最大值 um=A(e?αtP?e?βtP)(3)um=A(e-αtΡ-e-βtΡ)(3) (2) 最大值时间tp和半最大值时间t2对si tP=ln (α/β) / (α-β) (4) t2=ln (0.5um/A) /α(5) (3) 求解非线性方程 t1=5 (t90-t30) /3 (6) t2=ln (0.5um/A) /α- [t30-0.5 (t90-t30) ] (7) 式中t30、t90分别是对应于波幅值30%、90%um的时刻, 可通过求解非线性方程求得。 1.2 剩余波形上振荡的lifo 定义数字记录仪的采样数据为原始波形, 双指数拟合算法应用于原始波形后得到的为拟合波形, 为了判断波形的类型将原始波形与拟合波表之差称为剩余波形。通过对原始波形、拟合波形、剩余波形及其相互关系的分析及大量实际波形的计算, 对雷电全波波形的类型判别给出以下准则: 光滑雷电全波波形较接近双指数波, 若剩余波形最大值与原始波形最大值之比2%, 则波为光滑雷电全波 (LI) , 见图1;否则, 认为波形带有过冲或振荡。 (图1～6中剩余波形均使用左坐标, 原始波形和拟合波形均使用右坐标) 。 对于仅在波前有振荡的波 (此振荡在接近波峰时应已基本消失) , 其振荡频率较高, 在剩余波形上原始波形峰值时间之前通常有3个以上的振荡峰值。据此, 可判断出LIFO, 例见图2。 IEC60-1以图例的形式给出的波顶带过冲和波带振荡的冲击波有很明确的区别, 前者叠加的振荡衰减很快, 因而出现的半周波个数不超过两个;而后者波形上出现的叠加振荡的半周波个数多于3个。但IEC1083-2提供的这两类波形并无这种特征, 其区别为前者由解析式通过计算机模拟产生, 而后者由实际试验回路产生。其关键是判断出其叠加振荡的频率或峰值左右过冲的持续时间。因此, 若剩余波形上出现第3个峰值的时间在原始波形峰值时间之后, 可判定为波顶带振荡的冲击波;若剩余波形上的振荡数3, 且振荡峰值发生在原始波形峰值左右, 可判断为波顶带过冲的冲击波 (IEC1083-2未给出此种波形) 。 波峰叠加过冲波形的过冲持续时间τ定义为与原始波形最大值处对应的剩余波形波峰两侧的两个过零点之间的时间段。若τ1μs, 判定为LILO, 见图3;否则, 判定为LISO, 见图4。 波峰叠加振荡波形的振荡频率f, 定义为与原始波形最大值处对应的剩余波形波峰与其下一个波峰之间时间长度的倒数。若f0.5 MHz, 则判定为LIFA, 见图5;否则, 判定为LISL, 见图6。 操作波的判别很简单, 根据双指数拟合算法得到的波前时间t110 μs, 可判定是SI。