首页/文章/ 详情

积鼎多相流仿真:石油及天然气工业瞬态多相流动模拟

1年前浏览4816

多相流存在于很多工业过程中,特别是石油工业,石油、天然气和水经常一起生产、运输。另外,当有外部热源,相间的热交换或由于压降而产生的闪蒸现象时,各相的相对体积分数可沿管道变化。某些情况下,这些流体动力学特征是不希望发生的,尤其是在烃类运输系统中,例如团状流可能对某些操作组件有害。这种多相流也存在于进出储层的油气管道中。此外,管道中多相流的复杂性随气体管道中固体颗粒(包括沙子、黑粉)的存在而增加。碳钢制成的石油和天然气管道中经常发生颗粒诱发的腐蚀,腐蚀后需要拆除受影响的管段,从而导致额外成本和生产中断。对于这些现象,需要高保真的CFD预测方法与实验结合起来为安全防范系统做准备。  


 多孔介质中的多相流动 

在CCS或EOR技术中,多孔介质流动通常涉及多个液相或气相,如:向油藏中注入水、蒸汽或二氧化碳可以提高石油的采收率。从实际角度来看,计算预测工具在选择合适的CCS地点、估算二氧化碳储存潜力、量化被捕获流体的数量方面是至关重要的。

VirtuaFlow有效精确的数值方法可以处理非线性现象,如相界面的传播及溶解。此外,可以对液体弯液面的演化过程、二氧化碳溶解到盐水相以及重力指进进行计算,可以实现在孔隙尺度上进行仿真。
图片
图1 通过柱状多孔介质阵列的两相流动:界面追踪方法

VirtualFlow使用level-set界面跟踪方法进行仿真,以捕获指进不稳定性。图1展示了VirtualFlow软件计算得到的二维重力指进的不稳定性,两种不同密度的液体在重力作用下通过一列柱体流动。红色的浓稠液体以手指形渗入较轻的蓝色相。大量的较轻液体残留在圆筒阵列内,限制了回收过程的有效性。


 液化天然气罐内的分层现象 

储罐内的液化天然气(LNG)通常因不同密度引起液体分层。而LNG分层可能会导致储罐翻覆和罐内压力急剧上升等问题,这会对LNG安全灌装作业造成危害,甚至破坏流动保障。为了防止出现LNG罐内的天然气泄漏导致的气体蒸发排放、燃料流量波动和组分变化等现象,LNG站需要进行合理的设计。

液化天然气储罐有两种类型的常用加注喷嘴,分别应用在储罐顶部和底部。加气点的选择由待罐装液化天然气密度与罐内(或根部)液化天然气密度相比的大小决定。一般来说,对于轻质天然气采用底部填充的方式,而重型液化天然气采用顶部填充的方式。分层过程会受到喷嘴、初始密度差、跟部的深度和填充速度的影响。

图片
图2 注入过程中轻质液化天然气的体积分数

注入处可以观察到大尺度的流体运动现象,而根据初始温差和注入速率的不同,该运动尺度可达容器尺度的一半。这些活跃的流体现象往往发生在混合剧烈的阶段,并推迟罐体液体分层现象的稳定时间。如图2所示,从VirtualFlow计算得到的速度场上可以看到这些大尺度运动。该流体运动有效地混合了不同密度的液化天然气,避免储罐内出现分层。


 气举 

在气举技术中,通常会在水和油流通的生产管道底部注入气体,以降低井内的重力压降,从而提高管道内的油流量。实际生产中,气体是通过管壁上的阀门注入的,容易在管道内产生大气泡。而有关水和空气的实验结果表明,注入小气泡可以提高气举效率。因为小气泡的上升速度较低,管中停留时间长,因此流体空隙率高,并且小气泡均匀地分布在管道的截面增加了气体空隙率,小延缓了气泡流向段塞流的转变,重力压降减少。因此产生大气泡对气举来说,是一个不利条件。

采用VirtualFlow软件对垂直立管底部的注入气泡(不考虑石油)进行模拟。气体以不同的流速从底部注入(如图3所示),然后根据气体在喷嘴处不同的表面速度值计算压降。
图片

图3 在垂直管道中注入气体,喷嘴处的表面速度分别为0.2、0.5和1m /s,黑色的等值线表示气泡

图片
图4 压降随喷嘴表面气速的变化曲线

图3所示的不同气体流量的结果表明,该模型能够预测不同大小的气泡在管道中以不同的速度、沿着特殊之字形路径上升。从图4可以看到CFD预测与实验结果是一致的——即使注入相对少量的气体,也会显著降低立管内的压降。


 垂直管道中的多相流 

管道中的气液流动在石油分离和运输中具有重要意义。气液混合物(油,固体颗粒,水合物,蜡,冷凝水和/或水的轻和重组分)会在各种流型(例如气泡,团状,环形,薄雾)下产生并随流体输运(图5)。另外,相含率会由于加热、相间热交换或由于减压而产生的闪蒸等沿着管道变化。在垂直管道流中(例如立管),流态识别(多达三个主相,当存在沙子和水合物时,相会更多)对于成功钻井和生产至关重要。管道中的多相流建模的主要任务是确定流态图。
图片
图5 垂直管内流态随气体表面速度的变化

目前对管道多相流和传热的预测工具都是基于双流体六方程模型,需要求解各相的守恒方程。在石油和天然气行业,这个模型被简化为一维,通常被称为“机械模型”。为了封闭机械模型方程,需要考虑流动特性如局部速度、壁面剪切应力、液体滞留量等的关系。这些关系在封闭模型中具有很大的不确定性,通常是经验模型,并且基于过度简化的假设。

  • 3D段塞流模拟
采用VirtualFlow软件对段塞流进行三维仿真,该计算在128个处理器的高性能计算机上进行。流型变化会影响时间步长,如当出现小气泡时,时间步长会减小至10-5秒。
图片
图6 三维LEIS模拟(Level Set+LES)过渡到段塞流
图6清楚地表明,采用VirtualFlow软件进可以得到实际中可能出现的流动图像。可以看到,不同大小和伸长的段塞自然形成而彼此没有碰撞,它们占据整个管道并向上移动。大量的气泡也在尾流中产生,聚集在泰勒气泡之间的区域。
  • 3D环状流模拟结果

对于环形流态,管道尺寸减小为3cm×100cm。计算域网格有80万个网格点,模拟计算需要使用超算的128核。采用VirtualFlow软件得到的环状流模拟结果如图7,展示了由密度等值线刻画的两相流动界面和横流速度矢量图。

图片
图7 环状流动状态下的三维LEIS模拟,三个不同的时间显示的流动分布
在不同的截面上,水膜厚度的变化如图8所示。
图片
图8 环状流动状态下的三维LEIS模拟(三个不同横截面的流动分布)
图8显示了管道的三个横流位置,表现出了液膜变形的特征。值得注意的是,这些流动具有一定的径向一致性。


仿真体系通用多相流积鼎 CFD
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2023-07-19
最近编辑:1年前
积鼎科技
联系我们13162025768
获赞 108粉丝 110文章 302课程 0
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈