采油工程课程设计.docx

1、采油工程课程设计采油工程课程设计课程设计姓 名：胡燕班 级： 1203学 号： 920956完成日期： 2014 年 4 月 29 日中国石油大学（北京）远程教育学院一、给定设计基础数据： 1二、设计计算步骤 12.1 油井流入动态计算 12.2 井筒多相流的计算 32.3 悬点载荷和抽油杆柱设计计算 102.4 抽油机校核 142.5 泵效计算 142.6 举升效率计算 17三、课程设计总结 17一、给定设计基础数据：井深： 2000+5610=2560m套管内径： 0.124m油层静压： 2560/1001.2 =30.72Pa油层温度： 70恒温层温度： 16地面脱气油粘度： 30mPa

2、.s油相对密度： 0.84气相对密度： 0.76水相对密度： 1.0油饱和压力： 10MPa含水率： 0.4套压： 0.5MPa油压： 1 MPa生产气油比： 50m3/m3测试产液量： 30t/d抽油机型号： CYJ10353HB电机额定功率： 37kw配产量： 50t/d管式泵径： 56mm冲程： 3m冲次： 6rpm沉没压力： 3MPa二、设计计算步骤2.1 油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系，它反映了油藏向该井供油的能力。从单井来讲， IPR 曲线表示了油层工作特性。因而，它既是确定油井合理工作方式的依据，也是分析油井动态的基础。本次设计油井流入动态计算采用

3、Petro bras方法 Petro bras 方法计算综合 IPR 曲线的实质是按含水率取纯油 IPR 曲线和水 IPR 曲线1的加权平均值。当已知测试点计算采液指数时，是按产量加权平均 ;预测产量时 ,按流压加权平均。(1)采液指数计算已知一个测试点： p wftest 、 qtxst 和饱和压力 Pb 及油藏压力 P 。如果 pwftestpb 则 j1qnest=30/（ 30.72-12） =1.6/(d.Mpa)p1pwftest如果 pwftestpb 则 jqtxest=1.11/(d.Mpa)pb A) f w ( p1 pwfxets )(1f w )( p1pbL8(2)

4、某一产量 qt 下的流压 pwfqt j ( p t pb ) =1.6*(30.72-10)=33.15t/dqomzx qb jp b =33.15+1.6*10/1.8=42.04t/d1.8qomzx -油 IPR 曲线的最大产油量。若 0q1qt 则，令 q t1=10 t/d，则 pwf p1qt =16.47 Mpaj若 qt q1qomzx令 q t 4=50 t/d，则按流压加权平均进行推导得；pwff w ( p1q1 )0.125(1 f w ) pb 181 80( q1qb ) jqo maxqb=6.34Mpa若 qomzx q1 ，则综合 IPR 曲线的斜率可近似

5、常数。pwff w ( p1qomzx )(q1qomzx )(8 f w 9)jj令 q t 6 =71t/d， P wf 6 =2.233 Mpa综上，井底流压与产量的关系列表如下：Pwf/Mpa16.4713.87312.010.06.342.233Q/(t/d)10203040.65350712得到油井的流入动态曲线如下图：IPR 曲线20) 15apM( 10fwP500 10 20 30 40 50 60 70 80 90Q（t/d ）图 1 油井 IPR曲线2.2 井筒多相流的计算井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增

6、加的动能和流体沿管路的摩阻损失，其数学表达式如下：dpdvmf mvm2dh m gsin+ m v m dh2 m /d*式中 m 为多相混合物的密度 ;v m 为多相混合物的流速 ;f m 为多相混合物流动时的摩擦阻力系数 ;d 为管径 ;p 为压力 ;h 为深度 ;g 为重力加速度 ; 为井斜角的余角。井筒多相管流计算包括两部分：（1）由井底向上计算至泵入口处；（2）油管内由井口向下计算至泵出口处。1）由井底向上计算至泵入口处，计算下泵深度 Lp 。采用深度增量迭代方法，首先估算迭代深度。在本设计中为了减小工作量，采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、 流体物性参

7、数， 然后利用 Orkiszewski 方法判断流型，进行压力梯度计算，最后计算出深度增量和下泵深度 Lp 。按深度增量迭代的步骤 :井底流压 12Mpa，假设压力降为 0.2 Mpa；估计一个对应的深度增量 h =40m，即深3度为 1960m 。由井温关系式可以计算得到该处的井温为： 89.96 。平均的压力和温度：T =（ 90+89.96） /2=89.98 。平均压力 P =11.9 Mpa 。由平均压力和平均温度计算的得到流体的物性参数为：溶解油气比RS =71.31 ；原油体积系数 B0 =1.25原油密度 P0 =739.00； 油水混合液的密度Pz =843.40 ； 死油粘

8、444度 od =6.537*10O水的粘度 w =3.263*10；活油粘度 =3.318*10 ；4液体的粘度 = 3.296*10；天然气的压缩因子 Z=0.9567； 天然气的密度g90.70 。以上单位均是标准单位。由以上的流体物性参数判断流型：不同流动型态下的 m 和 f 的计算方法不同，为此，计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表 1 所示。表 1 流型界限流动型态界限泡流qgLqtB段 塞 流qgLB ,vgLSqt过 渡 流LMvg LS雾流vgLM其中 LB =1.071-0.7277 n t2 / D 且 LB 0.13 （如果 LB 0.13 ，则

9、取 LB =0.13 ）；qgLS=50+36vgqt；qg0.75LMvg qt=75+84（） 。由计算得到，由于该段的压力大于饱和压力的值，所以该段的流型为纯液流。dPdP f dh g mdh mvmdvm计算该段的压力梯度dh 。由压力梯度的计算公式：4m=843.40 ； f =计算对应于 P 的该段管长 ( 深度差 ) h计 。将第 步计算得的 h计 与第步估计的 h 进行比较，两者之差超过允许范围，则以新的 h 作为估算值， 重复的计算， 使计算的与估计的 h 之差在允许范围 内为止。该过程之中只迭代一次。2）由井口向下计算至泵出口处，计算泵排出口压力 PZ。采用压力增量迭代方

10、法，首先估算迭代压力。同样为了减小工作量，也采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用 Orkiszewski 方法判断流型，进行压力梯度计算，最后计算出压力增量和泵排出口压力 PZ。按压力增量迭代的步骤已知任一点 (井底或井口 )的压力 P0 ， 选取合适的深度间隔 h (可将管 L 等分为 n段)。估计一个对应于计算间隔 h 的压力增量 P 。计算该段的 T 和 P ，以及 P 、 T 下的流体性质参数。dP计算该段压力梯度dhoPih dP计算对应于 h 的压力增量dho比较压力增量的估计量P 与计算值Pi，若二者之差不在允许范围内，则以计算值作

11、为新的估计值，重复第步，使两者之差在允许范围 o 之内为止。计算该段下端对应的深度Li 和压力 PiiLi ihPi PoPi1以 Li 处的压力 Pi为起点压力重复第步，计算下一段的深度 Li 1 和压力 Pi 1，直到各段累加深度等于或大于管长L 时为止。3）计算气 - 液两相垂直管流的Orkiszewski 方法本设计井筒多相流计算采用 Orkiszewski 方法。Orkiszewski 法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图1 所5示。在处理过渡性流型时，采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。

12、图 1 气液混合物流动型态 (Orkiszewski) 1.压力降公式及流动型态划分界限由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。由式 (2-36)可直接写出多项垂直管流的压力降公式：dPf dhg mdh mvm dvm(26)式中P 压力， Pa；f 摩擦损失梯度， Pa/m；h 深度， m；g 重力加速度， m/s2；m混合物密度， kg/m3；vm 混合物流速， m/s。动能项只是在雾流情况下才有明显的意义。出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：W qm vm dvmtg2dpAp P(27)式中Ap 管子流通截面

13、积， m2；Wt 流体总质量流量， kg/s；qg 气体体积流量， m3/s。6将式(27)代入式k，dPpk，mm， P P 经过整理后可得：(26)，并取 dhhm gf hk1Wt qg2(28)PkAp P式中Pk 计算管段压力降， Pa；hk 计算管段的深度差， m；P计算管段的平均压力， Pa。不同流动型态下的 m 和 f 的计算方法不同 ,下面按流型分别介绍。(1)泡流平均密度m H L L H g g 1 H g L H g gH L H g 1式中 HgHL气相存容比 (含气率 )，计算管段中气相体积与管段容积之比值；液相存容比 (持液率 )，计算管段中液相体积与管段容积之比

14、值；g、 、 m 在 P 、 T 下气、液和混合物的密度， kg/m3。气相存容比由滑脱速度 Vs 来计算。滑脱速度定义为：气相流速与液相流速之差。vsgvsLqgqqgtvs H g1 H gAp H gAp (1 H g )11t(1qt)24qgq可解出 Hg ： H g2vs Apvs Apvs Ap式中 vs 滑脱速度，由实验确定， m/s；vsg 、 vsL 气相和液相的表观流速， m/s。泡流摩擦损失梯度按液相进行计算：2qLL vLHvLHtf D 2Ap (1 H g )式中 f 摩擦阻力系数；vLH 液相真实流速， m/s。7摩擦阻力系数 f 可根据管壁相对粗造度 / D

15、和液相雷诺数 N Re 查图 2。液相雷诺数 :NReDvsL LL式中 L 在 P 、 T 下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值， Pa.s。图 2(2)段塞流混合物平均密度WtL vs ApmL(34)qtvs Ap式中 液体分布系数；vs 滑脱速度， m/s。8滑脱速度可用 Griffith 和 Wallis 提出的公式计算：vs C1C2 gD(35)(3)过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。LMvgvgLsmLMLSSLLMMiLSLMvgvgLgtLMLSSLLMMiLS(36)(37)式中

16、的 SL 、 SL 及 Mi 、 Mi 为分别按段塞流和雾流计算的混合物密度及摩擦梯度。(4)雾流雾流混合物密度计算公式与泡流相同：m H L L H g g (1 H g ) L H g g由于雾流的气液无相对运动速度，即滑脱速度接近于雾，基本上没有滑脱。所以qH ggqL qg(38)摩擦梯度则按连续的气相进行计算，即2ffgvsg2D(39)式中vsg 气体表观流速， vsgq g / A p，m/s。雾流摩擦系数可根据气体雷诺数N Reg 和液膜相对粗糙度由图2 查得。按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同，只是在计算混合物密度及摩擦之前需要根据流动型态界

17、限确定其流动型态。图3 为Orkiszewski 方法的计算流程框图。9以井口油压或井底流压为起点，选择合适的压力间隔 P ，假设 h计算平均 P 和 T，并求得在此 P 和 T 下的流体性质参数和流动参数，以及相应的流动型态界限 LB、 Lg 和 LM确定流动型态雾流过渡流段塞流泡流计算气相存容分别按段塞流和雾流计算滑脱速度、液计算气相存容比、平均密度计算平均密度及摩擦体分布系数、平均比、平均密度及摩擦梯度梯度，并进行内插密度和摩擦梯度和摩擦梯度计算并比较 h ，重复上述计算使 h 的计算值与假设值相等或在允许的误差范围内h重复上述步骤，直到 的等于或大于油层深度为止图 3 Orkiszew

18、ski 方法计算流程框图2.3 悬点载荷和抽油杆柱设计计算抽油杆柱设计的一般方法见采油工程设计与原理 。之所以设计方法较复杂，原因之一是因为杆柱的最大、最小载荷与杆长不是线性关系。例如在考虑抽油杆弹性时的悬点载荷、在考虑杆柱摩擦时的悬点载荷公式与杆长不是线性关系。原因之二是因为杆、管环空中的压力分布取决于杆径，而杆柱的设计有用到杆、管环空中的压力分布。由于综合课程设计时间较少，所以这里提供一种简化杆柱设计方法。暂将杆、管环空中的压力分布给定（按油水两相、不考虑摩擦时的压力分布） ，杆柱的最大、最小载荷公式采用与杆长成线性关系的下面公式。它是针对液体粘度较低、直井、游梁抽油机的杆柱载荷公式。悬点

19、最大、最小载荷的计算公式：iSN2Pmax(Wrj WL )(1)j11790（ 40）10i iWrj qrj Lrj gj 1 j 1WL f p (PZ PN )式中： qr i 第 i 级杆每米杆在空气中的质量， Kg/mLri 第 i 级杆杆长， m；i 抽油杆级数，从下向上计数；PZ泵排出口压力， Pa；PN泵的沉没压力， Pa；N冲次， rpm;S光杆冲程， m；fP活塞截面积， m2；g重力加速度， m/s2；iSN2iPminWrjWrjj 11790 j 1iiiWrjWrjPj( f rjf r1 j 1 )j 1j 1j 1式中：令 fr0=0Pj第 j 级抽油杆底部断

20、面处压力 ,Pa：j 1Pj Pt 0 (1 f w ) w f w g(L L t )t 1Pt井口套压， Pa;0地面油密度， kg/m3;fw 体积含水率，小数；L- 抽油杆总长度， mfr 抽油杆横截面积， m2；应力范围比 pL 计算公式：PLmax minall min11P max P minmax minfr frall ( T 0.5625 min) SF抽油杆柱的许用最大应力的计算公式： 4式中： all 抽油杆许用最大应力， Pa；T抽油杆最小抗张强度，对 C 级杆， T=6.3*108Pa,对 D 级杆 T=8.1*108Pa;m i n抽油杆最小应力， Pa；SF 使

21、用系数，考虑到流体腐蚀性等因素而附加的系数（小于或等于 1.0），使用时可考表 2 来选值。表 2 抽油杆的使用系数使用介质API D 级杆API C 级杆无腐蚀性1.001.00矿化水0.900.65含硫化氢0.700.50若抽油杆的应力范围比小于 pL 则认为抽油杆满足强度要求，此时杆组长度可根据 pL 直接推导出杆柱长度的显示公式。对于液体粘度低的油井可不考虑采用加重杆，抽油杆自下而上依次增粗，所以应先给定最小杆径（ 19mm）然后自下而上依次设计。有应力范围比的计算公式即给定的应力范围比（ pL 0.85）计算第一级杆长 L1，若 L1 大于等于泵深 L ，则抽油杆为单级杆，杆长为 L

22、 ，并计算相应的应力范围比，若 L1 小于泵深 L,则由应力范围比的计算公式及给定的应力范围比计算第二级杆长 L2，若 L2 大于等于（ L-L1 ） ,则第二级杆长为 L2 ，并计算相应的应力范围比，若 L2 小于（ L-L1 ）,则同理进行设计。在设计中若杆径为 25mm 仍不能满足强度要求，则需改变抽汲参数。在设计中若杆径小于或等于25mm 并满足强度要求， 则杆柱设计结束。 此为杆柱非等强度设计方法。 若采用等强度设计方法，则需降低 pL 重新设计杆的长度。在设计抽油杆的过程中油管直径一般取 2 12 （外径 73mm，内径 62mm）。若泵径12大于或等于 70mm，则油管全用 3

23、（外径 89mm, 内径 76mm），原因是作业时大柱塞不能下如小直径油管中；若采用 25mm 抽油杆，则相应油管直径应用 3 ，原因是 25mm抽油杆节箍为 55mm，与 62mm 油管间隙太小。当采用多级杆时 3 油管长度比 25mm 杆长多 10m。为了减小计算工作量，在本次课程设计中杆柱设计简化处理，采用单级杆设计（19mm）。设计内容如下：由于采用单级杆设计，且杆径为 19mm，所以选用油管的直径为： 62mm。计算内容和步骤：1、最大载荷：W Lf p( PZPN)6=0.0014999110=1499.9NiiWrjqrj Lrj gj 1j1；由于是单级的计算，所以简化为：WrqrLr g19210 64 =26174.24N=78509.81200iSN2Pmax(WrjWL )(1)j11790