摘要
泥浆举升离心泵作为无隔水管泥浆循环钻井的核心功能单元,作用是驱动储存在吸入模块中的泥浆携带岩屑沿上返管线回流至钻井平台。以往的无隔水管泥浆循环系统以圆盘泵作为举升单元,围绕圆盘泵做了具体的结构设计和理论分析,对于水力性能更好的多级离心泵研究较少。以多级离心泵作为研究对象,介绍其设计的理论依据及具体的参数,采用CFD-DEM耦合计算的方法构建离心泵的计算模型,进行了全尺度流场模拟,对离心泵的叶轮设计、流场和内颗粒分布规律、泥浆携岩运移规律进行了探讨。计算结果验证了离心泵在设计工况下具有较好的颗粒通过性。基于泥浆举升试验台,对离心泵样机进行了水力性能试验,测试结果表明该泵具有较好的工作特性,其水力性能达到设计要求,可以应用于深海无隔水管钻井作业。
使用无隔水管泥浆循环钻井(Riserless Mud Return,以下简称RMR)工艺施工时,钻井平台上的钻井泥浆泵通过钻杆将泥浆输送到钻头,钻进时产生的岩屑被泥浆携带通过环空返回到井口吸入模
国内外对深海环境的流体举升泵进行了大量的理论研究,秦如雷

图1 井底泥浆举升循环系统
Fig.1 Bottom hole mud lifting and circulation system
为了设计出具有优异水力性能和颗粒通过性能的离心式泥浆举升泵,首先需要针对离心泵研究其内部复杂流场作用下固液两相流的分布特性,特别是要深入分析固相颗粒与流体的相互作用机理,研究固液两相流输送泵中的颗粒运动规律,提高颗粒在泵中的通过性,降低泵输送固液两相流时叶轮堵塞的概
固液两相流是指在流体流动中同时存在固体颗粒的现象。在进行固液两相流的数值模拟时,可以采用多种计算模型,这些模型根据不同的应用场景和计算需求有所区别。在欧拉-欧拉模型(Eulerian-Eulerian Model)中,固体和液体都被视为连续相,各自有其速度和压力场,适用于固体颗粒浓度较高、颗粒间相互作用显著的情况;此时需要为每个相分别求解Navier-Stokes方程。在欧拉-拉格朗日模型(Eulerian-Lagrangian Model)中,液体作为连续相采用欧拉方法,而固体颗粒作为离散相采用拉格朗日方法,适用于颗粒浓度不高,颗粒可以在流体中自由运动的情况。
根据Pagalthivarth
流体连续相的运动特性符合雷诺平均Navier-Stokes方程(RANS),可选用基于欧拉坐标系的SST k-ε湍流模型用于流场数值分析,其流体动力学特征可表示为
(1) |
式中:——流体密度,kg/
离散相的求解主要是通过接触模型求解粒子碰撞过程中的力,利用牛顿第二定律计算粒子的加速度,然后更新粒子的速度和位移。粒子上的力在

图2 粒子碰撞的受力分析
Fig.2 Force analysis of particle collision
(2) |
(3) |
式中:——质量,kg;——颗粒的速度,m/s;——重力加速度,;——颗粒之间的法向接触力,;——切向接触力,;——粒子的惯性矩,kg·2;——粒子的角速度,rad/s;——粒子的半径,m;——液相施加在颗粒上的力,;——滚动摩擦扭矩,。
采用CFD-DEM耦合的固液两相流理论,对模型泵进行水-颗粒两相流数值计算,旨在掌握颗粒在泵流场内的运动和迁移规律,以便分析泵内颗粒的运动和堵塞情况,为样泵的设计研发提供技术支撑。
采用CFD-DEM双向同时耦合方法,求解连续相和离散相,不同相之间的边界满足相容性方程。采用欧拉方法计算连续相,根据连续相流场的计算结果,在悬浮于连续相中的粗颗粒上施加局部载荷。采用拉格朗日跟踪法计算离散相位,根据力平衡和动量守恒原理,通过积分粒子在轨道上的力平衡方程,可以获得粒子轨迹拉格朗日方程坐标。轨迹方程通过在独立的时间步长内逐步积分求解,然后重新计算连续相的流场。通过这种方式,反复使用迭代方法,直到连续相和离散相都达到运动平衡。最后,实现双向耦合计算。
根据前述的多级泵设计思路,结合常规钻进的环境条件,设计工况条件:举升泥浆岩屑粒径10 mm、固相体积稠度为10%、密度为1.2 g/c
(4) |
(5) |
式中:——电机轴功率,等于泵输入功率,W;——泵输出功率,W;Q——泵的排量,m³/s;——泵的扬程,m;——流体的密度,kg/m³;——重力加速度,m/
经过计算,得到,。
在此基础上,确定泵的水力参数,具体参数见
叶轮几何参数 | |
---|---|
泵进口直径/mm | 140 |
泵进口速度/(m· | 2.83 |
泵出口直径/mm | 140 |
泵出口速度/(m· | 3.77 |
额定转速n/(r·mi | 1450 |
流量Q/( | 120 |
数量 | 4 |
导叶几何参数 | |
叶轮外径D2/mm | 480 |
出口宽度b2/mm | 30 |
叶片数Z/枚 | 4 |
叶片出口角β2/(°) | 5.1~5.3 |
额定转速n/(r·mi | 1450 |
流量Q/( | 120 |
数量 | 4 |

图3 四叶泵叶轮三维模型
Fig.3 3D model of four‑blade pump wheel
该结构特征为:4片大叶片再加上4个导流叶片。设计采用的流道相对较宽,有利于颗粒的运动,降低了泵内堵塞的风险。而且叶片较厚,尤其是靠近入口的叶片,有利于提高叶片的强度,对抗颗粒碰撞和磨损的能力更强。增加小叶片的目的是确保叶轮流道的体积变化不会太大,流场更加均匀,这有利于提高颗粒的输送能力。
采用EDEM-FLUENT软件耦合计算泵中的三维非定常固液两相流。数值分析的边界条件如
湍流模型 | SST k-ω |
---|---|
入口边界 | 质量流量入口 |
出口边界 | 自由出口 |
壁面边界 | 无滑移壁面 |
壁面函数 | 可放缩壁面函数 |
接触力 | Hertz-Mindlin无滑移模型 |
流体是4 °C海水与密度为2000 kg/
考虑到计算的准确性、计算机的计算能力以及粒径,确定4叶片泵的网格总数为385315。网格独立性检查结果如

图4 泵网格独立性检查
Fig.4 Mesh independent check of pump

图5 计算域三维图和网格图
Fig.5 3D model and mesh grid of computational domain
根据泵头流线分布图(

图6 泵叶轮流线分布
Fig.6 Streamline distribution of pump impeller.

图7 泵叶轮中颗粒的速度分布
Fig.7 Velocity distribution of particles in the pump impellers
区域 | 颗粒数量 | 颗粒平均速度/(m· | 颗粒停留时间/s | 颗粒碰撞次数 | 颗粒体积浓度/% |
---|---|---|---|---|---|
区域1 | 122 | 3.13 | 0.0937 | 11 | 12.6 |
区域2 | 441 | 9.48 | 0.1396 | 397 | |
区域3 | 1094 | 3.62 | 0.3539 | 526 | |
区域4 | 96 | 2.61 | 0.3724 | 23 | 9.93 |
在数值模拟的基础上,试制了两级模型泵用于开展水力性能试验,以验证泵型设计的合理性。将模型泵与现场管线系统连接,开展了清水工况下的泵性能测试。

图8 举升试验台原理
Fig.8 Principle of lifting test bench
本次试验将对模型泵在不同介质、不同转速下的泵扬程、轴功率和效率,以及颗粒通过性进行测试。两级泵的扬程110 m。其中关键的物理量测量依据如下:
扬程和效率的测量。泵的进口压力和出口压力分别由进口压力变送器和出口压力变送器测量,其测量误差0.1%。流量由电磁流量计记录,其测量误差0.5%,实际在颗粒工况下,通过标定装置对管路系统内的颗粒量进行物理标定,确定两相流量数值。驱动泵的电机最大功率为160 kW,泵的功率将通过电测法进行测量。轴功率可通过在电机与两级泥浆泵的连接处设置扭矩仪进行测量。
在背压工况下,介质为清水,电机转速为1450 r/min,流量分别为90、100、110、120和130
对监测的泵扬程和功率、效率数据进行分析,发现该泵性能曲线较为平稳,设计点120

图9 泵的扬程-流量性能特征
Fig.9 Lift-flow performance of pump

图10 泵的轴功率-流量性能特征
Fig.10 Shaft power-flow performance of pump

图11 泵的效率-流量性能特征
Fig.11 Efficiency-flow performance of pump
(1)扬程曲线比较平坦,尤其在设计运行区间(90~120 m³/h),泵的扬程变化较小,可使水下钻井获得恒定的出口压力。
(2)轴功率曲线起点较高,关死点扬程很大,在运行区间(90~120 m³/h),泵的轴功率变化较小,因此泵的轴功率不容易突增,可保障电机运行可靠。
(3)泵的效率曲线是典型的放大设计效果,泵在大流量点的效率显著高于小流量点。设计工况120 m³时,泵的效率36.5%,90 m³/h时,泵的效率28.8%。
本文基于欧拉坐标系的SST k-ε湍流模型和CFD-DEM双向流固耦合分析,进行了用于深海泥浆举升的离心泵内部流场和岩屑颗粒运移特性的研究,验证了离心泵叶轮结构设计方案及参数的合理性和可行性。
应用水力提升模拟试验台进行了离心式泥浆举升泵样机的水力特性测试,通过对试验数据和原始数据的对比分析,清水工况下,泵的效率曲线是典型的放大设计效果,大流量点泵的效率显著高于小流量点。验证了该泵具有更好的工作特性,其水力性能达到了设计的要求。
对照深海RMR系统的设计标准,离心式泥浆举升泵的数值模拟和样机试验结果均表明所设计的离心泵可应用于深海无隔水管钻井作业中,并且该泵的设计方法和试验过程是科学合理的。
参考文献(References)
陈浩文,于彦江,王艳丽,等.气举技术应用于深海无隔水管泥浆回收钻井工艺可行性分析[J].钻探工程,2022,49(2):9-15. [百度学术]
CHEN Haowen, YU Yanjiang, WANG Yanli, et al. Feasibility analysis of gas lift technology for application in deep‑sea riserless mud recovery drilling[J]. Drilling Engineering, 2022,49(2):9-15. [百度学术]
陈浩文,刘晓林,王林清,等.无隔水管泥浆回收钻井技术控制系统功能设计[J].钻探工程,2021,48(S1):375-380. [百度学术]
CHEN Haowen, LIU Xiaolin, WANG Linqing, et al. Functional design of the control system of riserless mud recovery drilling technology[J]. Drilling Engineering, 2021,48(S1):375-380. [百度学术]
高本金,陈国明,殷志明,等.深水无隔水管钻井液回收钻井技术[J].石油钻采工艺,2009,31(2):44-47. [百度学术]
GAO Benjin, CHEN Guoming, YIN Zhiming, et al. Deepwater riserless mud recovery drilling technology[J]. Oil Drilling Production Technology, 2009,31(2):44-47. [百度学术]
Cohen J H, Kleppe J, Grønås T, et al. Gulf of Mexico’s first application of riserless mud recovery for top‑hole drilling—A case study[C]//Offshore Technology Conference. Houston, Texas, USA, 2010: OTC-20939-MS.. [百度学术]
王偲,谢文卫,张伟,等.RMR技术在海域天然气水合物钻探中的适应性分析[J].探矿工程(岩土钻掘工程),2020,47(2):17-23. [百度学术]
WANG Cai, XIE Wenwei, ZHANG Wei, et al. Adaptability of RMR for Marine gas hydrate drilling[J]. Exploration Engineering (Rock Soil Drilling and Tunneling), 2020,47(2):17-23. [百度学术]
Aird P. Deepwater well planning and Design-Science Direct[M]//Deepwater Drilling, 2019:257-283. [百度学术]
王江帅,李军,柳贡慧,等.基于井下分离的深水双梯度钻井参数优化[J].石油勘探与开发,2019,46(4):776-781. [百度学术]
WANG Jiangshuai, LI Jun, LIU Gonghui, et al. Parameters optimization in deepwater dual‑gradient drilling based on downhole separation[J]. Petroleum Exploration and Development, 2019,46(4):776-781. [百度学术]
于洋,钱则刚,张和平,等.多列往复式压缩机轴系扭振特性研究[J].流体机械,2013,41(7):14-18. [百度学术]
YU Yang, QIAN Zegang, ZHANG Heping, et al. Research on the crankshaft torsional vibration of multi‑row reciprocating compressor[J]. Fluid Machinery, 2013,41(7):14-18. [百度学术]
唐斌,刘广彬,许海平,等.大型工艺往复压缩机系统振动分析[J].中国机械工程,2014,25(7):873-877. [百度学术]
TANG Bin, LIU Guangbin, XU Haiping, et al. Vibration analysis of large scale process reciprocating compressor system[J]. China Mechanical Engineering, 2014,25(7):873-877. [百度学术]
秦如雷,于彦江,陈浩文,等.无隔水管泥浆循环举升泵选型及性能参数计算方法[J].钻探工程,2021,48(S1):381-385. [百度学术]
QIN Rulei, YU Yanjiang, CHEN Haowen, et al. Selection and calculation of parameters of the mud lifting pump for RMR[J]. Drilling Engineering, 2021,48(S1):381-385. [百度学术]
解永超,陈国明,张叶,等.叶片圆盘泵压力脉动及叶轮径向力特性分析[J].石油机械,2015(10):63-67. [百度学术]
XIE Yongchao, CHEN Guoming, ZHANG Ye, et al. Numerical analysis on pressure fluctuation and radial force of disc pump with straight blade[J]. China Petroleum Machinery, 2015(10):63-67. [百度学术]
陈永超,陈国明,周昌静,等.圆盘泵叶片结构形式对泵外特性影响分析[J].机械设计,2014,31(9):93-96. [百度学术]
CHEN Yongchao, CHEN Guoming, ZHOU Changjing, et al. Analysis on the influence of blade structure on exterior characteristics of disc pump[J]. Journal of Machine Design, 2014,31(9):93-96. [百度学术]
邹伟生,刘瑞仙,刘少军.粗颗粒海底矿石浆体提升电泵研究[J].中国机械工程,2019,30(24):2939-2944. [百度学术]
ZOU Weisheng, LIU Ruixian, LIU Shaojun. Study on lifting motor pumps for coarse particle slurry in sea bed mining[J]. China Mechanical Engineering, 2019,30(24):2939-2944. [百度学术]
邹伟生,卢勇,李哲奂.深海采矿提升泵的数值模拟分析[J].湖南大学学报(自然科学版),2013,40(6):59-63. [百度学术]
ZOU Weisheng, LU Yong, LI Zhehuan. Numerical simulation and analyses of lift pump in deep sea mining[J]. Journal of Hunan University(Natural Sciences), 2013,40(6):59-63. [百度学术]
Li X, Zhang J, Tang X, et al. Study on wellbore temperature of riserless mud recovery system by CFD approach and numerical calculation[J]. Petroleum, 2020,6(2):163-169. [百度学术]
黄思,杨富翔,宿向辉.运用CFD-DEM耦合模拟计算离心泵内非稳态固液两相流动[J].科技导报,2014,32(27):28-31. [百度学术]
HUANG Si, YANG Fuxiang, SU Xianghui. Unsteady numerical simulation for solid-liquid two‑phase flow in centrifugal pump by CFD-DEM coupling[J]. Science Technology Review, 2014,32(27):28-31. [百度学术]
Zhou M M, Wang S, Kuang S B, et al. CFD-DEM modelling of hydraulic conveying of solid particles in a vertical pipe[J]. Powder Technology, 2019,354:893-905. [百度学术]
夏铖,赵睿杰,施卫东,等.基于DEM-CFD耦合的泵内颗粒流动特性研究[J].工程热物理学报,2021,42(2):357-369. [百度学术]
XIA Cheng, ZHAO Ruijie, SHI Weidong, et al. Numerical investigation of particle induced erosion in a mixed pump by CFD-DEM coupled method[J]. Journal of Engineering Thermophysics, 2021,42(2):357-369. [百度学术]
汪家琼,蒋万明,孔繁余,等.固液两相流离心泵内部流场数值模拟与磨损特性[J].农业机械学报,2013,44(11):53-60. [百度学术]
WANG Jiaqiong, JIANG Wanming, KONG Fanyu, et al. Numerical simulation of solid-liquid two‑phase turbulent flow and wear characteristics of centrifugal pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013,44(11):53-60. [百度学术]
Zhu Z C, Xie P, Ou G F, et al. Design and experimental analyses of small‑flow high‑head centrifugal‑vortex pump for gas-liquid two‑phase mixture[J]. Chinese Journal of Chemical Engineering, 2008,16(4):528-534. [百度学术]
Ariely S, Khentov A. Erosion corrosion of pump impeller of cyclic cooling water system[J]. Engineering Failure Analysis, 2006,13(6):925-932. [百度学术]
Llewellyn R J, Yick S K, Dolman K F. Scouring erosion resistance of metallic materials used in slurry pump service[J]. Wear, 2004,256(6):592-599. [百度学术]
Morsi S A, Alexander A J. An investigation of particle trajectories in two‑phase flow systems[J]. Journal of Fluid Mechanics, 1972,55(2):193-208. [百度学术]
李昳.离心泵内部固液两相流动数值模拟与磨损特性研究[D].杭州:浙江理工大学,2014. [百度学术]
LI Die. The research on numerical simulation and abrasion property of solid-liquid two‑phase‑flow centrifugal pump[D]. Hangzhou: Zhejiang Sci‑tech University, 2014. [百度学术]
Zhang Z, Sun B J, Wang Z Y, et al. Intelligent well killing control method driven by coupling multiphase flow simulation and real‑time data[J]. Journal of Petroleum Science and Engineering, 2022,213:110337. [百度学术]
Pagalthivarthi K V, Gupta p K, Tyagi V, et al. CFD prediction of erosion wear in centrifugal slurry pumps for dilute slurry flows[J]. The Journal of Computational Multiphase Flows, 2011,3(4):225-245. [百度学术]