Processing math: 85%

气侵后井底初始气泡平均直径预测模型实验研究

葛磊, 杨春旭, 郭兵, 王志远, 王子毓

葛磊,杨春旭,郭兵,等. 气侵后井底初始气泡平均直径预测模型实验研究[J]. 石油钻探技术,2023, 51(2):46-53. DOI: 10.11911/syztjs.2022096
引用本文: 葛磊,杨春旭,郭兵,等. 气侵后井底初始气泡平均直径预测模型实验研究[J]. 石油钻探技术,2023, 51(2):46-53. DOI: 10.11911/syztjs.2022096
GE Lei, YANG Chunxu, GUO Bing, et al. Experimental study on initial bubble size distribution after gas invading [J]. Petroleum Drilling Techniques,2023, 51(2):46-53. DOI: 10.11911/syztjs.2022096
Citation: GE Lei, YANG Chunxu, GUO Bing, et al. Experimental study on initial bubble size distribution after gas invading [J]. Petroleum Drilling Techniques,2023, 51(2):46-53. DOI: 10.11911/syztjs.2022096

气侵后井底初始气泡平均直径预测模型实验研究

基金项目: 山东省重大科技创新工程项目“一体化智能石油钻井装备技术研究及产业化”(编号:2019JZZY010446)资助
详细信息
    作者简介:

    葛磊(1981—),男,山东烟台人,2004年毕业于石油大学(华东)石油工程专业,高级工程师,主要从事钻井技术与工艺方面的研究工作。E-mail:glhuanying@163.com

  • 中图分类号: TE21

Experimental Study on a Prediction Model for the Average Diameter of Initial Bubbles at Bottom Hole after Gas Cut

  • 摘要:

    为了提高气侵后井筒气液两相流动计算结果的准确性,实验分析了气侵后井底初始气泡直径分布特征,建立了初始气泡平均直径预测模型。用不同质量分数黄原胶溶液模拟钻井液,用多孔介质模拟地层,实验分析了不同液相流变性、地层平均孔隙直径和不同气侵速度下井筒底部气泡群的直径分布特征。实验结果表明:模拟钻井液切力越大、气侵速度越大,生成的初始气泡直径范围越大,出现频率最高的气泡直径和最大气泡直径均增大;地层孔隙直径对初始气泡直径影响不明显。基于实验结果,综合考虑钻井液黏度、气体流量和表面张力等因素的影响,得到了侵入井底初始气泡平均直径实验预测模型;并考虑实际钻井过程中井壁处气体径向侵入特征和井斜角的影响,建立了侵入井底初始气泡平均直径预测模型。井底初始气泡直径预测模型的建立,为气侵后井筒气液两相流动精确计算提供了理论支撑。

    Abstract:

    In order to improve the accuracy of the calculation results of the gas-liquid two-phase flow along the bottom hole after gas cut, the diameter distribution characteristics of initial bubbles at the bottom hole after gas cut were experimentally studied, and a prediction model for the average diameter of the initial bubbles was established. Xanthan gum solutions with different mass fractions were used to simulate drilling fluids, and porous media were utilized to simulate the formations. The bottom hole diameter distribution characteristics of the bubble groups were observed experimentally under different liquid phase rheology, average pore diameter of formation, and gas cut rates. The experimental results showed that a larger yield strength of the simulated drilling fluid and a higher gas cut rate resulted in a larger distribution of diameter range of the initial bubbles generated. In addition, the diameter of bubbles with the highest frequency of occurences and that of the biggest bubbles both increased. However, the formation’s pore diameter had no obvious effect on the initial bubble diameter. According to the experimental results, an experimental prediction model for the average diameter of the initial bubbles invading the bottom hole was established, which comprehensively considered the influence factors such as drilling fluid viscosity, gas flow rate, and surface tension. Furthermore, an prediction model for the average diameter of the initial bubbles invading the well bottom was established by considering the influence of the radial intrusion characteristics of the gas at the borehole wall and the hole deviation angle during the actual drilling. The establishment of the prediction model for the diameter of the initial bubbles at the bottom hole provides theoretical support for the accurate calculation of the gas-liquid two-phase flow along the wellbore after gas cut.

  • 井筒压力安全控制是优快钻井的重要保障技术[1-5],钻井过程中发生气侵后气体通常以气泡形式进入井筒,在井筒内形成气液两相流动体系[6-7]。气泡在钻井液中会发生滑脱上升、膨胀、聚并和悬浮等,导致气液两相流动规律复杂[8-10]。气侵后井底初始气泡直径作为井筒内气体运移的初始条件,是影响井筒气液两相流动精确计算的重要因素之一[11-12]

    国内外对初始气泡直径的研究多集中于孔口单气泡生成系统,多孔介质生成气泡的相关研究还很少。周强等人[13]通过理论推导和实验结果建立了孔口生成气泡脱离周期的演化理论模型。李旭光等人[14]根据实验和量纲分析提出了多孔介质生成气泡初始直径的预测公式,并且定义了多孔介质生成气泡的平均直径计算方法。尹浚羽等人[15]提出了黏弹性流体中孔口生成非球形单气泡的形状预测模型。Zhang Lei等人[16]基于实验结果得到了倾斜多孔介质生成初始气泡的描述模型。N. A. Kazakis等人[17]考虑黏度和表面张力的影响,提出了一种新的预测均质多孔介质气泡初始平均直径的方法。Wang Zhiyuan等人[18-19]考虑液相非牛顿流体特性,建立了屈服应力流体中悬浮气泡形状预测模型和运移模型。

    目前,计算井筒多相流动时大多采用孔口生成单气泡初始直径模型作为井底气泡初始条件[20-22]。实际地层是多孔介质,气侵后井底气体以气泡群形式分散在钻井液中,与孔口生成气泡差异较大。多孔介质中所生成气泡的初始直径研究多基于气−水两相流动,很少考虑钻井液非牛顿流变特性的影响[14-17]。为此,笔者设计并开展了气体通过多孔介质侵入井底实验,分析了不同孔隙直径、气侵速度和钻井液流变性下的井底初始气泡直径特征,并建立了初始气泡平均直径预测模型,为提高气侵后井筒气液两相流动计算精度提供了理论支撑。

    气体侵入井筒后,通常以气泡形式分布于钻井液中,由于气液相间存在密度差,气泡会在浮力作用下沿井筒滑脱上升;气泡上升过程中导致周围液相发生塑性流动,液相质点间切力会阻碍其流动,形成气泡上升过程中的阻力。此外,气体在压差下通过地层多孔介质通道进入井底时存在一个初始动量,也会影响气体的运动状态。

    气泡所受浮力由气体与钻井液之间密度差引起:

    Fb=(ρlρg)gVb (1)
    Vb = 16πD3b (2)

    式中:Fb为气泡所受浮力,N;ρl为钻井液密度,kg/m3ρg为井筒内气体密度,kg/m3g为重力加速度,m/s2Vb为气泡体积,m3Db为气泡沿水平方向的最大宽度,m。

    液相切力是由液相质点间黏性的阻力引起的,由液相本构模型切应力和气泡运动速度决定。液相切力对运动气泡的阻力为:

    FD = π4D2bτ (3)

    式中:FD为气泡所受阻力,N;τ为液相切应力,Pa。

    钻井液切应力与其本构模型相关,通常可以用Herschel-Bulkley模型进行描述:

    τ = τ0+K˙γn (4)

    式中:τ0为钻井液屈服应力,Pa;K为稠度系数,Pa·sn˙γ为剪切速率,s−1n为流性指数。

    τ0=0时,为幂律流体模型;当n=1时,为宾汉流体模型;当τ0=0、且n =1时,为牛顿流体本构模型。

    孔口生成气泡时,气泡脱离前除受浮力和阻力外,还受到孔口与气泡间表面张力和气泡初始动量的影响,其中表面张力作用在气泡上的力为:

    Fs = πdoσ (5)

    式中:Fs为表面张力作用在气泡上的力,N;do为孔口直径,m;σ为溶液表面张力,N/m。

    生成气泡的初始动量与气体注入速度和气泡体积相关:

    Fmt=Vbρg4qgπd2o (6)

    式中:Fm为注入气体动量作用在气泡上的力,N;t为注气时间,s;qg为气体注入流量,m3/s。

    注气时间可由单个气泡体积和注气流量计算得到:

    t=Vbqg (7)

    将式(7)代入式(6),可得注入气体动量作用在气泡上的力为:

    Fm=ρg4q2gπd2o (8)

    由式(1)、式(3)、式(5)和式(8)可知,气泡从地层多孔介质进入井筒底部时所受的力与气泡直径、气液相密度、钻井液流变性、注气孔直径和气侵速度等因素相关。此外,与单个孔口生成气泡相比,多孔介质生成气泡时多个微孔隙渗流出来的气体会汇入同一个气泡中,气泡脱离前尾部存在一个类似孔口的缩径段,表面张力引起的上升阻力存在一定差异,可能与孔隙直径相关。基于上述分析,设计了不同钻井液流变性、不同液相流变性和不同孔隙直径下的气侵后井筒底部气泡初始直径实验。

    侵入井筒初始气泡直径实验装置由支架、实验井筒、图像收集装置、气泡发生装置、减压装置和溶液注入装置等部分组成(见图1)。实验井筒由有机玻璃制成,高1.20 m,内径0.15 m,可承压5 MPa。井筒底部由多孔介质与气源连接,用于模拟地层。气体从井筒底部通过多孔介质注入黄原胶溶液中,来模拟不同液相流变性和气侵速度下的气侵过程,并利用高速摄像机记录上升气泡的大小和形状。实验井筒内置刻度尺,相机焦距定在井筒中的刻度尺上,用于校正因井筒“凸透镜”效应导致水平方向气泡宽度放大引起的测量误差。井筒外部同样有刻度尺,用于读出溶液的液面高度。

    图  1  侵入井筒初始气泡直径实验装置
    Figure  1.  Experimental device for the diameter of initial bubbles invading the wellbore

    图像收集装置由高速摄像机及放置于井筒背面正对镜头的背光灯组成,高速摄像机镜头、背光灯及井筒内刻度尺在同一水平线上,通过调节光照强度可以清晰显示出气液界面,即可得到气泡轮廓。

    气泡发生装置用氮气作为气源,气体通过与注气管线连接的多孔介质进入实验井筒,多孔介质由金属粉末压制而成,直径为10 cm,孔隙直径分别为60和10μm,孔隙直径为多孔介质的平均孔隙直径。假设所有孔隙都是独立的,相互之间无影响,每个孔隙内的气体压力相同。气泡的大小由注入气体的速度及多孔介质的孔隙直径控制[19-20]。为防止注气过程中阀门的开启或关闭导致实验管路底部注气压力的波动,保证注气过程平稳,在注气阀与实验管路之间设置一个阻尼室。

    减压装置主要由实验井筒顶部低压管线和真空泵组成,可以去除注入溶液过程中可能进入的气泡。真空泵最低压力为16 kPa,可以将实验井筒最底部处溶液压力降至负压,去除井筒内全部气泡,减少无关气泡的影响。溶液注入装置的齿轮泵与井筒之间的管路上装有球阀,用于防止溶液回流。

    黄原胶是一种粉末状高分子聚合物,常用作钻井液的增黏剂。黄原胶溶液是典型的非牛顿流体[23],低质量分数状态下黄原胶溶液是无色透明液体,适合拍摄溶液中的气泡形状与运移状态,因此,实验采用黄原胶溶液模拟钻井液。

    配置黄原胶溶液时,将黄原胶粉末缓慢加入去离子水中,加入过程中需注意搅拌机转速,既不能过高导致气体混入,也不能过低导致粉末在水中结成块团,影响后续水解过程。加入黄原胶后,搅拌机保持400 r/min的转速持续搅拌6~8 h,使黄原胶粉末完全分散在水中并充分水解;利用Physica MCR 301流变仪测量溶液的流变性,部分黄原胶溶液的流变曲线如图2所示。

    图  2  不同质量分数黄原胶溶液切应力随剪切速率的变化
    Figure  2.  Yield strength variation of xanthan gum solution of different mass fractions with shearing rates

    实验设计黄原胶溶液质量分数为0.20%~0.50%,剪切速率为1~100 s−1,测试温度与配置溶液温度及实验温度均为22 ℃。分别利用Herschel-Bulkley模型和幂律流体模型对各组溶液流变曲线进行拟合,得到不同质量分数黄原胶溶液的流变参数(见表1)。低质量分数黄原胶溶液的流变性符合幂律流体本构方程,而高质量分数黄原胶溶液流变性则符合Herschel-Bulkley流体本构方程。

    表  1  不同质量分数的黄原胶溶液流变参数
    Table  1.  Rheological parameters of xanthan gum solution with different mass fractions
    黄原胶质量分数,%密度/(kg·m−3)屈服应力/Pa稠度系数/(Pa·sn)幂律指数表面张力/(N·m−1)
    0.201 00401.680.310.064
    0.241 00400.790.310.065
    0.251 0040.150.910.290.065
    0.281 0040.461.380.280.066
    0.301 0041.321.290.270.065
    0.351 0051.961.440.260.066
    0.401 0052.431.280.260.067
    0.451 0052.911.530.240.068
    0.481 0053.241.360.230.067
    0.501 0053.321.770.230.068
    下载: 导出CSV 
    | 显示表格

    1)实验前,将已测好物性参数的黄原胶溶液从实验井筒底部注入井筒中,液面高度达到50 cm时停止注入溶液,检查实验装置气路的密闭性;

    2)溶液静置0.5 h,待气泡自由逸出后,打开氮气瓶的阀门,调节气体流量为0.2 L/min;

    3)使用高速摄像机拍摄记录多孔介质至多孔介质上方5 cm高度范围内的气泡形态;

    4)调节测试气体流量分别为0.4,0.6,0.8和1.0 L/min,重复试验步骤3);

    5)重复实验步骤2)—步骤4),每个气体流量条件实验3次;

    6)排出实验溶液,用去离子水清洗实验井筒后,通氮气干燥实验井筒,进行下一组溶液实验。

    气泡在生成阶段的受力十分复杂,生成过程具有一定随机性,因此,大多数学者关注的是远离多孔介质位置处的气泡直径分布,而忽略了气泡在上升过程中聚并导致的气泡直径分布的变化。绝大多数学者都认为气泡的破裂与聚并实际上发生在多孔介质生成气泡阶段或多孔介质的表面[23-25],远离多孔介质的气泡直径分布并不能代表多孔介质表面生成气泡的直径分布。因此,实验采用高速摄像机拍摄,聚焦在脱离多孔介质瞬间的气泡直径分布上。

    实验用多孔介质的微孔隙平均直径分别为10和60 μm ,黄原胶溶液质量分数分别为0.20%和0.24%、侵入气体流量为0.2 L/min时,井底处初始气泡直径分布如图3图4所示;黄原胶溶液质量分数为0.28%、侵入气体流量为0.8 L/min时,井底处初始气泡直径分布如图5所示。

    图  3  气体流量0.2 L/min、黄原胶溶液质量分数0.20%时的初始气泡直径分布
    Figure  3.  Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.20% at a gas flow rate of 0.2 L/min
    图  4  气体流量0.2 L/min、黄原胶溶液质量分数0.24%时的初始气泡直径分布
    Figure  4.  Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.24% at a gas flow rate of 0.2 L/min
    图  5  气体流量0.8 L/min、黄原胶溶液质量分数0.28%时的初始气泡直径分布
    Figure  5.  Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.28% at a gas flow rate of 0.8 L/min

    图3图4可以看出,气体流量相同时,初始气泡直径分布范围和各直径值对应的概率分布基本相同。从图3图5可以看出,初始气泡直径均呈单峰分布,与对数正态分布类似。黄原胶质量分数增大,气泡直径的上限增大,这是因为黄原胶溶液的表面张力增大;同时,由于屈服应力增大,气泡并不是生成之后直接上升,而是需要增大到一定体积,直至浮力大于屈服应力才会脱离多孔介质上升,但分布方式仍为单峰分布,与牛顿流体中多孔介质孔隙直径大于45 μm时气泡直径呈双峰分布[20]存在较大差异。

    分析认为,造成黄原胶溶液中井底初始气泡直径分布特征与牛顿流体中不同的主要原因是:

    1)黄原胶溶液剪切稀释性的影响。多孔介质生成气泡过程中,气体经多孔介质进入井底黄原胶溶液时不是所有孔隙都会有气体通过,只有部分“有效孔隙”会成为气体通道并参与气泡生成[25],而剪切稀释性会降低气体从孔隙进入溶液时的阻力,使更多气体从“有效孔隙”进入黄原胶溶液中,导致气体难以“开辟”出新的孔隙通道。

    2)黄原胶溶液高黏度的影响。气泡在低质量分数黄原胶溶液中的直径分布范围较小;但随着黄原胶溶液质量分数增大,其黏度增大,气泡脱离多孔介质的时间延长,生成的初始气泡更容易与后续生成的气泡聚并,从而造成气泡的最大直径和直径分布范围均较大。

    3)多孔介质生成气泡的过程的影响。多孔介质生成的气泡直径最小为1.25 mm,比其孔隙直径高出至少2个数量级,而多孔介质的孔隙十分密集,当气体经过多孔介质进入溶液时,会有多个相邻孔隙的气体进入同一个气泡中,形成较大的气泡之后脱离多孔介质并上升。

    综上所述,多孔介质在具有剪切稀释性的非牛顿流体中的生成气泡直径分布主要受黄原胶溶液流变性和多孔介质孔隙直径分布情况影响;孔隙直径增大,会导致气泡直径范围整体右移。

    侵入气体流量影响气泡初始动量,是影响初始气泡直径的另一因素。黄原胶质量分数分别为0.25%和0.30%的溶液在不同侵入气体流量下的侵入井筒初始气泡直径分布如图6图7所示。

    图  6  黄原胶溶液质量分数0.25%、平均孔隙直径60 μm时的初始气泡直径分布
    Figure  6.  Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.25% and an average pore diameter of 60 μm
    图  7  黄原胶溶液质量分数0.30%、平均孔隙直径60 μm时的初始气泡直径分布
    Figure  7.  Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.30% and an average pore diameter of 60 μm

    图6图7可知,随着侵入气体流量增大,侵入井筒初始气泡的直径范围上下限均增大,初始气泡直径分布范围和初始气泡最大概率直径增大。侵入气体流量增大,导致通过多孔介质后初始气泡直径增大,其规律与牛顿流体中相同[23]

    多孔介质孔隙直径与黄原胶溶液质量分数都不变的条件下,多孔介质生成气泡过程中“有效孔隙”的数量并没有增加,侵入气体流量增大后,初始气泡动量增大,上升速度也随之增大,导致气泡脱离多孔介质的频率增大。气泡脱离多孔介质后在黄原胶溶液中上升时与周围气泡发生碰撞概率增加,即气泡聚并的概率增大,最终导致初始气泡直径分布范围和初始气泡最大概率直径增大。

    液相流变性是影响初始气泡运动阻力的主要因素,实验中利用不同质量分数的黄原胶溶液模拟不同流变参数的钻井液。控制侵入气体流量为0.4 L/min,黄原胶溶液质量分数分别为0.20%,0.25%,0.30%,0.35%和0.40%时,溶液中初始气泡直径分布情况如图8所示。

    图  8  不同质量分数黄原胶溶液中的初始气泡直径分布
    Figure  8.  Diameter distribution of initial bubbles in xanthan gum solution with different mass fractions

    图8可以看出,黄原胶质量分数由0.20%增大至0.40%,初始气泡直径分布由1.25~3.25 mm增大为1.75~6.25 mm,初始气泡最大概率直径由1.75 mm增大至3.25 mm;随着黄原胶质量分数增大,侵入井筒初始气泡直径分布范围的上下限均增大,初始气泡直径分布范围和初始气泡最大概率直径增大。

    分析认为,液相流变性对气泡直径分布的影响主要是由切力变化导致的,即相同剪切速率下切力越大,表观黏度越大,井筒中初始气泡运动时产生的阻力也越大。表观黏度对气泡直径有2个方面的影响:1)对最小气泡直径的影响,表观黏度增大,会使生成气泡脱离多孔介质的时间增大,导致初始气泡直径增大;二是对气泡直径分布范围的影响,黏度增大,气泡越容易聚并形成大气泡,使初始气泡直径范围上限和最大概率直径增大。

    分析井筒内生成气泡受力可知,初始气泡直径和分布可能受气液相密度、黄原胶溶液流变性、液相表面张力、侵入气体流量和多孔介质平均孔隙直径等因素影响。实验结果分析可知,多孔介质平均孔隙直径对初始气泡直径的影响最大。为建立各因素对初始气泡平均直径的影响,对实验结果进行量化分析,构建了反映上述物理参数的无量纲参数:

    We=ρlv2sgdwσ (9)
    Re=ρlvsgdwμsl (10)
    Fr=ρlv2sg(ρlρg)gdw (11)
    vsg=4qgπdw (12)
    μsl=τ0+K˙γn˙γ (13)

    式中:We为表征初始动量和表面张力相对大小的无量纲数;Re为表征气泡惯性力和黏性阻力相对大小的无量纲数;Fr为表征惯性力与重力的无量纲数;vsg为井筒中气体表观速度,m/s;dw为实验井筒内径,m;μsl为液相表观黏度,Pa·s。

    利用井筒直径dw为长度标尺对各组实验的初始气泡平均直径进行无量纲化,将气体侵入流量、黄原胶溶液流变性参数代入式(9)、式(11)和式(13),得到各组实验的无量纲数,通过拟合得到侵入井底初始气泡平均直径分布特征实验预测模型:

    dbds=AWeBReCFrD(dofdw)E (14)

    式中:db为侵入井底初期气泡的平均直径,m;ds为井筒实际直径,m;ABCDE为模型系数;dof为实验井筒多孔介质平均孔隙直径,m。

    对实验数据拟合,得到式(14)中各系数,从而得到了无量纲初始气泡平均直径预测模型:

    dbds=3.45×107We0.943Re0.0135Fr0.07866(dofdw)0.0225 (15)

    将模型预测的无量纲初始气泡直径与实验结果进行对比,结果如图9所示,可以看出误差基本在10%以内,表明模型预测结果精确。

    图  9  气泡索特平均直径与无量纲数的关系曲线
    Figure  9.  Relationship between bubble Sauter average diameter and dimensionless number

    实际钻井过程中,地层进入井筒的初始气泡与实验相比存在如下不同之处:1)井底高压条件下,气体密度大,与实验过程中气体密度值相差大;2)井底大部分气泡是井壁径向渗流进入井筒中,与实验过程中气体从井底轴向渗流进入井筒不同;3)井壁近似于圆筒状,而实验过程中多孔介质为圆盘状。因此,需要对基于实验建立的平均气泡直径预测模型进行修正。

    根据井筒中气泡受力分析可知,渗流方向主要影响井筒内气泡的初始动量和表面张力的方向。因此,考虑井斜角度的影响对式(15)进行修正,使其适用于实际钻井中初始气泡平均直径的预测。

    参照包立炯[25]建立的孔口生成初始气泡直径随孔口倾角变化规律描述方法,建立了实际钻井中初始气泡平均直径的预测模型:

    db(β)=f(β)db(0o) (16)
    其中f(β)={1β250,0<β (17)

    式中:dbβ)为井斜角β时生成气泡直径,m;β为井斜角,(°);fβ)为井斜角β时的修正因子;db(0°)为垂直方向生成气泡直径,m。

    井底无钻具时,dw为实验井筒直径;井底有钻具时,需要将井筒直径换算为水力直径dd

    {d_{\text{d}}} = {d_{\text{w}}} - {d_{\text{c}}} (18)

    式中:dd为水力直径,m;dc为井底钻具直径,m。

    假设地层孔隙直径等效于上述实验的多孔介质直径,修正后的井底生成初始气泡的预测模型为:

    \begin{split} \frac{{{d_{\text{b}}}}}{{{d_{\text{s}}}}} = &3.45 \times {10^{{{ - }}7}}f\left( \beta \right)W{e^{0.943}}R{e^{ - 0.013\;5}}\\&F{r^{ - 0.078\;66}}{\left( {\frac{{{d_{\text{p}}}}}{{{d_{\text{w}}}}}} \right)^{ - 0.022\;5}} \end{split} (19)

    式中:dp为地层孔隙直径,m。

    1)多孔介质气侵实验表明,侵入井筒初始气泡直径范围、直径范围上下限和最大概率半径均随气侵速率和钻井液表观黏度增大而增大,与地层孔隙平均直径关系不大。

    2)建立了考虑液相流变性、气侵速率、液相表面张力等因素的初始气泡平均直径预测模型,与实验数据相比模型预测误差小于10%;参照倾角对孔口生成气泡直径的影响,将模型推广到实际钻井工况中。

    3)实际钻井过程中地层孔隙存在非均质性,地层流体在井筒中有溶解、悬浮等现象,都会对侵入井筒初始气泡直径产生影响,建议后续研究中考虑上述因素,建立更贴近实际工况的侵入井筒初始气泡直径预测模型。

  • 图  1   侵入井筒初始气泡直径实验装置

    Figure  1.   Experimental device for the diameter of initial bubbles invading the wellbore

    图  2   不同质量分数黄原胶溶液切应力随剪切速率的变化

    Figure  2.   Yield strength variation of xanthan gum solution of different mass fractions with shearing rates

    图  3   气体流量0.2 L/min、黄原胶溶液质量分数0.20%时的初始气泡直径分布

    Figure  3.   Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.20% at a gas flow rate of 0.2 L/min

    图  4   气体流量0.2 L/min、黄原胶溶液质量分数0.24%时的初始气泡直径分布

    Figure  4.   Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.24% at a gas flow rate of 0.2 L/min

    图  5   气体流量0.8 L/min、黄原胶溶液质量分数0.28%时的初始气泡直径分布

    Figure  5.   Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.28% at a gas flow rate of 0.8 L/min

    图  6   黄原胶溶液质量分数0.25%、平均孔隙直径60 μm时的初始气泡直径分布

    Figure  6.   Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.25% and an average pore diameter of 60 μm

    图  7   黄原胶溶液质量分数0.30%、平均孔隙直径60 μm时的初始气泡直径分布

    Figure  7.   Diameter distribution of initial bubbles in xanthan gum solution with a mass fraction of 0.30% and an average pore diameter of 60 μm

    图  8   不同质量分数黄原胶溶液中的初始气泡直径分布

    Figure  8.   Diameter distribution of initial bubbles in xanthan gum solution with different mass fractions

    图  9   气泡索特平均直径与无量纲数的关系曲线

    Figure  9.   Relationship between bubble Sauter average diameter and dimensionless number

    表  1   不同质量分数的黄原胶溶液流变参数

    Table  1   Rheological parameters of xanthan gum solution with different mass fractions

    黄原胶质量分数,%密度/(kg·m−3)屈服应力/Pa稠度系数/(Pa·sn)幂律指数表面张力/(N·m−1)
    0.201 00401.680.310.064
    0.241 00400.790.310.065
    0.251 0040.150.910.290.065
    0.281 0040.461.380.280.066
    0.301 0041.321.290.270.065
    0.351 0051.961.440.260.066
    0.401 0052.431.280.260.067
    0.451 0052.911.530.240.068
    0.481 0053.241.360.230.067
    0.501 0053.321.770.230.068
    下载: 导出CSV
  • [1] 王果. 基于线性节流阀的MPD井筒压力调控方法[J]. 石油钻采工艺,2021,43(2):197–202. doi: 10.13639/j.odpt.2021.02.010

    WANG Guo. MPD wellbore pressure control method based on linear throttle valve[J]. Oil Drilling & Production Technology, 2021, 43(2): 197–202. doi: 10.13639/j.odpt.2021.02.010

    [2] 史殊哲,吴晓东,韩国庆,等. 内边界旋转的垂直井筒气液两相流型和压降梯度[J]. 石油钻采工艺,2019,41(1):89–95. doi: 10.13639/j.odpt.2019.01.015

    SHI Shuzhe, WU Xiaodong, HAN Guoqing, et al. Flow pattern and pressure gradient of the gas-liquid two-phase flow in the vertical wellbore with rotary inner boundary[J]. Oil Drilling & Production Technology, 2019, 41(1): 89–95. doi: 10.13639/j.odpt.2019.01.015

    [3] 孙宝江,王雪瑞,王志远,等. 控制压力固井技术研究进展及展望[J]. 石油钻探技术,2019,47(3):56–61. doi: 10.11911/syztjs.2019066

    SUN Baojiang, WANG Xuerui, WANG Zhiyuan, et al. Research development and outlook for managed pressure cementing technology[J]. Petroleum Drilling Techniques, 2019, 47(3): 56–61. doi: 10.11911/syztjs.2019066

    [4] 彭明佳,周英操,郭庆丰,等. 窄密度窗口精细控压钻井重浆帽优化技术[J]. 石油钻探技术,2015,43(6):24–28. doi: 10.11911/syztjs.201506005

    PENG Mingjia, ZHOU Yingcao, GUO Qingfeng, et al. Optimization of heavy mud cap in narrow density window precise managed pressure drilling[J]. Petroleum Drilling Techniques, 2015, 43(6): 24–28. doi: 10.11911/syztjs.201506005

    [5] 余意,王雪瑞,柯珂,等. 极地钻井井筒温度压力预测模型及分布规律研究[J]. 石油钻探技术,2021,49(3):11–20. doi: 10.11911/syztjs.2021047

    YU Yi, WANG Xuerui, KE Ke, et al. Prediction model and distribution law study of temperature and pressure of the wellbore in drilling in arctic region[J]. Petroleum Drilling Techniques, 2021, 49(3): 11–20. doi: 10.11911/syztjs.2021047

    [6] 张玉山,赵维青,张星星,等. 深水井气侵在不同类型钻井液中运动特征[J]. 石油钻采工艺,2016,38(3):310–314. doi: 10.13639/j.odpt.2016.03.007

    ZHANG Yushan, ZHAO Weiqing, ZHANG Xingxing, et al. Movement features of gas in different types of drilling fluid during kicking in deepwater wells[J]. Oil Drilling & Production Technology, 2016, 38(3): 310–314. doi: 10.13639/j.odpt.2016.03.007

    [7] 孙士慧,于晓文,于国庆,等. 侵入气体沿井筒的上升规律研究[J]. 数学的实践与认识,2017,47(10):104–111.

    SUN Shihui, YU Xiaowen, YU Guoqing, et al. Rising migration rule of gas kick in wellbore[J]. Mathematics in Practice and Theory, 2017, 47(10): 104–111.

    [8] 孙宝江,王雪瑞,孙小辉,等. 井筒四相流动理论在深水钻完井工程与测试领域的应用与展望[J]. 天然气工业,2020,40(12):95–105. doi: 10.3787/j.issn.1000-0976.2020.12.011

    SUN Baojiang, WANG Xuerui, SUN Xiaohui, et al. Application and prospect of the wellbore four-phase flow theory in the field of deepwater drilling and completion engineering and testing[J]. Natural Gas Industry, 2020, 40(12): 95–105. doi: 10.3787/j.issn.1000-0976.2020.12.011

    [9] 孙士慧. 井底恒压钻井井筒流动模型研究[D]. 大庆: 东北石油大学, 2014.

    SUN Shihui. Research on wellbore flow model for managed pressure drilling at constant bottomhole pressure[D]. Daqing: Northeast Petroleum University, 2014.

    [10] 孟也,李相方,何敏侠,等. 气泡卡断过程中的喉道液领形态与聚并模型[J]. 断块油气田,2019,26(5):632–637.

    MENG Ye, LI Xiangfang, HE Minxia, et al. Shape and coalescence model of liquid collar in pore-throat during snap-off[J]. Fault-Block Oil & Gas Field, 2019, 26(5): 632–637.

    [11] 韦红术,杜庆杰,曹波波,等. 深水油气井关井期间井筒含天然气水合物相变的气泡上升规律研究[J]. 石油钻探技术,2019,47(2):42–49. doi: 10.11911/syztjs.2019035

    WEI Hongshu, DU Qingjie, CAO Bobo, et al. The ascending law of gas bubbles in a wellbore considering the phase change of natural gas hydrates during deepwater well shut-in[J]. Petroleum Drilling Techniques, 2019, 47(2): 42–49. doi: 10.11911/syztjs.2019035

    [12] 段文广,孙宝江,潘登,等. 考虑气体悬浮的关井井筒压力计算方法[J]. 中国石油大学学报(自然科学版),2021,45(5):88–96.

    DUAN Wenguang, SUN Baojiang, PAN Deng, et al. A wellbore pressure calculation method considering gas entrapment in wellbore shut-in condition[J]. Journal of China University of Petroleum(Edition of Natural Science), 2021, 45(5): 88–96.

    [13] 周强,张志东,李淼,等. 早期溢流监测方法[J]. 油气田地面工程,2012,31(3):72–73. doi: 10.3969/j.issn.1006-6896.2012.3.036

    ZHOU Qiang, ZHANG Zhidong, LI Miao, et al. Early overflow monitoring method[J]. Oil-Gas Field Surface Engineering, 2012, 31(3): 72–73. doi: 10.3969/j.issn.1006-6896.2012.3.036

    [14] 李旭光,杨波,孙晓峰. 关井侵入气体运移规律探讨[J]. 长江大学学报(自科版),2013,10(14):65–66. doi: 10.16772/j.cnki.1673-1409.2013.14.023

    LI Xuguang, YANG Bo, SUN Xiaofeng. Discussion on migration law of shut in intrusive gas[J]. Journal of Yangtze University (Natural Science Edition), 2013, 10(14): 65–66. doi: 10.16772/j.cnki.1673-1409.2013.14.023

    [15] 尹浚羽,蒋宏伟,李磊,等. 井底气泡的形成过程及其平均初始直径计算[J]. 科学技术与工程,2014,14(19):30–33. doi: 10.3969/j.issn.1671-1815.2014.19.006

    YIN Junyu, JIANG Hongwei, LI Lei, et al. Bubbles formation from the Bottom region of borehole and their initial average diameter calculation[J]. Science Technology and Engineering, 2014, 14(19): 30–33. doi: 10.3969/j.issn.1671-1815.2014.19.006

    [16]

    ZHANG Lei, SHOJI M. Aperiodic bubble formation from a submerged orifice[J]. Chemical Engineering Science, 2001, 56(18): 5371–5381. doi: 10.1016/S0009-2509(01)00241-X

    [17]

    KAZAKIS N A, MOUZA A A, PARAS S V. Experimental study of bubble formation at metal porous spargers: Effect of liquid properties and sparger characteristics on the initial bubble size distribution[J]. Chemical Engineering Journal, 2008, 137(2): 265–281. doi: 10.1016/j.cej.2007.04.040

    [18]

    SUN Baojiang, PAN Shaowei, ZHANG Jianbo, et al. A dynamic model for predicting the geometry of bubble entrapped in yield stress fluid[J]. Chemical Engineering Journal, 2020, 391: 123569. doi: 10.1016/j.cej.2019.123569

    [19]

    WANG Zhiyuan, LOU Wenqiang, SUN Baojiang, et al. A model for predicting bubble velocity in yield stress fluid at low Reynolds number[J]. Chemical Engineering Science, 2019, 201: 325–338. doi: 10.1016/j.ces.2019.02.035

    [20]

    LOU Wenqiang, WANG Zhiyuan, PAN Shaowei, et al. Prediction model and energy dissipation analysis of Taylor bubble rise velocity in yield stress fluid[J]. Chemical Engineering Journal, 2020, 396: 125261. doi: 10.1016/j.cej.2020.125261

    [21] 尹浚羽,周英操,张辉,等. 钻井环空中赫-巴流体内气泡的上升速度[J]. 西南石油大学学报(自然科学版),2016,38(3):135–143.

    YIN Junyu, ZHOU Yingcao, ZHANG Hui, et al. Bubble rise velocity in Herchel-Bulkley fluid of drilling annulus[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2016, 38(3): 135–143.

    [22] 艾池,张东,任帅勤,等. 环空气泡融合特性对流体流动规律影响研究[J]. 石油地质与工程,2011,25(3):102–105. doi: 10.3969/j.issn.1673-8217.2011.03.032

    AI Chi, ZHANG Dong, REN Shuaiqin, et al. Influence of bubble fusion characteristics in annulus on fluid flow law[J]. Petroleum Geology and Engineering, 2011, 25(3): 102–105. doi: 10.3969/j.issn.1673-8217.2011.03.032

    [23]

    ONG E E S, O’BYRNE S, LIOW J L. Yield stress measurement of a thixotropic colloid[J]. Rheologica Acta, 2019, 58(6): 383–401.

    [24]

    MOUZA A A, DALAKOGLOU G K, PARAS S V. Effect of liquid properties on the performance of bubble column reactors with fine pore spargers[J]. Chemical Engineering Science, 2005, 60(5): 1465–1475. doi: 10.1016/j.ces.2004.10.013

    [25] 包立炯. 毛细管管口气泡生长及脱离特性研究[D]. 重庆: 重庆大学, 2007.

    BAO Lijiong. The study on the characteristics of air bubble growth and departure at the orifice of capillary tubes[D]. Chongqing: Chongqing University, 2007.

  • 期刊类型引用(1)

    1. 蒋振新,李军,郭勇,吴德胜,时培忠,杨宏伟,张更. 井下双梯度控压钻井井筒多相流动规律. 断块油气田. 2024(05): 936-944 . 百度学术

    其他类型引用(0)

图(9)  /  表(1)
计量
  • 文章访问数:  310
  • HTML全文浏览量:  153
  • PDF下载量:  58
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-02-15
  • 修回日期:  2022-11-02
  • 网络出版日期:  2022-12-20
  • 刊出日期:  2023-03-24

目录

/

返回文章
返回