电流互感器暂态过程的仿真研究.docxVIP

电流互感器暂态过程的仿真研究 1 ct二次电流过节拍的保护 电力交换（ct）的暂时过程对电气系统的继电保护具有重要意义。现在，最常用的微型计算机保护的原理是保护电流差动。为了提高保护速度，许多采用了瞬时值差动的原则。因此，为了对ct的设计和性能进行更精确的临时模拟，有助于ct的设计和性能分析。本文在用7次多项式拟合铁芯基本磁化曲线的基础上,比较精确地对CT的暂态过程进行了仿真。 当CT通过较大的一次暂态电流(特别对于母线区外故障)时,其中会含有大量衰减的非周期分量,使得铁芯进入饱和区,CT二次电流将会出现较为严重的缺省,电流波形发生畸变,由此会造成很大的差电流,这样就容易引起继电器误动作。如何正确判断CT饱和的时刻,测定其线性传变区一直以来是此类保护的难点所在。 近年来提出了CT线性区开放母差保护的方法。由于在母线发生区外故障后,每周波前1/4周期CT铁芯不会发生饱和,可以将一个周波内二次电流波形第一次过零时刻与CT铁芯饱和时刻之间的时间作为CT的线性传变区。利用CT二次电流过零时的线性传变区投入差动保护,避开CT的饱和区。此类差动保护的关键技术是怎样精确检测铁芯饱和发生的时刻。由于当CT铁芯进入饱和时,CT二次电流波形会发生较大的畸变,出现奇异点(或突变点),近年来,从傅里叶变换发展而来的小波变换具有多尺度分析功能与良好的时频局部化特性,可以较为准确地捕捉突变信号的发生时刻。从而较好地检测出CT进入饱和的时刻。 本文基于小波变换的信号奇异性检测理论,对CT二次电流信号首先进行小波变换,以确定奇异点的位置,然后利用Hilbert变换提取信号包络,使得电流信号奇异点的定位更加精确,从而更加有效地防止由于CT铁芯饱和引起的电流瞬时值差动保护误动作。 2 ct临时模拟 2.1 硅钢片的拟合曲线 CT铁芯的基本磁化曲线有很强的非线性,选择适当的函数形式来拟合磁化曲线是十分必要的。本文中选用的铁芯材料是Z11-0.35型冷轧硅钢片,采用7次多项式对其进行了拟合,其磁化曲线及拟合曲线见图1。 拟合曲线与原磁化曲线各点误差的均值为0.8006,可见采用7次多项式已经可以取得很好的拟合效果,得到的拟合多项式方程为 H=f(B)=1253B7-6627B6+13761B5- 14100B4+7343B3-1802B2+191.0B-2 (1) 2.2 次电流cos 本文以一个保护用CT作为分析例子,其有关参数为:铁芯材料为Z11,一次绕组N1n=2,二次绕组N2n=120,二次负载为30 VA,功率因数为cosφ=0.8。正常时,二次电流为5 A,取最为严重的情况进行计算分析(即一次电流非周期衰减直流分量引起二次电流发生100%的偏移),一次时间常数为10 ms。本论文中假定故障均发生在0s时刻。图2(a)为正常时的励磁电流波形,图2(b)为二次电流波形。图3(a)为故障短路电流倍数为20时的励磁电流波形,图3(b)为二次电流波形。从图3中可以很明显地看出CT铁芯已经发生程度较深的饱和,二次电流波形明显畸变,出现了奇异点。 3 小波变换离散化 小波变换是指函数f(x)∈L2(R)与展缩小波的卷积,其实质是将f(x)分解为不同频带(s取不同值)的信号。连续小波变换定义为 Wsf(x)=f*Ψs(x)=1s∫∞?∞f(τ)Ψ(x?τs)dτ(2)Wsf(x)=f*Ψs(x)=1s∫-∞∞f(τ)Ψ(x-τs)dτ(2) 在实际应用中,为能在计算机上实现,通常要对小波变换离散化。取s=2j(j=1,2,…J),x=n(n=0,1,…,N-1),可得离散二进小波变换。为方便计算引入尺度函数,其展缩形式为 φ2j(x)=2-jφ(2-jx) (3) φ(x)满足 |φ?(ω)|2=∑j=1∞Ψ?(2Jω)χ?(2Jω)(4)|φ^(ω)|2=∑j=1∞Ψ^(2Jω)χ^(2Jω)(4) ∑j=?∞∞Ψ?(2Jω)χ?(2Jω)=1(5)∑j=-∞∞Ψ^(2Jω)χ^(2Jω)=1(5) 其中χ(x)为重构小波。 定义S2jf(x)=f*φ2j(x)为平滑算子,S2jf(n)=f*φ2j(n)为平滑系数序列,取S20f(n)为输入数字信号,则可得到离散二进小波变换的快速递推分解算法: S2jf(n)=∑k∈ZhkS2j?1f(n?2j?1k)W2jf(n)=∑k∈ZgkS2j?1f(n?2j?1k)(6)S2jf(n)=∑k∈ΖhkS2j-1f(n-2j-1k)W2jf(n)=∑k∈ΖgkS2j-1f(n-2j-1k)(6) 其中hk为低通滤波器系数,gk为高通滤波器系数。 4 信号包络提取方法 信号的突变信息往往体现在信号的包络里,而包络的波峰可以认为是信号变化最剧烈的时刻,在某些场合即是信号的突变点。用信号的包络来进行CT