
Citation: | ZHANG Xuemin, ZHANG Xueru, LI Houbu, et al. Simulation study on mechanical behavior of a nonmetallic composite coiled tubing with cable laying under tension load [J]. Petroleum Drilling Techniques, 2025, 53(1):94−101. DOI: 10.11911/syztjs.2025003 |
Nonmetallic composite coiled tubing with cable laying is subjected to tension load due to its self-weight in the process of frequent lifting and lowering of wells for oil extraction, and clarifying the mechanical behavior of the pipe under this load can provide guidance for the safe service of the pipe. A three-dimensional numerical model of nonmetallic composite coiled tubing with cable laying was constructed by finite element software, and the mechanical behavior of the pipe under tension load and the mechanical response of each structural layer were analyzed. The influence of cable laying process parameters, such as cable winding and distribution angles, on the mechanical properties of the pipe was explored. The results indicate that under tension load, the stresses in all structural layers of the pipe with cable laying exhibit a spiral distribution because of cable winding. When the pipe is stretched to failure, it undergoes three stages: elastic deformation, transition stage, and yield deformation. Meanwhile, the cables are in a state of small plastic uniform deformation. Reducing the cable winding angle can enhance the elastic modulus and axial load-bearing capacity of the pipe with cable laying. However, it may cause the pipe to enter into the transition stage prematurely and then yield in advance. The cable distribution angle has a minimal impact on the mechanical properties of the pipe. Therefore, when this type of pipe is manufactured, special emphasis should be placed on the cable winding angle, as this parameter is highly correlated with the mechanical properties of the pipe under tension loads.
The flow of non-Newtonian flow through annular section has received much attention of researchers to provide optimized solutions for industry challenges. Mud displacement in the annular section between drilling string and casing is one of its industry applications. The drilling fluid is used to carry out different functions, such as maintaining borehole stability, transporting drilled cuttings to the surface and cooling the bit[1], and preventing costly problems when proper hydraulics design of the well is needed during the process. Among models that correctly predict the rheological behavior of drilling muds, Herschel-Bulkley, also known as yield power-law, fits for a large range of shear rates[2]. C. Nouar et al[3], carried out a numerical and experimental study of Herschel-Bulkley fluid flow through annulus with rotation of the inner pipe, they found that rotation of the inner pipe causes a decrease of the axial velocity gradient
In the present paper, the influence of rheological parameters (yield stress τ0, consistency index K and behavior index n) and inner pipe rotation on the distribution of the axial and tangential velocity profile of Herschel-Bulkley fluid in the wide and narrow regions of eccentric annulus are studied. Effects of these parameters on pressure drop gradient are evaluated utilizing a CFD commercial code ANSYS Fluent 17.0 as well.
To study the impact of the inner pipe rotation and rheological parameters (yield stress τ0, consistency index K and behavior index n) on the axial and tangential velocity profiles of Herschel-Bulkley fluid in annular section, the flow is assumed to be fully developed, incompressible, steady and isothermal in laminar regime.
The continuity equation is expressed as:
∇⋅v=0 |
(1) |
where v is flow rate.
The momentum equation governing the flow in annulus is expressed in terms of τ in cylindrical coordinates
ρ(∂vr∂t+vr∂vr∂r+vθr∂vr∂θ+vz∂vr∂z−v2θr)=−∂p∂r−[1r∂(rτrr)∂r+1r∂τθr∂θ+∂τzr∂z−τθθr]+ρgr |
(2) |
ρ(∂vθ∂t+vr∂vθ∂r+vθr∂vθ∂θ+vz∂vθ∂z+vrvθr)=−1r∂p∂θ−[1r2∂(r2τrθ)∂r+1r∂τθθ∂θ+∂τzθ∂z+τθr−τrθr]+ρgθ |
(3) |
ρ(∂vz∂t+vr∂vz∂r+vθr∂vz∂θ+vz∂vz∂z)=−∂p∂θ−[1r∂(rτrz)∂r+1r∂τθz∂θ+∂τzz∂z]+ρgz |
(4) |
where
The rheological behavior of the non-Newtonian fluid is modeled as:
τ=τ0+K(˙γ)n |
(5) |
where
In the present study, Herschel-Bulkley fluid flows through eccentric annulus (E=0.5) to simulate mud pattern during the drilling of horizontal well, in which the eccentricity is defined as:
E=2eDo−Di |
(6) |
where E is eccentricity of the inner cylinder, e is the distance between the centers of the inner and outer pipes(m), Do is diameter of the outer cylinder(m), Di is diameter of the inner cylinder(m).
To predict the flow regime of yield power-law fluid in the annulus, Reynolds number is calculated using the relationship presented by K. Madlener et al[11]:
ReYPL=ρDnhu2−nτ08(Dhu)n+K(3m+14m)n8n−1 |
(7) |
where ReYPL is the Reynolds number of yield power-law fluid, u is the bulk flow velocity(m/s), m is the local shear rate(s–1), Dh is the hydraulic diameter(m) which are calculated as:
Dh=Do−Di |
(8) |
m=nK(8uDh)nτ0+K(8uDh)n |
(9) |
To prevent the entrance effects, the length of the cylinders is selected to be longer than hydrodynamic entry length, which is given as[12]:
Lh,laminar=0.05(Do−Di)Re |
(10) |
where Lh, laminar is length of the hydrodynamic entry.
The rheological parameters of Herschel-Bulkley fluid and geometry characteristics are:
The domain of the fluid flow is meshed into 800 000 elements (20 radial divisions, 80 circumferential divisions and 500 axial divisions), where the number of elements is selected to ensure the results independence of mesh model used as well as keeping the number of elements as low as possible to save computational time, as shown in the Fig. 1 and Fig. 2. The commercial code ANSYS Fluent 17.0 based on finite volume method is used to solve the differential equations, where the flow equations are integrated over each control volume. The solution is assumed to be converged, when the convergence criteria 10–4 is reached.
Results obtained from the numerical model are compared to those of experimental work of C. Nouar et al[3] are shown in Fig. 3. The Simulation and experimental results are in good agreement in which the mean error are of 7.7% and 6.4% for dimensionless axial and tangential velocity, respectively.
Comparison of the experimental and numerical results reveals the ability of numerical simulation to provide accurate results.
Fig. 4 exhibits that the increase of the inner pipe rotation from 100 to 400 r/min causes an increase of 120 % of the maximum axial velocity in the narrow region of the annulus, which could improve the well cleaning process by the transportation of cuttings from the lower side of the horizontal annulus preventing the formation of cuttings bed. However, the increase of the inner pipe rotation, induces a reduction of the velocity gradient (∂u/∂r) near the outer pipe in the wide region of the annulus. Similar to the conclusion reported by C. Nouar et al[3].
Fig. 5 shows that the increase of yield stress from 8 to 32 Pa causes a decrease of 6.1 % and 3.6 % of the maximum axial velocity in the wide and narrow regions, respectively.
When the consistency index diminishes from 16 to 4 Pa·sn, the axial velocity profile becomes more flat in the wide region of the annulus, however, a negligible effect is observed for the narrow region, as exhibited in the Fig. 6.
As the behavior index n decreases from 0.8 to 0.2, the axial velocity profile becomes flat in the wide region of the eccentric annulus, as shown in the Fig. 7. While a negligible effect is observed in the narrow region except for n=0.8, where the negative values of the axial velocity profile could be attributed the rotation of the inner pipe which probably increase pressure drop in the annulus.
For all range of the inner pipe rotation, the tangential velocity decreases dramatically near the inner pipe rotation in the wide region of the annulus. As the fluid moved away from the inner pipe, this decrease becomes gradual except for ω=100 r/min, the decrease becomes more gradual in the center of annulus. After that, the tangential velocity begins to decease sharply again until the outer pipe, as presented in the Fig. 4.
As the yield stress diminishes, the tangential velocity decreases more gradually in the center of annulus in the wide region, while a similar variation is observed near the outer and inner pipes, as showed in the Fig. 5.
It can be seen from Fig. 7 that tangential velocity profile in the center of annulus of the wide region decreases more gradually as the behavior index decreases from 0.8 to 0.2 where the tangential velocity profile presents negative values from center of the annulus until the outer pipe. The negative values of the tangential velocity explained by the presence of a secondary flow (also called counter rotating swirl) which rotates in the opposite direction of the inner pipe. The appearance of the secondary flow could decrease carrying capacity of the drilling mud, which affects the whole cleaning process.
In the narrow region of the annulus, the tangential velocity profile decreases dramatically as the fluid moves away from the inner pipe. Fig. 4 to Fig. 7 show that the variation of the inner pipe rotation and rheological parameters have a slight effect on the tangential velocity.
As can be seen in Fig. 8, the increase of the inner pipe rotation from 0 r/min to 400 r/min causes a decrease of pressure drop gradient of yield power-law fluid, this decrease is estimated around 10% for all eccentric annulus (E=0.2, 0.4, 0.6 and 0.8). Since shear thinning phenomenon tends to reduce the pressure drop of yield power-law fluid, inertial effects induced by rotation of the inner pipe are dominated by shear thinning phenomenon in the annulus. This trend is also depicted by R. M. Ahmed et al[7].
Fig. 9 depicts the effect of the yield stress of yield power-law fluid on pressure drop gradient, as shown in the figure, the increase of the yield stress from 8 Pa to 32 Pa induces an increment of 15% of pressure drop gradient for all eccentricities. This increase could be attributed to additional required stress on the fluid to initiate the flow which enhance inertial effects of yield power-law fluid in the annulus.
Fig.10 displays variation of pressure drop gradient of yield power-law fluid in the annulus when the flow consistency index increases from 4 Pa·sn to 16 Pa·sn. As seen in the figure, an increase of 150% of pressure drop gradient for all eccentricities caused by the increment of the flow consistency index. This considerable increase could be explained by the increase of the flow resistance in the annulus due to the increasing fluid viscosity.
Moreover, it can be stated that as the flow consistency index increases, the effect of the eccentricity of cylinders on pressure drop gradient becomes more pronounced.
Fig. 11 exhibits impact of the flow behavior index on pressure drop gradient of yield power-law fluid. As can be seen, an exponential increase of pressure drop gradient as the flow behavior index gets closer from Newtonian behavior where the fluid is less affected by shear thinning phenomenon, which makes inertial effects dominate the flow of yield power-law fluid in the annulus. It was also observed that for low values of the flow behavior index, the eccentricity of cylinders has a slight effect, however, as the flow behavior index gets greater, pressure drop gradient of yield power-law fluid decreases with the increase of the eccentricity.
For all rheological parameters and inner pipe rotation, pressure drop gradient of yield power-law fluid decreases with the increase of the eccentricity.
1) The increase of the inner pipe rotation from 100 to 400 r/min increases the axial velocity in the narrow region of the eccentric annulus, which could enhance the cleaning process in the lower part of the eccentric annulus. However, a slight effect is observed for the maximum axial velocity in the wide region of the annulus.
2) For the behavior index
3) Low value of the behavior index (
4) The variation of the inner pipe rotation and rheological parameters of the Herschel-Bulkley fluid have a negligible effect on the tangential velocity profile in the wide region for an annulus of
5) The increase of the inner pipe rotation from 0 rpm to 400 rpm causes a decrease of 10% of pressure drop gradient of yield power-law fluid for all eccentric annulus (E=0.2, 0.4, 0.6 and 0.8).
6) As the flow behavior index gets closer from the Newtonian behavior, the eccentricity of cylinders pressure drop gradient of yield power-law fluid decreases with the increase of the eccentricity.
[1] |
刘耀民,孙文盛,鲜林云. 敷缆连续管在油田电潜泵采油技术中的应用[J]. 焊管,2022,45(3):61–64.
LIU Yaomin, SUN Wensheng, XIAN Linyun. Application of coiled tubing with cables in oilfield electric submersible pump oil production technology[J]. Welded Pipe and Tube, 2022, 45(3): 61–64.
|
[2] |
李昂,杨万有,郑春峰,等. 海上油田采油技术创新实践及发展方向[J]. 石油钻探技术,2024,52(6):75–85.
LI Ang, YANG Wanyou, ZHENG Chunfeng, et al. Innovation practice and prospect of oil production technologies in offshore oil-fields[J]. Petroleum Drilling Techniques, 2024, 52(6): 75–85.
|
[3] |
李厚补,张学敏,马相阳,等. 井下用非金属复合材料连续管研究进展[J]. 石油管材与仪器,2021,7(2):9–14.
LI Houbu, ZHANG Xuemin, MA Xiangyang, et al. Research progress on non-metallic coiled composite pipes[J]. Petroleum Tubular Goods & Instruments, 2021, 7(2): 9–14.
|
[4] |
CHEN Wei, XIONG Haichao, BAI Yong. Failure behavior analysis of steel strip–reinforced flexible pipe under combined tension and internal pressure[J]. Journal of Thermoplastic Composite Materials, 2020, 33(6): 727–753. doi: 10.1177/0892705718811405
|
[5] |
祖世强,熊巍,张宝辉,等. 敷缆连续油管无杆采油试验及分析[J]. 内蒙古石油化工,2019,45(9):10–12.
ZU Shiqiang, XIONG Wei, ZHANG Baohui, et al. Test and analysis of rodless oil recovery of coiled tubing[J]. Inner Mongolia Petrochemical Industry, 2019, 45(9): 10–12.
|
[6] |
侯伟,刘传友,浮昀,等. 复合连续敷缆管采油技术探索实践[C]//第十七届宁夏青年科学家论坛石油石化专题论坛论文集. 银川:宁夏回族自治区科学技术协会,2021:197-199.
HOU Wei, LIU Chuanyou, FU Yun, et al. Exploration and practice of composite continuous cable pipe oil production technology[C]//Proceedings of the 17th Ningxia Young Scientists Forum Petroleum and Petrochemical Special Forum. Yinchuan: Ningxia Association for Science and Technology, 2021: 197-199.
|
[7] |
WANG Baodong, LIU Xiaoben, ZHANG Hong, et al. A combined experimental and numerical simulation approach for burst pressure analysis of fiber-reinforced thermoplastic pipes[J]. Ocean Engineering, 2021, 236: 109517. doi: 10.1016/j.oceaneng.2021.109517
|
[8] |
WANG Yangyang, LOU Min, TONG Bing, et al. Mechanical properties study of reinforced thermoplastic pipes under a tensile load[J]. China Ocean Engineering, 2020, 34(6): 806–816. doi: 10.1007/s13344-020-0073-x
|
[9] |
KAMP G P, BETTS M. Development of a power and data transmission thermoplastic composite coiled tubing for electric drilling[R]. SPE 60730, 2000.
|
[10] |
李帅,袁文才,张友军,等. 智能非金属敷缆连续管的夹持力学性能研究[J]. 化学工程与装备,2022(12):28–30.
LI Shuai, YUAN Wencai, ZHANG Youjun, et al. Study on the clamping mechanical properties of non-metallic intelligent coiled tubing[J]. Chemical Engineering & Equipment, 2022(12): 28–30.
|
[11] |
宿振国. 敷缆复合材料连续管结构设计与性能评测研究[D]. 青岛:中国石油大学(华东),2014.
SU Zhenguo. Study on structural design and performance evaluation of composite coiled tubing with cables[D]. Qingdao: China University of Petroleum(East China), 2014.
|
[12] |
丁楠,李厚补,古兴瑾,等. 非金属智能连续管拉伸层力学特性研究[J]. 石油机械,2020,48(11):119–125.
DING Nan, LI Houbu, GU Xingjin, et al. Study on the mechanical properties of the tensile layer of non-metallic intelligent coiled tubing[J]. China Petroleum Machinery, 2020, 48(11): 119–125.
|
[13] |
夏和萍. 热塑性复合材料海洋柔性管内压承载力研究[D]. 青岛:中国石油大学(华东),2021.
XIA Heping. On internal pressure capacity of thermoplastic composite pipes[D]. Qingdao: China University of Petroleum(East China), 2021.
|
[14] |
KNAPP R H. Derivation of a new stiffness matrix for helically armoured cables considering tension and torsion[J]. International Journal for Numerical Methods in Engineering, 1979, 14(4): 515–529. doi: 10.1002/nme.1620140405
|
[15] |
白勇,张尹,杨红刚. 纤维绳缠绕增强复合管拉伸有限元分析[J]. 低温建筑技术,2015,37(1):61–63.
BAI Yong, ZHANG Yin, YANG Honggang. Finite element analysis of fiber rope winding reinforced composite pipe tensile[J]. Low Temperature Architecture Technology, 2015, 37(1): 61–63.
|
[16] |
张学敏,周腾,李厚补,等. 涤纶纤维增强聚乙烯复合管承压性能模拟[J]. 工程塑料应用,2020,48(8):118–122.
ZHANG Xuemin, ZHOU Teng, LI Houbu, et al. Simulation of pressure-bearing performance of polyester fiber reinforced polyethylene composite pipes[J]. Engineering Plastics Application, 2020, 48(8): 118–122.
|
[17] |
YUE Qianjin, LU Qingzhen, YAN Jun, et al. Tension behavior prediction of flexible pipelines in shallow water[J]. Ocean Engineering, 2013, 58: 201–207. doi: 10.1016/j.oceaneng.2012.11.002
|
[18] |
周腾. 纤维增强塑料复合管关键承载性能仿真模拟[D]. 西安:长安大学,2021.
ZHOU Teng. Simulation of the key bearing performance of fiber reinforced plastic composite pipes[D]. Xi’an: Chang’an University, 2021.
|
[19] |
张尹. 纤维缠绕增强复合管在轴对称荷载下的力学行为研究[D]. 杭州:浙江大学,2015.
ZHANG Yin. Investigation on mechanisms of filament-wound fiber-reinforced composite pipe under axisymmetric loads[D]. Hangzhou: Zhejiang University, 2015.
|
[20] |
张杰,梁博丰. 复杂载荷下钢带缠绕增强复合管力学特性[J]. 复合材料学报,2021,38(1):246–254.
ZHANG Jie, LIANG Bofeng. Mechanical properties of reinforced composite pipe wound with steel strip under complex loads[J]. Acta Materiae Compositae Sinica, 2021, 38(1): 246–254.
|
[21] |
邹宵. 海洋复合软管结构设计的关键力学问题研究[D]. 大连:大连理工大学,2021.
ZOU Xiao. Research on key mechanical problems of marine composite hose structure design[D]. Dalian: Dalian University of Technology, 2021.
|
[22] |
黄婷. 钢丝缠绕增强塑料复合管的力学性能分析与研究[D]. 杭州:浙江大学,2014.
HUANG Ting. Study on mechanisms of plastic pipe reinforced by helically cross-winding steel wire[D]. Hangzhou: Zhejiang University, 2014.
|
[23] |
刘婷. 钢带缠绕复合管力学性能及可靠性分析[D]. 杭州:浙江大学,2018.
LIU Ting. Mechanical behaviors and reliability analysis of steel strip reinforced flexible pipes[D]. Hangzhou: Zhejiang University, 2018.
|
[24] |
XU Yuxin, BAI Yong, FANG Pan, et al. Structural analysis of fibreglass reinforced bonded flexible pipe subjected to tension[J]. Ships and Offshore Structures, 2019, 14(7): 777–787. doi: 10.1080/17445302.2018.1564534
|
[25] |
冯德华,綦耀光,余焱群. 海洋纤维增强复合柔性管拉伸性能[J]. 中国石油大学学报(自然科学版),2021,45(4):146–152. doi: 10.3969/j.issn.1673-5005.2021.04.018
FENG Dehua, QI Yaoguang, YU Yanqun. Tensile properties of marine fiber reinforced composite flexible pipes[J]. Journal of China University of Petroleum (Edition of Natural Science), 2021, 45(4): 146–152. doi: 10.3969/j.issn.1673-5005.2021.04.018
|
[26] |
陈伟. 涉内压荷载下钢带缠绕增强复合管力学性能分析[D]. 杭州:浙江大学,2018.
CHEN Wei. Research on the mechanical properties of steel strip reinforced flexible pipe under internal pressure related load[D]. Hangzhou: Zhejiang University, 2018.
|
1. |
龚勇,蒋凯,闫许峰,康建平,林黎明. 钻具旋转时同心窄间隙环空流动的数值模拟. 钻探工程. 2023(06): 56-62 .
![]() | |
2. |
田野,蒋东雷,马传华,徐一龙,于晓东,宋洵成. 钻柱偏心旋转对环空摩阻压降影响的数值模拟研究. 石油钻探技术. 2022(05): 42-49 .
![]() |