Numerical Simulation of the Effects of Eccentric Rotation of the Drill String on Annular Frictional Pressure Drop
-
摘要:
准确预测钻柱偏心旋转工况下的环空摩阻压降是复杂结构井控压钻井的重要理论基础，但常规钻井液环空摩阻压降计算方法无法直接计算复杂结构井的环空摩阻压降。为此，应用数值模拟方法，分析了偏心度（0～67.42%）和钻柱转速(0～114.65 r/min)对典型环空（ϕ127.0 mm钻杆和ϕ215.9 mm井眼）中摩阻压降梯度的影响。分析结果表明：偏心度小于45.00%时，转速和偏心度对摩阻压降梯度影响较弱，摩阻压降梯度随转速增大略有降低，随偏心度增大而增大；偏心度大于45.00%时，低转速（<60 r/min）下摩阻压降梯度随偏心度增大而降低，高转速（≥60 r/min）下摩阻压降梯度随偏心度增大而略有增大。基于数值模拟结果，建立了偏心度分类的无因次偏心环空摩阻压降梯度预测模型，计算了南海某水平井ϕ215.9 mm井段的ECD，并与PWD测试结果进行了对比，平均相对误差为0.45%，表明该模型具有较好的准确性。研究结果表明，无因次偏心旋转环空摩阻压降计算模型可以精细描述环空压力场和准确计算ECD，为控压钻井水力参数优化提供指导。
Abstract:Accurate prediction of the annular frictional pressure drop under eccentric rotation of drill string is an important theoretical basis for managed pressure drilling (MPD) in complex structure wells. However, the conventional calculation methods for the annular frictional pressure drop of drilling fluid cannot be directly applied to calculating the annular frictional pressure drop in complex structure wells. For this reason, the influences of eccentricity (0−67.42%) and the rotational speed (0−114.65 r/min)of the drill string on the frictional pressure drop gradient in a typical annulus (created by a ϕ127.0 mm drill pipe and a ϕ215.9 mm wellbore) were analyzed. The results show that when the eccentricity is lower than 45.00%, the rotational speed and eccentricity have a weak influence on the frictional pressure drop gradient. Specifically, the frictional pressure drop gradient decreases slightly with the increase of rotational speed but increases with the increase of eccentricity; when the eccentricity is higher than 45.00%, the frictional pressure drop gradient decreases with the increase of eccentricity at low rotational speed（<60 r/min）, and the friction pressure drop gradient increases slightly with the increase of eccentricity at high rotational speed (>60 r/min). According to the numerical simulation results, a dimensionless frictional pressure drop gradient prediction model with eccentricity classification was built. The equivalent circulating density (ECD) of ϕ215.9 mm section of a horizontal well in the South China Sea was calculated by the proposed model. The results were then compared with the pressure-while-drilling (PWD) test results, with an average relative error of 0.45%, indicating that the proposed model has favorable accuracy. This study concludes that the proposed calculation model of the dimensionless annular frictional pressure drop under eccentric rotation can precisely describe the annular pressure field and ECD, and provide guidance for the hydraulic parameter optimization of MPD.
-
斜井、大位移井和水平井等复杂结构井钻井过程中，钻柱在井眼中会出现偏心旋转，影响了环空钻井液的流动和循环压耗，进而使井底钻井液当量循环密度（ECD）发生变化，而准确预测环空摩阻压降是复杂结构井控压钻井的重要理论基础[1-2]。钻井液环空摩阻压降常规计算方法假设钻柱在井眼内居中，且未考虑钻柱旋转，因而无法直接计算复杂结构井的环空摩阻压降。为此，国内外学者开展了偏心环空摩阻压降计算理论和试验研究[3-25]，建立了各种偏心环空流动摩阻压降计算模型，但均存在一定局限性，如解析模型中只考虑了偏心，没有考虑钻柱旋转，尚未建立符合现场工况的偏心环空摩阻压降计算模型。此外，环空层流摩阻压降的数值本身并不高，且偏心和旋转对环空摩阻压降的影响一般小于20%，因此，试验研究对压力计精度和操作精度的要求更加苛刻，导致文献公布的试验数据不系统，规律认识不一致。所以，基于试验数据建立的回归模型适用性差，现场应用时存在计算结果不合理现象。为此，笔者应用数值模拟方法，研究了正常钻井参数条件下，转速和偏心对复杂结构井典型环空（ϕ127.0 mm钻杆和ϕ215.9 mm井眼）摩阻压降的影响，进而建立了偏心旋转环空无因次摩阻压降计算模型，可以精细描述环空压力场和准确计算ECD，为控压钻井水力参数优化提供指导。
1. 三维偏心环空流动数值模拟
1.1 数值模拟参数
复杂结构井的生产井段井径大多为215.9 mm，通常为控压钻井的主要井段。该井段钻进时，钻柱转速一般不超过120 r/min，钻井液排量不超过40 L/s。ϕ127.0 mm钻杆一般应用ϕ168.0 mm接头，考虑井眼10%的扩径率，钻杆接头紧贴井壁情况下的钻柱偏心度为62.89%，偏心度计算公式可表示为：
E=δrw−rd (1) 式中：E为偏心度；
δ 为偏心距，即沿井筒轴线视图中钻杆和井筒2个圆心间的距离，m；rw 为井筒半径，m；rd 为钻杆半径，m。数值模拟时，假设钻井液性能和排量不变，改变钻柱转速和偏心度，主要模拟参数设定为：钻井液流性指数n为0.615，稠度系数K为0.518 Pa·sn，钻井液密度1 300 kg/m3、入口速度1.50 m/s，钻柱偏心度0～67.42%，转速0～114.65 r/min。由于偏心度和转速模拟矩阵划分得越细密，建立的回归模型预测结果精度越高，但模型也越复杂。综合考虑这2方面因素和偏心度对摩阻压降梯度的影响规律，不均匀设置了12个偏心度：0（偏心距0），22.47%（偏心距10 mm），26.97%（偏心距12 mm），31.46%（偏心距14 mm），35.96%（偏心距16 mm），40.45%（偏心距18 mm），44.94%（偏心距20 mm），49.44%（偏心距22 mm），53.93%（偏心距24 mm），58.43%（偏心距26 mm），62.92%（偏心距28 mm）和67.42%（偏心距30 mm）；均匀设置了7个转速：0，19.11，38.22，57.32，76.43，95.54和114.65 r/min，构建了84个数值模拟点。
1.2 物理模型及网格划分
为消除入口段对模拟结果的影响，设定环空长度为5 m（大于20倍环空直径），重力方向为X轴，环空轴线方向为Z轴，应用Design Modeler建立了12个偏心度的环空流道几何模型，并进行了结构化网格划分以加速收敛。以22.47%偏心度模型为例，网格划分结果如图1所示，共划分了640 000个控制体网格，673 680个节点。网格最大扭曲度为0.28，最小扭曲度为1.54×10−2，平均扭曲度为0.13，网格质量优良，满足FLUENT对网格质量的要求，有利于计算收敛。
![]()
图 1 22.47%偏心度下偏心环空物理模型与网格划分Figure 1. Physical model and meshing of eccentric annulus with the eccentricity of 22.47%1.3 数值模拟方法有效性分析及参数设置
应用FLUENT进行三维偏心环空流动数值模拟分析时，需要确定的主要问题有：流态模型（层流或湍流模型）、离散格式、算法、非牛顿流体最大黏度及边界条件。笔者基于可以获取准确解析解的同心环空不旋转工况，确定上述参数。流速为1.50 m/s时，应用摩阻压降解析计算方法，得到雷诺数为2 300，层流摩阻压降计算结果（894.86 Pa/m）略大于紊流摩阻压降计算结果（853.00 Pa/m），故数值模拟研究时选择层流模型。为提高计算精度，选择二阶迎风格式，应用SIMPLE算法加速收敛。入口边界设定为速度入口（1.50 m/s），出口边界为压力出口，钻杆壁面和环空壁面为固体静止壁面，壁面处无滑移；在模拟钻柱旋转时，钻杆壁面处设定为对应的旋转角速度。调整钻井液黏度为最大值0.055 Pa·s时，数值模拟结果为896.91 Pa/m，与解析结果（894.86 Pa/m）相对误差为0.23%，故应用上述条件进行偏心旋转工况下的数值模拟分析。
2. 数值模拟结果分析
2.1 偏心环空流动特征分析
以偏心度53.93%和转速114.65 r/min的情况为例，在模型轴线4 m处建立一个截面，该界面处轴向、切向速度剖面如图2、图3所示。为进一步分析偏心和旋转对流速的影响，在该截面处建立一条自窄间隙到宽间隙的高边方向线，获取并绘制偏心度53.93%、转速114.65 r/min，偏心度53.93%、转速19.11 r/min，偏心度22.47%、转速114.65 r/min和偏心度22.47%、转速19.11 r/min等4种工况下的轴向和切向流速沿该线的分布，结果见图4和图5。为方便分析，定义X轴正方向为井筒高边方向(自钻杆下方窄环空间隙指向钻杆上方的宽环空间隙)，Y轴为截面中垂直于X轴的方向。
![]()
图 4 轴向流速沿高边方向线分布Figure 4. Axial velocity distribution along the high-sidedirection line![]()
图 5 切向流速沿高边方向线分布Figure 5. Tangential velocity distribution along the high-side direction line从图2可以看出：1）钻杆上方环空中间区带的轴向流速最高，钻杆下方窄环空区带的轴向流速最低；2）钻杆顺时针旋转时，Y轴右侧环空的轴向流速高于左侧环空；3）钻杆壁面和井壁处的轴向流速为0，与数值模拟中2个壁面处的无滑移边界条件一致。
从图3可以看出：1）距离钻杆壁面越近，钻井液的切向流速越大；2）钻杆壁面下部区域的钻井液切向流速最大，上部宽环空区带的切向流速最小。
从图4可以看出：1）钻杆上方宽环空区带和钻杆下方窄环空区带的钻井液轴向流速均基本呈抛物线分布，壁面处轴线流速为0，中间位置轴向流速最大；2）偏心度越大，宽环空区带钻井液轴向流速越大，窄环空区带钻井液轴向流速越小；3）偏心度越大，转速对窄环空区带的钻井液轴向流速分布影响越大；4）相同偏心度条件下，不同转速下的宽环空区带钻井液轴向流速分布基本重合，即转速对宽环空区带的钻井液轴向流速影响极小。
从图5可以看出：1）窄环空区带的钻井液切向流速大于宽环空区带；2）偏心度越大，宽环空区带钻井液切向流速越低，窄环空区带钻井液切向流速越大，这2个区域钻井液切向流速差异也越大；3）转速比偏心度对钻井液切向流速的影响更大。
钻柱偏心和旋转对环空钻井液流速分布的影响复杂，难以直接得到二者对摩阻压降的影响规律。但从岩屑运移的角度来说，钻柱偏心度越大，转速对窄环空区带钻井液轴向流速和切向流速的影响越大，而钻杆下方恰恰是岩屑堆积的重点区域，因此可以推断，提高转速可显著改善大斜度井段岩屑的运移效果。
2.2 摩阻压降梯度数值模拟计算结果
流体在流道内的流动通常分为流动发展阶段和稳定阶段，为了准确获取摩阻压降梯度数值模拟结果，沿模型中钻柱上方流域中心点建立一条轴线，分析压力沿该轴线的变化情况，以确定数值模拟的流动稳定阶段。钻井液静压力沿轴线变化的模拟结果见图6。
![]()
图 6 钻井液静压力沿宽流域中心线的分布Figure 6. Static drilling fluid pressure distribution along the center line of wide flow domain从图6可以看出，流道入口附近压降梯度较大，距流道入口1.50 m之后，钻井液静压力与Z轴位置呈理想线性关系，这表明流道入口1.50 m内为流动发展阶段，之后进入流动稳定阶段。基于此，在模型轴线2.00和5.00 m处建立2个截面，分析2个截面之间的压差，求取摩阻压降梯度。不同偏心度和钻杆转速下的摩阻压降梯度计算结果见表1。
表 1 偏心环空摩阻压降梯度数值模拟计算结果Table 1. Numerical computation results of frictional pressure drop gradient in eccentric annulus偏心度，% 不同钻柱转速下的摩阻压降梯度/（Pa·m−1） 0 r/min 19.11 r/min 38.22 r/min 57.32 r/min 76.43 r/min 95.54 r/min 114.65 r/min 0 898.91 898.69 898.73 898.80 898.91 899.03 899.18 22.47 921.99 919.85 916.67 915.89 915.74 914.72 912.90 26.97 930.93 927.63 922.68 921.31 921.12 919.76 917.43 31.46 939.65 934.85 927.73 925.76 925.82 924.38 921.85 35.96 946.11 939.65 930.41 928.23 929.19 928.37 926.25 40.45 947.04 939.42 929.21 927.89 931.11 932.09 931.20 44.94 938.21 931.38 923.14 924.87 932.38 936.60 937.83 49.44 916.31 914.30 912.82 920.49 934.65 943.11 946.98 53.93 882.02 890.43 901.17 917.22 939.97 952.41 958.70 58.43 841.34 865.33 892.13 916.90 948.91 964.11 971.14 62.92 800.91 842.91 886.76 917.93 957.17 974.83 977.08 67.42 763.60 820.64 879.73 913.85 956.12 977.72 968.23 2.3 偏心度对偏心环空摩阻压降梯度的影响
基于表1数据，绘制不同钻杆转速下环空摩阻压降梯度与偏心度的关系曲线，如图7所示（图中红色加粗直线为根据解析式计算出的同心不旋转摩阻压降梯度）。从图7可以看出：1）不同转速下，偏心度对环空摩阻压降梯度的影响规律不同，且没有单调关系；2）大多数情况下，偏心旋转环空摩阻压降梯度大于同心不旋转环空摩阻压降梯度，但在偏心度大于53.93%且转速小于38 r/min时，偏心旋转环空摩阻压降梯度低于同心不旋转环空摩阻压降梯度，也就是说，大斜度井段若采取滑动钻进方式且不考虑岩屑堆积造成的流动截面积减小，环空摩阻压降梯度会低于旋转钻进工况；3）不同转速下，偏心环空摩阻压降梯度与偏心度的关系曲线在偏心度约45%处存在一个交点，偏心度小于45%时，随偏心度增大，偏心环空摩阻压降梯度均略有增大；4）偏心度大于45%且转速高于57 r/min时，偏心环空摩阻压降梯度随偏心度增大而增大；5）偏心度大于45%且转速小于57 r/min时，偏心环空摩阻压降梯度随偏心度增大而减小。
![]()
图 7 偏心环空摩阻压降梯度与钻柱偏心度的关系曲线Figure 7. Variation of frictional pressure drop gradient ineccentric annulus with eccentricity of drill string2.4 转速对偏心环空摩阻压降梯度的影响
基于表1数据，绘制不同钻杆偏心度下环空摩阻压降梯度与转速的关系曲线（见图8）。从图8可以看出：1）钻杆偏心度小于53.93%时，转速对环空摩阻压降梯度的影响较小，随着转速增大，偏心环空摩阻压降梯度略有降低；2）钻杆偏心度大于53.93%时，偏心环空摩阻压降梯度随转速增大而增大，但转速大于60 r/min后，转速对偏心环空摩阻压降梯度的影响减弱，也就是说，低转速下偏心度对摩阻压降梯度的影响远大于高转速情况。
![]()
图 8 偏心环空摩阻压降梯度与钻柱转速的关系曲线Figure 8. Variation frictional pressure drop gradient in eccentric annulus with rotational speed of drill string3. 偏心环空摩阻压降梯度预测模型的建立
将表1中的摩阻压降梯度数据除以同心不旋转工况下的摩阻压降梯度，得到无因次偏心环空摩阻压降梯度（见表2）。由于转速和偏心度对环空摩阻压降梯度的影响规律复杂，不具有单调特征，若强行建立一个综合考虑转速和偏心度的二元关系模型R=f（N，E），则预测精度会比较低。为此，按照偏心度和转速进行分类，分别建立偏心环空摩阻压降梯度预测模型。
表 2 无因次偏心环空摩阻压降梯度数值模拟结果Table 2. Numerical simulation results of dimensionless eccentric frictional pressure drop gradient偏心度，% 不同钻柱转速下的无因次偏心环空摩阻压降梯度 0 r/min 19.11 r/min 38.22 r/min 57.32 r/min 76.43 r/min 95.54 r/min 114.65 r/min 0 1.0000 0.9998 0.9998 0.9999 1.0000 1.0001 1.0003 22.47 1.0257 1.0233 1.0198 1.0189 1.0187 1.0176 1.0156 26.97 1.0356 1.0320 1.0264 1.0249 1.0247 1.0232 1.0206 31.46 1.0453 1.0400 1.0321 1.0299 1.0299 1.0283 1.0255 35.96 1.0525 1.0453 1.0351 1.0326 1.0337 1.0328 1.0304 40.45 1.0535 1.0451 1.0337 1.0322 1.0358 1.0369 1.0359 44.94 1.0437 1.0361 1.0270 1.0289 1.0372 1.0419 1.0433 49.44 1.0194 1.0171 1.0155 1.0240 1.0398 1.0492 1.0535 53.93 0.9812 0.9906 1.0025 1.0204 1.0457 1.0595 1.0665 58.43 0.9360 0.9626 0.9925 1.0200 1.0556 1.0725 1.0804 62.92 0.8910 0.9377 0.9865 1.0212 1.0648 1.0845 1.0870 67.42 0.8495 0.9129 0.9787 1.0166 1.0636 1.0877 1.0771 建立心环空摩阻压降梯度预测模型时，定义了无因次旋转因子
RN 和无因次偏心因子RE ：RN=Δpf,NΔpf,0 (2) RE=Δpf,EΔpf,0 (3) 式中：RN为无因次旋转因子，即相同偏心度下某转速时的摩阻压降梯度与基准摩阻压降梯度（即转速为0、偏心度为0时的摩阻压降梯度）的比值；RE为无因次偏心因子，即相同转速下某偏心度时的摩阻压降梯度与基准摩阻压降梯度的比值；
Δpf,0 为同心不旋转工况下的幂律流体摩阻压降梯度，Pa/m；Δpf,N 为相同偏心度不同转速时的摩阻压降梯度，Pa/m；Δpf,E 相同转速不同偏心度时的摩阻压降梯度，Pa/m。将2类偏心环空摩阻压降梯度预测模型的预测结果与数值模拟结果进行比较发现，偏心度分类预测模型的残差平方和为3.61×10−4，而转速分类预测模型的残差平方和为2.25×10−3，前者的预测精度比后者高一个数量级，所以最终确定将表2中的数据进行回归，建立偏心度分类环空摩阻压降梯度预测模型。
1）E=0时的预测模型为：
RN=7.254×10−8N2−3.991×10−6N+0.9999 (4) 2）E=22.47%时的预测模型为：
RN=3.965×10−7N2−2.369×10−4N+1.0250 (5) 3）E=26.97%时的预测模型为：
RN=8.114×10−7N2−1.247×10−4N+1.0350 (6) 4）E=31.46%时的预测模型为：
RN=1.588×10−6N2−2.405×10−4N+1.0450 (7) 5）E=35.96%时的预测模型为：
RN=2.539×10−6N2−3.645×10−4N+1.0520 (8) 6）E=40.45%时的预测模型为：
RN=3.579×10−6N2−4.356×10−4N+1.0530 (9) 7）E=44.94%时的预测模型为：
RN=−6.572×10−8N3+1.544×10−5N2−9.16×10−4N+1.0450 (10) 8）E=49.44%时的预测模型为：
RN=−8.842×10−8N3+1.855×10−5N2−5.72×10−4N+1.0210 (11) 9）E=53.93%时的预测模型为：
RN=−1.067×10−7N3+1.875×10−5N2−9.066×10−6N+0.9822 (12) 10）E=58.43%时的预测模型为：
RN=−1.141×10−7N3+1.496×10−5N2+1.039×10−3N+0.9366 (13) 11）E=62.92%时的预测模型为：
RN=−1.158×10−7N3+8.551×10−6N2+2.252×10−3N+0.8916 (14) 12）E=67.42%时的预测模型为：
RN=−1.277×10−7N3+3.688×10−6N2+3.251×10−3N+0.8503 (15) 式中：N为转速，r/min。
利用上述预测模型计算偏心环空摩阻压降梯度的主要步骤为：
1）基于钻井液流变性和流道尺寸，根据解析式[26]计算同心不旋转工况下的幂律流体摩阻压降梯度
Δpf,0 ；2）根据计算点所在井段特征（增斜段、降斜段还是稳斜段），以及井斜角、造斜率、钻柱尺寸、井眼尺寸和稳定器尺寸等参数，依据钻柱屈曲理论[17,27-28]，确定偏心度E；
3）查找偏心度E在12个偏心度的区间[E1，E2]，选取相应的预测模型计算偏心度E1和E2时的无因次旋转因子
RN,1 和RN,2 ，应用插值法计算偏心度E时的无因次旋转因子RN ；4）根据式（2）计算偏心环空摩阻压降梯度
Δpf,N 。4. 模型预测结果验证
南海某水平井水平段的钻具组合为：ϕ215.9 mm钻头+ ϕ171.0 mm旋转导向工具+ϕ206.0 mm倒划眼稳定器+ϕ165.1 mm液压震击器+ϕ127.0 mm加重钻杆8根+ϕ127.0 mm钻杆。钻井液密度为1 280 kg/m3，稠度系数为0.72 Pa·sn，流性指数为0.66；钻井液排量为0.033 m3/s，钻柱转速为110 r/min。
基于上述数据，首先，应用Haciislamoglu方法[17]计算了偏心度为20%、30%、40%和50%时文中提出的无因次摩阻压降梯度系数
Rlam （对应文中的RN），分别为0.94、0.98、0.98和0.98，均小于1.0，与现场工况不符，进一步印证了基于有限室内试验数据建立的回归模型的实用性不强。之后，应用建立的偏心环空摩阻压降梯度预测模型，以10 m为距离步长，计算了该井钻井液循环工况下ϕ215.9 mm井段偏心环空摩阻压降。为方便与现场PWD实测数据进行对比，按照已有公式计算得到了该井段的当量循环密度，计算值与实测值对比见图9，平均相对误差为0.45%，验证了模型的准确性。![]()
图 9 南海某水平井ϕ215.9 mm井段当量循环密度计算值与实测值对比Figure 9. Comparison between calculated and measured values of ECD in ϕ215.9 mm section of a horizontal well in the South China Sea钻井液循环工况下某井深处当量循环密度的计算公式为[29]：
ρce=ph+pfDv (16) 其中pf=N∑i=1ΔLiΔpf,i (17) 式中：ρce为当量循环密度，kg/m3；
ph 为静液柱压力，Pa；pf 为摩阻压降，根据井眼轨迹分段计算并累计求和，Pa；Dv 为目标点垂深，m。5. 结论与认识
1）三维环空流动数值模拟结果显示：偏心度越大，转速对窄环空区带的钻井液轴向流速分布影响越大，但转速对宽环空区带的钻井液轴向流速影响极小；窄环空区带的钻井液切向流速大于宽环空区带，偏心度越大，宽环空区带钻井液切向流速越低，窄环空区带的钻井液切向流速越大，这2个区域的钻井液切向流速差异也越大，转速比偏心度对钻井液切向流速的影响更大。
2）模拟计算结果显示，偏心环空摩阻压降梯度与偏心度、钻柱转速之间的关系复杂，不存在单调关系，难以建立综合考虑转速和偏心度的具有较高预测精度的二元关系模型；而按照偏心度和转速进行分类，分别建立偏心环空摩阻压降梯度预测模型，更有利于工程应用。
3）对比预测模型的预测结果与数值模拟结果发现，偏心度分类预测模型的预测精度比转速分类预测模型的预测精度高一个数量级，为此，通过对数值模拟数据回归分析，建立了偏心度分类环空摩阻压降梯度预测模型。
4）应用偏心度分类环空摩阻压降梯度预测模型，计算了南海某水平井的ϕ215.9 mm井段当量循环密度，并与PWD测试结果进行了对比，平均相对误差为0.45%，验证了模型的准确性，可用于精细描述钻井环空循环压力场和准确计算当量循环密度。
-
![]()
图 1 22.47%偏心度下偏心环空物理模型与网格划分
Figure 1. Physical model and meshing of eccentric annulus with the eccentricity of 22.47%
![]()
图 4 轴向流速沿高边方向线分布
Figure 4. Axial velocity distribution along the high-sidedirection line
![]()
图 5 切向流速沿高边方向线分布
Figure 5. Tangential velocity distribution along the high-side direction line
![]()
图 6 钻井液静压力沿宽流域中心线的分布
Figure 6. Static drilling fluid pressure distribution along the center line of wide flow domain
![]()
图 7 偏心环空摩阻压降梯度与钻柱偏心度的关系曲线
Figure 7. Variation of frictional pressure drop gradient ineccentric annulus with eccentricity of drill string
![]()
图 8 偏心环空摩阻压降梯度与钻柱转速的关系曲线
Figure 8. Variation frictional pressure drop gradient in eccentric annulus with rotational speed of drill string
![]()
图 9 南海某水平井ϕ215.9 mm井段当量循环密度计算值与实测值对比
Figure 9. Comparison between calculated and measured values of ECD in ϕ215.9 mm section of a horizontal well in the South China Sea
表 1 偏心环空摩阻压降梯度数值模拟计算结果
Table 1 Numerical computation results of frictional pressure drop gradient in eccentric annulus
偏心度，% 不同钻柱转速下的摩阻压降梯度/（Pa·m−1） 0 r/min 19.11 r/min 38.22 r/min 57.32 r/min 76.43 r/min 95.54 r/min 114.65 r/min 0 898.91 898.69 898.73 898.80 898.91 899.03 899.18 22.47 921.99 919.85 916.67 915.89 915.74 914.72 912.90 26.97 930.93 927.63 922.68 921.31 921.12 919.76 917.43 31.46 939.65 934.85 927.73 925.76 925.82 924.38 921.85 35.96 946.11 939.65 930.41 928.23 929.19 928.37 926.25 40.45 947.04 939.42 929.21 927.89 931.11 932.09 931.20 44.94 938.21 931.38 923.14 924.87 932.38 936.60 937.83 49.44 916.31 914.30 912.82 920.49 934.65 943.11 946.98 53.93 882.02 890.43 901.17 917.22 939.97 952.41 958.70 58.43 841.34 865.33 892.13 916.90 948.91 964.11 971.14 62.92 800.91 842.91 886.76 917.93 957.17 974.83 977.08 67.42 763.60 820.64 879.73 913.85 956.12 977.72 968.23 
下载: 导出CSV
表 2 无因次偏心环空摩阻压降梯度数值模拟结果
Table 2 Numerical simulation results of dimensionless eccentric frictional pressure drop gradient
偏心度，% 不同钻柱转速下的无因次偏心环空摩阻压降梯度 0 r/min 19.11 r/min 38.22 r/min 57.32 r/min 76.43 r/min 95.54 r/min 114.65 r/min 0 1.0000 0.9998 0.9998 0.9999 1.0000 1.0001 1.0003 22.47 1.0257 1.0233 1.0198 1.0189 1.0187 1.0176 1.0156 26.97 1.0356 1.0320 1.0264 1.0249 1.0247 1.0232 1.0206 31.46 1.0453 1.0400 1.0321 1.0299 1.0299 1.0283 1.0255 35.96 1.0525 1.0453 1.0351 1.0326 1.0337 1.0328 1.0304 40.45 1.0535 1.0451 1.0337 1.0322 1.0358 1.0369 1.0359 44.94 1.0437 1.0361 1.0270 1.0289 1.0372 1.0419 1.0433 49.44 1.0194 1.0171 1.0155 1.0240 1.0398 1.0492 1.0535 53.93 0.9812 0.9906 1.0025 1.0204 1.0457 1.0595 1.0665 58.43 0.9360 0.9626 0.9925 1.0200 1.0556 1.0725 1.0804 62.92 0.8910 0.9377 0.9865 1.0212 1.0648 1.0845 1.0870 67.42 0.8495 0.9129 0.9787 1.0166 1.0636 1.0877 1.0771
下载: 导出CSV
-
[1] 朱丽华,向雪琳,邓玉涵,等. 控制压力钻井设备概述[J]. 钻采工艺,2010,33(3):43–47. ZHU Lihua, XIANG Xuelin, DENG Yuhan, et al. Summarization of managed pressure drilling equipment[J]. Drilling & Production Technology, 2010, 33(3): 43–47.
[2] 刘伟,周英操,石希天,等. 塔里木油田库车山前超高压盐水层精细控压钻井技术[J]. 石油钻探技术,2020,48(2):23–28. doi: 10.11911/syztjs.2020034 LIU Wei, ZHOU Yingcao, SHI Xitian, et al. Precise managed pressure drilling technology for ultra-high pressure brine layer in the Kuqa Piedmont of the Tarim Oilfield[J]. Petroleum Drilling Techniques, 2020, 48(2): 23–28. doi: 10.11911/syztjs.2020034
[3] 张洁,汤明,蒋振新,等. 椭圆井眼同心环空赫巴流体流动规律研究及压降计算简化模型[J]. 特种油气藏,2021,28(2):156–162. doi: 10.3969/j.issn.1006-6535.2021.02.024 ZHANG Jie, TANG Ming, JIANG Zhengxin, et al. Study on flow rules of Herschel-Bulkley fluid in concentric annulus of elliptical wellbore and simplified model for pressure drop calculation[J]. Special Oil & Gas Reservoirs, 2021, 28(2): 156–162. doi: 10.3969/j.issn.1006-6535.2021.02.024
[4] 蔡萌. 幂律流体偏心环空螺旋流压力梯度的数值计算[J]. 石油钻采工艺,2010,32(2):11–14. doi: 10.3969/j.issn.1000-7393.2010.02.003 CAI Meng. Numerical calculation of pressure gradient of helical flow of power-law fluid in eccentric annulus[J]. Oil Drilling & Production Technology, 2010, 32(2): 11–14. doi: 10.3969/j.issn.1000-7393.2010.02.003
[5] BAIRSTOW L, BERRY A. Two-dimensional solutions of Poisson’s and Laplace’s equations[J]. Proceedings of the Royal Society A, 1919, 95(672): 457–475.
[6] HEYDA J F. A green’s function solution for the case of laminar incompressible flow between non-concentric circular cylinders[J]. Journal of the Franklin Institute, 1959, 267(1): 25–34. doi: 10.1016/0016-0032(59)90034-1
[7] REDBERGER P J, CHARLES M E. Axial laminar flow in a circular pipe containing a fixed eccentric core[J]. The Canadian Journal of Chemical Engineering, 1962, 40(4): 148–151. doi: 10.1002/cjce.5450400405
[8] PIERCY N A V, HOOPER M S, WINNY H F. Viscous flow through pipes with cores[J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1933, 15(99): 647–676. doi: 10.1080/14786443309462212
[9] TAO L N, DONOVAN W F. Through-flow in concentric and eccentric annuli of fine clearance with and without relative motion of the boundaries[J]. Transactions of the ASME, 1955, 77(8): 1291–1299.
[10] DRYDEN H L, MURNAGHAN F D, BATEMAN H. Hydrodynamics[M]. New York: Dover Publications, 1956.
[11] VAUGHN R D. Axial laminar flow of non-Newtonian fluids in narrow eccentric annuli[J]. SPE Journal, 1965, 5(4): 277–280.
[12] FERROUDJI H,HADJADJ A,NOFEI T,等. 偏心环空中幂率流体层流流动特性数值模拟研究[J]. 石油钻探技术,2020,48(4):37–42. doi: 10.11911/syztjs.2020066 FERROUDJI H, HADJADJ A, NOFEI T, et al. Effects of the inner pipe rotation and rheological parameters on the axial and tangential velocity profiles and pressure drop of yield Power-Law fluid in eccentric annulus[J]. Petroleum Drilling Techniques, 2020, 48(4): 37–42. doi: 10.11911/syztjs.2020066
[13] MITSUISHI N, AOYAGI Y. Non-Newtonian fluid flow in an eccentric annulus[J]. Journal of Chemical Engineering of Japan, 1973, 6(5): 402–408.
[14] IYOHO A W, AZAR J J. An accurate slot-flow model for non-Newtonian fluid flow through eccentric annuli[J]. SPE Journal, 1981, 21(5): 565–572.
[15] LUO Yuejin, PEDEN J M. Flow of non-Newtonian fluids through eccentric annuli[J]. SPE Production Engineering, 1990, 5(1): 91–96. doi: 10.2118/16692-PA
[16] 刘希圣, 樊洪海, 丁岗. 幂律流体在定向井偏心环空内流动规律的研究[J]. 石油大学学报(自然科学版), 1988, 12(增刊1): 34−45. LIU Xisheng, FAN Honghai, DING Gang. A study on the flow properties of power law fluid in annuli of directional well[J]. Journal of the University of Petroleum, China (Edition of Natural Science), 1988, 12(supplement 1): 34−45.
[17] 刘希圣,崔海清. 幂律流体在倾斜旋转内管的偏心环空中层流流动近似解法[J]. 石油大学学报(自然科学版),1992,16(6):29–34. LIU Xisheng, CUI Haiqing. Approximate solution for laminar flow of power law fluid in inclined eccentric annulus with rotating inner pipe[J]. Journal of the University of Petroleum, China(Edition of Natural Science), 1992, 16(6): 29–34.
[18] 汪海阁,朱明亮. 屈服假塑性流体偏心环空流动的基本特征[J]. 钻采工艺,1997,20(6):5–12. WANG Haige, ZHU Mingliang. Basic characteristics of eccentric annular flow of yield pseudoplastic fluid[J]. Drilling & Production Technology, 1997, 20(6): 5–12.
[19] 汪海阁,苏义脑. Robertson-Stiff流体在偏心环空中的流动[J]. 应用数学和力学,1998,19(10):931–940. WANG Haige, SU Yinao. Flow of Robertson-Stiff fluids through an eccentric annulus[J]. Applied Mathematics and Mechanics, 1998, 19(10): 931–940.
[20] 吴疆. 偏心环空中非牛顿液轴向层流流动规律[J]. 石油钻采工艺,1985,7(2):1–14. doi: 10.13639/j.odpt.1985.02.001 WU Jiang. Law of non-Newtonian fluid axial laminar flow in eccentric annulus[J]. Oil Drilling & Production Technology, 1985, 7(2): 1–14. doi: 10.13639/j.odpt.1985.02.001
[21] HACIISLAMOGLU M. Practical pressure loss predictions in realistic annular geometries[R]. SPE 28304, 1994.
[22] OOMS G, KAMPMAN-REINHARTZ B E. Influence of drillpipe rotation and eccentricity on pressure drop over borehole with Newtonian liquid during drilling[J]. SPE Drilling & Completion, 2000, 15(4): 249–253.
[23] ESCUDIERA M P, GOULDSONA I W, OLIVEIRAB P J, et al. Effects of inner cylinder rotation on laminar flow of a Newtonian fluid through an eccentric annulus[J]. International Journal of Heat and Fluid Flow, 2000, 21(1): 92–103. doi: 10.1016/S0142-727X(99)00059-4
[24] OZBAYOGLU E M, SORGUN M. Frictional pressure loss estimation of non-Newtonian fluids in realistic annulus with pipe rotation[J]. Journal of Canadian Petroleum Technology, 2010, 49(12): 57–64. doi: 10.2118/141518-PA
[25] KHATIBI M, WIKTORSKI E, SUI Dan, et al. Experimental study of frictional pressure loss for eccentric drillpipe in horizontal wells[R]. SPE 191046, 2018.
[26] 樊洪海. 实用钻井流体力学[M]. 北京: 石油工业出版社, 2014. FAN Honghai. Practical drilling fluid mechanics[M]. Beijing: Petroleum Industry Press, 2014.
[27] 韩志勇. 关于钻柱一次弯曲临界条件的考证[J]. 石油钻探技术,2009,37(2):1–4. doi: 10.3969/j.issn.1001-0890.2009.02.001 HAN Zhiyong. An investigation of critical condition of drill string first order buckling[J]. Petroleum Drilling Techniques, 2009, 37(2): 1–4. doi: 10.3969/j.issn.1001-0890.2009.02.001
[28] 黄根炉,韩志勇. 设计井眼中钻柱轴向变形分析与计算[J]. 石油钻探技术,1999,27(6):4–6. doi: 10.3969/j.issn.1001-0890.1999.06.001 HUANG Genlu, HAN Zhiyong. Analysis and calculations of axial deformation of drill string in borehole plan[J]. Petroleum Drilling Techniques, 1999, 27(6): 4–6. doi: 10.3969/j.issn.1001-0890.1999.06.001
[29] 张更,李军,柳贡慧,等. 海上高温高压井环空ECD精细预测模型[J]. 钻井液与完井液,2021,38(6):698–704. ZHANG Geng, LI Jun, LIU Gonghui, et al. A precise model for prediction of annular ECD in offshore HTHP wells[J]. Drilling Fluid & Completion Fluid, 2021, 38(6): 698–704.
-
期刊类型引用(6)
1. 沈吉阳，孙晨皓，王建民，张木楠，邹树江，彭勇. 基于钻柱动力学大位移井钻具接头可靠性研究. 石油矿场机械. 2024(02): 1-9 . 
百度学术
2. 陈雨飞，安锦涛，张辉，李军，周英操，路宗羽. 考虑钻杆接头的小井眼环空压耗计算模型. 钻井液与完井液. 2024(03): 296-304 . 百度学术
3. 王彪，李军，杨宏伟，詹家豪，张更，龙震宇. 基于工程参数变化趋势的溢流早期智能检测方法. 石油钻探技术. 2024(05): 145-153 . 本站查看
4. 赵金，赵星，蔡鹏，彭齐，饶嘉骐. 连续管多孔喷嘴射流工具研制与清洗参数优化. 石油钻探技术. 2023(03): 83-89 . 本站查看
5. 李文拓，罗鸣，黄洪林，李军，肖平. 高温高压小井眼水平井环空ECD综合计算模型. 石油钻采工艺. 2023(03): 259-268 . 百度学术
6. 龚勇，蒋凯，闫许峰，康建平，林黎明. 钻具旋转时同心窄间隙环空流动的数值模拟. 钻探工程. 2023(06): 56-62 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 231
- HTML全文浏览量: 168
- PDF下载量: 62
- 被引次数: 6
