抽油机设计计算

抽油机设计计算第1页，共54页，2023年，2月20日，星期一1.1简介游梁式抽油机主要由游梁—连杆—曲柄机构、减速箱、动力设备和辅助装备等四大部分组成。工作时，电动机的传动经变速箱、曲柄连杆机构变成驴头的上下运动，驴头经光杆、抽油杆带动井下深井泵的柱塞作上下运动，从而不断地把井中的原油抽出井筒。游梁式抽油机具有性能可靠、结构简单、操作维修方便等特点。技术参数符合中华人民共和国行业标准SY/T5044《游梁式抽油机》和美国石油协会API标准，技术成熟。主要特点：

1、整机结构合理、工作平稳、噪音小、操作维护方便；

2、游梁选用箱式或工字钢结构，强度高、刚性好、承载能力大；

3、减速器采用人字型渐开线或双圆弧齿形齿轮，加工精度高、承载能力强，使用寿命长；

4、驴头可采用上翻、上挂或侧转三种形式之一；

5、刹车采用外抱式结构，配有保险装置，操作灵活、制动迅速、安全可靠；

6、底座采用地脚螺栓连接或压杠连接两种方式之一。第2页，共54页，2023年，2月20日，星期一分类常规游梁式抽油机

以常规曲柄平衡抽油机为例，结构示意图如下图所示。

图中，1—底座；2—支架；3—悬绳器；4—驴头；5—游梁；6—横梁轴承座；7—横梁；8—连杆；9—曲柄销装置；10—曲柄装置；11—减速器；12—刹车保险装置；13—刹车装置；14—电动机；15—配电箱。

游梁式抽油机是一种变形的四连杆机构，其整机结构特点像一架天平，一端是抽油载荷，另一端是平衡配重载荷。对于支架来说，如果抽油载荷和平衡载荷形成的扭矩相等或变化一致，那么用很小的动力就可以使抽油机连续不间断地工作。

第3页，共54页，2023年，2月20日，星期一异形游梁式抽油机

异形游梁式系列抽油机包括双驴头抽油机、弯游梁抽油机和对常规游梁抽油机节能改造的机型等机型。异形游梁式抽油机是以常规抽油机为基础模式而研制出的新机型，它采用变径圆弧状的游梁后臂，游梁与横梁之间采用柔性件连接结构，在主要结构上具有常规游梁式抽油机简单，牢靠，耐用等特点，在性能上易于实现长冲程，并且具有突出的节能特点。平衡方式

1．游梁平衡：在游梁的尾部装设一定重量的平衡板，这是一种简单的平衡方式，适用于3吨以下的轻型抽油机。第4页，共54页，2023年，2月20日，星期一2．曲柄平衡：这是一种在油田上常用的平衡方式。顾名思义是将平衡块装在曲柄上，适用于重型抽油机。这种平衡方式减少了游梁平衡引起的抽油机摆动，调整比较方便，但是，曲柄上有很大的负荷和离心力。

3．复合平衡：在一台抽油机上同时使用游梁平衡和曲柄平衡。特点：小范围调整时，可以调整游梁平衡：大范围调整时，则调整曲柄平衡。这种平衡方式适用于中深井。4．气动平衡：利用气体的可压缩性来储存和释放能量达到平衡的目的，可用于10吨以上重型抽油机。这种平衡方式减少了抽油机的动负荷及震动，但其装置精度要求高，加工复杂。第5页，共54页，2023年，2月20日，星期一1.2我国石油开采技术概况

目前，我国石油开采技术仍落后西方。中国油田的种类比较复杂，主要有中高渗透多层砂岩油田、复杂断块油田，低渗透砂岩油田以及稠油油田。根据不同的地质条件和油田类型，中国形成了具有特色的油气开发技术。三元化学负荷驳对提高中国可采储量具有明显的效益，约可以提高5亿动左右石油开采量。在大庆、胜利等油田进行的先导性试验和矿场实验，都取得了良好的效果，原油采收率平均可以提高15%—20%。此外，中国对其他提高原油采收率的新技术，如微生物采油、纳米膜驱油等也进行了研究与试验。尽管中国石油开采技术经过50多年开展已经取得很大发展，但是同西方先进国家相比，总体上应该还有5至10年的差距。而且从具体的技术上看，完全属于中国原始创新的技术比较少，一部分具有中国知识产权的国际先进技术，如聚合物驱、化学复合驱技术等，最早也是在模仿学习国外理论和方法上发展起来的。第6页，共54页，2023年，2月20日，星期一2.1系统总成第7页，共54页，2023年，2月20日，星期一游梁式抽油系统示意图1－电动机；2－皮带减速箱传动装置；3－抽油机主体机构；4-悬绳器；5－抽油杆柱；6－油管柱；7－抽油泵第8页，共54页，2023年，2月20日，星期一2.2电动机已知参数及机械特性(1)电动机参数

第9页，共54页，2023年，2月20日，星期一电动机机械特性

第10页，共54页，2023年，2月20日，星期一电机特性曲线Md随n的变化Med=132.789*Md;>>n=0:0.005:1000;Md=(2*2*360*0.0847*104.7198*(104.7198-pi*n/30))./(0.0847^2*104.7198^2+(104.7198-pi*n/30).^2);第11页，共54页，2023年，2月20日，星期一2.3抽油机已知参数(1)抽油机型号根据抽油机型号，可以确定抽油机传动系统的机构尺寸、运动件质量、转动惯量、质心位置等参数。具体参数包括：R——曲柄半径，m；P——连杆长度，m；C——游梁后臂长度，m；K——基杆长度，m；A——游梁前臂长度，m；I——基杆的水平投影，m；τ——曲柄平衡重滞后角，rad；JB1——小皮带轮转动惯量，kg.m2；JB2——大皮带轮转动惯量，kg.m2；JR1——二级减速箱高速轴（输入轴）系转动惯量，kg.m2；第12页，共54页，2023年，2月20日，星期一JR2——二级减速箱中间轴系转动惯量，kg.m2；JR3——二级减速箱底速轴（输出轴）系转动惯量，kg.m2；iB——皮带传动比；iR1——二级减速箱的第一级传动比；iR2——二级减速箱的第二级传动比；mQ——单曲柄质量，kg；

rQ——曲柄质心半径，m；MC——单侧曲柄平衡块质量，kg；

RC——曲柄平衡块质心半径，m；JC——单侧曲柄平衡块相对其质心的转动惯量，kg.m2；ML——连杆质量，kg；JL——连杆相对其质心的转动惯量，kg.m2；JB——游梁组建相对其回转中心的转动惯量，kg.m2。第13页，共54页，2023年，2月20日，星期一(2)油井参数油井参数包括油藏参数与生产参数。根据假设，和数学建模相关的几个参数包括：——抽油杆柱质量，kg；——抽油杆柱在液柱中的重量，N；——油井动液面以上整个柱塞面积上的液柱重量，N；——抽油杆柱在液体载荷

作用下的静变形，m。各参数具体数值可取：R＝1.060m；P＝3.390m；C＝2.820m；K=4.280m；A＝3.935m；I＝2.82m；第14页，共54页，2023年，2月20日，星期一τ＝0；JB1＝2.0kg.m2；JB2＝15.0kg.m2；JR1＝3.0kg.m2；JR2＝8.0kg.m2；JR3＝15.0kg.m2；iB=5.0；iR1=5.5；iR2=5.45；mQ=500kg；

rQ=1.0m；MC==1000kg；

RC=1.2m；JC=100kg.m2；ML=300kg；JL==300kg.m2；JB=1000kg.m2;MR=4000kg；=30000N；=36000N；=0.2m。第15页，共54页，2023年，2月20日，星期一3.1假设条件为了对异相曲柄平衡抽油机的优化设计，以放映抽油机动力性能的最大扭矩因素为最小目标函数优化抽油机的机构尺寸，将整个抽油机装置简化为以曲柄为等效构件的运动情况，我们做了以下假设条件：第16页，共54页，2023年，2月20日，星期一1、电网供电电压与供电频率是常数；2、不考虑电动机转子到减速箱曲柄轴各传动副的间隙与传动件的弹性变形；3、在分析抽油机主体机构的运动规律时，不考虑各传动副的间隙与传动件的弹性变形；4、悬绳与悬绳器在垂直方向上的等效弹簧常数为定值；5、油管全长锚定或不锚定，当油管不锚定时，考虑液体静负荷引起的油管柱弹性变形；6、不考虑油管内液柱的振动；7、假设抽油杆柱与油管同心；8、油井是铅直的；9、连杆与游梁之间的轴承座质量忽略；10、游梁与支架之间的轴承座质量忽略；11、所有构件为刚性构件；12、曲柄为等效构件；13、求等效力矩时不考虑全部构件的质量；14、求等效转动惯量时考虑连杆质量；15、悬点载荷简化为F16、减速箱效率为0.8806第17页，共54页，2023年，2月20日，星期一第18页，共54页，2023年，2月20日，星期一将图中各杆件看成矢量，可以得到如下的矢量方程

上述矢量方程可以用复变矢量表示为将上式两边对时间t求导，可得到各杆件的角速度为：第19页，共54页，2023年，2月20日，星期一将上式对时间求导，可得各杆件运动的角加速度：

第20页，共54页，2023年，2月20日，星期一

图为Theta2随theta变化曲线theta=0:0.005:2*pitheta2=2.219*pi-theta第21页，共54页，2023年，2月20日，星期一

图为Theta3随theta变化曲线beta=asin((1.064./L).*sin(2.219*pi-theta));theta3=acos((1.2043+L.^2)./(6.052*L))-beta;第22页，共54页，2023年，2月20日，星期一

图为Theta4随theta变化曲线theta4=acos((1.2043-L.^2)./(6.052*L))-beta;第23页，共54页，2023年，2月20日，星期一（2）下面讨论悬点位移，速度与加速度的计算公式：A.悬点中程长度SPR:

第24页，共54页，2023年，2月20日，星期一B.悬点位移PR以下死点为位移零点，向上为位移的正方向，则任意时刻的悬点位移PR为：另外以上死点为位移零点，向下为位移的正方则悬点位移

为C.悬点速度

与加速度

为：悬点向上运动为位移的正方向，则

与

为：第25页，共54页，2023年，2月20日，星期一

图为悬点pr随theta变化x=acos((L.^2-1.2043)./(5.64*L));fi=x+beta;pr=3.935.*(0.372*pi-fi);第26页，共54页，2023年，2月20日，星期一（3）扭矩因数

图为扭矩因数tf随theta变化tf=-1.064*3.935.*sin(theta3-theta2)./(2.82*sin(theta3-theta4));第27页，共54页，2023年，2月20日，星期一3.3等效驱动力矩第28页，共54页，2023年，2月20日，星期一等效驱动力矩第29页，共54页，2023年，2月20日，星期一

图为等效驱动力矩Med随n的变化Med=151.3508.*Md;第30页，共54页，2023年，2月20日，星期一3.4悬点载荷分析1,2曲线是抽油机上冲程3,4曲线是抽油机下冲程第31页，共54页，2023年，2月20日，星期一图中—抽油杆在液柱中的重量36000N—抽油杆液面以上整个柱塞面积上的液柱重量30000N—悬点冲程长度m—抽油杆在液体载荷

作用下的静变形，0.2m图中四点左边如下第32页，共54页，2023年，2月20日，星期一由曲线形状得1,2,3,4段曲线的方程为：第33页，共54页，2023年，2月20日，星期一特s，

，

，

带入公式①得如下结果：其中以下死点为位移零点，向上为位移正方向，任意时间悬点的位移

为：已知如下参量：第34页，共54页，2023年，2月20日，星期一将

带入③得

第35页，共54页，2023年，2月20日，星期一

图为悬点载荷F随theta的分段变化F=(180000*pr+30000).*(theta>=0&theta<=0.1378*pi)+(66000).*(theta<=0.9631*pi&theta>=0.1378*pi)+(180000*pr-489025.17).*(theta>=0.9631*pi&theta<=1.2136*pi)+(30000).*(theta>=1.2136*pi&theta<=2*pi);第36页，共54页，2023年，2月20日，星期一

图为L随theta的变化曲线theta=0:0.005:2*piL=sqrt(19.3223-9.0759*cos(2.219*pi-theta))第37页，共54页，2023年，2月20日，星期一3.6等效力矩、等效阻力矩因为构件重量对结果影响不大，故为了简化方便计算可以忽略全部构件的重量:又因为①AB

段时第38页，共54页，2023年，2月20日，星期一②BC段时③CD段时④DA段时第39页，共54页，2023年，2月20日，星期一

图为等效阻力矩Mef与theta的分段关系Mef=((180000.*pr+30000)*3.935*1.064/2.82.*sin(theta3-theta2)./sin(theta3-theta4)).*(theta>=0&theta<=0.1378*pi)+(66000.*sin(theta3-theta2)./sin(theta3-theta4)*3.935).*(theta<=0.9631*pi&theta>=0.1378*pi)+((180000*pr-489025.17)*3.935.*sin(theta3-theta2)./sin(theta3-theta4)).*(theta>=0.9631*pi&theta<=1.2136*pi)+(30000*3.935.*sin(theta3-theta2)./sin(theta3-theta4)).*(theta>=1.2136*pi&theta<=2*pi);第40页，共54页，2023年，2月20日，星期一1）简化后模型曲柄摇杆机构简图如右图所示2）等效转动惯量Je将抽油机中所有运动件质量（从电动机转子到抽油机悬点）向曲柄处简化，用作用于曲柄轴上的转动惯量代替系统中的所有质量的转动惯量，这样简化要满足动能相等条件，即：第41页，共54页，2023年，2月20日，星期一:第42页，共54页，2023年，2月20日，星期一:第43页，共54页，2023年，2月20日，星期一

图为Je随theta的变化Je=107549.7594+102195.885.*(1.064/2.82.*sin(theta3-theta2)./sin(theta3-theta4)).^2;第44页，共54页，2023年，2月20日，星期一3.8系统运动微分方程非保守系统lagrange方程注意：角速度只是时间的函数。第45页，共54页，2023年，2月20日，星期一整理的等效构件运动微分方程为将:第46页，共54页，2023年，2月20日，星期一②ＢＣ段时:③ＣＤ段时：第47页，共54页，2023年，2月20日，星期一④ＤＡ段时:

图为

关于t的运动曲线第48页，共54页，2023年，2月20日，星期一

图为

关于

的运动曲线第49页，共54页，2023年，2月20日，星期一p=(6*2*pi/10+3.55)/2h=0.0001n=0p_rec(1)=p;x_rec(1)=0;x=0;whilex<0.448n=n+1;k1=h*f_x10_9(p,x);k2=h*f_x10_9(p+k1/2,x+h/2)k3=h*f_x10_9(p+k2/2,x+h/2)k4=h*f_x10_9(p+k3,x+h)p=p+(k1+2*k2+2*k3+k4)/6x=n*h;p_rec(n+1)=p;x_rec(n+1)=x;end;whilex<3.501n=n+1;k1=h*f_x10_10(p,x);k2=h*f_x10_10(p+k1/2,x+h/2)k3=h*f_x10_10(p+k2/2,x+h/2)k4=h*f_x10_10(p+k3,x+h)p=p+(k1+2*k2+2*k3+k4)/6;x=n*h;p_rec(n+1)=p;x_rec(n+1)=x;end;whilex<3.72n=n+1;k1=h*f_x10_11(p,x);k2=h*f_x10_11(p+k1/2,x+h/2)k3=h*f_x10_11(p+k2/2,x+h/2)k4=h*f_x10_11(p+k3,x+h)p=p+(k1+2*k2+2*k3+k4)/6;x=n*h;p_rec(n+1)=p;x_rec(n+1)=x;end第50页，共54页，2023年，2月20日，星期一whilex<2*pin=n+1;k1=h*f_x10_12(p,x);k2=h*f_x10_12(p+k1/2,x+h/2)k3=h*f_x10_12(p+k2/2,x+h/2)k4=h*f_x10_12(p+k3,x+h)p=p+(k1+2*k2+2*k3+k4)/6;x=n*