基于电动机功率间接测量游梁式抽油机悬点力的方法与流程

本发明是一种全新的游梁式抽油机悬点运动测力方法，具体地是以曲柄转角(θ)和电动机功率(p)为自变量而间接测量获取悬点的实时受力值变化，属于石化行业油田施工生产与计算机信息化领域。

背景技术：

目前国内油田勘探施工领域，较普遍地采用游梁式抽油机。现有技术针对悬点运动参数的求解，基本上是通过示功图或电功图以揭示机井光杆的悬点载荷与冲程之间的变化关系。

如以下在先公开文献，如《石油学报》第19卷第2期第107-110页，游梁式抽油机运动参数的精确解，1998年4月刊。通过悬点运动参数建立游梁摆动方程，从而得出悬点运动的位移、速度和加速度计算公式，最终给出求游梁式抽油机悬点运动参数的精确解。

再如《石油机械》第34卷第5期第22-24页，游梁式抽油机运动规律的多体动力学分析，2006年3月刊。采用多体动力学方法进行分析，以动力学方程值积分的结果对位移和速度级约束方程进行求解和违约修正。

上述现有技术实现的测力方法，计算过程复杂且具有多维的特点，计算结果因此也就累积较大的误差，仍非真正意义上的精确解。另外，计算过程过于复杂而具有较低的实地使用参考价值。

有鉴于此，特提出本专利申请。

技术实现要素：

本发明基于电动机功率间接测量游梁式抽油机悬点力的方法，在于解决现有技术存在的问题而以游梁式抽油机的曲柄转角(θ)、电动机功率(p)为自变量，基于游梁式抽油机四连杆机构的连接与运动关系，建立起悬点的实时参数动态模型，即以曲柄转角的变化可以实时地计算出悬点位置、速度与加速度的对应变化，加之电动机功率(p)、曲柄转角(θ)以间接地测量出悬点实时力的大小，从而有效地提高油井监控的实时性与可靠性、降低设计与生产成本。

为实现上述设计目的，所述基于电动机功率间接测量游梁式抽油机悬点力的方法，在游梁式抽油机四连杆机构中，悬点动态受力符合以下表达式：

即在已知电动机功率(p)和曲柄转角(θ)的前提下，间接地测量求得悬点的实时动态受力值。

如上所述，基于电动机功率间接测量游梁式抽油机悬点力的方法具有的优点是：

1、实现了一种新的间接式测力方法，相较于现有技术具有实时性强、可靠性高、成本低、准确性高的特点。

2、在此基础上能显著地提高油井状态监控的实时性，及时反映悬点载荷、位移、加速度的动态关系，在游梁式抽油机自动化、信息化升级改造过程中具有使实用价值。

附图说明

图1是游梁式抽油机四连杆机构的工作原理图；

图2是电动机输入转矩与减速机输出扭矩之间的关系图；

图3是基于本申请和抽油机实际参数所获取悬点运动位移变化周期的示意图；

图4是基于本申请和抽油机实际参数所获取悬点运动速度变化周期的示意图；

图5是基于本申请和抽油机实际参数所获取悬点运动加速度变化周期的示意图。

具体实施方式

下面结合附图和实施例对本发明作进一步说明。

实施例1，一种基于电动机功率间接测量游梁式抽油机悬点力的方法，是已知游梁式抽油机电动机功率p数值的前提下，通过依次计算出电动机输入功率与转矩、减速机输出扭矩之间的关系，通过曲柄转角θ间接地测量出在四连杆机构运动模型中悬点实时的受力变化，即图1中悬点c受到的实时拉力变化情况。

如图1所示，游梁式抽油机四连杆机构的直角坐标运动模型原理示意图。

在四连杆机构中，o点为游梁式抽油机支架轴的支点，o1点为游梁式抽油机动力输入轴的中点，c点为游梁式抽油机的悬点，b点为游梁式抽油机横梁轴(游梁后臂)的端点，a点为游梁式抽油机曲柄的端点。

工作过程是，电动机的动力输入轴通过曲柄将力矩传递到a点，a点绕o1点以半径为r、角速度为ω做圆周运动，a点与b点通过连杆p连接，使b点以半径为m在a点的带动下做圆弧往复运动；c点与b点为同一梁(游梁)的两端，以o点为支点做圆弧往复运动，且c点的运动方向与b点相反。

其中，r是曲柄长度(单位：mm)；p是连杆长度(单位：mm)；i为支架轴中点到动力输入轴中点的水平距离(单位：mm)；f为支架轴中点到动力输入轴中点的垂直距离(单位：mm)；m为支架轴中点o到横梁轴端点b的长度(单位：mm)；n为支架轴中点o到悬点c的长度(单位：mm)。

曲柄的转角θ，是在电动机通过减速机出扭矩的驱动下，在直角坐标系中的实时转角。

具体的计算过程如下：

设定电动机的输出功率p与电压u、电流i之间的关系表达式如下：

p＝ui(1)

其中，功率p的单位为w，电压u的单位为v，电流i的单位为a。

电动机的输出转矩t，是使抽油机曲柄转动的力矩，简称转矩。转矩t与功率p、转速n的关系表达式如下：t＝9550p/n(2)

由此也可推导出：p＝t·n/9550(3)

其中，功率p的单位为kw，转速n的单位为r/min，转矩t的单位为n·m，9550是计算得出的常数值系数。

以下是常数值系数9550的推到过程：

通常可知，功率p＝扭力f*线速度v(1-1)

由转矩t＝扭力f*作用半径r，即t＝f*r，推出f＝t/r(1-2)

而线速度v＝2πr*每秒转速(n秒)＝2πr*每分转速(n分)/60＝πr*n分/30。

即得出，v＝πr*n/30(1-3)

将公式(1-2)、(1-3)代入公式(1-1)得：

p＝f*v＝t/r*πr*n/30＝π*t*n/30

如果将p的单位换成kw，那么就是如下公式：p*1000＝π/30*t*n

即，30000/3.1415926*p＝t*n

9549.297*p＝t*n→t＝9550p/n

由以上公式(1)、(3)，电功率p(单位为kw)可由下式表示：

p＝tn/9550＝ui/1000(4)

可变形为tn/9.55＝ui(5)

也可变形为t＝9.55ui/n(6)

如图2所示的是电动机输入转矩与减速机输出扭矩之间的关系示意图。

电动机与减速机之间角速度的比值为传动比，也称速比。

即，ε＝n/n1(7)

式中ε为传动比，无量纲；n为电动机向减速机输入的转数，单位r/min；n1为减速机输出的转数，单位r/min；

减速机输入扭矩t与输出扭矩t1具有如下关系：

即减速机输入轴的扭矩t与其转速n的乘积，与输出轴的扭矩t1与其转速n1的乘积是相等的。即，t·n＝t1·n1→t1＝εt(8)

式中，ε为传动比，无量纲；t为电动机输入扭矩，单位n·m；t1为减速机输出扭矩，单位n·m；

如图1所示，以o点为零点建立直角坐标系，即o点直角坐标为(0,0)，则o1点的坐标为o1(-i,-f)，a点以o1为圆心作圆周运动，运动角速度为ω。则有，(xa+i)²+(ya+f)²＝r²(9)

θ＝ωt(10)

xa＝rcosθ-i(11)

ya＝rsinθ-f(12)

则b点的运动方程：xb²+yb²＝m²(13)

由a到b点的距离公式得：

(xb-xa)²+(yb-ya)²＝p²(14)

将式(14)展开：

xb²+xa²-2xbxa+yb²+ya²-2y^by^a＝p²(15)

式(15)﹣(9)﹣(13)得到：

-2xbxa-2ixa-2ybya-2fya＝p²-m²-r²+f²+i²

2xbxa+2ybya＝m²+r²-p²-f²-i²-2ixa-2fya

将式(16)代入式(13)中：

设q＝(m²+r²-p²-f²-i²-2(xa-2fya)/2，则有

由一元二次方程根与分数的关系解得：

由图中可知，|yb|≤r，即-r≤yb≤r(18)

由式(9)得：

xa²+ya²＝r²-2ixa-2fya(19)

将(19)代入式(17)：

依据线性方程，悬点c位置：

如图3所示的悬点c位置曲线图，通过悬点c位置坐标值yc的变化，可求出悬点c的速度vc，即vc是yc的一阶导数：

据此可绘制出抽油机悬点c基于曲柄转角θ的速度曲线图，如图4所示。

基于相同的原理，由悬点c的坐标值yc得到悬点c的加速度ac，即ac是yc的二阶导数：

据此绘制抽油机悬点(c)基于曲柄转角θ的加速度曲线图，如图5所示。

如图1所示，在a点的力分解为以下表达式：

fx＝facosθ(24)

fy＝fasinθ(25)

悬点c悬挂的总重量w，则力矩tc

tc＝w(g+a)n(26)

同理，在b点的力矩为：

由以上公式(8)可知，在a点的力矩符合以下表达式：

ta＝εt

其中t为电动机的转矩；则在a点的受力是：

fa＝ta/r(28)

由于b点与c点，是以o为支点且力矩达到平衡的，则有：

tb＝tc(29)

即：

因此，悬点c实时动态受到的力可表示为：

将以上公式(8)和(6)t＝9.55p/n代入公式(30)中则悬点c实时受力是：

其中，xb由公式(16)计算得出，则p为电动机输入功率，θ为减速机驱动曲柄的转角，n为电动机转数，电动机功率单位为kw。

通过以上公式(31)可知，在已知电动机输入功率p的基础上，再通过曲柄的转角θ应可间接地测量出悬点c实时受到拉力的变化情况。

如上所述，结合附图和描述给出的方案内容，可以衍生出类似的技术方案。但凡是未脱离本发明的结构的方案内容，均仍属于本发明技术方案的权利范围。