页岩气水平井暂堵球运移特性数值模拟研究

邹龙庆，何怀银，杨亚东，龚新伟，肖剑锋，苌 北

（1.中国石油集团川庆钻探工程有限公司页岩气勘探开发项目经理部,四川成都 610051；2.油气藏地质及开发工程全国重点实验室（西南石油大学）,四川成都 610500）

我国油气资源开发已经从常规油气转向了页岩、致密砂岩油气等非常规资源，同时基于国家“十四五”发展规划，非常规油气将成为我国油气增储上产的主力。对于非常规油气特别是页岩油气的开发，大规模体积压裂是提高储层改造效果的关键[1-3]。但传统压裂方法成本高、效果差，因此多采用水平井分段多簇压裂技术改造页岩油气储层，实现页岩油气的经济开发[4]。页岩油气储层一般埋藏深、物性差，平面及纵向非均质性强，采用水平井分段压裂技术改造时，由于裂缝间应力相互干扰使多簇裂缝间扩展不均，储层增产改造效果不佳[5]。在此背景下，井筒暂堵压裂技术应运而生，该技术能够有效调控裂缝扩展，从而促使裂缝均衡扩展，有利于提高页岩油气储层的压裂改造效果，许多油田应用该技术改造油气层，并取得了显著成效[6-7]。井筒暂堵是指在非均质性较为严重的储层中投入暂堵球，对优势裂缝孔眼进行暂堵，从而使后续压裂液流向未改造储层，促使劣势裂缝进行扩展，达到分层分段改造、裂缝均衡起裂扩展的目的，从而提高油气产量[8-10]。颗粒类暂堵球是目前封堵裂缝及孔眼常用的材料，其承压能力强，对储层的伤害低，并且在一段时间后能够自行降解，满足施工需求[11]。

页岩储层水平井套管射孔段多簇孔眼螺旋分布，通过压力检测只能粗略判断部分孔眼是否封堵，对暂堵球在水平射孔段运移及坐封规律的研究不足，主要依靠现场经验进行暂堵压裂施工，不能精确预判暂堵球的运移与封堵行为，使投球暂堵后裂缝扩展不均，导致储层改造效果不佳[6]。对于暂堵球在井筒中的运动规律，学者们进行了相关研究，肖辉等人[12]研究发现暂堵球在水平井中的运动主要分为前期的加速运动和后期的匀速运动；蔡华等人[13]发现暂堵球坐封在孔眼上时主要受到来自于射孔压差的持球力和管内流体作用在暂堵球上的脱落力，当持球力大于脱落力时，暂堵球才会坐封在孔眼上；许江文等人[10]通过试验分析了暂堵球形状影响裂缝封堵的规律；刘明明[14]通过研究暂堵球在水平井中的受力得到了暂堵球的运动方程；郑志兵[15]研究了暂堵球的受力情况，以及泵注排量、暂堵球密度、孔眼数量对坐封效率的影响；李春月等人[16]改进了三轴压裂物模试验装置，利用其进行了暂堵转向压裂试验，研究了碳酸盐岩裂缝暂堵转向规律和实现缝内暂堵转向的条件；夏海帮[17]针对页岩气井暂堵后压力升高不明显的问题，优选了暂堵工艺，优化了暂堵参数，形成了页岩气井双暂堵压裂技术，南山页岩气田应用后取得了显著效果；周丹等人[18]基于多级暂堵转向压裂机理研发了新型暂堵材料，形成了多级暂堵转向压裂技术，新疆低渗透油田应用该技术后，产量得到大幅提高；Tan Xuebao 等人[19]对暂堵球运移模型进行了改进，使其适用于直井、斜井和水平井，并建立了暂堵球坐封效率的计算公式。

虽然暂堵球的受力情况及其在井筒中的运移规律得到了一定的研究，但关于暂堵球在水平井中的运移及封堵行为特征的研究较少。目前工程中对暂堵球的投放依旧以现场试验为主，而数值模拟能够有效弥补现场试验的缺陷。一些学者也针对暂堵球运移开展了数值模拟研究，Li Xiaohe 等人[20]建立了暂堵球综合运移模型，研究了影响井筒暂堵的因素，并验证了水平井暂堵的可能性；张锋等人[21]基于CFD-DPM 建立了暂堵球运移模型，分析了暂堵球的运移机理和坐封规律；张雄等人[22]基于复合暂堵压裂机理开展了暂堵压裂数值模拟，优选了暂堵材料，优化了暂堵压裂施工参数，并在塔河油田进行了应用，取得了良好成效；R.W.Brown 等人[23]建立了暂堵球在井筒中运移的模型，探究了泵注排量与流体黏度对暂堵球坐封的影响；周彤等人[24]利用边界元方法，建立了水平井“井筒-射孔炮眼-裂缝扩展”耦合模型，提出了压裂过程中暂堵球分配的计算方法，模拟了非均质场条件下暂堵球数量、投球时机和投球次数对于裂缝扩展的影响；达引朋等人[25]基于三维地质建模分析了含水井调堵压裂增产机理，根据暂堵压裂数值模拟结果优化了暂堵施工参数，并在长庆油田进行了应用，暂堵后缝内压力明显提升；Cheng Wan 等人[26]采用数值模拟方法追踪了暂堵球在螺旋射孔簇水平井中的运移轨迹，分析了暂堵球密度、注入速度及射孔效率对暂堵球坐封效率的影响。

对比各种模拟方法的特征发现，常规离散相模型（DPM）将暂堵球视为质点，在求解过程中忽略了暂堵球的体积及暂堵球自身旋转运动的影响[27]，导致模拟结果误差较大，而CFD-DEM 耦合模型能有效解决上述问题。考虑暂堵球颗粒直径相对较大，笔者基于CFD-DEM 耦合模型，建立了页岩气水平井井筒暂堵模型，模拟了暂堵球在页岩气水平井中的运移过程，研究了暂堵球的运移规律及暂堵球坐封影响因素，从而明确暂堵球坐封孔眼的条件，提升暂堵球坐封效率，为页岩气水平井暂堵压裂施工提供了理论支撑。

1 页岩气水平井井筒暂堵模型的建立

1.1 数学模型

采用CFD-DEM 耦合模型模拟暂堵球在页岩气水平井中的运移时，将液相视为连续相，在欧拉坐标系下求解纳维-斯托克斯方程；将暂堵球视为离散相，在拉格朗日坐标系下求解运动方程。

1.1.1 流体控制方程

由于在孔眼附近压裂液呈紊流状态，因此选择湍流模型作为流体流动模型，故压裂液流动由连续性方程及N-S 方程控制。连续性方程为：

动量守恒方程为：

根据标准k-ε湍流模型，不可压缩湍流的附加方程为：

式中：ρ为流体密度，kg/m3；t为时间，s；u为瞬时速度，m/s；p为压力，Pa；τ为应力张量，kg/（m·s2）；g为重力加速度，m/s2；SM为固液动量交换损失，kg·m/s；µ为层流黏度，Pa·s；µt为湍流黏度，Pa·s；k为紊流脉动动能，J；Gk为平均速度梯度引起的湍流动能，J；Gb为由浮力产生的湍流动能，J；ε为紊流脉动动能的耗散率，%；YM为可压缩湍流中波动膨胀对总耗散率的贡献；C1ε，C2ε和C3ε为常数，C1ε=1.44，C2ε=1.92；σk和σε分别为k和ε的湍流普朗特数，σk=1.0，σε=1.3；Sk和Sε为用户定义的源项。

1.1.2 离散相控制方程

暂堵球作为离散相，遵守牛顿第二定律，可以将暂堵球的运动分解为平滑运动与翻转运动。暂堵球的运动方程[28]为：

式中：mb为暂堵球的质量，kg；ub为暂堵球的平移速度，m/s；Ib为暂堵球的转动动量，kg·m2；ωb为暂堵球的角速度，rad/s；Fw-b为壁面对暂堵球的作用力，N；Fb-b为暂堵球之间的作用力，N；Ff为流体对暂堵球的作用力；Mb为暂堵球所受力矩，N·m。

1.1.3 连续相与离散相的耦合

计算模型中，暂堵球与流体之间的动量交换方程为：

式中：S为连续相到离散相的动量传递量，kg·m/s；μf为流体的黏度，Pa·s；ρb为暂堵球的密度kg/m3；db为暂堵球的直径，m；Re为相对雷诺数；u为暂堵球和流体的速度，m/s；CD为阻力系数；qb为暂堵球的质量流量，kg/s；Δt为时间步长，s；Fother代表单位质量其他作用力，N/kg。

暂堵球与流体之间的质量交换方程为：

式中：V为离散相到连续相的质量传递量，kg/s2；Δqb、qb，0分别为暂堵球的质量流量改变量、暂堵球注入的初始质量流量，kg/s；mb,0为暂堵球的初始质量，kg。

1.2 物理模型

页岩储层等非常规油气藏压裂时多采用水平井进行套管注入，基于现场施工数据以ϕ139.7 mm 套管为研究对象。由于建立完整的井筒模型，需要划分的网格数量过多，计算速度慢，为了提升效率，节约计算资源，建立了典型的单段射孔段套管模型，对单簇射孔段进行模拟。模型总长度1.00 m，套管入口距离射孔段0.40 m，射孔段长度0.50 m，套管外径为139.7 mm，内径为114.3 mm，从左至右8 个孔眼的编号为①，②，…，⑧，射孔沿套管表面呈螺旋分布，孔眼直径为9.5 mm，孔眼相位角为60°，孔眼密度为16 孔/m，如图1 所示。

图1 套管几何模型及孔眼分布示意Fig.1 Geometric model of casing and perforation distribution

在CFD-DEM 耦合的页岩气水平井井筒暂堵物理模型中建立暂堵球球体模型时，考虑了暂堵球自身的体积及旋转运动对暂堵球运移坐封的影响，同时将暂堵球与井筒壁面的接触形式设置为反弹，并考虑了暂堵球碰击孔眼后从孔眼脱落的情况。而以往的研究则将暂堵球视为质点，忽略暂堵球体积以及自身的旋转运动对暂堵球运移的影响，认为从孔眼流出的质点便是坐封在孔眼上的暂堵球，与之相比，该模型能更准确地反映暂堵球在水平射孔段的运移坐封情况。模型边界条件设置如图2 所示，套管入口边界条件设置为“速度入口”，孔眼出口边界条件设置为“压力出口”，泵注排量设置为1～13 m3/min，暂堵球初始注入速度与压裂液速度一致，压裂液密度为1 000 kg/m3，压裂液黏度为1 mPa·s，暂堵球密度分别为900，1 020 和1 200 kg/m3，暂堵球数量为1/10 个。

图2 边界条件设置Fig.2 Model boundary condition setting

2 模拟结果分析与讨论

2.1 暂堵球运动过程分析

井筒暂堵是指储层改造过程中向井筒中泵注含有暂堵球的压裂液，暂堵球随压裂液运移到水平井射孔段，并坐封在孔眼上封堵孔眼，从而使压裂液转向未改造的储层段[29]。暂堵球在水平井中运移直至坐封在孔眼上主要分为3 个阶段（见图3）：非射孔段运移、射孔段运移、孔眼坐封。暂堵球在孔眼处的受力主要为拖曳力、惯性力、持球力、脱落力（见图4）。当拖曳力大于惯性力时，暂堵球被流体携带向孔眼方向运移，而暂堵球坐封在孔眼上时主要受到来自于射孔压差的持球力和管内流体作用在暂堵球上的脱落力，当持球力大于脱落力时，暂堵球才会坐封在孔眼上[30]。

图3 暂堵球坐封孔眼过程示意Fig.3 Sealer ball setting perforation process

图4 暂堵球在孔眼处的受力分析Fig.4 Force analysis of the ball sealer at perforation

惯性力的计算公式为：

拖曳力的计算公式为：

持球力的计算公式为：

脱落力的计算公式为：

式中：Δp为孔眼压差，MPa；dp为孔眼直径，m；fd为阻力因子；FG为暂堵球受到的惯性力，N；Fr为暂堵球受到流体的拖曳力，N；Kd为拖曳力系数；FH为暂堵球受到的孔眼压差形成的持球力，N；Fu为暂堵球受到来流冲击的脱落力，N。

2.1.1 单个暂堵球运移封堵特征分析

采用CFD-DEM 耦合方法，可以对暂堵球在页岩气水平井中的运移过程进行可视化数值模拟，泵注排量设置为2 m3/min，暂堵球初始运移速度设置为3.05 m/s（与流体速度一致），对单个高密度暂堵球的运动过程进行分析，其中暂堵球运移速度变化及暂堵球运移轨迹如图5 所示，某时刻暂堵球运移位置如图6 所示，流场压力流线如图7 所示。

信息时代，学生每天都会浏览大量信息。通过时事新闻评论这一环节，教师可以了解有哪些时事新闻能引起学生的兴趣，并通过学生的分析评价来判断学生的价值观。“对于一些比较敏感、学生争议比较大的问题，教师要坚持正确的政治导向，对播报中原则性的错误，特别是学生认识上的偏差一定要给予纠正，以便充分发挥思想政治理论课主渠道、主阵地的作用，帮助大学生树立正确的世界观、人生观和价值观。”[2]因此，“概论”课不仅要提高大学生的知识水平和理论素养，更重要的是增加其理论自信、道路自信和制度自信。

图5 单个暂堵球的运移速度和运移轨迹Fig.5 Migration speed and trajectory of a ball sealer

图6 特定时刻单个暂堵球的运移位置示意Fig.6 Movement position of a ball sealer at a specific moment

图7 单个暂堵球时的流场压力流线Fig.7 Pressure streamline of flow field of a ball sealer

从图5（a）和图6（a）可以看出：暂堵球进入射孔段之前（t＜0.1 s），其运移速度略微增大，因为在未进入射孔段之前，暂堵球与流体流动方向一致，流体对于暂堵球存在一定的推动作用，使暂堵球运移速度升高，暂堵球能更快进入射孔段；暂堵球在进入射孔段后，其运移速度略微降低，因为暂堵球运移速度大于流体速度后，流体对暂堵球存在一定的阻力作用，使暂堵球运移速度降低。从图5（a）和图6（b）可以看出，暂堵球进入射孔段后（t＞0.14 s），由于流体不断从孔眼流出，井筒内流量减小，暂堵球运移速度显著降低，导致惯性力降低。同时，结合图7 可以看出，由于流体不断从孔眼流向储层，在孔眼处形成负压差，从而对暂堵球产生拖曳力，使暂堵球在流体携带作用下向孔眼附近运移（见图5（b））。从图5（a）和图6（c）可以看出，因为孔眼截面积较小，暂堵球被流体携带至孔眼附近时，暂堵球运移速度有所升高，使暂堵球因井筒流量降低而导致其运移速度降低的趋势减缓，但在离开孔眼拖曳力作用范围后暂堵球运移速度显著降低。从图5（a）和图6（d）可以看出，当暂堵球运移至⑥号孔眼处时，暂堵球速度因为孔眼截面减小而增大，并且以较大的速度碰击孔眼⑥，因为受到孔眼压差的束缚作用，暂堵球运移速度降低至0，从而坐封在孔眼⑥上。

2.1.2 多个暂堵球运移封堵特征分析

页岩储层在进行暂堵转向压裂施工时，一般采用多个暂堵球对多个孔眼坐封，因此，采用数值模拟分析多个暂堵球的运移坐封过程。泵注排量设置为6 m3/min，暂堵球速度设置为9.14 m/s（与流体速度一致），对多个低密度暂堵球的运动过程进行模拟分析，其中暂堵球的运移轨迹与坐封孔眼位置如图8 所示，坐封在不同孔眼上暂堵球的运移速度如图9 所示，某时刻暂堵球的运移位置如图10 所示，流场压力流线如图11 所示。

图8 多个暂堵球的运动轨迹与坐封孔眼位置Fig.8 Migration trajectory of ball sealers and sealed perforation position

图9 坐封在不同孔眼上暂堵球的运移速度Fig.9 Migration speed of ball sealers on different perforations

图10 特定时刻多个暂堵球的运移位置Fig.10 Migration position of ball sealers at a specific time

图11 多个暂堵球时的流场压力流线Fig.11 Pressure streamline of flow field of ball sealers

从图8（b）可以看出，排量为6 m3/min 时，暂堵球能坐封在孔眼⑥、⑦和⑧上，因此适当增加暂堵球数量，有利于增加坐封孔眼数量，提升井筒暂堵效果。结合图9、图10（a）和图10（b）可知，暂堵球进入射孔段前（t＜0.04 s），暂堵球的运移速度有所增大，因为在该时间段内暂堵球与流体流动方向一致，流体对暂堵球有一定推动作用，使暂堵球运移速度有所升高，随后由于流体速度比暂堵球运移速度低，对暂堵球产生阻力作用，暂堵球运移速度略微降低；暂堵球进入射孔段后（t＞0.05 s），由于流体不断从孔眼流向储层，井筒内流量降低，暂堵球运移速度降低，惯性力也降低。从图8（a）和图11可以看出，流体从孔眼流出，在孔眼处形成了负压差，产生了拖曳力，使暂堵球在流体携带作用下向孔眼运移，便于暂堵球坐封在孔眼上。从图10（c）、图10（d）、图10（e）可以看出，t=0.10，0.15 和0.18 s时，暂堵球逐渐坐封在孔眼⑥、⑦和⑧上，结合图9分析得到暂堵球被流体携带至孔眼⑥、⑦和⑧附近时，由于孔眼截面较小，暂堵球运移速度有所升高，使暂堵球因井筒流量降低，导致其运移速度降低的趋势有所减缓，同时能以较大的运移速度碰击孔眼⑥、⑦和⑧，并且在孔眼压差束缚下暂堵球运移速度降低至0，从而坐封在孔眼⑥、⑦和⑧上。

2.2 暂堵球运移坐封影响因素分析

2.2.1 暂堵球直径对暂堵球运移以及坐封的影响

暂堵球在孔眼处的坐封一般经历封堵、溶胀和变形3 个过程，因此暂堵球直径不能小于孔眼直径，否则暂堵球直接从孔眼流出而无法坐封在孔眼上[31]。暂堵球的直径对暂堵球受到的重力、拖曳力、惯性力等均有影响，因此基于数值模拟研究暂堵球直径对于暂堵球运移和坐封的影响。暂堵球直径设置为9.5～19.0 mm，泵注排量为4 m3/min，模拟孔眼直径与暂堵球直径之比在0.5～1.0 条件下暂堵球运移坐封情况，暂堵球坐封孔眼模拟结果见表1 和图12，坐封效率（坐封孔眼数量与孔眼总数量之比）和孔眼与暂堵球直径之比的关系如图13 所示。

表1 暂堵球直径对暂堵球运移及封堵行为的影响Table 1 Influence of ball sealer diameter on migration and sealing behavior of ball sealer

图12 暂堵球直径对暂堵球运移及封堵行为的影响Fig.12 Influence of ball sealer diameter on migration and sealing behavior of ball sealer

图13 坐封效率和孔眼与暂堵球直径之比的关系Fig.13 Relationship between sealing efficiency and ratio of perforation diameter to ball sealer diameter

结合表1、图12 和图13 可知：孔眼与暂堵球直径之比大于0.83 时，暂堵球才能坐封在孔眼上；孔眼与暂堵球直径之比小于0.97 时，由于暂堵球受到的阻力与暂堵球直径的平方成正比，此时暂堵球直径增大，暂堵球所受阻力也随之增加，且增加得更快，使暂堵球坐封难度增大，坐封效率降低；孔眼与暂堵球直径之比为0.97 时，暂堵球坐封孔眼数量最多，坐封效率最高；孔眼与暂堵球直径之比大于0.97 时，随着暂堵球直径增大，流体对暂堵球的拖曳力增大，使暂堵球更容易向孔眼处运移，暂堵球坐封效果更好。

2.2.2 泵注排量对暂堵球运移以及坐封的影响

暂堵球在孔眼处坐封主要受惯性力、拖曳力、持球力和脱落力的影响[26]。当拖曳力大于惯性力时，暂堵球被流体携带向孔眼附近运移，只有当持球力大于脱落力时，暂堵球才能坐封在孔眼上，不会在来流的冲击下从孔眼脱落[13]；而惯性力、拖曳力、持球力及脱落力均与泵注排量有关，因此基于数值模拟分析泵注排量对暂堵球运移及坐封的影响。泵注排量设置为1～13 m3/min，选择低密度暂堵球（900 kg/m3）、中密度暂堵球（1 020 kg/m3）、高密度暂堵球（1 200 kg/m3），模拟不同密度暂堵球在不同泵注排量下的运移封堵情况，暂堵球坐封孔眼模拟结果见表2，暂堵球坐封效率与泵注排量的关系如图14 所示。

表2 泵注排量对暂堵球运移及封堵行为的影响Table 2 Influence of pump rate on migration and sealing behavior of ball sealer

图14 暂堵球坐封效率与泵注排量的关系Fig.14 Relationship between sealing efficiency of ball sealer and pump rate

从表2 可以看出：暂堵球在水平井运移过程中，暂堵球坐封的孔眼大多位于水平段下半部分，表明重力对于暂堵球运移及封堵的影响较大，暂堵球投入井筒后，其在重力作用下沿着套管下半部分运移，并坐封在套管下半部分的孔眼上；只有采用低密度和中密度暂堵球才坐封在位于套管中上部的孔眼⑤上，因为此时暂堵球密度与流体密度相差较小，暂堵球在浮力作用下能在套管中部运移，从而能坐封在孔眼⑤上。从图14 可以看出：暂堵球坐封效率随着泵注排量增大呈先升高后降低的趋势；在低泵注排量范围内（1～6 m3/min），逐渐增大泵注排量，有利于提升暂堵球坐封效率。这是因为在低泵注排量范围内，暂堵球所受惯性力较小，逐渐增大泵注排量，则孔眼流量增加，导致孔眼压差增大，拖曳力也会相应增大，流体更容易携带暂堵球向孔眼运移，并且在孔眼处，因为孔眼流量增大，持球力也会增大，因此暂堵球不容易从孔眼处脱落，坐封效率更高；当泵注排量持续增大时，暂堵球所受惯性力增大，并逐渐占主导地位，此时惯性力对于暂堵球运移的影响更大，因此暂堵球不易向孔眼运移，并且孔眼处来流的冲击力增强，脱落力增大，使暂堵球更容易从孔眼上脱落，导致坐封效率降低。从图14 还可以看出，泵注排量在4～7 m3/min 时，暂堵球的坐封效率相对较高。

2.2.3 暂堵球密度对暂堵球运移以及坐封的影响

水平井进行分段多簇压裂时，不同水平井段所采用的暂堵球密度可能不一致[32]，而暂堵球密度对于暂堵球所受重力和惯性力均有影响，从而影响暂堵球在水平段的运移。采用数值模拟方法，分析暂堵球密度对暂堵球运移及坐封的影响，结果见表3。

表3 不同密度暂堵球在不同泵排量下的坐封效率Table 3 Sealing efficiency of ball sealers with different densities under different pump rates

从表3 可以看出：泵注排量为1～2 m3/min时，低、中、高密度暂堵球均能坐封在孔眼上，但坐封效率低，这是因为泵注排量小，暂堵球运移速度小，惯性力小，容易被流体携带向孔眼运移，但由于泵注排量小，各个孔眼流量低，孔眼压差小，暂堵球难以坐封在孔眼上，因此坐封效率低；泵注排量为3～8 m3/min 时，暂堵球坐封效率明显升高，因为泵注排量增大后，孔眼流量增大，拖曳力增大，暂堵球更容易被流体携带至孔眼处，并且孔眼压差增大后各孔眼持球力增大，使坐封效率升高；泵注排量为9～12 m3/min 时，暂堵球运移速度增大，惯性力增加得更快，暂堵球运移轨迹难以发生改变，流体难以携带暂堵球至孔眼处，此外泵注排量增大后，井筒内流体对孔眼处暂堵球产生的脱落力增大，使坐封效率降低；泵注排量大于13 m3/min 时，暂堵球运移速度极大，在惯性力的作用下流体几乎对暂堵球没有影响，暂堵球直接经过射孔段流出井筒而无法坐封在孔眼上。从表3 还可以看出：低密度暂堵球坐封在孔眼上的泵注排量范围最大，中密度暂堵球次之，高密度暂堵球最小；在中等排量范围内，3 种密度暂堵球的坐封效率均较高，但低密度暂堵球的坐封效率普遍比中密度暂堵球和高密度暂堵球高，综合分析认为，暂堵球的密度越大，暂堵球的重力越大，则暂堵球受到的惯性力越大，从而导致暂堵球在运移过程中流体对暂堵球的携带能力减弱，其运动轨迹难以发生改变，暂堵球不易向孔眼运移，坐封效率不高。

2.3 暂堵球坐封孔眼位置特征分析

对暂堵球运移轨迹进行分析发现，暂堵球在孔眼处的运动存在2 种情况：1）暂堵球碰击孔眼之后恰好坐封在孔眼上；2）暂堵球碰击孔眼之后在孔眼处短暂停留，随后在来流的冲击作用下从孔眼脱离继续流动，导致暂堵球无法坐封在孔眼上。对孔眼坐封次数进行统计，结果如图15 所示。从图15 可以看出，暂堵球在页岩气水平井射孔段运移时最容易坐封在射孔簇最后3 个孔眼处（见图2），因为流体不断从孔眼中流出，使射孔段后半部分井筒的流量和流速降低，暂堵球运移速度及其所受惯性力更小，暂堵球容易被流体携带向孔眼运移，因此更容易在射孔簇后半段的孔眼坐封。此外，暂堵球易坐封在第1 个孔眼上，其原因是暂堵球刚进入射孔簇时，惯性力较大，并且在拖曳力作用下，暂堵球以较高的速度碰击第1 个孔眼，由于其直径较大，所以容易被该孔眼卡住从而坐封在其上。

图15 孔眼暂堵球坐封次数统计结果Fig.15 Statistical results of sealing times of perforation by ball sealer

3 结论

1）暂堵球在页岩气水平井运移过程中，进入射孔段之前，暂堵球在流体的推动作用下运移速度增大，暂堵球加速进入射孔段。进入射孔段之后，由于流体不断从孔眼流向储层，在孔眼处产生拖曳力，而井筒内流量降低，使暂堵球运移速度逐渐降低，惯性力随之降低，因此在拖曳力作用下，暂堵球更容易被流体携带至孔眼附近；由于孔眼截面较小，孔眼附近的暂堵球运移速度升高，减缓了暂堵球因为井筒流量降低而导致运移速度降低的趋势，暂堵球运移速度变化幅度变小，同时也有利于暂堵球以较大的运移速度碰击孔眼从而坐封在孔眼上。

2）孔眼与暂堵球直径之比小于0.83 时，暂堵球在来流冲击作用下无法坐封在孔眼上；孔眼与暂堵球直径之比大于0.83 时，随着孔眼与暂堵球直径之比增大，暂堵球的坐封效率呈现先升高后降低的趋势；孔眼与暂堵球直径之比为0.97 时，暂堵球坐封孔眼的数量最多，暂堵球的坐封效率最高。

3）暂堵球坐封效率随着泵注排量增加呈先升高后降低的趋势。泵注排量较低时（1～6 m3/min），逐渐增大泵注排量，暂堵球坐封效率升高；当排量持续增加时，暂堵球的惯性力增加，此时暂堵球的运移轨迹不易发生改变，坐封效率降低；泵注排量为4～7 m3/min 时，暂堵球坐封效率相对较高。

4）低密度暂堵球坐封效果最好，高密度暂堵球坐封效果最差。

5）暂堵球最容易坐封在射孔簇后半段孔眼处，也能坐封在第1 个孔眼上。