0 引言
在油气井钻完井过程中,经常需要进行注水泥施工。以完井试油为例,利用射孔技术射开一层取得合格的试油资料后,若该层暂不生产或无工业油流,需要注水泥塞将该层封闭,上返至第二层进行试油。在注水泥塞施工过程中,若施工时间超过水泥浆的初凝时间,水泥浆还没有完全注入井内就失去了流动性,容易发生“插旗杆”等事故,导致施工失败[1]。一般采取两种方式解决这个问题:①提高水泥浆的性能,延长初凝时间[2-3];②提高施工效率,缩短施工时间[4]。
在整个注水泥塞施工中,地面调配水泥浆的环节费时费力、施工难度大、安全风险高,是提高施工效率的主要障碍。目前,井场需要人力将袋装干水泥粉倒入1.5 m3的水泥浆罐,同时泵车向刺枪供水冲击干水泥粉,施工人员利用铁锹搅拌调配水泥浆。该方法存在以下缺点:①配灰时间长,用G级水泥配制1 m3密度1.85 g/cm3的水泥浆耗时约40 min,占总施工时间的
目前针对提高水泥浆性能的研究较多,但对提高注水泥塞施工效率的研究较少,尤其是对水泥浆调配装置的研究鲜有报道。为此,本文根据射流抽吸原理设计了射流抽吸式水泥浆调配装置,基于离散颗粒法对注入漏斗内干水泥粉的运动规律进行了数值模拟,优化设计了注入漏斗的半顶角和调制开关塞的安装高度,以达到防堵塞的目的,最后通过现场应用验证了水泥浆调配装置的使用效果。
1 装置整体设计射流抽吸原理即具有一定压力的工作流体通过喷嘴高速喷出,使压力能转化为速度能,在喷嘴出口区域形成真空,从而将被抽介质吸引入工作流体,完成两种介质的混合后排出[5-6]。射流抽吸式水泥浆调配装置主要包括注入漏斗、开关搅拌器、绞龙输送机、射流抽吸器和水泥浆罐等,如图 1所示。其关键部件结构如图 2所示。
开关搅拌器上等距离120°螺旋焊接三根搅拌杆,底部为锥形调制开关塞。在调配水泥浆之前,绞龙输送机将干水泥粉运送至注入漏斗内。当需要打开调制开关塞注入干水泥粉时,逆时针旋转手轮,通过螺纹传动使开关塞向上移动,直到螺纹杆旋出,继续旋转手轮,搅拌杆旋转搅拌注入漏斗内的干水泥粉,同时调制开关塞改善干水泥粉的流动状态,保证干水泥粉快速顺利地被抽入射流抽吸器内。在调配水泥浆时,水流过射流抽吸器内的节流喷嘴,在节流作用下水流速度升高,在注入漏斗下部形成负压,产生射流抽吸作用,干水泥粉被抽吸后与水充分混合,以45°的倾角喷入水泥浆罐,形成冲刺力并带动水泥浆罐内的液体旋流,进一步调配水泥浆。在同样泵压和排量条件下,可通过更换不同内径的节流喷嘴调整注入干水泥粉的速度。
2 干水泥粉运动规律数值模拟在干水泥粉注入过程中,关键问题是防止干水泥粉堵塞注入漏斗口。为此,在开关搅拌器上设计了搅拌杆。注入漏斗的结构和调制开关塞的安装同样影响干水泥粉的流动性。因此,通过对注入漏斗内干水泥粉运动规律的数值模拟,分析干水泥粉运动的动态变化过程和在注入漏斗内的分布特征,据此优化设计注入漏斗的半顶角和调制开关塞的安装高度,达到防堵塞的目的。
2.1 理论计算模型为了更贴近于干水泥粉的运动状态,基于离散颗粒方法(Discrete Element Method,DEM)[7-8]将干水泥粉离散为三维受力球形颗粒,采用满足高浓度颗粒流的软球模型进行研究[9-10]。它的基本原理是将颗粒堆积体视为离散质元的集合,利用牛顿第二定律建立每个质元的运动方程,用动态松弛法迭代求解,从而求得颗粒堆积体的整体运动形态[11]。
2.1.1 力的计算模型干水泥粉之间的接触力采用弹簧-阻尼器-摩擦器的表达模式,合力由弹性力、黏滞阻尼力和摩擦力组成。DEM单元接触模型如图 3所示。计算过程中,法向弹性力采用Hertz接触理论模型确定,切向弹性力采用R.D.MINDLIN[12]理论确定。
(1) |
(2) |
式中:
黏滞阻尼力采用M.J.ADAMS等[13]提出的液桥黏性力的解析解,其表达式为:
(3) |
(4) |
式中:fvn为法向黏滞阻尼力,N;fvt为切向黏滞阻尼力,N;h为颗粒i与颗粒j之间的液膜宽度,m;vn为法向相对速度,m/s;vt为切向相对速度,m/s;η为流体黏滞系数,Pa·s;R为粒子平均半径,m。
2.1.2 运动方程的求解根据牛顿第二定律,结合颗粒的受力分析,建立颗粒的运动方程:
(5) |
式中:i=1, 2, 3,……,n;xi表示位移量,m;mi为颗粒i的质量,kg;∑Fi为粒子i所受合力,N。
在运动方程求解过程中,力与位移的计算在同一时间步内进行,即通过对某一时刻t的加速度用中心差分形式导出下一时刻t+Δt的位移,此时,颗粒重新承受接触力。根据t+Δt时刻新的受力和位移状态进行循环迭代,可获得颗粒群中每个颗粒的运动形态。
2.2 物理模型由于研究的重点是注入漏斗内干水泥粉的运动和注入规律,侧重于解决干水泥粉的堵塞问题,所以只建立注入漏斗内的物理模型,如图 4所示。半顶角的大小是影响干水泥粉运动规律的主要参数,需进行优化设计,设置为30°~90°。注入漏斗直径1.5 m,注入口直径0.05 m,高度0.8 m。
在数值模拟时,射流在注入漏斗口产生的抽吸压力可根据伯努利方程计算[14],工作压力为2 MPa,节流喷嘴直径为9.5 mm。
2.3 模拟结果分析 2.3.1 干水泥粉的运动规律由数值模拟显示不同时刻注入漏斗内干水泥粉的位置分布。为追踪干水泥粉的运动轨迹和位置,在4个不同高度处设置示踪层颗粒(蓝色颗粒)。图 5为注入漏斗半顶角45°时的数值模拟结果。
注入漏斗内的干水泥粉在注入过程中,漏斗中间部位的干水泥粉流动较快,靠近注入漏斗壁面的干水泥粉流动较慢,示踪层呈现出抛物线的形状,顶部呈现出“中间低、边部高”的状态,表现出典型的漏斗流的流动形态。注入漏斗中间部位的流动可称为垂直流动区。垂直流动区两侧的干水泥粉在滑移作用下容易进入漏斗流动区域,可称之为快速流动区。注入漏斗上半部分靠近注入漏斗壁的流动称之为慢速流动区。在慢速流动区的下部靠近注入漏斗壁处的干水泥粉,受到半顶角的影响,向注入口滑移的坡度不同,加之干水泥粉间的摩擦力和黏性力,使得速度基本为0,称之为停滞区。
在漏斗流状态下,干水泥粉注入过程中,由于干水泥粉间的相对运动,引起切向摩擦,并且干水泥粉之间挤压错动,使得摩擦力增大。在慢速流动区和停滞区,干水泥粉间黏性作用力增强,容易聚合成团,引起干水泥粉间更强的相互挤压,产生压实作用达到力平衡,形成难以破开的机械拱,从而阻碍干水泥粉的流动,最终堵塞于注入漏斗口。
由此可见,漏斗流通过摩擦力和黏性力作用,使得干水泥粉聚结成拱是导致注入漏斗口堵塞的主要原因。因此,需通过优化设计注入漏斗的半顶角和调制开关塞的安装高度,改善干水泥粉的流动状态,由漏斗流向整体流转变,达到防堵塞的目的。
2.3.2 注入漏斗半顶角的影响规律图 6是半顶角为30°、45°和60°时注入漏斗内干水泥粉的流动状态。由图 6可知,随着半顶角的增大,干水泥粉的停滞区域增大。这主要是因为当半顶角较大时,干水泥粉整体的流动速度较小,易搭接架桥,黏性力、干水泥粉间摩擦力以及干水泥粉与注入漏斗壁面的摩擦力使得干水泥粉的注入动能耗散,尤其在注入漏斗壁面附近的干水泥粉,本身流动速度较低,更易造成流动停滞。
提取半顶角为30°、45°、60°、75°和90°时的干水泥粉体积流率,结果如图 7所示。对于半顶角较小的注入漏斗,干水泥粉的体积流率较大,完全注入所需的时间较短。但是,当半顶角为30°时,干水泥粉的注入出现明显异常。当干水泥粉刚注入时,体积流率迅速上升,到达40 s时,体积流率突然下降,并在之后出现转折,达到新的峰值后体积流率又突然下降为0。这说明在注入过程中,注入漏斗口附近结拱堵塞,最终导致干水泥粉无法注入。分析认为,由于注入漏斗的锥体部分逐渐收缩,导致干水泥粉注入时重新排布,产生严重的挤压错动,摩擦力增大,导致锥体部分干水泥粉流动困难,而上部干水泥粉受到的影响较小,注入速度较快,从而压实锥体部分干水泥粉,最终导致干水泥粉紧密接触而堵塞注入漏斗口。
对比不同注入漏斗半顶角时的体积流率发现:随着半顶角的增大,体积流率逐渐减小;当半顶角为45°时,体积流率最大,表明此时干水泥粉注入的流动过程相对比较稳定。因此,优选注入漏斗的半顶角为45°。
2.3.3 调制开关塞的影响规律在仓体中引入圆锥形调制元件,能够降低颗粒流注入的中间流速,使得整体流速均匀[15-16]。对于干水泥粉的注入,在注入漏斗中设置圆锥形调制开关塞,建立如图 8所示的物理模型。安装高度h1是影响干水泥粉注入的主要因素,需优化设计,设置为70~100 mm。调制开关塞的顶角为90°,侧翼长度a为40 mm。
在注入漏斗设置调制开关塞后,干水泥粉的流动形态如图 9所示。与图 5对比可见,设置调制开关塞后,注入漏斗内的干水泥粉呈现明显的整体流动形态,中间部位干水泥粉的快速流动得到了较好的控制,无明显的沉降漏斗出现,示踪层的干水泥粉基本呈水平分布均匀下降,干水泥粉间的相对运动减少,摩擦力减小,注入漏斗壁面附近原处于静止状态的干水泥粉开始运动,不存在死料的现象,流动停滞区得到明显改善。由此可见,调制开关塞能够实现干水泥粉流动形态由漏斗流向整体流转变,从而有效解决堵塞问题。
调制开关塞的安装高度对干水泥粉流动状态的改善效果至关重要。提取安装高度h1为70、80、90和100 mm时的干水泥粉体积流率,如图 10所示。
由图 10可见,随着安装高度增加,干水泥粉的体积流率逐渐减小,即干水泥粉的注入速度逐渐减小。因此,优选调制开关塞的安装高度为70 mm。
3 现场应用 3.1 现场布置及操作射流抽吸式水泥浆调配装置的现场布置如图 11所示。工作过程为:①启动绞龙输送机,将干水泥粉输送至注入漏斗内备用;此时阀门3和阀门1处于关闭状态,阀门2和阀门4处于开启状态。②开水泥浆泵,同时旋转开关搅拌器,打开调制开关塞。水由清水罐经过阀门4,由泵入口4进泥浆泵,由泵出口泵入射流抽吸器,将注入漏斗内的干水泥粉抽吸入水中,混合后喷入水泥浆罐。③达到设计清水量时,先打开阀门3,然后关闭阀门4,水泥浆罐内的水泥浆被水泥浆泵吸入泵内后,经阀门2、射流抽吸器后重新喷入水泥浆罐,形成闭循环,进一步调配水泥浆。④水泥浆完全调配均匀后,先打开阀门3,然后关闭阀门2,水泥浆罐内的水泥浆经过井口装置泵入井内,进行注水泥施工,直到所有水泥浆注入完毕。⑤停水泥浆泵,施工结束。
3.2 应用实例
大港油田女1701井完钻井深3 725 m,注水泥封Mz层位135号层,上返Ek2层位111-112、114-115号层。使用ø73.0 mm(
利用射流抽吸式水泥浆调配装置在井场调配水泥浆,现场只需水泥浆泵操作工和配浆工两人配合施工,20 min后水泥浆调配均匀。调配过程中,整套装置工作状态良好,干水泥粉注入稳定,未出现堵塞现象。水泥浆注入井内后关井候凝72 h后,清水对水泥塞试压15 MPa,30 min压降不超过0.5 MPa,试压合格。
射流抽吸式水泥浆调配装置在大港油田应用30多井次,调配水泥浆用时节省50%以上,水、灰混合均匀,极大降低了“插旗杆”等施工事故发生的概率,降低了施工劳动强度和安全风险,彻底解决调配水泥浆施工过程中存在的费时费力的问题,从而促进施工作业安全、高效地进行。
4 结论(1) 根据射流抽吸原理设计了射流抽吸式水泥浆调配装置,通过射流抽吸器内的节流喷嘴在注入漏斗下部产生抽吸作用,将干水泥粉注入并与水充分混合,可快速调配出合格密度的水泥浆。
(2) 基于离散颗粒法对注入漏斗内干水泥粉的运动规律进行数值模拟研究,优化设计了注入漏斗的半顶角和调制开关塞的安装高度,以达到防堵塞的目的。
(3) 优化后的注入漏斗半顶角为45°,调制开关塞的安装高度为70 mm。
(4) 射流抽吸式水泥浆调配装置可节省调配水泥浆用时50%以上,水、灰混合均匀,极大降低了“插旗杆”等施工事故发生的概率,保证了施工作业安全、高效地进行。
[1] |
王兴武, 王时成, 马小龙, 等. 超深井注水泥塞技术[J]. 钻采工艺, 2012, 35(4): 32-34. WANG X W, WANG S C, MA X L, et al. Cementing plug technique in ultra-deep well of western region[J]. Drilling & Production Technology, 2012, 35(4): 32-34. |
[2] |
徐璧华, 冯青豪, 谢应权, 等. 考虑循环温度影响下注水泥流变性能计算方法[J]. 钻井液与完井液, 2017, 34(6): 79-82. XU B H, FENG Q H, XIE Y Q, et al. A method of calculating rheology of cement slurry affected by circulation temperature[J]. Drilling Fluid & Completion Fluid, 2017, 34(6): 79-82. |
[3] |
陶谦, 周仕明, 张晋凯, 等. 水泥浆流变性对水平井固井顶替界面的影响:基于天河一号大规模集群计算平台的数值模拟[J]. 石油钻采工艺, 2017, 39(2): 185-191. TAO Q, ZHOU S M, ZHANG J K, et al. Effect of rheological property of slurry on cementing displacement interface of horizontal well:the numerical simulation based on large-scale cluster computing platform Tianhe-1[J]. Oil Drilling & Production Technology, 2017, 39(2): 185-191. |
[4] |
南树歧, 刘小卫, 王峰琳, 等. 高温高压生产气井的连续管注水泥作业[J]. 石油机械, 2012, 40(11): 31-33. NAN S Q, LIU X W, WANG F L, et al. Cementing through high pressure CT for HTHP gas well[J]. China Petroleum Machinery, 2012, 40(11): 31-33. |
[5] |
黄涛, 李军, 宋学锋, 等. ReelWell钻井反向射流抽吸技术及工具研究[J]. 石油机械, 2019, 47(2): 8-13. HUANG T, LI J, SONG X F, et al. Research on reelwell drilling reverse jet technology and tools[J]. China Petroleum Machinery, 2019, 47(2): 8-13. |
[6] |
王方祥, 王瑞和, 周卫东, 等. 前混式磨料射流发生系统混合腔的优化设计[J]. 石油机械, 2015, 43(12): 19-24. WANG F X, WANG R H, ZHOU W D, et al. Optimization design of mixing chamber of the pre-mixed abrasive water jet generating system[J]. China Petroleum Machinery, 2015, 43(12): 19-24. |
[7] |
CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Geotechnique, 1979, 29(1): 47-65. |
[8] |
MARKAUSKAS D, RAMIREZ G, KACIANAUSKAS R, et al. Maize grain shape approaches for DEM modeling[J]. Computers and Electronics in Agriculture, 2015, 118: 247-258. |
[9] |
WALTON O R, BRAUN R L. Stress calculations for assemblies of inelastic spheres in uniform shear[J]. Acta Mechanica, 1986, 63(1): 73-86. |
[10] |
谭骏华, 罗坤, 樊建人. 软球模型在颗粒流全尺度模拟中的验证和分析[J]. 浙江大学学报(工学版), 2015, 49(2): 344-350. TAN J H, LUO K, FAN J R. Verification and analysis of soft-sphere model in fully-resolved simulation of particulate flow[J]. Journal of Zhejiang University (Engineering Science), 2015, 49(2): 344-350. |
[11] |
陈磐, 边彩霞, 崔晓波, 等. 料仓中湿颗粒的卸料流型[J]. 中国粉体技术, 2019, 25(1): 1-6. CHEN P, BIAN C X, CUI X B, et al. Flow pattern of wet particles discharging from a hopper[J]. China Powder Science and Technology, 2019, 25(1): 1-6. |
[12] |
MINDLIN R D, DERESIEWICZ H. Elastic spheres in contact under varying oblique forces[J]. Journal of Applied Mechanics, 1953, 20(3): 327-344. |
[13] |
ADAMS M J, PERCHARD V. The cohesive forces between particles with interstitial liquid[J]. Institution of Chemical Engineering Symposium, 1985, 91: 147-160. |
[14] |
张兆顺, 崔桂香. 流体力学[M]. 3版. 北京: 清华大学出版社, 2015. ZHANG Z S, CUI G X. Fluid mechanics[M]. 3rd edition. Beijing: Tsinghua University Press, 2015. |
[15] |
谭援强, 肖湘武, 郑军辉, 等. 锥形改流体下部孔径对筒仓卸料流态的影响[J]. 农业工程学报, 2016, 32(19): 82-87. TAN Y Q, XIAO X W, ZHENG J H, et al. Effect of outlet diameter of cone-in-cone insert on silo flow pattern[J]. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(19): 82-87. |
[16] |
谢宇, 韩保红, 程兆刚, 等. 料斗结拱及改流体的数值仿真研究[J]. 计算机仿真, 2013, 30(8): 221-224. XIE Y, HAN B H, CHENG Z G, et al. Numerical simulation study on arch in hopper and cone-in-cone insert[J]. Computer Simulation, 2013, 30(8): 221-224. |