地面示功图转化为泵示功图可用Gibbs提出的带黏滞阻尼的波动方程【见公式（A.1）】描述：

SY/T 74922021

槽钢标准.ou(x,t) at? ax2 at

2 =1 U(t) =. 2 r= ≤t≤T

立移变化的傅氏系数由实测D（t）和U(t）曲线数

c =2ZD(p)cos (n=0,1， ki t, =D(p)sin(3 2n元 (n=l,，n) 2,U(p)cos (n=0,1,..，n kp 2,U(p)sin| 2n元 (n=l,..,n)

U(x, t) = 2EA. 2

根据胡克定律，抽油杆柱任意深度x断面上的动载荷随时间的变化F（x，t）为【见公式（A.5） 公式（A.6)】：

式（A.5）和公式（A.6）中的函数和系数如下

SY/T 74922020

F(x, t) = 0o aP,(x) cosnot + ax sinnot

SY/T 74922021

L一一抽油杆柱长度，单位为米（m）。 以上公式只适用于单级抽油杆体，在实际计算时应附加接箍引起的阻尼。 在t时间，x断面上的总载荷等于动载荷F（x，t）加x断面以下的抽油杆柱在并液中的重力。 对于多级杆柱，根据杆柱上载荷的连续性原理，第二级杆柱顶部载荷和位移的傅氏系数与第一级 顶部（地面）的系数的关系为：

式中： L第一级杆柱长度，单位为米(m)。

002=00；0n2=EA.O,；Tn2=EA,P n2 EA

型包括抽油杆柱质量（Mass） 簧（Spring）一阻尼（Damping MSD连续模型

A.2.2.2MSD单元模型

将抽油杆柱离散化，分为m个单元（m趋近于无穷大且最长单元的长度趋近于0时，最后的解 答趋近于问题的精确解），然后对每个单元作如下的等效处理：去掉单元的弯曲刚度（弯曲刚度所引 起井身各处的杆管侧压力N,单独考虑），将单元的质量分配到单元下端结点，形成一些集中质量m， 可各个集中质量之间用拉压刚度为k：的弹簧代替，杆管侧压力N就作用在集中质量m；上，把杆液之 的黏性摩阻等效成阻尼系数为C，的阻尼器。由于弹簧没有质量，所以它两端的作用力相同，而集中 质量可认为是一个质点，所以它两端的位移可看作是相同的。 如图A.2所示，抽油杆柱的离散化方程为公式（A.8）

SY/T 74922020

SY/T 74922020

一位移（m）；t一时间（s）9一杆管间分布侧压力（N/m)，F一杆柱所受载荷（N）g一重力加速度（m/s"） 图A.3抽油杆柱MSD连续模型

图A.3抽油杆柱MSD连续模型

P 抽油杆柱材料的密度，单位为千克每立方米（kg/m）； E一一抽油杆柱材料的弹性模量，单位为帕斯卡（Pa）； 从 井液的动力黏度，单位为帕斯卡秒（Pa·s）； A一抽油杆柱的横截面面积，单位为平方米（m）： C 和杆管环形空间横截面形状尺寸有关的系数； 一 一井斜角，单位为弧度（rad）； 一杆管间库仑动摩擦系数： 杆管间分布侧压力，单位为牛每米（N/m）； 一抽油杆柱的单位长度杆重（质量），单位为千克每米（kg/m）。 抽油杆柱的运动是有阻尼的强迫振动，运动微分方程的解由两部分组成：自由振动和强迫振动， 前者由于阻尼的存在，会很快衰减消失，后者由于外部激扰的周期作用，会一直存在，为问题的稳 态解。 吸多 抽油杆柱的上部边界条件为【见公式（A.12)】：

SY/T 74922020

E(a) cosjo+be) sinjo) ..... (A.13)

A.3.1泵阀启闭判断法

SY/T 74922021

SY/T 74922020

则P点的曲率计算见公式（A.14）：

弧长△I的近似计算见公式（A.16）：

图A.10离散点曲率计算示意图

Al,= P.,P. + P.P+ PP. + P.P.?

则点P处的曲率变化量SK,根据P、Pi+I两点的曲率K、KI由公式（A.17）计算：

泵示功图曲线是通过数值方法计算得到的，曲线中含有复杂的频率成分。在实际计算中，为了 降低或消除其引起的曲率变化量的波动，一般采用相邻5点取平均值的方法计算中间点的曲率变化量 SK/，以提高计算精度，见公式（A.18）：

在泵示功图的高载荷段找出曲率变化最天的两个尖峰点，位移小者为固定阀开后点，反之为 闵关闭点；在低载荷段找出曲率变化最大的两个点，位移小者为游动阀关闭点，反之为游动阀开 ，分别计算固定阀启闭点的位置差和游动阀启闭点的位置差，则有效冲程为其中较小者。

A.3.2多边形逼近和矢量特征法

采用多边形逼近法对泵示功图进行预处理，再采用特征矢量法进行故障识别，充分考虑其他影响 因素判断泵有效冲程，然后根据判断的有效冲程计算产液量。具体步骤如下： a）泵示功图多边形逼近：由于计算得到的泵示功图数据点很多，而泵示功图的几何特征仅仅集 中在某些点上，其他点对泵示功图的几何特征影响不大，因此采用多边形逼近法对泵示功图 数据点进行预处理，过滤掉某些对泵示功图的几何特征影响不大的数据点 b）矢量特征法：将处理过的泵示功图标准化（无因次化）建筑常用表格，再用一系列连续的矢量把它描述出

SY/T 74922020

来，即建立该泵示功图的矢量链，把它和标准故障量链库中矢量链对比，来判别出一个或 多个故障，剔除抽油杆断脱、阀失灵、卡泵、气锁、脱筒等特殊工况井，对其余工作正常油 并结合生产参数，确定有效冲程，如图A.11所示。

A.4.1单张示功图计算日产液量

A.4.2日产液量平均值

Q实m一采集第m张示功图工况不变，持续T时间产液量，单位为立方米每天（m/d）； 一每采集两张示功图时间间隔，单位为分钟（min）

污水处理标准规范范本SY/T74922020附录B（资料性附录）油井基础数据油并基础数据见表B.1至表B.3。表B.1×××井基础数据表井型序号井深斜深井斜角方位角mm(°)(°)12注：井型分直井、定向井、斜直井、水平井等。表B.2×××井生产数据表杆柱组合（自上而下）层抽油机泵径泵挂冲程含水动液面油压套压当前修井位型号mmm杆径长度m材料%mMPaMPa完井日期mmm等级表B.3×××井流体物性数据表50℃原油黏度业标准信息服务平台原油密度mPa·skg/m312

中华人民共和国 石油天然气行业标准 抽油机井示功图法产液量计算技术规范 SY/T7492—2020 石油工业出版社出版 （北京安定门外安华里二区一号楼） 北京中石油彩色印刷有限责任公司排版印刷 新华书店北京发行所发行 880×1230毫米16开本1.25印张30千字印1—400 2021年1月北京第1版2021年1月北京第1次印刷 书号：155021·8193定价：25.00元 版权专有不得翻印