使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

定向井抽油杆柱纵横扭耦合振动仿真研究  PDF

  • 王宏博 1
  • 孙秀荣 2
  • 李凤鑫 2
  • 李伟成 3
  • 辛顺 1
  • 王长亮 1
  • 李东芳 1
1. 海洋石油工程股份有限公司,天津 300451; 2. 河北环境工程学院,河北 秦皇岛 066102; 3. 燕山大学机械工程学院,河北 秦皇岛 066004

中图分类号: O327TE933+.2

最近更新:2024-12-03

DOI:10.16385/j.cnki.issn.1004-4523.2024.11.014

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

抽油杆柱是由接箍与抽油杆通过螺纹连接组成的细长杆柱。对抽油杆柱动力学行为的研究是预防抽油杆柱失效破坏,降低油井运行成本的基础。针对抽油杆柱的纵横扭耦合振动,提出杆管摩擦力产生的绕轴扭矩是扭转振动的激励这一观点。考虑摩擦力对杆柱纵横扭耦合振动的影响,建立了具有初弯曲的定向井抽油杆柱在顶端位移激励、底端载荷激励、弯曲井眼轨道激励与油管约束条件下的纵横扭耦合振动仿真模型。采用四阶龙格‑库塔法实现了对仿真模型的求解和对杆柱动力学行为的仿真模拟。仿真结果表明:杆柱的纵横扭耦合振动对杆柱纵向、横向振动仿真结果均有影响,相对来说对纵向振动的影响并不明显;油井造斜段的杆管接触力大,杆柱受压段的杆管碰撞严重,因此油井的造斜段与杆柱的受压段为杆管偏磨的危险点;正常预紧情况下,杆柱的扭转振动不足以导致抽油杆柱脱扣,但当接箍螺纹预紧力降至杆柱的最大扭矩以下后,扭转振动即可能导致抽油杆柱脱扣。

引 言

抽油杆柱作为将地面抽油机能量传递至井底抽油泵的主要载体,是有杆泵抽油系统的重要组成部

1。在油井生产过程中,抽油杆柱的失效破坏是油井停井作业的主要原因之一。油井停井作业不仅增加了油井的生产成本,还降低了油井的产量,进而导致巨大的经济损失。油井运行过程中杆柱的振动是杆柱失效破坏的主要原因,因此对抽油杆柱振动的研究是预防杆柱失效破坏的基础。

目前大量学者对抽油杆柱的纵向振动与横向振动进行了研究,实现了对油井悬点示功图预测仿真与杆管偏磨现象的分析。1963年,GIBBS

2采用一维有阻尼波动方程建立了描述抽油杆柱纵向振动的仿真模型,并采用有限差分法实现了仿真求解,为抽油杆柱动力学研究提供了思路。后来有研究人员分别考虑柱塞磨损、动静摩擦、温度变化等不同的影响因素,改进了抽油杆柱纵向振动的仿真模3‑7。孙秀荣8考虑屈曲变形激励建立了直井抽油杆柱横向振动仿真模型。

LUKASIEWICZ

9基于纵横弯曲梁理论提出了定向井抽油杆柱在井斜平面的纵横耦合振动仿真模型。XU10‑11采用虚位移原理建立了抽油杆柱在弯曲井眼内的纵横耦合振动仿真模型。WANG12‑13考虑弯曲井眼轨道对轴向运动抽油杆柱横向振动的激励,建立了考虑弯曲井眼轨道激励的抽油杆柱纵横耦合振动仿真模型。在以上杆柱纵向与横向仿真模型中均仅考虑了杆管库仑摩擦力对杆柱纵向振动的影响,而忽略了其对横向振动的影响。

抽油杆柱失效破坏的形式除了杆管偏磨外还包括抽油杆柱的脱扣现象。抽油杆柱的脱扣现象表明:除了纵向和横向振动以外,抽油杆柱还存在扭转振动,即抽油杆柱的动力学行为是纵横扭的耦合振动。与钻柱、螺杆泵不同,抽油机井在运行过程中,抽油杆柱不受系统外部扭矩,因此对抽油杆柱扭转激励的研究也就成为了抽油杆柱扭转振动研究的难点。王文昌

14曾采用有限元法建立了定向井抽油杆柱的纵横扭耦合振动仿真模型,但其模型中未指出抽油杆柱扭转振动的激励。针对这一问题,本文作者结合抽油杆柱纵横耦合振动现象认识到:抽油机井运行过程中弯曲变形的抽油杆柱与油管之间存在接触力与相对运动,相对运动又将导致杆管接触点出现摩擦力,杆管之间的摩擦力作用点分别为油管的内壁以及抽油杆柱表面,作用于抽油杆柱表面上的摩擦力在与抽油杆柱轴线不同向时将产生绕其轴线的扭矩,进而激励抽油杆柱产生扭转振动。因此,杆管摩擦力产生的绕轴扭矩即为抽油杆柱扭转振动的真实激励。

抽油杆柱的摩擦扭矩是导致杆柱扭转振动的主要原因,因此对摩擦扭矩的仿真建模是抽油杆柱扭转振动建模的基础。目前大量学者对钻柱在井眼内的摩擦扭矩进行了分析。JOHANCSIK

15将钻柱简化为没有刚度的绳索,假设钻柱轴线与井眼重合,考虑轴向受力与钻柱自重推导出了钻柱与井壁的接触力模型;在不考虑摩擦力方向的前提下,基于钻柱与井壁接触力以及摩擦系数建立了钻柱摩擦扭矩的仿真模型。HO16考虑弯曲井眼轨道对钻柱弯曲变形的影响,以及弯曲钻柱刚度对钻柱与井壁接触力的影响,建立了钻柱摩擦扭矩的仿真模型。李子丰17考虑钻井液等对钻柱阻力的影响,考虑轴向运动与钻柱自身旋转运动对阻力方向的影响、阻力方向对摩擦扭矩的影响提出了等效摩擦系数的概念,建立了钻柱稳态拉力与扭矩仿真模型。此外,有的学者考虑不同的影响因素建立了描述钻柱动力学行为的有限元模型,但其在计算摩擦扭矩时所采用的模型与上述两类模型并无本质差18‑21。祝效华等22‑24基于Hamilton原理与有限元原理建立了描述钻柱在弯曲井眼轨道内纵横耦合振动的动力学仿真模型,并据此得到了钻柱与井壁接触力的仿真模型,然后以钻柱纵横扭耦合振动的仿真模型为基础,考虑纵向振动与扭转振动对杆柱摩擦力方向的影响建立摩擦扭矩的仿真模型。此外在部分钻柱纵横扭耦合仿真模型中,部分学者还考虑了杆柱横向振动导致的相对速度改变对摩擦扭矩的影25‑26

以上钻柱摩擦扭矩的研究为抽油杆柱摩擦扭矩的研究提供了思路。但针对抽油杆柱这一研究对象,易忽略如下问题:(1)由于钻柱工作特性的影响,自身的旋转运动是影响相对运动速度横向分量的主要原因,因此忽略杆柱横向振动对杆管相对运动速度的影响不大。但抽油杆柱工作过程中不存在主动的旋转运动,即横向振动是影响杆管之间的相对运动的主要因素,在仿真建模时应特别注意。(2)除杆管相对运动速度外,杆柱轴线方向也是影响摩擦扭矩的原因。对于钻柱来说,其轴向速度比较小、横向速度比较大,因此杆管相对运动速度对摩擦扭矩的影响要远大于杆柱轴线方向的影响,以上模型中均忽略了杆柱弯曲形态的影响。但对于抽油杆柱,其工作过程中存在轴向往复运动,轴向速度比较大、横向速度比较小,杆柱弯曲形态对摩擦扭矩的影响明显,因此在建立仿真模型时需要考虑杆柱弯曲变形的影响。

通过对以上研究内容的分析可知,由于纵向振动、横向振动与弯曲井眼轨道的共同作用,抽油杆柱受油管库仑摩擦力的方向不再与杆柱轴线方向一致,进而产生绕杆柱轴线的摩擦扭矩,导致抽油杆柱产生扭转振动。同时杆柱的扭转振动将影响杆柱弯曲剪切力与杆管相对运动,进而影响杆柱振动。

目前已有部分学者对抽油杆柱的纵横扭耦合振动进行了研究,总结其研究,存在如下不足:(1)没有研究弯曲井眼轨道对轴向往复运动抽油杆柱纵横扭耦合振动的影响;(2)未明确指出抽油杆柱扭转振动的激励。

本文针对上述不足,将定向井抽油杆柱的动力学问题概括为:具有初弯曲的抽油杆柱在顶端周期性位移激励、底端载荷激励与弯曲井眼轨道激励下,考虑油管壁与抽油杆柱的弹性碰撞、杆管之间的库仑摩擦力与液体黏滞阻力的纵横扭耦合振动问题,进而建立了定向井抽油杆柱纵横扭耦合振动的仿真模型。采用有限差分法对杆柱长度进行离散,得到常微分方程组,并进一步采用四阶龙格‑库塔法进行求解计算,实现了对定向井抽油杆柱动力行为的仿真模拟。

1 井眼轨道及杆柱位置矢量描述模型

图1为定向井井眼轨道及抽油杆柱空间形态示意图。图中抽油杆柱由一定数量的抽油杆、接箍、扶正器组成。为描述井眼轨道与抽油杆柱的轴向位置,分别建立了固定坐标系和运动坐标系。固定坐标系原点为井口;运动坐标系原点在悬点处,且随悬点运动。以固定坐标系下sw描述井眼轨道上任意点相对于井口的曲线长度,以运动坐标系下s描述抽油杆柱任意点相对于悬点的曲线长度。

图1  井眼轨道及抽油杆柱空间形态示意图

Fig.1  Schematic diagram of the borehole trajectory and space configuration of the rod string

1.1 井眼轨道位置矢量模型

定向井的井眼轨道为一条空间曲线,现以直角坐标系进行描述。以井口为坐标原点,建立坐标系o‑xyz,坐标系中xyz的方向分别为单位向量ijk的方向。结合静位置曲线坐标sw即可得到井眼轨道上任意一点A的位置矢量表示方法:

rAsw=xswi+yswj+zswk (1)

根据井眼轨道测量数据可以得到油井的井斜角α与方位角ϕ。据此可以计算井眼轴线上任意点A处的切向量τA、主法线向量nA、副法线向量bA、井眼轴线曲率kA与井眼轴线挠率TA,如下式表示:

τA=sinαcosϕi+sinαsinϕj+cosαknA=1kAτAswbA=τA×nAkA=α˙2+ϕ˙2sin2αTA=kA-2      sinαα˙ϕ¨-α¨ϕ˙+cosαϕ˙3sin2α+2α˙2ϕ˙ (2)

式中 α˙为井斜角对斜深sw的一阶求导;α¨为井斜角对斜深sw的二阶求导;ϕ˙为方位角对斜深sw的一阶求导;ϕ¨为方位角对斜深sw的二阶求导。

1.2 抽油杆位置矢量模型

以井眼轨道位置矢量模型为基础,通过动位置曲线坐标s描述抽油杆柱上任意点C的轴向位置;以与C点对应的井眼轨道轴线上A点的切线向量τA描述C点的轴向位移;以A点的主法线向量nA与副法线向量bA来描述C点的横向位移。抽油杆柱上任意点C的位置矢量即可表示为:

rCs=rAsw+es,tτA+us,tnA+ws,tbA (3)

式中 est)为t时刻抽油杆柱上任意截面s在井眼轴线τA方向上的位移;ust)为t时刻抽油杆柱上任意截面s在井眼轴线nA方向上的位移;wst)为t时刻抽油杆柱上任意截面s在井眼轴线bA方向上的位移。

根据Frenet公式可求得在抽油杆柱轴线C点处的切向量τCs)、主法线向量nCs)、副法线向量bCs)、杆柱轴线曲率kCs)与杆柱轴线挠率TCs),如下式所示:

τC=XτA+YnA+ZbAnC=kC-1Xs-kAYτA+kAX+Ys-TAZnA+          TAY+ZsbAsscbC=τC×nCkC=sscXs-kAY2+kAX+Ys-TAZ2+          TAY+Zs20.5TC=kC-2τC,τCsc,2τCsc2 (4)

其中:

X=1+es-kAusscY=kAe+us-TAwsscZ=TAu+wssscssc=1+es-kAu2+kAe+us-TAw2+          TAu+ws2-0.5 (5)

式中 sc为变形后杆柱的长度;eses求一阶偏导;usus求一阶偏导;wsws求一阶偏导;XsXs求一阶偏导;YsYs求一阶偏导;ZsZs求一阶偏导。

2 杆柱纵横扭耦合振动仿真模型

为便于研究,同时也为突出本文的研究重点,做如下简化和假设:(1)不考虑抽油机负载对电机转速的影响,假设电动机匀速转动;(2)油管柱的轴线与井眼轴线重合,且不考虑油管柱的振动;(3)忽略油管内液柱的振动;(4)忽略顶端钢丝绳的弹性;(5)抽油杆柱为均质弹性体;(6)杆柱各横截面轴力与杆柱轴线一致,杆柱各截面剪切力与杆柱轴线垂直。

2.1 抽油杆柱纵横扭耦合振动力学模型

根据油田实际生产情况可知,抽油杆柱在弯曲油管内的运动包括在油管约束下的空间位移与沿自身轴线的扭转。

对于抽油杆柱的纵向振动,杆柱顶端可简化为位移的边界条

27‑28;对于横向振动,抽油机悬点至井口盘根处抽油杆柱可以简化为可滑动的固支端;对于扭转振动,杆柱顶端扭转受钢丝绳的约束,因此可以简化为扭矩的边界条件。对于抽油杆柱的纵向振动,杆柱底端受柱塞载荷的影响为力的边界条件;对于横向振动,杆柱的底端可以简化为可滑动的固支端;对于扭转振动,杆柱底端的柱塞可以在油管内自由旋转,因此可以简化为自由端。图2为抽油杆柱纵横扭耦合振动的力学模型图与微元受力图。

图2  抽油杆柱纵横扭耦合振动力学模型

Fig.2  Mechanical model of the longitudinal transverse and torsional coupling vibration of the rod string

图2中:S为抽油机冲程;Pt)为t时刻作用于抽油杆柱底端的柱塞载荷,以拉力为正;qst)为t时刻轴线位置s处抽油杆柱所受的轴向分布载荷,包括:杆柱自重在轴向的分力、液体黏滞阻力、杆管库仑摩擦力在轴向的分量;qubst)为t时刻轴线位置s处抽油杆柱所受的横向分布载荷,包括:杆柱浮重在横向的分量、杆管库仑摩擦力在横向的分量与液体黏滞阻力;mfst)为t时刻轴线s处单位长度抽油杆柱绕杆柱轴线的外力矩,包括:液体黏滞阻力矩与杆管库仑摩擦扭矩;uSt)为任意t时刻悬点相对于上死点的位移,悬点位移向下为正;Fest)为作用于微元两端轴心上的集中力;Mst)为作用于微元两端的弯矩;fest)为作用于微元段上的分布外力;mτst)为单位长度抽油杆柱绕微元段轴线的动量矩。

2.2 抽油杆柱纵横扭耦合振动数学模型

对微元受力进行分析,根据力与力矩的平衡关系得到描述微元段纵向振动、横向振动与扭转振动的动力学平衡方程:

Fes+ds,t-Fes,t+fes,tds=        ρkAkd2Us,tdt2ds mτs,tds=ρkJkd2θs,tdt2ds;        k=1,2,,K (6)

其中:

Us,t=es,tτA+us,tnA+ws,tbAFes,t=Feτs,tτA+Fens,tnA+Febs,tbA=           Faτs,tτC+Fmns,tnC+Fmbs,tbCFeτs,t=Faeτs,t+Fmeτs,tFens,t=Faens,t+Fmens,tFebs,t=Faebs,t+Fmebs,tFes+ds,t=Fes,t+Fes,tsdsFes,ts=Feτs,ts-kAFens,tτA+           Feτs,tkA+Fens,ts-Febs,tTAnA+           Fens,tTA+Febs,tsbAfes,t=feτs,tτA+fens,tnA+febs,tbA (7)

式中 ρk为第k级抽油杆柱密度;Ak为第k级抽油杆柱横截面积;K为抽油杆柱级数;θst)为t时刻抽油杆柱任意截面s处的扭转角;Jk为第k级抽油杆柱绕自身轴线的极惯性矩;Fst)为t时刻抽油杆柱任意截面s处井眼轴线切线方向的内力;Fenst)为t时刻抽油杆柱任意截面s处井眼轴线主法线方向的内力;Febst)为t时刻抽油杆柱任意截面s处井眼轴线副法线方向的内力;Fst)为t时刻抽油杆柱任意截面s处杆柱轴线切线方向的内力;Fmnst)为t时刻抽油杆柱任意截面s处杆柱轴线主法线方向的内力;Fmbst)为t时刻抽油杆柱任意截面s处杆柱轴线副法线方向的内力;Faeτst)为Fst)在井眼轴线切线方向上的分力;Faenst)为Fst)在井眼轴线主法线方向上的分力;Faebst)为Fst)在井眼轴线副法线方向上的分力;Fmeτst)为[Fmnst)+ Fmbst)]在井眼轴线切线方向上的分力;Fmenst)为[Fmnst)+Fmbst)]在井眼轴线主法线方向上的分力;Fmebst)为[Fmnst)+Fmbst)]在井眼轴线副法线方向上的分力;(FsFs求一阶偏导;(FensFens求一阶偏导;(FebsFebs求一阶偏导;fst)为t时刻抽油杆柱任意截面s处井眼轴线切线方向所受外力;fenst)为t时刻抽油杆柱任意截面s处井眼轴线主法线方向所受外力;febst)为t时刻抽油杆柱任意截面s处井眼轴线副法线法向所受外力。

式(6)进行分解与化简,可以得到杆柱在弯曲井眼轨道内的纵横扭耦合振动平衡微分方程:

Feτs,ts-Fens,tkAsw+feτs,t=           ρkAkd2es,tdt2Feτs,tkAsw+Fens,ts-Febs,tTAsw+           fens,t=ρkAkd2us,tdt2Fens,tTAsw+Febs,ts+febs,t=           ρkAkd2ws,tdt2mτs,t=ρkJkd2θs,tdt2;k=1,2,,K (8)

式(8)中,轴向运动抽油杆静位置曲线坐标与动位置曲线坐标的关系为:

sw=s-S-uSt (9)

2.2.1 内力与动量矩计算方法

抽油杆柱长度的变化属于小变形,杆柱变形依然处于弹性变形范围之内,因此由弹性力学基本知识可以得到抽油杆柱轴线切线方向上的内力与其在井眼轴线切线方向、主法线方向与副法线方向分力的表达式:

Faτ=E1A1es-kAuFaττC=FaeττA+FaennA+FaebbA=           FaτXτA+FaτYnA+FaτZbA (10)

式中 E1为第1级抽油杆柱材料弹性模量。

根据微元受力图,对微元段建立如下力矩平衡方程:

Ms+ds,t-Ms,t+rCs+ds,t-rCs,t·Fes+ds,t+12rCs+ds,t-rCs,t·fes,tds-mfcs,tτCds-mfls,tτCds=mτs,tτCds (11)

化简式(11)得到如下公式:

Ms-FmbnC+FmnbC-mfc+mfl+mττC=0 (12)

其中:

M=EkIkkCbC+GkJkθsτCMs=GkJkθssτC+GkJkθskC-           EkIkkCTCnC+EkIkkCsbC (13)

式中 MsMs求一阶偏导;Ek为第k级抽油杆柱材料弹性模量;mfcst)为t时刻抽油杆柱任意截面s处的杆管库仑摩擦力导致的绕轴分布摩擦扭矩;mflst)为t时刻抽油杆柱任意截面s处的液体阻尼力导致的绕轴分布阻尼力矩;Ik为第k级抽油杆柱的抗弯惯性矩;Gk为第k级抽油杆柱材料的剪切弹性模量;θsθs求一阶偏导;θssθs求二阶偏导;(kCskCs求一阶偏导。

式(13)代入式(12),求得FmbFmnmτ的表达式:

mτ=GkJkθss-mfc-mflFmn=-EkIkkCsFmb=GkJkθs-EkIkTCkC (14)

通过坐标转换求得FmeτFmenFmeb的表达式:

Fmeτ=FmnXs-kAY+FmbTAY2+YZs-           kAXZ-YsZ+TAZ2kC-1Fmen=FmnkAX+Ys-TAZ+           FmbXsZ-kAYZ-TAXY-XZskC-1Fmen=FmnTAY+Zs+FmbkAX2+           YsX-TAXZ-XsY+kAY2kC-1 (15)

结合式(7),(10)和(15)即可求得Fe表达式:

Fes,t=Feτs,tτA+Fens,tnA+Febs,tbAFeτ=FaτX+FmeτFen=FaτY+FmenFeb=FaτZ+Fmeb (16)

杆柱各方向内力FFenFeb关于s的偏导均可根据式(16)直接求得,限于篇幅本文不再赘述。

2.2.2 外力与外力矩计算方法

杆柱所受分布外力可以由下式表示:

feτ=-μdedt-ffτ-ρkAkgcosτA,kfen=-μdudt-ffn-ρk-ρlAkgcosnA,kfeb=-μdwdt-ffb-ρk-ρlAkgcosbA,k (17)

式中 ρl为管柱内液体密度;μ为液体黏性阻尼系数;g为重力加速度;f为单位长度抽油杆柱沿井眼轴线切线方向的杆管库仑摩擦力;ffn为单位长度抽油杆柱沿井眼轴线主法线方向的杆管库仑摩擦力;ffb为单位长度抽油杆柱沿井眼轴线副法线方向的杆管库仑摩擦力。

杆管接触与相对运动导致的库仑摩擦力是杆柱受到的外力之一。图3为杆管库仑摩擦力示意图。

图3  杆管库仑摩擦力示意图

Fig.3  Diagram of Coulomb friction force between rod string and tubing

图3(a)为杆柱与管柱空间形态示意图,图中标注出了杆柱变形后的切线方向、油管的切线方向与杆管接触点处油管的切平面。在切平面内建立直角坐标系DPQ,其中DQ与油管切线方向相同,DP为油管切平面与法平面在接触点的交线,DR为杆柱轴线的切线在油管切平面内的投影。图3(b)为切平面内杆管相对速度与摩擦力示意图,图中,vτD为杆柱横向振动导致的杆柱在DP方向上产生的相对运动速度;vω为由于杆柱扭转导致的杆管相对运动速度;vτ为由于杆柱纵向振动导致的杆管相对运动速度;vr为抽油杆柱相对于油管在切平面内的运动速度,为vτDvωvτ的速度矢量和;vrnv分别表示vr在不同方向上的速度分量;vf为杆管接触点在管柱切平面内的相对运动速度;θDn为DP和DR方向的夹角;θfvf和DP的夹角;θmvωvf之间的夹角;ff为抽油杆柱与油管柱间的库仑摩擦力;fffmffr分别代表摩擦力在不同方向上的分量。图3(c)为接触点法平面内杆管相对速度与库仑摩擦力的示意图,图中,vnvb分别为du/dt,dw/dtvτC为抽油杆柱轴线C点处的切向速度;θC为杆管位置方向与vn之间的夹角;ffnffb分别为ffr在井眼轴线主法线方向、副法线方向的分量。

根据图3中所呈现的杆管相对运动关系,可以求得各方向上的杆管库仑摩擦力:

ffτ=μcFtsinθfffr=μcFtcosθfffn=ffrsinθCffb=ffrcosθC (18)

其中:

θf=0.5π,  vrn=0θf=arccosvrnvrτ2+vrn2,  vrn0vrτ=vτ+ωricosθDnvrn=vτCrtrt-ri+ωrisinθDnvτC=-sgnwvnsinθC-sgnuvbcosθCθDn=arccosxtX2+xt2xt=-sgnwYsinθC-sgnuZcosθCθC=0.5π,  u=0θC=arctanabsw/u,  u0i=r,c,s (19)

式中 μc为库仑摩擦阻尼系数;Ft为单位长度抽油杆所受杆管接触力;ω为抽油杆柱扭转振动角速度;当i取r时,表示杆柱节点;当i取c时,表示接箍节点;当i取s时,表示扶正器节点;rr为抽油杆半径;rc为接箍半径;rs为扶正器半径,rt为油管半径。

杆管摩擦力作用在抽油杆柱的表面上,当杆柱的摩擦力方向与DR不重合时,摩擦力ff在垂直于DR的方向上产生摩擦分力ffm。摩擦分力ffm的存在将产生绕杆柱轴线的摩擦扭矩mfc,该绕轴扭矩将作为杆柱扭转振动的外力矩驱动或阻碍杆柱的扭转振动。绕轴摩擦扭矩mfc的表达式如下式所示:

mfc=κmffmri;i=r,c,sffm=μcFtcosθm (20)

其中:

κm=sgncosθDn-sgnvrτcosθfcosθm=1-Xvrτ+xtvrnX2+xt2vrτ2+vrn2 (21)

杆柱在油液中扭转运动时还受到液体阻力矩mfl,其表达式为:

mfl=vωθt (22)

观察式(18),(20)可以发现,绕杆柱轴线的外界分布扭矩受杆柱形态、杆管相对速度与杆管接触力的影响,即受杆柱纵横耦合振动的影响。因此对杆柱纵横耦合振动的仿真计算是实现杆柱扭转振动仿真的基础。同时,杆柱的扭转振动将影响杆柱的内力与杆管摩擦力的方向,即杆柱的扭转振动也是杆柱纵横耦合振动的影响因素。因此,杆柱的动力学模型为具有初弯曲的细长杆柱的纵横扭耦合振动模型。

2.3 边界条件

在研究纵向振动时,因为杆柱的顶端受悬点位移的约束,因此可以将顶端简化为位移的边界条件;杆柱底端受抽油泵的约束,因此可以将底端简化为力的边界条件。杆柱纵向振动的边界条件可表示为:

es=0=uStFaτs=L=Pt (23)

式中 L为抽油杆柱长度。

在研究横向振动时,井口盘根约束了杆柱的横向位移和横截面转角,因此悬点至井口盘根处抽油杆柱可以简化为可轴向滑动的固定约束;杆柱底端的横向位移与横截面转角受到泵筒的约束,因此可以将杆柱底端简化为可滑动的固定约束。抽油杆柱横向振动两端的边界条件可表示为:

u0suSt=0,us=L=0us0suSt=0,uss=L=0w0suSt=0,ws=L=0ws0suSt=0,wss=L=0 (24)

抽油机与抽油杆柱通过钢丝绳与悬绳器连接在一起,两平行钢丝绳上部与驴头弧面贴合,下端与悬绳器连接,悬绳器与抽油杆柱固接在一起。在抽油机运行过程中,杆柱的扭转会带动悬绳器绕其轴线(与抽油杆柱轴线重合,与驴头圆弧面的中心线处于同一平面)旋转,导致两钢丝绳偏离其原始平面,进而使钢丝绳拉力绕悬绳器轴线产生一力矩,此力矩与杆柱顶端的扭转角相关。图4为悬绳器扭矩示意图。图4中,H1H2与H3H4表示两钢丝绳平行时的轴间距;H1H4与H2H3表示悬空段钢丝绳长度(钢丝绳与驴头弧面结合点到悬绳器连接点的长度,随抽油机的运转不断改变);pv表示悬点载荷的投影。抽油杆柱扭转振动顶部边界条件即可表示为:

G1J1θ0,ts=prlt·Lh2sinθ0,t4l0+uSt2-Lhsin0.5θ0,t2 (25)

式中 prlt)为悬点载荷;l0为悬点处于上死点时H1H4或H2H3的长度。

图4  悬绳器扭矩示意图

Fig.4  Schematic diagram of the beam hanger torque

抽油杆柱的底端与抽油泵柱塞相连,可以将底端简化为自由端。杆柱扭转振动的底部边界条件可表示为:

GKJKθL,ts=0 (26)

2.4 连续性条件

对于多级组合抽油杆柱的纵横扭耦合振动,各级杆柱连接点处存在位移相等的连续性条件:

us,ts=Lk-=us,ts=Lk+uns,ts=Lk-=uns,ts=Lk+ubs,ts=Lk-=ubs,ts=Lk+ (27)

式中 Lk为第k级抽油杆柱的长度;unub分别为井眼轴线主法线和副法线方向的位移。

转角相等的连续性条件:

θs,ts=Lk-=θs,ts=Lk+ (28)

轴力(大小和方向)相等的连续性条件:

EkAkscs-1s=Lk-=Ek+1Ak+1scs-1s=Lk+-Ak-Ak+1pckXs=Lk-=Xs=Lk+Ys=Lk-=Ys=Lk+Zs=Lk-=Zs=Lk+ (29)

式中 pck为连接点处的液体压力。

扭矩与弯矩(大小和方向)相等的连续性条件:

GkJkθs,tss=Lk-=Gk+1Jk+1θs,tss=Lk+EkIkkCs,ts=Lk-=Ek+1Ik+1kCs,ts=Lk+Xs-kAYkCs,ts=Lk-=Xs-kAYkCs,ts=Lk+kAX+Ys-TAZkCs,ts=Lk-=kAX+Ys-TAZkCs,ts=Lk+ (30)

剪切力(大小和方向)相等的连续性条件:

EkIkkCs,tTCs,ts=Lk-=Ek+1Ik+1kCs,tTCs,ts=Lk+EkIkkCs,tss=Lk-=Ek+1Ik+1kCs,tss=Lk+ (31)

2.5 初始条件

假设抽油杆柱在初始状态时,杆柱轴线与井眼轴线重合,悬点位于下死点。初始条件可表示为:

θs,0=0,θs,0t=0Us,0=0,Us,0t=0uS0=S (32)

式中 Us,0为初始时刻抽油杆柱任意截面s处的变形量。

2.6 井筒约束

抽油杆柱与油管接触后产生杆管接触力,杆管接触力由杆管碰撞产生的碰撞力Fs和侧向压力Fn组成,其中Fn由杆柱弯曲导致的内力、杆柱弯曲与轴向拉力产生的侧向力组成。

杆管接触条件可以表示为:

u2+w2R-ri;    i=r,c,s (33)

式中 R为油管内圆半径。

碰撞后相应节点的位置为:

ut+=sgnuR-ricosθCt-wt+=sgnwR-risinθCt-;i=r,c,s (34)

碰撞后相应节点的速度为:

vnt+=sgnuvnCcosθC+sgnwvτCsinθCt+vbt+=sgnwvnCsinθC-sgnuvτCcosθCt+ (35)

其中:

vτCt+=-γsvτCt-vnCt+=vnCt-vτCt-=sgnwvnsinθC-sgnuvbcosθCt-vnCt-=sgnuvncosθC+sgnwvbsinθCt-vnt-=utt-vbt-=wtt- (36)

式中 t+表示碰撞后,t-表示碰撞前;γs为碰撞恢复系数,其值取决于碰撞体的材

24

由于较难确定碰撞力的瞬时值,因此本文通过冲量的改变来描述碰撞的强弱。假设在ts

25时间段内杆管完成碰撞,此时间段内由于杆管碰撞产生的分布力可由下式计算:

Fs=-ρAts1+γsvτCt- (37)

杆柱侧向压力的计算公式为:

Fn=κnFeτkA+Fens-FebTA-ρ-ρlgAcosnAsw,kcosθC+κbFenTA+Febs-ρ-ρlgAcosbAsw,ksinθC (38)

杆管接触力可以表示为:

Ft=Fs+Fn (39)

3 纵横扭耦合振动数值仿真模型

将抽油杆柱离散为N个节点,各节点的间距为Δs。为实现对抽油杆柱振动的仿真计算,首先需要对描述杆柱纵横扭耦合振动的二阶偏微分方程进行降阶处理。对扭转振动的二阶偏微分仿真模型进行降阶处理,得到其一阶偏微分形式,如下式所示:

θt=ωωt=1ρkJkGkJk2θs2-mfc-mfl (40)

对纵向振动的二阶偏微分仿真模型进行降阶处理,得到其一阶偏微分形式,如下式所示:

et=vτvτt=1ρkAkFeτs-FenkA+feτ (41)

对横向振动的二阶偏微分仿真模型进行降阶处理,得到其一阶偏微分形式,如下式所示:

ut=vnvnt=1ρkAkFeτkA+Fens-FenTA+fenwt=vbvbt=1ρkAkFenTA+Febs+feb (42)

建立变量数组[y]用以分别代表不同的变量。数组[y]中各元素与变量的对应关系如下式所示:

yi=θi;  i=1,2,3,,NyN+i=ωi;  i=1,2,3,,Ny2N+i=ei;  i=1,2,3,,Ny3N+i=vτ,i;  i=1,2,3,,Ny4N+i=ui;  i=1,2,3,,Ny5N+i=vn,i;  i=1,2,3,,Ny6N+i=wi;  i=1,2,3,,Ny7N+i=vb,i;  i=1,2,3,,N (43)

应用四阶龙格‑库塔法求解数学方程需要知道数组[y]对t的微分形式。数组[y]对t的微分形式中各元素与变量的对应关系如下式所示:

y˙i=ωi;  i=1,2,3,,Ny˙N+i=ωit;  i=1,2,3,,Ny˙2N+i=vτ,i;  i=1,2,3,,Ny˙3N+i=vτ,it;  i=1,2,3,,Ny˙4N+i=vn,i;  i=1,2,3,,Ny˙5N+i=vn,it;  i=1,2,3,,Ny˙6N+i=vb,i;  i=1,2,3,,Ny˙7N+i=vb,it;  i=1,2,3,,N (44)

式中 θeuws的偏微分项可通过牛顿中心差分确定。式中大部分杆柱节点所对应的变量可以由式(40)~(42)确定。

杆柱上第一个仿真节点与最后一个仿真节点的动力学行为受边界条件的影响,第一个仿真节点在计算θs的偏导时采用前差公式,最后一个节点在计算θes的偏导时采用后差公式。根据顶部边界条件得到第一个仿真节点的运动行为,如下式所示:

y˙N+1=G1J1Δsθ2s-θ1s-mfc-mflρ1J1y2N+1=uSty4N+1=0y6N+1=0 (45)

根据杆柱振动的底部边界条件得到最后一个仿真节点的运动行为:

y˙2N=GKJKΔsθNs-θN-1s-mfc-mflρKJKy˙4N=PpΔs-EKAKΔseN-1s-FenkA+feτAKρKy5N=0y7N=0 (46)

式中 Pp为杆柱底端所受柱塞载荷。

多级抽油杆柱的连接点ik的位移和转角与其前、后各一个仿真节点(ik-1,ik+1)的位移可以由模型的连续性条件确定。通过泰勒级数展开求得连接点处轴向位移、横向位移与扭转角对s的偏导数。将连接点处偏导数的表达式代入式(27)~(31)即可得到K-1组10个关于θikeik-1eikeik+1uik-1uikuik+1wik-1wikwik+1的等式组成的方程组。求解方程组即可得到连接点的位移和转角与其前、后各一个仿真节点的位移。

最后根据式(40)~(46)即可完成对抽油杆柱纵横扭耦合振动的仿真计算。仿真计算流程图如图5所示。

图5  仿真计算流程图

Fig.5  Simulation calculation flowchart

4 仿真分析

4.1 基本参数

本文算例仿真参数如表1所示。油井全井布置扶正器,每根杆上均匀布置两个扶正器。查找油井井斜数据得到井眼轨道示意图如图6所示。

表1  算例仿真参数
Tab.1  Simulation parameters of the example
参数型号/数值参数型号/数值
抽油机型号 CYJ10-3-53HB 含水率/% 19
电机型号 Y280S-6 动液面/m 1570
l0/m 0.5 泵径/mm 44
Lh/m 0.21 油管内径/mm 62
下泵深度/m 2100 碰撞恢复系数γs 0.56
冲程/m 3 碰撞时间ts /s 0.03
冲次/min-1 4 库仑摩擦系数 0.07
油压/MPa 0.3 抽油杆材质 20CrMo
套压/MPa 0.5 钢密度/(kg·m‒3) 7800
气油比/(m3·m‒3) 21 弹性模量/GPa 210
杆柱组合/(mm×m)

25×400+22×750+

19×950

各级杆柱对应接箍外径/mm 55, 46, 41.3
各级杆柱对应扶正器外径/mm 58,58,58

图6  井眼轨道示意图

Fig.6  Diagram of wellbore trajectory

4.2 纵向振动仿真结果

分别采用本文建立的抽油杆柱纵横扭耦合振动仿真模型与文献[

29]建立的定向井抽油杆柱纵向振动仿真模型对杆柱的振动进行仿真,得到了如图7(图中悬点位移以下死点为0点,向上为坐标轴正方向)和8所示的杆柱纵向振动仿真结果。

图7  悬点示功图仿真结果

Fig.7  Simulation results of suspension indicator diagram

图8  泵功图仿真结果

Fig.8  Simulation results of pump power diagram

观察杆柱纵向振动仿真结果发现:(1)杆柱纵横扭耦合振动仿真模型所得杆柱纵向振动仿真结果与杆柱纵向振动仿真模型所得杆柱纵向振动仿真结果近似。(2)与杆柱纵向振动仿真结果相比,在抽油机的下冲程纵横扭耦合振动所得到的悬点载荷更小。

以上现象说明:(1)杆柱的纵横扭耦合振动对杆柱纵向振动的仿真结果有影响,但影响不明显。(2)考虑杆柱的纵横扭耦合振动后杆柱的下行阻力增大,进而导致下冲程悬点载荷降低。

4.3 横向振动仿真结果

分别采用本文建立的抽油杆柱纵横扭耦合振动仿真模型与文献[

30]建立的定向井抽油杆柱横向振动仿真模型对杆柱的振动进行仿真,得到如图910所示的相同时刻杆柱横向振动仿真结果对比。

图9  13.5 s时横向位移仿真结果

Fig.9  Simulation results of the transverse displacementat 13.5 s

图10  13.5 s时杆管接触力仿真结果

Fig.10  Simulation results of the contact force betwen rod string and tubing at 13.5 s

观察杆柱横向振动仿真结果对比发现:(1)相同时刻两种模型的杆柱横向位移仿真结果有明显区别,在杆柱的受压段(约1800~2100 m)尤为不同。(2)相同时刻两种模型的杆管接触力的仿真结果也有明显区别。由于考虑纵横扭耦合振动后杆柱下行阻力更大,因此在杆柱的造斜段(约400~800 m)纵横扭耦合振动仿真模型所得杆管接触力仿真值相对较小,在杆柱的受压段,两种模型所得杆管接触力仿真值明显不同。

以上现象说明:杆柱的纵横扭耦合振动对杆柱的横向振动仿真结果有明显影响,在进行横向振动仿真时有必要考虑纵横扭耦合振动的影响。

图11为一完整冲程不同节点杆柱横向振动仿真结果;图12为杆管最大接触力沿杆长的分布曲线。

图11  单次全冲程不同节点横向振动仿真结果

Fig.11  The simulation results of the rod string transversevibration of different nodes in a complete stroke

图12  最大杆管接触力仿真结果

Fig.12  Simulation results of maximum contact force between rod string and tubing

观察图1112仿真结果可以发现:(1)在杆柱的受拉段,接箍节点与扶正器节点均存在与油管的接触(接触力不为0),杆体节点不与油管接触(接触力为0)。在杆柱的受压段,接箍节点、扶正器节点与杆体节点均与油管接触。(2)油井造斜段的杆管接触力较杆柱受压段的接触力大;杆柱受压段的横向振动与杆管碰撞较油井造斜段剧烈。(3)相较于杆体节点与接箍节点,扶正器节点的杆管接触力更大。

以上现象说明:(1)油井的造斜段与杆柱的受压段为杆管偏磨的危险点。(2)扶正器有降低杆管偏磨现象的作用。

4.4 扭转振动仿真结果

采用本文建立的纵横扭耦合振动仿真模型对杆柱进行仿真,得到如图13所示的扭转振动仿真结果。

  

  

图13  杆柱扭转振动仿真结果

Fig.13  Simulation results of the rod string torsional vibration

根据图13发现:(1)杆柱上下冲程转化时杆柱轴向速度较慢,杆管库仑摩擦力在垂直于轴线方向上的分量更大,因此杆柱扭转角在上下冲程转化时变化最为明显。(2)整体来说杆柱扭转振动的幅值沿杆长逐渐增大。在油井的造斜段由于杆管接触力较大,因此杆柱扭转振动的幅值在此部分有突变;最后一级杆柱杆径最小、扭转刚度最低,且杆柱受压段杆管变形与横向振动更明显,因此最后一级杆柱扭转振动幅值增幅显著。(3)算例油井杆柱扭转振动所产生的扭矩远小于接箍的上扣预紧扭矩,即正常预紧情况下杆柱的扭转振动不足以导致抽油杆柱的脱扣现象。但在实际油井中,抽油杆柱轴向受较大交变轴向载荷作用,在循环载荷的作用下螺纹连接的预紧力有可能逐渐降低,当其预紧力降低到杆柱的最大扭矩后即可能导致抽油杆柱产生脱扣现象。

5 结 论

(1)指出了弯曲变形抽油杆柱与油管之间的摩擦力是影响杆柱横向振动的一项主要因素;杆管摩擦力导致的绕轴扭矩是抽油杆柱扭转振动的根本原因。考虑杆柱纵、横、扭振动之间的相互影响,建立了抽油杆柱的纵横扭耦合振动仿真模型。采用四阶龙格‑库塔法实现对仿真模型的求解。

(2)仿真结果表明:杆柱的耦合振动对杆柱纵向振动、横向振动仿真结果均有影响。油井造斜段的杆管接触力相对较大,杆柱受压段的横向振动更为剧烈,因此油井的造斜段与杆柱的受压段均为杆管偏磨的危险点。定向井中,正常预紧情况下杆柱的扭转振动不足以导致抽油杆柱的脱扣现象,但当其螺纹预紧力降低到杆柱的最大扭矩后即可能导致抽油杆柱脱扣。

(3)本文研究工作指出了抽油杆柱扭转振动的机理,实现了对定向井抽油杆柱纵横扭耦合振动的仿真求解,对指导油田生产具有理论与实际意义。除了杆柱偏磨外,杆柱的疲劳破坏也是油井常见故障,因此下一步有必要基于本文仿真结果分析抽油杆柱在多轴应力下的疲劳强度。

参考文献

1

刘合郝忠献王连刚. 人工举升技术现状与发展趋势[J]. 石油学报20153611):1441-1448. [百度学术] 

LIU HeHAO ZhongxianWANG Lianganget al. Current technical status and development trend of artificial lift[J]. Acta Petrolei Sinica20153611): 14411448. [百度学术] 

2

GIBBS S G. Predicting the behavior of sucker-rod pumping system[J]. Journal of Petroleum Technology1963157):769-778. [百度学术] 

3

WANG D YLIU H Z. Prediction and analysis of polished rod dynamometer card in sucker rod pumping system with wear[J]. Shock and Vibration201820184979405. [百度学术] 

4

MORENO G AGARRIZ A E. Sucker rod string dynamics in deviated wells[J]. Journal of Petroleum Science and Engineering2020184106534. [百度学术] 

5

LI QCHEN BHUANG Z Qet al. Study on equivalent viscous damping coefficient of sucker rod based on the principle of equal friction loss[J]. Mathematical Problems in Engineering:Theory, Methods and Applications201920199272751. [百度学术] 

6

WANG X BLÜ L GLI Set al. Longitudinal vibration analysis of sucker rod based on a simplified thermo-solid model[J]. Journal of Petroleum Science and Engineering2021196107951. [百度学术] 

7

YIN J JSUN DYANG Y S. Predicting multi-tapered sucker-rod pumping systems with the analytical solution[J]. Journal of Petroleum Science and Engineering2021197108115. [百度学术] 

8

孙秀荣江少波黄有泉. 直井抽油杆柱在屈曲位移激励下横向振动的仿真研究[J].振动工程学报2018315):854-861. [百度学术] 

SUN XiurongJIANG ShaoboHUANG Youquanet al. The simulation research of sucker rod string transverse vibration under buckling deformation excitation in vertical wells[J]. Journal of Vibration Engineering2018315):854-861. [百度学术] 

9

LUKASIEWICZ S A. Dynamic behavior of the sucker rod string in the inclined well[C]//SPE Production Operations Symposium. Oklahoma1991. [百度学术] 

10

XU J. A new approach to the analysis of deviated rod-pumped wells[C]//SPE International Petroleum Conference & Exhibition of Mexico. Veracruz, Mexico1994. [百度学术] 

11

XU JDOTY D RBLAIS Ret al. A comprehensive rod-pumping model and its application to vertical and deviated wells[C]//SPE Mid-Continent Operations Symposium. Oklahoma1999. [百度学术] 

12

王宏博董世民. 轴向往复运动抽油杆柱在弯曲井眼内横向振动的仿真模型[J]. 工程力学20203710):228-237. [百度学术] 

WANG HongboDONG Shimin. A model for the transverse vibration simulation of sucker rod strings with axial reciprocating motion in curved wellbores[J]. Engineering Mechanics20203710):228-237. [百度学术] 

13

WANG H BDONG S M. Spatial transverse vibration simulation model of axially moving sucker rod string under the excitation of curved borehole[J]. Shock and Vibration202020203108718. [百度学术] 

14

王文昌狄勤丰姚建林. 三维定向井抽油杆柱力学特性有限元分析新方法[J]. 石油学报2010316): 1018-1023. [百度学术] 

WANG WenchangDI QinfengYAO Jianlinet al. A new finite element method in analyzing mechanical properties of sucker rod strings in 3D directional wellbores[J]. Acta Petrolei Sinica2010316):1018-1023. [百度学术] 

15

JOHANCSIK C AFRIESEN D BDAWSON R. Torque and drag in directional wells-prediction and measurement[J]. Journal of Petroleum Technology1984366): 987-992. [百度学术] 

16

HO H S. An improved modeling program for computing the torque and drag in directional and deep wells[C]//SPE Annual Technical Conference and Exhibition. Houston, Texas1988407-418. [百度学术] 

17

李子丰李相方单志刚. 海底大位移井井眼净化程度和卡钻可能性监测[J]. 石油钻探技术2000285): 4-6. [百度学术] 

LI ZifengLI XiangfangSHAN Zhigang. Monitoring wellbore cleaning degree and sticking possibility in offshore extended wells[J]. Petroleum Drilling Techniques2000285): 4-6. [百度学术] 

18

张学鸿陈抡元刘巨保. 整体钻柱力学接触有限元分析[J]. 石油学报1992133): 102-108. [百度学术] 

ZHANG XuehongCHEN LunyuanLIU Jubao. Finite element analysis of the contact and friction of a whole drillpipe string[J]. Acta Petrolei Sinica1992133): 102-108. [百度学术] 

19

帅健吕英民蔡强康. 全井中钻柱的有限元模型及应用[J]. 石油学报1995161): 118-126. [百度学术] 

SHUAI JianLÜ YingminCAI Qiangkang. Finite element model of the drill string in a whole well bore and its application[J]. Acta Petrolei Sinica1995161): 118-126. [百度学术] 

20

郭永峰金晓剑谭树人. 非线性有限元法预测近海水平井钻进摩阻[J]. 中国造船2003442): 61-66. [百度学术] 

GUO YongfengJIN XiaojianTAN Shuren. Prediction of drilling drag and torque in horizontal wells of offshore oilfield by nonlinear FEM[J]. Shipbuilding of China2003442): 61-66. [百度学术] 

21

付建红龚龙祥胡顺渠. 基于ANSYS的水平井下套管摩阻分析计算[J]. 石油钻采工艺2007294): 32-35. [百度学术] 

FU JianhongGONG LongxiangHU Shunquet al. Calculation of frictional drag of casing running in horizontal well based on ANSYS[J]. Oil Drilling & Production Technology2007294): 32-35. [百度学术] 

22

祝效华童华刘清友. 基于钻柱系统动力学的摩阻扭矩计算新方法[J]. 系统仿真学报20072119): 4853-4856. [百度学术] 

ZHU XiaohuaTONG HuaLIU Qingyouet al. New method on calculation of torque and drag based on drilling string system dynamics[J]. Journal of System Simulation20072119): 4853-4856. [百度学术] 

23

ZHU X HLI BLIU Q Yet al. New analysis theory and method for drag and torque based on full-hole system dynamics in highly deviated well[J]. Mathematical Problems in Engineering20152015535830. [百度学术] 

24

祝效华童华刘清友. 旋转钻柱与井壁的碰撞摩擦边界问题研究[J]. 中国机械工程20171815):1833-1837. [百度学术] 

ZHU XiaohuaTONG HuaLIU Qingyouet al. Research on the dynamic boundary condition between revolving drill string and borehole wall[J]. China Mechanical Engineering20171815):1833-1837. [百度学术] 

25

FIRU L SCHELU TMILITARU-PETRE C. A modern approach to the optimum design of sucker-rod pumping system[C]// SPE Annual Technical Conference and Exhibition. Denver, Colorado20031-9. [百度学术] 

26

VOLPI L PLOBO D MRITTO T G. A Stochastic analysis of the coupled lateral-torsional drill string vibration[J]. Nonlinear Dynamic202110349-62. [百度学术] 

27

王宏博董世民张洋. 基于气液分离的泵内压力与抽油杆柱纵向振动耦合仿真[J]. 石油学报2023442): 394-404. [百度学术] 

WANG HongboDONG ShiminZHANG Yanget al. Coupling simulation of the pressure in pump and the longitudinal vibration of sucker rod string based on gas liquid separation[J]. Acta Petrolei Sinica2023442): 394-404. [百度学术] 

28

董世民姚春东. 抽油杆柱纵向振动共振条件的质疑与力学模型的修正[J]. 振动工程学报2003163): 135-138. [百度学术] 

DONG ShiminYAO Chundong. Doubting the resonance condition of the longitudinal vibration of the sucker rod string introduced by APIRP11L and improving of the mechanical modeling[J]. Journal of Vibration Engineering2003163): 135-138. [百度学术] 

29

董世民李宝生. 水平井有杆抽油系统设计[M]. 北京石油工业出版社199619-21. [百度学术] 

30

董世民王宏博. 定向井抽油杆柱横向振动仿真模型及扶正器布点优化[J].石油学报20204112):1686-1696. [百度学术] 

DONG ShiminWANG Hongbo. Simulation model of lateral vibration of sucker rod string in directional wells and point arrangement optimization of centralizer[J]. Acta Petrolei Sinica20204112):1686-1696. [百度学术]