水平井钻井过程中井底钻压预测及应用

吴泽兵 郭龙龙 潘玉杰

西安石油大学机械工程学院

水平井摩阻大，不容易施加井底钻压，使得地面钻压与井底钻压差别较大。然而，井底钻压是影响钻头性能的主要因素，因此研究水平井钻井的井底钻压非常有必要。

国内已有学者对井底钻压进行了相关研究。苏义脑［1］导出了造斜、稳斜等不同井段井底钻压的计算公式。杜仕强在文献［2］中提及钻压，但没有给出具体的计算。岳钦鹏［3］在“关于对定向井井底钻压分析的研究”一文中也只是给出了钻压的计算公式。邵冬冬等［4］利用实验装置进行了水平井井底钻压波动实验，认为水平井井底钻压呈正弦波动，井底的实际钻压与名义钻压近似。郑双进等［5］考虑钻井液对钻柱轴向力的影响和定向井钻柱结构，提出钻柱和钻具组合轴向力的计算方法。孙召红等［6］在分析靠近钻头处测量短节受力的基础上，应用应变原理建立了钻压、扭矩测量模型。肖建波等［7］采用弹性梁的变形微分方程，建立了钻进时轴向力的二维计算模式，得出DTU钻具组合的井底钻压计算模式，并给出钻进时的钻重和视压力。陈永明等［8］根据力的叠加原理，对井底钻压在各个井段产生的钻柱内力进行了分析，导出井口钻压与井底钻压之间的关系，认为井口钻压与井底钻压的比值与摩擦因数、井底井斜角有关；在水平井段，井口钻压与井底钻压其比值达高1.60。以上学者尽管针对井底钻压做了相关研究，但是没有给出摩擦因数的计算过程，也没有对模型进行验证。

国外，G.Hareland等［9-16］利用软杆模型计算了大钩载荷、摩擦因数及井底钻压，并利用有限元法进行了类似的探讨和计算。Johancsik等［17］首次开发了用于预测钻柱摩阻扭矩的计算机模型。Aadnoy等［18］建立了适用于直井段、增斜段、降斜段和侧向弯曲段的摩阻扭矩模型。Johancsik模型和Aadnoy模型都基于软杆假设，忽略了钻柱的刚度。随后，部分研究者在解析模型中考虑了刚度，但假设钻柱与井壁是连续接触。此外，有学者利用数值求解的方法对钻柱的摩阻扭矩进行分析。Yang D等［19］提出了一种三维有限差分方法用于静态载荷下底部钻具组合的分析，然而该模型仅仅分析了底部钻具组合部分，没有涉及整个钻柱、摩阻和扭矩的计算。

笔者着重介绍了钻柱与井壁间摩擦因数、井底钻压的计算模型和方法。同时，使用开发的程序计算摩擦因数和井底钻压，并将计算值与现场测量值进行对比，以验证模型和程序的准确性。本文开发的模型和程序可用于井后分析，未来也可将其与常规自动送钻系统集成，从而提高钻头的性能和整体钻井效率、降低钻井成本。

1 摩阻扭矩建模

摩阻和扭矩建模对井身设计相当有用，其有助于预测和防止出现钻井事故。通常，井的垂直段几乎没有摩阻和扭矩，扭矩和摩阻出现在曲线段、斜直井段，也必然出现在水平段。常用的Johancsik模型和Aadnoy 3D模型都忽略了钻柱的刚度，认为钻柱是有重量的软绳。以下建立的摩阻扭矩模型是基于前面2个模型的综合模型，表达式见式（1）～（4）。

斜直井段

弯曲井段钻柱受拉时

弯曲井段钻柱受压时

其中

式中，Ft为钻柱上端受力，N；β为浮力系数；w为单位长度钻柱的重量，N/m；ΔL为钻柱长度，m；α为井斜角，rad；μ为摩擦因数；Fb为钻柱下端受力，N；θ为井偏角，rad；θ为方位角，rad；Fn为钻柱法向力，N。

根据上述模型，从钻头开始计算轴向力一直到大钩，最后可以得到大钩载荷。井底钻压的计算分2步，先用钻头空钻数据计算钻柱与井壁间的摩擦因数，然后用所得摩擦因数预测钻井时的井底钻压。

2 摩擦因数计算

根据上述摩阻模型，在Visual C++ 2013的集成开发环境下，利用C#开发计算井底钻压的程序。钻柱与井壁间摩擦因数的计算流程如图1所示。

图1 钻柱与井壁间摩擦因数的计算流程Fig.1 Calculation flow chart of friction coefficient between drilling string and borehole wall

上述计算大钩载荷时假定摩擦因数已知，但真实的摩擦因数是未知的。因此，根据已知的大钩载荷反算摩擦因数，具体步骤是：假定一个初始摩擦因数μ，利用摩阻扭矩模型计算大钩载荷HLc；比较计算值HLc与已知大钩载荷HLa（已知的大钩载荷是对直接测出的大钩载荷去除钻井绳与滑轮的摩擦修正得到的），若计算值与修正的测量值很接近，即小于给定的容差ε，则认为假定的摩擦因数为所求；否则重新假设摩擦因数，再次计算大钩载荷、直到满足要求。计算所得加拿大西部某井的摩擦因数如图2所示，可以看出摩擦因数的值不是常数。

图2 摩擦因数随井深变化曲线（井A）Fig.2 Relationship of friction coefficient vs.well depth (Well A)

3 井底压力计算

求出摩擦因数后，便可利用其值计算井底钻压。步骤与摩擦因数的计算相似：假定一个初始钻压，利用所建立的摩阻扭矩模型计算出大钩载荷；比较计算的大钩载荷与已知的修正过的大钩载荷，如果两者很接近，即小于给定的容差，则认为假定的钻压为所求值；否则重新假设钻压，重复前面的步骤计算大钩载荷，直到满足要求为止。

图3为井底钻压预测值与测量值的比较，钻井数据来自加拿大西部。可以看出，无论是趋势还是数值，计算值和测量值匹配得都比较好。然而，有些地方也存在差异。这些差异主要与测量工具、模型、钻柱振动、工具的磨损等因素有关。

图3 井底钻压预测值与测量值对比（井A）Fig.3 Comparison between predicted bottom hole WOB and the measured value (Well A)

正常情况下，井底钻压比地面钻压小，但具体小多少不可而知。利用上述模型和程序可以得到井底钻压与地面钻压的准确数值差异，如图4所示。

图4 预测的井底钻压与地面钻压之间的对比（井B）Fig.4 Comparison between predicted bottom hole WOB and ground WOB (Well B)

4 应用

4.1 钻井事后分析

对于有现存钻井数据的情况，开发的模型和程序可用于钻井事后分析，分析结果可以用于指导钻相似条件的井，例如具有相似钻机、相似地层、相似井眼几何参数等。显著的例子就是指导钻偏移井（offset well）。

4.2 实时井底钻压监控

开发的模型和程序可用于实时钻井操作，实时预测和显示井底钻压、钻头扭矩，也可以显示其他重要钻井参数，如：钻柱受的最大压力，从而监控钻柱是否接近失稳状态。

4.3 增强自动送钻功能

目前部分钻机已安装了自动钻井设备，如自动钻井仪或自动送钻系统，但这种设备有不足之处。一方面，其难以准确地控制井底钻压，也无法获得井底钻压的精确值，除非实际井下测量，因此使用了过大的余量；另一方面，自动送钻系统主要控制表面钻压，而井底钻压是小于表面钻压的，具体小多少无法定量。若将本研究的井底钻压计算模型和程序与基于地面钻压的自动送钻系统集成，形成一种新型自动送钻系统，称其为“基于井底钻压的自动送钻系统”。此新系统不仅具有普通自动系统的优点，而且能够减小甚至消除这个余量，从而施加准确的井底钻压，使钻头性能达到最优。图5为井底钻压为170 kN时，另外3口井地面钻压的变化情况。结果表明：若给自动送钻系统输入图示地面钻压值，就可以在井底获得170 kN的常钻压值。

图5 给定常井底钻压时的地面钻压Fig.5 Ground WOB at the given constant bottom hole WOB

5 结论

（1） 摩擦阻力和井底钻压预测程序能依据地面钻井数据比较准确地预测摩擦因数、井底钻压，井底钻压预测值与现场测量值在变化趋势和数值方面都比较吻合。

（2） 如果有现存的钻井数据，可利用程序进行钻井事后分析，以便指导同类情况下的钻井作业。

（3） 该程序也可以实时预测和显示井底钻压以及井底钻头扭矩，确保安全钻进。

（4） 将该程序与常规自动送钻系统集成，可以实现准确控制井底钻压的自动钻井操作，从而提高钻头性能和整体钻井效率、降低钻井成本。

参考文献：

［1］苏义脑.求定向井井底真实钻压值的理论分析和初步计算［J］.石油钻采工艺，1991，13（2）：31-42.SU Yinao.Theoretical analysis and initial calculation of downhole weight on bit in directional drilling［J］.Oil Drilling & Production Technology, 1991, 13（2）: 31-42.

［2］杜仕强.定向井、水平井摩阻和钻压的理论分析与研究［J］.中国石油和化工标准与质量，2012（3）：141，159.DU Shiqiang.Theoretical analysis and study of drag and weight on bit in directional and horizontal drilling［J］.China Petroleum and Chemical Standard and Quality,2012（3）: 141, 159.

［3］岳钦鹏.关于对定向井井底钻压分析的研究［J］.中国石油和化工标准与质量，2012（13）：91-92.YUE Qinpeng.Study of downhole weight on bit in directional drilling［J］.China Petroleum and Chemical Standard and Quality, 2012（13）: 91-92.

［4］邵冬冬，管志川，张文斌，温欣，史玉才.水平井井底钻压波动规律的实验研究［J］.科学技术与工程，2012，13（15）：4347-4351.SHAO Dongdong, GUAN Zhichuan, ZHANG Wenbin,WEN Xin, SHI Yucai.Experimental research on weight on bit fluctuation in horizontal well［J］.Science Technology and Engineering, 2012, 13（15）: 4347-4351.

［5］郑双进，黄志强，陈彬，周吉祥.定向井钻压传导计算方法［J］.断块油气田，2011，18（3）：403-405.ZHENG Shuangjin, HUANG Zhiqiang, CHEN Bin,ZHOU Xiangji.Calculating method for WOB conducting of directional well［J］.Fault-Block Oil & Gas Field,2011, 18（3）: 403-405.

［6］孙召红，房军，盛利民，窦秀荣.井底钻压扭矩传感器设计与分析［J］.石油矿场机械，2010，39（4）：65-68.SUN Zhaohong, FANG Jun, SHENG Limin, DOU Xiurong.Deisgn and analysis of torque and bit weight sensor［J］.Oil Field Equipment, 2010, 39（4）: 65-68.

［7］肖建波，雷宗明.水平井井底钻压的研究［J］.天然气工业，2001，21（3）：41-44.XIAO Jianbo, LEI Zongming.A research on bottom hole bit weight in a horizontal well［J］.Natural Gas Industry,2001, 21（3）: 41-44.

［8］陈永明，闫振来.定向井井底真实钻压的计算方法［J］.钻采工艺，1996，19（4）：1-5.CHEN Yongmin, YAN Zhenlai.Calculation method of real weight on bit in directional drilling［J］.Drilling &Production Technology, 1996, 19（4）: 1-5.

［9］WU A, HARELAND G, LEI L.Real-time downhole weight on bit （DWOB） automation in directional drilling［C］.Proceedings of the ASME 2015 34th International Conference on Ocean, Offshore and Arctic Engineering,Newfoundland, Canada, May 31-June 5, 2015.

［10］HARELAND G, WU A, LEI L.The field tests for measurement of downhole weight on bit （DWOB） and the calibration of a real-time dwob model［C］.7th International Petroleum Technology Conference, Doha,Qatar, January 19-22, 2014.

［11］HARELAND G, WU A, FAZAELIZADEH M.Finite element analysis of drillstring and its application on torque and drag calculation［J］.The International Journal of Engineering and Science, 2013, 2（1）: 9-16.

［12］HARELAND G, WU A, LEI L, HUSEIN M M,ZAKARIA M F.Innovative nanoparticle drilling fluid and its benefits to horizontal or extended reach drilling［C］.SPE Canadian Unconventional Resources Conference,Alberta, Canada, October 30-November 1, 2012.

［13］WU A, HARELAND G.Calculation of friction coefficient and downhole weight on bit with finite element analysis of drillstring［C］.The 46th US Rock Mechanics/Geomechanics Symposium, Chicago, IL,June 24–27, 2012.

［14］WU A, HARELAND G, FAZAELIZADEH M.Comparison of analytical and numerical models on torque and hookload calculation［J］.IOSR Journal of Engineering, 2012, 2（1）: 77-83.

［15］WU A, HARELAND G, FAZAELIZADEH M.Torque &drag analysis using finite element method［J］.Modern Applied Science, 2011, 5（6）: 13-27.

［16］FAZAELIZADEH M, HARELAND G, WU A.Real time wellbore friction analysis to detect onset of pipe sticking during extended reach drilling: case study［C］.Brazil Offshore Conference, Macae, Brazil, June 14-17, 2011.

［17］JOHANCSIK C A, FRIESEN D B, DAWSON R.Torque and drag in directional wells prediction and measurement［J］.Journal of Petroleum Technology, 1984, 36（6）:987-992.

［18］AADNOY B S, FAZAELIZADEH M, HARELAND G.A 3-dimentional analytical model for wellbore friction［J］.Journal of Canadian Petroleum Technology, 2010,49（10）: 25-36.

［19］YANG D, RAHMAN M K, CHEN Y.Bottom hole assembly analysis by finite difference differential method［J］.International Journal for Numerical Methods in Engineering, 2008, 74（9）: 1495-1517.