湿汽在孔板和文丘里管内流动压降特性分析.docVIP

湿汽在孔板和文丘里管内流动压降特性分析 　　摘 要：该文采用CFD两相流模拟软件对蒸汽干度较高的湿汽在孔板和文丘里内流动过程进行了数值计算，分析了压力、流量、干度等参数对孔板和文丘里管内湿汽压降特性的影响，并比较了两种相同喉径比孔板和文丘里管对于湿汽测量的虚高系数的影响，为优化压差式湿汽测量装置结构奠定了基础。 　　关键词：孔板 文丘里管 湿汽 CFD数值模拟 　　中图分类号：TK262 文献标识码：A 文章编号：1674-098X（2015）08（a）-0056-03 　　向井下注入高温高压蒸汽是石油开采工业中稠油的主要开采方法[1]。油田注汽锅炉生产的蒸汽，经过长距离的注汽管线到达油田井口后，通常是含一定量液体的湿蒸汽，而稠油开采中起主要作用的是湿蒸汽中能量品质较高的干蒸汽。准确测量井口注入的干蒸汽流量，根据现场油井的生产状况，获得每口井每次注入的蒸汽量，确定整个油区蒸汽的总需求量，建立在线湿蒸汽流量控制对于减少采油过程中蒸汽的损失，提高油田节能减排效率具有重要的意义。 　　现有的油田开采过程中，井口的蒸汽普遍采用采油孔板或文丘里等压差式流量计测量，由于汽液两相湿蒸汽的两相流动特性导致测量结果高于实际干蒸汽流量，产生虚高。湿气测量虚高模型的研究是两相流基础研究重要内容之一。国内外许多学者研究得到了各种不同的虚高模型，比如Murdock模型[2]、Bizon模型[3]、林宗虎模型[4]等。上述虚高模型受被测介质、工况等因素影响，虚高模型的可移植性较差。 　　因此，为了获得精度较高的湿汽测量结果，需要对一定工况条件下测量虚高的影响因素进行研究，得到较为准确的虚高模型。该文采用CFD软件对于汽液两相湿蒸汽在孔板和文丘里管内流动过程数值模拟，对比分析了压力、流量、干度等参数对喉径比相同的孔板和文丘里管内湿汽压降的影响特性，分析了两种结构对于湿汽测量的虚高系数的影响，为优化压差式湿汽测量装置结构奠定了基础。 　　1 孔板流量计内是汽液两相蒸汽流动过程CFD数值分析 　　现有的CFD软件比如FLUENT、CFX等集成了大量的两相流模型，能够对复杂几何结构的各种流型条件进行数值仿真模拟，能够较清晰地分析流动过程。CFD方法已经广泛应用于航空航天、机械制造等领域的设计中。 　　汽相为连续相，液相为弥散相且具有稳定流动结构的汽液两相混合物在孔板内湿汽流动过程是十分复杂。精确分析孔板内湿汽流动过程，获得影响湿汽虚高系数的影响因素，是建立湿汽测量虚高模型的基础。 　　汽液两相蒸汽总流量G=5-10t/h，干度x=0.5-0.8，压力p=10-16.5MPa。由此可以计算出该工况参数条件下汽相的体积分数在0.80～0.96，液滴的体积分数在0.05～0.20之间，汽液两相蒸汽呈现出雾状流的流动形态。 　　1.1 离散相模型（DPM） 　　DPM模型[5]采用欧拉-拉格朗日方法。体积份数较大的主流体相被处理为连续相，直接求解时N-S方程，而离散相是通过计算流场中大量的粒子（气泡或是液滴）运动得到的。DPM模型可以计算离子浓度小于20%的两相流动。根据上文的计算，该文采用离散相（DMP）模型对孔板内湿汽的流动过程进行数值模拟，分析压力、流量、干度等参数对虚高系数的影响。 　　1.1.1 控制方程 　　（1）动量方程。 　　（1） 　　式中，αq为体积份数；u为流速（m/s）；ρ为密度（kg/m3）；p为压力（Pa）；SMA为体积力（Pa/m）；Ma为相间力（Pa/m）；下标：q为相，上标：T为u的逆矩阵。 　　（2）连续性方程。 　　（2） 　　式中，SM为质量源项（kg/m3s）。 　　1.1.2 几何模型及边界条件 　　如图1（a）所示文丘里管（WQL）几何结构，文丘里内径D=60 mm，喉部前端的渐缩角21°，后端的渐扩角12°；如图1（b）所示孔板（KB）几何结构，孔板内径D=60 mm。文丘里管和孔板喉径比两种结构β1=0.4，β2=0.7。 　　计算网格采用四边形结构网格，网格总数4.3万。进口设定为速度进口边界条件，出口设定为压力出口边界条件，中心采用轴对称边界条件。两相流模型采用DPM模型，湍流模型采用k-ε模型。计算工况压力p=7-10 MPa，流量G=5-10 t/h，干度x=0.5-0.8。 　　1.2 湿汽在文丘里管和孔板内流动特性分析 　　图2～图3的工况参数范围压力p=8 MPa，湿汽总流量G=10 t/h，干度x=0.5文丘里管和孔板内的速度分布计算结果，如图2所示湿汽在文丘里管内速度分布，可以发现文丘里管入口段的速度场较均匀，在喉部段速度增大，在渐扩段速度减小；喉部前后的速度变化较大，而且一直影响到渐扩段后部的稳定段。如图3所示，湿汽在孔板孔口前端速度降低区域较小