核电站汽轮发电机组周期性同步振动回顾性振动矢量分析

《Sensors》：Retrospective Vibration Vector Analysis of Periodic Synchronous Vibration in Nuclear Steam Turbine Generator Units

　　

摘要 这项回顾性研究重新分析了来自两台半速核电站汽轮发电机组的59个归档同步（1X）幅值-相位观测数据：25个A型6H测点和34个B型8H测点。研究人员采用Powell优化方法，对重构的实部和虚部联合拟合了一个恒定复响应和一个旋转慢变分量。基线幅值分别为27.96 μm和82.04 μm，慢变分量幅值分别为8.11 μm和6.21 μm，周期分别为340.27 min和365.02 min。全记录幅值R2值分别为0.918和0.883。当每段记录最后约30%的数据被保留不参与拟合时，幅值均方根误差（RMSE）分别为4.19 μm和2.85 μm，R2值分别为?0.782和0.185，相比简单周期幅值拟合仅有适度改善。残差结构和时间窗口敏感性限制了预测能力。归档密封检查记录和冷却趋势为相互竞争物理解释提供了背景信息。无论是分解方法还是有条件的油膜厚度示意，均未识别出唯一热源，也未验证密封到轴弯曲的因果链。研究结果支持对所保留响应模式的描述和比较。

核电站汽轮发电机组周期性同步振动回顾性研究解读

一、研究背景与意义

核电站汽轮发电机组是核电站二回路中的关键能量转换设备，轴系运行稳定性直接关系到机组安全、经济性能和计划性维护。研究对象属于ARABELLE型半速核电汽轮机平台，具有高中压联合模块和多个低压模块，其长跨距、多支承轴系特性促使研究人员关注局部激励、转子响应与支承条件之间的相互作用。研究缓慢变化的同步振动（1X分量）对该类机组的运行维护具有实际意义。

既有研究已分析了汽轮机蒸汽流激励、基础隔振、阀门组失稳、振动探头支架共振、转子-定子碰磨及热致振动、低压缸热变形以及动平衡处理等机制。这些机制即使测量响应中含有强1X分量，也需要不同的诊断证据加以区分。振动评估还取决于测量量和运行状态。ISO 20816-2针对大型汽轮发电机组给出了宽带轴和轴承座振动准则，但其判据不能在没有匹配通道、幅值约定和电厂保护设置的情况下直接应用于分解后的1X分量。本研究关注的监测问题是对稀疏趋势记录中缓慢变化的1X幅值和相位的解释。

对于热致转子振动，已有研究将发电机密封环和轴承油膜的周向轴颈加热与转子热弯曲联系起来。国际研究提供了更直接的反馈证据：Tong和Palazzolo在可倾瓦轴承试验台上用20个传感器测量了周向轴颈温度，识别出温度最大值与平均油膜厚度最小位置之间存在与转速相关的相位分离。后续研究显示枢轴柔性和摩擦影响预测不稳定范围，而挠性枢轴模型将轴承变形与完整转子的热弯曲耦合，并经过实验基准验证。Yang和Palazzolo还证明了瓦块入口温度分布和注油布置对Morton效应预测的重要性。

热弯曲可产生不同的振动轨迹。受控感应加热实验再现了稳定、不稳定和极限环螺旋响应，其周期取决于加热幅值和相位、模态阻尼以及接近临界转速的程度。因此，缓慢变化的1X矢量与热反馈机制相容，但不能识别其热源。迷宫密封研究强调流动引起的交叉耦合力及其对周围流动的敏感性，而主动磁轴承的近期模拟将电磁损耗与转子加热和弯曲耦合。这些机制不同于油密封中的黏性加热，不能在不匹配转子、轴承和密封条件的情况下定量移植到1500 r/min的多支承核电汽轮发电机组。

近期研究还涉及发电机特定油膜。Li等人将剪切流加热、热弯曲和转子响应耦合，并在核电机组上进行了油温调整试验。Zhang等人发现偏心间隙、温度相关黏度和流固换热控制双流氢密封环中的非均匀加热和浮动行为。后者配置不同于本研究的单流密封，但识别了相关耦合机制和参数。

本文回顾性检验了两台ARABELLE型机组记录中的周期性同步振动。线性响应框架提供物理背景，而复矢量分解直接应用于59个保留的1X幅值和相位观测。研究贡献在于可复现的信号分析程序，该程序区分恒定响应基线与慢变分量，检验幅值和相位一致性，并评估残差、时间保持性能和报告测量精度的敏感性。归档密封检查、冷却调整趋势和独立动平衡记录为竞争性解释提供了背景。响应拟合不依赖归档有限元计算，也未识别唯一激励源。

二、关键技术方法

研究人员对两台ARABELLE型机组的历史1X幅值和相位记录进行了回顾性分析。A型机组使用6H通道，包含25个观测点，时间跨度0–693 min，采样间隔10–65 min，中位数26.5 min；B型机组使用8H通道，包含34个观测点，时间跨度0–1008 min，采样间隔15–78 min，中位数28 min。测量系统采用涡流探头测量轴相对位移，标称灵敏度8 mV/μm，幅值以峰峰值位移报告。相位采用键相器参考。

核心方法为复响应矢量分解：将记录的幅值和相位重构为复相量，联合拟合恒定复基线C和恒定幅值的慢变项D exp(iωt)，采用Powell无导数优化方法，在120–1000 min周期范围内以0.5 min网格搜索两种旋转方向，再围绕最佳网格候选进行有界细化。评估指标包括幅值MSE、RMSE、Pearson r、R²、复RMSE和相位RMSE。时间评估保留最后约30%观测不参与拟合，仅用前70%数据进行参数估计。敏感性分析包括首尾80%窗口、奇偶记录、十种初始周期以及1000次循环移动块自举。报告精度敏感性采用±0.5 μm和±0.5°均匀扰动假设，通过2000次蒙特卡洛重拟合传播到拟合参数。

三、研究结果

A型机组数据分析： A型6H测点25个不均匀间隔观测，前17个观测（0–458 min）用于训练，后8个（484–693 min）用于评估。联合拟合重构复相量得到幅值MSE为2.979 μm²，RMSE为1.726 μm，Pearson r为0.961，R²为0.918。复RMSE为3.438 μm，相位RMSE为6.53°。基线幅值27.96 μm，旋转分量幅值8.11 μm，周期340.27 min。最大绝对相对偏差18.15%，出现在低幅值观测点t=390 min。保持数据集的幅值RMSE为4.192 μm，R²为?0.782，表明固定周期模型在该后期区间不能提供可靠外推。

B型机组数据分析： B型8H测点34个原始观测全部保留。相同复响应拟合程序得到幅值MSE为1.473 μm²，RMSE为1.214 μm，Pearson r为0.951，R²为0.883。复RMSE为2.089 μm，相位RMSE为1.18°。基线幅值82.04 μm，旋转分量幅值6.21 μm，拟合周期365.02 min。最大绝对相对偏差3.28%，低于A型，反映较小绝对残差和较高测量幅值，但不证明共同物理机制被独立验证。

两型机组比较： 两种情况下全记录基线幅值均超过旋转分量幅值，但该观察涉及响应矢量而非独立识别的质量不平衡。A型幅值残差范围?3.50至3.08 μm，中位数0.06 μm；B型对应值?2.26至2.75 μm，中位数?0.03 μm。连续记录时间残差间相关性分别为0.768和0.805，表明存在未解析的时间结构。条件95%周期区间A型为325.73–355.74 min，B型为356.85–373.95 min。A型周期在首尾80%窗口间从328.88 min变化到366.11 min，提示不应将单一全记录周期视为固定机器属性。

报告精度影响： 在指定舍入场景下，单观测组合正交贡献A型为0.301–0.344 μm，B型为0.480–0.533 μm，均小于全记录复RMSE值3.438 μm和2.089 μm。因此打印整数精度本身不能解释观测拟合残差。全记录P₂/P₁比值分别为0.290和0.076。时间原点平移和恒定相位偏移后重拟合参数变化小于2×10^?9相对量级，增益0.5和2给出预期比例响应缩放。

辅助现场观测： 2018年1月26日同步快照显示5H=11 μm∠218°、6H=30 μm∠280°、7H=22 μm∠216°、8H=7 μm∠57°，表明同时刻轴响应随轴承位置变化。归档冷却试验显示冷氢温度从约44 °C下降时，两通道幅值后期慢变幅度减小。密封环和座室再加工后监测窗口未见早期长周期波动明显复发，但同期包含负荷降低，非相同运行条件比较。独立动平衡记录报告励磁机侧2.15 kg@80°修正后三通道幅值下降而7V从11 μm增至17 μm，说明平衡评估必须检查所有监测方向。

四、讨论与结论

物理解释方面，在准稳态线性解释下，拟合慢响应系数D_s=H_s(Ω)q，无量纲H_s包含位移约定。对任意正无量纲a，(aH_s, q/a)给出相同D_s，因此无标定传递时P₂=|D_s|识别响应幅值但不能唯一恢复物理弯曲幅值。理想全周期幅值范围A型为19.85–36.07 μm，B型为75.83–88.25 μm，对应观测范围17–37 μm和76–89 μm。理想相位摆动A型33.71°、B型8.69°，记录摆动55°和12°，B型更接近幅值包络，A型更大相位摆动暴露了与圆路径的偏离。相关Arabelle轴系研究发现发电机密封处建模热激励在发电机轴承产生比相邻汽轮机轴承大得多的响应，且未再现整个轴系的观测波动，限制了局部密封加热单独解释每个测点的论断。多种热机制可产生缓慢旋转1X矢量，但现有研究确立物理可能性和判别测量，未识别本机组热源。密封接触痕迹、冷却相关变化和维护后监测窗口与密封相关热效应相容，但负荷和冷却的并发变化阻止因果比较。

工程意义方面，分解的实际用途始于验证振动通道和运行背景，然后检查基线、慢变分量和残差在可比时间窗口内是否保持一致。7V在独立动平衡修正后增加说明需验证所有可用方向上的响应。案例未建立P₁、P₂、其比值或密封油温的通用数值阈值，报警和跳机决策仍由经批准的工厂设置决定。

研究局限性方面，分析限于两段短而不均匀间隔的趋势记录和固定圆路径近似。残差结构、窗口敏感性和A型较弱的保持性能限制了预测。测量范围由保留材料固定。恒定参考不变性和舍入敏感性不建立完整测量不确定度或与过程事件的对齐。原始运行配置（包括热输入站和输出幅值约定）不可用，这些材料仅提供历史背景，不进入59点相量拟合。

研究结论：这项回顾性研究用一致的复响应矢量模型表示两台核汽轮发电机组记录中的缓慢变化1X振动。联合拟合重构实部和虚部分离恒定响应基线与旋转慢变分量。A型和B型全记录拟合中，基线幅值分别为27.96 μm和82.04 μm，旋转分量幅值分别为8.11 μm和6.21 μm，周期分别为340.27 min和365.02 min。幅值R²值分别为0.918和0.883，最大绝对相对偏差分别为18.15%和3.28%。保持后期观测给出幅值RMSE 4.19 μm和2.85 μm，R²值?0.782和0.185。显式整数舍入场景传播给出中心95%周期范围338.99–341.53 min和363.38–366.60 min。该分析提供了所保留振动趋势及其解释限度的透明描述。独立动平衡记录降低三个通道幅值但增加7V，说明需考虑所有监测方向。研究贡献是回顾性响应表征，案例未建立验证的因果诊断或通用运行阈值。