1、采油工程课程设计姓名: 李健星 班级: 1班学号: 915463中国石油大学(北京)二O一二年四月目 录1、设计基础数据:12、具体设计及计算步骤2(1) 油井流入动态计算2(2) 流体物性参数计算方法4(3) 井筒温度场的计算6(4) 井筒多相流的计算7(5) 悬点载荷和抽油杆柱设计计算16(6) 抽油机校核21(7) 泵效计算21(8) 举升效率计算243、设计计算总结果26有杆抽油系统包括油层,井筒流体、油管、抽油杆、泵、抽油机、电动机、地面出油管线直到油气分离器。有杆抽油系统设计就是选择合理的机,杆,泵,管以及相应的抽汲参数,目的是挖掘油井潜力,使生产压力差合理,抽油设备工作安全、高效
2、及达到较好的经济效益。本次采油工程课程设计的主要内容是进行有杆抽油生产系统设计,通过设计计算,让学生了解有杆抽油生产系统的组成、设计原理及设计思路。1、设计基础数据: 井深:2000+学号末两位6310m=2630m 套管内径:0.124m 油层静压:给定地层压力系数为1.2MPa/100m,即油层静压为井深2630m/100m1.2MPa=31.56MPa油层温度:90 恒温层温度:16 地面脱气油粘度:30mPa.s油相对密度:0.84 气相对密度:0.76 水相对密度:1.0 油饱和压力:10MPa含水率:0.4 套压:0.5MPa油压:1 MPa 生产气油比:50m3/m3 原产液量(
3、测试点):30t/d 原井底流压(测试点):12MPa(根据测试液面计算得到)抽油机型号:CYJ10353HB 配产量:50t/d 泵径:44mm(如果产量低泵径可改为56mm,70mm) 冲程:3m 冲次:6rpm 沉没压力:3MPa 电机额定功率:37kw2、具体设计及计算步骤(1) 油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系,它反映了油藏向该井供油的能力,从单井来讲,IPR曲线表示了油层工作特性。因而,他既是确定油井合理工作方式的依据,也是分析油井动态的基础。本次设计油井流入动态计算采用Petrobras方法。Petrobras方法计算综合IPR曲线的实质是按含水率取纯
4、油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算采液指数时,是按产量加权平均;当预测产量或流压加权求平均值。采液指数计算已知一个测试点:、和饱和压力及油藏压力。 因为,= =2.083t/(d.MPa)某一产量下的流压Pwf=j()=2.083 x(26.410)=34.161t/d=+=34.161+=45.733t/d-油IPR曲线的最大产油量。 当0q时,令q=10 t/d,则p=21.599MPa同理,q=20 t/d,P=16.798MPa q=30 t/d,P=11.998 MPa 当qq时,令q=40 t/d,则按流压加权平均进行推导得:P=f+0.125(1-f)P-1
5、+=0.4+0.125(1-0.4)10-1+=6.910MPa 当qq时,令q=46t/dP=0.4=2.521MPa综上,井底流压与产量的关系列表如下:Pwf/ MPa21.59916.79811.9986.9102.521Q/(t/d)1020304046得到油井的流入动态曲线如下图:图1 油井IPR曲线Pwf(MPa)Q/(t/d)(2) 流体物性参数计算方法在地层的压力26.4Mpa和温度90时。 原油的API度 y= 36.95地面条件下的原油相对密度:0.84。 溶解油气比的计算因为15y=36.95,使用Lastater的相关式 (12)式中,地面脱气原油的有效分子量; 天然气
6、的摩尔分数。 的计算y0.13(如果0.13,则取=0.13); =50+36 ; =75+84 ()。由计算得到,由于该段的压力大于饱和压力的值,所以该段的流型为纯液流。 计算该段的压力梯度。由压力梯度的计算公式: =843.40;摩擦损失系数可由图二查得。可据此计算对应于的该段管长(深度差)。 将第步计算得的与第步估计的进行比较,两者之差超过允许范围,则以新的作为估算值,重复的计算,使计算的与估计的之差在允许范围内为止。该过程之中只迭代一次。2)由井口向下计算至泵出口处,计算泵排出口压力PZ。采用压力增量迭代方法,首先估算迭代压力。同样为了减小工作量,也采用只迭代一次的方法。计算井筒多相管
7、流时,首先计算井筒温度场、流体物性参数,然后利用Orkiszewski方法判断流型,进行压力梯度计算,最后计算出压力增量和泵排出口压力PZ。按压力增量迭代的步骤已知任一点(井底或井口)的压力,选取合适的深度间隔(可将管等分为n段)。估计一个对应于计算间隔的压力增量。计算该段的和 ,以及、下的流体性质参数。计算该段压力梯度计算对应于的压力增量比较压力增量的估计量与计算值 ,若二者之差不在允许范围内,则以计算值作为新的估计值,重复第步,使两者之差在允许范围之内为止。计算该段下端对应的深度和压力 以处的压力为起点压力重复第步,计算下一段的深度和压力 ,直到各段累加深度等于或大于管长时为止。1.5.3
8、计算气-液两相垂直管流的Orkiszewski方法本设计井筒多相流计算采用Orkiszewski方法。Orkiszewski法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图1所示。在处理过渡性流型时,采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。图1 气液混合物流动型态(Orkiszewski) 1.压力降公式及流动型态划分界限由前面垂直管流能量方程可知,其压力降是摩擦能量损失、势能变化和动能变化之和。由式(2-36)可直接写出多项垂直管流的压力降公式: (26)式中 压力,Pa; 摩擦损失梯度,Pa/m; 深度,m;
9、重力加速度,m/s2; 混合物密度,kg/m3; 混合物流速,m/s。动能项只是在雾流情况下才有明显的意义。出现雾流时,气体体积流量远大于液体体积流量。根据气体定律,动能变化可表示为: (27) 式中 管子流通截面积,m2; 流体总质量流量,kg/s; 气体体积流量,m3/s。将式(27)代入式(26),并取,经过整理后可得: (28) 式中 计算管段压力降,Pa; 计算管段的深度差,m; 计算管段的平均压力,Pa。不同流动型态下的和的计算方法不同,下面按流型分别介绍。(1)泡流平均密度 式中 气相存容比(含气率),计算管段中气相体积与管段容积之比值; 液相存容比(持液率),计算管段中液相体积
10、与管段容积之比值; 在下气、液和混合物的密度,kg/m3。气相存容比由滑脱速度来计算。滑脱速度定义为:气相流速与液相流速之差。 可解出: 式中 滑脱速度,由实验确定,m/s; 、气相和液相的表观流速,m/s。泡流摩擦损失梯度按液相进行计算: 式中 摩擦阻力系数; 液相真实流速,m/s。摩擦阻力系数可根据管壁相对粗造度和液相雷诺数查图2。液相雷诺数: 式中 在下的液体粘度,油、水混合物在未乳化的情况下可取其体积加权平均值,Pa.s。图 2 (2)段塞流混合物平均密度 (34)式中 液体分布系数; 滑脱速度,m/s。滑脱速度可用Griffith和Wallis提出的公式计算: (35) (3)过渡流
11、过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算,然后用内插方法来确定相应的数值。 (36) (37)式中的、及、为分别按段塞流和雾流计算的混合物密度及摩擦梯度。(4)雾流雾流混合物密度计算公式与泡流相同: 由于雾流的气液无相对运动速度,即滑脱速度接近于雾,基本上没有滑脱。所以 (38) 摩擦梯度则按连续的气相进行计算,即 (39) 式中 气体表观流速, ,m/s。雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度由图2查得。按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同,只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。图3为Orkiszewsk
12、i方法的计算流程框图。图3 Orkiszewski方法计算流程框图(5)悬点载荷和抽油杆柱设计计算 抽油杆柱设计的主要应用采油工程设计与原理。因为杆柱的最大、最小载荷与杆长不是线性关系,设计相对复杂。例如在考虑抽油杆弹性时的悬点载荷、在考虑杆柱摩擦时的悬点载荷公式与杆长不是线性关系。原因之二是因为杆、管环空中的压力分布取决于杆径,而杆柱的设计又用到杆、管环空中的压力分布。 本次课程设计仅将杆、管环空中的压力分布给定(按油水两相、不考虑摩擦时的压力分布),杆柱的最大、最小载荷公式采用与杆长成线性关系的公式。它是针对液体粘度较低、直井、游梁抽油机的杆柱载荷公式。 悬点最大、最小载荷的计算公式: (
13、40) (41) (42) 式中:第i级杆每米杆在空气中的质量,Kg/m 第i级杆杆长,m; i 抽油杆级数,从下向上计数; PZ泵排出口压力,Pa;PN泵的沉没压力,Pa;N冲次,rpm;S光杆冲程,m;fP活塞截面积,m2;g重力加速度,m/s2; (43) (44) 式中:令fr0=0 Pj第j级抽油杆底部断面处压力,Pa: (45) Pt井口油压,Pa; 0地面油密度,kg/m3; fw体积含水率,小数;应力范围比计算公式: (46) (47) 抽油杆柱的许用最大应力的计算公式:式中:抽油杆许用最大应力,Pa; T抽油杆最小抗张强度,对C级杆,T=6.3*108Pa,对D级杆T=8.1
14、*108Pa; 抽油杆最小应力,Pa; 使用系数,考虑到流体腐蚀性等因素而附加的系数(小于或等于1.0),使用时可考表2来选值。表2 抽油杆的使用系数使用介质API D级杆API C级杆无腐蚀性1.001.00矿化水0.900.65含硫化氢0.700.50若抽油杆的应力范围比小于则认为抽油杆满足强度要求,此时杆组长度可根据直接推导出杆柱长度的显示公式。对于液体粘度低的油井可不考虑采用加重杆,抽油杆自下而上依次增粗,所以应先给定最小杆径(19mm)然后自下而上依次设计。有应力范围比的计算公式即给定的应力范围比(0.85)计算第一级杆长L1,若L1大于等于泵深L,则抽油杆为单级杆,杆长为L,并计算
15、相应的应力范围比,若L1小于泵深L,则由应力范围比的计算公式及给定的应力范围比计算第二级杆长L2,若L2大于等于(L-L1),则第二级杆长为L2,并计算相应的应力范围比,若L2小于(L-L1),则同理进行设计。在设计中若杆径为25mm仍不能满足强度要求,则需改变抽汲参数。在设计中若杆径小于或等于25mm并满足强度要求,则杆柱设计结束。此为杆柱非等强度设计方法。若采用等强度设计方法,则需降低重新设计杆的长度。在设计抽油杆的过程中油管直径一般取(外径73mm,内径62mm)。若泵径大于或等于70mm,则油管全用(外径89mm, 内径76mm),原因是作业时大柱塞不能下入小直径油管中;若采用25mm
16、抽油杆,则相应油管直径应用,原因是25mm抽油杆节箍为55mm,与62mm油管间隙太小。当采用多级杆时油管长度比25mm杆长多10m。在本次课程设计中采用单级杆设计(19mm)。设计内容如下:由于采用单级杆设计,且杆径为19mm,所以选用油管的直径为:62mm。计算内容和步骤:1、 最大载荷: =0.0014999110=1499.9N;由于是单级的计算,所以简化为: =78509.81200=26174.24N (1499.9+26174.24)(1+) =29343.66N 2、 最小载荷: 式中:令fr0=0. 由于,在该设计过程之中,只有一级杆,所以公式变为: =1+10.631= 1
17、1.631 Mpa =26174.24N 11.631 (0.0014999-0) 10=8728.9N =8728.9 - =7149.68N (6) 抽油机校核1)最大扭矩计算公式 =1800 3+0.202 3 (29343.66 - 7149.68)= 18849.55N.m 2)电动机功率计算, =7860.53W 所以,得到电机的计算功率小于电机的额定功率,因而符合要求。(7) 泵效计算 (1)泵效及其影响因素在抽油井生产过程中,实际产量Q一般都比理论产量Qt要低,两者的比值叫泵效,表示, (50) (2)产量计算 根据影响泵效的三方面的因素,实际产量的计算公式为 (51)式中:Q
18、实际产量,m3/d; Qt理论产量,m3/d; Sp柱塞冲程,m; S光杆冲程,m; 抽油杆柱和油管柱弹性伸缩引起冲程损失系数; Bl泵内液体的体积系数; 泵的充满系数; qleak检泵初期的漏失量,m3/d;1) 理论排量计算 =1400 0.001499936=37.80 m3/d 2) 冲程损失系数的计算根据静载荷和惯性载荷对光杆冲程的影响计算当油管未锚定时; 由于只有一级抽油杆柱,所以公式简化为:=1.018式中:uL/a=0.1478 曲柄角速度,rad/s;N/30=6/30=0.6283; a声波在抽油杆柱中的传播速度,5100m/s; =1 10 0.0014999=1499.
19、9NPZ泵排出口压力,Pa;Pin泵内压力,Pa;当液体粘度较低时,可忽略泵吸入口压力,故PinPN;PN泵的沉没压力,Pa; fp、fr、ft活塞、抽油杆及油管金属截面积,m2; L抽油杆柱总长度,m; l液体密度,kg/m3; E钢的弹性模数,2.061011Pa; Lf动液面深度,m; L1、L2、L3每级抽油杆的长度,m; fr1、fr2 、fr3每级抽油杆的截面积,m23) 充满系数的计算 = 0.4814 式中:K泵内余隙比;取0.1. R泵内气液比; = =0.892 =50,m3/m3;=10m3/m3;=3M Pa;=0.4;P0=105Pa; T0=293K;273t=35
20、1.66;Z=0.96 4) 泵内液体的体积系数Bl =1.0462 5)漏失量的计算 检泵初期的漏失量为=0m3/d; D=0.044m;=0.00053Pas;l= 1.5m;PPZPN=10Pa;g=9.8m/s2;e= 0.00005m;=0.6m/s;所以最终算出泵的效率: =46.82(8)举升效率计算 光杆功率:P光= SN/60 = 8728.936/60=2618.67kw水力功率:P水力Q实际(PZPN)/86.4=17.7110/86400=204.94 井下效率:井下P水力 / P光 =0.0783 地面效率:地P光/ P电机 =0.5831 系统效率:总P地* P井下
21、 =0.4351 3、设计计算总结果基础数据设计结果班级 1班配产量50t/d班级内序号 63采油指数5.338t/(d.MPa)井深2630米井底流压12MPa静压31.56MPa下泵深度1200m油层温度90泵排出口压力3.6 MPa 含水率0.4悬点最大载荷29343.66N套压0.5MPa悬点最小载荷9835.26油压1MPa液柱载荷1499.9N生产气油比50m3/m3泵效0.4682抽油机型号CYJ10353HB最大扭矩18849.55N.m泵径44mm光杆功率2618.67kw冲程3m水力功率204.94冲次6rpm系统效率0.4351抽油杆是否满足生产满足抽油机是否满足生产满足
22、 4.课程设计总结在这次采油工程的课程设计过程之中:1、熟悉并掌握了。IPR曲线的画法,了解了IPR曲线的应用;2、多相管流的计算是一个难点,首先是计算公式很多,计算重复性大,其次是过程复杂,总来说,掌握了多相管流的Orkiszewski方法计算步骤,能运用压力迭代法和深度迭代法;3、掌握了悬点载荷及抽油管柱的设计计算泵效计算。这次实习的收获很多,通过实习对书本上的知识加深了理解,在实习的过程中也培养了自己独立思考和解决问题的能力。比如在实习过程中对一些数据的计算,我应用了电子表格的公式计算,降低了计算出错的概率。但在实习中我也感到一些不足,由于自己对采油工程设计做得少,很多地方不能得心应手,
23、以后只有加大练习才能够做的更好。二一二年十一月1. 80196单片机IP研究与实现,TN914.42 2. AT89S52单片机实验系统的开发与应用,TG155.1 F406 3. 基于单片机的LED三维动态信息显示系统,O536 TG174.444 4. 基于单片机的IGBT光伏充电控制器的研究,TV732.1 TV312 5. 基于89C52单片机的印刷品色彩质量检测系统的研究,TP391.41 6. 基于单片机+CPLD体系结构的信标机设计,TU858.3 TN915.627. 基于单片机SPCE061A的汽车空调控制系统,TM774 TM621.3 8. 带有IEEE488接口的通用单
24、片机系统方案设计与研究,TN015 9. 基于VC的单片机软件式开发平台,TG155.1 F406 10. 基于VB的单片机虚拟实验软件的研究与开发,TG155.1 F406 11. 采用单片机的电阻点焊智能控制器开发,TG155.1 F406 12. 基于51系列单片机的PROFIBUS-DP智能从站研究,TG155.1 F406 13. 八位单片机以太网接入研究与实现,TG155.1 F406 14. 基于单片机与Internet的数控机床远程监控系统的研发,R319 TP319 15. 基于单片机和DSP控制的医用输液泵的研究,U467.11 16. 基于单片机控制新型逆变稳压电源的设计
25、与仿真,F426.22 TP311.52 17. 基于8位单片机的摩托车发动机电控单元软硬件的开发,TB61 18. 基于430单片机的变压器监控终端的研究,TG155.1 F406 19. 逆变点焊单片机控制系统研究,TG131 TG113.14 20. 单片机控制数字变量柱塞泵的研究,F426.22 TP311.52 21. 基于单片机控制的高通量药物筛选及检测系统开发,R730.55 R734.2 22. MCS8051以及DS80C320单片机软核的设计,TP391 23. 基于AVR单片机的应用设计实践,TN015 24. LPC2210单片机的KGW脉冲固体激光掩膜加工控制系统研究
26、,TG131 TG113.14 25. 基于单片机控制的交流伺服系统的多梳栉经编机的研究,TN916 TP317 26. 80C196单片机在铁路客车发电机控制系统中的应用研究,TP368.1 TP393 27. 基于单片机的工程车辆3参数自动换档技术研究,F426.22 TP311.52 28. 削方制材机摇尺机构单片机控制装置的研制,TH213.6 29. 8XC196单片机集成开发环境的研制,F426.22 TP311.52 30. 基于单片机与PC的光电靶测试系统研究,O536 TG174.444 31. 手机和单片机控制系统的理论与应用研究,TG155.1 F406 32. 基于单片
27、机数控实验教学绘图仪研究,TN916 TP317 33. 基于单片机控制的脉冲电化学齿轮修形研究,R319 TP319 34. 基于AT89S52单片机的三相电度表研究,TP274.2 35. 基于MSP430单片机的嵌入式网络终端,TH812 TP368.1 36. 基于MSP430单片机地下车库通风控制系统设计,TP273.2 TG333.26 37. 基于PIC单片机的血压测量计设计,TP311.5 TM910.6 38. 基于单片机的标记打印机的研究与开发,TP277 TG156.82 39. 基于C8051F单片机和CANbus的航空三轴伺服转台控制系统的设计与研究,TP368.1
28、V217.2 40. 基于单片机的全位置自动焊接控制系统的研究 项目论证,项目可行性研究报告,可行性研究报告,项目推广,项目研究报告,项目设计,项目建议书,项目可研报告,本文档支持完整下载,支持任意编辑!选择我们,选择成功!项目论证,项目可行性研究报告,可行性研究报告,项目推广,项目研究报告,项目设计,项目建议书,项目可研报告,本文档支持完整下载,支持任意编辑!选择我们,选择成功!毕业论文,毕业设计,毕业论文设计,商业计划、商业策划、大学生商业计划书、大学生商业策划书、大学生创业计划书,毕业论文,毕业设计,毕业论文设计,商业计划、商业策划、大学生商业计划书、大学生商业策划书、大学生创业计划书毕
29、业论文,毕业设计,毕业论文设计,商业计划、商业策划、大学生商业计划书、大学生商业策划书、大学生创业计划书毕业论文,毕业设计,毕业论文设计,商业计划、商业策划、大学生商业计划书、大学生商业策划书、大学生创业计划书项目论证,项目可行性研究报告,可行性研究报告,项目推广,项目研究报告,项目设计,项目建议书,项目可研报告,本文档支持完整下载,支持任意编辑!选择我们,选择成功!项目论证,项目可行性研究报告,可行性研究报告,项目推广,项目研究报告,项目设计,项目建议书,项目可研报告,本文档支持完整下载,支持任意编辑!选择我们,选择成功!单片机论文,毕业设计,毕业论文,单片机设计,硕士论文,研究生论文,单片机研究论文,单片机设计论文,优秀毕业论文,毕业论文设计,毕业过关论文,毕业设计,毕业设计说明,毕业论文,单片机论文,基于单片机论文,毕业论文终稿,毕业论文初稿,本文档支持完整下载,支持任意编辑!本文档全网独一无二,放心使用,下载这篇文档,定会成功!27