TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
题目有杆抽油系统的数学建模及诊断: [& c% h6 D P: c% W5 G+ g& v+ @* R
摘要:
?; v, P8 N0 o, H2 O有杆抽油系统目前广泛的的使用在开采原油活动中,它通过电机的旋转运动
- A ~6 n4 r: v* X, W转化为抽油杆上下往返周期运动,带动设置在杆下端泵的两个阀门相继开闭,从, s, d! h7 l5 I" T/ b U
而将地下上千米深处蕴藏的原油抽到地面上来。示功图通常用来描述抽油杆中任
m/ [& e1 |' P* G- A+ x& M意一水平截面处的基本信息,而悬点示功图可以初步诊断该井的工作状况,如产, Z% s' a+ R7 Q% _3 A- t6 }# _. r
量、气体影响、阀门漏液、沙堵等等。要精确诊断油井的工作状况,最好采用泵
; @0 f1 \4 G7 G3 u" U; R功图。然而泵在地下深处,使用仪器测试其示功数据实现困难大、成本高。通过) W' K. ~9 ?, K6 n
数学建模,把悬点示功图转化为杆上任意点的示功图(统称为地下示功图)并最
% o( o/ e/ f* b2 x( P& ~% q终确定泵功图,以准确诊断该井的工作状况具有很高的实际应用价值。
2 U6 [1 ~; Y+ I6 A+ A! ]通过查阅文献,我们在Gibbs 模型的基础上,运用Excel、SPSS 等软件进行
* ~. d! v5 S& i1 [4 f h数据分析,并通过MATLAB 软件编程求解,建立悬点运动规律、悬点示功图及泵5 K# H: F9 a+ ^- U2 G2 N" X
功图的数学模型。我们利用建立的模型对附件中所给出的两口油井的日产液量进: D3 _" Y/ s$ D5 \4 ?
行预测,同时对Gibbs 模型中粘滞阻尼系数的确定方法进行了研究与改进,大部
5 P2 {, _) g; s: v8 T2 V分模型都给出了具体的算法并附上程序,我们的工作主要包括:, A; t. c" P/ f7 `2 X) o* ?
(1)建立悬点运动规律模型,得到位移、速度和加速度表达式,并计算出0 z8 t) R# e2 _& s. X+ R& z
位移与附件中所给数据的平均误差e 为0.3806,拟合优度R 2 为0.7066,该位移- U# o G- o. k
模型可以一定程度的拟合悬点实际运动规律;' J" A' K) L* T- L! ^$ r( F
(2)对Gibbs 模型的边界条件和初始条件进行分析,利用附件1 和附件2
/ f: Z! A# ^0 D+ q& {; j的数据分别求出两口油井一级抽油杆和三级抽油杆的泵功图,由悬点示功图转换' ^. V) r* p- ]' i- V
23 L- ]! u* L2 |+ e/ y6 d) J
之后的泵功图有效地减弱了在上下冲程过程中抽油杆的波动,消除了摩擦和原油9 V8 W/ X& Y. q' T; k4 g/ [
稠度对荷载和位移的干扰,为依照示功图诊断油井工作状态提供了有效凭证;
+ z/ R- u. v( a7 a2 S9 b, @7 R(3)依照求得的泵功图,分别采用有效冲程法和面积法估算出两口油井的
& Q' M5 U* p6 d1 h. E0 l日产液量,有效冲程法的估算误差分别为32.25%和18.6%,面积法的估算误差
& _$ i0 z4 n/ @分别为9.21%和3.71%;第二小问,我们给出判断泵内是否充气的一种算法模型,4 m7 S4 a: R/ b
该模型将泵功图进行划分,通过计算划分后的区域面积并结合泵功图的曲率对泵! ^) H9 ~+ v4 c( z
内的气体影响进行判定;+ j8 ]. t( s1 @+ u+ |
(4)在一维Gibbs 振动模型的基础上使用有限元分析,加入抽油管、液柱# k( C/ G5 i. |5 E( s* |1 R
和油管振动这些因素,得到一种改进的Gibbs 模型;并使用迭代法从振动方程本& ^% F' R- C* E* P' w! F' \5 @% q
身推导出了粘滞阻尼系数的一种计算方法。
- ^! W" `' [2 [关键词: 泵功图有限元分析Gibbs 模型粘滞系数迭代法
G. ?% v- J6 G! G. O& x$ z0 E$ Q3 D5 _) G3 E- B8 W2 x+ ^; D1 o
# R8 [, U) d7 W) |! H) A
C19422003厉李王.rar
(714.1 KB, 下载次数: 27)
& Q y2 R3 b) \) h7 E6 A |
zan
|