4、谐波检测与测量技术:傅里叶变换(FFT)原理、加窗FFT与频谱泄漏、谐波测量仪器(电能质量分析仪)使用、谐波测量点选择

各位工程师朋友,大家好。今天我们来聊聊谐波检测与测量技术。说实话,这块内容我当年刚入行时也踩过不少坑。你想想看,谐波看不见摸不着,要是测不准,后面所有的治理方案都是白搭。所以这一节,我把自己十几年在现场摸爬滚打的经验,结合理论,给大家掰开揉碎了讲清楚。

4.1 傅里叶变换(FFT)原理——从时域到频域的“翻译官”

谐波分析的核心工具,就是傅里叶变换。说白了,它能把一个复杂的时域波形,拆解成不同频率的正弦波分量。我们平时用示波器看到的电压电流波形是“时域”的,而谐波分析需要的是“频域”信息——哪个频率有多少幅值、多少相位。

FFT(快速傅里叶变换)是DFT(离散傅里叶变换)的快速算法。我在项目中遇到过一位刚毕业的同事,拿着MATLAB直接对采样数据做FFT,结果出来的频谱图乱七八糟。为什么?因为他没搞懂几个关键参数。

FFT的核心参数:

  • 采样频率 fs:必须大于最高谐波频率的2倍(奈奎斯特定理)。我建议至少取4~5倍,留足余量。
  • 采样点数 N:决定了频率分辨率 Δf = fs / N。比如 fs=10kHz,N=1024,那么Δf≈9.77Hz。
  • 频率分辨率:你想分辨出49.9Hz和50.1Hz吗?那Δf必须小于0.2Hz,对应的N就要很大。

举个例子,一个50Hz基波叠加了5次谐波(250Hz)的信号:

// 伪代码示例:FFT分析
采样率 fs = 6400 Hz  // 128倍基波,够用
采样点数 N = 1024
频率分辨率 Δf = 6400/1024 = 6.25 Hz

// 对采样数据做FFT
X = fft(signal, N)
// 取前N/2个点,对应频率 0 ~ fs/2
f = (0:N/2-1) * fs / N
// 幅值计算
amp = abs(X(1:N/2)) * 2 / N

嗯,这里要注意:FFT出来的结果是对称的,我们只取前一半。幅值要乘以2/N才是真实幅值(直流分量除外)。

4.2 加窗FFT与频谱泄漏——一个让我头疼很久的问题

实际工程中,我们不可能采集无限长的信号。截断信号时,如果截取的不是整周期,就会发生频谱泄漏。什么意思呢?就是本来应该集中在50Hz一根线上的能量,泄漏到了旁边的频率上,看起来像“裙边”一样。

为什么会这样? 因为截断相当于在时域乘了一个矩形窗,矩形窗的频域是sinc函数,旁瓣很高。我刚开始做风电并网测试时,用矩形窗分析谐波,结果5次谐波旁边冒出来一堆莫名其妙的频率分量,查了半天才发现是泄漏造成的。

我的经验: 风电并网场景下,基波频率会波动(49.5~50.5Hz),很难做到整周期采样。所以加窗是必须的。我个人习惯用汉宁窗(Hanning),旁瓣衰减快,主瓣宽度适中。如果对幅值精度要求高,可以用布莱克曼窗(Blackman)。

加窗FFT的步骤很简单:

  1. 选择一个窗函数 w(n),长度等于采样点数 N
  2. 将采样数据 x(n) 与窗函数相乘:y(n) = x(n) * w(n)
  3. 对 y(n) 做FFT
  4. 注意:加窗后幅值会衰减,需要做幅值修正(乘以修正系数)
// 加窗FFT示例(Python风格)
import numpy as np

fs = 6400
N = 1024
t = np.arange(N) / fs
signal = np.sin(2*np.pi*50*t) + 0.3*np.sin(2*np.pi*250*t)

# 汉宁窗
window = np.hanning(N)
# 加窗
windowed_signal = signal * window
# FFT
X = np.fft.fft(windowed_signal, N)
# 幅值修正(汉宁窗的修正系数约为2/1.5)
amp = np.abs(X[:N//2]) * 2 / N / (np.sum(window)/N)

避坑指南: 我曾经犯过一个低级错误——加窗后忘了做幅值修正,结果谐波幅值偏小了一大截。另外,如果信号中有间谐波(非整数次谐波),加窗也无法完全消除泄漏,这时候需要用更高级的方法,比如插值FFT或锁相环同步采样。

4.3 谐波测量仪器(电能质量分析仪)使用——工欲善其事,必先利其器

现场用的谐波测量仪器,主流的有Fluke 43B/435、Dranetz、日置PQ系列等。我最早用的是Fluke 43B,那家伙皮实耐用,就是屏幕小了点。现在新出的仪器功能强大多了,但核心操作逻辑差不多。

使用要点:

  • 接线方式:三相四线制用Y型接法,三相三线制用Δ型接法。我见过有人把电压钳位接反了,测出来的谐波相位全错。
  • 量程设置:先预估电压电流大小,别让信号饱和或太小。风电并网出口电压一般是690V或10kV,电流几百到上千安,注意配CT/PT变比。
  • 采样设置:一般仪器默认10周波(200ms)做一次FFT,符合IEC 61000-4-30标准。但如果你要分析快速变化的谐波,可以设成1周波。
  • 数据记录:别只记总谐波畸变率(THD),要记录各次谐波的幅值和相位,特别是5、7、11、13次。我习惯同时录波形,方便事后回放分析。

现场小技巧: 测量前先做一次“背景噪声”测试——断开被测设备,看看电网本身的谐波水平。有一次我在风电场测谐波,发现5次谐波特别大,结果一查是电网背景谐波,跟风机没关系。差点冤枉了风机厂家。

4.4 谐波测量点选择——位置不对,数据白费

测量点选在哪里,直接决定了数据的有效性。我总结了几条原则:

测量点位置 适用场景 注意事项
风电机组出口(690V侧) 评估单台机组谐波发射特性 注意避开机组自身滤波器的影响
集电线路并网点(10kV/35kV侧) 评估多台机组叠加后的谐波水平 考虑线路阻抗和电容效应
公共连接点(PCC) 考核并网电能质量是否达标 按国标GB/T 14549要求测量
滤波装置前后 评估滤波器效果 同时测滤波前后数据,对比分析

我个人习惯的做法: 先在PCC点做长期监测(至少一周),摸清谐波的时间分布规律。然后针对谐波大的时段,到具体机组出口做短时测量,定位谐波源。最后在滤波装置前后测,验证治理效果。三步走,稳得很。

注意: 测量点要选在CT/PT的二次侧,注意二次回路不能开路(CT)或短路(PT)。我曾经遇到一个现场,PT二次侧保险丝烧了,仪器显示电压为零,折腾了半天才发现。所以测量前一定要检查二次回路是否正常。

4.5 本章知识体系总览

下面这张图,是我梳理的谐波检测与测量技术的核心逻辑。从信号采集到FFT分析,再到仪器使用和测量点选择,每一步都有坑,每一步也都有技巧。

谐波检测与测量技术知识体系 信号采集 采样率、采样点数、量程 FFT分析 频率分辨率、幅值计算 加窗处理 汉宁窗、幅值修正 电能质量分析仪使用 接线、量程、采样设置、数据记录 测量点选择 机组出口 → 集电线路 → PCC → 滤波装置前后 核心原则:从源头到并网点,逐级排查,定位谐波源

好了,这一节的内容就到这里。谐波检测是谐波治理的第一步,也是最重要的一步。测不准,后面全是白干。下一节我们会聊到谐波标准与限值,看看国标到底是怎么规定的。


公众号:蓝海资料掘金营,微信deep3321