-
油气资源勘探开发过程中,井壁坍塌、钻头泥包、落物入井等多种原因均能造成卡钻事故,卡钻后需要迅速解卡,避免出现更复杂的井下事故。处理各种卡钻事故时,相对准确的判断卡点深度的方法是利用测卡仪测量,但因使用仪器存在需专业技术人员配合、费用高、受地磁影响以及到井时间长等限制条件,现场多采用提拉法计算[1-4]。该方法需要运用卡点计算公式(自由管柱在一定拉力下的弹性伸长量)求解卡点深度。目前,该公式在计算卡点时并未考虑摩阻影响,虽有误差但对低摩阻直井基本有效,可对定向井、大斜度井等复杂井身结构井而言,较大的摩阻将导致计算公式误差加大,甚至出现错误结果[5-7]。因此,有必要在考虑摩阻因素条件下对现有方法进行修正。通过对上提和下放过程中管柱受力进行分析,并基于卡点深度计算原理,运用胡克定律计算了管柱整体受力与伸长量的关系,指定了特定的操作方式、实施步骤以及取值,为卡点深度计算提供参考。
-
依据材料力学中胡克定律可知,自由管柱(卡点以上管柱)的伸长量与拉力成正比,与自由管柱的长度成正比,与管柱的横截面积成反比,与钢材的弹性系数成反比,由此推导出目前常用卡点计算公式为
$$ L=\frac{\Delta {L}_{1-2}AE}{\Delta f} $$ (1) 其中
$\Delta {L}_{1-2}={\Delta L}_{1}-{\Delta L}_{2}$ $$ \Delta f={F}_{2}-{F}_{1} $$ 式中,L为卡点深度,m;Δf为两次轴向拉力差值 kN;ΔL1−2为两次上提管柱伸长量差值,m;ΔL1、ΔL2分别为轴向拉力F1、F2时管柱的伸长量m;F1、F2分别为超过自由管柱悬重G且小于管柱抗拉强度的任意拉力,且F1<F2,kN;E为管柱弹性模量,2.1×108 kPa;A为管柱横截面积,m2。
现场进行操作时在井下被卡管柱的抗拉强度范围内(即弹性形变范围内)用一定的上提力上提管柱,测得管柱在该上提力下的伸长量,然后计算被卡管柱的截面积,最后计算出卡点的深度,确定卡点位置。具体实施方式如图1所示,将卡钻后的管柱以大于卡钻前悬重的轴向拉力F1缓慢连续上提至F2,分别记录从F1到F2的管柱伸长量ΔL1−2和轴向拉力差值Δf。如果是单一管柱,直接利用式(1)即可求出卡点深度。
对于复合管柱,利用式(1)分别求出n段不同外径和壁厚的管柱伸长量
$\Delta L_{1 - 2}'$ 、${\Delta L}_{1-2}''$ 、···、$\Delta L_{1 - 2}^n $ ,对比实测ΔL1−2,当ΔL1−2<$\Delta L_{1 - 2}' $ 时,则卡点位于第1段管柱内,卡点深度采用式(1)计算;当$\Delta L_{1 - 2}' $ <ΔL1−2<$\Delta L_{1 - 2}'' $ 时,则卡点位于第2段管柱内,卡点深度采用式(2)计算:当$\Delta L_{1 - 2}' $ +···+$\Delta L_{1 - 2}^{n - 1} $ <ΔL1−2<$\Delta L_{1 - 2}' $ +···+$\Delta L_{1 - 2}^n $ 时,则卡点位于第n段管柱内,卡点深度采用式(3)计算。$$ L={L}'+\frac{\left(\Delta {L}_{1-2}-{\Delta L}_{1-2}'\right) {A}''E}{\Delta f} $$ (2) $$ L={L}'+{L}''+\cdots +{L}^{n-1}+\frac{\left(\Delta {L}_{1-2}-{\Delta L}_{1-2}'-\cdots -{\Delta L}_{1-2}^{n-1}\right){A}^{n}E}{\Delta f} $$ (3) 式中,
$ {\Delta L}_{1-2}^{n} $ 为长度为${L}_{}^{n}$ 的第n段管柱在两次轴向拉力下的伸长量差值,m。 -
目前现场采用的提拉法卡点计算公式均基于胡克定律,但对公式中关键参数拉力F、伸长量ΔL的测量和取值方式不同以及摩擦阻力的影响导致了不同的计算结果,精度有所差别。通过在计算公式中加入修正系数修正伸长量ΔL可有效消除接头及加厚部分的影响。此外,在假设摩擦因数为某经验值的基础上推导的考虑摩阻的卡点深度计算公式,也能适度提高卡点计算精度[8-12]。
由于摩擦因数受井况、地层等诸多因素影响,取值不同导致计算结果差异较大。通过深入研究,并对现有方法归纳总结、优化,利用上提下放过程中,相同钩载条件下管柱所受摩擦阻力大小相等、方向相反的特点,提出了一种可消除摩阻影响的卡点深度计算公式,其具体操作方法和实施步骤为:
(1)缓慢连续上提管柱到钩载等于F1,G<F1,在转盘面以上高度h处管柱部位做一标记a (图2);
图 2 消除摩擦阻力的卡点深度计算方法实施步骤
Figure 2. Implementation step of the free point depth calculation method that can eliminate the influence of friction drag
(2)缓慢连续上提管柱到钩载等于F2,G<F1<F2<Fmax,在转盘面以上高度h处的管柱部位做第2个标记b;
(3)缓慢连续上提管柱到钩载等于F3,F2<F3<Fmax;
(4)缓慢连续下放管柱减小管柱钩载至F2,在转盘面以上高度h处的管柱部位做第3个标记c;
(5)缓慢连续下放管柱减小管柱钩载至F1,同样在转盘面以上高度h处的管柱部位做第4个标记d;
(6)分别测量ab、ac和ad,令F2−F1=Δf,ab=ΔL1−2,ac−ad=ΔL2−1;另外考虑管柱接头及加厚部分的影响以提高计算结果的准确性,引入修正系数S (钻杆S值取值见参考文献[11],常用管柱S值取值见表1),则精确的卡点深度L为
$$ L=\frac{{\Delta L}_{1-2}+{\Delta L}_{2-1}}{2\Delta f}EAS $$ (4) 表 1 常用管柱卡点深度修正系数S
Table 1. Free point depth correction coefficient S of commonly used pipe strings
外径/mm S值 外径/mm S值 73 0.9923 127 1.0405 88.9 0.9920 139.7 1.0404 图3为管柱保持拉伸状态时的受力示意图。当上提管柱时,接触点i的摩擦阻力方向沿切线向下,轴向拉力、浮力与重力、摩擦阻力之差为胡克定律中的拉力F。设采用悬重为Fx轴向拉力提升管柱时,卡点以上管柱与井壁接触点i的摩擦因数为
$ {\mu }_{i} $ ,井斜角为$ {\theta }_{i} $ ,管柱自重为G,管柱所受浮力为Fxf,则摩擦力为$ {f}_{xi}=\left({F}_{x}-G+{F}_{x{\rm{f}}}\right)\sin{\theta }_{i}{\mu }_{i} $ 。上提管柱时考虑摩阻的胡克定律为
$$ \begin{split} &E = \frac{{\left[ {{F_1} - G + {F_{1{\rm{f}}}} - \displaystyle\sum\nolimits_{i = 1}^n {\left( {{f_{1i}} \cdot {\rm{cos}}{\theta _i}} \right)} } \right]\Big/A}}{{\Delta {L_1}/L}} \\ & =\frac{{\left[ {{F_2} - G + {F_{2{\rm{f}}}} - \displaystyle\sum\nolimits_{i = 1}^n {\left( {{f_{2i}} \cdot {\rm{cos}}{\theta _i}} \right)} } \right]\Big/A}}{{\Delta {L_2}/L}} \\ & = \frac{{ {\left[ {{F_2}\!\!+\!\!{F_{2\rm{f}}}\!\!-\!\!\displaystyle\sum\nolimits_{i = 1}^n {\left( {{f_{2i}} \cdot {\rm{cos}}{\theta _i}} \right)} } \right]\!\!-\!\!\left[ {{F_1}\!\!+\!\!{F_{1\rm{f}}}\!\!-\!\!\displaystyle\sum\nolimits_{i = 1}^n {\left( {{f_{1i}} \cdot {\rm{cos}}{\theta _i}} \right)} } \right]} }}{A{\left( {\Delta {L_2}\!\!-\!\!\Delta {L_1}} \right)/L}} \end{split} $$ (5) 图3b为管柱保持拉伸状态时下放过程受力示意图。与图3a不同,当管柱下放时,摩擦阻力的方向向上,因此轴向拉力、浮力、摩擦阻力之和与重力之差即为胡克定律中的拉力F。下放管柱时考虑摩阻的胡克定律为
$$ \begin{split} &E = \frac{{\left[ {{F_2} - G + {F_{2\rm{f}}} + \displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{2i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]\Big/A}}{{\Delta {L_3}/L}} \\ &=\frac{{\left[ {{F_1} - G + {F_{1\rm{f}}} + \displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{1i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]\Big/A}}{{\Delta {L_4}/L}} \\ & = \frac{{{\left[ {{F_2}\!\!+\!\!{F_{2\rm{f}}}\!\!+\!\!\displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{2i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]\!\!-\!\!\left[ {{F_1}\!\!+\!\!{F_{1\rm{f}}}\!\!+\!\!\displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{1i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]} }}{A{\left( {\Delta {L_3}\!\!-\!\!\Delta {L_4}} \right)/L}} \end{split} $$ (6) 将式(5)和式(6)两边分别相加后变换可得
$$ \begin{split} &\left( {\Delta {L_{1-2}}+\Delta {L_{2-1}}} \right)E/L = \\ & \left\{ {\left[ {{F_2}\!\!+\!\!{F_{2{\rm{f}}}}\!\!-\!\!\displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{2i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]\!\!-\!\!\left[ {{F_1}\!\!+\!\!{F_{1{\rm{f}}}}\!\!-\!\!\displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{1i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]} \right\}\Big/A +\!\!\\ & \left\{ {\left[ {{F_2}\!\!+\!\!{F_{2{\rm{f}}}}\!\!+\!\!\displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{2i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]\!\!-\!\!\left[ {{F_1}\!\!+\!\!{F_{1{\rm{f}}}}\!\!+\!\!\displaystyle\sum\nolimits_{i = 1}^n \left( {{f_{1i}} \cdot {\rm{cos}}{\theta _i}} \right)} \right]} \right\}\Big/A \end{split} $$ 即
$$ \left({\Delta L}_{1-2}+{\Delta L}_{2-1}\right)E/L=2\left({F}_{2}-{F}_{1}\right)/A $$ 其中
$ \Delta {L}_{2-1}={\Delta L}_{3}-{\Delta L}_{4}=ac-ad $ 式中,n为接触点数量;ΔL3、ΔL4分别为采用轴向拉力F2、F1下放管柱时管柱的伸长量,m;ΔL2−1为两次下放管柱伸长量差值,m。
因此,考虑摩擦因数的精确卡点深度计算公式即证明为式(4)。
-
(1)精确测量ΔL是准确计算卡点的前提。首先要校正指重表灵敏针,并通过灵敏针读取F值,不要通过悬重针读值。测量ΔL时,先用直角钢板尺放置在转盘面上,在管柱上划出初始起拉点记号(位置与转盘面平齐,其初始提拉力大于原悬重),提拉管柱一定载荷后再采用直角钢板尺划出与转盘面平齐的记号,并做上提拉次数序号,然后仔细测量两记号之间的距离。注意不要将记号画的过粗,以免造成测量误差过大。
(2)应确保初始记号所提拉载荷大于卡点以上悬重,即保证卡点以上管柱处于受拉状态。做好初始记号然后进行提拉卡点工作,就消除了钻具可能在压缩状态下受拉伸长计入提拉长度内,避免影响测量数据,进而影响计算的卡点深度误差过大。
(3)管柱与井壁之间存在一定的摩擦力,导致测量的钻具弹性伸长比实际弹性伸长小,因此在进行提拉法测卡点时,要先稍微上提钻具使其具备一个初始伸长量后再以一定拉力F连续上提,测量拉力F作用下的管柱伸长量ΔL1;再下放管柱直至拉力变为F,测钻具收缩量计入本文卡点计算公式进行计算,即可消除井壁摩擦力造成的影响,特别是对于摩阻大的定向井水平井更为重要。
(4)计算公式中的S值是经过理论计算得出的,在有条件时可以在拉力机上测得实际值对其进行实际校核(将钻具公母接头按规定扭矩上紧后进行拉力机实验测量),这样会使本文推算的卡点公式计算更加准确。
(5)对于复合管柱组合,绝对不能盲目代入本文公式进行计算,一定要按复合管柱测量计算方法计算,根据计算结果分析出卡点在哪一段尺寸类型的钻具上。否则计算出的卡点深度将是错误的,影响后续的故障处理方案和措施,有可能导致井下情况更加复杂。对于复合管柱的卡点计算,若计算出 L < L 1,则卡点在复合管柱中的上面第1组管柱上,通过应用单一管柱卡点计算公式即可计算出卡点深度;若计算出 L = L 1,则卡点在第1和第2种管柱连接处;若计算出 L > L 1,则卡点在第1种管柱之下,依次计算最终确认卡点深度。
(6)为了更加准确计算卡点深度,减少人为测量数据误差,现场实际操作中可以多做几次提拉卡点试验,分别计算出每次的卡点深度,再采用平均值作为最终的卡点深度。
-
该井完钻井深4 248 m,造斜点3 350 m,最大井斜为4 181.79 m处的 92.1°,Ø139.7 mm套管下深4 238.77 m,井下管柱为Ø73 mm J55。作业施工过程中,起原井管柱2.3 m遇卡,经过2次测量拉力与悬重,进行卡点计算,结果见表2。
表 2 Q69-P1井测卡数据
Table 2. Free point measurement data of Well Q69-P1
拉力/kN 伸长量/m 卡点位置/m 上提①:400→600 2.50 3 627 下放①:600→400 2.40 上提②:400→600 2.51 3 649 下放②:600→400 2.42 通过测卡仪测得卡点深度3 645 m。,对比表2中2次操作计算的卡点深度,误差均小于0.5%。
-
该井因采取固结堵漏导致卡钻,应用本文方法计算卡点深度,同时与注磁测卡仪测量结果进行对比验证。该井Ø215.9 mm井眼完钻,井深2 923.00 m,Ø311.15 mm井眼2018.00 m,井下管柱Ø127 mm 钻杆+ Ø139.7 mm钻杆。堵漏作业后,起钻至2 680.00 m静止堵漏,6.5 h后上提管柱遇卡。后经过3次测量拉力与悬重计算卡点深度,结果见表3。
表 3 JY143-5HF井测卡数据
Table 3. Free point measurement data of Well JY143-5HF
拉力/kN 伸长量/m 卡点位置/m 上提①:1100→1200 0.22 1 993 下放①:1200→1100 0.21 上提②:1100→1400 0.64 1 992 下放②:1400→1100 0.62 上提②:1000→1160 0.36 1 995 下放②:1160→1000 0.34 在爆破松扣前使用注磁测卡仪测得卡点深度1 997 m,本文方法计算结果与实测卡点深度对比,误差均小于0.5%。
值得注意的是:(1)由于管柱在提拉过程中伸长量远远小于管柱自身入井长度,因此由于伸长量导致的管柱浮力变化在计算过程中可忽略不计,即F1f=F2f ;(2)实施过程中应尽量保持管柱缓慢连续作业,有利于减小在同一拉力条件下提拉和下放管柱时接触点i的摩擦阻力波动;(3) 3次提拉力的差值应尽可能大,且F2−F1>F3−F2,同时Fx必须小于管柱抗拉强度,如此可减小操作误差,提高计算精度;(4)该方法虽然与文献[7] (准确计算管柱卡点深度的方法)实施方式相似,但细化的步骤3和步骤4考虑了上提和下放过程中,轴向拉力为F2时管柱摩擦阻力方向相反的客观事实,
$ {\Delta L}_{2-1}\ne \Delta {L}_{2} $ ,因此卡点深度计算结果更为准确。 -
(1)利用传统提拉法计算卡点深度时,提拉方法和实施步骤的不同及管柱伸长量取值不同,导致利用胡克定律进行卡点深度计算的结果不准确,而摩阻的存在加大了计算结果的误差,并且随着摩阻的增加误差加大。
(2)忽略上提下放过程中管柱浮力变化影响,通过特定的提拉方法、实施步骤和精确测量数据,可基本消除管柱摩擦阻力对卡点计算深度的影响。
Accurate free point calculation method and its application
-
摘要: 钻井过程中由于地层或操作不当等因素导致偶发的卡钻,不仅延长了钻井周期,而且可能导致次生故障甚至埋钻等事故,造成严重经济损失。目前,现场采用的提拉法计算卡点方法没有考虑井内管柱所受摩擦阻力影响,导致计算出的卡点深度误差随着摩阻增加而增加,难以适用于高摩阻直井及大斜度井。基于卡点深度计算原理,运用胡克定律,利用上提、下放过程中相同钩载条件下管柱在井眼内所受摩擦阻力大小相等、方向相反的特点,推导了可消除套管内管柱摩擦阻力影响的卡点深度计算公式,同时提出了该计算公式配套的操作方法和实施步骤。实践证明,该方法计算卡点深度结果准确,具有较强的实用价值。Abstract: In the process of drilling, formation factors or misoperation can lead to incidental pipe sticking, which not only extends the drilling cycle, but also may cause secondary failure and even pipe burying and other accidents, resulting in great economic loss. At present, the pulling method adopted on site to calculate free points doesn’t consider the influence of friction drag applied on downhole pipe strings, which increases the calculation error of free point depth with the increase of the friction drag, so it is hardly applicable to high-friction vertical wells and high-angle deviated wells. In this paper, the free point depth calculation formula which can eliminate the influence of the friction drag applied on the pipe string inside the casing was derived by making use of the performance that the friction drag applied on the downhole pipe string under the same hook load in the process of pulling and lowering is equal in value and opposite in direction, according to the free point depth calculation principle and the Hooke’s Law. What’s more, the operation method and implementation step corresponding to this calculation formula were put forward. The practice proves that this method can calculate the free point depth accurately and has a great practical value.
-
Key words:
- drilling /
- pipe sticking /
- free point depth /
- Hooke’s Law /
- pulling method /
- friction drag
-
表 1 常用管柱卡点深度修正系数S
Table 1. Free point depth correction coefficient S of commonly used pipe strings
外径/mm S值 外径/mm S值 73 0.9923 127 1.0405 88.9 0.9920 139.7 1.0404 表 2 Q69-P1井测卡数据
Table 2. Free point measurement data of Well Q69-P1
拉力/kN 伸长量/m 卡点位置/m 上提①:400→600 2.50 3 627 下放①:600→400 2.40 上提②:400→600 2.51 3 649 下放②:600→400 2.42 表 3 JY143-5HF井测卡数据
Table 3. Free point measurement data of Well JY143-5HF
拉力/kN 伸长量/m 卡点位置/m 上提①:1100→1200 0.22 1 993 下放①:1200→1100 0.21 上提②:1100→1400 0.64 1 992 下放②:1400→1100 0.62 上提②:1000→1160 0.36 1 995 下放②:1160→1000 0.34 -
[1] 刘琮洁, 周祥易, 张燕, 等. 管柱卡点计算公式及应用[J]. 江汉石油学院学报, 2001, 23(4):74-76. doi: 10.3969/j.issn.1000-9752.2001.04.027 LIU Congjie, ZHOU Xiangyi, ZHANG Yan, et al. Pipe string stuck-point calculating equations and their application[J]. Journal of Jianghan Petroleum Institute, 2001, 23(4): 74-76. doi: 10.3969/j.issn.1000-9752.2001.04.027 [2] 周劲辉, 高德利, 王宇新. 钻柱卡点预测实验[J]. 实验力学, 2010, 25(5):575-580. ZHOU Jinhui, GAO Deli, WANG Yuxin. Experiment on the stuck-point prediction of drill string[J]. Journal of Experimental Mechanics, 2010, 25(5): 575-580. [3] ADNOY B S, LARSEN K, BERG P C. Analysis of stuck pipe in deviated boreholes[J]. Journal of Petroleum Science and Engineering, 2003, 37(3/4): 195-212. doi: 10.1016/S0920-4105(02)00353-4 [4] 刘同茹, 魏群涛, 吴思琼, 等. 测卡松扣仪及其在卡钻事故检测中的应用[J]. 石油钻采工艺, 2006, 28(6):18-20. doi: 10.3969/j.issn.1000-7393.2006.06.006 LIU Tongru, WEI Quntao, WU Saiqiong, et al. Free point back-off instrument and its application for releasing downhole stuck pipe[J]. Oil Drilling & Production Technology, 2006, 28(6): 18-20. doi: 10.3969/j.issn.1000-7393.2006.06.006 [5] 郭凤, 刘合, 张劲. 卡点计算的有限元分析[J]. 机械工程学报, 2007, 43(1):158-163. doi: 10.3321/j.issn:0577-6686.2007.01.026 GUO Feng, LIU Ge, ZHANG Jin. FEM analysis of the free point of drill string[J]. Chinese Journal of Mechanical Engineering, 2007, 43(1): 158-163. doi: 10.3321/j.issn:0577-6686.2007.01.026 [6] 周劲辉, 张勇, 高德利. 大位移井卡点预测实验[J]. 实验室研究与探索, 2015, 34(5):11-15. doi: 10.3969/j.issn.1006-7167.2015.05.003 ZHOU Jinhui, ZHANG Yong, GAO Deli. Experimental research on prediction of stuck-point of downhole tubular string in extended reach well[J]. Research and Exploration in Laboratory, 2015, 34(5): 11-15. doi: 10.3969/j.issn.1006-7167.2015.05.003 [7] 刘东明. 井下被卡管柱的卡点公式计算法优化探讨[J]. 中国海上油气, 2004, 16(1):47-50. doi: 10.3969/j.issn.1673-1506.2004.01.009 LIU Dongming. The discussion of the formula calculating method for the free point of downhole string stuck[J]. China Offshore Oil and Gas, 2004, 16(1): 47-50. doi: 10.3969/j.issn.1673-1506.2004.01.009 [8] 许清海, 胡晋阳, 熊燃, 等. 考虑摩阻的圆弧井卡点深度计算模型研究[J]. 石油机械, 2018, 46(3):84-88, 94. doi: 10.16082/j.cnki.issn.1001-4578.2018.03.015 XU Qinghai, HU Jinyang, XIONG Ran, et al. Study on the calculation model of free point depth of circular arc well[J]. China Petroleum Machinery, 2018, 46(3): 84-88, 94. doi: 10.16082/j.cnki.issn.1001-4578.2018.03.015 [9] 况雨春, 熊威, 张立民, 等. 大斜度井修井作业卡点预测数值模拟方法[J]. 石油机械, 2014, 42(4):57-61. doi: 10.3969/j.issn.1001-4578.2014.04.014 KUANG Yuchun, XIONG Wei, ZHANG Limin, et al. Numerical simulation method for predicting the sticking point in highly deviated well workover[J]. China Petroleum Machinery, 2014, 42(4): 57-61. doi: 10.3969/j.issn.1001-4578.2014.04.014 [10] 韩志勇. 卡点深度计算中消除摩阻力影响的方法[J]. 石油钻探技术, 2010, 38(1):1-3. doi: 10.3969/j.issn.1001-0890.2010.01.001 HAN Zhiyong. Eliminating disturbance of drag to sticking depth calculations[J]. Petroleum Drilling Techniques, 2010, 38(1): 1-3. doi: 10.3969/j.issn.1001-0890.2010.01.001 [11] 蒋希文. 钻井事故与复杂问题[M]. 北京: 石油工业出版社, 2006. JIANG Xiwen. Drilling accidents and complex problems[M]. Beijing: Petroleum Industry Press, 2006. [12] 冉津津, 李长健, 唐义祥, 等. 油水井测卡点技术研究与应用[J]. 石油天然气学报, 2013, 35(5):134-137. doi: 10.3969/j.issn.1000-9752.2013.05.028 RAN Jinjin, LI Changjian, TANG Yixiang, et al. Research and application on prediction of stuck-point[J]. Journal of Oil and Gas Technology, 2013, 35(5): 134-137. doi: 10.3969/j.issn.1000-9752.2013.05.028 -