第6讲 计算机诊断模型有限差分法 、水平井诊断模型、现场应用.pptVIP

（一）等步长有限差分解 （一）等步长有限差分解 （一）等步长有限差分解 （二）变步长有限差分解 人工分析：该井从生产动态看,油井生产不正常,表现为动液面高,沉没度389 m ,从地面示功图看,深井泵受供液不足的影响,但动液面高,怀疑是凡尔堵塞。 软件分析：经采用故障诊断技术软件分析为气体影响、轻微结蜡。 实际情况：2000 年8 月21 日采用防气泵+ 螺旋式高效气锚防气工艺,下后动液面下降367 m ,证明诊断结果与实际完全吻合。 现场实例 3、凡尔漏失 实例3 :华77 - 71 井(见表3、图3) 。 现场实例 人工分析：该井从生产动态看,油井生产不正常,表现为动液面上升,沉没度710 m ,从地面示功图看,深井泵主要受固定凡尔漏失的影响。 软件分析：经采用故障诊断技术软件分析为固漏、严重结蜡。 实际情况：该井在1997 年9 月12日修井验证,固定凡尔失灵,油管结蜡严重,井口以下30 多根油管几乎堵死,蜡含砂、质硬,证明诊断结果与实际完全吻合。 现场实例 4、供液不足 实例4 :121 - 22 井(见表4、图4） 现场实例 有杆抽油系统的计算机 诊断模型 中国石油大学 吴晓东 诊断模型的有限差分解 傅氏级数解存在的问题是，由于傅氏系数数目太低就不能保证精度，因此计算时间较长。此外，傅氏级数法将使地面示功图平滑化。有限差分解可以克服这个问题。 可以用台劳级数推导出波动方程的有限差分解。设驴头下死点为x坐标原点，向下为正。u(x,t)也以向下为正，Δx为x的步长，Δt为时间步长。足标i表示位置，j表示时间，则 (a) (b) (c) 将（a）、(b)、(c)式代入一维波动方程，并经整理就得 ： 上式即为诊断模型的等步长有限差分解 已知地面光杆位移为u1,u2, ……uk，光杆动载荷为Fl，F2…… Fk,则边界条件为 (d) …………………… …………………… (d)和（e）即有限差分方程的边界条件 (e) 由上式有： （一）等步长有限差分解 差分格式，见左图。这种十字形差分格式结果会形成三角形，见右图。其中x为无法得出元点，但我们知道示功图是一个周期函数，根据这个概念得出 差分格式: （一）等步长有限差分解 上式实际是阻尼波动方程的两个初始条件。这样我们可以采取补格的办法求出全部未知点，见图 （一）等步长有限差分解 收敛条件：我们使用无阻尼波动方程推导收敛条件。无阻尼波动方程可写成 将(b)(c)代入（f）得出 (f) 上式是无阻尼波动方程的有限差分解 （一）等步长有限差分解 有限差分解的ui,j项的系数如果是负值的话，则其解是不稳定的，所以为了使其解是稳定的，必须满足以下条件： （一）等步长有限差分解 J.F.李1998耐庵提出了不同材料和杆径的诊断模型的变步长有限差分解如下. 对Δx与Δt的优选进行了研究的原则是在保证足够精度的前提下使计算时间尽量减少．通过 大量试验，得出以下结论： 单元长度Δx应为总长的15％至20％，但每级杆至少要有2至3个单元．对于变步长单元数，建议按下式计算： 式中 n1、n2：——分别为杆1与杆2单元数； l1、l2：——分别为杆1与杆2长度． （三）Δx与Δt的优选 时间步长Δt应桉下式进行选择，即： 显然在满足上式条件下， Δt愈小，精度愈高，计算时问愈长． 在满足上式的条件下，为了计算方便，最好使各杆杆Δt 一致，且采样点K取整数． （三）Δx与Δt的优选 水平井有杆泵诊断 模型 有杆泵诊断技术已有三十多年的历史。 60 年代的工作, 有S.G.Gibb s 提出的带阻尼的一维波动方程。在国内有张少南等的工作。 80 年代,D.R .Doty提出了一种改进模型, 是一组考虑了杆柱、液柱振动的双曲型偏微分方程组。在国内余国安等提出了同时考虑杆柱、液柱、管柱振动的数学模型。 大庆石油学院刘庆吉、宋考平于1997年给出了用于水平井有杆泵泵况诊断的, 同时考虑杆柱、液柱、管柱振动的数学模型。 该模型为六个方程组成的双曲型偏微分方程组,带有地面边界条件和周期条件。 1、运动和状态方程 取地面井口处为坐标原点, 铅直向下为y 轴正方向。设井轴曲线弧长为l 处, 切线与y 轴夹角为A(y ) , 其中y 为l 处在y 轴上投影点的坐标。经对井轴微元分析, 可得抽油杆的运动方程。 (1) 其中　v r = v r (y , t) 为抽油杆速度沿y 方向分量,m/s。f r = f r (y , t)