4. 风机振动信号预处理

从风机振动采集系统拿到的原始数据，几乎不可能直接送进频谱分析或预警模型。传感器零点漂移、温度引起的低频漂移、工频干扰、电磁脉冲、丢包、采样率不统一、启停阶段混入的转速变化，都会让频谱图上出现虚假的峰或把真实故障特征淹没。预处理的目标不是"把数据变好看"，而是让后续的时域指标、包络谱和阶次分析建立在物理上可信的信号之上。本节按工程清洗的实际顺序，把八个环节串成一条可复现的处理链，并给出可参考的参数与判断依据。

一、去趋势与去均值：先固定"零线"

去趋势（detrend）和去均值（remove mean）解决的是同一个问题的两个层次。去均值只减去整段数据的算术平均值，消除直流偏置；去趋势则进一步拟合并扣除缓慢变化的基线，通常用最小二乘拟合一条直线或低阶多项式。对于风机振动信号，加速度传感器常因安装应力或电荷放大器漂移产生线性趋势，速度信号经积分后还会叠加二次趋势，因此去趋势比单纯去均值更彻底。

工程上建议的顺序是：先去趋势，再去均值。以一段 10 s、采样率 25.6 kHz 的轴承座加速度信号为例，若原始均值约为 0.08 g，说明存在明显直流分量；用一次多项式去趋势后，均值通常降到 1e-4 g 量级。判断是否处理到位，可以看处理前后时域波形的中心线是否回到零附近，以及 0 Hz 附近谱峰是否被压下去。需要注意：去均值会移除信号的直流分量，如果后续要做积分求速度或位移，必须在积分前完成去趋势，否则趋势项会被积分放大成严重的低频漂移。

二、滤波：按频带保留有效成分

滤波决定"哪些频率进入分析"。风机振动分析常用的滤波有三类：低通用于抗混叠和去除高频噪声；高通用于去掉趋势和低频晃动；带通用于锁定关注频带，例如滚动轴承故障常集中在 1 kHz～10 kHz 的高频共振带。滤波器类型上，FIR 线性相位适合需要保持波形形状的场合，IIR（如 Butterworth）阶数低、计算省，但会引入相位失真。

参数选择要结合风机转速和特征频率。假设某风机转频 f_r = 25 Hz（1500 r/min），滚动轴承外圈故障特征频率 BPFO ≈ 3.5×f_r = 87.5 Hz，若只关心轴承早期故障，可设置带通 500 Hz～8000 Hz 做包络分析；若关心不平衡和不对中，则应保留 1×、2×、3× 转频，带通下限不能高于 25 Hz。截止频率不能贴着关注频率设置，一般留 10%～20% 过渡带。滤波后应检查通带内幅值是否基本不变、阻带衰减是否达到预期（如 40 dB 以上），并确认没有因阶数过高产生振铃。

三、重采样：统一时间基准

风机常采用变转速运行，若按固定时间间隔采样，转频会随转速漂移，频谱上 1× 峰变宽甚至分裂。重采样分两种：等时间重采样用于统一不同传感器的采样率；等角度重采样（阶次重采样）则依据转速脉冲，把时间域信号转换为角度域信号，使转频恒定为 1 阶，故障特征以固定阶次出现。工程上常用转速键相（tacho）信号做参考，通过插值把非均匀角度采样点映射到均匀角度网格。

重采样前必须先做抗混叠滤波，否则降采样会把高频噪声折叠进低频。例如把 25.6 kHz 降到 2.56 kHz，需先低通到 1.28 kHz 以下。重采样后要核对：转速稳定段内 1× 幅值是否集中、阶次谱中 BPFO 是否落在理论阶次附近。若没有键相信号，也可用振动信号自身的转速谐波做参考，但精度会下降。

四、异常值剔除与缺失值处理

异常值多来自电磁脉冲、传感器松动或采集卡丢帧，表现为孤立的尖峰或整段饱和。常用判断方法是 3σ 准则或中位数绝对偏差（MAD）：|x_i − median| > k·MAD 时判为异常，k 取 3～5。对已确认的异常点，不能简单删除，否则会破坏时间轴连续性；推荐用邻近点线性插值或样条插值替换，并记录替换位置和数量。若某段数据异常比例超过 5%，应标记该段为低质量数据，不参与预警建模。

缺失值处理与异常值类似：短时丢包（如连续少于 10 个点）用线性插值补齐；长时丢包应分段处理，把缺失段前后各自作为独立数据段，避免插值引入虚假低频。处理完成后应统计缺失率和异常率，作为数据质量报告的一部分。

五、数据分段与加窗

频谱分析要求信号近似平稳，因此需要把长信号切成短段。段长由频率分辨率决定：Δf = f_s/N，若希望分辨 1 Hz 且 f_s = 25.6 kHz，则 N = 25600 点，对应 1 s。工程上常取 0.5～2 s 一段，并设置 50%～75% 重叠，以减少信息损失。

直接截断会产生频谱泄漏，因此每段要加窗。常用窗函数及适用场景：Hanning 窗主瓣窄、旁瓣衰减快，适合一般频谱分析；Hamming 窗旁瓣更低但主瓣略宽；Flat-top 窗幅值精度高，适合校准和幅值测量；矩形窗仅适合整周期采样或瞬态冲击。加窗后需做幅值修正，Hanning 窗的幅值恢复系数约为 2.0（单边谱），否则测得的振动幅值会偏小。分段加窗后，可对多段做平均，得到更稳定的自功率谱，再进入包络解调或阶次分析。

六、工程落地建议

把上述环节串成一条可复用的处理链时，建议固定顺序：去趋势 → 去均值 → 抗混叠滤波 → 重采样 → 异常值/缺失值处理 → 分段加窗。顺序不能随意调换：先重采样后滤波会导致混叠；先分段后去趋势会让每段趋势不一致；先加窗后去均值会破坏窗函数对称性。每一步都应保留中间结果和参数记录，便于回溯。

对于在线预警系统，预处理还要考虑实时性。FIR 滤波阶数不宜过高，重采样尽量用查表插值，分段长度与模型输入长度对齐。离线分析则可以放宽计算量，采用更精细的样条插值和多次平均。无论在线离线，最终都要输出一份数据质量摘要，至少包含：原始点数、有效点数、缺失率、异常率、处理参数、有效数据段数。只有预处理链稳定、参数可追溯，后续的频谱、包络和阶次分析结果才具备可比性和可信度。