CN106772036A - 基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法 - Google Patents

Abstract

一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，包括如下步骤：1)机组并网运行时，实时采集发电机转子三相电流i_A、i_B、i_C，利用Morlet小波，分别对其进行变换；2)提取小波脊线；3)根据小波脊线得出瞬时频率信息，检测瞬时频率中是否存在特征频率，即50Hz；4)对转子三相电流信号进行滤波、加窗、插值之后，进行FFT变换，分别比较当前在转子基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异，如果两处幅值都有减小的趋势，判定为发生打火。本发明直接采集转子三相电流，数据获取方便，对双馈风电机组碳刷和滑环间打火的早期诊断具有一定的工程应用价值。

Description

基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测 方法

技术领域

本发明属于风电机组的在线故障诊断技术领域，涉及双馈风力发电机碳刷和滑环火花监测方法。

背景技术

对于MW级的双馈风电机组，励磁电流较大，通常同一相分布着多个碳刷，也使相互之间的均流效果难以得到保证。实际运行中，碳刷和滑环间电阻的变化，极易导致电刷电流的二次分配，使部分碳刷和滑环高温过热，再加上周围空气湿度大，两者之间发生打火的几率会显著增加，轻则损坏碳刷，重则造成发电机滑环及转轴的损坏，影响发电机的正常运行。

实际中导致碳刷与滑环间出现电火花的原因很多，比如：整机的振动过大，引起电刷振动过大；滑环表面有灼伤或表面粗糙使两者接触不良；碳刷运行过热和周围环境湿度大；恒压弹簧压力不够导致碳刷在刷握内摇摆等。若能够对火花进行诊断，及时更换碳刷，将有非常重要的过程意义。

目前双馈风力发电机组碳刷和滑环火花的监测可通过光电管监测法、电刷电流、接触电压的峰值等手段，但是这些方法大都需要增加设备，不仅成本大，而且可靠性不高、适应性较差。近来各研究机构及风机制造商也对两者间的火花监测进行了相关研究，以及公布了一些文章和专利，但大都从电刷电流的角度出发。如“双馈异步发电机电刷滑环电阻变化理论与仿真研究”，该方法通过对电刷与滑环的等效电阻进行分析，测量每个碳刷的电流并进行傅里叶变换，提取特征分量。该方法需要对电刷电流进行处理，在现有风机硬件基础上难以实现。

发明内容

为了克服现有的双馈风力发电机碳刷和滑环火花监测方式的成本大、可靠性较差、适应性较差的不足，本发明提供一种简单易行、准确性良好、适用性较好的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法。

本发明解决其技术问题所采用的技术方案是：

一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，所述监测方法包括如下步骤：

1)机组并网运行时，实时采集发电机转子三相电流i_A、i_B、i_C。利用Morlet小波，分别对其进行变换；

2)提取小波脊线；

3)根据小波脊线得出瞬时频率信息，检测瞬时频率中是否存在特征频率，即50Hz；

4)对转子三相电流信号进行滤波、加窗、插值之后，进行FFT变换，分别比较当前在转子基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异，如果两处幅值都有减小的趋势，判定为发生打火。

进一步，所述步骤1)中，需实时采集发电机转子侧的电流，实际中可采用变流器采集的三相转子电流。

再进一步，在步骤1)中，对转子三相电流进行处理的过程如下：

1.1)设ψ(t)为平方可积函数，即ψ(t)∈L²(R)，如果Fourier变换满足容许性条件：

则ψ(t)为母小波，将其按照尺度参数a和平移参数b进行伸缩和平移，生成的函数族{ψ_a,b(t)}为分析小波。

1.2)变换中所用的Morlet小波时域形式如下。

1.3)根据步骤1.1)、1.2)中，可以得到信号各个频率的相位和幅值信息，对任意的(a,b)，相位可表示为：

对于式(4)，任意的(a,b)，其在t上只有一个驻点，满足：

φ'_a,b(t_s)＝0&φ”_a，b(t_s)≠0(5)。

在步骤3)中，由驻点的性质知φ'_a，b(t_s)＝0，而小波脊线是由相平面上满足t_s(a,b)＝b的点，两者结合即可得到脊线。

因此知道小波脊线a_r(b)，即可求出瞬时频率其中小波脊线的提取可借助小波曲线。

在步骤3)中，小波脊线与瞬时频率的关系如下：

式(2)中，a_r(b)为小波脊线，为瞬时频率。

在步骤4)中，对转子三相电流进行如下处理：

4.1)对转子相电流进行低通滤波、加窗函数，Hanning窗的时域形式如下：

4.2)进行分段线性插值，各插值点用折线连接，再将插值之后的信号做FFT变换；

4.3)将打火和正常运行时的处理结果作比较，观测两者在转子基频和50Hz处频谱的差异。

本发明的有益效果主要表现在：1)本发明涉及的信息的提取无需增加新的设备，硬件成本较低；

2)通过发明所述的方法可以有效的提取特征量，且与正常运行时的数据对比，可以明显观察到打火时的区别。

附图说明

图1为本发明流程图。

图2为滑环和碳刷电阻理论示意图，其中，(a)是等效电路，(b)是碳刷与滑环滑动接触示意图。

图3为实验机组滑环电阻周期性变化图。

图4为打火时转子A相和B相电流。

图5打火相瞬时频率信息。

图6为打火时转子A相电流FFT分析。

图7为正常运行时转子A相电流FFT分析。

具体实施方式

下面结合附图对本发明作进一步描述。

参照图1～图7，一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，包括如下步骤：

1)机组并网运行时，实时采集发电机转子三相电流i_A、i_B、i_C。碳刷和滑环之间发生打火时，会在转子电流中有所体现，现从两者间的电阻理论出发具体介绍，主要考虑接触电阻和滑环的周期性电阻。

碳刷和滑环的接触面存在附件电阻，表达式如下：

式中，ρ_l、ρ_r为碳刷、滑环的电阻率；ε为材料的变形系数；H_l、H_r为碳刷、滑环接触材料的布氏硬度；n、F分别是接触点数和接触压力；ρ_a膜层的电阻率。

滑环的引线触头与碳刷接触部分时刻变化，故滑环电阻也做周期性变化，分析原理如图2所示。从图中可以看到电刷电流经2路流向滑环，一路是L1，另一路是L2，两段电阻可表示如下：

式中，ρ为电阻率，A为滑环径向截面积。

式中ω可由频率f来表示，即ω(t)＝f×t×360。再有R1和R2的并联关系可得总的滑动电阻，其周期性变化如图3所示。

R_tot(t)＝R_ring×f×(t-f×t²)R_ring＝2×π×r×ρ/A(13)

进一步设转子A相电压u(t)＝cos(ωt)，对(4)式按泰勒变换转化成正余弦形式，可求出转子A相电流与碳刷和滑环电阻的关系，两者电阻的变化会引起转子电流在ω₀、ω₀+ω_m频率处的变化。

利用Morlet小波对其变换；

2)提取小波脊线；

3)根据小波脊线得出瞬时频率的变化，检测瞬时频率图中是否存在特征频率，即50Hz。

4)为进一步区分，打火与正常运行时，两者在特征频率及基波频率处的区别，对转子三相电流信号进行滤波、加窗、插值之后，进行FFT变换，分别比较当前在转子基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异，如果两处幅值都有减小的趋势，判定为发生打火。

即当前转子基频和50Hz处频谱处的幅值，与正常工况时转子基频和50Hz处频谱处的幅值相比较，如果同时满足当前转子基频的幅值都比正常工况时基频的幅值小，当前转子50Hz的幅值都比正常工况时50Hz的幅值小，判定为发生打火。

进一步，在步骤1)中，对转子三相电流进行处理的过程如下：

1.1)设ψ(t)为平方可积函数，即ψ(t)∈L²(R)，如果Fourier变换满足容许性条件：

则ψ(t)为母小波，将其按照尺度参数a和平移参数b进行伸缩和平移，生成的函数族{ψ_a,b(t)}为分析小波。

1.2)变换中所用的Morlet小波时域形式如下。

1.3)根据步骤1.1)、1.2)中，可以得到信号各个频率的相位和幅值信息，对任意的(a,b)，相位可表示为：

对于式(9)，任意的(a,b)，其在t上只有一个驻点，满足：

φ'_a，b(t_s)＝0&φ”_a，b(t_s)≠0(5)

1.4)由驻点的性质知φ'_a，b(t_s)＝0，而小波脊线是由相平面上满足t_s(a,b)＝b的点。两者结合即可得到脊线。

因此知道小波脊线a_r(b)，即可求出瞬时频率其中小波脊线的提取可借助小波曲线。

在步骤3)中，小波脊线与瞬时频率的关系如下：

式(2)中，a_r(b)为小波脊线，为瞬时频率。

再进一步，在步骤4)中，对转子三相电流进行如下处理：

4.1)对转子相电流进行低通滤波、加窗函数，Hanning窗的时域形式如下：

4.2)进行分段线性插值，各插值点用折线连接，再将插值之后的信号做FFT变换；

4.3)将打火和正常运行时的处理结果作比较，观测两者在转子基频和50Hz处频谱的差异。

本发明涉及的信息的提取无需增加新的设备，硬件成本较低；

因故障特征量本身很小，提取困难。文中采用瞬时频率的方法进行提取。并且通过分析及实验结果可以明显看到打火时，转子电流基波幅值及特征频率处的幅值较正常运行时偏小。

以某2.0MW双馈风力发电机组为研究对象，转子每相6个碳刷，为使打火易于发生使转子K相2个碳刷运行，实验时转速为1750rpm附近，输出有功1.5MW，转子相电流有效值410A。图2和图3分别是滑环和碳刷间电阻理论的等效示意图和实验机组的滑环周期性电阻，实验时在滑环室观察到的火花现象以及实验结束后，转子K相碳刷表面留下的电蚀现象。

需要说明的是，实际实验时打火一般难以实现，故将转子K相2个碳刷运行，以增加电刷电流密度，增大打火的几率。

记录打火前后转子电流数据，数据包时间长度为10s，图4记录的是该时间段内转子A和B相的电流，其中A相对应K相，B相对应L相，期间A相发生打火。采用Morlet小波分析，其瞬时频率如图5所示，图中平均频率对应转子基频8.3Hz附近，在6.3s左右出现50Hz分量。

为进一步分析，A相打火与正常运行时的区别，分别记录两种工况的数据进行分析。图6和图7分别对应打火时和正常运行时A相电流的FFT分析结果，所用方法如上述所述。打火期间，对应转子基频8.3Hz频谱幅值523.3，特征频率50Hz处幅值8.014；而正常工况时，转子基频处幅值532.6，50Hz处幅值8.394。从结果看出不管是转子基频还是50Hz处，在打火发生时两者幅值明显减小，与根据滑环和碳刷电阻理论分析的结果相一致。

根据实验结果，可以根据50Hz处幅值和转子基频幅值的变化判断碳刷和滑环轻微打火的发生，并且本发明的判定条件是基于大量实验数据的，并不局限于该实例。

本发明通过双馈风电机组转子电流来判断碳刷和滑环之间的打火，简单有效，成本相对较低，是一种有效可靠的判断方法。

最后需要说明的是，以上实施例仅用以说明本发明的技术方案而非限制，本领域技术人员应当理解，可以在形式上和细节上对本发明做出改变，而不偏离发明权利要求书所限定的范围。

还需要说明的是，本文中所公开的实例描述的各参数，仅是为了更好的描述本发明，专业人员可以意识到，通过修改本发明的参数值可达到同样的诊断效果，但是这种实现不应认为超出本发明的范围。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (5)

1.一种基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，其特征在于：所述监测方法包括如下步骤：

1)机组并网运行时，实时采集发电机转子三相电流i_A、i_B、i_C。利用Morlet小波，分别对其进行变换；

2)提取小波脊线；

3)根据小波脊线得出瞬时频率信息，检测瞬时频率中是否存在特征频率，即50Hz；

4)对转子三相电流信号进行滤波、加窗、插值之后，进行FFT变换，分别比较当前在转子基频和50Hz处频谱幅值与正常工况下在转子基频和50Hz处频谱幅值的差异，如果两处幅值都有减小的趋势，判定为发生打火。

2.如权利要求1所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，其特征在于：在步骤1)中，对转子三相电流进行处理的过程如下：

1.1)设ψ(t)为平方可积函数，即ψ(t)∈L²(R)，如果Fourier变换满足容许性条件：

$C_{\mathrm{\&psi;}}=\underset{R}{\&Integral;}\frac{\widetilde{{\mathrm{\&psi;}}^2\left(\mathrm{\&omega;}\right)}}{\left|\mathrm{\&omega;}\right|}d\mathrm{\&omega;}<\mathrm{\&infin;}---\left(1\right)$

则ψ(t)为母小波，将其按照尺度参数a和平移参数b进行伸缩和平移，生成的函数族{ψ_a,b(t)}为分析小波。

${\mathrm{\&psi;}}_{a,b}\left(t\right)=|a{|}^{-0.5}\mathrm{\&psi;}\left(\frac{t-b}{a}\right),a,b\&Element;R,a\&NotEqual;0---\left(2\right)$

1.2)变换中所用的Morlet小波时域形式如下。

$\mathrm{\&psi;}\left(t\right)={\mathrm{\&pi;}}^{-1/4}(e^{-{\mathrm{i\&omega;}}_{0}t}-e^{-{\mathrm{\&omega;}}_0^2/2})e^{-t^2/2}---\left(3\right)$

1.3)根据步骤1.1)、1.2)中，可以得到信号各个频率的相位和幅值信息，对任意的(a,b)，相位可表示为：

对于式(4)，任意的(a,b)，其在t上只有一个驻点，满足：

3.如权利要求2所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，其特征在于：在步骤3)中，由驻点的性质知φ′_a,b(t_s)＝0，而小波脊线是由相平面上满足t_s(a,b)＝b的点，两者结合即可得到脊线。

因此知道小波脊线a_r(b)，即可求出瞬时频率其中小波脊线的提取可借助小波曲线。

4.如权利要求1～3之一所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，其特征在于：在步骤3)中，小波脊线与瞬时频率的关系如下：

式(2)中，a_r(b)为小波脊线，为瞬时频率。

5.如权利要求1～3之一所述的基于双馈风力发电机转子侧电量信息的碳刷和滑环火花监测方法，其特征在于：在步骤4)中，对转子三相电流进行如下处理：

4.1)对转子相电流进行低通滤波、加窗函数，Hanning窗的时域形式如下：

$w\left(t\right)=\left\{\begin{array}{cc}(\frac{1}{2}+\frac{1}{2}cos\frac{\mathrm{\&pi;}t}{T}),& \left|t\right|\&le;T\\ 0,& \left|t\right|>T\end{array}\right.---\left(9\right)$

4.2)进行分段线性插值，各插值点用折线连接，再将插值之后的信号做FFT变换；

4.3)将打火和正常运行时的处理结果作比较，观测两者在转子基频和50Hz处频谱的差异。