基金项目:国家重点研发计划“中低压纯氢与掺氢燃气管道输送及其应用关键技术”,2021YFB4001604;国家自然科学基金“超高压压缩机气缸-管路气流脉动双向传播耦合规律与解耦抑制机理研究”,52076166;压缩机技术国家重点实验室(压缩机技术安徽省实验室)开放基金项目“基于数字图像处理的往复式压缩机智能故障诊断方法研究”,No. SKL-YSJ202111。
(收稿日期:2023-06-19;修回日期:2023-07-06;编辑:张雪琴)
1.西安交通大学能源与动力工程学院;2.西安交通大学化学工程与技术学院;3.合肥通用机械研究院有限公司·压缩机技术国家重点实验室(压缩机技术安徽省实验室)
1.School of Energy and Power Engineering, Xi'an Jiaotong University; 2.School of Chemical Engineering and Technology, Xi'an Jiaotong University; 3.Hefei General Machinery Research Institute//State Key Laboratory of Compressor Technology (Anhui Province Key Laboratory of Compressor Technology)
hydrogen compressor, gas pulsation, modal analysis, vibration control
DOI: 10.6047/j.issn.1000-8241.2023.08.012
氢气作为一种公认的清洁能源和重要的化工原料,在各个行业均有广泛应用。氢能在应用时,一般需要在较高的压力下进行,因此氢气压缩机对于氢能的使用具有不可替代的作用。往复式压缩机在各种压缩机中密封性能较好、热效率较高,且能提供最高的压缩比,因此在石化行业被大量使用,其正常、高效运行至关重要[1]。管道系统是氢气压缩机成套装备中重要的组成部分,管道的振动直接影响压缩机装备的正常、安全运行,尤其是对于氢气这种易燃易爆的工艺介质,严重的管道振动会带来极大的事故风险[2]。据统计, 20世纪末全世界范围内石化工业事故中31%的事故是由管道系统故障造成的,带来了巨大的经济损失[3]。因此,对氢气压缩机进行管道振动分析,减轻其振动具有重要意义。
国外对于往复式压缩机管路气流脉动以及管路振动的研究较早,在20世纪50年代,美国开展了管路气流脉动的相关研究,提出了平面波动理论。Maclaren等[4]综合考虑了管道系统与外界热交换的影响,建立了更加完善的特征方程组,适用于计算一维非定常流动,更加接近管道内实际的流动情况。Wiggert等[5]在不考虑重力、压力及阻尼等外载荷影响的情况下,将等直长管道视为梁单元,得到了管道振动方程。文献[6]给出了组合管道的气柱固有频率的计算方法,但该计算方法仅对若干简单管道组合有效,与实际情况相差较远,并不能用于解决工程实际问题。中国对于往复式压缩机的管道振动研究较早的是西安交通大学,其振动科研小组于1976年总结了两种计算复杂管道系统内气流脉动的方法,根据计算结果,决定在某氮氢压缩机和某空压机上通过使用孔板来实现管道的减振[7]。徐稼轩等[8]构建了专用程序,应用有限单元法,为压缩机管系组成部分建立了正确的力学模型进行振动计算,并应用该程序设计了某化肥厂的压缩机管系,效果良好。鲁录义等[9]提出了一种管道气流脉动-机械振动一体化模型,通过ANSYS多物理场实现了压缩机管道气流脉动与机械振动的高度耦合。徐斌等[10-11]研究指出,以大容器为分段点,将复杂管路分成简单管路进行气流脉动计算,其准确性可以满足工业要求。韩文龙等[12]研究了两台往复式压缩机并联运行时管道系统中关键位置的气流脉动,应用CFD方法分析了湍流和层流模型下管路系统的气流脉动,对比得出湍流模型计算管路气流脉动的准确性较之层流模型更高的结论。翟延科等[13]运用ANSYS软件对出口处管道剧烈振动的某往复式压缩机管道模型进行模态分析,得到其前5阶固有频率和振型,并提出了相应的管路振动治理方案。江志农等[14]提出应用LMS Test.lab与Bentley AutoPIPE相结合对往复式压缩机管道系统进行振动测试的新方法,可以提取压缩机管道系统的固有频率、振型等模态测试结果。袁伟[15]对往复式压缩机管道振动特性、减振技术进行了全方面的研究,搭建了往复式压缩机管道振动实验系统,对一段实验管路采用ANSYS软件进行数值模拟研究,并对某石化炼油企业压缩机系统的振动治理案例进行了详细说明。目前,对于管路振动的研究已较为全面,但大部分对于压缩机振动的研究缺少系统的流程方法,且大多通过商业软件进行仿真分析,缺少现场的实际测试。在此,针对变工况氢气压缩机一级进气管路振动剧烈的现象,结合现场测试,提出了一套完整的管路振动分析流程,并制定了有效的治理方案。
往复式压缩机管路振动的主要原因集中于两个方面,即管道内气流脉动和管道机械共振。气流脉动是由往复式压缩机间歇性吸排气造成的,是往复式压缩机的一个固有特性。脉动气流在经过阀门、弯头、异径管、三通等管道元件时,其流速和压力变化会产生激振力,从而引起管道振动,尤其是在发生气柱共振时,管道内产生的压力波动最为剧烈,引起的管道振动也最为剧烈[16]。孔板、缓冲罐、亥姆霍兹共鸣器等元件均可用来衰减管道内的气流脉动[17];同时,当管道的某阶固有频率与压缩机的某阶激发频率相等或接近时,也会使管道发生机械共振[18]。按照一般规律,管道的固有频率越低,发生机械共振的可能性越大,根据API-618-2007《石油化工和天然气工业用往复式压缩机》规定,往复式压缩机管路系统元件的最小机械固有频率应大于压缩机最大额定转速的2.4倍。对于发生机械共振的管路,可以通过改变管路的走向、添加支撑或改变支撑位置的方式来改变管路的固有频率,从而达到防止机械共振的效果[19]。此外,管道振动还可能由外力及压缩机主机振动引起[20]。
氢气压缩机属于往复式压缩机,其吸排气具有间歇性,相比于其他形式的压缩机会有较大的管道内气流脉动,可因此引起管路振动;同时,氢气与同温度、压力下其他常见工质相比声速高、密度极小,导致氢气压缩机管路内气柱固有频率、气流脉动与其他工质的往复式压缩机不同,造成其管道设计、布置与其他工质往复式压缩机不同。管路振动易导致管道破裂,由于氢气易燃易爆的特性,一旦发生泄漏将会造成严重事故。因此,需要对氢气压缩机管路的振动原因进行分析。
以某石化工厂的氢气压缩机为例,其安装在工作环境较好的压缩机厂房,且经过现场监测,压缩机主机振动较小,在可接受范围内,因此,排除由外力及压缩机主机振动引起的管路振动。在此,通过现场测试及仿真研究,从气流脉动和管道机械共振两个方面进行氢气压缩机管路振动原因分析,探究氢气压缩机一级进气管路振动剧烈的原因。
该氢气压缩机为两级两列对动式压缩机(表1),气量可从40%调节到100%,通过部分行程压开进气阀实现无级气量调节。该压缩机一级进气管路振动剧烈,通过手持式测振仪对一级进气缓冲罐与过滤器之间管路的振动进行测量(图1),设置了A、B两个监测点:将一级进气管路(即缓冲罐与过滤器之间管道上某点)的测试位置设为监测点A,再将缓冲罐与气缸之间接管上的测试位置设为监测点B。由此,可得到监测点A处在各个气量负荷工况下x、y方向上的振动值(表2,定义一级进气管路中重力方向为z方向,缓冲罐轴向为y方向,垂直yz平面方向为x方向)。由于x 、 y方向振动值较大,z方向振动值很小,只记录了各个气量负荷工况下x、y方向上的振动值。可见,气量从40%调节到100%各工况下,y方向振动值均较大,综合两个方向上的振动值,在60%气量负荷工况下一级进气管路振动最为强烈,其振动值严重超过ISO 20816-8-2018《机械振动 机械振动的测量和评估 第8部分:往复式压缩机系统》规定的限定值(12.7 mm/s)。
表1 某氢气压缩机运行参数及结构参数表Table 1 Operating and structural parameters of a hydrogen compressor
图1 某氢气压缩机一级进气管路气体流程及振动测试位置示意图Fig. 1 Gas flow process and vibration test position of primary gas intake pipeline of a hydrogen compressor
表2 某氢气压缩机在不同气量负荷工况下监测点A处振动值数据表Table 2 Vibration data of a hydrogen compressor at measuring point A under various gas loads
对该氢气压缩机进行现场气流脉动测试,从多方面分析判断压缩机异常振动的原因。压缩机转速为300 r/min,现场检测得到该氢气压缩机的监测点A处在气量负荷为60%时振动最为强烈。针对此现象,对压缩机在各个气量负荷工况下的运行状况分别进行测试:先将气量负荷从40%阶段性每隔10%升至100%,并进行测试;再将气量负荷从100%阶段性每隔10%降至40%,进行测试;最后再将气量负荷从40%阶段性每隔10%升至100%进行测试,从而确保数据的真实性及准确性。
在监测点A、监测点B处分别安装高频动态压力传感器以在每种气量负荷工况下进行脉动测试,取稳定运行后的压力数据进行分析,经过时域、频域分析得到各个频率分量下压力数据(表3、表4)。可见,压缩机在60%气量负荷工况下的气流脉动确实较大,其中测点A处的脉动最大。
为了进一步将压力数据可视化,取60%气量负荷工况下测得的监测点A处压力脉动时域结果(图2,取两个周期内的结果),经过计算其压力不均匀度为2.88%。
为更好地分析其脉动情况,对时域数据进行FFT(Fast Fourier Transform)分析,得到其各阶频率分量下的脉动结果,压缩机转速为300 r/min,对应基频为5 Hz,其整数倍分别对应各阶脉动分量下的频率。API-618-2007规定了缓冲装置管侧接管处或接管以外的最大许用脉动范围,即API-618-2007限定值,其计算式为:
式中:p1为各脉动分量对应于基本谐波频率的最大许用峰峰值,以管路平均绝对压力的百分比表示;a为气体的声速,m/s;pL为管路的平均绝对压力,MPa;DI为管道内径,mm;f为脉动频率,Hz。
表3 气量负荷从40%升至100%时一级进气缓冲罐接管和一级进气管路的时域、频域信息表Table 3 Time and frequency domain data of primary gas intake buffer tank nozzle and primary gas intake pipeline with the gas load rising from 40% to 100%
表4 气量负荷从100%降至40%时一级进气缓冲罐接管和一级进气管路的时域、频域信息表Table 4 Time and frequency domain data of primary gas intake buffer tank nozzle and primary gas intake pipeline with the gas load dropping from 100% to 40%
图2 某氢气压缩机监测点A处在60%气量负荷工况下的压力脉动时域结果图Fig. 2 Time domain results of pressure pulsation at measuring point A of a hydrogen compressor under 60% gas load
将FFT分析得到的各频率分量下脉动幅值与API-618-2007限定值进行对比,即可得到60%气量负荷工况下监测点A处的压力脉动频域结果图(图3)。
图3 某氢气压缩机监测点A处在60%气量负荷工况下的压力脉动频域结果图Fig. 3 Frequency domain results of pressure pulsation at measuring point A of a hydrogen compressor under 60% gas load
从图3可知,该氢气压缩机压力脉动2阶分量的脉动幅值最大,是典型的双作用压缩机管路脉动情况。但经过计算并与API-618-2007限定值对比,管道内的脉动分量幅值远小于限定值,均满足API-618-2007标准要求,再结合表3、表4,初步分析气流脉动对压缩机一级进气管路异常振动的影响较小。因此,判断管路振动较大的原因并非是气流脉动。
脉动模拟基于平面波动理论进行,可以用来分析稳定流动时管道内的气流脉动情况,可对任何工况下往复式压缩机管道系统的压力脉动及激振力进行预测。为相互验证测试与仿真分析的准确性,对100%气量负荷工况下一级进气管路进行气流脉动仿真分析,并与测试结果进行对比。结合管路设计参数与现场测绘,按照实际管道尺寸与走向建立一级进气管路气流脉动分析模型(图4):管道内径为154.06 mm、壁厚为7.11 mm,缓冲罐筒体内径为800 mm、壁厚为16 mm,缓冲罐长度为2 565 mm,图中绿线处表示对应管道在y方向、z方向长度均为995 mm。考虑压缩机阀腔对气流脉动的影响[21],在管道与分离器以及压缩机连接处分别设置开口边界条件、压缩机边界条件,气流从前分离器中进入管道,流经缓冲罐后进入压缩机气缸内进行压缩。设置好边界条件以及输入压缩机结构参数后,通过分析计算可得到监测点A处在100%气量负荷工况下的气流脉动仿真计算时域及频域结果,并与测试得到的100%气量负荷工况下的频域结果进行对比(图5)。可见,各阶脉动分量的结果基本相同,在100%气量负荷工况下相互验证了测试与仿真分析的准确性,也验证了一级进气管路振动较大的原因并非是管道内的气流脉动。
图4 某氢气压缩机一级进气管路气流脉动分析模型图(mm) Fig. 4 Gas pulsation analysis model of the primary gas intake pipeline of a hydrogen compressor (mm)
图5 某氢气压缩机监测点A处在100%气量负荷工况下气流脉动频域仿真结果与测试结果对比图Fig. 5 Comparison between simulation and test results of gas pulsation frequency domain at measuring point A of a hydrogen compressor under 100% gas load
管道机械共振也是造成管道振动的重要原因,因此需要对该压缩机一级进气管路进行模态测试,分析其振动特性。由于整个管路较长,在整条管道上布置测点较为困难,而振动较大的位置靠近一级进气缓冲罐(过滤器与一级进气缓冲罐之间),且管道在过滤器处设置有固定支架,因此在过滤器及一级进气缓冲罐出口之间布置12个测点(图6)进行模态测试,得到一级进气管路前两阶振型(图7,图中白色线为管道初始位置、彩色线为振型,颜色与3个坐标轴颜色相对应,与某一方向坐标轴颜色相同代表向该方向偏离,偏离初始位置越多表示相对振幅越大)。可见,一级进气管路系统的1阶、2阶的固有频率分别为13.4 Hz、24.6 Hz,其中1阶主要沿x方向振动,2阶则主要沿y方向振动。该压缩机的转速为300 r/min,对应的激发频率为5 Hz及其倍频。一般以管道固有频率的10%计算共振的频率区域,测试得到管路的2阶固有频率为24.6 Hz,该压缩机激发频率的3倍频恰好处于管道2阶固有频率的共振区域内。
图6 某氢气压缩机模态测试以及动态振动测试测点位置示意图Fig. 6 Measuring points for modal test and dynamic vibration test of a hydrogen compressor
图7 某氢气压缩机一级进气管路前两阶振型测试图Fig. 7 First two orders of vibration modes of the primary gas intake pipeline of a hydrogen compressor
经过上述分析,考虑到管路振动可能是在管道2阶固有频率处发生共振,因此需要对振动较大的位置(即图6中的测点5)进行动态振动测试,通过在压缩机管路上放置加速度传感器,采集其各个气量负荷工况下稳定运行时的振动加速度数据进行分析可知,在60%气量负荷工况下,振动分量达到最大值。因此,在测点5处选取x、y、z共3个方向测量其振动加速度,并记录60%气量负荷工况下加速度的时域结果(图8a),并对其进行频域分析(图8b)。可见,3个方向下测点5在25 Hz均有较大的振动分量,且y方向的振动最大,结合模态测试结果可以判断,一级进气管路由于其自身结构的固有频率在25 Hz处产生了共振现象。
图8 某氢气压缩机在60%气量负荷工况下稳定运行时测点5的振动加速度的时域和频域图Fig. 8 Time and frequency domains of vibration acceleration during stable running at measuring point 5 of a hydrogen compressor under 60% gas load
利用数值模拟进行模态分析是开展振动分析时常用的一种有效分析方法。ANSYS Workbench结合了结构、流体、磁场等多域多物理场及其耦合仿真,能够满足多个行业的需求,在此,采用ANSYS Workbench中Modal模块对一级进气管路系统建立有限元模型进行模态分析。选择与模态测试相同的管路,网格划分采用六面体网格,过滤器出口、气缸入口以及缓冲罐支撑位置采用固定约束,管夹1~管夹3位置(图9)处采用弹性约束。根据管道、缓冲罐的实际走向和结构尺寸(图4)进行适当简化(忽略法兰中的螺栓螺纹),建立三维仿真模型(图9)。根据模态测试得到的固有频率及振型调整管夹的约束大小,最后得到与实际情况基本相同的结果,前两阶固有频率为13.587 Hz、24.539 Hz。从模拟得到的振型图(图 10)可知,1阶振型主要是x方向振动位移较大, 2阶振型主要是y方向振动位移较大,符合实际模态测试结果。因此,可以在此基础上进行振动治理方案的模拟。
图9 某氢气压缩机一级进气管路三维模型图(mm) Fig. 9 3D model of primary gas intake pipeline of a hydrogen compressor (mm)
图 10 某氢气压缩机一级进气管路前两阶振型模拟云图Fig. 10 Simulation results of the first two orders of vibration modes of primary gas intake pipeline of a hydrogen compressor
由于共振导致管路振动,可以通过改变管道固有频率来减小振动。对于管道的约束越多,其各阶固有频率越大,因此,可以通过增加管夹来改变管道的固有频率。管路振动大的位置处于靠近缓冲罐的两个弯头处,考虑到缓冲罐轴向管道部分已有一个管夹,且倾斜的管道部分不便加装管夹,故考虑将新增的管夹安装在竖直管道靠近上方弯头部位(图 11)。由于现场缓冲罐位于二楼,故将管夹的支架设置在二楼的框架上,支架高度为1 700 mm,仅在图中标注出新增管夹及支架,原有管夹及支架不再标注。支架的结构采用H型钢的形式,新增管夹主要限制y方向的振动,因此支架竖直部分H型钢的腹板与y轴平行。经过数值模拟可知,添加支架后一级进气管路前两阶固有频率分别变为16.10 Hz、32.383 Hz,刚好可以避开压缩机激振频率的倍频,从而避免共振。现场按照该治理方案实施后,测试得到一级进气管路的振动值完全降至可接受的范围内。
图 11 某氢气压缩机一级进气管路新增支架位置示意图(mm) Fig. 11 Position of new supports of primary gas intake pipeline of a hydrogen compressor (mm)
针对某一级进气管路振动剧烈的氢气压缩机进行了气流脉动及管路振动分析,提出了一套完整的管路振动原因及治理的分析方法:首先,对其一级进气管路进行气流脉动分析,通过现场测试以及仿真分析进行相互验证,证明了管道内气流脉动符合API-618-2007标准限定值,即管路振动并非是由气流脉动造成的;然后,进行振动特性分析,对其一级进气管路进行模态测试以及动态振动测试,通过数据分析可以发现压缩机激发频率的3倍频刚好处在管路的2阶固有频率共振区域内,因此判断管路振动过大的原因是发生了机械共振;最后,根据模态测试的结果建立三维仿真模型,并在此基础上提出振动治理方案,将管道的固有频率避开压缩机的激振频率,从而避免管道机械共振。现场实施后减振效果显著,验证了该方法的可行性。
对于非共振原因的管路振动,仍需进一步对管道进行强迫振动响应分析,找到其振动原因。此外,可以优化气流脉动分析模型,如将气缸工作过程以及气阀运动影响考虑在内,从而得到更精确的气流脉动计算结果。
主管单位:中国石油天然气集团公司
主办单位:中国石油天然气股份有限公司
管道分公司
编辑出版:《油气储运》杂志社
通信地址:河北省廊坊市金光道51号
邮政编码:065000
主 编:艾慕阳
副 主 编:关中原
电 话:0316-2177193
0316-2176753
传 真:0316-2177392
发行范围:公开发行
国内发行:河北省廊坊市邮政局
国外发行:中国出版对外贸易总公司
(北京782信箱)100011
国际刊号:ISSN 1000-8241
国内刊号:CN 13-1093/TE
邮发代号:18-89
广告许可:1310004000002