TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
题目有杆抽油系统的数学建模及诊断7 m( `2 e" Q3 E( `& M
摘要:
- } F# |9 V) c8 H. P( P有杆抽油系统目前广泛的的使用在开采原油活动中,它通过电机的旋转运动
& y3 _3 V7 \8 W. Z. Q转化为抽油杆上下往返周期运动,带动设置在杆下端泵的两个阀门相继开闭,从. C* F# f2 j7 L5 m! o
而将地下上千米深处蕴藏的原油抽到地面上来。示功图通常用来描述抽油杆中任. p* f$ f* \; _- O# A
意一水平截面处的基本信息,而悬点示功图可以初步诊断该井的工作状况,如产2 ?2 l9 m$ K B5 R2 u" V
量、气体影响、阀门漏液、沙堵等等。要精确诊断油井的工作状况,最好采用泵$ S+ f* U& c, Y2 I# X* \
功图。然而泵在地下深处,使用仪器测试其示功数据实现困难大、成本高。通过 r& l1 q& B; |2 }' J
数学建模,把悬点示功图转化为杆上任意点的示功图(统称为地下示功图)并最
' O P( V1 ^9 R+ H6 }+ F B终确定泵功图,以准确诊断该井的工作状况具有很高的实际应用价值。
0 k2 a/ Y, g" F1 `通过查阅文献,我们在Gibbs 模型的基础上,运用Excel、SPSS 等软件进行! K% A j" z2 l2 u
数据分析,并通过MATLAB 软件编程求解,建立悬点运动规律、悬点示功图及泵
7 g5 v; t( V4 O- s0 f! L/ p功图的数学模型。我们利用建立的模型对附件中所给出的两口油井的日产液量进
( ~" _2 b& P4 S9 R Z3 K行预测,同时对Gibbs 模型中粘滞阻尼系数的确定方法进行了研究与改进,大部
: T D- a( @! h" O6 n0 c分模型都给出了具体的算法并附上程序,我们的工作主要包括:
9 \0 g7 K! X: O" f0 \(1)建立悬点运动规律模型,得到位移、速度和加速度表达式,并计算出
! O, V; k2 t" H$ F; M9 {# R" @位移与附件中所给数据的平均误差e 为0.3806,拟合优度R 2 为0.7066,该位移
7 D: }: b4 f7 P8 J) Y2 {6 x模型可以一定程度的拟合悬点实际运动规律;
5 H1 S I4 E* Z: z, \$ b(2)对Gibbs 模型的边界条件和初始条件进行分析,利用附件1 和附件2
, p" Y: M, [8 \) w$ S$ y" R7 r+ O的数据分别求出两口油井一级抽油杆和三级抽油杆的泵功图,由悬点示功图转换( X- G1 V$ q2 q9 I* D
2
- ~: D2 l: R' _9 H6 ]之后的泵功图有效地减弱了在上下冲程过程中抽油杆的波动,消除了摩擦和原油+ D$ F( L8 Z! q- g; W" G$ c
稠度对荷载和位移的干扰,为依照示功图诊断油井工作状态提供了有效凭证;( ^: F8 v& Y, q4 Y Z4 ?# y' O
(3)依照求得的泵功图,分别采用有效冲程法和面积法估算出两口油井的9 [3 Z; W E x1 a* G5 r
日产液量,有效冲程法的估算误差分别为32.25%和18.6%,面积法的估算误差" P8 z: ~. g. z$ W: f
分别为9.21%和3.71%;第二小问,我们给出判断泵内是否充气的一种算法模型,
+ [6 @3 G, U- i' N" p4 c9 B1 r该模型将泵功图进行划分,通过计算划分后的区域面积并结合泵功图的曲率对泵2 V" Q; i2 K" A1 I& z7 f
内的气体影响进行判定;# c: U5 q& m7 ]$ R5 l0 ]8 o/ }
(4)在一维Gibbs 振动模型的基础上使用有限元分析,加入抽油管、液柱; ~& M8 B" C/ X
和油管振动这些因素,得到一种改进的Gibbs 模型;并使用迭代法从振动方程本4 j% f, \: M1 I. Y4 ~
身推导出了粘滞阻尼系数的一种计算方法。
0 n7 }8 i8 V5 B% k0 l) Z/ M! d! i关键词: 泵功图有限元分析Gibbs 模型粘滞系数迭代法1 k3 j- D8 n) ?- M
3 M: B( u) f+ S/ h$ ?9 e6 E% c" d9 _7 x" |0 s
C19422003厉李王.rar
(714.1 KB, 下载次数: 27)
2 L# u' h/ Q- t, P a7 F( R- {5 K |
zan
|