基于VOF模型的水平管道油水两相流动数值模拟.docVIP

中国工程热物理学会 多相流 学术会议论文 编号：116086 基于VOF模型的水平管道油水两相流动数值模拟 刘晓燕 陈伟 刘立君 刘殿玮 东北石油大学土木建筑学院，黑龙江大庆，163318 （Tel：0459-6503441 E-mail：Liu_xydq@163.com） 摘要：由于油水两相间存在相界面，在不同的条件下可形成分层流、泡状流、段塞流、环状流等多种流型。本文采用VOF多相流模型对水平管道中高含水期油水两相流动进行数值模拟，得出不同体积含油率、流速及粘度条件下油水两相在水平管道中流动相分布图。结果表明，当体积含油率在0.1~0.2、流速在0.3~0.5m/s，粘度在0.02kg/(m·s)~0.06kg/(m·s)时，油水两相流动趋于分层流动，当流速大于0.5m/s时，油水两相由分层流动向分散流动转变；数值模拟结果与实验结果基本一致。 0 前言 在油田生产过程中，油田地面集输系统广泛采用油水混输管路。采用多相混输工艺技术具有降低工程投资、减少运行费用等单相管路不可比拟的优势。随着陆上油田的开采，大多数油田进入了特高含水采油期，集输管路中的油水两相流动形态发生了变化研究油水两相输送过程中的形态，对于集输管路的水力、热力计算，以及管路优化节能都具有重要的现实意义。 传统的油水两相流型的研究是以实验为基础的研究，通过配比特定比例的油水混合液，在不同的参数范围内测试两相的流动形态。这种方法运行周期长，搭建实验台复杂，而且运行费用高。近些年随着计算机技术的发展，数值模拟开始广泛地应用于流场分析研究中。计算流体力学(computational fluid dynamics．CFD)就是一种有效的研究流体动力学的数值模拟方法。顾汉洋等采用VOF模型对水平管道油水两相分层流动进行了数值模拟，分析了两相表观速度以及粘度、密度等物性参数对截面相含率和压降的影响。本文油水两相流动进行研究，利用VOF模型对水平管道中的油水两相流动进行模拟模型1.1 VOF模型 F模型是一种在固定的欧拉网格下的表面跟踪方法。当需要得到一种或多种互不相溶流体间的交界面时，可以采用这种模型。 (1) (2) 体积分数连续性方程 (2) 其中: (3) (3) 动量方程在VOF模型中，不同的流体组分共用着一套动量方程，计算时在流场的每个计算单元内，都记录下各流体组分所占有的体积率，从而确定界面的运动位置。 (4) 式中Fs为表面张力和壁面粘附作用引起的动量方程的源项： (5) 式中：为表面张力常数；为速度矢量；μ为流体动力粘度；ai为第i相的体积分数；k为表面曲率。 1.2 假设条件 建立能够真实反映水平管道中油水两相流动且具有足够计算精度的数学模型是数值模拟的关键。油水两相在管道截面上的许多参数分布不均，因此，油水两相在水平管道中的流动属于三元流动问题。但按三元流动对油水两相流型进行分析是非常困难的，模型的假设如下： 管道中的流体全部视为不可压缩流体； 不计流动过程中的能量损失；忽略温度变化对原油与水粘度、密度的影响； 4)油水两相在出口处不产生回流，即出口处表压为0.3计算模型及边界条件 本文采用管道尺寸为1000mm×mm的水平直管段对油水两相流动进行模拟，网格结构放大如图1所示，均为结构化网格。 图1 管道网格划分图 自由面重构格式采用较为精确的追踪自由表面的计算格式GeoReconstruct（几何重构)格式。考虑流动为湍流的情形，采用k湍流模型，近壁区流动采用壁面函数法进行处理。对流项按如下格式离散压力速度方式采用PISO(压力隐式算子分割)算法。压力方程采用PRESTO线性插值格式湍动能、湍动耗散率方程均采用阶迎风格式。 入口处设为油水混合的速度入口；出口处设定为压力出口边界条件；壁为无滑移。 计算结果及分析 不同含油率条件下流型变化规律 油水混合入口速度0.m/s（50t/d）、原油粘度0.0kg/(m·s)，含油率0.1、0.和0.进行不同含油率条件下流型变化规律的数值模拟，图2所示。 图2 随含油率变化的油水两相轴向体积相分布图 表1 随含油率变化的油水两相径向体积相分布图 距入口距离/m 体积含油率 0.2 0.5 0.8