冯晓丽, 王智勇, 杨紫明, 肖 毅, 张宏波, 方志朋
(1. 辽宁工程技术大学电气与控制工程学院, 辽宁 葫芦岛 125105; 2. 国家电网公司葫芦岛供电公司,辽宁 葫芦岛 125000; 3. 国家电网公司北辰新汇公司, 辽宁 沈阳 110013)
滑动电接触广泛应用于高速电气化铁路、航空航天等领域。在操作过程中,接触不良会产生电弧而引发强烈的电磁噪声,其中辐射电磁噪声对通信系统产生干扰,严重时会导致设备无法正常运行。电弧辐射电磁噪声具有持续时间短、突变快等特点,其频段分布因外界环境的影响而不同,对测试接受设备要求高,给现场测试造成诸多困难。因此,研究电弧辐射电磁噪声的频率特性具有重要意义。
国内外学者对电弧辐射电磁噪声频率范围进行了大量研究并获得了可观的结论。Keiichi Uchimura[1]研究了稳定燃烧电弧和喷弧的辐射电磁噪声的频段及特性。郑宁敏等[2]获得了刀闸操作电弧的幅频特性,并指出电弧辐射的能量集中分布在50~150kHz与300~350kHz之间。蓝会立等[3]分析了开关柜内故障电弧的频谱特性并指出故障电弧的频率主要分布在5~10kHz之间。Pan Tao等[4]认为矿用直流电机车弓网电弧的辐射噪声的频率低于1900MHz。郭凤仪等[5,6]获得了弓网电弧辐射电磁噪声幅值与滑动速度、接触电流、接触压力等因素之间的关系,并指出弓网电弧主要分布在30~500MHz之间。马云双等[7]对不同车速条件下弓网离线电弧放电电磁骚扰发生的幅值及频率进行了仿真研究。王莉等[8]研究了航空交流故障电弧特性,并指出其高频分量分布在10~50kHz之间。郭云梅[9]分析了航空直流电弧时频特性并获得了航空直流电弧的频率在2kHz以内的结论。
现有的电弧辐射电磁噪声的频谱特性大多通过干扰接收机获得,干扰接收机的工作方式为扫频模式,完成一次测试需要一定的时间,因而不能实时地获得电磁噪声的频谱信号,获得电弧辐射电磁噪声的瞬时频率更为困难。
为解决电弧辐射电磁噪声瞬时频谱特性难于测量的问题,并获得电磁噪声的瞬时频率,本文提出了一种基于Hilbert-Huang变换(HHT)的电弧辐射电磁噪声瞬时频率分析方法,并对其有效性进行实验验证。
利用自行研制的电弧辐射电磁噪声实验系统开展实验,该系统由电波暗室、电弧放电回路和信号检测系统组成,如图1所示。
图1 电弧辐射电磁噪声实验系统Fig.1 Electromagnetic noise experimental system
电弧放电回路用于产生特定电流的单次电弧,由调压回路和电弧发生器组成,如图2所示。
图2 电弧放电回路Fig.2 Arc discharge circuit
电弧发生器由平板状的浸铜碳电极、运动的铜电极和电极运动控制系统组成,如图1(a)所示。电极运动控制系统由步进电机、滑台和电机控制器组成。通过控制步进电机使安装在滑台上的铜电极实现往复运动,与静止的碳电极配合获得电弧放电。为了在每个运动周期获得一系列的单次电弧并提高实验效率,在碳电极上刻制了特定参数的V型槽,如图3所示。
图3 碳电极V槽Fig.3 V-groove structure in carbon electrode
信号检测系统由HFH2-Z2 型环形天线、HL562-Z1 型对数周期天线和Tektronix MSO4104B型示波器组成,如图1(b)所示。环形天线和对数天线分别检测电弧辐射的高频磁场信号和电场信号,环形天线的测量频段为9kHz~30MHz,对数周期天线的测量频段为30MHz~3GHz。天线的信号通过屏蔽电缆连接至示波器。示波器的采样率为5GS/s。
碳电极的尺寸为 250mm×25mm×10mm,铜电极的截面积为120mm2。在20℃时其物理参数如表1所示。按照表2的实验条件进行了单次电弧辐射电磁噪声实验。
表1 碳电极和铜电极的物理特性Tab.1 Phycical parameters of slide and contact wire
表2 实验条件Tab.2 Experimental conditions
HHT是一种较新的时频分析技术,由经验模态分解(EMD)和Hilbert变换两部分组成,其中EMD为核心部分,其功能是对信号进行分解;Hilbert变换求解信号的Hilbert谱。
经验模态分解方法针对数据固有的时间尺度进行信号分解,是一种自适应的信号分析方法,其目的是为了提取固有模态函数。利用EMD求取信号的固有模态函数主要包括三个步骤[10]:
(1)取出信号x(t)的各个局部极值,包括极大值和极小值,然后利用三阶样条函数对局部极值进行插值,分别获得信号x(t)的上包络序列xmax(t)和下包络序列xmin(t)。
(2)对每个时刻的xmax(t)和xmin(t)求平均值,得到上下包络的平均值m(t):
(1)
(3)用信号x(t)减去m(t),得到新数据列h(t):
h(t)=x(t)-m(t)
(2)
判断h(t)是否为固有模态函数,若不满足定义,重复步骤(1)~步骤(3),直至求出固有模态函数为止。
求出第一个固有模态函数I1(t),即从信号中分解出了第一个分量,然后用信号x(t)减去I1(t),得到余项r1(t):
r1(t)=x(t)-I1(t)
(3)
至此,第一个固有模态函数提取完成。然后对r1(t)重复步骤(1)~步骤(3)提取n个固有模态函数,直至rn(t)单调。
信号x(t)被分解为多个固有模态函数和一个余项:
(4)
基于Hilbert变换的瞬时频率定义为:
(5)
式中
H[x(t)]为信号x(t)的Hilbert变换。
对于单个固有模态函数Ii(t),通过Hilbert变换获得其瞬时频率fi(t):
(6)
信号的瞬时频率表征了信号的频率随时间的变化,通过HHT获得的解析信号与原信号的频谱完全相同[11],且不会出现没有物理意义的负频率。
测量距离为1m、接触电流为10A条件下的电弧辐射电磁噪声波形如图4所示。单次电弧辐射电磁噪声时域波形是一个短暂的阻尼振荡波。在不同测量距离、不同接触电流下,电弧辐射电磁噪声幅值有所不同,但波形大体轮廓基本保持一致。
图4 电弧辐射电磁噪声时域波形Fig.4 Time domain waveform of electromagnetic radiation noise produced by single arc
电弧辐射电磁噪声经过EMD分解后,获得了一系列固有模态函数(IMF)分量和一个余项。以图4(b)的电场噪声为例,其EMD分解结果如图5所示。
图5 电场噪声的EMD分解Fig.5 EMD results of electromagnetic noise
EMD将辐射电磁噪声分解成若干IMF分量,各个IMF反映了电磁噪声不同的频段范围,因此需要对IMF进行筛选,根据IMF与电磁噪声的相关系数及其方差贡献率进行选择。
相关系数及方差贡献率的定义为:
(7)
(8)
以图4(b)的电场噪声为例,计算后的相关系数及方差贡献率如表3所示。
表3 IMF分量的ρ和KiTab.3 ρ and Ki of IMF components
选择ρ和Ki值最大的IMF4,根据式(6)计算IMF4的瞬时频率,如图6所示。图6中,IMF的瞬时频率高达2.5GHz,而示波器的采样频率为5GHz,说明2.5GHz的瞬时频率不一定具有参考价值。出现如此高的频率是因为计算误差被放大[12],因此对瞬时频率进行平滑处理,选择Matlab中一维数字滤波器filter,处理后效果如图7所示。
图6 IMF的瞬时频率Fig.6 Instantaneous frequency of IMF
图7 平滑处理后IMF的瞬时频率Fig.7 Instantaneous frequency of smoothing-processed IMF
一维数字滤波器filter按式(9)计算:
(9)
式中,Y(z)为输出信号;a、b为滤波器参数,其值分别为0.15和0.8;X(z)为输入信号。
平滑处理后,IMF的瞬时频率的畸变成分被滤除,更能表征电磁噪声信号的频率随时间的变化趋势。图7中,IMF4的瞬时频率的均值为236.51MHz,不同测量距离、不同接触电流条件对瞬时频率有一定影响。
为验证本文方法的有效性,对EMD分解后的各个IMF分别求取功率谱密度(PSD)。功率谱密度表征了信号的能量与其频率之间的关系,对于离散随机信号x(n)其功率谱密度为[13]:
(10)
式中,rxx(k)=E[x(n)x*(n-k)]。
功率谱密度的计算方法有直接法、间接法和Welch法等。Welch法计算后的PSD谱曲线更为光滑,避免当数据较长时,直接法计算的谱曲线起伏剧烈的缺点,同时提高了PSD谱的分辨率。采用Welch法计算电弧辐射电磁噪声的功率谱密度。以图4为例,计算结果如图8所示。
图8 电弧辐射电磁噪声功率谱密度Fig.8 PSD of electromagnetic radiation noise produced by single arc
EMD分解的IMF分量分别表征了电磁噪声低频和高频部分。磁场噪声和电场噪声在低频部分幅度较高,因而在低频部分能量较高。分别对各IMF分量求PSD并对比发现,以电场噪声为例,IMF1分量的PSD是电磁噪声最高频部分,IMF4的PSD是电磁噪声的低频部分,筛选后IMF分量与电磁噪声的PSD对比如图9所示。
图9 IMF分量与电场噪声的PSD对比Fig.9 Comparison results between IMF component and PSD of electric field noise
图9中,IMF1在高频部分与电磁噪声的PSD峰值点基本一致;IMF4在低频部分与电磁噪声的PSD峰值点基本一致,说明本文提出的HHT方法分解的IMF可以表征电磁噪声不同频段的特性[14]。由于电磁噪声低频部分PSD幅度远高于高频部分,说明电磁噪声低频部分能量分布较为密集,故选择能表征电磁噪声能量分布密集的IMF分量进行研究更具参考价值。
均方根表征了随机数据的一般强度,方差表征了信号的变化程度。不同测量距离、不同接触电流条件下,按照图10所示流程对电磁噪声进行处理。
图10 电磁噪声分析流程图Fig.10 Flowchart of analyzing electromagnetic noise
对每组实验条件下电弧辐射电磁噪声的瞬时频率统计了均方根(RMS)和方差(VAR)并对统计值取均值,获得了不同工况下,电弧辐射电磁噪声瞬时频率的统计图,如图11所示。
图11 不同工况下电磁噪声瞬时频率的均方根和方差Fig.11 RSM and VAR of instantaneous frequency of electromagnetic noise under different experimental conditions
图11(a)和图11(b)中,磁场噪声瞬时频率的均方根在60MHz附近,随着接触电流、测量距离的增加,磁场噪声瞬时频率的均方根稍有增加。电场噪声瞬时频率的均方根在96~165MHz,随着接触电流、测量距离的增加,电场噪声瞬时频率的均方根明显减小。
图11(c)和图11(d)中,随着实验条件的改变,磁场噪声的瞬时频率的方差轻微波动,说明磁场噪声的频率随时间的变化较小。随着实验条件的改变,电场噪声瞬时频率的方差波动较大,说明电场噪声的频率随时间的变化较大。相同实验条件下,电场噪声瞬时频率的方差远大于磁场噪声,说明电场噪声的变化快于磁场噪声。
本文分析了不同测量距离、接触电流下单次电弧辐射电磁噪声的瞬时频率特性,通过使用HHT法获得的IMF分量与实测电磁噪声的功率谱密度对比,发现二者的谱曲线峰值基本重合,表明HHT法可以用于电弧辐射电磁噪声的时频分析。与磁场噪声相比,电场噪声瞬时频率的均方根更大,电场噪声的瞬时频率变化更快。电场噪声与磁场噪声的瞬时频率均方根具有相反的变化趋势。
参考文献(References):
[1] Uchimura K. Electromagnetic interference from discharge phenomena of electric contacts[J]. IEEE Transactions on Electromagnetic Compatibility, 1990, 32(2):81-86.
[2] 郑宁敏, 林匹(Zheng Ningmin, Lin Pi). 刀闸操作电弧对载波通道的干扰(Disturbance caused by the arcing during operation of isolating switch for power carrier channel) [J]. 电力系统保护与控制(Power System Protection and Control), 2009, 37(17):128-131.
[3] 蓝会立, 张认成(Lan Huili, Zhang Rencheng). 开关柜内部弧光短路伴生弧声的频谱特性研究(Study on arc sound signature of fault arc in switchgear) [J]. 高压电器(High Voltage Apparatus), 2008, 44(2):107-109.
[4] Pan Tao, Liang Hong. Analysis of electromagnetic interference from pantograph arcing with different trolley locomotive speed in coal mine tunnel [A]. Proceedings of 8th International Symposium on Antennas, Propagation and EM Theory [C]. 2008. 1096-1098.
[5] 郭凤仪, 王喜利, 王智勇,等(Guo Fengyi, Wang Xili, Wang Zhiyong, et al.). 弓网电弧辐射电场噪声实验研究(Research on radiated electric field noise of pantograph arc) [J]. 电工技术学报(Transactions of China Electrotechnical Society), 2015, 30(14):220-225.
[6] 郭凤仪, 张艳立, 王智勇,等(Guo Fengyi, Zhang Yanli, Wang Zhiyong, et al.). 弓网电弧实验系统与辐射噪声实验研究(Design of pantograph arc experiment system and study on radiated noise) [J]. 电工电能新技术(Advanced Technology of Electrical Engineering and Energy), 2015, 34(12):49-53.
[7] 马云双, 刘志刚, 马岚,等(Ma Yunshuang, Liu Zhigang, Ma Lan, et al.). 动车组弓网离线放电电磁骚扰源模型研究(Electromagnetic disturbance source model of discharge upon pantograph-catenary disconnection of EMU) [J]. 中国铁道科学(China Railway Science), 2013, 34(5):76-81.
[8] 王莉, 曹璐, 严仰光(Wang Li, Cao Lu, Yan Yangguang). 航空交流故障电弧特性的研究(Research on arc fault characteristic for AC power system in aircraft) [J]. 电器与能效管理技术(Electrical & Energy Management Technology), 2011, (2):19-23.
[9] 郭云梅(Guo Yunmei). 航空直流电弧故障检测及保护技术研究(Research on aerospace DC arc fault detection and protection) [D]. 南京:南京航空航天大学(Nanjing: Nanjing University of Aeronautics & Astronautics), 2010.
[10] 熊学军, 郭炳火, 胡筱敏,等(Xiong Xuejun, Guo Binghuo, Hu Xiaomin, et al.). EMD方法和Hilbert谱分析法的应用与探讨(Application and discussion of empirical mode decomposition method and Hilbert spectral analysis method) [J]. 黄渤海海洋(Journal of Oceanography of Huanghai & Bohai Seas), 2002, 20(2):12-21.
[11] 陈平, 李庆民, 赵彤(Chen Ping, Li Qingmin, Zhao Tong). 瞬时频率估计算法研究进展综述(Advances and trends in instantaneous frequency estimation methodology) [J]. 电测与仪表(Electrical Measurement & Instrumentation), 2006, 43(7):1-7.
[12] 任达千, 杨世锡, 吴昭同,等(Ren Daqian, Yang Shixi, Wu Zhaotong, et al.). 信号瞬时频率直接计算法与Hilbert变换及Teager能量法比较(Comparison of instantaneous frequency directed computing method and Hilbert transform and Teager energy method) [J]. 机械工程学报(Journal of Mechanical Engineering), 2013, 49(9):42-48.
[13] 张蓉竹, 蔡邦维, 杨春林,等(Zhang Rongzhu, Cai Bangwei, Yang Chunlin, et al.). 功率谱密度的数值计算方法(Numerical method of the power spectral density) [J]. 强激光与粒子束(High Power Laser and Particle Beams), 2000, 12(6):661-664.
[14] 颜诗源, 李新俊, 吴勋, 等(Yan Shiyuan, Li Xinjun, Wu Xun, et al.). 基于HHT方法的公路运输振动信号时频分析(Frequency analysis for vibration signal from road transportation using HHT) [J]. 装备环境工程(Equipment Environmental Engineering), 2014, (5):23-26.