设备振动分析仪：从总振值到轴承故障定位的旋转机械诊断方案

简介：《设备振动分析仪交付大庆油田》是一份面向设备管理、状态监测与故障诊断技术人员及石油行业运维人员的案例型文档，围绕VAST-DC23振动分析仪在大庆油田的交付应用展开，可用于了解旋转机械振动监测、故障诊断与油田设备管理场景。压缩包内共1个doc文件，约138KB，内容集中呈现仪器特点、大庆油田概况及现场应用价值，便于快速查阅和归档。已有60人学习，说明该案例对相关从业者具有一定参考意义。读者可从中获取高性价比诊断仪器的功能定位、多参数多频谱监测思路，以及提高设备可靠性、降低维护成本、提升安全性和促进技术创新的应用逻辑；文档以交付案例为线索串联设备选型、现场应用与运维价值，适合作为振动监测方案选型及油田设备管理的参考资料。

1. 从一台「总振值合格却抱死」的注水泵说起

某特大型砂岩油田一个注水站里，值班人员用手持测振表测外输泵驱动端，总振值 4.2 mm/s，落在合格区间；两周后同一台泵的轴承保持架碎裂，非计划停机十一天。事后翻数据才发现，总振值从来不是答案，它把三四个方向的能量压成一个数字，也把早期轴承缺陷的冲击成分抹平在平均值里。这正是设备振动分析仪要解决的事：在旋转机械的振动故障还没长成事故之前，把「哪一根轴、哪一个轴承、哪一类故障」讲清楚。这套方案面向两类人——做设备状态监测课程设计、毕业设计的学生，要给出一套完整可演示的系统；以及要给厂里出一份点检与诊断方案的一线工程师，从测点怎么布、参数怎么选，到频谱怎么读、报告怎么签字，一条线走完。

2. 需求拆解与技术选型：把「专家经验」翻译成四层可验证的能力

「基于专家经验」这句话写在宣传页上很好看，落到设计文档里就必须拆成能测、能算、能判、能签字的四层。否则答辩时被问一句「你的专家经验存在哪儿」，整套方案就悬空了。这一章先把现场需求翻成技术指标，再解释参数和采样为什么这么定。

2.1 现场真正要回答的三个问题

设备管理人员不关心你用的是什么窗函数，他只问三件事：这台机器还能不能继续转、坏在哪个部件、大概还能撑多久。这三个问题对应三类输出，缺一不可。

第一类是状态分级，回答「能不能转」。输出不是一个孤零零的数字，而是按振动烈度四分区分级（A 区良好、B 区可接受、C 区报警、D 区危险）给出的等级，同时标注测点、方向、工况。第二类是故障定位，回答「坏在哪」。这就需要频谱，需要看 1X、2X、谐波族、边带、包络谱里的轴承特征频率，把能量分布翻译成部件名称。第三类是趋势，回答「还能撑多久」。单次测量只能定性，连续多次测量才能画出劣化斜率，给出剩余时间的粗估区间。

很多课程设计只做了第一类，交上去是一张振动值表格，看着完整其实没回答任何问题。我一般的做法是把三类输出做成一张报告的三段：结论区、图谱区、趋势区，后面第 4 章会给出组装逻辑。

2.2 三参数并存：加速度、速度、位移各管哪一段

只用一种参数是这类方案最常见的先天缺陷。旋转机械的故障能量分布在不同频段，单一参数注定只能覆盖一段。

位移（电涡流或低频速度积分）负责 10 Hz 以下的低频段，管的是不平衡、油膜涡动、喘振、机体松动这类「大行程、慢晃悠」的问题。速度负责 10 Hz 到 1000 Hz 的中频段，这一段是人耳和大多数振动烈度评价标准的公共区间，也是判定「整体振得厉害不厉害」的主依据。加速度负责 10 Hz 到 10 kHz 的高频段，滚动轴承的早期点蚀、齿轮的啮合冲击都藏在这里，幅值很小但峭度很高。

三者的关系不是简单积分微分就能互换。加速度积分成速度时，低频噪声会被放大，必须加 10 Hz 高通；速度微分成加速度则会把高频噪声抬高。设计里如果偷懒直接积分，谱图低频端会翘起来，读谱的人很容易把噪声当成 1X 不平衡。通常的做法是：硬件层分通道采集，软件层各自独立分析，只在报告层做交叉印证，不做数据层的互相转换。

2.3 采样率与频率分辨率：决定「看得见」还是「看得清」

采样定理只保证不混叠，不保证分得开。两个参数必须一起定：采样率 fs 和分析点数 N，频率分辨率 Δf = fs / N。

以常见配置为例，加速度通道 fs 取 25.6 kHz，分析点数 8192，Δf 约 3.1 Hz。对一台 2980 r/min 的电机（基频约 49.7 Hz）来说，3.1 Hz 的分辨率足够把 1X 和 2X 分开。但换到低速设备就翻车了：一台 300 r/min 的螺杆泵，基频只有 5 Hz，Δf 3.1 Hz 意味着 1X 和 2X 挤在同一根谱线上，根本没法判断是不平衡还是不对中。解决办法有两个，一是加长采样时间把 N 提到 32768，二是切到低速档降低 fs。

线数（谱线数量）跟 N 是同一件事的两种说法，400 线对应 1024 点、800 线对应 2048 点、1600 线对应 4096 点，工程上习惯按线数报参数。选型时把这张表直接写进设计文档，比写「高分辨率」四个字有用得多。

2.4 仪器形态对比：为什么不是「买个测振笔就行」

测振笔不是没用，它适合高频次的粗筛。问题在于把粗筛结果当诊断结论用，相当于拿体温计判断是肺炎还是骨折。方案设计里我会明确定位：测振笔做日常巡检，分析仪做异常复核和精密诊断，两者是漏斗关系不是替代关系。

这一层想清楚之后，「高性价比」才有落点——不是便宜，是同样的预算覆盖了筛查和诊断两段，不用买两台设备。

3. 总体架构设计：三层结构与四条数据通路

架构图画得再漂亮，落不到测点编码和数据库字段上就是装饰。这一章把系统切成采集层、分析层、诊断层，再把支撑它们的数据通路和存储模型定下来，保证后面写代码时有地方放。

3.1 采集层：测点布置与编码规则

测点布不好，后面算法再精细也是垃圾进垃圾出。滚动轴承座的测点通常布三个方向：水平（H）、垂直（V）、轴向（A）。卧式泵组在驱动端和非驱动端各取两到三个方向；立式设备以水平两个正交方向为主。轴向测点对判断不对中特别敏感，不能省。

测点编码必须有规则，否则几百个测点做完两次巡检就乱了。我一般用「区域-装置-设备类型-设备编号-端位-方向」六段式，段间用短横线连接，全大写定长，方便排序和模糊检索。

# 测点编码生成：六段式，全大写定长，便于排序与检索 def build_point_code(area: str, unit: str, dev_type: str, dev_no: int, end: str, direction: str) -> str: """ area 区域码，如 P01 表示某采油厂一号站 unit 装置码，如 WIP 表示注水 dev_type 设备类型，PMP 泵 / MOT 电机 / CMP 压缩机 / FAN 风机 dev_no 设备序号，三位补零 end 端位，DE 驱动端 / NDE 非驱动端 direction 方向，H 水平 / V 垂直 / A 轴向 """ return f"{area.upper()}-{unit.upper()}-{dev_type.upper()}-{dev_no:03d}-{end.upper()}-{direction.upper()}" # 示例：某站注水装置 7 号泵驱动端水平测点 print(build_point_code("p01", "wip", "pmp", 7, "de", "h")) # 输出 P01-WIP-PMP-007-DE-H

这段看起来简单，价值在约束。定长六段让数据库索引稳定，让报表按字符串排序就等于按区域-装置-设备排序，也让后来接手的人不用查表就知道测点在哪。设备编号补零是为了避免 10 号设备排在 2 号前面这种低级问题，血泪经验，栽过一次就记住了。

3.2 分析层的三条并行通道

同一段时域波形进来，要同时走三条通道，各自出结果，最后在诊断层汇总。

时域通道算有效值、峰值、峰峰值、峭度、峰值因子、裕度因子，输出的是「现在振得有多厉害」。频域通道做 FFT，加汉宁窗，输出幅值谱和相位谱，看的是「能量分布在哪些频率上」。包络通道先带通滤波再 Hilbert 解调，输出包络谱，专门盯轴承和齿轮的早期冲击。

三条通道并行而不是串行，是因为它们服务的判断不同。总振值合格但峭度突然从 3 涨到 6，说明有冲击成分在发育，这时候必须看包络谱。反过来，包络谱干净但 1X 幅值持续上涨，那是典型的不平衡发展过程，跟轴承没关系。只看一条通道就会误判。

3.3 诊断层：规则匹配与置信度打分

诊断层不做「人工智能」，做的是规则匹配加加权打分。每条规则由主判据、辅助判据、排除条件三部分组成，命中主判据得基础分，命中辅助判据加分，命中排除条件减分，最后归一化成 0 到 1 的置信度。

比如不平衡的主判据是 1X 幅值占振动总量 70% 以上，辅助判据是水平与垂直方向相位差接近 90°、轴向振动小，排除条件是 2X 及高次谐波同样显著（那更可能是松动或不对中）。这套打分逻辑后面第 4 章会给完整规则表。

规则库用配置文件存，不硬编码在代码里。现场工程师发现某类设备误报多，改配置就能调，不用重新编译。

3.4 数据模型：三张表撑起全部查询

测点、测次、特征值三张表，加上一张诊断结论表，够用了。不要一上来就设计十几张表，课程设计里表越多越讲不清。

-- 测点档案表：一台设备的每个测点一行 CREATE TABLE t_point ( point_code VARCHAR(32) PRIMARY KEY, -- 六段式测点编码 area_name VARCHAR(64) NOT NULL, dev_name VARCHAR(64) NOT NULL, dev_type VARCHAR(16) NOT NULL, rated_speed INT, -- 铭牌转速 r/min bearing_spec VARCHAR(64), -- 轴承型号，用于算特征频率 created_at TIMESTAMP DEFAULT CURRENT_TIMESTAMP ); -- 测次表：一次测量一行，记录工况 CREATE TABLE t_measure ( measure_id BIGSERIAL PRIMARY KEY, point_code VARCHAR(32) REFERENCES t_point(point_code), measured_at TIMESTAMP NOT NULL, speed_rpm INT, -- 实测转速，优先取键相 load_ratio NUMERIC(5,2), -- 负载率，用于分工况建档 sensor_type VARCHAR(16) -- ACC / VEL / DISP ); -- 特征值表：一次测次对多条通道算出多个特征 CREATE TABLE t_feature ( feature_id BIGSERIAL PRIMARY KEY, measure_id BIGINT REFERENCES t_measure(measure_id), channel VARCHAR(16), -- 1X / 2X / BPFO / RMS / KURTOSIS value NUMERIC(12,4), unit VARCHAR(16) ); -- 诊断结论表：一次测次可给出多条结论 CREATE TABLE t_diagnosis ( diag_id BIGSERIAL PRIMARY KEY, measure_id BIGINT REFERENCES t_measure(measure_id), fault_type VARCHAR(32), -- 不平衡 / 不对中 / 松动 / 轴承损伤 confidence NUMERIC(4,3), -- 0.000 ~ 1.000 suggestion TEXT );

分表的意义在于查询路径清晰：查某台设备的历史趋势，从 t_point 拿编码，连 t_measure 按时间排序，再连 t_feature 取 RMS 通道；查某类故障的历史案例，直接扫 t_diagnosis 按 fault_type 分组。单表堆所有字段看着省事，实际每次查询都要扫全表，测点上到几百个以后现场演示会卡。

4. 核心算法实现：从原始波形到一张能签字的报告

前面定的是骨架，这一章填肉。三个部分：时域指标怎么算得可信、轴承特征频率怎么算得准、规则怎么组合成结论。每段代码后面都要说清楚参数怎么改、失败时看什么。

4.1 时域特征：六个指标各有各的用处

import numpy as np def time_features(x: np.ndarray, fs: float) -> dict: """ x 一维时域波形（已去直流） fs 采样率 Hz 返回六个常用时域指标 """ x = np.asarray(x, dtype=float) peak = np.max(np.abs(x)) # 峰值 rms = np.sqrt(np.mean(x ** 2)) # 有效值 p2p = np.max(x) - np.min(x) # 峰峰值 mean_abs = np.mean(np.abs(x)) # 峭度：对冲击最敏感，正常轴承接近 3 kurt = np.mean((x - np.mean(x)) ** 4) / (np.std(x) ** 4) crest = peak / rms if rms > 0 else 0.0 # 峰值因子 clearance = peak / (mean_abs ** 2) if mean_abs > 0 else 0.0 # 裕度因子 return { "peak": round(peak, 4), "rms": round(rms, 4), "p2p": round(p2p, 4), "kurtosis": round(kurt, 3), "crest": round(crest, 3), "clearance": round(clearance, 3), }

有效值（RMS）对应振动能量，是状态分级的主依据，换算成速度后直接对四分区查表。峰值和峰峰值反映瞬时最大载荷，做强度校核时用。峭度是这套方案里最有价值的一个指标：健康轴承的振动接近高斯分布，峭度在 3 附近；一旦出现点蚀，冲击把尾部拉长，峭度会先于总振值上涨，通常从 3 涨到 5 以上时总振值还没动。这一条就是「早期预警」的技术来源，比任何宣传语都实在。

峰值因子和裕度因子是峭度的补充，对早期微弱冲击更敏感一些，但抗干扰差，一般只在趋势图里做辅助参考，不单独拿来报警。参数上要注意，峭度对单个野点极其敏感，采集时如果混进一次电磁干扰脉冲，峭度能从 3 跳到 15。工程做法是先做 3σ 截断或者中值滤波再算，代码里没写是为了不遮蔽问题，实际部署时这一层要加。

4.2 包络解调与轴承特征频率

轴承的早期损伤产生的冲击很弱，藏在噪声里，直接看频谱看不见。做法是先用结构共振频带做带通滤波，把冲击激起的高频共振放大，再用 Hilbert 变换取包络，最后对包络做 FFT。这样出来的谱，横轴是故障特征频率，不是共振频率。

特征频率靠几何参数算，公式不复杂，但参数输错一个数结果全废。

import numpy as np def bearing_fault_freqs(n: int, d: float, D: float, alpha_deg: float, fr: float) -> dict: """ n 滚动体数量 d 滚动体直径 mm D 节圆直径 mm alpha_deg 接触角，度 fr 轴转频 Hz（注意不是转速 r/min） 返回四个特征频率，单位 Hz """ a = np.radians(alpha_deg) ratio = (d / D) * np.cos(a) bpfo = (n / 2) * fr * (1 - ratio) # 外圈通过频率 bpfi = (n / 2) * fr * (1 + ratio) # 内圈通过频率 bsf = (D / (2 * d)) * fr * (1 - ratio ** 2) # 滚动体自转频率 ftf = (fr / 2) * (1 - ratio) # 保持架频率 return { "BPFO": round(bpfo, 2), "BPFI": round(bpfi, 2), "BSF": round(bsf, 2), "FTF": round(ftf, 2), } # 示例：某 6308 轴承，8 个滚动体，滚动体直径 15.08 mm， # 节圆直径 65 mm，接触角 0 度，轴转频 49.7 Hz print(bearing_fault_freqs(8, 15.08, 65.0, 0.0, 49.7)) # BPFO 约 152.6 Hz，BPFI 约 245.0 Hz

三个参数最容易错：fr 是转频不是转速，49.7 Hz 对应 2980 r/min，直接把 2980 代进去会得到离谱的结果；d 和 D 不能混，滚动体直径和节圆直径差一个量级；接触角对深沟球轴承是 0，对角接触轴承可能是 15° 或 25°，查不到就按 0 算，误差可接受。

判读上有个经验：外圈故障的包络谱峰通常又高又稳，还会出现 2 倍、3 倍谐波；内圈故障因为随轴转动，受载荷区调制，峰两侧会出现转频边带，也就是 BPFI ± fr。看到边带基本可以确认是内圈。保持架故障的 FTF 峰往往很低，容易被噪声淹没，这种要靠峭度趋势而不是单次谱图判断，玄学的地方就在这儿。

带通频带的选取也有讲究。常见的做法是看原始频谱，找 2 kHz 到 8 kHz 之间能量抬升的那一段，把带宽取在共振峰两侧。频率选低了包络谱会很脏，选高了信噪比不够。实在判断不了就多试两三个频带，哪个出来的 BPFO 峰最干净用哪个。

4.3 故障规则库：五类常见故障的判别表

规则写成配置，不写死在代码里：

{ "rule_id": "UNBALANCE_01", "fault_type": "不平衡", "device_scope": ["PMP", "FAN"], "main": {"feature": "1X_RATIO", "op": ">=", "value": 0.70, "score": 0.6}, "aux": [ {"feature": "PHASE_HV_DIFF", "op": "between", "value": [70, 110], "score": 0.2}, {"feature": "AXIAL_RATIO", "op": "<=", "value": 0.30, "score": 0.2} ], "exclude": [ {"feature": "HARMONIC_RATIO", "op": ">=", "value": 0.40, "penalty": 0.4} ], "suggestion": "安排动平衡校正，校正前复测一次确认趋势" }

这套结构的价值是可解释。现场问「凭什么说不平衡」，能直接指出来 1X 占比 0.78、相位差 92°、轴向占比 0.21，三条都是可复现的。用神经网络黑匣子做诊断，在设备科这种要签字担责的场景里推不动，这是选型时就要想清楚的。

阈值不是拍脑袋定的。1X 占比 70% 这条线来自公开的振动诊断工程经验，不同行业会微调，比如风机可以放到 65%，往复机械不适用。设计文档里要写明阈值来源和适用范围，答辩时被追问才有话说。

4.4 报告组装：结论、置信度、建议三段式

一次测量的输出不止一张谱图。我一般组装成三段：第一段是结论，一句话说清状态等级和主要故障；第二段是证据，附上时域指标表、频谱截图、包络谱截图；第三段是建议，给出复测周期或者工单建议。

置信度按规则得分归一化，低于 0.5 的不出结论，只标「待观察」，避免误报把现场折腾疲了。这一条在方案里写清楚，比事后解释说「算法有误差」体面得多。

5. 现场落地避坑：五类最容易翻车的点

实验室里跑通和现场能用是两件事。这一章的五条都是从返工里攒出来的，每条按现象、原因、解决写。

5.1 磁座吸上去，加速度值虚低一半

现象：同一测点，磁座吸附测得加速度 2.1 m/s²，螺栓固定测得 4.6 m/s²，差了一倍多，包络谱上什么也看不到。

原因：磁座本质是一个弹簧-质量系统，它的固有频率会落在 3 kHz 到 7 kHz 之间，正好压在轴承共振带上，把要测的高频成分滤掉了。加上油漆层、锈层、弧面接触不良，耦合刚度更差。

解决：加速度测点优先用螺栓固定或者粘接，磁座只用于粗筛和低频段测量。现场不能打孔的话，用平面磨平一小块漆面再吸，或者选带弧形底面的专用磁座。报告里要标注安装方式，不同安装方式的数据不能放进同一条趋势线。

5.2 转速没测准，整张谱都偏

现象：一台泵的 1X 峰算出来对不上铭牌转速，偏差 2%，导致 2X、3X 定位全偏，包络谱上算出来的 BPFO 也找不到对应峰。

原因：异步电机有滑差，额定 2980 r/min 实际可能 2965 r/min；变频工况下转速随频率变化；铭牌参数跟实际叶轮改了也可能不一致。

解决：有键相就用键相，直接测脉冲间隔算转频；没有键相就在频谱上取 1X 峰的真实位置反推，不要用铭牌数。变频设备必须把实际频率和转速一起记录，后面做阶次分析时用。

5.3 没有基线，只能拿绝对值硬比

现象：新接手的设备，第一次测出来速度有效值 3.8 mm/s，按标准判 B 区可接受，三个月后涨到 6.2 mm/s 才发现问题，中间没有数据支撑趋势判断。

原因：绝对值只能跟标准比，标准是按设备类型和功率分级的粗粒度阈值，反映不出这台设备自己的「健康状态」。

解决：新设备投运一个月内建立基线，同一工况下连续测三次取均值和离散度。之后每次测量跟基线比，涨幅超过 30% 就进观察名单，跟标准阈值形成双判据。基线缺失的情况下，宁可把首次测量标为「基线建立」，不要直接出结论。

5.4 变频工况下阈值误报，一天弹几十条

现象：一台变频泵在 30 Hz 和 50 Hz 两个工况下运行，用同一套阈值，30 Hz 时总振值偏低不报警，50 Hz 时又因为共振区被激起来连续报 C 区，值班员直接把报警关了。

原因：变频设备的激振频率随转速移动，结构的共振点在某些转速下被激发，振动总量跟转速强相关，单一阈值不成立。

解决：按转速分档建档，每档单独设阈值和基线。转速波动超过 5% 的测次不参与趋势外推。共振区识别出来之后，在运行规程里建议避开，或者在结构上加减振处理。

5.5 阈值照搬标准，设备重要度不分

现象：一台备用风机和一台主进料泵用同一条报警线，风机报 C 区没人管，泵报 C 区已经快出事。

原因：振动烈度标准的阈值是针对设备类型和功率的，不区分设备在流程里的关键程度。备机可以带着 C 区等到检修窗口，主泵不行。

解决：在测点档案里加一个关键度字段，分 A（流程关键，停机即减产）、B（重要）、C（一般备用）三级，A 类设备阈值收紧 20%、采样周期缩短一半，C 类放宽。这一条写进方案，评审时最容易被认可，因为它直接对应生产损失。

6. 验收与答辩：用一组对照数据把方案讲透

方案写完不等于交付。我一般按「仪器精度、诊断命中、报告完整」三条线做验收，每条都有可复现的方法和判定值。

这张表的价值在于每一条都能当场演示。幅值线性度用振动台打三个点，频率误差用信号源输标准正弦，轴承识别用带缺陷的试验轴承，重复性在同一点连测三次。评审或答辩现场，与其讲算法多先进，不如把这六行数据摆出来。

答辩时被问得最多的是三个问题，我提前把回答准备好。

第一个是「为什么不直接用总振值加固定阈值」。回答要落到频段：总振值把不同频率的能量压成一个数，轴承早期缺陷的冲击能量小、频率高，混在总量里被低频不平衡成分掩盖，峭度和包络谱才能在总量没动的时候先报警。

第二个是「诊断结论的可靠性怎么量化」。回答置信度的打分机制：主判据、辅助判据、排除条件三部分加权，0 到 1 归一化，低于 0.5 只出「待观察」不出结论。再补一句误报率的统计口径，比空说「准确率高」有说服力。

第三个是「这套东西跟厂里现有的点检制度怎么衔接」。回答漏斗结构：测振笔做日常巡检，筛出异常测点送分析仪做精密诊断，诊断结论进工单系统。不推翻现有制度，只补上诊断这一段。

还有一个容易被忽略的技巧：先把「基线—趋势—报警」这条时间线画出来。评审看方案看多了，对功能列表已经麻木，一条能看清「什么时候测、测了存哪儿、什么时候报警、报警之后谁处理」的流程图，比十页功能介绍管用。

从那以后我每次交付前都强制走一遍验收表，哪怕工期紧也至少把幅值线性度和重复性两项跑掉——这两项不过，后面所有诊断都是建立在错误数据上的空中楼阁。设备诊断这行，测得准永远排在判得准前面。希望这套结构能帮到你，无论是交一份毕设，还是给厂里落一套能签字的点检方案。

本文还有配套的精品资源，点击获取