0 引言
旋切式PDC钻头作为一种新型PDC钻头,它为页岩气开采过程中遇到的问题提供了一种新的解决方法。钻头结构包括切削结构和水力结构,其水力结构影响钻头的水力性能。
国内外不少学者已经对PDC钻头的水力结构开展了研究。管志川等[1-2]采用试验方法对安装矩形喷嘴及圆形喷嘴的PDC钻头的井底流场进行研究,且用数值模拟方法对圆形喷嘴射流冲击井底的问题进行研究和试验验证。刘刚等[3]利用试验装置研究了0°~20°喷嘴角度对PDC钻头井底清洗能力的影响。谢翠丽等[4]利用数值方法对喷嘴垂直和倾斜射流两种情况进行了研究。况雨春等[5-6]利用数值模拟方法研究了喷嘴的方位角、喷距和倾斜角等因素对PDC钻头井底流场的影响。R.CROUSE等[7]首次采用数值计算方法对井底流场进行分析,并对PDC钻头水力结构进行了优化。T.NAZARI等[8-9]运用DPM模型来研究井底流场岩屑运移比和岩屑尺寸与喷嘴尺寸的线性关系。拉夫堡大学的G.R.WATSON等[10]使用数值模拟方法对直径215.9 mm的PDC钻头进行了水力结构优化设计研究。
旋切式PDC钻头比常规PDC钻头的结构更复杂,工作原理也比较特殊,常规的数值模拟方法不能全面描述旋切式PDC钻头水力结构及井底流场特性。本文基于CFD软件Fluent对旋切式PDC钻头进行有限元计算,分析了刀盘结构形状、喷嘴结构、喷嘴弯曲曲率及中心喷嘴和边缘喷嘴直径组合参数对井底岩屑颗粒质量浓度的影响,同时分析了优化后的刀盘和喷嘴结构及中心喷嘴和边缘喷嘴直径组合对井底清洗效果的影响。所得结果对未来页岩气开采方式的研究具有一定的指导意义。
1 旋切式PDC钻头结构及仿真分析 1.1 三维结构旋切式PDC钻头主刀翼在破岩的同时旋切刀盘反向旋切破岩,实现对井下地层机械交叉式切削破岩,其三维结构如图 1所示。
1.2 内部结构
旋切式PDC钻头主要由上部壳体、卡环装置、内齿轮部、小齿轮组、齿轮轴、密封装置、下部壳体、钻头主刀翼、旋切刀盘、连接装置和钻具接头等部分组成,其内部结构如图 2所示。
1.3 旋切式PDC钻头刀盘结构优化
优化前、后刀盘厚度和直径保持一致,其三维结构如图 3所示。
1.4 旋切式PDC钻头喷嘴结构优化
在等变速形变截面喷嘴的出口处增加了一个弯曲喷口,其高度为7 mm。中心喷嘴优化前、后三维结构如图 4所示。
1.5 仿真模型 1.5.1 物理模型和网格划分
利用SolidWorks软件导出旋切式PDC钻头井底流体域,并进行数值分析。采用ICEM软件对井底流体域三维模型划分网格,考虑管壁处黏性底层的影响,使用Inflation将边界层细化为5层,设生长因子为1.1。井底流体域包括流体入口、环空出口和井底流场以及各结构对应的网格模型,如图 5所示。
1.5.2 边界条件
液相为连续相,砂粒为离散相[11]。连续相采用标准k-ε湍流模型,近壁区域采用标准壁面函数处理。采用50 L/s速度入口边界及30 MPa压力出口边界,采用标准壁面边界。离散相通过Injection注射,粒子喷射选为表面Surface喷射,设置入射面为Inlet,粒子为尺寸均匀的石英石。分别设置粒子注入速度、颗粒直径和质量流量,其中粒子注入速度与连续相流入速度相同。在DPM选项卡中入口和出口采用Escape模型,固体颗粒与壁面碰撞时采用Reflect壁面反弹模型[12]。
动量和湍动能采用二阶迎风离散格式,选半隐式Simple算法,用该算法求解湍流问题可以得到较好的收敛计算结果。选标准初始化方法,并设置从入口初始化。
1.5.3 网格适应性分析网格形式和质量直接影响计算精度和计算量。进行网格适应性分析能够在减少计算量的同时获得较高的计算精度。分析入口速度50 L/s,质量流量0.01 kg/s,出口压力30 MPa时不同网格密度下井底同一点速度的变化情况,计算结果如表 1所示。
网格尺 寸/mm |
网格数量 | 网格变 化/% |
井底流速/ (m·s-1) |
速度变 化/% |
8 | 1 193 642 | — | 12.6 | — |
7 | 11 125 760 | 8.32 | 11.5 | 0.087 |
6 | 15 192 324 | 3.70 | 10.8 | 0.061 |
5 | 20 211 202 | 33.00 | 10.1 | 0.064 |
从表 1可以看出:随着网格的加密,井底流速不断变化,但是基本在12.6~10.1 m/s之间;网格每细化1 mm,网格数量逐渐增大,计算量逐渐增加,当网格尺寸从8 mm依次变到5 mm时,速度变化量依次为0.087%、0.061%和0.064%;当网格尺寸从6 mm变为5 mm时,网格数量增加了33.0%,但井底速度变化趋于稳定。为减少计算量,选用网格尺寸6 mm进行分析。
2 旋切式PDC钻头水力结构优化设计及模拟结果分析 2.1 旋切刀盘未优化井底流场分析当液固两相流速度为50 L/s、颗粒直径为100 μm、颗粒质量流量为0.005 kg/s、网格尺寸为8 mm、刀盘结构未优化时的仿真结果如图 6~图 8所示。选取1号、2号和3号喷嘴为研究对象,分析井底流动特性和流动规律。射流区指工作流体从钻头喷嘴喷出后形成的高速流体区域。漫流区指流体离开射流区,冲击方向与井底平行,形成层状且具有携岩和冷却作用的流体区域。上返区指井底漫流携岩到达井壁附近后,在井壁空间的限制下,从钻头排屑槽向上运动到达环空的区域。钻头喷嘴轴向剖面速度如图 6所示。
随着喷嘴倾斜角度的增大,流体获得水力能量分量增加,流量增大,而冲击方向的反方向流量较少,在高速射流带动下会形成漩涡区,如1号和3号喷嘴红色框选区域。在旋切刀盘处,刀盘附近淹没流体的射出和流出方向受力不均衡,会产生漩涡区,如2号喷嘴红色框选区域。大量漩涡区的存在会影响井底岩屑排出,导致钻头泥包,因此应尽量减小漩涡区域的面积,降低井底岩屑颗粒的质量浓度。
未优化旋切刀盘结构对井底影响如下。
(1) 刀盘上方淹没流体由于受力不均衡会产生比较多的漩涡区,导致井底岩屑不能及时排出,聚集在刀盘附近,引起井底岩屑的重复切削,产生钻头泥包。
(2) 部分喷嘴射流冲击在旋切刀盘上,导致射流的水力能量被分散,影响流体在漫流区的速度,同时射流在刀盘的边缘发生偏折,刀盘上方区域产生了一些漩涡区,使射流无法有效地清洗旋切刀盘区域的岩屑,也会增加钻头泥包。
(3) 高速射流会不断冲击旋切刀盘边缘,使刀盘受较大冲击力,刀盘安装轴受到弯矩作用,导致刀盘与轴的连接部分发生松动,使旋切刀盘不能有效地对井底岩层进行切削,旋切刀盘处速度如图 7所示。
利用DPM模型进行数值模拟。为了便于观察结果,选取1条典型的岩屑颗粒运动轨迹进行分析。岩屑颗粒的运动轨迹并不是二维运动轨迹,而是在流体冲击下与井壁和钻头体之间产生碰撞和反弹形成的不规则三维运动轨迹。运动开始时岩屑颗粒的速度较小,经过漫流区时速度逐渐增大,但由于在岩屑运移过程中会受到井壁和未优化旋切刀盘的作用而发生反弹,导致运动速度减小,井底岩屑颗粒输运轨迹如图 8所示。
2.2 旋切刀盘结构优化前、后井底流场分析旋切刀盘结构优化前、后钻头表面岩屑颗粒质量浓度如图 9所示。
从图 9可以看出,优化后的旋切式PDC钻头体表面的岩屑颗粒质量浓度大大降低,表明此区域流体的流动比较充分。刀盘轴处的岩屑黏附堆积情况有所改善,这是因为刀盘结构改变,有更多的流体可以通过刀翼之间的空间到达刀盘轴,同时减少了喷嘴射流与刀盘之间的干涉,井底漫流速度有所提高,有利于提高刀盘区域的冲洗能力。
在主刀翼上增大排屑槽区域面积,岩屑颗粒质量浓度也有所降低,这是因为井底漫流速度的提高,使得上返区流体速度也随之提高,从而增强了工作流体的携岩能力,可使岩屑快速离开排屑槽,进入环空区域。
旋切刀盘结构优化前、后钻头井底压力云图如图 10所示。图 10a中井底负压区面积较小,主要集中在钻头的中心部分,可能会造成刀盘区域岩屑颗粒质量浓度升高;而图 10b中的负压区域面积明显增大,在以射流冲击点为中心的区域形成了完整的压力梯度。井底压力的降低趋势代表了流体携岩的流动方向,有利于射流冲击井底时,井底漫流将岩屑带入环空处。
中心喷嘴结构优化前轴向剖面速度云图如图 11所示。从图 11可以发现,虽然井底漫流区的速度增加,但是上返区的速度却减小了,导致漫流携岩能力减弱。这是因为优化后刀盘的体积减小,井底漫流在经过刀盘下表面和井底平面形成的空间时,流体扩散,导致井底漫流水力能量被分散,漫流速度整体下降。为了改善这种情况,下面将对中心喷嘴(2号和5号喷嘴)进行相应的结构优化。
2.3 喷嘴结构优化后井底流场分析
中心喷嘴结构优化后轴向剖面速度云图如图 12所示。从图 12可以看出,弯曲喷口使喷嘴射流有了一个喷射角度,增大了在井底漫流的水力能量,提高了漫流速度,使漫流以较高速度进入上返区,提高了上返区举升高度,增强了该区域的携岩能力,降低了井底岩屑颗粒的质量浓度。
中心喷嘴优化前、后钻头井底压力云图如图 13所示。由图 13可以看出,中心喷嘴优化前、后井底平面压力保持不变,但弯曲喷嘴射流在井底平面产生的负压区域面积较大,同时可以将负压区扩展到井底平面边缘,更有利于该区域岩屑的排出。弯曲喷嘴在钻头的中心区域形成了一定的压力梯度,有利于钻头中心的岩屑清除,降低井底岩屑颗粒的质量浓度。
中心喷嘴优化前、后钻头表面岩屑颗粒质量浓度云图如图 14所示。由图 14可知,安装有弯曲喷嘴的旋切式PDC钻头,其旋切刀盘区域的岩屑颗粒浓度降低,说明弯曲喷嘴射流产生的井底漫流可以有效携带岩屑排出,缩短岩屑与旋切刀盘和钻头体表面的接触时间,从而减少了岩屑的黏附堆积。
为了更好地描述井底速度变化,在井底附近和旋切刀盘下表面附近各选一个位置,观察喷嘴射流从钻头中心到井壁的速度变化,结果如图 15和图 16所示。
从图 15和图 16可以看出,弯曲喷嘴的水力性能明显优于普通喷嘴,但两个喷嘴的速度变化趋势大致相似。带有弯曲喷口的喷嘴射流具有较大的水力能量,可以保持一个较大的井底漫流速度,在井底附近的最高速度为58.8 m/s,在旋切刀盘附近的最高速度为62.7 m/s。较高的井底漫流速度有利于岩屑及时排出,可降低井底岩屑颗粒的质量浓度。
3 改进后钻头的水力结构参数优化 3.1 中心喷嘴曲率优化由上述分析可知,弯曲喷嘴可以使流体改道,加快了井底流体上返速度。为了研究弯曲喷嘴曲率ρ对井底流场的影响,设计了中心喷嘴曲率分析方案,如表 2所示。
保持旋切式PDC钻头其他结构参数不变,根据钻头中心喷嘴曲率分析方案进行数值模拟分析,选出最优的中心喷嘴曲率。
不同曲率中心喷嘴轴向剖面速度云图如图 17所示。由图 17可见,随着喷嘴曲率ρ增大,漫流层厚度逐渐增加,漫流的水力能量也逐渐增大。漫流层在排屑槽中上返的速度越快,漫流层的举升高度越高。当曲率ρ=1/50时,速度达到最大值;当曲率ρ=1/30时,速度略有衰减。较大的井底漫流速度有利于降低井底岩屑颗粒的质量浓度。
为了更好地描述井底速度的变化趋势,在井底附近选取一个位置,统计喷嘴射流从钻头中心到井壁的速度变化,如图 18所示。从图 18可见,4条曲线整体变化趋势大致相似,当ρ=1/50时,可以保持一个较大的井底漫流速度,在井底附近的最高速度为69.4 m/s,在旋切刀盘附近的最高速度为68.9 m/s。因此,当ρ=1/50时,可以获得较为优良的井底流场特性,既保持了射流冲击方向的漫流速度,又使得射流冲击反方向的漫流速度不会过小,保证了钻头中心区域和旋切刀盘区域的岩屑有效排出,降低了井底岩屑颗粒的质量浓度。
不同曲率下钻头表面岩屑颗粒质量浓度如图 19所示。从图 19可见,当喷嘴曲率为1/150时,钻头中心区域岩屑颗粒质量浓度较低,但是外部区域岩屑颗粒质量浓度较高。当喷嘴曲率为1/30时,钻头外部区域岩屑颗粒质量浓度较低,但是中心区域岩屑颗粒质量浓度较高。当喷嘴曲率为1/50时,钻头中心区域和外部区域整体岩屑颗粒质量浓度适中,能够保证钻头中心区域和旋切刀盘区域的岩屑有效排出。
3.2 中心与边缘喷嘴直径组合优化
中心喷嘴(2号、5号喷嘴)分布在钻头中心区域,喷嘴射流从井底携岩到上返排屑槽的运移距离较长,井底清洗区域大;而钻头的边缘喷嘴(1号、3号、4号、6号喷嘴)分布在钻头主刀翼附近,清洗区域较小。因此,根据中心喷嘴和边缘喷嘴所需流量设计不同喷嘴直径组合。
喷嘴当量面积是指钻头上所有喷嘴出口面积之和,它与钻头上每个喷嘴直径有关。不同的喷嘴直径组合方案会得到不同的喷嘴当量面积。中心区域的喷嘴有两个,采用曲率ρ=1/50的弯曲喷嘴;位于边缘区域的喷嘴有4个,采用等变速形变截面喷嘴。在喷嘴当量面积相等或相近情况下,设计出4种喷嘴的组合方案(见表 3),用于钻头井底流场研究。
方案 | 中心喷嘴 数量/个 |
中心喷嘴 直径/mm |
边缘喷嘴 数量/个 |
边缘喷嘴 直径/mm |
当量直 径/mm |
1 | 2 | 10 | 4 | 10.0 | 24.5 |
2 | 2 | 11 | 4 | 9.5 | 24.6 |
3 | 2 | 12 | 4 | 8.8 | 24.4 |
4 | 2 | 13 | 4 | 8.1 | 24.5 |
保持旋切式PDC钻头其他结构参数不变,根据钻头喷嘴直径组合分析方案进行数值模拟,选取最优的喷嘴直径组合方案。
不同直径(D1)中心喷嘴轴向剖面速度云图如图 20所示。从图 20可见,当中心喷嘴(2号喷嘴)的直径由方案1中的10 mm增大到方案3的12 mm时,井底漫流层的上返高度逐渐增加;当中心喷嘴直径增大到方案4的13 mm时,受到上返区涡漩影响,漫流层的上返高度有所下降。
为了更好地描述井底速度的变化趋势,分别在井底附近选取一个位置,观察4种方案下喷嘴射流从钻头中心到井壁的速度变化,如图 21所示。从图 21可见,从方案1到方案4,随着中心喷嘴直径D1逐渐增大,射流冲击方向的井底漫流速度逐渐增大。说明适当增大喷嘴直径可以提高井底漫流速度及漫流携岩能力,降低井底岩屑颗粒的质量浓度。方案3中的速度总体变化趋势优于方案4。
不同直径(D2)边缘喷嘴轴向剖面速度云图如图 22所示。从图 22可以看出,在边缘喷嘴(1号喷嘴)直径从方案1的10.0 mm逐渐减小到方案4的8.1 mm的过程中,漫流层在排屑槽的上返高度逐渐增大,当直径为8.1 mm时上返高度最大,其次为8.8 mm的喷嘴。
井底附近流场速度的变化曲线如图 23所示。从图 23可见,从方案1到方案4,随着边缘喷嘴(1号喷嘴)直径逐渐减小,射流冲击方向的井底漫流速度逐渐增大。这说明适当减小边缘喷嘴直径可以提高井底漫流速度和漫流携岩能力,降低井底岩屑颗粒的质量浓度。方案3中的速度总体变化趋势优于方案4。
在CFD-post中,以旋切式PDC钻头的保径块上表面为基准,建立一个平行于井底的平面。以此平面沿Z轴的负方向建立一个圆柱体,用来分析钻头体表面岩屑颗粒的平均质量浓度。
如图 24所示,采用方案3时,钻头体上的岩屑颗粒平均质量浓度最低,为22.73 kg/m3。这与上述分析结果一致。
4 结论
(1) 刀盘形状的改变有利于工作流体的流动,可以降低刀盘附近岩屑颗粒的质量浓度。
(2) 采用带有弯曲喷口的喷嘴可以加强对旋切式PDC钻头刀盘区域的井底清洗效果。
(3) 当采用曲率为1/50的喷嘴时,刀盘区域的井底清洗效果最好。
(4) 采用中心喷嘴直径12 mm和边缘喷嘴直径8.8 mm的组合方案时,旋切式PDC钻头的清洗和携岩能力达到最优。
[1] |
管志川, 陈庭根, 刘希圣. PDC钻头水力结构研究进展[J]. 石油大学学报(自然科学版), 1994, 18(6): 136-142. GUAN Z C, CHEN T G, LIU X S. Review of studies of PDC bit hydraulic configurations[J]. Journal of the University of Petroleum, China, 1994, 18(6): 136-142. |
[2] |
管志川, 陈庭根, 蔡镜仑. 用矩形出口喷嘴改善PDC钻头水力结构的研究[J]. 石油大学学报(自然科学版), 1992, 16(2): 20-28. GUAN Z C, CHEN T G, CAI J L. Experimental study on improving PDC bit hydraulics by using rectangular nozzles[J]. Journal of the University of Petroleum, China, 1992, 16(2): 20-28. |
[3] |
刘刚, 陈庭根, 管志川. 喷嘴倾角对PDC钻头齿水射流清洗力影响的实验研究[J]. 石油大学学报(自然科学版), 1996, 20(4): 30-33. LIU G, CHEN T G, GUAN Z C. Effect of nozzel inclination angle on cleaning force acting on PDC bit teeth[J]. Journal of the University of Petroleum, China, 1996, 20(4): 30-33. |
[4] |
谢翠丽, 杨爱玲, 陈康民. 非对称多喷嘴平底钻头井底三维流场数值模拟[J]. 石油学报, 2002, 23(6): 77-80. XIE C L, YANG A L, CHEN K M. Numerical simulation for downhole flow field of jet on PDC bits with asymmetric multi-nozzle distribution[J]. Acta Petrolei Sinica, 2002, 23(6): 77-80. DOI:10.3321/j.issn:0253-2697.2002.06.017 |
[5] |
况雨春, 曾恒, 周学军, 等. PDC钻头水力结构优化设计研究[J]. 天然气工业, 2006, 26(4): 60-61, 67. KUANG Y C, ZENG H, ZHOU X J, et al. Study on optimized design of hydraulic structure for PDC bits[J]. Natural Gas Industry, 2006, 26(4): 60-61, 67. DOI:10.3321/j.issn:1000-0976.2006.04.020 |
[6] |
况雨春, 杨云山, 杨迎新, 等. PDC钻头井底流场CFD仿真及二次开发[J]. 钻采工艺, 2007, 30(1): 77-79. KUANG Y C, YANG Y S, YANG Y X, et al. Advanced development and smulation of bottom fluid of PDC bits based on CFD[J]. Drilling & Production Technology, 2007, 30(1): 77-79. DOI:10.3969/j.issn.1006-768X.2007.01.026 |
[7] |
CROUSE R, CHIA R. Optimization of PDC bit hydraulics by fluid simulation[R]. SPE 14221, 1985.
|
[8] |
NAZARI T, HARELAND G, AZAR J. Review of cuttings transport in directional well drilling: systematic approach[R]. SPE 132372, 2010.
|
[9] |
MOSLEMI A, AHMADI G. Study of the hydraulic performance of drill bits using a computational particle-tracking method[J]. SPE Drilling & Completion, 2014, 29(1): 28-35. |
[10] |
WATSON G R, BARTON G R, HARGRAVE G K. Using new computation fluid dynamics techniques to improve PDC bit performance[R]. SPE 37580, 1997.
|
[11] |
MORSI S A, ALEXANDER A J. An investigation of particle trajectories in two-phase flow systems[J]. Journal of Fluid Mechanics Digital Archive, 1972, 55(2): 193-208. |
[12] |
FORDER A, THEW M, HARRISON D. A numerical investigation of solid particle erosion experienced within oilfield control valves[J]. Wear, 1998, 216(2): 184-193. DOI:10.1016/S0043-1648(97)00217-2 |