基于Fluent的天然气运...冲蚀模拟分析与防控措施研究_吕超.pdf

试验与研究第46卷 第1期2023年1月Vol.46 No.1Jan.2023HAN GUAN 焊管WELDED PIPE AND TUBE基于Fluent的天然气运输管道弯头冲蚀模拟分析与防控措施研究*吕超1，陈绪鑫1，刘艳龙1，殷宏鑫1，于林辉2（1.东北大学秦皇岛分校 控制工程学院，河北 秦皇岛 066000；2.营口东邦环保科技股份有限公司，辽宁 营口 115000）摘 要：针对天然气运输管道的典型弯道冲蚀问题，依据现场实际生产情况，基于Fluent流体模拟软件建立相应的物理模型，选用standard k-模型、DPM离散相模型和冲蚀速率方程，分析讨论了固体颗粒直径、气体流速以及弯管角度对弯管冲蚀磨损的影响规律。结果显示：流体流速越快，颗粒对弯道的冲蚀越明显，流速与冲蚀率呈正相关；在同一流速下，随着颗粒粒径的增大，冲蚀率逐渐减小，最后趋于稳定；不同角度的弯头受到的侵蚀是不同的，应根据现场操作环境的差异，设计合适角度的弯头。模拟分析结果表明，天然气流速和固体颗粒直径是影响管道弯头冲蚀率的主要因素，并根据模拟结果提出了冲蚀防控措施。关键词：天然气管道；弯头；冲蚀；Fluent仿真模拟；防控措施中图分类号：TE832 文献标识码：A DOI:10.19291/ki.1001-3938.2023.01.003Simulation Analysis and Prevention and Control Measures of Elbow Erosion of Natural Gas Transportation Pipeline based on FluentLYU Chao1,CHEN Xuxin1,LIU Yanlong1,YIN Hongxin1,YU Linhui2(1.College of Control Engineering,Qinhuangdao Branch of Northeastern University,Qinhuangdao 066000,Hebei,China;2.Yingkou Dongbang Environmental Protection Technology Co.,Ltd.,Yingkou 115000,Liaoning,China)Abstract:Aiming at the typical curve erosion problem of natural gas transportation pipeline,according to the actual production situation of the site,the corresponding physical model is established based on Fluent fluid simulation software,and the standard k-model,DPM discrete phase model and erosion rate equation are selected to analyze and discuss the impact of solid particle diameter,gas flow rate and bend angle on bend erosion wear.The results show that the faster the fluid velocity is,the more obvious the erosion of particles on the bend is,and the velocity is positively correlated with the erosion rate.At the same flow rate,with the increase of particle size,the erosion rate gradually decreases and finally tends to be stable.Elbows with different angles are subject to different erosion.Elbows with appropriate angles should be designed according to the differences in the field operating environment.The simulation results show that the flow rate of natural gas and the diameter of solid particles are the main factors affecting the erosion rate of pipeline elbow,and the erosion prevention and control measures are proposed according to the simulation results.Key words:natural gas pipeline;elbow;erosion;fluent simulation;prevention and control measure目前，我国在进行天然气探勘开发及生产利用过程中，从井底直接开采出的天然气往往不是单一的气体，而是天然气、固体颗粒以及一些液体的混合物。其中，固体颗粒和流体介质产生了*基金项目：国家自然科学基金青年基金项目“喷雾热解过程中微纳多孔结构氧化铈颗粒形成机理及模拟研究”（项目编号51904069）。132023年 第 46 卷焊 管良好的耦合性，从而带有一定的速度，具备了动量。在流体流经弯管时，较大的离心力会使颗粒在弯头处发生碰撞，从而对管体形成冲蚀。冲蚀现象严重影响运送管道的安全性和可靠性。由于弯头结构的特殊性，与管道其他部位相比，弯头更易受到颗粒的冲蚀磨损1。因此，研究弯头的冲蚀规律对工程实践具有重要意义。相关研究表明，弯管冲蚀与管道材料和结构、内部流体流速、弯径比及固体颗粒直径都有关系2-9。鲁剑啸10采用正交试验方法分析了 90弯管的冲蚀影响规律，确定了影响弯管冲蚀的诸多因素。现有的研究大多是基于理想输油管道11-12的冲蚀影响规律，对于以天然气为流体介质的研究依然较为缺乏，而且与现实工况结合的研究更少。本研究以黄土高原抽取天然气的现实工况为背景，运用 Fluent 仿真软件，以天然气作为流体介质，设置工况中真实的气体速度与颗粒直径进行数值模拟。对于复杂的地下形貌，通常运用不同角度的弯头13可以有效解决天然气运输问题，因此分析不同角度弯头的冲 蚀 规 律14，对 管 道 设 计 具 有 一 定 的 积 极作用。1理论模型1.1离散相（DPM）模型一般来说，在气固两相流体系中，假设流体相介质是连续的，在具体分析流体相的规律时，将内部颗粒作离散化处理，因此该模型被称为离散相（DPM）模型15。在Fluent模拟中，一般采用欧拉-拉格朗日模型，采用欧拉描述连续相，采用拉格朗日方法描述离散相，其轨迹数学模型为dupdt=FD(u-up)+gx()p-p+FX（1）FD=18pd2pCDRe24（2）Re=dp|up-u|（3）式中：uP颗粒速度，m/s；u流体速度，m/s；p颗粒密度，kg/m3；流体密度，kg/m3；FD方向作用力，N；FX其他作用力，N；g重力加速度，m/s2；Re雷诺数；流体粘度，Pas；dp颗粒直径，m；CD拖曳力系数，CD=1+2Re+3R2e。1.2冲蚀模型根据管道内固液两相流的耦合情况以及管内流体的湍流动能方程，建立冲蚀磨损模型16，模型主要涉及颗粒直径、颗粒质量流量以及颗粒数量。具体的数学模型为Rero=p=1NmpC()dpf()vb()vpAf（4）f()=2.69+1.612-8.843+7.334-1.855（5）式中：Rero管道内侧的冲蚀率，kg()m2s；N碰撞颗粒的数量；mp颗粒质量流量，kg s；C(dp)颗粒粒径函数，取1.810-9；vp颗粒冲击管道壁面的速率，m s；Af管道壁面有限元的计算面积，m2；b(v)相对颗粒速度的函数，取2.6；颗粒与管道壁面的冲击角度，表 1 为冲击角度的具体值及f()。2管道弯头冲蚀模拟2.1管道弯头几何模型建立及有限元处理根据实际输送天然气管道规定表，选择直径 245 mm、壁厚 6 mm 的弯管进行建模。为保证管内流体充分，取进出口长度 L1=L2=8D。将表1冲击角线性分段函数值序号12345/（）020304590f（）00.810.50.4 14第1期吕超等：基于Fluent的天然气运输管道弯头冲蚀模拟分析与防控措施研究HAN GUAN 模型导入 ICEM 中进行网格划分，进出口网格采用“O形”划分，网格单元总数约为30万个，90弯管的几何模型和进出口网格划分如图 1所示。2.2边界条件设置2.2.1物性参数天然气的成分比较复杂，主要是由甲烷、氮气以及烃类等物质组成，具体组分见表2。在Fluent软件中设置多组分运输，通过组分的比例，可以计算出此状态下天然气的密度为0.762 6 kg/m3，黏度为1.12105 Pas。在进行数值模拟时，将天然气看作理想的不可压缩的流体。2.2.2边界条件使用Fluent软件模拟计算天然气对管道内壁冲蚀磨损率。首先，在Fluent中选用standard k湍流模型，湍流强度设置为5%，将模型的入口位置设置为速度进口型，为了保证颗粒受反弹条件的限制，将模型出口处的边界条件设置为自由流出型，并且要满足足够的逃逸条件。为了使仿真模拟结果更加贴近现实，对模型的壁面进行合理的参数设置尤为重要，需要在壁面上设置反弹系数、冲击角函数、粒径函数和速度指数函数，其中粒径函数设置为常数1.810-9，速度指数函数设置为2.6。设置的弯头模拟关键参数见表3。3数值模拟结果分析3.1弯头内流场的数值分析在进行数值计算时，将进口设置为速度入口，出口条件为自由型。其中，设置颗粒直径为200 m，颗粒流量为0.01 kg/s。进行弯头内流场的数值分析是为了预测颗粒冲蚀管道的大体位置，所以只需要模拟一组数据即可。图2为弯管的速度云图和压力云图，从图2可以看出，流体在进入弯头之前的一段距离时，管道内的压力比较均匀，流体在流入弯头时，压力和速度都发生了较大的变化。弯管内的压力变化如图 3 所示，在靠近弯头内侧，压力先减小后增大，速度却是先增大后减小；在图190弯管的几何模型以及“O”形网格划分表2天然气组分表天然气组分CH4CO2N2H2SC2H6H2He比例/%92.44.541.21.20.40.040.02表3管道弯头模型模拟关键参数天然气密度/（kgm-3）0.756 2 气体流速/（ms-1）816颗粒密度/（kgm-3）1 800 颗粒直径/m50250颗粒质量流量/(kgs-1)0.010.1重力加速度/(ms-2）9.81 管道直径/mm245水力直径/mm233 152023年 第 46 卷焊 管靠近弯头外侧，压力先增大后减小，速度是先减小后增大，即速度与压力的分布呈负相关。这是由于流体进入弯头时，流体流向突然发生变化，从而产生较大的离心力，这个离心力对弯管外侧进行挤压、碰撞，压力沿着离心力的方向慢慢增加，从而预测弯管的外侧是受侵蚀最严重的地方，且受到冲蚀的位置一般在弯头的外侧7090之间。3.2流体速度对冲蚀率的影响为了研究颗粒在不同流速下对弯管的冲蚀规律，在模拟过程中保持其他条件不变。管道直径为245 mm，弯径比1.5，颗粒直径200 m，颗粒质量流量分别为 0.01 kg/s、0.05 kg/s、0.1 kg/s，弯管入口流速分别为 8 m s、10 m s、12 m s、14 m s和16 m s，在以上条件下，对弯头冲蚀情况进行数值模拟。图4所示为不同颗粒质量流量以及不同流速下的冲蚀率曲线，图5为8 m/s和12 m/s流速下的冲蚀