钻井环空中赫—巴流体内气泡的上升速度

尹浚羽*，周英操，张辉，蒋宏伟，朱磊

1.中国石油集团钻井工程技术研究院，北京昌平102206；2.中国石油大学（北京）石油工程学院，北京昌平102249

钻井环空中赫—巴流体内气泡的上升速度

尹浚羽1，2*，周英操1，张辉2，蒋宏伟1，朱磊1，2

1.中国石油集团钻井工程技术研究院，北京昌平102206；2.中国石油大学（北京）石油工程学院，北京昌平102249

气泡在钻井环空中的运动规律是井筒多相流的一个研究重点，气泡上升速度则是其中的一个关键参数。井筒内各流型内气泡/气体上升速度的准确性很大程度上影响了井筒多相流描述的精确性。为此，按照含气率对环空流型进行了划分，对各流型气泡的运动行为进行了分析。系统地总结分析了单气泡在无限流域中的滑脱速度及气泡纵横比的计算方法。利用赫-巴流变模式对气泡周围表观黏度进行了修正，通过实例对比分析单气泡在赫-巴流体内滑脱速度几个计算公式，认为Rodrigue公式相对于Harmathy公式来说是极大的提高，能适应不同形状的气泡，可作为泡状流流域气泡滑脱速度计算的首选公式。通过对单气泡滑脱速度进行含气率、尺寸和井斜修正，得出了环空中各流型的气体滑脱速度计算方法，过渡流型的气泡滑脱速度则采用线性处理，并对环空气泡滑脱速度随含气率变化关系进行了实例分析。最终通过漂移流模型给出了计算环空气泡上升速度的方法。

钻井环空流型；赫-巴流体；气泡滑脱速度；漂移流模型；气泡上升速度

引言

气泡上升运动规律一直是热点研究问题。初期的研究重点是单个气泡的上升过程。单气泡在上升过程中受到浮力、黏滞力、表面张力、惯性力、气泡尺寸、气泡变形率、不稳定扰动量等因素影响。前人通过大量的实验研究和理论分析发展了众多单气泡滑脱速度计算公式。目前钻井环空中泡状流流域的气泡滑脱速度大多采用Harmathy无限流域中单气泡上升速度公式［1］。该公式没有考虑气泡尺寸等因素对上升速度的影响，仅对表面张力占主导地位的流域有较高的精度。若将该公式用于整个泡状流域的气泡滑脱速度计算会造成较大误差。因此，寻找考虑因素较完整、适用范围较广的单气泡上升速度公式是本文的一个重点。

除了以上因素以外，井筒中气泡上升过程会受到含气率的影响。气泡群在上升过程中可能相互碰撞、挤压或融合为一个气泡从而改变气泡尺寸形状，前行气泡形成的尾流会使得后续气泡的速度增大。对于大尺寸气泡（Taylor气泡），井筒尺寸也是一个重要的影响因素。这两个方面的影响，本文会给出相应的影响因子对单气泡上升速度公式进行修正。

值得一提的是大部分钻井液是非牛顿流体，流变模式考虑为赫-巴流变模式。该模式的表观黏度是随着剪切速率变化的。而大部分气泡上升速度计算方法认为气泡周围黏度是常数，故需要进行黏度修正。井筒气液混合流动过程中，针对同一截面液流速度剖面的非线性和气泡分布不均等造成的气泡上升速度变化，本文会采用漂移流模型进行讨论。

综上所述，气泡上升速度是一个多因素影响的函数。本文给出气泡上升速度计算的一般思路是：首先计算出单个气泡在无限流域中的滑脱上升速度，再通过与含气率、流道尺寸和倾角相关的修正因子改进，最后通过漂移流模型得出气泡上升速度。

1 环空气泡形状和流型

图1展示了钻井环空可能出现的气泡形状和流型［2-3］。同一口气井可能自下而上在不同井段分别经历这4个过程。下面将根据Julia JE等［2-3］及Waltrich P J等［4-7］对环空流型的研究介绍各流型特性。

图1 气体形状和流型Fig.1 Flow regimesdefinition

1.1 泡状流

在体积含气率α较小的情况下（一般α＜15%或者α＜20%），液体为连续相，气体以气泡形式分散在液相中。在泡状流型下，气泡有3个基本形状：球形，椭球和球冠状。Clift R［5］利用气泡的Reynolds（Re）数、Eötvös（Eo）数和Morton（M）数绘制了气泡的形状版图以方便使用。

式中：Eo Eötvös数，无因次；

ReReynolds数，无因次；

M-Morton数，无因次；

g重力加速度，g=9.8m/s2；

∆ρ气液密度差，kg/m3；

de气泡等效直径（当量直径），m；

σ气液表面张力，N/m；

ρl液体密度，kg/m3；

U∞气泡在无限流域中稳定滑脱速度，m/s；

µa气泡界面附近的黏度，Pa·s。

当气泡Re＜1时，气泡近似为球形；随着Re和Eo（1＜Eo＜40）增大，气泡横纵比开始增大，气泡形状为椭球形和类椭球形；当气泡Eo＞40时，气泡变形继续增强，气泡形状为球冠形或类球冠形。

1.2 冠状-断塞流

气泡在环空上升过程中，液柱压力逐渐减小，气泡体积膨胀增大，使得气泡含气率增大（20%＜α＜52%）。气泡之间的相互作用也逐渐加强，小气泡开始融合成大气泡。由于管壁的限制，气泡横向膨胀受到限制，纵向则变长。当体积足够大时，气泡便成为了断塞状（图1b）。值得一提的是，环空断塞气泡形状又与圆形截面流道不同：圆形截面流道形成的断塞气泡顶部为球冠状，底部类似为圆柱；而环空中的断塞气泡横向左侧膨胀受到内管壁阻碍，右侧受到外管壁限制，气泡只能沿着环空流道周向膨胀变长，使得气泡截面变为月牙状或者圆弧环状（图1b）。

1.3 搅拌流

气泡继续上升膨胀，气泡含气率进一步加大（52%＜α＜90%），气泡之间的相互作用也更加剧烈，小气泡时而融合形成断塞，断塞时而被湍流震荡冲散为小气泡（图1c）。在搅拌流域，气液交互也十分混乱，液体时而形成段桥阻碍气泡的运动，又被气泡抬升驱散。

1.4 环雾流

气泡含气率达到一定程度（α＞90%），气体流速会非常大，气体和液体运动变为分层流动状态（气泡处于环空中心，液体则在内外管壁上形成液膜），部分液体则以液珠的形式分散在气相中（图1d）。

2 单气泡在赫-巴流体无限流域中的滑脱速度

2.1 气泡表观黏度

钻井液的常用体系绝大部分都是非牛顿流体，而应用范围较广的钻井液流变模式是赫-巴流变模式

式中：τ0屈服应力，Pa；

k 黏度系数，Pa·sn；

γ-剪切速率，s−1；

n-流性指数，无因次。

参照Mendelson H D等对气泡周围幂律流体表观黏度的定义［8］，赫-巴流体中气泡表面的平均剪切速率等于µa/a（a-气泡的横向最大尺寸，m），气泡附近流体表观黏度可写为

2.2 球形气泡滑脱速度和阻力系数

当气泡尺寸较小和气泡雷诺数Re＜1时，表面张力对气泡形状影响占主导地位，气泡形状成球形，而黏性力是气泡最主要的阻力。CliftR等［5］总结了两个适用于球形气泡的两种计算方法。

第一种是Stokes方法，该方法认为小尺寸气泡的气液界面表现出固壁特性。基于这个前提，球形气泡的滑脱速度与球形固体的滑落速度类似

式中：db气泡直径，m。

第二种方法是Hadamard-Rybczynski方法，该方法考虑了气泡上升过程中，黏性切力引起的气泡内部气体循环对气泡滑脱速度的影响。计算公式为

对于稳定上升的气泡来说，气泡处于浮力和上升阻力受力平衡状态，气泡阻力系数CD表示为

对于球形气泡来说，db=de=a，将式（6）代入式（4），得到Stokes阻力系数

将式（6）代入式（5），得到Hadamard-Rybczynski阻力系数

2.3 椭球形气泡滑脱速度、纵横比和阻力系数

（1）椭球形气泡滑脱速度

椭球形气泡一般处于Re＞0.1和Eo＜40的区域。该区域气泡的上升阻力中黏性力比例减弱，表面张力占主导地位，惯性力比例开始增加。忽略黏性阻力的作用，Medelsonb H D［8］创新性地提出了预测气泡上升速度的波动理论，他认为气泡界面扰动的运动特性和理想流体中波的运动特性十分相似，从而导出了气泡的上升速度公式

文献［9］、［10］认为该公式特别适合椭球形气泡滑脱速度计算。Lahrer IH认为上升过程中气泡势能会转换为动能，随后又在波动中耗散掉一部分，这刚好跟Medelson波动理论契合，并对Medelson气泡滑脱速度公式进行了改进［11］

式（10）根号下第一项表征表面张力对气泡滑脱速度的作用，第二项则表征惯性力对气泡滑脱速度的作用。式（10）相对于式（9）在第一项有所强化。

除此之外，CliftR等［5］在前人研究成果基础上，给出了适用于M＜10−3，Eo＜40，Re＞0.1条件的椭球形气泡滑脱速度计算公式

CliftR等［5］总结分析了大量的前人实验数据并得出了气泡纵横比版图，该版图表明，气泡的纵横比不但随着Eo增大而逐渐增大，还跟Morton数密切相关。大Morton数下，黏度大，纵横比小，气泡更不易变形。Bozzano G［12］提出了类似的气泡变形率（a/de）2计算方式（式（17）），从而得到Bozzano和Dente气泡纵横比版图（图2）。

Tadaki和Maeda则认为气泡变形率是与Re和Eo有关的Ta=ReMo0.23的函数，Vakrushev和Efremov扩展了这一方法并得到了适用于M＜10−3的气泡横纵比的计算公式

Vakrushev和Efremov气泡纵横比版图见图3。

图2 Bozzano和Dente气泡纵横比版图Fig.2 Bozzano and Dentebubb le aspect ratio E Vs.Eo

图3 Vakrushev和Efremov气泡纵横比版图Fig.3 Vakrushevand Efremov bubb le aspect ratio E Vs.Re

2.4 球冠形气泡滑脱速度

球冠形气泡一般处于Eo＞40的区域。在该区域，气泡所受的上升阻力中，惯性力占主导地位，黏性力十分微小，表面张力也较小。忽虑黏性力和表面张力的作用，球冠形气泡滑脱速度可以写为［5］

2.5 单气泡滑脱速度的几个通用公式

（1）Karamanev公式

对于小气泡，气泡的滑脱速度遵从Stokes公式。但是随着气泡形状的变化，气泡的阻力系数也随之改变。Karamanev阻力系数公式［13］是利用Reynolds数对Stokes公式进行修正的。

纯净液体中的气泡在上升过程中气泡内部还可能出现内部循环，气泡的阻力系数可利用Hadamard-Rybczynski公式进行修正。而井筒内的钻井液包含了很多表面活性剂和其他杂质，这大大地阻碍了气泡的内部循环，钻井液中的气泡表面更加倾向为“固壁”。阻力系数的修正也应针对于Stokes公式

将式（20）带入式（16）可得到气泡的滑脱速度

（2）Jamaialahamdi公式

Jamaialahamdi等综合考虑Stokes公式和Medelson公式，得出了无限流域中气泡滑脱速度的计算公式［14］

气泡尺寸较小时，Stokes公式占主导地位；当气泡处于大尺寸时，Medelson公式占主导地位。

（3）Bozzano和Dente公式

Bozzano G和DenteM认为气泡阻力系数是摩擦因子和气泡变形率的函数［12］

式中：f 摩擦因子，无因次。

摩擦因子f是Eo和M的函数

（4）Rodrigue公式

Rodrigue［15］通过无量纲分析方法提出了一个很好的单气泡上升速度计算公式

从式（27）可见，该方法考虑到了气泡尺寸、表面张力、黏性和密度等因素，能够在适应非常宽的气泡变形区域和不同的流体介质。

3 不同含气率区域内气泡/气体滑脱速度和井斜校正

3.1 泡状流流域气泡滑脱速度

气泡之间的相互作用会影响气泡的滑脱速度，在泡状流流域（α＜20%），含气率对气泡滑脱速度的影响可用含气率修正因子来表征。根据Richardson和Zaki的启发，气泡的修正因子可写为［14］

式中：Fα含气率修正因子，无因次；

n′系数，无因次。

n′在1.5～2.0取值［16］，本文取为1.75。从而泡状流流域的气泡滑脱速度经含气率修正因子修正后可表示为式中：US滑脱速度，m/s。

3.2 环空Taylor气泡滑脱速度

气泡尺寸和流道尺寸的比例关系也会对气泡的滑脱速度产生影响，气泡尺寸越大，管壁对气泡的上升阻力越大。类似于含气率，流道尺寸对气泡滑脱速度的影响可以流道尺寸修正系数来表征。Krishna R给出了Collins圆形流道尺寸修正系数［17］

式中：Fp流道尺寸修正系数，无因次；

DT圆管内径，m。

当de趋近于DT时，气泡变为Taylor气泡，此时Fp≈0.496。Taylor气泡的前端为球冠状，该结构是影响气泡滑脱速度最主要的因素。因此，圆管内其滑脱速度可以利用流道尺寸修正系数0.496乘以球冠状气泡进行修正得到

Hasan A R等［18］发现环空中的Taylor气泡滑脱速度也和环空几何有关，并扩展得出了环空中Taylor气泡的滑脱速度计算公式（40%＜α＜52%）

式中：Di环空的内径，m；

Do环空的外径，m。

断塞流流域的气泡不仅仅是Taylor气泡，还有一些体积较小的气泡。因此式（32）仅可作为以Taylor气泡为主的断塞流流域（40%＜α＜52%）气泡滑脱速度。

3.3 搅拌流流域的气泡

搅拌流流域的气泡运动则较为不稳定，湍流波动使得气泡在大量破碎的同时又大量融合。对于该流域气泡的滑脱速度，Zuber认为它和Harmathy气泡上升速度公式相同［1，16］

3.4 环雾流流域气体滑脱速度

当环空流域为环雾流流域时，气体为连续相，当含气率趋进于100%时，滑脱速度接近于0。

3.5 过渡含气率区域的气泡滑脱速度

过渡含气率区域为式（31）～式（34）没有覆盖所有的含气率区域。针对这一问题，参考ShiH等［19］的处理方法，对处在过渡区域内的气泡滑脱速度看作是含气率的线性函数，函数的前后两端点分别服从前一含气率区域的气泡滑脱速度公式和后一含气率区域的滑脱速度公式。

3.6 气泡/气体滑脱速度的井斜校正

钻井技术的飞速发展，使得井眼形状多样，钻井环空不再局限于垂直状态，定向井、水平井等迅速发展。这些井型的井斜角必然会对气泡的滑脱速度造成影响。基于实验数据的观察，Hasan A R等［18］得到了气泡/气体滑脱速度的井斜修正因子，修正后气泡/气体的滑脱速度可表示为

式中：θ-井斜角，（◦）。

4 漂移流模型

两相流漂移流模型认为气体和液体之间的滑脱是两个机理联合作用的结果。第一个机理是流体速度剖面非均匀性和流道剖面上气体的非均匀分布使得流道中心处的气体/液体流量最高，混合物的局部速度也最大，对流道剖面积分后得出气体的平均速度会比液体大。第二个机理是由于浮力等作用使得气体超越液体流动。

ZuberN和Findlay JA［20］给出了描述气泡/气体上升速度的漂移流模型

式中：Ug气泡/气体上升速度，m/s；

Co速度分布因子，无因次；

Um气液混合物流速，m/s。

气液混合物流速计算公式为

式中：Ql液体流量，m3/s；

Qg气体流量，m3/s；

A-横截面积，m2。

IshiiM给出了圆管的速度分布系数［21］

式中：ρg气体密度，kg/m3。

环空中的速度分布系数则要考虑环空尺寸。对式（38）进行扩展，高永海［6］给出了环空中各流域气泡的速度分布系数

泡状流域

段塞流域

5 实例分析

根据实际情况，计算钻井过程中的滑脱速度。计算使用的参数见表1。

5.1 泡状流流域中赫-巴流体内单气泡滑脱速度的几个通用公式比较

钻井井筒泡状流流域的气泡滑脱速度多采用Harmathy公式［1，22-23］，在基本参数不变的情况下，该公式计算出的滑脱速度为一常数，如图4虚线所示。与其他气泡计算公式相比较，Harmathy气泡计算公式表征的仅仅是中等尺寸气泡的滑脱速度。因此，利用Harmathy公式计算整个泡状流流域气泡的滑脱速度必然会产生很大的误差。

表1 基本参数列表Tab.1__ Basic parameter list

图4 泡状流流域中单气泡滑脱速度Fig.4 Single bubble slip velocity of bubbly flow pattern

图4展示了4个通用公式计算的滑脱速度随气泡当量直径的变化关系。在气泡当量直径很小（de＜1～2mm）的情况下，气泡的滑脱速度接近于0。主要原因在于赫-巴钻井液特殊的流变模式，赫-巴流体的流变曲线具有较大的屈服值，使得小尺寸气泡周围的表观黏度很高（µa＞100 Pa·s），此时气泡所受到的黏滞阻力会很大，气泡很难向上运动。随着气泡膨胀或气泡之间的融合，气泡当量直径会逐渐增大，赫-巴流体的剪切稀释性显著增强，气泡周围的表观黏度迅速减小（µa＜1Pa·s），气泡在浮力的作用下克服黏滞力向上运动，速度增大幅度也较强。对于小气泡（0～7mm）来说，气泡形状近似为圆球体，黏滞力是最主要的滑脱阻力，气泡滑脱速度遵从Stokes定律（图4中Stokes曲线）。在这一区域与Stokes公式符合最好的通用公式是Rodrigue公式，其次是Bozzano&Dente公式和Jam ialahmadi公式，最差的是Karamaev公式。随着气泡的当量直径继续增大，气泡不再为圆球形，气泡形变增大，形状近似为椭球体，此时，气泡不但受到黏滞力的作用，还受到表面张力和惯性力的影响，气泡滑脱速度的增大幅度减弱。在这区域，Medelson&Lahrer方法计算出的椭球形气泡滑脱速度较之4个通用公式都偏大，原因在于Medelson&Lahrer方法没有考虑黏滞阻力对气泡滑脱速度的影响。这一区域，根据实例计算出的气泡周围的表观黏度量级大于等于10−1，M的量级则大于等于10−2，黏滞力对气泡滑脱运动的阻力不能忽略，Medelson&Lahrer方法计算值偏高。而这种情况不满足Clift关于椭球形气泡滑脱速度计算公式的应用条件（M＜10−3），比较符合的是Rodrigue公式和Bozzano&Dente公式。当气泡体积增大到一定程度，气泡变形率趋于恒定，气泡形状类似于球冠状，黏滞力和表面张力对气泡滑脱运动影响变小，惯性力成为主要的影响因素。气泡滑脱速度接近球冠形气泡的滑脱速度公式计算值（如图4中球冠形气泡滑脱速度公式）。在这一区域符合最好的通用公式是Jam ialahmadi公式，其次是Rodrigue公式和Bozzano&Dente公式。而在各区域表现较差的是Karamaev公式。

综合考虑，本文建议泡状流流域赫-巴流体内单气泡滑脱速度计算使用Rodrigue公式。该公式在各区域的表现都较好，特别是在小气泡区域内计算值和Stokes公式计算值重合度极高。

5.2 钻井环空中赫-巴流体内气泡/气体滑脱速度和上升速度

泡状流流域的气泡运动较为稳定，气泡之间的融合率也很低。气泡尺寸的增大主要是由于气泡的上升过程中周围压力降低引起的膨胀作用。因此，在泡状流流域可以认为气泡含气率和气泡平均体积成正比关系，即气泡的等效直径与含气率的1/3次幂成正线性关系。假设气泡在形成后经过一个不稳定过程变为5mm左右的气泡，此时的含气率为1%。而当含气率1%～20%时，气泡滑脱速度可使用Rodrigue公式，再经含气率校正得到。再分别计算出其他不同含气率区域环空气泡/气体滑脱速度最终得到了气泡含气率与环空气泡滑脱速度的关系曲线（图5）。

从图5中可以看出整个泡状流流域的气泡滑脱速度均小于0.17m/s，如果采用Harmathy公式则计算值在0.25m/s左右。当含气率大于20%时，Taylor气泡开始形成，直到含气率达到40%左右时，Taylor气泡成为环空中气泡的主要组成形式。此时气泡滑脱速度计算采用环空Taylor气泡计算公式较为合理。在含气率大于20%和40%的区间，滑脱速度则认为随含气率线性变化。而一些文献视含气率一旦达到泡状流向段塞流转换边界（α=20%），滑脱速度就可采用环空Taylor气泡公式计算的做法欠考究。当含气率达到52%时，湍流扰动突然加强，Taylor气泡的破碎概率增大，同时伴随着小气泡融合成Taylor气泡过程。当含气率达到60%左右时，这一种扰动现象成为主流，环空流域从段塞流真正过渡到搅拌流。对于过渡区域的处理方法与泡状流向段塞流转换的处理类似。当含气率达到90%时，环空中气体为连续相，部分液体则靠近环空外壁形成液膜，部分液体成为液滴分散在气体中。当含气率达到趋近于1时，环空中为气体单相流动，气体滑脱速度为0。这一区域内，气泡的滑脱速度是随着含气率的升高逐步减小的（图5环雾流区域）。

图5 气体含气率与环空气泡滑脱速度的关系Fig.5 Void fraction vs.bubble slip velocity in annulus

6 结语

（1）根据含气率把环空流型进行了划分：含气率在0～20%为泡状流，20%～52%为断塞流，52%～90%为搅拌流，90%～100%为环雾流。并对各流型的气泡/气体运动行为进行了分析。

（2）对单气泡在赫-巴流体无限流域中的滑脱速度进行了广泛深入的研究。根据赫-巴流变模式对气泡周围的表观黏度进行了修正，分别给出了适合球形气泡、椭球形气泡和球冠形气泡滑脱速度的计算公式以及几个通用公式。通过实例对比分析发现泡状流流域赫-巴流体内单气泡滑脱速度使用Rodrigue公式计算较好。该公式在各区域的表现都较优，特别是在小气泡区域内计算值和Stokes公式计算值重合度极高。该公式相对于Harmathy公式来说有了长足的进步，能作为精确描述井筒泡状流的一个有力工具。

（3）给出了不同含气率区域内气泡/气体滑脱速度的计算公式，对于过度区域的气泡/气体滑脱速度采用线性处理。通过环空气泡滑脱速度随含气率变化关系的实例分析发现这种处理方式能改善一些文献视含气率一旦达到泡状流向段塞流转换边界就采用Taylor气泡计算公式计算滑脱速度等处理方法的缺陷。在得到环空中气泡滑脱速度以后，如果存在井斜，首先利用（35）式进行井斜校正，最后通漂移流模型便可求出气泡的上升速度。

（4）钻井环空中流型的转换规则还需进一步研究。

［1］HARMATHY T Z.Velocity of large dropsand bubbles in media of infinite or restricted extent［J］.AIChE Journal，1960，6（2）：281-288.

［2］JULIA JE，OZARB，JEONG JJ，etal.Flow regime development analysis in adiabatic upward two-phase flow in a verticalannulus［J］.International Journalof Heatand Fluid Flow，2011，32（1）：164-175.

［3］JULIA JE，HIBIKIT.Flow regime transition criteria for two-phase flow in averticalannulus［J］.International Journalof Heatand Fluid Flow，2011，32（5）：993-1004.

［4］WALTRICH P J，FALCONEG，BARBOSA JR.Axial development of annular，churn and slug flows in a long verticaltube［J］.International JournalofMultiphaseFlow，2013，57：38-48.

［5］CLIFT R，GRACE JR，WEBERM E.Bubbles，drops，and particles［M］.New York：ACademic Press，Inc，1978.

［6］高永海.深水钻探井筒多相流流动与井控研究［D］.青岛：中国石油大学（华东），2007.

［7］DRANCHUK P M，PURVIS R A，ROBINSON D B. Computer calculation of natural gas compressibility factors using the Standing and Katz correlation［C］//Annual TechnicalMeeting，Petroleum Society of Canada，1973.

［8］MENDELSON H D.The prediction of bubble term inal velocities from wave theory［J］.AIChE Journal，1967，13（2）：250- 253.

［9］Hassan N M S，Khan M M K，RASUL M G.Bubble risevelocity and drag co-efficientathigh reynoldsnumber in power-law fluids［C］.InternationalConferenceon Fluid Mechanicsand Aerodynamics，Athens，Greece，2007.

［10］NEDELTCHEV S，SCHUMPE A.New approaches for theoreticalestimation ofmass transfer parameters in both gas-liquid and slurry bubble columns［M］.InTech Open Access Publisher，2011.

［11］LEHRER IH.A rational terminal velocity equation for bubblesand dropsatintermediateand high Reynoldsnumbers［J］.Journalof Chem icalEngineering of Japan，1976，9（3）：237-240.

［12］BOZZANO G，DENTEM.Shape and terminal velocity of singlebubblemotion：A novelapproach［J］.Computers &Chem ical Engineering，2001，25（4）：571-576.

［13］KARAMANEVDG，NIKOLOV LN.Free rising spheres donotobey New ton’s law for freesettling［J］.AIChEJournal，1992，38（11）：1843-1846.

［14］SIMONNETM，GENTRICC，OLMOSE，etal.Experimentaldeterm ination of the drag coefficientin a swarm of bubbles［J］.Chem icalEngineering Science，2007，62（3）：858-866.

［15］RODRIGUED.A generalcorrelation for the risevelocity of single gasbubbles［J］.The Canadian Journalof Chem ical Engineering，2004，82（2）：382-386.

［16］WALLIS G B.One-dimensional two-phase flow［M］. America：McGraw-Hill，1969.

［17］KRISHNA R，URSEANUM I，VAN BATEN JM，etal. Risevelocity ofaswarm of largegasbubblesin liquids［J］. Chem icalEngineering Science，1999，54（2）：171-183.

［18］HASAN A R，KABIR C S.Two-phase flow in vertical and inclined annuli［J］.International JournalofMultiphase Flow，1992，18（2）：279-293.

［19］SHIH，HOLMESJA，DURLOFSKY L J，etal.Drift-flux modeling of two-phase flow inwellbores［J］.SPEJournal，2005，10（1）：24-33.

［20］ZUBER N，FINDLAY JA.Average volumetric concentration in two-phase flow systems［J］.Journal of Heat Transfer，1965，87（4）：453-468.

［21］SUN X，KURAN S，ISHIIM.Cap bubbly-to-slug flow regime transition in a vertical annulus［J］.Experiments in Fluids，2004，37（3）：458-464.

［22］KULKARNIA A，JOSHIJB.Bubble formation and bubble risevelocity in gas-liquid systems：A review［J］.Industrial&Engineering Chem istry Research，2005，44（16）：5873-5931.

［23］WANG Z，PEDEN J M，LEMANCZYK R Z.Gas kick simulation study forhorizontalwells［C］.IADC/SPE 27498，1994.

张辉，1971年生，女，汉族，山东平原人，研究员，博士生导师，主要从事油气井信息开发与应用方面的研究。E-mai：zhanghui3702@163.com

蒋宏伟，1978年生，男，汉族，河南鄢陵人，高级工程师，博士，主要从事钻井技术、钻井软件研究。E-mai：jianghwdri@cnpc.com.cn

朱磊，1988年生，男，汉族，河南南阳人，硕士研究生，主要从事油气井工程、钻井软件开发、钻井水力学研究。E-mai：982846086@qq.com

编辑：王旭东

编辑部网址：http：//zk.swpuxb.com

Bubble Rise Velocity in Herchel-Bulk ley Fluid of Drilling Annulus

YIN Junyu1，2*，ZHOU Yingcao1，ZHANGHui2，JIANGHongwei1，ZHU Lei1，2
1.CNPCDrilling Research Institute，Changping，Beijing 102206，China
2.Schoolof Petroleum Engineering，China University of Petroleum（Beijing），Changping，Beijing 102249，China

Bubblemotion laws in drilling annulus is an importantpartof boreholemultiphase flow research.Among which bubble rise velocity isa key parameter.The accuracy of boreholemultiphase flow description highly dependson the accuracy of bubble/gas rise velocity of each flow pattern.In the view of this，according to the void fraction，the drilling annulus flow is divided into four patterns and then the analysis of bubble or gasmotion characteristic of every pattern is carried out.The calculationmethodsof single bubble slip velocity and aspect ratio in Infinitemedia are summarized.Apparentviscosity round bubbles has been corrected through the Herchel-Bulkley reheology model.Through the contrastive analysis of case study among several formulas of single bubble′s slip velocity in Herchel-Bulkley fluid，the Rodrigue′s formula is considered as an improvement compared to the Harmathy′s formula and can be the optimal formula of bubble rise velocity in bubbly flow pattern.Themodificationof calculation formulasofbubbleslip velocitiesw ith regard to void fraction，bubblesizeand theangle ofboreholehasbeen adopted to correct theslip velocities in four flow patterns.For thebubbleslip velocity in the transforming pattern，the linermethod isproposed.Thecasestudy isapplied in the relationship betweenbubble risevelocity and void fraction. Finally，bubble rise velocity can be acquired through themethod of driftmodel

flow pattern in driling annulus；Herchel-Bulkley fluid；bubble slip velocity；driftmodel；bubble rise velocity

尹浚羽，1989年生，男，汉族，四川会东人，硕士，主要从事井筒多相流和钻井软件开发研究工作。E-mail：CUPyinjunyu@163.com

周英操，1962年生，男，汉族，黑龙江大庆人，教授级高级工程师，博士生导师，主要从事钻井工程方面的科研与管理工作。E-mail：zhouyingcaodri@cnpc.com.cn

10.11885/j.issn.1674-5086.2014.02.20.01

1674-5086（2016）03- 0135-09

TE254；O359.1

A

http：//www.cnki.net/kcms/detail/51.1718.TE.20160527.1625.012.htm l

尹浚羽，周英操，张辉，等．钻井环空中赫-巴流体内气泡的上升速度［J］．西南石油大学学报（自然科学版），2016，38（3）：135-143．

YIN Junyu，ZHOU Yingcao，ZHANGHui，etal．BubbleRiseVelocity in Herchel-Bulkley Fluid of Drilling Annulus［J］．JournalofSouthwestPetroleum University（Science&Technology Edition），2016，38（3）：135- 143．

2014- 02-20网络出版时间：2016- 05-27

尹浚羽，E-mail：CUPyinjunyu@163.com

国家科技重大专项（2016ZX05020- 006）。