井下泡沫发生器内流场数值模拟及结构优选

井下泡沫发生器内流场数值模拟及结构优选
闫月娟曹宇航桑小娜李默李森
(东北石油大学机械科学与工程学院)
摘要井下泡沫发生器的应用可以避免泡沫从地面注入井下过程中的消泡问题，改善泡沫驱油的效果。

但是由于井下径向空间有限，其内部泡沫发生结构的设计受到诸多因素的限制。

应用Fluent软件，对射流挡板式、螺旋挡板式和螺旋搅拌式3种不同的井下泡沫发生结构进行了内部流场的数值模拟，对比分析了流场的迹线、速度分布、气液相分布和压力分布情况。

结果表明：螺旋搅拌式发泡结构对气液相的掺混效果较好，压力损失小，为所研究的最佳井下泡沫发生结构。

关键词泡沫发生器气液两相流数值模拟结构优选
中图分类号TQ051.69文献标识码A 泡沫驱因具有波及体积大、驱油效率高的优点，已被我国各大油田普遍认可、广泛应用#T。

泡沫发生器是泡沫驱作业的重要部件，其性能好坏直接影响着泡沫驱的应用效果［4，&$。

目前，我国油田实施泡沫驱油作业时普遍采用地面发泡器发泡，将生成的泡沫通过管柱注入地层％这种注入方式,对于一些深井和超深井来说,存在严重的泡沫消泡问题，不能保证井下驱油泡沫的质量%井下泡沫发生器能够实现在井下产生泡沫，将所生成的泡沫直接注入地层，因而有效减少了泡沫的消泡，很大程度上改善了泡沫作业效果，具有良好的研究和开发价值％但是，由于井下径向空间狭小、操作局限，井下泡沫发生结构既要保证良好的发泡性能，同时还要保证径向尺寸满足要求且具有不容易发生堵塞的特点％
我国油田用泡沫发生器的结构形式多种多样，但大部分为地面发泡器％典型的地面泡沫发生器结构形式有涡轮式、螺旋式、同心管式、孔隙式和挡板式近些年来，一些专家学者也设计出了新型的泡沫发生器结构％如，林日亿设计了同心管与变向螺旋块组合形式的耐压耐酸泡沫发生器［7］；刘承婷等设计了螺旋挡板式冲砂洗井泡沫发生器相比于地面发泡器，对井下泡沫发生器的报道较少，仅在国内外的一些专利中有所涉文章编号0254-6094(2021)01-0089-06
及％为此，笔者在现有的泡沫发生结构的基础上，初选了3种结构简单、不容易发生堵塞的组合式泡沫发生结构，即射流挡板式、螺旋挡板式和螺旋搅拌式％应用数值模拟方法对三者的流体域进行数值模拟，对比分析不同泡沫发生结构内部流场的流动特性，优选出适于井下发泡的泡沫发生结构。

1物理模型及计算前处理
1.1模型建立及网格划分
根据笔者所选择的3种泡沫发生结构，结合目前已经优化好的泡沫发生器结构尺寸参数［9$11］，建立了相应的流体域模型，流体域模型的三维结构如图1所示％
射流挡板式泡沫发生结构由喷嘴、接受室、喉管、扩散管和4块等距放置的半圆形挡板组成；螺旋挡板式泡沫发生结构由焊接在固定杆上的螺旋叶片和4块等距放置的半圆形挡板组成；螺旋搅拌式泡沫发生结构由焊接在固定杆上的等距螺旋叶片和螺旋状布置且与固定杆横截面呈7.5。

左右升角的8对搅拌桨片组成％3种流体域模型的长度均为300mm,最大外径均为32mm；气液相入口形状和尺寸相同，气相入口直径为20mm，液相入口共8个，直径为5mm%各模型具体尺寸见表1~3。

作者简介：闫月娟(1971-)，教授,从事流体机械设计及多相流模拟仿真的研究和教学工作，******************
图1泡沫发生器流体域几何模型
表1射流挡板式流体域结构尺寸参数表
喷嘴直径接受室收缩角喉管长度喉管直径扩散管入口
扩散角
(°)
板间间距
挡板个数
挡板半径
mm(°)mm mm mm mm 10257020725415
表2螺旋挡板式流体域结构尺寸参数表
螺距mm 圈数
固定杆长度
mm
固定杆直径
mm
板间间距
mm
挡板个数
挡板半径
mm
253圈77.5925415
表3螺旋搅拌式流体域结构尺寸参数表
螺距mm 圈数
固定杆长度
mm
固定杆直径
mm
桨片对数
桨片宽度
mm
桨片升角
(°)
桨片间距
mm
相邻桨片
方向夹角
(°)
253圈1609897.51090
由于3种流体域几何模型结构较为复杂，因此利用ANSYS Meshing分别对3种流体域结构进行四面体非结构化网格划分。

为了保证网格质量,分别进行了网格无关性验证,当3种模型的网格数分别达到32.8万#33.4万和35.9万时,模型进出口的静压差值基本保持不变,网格质量满足数值模拟要求。

1.2边界条件及计算模型选择
以塔河油田某井实施泡沫驱的实际工况为例：井深40007,井下压力40MPa,井下温度140!。

向该井注入泡沫时,发泡液为水与耐高温起泡剂的混合基液(起泡剂质量浓度为3%),发泡液的注入流量为2.76m3/d,井下氮气的体积为发泡液的2倍。

笔者采用气液两相流的混合模型对发泡器内部的气液流场进行模拟'12(。

由于发泡液中发泡剂的质量浓度较小，因此以水相代替混合基液作为主相,次级相选择为氮气,设定气体粒径为100从7'13(。

通过查阅AP1700物质物性平台，该工况下水的密度为946.2037kg/m3,黏度为0.206547mPa・s,氮气的密度为263.8372kg/m3,黏度为3.00963"10B5Pa・s。

气液相的入口边界条件选择速度入口 ,入口速度根据实际工况下氮气与发泡液的流量均设置为4.88m/s，出口设定为自由出流。

考虑重力对发泡效果的影响，设置重力
加速度为9.81m/s2,方向为！正方向°
选用双精度压力基准算法隐式求解器稳态求解，湍流计算模型选择标准模型,SIMPLEC 算法用于进行速度压力耦合,采用有限体积法对控制方程进行离散，压力项采用PRESTO!方法,体积百分含量采用QUICK方法,动量、湍动能和湍流耗散率采用二阶迎风格式进行离散,壁面区域采用标准壁面函数处理，为无滑移边界条件,收敛精度设为10-6$14%0
2数值模拟结果分析
在泡沫发生器的设计过程中，应考虑以下几方面的因素:气液两相呈高速、紊态流动;尽量使得液体分散,以增大气液接触面积;适当地加大机械搅拌作用，使得泡沫的粒径更均匀。

因此，在泡沫发生器内部不仅要使得气液两相获得较大的流速，还要保证气液相之间获得足够的扰动,增加气液相之间的接触面积，以利于泡沫的形成。

2.1迹线分布
泡沫发生器内部流场的迹线分布如图2所示。

为了能够较为直观地展现迹线分布情况,分别取流体域全局、中分面和板间与桨叶间横截面上的迹线分布进行说明。

图2泡沫发生器内部流场的迹线分布
从射流挡板式的迹线分布（图2；）可以看出,气液两相经过挡板时,在挡板后部产生了较为明显的涡旋，但是涡旋内部的流线相对较少,流体大部分呈波浪状流态在板间流动,板间横截面上没有形成横向涡旋；从螺旋挡板式的迹线分布（图2b）可以看出,气液相经过螺旋叶片时形成了较为稳定的螺旋流动，螺旋流继续经过挡板时,流体呈螺旋流态在挡板间流动,从流体域中分面和板间横截面上的局部放大图可以看出,流体在挡板后部产生了较大的涡旋,涡旋内部的流线较多,说明大部分流体经过该涡旋实现混合,但在横截面方向仅有螺旋流动并没有形成多涡旋流动;从螺旋搅拌式的迹线分布（图2=）可以看出,与螺旋挡板式类似,气液两相在螺旋叶片位置形成了较为稳定的螺旋流,螺旋流继续经过搅拌桨片位置时,螺旋在多对桨片的分流作用下逐渐被削弱,在流体域的泡沫出口段,气液相呈微弱的螺旋流动。

从流体域中分面和桨叶间横截面上可以看出,气液相在该结构的桨叶前后轴向位置和径向位置上均形成了多个不规则的小尺度涡旋,这些旋涡造成的宏观不稳定性,可以加速气液间的混合,促进泡沫的形成。

2.2速度分布
泡沫发生器内部流场的速度分布如图3所示。

射流挡板式的速度分布如图3；所示,气体经过喷嘴时形成高速流,最大速度达到23.36m/s,带动接受室内的液体共同流入喉管。

流体在喉管位置由于存在紊流扩散作用,使得混合相的流速降低。

流体继续以垂直于挡板速度流经挡板时,由于受到挡板的阻碍作用，流体在挡板前的流速为11〜14m/s,在挡板后的流速为1〜2m/s,且低速区占到流体域的一半，挡板对流体的阻碍作用明显。

气液相在该结构的喉管和挡板位置的流速较大,因此湍流强度较高，但是气液两相在该结构内部并没有得到充分的扰动，速度分布极不均匀，因此不利于气液充分接触生成泡沫。

螺旋挡板式的速度分布如图3b所示，气液两相经过螺旋叶片后形成强螺旋流,速度最高达到16.25m/s,
且在螺旋流道内保持良好,使得螺旋流经过挡板时维持了一定的螺旋流速,相较射流挡板结构减小了挡板对流体的阻碍，挡板后低速区范围减小,流速影响范围更广,气液相在该位置形成了较为强烈的涡旋,有利于气液相间的接触。

螺旋搅拌式的速度分布如图3c所示，气液相在螺旋叶片位置的速度分布与螺旋挡板式基本相同。

螺旋流经过桨片位置时,由于桨片面积较小,且呈螺旋状多对错落布置,进一步减小了对流体的阻碍作用,使得流体保持了一定的螺旋流速,相邻两对桨叶前后流体最大速度的差值仅为1.5m/s 左右。

桨片对流体的分流作用更加明显,促进了桨叶间小尺度涡旋的形成。

同时,该结构出口处的速度分布相比于其他两种分布的更加均匀，因此有利于泡沫在管柱中的输送。

c.螺旋搅拌式
图3泡沫发生器内部速度矢量分布
2.3气液相分布
泡沫发生器内部气液相分布的均匀程度可以反映泡沫发生结构对气液的掺混效果。

泡沫发生器内部流场中的气相体积分数分布如图4所示,其分布规律验证了迹线和速度场分析的正确
性。

Vblume fraction(phase-2)
0.000.150.300.450.600.750.90 1.00
a.射流挡板式
b.螺旋挡板式
c.螺旋搅拌式
图4模型中分面上的气相体积分数分布
从图4可以看出射流挡板式的流场中存在较为明显的气液交界面，气液相没有很好地掺混。

气相射流后受挡板扰流呈波浪式流动,在挡板后因存在涡旋流动气液相混合加强,但因为流速和气量较少对混合的贡献较少,气液相混合效果不佳。

螺旋挡板式对气液相的掺混效果较好,流体在螺旋叶片位置由于螺旋流的流速较高,湍流脉动较强,具有一定混合效果。

流体继续经过挡板位置时,由于螺旋流的存在和挡板间所形成的涡旋促进了气液相间的混合,气液分布逐渐均匀,较射流挡板式有了很大改善,在尾部挡板的后面流体的气相体积分数基本不发生变化。

相比于螺旋挡板式,螺旋搅拌式对气液两相的掺混效果更好，流体经过螺旋叶片之后,相较于螺旋挡板式结构,在多组搅拌桨片的剪切搅拌作用下形成了更多不规则的小尺度涡旋，且流动更加均匀,促进了气液混合,使得气液相在前4对桨叶位置已经掺混得比较均匀。

后4对桨叶可以对所形成的泡沫进一步剪切细化,提高泡沫质量。

3种泡沫发生器在出口位置的气相体积分数分布对比曲线如图5所示。

从图中可以看出,射流挡板式出口位置气相体积分数波动最大,螺旋 搅拌式出口截面的气相体积分数波动最小,更加接近66.7%(气液比2：1)。

表明螺旋搅拌式相比于其他两种结构的发泡效果更好。

图53种泡沫发生器在出口位置的气相
体积分数分布
2.4压力分布
泡沫发生结构内部流场的静压分布如图6所示。

从图中可以看出,射流挡板式泡沫发生结构在喷嘴位置的节流效应较强,气体的压力损失
较大。

在喷嘴出口位置，由于气体的流速高,使得喉管和液相入口的静压值相对较低,更加有利于液体的进入。

在扩散管至第1块挡板位置,由于流道的扩大和挡板的阻挡作用,使得流体的流速降低,静压值增大。

流体继续经过其他挡板时，静压力逐渐降低。

螺旋挡板式和螺旋搅拌式由于在螺旋叶片位置的结构相同，因此压力变化相似!但是，由于螺旋搅拌式的搅拌桨叶呈螺旋状布置,且桨片尺寸较小,相比于螺旋挡板式的挡板结构对螺旋流所产生的阻碍作用更小，因此压力损失低。

图6模型中分面上的静压分布云图
泡沫发生器沿轴向方向上静压值的变化曲线如图%所示。

从图中可以看出,射流段压降为0.06MP*,射流结构每级挡板压降约为0.02MP*,总体压力损失最大，为0.12MP*；螺旋挡板式和螺旋搅拌式由于螺旋叶片结构相同，因此气液相在%=100mm以前静压值的变化趋势相似，压降约为0.035MP*。

但是在%"100mm以后,由于螺旋挡板相较于搅拌桨对流体的阻碍作用更强,螺旋挡板式的压力损失较大,螺旋搅拌式的压力损失较
图7泡沫发生器沿轴向方向上静压值的变化曲线小。

螺旋挡板式结构的整体压力损失为0.11MP*,螺旋搅拌式结构的整体压力损失为0.09MP*。

3结论
3.1单一的射流或者螺旋结构都不能使得气液两相得到有效混合。

气体介质在射流结构的喷嘴位置节流效应较强,高速的气体射流使得气液之间存在明显的交界面，不利于气液相间的掺混。

螺旋结构可以使得气液相形成稳定的螺旋流态,但由于螺旋流道内没有其他的扰流部件，因此气液相间没有得到有效的掺混。

3.2射流挡板式对气液相的掺混效果较差，流体在挡板间主要呈波浪状流态流动,在挡板后部涡旋位置的流体较少，流体没有得到充分的混合。

螺旋挡板式和螺旋搅拌式对气液相的掺混效果较好。

气液相经过螺旋叶片之后形成高速螺旋流,螺旋流继续经过挡板位置时产生了较大的涡旋,涡旋内部的流线较为密集,说明大部分流体经过该涡旋实现混合;高速螺旋流经过搅拌桨片时保持了一定的螺旋流速，同时气液相在桨片位置形成了多尺度、不规则的涡旋,促进了气液相间的混合。

3.3螺旋搅拌式相比于其他两种泡沫发生结构,在具有良好的发泡效果的同时，压力损失小。

且在前4对桨片就已经将气液掺混得比较均匀,其余多对桨叶的搅拌剪切可以使得所形成的泡沫细小致密以提高泡沫质量，因此更适合作为井下泡沫发生器的泡沫发生结构。

参考文献
[1]李兆敏,徐正晓,李宾飞,等.泡沫驱技术研究与应用
进展[J].中国石油大学学报（自然科学版）,2019,43
（5）：118〜127.
[2]王其伟.泡沫驱油发展现状及前景展望[J].石油钻采
工艺,2013,35（2）：94〜97.
[3]Sun L,Wei P,Pu W F,et al.The Oil Recovery Enhan­
cement by Nitrogen Foam in High-Temperature and
High-Salinity Environments[J].Journal of Petroleum
Science and Engineering,2016,147：485〜494.
[4]李宾飞.氮气泡沫调驱技术及其适应性研究[D].青
岛：中国石油大学（华东）,2007.
[5]何金钢.泡沫物理性能表征和泡沫驱油效果研究
[D].大庆：东北石油大学,2015.
（下转第118页
）
Leakage[J].Applied Mechanics and Materials,2013,
423-426：2078〜2081.
[16]Kimura T,Tanaka T,Fujimoto K,et al.Hydrodynamic
Characteristics of a Butterfly Valve—Prediction of Pressure Loss Characteristics[J].ISA Transactions, 1995,34(4)：319〜326.
(收稿日期$2020-04-01,修回日期$2021-01-13)
Interference Analysis of Truss Over-flow Triple
Eccentric Butterfly Valves
WEI Yan-qiang,HE Shi-quan,WU Qian-qian,LIU Shuai-shuai
(School of Mechanical and Power Engineering,Nanjing University of Technology)
Abstract Basing on the relationship of various eccentric parameters,a formula for calculating the area where three eccentric butterfly valve*s center of rotation stays was derived;and based on the MATLAB GUI module, a program was specially developed to check whether the three-eccentric butterfly valve interferes so as to pro­vide a method for quickly checking whether various types of three-eccentric butterfly valves interfere.The re­sults showed that,this software can accurately confirm whether the three-eccentric butterfly valve interferes. Key words triple eccentric butterfly valve,center of rotation,eccentricity parameter,MATLAB,interfer­ence analysis
（上接第93页）
[6]申瑞臣.泡沫发生器结构设计综述[J].石油机械,
1993,21（5）：52〜55.
[7]林日亿.蒸汽-氮气泡沫体系流动规律研究及应用
[D].青岛冲国石油大学（华东），2007.
[8]刘承婷，张维薇，刘钢，等.螺旋挡板式泡沫发生器的
设计及内部流动特性研究[J].数学的实践与认识,
2018,48（8）：120〜127.
[9]惠恒雷.射流发泡制造微气泡技术试验研究[D].青
岛：中国石油大学（华东）,2011.
[10]张维薇.多级挡板扰流泡沫发生器流场分析及实验
研究[D].青岛：中国石油大学（华东）,2019.
[11]Chen J S.Property Experiments on the Foam Genera­
tor and Its Influencing Factors during Down-the-Hole
Drilling[J].Process Safety and Environmental Protec-
tion,2018,114：169〜178.
[12]任万兴.煤矿井下泡沫除尘理论与技术研究[D].徐
州：中国矿业大学,2009.
[13]黄浩,李华斌,牛忠晓,等.高压下泡沫性能参数的
研究[J].科学技术与工程,2013,13（3）：694〜696,
710.
[14]闫月娟，王尊策,徐艳,等.采出液粘度对井下旋流
除砂器内流特性及分离性能的影响[J].化工机械,
2014,41（4）：505〜509.
（收稿日期:2020-05-07，修回日期:2021-01-11）
Numerical Simulation and Structural Optimization of Flow
Field in Downhole Foam Generator
YAN Yue-juan,CAO Yu-hang,SANG Xiao-na,LI Mo,LI Sen
(School of Mechanical Science and Engineering,Northeast Petroleum University)
Abstract Application of downhole foam generator can solve the defoaming while injecting the foam downhole from the surface and improve the effect of foam flooding.The downhole*s limited radial space and other factors jointly influence the design of internal foam generator.In this paper,having Fluent software used to simulate internal flow field of three types of downhole foam generators such as jet baffle type,screw baffle type and screw stirring type and having trace line,velocity distribution,gas-liquid distribution and pressure distribu­tion of the flow field compared and analyzed to show that,the spiral-stirred foaming generator has better mix­ing effect on gas and liquid and less pressure loss.
Key words foam generator,gas-liquid flow,numerical simulation,structure optimization。