甲醇制烯烃(MTO)联合装置是指以煤基或天然气基合成的甲醇为原料,借助类似催化裂化装置的流化床反应形式主要生产烯烃(聚乙烯、聚丙烯等)的工艺技术装备,属于商业性运行的大型煤基甲醇制烯烃工业装置[1-3]。MTO装置中的急冷塔是对甲醇-反应气降温、催化剂等杂质去除以及热量回收的关键设备[4-5],但是经常存在管线振动问题。
目前石化管线常用的减振方法有改变管线结构布置、适量增加支撑刚度、设置孔板、设置缓冲器及安装活塞液压阻尼器等[6-7]。改变管线结构布置、设置孔板及设置缓冲器几类方法均需停车施工,会带来一定经济损失,减振效果也不理想;增加管道支撑刚度法不能耗散振动能量,也容易造成振动的转移[6-7];活塞液压阻尼器只能单自由度减振[8],另外还存在空行程、不能有效解决复杂振动问题。蜂窝型黏滞阻尼器是本科研团队研发的产品,主要由蜂窝活塞、外缸体、传力管夹、高分子阻尼液等组成,可以多自由度减振,不增加管线附加应力,可以不停机进行管线设备减振。
新疆某大型煤化工企业MTO装置中的甲醇-反应气进急冷塔的DN1800管线存在剧烈振动问题,尤其是进塔口附近管线上的大型拉杆横向补偿器,其振动位移峰峰值可达2 650 μm,振动疲劳极易引起补偿器波纹管断裂及管道焊缝的开裂,无法正常开车运行。本文根据急冷塔进塔管线空间结构复杂及拉杆补偿器刚性小等特性,通过计算管线模态和热位移,对阻尼减振方案进行了优化模拟计算,并结合现场空间位置等情况,提出了黏滞阻尼减振方案。该方案不仅可以不停车施工,而且适用大型拉杆横向补偿器刚度小特点,不妨碍其发挥功能作用,有效解决了特大型管线(DN1800)振动问题。
1 急冷塔进塔管线及大型拉杆横向补偿器参数及振动情况MTO装置甲醇制烯烃工艺段,甲醇在反应器内与催化剂接触迅速发生放热催化反应,生成的反应气(乙烯等)与残余甲醇换热至277 ℃后,经DN1800管道送至急冷塔。急冷塔冷却介质为60 ℃水,管线中介质甲醇-反应气的压力为0.2 MPa。进塔管道在进塔口附近存在两个90°弯头,构成Z字形管线,管道温度与急冷塔温度存在180 ℃左右的温度差,管线会产生较大的热膨胀位移。在两个90°弯头之间安装有大型拉杆横向位移补偿器,该补偿器由接管与两个波纹管以及大拉杆等零件构成,可以降低管道热膨胀位移产生的应力[9-11]。现场管线如图 1所示,急冷塔附近的DN1800管线振动剧烈,其中大型拉杆横向位移补偿器振动最大,振动位移峰峰值达2 650 μm。剧烈振动易造成补偿器波纹管疲劳开裂并引起管道弯头焊缝开裂,最终导致介质泄漏和爆炸。
运用ANSYS Workbench 16.0有限元软件对急冷塔的进塔管线进行实体建模和模态分析。由于拉杆补偿器结构复杂,存在拉杆、端板、铰链等结构,计算模态时会出现拉杆、端板等结构的模态,使得计算量大增;但影响管道振动的主要因素与拉杆补偿器辅助构件的模态相关性不大,为了提高计算效率,实体建模时省略了这些辅助构件。计算得到管线前8阶振动频率如表 1所示,第3阶和第7阶模态振型位移云图如图 2所示。
由表 1中的计算结果可见,管线前8阶固有频率分布在7.9~36 Hz范围内,相邻的固有频率之间相差3 Hz左右。现场测量的管道频率为15.5 Hz,非常接近第3阶模态频率(16.8 Hz),可见管线由于介质的冲击力发生了基频振动。当介质流量变化时,激振力频率变化,由于相邻的固有频率比较接近,极易导致管道基频振动。另外计算得到管线前8阶模态位移云图中振动位移最大位置均发生在拉杆补偿器处,图 2只列举了第3阶和第7阶模态位移云图。
2.2 管道热膨胀位移急冷塔冷却水温度为60 ℃,进急冷塔管线(DN1800)中的介质温度为277 ℃,由于管线前后温差较大,会产生较大的热应力和热位移,因此在管线进塔处安装了拉杆横向位移补偿器来补偿管线的横向热位移,降低管线热应力。采用SolidWorks2014软件对进塔管线及拉杆补偿器进行实体建模,再运用ANSYS Workbench 16.0有限元软件对模型进行温度场和热位移变形数值计算,得到管线的温度场分布如图 3所示,管线热位移云图如图 4所示。
由图 3可见,在大型拉杆横向位移补偿器处存在较大的温度梯度,补偿器中间位置的温度为189 ℃左右;由图 4可见,管线弯头处产生的热位移最大,最大位移达39 mm,补偿器处的热位移达28~39 mm。而该拉杆横向位移补偿器设计的补偿管线热位移的范围是0~80 mm,可见补偿器段管线刚度很小,受到流体冲击后更容易产生振动,这也是整条进急冷塔管线的最大振动发生在该处的原因。
3 管线振动原因分析结合第2节有限元数值计算结果和现场测量的相关数据(振动、温度等)可以得到DN1800进急冷塔管线振动的主要原因如下。
① 进急冷塔的管线空间走向复杂,在进塔口附近有两处90°弯头、两处三通结构,介质在弯头以及三通处速度及流向均发生急剧变化,流体对弯头和三通处管道产生脉动冲击[6, 12-14],造成管线振动。
② 由管线模态计算结果和现场测量的管线振动频率(15.5 Hz)可知,管线发生了第3阶基频振动,原因是介质流动激振力频率造成了管线共振。
③ 急冷塔进塔管线振动管段长约20 m,只在3处设置了弹簧支架支撑管道和介质的重量,对管线的约束作用很小,造成管线的刚度不足[6, 12-13],在介质的脉动激振下管道极易发生振动。
④ 在进急冷塔管线第二个弯头前的DN1800管道周向布置了很多DN50的小管(图 5),由于小管向DN1800管道高速喷入冷水对甲醇-反应气进行初步冷却,部分气相介质液化导致管道内压力脉动增大,加剧了管道振动[12-13]。
⑤ 在进塔管线安装的拉杆横向位移补偿器是由两段相同的波纹管及端管、端板、4个拉杆和铰链关节组成的挠性部件,结构见图 6。波纹管的角位移可以吸收管道多平面内的位移,具有补偿量大、刚性小的特点[9-11],而拉杆补偿器的刚度不足,造成管线在该处振动加剧。
设阻尼器提供给振动设备及振动管道的阻尼力为F阻=cVα,其中F阻为阻尼器产生的阻尼力,c为阻尼液的阻尼系数,V为阻尼液蜂窝型活塞运动速度(即管道振动速度),α为阻尼指数[15]。
将振动管道、支撑、黏滞性阻尼器等整个系统进行模型简化分析,得到简化模型如图 7所示。
根据简化物理模型得到管道运动微分方程为
$ m\ddot x + c\dot x + kx = F\sin \left( {\omega t} \right) $ | (1) |
式中,ω为角频率,t为时间,Fsin(ωt)为管道激振力,m为管道质量,c为阻尼系数,k为支架支撑刚度,x为管道位移。求解微分方程得
$ x = B\sin \left( {\omega t - \varphi } \right) = \left( {F/k} \right)\beta \sin \left( {\omega t - \varphi } \right) $ | (2) |
式中,B为受迫振动的振幅,φ为初相位,β为无量纲的振幅放大因子,β越大,振幅越大,且有
$ \beta = \frac{B}{{{B_0}}} = \frac{1}{{\sqrt {{{(1 - {\eta ^2})}^2} + {{\left( {2\zeta \eta } \right)}^2}} }} $ | (3) |
式中,η=ω/ωn为频率比,ζ=c/4(πmf)为阻尼比,B0为静変位移。
绘制振幅放大因子、频率比以及阻尼比之间的关系曲线如图 8所示[15-16]。由图可见,在η=1附近区域,增大阻尼可以明显降低振幅。
黏滞性阻尼器属于速度相关型减振设备[6],可以提供较大的阻尼比,减振效果较好。该方法具有以下优点。
① 响应具有即时性,只要有振动速度,就会提供阻尼力(F阻=cVα)进行减振。
② 可以实现6个自由度减振,没有特定的安装方向;只提供阻尼力,不提供支撑刚度,避免了应力集中和振动转移。
③ 没有机械摩擦副,阻尼液耐高温,使用寿命长,维护简单。
④ 有隔热部件,可对高温管道(350 ℃以下)减振,对于数十微米至几十毫米的振动峰峰值均能有效减振。
黏滞性阻尼器已有数十个石化设备减振的成功案例,减振效果达60%~99%[6, 15]。黏滞性阻尼器主要由缸体、蜂窝型运动活塞、高分子阻尼液、振动能量传递管夹等部件构成[6],其安装形式见图 9。
急冷塔进塔管线(DN1800)的最大振动位置位于拉杆补偿器处,其次在拉杆补偿器前的第二个三通处(图 1(b)),振动峰峰值分别为2 650 μm、1 620 μm。根据振动测量值和现场空间因素,选择在上述位置安装黏滞性阻尼器。
为了预测减振效果和优化安装阻尼器的个数,运用SAP2000 15有限元软件进行阻尼减振模拟计算。建立管道模型,节点数共16 375个,网格边长100 mm,计算材料选择Q235。计算中施加的黏滞性阻尼器刚度13 700 N/mm,阻尼系数137 N·s/mm,阻尼指数0.3,阻尼器作用方向属性设置为3个振动方向和3个扭转方向[17]。在两个弯头及三通处各选择1 300个节点施加载荷,设定载荷时程函数为sine函数;因管道共振频率为16.8 Hz,设定函数周期为0.06 s;每个节点施加力的峰峰值为5 N,各弯头及三通处作用力的峰峰值为6 500 N。
分别计算施加2、3、4、5、6个阻尼器及不加阻尼器时的管道振动位移值,阻尼减振模拟计算结果见图 10。读取测点1(节点14 136)、测点2(节点12 342)、测点3(节点9 154)、测点4(节点4 298)的振动位移值,结果如表 2所示。将施加的阻尼器进行编号并设定施加位置为:阻尼器①(节点11 897)、阻尼器②(节点11 952)、阻尼器③(节点11 001)、阻尼器④(节点11 056)、阻尼器⑤(节点3 035)、阻尼器⑥(节点2 586)。
模拟计算结果整理得到表 2。由表 2可见,施加2个阻尼器时振动降幅为41%~51%,施加3个阻尼器时振动降幅为51%~65%,施加4个阻尼器时振动降幅为65%~80%,施加6个阻尼器时振动降幅为77%~86%。为更直观地表达减振模拟计算结果,将计算结果数据绘制成曲线如图 11所示。根据某设计院管道运行经验,振动频率16.5 Hz时管道振动峰峰值不超过220 μm[17]。在本文模拟减振计算中,施加5个或者6个阻尼器后各测点的振动峰峰值小于220 μm,满足要求。
综合考虑以上管道热位移、管道模态的计算分析,振动原因分析,阻尼减振效果数值模拟计算以及阻尼器安装空间位置和厂方成本要求,最后确定的减振方案为:在拉杆补偿器处安装4个黏滞型阻尼器,在第二个三通处安装2个黏滞性阻尼器。现场阻尼器安装情况如图 12所示。
对管道4个测点的振动峰峰值在安装阻尼器前后进行了测量,改造结果见表 3和图 13。安装黏滞性阻尼器后振动降幅达87.8%~92.8%,最大振动由2 650 μm降为192 μm,保障了管道安全运行。
针对煤基制烯烃MTO装置中急冷塔进塔大型管线(DN1800管段和拉杆补偿器)振动问题,分析了管道热位移、管道模态和管道振动原因,针对拉杆补偿器刚度小、易振动的特性,运用黏滞性阻尼减振技术,在不停车、不改造管道结构、不转移振动、不产生附加应力的情况下,通过施加合适的阻尼,使大型拉杆横向补偿器及特大型管道振动降幅达88%以上,有效控制了急冷塔进塔管道振动,避免了管道因振动过大而发生开裂甚至爆炸等重大事故,保障了生产的安全运行。
[1] |
宗言恭. 甲醇制烯烃工艺(MTO)[J]. 化工技术经济, 1999(4): 18. ZONG Y G. Process of making olefins with methanol(MTO)[J]. Chemical Techno-Economics, 1999(4): 18. (in Chinese) |
[2] |
钱伯章. 陕西延长石油集团MTO联合装置一次投产成功[J]. 石油炼制与化工, 2014, 45(9): 15. QIAN B Z. Shanxi Yanchang oil group MTO joint unit was successfully put into production once[J]. China Petroleum Processing Petrochemical Technology, 2014, 45(9): 15. (in Chinese) |
[3] |
世界首个煤制烯烃项目成功产出MTO级甲醇[J].精细化工原料及中间体, 2010(8): 52. The world's first coal olefin project has successfully produced MTO grade methanol[J]. Fine Chemical Industrial Raw Materials & Intermediates, 2010(8): 52. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-JXHY201008030.htm |
[4] |
叶帅, 袁春亮, 闫程. MTO装置急冷塔和水洗塔塔底搅拌蒸汽的技改分析[J]. 内蒙古石油化工, 2016, 42(5): 36-37. YE S, YUAN C L, YAN C. The technical modification analysis of the steam stirring of the quench tower of the MTO unit and the bottom of the water washing tower[J]. Inner Mongolia Petrochemical Industry, 2016, 42(5): 36-37. (in Chinese) DOI:10.3969/j.issn.1006-7981.2016.05.014 |
[5] |
赵彬. 60万吨/年S-MTO装置分离-急冷塔技术改造[J]. 工程技术(全文版), 2017(3): 124. ZHAO B. Technical transformation of 600 thousand ton/year S-MTO unit separation-quench tower[J]. Engineering Technology (Full Text Edition), 2017(3): 124. (in Chinese) |
[6] |
黄文超, 何立东, 王晨阳, 等. 焦炭塔进料管道振动分析及阻尼减振技术应用[J]. 北京化工大学学报(自然科学版), 2016, 43(6): 83-88. HUANG W C, HE L D, WANG C Y, et al. Vibration analysis of feed pipe for coke tower and application of damping technology[J]. Journal of Beijing University of Chemical Technology(Natural Science), 2016, 43(6): 83-88. (in Chinese) |
[7] |
SOROKIN S V, OLHOFF N, ERSHOVA O A. Analysis of the energy transmission in spatial piping systems with heavy internal fluid loading[J]. Journal of Sound and Vibration, 2008, 310(4): 1141-1166. |
[8] |
蔡文军, 王平, 祝远征, 等. 一种液压阻尼器的结构及阻尼性能分析[J]. 机床与液压, 2006(6): 149-150, 153. CAI W J, WANG P, ZHU Y Z, et al. Structure and damping performance analysis of a hydraulic damper[J]. Machine Tool & Hydraulics, 2006(6): 149-150, 153. (in Chinese) DOI:10.3969/j.issn.1001-3881.2006.06.054 |
[9] |
张丽霞, 冯朝晖. 大拉杆波纹补偿器在供热管网中的应用[J]. 设备管理与维修, 2005(5): 27-28. ZHANG L X, FENG C H. Application of large pull rod ripple compensator in heating pipe network[J]. Plant Maintenance Engineering, 2005(5): 27-28. (in Chinese) DOI:10.3969/j.issn.1001-0599.2005.05.019 |
[10] |
袁立方. 基于CAESARⅡ的大拉杆横向型波纹补偿器在供热管道设计中的应用分析[J]. 轻工科技, 2013, 29(1): 36-37, 60. YUAN L F. Application and analysis of large pull rod transverse corrugated compensator based on CAESAR Ⅱ in the design of heating pipeline[J]. Light Industry Science and Technology, 2013, 29(1): 36-37, 60. (in Chinese) |
[11] |
方明宇.直管压力平衡复式拉杆型补偿器[C]//第十三届全国膨胀节学术会议论文集——膨胀节技术进展.石家庄, 2014: 4. FANG M Y. Direct tube pressure balance compound tension rod type compensator[C]//Proceedings of the 13th National Expansion Joint Academic Conference-Progress of Expansion Joint Technology. Shijiazhuang, 2014: 4. (in Chinese) http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGVE201410002025.htm |
[12] |
姜杨, 何立东, 伍伟. 丙烷塔空冷器集合管管道阻尼减振技术研究[J]. 石油化工设备技术, 2011, 32(2): 19-21, 25. JIANG Y, HE L D, WU W. Study on vibration attenuation of damping technology in pipe net of air cooler of propane tower[J]. Petro-Chemical Equipment Technology, 2011, 32(2): 19-21, 25. (in Chinese) DOI:10.3969/j.issn.1006-8805.2011.02.006 |
[13] |
刘明, 何立东, 裴正武. 丁二烯换热器出口管道振动原因分析及阻尼减振技术研究[J]. 石油化工设备技术, 2012, 33(5): 20-23, 69-70. LIU M, HE L D, PEI Z W. Analysis of vibration causes of outlet pipe of butadiene heat exchanger and research on damping technology[J]. Petro-Chemical Equipment Technology, 2012, 33(5): 20-23, 69-70. (in Chinese) DOI:10.3969/j.issn.1006-8805.2012.05.007 |
[14] |
胡朋, 何立东, 张震坤, 等. 基于阻尼减振技术的热电厂减温减压器管道研究[J]. 机电工程, 2014, 31(1): 43-46. HU P, HE L D, ZHANG Z K, et al. Study on the pipe of reducing temperature and pressure reducer in thermal power plant based on vibration attenuation of damping technology[J]. Journal of Mechanical & Electrical Engineering, 2014, 31(1): 43-46. (in Chinese) |
[15] |
吕江.管道及旋转机械减振技术研究[D].北京: 北京化工大学, 2015. LV J. Study on vibration attenuation technology of pipe and rotating machinery[D]. Beijing: Beijing University of Chemical Technology, 2015. (in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10010-1015726072.htm |
[16] |
倪振华. 振动力学[M]. 西安: 西安交通大学出版社, 1988: 79-84. NI Z H. Vibration mechanics[M]. Xi'an: Xi'an Jiao Tong University Press, 1988: 79-84. (in Chinese) |
[17] |
胡朋.管道阻尼减振技术研究及其应用[D].北京: 北京化工大学, 2014. HU P. Research and application of damping technology for pipe attenuation[D]. Beijing: Beijing University of Chemical Technology, 2014. (in Chinese) |