1、1有杆抽油井工程诊断方法及程序2第一章 绪论1.1 研究目的和意义有杆泵采油是世界石油工业传统的采油方式之一，也是迄今在采油工程中一直占主导地位的人工举升方式。在我国，采油生产井中大约有 90%采用有杆抽油技术，全国各油田产液量的 60%，产液量的 75%是靠有杆抽油技术采出的 1。然而由于抽油泵是在近千米到数千米的井下，工况十分复杂，工作环境极其恶劣，不但受“机、杆、泵”抽油设备的影响，而且直接受到“砂、蜡、气、水”的影响，故障发生率很高，1998 年的统计结果表明：我国平均每口有杆抽油井年作业 1.25 次，严重影响油田的正常生产。因此，及时、准确地掌握有杆抽油系统井下设备的工作状况，诊断
2、油井所存在的故障问题，制定合理的技术措施，使油井及时恢复正常生产，最大限度的提高原油产量、降低生产成本和提高抽油效率，对石油工业的发展和提高经济效益都具有非常重大的理论意义和现实意义 1。石油工业的发展对抽油丼故障诊断技术的研究提出强烈要求。几十年来，抽油丼故障诊断技术一直是国内外采油工程技术人员的一个重要研究课题，经过长期的研究与实践，抽油丼故障诊断技术取得了相当大的发展。1.2 国内外研究现状与发展趋势自从有杆泵在油田服役以来，有杆抽油井的故障诊断先后经历了从靠感觉分析到靠仪器测量分析，从仪器测量分析到计算机处理分析，并向智能化诊断不断迈进的发展过程。近几十年来诊断技术得到了飞速的发展，特
3、别是计算机诊断技术的出现，使有杆抽油系统故障诊断技术进入了一个崭新的阶段，即从定性分析到部分定量分析、从人工到初级智能。具体发展历程如下：1.2.1 国外发展概况早期的油井故障诊断仅靠手感，工作人员用手握住光杆，上下运动几个冲程，凭感觉来判断抽油泵的某些故障，这种方法只适用3于浅井，并且误差比较大 3。到了二十年代，1927 年发明地面光杆动力仪，利用光杆动力仪绘制光杆载荷与位移的关系曲线，即光杆示功图，然后对光杆示功图进行解释，以判断油井与设备故障，几十年来，许多国家进行了大量的研究工作：一方面，不断改进动力仪，提高检测精度；另一方面，不断改进示功图的解释方法，扩大解释范围 4。1936 年
4、美国的 Gilbert 和 Surgent 发明井下动力仪。这种方法是将井下动力仪随同抽油泵一起下入丼内，用其直接测量泵示功图，因这种方法耗资甚巨，工艺也较复杂，没有得到推广应用 5。1966 年美国壳牌石油公司的 Gibbs 建立了带阻尼的波动方程作为描述应力波在抽油杆柱中过程的基本微分方程，通过求解方程，可以得到抽油杆柱任意截面及泵处的示功图，随后很多人在模型完善及求解方面做了大量的工作 1。20 世纪 70 年代末期，Schafer 和 Jennings 在使用有限差分求解波动方程时，使用等步长差分形式，在抽油杆截面积和性质发生变化时，采用了等效值的概念，解决了多级杆及混合杆的波动方程求
5、解问题。随后 Nikea 采用有限元计算求解波动方程，取得了较好的效果 17。1969 年 8 月，美国有杆泵研究股份公司与中西部研究所经美国石油学会采油设备标准委员会批准，发表 API 标准示功图。1981 年美国塔尔萨大学的 Doty 和 Schmidt 建立了考虑液柱振动的二维预测数学模型。1988 年 Daek 等在走访许多著名专家后研制处有杆抽油井故障诊断专家系统，它是将地面实测的示功图转换成井下示功图，然后与标准示功图进行比较以判断故障类型。同年，Svinos 等，推出了一种由 Basic 语言编译的有杆泵诊断专家系统，该系统有 5 个模块，用产生式法则建立规则库，运用反向推理机建
6、立了一个可以识别典型示功图并计算出有关数据的专家4系统，然后利用这些数据诊断有杆抽油系统的故障。随后出现了各种诊断系统。美国 Arc 油气公司（DES）委内瑞拉Corpoven 公司和 Central de Venezuela（有杆泵抽油专家系统） 。1.2.2 国内发展概况1989 年西安石油学院的余国安和乌亦炯建立了综合考虑抽油杆柱、液柱和油管振动的三维数学模型，并进行了求解。这种三维振动的数学模型比起前两种数学模型来，在理论上显然更加符合实际一些 3。1990 年，石油大学用专家系统建造工具 M.1 开发有杆泵抽油丼故障诊断专家系统 ESROFD。1991 年，张润伟运用故障树分析技术也
7、建立了一个有杆泵抽油诊断专家系统 5。1993 年，天津大学与大港油田共同开发了抽油机井集成化智能诊断，它集成多种诊断方法进行综合诊断（但不是多种方法的有机融合，而是几种方法的简单堆砌） 。该系统中的元系统由 C 语言写成，子系统均用建造工具 M.1 专家系统外壳建立 5。上述研究成果有力地推进了有杆抽油系统故障诊断技术的进一步发展，但是，由于抽油泵工况诊断极其复杂，我国油田的抽油丼故障诊断还没有实现智能化，还处在人工诊断阶段，有待于进一步深入研究 5。近年来我国学者的研究与各石油单位联合，都具有一定的针对性，主要集中在：（1）有杆抽油系统预测和数值模拟的动力学模型建立；（2）预测和数值模拟数
8、学模型的计算求解方法；（3）井下工况诊断；（4）地面设备的仿真；（5）有杆抽油系统的优化设计 4。1.3 本文的主要工作本设计以有杆抽油丼抽油杆柱力学行为分析为基础，以有杆抽油机的诊断为研究对象，建立了基于波动方程的诊断数学模型，并5仔细推导了波动方程的数值解法。本文主要完成了以下几方面的工作：（1） 分析抽油机的运动特性，建立常规型抽油机的运动规律模型，得出了抽油机悬点位移、速度、加速度的运动规律曲线图；（2） 分析抽油杆柱的动力学特性，根据描述抽油杆动态特性的波动方程，有限差分求解波动方程；（3） 利用 Visual Basic 语言根据上述理论建立的模型编写程序，实现有杆抽油系统工况的计
9、算机诊断；（4） 根据所编程序对实际生产井进行实例计算，绘出井下泵功图，诊断油井工况。6第二章 游梁式抽油机动力学特性分析游梁式抽油机可以分为常规型、前置型、异相型等几种类型。抽油机作为主要动力传输单元是最重要的地面设备之一，其悬点运动规律是有杆抽油系统诊断数学模型的边界条件，其运行状况直接影响有杆抽油系统效率的高低，本文针对常规游梁抽油机运动特性的研究，建立了悬点运动规律数学模型，并通过模型使用 VB 程序来绘出悬点的位移、速度、加速度曲线图 2。2.1 常规游梁式抽油机简介抽油机是有杆抽油的地面驱动设备。按其基本结构，抽油机可分为游梁式和无游梁式，目前我国应用最为广泛的是游梁式抽油机。游梁
10、式抽油机主要由游梁连杆曲柄机构、减速机构、动力设备和辅助装置等四部分组成，游梁式抽油机工作时，电动机通过传动皮带将高速旋转运动传递给减速器的输出轴，经减速后由低速旋转的曲柄通过四连杆机构带动游梁作上下往复摆动。游梁前端圆弧状的驴头经悬绳器带动抽油杆作上下往复直线运动 。2图 2-1 常规型游梁式抽油机结构1刹车装置；2电动车；3减速器皮带轮；4减速器；5输入轴；6中间轴；7输出轴；8曲柄；9连杆轴；10支架；11曲柄平衡块；12连杆；13横船轴；147横船；15游梁平衡块；16游梁；17支架轴；18驴头；19悬绳器；20底座。常规型游梁式抽油机如图 2-1 所示。它是目前使用最广的一种抽油机。
11、其结构特点是：支架位于游梁的中部，驴头和曲柄连杆分别位于游梁的两端，曲柄轴中心基本位于游梁尾轴承的正下方，上下冲程运行时间相等。2.2 抽油机悬点运动分析掌握抽油机驴头悬点的位移、速度和加速度的变化规律是研究抽油装置动力学和进行抽油动态分析的基础。游梁式抽油机是以游梁支点和曲柄轴中心的连线作为固定杆，以曲柄、连杆和游梁后臂为三个运动杆所构成的四连杆机构（图 2-2） 。为了便于分析可简化为简谐运动和曲柄滑块机构。因为简谐运动模型只能用于粗略估算和简单分析，所以本文主要研究抽油机悬点运动的精确模型分析 2。图 2-2 常规型抽油机四连杆机构简图2.2.1 简化分析2.2.1.1 简谐运动视 r/
12、l0 及 r/b0 时，可将游梁与连杆的连接点 B 的运动简化为简谐运动，即认为 B 点的运动规律和 D 点做圆周运动时在垂直中心线上的投影（C 点）的运动规律相同，其悬点位移 SA、速度 VA8和加速度 aA分别为（2-1）)cos1(rbaSA（2-2）inv（2-3）cos2rbaA式中 曲柄转角 ；)(t曲柄角速度；t时间。2.2.1.2 曲柄滑块机构假设曲柄半径 r 与连杆长度 l 的比值范围为 0r/l1/4,抽油机的运动可简化为曲柄滑块机构的运动。令 =r/l，悬点的运动规律为 6(2-4)sin2co1(rbaSA(2-5)isiv(2-6)2cos(2rbaA由式（2-6）对
13、 求导并令其等于 0，可求得在 和 0（上、下死点）处悬点的最大加速度。180曲柄滑块机构模型是常用的模型，可用于一般的计算和分析。而简谐运动模型只能用于粗略估算和简单分析。2.2.2 精确分析常规型游梁式抽油机的几何关系于图 2-3 所示.采用以下符号表示抽油机的几何关系:9图 2-3 常规抽油机几何结构图中：A游梁前臂长度；C游梁后臂长度；P连线长度；R曲柄半径；I游梁轴中心到曲柄轴中心的水平距离；XC 与 J 的夹角; H游梁轴中心到底座底部的高度;K 与 J 的夹角；G曲柄轴中心到底座底部的高度;C 与 K 的夹角；K游梁轴中心到曲柄轴中心的距离； 驴头在下死点位置的 角;J曲柄-连杆
14、轴径b到游梁轴中心的距离; 驴头在上死点位置的 角;曲柄半径 R 与连杆 P 之t间的夹角；游梁后臂 C 与连杆 P 之间的夹角;R 与零度线的夹角即曲柄转角；零度线与 K 的夹角。运动分析：根据常规型抽油机几何结构参数，由常规行抽油机构成的四边形的角度关系得 ：2)(2arcos2RPKCcos2J)arcs(J22GHIK)sin(arcJRCKRPb2arcos210CKRPt 2)(arcos2根据前述常规型抽油机的几何参数间的关系式,可得:悬点位移方程：（2-7）)(bAS位移比：（2-8）tbPR悬点速度方程：（2-9）sinCAvA悬点加速度方程：（2-10）aA2式（2-9） 、 （2-10）中的 和 表示 对 的一阶和二阶导数，由复合函数直接求导得： 2)/()( JZJBJZ2/yURyJZB/yRU/JKJ/)sin(D/co2以上公式的中间变量 B、Z、U、y 及其导数分别为：
