1. 课程导入与船舶设备维护场景认知

船舶是一个高度集成的移动工业系统，主机、舵机、泵组、风机等设备在盐雾、振动、变负荷的恶劣环境中长期连续运行。一旦关键设备在航行途中失效，轻则延误船期，重则危及航行安全。因此，船舶设备维护从来不是"坏了再修"这么简单，而是围绕"如何在故障造成后果之前介入"展开的一整套工程决策。本课程要做的，就是用 Python 把这种决策从经验判断推进到数据驱动。

从三种维护策略说起

船舶维护策略大致经历三个阶段。事后维护（ corrective maintenance）在设备失效后更换或修复，成本最低但风险最高，只适合非关键、易更换的辅件。定期维护（preventive maintenance）按运行小时或日历周期强制检修，例如主机每 8000 小时检查缸套、泵组每季度更换密封，它降低了突发故障率，但存在"过度维护"——状态良好的部件被提前拆换，以及"维护不足"——周期未到却已劣化的矛盾。

预测性维护（predictive maintenance, PdM）则依据设备实际状态决定何时维护：通过传感器采集振动、温度、压力、电流等信号，提取特征并建立健康指标，在劣化趋势达到阈值前安排检修。它的核心价值不是"更频繁地监测"，而是把维护时机从固定周期转移到设备真实状态上，从而同时降低非计划停机与不必要的备件消耗。对船舶而言，这意味着可以结合航次计划，把需要停机的检修安排在靠港期间，而不是在海上被动抢修。

四类典型设备的故障模式

不同设备的失效机理决定了监测信号与特征的选择，这是后续建模的物理基础。

• 离心泵：常见故障为轴承磨损、叶轮不平衡与汽蚀。轴承外圈剥落会在振动频谱上产生特征频率及其谐波；汽蚀则表现为高频宽带能量上升并伴随流量、出口压力波动。监测重点通常是泵体轴承座振动、进出口压力和电机电流。
• 风机：以叶片积灰不平衡、轴承故障和皮带打滑为主。积灰导致 1 倍转频振动幅值缓慢上升，皮带打滑则表现为转速比异常与传动侧温升。监测重点为轴承振动、转速与电机电流。
• 主机：柴油主机故障模式复杂，包括缸套磨损、活塞环断裂、喷油器堵塞、排气阀烧蚀等。工程上常用缸压、排温、滑油压力与金属磨粒浓度作为健康指标，其中各缸排温偏差是判断燃烧异常最直接的信号之一。
• 舵机：液压舵机的典型问题是液压油泄漏、油泵磨损与阀件卡滞，表现为转舵响应变慢、油压建立不足和油温升高。监测重点为液压油压力、油温、转舵角度与响应时间。

可以看到，四类设备的共同点是：故障都有一个从轻微到严重的渐进过程，且这个过程会在可测量的信号上留下痕迹。预测性维护的任务，就是把这个"痕迹"变成可比较的数字。

上图把本节的核心逻辑串成一条信号链：设备故障模式决定该测什么信号，信号经特征提取变成健康指标，再由 Python 完成趋势建模与阈值判断，最终输出维护决策。底部色条强调一个关键事实——预测性维护真正有价值的窗口，是"早期劣化"到"明显劣化"之间那段可检出但尚未失效的区间。如果等到功能失效才报警，那只是事后维护的数字化版本。

Python 在其中的角色

船舶设备的数据往往来自不同厂商的采集系统，格式不统一、采样率不一致、存在缺失与漂移。Python 的价值在于用一套轻量、可复现的代码把这条链路打通。具体来说，它承担四件事：用 pandas 完成多源数据的对齐与清洗；用 numpy 与 scipy 做重采样、滤波和 FFT 频谱分析；用 matplotlib 绘制趋势曲线与频谱图，让劣化过程"看得见"；用 scikit-learn 或简单的统计方法建立基线模型与异常判据。相比专用诊断软件，Python 方案的优势是透明、可修改、便于把领域知识（例如某型泵的轴承特征频率）直接写进代码。

下面是一段最小可运行的示例，演示如何从一段振动信号中提取 RMS 与主频幅值，作为健康指标的候选特征。数据为教学构造的模拟信号，不代表任何实测结果：

import numpy as np
from scipy.fft import rfft, rfftfreq

fs = 2000                      # 采样频率 Hz（示例）
t = np.arange(0, 2, 1/fs)      # 2 秒数据
f_rot = 25                     # 转频 25 Hz，对应 1500 r/min
# 模拟信号：转频分量 + 2 倍频 + 噪声
sig = (0.8*np.sin(2*np.pi*f_rot*t)
       + 0.3*np.sin(2*np.pi*2*f_rot*t)
       + 0.05*np.random.randn(len(t)))

rms = np.sqrt(np.mean(sig**2))          # 时域有效值
spec = np.abs(rfft(sig)) / len(sig)     # 幅值谱
freqs = rfftfreq(len(sig), 1/fs)
peak_amp = spec[np.argmin(np.abs(freqs - f_rot))]

print(f"RMS={rms:.3f}, 转频幅值={peak_amp:.3f}")

运行后应得到约 RMS≈0.60、转频幅值约 0.40 的结果（因随机噪声每次略有差异）。这两个数字本身没有意义，意义在于把它们按时间排列成序列：当轴承逐渐磨损时，转频及其谐波幅值会持续上升，RMS 随之抬升。工程上通常先采集设备健康期的数据建立基线，再设定如"连续 3 个采样周期超过基线均值加 3 倍标准差"的判据，避免单点噪声造成误报。

课程学习路径与实战项目

本课程按"认知—数据—特征—模型—决策"的顺序推进。前段建立船舶设备故障模式与信号对应关系的认知；中段集中训练 Python 数据处理能力，包括时间序列清洗、重采样、滤波与频谱分析；后段进入建模，从简单的阈值与趋势外推，逐步过渡到异常检测与剩余寿命估计；最后落到维护决策，把模型输出转化为可执行的检修建议。

贯穿全程的实战项目是"船用离心泵健康监测与预警"。项目使用一组按真实工况特征构造的泵组运行数据（含振动、出口压力、电机电流与轴承温度），要求完成：数据清洗与统一时间轴、特征提取、健康指标构建、劣化趋势识别，并在指标越过预警线时输出包含设备编号、异常特征与建议检修窗口的预警记录。完成这个项目后，你将拥有一套可迁移到风机、舵机等设备的分析框架——换的是信号与特征，不变的是从物理机理到数据判据的这条主线。