采油工程课程设计报告

1、页眉采油工程课程设计指导书页脚页眉中国石油大学（北京）石油天然气工程学院页脚页眉本次采油工程课程设计的主要内容是进行有杆抽油生产系统设计，通过设计计算，让学生了解有杆抽油生产系统的组成、设计原理及设计思路。1有杆泵抽油生产系统设计11 有杆抽油生产系统设计原理有杆抽油系统包括油层，井筒流体、泵、油管、抽油杆、抽油机、电动机、地面出油管线直到油气分离器。有杆抽油系统设计就是选择合理的机，杆，泵，管以及相应的抽汲参数，目的是挖掘油井潜力，使生产压差合理，抽油设备工作安全、高效及达到较好的经济效益。在生产过程中，井口回压ph 基本保持不变，可取为常数。它与出油管线的长度、分离器的入口压力有关，此处取

2、ph1.0MPa 。抽油井井底流压为pwf 向上为多相管流 ,至泵下压力降至泵的沉没压力(或吸入口压力 ) pn ,抽油泵为增压设备 ,故泵出口压力增至 pz , 称为泵的排出口压力 .在向上 ,为抽油杆油管间的环空流动 .至井口 ,压力降至井口回压 ph 。(1) 设计内容对刚转为有杆泵抽油的井和少量需调整抽油机机型的有杆抽油井可初选抽油机机型。对大部分有杆抽油油井。 抽油机不变， 为己知。对于某一抽油机型号， 设计内容有：泵径、冲程、冲次、泵深及相应的泵径、杆长，并求载荷、应力、扭矩、功率、产量等技术指标。(2) 需要数据井：井深，套管直径，油层静压，油层温度混合物：油、气、水比重，饱和压

3、力页脚页眉生产数据：含水率，套压，油压，生产气油比，原产量，原流压（或原动液面）。(3) 设计方法这里介绍给定配产时有杆抽油系统的设计方法。首先需要获得油层的IPR 曲线。若没有井底流压的测试值，可根据测试液面和套压计算得井底流压，从而计算出采液指数及IPR 曲线。1） 根据测试液面计算测试点流压从井口到井底可分为三段。从井口到动液面为气柱段，若忽略气柱压力，则动液面顶端压力仍为套压。从动液面到吸入口为纯油柱段，可以将这一段分为许多小段，采用迭代压力方法可求出每小段油的密度，最后求出吸口处的压力。从吸入至油层中部分多相管流段。通过分小段计算多相管流压力分布，可求得测试点流压。2） 根据测试点流

4、压和产量计算IPR 曲线3） 给定配产量时有杆泵油井设计步骤（简化设计方法）a 利用 IPR 曲线，由给定产量 Q 计算流压。b 按 Q 由流压向上进行多相管流计算，得不同深度处的压力分布。一般分若干小段进行压力分布计算。为了计算简便，此处可按深度增量迭代方法分两段计算。若井底流压 pwf 高于饱和压力 pb ，则以饱和压力点 pb 为分界线分为两段， pwf 从到 pb 为一段，从 pb 到零为一段。若井底流压 pwf 低于饱和压力 pb ，则以 pwf / 2 为分界线分为两段，从pwf 到 pwf / 2 为一段，从 pwf / 2 到零为一段。c 根据泵沉没压力内插确定泵深；d 初选杆

5、、管直径，按Q 由井口向下进行杆、管环空压力分布计算，得不同深度处的压力分布，为了简化计算，给定压力分布；e 对某一抽汲参数组合：泵径、冲程、冲次、泵沉没压力，计算液柱载荷，设计抽油杆柱；页脚页眉f 计算扭矩和需要电机功率等校核抽油机：g 计算泵效：从而计算出产量Q h 判断 Q Q 。若不成立，则换另一组抽汲参数，转第 e 步；若成立转第 iQ步。i 计算举升效率。j 通过计算多组抽汲参数的产量，最后得到产量比配产高但最接近且经济、技术指标较好的抽汲参数组合。1.2 油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系，它反映了油藏向该井供油的能力，从单井来讲， IPR 曲线表示了油

6、层工作特性。因而，它既是确定油井合理工作方式的依据，也是分析油井动态的基础。本次设计油井流入动态计算采用Petrobras方法。Petrobras方法计算综合 IPR 曲线的实质是按含水率取纯油IPR 曲线和水 IPR 曲线的加权平均值。当已知测试点计算采液指数时，按产量加权平均。采液指数计算已知一个测试点；pwftest 、 qtxst 和饱和压力 pb 及油藏压力p 。如果 pwftestpb则 Jqnest（1）p1pwftest如果 pwftestpb采液指数Jqtest（ 2）pb A1 f w p1 pbf w p1pwftest1.8页脚页眉式中，A10.2( pwftest )

7、 0.8( pwftets ) 2pbpbqtest 对应流压 pwfxets 时总产液量；f w 含水率，小数：qomzx 油 IPR 曲线的最大产油量。某一产量qt 下的流压 pwf计算qtJ ( ptpb )qomzxqbJpb1.8 若 0q1qt 则pwfp1qtJ 若 qtq1qomzx 则按流压加权平均进行推导得；pwffw ( p1q1 )0.125(1fw ) pb 181 80(q1qb )Jqo maxqb 若 qomzxq1 ，则综合 IPR 曲线的斜率可近似常数。pwffw ( p1qomzx )(q1qomzx )(8 fw9)JJ13 流体物性参数计算方法(1)

8、原油密度计算1000( o1.20610 3 Rsg )oBo式中，o 在压力 P 及温度 t 下的原油密度， kgm3 ；o 地面条件下的原油相对密度；g 地面条件下的气体相对密度；Rs 在压力 p 及温度下的溶解油汽比，m3 m3 ；（3）（4）（5）（6）（7）（ 8）页脚页眉Bo 在压力 p 及温度 T 下的原油体积系数。(2)原油的 API 度141.5131.5yAPIo式中，y API 原油的 API 度。(3) 原油体积系数的计算Bo0.9720.000147F 1.175式中， F 5.615Rsg2.25t40o(4) 溶解油气比的计算当 yAPI15 时，使用 stand

9、ing的相关式RS0.17812 g(8.0558p 10A )1.2048式中， A 0.0125yAPI 0.00091(1.8t 32)t 温度，；p 泡点压力（在多相管流中取计算段的平均压力p ），Pa 。（ 9）（10）（ 11）当 yAPI 15 时，使用 Lastater的相关式RS23650oyng（12）mo 1yng式中， mo 地面脱气原油的有效分子量；yng 天然气的摩尔分数。其中，mo 和 yng 可以通过差图来获得。为便于计算，可以采用以下公式计算mo 和 yng 。mo 的计算当API38.3 时0.6631 ln (1346 )（13）m0 10yAPI 2.1

10、o页脚页眉当 o API38.3时61.933y API（14）mo0.0943yng 的计算；首先计算泡点压力系数；8.0558p g（15）xg5273.15)10(t当 x g3.448 时yng 0.3531 lnxg（16）0.5967当0.7 xg3.448时yngxg0.2401 ln0.27（17）当 xg0.7 时yng0.1236lnxg（17-1 ）0.1223如果计算出来的溶解油气比大于生产油气比，则等于生产油气比。(5) 油水混合液体的密度lo (1 f w) w fw（18）式中， f w 体积含水，小数。(6) 液体粘度原油粘度“死油”（脱气油）粘度：10 x1o

11、d（19）1000xy (32 1.8t) 1.163式中， y10 Zz3.03240.02023 yAPI页脚页眉“活油”（饱和油）粘度：A(1000OD ) BO1000式中， A10.715(5.615RS100) 0 .515B5.44(5.615RS150) 0.338OD ,O 原油死油与活油粘度，Pa s 。水的粘度1.0031.479 10 2 (32 1.8 t ) 1.98210 5 (32 1.8t )2ew1000式中，w 水的粘度， Pa s 。液体的粘度lo (1 f w )wfw(7)油、天然气的表面张力742.40.047(1.8t32)0.267( yAPI

12、 )e 1.015 10 Pog1000式中，og 油、气的表面张力，N m ；P压力 ,Pa 。(8) 水、天然气的表面张力wg 248 1.8t (23.33)137.78 (137.78) 206其中776e 3.62575 10p( 23.33)100052.58.701810 7p(137.78 )1000（20）（ 21）（22）（23）页脚页眉式中，wg 温度为 t时水、气的表面张力，N m ；P压力 ,Pa 。(9) 油水混合物和天然气的表面张力log(1f w)wgf w(10) 天然气的压缩因子 ZTc92.22176.67gPc1064.880.39ggPc1064.78

13、0.25gg273.15tTrTcPPrPc0.70.70.27 Pr（迭代 -1）PRZTr1.04670.57830.02320.68152）Z 1 (0.315063 )PR (0.53531.55) PRPR（迭代 -21.551.551.553此处迭代过程 Z 取初值 1，迭代 1 次。式中， P压力 ,Pa 。(11) 天然气的密度g3.4844 10 3g PKg / m3 ;Z (t273.15)式中， P压力 Pa 。(12) 天然气的粘度天然气粘度取g1.22 10 2 mPas页脚页眉14 井筒温度场计算根据经验公式计算沿井筒的温度分布：t t otrtoBATA L 1

14、e B ATA ( HL )（24）BATAHBATA2K P(1FW )GK P1G1.15735.4246e 1000QL1000G24式中， QL 油井产液量， t/d ；Fw 质量含水率，小数，体积含水率与质量含水率的换算公式为：Fww f wo (1 f w )w f wt o 恒温层温度，；t r 油层温度，；H油层中部深度， m ；L井筒中任意点深度， m 。1.5 井筒多相流计算井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：dpmvmdvmf mmvm2m g sindhd(

15、25)dh2页脚页眉式中m 为多相混合物的密度； vm 为多相混合物的流速；f m 为多相混合物流动时的摩擦阻力系数； d 为管径； p 为压力； h 为深度； g 为重力加速度；为井斜角的余角。多相垂直管流压力分布计算步骤根据多相管流的压力梯度就可计算出沿程压力分布。由于多相管流中每相流体影响流动的物理参数( 密度、粘度等 )及混合物密度和流速都随压力和温度而变，沿程压力梯度并不是常数。因此，多相管流需要分段计算，并要预先求得相应段的流体性质参数。然而，这些参数又是压力和温度的函数，压力却又是计算中需要求得的未知数。所以，多相管流通常采用迭代法进行计算。有两种不同的迭代途径：按深度增量迭代和

16、按压力增量迭代。按深度增量迭代的步骤已知任一点 ( 井口或井底 ) 的压力 P 作为起点，任选一个合适的压力降P作o为计算的压力间隔。具体要根据流体流量( 油井的气、液产量 ) 、管长 ( 井深 ) 及流体性质来定。估计一个对应P 的深度增量h ，以便根据温度梯度估算该段下端的温度T1 。计算出该管段的平均温度T 及平均压力P ，并确定在该 T 和 P 下的全部流体性质参数 ( 溶解气油比Rg 、 原油体积系数Bo 和粘度 o 、 气体密度g和粘度g ，混合物粘度m 及表面张力 等 )。计算该段的压力梯度dP 。dh计算对应于P的该段管长 ( 深度差 )h。计将第步计算得的计与第步估计的h 进

17、行比较，两者之差超过允许范h围，则以新的h 作为估算值，重复的计算，使计算的与估计的h 之差在允页脚页眉许范围内为止。计算该段下端对应的深度iiL及压力 PiLihjPiPo i Pj1i=1,2,3,n以i处的压力为起点，重复步，计算下一段的深度i 1i 1，直LL和压力 P到各段的累加深度等于或大于管长( LnL ) 时为止。按压力增量迭代的步骤已知任一点( 井底或井口)的压力 P0 ， 选取合适的深度间隔h ( 可将管L 等分为 n 段 ) 。估计一个对应于计算间隔h 的压力增量P 。计算该段的T 和 P，以及 P 、 T 下的流体性质参数。计算该段压力梯度dPdh o计算对应于h的压力

18、增量idPPh dho比较压力增量的估计量P与计算值P，若二者之差不在允许范围内，则i以计算值作为新的估计值，重复第步，使两者之差在允许范围o 之内为止。计算该段下端对应的深度Li 和压力 PiiPiPoiPiihL1以 Li 处的压力 Pi为起点压力重复第步，计算下一段的深度Li 1 和压力P，直到各段累加深度等于或大于管长L时为止。i1为了简化计算，通常对各段选取同样的增量间隔。而在有些情况下，各段的增量间隔可以不同，这样既能节约计算时间，而又能较好地反映出压力分布。计算气 -液两相垂直管流的Orkiszewski方法页脚页眉本设计井筒多相流计算采用Orkiszewski方法。Orkisz

19、ewski法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图1 所示。在处理过渡性流型时，采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。图 1 气液混合物流动型态 (Orkiszewski)1. 压力降公式及流动型态划分界限由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。由式(2-36) 可直接写出多项垂直管流的压力降公式：dP f dh g mdh mvmdvm(26)式中，P 压力，Pa ；f 摩擦损失梯度，Pa/m ；h 深度， m ；g 重力加速度，m/s 2 ；m 混合物密度，kg/m

20、3 ；vm 混合物流速，m/s 。动能项只是在雾流情况下才有明显的意义。出现雾流时，气体体积流量远大页脚页眉于液体体积流量。根据气体定律，动能变化可表示为：mvmdvmWt qgdp(27)2Ap P式中，pm2 ；A 管子流通截面积，Wt 流体总质量流量，kg/s ；qg 气体体积流量，m 3 /s 。将式 (27) 代入式 (26)，并取 dhhk ， dPpk ， mm ， P P 经过整理后可得：kmgf hk(28)PWt qg12Ap P式中，kPa ；P 计算管段压力降，hk 计算管段的深度差，m ；P 计算管段的平均压力，Pa 。表 1流型界限流动型态界限泡流qgLBqt段 塞

21、 流qgLSLB , vgqt过 渡 流LMvgLS雾流vgLM不同流动型态下的m和f 的计算方法不同，为此，计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表1 所示。2. 平均密度及摩擦损失梯度的计算由于不同流动型态下各种参数的计算方法不同，下面按流型分别介绍。(1) 泡流页脚页眉平均密度mH L LH g g 1 H g L Hg gH LH g 1式中，H g 气相存容比(含气率 ) ，计算管段中气相体积与管段容积之比值；H L 液相存容比( 持液率 ) ，计算管段中液相体积与管段容积之比值；g、m 在P、 T 下气、液和混合物的密度，kg/m 3 。气相存容比由滑脱速度Vs

22、 来计算。滑脱速度定义为：气相流速与液相流速之差。vsgvsLqgqt qgvs1 H gAp H gAp (1 H g )H g可解出H g ：H gt(1t) 24qg11 qqvs Ap2 vs Apvs Ap式中， vs 滑脱速度，由实验确定，m/s ；vsg 、 vsL 气相和液相的表观流速，m/s 。泡流摩擦损失梯度按液相进行计算：(29)(30)(31)tfLvLH2D2vLHqL(32)Ap(1H g )式中， f 摩擦阻力系数；vLH 液相真实流速，m/s 。摩擦阻力系数 f可根据管壁相对粗造度/ D 和液相雷诺数NRe 查图 2 。对于普通油管，其管壁绝对粗糙度，一般取4.

23、57 10 5 m 。液相雷诺数页脚页眉Dv sLL(33)N ReL式中，L 在P、T 下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值，Pa.s 。(2) 段塞流混合物平均密度WtL vs Ap(34)mLqtvs Ap式中，液体分布系数；v s 滑脱速度，m/s 。滑脱速度可用Griffith和 Wallis 提出的公式计算：页脚图 2摩擦阻力系数曲线页眉vs C1C2 gD(35)(3) 过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。LMvgvgLsm LMLSSLLMLSMitLMvgSLvgLgMiLMLSLML

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

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

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

27、、最小载荷的计算公式：页脚页眉iSN2Pmax(WrjWL )(1)j 11790iiWrjqrj Lrj gj1j1WLf p (PZ PN )式中， qr i第 i 级杆每米杆在空气中的质量， Kg/m ；Lri第 i 级杆杆长， m ；i 抽油杆级数，从下向上计数；PZ泵排出口压力， Pa；PN泵的沉没压力， Pa ；N冲次， rpm ；S光杆冲程， m ；fP活塞截面积， m 2；g 重力加速度， m/s 2。iSN2iPminWrjWrj1790 jj11iiiWrjWrjPj( f rjf r 1 j 1 )j 1j 1j 1令 f r0 =0式中， Pj第 j 级抽油杆底部断面处

28、压力，Pa ：j1PjPt0(1f w )w f w ? g(LLt )t 1Pt井口油压， Pa;地面油密度， kg/m 3；0fw 体积含水率，小数。（40）（41）（42）（43）（44）（45）页脚页眉应力范围比 pL 计算公式：PLmaxmin（46）allminPm axPm in（47）m axm infrfr抽油杆柱的许用最大应力的计算公式：all ( T0.5625min) SF4式中，all 抽油杆许用最大应力，Pa ；T抽油杆最小抗张强度，对C 级杆， T=6.3*10 8 Pa, 对 D 级杆 T=8.1*10 8Pa;min 抽油杆最小应力， Pa ；SF 使用系数，

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

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

31、 76mm ），原因是作业时大柱塞不能下如小直径油管中；若采用25mm 抽油杆，则相应油管直径应用3 ，原因是25mm 抽油杆节箍为55mm ，与 62mm 油管间隙太小。当采用多级杆时3 油管长度比25mm 杆长多 10m 。为了减小计算工作量，在本次课程设计中杆柱设计简化处理，采用单级杆设计（ 19mm ）。1.7 抽油机校核最大扭矩计算公式M max 1800S 0.202S(Pmax P min)（48）式中， M max 最大扭矩， Nm ；Pm ax悬点最大载荷， N；P min 悬点最小载荷， N ；S冲程， m 。电动机功率计算，Nt1000M maxn14388（49）式中，

32、 N t 需要的电动机功率，W；页脚页眉n 冲数， rpm ；如果计算的最大扭矩超过抽油机所配减速箱允许的最大扭矩或计算电动机功率超过电动机额定功率则需改变抽油参数（fp ，s，N 及 L）重行进行设计计算。1.8 泵效计算泵效及其影响因素在抽油井生产过程中， 实际产量 Q 一般都比理论产量Qt 要低，两者的比值叫泵效，表示，Q（50）Qt产量计算根据影响泵效的三方面的因素，实际产量的计算公式为SPqleak（51）Q QtBlS Bl式中， Q 实际产量， m 3/d;Qt 理论产量， m 3/d;Sp柱塞冲程， m ；S光杆冲程， m ；Sp S 抽油杆柱和油管柱弹性伸缩引起冲程损失系数；

33、Bl泵内液体的体积系数；泵的充满系数；q leak 检泵初期的漏失量， m 3/d;理论排量计算Ql 1440 f p SN（52）页脚页眉式中， Qt 泵的理论产量， m 3/d 。冲程损失系数 Sp S 的计算根据静载荷和惯性载荷对光杆冲程的影响计算本设计按照油管未锚定计算。当油管未锚定时：SP(1u2WlL1L2L3L（53）S)SE(f r 2f r 3)2f r 1f t当油管锚定时：SP(1u 2WlL1L2L3)（54）S)(f r 2f r 32SE f r1其中u L/a式中，曲柄角速度， rad/s ；N/30 ；a 声波在抽油杆柱中的传播速度， 5100m/s ；Wl 考虑沉没度影响后的液柱载荷为上下冲程中静载荷之差，N；Wl (PZPin ) ? f plL f gf p（55）PZ泵排出口压力， Pa ；P 泵内压力， Pa ；当液体粘度较低时，可忽略泵吸入口压力，故P P ；ininNPN泵的沉没压力， Pa ；fp 、fr、ft 活塞、抽油杆及油管金属截面积，m 2 ；L抽油杆柱总长度， m ； l 液体密度，kg/m 3;E钢的弹性模数，112.06 Pa10；页脚页眉Lf动液面深度， m