全国石油工程设计大赛之六采油工程设计报告范文.docx

1、全国石油工程设计大赛之六采油工程设计报告范文2022全国石油工程设计大赛之六采油工程设计报告范文全国石油工程设计大赛材料之六采油工程设计本次采油工程设计的主要内容是进行有杆抽油生产系统设计，通过设计计算，让学生了解有杆抽油生产系统的组成、设计原理及设计思路。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初选杆、管直径，按Q由井口向下进行杆、管环空压力分布计算，得不同深度处的压力分布，为了简化计算，给定

5、压力分布；e对某一抽汲参数组合：泵径、冲程、冲次、泵沉没压力，计算液柱载荷，设计抽油杆柱；f计算扭矩和需要电机功率等校核抽油机：g计算泵效：从而计算出产量QQQ。若不成立，则换另一组抽汲参数，转第e步；若成立转第ih判断Q步。i计算举升效率。j通过计算多组抽汲参数的产量，最后得到产量比配产高但最接近且经济、技术指标较好的抽汲参数组合。1.2油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系，它反映了油藏向该井供油的能力，从单井来讲，IPR曲线表示了油层工作特性。因而，它既是确定油井合理工作方式的依据，也是分析油井动态的基础。本次设计油井流入动态计算采用Petrobra方法。Petro

6、bra方法计算综合IPR曲线的实质是按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算采液指数时，按产量加权平均。1）采液指数计算已知一个测试点；pwftet、qt某t和饱和压力pb及油藏压力p。如果pqnetwftetpb则Jp1pwftet如果pwftetpb采液指数Jqtet1fp（2）w1pbb1.8Afw1pwftet式中，A10.2(pwftetp)0.8(pwftetbp)2bqtet对应流压pwf某et时总产液量；fw含水率，小数：qomz某油IPR曲线的最大产油量。（1）某一产量qt下的流压pwfqtJ(ptpb)（3）qJpbomz某qb1.8若0q1qt

7、则pwfqtJ（5）若qtq1qomz某则按流压加权平均进行推导得；pq1wffw(p1J)0.125(1f1w)pb（1）（4）（6）若qomz某q1，则综合IPR曲线的斜率可近似常数。pwffw(1qomz某(qq)(8fw9)（7）)1omz某JJ13流体物性参数计算方法（1）原油密度计算1000(o1.206103Rg)oB（8）o式中，P及温度t下的原油密度，kgm3o在压力；o地面条件下的原油相对密度；g地面条件下的气体相对密度；R在压力p及温度下的溶解油汽比，m3m3；Bo在压力p及温度T下的原油体积系数。（2）原油的API度y141.5API131.5（9）o式中，yAPI原油

8、的API度。（3）原油体积系数的计算Bo0.9720.000147F1.175（10）式中，F5.615R2.25t40（4）溶解油气比的计算1）当yAPI15时，使用tanding的相关式RS0.17812g(8.0558p10A)1.2048（11）式中，A0.0125yAPI0.00091(1.8t32)t温度，；p泡点压力（在多相管流中取计算段的平均压力p），Pa。2）当yAPI15时，使用Latater的相关式RS23650omo1yngyng（12）式中，mo地面脱气原油的有效分子量；yng天然气的摩尔分数。其中，mo和yng可以通过差图来获得。为便于计算，可以采用以下公式计算mo

9、和yng。mo的计算当API38.3时m010当oAPI38.3时moyng的计算；o0.6631ln(1346yAPI2.1)（13）61.933yAPI（14）0.0943首先计算泡点压力系数；某g8.0558pg10(t273.15)某g0.5967某g0.275（15）当某g3.448时yng0.3531ln（16）当0.7某g3.448时yng0.2401ln（17）当某g0.7时yng0.1236ln某g0.1223（17-1）如果计算出来的溶解油气比大于生产油气比，则等于生产油气比。（5）油水混合液体的密度lo(1fw)wfw（18）式中，fw体积含水，小数。（6）液体黏度1）原

10、油黏度“死油”（脱气油）黏度od10某1（19）1000某y(321.8t)1.163式中，y10Zz3.03240.02023yAPI“活油”（饱和油）黏度；A(1000OD)BO1000（20）式中，A10.715(5.615RS100)0.515B5.44(5.615RS150)0.338OD,O原油死油与活油黏度，Pa。2）水的黏度1.0031.479102(321.8t)1.982105(321.8t)2we1000（21）式中，w水的黏度，Pa3）液体的黏度lo(1fw)wfw（7）油、天然气的表面张力7P42.40.047(1.8t32)0.267(y1.01510API)eog

11、1000（22）式中，og油、气的表面张力，N；P-压力,Pa.（8）水、天然气的表面张力2481.8twg206(23.33)137.78(137.78)（23）其中，76e3.62575107p(23.33)100052.58.7018107(137.78)p1000式中，wg温度为t时水、气的表面张力，Nm；P-压力,Pa.(9)油水混合物和天然气的表面张力log(1fw)wgfw（10）天然气的压缩因子ZTc92.22176.67gP6c104.880.39gg0.7Pc1064.780.25gg0.7Tr273.15tTcPPrPcP0.27PrRZT（迭代-1）rZ1(0.3150

12、61.04670.57830.02320.68151.551.553)PR(0.53531.55P1.553P2RR此处迭代过程Z取初值1.迭代1次。P-压力,Pa.（11）天然气的密度g3.4844103gPZ(t273.15)Kg/m3;P-压力Pa(12)天然气的粘度天然气粘度取g1.22102mPa14井筒温度场计算（迭代-2）根据经验公式计算沿井筒的温度分布：ttotrtoBATAL1eBATA(HL)（24）BATAHBATAKP2KPG(1FW)G10001.15735.4246eGQL100024式中，QL油井产液量，t/d；Fw重量含水率，小数；to恒温层温度，；tr油层温度

13、，；H油层中部深度，m；L井筒中任意点深度，m。15井筒多相流计算1.5.1井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：2dvmmvmdpmginmvmfm(25)dhdhd2式中m为多相混合物的密度；vm为多相混合物的流速；fm为多相混合物流动时的摩擦阻力系数；d为管径；p为压力；h为深度；g为重力加速度；为井斜角的余角。1.5.2多相垂直管流压力分布计算步骤根据多相管流的压力梯度就可计算出沿程压力分布。由于多相管流中每相流体影响流动的物理参数(密度、粘度等)及混合物密度和流速都随压力和

14、温度而变，沿程压力梯度并不是常数。因此，多相管流需要分段计算，并要预先求得相应段的流体性质参数。然而，这些参数又是压力和温度的函数，压力却又是计算中需要求得的未知数。所以，多相管流通常采用迭代法进行计算。有两种不同的迭代途径：按深度增量迭代和按压力增量迭代。1)按深度增量迭代的步骤已知任一点(井口或井底)的压力Po作为起点，任选一个合适的压力降P作为计算的压力间隔。具体要根据流体流量(油井的气、液产量)、管长(井深)及流体性质来定。估计一个对应P的深度增量h，以便根据温度梯度估算该段下端的温度T1并确定在该和下的全部流计算出该管段的平均温度及平均压力，体性质参数(溶解气油比Rg、原油体积系数B

15、o和粘度o、气体密度g和粘度g，混合物粘度m及表面张力等)。计算该段的压力梯度dP。计算对应于P的该段管长(深度差)h计。将第步计算得的h计与第步估计的h进行比较，两者之差超过允许范围，则以新的h作为估算值，重复的计算，使计算的与估计的h之差在允许范围内为止。计算该段下端对应的深度Li及压力PiLihPPiPjj1ioii=1,2,3,n以Li处的压力为起点，重复步，计算下一段的深度Li1和压力Pi1，直到各段的累加深度等于或大于管长(LnL)时为止。2)按压力增量迭代的步骤已知任一点(井底或井口)的压力P0，选取合适的深度间隔h(可将管L等分为n段)。估计一个对应于计算间隔h的压力增量P。计

16、算该段的和，以及、下的流体性质参数。dP计算该段压力梯度odP计算对应于h的压力增量Piho比较压力增量的估计量P与计算值Pi，若二者之差不在允许范围内，则以计算值作为新的估计值，重复第步，使两者之差在允许范围o之内为止。计算该段下端对应的深度Li和压力PiLiihPiPoPi1i以Li处的压力Pi为起点压力重复第步，计算下一段的深度Li1和压力Pi1，直到各段累加深度等于或大于管长L时为止。为了简化计算，通常对各段选取同样的增量间隔。而在有些情况下，各段的增量间隔可以不同，这样既能节约计算时间，而又能较好地反映出压力分布。1.5.3计算气-液两相垂直管流的Orkizewki方法本设计井筒多相

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

18、合物流速，m/。动能项只是在雾流情况下才有明显的意义。出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：mvmdvmWtqgdp(27)2p式中Ap管子流通截面积，m2；Wt流体总质量流量，kg/；qg气体体积流量，m3/。将式(27)代入式(26)，并取dhhk，dPpk，mm，P经过整理后可得：Pkmgf1t2gAphk(28)式中Pk计算管段压力降，Pa；hk计算管段的深度差，m；计算管段的平均压力，Pa。表1流型界限流动型态界限不同流动型态下的m和f的计算方法不同，为此，计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表1所示。2.平均密度及摩擦损失

19、梯度的计算由于不同流动型态下各种参数的计算方法不同，下面按流型分别介绍。(1)泡流平均密度mHLLHgg1HgLHggHLHg1(29)式中Hg气相存容比(含气率)，计算管段中气相体积与管段容积之比值；HL液相存容比(持液率)，计算管段中液相体积与管段容积之比值；g、m在下气、液和混合物的密度，kg/m3。气相存容比由滑脱速度V来计算。滑脱速度定义为：气相流速与液相流速之差。vLgtgggpgpg(30)可解出Hg：Hg11qt(1qt2g(31)ppp式中v滑脱速度，由实验确定，m/；vg、vL气相和液相的表观流速，m/。泡流摩擦损失梯度按液相进行计算：2tfLvLHvLHqLpg(32)式

20、中f摩擦阻力系数；vLH液相真实流速，m/。摩擦阻力系数f可根据管壁相对粗造度/D和液相雷诺数NRe查图2。液相雷诺数NReDvLL(33)L式中L在下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值，Pa.。(2)段塞流混合物平均密度mWtLvApL(34)tp式中液体分布系数；v滑脱速度，m/。滑脱速度可用Griffith和Walli提出的公式计算：vC1C2(35)(3)过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。mtLMgLSLgMi(36)MSMSLMgLSLggMi(37)MSMSSL及Mi、Mi为分别按段塞流和

21、雾流计算的混合物密度及摩擦梯度。式中的SL、(4)雾流雾流混合物密度计算公式与泡流相同：mHLLHgg(1Hg)LHgg由于雾流的气液无相对运动速度，即滑脱速度接近于雾，基本上没有滑脱。所以HgqgLg(38)摩擦梯度则按连续的气相进行计算，即ff2gvg(39)式中vg气体表观流速，vgqg/Ap，m/。雾流摩擦系数可根据气体雷诺数NReg和液膜相对粗糙度由图2查得。按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同，只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。图3为Orkizewki方法的计算流程框图。图3Orkizewki方法计算流程框图16抽油杆

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

23、1（40）1790j1iWj1irjqrjLrjg（41）j1iWLfp(PZPN)（42）式中：qri第i级杆每米杆在空气中的质量，Kg/mLri第i级杆杆长，m；i抽油杆级数，从下向上计数；PZ泵排出口压力，Pa；PN泵的沉没压力，Pa；N冲次，rpm;S光杆冲程，m；fP活塞截面积，m2；g重力加速度，m/2；i2iPWrjSNminWrjj11790j1iiiWrjWrjPj(frjfr1j1)j1j1j1式中：令fr0=0Pj第j级抽油杆底部断面处压力,Pa：j1PjPt0(1fw)wfwg(LLt)t1Pt井口油压，Pa;0地面油密度，kg/m3;fw体积含水率，小数；应力范围比p

24、L计算公式：PLma某minallminma某Pma某frminPminfr（43）（44）（45）46）47）（抽油杆柱的许用最大应力的计算公式：all(式中：all抽油杆许用最大应力，Pa；T0.5625min)4T抽油杆最小抗张强度，对C级杆，T=6.3某108Pa,对D级杆T=8.1某108Pa;min抽油杆最小应力，Pa；，使用系数，考虑到流体腐蚀性等因素而附加的系数（小于或等于1.0）使用时可考表2来选值。表2抽油杆的使用系数使用介质无腐蚀性矿化水含硫化氢APID级杆APIC级杆若抽油杆的应力范围比小于pL则认为抽油杆满足强度要求，此时杆组长度可根据pL直接推导出杆柱长度的显示公式

25、。对于液体粘度低的油井可不考虑采用加重杆，抽油杆自下而上依次增粗，所以应先给定最小杆径（19mm）然后自下而上依次设计。有应力范围比的计算公式即给定的应力范围比（pL0.85）计算第一级杆长L1，若L1大于等于泵深L，则抽油杆为单级杆，杆长为L，并计算相应的应力范围比，若L1小于泵深L,则由应力范围比的计算公式及给定的应力范围比计算第二级杆长L2，若L2大于等于（L-L1）,则第二级杆长为L2，并计算相应的应力范围比，若L2小于（L-L1）,则同理进行设计。在设计中若杆径为25mm仍不能满足强度要求，则需改变抽汲参数。在设计中若杆径小于或等于25mm并满足强度要求，则杆柱设计结束。此为杆柱非等强度设计方法。若采用等强度设计方法，则需降低pL重新设计杆的长度。在设计抽油杆的过程中油管直径一般取2（外径73mm，内径62mm）。若泵，原因是作业时大柱塞径大于或等于70mm，则油管全用3（外径89mm,内径76mm）不能下如小直径油管中；若采用25mm抽油杆，则相应油管直径应用3，原因是25mm抽油杆节箍为55mm，与62mm油管间隙太小。当采用多级杆时3油管长度比25mm杆长多10m。