国际化学工程杂志》上

PDF
国际化学工程杂志》上/2019年/文章

研究文章|开放获取

体积 2019年 |文章的ID 7068989 | https://doi.org/10.1155/2019/7068989

Rasel a .苏丹穆罕默德Azizur拉赫曼细哔叽·鲁世德,索拉博Zendehboudi, Vassilios Kelessidis, CFD分析的压力损失和沉积速度的水平环形”,国际化学工程杂志》上, 卷。2019年, 文章的ID7068989, 17 页面, 2019年 https://doi.org/10.1155/2019/7068989

CFD分析的压力损失和沉积速度的水平环形

学术编辑器:Evangelos Tsotsas
收到了 2018年10月31日
接受 2019年1月05
发表 2019年2月3日

文摘

估计的压力损失和沉积速度是至关重要的水力设计的环形钻洞的石油工业。本研究调查的影响,流体速度、流体类型、粒子大小、粒子浓度、钻柱转速和偏心压力损失和沉降条件下利用计算流体动力学(CFD)。古怪的钻杆在0 - 75%的范围不同,它绕自己的轴旋转在0 - 150 rpm。模拟钻孔的直径比是0.56。实验数据证实了当前CFD模型的有效性使用ANSYS 16.2开发平台。

1。介绍

在石油行业,预测摩擦压力损失和沉淀条件的运输钻井液在轮对钻井作业很重要。不准确的预测会导致昂贵的钻井问题。几个例子的问题是循环,好玩,堵塞,磨损,磨损和钻机功率选择不当。现有的实证模型变得不那么准确的涉及许多简化的假设。CFD模拟有助于减少这些假设通过使用基于物理的n - s方程模型流的流体动力学系统。当前的工作重点是开发一个全面的CFD模型,能够考虑到所有重要的钻井参数的影响,如流体速度、流体类型、粒子大小、粒子浓度、钻杆转速和钻杆偏心。

2。文献综述

在环空压力损失的估计更困难与管流相比,由于产生的复杂性在水力学复杂几何(1,2]。从实证的角度来看,这个问题通常是通过更换管直径的管流模型的有效直径环。大量的“有效直径”的定义提出了至今为止。然而,很难选择一个字段定义应用程序的开发和/或应用经验。比较预测压力损失是由多个定义Anifowoshe和Osisanya3]。其他问题,使压力损失的估计困难是偏心和钻孔内钻杆的旋转速度。许多研究已经完成在环形非牛顿流体的流动引入经验/分析模型,允许考虑这些影响(1,3- - - - - -10]。先前的研究结果表明,环形压力损失非牛顿幂律)流体流动在一个钻孔取决于钻杆旋转速度、流体性质、流动机制(层流过渡/紊流),直径比、偏心度和等效水力粗糙。

使用ANSYS等商用CFD软件包流利预测钻井液的环形交通压力损失是相对较新的方法。Sorgun和Ozbayoglu9]证明了CFD模型的更好的性能比现有的实证模型在预测摩擦压力损失。Sorgun [8]调查管偏心度对压力损失的影响,切向速度、轴向速度,和有效利用CFD粘度。Erge et al。6]介绍了CFD建模方法,适用于摩擦压力的损失估计在偏心环空内管旋转循环时产生幂律流体。然而,大多数这些CFD研究仅限于一个单相层流的动态水力光滑轮。

环状流含岩屑的钻井液。泥浆,没有足够详细地研究。领域的工作的例子环形泥浆流在引用(11- - - - - -14]。这些研究的重点是理解在环空泥浆流的流体动力学实时实验和生产基于数据分析的经验模型。最近,不同的研究人员15- - - - - -17)利用CFD在环形研究泥浆的运输。Ofei [15)检查携带液体的流变参数的影响速度的固体。Sorgun和Ulker16)而获得的压力损失的预测使用人工神经网络(ANN)和计算流体动力学。这两种方法产生了类似的结果。太阳et al。17]研究倾向的影响,转速和流量的分布固体浓度和摩擦压力损失。类似于单相环状流的作品,其中大部分CFD研究仅限于层流浆流条件。

3所示。方法

3.1。CFD模拟

在当前的工作中,环形泥浆流动的CFD仿真模型是使用ANSYS 16.2流利的开发平台。以前作品后(18- - - - - -20.),multi-fluid颗粒模型用于描述流体混合物的流动行为。细粒度的版本的欧拉模型选为多相模型(附录C)。这是因为高固相体积分数将被用于这项研究和细粒度的版本捕获的流体动力学高浓度料浆组成不同的晶粒大小。它允许多个单独的建模,但互动阶段。阶段可以是液体、气体或固体在几乎任何组合。欧拉治疗用于每个阶段,拉格朗日治疗用于与离散相模型。

多相流的描述贯穿连续包含相位的体积分数的概念,它代表了每个阶段占据的空间。每个阶段满足质量守恒定律和动量。平均守恒方程被修改的地方瞬时平衡的每个阶段(21)或通过使用混合理论方法(22]。详细描述可在附录一个

欧拉多相计算,态势的简单(PC-SIMPLE)算法用于压力速度耦合。PC-SIMPLE是多相流的简单算法的扩展(23,24]。阶段的速度耦合隔离的方式来解决。使用的块代数多重网格方案density-based解算器是用于解决一个矢量方程形成的所有阶段的速度组件同时[25]。然后,压力修正方程建立基于总量的连续性。压力和速度然后修正,以满足连续性约束。

确保稳定和收敛的迭代过程,二阶逆风离散化是用于动量方程,和第一逆风离散化采用体积分数、湍动能及其耗散。逆风是指票面价值来源于上游细胞数量,或“逆风”,相对于正常速度的方向。当一阶精度,在细胞数量面临取决于假设任何领域的cell-center值变量代表cell-average价值和保存在整个细胞;面对数量是相同的细胞数量。因此,面对值设置等于cell-center值选择一阶逆风时上游的细胞。相比之下,当二阶精度或二阶逆风,细胞数量在面临计算使用多维线性重建方法(26]。在这种方法中,高阶精度达到细胞面临通过泰勒级数展开的cell-centered解决方案单元形心。整个模拟过程如图1

3.2。湍流模型的选择

流体湍流量计算采用雷诺应力模型(27- - - - - -29日(附录D)。放弃各向同性涡粘性假设,RSM关闭Reynolds-averaged n - s方程通过求解运输雷诺应力方程,结合耗散率的方程(30.]。这里,需要五个额外的传输方程在二维流动和七个额外的传输方程必须解决3 d(请见附件B为进一步的细节)。湍流模型被选为当前工作,分析不同模型的相对性能。一个典型的例子的分析呈现在图2。在这个图中,%的差异指的是百分比区别实验测量压力损失和相应的CFD预测。在大多数情况下,当采用RSM差异小于10%。一个重要的结果是在表的列表1


参考 实验条件 %的差异[100(模拟实验−)/实验)
RSM 风场k- - - - - -ω 标准k- - - - - -ω 标准k- - - - - -ε

Kaushal et al。(2005)31日] 管方向:水平 8 15 17 17.48
管印度。:0.0549 m
管材的施工:不锈钢(密度8030 kg / m3光滑的墙)
液体:单相水(密度:998.2公斤/米3粘度:0.001003公斤/ m)
流体速度:2 m / s

Camci (2003) (32] 环方向:水平和同心 9 17 20. 23
环尺寸:0.0432米ID、0.123 m OD
施工材料:铝(光滑的墙)
液体:单相水
流体速度:0.2米/秒

Kelessidis et al。(2011)33] 环方向:水平和同心 8 13 16 22.2
环尺寸:0.04米ID、0.07 m OD
施工材料:有机玻璃(光滑的墙)
液体:单相水
流体速度:1.12米/秒

Skudarnov et al。(2004)34] 管方向:水平 −2 −2.8 −3 −3.2
管直径:0.023米
管材的施工:不锈钢(壁粗糙度:32µ米)
液体:两阶段固液浆
液体:水(密度:998.2公斤/米3粘度:0.001003公斤/ m)
固体:玻璃球体(双物种密度2490 kg / m3和4200公斤/立方米,50% 50%体积的混合物,粒子直径(d):140µ米,体积浓度(Cv):15%)
流体速度:1.724米/秒

Ozbelge和Beyaz (2001) [35] 环方向:垂直和同心 −6 14 15 −15.1
OD环尺寸:0.125米,0.025米的ID
施工材料:不锈钢(密度8030 kg / m3,光滑的墙)
液体:两阶段固液浆
液体:水(密度998.2公斤/立方米,粘度0.001003公斤/ m)
固体:长石(平均颗粒直径0.23毫米,平均密度2500 kg / m3,体积浓度,简历= 1.80%)
流体速度:0.135米/秒

Ozbelge和Beyaz (2001) [35] 环方向:垂直和同心 −4 8 13 16
OD环尺寸:0.125米,0.025米的ID
施工材料:不锈钢(密度8030 kg / m3,光滑的墙)
液体:两阶段固液浆
液体:水(密度998.2公斤/立方米,粘度0.001003公斤/ m)
固体:长石(平均颗粒直径0.23毫米,平均密度2500 kg / m3,体积浓度,简历= 1.80%)
流体速度:0.197米/秒

3.3。长度独立研究

流的长度域被认为是足够长的时间来达到充分发展流动。最小入口长度考虑流开发是50Dh,水力直径Dh= OD-ID [36,37]。独立的测试长度的一个例子是图所示3。仿真结果不依赖于从入口长度在3 m。

3.4。网孔分析法

环形截面的计算网格生成使用ANSYS流畅,啮合是最终确定的基础上适当的网格独立性检查。多层壁附近的通货膨胀被添加来自内部和外墙墙附近计算不同参数的特点,更准确地说。墙面和流体之间的剪切应力更高,这通货膨胀有助于创建墙附近的密集的啮合。计算网格的一个例子分布和网格独立性测试数据所示45。网格独立结果可以产生超过800000的节点数量。所有的结果提出了在当前工作得到使用大约900000个节点。

无因次壁距离的值(y+)检查在墙附近的细胞生成的融合要求的考虑y+为细胞相邻的墙上。的价值y+取决于壁剪切应力、流体密度、水力直径,和分子粘度: 在哪里 是细胞墙的距离中心; ,分子的粘度; ,流体密度;和 ,墙剪切应力。最终,y+取决于网格的分辨率和流雷诺数。默认标准壁函数通常适用如果第一个细胞中心附近墙上有一个y+值大于30 (38]。的最低要求(y+> 30)的价值y+在我们的研究保持45以上。

3.4.1。收敛速度分析

10的最优收敛速度−5被选为终止迭代。图6显示了一个示例的分析找出最优收敛速度在10−6-10年−4。仿真结果收敛值大于10时不同−5。但是,结果并没有改变的值小于10−5

3.5。CFD模型验证

作为初步的步骤,验证了CFD建模方法对数据在开放的文学。验证提出的一些例子如下。

3.6。单相流过环形

实验数据集Kelessidis et al。33]和Camci [32比较图7拟议的CFD模型。几何是取自Kelessidis et al。(2011),在内径(ID)是0.04米,外径(OD)是0.07米,长度是5米(水平同心轮)。流体为水和墙材料是树脂玻璃(动态水力光滑的墙,ε一个= 0)。在Camci(2003),内径是0.0432米,外径0.123米,长度是5米。墙材料是铝(光滑的墙)。

在图7图,代表对数刻度。双对数尺度的压力梯度增加率几乎是线性的。仿真结果的平均百分比误差Kelessidis et al。(2011)和Camci(2003)分别为9.88%和8.46%,分别,这表明一个非常好的协议(据估计,实验数据的误差为±10%)。

3.7。两相(固液)流过环形

压力梯度(Pa / m)的水沙泥浆流经垂直同心环形相比Ozbelge和Beyaz [35实验数据图8。CFD模拟,液相密度是水(9982公斤/米3和粘度0.001003公斤/ m)和固相为长石(平均颗粒直径0.23毫米,平均密度2500 kg / m3)。长度5米,外径0.125米,内径0.025米,入口速度范围0.0738 - -0.197 m / s,和整体泥浆体积浓度与0.23毫米粒径范围1.0% - -1.8% (dp)被认为是边界条件。光滑的不锈钢管(密度8030 kg / m3)用于模拟。管子被认为是垂直的。,gravity effect is included, and gravity acceleration is directed opposite to the outlet. No slip condition for liquid and solid phases is used at the walls. Figure9显示了模拟和实验的比较两相摩擦压降通过垂直轮在不同速度和不同体积浓度的混合浆为0.23毫米的意思是砂粒径(dp)。模拟结果与实验值吻合较好,平均误差为2.62%。

4所示。结果与讨论

后达到良好的验证模型与实验数据,参数分析了观察流量变化的影响,钻杆旋转,离心率和颗粒大小(附录中给出详细的数据表(一个B))。中给出的参数用于分析表2


参数 价值/范围 单位

外直径 0.1143
内径 0.0635
长度 5
壁粗糙度 0 毫米
砂浓度 5 - 20 %
砂粒径 0.1 - -0.005 毫米
内管旋转 0 - 150 rpm
内管偏心 0 - 75 %
水当量雷诺数 - - - - - - 单位少

4.1。流量的影响

流体流动床上最大浓度的影响分析了数据10- - - - - -13。水与沙粒子混合物(浆)作为工作液。四个不同的条件考虑与解决砂进气浓度(20%)和沙颗粒大小(0.1毫米)。条件如下:(我)同心环形固定内管(2)与150 rpm同心环形旋转内管(3)50%偏心环形固定内管(iv)50%偏心环形150 rpm内管

每个案例的分析,可以看出床底部墙附近的浓度增加时,流量减少。由于重力和水平方向,砂颗粒倾向于聚集下盘附近并创建一个流堵塞。减少流动提高了过程。下面的分析,发现 雷诺数床上下盘附近浓度超过25%在所有情况下,当它下降 ,比例是50%以上。在图14,砂颗粒样本轮廓分布在不同的流率。操作条件从图10。等高线分布的环形截面3米距离入口显示流体流动和重力对浓度分布的影响。

4.2。钻杆旋转和偏心的影响

内管轮可以影响压力损失通过改变其偏心和旋转。内管旋转和偏心率的影响提出了数字915。在这个分析单相水用作液体。从图15它可以可视化,压力损失随转速的增加,这一趋势是相同的,不同的流率。静止不动的,50 rpm, 100 rpm, 150 rpm的旋转内管是考虑。增加的趋势在高流率高。由于转速粒子与粒子之间的和particle-wall碰撞和反弹上升,从而导致更高的压力损失。

9显示了内管偏心对压力损失的影响。同心,25%的偏心,偏心度分析了离心率50%,75%这个数字。在一个固定的流量(Re = 50000),与内管偏心距的增加,压力损失减少。这一趋势适用期间内管固定条件。与内管旋转的趋势是相反的。因为额外的碰撞增加管旋转,这种变化发生。

4.3。颗粒大小的影响

在泥浆流中,砂粒径粒子堵塞底部附近的墙上有一个有效的作用水平环形管。图16分析砂粒径的影响在同心和固定内管进口砂浓度为5%。与砂粒径的增加,粒子沉积下盘附近增加因为单个粒子重量增加,最终导致堵塞。从之前的分析中,我们发现最大床浓度随流量的增加而减小(数字8,10- - - - - -13),但对于粒径0.005毫米,最大床浓度下盘附近几乎是常数在不同的流率。这意味着更小的粒度(< 0.01毫米),粒子沉积的流率的影响可以忽略不计。

注意,这些粒子大小是选择从粒度分布(PSD)图表在CFD分析考虑。PSD的一个图表,如图17与平均粒径0.1毫米,0.1毫米粒度从八个不同的大小选择基于累积每个粒子的重量(%)。为了简化参数分析过程,只显示平均粒子尺寸。

5。结论

为了开发一个被广泛接受的CFD模型,通过钻井多相流轮,目前的工作涉及项目的总体规划和最初的进步。总之,本研究可以叙述如下:(我)CFD建模方法预测摩擦压力损失和沉淀条件进行了验证。通过示例验证了展示其在复杂钻井条件下的适用性。(2)以下重要的钻井参数对压力损失的影响和结算条件测试:流体流量、转速、偏心钻杆和固体粒子的大小。(3)流体流量的最大影响粒子沉积。在所有分析条件下,与流体流动的减少,粒子沉积在底部附近墙上增加。特定的沉积速度取决于特定的需求最大壁附近沉积,可以使用我们的方法计算。(iv)随着转速的增加和内管偏心,流体流动的能量损失(压力损失)增加。然而,在静止的条件下,压力损失随偏心率的增加而减小。(v)提出了当前研究项目的初步结果与牛顿流体流动。项目预计将产生一个综合CFD模型能够考虑到所有重要的钻井参数与牛顿和非牛顿流体。与非牛顿流体流动分析正在进行中。

附录

答:参数与单相流体研究图表

参数研究的结果图给出了单相流体在表3


控制变量 CFD结果
雷诺数(Re =dH ) 转速(rpm) 离心率(%) 压力梯度(Pa /米)

20000年 0 0 82年
40000年 0 0 150年
50000年 0 0 230年
60000年 0 0 320年
80000年 0 0 557.575
100000年 0 0 920.161
20000年 50 0 103年
40000年 50 0 158年
50000年 50 0 234年
60000年 50 0 324年
80000年 50 0 562.112
100000年 50 0 923.5
20000年 One hundred. 0 165年
40000年 One hundred. 0 209年
50000年 One hundred. 0 269年
60000年 One hundred. 0 345年
80000年 One hundred. 0 595.795
100000年 One hundred. 0 933.455
20000年 150年 0 275年
40000年 150年 0 293年
50000年 150年 0 343年
60000年 150年 0 407年
80000年 150年 0 645.529
100000年 150年 0 966.1
20000年 0 25 77年
40000年 0 25 140年
50000年 0 25 214年
60000年 0 25 299年
80000年 0 25 555.578
100000年 0 25 902.474
20000年 50 25 97年
40000年 50 25 150年
50000年 50 25 221年
60000年 50 25 304年
80000年 50 25 559.412
100000年 50 25 906.422
20000年 One hundred. 25 169年
40000年 One hundred. 25 206年
50000年 One hundred. 25 259年
60000年 One hundred. 25 327年
80000年 One hundred. 25 542.354
100000年 One hundred. 25 916.803
20000年 150年 25 291年
40000年 150年 25 298年
50000年 150年 25 340年
60000年 150年 25 395年
80000年 150年 25 582.623
100000年 150年 25 951.93
20000年 0 50 80年
40000年 0 50 150年
50000年 0 50 238年
60000年 0 50 342年
80000年 0 50 639.231
100000年 0 50 955.782
20000年 50 50 112年
40000年 50 50 178年
50000年 50 50 263年
60000年 50 50 364年
80000年 50 50 644.75
100000年 50 50 965.534
20000年 One hundred. 50 203年
40000年 One hundred. 50 251年
50000年 One hundred. 50 323年
60000年 One hundred. 50 416年
80000年 One hundred. 50 689.054
100000年 One hundred. 50 1003.29
20000年 150年 50 355年
40000年 150年 50 373年
50000年 150年 50 428年
60000年 150年 50 505年
80000年 150年 50 750.208
100000年 150年 50 1056.07
20000年 0 75年 71年
40000年 0 75年 135年
50000年 0 75年 215年
60000年 0 75年 310年
80000年 0 75年 574.076
100000年 0 75年 907.722
20000年 50 75年 123年
40000年 50 75年 197年
50000年 50 75年 288年
60000年 50 75年 395年
80000年 50 75年 680.247
100000年 50 75年 1013.64
20000年 One hundred. 75年 231年
40000年 One hundred. 75年 296年
50000年 One hundred. 75年 388年
60000年 One hundred. 75年 501年
80000年 One hundred. 75年 810.132
100000年 One hundred. 75年 1166.54
20000年 150年 75年 401年
40000年 150年 75年 449年
50000年 150年 75年 528年
60000年 150年 75年 631年
80000年 150年 75年 948.062
100000年 150年 75年 1313.03

b参数与两相流体研究图表

参数研究的结果图给出了两相流体在桌子上4


控制变量 CFD结果
雷诺数(Re =dH ) 转速(rpm) 离心率(%) 粒子输入浓度(%) 颗粒大小(毫米) 压力梯度(Pa /米) 下盘附近最大床浓度(%)

40000年 静止的 同心 5 0.05 766.966 58.12
70000年 静止的 同心 5 0.05 848.724 22.09
100000年 静止的 同心 5 0.05 1120.29 12.38
200000年 静止的 同心 5 0.05 3126.42 7.49
40000年 静止的 同心 20. 0.05 767.134 58.14
70000年 静止的 同心 20. 0.05 848.557 22.04
100000年 静止的 同心 20. 0.05 1120.29 12.36
200000年 静止的 同心 20. 0.05 3118.21 7.47
40000年 静止的 同心 5 0.1 767.674 58
70000年 静止的 同心 5 0.1 837.6 42.1
100000年 静止的 同心 5 0.1 1112.26 25.91
200000年 静止的 同心 5 0.1 3120.31 10.15
40000年 静止的 同心 20. 0.1 761.902 52.73
70000年 静止的 同心 20. 0.1 837.166 42.15
100000年 静止的 同心 20. 0.1 1111.41 26.02
200000年 静止的 同心 20. 0.1 3129.55 10.16
40000年 静止的 同心 5 0.005 759.968 5.27
70000年 静止的 同心 5 0.005 848.862 5.14
100000年 静止的 同心 5 0.005 1121.32 5.1
200000年 静止的 同心 5 0.005 3143.08 5.06
40000年 静止的 同心 20. 0.005 759.924 5.27
70000年 静止的 同心 20. 0.005 848.871 5.14
100000年 静止的 同心 20. 0.005 1121.6 5.1
200000年 静止的 同心 20. 0.005 3119.51 5.06
40000年 150年 同心 5 0.05 812.232 57.79
70000年 150年 同心 5 0.05 888.334 21.84
100000年 150年 同心 5 0.05 1162.98 12.33
200000年 150年 同心 5 0.05 3151.97 7.48
40000年 150年 同心 20. 0.05 812.265 57.74
70000年 150年 同心 20. 0.05 888.189 21.65
100000年 150年 同心 20. 0.05 1162.52 12.32
200000年 150年 同心 20. 0.05 3151.46 7.47
40000年 150年 同心 5 0.1 808.946 53.66
70000年 150年 同心 5 0.1 886.414 41.42
100000年 150年 同心 5 0.1 1155.53 25.81
200000年 150年 同心 5 0.1 3160.58 10.17
40000年 150年 同心 20. 0.1 810.85 51.43
70000年 150年 同心 20. 0.1 886.051 41.47
100000年 150年 同心 20. 0.1 1156.06 25.85
200000年 150年 同心 20. 0.1 3149.8 10.16
40000年 150年 同心 5 0.005 799.978 5.27
70000年 150年 同心 5 0.005 882.585 5.14
100000年 150年 同心 5 0.005 1162.34 5.1
200000年 150年 同心 5 0.005 3166.04 5.06
40000年 150年 同心 20. 0.005 799.909 5.27
70000年 150年 同心 20. 0.005 882.953 5.14
100000年 150年 同心 20. 0.005 1162.38 5.1
200000年 150年 同心 20. 0.005 3146.69 5.06
85000年 静止的 50 5 0.05 1042.12 33.71
100000年 静止的 50 5 0.05 1198.23 26.03
200000年 静止的 50 5 0.05 3255.59 9.96
300000年 静止的 50 5 0.05 6848.14 7.76
85000年 静止的 50 20. 0.05 1055.6 33.46
100000年 静止的 50 20. 0.05 1206.28 24.9
200000年 静止的 50 20. 0.05 3278.31 9.9
300000年 静止的 50 20. 0.05 6693.65 7.78
85000年 静止的 50 5 0.1 1011.03 26.39
100000年 静止的 50 5 0.1 1263.45 23.56
200000年 静止的 50 5 0.1 3170.89 9.91
300000年 静止的 50 5 0.1 6479.87 6.83
85000年 静止的 50 20. 0.1 1016.64 30.3
100000年 静止的 50 20. 0.1 1179.16 23.6
200000年 静止的 50 20. 0.1 3180.33 10.62
300000年 静止的 50 20. 0.1 6636.22 7.2
85000年 静止的 50 5 0.005 1036.05 5.35
100000年 静止的 50 5 0.005 1192.55 5.3
200000年 静止的 50 5 0.005 3184.09 5.15
300000年 静止的 50 5 0.005 6544.09 5.1
85000年 静止的 50 20. 0.005 1035.73 5.36
100000年 静止的 50 20. 0.005 1195.77 5.3
200000年 静止的 50 20. 0.005 3186.63 5.15
300000年 静止的 50 20. 0.005 6519.47 5.1
85000年 150年 50 5 0.05 1095.77 32.5
100000年 150年 50 5 0.05 1259.55 24.06
200000年 150年 50 5 0.05 3342.23 10.21
300000年 150年 50 5 0.05 6794.4 7.72
85000年 150年 50 20. 0.05 1096.93 29.37
100000年 150年 50 20. 0.05 1277.49 22.62
200000年 150年 50 20. 0.05 3352.92 9.85
300000年 150年 50 20. 0.05 6849.01 7.55
85000年 150年 50 5 0.1 1075.24 25.28
100000年 150年 50 5 0.1 1226.04 24.57
200000年 150年 50 5 0.1 3196.13 9.8
300000年 150年 50 5 0.1 6533.21 6.94
85000年 150年 50 20. 0.1 1071.57 25.25
100000年 150年 50 20. 0.1 1244.69 22.86
200000年 150年 50 20. 0.1 3216.8 10.19
300000年 150年 50 20. 0.1 6656.04 7.32
85000年 150年 50 5 0.005 1092.88 5.35
100000年 150年 50 5 0.005 1257.26 5.29
200000年 150年 50 5 0.005 3246.63 5.15
300000年 150年 50 5 0.005 6567.77 5.1
85000年 150年 50 20. 0.005 1091.97 5.35
100000年 150年 50 20. 0.005 1261.21 5.29
200000年 150年 50 20. 0.005 3243.02 5.15
300000年 150年 50 20. 0.005 6561.63 5.1

c .多相模型的描述

相的体积 , ,被定义为 在哪里

的有效密度的阶段 在哪里 阶段的物理密度吗

油和颗粒多相流动的方程是这里介绍的一般情况n步流。每个阶段的体积分数从连续性方程计算如下。 在哪里 相位参考密度或体积平均密度的 阶段的解决方案域, 特征的传质 阶段, 特征的传质 阶段。

流体相的动量守恒 在哪里 重力加速度, 阶段应力-应变张量, 是一个外部的身体力量, 是一个升力, 是一个虚拟质量力。

固相的动量守恒 在哪里 固体压力, 是液体或固相之间的动量交换系数 和固相 , 阶段的总数。

d .湍流模型的描述

确切的雷诺应力输运方程( )如下: 在哪里 率变化的总和吗 和运输 通过对流。 生产的速度吗 的扩散运输吗 耗散率。 是压力应变相关术语。 是旋转项。

模拟中使用的扩散项如下: 在哪里 ,

产率 可以表示为

pressure-strain项 或减缓pressure-strain项也被称为return-to-isotropy项, 或快速pressure-strain项, wall-reflection术语。它可以表示为 在哪里 在哪里

wall-reflection项 负责正常的墙附近的应力分布。这一项是建模为 在哪里 组件的水平分量是正常的在墙上, 是最短的距离在墙上, ,在哪里 卡门常数(= 0.4187)。

耗散率 或破坏的速度 被建模为 在哪里 =波动变形率。

旋转术语表示为 在哪里 =旋转矢量。 =交替的象征,+ 1−1或0根据、j和k

命名法

CFD: 计算流体动力学
RSM: 雷诺应力模型
风场: 剪切应力运输
跑: 雷诺平均n - s
简单: Semi-implicit方法有关的压力方程
d: 砂粒径
Cv: 砂体积浓度
ID: 管内径
OD: 管外径
DH: 水力直径
v: 流体速度
: 水速度
: 砂速度
: 壁粗糙度
: 雷诺数
: 摩擦系数。

数据可用性

ANSYS模拟的数据用于支持本研究的发现可以从相应的作者。

信息披露

会议论文(39)由同一作者提出了与主数据集的更新和修改。

的利益冲突

作者还宣布没有利益冲突有关的出版。

确认

作者感谢工学院和纽芬兰纪念大学的应用科学,卡塔尔国家研究基金(格兰特nprp10 - 0101 - 170091),德克萨斯A&M大学卡塔尔,和ANSYS Inc .)帮助和资助这项研究工作。

引用

  1. a . Pilehvari和r . Serth广义轴流的水力计算方法非牛顿流体在偏心环形,”SPE钻和完成,24卷,不。4、553 - 563年,2013页。视图:出版商的网站|谷歌学术搜索
  2. t·里德和A . Pilehvari”新模式层流、过渡和湍流的钻井泥浆,”学报SPE生产业务研讨会石油工程师学会,俄克拉何马州,美国,1993年。视图:谷歌学术搜索
  3. o . l . Anifowoshe和s . o . Osisanya当量直径的影响定义在摩擦压力损失估计管旋转的环”学报SPE深水钻井和完成会议石油工程师学会,加尔维斯顿TX,美国,2012年。视图:谷歌学术搜索
  4. r·m·艾哈迈德·m . s . Enfis h . m . El Kheir m . Laget和a . Saasen”钻旋转等效循环密度的影响:实地测量的建模和分析”SPE学报》年度技术会议和展览石油工程师学会,佛罗伦萨,意大利,2010年。视图:谷歌学术搜索
  5. i t Dosunmu和s . n .沙”,评价摩擦系数相关性和当量直径定义为管和环形非牛顿流体的流动,”石油科学与工程》杂志上卷,109年,第86 - 80页,2013年。视图:出版商的网站|谷歌学术搜索
  6. o . Erge e . m . Ozbayoglu s . z Miska et al .,“等效循环密度收益率幂律流体建模验证使用CFD方法,”石油科学与工程》杂志上卷。140年,16-27,2016页。视图:出版商的网站|谷歌学术搜索
  7. s a·汉森和n . Sterri”钻杆旋转对苗条的环形摩擦压力损失的影响,”学报SPE年度技术会议和展览,石油工程师学会,1995年10月。视图:谷歌学术搜索
  8. m . Sorgun”计算流体动力学建模的管偏心影响牛顿和非牛顿流体的流动特征,“能源,一个部分:复苏,利用率和环境影响,33卷,不。12日,第1208 - 1196页,2011年。视图:出版商的网站|谷歌学术搜索
  9. m . Sorgun和m . Ozbayoglu”预测摩擦压力损失为非牛顿流体在水平钻井,”能源部分复苏利用率和环境影响,33卷,不。11日,第1106 - 1104页,2011年。视图:出版商的网站|谷歌学术搜索
  10. r .萨勃拉曼尼亚和j .阿扎尔的实验研究对摩擦压降nonnewtonian钻井液管和环状流”学报》国际石油和天然气在中国会展,石油工程师学会,北京,中国,2000。视图:谷歌学术搜索
  11. g . Camci和t . a . Ozbelge稀释测定泥浆密度在垂直环用等速采样,”化学工程通信,卷193,不。11日,第1501 - 1482页,2007年。视图:出版商的网站|谷歌学术搜索
  12. s . k . Dewangan和s . l . Sinha”,探索水平井的井眼净化参数使用两阶段欧拉CFD方法,”《计算多相流,8卷,不。1、15 - 39,2016页。视图:出版商的网站|谷歌学术搜索
  13. t . A . Ozbelge和A . n . Eraslan”,计算流体动力学和传热研究湍流稀释的过程通过同心环空泥浆,”土耳其工程和环境科学》杂志上,30卷,不。1,1-13,2006页。视图:谷歌学术搜索
  14. t . A . Ozbelge和g·c·Unal”,一种新的两阶段相关压力下降充分发展稀浆过程通过一个环,“化学工程通信,卷196,不。4、491 - 498年,2008页。视图:出版商的网站|谷歌学术搜索
  15. t . n . Ofei”效应的收益率幂律流体流变特性在偏心岩屑运输水平狭窄环,“杂志的液体,2016年,页1 - 10。视图:出版商的网站|谷歌学术搜索
  16. m . Sorgun和大肠Ulker建模和实验研究固液两相压降在水平井筒,管旋转,“能源技术交易的ASME杂志》上,卷138,不。2、2016。视图:谷歌学术搜索
  17. t . x的太阳,k . Wang Yan, s .邵和j .焦”钻杆旋转对岩屑的影响运输使用计算流体动力学(CFD)在复杂结构井,“石油勘探和生产技术杂志》上,4卷,不。3、255 - 261年,2014页。视图:出版商的网站|谷歌学术搜索
  18. b . j .桤木和t·e·温赖特在分子动力学研究。二世。少量的弹性球体的行为。”化学物理学报,33卷,不。5,1439 - 1451年,1960页。视图:出版商的网站|谷歌学术搜索
  19. 查普曼和t . g .整流罩非均匀气体的数学理论:一个帐户的粘度、分子运动论热传导和扩散的气体英国剑桥,剑桥大学出版社,1970年。
  20. m . Syamlal w·罗杰斯,t . j . O ' brienMFIX文档:理论指导国家能源技术实验室,能源部,1993年,Technical Note DOE / METC-95/1013和NTIS / DE95000031。
  21. t·b·安德森和r·杰克逊,“流化床的流体力学描述。运动方程。”工业化学与工程化学基础》第六卷,没有。4、527 - 539年,1967页。视图:出版商的网站|谷歌学术搜索
  22. r·m·鲍恩“混合物理论”,连续物理学,3卷,不。第一部分,1976年。视图:谷歌学术搜索
  23. s . Patankar数值传热和流体流动美国佛罗里达州波卡拉顿,CRC新闻,1980年。
  24. 美国Vasquez诉伊万诺夫,“相耦合的方法解决多相非结构化网格上的问题,”ASME 2000流体工程部门夏季会议论文集(ASME FEDSM ' 00),纸fedsm2000 - 11281美国,页1 - 6,波士顿,MA, 2000。视图:谷歌学术搜索
  25. j·m·维斯,j.p. Maruszewski和w·a·史密斯,”隐式解决方案上的n - s方程使用代数多重网格,”张仁杂志,37卷,不。1,29-36,1999页。视图:出版商的网站|谷歌学术搜索
  26. t·巴斯和d . Jespersen逆风方案的设计和应用程序在非结构化网格,”27日航空学报》美国航空航天学会科学会议1989年美国NV,雷诺。视图:谷歌学术搜索
  27. m·m·吉布森和b . e .槽”,地面大气边界层对压力波动的影响,“流体力学杂志,卷86,不。3、491 - 511年,2006页。视图:出版商的网站|谷歌学术搜索
  28. b . e .槽、g·j·莉丝和w·Rodi“进步发展的雷诺应力湍流闭合,”流体力学杂志,卷68,不。03年,537 - 566年,2006页。视图:出版商的网站|谷歌学术搜索
  29. b . e .槽”,二次矩闭包和它的使用建模的工业流,”国际期刊液体中的数值方法,9卷,不。8,963 - 985年,1989页。视图:出版商的网站|谷歌学术搜索
  30. a .科林j . Chorin navier - stokes方程的数值解,“数学的计算,22卷,不。104年,745页,1968。视图:出版商的网站|谷歌学术搜索
  31. 佐藤d·r·Kaushal k, t .丰田,k . Funatsu y获利,”粒度分布对压降的影响和集中在管道的高度浓缩浆流,”国际多相流杂志》上没有,卷。31日。7,809 - 823年,2005页。视图:出版商的网站|谷歌学术搜索
  32. g . Camcı“等速采样技术的应用为当地固体密度上升固液流经一个环,“METU, 2003,博士论文。视图:谷歌学术搜索
  33. v . c . Kelessidis p Dalamarinis, r . Maglione”实验研究和预测的压力损失的液体建模为Herschel-Bulkley同心和偏心环形层流、过渡和湍流流动,”石油科学与工程》杂志上,卷77,不。3、305 - 312年,2011页。视图:出版商的网站|谷歌学术搜索
  34. p . v . Skudarnov、林c . x和m . a . Ebadian“Double-species浆流在水平管道,”流体工程,卷126,不。1,第132 - 125页,2004。视图:出版商的网站|谷歌学术搜索
  35. t . a . Ozbelge和a . Beyaz稀释固液向上流经垂直环在闭环系统中,“国际多相流杂志》上,27卷,不。4、737 - 752年,2001页。视图:出版商的网站|谷歌学术搜索
  36. n·p·布朗和n . i .海伍德泥浆处理:固液系统的设计施普林格科学与商业媒体,柏林,德国,1991年。
  37. e . j .黄蜂j.p.肯尼,r·l·甘地“固体——液体流:煤浆管道运输。(泵、阀门、机械设备、经济学),“散装材料处理,1卷,不。4、1977。视图:谷歌学术搜索
  38. 有限元分析软件,12.0/12.1的文档。用户指南手册Ansys Inc . Canonsburg, PA,美国,2009年。
  39. s·鲁世德r·a·苏丹a·拉赫曼和诉Kelessidis”在偏心斜轮调查的压力损失ASME学报》2017年第36届国际海洋会议,离岸和北极工程美国机械工程师协会,特隆赫姆挪威,2017。视图:谷歌学术搜索

版权©2019 Rasel苏丹等。这是一个开放的分布式下文章知识共享归属许可,它允许无限制的使用、分配和复制在任何媒介,提供最初的工作是正确引用。


更多相关文章

PDF 下载引用 引用
下载其他格式更多的
订单打印副本订单
的观点5468年
下载1253年
引用

相关文章

文章奖:2020年杰出的研究贡献,选择由我们的首席编辑。获奖的文章阅读