天然气节流过程数值模拟研究.pdfVIP

中困丁程热物理尝-／-厶z。 传热传质学 学术会议论文 编号：063064 天然气节流过程数值模拟 许 中 徐明海 (中固石油人学(华东)热能_T程研究所山东东营257061) Tel：0546—8393759 Email：xuzhong_upc@163．com 摘要：天然气管道运行过程中，经常在节流管件处因节流降温导致生成水合物，造成生产事故。本文 采用流动‘j传热数值汁算技术模拟研究了天然气在节流过程中爪力、温度及速度的变化关系，．’i天然 气水合物形成的压／J温度关系结合，可以判断柏：节流孔处足杏町能形成水合物，为防止水合物堵采措 ’ 施的制定提供技术依据。 关键字：节流；数值模拟；水合物 1．引言 天然气在管道内流经节流管什时，会出现节流降压和降温现象，若降温幅度过火， 会导致天然气水合物生成，导致管路堵塞。目前通常采用图版法或简化模型求解节流温 降，其主要理论依据是平衡态下的焓相等，主要研究平衡状态下的气体温度变化【l圳。 ：1了流过程中的气体压力、温度和速度随位置的变化较少有研究。 本文利用流动与传热的数值计算方法研究了管道内天然气流经节流件的速度、压力 以及温度的变化情况，，结合水合物形成的压力温度关系，判断出节流孔附近是否可能形 成水合物。 2．物理模型与控制方程 本文研究天然气流经一标准节流孔板的流场和温度场，孔板尺寸符合国家标准 GB／T2624--1993的规定，物理模型如图1所示。 图中划剖面线的部分为标准节流孔板。管道直径D=0．1m，，孔板直径d=0．04m，节 体流动是否充分发展。这将通过后面的计算加以验证。 节流过程是一不可逆的绝热过程，在此过程中假设： 1)不考虑体积力；2)不考虑热辐射的影响：3)流动过程中无化学反应；4)流动 为定常流动。5)流体与外界无热交换。 ’ nlet, 荔 Outlet 占 j惑￡。6r” 荔。一 图i含节流乱板的管道，J：意图 坐标系设置如卜．：原点住节流孔入口的中心，沿管道流动方向为X轴正向，YOZ 878 平面为管道的横截面。管道内大然气通过1，流装置的流动是三维、稳态可乐缩流动问题， 连续性方程、动餐方程和能最方程的统一形式为：． 这里0为，“义源项，名、S≯分别代表由丁压力及速度引起的源项。通过给源项和扩导 赋值，可以分别得到动餐方程和能量方程。 3．网格划分 本文利用分块划分网格，对。1了流孔板选较小的网格尺寸，而对节流孔扳前、后区域 均选较人的网格尺寸；然后将此二块区域迕接为一个整体，生成的整个计算区域网格， 结果如图2所示。 图2计算区域删恪 4．求解模型和边界条件的设定 根据前面描述的物理模型，选择三维稳态分离隐式解算器，压力一速度的耦台方式 选择SIMPLEC算法，密度、动量、紊流动能、紊流耗散率以及能量方程均选择二阶迎 风格式。湍流模型选择标准的k-e模型。入口和出口都给定压强，管壁绝热肇。由于天 然气主要成分为甲烷，故计算了纯甲烷气体在通过节流孔板时，流场和温度场的变化情 况。假定甲烷气体为理想气体，给定入口压力为3MPa，温度为310K，流速为10m／s： 给定节流后出口压力为1．4MPa。 ． 5．计算结果及分析 由丁几何结构和流动关于管道中心轴对称，在中心轴上的流动最能反映出流