环空涡动流场特性分析范文

本站小编为你精心准备了环空涡动流场特性分析参考范文，愿这些范文能点燃您思维的火花，激发您的写作灵感。欢迎深入阅读并收藏。

《石油钻采工艺杂志》2016年第二期

摘要：

为了深入研究环空涡动时的流场特性，根据流体力学理论，以连续性方程和N-S方程为控制方程，利用计算流体力学技术对钻柱涡动时环空赫巴流体的流动进行了系统的数值模拟，研究了不同公转方向下钻柱偏心度的变化对环空切向速度剖面与合速度剖面的影响。通过对模拟数据的对比分析，发现不同公转方向下环空流场分布截然不同，正向公转时切向速度在环空宽间隙处随着公转速度、自转速度和偏心度增大正向增大；反向公转时会出现二次流，切向速度在环空宽间隙处随偏心度的减小整体反向减小，同时二次流趋势越明显，摩阻压耗越大。合理应用这些规律有助于完善现有钻井水力学理论，更好地揭示井下环空流场特性，并为钻井水力参数优化设计提供有效的理论指导。

关键词：

涡动；环空；赫巴流体；摩阻压耗；数值模拟；流场

钻进过程中钻柱在井筒内一般会出现涡动，即在地面转盘带动下顺时针自转的同时还绕井眼轴线公转。涡动的出现以及公转的半径（即涡动时钻柱的偏心度）和方向（顺时针方向、逆时针方向）就目前钻井水平而言无法控制且很难预测，通过实验难以揭示其内部规律，最可行的方法是借助数值模拟技术进行研究［1-5］。因此研究涡动时的偏心度对由自转旋流场、公转旋流场和轴向流场叠加而成的环空真实流场［6-12］的影响有助于完善现有钻井水力学理论，更好地揭示井下环空流场特性，并为钻井水力参数优化设计提供理论指导［13-16］。

1环空流动的建模与模拟

1.1钻井液模型为更好地描述非牛顿不可压缩流体的流变特性，采用了更为合理的赫巴流变模式，该模式也可看作修正的幂律模式。

1.2几何模型数值模拟所涉及到的几何尺寸均采用现场实际尺寸，钻柱外径139.7mm，套管内径244.0mm，环空长度15m，其他参数见表1。选取15m长的环空钻井液作为研究对象，采用基于交错网格的SIMPLEC算法对连续性方程和N-S方程进行离散化，动量方程的对流项进行离散化时为保证计算精度，采用二阶格式；建立的网格几何模型如图1。设置3类边界类型：速度入口边界，压力出口边界以及旋转壁面边界，设定该边界绕其自身轴线自转的同时绕井壁轴线公转，井壁为静止壁面，所有壁面都简化为无滑移水力光滑的壁面，且近壁面处的网格加密划分，采用柱坐标系，对称轴为z轴，按照速度分量方式设定速度入口值［20］。

1.3流动模型将流体视为不可压缩流体，且具有一定的黏度值和密度值，不考虑传热和重力的影响［21］。用连续性方程和N-S方程共同描述环空流动，基于有限体积法的思想进行离散。由于涡动时的环空流场属于复杂的旋流场，流体质点脉动程度较高，因此采用湍流两方程模型中的Realizable模型，该模型能够有效地将湍动黏度模型中的系数与应变率耦合起来，能够更为精确地描述包括旋转均匀剪切流、包含有射流和混合流的自由流、管道内流动、边界层流动以及带有分离的流动。因此文中采用Realizable模型［22］。

1.4环空流动的数值模拟涡动时自转和公转坐标系都遵循右手准则，规定转盘自转为顺时针方向，以此作为旋转正方向，钻柱在井眼内的公转根据已有研究的描述有顺时针和逆时针2个方向，当公转为顺时针时，与自转方向一致，根据右手准则，为正向涡动，反之为反向涡动。为更好地研究不同偏心度下流场的变化规律，设置其余钻井参数在合理范围，如表2。取15m长偏心环空钻井液距入口9m处横截面上沿对称轴的剖面进行分析，横坐标0表示井筒中轴线。

2钻柱涡动时环空流场特性分析

设井筒圆心为O1，钻柱圆心为O2，井筒半径为R，钻柱半径为r，偏心距为h，公转速度为ω1，自转速度为ω2，P1为最宽间隙处对应的钻柱外壁面位置，P2为最窄间隙处对应的钻柱外壁面位置。设图2中时刻的自转线速度方向为正方向，根据分析可知：

2.1偏心度对流速剖面的影响分析模拟数据发现，其他条件不变时，随轴向流速增大，环空内宽窄间隙处合速度剖面变化规律与一般环空流动相同，由于径向速度较小且变化不明显，本研究未作讨论。下面着重分析涡动时的偏心度及公转正反方向的变化对环空气场的影响。设定环空钻井液密度为1.0g/cm3，轴向流速为2.0m/s，公转速度为18.84rad/s，自转速度为18.84rad/s，考查偏心度分别为0.32、0.64、0.95时各流动参数剖面的变化规律。切向速度剖面：偏心度增大，宽间隙处切向速度剖面变大变宽，窄间隙处的切向速度剖面整体增大，钻柱外壁面在宽间隙处根据式（3）的推导切向速度越小，在窄间隙处根据式（4）的推导切向速度越大，如图3a。合速度剖面：偏心度增大，宽间隙处合速度剖面整体增大，窄间隙处合速度剖面整体减小，钻柱外壁面在宽窄间隙处的变化规律同切向速度的变化规律，如图3b。设定环空钻井液密度为1000kg/m3，轴向流速为2.0m/s，公转速度为–18.84rad/s，自转速度为18.84rad/s，考查偏心度分别为0.32、0.64、0.95时各流动参数剖面的变化规律。切向速度剖面：随偏心度增大，宽间隙处的切向速度剖面整体反向增大，窄间隙处的切向速度剖面先增大后减小，钻柱外壁面在宽间隙处根据式（5）的推导切向速度正向增大，在窄间隙处根据式（6）的推导切向速度减小，井壁处切向速度为0，如图3c。合速度剖面：随偏心度增大，宽间隙处合速度剖面整体增大，窄间隙处的合速度剖面整体减小，钻柱外壁面宽窄间隙处的变化规律同切向速度的变化规律，如图3d。

2.2二次流的分析反向公转时，宽间隙钻柱外壁面附近的流体质点跟随钻柱自转呈正向旋转，此处的切向速度为正值，但在远离钻柱外壁面的区域，由于钻柱正向自转的趋势弱于反向公转的趋势，故此区域的流体质点跟随钻柱反向公转，所以此处切向速度为负值，如图3c所示。在宽间隙区域，同时出现正反2个方向的切向速度，两股不同方向的旋流在环空窄间隙的两侧分离与汇合，使得环空宽间隙处出现一个明显的大涡，环空窄间隙处由于钻柱旋转壁面带动的旋流将本应对应出现的另一个不对称的小涡淹没，最终形成了不对称的二次流，如图4所示。反向公转时，其他参数不变，偏心度越大，越不容易出现二次流。宽间隙处二次流产生的大涡越靠近旋转钻柱外壁面，同时二次流的分离与汇合点越靠近窄间隙处，当偏心度大到一定程度后，钻柱的正向自转趋势被二次流在宽窄间隙处的反向流动趋势完全包裹且基本抵消掉，此时状态与偏心环空自转的螺旋流场（自转方向为反向，自转速度大小同公转最为接近。二次流的明显出现使得环空流体质点的脉动趋势进一步加强，随着偏心度减小，钻柱正向自转趋势越不容易被抵消，二次流趋势越明显，偏心环空流动的摩阻压耗将进一步增大。

3结论

（1）偏心度对涡动流场的影响较为复杂，不能轻易忽略。宏观上不论公转方向如何，随偏心度增大，宽间隙处的合速度剖面整体增大，窄间隙处的合速度剖面整体减小，但当具体分析不同公转方向的环空流场时，又有各自截然不同的流动特性。（2）正向公转时，钻柱的自转和公转共同增强了宽窄间隙处的正向旋流趋势，此时切向速度整体为正值，无二次流出现，偏心度增大，宽间隙处切向速度增大。（3）反向公转时，切向速度同时出现正反向旋流趋势，钻柱的公转流场与自转流场出现相互抵消状态，环空截面呈现不对称的二次流。偏心度减小，宽间隙处切向速度反向越小，公转流场与自转流场相互抵消越小，二次流趋势越明显，摩阻压耗越大。

作者：张晋凯 李根生 黄中伟 田守嶒 宋先知 王海柱 单位：油气资源与探测国家重点实验室 中国石化石油工程技术研究院